WRERAQ 0043-1397 Crossref Web of Science Google Scholar Barrash, 2009, Estimating porosity with ground-penetrating radar reflection tomography: A controlled 3-D experiment at the Boise hydrogeophysical research site: Water Resources Research, 45, W00D26, doi: 10.1029/2008WR006960. Knoll, 1999, Boise hydrogeophysical research site (BHRS): Objectives, design, initial geostatistical results: Symposium on the Application of Geophysics to Engineering and Environmental Problems, 389–398. Clemo, 2002, Hierarchical geostatistics and multifacies systems: Boise hydrogeophysical research site, Boise, Idaho: Water Resources Research, 38, 14-1–14-18, doi: 10.1029/2002WR001436. E., 1942, The electrical resistivity log as an aid in determining some reservoir characteristics: Transactions of the AIME, 146, 54–62, doi: 10.2118/942054-G. First, on a 2.5D synthetic example, and then on field data acquired in a controlled alluvial aquifer, where we were able to match the recovered conductivity to borehole observations. We evaluate two case studies to assess the validity of our algorithm. We physically appraise the domain of the recovered conductivity using a cutoff of the electric current density present in our survey. Moreover, our method linearly outperforms operation counts when compared with 3D Gauss-Newton conjugate-gradient schemes, which scales cubically in our favor with respect to the thickness of the 3D domain. By assuming a 2.5D subsurface, we are able to linearly reduce memory requirements when compared to a 3D gradient descent inversion, and by a power of two when compared to storing a 2D Hessian. Rather, we exploit matrix-vector multiplication of sparse matrices and find successful convergence using gradient descent for our inversion routine without having to resort to the Hessian of the objective function. In doing so, we make no finite-difference approximation on the Jacobian of the data and avoid the need to store large and dense matrices. We have used the adjoint method directly in the discretized Maxwell’s steady state equation to compute the data sensitivity to the conductivity. Our algorithm enables the ER data sensitivities to be directly joined with other geophysical data without the need of interpolating or coarsening the discretization. Due to numerical stability criteria and available computational memory, joint inversion of different types of geophysical data can impose different grid discretization constraints on the model parameters. We have developed a memory and operation-count efficient 2.5D inversion algorithm of electrical resistivity (ER) data that can handle fine discretization domains imposed by other geophysical (e.g, ground penetrating radar or seismic) data.
0 kommentar(er)
