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

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

MT 법의 3차원 모델링 개관 (Review on the Three-Dimensional Magnetotelluric Modeling)

  • 김희준;남명진;송윤호;서정희
    • 지구물리와물리탐사
    • /
    • 제7권2호
    • /
    • pp.148-154
    • /
    • 2004
  • 자기지전류(MT)법의 3차원 모델링에 대해 소개한다. 3차원 MT 모델링은 MT 반응의 물리적 특성의 이해뿐만 아니라 지하의 3차원적 전기비저항 구조를 재구성하기 위한 역산법의 개발에도 필수적이다. 지난 20년 동안 3차원 모델링에 관한 여러 수치기법들이 개발되었으나 그 실용성에는 많은 한계가 있었다. 그러나 최근에는 컴퓨터의 급속한 발전과 대형 연립방정식에 대한 반복해법의 발전에 힘입어 이전에는 어려웠던 복잡한 3차원 구조에 대한 MT 반응을 효율적으로 모델링할 수 있게 되었다. 유한차분법에서는 자기 flux와 전류의 보존법칙을 만족하면서 전기장의 불연속을 표현할 수 있는 staggered 격자의 사용이 보편화되었다. 대형 연립방정식에 대한 수치해의 수렴성은 Krylov 부분공간법, 적당한 전처리 기술 및 정적 발산보정법을 채택함으로써 크게 향상된다. 변요소를 사용하는 벡터 유한요소법으로도 전기장의 불연속 문제를 해결할 수 있으며 이 방법이 가진 기하학적 유연성은 불규칙한 지표기복을 포함한 복잡한 구조를 모델화할 때 특히 유용하다.

3차원 불균질 횡등방성 매질에 대한 탄성파 초동 주시 모델링 (Seismic First Arrival Time Computation in 3D Inhomogeneous Tilted Transversely Isotropic Media)

  • 정창호;서정희
    • 지구물리와물리탐사
    • /
    • 제9권3호
    • /
    • pp.241-249
    • /
    • 2006
  • 한국의 지질 환경은 암석 분포가 매우 다양하고 복잡한 구조 활동의 영향을 받아 지하매질의 이방성 특성이 국부적으로 심하게 변화한다. 기존의 이방성 주시 모델링의 경우 지질 모델을 2차원으로 단순화시킴으로써 이러한 복잡한 지질 환경을 제대로 고려할 수 없었다. 또한 약 이방성 가정을 사용하여 실제로 나타날 수 있는 지하 매질의 심각한 이방성 영향을 주시 모델링에서 고려할 수 없었다. 이에 이 연구에서는 보다 실제적이고 복잡한 3차원 횡등방성 매질(transversely isotropic media)에서 q-P파의 초동 주시 양상을 모사할 수 있는 주시 모델링 알고리듬을 개발하였다. 이 알고리듬에서는 2차원 비선형 주시 내삽(2D nonlinear traveltime interpolation) 기법과 주시의 3차원 격자 채움법(mapping)을 이용한 직접 전파법(direct calculation)을 통해 급격한 물성의 변화에도 주시 계산이 가능하도록 하였다. 또한, 최소 주시 계산과정에서 수치 미분을 통한 최대 경사법(steepest descent method)을 사용하여 약 이방성 가정을 극복하였다. 개발된 알고리듬은 해석해와 비교하여 그 타당성을 검증하였고 3차원 2층구조에 대한 주시 계산을 수행하여 물성이 급격히 변화하는 모델에 대해서도 안정적으로 주시 계산이 이루어짐을 확인하였다. 이 연구에서 개발한 3차원 주시 모델링 알고리듬은 향후 구조보정이나 토모그래피 알고리듬 개발에 사용될 수 있을 것으로 기대한다.

HWAW 기법의 비수평 출상구조지반 적용에 대한 고찰 (A Study on Application of HWAW Method to the Non-horizontally Layered Soil Structure)

  • 방은석;박헌준;박형춘;김동수
    • 한국지반공학회논문집
    • /
    • 제25권4호
    • /
    • pp.5-17
    • /
    • 2009
  • HWAW 기법은 시간-주파수 해석을 통해 신호/잡음비가 최대가 되는 부분을 이용하여 실험 분산곡선을 획득하고 full wavefield를 고려하는 정모델링 기법을 통해 역산을 수행한다. 따라서 다른 표면파 기법에 비해 짧은 측선 설정이 가능하며 2차원 전단파 속도 분포를 도출하기에 유리하다. 수평층 가정이 성립되지 않는 지반에서의 HWAW 기법 적용성을 평가하기 위해서 수치해석 검증 연구를 수행하였다. HWAW 기법을 통해 각 거 리별로 획득한 실험 분산곡선들은 전반적으로 full wavefield를 고려한 이론 분산곡선과 유사한 형태를 보이는 것을 확인할 수 있었다. 경사면 반사파에 의한 분산곡선의 왜곡현상은 적절한 감지기 간격 조절을 통해 줄여줄 수 있음을 확인할 수 있었다. 도출된 실험 분산곡선은 가진 원으로부터 감지기까지의 영역보다는 감지기 사이 영역에 주되게 영향을 받는 것으로 나타났다. 이로써 HWAW 기법으로 합리적인 전단파 속도 분포의 도출이 가능한 것으로 판단되었으며 고찰한 결과를 바탕으로 실제 현장에서 지반 영상화를 시도하였다.

전기비저항 토모그래피에 의한 지하구조의 3차원 영상화 (Three-dimensional Imaging of Subsurface Structures by Resistivity Tomography)

  • 이명종;김정호;정승환;서정희
    • 지구물리와물리탐사
    • /
    • 제5권4호
    • /
    • pp.236-249
    • /
    • 2002
  • 본 연구에서는 종래의 지표 전기비저항 탐사자료의 3차원 영상화 알고리듬을 확장하여 여러 시추공을 이용한 3차원 전기비저항 토모그래피 영상화 알고리듬을 개발하였다. 개발된 알고리듬은 유한요소법을 이용한 3차원 전기비저항 모델링과 최소자승 역산 알고리듬으로 구성되며, 역산의 분해능을 높이기 위하여 ACB(Active Constraint Balancing)법을 채용하였다. 수치모형실험 및 분해능 분석을 통하여 3차원 전기비저항 토모그래피가 효용성이 높은 탐사법임을 알 수 있었으며, 지하구조에 대한 고분해능의 입체 영상을 획득할 수 있는 방법임을 입증하였다. 또한, 전기비저항 토모그래피에서 지형효과를 정확히 반영함으로써, 영상의 왜곡을 배제하여 더욱 정확한 지하구조 영상을 획득할 수 있었다. 현장자료에 대한 적용성 검토를 위하여 화강암 석산지역에서 획득한 3차원 탐사자료에 3차원 전기비저항 토모그래피를 수행한 결과 신뢰도 높은 3차원 전기비저항 영상을 획득할 수 있었다.

비저항 측선 근처 철케이싱 시추공 한개에 의한 $-{\rho}a$ ($-{\rho}a$ by One Steel Casing Borehole near Resistivity Survey Line)

  • 정현기
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2006년도 공동학술대회 논문집
    • /
    • pp.83-86
    • /
    • 2006
  • 쌍극자 전기탐사현장에서 드물지않은 음의 겉보기비저항에 대해 선진 장비회사의 메뉴얼들에서 언급되되 정확한 현상이 규명 안되었으나, 국내에서 계측 측면 및 가능한 지질모델 제시와 함께 전산모델링 결과와 $-{\rho}a$ 포함유무 3차원 역산결과가 일부 보고된 바 있다. 전도성 탄층, 전봇대 접지, 신호등, 철케이싱 시추공 등과 관련한 현장문제가 계속되어 해석시 딜레마 극복 지침이 필요하다. 여기서는 단순 철케이싱 한개의 효과 전산 모델링 결과 및 일부 측선 설계요령과 해석시 유의사항을 고찰한다. $-{\rho}a$는 분명한 지하정보이지만 잘못 해석시 진실과 먼 결과 오류 가능성이 높으며, 향후 역산 알고리즘의 유용한 발전이 절실하다.

  • PDF

근거리장에서 HWAW 기법의 수치해석적 검증 (Numerical Verification of HWAW Method in the Near Field)

  • 방은석;박형춘;김동수
    • 한국지반공학회논문집
    • /
    • 제23권2호
    • /
    • pp.5-17
    • /
    • 2007
  • 기존의 표면파 기법들은 근거리장 효과를 피하기 위하여 적절한 실험 방법 및 필터 기준을 제시하고 있다. HWAW 기법은 다른 표면파와는 다르게 이러한 근거리장 영역에서의 파의 거동을 적극적으로 활용하여 지반의 전단파 속도 주상도를 도출한다. 에너지가 작은 저주파 신호를 시간-주파수 해석을 통해 획득하며 체적파의 영향을 고려한 단일 배열 역산을 통해 가능하다. 이러한 사항을 검증하기 위해서 근거리장에서의 탄성파의 전파 형상을 수치해석을 통하여 모사하였으며 다양한 지반 형태를 대표하기 위해 5개의 지반모델을 구성하였다. 수치모델링을 통해 획득한 시간영역 신호를 이용하여 HWAW 기법으로 실험분산곡선을 도출하였다. 도출된 실험 분산곡선과 3차원 정모델링(Kausel의 방법)을 통해 계산된 이론 분산곡선의 형상이 각 위치에서 저의 일치하였다. 그러므로 HWAW 기법을 수행하는데 있어 근거리장 영역이 포함된 실험 분산곡선을 획득하여도 Kausel의 방법을 이용하여 역산을 수행할 경우 신뢰성 있는 지반의 전단파 속도를 도출하는 것이 가능하다는 것을 확인할 수 있었다.

국내 금속광 탐사를 위한 쌍극자-쌍극자 전기탐사의 적용성 연구 (A Study of Feasibility of Dipole-dipole Electric Method to Metallic Ore-deposit Exploration in Korea)

  • 민동주;정현기;박삼규;전효택;곽나은
    • 지구물리와물리탐사
    • /
    • 제11권3호
    • /
    • pp.250-262
    • /
    • 2008
  • 금속광산탐사에서 많이 이용되는 물리탐사 방법 중의 하나인 쌍극자-쌍극자 배열 전기비저항탐사의 실제 광산탐사에 대한 적용성을 평가하기 위하여 단순화된 광상구조에 대하여 현장자료 시뮬레이션 및 역산해석을 수행하였다. 실제 광상모형으로는 국내 여러 금속광상 중에서 가장 분포 빈도가 높은 맥상광상을 대상으로 하였다. 국내 맥상광상의 맥폭이 대부분 수십 cm에서 2m이내라는 점에 기초하여 맥폭을 변화시키고, 광맥의 물성을 변화시켰으며, 쌍극자-쌍극자 배열 전기비저항탐사의 전극간격은 40m로 가정하였다. 맥폭이 1m 이하일 경우 주변암과 광맥의 전기적 물성차가 300배 이상 되는 경우에도 저비저항대가 나타나지 않았다. 맥폭이 커질수록 그리고 비저항의 차가 클수록 저비저항대가 뚜렷하게 나타났으나, 넓은 전극간격과 역산에 이용된 셀 크기가 심부로 갈수록 커져서 대체로 저 비저항대의 폭이 과장되게 나타났고, 깊이 또한 실제 깊이보다 깊게 나타났다. 실제 천열수광상 모형을 가정하고 현재 가행 중인 광산에서 측정한 전기비저항을 근거로 하여 현장자료 시뮬레이션 및 역산을 수행한 결과 광맥의 전기비저항 값이 높은 경우에도 광맥에 의한 영향이 잘 나타나지 않는다는 것을 알 수 있었다. 이와 같이 구조지질학적으로 광맥에 수반되는 열수변질대가 전기비저항탐사 자료에서 이상대로 포착되므로 천열수 광산의 경우 이러한 주변 열수변질대를 탐사대상으로 선정하는 것이 바람직하다. 이상의 결과들을 종합해 볼 때 광맥이 지하 깊은 곳에 존재하며 광맥과 주변암과의 전기적 물성차가 크지 않은 경우, 일반적으로 수행되는 지표면 전기비저항탐사로는 광맥의 존재여부조차 파악할 수 없는 경우가 있으며, 저비저항대가 나타난다 하더라도 실제 구조와 다른 과장된 양상으로 나타난다는 것을 알 수 있었다. 이 연구에서 제시한 모델링 및 역산자료는 기본적이고 단순한 모형실험이지만 일선 현장 자료를 바탕으로 수행한 실험연구로서 향후 보다 다양한 모형에 대하여 모델링을 수행하고 표준 및 정량적 지침을 제시해 나간다면 현장 측선 설계나 탐사결과의 해석시에 유용한 지침으로 활용될 수 있을 것이다.

단일 UAV를 이용한 해안 지형 측량 (Survey of coastal topography using images from a single UAV)

  • 노효섭;김병욱;이민재;박용성;방기영;유호준
    • 한국수자원학회논문집
    • /
    • 제56권spc1호
    • /
    • pp.1027-1036
    • /
    • 2023
  • 해안 환경에서 정확한 지형 조사는 필수적이나 지점식 조사 기법이 일반적이며, 이마저도 육상과 해저면을 독립적으로 계측한다. 본 연구에서는 단일 UAV만을 이용해 육상 및 해저 지형을 측량하는 방법을 소개한다. 세부적으로, UAV 영상을 활용해 지형 및 수심을 계측하는 두 알고리즘을 각각 적용한 뒤 결과물을 정합하여 수행된다. 해빈 지형의 취득은 공간 스캔 영상을 이용하는 Structure-from-Motion Multi-View Stereo 기술이 적용된다. 해저 지형 측량을 위해서는 고정비행으로 취득된 시계열 파랑 영상을 이용하는 분산관계식 기반 수심 역산 기법이 적용된다. 두 요소기술로 산정한 수치 표고모형을 좌표에 따라 정합한 후, 쇄파대 및 포말대와 같이 두 요소기술 적용이 불가능한 부분을 내삽하여 최종적으로 연속된 근해역 지형을 취득할 수 있다. 본 UAV 기반 지형 측량 기법을 경상북도 포항시의 장사해수욕장에 적용한 결과 세부적인 지형적 특징을 재현해낼 수 있었다. 본 연구에서 제안되는 통합 모니터링 방법은 기존 방법들에 비해 시간, 비용, 안전 측면에서 이점이 있어 해안지역의 침퇴적 분석에 효과적으로 적용될 수 있을 것으로 기대된다.

MT 탐사 반응에서 정적효과: 적분방정식을 통한 고찰 (Static Effect in Magnetotelluric Responses: An Implication from the EM Integral Equation)

  • 송윤호
    • 지구물리와물리탐사
    • /
    • 제27권3호
    • /
    • pp.181-195
    • /
    • 2024
  • 이 해설 논문에서는 자기지전류(magnetotelluric; MT) 탐사에서의 정적효과가 지하 불균질체의 경계면에서 축적되는 전하 축적에 의해 발생하는 물리적 현상임을 설명하였다. 먼저 정적효과를 일으키는 물리적 현상에 대한 이해를 위해 정전유도와 경계면에서의 전하 축적의 차이에 대해 모식도를 이용해 설명하고 분석하였다. 그리고 적분방정식을 통해 지하에 있는 불균질체에 1차 전기장이 가해졌을 때, 전기전도도 경계면에서의 전하 축적에 의한 2차 전기장이 정적효과로 나타남을 명확히 하였다. 따라서 MT 탐사에서 1차원 층서구조나 2차원 TM (transverse magnetic) 모드의 경우를 제외하고 정적효과는 항상 존재하고 또한 그 용어와는 달리 '정적'이 아닌 주파수 의존적임을 보였다. 정적효과가 지하 불균질체에 의한 2차장임에도 MT 탐사의 역산에서 수식적으로 해결되지 않고 심부 구조의 해석에 장애 요인으로 작용하는 이유는 현장 탐사에서 피할 수 없는 '주파수 및 공간 영역에서의 과소 샘플링'이기 때문이다. 따라서 광대역 주파수와 전체 탐사 대상 면적에서의 연속 측정이 아닌 실제 MT 탐사를 고려할 때, 3차원 역산시에 정적효과를 전체 탐사 영역에서 정규분포를 갖는 변수로 가정하고 이를 역산의 제약 조건으로 포함하는 시도에 대해 알아보았다. 적분방정식 전자탐사 모델링의 이해를 돕기 위해 3차원 적분방정식의 유도, Green 텐서 및 산란전류의 수학적 분석을 부록에 자세히 서술하였다.

2차원 전기비저항 탐사를 위한 변형된 전극배열법에 관한 연구 (A Study on the Modified Electrode Arrays in Two-Dimensional Resistivity Survey)

  • 김정호;이명종;송윤호;정승환
    • 지구물리와물리탐사
    • /
    • 제4권3호
    • /
    • pp.59-69
    • /
    • 2001
  • 2차원 전기비저항 탐사에 일반적으로 사용되는 전극배열들의 단점을 보완하기 위해 단극과 쌍극자를 사용하는 기존의 전극배열로부터 5 종류의 변형된 전극배열법을 제안하였다. 변형된 단극 배열법은 원거리 접지전극이 필요한 단극 배열법의 비효율성을 보완하기 위한 전극배열법이다. 4 종류의 쌍극자를 이용하는 변형된 전극배열법은 쌍극자 또는 단극-쌍극자 배열이 갖는 문제점, 즉 낮은 측정값에 의한 낮은 신호대 잡음비를 개선하여 높은 측정값의 획득이 가능하도록 설계하였다. 수치 모델링과 2차원 역산을 이용한 수치 실험을 수행하였으며, 전기적 잡음에 의한 영향과 분해능에 대하여 고찰하여 본 논문에서 제안된 전극배열의 효용성을 입증하였다.

  • PDF