南京大学学报(自然科学版) ›› 2016, Vol. 52 ›› Issue (3): 438–447.

• • 上一篇    下一篇

 向导点-正则化方法在水文地质参数反演中的应用

 姜蓓蕾1*,冯燕生2,3,施小清2,吴吉春2,孙媛媛2   

  • 出版日期:2016-07-02 发布日期:2016-07-02
  • 作者简介:1.南京水利科学研究院 水文水资源与水利工程科学国家重点实验室,南京,210029;
    2.表生地球化学教育部重点实验室,南京大学地球科学与工程学院,南京,210023;
    3.华北地质勘查局综合普查大队,廊坊,065201
  • 基金资助:
     国家自然科学基金(41172206,41302201)

 Inversion of hydrogeological parameters using pilot­points and regularization method

 Jiang Beilei1*,Feng Yansheng2,3,Shi Xiaoqing2,Wu Jichun2,Sun Yuanyuan2   

  • Online:2016-07-02 Published:2016-07-02
  • About author: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

摘要:  水文地质参数反演是地下水数值模拟中的一个重要问题.基于向导点方法,通过设定理想算例,引入SVD­Assist正则化方法,反演求解向导点个数大于观测数据个数的高维不适定问题.与传统分区方法的反演结果对比结果表明:向导点-正则化法是一种有效的高维水文地质参数的反演方法,可不依赖于人为的概化分区有效地反演推估非均质渗透系数场的空间分布.还讨论了向导点个数对反演结果的影响,结果表明:向导点个数并非越多越好,随着向导点个数的逐步增加,反演结果精度逐渐提高,但达到一定数量后,向导点个数的增加不仅会增加计算工作量,同时反而会降低反演结果精度.

Abstract:  Inversion of hydrogeological parameters is an important issue for the groundwater numerical simulation.The pilot­points method combined with SVD­Assist regularization is adapted in this paper to solve the high­dimensional ill­posed problems by setting a number of pilot points in the study area.Base on a synthetic case,the inversion result shows that the pilot­points method combined with SVD­Assist 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 pilot­point number on the parameter inversion results is also discussed and the results show the accuracy improves with the increase of the pilot­point 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.

 [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 multiple­objective decision process under uncertainty.Water Resources Research,1973,9(4):1006-1021.
[7] Yeh W W G.Review of parameter­identification 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 hydro­geological para­meters.Advances in Water Science,2009,20(3):398-402.)
[9] 周海燕.基于集合卡尔曼滤波法的非高斯含水层参数识别.博士论文.北京:中国地质大学,2011.(Zhou H Y.Charactering non­Gaussian 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:Pilot­point 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 pilot­points 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.Pilot­point 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.Pilot­point 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.Three­dimensional interference test interpretation in a fractured aquifer using the pilot­point inverse method.Water Resources Research,2001,37(11):2659-2675.
[17] Doherty J.Groundwater model calibration using pilot­points 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 ill­posed 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,Model­independent 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.MODFLOW­2000,the U.S.geological survey modular ground­water model—User guide to modularization concepts and the ground­water 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.
No related articles found!
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!