• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.207-215
• /
• 2017

To develope shale play which is one of unconventional energy resources, horizontal drilling and hydraulic fracturing are necessary and those are applied to the place where the differential horizontal stress ratio (DHSR) is low. The differential horizontal stress ratio is generally calculated by the minimum and maximum horizontal stress, but it is also calculated from dynamic elastic constants and anisotropic parameters. In this study we analyzed anisotropic properties through the core samples from Danyang limestone and calculated DHSR. The three types of core samples shaped in three directions (vertical, parallel and 45 degree to bedding) were used for laboratory test. We measured P-, S-wave velocities, and density and then calculated dynamic elastic constants, compliance and DHSR. According to the results of the core sample analysis the calculated DHSR is 0.185. Thomsen parameters of the Danyang limestone used in this study are characterized by the P- and S-wave velocities varying along the bedding symmetry axis. It is observed that the DHSR value is more affected by the change in compliance value than the Poisson's ratio. It is necessary to measure SH-wave velocity for more correct petrophysical properties.

• Geophysics and Geophysical Exploration
• /
• v.20 no.4
• /
• pp.216-225
• /
• 2017

Recently, steel casing became an interesting issue when applying controlled-source electromagnetic (EM) method to various fields because effects of steel casing on EM responses are not negligible. This study employed an approach that approximates the steel casing as a series of electric dipole sources in order to develop the numerical algorithm for the efficient simulation of EM responses in the presence of steel casing. After verifying the validity of the developed algorithm, we analyze effects of steel casing on EM responses with the synthetic model simulating geothermal reservoir environment. The analysis showed that the effects of steel casing on EM responses are localized near the casing and increase as the transmitter becomes close to the casing. In addition, through the analysis on the EM responses by the injection of clean water, we confirm that the effects of casing are negligible when interpreting the after-injection data acquired using the transmitter located far enough from the casing. Considering the difference in EM responses between before and after injection in inversion, the effects of the casing can be neglected although after-injection data shows considerable difference due to the close distance between the transmitter and casing. To investigate this kind steel casing effect, the precise analysis on EM responses should be preceded. The algorithm introduced in this study will contribute to the reliable calculations of EM responses distorted by the conductive steel casing.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.210-224
• /
• 2019

Surface wave (SW) surveys, which have been applied to numerous application fields ranging from micro-scale ultrasonic analysis to geological scale analysis, are widely used to monitor near-surface stability. The survey method is basically made through analysis on dispersion of SW propagating along the earth surface, in order to delineate shear velocity structure of subsurface. SW survey data are inverted with assuming one-dimensional (1D) layered-earth in order to recover shear wave velocities of each layer, after being analyzed to make the dispersion curve that shows phase velocity of SW with respect to frequency. This study reviews surface wave surveys with explaining the basic theory including the characteristics of dispersion and the procedure of general data processing. Even though surface wave surveys can be categorized into active and passive methods, this paper focuses only on active surface wave methods which includes continuous SW (CSW), spectral analysis of SW (SASW) and multichannel analysis of SW (MASW). Passive method will be reviewed in the subsequent paper.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.225-238
• /
• 2019

Migration velocity analysis (MVA) for creating optimum depth-domain velocities in seismic imaging was applied to marine long-offset multi-channel data, and the effectiveness of the MVA approach was demonstrated by the combinations of conventional data processing procedures. The time-domain images generated by conventional time-processing scheme has been considered to be sufficient so far for the seismic stratigraphic interpretation. However, when the purpose of the seismic imaging moves to the hydrocarbon exploration, especially in the geologic modeling of the oil and gas play or lead area, drilling prognosis, in-place hydrocarbon volume estimation, the seismic images should be converted into depth domain or depth processing should be applied in the processing phase. CMP-based velocity analysis, which is mainly based on several approximations in the data domain, inherently contains errors and thus has high uncertainties. On the other hand, the MVA provides efficient and somewhat real-scale (in depth) images even if there are no logging data available. In this study, marine long-offset multi-channel seismic data were optimally processed in time domain to establish the most qualified dataset for the usage of the iterative MVA. Then, the depth-domain velocity profile was updated several times and the final velocity-in-depth was used for generating depth images (CRP gather and stack) and compared with the images obtained from the velocity-in-time. From the results, we were able to confirm the depth-domain results are more reasonable than the time-domain results. The spurious local minima, which can be occurred during the implementation of full waveform inversion, can be reduced when the result of MVA is used as an initial velocity model.

• Journal of the Mineralogical Society of Korea
• /
• v.1 no.2
• /
• pp.94-103
• /
• 1988

Diamond anvil cell (DAC) interfaced with a YAG laser heating system has been used to study the phase transformations on perovskite structured titanates ($BaTiO_3$, and $PbTiO_3$) and to synthesize the silicate perovskite phase from the orthopyroxenes of $MgSiO_3$ and $(Mg_{0.87},\;Fe_{0.13})SiO_3$. $BaTiO_3$ and $PbTiO_3$ transform from tetragonal phase to cubic at the pressures of approximately 2.6 GPa and 4.0 GPa at room temperature, respectively. Cubic phases of the both show wide range of stability in the extended in-situ high pressures and high temperature regions. Starting orthoenstatite of $MgSiO_3$ has yielded the perovskite phase as the major structure with ilmenite, gamma-spinel, betta-spinel and stishovite phases at ~38 GPa and ${\sim}1,000^{\circ}C$. $(Mg_{0.87},\;Fe_{0.13})SiO_3$ has shown the perovskite as the major phase with betta-spinel, stishovite and enstatite phases at ~35 GPa and ${\sim}1,000^{\circ}C$. The ilmenite phase does not occur at this condition.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.56-63
• /
• 2002

Tracing the history of study, problems of seawater intrusion are commonly investigated with electrical techniques because seawater saturated zone is indicative of the low resistivity anomaly. There we, however, silt and mud layers in the western and southern coastal areas of Korea, so we may make a mistake in case we determine seawater intrusion only with resistivity survey. Hence, reference IP survey was carried out in Kimje, Jeollabuk-Do and Youngkwang, Jeollanam-Do in order to decide whether or not the area is under the influence of seawater intrusion. With the use of a electric field cable to minimize EM coupling, we obtained more accurate results by appling reference If technique measuring simultaneously wavelet of current as well as potential. With the aid of reference IP technique, it is possible that we can exactly evaluate seawater intrusion by discriminating seawater saturated area (no IP effect) from very highly conductive layer composed of clay mineral (high IP effect).

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.23-32
• /
• 2002

Capacitively-coupled resistivity (CCR) system is known to be very useful where galvanic contact to earth is impossible, such as the area covered with thick ice, snow, concrete or asphalt. This system injects current non-galvanically, i.e., capacitively to earth through line antenna and measures potential difference in a same manner. We derived geometric factor for two types of antenna configuration and presented the method of processing and converting the data obtained with CCR system suitable to conventional resistivity inversion analysis. The CCR system, however, has limitations on use at conductive area or electrically noisy area since it is very difficult to inject sufficient current to earth with this system as with conventional resistivity system. This causes low SM ratio when acquiring data with CCR system and great care must be taken in acquiring data with this system. Additionally the uniform contact between line antennas and earth is also crucial factor to obtain good S/N ratio data. The CCR method, however, enables one to perform continuous profiling over a survey line by dragging entire system and thus will be useful in rapid investigation of conductivity distribution in shallow subsurface.

• 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.2 no.4
• /
• pp.180-183
• /
• 1999

A proto-type seismic measurement system with a reference was developed to improve S/N (signal-to-noise ratio) of seismic data, especially in noisy urban areas. Two pairs of correlation measurements (the one for microphone and geophone, and another for electromagnetic (EM) loop and geophone) were carried out near Kimpo Airport and at Kimje. The spectrum analyses were also performed to investigate the correlation of two pairs of time series; one for microphone and geophone, and another for EM loop and geophone. The sound waves measured with the microphone and the geophone are highly correlated. However, differences in the reponses are readily identifiable across 200 Hz; in the vicinity of 100 Hz, the spectral energy for geophone is 20 dB higher than that for microphone, and at near 500 Hz, the spectral energy for microphone is 30 dB higher than that for geophone. Overall, the spectral energy appears concentrated on the frequency window below 600 Hz for geophone. It contrasts with the observation of dominant frequency at the range of above 200 Hz for microphone. The wave forms of EM noise (due to an ACDC inverter) measured with EM loop and geophone are consistently and highly correlated each other. The power spectrum of the EM noise for EM loop shows that the spectral energies at odd harmonic frequencies of 60 Hz are higher than those at even harmonic frequencies of 60 Hz. It is compared to the power spectrum for geophone; the spectral energies at odd harmonics are nearly same as those at even harmonic frequencies.