• 제목/요약/키워드: Inversion method

검색결과 876건 처리시간 0.021초

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

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

분해능 모델 제한자를 사용하는 자력탐사자료의 2차원 역산 (2D Inversion of Magnetic Data using Resolution Model Constraint)

  • 조인기;강혜진;이근수;고광범;김종남;유영준;한경수;신홍준
    • 지구물리와물리탐사
    • /
    • 제16권3호
    • /
    • pp.131-138
    • /
    • 2013
  • 새로운 2차원 자력탐사자료 역산 방법을 개발하였다. 중,자력탐사와 같은 포텐셜 자료의 역산에서 가장 문제가 되는 점은 비유일해 문제이다. 일반적으로 자력탐사 자료의 역산은 모델변수의 수가 자료의 수보다 훨씬 많은 불충분 문제이며, 이는 비유일해 문제를 더욱 심화시키게 된다. 일반적인 최소자승법을 자력탐사자료의 역산에 적용하게 되면, 이 상대가 지표면에 집중되는 결과를 초래한다. 본 연구에서는 이러한 비유일해 문제를 극복하기 위하여 모델분해능에 근거한 새로운 모델제한자를 제안하였다. 이 모델제한자는 분해능이 높은 모델변수에는 큰 제한을 가하고, 작은 모델변수에는 약한 제한을 가하게 된다. 따라서 분해능이 낮은 심부의 모델변수도 효과적으로 추정할 수 있다. 개발된 역산 알고리듬을 이용하여, 전형적인 모델에 대한 이론자료의 역산에 적용하였다. 또한 옥천대에서 얻어진 항공자력탐사자료 역산에 적용하였다.

지구물리 자료의 고속 베이지안 역산 (Fast Bayesian Inversion of Geophysical Data)

  • 오석훈;권병두;남재철;이덕기
    • 지구물리
    • /
    • 제3권3호
    • /
    • pp.161-174
    • /
    • 2000
  • 베이지안 역산(Bayesian inversion)은 불충분한 자료를 가지고 지하구조를 추정해야 하는 지구물리자료의 해석에 있어서 안정적이고 신뢰를 줄 수 있는 방법 중의 하나이다. 관측 자료가 측정 과정부터 불확실성을 함유하고 있으며, 역산에 이용되는 이론 자료 또한 모델의 매개변수화에 따른 각종 불확실성을 포함하고 있다. 따라서 지구물리 자료의 역산은 확률적으로 접근하는 것이 가장 바람직하며 베이지안 역산은 이에 대한 처리뿐만 아니라, 추정에 대한 신뢰도와 불확실성에 대한 이론적 근거를 제공한다. 그러나 대부분의 베이지안 역산이 고차원의 적분을 필요로 하므로 몬테 카를로 방법과 같은 대규모의 계산이 요구되는 방법에 의해 사후 확률분포가 구해지는 경우가 많다. 이는 특히 지구물리 자료와 같이 고도의 비선형 자료에 대하여 매우 적합한 접근 방법이기는 하지만, 점차 현장화, 고속화되어가는 자료의 해석 경향에 맞추어 간략하게 사후 확률분포를 근사한 수 있는 기법의 연구 또한 필요하다. 따라서 이 연구에서는 관측자료와 사전 확률분포가 정규분포에 의해 근사 될 수 있는 지구물리자료에 대한 베이지안 역산에 대해 논의 하고자 한다. 사전 확률분포의 작성을 위해 지구통계학적 기법이 이용되었으며, 관측자료의 통계적 불화실성을 추정하기 위해 교차 검사(cross-validation) 방법을 이용하여 공분산(covariance)을 유도하고 그것에 의한 우도 함수(likelihood function)를 작성하였다. 베이지안 해석을 위해 두 확률분포를 곱하여 근사적인 사후 확률분포를 얻을 수 있었으며, 이에 대해 최적화(optimization) 기법을 이용하여 최대 사후 확률(Maximum a Posterior)을 따르는 지하 구조를 얻을 수 있었다. 또한 사후 확률 분포의 공분산 항을 이용하여 지하 비저항 구조를 시뮬레이션 하여 불확실성분석을 수행하였다.

  • PDF

MT 자료의 3차원 역산 개관 (Review on the Three-Dimensional Inversion of Magnetotelluric Date)

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

불균질 횡등방성 매질에서의 탄성파 주시토모그래피 (Seismic Traveltime Tomography in Inhomogeneous Tilted Transversely Isotropic Media)

  • 정창호;서정희
    • 지구물리와물리탐사
    • /
    • 제10권4호
    • /
    • pp.229-240
    • /
    • 2007
  • 이 연구에서는 지하 매질의 속도 이방성을 규명하고 이를 고려하여 지하 단면을 해석하는 탄성파 이방성 토모그래피(seismic anisotropic tomography) 알고리듬을 개발하였다. 이는 지하 매질을 대칭축이 기울어진 횡등방성(Tilted Transversely Isotropic) 매질로 가정하여 지하의 속도 구조와 더불어 이방성 구조를 함께 파악할 수 있도록 고안되었다. 개발된 알고리듬은 이방성 토모그래피에서 흔히 나타날 수 있는 역산의 안정성 문제를 해결하기위해여 역산 매개 변수를 역속도비, 역속도 그리고 이방성 대칭축 방향으로 설정하였고, 의사 베타 변환(pseudo-beta transform)을 통해 이방성 계수의 범위를 제한하였으며, ACB (Active Constraint Balancing) 법을 이용하여 효율적으로 역산 제한조건을 적용하였다. 특히, 기존의 등방성 토모그래피에서 사용하던 프레넬대(Fresnel volume)를 고려한 역산 기법을 이방성 토모그래피로 확장하여 파선 토모그래피에 비해 넓은 범위의 전파각 영역을 확보함으로써 이방성 토모그래피의 안정성을 높일 수 있었다. 개발된 이방성 토모그래피 알고리듬을 복잡한 수치 모델에 적용하여 알고리듬의 타당성을 검증하였으며 현장 자료에 적용하여 시추 결과 및 지질 정보들과 일치하는 향상된 결과를 얻었다. 개발된 이방성 토모그래피 알고리듬을 통해 지하매질의 속도 정보와 이방성 정보를 함께 고려함으로써 보다 정확한 지하 구조의 해석이 가능할 것으로 기대된다.

퇴적층들의 전단파 속도 평가를 위한 수신함수와 표면파 위상 속도의 통합 역산 (Joint inversion of receiver function and surface-wave phase velocity for estimation of shear-wave velocity of sedimentary layers)

  • ;산중호명
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.93-101
    • /
    • 2006
  • 이 연구에서는 유전자 알고리듬을 이용하여 수진기 함수와 표면파 위상 속도로부터 심부 퇴적층의 전단파 속도 구조를 결정하는 복합역산 방법을 제시하였다. 합성 탄성파 자료를 이용한 수치모형실험은 제시된 방법이 단지 수진기 함수만을 역산했을 때 발생하는 전단파 속도와 층 두께 사이의 trade-off와, 표면파 위상 속도만을 역산했을 때 심부구조에서의 불확실성을 피할 수 있음을 보여주고 있다 이 방법은 진앙거리 100km의 지진기록으로부터 얻은 수진기 함수들과 일본 Kanto 평원의 상시진동 배열 탐사로부터 얻은 레일리파의 위상속도에 적용되었으며, 추정된 심부구조는 선행된 굴절법 탐사 결과 및 심부 시추공 자료와 잘 일치하였다.

천부지각 2차원 속도구조를 위한 레일리파의 군속도와 위상속도 역산의 비교 연구 (A Study of the comparison of Inversion of Rayleigh wave Group and Phase Velocities for Regional Near-Surface 2-Dimensional Velocity Structure)

  • 이보라;정희옥
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2006년도 공동학술대회 논문집
    • /
    • pp.51-59
    • /
    • 2006
  • 서해안 조간대에서 24 channel 탐사기로 획득한 표면파 자료의 위상속도와 군속도를 구한다음, 이를 역산하여 그 결과를 비교하였다. 위상속도의 분산곡선은 tau-p stacking 방법에 의하여, 군속도의 분산곡선은 wavelet analysis와 Multiple Filtering Technique의 두가지 방법을 사용하여 구하였다. 위상속도의 오차가 군속도의 오차보다 더 큰 것을 확인하였다. 군속도의 경우, wavelet analysis가 Multiple Filtering Technique 보다 fundamental mode와 higher mode를 구분하는데 더 효과적이었다. 역산결과, 군속도의 fundamental mode와 1st higher mode 를 동시에 사용했을 때, 공간적 해상도가 가장 좋았다. 이연구는 천부 지반의 S파 속도 구조를 구하는데, 군속도의 higher mode를 포함한 군속도 분산곡선을 사용하는 것이 효과적임을 시사한다.

  • PDF

자료 가중을 통한 전기비저항 탐사 자료의 역산 (Inversion of Resistivity Data using Data-weighting)

  • 조인기;이근수;김연정;윤대성
    • 지구물리와물리탐사
    • /
    • 제18권1호
    • /
    • pp.9-13
    • /
    • 2015
  • 전기비저항 탐사 자료는 다양한 잡음을 포함하고 있다. 즉 전기비저항 자료는 높은 접촉저항, 장비의 측정 오차 및 주변의 불규칙한 전기적 잡음에 의해 영향을 받는다. 전기비저항 탐사 자료의 올바른 해석을 위해서는 이들 잡음의 정확한 추정이 요구된다. 이 연구에서는 상반성 시험을 통하여 추정된 잡음을 역산시 자료 가중에 반영하는 방법론을 제안하였다. 또한 역산시 현장 자료와 이론 자료 사이의 적합 오차와 상반성 오차를 분석하고, 상반성 오차와 적합 오차를 모두 이용하는 자료 가중법을 제안하였다. 현장 자료에 제안된 자료 가중법을 적용한 결과 통상적인 역산 결과에 비하여 국지적 이상대의 출현 빈도가 감소하는 것을 확인할 수 있었다.

가스 하이드레이트 부존층의 구조파악을 위한 탄성파 AVO 분석 AVO모델링, AVO역산 (Seismic AVO Analysis, AVO Modeling, AVO Inversion for understanding the gas-hydrate structure)

  • 김건득;정부흥
    • 한국신재생에너지학회:학술대회논문집
    • /
    • 한국신재생에너지학회 2005년도 춘계학술대회
    • /
    • pp.643-646
    • /
    • 2005
  • The gas hydrate exploration using seismic reflection data, the detection of BSR(Bottom Simulating Reflector) on the seismic section is the most important work flow because the BSR have been interpreted as being formed at the base of a gas hydrate zone. Usually, BSR has some dominant qualitative characteristics on seismic section i.e. Wavelet phase reversal compare to sea bottom signal, Parallel layer with sea bottom, Strong amplitude, Masking phenomenon above the BSR, Cross bedding with other geological layer. Even though a BSR can be selected on seismic section with these guidance, it is not enough to conform as being true BSR. Some other available methods for verifying the BSR with reliable analysis quantitatively i.e. Interval velocity analysis, AVO(Amplitude Variation with Offset)analysis etc. Usually, AVO analysis can be divided by three main parts. The first part is AVO analysis, the second is AVO modeling and the last is AVO inversion. AVO analysis is unique method for detecting the free gas zone on seismic section directly. Therefore it can be a kind of useful analysis method for discriminating true BSR, which might arise from an Possion ratio contrast between high velocity layer, partially hydrated sediment and low velocity layer, water saturated gas sediment. During the AVO interpretation, as the AVO response can be changed depend upon the water saturation ratio, it is confused to discriminate the AVO response of gas layer from dry layer. In that case, the AVO modeling is necessary to generate synthetic seismogram comparing with real data. It can be available to make conclusions from correspondence or lack of correspondence between the two seismograms. AVO inversion process is the method for driving a geological model by iterative operation that the result ing synthetic seismogram matches to real data seismogram wi thin some tolerance level. AVO inversion is a topic of current research and for now there is no general consensus on how the process should be done or even whether is valid for standard seismic data. Unfortunately, there are no well log data acquired from gas hydrate exploration area in Korea. Instead of that data, well log data and seismic data acquired from gas sand area located nearby the gas hydrate exploration area is used to AVO analysis, As the results of AVO modeling, type III AVO anomaly confirmed on the gas sand layer. The Castagna's equation constant value for estimating the S-wave velocity are evaluated as A=0.86190, B=-3845.14431 respectively and water saturation ratio is $50\%$. To calculate the reflection coefficient of synthetic seismogram, the Zoeppritz equation is used. For AVO inversion process, the dataset provided by Hampson-Rushell CO. is used.

  • PDF

Time domain and frequency domain interpretation of safety diagnosis for concrete structure

  • Suh Baeksoo;An Jehun;Kim Hyoungjun;Kim Yongin
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 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.

  • PDF