一种天然气管道多状态故障模型及电-气耦合连锁故障模拟方法技术

本发明专利技术所提供的一种天然气管道多状态故障模型及电

【技术实现步骤摘要】
一种天然气管道多状态故障模型及电
‑
气耦合连锁故障模拟方法


[0001]本专利技术涉及综合能源系统安全性及可靠性领域，尤其涉及系统故障在电、气两系统间的传播以及初始故障导致的一系列连锁故障。

技术介绍

[0002]现如今，新能源发展迅速，世界各国都在大力发展新能源以期转变能源发展方式、减轻环境污染以及优化传统能源的结构，未来新能源将发挥越来越多的作用。社会发展对清洁能源的需求大大增加了天然气系统与电力系统的耦合程度，这就会使天然气系统的故障影响到电力系统，这对综合能源系统的安全运行产生了巨大威胁。
[0003]Y ao R等人(2017)主要研究了电力系统的连锁故障，而忽略了天然气系统带给耦合系统的影响。目前，对于电力系统的故障分析已经较为成熟。在新能源发展背景下，其它系统对于电力系统的影响亟待研究。王伟亮等人(2016)对典型电
‑
气
‑
热综合能源系统划分了层次，建立了系统的稳态模型，分析了小干扰事件对综合能源系统安全性及稳定性的影响。然而，耦合系统的运行是一个复杂的动态过程，使用本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种天然气管道多状态故障模型及电
‑
气耦合连锁故障模拟方法，其特征在于，包括以下步骤：S1：统计某地区近几年来导致天然气管道故障的因素：天然气管道泄漏的因素有老化、极端事件、低温断裂、腐蚀以及其他原因(第三方破坏、违章操作、施工不当等)；天然气管道堵塞的原因主要有积水堵塞以及冰霜堵塞；因为管道脆性断裂发生泄漏和管道冰堵多发生在温度低的春冬季，具有明显的季节性，所以这里采用月故障率，进一步的，计算故障因素权重：对于管道泄漏，有前述5种因素，在泄漏中所占权重为ω
j
(1≤j≤5)；N＝[N
1 N
2 N
3 N
4 N5]N
j
为第j种因素引起的管道泄漏事件数量，同理可用K
j
表示第j种因素引起的管道堵塞时间数量,不同因素在管道堵塞的权重为σ
j
(1≤j≤2)；进一步的，归一化故障因素，得到不同故障因素对管道故障的权重：进一步的，归一化故障因素，得到不同故障因素对管道故障的权重：S2：计算故障因素修正系数：天然气管道的故障率随使用年限的变化符合浴盆曲线，利用改进的Weibull分布来拟合：式中：T为管道使用时间，t1为管道使用初期磨合中故障偶发时段，t1～t2表示管道使用中期，管道故障率趋于稳定，t2为管道使用后期故障多发，β可以用极大似然估计法求解，t1、t2由统计数据给出；极端事件修正模型：因为极端事件突发性强，只能根据某地连续n年的给出的极端事件数据，比较本月所处情况与历史情况来进行拟合，修正值如式所示：式中，N表示待评估月发生极端事件的次数，N
i
表示统计年中第m个月发生极端事件的次数，同样的，低温断裂的修正模型以同样的方法得到；由于管道腐蚀以及其他原因导致的管道故障难以统计且难寻其规律，故只考虑其权重；管道水堵修正模型：统计连续n年管道水堵次数，拟合水堵故障率修正值：s为待评估月管道水堵次数，S
m
为当年第m月发生管道水堵的次数；通过连续几年的统计数据发现，管道冰堵通常发生在寒冷的冬季，且管道的故障率基
其中f
i
为状态i的失效频率，将S5中的式子代入上式，可以得到式子中的值；S7：采用双层蒙特卡罗方法对天然气管道故障过程进行模拟：在第一层蒙特卡罗方法中，确定现在的状态是管道泄漏还是管道堵塞；如果选中状态为管道泄漏，则继续采样两个独立的随机数来确定泄漏孔的位置和大小，设管道长度为L，泄漏点与管道起点的距离为x，统一转化为圆孔计算，孔径为d，管道的直径为D，假设泄漏前管道初始点的压力为P1，末端压力为P4，泄漏点压力为P3，对应位置管道中线上的压力为P2，每个点的温度为T
i
，泄漏量为q，原来的管道流量为Q0，全部转化为圆孔来进行计算；对0到1之间均匀分布的随机数U、U
’
进行采样x＝LUd＝DU
’
如果d/D≤0.2，此时为小孔模型，这种模型下泄漏对管道内的压力与流量的影响极小，故认为泄漏前后参数不变；如果0.2＜d/D＜0.5，此时为大孔模型，采用下式计算泄漏点后的管道流量量管道上游流量Q
so
＝Q+0.5q，管道下游流量Q
or
＝Q
‑
0.5q；A为泄漏孔面积，C
D
为流量系数，取0.9～0.98，k为绝热指数，一般取1.3，P
a
为外界环境压力，是已知量，M是天然气的摩尔质量(kg/mol)，G为流过单位面积管道的质量流量(kg/(m2.s))，通用气体常数是R(J/(mol.K))，水力摩擦系数ξ；如果d/D≥0.5，认为管道断裂，管道内流量为0；进一步的，在第二层蒙特卡罗算法中采样当前状态和过渡到下一个状态的持续时间，根据从当前状态到下一个可能状态的转换速率，对所有可能的当前状态持续时间进行采样，选择持续时间最小的状态确定下一个过渡状态，并且将此时间作为当前状态的持续时间：进一步的，如果在第一层采样中选中的状态为管道堵塞，继续对管道堵塞后实际内径D1随机采样：D1＝Dln U采用下式计算堵塞后管道内流量：
A为常数，ρ
a
为正常工况下天然气密度，T
a
为环境温度，Qo为堵塞前管道流量，Q
D
为堵塞后管道内流量，第二层采样过程与管道堵塞情况相同；S8：通过蒙特卡罗模拟天然气管道故障，计算故障后的能流，若是能流越限，天然气系统将会调整气源出力以及切负荷，切负荷目标函数需使系统运行成本及切负荷成本最小：M
G
作为天然气系统总的节点数量；C
i,G
与C
i,G,cut
分别为天然气系统在节点i处的供气成本以及削负荷成本；q
i,G
与Δq
i,G,cut
分别为天然气系统在节点i处的供气量以及天然气削减量；天然气系统故障对电力系统的影响主要是对耦合节点的供气量减少，进行最优能流计算确定Δq
i,G,cut
，正常运行时天然气系统供给燃气发电机节点的供气量为q
i,G,L
，故障后该节点的发电功率为P
i,G
＝(q
i,G,L
‑
Δq
i,G,cut
)
·
η
g
P
i,G
为燃气发电机组的有功功率，η
g
为天然气转换效率；S9：天然气系统故障将是燃气机组的供气量减少，由此其发电功率相应减少，电力系统发用不平衡，电力系统需调节发电功率及电负荷来实现功率平衡，如果调整发用已经不能平衡功率，则根据电力系统备用容量RI来确定发电量和切负荷量：式中：为电力系统节点m处传统发电机组f的最大功率，表示初始负荷，N
s
表示电力系统节点个数；若RI＞0，说明孤岛I内的备用容量能够满足负荷需求，这时只需调整发电功率与负荷达到发用平衡，若RI＜0，说明孤岛I内的备用容量无法满足负荷需求，此时就需要进行切负荷处理，负荷削减量通过故障后的供电缺额求得：K
Ic
＝min{ΔP
I
/∑L
dr,1
}K
Ic
为某时刻孤岛I的规定削减比例；ΔP
I
为孤岛I的供电缺额；∑L
dr,1
为可削减负荷点的负荷总量，若故障后的供电缺额比可削减负荷之和大，所有的可中断负荷就要以比例1完全削减，其他重要负荷的供电才能得到保证；S10：用系统负荷损失率α
g
评估天然气系统故障对两系统的影响：P
i
为正常运行时系统在节点i处的供电负荷/供气量，P
ks，i
为孤岛s第k次故障后系统在
节点i处的供电负荷/供气量，n...

【专利技术属性】
技术研发人员：黄镜月，韩璐，汪永祥，陈君祥，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：