专利名称：一种计算单井回灌量的有限差分方法
技术领域：
本发明属于工程水文地质、水文水资源、地下水科学领域，具体涉及一种计算单井回灌量的有限差分方法。
背景技术：
近年来，我国的地下水人工补给工作已经有了较大发展，特别是在控制地面沉降、扩大地下水资源开采、利用含水层贮能等方面起了明显作用。所谓地下水人工补给，即通过各种人工措施，把各种地表水资源补充 到含水层内可利用的地下水资源增加。地下水补给主要包括直接和间接两类方法，其中管井注入法是国内外大量使用的直接方法之一。回灌井的工作效率与单井回灌量(回灌井工作时单位时间流入含水层的水量)密切相关，成为工程中重点关注的问题之一。单井回灌量通常可以通过现场试验和理论计算获得。然而回灌井的成井过程中每道工序(冲孔、换浆、下管、填砾、止水、洗井)都与回灌量有直接关系，这造成了现场试验(试验井、试验工况较少情况下)获得的回灌量人为误差较大，不能作为代表值。其次，理论计算由于涉及经验公式、水位抬升量等因素，而大量现场回灌试验表明，常见的回灌量(P)与水位抬升量(S)的关系有直线型、抛物型、幂函数型、对数曲线型等多种形式，因此脱离大量的现场试验也难以直接通过理论计算确定回灌量。在这种情况下，按照工程实际情况建立数值模型，通过数值计算获得单井回灌量已经成为重要的辅助手段。
发明内容
本发明的目的在于提供一种计算单井回灌量的有限差分方法。本发明提出的计算单井回灌量的有限差分方法，包括下述步骤
(1)对井滤管进行简化，使得有限差分模型更为简洁；
(2)对井滤管计算空间进行有限差分划分；
(3)对步骤(2)划分得到的各网格建立节点水均衡有限差分关系；
(4)根据实际工程对简化井所在地层进行水文地质参数赋值，并通过设置已知水头边界或定水头边界、隔水边界来区分滤管、滤管周围填料以及实际地层；然后根据地下水三维非稳定流运动方程建立水均衡线性方程组，并运用Gauss-Seide迭代求解方程；
(5)通过累加流经简化滤管外围相邻网格的水量差来计算单井回灌量。具体步骤如下
(1)、对管井结构进行等效简化
①按照周长相等原则，将井截面由圆形简化为正方形；
②根据滤管开孔率相等原则，对滤管开孔及布置方式进行如下简化将圆形开孔方式简化为方形开孔方式，并使其开孔面积相等，各开孔之间的间距依据开孔率相等原则进行计算确定；
(2)、对管井结构进行网格划分，网格划分的密度不应小于滤管上的方形开孔边长及间距的最小值；
(3)、在步骤(2)得到的网格上进行井滤管骨架管结构的精细模拟井底(沉砂管底端)、滤管壁、滤管顶均定义为死单元，滤管开孔位置设置为定水头边界，滤管周围滤料段及地层根据岩土工程勘察报告提供的水文地质参数进行赋值(渗透系数、储水系数等)；
(4)、在步骤(2)划分得到的各单个网格建立水均衡的有限差分方程及边界条件与初始条件方程，见式I ;
权利要求
1.计算单井回灌量的有限差分方法，其特征在于具体步骤如下 (1)、对管井结构进行等效简化 ①按照周长相等原则，将井截面由圆形简化为正方形； ②根据滤管开孔率相等原则，对滤管开孔及布置方式进行如下简化将圆形开孔方式简化为方形开孔方式，并使其开孔面积相等，各开孔之间的间距依据开孔率相等原则进行计算确定； (2)、对管井结构进行网格划分，网格剖分的密度不应小于滤管上方形开孔边长及间距的最小值； (3)、在步骤(2)得到的网格上进行井滤管骨架管结构的精细模拟井底或沉砂管底端、滤管壁、滤管顶均定义为死单元，滤管开孔位置设置为定水头边界，滤管周围滤料段及地层根据岩土工程勘察报告提供的水文地质参数进行赋值，所述水文地质参数为渗透系数、储水系数； (4)、在步骤(2)划分得到的各单个网格建立水均衡的有限差分方程及边界条件与初始条件方程，见式I ; ■ d {, dh'\ 5 ,, dh) d {, dh' Bh , 、+， _____ I --- . I ] _____ Ir" ______ { -- . Lr· _______ I __— n/ __f y % 5 々、id f I Qx \ ox J φ K ay} az \ m}Bi Ηχ,γ,ζ, ) ^^(χ,γ,ζ, ..................................................GwJeri(I) k( x, yr z, t) |Mij= ^ (xr y, z)..................................................e Ω 式中是Λ ,分别为各向异性主方向渗透系数(m/d)为点(Hz)在f时刻的水头值(m) ψ为源汇项(1/d) 'Ss为储水系数(Ι/m) ; 为时间(d); P为计算域为第一类边界条件为第一类边界条件I]上的水头值(m)； (5)、采用差分近似方法对步骤(4)中式I进行离散化处理，三维流某节点与周围节点水均衡差分关系见式(2)所示 Bi K ,t+ K χ—, [(Mr4办(Ak)W) k^*l\ -~ MX*IK I (Aj)3(Δζ J f y*4；i φ-m timn _ Li^n d^n — 一 +(1^ r t-K ,...........................-Γ , ^^ (2)I 4# 御 1 如 _-) ;冲(_ ·Φ (Aj)* _ &(* )g**n _ I v'Jwi {hzf(ΔζΥ J Q- 画知一^ 式中4 表示点(UA)上々时刻的水头；OWSl , 5常取三种情况当5=0时称有限差分显不差分格式，当5 = 0.5时称有限差分中心差分格式，当Θ= I时称有限差分隐式差分格式；在离散过程中采取中心差分格式，解是无条件稳定的，对离散的线性方程进行Gauss-Seide迭代加快收敛速度； (6)、 对滤管上的所有定水头边界外围网格中流出的水量进行累加，累计结果即为单井单井回灌量。
全文摘要
本发明涉及一种计算单井回灌量的有限差分方法，具体为对井滤管进行等效简化处理，对井滤管计算空间进行有限差分划分；划分得到的各网格建立节点水均衡有限差分关系；根据实际工程对简化井所在地层进行水文地质参数赋值，并通过设置已知水头边界或定水头边界、隔水边界来区分简化滤管、滤管周围填料以及实际地层；然后根据地下水三维非稳定流运动方程建立水均衡线性方程组，运用Gauss-Seide迭代求解方程；通过累加流经简化滤管外围相邻网格的水量差来计算回灌量。本发明针对以往基坑工程中单井回灌量主要依靠工程类比和成井后回灌试验确定，事前计算预测困难的问题，依据水文地质勘察成果和回灌要求，采用有限差分方法预测不同滤管长度和不同滤管位置回灌井的单井回灌量，为整体回灌方案的设计及优化提供低成本、快速和定量化的依据。
文档编号G06F19/00GK102855390SQ20121028094
公开日2013年1月2日 申请日期2012年8月9日 优先权日2012年8月9日
发明者王建秀, 胡蒙蒙, 吴远斌, 隋东昌, 田普卓, 汪沛, 冯波 申请人:同济大学