Case00 Scaling Salt Viscosity

 

In this tutorial page the main concept and approach for scaling salt viscosity in a MEM workflow will be demonstrated. To this end a single axisymmetric element test will be used. A top surface stress is applied and the response of the salt will be monitored via a history point (see figure below).

 

 

MEM_002_013

Schematic of model geometry and boundary conditions

 

 

Salt rocks have a visco-plastic rheological behaviour that is associated with creep driven by deviatoric stress. The creep deformation facilitates deviatoric stress relaxation (decrease in stresses) and consequently the stress within salt rocks is often close to isotropic / hydrostatic. The strain rate at which salt creeps is controlled by salt viscosity, which defines the resistance of a given material to flow (creep) and is expressed in units of [ stress · time ]. Thus for a given viscosity and a given load generating deviatoric stress the salt will be able to creep at a given rate. Consequently, the length of the time frame during which salt creeps will determine the amount of strain and hence the amount of deviatoric stress relaxation. It should be noted that to model salt properly within MEM workflows the targets are:

 

1.To achieve a quasi isotropic stress state within salt during initialisation

2.To capture an appropriate/realistic salt response during production which is elastic with creep at production time scales

 

To achieve the above target #2, appropriate salt properties must be defined in the material model (e.g. properties that may be measured in the lab). However, a quasi isotropic stress state during initialisation may not be achieved with these properties due to the relative small time frame. Initialisation over a larger time frame would facilitate salt creep and deviatoric stress relaxation. However this would require either; too many time steps (resulting in excessive CPU time) or a sudden change in time step size between initialisation and production (which is not recommended to avoid stability issues). Hence the recommended approach is to initially define salt viscosity in time units that will facilitate creep ("dka" which corresponds to 1 simulation time unit = 10000 years), perform the initial initialisation stages and then during subsequent initialisation stages gradually scale the time units to achieve production time units (days).

 

To demonstrate the effect of such approach in salt response 2 single element tests are used. In the fist one the time units will be gradually updated as recommended for MEM initialisation whereas in the second one the same time units will be maintained during the simulation. The files are provided in: MEM_002\Creep_Test\Data

 

The five stages for the test have a duration of 1.0 simulation time units and are defined as follows:

 

1.Load application (1·106 Pa)

2.Creep with initial time units (dKa)

3.Creep with time units of cyears

4.Creep with time units of years

5.Creep with time units of days

 

Note that a case in which all stages occur at dka time units is provided to allow comparison of the results.

 

 

Initial Salt Properties

 

Here the salt properties are discussed. Note that only visco-plastic properties will be discussed.

 

Data File


 

* Material_data   NUM=1                

! ----------------------------------------

 Material_name                  "Salt_Pa_dKa"

! ----------------------------------------  

 Units IDM=4

  /Stress/             "Pa"

  /Length/              "m"

  /Time/              "dka"

  /Temperature/   "Celcius"

! --------------------------------------------------------        

 Grain_stiffness                      70000E6   ! Pa

 Grain_density                        2200.00   ! Kg/m3

! --------------------------------------------------------  

 Porosity_model_type                        1  

 Porosity                                0.02

! --------------------------------------------------------

 Elastic_model_type                         1  

 Elastic_properties  IDM=2

  /Young's Modulus (E)/              18.1E+9   ! Pa

  /Poisson's Ratio (v)/                0.265

! --------------------------------------------------------        

 Viscoplastic_material_type                 1

 Viscoplastic_properties  IDM=3

  /Yield Strength/                    100000   ! Pa

  /Initial K (viscosity)/             3.2E+6  

  /Exponent n/                           1.0

! --------------------------------------------------------        

 Singlephase_fluid_name            "Water"  

 

1.The salt material properties are input in Pa and dka units.  It should be noted that the Units listed here are for reference only and not used by the code.

2.The viscoplastic model used is set to 1 which corresponds to Herschel-Bulkley model (Herschel and Bulkley (1926)).

3.The properties are defined so that:

a.Viscoplastic yield occurs at 105 Pa (0.1 MPa)

b.The viscosity as a function of deviatoric stress is asymptotic to 1·1018 Pa·s at high stresses

4.This requires a viscosity constant K=3.2·106 Pa·dka and an exponent n = 1.

 

MEM_002_012

Viscosity as a function of deviatoric stress for the input properties

 

 

 

 

Change of time units

 

The time scale can be changed during the simulation by:

Definition of initial time units.

Re-specification of Units with Time_update in later stages as required.

 

Base on the initial time units and updated time units defined, the code will apply the appropriate scaling factor to the salt viscosity value input in the material data.

 

Definition Initial Units

 

Data File


 

* Units

! ---------------------------------

 Stress             "Pa"

 Length              "m"

 Time              "dka"

 Temperature   "Celsius"

 

1.Definition of Time in Units is compulsory when Time_update is going to be used to change the time scale during the simulation as it will affect the salt viscosity.

2.The time units defined must be "dka", consistent with the units used in the input value for viscosity

 

 

 

Update of time units

 

Data File


 

* Units

! ---------------------------------

 Time_update       "cyears"

 

* Control_data

! ========================================

 Control_title                   "Test"

 

 

 

 

* Units

! ---------------------------------

 Time_update       "years"        

 

* Control_data

! ========================================

 Control_title                   "Test"

 

 

 

 

* Units

! ---------------------------------

 Time_update       "days"  

 

* Control_data

! ========================================

 Control_title                   "Test"

 

1.In stages 3, 4 and 5 update of the time units is performed so that gradually the time units are changed from dka to days.

 

 

 

 

Results

 

The results for the project are located in the MEM_002\Creep_Test\Results folder, including the spreadsheet 00_MEM_002_Results_Creep_Test.xlsx which contains the history results output from the simulations.

 

The results below show that:

 

1.When the time units are updated at time 2, there is a drastic reduction in further development of viscoplastic strain due to creep.

 

2.As the time units are updated the viscosity increases accordingly.

 

3.The viscosity for both cases is scaled to Pa·s for every output data point and then both cases are compared.  Both cases  show a consistent viscosity as a function of stress which agrees with the input model, thus validating the approach.

 

 

MEM_002_014

Comparison of history results for case with and without update of time units.