Uchida Toshihiro;Song Yoonho;Mitsuhata Yuji;Lee Seong Kon
한국지구물리탐사학회:학술대회논문집
/
한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
/
pp.390-397
/
2003
A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.
본 논문에서는 지진파 역행 전파를 이용하여 지진원을 영상화 하는 방법에 대한 연구를 수행하였다. 시간 역행 전파는 파동 방정식의 시간 반전의 불변성과 상반정리로 설명할 수 있으며 파동의 시간 역행 전파는 의학적인 목적을 비롯한 비파괴 검사, 그리고 탄성파 파형 역산에 사용되어 매질을 영상화 하는데 활용되고 있다. 지진 관측소에서 관측된 지진파형은 지진원으로부터 방출된 에너지를 기록한 것이다. 이를 시간의 역순으로 매질로 전파시키면, 지진으로부터 전파되었던 에너지가 진원시각의 지진원 위치에 집중하게 되어 지진원을 영상화 시키는 원리이다. 본 연구에서는 3차원 지진파 전파의 수치 모의를 위해 엇갈린 격자 유한 차분법과 병렬 처리 기법을 이용하여 탄성 매질에서의 파동방정식의 해를 구하였다. 수치모의 결과로부터 지진파 역행 전파를 이용해 지진의 시공간적 발생양상을 살펴볼 수 있음을 확인하였다.
곡면의 경계를 가지는 수로에서 비선형 파랑의 상호작용을 모의하기 위한 비정수압 자유수면 모형이 개발되었다. 제안된 모형은 비선형의 3차원 Navier-Stokes 방정식을 곡선좌표 영역에서 계산단계 분리법의 일종인 압력수정법에 의하여 수치적으로 해석된다. 특히, 연직방향으로 변형된 형태의 엇갈린 격자를 이용하여 상대적으로 간단하게 압력방정식과 자유수면 경계조건을 구성하였다. 개발된 모형의 수치해석 정확도는 2차원의 수치 파수조에서 파랑의 비선형 정도에 대하여 5차의 스토스우크스 해석해와 비교하였다. 본 모형의 실제적 적용은 원형의 수로에서 회절과 반사에 의해 변형되는 비선형 파의 변형에 초점을 맞추어 수행하였다. 두개의 파를 중첩한 고비선형의 파에 대한 경우를 제외하고 수치해석 결과는 비선형적인 영향을 고려하지 않은 해석해의 선형적인 중첩과 일치하였다. 두개의 파를 중첩한 고비선형의 파에 대한 모의를 통하여 본 모형은 원형의 수로에서 비선형 군파의 변형에관한 수치적인 모의 가능성을 제시하였다.
지표 레이다(Ground Penetrating Radar; GPR) 탐사를 모사하기 위해서, 시간 영역 유한 차분법을 이용하여 송수신 안테나를 고려한 3차원 모델링 알고리듬을 개발하였다. 국내에서 주로 사용되는 GPR탐사기기인 쌍극자 안테나를 모사하고, 안테나 끝점에서의 내부반사(ringing) 현상을 줄이기 위해 subcellular법을 사용하였다. 반무한 매질에서 종단 저항에 따른 출력 전압의 변화를 비교하여, 종단 저항의 수가 많아지면 내부반사가 줄어드는 반면 출력 전압이 작아진다는 사실을 알 수 있었고, 이 결과를 토대로 내부반사를 줄이는 쌍극자 안테나를 구현하였다. 각도에 따른 전기장 분포를 알기 위해 안테나에 직교하는 면과 안테나를 포함하는 면에서 방사 패턴을 구하였다. 안테나를 포함하는 면의 전기장 에너지가 안테나에 직교하는 면보다 안테나 직하부로 집중되는 현상을 볼 수 있었다.
본 연구는 C형 CIP법을 이용한 3차원 시간영역 음장해석법의 정밀도에 대해서 세밀한 검토를 하였다. 3차원 C형 CIP 음장해석의 위상오차 특성과 전파방향에 따른 오차를 명확히 하고, 본 수치 해석법의 유효성을 나타내었다. 또 종래의 수치해석법으로 staggered-grid 모델을 이용한 FDTD법에 따른 계산결과와 M형 CIP법에 의한 계산결과와 비교를 하였다. 본 논문의 검토에 의해 같은 이산화조건에서는 C형 CIP법은 M형 CIP법 및 FDTD법보다 해석법이 가진 분산성이 적고, C형 CIP법으로 계산된 음압 파형은 변형이 적은 결과가 되는 것을 알 수 있었다. 그리고 C-1형 CIP법, C-2형 CIP법, M형 CIP법 및 FDTD법의 메모리와 계산시간을 비교하였다. 그 결과 C형 CIP해석은 FDTD해석에 비하여 수치분산성이 적지만 많은 사용 메모리와 계산시간이 필요하였다. C-1형 과 C-2형 CIP법은 입방체의 격자 대각 방향에서는 축 방향에 비하여 정밀도가 약간 저하하는 것을 알 수 있었다. 그리고 C-2형 CIP법은 C-1형 CIP법보다 사용 메모리와 계산시간이 적고, 계산정밀도도 고려해서 보면 유효한 해석법이라는 것이 명확하였다.
자연 전자기장을 이용하여 지하 매질의 전기적 구조를 규명하는 자기지전류(magnetotelluric; MT) 탐사의 정확한 해석을 위해서는 특정 전기적 구조에 대한 정확한 수치적 반응을 구할 수 있는 3차원 모델링이 필수적이다. 특히, 매질내에 전기적 이방성이 있을 때는 MT 반응이 달라지므로 전기적 이방성의 영향을 고려한 MT 탐사 모델링이 필요하다. 특히, MT 탐사기법을 이용한 지열저류층의 모니터링과 같이 MT 반응의 작은 변화를 분석해야 하는 시간경과 자료의 해석의 경우, 대상 지역에 이방성이 존재할 경우 이를 고려할 수 있는 정확한 모델링이 필수적이다. 이 연구에서는 기존의 등방성만을 고려하던 유한차분법 MT 모델링 알고리듬을 수직 혹은 수평 횡등방성 이방성을 고려할 수 있도록 개선하였다. 개발한 알고리듬을 박리층 모델을 이용하여 검증한 후, 수직횡등방성 이방성이 MT 반응에 미치는 영향에 대해서 분석하였다. 향후에는 수평 횡등방성 이방성이 MT 반응에 미치는 영향에 대해서도 분석하고자 하며, 알고리듬을 더욱 발전시켜 경사 횡등방성 이방성까지 고려할 수 있도록 발전시키고자 한다.
해상에 누출된 기름의 거동은 크게 초기확산, 수평이동 그리고 풍화과정으로 나눌 수 있다. 본 논문에서는 초기확산과 수평이동을 고려하는 단기예측모델을 개발하여 경기만에 적용하였다. 누출된 기름의 초기확산과정에서는 Fay의 3단계 분류법에 근거한 경험식으로 확산반경을 추정하고, 수평이동과정에서는 바람과 $M_2$조석을 고려하여 기름의 이동을 추정 하였다. 조석 수치모델을 위해서 깊이방향으로 적분된 운동방정식을 사용하였고, 수치 계산을 위해 교차 격자체계를 채택하였다. 시간전개를 하는 방법으로 양해법과 중간시간전개법을 이용하여 그 계산결과를 서로 비교하였다. 수치해의 정확도를 검증하기 위하여 간단한 직사각형 유체영역에 먼저 수치모델을 적용한 다음, 실제 계산영역인 경기만에 적용하였다. 또한, 사고 발생시간과 누출장소, 누출된 양과 누출 지속시간, 바람의 크기와 방향을 입력하면 기름의 이동방향과 오염지역 그리고 해안선의 오염정도를 가시화할 수 있는 기능을 단기 예측모델에 포함 시켰다.
MT 모델링에서는 송신원을 고려하지 않으므로, 평면파에 대한 배경 매질의 반응을 경계값으로 설정하는 Dirichlet 경계조건을 이용할 때에 그 경계값의 정확한 계산이 매우 중요하다. 이 연구에서는 1차원 배경 매질만을 가정하던 기존의 모델링 알고리듬을 2차원 배경 매질을 고려할 수 있도록 발전시켰다. 1차원 배경매질의 경우 경계값은 해석적으로 계산할 수 있으나, 2차원 구조가 존재하는 경우에는 2차원 모델링을 통해 경계값을 계산하여야 한다. 2차원 모델링의 TM(transverse magnetic) 및 TE (transverse electric) 모드는 3차원 모델링의 입사 전기장의 분극 방향과 2차원 구조의 주향에 따라서 결정된다. 전기장을 셀 모서리에 정의하는 기존의 3차원 모델링 알고리듬과 잘 부합하도록 2차원 모델링에서도 모서리에서 전기장을 계산하였다. 2차원 모델링을 통해 계산된 값을 3차원 모델링의 경계값으로 활용한 결과, 단층 모형 혹은 한 면에 바다를 포함한 모형에 대해 보다 정확한 겉보기비저항 및 위상을 얻을 수 있었다.
탄화수소 생산에서 화학공법은 폴리머 용액을 저류층에 주입함으로써 탄화수소의 생산량을 증대시키는 방법이다. 화학공법의 유용성을 파악하기 위해서는 저류층 내 폴리머 용액의 거동 양상을 모니터링하는 것이 매우 중요하다. 저류층 내 폴리머 용액과 유체 변화에 대한 모니터링의 가능성을 확인하기 위해, 저류층을 구성하고 있는 암석 및 유체의 종류, 포화도 등의 매개변수에 따른 탄성파 및 전기비저항 암석물리모델을 구성하고, 폴리머 주입에 따른 탄성파 및 EM 토모그래피의 반응을 계산하였다. 특히 전기비저항 암석물리모델 구성에서는 순수사암, 셰일질 사암, 사암셰일박리층 3가지의 서로 다른 저류층 암석을 고려하였다. 폴리머 용액의 농도는 가장 일반적으로 이용되는 2 wt%로 하고 담수를 이용한 경우와 염수를 이용한 경우로 나누어 분석하였다. 더 나아가, 화학공법에 대한 중성자검층 민감도 분석을 수행하였다. 이 논문에서 다루는 기법들은 탄화수소 저류층뿐만 아니라 지열 저류층 등의 모니터링에 있어서도 중요하게 적용될 것이다.
탐사 지구물리학에서 수치 모사는 지하매질에서의 탄성파 전파 현상을 이해하는데 중요한 통찰력을 제공한다. 탄성파 모사는 음향파 근사에 의한 수치 모사보다 계산시간이 많이 소요되지만 전단응력 성분을 포함하여 보다 현실적인 파동의 모사를 가능하게 한다. 그러므로 탄성파 모사는 탄성체의 반응을 탐사하는데 적합하다고 할 수 있다. 계산 시간이 길다는 단점을 극복하기 위해 본 논문에서는 그래픽 프로세서(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배 향상시킬 수 있었다.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.