CN109869132A - 一种下套管摩阻系数计算方法 - Google Patents
一种下套管摩阻系数计算方法 Download PDFInfo
- Publication number
- CN109869132A CN109869132A CN201910056716.6A CN201910056716A CN109869132A CN 109869132 A CN109869132 A CN 109869132A CN 201910056716 A CN201910056716 A CN 201910056716A CN 109869132 A CN109869132 A CN 109869132A
- Authority
- CN
- China
- Prior art keywords
- casing
- coefficient
- friction resistance
- hook load
- centralizer
- 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
Landscapes
- Earth Drilling (AREA)
Abstract
本发明公开了一种下套管摩阻系数计算方法,属于石油钻采固井技术领域。本发明提供的方法包括:S1、在套管下至预定井深时,获取套管串中每根套管长度、井斜角、套管线重、钢材密度及钻井液密度;S2、实时监测记录套管下至预定井深过程的大钩载荷数据;S3、根据套管串中扶正器类型及所处位置,确定摩阻系数类型及个数,并选取与摩阻系数个数相同的时刻数,采集每个时刻的大钩载荷数据及对应的套管串数据;S4、基于每个时刻的大钩载荷数据及对应的套管串数据,列举大钩载荷计算公式得到方程组;S5、解方程组,得到摩阻系数值。本发明可以准确计算不同类型扶正器的摩阻系数,简单易行,能为下套管施工提供可靠的技术指导。
Description
技术领域
本发明涉及石油钻采固井技术领域,尤其涉及一种下套管摩阻系数计算方法。
背景技术
在非常规油气资源开发过程中,经常会遇到长水平段水平井。在这类水平井固井时,为保障套管的顺利下入,准确计算下套管的摩阻系数十分重要,为了提高套管居中度,经常会在套管串中安放不同类型的扶正器,不同类型扶正器与井壁及上层套管的摩阻系数还不相同。现有的摩阻系数计算方法大多是取光套管与井壁及上层套管的摩阻系数经验值,从而导致不同类型扶正器摩阻系数计算不准确,难以为下套管施工提供可靠的技术参考。
发明内容
本发明实施例提供了一种下套管摩阻系数计算方法,能够准确计算不同类型扶正器在井壁或上层套管中的摩阻系数。
结合本发明实施例的第一方面,提供了一种下套管摩阻系数计算方法,包括:
S1、在套管下至预定井深时,获取套管串中每根套管长度、每根套管所处位置的井斜角、套管线重、套管钢材密度及钻井液密度;
S2、实时监测记录套管下至预定井深过程的大钩载荷数据;
S3、根据套管串中扶正器类型及所处位置,确定摩阻系数类型及个数,并在下套管过程中选取与摩阻系数个数相同的时刻数,采集每个时刻的大钩载荷数据及对应的套管串数据;其中,所述套管串数据至少包括当前套管序号、套管长度、套管线重及套管所处位置的井斜角。
S4、基于所述每个时刻的大钩载荷数据及对应的套管串数据,列举大钩载荷计算公式(1),得到方程组,
其中,qc为套管线重,ρc为套管钢材密度,ρd为钻井液密度,Li为每根套管长度,i表示套管序号,θi为套管所在井深的井斜角,μk表示不同类型扶正器的摩阻系数,P1为大钩载荷值;
S5、解所述方程组,得到摩阻系数值。
从以上技术方案可以看出,本发明实施例具有以下优点:
本发明实施例中,通过记录套管下至一定井深时的套管串相关数据以及大钩载荷数据,根据扶正器类型及所处位置,确定摩阻系数类型及个数后,在下套管过程中选取与摩阻系数个数相同的时刻数,采集每个时刻的大钩载荷数据及对应的套管串数据,列举大钩载荷计算公式形成方程组并求解,即可获得不同类型扶正器的摩阻系数值。通过本发明实施例提供的技术方案可以准确计算不同类型扶正器的摩阻系数,从而优化扶正器类型、数量及位置设计,为下套管施工提供可靠的技术指导及参考。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的摩阻系数计算方法一个实施例流程图。
具体实施方式
本发明实施例提供了一种下套管摩阻系数计算方法,用于准确计算下套管过程中不同类型扶正器对应的摩阻系数。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1,本发明实施例中下套管摩阻系数计算方法一个实施例包括:
S101、在套管下至预定井深时,获取套管串中每根套管长度、每根套管所处位置的井斜角、套管线重、套管钢材密度及钻井液密度;
示例性的,在套管下至预定井深L时刻,获取下入套管串中每根套管长度Li,其中,i为套管段数,套管段数可上到下依次编号为1,2…n,记录每根套管所在井深的井斜角θi,其中,i等于套管段数,井斜角可从上到下依次编号为1,2…n,并获取套管线重qc,套管钢材密度ρc,钻井液密度ρd。
S102、实时监测记录套管下至预定井深过程的大钩载荷数据;
所述大钩载荷数据是指在下套管过程中大钩上的载荷值。具体的,从下套管大钩载荷曲线上提取套管下至某井深时的大钩载荷数据P1。
S103、根据套管中扶正器类型及所处位置,确定摩阻系数类型及个数,并在下套管过程中选取与摩阻系数个数相同的时刻数,采集每个时刻的大钩载荷数据及对应的套管串数据;
其中,所述套管串数据至少包括当前套管序号、每根套管长度、套管线重及套管所处位置的井斜角。
可选的,所述摩阻系数为不同类型扶正器及光套管在上层套管内或裸眼段内的摩擦系数。由于扶正器在上层套管内和裸眼段的摩阻系数不一样,且不同类型的扶正器其摩阻系数又不一样,示例性的,假设滚轮扶正器在上层套管内的摩阻系数为μ1,刚性扶正器在上层套管内的摩阻系数为μ2,弹性扶正器在上层套管内的摩阻系数为μ3,光套管在上层套管内的摩阻系数为μ4;滚轮扶正器在裸眼段内的摩阻系数为μ5,刚性扶正器在裸眼段内的摩阻系数为μ6,弹性扶正器在裸眼段内的摩阻系数为μ7,光套管在裸眼段内的摩阻系数为μ8。
进一步的,在每次下套管过程中套管中的扶正器类型会存在差异,故而,不同类型扶正器及光套管与上层套管或裸眼段的摩阻也会不同,即存在不同类型的摩阻系数,摩阻系数类型有多少种,就需要采集几次不同时刻的样本,或采集不同深度的样本数据,具体的,所述样本数据至少包括大钩载荷值、套管编号及套管所述位置的井斜角等,以便于求解方程。
需要注意的是,所述套管数据(或样本数据)还可以包括套管钢材密度等,具体以实际下套管过程及方程求解需要为准,在此不做限定。
优选的,根据套管串中扶正器类型及所处位置,确定摩阻系数类型及个数。
S104、基于每个时刻的大钩载荷数据及对应的套管串数据,列举大钩载荷计算公式(1),得到方程组,
其中,qc为套管线重,ρc为套管钢材密度,ρd为钻井液密度,Li为每根套管长度,i表示套管序号,θi为套管所在井深的井斜角,μk表示不同类型扶正器的摩阻系数,P1为大钩载荷值;
示例性的,在本发明实施例中k可能为1、2、3、4、5、6、7、8中的任何一个数,具体由扶正器类型和安放位置决定:
如第i段套管处于裸眼段内,且安放有弹性扶正器,则k的取值为7。该方程最多具有8个未知数,即μ1、μ2、μ3、μ4,μ5、μ6、μ7、μ8。
假设在下套管的某一时刻,上层套管段内和裸眼段内同时包括刚性扶正器和弹性扶正器,且并非每根套管上都安放有扶正器,即存在光套管接触井壁的现象,则存在上层套管内的摩阻系数μ2、μ3、μ4,裸眼段内的摩阻系数μ6、μ7、μ8,此时刻的大钩载荷方程演变为一个包括μ2、μ3、μ4,μ6、μ7、μ8的六元一次方程,六个不同时刻的大钩载荷方程形成的方程组如下:
假设在下套管的另一时刻,刚性扶正器和弹性扶正器已全部下入到裸眼段内,则存在上层套管内的摩阻系数μ4,裸眼段内的摩阻系数μ6、μ7、μ8,此时步骤三的方程演变为一个包括μ4,μ6、μ7、μ8的四元一次方程。
当为六元一次方程时,至少需要采集6个不同时刻的样本数据求解方程组,当为四元一次方程时,至少需要采集4个不同时刻的样本数据求解方程组。
S105、解所述方程组,得到摩阻系数值。
当所述方程组满足求解条件时,即可求得摩阻系数值。
可选的,通过采集多个下套管样本数据验证所述摩阻系数值,所述样本数据指定是按上述步骤S101和S102采集的下套管数据及大钩载荷数据,如当前下套管过程中其他不同深度的样本数据,或同等条件下,其他下套管过程的样本数据,可以用来验证所述摩阻系数值,进一步保证其准确性。
本实施例中,通过确定摩阻系数类型,进而采样求解方程组,得到摩阻系数值,可以准确的计算不同类型扶正器的摩阻系数,方法简单易操作,可行性高。
为便于理解,在图1所描述的实施例基础上,下面以一个实际应用场景对本发明实施例中的一种下套管摩阻系数计算方法进行描述,具体数据计算如下:
假设某井上层套管深度为500米,500米以上为直径段,下入套管串中每根套管长度Li,每根套管所在井深的井斜角θi,i表示套管序号,套管长度Li以及井斜角θi。
安放扶正器类型如表1所示,钻井液密度值ρd=1.14g/cm3,套管密度ρc=7.85g/cm3,套管线重qc=34.23kg/m。
表1下套管至900米处的相关数据表
当套管下至该井深900米时,监测记录大钩载荷数据或从下套管大钩载荷曲线上提取套管下至该井深时的大钩载荷数据P1=237.49kN。由于该井上层套管500米以内井斜为0,套管串与上层套管不接触,理论上摩阻为0。裸眼段内安放有刚性扶正器、弹性扶正器,所以裸眼段内存在刚性扶正器在裸眼段内的摩阻系数为μ6,弹性扶正器在裸眼段内的摩阻系数为μ7,光套管在裸眼段内的摩阻系数为μ8。
套管下至该井深900米时,根据下式计算此时的大钩载荷数据,令其等于实测的大钩载荷数据P1,得到以下方程:
其中k可能为4、6、7、8中的任何一个数,具体由扶正器类型和安放位置决定。比如i=41时,深度900米处的管串上安放有弹性扶正器,k取值为7;i=40时,井深890米处的管串上未安放任何扶正器,k取值为8;i=20时,井深690米处的管串上安放有刚性扶正器,k取值为6;i=1时,井深500米处及以上井斜为0,套管串与上层套管不接触,摩阻系数不存在。所以该方程包括μ6、μ7、μ8三个未知数,即为三元一次方程。
经计算整理可得:
250.282-(2.668μ6+16.313μ7+22.249μ8)=237.49
进一步整理可得:
2.668μ6+16.313μ7+22.249μ8=12.792
同理,在该井下套管过程中取另外一个不同的时刻(套管下至井深700米处),相关数据表如表2所示。套管下至该井深700米时,监测的大钩载荷为196.103kN,根据下式计算此时的大钩载荷数据,令其等于实测的大钩载荷数据P2,得到以下方程。
井深500米处及以上井斜为0,套管串与上层套管不接触,摩阻系数不存在,500米以下只有弹性扶正器和光套管接触井壁,所以不存在刚性扶正器在裸眼段内的摩阻系数μ6,只有弹性扶正器在裸眼段内的摩阻系数μ7和光套管在裸眼段内的摩阻系数μ8。
经计算整理得到此时的大钩载荷方程:
199.535-(5.442μ7+4.953μ8)=196.103
进一步整理可得:
5.442μ7+4.953μ8=3.432
表2下套管至700米处的相关数据表
按以上思路,在该井下套管过程中取另外一个不同的时刻(套管下至井深800米处),相关数据表如表3所示。套管下至该井深800米时,,监测的大钩载荷为218.134kN,根据下式计算此时的大钩载荷数据,令其等于实测的大钩载荷数据P2,得到以下方程。
井深500米处及以上井斜为0,套管串与上层套管不接触,摩阻系数不存在,裸眼段套管串上安放有刚性扶正器和弹性扶正器,所以存在刚性扶正器在裸眼段内的摩阻系数μ6、弹性扶正器在裸眼段内的摩阻系数μ7和光套管在裸眼段内的摩阻系数μ8。
经计算整理得到此时的大钩载荷方程:
225.382-(0.898μ6+10.222μ7+11.601μ8)=218.134
进一步整理可得:
0.898μ6+10.222μ7+11.601μ8=7.248
表3下套管至800米处的相关数据表
最后整理得到该井在套管下至700m,800m,900m时刻的大钩载荷计算方程组为:
针对以上三元一次方程组,首先将其写成矩阵的求解形式:
针对以上矩阵方程,运用雅可比迭代法,求解可得到裸眼段内的摩阻系数
μ6=0.3427、μ7=0.4354、μ8=0.2146。
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种下套管摩阻系数计算方法,其特征在于,包括:
S1、在套管下至预定井深时,获取套管串中每根套管长度、每根套管所处位置的井斜角、套管线重、套管钢材密度及钻井液密度;
S2、实时监测记录套管下至预定井深过程的大钩载荷数据;
S3、根据套管串中扶正器类型及所处位置,确定摩阻系数类型及个数,并在下套管过程中选取与摩阻系数个数相同的时刻数,采集每个时刻的大钩载荷数据及对应的套管串数据;其中,所述套管串数据至少包括当前套管序号、套管长度、套管线重及套管所处位置的井斜角。
S4、基于所述每个时刻的大钩载荷数据及对应的套管串数据,列举大钩载荷计算公式(1),得到方程组,
其中,qc为套管线重,ρc为套管钢材密度,ρd为钻井液密度,Li为每根套管长度,i表示套管序号,θi为套管所在井深的井斜角,μk表示不同类型扶正器的摩阻系数,P1为大钩载荷值;
S5、解所述方程组,得到摩阻系数值。
2.根据权利要求1所述的方法,其特征在于,所述摩阻系数类型为不同类型扶正器及光套管在上层套管内或裸眼段内的摩擦系数。
3.根据权利要求1所述的方法,其特征在于所述步骤S5还包括:
S6、采集多个不同时刻下套管样本数据,验证所述摩阻系数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056716.6A CN109869132B (zh) | 2019-01-22 | 2019-01-22 | 一种下套管摩阻系数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056716.6A CN109869132B (zh) | 2019-01-22 | 2019-01-22 | 一种下套管摩阻系数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109869132A true CN109869132A (zh) | 2019-06-11 |
CN109869132B CN109869132B (zh) | 2022-10-25 |
Family
ID=66917921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910056716.6A Active CN109869132B (zh) | 2019-01-22 | 2019-01-22 | 一种下套管摩阻系数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109869132B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110566117A (zh) * | 2019-09-06 | 2019-12-13 | 新疆贝肯能源工程股份有限公司 | 一种水平井导向钻具综合阻力系数确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103616117A (zh) * | 2013-11-29 | 2014-03-05 | 中国石油天然气股份有限公司 | 测定三维井眼中抽油杆柱轴向力和侧向力的方法及设备 |
CN103883267A (zh) * | 2014-03-11 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种钻杆用扶正器的布置方法及装置 |
CN105545260A (zh) * | 2015-12-30 | 2016-05-04 | 中国石油天然气股份有限公司 | 一种套管柱下入方法及装置 |
RU2583796C1 (ru) * | 2013-02-27 | 2016-05-10 | Лэндмарк Графикс Корпорейшн | Способ и система для осуществления калибровки коэффициента трения |
CN105825028A (zh) * | 2016-03-31 | 2016-08-03 | 中国石油集团川庆钻探工程有限公司工程技术研究院 | 一种石油钻机大钩载荷的计算方法 |
US20170152737A1 (en) * | 2015-11-30 | 2017-06-01 | Weatherford Technology Holdings, Llc | Calculating Downhole Card in Deviated Wellbore Using Parameterized Segment Calculations |
-
2019
- 2019-01-22 CN CN201910056716.6A patent/CN109869132B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2583796C1 (ru) * | 2013-02-27 | 2016-05-10 | Лэндмарк Графикс Корпорейшн | Способ и система для осуществления калибровки коэффициента трения |
CN103616117A (zh) * | 2013-11-29 | 2014-03-05 | 中国石油天然气股份有限公司 | 测定三维井眼中抽油杆柱轴向力和侧向力的方法及设备 |
CN103883267A (zh) * | 2014-03-11 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种钻杆用扶正器的布置方法及装置 |
US20170152737A1 (en) * | 2015-11-30 | 2017-06-01 | Weatherford Technology Holdings, Llc | Calculating Downhole Card in Deviated Wellbore Using Parameterized Segment Calculations |
CN105545260A (zh) * | 2015-12-30 | 2016-05-04 | 中国石油天然气股份有限公司 | 一种套管柱下入方法及装置 |
CN105825028A (zh) * | 2016-03-31 | 2016-08-03 | 中国石油集团川庆钻探工程有限公司工程技术研究院 | 一种石油钻机大钩载荷的计算方法 |
Non-Patent Citations (2)
Title |
---|
付建红 等: "水平井套管扶正器合理安放位置计算", 《西南石油大学学报(自然科学版)》 * |
*** 等: "多扶正器套管串在弯曲井段下套中的摩阻分析", 《西安石油学院学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110566117A (zh) * | 2019-09-06 | 2019-12-13 | 新疆贝肯能源工程股份有限公司 | 一种水平井导向钻具综合阻力系数确定方法 |
CN110566117B (zh) * | 2019-09-06 | 2021-09-07 | 新疆贝肯能源工程股份有限公司 | 一种水平井导向钻具综合阻力系数确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109869132B (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Izgec et al. | Transient fluid and heat flow modeling in coupled wellbore/reservoir systems | |
Williamson et al. | Representing wells in numerical reservoir simulation: Part 1–Theory | |
EP2656022A2 (en) | Methods and systems for interpreting multiphase fluid flow in a conduit | |
BRPI0706805A2 (pt) | métodos, sistemas e meio legìvel por computador para atualização rápida de modelos de produção para campo de gás e óleo com simuladores proxy e fìsico | |
WO2016134018A2 (en) | Integrated well completions | |
CN105808913B (zh) | 一种页岩气产量递减分析方法 | |
Chokshi et al. | Experimental study and the development of a mechanistic model for two-phase flow through vertical tubing | |
Beirute | The phenomenon of free fall during primary cementing | |
CN109869132A (zh) | 一种下套管摩阻系数计算方法 | |
US20190017378A1 (en) | Natural Resource Reservoir Charging | |
CN109779607A (zh) | 一种基于井径不规则度的下套管摩阻系数确定方法及*** | |
Dejia et al. | Analysis and application of production logging data in shale gas well | |
CN105825028B (zh) | 一种石油钻机大钩载荷的计算方法 | |
CN111444637A (zh) | 一种页岩气长段水平井套管下入的安全性评估方法及*** | |
CN109653730B (zh) | 一种用于钻杆地层测试作业井下射孔井段深度校准方法 | |
BR112019014495A2 (pt) | Método e sistema para modelagem de uma formação subterrânea, e, meio legível por computador não transitório. | |
Onwunyili et al. | Coupled model for analysis of multilayer reservoir in commingled production | |
CN115059454A (zh) | 固井用钻井液顶替效率及冲洗液清洗效率评价装置和方法 | |
CN105089648A (zh) | 一种同位素示踪剂的嵌套式分层注入装置及方法 | |
Semenov et al. | History matching of integrated reservoir simulation model for green field offshore Abu Dhabi | |
CN204347020U (zh) | 一种新型油井钨合金镀层防蜡效果检测装置 | |
US10385656B2 (en) | Methods and systems for determining manufacturing and operating parameters for a deviated downhole well component | |
Al-Hajri et al. | Calculating Casing Leak Cross Flow Rate Using Surface Injection Data | |
Kawaguchi et al. | Optimal sensor placement for multi-phase flow rate estimation using pressure and temperature measurements | |
US10331288B2 (en) | Method and system for generating oilfield objects |
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 |