• Economic and Environmental Geology
• /
• v.14 no.3
• /
• pp.193-201
• /
• 1981

The problem of automatic inversion of apparent resistivity sounding curves resulting from horizontally layered earth models is solved using the least-squares technique. This method, which makes use of damped least-squares algorithm in conjunction with digital filtering technique, is found to be speedier and more accurate than the conventional curve-matching method. Four sounding curves were chosen to test the inversion scheme. The analysis of the theoretical sounding data associated with a three-layer model illustrates clear advantages over the conventional curve-matching method. The usefulness of the inversion method is also shown when applied to the actual field data. It was found that the best fit earth models coincide with the subsurface structures confirmed by drilling.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.129-136
• /
• 2005

The damped least-squares inversion has become a most popular method in finding the solution in geophysical problems. Generally, the least-squares inversion is to minimize the object function which consists of data misfits and model constraints. Although both the data misfit and the model constraint take an important part in the least-squares inversion, most of the studies are concentrated on what kind of model constraint is imposed and how to select an optimum regularization parameter. Despite that each datum is recommended to be weighted according to its uncertainty or error in the data acquisition, the uncertainty is usually not available. Thus, the data weighting matrix is inevitably regarded as the identity matrix in the inversion. We present a new inversion scheme, in which the data weighting matrix is automatically obtained from the analysis of the data resolution matrix and its spread function. This approach, named 'extended active constraint balancing (EACB)', assigns a great weighting on the datum having a high resolution and vice versa. We demonstrate that by applying EACB to a two-dimensional resistivity tomography problem, the EACB approach helps to enhance both the resolution and the stability of the inversion process.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.479-489
• /
• 1994

In this paper, the iterative least-squares inversion method is used to determine shapes and density contrasts of 2-D structures from the gravity data. The 2-D structures are represented by their cross-sections of N-sided polygons with density contrasts which are constant or varying with depth. Gravity data are calculated by theoretical formulas for the above structure models. The data are considered as observed ones and used for inversions. The inversions are performed by the following processes: I) polygon's vertices and density contrast are initially assumed, 2) gravity are calculated for the assumed model and error between the true (observed) and calculated gravity are determined, 3) new vertices and density contrast are determined from the error by using the damped least-squares inversion method, and 4) final model is determined when the error is very small. Results of this study show that the shape and density contrast of each model are accurately determined when the density contrast is constant or vertical density gradient is known. In case where the density gradient is unknown, the inversion gives incorrect results. But the shape and density gradient of the model are determined when the surface density contrast is known.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.459-463
• /
• 2003

Waveform inversion requires extracting a reliable low frequency content of seismic data for estimating of the low wave number velocity model. The low frequency content of the seismic data is usually discarded or neglected because of the band-limited response of the source and the receivers. In this study, however small the spectral of the low frequency seismic data is, we assume that it is possible to extract a reliable phase information of the low frequency from the seismic data and use it in waveform inversion. To this end, we exploit the frequency domain finite element modeling and source-receiver reciprocity to calculate the $Frech\`{e}t$ derivative of the phase of the seismic data with respect to the earth model parameter such as velocity, and then apply a damped least squares method to invert the phase of the seismic data. Through numerical example, we will attempt to demonstrate the feasibility of our method in estimating the correct velocity model for prestack depth migration.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.481-485
• /
• 1995

This paper presents a new approaching method to determine the velocity and geometry of shallow subsurface from seismic refraction events. After picking the first breaks from seismic refraction data, we assume that field refraction seismogram can be replaced by the unit delta function having time shift of first break. Time curves are generated by shooting ray tracing. The partial derivatives seismogram for a damped least squares method is computed analytically at each step of the forward ray tracing. The technique is successfully tested on synthetic and real data. It has the advantage of real full waveform inversion, which is robust at low frequency band even if the initial guess is far from the true model.

• Geophysics and Geophysical Exploration
• /
• v.17 no.1
• /
• pp.11-20
• /
• 2014

To reveal shear-wave velocities ($v_s$) and site characterization of Wonju, Korea, Rayleigh waves were recorded at 78 sites of lower altitude using 12 to 24 4.5-Hz vertical geophones for 20 days during the period of February to September 2013. Dispersion curves of the Rayleigh waves obtained by the extended spatial autocorrelation method were inverted using the damped least-squares method to derive $v_s$ models. From these 1-D models, the average $v_s$ to a depth of 30 m ($v_s30$), $v_s$ of weathered rocks, depths to these basement rocks, and average $v_s$ of the overburden layer were derived to be $16.3{\pm}0.7m$, $576{\pm}8m/s$, $290{\pm}7m/s$, and $418{\pm}13m/s$, respectively, in the 95% confidence range. To determine adequate proxies for $v_s30$, we computed correlation coefficients of $v_s30$ with topographic slope (r = 0.46) and elevation (r = 0.43). An empirical linear relationship is presented as a combination of individually estimated $v_s30$ with weighting factors of 0.45, 0.45, and 0.1 for topographic slope, elevation, and mapped lithology, respectively. Due to a weak correlation between $v_s30$ obtained from inversion of dispersion curves and the proxy-based estimation (r = 0.50), however, the relatively large error range should be considered for applications of this relationship.