CN101278209B - 钻前孔隙压力预测方法和系统 - Google Patents
钻前孔隙压力预测方法和系统 Download PDFInfo
- Publication number
- CN101278209B CN101278209B CN2006800366219A CN200680036621A CN101278209B CN 101278209 B CN101278209 B CN 101278209B CN 2006800366219 A CN2006800366219 A CN 2006800366219A CN 200680036621 A CN200680036621 A CN 200680036621A CN 101278209 B CN101278209 B CN 101278209B
- Authority
- CN
- China
- Prior art keywords
- sigma
- stress
- relevant
- pore pressure
- delta
- 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.)
- Expired - Fee Related
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 173
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000035945 sensitivity Effects 0.000 claims abstract description 87
- 238000005553 drilling Methods 0.000 claims abstract description 22
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 35
- 239000000463 material Substances 0.000 claims description 26
- 230000006835 compression Effects 0.000 claims description 24
- 238000007906 compression Methods 0.000 claims description 24
- 239000012530 fluid Substances 0.000 claims description 15
- 238000010008 shearing Methods 0.000 claims description 9
- 230000000704 physical effect Effects 0.000 claims description 6
- 150000003839 salts Chemical class 0.000 claims description 6
- 230000002706 hydrostatic effect Effects 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims 4
- 239000011435 rock Substances 0.000 description 18
- 230000008859 change Effects 0.000 description 12
- 230000015572 biosynthetic process Effects 0.000 description 9
- 238000005755 formation reaction Methods 0.000 description 9
- 239000004576 sand Substances 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000013076 target substance Substances 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000010287 polarization Effects 0.000 description 4
- 241000425571 Trepanes Species 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000003252 repetitive effect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910001369 Brass Inorganic materials 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 235000015076 Shorea robusta Nutrition 0.000 description 1
- 244000166071 Shorea robusta Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000010951 brass Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/284—Application of the shear wave component and/or several components of the seismic signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Measuring Fluid Pressure (AREA)
- Earth Drilling (AREA)
Abstract
一种孔隙压力预测方法,包括获得应力灵敏度系数,获得针对钻前位置的压缩波(P-波)速度和剪切波(S-波)速度,以及获得第一预测孔隙压力。另外,所述方法包括重复执行:使用与钻前位置有关的第一预测孔隙压力,计算与钻前位置有关的总应力值;以及使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从针对钻前位置的P-波速度和S-波速度组成的组中选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,并工具第二预测孔隙压力,调整与钻前位置有关的钻井操作。
Description
相关申请的交叉参考
本申请要求Colin Michael Sayers和Richard Anthony Birchwood于2005年8月8日申请的题为“Sub-Salt Pore Pressure Prediction”的美国临时专利申请No.60/595,814的优先权,本文将该申请全文引为参考文献。
背景技术
有些物质(如盐等)具有低渗透性,即不容易透过流体。在有些情况下,低渗透性物质可能处于钻井操作目标物质的上面。
图1示出钻井操作的示意图。本领域技术人员可以理解,图1中所示之钻井操作仅用于例举之目的,因而不应当被理解为限制本发明的范围。例如图1中所示的钻井操作是海底钻井操作,而钻井操作也可以是陆地钻井操作。
如图1中所示,将钻井机(105)设计成使用与钻杆柱(125)末端连接的钻头(未示出)钻入岩层(例如海底(110)下面的岩层)中。具体地说,使用钻头钻出延伸至目标物质(120)的钻孔(130)。目标物质(120)可以是钻井操作目的所在的碳氢化合物或矿物资源。另外,如图1中所示,低渗透性物质(115)可能处于目标物质(120)的上面。因而,为了到达目标物质(120),要求钻孔(130)穿过低渗透性物质(115)。
当低渗透性物质的沉积物被埋藏或者压紧时,有可能将流体挤压到所生成的结构内(即低渗透性物质本身内,和/或低渗透性物质下面的物质(如沙子等)内)的孔隙中。按照这种方式被挤压的流体向周围岩层施加被称做是孔隙压力的压力。在给定深度处孔隙压力超过静水压力的岩层,被称做过压。
当钻入过压岩层时,泥浆重量(即输送到钻孔的钻井液的重量)必须足够大,以防止孔隙压力将岩层流体移动到钻孔中。在比较糟糕的情况下, 进入钻孔中的岩层流体可能造成钻井失败,和/或伤及操作钻井机的工作人员。因此,为了进行安全和经济的钻井,必须足够精确地预测孔隙压力。特别是,这样有利于在钻前预测孔隙压力,即在开始任何钻井操作之前和/或在钻头还没有到达的任何位置处预测孔隙压力。
一般地说,钻前孔隙压力预测是基于使用钻前地震速度,以及使用偏移井数据(即来自钻井位置附近其他井的数据)校准的速度-孔隙压力转换。不过,在有些情况下(比如钻到盐下面时),传统的钻前孔隙压力预测不够精确。在Sayers CM,Johnson GM和Denyer G.,2002,“Pre-drillPore Pressure Prediction Using Seismic Data”Geophysics,67,1286-1292中可以找到对传统钻前孔隙压力预测技术的进一步讨论。
在钻井操作中使用泥浆来冷却钻头,将钻井操作产生的切削屑输送到地面,防止岩层流体流入钻孔,并使钻孔稳定。为了防止岩层流体流入,钻井操作员必须保持泥浆重量等于或者高于孔隙压力。为了使钻孔稳定,钻井操作员调节泥浆重量(即所用泥浆的密度)),抵制钻孔塌陷的趋势。然而,钻孔操作员必须很小心,以至于不会由于使用过高泥浆重量而使岩层断裂。
另外,太高的泥浆重量必然导致不能接受的低钻井速度。因而,泥浆重量必须较低,以至足以保持可接受的钻井速度,并且避免使岩层断裂。在这种情况下,当钻入过压岩层中,允许的泥浆重量窗口(即容许泥浆重量的范围)可能比较小。特别是,泥浆施加的力必须处于孔隙压力(或者防止塌陷的压力,如果高于孔隙压力的话)和使岩层断裂需要的压力之间的范围以内。
另外,当钻入过压岩层时,可能会增大所需套管柱(即插入钻孔内的结构支撑体)的数量。特别是,如果没有足够精确的钻前孔隙压力预测,为了防止钻井控制问题(如岩层流体的流入)和/或钻孔失败,有可能过早地插入附加的套管柱。过早插入套管柱会延缓钻井操作和/或减小钻孔的尺寸,造成经济损失。
发明内容
概括而言,按照一种方案,本发明涉及一种进行孔隙压力预测的方法。 所述方法包括获得应力灵敏度系数,获得钻前位置的压缩波(P-波)速度和剪切波(S-波)速度,获得第一预测孔隙压力,重复执行:使用与钻前位置有关的第一预测孔隙压力计算与钻前位置有关的总应力值;并使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从钻前位置的P-波速度和S-波速度组成的组选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,其中,如果第一预测孔隙压力与第二预测孔隙压力之差小于或等于阈值差值,则完成重复过程;而若所述差值大于阈值差值,则使用第二预测孔隙压力来更新第一预测孔隙压力,并根据第二预测孔隙压力调节与钻前位置有关的钻井操作。
概括而言,按照一种方案,本发明涉及一种计算机可读介质,所述介质包含执行进行孔隙压力预测的方法的指令。所述方法包括:获得应力灵敏度系数,获得钻前位置的压缩波(P-波)速度和剪切波(S-波)速度,获得第一预测孔隙压力,重复执行:使用与钻前位置有关的第一预测孔隙压力计算与钻前位置有关的总应力值;并使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从钻前位置的P-波速度和S-波速度组成的组选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,其中,如果第一预测孔隙压力与第二预测孔隙压力之差小于或等于阈值差值,则完成重复过程;而若所述差值大于阈值差值,则使用第二预测孔隙压力来更新第一预测孔隙压力,并根据第二预测孔隙压力调节与钻前位置有关的钻井操作。
概括而言,按照一种方案,本发明涉及一种用于执行孔隙压力预测方法的井下钻进工具。所述方法包括:获得应力灵敏度系数,获得钻前位置的压缩波(P-波)速度和剪切波(S-波)速度,获得第一预测孔隙压力,重复执行:使用与钻前位置有关的第一预测孔隙压力计算与钻前位置有关的总应力值;并使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从钻前位置的P-波速度和S-波速度组成的组选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,其中,如果第一预测孔隙压力与第二预测孔隙压力之差小于或等于阈值差值,则完成重复过程;而若所述差值大于阈值差值,则使用第二预测孔隙压力来更新第一预测孔隙压力,并根据第二预测孔隙压力调节与钻前位置有关的钻井操 作。
通过下面的描述和所附权利要求,将使本发明的其他方面愈为显而易见的。
附图说明
图1表示钻井操作的示意图;
图2表示本发明一种或多种实施例系统的视图;
图3表示本发明一种或多种实施例的流程图;
图4表示本发明一种或多种实施例的计算机系统的视图。
具体实施方式
以下参照附图详细描述本发明的具体实施例。各附图中为保持一致性,使用相同的附图标记表示相同元件。
在下面对本发明实施例的详细描述中,为了提供对本发明更充分的理解给出多个具体细节。然而,本领域普通技术人员显然应能理解,可以在不具有这些具体细节的情况下实施本发明。在其他情况下,为了避免使描述不必要地变得复杂,并未详细描述公知特征。
概括地说,本发明的实施例提供了钻前孔隙压力预测方法和系统。具体而言,使用参考数据、与钻前位置有关的声波或地震波速度、至少一种应力-速度关系公式,以及应力-压力关系公式预测钻前孔隙压力。按照本发明的一种或多种实施例,反复地预测孔隙压力,直至满足阈值条件为止。基于预测出的孔隙压力调节与钻前位置有关的钻井操作。
图2表示本发明一种或多种实施例系统的视图。具体地说,图2表示本发明一种或多种实施例计算环境(205)的视图。
按照本发明的一种或多种实施例,所述计算环境(205)包括一个或多个计算机系统(例如计算机系统A(210)、计算机系统N(215)),用以执行与钻井有关的任务。在本发明的一种或多种实施例中,所述计算机系统(如210,215)可以是网络服务器、嵌入式系统(比如位于井下钻进工具中的计算机)、桌面计算机、膝上型计算机、个人数字助理,以及任何其他同类计算机系统或者他们的任意组合。
具体而言,按照本发明的一种或多种实施例,所述一个或多个计算机系统(如210,215)包括应力灵敏度系数计算器(220)、总应力计算器(230)和/或孔隙压力计算器(235)。在本发明的一种或多种实施例中,可将上述部件(即220,230,235)设置于单个计算机系统(如210,215)中,分别设在多个计算机系统(如210,215)上,或者它们的任意组合。在本发明的一种或多种实施例中,上述每个部件(即220,230,235)可以包括一个或多个软件模块、一个或多个硬件模块,或者它们的组合。另外,按照本发明的一种或多种实施例,可以将上述部件(即220,230,235)配置成通过函数调用、应用程序接口(API)、网络协议(即有线或无线网络协议)、电子电路,以及任何其他同类通信和/或通信协议,或者它们的任意组合而彼此通信。
按照本发明的一种或多种实施例,可以将应力灵敏度系数计算器(220)配置成:利用参考数据,计算至少一种应力灵敏度系数。下面参照图3描述使用参考数据计算应力灵敏度系数。在本发明的一种或多种实施例中,可以将总应力计算器(230)配置成:利用预测孔隙压力,计算与钻前位置有关的总应力值。下面参照图3描述对总应力值的计算。
按照本发明的一种或多种实施例,可以将孔隙压力计算器(235)配置成:由使用总应力计算器计算出的总应力值计算预测孔隙压力。另外,按照本发明的一种或多种实施例,计算预测孔隙压力包括使用由应力灵敏度系数计算器(220)计算出的应力灵敏度系数。下面参照图3描述对预测孔隙压力的计算。
另外,按照本发明的一种或多种实施例,一旦孔隙压力计算器(235)已经计算出预测孔隙压力,则总应力计算器(230)可以使用预测孔隙压力重新计算总应力值。然后,孔隙压力计算器(235)使用重新计算出的总应力值,重新计算预测孔隙压力。按照本发明的一种或多种实施例,可以将孔隙压力计算器(235)和总应力计算器(230)配置成重复地重新计算总应力值和预测孔隙压力,直至预测孔隙压力满足阈值条件为止。下面参照图3描述对总应力值和预测孔隙压力的重复的重新计算。
图3表示本发明一种或多种实施例的流程图。具体而言,图3表示本发明一种或多种实施例的钻前孔隙压力预测方法流程图。
按照本发明的一种实施例,钻前位置相当于还没有进行钻井的位置;另一种说法是,钻头还没有到达钻前位置,不过,钻前位置处于钻头的预定路径中,除非改变钻孔的轨迹,钻头最终必将到达该钻前位置。在本发明的一个实施例中,可以在钻井中,钻井前(即开始钻出钻孔之前)或其任意组合中执行图3中所述的方法。
回到图3,可以使用参考数据计算一个或多个应力灵敏度系数(步骤305)。具体来说,按照本发明的一种或多种实施例,所述参考数据可以包括与一个或多个参考位置有关的值。按照本发明的一种或多种实施例,所述参考位置在一个或多个偏移井中可包括模拟地质岩层中的位置、任何其他同类的位置,或者它们的任意组合。
按照本发明的一种或多种实施例,所述参考数据包括—但不限于—总应力、有效应力、P-波速度、S-波速度,以及孔隙压力。可以利用声波岩芯记录测量(sonic log measurements),在实验室中进行的超声测量(比如参见Winkler,K.W.&Xingzhou,L.,1996,“Measurements ofthird-order Elastic Constants in Rocks”,J.Acoust.Soc.Am.,Vol.100,pp.1392-1398)、钻孔地震测量、垂直地震分布(VSP)、计算机模拟,任何其他方法或者它们的组合获得所述参考数据。
譬如,可以利用Modular Formation Dynamics Tester(MDT)工具或Repeat Formation Tester(RFT)工具获得偏移井中某一位置处的孔隙压力。作为可供选择的方式,可以使用比如Malinverno,A.等人的2004,“Integrating Diverse Measurements to Predict Pore Pressure withUncertainties While Drilling”SPE Paper 90001中所述的数值模块获得孔隙压力。
另外,可以使用应力模块获得参考位置处的总应力。下面的参考文献描述了本发明可使用的多种应力模块:Lee,D.,等人,2003,“Using aDynamic Mechanical Earth Model and Integrated Drilling Team toReduce Well Costs and Drilling Risks in San Martin Field”SPE PaperNo.84557和Moos,D.&Castillo,D.A.,2000,“Reservoir GeomechanicsApplied to Drilling and Completions Programs in challengingformations:North West shelf Timor Sea,North Sea and Columbia” APPEA Journal,February 16,pp.509-521。
按照本发明的一种或多种实施例,在步骤305中算出的应力灵敏度系数是,涉及与岩层中某一位置有关的压缩波(P-波)和剪切波(S-波)速度,以及该位置处的有效应力的系数。因而,在本发明的一种或多种实施例中,可以利用一个或多个应力-速度关系公式计算应力灵敏度系数。所述应力-速度关系公式可以是一阶、二阶、n阶(n为任意正整数),或者它们的任意组合。
在本发明的一个实施例中,所述应力-速度关系公式可包括下列公式中的一个或多个:
所述 作为一个系数,它使压缩波(P-波)速度的改变与σV(即有效垂直应力)的改变相关联。所述 作为一个系数,它使压缩波(P-波)速度的改变与σh(即有效水平应力)的改变相关联。所述 作为一个系数,它使剪切波(S-波)速度的改变与σV(即有效垂直应力)的改变相关联。所述 作为一个系数,它使剪切波(S-波)速度的改变与σh(即有效水平应力)的改变相关联。
上面所述的变量定义如下: 是与公共参考位置有关的有效垂直应力值, 是与公共参考位置有关的有效水平应力值,σV是与位置有关的有效垂直应力值,σh是与该位置有关的有效水平应力值,ΔσV是公共参考 位置处的有效垂直应力值与该位置处的有效垂直应力值之间的差值,Δσh是公共参考位置处的有效水平应力值与该位置处的有效水平应力值之间的差值, 是与公共参考位置有关的参考P-波速度, 是与公共参考位置有关的S-波速度,vP(σV,σh)是与位置有关的P-波速度,vS(σV,σh)是与位置有关的S-波速度。
利用公式(1)-(6),可以使用下述方法得到应力灵敏度系数(即 以及 )。最初,选择公共参考位置(L0)、第一位置(L1)和第二位置(L2)。本领域技术人员可以理解,这里使用的术语“第一”和“第二”仅用于将与一个位置有关的值同与另一位置有关的值区分开。因而,使用这些术语并不意味着顺序。继续该示例,L0,L1和L2可以位于偏移井内的不同深度处。在本发明的一个实施例中,L0位于偏移井中岩层的顶部。作为可供选择的方式,L0,L1和L2可以位于同一深度处,不过距离偏移井钻孔的轴不同距离处。譬如,L0,L1和L2可以位于海底下面10,000英尺处,其中L0位于钻孔壁上,L1位于沿垂直于钻孔轴的线测量时岩层内部3英寸的位置处,L2位于沿垂直于钻孔轴的线测量时岩层内6英寸的位置处。
在本发明的一个实施例中,所述L0,L1和L2处的岩层包含基本相同的物理性质(如岩性,多孔性等)。可以使用磁芯数据或岩芯记录(如电阻系数岩芯记录、中子多孔性岩芯记录,密度岩芯记录,核磁谐振岩芯记录等)检验物理性质的相似性。
一旦已经选择了L0,L1和L2,则得出针对每个位置(即L0,L1和L2)的vP(σV,σh)和vS(σV,σh)。按照本发明的一种或多种实施例,可以使用预测变井源距的垂直地震分布(VSP)得到P-波速度和S-波速度。在本发明的一个实施例中,所述VSP是多方位变井源距VSP。在Leaney,W.S.,Sayers,C.M.,和Miller,D.E.,1999,“Analysis of Multiazimuthal VSP Datafor Anisotropy and AVO”,Geophysics,64,1172-1180中可以寻找到对于多方位VSP的进一步讨论。
作为另一种可供选择的方式,可以使用Harrison,A.R.,1990,“Acquisition and Analysis of Sonic Waveforms from a BoreholeMonopole and Dipole Source and the Determination of Compressional and Shear Speads and their Relation to Rock Mechanical Propertiesand Surface Seismic Data”,SPE Paper No.20557,267-282中所描述的钻孔声波工具获得P-波速度和S-波速度。在美国专利No.5,838,633和6,904,365中描述了获得以距钻孔轴的距离为函数的P-波速度和S-波速度的其他方法。前面所述的所有专利都是本申请的受让人共同拥有的。
作为又一种可供选择的方式,按照本发明的一种或多种实施例,可以使用多分量地震勘测(如3C勘测、4C勘测等)获得P-波速度和S-波速度。譬如,可以使用P-波和S-波振荡器、P-波和S-波接收器,或者它们的任何组合,在地面上执行多分量地震勘测。作为另一个示例,可以使用位于海底上的P-波和S-波传感器,在海洋环境中进行多分量地震勘测。
作为再一种可供选择的方式,按照本发明的一种或多种实施例,可以使用PP和SS反射、PP和PS反射(比如在S-波接收器可变,但S-波源不可变的情况下)、任何其他同种类型的方法,或者它们的任何组合,获得P-波速度和S-波速度。另外,按照本发明的一种或多种实施例,可以通过转换与PP反射有关的振幅-偏移(AVO)数据,获得P-波速度和S-波速度。同样地,按照本发明的一种或多种实施例,可以使用与PS反射和SS反射有关的AVO数据。
回过来再来讨论使用公式(1)-(6)确定应力灵敏度系数,得出每个位置(即L0,L1,L2)处的孔隙压力。如上所述,可以使用多种向下钻进工具或多种方法获得孔隙压力。另外,还对每个位置(即L0,L1,L2)获得也被表示为张量的总应力。
利用每个位置处的孔隙压力和总应力,确定每个位置处的有效应力。在本发明的一个实施例中,将有效应力(σ)定义为总应力与表示成张量的孔隙压力之间的差值。然后,利用上面所述的有效应力(σ)与公式(5)和(6),确定L1和L2处的ΔσV和Δσh。利用针对L1和L2的vP(σV,σh),vS(σV,σh),ΔσV以及Δσh的值,创建四个公式(两个公式与公式(1)相应,两个公式与公式(2)相应)。在这个阶段,有四个公式和四个未知数(即 和 ),因而,可重新整理这些公式,求解,并获得 和 的值。
在本发明的一个实施例中,在最小水平应力和最大水平应力相等或者基本相等的情况下,可以使用公式(1)和(2)。不过,在最小水平应力与最 大水平应力不相等或不基本相等的情况下,可使用下面公式(1)和(2)的展开:
本领域技术人员可以理解,当SH与Sh不相等时,存在具有不同速度的具有不同偏振的剪切波。
按照本发明的一种或多种实施例,通过将微机械模型应用于参考数据,得到应力灵敏度系数。譬如,可以使用考虑了应力依赖性及颗粒(如砂岩中)之间的接触取向分布的微机械模型。在Sayers,C.M.,2002,“Stress-dependent elastic anisotropy of sandstones”,GeophysicalProspecting,50,85-95和Sayers,C.M.(1999)Stress-dependentseismic anisotropy of shales,Geophysics,64,93-98中可以寻找到对于微机械模型的进一步讨论。
有些情况下的垂直应力并非主要的应力。在这种情况下,可以使用将有效应力与P-波速度和S-波速度关联起来的公式(1)和(2)的更为通用的形式。例如,公式(1)和(2)可包括并非平行于或垂直于垂直轴(关于重力矢量)取向的剪切应力分量或主应力分量。
美国专利US 7,042,802;US 6,351,991;US 5,838,633;US 5,619,475;US 5,475,650;以及US 5,398,215中描述了用于获得应力灵敏度系数的其他方法。上述所有专利都为本申请的受让人共同拥有。
本领域技术人员可以理解,公式(1)和(2)中示例性的应力-压力关系依赖于以下的事实-垂直传播的P-波对于平行于传播方向作用的有效应力分量的改变最为敏感,S-波对于平行于传播方向作用的有效应力分量和平行于偏振方向作用的有效应力分量都敏感。另外,在Sayers,C.M.,2004,“Monitoring production-induced stress changes using seismicwaves”,第74届SEG年会增刊摘要中可以寻找到对于这些基本事实的进一步讨论。特别是由于对于P-波而言,传播方向与偏振方向基本平行,而对于S-波来说,传播方向与偏振方向垂直,P-波与S-波速度的组合使 得能够分别确定有效应力的垂直分量,从而能够确定孔隙压力以及主应力。
继续讨论图3,按照本发明的一种或多种实施例,在步骤310中,可以使用上面所述的一种或多种方法得出与钻前位置有关的P-波速度和/或S-波速度。在本发明的一个实施例中,可以使用通常被称做变井源距VSP的多偏移垂直地震分布(VSP),得到P-波速度和/或S-波速度。在这种情况下,可以使用沿一条线排列的用于最大化关于主应力信息提取的小球(shot)获得VSP。还可以使用多偏移、多方位VSP。
继续讨论图3,按照本发明的一种或多种实施例,在步骤315中,可以使用上面所述的一个或多个方法获得与参考位置有关的P-波速度和/S-波速度。
在步骤320中,可以使用预测孔隙压力计算与钻前位置有关的总应力值。特别是,按照本发明的一种或多种实施例,第一次执行步骤320时,可以从偏移井、计算机模拟或者任何其他来源获得预测孔隙压力。例如,可以使用偏移井数据上的速度-孔隙压力变换得到预测孔隙压力。在Eaton,“The Equation for Geopressure Prediction from Well Logs”,SPE 5544(Society of Petroleum Engineers of AIME,1975)和Bowers,“Pore Pressure Estimation from Velocity Data:Accounting forPore-Pressure Mechanisms besides Undercompaction”,SPE Drillingand Completion(1995年6月)89-95中可以寻找到对于速度-孔隙压力变换的进一步讨论。
按照本发明的一种或多种实施例,计算总应力值可以包括包含钻前位置和周围岩层的有限元分析。例如,有限元分析中使用的有限元模型可以包括钻前位置,以及一个或多个偏移井中的一个或多个位置。执行有限元分析还可以包括使用预测孔隙压力计算一个或多个深度上(如低渗透性物质上面、下面,和/或周围)的总应力改变。例如,可以由偏移井获得边界应力条件,并用于估计钻前位置处的总应力值。在Fredrich,J.T.等人,2003,“Stress Perturbations Adjacent to Salt Bodies in DeepwaterGulf of Mexico”,SPE Paper No.84554中可以寻找到对于有限元分析的进一步讨论。
作为可供选择的方式,可以使用有限差或任何其他类似的数值方法。下面的文章讨论了可以取代上述有限元方法使用的数值方法:Poliakov,A.N.B.等人,1993,“An Explicit Intertial Method for the Simulationof Viscoelastic Flow:An Evaluation of Elastic Effects on DiapiricFlow in Two-and Three-layer Models”,Stone,D.B.&Runcorn,S.K.(eds.)Flow and Creep in the Solar System:Observations,Modeling and Theory,pp.175-195。
继续讨论图3,按照本发明的一种或多种实施例,在步骤325,利用在步骤320中得到的总应力值计算第二预测孔隙压力。下面详细讨论了获得第二预测孔隙压力的实施例。
最初,选择一参考位置。步骤325中使用的参考位置可以是步骤305中使用的公共参考位置,或者可以为处于一个偏移井或者目标钻孔(即当前轨迹与钻前位置相交的钻孔)中的另一位置。在本发明的一个实施例中,所述参考位置和钻前位置包括基本相同的物理性质(如岩性,多孔性等)。可使用岩芯记录(如电阻系数岩芯记录,多孔性岩芯记录,密度岩芯记录,磁共振岩芯记录等)检验物理性质的相似性。
下式给出有效应力、总应力与孔隙压力之间的关系:
σij=Sij-αpδij (10)
式中的i和j表示张量的分量,σij是有效应力分量,Sij是总应力分量,α是多孔性弹性系数。如果i=j,则δij为1;如果i≠j,则δij为0。可以使用多种方法得到α,比如在Wang HF.,2000,“Theory of LinearPoroelasticity-with Applications to Geomechanics andHydrogeology”Princeton University,287pp中所描述的那样。本领域技术人员可以理解,上述参考文献中所述的方法并无意于限制本发明的范围。继续讨论步骤325,利用公式(8)取代公式(1)或(2)中的有效应力,得到:
对于这份参考文献而言,公式(11)-(12)中的变量定义如下: vP(Sv,Sh,p)是与钻前位置有关的P-波速度;vS(Sv,Sh,p)是与钻前位置有关的S-波速度;vP (0)是与参考位置有关的P-波速度;vS (0)是与参考位置有关的S-波速度;αP v是在步骤305中获得的垂直压缩应力灵敏度系数;αP h是在步骤305中获得的水平压缩应力灵敏度系数;αS v是在步骤305中获得的垂直剪切应力灵敏度系数;αS h是在步骤305中获得的水平剪切应力灵敏度系数;Sv是总垂直应力;Sh是总水平应力;p是孔隙压力;前面带有Δ的项表示钻前位置与参考位置处的值之间的差值;α为多孔弹性系数。
利用公式(11)和/或公式(12),以前述步骤中得到的数值确定Δp。然后,可将Δp与参考位置处的孔隙压力相加,得出第二预测孔隙压力。
继续讨论图3,按照本发明的一种或多种实施例,在步骤330,可以比较步骤320中使用的预测孔隙压力(下称“第一”预测孔隙压力)与步骤325中计算出的预测孔隙压力(下称“第二”预测孔隙压力)。具体而言,可以将第一和第二预测孔隙压力之间的一个或多个差值度量与一个或多个阈值差值进行比较。按照本发明的一种或多种实施例,阈值差值可以是差值度量的最大值,表示第二预测孔隙压力足够精确。在本发明的一个实施例中,当|(PP1-PP2)/PP1|≤0.001时,判断第二预测孔隙压力足够精确,其中PP1表示第一预测孔隙压力,PP2表示第二预测孔隙压力。本领域技术人员可以理解,阈值差值可使用其他数值,并且差值度量可使用其他数学结构。阈值差值可以为小于1的任何其他数值。
按照本发明的一种或多种实施例,如果在步骤330中第一与第二预测孔隙压力之间的差值大于阈值差值,则在步骤335中使用第二预测孔隙压力更新第一预测孔隙压力。在这种情况下,使用第一更新后的预测孔隙压力重复步骤320。按照本发明的一种或多种实施例,如果重复步骤320-325,则可以使用钻进时采集的当前数据(如地震波速度、应力数据等)。譬如,可以使用钻进的同时所进行的测量(MWD),钻进时的岩芯记录(LWD),任何其他类似种类的钻井数据,或者它们的任何组合采集当前数据。本领域技术人员可以理解,使用多个当前数据和/或接近钻前位置的某一位置的数据,可以提高所计算出的总应力值的精度。另外,本领域技术人员可以理解,“更新”第一预测孔隙压力可能简单地表明,使用在步骤325中计算出的预测孔隙压力重复步骤320。
按照本发明的一种或多种实施例,可以重复执行步骤320和325(即经由步骤330和335),直到第一与第二预测孔隙压力之间的差值小于或等于步骤330中的阈值差值为止。所述的另一种方法,可以重复步骤320和325,直至第一和第二预测孔隙压力充分收敛为止(比如第一与第二预测孔隙压力之差≤0.1%)。
按照本发明的一种或多种实施例,如果在步骤330中第一与第二预测孔隙压力之差小于或等于阈值差值,则在步骤340中可以根据第二预测孔隙压力调整与钻前位置有关的钻井操作。具体来说,按照本发明的一种或多种实施例,调整钻井操作可能涉及到调整钻井流体密度(即适当增大或减小钻井流体密度),调整钻井轨迹(比如避开过压区域,通过低压区域等),优化钻孔中套管柱的数量(即增加套管柱,延迟增加套管柱等),或者任何其他相似类型的调整。
在本发明的一个实施例中,如果所提议的钻井轨迹位于干扰周围应力场的复杂结构(比如盐)的附近,并且所提议的钻井轨迹穿过高渗透性岩层(比如沙),则可使用任何可用于确定此岩层中孔隙压力的技术来确定高渗透岩层中的孔隙压力。
例如,如果页岩(低渗透性物质)与沙子(高渗透性)接界,则可以在页岩中钻一口偏移井,并且可使用例如声波或地震速度来确定页岩中的孔隙压力。然后,可使用页岩中的孔隙压力,利用例如质心方法获得沙子中的孔隙压力。在序列号为WO2004/022905、被转让给本申请受让人的专利合作条约(PCT)申请中描述了一种具体的质心方法。质心方法的替代包括数字盆地建模。一旦确定了沙子中的孔隙压力,则可使用孔隙压力来验证多种应力灵敏度系数和/或用于产生钻前位置处的总应力的模型。
在本发明的一个实施例中,在钻井的同时使用电线和/或测井,得到偏移井中某一位置的孔隙压力测量值,其中,所述的位置处于沙子中,而且沙子延伸到钻前位置。在这种情况下,可以使用在偏移井中该位置处所得到的孔隙压力,以及偏移井中该位置与钻前位置之间的流体静力学梯度(由流体密度定义)的信息,估计钻前位置处的孔隙压力。一旦确定了钻前位置处沙子中的孔隙压力,就可以使用孔隙压力验证多种应力灵敏度系数和/或用于产生钻前位置处的总应力的模型。
在本发明的一个实施例中,可以使用在目标钻孔(即轨迹与钻前位置相交的钻孔)中钻进时获得的井下钻进测量,更新应力灵敏度系数和/或用于生成钻前位置处的总应力的模型。然后,在使用多偏移VSP获得的地震速度之前稍许时间,可使用这些系数预测孔隙压力。
按照本发明的一个或多个实施例,提供了一种用于精确预测钻前孔隙压力的方法。因而,采用本发明的一个或多个实施例,可以防止岩层流体进入钻孔,从而防止破坏井和/或伤害操作钻井机的工作人员。另外,按照本发明的一个或多个实施例,可以防止在财政上过早地插入套管柱。
实际上可以在任何种类的计算机上实现本发明,而与所使用的平台无关。例如,有如图4中所示那样,计算机系统(400)包括处理器(402)、相关的存储器(404)、存储装置(406),以及多个其他元件和当今计算机(未示出)的典型功能。计算机(400)还可以包括输入装置,如键盘(408)和鼠标(410),以及输出装置,如监视器(412)。计算机系统(400)可通过网络接口连接(未示出)与网络(414)(比如局域网(LAN),比如Internet的广域网(WAN),或者任何其他类似种类的网络)连接。本领域技术人员可以理解,这些输入和输出装置可采取任何其他形式。
另外,本领域技术人员可以理解,上述计算机系统(400)的一个或多个元件可以位于远程位置处,并通过网络与其他元件连接。另外,可以在具有多个节点的分布式系统上实现本发明,其中,本发明的每一部分(如应力灵敏度系数计算器、总应力计算器、孔隙压力计算器等)可以位于分布式系统内的不同节点处。在本发明的一个实施例中,节点与计算机系统相应。作为选择,所述节点可以相当于与物理存储器相关的处理器。所述节点还可以与具有共享存储器和/或资源的处理器相应。另外,可以将执行本发明实施例的软件指令存储到比如光盘(CD)、磁盘、磁带,文件的计算机可读介质上,或者任何其他计算机可读存储装置上。此外,在本发明的一个实施例中,可以通过图形用户界面(比如显示装置),向用户显示预测孔隙压力(包括使用图3中所示方法计算出的所有孔隙压力)。
尽管针对有限数量的实施例描述了本发明,但本领域技术人员得益于本说明书,可以想到不偏离本文所披露的本发明范围的其他实施例。因而,本发明的范围应当仅由所附权利要求限定。
Claims (38)
1.一种孔隙压力预测方法,包括如下步骤:
获得应力灵敏度系数;
获得钻前位置的压缩波速度和剪切波速度;
获得第一预测孔隙压力;
重复执行:
使用与钻前位置有关的第一预测孔隙压力,计算与钻前位置有关的总应力值;和
使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从钻前位置的压缩波速度和剪切波速度组成的组选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,
其中,若第一预测孔隙压力与第二预测孔隙压力之差小于或等于阈值差值,则完成重复,并且
若所述差值大于阈值差值,则使用第二预测孔隙压力更新第一预测孔隙压力;以及
根据第二预测孔隙压力,调整与钻前位置有关的钻井操作。
2.如权利要求1所述的方法,其中,使用应力-速度关系公式、偏移井中的第一位置、偏移井中的第二位置以及公共参考位置,执行获得应力灵敏度系数的操作。
3.如权利要求2所述的方法,其中,所述公共参考位置位于偏移井中。
4.如权利要求2所述的方法,其中,所述第一位置包括处于偏移井中的第一深度,所述第二位置包括处于偏移井中的第二深度。
5.如权利要求2所述的方法,其中,所述应力-速度关系公式包括选自由:
所组成的组中的至少一个,其中,
并且,为了计算应力灵敏度系数:
是与公共参考位置有关的有效垂直应力值,是与公共参考位置有关的第一有效最小水平应力值,σV是与从第一位置和第二位置组成的组选择出的一个位置有关的有效垂直应力值,σh是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最小水平应力值,ΔσV是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效垂直应力值之间的差值,Δσh是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最小水平应力值之间的差值,是与公共参考位置有关的参考压缩波速度,是与公共参考位置有关的参考剪切波速度,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是垂直剪切应力灵敏度系数,是最小水平剪切应力灵敏度系数,vP(σV,σh)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考压缩波速度,vS(σV,σh)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考剪切波速度,以及
6.如权利要求2所述的方法,其中,所述应力-速度关系公式包括选自由
组成的组中的至少一个,其中:
并且,为了计算应力灵敏度系数:
是与公共参考位置有关的有效垂直应力值,是与公共参考位置有关的第一有效最小水平应力值,是与公共参考位置有关的第一有效最大水平应力值,σV是与从第一位置和第二位置组成的组选择出的一个位置有关的有效垂直应力值,σh是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最小水平应力值,σH是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最大水平应力值,ΔσV是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效垂直应力值之间的差值,Δσh是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最小水平应力值之间的差值,ΔσH是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最大水平应力值之间的差值,是与公共参考位置有关的参考压缩波速度,是与公共参考位置有关的第一参考剪切波速度,是与公共参考位置有关的第二参考剪切波速度,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是最大水平压缩应力灵敏度系数,是第一垂直剪切应力灵敏度系数,是第二垂直剪切应力灵敏度系数,是第一最小水平剪切应力灵敏度系数,是第二最小水平剪切应力灵敏度系数,是第一最大水平剪切应力灵敏度系数,是第二最大水平剪切应力灵敏度系数,vP(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考压缩波速度,vS1(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的第一参考剪切波速度,vS2(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的第二参考剪切波速度,
其中,剪切波包括第一参考剪切波速度和第二参考剪切波速度,并且
7.如权利要求2所述的方法,其中,从由公共参考位置、偏移井中的第一位置和偏移井中的第二位置组成的组中选择的至少一个,包含与钻前位置基本相同的物理性质。
8.如权利要求2所述的方法,其中,计算总应力值包括使用数学模型,并且数学模型计算包括在一个体积中的总应力值,所述体积包含钻前位置以及从由公共参考位置、偏移井中的第一位置和偏移井中的第二位置组成的组中选择的至少一个。
9.如权利要求1所述的方法,其中,所述钻前位置位于低渗透性物质的下面。
10.如权利要求9所述的方法,其中,所述低渗透性物质为盐。
11.如权利要求1所述的方法,其中,所述计算应力灵敏度系数包括使用微机械模型。
12.如权利要求1所述的方法,其中,所述获得压缩波速度和剪切波速度包括执行提前变井源距垂直地震剖面(VSP)。
13.如权利要求1所述的方法,其中,通过使用偏移井中某一位置处的孔隙压力,偏移井中该位置与钻前位置之间的流体静力学梯度的信息,获得钻前位置处的替代孔隙压力,从而验证应力灵敏度系数;并且,偏移井中的该位置处于高渗透性物质中,而且所述高渗透性物质延伸到钻前位置。
14.如权利要求1所述的方法,其中,通过获得钻前位置处的替代孔隙压力验证应力灵敏度系数;使用偏移井中的孔隙压力和质心方法来确定替代孔隙压力;并且,所述偏移井中的位置处于低渗透性物质中,钻前位置处于高渗透性物质中,以及所述低渗透性物质与高渗透性物质接触。
15.如权利要求1所述的方法,其中,所述应力-压力关系公式为:
σij=Sij-αpδij
式中的i和j表示张量的分量,σij是有效应力分量,Sij是总应力分量,α是多孔性弹性系数,p与选自第一预测孔隙压力、第二预测孔隙压力和参考位置处的孔隙压力组成的组中的一个相对应;若i=j,则δij为1;若i≠j,则δij为0。
16.如权利要求1所述的方法,其中,使用选自由:
组成的组中的至少一个公式计算第二预测孔隙压力,式中的vP(Sv,Sh,p)是与钻前位置有关的压缩波速度,vS(Sv,Sh,p)是与钻前位置有关的剪切波速度,是与参考位置有关的参考压缩波速度,是与参考位置有关的参考剪切波速度,Sv是从总应力值得出的总垂直应力,Sh是从总应力值得出的总最小水平应力,p与选自第一预测孔隙压力、第二预测孔隙压力和参考位置处的孔隙压力组成的组中的一个相应,α为多孔性弹性系数,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是垂直剪切应力灵敏度系数,是最小水平剪切应力灵敏度系数,并且
前面带有Δ的项表示钻前位置与参考位置处的值之间的差值。
17.如权利要求1所述的方法,其中,调整钻井操作包括选自由调整钻井流体密度、调整钻井轨迹,以及优化钻孔中套管柱数量组成的组中的至少一个。
18.如权利要求1所述的方法,其中,所述钻前位置包括在钻孔中正在操作的钻头下面的位置。
19.如权利要求1所述的方法,其中,在图形用户界面上显示第一预测孔隙压力和第二预测孔隙压力。
20.一种井下钻进工具,它被构造成执行孔隙压力预测方法,所述方法包括:
获得应力灵敏度系数;
获得钻前位置的压缩波速度和剪切波速度;
获得第一预测孔隙压力;
重复执行:
使用与钻前位置有关的第一预测孔隙压力,计算与钻前位置有关的总应力值;和
使用应力-压力关系公式、应力-速度关系公式、应力灵敏度系数、参考位置,以及从钻前位置的压缩波速度和剪切波速度组成的组选择出的至少一个,计算与钻前位置有关的第二预测孔隙压力,
其中,若第一预测孔隙压力与第二预测孔隙压力之差小于或等于阈值差值,则完成重复,并且
若所述差值大于阈值差值,则使用第二预测孔隙压力更新第一预测孔隙压力;以及
根据第二预测孔隙压力,调整与钻前位置有关的钻井操作。
21.如权利要求20所述的井下钻进工具,其中,使用应力-速度关系公式、偏移井中的第一位置、偏移井中的第二位置以及公共参考位置,执行获得应力灵敏度系数的操作。
22.如权利要求21所述的井下钻进工具,其中,所述公共参考位置位于偏移井中。
23.如权利要求21所述的井下钻进工具,其中,所述第一位置包括处于偏移井中的第一深度,所述第二位置包括处于偏移井中的第二深度。
24.如权利要求21所述的井下钻进工具,其中,所述应力-速度关系公式包括选自由:
组成的组中的至少一个,
式中的
并且,为了计算应力灵敏度系数:
是与公共参考位置有关的有效垂直应力值,是与公共参考位置有关的第一有效最小水平应力值,σV是与从第一位置和第二位置组成的组选择出的一个位置有关的有效垂直应力值,σh是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最小水平应力值,ΔσV是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效垂直应力值之间的差值,Δσh是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最小水平应力值之间的差值,是与公共参考位置有关的参考压缩波速度,是与公共参考位置有关的参考剪切波速度,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是垂直剪切应力灵敏度系数,是最小水平剪切应力灵敏度系数,vP(σV,σh)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考压缩波速度,vS(σV,σh)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考剪切波速度,以及
25.如权利要求21所述的井下钻进工具,其中,所述应力-速度关系公式包括选自由:
组成的组中的至少一个,
式中的
并且,为了计算应力灵敏度系数:
σV是与从第一位置和第二位置组成的组选择出的一个位置有关的有效垂直应力值,σh是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最小水平应力值,σH是与从第一位置和第二位置组成的组选择出的一个位置有关的有效最大水平应力值,ΔσV是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效垂直应力值之间的差值,Δσh是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最小水平应力值之间的差值,ΔσH是从第一位置和第二位置组成的组选择出的一个位置与公共参考位置处的有效最大水平应力值之间的差值,是与公共参考位置有关的参考压缩波速度,是与公共参考位置有关的第一参考剪切波速度,是与公共参考位置有关的第二参考剪切波速度,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是最大水平压缩应力灵敏度系数,是第一垂直剪切应力灵敏度系数,是第二垂直剪切应力灵敏度系数,是第一最小水平剪切应力灵敏度系数,是第二最小水平剪切应力灵敏度系数,是第一最大水平剪切应力灵敏度系数,是第二最大水平剪切应力灵敏度系数,vP(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的参考压缩波速度,vS1(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的第一参考剪切波速度,vS2(σV,σh,σH)是与从第一位置和第二位置组成的组选择出的一个位置有关的第二参考剪切波速度,
所述剪切波包括第一参考剪切波速度和第二参考剪切波速度,并且
26.如权利要求21所述的井下钻进工具,其中,从由公共参考位置、偏移井中的第一位置和偏移井中的第二位置组成的组的中选择的至少一个,包括与所述钻前位置基本相同的物理性质。
27.如权利要求21所述的井下钻进工具,其中,计算总应力值包括使用数学模型;数学模型计算包括在一个体积中的总应力值,所述体积包含所述钻前位置以及从由公共参考位置、偏移井中的第一位置和偏移井中的第二位置组成的组中选择的至少一个。
28.如权利要求20所述的井下钻进工具,其中,所述钻前位置位于低渗透性物质的下面。
29.如权利要求28所述的井下钻进工具,其中,所述低渗透性物质为盐。
30.如权利要求20所述的井下钻进工具,其中,所述计算应力灵敏度系数包括使用微机械模型。
31.权利要求20所述的井下钻进工具,其中,所述获得压缩波速度和剪切波速度包括执行提前变井源距垂直地震剖面(VSP)。
32.如权利要求20所述的井下钻进,其中,通过使用偏移井中某一位置处的孔隙压力、有关偏移井中该位置与钻前位置之间的流体静力学梯度的信息,获得钻前位置处的替代孔隙压力,从而验证应力灵敏度系数;并且,所述偏移井中的该位置处于高渗透性物质中,而且所述高渗透性物质延伸到钻前位置。
33.如权利要求20所述的井下钻进工具,其中,通过获得钻前位置处的替代孔隙压力来验证应力灵敏度系数;使用偏移井中的孔隙压力和质心方法来确定替代孔隙压力;并且,所述偏移井中该位置处于低渗透性物质中,钻前位置处于高渗透性物质中,而且所述低渗透性物质与高渗透性物质接触。
34.如权利要求20所述的井下钻进工具,其中,所述应力-压力关系公式为:
σij=Sij-αpδij
式中的i和j表示张量的分量,σij是有效应力分量,Sij是总应力分量,α是多孔性弹性系数,p与选自第一预测孔隙压力、第二预测孔隙压力和参考位置处的孔隙压力组成的组中的一个相对应;若i=j,则δij为1;若i≠j,则δij为0。
35.如权利要求20所述的井下钻进工具,其中,使用选自由:
组成的组中的至少一个公式计算第二预测孔隙压力,式中的vP(Sv,Sh,p)是与钻前位置有关的压缩波速度,vS(Sv,Sh,p)是与钻前位置有关的剪切波速度,是与参考位置有关的参考压缩波速度,是与参考位置有关的参考剪切波速度,Sv是从总应力值得出的总垂直应力,Sh是从总应力值得出的总最小水平应力,p与选自第一预测孔隙压力、第二预测孔隙压力和参考位置处的孔隙压力组成的组中的一个相应,α为多孔性弹性系数,是垂直压缩应力灵敏度系数,是最小水平压缩应力灵敏度系数,是垂直剪切应力灵敏度系数,是最小水平剪切应力灵敏度系数,并且
前面带有Δ的项表示钻前位置与参考位置处的值之间的差值。
36.如权利要求20所述的井下钻进工具,其中,调整钻井操作包括选自由调整钻井流体密度、调整钻井轨迹,以及优化钻孔中套管柱的数量组成的组中的至少一个。
37.如权利要求20所述的井下钻进工具,其中,所述钻前位置包括在钻孔中正在操作的钻头下面的位置。
38.如权利要求20所述的井下钻进工具,其中,在图形用户界面上显示第一预测孔隙压力和第二预测孔隙压力。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US59581405P | 2005-08-08 | 2005-08-08 | |
US60/595,814 | 2005-08-08 | ||
US11/499,931 | 2006-08-07 | ||
US11/499,931 US7299132B2 (en) | 2005-08-08 | 2006-08-07 | Method and system for pre-drill pore pressure prediction |
PCT/US2006/030885 WO2007019497A2 (en) | 2005-08-08 | 2006-08-08 | Method and system for pre-drill pore pressure prediction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101278209A CN101278209A (zh) | 2008-10-01 |
CN101278209B true CN101278209B (zh) | 2011-11-16 |
Family
ID=37743593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006800366219A Expired - Fee Related CN101278209B (zh) | 2005-08-08 | 2006-08-08 | 钻前孔隙压力预测方法和系统 |
Country Status (9)
Country | Link |
---|---|
US (2) | US7299132B2 (zh) |
EP (1) | EP1913423A2 (zh) |
CN (1) | CN101278209B (zh) |
AU (1) | AU2006278297B2 (zh) |
CA (1) | CA2618555C (zh) |
EG (1) | EG24897A (zh) |
MX (1) | MX2008001875A (zh) |
NO (1) | NO20081227L (zh) |
WO (1) | WO2007019497A2 (zh) |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040034304A1 (en) * | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
US7486589B2 (en) * | 2006-02-09 | 2009-02-03 | Schlumberger Technology Corporation | Methods and apparatus for predicting the hydrocarbon production of a well location |
US8885440B2 (en) * | 2007-01-05 | 2014-11-11 | Madhumita Sengupta | Constructing velocity models near salt bodies |
US8098543B2 (en) * | 2007-01-05 | 2012-01-17 | Westerngeco L.L.C. | Estimation of stress and elastic parameters |
CN101685314B (zh) * | 2008-09-23 | 2014-07-09 | 成都芯源系统有限公司 | 自适应电压定位直流稳压器及其控制方法 |
US8522611B2 (en) * | 2009-02-19 | 2013-09-03 | Baker Hughes Incorporated | Method and apparatus for measuring pore pressure beyond the casing |
US8341984B2 (en) * | 2009-05-27 | 2013-01-01 | Westerngeco L.L.C. | Estimating velocities with uncertainty |
US9140816B2 (en) | 2009-08-18 | 2015-09-22 | Halliburton Energy Services, Inc. | Apparatus and method for determining formation anisotropy |
WO2011123774A2 (en) * | 2010-04-02 | 2011-10-06 | Geco Technology B.V. | Method and apparatus to build a three-dimensional mechanical earth model |
WO2012144922A1 (en) * | 2011-04-22 | 2012-10-26 | Baker Hughes Incorporated | Increasing the resolution of vsp ava analysis through using borehole gravity information |
CA2850482C (en) * | 2011-10-03 | 2018-10-30 | Landmark Graphics Corporation | Enhanced 1-d method for prediction of mud weight window for subsalt well sections |
US8725479B2 (en) | 2011-11-30 | 2014-05-13 | Baker Hughes Incorporated | Use of monte carlo simulations for establishing correlations and their uncertainties |
US9593567B2 (en) | 2011-12-01 | 2017-03-14 | National Oilwell Varco, L.P. | Automated drilling system |
AU2012376842B2 (en) * | 2012-04-09 | 2015-04-02 | Landmark Graphics Corporation | Compressional velocity correction apparatus, methods, and systems |
US9183182B2 (en) * | 2012-08-31 | 2015-11-10 | Chevron U.S.A. Inc. | System and method for determining a probability of well success using stochastic inversion |
US20140076632A1 (en) * | 2012-09-20 | 2014-03-20 | Baker Hughes Incoroporated | Method to predict overpressure uncertainty from normal compaction trendline uncertainty |
CN103089253B (zh) * | 2013-01-22 | 2015-06-03 | 中国石油大学(北京) | 一种利用小波变换计算地层孔隙压力的方法 |
EP3052970A1 (en) | 2013-12-19 | 2016-08-10 | Halliburton Energy Services, Inc. | Pore size classification in subterranean formations based on nuclear magnetic resonance (nmr) relaxation distributions |
CA2943346A1 (en) * | 2014-04-22 | 2015-10-29 | Landmark Graphics Corporation | Integrated 3d method for prediction of mud weight window for complex well sections |
CN104453879B (zh) * | 2014-11-14 | 2017-04-05 | 中国海洋石油总公司 | 钻前压力的预测方法 |
ES2640378T3 (es) | 2014-11-17 | 2017-11-02 | Repsol, S.A. | Método para gestionar la producción de un yacimiento petroquímico y producto de programa para el mismo |
US11022709B2 (en) | 2015-02-20 | 2021-06-01 | Board Of Regents, The University Of Texas System | Pore-pressure prediction based on velocities coupled with geomechanical modeling |
US11789170B2 (en) * | 2016-06-15 | 2023-10-17 | Schlumberger Technology Corporation | Induced seismicity |
CN106321090B (zh) * | 2016-08-25 | 2019-10-29 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 一种盐间地层孔隙压力的预测方法 |
US10754050B2 (en) | 2017-03-31 | 2020-08-25 | Chevron U.S.A. Inc. | Pore pressure prediction based on an integrated seismic and basin modeling approach |
US11035966B2 (en) * | 2017-10-17 | 2021-06-15 | Ha Dinh Nguyen | Method for estimating reservoir formation quality using seismic velocities |
CN108387931B (zh) * | 2018-01-22 | 2019-09-10 | 中国石油天然气股份有限公司 | 储层孔隙压力的确定方法和装置 |
CN108301825B (zh) * | 2018-01-23 | 2021-03-30 | 中国海洋石油集团有限公司 | 一种高温高压储层孔隙压力分析方法 |
CN109667573B (zh) * | 2018-12-12 | 2022-07-15 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 三维页岩储层孔隙压力预测方法、装置和电子设备 |
CN109931055B (zh) * | 2019-01-31 | 2020-09-25 | 西北大学 | 盆地深层复合成因的流体压力预测方法 |
CN110133725B (zh) * | 2019-05-08 | 2021-05-14 | 中国石油大学(华东) | 地震岩石横波速度预测方法及装置 |
CN111025386B (zh) * | 2019-12-13 | 2020-11-17 | 成都理工大学 | 一种无分离假象的纵横波分离方法 |
CN114076976B (zh) * | 2020-08-14 | 2024-03-26 | 中国石油天然气股份有限公司 | 一种利用纵横波速度联合预测有效应力系数的方法及系统 |
US11767751B2 (en) * | 2020-11-04 | 2023-09-26 | Reveal Energy Services, Inc. | Determining a characteristic associated with a reservoir |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5130949A (en) * | 1991-06-28 | 1992-07-14 | Atlantic Richfield Company | Geopressure analysis system |
US6694261B1 (en) * | 1999-06-07 | 2004-02-17 | Conoco Inc. | Method for identification of shallow water flow hazards using marine seismic data |
EP1410072A1 (en) * | 2000-10-10 | 2004-04-21 | ExxonMobil Upstream Research Company | Method for borehole measurement of formation properties |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5398215A (en) | 1993-11-19 | 1995-03-14 | Schlumberger Technology Corporation | Identification of stress induced anisotropy in formations |
US5475650A (en) | 1993-11-19 | 1995-12-12 | Schlumberger Technology Corporation | Measurement of nonlinear properties of formation using sonic borehole tool while changing pressure in borehole |
US5619475A (en) | 1994-03-30 | 1997-04-08 | Schlumberger Technology Corportion | Method of predicting mechanical failure in formation utilizing stress derivatives which measure formation nonlinearity |
US5838633A (en) | 1997-01-27 | 1998-11-17 | Schlumberger Technology Corporation | Method for estimating formation in-situ stress magnitudes using a sonic borehole tool |
US6430507B1 (en) * | 1999-04-02 | 2002-08-06 | Conoco Inc. | Method for integrating gravity and magnetic inversion with geopressure prediction for oil, gas and mineral exploration and production |
US6351991B1 (en) * | 2000-06-05 | 2002-03-05 | Schlumberger Technology Corporation | Determining stress parameters of formations from multi-mode velocity data |
US6751558B2 (en) | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
WO2003036044A1 (en) * | 2001-10-24 | 2003-05-01 | Shell Internationale Research Maatschappij B.V. | Use of cutting velocities for real time pore pressure and fracture gradient prediction |
US6714873B2 (en) * | 2001-12-17 | 2004-03-30 | Schlumberger Technology Corporation | System and method for estimating subsurface principal stresses from seismic reflection data |
US7042802B2 (en) * | 2003-09-18 | 2006-05-09 | Schlumberger Technology Corporation | Determination of stress characteristics of earth formations |
BRPI0415662A (pt) * | 2003-10-24 | 2006-12-19 | Shell Int Research | método de investigar uma região de reservatório em uma formação subterránea |
US7349807B2 (en) * | 2004-03-08 | 2008-03-25 | Geomechanics International, Inc. | Quantitative risk assessment applied to pore pressure prediction |
-
2006
- 2006-08-07 US US11/499,931 patent/US7299132B2/en active Active
- 2006-08-08 AU AU2006278297A patent/AU2006278297B2/en not_active Ceased
- 2006-08-08 MX MX2008001875A patent/MX2008001875A/es active IP Right Grant
- 2006-08-08 CN CN2006800366219A patent/CN101278209B/zh not_active Expired - Fee Related
- 2006-08-08 WO PCT/US2006/030885 patent/WO2007019497A2/en active Application Filing
- 2006-08-08 CA CA2618555A patent/CA2618555C/en not_active Expired - Fee Related
- 2006-08-08 EP EP06800965A patent/EP1913423A2/en not_active Withdrawn
-
2007
- 2007-08-23 US US11/844,303 patent/US7617051B2/en active Active
-
2008
- 2008-02-07 EG EG2008020232A patent/EG24897A/xx active
- 2008-03-07 NO NO20081227A patent/NO20081227L/no not_active Application Discontinuation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5130949A (en) * | 1991-06-28 | 1992-07-14 | Atlantic Richfield Company | Geopressure analysis system |
US6694261B1 (en) * | 1999-06-07 | 2004-02-17 | Conoco Inc. | Method for identification of shallow water flow hazards using marine seismic data |
EP1410072A1 (en) * | 2000-10-10 | 2004-04-21 | ExxonMobil Upstream Research Company | Method for borehole measurement of formation properties |
Non-Patent Citations (3)
Title |
---|
EATON B A ET AL.THE EQUATION FOR GEOPRESSURE PREDICTION FROM WELL LOGS.《THE EQUATION FOR GEOPRESSURE PREDICTION FROM WELL LOGS》.1975,全文. * |
FREDRICH J T ET AL.Stress Perturbations Adjacent to Salt Bodies in the Deepwater Gulf of Mexico.《SPE ANNUAL ECHNICAL CONFERENCE AND EXHIBITION》.2003,全文. * |
SAYERS C M ET AL.Predrill pore-pressure prediction using seismic data.《GEOPHYSICS SOC. EXPLORATION GEOPHYSICISTS USA》.2002,第67卷(第4期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
NO20081227L (no) | 2008-05-08 |
AU2006278297B2 (en) | 2010-06-03 |
AU2006278297A1 (en) | 2007-02-15 |
CA2618555A1 (en) | 2007-02-15 |
CA2618555C (en) | 2013-04-23 |
WO2007019497A2 (en) | 2007-02-15 |
WO2007019497A3 (en) | 2007-07-05 |
EG24897A (en) | 2010-12-13 |
US20070288170A1 (en) | 2007-12-13 |
US7299132B2 (en) | 2007-11-20 |
US7617051B2 (en) | 2009-11-10 |
MX2008001875A (es) | 2008-04-09 |
EP1913423A2 (en) | 2008-04-23 |
CN101278209A (zh) | 2008-10-01 |
US20070038377A1 (en) | 2007-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101278209B (zh) | 钻前孔隙压力预测方法和系统 | |
CN101512100B (zh) | 用于孔隙压力预测的方法及系统 | |
US10788604B2 (en) | Fracturing and reactivated fracture volumes | |
Moss | Quantifying measurement uncertainty of thirty-meter shear-wave velocity | |
RU2687668C1 (ru) | Способ и система комбинированного сопровождения процесса бурения скважины | |
US10528681B2 (en) | Methods and systems for constructing and using a subterranean geomechanics model spanning local to zonal scale in complex geological environments | |
RU2602409C2 (ru) | Инверсия анизотропии многоскважинной системы | |
Roecker et al. | Joint inversion of gravity and arrival time data from Parkfield: New constraints on structure and hypocenter locations near the SAFOD drill site | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
US8243549B2 (en) | Estimating seismic anisotropy of shales | |
Lepillier et al. | Variational Phase‐field modeling of hydraulic fracture interaction with natural fractures and application to enhanced geothermal systems | |
Sultan et al. | Analysis of slope failures in submarine canyon heads: an example from the Gulf of Lions | |
US9051815B2 (en) | Apparatus and method for predicting vertical stress fields | |
Wenning et al. | Image log analysis of in situ stress orientation, breakout growth, and natural geologic structures to 2.5 km depth in central Scandinavian Caledonides: Results from the COSC‐1 borehole | |
Hui et al. | Comprehensive characterization and mitigation of hydraulic fracturing-induced seismicity in Fox Creek, Alberta | |
Gu et al. | Geomechanical-data acquisition, monitoring, and applications in SAGD | |
US20190227192A1 (en) | Drilling Geomechanics Salt Creep Monitoring | |
US20220050224A1 (en) | Evaluating anisotropic effective permeability in rock formations having natural fracture networks | |
Blouin et al. | Evolution model for the Absheron Mud Volcano: from stratified sediments to fluid mud generation | |
GB2523460A (en) | Singularity spectrum analysis of microseismic data | |
Khaksar Manshad et al. | Sand production onset using 3D Hoek–Brown criterion and petro-physical logs: a case study | |
Savage et al. | Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses | |
Perchikolaee et al. | Building a Precise Mechanical Earth Model and its Application in Drilling Operation Optimization: A Case Study of Asmari Formation in Mansuri Oil Field. | |
Blake | Stress analysis for boreholes on Department of Defense lands in the western United States: A study in stress heterogeneity | |
Saberi | Rock physics integration: From petrophysics to simulation |
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: 20111116 Termination date: 20140808 |
|
EXPY | Termination of patent right or utility model |