一种压力管道泄露过程的仿真建模方法技术

本发明专利技术公开了一种压力管道泄露过程的仿真建模方法，其步骤如下：S1：建立流体性质的计算方法；S2：将压力管道视为压力容器，建立相应泄露过程的守恒方程；S3：建立管道初始条件的计算方法；S4：建立管道泄露处边界条件的计算方法；S5：将步骤S1

【技术实现步骤摘要】
一种压力管道泄露过程的仿真建模方法


[0001]本专利技术涉及模拟仿真领域，尤其涉及一种压力管道泄露过程的仿真建模方法。

技术介绍

[0002]压力管道被广泛地应用于运输工业流体，包括石油、天然气、二氧化碳、碳氢化合物等，所涉及的高压，加上流体易燃、易爆或者有毒的性质，一旦发生管道泄露事故，会对生命、财产和环境构成重大风险。因此，人们必须对压力管道进行相应的风险评估，从而制定相应的事故应急预案，提高其安全性能，保证其安全运行。上述的核心是对压力管道的泄露过程进行仿真模拟，对过程中所涉及的流体性质（例如，流量、流速、压力、温度、相态等）随时间的变化进行准确地预测。目前，用于上述目的的预测模型主要为基于求解偏微分方程的数值模型，此类模型的优点是数学描述符合实际物理过程，计算精度高，适用性高，缺点是计算资源消耗大，模拟耗时长（尤其是在模拟长距离运输管道泄露或者多根管道同时泄露的情况下），严重影响进行风险评估以及制定相应事故应急预案的效率。
[0003]针对以上，我们提出一种压力管道泄露过程的仿真建模方法对上述缺点进行改进。

技术实现思路

[0004]本专利技术旨在解决现有技术中存在的关键问题，特别创新地提出了一种压力管道泄露过程的仿真建模方法。
[0005]本专利技术的技术方案如下：
[0006]一种压力管道泄露过程的仿真建模方法，包括以下建模步骤：
[0007]S1：建立流体性质的计算方法；
[0008]S2：将压力管道视为压力容器，建立相应泄露过程的守恒方程；r/>[0009]S3：建立管道初始条件的计算方法；
[0010]S4：建立管道泄露处边界条件的计算方法；
[0011]S5：将步骤S1
‑
S4中所涉及的计算合理串联，建立算法用以模拟实际的泄露过程。
[0012]作为本专利技术的进一步方案，所述步骤S1具体涉及：首先判断管道内流体的相态，再采用合适的状态方程对其性质（主要包括压力、温度、密度、焓、熵、声速、相态等）进行计算。对于单相流，流体性质由状态方程直接求得。对于多相流，本专利技术采用经过大量实验数据验证过的均相流模型来处理其性质，即流体中每个相态的压力、温度和速度分别相等；多相流的其他性质则根据以下公式计算：
[0013] (1)
[0014]其中，指代多相流的性质，指代气相的性质，指代液相的性质，为蒸气量。
[0015]作为本专利技术的再进一步方案，所述步骤S2具体涉及将压力管道假设为压力容器，并根据管道流体的特征，对压力容器泄漏模型的守恒方程进行改写，包括质量和能量守恒
方程。改写后的质量和能量守恒方程分别为：。
[0016] (2)
[0017] (3)
[0018]其中，为管道内流体的质量，为管道内流体的能量，为时间，为管道上游的质量流量，为管道泄漏处的质量流量，为管道上游流体的总焓，为管道泄漏点流体的总焓，为管道与环境的热交换。
[0019]所得守恒方程(1)和(2)即作为描述压力管道泄露过程的全部守恒方程。
[0020]进一步地，管道上游和泄漏处的质量流量和分别由以下公式计算：
[0021] (4)
[0022] (5)
[0023]其中，为管道上游的流体密度，为管道上游的流体流速，为管道的内横截面积，为泄露系数，为管道泄露口的流体密度，为管道泄露口的流体流速，为管道泄露口的面积。
[0024]管道上游和泄漏口流体的总焓和分别由以下公式计算：
[0025] (6)
[0026] (7)
[0027]其中，为管道上游流体的焓，为泄漏口流体的焓。
[0028]作为本专利技术的再进一步方案，所述步骤S3具体涉及对管道发生泄漏前的管道压降进行计算，并由此计算出流体的其他性质。本专利技术具体采用以下方法来计算管道压降：
[0029]假设管道在发生泄漏前，压力管道内流体处于恒温稳态，其对应的压降计算公式为：
[0030] (8)
[0031]其中，为压力，为压缩因子，为气体常数，为温度，为速度，为距离，下标1和2分别代表管道上不同的位置，代表管道内壁与流体的摩擦力项。
[0032]作为本专利技术的再进一步方案，所述步骤S4具体涉及设置相应的边界条件对管道泄漏口的非稳态流动进行计算。具体方法为：
[0033]针对管道泄露处建立相应的能量守恒方程：
[0034] (9)
[0035]其中，为泄漏口上游处流体的焓，为泄漏口上游处流体的速度。
[0036]基于求解公式(9)，计算出泄漏口的流体性质。
[0037]作为本专利技术的再进一步方案，所述步骤S5具体涉及将步骤S1
‑
S4中所涉及的计算合理串联，建立相应的循环迭代算法来模拟压力管道泄漏在时间上的推进过程，从而得到管道内和泄漏口的流体性质随时间的变化。
[0038]与现有技术相比，本专利技术的有益效果是：
Peng
‑
Robinson状态方程（现有技术）。
[0056]所述步骤S2具体涉及将压力管道假设为压力容器，并根据管道流体的特征，对压力容器泄漏模型的守恒方程进行改写，包括质量和能量守恒方程。改写后的质量和能量守恒方程分别为：
[0057] (2)
[0058] (3)
[0059]其中，为管道内流体的质量，为管道内流体的能量，为时间，为管道上游的质量流量，为管道泄漏处的质量流量，为管道上游流体的总焓，为管道泄漏点流体的总焓，为管道与环境的热交换。
[0060]所得守恒方程(1)和(2)即作为描述压力管道泄露过程的全部守恒方程。
[0061]进一步地，管道上游和泄漏处的质量流量和分别由以下公式计算：。
[0062] (4)
[0063] (5)
[0064]其中，
[0065]为管道上游的流体密度，可由步骤S1所述方法求得；
[0066]为管道上游的流体流速，其计算方法取决于管道上游的类型。根据工程实际，压力管道的上游可有多种类型，例如，泵、压缩机等工程设备或者油气藏等天然物质聚集。针对不同的类型，管道上游的流速可有不同的计算方法，例如，若管道上游为泵，通常为时间的函数，即；再例如，若管道为隔离管道，则管道内无物质流动，即等于0。
[0067]为管道的内横截面积，由管道的内径决定；
[0068]为泄露系数，取值为0到1。此系数用以量化流体通过泄露口时的摩擦和传热效应对流量的影响，其通常由泄露口的大小决定。对于小孔，摩擦和传热效应对流量的影响通常较大，取值一般较小，对于大孔，摩擦和传热效应对流量的影响通常较小，取值一般较大；例如，对于压力管道的全口径破裂，由于泄露口面积为理论上最大值，摩擦和传热效应对流量的影响微乎其微，故一般取1。
[0069]为管道泄露口的流体密度，可由步骤S1所述方法求得；
[0070]为管道泄露口的流体流速，可由步骤S4所述方法求得；<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种压力管道泄露过程的仿真建模方法，其特征在于：其建模步骤如下：S1：建立流体性质的计算方法；S2：将压力管道视为压力容器，建立相应泄露过程的守恒方程；S3：建立管道初始条件的计算方法；S4：建立管道泄露处边界条件的计算方法；S5：将步骤S1
‑
S4中所涉及的计算合理串联，建立算法用以模拟实际的泄露过程。2.根据权利要求1所述的压力管道泄露仿真建模方法，其特征在于：所述步骤S1具体涉及：首先判断管道内流体的相态，再采用合适的状态方程对其性质（主要包括压力、温度、密度、焓、熵、声速、相态等）进行计算；对于单相流，流体性质由状态方程直接求得；对于多相流，本发明采用经过大量实验数据验证过的均相流模型来处理其性质，即流体中每个相态的压力、温度和速度分别相等；多相流的其他性质则根据以下公式计算： (1)其中，指代多相流的性质，指代气相的性质，指代液相的性质，为蒸气量。3.根据权利要求1所述的压力管道泄露仿真建模方法，其特征在于：所述步骤S2具体涉及将压力管道假设为压力容器，并根据管道流体的特征，对压力容器泄漏模型的守恒方程进行改写，包括质量和能量守恒方程；改写后的质量和能量守恒方程分别为： (2) (3)其中，为管道内流体的质量，为管道内流体的能量，为时间，为管道上游的质量流量，为管道泄漏处的质量流量，为管道上游流体的总焓，为管道泄漏点流体的总焓，为管道与环境的热交换；所得守恒方程(1)和(2)即作为描述压力管道泄露过程的全部守恒方程；进一步地...

【专利技术属性】
技术研发人员：易家欢，
申请(专利权)人：易家欢，
类型：发明
国别省市：