• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2010년도 학술대회 초록집
• /
• pp.83-86
• /
• 2010

• 자원환경지질
• /
• 제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 깊이영역에서 나타남을 알 수 있다. 또한 구조보정 결과 반사파 에너지가 집중되는 지점에서 영상화가 잘 이루어지므로 관심대상 지역에 에너지를 많이 보낼 수 있는 자료취득변수를 결정해야 함을 알 수 있다.

• 지구물리와물리탐사
• /
• 제16권2호
• /
• pp.71-78
• /
• 2013

일방향 파동방정식(one-way wave equation)을 이용한 참반사 보정은 거꿀 참반사 보정(reverse time migration)에 비해서는 계산시간의 효율성을 향상시킬 수 있고 파선이론에 기초한 참반사 보정에 비해서는 상대적으로 정확한 영상화를 수행할 수 있다는 장점이 있다. 이 연구에서는 이러한 일방향 파동방정식을 이용하는 참반사 보정 알고리즘 중 Generalized-Screen (GS) 전파자를 이용하여 이방성 매질 중 횡적등방성(Vertical Transversely Isotropy, VTI) 매질에 적용가능한 겹쌓기 전 심도 참반사 보정(PreStack Depth Migration, PSDM) 모듈을 개발하였다. GS 전파자는 위상막(phase screen) 참반사 보정의 얇은 막(thin slab) 전파자에서 수직느리기항(vertical slowness)의 Taylor 급수 전개 시 고차항을 고려했기 때문에 수평적 속도변화가 크거나 급격한 경사를 지니는 구조에 대해서도 정확한 영상을 얻을 수 있다. 개발된 GS 참반사 보정 모듈을 검증하기 위해서 VTI 매질에 적용되는 전파자(GSVTI propagator)를 고차항의 차수에 따라 분석한 결과 차수가 증가할수록 넓은 각으로 전파하는 파동장의 정확도가 향상되었다. 그리고 수치모형실험을 통해서 이방성 매질에서 얻은 합성탄성파탐사 자료에 본 연구에서 개발된 이방성 GS 참반사 보정과 등방성 GS 참반사 보정을 각각 적용한 결과 이방성 GS 참반사 보정이 급경사나 이방성이 강한 지역에서의 구조에 대해 더욱 정확한 영상을 제공하는 것을 확인할 수 있었다.

• 대한자원환경지질학회:학술대회논문집
• /
• 대한자원환경지질학회 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.

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

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 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권1호
• /
• pp.18-22
• /
• 2002

단일 도달주시로서 최대에너지 도달주시는 카르히호프 구조보정을 위한 가장 적합한 연산자로 알려져 있다. 단순구조모델의 경우에 초동주시나 최대에너지 도달주시는 모두 정확한 카르히호프 구조보정 연산자로 사용될 수 있다. 그러나 구조가 복잡하고 속도대비가 큰 모델에 대해서 초동주시를 이용한 구조보정은 정확한 심도단면도를 얻을 수 없다. 즉 카르히호프 구조보정에서 요구되는 주시는 최대에너지 도달주시이지만 이 주시는 초동주시에 비해 상당히 많은 계산시간을 요한다. 본 논문에서는 일방향 파동방정식을 이용하여 최대에너지 도달주시에 근사한 주시를 계산하는 방법을 제안한다. 일방향 파동 방정식을 이용해 초동주시를 계산할 때 사용되는 WAS(Wrap Around Suppression, 두루마리현상억제)계수를 수평방향의 격자간격과 깊이에 대한 함수로 정의하고 음원 파형의 지연시간을 고려하여 최대 에너지 도달 시간에 근사한 주시를 구하였다. 이 주시의 타당성을 검증하고자 계산된 주시를 단순구조 및 복잡한 구조의 구조보정에 적용하고 그 심도단면도를 초동주시를 이용해 구한 심도단면도와 비교하였다.

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