CN107591806B - 一种考虑主动配电网有功调节能力的主网日前调度方法 - Google Patents
一种考虑主动配电网有功调节能力的主网日前调度方法 Download PDFInfo
- Publication number
- CN107591806B CN107591806B CN201710927881.5A CN201710927881A CN107591806B CN 107591806 B CN107591806 B CN 107591806B CN 201710927881 A CN201710927881 A CN 201710927881A CN 107591806 B CN107591806 B CN 107591806B
- Authority
- CN
- China
- Prior art keywords
- network
- power
- power distribution
- distribution network
- major network
- 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.)
- Active
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑主动配电网有功调节能力的主网日前调度方法，属于电气工程技术领域，考虑到主网调度与配电网调度的信息收集限制，该方法通过分别收集主网与配电网常规技术参数，获取两级电网的负荷预测曲线，对主网与配电网分别进行详细的建模，同时基于主网与配网的信息交互提出了一种实现该模型高效求解的方法。本发明充分考虑了含储能的主动配电网的有功调节能力对日前调度产生的影响，经优化后的调度方案能够实现整个电网的全局最优，最大限度的降低电网发电机组的发电成本，减少资源的浪费。
Description
技术领域
本发明属于电气工程技术领域，更具体地，涉及一种考虑主动配电网有功调节能力的主网日前调度方法。
背景技术
随着能源危机与环境污染问题的日益严重，光伏风电等新能源在电网中的渗透率逐渐增加，尤其是在配电网中，新能源的渗透率不断增加，与此同时，随着电力电子技术的不断发展，储能在配电网中的分布也越来越多，原有的配电网内的可控设备逐渐增多，逐渐变为具有众多可调可控资源的主动配电网，将拥有更多的自主权，能够对储能等可控设备进行调度。另一方面在现有的日前调度体系中，一般是主网与配电网分开进行优化，首先按照一定的指标(如网络损耗，电压质量等)对配网进行优化，之后由配网向主网传输一条负荷曲线，而从主网的角度来看，配网此时即等效为一条负荷曲线来进行机组启停的优化，而当传统配电网逐渐向主动配电网转化时，就意味着配电网的调度反过来也会对主网产生一定的影响，即该负荷曲线变为一个可控的曲线，如果调节适当，能够达到削峰填谷的效果，进而降低发电成本，而现有的方法仅仅从主网侧或配网侧研究，往往不能得到全局最优的决策方案。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种考虑主动配电网有功调节能力的主网日前调度方法，由此解决现有调度方法不能取得全局最优解的技术问题。
为实现上述目的，本发明提供了一种考虑主动配电网有功调节能力的主网日前调度方法，包括：
(1)在主网与各配电网之间没有信息交换时，由主网与各配电网的目标技术参数得到主网的日前调度模型；
(2)对主网的日前调度模型进行求解得到由主网传输至各配电网的有功功率；
(3)对于任意配电网m，由引入的正向松弛变量与负向松弛变量构造可行性校验模型，并对可行性校验模型进行求解得到配电网m的目标函数值以及配电网m中储能的充放电状态值；
(4)若目标函数值为0，则表示由主网向配电网m传输的有功功率满足配电网m的功率需求，并执行步骤(5)，若目标函数值不为0，则表示由主网向配电网m传输的有功功率不满足配电网m的功率需求，并执行步骤(6)；
(5)主网向配电网m输出调度方案；
(6)将配电网m中储能的充放电状态值替换可行性校验模型中原有的充放电状态值，并根据替换后的可行性校验模型得到对偶乘子，由对偶乘子以及目标函数值得到配电网m向主网反馈的约束，利用配电网m向主网反馈的约束对主网的日前调度模型进行重新求解得到主网传输至配电网m 的新的有功功率，并执行步骤(3)。
优选地，在步骤(1)中，所述主网的日前调度模型为：
其中，为主网发电机个数，T为调度时间段，为主网发电机组i在t时刻的开机费用，为主网发电机组 i在t时刻的停机费用，为主网发电机组i在t时刻所产生的煤耗，为主网发电机组i在t时刻所发出的有功功率。
优选地，所述主网的日前调度模型满足的约束条件为：
由主网发电机集合、主网发电机个数、主网节点集合、主网节点个数、主网支路集合、主网支路个数、主网负荷节点集合、主网负荷节点个数、主网连接配电网的节点集合、连接配电网个数、主网负荷功率、主网支路允许通过的最大功率、主网发电机的最大/小出力值、主网发电机的最小开机时间、主网发电机的最小关机时间、主网发电机运行时单位时段内最大上/下爬坡功率值、主网发电机开机后单位时段内最大上爬坡功率值、主网发电机关机前单位时段内最大下爬坡功率值以及主网发电机开机/关机费用构成的功率平衡约束、潮流上限约束、发电机出力上下限约束、发电机状态约束、发电机最小起停时间约束以及发电机爬坡约束。
优选地，在步骤(3)中，所述可行性校验模型包括：
目标函数：以及由配电网m发电机集合、配电网m 发电机个数、配电网m节点集合、配电网m节点个数、配电网m支路集合、配电网m支路个数、配电网m负荷节点集合、配电网m负荷节点个数、配电网m节点的父/子节点集合、配电网m支路电阻/电抗值、配电网m支路允许通过的最大电流值、配电网m节点电压允许最大/小值、配电网m负荷的有功/无功功率以及主网能够向配电网m传输的最大/最小有功功率构成的潮流平衡约束、电压电流上下限约束以及储能功率与能量约束。
优选地，在步骤(6)中，所述由对偶乘子以及目标函数值得到配电网 m向主网反馈的约束为：
其中，表示配电网m在第Ite次的目标函数，Ite表示迭代次数，PTDm,t表示主网向配电网m在t时刻传输的有功功率，为待求解变量，完成求解后得到的值为第Ite+1次的迭代值，表示为第Ite次迭代主网向配电网m在t时刻传输的有功功率，与表示第Ite次时配电网m在t时刻的对偶乘子，T为调度时间段。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：本发明充分考虑了含储能的主动配电网的有功调节能力对日前调度产生的影响，使得在主网做日前调度时充分考虑主动配电网的调节能力，经优化后的调度方案能够实现整个电网的全局最优，最大限度的降低电网发电机组的发电成本，减少资源的浪费。
附图说明
图1是本发明实施例提供的一种case6主网算例拓扑示意图；
图2是本发明实施例提供的一种case69配电网1算例拓扑示意图；
图3是本发明实施例提供的一种case24配电网2算例拓扑示意图；
图4是本发明实施例提供的一种主网与配电网接线示意图；
图5是本发明实施例提供的一种主网与配电网负荷曲线示意图；
图6是本发明实施例提供的一种配电网储能功率及电量优化结果；
图7是本发明实施例提供的一种主网发电机输出功率优化结果；
图8是本发明实施例提供的一种考虑主动配电网有功调节能力的主网日前调度方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提出了一种考虑主动配电网有功调节能力的主网日前调度方法，使得在主网做日前调度时充分考虑主动配电网的调节能力，达到利用配电网的调节能力实现削峰填谷的目的，解决了现有调度方法不能取得全局最优解的弊端，与此同时，基于主网与配网的信息交互提出了一种实现该模型高效求解的方法。
如图8所示为本发明实施例提供的一种考虑主动配电网有功调节能力的主网日前调度方法的流程示意图，在图8所示的方法中具体包括以下步骤：
(1)在主网与各配电网之间没有信息交换时，由主网与各配电网的目标技术参数得到主网的日前调度模型；
其中，主网与各配电网的目标技术参数中，主网的目标技术参数包括：
主网发电机集合Gt，主网发电机个数主网节点集合Bt，主网节点个数主网支路集合Lt，主网支路个数主网负荷节点集合LDt，主网负荷节点个数主网连接配电网的节点集合TD，连接配电网个数Nd，主网负荷i在t时刻功率主网支路i允许通过的最大功率主网发电机i最大/小出力值主网发电机i最小开机时间主网发电机i最小关机时间主网发电机i运行时单位时段内最大上/下爬坡功率值主网发电机i开机后单位时段内最大上爬坡功率值主网发电机i关机前单位时段内最大下爬坡功率值主网发电机i开机/ 关机费用
配电网的目标技术参数包括：配电网m发电机集合Gd,m，配电网m发电机个数配电网m节点集合Bd,m，配电网m节点个数配电网m支路集合Ldis,m，配电网m支路个数配电网m负荷节点集合LDd,m，配电网m负荷节点个数配电网m节点j的父/子节点集合πd,m(j)/δd,m(j)，配电网m支路(i,j)电阻/电抗值配电网m支路(i,j)允许通过的最大电流值配电网m节点j电压允许最大/小值配电网m负荷j在t时刻有功/无功功率主网能够向配网m传输的最大/最小有功功率
其中电网负荷预测曲线可以从调度中心获取。
(2)对主网的日前调度模型进行求解得到由主网传输至各配电网的有功功率；
(3)对于任意配电网m，由引入的正向松弛变量与负向松弛变量构造可行性校验模型，并对可行性校验模型进行求解得到配电网m的目标函数值以及配电网m中储能的充放电状态值；
(4)若目标函数值为0，则表示由主网向配电网m传输的有功功率满足配电网m的功率需求，并执行步骤(5)，若目标函数值不为0，则表示由主网向配电网m传输的有功功率不满足配电网m的功率需求，并执行步骤(6)；
(5)主网向配电网m输出调度方案；
(6)将配电网m中储能的充放电状态值替换可行性校验模型中原有的充放电状态值，并根据替换后的可行性校验模型得到对偶乘子，由对偶乘子以及目标函数值得到配电网m向主网反馈的约束，利用配电网m向主网反馈的约束对主网的日前调度模型进行重新求解得到主网传输至配电网m 的新的有功功率，并执行步骤(3)。
其中，在步骤(6)中，所述由对偶乘子以及目标函数值得到配电网m 向主网反馈的约束为：
其中，表示配电网m在第Ite次的目标函数，Ite表示迭代次数，PTDm,t表示主网向配电网m在t时刻传输的有功功率，为待求解变量，完成求解后得到的值为第Ite+1次的迭代值，表示为第Ite次迭代主网向配电网m在t时刻传输的有功功率，与表示第Ite次时配电网m在t时刻的对偶乘子，T为调度时间段。
在一个可选的实施方式中，在步骤(1)中，所述主网的日前调度模型为：
其中，为主网发电机个数，T为调度时间段，为主网发电机组i 在t时刻的开机费用，为主网发电机组i在t时刻的停机费用，为主网发电机组i在t时刻所产生的煤耗，为主网发电机组i在t时刻所发出的有功功率。
其中，各变量所需满足的约束如下所示：
其中，为主网发电机i在t时刻开机动作，为主网发电机i在t时刻关机动作，均为0/1变量，以及分别表示为煤耗成本的常数项、一次项系数以及二次项系数。对机组费用做线性化处理可得：
式中：AGi,p/BGi,p为机组i的第p分段对应的燃料系数。
其中，主网的日前调度模型满足的约束条件为：由主网发电机集合、主网发电机个数、主网节点集合、主网节点个数、主网支路集合、主网支路个数、主网负荷节点集合、主网负荷节点个数、主网连接配电网的节点集合、连接配电网个数、主网负荷功率、主网支路允许通过的最大功率、主网发电机的最大/小出力值、主网发电机的最小开机时间、主网发电机的最小关机时间、主网发电机运行时单位时段内最大上/下爬坡功率值、主网发电机开机后单位时段内最大上爬坡功率值、主网发电机关机前单位时段内最大下爬坡功率值以及主网发电机开机/关机费用构成的功率平衡约束、潮流上限约束、发电机出力上下限约束、发电机状态约束、发电机最小起停时间约束以及发电机爬坡约束。
具体的，日前调度模型满足以下约束：
其中，约束(5)为功率平衡约束，(6)为潮流上限约束，(7)为发电机出力上下限约束，(8)(9)为发电机状态约束，(10)～(13)为发电机最小起停时间约束，(14)(15)为发电机爬坡约束。其中，PTDm,t为主网向配网m在t时刻传输的有功功率；均为0/1变量，主网发电机i在t时刻开机动作，0为在t时刻不动作，1为在t时刻开机；主网发电机i在t时刻关机动作，0为在t时刻不动作，1为在t时刻关机；为主网发电机i在t时刻状态，0为在t时刻处于关机状态，1为在t时刻处于开机状态；Nd以及分别为负荷节点j在t时刻的有功功率、配电网个数、在第Ite次迭代中由主网传送至配电网的有功功率，b为参考节点、为直流潮流中的传输矩阵(由网架结构决定)，Gt(b)为连接在b节点上的发电机节点集合，LDt(b)为连接在b节点上的负荷节点集合，为负荷节点j在 t时刻的有功功率，TD(b)为连接在b节点上的配电网集合，PTDk,t为主网传送至配电网的有功功率，为线路允许通过的最大功率值，为主网发电机i在t时刻关机动作，为主网发电机i在t时刻开机动作，为主网发电机i在t时刻状态，为主网发电机i在t+1时刻状态，为主网发电机i在 w时刻状态，为主网发电机i在t-1时刻状态，为发电机i在t+1时刻的有功出力值。
考虑到主网在后续会与配电网进行信息交换，在做本次优化时，记迭代次数Ite＝1，在后续迭代中，上述变量的数据会不断更新。
在主网完成主网的日前调度模型的求解后，会向下属的各个配电网传递信息——“主网向配电网传输的功率值”，由于主网的信息量有限，在对主网进行优化时未考虑配电网的功率需求，因此在配电网优化时需要验证由主网传输来的“信息”是否可行，因此称本次优化为配电网m的可行性校验。在下述优化中，为了保证优化有解，引入了正向松弛变量与负向松弛变量两个松弛变量均为正数。
另外为了方便描述，认为主网与配电网的连接点为配电网的1号节点，且认为配电网中不含可控的发电机组，即认为配电网的有功功率均由主网输送，配电网的无功需求由配电网变电站内的无功补偿满足，1号节点处的电压值在存在变电站变压器的情况下认为是可调节的。基于此假设，对于任意配电网m，在步骤(3)中，构建的可行性校验模型包括：目标函数：
以及由配电网m发电机集合、配电网m发电机个数、配电网m节点集合、配电网m节点个数、配电网m支路集合、配电网m支路个数、配电网 m负荷节点集合、配电网m负荷节点个数、配电网m节点的父/子节点集合、配电网m支路电阻/电抗值、配电网m支路允许通过的最大电流值、配电网 m节点电压允许最大/小值、配电网m负荷的有功/无功功率以及主网能够向配电网m传输的最大/最小有功功率构成的潮流平衡约束、电压电流上下限约束以及储能功率与能量约束。
具体地，各约束为：
其中，约束(17)～(24)为潮流平衡约束，约束(25)(26)为电压电流上下限约束，约束(27)～(32)为储能功率与能量约束。其中为配网m发电机j在t时刻有功/无功功率，为配网m支路(i,j)在t时刻传输的有功/无功功率，为配网m节点j在t时刻电压值，为配网m支路(i,j)在t时刻通过电流值，为配网m中储能i在t时刻功率/电能，为配网m中储能i在t时刻放电功率/充电功率，为配网m中储能i在t时刻放电状态/充电状态，为0/1变量，当为1 时代表储能处于放电/充电状态，为0时代表此时储能不放电/充电。ESd,m(j) 为连接在j节点上的储能节点的集合，为储能发出的有功功率，为线路j～k上流过的有功功率，为线路i～j电阻，为线路j～k上流过的无功功率，为线路i～j电抗，为1节点处的负荷有功功率，δd,m(1)为1 节点的子节点集合，为线路j～k上流过的有功功率，为1节点处的无功补偿器的无功功率，为1节点处的负荷无功功率，为线路j～k上流过的无功功率，为i节点在t时刻的电压值，Ld,m为第m个配电网的线路集合，ESd,m为储能节点的集合，分别为i节点处的储能放电功率的最小值与最大值，分别为i节点处的储能充电功率的最小值与最大值，分别为i节点处的储能电量的最小值与最大值，为i节点处的储能在t-1时刻的电量值，分别为i节点处的储能放电、充电状态。
需要注意的是，当存在Nd个配电网时，需要求解Nd个优化模型，分别记为子问题1～子问题Nd。
在本发明实施例中，若目标函数值不为0，则表示由主网向配电网m 传输的有功功率不满足配电网m的功率需求，配电网需要向主网进行信息的交换(即反馈)，即告诉主网需要将某时刻的传输功率增大或减少，反映在数学形式上即为对传输功率PTDm,t的一条数学约束。
假设目标函数不为0的可行性校验问题为第k个子问题(即第k个配电网的可行性校验问题不满足)，则需要向主网返回相应的约束。约束形成过程如下：
首先将储能在T时段内的出力状态保存，记为(上标Ite表示迭代次数)，再次求解(16)～(32)，在本次优化中储能的充放电状态不再是0/1变量，而是用代替。完成求解后，采用YALMIP工具包中的dual函数求解约束(19)(20)所对应的对偶乘子，记为记本次优化得到的目标函数为完成对偶乘子与目标函数的求解后，该配电网向主网反馈的约束形式如下所示：
其中，PTDk,t表示主网向配网k在t时刻传输的有功功率，为待求解变量，完成求解后得到的值为第Ite+1次的迭代值。，表示为第Ite次迭代主网向配网k在t时刻传输的有功功率。完成约束数学表达式的求解后，再次返回步骤(2)，令迭代次数Ite加1，并求解优化问题(1)～(15)(33)。
需要注意的是，若在步骤(6)中有Nk个子问题的目标函数不为0，则子问题向主问题返回的约束共有Nk条，每一条约束均按照以上约束生成的方式构造。
以下结合附图及实施例对发明方法进行进一步详细说明。
1、收集主网与配电网m常规技术参数，获取该电网次日负荷预测曲线。所述主网与Nd个配网的常规技术参数。主网常规的技术参数包括：
主网发电机集合Gt，主网发电机个数主网节点集合Bt，主网节点个数主网支路集合Lt，主网支路个数主网负荷节点集合LDt，主网负荷节点个数主网连接配电网的节点集合TD，连接配电网个数Nd，主网负荷i在t时刻功率主网支路i允许通过的最大功率主网发电机i最大/小出力值主网发电机i最小开机时间主网发电机i最小关机时间主网发电机i运行时单位时段内最大上/下爬坡功率值主网发电机i开机后单位时段内最大上爬坡功率值主网发电机i关机前单位时段内最大下爬坡功率值主网发电机i开机/ 关机费用
配电网常规的技术参数包括：配网m发电机集合Gd,m，配网m发电机个数配网m节点集合Bd,m，配网m节点个数配网m支路集合Ldis,m，配网m支路个数配网m负荷节点集合LDd,m，配网m负荷节点个数配网m节点j的父/子节点集合πd,m(j)/δd,m(j)，配网m支路(i,j)电阻/ 电抗值配网m支路(i,j)允许通过的最大电流值配网m节点j电压允许最大/小值配网m负荷j在t时刻有功/无功功率主网能够向配网m传输的最大/最小有功功率
其中所述电网负荷预测曲线可以从调度中心获取，所述电网的负荷曲线如图5所示。
在本实例中，主网算例采用6节点标准算例，如图1所示是本发明实施例提供的一种case6主网算例拓扑示意图，下面接入两个配电网(Nd＝2， 69节点标准算例和24节点算例，分别记为配电网1和2，主网算例与配电网算例连接示意图如图4所示。
在主网中，包含发电机三台，发电机参数如下表1所示：
表1
| 编号 | 出力上限 | 出力下限 | 上爬坡速度 | 下爬坡速度 | 最小开机时间 | 最小关机时间 |
| 1 | 100MW | 50MW | 60MW/h | 60MW/h | 6h | 6h |
| 2 | 75MW | 18.75MW | 45MW/h | 45MW/h | 6h | 6h |
| 3 | 90MW | 22.5MW | 54MW/h | 54MW/h | 6h | 6h |
在配电网1(case69算例)中，含有一台储能位于12号节点，如图2 所示是本发明实施例提供的一种case69配电网1算例拓扑示意图，其参数如下表所示(编号1)，配电网2(case24算例)中，含有一台储能位于6 节点处，如图3所示是本发明实施例提供的一种case24配电网2算例拓扑示意图。其储能参数如下表2所示(编号2)。
表2
2、初始化迭代次数Ite＝1，计算日前调度模型(1)～(15)，计算所得得到第一次迭代传输功率如下表3所示。
表3
3、将Ite＝1带入到优化(16)～(32) 中，求解该优化，得到第一次迭代松弛变量的值如下表4所示：
表4
两个子问题的目标函数分别为
4、完成步骤3的求解后，两个子问题的目标函数均不为0，转至步骤 5。
5、两个配电网子问题的目标函数均不为0，因此两个配电网均需要向主网返回约束，步骤如下：首先确定每个时刻储能的动作状态，得到第一次迭代储能充放电状态如下表5所示：
表5
将上表5保存，并再次求解优化(16)～(32)，在求解过程中，将优化中的变量用上表5代替，完成求解后，采用YALMIP工具包中的dual函数求解约束(19)(20)的对偶乘子，可以得到第一次迭代对偶乘子如下表6：
表6
保存本次求解中的目标函数，并向主网中返回约束：
同时，令Ite＝2，重新返回步骤2。
6、重新计算日前调度模型(1)～(15)，计算所得得到第二次迭代传输功率如下表7所示：
表7
7、将上表带入至优化(16)～(32)中，求解该优化，得到第二次迭代松弛变量的值如下表8所示：
表8
两个子问题的目标函数分别为
8、完成步骤3的求解后，两个子问题的目标函数均不为0，转至步骤 9。
9、两个配电网子问题的目标函数均不为0，因此两个配电网均需要向主网返回约束，步骤如下：首先确定每个时刻储能的动作状态，得到第二次迭代储能充放电状态如下表9所示：
表9
将上表保存，并再次求解优化(16)～(32)，在求解过程中，将优化中的变量用上表代替，完成求解后，采用YALMIP工具包中的dual函数求解约束(19)(20)的对偶乘子，可以得到第二次迭代对偶乘子下表10所示：
表10
保存本次求解中的目标函数，并向主问题中返回约束：
同时，令Ite＝3，重新返回步骤2。
10、重新计算日前调度模型(1)～(15)，计算所得得到第三次迭代传输功率如下表11所示：
表11
11、将上表带入至优化(16)～(32)中，求解该优化，得到第三次迭代松弛变量的值如下表12所示：
表12
两个子问题的目标函数分别为
12、完成步骤3的求解后，两个子问题的目标函数均不为0，转至步骤 13。
13、两个配电网子问题的目标函数均不为0，因此两个配电网均需要向主网返回约束，步骤如下：首先确定每个时刻储能的动作状态，得到第三次迭代储能充放电状态如下表13所示：
表13
将上表保存，并再次求解优化(16)～(32)，在求解过程中，将优化中的变量用上表代替，完成求解后，采用YALMIP工具包中的dual函数求解约束(19)(20)的对偶乘子，可以得到第三次迭代对偶乘子如下表14：
表14
保存本次求解中的目标函数，并向主网中返回约束：
同时，令Ite＝4，重新返回步骤2。
14、重新计算日前调度模型(1)～(15)，计算所得得到第四次迭代传输功率结果如下表15所示：
表15
15、将上表带入至优化(16)～(32)中，求解该优化，得到第四次迭代松弛变量的值如下表16所示：
表16
两个子问题的目标函数分别为
16、完成步骤3的求解后，在误差允许范围内可以认为两个子问题的目标函数为0，可以输出此时的调度方案，机组开机方式如下表17所示：
表17
在最终结果中，储能的各时段的充放电功率值与电量值如图6所示，各时段发电机组的输出功率如图7所示。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种考虑主动配电网有功调节能力的主网日前调度方法，其特征在于，包括：
(1)在主网与各配电网之间没有信息交换时，由主网与各配电网的目标技术参数得到主网的日前调度模型；
(2)对主网的日前调度模型进行求解得到由主网传输至各配电网的有功功率；
(3)对于任意配电网m，由引入的正向松弛变量与负向松弛变量构造可行性校验模型，并对可行性校验模型进行求解得到配电网m的目标函数值以及配电网m中储能的充放电状态值；
(4)若目标函数值为0，则表示由主网向配电网m传输的有功功率满足配电网m的功率需求，并执行步骤(5)，若目标函数值不为0，则表示由主网向配电网m传输的有功功率不满足配电网m的功率需求，并执行步骤(6)；
(5)主网向配电网m输出调度方案；
(6)将配电网m中储能的充放电状态值替换可行性校验模型中原有的充放电状态值，并根据替换后的可行性校验模型得到对偶乘子，由对偶乘子以及目标函数值得到配电网m向主网反馈的约束，利用配电网m向主网反馈的约束对主网的日前调度模型进行重新求解得到主网传输至配电网m的新的有功功率，并返回执行步骤(3)；
在步骤(1)中，所述主网的日前调度模型为：
其中，为主网发电机个数，T为调度时间段，为主网发电机组i在t时刻的开机费用，为主网发电机组i在t时刻的停机费用，为主网发电机组i在t时刻所产生的煤耗，为主网发电机组i在t时刻所发出的有功功率。
2.根据权利要求1所述的方法，其特征在于，所述主网的日前调度模型满足的约束条件为：
由主网发电机集合、主网发电机个数、主网节点集合、主网节点个数、主网支路集合、主网支路个数、主网负荷节点集合、主网负荷节点个数、主网连接配电网的节点集合、连接配电网个数、主网负荷功率、主网支路允许通过的最大功率、主网发电机的最大/小出力值、主网发电机的最小开机时间、主网发电机的最小关机时间、主网发电机运行时单位时段内最大上/下爬坡功率值、主网发电机开机后单位时段内最大上爬坡功率值、主网发电机关机前单位时段内最大下爬坡功率值以及主网发电机开机/关机费用构成的功率平衡约束、潮流上限约束、发电机出力上下限约束、发电机状态约束、发电机最小起停时间约束以及发电机爬坡约束。
3.根据权利要求2所述的方法，其特征在于，在步骤(3)中，所述可行性校验模型包括：
目标函数：以及由配电网m发电机集合、配电网m发电机个数、配电网m节点集合、配电网m节点个数、配电网m支路集合、配电网m支路个数、配电网m负荷节点集合、配电网m负荷节点个数、配电网m节点的父/子节点集合、配电网m支路电阻/电抗值、配电网m支路允许通过的最大电流值、配电网m节点电压允许最大/小值、配电网m负荷的有功/无功功率以及主网能够向配电网m传输的最大/最小有功功率构成的潮流平衡约束、电压电流上下限约束以及储能功率与能量约束。
4.根据权利要求3所述的方法，其特征在于，在步骤(6)中，所述由对偶乘子以及目标函数值得到配电网m向主网反馈的约束为：
其中，表示配电网m在第Ite次的目标函数，Ite表示迭代次数，PTDm,t表示主网向配电网m在t时刻传输的有功功率，为待求解变量，完成求解后得到的值为第Ite+1次的迭代值，表示为第Ite次迭代主网向配电网m在t时刻传输的有功功率，与表示第Ite次时配电网m在t时刻的对偶乘子，T为调度时间段。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710927881.5A CN107591806B (zh) | 2017-10-09 | 2017-10-09 | 一种考虑主动配电网有功调节能力的主网日前调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710927881.5A CN107591806B (zh) | 2017-10-09 | 2017-10-09 | 一种考虑主动配电网有功调节能力的主网日前调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107591806A CN107591806A (zh) | 2018-01-16 |
| CN107591806B true CN107591806B (zh) | 2019-10-25 |
Family
ID=61052883
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710927881.5A Active CN107591806B (zh) | 2017-10-09 | 2017-10-09 | 一种考虑主动配电网有功调节能力的主网日前调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107591806B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109165865B (zh) * | 2018-09-07 | 2021-07-06 | 国电南瑞科技股份有限公司 | 基于边界开关连通性分析的主配一体化模型校验方法及系统 |
| CN109586299B (zh) * | 2018-12-18 | 2023-11-03 | 中国电力科学研究院有限公司 | 一种配电网有功功率分区优化控制方法及系统 |
| CN111899501B (zh) * | 2020-08-04 | 2021-08-06 | 广东电网有限责任公司 | 一种配网自动化主站变电站内开关的遥控方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102694391B (zh) * | 2012-05-31 | 2014-09-03 | 国电南瑞科技股份有限公司 | 风光储联合发电系统日前优化调度方法 |
| CN105703369B (zh) * | 2016-02-04 | 2020-06-09 | 马瑞 | 一种多能耦合输配网最优能量流建模及求解方法 |
-
2017
- 2017-10-09 CN CN201710927881.5A patent/CN107591806B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107591806A (zh) | 2018-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Sizing a stand-alone solar-wind-hydrogen energy system using weather forecasting and a hybrid search optimization algorithm | |
| Armghan et al. | Nonlinear integral backstepping based control of a DC microgrid with renewable generation and energy storage systems | |
| Lu et al. | Short-term scheduling of battery in a grid-connected PV/battery system | |
| Machlev et al. | A review of optimal control methods for energy storage systems-energy trading, energy balancing and electric vehicles | |
| CN107330568A (zh) | 基于Benders解耦的储能、分布式电源与配电网协调规划方法 | |
| CN107591806B (zh) | 一种考虑主动配电网有功调节能力的主网日前调度方法 | |
| CN110289622B (zh) | 一种光储充能量路由器的日前经济优化调度方法 | |
| CN103824123A (zh) | 一种新的配电网蓄电池储能系统的优化配置算法 | |
| CN112202205A (zh) | 一种多能源三级自律协同控制方法及装置 | |
| CN111555370A (zh) | 一种基于云边协同的配电网分层协调调度方法及装置 | |
| Ezhilarasan et al. | Design and development of energy management system for DG source allocation in a micro grid with energy storage system | |
| Haddad et al. | Stand-alone renewable power system scheduling for a green data center using integer linear programming | |
| Raza et al. | Robust nonlinear control of regenerative fuel cell, supercapacitor, battery and wind based direct current microgrid | |
| CN106786550A (zh) | 一种微电网成本优化的分布式控制方法及装置 | |
| CN116865270A (zh) | 一种含嵌入式直流的柔性互联配电网优化调度方法及系统 | |
| CN116805192A (zh) | 考虑最优弃能率的综合能源系统双层规划优化方法及应用 | |
| Alvarez-Mendoza et al. | Semi-dispatchable generation with wind-photovoltaic-fuel cell hybrid system to mitigate frequency disturbance | |
| Phan-Van et al. | A comparison of different metaheuristic optimization algorithms on hydrogen storage-based microgrid sizing | |
| Ignat et al. | Renewable Energy Microgrid Model using MATLAB—Simulink | |
| Samy et al. | Fuzzy logic based battery power management for PV and wind hybrid power system | |
| Yu et al. | A fuzzy Q-learning algorithm for storage optimization in islanding microgrid | |
| CN111082442B (zh) | 一种基于改进fpa的储能容量优化配置方法 | |
| CN114970096A (zh) | 基于概率最优潮流的电力系统灵活性量化评估方法及系统 | |
| Franco et al. | MIMO auto-regressive modeling-based generalized predictive control for grid-connected hybrid systems | |
| Ignat et al. | Islanded microgrid simulation and cost optimisation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190425 Address after: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037 Applicant after: Huazhong University of Science and Technology Applicant after: Power Science Research Institute of Anhui Electric Power Co., Ltd. Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037 Applicant before: Huazhong University of Science and Technology |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |