CN86104849A - 避免钻井设备卡在井中的方法 - Google Patents

避免钻井设备卡在井中的方法 Download PDF

Info

Publication number
CN86104849A
CN86104849A CN86104849.0A CN86104849A CN86104849A CN 86104849 A CN86104849 A CN 86104849A CN 86104849 A CN86104849 A CN 86104849A CN 86104849 A CN86104849 A CN 86104849A
Authority
CN
China
Prior art keywords
well
class
drilling
wells
matrix
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.)
Withdrawn
Application number
CN86104849.0A
Other languages
English (en)
Other versions
CN1011429B (zh
Inventor
W·布伦特·海姆普金斯
罗格·H·金格斯伯罗夫
维斯利·E·洛海克
康罗·丁·尼尼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Original Assignee
Chevron Research Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chevron Research Co filed Critical Chevron Research Co
Publication of CN86104849A publication Critical patent/CN86104849A/zh
Publication of CN1011429B publication Critical patent/CN1011429B/zh
Expired legal-status Critical Current

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B31/00Fishing for or freeing objects in boreholes or wells
    • E21B31/035Fishing for or freeing objects in boreholes or wells controlling differential pipe sticking
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B21/00Methods or apparatus for flushing boreholes, e.g. by use of exhaust air from motor
    • E21B21/08Controlling or monitoring pressure or flow of drilling fluid, e.g. automatic filling of boreholes, automatic control of bottom pressure
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/30Specific pattern of wells, e.g. optimising the spacing of wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B44/00Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mechanical Engineering (AREA)
  • Earth Drilling (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明介绍在同一个地质区中,类似井中相同设备已发生卡钻的深度范围内,避免钻井设备卡钻的方法。在所知深度下,对许多井同时测得许多钻井变量。这些钻井包括机械卡钻、差压卡钻、或两者都有,以及未卡钻等类井对每类井的所有井的所有变量进行多变量统计分析,然后在允许值内改变变量,使钻井的点标位置朝向不使钻柱卡钻的井的平均值移动。

Description

本发明介绍了,在已知钻杆会卡住的地质区域中,测定钻杆被卡的概率的方法。更确切地说,本发明介绍控制或改进钻井条件的方法,以避免以下原因造成的钻杆卡钻:即,钻柱的机械状况和在井筒中诸如高井孔角,过大的钻铤等一类的情况,或在钻杆反抗井筒周围的低压地层时产生过大的静水压差的结果造成差压卡钻。
本发明的特定目的是,利用统计计算钻井时钻杆卡在井筒中的概率和点绘该概率,或两种方法都用,来控制钻井并修正钻井条件,以避免卡钻。该种概率由大量的独立和相关的变量或代表井中机械、化学和水力钻井条件的物理量计算得出。大量井的同一些物理量是在钻柱由于机械状况或差压卡钻的浓度下测得的,或在钻柱未被卡的许多相似的井的相应浓度下测得的。然后,以在发生过卡钻的给定的地质区域的多口井中类似地测得的量,计算统计概率。在此所用的“地质区域”一词,包括钻出了许多口井的沉积盆地的地理区域,和相似地层次序的地区,在所知井深的范围内常常遇到的不同成分组成的页岩-砂地质体。由大量的因井筒内的机械状况和差压条件而造成卡钻的实例井的测定,以及类似的大量钻杆未被卡的实例井的测定,通过依次控制与钻井条件有关的这种测定量,避免由于机械或是差压或两种原因同时造成钻井过程中卡钻的概率是增加了。
按照本发明监测和校正在钻井时测得的机械的和水力学的量,是利用三类这种数据的称之为多变量分析的统计方法来完成的。该种分析依赖于矩阵代数,产生一个在给定深度范围内,代表每一类中所有井的条件的每口井的单一向量。然后,该代数值在所选出的最佳区分三类井的二维平面内,标绘为相应井向量的交点。接着,该多个相关和不相关(但是已测定的和可测定的)变量的统计概率可产生在给定井当前钻井条件的相似向量,用此向量就可测定该井与三类井中各类的相关位置。然后,就改变变量来修改单个井的钻井控制条件,如钻井泥浆特性,井孔角,钻杆组成等,具体修改则取决于这些变量对于井向量相对于三类井中各类的三个空间区代表量的标绘点位置的影响是正还是负。
在用水基钻井液和在未安放防钻杆卡钻的油井套管的条件下钻深井,比如说超过12,000英尺的井,是一个长期存在的问题。特别在海上钻井时,通常都是从工作面积不到1/4英亩的单一固定平台上钻多口深井。因此,必须以偏离垂直线相当大的角度定向钻井以达到离单一平台较大距离之处(“造斜”或“喷射造斜”)。用这种方法,就可以从复盖了包含多个生产间距的相当大的地下区域的地层中产油。
一般来说,使用水基钻井液(能润滑和冲洗井眼中的钻屑)钻这样的井是最经济的,而更准确地说,它可在井筒中产生静水压或静水压头,从而控制在含油层中可能遇到的压力。在钻井时,这种静水压头可防止油或气“喷出”或漏失于井中。进而,钻井液中含有固体物质,该固体物质将在井筒壁上形成薄泥浆饼,从而密封任何被井穿透的可渗透地层。这种水基钻井液体,包括海水,从原始成本,维修保养和保护海洋环境的观点看,实质上比任何一种油基钻井液都便宜。
早已知道,钻柱“卡钻”的一个主要原因是在井筒中的静水压头和钻柱通过的任何多孔低压地层之间存在的差压的影响。在这种情况下,压差以充足的力把钻杆推向井壁以防止钻杆移动。之所以产生这种情况,是因为井筒中钻井液的密度或重量对钻杆造成的静水压力甚大于井筒所通过的多孔地层中产生的静水压力。这是由于流过井筒壁的过滤物(钻井液中的水)和所希望的“泥浆饼”进入低压地层。这种状况可能出现在用于直接从钻头上对钻头施加重量的钻柱的钻铤部分,但似乎更经常的是出现在较浅的深度上,在此深度下,围绕较小直径钻柱流动的返回泥浆几乎不产生涡动,而相对地处于层流状态。因此,在斜井的情况下,钻杆靠近井筒一侧的地方,横过钻杆的较高差压增加了它对井筒壁的附着力。在最坏的情况下,这种差压的结果使钻柱卡住。
纠正钻柱卡钻条件通常需要降低井中钻井液的压力,这可用减少钻井液的静水压头或增加钻井液的固体含量的方法达到,因为这样做就可减少过滤物损失,随之积累起较厚的滤饼从而增加钻杆的接触面积。另外,有时使用较小直径的钻杆或在钻头的权物部件内使用较少的钻铤也可避免卡钻。当在过压地层钻井时,差压卡钻问题常常更严重。在这种井里,地层压力通常超过单由静水压头在该深度上造成的预计压力。在通过过压地层的这种井中,井中平衡的静水压在更深的深度下就不能安全地减小。然而,更深地层上的该种较大压力可能大大地增加地层断裂、随之出现钻井液从井中落入裂缝、产生可怕井喷的危险。
我们还知道,钻柱卡在钻井里常常还由于钻柱和井筒本身之间的机械问题。该种状况有时可能出现在称之为“键槽效应”时,即当钻杆接箍或管子连接处侵蚀出一个环形槽时,钻杆管的大小或工具连接的尺寸超过较大环形井眼一侧的直径,如最初钻头切割时那样,就会产生了键槽。该种槽可能造成钻柱和地层之间的摩擦或阻力大为增加,当试图从井眼中起钻时,造成钻铤卡住,钻铤卡在键槽内。该种问题也可能由于钻杆过重而产生,这时在钻柱的较下部分弯曲,特别是当井眼处于高角度,比如说与垂直方向成60°时,或在钻井平台和钻头之间井眼有不止一次方向变化,如呈S-曲线或形成一个或多个“狗腿”时也会造成。我们还知道,在钻杆的机械卡钻中,井周围的地层可能十分不稳定,以致侧壁坍塌在井筒中,造成卡钻。
据估计,石油工业为钻井时卡钻的花费每年为一至五亿美元,纠正每个事故的花费可能为五十万美元。每个卡钻问题的程度一般依赖于操作人员乐意“冲洗”钻杆卡住部分(松开并取出未卡部分后)的时间量,或者是操纵钻杆打捞“落鱼”的时间量。纷正方法也可包括测定点位或完全用油基钻井液代替水基钻井液。如果无法松开钻柱,就只有弃去井筒或者在被卡点之上侧钻新孔。这样就可能损失钻头,钻铤和井眼内被卡的钻杆。
卡钻问题在许多文献中已有叙述,特别是涉及井筒的差压卡钻,即钻柱贴附于多孔地层以致钻柱一侧周围无钻液循环的文献也很多。如上所述,该种卡钻一般发生在下列情况下:钻井液所含的固体太少,或钻柱和井筒侧之间由于液体从钻井液中漏入到多孔地层中而使液体失去了使泥浆或滤饼稠度增大的控制剂。该类文献主要介绍通过确保钻井液能适合井筒穿过的地层的需要而避免差压卡钻的方法。
在钻深井时,在得不到关于地层的知识、特别是可能遇到低压地层的地方,在卡钻之前就难于预测和采取纠正或预防的行动。进而,当这些问题可用钻柱周围加更深的井眼套管来避免时,这种套管又昂贵,且在一般情况下是不合乎需要的,因为它限制了用通常的测井工具进行地层评价。这也是除非钻井作业必要,而一般不希望用油基钻井液的主要原因。许多地层估计设备或测井设备都依赖于使用水基钻井液,因为水基钻井液通过地层时是导电的,而不是象油基钻井液那样是绝缘的。因为与一般的钻井***相比,预防措施的花费可能过高高,所以,在一切可能的情况下,都希望用普通的水基钻井液钻井,同时仍要避免钻杆卡住。
披露避免或补救卡钻的方法和装置的专利的例子有如下几种:
专利4,428,441中,Dellinger提出使用非圆形或方形工具接头或钻铤,特别是在钻头正上方的钻杆上使用。该种形状确保钻杆周围保持循环,减少差压可能起作用的钻杆与侧壁之间的密封面积。然而,该种工具昂贵,一般也买不到。另外,在比较软的地层它们倾向于使键槽问题加重,因为该种箍的方边在高角度井中可能会切割井的侧壁。
专利4,298,078中,Lawrence提出钻头正上边使用特定的钻进段,当钻杆趋于卡钻时,能用振击器振击钻头。并可以起动工具中附加的阀门来释放钻杆周围的钻井液,从而有助于防止或减缓卡钻条件。
专利4,427,080中,Steiger提出在钻柱的外部粘上一层多孔层,据说,该种涂层由于增加了钻柱周围的液流而防止钻杆的差压卡钻。
专利4,423,791中,Moses披露了在钻井液中加玻璃珠,抑制滤饼在邻近低压区的钻柱和井筒之间形成隔离层,从而避免差压卡钻的方法。
至今,已提出了用统计方法研究解脱钻杆差压卡钻的概率,这种统计分析目标集中解决估计最小浸透时间和最大打捞时间的问题,而此问题的解决可以经济地用于使被卡钻杆松卡。叙述该种方法的文章已在1984年召开的海洋技术会议上发表了,题为“在打捞作业时测定点位液体的时间限制的经济和统计分析”,作者为P.S.Keller等。“卡钻性因子-检查被卡钻杆的新方法”,1983年钻井会议,IADC/SPE论文11383,225-231页,作者为T.E.Love,他使用估计钻井中的一些明显变量,如裸眼的长度,泥浆重量,钻井液的损失和井底装置的长度等的经验公式,进行了评价被卡钻杆松动概率的“卡钻性因子”的统计研究。已利用在墨西哥海弯给定地区的钻井处测得的14个主要参数的互相关方法,从钻杆已被卡的井中和钻杆未被卡的井中导出了该公式。公式的主要用途是测定卡管松动的机会,并仅通过控制经验式中所用变量的选择来指导钻井。还没有人建议对压差卡钻井和机械卡钻井进行统计分析,也没有人建议测定仅仅改变某些测井变量,使上述两种卡钻井状况改变成非卡钻井状况的概率。
M.Stewart已经报导了(1984年在新奥尔良实验室对石油工程师协会,新奥尔良分会的演讲)关于用差压卡钻的统计研究,对在井中、特别是在墨西哥湾遇到的过压地层的井中的特定深度上安放套管,以避免这种地层中不适当的井孔静压水头,或如上所述,避免低压地层破裂的研究。
本发明详细介绍在得不到所遇地层的准确知识的情况下,评价在已知的地质区域中(如上所述)正确地划分正在钻的井或打算钻的井的当前状况和未来状况的概率的方法,然后,控制所选择的任何大量的可变条件或测定量的一个或多个量,如测定钻井液物理和化学特性性,钻栓构型,井眼物理尺寸和井筒穿透的地层。按照本方法,所计算的该种概率再用于更正钻井条件以避免卡钻。然而,如果钻杆已卡住,卡住原因的概率就可测定,直接通过消除该种原因而解脱钻杆,而不是如现有技术中那样只假定钻杆是差压卡钻。
根据本发明,井筒中钻杆卡钻的概率的统计分析指出,卡钻不是象本领域以前的研究人员所指出的那样只是由差压所引起的,而且也由于实质上与差压无关的机械或物理卡钻所引起。在避免钻杆卡钻时,这些条件已发现是同等重要的。特别是,通过对三类井,即出现差压和机械卡钻的井,已完钻而钻杆未卡的井的统计分析,本发明将在指导未来钻井作业时可能作出明显的改进。
为了对钻井进行统计控制,并且有了足够数量的所有三类井,就从给定地质区域的众多口井中,每口井的和井中一定水平上钻柱参数的大量测量数据建立起数据库。这三类井包括钻杆因为(1)机械原因或(2)差压原因而卡钻的井,以及(3)井钻穿过在(1)或(2)类中的深度间隔而没有卡钻的井。选用一种较好的形式,用打点标绘或记录表示每口井数据矩阵解的向量的方式建立了这种概率图。三类井之每一类形成数据矩阵,每类井中每个测定的变量是三个矩阵中每一个数阵(行或列)的一个元素Xij。每个矩阵的大小或是阶数等于每个矩阵中所选的记录变量数V。每个矩阵的补充行或列的大小或阶数则是该矩阵类中所包含的井数N。根据每个矩阵,求出每个变量是相对于所有本类中其它井的同一个变量的标准平均偏差矩阵。根据这些矩阵,就可求出每类井的皮尔逊-积-矩相关系数矩阵,其中所有的系数值介于-1和+1之间。然后,通过称之为多变量判别分析法,解出每个矩阵的这些相关系数的特征据或特征向量。该种分析将这些向量分解成三个明显不同的组,这些组在图形显示时空间上可分离并又能代表在给定地质区域采样的所有的井。
在实行本发明的优先选用的方法中,该种数据矩阵的多变量判别分析,包括找出一个能最佳区分三组中的两组的数学平面。而第三组就被垂直于另一个分隔面的平面分隔出来。因此,两个平面就把三个组彼此分隔开。然后,代表在单个井中许多测定量的完整的一套测定量的每个向量,被投影到垂直于两个面的一个单一平面上,这样,每个井向量表现为一个点,这些点在绘图平面上的坐标与三个向量的间隔相关。根据这些点,可以由每组的重心计算出组间距离,然后,所有所测定的这些值的总的重心,就画在绘图平面上或点在该平面上。根据每口井被正确归类于适当组的计算的概率,可以勾画出正确性概率的等值线。在概率接近相等的地方,这口井属于两组中的任一组,则这个向量一般就处于靠近这两个平面交线之处。因此,一个点离交线越远,这口井被正确地分类的概率就越大。
从概率“图”,然后就可根据同样测定的多个变量,点绘出钻井的进程。通过每个变量元素的系数值的计算,并总合该值以按它的当前钻井深度定出图面上井数据向量的相交位置。然后,按测定的变量的条件改变钻井“概率向量”的控制,以使概率井向量投影的座标移向或超过“从未卡钻”概率的重心。
例如,许多测定变量产生了把当前钻井过程与钻杆的机械卡钻联系起来的井向量时,该种条件很大程度上依赖于偏离于垂直线的角度,井眼直径,钻铤尺寸和井眼的总深度以及钻杆上的摩擦力(阻力)和力矩,但它们也与钻井液的流体力学性质和化学性质有关。当该种向量投影处于基本上相应于钻杆差压钻的高概率的向量空间时,该向量就极大地依赖于钻井液特性,诸如密度(每加仑的重量),粘度,胶体强度,失去性和流速;但也可能与井眼的深度和偏角有关。可引起差压卡钻或是机械问题或能使两者同时出现的其它的测定钻井***变量,诸如真垂直深度,钻井液的PH值和钻井气,用本法估价也是有希望的。在每种情况下,每个测定的变量仅可在它们可用值许可的范围内调节。
因为每口井中多个测定的参数都适当而清楚地描绘该种概率,即在选出的深度范围内,任何井的钻进中落入三类之中,正确的一类的概率,所以对任何要钻的井或正在钻的井,都可控制,以把它们的钻井条件从卡钻的危险处“拉开”,并使之倾向于不卡钻的概率。
在实行本发明的优先选用方法中,每口井产生由测定的多个变量的每一个相对贡献所组成的特性井向量,这个向量可从多维空间投影为一个单值的量,并可以用两个座标点画在选择的两维制图空间上。然后,把它相对于三组井或三组井之每一类的许多口井的位置表示出来。因此,在任何给定的深度钻井时,可以通过多口井在同一个制图空间的向量投影,对这口井作相似的估价。投影于图上的向量的两个座标,最好是同一些变量和相应于图上所有井的相同变量的系数的乘积的总和。然后就可采用改正行动以确保井向量脱离差压卡钻或机械卡钻或两者原因皆有的高概率区,而趋向于井具有高概率不卡钻的绘图区域中的“安全”值。
根据实行本发明方法的最佳形式,在地质区域的每一个单独井的选择深度下测定大量的井变量,以建立数据库。最好的深度应选择实际上发生了或是差压卡钻或是机械原因卡钻的深度。对于没有卡钻的井,就选择被卡井范围内的同一个深度。然后以相应于三类中每一类井的三个独立矩阵排列数据库。在每个矩阵中,一行(或一列)的每个元素相应于一个井内按选择的深度测定的变量。然后,计算每个井中每个数据元素的标准均偏差,以产生三类井中每一类井的标准正常变量矩阵。根据标准正常变量矩阵,通过相应的测定变量的交叉相乘乘,再加上每个矩阵中所有可能井对的交叉乘积,就产生皮尔逊-积-矩相关系数矩阵。许多井的许多这种井向量,就形成为适用于整个地质区域的同样大小的概率矩阵。因此,在该矩阵内的元素包括以下井中产生的元素:
(1)已知为差压卡钻的井,(2)已知为机械问题卡钻的井,及(3)钻柱未卡的井。然后利用在统计学上称之为该类矩阵的“多变量判别分析”方法隔开三个组,该项方法中,由一对彼此正交的数学面分隔开三个组。然后,多维空间中的每个井向量分解为一对系数,这一对系数可在正交于两个面的制图平面上表示为一个点。这就使从多维空间投射向量被最大程度地分隔开,向量则与用两维来点画出的绘图平面相交。通过在绘图平面上勾画出每口井由其向量系数所代表的概率的等值线,就可以区分开差压卡钻井和机械卡钻井,并能区分出两种卡钻向量与“从未卡钻”钻栓向量。然后,用从正在钻的井筒的任一水平产生的同种变量的单个测定量,每个该种变量的系数值用于计算向量系数乘以现时变量值之乘积的和。这些总数产生被控井在图面上的向量座标,并显示出正钻进的井相对于三组井的现在的概率。从该种计算得到的位置看,就可正确地估价并改进钻井中的可控变量,诸如泥浆重量,固体,钻铤尺寸等,从而使钻井的概率移向图上表示井处于我们所希望的高概率“非卡钻区”的座标。该种方法能对钻井进行分析并能定向控制,从而避免钻井中钻杆出现或是机械卡钻或是差压卡钻的问题。
从下列附图的详述和本发明最佳方案的叙述中,就将显示出本发明的进一步目的和优点。
图1是表示从单个海洋平台所钻的多个井的透视横截面立视图,表示出几类深而高偏离的井,本发明的钻井法特别可应用于这些井,以提高避免钻杆或由于差压或由于机械问题卡在井筒中的概率。
图2是部分井筒的透视立视图,说明机械卡钻中涉及的一类问题,即钻杆在井筒侧面形成了小直径的键槽。
图3是部分井筒的透视立视图,说明由于差压钻柱卡在低压地层上。
图4是按图2的4-4箭头方向通过钻柱和井筒的横截面图,指示出键槽中的钻杆。
图5是对给定地质区域由于机械和差压问题卡钻的大量井的调查的条线图。
图6是按图5试样中井的测定深度范围与所有井总发生率的百分比绘制的条线图,总发生率百分比包括机械和差压卡钻的发生率和不卡钻的发生率的百分比。
图7是相似于图5和图6的条线图,表示按机械和差压卡钻的总发生率的百分比绘出的井眼尺寸范围。
图8是卡钻概率“图”,在此图中,每个井的向量表示为多维空间中该向量与二维表面的交点。二维表面正交于分隔代表三类井的三个空间向量群的两个面,这三类井是(1)机械卡钻的井,(2)差压卡钻的井和(3)未卡钻的井。
图9是卡钻概率“图”,在此图中,按每口井被正确地分类入正确组的概率勾画出等值线。
图10是单口井的进程图,井是按规则的深度间隔通过取样的变量进行分析的,该井已差压卡钻。图显示井进行的过程是:从未卡钻概率,通过机械或差压卡钻的概率,一直到钻柱应当差压卡钻的高概率结尾状态,而实际上也就发生了差压卡钻。
图11是图9所示的井向量的三角图。
图12是三类井中,每类的三口井的四个可测定变量的说明性的例子所产生的,通过计算机程序计算的井向量图。
图1分别用立视和部分透视的方法,表示一个固定的海洋钻井平台10,通常,这类平台作为开发一个或多个水下生产地层的主要部分。本发明的钻井控制***特别应用于这种钻井,因为,有许多口井,比如说10-30口井,象11、12、13、14和15那样由单个平台10上以对垂直线的高偏离角钻的井,以开发从平台位置横向延伸数千英尺的水下贮油层16。如图所指,11至15号井是按不同的角度选择性钻的,可包含一个或多个“狗腿”17(与垂直线偏离角不同),在钻到所需深度时,它们甚至可以采取S-型的构型,如14号井的情况。该种构型或因地质条件而事先设计的,或在钻井时偶然出现的。
我们早已知道,高角度井有卡钻的趋向。特别在钻井深度超过12,000英尺时尤其是这样。一般都假设,这种卡钻是由于井筒和作用于钻杆上的地层之间存在差压;该种差压是由于井筒中的压力高于井筒所穿过的地层的压力。在某些地质区域,包括墨西哥海弯的海上井,在相当线的深度常常遇到高压,也即是,该地层的压力超过在这个深度下预料的静水压头或地层静压头的标称垂直梯度(标称井压实质上是在给定深度下井筒的水压)。为了控制过压地层,由井筒中的钻井液或泥浆施加于井眼的井压必须超过地层的压力。然而,在井的更大深度上,地层压力可能接近于该种深度的标称值。因此,为了保持对应于上部的超高压的适当井压,对较低地层的静水压可能过压。该种过高的井压可能使地层破裂,结果,钻井液漏失到地层中,产生井喷的危险。
利用水基钻井液,通过低压可渗透地层钻过压井时,水可能流入地层。该种流动通过井筒21周围的井壁泥浆或滤饼20,这种滤饼一般是堵塞可渗透地层23的薄层胶化固体。这种流动可引起滤饼内固体的过多沉淀。这种状况在图2和3的22处所示。液体连续流入地层就增加滤饼的厚度和增加钻杆17的接触面积,从而使钻杆堵在或卡在井筒壁17处。滤饼厚度的增加使钻井液恢复在钻杆和井筒之间循环的趋势就变得困难。进而,摇溶性的钻井液从钻头回流到表面。再流过井眼21的剩余区,流动会变得相当平稳,从而导致钻井液固结或结胶。如钻井工艺中熟知的那样,该种差压卡钻的准确原因常常难于确定,因此,改正这种状况一般都是通过反复试验、不断摸索进行的。
进而,改正卡钻状况的前景可以确定相对于放弃井筒这一部分的花费而言,操作员有多少钻机时间能花费在“打捞”上。放弃井筒常常需要在未卡钻的最后管段周围井壁上另钻新眼。这就要求安置一个柱塞,这会带来设备的损失,然后再钻到同样深度。因此,知道避免卡钻或解除差压卡钻的概率,以及知道钻柱是机械卡钻而不是差压卡钻的概率,具有很高的经济价值。这特别对于海上钻井时、钻机每小时的花费在数千美元的情况更是如此。
图2和4表示钻铤25和钻头27之上的钻杆17部分。如图所示,整个钻杆17的直径基本上都小于钻头27开始钻切的井眼21的直径。总的来说,钻杆本身其挠性大于包括钻铤25和钻头27的井底装置。因此,在高角度下,钻杆可能倾向于向井筒壁一侧下垂。在这种情况下,钻柱可能机械切割井筒侧壁,如图2和4的29所示,形成所说的“键槽”。在这种情况下,钻杆17的直径或管段间的连接段的直径小于钻铤部分或钻头。当钻杆上下移动时(如在钻柱的“一次起下钻”后改变钻头),钻杆和连接段就可能由于机械原因而卡在井眼内。
低压地层坍塌于井筒中可以产生其它的机械问题。虽然都知道,钻柱卡钻可以由于差压条件也可由机械问题而发生,但一般都假定最大的危险是差压卡钻,并且,过去的实践一般也都假定,任何卡钻井都是差压卡钻。
从我们对许多钻杆卡钻情况的统计研究发现,上面的假定未必正确。因此,试图使钻杆解脱卡钻状态的方法就可能是专门针对一种最可能的原因,如或是机械的,或是差压的,或两者皆有的原因。因此,测定钻杆为何卡住或可能要卡住的概率的方法,和如何避免钻井卡钻的方法是钻井中所盼望解决的。
我们的研究包括从数百口井测得的钻井变量,其中一些井已知是差压造成卡钻的,另一些所知或怀疑是机械问题造成卡钻的井。然而,在同样地质区域中,还有大量井在钻进过程中未出现钻柱卡钻。所有井都是在墨西哥海弯重要的地质区域钻的。大体说来,在该地质区域采样的井包含在一般具有相似地质结构的盆地中钻进深度超过12,000英尺的井。该类井钻穿过形成贮油层的砂和页岩地层,诸如包围盐丘或由于断层中断的地层。
如以下要更充分地解释的,在每口井中测量钻井变量。在三类井中每类井有许多个井的每一个井中的选择深度上,记录数十个测定的和可测定的变量中的20个。在图5、6和7中,表示出三类井中每类井的相对数。图5用条线图的方式表示在与垂直线偏离0°-75°范围内的采样井中机械或差压卡钻的百分数。图6用条线图的方式绘出形成数据矩阵的三类井的分布,都作为井深度的函数而绘出的。图7也是相似的条线图,按取样井中的孔径范围绘制。
图8、9和10是三类井中每类井的各口井的、在单一平面的向量投影图或分布概率图。这些图是利用本发明的方法,通过三类井中每类井全部测定变量的多变量分析而绘制的。这些图显示,这三类井容易以充分高的概率区分开,所以,通过在任何给定深度上测量同样一行待测变量,就可以绘制单一钻井的钻井状况,从而控制正被开钻的井。该种控制可以是由预先设计钻井计划来进行,或是在钻井时执行改正行动而进行。在图10的二维图上,绘制了显示钻井时该井与相对于三类条件的进展。
如图8、9和10所示,绘制的图上的点子对于控制是有用的,这些图是利用称之为多变量判别分析法对概率进行统计分析而绘制的。在给定的地质区域,三类井每一类的大量井用于形成统计分析的可靠样本。然后利用在赋值的矩阵中的每口井的同样多的变量,形成可比较的数据矩阵。对于本领域的熟练人员来说,显然也能按照钻井人员的需要,根据大量有明显差别的所测定的钻井变量,对其它地质区绘制相似的概率图。
在图8中,用在图中心相交的三条线来表示三个组被二个彼此正交的平面分隔开。两个平面正交于绘图平面。
图9相似于图8,表明三个组内每组的等值线,这些等值线代表每个井向量被正确的绘制在指定组内的概率。图10为向量系数图,图中所点画的井同于图8和图9中点画的井。
图11用三角座标图的形式,表示图9中所示的三类井中每类井的概率的另一种绘图方法,如所指,每口井越接近每类井的顶点,为通过改变起作用的变量而进行改正所作的正确分类的概率就越大。
例子
为了说明本发明的方法的工作进程,下面计算一个简化的例子。共有三组或三类中,每组三口井,每口井总共四个测得的变量。很显然,在实际应用时,对每个矩阵中所有的井,比如说是40到100口井,及每口井的所有变量,比如说20个测得的变量,所用的也就是这个计算方法。
如上所指,为了区分出三个组,井的选择是基于每口井的已知深度上的一组20个变量。在每个卡钻的情况下,这一组数最好是这种变量的最后一组,即钻柱因机械原因或因差压而卡住的深度上的一组变量。然而,还要用钻柱刚刚卡住之前的这口井的测定的条件。每个没有卡钻的井的单组20个变量,可以在随机选定的深度上选取,这个深度要在差压卡钻或机械卡钻的典型深度范围之内。
然后,用V个变量和N口井组成矩阵x,组成方式如下面三个矩阵,每个矩阵四个变量三口井的例子:
三组中第一组,三口井,四个变量:
变量V=
井N=    i=1    i=2    i=3    i=4
J=1 〔x11=〕9750 13.7 4750 70.0
J=2    9500    14.5    5000    60.0
J=3 10000 13.1 4500 〔xij=〕50.0
式中从i=1到i=4各列中的V个变量:
i=1    是总深度(英尺)
i=2    是泥浆比重(磅/加仑)
i=3    是井底处的钻探重量(磅)
i=4    是钻孔对垂直线的偏离角度(度)
然后,从每个元素中,如从x11中,减去列的平均值xi,则得到每一列的零均值。
在本例子中,每一列的列均值由下式确定:
Xi= 1/(N) Σ i = 1 N x i 或者
Xi=1= 1/3 (9750+9500+10000)=9750
对其余各列做类似的计算,则本组的几个列均值是
9750.0000    13.7666626    4750.0000    60.0000000
然后来计算每一列的标准偏差,算法是先算出每列中每个元素对列均值的偏差的平方,再于一列中把这些平方相加,并被变量数减一来除。每一列的这个和的平方根就是标准偏差Si
上述例子中计算标准偏差的过程如下:
第一列数据的方差计算如下:
方差= 1/(N-1) (9750-9750)2+(9500-9750)2+
(10000-9750)2=62500
(如下表中所用的,62500就是0.625×105,表示为0.625E+5)。标准偏差就是方差的平方根,为250.99)。这个值由计算机算出时表示为:249.927994,在数据精度范围内,和250.00是一样的。类似的,各列的标准偏差为:
249.927994    0.7024302    250.007996    10.0000000
为了表示变量之间的任一线性关系,协方差按下式计算:
1/(N-1) Σ i (X i j - Xj)(Xik- Xk)
式中i代表井号,j、k从1到4,表示变量。当j=k时,该乘积就是方差。
则方差-协方差矩阵如下:
变量→
↓1    2    3    4
1    0.625E+05    -0.175E+03    -0.625E+05    -0.125E+4
2    -0.175E+03    0.493E+00    0.175E+03    0.300E+01
3    -0.625E+05    0.175E+03    0.625E+05    0.125E+04
4    -0.125E+04    0.300E+01    0.125E+04    0.100E+03
当对角线各元素被该变量的方差除后,该值就恒定为1。不在对角线上的元素被行一列交点所代表的变量的两个标准偏差的乘积来除,即第一行与第二列相交,就被变量1和变量2的标准偏差来除,这就给出了相关矩阵。
相关矩阵为:
变量→
↓1    2    3    4
1    0.100E+01    -0.996E+00    -0.100E+01    -0.500E+00
2    -0.996E+00    0.100E+01    0.997E+00    0.427E+00
3    -0.100E+01    0.997E+00    0.100E+01    0.500E+00
4    -0.500E+00    0.427E+00    0.500E+00    0.100E+01
这个矩阵是为对角线对称的,即第一行与第二列的交叉点和第二行与第一列的交叉点是一样的。相关矩阵有特殊的性质,它是正的,半确定的(即所有的特征根都是非负的)。
另外两组有以下的统计特征值:
三组中的第二组,三口井,四个变量
原始数据:
井    1    2    3    4
1    5500.0000    10.80000    3700.00000    21.00000
2    5000.0000    10.40000    3500.00000    25.00000
3    6000.0000    11.20000    3250.00000    30.00000
这一组的平均值是:
5500.0000    10.7999973    3483.33325    25.3333282
这一组的标准偏差是:
500.023926    0.4000427    225.459534    4.5092545
方差-协方差矩阵:
变量→    1    2    3    4
↓1    0.250E+06    0.200E+03    -0.625E+05    0.125E+04
2    0.200E+03    0.160E+00    -0.500E+02    0.100E+01
3    -0.625E+05    -0.500E+02    0.508E+05    -0.102E+04
4    0.125E+04    0.100E+01    -0.102E+04    0.203E+02
相关矩阵:
变量→
↓1    2    3    4
1    0.100E+01    0.100E+01    -0.554E+00    0.554E+00
2    0.100E+01    0.100E+01    -0.555E+00    0.554E+00
3    -0.554E+00    -0.555E+00    0.100E+01    -0.100E+01
4    0.554E+00    0.554E+00    -0.100E+01    0.100E+01
三组中第三组,三口井,四个变量
原始数据:
1    7000.00000    12.10000    3875.00000    35.00000
2    7250.00000    12.00000    3000.00000    48.00000
3    8000.00000    12.80000    3950.00000    40.00000
本组的平均值:
7416.66406    12.2999926    3941.66650    41.0000000
本组的标准偏差:
520.453613    0.4361027    62.8649292    6.5574389
变量→1    2    3    4
↓1    0.271E+06    0.213E+03    0.115E+05    0.375E+03
2    0.213E+03    0.190E+00    0.625E+01    -0.699E+00
3    0.115E+05    0.625E-01    0.395E+04    0.400E+03
4    0.375E+03    -0.699E+00    0.400E+03    0.430E+02
相关矩阵:
变量→
↓1    2    3    4
1    0.100E+01    0.937E+00    0.350E+00    0.110E+00
2    0.937E+00    0.100E+01    0.228E-02    -0.245E+00
3    0.350E+00    0.228E-02    0.100E+01    0.970E+00
4    0.110E+00    -0.245E+00    0.970E+00    0.100E+01
这些矩阵再相加到一起,得到各组所有井的(W)组内联合矩阵:
联合矩阵:
W    MAT部分1
变量→
↓1    2    3    4
1    0.117E+07    0.475E+03    -0.227E+06    0.750E+03
2    0.475E+03    0.169E+01    0.250E+03    0.660E+01
3    -0.227E+06    0.250E+03    0.235E+06    0.127E+04
4    0.750E+03    0.660E+01    0.127E+04    0.327E+03
井的总数=9
各组所有井的全部统计数据综合在一起,就是:
总样本的均值
7555.5547    12.2889    4058.3333    42.1111
总样本的标准偏差:
1882.3816    1.3643    581.2178    16.3359
总相关矩阵
TOT    R部分1
变量→
↓1    2    3    4
1    0.100E+01    0.943E+00    0.905E+00    0.904E+00
2    0.943E+00    0.100E+01    0.927E+00    0.902E+00
3    0.905E+00    0.927E+00    0.100E+01    0.892E+00
4    0.904E+00    0.902E+00    0.892E+00    0.100E+01
所有井的总均值周围的组间距离计算出来为:
变量→
↓    1    2    3    4
1    0.272E+08    0.189E+05    0.815E+07    0.222E+06
2    0.189E+05    0.132E+02    0.563E+04    0.154E+03
3    0.815E+07    0.563E+04    0.247E+07    0.665E+05
4    0.222E+06    0.154E+03    0.665E+05    0.181E+04
求得总相关矩阵的特征向量:
特征值1    73.3556061
特征值2    0.2083998
并进行检验以确定结果精度(所有检验都应是同样的值):
特征值之和=73.5640259
B-1/2′AB-1/2的迹=73.5639648
W-逆矩阵A的根:
73.3556    0.2084
W-逆矩阵A的迹=73.56403
并且,每个特征值所说明的数据中变化的百分比,相加应为100%:
每个根的百分比:
99.7167    0.2833
判别函数的计算为:
W-逆矩阵A的向量,作为例:
向量部分1.
变量→
↓    1    2
1    0.244E-02    0.139E-03
2-0.100E+01    -0.100E+01
3    0.492E-02    0.206E-02
4    0.274E-01    -0.809E-02
下面可以给出特征值和判别函数的偏差的筒单解释:
取某个矩阵Q,解行列式方程:
|Q-λI|=0
式中I是单位矩阵,λ是特征值。
求特征值和特征向量:
|(Q)-(λI)|=0
|(1 2 3 2)-(λ 0 0 λ)|=0
1.求得特征值:|1-λ 2 3 2-λ|=0=(1-λ)(2-λ)-6
或:
=λ2-3λ-4=0
因之我们得到:λ1=4
λ2=-1
2.把特征值代入,就求得有关的特征向量:
a.λ1=4
Figure 86104849_IMG4
注意,系数矩阵的秩为1,这就意味着存在一个线性独立的解向量,其他解都是这个解向量的倍数。
经检验,C11 1)是这个向量。
b.λ2=-1
只存在一个解向量:C23 2)。因之特征值是4和-1,特征向量分别是C11 1)和C23 2)。
特征向量可以考虑为判别函数,并且是适当正规化的判别函数。
这个例子不具有同相关矩阵一样的性质,因为有一个特征值是负的。之所以选择一个简单的矩阵例子,是因为原来给出的三组的例子的矩阵太复杂,不能用袖珍计算器快速解出。
在得到特征向量之后,要对特征向量进行定标,以表明每个变量对判别函数的相对重要性。
定标的向量
定标部分1
变量→
↓    1    2
1    0.264E+01    0.150E+00
2    -0.130E+01    -0.130E+01
3    0.238E+01    0.996E+00
4    0.495E+00    -0.146E+00
用Wilk的Lambda准则(λ准则)和F比值来对显著性进行统计检验。
H2=0.0111295检验的λ
F1=8.0000000
F2=6.0000000
对于H2的检验,F=6.3592415
这是在0.01概率水平下的显著性。每口井的判别值的计算是把原始数据乘以每个变量的判别系数,再加上每组中每口井四个变量的结果:
原始数据乘以特征向量一第一组井
N    1    2
1    35.370758    -3.142392
2    34.916916    -3.382162
3    34.803467    -2.860050
原始数据乘以特征向量一第二组井
N    1    2
4    21.395081    -2.596268
5    19.700882    -2.709401
6    20.248352    -3.924898
原始数据乘以特征向量一第三组井
N    1    2
7    24.999207    -3.441051
8    26.679733    -3.154366
9    27.245026    -3.888223
这就完成了主要的判别分析。
由以下来计算正确分类的概率:
检验空间中各组的均值
9750.00000    13.7666626    4750.00000    60.0000000
5500.00000    10.7999973    3483.33325    25.3333282
7416.66406    12.2999926    3941.66650    41.0000000
判别空间中各组的重心,按行给出
35.0303802    -3.1281977
20.4481049    -3.0768538(三组的联合
36.3079834    -3.4945393    判别均值)
第一组在判别空间中的离散或标准偏差:
0.0901396    -0.0185371
-0.0185374    0.0683792
第二组在判别空间中的离散
0.7482136    0.1753250
0.1753258    0.5427456
第三组在判别空间中的离散
1.3636608    -0.1567893
-0.1567892    0.1366703
利用X2近似贝叶斯统计,求得概率值:
各组的X2值 正确分类的概率
1    2    3    1    2    3
1    1.334    322.918    76.613    1.000    0.0    0.0
2    1.334    307.021    64.589    1.000    0.0    0.0
3    1.331    295.166    74.808    1.000    0.0    0.0
4    2142.553    1.332    18.637    0.0    1.000    0.0
5    2722.738    1.333    32.018    0.0    1.000    0.0
6    2652.085    1.333    37.634    0.0    1.000    0.0
7    1203.734    31.762    1.335    0.0    0.0    1.000
8    820.693    56.615    1.337    0.0    0.0    1.000
9    758.760    73.265    1.333    0.0    0.0    1.000
图12表示了与其特征值相一致的这些组的结果,共9口井,都按其特征向量的座标点上去。三组的分离是明显的。
最佳方式
由前面的例子可以看到,三组中每组有40到100口井,每口井中每一个深度上有20个或20个以上的测量变量,每口井的分类计算和图形表示最好用计算机来完成。
井和变量都不多时,用HP35(Hewlett    Pa-ckard公司产品)袖珍计算机来计算无量纲的矩阵系数。然而,对大的数据组,比如说三个矩阵,每个矩阵是20个变量和80口井,用称为SAS的程序可以进行上述统计计算,这个程序可以从SAS研究所Raleigh,N.C.处得到。这个程序能执行多参数分析的每一步,包括主元素的矩阵计算,因子,回归和判别分析。此外,W.W.Cooley和P.R.Lohnes写的教科书“行为科学的多元参数过程”(John    Wiley    and    Sons,New    York    NY,1962),给出了统计分析的FORTRAN编码。三类井的图形表示和每个井向量的位置可以用Lotus    1-2-3程序点出,这是个可以从剑桥MA,Lotus    Development处买到的程序,这个程序要和dBASEⅡ程序一起用,以便管理数据文档,dBASE    Ⅱ可以从Ashton-Tate,Culver城,CA处得到。计算每个单独的井向量,以便点图并控制钻井的线性程序,可以用程序OMNI(可以从Harerly***公司,Denville,N.J.得到)来执行,程序MPSX(可以从IBM公司,White    Planins,NY处得到)也可以用。
本发明的方法在野外现场使用时,一般测量下列的井变量或井参数。
(1)测井深度
(2)实际垂直井深
(3)裸孔(未下套管)深度
(4)旋转钻杆钻进转矩
(5)旋转钻杆阻力
(6)测量孔角度(由垂直线起算)
(7)钻井液(泥浆)重量
(8)钻井液塑性粘度
(9)钻井液屈服点
(10)钻井液19秒胶粘力
(11)钻井液10分钟胶粘力
(12)API(美国石油学会)标准钻井液失水量(滤液)
(13)钻井液pH值
(14)钻井液绿泥石含量
(15)井孔尺寸(直径)
(16)钻井液固体百分比
(17)钻井液含水百分比
(18)钻井液流量(泵送率)
(19)钻铤外径
(20)钻杆上钻铤段的垂直长度
还成功地使用过钻井液的气体含量,气体类型等各种测定变量。
通过上面叙述,很清楚的是,最好是用三组井确定卡钻的概率,本方法很明确地适用于只区分成两组。这两组可以是卡钻的井和未卡钻的井,或是不加控制的井或加以控制的的井。另外,该分析还可以只用来区分机械卡钻和差压卡钻。井变量既对特定深度的井向量起作用,又和整个井的状态有关,井变量的校正作用由每个变量单独的系数表示。

Claims (11)

1、当用一种水基钻井液钻井时,为避免或机械卡钻或差压卡钻而改变钻井条件的方法,其特征在于,
(a)在许多(N口)不同的井中,在每口井的井筒中给定深度上,测量相应的许多(V个)有关的钻井变量,而井筒中钻杆处于以下三种状况:
(ⅰ)钻井时发生了机械卡钻,或
(ⅱ)由于井筒和被井穿透的可渗透地层之间的压力差而卡钻,或
(ⅲ)已钻过了(ⅰ)和(ⅱ)所选定的井的深度间隔,而没有发生卡钻。
(b)把(a)步中三组N口井每一组形成一个矩阵,在矩降中,每个所说的测量变量V都是一个公用数阵(行或列)中 Xij的元素,这个矩阵包括所说的N口选出做为代表组的成员井的补充数阵(行或列)。
(c)首先形成零均值偏差矩阵,其每个元素都是各个被测变量元素 Xij和零均值 Xj之差,再将每个偏差矩阵元素转换成为它相应的标准平均偏差Sj,从而形成每个变量数阵中每个元素的标准平均偏差矩阵:
所说的平均变量 Xj是: Xj= 1/N Σ i = 1 N X i j
式中j=1,2,……V(变量)
    i=1,2,……N(井数)
平均方差= 1/((N-1)) Σ i = 1 N (X i j - Xj)2
所说的每个元素的标准平均偏差Sj是:
Sj=( 1/((N-1)) Σ i = 1 N (X i j - Xj)2 ) 1 2
(d)由所说的标准平均偏差系数矩阵形成皮尔逊积一矩相关系数,任何两个X之间的值,如Xij,j=1,2,由下式定义:
它用来表示所说的一对X(如j=1,k=2)的线性关系或线性相关,使每个所说的、用一个方形对称阵表示的系数R12来计算一个组的总体RT,如:
Figure 86104849_IMG3
式中各个j和各个k相应于总体每个变量的均值,每个组相关也类似地定义,使各个i只相应于该组的元素,各个 X则只相对于该组的平均值。
(e)然后,再形成所说的系数矩阵R,该矩阵中每个元素都是所说的标准平均偏差矩阵的一行与一列的交叉相乘之和,所说的矩阵是对称的,万形的,是半确定的。
(f)解三个矩阵RT的每一个,求在下式中的矩阵尺度向量系数γi
(RTiI)i=0
式中λi是特征值(特征根),I是单位矩阵。
(g)对在(b)步骤中形成的原始矩阵中每一个原始的测量变量元素乘以该元素相应的特征向量系数γi,再用λi来定标,并分别把每个变量数阵的乘积相加,并且
(h)对每个数阵,把所说的带上γi值并被λi标定过的乘积之和点图,用以表示井已被正确地放入它应属类别的概率。
(i)然后,对在所说的地质区域及深度范围内要钻的、并且要确定其卡钻概率的另一口井中测得每个变量的γi和λi乘积并相加,且
(j)把所说的井的所说的座标的点子用可见的方式显示出来,以表示它相对于所说的三类井中每一组的组均值的卡钻的概率。
2、根据权利要求1的方法,其特点在于,在所说的另一口井中所说的大量测定的参数,按照所说的变量通过特征向量系数所起的作用而加以修改,以改变所说的另一口井相对于所说的组均值的位置。
3、按照避免或机械原因,或因井筒与井筒所穿过的可渗透地层之间的压力差的原因而卡钻的概率,来予测钻井的方向的方法,其特征在于,这个方法包括:
对在类似的地质区域中的至少两类井,形成相关系数矩阵,所说的两类井包括由各组中选出的井,各个组由机械卡钻的井(1),差压卡钻的井(2),以及没有卡钻的井(3)组成。
对于在各个类里的每一口井,所说的每一类的相关矩阵,包括许多个相关于所选深度间隔的基本相同的变量,
每一类井形成许多个单值向量,每个向量在其各自的矩阵中表示一口井,
确定把所说的至少两组向量分开的第一个平面,从而第二个平面垂直于该第一个平面,并限定一个二维的绘图表面,该表面一般地以总均值为中心,表面用于点绘出由所说的两类到每一类的重心或均值的向量的投影,所说的重心建立了每口井向量被正确分类的概率,
然后在正在钻进的井的所选深度上测定同一些变量,
产生一个表示每个所说的变量与所说的垂直面和与重心投影之间关系的一个单一向量,
所说的相关于所说的垂直平面上的所说的总均值的位置,由每个测得的变量的系数值的总和而定,并且
修改所说的井中一些所选的测量变量,以使所说的单一向量离开卡钻的概率重心。
4、根据权利要求3的方法,其特点在于,所说的至少两类井的每一类井被一个垂直于所说的绘图表面的平面所分开。
5、根据权利要求4的方法,其特点在于,三类井由一个同时垂直于所说的绘图表面和两类的分隔平面的附加平面所分开。
6、一个确定钻进过程中钻杆卡在井筒里的概率的方法,所说的卡钻概率由代表所说的钻杆在所说的井筒中的基本上所有钻井条件的许多可测变量的测量所建立,这些钻井条件包括机械的和钻井液体的,与该钻进有关的量,其特征在于,本方法包括:
由多口所钻的井建立一个地质区域的数据库,这些井至少包括三类井,即机械卡钻或压力差卡钻的井,和没有卡钻的井,
所说的数据库,是在所有这些井里,在给定的深度上基本上同时测量的多个变量的组合矩阵中每口井的井向量解。
把所说的向量的每口井点绘为平面上一个点的座标,每个所说的向量是多个变量中每个变量相对于数据库中其它井向量的相对贡献的总和。
7、根据权利要求6的方法,其特点在于,每个所说的三类井中的井向量,在所说的平面上,用多变量分析法分隔开,要用在所说的平面表面上正交的两个平面把所说的组最佳地分隔开。
8、一种对多个测定的钻井参数进行多变量分析的方法,用以控制所说的变量,以减少钻井过程中卡钻的概率,其特征在于,这个方法包括:
以矩阵形式记录在许多口井的给定深度下测定的同样的许多变量,所测的井包括对本组井选出来的至少两类井,本组井包括:钻柱(1)没有卡住,(2)因压差而卡住,(3)机械卡住的井,
对所说的矩阵之内的每口井,确定所说的多个变量中每一个的特征向量的贡献,
把所说的对于每个特征向量的贡献,与在钻进的井中所说的许多变量的测量值相乘,把乘积相加,从而产生当前井向量相对于所说的两类井的平均值的座标,
把井向量相对于所说的两类井的平均值点画出来,以表示所说的井相对于所说的至少两类井的平均值的当前井条件。
9、根据权利要求8的方法,其特征在于,三类井的每一类都分别记录为分离的矩阵,所说的井向量就是相对于三类井的每一类的平均值和所说的这三类井的总平均值而点出的。
10、在给定的地质区域内定向钻井时,避免钻柱卡在井筒中的方法,其特征在于,该方法包括:
对在类似的地质区域中至少两类井形成相关系数矩阵,这两类井包括由一组井中选出的井,这一组井包括钻杆(1)机械卡钻的、(2)差压卡钻的和(3)没有卡钻的井。
所说的每一类井的相关矩阵包括该类中每口井的,相应于所选深度间隔的大量基本相同的变量。
每一类井形成许多单值向量,每个向量在其各自的矩阵中代表一口井,
先确定把至少两组向量分开的第一个平面,从而第二个平面垂直于头一个平面,因之就限定了一个两维的绘图平面,一般以总的平均值为中心,平面用以点绘出从两类中的每一口井到每类的重心或平均值的向量的投影,所说的重心就确立了每口井的井向量被正确归类的概率,
然后在一个正在钻进的井中所选的深度上测量同一些向量,
产生一个单一向量,用来表示所说的每个变量和垂直平面及和重心的关系,
所说的在垂直平面上相对于总平均值的位置,由每个测定的变量的系数值之和来确定,
修改所说的井中的一些所说的测定变量,以使所说的单一向量离开卡钻的概率重心,
11、一种测定钻井时卡住钻杆的统计概率,和按所说的概率修改钻井条件,以求避免钻杆卡在井中的方法,其特征在于,该方法包括:
用在一个地质区域钻的多口井,来计算钻柱(包括钻头)和井筒之间的许多测定的机械条件和井筒中所用的钻井液的测定性质之间的统计关系,
所说的多口井包括第一类是机械卡钻的井,第二类是井筒的钻井液的压力和井筒穿透的地层的压力之差所造成卡钻的井,第三类是设有卡钻的井,
在所说的第一,第二和第三类井的每一类井中,分别计算测得的同样多个的机械条件与钻井液特性的统计关系,
对三类井中每一类井的基本上所有的井的测定的条件,进行多变量统计分析,来确定一个使三类井要彼此分离开的一个绘图平面,
把三类井中每一类的每一口井的井向量,投影到绘图平面上,然后,在被钻井的给定深度上,测量同一些变量条件,
根据整个一组井的每个变量的系数乘以相应的测定的条件的值的乘积之和,把每口井点到图上,以确定根据正在钻的井的测定的条件,把这口井归于三类井中某一类的概率。
CN86104849A 1985-07-15 1986-07-15 避免钻井设备卡在并中的方法 Expired CN1011429B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US75630785A 1985-07-15 1985-07-15
US756,307 1985-07-15

Publications (2)

Publication Number Publication Date
CN86104849A true CN86104849A (zh) 1987-01-14
CN1011429B CN1011429B (zh) 1991-01-30

Family

ID=25042916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN86104849A Expired CN1011429B (zh) 1985-07-15 1986-07-15 避免钻井设备卡在并中的方法

Country Status (9)

Country Link
US (1) US4791998A (zh)
EP (1) EP0209343B1 (zh)
CN (1) CN1011429B (zh)
AU (1) AU608503B2 (zh)
CA (1) CA1257701A (zh)
DE (2) DE209343T1 (zh)
DK (1) DK334286A (zh)
ES (1) ES2000508A6 (zh)
NO (1) NO862850L (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101353959B (zh) * 2008-09-10 2012-08-29 杜书东 防溜钻、遇阻卡自动报警保护装置
CN102979501A (zh) * 2010-04-12 2013-03-20 国际壳牌研究有限公司 自动地提升钻头离开地下地层中的开孔底面的方法
CN106150476A (zh) * 2015-04-09 2016-11-23 中国石油化工股份有限公司 一种预测钻柱的粘吸卡钻风险的***

Families Citing this family (52)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU608503B2 (en) * 1985-07-15 1991-04-11 Chevron Research And Technology Company Method of avoiding stuck drilling equipment
NO893113L (no) * 1988-08-03 1990-02-05 Chevron Res Fremgangsmaate og anlegg for aa unngaa fastkiling under broennboring, saerlig basert paa sannsynlighetsberegninger.
US4972703A (en) * 1988-10-03 1990-11-27 Baroid Technology, Inc. Method of predicting the torque and drag in directional wells
US4848144A (en) * 1988-10-03 1989-07-18 Nl Sperry-Sun, Inc. Method of predicting the torque and drag in directional wells
US5044198A (en) * 1988-10-03 1991-09-03 Baroid Technology, Inc. Method of predicting the torque and drag in directional wells
US5660239A (en) * 1989-08-31 1997-08-26 Union Oil Company Of California Drag analysis method
US5181172A (en) * 1989-11-14 1993-01-19 Teleco Oilfield Services Inc. Method for predicting drillstring sticking
JPH03201067A (ja) * 1989-12-28 1991-09-02 Nissan Motor Co Ltd デザイン装置
IE910209A1 (en) * 1990-02-28 1991-09-11 Union Oil Co Drag analysis method
US5508915A (en) * 1990-09-11 1996-04-16 Exxon Production Research Company Method to combine statistical and engineering techniques for stuck pipe data analysis
US5861362A (en) * 1992-01-06 1999-01-19 Blue Diamond Growers Almond shell additive and method of inhibiting sticking in wells
US5448911A (en) * 1993-02-18 1995-09-12 Baker Hughes Incorporated Method and apparatus for detecting impending sticking of a drillstring
US5316091A (en) * 1993-03-17 1994-05-31 Exxon Production Research Company Method for reducing occurrences of stuck drill pipe
US5327984A (en) * 1993-03-17 1994-07-12 Exxon Production Research Company Method of controlling cuttings accumulation in high-angle wells
GB2279381B (en) * 1993-06-25 1996-08-21 Schlumberger Services Petrol Method of warning of pipe sticking during drilling operations
FR2732403B1 (fr) * 1995-03-31 1997-05-09 Inst Francais Du Petrole Methode et systeme de prediction de l'apparition d'un dysfonctionnement en cours de forage
FR2768818B1 (fr) * 1997-09-22 1999-12-03 Inst Francais Du Petrole Methode statistique de classement d'evenements lies au proprietes physiques d'un milieu complexe tel que le sous-sol
GB2357097A (en) * 1999-12-08 2001-06-13 Norske Stats Oljeselskap Method of assessing positional uncertainty in drilling a well
US6401838B1 (en) 2000-11-13 2002-06-11 Schlumberger Technology Corporation Method for detecting stuck pipe or poor hole cleaning
GB2404988B (en) * 2002-04-10 2006-04-12 Schlumberger Technology Corp Method,apparatus and system for pore pressure prediction in presence of dipping formations
GB2385923B (en) * 2002-05-24 2004-07-28 Statoil Asa System and method for electromagnetic wavefield resolution
DE10317065A1 (de) * 2002-12-16 2004-07-22 Koenig & Bauer Ag Verfahren und Vorrichtung zur Steuerung und Verfahren zum Konfigurieren einer Anlage
GB2399640B (en) 2003-03-17 2007-02-21 Statoil Asa Method and apparatus for determining the nature of submarine reservoirs
GB2409900B (en) 2004-01-09 2006-05-24 Statoil Asa Processing seismic data representing a physical system
US20060100836A1 (en) * 2004-11-09 2006-05-11 Amardeep Singh Performance forecasting and bit selection tool for drill bits
US20060166224A1 (en) * 2005-01-24 2006-07-27 Norviel Vernon A Associations using genotypes and phenotypes
KR100735555B1 (ko) * 2005-09-15 2007-07-04 삼성전자주식회사 동작에 따른 기능을 수행하는 장치 및 방법
GB2434868B (en) 2006-02-06 2010-05-12 Statoil Asa Method of conducting a seismic survey
GB2435693A (en) 2006-02-09 2007-09-05 Electromagnetic Geoservices As Seabed electromagnetic surveying
GB2439378B (en) 2006-06-09 2011-03-16 Electromagnetic Geoservices As Instrument for measuring electromagnetic signals
GB2442749B (en) 2006-10-12 2010-05-19 Electromagnetic Geoservices As Positioning system
GB2445582A (en) 2007-01-09 2008-07-16 Statoil Asa Method for analysing data from an electromagnetic survey
CN104024572A (zh) 2011-11-02 2014-09-03 兰德马克绘图国际公司 用于预测钻柱卡钻事件的方法和***
AU2013258157B2 (en) * 2012-05-08 2016-02-25 Shell Internationale Research Maatschappij B.V. Method and system for sealing an annulus enclosing a tubular element
RU2015123680A (ru) * 2013-01-03 2017-02-08 Лэндмарк Графикс Корпорейшн Система и способ прогнозирования и визуализации событий в процессе бурения
MX2015009661A (es) * 2013-02-27 2016-03-21 Landmark Graphics Corp Metodo y sistema para predecir eventos de perforacion.
GB201406131D0 (en) * 2014-04-04 2014-05-21 Epidote Holdings Ltd System and method for determining deformed pipe geometry
US11598195B2 (en) * 2014-10-27 2023-03-07 Baker Hughes, A Ge Company, Llc Statistical approach to incorporate uncertainties of parameters in simulation results and stability analysis for earth drilling
CN104632076B (zh) * 2014-12-22 2017-04-05 中国石油天然气股份有限公司 一种丛式井组的钻井方法
CN105842152B (zh) * 2015-01-15 2018-11-16 中国石油天然气股份有限公司 泥饼力学性质测量仪
US10513920B2 (en) 2015-06-19 2019-12-24 Weatherford Technology Holdings, Llc Real-time stuck pipe warning system for downhole operations
CN105350932B (zh) * 2015-11-03 2017-10-03 辽河石油勘探局 一种气井带压诱喷解堵排液工艺
US10947811B2 (en) 2017-12-01 2021-03-16 Saudi Arabian Oil Company Systems and methods for pipe concentricity, zonal isolation, and stuck pipe prevention
US10557326B2 (en) 2017-12-01 2020-02-11 Saudi Arabian Oil Company Systems and methods for stuck pipe mitigation
US10557317B2 (en) 2017-12-01 2020-02-11 Saudi Arabian Oil Company Systems and methods for pipe concentricity, zonal isolation, and stuck pipe prevention
US10612360B2 (en) 2017-12-01 2020-04-07 Saudi Arabian Oil Company Ring assembly for measurement while drilling, logging while drilling and well intervention
CN108343424B (zh) * 2017-12-19 2021-08-03 中国石油天然气股份有限公司 钻井位置的确定方法和装置
US11286766B2 (en) 2017-12-23 2022-03-29 Noetic Technologies Inc. System and method for optimizing tubular running operations using real-time measurements and modelling
CN108894744B (zh) * 2018-08-29 2024-05-31 中煤特殊凿井有限责任公司 一种破除孤石卡钻的装置
NO20231246A1 (en) * 2021-06-29 2023-11-14 Landmark Graphics Corp Calculating pull for a stuck drill string
US20230121791A1 (en) * 2021-10-18 2023-04-20 Saudi Arabian Oil Company Pre-emptive jarring apparatus and methods of use thereof
CN114482856B (zh) * 2021-12-22 2023-02-28 中煤科工集团西安研究院有限公司 复杂破碎地层近水平定向钻进卡钻处理钻具组合及方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3152642A (en) * 1961-01-30 1964-10-13 Jr Albert G Bodine Acoustic method and apparatus for loosening and/or longitudinally moving stuck objects
US3552502A (en) * 1967-12-21 1971-01-05 Dresser Ind Apparatus for automatically controlling the killing of oil and gas wells
NL7209281A (zh) * 1971-09-15 1973-03-19
US3811509A (en) * 1972-08-31 1974-05-21 Cities Service Oil Co Well drilling method
US4270620A (en) * 1979-01-12 1981-06-02 Dailey Oil Tools, Inc. Constant bottom contact tool
US4428441A (en) * 1979-04-04 1984-01-31 Mobil Oil Corporation Method and apparatus for reducing the differential pressure sticking tendency of a drill string
US4384625A (en) * 1980-11-28 1983-05-24 Mobil Oil Corporation Reduction of the frictional coefficient in a borehole by the use of vibration
US4427080A (en) * 1980-12-11 1984-01-24 Exxon Production Research Co. Apparatus for preventing differential sticking in wells
US4423791A (en) * 1982-01-21 1984-01-03 Moses Jr Edward L Method of inhibiting differential wall sticking in the rotary drilling of hydrocarbon wells
AU1632283A (en) * 1983-06-28 1985-01-03 Moses, E.L. Inhibiting sticking
GB8411361D0 (en) * 1984-05-03 1984-06-06 Schlumberger Cambridge Researc Assessment of drilling conditions
AU608503B2 (en) * 1985-07-15 1991-04-11 Chevron Research And Technology Company Method of avoiding stuck drilling equipment

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101353959B (zh) * 2008-09-10 2012-08-29 杜书东 防溜钻、遇阻卡自动报警保护装置
CN102979501A (zh) * 2010-04-12 2013-03-20 国际壳牌研究有限公司 自动地提升钻头离开地下地层中的开孔底面的方法
CN102979501B (zh) * 2010-04-12 2015-11-18 国际壳牌研究有限公司 自动地提升钻头离开地下地层中的开孔底面的方法
US10415365B2 (en) 2010-04-12 2019-09-17 Shell Oil Company Methods and systems for drilling
CN106150476A (zh) * 2015-04-09 2016-11-23 中国石油化工股份有限公司 一种预测钻柱的粘吸卡钻风险的***
CN106150476B (zh) * 2015-04-09 2019-04-30 中国石油化工股份有限公司 一种预测钻柱的粘吸卡钻风险的***

Also Published As

Publication number Publication date
DE3688571T2 (de) 1993-10-07
EP0209343B1 (en) 1993-06-16
EP0209343A2 (en) 1987-01-21
DE209343T1 (de) 1990-04-12
CA1257701A (en) 1989-07-18
NO862850D0 (no) 1986-07-14
ES2000508A6 (es) 1988-03-01
DE3688571D1 (de) 1993-07-22
AU5944586A (en) 1987-01-22
AU608503B2 (en) 1991-04-11
US4791998A (en) 1988-12-20
DK334286D0 (da) 1986-07-14
NO862850L (no) 1987-01-16
CN1011429B (zh) 1991-01-30
DK334286A (da) 1987-01-16
EP0209343A3 (en) 1989-03-22

Similar Documents

Publication Publication Date Title
CN86104849A (zh) 避免钻井设备卡在井中的方法
CN1040653A (zh) 钻机避卡装置反压差卡钻解卡概率的确定方法
CN1081721C (zh) 测定岩石抗压强度的方法
CN1966934A (zh) 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法
CN1038327A (zh) 从油藏中采油
CN1300439C (zh) 用于确定到达定向目标的钻探路线的方法和装置
CN1791808A (zh) 用于钻井和测井的重力技术
CN101038680A (zh) 基于三维建模的立方体预测模型找矿方法
CN1628207A (zh) 钻井***
CN101253402A (zh) 评价欠饱和煤层气储层的方法
CN1573013A (zh) 定向电磁波电阻率装置和方法
Oloruntobi et al. Application of specific energy for lithology identification
CN1206546C (zh) 检测地壳中最大水平应力方向和相对大小的方法
Land et al. Louann salt geochemistry (Gulf of Mexico sedimentary basin, USA): a preliminary synthesis
CN1050088A (zh) 为绘制地下岩石夹杂物组成图而获取聚集流体夹杂挥发物
CN1043763A (zh) 基桩的性能与质量的判定和设计方法,地基性能的计测装置,以及就地灌注桩等基桩的掘孔方法及其装置
CN1302386C (zh) 低浓度表面活性剂与相态结合的三元复合驱计算机仿真方法
CN1896458A (zh) 导出井下流体的流体性质及其不确定度的***和方法
CN1271310C (zh) 确定钻井液密度的方法及控制钻井液密度的设备
Azi Casing Setting Depth using Bottom-Up Method for Development Well in the Offshore
CN1194233C (zh) 一种采用单井化学示踪剂测剩余油饱和度的方法
CN1019520B (zh) 估计一个储油层中至少两层地层的渗透性及表层系数的方法
Browne Mineralogical guides to interpreting the shallow paleohydrology of epithermal mineral depositing environments
RU2540701C2 (ru) Способ вскрытия продуктивного пласта на депрессии
Coker Managed pressure drilling applications index

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C13 Decision
GR02 Examined patent application
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee