一种油气输送管道管内流体密度随机变化的模拟方法技术

本发明专利技术公开了一种油气输送管道管内流体密度随机变化的模拟方法，包括：根据管道内输送石油和天然气的相关随机参数，构建管内流体密度随机变化模型；通过气液两相流物理模型实验，获得所述相关随机参数的概率分布，并进行随机取值；根据所述管内流体密度随机变化模型以及所述相关随机参数的随机取值，耦合柔性管道结构振动与密度随机变化流体的动力控制方程；采用有限元法或有限差分法，对所述动力控制方程进行求解，求得管道结构振动的动力响应结果及结构内部应力应变的时域响应结果。该模拟方法使得管内气液两相流随机性分析结果更贴近实际情况，提高了海洋管道结构疲劳使用寿命预报的可靠性和有效性。命预报的可靠性和有效性。命预报的可靠性和有效性。

【技术实现步骤摘要】
一种油气输送管道管内流体密度随机变化的模拟方法


[0001]本专利技术涉及油气输送
，特别涉及一种油气输送管道管内流体密度随机变化的模拟方法。

技术介绍

[0002]地球上71％的面积为海洋，在海洋的底部蕴藏着丰富的石油和天然气。随着陆地上油气资源的开采逐渐枯竭，产量日益减少，人们不断加大了对海洋油气资源的开采。在海洋油气资源的开采中，管道常用于输送石油和天然气，因而得到了广泛的应用，比如海底管线、顶张力立管、钢悬链线立管和缓波型立管等。
[0003]由于海底地质、化学反应和历史变迁等的不同，在不同海域，海底油气资源的储存可能不同。有些海域储存着大量的天然气，有些海域储存着大量的石油，而有些海域同时储存着天然气和石油。这些油气资源一般采用管道进行输送，管道输送具有效率高、连续性好、输送量大和成本较低等优点。管道对海底产出的天然气进行输送时，由于周围环境温度、压强等的变化，某些重烃将会液化为液体，此时管道内部为气体和液体组成的气液两相流。在另外一方面，管道对海底产出的石油进行输送时，由于周围环境温度、压强等的变化，某些轻烃将会气化为气体，管道内部为气液两相流体。对于同时产出石油和天然气的油气井，管道内部通常输送着气体和液体。
[0004]由此可见，由气体和液体组成的气液两相流在海洋油气资源输送管道中较为常见。气体和液体在管道内部进行流动时，容易发生变形、分离和混合，管道内部可具有多种流型，比如气泡流、段塞流、泡沫流、环形雾流和雾状流等，而且这些流型极不稳定，容易受到外界干扰而发生变化。管道内部气液两相流流体的总密度、流速和压强等参数将会发生随机性的变化。海洋油气输送管道长达数百米，甚至上千米，具有较强的柔性。在内部流体随机性激励下将会发生振动，管道结构内部的应力应变将会循环变化。在长时间激励下，管道结构容易出现疲劳裂纹，甚至是疲劳破坏。一旦管道发生疲劳破坏，管内石油和天然气就会泄露，从而造成重大的安全事故、严重的环境污染和巨额的经济损失。
[0005]现有技术对于管内石油和天然气组成的气液两相流做了大量的简化，一般采用单相流模型描述管内流体的流动，亦或者采用稳定的段塞流模型模拟气液段塞流对管道的作用，并未考虑管道内部气体和液体混合流动过程中的随机性变化、管内流体总密度的随机性变化、流体密度随机性变化与管道结构振动的耦合作用等。

技术实现思路

[0006]本专利技术的目的在于克服现有技术的缺陷，并为此提供一种油气输送管道管内流体密度随机变化的模拟方法，该方法考虑了天然气和石油在管道内部混合流动时流体总密度随机性变化的幅值、随机性变化的频率、随机性变化的流速和随机性变化的初相位，并将密度随机性变化流体的流动与管道结构的运动进行耦合，从而改进了现有技术对管内气液两相流随机性的模拟方法，使得管内气液两相流随机性分析结果更贴近实际情况，提高了海
洋管道结构疲劳使用寿命预报的可靠性和有效性。
[0007]为实现上述目的，本专利技术采取的技术方案为：
[0008]本专利技术提供一种油气输送管道管内流体密度随机变化的模拟方法，包括以下步骤：
[0009]根据管道内输送石油和天然气的相关随机参数，构建管内流体密度随机变化模型；
[0010]通过气液两相流物理模型实验，获得所述相关随机参数的概率分布，并进行随机取值；
[0011]根据所述管内流体密度随机变化模型以及所述相关随机参数的随机取值，耦合柔性管道结构振动与密度随机变化流体的动力控制方程；
[0012]采用有限元法或有限差分法，对所述动力控制方程进行求解，求得管道结构振动的动力响应结果及结构内部应力应变的时域响应结果。
[0013]进一步地，所述管道内输送石油和天然气的相关随机参数包括：
[0014]石油和天然气在管道内部混合流动时流体总密度随机性变化的幅值、随机性变化的频率、随机性变化的流速和随机性变化的初相位。
[0015]进一步地，所述管内流体密度随机变化模型为：
[0016][0017][0018]式中：ρ(s,t)表示管内流体的总密度，s表示拉格朗日坐标，为当前点沿着管道到原点的距离，t表示时间；
[0019]表示管内流体的平均密度；
[0020]ε
j
表示第j项管内流体密度变化的幅值；
[0021]exp[]表示以自然指数e为底的指数函数；
[0022]i表示虚数单位；
[0023]k
j
表示第j项管内流体密度变化的波速，其值与圆频率ω
j
和流速u
j
有关；
[0024]ω
j
表示第j项管内流体密度变化的圆频率；
[0025]θ
j
表示第j项管内流体密度变化的初相位；
[0026]u
j
表示第j项管内流体密度变化的流速。
[0027]进一步地，通过气液两相流物理模型实验，获得所述相关随机参数的概率分布，并进行随机取值；包括：
[0028]通过气液两相流物理模型实验，测量统计石油和天然气在管道内部混合流动时管内流体密度随机性变化的幅值、管内流体随机性变化的频率、管内流体随机性变化的流速和随机性变化的初相位；
[0029]根据测量统计的概率分布，进行随机取值。
[0030]进一步地，耦合柔性管道结构振动与密度随机变化流体的动力控制方程；包括：
[0031][0032]式中：m
s
表示管道结构单位长度的质量；
[0033]y表示管道结构振动的横向位移；
[0034]c
s
表示管道结构振动的阻尼系数；
[0035]A
in
表示管道结构内部横截面的面积；
[0036]T表示管道结构的轴向力；
[0037]EI表示管道结构的弯曲刚度。
[0038]与现有技术相比，本专利技术具有如下有益效果：
[0039]该方法考虑了海洋管道内部天然气和石油混合输送时，气液两相流流体总密度随机变化的幅值、随机变化的频率、随机变化的流速和随机变化初相位等因素的影响，采用物理模型实验确定各随机参数的概率分布，耦合管道结构振动与管内流体运动，从而建立了一个新的柔性管道输送密度随机变化流体的模拟方法，为海洋油气输送管道的安全性、可靠性和耐久性分析提供了理论依据和途径。
附图说明
[0040]图1为本专利技术实施例提供的油气输送管道管内流体密度随机变化的模拟方法流程图；
[0041]图2为本专利技术实施例提供的油气输送管道管内流体密度随机变化的模拟方法原理图。
具体实施方式
[0042]为使本专利技术实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本专利技术。
[0043]在本专利技术的描述中，需要说明的是，术语“上”、“下”、“内”、“外”“前端”、“后端”、“两端”、“一端”、“另一端”等指示的方位或位置关系为基于附图所示的方本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种油气输送管道管内流体密度随机变化的模拟方法，其特征在于，包括以下步骤：根据管道内输送石油和天然气的相关随机参数，构建管内流体密度随机变化模型；通过气液两相流物理模型实验，获得所述相关随机参数的概率分布，并进行随机取值；根据所述管内流体密度随机变化模型以及所述相关随机参数的随机取值，耦合柔性管道结构振动与密度随机变化流体的动力控制方程；采用有限元法或有限差分法，对所述动力控制方程进行求解，求得管道结构振动的动力响应结果及结构内部应力应变的时域响应结果。2.根据权利要求1所述的一种油气输送管道管内流体密度随机变化的模拟方法，其特征在于，所述管道内输送石油和天然气的相关随机参数包括：石油和天然气在管道内部混合流动时流体总密度随机性变化的幅值、随机性变化的频率、随机性变化的流速和随机性变化的初相位。3.根据权利要求2所述的一种油气输送管道管内流体密度随机变化的模拟方法，其特征在于，所述管内流体密度随机变化模型为：征在于，所述管内流体密度随机变化模型为：式中：ρ(s,t)表示管内流体的总密度，s表示拉格朗日坐标，为当前点沿着管道到原点的距离，t表示时间；表示管内流体的平均密度；ε
j
表示第j项管内流体密度变化的幅值；exp[]表示以自然指数e为底的指数函数；i表示虚数单位；k
j
...

【专利技术属性】
技术研发人员：谢武德，梁振林，胡中元，潘加辉，
申请(专利权)人：山东大学，
类型：发明
国别省市：