• 한국신재생에너지학회:학술대회논문집
• /
• 2007.06a
• /
• pp.497-500
• /
• 2007

탄성파 코다 파는 두 수진기에서 기록된 탄성파 자료의 상호상관으로부터 두 신호에 대한 순간응답을 구하고 이로부터 지층정보를 구하는데 이용된다. 여기에서는 인공합성 탄성파 자료와 가스 하이드레이트 현장자료에 적용하여 상호상관 모음도와 가상음원 모음도 (virtual source)를 구하고자 하였다. 인공합성자료는 해저면 탄성파 탐사법 (ocean bottom seismic)을 모델로 이용하여 인공합성 탄성파 단면도를 제작하였으며, 탄성파 코다 파를 살펴보기 위해 인공 OBS 자료 중 첫 번째 트레이스를 가상음원으로 정하고 모든 음원 모음도와 상호상관으로 가상응원 단면도를 제작하였다. 현장자료 적용으로는 해저면 기인 고진폭 반사파인 BSR (bottom simulating reflection)을 포함하고 있는 자료를 선정하여 상호상관 단면도와 가상음원 단면도를 제작하였다. 중합단면도상에 나타난 가스 분출지역은 상호상관 단면도에서도 나타났으며, 중합단면도상 BSR부분은 vs 단면도에서 강한 반사파를 보여줌을 알 수 있었다.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.457-465
• /
• 2013

We integrated and correlated datasets from surface and subsurface geophysics, drilling cores, and engineering geology to identify geological interfaces and characterize the joints and fracture zones within the rock mass. The regional geometry of a geologically weak zone was investigated via a fence projection of electrical resistivity data and a borehole image-processing system. Subsurface discontinuities and intensive fracture zones within the rock mass are delineated by cross-hole seismic tomography and analyses of dip directions in rose diagrams. The dynamic elastic modulus is studied in terms of the P-wave velocity and Poisson's ratio. Subsurface discontinuities, which are conventionally identified using the N value and from core samples, can now be identified from anomalous reflection coefficients (i.e., acoustic impedance contrast) calculated using a pair of well logs, comprising seismic velocity from suspension-PS logging and density from logging. Intensive fracture zones identified in the synthetic seismogram are matched to core loss zones in the drilling core data and to a high concentration of joints in the borehole imaging system. The upper boundaries of fracture zones are correlated to strongly negative amplitude in the synthetic trace, which is constructed by convolution of the optimal Ricker wavelet with a reflection coefficient. The standard deviations of dynamic elastic moduli are higher for fracture zones than for acompact rock mass, due to the wide range of velocities resulting from the large numbers of joints and fractures within the zone.

• Geophysics and Geophysical Exploration
• /
• v.25 no.2
• /
• pp.71-84
• /
• 2022

When acquiring seismic data, various types of simultaneously recorded seismic noise hinder accurate interpretation. Therefore, it is essential to attenuate this noise during the processing of seismic data and research on seismic noise attenuation. For this purpose, machine learning is extensively used. This study attempts to attenuate noise in prestack seismic data using unsupervised machine learning. Three unsupervised machine learning models, N2NUNET, PATCHUNET, and DDUL, are trained and applied to synthetic and field prestack seismic data to attenuate the noise and leave clean seismic data. The results are qualitatively and quantitatively analyzed and demonstrated that all three unsupervised learning models succeeded in removing seismic noise from both synthetic and field data. Of the three, the N2NUNET model performed the worst, and the PATCHUNET and DDUL models produced almost identical results, although the DDUL model performed slightly better.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.

• Geophysics and Geophysical Exploration
• /
• v.27 no.3
• /
• pp.171-180
• /
• 2024

The processing of seismic data involves analyzing earthquake wave data to understand the internal structure and characteristics of the Earth, which requires high computational power. Recently, machine learning (ML) techniques have been introduced to address these challenges and have been utilized in various tasks such as noise reduction and velocity model construction. However, most studies have focused on specific seismic data processing tasks, limiting the full utilization of similar features and structures inherent in the datasets. In this study, we compared the efficacy of using receiver-wise time-series data ("receiver array") and synchronized receiver signals ("time array") from shotgathers for pretraining a Bidirectional Encoder Representations from Transformers (BERT) model. To this end, shotgather data generated from a synthetic model containing faults was used to perform noise reduction, velocity prediction, and fault detection tasks. In the task of random noise reduction, both the receiver and time arrays showed good performance. However, for tasks requiring the identification of spatial distributions, such as velocity estimation and fault detection, the results from the time array were superior.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.50-57
• /
• 2011

For effective and accurate prediction of overpressure in the Efomeh field, located in the Niger delta basin of Nigeria, integrated seismic and borehole analyses were undertaken. Normal and abnormal pore pressure zones were delineated based on the principle of normal and deviation from normal velocity trends. The transition between the two trends signifies the top of overpressure. The overpressure tops were picked at regular intervals from seismic data using interval velocities obtained by applying Dix's approximation. The accuracy of the predicted overpressure zone was confirmed from the sonic velocity data of the Efomeh 01 well. The variation to the depth of overpressure between the predicted and observed values was less than 10mat the Efomeh 01 well location, with confidence of over 99 per cent. The depth map generated shows that the depth distribution to the top of the overpressure zone of the Efomeh field falls within the sub-sea depth range of 2655${\pm}$2m (2550 ms) to 3720${\pm}$2m (2900 ms). This depth conforms to thick marine shales using the Efomeh 01 composite log. The lower part of the Agbada Formation within the Efomeh field is overpressured and the depth of the top of the overpressure does not follow any time-stratigraphic boundary across the field. Prediction of the top of the overpressure zone within the Efomeh field for potential wells that will total depth beyond 2440m sub-sea is very important for safer drilling practice as well as the prevention of lost circulation.

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

Polarization filtering based on the multicomponent complex trace analysis method is performed to reject Rayleigh waves and to enhance S waves. To test the polarization filter, synthetic seismic data were constructed for a simple two-layer model based on the finite difference method. Rayleigh waves with elliptic motion are eliminated effectively and P and S waves with linear motions are well separated each other.

• Journal of the Korean Geophysical Society
• /
• v.8 no.1
• /
• pp.35-38
• /
• 2005

Polarization filtering based on the multicomponent complex trace analysis method is performed to reject Rayleigh waves and to enhance S waves. To test the polarization filter, synthetic seismic data were constructed for a simple two-layer model based on the finite difference method. Rayleigh waves with elliptic motion are eliminated effectively and P and S waves with linear motions are well separated each other.

• Economic and Environmental Geology
• /
• v.27 no.6
• /
• pp.557-561
• /
• 1994

Q estimates are made by direct measurements of energy loss per cycle from primary P and S waves, as a function of frequency. Assuming that intrinsic Q is frequency independent and scattering Q is frequency dependent over the frequencies of interest, the relative contributions of each, to a total observed Q, may be estimated. Test examples are produced by computing viscoelastic synthetic seismograms using a pseudo spectral solution with inclusion of relaxation mechanisms (for intrinsic Q) and a fractal distribution of scatterers (for scattering Q). The composite theory implies that when the total Q for S-waves is smaller than that for P-waves (the usual situation), intrinsic Q is dominating; when it is larger, scattering Q is dominating. In the inverse problem, performed by a global least squares search, intrinsic $Q_p$ and $Q_s$ estimates are reliable and unique when their absolute values are sufficiently low that their effects are measurable in the data. Large $Q_p$ and $Q_s$ have no measurable effect and hence are not resolvable. Standard deviation of velocity $({\sigma})$ and scatterer size (A) are less unique as they exhibit a tradeoff as predicted by Blair's equation. For the P-waves, intrinsic and scattering contributions are of approximately the same importance, for S-waves, the intrinsic contributions dominate.

• Journal of the Korean Geophysical Society
• /
• v.3 no.1
• /
• pp.37-48
• /
• 2000

In order to investigate in-line ground motions caused by earthquakes, we examine the multicomponent complex trace analysis method (MCTAM) for the synthetic data and apply it to real earthquake data. An experimental result for synthetic data gives correct information on the arrival times, duration of individual phases, and approaching angles for body waves. Rayleigh waves are also easily identified with the MCTAM. A deep earthquake with magnitude of 7.3 was chosen to test various polarization attributes of ground motions. For P waves, instantaneous phase difference between the vertical and the in-line horizontal components ${\phi}(t)$, instantaneous reciprocal ellipticity ${\rho}(t)$, and approaching angle ${\tau}(t)$ are computed to be ${\pm}180^{\circ},\;0{\sim}0.25,\;and\;-30^{\circ}{\sim}-45^{\circ}$, respectively. For S waves, ${\phi}(t)$ tends to vary while ${\rho}(t)$ have values of $0{\sim}0.3\;and\;{\tau}(t)$ remains near vertical, respectively. A relatively low frequency signal registered just prior to the S wave event is interpreted as a P-wave phase based on its polarization characteristics. Velocities of P and S waves are computed to be 8.633 km/s and 4.762 km/s, and their raypath parameters 0.074 s/km and 0.197 s/km. Dynamic Poisson's ratio is obtained as 0.281 from the velocities of P and S waves.