Residual and solubility trapping
In the result we have
1-sg which represents the saturation of CO2 in gas phase.
2- sliq which is saturation of Co2 in liquid phase
3- saq represent the water saturation
4- Xco2N is it represent the mole fraction of disloved CO2 in water?
If I am right? How I can calculate the amount of residual trapping and solubility trapping..
I have not used ECO2M in a while, so if anyone else has better information, please chime in. But here is what I suggest:
For Solubility Trapping (mass of CO2 dissolved in the aqueous phase) there are two possibilities and they are the same for ECO2N and ECO2M.
1. Look at the VOLUME- AND MASS_BALANCES section of the output file. There is a table showing the mass of each component in each phase, for individual materials, and for the problem as a whole. The mass of dissolved CO2 in the aqueous phase is in the AQUEOUS row, in the CO2 column.
2. In the main output (OUTPUT DATA), for each grid block, mass fraction of CO2 is given as XCO2AQ. If you multiply this by aqueous saturation and aqueous density, you will get the mass of dissolved CO2 in that grid block. In ECO2N, aqueous saturation is 1-SG-SS and aqueous density is DL. In ECO2M, aqueous saturation is SAQ and aqueous density is Daq (I think you need to specify KDATA=3 to get Daq to appear in a separate section of the output file).
For Residual Trapping, you need to use the main output and check for each grid block whether gas saturation is less than residual gas saturation. In ECO2N, gas saturation is SG. In ECO2M, gas saturation is SGAS. For a non-hysteretic simulation, residual gas saturation is an input parameter in the relative permeability function (could be different for different materials). For a hysteretic simulation (ECO2N only), residual gas saturation is calculated dynamically as part of the TOUGH simulation, and is variable SOR in the main output.
Hope this helps, Christine