第六章(1)_地下水数值模拟模型应用实例xiugai
地下水数值模拟
对海水入侵问题的认识摘要：随着地下水的超采，我国许多沿海城市都出现了海水入侵的现象，严重影响到了自然生态环境及人们的生活，本文介绍了关于海水入侵的一些概念，如何运用数值模拟的方法对该现象进行预测分析，以及据此得到的防治对策。
关键字：海水入侵，数值模拟，防治对策一、海水入侵的概念及现状海水入侵是由于滨海地区水动力条件发生了变化，引起高盐度的海水或者咸水向淡水含水层运动，从而发生水体入侵，造成水质恶化，土壤次生盐渍化的现象。
海水入侵现象与岩相及构造、地理环境、气象条件，以及人类的活动密不可分。
自然条件固然是造成这一现象的基本因素，但是在许多地区，人类对地下水的过量开采是最直接的诱因。
沿海地区，海水与淡水之间存在一个咸淡水接触面，当淡水位高于海水位时，淡水向海中流动，淡水不会被咸化。
但是人类过量开采地下水，直接导致地下水位下降，水压降低，海水向淡水含水层进行入侵。
海水入侵在中国沿海城市是一个非常普遍且严重的问题，由于经济的发展和人类生活的需求，许多沿海城市都存在地下水超采的问题，地下水的补给速度远小于开采速度，水位不断下降，海水入侵问题也愈演愈烈。
海水入侵对自然环境和人类生活都造成了非常巨大的危害，海水入侵使得地下水盐分增加，在入侵区域的工业设备容易受到腐蚀，使用寿命减少，产品质量降低，因水质处理或远距离调水增加了生产成本，甚至导致工厂不堪重负，搬迁或者倒闭，减缓了经济的发展；长期用咸水灌溉，会导致土壤发生盐渍化甚至板结，肥力下降，粮食减产；海水入侵还会导致居民的饮用水咸度增加，水质变差，包含更多对人体有害的元素，对人类的身体健康造成危害，许多淡水源地更因此受到影响，严重者甚至被废弃，大大干扰了人们正常的生活饮水，造成用水紧缺，降低了人们的生活质量。
20世纪60年代初期，中国的大部分沿海城市还处于天然的水动力条件下，很少发生海水入侵的现象，水质良好。
20世纪70年代开始，由于国家经济的发展，地下水开采量渐渐增大，咸淡水平衡遭到破坏，海水入侵的现象变得越来越普遍，辽宁、河北、天津、山东、江苏、上海、浙江、海南、广西9个省份的沿海地区都有不同程度的海水入侵现象发生，其中最严重的是山东、辽宁两省，入侵总面积已超过2000平方公里。
第三讲 地下水数值模拟原理及建模方法和步骤_xiugai
矩形网格 多边形网格
方法一：差商代替微商
1、网格划分的基本类型
（1）先划格线，格点位
均
于网格中心
衡
网
格
（2）先规定格点位置，
再垂直平分两相邻结点的连
节
线作格线，形成的网格即为 水均衡区
H (x0
x, t0 )
2H (x0 , t0 ) (x)2
H (x0
x, t0 )
一维控制方程差分格式
方法一
控制方程
T
2h( x, t ) x2
h( x, t ) t
网格剖分nx个
显式差分格式
H (x0 x, t0 ) 2H (x0 , t0 ) H (x0 x, t0 ) H (x0 , t0 t) H (x0 ,t0 )
(x) 2
t
Hn i 1
2
H
n i
Hn i 1
(x)2
H
n1 i
H
n i
t
i 2,3,, nx 1
隐式差分格式
H (x0 x, t0 t) 2H (x0 , t0 t) H (x0 x, t0 t) H (x0 , t0 t) H (x0 , t0 )
(x) 2
t
H n1 iБайду номын сангаас1
数值方法很多，但是最简单实用的是有限差分法：
✓ 有限差分法
✓ 有限单元法
✓ 积分有限差分法
✓ 半解析半数值法
✓ 边界元法
✓ 有限体积法
只讲有限差分法
一、有限差分法的基本原理
第五章 地下水流模型MODFLOW简介xiugai(1)
2.地下水数值模拟发展趋势
（5）尚未建立地下水动态模拟体系
近几年来，随着GIS技术的发展和普及，GIS已被广泛接受为管理、 存储、查询分析有关地理空间分布信息的成熟的工具和技术。随着其 本身空间分析功能和相关学科的迅猛发展，应用领域不在拓宽。随着 信息化成为当今社会发展的潮流和趋势，地下水动态模拟评价信息化 已成为水文地质工作的工作方向。如何充分利用GIS在数据采集、存 储、处理和可视化表达等地下水模拟重要环节上的优势，把GIS与地
2.模拟软件发展趋势
(2)遥感、地理信息系统和Leabharlann 球定位系统与地 下水数值模拟软件的耦合
遥感、地理信息系统和全球定位系统在地下水数 值模拟中的进一步应用，地下水模型需要一个专业的G 璐数据结构，当前国际上的主流地下水模拟软件都提 供了与GIS(如ARCYm2FO)的数据接口，随着GIS 在水文地质领域中的应用日益广泛和深入，它们之间 的联系会更加密切。GIS与地下水模拟模型紧密集成 技术将是将来的主要研究方向。
Visual MODFLOW主要处理模块：
Visual MODFLOW 界面设计包括三大彼此 联系但又相当独立的模块:
a)输入模块 b)运行模块
c)输出模块
3.1 Visual MODFLOW
输入模块
允许用户直接在计算机上赋值所有必要的输入参数以便自动生成 一个新的三维渗流模型。输入菜单把MODFLOW 、MODPATH 和MT3D 的数 据输入作为一个基本建模块,这些菜单以逻辑顺序排列并显示, 指导用 户逐步完成建模和数据输入工作。允许用户直接在计算机上定义和剖 分模拟区域, 用户可随意增减剖分网格和模拟层数, 确定边界几何形 态和边界性质,定义抽(排) 水井的空间位置和出水层位以及非稳定抽 排水量。参数菜单允许用户直接圈定各个水文地质参数的分区范围并 赋值相应参数,同时上、下层所有参数可相互拷贝。可预先定义水位校 正观测孔的具体空间位置和观测层位, 并输入其观测数据, 以便在后 续的模型识别工作中模拟使用。
地下水数值模拟01
基础知识
专业基础
水文学原理 水文地质学基础 地下水动力学 工程地质学 环境地质学 高等数学 线性代数 数学物理方法 概率论
数学基础
计算机基础
编程语言 计算机绘图 应用软件
数学模型
•利用模型与原型之间在数学上 的相似为基础 •例：达西定律、溶质运移方程
基本概念2
地下水数学模型的分类
确定性模型与随机模型
根据模型中变量的取值性质划分 •确定性模型：变量取确定值 •随机模型：含有随机变量，只知道 取值概率
集中参数模型与分布参数模型
线性模型与非线性模型
………
根据模型中是否含有空间变量划分 •集中参数模型：流量与降深经验公式 •分布参数模型：含有偏微分方程，裘布 依模型、泰斯模型
水文地质计算的数值方法 孙峰根 中国矿业大学出版社，1995
授课内容
• • • • • • 第一章 第二章 第三章 第四章 第五章 第六章 绪论 地下水运动的数学模型 有限差分法 有限单元法 其他数值方法 反求参数的数值方法
第一章 绪论
主要内容
地下水数值模拟及其基本任务 常用的数值方法简介 数值模拟的基本过程
地下水数值模拟
2012.5
参考教材
地下水数值模拟 薛禹群，谢春红 科学出版社，2007 水文地质学的数值法 薛禹群，谢春红 煤炭工业出版社，1980
地下水流动问题数值方法 陈崇希，唐仲华 中国地质大学出版社，1990
地下水流数值模拟 李俊亭 地质出版社，1989 地下水流的数学模型 和数值方法 孙纳正 地质出版社，1981
根据模型中变量的阶次划分 •线性模型：承压水运动 •非线性模型：潜水运动
基本概念3
数学模型的求解方法
解析法
地下水数值模型在工程上的应用
地下水数值模型在工程上的应用【摘要】本文以某水利枢纽工程库区为研究对象，在前期水文地质测绘、勘探工作的基础上，对上述均衡计算区建立地下水数学模型，以定量评价羌河一级水电站、二级水电站及三级水电站建成运行以后，对规划河段河谷地下水环境的影响，进行定量或定性的分析，并做出预测评价。
【关键词】水文地质；数值模型；水资源量1、含水层结构概化评价区地处羌河谷中，与周边中低山区对照比较，地形地貌及堆积物有明显的区别，可综合归属为一个独立的水文地质单元。
区内含水层在河床中自上到下为第四系全新统-上更新统冲积卵砾石层；在河谷两侧分布冲洪积扇，由上到下依次为第四系全新统冲洪积碎石土层、第四系全新统冲积卵砾石层及第四系上更新统冲积卵砾石层。
含水层的岩性较不均一，为非均质含水层，其非均质性用含水层参数（导水系数T、给水度μ）分区概化处理。
根据勘探试验获取的参数系列值做成含水层参数分区图，给出各区的参数均值作为数值计算的初值，经模型调试和识别，最终将试验参数系列转化为模型参数系统。
由于浅层与深部含水层之间存在较好的水力联系，所以将计算区含水层概化为与地表水有密切水力联系的单层各向同性平面二维渗流的潜水含水层。
根据模型范围内地质、水文地质条件分析，现有条件下，模拟区主要接受侧向补给（包括基岩裂隙水和河谷潜流补给量）和田间灌溉水入渗补给，以及少量降水入渗补给，对概念模型可概化。
2、地下水模型中数据选用在1∶1000的比例尺的地形图上获取地面高程数据及羌河水位，用surfer软件进行对插值，得到地面等高线及地下水等水位线图，再利用MODFLOW自带的插值模块Field Interpolator，确定了模拟区各个单元的高程和河流单元水位。
根据地形图上地形和河水水位，以及现场调查的资料，绘出地下水初步流场图，作为稳定流模型的地下水初步水头。
各单元的河水水位确定主要依据尾部水文站断面实测流量—水位关系曲线及规划报告水文章节中相关断面的水位—流量关系曲线来模拟确定。
地下水数值模拟gms应用基础和实例
地下水数值模拟gms应用基础和实例
《地下水数值模拟: GMS应用基础与实例》这本书介绍了地下水数值模拟的基本原理和应用。
以下是关于这本书的详细信息：
首先，这本书介绍了地下水运动的基本方程，包括水流运动方程、非水相流体（NAPL）运动方程和地下水中化学组分迁移方程。
其中，非水溶性液体运动方程与GMS采用的方程一致。
其次，这本书以实例为基础，详细介绍了在GMS环境下进行地下水数值模拟的全过程。
这个过程包括水文地质模型概化、几何模型建立、边界条件设定、参数赋值和反演，以及计算模拟与结果分析。
通过这本书的学习，读者可以较为深入地掌握地下水数值模拟的技术方法。
此外，这本书还介绍了地下水数值模拟中常用的模块，如TIN、Solids、Modflow2000/2005、MT3DMS、MODPATH、PEST、SEAWAT等，并详细阐述了这些模块在模拟地下水流动、地下水溶质运移、质点运移和海水入侵等过程中的应用。
最后，这本书强调了三维地质结构建模、水文地质模型概化、边界条件设定、参数反演和模型校核等关键环节在地下水数值模拟中的重要性。
此外，还介绍了如何处理GMS模型输出数据、编制相关图件以及如何进行模拟结果的三维可视化展示等内容。
总的来说，《地下水数值模拟: GMS应用基础与实例》是一本介绍地下水数值模拟
的综合性教材，适合环境、地矿、城建、水利等部门从事地下水研究的技术人员阅读，也可作为环境、地矿、城建、水利等专业的高年级本科生及研究生的教学用书。
第六章(3) 地下水数值模拟模型的应用实例xiugai
图5-51 距垃圾场中心约150m（cell Id：4227） 溶解氧浓度变化曲线
3、硫酸根在地下水中的浓度分布规律
图5-60 垃圾场中心（cell Id：4867）硫酸根浓度变化曲线
图5-63 距垃圾场中心约200m污染羽边缘（cell Id：4147） 硫酸根浓度变化曲线
在所有点中硫酸根都随季节有轻微的 波动，这可能是由于雨季降雨对地下 水中随酸根的稀释作用。 在填埋场中心，由于硫酸根是补给源， 硫酸根的浓度是持续上升的，并在 6000d左右时基本稳定。 随着距填埋场的距离增加，硫酸根上 升的速率逐渐降低。 在污染羽边缘（约200m）处，硫酸根 先是降低然后增加，这是因为一开始 该区域固有的硫酸根先与BTEX反应， 而垃圾渗滤液中的硫酸根还没有补给 到这一区域，所以硫酸根离子先降低， 随后硫酸根逐渐补给，故其浓度又开 始逐渐增加。 在污染羽范围之外，硫酸根离子则是 逐渐降低。
a：概念模型
场区为三面环山地形，与单斜地层基本构成一独立的水文地质单元。 填埋场地下方分布有新生界第四系及中生界侏罗系、古生界二叠系地层。 本区第四系地层分布不甚规律，缺少全新统。上更新统下段洪积层（莲花组 PLQ3）以黄色粘土、亚粘土组成，缺少砂砾层位。平均厚度约3米左右；中更 新统洪积层（之江组PLQ2）以棕红色粘土、亚粘土组成，结构紧密，层厚一 般为7米左右。 第四系下覆有埋藏型裂隙岩溶水，含水层由古生界二叠系下统栖霞组 （P3q）、石炭系上统船山组（C3c）、中统黄龙组（Czh）灰岩、白云质灰岩 组成。
（2）厌氧状态下在反硝化细菌作用下，NO3-作为电子接受体，接 受电子氧化BTEX。 （3）厌氧状态下在铁还原菌的作用下，Fe(Ⅲ)还原为Fe(Ⅱ)接受 电子，氧化BTEX。 （4）厌氧状态下在硫酸根还原菌的作用下，S6+还原为S2-接受电 子，氧化BTEX。
地下水流数值模拟的基本理论及应用
地下水流数值模拟的基本理论及应用王旭升博士目录1.地下水系统及其概念模型2.地下水流的数学描述与参数3.三维有限差分模型与MODFLOW4.模块及其作用5.VMODFLOW的应用23McWhorter and Sunada (1977)Q 1Q 3Q 2承压含水层地下水系统含水层概念模型：潜水含水层潜水面底板底板61. 地下水系统及其概念模型含水层概念模型：承压含水层顶板底板底板只有顺层流动测压水位面顶板含水层概念模型：多个含水层123模型范围含水层概念模型：多个含水层底板13425弱透水层：越流VMODFLOW4.1中文版培训 北京 2008年12月2. 地下水流的数学描述与参数* 承压水运动方程surface方 向 渗 透 系 数f lo w含水层 厚度y C MεhABεx补给强度 贮水系数∂H =0稳定流 ∂H ⎤ ∂ ⎡ ∂H ⎤ ∂ ⎡ ⎢ K xx M ∂x ⎥ + ∂y ⎢ K yy M ∂y ⎥ + ε = S ∂t ∂x ⎣ ⎦ ⎣ ⎦11VMODFLOW4.1中文版培训 北京 2008年12月2. 地下水流的数学描述与参数* 承压水运动方程 导水系数 Txx=KxxM Tyy=KyyMsurface f lo w贮水系数yεhAB C MS=SsMxε贮水率∂H ∂ ⎡ ∂H ⎤ ∂ ⎡ ∂H ⎤ ⎢Txx ∂x ⎥ + ∂y ⎢Tyy ∂y ⎥ + ε = S s M ∂t ∂x ⎣ ⎦ ⎦ ⎣12VMODFLOW4.1中文版培训 北京 2008年12月2. 地下水流的数学描述与参数* 三维渗流方程∂H ⎤ ∂ ⎡ ∂H ⎤ ∂ ⎡ ⎢ K xx ∂x ⎥ + ∂y ⎢ K yy ∂y ⎥ ∂x ⎣ ⎦ ⎣ ⎦ ∂ ⎡ ∂H ⎤ + ⎢ K zz ⎥ +W ∂z ⎣ ∂z ⎦ ∂H = Ss ∂t 体积源汇项h=500 500 700 1400 xh=0400 z贮水率13VMODFLOW4.1中文版培训 北京 2008年12月2. 地下水流的数学描述与参数定解条件 A h1 Q 边界条件 模型边界的水头、流量。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
6 257.5 302.5 272.9 236.5 240.5 268.6 244.6 249.3 241.5 270.3 256.5 305.4 262.2
7 213.80 255.1 219.70 193.4 218.1
211 234.7 202.7 204.22 227.4 211.8 228.60 218.4
0.1
10
0.08
8
0.08
9
0.05
5
0.06
7
0.02
3
河流入渗
（1）汾河入渗系数
（2）晋阳湖渗漏模数
按越流公式计算出湖水渗漏量，再除以其面积，则得晋阳 湖渗漏模数。
灌溉回渗补给地下水系数
受岩性、水位埋深土壤含水量、灌溉定额的影响。其值可近 似用次降雨量入渗补给地下水系数（α次）代替。代替的条 件是前十天以上时间无降雨，选取的次降雨量要接近灌溉定 额。本区稻田及菜地较多，灌溉量大，故取次降雨量为70～ 150毫米。
●把侧向补给量分配到相应 的层面上。
●在模拟的过程中有小部分 的流量边界，利用水均 衡原理做了适当的调节， 以便于后期的模拟。
参数分区及参数取值
模型参数有：
渗透系数及储水系数（分层（10）） 降雨入渗系数（列出公式、岩性） 河流入渗系数（沿河流分段，岩性、水位与河） 灌溉入渗系数（平面分出） 土面蒸发系数
10 144.72 129.2 141.40 93.7 108.7
109 114.2 95.7 99.09 107.8 105.1 114.37 113.6
11 63.73 76.3 66.70 50.9 62.2 82.3 78.1 62.5 61.13 75.3 71.3 62.15 67.7
12 42.54 49.7 42.90 30.6 34.8
地下水开采与地面沉降
图3-4 不同地点地下水开采量、水位降深与地面沉降对比曲线图
图3-3 太原市1956-2000年地面沉降等值线与2002年中深层水等水位线关系图
开始
确定计算区范围，剖分 网格，确定边界类型
整理输入基本数据：初 始水头，抽水井开采量， 降雨量，蒸发量及河流
数据等
第一次给定参数的初始估计值
图4-24 第六压缩层弹性释水因子
太原盆地观测井的拟合
观测孔分布图
水头高程（米）
790
780
770
760
750
740
730
720
710
700
Xa14
Ya9
D47
D10
Xa7
D57
D9
D13
D56
D36
87
Ja6
Pa29
井号
实测值 计算值
52个观测孔2003年12月16日计算值 与观测值的比较图
有限差分法计算井水位的校正
hw
hi, j,k
渗透系数
太原盆地参数分区
分区
地点
1
太原北郊汾河冲积扇
2
文峪河洪积扇
3
汾阳孝义冲积、洪积扇裙
4
龙凤河洪积扇及冲积平原
5
潇河洪积扇
6
昌源河洪积扇
7
太原西山倾斜平原
8
太原、清徐、交城西边山
9
汾阳、孝义倾斜平原
10
太原东山
11
榆次太谷祁县倾斜平原
12
祁县、平遥、介休倾斜平原
13
太原南郊汾河、潇河冲积平原
14
50 52.1 44 43.69 48.9 51.6 39.08 44.2
人工开采量
模型中：
工业及城市生活用水， 采用深井取水。
农村生活用水均采用浅 层面源取水。
乡村地区的农业用水， 按春灌、夏灌、冬 浇把取水量分配到 整个区域面积上。
模型识别校正
方法：采用间接法 原则及标准：要求
①观测孔水位动态趋势基本一致 ②地下水流场特征相同 ③水均衡计算相同 对比分析： ①观测孔水位动态对比结果 ②流场水位分布对比结果 ③水均衡对比结果
太原市地面沉降
太原市地面沉降范围为：北 纬37°40′—38°00′，东经 112°25′—112°45′；北起 上兰镇、南至西草寨；西起 金胜、东至武宿；包括太原 市辖区六个区（尖草坪区、 万柏林区、杏花岭区、迎泽 区、晋源区及小店区），面 积约585km2（图1- 140
19660
2003年12月第3层实测
2003年12月第3层计算
利用2003年8月和12月分别在研究区做过的两次水位统测.
水均衡分析对比
表 4-3 2004 年太原市孔隙水补给量（单位：亿米 3）
时间
降雨补给 侧向补给 汾河补给 渠道渗漏
春季 3-5 0.112
0.195
0.176
0.105
1 46.88 53.7 47.10 36.7 39.2 50.2 52.6 47.4 47.44 53.1 52.9 41.83 47.40
2 65.36 76.3 68.33 62.6 60.9 84.6 65.7 62.1 64.57 69.2
68 61.71 67.4
3 129.2 142.4 132.5 115 122 128.6 127.2 114.5 126.5 123.7 127.4 130.6 126.6
太原盆地结构模型
太原盆地实际钻孔位置分布图
太原盆地立体图
太原盆地结构模型
太原盆地孔隙介质三维模型
太原盆地孔隙介质结构三维可视化模型
太原盆地结构模型
盆地孔隙介质简化三维模型
盆地孔隙介质简化三维模型剖面图
二、太原盆地数值模型
太原盆 地数值 模型研 究区域 范围
模型范围： 周边界： 盆地与基岩山区接触带 上界：地表 下界：人为边界300m-350m深。 地下水流动特征： 边山岩溶水补给盆地孔隙水； 没有稳定隔水层（参考结构图） 盆地内部构成统一流动系统； 地下水混合开采； 因此属于：三维流系统。
Q
2 T
(In
x rw
)
2
H（m） 20.0
19.0
18.0 17.0 16.0
GMS（x=50） hw GMS（x=10）
15.0
14.0
13.0
1 2 3 4 5 6 7 8 9 10 T（d）
井水位校正对比图
图3-5（a）初始网格 图3-5（b）加密网格
地下水流场对比
拟合结果
1：流场趋 势基本一致
二、太原盆地数值模型
模型类型：三维流模型 构建思路： 1、太原盆地地层结构突出，不存在相对稳定的隔水层，难以划分含
水层； 2、三维建模将模拟区分为若干层（有人称为物理层），各层可以按
岩石含水渗透性确定其含水岩组、或隔水岩组分布； 3、边界的处理问题
底边界：基岩为不透水边界；第四系人为边界为流量边界，模型识 别时确定。
地下水数值模拟模型 的应用实例
地下水数值模拟模型的应用实例
7.1 区域地下水流模型及其应用-太原盆地 地下水数值模型
7.2 地下水水源地可开采资源量评价—宿 州市水源地数值模型
7.3 地下水污染预测研究—嘉兴垃圾填埋 场地下水溶质运移模型
7.1
太原盆地 地下水
数值模型
太原盆地地下水数值模型
4 216.10 251.9 237.40 199.5 212.3 225.6 231.7 199.4 202.80 206.0 208.8 240.10 219.3
5 272.49 309.6 294.50 249.1 240.7 282.2 251.9 280.5 245.47 264.9 251.7 327.10 272.5
划分调整各参数分区，确定各 参数上下限值
确定优选参数
数值模拟求逆问题
试
估
、
求目标函数G值，选出参数值
校
正
调
参数是否合理
整
是
参
数否
观测孔水位和地面
沉降量是否拟合好
是 最终确定参数
输入预测期有关数据进行预测
结束
图4-1地面沉降求解的计算流程图
图4-6 潜水水平传导系数分区图（K） 图4-7 潜水垂直传导系数分区图（Kz ）
我们依据1981-2000年晋 祠地震基准台跨断层I等短 水准高程实测资料（BM1BM3水准点），推算出晋 祠大断裂下降盘累计位移 量为21.88mm，平均为 1.15mm/a。这与吴家堡 沉降中心的96.18mm/a的 沉降速率相比，似乎说明
该地的构造运动与地面沉 降关系不大。
图3-2 太原地区构造图 1.活动断裂 2.地垒与凹陷边界
2：漏斗等 重点开采区 能反映实际 流场状况
3：在盆地 边缘的山地 还有某些区 域水头有出 入
4220
4200
4180
4160
4140
4120
4100
19560
19580
19600
19620
19640
19660
4220
4200
4180
4160
4140
4120
4100
19560
19580
19600
19620
为解决太原盆地出现的各种环境问题，1992年，研 究工作者们初步解决了太原市和太原盆地岩溶水与孔隙 水统一管理的问题，但由于当时的勘探目标是为农田和 城市供水，勘探技术、方法和指导思想相对落后，在勘 探过程中缺乏地下水系统指导思想和动态观念，大部分 勘察项目以行政区为单元，而不是以盆地地下水系统为 单元，导致边界不好控制、水量重复计算、忽略了相邻 地区的影响和边山地带裂隙水、岩溶水及地表水-地下水 相互作用的调查；同时由于缺乏现代先进的数据处理手 段，对于钻孔数据和水土样分析数据、地下水监测数据 等资料开发、利用程度很低，地下水系统三维空间结构 并没有弄清。