• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.105-115
• /
• 2011

Although there are many possible mechanisms for the intrinsic seismic attenuation in composite materials that include fluids, relative motion between solids and fluids during seismic wave propagation is one of the most important attenuation mechanisms. In our previous study, we conducted ultrasonic wave transmission measurements on an ice-brine coexisting system to examine the influence on ultrasonic waves of the unfrozen brine in the pore microstructure of ice. In order to elucidate the physical mechanism responsible for ultrasonic wave attenuation in the frequency range of 350.600 kHz, measured at different temperatures in partially frozen brines, we employed a poroelastic model based on the Biot theory to describe the propagation of ultrasonic waves through partially frozen brines. By assuming that the solid phase is ice and the liquid phase is the unfrozen brine, fluid properties measured by a pulsed nuclear magnetic resonance technique were used to calculate porosities at different temperatures. The computed intrinsic attenuation at 500 kHz cannot completely predict the measured attenuation results from the experimental study in an ice-brine coexisting system, which suggests that other attenuation mechanisms such as the squirt-flow mechanism and wave scattering effect should be taken into account.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.18-24
• /
• 2011

A multi-channel seismic reflection (MCS) survey was conducted in 2009 to explore the deep crustal structure of the Pacific Plate south of Hokkaido. The survey line happened to traverse a 250-km-wide Warm Core Ring (WCR), a current eddy that had been generated by the Kuroshio Extension. We attempted to use these MCS data to delineate the WCR fine structure. The survey line consists of two profiles: one with a shot interval of 200m and the other with a shot interval of 50 m. Records from the denser shot point line show much higher background noise than the records from the sparser shot point line. We identified the origin of this noise as acoustic reverberations between the sea surface, seafloor and subsurface discontinuities, from previous shots. Results showed that a prestack migration technique could enhance the signal buried in this background noise efficiently, if the sound speed information acquired from concurrent temperature measurements is available. The WCR is acoustically an assemblage of concave reflectors dipping inward, with steeper slopes (${\sim}2^{\circ}$) on th ocean side and gentler slopes (${\sim}1^{\circ}$) on the coastal side. Within the WCR, we recognised a 30-km-wide lens-shaped structure with reflectors on the perimeter.

• 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.20 no.1
• /
• pp.1-11
• /
• 2017

In order to improve the accuracy of microseismic epicenter location through the inversion techniques using P and S wave first arrivals, field experiments of microseismic monitoring were performed using borehole 3-component geophones. The direction of epicenter was estimated from the hodograms of P-wave first arrivals through the weight drop experiments in which the $\times$ component of 3-component geophone was aligned to the magnetic north. The picking of S wave first arrival was possible in the polarization filtered data even if S waves are difficult to be identified in raw data. The inversion technique using only P wave first arrival times can often converge to the local minimum when the initial values for epicenter are largely apart from the true epicenter, so that the correct solution can not be found. To solve this problem, the epicenter determination method using differences between P and S wave arrival times was used to estimate proper initial values of epicenter. The inversion result using only P-wave first arrival times which started from the estimated initial values showed the improved accuracy of the epicenter location.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.33-42
• /
• 2017

Recently, the grounded electrical-source airborne transient electromagnetic (GREATEM) system with high power source was introduced to achieve deeper investigation depth and to overcome high noise level. Although the GREATEM is a transient electromagnetic system using a long grounded wire as the transmitter, GREATEM data have been interpreted with 1D earth models because 2D or 3D modeling and inversion of vast airborne data are complicated and expensive to calculate. Generally, 1D inversion is subsequently applied to every survey point and combining 1D images together forms the stitched conductivity-depth image. However, the stitched models often result in abrupt variations in neighboring models. To overcome this problem, laterally constrained inversion (LCI) has been developed in inversion of ATEM data, which can yield layered sections with lateral smooth transitions. In this study, we analysed the GREATEM data through 1D numerical modeling for a curved grounded wire source. Furthermore, we developed a laterally constrained inversion scheme for continuous GREATEM data based on a layered earth model. All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. Applying the developed LCI technique to the GREATEM data, it was confirmed that the laterally constrained inversion can provide laterally smooth model sections that reflect the layering of the survey area effectively.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.25-32
• /
• 2017

We present a 3D SH-wave velocity model of the crust and uppermost mantle and seismic radial anisotropy beneath East Asia. The SH-wave velocity structure model was built using Love-wave group-velocity dispersion data from earthquake data recorded at broadband seismic networks of Korea, Japan, and China. Love-wave group-velocity dispersion curves were obtained by using the multiple filtering technique in the period range of 3 to 150 s for 3,369 event-station pairs. The inverted model using these data sets provides a crust and upper mantle SH-wave velocity structure down to 100 km depth. At 10 ~ 40 km depths SH-wave velocity beneath the East Sea is higher than beneath the Japanese island region. We estimated the Moho beneath the East Sea to be between 10 ~ 20 km depth, while Moho beneath the Korean Peninsula at around 35 km based on the depth where high-velocity anomalies are detected. We estimated the lithosphere-asthenosphere boundary beneath the East Sea to be at around 50 km based on the depth where strong low-velocity anomalies are observed. Widespread low-velocity anomalies are found between 50 ~ 100 km depth in the study region. Positive radial anisotropy ($V_{SV}$ > $V _{SH}$) is observed down to 35 km depth, while negative radial anisotropy ($V_{SV}$ > $V _{SH}$) is observed for deeper depth.

• Geophysics and Geophysical Exploration
• /
• v.23 no.3
• /
• pp.168-177
• /
• 2020

Seabed multiple reflections (seabed multiples) are the main cause of misinterpretations of primary reflections in both shot gathers and stack sections. Accordingly, seabed multiples need to be suppressed throughout data processing. Conventional model-driven methods, such as prediction-error deconvolution, Radon filtering, and data-driven methods, such as the surface-related multiple elimination technique, have been used to attenuate multiple reflections. However, the vast majority of processing workflows require time-consuming steps when testing and selecting the processing parameters in addition to computational power and skilled data-processing techniques. To attenuate seabed multiples in seismic reflection data, input gathers with seabed multiples and label gathers without seabed multiples were generated via numerical modeling using the Marmousi2 velocity structure. The training data consisted of normal-moveout-corrected common midpoint gathers fed into a U-Net neural network. The well-trained model was found to effectively attenuate the seabed multiples according to the image similarity between the prediction result and the target data, and demonstrated good applicability to field data.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.157-164
• /
• 2005

Urban conditions such as underground facilities and ambient noises due to cultural activity restrict the application of conventional geophysical techniques in general. We used the refraction microtremor (REMI) technique as an alternative way to get the geotechnical information, in particular shear-wave (S-wave) velocity information, at a site along an existing rail road. The REMI method uses ambient noises recorded using standard refraction equipment to derived shear-wave velocity information at a site. It does a wavefield transformation on the recorded wavefield to produce Rayleigh wave dispersion curve, which are then picked and modeled to get the shear-wave velocity structure. At this site the vibrations from the running trains provided strong noise sources that allowed REMI to be very effective. REMI was performed along the planned new underground rail tunnel. In addition, Suspension PS logging (SPS) were carried out at selected boreholes along the profile in order to draw out the quantitative relation between the shear wave velocity from the PS logging and the rock mass rating (RMR) determined from the inspection of the cores recovered from the same boreholes, These correlations were then used to relate the shear-wave velocity derived from REMI to RMR along the entire profile. The correlation between shear wave velocity and RMR was very good and so it was possible to estimate the RMR of the total zone of interest for the design of underground tunnel,

• Geophysics and Geophysical Exploration
• /
• v.6 no.1
• /
• pp.13-22
• /
• 2003

TSP (Tunnel Seismic Profiling) survey is a technique for imaging and characterizing geological structures ahead of a tunnel face. The seismic modeling algorithm and the synthetic data could be helpful for TSP surveys. However, there is few algorithm to describe the propagation of the elastic waves around the tunnel. In this study, existing 2-dimensional seismic modeling algorithm using finite element method was modified to make a suitable algorithm for TSP modeling. Using this algorithm, TSP modeling was practiced in some models. And the synthetic data was analyzed to examine the propagation characteristics of the elastic waves. First of all, the modeling for the homogeneous tunnel model was practiced to examine the propagation characteristics of the direct waves in the vicinity of the tunnel. And the algorithm was applied to some models having reflector which is perpendicular or parallel to the excavation direction. From these, the propagation characteristics of the reflected waves were examined. Furthermore, two source-receiver arrays were used in respective models to investigate the properties of the two arrays. These modeling algorithm and synthetic data could be helpful in interpreting TSP survey data, developing inversion algorithm and designing new source-receiver arrays.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.60-68
• /
• 2007

Conventional analysis of spectral induced polarization (SIP) data consists of measuring impedances over a range of frequencies, followed by spectral analysis to estimate spectral parameters. For the quantitative and accurate estimation of subsurface SIP parameter distribution, however, a sophisticated and stable inversion technique is required. In this study, we have developed a two-step inversion approach to obtain the two-dimensional distribution of SIP parameters. In the first inversion step, all the SIP data measured over a range of frequencies are simultaneously inverted, adopting cross regularisation of model complex resistivities at each frequency. The cross regularisation makes it possible to enhance the noise characteristics of the inversion by imposing a strong assumption, that complex resistivities should show similar characteristics over a range of frequencies. In numerical experiments, we could verify that our inversion approach successfully reduced inversion artefacts. As a second step, we have also developed an inversion algorithm to obtain SIP parameters based on the Cole-Cole model, in which frequency-dependent complex resistivities from the first step are inverted to obtain a two-dimensional distribution of SIP parameters. In numerical tests, the SIP parameter images showed a fairly good match with the exact model, which suggests that SIP imaging can provide a very useful subsurface image to complement resistivity.