• 지구물리와물리탐사
• /
• 제24권4호
• /
• pp.139-148
• /
• 2021

많은 경우 시간영역 유도분극 탐사에서 음의 겉보기 충전성은 잡음으로 간주되어 해석에 사용되지 않는다. 그러나 음의 겉보기 충전성은 실제 지하의 충전성 분포에 의해 발생할 수 있으며, 특정 모델의 경우 음의 겉보기 충전성이 획득된 자료의 대부분을 차지하기도 한다. 음의 겉보기 충전성은 송수신 전극과 충전성 이상대의 기하학적 위치에 의해 결정되며, 특히 쌍극자 배열을 사용할 경우 출현 빈도가 높다. 이 연구에서는 음의 겉보기 충전성 자료를 포함하는 시간영역 유도분극 탐사 2차원 모델링 및 비선형 역산에 관한 수치 실험을 통하여 음의 겉보기 충전성 자료가 역산결과에 미치는 영향을 분석하였다. 분석 결과, 음의 겉보기 충전성 자료도 지하의 충전성 분포에 대한 정보를 가지고 있으므로 역산에 포함되어야 하는 것으로 확인되었다.

• 한국전산구조공학회논문집
• /
• 제26권4호
• /
• pp.223-231
• /
• 2013

본 논문에서는 물성이 균일하지 않은 반무한 고체영역의 탄성파속도 분포를 재구성하기 위한 시간영역 Gauss-Newton 전체파형 역해석 기법을 소개한다. 반무한 영역을 유한 계산영역으로 치환하기 위하여 유한영역의 경계에 수치적 파동흡수 경계조건인 perfectly-matched-layers(PMLs)를 도입하였다. 이 역해석 문제는 PML을 경계로 하는 영역에서의 탄성파동방정식을 구속조건으로 하는 최적화 문제로 성립되며, 표면에서 측정된 변위응답과 혼합유한요소법에 의해 계산된 응답간의 차이를 최소화함으로써 미지의 탄성파속도 분포를 결정한다. 이 과정에서 Gauss-Newton-Krylov 최적화 알고리즘과 정규화기법을 사용하여 탄성파속도의 분포를 반복적으로 업데이트하였다. 1차원 수치예제들을 통해 Gauss-Newton 역해석으로 부터 재구성된 탄성파속도의 분포가 목표값에 충분히 근사함을 보였으며, Fletcher Reeves 최적화 알고리즘을 사용한 기존의 역해석 결과에 비해 수렴율이 현저히 개선되고 계산 소요시간이 단축됨을 확인할 수 있었다.

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

파형 역산에 사용하는 로그 목적함수는 관측 자료와 모델링 자료의 로그값의 차이를 최소화하는 목적함수이다. 라플라스 영역 파형 역산에서는 주로 로그 목적함수와 유사 헤시안의 대각 성분을 이용하여 최적화를 수행한다. 이 때 유사 헤시안의 대각 성분이 0 또는 0에 가까운 값이 되는 것을 막기 위해 레벤버그-마쿼트 알고리듬을 적용한다. 본 연구에서는 로그 목적함수의 유사 헤시안의 대각 성분을 분석하여 음향파 라플라스 영역 파형 역산에서는 유사 헤시안의 대각 성분이 0 또는 0에 가까운 값을 가지지 않음을 보였다. 따라서 로그 목적함수의 유사 헤시안을 이용한 경사 방향 정규화시 레벤버그-마쿼트 알고리듬을 적용할 필요가 없다. 수치 예제에서 인공합성 자료와 현장 자료를 이용해 레벤버그-마쿼트 기법 없이도 역산 결과를 얻을 수 있음을 보였다.

• 지구물리와물리탐사
• /
• 제24권4호
• /
• pp.171-179
• /
• 2021

시간영역 유도분극 탐사 자료로부터 Cole-Cole 변수를 추정하는 2차원 역산법을 개발하였다. 모든 유도분극 과도 전위 자료를 역산하여 전기비저항, 충전성, 완화 시간 및 주파수 승수 등의 2D Cole-Cole 변수를 추정하였다. 개발된 역산법은 2단계로 구성된다. 우선 음의 겉보기 충전성 문제를 피하기 위하여 측정된 유도분극 반응을 전류 주입 중 겉보기 전기비저항으로 변환하였다. 1단계 역산에서는 시간에 따라 항상 증가하는 전기비저항을 추정하는 4차원 역산을 통하여 각 역산 블록에서의 전기비저항 시계열 모델을 구축하였다. 2단계 역산에서는 4차원 역산에서 얻어진 전기비저항 시계열 자료를 역산하여 Cole-Cole 변수를 추정하였다. 이때 격자 탐색법을 통하여 참값에 근접한 초기 모델을 설정하는 방법을 통하여 신속한 역산이 가능하였다. 마지막으로 수치 자료에 대한 역산 실험을 통해 개발된 알고리즘이 Cole-Cole 지하 모델을 효과적으로 영상화할 수 있음을 확인하였다.

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

The traditional and still most widely used, test methods for concrete structures are destructive method, such as coring, drilling or otherwise removing part of the structure to permit visual inspection of the interior. While these methods are highly reliable, they are also time consuming and expensive, and the defects they leave behind often become focal point for deterioration. In this study, tomography by theoretical inversion method in case of elastic wave using impact-echo method among concrete non-destruction test method was made. Taken model experiments are theoretical inversion method and time domain and frequency domain test on pier test model at laboratory level. Also experiment concerning frequency domain on 3 kinds of tunnel model with I-dimension form was carried out.

• 지구물리와물리탐사
• /
• 제8권2호
• /
• pp.163-169
• /
• 2005

속도중합을 역산을 이용하면 탄성파 자료처리에서 있어서 다양한 처리가 가능하므로, 이 분야는 최근에 들어 매우 유용한 영역으로 주목을 받고 있다. 이러한 다양한 처리에 속도중합 역산을 응용하기 위해서는 사용하는 역산이 잡음에 강하면서도 고해상도의 속도중합 결과를 얻을 수 있어야 한다. 이러한 성질을 갖는 역산 방법들 중에서 가장 성공적인 방법 중의 하나라고 볼 수 있는 반복적 가중의 최소자승법(Iteratively Reweighted Least-Squares: IRLS)의 이론적 배경과 구현 방법을 소개하고, 기존 기술 특성과 한계성을 살펴보았다.

• 지구물리와물리탐사
• /
• 제14권3호
• /
• pp.220-226
• /
• 2011

탄성파동방정식에서는 변위가 수직 및 수평방향으로 정의된다. 실제 탐사에서는 수직변위와 수평변위를 모두 측정할 수 있기 때문에 이를 이용하여 방향성을 갖는 변위벡터를 구성할 수 있다. 본 연구에서는 이러한 변위벡터의 크기를 목적함수로 이용하는 주파수 영역 탄성파 파형역산 기법을 제안하고자 한다. 변위벡터 목적함수는 주파수 영역 파형역산 알고리듬에 적용할 경우 기존의 역전파 알고리듬과 동일한 방식으로 역산을 수행할 수 있다. 변위벡터 목적함수를 이용하여 Marmousi 모델과 SEG/EAGE 암염 모델의 합성탄성파 자료를 역산한 결과, 기존의 역산기법에 비해 RMS 오차가 안정적으로 감소하였다. 특히, Marmousi 모델의 밀도와 SEG/EAGE 암염 모델의 암염 하부의 저속도층을 실제 모델에 더 가깝게 구현할 수 있었다. 변위벡터의 크기를 목적함수로 사용할 경우 경사방향이 수치적으로 불안정한 형태로 정의되므로 이를 안정화시키기 위한 추가적인 연구가 필요할 것이다. 또한 본 논문에서 제안한 변위벡터 목적함수를 이용한 파형역산을 수행하기 위해서는 다성분 탐사자료 획득이 필수적이므로 육상탐사에서의 다성분 탐사나 해저면 다성분탐사(OBC, Ocean Bottom Cable) 등의 연구와 병행되어야 할 것이다.

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

주파수 영역 항공 전자탐사 자료의 holistic 역산은 사전정보가 충분한 경우에 측정 자료의 보정, 처리 및 역산을 동시에 진행하는 역산법으로 개발되었다. 이 연구에서는 사전정보가 없는 경우에도 적응 가능하도록 holistic 역산을 발전시킨다. 이 역산에서는 수직방향으로 전기적 물성이 부드럽게 변화하는 다층 구조를 가정하고, 시스템 오차로 인한 편향 매개변수를 단순화하며. 고도에 따른 수평 방향의 평활화 조건을 적용하고, 클러스터를 이용하여 병렬 계산을 수행한다. 사전정보를 전혀 이용하지 않고 holistic 역산법으로 800만 개의 자료를 340만 개의 매개댄수에 대해 역산한 결과, 이와는 별개로 계산된 보정 계수, 다운홀 로깅 자료 및 지하수위 자료와 잘 일치하는 결과를 굴을 수 있었다. 이로부터 holistic 역산의 성공 여부가 정교하게 작성된 초기 모형이나 탐사 지역의 특별한 사전정보와는 무관함을 알 수 있었다. 또한. holistic 역산은 높은 고도에서의 원점 측정을 필요로 하지 않으므로 자료 획득에 필요한 비용을 상당히 절감할 수 있을 것이다.

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

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

Frequency-space domain full-wave tomography is a promising technique for delineating detailed subsurface structure with high resolution. However, this method requires criteria for the selection of a set of optimal temporal frequency components, to achieve stability in the sequence of inversion processes together with computational efficiency. We propose a method of selecting optimal temporal frequencies, based on wavenumber continuity. The proposed method is tested numerically and is shown to be able to select an optimal set of frequency components that are sufficient to image the anomalies.