1.Nanjing Hydraulic Research Institute,National Key Laboratory of Water Resources and Hydraulic Engineering, Nanjing,210029,China; 2.Key Laboratory of Surficial Geochemistry,Ministry of Education,School of Earth Sciences and Engineering,Nanjing University,Nanjing,210023,China; 3.General Investigation Group of North China Geological Prospecting Bureau,Langfang,065201,China
Inversion of hydrogeological parameters is an important issue for the groundwater numerical simulation.The pilotpoints method combined with SVDAssist regularization is adapted in this paper to solve the highdimensional illposed problems by setting a number of pilot points in the study area.Base on a synthetic case,the inversion result shows that the pilotpoints method combined with SVDAssist regularization is efficient to estimate the spatial distribution of the hydraulic conductivity for the heterogeneous media,which does not rely on the artificially zone.The effect of the pilotpoint number on the parameter inversion results is also discussed and the results show the accuracy improves with the increase of the pilotpoint number,but deteriorates when the number exceeds the critical number.The computational cost increases with the number of pilot points.The number of pilot points should be controlled within a reasonable range.
Jiang Beilei1*,Feng Yansheng2,3,Shi Xiaoqing2,Wu Jichun2,Sun Yuanyuan2.
Inversion of hydrogeological parameters using pilotpoints and regularization method [J]. Journal of Nanjing University(Natural Sciences), 2016, 52(3): 438-447
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] 薛禹群,谢春红.地下水数值模拟.北京:科学出版社,2007,1-5.(Xue Y Q,Xie C H.Numerical simulation of groundwater.Beijing:Science Press,2007,1-5.) [2] 陈 彦,吴吉春.含水层渗透系数空间变异性对地下水数值模拟的影响.水科学进展,2005,16(4):482-487.(Chen Y,Wu J C.Effect of the spatial variability of hydraulic conductivity in aquifer on the numerical simulation of groundwater.Advances in Water Science,2005,16(4):482-487.) [3] 施小清,吴吉春,袁永生.渗透系数空间变异性研究.水科学进展,2005,16(2):210-215.(Shi X Q,Wu J C,Yuan Y S.Study on the spatial variability of hydraulic conductivity.Advances in Water Science,2005,16(2):210-215.) [4] 陆 乐,吴吉春.地下水数值模拟不确定性的贝叶斯分析.水利学报,2010,41(3):264-271.(Lu L,Wu J C.Bayesian analysis of uncertainties in groundwater numerical simulation.Journal of Hydraulic Engineering,2010,41(3):264-271.) [5] 吴吉春,陆 乐.地下水模拟不确定性分析.南京大学学报(自然科学),2011,47(3):227-234.(Wu J C,Lu L.Uncertainty analysis for groundwater modeling.Journal of Nanjing University(Natural Sciences),2011,47(3):227-234.) [6] Neuman S P.Calibration of distributed parameter groundwater flow models viewed as a multipleobjective decision process under uncertainty.Water Resources Research,1973,9(4):1006-1021. [7] Yeh W W G.Review of parameteridentification procedures in groundwater hydrology:The inverse problem.Water Resources Research,1986,22(2):95-108. [8] 刘玉珍,王本德,姜英震.基于可变水文地质参数的地下水系统数学模型.水科学进展,2009,20(3):398-402.(Liu Y Z,Wang B D,Jiang Y Z.Mathematical model for groundwater system based on the variable hydrogeological parameters.Advances in Water Science,2009,20(3):398-402.) [9] 周海燕.基于集合卡尔曼滤波法的非高斯含水层参数识别.博士论文.北京:中国地质大学,2011.(Zhou H Y.Charactering nonGaussian aquifer model parameters based on the ensemble Kalman Filter.Ph.D.Dissertation.Beijing:China University of Geosciences,2011.) [10] Doherty J,Fienen M,Hunt R.Approaches to highly parameterized inversion:Pilotpoint theory,guidelines,and research directions.US Geological Survey Scientific Investigations Report,2010-5168,1-36. [11] De Marsily G,Lavedan C,Boucher M,et al.Interpretation of interference tests in a well field using geostatistical techniques to fit the permeability distribution in a reservoir model,Geostatistics for Natural Resources Characterization:Part 2,1984,831-849. [12] Certes C,de Marsily G.Application of the pilotpoints method to the identification of aquifer transmissivities.Advances in Water Resources,1991,14(5):284-300. [13] LaVenue A M,Pickens J F.Application of a coupled adjoint sensitivity and Kriging approach to calibrate a groundwater flow model.Water Resources Research,1992,28(6):1543-1569. [14] RamaRao B S,LaVenue A M,de Marsily G,et al.Pilotpoint methodology for automated calibration of an ensemble of conditionally simulated transmissivity fields—1,Theory and computational experiments.Water Resources Research,1995,31(3):475-494. [15] LaVenue A M,RamaRao B S,de Marsily G,et al.Pilotpoint methodology for automated calibration of an ensemble of conditionally simulated transmissivity fields—2.Application:Water Resources Research,1995,31(3):495-516. [16] LaVenue M,de Marsily G.Threedimensional interference test interpretation in a fractured aquifer using the pilotpoint inverse method.Water Resources Research,2001,37(11):2659-2675. [17] Doherty J.Groundwater model calibration using pilotpoints and regularization.Ground Water,2003,41(2):170-177. [18] Burrows W,Doherty J.Efficient calibration/uncertainty analysis using paired complex/surrogate models.Ground Water,2015,53(4):531-541. [19] Moeck C,Hunkeler D,Brunner P.Tutorials as a flexible alternative to GUIs:An example for advanced model calibration using Pilot Points.Environmental Modelling & Software,2015,66:78-86. [20] Tikhonov A N.Regularization of incorrectly posed problems.Soviet Mathematics Doklady,1963,4:1624-1627. [21] Tikhonov A N,Arsenin V I A.Solution of illposed problems.New York,Halsted,1977,1-258. [22] Lawson C L,Hanson J R.Solving least squares problems.Classics in Applied Mathematics,Society for Industrial and Applied Mathematics,Philadelphia,Pennsylvania,1995,18-22. [23] Tonkin M,Doherty J.A hybrid regularized inversion methodology for highly parameterized models.Water Resources Research,2005,41,W10412,doi:10.1029/2005WR003995. [24] Doherty J.PEST,Modelindependent parameter estimation—User manual.The 5th Edition with slight additions.Brisbane,Australia,Watermark Numerical Computing,2010,1-336. [25] Doherty J.Guidelines for groundwater model calibration using regularized inversion.Brisbane,Australia,Watermark Numerical Computing,2008,1-41. [26] Seber G A F,Wild C J.Nonlinear regression.John Wiley,New York,1989,680-681. [27] Harbaugh A W,Banta E R,Hill M C,et al.MODFLOW2000,the U.S.geological survey modular groundwater model—User guide to modularization concepts and the groundwater flow process.US Geological Survey Open File Report 00-92,2000,1-121. [28] Deutsch C V,Journel A G.GSLIB:Geostatistical software library and users guide.The 2nd Edition.New York:Oxford University Press,1997,1-369.