A Study of the Dynamic Analysis and Design Method of Deepwater Drilling String Systems
-
摘要:
深水钻井过程中，从平台延伸到井底的钻柱会在不同深度处与隔水管或井筒发生多点碰撞和摩擦，呈现出接触非线性特征。为了准确掌握深水钻井管柱系统的非线性动力学特性，将其简化为管中管结构，并提出了对应的管柱动力学模型。采用Abaqus有限元软件，对建立的管柱动力学模型进行动态响应模拟，并将模拟结果导入Isight优化软件，进行基于可靠度分析的多目标优化设计，确定出在工程可行性和安全可靠性方面都满足要求的设计参数组合。研究发现，相比于单独考虑隔水管的模型，提出的管中管模型所模拟的管柱系统整体变形程度较小，说明内外管柱之间的相互作用对深水钻井管柱系统的整体偏移有抑制作用；此外，提出的基于可靠度分析的多目标优化设计方法，可以避免优化设计结果因靠近约束边界而在参数波动情况下失效的问题。
Abstract:During deepwater drilling operations, the drilling string extending from a platform to the bottom hole will collide and rub with the riser or wellbore at multiple points at different depths, presenting the characteristics of nonlinear contact. In order to accurately understand the nonlinear dynamic characteristics of deepwater drilling string systems, we proposed a dynamic model that based on the structure of pipe-in-pipe. The dynamic response simulation was conducted by using Abaqus finite element software, and then the simulation results were imported to Isight optimization software to carry out the multi-objective optimization design based on reliability analysis to determine the combination of design parameters that meet the requirement of engineering feasibility and safety reliability. The study suggests that the overall deformation of the string system from the proposed pipe-in-pipe model is smaller when compared with a model that only considers the riser, indicating that the overall offset of deepwater drilling string system is restrained by the interaction between the inner and outer strings. In addition, the proposed multi-objective optimization design method based on reliability analysis can effectively avoid the failed optimization design caused by parameter fluctuations in the vicinity of the constraint boundary.
-
深水钻井与陆地钻井最大的区别是，钻具和钻井液要通过几百甚至数千米长的隔水管系统到达海底实现钻进[1]。隔水管系统上端连接钻井平台，下端通过海底防喷器与井筒相连。钻井过程中，隔水管和井筒与钻柱不断发生碰撞、摩擦等相互作用。因此，对这一管中管系统进行准确的动力学分析和优化设计是确保深水钻井安全的关键，具有重要的工程实际意义。文献调研发现：相关研究大多侧重于对隔水管本身的力学分析[2]，一般采用时域分析[3]或频域分析[4–6]来展开。在建立物理模型时，将隔水管系统视为一个上端自由、下端铰接的梁[7]，然后用Morison方程[2]来描述其外部所受的海浪、海流等作用。部分模型考虑了隔水管内部钻井液流动的影响[8]，有些模型也考虑了其内部钻柱振动的影响[9]。建立的模型大多采用有限差分和有限元法进行求解，并主要分析了顶张力、海水深度、平台漂移、浮力块配置和海况条件等因素对管柱系统横向弯曲变形的影响[10–11]。总体而言，现有研究对深水钻井管柱系统内外管柱相互作用的接触非线性的影响考虑较少，对设计参数组合优化以及优化结果可靠性的分析也不够。为此，笔者建立了管中管结构模型，分析了深水钻井管柱系统的动力学特性，并在此基础上结合可靠度分析进行多目标优化，提出了深水钻井管柱系统的优化设计方法。
1. 深水钻井管柱系统动力学模型
深水钻井管柱系统动力学模型主要包括管中管模型和内外管柱相互作用模型（见图1，图中：μ为摩擦系数；K为碰撞刚度，N/m；N－i和N－i+1为管柱分段编号），前者主要用于描述隔水管和钻柱的弯曲变形以及隔水管外部受到的海洋环境载荷，后者主要用于描述内外管柱之间的相互作用。这2部分统筹起来形成的整体模型包含2个独立坐标系（即O-xyz和O′-x′y′z′），分别用于描述外管和内管系统。
在对深水钻井管柱系统进行动力学建模时，首先建立隔水管和钻柱微元段的运动微分方程，然后将离散的微元段向整体坐标系转化，经拼装得到整个深水钻井管柱系统的有限元模型。就隔水管微元段而言，需对其外部的海洋环境载荷进行描述，采用Ekman漂流理论[1]计算海流流速，然后采用Airy波理论[1]计算海浪的速度和加速度，根据计算的海流和海浪参数，采用Morison方程[2]计算隔水管受到的由海洋环境载荷引起的拖曳力和惯性力。对于隔水管的弯曲变形，采用Euler-Bernoulli梁[2]进行描述。假设管柱系统在y方向发生弯曲，考虑弯曲截面剪切力、张紧力和重力三者在y方向的分量，同时根据达朗伯原理[2]引入惯性力。此外，为进一步考虑内外管柱相互作用的影响，当内外管柱发生接触时，它们之间的相互作用由沿管柱系统引入的N个弹簧–摩擦单元来模拟，并通过设置Heaviside方程来判断是否触发接触位置的弹性力和摩擦力。由此得到隔水管微元段的运动微分方程为：
12CdDρ(uw+uc−∂y∂t)|uw+uc−∂y∂t|+14D2πρ[Cmaw−(Cm−1)∂2y∂t2]+∂2∂z2[EI(z)∂2y∂z2]+p(z)∂2y∂z2+(1−Bf)m(z)g∂y∂z+m(z)∂2y∂t2+H(y−y′−G)[K(z)(y−y′−G)+μ(z)K(z)(y−y′−G)]=0 (1) 其中
H(y−y′−G)={0(y−y′<G)1(y−y′⩾ (2) 式中：Cd为拖曳力系数；D为隔水管外径，m；ρ为海水密度，kg/m3；uw为海浪速度，m/s；uc为海流流速，m/s；Cm为惯性系数；aw为海浪加速度，m/s2； EI为隔水管抗弯刚度，N·m2；p为张紧力，N；m为隔水管的线质量，kg/m；Bf为浮力块浮力系数；g为重力加速度，m/s2；H为Heaviside方程；G为初始间距，m。
同理可得钻柱微元段的运动微分方程为：
\begin{array}{l} \dfrac{{{\partial ^2}}}{{\partial z{{\rm{'}}^2}}}\left[\! {EI{\rm{'}}\left(\! {z{\rm{'}}} \!\right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial z{{\rm{'}}^2}}}} \right] \!+\! p{\rm{'}}\left(\! {z{\rm{'}}}\! \right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial z{{\rm{'}}^2}}} \!+\! m{\rm{'}}\left(\! {z{\rm{'}}}\!\right)g\dfrac{{\partial y{\rm{'}}}}{{\partial z{\rm{'}}}} + \\ m{\rm{'}}\left( {z{\rm{'}}} \right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial {t^2}}} + H\left( {y - y' - G} \right) \left[ {K\left( z \right)}\right. \left( {y' - y - G} \right) +\\\quad\qquad\qquad\left. {\mu \left( z \right)K\left( z \right)\left( {y' - y - G} \right)} \right] = 0 \end{array} (3) 式中：
EI' 为钻柱抗弯刚度，N·m2；m {\rm{'}} 为钻柱的线质量，kg/m；p' 为钻柱悬重，N。设置初始条件时，隔水管与钻柱的初始位移与速度均设为0，即：
\left\{ {\begin{array}{l} \!\!\!{y\left( {z,0} \right) = 0}\\ \!\!\!{\dfrac{{\partial y\left( {z,0} \right)}}{{\partial t}} = 0}\\ {\begin{array}{l} \!\!\!\!\!\!\!{y'\left( {z',0} \right) = 0}\\ \!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {z',0} \right)}}{{\partial t}} = 0} \end{array}} \end{array}} \right. (4) 设置边界条件时，隔水管上端设置为自由、下端固定于井口，上下两端由于存在挠性接头均可旋转；钻柱设置为上端随隔水管偏移、下端受岩石反作用力影响进行轴向振动，上下两端均不能旋转，即：
\left\{ {\begin{array}{l} {y\left( {0,t} \right) = 0}\\ {EI\left( 0 \right)\dfrac{{{\partial ^2}y\left( {0,t} \right)}}{{\partial {z^2}}} = {k_{\rm{b}}}\dfrac{{\partial y\left( {0,t} \right)}}{{\partial z}}}\\ {\begin{array}{l} \!\!\!\!{y\left( {L,t} \right) = S}\\ \!\!\!\!\!{EI\left( L \right)\dfrac{{{\partial ^2}y\left( {L,t} \right)}}{{\partial {z^2}}} = {k_{\rm{t}}}\dfrac{{\partial y\left( {L,t} \right)}}{{\partial z}}}\\ {\begin{array}{l} \!\!\!\!\!\!\!\!\!{y'\left( {0,t} \right) = A\sin \left( {\varOmega t} \right)}\\ \!\!\!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {0,t} \right)}}{{\partial z'}} = 0}\\ {\begin{array}{l} \!\!\!\!\!\!\!\!\!\!\!\!\!{y'\left( {L,t} \right) = S}\\ \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {L,t} \right)}}{{\partial z'}} = 0} \end{array}} \end{array}} \end{array}} \end{array}} \right. (5) 式中：A为钻头振幅，m；Ω为钻头振动角速度，rad/s；S为隔水管顶部偏移，m；kt为上挠性接头旋转刚度，（N·m）/rad；kb为下挠性接头旋转刚度，（N·m）/rad。
2. 动力学模型的有限元分析
由于建立的动力学模型考虑了内外管柱之间的相互作用，模型的求解难度和计算时间显著增加。为了更加高效，利用有限元软件Abaqus进行动态模拟分析。在Abaqus中，采用B31型梁单元描述内外管柱的弯曲变形。将井筒简化为一个与隔水管外径相等且厚度为0.10 m的梁，认为当其与钻柱接触时发生弹塑性变形。内外管柱之间的相互作用使用general-contact模块中的edge-to-edge模型来设定。隔水管顶部和底部分别与钻井平台和井口铰接，采用Abaqus中的CONN3D2型铰链模拟顶部和底部的铰链，铰链的连接方式选用Joint-Rotation组合。隔水管外部受到的海洋环境载荷利用Aqua模块施加，Aqua模块算法的基础是Morison方程[2]。Aqua模块中的海流为分层运动，南海海洋监测数据分析表明，从海面到海底的海流流速从1.35 m/s到0 m/s呈指数式递减，符合Ekman漂流理论[2]。建立有限元模型后，首先进行静态分析，然后进行时长100个海浪周期的动态分析，以确保管柱系统达到动态稳定状态。
有限元模拟用到的主要参数：隔水管和钻柱的管材密度均为7 850 kg/m3，弹性模量均为210 GPa；钻柱的外径、壁厚分别为149.20和9.17 mm，钻柱的悬重比为0.9；钻头的振幅为0.02 m，隔水管的外径、壁厚分别为533.4和19.1 mm，隔水管的张力比为1.2；内外管柱间的碰撞刚度为18 000 kN/m，摩擦系数为0.05；上、下挠性接头的旋转刚度分别为5 000和7 300 kN·m/rad，浮力块的浮力系数为0.9；海水密度为1 050 kg/m3，钻井液密度为1 200 kg/m3；拖曳力系数为0.7，惯性力系数为2.0；海浪振幅为3.0 m，海浪周期为9.0 s，海面风速为5.0 m/s，风速因子为0.03。基于以上参数，使用Abaqus软件进行数值模拟，结果见图2。其中，图2（a）所示为深水钻井管柱系统达到动态稳定时的空间姿态，为便于观察，管柱横向变形放大了10倍，图中红色部分是管柱横向偏移较大的位置，其主要位于管柱中上部的1/3处。图2（b）为图2（a）的局部截图，展示了内外管柱之间发生碰撞、摩擦时的情况。由于在建立模型时考虑了内外管柱间的相互作用，因此首先比较了这种相互作用对泥线以上管柱整体偏移的影响，以2 000 m水深、1 000 m钻进深度为例进行数值模拟，分别模拟管中管系统和不考虑钻柱的隔水管系统的动态响应，结果见图2（c）。
由图2（c）可知，相比于单独考虑隔水管的情况，将内外管柱的碰撞和摩擦等相互作用纳入模型考虑时，泥线以上管柱系统的整体偏移量受到了明显抑制。这说明对于存在管中管结构的钻井管柱系统，在给定的安全偏移范围内能够承受更大的海洋环境载荷。
3. 参数敏感性分析
由上文所建有限元模型可知，主要设计参数包括隔水管外径与壁厚、顶部偏移量、钻压和顶张力。为了分析这些设计参数对钻井管柱系统动态响应的影响，将图2所示模拟结果作为标准状态，每次只改变1个设计参数重复进行数值模拟，当管柱系统达到动态稳定时，导出计算结果进行分析。5个模拟实例的计算结果与标准状态的对比如图3所示。
![]()
图 3 主要设计参数对钻井管柱系统动态响应的影响Figure 3. Effect of main design parameters on the dynamic response of drilling string system由图3可知，隔水管外径与壁厚、顶部偏移、钻压和顶张力等参数对管柱系统的动态响应都有不同程度的影响。增大隔水管壁厚、减小隔水管顶部偏移量、钻压及增大顶张力，都可以减小管柱系统的整体弯曲，并能同时减小隔水管及上、下挠性接头的转角。特别是，由于隔水管外部受到较强海洋环境载荷的影响，其应力分布呈现明显的弯曲效应；而钻柱除了在井口处由于与井壁接触产生应力集中之外，其余主要体现为轴向应力，因此呈现出近似线性分布的特征。
除了前面提到的4项主要设计参数，上、下挠性接头的旋转刚度也是深水钻井管柱系统配置的关键。因此，需分析上、下挠性接头旋转刚度组合对钻井管柱系统动态响应的影响。旋转刚度组合从（10 kN·m/rad，150 kN·m/rad）到（5 000 kN·m/rad，7 300 kN·m/rad）选取4组，模拟其对钻井管柱系统动态响应的影响，结果如图4所示（为了简便，各图例中省去了数值的单位kN·m/rad）。
![]()
图 4 上、下挠性接头旋转刚度的不同组合对钻井管柱系统动态响应的影响Figure 4. Effect of different combinations of the rotational stiffness of upper and lower flexible joints on the dynamic response of drilling string system由图4（a）可知，上、下挠性接头的旋转刚度对管柱系统整体偏移的影响很小，可以忽略；但是，其不同组合对隔水管顶、底两端的受力状态有明显影响（见图4（b）和图4（c））；此外，随着上、下挠性接头的旋转刚度增大，其转角在较小范围内波动（见图4（h）和图4（i）），这有助于在钻井过程中将挠性接头的转角控制在允许的最大转角范围之内。因此，在设计深水钻井管柱系统时，需要将上、下挠性接头的旋转刚度组合跟前述主要设计参数一起优化。
4. 基于可靠度分析的多目标优化设计
根据参数敏感性分析结果，探讨隔水管壁厚、顶部偏移量、钻柱悬重、顶张力及上、下挠性接头旋转刚度的优化设计问题。
4.1 优化设计方法
根据深水钻井工程实践，进行管柱设计时应考虑的主要约束条件有：1）管柱的最大Mises应力小于其材料屈服应力的67%，如X80钢屈服应力的67%为370 MPa；2）管柱的最大偏移量小于海水深度的2%，如2 000 m水深下管柱的偏移极限为40 m；3）上挠性接头的最大转角小于5°（即0.087 rad），下挠性接头的最大转角小于2°（即0.035 rad）。
为了提高效率，采用Isight优化分析软件进行多目标优化设计。Isight软件集成了多项成熟的优化分析算法，并且可以将在Abaqus中建立的有限元模型直接嵌入其优化分析的流程之中，并通过识别Abaqus模型中的设计参数和动态响应设置参数范围和动态响应约束条件。基于此，建立了基于可靠度分析的深水钻井管柱系统多目标优化设计方法，设计流程如图5所示。
![]()
图 5 基于可靠度分析的多目标优化设计流程Figure 5. Workflow for the multi-objective optimization design based on reliability analysis该设计方法的设计思路是：1）选取邻域培植遗传算法（NCGA）作为多目标优化算法；2）海洋深水环境具有极不稳定的气候特性，需要分析由于系统参数波动所引起的优化方案的可靠性降低问题，因此选取Monte Carlo算法进行可靠度分析；3）为了将多目标优化设计和可靠度分析贯穿起来，采用可靠度优化算法Six Sigma进行算法整合。具体设计步骤为：用Six Sigma法驱动NCGA算法进行多目标优化设计，并将得到的非劣解传递给Monte Carlo算法；然后，用Monte Carlo算法计算该非劣解在参数波动范围内的可靠度，如果该可靠度满足Six Sigma算法设定的最低可靠度要求，则输出该非劣解；反之，则将可靠度不满足要求的情况反馈给NCGA算法，以帮助它调整并优化策略。
4.2 可靠度验证
采用图5所示流程进行了优化设计测试，设计参数取值：隔水管壁厚9.0～28.0 mm，顶部偏移量－20～20 m，隔水管顶张力400～1 500 kN，钻柱悬重600～1 000 kN，上挠性接头旋转刚度100～7 000 kN·m/rad，下挠性接头旋转刚度100～10 000 kN·m/rad。动态响应约束条件：最大偏移量40 m，最大Mises应力370 MPa，上挠性接头最大转角5.0°，下挠性接头最大转角2.0°。
采用Six Sigma法进行可靠度优化运算，将各项动态响应的可靠度最低要求都设为98%，得到的计算结果见表1（不收敛或违反约束的算例标为红色，寻优过程中间解标为白色，非劣解标为蓝色，推荐的最优解标为绿色）。
表 1 基于可靠度分析的多目标优化设计结果Table 1. Results of the multi-objective optimization design based on reliability analysis算例
序号旋转刚度/(kN·m·rad–1) 隔水管
壁厚/mm顶部
偏移/m隔水管
顶张力/kN钻柱
悬重/kN最大
偏移量/m最大Mises
应力/MPa最大转角/rad 下挠性接头 上挠性接头 下挠性接头 上挠性接头 1 5 861. 9 3 005.0 18.8 –8.18 930.1 977.9 15.40 242.4 –0.011 3 0.077 4 2 4 708.8 214.0 16.9 6.01 848.0 914.6 27.94 236.9 –0.024 8 0.090 1 3 4 383. 4 3 542.9 25.0 –12.54 1 391.7 960.8 5.90 226.6 –0.003 1 0.055 4 4 7 194.9 5 758.1 11.2 18.56 1 059.6 805.1 32.06 199.4 –0.018 2 0.059 5 5 2 360.0 1 198.2 10.6 6.25 1 339.5 850.5 21.27 212.5 –0.017 1 0.069 5 6 541.6 1 195.5 23.7 9.35 486.7 966.4 34.55 251.8 –0.038 6 0.096 4 7 10 009.6 4 037.1 15.1 12.03 1 109.4 798.1 27.70 198.4 –0.012 8 0.065 6 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 887 10 314.4 6 971.1 25.1 –16.75 1 468.0 735.7 2.48 163.1 0.009 7 0.048 3 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 995 2 170.9 6 401.2 25.1 –17.35 1 468.4 917.9 0.75 207.8 0.000 4 0.044 1 996 735.7 6 407.1 25.1 –16.76 1 468.7 918.9 1.15 208.0 –0.001 7 0.043 9 997 7 560.4 7 063.5 9.7 –16.75 1 471.6 983.8 0.05 227.9 0.011 5 0.043 4 998 9 299.3 6 880.8 9.7 –16.33 1 474.8 834.2 1.69 191.8 0.012 1 0.049 7 999 7 446.5 6 401.2 27.2 –17.35 1 468.4 754.2 2.18 167.7 0.005 5 0.049 5 1 000 9 375.0 6 554.6 9.7 –16.75 1 471.6 983.8 0.25 229.1 0.013 3 0.045 4 为了验证表1中各个非劣解在设计参数波动情况下的可靠性，对推荐的最优解（第887号算例的解）进行了Monte Carlo可靠度分析，允许各项设计参数在±10%范围内扰动，并将样本点设置为正态分布，最终得到的可靠度分析结果如图6所示。
![]()
图 6 最优解在设计参数±10%范围内波动时的可靠度分析结果Figure 6. Reliability analysis results of the optimal solution when fluctuations occur within ±10% of the design parameters由图6可以看出，最大偏移量的变化范围为8.35～16.85 m，最大Mises应力的变化范围为148～240 MPa，上挠性接头转角的变化范围为3.14°～4.50°（即0.054 8～0.078 5 rad），下挠性接头转角的变化范围为–0.43°～0.10°（即–0.007 5～0.001 7 rad）。与动态响应约束条件相比，各项动态响应值均符合设置的约束条件，即设计参数在±10%范围内扰动情况下，该设计的可靠度仍达到了100%。由此可知，设计的该参数组合满足要求。
进一步研究可知，多目标优化设计方法可以避免由于设计结果靠近约束边界而在参数波动情况下可靠度达不到要求的问题，因此所得到的非劣解集能够同时满足工程可行性和安全可靠性2方面的要求。
5. 结论与建议
1）相比于仅考虑隔水管的模型，根据管中管模型计算得到的管柱整体弯曲程度更小，即内外管柱的相互作用对管柱系统的整体偏移有抑制作用。这说明对于实际存在管中管结构的深水钻井管柱系统而言，在给定的管柱安全偏移范围内能够承受更大的海洋环境载荷。
2）建立了基于可靠度分析的多目标优化设计方法，采用该设计方法得到的非劣解集能够同时满足工程可行性和安全可靠性2方面的要求。
3）建立管中管模型时，没有考虑钻井液的影响，而钻井液可以等效为一个贯穿于管中管模型的弹塑性减振装置，其会在一定程度上减轻内外管柱之间的相互作用，从而影响管柱系统的整体动态响应，因此需要在接下来的研究中予以考虑。
-
![]()
图 3 主要设计参数对钻井管柱系统动态响应的影响
Figure 3. Effect of main design parameters on the dynamic response of drilling string system
![]()
图 4 上、下挠性接头旋转刚度的不同组合对钻井管柱系统动态响应的影响
Figure 4. Effect of different combinations of the rotational stiffness of upper and lower flexible joints on the dynamic response of drilling string system
![]()
图 5 基于可靠度分析的多目标优化设计流程
Figure 5. Workflow for the multi-objective optimization design based on reliability analysis
![]()
图 6 最优解在设计参数±10%范围内波动时的可靠度分析结果
Figure 6. Reliability analysis results of the optimal solution when fluctuations occur within ±10% of the design parameters
表 1 基于可靠度分析的多目标优化设计结果
Table 1 Results of the multi-objective optimization design based on reliability analysis
算例
序号旋转刚度/(kN·m·rad–1) 隔水管
壁厚/mm顶部
偏移/m隔水管
顶张力/kN钻柱
悬重/kN最大
偏移量/m最大Mises
应力/MPa最大转角/rad 下挠性接头 上挠性接头 下挠性接头 上挠性接头 1 5 861. 9 3 005.0 18.8 –8.18 930.1 977.9 15.40 242.4 –0.011 3 0.077 4 2 4 708.8 214.0 16.9 6.01 848.0 914.6 27.94 236.9 –0.024 8 0.090 1 3 4 383. 4 3 542.9 25.0 –12.54 1 391.7 960.8 5.90 226.6 –0.003 1 0.055 4 4 7 194.9 5 758.1 11.2 18.56 1 059.6 805.1 32.06 199.4 –0.018 2 0.059 5 5 2 360.0 1 198.2 10.6 6.25 1 339.5 850.5 21.27 212.5 –0.017 1 0.069 5 6 541.6 1 195.5 23.7 9.35 486.7 966.4 34.55 251.8 –0.038 6 0.096 4 7 10 009.6 4 037.1 15.1 12.03 1 109.4 798.1 27.70 198.4 –0.012 8 0.065 6 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 887 10 314.4 6 971.1 25.1 –16.75 1 468.0 735.7 2.48 163.1 0.009 7 0.048 3 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 995 2 170.9 6 401.2 25.1 –17.35 1 468.4 917.9 0.75 207.8 0.000 4 0.044 1 996 735.7 6 407.1 25.1 –16.76 1 468.7 918.9 1.15 208.0 –0.001 7 0.043 9 997 7 560.4 7 063.5 9.7 –16.75 1 471.6 983.8 0.05 227.9 0.011 5 0.043 4 998 9 299.3 6 880.8 9.7 –16.33 1 474.8 834.2 1.69 191.8 0.012 1 0.049 7 999 7 446.5 6 401.2 27.2 –17.35 1 468.4 754.2 2.18 167.7 0.005 5 0.049 5 1 000 9 375.0 6 554.6 9.7 –16.75 1 471.6 983.8 0.25 229.1 0.013 3 0.045 4 
下载: 导出CSV
-
[1] 周守为, 刘清友．深水钻井隔水管系统力学行为理论及应用研究[M]．北京: 科学出版社, 2016: 1-8． ZHOU Shouwei, LIU Qinyou. Theory and application research of deepwater drilling riser system mechanical behavior[M]. Beijing: Science Press, 2016: 1-8.
[2] BURKE B G. An analysis of marine risers for deep water[J]. Journal of Petroleum Technology, 1974, 26(4): 455–465.
[3] SIMMONDS D G. Dynamic analysis of the marine riser[R]. SPE 9735, 1980.
[4] KROLIKOWSKI L P, GAY T A. An improved linearization technique for frequency domain riser analysis[R]. OTC 3777, 1980.
[5] DAREING D W, HUANG T. Natural frequencies of marine drilling risers[J]. Journal of Petroleum Technology, 1976, 28(7): 813–818.
[6] KIRK C L, ETOK E U, COOPER M T. Dynamic and static analysis of a marine riser[J]. Applied Ocean Research, 1979, 1(3): 125–135.
[7] PATEL M H, SAROHIA S, NG K F. Finite-element analysis of the marine riser[J]. Engineering Structures, 1984, 6(3): 175–184.
[8] 王宴滨, 高德利, 房军. 海洋钻井隔水管–钻井液横向耦合振动特性[J]. 石油钻采工艺, 2015, 13(1): 25–29. WANG Yanbin, GAO Deli, FANG Jun. Characteristic analysis on transverse coupling vibration of offshore drilling riser-drilling fluid[J]. Oil Drilling & Production Technology, 2015, 13(1): 25–29.
[9] 闫铁, 刘珊珊, 毕雪亮, 等. 考虑钻柱作用的深海隔水管柱力学分析[J]. 石油矿场机械, 2014, 43(8): 13–17. YAN Tie, LIU Shanshan, BI Xueliang, et al. Mechanical analysis of deep sea riser considering the effects of drill string[J]. Oil Field Equipment, 2014, 43(8): 13–17.
[10] YAZDCHI M, CRISFELD M A. Non-linear dynamic behaviour of flexible marine pipes and risers[J]. International Journal for Numerical Methods in Engineering, 2002, 54(9): 1265–1308.
[11] RAMAN-NAIR W R, BADDOUR R E. Three-dimensional dynamics of a flexible marine riser undergoing large elastic deformations[J]. Multibody System Dynamics, 2003, 10(4): 393–423.
-
期刊类型引用(13)
1. 马永乾，赵鹏，畅元江，王仕超，张晏铭. 基于LAGRANGE方程的深水钻井隔水管–水下井口系统动力分析. 应用科技. 2024(01): 151-157 . 
百度学术
2. 王文昌，徐祖凯，周星，王昭彬，陈锋. 超深井钻柱动态疲劳失效特征及参数优选. 石油钻探技术. 2024(02): 118-125 . 本站查看
3. 畅元江，王仕超，李健，刘秀全. 深水隔水管-防喷器组耦合系统动力学建模与仿真验证. 振动与冲击. 2024(11): 1-8 . 百度学术
4. 狄勤丰，骆大坤，秦垦，王文昌，陈锋. 考虑“筒中筒”耦合作用的深水隔水管柱与钻柱碰摩运动规律. 石油钻采工艺. 2024(01): 25-32 . 百度学术
5. 王高伟，廖茂林，蒋宏伟，郜志英，房超. 基于钻井液上返流动作用下的隔水管动力学行为研究. 船舶力学. 2023(05): 719-730 . 百度学术
6. 孙巧雷，张辉易，姚力萍，李乐勤，张圆圆，冯定. 海上管中管结构力学研究现状与发展. 科学技术与工程. 2023(21): 8907-8915 . 百度学术
7. 曹峰，胡波，王安伟，刘永亮. 压裂泵的数值分析及仿真标定验证. 机械设计与研究. 2023(06): 205-210 . 百度学术
8. 李旭，朱志强，冯颖韬，赵琥，温达洋，山永林，林黎明. 海洋钻井隔水管系统动力学行为研究. 钻采工艺. 2022(03): 31-36 . 百度学术
9. 张昊，雷雨辰，高云飞，张维维，闫启方，刘林超. 分数黏弹性可压缩油气井管杆准静态响应分析. 信阳师范学院学报(自然科学版). 2021(02): 331-338+344 . 百度学术
10. 李泉新，褚志伟，许超，杨冬冬. 煤矿井下双动力复合定向钻进轨迹调控研究. 工矿自动化. 2021(12): 25-31 . 百度学术
11. 王江帅，李军，柳贡慧，罗晓坤. 气侵条件下新型双梯度钻井环空出口流量变化规律研究. 石油钻探技术. 2020(04): 43-49 . 本站查看
12. 张光伟，乔阳，田帆，高嗣土，尹福来. 全旋转内置式导向钻井工具动力学分析及仿真. 机械强度. 2020(06): 1417-1423 . 百度学术
13. 董学莲，杨芯萍，赵雄，张聪林，郑亮，费凌. 液缸式隔水管张紧系统建模及性能分析. 液压与气动. 2019(08): 39-48 . 百度学术
其他类型引用(9)
计量
- 文章访问数: 1260
- HTML全文浏览量: 753
- PDF下载量: 74
- 被引次数: 22
