Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

被引:132
作者
Doetsch, Joseph [1 ]
Linde, Niklas [2 ]
Coscia, Ilaria [1 ]
Greenhalgh, Stewart A. [1 ,3 ]
Green, Alan G. [1 ]
机构
[1] ETH, Inst Geophys, CH-8093 Zurich, Switzerland
[2] Univ Lausanne, Inst Geophys, CH-1015 Lausanne, Switzerland
[3] Univ Adelaide, Dept Phys, Adelaide, SA 5005, Australia
基金
瑞士国家科学基金会;
关键词
TRAVEL-TIME INVERSION; DC RESISTIVITY; ELECTRICAL-CONDUCTIVITY; P-WAVE; VELOCITY; CONSTRAINTS; TOMOGRAPHY; EQUATIONS; MODELS;
D O I
10.1190/1.3496476
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.
引用
收藏
页码:G53 / G64
页数:12
相关论文
共 42 条
[1]  
[Anonymous], 1998, GSLIB Geostatistical software library and users guide
[2]  
[Anonymous], 1997, MACHINE LEARNING, MCGRAW-HILL SCIENCE/ENGINEERING/MATH
[3]   The electrical resistivity log as an aid in determining some reservoir characteristics [J].
Archie, GE .
TRANSACTIONS OF THE AMERICAN INSTITUTE OF MINING AND METALLURGICAL ENGINEERS, 1942, 146 :54-61
[4]   Seismic reservoir mapping from 3-D AVO in a North Sea turbidite system [J].
Avseth, P ;
Mukerji, T ;
Jorstad, A ;
Mavko, G ;
Veggeland, T .
GEOPHYSICS, 2001, 66 (04) :1157-1176
[5]  
Bedrosian PA, 2007, GEOPHYS J INT, V170, P737, DOI [10.1111/j.1365-246X.2007.03440.x, 10.1111/J.1365-246X.2007.03440.x]
[6]   Cross-hole resistivity tomography using different electrode configurations [J].
Bing, Z ;
Greenhalgh, SA .
GEOPHYSICAL PROSPECTING, 2000, 48 (05) :887-912
[7]  
Bouman C.A., 1997, Cluster: An unsupervised algorithm for modeling Gaussian mixtures
[8]   Cross-property relations between electrical conductivity and the seismic velocity of rocks [J].
Carcione, Jose M. ;
Ursin, Bjorn ;
Nordskag, Janniche I. .
GEOPHYSICS, 2007, 72 (05) :E193-E204
[9]   Assessing the resolution-dependent utility of tomograms for geostatistics [J].
Day-Lewis, FD ;
Lane, JW .
GEOPHYSICAL RESEARCH LETTERS, 2004, 31 (07) :L075031-4
[10]   MAXIMUM LIKELIHOOD FROM INCOMPLETE DATA VIA EM ALGORITHM [J].
DEMPSTER, AP ;
LAIRD, NM ;
RUBIN, DB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-METHODOLOGICAL, 1977, 39 (01) :1-38