Experimental Study on Plugging Mechanisms and Influencing Factors of a New Foam Metal Screen
-
摘要: 为了解新型泡沫金属筛管的堵塞机理和优化筛管结构,利用全尺寸防砂试验装置模拟泡沫金属筛管的堵塞过程,采用比采油指数评价筛管堵塞程度,分析不同泥质含量、筛管精度和泡沫金属层厚度下的堵塞程度及堵塞机理。结果表明,表层颗粒架桥是泡沫金属筛管堵塞的主要原因,筛管的比采油指数随泥质含量增大呈指数函数衰减,最终约稳定在0.075 m3/(d·MPa·m),与不含泥质时相比,下降了约78%;泡沫金属层孔径增大,比采油指数随之增大;泡沫金属层厚度增大,比采油指数的变化很小。新型泡沫金属筛管堵塞机理及影响因素的研究结果,为进一步优化筛管结构提供了理论依据。Abstract: In order to study the plugging mechanisms and optimize the structure of foam metal screen,the plugging process for foam metal screens was simulated by a large-scale sand control simulation device. With the screen plugging degree evaluating the specific productivity index, the plugging degree and mechanism under the different shale content, the screen precision and the thickness of the foam metal were analyzed. The experimental results showed that surface particle bridging was the main cause of foam metal screen plugging. The specific productivity index exponentially decreases with the increase of shale content where it finally stabilized at 0.075 m3/(d·MPa·m), which decreased by 78% compared with the index where there was no shale content. When the pore size of foam metal is increased, the specific productivity index increases. The increase of foam metal thickness had slight impact on the specific productivity index. The experimental study of plugging mechanism and influencing factors of foam metal screen provided theoretical basis for the optimization of screen structure.
-
Keywords:
- metal screen with foam /
- plugging mechanism /
- surface particle /
- bridging /
- shale content /
- sand control
-
深层干热岩储层具有温度高、硬度大和基质渗透率低等特点，需采用水力压裂等方法改造致密干热岩储层，构建复杂的裂缝网络体系，提高循环工质换热体积，实现高效循环取热[1–6]。干热岩水力压裂是一个多物理场耦合问题，需要考虑位移、流体压力、温度等多个物理量之间的相互影响[7–17]。H. Hofmann等人[18]通过研究发现，复杂裂缝网络有助于提高干热岩储层开采时的经济效益。Zhang Yanjun等人[19]提出，水平井分段压裂技术可提高循环工质的换热效率，从而提高裂缝型干热岩储层的开采效益。谢紫霄等人[20]通过干热岩压裂物理模拟试验，发现水力裂缝沟通天然裂缝后能显著提高裂缝渗流系数，实现地热能的高效开发。干热岩储层裂缝扩展过程中受多因素综合影响，且温度影响不容忽略。
相场法是一种基于系统能量最小化原理判断裂缝扩展的有限元方法，在处理多物理场耦合问题时具有网格处理简单、结果准确性高和计算量小等优点，近年来被广泛应用于水力压裂数值模拟领域[21–26]。侯冰等人[23]基于相场法建立了流−固−化耦合酸压裂缝扩展模型，分析了孔洞结构对酸压裂缝的影响。N. Noii等人[25]提出适用于正交各向异性的相场裂缝描述方法，对相场裂缝建模进行了初步分析。
相场法在处理复杂地层多物理场耦合裂缝扩展模拟方面具有显著优势。基于相场法建立了流−固−热耦合裂缝扩展数值模型，采用自适应网格方法平衡计算成本与计算精度，通过解析解验证模拟结果的准确性；在此基础上，研究了干热岩储层水力裂缝的起裂与扩展规律，分析了地质与工程因素对干热岩储层水力裂缝扩展的影响，以期为深层干热岩储层水力压裂提供参考。
1. 深层干热岩储层水力压裂数值模型
基于相场法构建流−固−热耦合裂缝扩展数值模型时，其核心思想是引入无量纲相场参数ϕ描述空间点的物理状态。ϕ的取值在 [0,1]区间，ϕ= 0 表示材料断裂，ϕ= 1 表示材料未断裂， 0 < ϕ < 1 表示材料处于裂缝过渡区域。该参数可以描述裂缝形态，同时可作为指示参数耦合多物理场方程。
1.1 控制方程
1.1.1 相场−位移控制方程
在热弹性多孔介质的相场裂缝模型能量泛函基础上[24–25]，将裂缝扩展问题转化为系统能量最小化问题并求解，可得相场−位移控制方程强形式：
−[(1−k)ϕε(u)−Gcε∇ϕ−Gcε1−ϕ−2(αB−1)ϕ(p−p0)∇⋅u+2ϕ∇(p−p0)u−2(αTKd+C)ϕ(θ−θ0)∇⋅u+2Cϕ∇(θ−θ0)u]∂tϕ=0 (1) −∇⋅[((1−k)ϕ2+k)σ+(u)+σ−(u)]+(α−1)∇(ϕ2p)+ϕ2∇p+(αTKd+CT)∇(ϕ2θ)+ϕ2∇θ=0 (2) 式中：u为位移，m；k为正则化参数；ε(u)为应变；ε为正则化参数；αB为Biot系数；αT为热膨胀系数，℃－1 ；p为流体压力，Pa；p0为初始压力，Pa；
θ 为温度，℃；θ0 为初始温度，℃；C为界面热应力常数，Pa/℃；σ+(u)为裂缝驱动力，Pa；σ－(u)为裂缝非驱动力，Pa；在二维问题中，Kd=μ+λ，μ，λ为Lame系数；Gc为临界能量释放率，Pa·m。1.1.2 温度场控制方程
基于相场参数的定义，可构建指示函数：
Xf=1−ϕ，Xr=ϕ (3) 基于局部热平衡假设，同一空间点处温度值具有唯一性。引入指示函数后，得到温度场控制方程的有限元强形式：
(XfCf+XrCr)∂θ∂t=(XfKf+XrKr)[∂∂x(∂θ∂x)+∂∂y(∂θ∂y)]+v∇θ+Qi (4) 式中：Cf为压裂液热容，kJ·m3/℃ ；Cr为储层热容，kJ·m3/℃ ；Kf为压裂液热传导系数，W/(m·℃)；Kr为地层岩石热传导系数，W/(m·℃)； v为流体流速，m/s；Qi为注入热源，J。
1.1.3 流体压力控制方程
对于流体场，其物质守恒方程为：
∂t(ρfφ∗)+∇⋅(ρfv)=q (5) 式中：ρf为液体密度，g/cm3；φ*为孔隙度；q为单位时间内注入的压裂液质量，g。
计算流体流速时忽略重力影响，表达式为：
v=−KηΔp (6) 式中：K为渗透率，D；η为液体黏度，mPa·s。
φ∗=φ∗0+αB∇⋅u+1M(p−p0)+KT(θ−θ0) (7) 式中：
φ∗0 为初始时刻孔隙度；M为Biot系数；KT为孔隙热膨胀系数，℃−1。流体压力控制方程的强形式为：
ρf∂t(αB∇⋅u+1Mp+KTθ)−ρfKη∇⋅Δp=q (8) 1.2 自适应网格
模型中以相场参数为指示参数，引入自适应网格技术，在保证模拟精度的前提下，减少模型的计算量。自适应网格算法流程如图1所示。
自适应网格可对同一时间步内的裂缝扩展状态进行多次求解，但与全局细分网格相比节点数量少、计算简单，计算精度仍保持良好[26]。
1.3 多物理场耦合求解
模型共包括固体场模块、温度场模块和流体场模块等3个模块，模型中包括4个解变量，分别为相场参数（ϕ）、位移（u）、温度（T）和流体压力（p）。在固体场模块中，采用全耦合方法求解相场参数与位移；利用温度场模块、流体场模块分别求解温度和压力值，并采用固定应力迭代的方式与固体场耦合。具体流程如下：
1）在任意时间步tn+1，给定上一个时间步的解θn，un，ϕn和pn。
2）基于固定应力法，构建迭代序列l=0，1，2…，定义初始解变量ul=0=un，ϕl=0=ϕn，pl=0=pn，θl=0=θn。
3）给定θn，un，ϕn，pn，ul，ϕl和pl，求解式（4），得到θl+1。
4）给定θn，un，ϕn，pn，ul，ϕl和θl，求解式（8）获得pl+1。
5）给定θn，un，ϕn，pn，pl和θl，求解式（1）和式（2），可得{ul+1，ϕ l+1}。
6）若max{||pl+1−pl||，||θl+1−|θl||，||ul+1−ul||，||ϕl+1−ϕl||}≤TOLfs，则迭代收敛（TOLfs为固定应力法容差，取10 －8），令pn+1=pl+1，θn+1=θl+1，un+1=ul+1，ϕn+1=ϕl+1；若迭代不收敛，则令ul+1=ul，ϕl+1=ϕl，pl+1=pl，θl+1=θl，重复步骤3）—5），直至满足收敛条件。
7）检查网格系统是否满足计算精度要求，若不满足，则更新自适应网格系统。
8）输出θn+1，un+1， ϕn+1和pn+1作为时间步tn+1的解变量，并进行下一个时间步的计算。
1.4 初始条件与边界条件
模拟二维条件下的裂缝扩展过程，忽略重力影响与垂向应变。模型初始状态及边界条件如图2所示。
模型长30 cm，宽30 cm。设置天然裂缝等结构时，将对应位置相场参数固定为0，其余岩石基质区域处相场参数设置为1，进行模拟分析，根据相场参数判断裂缝扩展情况。设定压裂液温度为20 ℃，岩石弹性模量为40 GPa，泊松比为0.20，岩石密度为2.50 g/cm3，Biot系数为0.80，渗透率为1 mD，地层岩石热容为2.0 MJ·m3/℃ ，压裂液热容4.2 MJ·m3/℃ ，地层岩石热传导系数为0.63 W/(m·℃)。
1.5 解析解验证
相场裂缝模型中，最小网格尺寸受网格细分次数控制。采用自适应网格方法，全局细分5次、局部细分3次时，最小网格尺寸为1.17 mm；采用全局网格细分方法，全局细分8次时，最小网格尺寸同样为1.17 mm。固定边界位移，分别采用全局细分网格与自适应网格对裂缝扩展进行模拟，模拟结果如图3所示。
![]()
图 3 固定边界位移条件下干热岩水力压裂数值模拟Figure 3. Numerical simulation of hydraulic fracturing in HDR reservoir under fixed boundary displacement无限大线弹性物体内有一条缝内压力均匀分布的裂缝，其裂缝宽度的计算公式为[27]：
w=4(1−ν2)pfE√c2−y2 (9) 式中：w为裂缝宽度，mm；pf 为裂缝内流体压力，Pa；c为裂缝半缝长，mm；y为裂缝内某点到裂缝中心点的距离，mm；
ν 为岩石泊松比；E为岩石弹性模量，GPa。将流体压力代入式（9），计算裂缝宽度解析解，并与不同网格方案下裂缝宽度对比，结果如图4所示。
模拟结果表明，裂缝宽度的数值解与解析解吻合程度较好，自适应网格与全局细分网格的模拟结果相接近。自适应网格方法节点少，算力消耗低，计算相场裂缝模型时有显著优势。
2. 地质与工程因素对裂缝扩展的影响
2.1 温度
温度场的变化会造成岩石体积变化，在系统内部互相约束，从而产生热应力[28–29]。为研究温度对水力裂缝起裂与扩展的影响，设置初始地层温度分别为20，100和200 ℃，压裂液注入温度均为20 ℃，排量均为20 mL/min，边界应力分别为10和7 MPa，裂缝长度变化趋势如图5所示。
![]()
图 5 不同地层温度下干热岩储层水力压裂缝长变化趋势Figure 5. Variation trend of hydraulic fracture length in HDR reservoir with different formation temperature模拟过程中，由于压裂液温度低于地层温度，热应力呈现为拉应力，且地层温度越高，热应力越大；由于水力裂缝以拉伸缝为主，热应力在裂缝扩展过程中起辅助作用。压裂液排量相同时，地层温度越高，裂缝扩展速度越快，起裂压力越低。真实干热岩储层中，岩石不同矿物热膨胀时存在各向异性，热应力方向与大小更为复杂，更容易产生局部应力集中；若地层中微裂缝发育，热应力易激活并诱导微裂缝起裂，提高水力裂缝复杂程度，且热应力越大，微裂缝激活程度越高。
2.2 压裂液排量
为探究不同压裂液排量对水力裂缝扩展的影响，设置压裂液排量分别为10，20和30 mL/min，地层初始温度均为200 ℃，压裂液温度均为20 ℃，边界应力分别为10和7 MPa，注液30 s后，模拟结果如图6所示。
![]()
图 6 不同压裂液排量下干热岩储层水力压裂模拟结果Figure 6. Simulation results of hydraulic fracturing in HDR reservoir under different pumping rates of fracturing fluids从图6可以看出，压裂液排量为10 mL/min时，缝长为11.63 cm；压裂液排量为20 mL/min时，缝长为20.26 cm；压裂液排量为30 mL/min时，缝长为27.48 cm。研究表明，随着压裂液排量增加，液体憋压速度更快，裂缝周边岩石温度降低更明显。干热岩储层水力裂缝扩展受流体压力和热应力的共同驱动，裂缝扩展速度更快，裂缝宽度也随之增大，更有利于后续循环介质的注入与热交换。
2.3 天然裂缝夹角
为探究天然裂缝对水力裂缝扩展的影响，设置天然裂缝与水平最大主应力方向的夹角分别为90°，45°和20°，地层初始温度均为200 ℃，压裂液温度均为20 ℃，压裂液排量为20 mL/min，边界应力分别为10和7 MPa，模拟结果如图7所示。
![]()
图 7 不同天然裂缝夹角下干热岩储层水力压裂模拟结果Figure 7. Simulation results of hydraulic fracturing in HDR reservoir with different angles of natural fractures水力压裂裂缝起裂后沿水平最大主应力方向扩展，遇到天然裂缝后发生转向并沟通天然裂缝。天然裂缝与水平最大主应力方向夹角为90°时，水力裂缝转向不明显，沟通天然裂缝后在天然裂缝两侧同时起裂并扩展，扩展角度与最大水平地应力夹角约为45°，2条裂缝呈分离趋势。岩石温度在水力压裂裂缝与天然裂缝沟通处降幅较大，主要原因为压裂液在该区域流动情况更复杂，传热效率高。
天然裂缝与水平最大主应力方向夹角为45°时，水力裂缝转向效果较明显，沟通天然裂缝后在天然裂缝两侧起裂并扩展，左右两侧扩展速度差异明显，左侧扩展速度更快，这是因为该侧与水平最大主应力方向夹角较小，扩展速度更快。另外，岩石温度在水力裂缝与天然裂缝右侧夹角处降低幅度较大，主要是该侧天然裂缝与水力压裂裂缝距离较近，传热效率更高。
天然裂缝与水平最大主应力方向夹角为20°时，天然裂缝诱导水力压裂裂缝产生转向现象。水力压裂裂缝沟通天然裂缝后在天然裂缝两侧起裂并扩展，左右两侧扩展速度差异明显，左侧扩展速度更快，裂缝沿水平最大主应力方向继续扩展。岩石温度在水力压裂裂缝与天然裂缝右侧夹角处降低幅度较大。
基于上述模拟结果可知，水力压裂裂缝遇到天然裂缝后易沟通天然裂缝，天然裂缝与水平最大主应力夹角越小，转向现象越明显。水力压裂裂缝与天然裂缝沟通位置流动情况复杂，传热效率高。水力压裂裂缝再次起裂方向受到地应力与天然裂缝方向的共同作用，在水平最大主应力方向夹角较小的方向扩展迅速，在水平最大主应力方向夹角较大的方向扩展较慢，且该侧岩石降温幅度更大。
2.4 多天然裂缝干扰下裂缝扩展
为探究多天然裂缝干扰下干热岩储层水力裂缝扩展力学行为，设置不同地应力与排量进行模拟。单侧天然裂缝条数为2条，天然裂缝与水平最大主应力方向夹角为45°，地层初始温度为200 ℃，压裂液温度为20 ℃，压裂液排量分别为20 和50 mL/min，边界应力差分别为3和8 MPa，模拟结果如图8所示。
![]()
图 8 不同天然裂缝密度下干热岩储层水力压裂模拟结果Figure 8. Simulation results of hydraulic fracturing in HDR reservoir with different densities of natural fractures地层中存在多条天然裂缝时，水力压裂裂缝优先与距离最近的天然裂缝发生相互作用。地应力差为3 MPa、压裂液排量为20 mL/min时，水力压裂裂缝在靠近天然裂缝时轻微转向并开启初条天然裂缝，随后在天然裂缝尖端二次起裂；地应力差为8 MPa、压裂液排量为20 mL/min时，水力压裂裂缝靠近天然裂缝无明显转向现象，开启初条天然裂缝后在裂缝尖端二次起裂；地应力差为8 MPa、压裂液排量增加至50 mL/min时，水力压裂裂缝贯穿初条天然裂缝并继续扩展，由于水力压裂裂缝贯穿天然裂缝后会消耗水力能量，因此仅能贯穿单条天然裂缝。
3. 结 论
1）基于相场法建立了流−固−热多场耦合裂缝扩展数值模型，根据能量最小化原理计算裂缝形态，无需额外添加裂缝扩展准则，模型搭建简单，可用于干热岩储层水力裂缝扩展模拟研究。
2）干热岩储层水力压裂过程中，热应力体现为拉应力，有助于增加裂缝开度，在裂缝扩展过程中起辅助作用。压裂液与地层温度差越高，排量越大，热应力影响效果越显著；若地层中微裂缝发育，热应力易激活并诱导微裂缝起裂，提高水力压裂裂缝的复杂程度。
3）地层中存在多条天然裂缝时，水力压裂裂缝优先与距离最近的天然裂缝发生相互作用。水力压裂裂缝遇到天然裂缝后易沟通天然裂缝，天然裂缝与水平最大主应力方向夹角越小，地应力差越小，转向现象越明显。水力压裂裂缝再次起裂时沿天然裂缝方向起裂，且向水平最大主应力方向偏转。
-
[1] 邓金根,郑勋,闫新江,等.深水气田防砂方案优选试验研究[J].断块油气田,2012,19(3):382-385. Deng Jingen,Zheng Xun,Yan Xinjiang,et al.Experimental study on optimization of sand control scheme in deep-water gas field[J].Fault-Block Oil Gas Field,2012,19(3):382-385. [2] 赵益忠,孙磉礅,高爱花,等.稠油油藏蒸汽吞吐井长效防砂技术[J].石油钻探技术,2014,42(3):90-94. Zhao Yizhong,Sun Sangdun,Gao Aihua,et al.Long-term sand control technology for multiple round steam huff and puff wells in heavy oil reservoirs[J].Petroleum Drilling Techniques,2014,42(3):90-94. [3] 匡韶华,石磊,于丽宏,等.防砂筛管测试技术现状及发展探讨[J].石油矿场机械,2013,42(4):83-88. Kuang Shaohua,Shi Lei,Yu Lihong,et al.Discussion on state and development of testing methods for sand control screen[J].Oil Field Equipment, 2013,42(4):83-88. [4] 王爱国.金属棉优质筛管适度防砂实验研究及应用[J].断块油气田,2013,20(4):535-538. Wang Aiguo.Experimental study and application of proper sand control of metallic cotton sand screen pipe[J].Fault-Block Oil Gas Field,2013,20(4):535-538. [5] 陈宇,邓金根,何宝生,等.砾石充填井工业砾石尺寸优选[J].石油钻探技术,2011,39(3):92-95. Chen Yu,Deng Jingen,He Baosheng,et al.Optimization of industrial gravel size in gravel packed wells[J].Petroleum Drilling Techniques, 2011,39(3):92-95. [6] 韩宗正.金属塞充填射孔防砂联作技术研究与应用[J].石油钻探技术,2013,41(3):119-122. Han Zongzheng.Research and application of combined perforation and sand control with metal plug[J].Petroleum Drilling Techniques,2013,41(3):119-122. [7] Larsen O S, Fjellstad V, Mathisen A M,et al.New sand retention test setup exhibits no plugging tendencies with various screen types using nonuniform test sand[R].SPE 151346,2012.
[8] 施进,李鹏,贾江鸿.筛网式滤砂管挡砂效果室内试验[J].石油钻探技术,2013,41(3):104-108. Shi Jin,Li Peng,Jia Jianghong.Laboratory tesing of sand control effect for mesh type screen[J].Petroleum Drilling Techniques, 2013,41(3):104-108. [9] 董长银,贾碧霞,刘春苗,等.机械防砂筛管挡砂介质堵塞机制及堵塞规律试验[J].中国石油大学学报:自然科学版,2011,35(5):82-88. Dong Changyin, Jia Bixia, Liu Chunmiao,et al.Blocking mechanism and blocking laws experiments of sandretention media in mechanical screens[J].Journal of China University of Petroleum:Edition of Natural Science,2011,35(5):82-88. [10] 徐新邦,刘培生,崔光,等.泡沫金属力学性能研究的分析概述[J].金属功能材料,2012,19(6):46-50. Xu Xinbang,Liu Peisheng,Cui Guang,et al.Analysis into study on mechanical properties of metal foams[J].Metallic Functional Materials,2012,19(6):46-50. [11] 左孝青,孙加林.泡沫金属制备技术研究进展[J].材料科学与工程学报,2004,22(3):452-456. Zuo Xiaoqing,Sun Jialin.Review on foam metal manufacture techniques[J].Journal of Materials Science Engineering,2004,22(3):452-456. [12] 邓金根,李萍,王利华,等.渤海湾地区适度防砂技术防砂方式优选[J].石油钻采工艺,2011,33(1):98-101. Deng Jingen, Li Ping, Wang Lihua,et al.The optimization of sand control method for moderate sand control technique application in Bohai Bay[J].Oil Drilling Production Technology,2011,33(1):98-101 [13] 邓金根,李萍,周建良,等.中国海上疏松砂岩适度出砂井防砂方式优选[J].石油学报,2012,33(4):676-679. Deng Jingen,Li Ping,Zhou Jianliang,et al.Sand control optimization applied to moderately sanding wells in offshore loose sandstone reservoirs[J].Acta Petrolei Sinica, 2012,33(4):676-679. [14] 付全荣,张铱鈖,段滋华,等.多孔泡沫金属及其在化工设备中的应用[J].化工机械,2010,37(6):805-809. Fu Quanrong,Zhang Yifen,Duan Zihua,et al.Applications of porous foam metals in chemical equipment[J].Chemical Equipments,2010,37(6):805-809. [15] 吴向红,何伶,方宏长.比采油指数曲线的分析和应用[J].石油勘探与开发,2007,34(6):745-749. Wu Xianghong,He Ling,Fang Hongchang.Analysis and application of specific productivity index[J].Petroleum Exploration and Development,2007,34(6):745-749. [16] 匡韶华,柳燕丽.油井滤砂管挡砂精度测试方法探讨[J].重庆科技学院学报:自然科学版,2012,14(6):61-64. Kuang Shaohua,Liu Yanli.Discussion on testing methods for size of aperture of sand control pipe[J].Journal of Chongqing University of Science Technology:Natural Science Edition, 2012,14(6):61-64. [17] Rajesh A Chanpura,Richard M hodge,Jamie Stuart Andrews,et al.State of the art screen selection for standalone screen applications[R].SPE 127931,2010.
-
期刊类型引用(4)
1. 陈德春，张文宣，阳成，常峰，邴绍强，张鹏，王亮亮，马硕. 基于井口流体温度的含水率计算模型. 断块油气田. 2025(03): 514-521 . 
百度学术
2. 魏锋，陈现，王迪，夏瑜. 海上少井条件下含水率计算方法研究及应用. 海洋石油. 2022(04): 63-66 . 百度学术
3. 王谦，谭茂金，石玉江，李高仁，程相志，罗伟平. 径向基函数神经网络法致密砂岩储层相对渗透率预测与含水率计算. 石油地球物理勘探. 2020(04): 864-872+704 . 百度学术
4. 王谦，石玉江，谭茂金，李高仁. 基于孔隙结构分类的致密砂岩含水率计算模型——以鄂尔多斯盆地陇东西部延长组长8_1储层为例. 石油物探. 2019(05): 669-680 . 百度学术
其他类型引用(5)
下载:
计量
- 文章访问数: 3244
- HTML全文浏览量: 68
- PDF下载量: 3935
- 被引次数: 9
