CN109787251B - 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 - Google Patents
一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 Download PDFInfo
- Publication number
- CN109787251B CN109787251B CN201910038801.XA CN201910038801A CN109787251B CN 109787251 B CN109787251 B CN 109787251B CN 201910038801 A CN201910038801 A CN 201910038801A CN 109787251 B CN109787251 B CN 109787251B
- Authority
- CN
- China
- Prior art keywords
- temperature
- cluster
- temperature control
- power
- loads
- 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
- Feedback Control In General (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种集群温控负荷的双线性离散模型。由于实际系统参数难以获取,本发明同时设计了基于模型分解的带遗忘因子递推最小二乘法的系统参数辨识方法,通过对系统参数进行辨识的方法获取系统建模参数。并且基于集群温控负荷双线性离散模型设计了基于反推控制的集群温控负荷聚合功率控制方法。从理论和仿真上证明了系统建模的准确性,以及系统参数辨识方法和控制方法的稳定性和有效性。相比其它模型及控制方法,本发明的新颖性在于建立了时间域和温度域均离散化的双线性离散模型,并将系统参数辨识方法作为系统参数获取的方式,并设计了性能良好的集群温控负荷功率的控制器,实现集群温控负荷对目标功率曲线的良好跟踪。
Description
技术领域
本发明涉及集群温控负荷的建模、系统参数辨识与控制方法,属于电力系统辅助服务和需求侧响应控制领域。
背景技术
随着新能源领域技术的发展,大量的以风电和光伏为代表的间歇式电源并入电力系统中,研究显示,间歇式电源输出功率呈现出分钟级或小时级的波动,波动幅度甚至可从最大输出骤降至零,这将对电网运行造成严重的冲击和影响。与此同时,温控负荷(Thermostatically Controlled Loads,TCL)占电网总负荷的比重日益增加。在我国大中城市,夏季空调负荷占总负荷的比例普遍高达30%-40%,更有个别城市超过50%。因TCL所具有的储能特性和灵活调度特性,需求侧“海量”的TCL,为平抑间歇式电源功率波动带来了丰富和优质的可调控资源。
当前,TCLs辅助电力系统调控的研究热点聚集在集群温控负荷(AggregatedThermostatically Controlled Loads,ATCLs)聚合功率的建模和控制方法。对于,ATCLs模型的建立。有文献提出了基于历史数据的回归模型、基于耦合Fokker-Plank方程的统计扩散模型、基于耦合Fokker-Plank方程的概率特性模型等,这些模型大都结构复杂、非常计算量大,难以直接用于实际控制程序;在控制方法上,有基于模糊理论的状态序列法、基于智能机器学习理论的状态序列法、基于温度延伸裕度的时变设定温度的控制方法等,这些控制方法大多针对着不同的场景,控制算法需较多的信息数据,控制的计算量大。因此,构建一种精确的ATCLs聚合功率模型,设计其计算量小的控制方法,对ATCLs参与电力系统运行有着重要的意义。
发明内容
针对上述问题,本发明的目的在于针对现有技术的不足,做出三点创新:一是提出了一种ATCLs双线性离散模型,推导出了具体的ATCLs聚合功率解析表达式以便于控制实施。二是引入了一种基于模型分解的带遗忘因子递推最小二乘法(Model Decompositionbased Recursive Least Squares algorithm with forgetting factor,MD-FRLS)的ATCLs模型参数辨识方法,以便于精确建模。三是基于反推控制原理设计了ATCLs聚合功率的控制方法,实现ATCLs功率对目标功率曲线的跟踪。
本发明所述的问题是由以下技术方案实现的:
首先,假定每个TCL具有均质性,建立ATCLs双线性离散数学模型:
式中:xj(k)、xj(k+1)分别为第k时刻、第k+1时刻第j段温度区间负荷的数量,xj-1(k)、xj+1(k)分别为第k时刻第j-1段、第j段温度区间负荷的数量,αon/off代表“开/关”状态下温度的平均变化率,单位“℃/h”,Q为划分的总离散温度区间,M为“ON”或“OFF”状态下的温度区间数量,Tset为温度设定值,为k时刻温度设定值的变化率,ΔT为离散化温度长度,k为采样点,Δτ为时间离散化步长。
边界小段的负荷数量变化表达式为:
式中:xQ(k)、xQ(k+1)分别为第k时刻、第k+1时刻第Q段温度区间负荷的数量,xM(k)、xM+1(k)为第k时刻第M段、第M+1段温度区间负荷的数量,x1(k)、x1(k+1)为第k时刻、第k+1时刻第1段温度区间负荷的数量。
ATCLs的聚合功率表达式为:
式中:P为温控负荷的等效功率,PTCL(k+1)为k+1时刻ATCLs的聚合功率,η为温控负荷的能量转换效率。
将温度区间负荷的数量代入ATCLs的聚合功率表达式中。最终,ATCLs聚合功率的离散表达式为:
式中:xM+1(k+1)为第k+1时刻第M+1段温度区间负荷的数量,Non(k+1)为第k+1时刻“开”状态下系统中TCLs个数。
其次,准确获取系统参数是ATCLs系统精确建模的基础,对于基于MD-FRLS的ATCLs模型参数辨识方法,所述辨识方法为:以MD-FRLS为基础,对ATCLs聚合功率的离散模型进行处理得到:
在前次估计结果的基础上,就新的数据根据递推规则对前次估计的结果进行修正,得出新的参数估计值,直至参数估计值收敛。
最后,根据反推控制原理,设计ATCLs聚合功率的控制方法。给定跟踪目标功率Pt,推导出ATCLs的控制量u(k),即ATCLs的功率控制器:
式中:l为控制增益,eP=e(k+1)-e(k),e(k+1)、e(k)分别为第k+1时刻、第k时刻的ATCLs功率跟踪误差。
附图说明
图1为常参数下系统参数αon的辨识结果;
图2为常参数下系统参数αoff的辨识结果;
图3为变参数下系统参数αon的辨识结果;
图4为变参数下系统参数αoff的辨识结果;
图5为目标功率跟踪结果;
图6为目标功率跟踪局部放大图;
图7为[TLow,THigh]区间内的离散化ATCLs温度域动态过程。
文中各符号为:xj(k)代表第j段温度区间负荷的数量,αon/off代表“开/关”状态下温度的平均变化率,单位“℃/h”,Q为划分的总离散温度区间,M为“ON”或“OFF”状态下的温度区间数量,Δ温度死区,Tset为温度设定值,THigh为温度设定值死区温度的上限,TLow为温度设定值死区温度的下限,为温度设定值的变化率,ΔT为离散化温度长度,k为采样点,Δτ为时间离散化步长,P代表温控负荷的等效功率,Non(k+1)代表k+1时刻“开”状态下系统中TCLs个数,Noff(k+1)分别代表k+1时刻“/关”状态下系统中TCLs个数,PTCL(k+1)代表k+1时刻ATCLs的聚合功率,η代表温控负荷的能量转换效率,A1和A2为分解后的待辨识参数,L1(k)和L2(k)为增益向量,P1(k)和P2(k)为协方差矩阵,和为信息向量,λ为遗忘因子,y(k)为系统的输出向量,k为采样点,Pt为踪目标功率,u(k)为ATCLs的控制量,l为控制增益。
具体实施方式
本发明由以下技术方案实现:
1.ATCLs聚合双线性时间离散模型的建立
假定假定每个TCL具有均质性,每个TCL都能在激励信号下正常响应。ATCLs双线性离散数学模型为:
其中:xj(k)、xj(k+1)分别为第k时刻、第k+1时刻第j段温度区间负荷的数量,xj-1(k)、xj+1(k)分别为第k时刻第j-1段、第j段温度区间负荷的数量,αon/off代表“开/关”状态下温度的平均变化率,单位“℃/h”,Q为划分的总离散温度区间,M为“on”或“off”状态下的温度区间数量,ΔT为离散化温度长度,k为采样点,Δτ为时间离散化步长,Tset为温度区间[TLow,THigh]中的温度设定值,为k时刻为温度设定值的变化率。Tset表达式为:
式中,Δ为温度上限与温度下限的差,即温度死区。
边界小段的负荷数量变化表达式为:
式中:xQ(k)、xQ(k+1)分别为第k时刻、第k+1时刻第Q段温度区间负荷的数量,xM(k)、xM+1(k)为第k时刻第M段、第M+1段温度区间负荷的数量,x1(k)、x1(k+1)为第k时刻、第k+1时刻第1段温度区间负荷的数量。
将式(1)的第一式与式(2)的第一式求和,可得:
将式(1)的第二式与式(2)的第二式求和,可得:
若Non(k+1)、Noff(k+1)分别代表k+1时刻“开/关”状态下系统中TCLs个数,将其代入式(3)、(4)中,得到:
Non(k+1)与Noff(k+1)相加即为k+1时刻集群系统中TCLs的总数。
假设本研究TCL具有均质性,每个TCL的等效功率和能量转换效率都为P和η,那么,在式(14)两边同乘以P/η,就可得到一个数量为Q的ATCLs系统的整体聚合功率。令PTCL(k+1)代表k+1时刻ATCLs的聚合功率,则ATCLs聚合功率的离散表达式可写为:
2.基于MD-FRLS的ATCLs系统参数辨识
准确的系统参数是ATCLs系统精确建模的基础,基于MD-FRLS的系统参数辨识是工业中常用的一种辨识方法:
基于以上方法,提出一种基于MD-FRLS的ATCLs系统参数辨识算法。假设ATCLs总功率PTCL已知,根据MD-FRLS算法架构将式(7)进行相应改写,可得:
其中:
A1=αon (19)
A2=αoff (20)
将式(16)至(20)代入式(8)至(13)中,即可辨识得到ATCLs的系统参数
3.基于反推控制原理的ATCLs聚合功率控制方法
本发明的控制目标是控制ATCLs对目标功率进行跟踪。ATCLs控制问题的本质是功率跟踪问题,通过控制算法实现ATCLs对参考功率的准确跟踪。
基于反推控制原理,利用改进的ATCLs聚合双线性离散模型及其参数辨识结果来设计控制器。假设ATCLs参数辨识结果为其中和分别代表αon和αoff的辨识值。同时,设Pt为跟踪信号,表征了间歇式电源的输出功率波动。令e=Pt(k)-PTCL(k),e为功率跟踪误差。根据反推控制原理,可得:
eP=e(k+1)-e(k)=Pt(k+1)-Pt(k)-(PTCL(k+1)-PTCL(k)) (21)
结合式(21)和(22),可得:
取Lyapunov函数V为:
根据李雅普诺夫稳定性第二定理,若控制系统稳定需取ΔV<0,即:
令:
可得到控制器为:
式中:l为控制增益。
此时:
即,当-1<l<1时ΔV<0,由此从理论上证明了该系统是稳定的。
最后,通过算法仿真验证分析,检验控制器的性能。在系统参数辨识仿真中设计了变参数识别和常参数识别仿真,并加入递推最小二乘辨识方法(Recursive Least Squaresalgorithm,RLS)与之对比。在功率控制仿真中给定了目标跟踪功率曲线,验证控制方法的有效性和准确性。
选取由10000个TCL构成的ATCLs系统,在Matlab/Simulink中搭建ATCLs参数辨识和控制系统并进行仿真。系统参数辨识仿真各项参数:Q=20,M=10,Δ=0.5℃,Δτ=10- 4s,ΔT=0.05℃,λ=0.9,R=2℃/kW,C=10(kW·h)/℃,P=14kW,η=2.5,Tset=20℃,初始状态下,每个温度区间开/关状态下的温控负荷初始值均设为500个。控制系统仿真各项参数:Q=20,M=10,Δ=0.5℃,Δτ=10-4s,ΔT=0.05℃,R=2℃/kW,C=10(kW·h)/℃,P=14kW,η=2.5, l=0.9,初始状态下每个温度区间开/关状态下的温控负荷初始值均设为500个。
结果如图1至图6所示。
图1至图4给出在系统参数为常参数下,MD-FRLS和RLS的辨识结果。其中实线(MD-FRLS)和点画线(RLS)分别代MD-FRLS算法和RLS算法的辨识结果,虚线(STD)代表系统实际参数。从图1至图2中可见,MD-FRLS和RLS都能准确的辨识出系统参数,但MD-FRLS算法的辨识速度快于RLS算法,辨识精度也更高。图3至图4给出了系统参数为变参数下,MD-FRLS和RLS的辨识结果。在系统参数发生变化时MD-FRLS算法仍能快速、准确地辨识出系统参数,而RLS算法则需要经过较长时间才能振荡收敛,并且其辨识结果还存在一定的误差。由此表明,MD-FRLS算法应用于变参数系统时亦具有较高的准确性,且与RLS算法相比,MD-FRLS算法所需时间更短、辨识结果更精确。图5至图6给出了将系统参数辨识结果用于建模后的ATCLs系统对目标功率跟踪仿真结果。图中点画线(PTCL)代表实时ATCLs负荷功率,实线(Pt)代表功率跟踪目标曲线,虚线(Dev)代表ATCLs对参考功率曲线的跟踪误差。图5表明,设计的控制方法能够使ATCLs聚合功率很好地跟踪时变的目标功率。图6为给功率跟踪结果的局部放大图。可见算法的跟踪效果良好。因此,提出的ATCLs聚合功率建模方法、系统参数辨识方法和ATCLs聚合功率控制方法具有良好实际意义。
Claims (3)
1.一种集群温控负荷聚合模型,其特征在于,对于集群温控负荷聚合的建模,所述建模为:
首先,假定每个温控负荷具有均质性,建立集群温控负荷双线性离散数学模型:
式中:xj(k)、xj(k+1)分别为第k时刻、第k+1时刻第j段温度区间负荷的数量,xj-1(k)、xj+1(k)分别为第k时刻第j-1段、第j段温度区间负荷的数量,αon/off为“开/关”状态下温度的平均变化率,单位“℃/h”,Q为划分的总离散温度区间,M为“ON”或“OFF”状态下的温度区间数量;Tset为温度设定值,为k时刻温度设定值的变化率,ΔT为离散化温度长度,k为采样点,Δτ为时间离散化步长;
边界小段的负荷数量变化表达式为:
式中:xQ(k)、xQ(k+1)分别为第k时刻、第k+1时刻第Q段温度区间负荷的数量,xM(k)、xM+1(k)为第k时刻第M段、第M+1段温度区间负荷的数量,x1(k)、x1(k+1)为第k时刻、第k+1时刻第1段温度区间负荷的数量;
集群温控负荷的聚合功率表达式为:
式中:P为温控负荷的等效功率,PTCL(k+1)为k+1时刻集群温控负荷的聚合功率,η为温控负荷的能量转换效率;
将温度区间负荷的数量带入集群温控负荷的聚合功率表达式中;最终,集群温控负荷聚合功率的离散表达式为:
式中:xM+1(k+1)为第k+1时刻第M+1段温度区间负荷的数量,Non(k+1)为第k+1时刻“开”状态下系统中温控负荷个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038801.XA CN109787251B (zh) | 2019-01-16 | 2019-01-16 | 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038801.XA CN109787251B (zh) | 2019-01-16 | 2019-01-16 | 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109787251A CN109787251A (zh) | 2019-05-21 |
CN109787251B true CN109787251B (zh) | 2022-11-18 |
Family
ID=66500698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910038801.XA Active CN109787251B (zh) | 2019-01-16 | 2019-01-16 | 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109787251B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112838580B (zh) * | 2019-11-22 | 2023-03-24 | 华北电力大学(保定) | 一种改进异质温控负荷聚合模型及多目标协调控制方法 |
CN113078629B (zh) * | 2020-01-06 | 2022-08-26 | 华北电力大学(保定) | 一种用于集群温控负荷聚合商功率调控的聚合体功率分配模型及分布式一致性控制方法 |
CN114069599B (zh) * | 2020-07-29 | 2023-12-19 | 华北电力大学(保定) | 预设性能条件下聚合温控负荷改进模型预测控制 |
CN114079293B (zh) * | 2020-08-18 | 2023-08-15 | 华北电力大学(保定) | 一种计及温控负荷与微型燃气轮机的电-气互联系统新能源功率波动平滑方法 |
CN112632800B (zh) * | 2021-01-07 | 2023-11-28 | 广东华中科技大学工业技术研究院 | 一种3d盖板玻璃模具高精度温度控制方法、系统及设备 |
CN113067340B (zh) * | 2021-03-25 | 2023-02-28 | 山东大学 | 一种面向恒温控制负荷系统的动态状态估计方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427011A (zh) * | 2015-10-27 | 2016-03-23 | 中国电力科学研究院 | 一种温控负荷聚合系统稳定性分析方法 |
CN107276069A (zh) * | 2017-06-08 | 2017-10-20 | 中国电力科学研究院 | 一种地区电网居民温控负荷的近似聚合建模方法及系统 |
-
2019
- 2019-01-16 CN CN201910038801.XA patent/CN109787251B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427011A (zh) * | 2015-10-27 | 2016-03-23 | 中国电力科学研究院 | 一种温控负荷聚合系统稳定性分析方法 |
CN107276069A (zh) * | 2017-06-08 | 2017-10-20 | 中国电力科学研究院 | 一种地区电网居民温控负荷的近似聚合建模方法及系统 |
Non-Patent Citations (1)
Title |
---|
Modeling and Control of Aggregate Air Conditioning Loads for Robust Renewable Power Management;Saeid Bashash等;《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》;20130731;第21卷(第4期);1318-1327 * |
Also Published As
Publication number | Publication date |
---|---|
CN109787251A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109787251B (zh) | 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法 | |
CN107203137B (zh) | 非线性抽水蓄能机组调速系统增益启发式自适应pid控制方法 | |
CN111563610A (zh) | 一种基于lstm神经网络的建筑物电负荷综合预测方法及系统 | |
CN103106544B (zh) | 一种基于t-s型模糊神经网络的光伏发电预测系统 | |
CN103773900B (zh) | 基于神经网络及粒子群算法的固态发酵控制方法 | |
CN104933483A (zh) | 一种基于天气过程划分的风电功率预测方法 | |
CN105354620A (zh) | 一种风机发电功率的预测方法 | |
CN110807554A (zh) | 基于风电/光伏经典场景集的生成方法及系统 | |
CN107276069A (zh) | 一种地区电网居民温控负荷的近似聚合建模方法及系统 | |
CN105262145B (zh) | 一种新能源混合系统控制参数的优选方法 | |
CN104036328A (zh) | 自适应风电功率预测系统及预测方法 | |
CN105787594A (zh) | 基于多元时间序列及回归分析的辐照预测方法 | |
CN112182709A (zh) | 大型水库叠梁门分层取水设施的下泄水温快速预测方法 | |
CN108038568A (zh) | 一种基于粒子群算法的变权组合电力负荷短期预测方法 | |
CN103823504A (zh) | 一种基于最小二乘支持向量机的最大功率跟踪控制方法 | |
CN105888971A (zh) | 一种大型风力机叶片主动降载控制系统及方法 | |
CN102509152A (zh) | 一种基于rbf神经网络的开关磁阻电机在线建模方法 | |
CN102645894A (zh) | 模糊自适应动态规划方法 | |
CN109408896B (zh) | 一种污水厌氧处理产气量多元智能实时监控方法 | |
CN116865343B (zh) | 分布式光伏配电网的无模型自适应控制方法、装置及介质 | |
CN103279030A (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
CN110555566B (zh) | 一种基于b样条分位数回归的光电概率密度预测方法 | |
CN104182854A (zh) | 一种面向企业能源管理系统的混合能耗测算的方法 | |
CN115796327A (zh) | 一种基于vmd和iwoa-f-gru模型的风电功率区间预测方法 | |
Vergini et al. | Advanced state fuzzy cognitive maps applied on nearly zero energy building model |
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 |