CN102866983A - 一种精细模拟管井结构的有限差分方法 - Google Patents

一种精细模拟管井结构的有限差分方法 Download PDF

Info

Publication number
CN102866983A
CN102866983A CN2012102811308A CN201210281130A CN102866983A CN 102866983 A CN102866983 A CN 102866983A CN 2012102811308 A CN2012102811308 A CN 2012102811308A CN 201210281130 A CN201210281130 A CN 201210281130A CN 102866983 A CN102866983 A CN 102866983A
Authority
CN
China
Prior art keywords
finite difference
well
chimney filter
difference
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012102811308A
Other languages
English (en)
Other versions
CN102866983B (zh
Inventor
王建秀
胡蒙蒙
吴远斌
汪沛
田普卓
隋东昌
冯波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN201210281130.8A priority Critical patent/CN102866983B/zh
Publication of CN102866983A publication Critical patent/CN102866983A/zh
Application granted granted Critical
Publication of CN102866983B publication Critical patent/CN102866983B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种精细模拟管井结构的有限差分方法，它首先根据研究需要对井计算空间进行有限差分网格划分，在网格节点上建立水均衡的有限差分关系；接着根据实际工程对井所在地层进行水文地质参数赋值，在井滤管处通过设置已知水头边界或定水头边界、隔水边界来区分井壁、滤管开孔、滤管周围填料以及实际地层；然后根据地下水三维非稳定流运动方程建立水均衡线性方程组，最后通过Gauss-Seide迭代求解方程，得到节点水头值。本发明简单有效可适用于任何井精细结构的有限差分数值模拟设计。

Description

一种精细模拟管井结构的有限差分方法
技术领域
本发明属于工程水文地质、水文水资源、环境地质工程领域，具体涉及一种精细模拟井结构（主要指井滤管的开孔方式、长度、填砾和开孔直径）的有限差分方法。
背景技术
    深基坑工程对承压含水层进行减压降水时，周围地层可能产生不均匀沉降，管井回灌是重要的主动控沉措施，回灌井由井口、井管、滤管和沉沙管组成，要正确评价管井的工作效率，必须设计回灌井的滤管。以往的工程经验认为，回灌井回灌效果除地层水文地质条件影响外主要与滤管的直径、长度及开孔率有关，特别是在滤管直径和长度都确定条件下，回灌效率仅与滤管开孔率相关，而忽略了开孔方式的影响。为了深入探究管井结构对工作效率的影响，若采用现场对比试验，由于成井工艺差异，少量的试验井难以达到试验效果，而精细模拟井结构就成为首选的经济有效的研究手段。
    就回灌井而言，在回灌过程中，由于回灌井中的回灌水位与地下水的静水位(承压的和非承压的)之间形成一个水头差，注入回灌井里的水才有可能向含水层渗流。当渗流稳定时，则回灌水位就不再继续上升而稳定下来，此时在回灌井周围形成一个水位降深反漏斗（即降深为负值）。反漏斗的形状是回灌井工作效率的直观体现，相对而言，平缓型回灌水位降深反漏斗对于工程更有意义。
发明内容
当水文地质条件、回灌压力相同情况下，直径相同的回灌井的工作效率不仅仅与滤管长度及其开孔率有关，还与滤管的开孔方式（如圆孔、条孔等）有关，为了利用数值模拟来研究井结构对井效率的影响，本发明的目的在于提供一种精细模拟管井结构的有限差分方法。
本发明提出的精细模拟管井结构的有限差分方法，具体步骤如下：
(1)、对管井结构进行等效简化
①按照周长相等原则，将井截面由圆形简化为正方形；
②根据滤管开孔率相等原则，对滤管开孔及布置方式进行如下简化：将圆形开孔方式简化为方形开孔方式，并使其开孔面积相等，各开孔之间的间距依据开孔率相等原则进行计算确定；
(2)、对管井结构进行网格划分，网格剖分的密度不应小于滤管上方形开孔边长及间距的最小值；
(3)、在步骤(2)得到的网格上进行井过滤器骨架管结构的精细模拟：井底（沉砂管底端）、滤管壁、滤管顶均定义为死单元，滤管开孔位置设置为定水头边界，滤管周围滤料段及地层根据岩土工程勘察报告提供的水文地质参数进行赋值（渗透系数、储水系数等）；
(4)、在步骤(3)划分得到的各单个网格建立水均衡的有限差分方程及边界条件与初始条件方程，见式(1)；
Figure 2012102811308100002DEST_PATH_IMAGE002
        (1)
式中：
Figure DEST_PATH_IMAGE004
， 
Figure DEST_PATH_IMAGE006
 ， 
Figure DEST_PATH_IMAGE008
 分别为各向异性主方向渗透系数（m/d）；
Figure DEST_PATH_IMAGE010
为点
Figure DEST_PATH_IMAGE012
在
Figure DEST_PATH_IMAGE014
时刻的水头值 (m)；
Figure DEST_PATH_IMAGE016
 为源汇项 (1/d)； 为储水系数 (1/m)；
Figure 253628DEST_PATH_IMAGE014
为时间 (d); 
Figure DEST_PATH_IMAGE020
为计算域；
Figure DEST_PATH_IMAGE022
为第一类边界条件；为第一类边界条件 
Figure DEST_PATH_IMAGE022A
 上的水头值 (m)。
(5)、采用差分近似方法对步骤(4)中所述式(1)进行离散化处理，三维流某节点与周围节点水均衡差分关系见式（2）所示：
Figure DEST_PATH_IMAGE027
  （2）
式中：
Figure DEST_PATH_IMAGE029
表示点
Figure DEST_PATH_IMAGE031
上
Figure DEST_PATH_IMAGE033
时刻的水头；
Figure DEST_PATH_IMAGE035
，常取三种情况：当
Figure DEST_PATH_IMAGE039
时称有限差分显示差分格式，当
Figure DEST_PATH_IMAGE041
时称有限差分中心差分格式，当时称有限差分隐式差分格式。在离散过程中采取中心差分格式，解是无条件稳定的，对离散的线性方程进行Gauss-Seide迭代加快收敛速度。
与现有发明相比，本发明的有益效果是：针对复杂地质条件下井结构选型及设计主要依靠经验和理想试验，对于现场变异地质条件适应性差、成本高、周期长的问题，发明了利用有限差分方法精细模拟管井工作效率的方法，可以针对特定地质条件模拟井过滤器的开孔方式、空隙率、滤料组合条件下的工作效率问题，为管井精细设计及优化提供了高效、低成本、快速实用的计算方法。
附图说明
图1是滤管开孔简化示意图。其中：(a)为井实际滤管开孔（圆形开孔）布置示意图，(b)为简化后的滤管开孔布置示意图。S 1、 S 2 是实际滤管上开孔横向间隔和竖向间隔，d是圆形开孔直径，S 3 是简化后滤管开孔横向间隔，而竖向间隔简化为0，b是方形开孔边长。
图2是滤管开孔简化示意图。其中：（a）为实际滤管开孔及布置方式，（b）为简化后滤管开孔及布置方式。
图3是简化滤管有限差分网格划分示意图。其中：（a）为横向剖面，（b）为竖向剖面。S 1是简化滤管以内区域， S 2是简化滤管以外滤料填充区域，S 3是在滤管顶及滤管底设置的隔水边界，p是简化滤管所在区域，滤管壁通过p区域各小均衡域的中心。
图4是滤管（相邻两层滤管孔）有限差分网格划分横向剖面详图。S 1 是井滤管以内区域，等压回灌过程中该区域设置为等压边界，q 1 为滤管开孔区域，q 2 为滤管封闭区域（管材部分）。
图5是滤管有限差分网格划分竖向剖面详图。S 1 是井滤管以内区域，等压回灌过程中该区域设置为等压边界，q 1 为滤管开孔区域，q 2 为滤管封闭区域（管材部分），S 4 是在滤管顶及滤管底设置的隔水边界。
具体实施方式
下面通过实施例结合附图进一步说明本发明。
实施例1：以下以滤管位于承压含水层完整井为例说明本发明方法，其中滤管开孔方式为圆形品字形布置、开孔率假设为25%，为了简易起见，假设地表标高为0m，地层一共设为四层，从上到下分别为：
    
Figure DEST_PATH_IMAGE045
层：潜水含水层，厚5m，水头为-2m；
    
Figure DEST_PATH_IMAGE047
层：隔水层，厚3m，水头为-2m；
    
Figure DEST_PATH_IMAGE049
层：承压含水层，厚10m，水头为-2m，滤管所在层；
    层；隔水层，厚2m，水头为-2m。
1、图1为滤管等效图，a图为等效前滤管开孔布置图，圆孔以品字形分布，相邻两横排开孔间存在狭小间隙S 2 ，为了便于有限差分网格划分，在开孔率相同情况下将圆孔等效为方孔，并对其布置稍作调整，即把间隙S 2 消除使得上下两层方孔紧密排列，如b图所示。
                              （1）
式中：
Figure DEST_PATH_IMAGE055
为圆孔个数，
Figure DEST_PATH_IMAGE057
为圆孔直径，
Figure DEST_PATH_IMAGE059
为方孔个数，
Figure DEST_PATH_IMAGE061
为方孔边长。
2、根据研究需要把计算空间划分为若干个单元（即有限差分网格划分），每个单元上有，一个计算节点，在节点上进行水均衡计算，在平面直角坐标系中，首先分别在x轴、y轴方向上进行均匀划分，然后在拟建井附近进行网格加密；在竖向上沿z轴均匀划分，划分层数与等效后滤管方孔层数相同，即以方孔边长b为间隔均匀划分，如图5所示。
3、若地下水运动为三维非稳定流且符合达西定律，承压含水层为非均质、各向异性，含水层及其中地下水视为弹性体，水头改变时弹性储量释放（或储存）是瞬时完成的，划分网格得到的各个小均衡域的水均衡符合下式：
Figure DEST_PATH_IMAGE063
      （2）
式中：
Figure DEST_PATH_IMAGE004A
，
Figure DEST_PATH_IMAGE006A
，
Figure DEST_PATH_IMAGE008A
分别为各向异性主方向渗透系数（m/d）；
Figure DEST_PATH_IMAGE010A
为点
Figure 663444DEST_PATH_IMAGE012
在
Figure 131597DEST_PATH_IMAGE014
时刻的水头值 (m)；
Figure DEST_PATH_IMAGE016A
为源汇项(1/d)；
Figure DEST_PATH_IMAGE018A
为储水系数 (1/m)；
Figure 212292DEST_PATH_IMAGE014
为时间(d)；为计算域。以上方程加上定解条件就构成非均质各向异性孔隙介质地下水数学模型，该数学模型是抛物型的，具有耗散效应，流场中的任何扰动都会影响全流场。采用差分近似方法对偏微分方程进行离散化处理，三维流某节点与周围节点水均衡差分关系如公式（3）所示：
Figure DEST_PATH_IMAGE065
                                                                                    （3）
式中：
Figure 964347DEST_PATH_IMAGE029
表示点（i，j，k）上t n 时刻的水头；，θ常取三种情况：当θ=0时称有限差分显示差分格式，当θ=0.5时称有限差分中心差分格式，当θ=1时称有限差分隐式差分格式。在离散过程中采取中心差分格式，解是无条件稳定的，对离散的线性方程进行Gauss-Seide迭代加快收敛速度。
4、模拟等压回灌，在式（1）中需要两类边界条件，一种为定水头边界，通过节点赋值来实现；另一种为隔水边界，通过设置死单元来实现。
5、对已划分好的网格进行定义水文地质参数。图2为井横向剖面示意图，图3为竖向剖面示意图，井、滤料填充区、地层区分别定义其水文地质参数，S 1 （井滤管以内区域）中水头赋值为回灌定水压，滤管开孔位置（图4、图5中的q 1区域）与S 1 设置相同的定水压，代表回灌水以等压方式向地层中渗流，滤管处的井壁部分（图4、图5中的q 2 区域）设置死单元，其中滤管部分详图见图4、图5。由于在回灌过程中井滤管顶部和底部不参与回灌，故在图3中设置死单元区S 4 ，结合图1，由于滤管开孔是交错分布，故在横剖面上表现为两种形式（图4），在纵剖面上表现为相邻交错分布（图5）。
6、完成了井结构（主要是滤管开孔方式）有限差分精细模拟，通过仿真计算反演求参，设计的结果满足设计条件并可用于指导实践。

Claims (1)

1.一种精细模拟管井结构的有限差分方法，其特征在于具体步骤如下：
(1)、对管井结构进行等效简化
①按照周长相等原则，将井截面由圆形简化为正方形；
②根据滤管开孔率相等原则，对滤管开孔及布置方式进行如下简化：将圆形开孔方式简化为方形开孔方式，并使其开孔面积相等，各开孔之间的间距依据开孔率相等原则进行计算确定；
(2)、对管井结构进行网格划分，网格剖分的密度不应小于滤管上方形开孔边长及间距的最小值；
(3)、在步骤(2)得到的网格上进行井过滤器骨架管结构的精细模拟：井底或沉砂管底端、滤管壁、滤管顶均定义为死单元，滤管开孔位置设置为定水头边界，滤管周围滤料段及地层根据岩土工程勘察报告提供的水文地质参数进行赋值，所述水文地质参数为渗透系数或储水系数；
(4)、在步骤(3)划分得到的各单个网格建立水均衡的有限差分方程及边界条件与初始条件方程，见式(1)；
Figure 2012102811308100001DEST_PATH_IMAGE001
        (1)
式中：
Figure DEST_PATH_IMAGE003
， 
Figure DEST_PATH_IMAGE005
 ，  分别为各向异性主方向渗透系数（m/d）；为点在
Figure DEST_PATH_IMAGE011
时刻的水头值 (m)； 为源汇项 (1/d)；
Figure DEST_PATH_IMAGE015
 为储水系数 (1/m)；为时间 (d); 
Figure DEST_PATH_IMAGE017
为计算域；为第一类边界条件；
Figure DEST_PATH_IMAGE021
为第一类边界条件 
Figure 657463DEST_PATH_IMAGE022
 上的水头值 (m)；
(5)、采用差分近似方法对步骤(4)中所述式(1)进行离散化处理，三维流某节点与周围节点水均衡差分关系见式（2）所示：
Figure DEST_PATH_IMAGE023
  （2）
式中：
Figure 37891DEST_PATH_IMAGE024
表示点上
Figure DEST_PATH_IMAGE026
时刻的水头；
Figure 522706DEST_PATH_IMAGE027
，常取三种情况：当时称有限差分显示差分格式，当
Figure DEST_PATH_IMAGE030
时称有限差分中心差分格式，当
Figure 226668DEST_PATH_IMAGE031
时称有限差分隐式差分格式；在离散过程中采取中心差分格式，解是无条件稳定的，对离散的线性方程进行Gauss-Seide迭代加快收敛速度。
CN201210281130.8A 2012-08-09 2012-08-09 一种精细模拟管井结构的有限差分方法 Expired - Fee Related CN102866983B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210281130.8A CN102866983B (zh) 2012-08-09 2012-08-09 一种精细模拟管井结构的有限差分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210281130.8A CN102866983B (zh) 2012-08-09 2012-08-09 一种精细模拟管井结构的有限差分方法

Publications (2)

Publication Number Publication Date
CN102866983A true CN102866983A (zh) 2013-01-09
CN102866983B CN102866983B (zh) 2016-06-29

Family

ID=47445861

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210281130.8A Expired - Fee Related CN102866983B (zh) 2012-08-09 2012-08-09 一种精细模拟管井结构的有限差分方法

Country Status (1)

Country Link
CN (1) CN102866983B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105386474A (zh) * 2015-10-22 2016-03-09 上海交通大学 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法
CN113050190A (zh) * 2021-03-03 2021-06-29 河北益坤岩土工程新技术有限公司 直线边界非稳定流抽水试验水文地质参数智能计算方法
CN117848423A (zh) * 2024-03-07 2024-04-09 南京中鑫智电科技有限公司 一种换流变阀侧套管壳体完整性的在线监测方法、***、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1470738A (zh) * 2002-07-26 2004-01-28 中国石油天然气股份有限公司 光油管井氮气隔热井筒传热模拟技术
US20090294122A1 (en) * 2006-05-24 2009-12-03 Jens Henrik Hansen Flow simulation in a well or pipe
CN102063577A (zh) * 2011-01-13 2011-05-18 黄河勘测规划设计有限公司 一种基于等势面的地下水渗流量计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1470738A (zh) * 2002-07-26 2004-01-28 中国石油天然气股份有限公司 光油管井氮气隔热井筒传热模拟技术
US20090294122A1 (en) * 2006-05-24 2009-12-03 Jens Henrik Hansen Flow simulation in a well or pipe
CN102063577A (zh) * 2011-01-13 2011-05-18 黄河勘测规划设计有限公司 一种基于等势面的地下水渗流量计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘持峰等: "渗渠产水量与滤管开孔率的关系研究", 《给水排水》, vol. 34, no. 12, 31 December 2008 (2008-12-31), pages 48 - 51 *
宋若龙等: "非轴对称套管井中声场的并行有限差分模拟", 《地球物理学报》, vol. 53, no. 11, 30 November 2010 (2010-11-30), pages 2767 - 2775 *
袁建伟等: "井结构参数对同井回灌地下水源热泵的影响", 《制冷与空调》, vol. 9, no. 3, 30 June 2009 (2009-06-30), pages 82 - 86 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105386474A (zh) * 2015-10-22 2016-03-09 上海交通大学 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法
CN105386474B (zh) * 2015-10-22 2017-05-10 上海交通大学 确定基坑开挖面上方止水帷幕渗漏对周边环境影响的方法
CN113050190A (zh) * 2021-03-03 2021-06-29 河北益坤岩土工程新技术有限公司 直线边界非稳定流抽水试验水文地质参数智能计算方法
CN117848423A (zh) * 2024-03-07 2024-04-09 南京中鑫智电科技有限公司 一种换流变阀侧套管壳体完整性的在线监测方法、***、设备及介质
CN117848423B (zh) * 2024-03-07 2024-05-17 南京中鑫智电科技有限公司 一种换流变阀侧套管壳体完整性的在线监测方法、***、设备及介质

Also Published As

Publication number Publication date
CN102866983B (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN103745022B (zh) 聚合物驱后变流线调整方法
CN105089582B (zh) 基于井下流量控制设备的油藏数值模拟方法及装置
CN102041995B (zh) 复杂油藏水淹状况监测***
CN104809277A (zh) 一种超低渗致密储层水平井地质建模方法
CN107664029B (zh) 砂岩热储地热资源可循环利用最优井网布局方法
US20160376885A1 (en) Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs
CN110838175B (zh) 注气开发油藏的地质模型建立方法
CN104632200A (zh) 河流相多层及厚层油藏精细建模方法
CN102930345A (zh) 一种基于梯度算法的自适应井网优化方法
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
CN112394404B (zh) 一种渐进式储层精细表征方法
Dong et al. Optimization of Mine Drainage Capacity Using FEFLOW for the No. 14 Coal Seam of China's Linnancang Coal Mine
CN114239904A (zh) 一种地下水管理方法及装置
CN103455667A (zh) 充气法治理承压含水层海水入侵的数值模拟方法
CN104533519B (zh) 立井井筒通过强含水厚岩层时涌水水害的治理方法
CN101915088A (zh) 一种油气运移路径生成方法与装置
CN102866983A (zh) 一种精细模拟管井结构的有限差分方法
CN116432546A (zh) 一种水平井网地浸采铀井储耦合模拟方法及***
CN102855390B (zh) 一种计算单井回灌量的有限差分方法
CN110687603B (zh) 一种海上油田储层内部渗流屏障的地质建模方法
Jiang et al. A new variational inequality formulation for unconfined seepage flow through fracture networks
CN111381277B (zh) 碳酸盐岩等时地层格架建立方法及装置
CN105275457B (zh) 整装老油田河流相储层水平井部署最佳模式研究方法
CN110232245B (zh) 悬挂式帷幕承压水基坑定流量抽水三维流场水位分布确定方法
Luo et al. 3-D variable parameter numerical model for evaluation of the planned exploitable groundwater resource in regional unconsolidated sediments

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160629

Termination date: 20180809

CF01 Termination of patent right due to non-payment of annual fee