Full-waveform inversion of Crosshole GPR data: Implications for porosity estimation in chalk

Abstract The Maastrichtian-Danian chalk is a widely distributed hydrocarbon and groundwater reservoir rock in north-western Europe. Knowledge of lateral and vertical heterogeneity and porosity variation in this type of rock is essential, since they critically determine the reservoir properties. We have collected a densely sampled crosshole ground-penetrating radar (GPR) dataset from a highly heterogeneous section of the chalk and inverted it with a full-waveform inversion (FWI) approach. To date, successful crosshole FWI has only been reported for a handful of GPR field data, none of which include strongly heterogeneous environments like the one considered in this study. Testing different starting models shows that all FWI results converge to very similar subsurface structures indicating that the results are robust with regard to local variations in the permittivity starting models and are not very sensitive to the conductivity starting models. Compared to their ray-based counterparts, the obtained FWI models show significantly higher resolution and improved localization of fine-scale heterogeneity. The final FWI permittivity tomogram was converted to a bulk porosity model using the Complex Refractive Index Model (CRIM) and comparisons with plug sample porosities and televiewer image logs verify that variations in the obtained permittivity are related to facies and lithology changes. The inferred porosity varies from 30 to 54%, which is consistent with values in the chalk cores from the investigated boreholes and in agreement with other studies conducted in similar rocks onshore. Moreover, porosities vary significantly over scales of less than a meter both laterally and vertically. The FWI constrains porosity variation with decimeter scale resolution in our 5 m (horizontally) by 10 m (vertically) model section bridging the gap between what is measured on the core sample scale and the scale typical of hydrogeophysical field experiments conducted to characterize fluid flow in the subsurface. The results provide complementary knowledge to traditional chalk reservoir characterization.

[1]  A. Krogsbøll,et al.  Chalk porosity and sonic velocity versus burial depth: Influence of fluid pressure, hydrocarbons, and mineralogy , 2008 .

[2]  P. Jakobsen,et al.  Infrared Thermography and Fracture Analysis of Preferential Flow in Chalk , 2005 .

[3]  H. Maurer,et al.  Ray-based amplitude tomography for crosshole georadar data: a numerical assessment , 2001 .

[4]  N. Bleistein.  Two-and-One-Half Dimensional In-Plane Wave Propagation. , 1984 .

[5]  Knud Skou Cordua,et al.  Estimation of Chalk Heterogeneity from Stochastic Modeling Conditioned by Crosshole GPR Traveltimes and Log Dat , 2010 .

[6]  S. Greenhalgh,et al.  High resolution imaging of the unsaturated and saturated zones of a gravel aquifer using full-waveform inversion , 2011, 2011 6th International Workshop on Advanced Ground Penetrating Radar (IWAGPR).

[7]  E. Sheldon,et al.  Upper Campanian–Maastrichtian holostratigraphy of the eastern Danish Basin , 2013 .

[8]  Harry Vereecken,et al.  Improved Characterization of Fine-Texture Soils Using On-Ground GPR Full-Waveform Inversion , 2014, IEEE Transactions on Geoscience and Remote Sensing.

[9]  P. Frykman.  Spatial variability in petrophysical properties in Upper Maastrichtian chalk outcrops at Stevns Klint, Denmark , 2001 .

[10]  E. Sheldon,et al.  3D seismic mapping and porosity variation of intra-chalk units in the southern Danish North Sea , 2010 .

[11]  F. Surlyk,et al.  The influence of depositional processes on the porosity of chalk , 2012, Journal of the Geological Society.

[12]  Andrew Binley,et al.  Monitoring Unsaturated Flow and Transport Using Cross‐Borehole Geophysical Methods , 2008 .

[13]  Stéphane Lanteri,et al.  Discontinuous Galerkin frequency domain forward modelling for the inversion of electric permittivity in the 2D case , 2011 .

[14]  I. Fabricius.  Chalk: composition, diagenesis and physical properties , 2007 .

[15]  L. Stemmerik,et al.  Integrated stratigraphy of the late Campanian – Maastrichtian in the Danish Basin: revision of the Boreal calcareous nannofossil zonation , 2016 .

[16]  Karsten H. Jensen,et al.  Monitoring CO2 gas-phase migration in a shallow sand aquifer using cross-borehole ground penetrating radar , 2015 .

[17]  D. Daniels.  Ground Penetrating Radar , 2005 .

[18]  Hansruedi Maurer,et al.  Taming the non-linearity problem in GPR full-waveform inversion for high contrast media , 2011 .

[19]  Jan Vanderborght,et al.  Imaging and characterization of facies heterogeneity in an alluvial aquifer using GPR full-waveform inversion and cone penetration tests , 2015 .

[20]  F. Surlyk,et al.  Late Maastrichtian chalk mounds, Stevns Klint, Denmark — Combined physical and biogenic structures , 2007 .

[21]  Knud Skou Cordua,et al.  Geostatistical inference using crosshole ground-penetrating radar , 2010 .

[22]  Andreas Becht,et al.  Inversion strategy in crosshole radar tomography using information of data subsets , 2004 .

[23]  B. Jacobsen,et al.  Resolution and error propagation analysis for tomographic data with correlated errors , 1996 .

[24]  Michel D'Heur.  Porosity and hydrocarbon distribution in the North Sea chalk reservoirs , 1984 .

[26]  P. Japsen.  Regional Velocity-Depth Anomalies, North Sea Chalk: A Record of Overpressure and Neogene Uplift and Erosion , 1998 .

[27]  Jens Tronicke,et al.  Effects of gas- and water-filled boreholes on the amplitudes of crosshole georadar data as inferred from experimental evidence , 2004 .

[28]  A. Neal.  Ground-penetrating radar and its use in sedimentology: principles, problems and progress , 2004 .

[29]  David L. Alumbaugh,et al.  Estimating moisture contents in the vadose zone using cross‐borehole ground penetrating radar: A study of accuracy and repeatability , 2002 .

[30]  Niklas Linde,et al.  3-D characterization of high-permeability zones in a gravel aquifer using 2-D crosshole GPR full-waveform inversion and waveguide detection , 2013 .

[31]  Florian A. Belina,et al.  Waveform Inversion of Crosshole Georadar Data: Influence of Source Wavelet Variability and the Suitability of a Single Wavelet Assumption , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[32]  T. Hansen,et al.  Monte Carlo full-waveform inversion of crosshole GPR data using multiple-point geostatistical a priori information , 2012 .

[33]  James Irving,et al.  Evaluation of the reconstruction limits of a frequency-independent crosshole georadar waveform inversion scheme in the presence of dispersion , 2012 .

[34]  Harry Vereecken,et al.  Detection of spatially limited high‐porosity layers using crosshole GPR signal analysis and full‐waveform inversion , 2014 .

[35]  G. Mount,et al.  Estimating porosity and solid dielectric permittivity in the Miami Limestone using high‐frequency ground penetrating radar (GPR) measurements at the laboratory scale , 2014 .

[36]  Jacques R. Ernst,et al.  A New Vector Waveform Inversion Algorithm for Simultaneous Updating of Conductivity and Permittivity Parameters From Combination Crosshole/Borehole-to-Surface GPR Data , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[37]  Ludovic Métivier,et al.  Two-dimensional permittivity and conductivity imaging by full waveform inversion of multioffset GPR data: a frequency-domain quasi-Newton approach , 2014 .

[38]  D. Robinson,et al.  The Dielectric Permittivity of Calcite and Arid Zone Soils with Carbonate Minerals , 2004 .

[39]  R. Knight,et al.  Effect of antennas on velocity estimates obtained from crosshole GPR data , 2005 .

[40]  G. McMechan,et al.  Modeling 3D porosity and permeability from GPR data in the Ellenburger Dolomite, central Texas , 2011 .

[41]  Harry Vereecken,et al.  Improvements in crosshole GPR full-waveform inversion and application on data measured at the Boise Hydrogeophysics Research Site , 2013 .

[42]  Jason R. McKenna,et al.  Phase and amplitude inversion of crosswell radar data , 2011 .

[43]  Harry Vereecken,et al.  Optimization of acquisition setup for cross-hole GPR full-waveform inversion using checkerboard analysis , 2013 .

[44]  Knud Skou Cordua,et al.  Accounting for Correlated Data Errors during Inversion of Cross‐Borehole Ground Penetrating Radar Data , 2008 .

[45]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[46]  Jacques R. Ernst,et al.  Application of a new 2D time-domain full-waveform inversion scheme to crosshole radar data , 2007 .

[47]  W. P. Clement,et al.  Crosshole Radar Tomography in a Fluvial Aquifer Near Boise, Idaho , 2006 .

[48]  Harry Vereecken,et al.  Frequency-domain Full-waveform Inversion of GPR Data , 2012 .

[49]  Jeffrey W. Roberts,et al.  Estimation of permeable pathways and water content using tomographic radar data , 1997 .

[50]  A. P. Annan,et al.  Ground-penetrating radar for high-resolution mapping of soil and rock stratigraphy , 1989 .

[51]  Niklas Linde,et al.  Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland , 2010 .

[52]  Seiichiro Kuroda,et al.  Full-waveform inversion algorithm for interpreting crosshole radar data: a theoretical approach , 2007 .

[53]  S. Friedman,et al.  Relationships between the Electrical and Hydrogeological Properties of Rocks and Soils , 2005 .

[54]  Harry Vereecken,et al.  Crosshole GPR full-waveform inversion of waveguides acting as preferential flow paths within aquifer systems , 2012 .

[55]  Andrew Binley,et al.  High‐resolution characterization of vadose zone dynamics using cross‐borehole radar , 2001 .

[56]  M. Erlström,et al.  Chalk depth structure maps, Central to Eastern North Sea, Denmark , 2007 .

[57]  L. Stemmerik,et al.  Production of Calcareous Nannofossil Ooze For Sedimentological Experiments , 2015 .

[58]  I. Fabricius,et al.  How depositional texture and diagenesis control petrophysical and elastic properties of samples from five North Sea chalk fields , 2007, Petroleum Geoscience.

[59]  Alan G. Green,et al.  Results of 3-D georadar surveying and trenching the San Andreas fault near its northern landward limit , 2003 .

[60]  M. Engkilde,et al.  Deep channels in the Cenomanian–Danian Chalk Group of the German North Sea sector: Evidence of strong constructional and erosional bottom currents and effect on reservoir quality distribution , 2008 .

[61]  Jacques R. Ernst,et al.  Full-Waveform Inversion of Crosshole Radar Data Based on 2-D Finite-Difference Time-Domain Solutions of Maxwell's Equations , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[62]  James Irving,et al.  Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole georadar data , 2012 .

[63]  J. J. Peterson.  Pre-inversion Corrections and Analysis of Radar Tomographic Data , 2001 .

[64]  L. Stemmerik,et al.  Diagenesis of Flint and Porcellanite in the Maastrichtian Chalk at Stevns Klint, Denmark , 2010 .