地下水资源评价实习报告

《地下水资源评价实习报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、中国地质大学（北京）地下水资源评价论文PMWIN数值分析实习班级：053082姓名：苏思慧学号：0530820615中国地质大学（北京）地下水资源评价论文指导老师：崔亚莉日期：2011年6月2日一、用均衡法和达西公式建立河间地块稳定流条件下的地下水平衡方程，并求解1、用单元中心法建立计算网格。15中国地质大学（北京）地下水资源评价论文15中国地质大学（北京）地下水资源评价论文2、建立单元流量平衡方程（1）建立1号计算单元的流量平衡方程（2）建立2号计算单元的流量平衡方程（3）建立3号计算单元的流量平衡方程（4）建立4号计算单元的流量平衡方程（5）建立5号计算单

2、元的流量平衡方程15中国地质大学（北京）地下水资源评价论文3、解单元流量平衡方程组（1）联立方程1号单元：2号单元：3号单元：4号单元：5号单元：（2）确定水文地质参数K＝10m/d；降雨量：445毫米/年；入渗系数：0.35；H0＝53.00m；H6＝52.00m；将水文地质参数带入方程组，得：（3）迭代求解将方程（2）中方程分别改写，其中表示单元i前一次计算水头。15中国地质大学（北京）地下水资源评价论文整理方程(a)~(e)，计算式如下：（4）利用EXCEL绘制迭代曲线15中国地质大学（北京）地下水资源评价论文4、求解地下水补给量和排泄量，并分析误差产生

3、的原因（1）补给量其中，W为降雨入渗量，A为河间地块面积。（2）排泄量其中，为河间地块的宽度。（3）误差分析有限差分法是以每个单元中心点的水头近似代表该单元的水头，然而真实情况中地下水水头是连续变化的，这一近似化与实际情况的区别进而导致了误差的产生。虽然有误差的产生，但只要单元划分的足够多，就可以使所模拟的地下水流场越近似于真实流场，即满足精度要求即可。二、用PMWIN建立地下水稳定流条件下的渗流模型（一）原理1、水量均衡其中，为计算期间内地下水系统各补给量总和；为计算期间内地下水系统各排泄量总和；为计算期间内地下水系统储存量的变化量。2、达西定律其中，Q为通

4、过过水断面的流量；A为过水断面的面积；15中国地质大学（北京）地下水资源评价论文和分别为左、右过水断面处的水头(m);为左、右过水断面之间的水头差(m);L为左、右过水断面之间的距离(m);为水力梯度(无量纲);K为过水断面之间含水介质的渗透系数(m/d)。（二）模型参数1、单元剖分：20列，20行，2层。2、水文地质参数：（1）水平渗透系数Layer1整个层K1=40m/d;Layer2左半部分K2=40m/d，右半部分K3=20m/d(2)垂直渗透系数K1=K2=1m/d3、含水层类型Layer1潜水，代码=1；Layer2承压/无压转化，代码=24、顶板

5、标高Layer1顶板标高58.85m；Layer2顶板标高48.85m5、底板标高Layer1底板标高48.85m；Layer2底板标高41.85m6、降雨补给量计算降雨量：445mm/a；入渗系数：0.35；7、左右侧水头;（三）模拟天然条件下的地下水流场1、天然条件下的地下水流场（如图1）流场分析：流场之所以呈现图1所示形态，原因是第二层左半部分渗透系数为40m/d，而右半部分渗透系数为20m/d。由于左右的渗透系数的差异，导致左半部分等水位线稀疏，而右半部分等水位线密集。2、用WaterBudget分析天然条件下的地下水均衡（见表1）第一层的补给项有两个

6、。一个是降雨入渗补给，补给量大小为1393.2021m^3/d，占总补给量的66%；另一个是下层的越流补给，补给量大小为717.8329m^3/d，占总补给量的34%。第一层的排泄项有两个。一个是侧向流出排泄，排泄量大小为1393.1968m^3/d，占总排泄量的66%；另一个是向下层的越流排泄，排泄量大小为717.8344m^3/d，占总排泄量的34%。补排差为0.0038m^3/d，即补给量近似等于排泄量，所以潜水含水层水量均衡。15中国地质大学（北京）地下水资源评价论文第二层的补给项有一个，是上层的越流补给，补给量大小为717.8344m^3/d。第一层

7、的排泄项有一个，是向上层的越流排泄，排泄量大小为717.8329m^3/d。补排差为0.0015m^3/d，即补给量近似等于排泄量，所以承压含水层水量均衡。图1天然条件下地下水流场等水位线图表1天然条件下地下水均衡表研究范围补给项（m^3/d)排泄项(m^3/d)第一层降水入渗量1393.2021侧向流出量1393.1968侧向流入量0.0000向下层越流排泄量717.8344下层越流补给量717.8329井开采量0.0000小计2111.0350小计2111.0312补排差0.0038第二层上层越流补给量717.8344向上层越流排泄量717.8329小计7

8、17.8344小计717.832915