基于有限元的下部钻具组合优化研究.docxVIP

基于有限元的下部钻具组合优化研究 调整部分排放系统的克氏分析,确定其是和其发展过程 在钻井过程中，适当的钻孔压力是提高钻井效率的有效手段。随着井眼的加深，钻柱变厚。当钻孔压超过下钻具组合（boleassembly，简称b.a）的临界失稳负荷时，钻头失去的稳定性，也被称为结构畸变。这将有可能导致钻柱的失效、钻井效率的降低以及造成井斜。下部钻具组合的失稳也许不会导致钻柱的立即失效,但由于其引起的钻具疲劳和磨损将导致钻柱的过早失效。钻柱失稳后将导致钻柱在井眼里成蛇形或螺旋形,使钻压不能加在钻头上,从而造成钻井效率的降低。 增强下部钻具的稳定性可以通过安放扶正器来实现,然而扶正器安放的位置将直接影响下部钻具组合所能承受临界失稳载荷的大小,在后面的有限元分析中我们可以发现,对于同样的钻具组合,通过调整扶正器的安放位置,可以实现尽量大的临界失稳载荷,从而达到防斜打快的目的。 钻柱的屈曲模型最先由Lubinski在50年代提出,1984年Bogy和Paslay提出斜井中的屈曲模型,其他的屈曲模型基于欧拉模型,大多数屈曲模型都是基于临界载荷假设。1977年,白家祉提出了纵横弯曲连续梁理论,用于转盘钻井下部钻具组合的受力分析;苏义脑院士将纵横弯曲连续梁理论发展至可求解井下动力钻具、带弯角、变截面组合的三维分析,拓宽并完善了纵横弯曲法的理论体系。Millheim K K自1978年来先后发表了多篇使用有限元法求解下部钻具组合的受力与变形分析的论文,他在对下部钻具组合静力分析的基础上引入了惯性力和摩擦力,建立了对BHA的动力分析有限元模型。 通过有限元软件,我们可以进行特征值屈曲分析,得到结构的理论临界屈曲载荷的上限值,通常,考虑非线性因素影响,这个值比实际临界屈曲载荷要大30%～50%。在结构受压情况下,当钻压增加时,钻柱的刚度弱化效应也随之增加,当达到某个载荷时,弱化效应超过钻柱结构的固有刚度,此时没有了净刚度,位移无限增加,结构发生失稳。基本思路是:采用有限元分析软件对包括钻铤和扶正器的整个中和点以下的下部钻具组合进行特征值屈曲分析,求得特征值屈曲荷载因子,从特征值屈曲模态可以分析不同扶正器及钻铤组合的稳定性。我们采用大型商用有限元分析软件ANSYS进行下部钻具组合在钻压作用下的失稳分析。 1 特征值屈曲分析 ANSYS的特征值屈曲分析考虑了应力刚化效应,这种效应会导致结构在承受应力后抵抗横向载荷能力降低,当压力增加时,结构抵抗横向变形能力减小。在某一载荷水平下,这种负的应力刚度超过线性结构刚度,造成结构屈曲。屈曲发生的点称为分叉点(或歧点),由于力-变形曲线达到该点后可能沿两条不同途径前进,所以当超过分叉点时,结构将屈曲或者在不稳定状态下承受载荷。 ANSYS的特征值屈曲分析使用特征值的公式计算造成负刚度的应力刚度矩阵的比例因子。公式如下 ([K]+λ[S]){ψ}=0 式中,[K]= 刚度矩阵;[S]= 应力刚度矩阵;{ψ}= 位移特征矢量;λ= 特征值(也叫做比例因子或载荷因子,即分析图形中的FREQ)。 利用上面的特征值公式可以决定结构的分叉点,具有分叉屈曲的结构在达到屈曲载荷之前其位移-变形曲线表现出线性关系,达到屈曲载荷之后,曲线将跟随另外的路线。 需要说明的是,在下部钻具组合的稳定性分析中,如果给定的载荷为单位载荷,特征值即是屈曲载荷,我们在分析中考虑钻柱的自重以及泥浆浮力的影响,而不是一个恒定的单位集中力,因此得到的特征值只是载荷因子,结构的临界屈曲载荷是载荷因子与钻压的乘积,但这只是一个理论解,可以利用特征值屈曲分析预测出屈曲载荷的上限值,其特点是计算速度快,因此在考虑实际情况时应有一个安全系数。 屈曲分析的基本步骤一般分以下三步进行:首先设置预应力选项进行静力分析;计算屈曲特征值;计算特征向量。具体如下: (1) 建立模型; (2) 完成静力分析:注意将求解预应力选项打开; (3) 分析模型特征值; (4) 分析模型特征向量:注意控制选项计算应力等参数; (5) 进入后处理菜单获得计算结果等。 (6) 评价分析结果。 2 防斜能力计算分析 对某油田目前正在使用的大尺寸井眼(17.5″)钻具组合,进行有限元钻柱受力和变形的计算分析,分析下部钻柱中心线产生最大位移和应力的位置,对其进行防斜机理评价分析。我们对该油田大井眼钻井使用的两种典型钻具结构进行防斜能力的力学计算分析。 计算使用的参数有 钻压:牙轮钻头20t; 泥浆比重:1.1g/cm3; 井斜角:0°; 钻挺内径: 11″ 为76.2mm, 9″ 为76.2mm,8″ 为71.44mm; 钢的比重:7.8g/cm3。 2.1 无排液钻具结构 17.5″钻头(0.5m)+接头(0.5m)+9″钻铤(9m)×3根+扶正器1(2m)+9″钻铤(9m)×1根+扶正器2(2m)