• 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.19 no.1
• /
• pp.1-10
• /
• 2016

The hypocenter distribution of microseismic events generated by hydraulic fracturing for shale gas development provides essential information for understanding characteristics of fracture network. In this study, we evaluate how inaccurate velocity models influence the inversion results of two widely used location programs, hypoellipse and hypoDD, which are developed based on an iterative linear inversion. We assume that 98 stations are densely located inside the circle with a radius of 4 km and 5 artificial hypocenter sets (S0 ~ S4) are located from the center of the network to the south with 1 km interval. Each hypocenter set contains 25 events placed on the plane. To quantify accuracies of the inversion results, we defined 6 parameters: difference between average hypocenters of assumed and inverted locations, $d_1$; ratio of assumed and inverted areas estimated by hypocenters, r; difference between dip of the reference plane and the best fitting plane for determined hypocenters, ${\theta}$; difference between strike of the reference plane and the best fitting plane for determined hypocenters, ${\phi}$; root-mean-square distance between hypocenters and the best fitting plane, $d_2$; root-mean-square error in horizontal direction on the best fitting plane, $d_3$. Synthetic travel times are calculated for the reference model having 1D layered structure and the inaccurate velocity model for the inversion is constructed by using normal distribution with standard deviations of 0.1, 0.2, and 0.3 km/s, respectively, with respect to the reference model. The parameters $d_1$, r, ${\theta}$, and $d_2$ show positive correlation with the level of velocity perturbations, but the others are not sensitive to the perturbations except S4, which is located at the outer boundary of the network. In cases of S0, S1, S2, and S3, hypoellipse and hypoDD provide similar results for $d_1$. However, for other parameters, hypoDD shows much better results and errors of locations can be reduced by about several meters regardless of the level of perturbations. In light of the purpose to understand the characteristics of hydraulic fracturing, $1{\sigma}$ error of velocity structure should be under 0.2 km/s in hypoellipse and 0.3 km/s in hypoDD.

• Economic and Environmental Geology
• /
• v.31 no.2
• /
• pp.127-139
• /
• 1998

A new technique of simultaneous inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the southeast part of the Korean Peninsula including Pohang Basin, Kyongsang Basin and Ryongnam Massif. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 44 events with 554 seismic rays are inverted for locations and crustal structure. $6{\times}6$ with $0.5^{\circ}$ and 8 layers (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from surface to Moho, ten profiles along latitude and longitude are analyzed. The results are as follows: 1) the average velocity and thickness of sediment are 5.04 km/s and 3-4 km, and the velocity of basement is 6.11 km/s. The shape of velocity in shallower layer is agreement with Bouguer gravity anomaly (Cho et al., 1997). 2) the velocities fluctuate strongly in the upper crust. The velocity distribution of the lower crust under Conrad appears basically horizontal. 3) the average depth of Moho is 30.4 km, and velocity is 8.01 km/s. 4) from the velocity and depth of the sediment, the thickness, velocity and form of the upper crust, and the depth and form of Moho, we can find the obvious differences among Ryongnam Massif, Kyongsang Basin and Pohang Basin. 5) the deep faults (a Ulsan series faults) near Kyongju and Pohang areas can be found to be normal and/or thrust faults with detachment extended to the bottom of the upper crust.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.117-122
• /
• 2006

In the geophysical monitoring to understand the change of subsurface material properties with time, the time-invariant static subsurface model is commonly adopted to reconstruct a time-lapse image. This assumption of static model, however, can be invalid particularly when fluid migrates very quickly in highly permeable medium in the brine injection experiment. In such case, the resultant subsurface images may be severely distorted. In order to alleviate this problem, we develop a new least-squares inversion algorithm under the assumption that the subsurface model will change continuously in time. Instead of sampling a time-space model into numerous space models with a regular time interval, a few reference models in space domain at different times pre-selected are used to describe the subsurface structure continuously changing in time; the material property at a certain space coordinate are assumed to change linearly in time. Consequently, finding a space-time model can be simplified into obtaining several reference space models. In order to stabilize iterative inversion and to calculate meaningful subsurface images varying with time, the regularization along time axis is introduced assuming that the subsurface model will not change significantly during the data acquisition. The performance of the proposed algorithm is demonstrated by the numerical experiments using the synthetic data of crosshole dc resistivity tomography.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.167-173
• /
• 1999

Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

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

• Economic and Environmental Geology
• /
• v.35 no.5
• /
• pp.445-454
• /
• 2002

3-D P-wave velocity model in the southern Korean Peninsula is investigated by using the earthquake tomography method. This velocity model would be used to locate the exact hypocenter position, and also useful for our understanding of the crustal structure. The simultaneous inversion is used to get the minimum 1-D model and hypo-center relocation, which are used as an initial 3-D velocity model. The velocities in the minimum 1-D model are 6.04 km/s, 6.45 km/s, and 7.78 km/s between the depth of 0-19 km, 19-32 km, and 32-55 km respectively. In the 3-D P-wave velocity model, Layer 1 (0~3 km) has high velocities in Kyongsang basin, Yonglam massif, and Okchon folded belt, and low velocities in Kyonggi massif. In layer 2 (3~19 km) high velocities are predominent around Kyonsang basin and Yongnam massif except Yonil basin, but low velocities exist around Kyonggi massif and Okchon folded belt. In Laye. 3 (19~32 km) high velocities prevail throughout the southern part of Korean Peninsula, but low velocity does throughout the middle except SNU, YIN station in Konggi massif. In Layer 4 (32 km), the maximum velocity is showed in the middle and southwestern part, while the minimum velocity in the southeastern and coastal area. The depth of the velocity boundary corresponds to the crustal structure of the southern Korean Peninsula which is calculated by gravity data.

• Journal of the Korean Geophysical Society
• /
• v.8 no.4
• /
• pp.211-214
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10 m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

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

• Economic and Environmental Geology
• /
• v.39 no.4 s.179
• /
• pp.473-483
• /
• 2006

Although application of electrical methods in Korea began with observation of self potentials before World War II, the methods were developed slowly by the beginning of 1980's when a major burst of development activity took place. DC resistivity methods are applied in Korea more to geotechnical problems rather than to environmental ones unlike other developed countries. As with every other branch of technology, the evolving speed of the silicon chip and of streaming data to hard disk has revolutionized data collection and noise reduction processing. The last two decades saw major advances in data collection, processing, and interpretation of electrical data. Development of smooth-model two-dimensional (2D) resistivity inversion is one of the most visible changes to geophysical interpretation of the last 40 years and is now routinely applied to apparent resistivity data. The ability to represent resistivities in section rather than pseudosection view has revolutionized interpretation. Although calculation of sensitivities for general electromagnetic problems require numerous forward modelings, DC resistivity methods can enjoy computational efficiencies if sources and receivers occupy the same position, and previously intractable 3D inversion is now becoming available.