Calculating Total Density from TOUGH2/ECO2M Output file
Dear all,
I am sorry maybe this is a simple questions, but I need to know it. When I run injection CO2 with eco2m module for my model with initial conditions as follow : salt mass fraction=0 and CO2 mass fraction=0. And then I got the result with the following format :
ELEM INDEX I P T SAQ SLIQ SGAS XCO2AQ DLIQ VISLIQ DGAS VISGAS
I know the porosity value of each element, then I need to calculate the fluid density of each element based on those data (SAQ, SLIQ, SGAS, XCO2AQ, DLIQ, and DGAS) assumes that there is no salinity in the reservoir. Thank you very much.
Here is what I've done, is that true or not?
1. I calculate Fluid density in Aqueous phase with a dissolved CO2 (as mentioned in ECO2M Manual book), let say it RHOaq.
2. Then I calculate the total fluid density by
RHOfluid = SAQ*RHOaq + SLIQ*DLIQ + SGAS*DGAS
3. Then I can calculate the total density by following expression:
RHOtot = RHOgrain*(1-porosity)+porosity*RHOfluid.
RHOgrain is material rock density that be assigned in DROK block in TOUGH2 input file.
I need to know whether the calculation for fluid density is true or not? Thank you very much in advance.
4 replies
-
Suhendi,
That looks right to me. When you say "salt mass fraction=0", are you referring to the absence of solid salt or the absence of any salt (i.e., solid and dissolved in aqueous phase)?
Thanks!
-Mikey
-
Hi Mikey,
Thank you for your reply. Yes, I mean in initial condition "salt mass fraction = 0" so I think that the rock is fully saturated by water, so any salt is absence. The main question is how to calculate the total density of each element using the TOUGH2/ECO2M output file with above format output file. I've calculated the total density of each element by those equation, but the result is less than initial density.
For your information, the permeability of lateral and bottom boundary of the model is not zero, so is it possible the fluid content go outside the model? Please do not hesitate to correct me if I'm wrong.
Thank you very much in advance,
Suhendi,
-
Suhendi,
Depending on how permeable your boundaries are, your conclusion about fluid loss from the domain (particularly through the lateral boundary) is entirely plausible. Especially considering you're displacing a denser fluid (water) with a lighter CO2.
- Mikey
-
Mikey,
Thank you very much for your answer.
Suhendi