[发明专利]一种计算单井回灌量的有限差分方法

权利要求书

1.计算单井回灌量的有限差分方法，其特征在于具体步骤如下：

(1)、对管井结构进行等效简化

①按照周长相等原则，将井截面由圆形简化为正方形；

②根据滤管开孔率相等原则，对滤管开孔及布置方式进行如下简化：将圆形开孔方式简化为方形开孔方式，并使其开孔面积相等，各开孔之间的间距依据开孔率相等原则进行计算确定；

(2)、对管井结构进行网格划分，网格剖分的密度不应小于滤管上方形开孔边长及间距的最小值；

(3)、在步骤(2)得到的网格上进行井滤管骨架管结构的精细模拟：井底或沉砂管底端、滤管壁、滤管顶均定义为死单元，滤管开孔位置设置为定水头边界，滤管周围滤料段及地层根据岩土工程勘察报告提供的水文地质参数进行赋值，所述水文地质参数为渗透系数、储水系数；

(4)、在步骤(2)划分得到的各单个网格建立水均衡的有限差分方程及边界条件与初始条件方程，见式1；

        (1)

式中：，  ，  分别为各向异性主方向渗透系数（m/d）；为点在时刻的水头值 (m)； 为源汇项 (1/d)； 为储水系数 (1/m)；为时间 (d); 为计算域；为第一类边界条件；为第一类边界条件  上的水头值 (m)；

(5)、采用差分近似方法对步骤(4)中式1进行离散化处理，三维流某节点与周围节点水均衡差分关系见式（2）所示：

  （2）

式中：表示点上时刻的水头；，常取三种情况：当时称有限差分显示差分格式，当时称有限差分中心差分格式，当时称有限差分隐式差分格式；在离散过程中采取中心差分格式，解是无条件稳定的，对离散的线性方程进行Gauss-Seide迭代加快收敛速度；

(6)、对滤管上的所有定水头边界外围网格中流出的水量进行累加，累计结果即为单井单井回灌量。