A Solute Transport Analysis around Underground Storage Cavern by using Eigenvalue Numerical Technique

고유치 수치기법을 이용한 지하저장공동 주위의 용질이동해석

  • Chung, Il-Moon (Hydrology Research Division, Korea Institute of Construction Technology) ;
  • Kim, Ji-Tae (Hydrology Research Division, Korea Institute of Construction Technology) ;
  • Cho, Won-Cheol (Department of Civil and Environmental Engineering, Yonsei University) ;
  • Kim, Nam-Won (Hydrology Research Division, Korea Institute of Construction Technology)
  • 정일문 (한국건설기술연구원 수문연구실) ;
  • 김지태 (한국건설기술연구원 수문연구실) ;
  • 조원철 (연세대학교 공과대학 토목환경공학과) ;
  • 김남원 (한국건설기술연구원 수문연구실)
  • Published : 2008.12.30

Abstract

The eigenvalue technique is introduced to overcome the problem of truncation errors caused by temporal discretization of numerical analysis. The eigenvalue technique is different from simulation in that only the space is discretized. The spatially discretized equation is diagonized and the linear dynamic system is then decoupled. The time integration can be done independently and continuously for any nodal point at any time. The results of eigenvalue technique are compared with the exact solution and FEM numerical solution. The eigenvalue technique is more efficient than the FEM to the computation time and the computer storage in the same conditions. This technique is applied to the solute transport analysis in nonuniform flow fields around underground storage caverns. This method can be very useful for time consuming simulations. So, a sensitivity analysis is carried out by using this method to analyze the safety of caverns from nearly located contaminant sources. According to the simulations, the reaching time from source to the nearest cavern may takes 50 years with longitudinal dispersivity of 50 m and transversal dispersivity of 5 m, respectively.

본 연구에서는 수치해석의 시간차분에서 발생하는 절단오차를 극복하는 방안으로 고유치 기법을 도입하였다. 고유치기법은 모의를 할 때 공간만을 이산화하는 특징을 가지며, 공간적으로 이산화된 방정식을 대각화시킴으로써 선형동력학적 시스템을 분리시킨 후 시간적분을 이용한 계산이 임의의 위치에서 임의의 시간에 대해 개별적으로 또 연속적으로 수행된다. 이러한 고유치기법을 이용하여 오염물 이동을 모의하고 이를 해석해와 비교 검증하였고, 동일한 조건에서 유한요소법을 이용한 수치모형과 고유치 기법을 이용한 용질이동의 예측을 실시한 결과 고유치기법을 이용할 경우 계산시간과 저장용량이 수치모형에 비해 절약됨을 확인할 수 있었다. 고유치 기법을 이용하여 지하유류저장 공동주위의 불균일 유속장에서 용질의 이동을 분석하였다 이 방법이 모의발생에 오랜 시간이 걸리는 문제에 유용하게 사용될 수 있으므로, 공동에 인접한 오염원으로부터 공동의 안전성을 평가하기 위한 민감도 분석에 이 방법을 적용하였으며, 모의결과에 의하면, 종분산지수와 횡분산지수가 각각 50 m, 5 m일 때 공동에 도달하는 시간은 약 50년으로 추정되었다.

Keywords

References

  1. 정일문, 조원철, 배덕효, 1997, 지하공동주변의 지하수 흐 름해석을 위한 수치모형의 확립. 한국수자원학회 논 문집, 30(1), 63-73
  2. 정일문, 조원철, 허준행, 1998, 지하저유공동의 기밀성을 위한 수막시설의 수리학적 분석, 대한토목학회 논문 집, 18(II-2) 131-139
  3. 조원철, 정일문, 1997, 지하공동의 기체누출방지를 위한 수막시설 설치에 따른 지하수 흐름해석, 97년 비축기 지 건설기술 세미나 논문집. 한국석유개발공사
  4. 한국석유개발공사, 1985, A사업조사설계용역 기본설계보 고서
  5. 한국석유개발공사 1991, A비축기지 기본설계보고서
  6. 한일영, 서일원, 오경택, 1996, 수치모의를 통한 지하 LPG저장시설에 인접한 폐기물 매립지에서의 침출수 이동예측 및 제어공법검토. 한국지하수환경학회지, 3(2), 51-59
  7. Bellman, R. E. 1960, Introduction to Matrix Analysis, McGraw Hill
  8. Chung, I. M., Cho, W. C., and Heo, J. H. 2003, Stochastic hydraulic safety factor for gas containment in underground storage caverns, Journal of Hydrology, Volume 284, Issues 1-4, 77-91 https://doi.org/10.1016/S0022-1694(03)00275-0
  9. Guymon, G. L. 1970, A finite element solution of the one-dimensional diffusion-convection equation. Water Resources Research, 6(1), 204-210 https://doi.org/10.1029/WR006i001p00204
  10. Hwang, J. C., Cho, W. C., and Yeh, G. T. 1984, An eigenvalue solution continuous in time to the spatially discretized solute transport equation in steady groundwater flow. Water Resources Research, 20(11), 1725-1732 https://doi.org/10.1029/WR020i011p01725
  11. Kim, J. T., Cho, W. C., Chung, I. M., and Heo, J. H. 2007, On the stochastic simulation procedure of estimating critical hydraulic gradient for gas storage in unlined rock caverns, Geosciences Journal, 11(3), 185-270 https://doi.org/10.1007/BF02913932
  12. Cleary, R. W. and Ungs, M. J. 1978, Analytical models for groundwater pollution and hydrology, Water Resources Program Report, 78-WR-15
  13. Gureghian, A. B., Ward, D. S., and Cleary, R. W., 1980, A finite element model for the mitigation of leachate from a sanitary landfill on Long Island, New York, Water Resources Bulletin, 16(5) 900-906 https://doi.org/10.1111/j.1752-1688.1980.tb02506.x
  14. Kuiper, L. K. 1973, Analytical solution of spatially discretized groundwater flow equations. Water Resources Research, 9(4), 1094-1097 https://doi.org/10.1029/WR009i004p01094
  15. Pickens, J. F. and Lennox, W. C. 1976, Numerical simulation of waste movement in steady groundwater flow systems. Water Research, 12(2), 171-180 https://doi.org/10.1029/WR012i002p00171
  16. Sahuquillo, A. 1983, An Eingenvalue numerical techique for solving unsteady linear groundwater models continously in time. Water Resources Research, 19(1),87-93 https://doi.org/10.1029/WR019i001p00087
  17. Smith, B. T. 1976, Matrix Eigensystem Routines-EISPACK Guide, 2nd Edition, Springer-Verlag
  18. Yeh, G. T. and Ward, D. S. 1981, FEMWASTE : A Finite- Element Model of Waste Transport through Saturated- Unsaturated Porous Media, Oak Ridge National Laboratory Environmental Science Division Publication ORNL-5601