Page 1 of 2

CO2 injection in Caprock- Aquifer

Posted: Sat Jan 23, 2016 2:59 pm
by Taiyeb Badr
I have tried modelling this tutorial problem on Code_bright as given in the tutorial available on the website of Code_bright. When I run the calculation, it ends within 3 seconds and pops up flashing "successful completion" . As expected, there is no result available in the post-process. I think the solver terminates within the first interval itself (The problem has 3 intervals). I have tried giving smaller time steps but it did not help. Where am I going wrong ?

Re: CO2 injection in Caprock- Aquifer

Posted: Mon Jan 25, 2016 8:39 pm
by code_bright
Hi!

Please find attached the tutorial project.

Hope it helps.

Regards,

Alfonso

Re: CO2 injection in Caprock- Aquifer

Posted: Thu Jan 28, 2016 10:52 am
by Taiyeb Badr
Thank you for the swift response but I am trying to model the hydro-mechanical analysis, the one with caprock . The attached file is for hydraulic analysis alone. My query is related to the next tutorial problem which has a caprock and solves stress equilibrium equations. It will be of great help if you could attach the correct file.

Re: CO2 injection in Caprock- Aquifer

Posted: Fri Jan 29, 2016 11:38 am
by code_bright
You are right. Here it is the right one.

Regards,

Alfonso

Re: CO2 injection in Caprock- Aquifer

Posted: Thu Mar 10, 2016 1:10 pm
by Taiyeb Badr
My question is related to the CO2 properties given in case of injection. The carbon dioxide properties are given at the flux bc in terms of gamma for gas. What exactly is this gamma? I guess it is not specific heat ratio as its value is 1000 in some cases in the tutorial problems. Also how does CO2 differ from air in terms of the input values for gas?
Further, how do I vary the CO2 properties?

Re: CO2 injection in Caprock- Aquifer

Posted: Mon Mar 14, 2016 6:29 pm
by code_bright
Please find attached a brief explanation for the gamma for gas.

Re: CO2 injection in Caprock- Aquifer

Posted: Tue Mar 15, 2016 12:41 pm
by Taiyeb Badr
Thank you for the explanation of gamma. But how does one change the carbon dioxide properties? Lets say that I wish to inject any other gas like methane or change the state of CO2 from supercritical to gaseous and vice versa. What input parameters determine this change?

Also, this question is for the code_bright tutorial team. How do you come up with the gas influx mass rate? How to calculate it ? Is there a factor used in multiplication ? In my calculations, I was off by about a factor of 10 in the kg/s rate of influx. Can you shed some light on this as well.

Re: CO2 injection in Caprock- Aquifer

Posted: Tue Mar 15, 2016 5:08 pm
by code_bright
You have that information in the user's guide.

For C02, you have several laws you can use:

ICL=15 (LIQUID DENSITY) --> You have one specific law for CO2: ITYCL=4
ICL=17 (GAS DENSITY) --> ITYCL=4
ICL=18 (GAS VISCOSITY) --> ITYCL=4
ICL=12 (DIFFUSIVE FLUXES) --> ITYCL=1

Then, if you want to use another gas instead CO2, you can use more general laws, like:

ICL=17 (GAS DENSITY) --> ITYCL=3
ICL=18 (GAS VISCOSITY) --> ITYCL=2

You have all the parameters that you need explained in the user's guide, and also some references where you can find more information about the laws.

Regards,

Alfonso

Re: CO2 injection in Caprock- Aquifer

Posted: Wed Mar 16, 2016 12:31 pm
by Taiyeb Badr
Also, this question is for the code_bright tutorial team. How do you come up with the gas influx mass rate? How to calculate it ? Is there a factor used in multiplication ? In my calculations, I was off by about a factor of 10 in the kg/s rate of influx. Can you shed some light on this as well.
example 0.84 kg/s is off by a factor of ten i.e. ~0.09 kg/s (in my calculations )to reach the per year injection requirement . Please refer the tutorial problems. Can you tell me how to convert this ton/year injection rate to kg/s ?

Re: CO2 injection in Caprock- Aquifer

Posted: Wed Mar 16, 2016 6:35 pm
by code_bright
Hi!

For the calculation of the gas flow, you have to take into account the surface where you are applying the boundary condition. Remember that it is an axisymmetric problem.

Surface = 100 * 2 * pi* 0.15 = 94,25 m2

10000 ton/week = 1.65 kg/s

gas flow = 1.65/94.25 = 0.175

Regards,

Alfonso