高压输气管线减压波分析预测方法技术

高压输气管线减压波分析预测方法，应用于高压天然气输送管道断裂控制和安全保障技术领域。在测得的管线某点开裂前初始压力Ｐ１、温度Ｔ１下，根据不同的气质选择相应的状态方程，即ＳＲＫ、ＰＲ、ＢＷＲＳ－ＰＲ、ＢＷＲＳ－ＳＲＫ四个状态方程及组合，计算等熵点压力温度下流体的密度、逸度、定容热容、定压热容、比热等流体特性，用于判断管线中流体所处相态，描述真实流体等熵过程的特性，准确得到减压波速。效果是：对重要输气管道的气体减压波曲线进行分析预测，进而预测管材的止裂韧性，分析预测结果，优化了管道设计和运行参数，使用方便。

【技术实现步骤摘要】

本专利技术涉及高压天然气输送管道断裂控制和安全保障
，提出了高压输气 管线减压波的分析和预测方法，开发含气液两相流的多组分天然气的减压波曲线计算分析 方法及软件，用于重大输气管道工程气体减压波分析预测及管材止裂韧性预测。
技术介绍
输气管道的安全运行一直是国内外油气行业关注的焦点。管道建设和运行过程中 由于制造缺陷、施工缺陷，及冲刷、腐蚀等可能产生各种缺陷，这些缺陷管道在运行过程中 一旦有强烈的波动或剧烈的外部干扰就有可能发生断裂或促进裂纹扩展。而高压输气管线 一旦发生断裂，往往引发不堪设想的重大事故。随着国内外输气管道逐渐向高压、高钢级管 道方向发展，对管道安全性提出了更新、更高的要求。因此输气管道止裂韧性预测及止裂控 制具有重要意义。目前，Battelle的双曲线法是一种常用的管道止裂韧性预测方法，也适用于重大 输气管道工程如西气东输二线的止裂韧性预测。钢管止裂韧性预测需要知道输送气体的减 压波曲线和材料的阻力曲线。其中气体减压波曲线是止裂韧性分析预测的核心和难点。而 管道气体减压波曲线受到管径、钢级、壁厚以及输送气体组分、运行温度、压力等的强烈影 响。国际上几十年来进行了大量的工作，发展了基于相应的气体减压波分析的止裂韧性预 测方法及软件，主要使用的有GasDecom、PipeDecom, HLP等，但是由于版权问题，我国还没 有一套类似软件，其中的气体减压波分析模型和方法也鲜有公开发表。因此，研究开发适合 于高压天然气管线减压波分析预测方法及软件已成当务之急。
技术实现思路
本专利技术的目的是提供一种，在对国际上高压 输气管道止裂韧性预测方法特别是气体减压波分析预测方法进行调查分析研究的基础上， 结合我国管道工程特点，研究建立天然气气质、温度、压力、管径等与压力波曲线的关系，开 发含气液两相流的多组分天然气的减压波曲线计算分析方法及软件，用于重大输气管道工 程止裂韧性预测。本专利技术采用的技术方案是步骤A 测得管线开裂点的压力及温度的信息，或者开裂点位置及管线入口处压 力、温度、管道基础资料，如管长、管径、管道粗糙度、传热系数、流量。其次，化验取得开裂前 管线内的气质组成的资料信息，并输入计算机；步骤B:在测得的管线某点开裂前初始压力P1、温度Tl下，根据不同的气质选 择相应的状态方程，即选用索阿韦-雷德利克-邝(Saove-Redlich-Kwong，简称SRK)方 程、彭-鲁宾逊(Peng-Robinson，简称PR)方程、本尼迪克特-韦布-鲁宾-斯塔林方程 和彭-鲁宾逊方程(Benedict-Webb-Rubin-Starling 方程和 Peng-Robinson 方程，简称 BffRS-PR方程)、本尼迪克特_韦布-鲁宾-斯塔林方程和索阿韦_雷德利克_邝方程(Benedict-Webb-Rubin-Starling 方程-Saove-Redlich-Kwong 方程，简称 BWRS-SRK 方程) 四个状态方程及组合，计算等熵点压力温度下流体的密度(P )、逸度(f)、定容热容(Cv)、定 压热容(Cp)、比热(Y)流体特性，用于判断管线中流体所处相态，描述真实流体等熵过程的 特性，以便准确计算减压波速(w)。SRK、PR、BffRS状态方程分别为 p_ RT ( α(Τ) 其中，状态方程计算中采用的混合物规则为经典的范德华(VanderWaals)单流体 混合规则，也称二次型混合规则，是目前最常用的混合规则，其形式如下在该混合准则中，交互作用参数采用以下的组合规则 式中P-—压力，T-—温度，R-—气体常数，ν-—摩尔体积，P -—密度，i-—组 分名称，j—组分名称，A0-特征参数，B0-特征参数，C0-特征参数，D0-特征参数， E0—特征参数，a—特征参数，b—特征参数，c—特征参数，d—特征参数，α ---特 征参数，Y —特征参数， -一交互作用参数，…一交互作用参数，bi—交互作用参数， br-交互作用参数，交互作用参数，交互作用参数，交互作用参数， Xi—i组分含量，Xj—j组分含量。步骤C 计算减压过程中的熵熵的计算直接决定着流体所处的状态，影响到气液 相的分布及组成，进而影响到减压波的计算。S = S0+(S-S0)对于压力为显函数的状态方程，混合物的熵值为 S0为混合物在系统温度T下的理想气体熵，表达式为 式中S—混合物的熵，S0—理想气体混合物的熵，P—压力，T—温度，R—气 体常数，V—体积，Cp-定压热容，η—物质的摩尔数。根据以上方程，计算出初始状态下气质的熵，记为Sp步骤D 根据步骤B，判断压力1 (P)、温度I(T1)下的相态。若为单相，则T2 = \-ΔΤ，ΔΤ取1开尔文。假设压力2 (P2)接近于压力I(P1),由步骤C求出压力2扒)、温度2 012)下的熵 2 (S2)0计算S2-S1 < ε，若是，则进入步骤Ε;若否，则减小压力(P2),重新计算，直至满足 条件。若为两相，则通过减小压力(P)JPP2 = P1-AP, ΔΡ取0.0001 O-Oi^P1。假设温度2(Τ2)接近于温度I(T1),由步骤C求出温度2(Τ2)，压力2(Ρ2)下的熵 2 (S2)。计算S2-S1) I < ε，若是，则进入步骤F;若否，则减小温度2 (T2)，重新计算，直至满 足条件。步骤E 计算压力2 (P2)、温度2 (T2)下的气体密度(见步骤B)和声速。气相声速计算模型为 通过分析得 c = aG式中，ae-—气体声速，P-—压力，P -—密度，T-—温度，R-—气体常数，Y -—气 体比热，M-气体摩尔质量，C—管线中流体声速。步骤F 两相流声速计算在查阅大量文献，以及对已有模型和新建声速模型进行 比较的基础上，最终确定采用中国石油大学(北京)研究的两相流声速计算，其中气体采用 经典的气相声速计算模型(见步骤E)，液相声速方程采用泰特(Tait)刚性气体模型。液相声速计算模型为 其中 两相流声速计算模型为宫-徐(Gong-Xu)模型(中国石油大学(北京)建立)，该 模型在双流体模型基础上，引入虚拟质量力系数，推导了不同流型工况气液两相管流压力 波速预测的统一模型 其中，Cvm为虚拟质量力系数。(1)当Cvm = 0时，适用于分层流(忽略相间阻力) (2)当Cvm = 0.51 + 27 , I-Rn时，适用于泡状流； (3)当6； ^ =3.3 + 1.7时，适用于段塞流；15 / 45 J(4)当Cvm —c 时，适用于均质小泡状流、雾状流。式中，Rl—液相体积分数，Rg—气相体积分数，常数η = 7. 15，B = 3. 3 X IO8Pa, P0l —对比液相(通常选用水)密度，P^—液相密度，Pe—气相密度，—液体声速， aG—气体声速，C—管线中流体声速。步骤G 根据减压波过程中各个等熵点的密度差和速度差，计算开裂处流体的流 出速度。2二2 (Au)sAu = c式中，u---流体从开裂处的流出速度，Au---不同等熵点的速度差，P---压力， P ---密度，Δ ρ ---等熵点的密度差，c—管线中流体声速。步骤H:由步骤E (F)和G，声速(c)和减压波速度(U)之差求得减压波传播速度 (w)。若减压波传播速度(w)接近本文档来自技高网...

【技术保护点】
一种高压输气管线减压波分析预测方法，其特征是：　　步骤Ａ：测得管线开裂点的压力及温度的信息，或者开裂点位置及管线入口处压力、温度、管道基础资料，如管长、管径、管道粗糙度、传热系数、流量，其次，化验取得开裂前管线内的气质组成的资料信息，并输入计算机；　　步骤Ｂ：在测得的管线某点开裂前初始压力Ｐ１、温度Ｔ１下，根据不同的气质选择相应的状态方程，即选用索阿韦－雷德利克－邝（Ｓａｏｖｅ－Ｒｅｄｌｉｃｈ－Ｋｗｏｎｇ，简称ＳＲＫ）方程、彭－鲁宾逊（Ｐｅｎｇ－Ｒｏｂｉｎｓｏｎ，简称ＰＲ）方－交互作用参数，ｋ↓［ｉｊ］－－－交互作用参数，ｘ↓［ｉ］－－－ｉ组分含量，ｘ↓［ｊ］－－－ｊ组分含量；　　步骤Ｃ：计算减压过程中的熵：熵的计算直接决定着流体所处的状态，影响到气液相的分布及组成，进而影响到减压波的计算，　　Ｓ＝Ｓ↑［０］＋（Ｓ－Ｓ↑［０］）　　对于压力为显函数的状态方程，混合物的熵值为　　Ｓ＝Ｓ↑［０］＋（Ｓ－Ｓ↑［０］）＝Ｓ↑［０］＋＊［（＊ｐ／＊Ｔ）Ｖ－Ｒ／Ｖ］ｄＶ＋Ｒ＊ｎ↓［ｉ］ｌｎＶ／ｎ↓［ｉ］ＲＴ　　Ｓ↑［０］为混合物在系统温度Ｔ下的理想气体熵，表达式为　　Ｓ↑［０］＝＊Ｃ↓［ｐ］ｄＴ／Ｔ　　式中：Ｓ－－－混合物的熵，Ｓ↑［０］－－－理想气体混合物的熵，Ｐ－－－压力，Ｔ－－－温度，Ｒ－－－气体常数，Ｖ－－－体积，Ｃ↓［Ｐ］－－－定压热容，ｎ－－－物质的摩尔数，　　根据以上方程，计算出初始状态下气质的熵，记为Ｓ↓［１］；　　步骤Ｄ：根据步骤Ｂ，判断压力１（Ｐ↓［１］）、温度１（Ｔ↓［１］）下的相态，　　若为单相，则Ｔ↓［２］＝Ｔ↓［１］－ΔＴ，ΔＴ取１开尔文，　　假设压力２（Ｐ↓［２］）接近于压力１（Ｐ↓［１］），由步骤Ｃ求出压力２（Ｐ↓［２］）、温度２（Ｔ↓［２］）下的熵２（Ｓ↓［２］），计算｜Ｓ↓［２］－Ｓ↓［１］｜＜ε，若是，则进入步骤Ｅ；若否，则减小压力（Ｐ↓［２］），重新计算，直至满足条件，　　若为两相，则通过减小压力（Ｐ），即Ｐ↓［２］＝Ｐ↓［１］－ΔＰ，ΔＰ取０．０００１～０．００２Ｐ↓［１］，　　假设温度２（Ｔ↓［２］）接近于温度１（Ｔ↓［１］），由步骤Ｃ求出温度２（Ｔ↓［２］），压力２（Ｐ↓［２］）下的熵２（Ｓ↓［２］）；计算｜Ｓ↓［２］－Ｓ↓［１］）｜＜ε，若是，则进入步骤Ｆ；若否，则减小温度２（Ｔ↓［２］），重新计算，直至满足条件，　　步骤Ｅ：计算压力２（Ｐ↓［２］）、温度２（Ｔ↓［２］）下的...

【技术特征摘要】

【专利技术属性】
技术研发人员：冯耀荣，宫敬，庄传晶，霍春勇，吉玲康，杨龙，赵建奎，邱伟伟，
申请(专利权)人：中国石油天然气集团公司，中国石油天然气集团公司管材研究所，
类型：发明
国别省市：11[中国|北京]