爬坡管路石蜡沉积层最大剪应力的数值模拟.PDF

工程技术 “爬坡”管路石蜡沉积层最大剪应力的数值模拟 1 2 （1. 天津市大港油田集团建设监理有限责任公司；2. 中国石油大学（华东）储建学院） 毕露 赵舒宇 摘要：由于长输管路的距离很长，所经地区地形复杂，部分热含蜡原油管 化为沿轴向方向限制的固定边界。石蜡沉积物的内表面设置为流固 路不可避免的形成了落差大、翻越点多的恶劣状况，管路出现了“爬坡”现象。 耦合界面。固体结构的被表面和流体结构的外表面都设置为耦合界 “爬坡”热含蜡原油管道在正常运行过程中，管道内的石蜡沉积物在油流剪切 面，进而实现结构和流体的耦合。施加的载荷是石蜡沉积物自身的重 应力作用下处于受力平衡状态。但是在管道停输后，由于油流流动形成的剪 力。固体结构有限元模型如图1- 1。 力消失，管壁处的石蜡沉积物处于失稳状态。于是在重力和未凝原油的共同 作用下，石蜡沉积物会发生剪切变形、蠕变，结构破坏后向管道低洼处运移， 进而造成管道低洼处凝管堵塞，使整条管线面临启动困难和全线原油凝管的 危险。本文利用ADINA 的流固耦合功能，对停输后石蜡沉积物的力学状况进 行模拟，探讨了管径等因素对石蜡沉积物的最大剪应力的大小以及分布的影 [1] 响，为工程实际提供了参考 。 关键词： 管路 石蜡沉积层 最大剪应力 数值模拟 1 有限元模型 1.1 材料计算模型 图1-1 石蜡沉积物结构的有限元模型 在管道停输的瞬间石蜡沉积物并不会发生破坏滑脱流动。这种 液态原油部分模型上、下端面和石蜡沉积物相应端面的边界条 情况下，在力学作用下，石蜡沉积物首先发生瞬时弹性应变，紧接着 件设置一致，外表面设置为流固耦合界面，耦合方式为双向直接耦 发生粘弹性蠕变，最后是塑性断裂，间隔时间短暂。为达到分析管道 合，收敛准则为力和位移以以使结果更精确。这部分有限元模型的建 停输后等效应力的分布规律的目的，数值计算时，石蜡沉积物选取多 立如图1- 2。 [2] 线性弹塑性中的lubby2 算法 ，此方法广泛用于模拟岩石，混凝土材 料的长期蠕变行为，并考虑了卸载等因素的影响。中心未凝原油采取 粘弹性流体模型。对于lubby2 模型，其蠕变应变率表达式如下： 1- 1 其中蠕变参数与应力的关系如下： 1- 2 1- 3 1- 4 图1-2 液态原油部分有限元模型 对于上述方程， 是六个材料相关参数，在 2 计算结果以及应力分析 ADINA 软件输入界面中分别用a0、a1、a2、a3、a4、a5 表示，通过试 2.1 管径对沉积层应力分布的影响 验结果来拟合上述六个蠕变参数，进而得到准确的材料模型。参数由 在长距离含蜡原油输送管路中，管道的管径是相差很大的。因此 石蜡沉积物的直剪实验，可得相关参数如下： 讨论在不同管径下的剪应力分布很有必要。相同物性的石蜡沉积物 - 1 在不同的管径的管道里，沉积层的剪切应力的大小与分布也是不同 m=- 0.01995MPa ；G=0.04991 MPad - 1 的。以下是当沉积物厚度为0.11m ，爬坡角度为十度，管长为4m 时 k1=0.009034MPa ； =100757.588MPa - 1 刻不同管径的剪应力分布。 k2=0.01625 Mpa ； =337MPad 1.2 单元选择与网格划分