• 지구물리와물리탐사
• /
• 제20권4호
• /
• pp.226-231
• /
• 2017

고정밀 중력 탐사는 천부의 소규모 밀도 이상 구조를 명확하게 탐지하는 것을 목적으로 하므로, 정밀한 측정과 보정뿐 아니라 소규모 이상 구조를 명확하게 분석해 낼 수 있는 역산법도 필요하다. 이것은 역산에서의 타당치 않은(ill-posed) 문제를 해결하기 위한 안정자를 밀도가 급격하게 변화하는 경계부를 분명하게 구분할 수 있도록 설계함으로써 가능해 진다. 이 논문은 다양한 탐사 목적에 대응할 수 있도록 여러 가지 역산 알고리즘을 포함하는 2차원 고정밀 중력 역산 패키지를 구성하였다. 패키지는 Matlab 기반으로 작성하였으며, 최선의 역산 결과를 얻을 수 있도록 대화형으로 만들었다. 몇 가지 대표적인 모델에 대한 수치모델링 자료에 역산 패키지에 포함되어 있는 여러 가지 역산 알고리듬을 적용하여 역산 방법들을 비교 및 검토하였다. 마지막으로 실제 탐사 자료에도 적용하였다.

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

본 연구에서는 지형을 포함한 2차원 MT 역산 알고리듬을 개발하였다. 역산 과정시 필요한 모델 반응 계산을 위하여 유한요소법을 이용하였다. 공기와 지표면의 경계를 기준으로 고도에 따라 각 요소들의 절점을 수직으로 이동시킴으로써 추가적인 계산시간의 증가없이 간편히 지형을 구현하였다. 역산에서는 공간적인 함수로서 라그랑지 곱수를 결정하는 알고리듬을 채택하여 역산의 분해능과 안정성을 높이고자 하였다. 수치모델 실험을 통하여 TM과 TE모드 자료의 지형효과를 고찰하였고, 수치 자료의 역산을 통하여 지형을 포함한 역산의 타당성을 살펴보았다. 또한 현장 자료에 대하여 적용하여 본 연구에서 개발된 지형을 포함한 MT자료 역산 알고리듬의 적용성을 확인하였다.

• 지구물리와물리탐사
• /
• 제8권1호
• /
• pp.26-33
• /
• 2005

최근 들어 토목관련 천부층 조사에 다중 모드 표면파 위상 속도의 역산이 많은 관심을 받고 있다. 감도 분석, 그리고 합성탄성파자료와 현장자료의 역산 결과는 이 방법이 기본 모드만을 이용하는 것에 비해 매우 효과적임을 보여주고 있다. 이중 모드 레일리 파의 위상속도들에서 층의 두께와 전단파 속도에서의 조그만 변화는 고차 모드의 감도들을 (a) 다른 주파수 대역들에 모이게 하고 (b) 심도가 깊어질수록 기본 모드보다도 더 크게 한다. 이 관찰을 통해 다중 모드 위상 속도 역산을 이용하면 기본 모드 자료들만의 역산에 비해 변수값들을 더 잘 구분해 낼 수 있고 깊은 구조, 특히 속도 역전이 일어난 구조에 대해 보다 나은 영상을 얻을 수 있음을 알 수 있다. 20 m 깊이에 저속도층이 존재하는 모델에서 이론적인 위상 속도들의 역산은 1차 모드만 첨가될 때 단지 연암층만을 영상화할 수 있다. 이 사실은 측정 가능한 가장 낮은 주파수가 단지 6 Hz 일 때 특히 중요하다. 현장시험들이 시추공과 PS 검층을 이용하여 조사된 지역들에서 행해졌다. 첫 번째 지역에서는 일본에서 심부 지질조사에 주로 이용되는 microtremor 배열 탐사가 35 m 깊이까지 토양층을 탐사하기 위해 사용되었다. 두 번째 지역에서는 12 m 깊이까지 조사하기 위해 sledgehammer 음원과 선형 다중 채널 수진기 전개를 이용하여 자료가 얻어졌다. 분산곡선 분석을 위해서 주파수-파수 파워 스펙트럼법이 사용되었고 각각의 시험에서 2차 모드의 속도까지 구해졌다. 다중 모드 역산 결과는 PS 검층기록과 잘 일치한다. 하지만 단지 기본 모드만을 이용하여 얻어진 결과는 매립지 아래의 천부 연암층까지의 깊이를 매우 작게 평가하였다.

• Journal of Forest and Environmental Science
• /
• 제28권1호
• /
• pp.24-29
• /
• 2012

This paper assesses the stmpage value of Korean white pine (Pinus koraiensis) stands in the Research Forest of Kangwon National University. Assessment was done by means of the cost value method (Faustmann formula) for age class I, Glaser method for age class II-V, and inversion market method for age class over VI. Generally the value of stumpage is calculated by the inversion market method in the market. However, immature trees and middle age class trees are not assessed in market, and the Korean forest is not old enough to harvest. So, when forests are damaged by attacks from forest fire, blight and harmful insects, the forest cannot be compensated for the loss from the government or insurance company. For this reason, the value of all-age class trees are calculated by using appropriate methods. As a result, the value of age I class stands (0.3 ha) is calculated as 1,786,305 (won), age II-Vclass stands (22.1 ha) 206,677,975 (won) and age VI class and over stands (24.8 ha) 523,789,603 (won).

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

Seismic exploration is divided by reflection and refraction method greatly, and reflection method can analyze complicated underground structure in the basis high resolution image, and refraction method can grasp the velocity structure of underground accurately. This thesis confirmed application of mixed exploration techniques using advantages of reflection and refraction. Reflection data processing applied conventional technique, and inversion of refraction data applied travel time tomographic technique that using SIRT method. Also, could establish initial information in model variable and improved the result of inversion by restricting model parameter value and dimension of area. Confirmed efficient fact in sequence and velocity structure grasping by utilizing accurate initial velocity model made out on the basis of marine reflection data, and mixed exploration technique using reflection and refraction have propriety that can trust in field application.

• 한국지구과학회지
• /
• 제26권7호
• /
• pp.691-697
• /
• 2005

최근에 수행되고 있는 중력 탐사는 고분해능의 중력계와 GPS(Global positioning system)를 통한 정밀한 측정과 측지가 이루어지고 있다. 중력탐사에서 모델링과 역산의 기술은 많은 발전이 있어왔지만, 중력자료처리는 거의 변화가 없었다. 통상적인 정밀한 중력 자료 보정을 통한 부우게이상은 측정점의 고도에서 기준면까지의 물질의 영향을 일정한 밀도를 이용해 제거해 버리기 때문에 측정점 바로 하부의 이상체에 의한 영향을 상당히 왜곡시키게 된다. 본 연구에서는 탐사 지역의 지형을 DEM(Digital Elevation Map) 자료와 Multiquadric equation을 이용하여 실제 지형과 유사한 Multiquadric surface를 자동적으로 구성하고, 이를 블록화 함으로써 보정의 대상이었던 기준면 상부에 대한 밀도를 탐사 지역의 지질 정보와 지형을 포함하는 역산을 통해 수치적으로 계산하였다. 이러한 지형을 포함한 역산 방법을 3차 원중력지형역산(3DGTI; 3-D Gravity Terrain Inversion)이라 한다. 이 연구의 효율성을 검증하기 위하여 주변암과 밀도차가 존재하는 관입지역에 대한 모델을 구성하고 적용한 결과 기존의 부게 보정 방법을 적용한 부게 이상도에 비해 자료의 왜곡이 감소하는 효과를 얻을 수 있었다. 이를 통하여 지형 역산을 통한 객관적인 부게 밀도의 결정과 부게 보정시 실제의 수평적인 밀도 변화를 반영함으로써 기존의 문제점을 보완하였다. 게다가, 3DGTI로부터 얻어진 밀도분포는 지형의 윤곽을 그대로 표현하고 있어서 보다 실질적인 지질을 보여준다고 하겠다. 이 방법을 화강암체가 관입하고 있는 마산$\cdot$창원 일대에서의 중력 탐사 자료에 적용해본 결과 기존 방법보다 관입 화강암체의 위치와 그 규모를 알아내는데 더 효과적이었다. 따라서, 수평적인 밀도 변화가 뚜렷하게 존재하는 지역의 경우, 새로운 중력 자료 처리 방법이 기존의 부게 보정에서 발생하였던 문제점을 해결함으로써 천부의 분해능을 높이고, 심부의 밀도 분포도 좀더 정확하게 계산할 수 있으리라 생각된다.

• 지구물리와물리탐사
• /
• 제18권2호
• /
• pp.54-63
• /
• 2015

매장량의 평가에 직접적인 정보를 제공해주는 저류층 변수들인 공극률과 유체포화율은 물리탐사방법을 통해 직접적으로 획득이 가능한 탄성파 속도나 전기비저항의 물성값을 암석물리모델 구성법에 적용하여 추정이 가능하다. 따라서, 정확한 저류층 변수들의 추정을 위해서는 우선적으로 정확한 속도나 전기비저항과 같은 물성값들의 추정이 필요하다. 이종의 물리탐사자료를 이용한 복합역산은 단일 물리탐사자료를 이용한 역산과 비교시, 역산의 불확실성을 줄일 수 있고, 두 탐사자료의 장점을 함께 이용할 수 있기 때문에 단일물리탐사자료에 비하여 보다 신뢰성있는 물성정보를 추정할 수 있다. 이 연구에서는 효율적인 공극률과 유체포화율의 분포를 추정하기 위하여, 탄성파탐사자료와 전자탐사자료의 복합역산을 통해서 지하의 속도 정보와 전기비저항 정보를 획득한 뒤, 이 두 물성을 모두 이용하여 암석물리모델을 구성하는 과정이 포함 된 격자탐색법(grid-search)을 제안하였다. 오일저류층 모델의 합성자료에 적용한 결과, 보다 신뢰성 있는 저류층 변수들을 추정하였으며, 이는 오일저류층의 정확한 매장 위치 추정과, 매장량 계산에 보다 정확한 정보를 제공해 줄 것으로 기대된다.

• 자원환경지질
• /
• 제14권3호
• /
• pp.193-201
• /
• 1981

수평다층(水平多層) 지질구조(地質構造)모델에 있어서 외견전기비저항곡선(外見電氣比抵抗曲線)의 자동해석(自動解析)에 최소자승법(最小自乘法)을 적용(適用)하여 보았다. 이 방법(方法)은 digital filtering 법(法)과 종합(綜合)된 damped least-squares algorithm으로 구성된 것으로서, 일반적(一般的)으로 사용(使用)되고 있는 curve-matching 법(法)보다 시간(時間)이 빠르고 정도(精度)도 높다는 것을 알게되었다. 이 반전법(反轉法)을 시험(試驗)하기 위해서 한개의 이론(理論)데이터와 3개의 현장탐사(現場探査) 결과(結果)를 선택(選擇)하여 해석(解析)하였다. 3층(層) 지질구조(地質構造)모델로부터 나오는 이론(理論)데이터의 해석(解析)을 통(通)하여 이 방법(方法)의 특징(特徵)을 파악(把握)할 수 있고, 또 종래(從來)의 Curve-matching법(法)에는 없는 명확(明確)한 특징(特徵)이 나타났다. 더우기 반전법(反轉法)의 유효성(有效性)은 현장탐사(現場探査) 결과(結果)의 해석(解析)에 의해서도 알 수 있었다. 그 때의 최적(最適) 지질구조(地質構造)모델은 시굴(試掘)로 확인(確認)된 지하구조(地下構造)와 일치(一致)함을 알 수 있었다.

• 지구물리와물리탐사
• /
• 제10권2호
• /
• pp.124-130
• /
• 2007

탄성파 역산에 있어서 가장 널리 사용되는 최소자승($l^2$ norm)해는 이상치(outlier)에 매우 민감하게 반응하는 경향이 있다. 이에 반해서 $l^1$ norm을 최소화하는 해는 이상치에 강인한 면을 보이나 일반적으로 좀 더 많은 계산이 필요하다. 반복적가중의 최소자승법(Iteratively reweighted least squares [IRLS] method)을 이용하면 이러한 $l^1$ norm 문제의 근사해(approximate solution)를 효율적으로 구할 수 있다. 본 논문에서는 작은 크기의 잔여분은 $l^2$ norm으로 처리하며, 큰 크기의 잔여분은 $l^1$ norm으로 처리하는 하이브리드 $l^1/l^2$ norm 최소화를 IRLS 방법에 쉽게 적용하는 구현 기법을 소개한다. 소개된 알고리즘은 특이치(singularity)처리를 위한 임계값의 결정에 민감하게 반응하는 기존의 $l^1$ norm IRLS 방법과는 달리 임계값 결정에 상관없이 늘 강인한 역산의 특성을 보여준다.

• 한국전자파학회논문지
• /
• 제12권7호
• /
• pp.1157-1166
• /
• 2001

본 논문에서는 분산 천이 광섬유를 전송로로 채택한 10 Gbps, 20 Gbps, 40 Gbps 전송 시스템에서의 전력 대칭 MSSI(mid-span spectral inversion)에 의한 보상 정도를 변조된 광 펄스의 다양한 첩 파라미터에 따라 분석하였다. 우선 각각의 전송 속도에서 입력 전력 변화에 따른 수신단에서의 EOP(eye-opening penalty)를 계산하여 수신 성능을 양호하게 유지할 수 있는 최대 입력 전력의 크기를 첩 파라미터에 따라 살펴보았다. 또한 MSSI의 장거리 광대역 WDM 전송 시스템에의 적용 가능성을 확인해 보기 위하여 송신단부터 광 위상 공액기(OPC; optical phase conjugator)가지의 첫 번째 광섬유의 분산 계수 D$_{11}$ 변동에 따른 EOP의 고찰을 통해 수신 성능이 양호하게 유지될 수 있는 송신 과장의 범위를 살펴보았다. 본 논문에서 제안된 최적 펌프 전력 조건을 유지하는 MSSI 방법은 이상 분산(anomalous dispersion) 영역에서 변조 과정을 통해 광 펄스에 인가된 초기 첩이 up-chirp 인 경우보다 down-chirp인 경우에서 더욱 효과적임을 확인할 수 있었고. 장거리 WDM 전송에서 비트율에 따라 3.5 dBm 이상의 비교적 높은 전력으로 수~수 십 nm 이상의 광대역 전송이 가능하다는 것을 확인할 수 있었다.다.