• 제목/요약/키워드: seismic inversion

검색결과 170건 처리시간 0.023초

지하공동 충전효과 평가를 위한 시차 공대공 탄성파 토모그래피 연구 (Time-Lapse Crosswell Seismic Study to Evaluate the Underground Cavity Filling)

  • 이두성
    • 지구물리와물리탐사
    • /
    • 제1권1호
    • /
    • pp.25-30
    • /
    • 1998
  • 공동충전 효과를 검증하기 위하여 실시한 시차 공대공 탄성파 탐사자료로부터 지하공동 부존 지역에서 충전 전과 후에 매질의 탄성파 전파속도의 변화를 확인하였다. 시차 공대공 탄성파 탐사자료에 나타난 반응과 시추조사 결과에 의하면 본 지역의 공동은 규모가 극히 소규모이거나 또는 폐석 등으로 충전된 것으로 보인다. 공동충진 효과는 토모그래피로부터 도출된 속도단면상의 탄성파 속도의 증가량을 분석함으로써 평가하였다. 시추공용 에어건을 진원으로 24-채널 하이드로폰을 수진기로 하여 자료를 취득하였다. 취득한 자료에는 무시할 수 없을 정도의 source statics를 확인할 수 있었다. 본 논문에서 제시한 보정방법은 2단계로; 1) 불규칙한 발파시점에 의한 영향 보정과 2) 잔여 정보정으로 이는 진원의 부정확한 위치에 대한 정보정이다. 본 논문에서는 고주파수 성분의 수치잡음이 억제되고 관심대상 부분에서 비교적 고분해능 영상을 도출할 수 있는 다단계 역산 방안을 제시하였다. 일반적으로 최소자승 주시토모그래피로는 평활화된 속도 영상을 얻을 수 있다. 따라서 이러한 역산으로는 비교적 소규모의 구간에서 발생한 적은 속도변화를 영상화하기에는 어려운 면이 있다. 본 논문에서는 속도모델의 파라메터를 변화시킨 2단계 제어 역산법으로 도출한 시차 토모그램으로부터 채굴 영향대에서 발생한 매질의 속도변화를 시각화 할 수 있었다. 2단계 역산법은 1-단계에서는 적정한 크기의 균일 격자로 구성된 모델을 사용하여 토모그램을 작성하고 이 토모그램에 2차원 중위수 필터를 적용하여 대략적인 속도구조 모델을 작성한다. 2-단계 역산시는 1-단계에서 작성한 속도모델을 수정하여 초기 모델로 한다. 모델 수정은 관심대상 부분만을 작은 크기의 균일격자로 재구성하는 것이다. 기준조사 토모그램을 2차 조사자료 역산의 초기 속도모델로 사용하였다. 속도변화는 공동대 부근에서만 예상되므로 그 이외 부분의 속도는 기준 토모그램과 동일하게 고정시키고 역산을 수행하였다.

  • PDF

호주 Victoria주 MT 탐사 자료 해석 (Analysis of MT Data Acquired in Victoria, Australia)

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

한반도의 소규모지진 모멘트 텐서 역산의 응용 (Application of Moment Tensor Inversion to Small Local Earthquakes in the Korean Peninsula)

  • 김소구;판티김반;이승규
    • 한국방재학회 논문집
    • /
    • 제1권3호
    • /
    • pp.123-136
    • /
    • 2001
  • 모멘트 텐서 역산의 응용 목적은 지진 메커니즘, 지진 모멘트, 한반도에서 국지적으로 발생하는 규모 5.5 이하의 소규모 지진의 진원깊이 등의 진원인자를 결정하는 것이다. 이와 같이 광범위한 지형에서 발생하는 지진들은 어느 지역에서의 활동적인 지체구조를 밝혀내는 실마리를 제공해 준다. 더욱이 소규모 지진의 지진 메커니즘 연구를 위해서는 반드시 필요하다. 본 논문에서는 한반도에서 최근에 발생한 규모 3.3 - 4.8의 14개 지진을 모멘트 역산방법 (TDMT_INV)을 이용하여 지진 메커니즘 분석 결과를 소개하고 있다. 단층들의 주향은 북동-남서, 북서-남동 방향으로 주향이동 또는 소규모 역단층의 주향이동 단층을 보이고 있다. 응력장의 압축방향은 몇몇 단층을 제외하고는 동-서 또는 동북동-서남서 방향이 우세한 것으로 나타났다. 그리고 영남 육괴와 원산 지역에 분포하는 지진들은 남-북 또는 북북서-남남동 방향이 우세한 것으로 나타났다.

  • PDF

퇴적층들의 전단파 속도 평가를 위한 수신함수와 표면파 위상 속도의 통합 역산 (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 평원의 상시진동 배열 탐사로부터 얻은 레일리파의 위상속도에 적용되었으며, 추정된 심부구조는 선행된 굴절법 탐사 결과 및 심부 시추공 자료와 잘 일치하였다.

강인한 역산으로서의 하이브리드 $l^1/l^2$ norm IRLS 방법의 효율적 구현기법 (An Efficient Implementation of Hybrid $l^1/l^2$ Norm IRLS Method as a Robust Inversion)

  • 지준
    • 지구물리와물리탐사
    • /
    • 제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 방법과는 달리 임계값 결정에 상관없이 늘 강인한 역산의 특성을 보여준다.

확장 탄성 임피던스 역산을 이용한 저류층 물성 예측 (Prediction of Reservoir Properties Using Extended Elastic Impedance Inversion)

  • 김현주;이광훈;문성훈
    • 자원환경지질
    • /
    • 제48권2호
    • /
    • pp.115-130
    • /
    • 2015
  • 확장 탄성 임피던스(extended elastic impedance, EEI)는 입사각에 따른 음향 임피던스(acoustic impedance, AI) 를 일반화한 탄성 임피던스(elastic impedance, EI)를 확장한 개념으로서 다양한 저류암 물성과 대비가 가능한 것으로 알려져 있다. 하지만 EEI 역산을 적용하여 예측한 저류암적 물성이 실제 물성을 얼마나 정확히 예측하는지를 검증한 사례는 거의 없다. 본 연구에서는 EEI 역산 기법을 이용하여 미국 와이오밍주 Teapot Dome 유전의 주요 저류층 중 하나인 Second Wall Creek 사암층의 P파속도($V_p$), S파속도($V_s$), P파속도-S파속도 비($V_p/V_s$), 포아송비(Poisson's ratio)와 같은 속도 물성들을 유추하고 실제 물리검층 자료와 비교하여 EEI 역산 기법의 정확도를 검증했다. 사용된 자료는 Teapot Dome 유전 남부의 3차원 공심점 모음 자료(CDP gather)와 많은 시추공에서 선택된 4개의 시추공 자료이다. $V_s$ 검층자료는 경험식을 통해 $V_p$ 검층자료로부터 계산되었다. 4개의 속도 물성 EEI 예측 %에러는 한 시추공에서의 포아송비를 제외하면 약 5%를 넘지 않는다. 그러나 속도로부터 직접적으로 계산되지 않는 공극률, 감마선 검층값, 밀도와 같은 물성들은 시추공에서의 EEI 역산 분석 결과가 만족스럽지 못하여 전체 자료에 EEI 역산을 적용할 수 없었다. 따라서 속도 물성의 경우 EEI 역산을 적용할 수 있지만 속도로부터 직접 계산이 되지 않는 물성의 경우는 EEI 역산 적용에 신중해야 할 것으로 판단된다.

${\ell}^1/{\ell}^2$ norm IRLS 방법을 사용한 강인한 탄성파자료역산 (Robust inversion of seismic data using ${\ell}^1/{\ell}^2$ norm IRLS method)

  • 지준
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2005년도 공동학술대회 논문집
    • /
    • pp.227-232
    • /
    • 2005
  • 탄성파 역산에 있어서 최소자승(${\ell}^2-norm$)해는 큰 오차에 매우 민감하게 반응하는 경향이 있다. 이에 반해서 ${\ell}^p-norm$ ($1{\le}p<2$)을 최소화하는 해는 잡음에 강인한 해를 보이나 보통은 좀 더 많은 계산이 요구된다. 반복적가중의 최소자승법(Iteratively reweighted least squares [IRLS] method)은 이러한 ${\ell}^p-norm$ 문제의 근사해를 효율적으로 구할 수 있도록 해준다. 본 논문에서는 작은 크기의 잔여분은 ${\ell}^2-norm$으로 큰 크기의 잔여분은 ${\ell}^2-norm$으로 적용되는 하이브리드 ${\ell}^1/{\ell}^2$최소화를 IRLS 방법에 쉽게 적용하는 기법을 소개한다. 모의 자료와 실제 현장자료에의 적용결과 큰 잡음이 포함된 경우 최소자승해보다 하이브리드 방법의 경우에 개선된 결과를 보임을 확인할 수 있었다.

  • PDF

한반도 남부의 지진파 토모그라피 연구 (Application of Seismic Tomography to the Region in and Near Southern Korean Peninsula)

  • 강익범;박정호
    • 자원환경지질
    • /
    • 제39권4호
    • /
    • pp.507-524
    • /
    • 2006
  • 한반도 남부 및 인근 지역에 대하여 3차원 토모그래피 연구를 수행하였다. 먼저 1차 연구지역은 지진계의 분포가 비교적 조밀하게 설치되어 있는 한반도 남동부지역으로 연구결과에 의하면 경상분지 지역의 지표에서는 퇴적층 지역에서는 상대적으로 저속도가 나타나며 7.5km 깊이에서는 상대적 고속도 분포가 넓게 나타났다. 이는 경상분지 내에서 화성기원의 암석들이 넓게 분포하는 것으로 보인다. 양산단층대 하부 약 $20{\sim}30km$ 깊이에서는 저속도 분포가 나타나는데 이는 양산단층대의 발달과 관계가 있을 것으로 여겨진다. 한반도 및 일본에서 관측된 자료를 이용하여 한반도 및 동해지역의 연구지역에서의 연구결과에 의하면 일본 동부에서 시작된 섭입대는 동해지역에까지 연장되는 것으로 분석되었다.

상시미동 표면파 분석에 의한 김해평야 퇴적층 심부 전단파 속도 결정 (Characterization of Deep Shear Wave Velocity Profiles in the Gimhae Plains Using the Microtremor Array Method)

  • 김재휘;정석호
    • 한국지반공학회논문집
    • /
    • 제38권8호
    • /
    • pp.17-27
    • /
    • 2022
  • 김해평야 퇴적층의 동적 특성을 결정하기 위해 6개의 위치에서 상시미동 수평-수직비(HVSR)법으로 지반 고유진동수(f0)를 도출하고 High frequency f-k방법과 Modified Spatial Autocorrelation 방법을 적용하여 상시미동 기반 표면파 분석을 수행하였다. 레일리파 분산곡선과 지반 고유진동수를 바탕으로 역해석을 실시하여 전단파 속도 주상도를 도출하였다. 분석 결과에 의하면 김해평야 퇴적층 지역은 1hz 이상의 진동수 대역에서 상당한 규모의 지반운동 증폭이 예상된다. 천부 퇴적층은 대체로 200m/s 내외의 전단파 속도를 보이고 60m~100m 사이의 깊이에서 기반암이 존재할 것으로 추정되며, 다수의 위치에서 천부 퇴적층과 기반암 사이에 Vs=400m/s 가량의 하성 퇴적층이 존재하는 것으로 나타났다. 역해석 결과를 바탕으로 김해평야 퇴적층의 전단파 속도-깊이 모델을 도출하였으며, 본 연구에 따르면 김해평야 퇴적층 지역에서 내진설계일반(KDS 17 10 00)의 S6지반에 해당하는 위치의 면적은 상당할 것으로 예상된다.

GRM과 RCS 방법을 이용한 굴절파 정적 시간차를 구하는 간단한 방법 (A simple approach to refraction statics with the Generalized Reciprocal Method and the Refraction Convolution Section)

  • Palmer Derecke;Jones Leonie
    • 지구물리와물리탐사
    • /
    • 제8권1호
    • /
    • pp.18-25
    • /
    • 2005
  • 풍화대의 심도가 급변하는 경암지형에서 얻은 자료로부터 굴절파의 정적 시간차(refraction statics)들이 계산되었다. 정보정값은 수신기 간격들의 12 배 정도되는 거리에 대해 보통 10 ms 이하에서 70 ms 이상의 값을 가진다. 이 논문에서는 자동 잔여 시간차가 항상 중요하지만은 않은 경암 지형에서 풍화대의 정확한 초기 굴절법 모델을 얻는 것이 얼마나 중요한 가에 대한 한 예를 보여주고 있다. GRM과 RCS(Refraction Convolution Section)법을 이용해 구한 간단한 풍화대 모델의 시간차 값들과 CG 법(Taner et al., 1998)을 이용한 최소 평균 제곱 역산에 의해 좀 복잡한 풍화대 모델에 대해 구해진 값들도 정확도 면에서 비교될 만하다. GRM 모델과 Taner 모델은 8.8 km 거리에 대해 체계적으로 평균 2 내지 4 ms 의 차이가 났다. 이들 두 모델들과 자동 잔여값을 포함한 최종 시간차 사이의 차이는 일반적으로 5 ms 이내이다. GRM 모델에서의 잔여값들은 때때로 Taner et al.(1998)의 모델의 잔여값들보다 작다. RCS법에서의 시간차들은 대략 10 ms 뒤에 발췌되지만 상대적인 정확도는 GRM 시간차들에 견줄만하다. 잔여 시간차값들은 굴절파 시간차값들과 일반적으로 상관관계를 보이며, 풍화대에서 더 낮은 탄성파 평균 속도를 이용함으로써 그 크기를 줄일 수 있다. 이들 결과들을 통해 풍화대에 적용된 부정확한 탄성파 평균 속도들은 주시로부터의 평균지연시간을 결정하는 역산 알고리듬의 어떤 문제점들보다도, 짧은 파장의 시간차의 원인이 될 수 있음을 알 수 있다.