一种气-电互联综合能源系统机组组合方法与流程

本发明属于电力信息技术领域，特别涉及了一种气-电互联综合能源系统机组组合方法。

背景技术：

随着化石能源的短缺和环境污染，天然气正迅速成为电网中新发电机组的最佳选择。在气-电互联综合能源系统中，能源综合利用系统机组是电力系统和天然气系统的传统耦合元件，能够在降低系统规划的投资成本和运行成本的同时提高系统的可靠性。

机组组合是指在调度周期内，在满足机组运行约束条件下，能够合理安排机组启停状态，调节各机组出力，从而达到经济性最优的目的。机组组合的优化在电力系统调度中占有重要地位，有利于电力系统经济运行。

然而，在气-电互联综合能源系统中存在多种不确定性，不仅有风电、电负荷的不确定性，而且还有气负荷的不确定性，这些不确定性因素会影响综合能源系统运行的稳定性、可靠性。随机规划是一种针对含有不确定性因素的优化问题的有效方法，通过多个不同场景来体现不确定性因素，用少量代表性的场景来刻画不确定性因素的波动特性，以此全面考虑不确定性因素对系统机组组合的影响。

技术实现要素：

为了解决上述背景技术提出的技术问题，本发明旨在提供一种气-电互联综合能源系统机组组合方法，能够通过随机规划的方法处理综合能源系统中气负荷、电负荷以及风电的不确定性，从而确定机组组合，提高系统运行经济性。

为了实现上述技术目的，本发明的技术方案为：

基于多场景随机规划模型的气-电互联综合能源系统机组组合方法，该方法包括构建不确定性因素的多场景和确定气-电互联综合能源系统的机组组合两个部分；所述构建不确定性因素的多场景是通过等效负荷综合考虑电负荷和风电的不确定性，同时考虑气负荷的不确定性，基于拉丁超立方抽样法生成大量场景，再通过场景削减得到所需的符合随机因素特征的场景个数；所述确定气-电互联综合能源系统的机组组合是以经济性最优为优化目标，确定电力网络约束条件、天然气网络约束条件以及二者耦合关系的约束，建立随机规划模型，求解得到气-电互联综合能源系统中最优机组组合。

进一步地，所述等效负荷定义为电负荷减去风电。

进一步地，构建不确定性因素的多场景的具体过程如下：

设风电、电负荷以及气负荷符合正态分布，采用拉丁超立方抽样，抽样规模为n，不确定性因素中随机变量个数为z，则第i个样本xi＝[xi1,.xi2,...,.xiz]，xiz为xi中的第z个随机变量，i＝1,2,...,n；

针对等效负荷、气负荷，分别采用拉丁超立方抽样产生样本空间：

(a)设随机变量xm状态为xmp，计算其数目为n×ρ(xmp)；其中ρ(xmp)为xm取xmp的概率，m,p＝1,2,...,z；

(b)生成n×z维样本矩阵xi，第m列的数据由n×ρ(xmp)个xmp构成，且每一列元素的概率和为1，即

(c)在第i次抽样中，依次从xi的每一列随机抽样得到样本值xim，i＝1,2,...,n；

对所得大量样本进行削减，过程如下：

(a)对于任意两个样本xi和xj，i≠j，其距离dij通过下式计算：

(b)去除概率密度pdi最小的样本xi：

上式中，ρi是样本xi出现的概率，ci是样本xi的密度距离，dil和dik是dij中最小的两个值，对应的样本xl和样本xk是距离样本xi最近的两个样本；

(c)通过下式更新样本xl和样本xk的概率ρl和ρk；

(d)重复步骤(a)-(c)，直至达到要求的样本数目；

通过以上步骤得到等效负荷以及气负荷个数分别为n^e和n^g，则构建的场景数个数ns＝n^en^g。

进一步地，所述经济性最优的优化目标函数如下：

上式中，i，j，t，s分别是火电机组、天然气机组、时间、场景数目的序号；np，ng，nt，ns分别是火电机组、天然气机组、时间、场景总数；ftotal是发电总成本；ρs是场景s下的概率；分别是第i台火电机组、第j台天然气机组在场景s下t时刻的出力；分别是火电机组、天然气机组的运行成本；分别是第i台火电机组t时刻的启停成本；分别是第j台天然气机组的t时刻启停成本。

进一步地，所述电力网络约束条件包括：

(1)功率平衡约束：

上式中，d是等效负荷的序号；nd是负荷的总数；是场景s下t时刻第d个等效负荷；是场景s下t时刻第d个负荷节点的切风；是场景s下t时刻第d个负荷节点的失负荷；

(2)机组出力约束：

上式中，分别是第i台火电机组与第j台天然气机组在t时刻的运行状态；分别是第i台火电机组的最大和最小出力；分别是第j台天然气机组的最大和最小出力；

(3)机组启停约束：

上式中，分别表示第i台火电机组在t时刻的开关状态；分别表示第j台天然气机组在t时刻的开关状态；

(4)机组爬坡约束：

上式中，分别是第i台火电机组的爬坡速率上下限；分别是第j台天然气机组的爬坡速率上下限；

(5)潮流约束：

上式中，sf是转移系数矩阵；kp，kg，kd分别是火电机组、天然气机组、等效负荷的关联矩阵；分别是火电机组和天然气机组在场景s下t时刻的出力矩阵；和分别是场景s下t时刻的等效负荷矩阵、切风矩阵以及失负荷矩阵；plmax是线路潮流上限矩阵。

进一步地，所述天然气网络约束条件包括：

(1)气源流量约束：

上式中，m为气源的个数；是场景s下t时刻第m个气源中天然气流量；分别是第m个气源中天然气流量的最大值和最小值；

(2)储气罐约束：

上式中，n为储气罐的个数；是场景s下t时刻第n个储气罐的容量；分别是第n个储气罐的容量的最大值和最小值；分别是场景s下t时刻第n个储气罐中流入与流出的天然气量；分别是第n个储气罐中流入的天然气量的最大值和最小值；分别是第n个储气罐中流出的天然气量的最大值和最小值；

(3)流量守恒约束：

上式中，a，b，c，d，e分别表示气源、储气罐、天然气负荷、天然气机组供电的负荷以及管道的关联矩阵；是场景s下t时刻气源中天然气流量矩阵；分别是场景s下t时刻储气罐中流入与流出的天然气量矩阵；是场景s下t时刻天然气负荷矩阵；是场景s下t时刻由天然气机组供电的负荷矩阵，由调度决定；是场景s下t时刻管道中天然气流量矩阵。

进一步地，所述电力网络和天然气网络的耦合约束：

上式中，场景s下t时刻由第j台天然气机组供电的负荷；分别是第j台天然气机组的能耗系数。

采用上述技术方案带来的有益效果：

(1)本发明用等效负荷代替电负荷和风电负荷，减少了不确定性因素个数，从而减少场景数目，提高了方法的效率；

(2)本发明不仅考虑了气-电互联综合能源系统中风电和电负荷的不确定性，而且考虑了气负荷的不确定性，有利于综合能源系统的稳定运行；

(3)本发明通过多场景随机规划，能够显著优化气-电互联综合能源系统的机组组合。

附图说明

图1是本发明的方法流程图；

图2是实施例中等效负荷的场景图；

图3是实施例中气负荷的场景图。

具体实施方式

以下将结合附图，对本发明的技术方案进行详细说明。

如图1所示，本发明设计的气-电互联综合能源系统机组组合方法，主要包括了构建不确定性因素的多场景，确定气-电互联综合能源系统的机组组合两个部分。

(1)构建不确定性因素的多场景

假设风电、电负荷以及气负荷符合正态分布，采用拉丁超立方抽样，抽样规模为n，不确定性因素中随机变量个数为z，样本表示为x＝[x1,.x2,...,xz]，则第i个样本xi＝[xi1,.xi2,...,.xiz]，i＝1,2,...,n。

为了减少场景数目，提高效率，采用等效负荷(即电力负荷减去风电负荷)，则针对等效负荷、气负荷，分别采用拉丁超立方抽样产生样本空间过程如下：

1)假设随机变量xm状态为xmp，计算其数目为n×ρ(xmp)；其中ρ(xmp)为xm取xmp的概率，m,p＝1,2,...,z；

2)生成n×z维样本矩阵xi，第m列的数据由n×ρ(xmp)个xmp构成，且每一列元素的概率和为1，即

3)在第i次抽样中，依次从xi的每一列随机抽样得到样本值xim，i＝1,2,...,n。

接下来对所得大量样本进行削减，过程如下：

1)对于任意两个样本xi和xj(i≠j)，其距离dij通过式(1)计算：

2)去除概率密度pdi最小的样本xi，如式(2)所示：

式中：ρi是样本xi出现的概率，ci是样本xi的密度距离，dil和dik是dij中最小的两个值，样本xl和样本xk是离样本xi最近的两个样本。

3)通过式(3)更新样本xl和样本xk的概率；

4)重复步骤1)～3)，直到达到要求的样本数目。

通过以上步骤可以得等效负荷以及气负荷个数分别为n^e，n^g，则构建的场景数个数ns如式(4)所示：

ns＝n^en^g(4)

注意每一个场景的概率之和为1。

假设所需等效负荷和气负荷个数都为3，通过实施上述步骤，可得到如图2、图3所示的等效负荷和气负荷的场景图。

(2)确定气-电互联综合能源系统的机组组合

在气-电互联综合能源系统中考虑混合机组，以经济性最优为目标函数：

式中：i，j，t，s分别是火电机组、天然气机组、时间、场景数目的序号；np，ng，nt，ns分别火电机组、天然气机组、时间、场景总数；ftotal是发电总成本；ρs是场景s下的概率；分别是第i台火电机组、第j台天然气机组在场景s下t时刻的出力；分别是火电机组、天然气机组的运行成本；分别是第i台火电机组t时刻的启停成本；分别是第j台天然气机组的t时刻启停成本。

电力网络约束条件如下：

1)功率平衡约束

式中：d是等效负荷的序号；nd是负荷的总数；是场景s下t时刻第d个等效负荷；是场景s下t时刻第d个负荷节点的切风；是场景s下t时刻第d个负荷节点的失负荷。考虑切风与失负荷是为了在系统负荷失去平衡时，维持系统平衡。

2)机组出力约束

式中：分别是第i台火电机组与第j台天然气机组在t时刻的运行状态；分别是第i台火电机组的最大和最小出力；分别是第j台天然气机组的最大和最小出力。

3)机组启停约束

式中：分别表示第i台火电机组在t时刻的开关状态；分别表示第j台天然气机组在t时刻的开关状态。

4)机组爬坡约束

式中：分别是第i台火电机组的爬坡速率上下限；分别是第j台天然气机组的爬坡速率上下限。

5)潮流约束

式中：sf是转移系数矩阵；kp，kg，kd分别是火电机组、天然气机组、等效负荷的关联矩阵；分别是火电机组和天然气机组在场景s下t时刻的出力矩阵；和分别是场景s下t时刻的等效负荷矩阵、切风矩阵以及失负荷矩阵；plmax是线路潮流上限矩阵。

天然气网络约束条件如下：

1)气源流量约束

式中：m为气源的个数；是场景s下t时刻第m个气源中天然气流量；分别是第m个气源中天然气流量的最大值和最小值。

2)储气罐约束

式中：n为储气罐的个数；是场景s下t时刻第n个储气罐的容量；分别是第n个储气罐的容量的最大值和最小值；分别是场景s下t时刻第n个储气罐中流入与流出的天然气量；分别是第n个储气罐中流入的天然气量的最大值和最小值；分别是第n个储气罐中流出的天然气量的最大值和最小值。

3)流量守恒约束

式中：a，b，c，d，e分别表示气源、储气罐、天然气负荷、天然气机组供电的负荷以及管道的关联矩阵；是场景s下t时刻气源中天然气流量矩阵；分别是场景s下t时刻储气罐中流入与流出的天然气量矩阵；是场景s下t时刻天然气负荷矩阵；是场景s下t时刻由天然气机组供电的负荷矩阵，由调度决定；是场景s下t时刻管道中天然气流量矩阵。

在气-电互联综合能源系统中，天然气机组作为耦合元件，消耗天然气来发电，其耦合关系约束如下：

式中：场景s下t时刻由第j台天然气机组供电的负荷；分别是第j台天然气机组的能耗系数。

最后通过优化工具即可求解得到最优机组组合。

假设系统中有一台火电机组g1，两台天然气机组g2和g3，通过求解式(5)～式(18)，即可得到如表1的最佳机组组合结果。

表1

实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。