一种电转气的综合能源系统多点线性概率能量流方法

文档序号:10562918阅读:386来源:国知局
一种电转气的综合能源系统多点线性概率能量流方法
【专利摘要】本发明公布了一种电转气的综合能源系统多点线性概率能量流方法,先提出一种含风电、热电联供以及电转气的能源集线器模型;接着定量分析电网与气网线性化能量流方程的精度,由此得到气网线性化能量流精度明显低于电网的结论;基于此,本发明提出基于多点线性化的能量流计算方法。同时,计及了电?气?热负荷与风速的相关性,基于不确定能量流注入的概率密度分布函数,采用蒙特卡罗模拟法生成随机样本;对于每个随机样本,采用多点线性法计算能量流。算例测试结果表明,本发明所提方法计算效率高,同时线性化误差在可接受范围内。
【专利说明】
一种电转气的综合能源系统多点线性概率能量流方法
技术领域
[0001] 发明涉及一种电转气的综合能源系统多点线性概率能量流方法,属于综合能源系 统不确定性分析与安全评估技术领域。
【背景技术】
[0002] 近年来,燃气轮机组装机容量与发电占比逐年增加,电力系统与天然气系统的耦 合将逐步加强。同时,燃气轮机快速的响应特性可用于平抑大规模可再生能源的波动,因此 从促进可再生能源并网的角度而言,应优先考虑安装燃气机组。当以燃气轮机平抑可再生 能源出力波动时,燃气轮机等效的气负荷(燃气轮机在天然气系统中相当于负荷)随之相应 波动,即可再生能源出力波动此时通过燃气轮机耦合传播至气网,但不同于电网,气网有良 好的灵活性应对这一波动,主要体现在两个方面:1)电能难以经济高效存储,而天然气能够 大规模存储于储气罐中;2)天然气管道linepack的存储功能,linepack的储量与节点压力 正相关,可在低谷时段存储一定量的天然气供高峰时段使用(例如,假定高峰、低谷时段负 荷分别为100MMCFD和50MMCFD,气源供应量为75MMCFD,则低谷时段25MMCFD天然气可存储于 管道中供高峰时段使用)。
[0003] 另一方面,电转气技术的发展与成熟为电能的大规模存储开辟了新的道路,相比 于蓄电池、压缩空气储能、飞轮储能等储能技术,电转气最大的优势在于能量的大规模存储 (容量达GW级)。电转气首先以电解水分解出H 2,可将出注入天然气网络,氢气/天然气混合气 可有效提高气燃烧速度,提高发电机热效率及降低排放,但混合气中氢气占比限制较大(如 英国要求不超过〇. 1 % ),电转气容量受限制;进一步的,可将H2甲烷化,利用出与0)2的化学 反应,生成大量的CH 4(天然气主要组成部分)注入至天然气网络中。
[0004] 可见,电转气技术有助于缓解输电阻塞,支撑新能源的并网,例如可将电力系统多 余的电能(比如无法消纳的新能源出力)经电转气技术转化为天然气,存储于储气罐中(天 然气易于大规模存储),然后经输气网络,由燃气轮机进一步转化为电能,此时天然气系统 充当了电能转化、存储的媒介。
[0005] 由此可知,电一气互联综合能源系统为能源领域的变革带来了新的机遇,但与此 同时,电力系统与天然气系统高度耦合运行也面临着一系列挑战,如缺乏统一的规划、运行 而带来的潜在运行风险。事实上,电网中的随机扰动可通过耦合传播至气网,进一步影响气 网的安全性,反之亦然。因此,如何定量评估不确定能量流注入对整个综合能源系统的影响 至关重要。本发明在此需求下提出一种电转气的综合能源系统多点线性概率能量流方法。

【发明内容】

[0006] 发明目的:本发明所要解决的技术问题是针对现有技术存在的不足而提供一种电 转气的综合能源系统多点线性概率能量流方法。
[0007] 技术方案:本发明提供以下技术方案:一种电转气的综合能源系统多点线性概率 能量流方法,在计算机中按以下步骤实现:
[0008] 1)基于各随机因素的概率密度分布函数,生成Ns个模特卡罗模拟的随机样本;
[0009] 2)计算各样本的等效气负荷;
[0010] 3)根据等效气负荷的概率密度分布函数,确定多个线性化点和线性化区间;
[0011] 4)对于各个线性化点,执行综合能源系统能量流计算,保留状态量及雅克比计算 结果;
[0012] 5)设置循环次数s = l;
[0013] 6)计算当前样本的等效气负荷,据此判定所属线性化区间;
[0014] 7)根据选定线性化区间,在所属区间内的线性化点附近泰勒展开,保留一阶项,得 当前运行状态并保存;
[0015] 8)s = s+l ;
[0016] 9)判断s是否大于Ns,是则转步骤10,否则转步骤6;
[0017] 10)输出运行状态概率计算结果,结束。
[0018] 作为优化,所述步骤3)中:
[0019] 气网的线性化潮流方程模型,nmn为管道mn的流量,nmn= nm-nn为节点m与n之间 压力的平方差(假定~>11 11),八(111",/^)为线性化点,5(!!:,7^)为当前运行点, 为点B的线性预测点;根据原始的非线性天然气潮流方程,行I可通过如下方 程计算得到:
(A-1)
[0021]在点A附近线性化,可得 mn
(A-2)
[0023]对比式(A-1)与(A-2),可知节点n压力值的线性化误差为:
[0026]式中:^与;rB为节点m与n在当前运行点B的压力值,;为节点n压力的线性预测 (A-3) (A - 4) (A-5) 值,为节点n压力线性预测的相对误差;[0027]电网电压降落与线路潮流的关系可表述为:
(A-6)
[0030]式中:Pij与Qij分别为线路有功与无功功率;rij与Xij分别为线路电阻与电抗;Vi为 节点i电压幅值;为线路ij电压降落矢量值;
[0031 ]引入多点线性化法含多个线性化点和线性化区间,2N+1个线性化点 {FDm, 、断,…, 定义为: FD>K,〇t-2NaDM, Fl〇rFl〇t+2NaD_10t
[0032] : : (A-7) F二-2a_,疗dm
[0033] 式中:Pro(F《i7二,)=.0,5 , PrQ(F < F^m) = 0.5 + iV / (2N + 2)., 、=《,-〇/(2 AO;
[0034] 基于式(A-7)生成的2N+1线性化点,2N+1个线性化区间他,R2,R3,…,R 2N,R2N+1}定 义为: Rl iFD,Tot I FDjot ~ ^Djot ~ FD,1o, - FD,to, + aD,to,) R2 -> [FDjo! I FlDm ~ aDjot < FDjo! < F{Djot + aLVor} rnno.l ^ ^ 1 F^ot ~ ^ ~ Fd^ ~ + |_ 0035 J 0. \A 0/ ^2 \ i i ^ \FI>.tot I ^D,t0r - ^Dio, +
[0036]作为优化,所述步骤4)中:
[0037] 建立计及电转气的能源集线器模型,能量流集线器含类型能源,其输入与输出间 的能量转化可表示为:
[0038] Le= (Pe+Pwind) ( l_Vl)Tlee+Pg,V2llCHP,e (A~9)
[0039] Lh=Pg,(v2ncHP,h+(l-V2)ngh) (A-io)
[0040] Pg,=Pg+(Pe+Pwind)vmeg (A-ll)
[0041] 式中:Pe、Pg分别为输入电功率和输入气功率;pwind为输入风电功率; V1为输入电功 率中用于电转气的比例;^为输入天然气功率中用于热电联供的比例;为变压器效率; ncHP, e与ncHP, h分别为热电联供的电效率与热效率;%h为热锅炉效率;neg为电转气效率;
[0042] 式(A-9)_(A-11)可表示为如下矩阵形式: _]卜](1-"队哪齡 v:h ir^+^i (A-12) 卜1"呀(別删 )+ (1 - )"s J W.h + G 」L Pg J
[0044] 作为优化,所述步骤7)中,对于任意运行点,计算其等效气负荷并判定所属线性化 区间,根据下述线性方程计算当前运行状态:
[0045] 万-& =-认.)1(歹-)巧:e 馬 (A-13)
[0046] 式中:Xk,Yk与Jk分别为第k个线性点能量流的状态量、已知量与雅克比矩阵,;f为 当前运行点能量流的已知量,;^为状态量的线性预测值,^^为当前运行点的等效气负 荷。
[0047] 有益效果:本发明与现有技术相比:本发明提供了一种电转气的综合能源系统多 点线性概率能量流方法,相对于现有技术而言,计及了电-气-热负荷与风速的相关性,基于 不确定能量流注入的概率密度分布函数,采用蒙特卡罗模拟法生成随机样本;对于每个随 机样本,采用多点线性法计算能量流。算例测试结果表明,本发明所提方法计算效率高,同 时线性化误差在可接受范围内。
【附图说明】:
[0048] 图1:本发明方法流程图;
[0049] 图2:含电转气的能源集线器示意图;
[0050] 图3:气网线性化潮流方程示意图;
[0051 ]图4:多点线性示意图。
【具体实施方式】:
[0052] 下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明 本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各 种等同形式的修改均落于本申请所附权利要求所限定的范围。
[0053] 1电-气互联综合能源系统能量流模型
[0054] 1.1电力潮流模型
[0055]连接节点i与j的支路有功、无功功率可表示为:
[0056] 1+ ^,)^2 - c〇s ^ sin)
[Qy ~ -(b,i + + C0S^/ ~ SiVVj Sin
[0057] 式中:Pij与Qij分别为支路ij的有功、无功功率;Vi与Vj分别为节点i与j的电压幅 值;9ij为节点i与j的电压相角差;gij与bij分别为支路ij的电导、电纳^~与1^4为节点对 地电导与电纳。
[0058] 各节点的有功、无功功率平衡方程为:
[0059] WlA (B-2)
[0060] QG',HQV (B-3)
[0061] 式中:Pc,i、Qc,i分别为节点i的发电机有功、无功功率;分别为节点i的有 功、无功负荷。
[0062] 当电网风电渗透率较高时,需要多台发电机参与频率控制,即分布式平衡节点:
[0063] PG = PG_set+APG (B-4)
[0064] APG = -KGAf (B-5)
[0065] Qg = QGsel + aQAPG + bgAP^ (B-6)
[0066] 式中:Pc_set与Qc_set发电机有功、无功功率设定值;A f为频率偏移;Kc为发电机单位 调节功率;_与1^为发电机无功功率调节因子。
[0067] 1.2天然气潮流模型
[0068]对于连接节点m与n的管道w,其流量方程为:
(B-7) f+1 Km > nn
[0070] sm=l (B-8^ l-1 mm<Wn
[0071] 式中:~为节点m压力;Fw为管道r流量。
[0072]为降低式(B-7)的非线性程度,以节点压力的平方项代替节点压力向量:
(B-9)
[0074]式中:n=ji2。
[0075]除天然气管道外,天然气网络需要配置一定数量的加压站以补偿天然气传输过程 中的压力损失。加压站消耗的能量可表示为:
(B-10)
[0077] 式中:H。为加压站消耗的能量。F。为流过加压站的流量。B。为取决于加压站工作特 性的常数,Zc = 0.239。
[0078] 若加压站由燃气轮机驱动,则需从加压站首端节点吸收额外的天然气:
[0079] r = a + [5H r + yH: :(B-11.)
[0080] 式中:T。为燃气轮机驱动加压站从气网吸收的流量。a与y为能量转换效率常 数。
[0081 ]气网的节点流量平衡方程可表示为:
[0082] TsFs-TdFd-TcTc=AwFw+AcFc (B-12)
[0083] 式中:Ts,Td与T。对应于节点-气源,节点-负荷,节点-加压站燃气轮机关联矩阵。Aw 与A。对应于节点-管道,节点-加压站关联矩阵。
[0084] 1.3电网与气网耦合模型
[0085]电网与气网之间的耦合包含三个层面:1)燃气轮机组;2)电力驱动加压站;3)计及 电转气的能源集线器模型。
[0086] 燃气轮机组消耗的天然气可通过式(B-13)表示,加压站消耗的电能即为H。。本小 节将详细介绍能源集线器模型。
[0087] Fg = Pg/% (B-13)
[0088] 式中:Fc为燃气轮机消耗的流量,%为燃气轮机转换效率。
[0089] 计及电转气的能源集线器模型如图2所示。图2中能源集线器输入能量包括风电功 率以及从电网、气网吸收的电功率与天然气功率,输出包含电负荷与热负荷。能量转换元件 包含变压器、电转气、热电联合以及热锅炉。能源集线器的输入-输出之间的能量转化关系 可表示为:
[0090] Le= (Pe+Pwind) ( l_Vl)Hee+Pg,V2HCHP,e (B~14)
[0091] Lh=Pg,(v2ncHP>h+(l-V2)ngh) (B-15)
[0092] Pg,=Pg+(Pe+Pwind)vmeg (B_ 16 )
[0093] 式中:pe、Pg分别为输入电、气功率;pwind为输入风电功率; vi为输入电功率中用于 电转气的比例;^为输入天然气功率中用于热电联供的比例;为变压器效率;批仰^与 ncHP,h分别为热电联供的电效率与热效率;%h为热锅炉效率;为电转气效率。
[0094] 式(B-14)-(B-16)可表示为如下矩阵形式: _4l「(L-V2>lcMP.e ~\\Pe+P^ f 、
[0095] = 、 ' (B-17) Lvi知a )知)々%郷+(i-&)~」[_ pg _
[0096] 1.4综合能源系统能量流模型
[0097] 计及电网与气网的耦合元素,则电力有功平衡方程(B-2)与天然气流量平衡方程 (B-12)需转换为
[0098] P(;j ~ Pu ~ Hc.i ~ Pcj = X P:i (B-18)
[0099] TsFs~T dF d~T c t c~T ggF g~T gPg=AWFW+ACF c (B-19)
[0100] 式中:Tg与Tgg分别为节点-能源集线器,节点-燃气轮机关联矩阵。
[0101] 2多点线性蒙特卡罗模拟模型
[0102] 2.1气网潮流线性化误差分析
[0103]气网的线性化潮流方程示意图是如3所示。图3中,nmn为管道mn的流量,n mn=nm-nn为节点m与n之间压力的平方差(假定nm> nn),A(nUl)为线性化点, 为当前运行点,713为点B的线性预测点。根据原始的非线性天然气潮流方程, 可通过如下方程计算得到:
[0105] 在点A附近线性化,可得:
[0107]对比式(B-20)与(B-21),可知节点n压力值的线性化误差为:
[0110] 式中:;rw与;^为节点m与n在当前运行点B的压力值,4为节点n压力的线性预测 值,見"为节点n压力线性预测的相对误差。
[0111] 由式(A-8)可知,当出现以下任何一种情况时,气网潮流线性误差较大:
[0112] ?管道流量不确定性很大(即A £_数值很大)
[0113] ?管道的长度很长、摩擦因子很大(即kmn很小)
[0114] ?在负荷处于峰值时段,气网末端的节点压力很低(即^很小)
[0115] 2.2电网潮流线性化误差分析
[0116] 电网电压降落与线路潮流的关系可表述为:
(B-24) (B-25)
[0119] 式中:Pi j与Qij分别为线路有功与无功功率;rij与Xij分别为线路电阻与电抗;Vi为 节点i电压幅值;A%为线路ij电压降落矢量值。
[0120] 由式(A-9)与(A-10)可知,若线路ij首端电压幅值I恒定,则与线路功率Pu与 Qij线性相关;即使Vi数值可变,由于其波动范围相对较小(0.9~l.lp.u.),与Pij (Qij)呈 现近似线性相关关系。
[0121] 2.3电网潮流与气网潮流线性化误差比较分析
[0122] 对比电力潮流与天然气潮流方程,可发现存在两方面的差异:
[0123] ?天然气管道潮流与节点压力平方差呈现二次非线性函数,然而线性电压降落与 线路潮流近似线性相关;由此可知,天然气潮流方程的非线性程度度明显高于电力潮流。
[0124] ?由于天然气节点压力波动范围较广(例如从300Psia到1500Psia),气网线性化 点与当前运行点"距离"可能较远;而电力节点电压波动范围较小(例如从〇.9p.u.至 1. lp.u.),电网线性化点与当前运行点"距离"相对较近。
[0125] 鉴于以上两点,可认为天然气潮流方程线性化精度低于电力潮流方程线性化精 度。传统单点线性化方法不再适用,因此本文提出一种多点线性化的方法。
[0126] 2.4多点线性化法
[0127] 本发明所提多点线性化法含多个线性化点和线性化区间,如图4所示。图4中,2N+1 个线性化点}定义为: ^2N<JD,tot^ Kjot^mot +2^aD,tot
[0128] : : (B-26) FDm=FD,,0t ^ ^UDjot > FDj〇t =FDjo, +
[0129] 式中:Pro(/^FU= 0.5 , Pro(F < F^Jot) = 0.S + N((2N + 1), 二)/(2A〇。
[0130] 基于式(A-ll)生成的2N+1线性化点,图4中2N+1个线性化区间{RhRhRs,…,R2N, R2N+i}可定义为: {FD,,〇r I Kjot ~ aDjot - FDjot ^ Km + ^DJO,} R2 iFD,rot I K,t〇r ^ ^Djot ^ FD,rot ^ K,tot + R3 ^ {FDtot | FD]tot - aDtot < FDtot < FDltot + aDtot}
[0131 ] (B-27) R2^{FDjot\Fm〇t>Fl〇t-uDjJ R2N-a Wojot I FDjot - + UD,tot}
[0132] 对于任意运行点,计算其等效气负荷并判定所属线性化区间,根据下述线性方程 计算当前运行状态:
[0133] JE-Xk^-{JI<)\Y-Yli) F^r^Rk (B-28)
[0134] 式中:Xk,Yk与Jk分别为第k个线性点能量流的状态量,已知量与雅克比矩阵。f为 当前运行点能量流的已知量,Y7为状态量的线性预测值。^;为当前运行点的等效气负 荷。
[0135] 3算例分析
[0136] 本发明测试的算例由IEEE-39节点系统与NGS-48节点系统耦合构成的综合能源系 统。两系统之间的耦合如表1所示。假定风速服从威布尔分布,其比例参数与形状参数分别 为10.7与3.97,各风电场间风速相关性为0.7。电/气/热负荷波动的标准差为其数学期望的 7%,互相关系数为0.7。
[0137] 蒙特卡罗样本总数为10000。概率结果的数学期望、标准差作为不确定性指标。以 非线性蒙特卡罗模拟结果作为参照,表2-3给出单点线性与多点线性蒙特卡罗模拟结果比 较。由表2-3可知,多点线性化精度明显高于单点线性化精度,且和非线性蒙特卡罗模拟结 果基本相近。
[0138] 另一方面,非线性蒙特卡罗模拟耗时496s,而多点线性蒙特卡罗模拟时间仅需 19s,因而多点线性化方法大幅提高了概率能量流的计算效率。
[0139] 表1 IEEE-39节点系统与NGS48节点系统耦合元素
[0142]表2气网蒙特卡罗模拟结果
[0144]注: <与 <分别为节点压力均值与标准差的相对误差。
[0145]表3电网蒙特卡罗模拟结果
[0148]注: < 与<分别为节点电压幅值均值与标准差的相对误差。
【主权项】
1. 一种电转气的综合能源系统多点线性概率能量流方法,其特征在于:在计算机中按 以下步骤实现: 1) 基于各随机因素的概率密度分布函数,生成Ns个模特卡罗模拟的随机样本; 2) 计算各样本的等效气负荷; 3) 根据等效气负荷的概率密度分布函数,确定多个线性化点和线性化区间; 4) 对于各个线性化点,执行综合能源系统能量流计算,保留状态量及雅克比计算结果; 5) 设置循环次数s = l; 6) 计算当前样本的等效气负荷,据此判定所属线性化区间; 7) 根据选定线性化区间,在所属区间内的线性化点附近泰勒展开,保留一阶项,得当前 运行状态并保存; 8. s = s+l ; 9) 判断s是否大于队,是则转步骤10,否则转步骤6; 10) 输出运行状态概率计算结果,结束。2. 根据权利要求1所述的电转气的综合能源系统多点线性概率能量流方法,其特征在 于,所述步骤3)中: 气网的线性化潮流方程模型,IImn为管道mn的流量,IImn= IIm-IIn为节点!!!与!!之间压力 的平方差(假定nm>nn),Α(Π 1, /1)为线性化点,TD为当前运行点, 为点B的线性预测点;根据原始的非线性天然气潮流方程,I可通过如下方 程计算得到:在点A附近线性化,可得?对比式(A-I)与(A-2),可知节点η压力值的线性化误差为:式中:Tw与%为节点m与η在当前运行点B的压力值,;r;f为节点η压力的线性预测值, 为节点η压力线性预测的相对误差; 电网电压降落与线路潮流的关系可表述为:式中:Pij与Qij分别为线路有功与无功功率;rij与Xij分别为线路电阻与电抗;Vi为节点i 电压幅值为线路ij电压降落矢量值; 弓丨入多点线性化法含多个线性化点和线性化区间,2N+1个线性化点 …,巧,L,,FiL,…,FXJ定义为:基于式(A-7)生成的2N+1线性化点,2N+1个线性化区间(R1,R2,R3,…,R 2n,R2N+1}定义为:3.根据权利要求1所述的电转气的综合能源系统多点线性概率能量流方法,其特征在 于,所述步骤4)中: 建立计及电转气的能源集线器模型,能量流集线器含类型能源,其输入与输出间的能 量转化可表示为: Le 一 (Pe+Pwind ) ( l_Vl ) Hee+P gV2HCHP,e ( A-9 ) Lh = P7 g(V2%HP,h+( l_V2)%h) (A-IO) P7 g = Pg+(Pe+Pwind) VlHeg (A-Il) 式中:Pe、Pg分别为输入电功率和输入气功率;Pwind为输入风电功率;^为输入电功率中 用于电转气的比例巧2为输入天然气功率中用于热电联供的比例;rw为变压器效率;%^^与 ncHP,h分别为热电联供的电效率与热效率;ngh为热锅炉效率;为电转气效率; 式(Α-9)-(Α-11)可表示为如下矩阵形式:4.根据权利要求1所述的电转气的综合能源系统多点线性概率能量流方法,其特征在 于:所述步骤7)中,对于任意运行点,计算其等效气负荷并判定所属线性化区间,根据下述 线性方程计算当前运行状态:式中:Xk,Yk与Jk分别为第k个线性点能量流的状态量、已知量与雅克比矩阵,歹为当前运 行点能量流的已知量,P为状态量的线性预测值,^为当前运行点的等效气负荷。
【文档编号】G06Q50/06GK105930980SQ201610300615
【公开日】2016年9月7日
【申请日】2016年6月8日
【发明人】卫志农, 陈 胜, 孙国强, 孙永辉, 臧海祥, 陈霜, 张思德
【申请人】河海大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1