CN107842908B - 一种基于环境参数补偿的实时供热负荷控制方法 - Google Patents
一种基于环境参数补偿的实时供热负荷控制方法 Download PDFInfo
- 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
Links
- 230000007613 environmental effect Effects 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000010438 heat treatment Methods 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 7
- 238000013277 forecasting method Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 230000004044 response Effects 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 10
- 241001123248 Arma Species 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 4
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 3
- 230000001364 causal effect Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 101150036841 minJ gene Proteins 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 230000033228 biological regulation Effects 0.000 abstract description 3
- 238000012546 transfer Methods 0.000 description 8
- 239000008400 supply water Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000008236 heating water Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24D—DOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
- F24D19/00—Details
- F24D19/10—Arrangement or mounting of control or safety devices
- F24D19/1006—Arrangement or mounting of control or safety devices for water heating systems
- F24D19/1009—Arrangement or mounting of control or safety devices for water heating systems for central heating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy 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为当前时刻的标度;
13)获取供热负荷预报值Q(k+j|k)为:
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值,P为预测步长。
所述预测控制步骤包含:
21)获取热网对象的控制模型,稳态时对对象进行阶跃响应试验;
22)建立控制系统的预测模型yPM(k)=yP0(k)+ADΔ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)为模型误差;
23)建立模型误差差分ω(k)的ARMA模型为:
ω(k)+a1ω(k-1)+…anω(k-n)=ξ(k)+c1ξ(k-1)+…+cnξ(k-n) (12)
式中:i、j代表相对k时刻的未来时点序号;
步骤21)中,对一级网流量阀和二级网循环水泵变频器频率加入10%的阶跃扰动信号,采样时间为10s,获得二级网供回水流量和温差的阶跃响应特征。
所述环境参数识别补偿步骤包含:
32)利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:
33)假设噪声模型C(z-1)参数未知,用迭代估计的方法求得;
a)令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次数;
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)
本发明基于环境参数补偿的实时供热负荷控制方法,主要采取以下措施: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时刻预报序列集合的转置N为模型长度。
θ=(ΦTΦ)-1ΦTQN (18)
求取估计值后,由最小方差预报原理,获取预报值:
获取供热负荷预报值Q(k+i|k)为:
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值;P为预测步长;
2)针对供热负荷系统大惯性的特点,采用预测控制进行控制调节。根据预测的负荷需求,提前进行操作,保证用户的用热品质。同时保证控制系统优良的跟踪性能和鲁棒性,增强了对干扰的抑制能力;
获取热网对象的控制模型,设热网供热负荷的传递函数可以由下式表示:
其中y1和y2分别为二次网供回水温差和二级网流量,u1和u2分别为调节一级网供水阀门和调节二级网循环水泵变频器频率。为供热过程耦合模型传递函数矩阵;g11(s)为调节一级网供水阀门对二次网供回水温差的传递函数;g12(s)为调节二级网循环水泵变频器频率对二次网供回水温差的传递函数;g21(s)为调节一级网供水阀门对二级网流量的传递函数;g22(s)为调节二级网循环水泵变频器频率对二级网流量的传递函数。各传递函数为:
稳态时对对象进行阶跃响应试验。即对一级网流量阀和二级网循环水泵变频器频率加入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的参数估计法进行估计。再利用最小方差预报原理,求取误差差分的预报值为:
3)将环境参数看作为扰动,仅仅考虑环境参数变化对供热温度的影响,认为环境参数变化对供热水量没有影响。设置前馈补偿回路。对可测环境参数的动态特性进行建模,提取对过程动态影响的因果性信息,利用广义最小二乘方法估计出可测环境参数与系统输出之间的辨识模型,求取辨识模型的阶跃响应特征,得到前馈补偿矩阵,减小环境参数波动对供热品质的影响。
利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:
令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次数;
依据前馈补偿修改动态矩阵控制算法预测模型为:
yPM(k)=yP0(k)+AΔuM(k)+bΔvP(k) (30)
式中ΔvP(k)为k时刻时未来P个时刻内的扰动增量;
相应的求解性能指标得到控制律;每一采样时刻进行一次优化计算,获得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)
设置相关参数,取采用时间为Ts=10s,预测步长P=40,控制步长M=5,模型长度N=100,某一个控制量的误差权矩阵为某一输出量的误差权矩阵为其中M为控制步长;P为预测步长。最后控制权矩阵输出误差权矩阵因此优化矩阵dT为:
计算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后,由最小方差预报原理,获取预报值:
13)获取供热负荷预报值Q(k+j|k)为:
其中,j为相对于当前时刻的未来时点序号;Q(k+j|k)为在k时刻对k+j时刻的负荷预报值,P为预测步长;
所述环境参数识别补偿步骤包含:
22)利用广义最小二乘方法辨识出环境参数扰动模型,考虑差分方程:式中:A(z-1),B(z-1)和C(z-1)是稳定的多项式;
23)假设噪声模型C(z-1)参数未知,用迭代估计的方法求得;
a)令Ci(z-1)=1,i=0作为下标,表示迭代次数为0,ξ0=100000为最大迭代次
数;
c)最小二乘法估计出中的参数,式中:Ai(z-1),Bi(z-1)是迭代i次的稳定的多项式;
e)计算估计残差平方和及估计残差平方和增量ξ=ξi-ξi-1,如果ξ小于增量限值,则结束辨识,否则继续下一步;
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)
所述预测控制步骤包含:
31)获取热网对象的控制模型,稳态时对对象进行阶跃响应试验;
32)建立控制系统的预测模型yPM(k)=yP0(k)+ADΔ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)为模型误差;
33)建立模型误差差分ω(k)的ARMA模型为:
ω(k)+a1ω(k-1)+…anω(k-n)=ξ(k)+c1ξ(k-1)+…+cnξ(k-n) (7)
式中:i、j代表相对k时刻的未来时点序号;
2.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:所述环境参数识别补偿步骤具体为:对可测环境参数的动态特性进行建模,提取对过程动态影响的因果性信息,利用广义最小二乘方法估计出可测环境参数与系统输出之间的辨识模型,求取辨识模型的阶跃响应特征,得到前馈补偿矩阵,减小环境参数波动对供热品质的影响。
3.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:所述供热负荷预报步骤中,以时间为预测对象的唯一影响因素,采用惯性原理对预测对象未来的状态进行推测,最后把预测供热负荷作为供热负荷控制系统的参考值。
4.根据权利要求1所述的基于环境参数补偿的实时供热负荷控制方法,其特征在于:步骤31)中,对一级网流量阀和二级网循环水泵变频器频率加入10%的阶跃扰动信号,采样时间为10s,获得二级网供回水流量和温差的阶跃响应特征。
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)
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)
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)
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 |
-
2017
- 2017-09-26 CN CN201710886793.5A patent/CN107842908B/zh not_active Expired - Fee Related
Patent Citations (6)
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 |