$\require{mediawiki-texvc}$

연합인증

연합인증 가입 기관의 연구자들은 소속기관의 인증정보(ID와 암호)를 이용해 다른 대학, 연구기관, 서비스 공급자의 다양한 온라인 자원과 연구 데이터를 이용할 수 있습니다.

이는 여행자가 자국에서 발행 받은 여권으로 세계 각국을 자유롭게 여행할 수 있는 것과 같습니다.

연합인증으로 이용이 가능한 서비스는 NTIS, DataON, Edison, Kafe, Webinar 등이 있습니다.

한번의 인증절차만으로 연합인증 가입 서비스에 추가 로그인 없이 이용이 가능합니다.

다만, 연합인증을 위해서는 최초 1회만 인증 절차가 필요합니다. (회원이 아닐 경우 회원 가입이 필요합니다.)

연합인증 절차는 다음과 같습니다.

최초이용시에는
ScienceON에 로그인 → 연합인증 서비스 접속 → 로그인 (본인 확인 또는 회원가입) → 서비스 이용

그 이후에는
ScienceON 로그인 → 연합인증 서비스 접속 → 서비스 이용

연합인증을 활용하시면 KISTI가 제공하는 다양한 서비스를 편리하게 이용하실 수 있습니다.

지구통계학적 진화전략 역산해석 기법의 소개 및 가상 대수층 수리전도도 분포 예측에의 적용
Introduction of Inverse Analysis Model Using Geostatistical Evolution Strategy and Estimation of Hydraulic Conductivity Distribution in Synthetic Aquifer 원문보기

자원환경지질 = Economic and environmental geology, v.53 no.6, 2020년, pp.703 - 713  

박은규 (경북대학교 지구시스템과학부 수리지질학 연구실)

초록
AI-Helper 아이콘AI-Helper

수리지질학을 포함한 많은 지질학 분야에 있어 지하 매질의 불균질성을 규명하는 것은 큰 중요성을 가진다. 본 연구에서는 최근 Park(2020)에 의해 소개된 바 있는 대수층을 구성하는 매질의 수리물성을 예측할 수 있는 기법의 개념 및 이론을 간단히 소개하고, 가상의 대수층에 해당 기법을 적용하여 기법에 의한 결과들의 다양한 시사점을 도출하였다. 소개하는 기법은 공분산행렬 적응 진화전략이라는 광역최적화 기법을 사용하며, 개념적으로 대수층에 가해지는 수리적 스트레스에 의한 지하수위 변화 자료를 동화하여 대수층 불균질성을 특성화하는 방법론이다. 가상의 대수층의 수리전도도 예측에 개발 기법을 적용한 결과, 총 40000개 미지의 값을 매우 빠른 시간 내에 예측함을 확인하였다. 또한, 예측의 결과는 레퍼런스 수리전도도와 수치적 및 구조적으로 큰 유사성을 보여 예측의 질적 수준이 높음을 확인하였다. 본 연구에서는 매우 제한적인 케이스에 대하여 적용을 실시하였으나, 기법의 추가개발을 통하여 보다 다양한 케이스에의 적용이 가능할 것으로 예상되며 현재 이를 위한 추가 개발이 이루어지고 있는 상황이다. 개발 기법은 수리지질학 분야 뿐만 아니라 다양한 지질학 및 지구물리 분야에 적용될 수 있는 잠재성을 갖추고 있다.

Abstract AI-Helper 아이콘AI-Helper

In many geological fields, including hydrogeology, it is of great importance to determine the heterogeneity of the subsurface media. This study briefly introduces the concept and theory of the method that can estimate the hydraulic properties of the media constituting the aquifer, which was recently...

주제어

표/그림 (7)

AI 본문요약
AI-Helper 아이콘 AI-Helper

* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.

문제 정의

  • 모리 및 연산 효율성의 획기적인 개선을 위해서는 B와 D를 각각 주성분(principal components) 행렬 및 특이값(singular values) 행렬로 대체할 수 있다. 그러나, 본 논문에서는 불필요한 수학적 복잡성을 최소화하기 위하여 보다 상세한 전개는 생략하였다. 초기 공분산 행렬을 구성하기 위해서는 다양한 공분산 함수가 이용될 수 있다.
  • 따라서 본 논문은 Park(2020)에 의해 소개된 광역최적화 기반 대수층 물성 예측 역산 방법론에 대하여 소개하고, 방법론을 가상 대수층에 적용하는 예제를 통해 예측 효율성과 정확성을 확인하며, 현재까지 개발된 2차원 모델의 소프트웨어를 공개하고자 하는 목적을 지니고 있다.
  • 본 논문에는 내용의 간결성 및 핵심성을 위하여 주요 이론만을 서술하였다. CMA-ES 기법을 이용한 대수층 수리물성 변화 예측 이론과 관련된 보다 상세한 내용은 Park(2020)에 수록되어 있으며, 관심 있는 독자들은 해당 연구를 참조하기 바란다.
  • 모사에 이용된 소프트웨어는 USGS의 MODFLOW-2005 (Harbaugh, 2005)이며, 전산시스템의 운영체제인 OS X에 적합하도록 컴파일 하였으며 컴파일러로는 gfortran을 이용하였다. 본 논문에서는 개발 기법의 검증에 목적을 두고 있다. 따라서, 3차원 도메인, 복잡한 경계조건, 저유계수의 불균질성 등을 반영한 도메인은 본 연구의 고려대상이 아니다.
  • 본 저자는 이 논문을 통해 역산해석 기법을 이용한 대수층 수리물성 특성화 및 최적화 알고리듬에 의한 모델 보정(model calibration) 방법이 보다 일반적인 방법으로 정착되어 실용적인 목적으로 빈번히 이용되기를 바라는 바이며, 이에 본 연구에 이용된 소프트웨어 패키지를 일반에 공개하고자 한다.

가설 설정

  • 관측정에서의 수리전도도 측정값이 가용한 경우 본 연구에서와 같이 이를 반영하여 예측의 질을 향상시킬 수 있으나, 수리전도도 측정값이 없을 경우에도 개발된 기법을 통하여 수리전도도 분포를 예측할 수 있다. 다만, 본 연구에서는 수리전도도가 있는 경우만을 상정하여 예측을 수행하였다.
  • 경계조건의 경우, 도메인의 남측을 제외한 모든 경계는 불투수 경계이며 남측은 100 m의 고정수두 경계로 설정하였다. 도메인 내 지하수위의 변화는 강우에 의한 지하수 함양에 기인하여 발생하며, 인위적인 지하 수위 변화는 발생하지 않는 것으로 가정하였다. Figure 1b는 해당 도메인 지표에서 가정된 강수량을 표현한 것이다.
  • 모사에서 가정된 대수층은 자유면 대수층으로 대수층의 두께는 도메인 전반에 걸쳐 균질하게 100 m를가정하였다. 또한 도메인에 설정된 저유계수(storagecoefficient, S)는 전체에 균질하게 5×10-3 m-1설정하였으며, 이는 파쇄암반의 특성에 가깝다.
  • 그러나, 해당 연구에서 개발된 주요 이론은 비단 2차원에만 한정되지 않으며, 1차원 혹은 3차원으로 확장되어 적용할 수 있다. 앞서 언급한 바와 같이, 가상의 수리전도도 필드의 공간적 상관성은 가우시안 프로세스를 따르므로, 부드러운 공간적 변화를 가정한다. 기법에 의한 진화는 여러 세대(g)에 걸쳐 발생하며 g는 0(최초 예측)에서 gmax(최종 예측) 까지 순차적으로 증가한다.
  • 대수층의 수리전도도 분포는 Figure 2a를 통해 확인할 수 있다. 이는 GSLIB(Deutsch and Journel, 1992)의 순차 가우시안 모사(sequential Gausian simulation)을 이용하여 재현된 것으로 남-북 방향으로의 상관거리는 약 10 km 이며 동-서 방향의 상관거리는 약 2 km를 가정하여 제작되었다. Figure 2a에서는 로그 화(log-transformd)가 적용된 수리전도도(K)를 보여주고 있으며, logK 범위는 대략 –1에서 6 사이(K=0.
  • Figure 1b는 해당 도메인 지표에서 가정된 강수량을 표현한 것이다. 지하수 함양율은 강수량의 20%로 가정하였다.
본문요약 정보가 도움이 되었나요?

참고문헌 (18)

  1. Cirpka, O.A. and Kitanidis, P.K. (2000) Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data. Advances in Water Resources, v.24(1), p.89-103. 

  2. Deutsch, C.V. and Journel, A.G. (1992) Geostatistical software library and user's guide. New York, 119(147). 

  3. Evensen, G. (2003) The ensemble Kalman filter: Theoretical formulation and practical implementation. Ocean Dynamics, v.53(4), p.343-367. 

  4. Hansen, N. and Ostermeier, A. (2001) Completely derandomized self-adaptation in evolution strategies. Evolutionary Computation, v.9(2), p.159-195. 

  5. Harbaugh, A.W. (2005) MODFLOW-2005, the U.S. Geological Survey modular ground-watermodel - the Ground-Water Flow Process: U.S. Geological Survey Techniques and Methods 6-A16. 

  6. Jeong, J. and Park, E. (2019) Theoretical development of the history matching method for subsurface characterizations based on simulated annealing algorithm. Journal of Petroleum Science and Engineering, v.180, p.545-558. 

  7. Kang, P.K., Lee, J., Fu, X., Lee, S., Kitanidis, P.K. and Juanes, R. (2017) Improved characterization of heterogeneous permeability in saline aquifers from transient pressure data during freshwater injection. Water Resources Research, v.53(5), p.4444-4458. 

  8. Kitanidis, P.K. (1995) Quasi-linear geostatistical theory for inversing. Water Resources Research, v.31(10), p.2411-2419. 

  9. Kitanidis, P.K. and Lee, J. (2014) Principal Component Geostatistical Approach forlarge-dimensional inverse problems. Water Resources Research, v.50, p.5410-5427. 

  10. Lee, J. and Kitanidis, P. (2014) Large-scale hydraulic tomography and joint inversion of head and tracer data using the Principal Component Geostatistical Approach(PCGA). Water Resources Research, v.50, p.5410-5427. 

  11. Lee, J., Kokkinaki, A. and Kitanidis, P.K. (2018) Fast large-scale joint inversion for deep aquifer characterization using pressure and heat tracer measurements. Transport Porous Med., v.123(3), p.533-543. 

  12. Nash, J.E. and Sutcliffe, J.V. (1970) River flow forecasting through conceptual models part I-A discussion of principles. J. Hydrol., v.10(3), p.282-290. 

  13. Nowak, W. and Cirpka, O.A. (2006) Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data. Water Resources Research, v.42(8). 

  14. Park, E. (2020) A geostatistical evolution strategy for subsurface characterization: Theory and validation through hypothetical two-dimensional hydraulic conductivity fields. Water Resources Research, v.56. 

  15. Sanchez-Leon, E., Erdal, D., Leven, C. and Cirpka, O.A. (2020) Comparison of Two Ensemble Kalman-Based Methods for Estimating Aquifer Parameters from Virtual2-D Hydraulic and Tracer Tomographic Tests. Geosciences, v.10(7), p.276. 

  16. Wang, Z., Bovik, A.C., Sheikh, H.R. and Simoncelli, E.P. (2004) Image quality assessment: from error visibility to structural similarity. IEEE Transaction on Image Processing, v.13(4), p.600-612. 

  17. Yeh, T.C.J., Jin, M. and Hanna, S. (1996) An iterative stochastic inversemethod: Conditional effective transmissivity and hydraulic head fields. Water Resources Research, v.32(1), p.85-92. 

  18. Zha, Y., Yeh, T.C.J., Illman, W.A., Zeng, W., Zhang, Y., Sun, F. and Shi, L. (2018) A reduced-order successive linear estimator for geostatistical inversion and its application in hydraulic tomography. Water Resources Research, v.54(3), p.1616-1632. 

저자의 다른 논문 :

관련 콘텐츠

오픈액세스(OA) 유형

GOLD

오픈액세스 학술지에 출판된 논문

섹션별 컨텐츠 바로가기

AI-Helper ※ AI-Helper는 오픈소스 모델을 사용합니다.

AI-Helper 아이콘
AI-Helper
안녕하세요, AI-Helper입니다. 좌측 "선택된 텍스트"에서 텍스트를 선택하여 요약, 번역, 용어설명을 실행하세요.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.

선택된 텍스트

맨위로