原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
地下水资源的最终评估
地下水资源评价大作业—有限差分法
一、用均衡和达西公式建立河间地块地下水平衡方程，并求解
1、用单元中心法建立5个计算网络（如图）。
2、建立单元流量均衡方程（达西公式）。
3、解单元流量均衡方程组（Excel，其中H0=11.15m，H6=10.15m），绘制迭代曲线。
４、求解地下水补给和排泄量，并分析误差产生的原因。
1， 图一 河间地块模型
2，如上图，为一河间地块及其概化模型。已知水文地质参数K＝10m/d；降雨量：445毫米/年；入渗系数：0.35；H0＝53.00－41.85＝11.15m；H6＝52.00－41.85＝10.15
建立概化模型，将该河间地块按照流线边界划分为五个小模块P1，P2，P3，P4，P5。然后对五个模块分别计算。
对P1单元格研究，建立水均衡方程：
左侧流入P1的水量 ：
对于P1来说，K为一定值，左侧过水断面为 ： ，
水力梯度： ，得：
根据达西定理：
右侧流出P1的水量 ：
对于P1来说，K为一定值，
右侧过水断面为： ，
水力梯度： ，
根据达西定律得：
则单元P1的均衡方程为：
（2）对P2单元格研究，建立水均衡方程：
左侧流入P2的水量 等于从P1右侧流出的水量 ：
右侧流出P2单元格的水量 ：
过水断面：
水力梯度：
根据达西定律得水量：
P2的均衡方程为：
（3）同理可得P3，P4，P5单元格的均衡方程如下：
P3均衡方程：
P4均衡方程：
P5均衡方程：
（4）确定相应的水温地质参数：K=10m/d，降雨量445mm/n，入渗系数：0.35；每天入渗量W=0.00043m/d；
H0＝53.00－41.85＝11.15m；H6＝52.00－41.85＝10.15m
将以上水文地质数据代入均衡方程式得到：
P1均衡方程：
P2均衡方程：
P3均衡方程：
P4均衡方程：
P5均衡方程：
改写以上五个均衡方程得：
P1：
P2：
P3：
P4：
P5：
计算上述五式的，得到水位迭代公式：
令计算个迭代初始水头为11.25，并代入以上方程，进行第一次迭代：
以下进行循环迭代（用EXCEL做）：可得如下数据。
表一 利用以上公式在EXCEL中作出的数据表格
h1
h2
h3
h4
h5
0
11.15
11.15
11.15
11.15
11.15
1
11.30161
11.42679
11.26958
11.29904
10.84576
2
11.38722
11.55574
11.47386
11.24542
10.90747
3
11.42789
11.70357
11.50104
11.3919
10.885
4
11.47498
11.73469
11.64724
11.39824
10.94759
5
11.48509
11.8352
11.66395
11.5104
10.95085
6
11.51744
11.84767
11.77023
11.52096
10.99933
7
11.5216
11.92019
11.78107
11.60373
11.00437
8
11.54509
11.92717
11.85878
11.61168
11.04041
9
11.54748
11.98023
11.86588
11.67249
11.04426
10
11.56474
11.98461
11.92289
11.67787
11.07087
11
11.56627
12.0236
11.92752
11.72257
11.07353
12
11.57899
12.0264
11.96942
11.72611
11.09316
13
11.57999
12.0551
11.97242
11.75901
11.09495
14
11.58938
12.05688
12.00326
11.7613
11.10943
15
11.59002
12.07804
12.00517
11.78554
11.11061
16
11.59696
12.07914
12.0279
11.78699
11.1213
17
11.59736
12.09475
12.02908
11.80487
11.12208
18
11.60249
12.09541
12.04585
11.80577
11.12997
19
11.60274
12.106
文档评论（0）