一种电-气互联综合能源系统静态安全域分析方法 |
|||||||
申请号 | CN201510707513.0 | 申请日 | 2015-10-27 | 公开(公告)号 | CN105356447A | 公开(公告)日 | 2016-02-24 |
申请人 | 河海大学; | 发明人 | 卫志农; 陈胜; 孙国强; 孙永辉; 臧海祥; 朱瑛; 陈霜; 张思德; | ||||
摘要 | 本 发明 公布了一种电-气互联综合 能源 系统静态安全域分析方法。首先建立了电-气互联综合能源系统稳态 能量 流模型;接着参照电 力 系统静态安全域的定义与模型,提出电-气互联综合能源系统静态安全域的概念;然后通过重复能量流计算,获取了安全域边界上的一组边界运行点;最后基于线性 超平面 拟合,可得电-气互联综合能源系统的静态安全域。本发明所提安全域模型将为电-气互联综合能源系统的实时运行调度及在线安全评估提供了理论 基础 。 | ||||||
权利要求 | 1.一种电-气互联综合能源系统静态安全域分析方法,其特征在于:包括以下步骤: |
||||||
说明书全文 | 一种电-气互联综合能源系统静态安全域分析方法技术领域背景技术[0002] 相比于火电机组,燃气轮机组碳排放量少、动态响应快速、建设周期短,因此近年来燃气轮机发电比重逐步增加。相应地,电力系统与天然气系统的耦合也不断增大。另一方面,电转气技术在大规模储能方面的优势,有望实现电网与气网能量流的双向流动。电力系统与天然气系统的耦合互联有望从两方面促进新能源的大规模并网:1)燃气轮机组快速的动态响应特性可用于平抑新能源出力的波动;2)对于电网无法消纳的新能源,可通过电转气技术将这部分能量转化为天然气,存储于天然气网络中。因此,电-气互联综合能源系统是人类社会向低碳、可持续能源系统转变的必经之路。 [0003] 然而,电网与气网之间的深度耦合给二者的安全运行带来了潜在的风险。例如2011年美国大停气(电)事件,其演变机理如附图2所示,该事件的根源在于两个方面:1)极端天气;2)电网与气网的深度耦合。由于燃气轮机一般签订可中断合同,当气网发生输气阻塞或故障时,燃气轮机组将首先被切除;这种情况下,当燃气轮机发电比重较大时,大量电源的缺失将导致大规模的停电事故,使得天然气加压站失去可靠电力供应,致使天然气网络输气能力进一步降低,即事故从气网蔓延至电网后,又反过来从电网蔓延至气网,导致事故愈演愈烈。 [0004] 事实上,当电网与气网深度耦合时,燃气轮机的发电比重不仅取决于自身装机容量,更取决于天然气网络的安全约束,而传统的电力系统静态安全域则难以体现这点,即电力系统静态安全域结果过于乐观。发明内容: [0005] 发明目的:本发明所要解决的技术问题是解决现有电力系统静态安全域难以体现天然气网络的安全约束,即电力系统静态安全域结果过于乐观的问题。 [0006] 技术方案:本发明为实现上述目的,采用如下技术方案: [0007] 一种电-气互联综合能源系统静态安全域分析方法,包括以下步骤: [0008] 1)建立电力系统稳态能量流模型: [0009] 对于连接节点i、j的支路l,其支路功率可表示为: [0010] [0012] 同时,每个节点流入功率必须等于流出功率: [0013] TGPG-TLPL=AePl; [0014] TGQG-TLQL=AeQl; [0015] 式中:TG为电网节点-发电机关联矩阵,TL为电网节点-负荷关联矩阵,Ae为电网节点-支路关联矩阵;PG、QG为发电机注入有功功率、无功功率,PL、QL为负荷吸收的有功功率、无功功率; [0016] 2)建立天然气系统稳态能量流模型: [0017] 对于连接节点m与n的管道l,流过该管道的流量为 [0018] [0019] [0020] 式中:Fl为管道l流量;π为节点压力;kl为与管道l压力损失相关的常数; [0021] 基于天然气系统各节点流量的平衡,可得: [0022] TSFS-TDFD=AgFl; [0023] 式中:TS为气网节点-气源关联矩阵,TD为气网节点-负荷关联矩阵,Ag为气网节点-管道关联矩阵;FS为气源注入流量,FD为气负荷流量; [0024] 3)燃气轮机耦合 [0025] 燃气轮机消耗的燃气量与其有功输出呈如下非线性关系: [0026] [0027] 式中:Gg为燃气轮机集合; 为燃气轮机输出有功, 为燃气轮机吸收的气流量; [0028] 4)电力系统静态安全域 [0029] 电力系统静态安全域为满足电力潮流约束以及网络静态安全约束的一组功率注入集: [0030] [0031] 式中:ye=(PG,PL,QL),xe=(VL,θ);Re为静态网络约束下(VL,θ)空间中的可行域;fe为非线性潮流函数; [0032] 5)电-气互联综合能源系统静态安全域 [0033] 参照电力系统静态安全域的定义,电-气互联综合能源系统静态安全域可定义为:满足电网与气网能量流约束以及静态安全约束的一组能量流注入集合: [0034] [0035] 式中:y=(PG,PL,QL,FS,FD)为能量流注入;x=(VL,θ,πD);R为电网与气网静态安全约束下x空间的可行域;f为非线性能量流方程; [0036] 6)安全域边界的运行点 [0037] 为保持已知量与未知量数量的相等,已知量 中某个元素(假定为 )设为未知量,相应地,电力系统能量流方程转化为: [0038] [0039] 式中:flimit对应于转化为等式约束的关键约束;P'G=(PG,1,L,PG,i-1,PG,i+1,L,PG,ng)T为修正后的功率注入向量; [0040] 通过求解修正后的能量流方程,可得安全域边界上的一个运行点;通过调整P'G的大小,可得一组边界运行点; [0041] 7)超平面拟合 [0042] 假定安全域边界为超平面,则电力系统安全域边界可表示为: [0043] [0044] 式中:CBe为电力系统关键约束集;η为CBe对应的超平面系数; [0045] 同理,天然气系统安全域边界为: [0046] [0047] 式中:CBg为天然气系统关键约束集;τ为CBg对应的超平面系数; [0048] 最终,电-气互联综合能源系统静态安全域Ω可表示为: [0049] [0050] 作为优化,所述步骤3)中燃气轮机在电网中相当于电源,同时在气网中相当于气负荷;因此燃气轮机连接了电网与气网。 [0051] 作为优化,所述步骤6)中安全域的边界对应于某个不等式关键约束转化为等式约束,该等式约束可加入到能量流方程中。 [0052] 作为优化,所述步骤7)中电-气互联综合能源系统静态安全域边界的超平面系数η与τ可采用最小二乘法拟合得到。 [0053] 有益效果:本发明与现有技术相比:当电网与气网深度耦合时,燃气轮机的发电比重不仅取决于自身装机容量,更取决于天然气网络的安全约束,而传统的电力系统静态安全域则难以体现这点,即电力系统静态安全域结果过于乐观。本发明提出电-气互联综合能源系统静态安全域模型,同时计及了电网与气网的静态安全约束,从而为电-气互联综合能源系统在线安全分析,事故前的预防控制以及事故后的校正控制提供理论基础。附图说明: [0054] 图1为本发明的流程示意图; [0055] 图2为2011年美国大停气(电)事件演变机理; [0056] 图3为NGS5节点系统示意图; [0057] 图4为静态安全域边界示意图; [0058] 图5为静态安全域可视化示意图。具体实施方式: [0059] 下面结合附图1~5对发明的技术流程进行详细说明: [0060] 建立电-气互联综合能源系统稳态能量流模型 [0061] 建立电力系统稳态能量流模型: [0062] 假定电力系统含ne+1个节点,be条线路;定义节点0为平衡节点,G={1,2,L,neg}为PV节点,L={neg+1,neg+2,L,ne}为PQ节点,Be为线路集。 [0063] 对于连接节点i、j的支路l,其支路功率可表示为: [0064] [0065] 式中:Pl与Ql分别为支路有功功率与无功功率;V与θ分别为节点电压幅值与相角,θij=θi-θj;gij与bij分别为支路电导与电纳;gsh,i与bsh,i分别为对地电导与电纳。 [0066] 同时,每个节点流入功率必须等于流出功率: [0067] TGPG-TLPL=AePl; [0068] TGQG-TLQL=AeQl; [0069] 式中:TG为电网节点-发电机关联矩阵,TL为电网节点-负荷关联矩阵,Ae为电网节点-支路关联矩阵。PG、QG为发电机注入有功功率、无功功率,PL、QL为负荷吸收的有功功率、无功功率。 [0070] 建立天然气系统稳态能量流模型 [0071] 首先需要说明的是,天然气系统可由暂态模型及稳态模型表述。其中暂态模型更为具体,但由分布参数与时变状态变量描述,计算复杂度相当大;相比而言,稳态模型复杂度低,且精度在可接受范围内。因此本文对天然气系统的建模选用稳态能量流模型。 [0072] 假定天然气系统含ng个节点和bg条线路;定义W={1,2,L,ngs}为压力已知节点集,H={ngs+1,ngs+2,L,ng}为流量注入已知节点集;对于集合H,定义S={ngs+1,ngs+2,L,ngs1}为气源节点集,D={ngs1+1,ngs1+2,L,ngs}为负荷节点集,H=S∪D。 [0073] 对于连接节点m与n的管道l,流过该管道的流量为 [0074] [0075] [0076] 式中:Fl为管道l流量;π为节点压力;kl为与管道l压力损失相关的常数。 [0077] 基于天然气系统各节点流量的平衡,可得: [0078] TSFS-TDFD=AgFl; [0079] 式中:TS为气网节点-气源关联矩阵,TD为气网节点-负荷关联矩阵,Ag为气网节点-管道关联矩阵;FS为气源注入流量,FD为气负荷流量。 [0080] 计算燃气轮机耦合 [0081] 燃气轮机在电网中相当于电源,同时在气网中相当于气负荷。因此燃气轮机连接了电网与气网。燃气轮机消耗的燃气量与其有功输出呈如下非线性关系: [0082] [0083] 式中:Gg为燃气轮机集合; 为燃气轮机输出有功, 为燃气轮机吸收的气流量。 [0084] 计算电-气互联综合能源系统静态安全域 [0085] 计算电力系统静态安全域: [0086] 电力系统静态安全域定义为满足电力潮流约束以及网络静态安全约束的一组功率注入集。 [0087] 在电力系统静态安全约束下,(VL,θ)空间的可行域为: [0088] [0089] [0090] [0091] [0092] [0093] 式中:RV、Rl、 及 分别为在电压幅值约束、支路容量约束、发电机有功与无功约束下(VL,θ)空间的可行域;Re为在所有电力系统静态安全约束下(VL,θ)空间的可行域;Sl为支路复功率, [0094] 接着,电力系统静态安全域可表示为: [0095] [0096] 式中:ye=(PG,PL,QL),xe=(VL,θ);fe为非线性潮流函数。 [0097] 计算电-气互联综合能源系统静态安全域 [0098] 参照电力系统静态安全域的定义,电-气互联综合能源系统静态安全域可定义为:满足电网与气网能量流约束以及静态安全约束的一组能量流注入集合。 [0099] 在天然气系统静态安全约束下,πD空间的可行域为 [0100] [0101] [0102] Rg=RF∩Rπ; [0103] 式中:RF与Rπ分别为气源容量约束、节点压力约束下的πD空间的可行域;Rg为所有天然气系统约束下πD空间的可行域。 [0104] 参照电力系统静态安全域模型,电-气互联综合能源系统静态安全域可表示为: [0105] [0106] 式中:y=(PG,PL,QL,FS,FD)为能量流注入;x=(VL,θ,πD);R=Re×Rg,为电网与气网静态安全约束下x空间的可行域;f为非线性能量流方程。 [0107] 计算决策空间的安全域 [0108] 系统调度人员的任务是在给定系统结构和负荷下,通过调整发电机与气源的出力达到确保系统安全、高效运行的目的。在这种应用背景下,通过固定电负荷与气负荷的吸收的能量,电-气互联综合能源系统的安全域可简化为决策空间中的安全域: [0109] [0110] 计算安全域边界的运行点 [0111] 安全域的边界对应于某个不等式关键约束转化为等式约束,该等式约束可加入到能量流方程中。以电力系统为例,为保持已知量与未知量数量的相等,已知量 中某个元素(假定为 )设为未知量,相应地,电力系统能量流方程转化为: [0112] [0113] 式中:flimit对应于转化为等式约束的关键约束;P'G=(PG,1,L,PG,i-1,PG,i+1,L,PG,ng)T为修正后的功率注入向量。 [0114] 通过求解修正后的能量流方程,可得安全域边界上的一个运行点。通过调整P'G的大小,可得一组边界运行点。 [0115] 进行超平面拟合: [0116] 假定安全域边界为超平面,则电力系统安全域边界可表示为: [0117] [0118] 式中:CBe为电力系统关键约束集;η为CBe对应的超平面系数。 [0119] 同理,天然气系统安全域边界为: [0120] [0121] 式中:CBg为天然气系统关键约束集;τ为CBg对应的超平面系数。 [0122] 需要说明的是,上述天然气安全域边界计及了燃气轮机组消耗的天然气 进一步的,将燃气轮机组的有功输出 代替其天然气消耗 则天然气系统静态安全域边界转化为: [0123] [0124] 同时,电-气互联综合能源系统静态安全域边界的超平面系数η与τ可采用最小二乘法拟合得到。 [0125] 基于电力系统与天然气系统静态安全域的边界,电-气互联综合能源系统静态安全域Ω可表示为: [0126] [0127] 算例分析 [0128] 本发明测试的算例为由电力系统算例IEEE39节点系统与天然气系统算例NGS5节点系统(如附录图3所示)构成。假定IEEE39节点中节点39和节点30所连的发电机为燃气轮机,且分别与NGS5节点中的节点2和节点5相连。 [0129] 基于最小二乘法拟合,IEEE39节点系统中节点37与30所连发电机以及NGS5节点系统中节点4所连储气罐的静态安全域边界如附录图4所示。定义最小二乘拟合误差ε: [0130] ε=||(PG,FS)EF-(PG,FS)LS||/||(PG,FS)EF||×100% [0131] 式中:(PG,FS)EF为根据重复潮流计算得到的边界运行点;(PG,FS)LS为最小二乘法拟合得到的边界运行点。 [0132] 误差ε的最大值和平均值分别为3.83%和0.82%,该误差精度能够满足工程实用需求。 [0133] 基于各静态安全域边界(超平面)的交集,静态安全域的可视化为附录图5所示。 [0134] 此外,超平面系数η体现了关键约束对发电机出力的灵敏度。对于天然气节点5压力约束,发电机30的超平面系数近似为发电机37的3倍。因此,当大量的电能需要由发电机30与37供应时,从天然气系统安全性角度出发,应首先考虑由发电机37供电。 |