[go: up one dir, main page]

CN108306617A - 一种理想太阳能电池最大功率点参数的求解方法 - Google Patents

一种理想太阳能电池最大功率点参数的求解方法 Download PDF

Info

Publication number
CN108306617A
CN108306617A CN201810064916.1A CN201810064916A CN108306617A CN 108306617 A CN108306617 A CN 108306617A CN 201810064916 A CN201810064916 A CN 201810064916A CN 108306617 A CN108306617 A CN 108306617A
Authority
CN
China
Prior art keywords
solar cell
maximum power
power point
ideal
expression
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
CN201810064916.1A
Other languages
English (en)
Other versions
CN108306617B (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.)
Anhui Polytechnic University
Original Assignee
Anhui Polytechnic 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 Anhui Polytechnic University filed Critical Anhui Polytechnic University
Priority to CN201810064916.1A priority Critical patent/CN108306617B/zh
Publication of CN108306617A publication Critical patent/CN108306617A/zh
Application granted granted Critical
Publication of CN108306617B publication Critical patent/CN108306617B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02SGENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
    • H02S50/00Monitoring or testing of PV systems, e.g. load balancing or fault identification
    • H02S50/10Testing of PV devices, e.g. of PV modules or single PV cells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • 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/50Photovoltaic [PV] energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Photovoltaic Devices (AREA)
  • Control Of Electrical Variables (AREA)

Abstract

本发明公开了一种理想太阳能电池最大功率点参数的求解方法,包括:建立理想太阳能电池的数学模型并根据数据模型求取输出特性表达式;根据最大功率理论结合输出特性表达式推导出最大功率点的初始表达式;采用电池温度、日常强度来修正太阳能电池最大功率点的初始表达式得到太阳能电池最大功率点的求解公式。本发明的优点在于:(1)本解析求解方法是建立在严格的理论推导基础之上,避开了循环迭代算法,变量关系明确,求解方便,克服了一般数值求解方法的缺陷。(2)本解析求解方法能直观反映太阳电池温度、日照强度等因素对最大功率点参数的影响,实际工程应用也比较方便,而参数求解精度比一般的解析求解方法要高。

Description

一种理想太阳能电池最大功率点参数的求解方法
技术领域
本发明涉及太阳能光伏发电技术领域,尤其是涉及理想太阳能电池功率预测和参数的解析求解方法。
背景技术
太阳能电池由于能直接光-电转换、清洁无污染和高功率密度等优势,使得我国太阳能发电占总发电量的比例呈现出逐年指数增长的趋势。但是太阳能电池在最大功率点跟踪的过程中需要获得相应的最大功率点参数。另外,太阳能电池长时间暴露于阳光和恶劣天气之下,效能会逐渐下降,同时太阳能电池板会出现一定比例的故障,因此太阳能发电系统需要周期性的维护和设备更换。而太阳能电池组件故障的判断和设备更换也同样需要依据最大功率点参数。
但是,太阳能电池输出特性的数学表达式是包含指数的超越方程,并且表达式中多个变量的取值和日照强度、电池温度有关,使得太阳能电池的电气参数计算过程复杂、求解比较困难。常用的以Lambert W函数、模式搜索和人工蜂群等算法为基础的数值求解方法计算量大、收敛困难、计算结果受初始值影响较大。而另外一些通过中值定理和线性最小二乘法推导的解析求解公式具有参数计算误差较大的缺点。
因此,现有的技术有待改进和发展。
发明内容
本发明的目的在于克服现有技术的不足,提供一种理想太阳能电池最大功率点参数的求解方法,该方法求解最大功率电参数时变量关系直接、求解方便、参数求解精度高、工程实现方便。
为了实现上述目的,本发明采用的技术方案为:一种理想太阳能电池最大功率点参数的求解方法,包括如下步骤:
建立理想太阳能电池的数学模型并根据数据模型求取输出特性表达式;
根据最大功率理论结合输出特性表达式推导出最大功率点的初始表达式;
采用电池温度、日照强度来修正太阳能电池最大功率点的初始表达式得到太阳能电池最大功率点的求解公式。
理想太阳能电池模型为一个电流源和一个二极管并联。
根据理想太阳能电池的数学模型求取的输出特性表达式为:
其中,I、V分别是理想太阳能电池的输出电流和输出电压,Iph是代表电流源的光生电流,Io是二极管的反向饱和电流,q是一个电子的电荷量,k是玻尔兹曼常数,Ns是串联电池片的数量,T是太阳能电池的温度,n是太阳能电池的理想因数。
太阳能电池的光生电流Iph的计算公式为:
其中,S是日照强度,Sr是参考日照强度,Tr是太阳能电池参考温度,Iph_r是太阳能电池生产厂家提供的标准测试条件下太阳能电池的短路电流,μI是太阳能电池生产厂家提供的短路电流温度系数。
引入太阳能电池生产厂家提供的短路电流温度系数μI和开路电压温度系数μV得到修正的反向饱和电流的计算公式为:
其中,Voc_r是太阳能电池生产厂家提供的标准测试条件下太阳能电池的开路电压,μV是太阳能电池生产厂家提供的开路电压温度系数,太阳能电池理想因数n的数值由太阳能电池所使用的材料来决定。
最大功率点初始表达式获取方法包括:
当理想太阳能电池处于最大功率点(Vmp,Imp)时,根据最大功率理论得:
对输出特性表达式求关于理想太阳能电池输出电压V的导数并结合最大功率理论得到:
将该表达式与输出特性曲线结合得到最大功率点初始表达式为:
最大功率点电压:
最大功率点电流:
考虑电池温度、日照强度对最大功率点得影响,引入电池温度、日照强度与最大功率点之间的线性关系式对最大功率点初始表达式进行简化,其中电池温度、日照强度与最大功率点电压之间的线性关系式为:
Vmp_r是太阳能电池生产厂家提供的标准测试条件下的最大功率点电压;将该线性关系式代入到最大功率点初始表达式即可得到太阳能电池最大功率点电压、最大功率点电流、最大功率得求解表达式分别如下:
Pmp=Vmp·Imp
本发明的优点在于:
(1)本解析求解方法是建立在严格的理论推导基础之上,避开了循环迭代算法,变量关系明确,求解方便,克服了一般数值求解方法的缺陷。
(2)本解析求解方法能直观反映太阳电池温度、日照强度等因素对最大功率点参数的影响,实际工程应用也比较方便,而参数求解精度比一般的解析求解方法要高。
附图说明
下面对本发明说明书各幅附图表达的内容及图中的标记作简要说明:
图1为本发明理想太阳能电池最大功率点参数求解流程图;
图2为本发明理想太阳能电池等效电路模型。
具体实施方式
下面对照附图,通过对最优实施例的描述,对本发明的具体实施方式作进一步详细的说明。
图1所示为本发明一个优选实施例,包括如下实现步骤:
步骤S1、建立理想太阳能电池的数学模型和输出特性表达式。
太阳能电池的内部结构是一个由半导体材料制作的PN结,能将太阳能直接转换成直流电流,其非常类似于一个二极管。因此,若忽略电池材料的电阻特性以及损耗等因素,太阳能电池可以等效为一个电流源和一个二极管的并联,即为理想太阳能电池,如图2所示。
另外,单个太阳能电池片的电压只有0.5V左右,功率一般不超过2W,生产厂商通常采用Ns块单个电池片串联来获得足够的电压和功率。因此,由图2的数学模型,可以推导出理想太阳能电池的输出特性表达式为:
其中,I、V分别是理想太阳能电池的输出电流和输出电压,Iph是代表电流源的光生电流,Io是二极管的反向饱和电流,q是一个电子的电荷量(1.6×10-19C),k是玻尔兹曼常数(1.38×10-23J/K),T是太阳能电池的温度,n是太阳能电池的理想因数。
步骤S2、求解理想太阳能电池的光生电流Iph、反向饱和电流Io和理想因数n。
当太阳能电池开路(I=0)和短路(V=0)时,通过式(1)可以分别得到太阳能电池短路电流Isc和开路电压Voc的表达式
Isc=Iph 式(2)
太阳能电池的光生电流的大小取决于日照强度和电池的工作温度,其计算公式一般取为:
其中,S是日照强度,Sr是参考日照强度(1000W/m2),Tr是太阳能电池参考温度(25℃),Iph_r是太阳能电池生产厂家提供的标准测试条件(日照强度为1000W/m2,温度为25℃)下的短路电流,μI是太阳能电池生产厂家提供的短路电流温度系数。
当太阳能电池开路时,由式(1)可得到在参考温度Tr下,反向饱和电流Io_r的表达式为:
其中,Voc_r是太阳能电池生产厂家提供的标准测试条件下太阳能电池的开路电压。
而反向饱和电流Io由太阳电池PN结的结构和材料属性决定,并且受电池温度影响较大,和日照强度S基本没有关系。因此,在一般温度下,引入太阳能电池生产厂家提供的短路电流温度系数μI和开路电压温度系数μV,把短路电流和开路电压表示成电池温度T的线性函数。因此反向饱和电流Io关于电池温度T的计算公式为:
另外,太阳能电池理想因数n的数值由太阳能电池所使用的材料来决定,本求解方法中,n的数值按电池所使用的材料来取值:单晶硅材料n取1.2,多晶硅和AsGa材料n取1.3,CdTe和CIS材料n取1.5,非晶硅材料n取1.8。
步骤S3、求解理想太阳能电池的最大功率点电流Imp、最大功率点电压Vmp和最大功率Pmp
当理想太阳能电池处于最大功率点(Vmp,Imp)时,根据最大功率理论有等式(7)的成立:
对式(1)求关于理想太阳能电池输出电压V的导数,然后代入式(7),有等式(8)的成立:
结合式(1)和(8),可以推导出最大功率点电压Vmp和最大功率点电流Imp的表达式:
从式(9)可以看出,其属于超越方程,没有办法对Vmp直接求解,进而导致Imp和Pmp没有办法直接求解。目前Vmp较通用的解析求解公式是式(9)右边的用开路电压Voc代替,Isc+Io用短路电流Isc代替,则Vmp的表达式变为等式(11):
从式(11)可以看出,这是Vmp的直接求解解析表达式,其和通过中值定理推导出的Vmp的解析表达式完全一致,被很多文献采用作为Vmp的计算公式。虽然式(11)的求解比式(9)方便,但是用开路电压Voc代替以及用Isc代替Isc+Io后,使得Vmp的参数误差明显变大,使得式(11)在对参数求解精度要求较高的场合下不能得到应用。
对式(9)的另一种处理方法是引入ε(0<ε<1)因子,用ε因子来表征开路电压Voc与最大功率点电压Vmp之间的比例关系,从而使得Vmp的求解精度变高,相应的Vmp计算公式变为:
但是实际太阳能电池的最大功率点电压Vmp数值是和电池温度、日照强度有关联的,而式(12)右边的εVoc并不能反映这一点,所以式(12)也存在较大的参数误差。
另一方面,最大功率点电压Vmp和电池温度、日照强度之间的关系可以用等式(13)的近似线性估算公式来表示:
其中,Vmp_r是太阳能电池生产厂家提供的标准测试条件下的最大功率点电压。
式(13)能较准确的反映最大功率点电压Vmp随电池温度T和日照强度S的变化关系。因此,在本发明中,对式(9)右边的Vmp用式(13)的线性估算公式代替,从而可以得到理想太阳能电池最大功率点电压、最大功率点电流和最大功率的精确解析求解公式,分别如式(14)、(15)和(16)所示:
Pmp=Vmp·Imp 式(16)
式(14)、(15)和(16)即为理想太阳能电池最大功率点参数的解析求解公式,可以看出这里变量关系直接、求解不需要迭代算法,并且能够取得较高的参数求解精度。式(14)、(15)和(16)相较于公式(9)、(10)没有超越方程,计算更为方便,式(14)、(15)和(16)由于采用了式(13)的近似线性估算公式对公式(9)、(10)进行了变换而得到,虽然引入了近似估算公式(13)的部分误差,但是通过工程实际计算,接近于精确计算公式(9)、(10),实际误差小于直接采用公式(13)的误差,误差在工程计算的允许范围内,而且可以达到快速计算的目的,实际工程应用也比较方便,参数求解精度比一般的解析求解方法(11)式、(12)式要高。
显然本发明具体实现并不受上述方式的限制,只要采用了本发明的方法构思和技术方案进行的各种非实质性的改进,均在本发明的保护范围之内。

Claims (7)

1.一种理想太阳能电池最大功率点参数的求解方法,其特征在于:包括
建立理想太阳能电池的数学模型并根据数据模型求取输出特性表达式;
根据最大功率理论结合输出特性表达式推导出最大功率点的初始表达式;
采用电池温度、日照强度来修正太阳能电池最大功率点的初始表达式得到理想太阳能电池最大功率点的求解公式。
2.如权利要求1所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:理想太阳能电池模型为一个电流源和一个二极管并联。
3.如权利要求2所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:根据理想太阳能电池的数学模型求取的输出特性表达式为:
其中,I、V分别是理想太阳能电池的输出电流和输出电压,Iph是代表电流源的光生电流,Io是二极管的反向饱和电流,q是一个电子的电荷量,k是玻尔兹曼常数,Ns是串联电池片的数量,T是太阳能电池的温度,n是太阳能电池的理想因数。
4.如权利要求3所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:太阳能电池的光生电流Iph的计算公式为:
其中,S是日照强度,Sr是参考日照强度,Tr是太阳能电池参考温度,Iph_r是太阳能电池生产厂家提供的标准测试条件下太阳能电池的短路电流,μI是太阳能电池生产厂家提供的短路电流温度系数。
5.如权利要求3所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:引入太阳能电池生产厂家提供的短路电流温度系数μI和开路电压温度系数μV得到修正的反向饱和电流的计算公式为:
其中,Voc_r是太阳能电池生产厂家提供的标准测试条件下太阳能电池的开路电压,μV是太阳能电池生产厂家提供的开路电压温度系数,太阳能电池理想因数n的数值由太阳能电池所使用的材料来决定。
6.如权利要求1所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:最大功率点初始表达式获取方法包括:
当理想太阳能电池处于最大功率点(Vmp,Imp)时,根据最大功率理论得:
对输出特性表达式求关于理想太阳能电池输出电压V的导数并结合最大功率理论得到:
将该表达式与输出特性曲线结合得到最大功率点初始表达式为:
最大功率点电压:
最大功率点电流:
7.如权利要求1-6任一所述的一种理想太阳能电池最大功率点参数的求解方法,其特征在于:考虑电池温度、日照强度对最大功率点得影响,引入电池温度、日照强度与最大功率点之间的线性关系式对最大功率点初始表达式进行简化,其中电池温度、日照强度与最大功率点电压之间的线性关系式为:
Vmp_r是太阳能电池生产厂家提供的标准测试条件下的最大功率点电压;将该线性关系式代入到最大功率点初始表达式即可得到太阳能电池最大功率点电压、最大功率点电流、最大功率得求解表达式分别如下:
Pmp=Vmp·Imp
CN201810064916.1A 2018-01-23 2018-01-23 一种理想太阳能电池最大功率点参数的求解方法 Active CN108306617B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810064916.1A CN108306617B (zh) 2018-01-23 2018-01-23 一种理想太阳能电池最大功率点参数的求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810064916.1A CN108306617B (zh) 2018-01-23 2018-01-23 一种理想太阳能电池最大功率点参数的求解方法

Publications (2)

Publication Number Publication Date
CN108306617A true CN108306617A (zh) 2018-07-20
CN108306617B CN108306617B (zh) 2020-04-10

Family

ID=62866008

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810064916.1A Active CN108306617B (zh) 2018-01-23 2018-01-23 一种理想太阳能电池最大功率点参数的求解方法

Country Status (1)

Country Link
CN (1) CN108306617B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109800515A (zh) * 2019-01-24 2019-05-24 杭州电子科技大学 一种基于头脑风暴算法的太阳能电池参数求解方法
CN111769802A (zh) * 2020-05-09 2020-10-13 安徽工程大学 一种获取光伏电池输出特性曲线的方法和系统
CN112098853A (zh) * 2020-08-21 2020-12-18 北京车和家信息技术有限公司 一种容量衰减电池放电功率确定方法及装置
CN113346842A (zh) * 2021-08-02 2021-09-03 国网江西省电力有限公司电力科学研究院 一种太阳能电池串联电阻阻值的在线计算方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104320043A (zh) * 2014-09-24 2015-01-28 中国人民解放军空军工程设计研究局 一种光伏单元及其实现方法
US20160126847A1 (en) * 2014-10-31 2016-05-05 Majid Pahlevaninezhad Current sensorless mppt for pv micro-inverters
CN106249042A (zh) * 2016-07-08 2016-12-21 国网山东省电力公司东营供电公司 一种基于数据分析的光伏发电监测系统及方法
CN107423257A (zh) * 2017-03-31 2017-12-01 安徽工程大学 一种太阳能电池最大功率点参数求解方法
US9966899B2 (en) * 2012-10-16 2018-05-08 Volterra Semiconductor, LLC Systems and methods for controlling maximum power point tracking controllers

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9966899B2 (en) * 2012-10-16 2018-05-08 Volterra Semiconductor, LLC Systems and methods for controlling maximum power point tracking controllers
CN104320043A (zh) * 2014-09-24 2015-01-28 中国人民解放军空军工程设计研究局 一种光伏单元及其实现方法
US20160126847A1 (en) * 2014-10-31 2016-05-05 Majid Pahlevaninezhad Current sensorless mppt for pv micro-inverters
CN106249042A (zh) * 2016-07-08 2016-12-21 国网山东省电力公司东营供电公司 一种基于数据分析的光伏发电监测系统及方法
CN107423257A (zh) * 2017-03-31 2017-12-01 安徽工程大学 一种太阳能电池最大功率点参数求解方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘海洋: "基于改进黄金分割法的MPPT算法研究及实现", 《CNKI硕博士学位论文库》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109800515A (zh) * 2019-01-24 2019-05-24 杭州电子科技大学 一种基于头脑风暴算法的太阳能电池参数求解方法
CN111769802A (zh) * 2020-05-09 2020-10-13 安徽工程大学 一种获取光伏电池输出特性曲线的方法和系统
CN112098853A (zh) * 2020-08-21 2020-12-18 北京车和家信息技术有限公司 一种容量衰减电池放电功率确定方法及装置
CN113346842A (zh) * 2021-08-02 2021-09-03 国网江西省电力有限公司电力科学研究院 一种太阳能电池串联电阻阻值的在线计算方法及装置
CN113346842B (zh) * 2021-08-02 2021-12-31 国网江西省电力有限公司电力科学研究院 一种太阳能电池串联电阻阻值的在线计算方法及装置

Also Published As

Publication number Publication date
CN108306617B (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
CN103593577B (zh) 一种光伏发电系统输出功率建模及估算方法
CN103838291B (zh) 一种太阳能电池的最大功率点跟踪方法
Zhang et al. Modeling and prediction of PV module performance under different operating conditions based on power-law I–V model
CN106295068B (zh) 一种光伏组件双二极管模型的参数简化及提取方法
TWI461882B (zh) 太陽能模組系統之多點直接預測最大功率點追蹤方法及太陽能模組陣列之控制裝置
CN108306617A (zh) 一种理想太阳能电池最大功率点参数的求解方法
CN108694276B (zh) 一种计算串并联光伏组件输出特性的方法
TWI499886B (zh) 估算電路的最大功率點功率的方法
Abdulkadir et al. A new approach of modelling, simulation of mppt for photovoltaic system in simulink model
CN108491025B (zh) 一种离线系数拟合的光伏发电系统主动减载方法
CN102893264A (zh) 功率点跟踪
CN103105573A (zh) 一种基于朗伯w函数的太阳电池参数提取方法
Kumar et al. Effect of temperature, irradiation, humidity and wind on ideal/double diode PV system performance
Hassan et al. Investigating parasitic resistance of mesoporous-based solar cells with respect to thin-film and conventional solar cells
CN106951024B (zh) 一种太阳能组件快速mppt方法
CN207720094U (zh) 一种聚光光伏电池的测试装置
CN103454502A (zh) 一种测量任意光强和温度下光伏电池串联内阻的方法
CN105068590B (zh) 光伏发电功率控制方法和系统
CN107423257B (zh) 一种太阳能电池最大功率点参数求解方法
Hosseini Modeling and Simulation of Silicon Solar Cell in MATLAB/SIMULINK for Optimization
CN111446732A (zh) 一种光伏发电系统最大功率估计方法
CN104242819A (zh) 基于两种不同温度下电池电性参量预测光伏发电的方法
Shannan et al. Two diode model for parameters extraction of PV module
CN106970677A (zh) 一种太阳能逆变器快速mppt方法
Shen et al. Mathematical model of a solar module for energy yield simulation in photovoltaic systems

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