• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.459-463
• /
• 2003

Waveform inversion requires extracting a reliable low frequency content of seismic data for estimating of the low wave number velocity model. The low frequency content of the seismic data is usually discarded or neglected because of the band-limited response of the source and the receivers. In this study, however small the spectral of the low frequency seismic data is, we assume that it is possible to extract a reliable phase information of the low frequency from the seismic data and use it in waveform inversion. To this end, we exploit the frequency domain finite element modeling and source-receiver reciprocity to calculate the $Frech\`{e}t$ derivative of the phase of the seismic data with respect to the earth model parameter such as velocity, and then apply a damped least squares method to invert the phase of the seismic data. Through numerical example, we will attempt to demonstrate the feasibility of our method in estimating the correct velocity model for prestack depth migration.

• 대한자원환경지질학회:학술대회논문집
• /
• 대한자원환경지질학회 2002년도 춘계 공동학술발표회
• /
• pp.159-162
• /
• 2002

We tried to obtain an initial velocity model for prestack depth migration via waveform inversion. For application of any field data we chose a smooth background layered velocity model (v=v0 + k x z) as an initial velocity model. Newton type waveform inversion needs to invert huge Hessian matrix. In order to compute full Hessian matrix arising from full aperture data and full illumination zone, we meet insurmountable difficulties of paying astronomical computing cost. For the layered media, approximate Hessian emerging from single shot aperture data can be used repeatedly for split spread source configuration. In our work of using this Hessian characteristic of layered media we attempted to obtain the approximate velocity model as close as possible to the true velocity model in first iteration.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2005년도 춘계학술대회
• /
• pp.638-642
• /
• 2005

For gas hydrate exploration, long offset multichannel seismic data acquired using by the 4km streamer length in Ulleung basin of the East Sea. The dataset was processed to define the BSRs (Bottom Simulating Reflectors) and to estimate the amount of gas hydrates. Confirmation of the presence of Bottom Simulating reflectors (BSR) and investigation of its physical properties from seismic section are important for gas hydrate detection. Specially, faster interval velocity overlying slower interval velocity indicates the likely presences of gas hydrate above BSR and free gas underneath BSR. In consequence, estimation of correct interval velocities and analysis of their spatial variations are critical processes for gas hydrate detection using seismic reflection data. Using Dix's equation, Root Mean Square (RMS) velocities can be converted into interval velocities. However, it is not a proper way to investigate interval velocities above and below BSR considering the fact that RMS velocities have poor resolution and correctness and the assumption that interval velocities increase along the depth. Therefore, we incorporated Migration Velocity Analysis (MVA) software produced by Landmark CO. to estimate correct interval velocities in detail. MVA is a process to yield velocities of sediments between layers using Common Mid Point (CMP) gathered seismic data. The CMP gathered data for MVA should be produced after basic processing steps to enhance the signal to noise ratio of the first reflections. Prestack depth migrated section is produced using interval velocities and interval velocities are key parameters governing qualities of prestack depth migration section. Correctness of interval velocities can be examined by the presence of Residual Move Out (RMO) on CMP gathered data. If there is no RMO, peaks of primary reflection events are flat in horizontal direction for all offsets of Common Reflection Point (CRP) gathers and it proves that prestack depth migration is done with correct velocity field. Used method in this study, Tomographic inversion needs two initial input data. One is the dataset obtained from the results of preprocessing by removing multiples and noise and stacked partially. The other is the depth domain velocity model build by smoothing and editing the interval velocity converted from RMS velocity. After the three times iteration of tomography inversion, Optimum interval velocity field can be fixed. The conclusion of this study as follow, the final Interval velocity around the BSR decreased to 1400 m/s from 2500 m/s abruptly. BSR is showed about 200m depth under the seabottom

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2006년도 공동학술대회 논문집
• /
• pp.131-135
• /
• 2006

본 논문에서는 가우스 뉴튼법을 이용한 중합전 탄성파 자료의 파형역산에 관한 연구를 수행하였다. 탄성파 파형역산에 가우스 뉴튼법을 적용하는 방법은 80년대에 제시되었으나 최근 들어서야 활발히 연구가 진행되고 있는데 이는 연산 능력과 기억용량의 한계에 기인한 것이다. 이를 극복하기 위해 본 연구에서는, 파동 전파 수치모의와 역산과정에서 각각 다른 크기의 격자간격을 사용하고, 필요한 시간영역의 파동전파 모사와 가상 진원의 근사를 통해 편미분 파형을 계산하였으며, 효과적으로 슈퍼컴퓨터를 활용하기 위해 병렬처리 기법을 사용하였다. 수치모의를 통해, 가우스 뉴튼법을 이용한 파형 역산의 수렴속도가 빠르고 정확한 것을 알 수 있었으며, 이를 통해 본 연구에서 제시한 방법의 실제 탄성파 자료를 이용한 역산에의 적용가능성을 확인하였다.

• 대한자원환경지질학회:학술대회논문집
• /
• 대한자원환경지질학회 2002년도 춘계 공동학술발표회
• /
• pp.163-166
• /
• 2002

Initial velocity model close to real velocity structure of the subsurface governs the quality of image of prestack depth migration. Geophysicists employ velocity estimation tools such as velocity analysis (curvature method, coherency inversion), tomography and waveform inversion. We present a reflection tomography that parameterizes the subsurface into the movable blocks. By interpreting the depth-migrated section or stacked section, we can design an initial constant velocity model having only impedance boundaries. We use shooting-raytracing method that allows us to calculate the Jacobian-matrix efficiently.

• 자원환경지질
• /
• 제49권3호
• /
• pp.225-242
• /
• 2016

암석 물리(rock physics)는 암석의 물성과 탄성파 자료의 정량적인 연결 고리로서 저류암 특성화와 모니터링 등에 적용되는 중요한 도구이다. 암석 물리를 기반으로 시추공 검층 자료로부터 공극 유체와 구성 성분이 다양한 암석의 성질을 대표하는 유효 탄성 계수(effective elastic constants)를 모델링하여 탄성파 자료에 적용함으로서 탄성파 자료전체에 대한 저류암 분포와 공극 유체를 유추할 수 있다. 본 보고에서는 먼저 Voigt, Reuss, Hashin-Shtrikman의 유효 탄성 계수 한계값과 Gassmann 방정식을 이용한 유체 치환에 대해서 설명한 후에 공극률이 비교적 높은 저류 사암에 널리 적용되는 접촉(contact) 암석 물리 모델들을 소개했다. 그리고 접촉 암석 물리 모델 중에서 분포 깊이가 어느 정도 일정하고 공극률이 비교적 높은 사암에 적합한 것으로 알려진 일정 교결량(constant-cement) 모델을 북해 유전의 검층 자료에 적용하여 암석 물리 견본(rock physics template)을 완성했다. 마지막으로 이 암석 물리 견본과 현장 2D 탄성파 자료의 중합전 역산 결과로부터 이 유전의 저류암과 덮개암, 공극 유체 분포를 예측해 보았다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.402-409
• /
• 2003

Inversion of traveltime requires an efficient algorithm for computing the traveltime as well as its $Frech\hat{e}t$ derivative. We compute the traveltime of the head waves using the damped wave solution in the Laplace domain and then present a new algorithm for calculating the $Frech\hat{e}t$ derivative of the head wave traveltimes by exploiting the numerical structure of the finite element method, the modem sparse matrix technology, and SWEET algorithm developed recently. Then, we use a properly regularized steepest descent method to invert the traveltime of the Marmousi-2 model. Through our numerical tests, we will demonstrate that the refraction tomography with large aperture data can be used to construct the initial velocity model for the prestack depth migration.