CN112665159B - 基于需求响应的负荷反弹量优化、负荷调节方法及*** - Google Patents
基于需求响应的负荷反弹量优化、负荷调节方法及*** Download PDFInfo
- Publication number
- CN112665159B CN112665159B CN202110016757.XA CN202110016757A CN112665159B CN 112665159 B CN112665159 B CN 112665159B CN 202110016757 A CN202110016757 A CN 202110016757A CN 112665159 B CN112665159 B CN 112665159B
- Authority
- CN
- China
- Prior art keywords
- load
- air conditioner
- model
- rebound
- indoor
- 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
Images
Landscapes
- Air Conditioning Control Device (AREA)
Abstract
本发明属于电力需求侧管理技术领域,公开了一种基于需求响应的负荷反弹量优化、负荷调节方法及***,在传统需求响应直接负荷控制方法上增加了优化层,优化后的参数以负反馈的形式返回直接负荷控制层抑制负荷反弹量;再利用迭代学***滑,使得负荷反弹量进一步降低,增强了***的抗干扰能力。从而可以使电网运行安全稳定。
Description
技术领域
本发明属于电力需求侧管理技术领域,具体涉及一种基于需求响应的负荷反弹量优化、负荷调节方法及***。
背景技术
随着需求的增加,各种新型发电方式主导的间歇性新能源大规模地接入到了电网之中,这大大增加了电力负荷的波动,也增加了电力潮流的不稳定性,对电网***稳健运行造成了巨大的威胁。同时,这种情况也促使着人们对电网的自适应能力提出更高的要求。另一方面,虽然电网整体的供需相对比较平衡,不会出现大规模崩溃现象,但是在某些地区的电力负荷高峰期,仍然会出现阶段性、局部性的紧张状态。
在实际需求响应负荷调控过程中,当用户负荷结束受控,再一次转换到运行状态时,出现的负荷反弹对用电峰值的调整起到了很大的影响作用,甚至可以造成新一轮的负荷高峰,造成更大的威胁,同时也对峰值负荷的削减提出了新的约束。
发明内容
本发明的目的在于提供一种基于需求响应的负荷反弹量优化、负荷调节方法及***,用以解决现有技术中的响应控制技术在控制期结束后出现的负荷反弹问题等问题。
为了实现上述任务,本发明采用以下技术方案:
基于需求响应的负荷反弹量优化方法,包括如下步骤:
步骤1:建立室内空调负荷模型;
步骤2:采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
步骤3:建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和步骤2得到的负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
步骤4:将步骤3得到的室内空调负荷反弹量通过PD控制算法进行优化,获得优化后的负荷反弹量。
进一步的,步骤1中先通过一阶等效热参数法建立单体空调负荷模型,然后通过蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
进一步的,步骤1建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
进一步的,步骤2中负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
进一步的,步骤3将室内温度上限、控制时长、参与负荷数作为约束条件包括:
控制时长的约束条件为式Ⅱ和式Ⅲ:
τon+τoff≤τ 式Ⅲ
室内温度上限的约束条件为式Ⅳ:
T(τoff)≤Tmax 式Ⅳ
其中,T(τoff)表示每个控制周期的停机期末室内温度,单位为℃;Tmax为当前控制周期内允许的室内最高温度上限,单位为℃。
一种负荷调节方法,用于抑制在需求响应控制时期出现的负荷反弹现象,包括如下步骤:
步骤A:根据负荷反弹量优化方法,获得优化后的负荷反弹量;
步骤B:根据获得的优化后的负荷反弹量进行负荷控制和调节。
基于需求响应的负荷反弹量优化***,包括模型建立模块、负荷控制模块、室内空调负荷反弹量计算模块和优化模块;
所述模型建立模块用于建立室内空调负荷模型;
所述负荷控制模块采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
所述室内空调负荷反弹量计算模块用于建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
所述优化模块用于将室内空调负荷反弹量计算模块得到的室内空调负荷反弹量通过PD控制算法进行优化,获得优化后的负荷反弹量。
进一步的,模型建立模块中先通过一阶等效热参数法建立单体空调负荷模型,然后通过蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
进一步的,模型建立模块建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
进一步的,负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
本发明与现有技术相比具有以下技术特点:
(1)本发明是通过MATLAB对建立的模型进行直接负荷控制不同参数下的仿真模拟,获取运行参数以及需求响应控制结束之后负荷反弹量数据。通过对数据进行分析处理,确定了直接负荷控制过程中直接影响负荷反弹量的三个指标:直接负荷控制温度上限、直接负荷控制时长和参与负荷数量。通过对这三个指标和负荷反弹量进行拟合,获取需求响应负荷反弹指标和负荷反弹量之间的关系,从而更好的描述了负荷反弹;
(2)本发明以直接负荷控制温度上限、直接负荷控制时长和参与负荷数量这三个指标所引起的负荷反弹量总和最小为目标,以空调最小启停时间、控制周期长度以及温度上限为约束,对传统直接负荷控制策略进行改进,提出基于直接负荷控制的负荷反弹抑制优化策略,以弥补传统控制策略在负荷反弹方面的不足;
(3)本发明针对基于直接负荷控制的负荷反弹抑制优化策略优化后的数据通过迭代学***滑,提出基于迭代学***滑的负荷曲线,使冷机、空调机组等设备逐步转变工作状态,可进一步减少***负荷反弹,可为后续需求响应工作更好地展开提供参考依据。
附图说明
图1为本发明的实施例中空调一阶等效热参数模型;
图2为本发明的实施例中室外温度逐时变化曲线图;
图3为本发明的实施例中一天内用户组1-10的自然负荷变化情况图;
图4为本发明的实施例中用户组6在执行DLC控制后与自然状态全天负荷对比图;
图5为本发明的实施例中控制前后负荷对比;
图6为本发明的实施例中采用PD型迭代学习控制算法优化后负荷曲线图;
图7为迭代学习算法示意图。
具体实施方式
以下给出本发明的具体实施例,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。以下各实施例的说明是参考附图,用以示例本发明可以用以实施的特定实施例。
一阶等效热参数方法:负荷建模是负荷参与需求响应控制的基础,主要描述的是建筑物室内温度与制冷(热)量之间的关系,是需求响应策略的依据。是需求响应策略的依据。如图1所示为现在比较常用的空调***一阶等效热参数模型。由图1可以得出:
其中,Tin表示空调所处建筑物内的室内温度,℃;Tout表示室外温度,℃;R表示房间的等效热阻,℃/kW;C表示房间的等效热容,kJ/℃;Qac表示空调的制冷量,kW。
NSGA2算法:即带有精英保留策略的快速非支配多目标优化算法,是一种基于Pareto最优解的多目标优化算法。
PD控制算法:一种比例(P)微分(D)控制算法,算法原理如图7所示。
在本实施例中公开了一种负荷反弹量优化方法,包括如下步骤:
步骤1:建立室内空调负荷模型;
步骤2:采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
本步骤进行模拟时,是对同一工况下不同用户群的数据分别进行模拟,每组用户群的负荷数量和大小均不同,获得多组运行参数,对所有的运行参数进行参数拟合后,得到拟合方程,该拟合方程适用于该同一工况下的任何用户;在本实施例中,工况为夏季用户活跃时间(09:00-21:00),本方法也适用于其他季节的工况;
步骤3:建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和步骤2得到的负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
步骤4:将步骤3得到的室内空调负荷反弹量通过PD控制算法进行优化,获得优化后的负荷反弹量。
具体的,建立出的空调负荷模型包括房间温度变化模型、空调启停模型和空调***制冷量模型,如下所示:
1)房间温度变化模型
其中,Δt表示时间间隔,min;s(t)表示空调的运行状态,等于1的时候表示开启状态,等于0表示关闭状态。
2)空调启停模型
空调房间内温度控制上下限分别为:
式中,Tmax、Tmin分别表示室内允许的温度范围上下限,℃;δ表示温度调节参数,℃,一般取2℃。
当室温波动范围为[Tmin,Tmax],可求得启停周期中停机期和制冷期的时长分别为:
式中,τon、τoff分别表示启停周期内制冷期和停机期的时长,min。
3)空调***制冷量模型
根据能量守恒定律,一段时间内,空调***制冷量应该等于围护结构冷负荷、室内热源散热引起的冷负荷、人体散热引起的冷负荷及新风负荷之和;
Qc=Qe+Qh+Qp+Qn
其中:Qc表示空调***制冷量,W;Qe表示围护结构的冷负荷,W;Qh表示室内热源散热引起的冷负荷,W;Qp表示人体散热形成的冷负荷,W;Qn表示新风负荷,W。
具体的,步骤1中先通过一阶等效热参数法建立单体空调负荷模型,然后通过能量守恒定律以及蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
空调负荷模型中的额定功率P、房间等效热阻R、房间等效热容C均服从对数正态分布。假设某区域中有M台空调,每台空调的额定运行功率为P,假设所有的空调均为定频空调,每台空调的运行状态为s(t),根据大数定律,该区域中空调负荷聚合模型即室内空调负荷模型的实时运行状态及功率如下。
具体的,步骤1建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
具体的,步骤2中负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
具体的,整体降负荷目标函数J2表示为:
J2=min(max(LDLC,i,j))
其中,LDLC,i,j表示在i时段,j组用户的总体负荷。该目标函数的含义为:降低第i时段内第j组用户的总体负荷量的最大值。
具体的,步骤3将室内温度上限、控制时长、参与负荷数作为约束条件包括:
控制时长的约束条件为式Ⅱ和式Ⅲ:
τon+τoff≤τ 式Ⅲ
室内温度上限的约束条件为式Ⅳ:
T(τoff)≤Tmax 式Ⅳ
其中,T(τoff)表示每个控制周期的停机期末室内温度,单位为℃;Tmax为当前控制周期内允许的室内最高温度上限,单位为℃。
本实施例中还公开了一种基于需求响应的负荷反弹抑制方法,用于抑制在需求响应控制时期出现的负荷反弹现象,包括如下步骤:
步骤A:根据任一种负荷反弹量优化方法,获得优化后的负荷反弹量;
步骤B:根据获得的优化后的负荷反弹量进行负荷控制和调节。
本实施例中还公开了一种负荷反弹量优化***,包括模型建立模块、负荷控制模块、室内空调负荷反弹量计算模块和优化模块;
所述模型建立模块用于建立室内空调负荷模型;
所述负荷控制模块采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
所述室内空调负荷反弹量计算模块用于建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
所述优化模块用于将室内空调负荷反弹量计算模块得到的室内空调负荷反弹量通过PD迭代学习算法进行优化,获得优化后的负荷反弹量。
具体的,模型建立模块中先通过一阶等效热参数法建立单体空调负荷模型,然后通过能量守恒定律以及蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
具体的,模型建立模块建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
具体的,负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
实施例1
本实施例中公开了一种负荷反弹量优化方法,在上述实施例的基础上还公开了如下技术特征:
优选的,目标函数J1为:
J1=minLPB
=min∑15.202399x1 2-273.266156x1-2659.039304-132.8800x2 2+747.5x2+5.3630-0.000323x3 2+2.282920x3-9.865600
基于上述聚合模型,对夏季用户活跃时间(09:00-21:00)内的12个时段负荷分布情况分别进行模拟。假定有10组用户群参与DLC,每组用户群包含500个控制对象;用户参与需求响应期间内实施了若干次DLC,其控制时段互不连续。该时段内逐时室外温度变化情况如图2所示。
用户的热力学参数如表1所示。
表1用户热力学参数表
将总数为5000的用户分类成为10组,每组含500个用户,为了便于计算,将每组内的500个用户理想化为同类型同参数。其中,房间热容、房间热阻、空调额定功率采用MATLAB从连续高斯分布中随机获取生成。空调能效比统一取为2.50。一天内用户组1-10的自然负荷变化情况如图3所示。
假定用户组1-10的参数不变,此时假定所有用户选取室内温度上下限为[24℃,28℃]。以组别6为例,全天在自然状态下的负荷情况及受控于DLC后的负荷对比如图4所示。
由图4可以看出,在11:00-17:00采用了DLC控制方法之后,用户的负荷有了明显降低,削峰效果明显,且负荷削减量达到全天负荷的9.312%。但在上述条件下,在结束DLC控制时期11:00-17:00之后,在17:00-18:00之间用户组1-10会产生一定量的反弹,具体负荷反弹量及负荷反弹比例如表2所示。
表2用户组1-10负荷反弹量及负荷反弹比例
如表2所示,随着DLC控制时期的结束,会伴随着一定量的负荷反弹,且最高占整个时段负荷的9.312%。针对这种情况,确定了三个与直接负荷控制相关的参数作为负荷反弹指标,分别为直接负荷控制温度上限、直接负荷控制时长和参与负荷数量。
对上述三个指标在不同工况下进行仿真模拟,仿真数据采用正交最小二乘法进行拟合,结果分别如式所示。
y1=15.202399x1 2-273.266156x1-2659.039304
y2=-132.8800x2 2+747.5x2+5.3630
y3=-0.000323x3 2+2.282920x3-9.865600
y1、y2、y3分别表示对应工况下引起的负荷反弹总量。
如图5为用户组6运行在自然状态、传统DLC以及优化方法下的负荷变化情况。当***运行在传统DLC控制下,相比于自然状态,***的负荷有了明显的削减,但在集中控制期结束后,在17:00-18:00之间的负荷又达到一个高峰,且该负荷反弹几乎与未受控时的负荷持平,会对电网造成较大冲击。采用本发明所用优化方法,可以看出,虽然负荷值也有波动,但***的负荷削减效果明显,且负荷反弹量均低于受控前的负荷,整体运行在节能状态,相比于自然状态下节能效果明显。
采用NSGA2算法对提出的基于直接负荷控制的反弹抑制优化模型求解。具体参数设置为:遗传的迭代次数100,交叉率为2,变异率为2,种群数为5,其中每个种群规模为500;各种群之间的迁移率设为0.2。选取控制时段为9:00-21:00,一共12小时,直接负荷控制每个阶段为60min,共有12个阶段。需求响应负荷反弹抑制优化仿真的负荷情况如图5所示。
用户组1-10受需求响应负荷反弹抑制优化策略控制下,全天各时段内的受控情况如表3所示。
表3用户组在全天的受控情况
如表3所示,表示用户组1-20在全天各个时段受控情况,其中,1表示该时段该用户受控于该负荷反弹抑制优化策略,0表示不受控。在这样的受控情况下,即为本文所描述的负荷反弹抑制优化后的状态。相比于传统DLC控制方法,需求响应负荷反弹抑制策略使得用户的受控状态分散于全天时段内,对用户的舒适度有了一定的提高。
进一步,基于MATLAB,采用PD型迭代学习控制算法,设置仿真采样时间ts=1s,比例学习增益Tp=0.09,积分学习增益Td=0.995。用户组6在需求响应负荷反弹抑制优化策略控制下、自然状态下与普通DLC控制下,全天负荷的对比如表4。将经负荷反弹抑制策略优化后的用户组6的负荷数据作为期望轨迹进行仿真,结果如图6所示。
表4用户组6在三种状态下全天负荷的对比
经过1600次迭代优化,如图6,在使用PD型迭代学***稳,几乎没有波动和超调量,而且具有更短的调节时间。并且当负荷值发生变化时,***能够从相对稳定的一个状态平滑过渡到另一个状态,这样会使得负荷反弹量进一步降低,能够使电网运行的安全性进一步加大。
本发明基于一阶等效热参数法建立的空调负荷聚合模型,结合直接负荷控制控制方法,对其进行了不同工况下的模拟仿真,确定了三个负荷反弹指标,并对其进行了敏感度分析以及拟合。然后,制定了负荷反弹抑制优化策略以及基于迭代学***稳,进一步减少负荷反弹对电网运行的影响。发明的研究对电网负荷调度优化以及需求响应更好的展开具有一定现实意义。
Claims (10)
1.基于需求响应的负荷反弹量优化方法,其特征在于,包括如下步骤:
步骤1:建立室内空调负荷模型;
所述室内空调负荷模型包括房间温度变化模型、空调启停模型和空调***制冷量模型,如下所示:
1)房间温度变化模型
其中,Δt表示时间间隔;s(t)表示空调的运行状态,等于1的时候表示开启状态,等于0表示关闭状态;
2)空调启停模型
空调房间内温度控制上下限分别为:
式中,Tmax、Tmin分别表示室内允许的温度范围上下限;δ表示温度调节参数,取2℃;
当室温波动范围为[Tmin,Tmax],可求得启停周期中停机期和制冷期的时长分别为:
式中,τon、τoff分别表示启停周期内制冷期和停机期的时长;
3)空调***制冷量模型
根据能量守恒定律,一段时间内,空调***制冷量应该等于围护结构冷负荷、室内热源散热引起的冷负荷、人体散热引起的冷负荷及新风负荷之和;
Qc=Qe+Qh+Qp+Qn
其中:Qc表示空调***制冷量;Qe表示围护结构的冷负荷;Qh表示室内热源散热引起的冷负荷;Qp表示人体散热形成的冷负荷;Qn表示新风负荷;
步骤2:采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
步骤3:建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和步骤2得到的负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
所述整体降负荷目标函数J2表示为:
J2=min(max(LDLC,i,j))
其中,LDLC,i,j表示在i时段,j组用户的总体负荷,该目标函数的含义为:降低第i时段内第j组用户的总体负荷量的最大值;
步骤4:将步骤3得到的室内空调负荷反弹量通过PD控制算法进行优化,获得优化后的负荷反弹量。
2.如权利要求1所述的基于需求响应的负荷反弹量优化方法,其特征在于,步骤1中先通过一阶等效热参数法建立单体空调负荷模型,然后通过蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
3.如权利要求1所述的基于需求响应的负荷反弹量优化方法,其特征在于,步骤1建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
4.如权利要求1所述的基于需求响应的负荷反弹量优化方法,其特征在于,步骤2中负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
6.一种负荷调节方法,用于抑制在需求响应控制时期出现的负荷反弹现象,其特征在于,包括如下步骤:
步骤A:根据如权利要求1-5中任一种负荷反弹量优化方法,获得优化后的负荷反弹量;
步骤B:根据获得的优化后的负荷反弹量进行负荷控制和调节。
7.基于需求响应的负荷反弹量优化***,其特征在于,包括模型建立模块、负荷控制模块、室内空调负荷反弹量计算模块和优化模块;
所述模型建立模块用于建立室内空调负荷模型;
所述室内空调负荷模型包括房间温度变化模型、空调启停模型和空调***制冷量模型,如下所示:
1)房间温度变化模型
其中,Δt表示时间间隔;s(t)表示空调的运行状态,等于1的时候表示开启状态,等于0表示关闭状态;
2)空调启停模型
空调房间内温度控制上下限分别为:
式中,Tmax、Tmin分别表示室内允许的温度范围上下限;δ表示温度调节参数,取2℃;
当室温波动范围为[Tmin,Tmax],可求得启停周期中停机期和制冷期的时长分别为:
式中,τon、τoff分别表示启停周期内制冷期和停机期的时长;
3)空调***制冷量模型
根据能量守恒定律,一段时间内,空调***制冷量应该等于围护结构冷负荷、室内热源散热引起的冷负荷、人体散热引起的冷负荷及新风负荷之和;
Qc=Qe+Qh+Qp+Qn
其中:Qc表示空调***制冷量;Qe表示围护结构的冷负荷;Qh表示室内热源散热引起的冷负荷;Qp表示人体散热形成的冷负荷;Qn表示新风负荷;
所述负荷控制模块采用直接负荷控制方法对室内空调负荷模型进行模拟,得到室内空调负荷模型的运行参数,采用最小二乘法对室内空调负荷模型的运行参数进行参数拟合,获得拟合方程,将拟合方程作为负荷反弹最小的目标函数;
所述室内空调负荷反弹量计算模块用于建立整体降负荷目标函数,将室内温度上限、控制时长、参与负荷数作为约束条件,采用NSGA2算法对整体降负荷目标函数和负荷反弹最小的目标函数求解,获得室内空调负荷反弹量;
所述整体降负荷目标函数J2表示为:
J2=min(max(LDLC,i,j))
其中,LDLC,i,j表示在i时段,j组用户的总体负荷,该目标函数的含义为:降低第i时段内第j组用户的总体负荷量的最大值;
所述优化模块用于将室内空调负荷反弹量计算模块得到的室内空调负荷反弹量通过PD控制算法进行优化,获得优化后的负荷反弹量。
8.如权利要求7所述的基于需求响应的负荷反弹量优化***,其特征在于,模型建立模块中先通过一阶等效热参数法建立单体空调负荷模型,然后通过蒙特卡洛模拟法获得室内空调负荷模型,所述室内空调负荷模型包括聚合的多个空调。
9.如权利要求7所述的基于需求响应的负荷反弹量优化***,其特征在于,模型建立模块建立室内空调负荷模型后还包括如下步骤:
步骤a1:获取空调负荷历史数据,所述空调负荷历史数据包括采集日历史数据和预测日历史数据;
步骤a2:根据将室内空调负荷模型中的模型参数建立基于Elman神经网络的空调负荷预测模型,将采集日历史数据和预测日历史数据输入空调负荷预测模型进行训练,得到训练好的空调负荷预测模型;
步骤a3:获取训练好的空调负荷预测模型的模型参数,用训练好的空调负荷预测模型的模型参数更新室内空调负荷模型的模型参数,获得更新后的室内空调负荷模型。
10.如权利要求7所述的基于需求响应的负荷反弹量优化***,其特征在于,负荷反弹最小的目标函数J1为式Ⅰ:
J1=min∑a1x1 2-a2x1-C1-b1x2 2+b2x2+C2-c1x3 2+c2x3-C3 式Ⅰ
其中,x1表示直接负荷控制温度上限,x2表示直接负荷控制时长,x3表示直接负荷控制参数负荷数,a1、a2、b1、b2、c1、c2、C1、C2和C3均为常数,其中a1∈(15.0,15.5),a2∈(273.0,273.5),C1∈(-2659.5,-2659.0),b1∈(132.5,-133.0),b2∈(747.5,748.0),C2∈(5.0,5.5),c1∈(-0.5,0),c2∈(2.0,2.5),C3∈(-10,-9.5)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110016757.XA CN112665159B (zh) | 2021-01-07 | 2021-01-07 | 基于需求响应的负荷反弹量优化、负荷调节方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110016757.XA CN112665159B (zh) | 2021-01-07 | 2021-01-07 | 基于需求响应的负荷反弹量优化、负荷调节方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112665159A CN112665159A (zh) | 2021-04-16 |
CN112665159B true CN112665159B (zh) | 2021-12-21 |
Family
ID=75413416
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110016757.XA Active CN112665159B (zh) | 2021-01-07 | 2021-01-07 | 基于需求响应的负荷反弹量优化、负荷调节方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112665159B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104319777A (zh) * | 2014-10-15 | 2015-01-28 | 国家电网公司 | 一种需求侧负荷柔性控制方法 |
WO2015118749A1 (ja) * | 2014-02-10 | 2015-08-13 | 株式会社東芝 | 熱負荷推定装置および空調制御システム |
CN105071378A (zh) * | 2015-07-20 | 2015-11-18 | 东南大学 | 含柔性负荷的配电公司日前优化调度方法 |
CN105444343A (zh) * | 2015-09-30 | 2016-03-30 | 江苏省电力公司南京供电公司 | 一种基于用电舒适度的空调负荷优先级中断方法 |
CN105811454A (zh) * | 2016-03-14 | 2016-07-27 | 东南大学 | 一种考虑风电接入的直接负荷控制资源优化方法 |
CN106886951A (zh) * | 2017-02-23 | 2017-06-23 | 国网天津节能服务有限公司 | 一种基于最大供电能力的自动需求响应策略 |
CN107101322A (zh) * | 2017-04-13 | 2017-08-29 | 东南大学 | 统一最大削减负荷持续时间的变频空调群组潜力评估方法 |
CN108131722A (zh) * | 2017-12-15 | 2018-06-08 | 西南大学 | 一种面向电网调峰的终端用户制冷行为自适应调控 |
CN109409741A (zh) * | 2018-10-26 | 2019-03-01 | 国网上海市电力公司 | 一种中央空调的柔性调峰方法 |
CN109458706A (zh) * | 2018-11-07 | 2019-03-12 | 国网江苏省电力有限公司南通供电分公司 | 一种稳定的降负荷聚合商空调调控策略 |
CN109737570A (zh) * | 2019-01-21 | 2019-05-10 | 中国电力科学研究院有限公司 | 一种多联机运行调控方法及*** |
CN110425706A (zh) * | 2019-07-29 | 2019-11-08 | 南京理工大学 | 面向电网削峰的聚合空调负荷调控方法 |
CN110535142A (zh) * | 2019-08-22 | 2019-12-03 | 广州水沐青华科技有限公司 | 基于改进离散型pso算法的用电智能控制方法及计算机可读存储介质 |
CN110543969A (zh) * | 2019-07-31 | 2019-12-06 | 国网江苏省电力有限公司电力科学研究院 | 一种家庭用电行为优化算法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006019726A (ja) * | 2004-06-29 | 2006-01-19 | Samsung Electronics Co Ltd | ウェーハ移送システム及びシステム内の圧力調整方法 |
US10354345B2 (en) * | 2012-01-23 | 2019-07-16 | Whisker Labs, Inc. | Optimizing and controlling the energy consumption of a building |
CN104019520B (zh) * | 2014-05-20 | 2017-01-18 | 天津大学 | 基于spsa的制冷***最小能耗的数据驱动控制方法 |
CN104807143B (zh) * | 2015-05-14 | 2017-05-10 | 南通大学 | 基于电能友好空调负荷侧主动需求策略 |
CN105135623B (zh) * | 2015-09-17 | 2018-03-30 | 国网天津节能服务有限公司 | 一种满足电网调峰和用户舒适度的中央空调减载控制方法 |
US10564612B2 (en) * | 2016-06-30 | 2020-02-18 | Johnson Controls Technology Company | Variable refrigerant flow system with multi-level model predictive control |
US10655872B2 (en) * | 2017-04-06 | 2020-05-19 | Eaton Intelligent Power Limited | Load management system and method utilizing breaker metering data and occupancy data |
-
2021
- 2021-01-07 CN CN202110016757.XA patent/CN112665159B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015118749A1 (ja) * | 2014-02-10 | 2015-08-13 | 株式会社東芝 | 熱負荷推定装置および空調制御システム |
CN104319777A (zh) * | 2014-10-15 | 2015-01-28 | 国家电网公司 | 一种需求侧负荷柔性控制方法 |
CN105071378A (zh) * | 2015-07-20 | 2015-11-18 | 东南大学 | 含柔性负荷的配电公司日前优化调度方法 |
CN105444343A (zh) * | 2015-09-30 | 2016-03-30 | 江苏省电力公司南京供电公司 | 一种基于用电舒适度的空调负荷优先级中断方法 |
CN105811454A (zh) * | 2016-03-14 | 2016-07-27 | 东南大学 | 一种考虑风电接入的直接负荷控制资源优化方法 |
CN106886951A (zh) * | 2017-02-23 | 2017-06-23 | 国网天津节能服务有限公司 | 一种基于最大供电能力的自动需求响应策略 |
CN107101322A (zh) * | 2017-04-13 | 2017-08-29 | 东南大学 | 统一最大削减负荷持续时间的变频空调群组潜力评估方法 |
CN108131722A (zh) * | 2017-12-15 | 2018-06-08 | 西南大学 | 一种面向电网调峰的终端用户制冷行为自适应调控 |
CN109409741A (zh) * | 2018-10-26 | 2019-03-01 | 国网上海市电力公司 | 一种中央空调的柔性调峰方法 |
CN109458706A (zh) * | 2018-11-07 | 2019-03-12 | 国网江苏省电力有限公司南通供电分公司 | 一种稳定的降负荷聚合商空调调控策略 |
CN109737570A (zh) * | 2019-01-21 | 2019-05-10 | 中国电力科学研究院有限公司 | 一种多联机运行调控方法及*** |
CN110425706A (zh) * | 2019-07-29 | 2019-11-08 | 南京理工大学 | 面向电网削峰的聚合空调负荷调控方法 |
CN110543969A (zh) * | 2019-07-31 | 2019-12-06 | 国网江苏省电力有限公司电力科学研究院 | 一种家庭用电行为优化算法 |
CN110535142A (zh) * | 2019-08-22 | 2019-12-03 | 广州水沐青华科技有限公司 | 基于改进离散型pso算法的用电智能控制方法及计算机可读存储介质 |
Non-Patent Citations (1)
Title |
---|
基于ABC_SA混合算法的群控电梯优化调度;闫秀英;《计算机测量与控制》;20200825;第107-111页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112665159A (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lu et al. | Global optimization for overall HVAC systems––Part II problem solution and simulations | |
Yang et al. | Optimal control strategy for HVAC system in building energy management | |
CN105444343A (zh) | 一种基于用电舒适度的空调负荷优先级中断方法 | |
CN110848895B (zh) | 一种非工空调柔性负荷控制方法和*** | |
CN105790286A (zh) | 一种中央空调负荷聚合及参与平抑分布式电源出力的策略 | |
CN109737570B (zh) | 一种多联机运行调控方法及*** | |
CN115882463A (zh) | 一种商业建筑空调负荷可调度潜力评估方法 | |
CN110535142B (zh) | 基于改进离散型pso算法的用电智能控制方法及计算机可读存储介质 | |
Xing et al. | Chiller–pump system optimisation method for minimum energy operation based on multi-objective evolutionary algorithm | |
TW201027014A (en) | Method for managing air conditioning power consumption | |
WO2022198734A1 (zh) | 基于响应优先级的公共建筑电力需求响应双重优化方法 | |
CN109343461A (zh) | 一种楼宇自动化控制***的节能控制方法 | |
CN112665159B (zh) | 基于需求响应的负荷反弹量优化、负荷调节方法及*** | |
CN116085953A (zh) | 基于动态空调运行数据的节能控制方法、***及介质 | |
CN116294085A (zh) | 一种数据中心空调***运行方法 | |
CN115526504A (zh) | 泵站供水***节能调度方法、***、电子设备及存储介质 | |
CN109447368A (zh) | 一种对中央空调***进行基线负荷预测的方法 | |
CN113435042B (zh) | 一种建筑空调***需求响应的强化学习建模方法 | |
CN115103569A (zh) | 一种机房中设备控制的方法和装置 | |
CN113028604A (zh) | 基于聚合空调的温度调节控制及调度方法 | |
CN110912163B (zh) | 一种温控负荷能效电厂的频率控制策略的获取方法 | |
CN113420413A (zh) | 一种基于负荷可塑性的柔性负荷可调节能力量化方法及*** | |
CN113158397A (zh) | 一种暖通负荷虚拟电厂频率调整策略 | |
CN111623491A (zh) | 一种基于协同优化策略的变速水泵运行调节方法 | |
Zhou et al. | Research on Optimization Strategy of Central Air-conditioning System for Power Demand Response |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |