如何计算地下水流模型中的流量边界

在进行流量边界的设置时，通常使用达西定律来计算流量。

基于使用Visual MODFLOW设置好的地层高程的基础上，进行接下来的计算。

达西定律中，

K表示渗透系数，m/d；A表示过水断面的面积，；I表示水力梯度，无量纲。

K通常使用岩性经验值，A面积可由网格的长乘高得出，I水力梯度可以参考上一篇博客计算。最后根据流量的流入（+）或者流出（-）来确定值得正负。

1、在ArcGIS中，根据边界位置，创建点文件，点出抽水井的位置（Visual MODFLOW中用虚拟井来表示流量边界）；

2、将Visual MODFLOW中对应流量边界所在层位的高程导出为shp面文件，在ArcGIS-分析工具-叠加分析-空间连接中，连接上一步的抽水井点文件和地层面文件，这一步获取网格的长和宽，即可计算出A面积；

3、参考上一篇博客计算水力梯度得到I；

4、利用ArcGIS的文件属性表的字段计算器就可以完成的计算。