• 제목/요약/키워드: 역산법

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

한국의 전자탐사 현황 (Electromagnetic Survey in Korea)

  • 조동행
    • 자원환경지질
    • /
    • 제39권4호
    • /
    • pp.427-440
    • /
    • 2006
  • 다른 물리탐사법과는 달리, 전자탐사의 실질적인 국내 활용과 연구는 1980년대부터이다. 이렇게 늦게나마 시작하게 된 것은 1970년대 미국 유학한 네 명의 연구자들의 영향이고, 이 들과의 공동연구와 유대는 현재까지 지속되고 있다. 그 동안 국내 연구진이 이룬 성과는 괄목할 만 하다. 이제 우리는 MT와 CSAMT탐사, 지자기수직탐사, 소형루프탐사, VLF탐사, GPR탐사, 시간영역 전자탐사 등 전 분야에 걸쳐 국제 수준의 기술력을 확보하게 되었고 신기술인 고주파수 전자탐사, 시간영역 전자탐사, 역산기법 분야에서 창의 있는 연구 성과를 내고 있다. 활용 분야도 지하자원탐사에서 지반조사, 환경조사, 고고학조사, 구조물 결함조사, 땅굴조사 등으로 엄청나게 확장되었다. 아직도 자료해석이 난삽하기 이를 데 없는 전자탐사법을 한 단계 높일 결정적인 돌파구가 국내 연구진에 의해서 나올 날이 멀지 않을 것이다.

상시미동 측정을 통한 천부 횡파속도 연구 (Shallow Shear-wave Velocities Using the Microtremor Survey Method)

  • 황윤구;김기영
    • 지질공학
    • /
    • 제16권4호
    • /
    • pp.381-392
    • /
    • 2006
  • 상시미동을 이용하는 수동적 표면파 탐사를 결정질암 기반암 지역에서 실시하여, 내진설계 지반분류 기준인 30 m 깊이까지의 평균 횡파속도(Vs30)를 구하고자 연구를 수행하였다. 시험자료는 조사지 A와 B에서 삼각형 배열과 'L' 자형 직각 배열로 2 ms 간격으로 4초 동안 기록하였다. 조사지 A에서 기록한 자료를 공간 자기상관 방법으로 분석하여 위상속도 스펙트럼을 작성하고, 수신기 배열의 길이 및 형태, 기록 횟수, 수신기 수 등의 주요인자들의 영향을 고찰하였다. 조사지 B에서 상시미동 기록자료의 역산을 통해서 구한 횡파속도와 근처 시추공에서 수직탄성파 방법으로 구한 결과가 대체로 일치함을 보여 상시미동을 이용한 표면파 탐사법이 Vs30을 구하는데 효율적인 방법으로 판단된다.

1차원 이외의 배경 전기전도도 구조에서 근사 적분방정식을 이용한 3차원 전자탐사 모델링 (3-D EM Modeling Using Approximate Integral Equation Method for the Models with Non 1-D Background Conductivity)

  • 이성곤
    • 지구물리와물리탐사
    • /
    • 제8권3호
    • /
    • pp.207-217
    • /
    • 2005
  • 이 연구에서는 배경 전기전도도 구조가 1차원 층서구조가 아닌 경우에도 적분방정식법을 이용하여 비교적 정확히 모델링 할 수 있는 근사적 방법에 대한 유도 과정을 제시하였으며, 이를 기반으로 전자탐사 수치 모델링을 구현하여 그 정확도 및 효율성에 대해 고찰하였다. 이 연구에서 제시된 방법은 통상적인 적분방정식법과 마찬가지로 1차원 층서구조에서의 Green함수를 그대로 이용할 수 있어 계산 측면에서 매우 효율적이며, 대체적으로 적분방정식을 구성하기 위한 행렬의 크기를 줄일 수 있다. 본 연구에서 개발된 근사 적분방정식 알고리듬을 간단한 수치 모델에 대하여 적용한 결과 정확도는 반응을 구하고자 하는 이상체 부근에서는 거의 차이가 없었으며, 또한 계산 시간에 대한 분석을 통해 이 모델링 알고리듬을 이용하여 역산을 수행하고자 할 때 계산 시간을 대폭 감소시킬 수 있음을 확인하였다. 본 연구에서 제안된 방법은 1차원 층서구조가 아닌 이미 알려진 공통된 구조론 갖는 여러 지질 모델에 대한 전자탐사 반응을 매우 효율적으로 구할 수 있을 것으로 판단된다.

해양 다성분 탄성파 자료를 이용한 가스하이드레이트 유망지역의 BSR 상하부 S파 속도 도출 (S-wave Velocity Derivation Near the BSR Depth of the Gas-hydrate Prospect Area Using Marine Multi-component Seismic Data)

  • 김병엽;변중무
    • 자원환경지질
    • /
    • 제44권3호
    • /
    • pp.229-238
    • /
    • 2011
  • 가스하이드레이트 부존량 평가에 있어서 해당 부존 지역의 S파 속도 정보는 암상과 공극유체의 정보를 파악하는데 결정적인 역할을 한다. 만일 퇴적층 내에 가스하이드레이트가 존재한다면 이 층에서의 P파 속도와 S파 속도는 동시에 증가하게 되며, 그 하부에 자유가스가 존재하는 경우 P파의 속도는 감소한다. 하지만 S파의 경우 공극을 채우고 있는 유체의 영향을 받지 않고 순수하게 매질을 통해서 진행하므로 하이드레이트 층의 하부에 자유가스층이 존재한다고 해도 그 속도가 변하지 않거나 오히려 매질의 영향으로 그 속도가 증가한다. 본 연구에서는 이러한 특성을 확인하기 위해 울릉분지의 가스하이드레이트 유망지역 중 탄성파 단면상에서 BSR(해저변 모방 반사면)이 강하게 분포하는 한 지점에서 한국지질자원연구원이 2009년 5월에 OBS(해저면 탄성파 기록계)를 이용하여 취득한 해저면 다성분 탄성파 자료를 이용하여 가스하이드레이트 부존 심도 부근의 P파 빛 모드전환 S파의 속도를 구하였다.OBS의 하이드로폰(hydrophone) 성분에 기록된 P파 자료를 이용하여 탄성파 주시 역산법을 수행하여 P파 속도 및 섬도 구조를 도출하였다. 해당지역에 취득한 2차원 반사법 탐사 자료는 기본 전산처리를 통해 구한 탐사지역의 기본 층서모델을 초기모델로 삼았다. 여기에 수평 2성분 지오폰(geophone)에 기록된 자료의 극성 분석을 통해 S파의 에너지가 최대로 모인 radial 성분 단면도를 생성하고 여기서 발췌한 주요 S파 이벤트의 주시를 이용해 포아송 비 정모델링을 수행하여 OBS가 위치한 지점에서의 포아송 비와 S파 속도구조를 최종적으로 도출하였다. 본 연구를 통해 탐사지역의 가스하이드레이트 존재로 인한 BSR 상하부 층의 P파 속도 역전 현상과 P파와는 달리 BSR 상부에서 히부로 갈수록 S파의 속도가 약간 증가하는 경향을 보여 결과적으로 자유가스층의 존재로 인한 BSR 하부에서 포아송 비 감소현상이 뚜렷함을 확인하였다.

천부 굴절법 탄성파 탐사 자료의 정보정 (Statics corrections for shallow seismic refraction data)

  • Palmer Derecke;Nikrouz Ramin;Spyrou Andreur
    • 지구물리와물리탐사
    • /
    • 제8권1호
    • /
    • pp.7-17
    • /
    • 2005
  • 천부 탄성파 굴절법 탐사를 이용하여 굴절이 발생하는 지층의 속도를 산출하는 것은 ill-posed 문제이다. 계산된 시간 변수들에서의 작은 변화들이 이로부터 산출된 속도들에 커다란 수평적 변화를 가져올 수 있으며 이는 종종 역산 알고리듬의 인위적인 오차를 유발한다. 이러한 인위적인 오차들은 모델링을 통해 인지되거나 보정되지 않는다. 그러므로 만약 모델에 근거한 역산을 통해 정밀한 지하 굴절 모델을 얻고자 한다면 정확한 초기 모델이 필요하다. 탄성파속도에서 인위적인 오차의 원인은 일반적으로 불규칙한 굴절면에 있다. 대부분의 경우에 GRM 방법을 이용하면 불규칙한 굴절면을 다룰 수 있고 굴절면의 정밀한 초기 모델을 만들 수 있다. 하지만 지표에 매우 가까운 극천부 지역 또한 불규칙하다면 GRM 방법의 효능은 감소하고 풍화대 보정이 필요하다. 천부 불균질대에 대한 일반적인 보정방법들은 수평적 확장이 제한된 극천부지역의 불균질대의 경우 효과적이지 못하다. 이럴 경우 GRM 평활화 통계적 방법(Smoothing Statics Method; SSM)이 지층의 속도를 좀 더 정확하게 평가할 수 있는 간단하고 실용적인 방법이다. GRM SSM 방법은 제로 XY 값을 가지고 계산된 시간-심도값들로부터 실제 XY 값을 가지고 얻어진 시간-심도값들의 평균값을 빼줌으로써 평활화 정보정을 수행한다. 심도가 깊어질수록 시간-심도값들이 XY 값에 따라 크게 변하지 않으므로 이들의 평균값은 최적값과 훨씬 더 같아진다. 그러나 극천부의 불균질대에 대해 시간-심도값들은 XY 값들이 증가함에 따라 수평적으로 이동하고 평균화 과정을 통해 대폭 감소한다. 결과적으로, XY값들에 대해 평균화된 시간-심도단면도는 천부의 불균질대에 대한 보정에 효과적이다. 또한 제로 XY 값을 가지고 계산된 시간-심도값들은 천부 불균질대의 영향과 대상 굴절면에 대한 시간-심도값들의 합으로 주어지므로 그들의 차는 정보정을 위해 주시로부터 빼주어야 할 대략적인 값들을 제공한다. GRM SSM 방법은 결정론적인 풍화대에 대한 보정법이라기 보다는 평활화 과정이다. 이 방법은 수평적으로 확장이 매우 제한된 천부 불균질대에 대해 가장 효과적이다. 모델과 현장 적용 결과들을 통해 GRM SSM 방법을 이용하여 불규칙한 굴절면을 가진 지층들에 대해 좀 더 신뢰할 수 있는 정밀한 탄성파 속도를 산출할 수 있음을 보여주고 있다.

유한차분법에 기초한 인공송신원 해양전자탐사 모델링 (Three-dimensional Modeling of Marine Controlled-source Electromagnetic Surveys Based on Finite Difference Method)

  • 한누리;남명진;구본진;김희준
    • 지구물리와물리탐사
    • /
    • 제15권2호
    • /
    • pp.66-74
    • /
    • 2012
  • 이 연구에서는 다양한 환경에서 인공송신원 해양전자탐사 자료의 정확한 해석을 위해 3차원 모델링 알고리듬을 개발하고, 복잡한 탄화수소자원 개발 환경에서 인공송신원 해양전자탐사의 적용 가능성을 파악하고자 암염에 인접한 저류층 모형을 구성하고 그 반응을 계산 및 분석하였다. 모델링 알고리듬은 유한차분법에 기초하고 있으며 이차 전기장 정식화를 도입함으로써 빠르고 정확한 반응 계산이 가능하도록 하였다. 이 알고리듬의 정확성은 2층구조 모형에 대해 해석해 및 기존의 3차원 모델링 알고리듬의 결과와의 비교를 통해 검증되었다. 암염에 인접한 저류층 모형에서는 암염과 저류층 및 배경 매질의 경계에 존재하는 세 군데 송신원에 대한 전기장 및 자기장 반응을 파악하고, 배경매질의 반응으로 정규화한 전기장 및 자기장 반응도 분석하였다. 그 결과 세 군데 송신원에 대한 반응을 종합하면 정규화된 전자기장 반응으로부터 지하의 고비저항 이상체의 개수와 각각의 위치를 파악할 수 있었으며 탄성파탐사를 통해 암염의 구조가 사전 정보로 알려졌다면 저류층의 존재를 밝힐 수 있었다. 이 연구에서 개발한 알고리듬은 향후 다양한 환경에서 인공송신원 해양전자탐사 자료의 반응 분석은 물론 역산 알고리듬 개발에 유용할 것으로 기대된다.

청천댐 주변의 천부 P파 및 S파 속도구조 (Near-surface P- and S-wave Velocity Structures in the Vicinity of the Cheongcheon Dam)

  • 박영환;김기영
    • 지구물리와물리탐사
    • /
    • 제16권3호
    • /
    • pp.109-118
    • /
    • 2013
  • 높이 23 m인 보령시 청천댐 상부와 인근에서 5 kg 해머를 이용하여 소규모 굴절법 및 표면파 탐사를 실시하였다. 인공 지진파의 수직 및 수평성분을 초동주시 토모그래피 및 레일리파 분산곡선 역산을 통하여 천부 P파속도(${\nu}_p$)와 S파속도(${\nu}_s$) 구조를 파악하였다. 중생대 퇴적암질 기반암의 평균 ${\nu}_p$${\nu}_s$는 댐마루 30 m 깊이에서 각각 1650 m/s와 950 m/s, 하류쪽 댐체 끝 지점 10 m 깊이에서 각각 1650 m/s와 940 m/s로 분석된다. 이 층들의 동포아송비는 0.24 ~ 0.25의 범위로, 고화된 퇴적층의 값과 일치한다. 댐체로부터 152 m 하류지점의 깊이 45 m 시추공 부근에서의 SH파 굴절법 토모그램은 10 ~ 12 m 깊이에 평균 vs가 870 m/s인 임계굴절면이 존재함을 보여준다. 이 지역에서는 덮개층의 ${\nu}_p$${\nu}_s$가 각각 500 m/s와 200 m/s인 상부층과 깊이에 따라 속도가 거의 선형으로 증가하는 하부층으로 구성되어 있다.

경관도로 등 신개념의 도로사업 개발에 관한 연구 (A Strategy for Developing New Road Projects)

  • 김응철
    • 한국도로학회논문집
    • /
    • 제9권2호
    • /
    • pp.115-127
    • /
    • 2007
  • 최근 도로선진국이라 할 수 있는 나라들에서는 다양한 도로사업의 발굴 및 장려정책이 도입되고 있다. 본 논문에서는 새로운 도로사업의 발굴을 위한 대안모형을 제시하였다. 새로운 도로사업 발굴을 위한 절차로는 근거법령의 모색 신개념 도로사업의 발굴, 평가 모형구축 및 운영지침 개발, 관련행정체계 개선 등 크게 4가지의 주요 절차가 필요한 것으로 제시되었다. 근거법령의 모색은 도로법 및 도로법 시행령, 시행규칙 등이 우선 고려되어야 할 것으로 분석되었으며 신개념의 도로사업 발굴절차는 크게 국가적 차원의 유일성 확보 도로, 다양한 가치기준에 의해 선정되는 도로 및 단일시설, 기술(설계 우수VE 및 LCCA사례 포함) 및 구조물 형태로서 발굴될 수 있는 사례로 세분류하였다. 평가모형 구축 및 운영지침 개발절차에서는 미국 NSBP프로그램의 절차 혹은 지속가능한 도시대상의 절차를 사례로 할 수 있으며 분석적 계층화법(AHP) 및 다지역산업연관모형(MRIO)의 적용가능성도 제시하였다. 또한, 관련 행정체계의 개선에서는 건설교통부의 도로기획관실 업무내역 분석을 통해 도로정책팀과 도로환경팀의 역할을 제시하였으며 평가위원회의 구성, 평가과정과 일체화된 예산지원절차를 추진방안으로 제시하였다.

  • PDF

유한 요소법을 이용한 2차원 TSP 모델링 (2 Dimensional TSP Modeling Using Finite Element Method)

  • 이홍;서정희;신창수
    • 지구물리와물리탐사
    • /
    • 제6권1호
    • /
    • pp.13-22
    • /
    • 2003
  • TSP탐사는 터널 및 지하공간의 건설 시 굴착 진행 중인 터널 주위의 지질구조를 파악하기 위하여 VSP탐사를 터널 내에 적용한 물리탐사 기법이다. 하지만 터널 굴착예정 지역의 지질구조에 따른 터널 주위에서 파의 전파특성을 보여줄 수 있는 인공 탄성파합성 알고리듬이 부족한 실정이다. 이에 본 연구에서는 기존의 2차원 유한 요소 탄성파 모델링 알고리듬을 발전시켜 터널 전방에 여러 형태의 지질구조가 존재하는 경우에 대하여 2차원 시간영역 탄성파 모델링을 수행하였다. 그리고 이로부터 얻은 인공합성 탄성파단면과 파면단면을 통하여 파의 전파특성을 규명하고자 하였다. 우선 TSP탐사와 동일한 형태의 송수신 배열을 한 후, 균질 매질에 대한 모델링을 수행하여 인공합성 탄성파단면과 파면단면을 동시에 구해 직접파의 전파특성을 분석하였다. 그리고 터널 굴착방향과 수직인 반사면과 평행인 반사면이 존재하는 경우에 대해 각각 모델링을 수행하여 균질 매질의 경우와의 비교, 분석을 통해 반사파들의 전파특성을 고찰하였다. 또한 송신원을 수진기의 앞쪽과 뒤쪽에 배열한 두 가지 경우에 대하여 모델링을 수행하여, 두 방법의 특성과 장단점을 고찰하였다. 이와 같은 모델링을 통하여 TSP탐사자료의 해석, 역산 알고리듬의 개발 및 송수신 개발 등에 도움을 줄 수 있을 것으로 사료된다.

한반도 남동부에서 부지효과를 고려한 스펙트럼 감쇠상수 χ 추정 및 강지진동의 추계학적 모사 (Estimation of Spectrum Decay Parameter χ and Stochastic Prediction of Strong Ground Motions in Southeastern Korea)

  • 조남대;박창업
    • 한국지진공학회논문집
    • /
    • 제7권6호
    • /
    • pp.59-70
    • /
    • 2003
  • 한반도 남동부에서 관측된 지진기록을 이용하여 스펙트럼 감쇠상수 $\chi$와 응력강하상수를 계산하였다. 특히, 스펙트럼 감쇠상수 $\chi$는 관측소별로 부지고유의 특성을 반영하여 새로운 방법으로 계산하였다. 이 방법은 2단계의 계산과정으로 구성된다. 첫 번째 단계에서는 코다(coda)파 정규화법을 이용하여 광역적 상수인 비탄성감쇠 Q를 역산하고, 이로부터 스펙트럼의 거리종속적인 감쇠효과를 일으키는 $\chi$$_{q}$값을 계산한다. 두 번째 단계에서는 앞서 계산된 $\chi$$_{q}$값을 관측소별로 동일하게 적용하여 거리와 무관한 스펙트럼 감쇠효과인 부지고유의 $\chi$$_{s}$값을 추정한다. 관측소에 관계없이 평균적으로 계산된 $\chi$값의 진원거리(R)에 따른 경향은 0.016+0.000157R이다. 응력강하상수는 1999년 경주근방에서 세 차례 발생한 지진의 관측자료를 이용하여 조남대와 박창업 방법으로 계산하였으며, 그 결과는 92-bar이다. 이들 상수를 이용하여 추계학적 모사법으로 한반도 남동부에서 발생 가능한 강지진동의 최대지반운동과 주파수에 따른 특성을 간접적으로 추정하였다. 이러한 추계학적 모사의 결과를 바탕으로 진원거리에 따른 지진동 감쇠공식을 유도하였으며 기존의 국내연구와 비교, 분석하였다.