• 제목/요약/키워드: Seismic data

검색결과 1,399건 처리시간 0.029초

원전 격납 건물의 실시간 모니터링을 위한 강건한 최적 센서배치 연구 (A Study on Robust Optimal Sensor Placement for Real-time Monitoring of Containment Buildings in Nuclear Power Plants)

  • 이찬우;김유진;정형조
    • 한국전산구조공학회논문집
    • /
    • 제36권3호
    • /
    • pp.155-163
    • /
    • 2023
  • 원전 구조물의 실시간 모니터링 기술이 요구되고 있지만, 현재 운영 중인 지진 감시계통으로는 동특성 추출 등 시스템 식별이 제한된다. 전역적인 거동 데이터 및 동특성 추출을 위해서는 다수의 센서를 최적 배치하여야 한다. 최적 센서배치 연구는 많이 진행되어 왔지만 주로 토목, 기계 구조물이 대상이었으며 원전 구조물 대상으로 수행된 연구는 없었다. 원전 구조물은 미미한 신호대잡음비에도 강건한 신호를 획득하여야 하며, 모드 기여도가 저차 모드에 집중되어 있어 모드별 잡음 영향을 고려해야 하는 등 구조물 특성을 고려해야 한다. 이에 본 연구에서는 잡음에 대한 강건도와 모드별 영향을 평가할 수 있는 최적 센서배치 방법론을 제시하였다. 활용한 지표로서 auto MAC(Modal Assurance Criterion), cross MAC, 노드별 모드형상 분포를 분석하였으며, 잡음에 대한 강건도 평가의 적합성을 수치해석으로 검증하였다.

$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$ 주입 시 현장에서의 탄성파 거동을 예상하는데 있어서 어느 정도 확신을 준다.

충청 동부지역 지각의 P, S파 감쇠 분석 (Attenuation of High-Frequency P and S Waves in the Crust of Eastern Part of Choongchung Provinces)

  • 경재복;김규동
    • 한국지구과학회지
    • /
    • 제24권8호
    • /
    • pp.684-690
    • /
    • 2003
  • 최근 새로운 행정수도의 이전이 예상된 남한 중부 충청지방에 대한 관심이 고조되고 있다. 우리는 이 지역에서 1996년 9월 이후 한국교원대학교 지진관측망을 설치하석 운영 중이며, 이 관측망 중 2개 관측소에서 기록된 60개의 지진자료를 바탕으로 확장 Coda 규격화법을 이용하여 Q$_P^{-1}$ 과 Q$_S^{-1}$을 동시에 측정하였다. 부족한 자료를 보충하기 위해서 한국지질자원연구원 관측망 중 한국교원대학교 관측망과 지리적으로 인접한 1개의 관측소에서 기록된 33개의 지진자료도 함께 분석하였다. 측정된 Q$_P^{-1}$은 (1.9${\pm}$3.0)${\times}$10$^{-3}$에서 (5.4${\pm}$1.5)${\times}$10$^{-4}$으로, Q$_S^{-1}$은 (2.4${\pm}$1.4)${\times}$10$^{-3}$에서 (6.3${\pm}$1.1)${\times}$10$^{-4}$으로 주파수가 3.0Hz에서 24Hz로 늘어남에 따라 줄어드는 주파수 의존적 특성을 보인다. 이 값을 주파수의 지수형태로 나타내면 Q$_P^{-1}$=0.003f^{-0.62}$, Q$_S^{-1}$=0.006f$^{-0.071}$이다. 이 값들은 지진학적으로 안정한 전 세계의 다른 지역의 값과 유사하며 남한 남동부 지역값에 비해서는 다소 약한 주파수 의존을 나타낸다. 이는 고주파수에서 높은 Q$^{-1}$값 때문이다. 남한 남동부와 충청 동부지역 간의 Q$^{-1}$의 차이를 평가하기 위해서는 보다 많은 관측이 요구된다.

남해 대륙붕 PZ-1 시추공 주변 현무암 대지 구조의 CO2 지중저장용량 평가 (Assessment of CO2 Geological Storage Capacity for Basalt Flow Structure around PZ-1 Exploration Well in the Southern Continental Shelf of Korea)

  • 신승용;강무희;신영재;정순홍
    • 자원환경지질
    • /
    • 제53권1호
    • /
    • pp.33-43
    • /
    • 2020
  • CO2 지중저장은 현재 온실가스 감축을 위한 CCS(Carbon Capture & Storage) 저장기술 중 가장 안정적이고 효과적인 기술로 평가되고 있다. 지중저장 대상 염대수층은 일반적으로 CO2가 초임계상태로 저장될 수 있는 지하 800 m 이상의 심도에 위치하고, 그 상부에 지표로의 CO2 누출을 막는 광역적인 불투수성 덮개암층이 분포해야 한다. 본 연구에서는 남해 대륙붕 탄성파 및 시추공 자료를 해석하여 덮개암으로 활용할 수 있는 현무암층과, 하부 염대수층이 발달하고 있는 현무암 대지 구조를 CO2 지중저장 유망구조로 제시하였다. 연구지역 주입대상 염대수층의 지중저장용량을 평가하기 위해 총 공극률은 시추공 감마 및 음파검층 자료를 이용하여 산출하였으며, 염대수층 구간의 온도/압력 조건을 고려하여 CO2 밀도를 계산하였다. 정확한 주입대상 염대수층의 체적을 산정하기 위해 x, y, z 축 방향으로 일정한 크기를 가지는 3차원 지질 격자 모델을 구성하였고, 염대수층 3차원 공간의 물성 분포를 지구통계학적 기법으로 예측하는 물성 모델링을 수행하여 총 공극률 값을 지질 격자에 할당하였다. U.S. DOE 방법을 이용하여 남해 대륙붕 현무암 대지 구조의 CO2 지중저장용량 평가 결과 평균 약 8,417만 CO2톤(최소 4,207만 ~ 최대 1억 4,379만 CO2톤)이 주입대상 염대수층 구간에 저장 가능한 것으로 예측되었다.

개정된 지반증폭계수의 Macro적 액상화 평가에 미치는 영향 분석 (Analysis of the Effect of the Revised Ground Amplification Factor on the Macro Liquefaction Assessment Method)

  • 백우현;최재순
    • 한국지반공학회논문집
    • /
    • 제36권2호
    • /
    • pp.5-15
    • /
    • 2020
  • 포항지진(ML=5.4) 시 발생한 액상화 현상은 국민들에게 지진으로 유발되는 액상화의 위험성을 새롭게 각인시켰고, 이에 대한 대비책으로 액상화 위험지도의 관심이 높아지고 있다. 현재 행정안전부가 보유하고 있는 액상화 위험지도는 2014년 제작된 것으로 전국 100,000개 이상의 시추 자료를 토대로 지하수위 0m인 조건으로 지반조건별 증폭계수를 사용하였으며 시추정보가 없는 지역은 보간법을 이용하여 2km × 2km 격자형식으로 제작된 것이 특징이다. 이러한 가운데, 2018년 행정안전부는 내진설계 공통기준의 새로운 지반분류법과 증폭계수를 공표하였다. 따라서 개정된 행정안전부의 증폭계수를 반영한 액상화 위험지도의 재작성이 필요하다. 본 연구는 내진설계 공통기준 개정 전·후 두 개의 기준으로 전 국토를 대상으로 지반분류를 수행하여 변동성을 분석하였으며, 지반조건별 증폭계수를 적용한 액상화 평가결과를 부산시 강서구를 대상으로 수행하였다. 이때 재현주기 500년과 1,000년에 해당하는 지반가속도를 적용하였으며 우리나라 평균 지하수위인 5m와 극한 조건인 0m로 구분하여 액상화 위험도를 평가하였다. 액상화 위험지도는 기존의 2km × 2km보다 높은 해상도를 확보하기 위해 500m × 500m 격자를 생성하여 위험지도를 작성하였다. 연구결과, 기존 지반분류 기준을 통해 SC, SD 지반으로 분류되었던 지반상태가 개정된 지반분류 기준을 통해 S2, S3, S4로 재분류되었다. 재현주기 500년과 1,000년으로 액상화 평가를 수행한 결과 개정 전 지반증폭계수 적용한 LPI가 상대적으로 과대평가되는 결과를 도출하였다. 본 연구결과는 증폭계수를 이용하는 광역지역 액상화 위험지도 작성의 근간인 액상화 평가에 큰 영향을 미치는 요소로써 향후 광역지역 액상화 위험지도 작성의 경우 반드시 고려될 사항으로 판단된다

한국 서남근해에 대한 항공자력탐사 해석 (Aeromagnetic Interpretation of the Southern and Western Offshore Korea)

  • 박창고;박창업
    • 한국석유지질학회지
    • /
    • 제2권2호
    • /
    • pp.51-57
    • /
    • 1994
  • 1969년 미해군이 황해지역에서 실시한 항공자력측정자료를 분석 연구한 결과, 한국 서남근해 제2광구 내에 있는 기존 군산분지(群山盆地) 남쪽에 흑산분지(黑山盆地)라고 명명하는 신퇴적분지가 존재함을 알게 되었다. 이 흑산분지는 서서 북-남남동으로 뻗은 3개의 소분지로 구성되었다. 이번 연구결과로 인하여 군산분지에 대해 독립적인 새로운 분석을 하게 했고 또한 그 두 분지뿐만 아니라 서근해 전역에 걸쳐 지질구조적 성인 및 그 생성역학(生成力學)에 대한 중요한 자료를 얻게 되었다. 분지생성 및 구조형태는 ‘두번-건너뛴-좌수향-렌치단층(斷層)’(Double overstepped left-lateral wrench faults)에 의한 롬보캐즘(rhombochasm)으로 해석된다. 한반도 내에 발달하고 있는 대단층들을 본 항공자력이상으로부터 얻은 중거에 의해서 본 연구지역 즉 근해까지 연장시켰다. 또한 이들 단층들의 성질을 좌수향렌치로 해석했다. 이와 같은 전 지역적인 단층 구조형태는 이동우의 영동분지에 대한 연구결과와 같다. 단층들의 운동방향은 필리핀이 에오세(Eocene) 이후 시계방향으로 약$80^{\circ}$회전했다는 최근 발표된 고지자기연구(古地磁氣硏究)에서 얻은 증거에 의하여 뒷받침되고 있다. 2$\frac{1}{2}$-차원의 모형 연구결과에 의하면 군산분지는 대략 $7,500{\cal}m$의 두터운 퇴적층을 갖는 것 같다. 흑산분지는 그 상부에 탄성파 속도가 빠른 층이 덮고 있어서 지금까지 발견되지 않고 있었다고 믿어진다. 흑산분지의 개략적인 형태가 석유부존에 적합한 듯 하기 때문에 두터운 퇴적층을 보유하고 있는 듯한 자력적 표현을 규명하기 위한 정밀항공 자력탐사를 실시할 것을 추천한다. 만약 고속도의 층이 확인된다면 자력탐사법에 대한 보충이 될 뿐 아니라 퇴적층에 대한 구조적인 묘사를 하는데 도움이 되는 중력탐사를 추가할 것을 추천한다. 정밀자력탐사 또는 중력탐사를 할 수 없다면 높은 에너지원의 탄성파 탐사를 실시하는 것이 더 효과적일 수도 있다.

  • PDF

Sentinel-1 InSAR 긴밀도 영상을 이용한 3월5일청년광산의 지표 변화 탐지 (Surface Change Detection in the March 5Youth Mine Using Sentinel-1 Interferometric SAR Coherence Imagery)

  • 문지현;김근영;이훈열
    • 대한원격탐사학회지
    • /
    • 제37권3호
    • /
    • pp.531-542
    • /
    • 2021
  • 노천 채광을 수행하는 광산은 지표 변화와 환경 교란을 발생시킬 수 있기 때문에 지속적인 모니터링이 필요하다. 노천 광산은 채광 작업장에 식생이 거의 분포하지 않아 InSAR 긴밀도 영상을 이용한 모니터링이 가능하다. 본 연구는 최근 개발된 InSAR 긴밀도 영상 기반의 Normalized Difference Activity Index(NDAI)를 적용하여 광산에서 발생하는 활동을 분석하였다. 3월5일청년광산은 2008년 이후 본격적으로 개발이 확장된 북한의 광산이다. 3월5일청년광산을 촬영한 12일 간격의 Sentinel-1 SAR 영상을 이용하여 획득된 InSAR 긴밀도 영상으로 NDAI 분석을 진행하였다. 우선 2000년부터 약 14년간 발생한 75.24 m의 고도 하강 지역과 약 9.85 m의 고도 상승 지역을 채광 작업장 및 광미 적치장으로 정의하였다. 이후 NDAI 영상을 이용하여 기간별 활동 분석을 진행하기 위해 전체 기간의 평균 영상, 1년 단위의 평균 영상, 및 4개월 단위의 평균 영상을 제작하였다. 2017년부터 2019년까지 광산 활동은 평균적으로 채광 작업장의 중심에서 비교적 활발하였다. 보다 자세한 광산의 활동 변화를 확인하기 위해 시간 간격을 좁혀 1년간의 활동을 알아보고자 하였다. 2017년은 지진파 자료의 정보와 NDAI 영상을 이용하여 인공 지진의 발생 시점과 그 전후에 대하여 RGB 합성 영상을 제작하고 채광 작업장의 활동 변화를 분석하였다. 2017년 4월 30일 발생한 대규모 발파 이후 채광 작업장의 서쪽에서 활발한 활동이 감지되었다. 9월 30일의 두 차례의 발파 이후에는 채광 작업장의 크기가 확장된 것으로 추정된다. 2018년 및 2019년의 활동 변화는 4개월 단위의 시간 평균 영상을 RGB 영상으로 합성하여 분석하였다. 연도별 활동을 분석한 결과, 2018년은 채광 작업장의 북동쪽에서 활발하게 활동하는 영역을 찾을 수 있었으며, 2019년은 광미 적치장에서 확장에 따른 특징적인 활동이 확인되었다. NDAI를 이용한 시계열 분석으로 광학 영상으로는 확인하기 어려운 노천 광산의 무작위적인 지표 변화를 탐지할 수 있었다. 특히 현장 자료를 획득할 수 없는 지역의 광산 활동을 원격 탐사를 이용하여 효과적으로 수행할 수 있었다.

황해 중동부 제4기 퇴적층의 지음향 특성 (Geoacoustic characteristics of Quaternary stratigraphic sequences in the mid-eastern Yellow Sea)

  • 진재화;장성형;김성필;김현태;이치원;장정해;최진혁;양우헌
    • 한국해양학회지:바다
    • /
    • 제6권2호
    • /
    • pp.81-92
    • /
    • 2001
  • 황해 중동부 해역 ($36^{\circ}00'N{\sim}36^{\circ}45'N$,$\;125^{\circ}00'E{\sim}125^{\circ}45'E$)에 분포하는 퇴적층을 대상으로 획득한 고해상 탄성파탐사자료(에어건, 스파커, SBP)와 심부 시추시료(YSDP 105, ${\sim}$64m 깊이)의 종합분석에 따라, 연구 해역의 층서모델을 설정하고 탄성 -음향 모델링을 수행하였다. 층서모델은 탄성파${\cdot}$${\cdot}$시 층서 단위들을 비교 분석하여 각 지층단위의 음향학적 특성, 퇴적과정, 생성시기를 규명하였다 또한 단위층 형성과정을 해수면 변동과 관련지어 순차층서를 복원하였다. 각 순차층은 육성 또는 천해성 기원이며 조립질의 저해수면 퇴적계와 조석 기원 세립질의 해침-고해수면 퇴적계로 구성되어 있다. 탄성-음향 모델링의 매개변수 산출을 위해, 0.S~90 cm 심도 간격으로 측정된 시추코아의 121 개 평균입도 자료를, Hamilton의 회귀식과 퇴적심도에 따른 변화구배를 이용하여, 전밀도(bulk density)값과 종파속도값으로 변환${\cdot}$계산하였다. 외삽${\cdot}$보정된 121 쌍의 전밀도${\cdot}$종파속도 물성값은 설정된 층서모델에서 제시된 단위층 별로 평균하여 각 단위층의 대표값으로 산출하였다. 퇴적층 내의 탄성-음향 전파 모델링을 위해, 층서모델의 각 단위층에 전밀도와 종파속도의 대표값을 입력 매개변수값으로 부여한 후, 층서 구성 단위층들을 유한요소격자로 분해하였다. 파선추적법의 컴퓨터 모델링 결과, 탄성-음향의 전파는 층서모델을 구성하는 단위층 경계면의 불규칙성과 저속도층의 다수 존재에 의하여 매우 복잡한 음경로를 보이며, 음원의 위치에 따라 탄성-음파의 암영대가 달리 존재함을 확인하였다. 것으로 보아 층 2는 고에너지환경에서 형성되었음을 시사한다. 또한 층 2에서 Al의 증가는 쇄설기원 물질의 유입에 의한 것으로 추정된다. 플라이오세 이후에 형성된 층 1은 치밀한 조직을 보이며, 형성 이후 속성작용을 받지 않았다.선택적으로 분해되는 것으로 여겨진다.두냔 tsuchigae, S. nigripinni morii, M. jeoni, C. splendidus, P. koreanus, C. lutheri, I. koreensis, C. herzi, R. brunneus 등으로 본 종들의 서식 상태는 하천의 오염 정도 및 하천개수와 직접적인 관계가 있으므로 어류서식 환경을 유지시키기 위해서 보다 적극적인 하천 수질관리의 대책과 방안이 수립되어야 될 것으로 사료된다.와 운동, H-R도상에서의 별의 특성과 진화 등이다. 탐구 과정에 대한 학생의 성취도는 비교적 높으므로 이해도가 상대적으로 낮은 영역에 대해 학생의 오답 유형을 참고하여 기본 개념 지도에 보다 관심을 기울일 필요가 있다.Cr은 FW에서 전반적으로 감소하는 경향을 보였지만 SW에서는 실험 초기에 감소하다 24시간 이후에는 증가 후 일정한 양상을 보였다. Pb은 FW에서 전반적으로 감소했지만 SW에서는 초기에 급격히 증가 후 다시 급격히 감소하는 양상을 보였다 Pb 또한 Cu, Cd, As와 마찬가지로 SW1&2에서 제거속도가 가장 빠르게 나타났다. FW 상층수 중 Hg는 시간에 따라 급격히 감소했고, 제거속도는 Fw5&6에서 가장 느렸다. 이러한 결과에 근거할 때 벼가 자라고 있고 이분해성 유기물이 풍부한 FW1&2, FW3&4 토양과 상층수에서는 유기물의 분해 활동이 활발하였지만, 벼가 경작되지 않는 FW5&6과 SW 에서는 유기물이 상대적으로 결핍되어 유기물의 분해활동이 적었을 것으로 판단된다. 한편, 수조에 인위적으로 유기물을 첨가한 경우 박테리아 세포수는 SW1에서 164시간 동안 4배 증가하였으나 SW3과 SW5에서는 각각 2.7배, 1.

  • PDF

북한 핵실험 모니터링 : 1, 2차 비교 (Monitoring North Korea Nuclear Tests: Comparison of 1st and 2nd Tests)

  • 지헌철;박정호;김근영;제일영;신동훈;신진수;조창수;이희일
    • 지구물리와물리탐사
    • /
    • 제13권3호
    • /
    • pp.243-248
    • /
    • 2010
  • 북한지역에서 핵실험으로 추정되는 두 번의 발파가 관측되었다. 한국지질자원연구원 관측소와 한중 공동관측소는 북한과 주변국간의 경계에 고르게 분포하고 있다. 본 연구에서는 북한 핵실험 장소로부터 200 km에서 550 km 거리에 있는 광대역 지진 관측소의 자료를 사용하여 북한의 2차례 핵실험을 비교 분석하였다. 관측소별 1차 실험과 2차 실험의 초동 Pn 도착 시간차를 비교함으로서 상대적인 위치이동을 계산할 수 있다. Pn 속도를 8 km/s로 가정하고, 실험 장소와 관측소간의 기하학적인 관계를 이용하여 계산한 결과, 2차 장소는 1차 장소로부터 서북서 방향으로 2 km 거리에 위치하는 것으로 추정된다. P 파로부터 계산된 2차 실험의 실체파 규모는 평균적으로 4.5이나, 관측소별로는 최대 5.2에서 최소 4.1로 아주 큰 차이를 보인다. 이에 비해 Lg 파로부터 계산한 2차 실험의 규모는 평균적으로 4.6이며, 관측소별로 최대 4.7에서 최소 4.3사이로 P 파에 의한 규모에 비해 관측소간의 차이가 작다. 1, 2차 실험의 이동 윈도우 주파수 스펙트럼은 매우 비슷한 패턴을 보여 주며 두 실험의 초동 P 파의 모서리 주파수는 거의 차이가 없다. 따라서 2차 실험의 깊이가 1차 때와 비슷한 것으로 추정된다. 2차 실험의 폭발력은 관측소별 1차와 2차의 지반속도비로부터 계산한 결과 1차에 비하여 8배 큰 것으로 추정된다.

그래픽 프로세서를 이용한 탄성파 수치모사의 계산속도 향상 (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배 향상시킬 수 있었다.