• 지구물리와물리탐사
• /
• 제10권4호
• /
• pp.268-274
• /
• 2007

2차원 측선을 따라 획득한 MT 탐사 자료를 3차원 역산하여 제주도 지전기 구조의 해석을 시도하였다. 이에 앞서 TM 모드 자료와 TE 모드 자료를 각각, 또는 두 모드의 자료를 동시에 이용하여 2차원 역산을 수행하였다. 이 세가지 2차원 역산을 통해 얻은 영상을 3차원 영상과 비교분석하기 위해 측선을 포함하는 3차원 모델을 구성하여 3차원 역산 결과를 도출하였다. 모든 역산 결과에서 전체적으로 동일한 층서구조를 확인할 수 있었으나 3차원 역산에서는 보이지 않는 전기비저항 이상대가 2차원 역산에는 나타난다. 이는 측선 직하부에 있지 않은 3차원 이상체의 영향으로 2차원 역산 결과가 왜곡된 것이라 유추할 수 있다. 역산 영상의 비교 분석 결과 측선 하부에서 제주도의 지전기 구조는 심도 5 km 이내에서는 전체적으로 고비저항.저비저항.고비저항의 3층 수평 층서 구조임을 알 수 있었고, 3차원 역산이 더욱 신뢰성 있는 영상을 보임을 확인하였다.

• 지구물리와물리탐사
• /
• 제6권4호
• /
• pp.199-206
• /
• 2003

MT 탐사자료의 역산에 있어서 지하의 전기비저항과 함께 정적효과를 파라미터로 설정하여 동시 역산을 수행하는 알고리듬을 하나의 지하구조 모델에 각기 다른 양의 정적효과를 포함시킨 4개의 자료에 대하여 적용시키고 이를 정적효과가 전혀 고려되지 않은 경우와 비교하여, 3차원 역산에서 정적효과가 미치는 영향 및 그 특성에 대하여 분석하였다. 일반적으로 현장자료에 정적효과가 어느 정도 포함되어 있는지에 대한 사전 정보가 전혀 없으므로 역산과정에서 이를 조절하는 trade-off 파라미터의 적절한 선택이 매우 중요하며, 본 연구에서는 모델의 smoothness와 static shift의 양을 조절하는 각각의 파라미터의 크기를 매 반복마다 구하는 알고리듬을 동시역산에 적용하였으며 4개의 이론자료에 적용한 결과 만족할 만한 결과를 얻었다. 정적효과가 포함된 자료에 대하여 정적효과를 고려하지 않은 역산(기존의 MT 역산)에서는 지표 block의 전기비저항을 바꿔 역산 스스로가 정적효과를 유발하려는 경향을 보였으며 이의 결과로 저주파수에서는 상당한 정적효과를 발생시켜 정적효과가 그리 크지 않은 경우 심부구조를 어느 정도 규명해 내는 것으로 나타났다. 그러나 고주파수에서는 이들 지표 block의 영향이 주파수에 무관하지 않게 되어 정적효과를 포함하는 자료의 겉보기 전기비저항과 위상을 동시에 만족시키지 못하게 된다. 그러나 정적효과를 파라미터로 하는 동시역산의 경우, 매우 심한 정적효과를 포함하는 자료에 대해서도 지하구조를 매우 정확히 영상화 하는 것이 가능하였다.

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

An efficient three-dimensional (3D) inversion of magnetotelluric (MT) data can be carried out by using approximate sensitivities or avoiding the calculation of a full sensitivity matrix. In this paper, we propose approximate sensitivities for efficient 3D MT inversion based on the Gauss-Newton method and test and compare four kinds of sensitivities. Applying the four sensitivities to both synthetic and field data shows that the effects of sensitivities are highly dependent on data and thus applying various combinations of sensitivities is recommended for efficient inversion and good images.

• 지구물리와물리탐사
• /
• 제10권4호
• /
• pp.259-267
• /
• 2007

정적효과가 포함된 MT 자료를 효율적으로 역산하기 위해 새로운 근사 감도(sensitivity)를 제시하였다. 근사 감도를 이용하여 감도 행렬을 구하는데 필요한 막대한 계산 시간을 합리적인 수준으로 줄이면서 정확한 감도를 이용한 경우와 비슷한 수준의 역산 결과를 얻을 수 있었다. 이 연구에서 제시한 근사 감도 외에 초기 균질 모형에 대한 감도, Broyden 법을 이용하여 업데이트한 감도, 정확한 해에 의존한 감도 등을 여러 가지 조합으로 3차원 역산을 수행할 수 있는 효율적인 방법에 대해 연구하였다. 인공합성자료 및 포항 지역에서 획득한 현장 자료에 대하여 3차원 역산한 결과, 네가지 감도의 조합으로 성공적인 역산이 가능함을 확인하였다. 이로부터 다양한 감도를 적용하여 역산을 수행하는 것이 보다 효율적이고 합리적인 역산 결과를 얻는데 도움이 됨을 알 수 있었다.

• 지구물리와물리탐사
• /
• 제17권2호
• /
• pp.95-103
• /
• 2014

MT 자료에서 3차원 전기비저항 구조와 정적효과를 동시에 구하기 위한 실용적인 역산법을 소개한다. 이 방법은 감도행렬이 필요한 Gauss-Newton법을 기반으로 하고 반복과정에서 Broyden의 방식으로 감도를 수정하는 것을 기본으로 하고 있다. 이 논문에서는 합성 MT 자료에 대한 역산실험을 통해 근사역산법의 성능과 정적효과에 대한 가중치에 대해 검토하였다. 해석적으로 구해지는 초기감도를 Broyden의 방식으로 수정하는 역산법은 초기감도를 끝까지 쓰는 역산법보다 자료오차를 줄이는데 효과적이었다. 그리고 완전한 감도행렬을 반복 중간에서 단 한번만 사용하는 근사역산법으로서는 반복 전반부에서 완전한 감도를 사용할 때 자료오차를 가장 많이 줄이는 것으로 나타났다. 정적효과에 대한 가중치는 어느 특정 한계값 이하로 선택하면 최종 자료오차에는 결정적인 영향을 주지 않는다. 합성 MT 자료에 대한 실험 결과 이 역산법은 정적효과가 포함된 MT 자료로부터 3차원 전기비저항 구조를 재현하는데 효과적임을 확인하였다.

• 지구물리와물리탐사
• /
• 제11권3호
• /
• pp.184-196
• /
• 2008

호주 Victoria주 Bendigo 북부 지역에서 심부 지질 구조를 조사하기 위하여 MT 탐사를 수행하였다. 이 탐사의 주 목적은 이 지역에 발달한 단층대의 연장 상황을 파악하는데 일차적으로 그 목적이 있다. 총 11일에 걸쳐 71개의 측점에서 MT시계열을 획득하였으며 이와 동시에 원거리기준점에서도 자료를 획득하여 자료 처리 및 해석의 신뢰도를 향상시키고자 하였다. 획득된 자료는 임피던스 추출, 원거리기준점 자기장을 이용한 로버스트 자료처리 등을 거쳐 임피던스를 추출하였으며, 3차원 역산 알고리듬을 이용하여 최종적으로 전기비저항 영상을 획득하였다. 이러한 MT 탐사 결과를 이 지역에 대한 지질 자료 및 동일한 측선에서의 탄성파 반사법 자료와 비교한 결과 상호 부합함을 알 수 있었으며 이 지역에 발달된 단층대의 연장 상황을 파악할 수 있었다.

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

A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.

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

자기지전류(MT) 자료의 3차원 역산에 대해 소개한다. MT 자료의 역산 문제는 기본적으로 악조건이므로 유일한 해가 존재하지 않는다. 이러한 비유일성을 줄이고 정확한 역산해를 구하기 위해서는 역산 시 사전정보를 추가하는 제약조건을 가해야 한다. 물리탐사 분야에서 비선형 역산에 사용되는 가장 일반적인 방법은 일련의 선형화된 역산문제를 푸는 Gauss-Newton법이다. 이 알고리듬은 수렴 시, 모델 공간에서 역산문제에 대한 목적함수를 최소화하는 최적해를 준다. 그러나 이러한 반복적 선형화기법은 3차원 MT 역산의 경우 Jacobian 행렬을 구하기 힘들기 때문에 그 유용성에 한계가 있다. 이러한 어려움은 CG법에 의해 완화할 수 있다. 선형 CG법은 Gauss-Newton 반복의 각 단계를 근사적으로 풀기 위해서 사용된다. 한편 비선형 CG법은 목적함수의 최소화에 직접적으로 적용된다. 이들 CG법은 Jacobian 행렬의 계산 및 대형 선형방정식의 해를 반복 당 세 번의 모델링으로 대치할 수 있어서 3차원 역산에 적합하다.

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

Uiseong subbasin belonging to Kyungsang basin resulted from volcanic activity in the late Cretaceous. In this study, we carry out MT and gravity survey at the Hwasan caldera, which was formed of volcanic and abyssal rocks complex, then analyze and identify geological substructure. Potential survey such as gravity and magnetic survey has been mainly carried out in former studies, so depth information for understanding substructure was not enough. To complement a potential survey, we use MT method, which has high vertical resolution. Moreover we make a simple 2D model comparing with former study. The result of MT and gravity 2D modeling shows that this area is roughly composed of 3 layers; The bottom layer is a basement. In the second layer, intrusive rocks having high resistivity is placed along the ring faults and the sedimentary layer of low resistivity is inside caldera. The highest layer is alluvium. To comprehend the 3D structure of the Hwasan caldera, we perform 3D gravity inversion, and construct the 3D model from the result of 3D gravity inversion. MT responses are calculated by using the constructed 3D model and the 3D model of the Hwasan caldera's structure is suggested after comparing the calculated values with the observed values at MT line.

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

해양과 인접한 지역에서 MT탐사 자료를 해석하는 경우, 해양은 심부 구조의 반응을 왜곡시키는 역할을 한다. 본 연구에서는 해저면 MT탐사에서 해저 지형의 영향을 제거하는 지형보정 기법을 바탕으로, 반복적으로 해양효과를 보정하는 기법을 개발하였다. 제안된 기법은 우선 관측 MT 반응에서 주변 해양의 영향을 보정하고, 그 후 해양이 없는 모델공간에서 역산을 수행한다. 주변 해양과 지하구조와의 상호결합 때문에 이 과정은 반복적으로 수행되며, 보정 결과의 변화가 미미할 때 반복과정이 종료된다. 제안된 기법의 검증을 위해 해양을 포함하는 3차원 순산 모델링을 통하여 합성 자료를 생성하였고, 1차원 및 2차원 구조에서 보정기법을 적용하였다. 대체적으로 제안된 해양효과 보정 기법은 $2{\sim}3$회의 반복단계를 거친 후 해양이 없는 경우의 지하구조를 성공적으로 복원하였다. 실제 MT 자료와 유사한 자료를 획득하기 위해 1차원 구조로 잘 알려져 있는 제주도에 대한 3차원 모델링을 수행하여 유사 현장 자료를 생성하였다. 제주도 모델의 경우, 해양효과는 약 1 Hz 이하에서 나타나기 시작하였으며 측선의 위치 때문에 해안선과 수직한 전기장 성분에서 상대 적으로 왜곡이 심하게 나타났다. 이러한 왜곡은 3회의 해양효과 보정과정을 통해 성공적으로 제거되었으며, 1차원 및 2 차원 역산은 모델링 시 가정한 제주도의 지하구조를 성공적으로 복원하였다.