TY - JOUR

T1 - Transient-state instability analysis of dissolution-timescale reactive infiltration in fluid-saturated porous rocks

T2 - Purely mathematical approach

AU - Zhao, Chong Bin

AU - Hobbs, Bruce

AU - Ord, Alison

PY - 2019/4/16

Y1 - 2019/4/16

N2 - This paper deals with how the purely mathematical approach can be used to solve transient-state instability problems of dissolution-timescale reactive infiltration (DTRI) in fluid-saturated porous rocks. Three key steps involved in such an approach are: (1) to mathematically derive an analytical solution (known as the base solution or conventional solution) for a quasi-steady state problem of the dissolution timescale, which is viewed as a frozen state of the original transient-state instability problem; (2) to mathematically deduce a group of first-order perturbation partial-differential equations (PDEs) for the quasi-steady state problem; (3) to mathematically derive an analytical solution (known as the perturbation solution or unconventional solution) for this group of first-order perturbation PDEs. Because of difficulty in mathematically solving a transient-state instability problem of DTRI in general cases, only a special case, in which some nonlinear coupling between governing PDEs of the problem can be decoupled, is considered to illustrate these three key steps in this study. The related theoretical results demonstrated that the transient chemical dissolution front can become unstable in the DTRI system of large Zh numbers when the long wavelength perturbations are applied to the system. This new finding may lay the theoretical foundation for developing innovative technique to exploit shale gas resources in the deep Earth.

AB - This paper deals with how the purely mathematical approach can be used to solve transient-state instability problems of dissolution-timescale reactive infiltration (DTRI) in fluid-saturated porous rocks. Three key steps involved in such an approach are: (1) to mathematically derive an analytical solution (known as the base solution or conventional solution) for a quasi-steady state problem of the dissolution timescale, which is viewed as a frozen state of the original transient-state instability problem; (2) to mathematically deduce a group of first-order perturbation partial-differential equations (PDEs) for the quasi-steady state problem; (3) to mathematically derive an analytical solution (known as the perturbation solution or unconventional solution) for this group of first-order perturbation PDEs. Because of difficulty in mathematically solving a transient-state instability problem of DTRI in general cases, only a special case, in which some nonlinear coupling between governing PDEs of the problem can be decoupled, is considered to illustrate these three key steps in this study. The related theoretical results demonstrated that the transient chemical dissolution front can become unstable in the DTRI system of large Zh numbers when the long wavelength perturbations are applied to the system. This new finding may lay the theoretical foundation for developing innovative technique to exploit shale gas resources in the deep Earth.

KW - chemical dissolution

KW - dissolution timescale

KW - mathematical study

KW - perturbation solution

KW - porous rocks

KW - reactive infiltration instability

UR - http://www.scopus.com/inward/record.url?scp=85064754761&partnerID=8YFLogxK

U2 - 10.1007/s11431-018-9448-8

DO - 10.1007/s11431-018-9448-8

M3 - Article

AN - SCOPUS:85064754761

JO - Science in China Series E Engineering and Materials Science

JF - Science in China Series E Engineering and Materials Science

SN - 1674-7321

ER -