一种天然气管道气压波动预测方法及系统技术方案

本发明专利技术涉及气

【技术实现步骤摘要】
一种天然气管道气压波动预测方法及系统


[0001]本专利技术属于气
‑
电耦合系统态势感知
，特别涉及一种天然气管道气压波动预测方法及系统。

技术介绍

[0002]随着新能源开发规模的扩大，如何提升新能源消纳水平，缓解由于其间歇性对电力系统安全稳定运行带来的影响，已成为新型电力系统构建面临的挑战。天然气发电具有负荷调节范围宽、响应快速、变负荷能力强的特点，与新能源形成互补，是电网调峰的理想选择。然而近年来多起事故说明了天然气系统与电力系统之间的耦合会带来一定的安全风险：极端恶劣天气导致气井和管道结冰，天然气供应减少，使得签订可中断供气合同的燃气轮机组停机，而其他机组也会因为快速爬坡会导致进气口管道压力低于阈值而停机。因此有必要对气
‑
电耦合系统进行实时状态感知，有助于对存在的安全风险进行预警。
[0003]除了恶劣天气之外，燃气轮机快速启停、阀门操作等现象也会对气
‑
电耦合系统带来安全风险。管道内天然气的压力、流量和管存会产生剧烈波动，当与燃气轮机连接的天然气管道气压在短时间尺度内波动程度较大时，燃气轮机将存在跳闸的可能，造成电力系统的不稳定。此外，当需要增大燃气轮机出力而快速操作阀门时会带来气锤效应，燃气轮机的输出功率会存在先减小后增大的现象，造成电力系统的功率不平衡，进而存在触发低频减载的可能。而气锤效应严重时则会产生管道爆裂等现象，进而给系统带来巨大的安全风险。
[0004]目前关于天然气管道气压波动特性的主要研究有如下方向：
[0005](1)将Weymouth稳态方程进行扩展近似刻画天然气管网的动态过程，一定程度上考虑了气网的慢动态特性，但精度难以支撑实时态势感知；
[0006](2)通过隐式差分法和特征线法实现天然气管网的动态仿真以及优化，通过逐点求解/差分高维非线性方程，计算量较大；
[0007](3)通过气路理论和状态空间法将气网动态过程转换成频率域上的二端口模型并优化求解，利用气路理论将天然气管道状态偏微分方程组转换成时域上的常微分方程组，再通过拉普拉斯变换转换为频域上的代数方程组，提高了模型精度和求解效率。已有部分研究基于这种思路提出了如下天然气管道气压可解析表达式：
[0008][0009][0010]但上述表达式无法体现流量调节时间，且并不能正确体现管道气压的波动特性。

技术实现思路

[0011]针对
技术介绍
存在的问题，本专利技术提供一种天然气管道气压波动预测方法，简单、
准确地判断出调度指令会使得天然气管道气压在何时产生多大程度的波动。
[0012]为解决上述技术问题，本专利技术采用如下技术方案：一种天然气管道气压波动预测方法，包括：
[0013]获取初始状态天然气管道参数；
[0014]设定用于生成天然气管道气压序列的时间参数；
[0015]计算流量调节开始后天然气管道首端气压在20秒内的序列；
[0016]按顺序分析天然气管道首端气压在20秒内序列，通过联机算法提取所有波动幅度不小于气压波动幅度阈值Th的子序列，得到波动预测结果。
[0017]在上述天然气管道气压波动预测方法中，所述获取初始状态天然气管道参数包括：通过SCADA系统获取初始状态下天然气管道首端气压p
s
、天然气管道末端气压p
e
、管道气流m1；获取管道的静态信息：管道长度L、管道半径R、天然气波速a；电网对燃气燃机发出的调度指令为Δm,T，从m1调节至m1+Δm。
[0018]在上述天然气管道气压波动预测方法中，所述设定用于生成天然气管道气压序列的时间参数包括：设定计算时间间隔Δt＝0.01s、时间间隔数量N，N为正整数，且NΔt≤20s，气压波动幅度阈值Th设置为0.05kPa。
[0019]在上述天然气管道气压波动预测方法中，所述计算流量调节开始后天然气管道首端气压在20秒内的序列包括：通过获取的始状态天然气管道参数和设定的用于生成天然气管道气压序列的时间参数，计算流量调节开始后，天然气管道首端气压在20秒内的序列：p0、p1、p2、
…
、p
N
；p
i
,，i＝1，2，
…
，N的可解析表达式如下：
[0020][0021]其中H为heaviside阶梯函数，erfc为误差互补函数，定义方式为：
[0022][0023]其中，c1、υ、c
l,k
、b
l,k
、e
l,k
、f
l,k
、g
l,k
，l的范围为1至10，k的范围为1至2，均为方便表达而引入的过程参数，无实际物理意义，ξ为积分操作中的被积分符号，无实际意义，过程参数定义如下：
[0024][0025][0026][0027][0028][0029][0030][0031]在上述天然气管道气压波动预测方法中，所述通过联机算法提取子序列，得到波动预测结果包括：按照顺序分析所得序列，通过联机算法以线性时间提取出所有波动幅度不小于Th的子序列p
i
,
…
,p
j
，定义如下：
[0032]|p
i
‑1‑
p
i
|<Th|p
j
‑
p
j+1
|<Th
[0033]|p
m
‑
p
m+1
|≥Th,m＝i,i+1,...,j
‑1[0034]子序列p
i
，
…
，p
j
的波动幅度为：max(p
m
)
‑
min(p
m
),m＝i,i+1,...,j
[0035]提取的所有子序列的起始时间、终止时间和波动幅度即为波动预测结果。
[0036]一种用于天然气管道气压波动预测方法的系统，包括参数获取模块，用于获取天然气管道初始参数；管道气压序列生成模块，用于生成管道气压序列；计算模块，用于计算流量调节后天然气管道气压序列；序列提取模块，用于提取子序列，得到波动预测结果。
[0037]一种计算机可读存储介质，所述计算机可读存介质用于存储程序代码，程序代码用于执行所述的天然气管道气压波动预测方法。
[0038]一种包括指令的计算机程序产品，当其在计算机上运行时，使得所述计算机执行所述的天然气管道气压波动预测方法。
[0039]与现有技术相比，本专利技术的有益效果：本专利技术基于一种可解析表达式计算出管道气压的波动幅度以波动出现的时刻为气
‑
电耦合系统的实时状态感知、安全预警提供支撑。
[0040]与差分法、特征线法等动态仿真方法比，本专利技术的计算效率大幅提高。
[0041]与其他管道气压可解析表达式的研究相比，本专利技术真实地模本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种天然气管道气压波动预测方法，其特征在于，包括：获取初始状态天然气管道参数；设定用于生成天然气管道气压序列的时间参数；计算流量调节开始后天然气管道首端气压在20秒内的序列；按顺序分析天然气管道首端气压在20秒内序列，通过联机算法提取所有波动幅度不小于气压波动幅度阈值Th的子序列，得到波动预测结果。2.根据权利要求1所述天然气管道气压波动预测方法，其特征在于，所述获取初始状态天然气管道参数包括：通过SCADA系统获取初始状态下天然气管道首端气压p
s
、天然气管道末端气压p
e
、管道气流m1；获取管道的静态信息：管道长度L、管道半径R、天然气波速a；电网对燃气燃机发出的调度指令为Δm,T，从m1调节至m1+Δm。3.根据权利要求1所述天然气管道气压波动预测方法，其特征在于，所述设定用于生成天然气管道气压序列的时间参数包括：设定时间间隔Δt＝0.01s、时间间隔数量N，N为正整数，且NΔt≤20s，气压波动幅度阈值Th设置为0.05kPa。4.根据权利要求3所述天然气管道气压波动预测方法，其特征在于，所述计算流量调节开始后天然气管道首端气压在20秒内的序列包括：通过获取的初始状态天然气管道参数和设定的用于生成天然气管道气压序列的时间参数，计算流量调节开始后，天然气管道首端气压在20秒内的序列：p0、p1、p2、
…
、p
N
；p
i
,，i＝1，2，
…
，N的可解析表达式如下：其中H为heaviside阶梯函数，erfc为误差互补函数，定义方式为：其中，c1、υ、c
l,k
、b
l,k
、e
l,k
、f
l,k
、g
l,k
，l的范围为1至10，k的范围为1至2，均为方便表达而引入的过程参数，无实际物理意义，ξ为积分操作中的被积分符号，过程参数定义如下：引入的过程参数，无实际物理意义，ξ为积分操作中的被积分符号，过程...

【专利技术属性】
技术研发人员：吴煜晖，徐箭，柯德平，廖思阳，姜新雄，陈睿，王俊，冯帅帅，
申请(专利权)人：武汉大学，
类型：发明
国别省市：