CN107842908B - 一种基于环境参数补偿的实时供热负荷控制方法 - Google Patents

一种基于环境参数补偿的实时供热负荷控制方法 Download PDF

Info

Publication number
CN107842908B
CN107842908B CN201710886793.5A CN201710886793A CN107842908B CN 107842908 B CN107842908 B CN 107842908B CN 201710886793 A CN201710886793 A CN 201710886793A CN 107842908 B CN107842908 B CN 107842908B
Authority
CN
China
Prior art keywords
control
model
load
heat supply
moment
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.)
Expired - Fee Related
Application number
CN201710886793.5A
Other languages
English (en)
Other versions
CN107842908A (zh
Inventor
刘西陲
沈炯
张俊礼
王硕
潘天尧
申鹂
王舒颦
颜虹
张如通
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Southeast University
Nari Technology Co Ltd
Hangzhou Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Southeast University
Nari Technology Co Ltd
Hangzhou Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Southeast University, Nari Technology Co Ltd, Hangzhou Power Supply Co of State Grid Zhejiang Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710886793.5A priority Critical patent/CN107842908B/zh
Publication of CN107842908A publication Critical patent/CN107842908A/zh
Application granted granted Critical
Publication of CN107842908B publication Critical patent/CN107842908B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24DDOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
    • F24D19/00Details
    • F24D19/10Arrangement or mounting of control or safety devices
    • F24D19/1006Arrangement or mounting of control or safety devices for water heating systems
    • F24D19/1009Arrangement or mounting of control or safety devices for water heating systems for central heating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Thermal Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于环境参数补偿的实时供热负荷控制方法,包括:供热负荷预报步骤:采用了基于时间序列的供热负荷预报方法,对供热系统负荷进行预报;环境参数识别补偿步骤:将环境参数作为扰动,设置前馈补偿回路消除环境参数变化对供热温度的影响;预测控制步骤:采用预测控制,根据预测的负荷需求,进行操作,提前对大惯性供热负荷系统进行控制调节操作。与现有技术相比,本方法采用预测控制对供热系统进行负荷控制,并设计前馈补偿回路消除环境参数波动对供热的影响。控制效果好,具有一定的抗扰动能力,有效的提高了供热品质。

Description

一种基于环境参数补偿的实时供热负荷控制方法
技术领域
本发明涉及一种基于环境参数补偿的实时供热负荷控制方法,属于热工控制领域。
背景技术
供热行业密切联系人们的日常生活,对热电厂经济性的发展有着一定的作用。
目前我国供热负荷自动化控制技术比较薄弱,主要依靠工作人员的经验进行调节,因而造成了集中供热效率低、供热品质差、系统热能效益低、用户不能自行设定和调整用热温度。这些问题影响了供热质量、增加了热网的运行费用、浪费了能源、降低了供热的经济性。
由于供热管网较长,供热温度被控对象存在较大的惯性和延迟,造成传热过程的滞后。供热系统存在阀门等设备的非线性特征和一、二次供热回路间的耦合特性,导致传统控制方法难以取得满意的控制效果。同时,温度、风速等环境参数的波动对供热系统供热品质有一定干扰作用,很难及时消除环境参数的影响,保证供热品质。
目前大多数的供热管网控制系统仍采用常规PID(比例-积分-微分)串级控制方案,难以有效应对被控对象的大延迟、非线性特征。本发明有效利用预测控制(动态矩阵控制)的思想,进行滚动优化求解每一采样时刻的最优控制量,同时设计前馈补偿回路减小环境参数波动对供热温度的影响,保证控制品质。
发明内容
发明目的:针对上述存在的问题和不足,本发明提供一种热网控制的有效方法,能够解决上述问题,提高系统抗扰动能力,保证供热系统供热品质。
技术方案:为实现上述发明目的,本发明采用的技术方案为:
一种基于环境参数补偿的实时供热负荷控制方法,其特征在于,包括:
供热负荷预报步骤:采用了基于时间序列的供热负荷预报方法,对供热系统负荷进行预报;
环境参数识别补偿步骤:将环境参数作为扰动,设置前馈补偿回路消除环境参数变化对供热温度的影响;
预测控制步骤:采用预测控制,根据预测的负荷需求,进行操作,提前对大惯性供热负荷系统进行控制调节操作。
所述环境参数识别补偿步骤具体为:对可测环境参数的动态特性进行建模,提取对过程动态影响的因果性信息,利用广义最小二乘方法估计出可测环境参数与系统输出之间的辨识模型,求取辨识模型的阶跃响应特征,得到前馈补偿矩阵,减小环境参数波动对供热品质的影响。
所述供热负荷预报步骤中,以时间为预测对象的唯一影响因素,采用惯性原理对预测对象未来的状态进行推测,最后把预测供热负荷作为供热负荷控制系统的参考值。
所述供热负荷预报步骤包含:
11)针对k时刻负荷序列Q(k),通过求取负荷差分ΔQ(k)获得平稳随机序列,建立相应的ARMA预报模型,预报模型描述为:
A(z-1)ΔQ(k)=C(z-1)ξ(k) (9)
式中:A(z-1)=1+a1z-1+…anz-n,C(z-1)=1+c1z-1+cnz-n;A(z-1)和C(z-1)是稳定的多项式,a1,a2,...,an和c1,c2,...,cn是多项式系数,Δ=1-z-1为差分算子,ξ(k)为均值为零、方差为σ2的白噪声序列,n为多项式阶次,k为当前时刻的标度;
12)式(1)中的参数采用减消记忆的递推最小二乘方法进行估计,求取多项式系数的估计值
Figure BDA0001419166660000021
Figure BDA0001419166660000022
后,由最小方差预报原理,获取预报值:
Figure BDA0001419166660000023
其中,为在k时刻对k+i时刻负荷差分的预报值;
Figure BDA0001419166660000025
Figure BDA0001419166660000026
为系数为估计值的多项式,下标1代表迭代次数为1;i为相对于当前k时刻的未来时点序号;
13)获取供热负荷预报值Q(k+j|k)为:
Figure BDA0001419166660000027
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值,P为预测步长。
所述预测控制步骤包含:
21)获取热网对象的控制模型,稳态时对对象进行阶跃响应试验;
22)建立控制系统的预测模型yPM(k)=yP0(k)+ADΔuM(k),通过求解性能指标
Figure BDA0001419166660000028
求出最优控制量Δu(k)进行滚动优化;利用反馈项e(k+1)=y(k+1)-y1(k+1|k)进行误差校正,保证控制精度;式中:yP0(k)为k时刻起控制作用保持不变时,未来P个时刻内对象的输出;AD为动态矩阵;J(k)为k时刻的代价函数;wP(k)为k时刻时未来P个时刻内的设定值序列;Q,R分别为输出误差权矩阵和控制权矩阵;yPM(k)为k时刻起对象受到M个连续控制增量ΔuM(k)作用后在未来P个时刻内的输出;M为控制步长;y(k+1)为对象的在k+1时刻的实际输出;y1(k+1|k)为对象在k时刻受到1个控制增量,经过一步采样时间后在k+1时刻的预测输出;e(k+1)为模型误差;
23)建立模型误差差分ω(k)的ARMA模型为:
ω(k)+a1ω(k-1)+…anω(k-n)=ξ(k)+c1ξ(k-1)+…+cnξ(k-n) (12)
其中,ω(k),...,ω(k-n)为误差差分序列,多项式参数a1,a2,…an和c1,c2,…cn采用ARMA的参数估计法进行估计;再利用最小方差预报原理,求取k时刻对未来的误差差分预报值
Figure BDA0001419166660000031
为:
Figure BDA0001419166660000032
式中:i、j代表相对k时刻的未来时点序号;
最后求取在k时刻对k+i时刻的预报误差
Figure BDA0001419166660000033
式中e(k)为k时刻的模型误差,用预报误差代替e(k)提高模型精度。
步骤21)中,对一级网流量阀和二级网循环水泵变频器频率加入10%的阶跃扰动信号,采样时间为10s,获得二级网供回水流量和温差的阶跃响应特征。
所述环境参数识别补偿步骤包含:
31)利用最终预报误差准则估计环境参数对供热温度扰动模型的阶次,求解下式获取扰动模型阶次
Figure BDA0001419166660000034
其中,
Figure BDA0001419166660000036
为模型输出残差的方差,L为数据长度。
32)利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:
Figure BDA0001419166660000041
式中:A(z-1),B(z-1)和C(z-1)是稳定的多项式;
33)假设噪声模型C(z-1)参数未知,用迭代估计的方法求得;
a)令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次数;
b)计算输出估计值
Figure BDA0001419166660000042
输入估计值再更新迭代次数i=i+1,式中y(k),u(k)分别为k时刻实际输入和实际输出;
c)最小二乘法估计出
Figure BDA0001419166660000044
中的参数,式中:Ai(z-1),Bi(z-1)是迭代i次的稳定的多项式;
d)用估计模型
Figure BDA0001419166660000045
以及各时刻的观测数据,
Figure BDA0001419166660000046
分别为迭代i次、参数为估计值且与Ai(z-1),Bi(z-1)同形的多项式,计算出估计残差
Figure BDA0001419166660000047
e)计算估计残差平方和
Figure BDA0001419166660000048
及估计残差平方和增量ξ=ξii-1,如果ξ小于增量限值,则结束辨识,否则继续下一步;
f)对于新得到的噪声模型
Figure BDA0001419166660000049
用最小二乘法估计出参数,得到更新i次后的噪声模型Ci(z-1),返回第b)步,继续按照步骤b)至f)顺序进行计算;
34)已知供热系统环境参数扰动模型,对该模型进行阶跃响应试验,每隔10s采集一次数据,共采集N组数据组成前馈补偿矩阵
Figure BDA00014191666600000410
N为模型长度;
35)依据前馈补偿修改动态矩阵控制算法预测模型为:
yPM(k)=yP0(k)+ADΔuM(k)+bΔvP(k) (15)
式中:ΔvP(k)为k时刻时未来P个时刻内的扰动增量;
相应地求解优化问题minJ(k)得到控制律为:
Δu(k)=dT[wP(k)-yP0(k)-bΔvP(k)] (16)
其中,优化矩阵
Figure BDA00014191666600000411
Q,R分别为输出误差权矩阵和控制权矩阵。
本发明基于环境参数补偿的实时供热负荷控制方法,主要采取以下措施:1)采用基于时间序列法的供热负荷短期预报方法,预测未来的供热负荷并作为控制系统的参考值;2)将环境参数看做扰动,认为该扰动对供热流量没有影响。采用带前馈的动态矩阵控制方法进行控制,设计前馈补偿回路,消除环境参数波动对供热温度的影响;3)考虑到环境参数的不可控因素,缺乏环境参数对输出动态影响的阶跃响应模型,因而采用输入-输出数据进行辨识,以获得较为准确的阶跃响应模型,方便前馈回路设计。
与现有技术相比,本方法采用预测控制(动态矩阵控制)对供热系统进行负荷控制,并设计前馈补偿回路消除环境参数波动对供热的影响。控制效果好,具有一定的抗扰动能力,有效的提高了供热品质。
有益效果:当外界温度等环境参数发生变化时,保持供热温度不变,减小环境参数扰动对供热品质的影响。此时,增加或减少供热出水流量,用户端水流量不变,响应减小或增加混水流量。预测供热负荷会根据负荷时间序列在线进行估计,也包含环境参数变化的影响,本质上也能提高抵抗环境参数扰动作用。
附图说明
图1为本发明的系统结构框图。
具体实施方式
下面结合附图和具体实施方案,进一步阐明本发明,应该理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定范围。
本发明的基于环境参数补偿的实时供热负荷控制方法,通过预测控制的方法,设置前馈补偿回路实现对供热温度和供热流量的控制。同时有效处理对象大迟延特性,提高供热系统对环境参数的抗扰动能力。
本发明实时供热负荷控制方法具体步骤如下:
1)采用了基于时间序列的供热负荷预报方法。测量中对象与外部因素之间有着密切的联系,时间序列中的数据也是外部因素综合作用的反映。为了简便研究过程,对外部因素复杂的作用过程进行简化处理,即以时间为预测对象的唯一影响因素,采用惯性原理对预测对象未来的状态进行推测。最后把预测供热负荷作为供热负荷控制系统的参考值;
已知k时刻及历史负荷序列Q(k),Q(k-1),...,通过求取负荷差分ΔQ(k)获得平稳随机序列,建立相应的ARMA预报模型。预报模型描述如式(17)所示:
A(z-1)ΔQ(k)=C(z-1)ξ(k) (17)
式中:A(z-1)=1+a1z-1+…anz-n,C(z-1)=1+c1z-1+cnz-n。A(z-1)和C(z-1)是稳定的多项式,a1,a2,...,an和c1,c2,...,cn是多项式系数,Δ=1-z-1为差分算子,ξ(k)为均值为零、方差为σ2的白噪声序列,n为多项式阶次,k为当前时刻的标度;
式(17)中的多项式系数采用最小二乘方法进行估计,记多项式系数向量θ=[a1,...,an,c1...,cn]T,k时刻负荷预报值序列-Q(k),...,-Q(k-n+1)]T,k时刻至N+k-1时刻预报序列集合的转置
Figure BDA0001419166660000062
N为模型长度。
估计值
Figure BDA0001419166660000063
Figure BDA0001419166660000064
可通过求解如式(18)的多项式系数向量表达式得到:
θ=(ΦTΦ)-1ΦTQN (18)
求取估计值后,由最小方差预报原理,获取预报值:
其中,
Figure BDA0001419166660000066
为在k时刻对k+i时刻负荷差分的预报值;为系数为估计值的多项式,下标1代表迭代次数为1;i为相对于当前k时刻的未来时点序号。
获取供热负荷预报值Q(k+i|k)为:
Figure BDA0001419166660000069
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值;P为预测步长;
2)针对供热负荷系统大惯性的特点,采用预测控制进行控制调节。根据预测的负荷需求,提前进行操作,保证用户的用热品质。同时保证控制系统优良的跟踪性能和鲁棒性,增强了对干扰的抑制能力;
获取热网对象的控制模型,设热网供热负荷的传递函数可以由下式表示:
Figure BDA00014191666600000610
其中y1和y2分别为二次网供回水温差和二级网流量,u1和u2分别为调节一级网供水阀门和调节二级网循环水泵变频器频率。
Figure BDA0001419166660000071
为供热过程耦合模型传递函数矩阵;g11(s)为调节一级网供水阀门对二次网供回水温差的传递函数;g12(s)为调节二级网循环水泵变频器频率对二次网供回水温差的传递函数;g21(s)为调节一级网供水阀门对二级网流量的传递函数;g22(s)为调节二级网循环水泵变频器频率对二级网流量的传递函数。各传递函数为:
Figure BDA0001419166660000072
Figure BDA0001419166660000074
Figure BDA0001419166660000075
稳态时对对象进行阶跃响应试验。即对一级网流量阀和二级网循环水泵变频器频率加入10%的阶跃扰动信号,采样时间为10s,获得二级网供回水流量和温差的阶跃响应特征。因为假设二级网供回水流量不变,即动态矩阵AD为:
其中,a11=[a11(1),…a11(N)]T,a12=[a12(1),…a12(N)]T,a21=[a21(1),…a21(N)]T
a22=[a22(1),…a22(N)]T为不同输入对不同输出的N步阶跃响应序列;
建立控制系统的预测模型yPM(k)=yP0(k)+AΔuM(k),通过求解性能指标求出最优控制量Δu(k)进行滚动优化,利用反馈项e(k+1)=y(k+1)-y1(k+1|k)进行误差校正,保证控制精度;式中:yP0(k)为k时刻起控制作用保持不变时,未来P个时刻内对象的输出;AD为动态矩阵;J(k)为k时刻的代价函数;wP(k)为k时刻时未来P个时刻内的设定值序列;Q,R分别为输出误差权矩阵和控制权矩阵;yPM(k)为k时刻起对象受到M个连续控制增量ΔuM(k)作用后在未来P个时刻内的输出;M为控制步长;y(k+1)为对象的在k+1时刻的实际输出;y1(k+1|k)为对象在k时刻受到1个控制增量,经过一步采样时间后在k+1时刻的预测输出;e(k+1)为模型误差。
建立模型误差差分ω(k)的ARMA模型为:
ω(k)+a1ω(k-1)+…anω(k-n)=ξ(k)+c1ξ(k-1)+…+cnξ(k-n) (27)
其中,ω(k),...,ω(k-n)为误差差分序列,多项式参数a1,a2,…an和c1,c2,…cn采用ARMA的参数估计法进行估计。再利用最小方差预报原理,求取误差差分的预报值为:
Figure BDA0001419166660000082
求取预报误差
Figure BDA0001419166660000083
式中e(k)为k时刻的模型误差,用预报误差代替e(k)提高模型精度。
3)将环境参数看作为扰动,仅仅考虑环境参数变化对供热温度的影响,认为环境参数变化对供热水量没有影响。设置前馈补偿回路。对可测环境参数的动态特性进行建模,提取对过程动态影响的因果性信息,利用广义最小二乘方法估计出可测环境参数与系统输出之间的辨识模型,求取辨识模型的阶跃响应特征,得到前馈补偿矩阵,减小环境参数波动对供热品质的影响。
通过求解环境参数扰动模型最终预报误差准则的最小值求解出模型扰动的阶次
Figure BDA0001419166660000084
最终预报误差
Figure BDA0001419166660000085
准则为:
Figure BDA0001419166660000086
其中,
Figure BDA0001419166660000087
为模型输出残差的方差,L为数据长度;
利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:
Figure BDA0001419166660000088
式中:A(z-1),B(z-1)和C(z-1)是稳定的多项式;假设噪声模型C(z-1)未知,考虑用迭代估计的方法求得;
令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次数;
计算输出估计值
Figure BDA0001419166660000091
输入估计值
Figure BDA0001419166660000092
再更新迭代次数i=i+1,式中y(k),u(k)分别为k时刻实际输入和实际输出;
最小二乘法估计出
Figure BDA0001419166660000093
中的参数,式中:Ai(z-1),Bi(z-1)是迭代i次的稳定的多项式;
用估计模型
Figure BDA0001419166660000094
以及各时刻的观测数据,分别为迭代i次、参数为估计值且与Ai(z-1),Bi(z-1)同形的多项式,计算出估计残差
Figure BDA0001419166660000096
计算估计残差平方和
Figure BDA0001419166660000097
及估计残差平方和增量ξ=ξii-1,如果ξ小于增量限值,则结束辨识,否则继续下一步;
对于新得到的噪声模型
Figure BDA0001419166660000098
用最小二乘法估计出参数,得到更新i次后的噪声模型Ci(z-1),返回第b)步,继续按照步骤b)至f)顺序进行计算。
已知供热系统环境参数对供热温度的扰动模型,对该模型进行阶跃响应试验,每隔10s采集一次数据,共采集N组数据组成前馈补偿矩阵
Figure BDA0001419166660000099
依据前馈补偿修改动态矩阵控制算法预测模型为:
yPM(k)=yP0(k)+AΔuM(k)+bΔvP(k) (30)
式中ΔvP(k)为k时刻时未来P个时刻内的扰动增量;
相应的求解性能指标
Figure BDA00014191666600000910
得到控制律;每一采样时刻进行一次优化计算,获得M时域内的最优控制增量序列[Δu(k),Δu(k+1),Δu(k+M-1)]T。只对当前时刻的控制输入量Δu(k),由此Δu(k)计算公式为:
Δu(k)=dT[wP(k)-yP0(k)-bΔvP(k)] (31)
其中,优化矩阵
Figure BDA00014191666600000911
Q,R分别为输出误差权矩阵和控制权矩阵;;
设置相关参数,取采用时间为Ts=10s,预测步长P=40,控制步长M=5,模型长度N=100,某一个控制量的误差权矩阵为
Figure BDA0001419166660000101
某一输出量的误差权矩阵为其中M为控制步长;P为预测步长。最后控制权矩阵
Figure BDA0001419166660000103
输出误差权矩阵因此优化矩阵dT为:
Figure BDA0001419166660000105
计算k时刻的控制量为u(k)=Δu(k)+u(k-1);
输出供热负荷的控制量u(k),计算供热负荷控制系统输出:供热水温和供热流量,再重复1)至3)步,直到系统趋于稳定或者到达计算上限。

Claims (4)

1.一种基于环境参数补偿的实时供热负荷控制方法,其特征在于,包括:
供热负荷预报步骤:采用了基于时间序列的供热负荷预报方法,对供热系统负荷进行预报;
环境参数识别补偿步骤:将环境参数作为扰动,设置前馈补偿回路消除环境参数变化对供热温度的影响;
预测控制步骤:采用预测控制,根据预测的负荷需求,进行操作,提前对大惯性供热负荷系统进行控制调节操作;
所述供热负荷预报步骤包含:
11)针对k时刻负荷序列Q(k),通过求取负荷差分ΔQ(k)获得平稳随机序列,建立相应的ARMA预报模型,预报模型描述为:
A(z-1)ΔQ(k)=C(z-1)ξ(k) (1)
式中:A(z-1)=1+a1z-1+…anz-n,C(z-1)=1+c1z-1+cnz-n;A(z-1)和C(z-1)是稳定的多项式,a1,a2,...,an和c1,c2,...,cn是多项式系数,Δ=1-z-1为差分算子,ξ(k)为均值为零、方差为σ2的白噪声序列,n为多项式阶次,k为当前时刻的标度;
12)式(1)中的参数采用减消记忆的递推最小二乘方法进行估计,求取多项式系数的估计值a1,a2,…an和c1,c2,…cn后,由最小方差预报原理,获取预报值:
Figure FDA0002192360200000011
其中,
Figure FDA0002192360200000012
为在k时刻对k+i时刻负荷差分的预报值;
Figure FDA0002192360200000013
Figure FDA0002192360200000014
为系数为估计值的多项式,下标1代表迭代次数为1;i为相对于当前k时刻的未来时点序号;
13)获取供热负荷预报值Q(k+j|k)为:
Figure FDA0002192360200000015
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值,P为预测步长;
所述环境参数识别补偿步骤包含:
21)利用最终预报误差准则估计环境参数对供热温度扰动模型的阶次,求解下式获取扰动模型阶次
Figure FDA0002192360200000021
Figure FDA0002192360200000022
其中,
Figure FDA0002192360200000023
为模型输出残差的方差,L为数据长度;
22)利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:式中:A(z-1),B(z-1)和C(z-1)是稳定的多项式;
23)假设噪声模型C(z-1)参数未知,用迭代估计的方法求得;
a)令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次
数;
b)计算输出估计值
Figure FDA0002192360200000025
输入估计值再更新迭代次数i=i+1,式中y(k),u(k)分别为k时刻实际输入和实际输出;
c)最小二乘法估计出中的参数,式中:Ai(z-1),Bi(z-1)是迭代i次的稳定的多项式;
d)用估计模型
Figure FDA0002192360200000028
以及各时刻的观测数据,
Figure FDA0002192360200000029
分别为迭代i次、参数为估计值且与Ai(z-1),Bi(z-1)同形的多项式,计算出估计残差
Figure FDA00021923602000000210
e)计算估计残差平方和及估计残差平方和增量ξ=ξii-1,如果ξ小于增量限值,则结束辨识,否则继续下一步;
f)对于新得到的噪声模型
Figure FDA00021923602000000212
用最小二乘法估计出参数,得到更新i次后的噪声模型Ci(z-1),返回第b)步,继续按照步骤b)至f)顺序进行计算;
24)已知供热系统环境参数扰动模型,对该模型进行阶跃响应试验,每隔10s采集一次数据,共采集N组数据组成前馈补偿矩阵N为模型长度;
25)依据前馈补偿修改动态矩阵控制算法预测模型为:
yPM(k)=yP0(k)+ADΔuM(k)+bΔvP(k) (5)
式中:ΔvP(k)为k时刻时未来P个时刻内的扰动增量;
相应地求解优化问题minJ(k)得到控制律为:
Δu(k)=dT[wP(k)-yP0(k)-bΔvP(k)] (6)
其中,优化矩阵
Figure FDA0002192360200000031
Q,R分别为输出误差权矩阵和控制权矩阵;
所述预测控制步骤包含:
31)获取热网对象的控制模型,稳态时对对象进行阶跃响应试验;
32)建立控制系统的预测模型yPM(k)=yP0(k)+ADΔuM(k),通过求解性能指标
Figure FDA0002192360200000032
求出最优控制量Δu(k)进行滚动优化;利用反馈项e(k+1)=y(k+1)-y1(k+1|k)进行误差校正,保证控制精度;式中:yP0(k)为k时刻起控制作用保持不变时,未来P个时刻内对象的输出;AD为动态矩阵;J(k)为k时刻的代价函数;wP(k)为k时刻时未来P个时刻内的设定值序列;Q,R分别为输出误差权矩阵和控制权矩阵;yPM(k)为k时刻起对象受到M个连续控制增量ΔuM(k)作用后在未来P个时刻内的输出;M为控制步长;y(k+1)为对象的在k+1时刻的实际输出;y1(k+1|k)为对象在k时刻受到1个控制增量,经过一步采样时间后在k+1时刻的预测输出;e(k+1)为模型误差;
33)建立模型误差差分ω(k)的ARMA模型为:
ω(k)+a1ω(k-1)+…anω(k-n)=ξ(k)+c1ξ(k-1)+…+cnξ(k-n) (7)
其中,ω(k),…,ω(k-n)为误差差分序列,多项式参数a1,a2,…an和c1,c2,…cn采用ARMA的参数估计法进行估计;再利用最小方差预报原理,求取k时刻对未来的误差差分预报值
Figure FDA0002192360200000033
为:
式中:i、j代表相对k时刻的未来时点序号;
最后求取在k时刻对k+i时刻的预报误差
Figure FDA0002192360200000041
式中e(k)为k时刻的模型误差,用预报误差代替e(k)提高模型精度。
2.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:所述环境参数识别补偿步骤具体为:对可测环境参数的动态特性进行建模,提取对过程动态影响的因果性信息,利用广义最小二乘方法估计出可测环境参数与系统输出之间的辨识模型,求取辨识模型的阶跃响应特征,得到前馈补偿矩阵,减小环境参数波动对供热品质的影响。
3.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:所述供热负荷预报步骤中,以时间为预测对象的唯一影响因素,采用惯性原理对预测对象未来的状态进行推测,最后把预测供热负荷作为供热负荷控制系统的参考值。
4.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:步骤31)中,对一级网流量阀和二级网循环水泵变频器频率加入10%的阶跃扰动信号,采样时间为10s,获得二级网供回水流量和温差的阶跃响应特征。
CN201710886793.5A 2017-09-26 2017-09-26 一种基于环境参数补偿的实时供热负荷控制方法 Expired - Fee Related CN107842908B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710886793.5A CN107842908B (zh) 2017-09-26 2017-09-26 一种基于环境参数补偿的实时供热负荷控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710886793.5A CN107842908B (zh) 2017-09-26 2017-09-26 一种基于环境参数补偿的实时供热负荷控制方法

Publications (2)

Publication Number Publication Date
CN107842908A CN107842908A (zh) 2018-03-27
CN107842908B true CN107842908B (zh) 2020-01-03

Family

ID=61661421

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710886793.5A Expired - Fee Related CN107842908B (zh) 2017-09-26 2017-09-26 一种基于环境参数补偿的实时供热负荷控制方法

Country Status (1)

Country Link
CN (1) CN107842908B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108983595B (zh) * 2018-07-18 2021-04-20 天津大学 一种前馈控制器参数的自动整定方法
CN109297086B (zh) * 2018-09-10 2020-10-09 常州英集动力科技有限公司 热力站负荷分时段滚动预测及自适应矫正方法及系统
CN109214948A (zh) * 2018-09-25 2019-01-15 新智数字科技有限公司 一种电力系统热负荷预测的方法和装置
CN110779075B (zh) * 2019-09-18 2021-09-07 新疆昌吉特变能源有限责任公司 一种热网加热器水位前馈控制方法
CN111998423A (zh) * 2020-08-27 2020-11-27 沈阳工业大学 电蓄热循环风量-水温控制系统及其预测控制方法
CN113739296B (zh) * 2021-09-08 2022-09-06 山东佐耀科技有限公司 基于模型预测控制的空气源热泵负荷水温控制方法及系统
CN113835342B (zh) * 2021-09-18 2024-04-16 国网河北能源技术服务有限公司 一种过热汽温系统的抗扰预测控制方法
CN113659246B (zh) * 2021-10-20 2022-01-25 中国气象科学研究院 适应于极地超低温环境的电池系统及其控温方法
CN114484557B (zh) * 2022-01-22 2022-10-11 天津大学 一种基于目标能耗管控的建筑群供热负荷调控方法
CN114608065B (zh) * 2022-03-29 2022-10-11 天津大学 一种新型热力入口阀位调控方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101021914A (zh) * 2006-03-22 2007-08-22 侯春海 暖通空调负荷预测方法和系统
CN102508466A (zh) * 2011-09-28 2012-06-20 哈尔滨工业大学 建筑供热计量收费与节能监控物联网系统
CN103177391A (zh) * 2013-01-31 2013-06-26 李生林 一种基于三维实景的设施设备监控系统及系统集成方法
CN103322553A (zh) * 2013-07-04 2013-09-25 东南大学 一种火电机组过热汽温多模型扰动估计预测控制方法
CN105570866A (zh) * 2015-12-29 2016-05-11 北京市燃气集团有限责任公司 一种供热用燃气锅炉负荷参数的调节方法及系统
CN105588177A (zh) * 2014-11-18 2016-05-18 天津银箭科技有限公司 基于物联网环境下的智能供热监控应用系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9395712B2 (en) * 2013-05-31 2016-07-19 Patrick Andrew Shiel Building energy usage reduction by automation of optimized plant operation times and sub-hourly building energy forecasting to determine plant faults
US9869481B2 (en) * 2015-01-27 2018-01-16 Patrick Andrew Shiel Method of controlling ventilation and chilling systems to conserve energy in commercial buildings

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101021914A (zh) * 2006-03-22 2007-08-22 侯春海 暖通空调负荷预测方法和系统
CN102508466A (zh) * 2011-09-28 2012-06-20 哈尔滨工业大学 建筑供热计量收费与节能监控物联网系统
CN103177391A (zh) * 2013-01-31 2013-06-26 李生林 一种基于三维实景的设施设备监控系统及系统集成方法
CN103322553A (zh) * 2013-07-04 2013-09-25 东南大学 一种火电机组过热汽温多模型扰动估计预测控制方法
CN105588177A (zh) * 2014-11-18 2016-05-18 天津银箭科技有限公司 基于物联网环境下的智能供热监控应用系统
CN105570866A (zh) * 2015-12-29 2016-05-11 北京市燃气集团有限责任公司 一种供热用燃气锅炉负荷参数的调节方法及系统

Also Published As

Publication number Publication date
CN107842908A (zh) 2018-03-27

Similar Documents

Publication Publication Date Title
CN107842908B (zh) 一种基于环境参数补偿的实时供热负荷控制方法
CN112128841B (zh) 一种基于负荷预测和室温反馈修正的全网平衡调节方法
CN108564204B (zh) 基于最大相关熵准则的最小二乘支持向量机电量预测方法
CN103730006B (zh) 一种短时交通流量的组合预测方法
CN103676651B (zh) 基于状态观测模型的锅炉汽温预测控制方法
CN102444784B (zh) 基于动态矩阵控制的钢铁企业蒸汽管网压力控制系统
CN111898816A (zh) 一种综合能源系统动态状态估计方法
CN107270283B (zh) 一种基于循环流化床机组的多变量约束预测控制方法
CN106707756B (zh) 融合扩张观测器的超临界火电机组机炉协调控制方法
CN103116283A (zh) 一种非自衡对象的动态矩阵控制方法
CN103322553A (zh) 一种火电机组过热汽温多模型扰动估计预测控制方法
CN104270055A (zh) 一种交流伺服系统位置环ip控制器的参数自校正方法
CN106200379B (zh) 一种非自衡对象的分布式动态矩阵控制方法
CN107783941B (zh) 一种基于谷电价格变化的太阳能补偿型电锅炉供热控制方法
CN113091088B (zh) 一种基于双阶段神经网络的锅炉燃烧广义预测控制方法
CN101221437B (zh) 网络信息交换模式下的工业生产全过程优化控制方法
CN106444362A (zh) 一种废塑料裂解炉炉膛温度的分布式pid型预测函数控制方法
CN106444388A (zh) 一种焦炭炉炉膛压力的分布式pid型动态矩阵控制方法
CN103439968A (zh) 存在控制阀粘滞特性的化工过程的控制性能测定方法
CN111998423A (zh) 电蓄热循环风量-水温控制系统及其预测控制方法
CN116430726A (zh) 基于减法聚类和模糊神经网络的涡轮机组控制方法
Wanfeng et al. Adaptive PID controller based on online LSSVM identification
CN110366232B (zh) 用于远程状态估计的传感器传输能量控制方法
CN115271194A (zh) 基于多变量预测控制的综合能源系统双层优化运行方法
Chen et al. Research on Wind and Solar Power Generation Forecasting Based on SARIMA-LSTM 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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200103

Termination date: 20200926