ResInsight computes several derived results. In this section we will explain what they are, and briefly how they are calculated.
ResInsight calculates several derived cell properties that is made available as Static or Dynamic cell properties. The derived results listed at the bottom of the Static result properties, are shown below.
The transmissibility for cells and Non-Neighbor Connections (NNCs) are dependent on both cell properties and geometry. ResInsight normalizes TRANX, TRANY and TRANZ with the overlapping flow area for both neighbor cells and NNC-cells. The results are named riTRANXbyArea, riTRANYbyArea and riTRANZbyArea respectively.
The normalized transmissibilities make it easier to compare and check the flow capacity visually. This can be useful when history matching pressure differences across a fault.
Transmissibility can be set or adjusted with multiple keywords in an Eclipse data deck. To visualize the adjustments made, ResInsight calculates a multiplicator for the overall change. First unadjusted transmissibilities for all neighbor cells and NNCs are evaluated based on geometry and permeabilities, similar to the NEWTRAN algorithm in Eclipse. For x- and y-directions, the NTG parameter is also included. The results are named riTRANX, riTRANY and riTRANZ respectively.
The TRANX, TRANY and TRANZ used in the simulation are divided by the ResInsight calculated transmissibilities and the resulting multiplicators are named riMULTX, riMULTY and riMULTZ respectively. The derived properties are listed under Static properties. The riMULT-properties are useful for quality checking consistence in user input for fault seal along a fault plane.
Cell properties with names ending in I, J, K, X, Y, or Z, and an optional “+” or “-” are combined into derived results post-fixed with IJK, or XYZ depending on their origin. (Eg. the static cell properties MULTX, MULTY, MULTZ, and their negatives are combined into the result MULTXYZ, while the dynamic cell properties FLRGASI, FLRGASJ, FLRGASK are combined to FLRGASIJK).
These combined cell properties visualize the property as a color in all directions combined when selected in as a Cell Result and Separate Fault Result.
The face of a cell is then colored based on the value associated with that particular face. The Positive I-face of the cell gets the cell X/I-value, while the J-face gets the Y/J-value etc. The negative faces, however, get the value from the neighbor cell on that side. The negative I-face gets the X-value of the IJK-neighbor in negative I direction, and so on for the J- and K-faces.
The directional combined parameters available are:
The dynamic cell property named Completion Type is calculated from the intersections between Completions and the grid cells. All grid cells intersected by a completion will be assigned a color based on the type of completion that intersects the cell.
If a cell is completed with multiple completions, the following priority is used : Fracture, Fishbones, and Perforation Interval.
In the process of normalizing transmissibility by the overlapping flow area, the NNCs in the model without any shared surface between two cells are identified. These NNCs are listed in the Faults/NNCs With No Common Area folder. These NNCs are questionable since flow normally is associated with a flow area.
Water Flooded PV, also called Number of flooded porevolumes shows the amount of flow from a selected set of simulation tracers into a particular cell, compared to the cells mobile pore volume. A value of 1.0 will tell that the tracers accumulated flow into the cell has reached a volume equal to the mobile pore volume in the cell.
ResInsight calculates several of the presented geomechanical results based on the native results present in the odb-files.
ResInsight can calculate and display relative results, sometimes also referred to as Time Lapse results. When enabled, every result variable is calculated as :
Value’(t) = Value(t) - Value(BaseTime)
Enable the Enable Relative Result option in the Relative Result Options group, and select the appropriate Base Time Step.
Each variable is then post-fixed with “_DTimeStepIndex” to distinguish them from the native variables.
Note: Relative Results calculated based on Gamma values are calculated slightly differently:
Gamma_Dn = ST_Dn / POR_Dn
The calculated result fields are:
In this text the label Sa and Ea will be used to denote the unchanged stress and strain tensor respectively from the odb file.
Components with one subscript denotes the principal values 1, 2, and 3 which refers to the maximum, middle, and minimum principals respectively.
Components with two subscripts however, refers to the global directions 1, 2, and 3 which corresponds to X, Y, and Z and thus also easting, northing, and depth.
Two constants can be assigned to a Geomechanical case:
In the following they are denoted s0 and fa respectively. Some of the derived results use these constants, that can be changed in the property panel of the Case.
Compaction is the difference in vertical displacement (U3) between a grid node and a specified reference K layer. The reference K layer is specified in the property editor.
For each node n in the grid, a node nref in the reference K layer is found by vertical intersection from the node n.
If Depthn <= Depthnref:
COMPACTIONn = -(U3n - U3nref)
COMPACTIONn = -(U3nref - U3n)
STii = -Saii + POR (i= 1,2,3)
STij = -Saij (i,j = 1,2,3 and i not equal j)
We use a value of POR=0.0 where it is not defined.
STi = Principal value i of ST
STM = (ST11 + ST22 + ST33)/3
Q = sqrt( (3⁄2) * ( (ST1 – STM)2 + (ST2 – STM)2 + (ST3 – STM)2 ))
Gammaii = STii/POR (i= 1,2,3)
Gammai = STi/POR
In these calculations we set Gamma to undefined if abs(POR) > 0.01 MPa.
SEij = -Saij (Where POR is defined)
SEij = undefined (Were POR is not defined)
SEi = Principal value i of SE
SEM = (SE11 + SE22 + SE33)/3
SFI = ( (s0/tan(fa) + 0.5*(SE1 + SE3))sin(fa) ) /(0.5(SE1-SE3) )
DSM = tan(rho)/tan(fa)
rho = 2 * (arctan (sqrt (( SE1 + a)/(SE3 + a)) ) – pi/4)
a = s0/tan(fa)
FOS = 1/DSM
Eij = -Eaij
EV = E11 + E22 + E33
ED = 2*(E1-E3)/3
For each face displayed, (might be an element face or an intersection/intersection box face), a coordinate system is established such that:
The stress tensors in that particular face are then transformed to that coordinate system. The following quantities are derived from the transformed tensor named TS in the following:
SN = TS33
TPH = TS31 = TSZX
TPQV = TS32 = TSZY
TP = sqrt(TPH2 + TPQV2)
Angle of the total in-plane shear relative to the Quasi Vertical direction
TPinc = acos(TPQV/TP)
These are the directional angles of the face-normal itself.