• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.252-265
• /
• 2022

Subsurface velocity is critical for the accurate resolution geological structures. The estimation of acoustic impedance is also critical, as it provides key information regarding the reservoir properties. Therefore, researchers have developed various inversion approaches for the estimation of reservoir properties. The Markov chain Monte Carlo method, which is a stochastic method, has advantages over the deterministic method, as the stochastic method enables us to attenuate the local minima problem and quantify the uncertainty of inversion results. Therefore, the Markov chain Monte Carlo inversion method has been applied to various kinds of geophysical inversion problems. However, studies on the Markov chain Monte Carlo inversion are still very few compared with deterministic approaches. In this study, we reviewed various types of reversible jump Markov chain Monte Carlo applications and explained the key concept of each application. Furthermore, we discussed future applications of the stochastic method.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.131-137
• /
• 2007

With consideration of the first arrival picking methodology and inherent errors in picking process, I propose, from the computerization point of view, a practical algorithm for picking and error computation. The proposed picking procedure consists of 2-step; 1) picking the first coherent peak or trough events, 2) derive a line which approximates the record in the interval prior to the pick, and set the intercept time of the line as the first break. The length of fitting interval used in experiment, is few samples less than 1/4 width of the arriving wavelet. A quantitative measure of the error involved in first arrival picking is defined as the time length that needed to determine if an event is the first arrival or not. The time length is expressed as a function of frequency bandwidth of the signal and the S/N ratio. For 3 sets of cross-well seismic data, first breaks are picked twice, by manually, and by the proposed method. And at the same time, the error bound for each trace is computed. Experiment results show that good performance of the proposed picking method, and the usefulness of the quantitative error measure in pick-quality evaluation.

• Journal of the Korean Geophysical Society
• /
• v.5 no.3
• /
• pp.175-198
• /
• 2002

The techniques using geophysical methods were adopted to obtain quantitative criteria for assessment of grouting effectiveness. Various surface and borehole geophysical surveys including seismic, GPR(ground penetrating radar), resistivity and electromagnetic methods were conducted in fractured rock pilot site before and after grouting execution. However, it is not enough that geophysical data provide criteria for field engineers to confirm the grouting effectiveness in that site even though there is somewhat difference before and after grouting. This study will be continued for the detailed criteria and assessment of grouting effectiveness in other sites.

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

The first-arrival traveltime calculation method based on the damped wave equation overcomes the shortcomings of ray-tracing methods. Since this algorithm needs to solve the damped wave equation, numerical modeling is essential. However, it is not desirable to use the finite-difference method (FDM), which has good computational efficiency, for simulating the land seismic data because of irregular topography. Thus, the finite-element method (FEM) which requires higher computational cost than FDM has been used to correctly describe the irregular topography. In this study, we computed first-arrival traveltimes in an irregular topographic model using FDM incorporating embedded boundary method (EBM) to overcome this problem. To verify the accuracy and efficiency of the proposed algorithm, we compared our results with those of FEM. As a result, the proposed method using EBM not only provided the same accuracy as the FEM but also showed the improved computational efficiency.

• The Journal of Engineering Geology
• /
• v.30 no.3
• /
• pp.289-302
• /
• 2020

To determine the shallow subsurface structure and sliding surface of land creeping in 2016 at Hadong-gun, Gyeongsangnam-do, geophysical surveys (electric resistivity, and refraction seismic methods, borehole televiewer) and slope stability analysis were conducted. The subsurface structure delineated with borehole lithologies and seismic velocity structures provided the information that the sediment layer on the top of the slope was rather as thick as 20 m and the underlying weathered rock (anorthosite) was thinner than 1 m. Based on the tension cracks observed during the geological mapping, televiewer scanning was performed at the borehole BH-2 and detected the intensive fracture zones at the ground-water level, associated with the slip weak zones mapped in dipole-dipole electrical resistivity section. Downslope sliding and slightly upward pushing at the apex of high resistive bedrock explains the curved slip plane of the land creeping. Such a convex structure might play a role of natural toe abutment for preventing the downward development of slip weak zones. In slope stability analysis, the safety factors of the slip weak zone are calculated with varying the groundwater levels for dry and rainy seasons and the downslope is founded to be unstable with safety factor of 0.89 due to fully saturated material in rainy season.

• Journal of the Korean earth science society
• /
• v.24 no.8
• /
• pp.657-667
• /
• 2003

Deep-drilled core sampling and high-resolution seismic survey were carried out to identify a Holocene-late Pleistocene boundary in Gyeonggi Bay, western coast of Korea. Analysis of core sections revealed the existence of an oxidized and semi-consolidated sediment layer, Iying immediately below a Holocene horizon (Unit I) and being developed at the top of a late Pleistocene deposit (Unit II). The oxidized sedimentary layer (uppermost part of Unit II) is characterized by semi-consolidated, yellowish sediments showing signs of desiccation and alteration such as high N value, low water content, periglacial cryogenic structure, depletion of smectite, and high geochemical weathering index (Ba/Sr ratio). This feature, together with radiocarbon ages, suggests that the layer has formed as a result of prolonged subaerial exposure of Unit II sediments during the late Wisconsin sea-level lowstand, producing a regional unconformity. Such unconformitic-bounding surface corresponds to a prominent near-surface reflector (R), which is observed in seismic profiles obtained across the drilled-core sections in the study area. Consequently, the buried oxidized-sedimentary layer associated with the seismic reflector possibly plays a key horizon for the understanding of late Quaternary environmental changes as well as evidence of the emergence of the Yellow Sea shelf during the late Wisconsin sea-level lowstand.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.211-216
• /
• 2013

In general, smoothing filters regularize functions by reducing differences between adjacent values. The smoothing filters, therefore, can regularize inverse solutions and produce more accurate subsurface structure when we apply it to full waveform inversion. If we apply a smoothing filter with a constant coefficient to subsurface image or velocity model, it will make layer interfaces and fault structures vague because it does not consider any information of geologic structures and variations of velocity. In this study, we develop a selective smoothing regularization technique, which adapts smoothing coefficients according to inversion iteration, to solve the weakness of smoothing regularization with a constant coefficient. First, we determine appropriate frequencies and analyze the corresponding wavenumber coverage. Then, we define effective maximum wavenumber as 99 percentile of wavenumber spectrum in order to choose smoothing coefficients which can effectively limit the wavenumber coverage. By adapting the chosen smoothing coefficients according to the iteration, we can implement multi-scale full waveform inversion while inverting multi-frequency components simultaneously. Through the successful inversion example on a salt model with high-contrast velocity structures, we can note that our method effectively regularizes the inverse solution. We also verify that our scheme is applicable to field data through the numerical example to the synthetic data containing random noise.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.323-329
• /
• 2010

Three deconvolution methods were applied to stacked seismic data obtained to investigate gas-hydrates in the Ulleung Basin, East Sea: (1) minimum-phase spiking deconvolution, (2) minimum-phase spiking deconvolution using an averaged wavelet from all traces, and (3) deterministic deconvolution using a wavelet with phases computed from well-logs. We analyzed the resolving property of these methods for lithological boundaries. The first deconvolution method increases temporal resolution but decreases lateral continuity. The second method shows, in an overall sense, similar results to the spiking deconvolution using a minimum phase wavelet for each trace; however, it results in a more consistent and continuous bottom-simulating reflector (BSR) and better resolved sub-BSR reflectors. The results from the third method reveal more detailed internal structures of debris-flow deposits and increased continuity of reflectors; in addition, the seafloor reflection and the BSR appear to have changed to a zero-phase waveform. These properties help more precisely estimate the distribution and reserves of gas hydrates in the exploration area by improving analysis of facies and amplitude of the BSR.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.148-152
• /
• 2008

Multichannel seismic profiles reveal a strong bottom simulating reflector (BSR) occurring below the seafloor in the plain of the Ulleung Basin, East Sea (Japan Sea). The essential characteristics of the BSR include its cross-cutting relationship to strata, strong amplitude, and reverse polarity with respect to the seafloor reflection, representing the base of the gas hydrate stability zone (BHSZ). The BSR reflection coefficient ranging from -0.23 to -0.26 is 1.5${\sim}$1.7 times that of the seafloor reflection and interval velocities decrease to less than 700 m/s below the BSR. These features indicate the existence of free gas beneath the GHSZ. Heat flow, estimated from the BSR depth as $95{\sim}98mW/m^2$, is in good agreement with measured values. Therefore, the BSR can be efficiently used to estimate regional distribution of heat flow in the Ulleung Basin.

• Journal of the Korean Geophysical Society
• /
• v.7 no.2
• /
• pp.99-112
• /
• 2004

The Basin, a marginal basin located between the Antarctic Peninsula and the South Shetland Islands, is consist of three small basins, the Central, Eastern, Western Basins. Seismic data obtained on December 1995 show well-defined spreading ridges, basement highs, faults, morphology of the basin, distribution of sediments, crustal and sedimentary deformation, diapirs, and contourites. The main spreading axis of the Central Bransfield Basin connecting Deception and Bridgeman Islands continues up to the central part of the Eastern Basin, whereas deep basin covered by thick sediments without any spreading structures develops in the northeastern part. This indicates that back-arc spreading along the axis of the Bransfield Basin has been taken place in the southwestern part of the Eastern Basin, not in the northeastern part. Many NW-SE trending faults perpendicular to the axis of the basin would be related with strike-slip movement of the Shackleton Fracture. Zone. Extensinal strutures like deep basin without any spreading structures in the northeastern part, normal faults and diapirs on both continental slopes of the Eastern Basin would be formed by extension as a consequence of the sinistral movement between Antarctic and the Scotia plates.