海洋油田水下集输-立管内两相流冲击力预测模型

2020-05-08李文升冯耀荣张西宁朱文军王广明

装备环境工程 2020年4期

李文升，冯耀荣，张西宁，朱文军，王广明

（1. 中国石油集团石油管工程技术研究院，西安 710077；2. 西安交通大学，西安 710049；3. 渤海装备（天津）新世纪机械制造有限公司抽油机制造厂，天津 300280）

海洋油田水下集输立管系统内两相流引发的压力波动会引发立管系统的振动，而立管是整个海底管道中约束最弱的部位。管道在内外流载荷的交互作用下的振动，严重威胁了立管系统的安全运行。当管内产生严重段塞流这类周期性作用力时，其频率若接近立管的固有频率，将引发流固耦合振动，加速立管的破坏，并直接威胁到中心平台或海上油轮的安全。另外，管道内壁所受冲击力直接影响到了其腐蚀速率，特别是管道弯头处。因此分析研究管内两相流对管道冲击力的特性，对于保障海洋油气的安全生产具有重要意义[1-3]。

国内外目前对严重段塞流引起的流固耦合振动研究相对较少，Fylling等[4]采用准稳态模型研究了段塞流经过悬链管弯管时引起的支撑力变化。Patel与Seyed[5]发现柔性管内段塞流可以引起轴向张力的变化，从而引发管道破坏。Valdivia[6]在2008年采用小型试验装置分析了管道内流对不锈钢悬链管动力学特性的影响。Ita[7]对塑料柔性管在严重段塞流工况下的位移特性进行了分析。Ortega等[8]基于线性有限元方法，构建了柔性立管在不同流态下的非稳态结构动力学预测模型，其研究发现，管内段塞流的出现明显影响了柔性立管的位移及应力变化[9-10]。而国内则鲜有严重段塞流对管道内壁冲击的相关研究。

文中针对严重段塞流这一典型不稳定流动，构建了集输立管系统的气液两相流冲击力预测模型，分析了管内段塞频率、管外涡街频率等流动特性与管道固有频率的关系，得到了管道不同部位所受的冲击力频率、大小及振幅特性。所分析的实验管道内径为50.8 mm、水平段管道长度为114 m、下倾段为19 m（倾角为-2°）、垂直立管高度为15.3 m的集输-立管系统[5]。

1 频率分析

在海底输油系统中，当立管道发生严重段塞流时，有一定的段塞频率，而外流绕过立管时，存在着漩涡释放频率。流体（含内流、外流）的流动是激发立管振动的动力之一，当流体流动产生的激振力频率与立管道固有频率相近时，将有可能发生立管系统流固耦合振动的现象，从而给生产安全带来危害。因此，必须对立管固有频率、涡街频率、严重段塞流频率进行分析，确定是否产生流固耦合。

1.1 立管固有频率

管线的固有频率，取决于圆管的特性和管跨两端的支撑条件（约束），以及管跨长度和管线单位长度的质量（包括附加质量）。以文中管道为例，根据结构力学，可得出管跨的一阶固有频率（即基频）为：

式中：C为系数，在真空中取1，在水中取0.7；K为管跨约束相关的系数，两端固支取4.73，两端铰支取3.14，一固一铰也可取3.14；E管材的弹性模量，取2.06×1011N/m2；I为管界面的惯性矩，可由管径得到，为4.4×10-5N/m2；L为管跨长，L=133 m；Mp为管线单位长度质量，含管材、管内流体、管外附加流体质量（即排出水体积的 1～2倍来考虑，文中为554.1 g）。由以上计算得出，管材的固有频率 fp=1.6×10-2Hz。

1.2 涡街释放频率

外流绕过立管道会释放出漩涡，其释放频率可采取式（2）估算：

式中：Sr为Strouhal数，它是雷诺数和管线界面形状所组成的数，可通过实验得出；u为海流速度；D为管线外径。若取海流速度u=0.4 m/s，D=0.0508 m，则通过查表得Sr≈0.21，代入到式（2）可得漩涡释放的频率fs=0.3 Hz。

1.3 严重段塞流频率分析

根据李文升[11]的实验结果，严重段塞流周期约为0.5～2 h，即频率约为10-4量级。由前述分析可知，漩涡释放频率约为10-1数量级，管道固有频率约为10-2数量级。因此，对于海底管道生产系统，当内流发生严重段塞流流动，且外流遭受涡激时，由于涡激力与严重段塞力的频率无法耦合在一起，一般不会发生同时诱发管道共振的可能。严重段塞流发生频率与管道固有频率相差较大，也不会因此而产生共振。

2 严重段塞流冲击力预测

2.1 控制方程

当立管产生严重段塞流时，在管道的不同界面上，管内流体速度、相含率、压力等参数均随时间而发生变化。在这些流体参数的共同作用下，管道上的受力也随之变化。从便于分析问题的角度出发，首先在任意一曲管单元来分析流体的作用力，在此基础上，再详细地分析严重段塞流对立管结构冲击力作用。管道单元结构如图1所示。

图1 管道单元Fig.1 Pipe unit

图 1 为任意一段曲管单元，内充有流体，其流速在两截面处的速度分别为 u1和 u2。采用固结于管道上的绝对坐标系，并取弯管的体积为控制体，以A0表示控制体的侧面，A1、A2表示其截面面积，n1、n2分别为截面外法线单位向量。若以F表示流体作用于管道壁的合力，则有：

式中：Pn为应力张量法向分量。由动量守恒方程可得：

将式（4）代入式（3）中，可以得到管道控制体A0所受合力的表达式为：

假设进出口截面上的速度均匀分布，即其切应力为0，在截面上只有压力P，同时在A0面上有u⊥n，因此得：

从式（8）可知，若已知流体参数在各截面的分布，则可通过式（6）得到流体对曲管单元的作用力，同时可以看出影响流体对结构作用力的因素有两类：流体参数和结构参数。其中流体参数有流体密度ρ、流体速度u和压力P；而结构参数有控制体体积V、截面面积A和形状参数n。式（8）中的6项可以分成重力动量变化所引起的力和压力变化所引起的力（ - p1A1n1- p2A2n2）等三部分。

2.2 预测结果

根据上述冲击力计算模型，对实验管道冲击力进行数值模拟，计算不同气相、液相折算速度下严重段塞流对管道冲击力影响，计算点布置如图2所示。选取典型严重段塞流工况进行管壁冲击力分析[11]，气相折算速度为0.4 m/s、液相折算速度为0.4 m/s。管道位置1、2、3、4、5处所受的冲击力如图3所示，管内工质分别为水和空气，管道单元长度为1 m。

图2 计算点布置Fig.2 Layout of calculation point

图 3a—f分别为五个典型管单元 x、y方向所受段塞流作用力大小变化规律。在x方向，管道单元受到的冲击力最大部位是水平管道和下倾管道连接处（位置2），最大值达到703.3 N，振幅波动为266.2 N。y方向管道单元受段塞流冲击力最大部位在下倾管底部与立管连接处（位置4），最大达到-993.5 N，振幅波动为450.7 N。值得注意的是，立管管道单元5受到最大冲击力为80.1 N，振幅波动为88.7 N。其余部分冲击力较小，见表1。

由于1、2、3、5四个管单元的长度一样，即1 m，而管道的截面面积沿线都不变，因此四个管单元的体积是相同的。对于结构参数来说，形状参数n是四个管单元唯一不同之处，它决定了流体的流动方向。因此，流动方向发生变化的部位，流体作用力很大。另一方面，流体参数的变化，也使得冲击力的幅值发生变化，结合这两方面的因素，在水平管和下倾管、下倾管和立管连接处，由于流动方向发生变化，力的大小和变化幅度都比其他位置的管道单元要大很多。在流动方向不发生改变的管单元作用力（如立管段和倾斜管段），完全取决于流体参数的变化，而流体参数的改变影响动量变化，如密度、速度等参数。由图3压力变化曲线对比可知，在液塞生长和液体排出过程中，各管道中冲击力基本保持恒定。当气体喷发时，各管道处压力、流动速度和密度快速变化，从而导致各管道处冲击力剧烈变化。

图3 不同管道节点位置的冲击力预测结果Fig.3 Prediction results of impact at different pipeline node positions: a) position 1; b) position 2; c) position 3;d) position 4; e) position 5; f) severe slug flow pressure fluctuations

表1 计算节点位置幅值统计Tab.1 Compute node location amplitude statistics N

需要注意的是，文中计算单元的长度为1 m，未充分考虑局部湍流等微观现象引起的流动对管壁冲击力的影响，但提供了一种快速评估严重段塞流工况对管壁冲击力特性的分析方法，有助于对严重段塞流引起的流固耦合振动现象进行深层次的模拟计算，同时对工程现场中水下柔性立管的设计也具有一定的借鉴意义。后续相关研究可将严重段塞流冲击力计算方法与瞬态双流体模型相结合，从而得到整个集输-立管系统内严重段塞流对壁面冲击力的分布特性。