• 제목/요약/키워드: Korean Geophysical Society

검색결과 1,538건 처리시간 0.028초

그래픽 프로세서를 이용한 탄성파 수치모사의 계산속도 향상 (Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit)

  • Nakata, Norimitsu;Tsuji, Takeshi;Matsuoka, Toshifumi
    • 지구물리와물리탐사
    • /
    • 제14권1호
    • /
    • pp.98-104
    • /
    • 2011
  • 탐사 지구물리학에서 수치 모사는 지하매질에서의 탄성파 전파 현상을 이해하는데 중요한 통찰력을 제공한다. 탄성파 모사는 음향파 근사에 의한 수치 모사보다 계산시간이 많이 소요되지만 전단응력 성분을 포함하여 보다 현실적인 파동의 모사를 가능하게 한다. 그러므로 탄성파 모사는 탄성체의 반응을 탐사하는데 적합하다고 할 수 있다. 계산 시간이 길다는 단점을 극복하기 위해 본 논문에서는 그래픽 프로세서(GPU)를 이용하여 탄성파 수치 모사 시간을 단축하고자 하였다. GPU는 많은 수의 프로세서와 광대역 메모리를 갖고 있기 때문에 병렬화된 계산 아카텍쳐에서 사용할 수 있는 장점이 있다. 본 연구에서 사용한 GPU 하드웨어는 NVIDIA Tesla C1060으로 240개의 프로세서로 구성되어 있으며 102 GB/s의 메모리 대역폭을 갖고 있다. NVIDIA에서 개발된 병렬계산 아카텍쳐인 CUDA를 사용할 수 있음에도 불구하고 계산효율을 상당히 향상시키기 위해서는 GPU 장치의 여러 가지 다양한 메모리의 사용과 계산 순서를 최적화해야만 한다. 본 연구에서는 GPU 시스템에서 시간영역 유한차분법을 이용하여 2차원과 3차원 탄성과 전파를 수치 모사하였다. 파동전파 모사에 가장 널리 사용되는 유한차분법 중의 하나인 엇갈린 격자기법을 채택하였다. 엇갈린 격자법은 지구물리학 분야에서 수치 모델링을 위해 사용하기에 충분한 정확도를 갖고 있는 것으로 알려져 있다. 본 논문에서 제안한 모델링기법은 자료 접근 시간을 단축하기 위해 GPU 장치를 메모리 사용을 최적화하여 가능한 더 빠른 메모리를 사용한다. 이점이 GPU를 이용한 계산의 핵심 요소이다. 하나의 GPU 장치를 사용하고 메모리 사용을 최적화함으로써 단일 CPU를 이용할 경우보다 2차원 모사에서는 14배 이상, 3차원에서는 6배 이상 계산시간을 단축할 수 있었다. 세 개의 GPU를 사용한 경우에는 3차원 모사에서 계산효율을 10배 향상시킬 수 있었다.

베트남 북부 네안 희토류 산출지의 지구화학탐사 (Geochemical exploration for REE occurrence in Nghe An Area within Northern Vietnam)

  • 허철호;;이재호
    • 자원환경지질
    • /
    • 제45권6호
    • /
    • pp.599-622
    • /
    • 2012
  • 베트남 북부 네안(Nghe An)성 뀌 차우(Quy Chau)에서 남동쪽으로 약 14 km 떨어진 쵸 빈(Chau Binh) 일대를 대상으로 1차 토양지구화학탐사를 수행하였다. 시료채취 간격은 가로 300 m 간격으로 14 line, 세로 500 m 간격으로 15 line이었으며 토양시료 194개를 채취하였다. 그리고, 1차 토양지구화학탐사 결과를 바탕으로 총희토류(TREO) 함량이 높은 격자점을 대상으로 2차 정밀 피트탐사를 수행하여 총 7개 피트에서 토양 56개를 채취하였다. 조사지역의 지질은 반 찌응(Ban Chieng) 흑운모 화강암 복합체와 다이 록(Dai Loc) 편마암질 화강암 복합체가 편암, 편마암, 석회암으로 구성된 부캉(Bukang)층을 관입하고 있다. 본 역의 주요 광화작용은 주석, 루비, 그리고 화강암 복합체의 주변부 변질대(albitization)에 분포하고 있는 것으로 사료되는 함희토류 모나자이트(약 300 g/t)와 제노타임(약 10 g/t)의 산출이 특징적이다. 중사시료에서 확인한 모나자이트와 제노타임의 근원암을 규명하기 위해 토양지구화학탐사를 수행하였다. 1차 토양지구화학탐사에서 채취한 토양시료에 대하여 희토류원소(La, Ce, Nd, Y 등) 및 미량원소(Th, Sr, Ba, Nb 등)의 화학성분을 ICP-MS를 이용하여 분석한 결과, 배경토양의 총희토류 함량이 지각 평균의 약 2배를 초과하고 있으며, 중희토류(약 2배)와 경희토류(약 1.84배)가 모두 부화되어 있다. 2차 정밀 피트탐사에서 채취한 토양시료에 대하여 희토류원소(La, Ce, Nd, Y 등)의 화학성분을 ICP-MS를 이용하여 분석한 결과, 격자점 4-7번 피트에서 광체로 추정되는 노두를 인지하였고 본 광체는 N40W 방향으로 발달되어 있다. 1차 토양지구화학탐사 및 2차 정밀 피트탐사를 종합적으로 고찰한 결과, 네안성 뀌차우 일대에서는 토양에서 이상대를 보이며 광체의 연장방향으로 추정되는 희토류 부존 유망 지역을 잠정적으로 선정하였으며, 향후 광체로 추정되는 구간에 대하여 지구물리탐사 및 시험시추를 수행할 예정이다.

포항분지 해상 중소규모 CO2 저장 실증연구 (Demonstration-scale Offshore CO2 Storage Project in the Pohang Basin, Korea)

  • 권이균
    • 지질공학
    • /
    • 제28권2호
    • /
    • pp.133-160
    • /
    • 2018
  • 신기후체제 출범 이후 온실가스 감축목표 달성을 위한 대규모 $CO_2$ 직접 감축의 거의 유일한 방안으로 알려져 있는 대용량 $CO_2$ 지중저장의 중요성이 강조되고 있다. 포항 분지 해상 중소규모 $CO_2$ 지중저장 실증사업은 대규모 $CO_2$ 저장 사업에 필요한 핵심 기술을 우선적으로 확보하기 위해 추진하고 있는 연구개발 사업이다. 국내 저장 실증을 위한 $CO_2$ 저장소 확보를 위해 2010년부터 3년간 육상-연안 $CO_2$ 저장소 탐사 연구 과제를 수행하여, 한반도 동남부 육상 및 연안에 걸쳐 발달하고 있는 포항분지를 중소규모 $CO_2$ 지중저장 실증 후보지역으로 선정한 바 있다. 2013년부터는 포항분지 해상 지층을 대상으로 저장소 탐사, 탐사시추, 저장소 특성화, 저장 설계, 플랫폼 설계 및 설치, 주입공 시추 및 완결, 주입설비 설치, 주입 및 모니터링 등의 개별 기술 개발을 수행하였으며, 2017년 초 국내 기술력을 기반으로 해저 실증 저장소에 국내 최초로 $CO_2$ 시험주입 실증에 성공하였다. 향후에는 $CO_2$ 저장소 연속 운전 실증을 위한 수송체계 구축과 저장소 운영체계 확립 연구를 진행할 계획이다. 안전하고 효율적인 저장소 운영을 위해 상시 주입 모니터링 및 정기적인 거동 누출 환경 모니터링을 수행할 계획이며, 이를 위한 주입 전 모니터링을 완수하고 모니터링 체계를 구축하였다. 포항분지 해상 중소규모 $CO_2$ 지중저장 실증 연구를 통해 국내 대규모 CCS 통합실증 및 상용화 사업에 필요한 핵심적 $CO_2$ 저장 기술을 자립화 하였고, 성공적인 시험 주입으로 대규모 $CO_2$ 지중저장 사업의 기술적 경험적 기초를 확보하였다. 또한 본 연구 사업은 국제적으로 점차 확대되고 있는 $CO_2$ 지중저장 기술 시장에서 우리나라의 기술 경쟁력을 확보하는데 기여할 것으로 기대된다.

고해상 해저 탄성파 탐사자료에 대한 너울영향 보정 (Swell Effect Correction for the High-resolution Marine Seismic Data)

  • 이호영;구남형;김원식;김병엽;정순홍;김영준
    • 지구물리와물리탐사
    • /
    • 제16권4호
    • /
    • pp.240-249
    • /
    • 2013
  • 해양지질 조사 또는 엔지니어링 목적의 해양 탄성파 탐사에서 해상의 너울로 인하여 탐사자료의 품질이 저하된다. 1 ~ 2 m의 너울은 탐사 중에 종종 발생하는데, 1 m 이내의 정밀도를 요하는 고해상 해저탐사자료에서 이를 보정하여 줌으로써 탄성파 단면도의 수평적 연속성을 높일 수 있다. 이 연구에서는 8채널 고해상 에어건 탐사자료와 3.5 kHz 천부지층탐사자료에 대하여 너울영향 보정을 적용하였다. 너울영향을 효율적으로 보정하기 위해서는 탐사자료에 나타나는 해저면의 위치를 정확하게 산출하는 것이 중요하다. 이 연구에서는 해저면 수심값 추출을 위하여 해저면 부근의 최대 진폭을 이용하거나 최대진폭값의 일정 기준을 넘어서는 지점을 해저면 위치로 추출하는 방법 등을 사용하였다. 품질이 낮은 자료에 대해서는 수심값 추출이 용이하도록 엔벨로프 또는 해저면 요소파(wavelet) 신호와 상호상관을 수행한 자료를 사용하였다. 이와 같은 방법으로 산출한 수심값으로부터 평균값을 구한 후, 그 차이를 보정하여 주었다. 시험 적용된 에어건 탐사자료에서는 약 0.8 m, 2종의 3.5 kHz 천부지층탐사자료에서는 각각 약 0.5 m와 약 2.0 m 범위의 너울영향을 보정하였다. 자료의 상태에 따라 적절한 해저면 수심값 추출 방법을 사용하여 탐사자료의 너울 영향을 보정함으로써 지층의 연속성이 향상된 고품질의 고해상 해저 탄성파 단면도를 제작할 수 있었다.

호주 Gippsland Basin에서 $CO_2$ 주입 중 단층 재활성화의 가능성 (Fault reactivation potential during $CO_2$ injection in the Gippsland Basin, Australia)

  • Ruth, Peter J. van;Nelson, Emma J.;Hillis, Richard R.
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.50-59
    • /
    • 2006
  • 현재의 응력장내에서 단층 재활성화를 야기하는데 필요한 공극압의 증가를 추정함으로써 재활성화 위험도를 결정하는 FAST(단층 분석 확인 기술)를 이용해 Gippsland Basin의 단층 재환성의 위험도가 계산되었다. Gippsland Basin의 응력 형태는 주향이동단층과 역단층의 경계부근으로서 즉, 최대 수평 압력$({\sim}40.5\;MPa/km)$ > 수직 압력(21 MPa/km) ${\sim}$ 최소 수평 압력(200 MPa/km)이다. 공극압은 Golden Beach Subgroup의 Campanian volcanics 상부에서 정수압이다. 여기에서 결정된 NW-SE 최대 수평 응력 방향$(139^{\circ}N)$은 이전의 측정값들과 대체로 일치하고 Gippsland Basin에서의 NW-SE 최대 수평 응력 방향을 입증한다. Gippsland Basin의 단층 재활성화 위험도는, cohesionless fault$(C=0;\;{\mu}=0.65)$와 healed fault$(C=5.4;\;{\mu}=0.78)$, 두 가지 단층 강도 시나리오를 이용해서 계산되었다. 상대적으로 높고 낮은 재활성화 가능성을 가진 단층들의 방향은 cohesionless fault 와 healed fault 모두에 대해 거의 동일하다. NE-SW 주향방향의 큰 각을 가진 단층들은 현재의 응력상태하에서는 재활성화 가능성이 거의 없다. SSE-NNW 과 ENE-WSW 방향의 큰 각을 가진 단층들이 단층 재활성화 위험도가 가장 높다. 부가적으로 NE-SW 주향 방향의 작은 각을 가진 단층(thrust 단층)은 상대적으로 높은 재활성화 위험도를 가지고 있다. 최적 방향 단층들에 대한 가장 높은 재활성화 위험도는 cohesionless fault에 대해서는 추정 공극압의 3.8MPa$({\sim}548psi)$ 증가(Delta P), healed fault에 대해서는 15.6MPa 증가에 해당된다. 이 논문에서 제시된 단층 재활성화 분석으로부터 얻은 공극압 증가의 절대값은 지구역학적인 모델(원위치 응력과 암석 강도 자료)에서의 불확실성으로 인해 큰 오차를 수반한다. 특히, 최대 수평 응력 강도와 단층 강도 자료는 좁은 범위에 한정되어 있지 않다. 그러므로 단풍 재활성화 분석은 저류층 내에서 최대로 허용할 수 있는 공극압 증가를 직접 측정하는데 사용될 수 없다. 이러한 종류의 단층 재활성화 분석은 단지 단층 재활성화의 상대적인 위험도의 평가에 사용될 수 있을 뿐이고, 재활성화에 앞서 단층이 견딜 수 있는 공극압 증가의 최대 허용치를 결정하는데는 사용할 수 없다고 주장하고자 한다.

$CO_2$로 포화된 사암의 속도-유효응력 반응 (Velocity-effective stress response of $CO_2$-saturated sandstones)

  • Siggins, Anthony F.
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.60-66
    • /
    • 2006
  • 세 개의 서로 다른 사암 샘플들 -두 개의 합성 샘플과 한 개의 현장 샘플- 에 대해 현장 저류층의 대표적인 구속압력과 공극압력하에서 초음파 시험이 수행되었다. 세가지 사암 샘플들은 (a) 캘사이트 시멘트(calcite intergranula. cement (CIPS))로 만든 합성사암, (b) 실리카 시멘트(silica intergranular cement)로 만든 합성 사암 (c) Otway Basin 의 Boggy Creek 1 시추공에서 시도되는 $CO_2$ 파일럿 프로젝트의 대상 암석층 중 Waarre 층으로부터 추출한 코아 샘플로 구성되어 있다. 공극률은 각각 32%, 33%, 26%이다. 초기시험은 실내건조(room-dried) 상태에 있는 코아들에 대해 구속응력을 5 MPa 씩 단계별로 65 MPa 까지 증가시키며 이루어졌다. 그리고 나서 모든 코아들에 처음에는 온도 $22^{\circ}C$에서 6 MPa 공극압력으로 기체상의 $CO_2$를, 그 다음에는 온도 $22^{\circ}C$ 에서 7 MPa 부터 17 MPa 까지 5 MPa 씩 증가시키면서 액체상의 $CO_2$를 주입하였다. 구속응력은 10MPa부터 65 MPa까지 달리 하였다. P와 S 초음파 파형들이 유효응력이 증가할 때마다 기록되었다. 속도-유효응력 반응들이 P 파와 S 파에 대해 실험 자료들로부터 계산되었으며, 감쇠(1/Qp)들은 스펙트럼 비 방법을 이용하여 파형들로부터 계산되었다. 각각의 사암들에 대한 이론적인 속도-유효응력 계산은 $CO_2$ 압력-밀도 와 $CO_2$ 체적계수-압력 상 다이어그램(phasediagram), Gassmann 유효 매질 이론(effective medium theory)을 이용하여 구하였다. 기체상의 $CO_2$ 주입은 속도-유효응력에서 건조상태(공기로 포화된 상태)에 비해 거의 무시할만한 변화를 가져왔다. 다양한 공극압력에서 액체상의 $CO_2$ 주입은 공기로 포화된 상태에 비해 속도-유효응력 반응을 평균 약 8% 정도 낮게 한다. 실험자료들은 높은 유효응력에서 Gassmann 계산들과 일치한다. 이러한 이론과 일치하는 "임계" 유효응력은 사암의 종류에 따라 달라진다. 이 차이는 각각의 사암 종류의 미세구조에서 미세 균열 수의 차이에 기인한 것이라 생각된다. 높은 유효응력에서의 이론과 의미있게 일치하였으며, $CO_2$ 주입 시 현장에서의 탄성파 거동을 예상하는데 있어서 어느 정도 확신을 준다.

$CO_2$ 격리 처리를 위한 암석물리학 모의실헝장치와 그 응용 (A rock physics simulator and its application for $CO_2$ sequestration process)

  • Li, Ruiping;Dodds, Kevin;Siggins, A.F.;Urosevic, Milovan
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.67-72
    • /
    • 2006
  • 지하 염수층의 $CO_2$ 주입은 큰 저장 능력으로 인하여 대기 중으로의 $CO_2$ 방출을 감소시키기 위한 가장 유망한 방법일 것이다. $CO_2$ 저장은 적어도 수 천년 간 $CO_2$가 지층 안에 안전하게 남아있도록 주의깊게 계획되고 모니터링되어야 한다. 특히 해양 저류층에 대한 탄성파 탐사 방법들은 알맞은 저류층특성이 제공된다면 $CO_2$의 주인공정과 분산을 모니터링하기 위한 일차적인 수단이다. 탄성파탐사 방법은 잠재적인 트랩, 저류층 특성, 저류층 저장능력의 규명에 또한 필수적이다. 따라서 $CO_2$ 저장에 대한 탄성파 반응의 변화에 대한 평가는 매우 초기 단계에 이루어져야 한다. 이것은 모암과 $CO_2$ 사이의 화학적 작용에 의해 일어날 수 있는 유체의 특성이나 광물 조성의 변화에 따른 탄성파 반응에서의 잠재적 변화를 평가하기 위해 나중 단계에 다시 고려될 필요가 있다. 따라서 저류층에 일정시간 이상의 $CO_2$ 주입에 의한 탄성파 반응 변화에 대해 섬세히 구축된 모형은 장기간의 모니터링 프로그램 설계에 도움을 준다. 그러한 목적으로 주입된 $CO_2$에 대한 단기간과 장기간의 4차원 탄성파 반응을 모델링하도록 설계된, 그래픽 사용자 인터페이스((GUI)를 채택한 암석물리학 모의실험장치를 개발했다. 적용분야는 $CO_2$ 위상 변화, 국부적인 압력과 온도 변화, 화학 반응 및 광물의 침전을 포함한다. 이방성 가스만(Gassmann) 식을 모의실험장치에 고려시킴으로써 단층과 파쇄대를 재활성화 시키는 $CO_2$의 탄성파 반응 또한 예측될 수 있다. 이 논문에서는 암석물리학 모의실험장치를 적용했던 현장(해상과 육상의 잠재적 $CO_2$ 격리 지역)의 사례를 보여주고 있다. 4차원 탄성파 반응들이 모니터링 프로그램의 설계를 돕기 위하여 만들어 졌다.

고해상도 항공자력탐사를 이용한 Italia Vulcano-Lipari 화산 복합체의 천부 지하 구조 (Shallow subsurface structure of the Vulcano-Lipari volcanic complex, Italy, constrained by helicopter-borne aeromagnetic surveys)

  • 대웅무웅;중총 정;구택정부;삼원광언;중야 준
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.129-138
    • /
    • 2006
  • 남부 Italia Aeolian 군도의 Vulcano-Lipari 화산 복합체의 천부 지하구조를 잘 이해하고 또한 이 지역의 화산활동을 모니터링하기 위해서 고해상도 항공자력탐사가 3 년간의 간격을 두고 두 번 수행되었다. 두 개의 서로 다른 자력탐사 자료가 화산활동의 변화를 지시하는 어떠한 의미있는 차이를 보이지 않기 때문에, 자료들은 서로 합쳐져서 단일 자료보다 넓은 영역에 대한 항공자력도로 만들어졌다. 지형보정된 자력이상으로부터 겉보기 자화강도 분포도가 만들어졌으며 이로부터 Fossa 원추구의 이질성을 제시하는 국부적인 고 자화 이상을 볼 수 있었다. 이중 세 개의 고자화 이상에 대해 자력 모델링이 수행되었다. 각 모델은 Fossa 화구의 화산쇄설류로 덮혀있는 화산생성물의 존재를 밝히는 데에 적용되었다. Fossa 화구 지역에 대한 모델로부터 현재 화구의 남쪽 가장자리에는 조면암질 용암류가 묻혀있다는 것이 제시되었다. Forgia Vecchia에서 적용한 자력모델은 수증기 폭발성 원추구가, Fossa 칼데라를 메운 레타이트질 용암류(현무암질조면암과 안산암질조면암을 통칭)에 덮혀버린 한 분출중심으로부터 형성되었다는 것을 제시해 준다. 하지만 용암류의 분포는 기존의 시추 결과들로부터 알려진 것보다도 적은 지역에 국한되는 것처럼 나타난다. 이는 Porto Levante에 인접한 지열지역에서 알 수 있듯이, 강렬한 열수활동으로 인한 용암류의 부분적인 변질에 기인한 것으로 설명될 수 있다. Fossa 원추구 북동부에서의 모델은 두꺼운 용암류가 Fossa 화산활동의 초기단계에 또 다른 분출중심에서 집적되어 있다는 것을 암시한다. 최근의 전기탐사는 마지막 두 자력모델 지역에서 고비저항대를 보여준다.를 재활성화 시키는 $CO_2$의 탄성파 반응 또한 예측될 수 있다. 이 논문에서는 암석물리학 모의실험장치를 적용했던 현장(해상과 육상의 잠재적 $CO_2$ 격리 지역)의 사례를 보여주고 있다. 4차원 탄성파 반응들이 모니터링 프로그램의 설계를 돕기 위하여 만들어 졌다. 액체상의 $CO_2$ 주입은 공기로 포화된 상태에 비해 속도-유효응력 반응을 평균 약 8% 정도 낮게 한다. 실험자료들은 높은 유효응력에서 Gassmann 계산들과 일치한다. 이러한 이론과 일치하는 "임계" 유효응력은 사암의 종류에 따라 달라진다. 이 차이는 각각의 사암 종류의 미세구조에서 미세 균열 수의 차이에 기인한 것이라 생각된다. 높은 유효응력에서의 이론과 의미있게 일치하였으며, $CO_2$ 주입 시 현장에서의 탄성파 거동을 예상하는데 있어서 어느 정도 확신을 준다.극압 증가의 최대 허용치를 결정하는데는 사용할 수 없다고 주장하고자 한다. 초기폐사율이 낮음을 확인할수 있었으며, 상기 결과를 토대로, 넙치 치어의 경우 ${\beta}-1,3$ 글루칸을 0.05% 이상 0.1% 미만을 사료에 첨가하는 것이 성장, 사료효율 증진, 항산화능 및 질병저항성에 가장 좋은 효과를 나타낼 수 있을 것을 사료된다./Cip1}(-)/p27^{kip1}(-)$인 경우는 미만형인 경우(87.0%)가 장형(54.9%)의 경우보다 많은 비율을 차지하였다(P<0.05). 5년 장기 생존율에 있어서는 각각의 $p21^{Waf1/Cip1}$$p27^{kip1}$의 발현 유무에 따른 통계적인 유의성은 없었고 복합

일본 Nagaoka의 육상 대수층에 주입된 $CO_2$의 관찰을 위한 시간차 시추공간 탄성파 토모그래피 (Time-lapse crosswell seismic tomography for monitoring injected $CO_2$ in an onshore aquifer, Nagaoka, Japan)

  • ;신강 대;동 굉행;설 자구
    • 지구물리와물리탐사
    • /
    • 제9권1호
    • /
    • pp.30-36
    • /
    • 2006
  • 일본의 첫 번째 파일럿 규모의 $CO_2$격리실험이 Nagaoka 에서 약 1,100m 깊이의 육상 대수층에 10,400 톤의 $CO_2$를 주입하면서 행해졌다. 주입된 $CO_2$ 의 모니터링을 위해서 그 지역에서 행해진 다양한 측정들 중, P 파 속도의 변화를 이용하여 2 개의 관측정 사이에서 대수층 내의 $CO_2$ 분포를 알기 위해 시간차 시추공간 탄성파 토모그래피가 실시되었다. 이 논문은 탄성파 토모그래피에 대한 예비적인 결과들이다. 시추공간 탄성파 토모그래피 측정이 3 회에 걸쳐 실시되었는데, 한 번은 기준 조사로서 주입 전에, 나머지 두 번은 모니터링 조사를 위해 주입 중에 이루어졌다. 모니터링 조사로부터의 속도 분포도를 기준 조사 속도 분포도와 비교하여 속도 차 분포도가 만들어졌다. 속도 차 분포도는 주입된 $CO_2$가 퍼져있다고 예상되는 주입정 주변 대수층 부분에서 소도 감소 지역을 보여준다. 비록 탄성파 토모그래피가 $CO_2$주입에 따른 속도 감소 지역을 제공할 수 있었지만 속도 차 분포도에서 다소 이상한 속도 분포 지역이 관측되었다. 그러한 현상의 발생을 조사하기 위해서 조사 지역을 모사하는 수치모형 실험을 수행하였다. 그 결과 우리는 속도 차 분포도에서 발생한 이상한 속도 분포는 인위적인 오차이거나 이 지역에서 행해진 시추공간 탐사의 송수신기 배열에 의한 고스트라고 확신할 수 있었다. 토모그래피에 의해서 얻어진 최대 속도 감소 비율은 음파검층에 의해 관측된 것보다 아주 작았다. 비록 탄성파 토모그래피에 대한 수치모형 실험에 의해서도 불일치가 확인되었으나 음파검층에 의해 획득한 속도 값에 비교될 수 있는 좀 더 정확한 속도 값을 얻기 위해 앞으로 더 많은 연구가 필요하다.

심부 지열자원 개발을 위한 원거리 기준점 MT 탐사자료의 2차원 역산 해석 (Two-Dimensional Interpretation of Ear-Remote Reference Magnetotelluric Data for Geothermal Application)

  • 이태종;송윤호
    • 지구물리와물리탐사
    • /
    • 제8권2호
    • /
    • pp.145-155
    • /
    • 2005
  • 경상북도 포항시 흥해읍 일대에서 지열자원 개발을 위한 심부 파쇄대 탐지를 목적으로 수행된 자기지전류 (MT) 탐사에서 경북 안동과 충북 청원시의 대청호 및 일본의 Kyushu에 원거리 기준점을 설치하고 이를 이용한 자료처리 결과의 상호 비교를 통하여 원거리 기준점에 따른 겉보기 비저항과 위상 자료의 질을 검토하였다. 대상지역에서 약 165 km떨어진 대청호 자료를 이용한 경우는 전력주파수인 60 Hz대역과 자연 전자기장 신호가 미약한 $10^{-1}Hz\~1\;Hz$ 대역을 제외한 주파수 대역에서 매우 좋은 자료를 얻을 수 있었으며, 약 480 km 떨어진 Kyushu 원거리 기준점 자료를 이용해서 거의 모든 주파수 대역에서 연속성이 매우 좋은 전기비저항과 위상 곡선을 얻을 수 있었다. 이를 이용하여 2차원 역산을 수행한 결과, 지표의 10ohm-m 이내의 전기비저항을 보이는 반고결 이암층은 대상지역의 남쪽에서는 약 500 m 이상, 북쪽에서는 200 m 이내의 두께로 분포하며 북에서 남으로 경사진 형태로 나타났다. 심도 $500\~1,500m$에서 저비저항(L-2)과 고비저항(H-2)의 경계면은 단층면으로 해석되며 측점 206, 112와 414를 지나는 대략 $N15^{\circ}E$의 주향을 보인다. 1 km 이내의 천부의 저비저항 이상(L-4)은 파쇄대로 해석되며 측점 105를 지나면서 데략 $N60^{\circ}W$의 주향을 보인다. 또한, 북쪽의 큰 하천을 따라 저비저항 이상대가 공통적으로 나타나 구조선일 가능성을 시사한다. 그리고, 대상지역의 서쪽과 서남쪽의 $2\~3km$ 하부에서 10 ohm-m이내의 저비저항 이상대(L-3)가 나타나는데 이 층에 대해서는 향후 추가적인 연구가 요구된다.