• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.66-73
• /
• 2014

To estimate the shear-velocity ($v_s$) structure beneath the WIZ station on White Island in New Zealand, we applied receiver function (RF) inversion and H-${\kappa}$ stacking methods to 362 teleseismic events (Mw > 5.5) recorded during April 20, 2007 to September 6, 2013. Using 71 RFs with errors less than 20% after 200 iterative computations, we determined that the depth to Moho of $v_s$ = 4.35 km/s is $24{\pm}1km$ within a 15 km radius of the station. In an 1-d $v_s$ model derived by RF inversions, a 4-km thick low-velocity layer (LVL) at depths of 18 ~ 22 km was identified in the lower crust. This LVL, which is 0.15 km/s slower than the rocks above and below it, may indicate the presence of a deep magma reservoir. The H-${\kappa}$ stacking method yielded an estimate of the depth to the Moho of 24.5 km, which agrees well with the depth determined by RF inversions. The low $v_p/v_s$ ratio of 1.64 may be due to the presence of gas-filled rock or hot crystallizing magma.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed, when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.152-157
• /
• 2011

The small-loop electromagnetic (EM) method is one of the rapid and non-destructive geophysical methods and has been used widely for many geophysical investigations, particularly for shallow engineering and environmental surveys. Especially in the shallow marine environment, the small-loop EM technique is very effective because of rapid and convenient data acquisition, large signal and low noise level. However, the method has been rarely applied in the very conductive marine environment since it's penetration or investigation depth might be considered too low. In this study, we demonstrated that the small-loop EM method can be effectively applied in the extremely conductive marine environment through the analysis of 1D small-loop EM data. Furthermore, we confirmed that the resistivity distribution under the sea bottom can be quantitatively predicted from the 1D inversion results of synthetic and field data.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.21-36
• /
• 1999

Recently resistivity survey is widely used for site investigations in the field of civil engineering. Since such application area requires accurate interpretation tools especially in the area of complicated geology and rough terrain topography, we developed a three-dimensional (3-D) resistivity inversion code, which can reconstruct real earth structures. Furthermore, the inversion code gives resolution-enhanced images by applying the ACB(Active Constraint Balancing) method. With the help of this inversion code, 3-D resistivity survey is now used as new techniques for site investigations in civil engineering problem. By imaging the 3-D resistivity distribution, we could get useful informations such as depth distribution of basement rock, distribution of weak zone, fractures and cavities which is crucial to civil engineers.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.124-130
• /
• 2007

Least squares ($l^2$ norm) solutions of seismic inversion tend to be very sensitive to data points with large errors. The $l^1$ norm minimization gives more robust solutions, but usually with higher computational cost. Iteratively reweighted least squares (IRLS) method gives efficient approximate solutions of these $l^1$ norm problems. I propose an efficient implementation of the IRLS method for a hybrid $l^1/l^2$ minimization problem that behaves as $l^2$ norm fit for small residual and $l^1$ norm fit for large residuals. The proposed algorithm shows more robust characteristics to the decision of the threshold value than the l1 norm IRLS inversion does with respect to the threshold value to avoid singularity.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.217-224
• /
• 2013

We developed a time-lapse inversion using new cross-model constraints based on change ratio and resolution of model parameters. The cross-model constraint based on change ratio imposes the same penalty on the model parameters with equal change ratio. This constraint can emphasize the model parameters with significant change regardless of their increase or decrease. The resolution cross-model constraint imposes a small penalty on the model parameters with poor resolution, but a large penalty on the model parameters with good resolution. Thus, the model parameter with poor resolution can be effectively identified in the inversion result if they are significantly changed with time. Through the numerical tests for 3D resistivity monitoring data sets, the performance of these two cross-model constraints was confirmed. Finally, for the safety estimation of a sea dyke, we applied the developed time-lapse inversion to the 3D resistivity monitoring data that were acquired at a sea dike located in western coastal area of Korea. The result of time-lapse inversion suggested that there were no significant changes at the sea dike during the monitoring period.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.1-11
• /
• 2019

Well logging technologies are used to measure the physical properties of reservoirs through boreholes. These technologies have been utilized to understand reservoir characteristics, such as porosity, fluid saturation, etc., using equations based on rock physics models. The analysis of well logs is performed by selecting a reliable rock physics model adequate for reservoir conditions or characteristics, comparing the results using the Archie's equation or simandoux method, and determining the most feasible reservoir properties. In this study, we developed a joint inversion algorithm to estimate physical properties in shaly sandstone reservoirs based on the pre-existing algorithm for sandstone reservoirs. For this purpose, we proposed a rock physics model with respect to shale volume, constructed the Jacobian matrix, and performed the sensitivity analysis for understanding the relationship between well-logging data and rock properties. The joint inversion algorithm was implemented by adopting the least-squares method using probabilistic approach. The developed algorithm was applied to the well-logging data obtained from the Colony gas sandstone reservoir. The results were compared with the simandox method and the joint inversion algorithms of sand stone reservoirs.

• The Journal of the Acoustical Society of Korea
• /
• v.23 no.4
• /
• pp.268-273
• /
• 2004

A method of reflection-transmission type ultrasonic inverse scattering image was presented using linearly arrayed transducers in inner surface of half-cylinder. In this method, to reduce the number of data, the mirror effect using a reflector behind object and pulse wave with finite frequency band, To verify the proposed method, a computer simulation was performed for organic phantom specimen, As the results. it was verified that the reconstructed image was satisfactory even when the limitation view angle was limited to around 30 deg.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.23 no.4
• /
• pp.524-532
• /
• 2012

In this paper, we used MOM-based BIM(Born Iterative Method) algorithm to implement the inverse scattering for the detection of cancer. We adopted two-dimensional breast structure, integral equations and two-dimensional Green's function is solved with MoM(Method of Moment) to analyzing electromagnetic scattering phenomena. In addition, verifying the calculation of developed inverse scattering algorithm and analyzing medical applicability and limitations of the algorithm.