地下水三维流多边形有限差分模拟软件开发研究及实例应用
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
地下水流数值模拟技术已成为地下水资源和地质环境评价问题分析的主要方法之一,而三维地下水流数值模拟软件本身存在两个主要的内容,一是可视化的前后处理,二是核心程序的仿真性。通用的地下水流模拟软件的形成已成为发展的主流趋势,它不仅要求能有效快速处理数据,而且强调核心程序对水文地质条件刻画的仿真性以及核心程序与前后处理程序良好的接口。目前国际上已推出多款地下水流模拟商业软件,而其对若干水文地质模拟要素的仿真技术有待改进,同时国内没有形成统一的模拟软件。本文的研究正是基于这个问题展开的。
    本文以地下水三维流多边形有限差分模拟理论为研究中心,采用了多学科技术综合应用研究的新技术和新方法,遵循着面向对象程序设计和系统设计的理论,在数字化和可视化思想的指导原则下,对三维模拟软件开发中的各个关键问题进行了探讨与改进,并且以河西走廊黑河干流中游地区为例,应用设计的三维模拟软件实际完成该区的三维地下水数值模拟研究。本文的主要研究内容、方法及工作有如下几点:
    1.全文围绕着三维数值模拟的前后处理及模拟系统集成的内容,对模拟中的各个关键问题,包括平面三角网络的剖分、平面等值线的形成、空间插值的实现、模拟数值科学计算可视化、数值模拟的各个水文地质要素处理问题,分别展开了较详细的论述及研究。
    2. 正确刻画水文地质模拟要素是水流数值模拟的基本出发点,因此本文详细介绍了地下水三维流任意多边形有限差分法数值模拟的基本理论,并对如何模拟经常遇到的大气降雨入渗的滞后性、潜水蒸发排泄、河流和水库及包括混合井在内的各种井孔、泉等问题作了较充分的分析和讨论。
    3.基于三维多边形有限差分地下水数值模拟的理论,对数值模拟系统进行了总体框架设计,采用面向对象语言C++中的类表示了不同的可视化研究内容,同时设计了模拟中的文件系统,认为在模拟系统中应首先把握住地下水的含水系统结构,整个系统必须有一个良好的总体设计和文件组织。
    4.在重点分析了前人平面三角剖分中波前推进法(advancing front method)和Delaunay 三角剖分方法的基础上,对Delaunay三角剖分方法中结点位置的设置作了改进,并实现了平面三角网络的剖分程序的编制,同时对网络中部分存在的钝角提出了网络编辑及优化的解决办法,其结果基本可满足水文地质数值模拟的要求。另外本文提出了平面三角网络的逻辑结构,并且在这个逻辑结构基础上很容易得出网络的侧边结点序号等信息。
    5.总结和介绍了目前存在的各种不同的离散化数据插值方法,以一个实际带有一定数
    
    
    量的采样点的函数分别对不同的插值方法进行了测试,分析认为MQS(改进的距离反比法)、径向基函数法(包括Multiquadirc和薄板样条法)、克里格法插值结果比较理想,而且插值结果同采样点的数量及分布有直接的关系。
    6.对数值模拟科学计算可视化的要素从一维、二维、三维数据场分别作了详细的讨论,其中简要分析了平面等值线绘制算法中经常用到的网络序列法和三角网格等值线法,同时编制了相应的程序,实现结果同国内外广泛应用的Surfer软件生成的结果一致;而且综合分析了已有的三维建模的各种开发环境及建模采用的方法,根据水流三维模型重建的要求及OpenGL本身方便独立的特点,详细介绍了利用OpenGL构建三维水流实体模型的“六结点单元法”及“八结点单元法”,提出的方法不仅简单,而且显示起来速度较快,可方便利用水流数值模拟前处理已生成的网络数据;同时从模拟计算中水文地质参数的意义出发,设计了地下水水流线的绘制算法,并予以实现。
    7.一个模型的好坏与水文地质模型的合理概化、水文地质模拟要素的仿真技术有着直接的关系,而相应的三维水流模型基础资料则是运行模拟必备的条件之一,因此本文详细归纳了建立三维水流模型所需的资料,这将为实际水文地质勘探工作提供参考和依据。在资料的获取基础之上,同时分析了数值模拟前后处理、核心程序及两者接口实现的内容,并以图件展示了前后处理实现的图形用户界面,探讨了不同格式的图形数据格式的定义及如何输入、输出及转化这些图形数据。
    8.以河西走廊黑河干流中游地区地下水三维流数值模拟为例,在充分分析研究区水文地质条件的基础之上,对水文地质参数进行了分区,进行了模型识别工作,从不同观测孔的拟合、泉的拟合及地下水水资源的均衡分析得出模拟结果可满足目前的水文地质情况,特别探讨和研究了地下水与地表水的转化关系,再以现状的开采动态进行了15年后地下水动态进行了预测,并对结果进行了分析。
Groundwater numerical simulation has been one main method for groundwater resources analysis and geology environmental evaluation. Three-dimensional groundwater modeling software has two basic research contents. One is pre-post processing with the technique of the visualization, and the other is model accuracy. Common groundwater modeling software is final trend of software development, and it not only requires effective and quick data process, but also the model accuracy of depicting hydro-geological condition, and the well-designed interface between core program and pro-post processing program.
    At present, several kinds of abroad commercial groundwater modeling softwares have turned on. However, many problems in depicting hydro-geological conditions of these softwares have to be modified. In the homeland uniformed modeling software hasn’t been formed. This paper is from this.
    Based on the theory of three-dimensional finite different groundwater modeling, making full use of new technique and methods in multi-subjects synthetical researches, importing the thought of object-oriented programming and software system design, and under the direction of the theory of digitization and visualization, this paper discussed and supplied the treatments for common key problems of modeling software development, and finally made one case study on middle reach of mainstream of the Heihe River in Hexi corridor. This case applied three-dimensional groundwater modeling programs designed by the theories in this paper. The main contents, methods and works are followed:
    1. Centered on pre-post processing and system integration of three-dimensional groundwater modeling, this paper makes researches on such key problems as automatic triangle mesh generation in arbitrary plane domain, plain contour map generation, and methods for space interpolations, scientific calculation visualization and how to depict the problems of many special hydro-geological conditions.
    2. Depicting hydro-geology elements correctly is the basic element of numerical simulation. So this paper introduced basic theory of three-dimensional finite different groundwater modeling, and discussed how to depict such general hydro-geology
    
    
    problems as the lag properties of infiltration,evaporation of phreatic water, river, reservoir, all kinds of wells including mixing well, springs and etc.
    3. From the view of the theory of groundwater flow modeling, this paper design the whole structure of modeling system, applied the class of C++ to express different objects of visualization theory, designed file system of modeling system. The author considers well-organized modeling structure will be the first-place content of the whole system.
    4. Comparing and analyzing advancing front method and Delaunay method in automatic triangle mesh generation of arbitrary plane domain, this paper made some improvements on the placement of node, and realized the program of the plain triangle mesh generation. This article also supplied the function for editing the obtuse triangle, and the results satisfied the requirement of modeling program successfully. Effective logic structure of triangle mesh is supplied, from which we can get adjacent side and node information easily.
    5. This paper summered current method of scatter data interpolations, and selected one factual function and compared the results from different inserting methods. By analysis and comparison, the author thinks the results by insertion of MQS, Radical Basis Function or Kriging method are perfect, and the results have direct relation with the numbers and distribution of sampling points.
    6. This paper discussed scientific calculation visualization of one, two, three-dimension data field of in groundwater numerical simulation, and analyzed the methods of the grid-sequence method and triangle mesh contour map generation, and at the same time programmed related programs, which are consistent with the results by surfer software. Also this article summered the methods of three-dimension solid model construction, and according
引文
Chunmiao Zheng; P Patrick Wang, A field demonstration of the simulation optimization approach for remediation system design ,Ground Water; May/Jun 2002
    Osiensky, James L; Williams, Roy E, A two-dimensional MODFLOW numerical approximation of mise-a-la-masse electrical flow through porous media ,Ground Water; Jul/Aug 1996
    Russell L Detwiler; Steffen Mehl; Harihar Rajaram; Wendy W Cheung, Comparison of an algebraic multigrid algorithm to two iterative solvers used for modeling ground water flow and transport ,Ground Water; May/Jun 2002
    Yongqiang Lan; James L Osiensky, Discussion of papers: Potential inaccuracies in MODFLOW simulations involving the SIP and SSOR methods for matrix solution ,Ground Water; Jan/Feb 1998
    Amar C Bumb; Judson T Mitchell; Seth K Gifford, Design of a ground-water extraction/reinjection system at a superfund site using MODFLOW ,Ground Water; May/Jun 1997
    Liang-Cheng Chang; Chin-Tsai Hsiao, Dynamic optimal ground water remediation including fixed and operation costs ,Ground Water, Sep/Oct 2002
    Keith J Halford,Effects of unsaturated zone on aquifer test analysis in a shallow-aquifer system ,Ground Water, May/Jun 1997
    Sophocleous, Marios; Koussis, Antonis; Martin, J L; Perkins, S P, Evaluation of simplified stream-aquifer depletion models for water rights administration ,Ground Water, Jul 1995
    Samuel P Perkins; Marios Sophocleous, Development of a comprehensive watershed model applied to study stream yield under drought conditions ,Ground Water, May/Jun 1999
    Norman K Jones; Trevor J Budge; Alan M Lemon; Alan K Zundel, Generating MODFLOW grids from boundary representation solid models ,Ground Water, Mar/Apr 2002
    Michael A Robinson; William G Reay, Ground water flow analysis of a Mid-Atlantic outer coastal plain watershed, Virginia, U.S.A. ,Ground Water, Mar/Apr 2002
    Richard B Winston, Ground water modeling software on the Internet ,Ground Water, Jul/Aug 2002
    Reynolds, Jeffrey W; Spruill, Richard K, Ground-water flow simulation for management of a regulated aquifer system: A case study in the North Carolina coastal plain ,Ground Water, Sep 1995
    Richard M Yager; Christopher J Neville, GFLOW 2000: An analytical element ground water flow modeling system ,Ground Water, Nov/Dec 2002
    
    Gary Miner, Groundwater Modeling Using Geographical Information Systems , American Water Works Association. Journal, Nov 2002
    Richard E Gerber; Ken Howard, Hydrogeology of the Oak Ridges Moraine aquifer system: Implications for protection and management from the Duffins Creek watershed ,Canadian Journal of Earth Sciences; Ottawa; Sep 2002
    Swain, Eric D, Implementation and use of direct-flow connections in a coupled ground-water and surface-water model ,Ground Water, Jan 1994
    Beverly L Herzog; David R Larson; Curtis Abert; Steven D Wilson; George S Roadcap, Hydrostratigraphic modeling of a complex, glacial-drift aquifer system for importation into MODFLOW ,Ground Water, Jan/Feb 2003
    John Doherty,Improved calculations for dewatered cells in MODFLOW ,Ground Water, Nov/Dec 2001
    William C Burton; L Niel Plummer; Eurybiades Busenberg; Bruce D Lindsey; William L Gburek, Influence of fracture anisotropy on ground water ages and chemistry, Valley and Ridge Province, Pennsylvania ,Ground Water, May/Jun 2002
    Michael P Kladias; Gregory J Ruskauff, Implementing spatially variable anisotropy in MODFLOW ,Ground Water, Mar/Apr 1997
    Warren Barrash; Martin E Dougherty, Modeling axially symmetric and nonsymmetric flow to a well with MODFLOW, and application to Goddard 2 well test, Boise, Idaho ,Ground Water, Jul/Aug 1997
    Michael A Jones, NCF: A finite-element computer program to simulate ground-water flow within the U.S. Geological Survey Modular Ground-Water Flow Model (MODFLOW) ,Ground Water, Jul/Aug 1997
    Daniel W Matthews; Joel Massmann; Stuart E Strand, Influence of aquifer properties on phytoremediation efffectiveness ,Ground Water, Jan/Feb 2003
    Rumbaugh, James O III ,New books -- A Manual of Instructional Problems for the U.S.G.S. MODFLOW Model by the U.S. Environmental Protection Agency ,Ground Water, Sep 1994
    M Wang; C Zheng, Optimal remediation policy selection under general conditions ,Ground Water, Sep/Oct 1997
    John Robert Hoaglund III; Gary Cecil Huffman; Norman Guy Grannemann, Michigan basin regional ground water flow discharge to three Great Lakes ,Ground Water, Jul/Aug 2002
    Mikko I Jyrkama; Jon F Sykes; Stefano D Normani, Recharge estimation for transient ground water modeling ,Ground Water, Nov/Dec 2002
    Jon P Fenske; S A Leake; David E Prudic, Simulating reservoir leakage in ground-water models ,Ground Water, Sep/Oct 1997
    David M Diodato, Software spotlight ,Ground Water, Jan/Feb 1998
    Adam D Festger; Gary R Walter, The capture efficiency map: The capture zone under time-varying flow ,Ground Water, Nov/Dec 2002
    
    John D Bredenhoeft, The water budget myth revisited: Why hydrogeologists model ,Ground Water, Jul/Aug 2002
    James L Osiensky; Roy E Williams, Potential inaccuracies in MODFLOW simulations involving the SIP and SSOR methods for matrix solution ,Ground Water, Mar/Apr 1997
    Hakan Basagaoglu; Timothy R Ginn; Christopher T Green; Benjamin J McCoy, Transport in heterogeneous media: Tracer dynamics in complex flow networks ,American Institute of Chemical Engineers. AIChE Journal; New York; May 2002
    http://www.visual-modflow.com/
    http://www.scisoftware.com/
    郭卫星,卢国平编译,朱学愚,刘建立校对,M.G.McDonald and A.W.Harbaugh原著,MODFLOW三维有限差分地下水流模型,南京大学地球科学系,1998
    董壮,三维水流数值模拟研究进展, 水利水运工程学报, 2002.9,第3期
    吴剑锋,朱学愚,由MODFLOW浅谈地下水流数值模拟软件的发展趋势,工程勘察,2000年第2期
    孙永堂,朱燕丽,空间插值方法及其在地下水数值模拟中的应用,吉林水利,2000年3月,总第207期
    张丽琼,崔广柏,杨珏,通量向量分裂格式及有限体积法在水流模拟中的应用,水利水电技术,2001年第8期,第32卷
    陈崇希,岩溶管道-裂隙-空隙三重空隙介质地下水流模型及模拟方法研究(国家自然科学基金资助),地球科学,1995,20(4),361-366
    陈崇希、裴顺平, 地下水开采-地面沉降数值模拟及防治对策研究[M] (国家“九五”重点科技攻关项目95-913-06-04). 武汉:中国地质大学出版社,2001
    陈崇希,滞后补给权系数-----降雨补给潜水滞后性处理方法,水文地质工程地质,1998,(6), 22-24
    McDonald M.G.and A.W.Harbaugh, A MODULAR THREE-DIMENSIONAL FINITE-DIFFERENCE GROUND-WATER FLOW MODEL, U.S. Geological Survey Techniques of Water Resources Investigations Book 6. United States Government Printing Office, Washington, D.C Chapter A1,1988.
    陈崇希、林敏等, 地下水混合井流理论及其应用(国家自然科学基金和地质行业科技发展基金资助),中国地质大学出版社[M].武汉:1998
    陈崇希、焦赳赳, Numerical simulation of pumping test in multilayer wells with non-Darcian flow in the wellbore , GROUND WATER ,1999,37(3),465-474
    陈崇希、唐仲华, 地下水流动问题数值方法, 中国地质大学出版社(武汉),1990年
    ftp:// 166.111.49.137资源:石教英-可视化技术正文.pdf
    网页http://www.cs.berkeley.edu/~jrs/
    孙文理,杨文硕,自动剖分生成有限元网格的算法研究,哈尔滨理工大学学报,2000,5(3):4-7
    程世才,林国成,王永洁,二维有限元网格全自动生成方法-AFM法,哈尔滨工业大
    
    
    学学报,1995,27(4):103-108
    郑伟,蔡永昌,张湘伟,多连通域三角剖分的网格推进法,汕头大学学报(自然科学版),2000,15(1):39-44
    丁永祥,约束Delaunay三角剖分与有限元网格自动生成,华中理工大学学报,1995,23(6):39-43
    闵卫东,唐泽圣,二维Delaunay三角划分的平均形态比最大性质,计算机学报,1994,17:20-25
    李国敏,陈崇希,克立格法在地下水数值模型中的应用,地质科技情报,1995,14(3):95-99
    胡立堂,周爱国,工程中三角网格等值线的实现与应用,工程地质计算机应用,2002,26(2):14-16
    唐泽圣等著,清华大学学术专著:三维数据场可视化, 1999年12月第1版
    胡友元,黄杏元主编,计算机地图制图,测绘出版社,1987年12月第一版
    张新长,曾广鸿,张青年,城市地理信息系统,科学出版社,2001年9月第一版
    彭群生,鲍虎军,金小刚编著,计算机真实感图形的算法基础,科学出版社,1999年6月第1版
    陈建春,Visual C++开发GIS系统—开发实例剖析,电子工业出版社,2000年9月第1版
    中地软件丛书,MapGis地理信息系统使用手册—空间分析篇,2001.5
    中地软件丛书,MapGis地理信息系统使用手册—数字制图篇,2001.5
    David J.Kruglinski著,潘爱民,王国印译,Visual C++技术内幕,第四版,MicroSoft Press,1999.4月第2次印刷
    廖朵朵,张华军,OpenGL三维图形设计,星球地图出版社,1996
    汪同庆主编,Fortran90程序设计,武汉大学出版社,2002.7
    曲焕林主编,中国干旱干旱干旱地区地下水资源评价,科学出版社,1991.7
    康凤举主编,现代仿真技术与应用,国防工业出版社,2001.9
    供水水文地质手册,地质出版社,1977
    张秀山等编著,虚拟现实技术及编程技巧,国防科技大学出版社,1999.9
    朱学愚,钱孝星,刘新仁编著,地下水资源评价,南京大学出版社,1987.11
    黄河水利委员会黑河流域管理局,黑河流域水资源评价利用及分配方案研究,1999.3
    甘肃省地质矿产局第二水文地质工程地质队,甘肃省黑河干流中游地区地下水资源及其合理开发利用勘察研究,1990.3
    姚足金,陈德华,姚洪,干旱区高山-盆地系统深部水资源的赋存和估算--以河西走廊西段为例,地球科学,2001.5,第26卷第3期
    徐中民,黑河流域水资源管理方式的初步探讨,地域研究与开发,2000.3
    张济世,康尔泗等,黑河中游水土资源开发利用现状及水资源生态环境安全分析,地球科学进展,2003.4
    龙爱华,徐中民等,基于边际效益的水资源空间动态优化配置研究,冰川冻土,2002.8
    张勃,李吉均,河西地区黑河流域水资源空间组合与优化利用研究,盐湖研究,2001.3
    许平,黑河流域中游灌区节水对策与措施,中国水利,2002.8
    
    郝庆凡,黑河调水及思考,人民黄河,2001.7
    曲耀光,樊胜岳,黑河流域水资源承载力分析计算与对策,中国沙漠,2000.3
    王宝玉,安新代,黑河流域水资源合理利用探讨,西北水资源与水工程,1998.4
    黄高宝,柴强,柴守玺,张掖市水资源利用现状与可持续开发,干旱地区农业研究,2001.9
    王根绪,黑河流域额济纳绿洲区水资源合理利用分析,兰州大学学报(自然科学版),1997.33
    曹琪,李开月,深化水利管理改革—实现黑河流域综合治理目标,甘肃水利水电技术,2003.9
    藏忠淑,邢俊兵,MapGis与ArcView之间的文件转换技巧,黄金地质,2003年3月,Vol.9,No.1
    卢文喜,地下水运动数值模拟过程中边界条件问题探讨,水利学报,2003年3月,第3期
    陈崇希,"防止模拟失真,提高仿真性"是数值模拟的核心,水文地质工程地质,2003年第2期
    罗振东,非饱和水流问题的混合元法及其数值模拟,计算数学,2003年2月,Vol.25,No.1
    潘世兵,朱雪芹,王忠静,河西走廊地下水系统数值模拟中的几个问题探讨-以酒泉盆地为例,工程勘察,2003年6月
    周卫,利用ArcView实现地面沉降数据的处理和可视化,南京建筑工程学院学报,2001年第3期,N0.59
    尹大凯,胡和平,惠士博,宁夏银北灌区井渠结合灌溉三维数值模拟与分析,灌溉排水学报,2003年2月,
    杜绍敏,宁冰,水文地质数值模拟的立足点,东北林业大学学报,2003年5月,Vol.31,NO.3
    王玮,水文地质数值模拟中节点地面标高的获取方法,长安大学学报(地球科学版),2003年6月,Vol.25,No.2
    王晓梅,朱国荣,余勤等,苏锡常地区主采层地下水流数值模拟,水文地质工程地质,2003年第1期
    郑文刚,离散元工程计算软件前后处理系统,大连理工大学博士学位论文,2002.5
    李硕,GIS和遥感辅助下流域模拟的空间离散化与参数化研究与应用,南京师范大学博士学位论文,2002.5
    秦绪佳,医学图学三维重建及可视化技术研究,大连理工大学博士学位论文,2001.6
    万华根,基于虚拟现实技术的CAD方法研究,浙江大学应用数学专业博士学位论文,1999.10
    阎丽霞,虚拟手术关键技术研究,浙江大学计算机应用技术博士学位论文,2001,2
    王文成,三维数据场体绘制技术的研究与实现,中国科学院软件所博士学位论文,1998.1
    马小虎,虚拟现实中多细节层次模型的研究,浙江大学博士学位论文,1997
    李桂清,细分曲面造型及应用,中国科学院博士学位研究生学位论文,2001.5
    刘慎权,CAD与图形学学术年会二十年回顾,计算机辅助设计与图形学学报,1999.5
    陈永华,WEBGIS三维可视化的研究,中国人民解放军信息工程大学测绘学院博士学
    
    
    位论文,2000.4
    刘贤喜,机械系统虚拟样机仿真软件的应用化研究,中国农业大学车辆工程专业博士学位论文,2001.10
    李智慧,基于可视化的三维放射治疗计算机模拟系统的研究,四川大学固体力学专业博士学位论文,2002.4
    李维诗,基于医学断层轮廓数据的反求CAD建模理论与方法研究,浙江大学机械制造及其自动化专业博士学位论文,2001.8
    费广正,三维复杂模型快速简化及交互方法的研究,中国科学院博士学位研究生学位论文,2001.1
    孙文武,三维数据场体绘制中几类问题的研究与实现,中国科学院博士学位研究生学位论文,2001.2
    杨晓松,三维有限元数据场可视化技术研究与实现,大连理工大学工程力学系博士学位论文,2000.3
    张文,矢量场可视化算法研究与系统设计,国防科技大学计算机科学与技术专业博士学位论文,2001.10
    吴良武,体素建模、绘制及其在医学图象中的应用,大连理工大学机械制造及自动化专业博士学位论文,2001.7
    汪惠芬,协同开发环境中CAD/CAM集成技术研究,南京理工大学控制理论与控制工程专业博士学位论文,2001.10
    吴镇,优化设计中的智能方法与可视化技术研究,武汉理工大学博士学位论文,2001.7
    申杰,水利土木工程领域中的扩充现实技术,四川大学博士学位论文,2001.5
    郭洪宇,区域水资源评价模型技术及其应用研究,中国农业大学博士学位论文,2001.6
    王惠中,浅海与湖泊三维环流与水质数值模拟研究和应用,河海大学博士学位论文,2001.7
    胡云进,裂隙非饱和渗流试验研究及有地表入渗的裂隙岩体渗流数值分析,河海大学博士学位论文,2001.7
    张新民,干旱区水资源量与质统一管理研究,西安理工大学博士学位论文,2000.10
    管群,VR-GIS技术在岩土工程中的应用,四川大学博士学位论文,2001.12

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700