• 자원환경지질
• /
• 제39권6호
• /
• pp.711-717
• /
• 2006

한국지질자원연구원은 1997년부터 새로운 에너지 자원으로 활용 가능성을 포함하고 있는 가스 하이드레이트를 조사하기 위해 동해 일원에서 탄성파탐사를 실시하고 있다. 탄성파 반사자료로부터 가스 하이드레이트 부존여부를 확인하는 방법은 해저면과 평행하면서 위상이 반대로 나타나는 고진폭 반사파 Bottom Simulating Reflector (BSR)과 BSR상부에서의 진폭감소, 하부에서 진폭증가와 구간속도 감소 둥을 들 수 있다. 대용량 탐사자료로 구성된 탄성파 반사자료에 깊이영역 구조보정을 적용하기 위해서는 고성능 컴퓨터와 병렬처리 기술이 필요하다. PSPI법은 적은 컴퓨터 계산량과 효율성 그리고 주파수 영역에서 구조적으로 병렬화가 용이한 특성을 지니고 있어 구조보정에 많이 이용되고 있다. 여기에서는 동해 가스 하이드레이트 탄성파 반사자료에 대한 일반자료처리와 함께 BSR로 여길 수 있는 구간에 대해 message passing interface_local area multicomputers(MPI_LAM)으로 병렬 코드화된 MPI PSPI를 이용하여 깊이영역 중합 전 구조보정에 적용하였다. 중합 전 깊이영역 구조보정 입력자료를 위한 속도모델은 자체 개발된 지오빗을 이용하여 중합 단면도로부터 지층경계면을 구하고 중합속도를 이용하여 제작하였다. BSR은 시간영역구조보정 된 중합 단면도상에서 음원모음도 3555-4162 사이와 왕복주시 2950 ms 부근에서 확인되지만 깊이영역 단면도에서는 해수면 6 km에서 17 km사이, 해저면에서 약 2.1km 깊이영역에서 나타남을 알 수 있다. 또한 구조보정 결과 반사파 에너지가 집중되는 지점에서 영상화가 잘 이루어지므로 관심대상 지역에 에너지를 많이 보낼 수 있는 자료취득변수를 결정해야 함을 알 수 있다.

• 지구물리와물리탐사
• /
• 제13권4호
• /
• pp.315-322
• /
• 2010

위상막 구조보정과 split-step Fourier 구조보정은 주파수-파수, 주파수-공간 영역에서 단방향 파동방정식을 이용하여 빠른 계산 속도로 수평적 속도변화를 고려할 수 있는 구조보정이다. 일반화된-막(generalized-screen) 구조보정은 주파수-파수영역에서 수직전파를 가정하는 위의 두 구조보정과는 달리 수직전파를 가정하지 않고, 지수함수의 무한급수 전개를 이용한다. 또한 수직느리기항의 테일러 급수전개를 일반화하여 고차항을 추가함으로써 급격한 속도변화를 갖는 지하구조에서 넓은 각으로 전파하는 파동장에 대한 정확도를 향상시켰다. 이 논문은 다양한 경사와 급격한 속도변화를 포함하는 복잡한 지하구조를 효율적으로 보다 정확하게 영상화하기 위하여 2차원 일반화된-막 구조보정에 대하여 연구하였다. 일정한 미소변량(constant perturbation)을 갖는 매질과 SEG/EAGE 암염돔을 모사한 모델에 대하여 일반화된-막 전파자와 위상막 전파자의 전파된 파동장을 비교한 결과, 일반화된-막 전파자가 파동장의 넓은각 전파에 대해 위상막 전파자보다 높은 정확도를 보였다. 또한 일반화된-막 전파자의 차수를 증가시킬수록 넓은 각으로 전파하는 파동장의 정확도가 향상되었다. 큰 수평적 속도변화와 급경사를 갖는 모델과 SEG/EAGE 암염돔 합성 탄성파탐사 자료에 대하여 일반화된-막 구조보정과 위상막 구조보정을 적용한 결과, 일반화된-막 구조보정이 속도변화가 크고 급격한 경사를 갖는 반사면을 보다 정확한 위치에 뚜렷하게 영상화하였다.

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

Refraction tomography requires an algorithm for efficiently computing the traveltimes and their $Fr\'{e}chet$ derivatives. We have attempted to solve the damped wave field using the frequency domain finite element model ing and then invoked the reciprocity theorem to calculate the $Fr\'{e}chet$ derivative of the traveltime with respect to the subsurface parameter. Then, we used a damped least square method to invert the traveltimes of the Marmousi 2 model. Numerical tests demonstrate that the refraction tomography with large aperture data can be used to estimate the smooth velocity model for the prestack depth migration.

• 지구물리와물리탐사
• /
• 제5권3호
• /
• pp.185-192
• /
• 2002

구조물의 안전진단을 위한 비파괴 검사 분야에 사용되는 지표레이다 자료를 정밀하게 영상화시켜줄 수 있는 3 차원 구조보정기법을 살펴보았다. 이 연구에서는 탄성파 자료처리를 위해 개발된 3차원 Kirchhoff 중합전 심도구조보정기법을 이용하여 구조물 안전진단시 중요대상체인 콘크리트 내의 철근, 균열 등을 3차원적으로 영상화시켰다. 철근과 공동대상체가 동시에 존재하는 경우 물성차이의 구분이 가능토록 영상화되었고, 여러 가지 크기의 철근이 존재하는 경우 XX배열법에 의한 영상이 YY배열법보다 더 높은 분해능을 보여주었다. 따라서 보다 정밀한 정보를 얻기 위해서는 일반적인 YY배열 외에 XX배열법의 영상화 결과를 연계하여 분석하는 것이 필요하다. 마지막으로 주향이 교차하는 두 개의 철근이 다른 심도에 존재하는 경우에는 상부의 철근은 정확하게 영상화되나 철근의 높은 전도도로 인하여 하부의 철근은 실제크기보다 좀 더 크게 영상화되는 한계를 보여주었다.

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

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

• 지구물리와물리탐사
• /
• 제25권2호
• /
• pp.71-84
• /
• 2022

탄성파 자료 취득 시 신호와 함께 기록되는 다양한 형태의 잡음은 탄성파 자료의 정확한 해석을 방해하는 요인으로 작용한다. 따라서 탄성파 자료의 잡음 제거는 탄성파 자료 처리 과정 중 필수적인 절차이므로 기계 학습을 포함한 다양한 방식의 잡음 제거 연구가 수행되고 있다. 본 연구에서는 비지도 학습 기반의 탄성파 잡음 제거 모델을 이용하여 중합 전 탄성파 자료의 잡음 제거를 수행하고자 하였으며 총 세 가지의 비지도 학습 기반 기계 학습 모델을 비교하였다. 세 가지의 비지도 학습 모델은 N2NUNET, PATCHUNET, DDUL로 각각 서로 다른 신경망 구조를 통해 정답 자료 없이 탄성파 잡음을 제거한다. 세 가지 모델들을 인공 합성 및 현장 중합 전 탄성파 자료에 적용하여 잡음을 제거한 후 그 결과를 정성적·정량적으로 분석하였으며, 분석 결과 세 가지 비지도 학습 모델 모두 인공 합성 및 현장 자료의 탄성파 잡음을 적절히 제거하였음을 확인하였다. 그 중 N2NUNET 모델이 가장 낮은 잡음 제거 성능을 보여주었으며, PATCHUNET과 DDUL은 거의 유사한 결과를 도출하였지만, DDUL이 정량적으로 근소한 우위를 보였다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2009년도 특별 심포지엄
• /
• pp.67-71
• /
• 2009

The high resolution seismic data processing essentially requires the latest data processing technology with a system which enables to calculate the large volume of input/output and computing of data. Nowadays, there is a high foreign dependency on the data processing and QC. In this research, The test processing had been carried out to deduce the best suited software by evaluating functionality, usability and user environment. And several attempts are made to select the most appropriated specification of hardware in order to maximize the software efficiency. The results showed that more numbers of multiple and noise were eliminated in existing prestack time migration sections, which effectively enhanced the resolution.

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

Reverse-time migration has no dip limitations and one of the most promising methods to preserve true amplitudes. We applied 3-D prestack reverse time migration based on a pseudo-spectral implementation of the acoustic wave equation to the SEG/EAGE salt dome synthetic data set. We were able to illuminate sub salt reflectors of the SEG/EAGE salt model that were barely observable in the Kirchhoff migration images. Using the pseudo-spectral modeling technique, we could implement reverse-time migration within the core memory, which could be equipped to a personal computer.