一种综合能源系统日前滚动优化方法与流程

本发明涉及一种综合能源系统日前滚动优化方法，属于综合能源系统技术领域。

背景技术：

综合能源系统融合了电、热、冷、气等不同类型的能源，为实现多能源的优化互补运行提供可能。采用互联网技术，将以电力资源为主的多种能源互联，实现综合能源系统中的横向多能互补、纵向需求响应联合优化，使能量生产与消耗双向互动。从源端和荷端两方面实现能量协同优化，将柔性负荷进行最佳匹配，对于用户用能经济性和整个综合能源系统的稳定优化运行都具有重要的现实意义。

现有的日前电量计算方法不合理，结算效率低下。

技术实现要素：

本发明要解决的技术问题是：提供一种综合能源系统日前滚动优化方法，以解决现有技术中存在的问题。

本发明采取的技术方案为：一种综合能源系统日前滚动优化方法，该方法包括以下步骤：

(1)获取电负荷、冷负荷、热负荷用能需求；

(2)获取实时电价、上网电价、天然气价；

(3)用户用能计划的确定、分布式电源出力和用户负荷预测；

(4)根据多目标优化指标计算用户一天的总运行成本、不可再生能源消耗量、碳排放减少量、用户舒适度指标；

(5)根据用户偏好确定各目标函数权重；

(6)将各个目标函数融合成单一目标函数，融合后的单一目标函数为：

ming＝λ1f1+λ2φcf2+λ3φef3-λ4φuf4(1)

式中

λ1，λ2，λ3，λ4——用户赋予经济性、环保性、节能性和舒适度目标函数的权系数；

φc——碳排放量经济系数，随co2排放交易市场的成交价格设定；

φe——一次能源经济系数，随一次能源价格设定；

φu——用户舒适度经济系数，随用户舒适度敏感系数设定；

(7)通过分支定界法进行求解，得出用户当天最佳用能计划，通过系统人机交互界面反馈给用户。

分支定界法求解单一目标函数：将整型变量约束松弛为实数连续变量，转化为二次规划qp子问题进行求解，然后得到其中符合整数约束的原单一目标函数的最优解。

qp子问题采取可行方向法求解，qp子问题的一般形式如下：

采用zoutendijk可行方向法求解qp子问题的步骤如下：

(1)取初始可行点x⁽¹⁾，令k＝1；

(2)进行第k次迭代，在x^(k)处将不等式约束分解为起作用约束和不起作用约束：和使得a1x^(k)＝b1，a2x^(k)≤b2，计算▽f(x^(k))＝qx^(k)；

(3)求解下列线性规划问题，得到最优解d^(k)；

(4)如果(qx^(k))^td^(k)＝0，则x^(k)为最优解，停止计算；否则，转至(5)；

(5)求步长λk。令则最大步长为：

求解：得到最优解λk，令x^(k+1)＝x^(k)+λkd^(k)；

(6)令k＝k+1，转至(2)。

本发明的有益效果：与现有技术相比，本发明考虑经济性、环保性、节能性及舒适度的多目标优化问题，采用基于b&b算法和zoutendijk可行方向法的miqp求解算法对问题进行求解，能够根据用户偏好对用户端能量流进行优化，达到经济、节能、环保、用户舒适度较高的综合效益，计算更科学，效率更快。

附图说明

图1为日前计划优化流程图；

图2为3个整型变量约束的miqp的完全二叉树图；

图3为基于b&b算法和zoutendijk可行方向法的miqp求解流程图；

图4为miqp分支定界求解日前计划优化过程与最优值曲线图。

具体实施方式

下面结合附图及具体的实施例对发明进行进一步介绍。

日前计划优化基于上级电网日前下发的实时电价和用户日前预安排的用能计划，由用户设定日前用能计划后触发。

日前计划优化基于分布式能源发电和用户负荷预测等数据，希望在日前多种用能方案中为用户制定最佳的用能计划，以用户为核心，构建综合多能互补与需求响应的联合优化模型，融合经济性、环保性、节能性、舒适性四方面的因素，同时考虑用户偏好。

在前述讨论的多目标优化性能评估指标中，对于经济性指标，运行成本越小表明系统越经济；对于环保性指标，碳排放量越少表明用户环保效益越高；对于节能性指标，不可再生能源消耗量越少表明系统越节能；对于舒适度指标，用户总舒适度越高意味着用户整体舒适度体验越好。

日前计划优化详细的目标函数如下：

1)对于综合能源系统经济性指标，优化目标主要是在满足系统用能需求的前提下，使系统总运行成本最小，目标函数可表示为：

式中

t——时域长度。

2)考虑系统环保性，优化目标尽可能满足最小的碳排放量，目标函数可表示为：

3)考虑系统节能性，优化目标尽可能使不可再生能源消耗量最少，目标函数可表示为：

4)考虑用户舒适度要求，优化目标尽可能使用户整体舒适度最高，目标函数可表示如下：

日前计划优化的约束条件包括：能量流供需平衡约束条件、分布式能源出力约束、负荷需求约束。

能量流供需平衡约束条件包括电平衡、热平衡和冷平衡三部分。

1)电平衡

egrid,t+egt,t+epv,t+ewp,t+eess,t＝eeload,t(8)

2)热平衡

3)冷平衡

分布式能源出力约束考虑光伏发电、风力发电和微型燃气轮机的出力限制。

式中

——光伏发电、风力发电在t时段的最大输出功率(w)；

——微型燃气轮机爬坡率上下限。

冷热负荷需求满足约束在夏季模式下考虑室内空气制冷冷负荷和生活热水热负荷需求，在冬季模式下考虑室内空气制热和生活热水热负荷需求。夏季模式下具体约束为：

冬季模式下室内空气制热热负荷需求具体约束为：

针对冷热负荷需求的约束范围是温度变量，而系统优化变量是空调、电热水器或者三联供出力，因此需要根据优化变量解出温度变量，然后将约束表示为针对优化变量的线性约束。

令λa＝caρava，ta＝[temp1temp2…tempt]^t，在系统优化初始时刻的室内温度temp0为已知量，室外温度曲线可根据天气预报数据获取，则式(12)可表示为：

式中：

令则式(15)可写为：

s^cata＝p^cae^ca-b^ca(16)

解出室内温度序列：

得到标准形式约束：

记

生活热水热负荷的温度约束以及电动汽车、储能电池的能量约束上下限与上述推导过程类似，上述矩阵变换过程可以推导得出其他用能需求的线性约束。

令λw＝cwρw，tw＝[θ1θ2…θt]^t，在系统优化初始时刻的热水温度为θ0。令

则式(13)可表示为：

s^hwtw＝p^hwe^hw+b^hw(19)

解出生活热水温度序列：

得到标准形式约束：

记

令oev＝qev[socev,1socev,2…socev,t]^t，在系统优化初始时刻的电动汽车soc为socev,0。令

则电动汽车的soc约束可表示为：

s^evoev＝p^eve^ev+b^ev(22)

解出生活热水温度序列：

得到标准形式约束：

令电动汽车的充放电功率约束可表示为：

记

令记储能电池的不等式约束为aesse^ess≤bess，储能电池的约束与电动汽车类似，aess，bess的推导过程可参考电动汽车约束推导过程。此外，记冷热电联供系统的不等式约束为acchpe^cchp≤bcchp；ae、ah、ac和be、bh、bc分别为电、热、冷能量平衡的等式约束，详细推导过程本发明在此不再赘述。

考虑经济性、环保性、节能性和用户舒适度的日前计划优化模型是一个多目标规划问题(multi-objectiveprogramming，mop)。多目标优化中的各个子目标往往是相互冲突的，即各个子目标不能同时达到最优解，只能进行折中处理，使各个子目标函数尽可能接近最优。

针对多目标优化问题，其求解算法分为传统优化算法和智能优化算法两类。具体来讲，传统优化算法通过将多目标函数转化为单目标函数进行求解，常用的有线性加权法、约束法和线性规划法等。而智能优化算法包括遗传算法和粒子群算法等。在众多求解算法中，线性加权法是最基本的解法，该方法对单个目标函数赋予一定权重，从而将各个目标函数线性组合成为一个目标函数：

式中

λ——权重向量，λ的各个分量λi称为权系数，为不失一般性，

x——优化变量向量，包括e^ca、e^hw、e^ev、e^ess、e^cchp；

f——目标函数向量，各个分量对应各目标函数值。

由于本发明的优化针对用户端能源的经济性、环保性、节能性和用户舒适度进行综合优化，而用户的个人偏好正好决定了各目标函数的权系数，因此本发明采用线性加权法将多目标优化问题转化为单目标优化问题求解。

实施例1：一种综合能源系统日前滚动优化方法，该方法包括以下步骤：

(1)获取电负荷、冷负荷、热负荷用能需求；

(2)获取实时电价、上网电价、天然气价；

(3)用户用能计划的确定、分布式电源出力和用户负荷预测；

(4)根据多目标优化指标计算用户一天的总运行成本、不可再生能源消耗量、碳排放减少量、用户舒适度指标；

(5)根据用户偏好确定各目标函数权重；

(6)将各个目标函数融合成单一目标函数，融合后的单一目标函数为：

ming＝λ1f1+λ2φcf2+λ3φef3-λ4φuf4(1)

式中

λ1，λ2，λ3，λ4——用户赋予经济性、环保性、节能性和舒适度目标函数的权系数；

φc——碳排放量经济系数，随co2排放交易市场的成交价格设定；

φe——一次能源经济系数，随一次能源价格设定；

φu——用户舒适度经济系数，随用户舒适度敏感系数设定；

(7)通过分支定界法进行求解，得出用户当天最佳用能计划，通过系统人机交互界面反馈给用户。

根据以上日前计划优化模型，用户端综合能源系统能量日前计划优化流程图如图1所示，实行日前计划优化需要完成电网数据的采集、用户用能计划的确定、分布式电源出力和用户负荷预测等数据，随后根据多目标优化指标计算用户一天的总运行成本、不可再生能源消耗量、碳排放减少量、用户舒适度指标，然后根据用户偏好确定各目标函数权重，将各个目标函数融合成单一目标函数，求解得出用户当天最佳用能计划，通过系统人机交互界面反馈给用户。

由于引入了电动汽车充放电标志σev、储能电池充放电标志σess，充放电标志σev、σess均为{0,1}型整数变量，使得整个优化问题混入了整数约束，成为混合整数规划问题(mixedintegerprogramming，mip)。各子目标函数中，经济性、环保性、节能性目标函数都是线性的，只有舒适度目标函数为二次目标函数。因此日前计划优化问题为混合整数二次规划(mixedintegerquadraticprogramming，miqp)问题，可以通过分支定界法(branch&bound，b&b)进行求解。

mip问题的决策变量中既包含了整型变量，也包含连续变量，一般来讲，mip问题是np-hard问题，算法时间复杂度为指数级。目前针对混合整数非线性规划问题的有效求解方法有分支定界法(branch&bound，b&b)、广义bender分解法(generalizedbendersdecomposition，gbd)、外逼近法(outerapproximation，oa)、遗传算法(genericalgorithm，ga)、粒子群算法以及模拟退火算法等。但miqp作为特殊的混合整数非线性规划问题，具有如下特点：二次规划问题相对容易求解，当加权矩阵为正定矩阵时，优化问题有唯一最优解，因此，在问题规模小于一定范围内，可以基于分支定界法有效求解miqp问题。

本发明的miqp问题可以写作如下形式：

其中：

通过前文建模过程可以发现，本发明中miqp问题具有两大特点：加权矩阵q为半正定矩阵；整型变量的取值为0或者1(整型变量上下限约束为1和0)。因此本发明讨论的b&b算法针对此类特殊问题的求解更高效。

分支定界法求解单一目标函数：将整型变量约束松弛为实数连续变量，转化为二次规划qp子问题进行求解，然后得到其中符合整数约束的原单一目标函数的最优解。

qp子问题采取可行方向法(如zoutendijk可行方向法、frank-wolfe方法等)求解，qp子问题的一般形式如下：

采用zoutendijk可行方向法求解qp子问题的步骤如下：

(1)取初始可行点x⁽¹⁾，令k＝1；

(2)进行第k次迭代，在x^(k)处将不等式约束分解为起作用约束和不起作用约束：和使得a1x^(k)＝b1，a2x^(k)≤b2，计算▽f(x^(k))＝qx^(k)；

(3)求解下列线性规划问题，得到最优解d^(k)；

(4)如果(qx^(k))^td^(k)＝0，则x^(k)为最优解，停止计算；否则，转至(5)；

(5)求步长λk。令则最大步长为：

求解：得到最优解λk，令x^(k+1)＝x^(k)+λkd^(k)；

(6)令k＝k+1，转至(2)。

基于b&b求解本发明中miqp问题的原理可以采用完全二叉树表示。设一组向量ξ∈{0,1,#}^|ω|，该向量维数为miqp问题中整型变量的个数，向量中元素值与各个整型变量的取值一一对应，其中“#”表示[0,1]区间内的任意实数。一个向量ξ对应了一个qp问题，ξ中元素的取值对应了qp问题中整型变量的取值，同时ξ向量也对应二叉树中的一个节点，二叉树的根节点对应ξ0＝[#,…#]，表示对miqp问题中的所有整型变量都松弛为[0,1]。通过将ξ向量中的某一元素设定值为0或1，可以在ξ向量基础上生成新的qp问题，例如ξ0向量可以生成ξ1＝[0,#,…#]和ξ2＝[1,#,…#]向量，ξ1和ξ2在二叉树中对应了ξ0向量节点的子节点，以ξ1向量为例，其对用的qp问题为：

上述qp问题可以转化为式(3-78)表示的一般qp问题采用zoutendijk可行方向法进行求解。具有3个整型变量约束的miqp的完全二叉树如图2-3所示，图中每个节点对应一个ξ向量，也对应一个qp问题。

根据以上二叉树分析，求解含n个整型变量的miqp问题最多需要求解2ⁿ⁺¹-1个qp子问题，算法复杂度为指数级，但在问题规模较小时，仍可以在短时间内求得最优解。

基于b&b算法和zoutendijk可行方向法求解miqp问题具体流程如图3所示。从图中可以看出，当某一个qp问题的最优解大于当前最优解时，其对应的子qp问题将无需求解，这也大大缩减了miqp问题的求解时间。采用图3所示的b&b和zoutendijk可行方向法miqp求解流程求解日前计划优化问题，随着求解子qp问题次数的增加，当前最优值不断更新，最优值随求解次数更新的曲线如图4所示，最终求出问题的最优解。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内，因此，本发明的保护范围应以所述权利要求的保护范围为准。