reactive transport modeling of carbonate …...reactive transport modeling of carbonate diagenesis...
TRANSCRIPT
![Page 1: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/1.jpg)
Reactive transport modeling of carbonate diagenesis on
unstructured grids
Alina Yapparova
Chair of Reservoir Engineering
Montanuniversität Leoben
![Page 2: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/2.jpg)
2
Dolomitization and geochemical modeling
“Hollywood may never make a movie about geochemical
modeling, but the field has its roots in top-secret efforts to formulate rocket fuels in the
1940s and 1950s. Anyone who reads cheap novels knows that
these efforts involved brilliant scientists endangered by spies, counterspies,hidden microfilm, and beautiful but treacherous
women.”
(Bethke, C. Geochemical and Biogeochemical Reaction Modeling,
Cambridge University Press 2008)
![Page 3: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/3.jpg)
3
Agenda
● Method description. Software design● Benchmark 1d: calcite dissolution –
dolomite precipitation ● Calcite dolomitization on a 2d cross-
section with realistic geology
![Page 4: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/4.jpg)
4
Chemical modeling
TOUGH, PHREEQC GEMS
Equilibrium constant methodMinimization of Gibbs
free energy
Two methods are computationally and conceptually equivalent (Zeleznik and Gordon (1960, 1968) and Brinkley (1960) )
Variation in free energy G with reaction progress for the reaction
bB+cC↔dD+eE.
The reaction’s equilibrium point is the minimum along the free energy curve
![Page 5: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/5.jpg)
5
Chemical reactions and Math
“Fast” reaction(aqueous complexation)
“Slow” reaction(dissolution-precipitation)
takes place at a significantly greater rate than the transport processes that redistribute mass
takes place at a significantly smaller rate than the transport processes that redistribute mass, requires a kinetic description because the flow system can remove products and reactants before reactions can proceed to equilibrium
A+B↔C
K eq=[ A ][B]
[C ]
[C ]=1
K eq
[ A ][B]d [C ]
dt=−r ([A ] , [B ] ,T , SA ,SI ...)
Algebraic system of nonlinear equations
System of ordinary differential equations (ODE)
![Page 6: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/6.jpg)
6
Simple carbonate system
= =
=
One reaction for each mineral and each secondary species
Secondary species concentrations are expressed using concentrations of primary species and equilibrium constants
![Page 7: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/7.jpg)
7
Mass conservation
For each primary species total concentration in mobile and immobile components remains constant
– mineral concentrations for calcite and dolomite,
![Page 8: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/8.jpg)
8
Chemical system of equations. Change of variables
- initial total aqueous concentrations of primary species
- initial mineral amounts
![Page 9: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/9.jpg)
9
Reactive transport modeling
ddt
T+L(C )=q
Φ(X ,T )=0C−C (X )=0
- transport equations written for each primary species,
- mass conservation equations
T – total concentration of chemical species (mobile and immobile), C – mobile concentration of chemical species, L – transport operator (advection-diffusion-dispersion), q – source term
Sequential non-iterative approach (SNIA)
Sequential iterative approach (SIA)
T n+1+Δ t L(Cn)−T n=Δ t qn
Φ(X ,T n+1)=0Cn+1−C (X )=0
T n+1+Δ t L(Cn+1
)−T n=Δ t qn
Φ(X ,T n+1)=0Cn+1−C (X )=0
(after C. de Dieuleveult, J. Erhel, M. Kern A global strategy for solving reactive transport equations, Journal of Computational Physics, 2009)
![Page 10: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/10.jpg)
10
Solution procedure
Initial equilibration
time looptn+1
=tn+Δt
Solve transport equations
Solve chemical equilibrium equations
tn+1
≥tmax
- for all finite volumes- for each basic species separately
- for each finite volume separately- for all species
yes
no
Transport equations are solved by operator splitting using FEFV method implemented in CSMP++
System of nonlinear algebraic equations for chemical equilibrium is solved using Newton-Raphson method using KINSOL library
![Page 11: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/11.jpg)
11
Numerical issues. Oscillations
Initial guess for Newton-Raphson:
● Always choose the most abundant species as primary species
● 95% of total species aqueous concentration is a good initial approximation for [Ca2+] and [Mg2+]
● For [CO3
2-] initial guess from the
previous time step
![Page 12: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/12.jpg)
12
Mineral swap procedure
After equilibration some mineral concentrations
are negative
Exclude the most negative mineral
Recalculate equilibrium
Yes
NoNext time step
![Page 13: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/13.jpg)
13
Agenda
● Method description. Software design● Benchmark 1d: calcite dissolution –
dolomite precipitation ● Calcite dolomitization on a 2d cross-
section with realistic geology
![Page 14: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/14.jpg)
14
Benchmark: calcite dissolution, dolomite precipitation
Parameter Value Unit
Column length
0.5 m
Porosity 0.32 -
Bulk density 1800 kg/m3
Pore velocity 9.375x10-6 m/s
Flow rate 3x10-6 m/s
Material properties
0.000122M Calcite0.001MMgCl
2
0.0m 0.5mModel domain
P. Engesgaard, K., Kipp, A geochemical transport model for redox-controlled movement of mineral fronts in groundwater flow systems: A case of nitrate removal by oxidation of pyrite, Water Resources Research, 28, pp. 2829-2843, 1992.
● Pore fluid in 1D column is initially equilibrated with calcite
● Column is flushed with MgCl2 solution
● Calcite dissolves, dolomite is formed temporarily as a moving zone
temperature 25oCpressure 1atmmesh size 0.01mtime step 200s
![Page 15: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/15.jpg)
15
Benchmark: calcite dissolution, dolomite precipitation
P. Engesgaard, K., Kipp, A geochemical transport model for redox-controlled movement of mineral fronts in groundwater flow systems: A case of nitrate removal by oxidation of pyrite, Water Resources Research, 28, pp. 2829-2843, 1992.
Boundary Initial
pH 7.06 9.91
Ca2+ 0.0 1.239e-4
CO32- 0.0 1.239e-4
Mg2+ 1.e-3 0.0
Cl- 2.e-3 0.0
CaCO3(s) 0.0 2.17e-5
CaMg(CO3)
20.0 0.0
CaCO3(s) = Ca2+ + CO
32- (-8.47)
CaMg(CO3)
2= Ca2+ + Mg2+ + 2CO
32- (-17.17)
Ca2+ + CO3
2- = CaCO3(aq) (3.23)
Mg2+ + CO3
2- = MgCO3(aq) (2.98)
H+ + OH- = H2O (-14.01)
H+ + CO3
2- = HCO3
(10.31)
2H+ + CO3
2- = H2O + CO
2(aq) (16.71)
Cl- = Cl-
Component and solid concentrations
![Page 16: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/16.jpg)
16
Benchmark: calcite dissolution, dolomite precipitation
![Page 17: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/17.jpg)
17
Benchmark: calcite dissolution, dolomite precipitation
(Engesgaard, Kipp, 1992)
Cl-
Mg2+
Ca2+
time = 21000s
Calcite
Dolomite
![Page 18: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/18.jpg)
18
Agenda
● Method description. Software design● Benchmark 1d: calcite dissolution –
dolomite precipitation ● Calcite dolomitization on a 2d cross-
section with realistic geology
![Page 19: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/19.jpg)
19
2D simulation:geological cross-section
![Page 20: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/20.jpg)
20
2D simulation: rock properties
![Page 21: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/21.jpg)
21
2D simulation: essential conditions
Initial conditions: water equilibrated with calciteBoundary conditions:
LEFT – constant inflow of MaCl2 solution,
RIGHT – constant pressure
![Page 22: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/22.jpg)
22
2D simulation: results
![Page 23: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/23.jpg)
23
Conclusions: 2D simulation
● Bioherms get partially dolomitized in the direction of fluid flow
● Calcite sands get partially dolomitized due to Mg provided from permeable sandstone layers
● Calcite mudstone layer gets partially dolomitized after calcite sand layer gets completely dolomitized in the middle part
![Page 24: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/24.jpg)
24
Conclusions: general
● RTM prototype module was implemented in CSMP
● 1D benchmark was conducted and results are in a good agreement
● 2D simulation was performed and the ability of combining chemistry with transport with realistic geometry was shown
● General framework for further software development was established
![Page 25: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/25.jpg)
25
Thank you for your attention!
![Page 26: Reactive transport modeling of carbonate …...Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität](https://reader033.vdocument.in/reader033/viewer/2022052406/5f0849de7e708231d421450f/html5/thumbnails/26.jpg)
26
Future work
● Include more minerals and aqueous species
● Pitzer activity model● Kinetic control reactions● New porosity/permeability correlation● 3d simulations