Modflow抽水井案例模拟1

Modflow抽水井案例模拟1

案例描述

题干：
一侧向无限分布的承压含水层，其导水系数T=200m²，贮水系数S=2*10^-3.有一完整生产井以抽水量Q=3140m³/d进行开采。试计算距井300m处，开采10、20、30时的降深。

Modflow抽水井案例模拟1M=10m(含水层厚度) 初始水位H0=20m(承压含水层)
渗透系数:K=20m/d

Modflow抽水井案例模拟1
s1=4.048m;s2=4.927m;s3=5.401m
h1=15.952m;h2=15.073;h3=14.599

软件版本：Visaul Modflow 2011
界面介绍：
File :文件
Input:输入信息
Run:运行
Output:输出
PS：路径中最好不要有中文出现

点击Next

由于仅需要10d、20d、30d，所以将渗流模拟时间设置为50d备用。

网格划分

注意：Columns和Rows设置为101是因为网格为奇数个时，可以找到一个中心点网格，左右两边网格数相同，毕竟对称比较舒服。

点击Finish

本模型较为简单，网格部分就此设置完成，其余教程再细说网格划分。
点击Wells设置井，此处我们需要两种井，一种是抽水井，一种是观测井，抽水井参与实际模型计算，观测井顾名思义只是用来观察的，可以认为观测井是用来观察水位并和计算结果进行对比的，后续计算完成后会再提到。

首先设置降水井
点击界面左侧Add Well, 然后在右侧网格中随便点击一个点，再输入横纵坐标确定井的准确位置，我这里输入0，0


设置完后弹出井的信息界面

点击OK后出现此图

再设置观测井 点击Add Obs.

其他操作类似


接下来分别设置渗透性、贮水性和水头参数

渗透性—database

贮水性—database

水头参数—database

井的信息设置完成，Input界面就完成了
接下来点击右下角Main Menu回到主菜单，再点击Run

首先进行时间设置(默认即可)

层设置
Layer：将层设置为承压水层

点击运行Run即可

运行结果显示:回到Main Menu—点击Output出现以下界面

点击左上角View Row并将剖面线放在观测井与抽水井连线上，如图：

将数轴扩大50倍方便观看

再左边Options选项中点击 Use Color Shading显示彩色云图


接下来绘制Time—head 图

再W1.1-1前打勾并点击Apply

从图中可见，观测井有三个点(较大方块)，曲线为数值结果，结果表明，数值模拟结果和观测井(解析解)较为接近，十分可信！！