You are here:
Time-Lapse Joint Inversion of Cross-Well DC Resistivity and Seismic Data: A Numerical Investigation
Karaoulis, M., A. Revil, J. Zhang, AND D. Werkema. Time-Lapse Joint Inversion of Cross-Well DC Resistivity and Seismic Data: A Numerical Investigation. GEOPHYSICS. Society of Exploration Geophysicists, 77(4):141-157, (2013).
The time-lapse joint inversion of geophysical data is required to solve a number of problems such as the management of oil and gas reservoirs, the sequestration of carbon dioxide, the leakage of water in earth dams and embankments through internal erosion, bioremediation, the production of geothermal reservoirs, and the monitoring of active faults and volcanoes (Lazaratos and Marion, 1997; McKenna et al., 2001; Kowalsky et al., 2006; Ajo-Franklin et al., 2007a, b; Miller et al., 2008; Doetch et al., 2010; Ayeni and Biondi, 2010; Liang et al., 2011).
Time-lapse joint inversion of geophysical data is required to image the evolution of oil reservoirs during production and enhanced oil recovery, CO2 sequestration, geothermal fields during production, and to monitor the evolution of contaminant plumes. Joint inversion schemes reduce space-related artifacts in filtering out noise that is spatially uncorrelated while time lapse inversion algorithms reduce time-related artifacts in filtering out noise that is uncorrelated over time. There are several approaches that are possible to perform the joint inverse problem. In this work, we investigate both the Structural Cross-Gradient (SCG) joint inversion approach and the Cross-Petrophysical (CP) approach, which are both justified for time-lapse problem by petrophysical models. In the first case, the inversion scheme looks for models with structural similarities. In second the case, we use a direct relationship between the geophysical parameters 38. Time-lapse inversion is performed with an actively time-constrained (ATC) approach. In this approach, the subsurface is defined as a space-time model. All the snapshots are inverted together assuming a regularization of the sequence of snapshots over time. First we show the advantage of combining the SCG or CP inversion approaches and the ATC inversion by using a synthetic problem corresponding to cross-hole seismic and DC-resistivity data and piecewise constant resistivity and seismic velocity. We show that the combined SCG/ATC approach reduces the presence of artifacts both with respect to individual inversion of the resistivity and seismic datasets as well as with respect to the joint inversion of both data sets at each time step. We also performed a synthetic study using a secondary oil recovery problem. The combined CP/ATC approach is successful in retrieving the position of the oil/water encroachment front.