• 지구물리와물리탐사
• /
• 제7권3호
• /
• pp.191-196
• /
• 2004

탄성파 역산에는 고전적인 Gauss-Newton 방법이 주로 사용된다. 이 방법은 Jacobian을 직접 계산하여 거대한 크기의 Hessian 행렬을 만드는 것을 필요로 한다. Hessian 행렬의 구성은 몇 가지의 요소들에 의해 결정되는데, 음원과 수진기의 위치, 영상화 구역(image zone), 음원 파형의 형태 등 다양한 형태의 모델링에 영향을 미치는 요소에 따라서 다른 모습으로 나타난다. 이 논문에서는 Gauss-Newton 방법에 나타나는 거대한 Hessian 행렬을 조절함으로써 Marmousi 탄성파 모델 자료를 역산하고자 한다. 또한 근사 Hessian행렬의 대안으로 두 가지의 유사 Hessian행렬들을 제시하고자 한다. 하나는 유한 폭을 갖는 Hessian행렬이고 다른 하나는 자동안정함수(automatic gain function, AGC)를 이용한 Hessian 행렬이다. 작은 크기의 모델에 대한 수치결과로부터 몇 가지의 사실을 알 수 있다. 하나는 유한 폭을 갖는 Hessian 행렬을 이용하여 얻어진 한번 근사된 속도모델은 원래의 Hessian 행렬을 이용하여 얻은 결과와 매우 유사하다는 것이고, 둘째로 자동안정함수를 이용한 근사 Hessian 행렬의 안정성이 많이 개선된다는 것이다.

• 지구물리와물리탐사
• /
• 제12권1호
• /
• pp.99-104
• /
• 2009

본 연구에서는 천부 토양층의 2차원 불균질 S파 단면을 결정하기 위해 천부 굴절법 탐사로부터 얻은 SH파 자료에 대한 파형역산 기법을 제안한다. 2차원 매질에서 SH파의 전파를 모사하기 위해 2.5차원 파동 방정식을 사용한다. 합성탄성파를 계산하기 위해 공간축으로 4차, 시간축으로 2차 근사한 staggered grid 유한차분법을 사용하여 파동 방정식을 푼다. 계산된 파형과 측정 파형의 잔여오차로 정의되는 목적함수를 하이브리드 발견적 탐색기법에 의해 최소화한다. 2차원 지하구조 모형은 각기 다른 심도 경계면을 갖는 블록과 블록 내부의 S파 속도에 의해 매개화한다. 수치실험은 암영층파 불규칙한 경계를 갖는 모델에 대해 백색잡음을 추가한 합성 SH파 자료를 이용하여 수행하였다. 지표 굴절법 자료로부터 암영층을 포함한 구조를 적절한 계산시간 내에 영상화할 수 있었다

• 한국전산구조공학회논문집
• /
• 제34권4호
• /
• pp.221-230
• /
• 2021

전역 최적화 문제의 해를 유전 알고리즘을 사용하여 얻어 완전파형역산을 수행하고 층상 반무한체의 물성치를 추정하는 기법을 제안한다. 조화 수직 하중이 작용하는 층상 반무한체의 동적 응답을 측정하고, 이를 추정 물성치를 사용하여 계산된 응답과 비교한다. 응답의 추정치는 mid-point integrated finite element와 perfectly matched discrete layer를 사용하여 구성된 thin-layer model로부터 얻는다. 전역 최적화 문제의 목적 함수는 응답의 관측치와 추정치의 차이에 대한 L2-norm으로 계산된다. 유전 알고리즘을 사용하여 전역 최적화 문제의 해를 구하여 완전파형역산을 수행한다. 제안된 기법을 기본 진동 모드 뿐만이 아니라 고차 진동 모드도 우세한 다양한 층상 반무한 매질에 적용하여, 측정치가 잡음을 포함하지 않는 경우와 포함하는 경우 모두에 대해서 제안된 완전파형역산 기법은 층상 반무한체의 재료 특성을 추정하는데 적합함을 확인할 수 있다.

• 지구물리와물리탐사
• /
• 제17권4호
• /
• pp.231-241
• /
• 2014

최근 엔지니어링 목적의 단일채널 탄성파 탐사가 많이 수행되고 있다. 단일채널 탄성파 탐사는 일반적으로 특별한 자료처리 없이 효율적으로 지하 지질구조를 파악할 수 있는 장점이 있지만, 복잡한 지질구조에 대한 정확한 영상화에는 한계를 가진다. 자원개발 목적의 다중채널 탄성파 탐사에서는 최근 파형역산에 대한 연구가 활발히 진행되고 있으며, 복잡한 지하구조에 대해서도 정확한 지하영상화 결과를 제시하고 있다. 이에 본 논문에서는 단일채널 탄성파 탐사 자료를 이용하여 지하 속도모델을 구하기 위한 탄성파 파형역산 알고리즘을 제안하였다. 제안한 알고리즘은 단일채널 탄성파 탐사를 고려하여 지하 매질을 1차원으로 가정하였으며, 벌림에 의한 지연시간을 제거하여 벌림에 의한 효과를 보정하였다. 파형역산은 안정적인 해의 계산이 가능한 가우스-뉴턴법을 이용하였다. 알고리즘은 수정된 Marmousi2 모델에 적용하여 검증하였으며, 부산항에서 얻은 현장자료에 적용해 보았다.

• 한국지반환경공학회 논문집
• /
• 제25권2호
• /
• pp.5-14
• /
• 2024

철근 콘크리트 교량 바닥판은 차량 하중과 우수침투 등으로 가장 먼저 손상이 발생하며 제설염화물 등으로 인한 철근 및 기타 금속 부재가 부식되면서 콘크리트 열화가 주로 발생한다. 교량 바닥판의 시공 상태 및 포장 내부 바닥판의 열화는 지표투과레이더(ground-penetrating radar, GPR) 탐사 자료를 이용하여 평가하고 있다. 교량의 콘크리트 열화 상태를 평가하기 위해서는 철근의 위치 및 열화 지점을 정확하게 파악하기 위한 GPR 자료 해석 기술 개발이 필요하다. GPR 탐사에서는 지반 매질의 레이더파 전파 속도 차이에 의한 반사 및 회절파 신호를 취득한다. 그러므로 이 연구에서는 GPR 탐사 자료를 이용해서 지반 매질의 레이더파 전파 속도를 추정하고 교량 바닥판의 열화를 평가하기 위한 완전파형역산(full-waveform inversion, FWI) 기술을 개발하였다. 개발된 FWI 기술의 적용성은 콘크리트 박리 및 철근 부식과 같은 교량 바닥판 열화현상을 보여주는 GPR 속도 모델을 만들어서 수치 실험을 수행하여 검증하였다. 합성 GPR 자료로 역산한 결과 교량 바닥판의 철근 위치와 열화 지점을 역산으로 계산된 속도 영상으로 확인할 수 있었다.

• 지구물리와물리탐사
• /
• 제23권1호
• /
• pp.38-49
• /
• 2020

음향파 완전파형역산은 탄성파 탐사를 통해 얻은 관측자료와 음향파 모델링자료를 맞춤으로써 지층의 속도모델을 고해상도로 구축하는 최적화 과정이다. 기존의 스트리머를 이용한 해양 탄성파 탐사 자료에 대한 음향파 완전파형역산에서는 압력자료만을 활용하여 P파 속도모델을 구축한다. 그러나 최근 다성분 해저면 탄성파 탐사기술의 발달로 다성분자료의 취득 사례가 늘고 있으며, 이에 따라 해양에서 얻어지는 다성분 자료를 활용한 음향파 완전파형역산 기법에 대한 연구가 필요하다. 이 연구에서는 수평 및 수직 입자가속도 자료를 활용한 효과적인 음향파 완전파형역산 전략을 제시하고자 한다. 이를 위해, 우선 음향파 모델링으로 제작된 압력 및 입자가속도 자료와 민감도커널을 분석하여 파형역산 과정에서 각 자료의 성분별 특성을 관찰하였다. 압력 자료에 함께 나타났던 직접파, 다이빙파 및 반사파가 수직 및 수평 입자가속도 자료에서 파동의 진행방향에 따라 분리되어 나타나는 것을 확인하였으며, 수평 입자가속도 자료는 상부의 장파장구조를, 수직 입자가속도 자료는 하부의 장파장구조와 전체 영역에서의 단파장구조를 구축하는 데 효과적임을 확인할 수 있었다. 이러한 분석 결과를 바탕으로 입자가속도 자료만을 활용해 음향파 완전파형역산을 수행하는 순차적 자료 활용전략을 제시하며, 압력자료를 얻지 못하였거나 품질이 낮은 경우에도 입자가속도 자료를 활용하는 음향파 완전파형역산을 통해 양호한 P파 속도모델을 구축할 수 있을 것으로 기대된다.

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

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

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

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

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