Poroelastic model-based time-lapse modeling of the Quest carbon storage project in Alberta

A finite-difference algorithm was developed based on the Biot’s equations of motion for modelling wave propagation in poroelastic media. In contrast with the elastic modelling, in the poroelastic approach the properties of the pore fluid are taken into account in the algorithm. Poroelastic modelling could be useful in cases where the fluid content of the rock is of interest, i.e. Carbon Capture and Storage (CCS) projects. We examined our program using a model based on the Quest CCS project in Alberta to investigate the detectability of CO2 after one year of injection. This was done by defining two models for the baseline and monitor scenarios that represented the subsurface before and after injecting CO2, respectively. The difference between the calculated seismic sections for the two scenarios shows that the residual amplitude is comparable with the signal amplitude. With this result, the injected CO2 in the Quest project over a year could be detected providing the data have good bandwidth and a high signal-to-noise ratio. Furthermore, a comparison between the poroelastic algorithm and the elastic algorithm shows that the time-lapse effect in the poroelastic case is smaller than the one in the elastic case. In the fluid saturated media some of the wave energy dissipates due to the fluid induced flow, and the poroelastic approach helps us to take this loss into account in the modeling process.

[1]  A. Levander Fourth-order finite-difference P-SV seismograms , 1988 .

[2]  J. Carcione,et al.  Computational poroelasticity — A review , 2010 .

[3]  Q. H. Liu,et al.  A staggered-grid finite-difference method with perfectly matched layers for poroelastic wave equations. , 2001, The Journal of the Acoustical Society of America.

[4]  Xiuming Wang,et al.  Modelling Seismic Wave Propagation in Heterogeneous Poroelastic Media Using a High‐Order Staggered Finite‐Difference Method , 2003 .

[5]  José M. Carcione,et al.  SOME ASPECTS OF THE PHYSICS AND NUMERICAL MODELING OF BIOT COMPRESSIONAL WAVES , 1995 .

[6]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[7]  Nanxun Dai,et al.  Wave propagation in heterogeneous, porous media: A velocity‐stress, finite‐difference method , 1995 .

[8]  M. Biot MECHANICS OF DEFORMATION AND ACOUSTIC PROPAGATION IN POROUS MEDIA , 1962 .

[9]  Kagan Tuncay,et al.  Parallel implementation of a velocity-stress staggered-grid finite-difference method for 2-D poroelastic wave propagation , 2006, Comput. Geosci..

[10]  Velocity-stress finite-difference modeling of poroelastic wave propagation , 2013 .

[11]  F. Gaßmann Uber die Elastizitat poroser Medien. , 1961 .

[12]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[13]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[14]  L. Bartel,et al.  Velocity‐stress‐pressure algorithm for 3D poroelastic wave propagation , 2004 .

[15]  G. McMechan,et al.  Numerical simulation of seismic responses of poroelastic reservoirs using Biot theory , 1991 .