CN112084707A - 基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 - Google Patents
基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 Download PDFInfo
- Publication number
- CN112084707A CN112084707A CN202010912061.0A CN202010912061A CN112084707A CN 112084707 A CN112084707 A CN 112084707A CN 202010912061 A CN202010912061 A CN 202010912061A CN 112084707 A CN112084707 A CN 112084707A
- Authority
- CN
- China
- Prior art keywords
- water
- cooling water
- model
- chilling unit
- machine room
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 223
- 239000000498 cooling water Substances 0.000 title claims abstract description 122
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 title claims abstract description 22
- 238000005057 refrigeration Methods 0.000 title claims description 40
- 238000005265 energy consumption Methods 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 30
- 238000010998 test method Methods 0.000 claims abstract description 6
- 238000013178 mathematical model Methods 0.000 claims description 18
- 238000012549 training Methods 0.000 claims description 18
- 239000013598 vector Substances 0.000 claims description 13
- 238000013528 artificial neural network Methods 0.000 claims description 12
- 238000001816 cooling Methods 0.000 claims description 11
- 238000010219 correlation analysis Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000011217 control strategy Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000009499 grossing Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000007710 freezing Methods 0.000 claims description 3
- 230000008014 freezing Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 14
- 230000009466 transformation Effects 0.000 abstract description 4
- 238000000611 regression analysis Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 9
- 238000004378 air conditioning Methods 0.000 description 8
- 238000004134 energy conservation Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000004615 ingredient Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000017525 heat dissipation Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000009467 reduction Effects 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Biophysics (AREA)
- Computing Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Air Conditioning Control Device (AREA)
Abstract
本发明公开了一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统,该方法在原冷水机组MP模型的基础上应用基于GRNN的冷却水进水温度建模方法对模型进行改进。通过设备建模、系统仿真、正交试验、回归分析,当建筑冷负荷确定时,使用了一种适用于制冷机房在冷冻水和冷却水变流量解耦的正交试验方法,来确定该条件下使制冷机房总能耗最小的冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度参数值,对优化前后的各设备及系统的能耗进行对比,并对各设备进行相应的调节,从而为冷水机组及制冷机房的节能改造提供理论依据和技术指导。
Description
技术领域
本发明涉及一种制冷机房节能优化方法,具体涉及一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统。
背景技术
大型公共建筑节能的最主要任务是空调系统的节能,实现空调系统节能必须解决其在部分负荷运行情况下能效偏低的问题。而变流量技术的应用可在一定程度上保证制冷系统高效地运行,因此,在空调水系统变流量研究基础上,对制冷系统进行有效的能源管理,对于公共建筑节能降耗具有重要意义。
目前随着对制冷机房节能技术的逐步深化研究,多种新型节能技术及调节措施被应用到制冷机房节能改造当中,但是制冷机房节能优化及调控仍存在部分设备不可调、不稳定、调控策略不科学、节能效果不理想等问题。因此如何在部分负荷运行情况下,确定合理的制冷机房系统设备运行参数,使系统在建筑运行期间都保持在一个相对节能的运行状态,是现在亟待解决的问题。
发明内容
针对现有制冷机房节能控制存在的不足,本发明提出一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统,该方法能够实现对各类模拟设备的有效控制,简化了制冷机房的控制方法,易于实现。
为了实现上述任务,本发明采用以下技术方案:
一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法,包括以下步骤:
首先采集冷源系统各设备运行参数,对冷水机组性能系数与影响因素进行相关性分析并选取冷水机组的基本数学模型;并利用广义回归神经网络的非线性映射特性,建立冷却水进水温度模型以此对冷水机组基本数学模型进行改进;
其次,在建立冷水机组能效模型和水泵功率模型的基础上,将冷冻水和冷却水变流量解耦控制策略植入平台,在平台上进行正交回归试验,采用多项式回归的方法得到制冷机房能耗与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
最后对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
作为本发明的进一步改进,所述冷水机组基本数学模型采用以下方法计算:
COP=D0+D1Qe+D2Tei+D3Tci+D4Qe 2+D5Tei 2+D6Tci 2+D7QeTei+D8QeTci+D9TeiTci;
式中,Qe为冷水机组实时制冷量;Tei为冷冻水回水温度;Tci为冷却水进水温度;D0~D9为与设备物理特性无关的回归系数。
作为本发明的进一步改进,建立冷却水进水温度模型的具体步骤为:
步骤3.1、数据预处理:
将所有数据进行归一化处理,具体计算方法如下:
式中:xi为输入数据;yi为输出数据;xi *为归一化后的输入向量;yi *为归一化后的输出向量;l为输入参数的个数;
步骤3.2、建立GRNN模型
将经过归一化处理后的Tco,Twb,ffan,Mc的学习样本作为网络输入层的输入向量,将Tci作为输出层的输出向量,训练样本数量用N表示;并将输入样本中包含的训练样本和测试样本的比例定为3:1;根据广义回归神经网络理论,构建的4层GRNN冷却水进水温度模型;
步骤3.3、确定网络参数
设定光滑因子σ值,从0.1开始,每次以增量0.01在[0.1,1]内递增;在同样训练样本和测试样本情况下,得到各光滑因子对应误差序列的均方根误差RMSE作为网络性能的评价指标,如式(3)所示,取得最小误差的σ值就是最优的参数;
步骤3.4、模型评价指标的确定
采用均方根误差RMSE和决定系数R2两个统计量对冷却水进水温度模型的预测精度进行评价如下:
步骤3.5、结果与分析
基于GRNN的冷却水进水温度参数模型为式(5),将式(5)的结果带入MP模型以此对模型进行改进,改进后的冷水机组能效数学模型如式(6)所示:
Tci=f(Tco,Twb,Mc,ffan) (5)
式中,Qe为冷水机组实时制冷量;Tei为冷冻水回水温度;Tco为冷却水出水温度;Twb为室外湿球温度;Mc为冷却水流量;ffan为冷却塔风机频率;D0~D9为与设备物理特性无关的回归系数;
作为本发明的进一步改进,采用多项式回归的方法得到水泵的功率随流量的关系曲线,具体如下:
步骤4.1、表达式的选取
对于变频水泵,满负荷工况下其功率计算公式如下:
η=ηp·ηm·ηVFD (8)
式中,Ppump为水泵功率;ρ为流体的密度;g为重力加速度;M为流体的体积流量;H为水泵的扬程;η为水泵的全效率;ηp为水泵效率;ηm为电机效率;ηVFD为变频器效率;
变频水泵大部分时间运行在部分负荷下,其功率模型表示为与流量的多项式关系,表达式如下:
其中Mrated为水泵的额定流量;a0~a3为模型系数;
步骤4.2、参数的辨识;
利用夏季和过渡季典型日采集到的冷冻水泵和冷却水泵的运行数据对水泵的功率模型进行辨识,得到水泵的功率随流量的关系曲线。
作为本发明的进一步改进,所述最佳运行参数是在建筑冷负荷确定的工况下,借助平台,采用正交回归试验法对冷水机组和水泵的功率进行计算,从而得到系统总能耗最小情况下的最佳参数设定值。
作为本发明的进一步改进,所述节能潜力对应的的总功率与冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度的拟合关系式如下:
式中,Me为冷冻水流量;Mc为冷却水流量;Teo为冷冻水供水温度;Tci为冷却水进水温度;
一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化系统,包括:
建立模型模块,用于采集冷源系统各设备运行参数,对冷水机组性能系数与影响因素进行相关性分析并选取冷水机组的基本数学模型;并利用广义回归神经网络的非线性映射特性,建立冷却水进水温度模型以此对冷水机组基本数学模型进行改进;
参数确定模块,用于在建立冷水机组能效模型和水泵功率模型的基础上,将冷冻水和冷却水变流量解耦控制策略在平台上进行正交回归试验,采用多项式回归的方法得到制冷机房能耗与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
节能潜力确定模块,用于对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
本发明与现有技术相比具有以下技术效果:
本发明基于GRNN的冷却水进水温度参数模型以对基本冷水机组MP模型进行改进,在建模研究上引入冷却水流量因素,可凸显变流量研究条件下流量因素在设备建模中的重要性。而引入室外湿球温度的变量因素,将室外气象条件的变化与冷源系统内部参数紧密结合起来,模型更直观,也更加具有现实意义。通过设备建模、系统仿真、正交试验、回归分析,当建筑冷负荷确定时,使用正交试验方法来确定该条件下使制冷机房总能耗最小的冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度参数值,对优化后的各设备及系统的能耗进行对比,并对各设备进行相应的调节,从而为冷水机组及制冷机房的节能改造提供理论依据和技术指导。
附图说明
图1为本发明的GRNN网络拓扑结构;
图2为本发明的不同平滑因子下GRNN的均方根误差结果;
图3为本发明的训练样本冷却水进水温度预测值与实际值对比结果;
图4为本发明的测试样本冷却水进水温度预测值与实际值对比结果;
图5为本发明的冷冻水泵功率随流量的关系曲线图;
图6为本发明的冷却水泵功率随流量的关系曲线图;
图7为本发明的冷源系统动态仿真流程图;
图8为本发明的中央空调冷源系统Trnsys仿真平台主界面。
图9为测试日优化前后能耗对比图。
具体实施方式
以下给出本发明的具体实施例,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换都属于本发明的保护范围。
针对研究对象实际空调冷源系统控制过程复杂和繁琐,且加入控制策略之后难以实现对相关设备的有效控制等问题,本发明提供的一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法,该方法基于Trnsys仿真软件高效灵活,部件之间的联系容易建立等特点,在该仿真平台上进行正交回归试验,通过数据处理得到使制冷机房系统功率最小的最佳运行参数,并将得到的系统能耗值与传统变流量控制方式下的能耗值进行比较,可得出在冷冻水和冷却水变流量解耦的研究背景下系统的节能潜力。
首先对本发明涉及的定义或概念内涵做以说明:
冷水机组性能系数(Coefficient Of Performance,COP):指在规定工况下,机组以同一单位表示的制冷(热)量除以总输入电功率的比值。
冷负荷:指为保持建筑物的热湿环境和所要求的室内温度,必须由空调系统从房间带走的热量叫空调房间冷负荷。
部分负荷率(Part Load Ratio,PLR):冷水机组的实际制冷量与额定制冷量的比值。
广义回归神经网络(General Regression Neural Network,GRNN):广义回归神经网络是一种基于one-pass学习算法的高度并行径向基网络。根据广义回归神经网络理论,GRNN总体的结构由四层构成,这四层分别是输入层,模式层,求和层和输出层。网络具有较强的非线性映射能力和柔性网络结构以及高度的容错性和鲁棒性,适用于解决非线性问题。
本发明的方法具体步骤为:
首先对冷水机组COP与影响因素进行相关性分析,选择冷水机组能效的基本数学模型。利用GRNN强大的非线性映射特性,以建立基于GRNN的冷却水进水温度参数模型对基本MP模型进行改进并建立水泵模型。
其次,在Trnsys冷源系统仿真平台上进行四因素二次回归正交试验,得到制冷机房与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
最后对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
以下结合具体实施例对本发明的进行详细说明。
实施例1
在本实施例中公开了一种基于GRNN的冷却水进水温度建模方法,为改进冷水机组能效模型提供可靠理论依据。
由于冷水机组COP受部分负荷率、冷冻水供回水温度、冷却水进出水温度、冷冻水流量和冷却水流量等诸多因素影响,且各种影响因素之间也存在着错综复杂的交互关系。
为了明确各运行参数对冷水机组运行能效的影响程度,本发明利用SPSS软件对冷水机组运行能效与多个影响因素进行相关性分析,通过Pearson相关系数、双侧显著性检验进行相关性分析,结果见表1。结果显示影响冷水机组COP的相关性大小依次为:部分负荷率PLR、冷却水进水温度Tci、冷冻水流量Me、冷冻水回水温度Tei、冷却水流量Mc、冷却水出水温度Tco、冷冻水供水温度Teo。通过相关性分析明确了哪些优化变量对降低冷水机组功耗贡献最大,从而可帮助工程师确定所选优化变量的优先级。
表1机组COP与影响因素相关性分析
影响因素 | 相关系数 |
部分负荷率 | 0.869 |
冷冻水供水温度 | 0.212 |
冷冻水回水温度 | 0.563 |
冷却水进水温度 | -0.693 |
冷却水出水温度 | -0.242 |
冷冻水流量 | 0.655 |
冷却水流量 | 0.262 |
运用测试数据对常见的基于经验的冷水机组数学模型进行参数辨识,采用CV值比较各模型的预测精度并考虑相关模型的训练样本数量及计算量,最终选用MP模型作为冷水机组的基本数学模型,表达式为:COP=D0+D1Qe+D2Tei+D3Tci+D4Qe 2+D5Tei 2+D6Tci 2+D7QeTei+D8QeTci+D9TeiTci。
式中,Qe为冷水机组实时制冷量,kW;Tei为冷冻水回水温度,℃;Tci为冷却水进水温度,℃;D0~D9为与设备物理特性无关的回归系数。
在本实施例中,以7月23日及10月15日分别作为夏季和过渡季的典型日,采集冷却水出水温度Tco、室外湿球温度Twb、室外干球温度T、冷却塔风机频率ffan和冷却水流量Mc的实时运行数据,采用SPSS统计软件对以上因素对冷却水进水温度Tci的影响程度进行分析,结果如下表所示。由相关性分析结果可知,冷却水出水温度Tco这一因素的相关系数最大,它是影响冷却水进水温度Tci的关键因素;其次是室外湿球温度Twb;室外干球温度T的相关系数最小,表明其数值的变化对冷却水进水温度Tci影响程度最小。因此本实施例选取冷却水出水温度Tco、室外湿球温度Twb、冷却塔风机频率ffan和冷却水流量Mc以上4种因素作为GRNN的输入向量建立冷却水进水温度模型。
表2冷却水进水温度与影响因素相关性分析
方法按照以下步骤执行:
步骤1、数据预处理。由于GRNN模型的输入输出数据具有不同的物理意义,为了消除量纲对神经网络的影响,需要将所有数据进行归一化处理,在本实施例中,采用离差标准化的方法将样本数据处理成[0,1]的无因次数,有效提高预测收敛速度并提升模型的精度,具体计算方法如下:
式中:xi为输入数据;yi为输出数据;xi *为归一化后的输入向量;yi *为归一化后的输出向量;l为输入参数的个数。
步骤2、建立GRNN模型。基于GRNN冷却水进水温度模型方法的主要思想是将经过归一化处理后的Tco,Twb,ffan,Mc的学习样本作为网络输入层的输入向量,将Tci作为输出层的输出向量,训练样本数量用N表示,如果输入样本中训练集和测试集的数量级差距过大,则很有可能会出现较大的预测误差,因此本发明将输入样本中包含的训练样本和测试样本的比例定为3:1。根据广义回归神经网络理论,构建的4层GRNN冷却水进水温度模型,如图1所示。
步骤3、确定网络参数。平滑因子σ的确定是GRNN网络训练的实质和难点,本发明设定光滑因子σ值,从0.1开始,每次以增量0.01在[0.1,1]内递增;在同样训练样本和测试样本情况下,得到各光滑因子对应误差序列的均方根误差RMSE作为网络性能的评价指标,见式(3)所示,取得最小误差的σ值就是最优的参数。
经过上述步骤计算,为了得到较高的模型泛化能力,最终确定最优光滑因子σ的数值为0.19,对应的均方根误差为0.6015,如图2所示。因此,将最优光滑因子σ=0.19应用于网络的输入层。
步骤4、模型评价指标的确定。在本实施例中,采用均方根误差RMSE和决定系数R2两个统计量对冷却水进水温度模型的预测精度进行评价,见式(4)所示。
步骤5、结果与分析。GRNN模型对训练样本的预测值与实际值对比见图3。总体来说,训练样本中绝大部分冷却水进水温度值均匀分布在实际值曲线两侧,少部分预测值出现不同程度的偏离。对于36组训练样本,模型预测值与实际值的拟合程度较高,决定系数R2为0.9797,RMSE为0.1114℃,可见GRNN训练网络的模拟精度符合要求。
用12组测试样本验证已建立的GRNN模型的泛化性能与模拟精度,其预测值与实际值对比如图4所示。模型对测试样本的预测值与实际值的拟合程度较高,决定系数R2为0.9536。预测值与实际值的RMSE较训练样本略有增加,为0.1388℃。分析结果可知GRNN模型对测试样本的预测精度较高,模型预测值与实际值的误差在10%的允许范围内,可用于冷却水进水温度模型预测。
由分析可知基于GRNN的冷却水进水温度参数模型可表达成式(5),因此本发明提出将式(5)的结果带入MP模型以此对模型进行改进,改进后的冷水机组能效数学模型如式(6)所示。
Tci=f(Tco,Twb,Mc,ffan) (5)
改进的冷水机组能效模型在建模研究上引入冷却水流量因素,可凸显变流量研究条件下流量因素在设备建模中的重要性。而引入室外湿球温度的变量因素,将室外气象条件的变化与冷源系统内部参数紧密结合起来,模型更直观,也更加具有现实意义。
实施例2
在本实施例中采用多项式回归的方法得到水泵的功率随流量的关系曲线,具体如下:
步骤1、表达式的选取。对于变频水泵,满负荷工况下其功率计算公式如下:
η=ηp·ηm·ηVFD (8)
式中,Ppump为水泵功率,kW;ρ为流体的密度,kg/m3;g为重力加速度,m/s2;M为流体的体积流量,m3/h;H为水泵的扬程,m;η为水泵的全效率;ηp为水泵效率;ηm为电机效率;ηVFD为变频器效率。
从实际分析看,变频水泵大部分时间运行在部分负荷下,其功率模型可表示为与流量的多项式关系,表达式如下:
其中Mrated为水泵的额定流量,m3/h;a0~a3为模型系数。
步骤2、参数的辨识。在本实施例中,利用夏季和过渡季典型日采集到的冷冻水泵和冷却水泵的运行数据对水泵的功率模型进行辨识,结果如图5和图6所示。
实施例3
本实施例为冷冻水和冷却水变流量解耦控制策略在Trnsys仿真平台上的应用。结合本发明的综合办公楼研究对象,空调冷源系统各环路负荷特性和压损相差较大,每层的房间内部需要有个性化设置等特性,选择定供回水干管压差控制作为冷冻水变流量控制方式能满足上述要求。在冷却水变流量侧,由于空调冷源系统的运行大部分处于部分负荷情况下,采用定温差控制方式相当于冷却水流量随负荷成比例变化,控制方法简单、有效,因此本实施例在冷却水变流量侧采用定温差控制方法。
具体步骤如下:
步骤1、基于工程案例,该冷源系统动态仿真平台主要由气象模块、冷水机组模块、水泵模块、冷却塔模块、混流器模块、PID控制模块和数据输出模块等组成,系统动态仿真流程如图7所示。
步骤2、仿真平台的搭建及控制策略的应用。仿真平台的主界面如图8所示。
实施例4
本发明借助建立的设备模型和Trnsys仿真平台,在该仿真平台上进行正交试验,通过数据处理得到使制冷机房系统功率最小的最佳运行参数,将得到的最小功率值与传统变流量控制方式下的功率值进行比较,得到在冷冻水和冷却水变流量解耦背景下系统的节能潜力。
对于确定设备型号的制冷机房,冷冻水流量、冷冻水供水温度和冷却水流量是可以通过设备控制的变量。而冷却水进水温度受冷却塔散热性能和冷却水流量影响,但冷却塔散热性能主要与冷却水流量和进塔风量有关,冷却水流量可以通过改变冷却水泵频率进行调节,进塔风量又可以通过改变风机频率调节,因此冷却水进水温度间接上也是可控的。综上所述,选择冷冻水供水温度Teo、冷冻水流量Me、冷却水进水温度Tci、冷却水流量Mc作为系统的优化变量;
为了量化制冷机房在冷冻水和冷却水变流量控制策略下的运行能耗,在建筑冷负荷确定的工况下,借助Trnsys仿真平台,采用正交回归试验法对冷水机组和水泵的功率进行计算,从而得到系统总能耗最小情况下的最佳参数设定值。选取的过渡季某日9时的冷负荷为987.35kW,考虑到冷冻水和冷却水流量均存在一个最小值,超过这个最小值将导致冷水机组效率大幅度下降,因此,在试验中保证水流量在50%~100%之间调节;冷冻水供水温度范围设定为4~12℃;冷却水进水温度范围设定为15~35℃。为了简化分析,选择只开启一台冷水机组、一台冷冻水泵和一台冷却水泵的工况作为本实施例的研究背景,不考虑空调末端设备及其能耗。具体各参数取值见表3。试验结果见表4。
表3正交试验水平编码
表4四因素二次回归正交试验结果
对上述试验结果进行数据处理,剔除不显著因子后得到制冷机房系统总功率与冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度的拟合关系式,如式(10)所示:
对上述回归模型进行F检验,式(11)可以看出回归模型极显著。
F=1110.13>F0.01(9.5)=14.7 (11)
根据回归模型进行分析,可得该模型的最小值为189.90kW,此时的冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度分别为164.83m3/h、290.97m3/h、8.14℃、15.00℃。将参数代入设备模型,计算所得的总功率值为192.71kW,与分析值非常接近,误差仅为1.57%,这说明通过正交试验寻找制冷机房总功率最小值具有一定的可靠性。图9所示为在优化前后水泵、冷水机组和制冷机房系统能耗对比。水泵的能耗增加了3.44kW,但是冷水机组能耗降低了29.88kW,制冷机房系统总的节能率达到了12.23%。从优化结果可以看出,冷水机组能耗是制冷机房系统节能的重点,系统参数优化后的结果是水泵能耗略有增加,冷水机组能耗大幅降低。
本发明第二个目的是提供一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化系统,包括:
建立模型模块,用于采集冷源系统各设备运行参数,对冷水机组性能系数与影响因素进行相关性分析并选取冷水机组的基本数学模型;并利用广义回归神经网络的非线性映射特性,建立冷却水进水温度模型以此对冷水机组基本数学模型进行改进;
参数确定模块,用于在建立冷水机组能效模型和水泵功率模型的基础上,将冷冻水和冷却水变流量解耦控制策略在平仿真平台上进行正交回归试验,采用多项式回归的方法得到制冷机房能耗与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
节能潜力确定模块,用于对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
总之,本发明一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统,在原冷水机组MP模型的基础上应用基于GRNN的冷却水进水温度建模方法对模型进行改进。通过设备建模、系统仿真、正交试验、回归分析,当建筑冷负荷确定时,使用了一种适用于制冷机房在冷冻水和冷却水变流量解耦的正交试验方法,来确定该条件下使制冷机房总能耗最小的冷冻水流量、冷却水流量、冷冻水供水温度和冷却水进水温度参数值,对优化后的各设备及系统的能耗进行对比,并对各设备进行相应的调节,从而为冷水机组及制冷机房的节能改造提供理论依据和技术指导。
以上披露的所有文章和参考资料,包括专利申请和出版物,出于各种目的通过援引结合于此。描述组合的术语“基本由…构成”应该包括所确定的元件、成分、部件或步骤以及实质上没有影响该组合的基本新颖特征的其他元件、成分、部件或步骤。使用术语“包含”或“包括”来描述这里的元件、成分、部件或步骤的组合也想到了基本由这些元件、成分、部件或步骤构成的实施方式。这里通过使用术语“可以”,旨在说明“可以”包括的所描述的任何属性都是可选的。
多个元件、成分、部件或步骤能够由单个集成元件、成分、部件或步骤来提供。另选地,单个集成元件、成分、部件或步骤可以被分成分离的多个元件、成分、部件或步骤。用来描述元件、成分、部件或步骤的公开“一”或“一个”并不说为了排除其他的元件、成分、部件或步骤。
应该理解,以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述,在所提供的示例之外的许多实施例和许多应用对本领域技术人员来说都将是显而易见的。因此,本教导的范围不应该参照上述描述来确定,而是应该参照前述权利要求以及这些权利要求所拥有的等价物的全部范围来确定。出于全面之目的,所有文章和参考包括专利申请和公告的公开都通过参考结合在本文中。在前述权利要求中省略这里公开的主题的任何方面并不是为了放弃该主体内容,也不应该认为申请人没有将该主题考虑为所公开的发明主题的一部分。
Claims (7)
1.一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法,其特征在于,包括以下步骤:
首先采集冷源系统各设备运行参数,对冷水机组性能系数与影响因素进行相关性分析并选取冷水机组的基本数学模型;并利用广义回归神经网络的非线性映射特性,建立冷却水进水温度模型以此对冷水机组基本数学模型进行改进;
其次,在建立冷水机组能效模型和水泵功率模型的基础上,将冷冻水和冷却水变流量解耦控制策略植入平台,在平台上进行正交回归试验,采用多项式回归的方法得到制冷机房能耗与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
最后对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
2.根据权利要求1所述的方法,其特征在于,所述冷水机组基本数学模型采用以下方法计算:
COP=D0+D1Qe+D2Tei+D3Tci+D4Qe 2+D5Tei 2+D6Tci 2+D7QeTei+D8QeTci+D9TeiTci;
式中,Qe为冷水机组实时制冷量;Tei为冷冻水回水温度;Tci为冷却水进水温度;D0~D9为与设备物理特性无关的回归系数。
3.根据权利要求1所述的方法,其特征在于,建立冷却水进水温度模型的具体步骤为:
步骤3.1、数据预处理:
将所有数据进行归一化处理,具体计算方法如下:
式中:xi为输入数据;yi为输出数据;xi *为归一化后的输入向量;yi *为归一化后的输出向量;l为输入参数的个数;
步骤3.2、建立GRNN模型
将经过归一化处理后的Tco,Twb,ffan,Mc的学习样本作为网络输入层的输入向量,将Tci作为输出层的输出向量,训练样本数量用N表示;并将输入样本中包含的训练样本和测试样本的比例定为3:1;根据广义回归神经网络理论,构建的4层GRNN冷却水进水温度模型;
步骤3.3、确定网络参数
设定光滑因子σ值,从0.1开始,每次以增量0.01在[0.1,1]内递增;在同样训练样本和测试样本情况下,得到各光滑因子对应误差序列的均方根误差RMSE作为网络性能的评价指标,如式(3)所示,取得最小误差的σ值就是最优的参数;
步骤3.4、模型评价指标的确定
采用均方根误差RMSE和决定系数R2两个统计量对冷却水进水温度模型的预测精度进行评价如下:
步骤3.5、结果与分析
基于GRNN的冷却水进水温度参数模型为式(5),将式(5)的结果带入MP模型以此对模型进行改进,改进后的冷水机组能效数学模型如式(6)所示:
Tci=f(Tco,Twb,Mc,ffan) (5)
式中,Qe为冷水机组实时制冷量;Tei为冷冻水回水温度;Tco为冷却水出水温度;Twb为室外湿球温度;Mc为冷却水流量;ffan为冷却塔风机频率;D0~D9为与设备物理特性无关的回归系数。
4.根据权利要求1所述的方法,其特征在于,采用多项式回归的方法得到水泵的功率随流量的关系曲线,具体如下:
步骤4.1、表达式的选取
对于变频水泵,满负荷工况下其功率计算公式如下:
η=ηp·ηm·ηVFD (8)
式中,Ppump为水泵功率;ρ为流体的密度;g为重力加速度;M为流体的体积流量;H为水泵的扬程;η为水泵的全效率;ηp为水泵效率;ηm为电机效率;ηVFD为变频器效率;
变频水泵大部分时间运行在部分负荷下,其功率模型表示为与流量的多项式关系,表达式如下:
其中Mrated为水泵的额定流量;a0~a3为模型系数;
步骤4.2、参数的辨识;
利用夏季和过渡季典型日采集到的冷冻水泵和冷却水泵的运行数据对水泵的功率模型进行辨识,得到水泵的功率随流量的关系曲线。
5.根据权利要求1所述的方法,其特征在于,所述最佳运行参数是在建筑冷负荷确定的工况下,借助平台,采用正交回归试验法对冷水机组和水泵的功率进行计算,从而得到系统总能耗最小情况下的最佳参数设定值。
7.一种基于冷冻水和冷却水变流量解耦的制冷机房节能优化系统,其特征在于,包括:
建立模型模块,用于采集冷源系统各设备运行参数,对冷水机组性能系数与影响因素进行相关性分析并选取冷水机组的基本数学模型;并利用广义回归神经网络的非线性映射特性,建立冷却水进水温度模型以此对冷水机组基本数学模型进行改进;
参数确定模块,用于在建立冷水机组能效模型和水泵功率模型的基础上,将冷冻水和冷却水变流量解耦控制策略在平台上进行正交回归试验,采用多项式回归的方法得到制冷机房能耗与各参数之间的经验关系式,通过数据处理得到使制冷机房系统功率最小的最佳运行参数;
节能潜力确定模块,用于对比优化前后水泵、冷水机组和制冷机房系统能耗,得到冷冻水和冷却水变流量解耦条件下系统的节能潜力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010912061.0A CN112084707B (zh) | 2020-09-02 | 2020-09-02 | 基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010912061.0A CN112084707B (zh) | 2020-09-02 | 2020-09-02 | 基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112084707A true CN112084707A (zh) | 2020-12-15 |
CN112084707B CN112084707B (zh) | 2024-08-20 |
Family
ID=73733108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010912061.0A Active CN112084707B (zh) | 2020-09-02 | 2020-09-02 | 基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112084707B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112665102A (zh) * | 2020-12-25 | 2021-04-16 | 江苏联宏智慧能源股份有限公司 | 一种基于冷量计算的风机盘管控制方法 |
CN113269365A (zh) * | 2021-06-02 | 2021-08-17 | 西安建筑科技大学 | 一种基于麻雀优化算法的短期空调负荷预测方法及系统 |
CN113739368A (zh) * | 2021-08-31 | 2021-12-03 | 广州汇电云联互联网科技有限公司 | 一种中央空调系统的冷站控制方法及系统 |
CN113790516A (zh) * | 2021-09-18 | 2021-12-14 | 深圳达实智能股份有限公司 | 中央空调制冷站全局优化节能控制方法、系统和电子设备 |
CN113806938A (zh) * | 2021-09-17 | 2021-12-17 | 上海美控智慧建筑有限公司 | 制冷机房能效计算方法、装置和电子设备 |
CN113821902A (zh) * | 2021-06-17 | 2021-12-21 | 李壮举 | 一种用于中央空调制冷站静态优化的自抗扰控制系统 |
CN114396731A (zh) * | 2021-12-30 | 2022-04-26 | 青岛海信日立空调系统有限公司 | 制冷机组 |
CN114754432A (zh) * | 2022-04-26 | 2022-07-15 | 西安科技大学 | 一种基于降温需冷量变化的冷冻水供给调控方法及系统 |
CN115355580A (zh) * | 2022-08-16 | 2022-11-18 | 浙江大冲能源科技股份有限公司 | 中央空调二次泵平衡管蓄冷系统 |
CN115796055A (zh) * | 2023-01-10 | 2023-03-14 | 北京云庐科技有限公司 | 一种基于完整空调系统仿真模型的优化运行调节方法 |
WO2024146529A1 (zh) * | 2023-01-04 | 2024-07-11 | 青岛海信日立空调系统有限公司 | 制冷机房及其能效确定方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130167560A1 (en) * | 2010-10-13 | 2013-07-04 | Weldtech Technology (Shanghai) Co., Ltd. | Energy-saving optimized control system and method for refrigeration plant room |
WO2015172560A1 (zh) * | 2014-05-16 | 2015-11-19 | 华南理工大学 | 基于bp神经网络的中央空调冷负荷的预测方法 |
WO2018004464A1 (en) * | 2016-06-29 | 2018-01-04 | Kirkham Group Pte Ltd | Large scale machine learning-based chiller plants modeling, optimization and diagnosis |
CN109855238A (zh) * | 2019-02-27 | 2019-06-07 | 四川泰立智汇科技有限公司 | 一种中央空调建模与能效优化方法及装置 |
CN110107989A (zh) * | 2019-04-30 | 2019-08-09 | 北京工业大学 | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 |
CN110288164A (zh) * | 2019-07-02 | 2019-09-27 | 广州市特沃能源管理有限公司 | 一种建筑空调制冷站系统预测控制方法 |
CN110688740A (zh) * | 2019-09-10 | 2020-01-14 | 天津大学 | 基于Modelica联合仿真优化的冷水机房模型校准方法 |
CN111256294A (zh) * | 2020-01-17 | 2020-06-09 | 深圳市得益节能科技股份有限公司 | 一种基于模型预测的冷水机组组合运行优化控制方法 |
-
2020
- 2020-09-02 CN CN202010912061.0A patent/CN112084707B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130167560A1 (en) * | 2010-10-13 | 2013-07-04 | Weldtech Technology (Shanghai) Co., Ltd. | Energy-saving optimized control system and method for refrigeration plant room |
WO2015172560A1 (zh) * | 2014-05-16 | 2015-11-19 | 华南理工大学 | 基于bp神经网络的中央空调冷负荷的预测方法 |
WO2018004464A1 (en) * | 2016-06-29 | 2018-01-04 | Kirkham Group Pte Ltd | Large scale machine learning-based chiller plants modeling, optimization and diagnosis |
CN109855238A (zh) * | 2019-02-27 | 2019-06-07 | 四川泰立智汇科技有限公司 | 一种中央空调建模与能效优化方法及装置 |
CN110107989A (zh) * | 2019-04-30 | 2019-08-09 | 北京工业大学 | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 |
CN110288164A (zh) * | 2019-07-02 | 2019-09-27 | 广州市特沃能源管理有限公司 | 一种建筑空调制冷站系统预测控制方法 |
CN110688740A (zh) * | 2019-09-10 | 2020-01-14 | 天津大学 | 基于Modelica联合仿真优化的冷水机房模型校准方法 |
CN111256294A (zh) * | 2020-01-17 | 2020-06-09 | 深圳市得益节能科技股份有限公司 | 一种基于模型预测的冷水机组组合运行优化控制方法 |
Non-Patent Citations (4)
Title |
---|
刘国强;刘志斌;闫军威;周璇;: "集中式空调冷冻水系统运行特性分析与节能优化运行方法研究", 建筑科学, no. 04, 30 April 2018 (2018-04-30) * |
张广丽;陈丽萍;: "制冷机房的运行能耗优化", 流体机械, no. 08, 31 August 2012 (2012-08-31) * |
蒋小强;龙惟定;李敏;: "制冷机房冷冻水和冷却水变流量解耦控制策略", 重庆大学学报, no. 06, 30 June 2010 (2010-06-30) * |
陈霈;牛洪海;管晓晨;李兵;: "冷水机组运行特性建模与优化控制研究", 节能, no. 03, 31 March 2019 (2019-03-31) * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112665102A (zh) * | 2020-12-25 | 2021-04-16 | 江苏联宏智慧能源股份有限公司 | 一种基于冷量计算的风机盘管控制方法 |
CN113269365B (zh) * | 2021-06-02 | 2023-06-16 | 西安建筑科技大学 | 一种基于麻雀优化算法的短期空调负荷预测方法及系统 |
CN113269365A (zh) * | 2021-06-02 | 2021-08-17 | 西安建筑科技大学 | 一种基于麻雀优化算法的短期空调负荷预测方法及系统 |
CN113821902A (zh) * | 2021-06-17 | 2021-12-21 | 李壮举 | 一种用于中央空调制冷站静态优化的自抗扰控制系统 |
CN113739368A (zh) * | 2021-08-31 | 2021-12-03 | 广州汇电云联互联网科技有限公司 | 一种中央空调系统的冷站控制方法及系统 |
CN113806938A (zh) * | 2021-09-17 | 2021-12-17 | 上海美控智慧建筑有限公司 | 制冷机房能效计算方法、装置和电子设备 |
CN113790516A (zh) * | 2021-09-18 | 2021-12-14 | 深圳达实智能股份有限公司 | 中央空调制冷站全局优化节能控制方法、系统和电子设备 |
CN114396731A (zh) * | 2021-12-30 | 2022-04-26 | 青岛海信日立空调系统有限公司 | 制冷机组 |
CN114754432A (zh) * | 2022-04-26 | 2022-07-15 | 西安科技大学 | 一种基于降温需冷量变化的冷冻水供给调控方法及系统 |
CN115355580A (zh) * | 2022-08-16 | 2022-11-18 | 浙江大冲能源科技股份有限公司 | 中央空调二次泵平衡管蓄冷系统 |
WO2024146529A1 (zh) * | 2023-01-04 | 2024-07-11 | 青岛海信日立空调系统有限公司 | 制冷机房及其能效确定方法 |
CN115796055A (zh) * | 2023-01-10 | 2023-03-14 | 北京云庐科技有限公司 | 一种基于完整空调系统仿真模型的优化运行调节方法 |
CN115796055B (zh) * | 2023-01-10 | 2023-04-07 | 北京云庐科技有限公司 | 一种基于完整空调系统仿真模型的优化运行调节方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112084707B (zh) | 2024-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112084707B (zh) | 基于冷冻水和冷却水变流量解耦的制冷机房节能优化方法及系统 | |
CN110285532B (zh) | 基于人工智能的机房空调控制方法、装置及系统 | |
CN106920006B (zh) | 一种基于isoa-lssvm的地铁站空调系统能耗预测方法 | |
CN109855238B (zh) | 一种中央空调建模与能效优化方法及装置 | |
CN103912966B (zh) | 一种地源热泵制冷系统优化控制方法 | |
CN103049612B (zh) | 一种基于模型降阶技术的建筑室内环境优化方法 | |
CN110107989A (zh) | 基于冷冻水回水温度最佳设定点的小型定频冷水机组变水温控制方法 | |
CN114383299B (zh) | 基于大数据与动态仿真的中央空调系统运行策略优化方法 | |
CN107940679B (zh) | 一种基于数据中心冷水机组性能曲线的群控方法 | |
CN111723533B (zh) | 一种地源热泵系统变频水泵节能计算方法 | |
CN110567101A (zh) | 一种基于支持向量机模型的冷水机组高能效控制方法 | |
CN111649459B (zh) | 基于专家pid的纺织空调节能自控方法 | |
CN113268913B (zh) | 一种基于pso-elm算法的智能建筑空调冷机系统运行优化方法 | |
CN115776795A (zh) | 数据中心空调系统诊断方法及装置 | |
Liu et al. | Optimization based integrated control of building HVAC system | |
CN113339941A (zh) | 一种变频空调器控制方法 | |
CN115115143A (zh) | 基于ai算法的冷却塔最优开启台数和最小运行能耗的计算方法 | |
Zhuang et al. | A decentralized method for energy conservation of an HVAC system | |
CN111288610A (zh) | 一种变风量空调系统变静压自适应模糊控制方法 | |
CN110688740A (zh) | 基于Modelica联合仿真优化的冷水机房模型校准方法 | |
Ma et al. | Test and evaluation of energy saving potentials in a complex building central chilling system using genetic algorithm | |
CN117592657A (zh) | 一种多台冷水机组负荷分配优化方法及系统 | |
CN117366799A (zh) | 建筑物冷源系统的控制方法 | |
CN114543273B (zh) | 一种集中空调供冷系统自适应深度学习优化节能控制算法 | |
Yang et al. | Hybrid artificial neural network− genetic algorithm technique for condensing temperature control of air-cooled chillers |
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 |