Pressure reporting: Tgh2-ECO2M
Hello ECO2M users,
It is about pressure reporting of Tough2 when running a 1D radial model initialized with only water, flooded by CO2 at Its center using ECO2M. I get the following curves.
While expecting the pressure at initial part to be at the same level as the picked values and be continuous.
Is it a reporting issue here? It is a bit difficult to make sense of it. Input and output files are attached. Running the same model with different rel perm and capillary behavior gave negative values for the same initial part, see output file itgh_Sg_P_T_S.col.
Other outputs like saturations, temperature, etc. don't seems to have such discontinuity- at first glance.
Appreciate your feedback.
I tried your input file with T2Well-ECO2M, of course with no wellbore flow.
The P vs log(radius) at 1 year is shown below.
In your input file I noticed that you have omitted the phase index, which for stated aqueous conditions should be 1.
You have also used X3=10.00001 which is a wrong initialization, as X3 should be the mass concentration of CO2 in the aqueous phase. As a guess I used X3=0.00001.
7000000.0 0.000001 10.00001 35.0
I need to clarify my previous post. The ECOM version I used is V.2 by Lehua Pan which can simulate H2O-CO2-NaCl mixtures up to 300°C. It seems it is not able to accept the ECO2N style INCON format that you used. This is why I had to include the phase index in the INCON block, which is any way the standard ECO2M input.
But you initialized your model in 2 phase (aqueous+CO2phase), not in single aqueous conditions as you wrote, by using X3=10.00001= 10 + Sc. So, I guess the run with TOUGH3-ECO2M with ECO2N INCON style runs with the reservoir initially in two-phase conditions.