CN103713613A - Propr模式下火电机组负荷优化控制方法 - Google Patents

Propr模式下火电机组负荷优化控制方法 Download PDF

Info

Publication number
CN103713613A
CN103713613A CN201410001773.1A CN201410001773A CN103713613A CN 103713613 A CN103713613 A CN 103713613A CN 201410001773 A CN201410001773 A CN 201410001773A CN 103713613 A CN103713613 A CN 103713613A
Authority
CN
China
Prior art keywords
load
unit
value
propr
changing rate
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
Application number
CN201410001773.1A
Other languages
English (en)
Other versions
CN103713613B (zh
Inventor
李军
周长来
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201410001773.1A priority Critical patent/CN103713613B/zh
Publication of CN103713613A publication Critical patent/CN103713613A/zh
Application granted granted Critical
Publication of CN103713613B publication Critical patent/CN103713613B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Control Of Eletrric Generators (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明公开了一种PROPR模式下火电机组负荷优化控制方法,具体步骤为:步骤一:对调度中心下发的AGC负荷指令值进行分析,通过函数对负荷指令值进行优化处理;步骤二:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,生成一汽机主控前馈量;步骤三:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,进行负荷变化率的切换。本发明通过对机组的负荷指令值和实际功率值进行分析处理,对机组的负荷实时变化情况进行判断,对机组的汽机主控进行前馈优化调整,并且对机组的负荷变化率进行折线优化处理,从而在保证机组在PROPR模式下调节速率指标的同时,有效提升调节精度和响应时间指标,确保了电网负荷的供需平衡,保证了电网频率的稳定。

Description

PROPR模式下火电机组负荷优化控制方法
技术领域
本发明涉及一种火电机组负荷优化控制方法,尤其是一种PROPR模式下火电机组负荷优化控制方法。
背景技术
随着我国电力系统的快速发展,电网容量越来越大,对电能的品质也要求越来越高,为了电网的安全稳定运行,各大型火电机组都要求投入AGC功能,尤其是要PROPR(等可调比例正常调节)模式,这种模式要求机组的负荷变化范围大(一般要求50%--100%额定负荷),在此范围内要求机组快速、准确、稳定的响应电网发出的负荷变化需求。
(1)AGC考核标准
如图1典型AGC机组设点控制过程所示,这是网内某台机组一次典型的AGC机组设点控制过程。图中,Pmin,i是该机组可调的下限出力,Pmax,i是其可调的上限出力,PNi是其额定出力,Pdi是其启停磨临界点功率。整个过程可以这样描述:T0时刻以前,T1时刻以前,该机组稳定运行在出力值P1附近,T0时刻,AGC控制程序对该机组下发功率为P2的设点命令,机组开始涨出力,到T1时刻可靠跨出P1的调节死区,然后到T2时刻进入启磨区间,一直到T3时刻,启磨过程结束,机组继续涨出力,至T4时刻第一次进入调节死区范围,然后在P2附近小幅振荡,并稳定运行于P2附近,直至T5时刻,AGC控制程序对该机组发出新的设点命令,功率值为P3,机组随后开始降出力的过程,T6时刻可靠跨出调节死区,至T7时刻进入P3的调节死区,并稳定运行于其附近。
AGC考核主要涉及调节速率、调节精度、响应时间三个主要指标:
1)调节速率K1:是指机组响应设点指令的速率,可分为上升速率和下降速率。
2)调节精度K2:是指机组响应稳定以后,实际出力和设点出力之间的差值。
3)响应时间K3:是指EMS(调度能量管理系统)发出指令之后,机组出力在原出力点的基础上,可靠地跨出与调节方向一致的调节死区所用的时间。
(2)火电机组常规负荷控制方案
图2中,T为切换器,进行支路切换;A为模拟量发生器,可设定模拟量数值;≮、≯为不小于和不大于判断器,分别接收来自模拟量发生器A的负荷下限和负荷上限的限定,不超过设定值则输出与输入相同;V≯为速率限制器,接收负荷变化率的限定,可控制输入量的变化速率。机组若要投入AGC,前提是机组工作在协调(CCS)方式下,此时的负荷指令LDSP是经过速率限制以及负荷下限、上限限制后的值,所述速率以及负荷上下限由电厂操作人员设定。当在CCS方式下投入AGC后,LDSP前的切换器切至AGC支路,LDSP的值即为中调指令;否则为电厂操作员设定的指令。
火电机组在投入协调控制方式(CCS)时,一般而言协调控制方式为汽轮机侧进行负荷控制;锅炉侧控制压力。而在目前电网调度模式下,电网考核的是机组的调节速率、调节精度和响应时间,因此火电机组对汽轮机侧的控制要求比较严格,控制的精度要求较高。常规火电机组协调方式下汽机主控逻辑如图3所示。其中,PID环节PV输入量为现场变送器采集来的机组实际功率值,SP输入量为一次调频前负荷指令值(LDC OUT)与一次调频产生的负荷增量的叠加值,PV与SP的差值经过PID运算最终生成汽机主控指令,进而控制汽轮机的调门开度,改变蒸汽流量,从而达到控制机组负荷的目的。有时为加快反映速度,可以引入负荷指令的比例微分前馈。
实际运行证明,这种控制方案可以通过人为增大机组的负荷变化率,来提高考核指标总的调节速率,但是负荷变化率的增大会造成机组的主汽压力波动较大,从而对机组负荷的控制精度造成负面影响,影响到调节精度和响应时间,最终会造成整体考核结果不理想,严重时运行控制结果在调度的考核评价中出现品质不达标情况。
发明内容
本发明的目的就是为解决上述问题,提供一种PROPR模式下火电机组负荷优化控制方法,它通过对机组的负荷实时变化情况进行判断,对机组的汽机主控进行前馈优化调整,并且对机组的负荷变化率进行折线优化处理,从而在保证机组在PROPR模式下调节速率指标的同时,有效提升调节精度和响应时间,进而确保机组的整体控制指标达到优秀。
为实现上述目的,本发明采用如下技术方案:
一种PROPR模式下火电机组负荷优化控制方法,具体步骤为:
步骤一:对调度中心下发的AGC负荷指令值进行求差分析,找出在指定负荷点下的偏差量,并将偏差量叠加到原指令值上,通过非线性修正函数对指令进行优化处理;
步骤二:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,求出偏差量,将偏差量送至一饱和非线性环节,生成一汽机主控前馈量;
步骤三:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,求出偏差量,根据偏差量是否超过设定限值,通过控制切换器模块的选择开关,进行负荷变化率的切换。
所述步骤一的具体步骤为:利用GPS同步对时系统,将同一时间点下的调度中心下发的AGC负荷指令值和CCS接收到的负荷指令值进行求差。
所述步骤二中,死区的设置原则为功率偏差值小幅震荡下不应产生汽机主控前馈量,一般设置为±0.1-0.3MW,饱和非线性环节的具体函数如下:
所述步骤三的具体步骤为:
将偏差量送至一高低限报警模块中,设置上下限,越限则输出“1”,反之为“0”;将高低限报警模块的输出作为一切换器T2的选择开关,修正后的负荷变化率值和原有的负荷变化率值同时送至切换器T2的两个切换输入端中,若为“1”,则将原有的负荷变化率值送出,为“0”,则切换器T2将修正后的负荷变化率值送出。
所述高低限报警模块的上下限H/L的设置原则为:按规定300MW机组在PROPR模式下单次负荷调节的幅值为±3MW,其计算调节速率的负荷取点在负荷实际变化至±2MW之前,因此将H和L分别设置为2和-2。
所述切换器T2的两个输入值设置原则为:一个为原有负荷变化率值,另一个为原有负荷变化率值的60%-70%,也就是说,当实际功率变化将至负荷指令值时,以较低的负荷变化率接近,进而提高调节精度指标。
本发明的有益效果:本发明通过对机组的负荷指令值和实际功率值进行分析处理,对机组的负荷实时变化情况进行判断,对机组的汽机主控进行前馈优化调整,并且对机组的负荷变化率进行折线优化处理,从而在保证机组在PROPR模式下调节速率指标的同时,有效提升调节精度和响应时间,切实有效的确保了火电机组控制品质和考核指标,进而确保了电网负荷的供需平衡,保证了电网频率的稳定。
附图说明
图1为本发明典型AGC机组设点控制过程;
图2为火电机组负荷产生逻辑图;
图3为火电机组协调方式下汽机主控逻辑图;
图4为优化后火电机组协调方式下汽机主控逻辑图;
图5为优化后火电机组负荷产生逻辑图;
图6为饱和非线性函数的折线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图4、图5所示,在某300MW亚临界、中间再热、凝汽式燃煤机组中加以应用。
(1)因为调度考核机组的负荷范围一般为50%-100%额定功率,因此只需对150-300MW之间的负荷点进行修正,利用GPS同步对时系统,取出同一时间点下调度中心下发的AGC负荷指令值和CCS接收到的负荷指令值送入减法器模块中,找出在指定负荷点下的偏差量,并将偏差量叠加到原指令值上,生成一非线性修正函数F(x),如表1所示。
表1
Figure BDA0000452306360000041
Figure BDA0000452306360000051
(2)将经过修正的AGC负荷指令值和机组实际功率值送至减法器模块中,对两者进行减法运算,求出偏差量,将偏差量送至一饱和非线性环节,分析运算后输出一数值作为汽机主控的前馈。其中,死区的设置原则为功率偏差值小幅震荡下不应产生汽机主控前馈量,一般可设置为±0.1-0.3MW。按规定300MW机组在PROPR模式下单次负荷调节的幅值为±3MW,此机组根据实际情况,设置具体函数如表2、图6所示,实际运行证明可有效提升响应时间指标K3。
表2饱和非线性函数
Figure BDA0000452306360000052
(3)将修正后的AGC负荷指令值和机组实际功率值送至减法器模块中,对两者进行减法运算,求出偏差量,将偏差量送至高低限报警模块中,越限则输出“1”,反之为“0”;将高低限报警模块的输出作为一切换器T2的选择开关,修正后的负荷变化率值和原有的负荷变化率值同时送至切换器T2的两个切换输入端中,若为“1”,则将原有的负荷变化率值送出,为“0”,则切换器T2将修正后的负荷变化率值送出。
其中,高低限报警模块的上下限(H/L)的设置原则为:按规定300MW机组在PROPR模式下单次负荷调节的幅值为±3MW,其计算调节速率的负荷取点在负荷实际变化至±2MW之前,因此将H和L分别设置为2和-2。
切换器T2两个输入值设置原则为:一个为原有负荷变化率值,另一个为原有负荷变化率值的60%-70%。也就是说,当实际功率变化将至负荷指令值时,以较低的负荷变化率接近,进而提高调节精度指标。
如表3所示,1-3行7月份考核指标数据为优化前,后3行8月份考核指标数据为优化后,通过对比可以看出,在保证调节速率的前提下,技改后机组的调节精度和响应时间都有了明显提高,综合性能指标KP(KP=K1*K2*K3)提高14%以上。
表3
Figure BDA0000452306360000061
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (6)

1.一种PROPR模式下火电机组负荷优化控制方法,其特征是,具体步骤为:
步骤一:对调度中心下发的AGC负荷指令值进行求差分析,找出在指定负荷点下的偏差量,并将偏差量叠加到原指令值上,通过非线性修正函数对指令进行优化处理;
步骤二:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,求出偏差量,将偏差量送至一饱和非线性环节,生成一汽机主控前馈量;
步骤三:对修正后的AGC负荷指令值和机组实际功率值进行求差、判断,求出偏差量,根据偏差量是否超过设定限值,通过控制切换器模块的选择开关,进行负荷变化率的切换。
2.如权利要求1所述一种PROPR模式下火电机组负荷优化控制方法,其特征是,所述步骤一的具体步骤为:利用GPS同步对时系统,将同一时间点下的调度中心下发的AGC负荷指令值和CCS接收到的负荷指令值进行求差。
3.如权利要求1所述一种PROPR模式下火电机组负荷优化控制方法,其特征是,所述步骤二中,死区的设置原则为功率偏差值小幅震荡下不应产生汽机主控前馈量,一般设置为±0.1-0.3MW,饱和非线性环节的具体函数如下:
Figure FDA0000452306350000011
4.如权利要求1所述一种PROPR模式下火电机组负荷优化控制方法,其特征是,所述步骤三的具体步骤为:
将偏差量送至一高低限报警模块中,设置上下限,越限则输出“1”,反之为“0”;将高低限报警模块的输出作为一切换器T2的选择开关,修正后的负荷变化率值和原有的负荷变化率值同时送至切换器T2的两个切换输入端中,若为“1”,则将原有的负荷变化率值送出,为“0”,则切换器T2将修正后的负荷变化率值送出。
5.如权利要求4所述的一种PROPR模式下火电机组负荷优化控制方法,其特征是,所述高低限报警模块的上下限H/L的设置原则为:按规定300MW机组在PROPR模式下单次负荷调节的幅值为±3MW,其计算调节速率的负荷取点在负荷实际变化至±2MW之前,因此将H和L分别设置为2和-2。
6.如权利要求5所述的一种PROPR模式下火电机组负荷优化控制方法,其特征是,所述切换器T2的两个输入值设置原则为:一个为原有负荷变化率值,另一个为原有负荷变化率值的60%-70%,也就是说,当实际功率变化将至负荷指令值时,以较低的负荷变化率接近,进而提高调节精度指标。
CN201410001773.1A 2014-01-02 2014-01-02 Propr模式下火电机组负荷优化控制方法 Active CN103713613B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410001773.1A CN103713613B (zh) 2014-01-02 2014-01-02 Propr模式下火电机组负荷优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410001773.1A CN103713613B (zh) 2014-01-02 2014-01-02 Propr模式下火电机组负荷优化控制方法

Publications (2)

Publication Number Publication Date
CN103713613A true CN103713613A (zh) 2014-04-09
CN103713613B CN103713613B (zh) 2017-01-18

Family

ID=50406664

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410001773.1A Active CN103713613B (zh) 2014-01-02 2014-01-02 Propr模式下火电机组负荷优化控制方法

Country Status (1)

Country Link
CN (1) CN103713613B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900821A (zh) * 2014-04-22 2014-07-02 国家电网公司 一种发电机组一次调频功能优化试验方法
CN104089270A (zh) * 2014-07-11 2014-10-08 国家电网公司 一种发电机组锅炉负荷控制优化调整试验方法
CN106094531A (zh) * 2016-07-22 2016-11-09 河北省电力建设调整试验所 一种适用于agc工况的燃料量智能前馈控制器
CN107591847A (zh) * 2017-08-04 2018-01-16 西安五常电力技术有限公司 一种利用变参数的方式调节水电机组agc的方法
CN107976900A (zh) * 2017-11-21 2018-05-01 哈尔滨沃华智能发电设备有限公司 一种火电机组agc系统的前馈补偿设计方法
CN108599190A (zh) * 2017-12-25 2018-09-28 赫普科技发展(北京)有限公司 一种对电网进行调频的方法、装置和系统
CN110513158A (zh) * 2019-09-04 2019-11-29 黑龙江苑博信息技术有限公司 汽轮机前馈的多级速率调节方法
CN110690711A (zh) * 2019-10-25 2020-01-14 国电宝鸡发电有限责任公司 一种基于agc方式的机组负荷指令逻辑调控方法
CN114151207A (zh) * 2021-11-03 2022-03-08 中山嘉明电力有限公司 一种燃机机组快速变负荷控制方法
CN114527659A (zh) * 2022-02-16 2022-05-24 华北电力科学研究院有限责任公司 一种火电机组主动模式识别控制方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1230817A (zh) * 1998-02-11 1999-10-06 唐寅生 电力系统无功功率自动调控方法
US20050055137A1 (en) * 2001-09-13 2005-03-10 Anders Andren Method and system to calculate a demand for energy
CN101206754A (zh) * 2006-12-21 2008-06-25 北京华电天仁电力控制技术有限公司 基于多种约束规则的火电厂厂级负荷优化分配方法
CN101788788A (zh) * 2010-01-05 2010-07-28 华北电力大学 基于电网状态识别的多目标机组组合智能优化系统
CN103345205A (zh) * 2013-05-21 2013-10-09 国家电网公司 一种基于数据校核的agc负荷优化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1230817A (zh) * 1998-02-11 1999-10-06 唐寅生 电力系统无功功率自动调控方法
US20050055137A1 (en) * 2001-09-13 2005-03-10 Anders Andren Method and system to calculate a demand for energy
CN101206754A (zh) * 2006-12-21 2008-06-25 北京华电天仁电力控制技术有限公司 基于多种约束规则的火电厂厂级负荷优化分配方法
CN101788788A (zh) * 2010-01-05 2010-07-28 华北电力大学 基于电网状态识别的多目标机组组合智能优化系统
CN103345205A (zh) * 2013-05-21 2013-10-09 国家电网公司 一种基于数据校核的agc负荷优化控制方法

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900821A (zh) * 2014-04-22 2014-07-02 国家电网公司 一种发电机组一次调频功能优化试验方法
CN104089270A (zh) * 2014-07-11 2014-10-08 国家电网公司 一种发电机组锅炉负荷控制优化调整试验方法
CN104089270B (zh) * 2014-07-11 2016-03-30 国家电网公司 一种发电机组锅炉负荷控制优化调整试验方法
CN106094531A (zh) * 2016-07-22 2016-11-09 河北省电力建设调整试验所 一种适用于agc工况的燃料量智能前馈控制器
CN107591847B (zh) * 2017-08-04 2020-05-01 西安五常电力技术有限公司 一种利用变参数的方式调节水电机组agc的方法
CN107591847A (zh) * 2017-08-04 2018-01-16 西安五常电力技术有限公司 一种利用变参数的方式调节水电机组agc的方法
CN107976900A (zh) * 2017-11-21 2018-05-01 哈尔滨沃华智能发电设备有限公司 一种火电机组agc系统的前馈补偿设计方法
CN107976900B (zh) * 2017-11-21 2020-08-07 哈尔滨沃华智能发电设备有限公司 一种火电机组agc系统的前馈补偿设计方法
CN108599190A (zh) * 2017-12-25 2018-09-28 赫普科技发展(北京)有限公司 一种对电网进行调频的方法、装置和系统
CN110513158A (zh) * 2019-09-04 2019-11-29 黑龙江苑博信息技术有限公司 汽轮机前馈的多级速率调节方法
CN110690711A (zh) * 2019-10-25 2020-01-14 国电宝鸡发电有限责任公司 一种基于agc方式的机组负荷指令逻辑调控方法
CN114151207A (zh) * 2021-11-03 2022-03-08 中山嘉明电力有限公司 一种燃机机组快速变负荷控制方法
CN114151207B (zh) * 2021-11-03 2023-12-22 中山嘉明电力有限公司 一种燃机机组快速变负荷控制方法
CN114527659A (zh) * 2022-02-16 2022-05-24 华北电力科学研究院有限责任公司 一种火电机组主动模式识别控制方法及系统
CN114527659B (zh) * 2022-02-16 2024-09-27 华北电力科学研究院有限责任公司 一种火电机组主动模式识别控制方法及系统

Also Published As

Publication number Publication date
CN103713613B (zh) 2017-01-18

Similar Documents

Publication Publication Date Title
CN103713613A (zh) Propr模式下火电机组负荷优化控制方法
CN103345205B (zh) 一种基于数据校核的agc负荷优化控制方法
CN103378601B (zh) 一种基于bang-bang控制的一次调频方法及装置
CN102541028B (zh) 一种煤质变化下超临界机组agc优化控制方法
CN103174471B (zh) 一种火电机组变参数负荷优化控制方法
CN102748080B (zh) 基于主蒸汽压力变化的火电机组负荷控制方法
CN102654751A (zh) 基于非线性控制和模糊控制的协调控制方法
CN104810851B (zh) 火电机组agc方式下负荷调节速率动态调整的方法及装置
CN103104931B (zh) 一种基于锅炉一次风压动态补偿的一次调频方法及系统
CN105135409B (zh) 基于一次调频动作幅值的超临界机组锅炉主控控制方法
CN103294030A (zh) 一种dcs控制方法及smith控制器
CN111399448B (zh) 一种汽机跟随方式下燃煤机组一次调频控制方法
CN106773644B (zh) 一种基于供热抽汽量变化的agc控制系统及其方法
CN109378833A (zh) 一种通过控制汽轮机抽汽量实现机组快速调频的方法
CN104269871B (zh) 提高火力发电机组一次调频正确动作合格率的方法
CN103378784B (zh) 一种负荷指令动态补偿方法及系统
CN104600708A (zh) 含svg的风电场自动电压控制分配方法
Boddepalli et al. Design and analysis of firefly algorithm based PID controller for automatic load frequency control problem
CN105048514B (zh) 一种基于调节性能一致性分析的机组协调控制方法
CN105375511A (zh) 一种一次调频转速不等率函数的设置方法
CN105974800A (zh) 一种反馈前馈控制系统
CN204741287U (zh) 一种火电机组agc方式下负荷调节速率动态调整的装置
CN103670539B (zh) 补偿发电机组动态特性联合调频控制方法、系统与装置
CN110134003A (zh) 自适应agc工况的火电机组变负荷动态前馈控制方法
Kumar et al. Deregulated AGC scheme using dynamic programming controller

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Li Jun

Inventor after: Zhou Changlai

Inventor after: Zhao Yan

Inventor before: Li Jun

Inventor before: Zhou Changlai

COR Change of bibliographic data
C14 Grant of patent or utility model
GR01 Patent grant