Modeling fluid flow and compaction in sedimentary basins using mixed finite elements

Kikinzon, Evgeny; Kuznetsov, Yuri A.; Maliassov, Serguei; Sumant, Prasad S.

doi:10.1007/s10596-015-9470-2

Cited by 2 publications

(1 citation statement)

References 15 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The coupled system is iteratively solved using the mixed finite element method. Kikinzon et al [17] follows the model set forth in [9], while modifying the algorithm for the timedependent grid structure. They clearly divide the problem into physically distinct steps which are treated in sequence, an operator splitting approach.…”

Section: Introductionmentioning

confidence: 99%

Novel basin modelling concept for simulating deformation from mechanical compaction using level sets

et al. 2017

View full text Add to dashboard Cite

As sedimentation progresses in the formation and evolution of a depositional geologic basin, the rock strata are subject to various stresses. With increasing lithostatic pressure, compressional forces act to compact the porous rock matrix, leading to overpressure buildup, changes in the fluid pore pressure and fluid flow. In the context of petroleum systems modelling, the present study concerns the geometry changes that a compacting basin experiences subject to deposition. The purpose is to track the positions of the rock layer interfaces as compaction occurs. To handle the challenge of potentially large geometry deformations, a new modelling concept is proposed that couples the pore pressure equation with a level set method to determine the movement of lithostratigraphic interfaces. The level set method propagates an interface according to a prescribed speed. The coupling term for the pore pressure and level-set equations consists of this speed function, which is dependent on the compaction law. The two primary features of this approach are the simplicity of the grid and the flexibility of the speed function. A first evaluation of the model concept is presented based on an implementation for one spatial dimension accounting for vertical effective stress. Isothermal conditions with a constant fluid density and viscosity were assumed. The accuracy of the implemented numerical solution for the case of a single stratigraphic unit with a linear compaction law was compared to the available analytical solution [38]. The multi-layer setup and the nonlinear case were tested for plausibility.

show abstract

Section: Introductionmentioning

confidence: 99%