一种实现集中供热管网系统动态水力优化的方法技术方案

本发明专利技术公开一种实现集中供热管网系统动态水力优化的方法，基于非线性规划算法，包括以下步骤：1)对集中供热管网系统的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定；其中，管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定；电动阀特性曲线通过在某一供水压力下，测量若干个流量和阀门开度对应值来获得拟合特性曲线；水泵特性曲线通过闭合电动阀调节管网阻力，测量若干个总供回水流量和压差来获得拟合特性曲线；2)建立水力平衡方程；3)以水泵输送能耗最小为目标，以电动阀阻力、热力入口流量、总供回水压差为约束变量。回水压差为约束变量。回水压差为约束变量。

【技术实现步骤摘要】
Network with Remote Measurement and Control System[J].Applied Mechanics&Materials,2013,385-386:927-930.
[0009][3]陈拥政,姚强,于大勇,等.基于回水温度的二次侧供热自动平衡调节方法[P]:
[0010][4]Wang N,You S,Wang Y,et al.Hydraulic resistance identification and optimal pressure control of district heating network[J].Energy and Buildings,2018,170:83-94.
[0011][5]Lu C Y,Ma J F,Wang X H,et al.Hydraulic Balance Optimization Calculation and Experimental Research on Model of Fluid Network System[J].Applied Mechanics and Materials,2013,300-301:681-688.
[0012][6]周璇,刘国强,王晓佩,等.基于动态水力平衡的集中供热系统二次管网节能优化调节方法研究[J].中南大学学报(自然科学版),2019,50(01):220-233.
[0013][7]Olcay G,Atil K,Murat Y D,et al.Circular eco-industrial park design inspired by nature:An integrated non-linear optimization,location,and food web analysis.[J].Journal of environmental management,2020,270.
[0014][8]Wu N,Zhan X,Zhu X,et al.Analysis of biomass polygeneration integrated energy system based on a mixed-integer nonlinear programming optimization method[J].Journal of Cleaner Production,2020,271.
[0015][9]苏小珂,肖人峰,陈瑞,等.电力工业中电源规划的研究[J].科学技术创新,2020(25):19-20.

技术实现思路

[0016]本专利技术的目的是为了克服现有技术中的不足，提供一种利用非线性规划算法实现集中供热管网系统动态水力优化的方法，以水泵能耗最小为目标，以调节阀开度、流量为变量，通过非线性规划寻优算法进行水力工况优化。该方法可基于供热系统各个末端的用热需求，实现复杂管网系统的水力动态优化，并降低能源站内运行能耗，且具有较高的计算效率。
[0017]本专利技术的目的是通过以下技术方案实现的：
[0018]一种实现集中供热管网系统动态水力优化的方法，基于非线性规划算法，包括以下步骤：
[0019]1)对集中供热管网系统的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定；其中，管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定；电动阀特性曲线通过在某一供水压力下，测量若干个流量和阀门开度对应值来获得拟合特性曲线；水泵特性曲线通过闭合电动阀调节管网阻力，测量若干个总供回水流量和压差来获得拟合特性曲线；
[0020]2)建立水力平衡方程；水力平衡方程需同时满足流量守恒和回路压降和为零；将集中供热管网系统中的供热管网类比电路结构，基于基尔霍夫电流定律，供热管网的支路流量等于该支路所有热力入口流量之和，写成矩阵的形式为：
[0021]G＝AQ
[0022]G：支路流量；A：系数矩阵；Q：热力入口流量；
[0023]基于基尔霍夫电压定律，在供热管网中应满足闭合回路压降的代数和为0，在供热管网中，经过总供回水管的管路压降方程的矩阵形式为：
[0024]P-S(G*G)-(R+R
v
)(Q*Q)＝0
[0025]P：总供回水压差；S：沿程阻力系数；R：热力入口阻力系数；R
v
：电动阀阻力系数；*：矩阵对应元素乘积运算；
[0026]3)以水泵输送能耗最小为目标，以电动阀阻力、热力入口流量、总供回水压差为约束变量，建立非线性规划求解模型如下：
[0027][0028]上式中，目标函数除了满足水力平衡方程外，电动阀阻力、热力入口流量和总供回水压差分别满足最小阀阻R
v0
、最小流量Q0和最小压差0约束；由于线性规划求解模型中含有流量约束变量，对于未安装电动阀进行控制的热力入口，线性规划求解模型同样能够进行求解；
[0029]4)优化求得水泵的频率及电动阀的开度，并通过监控平台对各水泵和电动阀门发送指令执行。
[0030]进一步的，步骤(3)中非线性规划求解模型通过迭代法进行求解，迭代法的基本求解方法为：首先给定一个初始点x0∈R
n
，安照迭代规则生成一个点列{x
k
}，当{x
k
}是有限点列时，最后一个点便是这个非线性规划的最优解；设x
k
∈R
n
是第k轮迭代点，x
k+1
∈R
n
是第k+1轮迭代点，则x
k
和x
k+1
的关系为：
[0031]x
k+1
＝x
k
+t
k
p
k
[0032]上式中，p
k
是x
k
沿x
k+1
方向的单位向量，|p
k
|＝1；t
k
是p
k
沿方向的步长因子，该式表明了x
k
的下一轮迭代点x
k+1
是沿p
k
方向的t
k
长度，即确定了下一轮迭代点的位置；
[0033]使用迭代法求解非线性规划的关键在于如何确定每一轮搜索的步长和方向，在求解非线性规划之前引入下降方向和可行方向两个概念，下降方向指函数在一点处的下降方向是使函数值减少或下降的方向；可行方向指函数在一点处关于某区域的可行方向是使这个方向上存在可行点的方向；
[0034]由此可知，使用基本迭代格式求解非线性规划问题的步骤如下：
[0035](301)选取初始数据；选择初始迭代点x0，令迭代次数k＝0；
[0036](302)构造搜索方向；对于待求解的非线性规划问题，构造目标函数在点x
k
处的可行下降方向作为搜索方向p
k
；
[0037](303)寻求搜索步本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种实现集中供热管网系统动态水力优化的方法，基于非线性规划算法，其特征在于，包括以下步骤：1)对集中供热管网系统的管网阻力、热用户阻力、电动阀特性曲线及水泵特性曲线进行标定；其中，管网及热用户阻力通过结合管网设计图、压力表和热量表实测参数进行标定；电动阀特性曲线通过在某一供水压力下，测量若干个流量和阀门开度对应值来获得拟合特性曲线；水泵特性曲线通过闭合电动阀调节管网阻力，测量若干个总供回水流量和压差来获得拟合特性曲线；2)建立水力平衡方程；水力平衡方程需同时满足流量守恒和回路压降和为零；将集中供热管网系统中的供热管网类比电路结构，基于基尔霍夫电流定律，供热管网的支路流量等于该支路所有热力入口流量之和，写成矩阵的形式为：G＝AQG：支路流量；A：系数矩阵；Q：热力入口流量；基于基尔霍夫电压定律，在供热管网中应满足闭合回路压降的代数和为0，在供热管网中，经过总供回水管的管路压降方程的矩阵形式为：P-S(G*G)-(R+R
v
)(Q*Q)＝0P：总供回水压差；S：沿程阻力系数；R：热力入口阻力系数；R
v
：电动阀阻力系数；*：矩阵对应元素乘积运算；3)以水泵输送能耗最小为目标，以电动阀阻力、热力入口流量、总供回水压差为约束变量，建立非线性规划求解模型如下：上式中，目标函数除了满足水力平衡方程外，电动阀阻力、热力入口流量和总供回水压差分别满足最小阀阻R
v0
、最小流量Q0和最小压差0约束；由于线性规划求解模型中含有流量约束变量，对于未安装电动阀进行控制的热力入口，线性规划求解模型同样能够进行求解；4)优化求得水泵的频率及电动阀的开度，并通过监控平台对各水泵和电动阀门发送指令执行。2.根据权利要求1所述一种实现集中供热管网系统动态水力优化的方法，其特征在于，步骤(3)中非线性规划求解模型通过迭代法进行求解，迭代法的基本求解方法为：首先给定一个初始点x0∈R
n
，安照迭代规则生成一个点列{x
k
}，当{x
k
}是有限点列时，最后一个点便是这个非线性规划的最优解；设x
k
∈R
n
是第k轮迭代点，x
k+1
∈R
n
是第k+1轮迭代点，则x
k
和x
k+1
的关系为：x
k+1
＝x
k
+t
k
p
...

【专利技术属性】
技术研发人员：史鑫钰，杨俊红，刘德朝，崔棉善，马睿杰，
申请(专利权)人：格物智控天津能源科技有限公司，
类型：发明
国别省市：