发明名称 一种基于节能的确定空气源热泵地暖最佳运行工况的方法
摘要 本发明公开了一种基于节能的确定空气源热泵地暖最佳运行工况的方法。目前的方法均不够节能,造成能源的浪费。本发明考虑到环境温度和供水温度、循环水质量流量对系统的影响。包括建立地暖供热的热平衡的模型;建立热泵的能效比和功耗模型;建立地暖水循环的能耗模型方程;建立总的优化模型,依据上述几个模型对其进行优化求解,得到最优控制参数,使得系统运行在最佳的工作点。本发明可降低系统的运行能耗,达到降低运行成本,节能减排的目的。
申请公布号 CN103020481B 申请公布日期 2016.03.30
申请号 CN201210588559.1 申请日期 2012.12.29
申请人 杭州电子科技大学 发明人 江爱朋;丁强;姜周曙;王剑;黄国辉
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州君度专利代理事务所(特殊普通合伙) 33240 代理人 杜军
主权项 一种基于节能的确定空气源热泵地暖最佳运行工况的方法,其特征在于该方法包括以下步骤:步骤1.建立地暖供热的热平衡的模型;房间的单位时间散热损失ΔQ<sub>1</sub>表示为:ΔQ<sub>1</sub>=k<sub>1</sub>A<sub>1</sub>Δt<sub>1</sub>=k<sub>1</sub>A<sub>1</sub>(T<sub>z</sub>‑T<sub>amb</sub>)       (1)k<sub>1</sub>表示房间的平均传热系数,A<sub>1</sub>表示房间的传热面积,Δt<sub>1</sub>表示房间平均温度与环境温度的温差,T<sub>z</sub>表示房间内部平均温度,T<sub>amb</sub>表示环境温度;地暖系统单位时间供给房间的热量ΔQ<sub>2</sub>表示为:ΔQ<sub>2</sub>=k<sub>2</sub>A<sub>2</sub>Δt<sub>2</sub>=k<sub>2</sub>A<sub>2</sub>((T<sub>ws</sub>+T<sub>wr</sub>)/2‑T<sub>z</sub>)       (2)k<sub>2</sub>表示地暖通过地板与房间空气的换热系数,A<sub>2</sub>表示地暖与房间的换热面积,Δt<sub>2</sub>表示房间平均温度与热泵循环水系统平均温度的温差;T<sub>ws</sub>表示热泵系统的出水温度,也是地暖管道的入口水温,T<sub>wr</sub>表示热泵系统的进水温度,也是地暖管道的回水温度;地暖循环水的单位时间热量损失ΔQ<sub>3</sub>为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;Q</mi><mn>3</mn></msub><mo>=</mo><mi>C</mi><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>w</mi><mi>s</mi></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>w</mi><mi>r</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000011.GIF" wi="1681" he="84" /></maths>其中C表示水的质量比热容,<img file="FDA0000824718820000012.GIF" wi="63" he="57" />表示水循环质量流量;要保持房间温度不变,房间的单位时间散热损失和地暖系统单位时间供给房间的热量应该相等,又应该等于地暖循环水的单位时间热量损失;因此,必然满足:ΔQ<sub>2</sub>=ΔQ<sub>3</sub>=ΔQ<sub>1</sub>         (4)地暖通过地板与房间空气的换热系数k<sub>2</sub>除了与换热介质有关外,还与供水循环系统的水的流量有关系,随着流量增大,换热系数增大,并存在增加变缓趋势;因此在一定的流量范围满足:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><msub><mi>b</mi><mn>1</mn></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><msup><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mrow><mn>1</mn><mo>/</mo><mi>t</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000013.GIF" wi="1389" he="85" /></maths>这里<img file="FDA0000824718820000014.GIF" wi="55" he="60" />表示水循环质量流量,单位为kg/s,t表示大于1的数,b<sub>1</sub>和b<sub>2</sub>表示结构系数;步骤2.建立热泵的能效比和功耗模型;热泵的性能采用能效比COP来表示;具体形式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>C</mi><mi>O</mi><mi>P</mi><mo>=</mo><msub><mover><mi>Q</mi><mo>&CenterDot;</mo></mover><mrow><mi>h</mi><mi>p</mi></mrow></msub><msub><mi>/P</mi><mrow><mi>h</mi><mi>p</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000021.GIF" wi="1430" he="94" /></maths>其中COP表示热泵性能系数,<img file="FDA0000824718820000022.GIF" wi="221" he="87" />表示热泵的输出功率,P<sub>hp</sub>表示热泵压缩机功率,单位为瓦;热泵的性能系数表示为环境温度和供水温度的非线性函数,由式(7)所示;COP=Cop<sub>0</sub>+α<sub>1</sub>T<sub>amb</sub>+α<sub>2</sub>T<sub>ws</sub>+α<sub>3</sub>T<sub>amb</sub><sup>2</sup>+α<sub>4</sub>T<sub>ws</sub><sup>2</sup>+α<sub>5</sub>T<sub>amb</sub>T<sub>ws</sub>   (7)Cop<sub>0</sub>以及参数α<sub>1</sub>、α<sub>2</sub>、α<sub>3</sub>、α<sub>4</sub>、α<sub>5</sub>分别表示与热泵性能系数相关的参数;其具体数值根据热泵出厂时候标定的参数确定;在获得了热泵的能效比方程后,根据方程式(8)得到热泵的能耗功率:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>h</mi><mi>p</mi></mrow></msub><mo>=</mo><msub><mover><mi>Q</mi><mo>&CenterDot;</mo></mover><mrow><mi>h</mi><mi>p</mi></mrow></msub><mo>/</mo><mi>C</mi><mi>O</mi><mi>P</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000023.GIF" wi="1060" he="92" /></maths>步骤3.建立地暖水循环的能耗模型方程;水循环能耗P<sub>w</sub>为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>P</mi><mi>w</mi></msub><mo>=</mo><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mo>/</mo><msub><mi>&rho;P</mi><mi>p</mi></msub><mo>/</mo><mi>&eta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000024.GIF" wi="1286" he="97" /></maths>ρ表示水循环系统水的平均密度,ΔP<sub>p</sub>表示循环泵进出口水压之差,η表示机械效率;而ΔP<sub>p</sub>除了与管道的长度、直径有关外,还与水循环系统的流速相关,因此表示为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mi>p</mi></msub><mo>=</mo><mi>&lambda;</mi><mfrac><mrow><msub><mi>F</mi><mi>a</mi></msub><msup><mi>&rho;lv</mi><mn>2</mn></msup></mrow><mi>D</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000025.GIF" wi="1276" he="139" /></maths>F<sub>a</sub>表示摩擦阻力系数,l表示管道长度,D表示管道内径,v表示水的平均流速,λ表示校正参数;根据质量守恒定律,则<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mi>&pi;</mi><mn>4</mn></mfrac><msup><mi>D</mi><mn>2</mn></msup><mi>v</mi><mi>&rho;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000026.GIF" wi="1269" he="120" /></maths>步骤4.建立总的优化模型热泵地暖系统的输出为给房间的供热,输入为热泵压缩机和水循环系统的能耗;在满足房间供暖温度为T<sub>z</sub>的情况下,为使得运行费用最低,总的优化目标为热泵的能耗与水循环泵的能耗最低;因此总的优化目标为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><msub><mi>T</mi><mrow><mi>w</mi><mi>s</mi></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>w</mi><mi>r</mi></mrow></msub><mo>,</mo><mover><mi>m</mi><mo>&CenterDot;</mo></mover></mrow></munder><msub><mi>P</mi><mrow><mi>h</mi><mi>p</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mi>w</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000027.GIF" wi="1361" he="103" /></maths>以<img file="FDA0000824718820000031.GIF" wi="299" he="95" />为优化目标,以建立的方程式(1)‑(11)为约束方程,以0≤P<sub>hp</sub>≤P<sub>hp,max</sub>,0≤P<sub>w</sub>≤P<sub>w,max</sub>,0≤ΔP<sub>p</sub>≤ΔP<sub>p,max</sub>,0≤T<sub>ws</sub>≤60为边界约束条件,以上方程形成总的优化模型,模型中的环境温度、进出口水温、水循环质量流量均为变量;步骤5.对优化命题进行求解,获得热泵地暖系统的进出口水温,循环水量参数;以上的优化模型表示为以下非线性优化命题形式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mfenced open = "" close = "}"><mtable><mtr><mtd><mrow><munder><mi>min</mi><mrow><mi>x</mi><mo>&Element;</mo><msup><mi>R</mi><mi>n</mi></msup></mrow></munder><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>c</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>x</mi><mi>L</mi></msup><mo>&le;</mo><mi>x</mi><mo>&le;</mo><msup><mi>x</mi><mi>U</mi></msup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000032.GIF" wi="1256" he="291" /></maths>其中f:R<sup>n</sup>→R,作为目标函数,c:R<sup>n</sup>→R<sup>m</sup>,表示约束方程,f和c均连续可微;x表示与热泵输出功率和环境温度有关的变量,x<sup>L</sup>和x<sup>U</sup>分别表示变量的上下界约束;采用SQP算法求解以上优化命题;对于式(13)通过求解一系列的QP子问题来逐步逼近其最优解,具体求解过程为:a)x在迭代点x<sub>k</sub>处对式(13)进行泰勒展开,并忽略高次项和目标函数中的常数项,则原命题转为求解以下QP子问题<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mfenced open = "" close = "}"><mtable><mtr><mtd><mrow><munder><mi>min</mi><mrow><msub><mi>d</mi><mi>k</mi></msub><mo>&Element;</mo><msup><mi>R</mi><mi>n</mi></msup></mrow></munder><msub><mi>g</mi><mi>k</mi></msub><msub><mi>d</mi><mi>k</mi></msub><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn><msup><msub><mi>d</mi><mi>k</mi></msub><mi>T</mi></msup><msub><mi>Wd</mi><mi>k</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>c</mi><mi>k</mi></msub><mo>+</mo><msup><msub><mi>A</mi><mi>k</mi></msub><mi>T</mi></msup><msub><mi>d</mi><mi>k</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>x</mi><mi>L</mi></msup><mo>&le;</mo><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msub><mi>d</mi><mi>k</mi></msub><mo>&le;</mo><msup><mi>x</mi><mi>U</mi></msup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824718820000033.GIF" wi="1718" he="526" /></maths>式中d<sub>k</sub>为搜索方向矢量,g<sub>k</sub>和A<sub>k</sub><sup>T</sup>分别表示在x<sub>k</sub>处目标函数f的导数和约束c的雅克比矩阵,W为拉格朗日函数的Hessian阵;通过求解式(14),获得搜索方向矢量的值d<sub>k</sub>;b)令x<sub>k+1</sub>=x<sub>k</sub>+αd<sub>k</sub>,这样就得到了下一个迭代点;α∈(0,1],通过一维搜索方法求取;c)通过计算优化命题(14)的一阶优化条件,判断x<sub>k+1</sub>是否达到最优;如果达到,求解停止,x<sub>k+1</sub>就是最优的结果;否则返回步骤a)继续进行计算;步骤6.根据步骤5获得热泵地暖系统的进出口水温、循环水量参数;从而根据方程式(3)、(7)、(8)、(9)得到热泵的输出功率、热泵能耗功率和循环水泵的功率,通过控制热泵压缩机的频率和循环水泵的频率使热泵和循环水泵达到相应功率,实现整个地暖系统的优化运行,达到节能降耗的目的。
地址 310018 浙江省杭州市下沙高教园区2号大街