• Economic and Environmental Geology
• /
• v.56 no.4
• /
• pp.365-384
• /
• 2023

We applied the partitioned waveform inversion to 2,026 event data recorded at 173 seismic stations from the Incorporated Research Institutions for Seismology Data Managing Center and the Ocean Hemisphere network Project to estimate S-wave velocity and radial anisotropy models beneath the Western Pacific. In the Philippine Sea plate, high-Vs anomalies reach deeper in the West Philippine basin than in the Parece-Vela basin. Low-Vs anomalies found at 80 km below the Parece-Vela basin extend deeper into the West Philippine Basin. This velocity contrast between the basins may be caused by differences in lithospheric age. Low-Vs anomalies are observed beneath the Caroline seamount chain and the Caroline plate. Overall positive radial anisotropy anomalies are observed in the Western Pacific, but negative radial anisotropy is found at > 220 km depth on the subducting plate along the Mariana trench and at ~50 km in the Parece-Vela basin. Positive radial anisotropy is found at > 200 km depth beneath the Caroline seamount chain, which may indicate the 'drag' between the plume and the moving Pacific plate. High-Vs anomalies are found at 40 ~ 180 km depth beneath the Ontong-Java plateau, which may indicate the presence of unusually thick lithosphere due to underplating of dehydrated plume material.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.214-216
• /
• 2022

Distributed acoustic sensing (DAS), an increasingly growing acquisition technique in the oil and gas exploration and seismology fields, has been used to record seismic signals using optical cables as receivers. With the development of imaging methods for DAS data, full waveform inversion (FWI) is been applied to DAS data to obtain high-resolution property models such as P- and S-velocity. However, because the DAS systems measure strain from the phase distortion between two points along optical cables, DAS data must be transformed from strain to particle velocity for FWI algorithms. In this study, a plane-wave FWI algorithm based on the relationship between strain and horizontal particle velocity in the plane-wave assumption is proposed to apply FWI to DAS data. Under the plane-wave assumption, strain equals the horizontal particle velocity, which is scaled by the velocity at the receiver position. This relationship was confirmed using a numerical experiment. Furthermore, 4-layer and modified Marmousi-2 velocity models were used to verify the applicability of the proposed FWI algorithm in various survey environments. The proposed FWI was implemented in land and marine survey environments and provided high-resolution P- and S-velocity models.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.25 no.6
• /
• pp.583-603
• /
• 2023

Recently, the advancement of mechanical tunnel boring machine (TBM) technology and the characteristics of subsea railway tunnels subjected to hydrostatic pressure have led to the widespread application of shield TBM methods in the design and construction of subsea railway tunnels. Subsea railway tunnels are exposed in a constant pore water pressure and are influenced by the amplification of seismic waves during earthquake. In particular, seismic loads acting on subsea railway tunnels under various ground conditions such as soft ground, soft soil-rock composite ground, and fractured zones can cause significant changes in tunnel displacement and stress, thereby affecting tunnel safety. Additionally, the dynamic response of the ground and tunnel varies based on seismic load parameters such as frequency characteristics, seismic waveform, and peak acceleration, adding complexity to the behavior of the ground-tunnel structure system. In this study, a finite difference method is employed to model the entire ground-tunnel structure system, considering hydrostatic pressure, for the investigation of dynamic behavior of subsea railway tunnel during earthquake. Since the key factors influencing the dynamic behavior during seismic events are ground conditions and seismic waves, six analysis cases are established based on virtual ground conditions: Case-1 with weathered soil, Case-2 with hard rock, Case-3 with a composite ground of soil and hard rock in the tunnel longitudinal direction, Case-4 with the tunnel passing through a narrow fault zone, Case-5 with a composite ground of soft soil and hard rock in the tunnel longitudinal direction, and Case-6 with the tunnel passing through a wide fractured zone. As a result, horizontal displacements due to earthquakes tend to increase with an increase in ground stiffness, however, the displacements tend to be restrained due to the confining effects of the ground and the rigid shield segments. On the contrary, peak compressive stress of segment significantly increases with weaker ground stiffness and the effects of displacement restrain contribute the increase of peak compressive stress of segment.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.18 no.4
• /
• pp.59-66
• /
• 2010

Aerospace structures need high stability and long life because many personal injuries can result from an accident and securing structural integrity for various external environments is more important than any other thing. So first of all we must prove the destruction properties for operating environment, have prediction technology about damage evolution and life, and develop an economical non-destructive technology capable of detecting structure damage. Acoustic emission (AE) have no need of artificial environment like ultrasonic inspection or radio fluoroscopy to emit a certain energy, is a testing technique using seismic signal resulting from interior changes of solids, and enables to observe if any fault is appeared and it grows seriously or not while running. In this study we suggest the method of structural integrity evaluation for aerospace structures through the acoustic emission technique, for which a model plane was manufactured and an actual operation test was conducted.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.136-144
• /
• 2016

Since seismic inversion is based on the wave equation, it is important to calculate the solution of wave equation exactly. In particular, full waveform inversion would produce reliable results only when the forward modeling is accurately performed because it uses full waveform. When we use finite-difference or finite-element method to solve the wave equation, the convergence of numerical scheme should be guaranteed. Although the general proof of convergence is provided theoretically, the consistency and stability of numerical schemes should be verified for practical applications. The implementation of source function is the most crucial factor for the consistency of modeling schemes. While we have to use the sinc function normalized by grid spacing to correctly describe the Dirac delta function in the finite-difference method, we can simply use the value of basis function, regardless of grid spacing, to implement the Dirac delta function in the finite-element method. If we use frequency-domain wave equation, we need to use a conservative criterion to determine both sampling interval and maximum frequency for the source wavelet generation. In addition, the source wavelet should be attenuated before applying it for modeling in order to make it obey damped wave equation in case of using complex angular frequency. With these conditions satisfied, we can develop reliable inversion algorithms.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.8-18
• /
• 1998

In exploration seismology, the Kirchhoff hyperbola has been successfully used to migrate reflection seismo-grams. The mathematical basis of Kirchhoff hyperbola has not been clearly defined and understood for the application of prestack or poststack migration. The travel time from the scatterer in the subsurface to the receivers (exploding reflector model) on the surface can be a kinematic approximation of Green's function when the source is excited at position of the scatterer. If we add the travel time from the source to the scatterer in the subsurface to the travel time of exploding reflector model, we can view this travel time as a kinematic approximation of the partial derivative wavefield with respect to the velocity or the density in the subsurface. The summation of reflection seismogram along the Kirchhoff hyperbola can be evaluated as an inner product between the partial derivative wavefield and the field reflection seismogram. In addition to this kinematic interpretation of Kirchhoff hyperbola, when we extend this concept to shallow refraction seismic data, the stacking of refraction data along the straight line can be interpreted as a measurement of an inner product between the first arrival waveform of the partial derivative wavefield and the field refraction data. We evaluated the Kirchhoff hyperbola and the straight line for stacking the refraction data in terms of the first arrival waveform of the partial derivative wavefield with respect to the velocity or the density in the subsurface. This evaluation provides a firm and solid basis for the conventional Kirchhoff migration and the straight line stacking of the refraction data.

• Geophysics and Geophysical Exploration
• /
• v.24 no.1
• /
• pp.6-19
• /
• 2021

Subsurface physical properties can be obtained and imaged by seismic exploration, and various algorithms have been developed for this purpose. In this regard, root mean square error (RMSE) has been widely used to quantitatively evaluate the accuracy of the developed algorithms. Although RMSE has the advantage of being numerically simple, it has limitations in assessing structural similarity. To supplement this, full-reference image quality assessment (FR-IQA) techniques, which reflect the human visual system, are being investigated. Therefore, we selected six FR-IQA techniques that could evaluate the obtained physical properties. In this paper, we used the full-waveform inversion, because the algorithm can provide the physical properties. The inversion results were applied to the six selected FR-IQA techniques using three benchmark models. Using salt models, it was confirmed that the inversion results were not satisfactory in some aspects, but the value of RMSE decreased. On the other hand, some FR-IQA techniques could definitely improve the evaluation.

• Journal of the Korean earth science society
• /
• v.38 no.7
• /
• pp.561-569
• /
• 2017

Fault plane solutions of the Iksan earthquake ($M_L=3.85$) and two aftershocks were obtained using the FOCMEC (FOCal MEChanism determination) program. The main event showed the characteristic of strike slip faulting with reverse component. It has the fault planes with NE-SW or NW-SE direction. This is similar to the fault characteristics of earthquake pattern in the inland area of the Korean Peninsula. In order to detect micro-earthquake events, continuous seismic waveform data of the thirteen seismic stations within a radius of 100km from epicenter were analyzed by PQLII program (PASSCAL, 2017) for the period from December 15, 2015 to January 22, 2016. The epicenters of nineteen micro-events were newly determined by Hypoinverse-2000 program. They are not concentrated along some lineaments or fault lines. The intensity of the Iksan earthquake was obtained by estimating the telephone inquiries, the degree of ground shaking or damage all around the southern peninsula. The instrumental intensity was also obtained using PGA (Peak Ground Acceleration) records. As a result, the maximum MM intensity was estimated to be V near the epicenter.

• 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.20 no.3
• /
• pp.129-136
• /
• 2017

3D seismic data processing methods such as full waveform inversion or reverse-time migration require 3D wave propagation modeling and heavy calculations. We compared efficiency and accuracy of a Xeon Phi coprocessor to those of a high-end server CPU using 3D frequency-domain wave propagation modeling. We adopted the OpenMP parallel programming to the time-domain finite difference algorithm by considering the characteristics of the Xeon Phi coprocessors. We applied the Fourier transform using a running-integration to obtain the frequency-domain wavefield. A numerical test on frequency-domain wavefield modeling was performed using the 3D SEG/EAGE salt velocity model. Consequently, we could obtain an accurate frequency-domain wavefield and attain a 1.44x speedup using the Xeon Phi coprocessor compared to the CPU.