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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.05a
• /
• pp.47-56
• /
• 2009

Proper reservoir characterization is an integral part of formation evaluation, reserve estimation and planning of field development. Seismic inversion is a widely employed reservoir characterization tool that provides various rock properties of reservoir intervals. This study presents results of the inversion studies including Geostatistical Inversion carried out on the gas fields, offshore Myanmar. Higher resolution and multiple models can be produced by Geostatistical Inversion using input data such as pre-stack seismic data, well logs, petrophysical relationships and geological inferences for example reservoir shape and lateral extent. Detailed reservoir characterization was required for the development plan of gas fields, and the Geostatistical Inversion studies served as a basis for integrated geological modeling and development well planning.

• Geotechnical Engineering
• /
• v.14 no.3
• /
• pp.95-108
• /
• 1998

Downhole testing and seismic CPT (SCPT) have been widely used to evaluate stiffness profiles of the subgrade. Advantages of downhole testing and SCPT such as low cost, easy operation and a simple seismic source have got these testings more frequently adopted in site investigation. For the automated analysis of downhole testing and SCPT, the concept of interval measurements has been practiced. In this paper. a new inversion procedure to deal tilth the interval measurements for the automated downhole testing and SCPT (including a newlydeveloped continuous SCPT) is proposed. The forward modeling in the new inversion procedure incorporates ray path theory based on Snell's law. The formulation for the inversion analysis is derived from the maximum likelihood approach, which estimates the maximum likelihood of obtaining a particular travel time from a source to a receiver. Verification of the new inversion procedure was performed with numerical simulations of SCPT using synthesized profiles. The results of the inversion analyses performed for the synthetic data show that the new inversion analysis is a valid procedure which enhances Va profiles determined by downhole testing and SCPT.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.33-45
• /
• 2002

It was necessary to devise new techniques to overcome and enhance the resolution limits of traveltime tomography. Waveform inversion has been one of the methods for giving very high resolution result. High resolution image could be acquired because waveform inversion used not only phase but amplitude. But waveform inversion was much time consuming Job because forward and backward modeling was needed at each iteration step. Velocity-stress method was used for effective modeling. Resolution limits of imaging methods such as travel time inversion, acoustic and elastic waveform inversion were investigated with numerical models. it was investigated that Resolution limit of waveform inversion was similar tn resolution limit of migration derived by Schuster. Horizontal resolution limit could be improved with increased coverage by adding VSP data in cross hole that had insufficient coverage. Also, waveform inversion was applied to realistic models to evaluate applicability and using initial guess of travel time tomograms to reduce non-linearity of waveform inversion showed that the better reconstructed image could be acquired.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.302-307
• /
• 2005

Love wave and Rayleigh wave are the major elastic waves belonging to the category of the surface wave. The fact that Love wave is not contaminated by P-wave which makes Love wave superior to Rayleigh wave and other body waves. Therefore, the information that Love wave carries is more distinct and clearer than the information of Rayleigh wave. Based on theoretical research, the joint inversion analysis which is used both Love wave dispersion information and Rayleigh wave dispersion information was proposed. Purpose of the joint inversion analysis is to improve accuracy and convergency of inversion results utilizing that frequency contribution of each wave is different. This analysis technique is consisted of the forward modeling using transfer matrix, the sensitivity matrix determined to the ground system and DLSS(Damped Least Square Solution) as a inversion technique. The application of this analysis was examined through the field test.

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.643-646
• /
• 2005

The gas hydrate exploration using seismic reflection data, the detection of BSR(Bottom Simulating Reflector) on the seismic section is the most important work flow because the BSR have been interpreted as being formed at the base of a gas hydrate zone. Usually, BSR has some dominant qualitative characteristics on seismic section i.e. Wavelet phase reversal compare to sea bottom signal, Parallel layer with sea bottom, Strong amplitude, Masking phenomenon above the BSR, Cross bedding with other geological layer. Even though a BSR can be selected on seismic section with these guidance, it is not enough to conform as being true BSR. Some other available methods for verifying the BSR with reliable analysis quantitatively i.e. Interval velocity analysis, AVO(Amplitude Variation with Offset)analysis etc. Usually, AVO analysis can be divided by three main parts. The first part is AVO analysis, the second is AVO modeling and the last is AVO inversion. AVO analysis is unique method for detecting the free gas zone on seismic section directly. Therefore it can be a kind of useful analysis method for discriminating true BSR, which might arise from an Possion ratio contrast between high velocity layer, partially hydrated sediment and low velocity layer, water saturated gas sediment. During the AVO interpretation, as the AVO response can be changed depend upon the water saturation ratio, it is confused to discriminate the AVO response of gas layer from dry layer. In that case, the AVO modeling is necessary to generate synthetic seismogram comparing with real data. It can be available to make conclusions from correspondence or lack of correspondence between the two seismograms. AVO inversion process is the method for driving a geological model by iterative operation that the result ing synthetic seismogram matches to real data seismogram wi thin some tolerance level. AVO inversion is a topic of current research and for now there is no general consensus on how the process should be done or even whether is valid for standard seismic data. Unfortunately, there are no well log data acquired from gas hydrate exploration area in Korea. Instead of that data, well log data and seismic data acquired from gas sand area located nearby the gas hydrate exploration area is used to AVO analysis, As the results of AVO modeling, type III AVO anomaly confirmed on the gas sand layer. The Castagna's equation constant value for estimating the S-wave velocity are evaluated as A=0.86190, B=-3845.14431 respectively and water saturation ratio is $50\%$. To calculate the reflection coefficient of synthetic seismogram, the Zoeppritz equation is used. For AVO inversion process, the dataset provided by Hampson-Rushell CO. is used.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.291-298
• /
• 2002

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.231-236
• /
• 2000

In order to predict depth of the pile forward modeling and inversion of magnetic logging data was conducted by using a finite line of dipoles model. The horizontal component as well as the vertical component of magnetic fields can be measured in the borehole, and the magnetic anomalies can be obtained by subtracting the Earth's magnetic field from the measurement. The magnetic anomalies of the pile are considered as vector sum of induced magnetization due to the Earth's magnetic field and remnant magnetization possessed by steel strings in the pile. The magnetic anomalies are used as input data for inversion from which the length, the magnetic moment per unit length, and the dip angle of the pile can be obtained. From the inversion of synthetic noisy data, and the data obtained from the field model test it is found that the driving depth of the pile can be determined as close to the order of measuring interval (5∼10㎝). It is also found that the resultant magnetic anomalies due to an individual steel string in the pile are almost same as those due to a group of steel strings located at the center of the pile. The magnetic logging method also can be used for locating reinforced bars, pipes, and steel casings.

• 제어로봇시스템학회:학술대회논문집
• /
• 1996.10b
• /
• pp.1500-1503
• /
• 1996

This paper combines the disturbance accommodating control(DAC) and nonlinear model inversion control for a skid-to-turn(STT) missile. The missile autopilot may be designed to be robust with respect to a variety of uncertainties. We proposes the two step control design method. Nonlinear model inversion control is used as the main design method. Due to the model uncertainties and external disturbances, the exact nonlinear model inversion can not be achieved. DAC is designed to detect, to identify, and to compensate these uncertainties. DAC's disturbance observer is linear. Thus it is easy to implement. It does not cause the convergence problem due to coexistence between the modeling uncertainties and external disturbances. 6 DOF simulation results show that the proposed method may improve the missile tracking performance.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.129-135
• /
• 2012

The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.