• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.163-169
• /
• 2005

Recently, the velocity stack domain is having an attention as a very useful domain for various processing in seismic data processing. In order to be used in many applications, the velocity stack should be obtained through an inversion method and the used inversion should have properties like the robustness to noise and the parsimony of velocity stack result. Iteratively Reweighted Least-Squares (IRLS) method is the one of the inversion methods that have such properties. This paper describes the theoretical background, implementation of the method, and examines the characteristics and limits of the IRLS method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.69-72
• /
• 2008

EZTOMO is a crosswell seismic tomography software system. The system has capability of event picking, raytracing, inversion, error analysis, and visualization of the processing results. Waveform of the first arrival signal has been utilized to select the event of the first motion, and uncertainty measured in estimation of the first breaks has been utilized to improve the inversion process.

• The Journal of Engineering Geology
• /
• v.3 no.3
• /
• pp.231-239
• /
• 1993

This paper presents a result on the determination of velocity distribution by a tomographic inversion of crosshole seismic traveltimes in transversely isotropic(aniso tropic) media. The crosshole traveltimes used in this study are synthetic ones computed by ray tracing for some models having isotropic and transversely isotropic velocity distributions. The traveltimes are inverted by a general ART and ansotropic ART which considers the transversely isotropic effect during inver sion. The aniotropic ART gives accurate velodty distributions of transversely isotropic and isotropic models, while the isotropic ART determines accurate velocities only for the isotropic model but inaccurate for the transversely isotropic one. Therefore, the anisotropic ART may be used in case where no information is known on the isotropy or transverse isotropy of a survey area.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.33-42
• /
• 2017

Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

• Journal of the Korean earth science society
• /
• v.30 no.7
• /
• pp.834-844
• /
• 2009

Two-dimensional magnetotelluric (MT) inversion algorithm using minimum support (MS) stabilizer functional was implemented in this study to enhance the contrast of inverted images. For this implementation, this study derived a formula in discrete form for creeping model updates in the least-squares linearized inversion. A spatially varying regularization parameter determination algorithm, which is known as ACB (Active Constraint Balancing), was also adopted to stabilize the inversion process when using MS stabilizer as a model constraint. Inversion experiments for a simple isolated body model show well the feature of MS stabilizer in concentrating the anomalous body compared with the second-order derivative model constraint. This study also compared MS stabilizer and the second-order derivative model constraints for a model having multiple anomalous bodies to show the applicability of the algorithm into field data.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.140-142
• /
• 2002

• Proceedings of the KSEEG Conference
• /
• 2000.04b
• /
• pp.158-159
• /
• 2000

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.299-303
• /
• 2006

Geologic noise, especially located at shallow depth, can be a great obstacle in the interpretation of geophysical data. Thus, it is important to suppress geologic noise in order to accurately detect major anomalous bodies in the survey area. In the inversion of geophysical data, model parameters at shallow depth, which have small size and low contrast of physical property, can be regarded as one of geologic noise. The least-squares method with smoothness constraint has been widely used in the inversion of geophysical data. The method imposes a big penalty on the large model parameter, while a small penalty on the small model parameter. Therefore, it is not easy to suppress small anomalous boies. In this study, we developed a new inversion scheme which can effectively suppress geologic noise by imposing a big penalty on the slowly varying model parameter and a small penalty on the largely varying model parameter. We call the method MTE (main-target emphasizing) inversion. Applying the method to the inversion of 2.5D small loop EM data, we can ensure that it is effective in suppressing small anomalous boies and emphasizing major anomalous bodies in the survey area.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.147-153
• /
• 2007

The conjugate gradient (CG) method is one of the most efficient algorithms for solving a linear system of equations. In addition to being used as a linear equation solver, it can be applied to a least-squares problem. When the CG method is applied to large-scale three-dimensional inversion of magnetotelluric data, two approaches have been pursued; one is the linear CG inversion in which each step of the Gauss-Newton iteration is incompletely solved using a truncated CG technique, and the other is referred to as the nonlinear CG inversion in which CG is directly applied to the minimization of objective functional for a nonlinear inverse problem. In each procedure we only need to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector, significantly reducing the computational requirements needed to do large-scale inversion.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.