• 제목/요약/키워드: numerical modeling and inversion

검색결과 48건 처리시간 0.024초

전기비저항 탐사자료의 3차원 역산 (Three-dimensional Inversion of Resistivity Data)

  • 이명종;김정호;조성준;정승환;송윤호
    • 지구물리와물리탐사
    • /
    • 제2권4호
    • /
    • pp.191-201
    • /
    • 1999
  • 현재까지 전기비저항 탐사자료의 해석은 대부분 지하구조가 1차원 또는 2차원으로 가정하는데 근간을 두고 있다. 그러나, 국내와 같이 지질구조가 복잡하며 지형의 기복이 심한 경우, 이와 같은 가정을 사용은 전기비저항 탐사자료의 해석에 오류를 야기할 가능성이 높다. 본 연구에서는 3차원적으로 복잡하게 변화하는 지하구조 및 지형의 기복에 의한 지형효과를 정확히 계산할 수 있는 유한요소법에 의한 3차원 전기비저항 모델링 알고리듬을 근간으로 하는 역산 알고리듬을 개발하였다. 이 알고리듬은 ACB(Active Constraint Balancing)법을 채용함으로써 역산의 분해능 향상을 기하였다. 간단한 3차원 지하구조 및 지형기복에 대한 수치모형 실험결과 지표의 불규칙한 기복을 포함하여 역산을 수행함으로써 정확한 지하구조 영상을 획득할 수 있음을 보였으며, 국내의 교량건설 예정부지에서 수행된 격자망을 이룬 전기비저항 탐사자료에 대한 3차원 역산 결과 매우 신뢰도 높은 지하구조에 대한 영상을 도출할 수 있었다.

  • PDF

SAR Imaging에서 Physical Array와 합성 Array 신호의 Subspace 비교를 통한 Wavefront Inversion 기법 입증 (Verification of Wavefront Inversion Scheme via Signal Subspace Comparison Between Physical and Synthesized Array Data in SAT Imaging)

  • 최정희
    • 전자공학회논문지D
    • /
    • 제36D권4호
    • /
    • pp.34-41
    • /
    • 1999
  • Synthetic Aperture Radar 시스템은 레이더 또는 목적물을 일정거리 또는 시간 간격동안 움직여가면서 반사되는 신호를 측정하여 동시에 처리해 줌으로써 안테나 크기를 크게 하는 효과를 가지게 되므로 목적물의 형체를 파악하는 Imaging 시스템의 일종이다. Imaging의 성패는 어떤 Modeling 과 Inversion 과정을 거치느냐에 달려 있는데 본 논문에서는 가장 Approximation이 적은 Wavefront reconstruction 기법을 적용하여 SAR 모델에서 획득한 데이터와 같은 크기의 실제 안테나가 있다고 가정하였을 때의 데이터와 이론적이고 실험적인 Signal Subspace 비교를 수행하였다. 이 연구를 통해 새로운 기법인 Wavefront reconstruction 알고리즘을 충분히 뒷받침해 줄 수 있는 결과를 얻을 수 있다.

  • PDF

Chirp 신호를 이용한 해저퇴적층의 음향학적 특성 역산 (Inversion of Acoustical Properties of Sedimentary Layers from Chirp Sonar Signals)

  • 박철수;성우제
    • 한국음향학회지
    • /
    • 제18권8호
    • /
    • pp.32-41
    • /
    • 1999
  • 본 논문은 chirp신호와 두 개의 근거리 청음기를 이용한 해저퇴적층의 음향학적 특성치 역산기법을 제시한다. 역산문제를 확률론적 모델로 정식화하고, 역산의 해를 역산인자의 a priori분포와 유사도함수의 곱으로 표현되는a posteriori 확률분포로 정의하였다. 퇴적층의 음속과 층두께의 a priori정보를 파형 매칭 기법으로 추정한 후 다수의 퇴적층이 존재하는 환경모델을 부분퇴적층모델로 치환하고, 계측신호와 모의신호의 L₂노음을 이용하여 정의된 목적함수에 대해 반복적인 유전자알고리즘 탐색을 수행하여 탐색공간의 축소로 인한 탐색효율과 결과의 향상을 얻었다. A posteriori 확률분포의 다중적분의 형태로 정의되는 인자의 주변확률분포와 평균의 추정은 유전자알고리즘의 탐색과정에서 선택된 탐색점들을 이용하여 수행되었다. 제시된 역산기법의 검증을 위해 두 가지 퇴적층 환경모델을 설정하고 잡음을 첨가한 합성신호에 대해 역산기법을 적용하여 역산해를 추정하였고 역산결과로부터 본 역산기법의 유용성을 확인하였다.

  • PDF

고정밀 중력탐사 자료 역산 패키지 (An Inversion Package for Interpretation of Microgravity Data)

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

지형효과를 포함한 3차원 전기비저항 역산 (Three-dimensional Resistivity Inversion Including Topographic Effect)

  • 박종오;김희준;송무영
    • 지질공학
    • /
    • 제14권1호
    • /
    • pp.21-28
    • /
    • 2004
  • 지형효과를 포함하는 3차원 전기비저항 역산법은 전기비저항 탐사자료의 해석에 있어서 실제 지하구조에 존재하는 3차원 전도성 이상체의 특성을 파악함으로써 원리적으로 정확한 영상을 획득할 수 있는 방법이라고 본다. 본 연구는 수치모델링에서 유한요소법을 이용하였으며, 역산에서 설정된 블록에 대하여 감도 분석을 통하여 자코비안 계산의 효율성을 극대화하였다. 또한 역산과정에서는 라그랑지 곱수의 값을 변화시키면서 최소자승근 오차의 최소가 되는 최적의 값을 선택하는 방법으로 역산의 분해능을 향상 시켰다. 전기비저항 탐사에서 지형효과의 해석으로는 지형 의 기복을 포함하여 얻어진 순차모델링 자료를 가지고 역산 방법에서 지형의 기복을 무시하여 얻어진 결과와 지형의 기복을 포함하여 얻어진 결과를 비교하여 불규칙한 지형효과에 의한 왜곡된 전기장 반응으로 인한 오해석을 파악하는데 있다.

파동방정식 수치해의 일관성에 관한 연구 (A Study on Consistency of Numerical Solutions for Wave Equation)

  • 편석준;박윤희
    • 지구물리와물리탐사
    • /
    • 제19권3호
    • /
    • pp.136-144
    • /
    • 2016
  • 탄성파 자료의 역산은 파동방정식에 기초하고 있으므로 파동방정식의 해를 정확하게 구하는 것이 가장 중요하다. 특히, 전파형역산은 파동장 전체를 이용하기 때문에 정문제에 해당하는 모델링이 정확하게 이루어져야 신뢰할 수 있는 결과를 얻게 된다. 파동방정식의 수치해를 구하는 대표적인 기법인 유한차분법과 유한요소법은 해의 수렴성을 보장할 수 있어야 하는데, 해의 수렴성은 이론적으로 일반화된 증명이 되어 있으나 실제 문제에 적용할 경우 일관성과 안정성을 분석해야 한다. 모델링 결과의 일관성은 송신원 함수의 구현이 매우 중요한 부분인데, 유한차분법은 디랙 델타 함수(Dirac delta function)를 나타낼 때 격자 간격으로 표준화된 싱크 함수(sinc function)를 사용해야 하는 반면 유한요소법은 격자 간격에 관계없이 기저함수 값을 사용하면 된다. 주파수 영역 파동방정식을 사용할 경우 송신 파형 함수의 스펙트럼을 정확하게 표현하기 위해 샘플링 이론으로 정의되는 시간 간격보다 더 조밀한 샘플링 간격을 사용하고 나이퀴스트(Nyquist) 주파수보다 더 높은 주파수를 최대 주파수로 사용해야 한다. 또한, 복소 각주파수를 사용하는 경우 감쇠 파동방정식을 만족하기 위해서는 송신 파형 함수를 먼저 감쇠한 후 사용해야 한다. 이러한 요건들이 모두 만족되었을 때 신뢰할 수 있는 역산 알고리즘 개발이 가능하다.

역 사상법에 의한 시간영역 유도분극 자료의 역산 (Inversion of Time-domain Induced Polarization Data by Inverse Mapping)

  • 조인기;김연정
    • 지구물리와물리탐사
    • /
    • 제24권4호
    • /
    • pp.149-157
    • /
    • 2021
  • 유도분극 탐사와 전기비저항 탐사는 자료획득이 유사하며, 대부분의 전기비저항 탐사 시스템에는 시간영역 유도분극 탐사 기능이 함께 탑재되어 있다. 또한 시간영역 유도분극 탐사 자료에는 전기비저항 자료가 내포되어 있다. 따라서 유도분극 탐사와 전기비저항 탐사와는 불가분의 관계가 있으며, 유도분극 자료의 역산도 전기비저항 탐사 자료의 역산에 근거한 2단계 역산법이 적용되고 있다. 그러나 유도분극 탐사는 효과적인 해석법의 부재로 인하여 전기비저항 탐사에 비하여 널리 적용되지 못하고 있다. 이 연구에서는 수치 모델링 및 역산실험을 통하여 시간영역 유도분극 자료의 역산해석에 사용되는 역 사상법의 문제점을 분석하였다. 또한 역 사상법 적용시 문제가 되는 역산잡음을 효과적으로 억제할 수 있는 수정된 역 사상법을 제시하였다. 마지막으로 수치자료에 대한 역산실험을 통하여 개발된 역 사상법의 효과를 검증하였다.

음의 겉보기 충전성 자료를 포함한 시간영역 유도분극 자료의 비선형 역산 (Nonlinear Inversion of Time-domain Induced Polarization Data with Negative Apparent Chargeability Data)

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

Transient Analysis of Hybrid Systems Composed of Lumped Elements and Frequency Dependent Lossy Disributed Interconnects

  • Ichikawa, Satoshi;Shimoda, Tomokazu
    • 대한전자공학회:학술대회논문집
    • /
    • 대한전자공학회 2000년도 ITC-CSCC -2
    • /
    • pp.1096-1099
    • /
    • 2000
  • A method to analyze the high speed inter-connects that are composed of frequency dependent lossy distributed lines is presented. Network modeling of hybrid systems is implemented by using the modified nodal admittance matrix in the Laplace transformation domain. The network response is computed by different two methods. One method Is the asymptotic waveform evaluation (AWE) method and other is numerical Laplace inversion method. The merits and demerits of two methods are discussed by applying to several concrete illustrative networks.

  • PDF

물리 기반 유한 단층 미끌림 역산을 위한 CPInterface (COMSOL-PyLith Interface) 개발 (Development of a CPInterface (COMSOL-PyLith Interface) for Finite Source Inversion using the Physics-based Green's Function Matrix)

  • 김민수;소병달
    • 지구물리와물리탐사
    • /
    • 제26권4호
    • /
    • pp.268-274
    • /
    • 2023
  • 유한 단층 미끌림 역산에는 지진 변위 측지 자료와 그린 함수 행렬(Green's function matrix)을 주로 사용한다. 그린 함수 행렬은 일반적으로 오카다 모형(Okada, 1985)을 기반으로 한다. 그러나 최근 물리 기반 지진 모델링을 활용하여 그린 함수 행렬을 제작하고 유한 단층 미끌림 역산을 수행하는 연구가 활발하다. 물리 기반 지진 모델링은 다양한 물성(탄성, 점탄성, 탄소성 등)을 고려하여 현실적인 환경에서 지진을 모사할 수 있다는 장점이 있다. 물리 기반 유한요소 소프트웨어 PyLith는 단층을 구성하는 절점을 두 개로 나누어 지진을 모사할 수 있으므로 지진 모사 모델링에 적합하다. 하지만 PyLith는 격자망 생성 기능을 자체 제공하지 않아, 모형 내부에 수십~수백 개의 소단층과 관측점을 설정해야 하는 유한 단층 미끌림 역산 수행에는 어려움이 있다. 본 연구에서는 소단층과 관측점을 포함한 수치 모형을 제작하고, 지진 모사 모델링을 수행하여 그린 함수 행렬을 제작하는 일련의 과정을 연계하여 유한 단층 미끌림 역산의 편리성을 높이기 위해 CPInterface (COMSOL-PyLith Interface)를 개발하였다. CPInterface는 COMSOL의 격자 생성 능력과 PyLith의 지진 모사 능력을 결합하여 그린 함수 행렬을 자동으로 생성할 수 있다. CPInterface는 간단한 변수들로 모형 및 단층 정보를 조절할 수 있고, 지하 탄성 이상체와 GPS 관측점을 자유롭게 배치할 수 있다. 또한, 그린 함수 행렬을 생성하는 복잡한 과정을 간소화하여 더욱 편리하게 유한 단층 미끌림 역산을 할 수 있게 한다.