• 지구물리와물리탐사
• /
• 제8권3호
• /
• pp.207-217
• /
• 2005

이 연구에서는 배경 전기전도도 구조가 1차원 층서구조가 아닌 경우에도 적분방정식법을 이용하여 비교적 정확히 모델링 할 수 있는 근사적 방법에 대한 유도 과정을 제시하였으며, 이를 기반으로 전자탐사 수치 모델링을 구현하여 그 정확도 및 효율성에 대해 고찰하였다. 이 연구에서 제시된 방법은 통상적인 적분방정식법과 마찬가지로 1차원 층서구조에서의 Green함수를 그대로 이용할 수 있어 계산 측면에서 매우 효율적이며, 대체적으로 적분방정식을 구성하기 위한 행렬의 크기를 줄일 수 있다. 본 연구에서 개발된 근사 적분방정식 알고리듬을 간단한 수치 모델에 대하여 적용한 결과 정확도는 반응을 구하고자 하는 이상체 부근에서는 거의 차이가 없었으며, 또한 계산 시간에 대한 분석을 통해 이 모델링 알고리듬을 이용하여 역산을 수행하고자 할 때 계산 시간을 대폭 감소시킬 수 있음을 확인하였다. 본 연구에서 제안된 방법은 1차원 층서구조가 아닌 이미 알려진 공통된 구조론 갖는 여러 지질 모델에 대한 전자탐사 반응을 매우 효율적으로 구할 수 있을 것으로 판단된다.

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

쌍극자 전기탐사현장에서 드물지않은 음의 겉보기비저항에 대해 선진 장비회사의 메뉴얼들에서 언급되되 정확한 현상이 규명 안되었으나, 국내에서 계측 측면 및 가능한 지질모델 제시와 함께 전산모델링 결과와 $-{\rho}a$ 포함유무 3차원 역산결과가 일부 보고된 바 있다. 전도성 탄층, 전봇대 접지, 신호등, 철케이싱 시추공 등과 관련한 현장문제가 계속되어 해석시 딜레마 극복 지침이 필요하다. 여기서는 단순 철케이싱 한개의 효과 전산 모델링 결과 및 일부 측선 설계요령과 해석시 유의사항을 고찰한다. $-{\rho}a$는 분명한 지하정보이지만 잘못 해석시 진실과 먼 결과 오류 가능성이 높으며, 향후 역산 알고리즘의 유용한 발전이 절실하다.

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

한반도 남부 지각 속도구조를 밝히기 위해서, 서북서-동남동 방향의 2002년도 측선 294 km와 북북서-남남동 방향의 2004년도 측선 335 km를 따라 인공적으로 발생시킨 지진파 자료를 각각 120초와 150초 기록하였다. 초동주시 역산과정을 통하여 속도단면을 작성하였으며, 역산 시 초기모델은 측선 주변의 고정관측소에서 기록한 원거리 지진자료의 수신함수역산으로 구한 1차원 속도구조자료를 활용하였다. 파선경로는 2.0 km와 7.1 km 깊이에 속도 6.0 km/s와 7.1 km/s를 갖는 굴절면이 존재하며, 굴절파 속도 $7.8{\sim}8.1\;km/s$의 모호면은 $30.8{\sim}36.1\;km$ 깊이에 존재함을 보인다. 속도단면은 옥천계 하부 $6{\sim}7\;km$ 깊이에 상당한 규모의 저속도층이 15 km 깊이의 속도 불연속면 상부에 존재하며, 영동단층은 10 km 이상 깊이까지 연장되어 있고, 최대 4.2 km 정도두께를 갖는 경상분지 하부에 고속도층이 얕게 분포하는 것으로 분석된다.

• 지구물리
• /
• 제8권1호
• /
• pp.45-48
• /
• 2005

한반도 남부 지각 속도구조를 밝히기 위해서, 서북서-동남동 방향의 2002년도 측선 294 km와 북북서-남남동 방향의 2004년도 측선 335 km를 따라 인공적으로 발생시킨 지진파 자료를 각각 120초와 150초 기록하였다. 초동주시 역산과정을 통하여 속도단면을 작성하였으며, 역산 시 초기모델은 측선 주변의 고정관측소에서 기록한 원거리 지진자료의 수신함수역산으로 구한 1차원 속도구조자료를 활용하였다. 파선경로는 2.0 km와 7.1 km 깊이에 속도 6.0 km/s와 7.1 km/s를 갖는 굴절면이 존재하며, 굴절파 속도 7.8∼8.1 km/s의 모호면은 30.8∼36.1 km 깊이에 존재함을 보인다. 속도단면은 옥천계 하부 6∼7 km 깊이에 상당한 규모의 저속도층이 15 km 깊이의 속도 불연속면 상부에 존재하며, 영동단층은 10 km 이상 깊이까지 연장되어 있고, 최대 4.2 km 정도두께를 갖는 경상분지 하부에 고속도층이 얕게 분포하는 것으로 분석된다.

• 지구물리와물리탐사
• /
• 제20권3호
• /
• pp.121-128
• /
• 2017

국내에서 시간영역 전자탐사(time-domain/transient electromagnetic, TEM) 자료의 해석은 1차원 알고리듬에 주로 의존하고 있는 실정이기 때문에 정밀한 해석을 위하여 3차원 모델링 및 역산 해석 프로그램의 개발이 필요한 상황이다. 이 연구에서는 엇갈린 격자를 이용한 시간영역 유한차분(staggered-grid finite-difference time-domain, FDTD)법에 기초하여 3차원 TEM 반응 모델링 알고리듬을 개발하였다. 시간영역 전자탐사의 모델링을 위해 맥스웰 방정식을 현시적 중앙점 FDTD법을 이용하여 이산화하였으며 수치 안정성을 높이기 위해 가상 변위전류항을 도입하였다. 일반적으로 많이 활용되는 소형 코일 송신원을 수치적으로 구현하여 균질 반무한 공간에서의 해석해와 비교 검증하고 3차원 이상체에 대한 반응을 분석하였다. 이 연구에서 개발된 모델링 프로그램은 향후 TEM 전자탐사 자료의 정밀 해석에 기초가 될 것으로 기대한다.

• 지구물리와물리탐사
• /
• 제3권1호
• /
• pp.7-12
• /
• 2000

본 논문은 3차원 탄성파자료로부터 저류층의 공극율 분포를 파악하기 위한 자료 전처리 및 역산법에 대한 연구이다. 전처리과정은 음파 및 밀도 검층기록에서 도출한 반사계수 시계열을 사용하여 탄성파기록의 진폭을 보정하는 공정이다. 연구지역에서 획득한 8개의 검층기록과 그 위치의 탄성파 기록으로부터 산출된 보정함수를 크리깅(kriging)하여 모든 탄성파기록의 보정함수를 산출하였다. 일정 대역폭을 보유한 탄성파기록으로부터 검층과 같은 광대역 임피던스 시계열을 도출하기 위해서는 한정된 대역폭의 탄성파 기록에서 손실된 저주파와 고주파 성분을 복구하여야 한다. 본 연구에서는 autoregressive(AR) 방법을 사용하였다.

• 지구물리와물리탐사
• /
• 제11권4호
• /
• pp.302-309
• /
• 2008

3차원 전기비저항 토모그래피는 토목 및 수리지질분야 등에서 매우 고분해능의 지하구조 영상을 획득할 수 있는 강력한 성능을 제공한다. 이 논문에서는 특히 3차원 이상대가 주탐사 대상인 경우에 대한 3차원 전기비저항 토모그래피의 두가지 성공적 응용사례를 소개한다. 첫번째 예는 수리지질분야에의 응용을 위한 파쇄대 영상화 사례로서, 여기서는 시추공의 공곡변화가 중요한 문제로 대두되었으며 이에 따라 시추공 공곡변화를 토모그래피 역산과정에 포함함으로써 파쇄대의 3차원 영상을 획득할 수 있었다. 두번째 예에서는 크기 $2m{\times}2m$의 내부가 비어있는 터널의 영상화 문제를 다루었으며, 3차원 토모그래피를 통하여 터널의 발달을 3차원 공간상에 영상화할 수 있었다. 이들 예로부터 3차원 전기비저항 토모그래피가 3차원 지하구조 영상화에 강력한 성능을 보임과 아울러 특히 3차원 이상대인 파쇄대, 공동 또는 터널의 영상화에 더욱 향상된 영상화 능력을 제공할 수 있음을 보였다.

• 지구물리와물리탐사
• /
• 제12권4호
• /
• pp.361-366
• /
• 2009

전자탐사에서 수평 층서구조에 대한 Green 함수의 계산은 전자기 반응의 모델링에서 핵심적인 부분을 담당한다. 해석적으로 구해진 핵함수의 Hankel 변환으로 계산되는 Green 함수는 핵함수의 대수적 등가 표현방식에 의해 그 정확도가 결정된다. 특히 3차원 모델링의 경우 Green 함수 계산 횟수가 매우 많아서 Hankel 변환 계산이 전체 계산시간의 상당량을 차지하므로, 빠르고 정확한 Hankel 변환의 계산을 위해서 선형 수치필터를 이용한다. 최근 많이 시도되는 3차원 역산을 위한 모델링에서는 송신점에서의 특이성 문제를 피하기 위해 전기장을 1차장과 2차장으로 나누어 계산하는 것이 보통이다. 이 연구에서는 균질 반무한공간에 대해 지표면에 놓인 다섯 종류의 송신원에 대한 지하 매질에서의 전기장 세성분을 Hankel 변환을 이용하여 정리하고, 그 계산 방법에 대해 고찰하였다. 그리고 2중 반공간에서 EM1D를 이용하여 공기와 바다의 영향을 모두 고려한 전자기장을 계산할 때, 보다 정확한 해의 계산을 위해 TE 및 TM 모드에서의 반사계수를 유도하였다. 여기서 정리한 해를 이용하면 MT 문제는 물론, 해양 전자탐사의 경우에도 1차장을 정확히 계산할 수 있으므로 3차원 역산에서 보다 정확하고 효율적인 감도 계산이 가능할 것이다.

• 지구물리와물리탐사
• /
• 제9권4호
• /
• pp.291-298
• /
• 2006

전기비저항 탐사법은 저수지 누수구간 탐지에 매우 효과적이고 실질적인 물리탐사법이다. 일반적으로 저수지 제체의 마루에서 수행되는 전기비저항 탐사는 저수지가 2차원 구조를 갖고 있다고 가정한다. 하지만 저수지 제체에서 수행되는 전기비저항 탐사 자료는 저수지의 3차원 지형에 의해 크게 영향을 받는다. 본 연구에서는 유한요소법을 사용하는 3차원 전기비저항탐사 모델링 프로그램을 통하여 저수지의 3차원 구조가 전기비저항 탐사자료에 미치는 영향을 평가하였다. 제체의 상류사면, 상류쪽 마루, 마루중앙, 하류쪽 마루와 하류사면에 위치한 측선에서 얻어진 전기비저항 탐사자료를 비교하였다. 3차원 효과는 측선의 위치에 따라 서로 다른 양상을 보였으며, 측선이 하류사면 쪽에 위치할 경우 그 영향이 더 크게 나타났다. 또한 지형효과보다는 저수지를 구성하는 요소들의 전기비저항 값이 전기비저항 탐사자료에 더 큰 영향을 미치는 것으로 확인되었다. 마지막으로 누수구간이 존재하는 제체에 대한 3차원 모델링 자료에 대하여 2.5차원 역산을 수행하여 2차원 전기비저항 탐사법을 이용한 누수구간 탐지의 가능성을 분석하였다.

• 지구물리와물리탐사
• /
• 제18권1호
• /
• pp.1-8
• /
• 2015

뉴질랜드 Ngauruhoe 화산에 설치된 OTVZ 지진관측소 하부의 S파속도(${\nu}_s$)를 규명하기 위해, 2011년 11월 11일에서 2013년 9월 11일까지 기록된 원거리 지진자료($Mw{\geq}5.5$) 중에서 신호 대 잡음비가 높은 127개의 자료만을 이용하여 수신함수를 계산하였다. 수신함수는 시간영역에서 반복적 곱풀기 방법으로 계산되었으며, 21개의 수신함수에 유전자 알고리즘을 적용하여 역산을 실시하였다. 역산으로부터 구한 관측소 하부의 1차원 속도모델은 ${\nu}_s$가 3.7 km/s에서 4.7 km/s로 급격히 변하는 모호면이 14 km 깊이에 존재하며, 지각의 평균 ${\nu}_s$는 3.4 km/s임을 보인다. 하부지각에는 평균 ${\bar}{\nu}_s$ 가 3.1 km/s인 저속도층이 두께 9 km 이내로 존재하며, 관측소에서 멀어질수록 두께가 얇아지는 것으로 분석된다. 또한, 상부맨틀에도 ${\nu}_s$가 4.3 km/s 이하인 저속도층이 10 km 이상의 두께로 존재하며, 이러한 하부지각과 상부맨틀 내에 존재하는 저속도층과 상대적으로 얇은 지각은 태평양판의 섭입에 따른 마그마 활동과 관련이 있을 것으로 추정된다.