CN103679282B - 风电功率爬坡的预测方法 - Google Patents
风电功率爬坡的预测方法 Download PDFInfo
- Publication number
- CN103679282B CN103679282B CN201310457433.5A CN201310457433A CN103679282B CN 103679282 B CN103679282 B CN 103679282B CN 201310457433 A CN201310457433 A CN 201310457433A CN 103679282 B CN103679282 B CN 103679282B
- Authority
- CN
- China
- Prior art keywords
- climbing
- wind power
- time
- ramp
- prediction
- 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.)
- Expired - Fee Related
Links
- 230000009194 climbing Effects 0.000 title claims abstract description 166
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000003203 everyday effect Effects 0.000 claims abstract description 7
- 238000009499 grossing Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013277 forecasting method Methods 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000714 time series forecasting Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000003828 downregulation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种风电功率爬坡的预测方法,主要包括以下步骤:据时序先后逐点计算超短期风电功率爬坡持续时间、爬坡幅值、爬坡方向及实时爬坡率,比较与阈值判断是否为风电功率爬坡事件;获取每日更新的短期风电功率预测序列,对其逐点进行平滑处理,根据超过设定阈值部分的极大值确定功率爬坡事件的爬坡中心时间,根据附近超过阈值的临界点计算爬坡持续时间、爬坡幅值及爬坡方向,实现短期风电功率极端爬坡事件预测;以及分别对超短期风电功率爬坡和短期风电功率爬坡的历史预测结果和历史观测结果进行统计,实现风电功率爬坡预测。
Description
技术领域
本发明属于电力***预测与控制技术领域。尤其涉及一种超短期/短期风电功率爬坡预测方法。
背景技术
受风资源变化特性支配,风电功率输出具有很强的随机性和波动性。风电功率爬坡事件指较短时间内风电出力发生较大变化,大规模风电集中并网后的风电功率爬坡可在短时间内变化上百兆瓦,是电力***调度运行所面临的重要挑战,一旦爬坡幅度或爬坡速率超过***备用可以承受的范围,将对***的安全稳定构成威胁。
对风电功率爬坡进行预报可以帮助***提前采取措施应对即将发生的功率爬坡事件。风电功率爬坡事件的基本要素是爬坡持续时间、爬坡幅值和爬坡方向(简称为“爬坡三要素”)。根据预报结果,对于功率上爬坡事件,调度机构可以调整发电计划以吸纳更多的清洁能源,或者在下调备用容量不足时采取限电等必要措施;对于功率下爬坡事件,可以按照功率爬坡的规模大小提前准备投入各级备用,必要时采取切负荷等措施保证***安全。
传统风电功率预测技术以输出规定时段的预测功率序列为目的,而风电功率爬坡预报侧重对爬坡事件的识别和预判。由于功率爬坡事件涉及风资源的快速变化,采用的预测方法主要基于数值模式对天气变化物理过程的滚动模拟。国内已初步建立风电功率预测***技术规范,规定预测***提供每15分钟滚动的超短期未来4小时功率预报和按日上报的短期日前功率预报,但不包含对功率爬坡预报的要求,因此国内现有的风电功率预测***并不具备该项功能。另一方面,国内风电功率预测***所使用的数值天气预报服务一般由气象机构按日提供1-2次,尚不具备滚动数值预报服务,难以依靠实时气象预报对功率爬坡进行模拟。因此,如何依据现有风电功率预测***运行模式在超短期/短期时间尺度实现对风电功率爬坡的预报具有十分重要的工程价值。
发明内容
综上所述,确有必要提供一种能够预测风电功率爬坡的预测方法。
一种风电功率爬坡的预测方法,主要包括以下步骤:获取实时更新的超短期风电功率预测序列,根据时序先后逐点计算超短期风电功率爬坡持续时间爬坡幅值爬坡方向dm及实时爬坡率根据设定阈值判断是否为风电功率爬坡事件;获取每日更新的短期风电功率预测序列,对其逐点进行平滑处理,根据平滑后序列超过设定阈值Pval部分的极大值确定功率爬坡事件的爬坡中心时间tc,m,根据附近超过阈值的临界点计算爬坡持续时间爬坡幅值及爬坡方向dm,实现短期风电功率极端爬坡事件预测;以及分别对超短期风电功率爬坡和短期风电功率爬坡的历史预测结果和历史观测结果进行统计,计算预报完成率、预报准确率和爬坡综合指标,实现风电功率爬坡预测。
相对于现有技术,本发明提供的风电功率爬坡预报方法,基于现行风电功率预测***运行模式实现超短期/短期风电功率爬坡预报,通过设定不同含义爬坡阈值,在不同时间尺度将时间序列型预报转化为事件型预报,该功率爬坡预报方法可作为其附加功能模块或外置模块,具有较高灵活性,且无需额外的输入数据支持,经济成本低。
附图说明
图1为本发明提供的超短期风电功率爬坡的预测方法的流程图。
图2为本发明提供的超短期风电功率爬坡计算流程图。
图3为本发明提供的短期风电功率极端爬坡事件预测流程图。
图4为本发明提供的短期风电功率极端爬坡判断示意图。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
请参阅图1,图1为本发明提供的风电功率爬坡的预测方法流程图,主要包括如下步骤:
步骤S10,获取实时更新的超短期风电功率预测序列,根据时序先后逐点计算超短期风电功率爬坡持续时间、爬坡幅值、爬坡方向,得到实时爬坡率,比较实时爬坡率与设定阈值,判断是否为风电功率爬坡事件;
步骤S20,获取每日更新的短期风电功率预测序列,对其逐点进行平滑处理,根据平滑后序列超过设定阈值部分的极大值确定功率爬坡事件的爬坡中心时间,根据附近超过阈值的临界点计算爬坡持续时间、爬坡幅值、爬坡方向及实时爬坡率,实现短期风电功率极端爬坡事件预测;
步骤S30,分别对超短期风电功率爬坡和短期风电功率爬坡的历史预测结果和历史观测结果进行统计,计算预报完成率、预报准确率和爬坡综合指标,实现风电功率爬坡预测。
在步骤S10中,请一并参阅图2,所述爬坡持续时间、爬坡幅值、及爬坡方向可通过以下方法计算:
对于第k个预测点,其对应时刻为tk,Δt为事先设定的爬坡时间,取Δt小于等于1小时。max(p[tk,tk+Δt])为[tk,tk+Δt]时段内风电功率预测最大值,对应时刻为min(p[tk,tk+Δt])为[tk,tk+Δt]时段内风电功率预测最小值,对应时刻为则爬坡持续时间为:
爬坡幅值为:
爬坡方向dk为:
则对应tk的未来给定时段爬坡率为:
如果则判定其为功率爬坡事件。其中,为风电爬坡率阈值,一般作为旋转备用的火电机组,其上/下功率爬坡上限一般为3-4MW/min,故根据实际备用情况,一般应取不大于备用机组功率爬坡上限的80%,即2-3MW/min。
在步骤S20中,请一并参阅图3,所述短期风电功率极端爬坡事件预测可通过以下步骤进行。
首先,对短期风电功率预测序列{pt}(一般覆盖未来24-48小时)逐点进行平滑处理:
平滑阶数n一般取5≤n≤10,是{pt}的平滑结果,t为时刻,h为一段时间。
其次,根据风电爬坡阈值Pval判断功率爬坡事件,并计算其爬坡持续时间、爬坡幅值、及爬坡方向:
就极端功率爬坡事件而言,Pval一般取20%-30%风电总开机容量。
请一并参阅图4,所述爬坡持续时间、爬坡幅值、及爬坡方向的计算包括以下步骤:
由的局部极大值确定第m次功率爬坡事件的爬坡中心时间tc,m;
通过所述极大值附近超过阈值的两个临界点和可对应确定功率爬坡起止时刻t1,m和t2,m(t1,m<t2,m),则爬坡持续时间为:
所述爬坡幅值为:
所述爬坡方向dm为:
在步骤S30中,超短期/短期风电功率爬坡预测都属于事件型预测,可采用列联表对其历史预测结果和历史观测结果进行统计。
对于超短期风电功率爬坡预测,tk时刻历史预测爬坡率和历史观测爬坡率的绝对值均超过且爬坡方向相同的情况为成功;tk时刻历史预测爬坡率的绝对值超过历史观测爬坡率的绝对值未超过或超过但爬坡方向不同的情况为错报;tk时刻历史观测爬坡率的绝对值超过历史预测爬坡率的绝对值未超过或超过但爬坡方向不同的情况为漏报。
对于短期风电功率爬坡预测,对于每一次爬坡事件,根据爬坡起止时刻定义爬坡区间[t1,m,t2,m],该区间包含爬坡中心时间tc,m,则一次历史预测爬坡事件爬坡中心时间若存在一次历史观测爬坡事件爬坡中心时间属于且属于则统计为成功;没有被统计为成功的其余历史预测爬坡事件统计为错报;相应地,没有被统计为成功的其余历史观测爬坡事件统计为漏报。
完成所有历史结果统计后,计算评价指标:
评价指标越接近1,功率爬坡预测效果越好。
实施例一
作为具体的实施例,所述的超短期/短期风电功率爬坡预报方法基于我国西北某风电场风电功率预测***历史数据进行虚拟预测,具体步骤如下:
(1)计算超短期风电功率实时爬坡率。
该风电场风电功率预测***每15分钟滚动超短期未来4小时全场功率预报,对于每次更新的预报序列,设定爬坡时间Δt为1小时,按照时序先后,依据公式(1)-(4),逐点计算其Δt时段内的爬坡持续时间、爬坡幅值、及爬坡方向和爬坡率。风电爬坡率阈值取为2MW/min。
(2)短期风电功率极端爬坡事件预测。
该风电场风电功率预测***每天上报未来3天的全场风电功率功率预报结果,对于每天更新的预报序列,设定平滑阶数n为9,得到平滑序列,设定风电爬坡阈值Pval为25%风电总开机容量,从而依据Pval判定极端功率爬坡事件。对于每次爬坡事件,根据公式(7)-(9)计算其爬坡持续时间、爬坡幅值、及爬坡方向。
(3)风电功率爬坡预测。
采用列联表对其半年的历史预测结果和历史观测结果进行统计,计算预报完成率,预报准确率和爬坡综合指标。该风电场2012年上半年的统计结果为:
超短期风电功率爬坡预报:预报完成率76%,预报准确率70%,爬坡综合指标73%;
短期风电功率爬坡预报:预报完成率83%,预报准确率74%,爬坡综合指标78%。
本发明提供的风电功率爬坡预报方法,基于现行风电功率预测***运行模式实现超短期/短期风电功率爬坡预报,通过不同含义爬坡阈值的设定,在不同时间尺度将时间序列型预报转化为事件型预报:在超短期时间尺度着重实时判断风电功率爬坡速率是否超过旋转备用机组承受能力;在短期时间尺度重点对大规模极端功率爬坡事件进行日前预测。对于各类风电功率预测***,该功率爬坡预报方法可作为其附加功能模块或外置模块,具有较高灵活性,且无需额外的输入数据支持,经济成本低。
另外,本领域技术人员还可在本发明精神内作其它变化,当然这些依据本发明精神所作的变化,都应包含在本发明所要求保护的范围内。
Claims (10)
1.一种风电功率爬坡的预测方法,主要包括以下步骤:
步骤S10,获取实时更新的超短期风电功率预测序列,根据时序先后逐点计算超短期风电功率爬坡持续时间爬坡幅值爬坡方向dm,得到实时爬坡率比较与阈值判断是否为风电功率爬坡事件;
步骤S20,获取每日更新的短期风电功率预测序列,对其逐点进行平滑处理,根据平滑后序列超过设定阈值Pval部分的极大值确定功率爬坡事件的爬坡中心时间tc,m,根据附近超过阈值的临界点计算爬坡持续时间爬坡幅值及爬坡方向dm,实现短期风电功率极端爬坡事件预测;以及
步骤S30,分别对超短期风电功率爬坡和短期风电功率爬坡的历史预测结果和历史观测结果进行统计,计算预报完成率、预报准确率和爬坡综合指标,实现风电功率爬坡预测。
2.如权利要求1所述的风电功率爬坡的预测方法,其特征在于,所述爬坡持续时间、爬坡幅值、及爬坡方向通过以下方法计算:
对于第k个预测点,其对应时刻为tk,Δt为事先设定的爬坡时间,取Δt小于等于1小时;
设max(p[tk,tk+Δt])为[tk,tk+Δt]时段内风电功率预测最大值,对应时刻为
设min(p[tk,tk+Δt])为[tk,tk+Δt]时段内风电功率预测最小值,对应时刻为则爬坡持续时间为:
爬坡幅值为:
爬坡方向dk为:
3.如权利要求2所述的风电功率爬坡的预测方法,其特征在于,对应tk的未来给定时段爬坡率为:
如果则判定其为功率爬坡事件。
4.如权利要求1为所述的风电功率爬坡的预测方法,其特征在于,所述阈值小于等于备用机组功率爬坡上限的80%。
5.如权利要求1所述的风电功率爬坡的预测方法,其特征在于,所述阈值Pval为风电总开机容量的20%-30%。
6.如权利要求1所述的风电功率爬坡的预测方法,其特征在于,所述短期风电功率极端爬坡事件预测包括以下子步骤:
首先,对短期风电功率预测序列{pt}逐点进行平滑处理:
其中,平滑阶数n取5≤n≤10,是{pt}的平滑结果,t为时刻,h为一段时间;
其次,根据风电爬坡阈值Pval判断功率爬坡事件,并计算其爬坡持续时间、爬坡幅值、及爬坡方向。
7.如权利要求6所述的风电功率爬坡的预测方法,其特征在于,所述爬坡持续时间、爬坡幅值、及爬坡方向的计算包括以下子步骤:
由的局部极大值确定第m次功率爬坡事件的爬坡中心时间tc,m;
通过所述极大值附近超过阈值的两个临界点和对应确定功率爬坡起止时刻t1,m和t2,m(t1,m<t2,m),则爬坡持续时间为:
所述爬坡幅值为:
所述爬坡方向dm为:
8.如权利要求1所述的风电功率爬坡的预测方法,其特征在于,采用列联表对风电功率爬坡的历史预测结果和历史观测结果进行统计。
9.如权利要求8所述的风电功率爬坡的预测方法,其特征在于,对于超短期风电功率爬坡预测,tk时刻历史预测爬坡率和历史观测爬坡率的绝对值均超过且爬坡方向相同的情况为成功;tk时刻历史预测爬坡率的绝对值超过历史观测爬坡率的绝对值未超过或超过但爬坡方向不同的情况为错报;tk时刻历史观测爬坡率的绝对值超过历史预测爬坡率的绝对值未超过或超过但爬坡方向不同的情况为漏报。
10.如权利要求8所述的风电功率爬坡的预测方法,其特征在于,对于短期风电功率爬坡预测,对于每一次爬坡事件,根据爬坡起止时刻定义爬坡区间[t1,m,t2,m],该区间包含爬坡中心时间tc,m,则对于一次历史预测爬坡事件爬坡中心时间若存在一次历史观测爬坡事件爬坡中心时间 属于且属于则统计为成功;没有被统计为成功的其余历史预测爬坡事件统计为错报;相应地,没有被统计为成功的其余历史观测爬坡事件统计为漏报。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310457433.5A CN103679282B (zh) | 2013-09-30 | 2013-09-30 | 风电功率爬坡的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310457433.5A CN103679282B (zh) | 2013-09-30 | 2013-09-30 | 风电功率爬坡的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103679282A CN103679282A (zh) | 2014-03-26 |
CN103679282B true CN103679282B (zh) | 2017-01-04 |
Family
ID=50316766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310457433.5A Expired - Fee Related CN103679282B (zh) | 2013-09-30 | 2013-09-30 | 风电功率爬坡的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103679282B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106897789B (zh) * | 2017-01-03 | 2020-11-27 | 内蒙古电力(集团)有限责任公司 | 一种风电有功功率爬坡事件检测方法及装置 |
CN106933778A (zh) * | 2017-01-22 | 2017-07-07 | 中国农业大学 | 一种基于爬坡事件特征识别的风电功率组合预测方法 |
CN107947200A (zh) * | 2017-12-04 | 2018-04-20 | 长沙理工大学 | 一种含风电的***频率偏差系数确定的方法 |
CN110264375A (zh) * | 2019-06-21 | 2019-09-20 | 哈尔滨工业大学 | 一种风电场群出力间歇性的多时间尺度定量刻画方法 |
CN110348637B (zh) * | 2019-07-12 | 2022-10-11 | 哈尔滨工业大学 | 一种考虑场-网因素的风电爬坡事件预警方法 |
CN111582557A (zh) * | 2020-04-20 | 2020-08-25 | 哈尔滨工业大学 | 基于变差函数的风电爬坡事件多级预警方法 |
CN112270439B (zh) * | 2020-10-28 | 2024-03-08 | 国能日新科技股份有限公司 | 超短期风电功率预测方法、装置、电子设备及存储介质 |
CN112529275B (zh) * | 2020-12-02 | 2023-09-22 | 中国矿业大学 | 基于特征提取和深度学习的风电爬坡事件预测方法 |
CN113298290A (zh) * | 2021-04-27 | 2021-08-24 | 宁波职业技术学院 | 一种基于图像生成的风电功率预测方法与*** |
CN114254805A (zh) * | 2021-11-22 | 2022-03-29 | 华北电力大学 | 爬坡事件的时间窗口识别方法、装置、设备及存储介质 |
CN117196122B (zh) * | 2023-11-02 | 2024-02-06 | 湖南赛能环测科技有限公司 | 一种基于风电功率爬坡时长的风电场调整方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7613548B2 (en) * | 2006-01-26 | 2009-11-03 | General Electric Company | Systems and methods for controlling a ramp rate of a wind farm |
CN103248050A (zh) * | 2013-05-22 | 2013-08-14 | 山东大学 | 一种基于竞争博弈的风电爬坡协同控制***及方法 |
CN103279804A (zh) * | 2013-04-29 | 2013-09-04 | 清华大学 | 超短期风电功率的预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8930299B2 (en) * | 2010-12-15 | 2015-01-06 | Vaisala, Inc. | Systems and methods for wind forecasting and grid management |
-
2013
- 2013-09-30 CN CN201310457433.5A patent/CN103679282B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7613548B2 (en) * | 2006-01-26 | 2009-11-03 | General Electric Company | Systems and methods for controlling a ramp rate of a wind farm |
CN103279804A (zh) * | 2013-04-29 | 2013-09-04 | 清华大学 | 超短期风电功率的预测方法 |
CN103248050A (zh) * | 2013-05-22 | 2013-08-14 | 山东大学 | 一种基于竞争博弈的风电爬坡协同控制***及方法 |
Non-Patent Citations (3)
Title |
---|
《A Survey of Wind Power Ramp Forecasting》;Tinghui Ouyang,Xiaoming Zha,Liang Qin;《Energy and Power Engineering》;20130731;第5卷(第4B期);全文 * |
《基于统计分析和多支持向量机的风电功率坡度事件分类预测》;李福东等;《上海交通大学学报》;20121231;第46卷(第12期);全文 * |
《抑制风电爬坡率的风储联合优化控制方法》;王颖等;《电力***自动化》;20130710;第37卷(第13期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103679282A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103679282B (zh) | 风电功率爬坡的预测方法 | |
Cui et al. | A novel event detection method using PMU data with high precision | |
CN103543403B (zh) | 电力***机组一次调频能力检测方法 | |
CN104463356A (zh) | 一种基于多维信息人工神经网络算法的光伏发电功率预测方法 | |
CN106203693A (zh) | 一种风电场输出功率爬坡事件预测的***及方法 | |
CN105140967A (zh) | 一种含新能源电力***调峰需求的评估方法 | |
CN106611243A (zh) | 一种基于garch模型的风速预测残差修正方法 | |
CN111191832A (zh) | 一种台风灾害配电网杆塔故障预测方法及*** | |
CN102904248B (zh) | 基于风电出力不确定性集合的电力***调度方法 | |
CN104484728B (zh) | 一种电网安全综合指标体系架构方法 | |
CN104574211A (zh) | 基于风险源的电网调度操作风险预警方法和*** | |
Li et al. | Generating wind power time series based on its persistence and variation characteristics | |
CN116307257B (zh) | 一种特定天气下的风电场的输出功率预测方法及装置 | |
Lijuan et al. | A novel model for wind power forecasting based on Markov residual correction | |
CN103473461B (zh) | 基于数据特征提取的风功率预测误差估计方法 | |
CN112182952B (zh) | 一种提升电力***弹性的多目标优化调度方法 | |
CN205091957U (zh) | 基于实时仿真的在线调度员培训*** | |
Li et al. | Proportion optimization of wind and solar power capacity for regional power network | |
Lu | Multi-step ahead ultra-short-term wind power forecasting based on time series analysis | |
Wu et al. | Short-term wind power prediction based on time series analysis model | |
Ma et al. | A robust interval economic dispatch method accommodating significant wind power integration | |
CN105184411A (zh) | 一种基于最大Lyapunov指数-校正的风电功率实时预测方法 | |
CN104574216A (zh) | 一种基于wams数据的风电出力特性分析方法 | |
Zhao et al. | Real-time Inertia Estimation Using Probing Method and PMU Measurement | |
CN104166939A (zh) | 基于情景分析的配电网运行风险预警方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170104 |
|
CF01 | Termination of patent right due to non-payment of annual fee |