TOUGHREACT1.2/ECO2N: About the inject rate and maximum pressure
I'm using TOUGHREACT1.2/ECO2N with Petrasim. I am simulating the injection of CO2(isothermal) at a constant rate of 6.5kg/s（just one GENER）using a rectangle mesh, but I find that it will be an error that shows after 30 days, but I want to simulate for years. And I finally found that the P around my GENER ELEMNT is about 59MPa which calculates 30 days, so I'm wondering if it's overpressurized. （The manual shows the maximum pressure is 600bar）
if it does， How should I set up my GENER items or change some settings to let my high inject rate run successfully？
Attached is my Petrasim file，I urgently need your help, thank you very much!
Your output file does not show any problem for the 30 days simulated. The P builds up to 596.33 bara, as you mention, for a relatively limited injection rate of 5 kg/s. And you have not linked permeability to the precipitation of solid salt, which actually occurs around the injection well. The error you mention after 30 years of simulation may be linked to P going above 600 bar. I have not checked that by running the model.
You inject in eleme " 1" belonging to domain "sand ", which as permeabilities kx, ky and kz of 5.8E-16, 5.8E-16, and 5.8E-18 m2, respectively.
In my view, you are simply trying to inject a too high rate into a very low permeability formation: less than 0.6 mD ! If the assigned permeability is representative of true formation properties, then you have to reduce your injection rate to levels compatible with the very poor formation properties.
Or considering to inject in a different formation...
You are using a 2D vertical slice just 1 m wide. Injecting 5 kg/s CO2 in a vertical slice is not comparable to inject the same rate in a radial system. I think you should either use a radial grid or scale the injection rate from the radial injection to that on the vertical slice.
The aquifer domain "sands" has a considerable thickness of some 550 m, but you inject only in the bottom element 50 m thick. I imagine you wanted to simulate the upward migration of CO2 under gravity driven flow.
But you assigned a very low vertical permeability, initially with a kx/kz ratio of 100. Thus, your injected CO2 is constrained into a 1 m wide 50 m thick volume with an almost linear flow. Under these conditions, there might be anything wrong in the simulation, but a too high CO2 injection rate for the system permeability and geometry.
Hello, can you share a TOUGHREACT1.2 source code?Thank you very much