国内外破裂压力计算方法：总结 计划 汇报 设计 可编辑

国内外破裂压力计算方法：总结 计划 汇报 设计 可编辑 破裂压力计算概述 1引言 1.1破裂压力概念 地层破裂压力(P)定义为使地层产生水力裂缝或张开原有裂缝时的井底压力,B 要实现水力加砂压裂的前提条件是应该有足够的地面泵压使井底目的层地层开裂。实际生产中通常用破裂压力梯度G(地层破裂压力P与地层深度H的比值)BB 表示破裂压力的大小,破裂压力梯度值G一般由压裂实践统计得出。地层破裂压B 力与岩石弹性性质、孔隙压力、天然裂缝发育情况以及该地区的地应力等因素有关。在压裂施工中的地层破裂压力还可以这样来理解就是裂缝即将开启而未开启时的井底压力;在压裂施工作业中，如果起泵初期压力有比较明显的降落时，那么我们就可以确定出破裂压力来这一数值可用下面这一关系式来描述:地层破裂压力=裂施工作业初期的最高套管压力+层中部的液柱压力 1.2破裂压力的获取途径 水力压裂是油气井最常用的一种增产措施，而地层破裂压力是压裂设计和施工的一项重要参数，确定该参数正确与否，将关系到能否保证压开地层等问题。 该参数的获取有两种途径:一是进行室内岩石力学实验或井场水力压裂施工;二是从测井资料中提取。目前，用测井资料估算砂泥岩剖面地层破裂压力的方法与技术较为成熟。由于碳酸盐岩地层原生孔隙很小，次生孔隙的发育使岩石的刚性大大减弱，并呈现出明显的非均质性与各向异性，同时不同的构造部位受构造应力作用的强度难以确定，最小水平主应力和岩体抗张强度的度量较难，造成用测井资料计算的地层破裂压力精度较低。碳酸盐岩地层破裂压力与测井响应具有密切的关系。利用能够反映碳酸盐岩地层基本特性和岩石力学性质的测井信息，预测碳酸盐岩地层的破裂压力是一种经济、简便的可靠途径。 1957年，Hubbert和Willis根据三轴压缩试验，首先提出了地层破裂压力预测模式即H-W模式。到目前为止，国内外提出了许多预测地层破裂压力的方法。比较常用的有Eaton法，Stephen法，黄荣樽法等。1997年Holbrook发表了适于预测张性盆地裂缝扩展压力的一种方法。现场应用表明，修正后的模型具有较高的精度。 以上方法需要确定地层的泊松比、地层的构造应力系数、抗拉强度、室内岩心三轴试验和现场典型的破裂压力试验。 1.3地层破裂力学模型 压裂作业时，地层破裂力学模型如图1.1所示。此时，地层裂隙受地应力与压裂液共同作用。考虑深层水力压裂主要是形成垂直裂缝，且裂缝穿透整个油层。地应力与压裂液应力的最终有效合应力在裂隙壁面上是拉应力，当其合成应力强度因子K达到临界值时，裂隙就开始失稳延伸。 图1.1 压裂施工地层破裂模型 1.4破裂压力的应用 破裂压力数据在油田上应用较为广泛，应用于钻井、修井、压裂、试油、井下测试等井下工艺技术，钻井大多数是在裸眼中进行的，所以破裂压力数据在钻井方面尤为重要，它是钻井之前的井身结构设计，套管强度计算、钻井液密度设计等钻井工程设计内容的关键参数，特别是在一个新的区块开发之前，破裂压力这一数据为就重中之重了。它决定着在这一新的区域内的所有钻井是否正确， 并能否顺利执行和能否顺利完成。 2 国外破裂压力计算模型总结 地层的破裂压力对钻井液密度确定、井身结构和压裂设计施工等有着重要的指导作用。从上世纪五六十年代，国外就开始对地层破裂压力进行了研究，并取得了一系列的成果。主要如下: 1967年马修斯和凯利(Matthew，Kelly)提出了一个预测模型 ,,,，,,,,,) (2.1) ,,,, 式中 P— 地层破裂压力; f P— 地层空隙压力: p S — 上覆岩层压力: K— 随井深而变化的应力系数。 i 由于马修斯和凯利认为上覆层压力梯度等于1.0磅/平方英尺•英尺，是不随深度变化的常数，因而不符合实际情况。而且K值需要实际压裂资料来确定，所i 以未得到推广应用。 1969年伊顿(Eaton)提出上覆岩层压力梯度不是常数而是深度的函数，可由密度测井曲线求得，并把(1)式中的K值具体化为μ/(1-μ)，μ为地层的泊松i 比。提出预测破裂压力模式为 ,,,,，,,,,) (2.2) ,,,,,, 伊顿认为(2)式中的μ值应由地层破裂试验数据求得，提出西德克萨斯储积砂岩层的泊松比是不随深度变化的常数，其值等于0.25。但又认为墨西哥湾沿海地区的砂岩泊松比是随深度变化的，其值大于0.25，并在约2000m是都达到或超过0.4。由于伊顿的所谓泊松比是按(2)式反算的，其中包括了伊顿模式中未加考虑的经验周围的应力集中，地质构造应力和岩层的强度特性等在内的许多因素的影响，所以反算而得到的μ值不是岩层本身的真实泊松比，其数值明显偏大，有时超过不可压缩泊松比的上限值0.5，达到0.8以上。 1973年安德森(Anderson)等探索从测井资料中获得足以确定地层破裂压力的系数，考虑了井壁上应力集中的影响，并根据特查希(Terzadhi)的试验结果对比奥特(Biot)弹性多孔介质的应力、应变关系式进行简化后到处了预测地层破裂压力的模式为: ,,,,,，,,,,) (2.3) ,,,,,, 安德森提出用测井资料确定砂岩泥质含量和孔隙度并找出它们与岩层泊松比的关系后才能确定(3)式中的μ值，而对非砂岩地层的破裂压力仍无法预测。由于导出(3)式时没有计入地下构造应力的影响，所以这个预测模式亦不具普遍意义。 1982年斯蒂芬(Stephen)提出了再预测破裂压力的模式中考虑构造应力的问题，但又做了均匀水平构造应力的假设，其预测模式为: ,,,,，(，,),,,,) (2.4) ,,,,,, 式中 ξ—均匀构造应力系数可由实测破裂压力推算。 可见，斯蒂芬公式只是伊顿公式的改进，多了一项均匀构造应力系数，但是在水平方向均匀构造应力的假设是不符合全世界多数地区的地应力状况的。斯蒂芬主张用在常压下测得的动弹模量推算的泊松比值而没有考虑地下岩层围压的作用以及动弹模量和静弹模量之间的差别所应进行的修正。 上述四个模式中，均采用了岩层抗张强度为零的假说，这也与实际情况不符，而岩层抗张强度对其破裂压力也是有明显影响的。 1997年Holbrook 发表了适于预测张性盆地裂缝扩展压力的一种方法: ,,,,,,,,,,,,，, (2.5) ,,,, (1-Ф)表示地层的压实程度，经现场验证该方法对于泥岩地层适用性较好,但对于砂岩地层预测值偏高。 1999年M.M.Hossain提出了新的破裂压力理论即裸眼斜井的破裂压力理论。 3 国内破裂压力计算模型总结 1986年黄蓉樽考虑到一般地应力是不均匀的，在三向应力的影响下，考虑井眼周围处于平面应力状态，利用弹性理论中kursh关于无限平板中的小圆孔周围应力的解，推导出了地层破裂压力公式: ,,,,,,，,,,,,，,，, (3.1) ,,,, ,,,, 式中参数如下: ,——地层泊松比; ,——地层破裂压力; , ,——地层上覆岩层压力; , ,——地层孔隙压力; , ,,,,,, 非均质地质构造应力系数; ,，, ——水平两个主应力方向构造应力系数; ,——地层抗拉强度。 , 在不考虑岩石抗张强度的基础上，公式可归结为 ,,,，,,,,,) (3.2) ,,,,, 该公式没有考虑岩石中孔隙压力的作用和滤液侵入岩石的影响，而且假设了岩石满足特查希的有效应力条件(即作用于岩石固体骨架上的有效应力等于正应力减去孔隙压力)。 2000年李传亮、孔祥言提出了裸眼井完井条件下破裂压力的计算公式: 图1 裸眼完井垂直裂缝示意图 在裸眼完井条件下对油井进行压裂,垂直裂缝将沿着最大水平主应力的平行方向延伸,当裂缝开始形成时,井底流压即地层岩石的破裂压力计算公式为: ,,,,，,,,,,,,,,,, (3.3) ,,，,,,,; 裸眼完井条件下地层产生垂直裂缝时的岩石破裂压力除了与岩石的性质参数Φ和η有关外,主要受水平地应力参数σh 和σH 的影响。按照上式的推导方法,推导出裸眼完井条件下地层产生水平裂缝时的岩石破裂压力计算公式: ,，,,,,,,,,,,, (3.4) ,,,,,; 2001年葛洪魁、林英松、马善洲等人提出了修正Holbrook地层破裂压力预测模型。指出虽然Holb roo k的破裂压力与孔隙度的关系较好地解释了井漏容易在疏松砂岩地层中发生的现象，但对高孔隙度地层不适用。基于“临界孔隙度”概念和岩石力学特性通用预测模型, 对Holbrook地层破裂压力预测模型进行了修正得出如下公式: ,,,,,,,,,,？,,，, (3.5) ,,,,,; 2002年邓金根、王金凤、周建良根据线性孔隙弹性理论，在考虑孔隙压力及库伦-摩尔准则有效应力的情况下，推导出了地层渗透和地层不渗透两种情况下的破裂压力计算公式: ,,,,,,,,,,,,,？,,, ,，,,，, (3.6),,,,,,,,,, ,,,,,,,,,,,,,,,,,,,，，？,,,，,,，, (3.7) ,,,,,,,,,,，,,,,,,, 式中:,为地层的静态泊松比; , ,为地层的静态杨氏模量(MPa); , ,， ,为构成应力系数。 ,, 2002年李传亮提出了射孔完井条件下破裂压力的计算方法: 图2 油井射孔孔眼示意图 图3 射孔孔眼垂直裂缝示意图 对于射孔完井,情况则完全不同。由于油层段下了套管,地层是通过射孔孔眼与井筒进行联系的。高压液体首先从井筒流入射孔孔眼,然后通过孔眼把地层岩石压开。每个孔眼就相当于裸眼完井条件下的1个井眼。在所有的孔眼中,与最小水平主应力垂直或与最大水平主应力平行的孔眼中最容易产生垂直裂缝,在孔眼中产生垂直裂缝时的岩石破裂压力计算公式为: ,,，？,,,,,？,,,,,,, ,,,,, ,,，,,,,, 此公式的推导方法与裸眼井垂直裂缝破裂压力公式的推导方法类似。由射孔井的破裂压力的公式可见:射孔完井条件下地层产生垂直裂缝时的岩石破裂压力除了与岩石的性质参数Φc 和η有关外, 还要受地应力条件参数σh 和Pob的影响,而与地层的最大水平主应力σH 没有关系。 射孔完井条件下地层产生水平裂缝时pb 为: ,？，？,,,,,,,,,,,,, ,,,,, ,,，,,,,, 由上式可以看出:射孔完井条件下地层产生水平裂缝时的岩石破裂压力除了 与岩石的性质参数Φc和η有关外，主要受地应力条件参数σH 和pob的影响，而与地层的最小水平主应力σh 没有关系。所以，射孔完井条件下产生水平裂缝时的岩石破裂压力计算公式与裸眼完井条件下产生水平裂缝时的岩石破裂压力计算公式完全不同。 2003年胡永全、赵金洲、曾庆坤等提出了计算射孔井水力压裂破裂压力的有限元方法(有限元方法是一种以能量原理和变分法为理论基础,以矩阵代数和计算方法为主要工具的近似数值分析方法。广泛用于结构分析,以确定研究对象的应力与应变变化规律。)。他们将套管和岩石视为具有不同性能的两种材料,应用引进的CAD/ CA E/ CAM Pro/ E有限元软件在SGI工作站上计算得到了射孔孔眼周围的应力分布。根据主应力方位确定了压裂裂缝方位，结合岩石破裂准则得到水力压裂施工的破裂压力为。 2004年边芳霞、林平、王力等考虑到地层岩石存在天然微裂隙这一客观事实，建立了地应力与压裂液共同作用时，地层裂隙的应力强度因子表达式;并根据断裂判据，推得了地层破裂压力计算模型如下: 表1 不同完井方式破裂压力计算模型 式中个参数如下: K1 —套管承载系数(通常取K1 = 4165) ; K2 —射孔套管系数; K3 —无孔套管系数; 2Ps —孔隙压力,kg/ cm; Rw—井半径,cm ; 2St —岩石抗张强度,kg/ cm; μ—岩石泊松比; z—压裂层深度,m ; 2σ—有效应力值,kg/ cm; 3γ—岩石重度,kg/ cm; ζ—侧压系数; K渗—渗透系数。 2004年郭凯俊、常培锋从地应力与岩石力学参数研究入手，分析井周应力场分布，根据起裂准则，建立了浅层破裂压力预测模型。 直井破裂压力预测模型: (1)垂直缝的破裂压力: ,,,,,,,,,？,,,,,,,，,,？,,,,,,,,,, ,,,,,,,,,,,,,,,,,, 若考虑地层低渗透，强度各向异性不明显，则 ,,,,,？,？,,,，, ,,,,, (2)水平缝的破裂压力: ,,,,,,，,？,,,,,, ，,,,,,,,,,,,,,，,,,, 若考虑地层低渗透，强度各向异性不明显，则 ,,,,？，, ,,, 破裂压力P为: f ,, ,,,,,,,,，, ,,, 式中:p为孔隙压力，v为泊松比，Ф为孔隙度。 P 2007年张旭东认为孔隙性地层的渗透性具有较强的应力敏感性，井筒流体的流体力学特性及泵入排量影响井周地层渗流压力和地层应力分布，从而决定着井筒流体渗入地层的流量和沿程摩阻。考虑应力敏感性渗流对井口的破裂压力的影响，建立了孔隙性地层破裂压力的解释模型，得出: 在水力压裂作业时，地表压力表记录地层破裂时的压力为P，地层实际破裂t压力P 则为: f ,,,，,,,，, ,,,, 式中:, 为井筒沿程摩阻。 ,, 2008年刘翔在国内外相关应力分析研究的基础上，运用解析方法研究了射孔后孔眼围岩的地应力分布，通过应力分布计算,结合抗拉破坏准则，对垂直射 孔井水力压裂破裂压力进行了分析，得到了孔眼围岩应力分布和垂直射孔井水力压裂破裂压力的定性认识。如下: 由弹性力学理论，最大拉伸应力为: ,,,,,，？，？,？？,,,,,,,？, ,,,, 结合岩石拉伸破坏准则，引入有效应力，得到破裂压力的计算公式: ？,,,,, ,,,,, T 岩石的抗拉强度,MPa; σ孔眼围岩周向应力,MPa; Ф σzz 孔眼围岩轴向应力,MPa; ηp 孔隙压力贡献系数,无量纲; p 原始孔隙压力,MPa。 p 2009年王河、樊洪海考虑到目前预测地层破裂压力的方法中大多数的参数需经过大量的岩心试验获取实际工作中我们常常遇到没有岩心实验数据或者数据不足的情况，提出的是一种利用现场实测数据的方法，通过实测数据反算计算公式中的参数，然后计算破裂压力。在资料不全的情况下用此方法简单实用，现场应用亦很方便。主要过程如下: ,,,,,,, ,,,,, K 即水平有效应力与垂直有效应力的比值,称为侧压力系数。上覆岩层压力可由密度测井资料获得，待定参数K值的求取，我们可以从以下方法得到: (1) 查找地层破裂压力实验和地层破裂压力测试的数据，从中可以计算出地层的破裂压力，带人上式就可反算出该深度的侧压力系数; (2 )查找同区域地应力的数据，或者同区域其他井已取岩心的实验数据。构造平缓的区域，同一海拔地应力值差不多，对于构造剧烈的区域情况就不一样了。柏松比的计算可以，可以算出该深度的侧压力系数。 通过以上方法可以得到一组侧压力系数K 值。中间没有侧压力系数的地层,其K 值按线性插值方法求取，结合已经计算出来的上覆岩石压力和孔隙压力就可以计算出破裂压力。 2009年任岚、赵金洲、胡永全等人以多孔介质流体渗流与岩石应力–应变耦合理论为基础，推导水力压裂时岩石破裂压力数值计算方程，提出了一种全新的岩石破裂压力计算方法。该方法考虑岩石变形与流体渗流的全耦合作用，采用有限元法数值计算技术，能模拟计算流体向地层渗滤情况下井眼周围地层有效应力的时空分布。克服了传统破裂压力解析计算方法的许多不足，如无法精确计算井眼周围孔隙压力的升高导致应力集中加剧对破裂压力的影响，无法计算地层破裂的确切时间等问题;实现了岩石破裂压力数值计算的突破和计算精度的提高， 为水力压裂时岩石破裂压力的计算找到了新的理论和方法。其计算流程图如下: 图 计算流程图 2010年刘欣、刘从箐、刘同斌认为斜井破裂压力是压裂设计的一个重要参数，一般只能从已知射孔方位角来预测该参数。但现场采用的是非定向射孔,缺乏相应破裂压力的预测方法,而传统的破裂压力计算方法并没有获得非定向射井破裂压力的变化规律。因此，他们研究了相位对斜井射孔破裂压力的影响。 4结论 通过调研近年来国内外对于破裂压力的计算方法，这些方法都有各自的适用性。都没有全面考虑影响破裂压力的各种因素，因此，都有一定的实用范围。现对其优缺点及实用性进行简要总结。 (1)H-W方程是最早研究水力压裂机理的成果，它实际是计算的岩石破裂压力的上限值，而H-F公式计算的是岩石破裂压力的下限值，岩石的实际破裂压力介于他们之间，其数值随岩石性质变化而有所不同。 (2)伊顿Eaton首次提出了上覆岩层压力梯度不是常数而是深度的函数，其值可由密度测井曲线求得。但伊顿公式没有考虑井壁应力集中的影响和地层破裂的真实起因，更没有考虑到不均匀的原地应力或地质构造力对地层压力的作用。因次伊顿法往往局限于墨西哥湾地区。 (3)安德森法人为可以从电测资料中掌握砂岩泥质含量，来间接确定它的 μ值。但是这个关系式式随砂岩的结构特点而变的，且除砂岩以外的其他地层的破裂压力仍然无法预测。而且也没有充分考虑地下构造应力的影响，所以该预测破裂压力的模式也不具有普遍意义。 (4)黄荣樽法计算公式从地下岩层中两个水平主方向上作用着不等的地质构造应力的条件出发，分析井壁上应力集中现象，同时考虑岩层本身强度性质，提出了一个比较全面第考虑了各种因素的计算方法，在我国各主要油田中有较大的实用性。 (5)刘翔在研究了垂直射孔井破裂压力，得出了破裂压力随孔眼方位角增大而增加的结论，对射孔井的破裂压力计算有一定指导意义。 (6)任岚在提出的采用有限元数值计算技术进行模拟计算的方法，推导了流体渗流和岩石应力-应变的全耦合破裂压力方程，克服了传统破裂压力计算的许多不足，如无法计算井眼周围孔隙压力的升高导致应力集中加剧对破裂压力的影响，无法计算地层破裂的确切时间等问题。 (7)樊洪海提出了利用现场实测数据反算侧压力系数，然后计算破裂压力的方法是针对资料不全的情况下的一种经验与理论相结合的方法，对于没有水平构造力或水平构造力较小的地层有一定适用性。 4.1存在的问题及改进的方向 压裂过程中，地层岩石的破裂压力与油井的完井方式有很大关系。射孔完井条件下地层岩石破裂压力的计算公式与裸眼完井条件下的公式有所不同。油井压裂设计过程中，应根据油井的完井条件选择相应的计算公式来确定岩石的破裂压力。 目前，射孔完井的应用日益广泛，水力压裂的对象绝大多数都是射孔完井方式。因此，射孔参数对于破裂压力有重要影响，这方面的研究还有待深化。 对于低渗油井，特殊类油气藏等，破裂压力的计算方法研究不多。因此，从以下方面加强破裂压力的研究。 (1)寻求新的、方便快捷的计算方法来求解破裂压力，增强实用性。 (2)进一步提高破裂压力计算公式中参数的精度，进而提高破裂压力数值的准确性。 (3)进一步加强对低渗油气藏，页岩气等非常规油气藏的破裂压力计算模型的研究。 参考文献 [1]刘翔.垂直射孔井地应力及破裂压力研究.钻采工艺， 2008年31(2):36-38 [2]黄荣樽.地层破裂压力预测模式的探讨.华东石油学院学报，1984(4):336-347 [3]王河，樊洪海.一种地层破裂压力的估算方法.新疆石油科技，2009，19(4): 11-13 [4]李传亮，孔祥言. 油气井压裂过程中岩石破裂压力计算公式的理论研究.石油钻采工艺，2000，22(2):54-56 [5]任岚，赵金洲，胡永全等.水力压裂时岩石破裂压力数值计算.岩石力学与工程学报.2009,9月，28(2):3418-3422 [6]李传亮.射孔完井条件下的岩石破裂压力计算公式.石油钻采工艺.2003，24(2):37-38 [7]胡永全，赵金洲，曾庆坤等.计算射孔井水力压裂破裂压力的有限元方法.天然气工业.2003,23(2):58-59