CN106089577A - 基于tmd及主动锚链结构控制的浮式风电机组减载方法 - Google Patents
基于tmd及主动锚链结构控制的浮式风电机组减载方法 Download PDFInfo
- Publication number
- CN106089577A CN106089577A CN201610594764.7A CN201610594764A CN106089577A CN 106089577 A CN106089577 A CN 106089577A CN 201610594764 A CN201610594764 A CN 201610594764A CN 106089577 A CN106089577 A CN 106089577A
- Authority
- CN
- China
- Prior art keywords
- theta
- floating platform
- anchor chain
- tmd
- centerdot
- 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
- 238000007667 floating Methods 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000006073 displacement reaction Methods 0.000 claims abstract description 15
- 230000033001 locomotion Effects 0.000 claims abstract description 7
- 230000001419 dependent effect Effects 0.000 claims abstract description 4
- 230000001133 acceleration Effects 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 6
- 230000008676 import Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 6
- 230000008859 change Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 230000003028 elevating effect Effects 0.000 abstract description 2
- 230000004044 response Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/331—Mechanical loads
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明公开了一种基于TMD及主动锚链结构控制的浮式风电机组减载方法,一)在锚链处于悬链状态:1)通过分析得到阻尼器力与活塞相对于缸体的位移之间的对应关系,2)分别建立浮台、塔架和TMD的运动方程;二)在锚链处于张紧状态:1)建立浮台及锚链的相关动力学方程,2)对控制锚链的电机进行主动控制,以收缩和释放连接浮台的锚链,从而使浮台减少纵荡,横荡和艏摇的现象。本发明通过采用被动TMD调节的同时,并针对不同海况改变悬链状态,对海上驳船式结构的漂浮式风力机组的载荷响应进行控制,从而可以有效降低浮式风力机在不同工况下的总体载荷和机舱、浮台的俯仰运动,进而提高漂浮式风电机组的输出电能质量并增加其使用寿命。
Description
技术领域
本发明专利涉及漂浮式风电机组减载控制方法,尤其涉及一种基于TMD及主动锚链结构的控制方法,它是一种可以有效降低浮式风电机组在不同海况下的总体载荷,在一定程度上保证浮台的平稳,进而提高风力机的使用寿命和输出电能质量的减载控制方法。
背景技术
漂浮式海上风力发电机组本身虽然能够较大程度利用深海的风能资源,但其外界载荷条件比陆地上固定式风力发电机组更加复杂,其原因是除了受到通常的风和浪载荷以外,还要考虑风、浪载荷的大小问题。同时相对固定式海上风力发电机来说,浮式平台的运动对系统发电也有很大的影响。由于风和浪具有随机性的特点,环境载荷波动会造成风电机组的载荷波动变化,使风电机组输出功率不稳定。特别是风和波浪引起的浮台和塔架的倾斜运行,由于大惯性和重力作用,大大增加风机的结构载荷,从而影响风力机的使用寿命和输出电能质量。
在载荷控制方面,目前更多的采用独立变桨距控制和在机舱中加大振动能量耗散方法对风电机组载荷进行控制研究。但采用以上方法对受到较大波浪载荷作用的浮式风电系统而言,其载荷控制无法满足要求,同时荷载控制也以牺牲部分功率为代价。
发明内容
有鉴于此,本发明的目的是提供一种基于TMD及主动锚链结构控制的浮式风电机组减载方法,以实现有效减小海上浮式风电机组的总体载荷,并在一定程度上保证浮台的平稳,提高风力机的使用寿命和输出电能质量。
本发明基于TMD及主动锚链结构控制的浮式风电机组减载方法,包括:
一)在锚链处于悬链状态
1)对浮台内被动TMD的结构进行分析,得到阻尼器力F与活塞相对于缸体的位移x之间的对应关系:式中,Ftmd为变阻尼力;dtmd为阻尼系数;为活塞相对缸体的速度;ktmd为弹簧系数;xtmd为弹簧伸缩位移;
2)利用Kane动力学方程,分别建立浮台、塔架和TMD的运动方程:
其中,θp表示浮台偏离垂直平面的角度,θt表示塔架偏离垂直平面的角度;kp表示浮台的弹性常数,kt表示塔架的弹性常数;dp表示浮台的阻尼系数,dt表示塔架的阻尼系数;Rtmd为到阻尼器质心到铰链的距离;Rt为塔架质心到铰链的距离;Rp为浮台质心到铰链的距离;Ip表示浮台转动惯量,It表示塔架转动惯量;g表示重力加速度;mtmd为阻尼器的质量;mt为塔架的质量;mp为浮台质量;xtmd表示阻尼器质心到铰链所在垂直面的距离;Ftmd为变阻尼力;
二)在锚链处于张紧状态
1)在锚链处于张紧状态时,建立海上漂浮式风机组浮台及锚链的相关动力学方程:
Tsinθ=wl=Tv,
Tcosθ=T′,
其中,T'是长为l的锚链在水深h处于张紧状态时受到的水平力,TV是长为l的锚链在水深h处于张紧状态时受到的垂向力,θ为链部上端切线方向与水平面的夹角,w为锚链在水中单位长度的重量;
在考虑浮台与锚链之间的关系时,平衡状态的情况如下:
式中,θi为第i根锚链与海底水平面的夹角;Ti为第i根锚链所受到的拉力;Ff为浮台所受浮力;Mz为漂浮式风机和浮台总质量;g为重力加速度;
2)将以上的平衡状态作为参考状态,在受到风和浪作用时,利用传感器实时检测浮台相对于参考状态的俯仰角加速度角速度和角位移θ,并将各量的检测值与期望值相比较;
通过比较后的结果制定电机的转向和转速大小的规则库;
利用上述规则库编写相应的程序,导入PLC控制器,从而对电机进行主动控制,以收缩和释放连接浮台的锚链,从而使浮台减少纵荡,横荡和艏摇的现象。
本发明的有益效果:
本发明一种基于TMD及主动锚链结构控制的浮式风电机组减载方法,通过采用被动TMD调节的同时,并针对不同海况改变悬链状态,对海上驳船式基础结构的漂浮式风力机组的载荷响应进行控制;针对实际中,当波浪和风作用引起浮台的仰俯运动较小时,通过在浮台内安装被动TMD阻尼器,并对其进行结构减振控制;当波浪和风作用引起浮台的仰俯运动很大时,通过改变悬链状态后,实时监测浮台的状态,利用拉紧后悬链的收放对浮台进行主动控制。从而可以有效降低浮式风力机在不同工况下的总体载荷和机舱、浮台的俯仰运动,进而提高漂浮式风电机组的输出电能质量并增加其使用寿命。
附图说明
图1为浮台中配置了被动TMD的漂浮式风力机结构示意图,图中:1-海平面,2-浮台,3-塔架,4-被动TMD;
图2为多阶段TMD及锚链控制示意图,图中海况1指锚链处于悬链状态,海况2指锚链处于张紧状态;
图3为被动TMD的力-相关参数结构图;
图4为张紧式锚泊线示意图;
图5为浮台受力示意图,图中:5-第一锚,6-第二锚,7-第三锚,8-第四锚;
图6为浮台受力分析俯视图,图中-第一锚链水平分力,-第二锚链水平分力,-第三锚链水平分力,-第四锚链水平分力;
图7为PLC控制流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步描述。
如图所示,本实施例基于TMD及主动锚链结构控制的浮式风电机组减载方法,包括:
一)在锚链处于悬链状态,此时为海况1一般海况下:
1)对浮台内被动TMD的结构进行分析,得到阻尼器力F与活塞相对于缸体的位移x之间的对应关系:式中,Ftmd为变阻尼力;dtmd为阻尼系数;为活塞相对缸体的速度;ktmd为弹簧系数;xtmd为弹簧伸缩位移;TMD即调谐质量阻尼器;
2)由于风电机组的浮台是承受整个风电机组的关键部分,如果浮台减少纵荡,横荡和艏摇的现象,那么风电机组的稳定性大大增加;因此,本发明只考虑浮台的俯仰模式。在海况1下,风电机组配置了被动TMD后,风机在原有模型中耦合了新的自由度,附图1给出了耦合纵向自由度的风机模型,耦合模型多了与被动TMD相关的驱动力和惯性力;
利用Kane动力学方程,分别建立浮台、塔架和TMD的运动方程:
其中,θp表示浮台偏离垂直平面的角度,θt表示塔架偏离垂直平面的角度;kp表示浮台的弹性常数,kt表示塔架的弹性常数;dp表示浮台的阻尼系数,dt表示塔架的阻尼系数;Rtmd为到阻尼器质心到铰链的距离;Rt为塔架质心到铰链的距离;Rp为浮台质心到铰链的距离;Ip表示浮台转动惯量,It表示塔架转动惯量;g表示重力加速度;xtmd表示阻尼器质心到铰链所在垂直面的距离;Ftmd为变阻尼力;mtmd为阻尼器的质量;mt为塔架的质量;mp为浮台质量;进而可得到附图2中选择海况1模式,在受到风和浪作用后,通过控制被动TMD产生一个与浮台振动方向相反的惯性力作用在浮台上,使浮台的振荡衰减并受到抑制,即能调节浮台的稳定性;
二)在锚链处于张紧状态,此时为海况2极端海况下:
先将整个风电机组做如附图4,5和6所示考虑,假设在没有风和浪作用时,作为最佳平衡状态(参考状态),受到风和浪影响时风电机组的状态都期望朝着平衡状态运动;
1)在锚链处于张紧状态时,先对平衡状态做数学分析,建立海上漂浮式风机组浮台及锚链的相关动力学方程:
Tsinθ=wl=Tv,
Tcosθ=T′,
其中,T'是长为l的锚链在水深h处于张紧状态时受到的水平力,TV是长为l的锚链在水深h处于张紧状态时受到的垂向力,θ为链部上端切线方向与水平面的夹角,w为锚链在水中单位长度的重量;本实施例以连接浮台的锚链为四根为例,一根锚链对应一个控制其收缩和释放的电机,在考虑浮台与锚链之间的关系时,平衡状态的情况如下:
式中,θi为第i根锚链与海底水平面的夹角;Ti为第i根锚链所受到的拉力;Ti′为第i根锚链水平分力;θi′为第i根锚链水平分力方向与浮台一边的夹角;Ff为浮台所受浮力;Mz为漂浮式风机和浮台总质量,i=1,2,3,4;
以上分析的量均为参考量,所有的不平衡状态都要将此状态作为参考状态;
2)电机的转向和大小与浮台俯仰角位移,角速度及角加速度均有相关的联系,故如附图5所示在风电机组受到风和浪作用时(考虑单一方向作用),浮台偏移参考位置,发生了俯仰运动,产生浮台俯仰角位移、角速度及角加速度;
利用传感器实时检测浮台相对于参考状态的俯仰角加速度角速度和角位移θ,并将各量的检测值与期望值相比较;
通过比较后的结果制定电机的转向和转速大小的规则库:
上述规则库中,通过对俯仰角位移的正负性,角速度以及角加速度变化方向确定,将各个变化量的组合分为6组不同的情况;
然后通过分析每种情况确定电机工作电压μ的方向和大小,从而控制电机的转向和大小,偏移量方程如下所示:
其中为位移观测量,θi为位移参考量;
电机工作电压的表达式如下:
其中,μi表示第i个电机的工作电压;为浮台俯仰角位移偏移量;为浮台俯仰角速度;为浮台俯仰角加速度;ai,bi和ci分别角位移偏移量,角速度,角加速度相关系数;ki为放大系数;
利用上述规则库编写相应的程序,导入PLC控制器,从而对控制锚链的电机的工作电压实时控制,进而使电机转向和大小得到确定,以收缩或释放连接浮台的锚链,从而使浮台减少纵荡,横荡和艏摇的现象。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种基于TMD及主动锚链结构控制的浮式风电机组减载方法,其特征在于:包括:
一)在锚链处于悬链状态
1)对浮台内被动TMD的结构进行分析,得到阻尼器力F与活塞相对于缸体的位移x之间的对应关系:式中,Ftmd为变阻尼力;dtmd为阻尼系数;为活塞相对缸体的速度;ktmd为弹簧系数;xtmd为弹簧伸缩位移;
2)利用Kane动力学方程,分别建立浮台、塔架和TMD的运动方程:
其中,θp表示浮台偏离垂直平面的角度,θt表示塔架偏离垂直平面的角度;kp表示浮台的弹性常数,kt表示塔架的弹性常数;dp表示浮台的阻尼系数,dt表示塔架的阻尼系数;Rtmd为到阻尼器质心到铰链的距离;Rt为塔架质心到铰链的距离;Rp为浮台质心到铰链的距离;Ip表示浮台转动惯量,It表示塔架转动惯量;g表示重力加速度;xtmd表示阻尼器质心到铰链所在垂直面的距离;Ftmd为变阻尼力;mtmd为阻尼器的质量;mt为塔架的质量;mp为浮台质量;
二)在锚链处于张紧状态
1)在锚链处于张紧状态时,建立海上漂浮式风机组浮台及锚链的相关动力学方程:
T sinθ=wl=Tv,
T cosθ=T′,
其中,T′是长为l的锚链在水深h处于张紧状态时受到的水平力,TV是长为l的锚链在水深h处于张紧状态时受到的垂向力,θ为链部上端切线方向与水平面的夹角,w为锚链在水中单位长度的重量;
在考虑浮台与锚链之间的关系时,平衡状态的情况如下:
式中,θi为第i根锚链与海底水平面的夹角;Ti为第i根锚链所受到的拉力;Ff为浮台所受浮力;Mz为漂浮式风机和浮台的质量;g为重力加速度;
2)将以上的平衡状态作为参考状态,在受到风和浪作用时,利用传感器实时检测浮台相对于参考状态的俯仰角加速度角速度和角位移θ,并将各量的检测值与期望值相比较;
通过比较后的结果制定电机的转向和转速大小的规则库;
利用规则库编写相应的程序,导入PLC控制器,从而对控制锚链的电机进行主动控制,以收缩和释放连接浮台的锚链,从而使浮台减少纵荡,横荡和艏摇的现象。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610594764.7A CN106089577B (zh) | 2016-07-26 | 2016-07-26 | 基于tmd及主动锚链结构控制的浮式风电机组减载方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610594764.7A CN106089577B (zh) | 2016-07-26 | 2016-07-26 | 基于tmd及主动锚链结构控制的浮式风电机组减载方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106089577A true CN106089577A (zh) | 2016-11-09 |
CN106089577B CN106089577B (zh) | 2018-07-10 |
Family
ID=57450319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610594764.7A Active CN106089577B (zh) | 2016-07-26 | 2016-07-26 | 基于tmd及主动锚链结构控制的浮式风电机组减载方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106089577B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112682256A (zh) * | 2020-12-09 | 2021-04-20 | 北京华能新锐控制技术有限公司 | 一种基于tmd与变桨优化控制的风机联合减载方法 |
CN114162263A (zh) * | 2021-12-17 | 2022-03-11 | 浙江大学 | 一种基于主动控制的漂浮式风力机系泊系统和控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291546A (zh) * | 2013-06-27 | 2013-09-11 | 上海交通大学 | 调液柱阻尼控制的近海漂浮式风电机 |
CN103541861A (zh) * | 2013-10-30 | 2014-01-29 | 新疆金风科技股份有限公司 | 浮动式风电机组塔架负阻尼抑制系统和方法 |
CN104986301A (zh) * | 2015-06-29 | 2015-10-21 | 武汉理工大学 | 一种组合式漂浮风力发电平台 |
WO2015158348A1 (en) * | 2014-04-14 | 2015-10-22 | Envision Energy (Denmark) Aps | Wind turbine with floating foundation and position regulating control system and method thereof |
-
2016
- 2016-07-26 CN CN201610594764.7A patent/CN106089577B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291546A (zh) * | 2013-06-27 | 2013-09-11 | 上海交通大学 | 调液柱阻尼控制的近海漂浮式风电机 |
CN103541861A (zh) * | 2013-10-30 | 2014-01-29 | 新疆金风科技股份有限公司 | 浮动式风电机组塔架负阻尼抑制系统和方法 |
WO2015158348A1 (en) * | 2014-04-14 | 2015-10-22 | Envision Energy (Denmark) Aps | Wind turbine with floating foundation and position regulating control system and method thereof |
CN104986301A (zh) * | 2015-06-29 | 2015-10-21 | 武汉理工大学 | 一种组合式漂浮风力发电平台 |
Non-Patent Citations (1)
Title |
---|
贺尔铭等: "基于TMD的海上浮动风力机结构振动控制研究", 《西北工业大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112682256A (zh) * | 2020-12-09 | 2021-04-20 | 北京华能新锐控制技术有限公司 | 一种基于tmd与变桨优化控制的风机联合减载方法 |
CN114162263A (zh) * | 2021-12-17 | 2022-03-11 | 浙江大学 | 一种基于主动控制的漂浮式风力机系泊系统和控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106089577B (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Si et al. | The influence of power-take-off control on the dynamic response and power output of combined semi-submersible floating wind turbine and point-absorber wave energy converters | |
Yang et al. | Coupled modeling and structural vibration control for floating offshore wind turbine | |
Stewart et al. | Offshore wind turbine load reduction employing optimal passive tuned mass damping systems | |
Cermelli et al. | WindFloat: a floating foundation for offshore wind turbines—part II: hydrodynamics analysis | |
CN103708004B (zh) | 减摇装置、浮动基础和海上风力发电机 | |
Nielsen et al. | Integrated dynamic analysis of floating offshore wind turbines | |
CN109441733B (zh) | 汲能-减振深海风力发电浮式半潜平台 | |
Bagbanci et al. | Review of offshore floating wind turbines concepts | |
Ormberg et al. | Global analysis of a floating wind turbine using an aero-hydro-elastic model: Part 1—code development and case study | |
JP2017515033A (ja) | 浮体式基礎及び位置調整制御システムを有する風力タービン並びにその位置調整制御方法 | |
Bagbanci | Dynamic analysis of offshore floating wind turbines | |
CN110182326A (zh) | 一种具有减隔振功能的半潜式风力发电机平台 | |
Crozier | Design and dynamic modeling of the support structure for a 10 MW offshore wind turbine | |
Battezzato et al. | Performance assessment of a 2 DOF gyroscopic wave energy converter | |
Zhao et al. | Experimental study on new multi-column tension-leg-type floating wind turbine | |
Dinh et al. | On the modeling of spar-type floating offshore wind turbines | |
Jahangiri et al. | Multibody dynamics of a floating wind turbine considering the flexibility between nacelle and tower | |
CN106089577A (zh) | 基于tmd及主动锚链结构控制的浮式风电机组减载方法 | |
Karimirad et al. | Effect of aerodynamic and hydrodynamic damping on dynamic response of a spar type floating wind turbine | |
Meng et al. | Research on dynamic response characteristics of 6MW spar-type floating offshore wind turbine | |
Chen et al. | Study on gyroscopic effect of floating offshore wind turbines | |
Tang et al. | Experiment study of dynamics response for wind turbine system of floating foundation | |
Nallayarasu et al. | Experimental and numerical investigation on hydrodynamic response of spar with wind turbine under regular waves | |
Ramachandran | A numerical model for a floating TLP wind turbine | |
Nallayarasu et al. | Experimental and numerical investigation of hydrodynamic response of spar with wind turbine under random waves |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |