一种简谐冲击下的钻柱动力学特性的计算方法技术

本发明专利技术涉及的是一种简谐冲击下的钻柱动力学特性的计算方法方法，它包括如下步骤：一、将钻杆、钻铤、底部钻具组合、钻头等效为集中质量的二自由度轴向与扭转相互耦合的物理模型；二、基于振动力学的相关理论，采用牛顿定理描述二自由度轴向与扭转相互耦合的数学模型，建立钻柱动力学数学模型；三、采用状态依赖时间迟滞的方法，对数学模型进行修正；四、根据所述钻柱动力学特性，调整外加简谐动载的幅值与角频率，分析确定钻柱系统动力学特性。本发明专利技术分析了外加简谐冲击幅值与角频率对钻柱动力学的影响，为控制钻柱运动特性、提速降本、确保井下安全提供了理论指导。下安全提供了理论指导。下安全提供了理论指导。

【技术实现步骤摘要】
一种简谐冲击下的钻柱动力学特性的计算方法

：
[0001]本专利技术涉及的是石油钻井
，具体涉及的是一种简谐冲击下的钻柱动力学特性的计算方法。

技术介绍
：
[0002]近年来，我国石油行业在产能与规模方面发展迅速，而钻井技术作为石油勘探与开发的重要手段，其发展与相应理论的研究对我国石油行业具有重要意义。
[0003]从目前的钻井工程发展来看，深井、超深井数量激增。随着深度的增加，钻柱不再能被视作刚体，更易产生粘滑、跳钻问题，甚至会出现混沌现象。扭转方向的粘滑振动，会导致井下钻具过早疲劳失效；轴向方向跳钻问题，会大大降低钻井效率，增加成本。因此，为了有效避免粘滑与跳钻问题，提高机械钻速，降低钻井成本，保证钻井作业的安全，有必要对钻柱动力学特性进行分析。
[0004]目前，对钻柱动力学的研究以扭转振动和轴
‑
扭复合振动为主。研究扭转振动主要是研究扭转方向的粘滑效应，该方法基于振动力学相关理论，采用哈密顿原理的拉格朗日方程建立动力学模型，分析相关参数对角速度、角位移的影响。研究轴
‑
扭复合振动是在扭转振动的基础上，考虑轴向振动，此时钻头所受的钻压与扭矩相互耦合。
[0005]然而，现有的研究方法仅是针对最原始的冲击钻井技术提出的。随着钻井技术的不断发展，许多学者提出了能够高效破岩的新方法和工具，如共振钻井和水力振荡器等，简谐冲击技术开始走入钻井从业人员的视野。因此，钻井动力学的相关研究也不能停滞不前，亟需一种能够符合当前高效破岩发展，对简谐冲击下的钻柱动力学进行分析的方法。

技术实现思路
：
[0006]本专利技术的目的是提供一种简谐冲击下的钻柱动力学特性的计算方法，这种简谐冲击下的钻柱动力学特性的计算方法用于解决目前扭转振动和轴
‑
扭复合振动方法破岩效率较低的问题。
[0007]本专利技术解决其技术问题所采用的技术方案是：这种简谐冲击下的钻柱动力学特性的计算方法：
[0008]步骤一、将钻杆、钻铤、底部钻具组合、钻头等效为集中质量的二自由度轴向与扭转相互耦合的物理模型；
[0009]步骤二、基于振动力学的相关理论，采用牛顿定理描述二自由度轴向与扭转相互耦合的数学模型，建立简谐冲击下的轴向
‑
扭转复合动力学数学模型为：
[0010][0011][0012]式中为角加速度、为角速度、Φ为角位移、C
t
为扭转阻尼系数、K
t
为扭转刚度系
数、Ω0为转盘转速、T为钻头所受扭矩、为轴向加速度、为轴向速度、X为轴向位移、C
a
为轴向阻尼系数、K
a
为轴向刚度系数、v0为大钩运动速度、W0为稳态钻压、W为钻头所受钻压、F0为外加简谐冲击幅值、w为外加简谐冲击频率；
[0013]步骤三、采用状态依赖时间迟滞的方法，对简谐冲击下的轴向
‑
扭转复合动力学数学模型进行修正；
[0014]步骤四、根据所述钻柱动力学特性，调整外加简谐动载的幅值与角频率，分析确定钻柱系统动力学特性。
[0015]上述方案中步骤二中的T和W求取方法：
[0016]在考虑顺应性的基础上，钻压和扭矩为：
[0017][0018][0019]式中，T
c
为扭矩切削分量；T
f
为扭矩摩擦分量；W
c
为钻压切削分量；W
f
为钻压摩擦分量；d为钻头切削深度；为关于角速度的单位阶跃函数、为关于轴向速度的单位阶跃函数、H(d)为关于切削深度的单位阶跃函数；为关于角速度的符号函数、为关于轴向速度的符号函数、sgn(d)为关于切削深度的符号函数，分别定义为：
[0020]T
c
＝0.5εR
b
nd
n
，T
f
＝0.5R
b2
γμlσ
[0021]W
c
＝ξεR
b
nd
n
，W
f
＝R
b
lσ
[0022][0023][0024]式中，ε为岩石的破碎比能，R
b
为钻头半径，n为钻头刀翼数量，d
n
为每个刀翼的垂直切削深度，ξ为与切削力方向有关的参数，γ为大于1的钻头几何特征参数，μ为摩擦系数，l为钻头磨损平面长度，σ为接触应力。
[0025]上述方案中每个刀翼的垂直切削深度，根据轴向位移X(t)与状态依赖时滞变量t
n
、t
n
前的轴向位移X(t
‑
t
n
)与时间变量t计算得出：
[0026]d
n
＝X(t)
‑
X(t
‑
t
n
)
[0027]上述方案中状态依赖时间迟滞的方法中状态依赖时滞变量t
n
，根据以下表达式求出：
[0028][0029]其中Φ(t)、Φ(t
‑
t
n
)分别为角位移与t
n
前的角位移。
[0030]本专利技术有益效果在于：
[0031]1.本专利技术将钻杆、钻铤、BHA、钻头等效为集中质量，基于振动力学相关理论，采用牛顿定律建立符合高效破岩发展的在简谐冲击作用下的轴向
‑
扭转二自由度动力学模型，考虑顺应性，对模型进行修正。为了更加贴合实际工程状况，引入状态依赖时间迟滞变量，采用更为精确的四阶龙格
‑
库塔方法，对钻柱动力学分析。
[0032]2.本专利技术分析了外加简谐冲击幅值与角频率对钻柱动力学的影响，为控制钻柱运动特性、提速降本、确保井下安全提供了理论指导。
[0033]3.采用扭转振动模型和轴
‑
扭复合模型不能很好的分析符合高效破岩发展的钻柱动力学特性问题，本专利技术提供一种在简谐冲击下对井下钻柱动力学复杂特性的分析方法，通过建立轴向与扭转二自由度状态依赖时间迟滞动力学模型，考虑外加简谐动载的幅值与频率，依靠四阶龙格
‑
库塔方法对问题进行分析。
附图说明
[0034]图1为钻柱动力学系统轴向振动模型图；
[0035]图2为钻柱动力学系统扭转振动模型图；
[0036]图3为不同角频率对角速度影响图；
[0037]图4为不同角频率对轴向速度影响图；
[0038]图5为不同幅值对角速度影响图；
[0039]图6为不同幅值对轴向速度影响图。
具体实施方式：
[0040]下面结合附图对本专利技术做进一步的说明：
[0041]这种简谐冲击下的钻柱动力学特性的计算方法：
[0042]步骤一、将钻杆、钻铤、底部钻具组合、钻头等效为集中质量的二自由度轴向与扭转相互耦合的物理模型。
[0043]步骤二、基于振动力学的相关本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种简谐冲击下的钻柱动力学特性的计算方法，其特征在于包括如下步骤：步骤一、将钻杆、钻铤、底部钻具组合、钻头等效为集中质量的二自由度轴向与扭转相互耦合的物理模型；步骤二、基于振动力学的相关理论，采用牛顿定理描述二自由度轴向与扭转相互耦合的数学模型，建立简谐冲击下的轴向
‑
扭转复合动力学数学模型为：扭转复合动力学数学模型为：式中为角加速度、为角速度、Φ为角位移、C
t
为扭转阻尼系数、K
t
为扭转刚度系数、Ω0为转盘转速、T为钻头所受扭矩、为轴向加速度、为轴向速度、X为轴向位移、C
a
为轴向阻尼系数、K
a
为轴向刚度系数、v0为大钩运动速度、W0为稳态钻压、W为钻头所受钻压、F0为外加简谐冲击幅值、w为外加简谐冲击频率；步骤三、采用状态依赖时间迟滞的方法，对简谐冲击下的轴向
‑
扭转复合动力学数学模型进行修正；步骤四、根据所述钻柱动力学特性，调整外加简谐动载的幅值与角频率，分析确定钻柱系统动力学特性。2.根据权利要求1所述的简谐冲击下的钻柱动力学特性的计算方法方法，其特征在于：步骤二中的T和W求取方法：步骤二中的T和W求取方法：式中，T
c
为扭矩切削分量；T
f
为扭矩摩擦分量；W
c
为钻压切削分量；W
f
为钻压摩擦分量；d为钻头切削深度；为关于角速度的单位阶跃函数、为关于轴向速度的单位阶跃函数、H(d)为关于切削深度的单位阶跃函数；为关于角速度的符号函数、为关于轴向速度的符号函数、...

【专利技术属性】
技术研发人员：李思琪，陈卓，孙文睿，童叶霜，李玮，赵欢，
申请(专利权)人：东北石油大学，
类型：发明
国别省市：