CN104899414A - 一种电-气互联能源系统可用输电能力计算方法 - Google Patents
一种电-气互联能源系统可用输电能力计算方法 Download PDFInfo
- Publication number
- CN104899414A CN104899414A CN201510174989.2A CN201510174989A CN104899414A CN 104899414 A CN104899414 A CN 104899414A CN 201510174989 A CN201510174989 A CN 201510174989A CN 104899414 A CN104899414 A CN 104899414A
- Authority
- CN
- China
- Prior art keywords
- atc
- natural gas
- electric
- constraint
- static security
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公布了一种电-气互联能源系统可用输电能力计算方法。电力系统可用输电能力定义了不同区域间的功率交换能力,ATC信息对于电力市场参与者而言至关重要。然而,随着燃气轮机组发电比重的显著提升,传统ATC计算并没有考虑NGFPP的一次能源供应以及天然气系统的运行约束,显然是不恰当的。为此,本发明研究了计及电-气互联能源系统静态安全约束的电力系统ATC计算,并采用线性预测法求解,线性预测法首先基于灵敏度因子辨识关键约束,然后联立求解关键约束与稳态能量流方程,得ATC的计算值。
Description
技术领域
本发明涉及一种电-气互联能源系统可用输电能力计算方法,属于电-气互联能源系统分析、优化、市场技术领域。
背景技术
以不可再生的化石燃料作为一次能源的传统发电机组环境污染大,难以适应现代智能电网、低碳社会的发展需求。相比而言,天然气储量丰富,且以天然气为一次能源的燃气轮机(natural-gas fired power plants,NGFPP)有如下优点:1)发电效率高;2)碳排放量小;3)动作响应块。因而近年来NGFPP发电比例显著提高,电力系统与天然气系统的耦合不断增加。同时,随着电转气技术的成熟,有望实现电能的大规模存储,从而为风能、太阳能等可再生能源的消纳提供了保障,更深层次地加深电力系统与天然气系统之间的耦合。
传统意义上,电力系统与天然气系统是相互独立、高度自治的系统,二者的规划、运行、控制等方面均由各自的运营机构独立执行。但随着二者之间的高度融合,相互独立的运营方式未必适用。例如,当电力系统可再生能源出力、电负荷功率波动时,部分功率波动需要由天然气系统平衡。另一方面,天然气的市场价格直接决定了NGFPP的运行成本,继而影响了电力系统的优化调度;同时,当天然气系统发生大扰动或故障时,直接影响了NGFPP的天然气供应,进一步威胁到电力系统的安全性与可靠性。因而,迫切需要从统一、协调的方式的视角去分析电力系统与天然气系统。
电力系统可用输电能力(available transfer capacity,ATC)衡量了在满足电力系统安全约束的基础上互联电网不同区域间的功率交换能力。ATC的计算结果对于电力市场的买卖双方至关重要,同时也可用于电力系统的规划、安全稳定评估。ATC的求解方法主要包括直流潮流法、重复潮流法、最优潮流法、连续潮流法。值得引起注意的是,现有ATC求解方法仅仅计及了电力系统的安全约束,并没有考虑一次能源系统-天然气系统的运行约束。特别是当NGFPP发电比重较大时,NGFPP输出功率的大额度增长可能会迫使天然气系统的运行状态超出安全约束(例如节点压力越限、管道流量越限),从而影响NGFPP的一次能源供应。因而NGFPP的出力约束不仅仅取决于自身的装机容量,天然气系统的安全约束也是不可忽略的。
本发明研究了计及电-气互联能源系统静态安全约束的ATC计算。并采用线性预测法求解,首先基于灵敏度因子辨识制约ATC的关键约束,然后联立求解能量流方程与关键约束,实现了ATC的快速、准确计算,同时线性预测提供的关键约束信息恰好为天然气系统运行状态的调整提供了参考。最后,实际系统的算例测试验证了本发明所提方法的 有效性。
发明内容:
发明目的:本发明提出一种电-气互联能源系统可用输电能力计算方法,利用线性预测法考虑了天然气系统的运行约束,计算快速准确,而且为天然气系统运行状态的调整提供了参考。
技术方案:本发明提出一种基于线性预测法的电-气互联能源系统可用输电能力计算方法,包括以下步骤:
1)输入电力系统的实时运行参数,包括网络参数、拓扑、发电机实时出力、电负荷功率以及电力系统静态安全约束,包括发电机有功出力约束、节点电压幅值约束、支路功率约束;
2)输入天然气系统的实时运行参数,包括网络参数、拓扑、储气点流量供应、气负荷流量以及天然气系统静态安全约束,包括气源流量供应约束、节点压力约束、管道流量约束;
3)根据当前运行状态,分别求解电力系统、天然气系统状态量对负荷参数的灵敏度因子;
4)基于灵敏度因子,辨识制约电-气互联能源系统ATC的关键约束;
5)联立求解关键约束与稳态能量流方程,得互联能源系统的运行状态;
6)判断当前运行状态是否满足互联能源系统的静态安全约束(包括电力系统、天然气系统的静态安全约束),是则转下一步,否则在当前运行状态重新计算灵敏度因子,重复线性预测法的步骤;
7)判断是否需要调整天然气系统的运行状态,是则调整天然气系统的初始运行状态,否则转下一步;
8)输出两区域间的ATC信息。
优选地,所述电力系统静态安全约束包括:
式中:V分别为节点电压幅值;|Sij|为支路视在功率;PG为发电机输出功率;Vmax、Vmin分别为节点电压幅值约束的上下限;为支路视在功率约束上限;分别为发电机有功出力约束的上下限。
优选地,所述天然气系统静态安全约束包括:
式中:Fmn为管道流量;π分别为节点压力;Fsup为流量供应量;分别为节点压力约束的上下限;为管道输送流量的最大值;分别为气源流量供应约束的上下限。
优选地,所述灵敏度因子包括电力系统状态量对负荷参数λ的偏导,即 以及天然气系统状态量对负荷参数λ的偏导,即
优选地,所述基于线性预测法的ATC计算方法的数学模型为:
定义λe对应的约束即为电力系统静态安全约束下制约ATC的关键约束,λg对应的约束即为天然气系统静态安全约束下制约ATC的关键约束。比较λe与λg的大小,基于线性预测的ATC计算可分为两种情况:
1)λe<λg
若λe<λg,则可不考虑天然气系统的静态安全约束,直接联立求解电力潮流方程与λe对应的关键约束:
式中:fe(V,θ,λ)=0为电力系统等式约束,θ为节点电压相角;he(V,θ,λ)=0为电力系统静态安全约束下制约ATC的关键约束。
2)λe>λg
在这类情况下,ATC计算有必要计及天然气系统的静态安全约束。即统一求解互联能源系统的能量流方程与λg对应的关键约束:
式中:f(V,θ,π,λ)=0互联能源系统稳态能量流方程;hg(π,λ)=0为天然气系统静态安全约束下制约ATC的关键约束。
优选地,所述天然气系统初始运行的调整为:
当制约ATC的关键约束为天然气系统静态安全约束时(即λe>λg),可通过调整天然气系统的初始运行状态,例如增加关键约束附近的储气量、提高加压站升压比,扩大天然气系统的静态安全域,以增强电力系统不同区域间的输电能力。
有益效果:1)当研究二次能源系统-电力系统不同区域间的输电能力时,有必要计及一次能源系统-天然气系统的静态安全约束,以保证NGFPP一次能源供应的安全性与可靠性;2)线性预测法既实现了ATC的快速、准确计算,同时又为天然气系统运行状态的调整提供了参考;3)ATC信息对于电力、天然气市场参与者同等重要。以统一的能源市场分析不同区域间的ATC,在保证互联能源系统安全性的基础上,有望给整个能源市场的参与者带来更大的经济收益。
附图说明:
图1:本发明流程图;
图2:修改的比利时输气系统结构图;
图3:IEEE39节点系统结构图。
具体实施方式:
下面结合附图对发明的技术流程进行详细说明:
1 电-气混联能源系统ATC计算数学模型
同时计及电力系统、天然气系统静态安全约束的ATC数学模型为:
obj.maxλ
s.t.f(V,θ,π,λ)=0
gmin≤g(V,θ,π,λ)≤gmax
式中:f(V,θ,π,λ)=0为互联能源系统稳态能量流方程;gmin≤g(V,θ,π,λ)≤gmax为不等式约束,包括电力系统静态安全约束,天然气系统静态安全约束。
基于连续潮流法计算ATC,从当前运行状态出发,逐步增大λ,求解稳态能量流,直至静态安全约束越限。连续潮流法的优点在于易于处理各种约束,计算精度高,但由于需要求解多次潮流,计算复杂度高,一般仅用于离线分析。
2 混联能源系统状态量对ATC的灵敏度因子计算
电力系统状态量对λ的灵敏度计算:
1)发电机有功出力:
式中:KG为发电机功率增长方向。
2)节点电压幅值:
3)支路功率:
天然气系统状态量对λ的灵敏度计算:
1)节点流量注入:
2)节点压力:
3)管道流量:
3 关键约束的辨识
基于电力系统状态量与λ的灵敏度关系,分别可得各约束下λ的预测值:
式中:λV、分别为发电机出力、节点电压幅值、支路视在功率约束下λ的预测值。
定义则λe对应的约束即为电力系统静态安全约束下制约ATC的关键约束。
同理,对于天然气系统,各约束下λ的预测值可表示为:
式中:λπ、分别为气源流量供应、节点压力、管道流量约束下λ的预测值。
定义则λg对应的约束即为天然气系统静态安全约束下制约ATC的关键约束。
4 关键约束与能量流的联立求解
比较λe与λg的大小,基于线性预测的ATC计算可分为两种情况:
1)λe<λg
若λe<λg,则说明制约ATC的关键约束为电力系统静态安全约束,可不考虑天然气系统的静态安全约束,直接联立求解电力潮流方程与λe对应的关键约束:
式中:fe(V,θ,λ)=0为电力系统等式约束;he(V,θ,λ)=0为电力系统静态安全约束下制约ATC的关键约束。
2)λe>λg
在这类情况下,制约ATC的关键约束为电力系统静态安全约束,因而ATC计算有必要计及天然气系统的静态安全约束。即统一求解互联能源系统的能量流方程与λg对应的关键约束:
式中:hg(π,λ)=0为天然气系统静态安全约束下制约ATC的关键约束。
线性预测法仅需求解2次稳态(能量流)方程,相比于传统连续潮流法,计算效率大大提高。且在关键约束辨识准确的情况下,ATC的计算结果是严格精确的。
5 此外,由于互联能源系统本质上是非线性的,基于线性预测法辨识关键约束,可能会存在误辨识的情况。因而求得λ后,需重新校验当前运行状态是否越限,若存在越限,则在当前运行点重复线性预测法的步骤,ATC计算结果对应的运行状态满足所有静态安全约束。
所述电力系统静态安全约束包括:
式中:V分别为节点电压幅值;|Sij|为支路视在功率;PG为发电机输出功率;Vmax、Vmin分别为节点电压幅值约束的上下限;为支路视在功率约束上限;分别为发电机有功出力约束的上下限,i为电力系统中第i个节点编号。
所述天然气系统静态安全约束包括:
式中:Fmn为管道流量;π分别为节点压力;Fsup为流量供应量;分别为节点压力约束的上下限;为管道输送流量的最大值;分别为气源流量供应约束的上下限,m为天然气系统中第m个节点编号。
6 天然气系统初始运行状态的调整
天然气系统与电力系统的显著差别在于:天然气是易于存储的,电能则难以大规模存储。天然气的储气量对于天然气系统的安全性、可靠性至关重要。此外,当网络中节点压力偏低时,通过提高加压站的升压比,可增强天然气系统的输气能力。
当制约ATC的关键约束为天然气系统静态安全约束时(即λe>λg),可通过调整天然气系统的初始运行状态,例如增加关键约束附近的储气量、提高加压站升压比,扩大天然气系统的静态安全域,以增强电力系统不同区域间的输电能力。
算例分析
为验证所提方法的有效性,本发明电-气互联能源系统由中比利时输气系统(如图2)以及Matpower5.0软件中IEEE39节点系统(如图3)构成。将所有发电机有功出力上限提高50%,且假定节点30、31、32、34、37所连发电机为NGFPP。
1 线性预测法有效性测试
相比于传统连续潮流,线性预测法以灵敏度因子辨识关键约束,然后联立求解能量流方程与关键约束,避免了能量流方程的连续计算,大大提高了计算效率。连续潮流与线性预测法ATC计算结果比较如表1所示。
由表1可知,线性预测法准确辨识了制约ATC的关键约束,线性预测法的ATC计算结果也与连续潮流法基本相同。因而线性预测法实现了ATC的快速、准确计算,确保了能源市场的参与者能够获知系统实时的ATC信息。
表1 连续潮流与线性预测法ATC计算结果比较
2 天然气系统初始运行状态的调整
若天然气系统的静态安全约束制约了电力系统的ATC,可考虑增加储气点的天然气供应、提高加压站升压比,以增强天然气系统的输气能力,扩大天然气系统的静态安全域。
对于区域2到区域1、3的ATC,天然气系统节点6压力约束为关键约束。增加节点5天然气的供应量,则区域2到区域1、3的ATC变化如表2所示。
表2 节点5天然气供应量对ATC的影响
由表2可知,通过增加节点5天然气的供应,能有效提高区域2到区域1、3的ATC。当节点5天然气供应增加至2MMm3时,场景2与场景1下区域2到区域1、3的输电能力完全一致。
同理,对于区域3到区域1、2的ATC,将节点17与18之间的加压站升压比由1.2提高至1.3,则场景2区域3到区域1、2的ATC也提高至表1中场景1的计算结果。
通过调整天然气系统的初始运行点,可将制约ATC的关键约束由天然气静态安全约束转化为电力系统静态安全约束,从而增强了不同区域间的输电能力。而线性预测法提供的关键约束信息恰好为天然气系统运行状态的调整提供了参考。
Claims (6)
1.一种电-气互联能源系统可用输电能力计算方法,包括输入电力系统和天然气系统的实时运行参数,其特征在于,包括以下步骤:
分别求解电力系统和天然气系统的各个状态量对负荷参数λ的灵敏度因子;
根据灵敏度因子来辨识制约电力系统的关键约束λe和制约天然气系统的关键约束λg;
当λe<λg时,联立求解电力系统静态安全约束与稳态能量流方程,得到互联能源系统的运行状态,当λe>λg时,联立求解天然气系统静态安全约束与稳态能量流方程,得到互联能源系统的运行状态;
判断当前运行状态是否同时满足电力系统和天然气系统的静态安全约束,是则转下一步,否则在当前运行状态重新计算灵敏度因子,重复线性预测法的步骤;
当λe>λg时,调整天然气系统的运行状态。
2.根据权利要求1所述的电-气互联能源系统可用输电能力计算方法,其特征在于,所述电力系统静态安全约束包括:
式中:V分别为节点电压幅值;|Sij|为支路视在功率;PG为发电机输出功率;分别为节点电压幅值约束的上下限;为支路视在功率约束上限;分别为发电机有功出力约束的上下限,i为电力系统中第i个节点编号。
3.根据权利要求1所述的电-气互联能源系统可用输电能力计算方法,其特征在于,所述天然气系统静态安全约束包括:
式中:Fmn为管道流量;π分别为节点压力;Fsup为流量供应量;分别为节点压力约束的上下限;为管道输送流量的最大值;分别为气源流量供应约束的上下限,m为天然气系统中第m个节点编号。
4.根据权利要求1所述的电-气互联能源系统可用输电能力计算方法,其特征在于,所述灵敏度因子包括电力系统状态量对负荷参数λ的偏导,即以及天然气系统状态量对负荷参数λ的偏导,即
5.根据权利要求1所述的电-气互联能源系统可用输电能力计算方法,其特征在于,所述基于线性预测法的ATC计算方法的数学模型为:
定义λe对应的约束即为电力系统静态安全约束下制约ATC的关键约束,λg对应的约束即为天然气系统静态安全约束下制约ATC的关键约束,比较λe与λg的大小,基于线性预测的ATC计算可分为两种情况:
1)λe<λg
若λe<λg,则可不考虑天然气系统的静态安全约束,直接联立求解电力潮流方程与λe对应的关键约束:
式中:fe(V,θ,λ)=0为电力系统等式约束;he(V,θ,λ)=0为电力系统静态安全约束下制约ATC的关键约束;
2)λe>λg
在这类情况下,ATC计算有必要计及天然气系统的静态安全约束,即统一求解互联能源系统的能量流方程与λg对应的关键约束:
式中:f(V,θ,π,λ)=0互联能源系统稳态能量流方程;hg(π,λ)=0为天然气系统静态安全约束下制约ATC的关键约束。
6.根据权利要求1所述的电-气互联能源系统可用输电能力计算方法,其特征在于,所述天然气系统初始运行的调整包括:增加关键约束附近的储气量或提高加压站升压比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510174989.2A CN104899414B (zh) | 2015-04-14 | 2015-04-14 | 一种电‑气互联能源系统可用输电能力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510174989.2A CN104899414B (zh) | 2015-04-14 | 2015-04-14 | 一种电‑气互联能源系统可用输电能力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104899414A true CN104899414A (zh) | 2015-09-09 |
CN104899414B CN104899414B (zh) | 2017-11-03 |
Family
ID=54032076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510174989.2A Active CN104899414B (zh) | 2015-04-14 | 2015-04-14 | 一种电‑气互联能源系统可用输电能力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104899414B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105939029A (zh) * | 2016-06-21 | 2016-09-14 | 中国南方电网有限责任公司电网技术研究中心 | 综合能源系统规划方案获取方法和系统 |
CN106532709A (zh) * | 2016-10-25 | 2017-03-22 | 天津大学 | 一种含配电网重构的区域综合能源系统最优潮流计算方法 |
CN111798046A (zh) * | 2020-07-01 | 2020-10-20 | 中国电力科学研究院有限公司 | 计及atc的省间中长期交易出清方法、装置及存储介质 |
CN112510680A (zh) * | 2020-10-26 | 2021-03-16 | 华能定边新能源发电有限公司 | 综合能源系统中分布式电源与电转气设备协同规划方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101895116A (zh) * | 2010-07-16 | 2010-11-24 | 国电南瑞科技股份有限公司 | 基于分布因子的在线可用输电能力计算方法 |
-
2015
- 2015-04-14 CN CN201510174989.2A patent/CN104899414B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101895116A (zh) * | 2010-07-16 | 2010-11-24 | 国电南瑞科技股份有限公司 | 基于分布因子的在线可用输电能力计算方法 |
Non-Patent Citations (2)
Title |
---|
CHAUDRY M ET AL.: "Multi-time period combined gas and electricity network optimization", 《ELECTRIC POWER SYSTEMS RESEARCH》 * |
李中成 等: "含大规模风电场的电力系统概率可用输电能力快速计算", 《中国电机工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105939029A (zh) * | 2016-06-21 | 2016-09-14 | 中国南方电网有限责任公司电网技术研究中心 | 综合能源系统规划方案获取方法和系统 |
CN105939029B (zh) * | 2016-06-21 | 2018-06-19 | 中国南方电网有限责任公司电网技术研究中心 | 综合能源系统规划方案获取方法和系统 |
CN106532709A (zh) * | 2016-10-25 | 2017-03-22 | 天津大学 | 一种含配电网重构的区域综合能源系统最优潮流计算方法 |
CN106532709B (zh) * | 2016-10-25 | 2020-01-03 | 天津大学 | 一种含配电网重构的区域综合能源系统最优潮流计算方法 |
CN111798046A (zh) * | 2020-07-01 | 2020-10-20 | 中国电力科学研究院有限公司 | 计及atc的省间中长期交易出清方法、装置及存储介质 |
CN111798046B (zh) * | 2020-07-01 | 2024-04-30 | 中国电力科学研究院有限公司 | 计及atc的省间中长期交易出清方法、装置及存储介质 |
CN112510680A (zh) * | 2020-10-26 | 2021-03-16 | 华能定边新能源发电有限公司 | 综合能源系统中分布式电源与电转气设备协同规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104899414B (zh) | 2017-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104734155B (zh) | 一种获取电‑气互联能源系统可用输电能力的方法 | |
CN105046369B (zh) | 一种基于能源中心的电气混联系统建模和优化调度方法 | |
CN102593829B (zh) | 计及关键支路的输电断面确定方法与装置 | |
CN103944507B (zh) | 基于逆变器型式试验的光伏电站低电压穿越性能评价方法 | |
CN105322535B (zh) | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 | |
CN102222919B (zh) | 基于改进差分进化算法的电力系统无功优化方法 | |
CN103427609B (zh) | 一种mmc的谐波特性解析方法 | |
CN103973203B (zh) | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 | |
CN105449713A (zh) | 考虑分布式电源特性的有源配电网智能软开关规划方法 | |
CN102545207B (zh) | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 | |
CN103218690A (zh) | 基于碳排放流的主动配电网用户用电碳排放量的测量方法 | |
CN105005940A (zh) | 计及相关性的电-气互联系统概率最优潮流计算方法 | |
CN104899414A (zh) | 一种电-气互联能源系统可用输电能力计算方法 | |
CN102856899B (zh) | 一种微电网降低网损的方法 | |
CN105633948A (zh) | 一种分布式风电、光伏发电等不确定能源接入电力系统随机模糊潮流算法 | |
CN104009484A (zh) | 一种基于潮流计算的电网降损评估方法 | |
CN106786625A (zh) | 基于分布式供能系统互动能力的配电网电压协调控制方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补系统随机稳定性分析方法 | |
CN107887908A (zh) | 基于vsc与高压配电网拓扑重构的输电网阻塞管控方法 | |
CN105958496A (zh) | 基于连续潮流的含vsc系统最大输电能力计算方法 | |
CN104617576B (zh) | 考虑直流控制特性的多直流馈出交流电网故障计算方法 | |
CN106208099A (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
CN104217251A (zh) | 基于k2算法的装备故障贝叶斯网络预测方法 | |
CN105490282A (zh) | 考虑微电源无功出力均衡度的微电网实时电压控制方法 | |
CN105281327B (zh) | 考虑离散和连续决策变量的大规模配电网最优潮流计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |