CN112417662B - 一种实现集中供热管网***动态水力优化的方法 - Google Patents
一种实现集中供热管网***动态水力优化的方法 Download PDFInfo
- Publication number
- CN112417662B CN112417662B CN202011269983.0A CN202011269983A CN112417662B CN 112417662 B CN112417662 B CN 112417662B CN 202011269983 A CN202011269983 A CN 202011269983A CN 112417662 B CN112417662 B CN 112417662B
- Authority
- CN
- China
- Prior art keywords
- pipe network
- heat
- resistance
- electric valve
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 78
- 238000010438 heat treatment Methods 0.000 claims abstract description 32
- 238000005265 energy consumption Methods 0.000 claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 19
- 238000013461 design Methods 0.000 claims abstract description 6
- 238000005259 measurement Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000000243 solution Substances 0.000 claims description 10
- 238000012544 monitoring process Methods 0.000 claims description 3
- 239000003637 basic solution Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 6
- 230000001105 regulatory effect Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012887 quadratic function Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Thermal Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种实现集中供热管网***动态水力优化的方法,基于非线性规划算法,包括以下步骤:1)对集中供热管网***的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定;其中,管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定;电动阀特性曲线通过在某一供水压力下,测量若干个流量和阀门开度对应值来获得拟合特性曲线;水泵特性曲线通过闭合电动阀调节管网阻力,测量若干个总供回水流量和压差来获得拟合特性曲线;2)建立水力平衡方程;3)以水泵输送能耗最小为目标,以电动阀阻力、热力入口流量、总供回水压差为约束变量。
Description
技术领域
本发明涉及集中供热领域,具体涉及供热管网的水力平衡优化,并且实现随末端负荷变化的动态调节,实现水泵输送能耗最小的目标。
背景技术
水力失调是集中供热***最常见的问题之一,由于实际的供热管网运行工况受到制造和施工、工作条件、环境等多方面的影响,传统的管网难以消除管网之间的水力失衡,使得局部末端偏离设计要求,造成冷热不均的现象,并浪费大量热耗。改善末端管网的水力失调问题,是保证供热***良好运行的重要条件。
供热***在运行中为了达到各个末端的水力及热力平衡,通常安装一些平衡阀来进行调节。史凯[1]通过对济南某实际供热***安装动态及静态压差平衡阀进行管网改造,该调节阀可以在外管网压差变化时仍能保持局部流量稳定,实际运行表明可以显著改善***的水力水平和垂直失调的问题;Jiang[2]采用了一种远程控制***进行水力平衡调节,采用远程调节控制阀和流量计反馈的方式,有效地保证了供暖***水力平衡,节约供热能耗;陈拥政[3]提出了一种通过二次侧的温度控制阀进行回水温度调节***,该温度控制阀可以在外管网水力工况变动时保证局部供回水温差恒定,从而实现二次侧供热自动平衡。以上进行水力平衡的手段都是采用静态调节的方法,这也是大多数供热***进行水力调节的方法,但缺点是无法根据负荷变化进行动态调节,局部的水力工况变化时往往会引起整体的水力不平衡。
随着智慧热网的发展,实现随用热负荷变化的动态调节对供热***的节能降耗有重要作用,为了能够使***在负荷变化时又能保证末端管网的水力平衡,需要进行供热管网的动态水力调节。Wang[4]基于水泵变频,提出一种供热***运行的最优压力控制策略,与传统的恒定压力运行相比,最优压力控制策略在整个供暖季可以降低水泵能耗14.6%;Lu[5]将流量控制阀与变频水泵相结合,在满足各个末端的流量需求调节下,利用线性规划的方法优化了水泵频率和阀门开度,与水泵定频运行相比节能效果显著;周旋[6]基于某一换热站下的多个热力入口进行了动态水力建模,通过遗传算法对水泵频率和调节阀开度进行了优化,结果表明动态调节法可以进一步降低输配***的运行能耗约12.27%。以上通过供热管网的动态水力调节,可以保证***的动态水力平衡,且优化水泵运行工况。
动态水力调节的核心算法是寻优算法,如规划算法、遗传算法、粒子群算法等,对于一般较为复杂的问题往往是非线性的,非线性规划算法可以很好地适应复杂问题的优化求解。非线性规划算法在目标和约束条件中含有非线性函数,该问题一般通过迭代法进行求解,该类算法在选址、电力***电源规划,多能源互补优化等方面有较多应用[7-9]。
参考文献:
[1]史凯,花博,侯家涛,等.基于动态压差平衡阀和静态平衡阀的供热二次管网水力平衡案例分析[J].区域供热,2018(06):81-88.
[2]Jiang J Y,Yin Y B,Lin N T.The Hydraulic Stability of HeatingNetwork with Remote Measurement and Control System[J].Applied Mechanics&Materials,2013,385-386:927-930.
[3]陈拥政,姚强,于大勇,等.基于回水温度的二次侧供热自动平衡调节方法[P]:
[4]Wang N,You S,Wang Y,et al.Hydraulic resistance identification andoptimal pressure control of district heating network[J].Energy and Buildings,2018,170:83-94.
[5]Lu C Y,Ma J F,Wang X H,et al.Hydraulic Balance OptimizationCalculation and Experimental Research on Model of Fluid Network System[J].Applied Mechanics and Materials,2013,300-301:681-688.
[6]周璇,刘国强,王晓佩,等.基于动态水力平衡的集中供热***二次管网节能优化调节方法研究[J].中南大学学报(自然科学版),2019,50(01):220-233.
[7]Olcay G,Atil K,Murat Y D,et al.Circular eco-industrial park designinspired by nature:An integrated non-linear optimization,location,and foodweb analysis.[J].Journal of environmental management,2020,270.
[8]Wu N,Zhan X,Zhu X,et al.Analysis of biomass polygenerationintegrated energy system based on a mixed-integer nonlinear programmingoptimization method[J].Journal of Cleaner Production,2020,271.
[9]苏小珂,肖人峰,陈瑞,等.电力工业中电源规划的研究[J].科学技术创新,2020(25):19-20.
发明内容
本发明的目的是为了克服现有技术中的不足,提供一种利用非线性规划算法实现集中供热管网***动态水力优化的方法,以水泵能耗最小为目标,以调节阀开度、流量为变量,通过非线性规划寻优算法进行水力工况优化。该方法可基于供热***各个末端的用热需求,实现复杂管网***的水力动态优化,并降低能源站内运行能耗,且具有较高的计算效率。
本发明的目的是通过以下技术方案实现的:
一种实现集中供热管网***动态水力优化的方法,基于非线性规划算法,包括以下步骤:
1)对集中供热管网***的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定;其中,管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定;电动阀特性曲线通过在某一供水压力下,测量若干个流量和阀门开度对应值来获得拟合特性曲线;水泵特性曲线通过闭合电动阀调节管网阻力,测量若干个总供回水流量和压差来获得拟合特性曲线;
2)建立水力平衡方程;水力平衡方程需同时满足流量守恒和回路压降和为零;将集中供热管网***中的供热管网类比电路结构,基于基尔霍夫电流定律,供热管网的支路流量等于该支路所有热力入口流量之和,写成矩阵的形式为:
G=AQ
G:支路流量;A:系数矩阵;Q:热力入口流量;
基于基尔霍夫电压定律,在供热管网中应满足闭合回路压降的代数和为0,在供热管网中,经过总供回水管的管路压降方程的矩阵形式为:
P-S(G*G)-(R+Rv)(Q*Q)=0
P:总供回水压差;S:沿程阻力系数;R:热力入口阻力系数;Rv:电动阀阻力系数;*:矩阵对应元素乘积运算;
3)以水泵输送能耗最小为目标,以电动阀阻力、热力入口流量、总供回水压差为约束变量,建立非线性规划求解模型如下:
上式中,目标函数除了满足水力平衡方程外,电动阀阻力、热力入口流量和总供回水压差分别满足最小阀阻Rv0、最小流量Q0和最小压差0约束;由于线性规划求解模型中含有流量约束变量,对于未安装电动阀进行控制的热力入口,线性规划求解模型同样能够进行求解;
4)优化求得水泵的频率及电动阀的开度,并通过监控平台对各水泵和电动阀门发送指令执行。
进一步的,步骤(3)中非线性规划求解模型通过迭代法进行求解,迭代法的基本求解方法为:首先给定一个初始点x0∈Rn,安照迭代规则生成一个点列{xk},当{xk}是有限点列时,最后一个点便是这个非线性规划的最优解;设xk∈Rn是第k轮迭代点,xk+1∈Rn是第k+1轮迭代点,则xk和xk+1的关系为:
xk+1=xk+tkpk
上式中,pk是xk沿xk+1方向的单位向量,|pk|=1;tk是pk沿方向的步长因子,该式表明了xk的下一轮迭代点xk+1是沿pk方向的tk长度,即确定了下一轮迭代点的位置;
使用迭代法求解非线性规划的关键在于如何确定每一轮搜索的步长和方向,在求解非线性规划之前引入下降方向和可行方向两个概念,下降方向指函数在一点处的下降方向是使函数值减少或下降的方向;可行方向指函数在一点处关于某区域的可行方向是使这个方向上存在可行点的方向;
由此可知,使用基本迭代格式求解非线性规划问题的步骤如下:
(301)选取初始数据;选择初始迭代点x0,令迭代次数k=0;
(302)构造搜索方向;对于待求解的非线性规划问题,构造目标函数在点xk处的可行下降方向作为搜索方向pk;
(303)寻求搜索步长;以xk为起点沿搜索方向pk寻求合适的步长tk,使规划问题的目标函数有意义的下降;
(304)求解下一个迭代点;下一个迭代点xk+1由xk+1和xk的关系式(xk+1=xk+tkpk)求出,并令k=k+1再次从第(302)步开始执行;直到xk+1满足给定的终止条件后停止迭代,输出近似最优解xk+1。
进一步的,步骤4)中的具体如下:水泵特性曲线可表示为二次函数的形式为:P=k1Q2+k2Q+k3,其中k1,k2,k3为拟合系数,能够通过实验测定水泵多组的流量和压头数据拟合获得;当压头P和流量Q确定后,通过求解水泵频率f得:
与现有技术相比,本发明的技术方案所带来的有益效果是:
1)将非线性规划优化算法应用到集中供热***的水力平衡优化过程中,可以保证***中水泵输送能耗最小,对于部分热力入口未安装电动阀的***也可进行优化求解,因此该方法适用于较复杂的管网***的水力工况优化。
2)将热力入口流量作为非线性规划算法的约束变量,对于具有未管控入口的***仍然可以进行水力优化,尽可能地降低水力失调程度。这样不仅可以降低水泵电耗,还可以进一步降低***的供热能耗。
3)将参数标定和非线性规划优化算法相结合。具体表现为:首先通过压力表和流量计测量的实际数据对供热管网进行标定,之后再利用非线性规划优化算法对***的水力平衡问题进行优化,保证了模型具有较高的计算精度。
4)将非线性规划优化算法应用到集中供热***的水力平衡优化过程中,可以实现水力的动态优化,即基于末端的负荷变化进行流量调节,尤其适用于智能热网***的智能管控。
5)工程量小,易于实施。本发明仅对能源站及热力入口管道安装热量表、电动阀和压力变送器即可实现动态水力调节,不涉及大规模的管道改造及土木工程等。
6)该方法适用性广,不仅适用于供热的直供管网和二次管网的水力优化调节,其他任意具有变频水泵和调节阀的管网***也同样适用。
附图说明
图1是本发明涉及的集中供热管网***的结构示意图。
图2是集中供热管网水力动态优化流程图。
图3是某高校供热***供热管网分布及简化模型。
图4是各入口模拟平均开度和实际平均开度图。
图5是模拟总流量值与实际总流量图。
图6是动态优化后的水泵频率变化图。
图7是动态优化前后的水泵能耗变化图。
图8是动态优化前后的总管流量变化图。
附图标记:1-变频水泵,2-调节阀,3-热量表,4-压力表,5-锅炉/换热器
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一种基于非线性规划算法实现集中供热管网***动态水力优化的方法,其中集中供热管网***的结构示意图见图1,包括变频水泵1,调节阀2,热量表3,压力表4和锅炉/换热器5。
通过以下四个步骤,当输入给定的各热力入口流量后,可以优化求得水泵的频率及电动阀的开度,并通过监控平台各水泵和电动阀门发送指令执行,具体流程如图2所示。
1)模型参数标定。在通过该方法进行水力工况优化之前,需要对***的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定。其中,管网及热用户阻力可结合管网设计图、压力表和热量表实测参数标定;电动阀特性曲线可在某一供水压力下,测量多个流量和阀门开度对应值来获得拟合特性曲线;水泵特性曲线可通过闭合电动阀调节管网阻力,测量多个总供回水流量和压差来获得拟合特性曲线。
2)建立水力平衡方程。供热管网和电路具有相似的结构,供热管网中的压力、流量、阻力系数等参数可类比于电路中的电压、电流、电阻等参数,电学中的基尔霍夫定律在供热管网中同样适用,因此水力平衡方程需要同时满足流量守恒和回路压降和为零。基于基尔霍夫电流定律,供热管网的支路流量等于该支路所有热力入口流量之和,写成矩阵的形式为:
G=AQ
G:支路流量;A:系数矩阵;Q:热力入口流量。
基于基尔霍夫电压定律,在供热管网中应满足闭合回路压降的代数和为0,在供热管网中,经过总供回水管的管路压降方程的矩阵形式为:
P-S(G*G)-(R+Rv)(Q*Q)=0
P:总供回水压差;S:沿程阻力系数;R:热力入口阻力系数;Rv:电动阀阻力系数;*:矩阵对应元素乘积运算。由于支路流量G是矩阵的形式,G*G可表示为一系列支路流量按照一定顺序排列的流量向量
3)优化求解。以水泵输送能耗最小为目标,以电动阀阻力、热力入口流量、总供回水压差为约束变量,可建立非线性规划求解模型如下:
上式中,目标函数除了满足水力平衡方程外,电动阀阻力、热力入口流量和总供回水压差分别满足最小阀阻Rv0、最小流量Q0和最小压差0约束。由于规划模型中含有流量约束变量,对于未安装电动阀进行控制的热力入口,该模型同样可进行求解。
4)确定水泵频率和阀门开度。水泵特性曲线可表示为二次函数的形式为:P=k1Q2+k2Q+k3,其中k1,k2,k3为拟合系数,能够通过实验测定水泵多组的流量和压头数据拟合获得;当压头P和流量Q确定后,通过求解水泵频率f得:
以天津市某高校锅炉集中供热管网***为例,该供热管网***通过燃气锅炉为直连的38个热力入口供热,经过供热改造后,该***已对28个热力入口安装电动阀进行了管控,并对能源站和28个热力入口进行了热计量。***改造后在2019-2020供暖季对管控的热力入口进行变流量运行,节能效果明显,但站内水泵依然进行定频运行,且各个管控入口的阀门开度不为全开,未管控入口的供热占比高达39%(设计负荷占比为19%)。为了解决以上问题,通过该方法对该***进行水力工况的优化,图3是该高校供热***的管网分布以及经过模型简化后的示意图,基于简化模型建立水力平衡方程:将未管控1-9#入口简化为一个未管控入口,22#入口为一个未管控入口,因此简化模型只有两个未管控入口;另外还有一个热力入口进行了分支路管控,因此看做增加一个管控入口。该模型的系数矩阵A为:
系数矩阵B为:
非线性规划模型为:
建立好模型后,可通过Matlab、Python、Excel等具有编程功能的软件进行计算。在优化计算前,首先基于2019-2020供暖季实际运行数据进行模型验证,模型验证选取了11.20、12.5、12.10、12.21、1.5、1.25和2.15共7天的各入口流量和阀门开度数据,对模型进行了定频水力计算,即输入当天水泵实际台数和运行频率,以及各个入口实际流量,输出为各个入口阀门开度和总流量。经过计算,各个热力入口的平均开度模拟值较实际值整体偏小,偏差在1.6~18.9%,如图4所示。另外对模拟的总流量和实际的总流量进行了比较,如图5所示,模拟值和实际值基本吻合,除了个别时段偏差较大外,整体误差一般不超过5%(置信度98.2%)。该模型能够较为合理地表征实际管网工况,并通过模型进行动态水力计算。
2019-2020供热季水泵调节范围大致为38-45Hz,若基于各个末端的逐时负荷变化进行逐时变频调节,水泵频率在18-45Hz之间。图6是优化后的水泵变频变化,图中可以看到,在供暖季中期(12月中旬至1月中旬)水泵频率变化范围较大,这是因为在该阶段各个入口不同时间变化较大;供暖季前期水泵变化范围较小,这是因为该阶段处于初期调试,各个入口时间流量变化小;寒假期间水泵频率变化基本趋于平稳,这是因为该阶段实际全天无调节。
图7是该调节方案下的水泵能耗变化(站内实际流量11月和12月底有缺失)。由水泵能耗曲线可知,采用变频调节水泵能耗较实际能耗显著降低,另外在12月30日至1月21日增加水泵后,实际水泵能耗会骤增,而采用变频调节后能耗仍平稳在较低水平,因此采用变频调节可以减少由水泵台数的增加导致的能耗骤增。通过模型优化后的水泵逐时平均能耗降低38.2%,因此站内水泵仍具有较大的节能潜力。
图8是该调节方案下的站内总流量变化,由于水泵频率降低,未管控入口减少。通过模型优化的总流量总体低于实际流量,逐时平均流量较实际降低9.0%,即供热能耗降低9.0%,因此通过水力优化可以降低未管控入口的水力失调,减少供热能耗损失。
本发明并不限于上文描述的实施方式。以上对具体实施方式的描述旨在描述和说明本发明的技术方案,上述的具体实施方式仅仅是示意性的,并不是限制性的。在不脱离本发明宗旨和权利要求所保护的范围情况下,本领域的普通技术人员在本发明的启示下还可做出很多形式的具体变换,这些均属于本发明的保护范围之内。
Claims (2)
1.一种实现集中供热管网***动态水力优化的方法,基于非线性规划算法,其特征在于,包括以下步骤:
1)对集中供热管网***的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定;其中,管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定;电动阀特性曲线通过在某一供水压力下,测量若干个流量和阀门开度对应值来获得拟合特性曲线;水泵特性曲线通过闭合电动阀调节管网阻力,测量若干个总供回水流量和压差来获得拟合特性曲线;
2)建立水力平衡方程;水力平衡方程需同时满足流量守恒和回路压降和为零;将集中供热管网***中的供热管网类比电路结构,基于基尔霍夫电流定律,供热管网的支路流量等于该支路所有热力入口流量之和,写成矩阵的形式为:
G=AQ
G:支路流量;A:系数矩阵;Q:热力入口流量;
基于基尔霍夫电压定律,在供热管网中应满足闭合回路压降的代数和为0,在供热管网中,经过总供回水管的管路压降方程的矩阵形式为:
P-S(G*G)-(R+Rv)(Q*Q)=0
P:总供回水压差;S:沿程阻力系数;R:热力入口阻力系数;Rv:电动阀阻力系数;*:矩阵对应元素乘积运算;
3)以水泵输送能耗最小为目标,以电动阀阻力、热力入口流量、总供回水压差为约束变量,建立非线性规划求解模型如下:
上式中,目标函数除了满足水力平衡方程外,电动阀阻力、热力入口流量和总供回水压差分别满足最小阀阻Rv0、最小流量Q0和最小压差0约束;由于线性规划求解模型中含有流量约束变量,对于未安装电动阀进行控制的热力入口,线性规划求解模型同样能够进行求解;
4)优化求得水泵的频率及电动阀的开度,并通过监控平台对各水泵和电动阀门发送指令执行。
2.根据权利要求1所述一种实现集中供热管网***动态水力优化的方法,其特征在于,步骤(3)中非线性规划求解模型通过迭代法进行求解,迭代法的基本求解方法为:首先给定一个初始点x0∈Rn,按照迭代规则生成一个点列{xk},当{xk}是有限点列时,最后一个点便是这个非线性规划的最优解;设xk∈Rn是第k轮迭代点,xk+1∈Rn是第k+1轮迭代点,则xk和xk +1的关系为:
xk+1=xk+tkpk
上式中,pk是xk沿xk+1方向的单位向量,|pk|=1;tk是pk沿方向的步长因子,该式表明了xk的下一轮迭代点xk+1是沿pk方向的tk长度,即确定了下一轮迭代点的位置;
使用迭代法求解非线性规划的关键在于如何确定每一轮搜索的步长和方向,在求解非线性规划之前引入下降方向和可行方向两个概念,下降方向指函数在一点处的下降方向是使函数值减少或下降的方向;可行方向指函数在一点处关于某区域的可行方向是使这个方向上存在可行点的方向;
由此可知,使用基本迭代格式求解非线性规划问题的步骤如下:
(301)选取初始数据;选择初始迭代点x0,令迭代次数k=0;
(302)构造搜索方向;对于待求解的非线性规划问题,构造目标函数在点xk处的可行下降方向作为搜索方向pk;
(303)寻求搜索步长;以xk为起点沿搜索方向pk寻求合适的步长tk,使规划问题的目标函数有意义的下降;
(304)求解下一个迭代点;下一个迭代点xk+1由xk+1和xk的关系式(xk+1=xk+tkpk)求出,并令k=k+1再次从第(302)步开始执行;直到xk+1满足给定的终止条件后停止迭代,输出近似最优解xk+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011269983.0A CN112417662B (zh) | 2020-11-13 | 2020-11-13 | 一种实现集中供热管网***动态水力优化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011269983.0A CN112417662B (zh) | 2020-11-13 | 2020-11-13 | 一种实现集中供热管网***动态水力优化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112417662A CN112417662A (zh) | 2021-02-26 |
CN112417662B true CN112417662B (zh) | 2022-07-15 |
Family
ID=74831035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011269983.0A Active CN112417662B (zh) | 2020-11-13 | 2020-11-13 | 一种实现集中供热管网***动态水力优化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112417662B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113076707A (zh) * | 2021-04-25 | 2021-07-06 | 西安建筑科技大学 | 一种风机盘管***的热能输入控制方法、介质及设备 |
CN113158476B (zh) * | 2021-04-27 | 2023-06-09 | 唐山晓世科技有限公司 | 供热管网水力平衡调节计算机方法 |
CN113251321A (zh) * | 2021-05-28 | 2021-08-13 | 华能(广东)能源开发有限公司海门电厂 | 一种识别蒸汽热网爆管位置的方法和*** |
CN113958991A (zh) * | 2021-09-29 | 2022-01-21 | 西安沣东华能热力有限公司 | 一种分布式变频水泵控制方法 |
CN114857743A (zh) * | 2022-01-25 | 2022-08-05 | 西安建筑科技大学 | 一种基于商场分区负荷预测的末端阀门优化控制方法及*** |
CN115423295B (zh) * | 2022-08-30 | 2023-12-05 | 青岛艾德森物联科技有限公司 | 一种空调水***的双平衡状态监测的评估方法及*** |
CN116717839A (zh) * | 2023-08-10 | 2023-09-08 | 陕西拓普索尔电子科技有限责任公司 | 供热控制方法、控制装置和供热*** |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102162443A (zh) * | 2011-03-24 | 2011-08-24 | 杭州电子科技大学 | 一种多泵站城市排水管网***的调速节能控制方法 |
RO127420A2 (ro) * | 2010-11-29 | 2012-05-30 | Universitatea Tehnică ''gheorghe Asachi'' Din Iaşi | Metodă de optimizare a sistemelor hidraulice pentru alimentări cu apă |
CN102542081A (zh) * | 2010-12-21 | 2012-07-04 | 北京工业大学 | 基于拓扑结构演化的给水管网分区方法 |
CN103064285A (zh) * | 2012-12-29 | 2013-04-24 | 杭州电子科技大学 | 一种基于模型的热泵供暖多目标优化控制方法 |
DE102012020750A1 (de) * | 2012-10-23 | 2014-04-24 | Michael Sparn | Verfahren zur Optimierung eines thermischen und eines hydraulischen Abgleichs in einer Temperieranlage |
KR101562490B1 (ko) * | 2015-04-14 | 2015-10-23 | 고려대학교 산학협력단 | 메타휴리스틱 수원 추적 기반 관망 설계 방법 |
CN105674390A (zh) * | 2016-01-22 | 2016-06-15 | 张凡 | 一种集中供热***动态水力平衡调节方法 |
CN111125938A (zh) * | 2020-01-15 | 2020-05-08 | 华南理工大学 | 基于次优算法的大型中央空调冷冻水管网优化设计方法 |
CN111210152A (zh) * | 2019-12-30 | 2020-05-29 | 清华大学 | 一种排水***调度方法及装置 |
CN111815068A (zh) * | 2020-07-22 | 2020-10-23 | 谢仕炜 | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107609680B (zh) * | 2017-08-22 | 2020-11-20 | 天津大学 | 一种多热源环状集中供热管网水力工况优化调度方法 |
CN109636037B (zh) * | 2018-12-12 | 2023-02-14 | 石家庄华电供热集团有限公司 | 一种基于pso算法的多热源环状热网水力优化调度方法 |
-
2020
- 2020-11-13 CN CN202011269983.0A patent/CN112417662B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RO127420A2 (ro) * | 2010-11-29 | 2012-05-30 | Universitatea Tehnică ''gheorghe Asachi'' Din Iaşi | Metodă de optimizare a sistemelor hidraulice pentru alimentări cu apă |
CN102542081A (zh) * | 2010-12-21 | 2012-07-04 | 北京工业大学 | 基于拓扑结构演化的给水管网分区方法 |
CN102162443A (zh) * | 2011-03-24 | 2011-08-24 | 杭州电子科技大学 | 一种多泵站城市排水管网***的调速节能控制方法 |
DE102012020750A1 (de) * | 2012-10-23 | 2014-04-24 | Michael Sparn | Verfahren zur Optimierung eines thermischen und eines hydraulischen Abgleichs in einer Temperieranlage |
CN103064285A (zh) * | 2012-12-29 | 2013-04-24 | 杭州电子科技大学 | 一种基于模型的热泵供暖多目标优化控制方法 |
KR101562490B1 (ko) * | 2015-04-14 | 2015-10-23 | 고려대학교 산학협력단 | 메타휴리스틱 수원 추적 기반 관망 설계 방법 |
CN105674390A (zh) * | 2016-01-22 | 2016-06-15 | 张凡 | 一种集中供热***动态水力平衡调节方法 |
CN111210152A (zh) * | 2019-12-30 | 2020-05-29 | 清华大学 | 一种排水***调度方法及装置 |
CN111125938A (zh) * | 2020-01-15 | 2020-05-08 | 华南理工大学 | 基于次优算法的大型中央空调冷冻水管网优化设计方法 |
CN111815068A (zh) * | 2020-07-22 | 2020-10-23 | 谢仕炜 | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112417662A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112417662B (zh) | 一种实现集中供热管网***动态水力优化的方法 | |
CN105910169B (zh) | 基于机理模型预测控制的城市供热***热网调节方法及*** | |
Yabanova et al. | Development of ANN model for geothermal district heating system and a novel PID-based control strategy | |
Wang et al. | Hydraulic resistance identification and optimal pressure control of district heating network | |
CN111222213B (zh) | 一种热力网络动态仿真方法及装置 | |
CN106958855B (zh) | 供热***水力平衡的模型预测控制方法及*** | |
Wang et al. | Optimization modeling for smart operation of multi-source district heating with distributed variable-speed pumps | |
CN110793094B (zh) | 一种城市集中供热***综合灵活调度方法和*** | |
Kuosa et al. | Static study of traditional and ring networks and the use of mass flow control in district heating applications | |
Zheng et al. | Thermo-hydraulic coupled simulation and analysis of a real large-scale complex district heating network in Tianjin | |
CN113190999A (zh) | 供热管网流量调节提升风电消纳的电热协调方法及装置 | |
Zhou et al. | Hydraulic modeling of double-source and ring-shaped heating networks | |
CN108954484B (zh) | 一种供热管网***及自适应调控方法 | |
Lazarević et al. | A novel approach to real-time modelling of the district heating substation system using LabVIEW | |
CN103235512B (zh) | 一种发电机组运行的方法 | |
CN108954491A (zh) | 一种光热中温补偿型电锅炉供热***的控制方法 | |
Xu et al. | Discussions on the real potential of district heating networks in improving wind power accommodation with temperature feedback as one consideration | |
Liu et al. | A thermal-hydraulic coupled simulation approach for the temperature and flow rate control strategy evaluation of the multi-room radiator heating system | |
CN110222362B (zh) | 一种通过隧洞和调压井的多机形式建立一管多机微分方程模型的方法 | |
Li et al. | Novel effective room temperature-based predictive feedback control method for large-scale district heating substation | |
Zhu et al. | Operation optimization research of circulating cooling water system based on superstructure and domain knowledge | |
CN111783309A (zh) | 基于内部守恒的蒸汽供热网络动态仿真方法 | |
CN116611706A (zh) | 基于多能源主体的动态碳排放因子测算方法 | |
Tarasov et al. | Mathematical modeling of the forecast and standby heating modes | |
Dobos et al. | Dynamic model and control of heat exchanger networks for district heating |
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 |