CN111677494A - 一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 - Google Patents
一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 Download PDFInfo
- Publication number
- CN111677494A CN111677494A CN202010081334.1A CN202010081334A CN111677494A CN 111677494 A CN111677494 A CN 111677494A CN 202010081334 A CN202010081334 A CN 202010081334A CN 111677494 A CN111677494 A CN 111677494A
- Authority
- CN
- China
- Prior art keywords
- logging
- well
- logging instrument
- cable
- borehole
- 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
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000003860 storage Methods 0.000 title claims description 9
- 239000012530 fluid Substances 0.000 claims abstract description 51
- 230000014509 gene expression Effects 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000005553 drilling Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 32
- 239000007788 liquid Substances 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 19
- 238000004458 analytical method Methods 0.000 claims description 17
- 230000005484 gravity Effects 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 230000005489 elastic deformation Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000007667 floating Methods 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/09—Locating or determining the position of objects in boreholes or wells, e.g. the position of an extending arm; Identifying the free or blocked portions of pipes
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种三维曲井中测井仪器压差遇卡判断与计算方法,包括如下步骤:获取测井轨迹参数生成井眼轨迹参数;建立受力平衡方程;结合测井电缆在井眼的受力表达式与力矩表达式生成测井电缆受力模型;带入井眼轨迹参数得到测井电缆在测井轨迹任意位置的下放摩阻;构建测井仪器压差遇卡的计算模型;输出遇卡时对应的井眼参数。本发明综合测井电缆与测井仪器的物理特性及力学性能,钻井液的触变性、密度、粘度因素、地层压力参数,建立一种测井电缆在井眼中的受力模型,并根据测定的井眼轨迹参数进行三维曲井中测井电缆下放摩阻的测定,对测井过程中的测井遇阻和测井遇卡现象进行判定、分析、预测与评价,为测井的正常进行夯实基础。
Description
技术领域
本发明涉及测井领域,尤其涉及一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质。
背景技术
测井是根据地质和地球物理条件,合理地选用综合测井方法,可以详细研究钻孔地质剖面、探测有用矿产、详细提供计算储量所必需的数据,如油层的有效厚度、孔隙度、含油气饱和度和渗透率等,以及研究钻孔技术情况等任务的勘探方法。测井作为勘探与开发油气田的重要方法技术,至今已近80年的历史。随着科技进步和测井技术本身的发展,它在油气勘探、开发和生产的全过程中发挥着更大的作用,为油气工业带来更高的经济效益。近十几年来的测井技术,特别是20世纪90年代后,取得了重大进展。
在面对油田储层埋藏深的油井,由于地层压力系数高,井型与井身结构发生较大变化,与此同时,高密度钻井液固相含量高、黏切值高、触变性大,对测井带来了巨大的困难,造成测井事故复杂多,成功率低。再由于压差、钻井液性能和井眼质量等原因,高密度定向常规测井遇阻、遇卡问题尤为突出,电缆测井被迫改为钻杆传输测井,导致地质资料无法取全。综上所述,酒东地区定向井测井存在较为突出的测井遇阻、遇卡问题。超低的常规测井成功率,一方面严重影响了钻井提速,另一方面由于无法常规测井导致地质测井资料采集不足,影响了油田对油藏储量的评价和勘探开发的布局。同时,间接影响了钻井公司的声誉和市场。
为了对测井过程中的测井遇阻和测井遇卡现象进行分析、预测、评价,需要对测井过程中测井电缆的下放摩阻进行测定,现有的测井电缆下放摩阻的方法对测井仪器的影响因素以及受力分析不全面,会导致测定的下放摩阻准确性不高,影响测井过程的正常运行。
发明内容
本发明的目的在于,针对上述问题,提出一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质。
本发明的第一个方面,提供一种三维曲井中测井仪器压差遇卡判断与计算方法,包括如下步骤:
S1:获取测井过程中测井电缆在自然坐标系下的测井轨迹参数生成井眼轨迹参数;
S2:根据测井电缆在井中的载荷情况建立受力平衡方程,得到测井电缆在井眼的受力表达式与力矩表达式;
S3:计算因测井电缆的弹性形变引起的力矩变化表达式,并结合得到的受力表达式与力矩表达式生成测井电缆在井眼中的受力模型;
S4:将井眼轨迹参数带入测井电缆在井眼中的受力模型得到测井电缆在测井轨迹任意位置的下放摩阻;
S5:根据测井仪器在井中的压差遇卡受力表达式构建测井仪器压差遇卡的计算模型;
S6:利用构建的测井仪器压差遇卡计算模型判断测井仪器在井中下放是否遇卡,若遇卡则输出遇卡时对应的井眼参数。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的步骤S2包括如下子步骤:
S201:计算测井电缆自重、支反力、摩擦阻力、两端内力和流体的粘滞摩阻的载荷矢量,
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的自重载荷矢量的计算包括:
步骤A1:根据测井电缆的参数计算单位长度测井电缆浮重qm,其中qm=q·Kf,井段Kf为浮力系数;q为测井电缆在空气中的单位重量,kN/m;Kf=1-ρm/ρs;ρm为井筒内液密度,ρs为测井电缆材料密度;
步骤A2:根据获得的测井电缆浮重qm计算井眼内测井工具串自重载荷矢量其中其中α为井斜角,rad;为方位角, rad;kα为井斜变化率,rad/m;为方位变化率,rad/m;k为井眼曲率rad/m,分别表示测井电缆轨迹的切线、主法线和副法线方向单位矢量。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的流体的粘滞摩阻载荷矢量其中,τf为流体结构力,N/m;μ为流体粘度,Ns/m2;ω为测井仪器串转动角速度,rad/s;Dw为井眼直径,m;R为测井电缆外半径, m;v为流体速度,m/s。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的两端内力载荷矢量的计算包括:
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的力矩变化表达式为Mb=EIK其中E为杨氏弹性模量,kN/m3;I为测井电缆的惯性矩,m4。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的测井电缆在井眼中的受力模型为:
式中:τf为流体结构力,N/m;μ为流体粘度,Ns/m2;ω为测井仪器串转动角速度,rad/s; Dw为井眼直径,m;R为测井电缆外半径,m;v为流体速度,m/s。E为杨氏弹性模量,kN/m3; I为测井电缆的惯性矩,m4。Nn为主法线方向的正压力,kN;Nb为副法线方向的正压力,kN;μt为周向方向的摩擦系数;μα为轴线方向的摩擦系数。Kf为浮力系数;q为测井电缆在空气中的单位重量,kN/m;Kf=1-ρm/ρs(ρm为井筒内液密度,ρs为测井电缆材料密度)。α为井斜角,rad;为井斜方位角,rad;Kα为井斜变化率,rad/m;Kφ为方位变化率,rad/m; K为井眼曲率(或称全角变化率),rad/m。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的步骤S4包括如下子步骤:
S401:将生成的井眼轨迹参数离散化并利用牛顿差分公式进行处理;
S402:将处理后的参数带入测井电缆在井眼中的受力模型计算得到测井电缆各位置对应的摩阻值。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的步骤S5包括如下子步骤:
S501:分别对测井仪器在井中受到的液体浮力、重力、液体压强、地层压力以及液体粘滞力进行受力分析获得对应受力的表达式;其中:
测井仪器串受到垂直向上的液体浮力为:Ffu=ρgV;Ffu为测井仪器在钻井液中受到的浮力,N;ρ为钻井液密度,kg/m3;V为测井仪器串排开水的体积,m3;
测井仪器串受到垂直向下的重力为:G=mg;G为测井仪器所受重力,N;m为测井仪器串的质量,kg;g为重力加速度,9.8N/kg;
测井仪器串受到垂直向下作用在测井仪器串侧面的液体压强:Py=ρgh;Py为测井仪器串一侧所受的液体压强,Pa;h为遇卡井段垂深,m;
测井仪器串受到与测井仪器运动方向相反的液体粘滞力:Fn为测井仪器串所受粘滞阻力,N;r为测井仪器半径,m;τ为流体的剪切力,N/m2;δ为流体的粘度,Pa·s;v为测井仪器串的下入速度,m/s;R为裸眼井段的井眼半径,m;
S502:根据得到的测井仪器对应受力的表达式计算测井仪器下放摩阻表达式:
f为测井仪器下放所收到的摩阻,N;μ为泥饼与测井仪器串之间的动摩擦系数;α为井斜角,rad;A为测井工具串与井壁泥饼的接触面积。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的测井工具串与井壁泥饼的接触面积其中式中l 为测井工具串与泥饼的接触长度,m;d1为井眼半径,m;d2为测井工具串半径,m;h1为测井工具串与泥饼的封闭处的泥饼厚度,m;h2为压实后泥饼厚度,m。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述判断测井仪器在井中下放是否遇卡通过比较测井仪器下放摩阻f与测井仪器串所受重力与浮力之差的余弦值(mg-ρgV)cosα的大小,当f>(mg-ρgV)cosα时,测井仪器下放遇卡。
本发明的第二个方面,提供一种三维曲井中测井仪器压差遇卡判断与计算系统,包括:井眼轨迹参数生成单元:获取测井过程中测井电缆在自然坐标系下的测井轨迹参数生成井眼轨迹参数;
井眼受力分析单元:根据测井电缆在井中的载荷情况建立受力平衡方程,得到测井电缆在井眼的受力表达式与力矩表达式;
受力模型生成单元:计算因测井电缆的弹性形变引起的力矩变化表达式,并结合得到的受力表达式与力矩表达式生成测井电缆在井眼中的受力模型;
下放摩阻计算单元:将井眼轨迹参数带入测井电缆在井眼中的受力模型得到测井电缆在测井轨迹任意位置的下放摩阻。
本发明的第三个方面,提供一种三维曲井中测井仪器压差遇卡判断与计算的终端设备,所述终端设备包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述三维曲井中测井仪器压差遇卡判断与计算。
本发明的第四个方面,提供一种计算机可读存储介质,所述计算机可读存储介质上存储有指令,该指令被处理器执行时实现上述三维曲井中测井仪器压差遇卡判断与计算方法。
本发明的有益效果:本发明综合测井电缆与测井仪器的物理特性及力学性能,钻井液的触变性、密度、粘度因素、地层压力参数,建立一种测井电缆在井眼中的受力模型,并根据测定的井眼轨迹参数进行三维曲井中测井电缆下放摩阻的测定,对测井过程中的测井遇阻和测井遇卡现象进行判定、分析、预测与评价,为测井的正常进行夯实基础。
附图说明
图1是三维井眼轨迹的空间几何关系示意图。
图2是测井电缆微元段受力示意图。
图3是本发明的测井仪器压差遇卡的受力状况。
图4是本发明的原理流程示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图说明本发明的具体实施方式。
(一)井眼轨迹的几何描述和计算方法
钻井最终形成的空间轨迹称之为井眼轨迹,实指井眼轴线。实钻井眼轨迹都是三维的,是一条复杂的三维空间曲线。井眼轨迹的几何描述是后续建立测井仪器串及测井电缆轴向载荷模型的必要前提。因此在建立测井仪器串及测井电缆轴向载荷模型之前,有必要对井眼轨迹进行几何描述和插值计算。
1.井眼轨迹的描述方法:
一般认为井眼轨迹是一条连续光滑的空间曲线,常用空间直角坐标系Oxyz和自然坐标系OsTNB两种坐标系描述井眼轨迹,如图1所示。
用于描述井眼轨迹空间挠曲形态的参数主要有:
(1)井深(s)
井眼轨迹上任意一点到井口的长度称为井深,它是一条曲线的长度,所以也称之为斜深。对于实钻井眼轨迹来讲,测点处的井深称为测量深度,它通常以钻柱长度或测量电缆的长度来测量的。
(2)井斜角(α)
过井眼轴线上某测点作井眼轴线的切线,该切线向井眼前进方向延伸的部分称为井眼方向线。井斜角为井眼轨迹曲线上任意一点井眼方向线与铅垂线的夹角,如图2所示的α角。井斜角表示了井眼轨迹在该测点处倾斜的大小。
2.井眼轨迹的插值计算:
由测井得到的井眼轨迹坐标数据是不连续的,是一系列测深离散点及对应的井斜角和方位角。只用这些数据不能表示井眼轨迹的实际形态,也不便于力学分析计算。只有借助于插值计算方法,对参数进行曲线拟合,才能获得更多点来绘制连续光滑的井眼轨迹曲线。在众多的插值计算方法中,三次样条插值既有分段插值精度高的优点,又能使曲线保持光滑连续,其结果高度逼真。本文采用三次样条插值来计算井眼轨迹。设从测深s0开始至测深sN止,共测得N+1个点的井深、井斜角和方位角。
井深是井眼轨迹的一个基本参数,也是测点位置的标志。因此,在井眼轨迹计算中将井深作为自变量,而将井斜角和方位角表达为井深的函数。根据三次样条函数的定义和性质,可以构造出区间[sk-1,sk](k=1,2,3,......,N)上井斜样条函数α(S)和方位样条函数的表达式。
Mk=αn(sk),Mk-1=αn(sk-1);
K——测点序号;
Lk——测电长度,Lk=sk-sk-1,m;
s——插值点处的井深,m;
N——测点个数。
数Mk和mk既与井斜角和方位角的测量值有关,又与井口和井底的边界条件有关。设井口和井底的井斜角和方位角的二阶倒数为常数,则有
M0=MN=m0=mN=0 (2-4)
对于全井的N+1个测点,可以得到两组含有N-1个未知数Mk和mk(k=1,2,......,N-1)的线性方程组
上述方程组是对角线方程组,采用追赶法求解。将解Mk、mk(k=0,1,2,......,N)分别代入α(S)和的表达式中,即可求出[sk-1,sk](k=1,2,3,......,N)井段上任意井深处的井斜角和方位角,从而确定三维空间中的平滑井眼轨迹曲线。
3.微分几何曲线论:
对于图2所示的自然坐标系,原点Os位于空间井眼轨迹曲线上的任意一点,轴分别指向该曲线的切线方向、主法线方向和副法线方向。用表示切线方向单位矢量,表示主法线方向单位矢量,表示副法线方向单位矢量,两两垂直且符合右手定则。微分几何中将所构成的右手标架称为曲线在点Os的Frenet标架。
由微分几何中曲率、挠率的定义可得:
根据Frenet标架定义可得
将式(2-11)两边求导并代入式(2-8)、(2-9)和(2-10)、(2-12)得
联立式(2-8)、(2-9)、(3-13)即可以将Frenet标价的三个单位向量的导向量写成下列公式
这就是曲线论的基本公式,称为Frenet公式。Frenet公式在曲线论中起着极为重要的作用,它的特点是基本向量关于弧长的微商可用基本向量线性表示,而且表示系数由曲线的曲率和扰率完全确定。为方便使用,通常将其写成矩阵形式
Frenet公式是本文进行轴向载荷分析的基础。
(二)三维曲井中测井电缆力学模型建立
1.基本假设:
为了建立科学合理的测井电缆及测井仪器串轴向载荷力学模型,必须对井眼和测井工具串作适当简化。在认真分析测井仪器串轴向载荷研究中常用假设条件的基础上,根据本文建模特点,在分析中作如下基本假设:
(1)测井电缆与测井仪器串轴线与井眼轴线重合,即测井工具串单元的曲率与井眼曲率相同;
(2)井壁(套管内壁)是刚性的;
(3)测井电缆单元所受重力、正压力、摩阻力均匀分布;
(4)摩擦系数在同一口井或同一井段中为常数;
(5)暂不考虑动载荷对测井电缆和测井工具串的影响。
2.坐标系及几何关系:
由假设条件可知测井工具串轴线与井眼轴线具有相同轨迹。采用图2的空间直角坐标系 Oxyz和自然坐标系OsTNB两种坐标系来描述三维曲井中的测井仪器串。在空间直角坐标系中,用分别表示沿坐标轴x、y和z的单位矢量;在自然坐标系中,原点Os位于测井工具串轴线上的任意一点,分别表示测井仪器串的切线、主法线和副法线方向单位矢量。将测井工具串轴线上任一点矢径记为其中s为弧长。这样,测井工具串轴线上的任一点O(x,y,z)即可用来描述
两边对弧长s求导,得
由空间几何关系可得
将式(2-18)代入式(2-17)得
利用Frenet公式有
将式(3-19)代入上式得
根据Frenet标架定义可得
将式(2-19)、(2-21)代入(2-22),整理得到
这样,即可得到自然坐标系与直角坐标系的关系
式中:α——井斜角,rad;
kα——井斜变化率,rad/m;
k——井眼曲率(或称全角变化率)rad/m。
3.三维钢杆力学模型建立
通常,大斜度井的眼轨迹存在井眼轨迹曲率变化较大的井段。由于测井工具串受到三维弯曲井眼的约束作用,测井电缆刚度及井眼曲率会影响测井工具串在作业过程中的力学行为。因此,在受力分析中,在弯曲井段需要充分考虑测井电缆刚度和井眼曲率对测井工具串受力的影响。为此,本项目将建立考虑测井电缆刚度的三维“刚杆”力学模型用以描述大斜度井弯曲井段的测井电缆及测井仪器串受力情况。
首先,按照轨迹描述方法,在测井工具串轴线上取自然坐标系,并在测井电缆上取任一弧长为ds的微元体作为受力分析对象,如图2所示。以A点为始点,其曲线坐标为,B点为终点,其曲线坐标为,综合考虑各分部载荷、集中力、力矩的影响,对其进行受力分析。
(1)几何方程
根据微分几何的基本原理可得:
式中:
式中:α——井斜角,rad;
Kα——井斜变化率,rad/m;
Kφ——方位变化率,rad/m;
K——井眼曲率(或称全角变化率),rad/m。
(2)平衡方程
根据测井的工况特征可知,测井工具串在井眼内受到的载荷包括测井工具自重、测井电缆及测井仪器串与井壁的支反力、摩擦阻力、内外流体压力、流体的粘滞摩阻等。因此,根据受力工况,对测井电缆微元段进行受力分析。各种载荷的具体分析如下:
1)两端内力及内力矩
根据式(2-32),利用Frenet公式有
也即为:
也即为:
2)测井工具串自重
自重是测井电缆及测井仪器串在井眼中受到的一个对整体受力影响较大的均布载荷。在直井中,载荷即为测井电缆及测井仪器串自身重量,而在弯曲井眼中,则需要向各个方向进行分解分析。同时,由于井筒内充满了密度较高的钻井液,测井工具串将受到液体的浮力,因而在计算时采用浮重进行分析。单位长度测井电缆浮重qm为:
qm=q·Kf (2-39)
式中:Kf——浮力系数;
q——测井电缆在空气中的单位重量,kN/m;
Kf=1-ρm/ρs(ρm为井筒内液密度,ρs为测井电缆材料密度)。
因此,在弯曲井眼内测井工具串自重载荷矢量可以表示为:
3)支反力
在自重的作用下,测井工具串会受到井壁的支反力:
4)与井壁间的摩擦力
测井工具串在弯曲井眼内向下运动,因而会在测井工具串与井壁之间产生摩擦力:
式中:Nn——主法线方向的正压力,kN;
Nb——副法线方向的正压力,kN;
μt——周向方向的摩擦系数;
μα——轴线方向的摩擦系数。
5)流体粘滞摩阻
测井电缆受到流体的粘滞摩擦阻力可表示为:
式中:τf——流体结构力,N/m;
μ——流体粘度,Ns/m2;
ω——测井仪器串转动角速度,rad/s;
Dw——井眼直径,m;
R——测井电缆外半径,m;
v——流体速度,m/s。
经过测井工具串在测井作业过程中主要载荷的受力分析,根据受力平衡条件建立测井工具串微元段ds的受力平衡方程。
摩擦阻力产生的扭矩矢量为
粘滞阻力产生的扭矩矢量为
联立式(2-35)及式(2-47),可得
联立式(2-40)及式(2-43),可得
又由式(2-37)
(3)物理方程
在井眼轨迹中,测井电缆横向变形量受到井眼的约束,其变形仍在弹性变形范围内,因此测井电缆变形的物理关系为:
Mb=EIK(2-56)
式中:E——杨氏弹性模量,kN/m3;
I——测井电缆的惯性矩,m4。
将式(2-54)至(2-58)整理可得测井电缆在井眼中的受力模型:
由于微元段所在井眼轴线相邻两测点为空间斜平面上的一段圆弧,因此井眼挠率始终位于密切面内,由密切面定义可知τ=0,则(2-59)简化为:
式中:τf为流体结构力,N/m;μ为流体粘度,Ns/m2;ω为测井仪器串转动角速度,rad/s;Dw为井眼直径,m;R为测井电缆外半径,m;v为流体速度,m/s。E为杨氏弹性模量,kN/m3; I为测井电缆的惯性矩,m4。Nn为主法线方向的正压力,kN;Nb为副法线方向的正压力,kN;μt为周向方向的摩擦系数;μα为轴线方向的摩擦系数。Kf为浮力系数;q为测井电缆在空气中的单位重量,kN/m;Kf=1-ρm/ρs(ρm为井筒内液密度,ρs为测井电缆材料密度)。α为井斜角,rad;为井斜方位角,rad;Kα为井斜变化率,rad/m;Kφ为方位变化率,rad/m; K为井眼曲率(或称全角变化率),rad/m。
式(2-60)为综合考虑了井眼轨迹、测井电缆、井壁对测井电缆及测井仪器串摩擦力、井筒内液体对测井电缆和测井仪器串浮力以及液体粘滞阻力作用的测井工具串三维钢杆力学模型。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的步骤S4包括如下子步骤:引入牛顿差分公式,首先将微分方程离散化再利用牛顿差分公式,结合三维曲井中测井电缆及测井仪器串轴向载荷力学模型计算得到测井工具串在裸眼井段下放过程中各井深对应的摩阻。
(三)三维曲井中构建测井仪器压差遇卡的计算模型
一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的步骤S5包括如下子步骤: S501:如图4所示,分别对测井仪器在井中受到的液体浮力、重力、液体压强、地层压力以及液体粘滞力进行受力分析获得对应受力的表达式;其中:
测井仪器串受到垂直向上的液体浮力为:Ffu=ρgV;Ffu为测井仪器在钻井液中受到的浮力,N;ρ为钻井液密度,kg/m3;V为测井仪器串排开水的体积,m3;
测井仪器串受到垂直向下的重力为:G=mg;G为测井仪器所受重力,N;m为测井仪器串的质量,kg;g为重力加速度,9.8N/kg;
测井仪器串受到垂直向下作用在测井仪器串侧面的液体压强:Py=ρgh;Py为测井仪器串一侧所受的液体压强,Pa;h为遇卡井段垂深,m;
测井仪器串受到与测井仪器运动方向相反的液体粘滞力:Fn为测井仪器串所受粘滞阻力,N;r为测井仪器半径,m;τ为流体的剪切力,N/m2;δ为流体的粘度,Pa·s;v为测井仪器串的下入速度,m/s;R为裸眼井段的井眼半径,m;
S502:根据得到的测井仪器对应受力的表达式计算测井仪器下放摩阻表达式:
f为测井仪器下放所收到的摩阻,N;μ为泥饼与测井仪器串之间的动摩擦系数;α为井斜角,rad;A为测井工具串与井壁泥饼的接触面积。
进一步的,一种三维曲井中测井仪器压差遇卡判断与计算方法,所述的测井工具串与井壁泥饼的接触面积其中式中l 为测井工具串与泥饼的接触长度,m;d1为井眼半径,m;d2为测井工具串半径,m;h1为测井工具串与泥饼的封闭处的泥饼厚度,m;h2为压实后泥饼厚度,m。
一种三维曲井中测井仪器压差遇卡判断与计算方法,所述判断测井仪器在井中下放是否遇卡通过比较测井仪器下放摩阻f与测井仪器串所受重力与浮力之差的余弦值(mg-ρgV)cosα的大小,当f>(mg-ρgV)cosα时,测井仪器下放遇卡,计算该井段摩阻。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (17)
1.一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,包括如下步骤:
S1:获取测井过程中测井电缆在自然坐标系下的测井轨迹参数生成井眼轨迹参数;
S2:根据测井电缆在井中的载荷情况建立受力平衡方程,得到测井电缆在井眼的受力表达式与力矩表达式;
S3:计算因测井电缆的弹性形变引起的力矩变化表达式,并结合得到的受力表达式与力矩表达式生成测井电缆在井眼中的受力模型;
S4:将井眼轨迹参数带入测井电缆在井眼中的受力模型得到测井电缆在测井轨迹任意位置的下放摩阻;
S5:根据测井仪器在井中的压差遇卡受力表达式构建测井仪器压差遇卡的计算模型;
S6:利用构建的测井仪器压差遇卡计算模型判断测井仪器在井中下放是否遇卡,若遇卡则输出遇卡时对应的井眼参数。
9.根据权利要求1所述的一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,所述的力矩变化表达式为Mb=EIK,其中E为杨氏弹性模量,kN/m3;I为测井电缆的惯性矩,m4。
10.根据权利要求3所述的一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,所述的测井电缆在井眼中的受力模型为:
式中:τf为流体结构力,N/m;μ为流体粘度,Ns/m2;ω为测井仪器串转动角速度,rad/s;Dw为井眼直径,m;R为测井电缆外半径,m;v为流体速度,m/s。E为杨氏弹性模量,kN/m3;I为测井电缆的惯性矩,m4。Nn为主法线方向的正压力,kN;Nb为副法线方向的正压力,kN;μt为周向方向的摩擦系数;μα为轴线方向的摩擦系数。Kf为浮力系数;q为测井电缆在空气中的单位重量,kN/m;Kf=1-ρm/ρs(ρm为井筒内液密度,ρs为测井电缆材料密度)。α为井斜角,rad;为井斜方位角,rad;Kα为井斜变化率,rad/m;Kφ为方位变化率,rad/m;K为井眼曲率(或称全角变化率),rad/m。
11.根据权利要求3所述的一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,所述的步骤S4包括如下子步骤:
S401:将生成的井眼轨迹参数离散化并利用牛顿差分公式进行处理;
S402:将处理后的参数带入测井电缆在井眼中的受力模型计算得到测井电缆各位置对应的摩阻值。
12.根据权利要求1所述的一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,所述的步骤S5包括如下子步骤:
S501:分别对测井仪器在井中受到的液体浮力、重力、液体压强、地层压力以及液体粘滞力进行受力分析获得对应受力的表达式;其中:
测井仪器串受到垂直向上的液体浮力为:Ffu=ρgV;Ffu为测井仪器在钻井液中受到的浮力,N;ρ为钻井液密度,kg/m3;V为测井仪器串排开水的体积,m3;
测井仪器串受到垂直向下的重力为:G=mg;G为测井仪器所受重力,N;m为测井仪器串的质量,kg;g为重力加速度,9.8N/kg;
测井仪器串受到垂直向下作用在测井仪器串侧面的液体压强:Py=ρgh;Py为测井仪器串一侧所受的液体压强,Pa;h为遇卡井段垂深,m;
测井仪器串受到与测井仪器运动方向相反的液体粘滞力:Fn为测井仪器串所受粘滞阻力,N;r为测井仪器半径,m;τ为流体的剪切力,N/m2;δ为流体的粘度,Pa·s;v为测井仪器串的下入速度,m/s;R为裸眼井段的井眼半径,m;
S502:根据得到的测井仪器对应受力的表达式计算测井仪器下放摩阻表达式:
f为测井仪器下放所收到的摩阻,N;μ为泥饼与测井仪器串之间的动摩擦系数;α为井斜角,rad;A为测井工具串与井壁泥饼的接触面积。
14.根据权利要求13所述的一种三维曲井中测井仪器压差遇卡判断与计算方法,其特征在于,所述判断测井仪器在井中下放是否遇卡通过比较测井仪器下放摩阻f与测井仪器串所受重力与浮力之差的余弦值(mg-ρgV)cosα的大小,当f>(mg-ρgV)cosα时,测井仪器下放遇卡。
15.一种三维曲井中测井仪器压差遇卡判断与计算系统,其特征在于,包括:
井眼轨迹参数生成单元:获取测井过程中测井电缆在自然坐标系下的测井轨迹参数生成井眼轨迹参数;
井眼受力分析单元:根据测井电缆在井中的载荷情况建立受力平衡方程,得到测井电缆在井眼的受力表达式与力矩表达式;
受力模型生成单元:计算因测井电缆的弹性形变引起的力矩变化表达式,并结合得到的受力表达式与力矩表达式生成测井电缆在井眼中的受力模型;
下放摩阻计算单元:将井眼轨迹参数带入测井电缆在井眼中的受力模型得到测井电缆在测井轨迹任意位置的下放摩阻;
压差遇卡模型生成单元:根据测井仪器在井中的压差遇卡受力表达式构建测井仪器压差遇卡的计算模型;
压差遇卡判定单元:利用构建的测井仪器压差遇卡计算模型判断测井仪器在井中下放是否遇卡,若遇卡则输出遇卡时对应的井眼参数。
16.一种三维曲井中测井仪器压差遇卡判断与计算的终端设备,其特征在于,所述终端设备包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如权利要求1-11中任一所述的三维曲井中测井仪器压差遇卡判断与计算方法。
17.一种计算机可读存储介质,所述计算机可读存储介质上存储有指令,其特征在于,该指令被处理器执行时实现如权利要求1-11任意所述的三维曲井中测井仪器压差遇卡判断与计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010081334.1A CN111677494B (zh) | 2020-02-06 | 2020-02-06 | 一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010081334.1A CN111677494B (zh) | 2020-02-06 | 2020-02-06 | 一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111677494A true CN111677494A (zh) | 2020-09-18 |
CN111677494B CN111677494B (zh) | 2022-03-11 |
Family
ID=72451201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010081334.1A Active CN111677494B (zh) | 2020-02-06 | 2020-02-06 | 一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111677494B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113032987A (zh) * | 2021-03-11 | 2021-06-25 | 西南石油大学 | 一种无隔水管钻井气侵特性动态分析方法 |
CN113343544A (zh) * | 2021-07-21 | 2021-09-03 | 西南石油大学 | 一种连续油管下入水平井的力学分析方法 |
CN113391373A (zh) * | 2021-05-10 | 2021-09-14 | 三明学院 | 一种杨氏模量的计算方法、装置、设备和存储介质 |
CN115354996A (zh) * | 2022-10-21 | 2022-11-18 | 吉辰智能设备集团有限公司 | 遇阻指示仪及测井遇阻操作方法 |
CN115880482A (zh) * | 2023-02-17 | 2023-03-31 | 中海油田服务股份有限公司 | 一种测井图像遇卡识别与倾角校正方法、装置及计算设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1710249A (zh) * | 2005-08-05 | 2005-12-21 | 余丹 | 特殊井快速测井平台技术 |
US20100307742A1 (en) * | 2007-11-12 | 2010-12-09 | Phillips Wayne J | Method of determining and utilizing high fidelity wellbore trajectory |
CN204941415U (zh) * | 2015-08-05 | 2016-01-06 | 中国石油天然气股份有限公司 | 电缆测井用解卡装置 |
CN105550448A (zh) * | 2015-12-15 | 2016-05-04 | 中国石油天然气股份有限公司 | 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置 |
CN105863514A (zh) * | 2016-04-11 | 2016-08-17 | 西南石油大学 | 一种液控式随钻扩眼工具 |
CN105986803A (zh) * | 2015-07-15 | 2016-10-05 | 中国石油天然气集团公司 | 井下振动监测方法 |
CN106156385A (zh) * | 2015-04-09 | 2016-11-23 | 中国石油化工股份有限公司 | 一种预测钻柱的粘吸卡钻风险的方法 |
CN207420526U (zh) * | 2017-10-16 | 2018-05-29 | 西南石油大学 | 一种高压防硫化氢井口压力设备 |
CN110107224A (zh) * | 2019-05-05 | 2019-08-09 | 长江大学 | 一种水平井钻磨管柱下入摩阻计算方法 |
-
2020
- 2020-02-06 CN CN202010081334.1A patent/CN111677494B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1710249A (zh) * | 2005-08-05 | 2005-12-21 | 余丹 | 特殊井快速测井平台技术 |
US20100307742A1 (en) * | 2007-11-12 | 2010-12-09 | Phillips Wayne J | Method of determining and utilizing high fidelity wellbore trajectory |
CN106156385A (zh) * | 2015-04-09 | 2016-11-23 | 中国石油化工股份有限公司 | 一种预测钻柱的粘吸卡钻风险的方法 |
CN105986803A (zh) * | 2015-07-15 | 2016-10-05 | 中国石油天然气集团公司 | 井下振动监测方法 |
CN204941415U (zh) * | 2015-08-05 | 2016-01-06 | 中国石油天然气股份有限公司 | 电缆测井用解卡装置 |
CN105550448A (zh) * | 2015-12-15 | 2016-05-04 | 中国石油天然气股份有限公司 | 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置 |
CN105863514A (zh) * | 2016-04-11 | 2016-08-17 | 西南石油大学 | 一种液控式随钻扩眼工具 |
CN207420526U (zh) * | 2017-10-16 | 2018-05-29 | 西南石油大学 | 一种高压防硫化氢井口压力设备 |
CN110107224A (zh) * | 2019-05-05 | 2019-08-09 | 长江大学 | 一种水平井钻磨管柱下入摩阻计算方法 |
Non-Patent Citations (4)
Title |
---|
J.VENANCIO: "Problem detection during tripping operations horizontal and directional well", 《ONEPETRO》 * |
LIANGJIE MAO: "Dynamic mechanical behavior analysis of deep water drilling riser under hard hang-off evacuation conditions", 《OCEAN ENGINEERING 》 * |
吴应凯: "低压易漏地层防漏堵漏机理探讨及现场应用", 《天然气工业》 * |
赵勇: "冀中地区测井施工遇阻卡原因分析及建议", 《石油工业技术监督》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113032987A (zh) * | 2021-03-11 | 2021-06-25 | 西南石油大学 | 一种无隔水管钻井气侵特性动态分析方法 |
CN113391373A (zh) * | 2021-05-10 | 2021-09-14 | 三明学院 | 一种杨氏模量的计算方法、装置、设备和存储介质 |
CN113391373B (zh) * | 2021-05-10 | 2023-12-12 | 三明学院 | 一种杨氏模量的计算方法、装置、设备和存储介质 |
CN113343544A (zh) * | 2021-07-21 | 2021-09-03 | 西南石油大学 | 一种连续油管下入水平井的力学分析方法 |
CN113343544B (zh) * | 2021-07-21 | 2022-03-01 | 西南石油大学 | 一种连续油管下入水平井的力学分析方法 |
CN115354996A (zh) * | 2022-10-21 | 2022-11-18 | 吉辰智能设备集团有限公司 | 遇阻指示仪及测井遇阻操作方法 |
CN115354996B (zh) * | 2022-10-21 | 2023-01-13 | 吉辰智能设备集团有限公司 | 遇阻指示仪及测井遇阻操作方法 |
CN115880482A (zh) * | 2023-02-17 | 2023-03-31 | 中海油田服务股份有限公司 | 一种测井图像遇卡识别与倾角校正方法、装置及计算设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111677494B (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111677494B (zh) | 一种三维曲井中测井仪器压差遇卡判断与计算方法、系统、终端、存储介质 | |
RU2640324C2 (ru) | Калибровка моделирования бурения, включая оценку растяжения и скручивания бурильной колонны | |
US8594987B2 (en) | Systems and methods for modeling wellbore trajectories | |
US8490693B2 (en) | Determining fracture orientation using wellbore acoustic radial profiles | |
RU2663653C1 (ru) | Улучшенная оценка искривления ствола скважины, основанная на результатах измерений изгибающего момента инструмента | |
WO2020080973A1 (ru) | Способ и система комбинированного сопровождения процесса бурения скважины | |
Mitchell et al. | Drillstring analysis with a discrete torque/drag model | |
CN105911589B (zh) | 运用地震资料确定三维空间内断层滑动参数的定量化方法 | |
Liu et al. | Method and apparatus for monitoring of downhole dynamic drag and torque of drill-string in horizontal wells | |
BR112014013553B1 (pt) | Método e sistema para otimizar medições de peso em operações de perfuração, e, mídia legível por computador | |
Fazaelizadeh | Real time torque and drag analysis during directional drilling | |
Mirhaj et al. | Improvement of torque-and-drag modeling in long-reach wells | |
Cayeux et al. | Challenges and solutions to the correct interpretation of drilling friction tests | |
US11321506B2 (en) | Fast algorithm to simulate the response of PDC bits | |
CN104074503B (zh) | 一种超声测井仪系统的测井方法 | |
CN109854237B (zh) | 一种井底压力预测的方法及系统 | |
Mitchell et al. | Lateral Buckling—The Key to Lockup | |
EP3486427B1 (en) | Method of determining maximum stress in a well drilled in a reservoir | |
US5388044A (en) | Dipmeter processing technique | |
Heisig et al. | Continuous borehole curvature estimates while drilling based on downhole bending moment measurements | |
CN114991744A (zh) | 井下测量数据的时空转换方法及装置 | |
Zhang et al. | A new real-time hole cleaning monitoring method based on downhole multi-point pressure measurement and data driven approach | |
ElGizawy et al. | Is While-Drilling True Vertical Depth True? | |
Ifeanyi et al. | Predicting drillstring buckling | |
Sawaryn | Simple Engineering Applications Recycled as Effective Training Aids |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |