Specifying initial conditions for layered reservoir
I am trying to model CO2 injection in stratified saline aquifer reservoir. I have to provide initial conditions for my model considering hydrostatic pressure gradient and geothermal temperature gradient. I can obtain the initial hydrostatic pressure values by running the model with single pressure value for all layers (without wells) to obtain the hydrostatic pressure distribution. However for initial temperature, I have to specify the initial temperature for every layer manually. I would like to know if there is function/code in python that can specify initial temperature for every layer directly.
Thanks in advance.
is there any pattern to the temperature distributions, or do you want to have the ability to impose arbitrary temperatures to initialize your simulation? If the former, there may be tricks to run TOUGH to get the effects you want. If the latter, I’d suggest one of the preprocessors, like toughio developed by Keurfon Luu
Do I need to specify large grid volume for the upper and lower layers to keep temperature there constant during initialization? I expect after initialization that temperature at upper layer will be higher than 50 degrees C and the temperature of the bottom layer is lower than 57.5 degrees C. Am I right?