TY - JOUR
T1 - A scalable collocated finite volume scheme for simulation of induced fault slip
AU - Novikov, A.
AU - Voskov, Denis
AU - Khait, Mark
AU - Hajibeygi, Hadi
AU - Jansen, Jan Dirk
PY - 2022
Y1 - 2022
N2 - We present a scalable collocated Finite Volume Method (FVM) to simulate induced seismicity as a result of pore pressure changes. A discrete system is obtained based on a fully-implicit fully-coupled description of flow, elastic deformation, and contact mechanics at fault surfaces on a flexible unstructured mesh. The cell-centered collocated scheme leads to a convenient integration of the different physical equations, as the unknowns share the same discrete locations on the mesh. Additionally, a generic multi-point flux approximation is formulated to treat heterogeneity, anisotropy, and cross-derivative terms for both flow and mechanics equations. The resulting system, though flexible and accurate, can lead to excessive computational costs for field-relevant applications. To resolve this limitation, a scalable processing algorithm is developed and presented. Several proof-of-concept numerical tests, including benchmark studies with analytical solutions, are investigated. It is found that the presented method is indeed accurate and efficient; and provides a promising framework for accurate and efficient simulation of induced seismicity in various geoscientific applications.
AB - We present a scalable collocated Finite Volume Method (FVM) to simulate induced seismicity as a result of pore pressure changes. A discrete system is obtained based on a fully-implicit fully-coupled description of flow, elastic deformation, and contact mechanics at fault surfaces on a flexible unstructured mesh. The cell-centered collocated scheme leads to a convenient integration of the different physical equations, as the unknowns share the same discrete locations on the mesh. Additionally, a generic multi-point flux approximation is formulated to treat heterogeneity, anisotropy, and cross-derivative terms for both flow and mechanics equations. The resulting system, though flexible and accurate, can lead to excessive computational costs for field-relevant applications. To resolve this limitation, a scalable processing algorithm is developed and presented. Several proof-of-concept numerical tests, including benchmark studies with analytical solutions, are investigated. It is found that the presented method is indeed accurate and efficient; and provides a promising framework for accurate and efficient simulation of induced seismicity in various geoscientific applications.
KW - Collocated grid
KW - Contact mechanics
KW - Finite volume
KW - Induced seismicity
KW - Poroelastic
UR - http://www.scopus.com/inward/record.url?scp=85138047091&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2022.111598
DO - 10.1016/j.jcp.2022.111598
M3 - Article
AN - SCOPUS:85138047091
SN - 0021-9991
VL - 469
JO - Journal of Computational Physics
JF - Journal of Computational Physics
M1 - 111598
ER -