CN107388599B - 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法 - Google Patents

一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法 Download PDF

Info

Publication number
CN107388599B
CN107388599B CN201710653340.8A CN201710653340A CN107388599B CN 107388599 B CN107388599 B CN 107388599B CN 201710653340 A CN201710653340 A CN 201710653340A CN 107388599 B CN107388599 B CN 107388599B
Authority
CN
China
Prior art keywords
shadow
column
point
shielding
mirror
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
Application number
CN201710653340.8A
Other languages
English (en)
Other versions
CN107388599A (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.)
Lanzhou Jiaotong University
Original Assignee
Lanzhou Jiaotong University
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 Lanzhou Jiaotong University filed Critical Lanzhou Jiaotong University
Priority to CN201710653340.8A priority Critical patent/CN107388599B/zh
Publication of CN107388599A publication Critical patent/CN107388599A/zh
Application granted granted Critical
Publication of CN107388599B publication Critical patent/CN107388599B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24SSOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
    • F24S23/00Arrangements for concentrating solar-rays for solar heat collectors
    • F24S23/30Arrangements for concentrating solar-rays for solar heat collectors with lenses
    • F24S23/31Arrangements for concentrating solar-rays for solar heat collectors with lenses having discontinuous faces, e.g. Fresnel lenses
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/40Solar thermal energy, e.g. solar towers

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Thermal Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)
  • Optical Elements Other Than Lenses (AREA)

Abstract

本发明涉及光伏发电技术领域,特别是一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法。包括如下步骤:A镜场坐标系确定、B遮挡阴影坐标系确定、C遮挡阴影分析计算、D偏移影响分析计算、E阴影与遮挡效率。该方法更加贴近实际使用环境,可以准确计算线性菲涅尔式聚光镜场的阴影与遮挡效率,为分析线性菲涅尔式聚光系统的效率提供了可靠依据。进而修正线性菲涅尔式聚光镜场反射镜的分布方向和角度,进一步提升镜场转换效率。

Description

一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法
技术领域
本发明涉及光伏发电技术领域,特别是一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法。
背景技术
太阳能作为一种理想的可再生能源已备受关注,但是能流密度低,提高利用效率的最佳技术路线是先聚光再利用。根据聚光形式,太阳能聚光系统主要有槽式、线性菲涅尔式、塔式和碟式等。线性菲涅尔式聚光系统主要由光学聚光器、接收器、太阳跟踪控制装置三部分构成,具有结构简单、风阻小、成本低、土地利用率高等优点,正逐渐在大规模电站中得到应用。
目前,针对已公开的专利(名称:线性菲涅尔式聚光系统镜场的优化布置方法、专利号:201410439026.6、授权公告号:CN104236122B)。本发明针对实际需要,发明了一种线性菲涅尔式聚光系统镜场的优化布置方法。由于南北方向的镜场总体光学效率优于东西方向,因此本发明以南北方向的镜场为研究对象,提供 了一种以CPC最大接受半角控制镜场高宽比进而进行镜场无阴影布置的方法。本发明的方法 可以针对实际情况快速准确地计算出确定CPC、反射镜宽度和列数的线性菲涅尔式聚光集热 系统在特定时间(系统无阴影工作时间)内无阴影影响的反射镜布置间距。
针对南北布置的采用单管接收器的线性菲涅尔式聚光系统,推导了镜场无阴影布置的数学表达式,并给出了数值计算方法。在实际工程应用中,为了便于加工安装并采用简单的驱动系统,反射镜一般采用等间距布置,这必将带来相邻反射镜之间的阴影与遮挡影响。而对于线性菲涅尔式聚光系统,相邻反射镜之间的阴影与遮挡效率直接影响着整个系统的光学效率,而系统光学效率的分析计算是设计电站聚光系统的依据。因此,准确分析相邻反射镜之间的阴影与遮挡至关重要。现有设备和技术中多根据镜场无阴影进行设定进而布设镜场,但实际场地中无阴影几乎不可能实现,使用现有设备和技术进行测试和布设镜场时均未考虑或很少考虑阴影遮挡问题,进而造成数据不准确影响使用效率。
发明内容
本发明解决其现有技术的不足提供一种分析数据精准、分析维度多样、布设更加合理的线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法。
本发明解决其技术问题所采用的技术方案为:
一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法,包括如下步骤:
A、镜场坐标系确定
反射镜列数为N的镜场,从左侧至右侧依次编号为1、2、...、N,第列反射镜中心与镜场中心的距离为,第列反射镜与水平面的倾斜夹角为;反射镜宽度为;接收器放置在距反射镜所在平面的高处;
某时刻太阳光入射角为,则该时刻的跟踪倾角由下式求得:
列反射镜的坐标方程为
式中,
列反射镜的左、右端点坐标分别为()和(),则
B、遮挡阴影坐标系确定
遮挡阴影是指太阳入射光被相邻的反射镜所阻挡,遮挡阴影分析针对太阳入射光;
令经过第列反射镜左端点的太阳入射光线坐标方程为:
式中,
相邻的反射镜中过东侧左端点的太阳入射光线坐标为
若该点不在第列反射镜上,则无阴影影响;若该点在第列反射镜上,则形成了阴影影响,阴影长度为该点与第列反射镜右端点之间的距离,即
C、遮挡阴影分析计算
C1、接收器西侧镜场
令经过第列反射镜左端点的太阳反射光线方程为
式中,
相邻的反射镜中过西侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜右端点之间的距离,即
若N为偶数,则, 若N为奇数,则
C2、接收器东侧镜场
令经过第列反射镜右端点的太阳反射光线方程为
式中,
相邻的反射镜中过东侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜左端点之间的距离,即
若N为偶数,则;若N为奇数,则
若N为偶数,则都为0,若N为奇数,则为0。
还包括步骤:
D、偏移影响分析计算
阴影随时间变化产生偏移,阴影偏移中会出现空白阴影区域,其中
其中,为太阳方位角,求得时刻点阴影区域的中位线长度
同理求得空白阴影区域的中位线长度
E、阴影与遮挡效率
多个时刻点阴影与遮挡区域的中位线长度之和为,多个时刻点未形成阴影与遮挡区域的中位线长度之和为,得出测试时段内阴影与遮挡效率为
,式中为镜场长度。
本发明的有益效果为:
本发明根据线性菲涅尔式聚光镜场反射镜单轴跟踪的特点,建立平面直角坐标系,推导出任意时刻每列反射镜的坐标方程,利用射线追踪法追踪某时刻入射至每列反射镜边缘和从边缘反射至接收器的太阳光线,得到阴影与遮挡的长度,利用几何光学可求得由于偏移未形成阴影与遮挡区域的计算公式,进而求出线性菲涅尔式聚光镜场的阴影与遮挡效率。该方法更加贴近实际使用环境,可以准确计算线性菲涅尔式聚光镜场的阴影与遮挡效率,为分析线性菲涅尔式聚光系统的效率提供了可靠依据。进而修正线性菲涅尔式聚光镜场反射镜的分布方向和角度,进一步提升镜场转换效率。
附图说明
图1为线性菲涅尔式聚光镜场示意图;
图2为线性菲涅尔式聚光镜场阴影示意图;
图3为实施例中反射镜阴影情况示意图;
图4为实施例中反射镜遮挡阴影情况示意图;
图5为日光阴影与遮挡效率变化示意图。
具体实施方式
一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法,包括如下步骤:
A、镜场坐标系确定
反射镜列数为N的镜场,从左侧至右侧依次编号为1、2、...、N,第 列反射镜中心与镜场中心的距离为,第列反射镜与水平面的倾斜夹角为;反射镜宽度为;接收器放置在距反射镜所在平面的高处;
某时刻太阳光入射角为,则该时刻的跟踪倾角由下式求得:
列反射镜的坐标方程为
式中,
列反射镜的左、右端点坐标分别为()和(),则
B、遮挡阴影坐标系确定
遮挡阴影是指太阳入射光被相邻的反射镜所阻挡,遮挡阴影分析针对太阳入射光;
令经过第列反射镜左端点的太阳入射光线坐标方程为:
式中,
相邻的反射镜中过东侧左端点的太阳入射光线坐标为
若该点不在第列反射镜上,则无阴影影响;若该点在第列反射镜上,则形成了阴影影响,阴影长度为该点与第列反射镜右端点之间的距离,即
C、遮挡阴影分析计算
C1、接收器西侧镜场
令经过第列反射镜左端点的太阳反射光线方程为
式中,
相邻的反射镜中过西侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜右端点之间的距离,即
若N为偶数,则, 若N为奇数,则
C2、接收器东侧镜场
令经过第列反射镜右端点的太阳反射光线方程为
式中,
相邻的反射镜中过东侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜左端点之间的距离,即
若N为偶数,则;若N为奇数,则
若N为偶数,则都为0,若N为奇数,则为0。
还包括步骤:
D、偏移影响分析计算
阴影随时间变化产生偏移,阴影偏移中会出现空白阴影区域,其中
其中,为太阳方位角,求得时刻点阴影区域的中位线长度
同理求得空白阴影区域的中位线长度
E、阴影与遮挡效率
多个时刻点阴影与遮挡区域的中位线长度之和为,多个时刻点未形成阴影与遮挡区域的中位线长度之和为,得出测试时段内阴影与遮挡效率为
,式中为镜场长度。
实施例参照Puerto Errado 1电站,相关参数分别为:反射镜场东西方向布置,长800m,反射镜与集热管圆心距离为7.4m,反射镜宽度为0.8m,反射镜列数为21列,以0.16m等间距布置。反射镜自南向北依次编号为1-21,利用Matlab软件编程可计算春分日6:30时刻镜场的阴影与遮挡情况,此时刻1-12号反射镜存在阴影情况,具体数据表1所示。
表1 春分日6:30时刻镜场的阴影与遮挡情况
为了验证计算的正确性,利用一款由美国 NREL( National Renewable EnergyLaboratory) 开发,基于射线追踪法的,用于太阳能发电光学系统建模与性能分析的软件工具 SolTrace 对上述系统进行仿真计算,图3a、3b、3c和3d 分别为反射镜2、5、8和11的阴影情况。图4a和4b分别为反射镜2和4的遮挡情况。
由图3和图4可以看出,反射镜2、5、8、11的阴影长度和反射镜2、4的遮挡长度与本文计算完全一致。
对该镜场进行阴影与遮挡效率全天动态分析,该系统分别在春分、夏至、冬至日的镜场阴影与遮挡效率随时间的变化如图5所示。
由图5可知,典型日阴影与遮挡效率随时间的变化与文献[14]的结果基本吻合,春分日全天变化较为平缓,夏至日在日出、日落时段较小,随着太阳高度角的增加逐渐平稳且与春分日接近,冬至日效率较小,随时间的变化也较大。
计算春分、夏至、冬至日由于偏移未形成阴影与遮挡区域的中位线长度的平均值分别为为14.81mm、10.70mm、29.85mm,分别占镜场全长的0.0019%、0.0013%、0.0037%。可见由于偏移未形成阴影与遮挡区域很小。实际中,线性菲涅式聚光镜场的长度都在百米以上,因此阴影与遮挡分析时由于偏移未形成阴影与遮挡区域影响可以忽略不计。
以Puerto Errado 1电站为算例,分析了春分、夏至、冬至日的镜场阴影与遮挡效率随时间的变化,结果给表明春分日全天变化较为平缓,夏至日在日出、日落时段较小,随着太阳高度角的增加逐渐平稳且与春分日接近,冬至日效率较小,随时间的变化也较大;阴影与遮挡分析时由于偏移未形成阴影与遮挡区域影响可以忽略不计。通过实施例的方法可以准确计算线性菲涅尔式聚光镜场的阴影与遮挡效率,为分析线性菲涅尔式聚光系统的效率提供了可靠依据。

Claims (1)

1.一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法,其特征在于包括如下步骤:
A、镜场坐标系确定
反射镜列数为N的镜场,从左侧至右侧依次编号为1、2、...、N,第列反射镜中心与镜场中心的距离为,第列反射镜与水平面的倾斜夹角为;反射镜宽度为;接收器放置在距反射镜所在平面的高处;
某时刻太阳光入射角为,则该时刻的跟踪倾角由下式求得:
列反射镜的坐标方程为
式中,
列反射镜的左、右端点坐标分别为()和(),则
B、阴影分析计算
阴影是指太阳入射光被相邻的反射镜所阻挡,阴影分析针对太阳入射光;
令经过第列反射镜左端点的太阳入射光线坐标方程为:
式中,
相邻的反射镜中过东侧左端点的太阳入射光线坐标为
若该点不在第列反射镜上,则无阴影影响;若该点在第列反射镜上,则形成了阴影影响,阴影长度为该点与第列反射镜右端点之间的距离,即
C、遮挡分析计算
C1、接收器西侧镜场
令经过第列反射镜左端点的太阳反射光线方程为
式中,
相邻的反射镜中过西侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜右端点之间的距离,即
若N为偶数,则, 若N为奇数,则
C2、接收器东侧镜场
令经过第列反射镜右端点的太阳反射光线方程为
式中,
相邻的反射镜中过东侧镜场的太阳入射光线坐标为
若N为偶数,则;若N为奇数,则
若该点不在第列反射镜上,则无遮挡影响,遮挡长度为0;若该点在第列反射镜上,则形成了遮挡影响,遮挡长度为该点与第列反射镜左端点之间的距离,即
若N为偶数,则;若N为奇数,则
若N为偶数,则都为0,若N为奇数,则为0;
D、偏移影响分析计算
阴影随时间变化产生偏移,阴影偏移中会出现未形成阴影区域,其中
其中,为太阳方位角,为未形成阴影区域的纵向中位线;为第列反射镜未形成阴影区域南侧边缘与第n+1列反射镜南侧边缘的垂直距离,是两个垂点,求得时刻点未形成阴影区域的中位线长度
同理求得未形成遮挡区域的中位线长度
E、阴影与遮挡效率
多个时刻点阴影长度与遮挡长度之和为,多个时刻点未形成阴影与遮挡区域的中位线长度之和为,得出测试时段内阴影与遮挡效率为
,式中为镜场长度。
CN201710653340.8A 2017-08-02 2017-08-02 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法 Expired - Fee Related CN107388599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710653340.8A CN107388599B (zh) 2017-08-02 2017-08-02 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710653340.8A CN107388599B (zh) 2017-08-02 2017-08-02 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法

Publications (2)

Publication Number Publication Date
CN107388599A CN107388599A (zh) 2017-11-24
CN107388599B true CN107388599B (zh) 2019-11-05

Family

ID=60344127

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710653340.8A Expired - Fee Related CN107388599B (zh) 2017-08-02 2017-08-02 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法

Country Status (1)

Country Link
CN (1) CN107388599B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108304626B (zh) * 2018-01-16 2021-08-13 兰州交大常州研究院有限公司 一种基于SolTrace的线性菲涅尔式聚光器的建模方法
CN108775718B (zh) * 2018-03-26 2020-03-17 河海大学常州校区 双轴跟踪器最大聚光比的线性菲涅尔反射镜设计方法
CN108761780B (zh) * 2018-06-20 2020-11-10 兰州交通大学 一种用于聚光集热系统中聚光镜构造的光学建模方法
CN110658857A (zh) * 2019-10-08 2020-01-07 浙江正泰新能源开发有限公司 一种验证光伏跟踪系统跟踪准确度的方法及装置
CN112214871B (zh) * 2020-09-09 2022-11-01 东南大学 一种聚焦离子束辅助沉积加工的结构高度计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09189603A (ja) * 1996-01-08 1997-07-22 Kansai Electric Power Co Inc:The 日影の評価方法、および反射光の評価方法
CN102930160A (zh) * 2012-11-02 2013-02-13 浙江大学 塔式太阳能热电系统镜场光学效率的计算方法
CN104236122A (zh) * 2014-10-14 2014-12-24 兰州大成科技股份有限公司 线性菲涅尔式聚光系统镜场的优化布置方法
US9112078B2 (en) * 2010-06-09 2015-08-18 Sma Solar Technology Ag Method of recognizing and assessing shadowing events
CN105956331A (zh) * 2016-05-30 2016-09-21 浙江大学 塔式太阳能热发电系统中阴影和遮挡的快速判定方法
CN106951642A (zh) * 2017-03-22 2017-07-14 中国电建集团西北勘测设计研究院有限公司 一种新型的仿生型太阳能光热镜场的排布方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09189603A (ja) * 1996-01-08 1997-07-22 Kansai Electric Power Co Inc:The 日影の評価方法、および反射光の評価方法
US9112078B2 (en) * 2010-06-09 2015-08-18 Sma Solar Technology Ag Method of recognizing and assessing shadowing events
CN102930160A (zh) * 2012-11-02 2013-02-13 浙江大学 塔式太阳能热电系统镜场光学效率的计算方法
CN104236122A (zh) * 2014-10-14 2014-12-24 兰州大成科技股份有限公司 线性菲涅尔式聚光系统镜场的优化布置方法
CN105956331A (zh) * 2016-05-30 2016-09-21 浙江大学 塔式太阳能热发电系统中阴影和遮挡的快速判定方法
CN106951642A (zh) * 2017-03-22 2017-07-14 中国电建集团西北勘测设计研究院有限公司 一种新型的仿生型太阳能光热镜场的排布方法

Also Published As

Publication number Publication date
CN107388599A (zh) 2017-11-24

Similar Documents

Publication Publication Date Title
CN107388599B (zh) 一种线性菲涅尔式聚光镜场的阴影与遮挡分析优化布设方法
Sánchez-González et al. Aiming factor to flatten the flux distribution on cylindrical receivers
CN104236122B (zh) 线性菲涅尔式聚光系统镜场的优化布置方法
Mustafa et al. Analytical study of an innovated solar power tower (PS10) in Aswan
Barale et al. Optical design of a linear Fresnel collector for Sicily
Ghodbane et al. Optical numerical investigation of a solar power plant of parabolic trough collectors
Sánchez-González et al. Beam-down linear Fresnel reflector: BDLFR
Vouros et al. On the optimal shape of secondary reflectors for linear Fresnel collectors
US11644219B2 (en) Secondary reflectors for solar collectors and methods of making the same
Eddhibi et al. Design and analysis of a heliostat field layout with reduced shading effect in southern Tunisia
Fartaria et al. Simulation and computation of shadow losses of direct normal, diffuse solar radiation and albedo in a photovoltaic field with multiple 2-axis trackers using ray tracing methods
Rungasamy et al. 3-D CFD modeling of a slanted receiver in a compact linear Fresnel plant with etendue-matched mirror field
Li et al. Performance analysis of azimuth tracking fixed mirror solar concentrator
Rungasamy et al. Comparative study of the optical and economic performance of etendue-conserving compact linear Fresnel reflector concepts
Kim et al. Power performance assessment of PV blinds system considering self-shading effects
Xu et al. Performance comparison of solar parabolic trough system with glass and film reflector
Lee et al. Calculation of optical efficiency for the first central-receiver solar concentrator system in Korea
Bouarfa et al. Developing optical and thermal models with experimental validation of parabolic trough collector for moroccan industrial heat applications
Hu et al. Performance analysis and optimization of an integrated azimuth tracking solar tower
Li et al. Effects of receiver parameters on the optical efficiency of a fixed linear-focus Fresnel lens solar system with sliding adjustment
Giovannetti et al. Comparison of process heat collectors with respect to technical and economic conditions
Serrano-Arrabal et al. Dual-tower CSP plants: optical assessment and optimization with a novel cone-tracing model
CN115968442A (zh) 用于基于阵列级地形的回溯的系统和方法
Slootweg et al. Investigation into central receiver design for optimal optical and thermal performance
Georgiou et al. Optical evaluation of heliostat mirrors using caustics

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20191105

Termination date: 20200802