• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2011년도 학술대회 및 정기총회
• /
• pp.41-44
• /
• 2011

• 자원환경지질
• /
• 제28권5호
• /
• pp.481-485
• /
• 1995

본 논문에서는 탄성파 굴절법 탐사자료를 이용하여 천부지층의 속도와 심도를 결정하기 위한 새로운 접근방법을 소개한다. 굴절법 자료로부터 초동을 발췌한 후 실제 합성단면도를 이러한 초동의 시간이동에 해당하는 단위 델타 함수로 대치할 수 있다고 가정하였다. 주시의 계산은 발사법 파선추척을 이용하였다. 감쇠 최소자승법의 적용을 위한 편미분치의 계산은 이론주시의 계산과 동시에 해석적으로 구하였다. 본 역산법은 합성자료와 현장자료에 적용하여 성공적인 결과를 가져왔으며, 초기 가정 모델이 실제 모델과 많이 다르더라도 저주파수 대역에서 매우 양호한 결과를 보여주는 장점을 지닌다.

• 한국방재학회 논문집
• /
• 제1권1호
• /
• pp.139-155
• /
• 2001

본 연구는 수신함수 역산 기법을 이용하여 새로 설치된 한국기상청(KMA)의 광대역 지진 관측망 하부의 수직 속도모델을 찾는 것이다. 본 연구에서는 Moho(지각-맨틀 경계) 불연속면과 퇴적층-기반암의 접촉면과 같은 임피던스 차가 큰 두 지층사이에서 발생되는 P파에서 S파로 전환파인 Ps 및 다중 반사파 위상들을 이용한다. 원격지진자료는 서울(SEO), 인천(INCN), 대전(TEJ), 서산(SOS/SES), 강릉(KAN), 울진(ULC/ULJ), 대구(TAG), 부산(PUS), 그리고 울릉도(ULL) 관측소의 원격지진 P파 파형자료를 이용하였다. 광주와 춘천 관측소에서는 Moho 전환파인 Ps 도착시간과 Radial 성분 수신함수 파형이 후방위각에 따라 일치하지 않고 파형도 명확하지 않음이 발견되었다. 수신함수 역산결과 지각두께는 경기육괴에 속해있는 인천, 서울, 서산관측소에서는 29 km, 강릉(KAN) 관측소에서는 28 km, 옥천습곡대에 속해있는 대전(TEJ) 관측소에서는 32 km, 경상분지에 속해있는 대구(TAG) 관측소에서 34 km, 부산(PUS) 관측소에서는 33 km, 영동-광주 함몰지대인 광주(KWJ) 관측소에서 32 km, 영남육괴의 동쪽 경계에 위치한 울진관측소에서는 28 km, 그리고 동해의 울릉도의 울릉도관측소에서는 17 km로 각각 나타났다. 인천, 서산, 광주 그리고 강릉 관측소의 Moho 불연속면의 속도구조 양상은 약 $3{\sim}5km$ 두께의 완만한 속도 전이대를 갖는 것으로 나타났다. 강릉, 울진, 부산 관측소의 상부지각(${\sim}5km$)은 고속도의 복잡한 지각구조를 보이고 있다. 경상분지에 속해있는 대구(TAG)와 부산(PUS) 관측소에서는 한반도 서부 지역 (INCN, SEO, SOS, TEJ, KWA 관측소)의 얇은($29{\sim}32km$) 지각 두께에 비해 매우 두꺼운 지각 두께를 갖는 것으로 나타났다. 울릉도(ULL) 관측소하부의 지각 두께는 17 km인 준해양성(suboceanic) 지각으로 상부지각에서 고속도를 나타내는 복잡한 속도 구조를 보이며, 서쪽 방향에서 들어오는 Ps파형의 진폭은 다른 방향에서 입사하는 파형의 진폭에 비해 상대적으로 큰 것으로 확인되었다.

• 지구물리와물리탐사
• /
• 제19권3호
• /
• pp.136-144
• /
• 2016

탄성파 자료의 역산은 파동방정식에 기초하고 있으므로 파동방정식의 해를 정확하게 구하는 것이 가장 중요하다. 특히, 전파형역산은 파동장 전체를 이용하기 때문에 정문제에 해당하는 모델링이 정확하게 이루어져야 신뢰할 수 있는 결과를 얻게 된다. 파동방정식의 수치해를 구하는 대표적인 기법인 유한차분법과 유한요소법은 해의 수렴성을 보장할 수 있어야 하는데, 해의 수렴성은 이론적으로 일반화된 증명이 되어 있으나 실제 문제에 적용할 경우 일관성과 안정성을 분석해야 한다. 모델링 결과의 일관성은 송신원 함수의 구현이 매우 중요한 부분인데, 유한차분법은 디랙 델타 함수(Dirac delta function)를 나타낼 때 격자 간격으로 표준화된 싱크 함수(sinc function)를 사용해야 하는 반면 유한요소법은 격자 간격에 관계없이 기저함수 값을 사용하면 된다. 주파수 영역 파동방정식을 사용할 경우 송신 파형 함수의 스펙트럼을 정확하게 표현하기 위해 샘플링 이론으로 정의되는 시간 간격보다 더 조밀한 샘플링 간격을 사용하고 나이퀴스트(Nyquist) 주파수보다 더 높은 주파수를 최대 주파수로 사용해야 한다. 또한, 복소 각주파수를 사용하는 경우 감쇠 파동방정식을 만족하기 위해서는 송신 파형 함수를 먼저 감쇠한 후 사용해야 한다. 이러한 요건들이 모두 만족되었을 때 신뢰할 수 있는 역산 알고리즘 개발이 가능하다.

• 지구물리와물리탐사
• /
• 제13권3호
• /
• pp.235-242
• /
• 2010

2003년 3월 23일 한반도 홍도 해역, 3월 30일 한반도 백령도 해역에서 규모 5.0 및 4.8의 중규모 지진이 잇달아 발생하였다. 파형역산법을 적용하여 두 지진의 발진기구를 분석한 결과 3월 23일 지진은 약간의 역단층 성분을 포함하는 주향이동단층의 특성을 보이며, 3월 30일 지진은 정단층 특성을 보인다. P축의 방향은 동북동-서남서 방향으로 한반도 주응력 방향에 대한 이전 연구결과와 부합한다. 파형역산에 의한 각 지진의 모멘트규모는 4.7과 4.5로 결정되었으며, 스펙트럼 분석에 의한 각 지진의 모멘트규모는 4.8과 4.6으로 결정되었다. 홍도 해역 지진은 국지적인 응력보다는 황해에 전체적으로 작용하는 판구조적 응력 분포를 반영한다. 그러나 백령도 해역 지진은 기존의 섭입대나 충돌대에 의해 발생한 연약대가 인장방향의 응력을 받아 재활되어 발생하는 것으로 추정되며, 충돌대의 한반도 연장 가능성에 대한 지진학적 관측 증거를 간접적으로 제시한다.

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

Classical full waveform inversion for velocity estimation defines the objective function as the $l^2$ -norm of differences between the modeled and the observed wavefields. Although widely used, the results of this method have been less than satisfactory. A moderate improvement of this method is to define the objective function as the $l^2$ -norm of differences between the logarithms of the modeled and observed wavefields. In this paper we propose new objective functions of waveform inversion. They produce better results in sub-salt imaging than those of the classical and the logarithmic objective functions. One objective function defines the residual as the difference between $L^{th}$ power of the modeled wavefields and that of the observed wavefields. Another defines the residual as the difference between the integral of the $L^{th}$ power of the modeled wavefields and that of the observed wavefields. We apply these new objective functions to the synthetic SEG/EAGE salt model, and show that our new waveform inversion algorithms provide more accurate results than those of the classical and logarithmic waveform inversion methods.

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

For scaling of the gradient of misfit function, we develop a new pseudo-Hessian matrix constructed by combining amplitude field and pseudo-Hessian matrix. Since pseudo- Hessian matrix neglects the calculation of the zero-lag auto-correlation of impulse responses in the approximate Hessian matrix, the pseudo-Hessian matrix has a limitation to scale the gradient of misfit function compared to the approximate Hessian matrix. To validate the new pseudo- Hessian matrix, we perform frequency-domain elastic full waveform inversion using this Hessian matrix. By synthetic experiments, we show that the new pseudo-Hessian matrix can give better convergence to the true model than the old one does. Furthermore, since the amplitude fields are intrinsically obtained in forward modeling procedure, we do not have to pay any extra cost to compute the new pseudo-Hessian. We think that the new pseudo-Hessian matrix can be used as an alternative of the approximate Hessian matrix of the Gauss-Newton method.

• 지구물리와물리탐사
• /
• 제22권4호
• /
• pp.225-238
• /
• 2019

최적의 심도 영역 속도를 도출하기 위한 구조보정 속도분석(MVA, migration velocity analysis) 기법을 해양에서 취득한 원거리 다중채널 탄성파 자료에 적용하여 그 효용성을 확인한다. 지금까지 통상적으로 수행된 시간 영역 자료처리 결과는 지질학적 층서해석에는 무리 없는 결과이나, 어느 정도 가능성이 있는 플레이나 리드 지역에서의 유가스 탐사에서는 저류층 지질모델 구축, 시추 설계, 매장량 계산에서 반드시 심도 영역 속도 구조와 영상이 필요하다. 데이터 영역에서 근사 방식을 사용한 공통 중간점 기반 속도 분석으로부터 도출한 속도는 처음부터 오차를 내재하여 불확실성이 높다. 반면에, 이를 보완해 줄 검층 자료가 없는 상황에서 실측 규모의 속도 구조를 도출하는데 있어 이미지 영역 구조보정 속도분석 기법은 상당히 효율적인 방법이다. 이 연구에서는 해양에서 취득한 다중 채널 탄성파자료에 대해 합리적인 결과를 도출하기 위해 시간 영역에서 신호의 품질을 최적화하고, 이 자료에 대하여 반복적으로 MVA 기법을 적용함으로써 심도영역 속도 및 구조보정 단면도를 도출하였다. 시간 영역 속도를 단순히 Dix 방정식에 의해 심도영역으로 변환한 속도를 이용하여 생성한 결과(공통 수신점 모음도 및 중합 단면도)와 MVA 기법을 이용한 심도영역 자료처리를 통해 도출된 속도를 이용하여 생성한 결과를 비교함으로써, 심도영역 결과가 보다 합리적임을 확인하였다. 심도 영역으로 도출된 속도는 중합전 심도 구조보정에 바로 사용될 수 있을 뿐만 아니라, 현장 자료의 파형역산 적용 시 초기 모델로 활용함으로써 역산 수행 과정에서 발생할 수 있는 국부 최소(local minima) 문제를 최소화할 수 있다.

• 지구물리와물리탐사
• /
• 제24권1호
• /
• pp.6-19
• /
• 2021

탄성파 탐사는 속도와 밀도 같은 지하 매질 물성 정보를 파악하고 지하 지층 구조를 영상화 할 수 있으며, 이를 위한 다양한 알고리즘 개발이 이루어지고 있다. 이러한 알고리즘의 성능 검증을 위해 다양한 기준 모델이 사용되며, 정확도의 경우 참 물성 데이터와의 평균 제곱근 오차(Root Mean Squre Error, RMSE)를 통해 정량적으로 평가할 수 있다. RMSE는 수치적으로 단순하다는 장점이 있지만 구조적인 품질과의 상관도가 높지 않다는 한계가 있다. 이러한 한계를 보완하기 위해 인간지각시스템을 반영한 FR-IQA (Full Reference Image Quality Assessment) 기법이 연구되고 있으며, 지하 물성 정보 데이터를 다룰 수 있는 FR-IQA 기법들을 선정하였다. 본 연구는 물성 정보 도출 알고리즘으로 완전 파형 역산을 선정하여 세 가지 기준 모델에서 수치예제 실험을 진행하였으며, 선정 된 FR-IQA 기법들을 이용하여 물성 정보 도출 알고리즘 정확성 평가를 수행하였다. 주요 구조 정확성 평가 시 암염모델 하부 구조의 경우 구조적으로 좋지 않음을 육안으로 확인할 수 있었으나 RMSE 값은 감소하며 결과의 부정확성을 표출하지 못하였다. 반면, 몇몇 FR-IQA의 경우 결과의 부정확성을 수치적으로 표출하는 것을 확인하였다.

• 한국음향학회지
• /
• 제36권5호
• /
• pp.293-299
• /
• 2017

수중에서 기동하는 소형무인 플랫폼은 미래의 핵심적인 해군전력이지만, 이와 동시에 아군세력을 위협하는 비대칭 세력이기도 하다. 이러한 소형무인 플랫폼을 탐지, 식별하기 위한 하나의 방법으로 반복적 수치해석 기법에 기반한 물성추정 기법을 소개한다. 물성추정 기법은 표적의 위치정보뿐만 아니라 표적의 물리적인 물성치를 추정할 수 있기 때문에 정지한 표적을 탐지 및 식별할 수 있으며, 다수의 수중표적도 동시에 탐지 가능하며, 표적과 주변 자연물체와의 구분이 가능한 장점을 가진다. 본 연구에서는 기존에 제안한 기법을 보다 확장하여, 실제 호수시험 획득자료를 이용하여 실제 수중 소형표적을 대상으로 물성추정을 수행하였다. 그 결과 단 한 개의 채널신호만으로도 소형플랫폼에 대한 물성추정이 가능하였으며, 현장자료에도 강인한 결과를 획득할 수 있었다.