CN113971378B - 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 - Google Patents
一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 Download PDFInfo
- Publication number
- CN113971378B CN113971378B CN202111257648.3A CN202111257648A CN113971378B CN 113971378 B CN113971378 B CN 113971378B CN 202111257648 A CN202111257648 A CN 202111257648A CN 113971378 B CN113971378 B CN 113971378B
- Authority
- CN
- China
- Prior art keywords
- perforation
- cluster
- fracturing
- fracture
- temporary plugging
- 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.)
- Active
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B33/00—Sealing or packing boreholes or wells
- E21B33/10—Sealing or packing boreholes or wells in the borehole
- E21B33/13—Methods or devices for cementing, for plugging holes, crevices, or the like
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B33/00—Sealing or packing boreholes or wells
- E21B33/10—Sealing or packing boreholes or wells in the borehole
- E21B33/13—Methods or devices for cementing, for plugging holes, crevices, or the like
- E21B33/138—Plastering the borehole wall; Injecting into the formation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/261—Separate steps of (1) cementing, plugging or consolidating and (2) fracturing or attacking the formation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开了一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,包括:获取储层地质参数、压裂施工参数、射孔布簇参数、暂堵转向参数;计算各条水力裂缝延伸长度、裂缝延伸高度、裂缝延伸开度以及裂缝缝内压力;计算水平井筒各簇射孔流量分配、暂堵球堵塞射孔概率、射孔簇暂堵球堵塞数量以及各簇射孔剩余的有效孔眼数量;计算各射孔孔眼流动速度与射孔孔眼直径;进行下一时步的计算,直到压裂时间结束,确定暂堵球粒径。本发明根据泵入暂堵剂时刻的孔眼直径大小,优选合适的暂堵球粒径,实现射孔孔眼的有效封堵,解决了深层页岩气水平井缝口转向压裂缺少暂堵球粒径定量优选手段的问题,进一步提高页岩气压裂设计的科学性和针对性。
Description
技术领域
本发明涉及一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,属于页岩气开发技术领域。
背景技术
深层页岩气水平井压裂具有射孔簇数多,簇间距小等特点,可在一定程度上扩大缝网体积,并提高压裂缝网复杂程度。但是,由于压裂时各条水力裂缝之间距离较近,延伸过程中存在显著的应力干扰效应,导致各条水力裂缝相互竞争,延伸速度快慢不一,延伸长度存在较大差异,严重影响深层页岩气水平井压裂增产效果。为此,深层页岩气水平井压裂时通常需要结合缝口转向工艺,即在压裂期间泵入暂堵球,利用暂堵球对延伸较快裂缝的射孔孔眼进行封堵,减小其缝口压裂液流入流量,减慢其延伸速度,最终实现各条裂缝均匀延伸的目的。
目前,国内外学者针对水平井缝口转向压裂开展了相关研究。龚蔚(2017)、丁邦春(2018)、夏海帮(2020)分别在高石梯-磨溪区块、苏里格气田、川南页岩气区块开展了水平井缝口转向压裂工艺的矿场应用,并结合测井、微地震监测等手段证实了一部分暂堵球可对射孔孔眼进行封堵;金智荣(2019)、吕瑞华(2020)、刘明明(2020)、韩慧芬(2021)建立了暂堵球在井筒内运移与坐封射孔孔眼数学模型,可对泵入暂堵球时的压裂排量进行优化设计;Feng Zhang(2020)、Ming Chen(2020)、周彤(2020)、陈钊(2021)将暂堵球运移坐封模型与页岩气缝网压裂裂缝延伸模型结合,定量分析了暂堵球数量、暂堵时机、暂堵次数对压裂效果的影响,进而形成了一套暂堵球参数优化方法。
综上所述,目前水平井缝口转向压裂相关研究主要集中于暂堵球封堵各簇射孔孔眼情况下的水力裂缝延伸模拟,可用于暂堵球数量、暂堵时机、暂堵次数、暂堵排量等参数进行优化设计,但均未考虑深层页岩气水平井缝口转向压裂过程中支撑剂持续对射孔冲蚀导致其孔眼直径逐渐变化的行为,因此不能确定暂堵时刻各簇射孔孔眼直径大小,从而无法优选粒径合适的暂堵球对射孔孔眼进行有效封堵,使得深层页岩气水平井缝口转向压裂设计和工艺优化仍然存在一定盲目性。
因此,亟需建立一种适用于深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,将有助于进一步提高深层页岩气水平井缝口转向压裂设计的科学性和针对性,改善深层页岩气压裂增产效果。
发明内容
为了克服现有技术中的问题,本发明提供一种解决深层页岩气水平井缝口转向压裂过程中考虑各簇射孔孔眼直径动态变化的情况下,缺少暂堵球粒径定量优选手段问题的深层页岩气水平井缝口转向压裂暂堵球粒径优选方法。
本发明解决上述技术问题所提供的技术方案是:一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,包括以下步骤:
步骤一、获取储层地质参数、压裂施工参数、射孔布簇参数、暂堵转向参数;
步骤二、运用流固耦合理论建立深层页岩气水平井缝口转向压裂多簇裂缝延伸模型,并计算各条水力裂缝延伸长度、裂缝延伸高度、裂缝延伸开度以及裂缝缝内压力;
步骤三、运用流体力学理论建立深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型,计算水平井筒各簇射孔流量分配、暂堵球堵塞射孔概率、射孔簇暂堵球堵塞数量以及各簇射孔剩余的有效孔眼数量;
步骤四、运用摩擦力学理论建立深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型,计算各射孔孔眼流动速度与射孔孔眼直径;
步骤五、将步骤三中计算得到的水平井筒各簇射孔流量分配、各簇射孔剩余的有效孔眼数量和步骤四中计算的得到的射孔孔眼直径作为下一时步的初始条件;从而重复步骤二而至步骤四进行下一时步的计算,直到压裂时间结束并分别绘制深层页岩气水平井缝口转向压裂多簇裂缝延伸图、各簇射孔有效孔眼数量随时间变化图、各簇射孔孔眼直径随时间变化图;
步骤六、根据各簇射孔孔眼直径随时间变化图确定暂堵球粒径。
进一步的技术方案是,所述层地质参数包括地层最小水平主应力、地层岩石断裂韧性、地层岩石杨氏模量、地层岩石泊松比;
所述压裂施工参数包括压裂排量、压裂液量、压裂时间、压裂液体粘度、压裂液滤失系数、压裂液密度、压裂注入支撑剂平均浓度;
所述射孔布簇参数包括射孔簇数、单簇射孔孔眼数量、射孔孔眼直径;
所述暂堵转向参数包括暂堵球密度、泵入暂堵球总数量、泵入暂堵球时间、射孔转向流动系数。
进一步的技术方案是,所述深层页岩气水平井缝口转向压裂多簇裂缝延伸模型包括单条裂缝内流体流动方程、单条裂缝延伸物质平衡方程、裂缝高度方程、裂缝开度方程、裂缝延伸边界条件与初始条件方程。
进一步的技术方案是,所述单条裂缝内流体流动方程为:
式中:pi(s,t)—t时刻第i条裂缝内s位置处流体压力,Pa;qi(s,t)—t时刻第i条裂缝内s位置处流量,m3/s;hi(s,t)—t时刻第i条裂缝s位置处高度,m;wi(s,t)—t时刻第i条裂缝s位置处开度,m;si—第i条裂缝长度方向坐标,m;μ—压裂液体粘度,Pa·s;
所述单条裂缝延伸物质平衡方程为:
其中:
式中:qL,i(s,t)—压裂液滤失速度,m/s;CL—压裂液滤失系数,m/s0.5;t—压裂时间,s;τi—第i条裂缝s位置开始滤失时间,s;
所述裂缝高度方程为:
式中:KIC—地层岩石断裂韧性,Pa·m0.5;σh—地层最小水平主应力,Pa;
所述裂缝开度方程为:
式中:ν—泊松比,无量纲;E—杨氏模量,Pa;
所述裂缝延伸边界条件与初始条件方程为:
式中:Qi(t)—t时刻第i条裂缝缝口流量,m3/min;Li(t)—t时刻第i条裂缝半长,m;QT—压裂排量,m3/min。
进一步的技术方案是,所述深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型包括水平井筒各射孔簇流量分配方程、射孔簇内各孔眼流量分配方程、暂堵球堵塞射孔概率方程、射孔簇暂堵球堵塞数量方程。
进一步的技术方案是,所述所述水平井筒各射孔簇流量分配方程为:
式中:pheel—水平井跟端压力,Pa;pfi,i—第i条裂缝首个单元内压力,Pa;Δppf,i—第i簇射孔孔眼处的摩阻压降,Pa;Δpw,k—第k簇射孔与第k-1簇射孔之间水平井段沿程压降,Pa;npf,i—第i簇射孔孔眼数量,个;dpf,i—第i簇射孔孔眼直径,m;αpf—孔眼流量系数,一般取0.85,无量纲;ρ—压裂液密度,kg/m3;Lw,k—第k簇射孔与第k-1簇射孔之间水平井段长度,m;qw,k—第k簇射孔与第k-1簇射孔之间水平井段内流量,m3/s;qcl,i—i号簇流量,即第i条裂缝缝口流量,m3/s;Ncl—射孔簇数,个;dw—水平井筒直径,m。
射孔簇内各孔眼流量分配方程:
式中:qpf,j|j∈i—i射孔簇内的j号射孔孔眼流量,m3/s;Npf,i—i号簇的射孔总数量,个;
暂堵球堵塞射孔概率方程:
式中:fblock,j—j号射孔被暂堵球封堵概率,无量纲;ξdivert,j—j号射孔转向流动系数,表征暂堵球在射孔处转向的难易程度,取值0~1,无量纲;qw,j—j号射孔下游处水平井段内流量,m3/s;ρdivert,—暂堵球密度,kg/m3;ρfluid,—压裂液密度,kg/m3;
射孔簇暂堵球堵塞数量方程:
式中:Mdivert,i—水平井筒内i号簇位置处剩余的暂堵球数量,个;Mblock,i—i号簇被封堵的射孔数量,个;Mtotal—泵入暂堵球总数量,个;Meffective,i—i号簇位置处剩余的暂堵球数量,个。
进一步的技术方案是,所述深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型包括射孔孔眼流动速度方程、射孔孔眼直径方程。
进一步的技术方案是,所述射孔孔眼流动速度方程为:
式中:vi—i号射孔孔眼流动速度,m/s;dpf,j—j号簇的射孔孔眼直径,m;
射孔孔眼直径方程为:
式中:C—压裂注入支撑剂平均浓度,kg/m3。
进一步的技术方案是,所述步骤六中:根据各簇射孔孔眼直径随时间变化图确定各簇最大射孔孔眼直径,再优选暂堵球粒径大于各簇最大射孔孔眼直径的1.2倍。
本发明具有以下有益效果:
1、该方法专门针对缝口转向压裂工艺特征,建立了多簇裂缝延伸模型、暂堵球封堵射孔模型、射孔孔眼动态磨蚀模型,从而提出了一种适用于深层页岩气水平井缝口转向压裂的暂堵球粒径优选方法;
2、由于该方法综合考虑了深层页岩气水平井缝口转向压裂过程中支撑剂持续对射孔冲蚀导致其孔眼直径逐渐变化的行为,为此,可根据泵入暂堵剂时刻的孔眼直径大小,优选合适的暂堵球粒径,实现射孔孔眼的有效封堵,解决了深层页岩气水平井缝口转向压裂缺少暂堵球粒径定量优选手段的问题。
附图说明
图1为本发明的计算流程示意图;
图2为本发明某深层页岩气水平井缝口转向压裂多簇裂缝延伸图;
图3为发明某深层页岩气水平井缝口转向压裂期间各簇射孔有效孔眼数量随时间变化图;
图4为发明某深层页岩气水平井缝口转向压裂期间各簇射孔孔眼直径随时间变化图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
除非另外定义,本公开使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同,而不排除其他元件或者物件。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
如图1所示,本发明的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,包括以下步骤:
步骤一、获取储层地质参数、压裂施工参数、射孔布簇参数、暂堵转向参数;
其中所述储层地质参数包括地层最小水平主应力、地层岩石断裂韧性、地层岩石杨氏模量、地层岩石泊松比;
压裂施工参数包括压裂排量、压裂液量、压裂时间、压裂液体粘度、压裂液滤失系数、压裂液密度、压裂注入支撑剂平均浓度;
射孔布簇参数包括射孔簇数、单簇射孔孔眼数量、射孔孔眼直径;
暂堵转向参数包括暂堵球密度、泵入暂堵球总数量、泵入暂堵球时间、射孔转向流动系数;
步骤二、运用流固耦合理论建立深层页岩气水平井缝口转向压裂多簇裂缝延伸模型,并计算各条水力裂缝延伸长度、裂缝延伸高度、裂缝延伸开度以及裂缝缝内压力;
其中深层页岩气水平井缝口转向压裂多簇裂缝延伸模型包括单条裂缝内流体流动方程、单条裂缝延伸物质平衡方程、裂缝高度方程、裂缝开度方程、裂缝延伸边界条件与初始条件方程;
所述单条裂缝内流体流动方程为:
式中:pi(s,t)—t时刻第i条裂缝内s位置处流体压力,Pa;qi(s,t)—t时刻第i条裂缝内s位置处流量,m3/s;hi(s,t)—t时刻第i条裂缝s位置处高度,m;wi(s,t)—t时刻第i条裂缝s位置处开度,m;si—第i条裂缝长度方向坐标,m;μ—压裂液体粘度,Pa·s;
所述单条裂缝延伸物质平衡方程为:
其中:
式中:
qL,i(s,t)—压裂液滤失速度,m/s;CL—压裂液滤失系数,m/s0.5;t—压裂时间,s;τi—第i条裂缝s位置开始滤失时间,s;
所述裂缝高度方程为:
式中:KIC—地层岩石断裂韧性,Pa·m0.5;σh—地层最小水平主应力,Pa;
所述裂缝开度方程为:
式中:ν—泊松比,无量纲;E—杨氏模量,Pa;
所述裂缝延伸边界条件与初始条件方程为:
式中:Qi(t)—t时刻第i条裂缝缝口流量,m3/min;Li(t)—t时刻第i条裂缝半长,m;QT—压裂排量,m3/min;
步骤三、运用流体力学理论建立深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型,计算水平井筒各簇射孔流量分配、暂堵球堵塞射孔概率、射孔簇暂堵球堵塞数量以及各簇射孔剩余的有效孔眼数量;
其中深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型包括水平井筒各射孔簇流量分配方程、射孔簇内各孔眼流量分配方程、暂堵球堵塞射孔概率方程、射孔簇暂堵球堵塞数量方程;
所述水平井筒各射孔簇流量分配方程为:
式中:pheel—水平井跟端压力,Pa;pfi,i—第i条裂缝首个单元内压力,Pa;Δppf,i—第i簇射孔孔眼处的摩阻压降,Pa;Δpw,k—第k簇射孔与第k-1簇射孔之间水平井段沿程压降,Pa;npf,i—第i簇射孔孔眼数量,个;dpf,i—第i簇射孔孔眼直径,m;αpf—孔眼流量系数,一般取0.85,无量纲;ρ—压裂液密度,kg/m3;Lw,k—第k簇射孔与第k-1簇射孔之间水平井段长度,m;qw,k—第k簇射孔与第k-1簇射孔之间水平井段内流量,m3/s;qcl,i—i号簇流量,即第i条裂缝缝口流量,m3/s;Ncl—射孔簇数,个;dw—水平井筒直径,m。
射孔簇内各孔眼流量分配方程:
式中:qpf,j|j∈i—i射孔簇内的j号射孔孔眼流量,m3/s;Npf,i—i号簇的射孔总数量,个。
暂堵球堵塞射孔概率方程:
式中:fblock,j—j号射孔被暂堵球封堵概率,无量纲;ξdivert,j—j号射孔转向流动系数,表征暂堵球在射孔处转向的难易程度,取值0~1,无量纲;qw,j—j号射孔下游处水平井段内流量,m3/s;ρdivert,—暂堵球密度,kg/m3;ρfluid,—压裂液密度,kg/m3;
射孔簇暂堵球堵塞数量方程:
式中:Mdivert,i—水平井筒内i号簇位置处剩余的暂堵球数量,个;Mblock,i—i号簇被封堵的射孔数量,个;Mtotal—泵入暂堵球总数量,个;Meffective,i—i号簇位置处剩余的暂堵球数量,个;
步骤四、运用摩擦力学理论建立深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型,计算各射孔孔眼流动速度与射孔孔眼直径;
其中深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型包括射孔孔眼流动速度方程、射孔孔眼直径方程;
所述射孔孔眼流动速度方程为:
式中:vi—i号射孔孔眼流动速度,m/s;dpf,j—j号簇的射孔孔眼直径,m;
射孔孔眼直径方程为:
式中:C—压裂注入支撑剂平均浓度,kg/m3;
步骤五、将步骤三中计算得到的水平井筒各簇射孔流量分配、各簇射孔剩余的有效孔眼数量和步骤四中计算的得到的射孔孔眼直径作为下一时步的初始条件;从而重复步骤二而至步骤四进行下一时步的计算,直到压裂时间结束并分别绘制深层页岩气水平井缝口转向压裂多簇裂缝延伸图、各簇射孔有效孔眼数量随时间变化图、各簇射孔孔眼直径随时间变化图;
步骤六、根据各簇射孔孔眼直径随时间变化图,确定泵入暂堵球时刻各簇射孔孔眼直径,优选暂堵球粒径,保证暂堵球粒径大于各簇最大射孔孔眼直径的1.2倍。
实施例
已知某一典型深层页岩气井矿场实际参数如表1所示,根据图1流程开展实例计算:
表1某深层页岩气井矿场实际参数表
首先,运用流固耦合理论,建立深层页岩气水平井缝口转向压裂多簇裂缝延伸模型,又结合以下公式计算实施例1中各条水力裂缝延伸情况:
具体步骤为:利用方程式(1)~(6),利用有限差分方法,计算深层页岩气水平井缝口转向压裂过程中各条水力裂缝延伸长度、高度、开度,以及裂缝内压力。
随后,运用流体力学理论,建立深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型,又结合以下公式计算实施例1中各簇射孔流量与剩余的有效孔眼数量变化情况:
具体步骤为:利用方程式(7)~(11),利用牛顿迭代方法,计算深层页岩气水平井缝口转向压裂过程中水平井筒各簇射孔流量分配、暂堵球堵塞射孔概率、射孔簇暂堵球堵塞数量,以及各簇射孔剩余的有效孔眼数量,更新各簇射孔孔眼数量;将水平井筒各簇射孔流量带入方程式(6)中,更新各条裂缝缝口流量,作为下一时步内多簇裂缝延伸模型计算的初始条件。
随后,运用摩擦力学理论建立深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型,又结合以下公式计算实施例1中各簇射孔孔眼直径变化情况:
具体步骤为:利用方程式(12)~(13),利用微分方程求解方法,计算深层页岩气水平井缝口转向压裂过程中各射孔孔眼流动速度与射孔孔眼直径,并将射孔孔眼直径带入方程式(7)中,更新各簇射孔孔眼直径,计算下一时步内各条裂缝缝口流量。
随后,运用如图1所示的本发明方法的数值计算流程框图开展实例计算,并根据计算结果,分别绘制深层页岩气水平井缝口转向压裂多簇裂缝延伸图(如图2所示)、各簇射孔有效孔眼数量随时间变化图(如图3所示)、各簇射孔孔眼直径随时间变化图(如图4所示)。
最后,根据各簇射孔孔眼直径随时间变化图,确定泵入暂堵球时刻各簇射孔孔眼直径,优选暂堵球粒径,保证暂堵球粒径大于各簇最大射孔孔眼直径的1.2倍,如表2所示。
表2泵入暂堵球时刻各簇射孔孔眼直径与暂堵球粒径优选值
基于上述计算流程和结果,该深层页岩气水平井缝口转向压裂施工过程中,可于第43分钟第一次泵入粒径为15mm的12颗暂堵球,再于第88分钟第二次泵入粒径为18mm的12颗暂堵球,可实现射孔孔眼的有效封堵,促进各条裂缝均匀延伸。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (7)
1.一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,包括以下步骤:
步骤一、获取储层地质参数、压裂施工参数、射孔布簇参数、暂堵转向参数;
步骤二、运用流固耦合理论建立深层页岩气水平井缝口转向压裂多簇裂缝延伸模型,并计算各条水力裂缝延伸长度、裂缝延伸高度、裂缝延伸开度以及裂缝缝内压力;
步骤三、运用流体力学理论建立深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型,计算水平井筒各簇射孔流量分配、暂堵球堵塞射孔概率、射孔簇暂堵球堵塞数量以及各簇射孔剩余的有效孔眼数量;
其中深层页岩气水平井缝口转向压裂暂堵球封堵射孔模型包括水平井筒各射孔簇流量分配方程、射孔簇内各孔眼流量分配方程、暂堵球堵塞射孔概率方程、射孔簇暂堵球堵塞数量方程;
所述水平井筒各射孔簇流量分配方程为:
式中:pheel—水平井跟端压力,Pa;pfi,i—第i条裂缝首个单元内压力,Pa;Δppf,i—第i簇射孔孔眼处的摩阻压降,Pa;Δpw,k—第k簇射孔与第k-1簇射孔之间水平井段沿程压降,Pa;npf,i—第i簇射孔孔眼数量,个;dpf,i—第i簇射孔孔眼直径,m;αpf—孔眼流量系数,取0.85,无量纲;ρ—压裂液密度,kg/m3;Lw,k—第k簇射孔与第k-1簇射孔之间水平井段长度,m;qw,k—第k簇射孔与第k-1簇射孔之间水平井段内流量,m3/s;qcl,i—i号簇流量,即第i条裂缝缝口流量,m3/s;Ncl—射孔簇数,个;dw—水平井筒直径,m;μ—压裂液体粘度,Pa·s;σh—地层最小水平主应力,Pa;QT—压裂排量,m3/min;
射孔簇内各孔眼流量分配方程:
式中:qpf,j|j∈i—i射孔簇内的j号射孔孔眼流量,m3/s;Npf,i—i号簇的射孔总数量,个;
暂堵球堵塞射孔概率方程:
式中:fblock,j—j号射孔被暂堵球封堵概率,无量纲;ξdivert,j—j号射孔转向流动系数,表征暂堵球在射孔处转向的难易程度,取值0~1,无量纲;qw,j—j号射孔下游处水平井段内流量,m3/s;ρdivert,—暂堵球密度,kg/m3;ρfluid,—压裂液密度,kg/m3;
射孔簇暂堵球堵塞数量方程:
式中:Mdivert,i—水平井筒内i号簇位置处剩余的暂堵球数量,个;Mblock,i—i号簇被封堵的射孔数量,个;Mtotal—泵入暂堵球总数量,个;Meffective,i—i号簇位置处剩余的暂堵球数量,个;
步骤四、运用摩擦力学理论建立深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型,计算各射孔孔眼流动速度与射孔孔眼直径;
步骤五、将步骤三中计算得到的水平井筒各簇射孔流量分配、各簇射孔剩余的有效孔眼数量和步骤四中计算的得到的射孔孔眼直径作为下一时步的初始条件;从而重复步骤二而至步骤四进行下一时步的计算,直到压裂时间结束并分别绘制深层页岩气水平井缝口转向压裂多簇裂缝延伸图、各簇射孔有效孔眼数量随时间变化图、各簇射孔孔眼直径随时间变化图;
步骤六、根据各簇射孔孔眼直径随时间变化图确定暂堵球粒径。
2.根据权利要求1所述的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,所述层地质参数包括地层最小水平主应力、地层岩石断裂韧性、地层岩石杨氏模量、地层岩石泊松比;
所述压裂施工参数包括压裂排量、压裂液量、压裂时间、压裂液体粘度、压裂液滤失系数、压裂液密度、压裂注入支撑剂平均浓度;
所述射孔布簇参数包括射孔簇数、单簇射孔孔眼数量、射孔孔眼直径;
所述暂堵转向参数包括暂堵球密度、泵入暂堵球总数量、泵入暂堵球时间、射孔转向流动系数。
3.根据权利要求1所述的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,所述深层页岩气水平井缝口转向压裂多簇裂缝延伸模型包括单条裂缝内流体流动方程、单条裂缝延伸物质平衡方程、裂缝高度方程、裂缝开度方程、裂缝延伸边界条件与初始条件方程。
4.根据权利要求3所述的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,所述单条裂缝内流体流动方程为:
式中:pi(s,t)—t时刻第i条裂缝内s位置处流体压力,Pa;qi(s,t)—t时刻第i条裂缝内s位置处流量,m3/s;hi(s,t)—t时刻第i条裂缝s位置处高度,m;wi(s,t)—t时刻第i条裂缝s位置处开度,m;si—第i条裂缝长度方向坐标,m;μ—压裂液体粘度,Pa·s;
所述单条裂缝延伸物质平衡方程为:
其中:
式中:qL,i(s,t)—压裂液滤失速度,m/s;CL—压裂液滤失系数,m/s0.5;t—压裂时间,s;τi—第i条裂缝s位置开始滤失时间,s;
所述裂缝高度方程为:
式中:KIC—地层岩石断裂韧性,Pa·m0.5;σh—地层最小水平主应力,Pa;
所述裂缝开度方程为:
式中:ν—泊松比,无量纲;E—杨氏模量,Pa;
所述裂缝延伸边界条件与初始条件方程为:
式中:Qi(t)—t时刻第i条裂缝缝口流量,m3/min;Li(t)—t时刻第i条裂缝半长,m;QT—压裂排量,m3/min。
5.根据权利要求1所述的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,所述深层页岩气水平井缝口转向压裂射孔孔眼动态磨蚀模型包括射孔孔眼流动速度方程、射孔孔眼直径方程。
7.根据权利要求1所述的一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法,其特征在于,所述步骤六中:根据各簇射孔孔眼直径随时间变化图确定各簇最大射孔孔眼直径,再优选暂堵球粒径大于各簇最大射孔孔眼直径的1.2倍。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111257648.3A CN113971378B (zh) | 2021-10-27 | 2021-10-27 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
GB2209351.2A GB2612399B (en) | 2021-10-27 | 2022-06-27 | Method for optimizing the particle size of temporary plugging balls for fracture-opening and diverting fracturing in deep shale gas horizontal well |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111257648.3A CN113971378B (zh) | 2021-10-27 | 2021-10-27 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113971378A CN113971378A (zh) | 2022-01-25 |
CN113971378B true CN113971378B (zh) | 2022-08-02 |
Family
ID=79588610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111257648.3A Active CN113971378B (zh) | 2021-10-27 | 2021-10-27 | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113971378B (zh) |
GB (1) | GB2612399B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115659736B (zh) * | 2022-10-19 | 2023-11-03 | 西南石油大学 | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 |
CN116220647A (zh) * | 2023-05-06 | 2023-06-06 | 北京石油化工学院 | 基于应力笼效应的平面密切割压裂方法、装置及存储介质 |
CN116696279B (zh) * | 2023-08-07 | 2024-01-26 | 新疆斐德莱布能源科技有限公司 | 一种长直井段储气库的暂堵作业方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109977612A (zh) * | 2019-04-19 | 2019-07-05 | 高东伟 | 一种适用于加密开发页岩气井的压裂新工艺 |
CN110905472A (zh) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | 确定基于复合暂堵体系的实时转向压裂参数的方法 |
CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
CN111639434A (zh) * | 2020-06-04 | 2020-09-08 | 华美孚泰油气增产技术服务有限责任公司 | 一种页岩气水平井投球暂堵转向压裂的优化方法 |
CN112949129A (zh) * | 2021-03-03 | 2021-06-11 | 西南石油大学 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110210144B (zh) * | 2019-06-05 | 2019-12-27 | 西南石油大学 | 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法 |
US11549046B2 (en) * | 2020-01-30 | 2023-01-10 | Kuraray Co., Ltd. | Particulate polyvinyl alcohol plugging agents |
CN111539101B (zh) * | 2020-04-17 | 2021-05-07 | 中国石油大学(华东) | 一种基于刚性颗粒堵漏剂粒径优选裂缝地层堵漏用刚性颗粒堵漏剂浓度的方法 |
CN114233261B (zh) * | 2021-12-23 | 2023-05-09 | 西南石油大学 | 一种分段压裂实现油气井均匀改造的方法 |
-
2021
- 2021-10-27 CN CN202111257648.3A patent/CN113971378B/zh active Active
-
2022
- 2022-06-27 GB GB2209351.2A patent/GB2612399B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109977612A (zh) * | 2019-04-19 | 2019-07-05 | 高东伟 | 一种适用于加密开发页岩气井的压裂新工艺 |
CN110905472A (zh) * | 2019-10-29 | 2020-03-24 | 中国石油集团川庆钻探工程有限公司 | 确定基于复合暂堵体系的实时转向压裂参数的方法 |
CN111322050A (zh) * | 2020-04-24 | 2020-06-23 | 西南石油大学 | 一种页岩水平井段内密切割暂堵压裂施工优化方法 |
CN111639434A (zh) * | 2020-06-04 | 2020-09-08 | 华美孚泰油气增产技术服务有限责任公司 | 一种页岩气水平井投球暂堵转向压裂的优化方法 |
CN112949129A (zh) * | 2021-03-03 | 2021-06-11 | 西南石油大学 | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 |
Non-Patent Citations (4)
Title |
---|
Formation of temporary plugging in acid-etched fracture with degradable diverters;YANG Chen;《Journal of Petroleum Science and Engineering》;20201231;全文 * |
涪陵焦石工区页岩气上部气层压裂技术探讨;冯德华;《江汉石油职工大学学报》;20210731;全文 * |
非均匀应力场影响下的裂缝扩展模拟及投球暂堵优化;周彤等;《天然气工业》;20200325(第03期);全文 * |
页岩气水平井投球暂堵压裂裂缝数模研究;李奎东;《天然气勘探与开发》;20210630;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113971378A (zh) | 2022-01-25 |
GB2612399B (en) | 2024-01-10 |
GB202209351D0 (en) | 2022-08-10 |
GB2612399A (en) | 2023-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113971378B (zh) | 一种深层页岩气水平井缝口转向压裂暂堵球粒径优选方法 | |
CN110056336B (zh) | 一种页岩气缝网压裂施工压力曲线自动诊断方法 | |
CN112836446B (zh) | 一种致密储层水平井分簇限流布孔优化设计方法 | |
CN112949129B (zh) | 一种深层页岩水平井压裂多簇裂缝异步起裂延伸计算方法 | |
CN113850029B (zh) | 一种页岩气水平井密切割压裂射孔参数优化设计方法 | |
CN110485986B (zh) | 碳酸盐岩缝洞型油藏深裸眼水平井分段酸压控水增产方法 | |
CN107366530A (zh) | 一种深层页岩气藏增产方法及其应用 | |
CN109800523B (zh) | 一种水力压裂暂堵条件下非正交天然裂缝开启判定方法 | |
WO2021004442A1 (zh) | 一种用于微地震解释压裂裂缝参数结果的校正方法 | |
CN110439519B (zh) | 一种基于极限限流设计的压裂方法及系统 | |
CN105386746A (zh) | 水平井水力压裂射孔方法 | |
CN111271042A (zh) | 一种针对天然裂缝发育储层水力压裂形成复杂缝网的排量优化方法 | |
CN109424348A (zh) | 水平井多簇射孔压裂裂缝条数的判断方法及系统 | |
US11692423B1 (en) | Method for realizing uniform stimulation for the oil and gas well by low-cost multi-stage fracturing | |
CN112761609A (zh) | 一种水力压裂作业中高效铺置支撑剂的优化方法 | |
CN113836753A (zh) | 一种水平井段内分簇射孔缝间暂堵转向投球优化方法 | |
CN111472730A (zh) | 大段多簇压裂的射孔方案确定方法 | |
CN113982568B (zh) | 致密油气多裂缝竞争起裂井底压力预测方法 | |
CN115659736A (zh) | 一种深层页岩气水平井缝内转向压裂缝网扩展计算方法 | |
CN111639434A (zh) | 一种页岩气水平井投球暂堵转向压裂的优化方法 | |
CN111810108B (zh) | 一种页岩气水平井压后返排油嘴动态调整系统及方法 | |
CN114372428A (zh) | 砂砾岩储层水平井段内多簇压裂裂缝延伸跨尺度模拟方法 | |
CN107605449B (zh) | 一种用于非均质储层的投球暂堵分层压裂方法 | |
CN111946321B (zh) | 填砂暂堵压裂的支撑剂参数设计方法 | |
CN112031728A (zh) | 一种水平井复合酸压方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230418 Address after: 610500, Xindu Avenue, Xindu District, Sichuan, Chengdu, 8 Patentee after: SOUTHWEST PETROLEUM University Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Address before: 610500, Xindu Avenue, Xindu District, Sichuan, Chengdu, 8 Patentee before: SOUTHWEST PETROLEUM University |
|
TR01 | Transfer of patent right |