CN115618767A - 一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 - Google Patents
一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 Download PDFInfo
- Publication number
- CN115618767A CN115618767A CN202211411593.1A CN202211411593A CN115618767A CN 115618767 A CN115618767 A CN 115618767A CN 202211411593 A CN202211411593 A CN 202211411593A CN 115618767 A CN115618767 A CN 115618767A
- Authority
- CN
- China
- Prior art keywords
- stress
- perforation
- azimuth angle
- rock
- mpa
- 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.)
- Pending
Links
Images
Classifications
-
- 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
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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
技术领域
本发明涉及一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,属于油气增产改造领域。
背景技术
和常规压裂技术相比,连续油管拖动压裂技术具有较高的储层改造准确度与作业效率,其作业时间为常规压裂技术的一半以上,有效地改善作业的效果,在大规模的低渗透油藏中的应用频率愈来愈高。同时连续油管拖动压裂技术具有较高储层改造准确度和效率,能够精准定位裂缝起裂位置,单缝规模可控,有利于实施针对性的储层改造措施,便于地质工程一体化压裂改造。
目前,相关学者对影响连续油管拖动压裂的水力喷射效果展开了实验与理论研究。付钢旦、夏强、谢刚儒、张晶等人通过室内试验研究了不同喷砂射孔参数下喷射速度、射孔时间、喷嘴尺寸以及流体性质的优选,表明了影响连续油管水力喷砂射孔能力的因素主要有喷嘴参数、磨料类型、磨料粒径和浓度、射流压力、施工排量、围压以及喷射时间等,其中喷嘴参数、磨料类型、射流压力、排量、围压是关键因素。欧阳梦迪、王博学、李朝阳、高荣兴、王尊策等人通过数值模拟仿真与计算流体力学的方法,结合CFD软件分析了射流流体参数、磨料参数、喷嘴参数以及岩石力学参数对水力喷射效果的影响。
综上所述,实验和理论研究多以喷砂射孔参数对喷射效果的影响为目标,对射孔时地层破裂压力的影响尚未开展研究。而连续油管拖动压裂作业时,其地层破裂压力要大于常规压裂作业的破裂压力,同时由于喷砂射孔数较少,射孔孔眼位置可能不在主应力位置上,进而影响地层破裂压力的大小。
因此,亟需建立一种考虑应力叠置下连续油管拖动压裂射孔方位角优化方法,将有助于提高连续油管拖动压裂设计的科学性和针对性,进一步提高低渗透储层的开发潜力。
发明内容
为了克服现有技术中的问题,本发明提供一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,该方法能够便捷地分析连续油管压裂射孔方位角对地层破裂压力的影响,并为连续油管拖动压裂执行的参数设计和优化提供了一些见解。
本发明解决上述技术问题所提供的技术方案是:一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,包括以下步骤:
步骤S1、获取储层岩石力学参数、储层物性参数、热力学参数以及井筒参数;
步骤S2、利用步骤S1获得参数,根据弹性力学理论,建立考虑原地应力、井筒内压、压裂液滤失、胶结套管以及围岩温度变化导致的应力叠置影响下的射孔孔道的总应力分布模型,并计算得到射孔井筒周围应力分布;
步骤S3、根据射孔井筒周围应力分布,并利用弹性力学方法和二分法计算得到储层岩石的最大拉应力;
步骤S4、将储层岩石的最大拉应力代入破裂准则,计算得到不同射孔方位角所对应喷砂射孔时的地层破裂压力。
进一步的技术方案是,所述步骤S1中的储层力学参数包括最大水平主应力、最小水平主应力、垂向应力、孔隙压力、岩石抗张强度,所述储层物性参数包括孔隙度、杨氏模量、泊松比、渗透性系数、Biot多孔弹性系数、Haimson修正系数,所述热力学参数包括线性热膨胀系数、岩石温度差,所述井筒参数包括井筒倾斜角、方位角、套管外径、内径。
进一步的技术方案是,所述射孔孔道的总应力分布模型为:
式中:σs、σθp、σzz为射孔孔眼中的径向、周向和轴向应力,MPa;τsθ、τzzθ、τszz为射孔孔眼中剪应力的三个坐标分量,MPa;θp为射孔的周向角度,°;rhs为射孔孔眼半径,m;s为射孔孔眼平面内某一点到孔眼轴向的距离,m;Pp为孔隙压力,MPa。
进一步的技术方案是,所述步骤S2中射孔孔道的总应力分布的计算过程如下:
步骤S21、根据弹性力学理论,计算出原地应力引起的应力场分布;
步骤S22、计算出井筒内压引起的应力场分布;
步骤S23、计算出压裂液滤失引起的应力场分布;
步骤S24、计算出胶结套管引起的应力场分布;
步骤S25、计算出围岩温度变化引起的应力场分布;
步骤S26、结合上述计算的不同应力场变化分布,得到射孔孔道的总应力场分布。
进一步的技术方案是,所述步骤S22中的计算公式为:
式中:rw为井筒半径,m;r为井眼轴线到地层中某一点的距离,m;θ为任意径向与x轴的极角,°;σr、σθ、σz为井筒坐标中的径向、环向和轴向应力,MPa;τθz、τrθ、τrz为井筒坐标中剪应力的三个分量,MPa。
进一步的技术方案是,所述步骤S23中的计算公式为:
式中:Pw为井筒内液柱压力,MPa;c为Haimson修正系数,0.9<c<1。
进一步的技术方案是,所述步骤S24中的计算公式为:
其边界条件:
式中:α为Biot多孔弹性系数;Cr、Cb分别为岩石的骨架压缩率和容积压缩率,%;φ为岩石孔隙度,%;δ为渗透性系数;Pn(r)为半径r处的地层净应力。
进一步的技术方案是,所述步骤S25中的计算公式为:
式中:为水泥套管引起的井筒周围的径向应力和切向应力,MPa;TF为传递系数,表示传递到地层中岩石的井筒压力;vc为水泥套管的泊松比,无量纲;Ec为水泥套管的杨氏模量,MPa;Ro、Ri为分别为套管的内径和外径,m。
进一步的技术方案是,所述步骤S3中储层岩石的最大拉应力计算公式为:
式中:α为Biot多孔弹性系数;σ3为储层岩石的最大拉应力,MPa。
进一步的技术方案是,所述破裂准则为:
σmax(θp)-αPp≥σt
式中:α为Biot多孔弹性系数;Pp为孔隙压力,MPa。
本发明具有以下有益效果:本发明提出了一种较为准确的计算方法优化连续油管拖动压裂喷砂射孔方位,揭示了连续油管拖动压裂喷砂射孔方位角对破裂压力的影响规律,将有助于提高低渗透储层油气资源的开发潜力。
附图说明
图1为本发明的流程示意图;
图2为本发明射孔孔道切向应力分布示意图;
图3为本发明裂缝起裂方位与起裂角示意图;
图4为A井的地层破裂压力随射孔方位角变化曲线图;
图5为B井的地层破裂压力随射孔方位角变化曲线图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明的一种考虑应力叠置下连续油管拖动压裂射孔方位角优化方法,包括以下步骤:
步骤S1、获取储层力学参数(最大最小水平主应力、垂向应力、孔隙压力、岩石抗张强度),储层物性参数(孔隙度、杨氏模量、泊松比、渗透性系数、Biot多孔弹性系数、Haimson修正系数),热力学参数(线性热膨胀系数、岩石温度差),井筒参数(井筒倾斜角、方位角、套管外径、内径);
步骤S2、利用步骤S1获得参数,根据弹性力学理论,建立考虑原地应力、井筒内压、压裂液滤失、胶结套管以及围岩温度变化导致的应力叠置影响下的射孔孔道的总应力分布模型,并计算得到射孔井筒周围应力分布;
步骤S21、根据弹性力学理论,计算出原地应力引起的应力场分布;
式中:rw为井筒半径,m;r为井眼轴线到地层中某一点的距离,m;θ为任意径向与x轴的极角,°;σr、σθ、σz为井筒坐标中的径向、环向和轴向应力,MPa;τθz、τrθ、τrz为井筒坐标中剪应力的三个分量,MPa;
步骤S22、计算出井筒内压引起的应力场分布;
式中:Pw为井筒内液柱压力,MPa;c为Haimson修正系数,0.9<c<1。
步骤S23、计算出压裂液滤失引起的应力场分布;
其边界条件:
式中:α为Biot多孔弹性系数,α=1-βv=1–Cr/Cb;Cr、Cb分别为岩石的骨架压缩率和容积压缩率,%;φ为岩石孔隙度,%;δ为渗透性系数(地层可渗透时δ=1,不可渗透时δ=0);Pn(r)为半径r处的地层净应力,Pn(r)=P(r)-P0;
步骤S24、计算出胶结套管引起的应力场分布;
式中:为水泥套管引起的井筒周围的径向应力和切向应力,MPa;TF为传递系数,表示传递到地层中岩石的井筒压力。vc为水泥套管的泊松比,无量纲;Ec为水泥套管的杨氏模量,MPa;Ro、Ri为分别为套管的内径和外径,m;
步骤S25、计算出围岩温度变化引起的应力场分布;
其中:
式中:σx、σy、σz分别为x,y,z方向上的正应力,MPa;v、vˊ分别为三个主平面内和垂直于这三个主平面方向上的泊松比,无量纲;E、Eˊ分别为三个主平面内和垂直于这三个主平面方向上的杨氏模量,MPa;α、αˊ分别为平面内和垂直于该平面方向上岩石的热膨胀系数,1/℃;T0为地层岩石初始温度,℃;T为地层岩石受到钻井及压裂液等作用后的温度,℃;αT为岩石热膨胀系数,1/℃。
步骤S27、结合式(1)-(7)计算的不同应力场变化分布,根据图2所示的几何模型,并由式(8)计算得到射孔孔道的总应力分布;
式中:σs、σθp、σzz为射孔孔眼中的径向、周向和轴向应力,MPa;τsθ、τzzθ、τszz为射孔孔眼中剪应力的三个坐标分量,MPa;θp为射孔的周向角度,°;rhs为射孔孔眼半径,m;s为射孔孔眼平面内某一点到孔眼轴向的距离,m;Pp为孔隙压力,MPa;
步骤S3、根据射孔井筒周围应力分布,并利用弹性力学方法和二分法计算得到储层岩石的最大拉应力;
由式(8)计算得到射孔孔道的总应力分布,可推导出任一点射孔孔壁上的主应力可表达式;
根据弹性力学理论,计算岩石的最大拉应力大小为:
对于某一射孔方位角θp,裂缝起裂方位角为:
步骤S4、如图3所示,根据破裂准则,当井壁z-θ平面内最大拉伸应力大于等于岩石自身的抗拉强度时,即:
σmax(θp)-αPp≥σt (12)
满足θp即为射孔方位角为θ时井壁发生拉伸破裂时的起裂压力。
其中对于任一给定的射孔方位角可以求出其所对应的围岩破裂压力及起裂方位角,起裂方位角计算方程如式(13)所示;
得到γ1、γ2;
根据函数极值的定义,当二阶函数值小于0时,原函数才存在极大值。围岩最大拉应力及其二阶导数表达式如式(15)所示;
将γ1、γ2分别代入式(15),如果数值小于0,原函数存在极大值,则其为裂缝起裂的真实角度。
实施例1
以国内某油田连续油管拖动水力喷砂射孔压裂井的基本参数,对射孔方位角与破裂压力之间的关系进行研究,本实施例中连续油管拖动压裂井的基础参数如表1所示。
表1连续油管拖动压裂井的基本参数表
按如前所述步骤开展实施例计算,输出计算结果如图4和5所示。
A井在射孔方位角为160°或340°左右时所对应的地层破裂压力最小,B井在射孔方位角在6°或180°左右时所对应的地层破裂压力最小。随着射孔方位角变化,破裂压力降低了7~20MPa不等。由于水平井眼沿水平最小主应力方向布置,井眼受到的主应力为水平最大主应力和垂向主应力,当水平最大主应力大于垂向应力时,喷嘴要调整到水平位置上,当水平最大主应力小于垂向应力时,喷嘴要调整到垂直位置上,可有效降低起裂压力。若喷砂射孔时破裂压力不是最优值,使得套管的完整性降低,抗剪切能力减弱,不利于井筒的稳定性。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (10)
1.一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,其特征在于,包括以下步骤:
步骤S1、获取储层岩石力学参数、储层物性参数、热力学参数以及井筒参数;
步骤S2、利用步骤S1获得参数,根据弹性力学理论,建立考虑原地应力、井筒内压、压裂液滤失、胶结套管以及围岩温度变化导致的应力叠置影响下的射孔孔道的总应力分布模型,并计算得到射孔井筒周围应力分布;
步骤S3、根据射孔井筒周围应力分布,并利用弹性力学方法和二分法计算得到储层岩石的最大拉应力;
步骤S4、将储层岩石的最大拉应力代入破裂准则,计算得到不同射孔方位角所对应喷砂射孔时的地层破裂压力。
2.根据权利要求1所述的一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,其特征在于,所述步骤S1中的储层力学参数包括最大水平主应力、最小水平主应力、垂向应力、孔隙压力、岩石抗张强度,所述储层物性参数包括孔隙度、杨氏模量、泊松比、渗透性系数、Biot多孔弹性系数、Haimson修正系数,所述热力学参数包括线性热膨胀系数、岩石温度差,所述井筒参数包括井筒倾斜角、方位角、套管外径、内径。
4.根据权利要求3所述的一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,其特征在于,所述步骤S2中射孔孔道的总应力分布的计算过程如下:
步骤S21、根据弹性力学理论,计算出原地应力引起的应力场分布;
步骤S22、计算出井筒内压引起的应力场分布;
步骤S23、计算出压裂液滤失引起的应力场分布;
步骤S24、计算出胶结套管引起的应力场分布;
步骤S25、计算出围岩温度变化引起的应力场分布;
步骤S26、结合上述计算的不同应力场变化分布,得到射孔孔道的总应力场分布。
10.根据权利要求1所述的一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法,其特征在于,所述破裂准则为:
σmax(θp)-αPp≥σt
式中:α为Biot多孔弹性系数;Pp为孔隙压力,MPa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211411593.1A CN115618767A (zh) | 2022-11-11 | 2022-11-11 | 一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211411593.1A CN115618767A (zh) | 2022-11-11 | 2022-11-11 | 一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115618767A true CN115618767A (zh) | 2023-01-17 |
Family
ID=84878984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211411593.1A Pending CN115618767A (zh) | 2022-11-11 | 2022-11-11 | 一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115618767A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454561A (zh) * | 2023-12-19 | 2024-01-26 | 成都信息工程大学 | 一种连续油管在水平井眼极限延伸距离的分析方法及系统 |
-
2022
- 2022-11-11 CN CN202211411593.1A patent/CN115618767A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454561A (zh) * | 2023-12-19 | 2024-01-26 | 成都信息工程大学 | 一种连续油管在水平井眼极限延伸距离的分析方法及系统 |
CN117454561B (zh) * | 2023-12-19 | 2024-03-08 | 成都信息工程大学 | 一种连续油管在水平井眼极限延伸距离的分析方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110222477B (zh) | 维持水平井分段压裂裂缝均衡扩展的射孔参数优化方法 | |
CN106855897B (zh) | 适用于压力衰竭地层的井壁稳定的研究方法 | |
CN106869892B (zh) | 一种重复压裂井缝内暂堵起裂的判断方法 | |
CN109209350A (zh) | 一种煤层斜井套管射孔完井破裂压力的预测方法 | |
CN109162701B (zh) | 一种煤层裸眼井破裂压力预测方法 | |
CN106545336A (zh) | 考虑致密气藏渗流机理的产能计算方法 | |
Fu et al. | A novel tool to improve the rate of penetration–hydraulic-pulsed cavitating-jet generator | |
CN110580401A (zh) | 判断定向井、水平井分段多簇压裂井筒暂堵次数的方法 | |
CN107705215A (zh) | 一种页岩储层重复压裂选井选段方法 | |
CN115618767A (zh) | 一种考虑应力叠置下连续油管拖动压裂的射孔方位角优化方法 | |
CN111271042A (zh) | 一种针对天然裂缝发育储层水力压裂形成复杂缝网的排量优化方法 | |
CN113239499A (zh) | 一种页岩气井最大油嘴尺寸优选方法 | |
CN113850029A (zh) | 一种页岩气水平井密切割压裂射孔参数优化设计方法 | |
CN104847314A (zh) | 高温高压油气直井单相流射孔完井参数优化方法 | |
CN114372428B (zh) | 砂砾岩储层水平井段内多簇压裂裂缝延伸跨尺度模拟方法 | |
CN108121844B (zh) | 水力波及半径的获得方法 | |
CN116717224A (zh) | 一种针对低渗致密储层复杂缝网的压裂产能预测方法 | |
CN115952620A (zh) | 一种页岩储层防套变压裂施工方法与关键裂缝控制材料 | |
Zhang et al. | Numerical simulation on the flow field of self-propelled multi-orifices nozzle for ultra-short radius radial jet drilling | |
CN115186524A (zh) | 低渗透长缝压裂井组co2驱临界注气参数的评测方法 | |
CN111946321B (zh) | 填砂暂堵压裂的支撑剂参数设计方法 | |
CN106227903B (zh) | 底水层钻开厚度的确定方法及装置 | |
CN115522918A (zh) | 深层砂岩油藏射孔井出砂压差剖面预测方法 | |
CN115705454A (zh) | 一种基于相场法的裂缝扩展模拟压裂设计优化方法 | |
CN107916917B (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 |