Title :
Simulation of porosity field using wavelet bayesian inversion of crosswell GPR and log data
Author :
Gloaguen, Erwan ; Dubreuil-Boisclair, Camille ; Simard, Patrick ; Giroux, Bernard ; Marcotte, Denis
Author_Institution :
Inst. Nat. de la Rech. Sci., Quebec City, QC, Canada
Abstract :
In this paper, we present a novel approach to simulate porosity fields constrained by borehole radar tomography images. The cornerstone of the method is the bayesian analysis of the approximation wavelet coefficients of a petro-physical analogue. The method is tested with a two-dimensional porosity field generated from a digital picture of a real sand deposit. The porosity field is translated into electrical properties and a cross-hole tomography synthetic survey is modeled using a finite-difference modeling algorithm. In parallel, an analogue deposit is created based on the geological knowledge of the area under study. The analogue porosity field is converted into electrical property fields using the same equations as previously. A synthetic GPR tomography is also computed from the latter. Wavelet decomposition of both measured and analogue tomograms and porosity analogue fields is then calculated. Based on the assumption that geophysical data carry only the large-scale information about the geological model, statistical analysis of the approximation coefficients of each variable is carried out. From the measured tomogram approximation coefficients and the cross statistics evaluated on the analogues, the approximation of the real porosity field is inferred using bayesian inference. Finally, based on the geostatistical relationships between wavelet coefficients across the different scales, all the porosity wavelet detail coefficients are simulated using a standard geostatistical simulation algorithm. The wavelet coefficients are then back transformed in the porosity space. The final simulated porosity fields contain the large wavelengths of the measured radar tomogram and the texture of the analogue porosity field.
Keywords :
Bayes methods; finite difference methods; geophysical techniques; sand; statistical analysis; tomography; 2D porosity held; analogue porosity field; analogue tomograms; approximation wavelet coefficients; bayesian analysis; borehole radar tomography images; cross-hole tomography synthetic survey; crosswell GPR; digital picture; electrical properties; electrical property fields; finite-difference modeling algorithm; geological model; log data; petro-physical analogue; porosity analogue fields; porosity space; porosity wavelet detail coefficients; radar tomogram; sand deposit; standard geostatistical simulation algorithm; statistical analysis; synthetic GPR tomography; tomogram approximation coefficients; wavelet Bayesian inversion; wavelet decomposition; Bayesian methods; Geologic measurements; Geology; Geophysical measurements; Geophysics computing; Ground penetrating radar; Radar imaging; Testing; Tomography; Wavelet coefficients;
Conference_Titel :
Ground Penetrating Radar (GPR), 2010 13th International Conference on
Conference_Location :
Lecce
Print_ISBN :
978-1-4244-4604-9
Electronic_ISBN :
978-1-4244-4605-6
DOI :
10.1109/ICGPR.2010.5550181