• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.50-59
• /
• 2007

The advent of towed geo-electrical array surveying on water and land has resulted in datasets of magnitude approaching that of airborne electromagnetic surveying and most suited to 1D inversion. Robustness and complete automation is essential if processing and reliable interpretation of such data is to be viable. Sharp boundaries such as river beds and the top of saline aquifers must be resolved so use of smoothness constraints must be minimised. Suitable inversion algorithms must intelligently handle low signal-to-noise ratio data if conductive basement, that attenuates signal, is not to be misrepresented. A noise-level aware inversion algorithm that operates with one elastic thickness layer per electrode configuration has been coded. The noise-level aware inversion identifies if conductive basement has attenuated signal levels so that they are below noise level, and models conductive basement where appropriate. Layers in the initial models are distributed to span the effective depths of each of the geo-electric array quadrupoles. The algorithm works optimally on data collected using geo-electric arrays with an approximately exponential distribution of quadrupole effective depths. Inversion of data from arrays with linear electrodes, used to reduce contact resistance, and capacitive-line antennae is plausible. This paper demonstrates the effectiveness of the algorithm using theoretical examples and an example from a salt interception scheme on the Murray River, Australia.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.99-104
• /
• 2009

We propose a waveform inversion method for SH-wave data obtained in a shallow seismic refraction survey, to determine a 2D inhomogeneous S-wave profile of shallow soils. In this method, a 2.5D equation is used to simulate SH-wave propagation in 2D media. The equation is solved with the staggered grid finite-difference approximation to the 4th-order in space and 2nd-order in time, to compute a synthetic wave. The misfit, defined using differences between calculated and observed waveforms, is minimised with a hybrid heuristic search method. We parameterise a 2D subsurface structural model with blocks with different depth boundaries, and S-wave velocities in each block. Numerical experiments were conducted using synthetic SH-wave data with white noise for a model having a blind layer and irregular interfaces. We could reconstruct a structure including a blind layer with reasonable computation time from surface seismic refraction data.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.09a
• /
• pp.3-7
• /
• 2003

The velocity structure beneath the CHNB broadband station is determined by receiver function analysis using by from teleseismic P waveforms. The detailed broadband receiver functions are obtained by stacking method for source-equalized vertical, radial and tangential components of teleseismic P waveforms. A time domain inversion uses the stacked radial receiver function to determine vertical P wave velocity structure beneath the station. The crustal velocity structures beneath the stations are estimated using the receiver function inversion method in the case at the crustal model parameterized by many thin, flat-lying, homogeneous layers. Events divide into 4 groups. four azimuths corresponding to events in group a(southwest), b(south), c(southeast), d(northeast). The result of crust at model inversion shows the crustal velocity structure beneath the CHNB station varies smoothly with increasing depth. The conard discontinuity lies around 18 km and moho discontinuity lies range from 30 to 34 km.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.261-270
• /
• 2006

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeju Island has been carried out. Broad band MT sounding curves with good quality could be gathered by performing audio-frequency magnetotelluric (AMT) survey during the MT survey and by operating the remote reference in Kyushu Island, Japan. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands for the field data as well as for the data from numerical 2-D modeling said that high frequency information from AMT survey can be useful for interpreting not only the shallow part but also the deep structures, especially when the formation is resistive. The 2-D inversion models of field data show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. It can be either the effect of the surrounding sea water, or the structures due to ancient volcanic events. But unfortunately by now, we do not have any further information about the anomaly.

• 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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.303-306
• /
• 2007

After Kyoto protocol took effect, many countries are making efforts to reduce $CO_2$ one of effective which is geosequestration. But a percentage of geosequestration in total research budget is very small and the priority order of research also is receded in Korea. As one of efforts to activate the research on geosequestration in field of geophysics, we proposed the plan to build up a measurement system for rock physical properties monitoring during $CO_2$ injection which will function as original technology. The system consists of two part, one of which is a data acquisition system based on seismic and complex resistivity tomographic measurement and the other of which is a tri-axial compressive system to realize the in-situ condition. And also developments of various inversion algorithms are proposed to interpret data qualitatively such as a inversion algorithm for confined cylindrical boundary, a joint inversion algorithm and a 4-D inversion algorithm.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.199-206
• /
• 2017

Field data obtained from marine exploration are influenced by various environmental factors such as wind, waves, tidal current and flow velocity of a background medium. Most environmental factors except for the flow velocity are properly corrected in the data processing stage. In this study, the wave equation modeling considering flow velocity is used to generate observation data, and numerical experiments using the observation data were conducted to analyze the effect of flow velocity on waveform inversion. The numerical examples include the results with unrealistic flow velocities. In addition, an algorithm is suggested to numerically extract flow velocity for waveform inversion. The proposed algorithm was applied to the modified Marmousi2 model to obtain the results depending on the flow velocity. The effect of flow velocity on updated physical properties was verified by comparing the inversion results without considering flow velocity and those obtained from the proposed algorithm.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.326-334
• /
• 2008

The resistivity method has been used to image the electrical properties of the subsurface. Especially, this method has become suitable for monitoring since data could be rapidly and automatically acquired. In this study, we developed a time-lapse inversion algorithm for the interpretation of resistivity monitoring data. The developed inversion algorithm imposes a big penalty on the model parameter with small change, while a minimal penalty on the model parameter with large change compared to the reference model. Through the numerical experiments, we can ensure that the time-lapse inversion result shows more accurate and focused image where model parameters have changed. Also, applying the timelapse inversion method to the leakage detection of an embankment dam, we can confirm that there are three major leakage zones, but they have not changed over time.

• Journal of the Korean Geophysical Society
• /
• v.4 no.1
• /
• pp.21-33
• /
• 2001

Simultaneous inversion of first-arrivals of local earthquakes recorded by the Korea Meteorological Administration(KMA) seismograph network from 1991 to 1998 is made to derive 1D crustal velocity structure of the Korean peninsula. Twenty-nine events with 178 observations are used in the inversion. Average crustal P-wave velocity turns out to be about 6.3 km/sec, and crustal thickness and upper mantle P-wave velocity are estimated as 33 km and 7.9 km/sec, respectively. Results of inversion indicate the possibility of the low velocity layer in the lower crust. Joint inversion is applied to estimate hypocenters, station delays, and velocities simultaneously. Relative station corrections for 11 stations range from zero to about 1.2 sec. Analysis of the synthetic data shows that estimates of hypocenter locations and station corrections as well as averaged crustal structure are reliable for the given data set..

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.13 no.7
• /
• pp.716-725
• /
• 2002

In this paper, we showed the applicability of power symmetric MSSI(Mid-Span Spectral Inversion) to the long-haul WDM system with the channel interference due to FWM(Four Wave Mixing). And we showed the degree of performance improvement. We used 1 dB EOP(Eye-Opening Penalty) criterion so as to evaluate the degree of compensation dependent on the variation of chirp parameter of optical pulse for the various input power in high speed tansmission system. And we evaluated the maximum input power of channel be able to be the signal to crosstalk noise (SNR) above 20 dB in the transmission link with the channel interference due to FWM. Consequently the proposed MSSI compensation method is capable to transmitting the total 68 WDM channels simultaneously with a 0.4 nm channel spacing and 5.3 dBm maximum input power in a 10 Gbps transmission link. Therefore the proposed power symmetric MSSI compensation method may be very useful for the implementation of long-haul wideband WDM transmission systems with relatively high power and improved performance.