原创力文档- 1、本文档共21页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
地下水资源评价期末作业精选
地下水资源评价大作业—有限差分法
一、用均衡和达西公式建立河间地块地下水平衡方程，并求解
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.15m。
建立概化模型，将该河间地块按照流线边界划分为五个小模块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.12
文档评论（0）