CN105718699A - 大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法 - Google Patents
大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法 Download PDFInfo
- Publication number
- CN105718699A CN105718699A CN201610105728.XA CN201610105728A CN105718699A CN 105718699 A CN105718699 A CN 105718699A CN 201610105728 A CN201610105728 A CN 201610105728A CN 105718699 A CN105718699 A CN 105718699A
- Authority
- CN
- China
- Prior art keywords
- annular
- bed
- detritus
- height
- bottom hole
- 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.)
- Granted
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 55
- 238000005520 cutting process Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 29
- 230000008859 change Effects 0.000 title claims abstract description 27
- 238000013508 migration Methods 0.000 claims abstract description 46
- 230000005012 migration Effects 0.000 claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000013461 design Methods 0.000 claims abstract description 7
- 239000011435 rock Substances 0.000 claims description 65
- 230000008569 process Effects 0.000 claims description 14
- 239000002245 particle Substances 0.000 claims description 10
- 239000002002 slurry Substances 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 239000000725 suspension Substances 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 2
- 230000035515 penetration Effects 0.000 claims description 2
- 238000004062 sedimentation Methods 0.000 claims 2
- 239000012530 fluid Substances 0.000 description 6
- 230000007547 defect Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Earth Drilling (AREA)
Abstract
本发明提供一种大位移钻井环空岩屑床高度分布和井底的压力实时变化计算方法,包括根据钻井基本参数,获取环空稳态岩屑床高度;获取环空岩屑床动态运移速度;根据环空岩屑床动态运移速度,获取环空岩屑床动态运移距离;根据环空岩屑床动态运移距离,获取环空岩屑床高度分布;根据钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗;根据不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力。本发明能快速、准确的计算钻进时环空岩屑床高度和井底压力的实时分布变化,为长斜段钻进水力参数设计及计算提供帮助,有助于减少井下事故、提高钻井效率。
Description
技术领域
本发明涉及石油与天然气开技术领域,更为具体地,涉及一种大位移钻井大斜度井段和水平井井段(以下以长斜段统称)钻井时环空岩屑床高度分布和井底压力实时变化的简便计算方法。
背景技术
大位移井本质上是水平位移较长的定向井,目前,其水力学计算现场多是直接使用定向井研究成果。这在大斜度井段和水平井段不太长时,误差可以接受。但在应用于高水垂比大位移井,尤其是海上超高水垂比大位移井时,将会产生较大误差。其原因如下:一方面,海上很多大位移井其大斜度井段和水平井段长达3000-5000米,在此井段钻井时,地层破裂压力基本不变,而井底压力则随着环空压耗的增加而逐渐增加,两者差值越来越小,极易诱发井下事故。这就要求对于井底压力(环空压耗)的计算和预测必须足够准确和精确。另一方面,研究表明,大斜度井段和水平井段岩屑易下沉形成岩屑床,岩屑床对环空压耗的影响很大,要精确计算环空压耗就必须首先了解环空岩屑床分布。随着大位移井井深的增加,大斜度井段越来越长,钻井时岩屑床从形成到铺满整个环空需要很长时间,再加上一些施工操作的影响,在此过程中,环空岩屑床分布是实时变化的,稳态模型仅能计算环空达到稳定后最终形成的岩屑床高度,而无法计算整个过程中岩屑床的分布变化。岩屑床分布未知,环空压降自然也无法精确计算,无法满足大位移井井眼清洁和水力参数设计的需要。目前,国内外已建立了多个岩屑床动态运移模型,这些模型均是复杂的非线性方程组,求解难度高、速度慢,不适合工程应用,很难为工程设计提供直接帮助。
发明内容
鉴于上述问题,本发明的目的是提供一种大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,以解决现有复杂非线性岩屑床动态移动模型的缺陷的问题,此方法快速、准确的计算钻井时环空岩屑床高度和井底压力的实时分布变化,为长斜段钻井水力参数设计及计算提供帮助,减少井下事故、提高钻井效率。
本发明提供一种大位移钻井环空岩屑床高度分布和井底的压力实时变化计算方法,包括:
根据获取的钻井基本参数,获取环空稳态岩屑床高度;
根据所述环空稳态岩屑床高度,获取岩屑床推移速度和悬浮岩屑沉降成床速度;
根据所述岩屑床推移速度和悬浮岩屑沉降成床速度,获取环空岩屑床动态运移速度;
根据所述环空岩屑床动态运移速度,获取环空岩屑床动态运移距离;
根据所述环空稳态岩屑床高度和所述环空岩屑床动态运移距离,获取环空岩屑床高度分布;
根据所述钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗;
根据所述不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力。
此外,优选的方案是,根据钻井设计和实测数据,获取钻井基本参数;其中,所述钻井基本参数包括段井眼内径、钻柱外径、泥浆排量、泥浆密度、泥浆流变参数、钻柱转速、钻柱偏心度、机械钻速、井斜角、岩屑密度、岩屑粒径和井底的垂直深度。
此外,优选的方案是,在根据获取的钻井基本参数,获取环空稳态岩屑床高度的过程中,采用如下公式获取:
其中,Hc为岩屑床高度;
dw、dpo和db分别为井眼内径、钻杆外径和钻头直径;
ρf和ρs分别为泥浆和岩屑密度;
μe为有效粘度;Ua为环空返速;ROP为机械钻速;
RPM为钻杆转速;θ为井斜角;ds为岩屑直径;e为偏心度。
此外,优选的方案是,在根据所述环空稳态岩屑床高度,获取岩屑床推移速度和悬浮岩屑沉降成床速度的过程中,其中,所述岩屑床推移速度计算公式为:
所述悬浮岩屑沉降成床速度计算公式为:
其中,Ucm为岩屑床推移速度;
Ucd为悬浮岩屑沉降成床速度;
fhd为悬浮层与移动层界面摩擦系数;
Cd为移动岩屑浓度;
φ为颗粒内摩擦角,tanφ相当于颗粒动摩擦系数,当Cd=0.3~0.52时,tanφ=0.7~0.8;
Qs为入口岩屑体积流量;
Qf为入口泥浆体积流量;
k1、k2为修正系数;k3为比例系数。
此外,优选的方案是,在根据所述岩屑床推移速度和悬浮岩屑沉降成床速度,获取环空岩屑床动态运移速度的过程中,其中,所述环空岩屑床动态运移速度的计算公式为:
其中,Uc为环空岩屑床动态运移速度;
k3为比例系数;
Ucm为岩屑床推移速度;
Ucd为悬浮岩屑沉降成床速度。
此外,优选的方案是,在根据所述环空岩屑床动态运移速度,获取环空岩屑床动态运移距离的过程中,其中获取所述环空岩屑床动态运移距离的公式如下:
Lc=tUc
其中,Uc为环空岩屑床动态运移速度;
t为岩屑床开始形成所经历的时间;
Lc为经过t时间后,稳态岩屑床运移的距离。
此外,优选的方案是,在根据所述环空稳态岩屑床高度和所述环空岩屑床动态运移距离,获取环空岩屑床高度分布的过程中,其中,
若x≥MD-Lc,则Hcx=Hc
若x<MD-Lc,则Hcx公式如下:
其中,Hc为岩屑床高度;
MD为井总测深;
x为欲求位置的测深,在井底为MD,在井口为0;
Hcx为x位置处岩屑床高度。
此外,优选的方案是,在根据所述钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗的过程中,获取所述环空压耗的公式如下:
其中,ΔPc为有岩屑床存在时的环空压耗;
ΔPa为无岩屑床存在时的环空压耗;
fb为达西-韦斯巴赫系数,层流时紊流时
Hc为岩屑床高度;
dw、dpo分别为井眼内径、钻杆外径;
ρf和ρs分别为泥浆和岩屑密度;
Qf为入口泥浆体积流量。
此外,优选的方案是,在根据所述不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力的过程中,获取所述井底压力的公式如下:
Pw=ρfgh+ΔPc
其中,Pw为井底压力;
ρf为泥浆密度;
h为井底的垂直深度;
ΔPc为有岩屑床存在时的环空压耗。
从上面的技术方案可知,本发明提供的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,能够取得以下有益效果:
(1)该方法考虑岩屑床的形成过程,能够实时计算不同时刻环空岩屑床高度分布,以及对应的井底压力,相比稳态模型更加贴近实际情况,计算误差更小。
(2)通过计算岩屑动态运移速度,并利用几何曲线描述环空岩屑床高度分布,将稳态岩屑床高度计算模型与岩屑床高度分布动态变化结合起来,扩展了现有稳态岩屑运行模型的应用范围。
(3)本发明的方法相比于现有的岩屑动态运移模型(复杂非线性方程组),计算简便、速度快,非常适合工程应用,填补了复杂岩屑动态运移模型无法直接应用于工程实际的缺陷,为工程设计提供一种有效手段。
为了实现上述以及相关目的,本发明的一个或多个方面包括后面将详细说明并在权利要求中特别指出的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。然而,这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外,本发明旨在包括所有这些方面以及它们的等同物。
附图说明
通过参考以下结合附图的说明及权利要求书的内容,并且随着对本发明的更全面理解,本发明的其它目的及结果将更加明白及易于理解。在附图中:
图1为根据本发明实施例的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法流程示意图。
在所有附图中相同的标号指示相似或相应的特征或功能。
具体实施方式
在下面的描述中,出于说明的目的,为了提供对一个或多个实施例的全面理解,阐述了许多具体细节。然而,很明显,也可以在没有这些具体细节的情况下实现这些实施例。
以下将结合附图对本发明的具体实施例进行详细描述。
为了说明本发明提供的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,图1示出了根据本发明实施例的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法流程。
如图1所示,本发明提供的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法包括:
S110:根据获取的钻井基本参数,获取环空稳态岩屑床高度;
S120:根据环空稳态岩屑床高度,获取岩屑床推移速度和悬浮岩屑沉降成床速度;
S130:根据岩屑床推移速度和悬浮岩屑沉降成床速度,获取环空岩屑床动态运移速度;
S140:根据环空岩屑床动态运移速度,获取环空岩屑床动态运移距离;
S150:根据环空稳态岩屑床高度和环空岩屑床动态运移距离,获取环空岩屑床高度分布;
S160:根据钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗;
S170:根据不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力。
上述为大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法的具体流程,在获取环空稳态岩屑床高度之前,首先要获取钻井基本参数,根据钻井设计和实测数据,获取钻井基本参数,钻井基本参数包括计算段井眼内径、钻柱外径、泥浆排量、泥浆密度、泥浆流变参数、钻柱转速、钻柱偏心度、机械钻速、井斜角、岩屑密度、岩屑粒径、井底的垂直深度。
在上述步骤S110中,环空稳态岩屑床高度可以根据以下公式(1)计算,也可根据其它稳态岩屑床高度模型计算。
其中,Hc为岩屑床高度,单位为mm;dw、dpo和db分别为井眼内径、钻杆外径和钻头直径,单位为mm;ρf和ρs分别为泥浆和岩屑密度,单位为g/cm3;μe为有效粘度,单位为cP;Ua为环空返速,单位为m/s;ROP为机械钻速,单位m/h;RPM为钻杆转速,单位为转/分;θ为井斜角,单位为°;ds为岩屑直径,单位为mm;e为偏心度,无因次。
对于幂律流体,有效粘度μe可以根据公式(2)计算:
其中,n为幂律流体流性指数,无因次;k为幂律流体稠度系数,Pa﹒sn。
在上述步骤S120和S130中,计算环空岩屑床动态运移速度的具体步骤为:计算悬浮层与岩屑床层界面摩擦系数;计算岩屑床推移速度;计算悬浮岩屑沉降成床速度;岩屑床推移速度和悬浮岩屑沉降成床速度按比例相加即为岩屑床动态运移速度;具体如下:
环空岩屑床动态运移速度计算可以根据以下公式(3)计算:
其中,岩屑床推移速度Ucm计算公式为:
悬浮岩屑沉降成床速度Ucd计算公式为:
其中,fhd为悬浮层与移动层界面摩擦系数,计算方法如下公式(7),也可采用其它界面间摩擦系数;Cd为移动岩屑浓度;φ为颗粒内摩擦角,tanφ相当于颗粒动摩擦系数,与颗粒浓度有关,浓度越大,值越大,Cd=0.3~0.52时,tanφ=0.7~0.8;Qs为入口岩屑体积流量,m3/s;Qf为入口泥浆体积流量,m3/s;k1、k2为修正系数,与环空几何形态、公式简化误差等有关,可根据实验或实测数据求得,可取1;k3为比例系数,与泥浆悬浮携岩强度有关,可根据实验或实测数据求得,可取0.5。
其中,ds为岩屑直径;dw为井眼内径。
在步骤S140中,根据岩屑床运移时间和速度,计算岩屑床动态运移距离。其中,环空岩屑床动态运移距离可以根据以下公式(8)计算:
Lc=tUc(8)
其中,t为岩屑床开始形成或冲蚀到计算时刻所经历的时间,单位为s;Lc为经过t时间后,稳态岩屑床运移的距离,单位为m。
在步骤S150中,计算不同时刻环空岩屑床高度分布的具体步骤为:从井底向上,将环空分成两段,一段为稳态高度的岩屑床段,一段为从稳态高度逐渐递减至零的岩屑床段,前者环空岩屑床高度等于稳态岩屑床高度,后者岩屑床高度分布可以用一几何曲线描述,并用相应的公式计算。
具体地,不同时刻环空岩屑床高度分布可以根据以下公式(9)计算:
若x≥MD-Lc,则Hcx=Hc;
若x<MD-Lc,则Hcx按式(9)计算。
其中,MD为井总测深,单位为m;x为欲求位置的测深,在井底为MD,在井口为0,单位为m;Hcx为x位置处岩屑床高度,单位为m。
在步骤S160和S170中,计算不同时刻环空压耗和井底压力的具体步骤为:对井筒环空进行分段,使用以上方法计算某一时刻环空岩屑床高度分布,即每一段环空岩屑床高度,再根据公式计算每一段有/无岩屑床环空压耗,所有段相加即为当前时刻总环空压耗;根据泥浆密度和垂深,计算井底静压力;两者相加计算井底压力。
具体地,对井筒环空进行分段(段长可取25m或30m),使用以上公式(9)计算某一时刻环空岩屑床高度分布,即每一段环空岩屑床高度,再根据以下公式(10)计算每一段环空压耗,所有段相加即为当前时刻总环空压耗:
其中,ΔPc为有岩屑床存在时的环空压耗,Pa;ΔPa为无岩屑床存在时的环空压耗,Pa;fb为达西-韦斯巴赫系数,层流时紊流时
无岩屑床存在时的环空压耗ΔPa可以根据以下公式(11)计算:
其中,l为井段长度,m;f为模组系数,对于幂律流体,可以根据以下公式(12)或(13)计算:
层流摩阻系数:环空:f=24/Re(12)
紊流摩阻系数:
其中,Re为环空流体流动雷诺数,对于幂律流体,可以根据以下公式(14)计算:
不同时刻井底压力可以根据以下公式(15)计算:
Pw=ρfgh+ΔPc(15)
其中,Pw为井底压力;ρf为泥浆密度;h为井底的垂直深度;ΔPc为有岩屑床存在时的环空压耗。
通过上述实施方式可以看出,本发明提供的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,能够解决现有复杂非线性岩屑床动态运移模型的缺陷,该方法基于对岩屑运移机理和早期建立的岩屑动态运移模型的分析和简化,通过计算岩屑动态运移速度,并利用几何曲线描述环空岩屑床高度分布,将稳态岩屑床高度计算模型与岩屑床高度分布动态变化结合起来,能够快速准确计算环空岩屑床高度和井底压力实时分布变化。
如上参照附图以示例的方式描述了根据本发明提出的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法。但是,本领域技术人员应当理解,对于上述本发明所提出的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,还可以在不脱离本发明内容的基础上做出各种改进。因此,本发明的保护范围应当由所附的权利要求书的内容确定。
Claims (9)
1.一种大位移钻井环空岩屑床高度分布和井底的压力实时变化计算方法,包括:
根据获取的钻井基本参数,获取环空稳态岩屑床高度;
根据所述环空稳态岩屑床高度,获取岩屑床推移速度和悬浮岩屑沉降成床速度;
根据所述岩屑床推移速度和悬浮岩屑沉降成床速度,获取环空岩屑床动态运移速度;
根据所述环空岩屑床动态运移速度,获取环空岩屑床动态运移距离;
根据所述环空稳态岩屑床高度和所述环空岩屑床动态运移距离,获取环空岩屑床高度分布;
根据所述钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗;
根据所述不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力。
2.如权利要求1所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
根据钻井设计和实测数据,获取钻井基本参数;其中,
所述钻井基本参数包括段井眼内径、钻柱外径、泥浆排量、泥浆密度、泥浆流变参数、钻柱转速、钻柱偏心度、机械钻速、井斜角、岩屑密度、岩屑粒径和井底的垂直深度。
3.如权利要求1所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据获取的钻井基本参数,获取环空稳态岩屑床高度的过程中,采用如下公式获取:
其中,Hc为岩屑床高度;
dw、dpo和db分别为井眼内径、钻杆外径和钻头直径;
ρf和ρs分别为泥浆和岩屑密度;
μe为有效粘度;Ua为环空返速;ROP为机械钻速;
RPM为钻杆转速;θ为井斜角;ds为岩屑直径;e为偏心度。
4.如权利要求3所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述环空稳态岩屑床高度,获取岩屑床推移速度和悬浮岩屑沉降成床速度的过程中,其中,所述岩屑床推移速度计算公式为:
所述悬浮岩屑沉降成床速度计算公式为:
其中,Ucm为岩屑床推移速度;
Ucd为悬浮岩屑沉降成床速度;
fhd为悬浮层与移动层界面摩擦系数;
Cd为移动岩屑浓度;
φ为颗粒内摩擦角,tanφ相当于颗粒动摩擦系数,当Cd=0.3~0.52时,tanφ=0.7~0.8;
Qs为入口岩屑体积流量;
Qf为入口泥浆体积流量;
k1、k2为修正系数;k3为比例系数。
5.如权利要求4所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述岩屑床推移速度和悬浮岩屑沉降成床速度,获取环空岩屑床动态运移速度的过程中,其中,所述环空岩屑床动态运移速度的计算公式为:
其中,Uc为环空岩屑床动态运移速度;
k3为比例系数;
Ucm为岩屑床推移速度;
Ucd为悬浮岩屑沉降成床速度。
6.如权利要求5所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述环空岩屑床动态运移速度,获取环空岩屑床动态运移距离的过程中,其中获取所述环空岩屑床动态运移距离的公式如下:
Lc=tUc
其中,Uc为环空岩屑床动态运移速度;
t为岩屑床开始形成所经历的时间;
Lc为经过t时间后,稳态岩屑床运移的距离。
7.如权利要求6所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述环空稳态岩屑床高度和所述环空岩屑床动态运移距离,获取环空岩屑床高度分布的过程中,其中,
若x≥MD-Lc,则Hcx=Hc
若x<MD-Lc,则Hcx公式如下:
其中,Hc为岩屑床高度;
MD为井总测深;
x为欲求位置的测深,在井底为MD,在井口为0;
Hcx为x位置处岩屑床高度。
8.如权利要求1所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述钻井基本参数和所述岩屑床高度,获取不同时刻的环空压耗的过程中,获取所述环空压耗的公式如下:
其中,ΔPc为有岩屑床存在时的环空压耗;
ΔPa为无岩屑床存在时的环空压耗;
fb为达西-韦斯巴赫系数,层流时紊流时
Hc为岩屑床高度;
dw、dpo分别为井眼内径、钻杆外径;
ρf和ρs分别为泥浆和岩屑密度;
Qf为入口泥浆体积流量。
9.如权利要求8所述的大位移钻井环空岩屑床高度分布和井底压力实时变化的计算方法,其中,
在根据所述不同时刻的环空压耗和井底的垂直深度,获取不同时刻的井底压力的过程中,获取所述井底压力的公式如下:
Pw=ρfgh+ΔPc
其中,Pw为井底压力;
ρf为泥浆密度;
h为井底的垂直深度;
ΔPc为有岩屑床存在时的环空压耗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610105728.XA CN105718699B (zh) | 2016-02-25 | 2016-02-25 | 钻井环空岩屑床高度分布和井底压力实时变化的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610105728.XA CN105718699B (zh) | 2016-02-25 | 2016-02-25 | 钻井环空岩屑床高度分布和井底压力实时变化的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105718699A true CN105718699A (zh) | 2016-06-29 |
CN105718699B CN105718699B (zh) | 2019-03-12 |
Family
ID=56156044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610105728.XA Active CN105718699B (zh) | 2016-02-25 | 2016-02-25 | 钻井环空岩屑床高度分布和井底压力实时变化的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105718699B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598024A (zh) * | 2018-11-02 | 2019-04-09 | 长江大学 | 基于井眼清洁与环空压力耦合的钻井水力优化方法及设备 |
CN109854237A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 一种井底压力预测的方法及*** |
CN111236869A (zh) * | 2020-01-10 | 2020-06-05 | 中国石油大学(北京) | 停泵工况下岩屑分布的确定方法、装置和设备 |
CN112145167A (zh) * | 2020-11-10 | 2020-12-29 | 西南石油大学 | 一种基于施工数据的压裂动态井底压力计算方法 |
CN114382427A (zh) * | 2021-12-31 | 2022-04-22 | 中国石油天然气集团有限公司 | 水平井岩屑床处理方法及装置 |
CN114635654A (zh) * | 2021-12-22 | 2022-06-17 | 中国石油天然气集团有限公司 | 井眼清洁度确定方法、***、设备和计算机可读存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105178943A (zh) * | 2015-09-08 | 2015-12-23 | 中国石油天然气集团公司 | 一种实时校正井筒压力的方法 |
-
2016
- 2016-02-25 CN CN201610105728.XA patent/CN105718699B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105178943A (zh) * | 2015-09-08 | 2015-12-23 | 中国石油天然气集团公司 | 一种实时校正井筒压力的方法 |
Non-Patent Citations (6)
Title |
---|
TITUS N. OFEI 等: "Computational fluid dynamics simulation of the effect of drill pipe rotation on cuttings transport in horizontal wellbores using a Newtonian fluid", 《2015 INTERNATIONAL FIELD EXPLORATION AND DEVELOPMENT CONFERENCE》 * |
周英操 等: "欠平衡钻井环空多相流井底压力计算模型", 《石油学报》 * |
徐小龙: "大位移井井眼清洁评价及水力参数优选研究", 《中国优秀硕士学位全文数据库工程科技Ⅰ辑》 * |
苏堪华: "气侵过程中水平井井底压力变化规律研究", 《2011 INTERNATIONAL CONFERENCE ON ENVIRONMENTAL SYSTEMS SCIENCE AND ENGINEERING》 * |
郭晓乐 等: "大位移井岩屑运移规律研究", 《第二十届全国水动力学研讨会文集》 * |
郭晓乐 等: "大位移钻井岩屑动态运移模型及其应用研究", 《2011 INTERNATIONAL CONFERENCE ON ENERGY AND ENVIRONMENT》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598024A (zh) * | 2018-11-02 | 2019-04-09 | 长江大学 | 基于井眼清洁与环空压力耦合的钻井水力优化方法及设备 |
CN109854237A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 一种井底压力预测的方法及*** |
CN109854237B (zh) * | 2019-03-04 | 2019-10-22 | 西南石油大学 | 一种井底压力预测的方法及*** |
CN111236869A (zh) * | 2020-01-10 | 2020-06-05 | 中国石油大学(北京) | 停泵工况下岩屑分布的确定方法、装置和设备 |
CN111236869B (zh) * | 2020-01-10 | 2021-08-10 | 中国石油大学(北京) | 停泵工况下岩屑分布的确定方法、装置和设备 |
CN112145167A (zh) * | 2020-11-10 | 2020-12-29 | 西南石油大学 | 一种基于施工数据的压裂动态井底压力计算方法 |
CN114635654A (zh) * | 2021-12-22 | 2022-06-17 | 中国石油天然气集团有限公司 | 井眼清洁度确定方法、***、设备和计算机可读存储介质 |
CN114382427A (zh) * | 2021-12-31 | 2022-04-22 | 中国石油天然气集团有限公司 | 水平井岩屑床处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105718699B (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105718699B (zh) | 钻井环空岩屑床高度分布和井底压力实时变化的计算方法 | |
Busahmin et al. | Review on hole cleaning for horizontal wells | |
Cho et al. | A three-segment hydraulic model for cuttings transport in coiled tubing horizontal and deviated drilling | |
Ogunrinde et al. | Hydraulics optimization for efficient hole cleaning in deviated and horizontal wells | |
Maleki et al. | Comparing laminar and turbulent primary cementing flows | |
Pandya et al. | Wellbore cleanout in inclined and horizontal wellbores: the effects of flow rate, fluid rheology, and solids density | |
Huque et al. | Investigation of cuttings transport in a horizontal well with high-speed visualization and electrical resistance tomography technique | |
CN109854237B (zh) | 一种井底压力预测的方法及*** | |
Ytrehus et al. | Cuttings bed removal in deviated wells | |
CN113266344B (zh) | 一种水平井钻井液携岩效率预测方法 | |
Noah | Optimizing drilling fluid properties and flow rates for effective hole cleaning at high-angle and horizontal wells | |
CN109812236B (zh) | 一种确定异形井眼中的井眼清洁效果的方法 | |
Ytrehus et al. | Hydraulic behavior in cased and open-hole sections in highly deviated wellbores | |
Egbue | Improving the Efficiency of Transportation of Cuttings in Wellbore-Experimental Investigations on Critical Rolling and Lifting Velocities | |
Al-Shargabi et al. | Hole-cleaning performance in non-vertical wellbores: A review of influences, models, drilling fluid types, and real-time applications | |
Deng | Assessment of drilling fluid hole cleaning capacity in horizontal directional drilling–a parametric study of the effects of drilling fluid additives | |
Mahon | Development of an optimised integrated underbalanced drilling strategy for cuttings transport in gas-liquid flow through wellbore annuli. | |
Lund et al. | Displacement Efficiency In Eccentric Annuli | |
Abimbola et al. | Cuttings transport evaluation in deviated wells | |
Reda et al. | A statistical model for hole cleaning and cuttings transport efficiency during drilling of wells | |
CN109812237B (zh) | 一种满足异形井眼清洁的钻井液排量确定方法 | |
Moradi | A probabilistic study on hole cleaning optimization | |
Khaled | Innovative Tool for Optimizing Hole Cleaning Efficiency in Deviated and Extended Reach Wells | |
Cayeux et al. | Estimation of weight and torque on bit: assessment of uncertainties, correction and calibration methods | |
Rathgeber | Horizontal Flow Loop Design for the Study of Drill String Rotation Effects on Cuttings Transport |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
OL01 | Intention to license declared |