CN115618695B - 载荷计算模型及建立方法、应用、分析方法、设备、介质 - Google Patents
载荷计算模型及建立方法、应用、分析方法、设备、介质 Download PDFInfo
- Publication number
- CN115618695B CN115618695B CN202211619054.7A CN202211619054A CN115618695B CN 115618695 B CN115618695 B CN 115618695B CN 202211619054 A CN202211619054 A CN 202211619054A CN 115618695 B CN115618695 B CN 115618695B
- Authority
- CN
- China
- Prior art keywords
- torque
- friction
- model
- drill
- establishing
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004458 analytical method Methods 0.000 title claims abstract description 25
- 238000005553 drilling Methods 0.000 claims abstract description 135
- 239000002131 composite material Substances 0.000 claims abstract description 70
- 238000013178 mathematical model Methods 0.000 claims abstract description 44
- 238000003860 storage Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 40
- 239000012530 fluid Substances 0.000 claims description 17
- 229910000831 Steel Inorganic materials 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 15
- 239000010959 steel Substances 0.000 claims description 15
- 239000003381 stabilizer Substances 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 abstract 1
- 239000011435 rock Substances 0.000 description 9
- 239000007789 gas Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000005452 bending Methods 0.000 description 6
- 238000005094 computer simulation Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000010223 real-time analysis Methods 0.000 description 2
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012821 model calculation Methods 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
- 235000012149 noodles Nutrition 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Animal Husbandry (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Mathematical Analysis (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Earth Drilling (AREA)
Abstract
本发明提供了一种载荷计算模型及建立方法、应用、分析方法、设备、介质。所述建立方法包括:对地层中的钻具组合进行分段有限元离散化处理,建立软杆、刚杆和连续梁的物理模型;根据三个物理模型分别建立软杆、刚杆和连续梁的数学模型;根据三个数学模型,建立所述复合计算模型。所述应用为在卡钻预测中的应用。所述分析方法基于所述建立的模型进行。所述设备包括至少一个处理器;存储有程序指令的存储器,程序指令包括执行上述建立方法或评价方法的指令。所述存储介质上存储有程序指令。本发明计算出的摩阻力、扭矩预测值相对误差小,复合模型及其分析方法在现场应用具有可靠性;本发明对卡钻趋势的预测应当具有良好效果。
Description
技术领域
本发明涉及石油钻井技术领域,具体来讲,涉及一种摩阻扭矩复合计算模型及其建立方法、模型的应用、分析方法、计算机设备和存储介质。
背景技术
非常规油气资源勘探开发是我国油气规模上产,保障国家能源安全的重要领域。我国非常规油气资源十分丰富,其中,四川盆地五峰组~龙马溪组海相页岩品质优,分布连续稳定,鄂尔多斯盆地页岩油的地质条件最好,资源总量与储量规模在国内名列前茅。随着国内外石油与天然气勘探开发的不断深入,一些油气藏的开发效果很不理想。而长水平段水平井技术能够更大限度地提高产量和采收率的技术优势,因此国内外各油气田广泛地实施长水平段水平井技术,以提高对特殊油气藏的勘探开发效益。但长水平段水平井的技术难点多、钻井工艺复杂、对装备和工具要求高,所以加大对长水平段水平井技术的研究力度已成国内外石油工业的必然趋势。
长水平段安全快速钻井是缩短钻完井周期、降低作业成本,实现非常规油气资源效益开发的关键所在。超长段水平井施工过程中,随着井深递进,井壁压力增加。基于油气井管柱结构与材料特性,钻具组合可能于深部地层产生“屈曲”现象;这使得管柱变得像面条一样柔软,严重削弱钻头钻进动力,改变钻进方向,增大管柱摩阻扭矩,增加卡钻风险,使后续下套管、固井等施工带来极大困难。
公开号为CN110107224A的发明专利申请公开了一种水平井钻磨管柱下入摩阻计算方法,该方法能够模拟计算出水平井钻磨管柱下入过程中的摩阻,并进行水平井钻磨管柱强度分析。其虽然公开了有限元思想可用于摩阻扭矩计算,但计算过程关联因素单一,尤其是钻具屈曲效应尚未被充分考虑到。因此本发明换一种新模型或新思路,对屈曲效应下的钻具组合进行精细描述,结合现场工程地质数据,得到更准确的摩阻扭矩预测值。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如,本发明的目的之一在于对屈曲效应下的钻具组合进行精细描述,结合现场工程地质数据,得到更准确的摩阻扭矩预测值。
为了实现上述目的,本发明一方面提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法。
所述方法包括以下步骤:
S1、对地层中的钻具组合进行分段有限元离散化处理,建立软杆物理模型、刚杆物理模型和连续梁物理模型;其中,软杆物理模型为忽略钻柱刚度影响的钻杆段的物理模型,刚杆模型为考虑钻柱刚度影响的加重钻杆段的物理模型;
S2、根据所述软杆物理模型,建立第一数学模型,第一数学模型包括软杆条件下的轴向载荷和扭矩载荷计算的模型;
根据所述刚杆模型,建立第二数学模型,第二数学模型包括摩擦力、摩擦扭矩和转盘扭矩计算的模型;
根据所述连续梁模型,建立第三数学模型,第三模型包括下部钻具组合段在连续梁条件下的扭矩载荷计算模型;
S3、根据所述第一、第二和第三数学模型,建立钻井屈曲影响下的摩阻扭矩复合计算模型。
可选择地,所述地层包括可深部地层。
可选择地,所述软杆物理模型的建立条件可包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;仅考虑承受轴向压力。
可选择地,所述刚杆物理模型的建立条件可包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;计算单元柱处在某一空间斜平面上。
可选择地,所述第一数学模型包括:
其中,W为软杆轴向载荷,N;F i+1为第i单元柱下端拉力,N;N g 为单元柱在钻井液中的重力,N;α为第i单元柱上下端井斜角均值,°;δ为第i单元柱上下两端拉力向量的夹角,°;F μ 为软杆单元柱上的摩擦力,N;n为软杆段所划分出的单元柱的总数量;
其中,T为软杆扭矩载荷,N·m;T i+1为第i单元柱下端的扭矩载荷;μ t 为软杆摩擦系数,无量纲;r为第i单元柱半径,m;N为单元柱在钻井液中的正压力,N。
可选择地,所述第二数学模型包括:
其中,F μ 为刚杆单元柱上的摩擦力,N;μ t 为软杆摩擦系数,无量纲;N为单元柱在钻井液中的正压力,N;S为刚杆物理模型中轴线终点的坐标值,用作积分上限,m;M t 为钻柱所受扭矩分布,kN/m;R为钻柱外半径,m;T g 为转盘扭矩,N;q m 为钻柱单位长度重量,kN/m;k f 为钻柱材料密度与钻井液材料密度之比,无量纲;α为第i单元柱上下端井斜角均值,°。
可选择地,所述第三数学模型包括:
其中,T L 为下部钻具组合段在连续梁条件下的扭矩载荷,N·m;E为弹性杨氏模量,kN/m2;I为钻柱惯性矩,m4;K i 为第i单元柱井斜角变化率,无刚量;D O 、D S 分别为井眼和稳定器的直径,m;L为梁长度,m。
可选择地,所述复合计算模型包括:
可选择地,所述方法针对的是定向井。
可选择地,所述方法针对的是大位移井。
本发明另一方面提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型。
所述复合计算模型由根据如上所述的建立方法得到,并包括:
本发明再一方面提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型在卡钻预测中的应用。
本发明又一方面提供了一种考虑钻具屈曲的摩阻扭矩分析方法。
所述方法可包括:利用所上所述的复合计算模型,对目标地层中钻具的摩阻力和/或扭矩进行分析。
可选择地,所述分析可包括:根据所述复合计算模型和目标地层的现场地质数据,得到钻具的总摩阻力和总扭矩。
可选择的,在计算所述总摩阻力和总扭矩的情况下,需满足以下边界条件:
其中,F 0 为第一个单元柱上端拉力;T 0 为第1单元柱上端的扭矩载荷;M为钻井平台转盘实测扭矩数值;WOB为钻压。
本发明又一方面提供了一种计算机设备。
所述设备包括:至少一个处理器;存储有程序指令的存储器,其中,所述程序指令被配置为由所述至少一个处理器执行,所述程序指令包括用于执行根据如上所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法的指令,或者执行如上所述考虑钻具屈曲的摩阻扭矩分析方法的指令。
本发明又一方面提供了一种计算机可读存储介质。
所述储存介质上存储有计算机程序,所述计算机程序指令被处理器执行时实现如上所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,或者实现如上所述考虑钻具屈曲的摩阻扭矩分析方法。
与现有技术相比,本发明的有益效果包括以下内容中的至少一项:
(1)本发明计算出的摩阻力、扭矩预测值分别与实测值对比,相对误差小,复合模型及其分析方法在现场应用具有可靠性。
(2)本发明对卡钻趋势的预测应当具有良好效果,便于现场实时调整钻井参数,有效规避卡钻风险,提高钻井效率。
附图说明
通过下面结合附图进行的描述,本发明的上述和其他目的和/或特点将会变得更加清楚,其中:
图1a示出了本发明的软杆微元段物理模型的受力分解示意图。
图1b示出了本发明的软杆微元段在平面上的受力分解示意图。
图2示出了本发明的刚杆微元段的受力分解示意图。
图3示出了本发明的各支点不在同一条直线上的连续梁受力状况模型。
图4示出了井眼滑动钻进轴向载荷分布图。
具体实施方式
在下文中,将结合示例性实施例来详细说明本发明的载荷计算模型及建立方法、应用、分析方法、设备、介质。其中,荷载计算模型又称为考虑钻具屈曲的摩阻扭矩复合计算模型。
需要说明的是,“第一”、“第二”、“第三”等仅仅是为了方便描述和便于区分,而不能理解为指示或暗示相对重要性。“上”、“下”、“前”、“后”、“左”、“右”、“内”、“外”等仅仅为了便于描述和构成相对的方位或位置关系,而并非指示或暗示所指的部件必须具有该特定方位或位置。
示例性实施例1
本示例性实施例提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法。所述方法可包括以下步骤:
S1、对地层中的钻具组合进行分段有限元离散化处理,建立软杆物理模型、刚杆物理模型和连续梁物理模型。其中,软杆物理模型可以为忽略钻柱刚度影响的钻杆段的物理模型,刚杆物理模型可以为考虑钻柱刚度影响的加重钻杆段的物理模型。
本发明基于分段模型,将地层中的钻具组合分为三个部分进行计算:(1)钻杆段忽略钻柱刚度影响所采用的软杆模型;(2)加重钻杆段考虑钻柱刚度影响所采用的刚杆模型;(3)下部钻具组合段所采用的连续梁模型。基于此,深部地层定向井井下实际情况“从悬挂于钻井平台的正常钻杆,再到加重钻杆段,最后到连接钻头的下部钻具组合段”可以正好对应本发明的“软杆→刚杆→连续梁”。
其中,加重钻杆具有以下几个特点:①其钻杆接头比普通钻杆长,可以提供较大的耐磨表面和重量,接头螺纹可以多次修复;②比同尺寸的钻杆重,壁厚比普通钻杆增加了2-3倍,管体和接头外径与普通钻杆一致,内孔是内平的,内孔直径至少等于钻铤的内径。③中部外加厚段起小型稳定器作用。增加了抗弯曲性能。④接头表面敷焊耐磨的硬质合金,寿命可以延长4倍。⑤钻柱疲劳失效常发生在钻铤以上数根钻杆上,因为从钻铤过渡到钻杆时断面急剧变化,弯曲应力集中在这部分钻杆上。如果在钻铤和钻杆之间加入15~30根加重钻杆,加重钻杆,则可以缓和断面的变化,减少应力集中,防止与钻铤连接的钻杆发生疲劳破坏;从而减少钻具事故。另外,与钻铤相比,加重钻杆还有打捞容易、可缩短起下钻时间、搬运方便等优点,并能保持定向井的方位,起到稳斜作用。在大钩负荷和钻压相同的情况下,使用加重钻杆,还可提高钻机的钻深能力。⑥用于定向井:定向井钻柱中接上加重钻杆后,可使钻铤和加重钻杆本体不紧贴在井壁上,从而减少发生粘吸卡钻的可能性,此外,滑动钻进过程中有利于定向井、水平井井眼轨迹控制。还能更好地控制井斜角和方位角,减少定向工作的次数,从而缩短定向井的建井周期。
在本实施例中,所述地层可包括深部地层。其中,深部地层是指垂深超过3500m的地层,例如垂深为3600m、4000m、6000m的地层。深部地层随岩石围压增加,岩石孔隙度逐渐减小、岩石强度逐渐增加,且二者变化基本趋于平缓;同时,岩石塑性也随之增大。
在本实施例中,所述软杆物理模型的建立条件可包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;仅考虑承受轴向压力。
在本实施例中,所述刚杆物理模型的建立条件可包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;计算单元柱处在某一空间斜平面上。
在本实施例中,连续梁模型可不考虑假设条件。定向井下部钻具组合段主要以钻头、钻铤、稳定器、螺杆等构成,各钻具连接处常以丝扣连接,当考虑实际钻井过程中的稳定器作用时,将各个稳定器视作支座,则下部钻具组合段即可被看为连续梁条件下的一个整体。
S2、根据所述软杆物理模型,建立第一数学模型,第一数学模型包括软杆条件下的轴向载荷和扭矩载荷计算的模型。根据所述刚杆模型,建立第二数学模型,第二数学模型包括摩擦力、摩擦扭矩和转盘扭矩计算的模型。根据所述连续梁模型,建立第三数学模型,第三模型包括下部钻具组合段在连续梁条件下的扭矩载荷计算模型。
在本实施例中,所述第一数学模型可包括以下两式:
其中,W为软杆轴向载荷,N;F i+1为第i单元柱下端拉力,N,其中,i以钻井平台为起始来计;N g 为单元柱在钻井液中的重力,N;α为第i单元柱上下端井斜角均值,°;δ为第i单元柱上下两端拉力向量的夹角,°;F μ 为软杆单元柱上的摩擦力,N;n为软杆段所划分出的单元柱的总数量。
其中,T为软杆扭矩载荷,N•m;T i+1为第i单元柱下端的扭矩载荷,N•m;μ t 为软杆摩擦系数,无量纲,其中,i以钻井平台为起始来计;r为第i单元柱半径,m;N为单元柱在钻井液中的正压力,N。
在本实施例中,所述第二数学模型包括:
其中,F μ 为加重钻杆上的摩擦力,N;
μ t 为软杆摩擦系数,无量纲;
N为单元柱在钻井液中的正压力,N,与岩石围压有关;
S为加重钻杆的总长度,用作积分上限,m;其中,实际的刚杆、软杆、连续梁三者的长度划分可参照邻井数据确定,进而确定S数值;
M t 为钻柱所受扭矩分布,kN/m;
R为钻柱外半径,m;
T g 为刚杆条件下的转盘扭矩,N;
q m 为钻柱单位长度重量,kN/m;
k f 为钻柱材料密度与钻井液材料密度之比,无量纲;
α为第i单元柱上下端井斜角均值,°。
在本实施例中,所述第三数学模型可包括:
其中,T L 为下部钻具组合段在连续梁条件下的扭矩载荷,N·m;E为弹性杨氏模量,kN/m2;I为钻柱惯性矩,m4;K i 为第i单元柱井斜角变化率,无刚量;D O 、D S 分别为井眼和稳定器的直径,m;L为梁长度,m。
S3、根据所述第一、第二和第三数学模型,建立钻井屈曲影响下的摩阻扭矩复合计算模型。
在本实施例中,所述复合计算模型可包括:
示例性实施例2
如图1a示出了软杆的物理模型示意图,该图中字头顶端所带箭头符号表示方向,例如表示单元柱在钻井液中的重力分布方向。图1b示出了软杆在Ri平面上的空间投影分布图,图1b中的希腊字母伽马γ表示:单元柱在钻井液中正压力分布方向与Ri平面上的y轴之间的夹角,仅仅用于展示水平井段受力分布和空间方向。
图2示出了本发明的刚杆微元段的受力分解示意图。图3示出了本发明的各支点不在同一条直线上的连续梁受力状况模型。图4示出了井眼滑动钻进轴向载荷分布图。
本示例性实施例提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法。
本发明提供的技术方案是基于分段模型,将钻具组合分为三个部分进行计算:(1)钻杆段忽略钻柱刚度影响所采用的软杆模型;(2)加重钻杆段考虑钻柱刚度影响所采用的刚杆模型;(3)下部钻具组合段所采用的连续梁模型。考虑钻具屈曲现象,建立摩阻扭矩数值计算复合模型。具体的,包括以下步骤:
S1、对深部地层中的钻具组合进行分段有限元离散化处理;
S2、建立理想化条件下的摩阻扭矩分段数学模型并迭代求解;
S3、建立钻具屈曲影响下的摩阻扭矩复合计算模型。
以下结合具体实施例,对本发明进行详细说明,但所举实例不作为对本发明的限定。
本实施例的技术方案是基于分段模型,将钻具组合分为三个部分进行计算:(1)钻杆段忽略钻柱刚度影响所采用的软杆模型;(2)加重钻杆段考虑钻柱刚度影响所采用的刚杆模型;(3)下部钻具组合段所采用的连续梁模型。在考虑钻具屈曲的影响下,对不同钻柱结构段采用不同计算模型,既较为准确的描述了各段钻具力学特性,又能简化参数,提高摩阻、扭矩的计算精度,满足现场要求。具体包括以下操作步骤:
S1、对深部地层中的钻具组合进行分段有限元离散化处理。
在钻井作业中,钻柱与井壁的摩阻分布于整个钻柱上,其影响主要体现在大钩载荷和扭矩载荷的变化上,为了便于分析,我们将整个钻柱视为若干单元柱的集合(有限微元脚标i=1,…, n),并将深部地层中的钻具组合离散化,基于一定井下实际工况,提出合理假设,将钻具组合三分段,建立软杆、刚杆、连续梁各自物理模型。
进一步方案是,在钻杆段忽略钻柱刚度影响下,建立软杆模型,基本假设为:
(1)计算单元柱井眼曲率是常数;
(2)管柱与井壁接触面上下侧的曲率变化一致;
(3)仅考虑承受轴向压力。
进一步方案是,在加重钻杆段考虑钻柱刚度影响下,建立刚杆模型,基本假设为:
(1)计算单元柱井眼曲率是常数;
(2)管柱与井壁接触面上下侧的曲率变化一致;
(3)计算单元柱处在某一空间斜平面上。
S2、建立理想化条件下的摩阻扭矩分段数学模型并迭代求解。
根据软杆基本假设,通过对每一个单元柱进行受力分析,叠加摩擦力的影响,从而求得大钩载荷和扭矩载荷。在空间斜平面Ri上,以集中力代替分布力,并以单元柱中点为原点建立笛卡儿xyz坐标系,如图1b所示,x轴为切线方向,xy平面与Ri平面重合,z轴与Ri平面垂直向下。在xyz坐标系中,根据力平衡原理有如下方程:
几何投影、迭代计算后,软杆条件下轴向载荷与扭矩载荷分别为:
根据刚杆基本假设,在井眼轴线坐标系上任取一弧长为ds的微元体AB,并对其进行受力分析,以A点为始点,其轴线坐标为s,B点为终点,其轴线坐标为s+ds,此单元体受力如图2所示,在增加弯矩和剪力作用下,可得大位移井全刚度钻柱摩阻扭矩数学模型:
式中:q m 为钻柱单位长度重量,kN/m;k f 为钻柱材料密度与钻井液材料密度之比,无量纲;R为钻柱外半径,m;M t 为钻柱所受扭矩分布,kN/m;M b 为钻柱微段上的均布接触力分布,kN/m;N n 、N b 分别为主法线和副法线方向的压力分布,N/m;α为第i单元柱上下端井斜角均值,°;μ t 为软杆摩擦系数,无量纲;N为单元柱在钻井液中的正压力,N;T为软杆扭矩载荷,N·m; s为刚杆假设中轴线坐标数值,用作积分上限(实际的刚杆、软杆、连续梁三者的长度划分参照邻井数据确定,进而确定s数值),m。
采用基于有限差分思想的拟牛顿迭代法进行求解,将常微分方程离散化,得出法线方向上的均布接触力后,在距钻头任意井深处,推导出其摩阻力、摩擦扭矩、转盘扭矩计算公式:
式中:“±”代表起下钻,起钻取“+”,下钻取“-”。
式中:F μ 注为刚杆摩擦力,N;T g 为刚杆条件下的转盘扭矩,N·m。
当考虑钻井过程中的稳定器作用时,将各个稳定器视作支座,如图3所示,分析受力模式,可建立弯矩方程如下:
式中:L为梁长度,m;、分别为梁左右两端弯矩,N·m;为第i+1单元柱均布载荷,N/m;为弹性杨氏模量,kN/m2;为钻柱惯性矩,m4;、、分别为第i+1单元柱放大因子,可由钻柱力学性能实验给出,无量纲;K为第i+1单元柱井斜与方位角综合变化参数,可由测井数据获得,rad/m;为梁长变化参数,由钻柱力学性能实验给出,无量纲。图3中的M代表弯矩;q代表单元柱所受岩石围压载荷,即钻柱周围岩石对钻柱施加的一个环绕钻柱的压力。
在石油钻井中,支撑井壁的稳定器往往不在同一条直线上,当各支座不在同一条直线上时,简化处理,则可得到下部钻具组合段在连续梁条件下的扭矩载荷:
式中:T L 为下部钻具组合段在连续梁条件下的扭矩载荷,N·m;D O 、D S 分别为井眼、稳定器直径,m。
S3、建立钻具屈曲影响下的摩阻扭矩复合计算模型。
屈曲条件下,钻柱与井壁之间产生的附加接触压力可表达为:
基于以上分析,考虑钻具屈曲现象,建立摩阻扭矩复合模型,其中摩阻力复合计算公式为:
具体求解摩阻扭矩相关数值解时,须满足以下边界条件:
由于钻柱不旋转,在滑动钻进过程中,轴向摩擦力较大。图4为考虑钻具屈曲前后的轴向载荷对比图,其中,曲线a为不考虑钻具屈曲的轴向载荷曲线,曲线b为考虑钻具屈曲的轴线载荷。根据复合模型计算结果,可以看出大部分钻柱处于受压状态,当受压钻柱超过临界载荷时,将发生正弦屈曲或螺旋屈曲。不考虑钻柱屈曲时,井口轴向载荷明显增大,两者最大相差56.62KN。
示例性实施例3
本示例性实施例提供了一种考虑钻具屈曲的摩阻扭矩复合计算模型。
摩阻扭矩复合计算模型可以由示例性实施例1或2中的建立方法得到。
摩阻扭矩复合计算模型可包括:摩阻力复合计算模型和总扭矩复合计算模型。
示例性实施例4
本示例性实施例提供了本发明的复合计算模型的应用。该复合计算模型可以是示例性实施例3所述的模型。
该应用可包括在卡钻预测中的应用。
钻柱摩阻扭矩的实时分析对提高钻井效率、规避钻井卡钻风险具有重要意义,但钻井过程中摩阻扭矩的实时分析尚不成熟,故当前摩阻扭矩分析仍以钻前预测为基础、并行提高实时测算精准度,以规避钻井卡钻风险。考虑到钻柱摩阻扭矩数值在一定程度上表征钻柱卡钻趋势,可以根据本发明的模型计算出的数值变化亦可进一步对钻井卡钻趋势进行预测。
例如,根据本发明的模型发现某井钻柱摩阻、扭矩数值在4000~4300m区间整体呈现逐渐增大的趋势,且在4150m处附近,钻柱摩阻系数从0.35 附近陡增至0.75,变化极为剧烈,判断即将发生卡钻。经过对该井施工日志查证,该井在4155m 处附近蹩停顶驱钻具卡死。说明合理利用摩阻扭矩计算方法及其计算结果,可以有效对卡钻趋势进行准确预判,这也便于现场实时调整钻井参数,有效规避卡钻风险,提高钻井效率。
示例性实施例5
本示例性实施例提供了一种考虑钻具屈曲的摩阻扭矩分析方法。
根据示例性实施例3中的复合计算模型对目标区域的钻具进行摩阻和扭矩的分析。例如,可根据所述复合计算模型以及现场地质条件,得到钻具的总摩阻力和总扭矩。
示例性实施例6
本示例性实施例提供了一种考虑钻具屈曲的摩阻扭矩分析方法。所述方法可包括:
与示例性实施例1或2中相同的步骤S1、S2和S3;
以及步骤S4:根据所述复合计算模型和现场地质数据,得到钻具的总摩阻力和总扭矩。
根据本发明的示例性实施例5或6,在具体求解摩阻扭矩相关数值解时,须满足以下边界条件:
其中,F 0 为第1个单元柱上端拉力;T 0 为第1单元柱上端的扭矩载荷;M为钻井平台转盘实测扭矩数值;WOB为钻压,kN。
根据本发明的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法或考虑钻具屈曲的摩阻扭矩分析方法,可以被编程为计算机程序并且相应的程序代码或指令可以被存储在计算机可读存储介质中,当程序代码或指令被处理器执行时使得处理器执行上述方法,上述处理器和存储器可以被包括在计算机设备中。
示例性实施例7
本示例性实施例提供了一种计算机设备,包括:
至少一个处理器;
存储有程序指令的存储器,其中,所述程序指令被配置为由所述至少一个处理器执行,所述程序指令包括用于执行根据示例性实施例1或2所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法的指令,或者执行根据示例性实施例5或6所述考虑钻具屈曲的摩阻扭矩分析方法的指令。
示例性实施例8
本示例性实施例提供了一种计算机可读存储介质。
所述储存介质上存储有计算机程序。所述计算机程序指令被处理器执行时实现如示例性实施例1或2所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,或者实现示例性实施例5或6所述考虑钻具屈曲的摩阻扭矩分析方法。
该计算机可读存储介质可以是任意数据存储装置,该数据存储装置中存储有能够被计算机系统读出的数据。例如,计算机可读存储介质的示例可包括:只读存储器、随机存取存储器、只读光盘、磁带、软盘、光数据存储装置和载波(诸如经有线或无线传输路径通过互联网的数据传输)。
尽管上面已经结合示例性实施例及附图描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。
Claims (16)
1.一种考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,其特征在于,所述建立方法包括以下步骤:
S1、对地层中的钻具组合进行分段有限元离散化处理,建立软杆物理模型、刚杆物理模型和连续梁物理模型;其中,软杆物理模型为忽略钻柱刚度影响的钻杆段的物理模型,刚杆模型为考虑钻柱刚度影响的加重钻杆段的物理模型;
S2、根据所述软杆物理模型,建立第一数学模型,第一数学模型包括软杆条件下的轴向载荷和扭矩载荷计算的模型;
根据所述刚杆物理模型,建立第二数学模型,第二数学模型包括摩擦力、摩擦扭矩和转盘扭矩计算的模型;
根据所述连续梁物理模型,建立第三数学模型,第三模型包括下部钻具组合段在连续梁条件下的扭矩载荷计算模型;
S3、根据所述第一、第二和第三数学模型,建立钻井屈曲影响下的摩阻扭矩复合计算模型。
2.根据权利要求1所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,其特征在于,所述地层包括深部地层。
3.根据权利要求1所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,其特征在于,所述软杆物理模型的建立条件包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;仅考虑承受轴向压力。
4.根据权利要求1所述的考虑钻具屈曲的摩阻扭矩复合计算模型的建立方法,其特征在于,所述刚杆物理模型的建立条件包括:计算单元柱井眼曲率为常数;管柱和井壁接触面上下侧的曲率变化一致;计算单元柱处在某一空间斜平面上。
11.权利要求10所述的包含考虑钻具屈曲的摩阻扭矩复合计算模型的装置在卡钻预测中的应用。
12.一种考虑钻具屈曲的摩阻扭矩分析方法,其特征在于,所述方法包括:利用权利要求10所述的包含考虑钻具屈曲的摩阻扭矩复合计算模型的装置,对目标地层中钻具的摩阻力和/或扭矩进行分析。
13.根据权利要求12所述的考虑钻具屈曲的摩阻扭矩分析方法,其特征在于,所述分析包括:根据所述复合计算模型和目标地层的现场地质数据,得到钻具的总摩阻力和总扭矩。
14.根据权利要求13所述的考虑钻具屈曲的摩阻扭矩分析方法,其特征在于,在计算所述总摩阻力和总扭矩的情况下,需满足以下边界条件:
在起下钻的情况下:F0=0,T0=0;
在空转的情况下:F0=0,T0=0;
在正常钻进的情况下:F0=-WOB,T0=M;
在滑动钻进的情况下:F0=-WOB,T0=M;
其中,F0为第一个单元柱上端拉力;T0为第1单元柱上端的扭矩载荷;M为钻井平台转盘实测扭矩数值;WOB为钻压。
15.一种计算机设备,其特征在于,包括:
至少一个处理器;
存储有程序指令的存储器,其中,所述程序指令被配置为由所述至少一个处理器执行,所述程序指令包括用于执行根据权利要求1-9中任一项所述方法的指令,或者执行根据权利要求12-14中任一项所述方法的指令。
16.一种计算机可读存储介质,其上存储有计算机程序指令,其特征在于,所述计算机程序指令被处理器执行时实现权利要求1-9中任一项所述的方法,或者实现根据权利要求12-14中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211619054.7A CN115618695B (zh) | 2022-12-16 | 2022-12-16 | 载荷计算模型及建立方法、应用、分析方法、设备、介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211619054.7A CN115618695B (zh) | 2022-12-16 | 2022-12-16 | 载荷计算模型及建立方法、应用、分析方法、设备、介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115618695A CN115618695A (zh) | 2023-01-17 |
CN115618695B true CN115618695B (zh) | 2023-03-28 |
Family
ID=84880181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211619054.7A Active CN115618695B (zh) | 2022-12-16 | 2022-12-16 | 载荷计算模型及建立方法、应用、分析方法、设备、介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115618695B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5431046A (en) * | 1994-02-14 | 1995-07-11 | Ho; Hwa-Shan | Compliance-based torque and drag monitoring system and method |
CN104564019A (zh) * | 2013-10-27 | 2015-04-29 | 中国石油化工集团公司 | 一种大位移井摩阻扭矩监测方法 |
CN105701299A (zh) * | 2016-01-15 | 2016-06-22 | 西南石油大学 | 一种动态摩阻扭矩计算方法 |
CN108804808A (zh) * | 2018-07-16 | 2018-11-13 | 西南石油大学 | 基于相似度理论的复杂结构井动态摩阻扭矩计算方法 |
CN109209336A (zh) * | 2018-11-07 | 2019-01-15 | 西安石油大学 | 一种基于井底钻压的自动送钻系统及控制方法 |
WO2021137866A1 (en) * | 2020-01-02 | 2021-07-08 | Landmark Graphics Corporation | Combined soft and stiff-string torque and drag model |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015094174A1 (en) * | 2013-12-17 | 2015-06-25 | Halliburton Energy Services Inc. | Drilling modeling calibration, including estimation of drill string stretch and twist |
-
2022
- 2022-12-16 CN CN202211619054.7A patent/CN115618695B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5431046A (en) * | 1994-02-14 | 1995-07-11 | Ho; Hwa-Shan | Compliance-based torque and drag monitoring system and method |
CN104564019A (zh) * | 2013-10-27 | 2015-04-29 | 中国石油化工集团公司 | 一种大位移井摩阻扭矩监测方法 |
CN105701299A (zh) * | 2016-01-15 | 2016-06-22 | 西南石油大学 | 一种动态摩阻扭矩计算方法 |
CN108804808A (zh) * | 2018-07-16 | 2018-11-13 | 西南石油大学 | 基于相似度理论的复杂结构井动态摩阻扭矩计算方法 |
CN109209336A (zh) * | 2018-11-07 | 2019-01-15 | 西安石油大学 | 一种基于井底钻压的自动送钻系统及控制方法 |
WO2021137866A1 (en) * | 2020-01-02 | 2021-07-08 | Landmark Graphics Corporation | Combined soft and stiff-string torque and drag model |
Also Published As
Publication number | Publication date |
---|---|
CN115618695A (zh) | 2023-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11286766B2 (en) | System and method for optimizing tubular running operations using real-time measurements and modelling | |
Fazaelizadeh | Real time torque and drag analysis during directional drilling | |
Vlasiy et al. | Improving the aluminum drill pipes stability by optimizing the shape of protector thickening | |
Guan et al. | Dynamic behavior analysis of push-the-bit rotary steerable bottom hole assembly | |
Mirhaj et al. | Improvement of torque-and-drag modeling in long-reach wells | |
Zhang et al. | Buckling configurations of tubular strings constrained in three-dimensional curved wellbores | |
CN106639879A (zh) | 一种底部钻具组合造斜率预测方法 | |
Mitchell et al. | Lateral Buckling—The Key to Lockup | |
CN115618695B (zh) | 载荷计算模型及建立方法、应用、分析方法、设备、介质 | |
Mahyari et al. | Drill string instability reduction by optimum positioning of stabilizers | |
Li et al. | Fundamental equations for dynamic analysis of rod and pipe string in oil-gas wells and application in static buckling analysis | |
CN110457866B (zh) | 滑动钻进全过程摩擦阻力预测方法、降低摩擦阻力方法 | |
Ding et al. | Determination of effective drag forces and torques for jam release during drilling and workover operations in directional wells | |
Fan et al. | Method of suspender line trajectory design | |
Zhang et al. | On the failure configuration of suspended tubular strings when helical buckling just occurs | |
CN115270345B (zh) | 基于连续梁绳索取心防斜钻具力学模型的计算方法 | |
Cayeux et al. | Estimation of weight and torque on bit: assessment of uncertainties, correction and calibration methods | |
Zhu et al. | Vibration analysis of a drillstring in horizontal well | |
Dareing et al. | Tubular bending and pull-out forces in high-curvature well bores | |
Ding et al. | Analysis of Casing String Running Characteristics in Negative-Displacement Horizontal Wells | |
CN117932795A (zh) | 一种水平井钻井延伸极限的预测方法 | |
Perneder et al. | Steady-state solutions of a propagating borehole: Helical trajectory | |
Lei | Downhole Weight on bit prediction with analytical model and finite element method | |
Ifeanyi et al. | Predicting drillstring buckling | |
CN118690887A (zh) | Bha的钻进造斜率预测方法、装置、设备及存储介质 |
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 |