For the case I'm working on, I've created and INCON file filled with cell-by-cell porosities, permeabilities (I,J,K), pressure and temperature.
I run it for a few years in order to reach a more precise equilibrium, and then want to use the generated SAVE file as a new starting point for my actual run.
But in the SAVE file, I miss the writing of permeabilities. How can I get them reproduced in the SAVE file ?
I have to clarify this to avoid (further) confusion:
iTOUGH2 does have element-by-element permeabilities or permeability modifiers (for all three directions) and other user-specified parameters in the INCON and SAVE files.
Andri asked you the right questions. Here are the possibilities:
(1) You are using TOUGHREACT (but then you would not have posted the question)
(2) You are using iTOUGH2 (but then Andri's suggestion to set MOP(13)=1 would have worked).
(3) You are using TOUGH2 (the most likely scenario)
If you send one of your output files, as Andri asked, we would see immediately from the header which TOUGH package you are using (you may also check your license agreement).
If you are using TOUGH2, then specifying permeabilities in INCON has no effect. You can easily check this by comparing the output with and without the extra information in INCON.
If your permeabilities do not change during the simulation (e.g., through salt precipitation*), and all you want is (static) heterogeneity, then simply specify permeability modifiers in block ELEME (Columns 41-50) rather than INCON, and add a material SEED (check TOUGH2 user's guide for details).
* You can have salt precipitation in certain TOUGH2 modules, but these dynamic permeability changes do not need to be saved in the SAVE file, as they are internally recalculated as a function of solid saturation. In fact, only TOUGHREACT and iTOUGH2 include processes that lead to dynamic permeability changes.
Glad you are charging ahead.
Just to complete the topic and answer Eric's question:
TOUGH3 inherited some (but not all) of the iTOUGH2 features documented in the report "Enhancements to the TOUGH2 Simulator Integrated in iTOUGH2".
So yes, TOUGH3 does indeed read user-specified parameters (specifically permeabilities / permeability modifiers) from block INCON using the same format and approach as in iTOUGH2. Note that the format used in TOUGHREACT is different from that used by iTOUGH2 and TOUGH3.
I've been trying to convert Petrel/Eclipse-type grid files to MESH/INCON files. I searched the Internet and found your answer. I learned that T2B is a new tool developed by BRGM that can solve this problem. I meet some difficulties in the process of converting the grid, and I need your guidance to solve these problems. I have implemented the conversion of Petrel to a finite element mesh, but how to realize the Petrel grid to Tough grid in the"Petrel2Tough" module? Specifically, Petrel is an irregular hexahedral grid. how do you calculate VOLX(element volume), AHTX(interface area), X, Y, Z(Cartesian coordinates of grid block centers.)，ISOT(set equal to 1, 2, or 3;), D1(distance (m) from first element to the interface.), D2(distance (m) from the second element to the interface), AREAX(interface area), BETAX(cosine of the angle between the gravitational acceleration vector and the line between the two elements)? Can you provide me with a method for this problem or how to get T2B? Your help is greatly appreciated.