SAVE file

Hello !

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 ?



18replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
    • Andri Arnaldsson
    • Senior Geothermal Scientist - Vatnaskil Consulting Engineers
    • Andri_Arnaldsson
    • 3 yrs ago
    • Reported - view

    Do you need the permeability values in the SAVE file? Are the changed from that is in the ROCKS block?

    • Andri Arnaldsson Yes actually I don't use the ROCKS keyword for permeability definition as it only allows an homogeneous value for the whole grid. So I define them through INCON

      • Andri Arnaldsson
      • Senior Geothermal Scientist - Vatnaskil Consulting Engineers
      • Andri_Arnaldsson
      • 3 yrs ago
      • Reported - view

       Are you using TOUGH2 or iTOUGH2? If you use iTOUGH2 you could try to set MOP(13)=1

    • I think it's TOUGH2 but it's worth a try

    • No, that didn't work. I'll probably have to reprocess the Save file then

      • Andri Arnaldsson
      • Senior Geothermal Scientist - Vatnaskil Consulting Engineers
      • Andri_Arnaldsson
      • 3 yrs ago
      • Reported - view

       Can you share your input/output files?

  • Hi all,

    Only TOUGHREACT has grid-block dependent permeabilities in the INCON and SAVE files. I know for sure that TOUGH2 and ITOUGH2 do not have this feature. Unless you have a version that someone else has  modified to have this same feature.



  • Hi Eric,

    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.


    Hi Florian,

    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).

    Good luck!



    * 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.

  • I didn’t know that Stefan! Is that option in TOUGH3?

  • Thanks for the enlightenments Stefan.

    Actually I've drawn this conclusion after running a few test. Permeability modification in my INCON did not have any impact, so I more or less guessed that it was not read. 

    I'll try to use the method you suggest via ELEME and SEED.



    • (and yes, I checked and I'm using Tough2)

  • Florian,

    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.



    • Stefan Finsterle Interesting. I check with the IT department, and apprently we have a ToughReact license available, so I might try to use it . Can you give me some advice on were I can find a desription of the appropriate formatting for ToughReact ?

      • Stefan Finsterle
      • Finsterle GeoConsulting
      • Stefan_Finsterle
      • 3 yrs ago
      • Reported - view

      Florian DELSANTE 


      Please check the TOUGHREACT manual or contact Eric Sonnenthal.

      I'm actually curious to learn where you got the idea that permeabilities can be specified in INCON.


    • Stefan Finsterle The idea came up because I use a pre-processing software to convert Petrel/Eclipse type grid files to MESH/INCON files. It's called T2B, developped by BRGM.

      It reads Petrel meshes and properties grids formatted in Eclipse keywords : 

      -- Format      : ECLIPSE keywords (grid geometry and properties) (ASCII)
      -- Exported by : Petrel 2017.1 Schlumberger
      -- User name   : TP2030
      -- Date        : Thursday, January 04 2018 15:27:08
      -- Project     :
      -- Grid        : Reservoir

      PINCH                                  -- Generated : Petrel

      NOECHO                                 -- Generated : Petrel

      MAPUNITS                               -- Generated : Petrel
        METRES /

      MAPAXES                                -- Generated : Petrel
        543748.7368 2786171.823 543748.7368 2787171.823 544748.7368 2787171.823 /

      GRIDUNIT                               -- Generated : Petrel
        METRES /

      SPECGRID                               -- Generated : Petrel
        19 23 100 1 F /

      COORDSYS                               -- Generated : Petrel
        1 100 /

      COORD                                  -- Generated : Petrel
        1200 1100 158.552063 1200 1100 1158.552002 1250 1100 147.3273926 1250 1100 1147.327393 1300 1100 136.137085 1300 1100 1136.137085
        1350 1100 125.2067719 1350 1100 1125.206787 1400 1100 113.8356628 1400 1100 1113.835693 1450 1100 102.4964981 1450 1100 1102.49646


      ZCORN                                  -- Generated : Petrel
        158.552 147.327 147.327 136.137 136.137 125.207 125.207 113.836 113.836 102.496 102.496 90.2352 90.2352 77.6534 77.6534 64.0504
        64.0504 49.6922 49.6922 34.3368 34.3368 18.0418 18.0418 1.0358 1.0358 -16.982 -16.982 -35.0938 -35.0938 -54.4178 -54.4178 -73.3326


      PORO                                   -- Generated : Petrel
      -- Property name in Petrel : Poro_m_3
        43700*0.03 /

      PERMX                                  -- Generated : Petrel
      -- Property name in Petrel : Perm_f_15
        43700*15 /

      PERMY                                  -- Generated : Petrel
      -- Property name in Petrel : Perm_f_15
        43700*15 /

      PERMZ                                  -- Generated : Petrel
      -- Property name in Petrel : Perm_f_15
        43700*15 /

      INITPRESS -- pressure or bottom pressure- depth or top pressure - density or 0. --

      INITTEMP -- bottom temperature - top temperature or -280. - 0. or temperature gradient --


      And that's the output :

        0 1          0.300000000E-010.148038500E-130.148038500E-130.148038500E-13
        0 2          0.300000000E-010.148038500E-130.148038500E-130.148038500E-13
        0 3          0.300000000E-010.148038500E-130.148038500E-130.148038500E-13


      Permeabilities have been written beside porosities on the first line for each element (converted for mD to m2)

  • Thanks, Florian, for this information.

    This is the TOUGHREACT format (4E15.9; iTOUGH2 and TOUGH3 use 5E10.4). Is there a T2B documentation? If so, it should indicate for which TOUGH package this interface was written (specifically that this does not work with TOUGH2).



    • That's fine. The doc says : 

      What is T2B ?

      T2B is a workflow for handling heterogeneous 3D models with the TOUGH2 family of codes. It is composed of separate FORTRAN subroutines (or modules) that can be run independently using input and output files in ASCII format for TOUGH2.

      These modules have been developed specifically for modeling of carbon dioxide geological storage and their use

      • with TOUGH2 and the Equation of State module ECO2N, dedicated to CO2-water-salt mixture systems,
      • with TOUGHREACT, which is an adaptation of TOUGH2 with ECO2N and geochemical fluid-rock interactions
      • and with TOUGH2 and the EOS7C module dedicated to CO2-CH4gas mixture

      is described in "Audigane, P., et al., A workflow for handling heterogeneous 3D models with the TOUGH2 family of codes: Applications to numerical modeling of CO2 geological storage. Computers and Geosciences(2011), doi:10.1016/j.cageo.2010.11.020"


      The objective is to save time for the pre-processing, execution and visualization of complex geometry for geological system representation.


      So Tough2 in mentioned. I guess people using Tough2 are aware that there's no in using permeability grid, and anyway, even if you do and end up with permeabilities in your INCON (like I did), it's just no taken into account. 


      Thanks again for helping. I'll adapt my MESH file, and maybe try ToughReact as well, even if I don't need chemical reactions in my case.

  • Florian DELSANTE


    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.

Like Follow
  • 2 wk agoLast active
  • 18Replies
  • 550Views
  • 5 Following