Application of different inverse methods for combination of vS and vGPR data to estimate porosity and water saturation


1 M.Sc. Graduate, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran

2 Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran

3 Assistant Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran


Inverse problem is one of the most important problems in geophysics as model parameters can be estimated from the measured data directly using inverse techniques. In this paper, applying different inverse methods on integration of S-wave and GPR velocities are investigated for estimation of porosity and water saturation. A combination of linear and nonlinear inverse problems are solved. Linear least-squares and conjugate gradient are used as linear techniques, whereas grid search and Newton methods are selected as nonlinear ones. It is understood that vS depends on density and Lame Constant (shear modulus) and vGPR on dielectric constant. This combination seems to be logical. Shear modulus is related to porosity using Bruggeman’s rule. Density and dielectric constant is also related to porosity and water saturation. This implies that vS and vGPR are bivariate functions of porosity and water saturation, which are our unknown model parameters. The model parameters are estimated to minimize the cost functional ora system of the equations. In order to convert the nonlinear problem into the linear form, taking logarithm and changing variables were used. The problem was convex, which was inferred from the linear form, so there was just one local minimum as the global minimum of the problem. The grid search method shows that porosity and water saturation cannot be estimated by vGPR or vS uniquely. The results of the four methods were compared with each other and a good agreement was observed.


Main Subjects

Bachrach, R., 2006, Joint porosity and saturation using stochastic rock- physics modeling, Geophysics, 71, 53-63.
Birchak, J. R., Gardner, L. G., Hipp, J. W. and Victor, J. M., 1974, High electric constant microwave probes for sensing soil moisture, Proc. IEEE, 62(1), 93-98.
Dannowski, G. and Yaramanci, U., 1999, Estimation of water content and porosity using combined radar and geoelectrical measurements, European Journal of Environmental and Engineering Geophysics, 4, 71-85.
Boulanger, O. and Chouteau, M., 2001, Constraints in 3D gravity inversion, Geophysical Prospecting, 49, 265-280.
Aster, R. C., Borchers, B. and Clifford, H. T., 2005, Parameter estimation and inverse problems, Elsevier Academic Press.
Ghose, R. and Slob, E. C., 2006, Quantitative integration of seismic and GPR reflections to derive unique estimates for water saturation and porosity in subsoil, Geophysical Research Letters, 33, L05404.
Menke, W., 2012, Geophysical data analysis, discrete inverse theory, Elsevier Academic Press.
Sihvola, A., 1999, Electromagnetic mixing formulas and applications, 284 pp., Inst. of Electr. and Electron. Eng., New York.
Wang, J. R. and Schmugge, T. J., 1980, An empirical model for the complex dielectric permittivity of soils as a function of water content, IEEE Trans. Geosci. Remote Sens., 18, 288-295.