<p>This study employs a lattice Boltzmann model to investigate CO<InlineEquation ID="IEq1"> <EquationSource Format="TEX">$_2$</EquationSource> <EquationSource Format="MATHML"><math> <mmultiscripts> <mrow /> <mn>2</mn> <mrow /> </mmultiscripts> </math></EquationSource> </InlineEquation> reactive transport in fractured porous media, focusing on the effects of Reynolds number (Re), porosity (<InlineEquation ID="IEq2"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation>), and fracture geometry. Results indicate that the steady-state average reaction rate exhibits coupled dependencies on these parameters. At low Re, the steady-state average reaction rate increases monotonically with <InlineEquation ID="IEq3"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation>. In contrast, a non-monotonic relationship is observed at high Re, governed by a critical porosity threshold (<InlineEquation ID="IEq4"> <EquationSource Format="TEX">$\varepsilon _c$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>ε</mi> <mi>c</mi> </msub> </math></EquationSource> </InlineEquation>). For <InlineEquation ID="IEq5"> <EquationSource Format="TEX">$\varepsilon < \varepsilon _c$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>ε</mi> <mo><</mo> <msub> <mi>ε</mi> <mi>c</mi> </msub> </mrow> </math></EquationSource> </InlineEquation>, the steady-state reaction rate increases with <InlineEquation ID="IEq6"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation>, while it remains constant or decreases for <InlineEquation ID="IEq7"> <EquationSource Format="TEX">$\varepsilon > \varepsilon _c$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>ε</mi> <mo>></mo> <msub> <mi>ε</mi> <mi>c</mi> </msub> </mrow> </math></EquationSource> </InlineEquation>. The critical porosity <InlineEquation ID="IEq8"> <EquationSource Format="TEX">$\varepsilon _c$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>ε</mi> <mi>c</mi> </msub> </math></EquationSource> </InlineEquation> decreases initially with Re before stabilizing. Similarly, the reaction rate’s dependence on Re is modulated by porosity. At low <InlineEquation ID="IEq9"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation>, the reaction rate rises with Re. At high <InlineEquation ID="IEq10"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation>, a critical Reynolds number (<InlineEquation ID="IEq11"> <EquationSource Format="TEX">$\textrm{Re}_\textrm{c}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mtext>Re</mtext> <mtext>c</mtext> </msub> </math></EquationSource> </InlineEquation>) emerges. Below <InlineEquation ID="IEq12"> <EquationSource Format="TEX">$\textrm{Re}_\textrm{c}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mtext>Re</mtext> <mtext>c</mtext> </msub> </math></EquationSource> </InlineEquation>, the steady-state reaction rate increases with Re, but it remains invariant when Re exceeds <InlineEquation ID="IEq13"> <EquationSource Format="TEX">$\textrm{Re}_\textrm{c}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mtext>Re</mtext> <mtext>c</mtext> </msub> </math></EquationSource> </InlineEquation>. The value of <InlineEquation ID="IEq14"> <EquationSource Format="TEX">$\textrm{Re}_c$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mtext>Re</mtext> <mi>c</mi> </msub> </math></EquationSource> </InlineEquation> decreases and then stabilizes as <InlineEquation ID="IEq15"> <EquationSource Format="TEX">$\varepsilon $</EquationSource> <EquationSource Format="MATHML"><math> <mi>ε</mi> </math></EquationSource> </InlineEquation> increases. Furthermore, the reactive transport process is governed by two distinct regimes: diffusion-dominated and convection-dominated regimes. For the diffusion-dominated regime, variations in fracture geometry parameters (including fracture width, inclination angle, number, and cross configuration) exert negligible influence on reactive transport. Under the convection-dominated regime, however, increasing fracture width enhances reactant transport toward the downstream side of the porous medium and intensifies the spatial heterogeneity of the dissolution rate. Among the geometric factors, horizontal fractures exert the strongest influence on CO<InlineEquation ID="IEq16"> <EquationSource Format="TEX">$_2$</EquationSource> <EquationSource Format="MATHML"><math> <mmultiscripts> <mrow /> <mn>2</mn> <mrow /> </mmultiscripts> </math></EquationSource> </InlineEquation> transport and generate the most pronounced dissolution heterogeneity, whereas vertical fractures have the weakest effect. An increase in fracture number correlates with a reduction in the dissolution rate. The effects of cross-fracture configurations on transport are highly similar to those of a single fracture.</p>

Mesoscopic Numerical Study of CO₂ Reactive Transport in Fractured Porous Rock Formations

期刊论文
发表日期：

2026年6月23日

查看全文

Jizong Duan,
Qin Lou

摘要

This study employs a lattice Boltzmann model to investigate CO $_2$ $_{2}$ reactive transport in fractured porous media, focusing on the effects of Reynolds number (Re), porosity ( $\varepsilon $ $ε$ ), and fracture geometry. Results indicate that the steady-state average reaction rate exhibits coupled dependencies on these parameters. At low Re, the steady-state average reaction rate increases monotonically with $\varepsilon $ $ε$ . In contrast, a non-monotonic relationship is observed at high Re, governed by a critical porosity threshold ( $\varepsilon _c$ $ε_{c}$ ). For $\varepsilon < \varepsilon _c$ $ε < ε_{c}$ , the steady-state reaction rate increases with $\varepsilon $ $ε$ , while it remains constant or decreases for $\varepsilon > \varepsilon _c$ $ε > ε_{c}$ . The critical porosity $\varepsilon _c$ $ε_{c}$ decreases initially with Re before stabilizing. Similarly, the reaction rate’s dependence on Re is modulated by porosity. At low $\varepsilon $ $ε$ , the reaction rate rises with Re. At high $\varepsilon $ $ε$ , a critical Reynolds number ( $\textrm{Re}_\textrm{c}$ ${Re}_{c}$ ) emerges. Below $\textrm{Re}_\textrm{c}$ ${Re}_{c}$ , the steady-state reaction rate increases with Re, but it remains invariant when Re exceeds $\textrm{Re}_\textrm{c}$ ${Re}_{c}$ . The value of $\textrm{Re}_c$ ${Re}_{c}$ decreases and then stabilizes as $\varepsilon $ $ε$ increases. Furthermore, the reactive transport process is governed by two distinct regimes: diffusion-dominated and convection-dominated regimes. For the diffusion-dominated regime, variations in fracture geometry parameters (including fracture width, inclination angle, number, and cross configuration) exert negligible influence on reactive transport. Under the convection-dominated regime, however, increasing fracture width enhances reactant transport toward the downstream side of the porous medium and intensifies the spatial heterogeneity of the dissolution rate. Among the geometric factors, horizontal fractures exert the strongest influence on CO $_2$ $_{2}$ transport and generate the most pronounced dissolution heterogeneity, whereas vertical fractures have the weakest effect. An increase in fracture number correlates with a reduction in the dissolution rate. The effects of cross-fracture configurations on transport are highly similar to those of a single fracture.

Mesoscopic Numerical Study of CO2 Reactive Transport in Fractured Porous Rock Formations

摘要

Mesoscopic Numerical Study of CO₂ Reactive Transport in Fractured Porous Rock Formations