CN110986315B - 基于集中空调系统的室内温度动态调控方法 - Google Patents

基于集中空调系统的室内温度动态调控方法 Download PDF

Info

Publication number
CN110986315B
CN110986315B CN201911241162.3A CN201911241162A CN110986315B CN 110986315 B CN110986315 B CN 110986315B CN 201911241162 A CN201911241162 A CN 201911241162A CN 110986315 B CN110986315 B CN 110986315B
Authority
CN
China
Prior art keywords
temperature
time
stage
indoor
indoor temperature
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.)
Active
Application number
CN201911241162.3A
Other languages
English (en)
Other versions
CN110986315A (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.)
Xian University of Architecture and Technology
Original Assignee
Xian University of Architecture and Technology
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 Xian University of Architecture and Technology filed Critical Xian University of Architecture and Technology
Priority to CN201911241162.3A priority Critical patent/CN110986315B/zh
Publication of CN110986315A publication Critical patent/CN110986315A/zh
Application granted granted Critical
Publication of CN110986315B publication Critical patent/CN110986315B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/62Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
    • F24F11/63Electronic processing
    • F24F11/64Electronic processing using pre-stored data
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/50Control or safety arrangements characterised by user interfaces or communication
    • F24F11/61Control or safety arrangements characterised by user interfaces or communication using timers

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Human Computer Interaction (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

本发明属于建筑热环境领域,公开了一种基于集中空调系统的室内温度动态调控方法,在对人体热舒适主要影响因素进行分析的基础上,结合室外温度变化规律,提出空调环境动态温度调控方法,该方法在满足室内热环境舒适度要求的前提下,使室内热环境动态化,从而增加室内环境对人员的热刺激,在一定程度上减少了由于长期处于稳态空调环境中从而导致的建筑综合症的发生,由于本发明使得室内温度变化规律与室外温度变化趋势相似,因此可有效降低空调系统负荷,促进建筑节能。

Description

基于集中空调系统的室内温度动态调控方法
技术领域
本发明属于建筑热环境领域,具体涉及一种基于集中空调系统的室内温度动态调控方法。
背景技术
由于社会的不断发展,人们对于生活环境的关注与要求越来越高,空调成为人们生活中不可或缺的一部分。因此,室内环境是否舒适显然已成为空调室内环境的重要评价指标。同时,舒适健康的室内热环境也是人们健康生活和高效工作的前提。据相关调查显示,适当的热环境不仅可以带给人们精神上的愉悦,还可以将工作效率提高到18%以上。
为了研究人体的热舒适性,20世纪70年代丹麦学者Fanger教授通过分析大量的实验数据,建立了评估热环境舒适性的指标,即PMV-PPD指标。由于该指标综合考虑了影响人体热舒适的不同因素,目前被广泛使用。
PMV舒适性指标公式描述如下:
Figure BDA0002306267010000011
其中,
Figure BDA0002306267010000012
Figure BDA0002306267010000021
Figure BDA0002306267010000022
式中:
M——人体的新陈代谢率(W/m2),常用的新陈代谢率的取值范围为59~69;
W——人体对外做的机械功(W/m2),优选的人体对外做功的取值范围为0~1;
Pa——水蒸气分压力(Pa);
ta——室内空气温度(℃);
tr——平均辐射温度(℃);常用取值范围为25~35℃,优选的取值范围为25~30℃;
fcl——服装的表面积系数,由式-2计算得出;
Icl——服装热阻(m2·K/W);优选的取值范围为0.045~0.225;
tcl——服装的外表面温度(℃),由式-3计算得出;
hc——对流传热系数(W/(m2·℃)),由式-4计算得出;
Va——空气流速(m/s),常用的取值范围为0.5~3,优选的取值范围为2~3。
PMV与PPD之间的关系,如式-5所示。
Figure BDA0002306267010000023
当PMV在-0.5~+0.5之间,PPD小于等于10%(即允许有10%的人感到不满意)时,室内环境即为舒适环境。
由于传统空调力图使室内保持稳态的环境状况,即人一直处于稳态的环境中,主观热感觉处于热中性状态。研究表明,如果长期处于稳态空调环境中,由于缺少适当的热刺激,减弱了人类的热适应能力,最终导致人体体温调节能力和抵抗力下降,影响人体健康,出现各种与热环境有关的建筑关联症和病态建筑综合症,如目前常见的“空调病”等。针对该问题,国内外学者提出了动态热舒适调控策略。由于动态热舒适调控策略可有效利用热环境的小幅规律性波动,在满足热环境控制的同时实现舒适、健康和节能的统一,现已成为集中空调系统控制领域研究的热点。
然而,目前的相关研究在温度、风速以及舒适度的动态调节过程中,对于室外自然环境的变化规律及人体生理规律等因素考虑不足。因此,如何将室内热环境的控制过程与自然环境变化规律、人体生理规律等因素综合考虑,构建合理的室内环境热舒适调控方法是该领域尚需研究的问题。
发明内容
本发明的目的在于提供一种基于集中空调系统的室内温度动态调控方法,用以解决现有技术中长时间处于稳态空调环境中会相对减少热刺激,导致人体体温调节能力和抵抗力下降,影响人体健康的问题。
为了实现上述任务,本发明采用以下技术方案:
该调控方法是将空调系统的工作期间划分为多个温度调控周期进行调控,每个温度调控周期时长大于等于2h;
对于当前调控周期T,对T的开始整点时刻、结束整点时刻以及与中间转换时刻最接近的整点时刻所对应的室外温度t进行判断,并对室内温度进行升温或降温调控,所述中间转换时刻为R*T对应的时刻,R表示升降温时间比;
若t∈[t'min,t'max],则不对当前t所属时刻进行温度调控;其中:t'min≥27,t'max≤29;
若t∈(-∞,t'min)∪(t'max,+∞),则对当前t所属时刻的室内温度进行升温或降温调控,并且,
对于当前T的开始整点时刻,若当前T为首周期,则当前T的开始整点时刻的室内温度设定值等于室内温度映射值g(t),
Figure BDA0002306267010000041
tmax表示历史室外温度最大值,tmin表示历史室外温度最小值;若当前T不是首周期,则当前T的开始整点时刻的室内温度设定值等于上一个T的结束整点时刻室内温度设定值;
当前T的中间转换时刻的室内温度设定值等于与该中间转换时刻最邻近的整点时刻的g(t);
当前T的结束整点时刻的室内温度设定值根据式Ⅰ获得:
Figure BDA0002306267010000042
t1表示待调控日升温阶段中当前T的结束整点时刻的温度设定值,所述升温阶段为t持续增大的阶段,η1表示历史升温阶段的室内温度映射值变化率,t2表示待调控日降温阶段中每个T的结束整点时刻的温度设定值,所述降温阶段为t持续减小的阶段,η2表示历史降温阶段的室内温度映射值变化率。
进一步的,所述
Figure BDA0002306267010000043
进一步的,所述的温度调控周期
Figure BDA0002306267010000044
Δt表示温度变化量且Δt=t'max-t'min,cph表示待调控日的室内温度映射值变化率且cph<1.15。
更进一步的,所述cph=1。
更进一步的,所述温度调控周期T=2h。
进一步的,所述待调控日的升温阶段和降温阶段对应的整点时刻为8:00至20:00,总时长为12h。
本发明与现有技术相比具有以下技术特点:
本发明在对人体热舒适主要影响因素进行分析的基础上,结合室外温度变化规律,提出空调环境动态温度调控方法。该方法在满足室内热环境舒适度要求的前提下,使室内热环境动态化,从而增加室内环境对人员的热刺激,在一定程度上减少了由于长期处于稳态空调环境中从而导致的建筑综合症的发生,由于本发明使得室内温度变化规律与室外温度变化趋势相似,因此可有效降低空调系统负荷,促进建筑节能。
附图说明
图1是温度映射示意图;
图2是不同环境下热舒适温度范围图;
图3是室外温度与室内动态温度策略下的温度设定值变化趋势图;
图4是实施例中温度调控后的PMV指数示意图;
图5是实施例中温度调控后的PPD指数示意图。
具体实施方式
以下给出本发明的具体实施方式,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。首先对文中出现的专有名词进行解释:
升/降温阶段:室外温度连续升高的阶段为升温阶段,温度连续降低的阶段为降温阶段。
中间转换时刻:在每个温度调控周期内的温度发生上下波动的时刻,即升温与降温或降温与升温进行转换的时刻。
历史室外温度最大值:待调控日前一天或去年今日采集到的室外温度最大值。
历史室外温度最小值:待调控日前一天或去年今日采集到的室外温度最小值。
历史升温阶段的室内温度映射值变化率:待调控日前一天或去年今日计算得到的在升温阶段的各调控周期内室内温度映射值的变化率。
历史降温阶段的室内温度映射值变化率:待调控日前一天或去年今日计算得到的在降温阶段的各调控周期内室内温度映射值的变化率。
实施例1
在本实施例中公开了一种基于集中空调系统的室内温度动态调控方法,该调控方法是将空调系统的工作期间划分为多个温度调控周期进行调控,每个温度调控周期时长大于等于2h;
对于当前调控周期T,对T的开始整点时刻、结束整点时刻以及与中间转换时刻最接近的整点时刻所对应的室外温度t进行判断,并对室内温度进行升温或降温调控,所述中间转换时刻为R*T对应的时刻,R表示升降温时间比;
若t∈[t'min,t'max],则不对当前t所属时刻进行温度调控;其中:t'min≥27,t'max≤29;
若t∈(-∞,t'min)∪(t'max,+∞),则对当前t所属时刻的室内温度进行升温或降温调控,如当室外温度不在27℃~29℃的范围内时,需映射后调控,例如9点时,27+(2/10)*(31-29)=27.4,并且,
对于当前T的开始整点时刻,若当前T为首周期,则当前T的开始整点时刻的室内温度设定值等于室内温度映射值g(t),
Figure BDA0002306267010000061
tmax表示历史室外温度最大值,tmin表示历史室外温度最小值;若当前T不是首周期,则当前T的开始整点时刻的室内温度设定值等于上一个T的结束整点时刻室内温度设定值;
当前T的中间转换时刻的室内温度设定值等于与该中间转换时刻最邻近的整点时刻的g(t);
当前T的结束整点时刻的室内温度设定值根据式Ⅰ获得:
Figure BDA0002306267010000071
t1表示待调控日升温阶段中当前T的结束整点时刻的温度设定值,所述升温阶段为t持续增大的阶段,η1表示历史升温阶段的室内温度映射值变化率,t2表示待调控日降温阶段中每个T的结束整点时刻的温度设定值,所述降温阶段为t持续减小的阶段,η2表示历史降温阶段的室内温度映射值变化率。
具体的,所述升温时间和降温时间比例为10:7且
Figure BDA0002306267010000072
即在动态温度调控周期内升温时间约为70min,降温时间约为50min。因为调控周期为2个小时,即120分钟,则每个调控周期内升温时间为(120÷17)×10≈70.59分钟,降温时间为(120÷17)×7≈49.41分钟,则取每个调控周期内升温时间为70分钟,降温时间为50分钟。
具体的,所述待调控日的升温阶段和降温阶段对应的整点时刻为8:00至20:00,总时长为12h。
进一步的,所述的温度调控周期
Figure BDA0002306267010000073
Δt表示温度变化量且Δt=t'max-t'min,cph表示待调控日的室内温度映射值变化率且cph<1.15;
更进一步的,所述cph=1。
更进一步的,所述温度调控周期T=2h。
进一步的,g(t)表示的含义为:对将当前采集到的室外环境温度序列记为a1,a2,a3……a8,他们构成集合A,即A={a1,a2,a3……a8},V是由空调室内环境初始温度设定值序列(即映射值序列)v1,v2,v3……v8所组成的集合,即V={v1,v2,v3……v8}。设g为集合A到集合V上的一个映射,记作g:A→V。集合A与集合V都是有限集,如图1。
实施例2
以某市2016年8月15日夏季典型工况工作日为例,对早上8:00至20:00的室外环境温度数据进行采集,并依据上述步骤完成对室内初始温度设定值的计算。
表1某市2016年8月15日室外环境温度与室内温度映射值(单位:℃)
Figure BDA0002306267010000081
Figure BDA0002306267010000091
由表1可知,室内温度设定值的上界为29℃,下界为27.4℃,根据温度变化量与温度变化周期的关系可得,温度下降期间的温度变化率为0.24℃/h,温度上升期间的温度变化率为0.34℃/h,则在对应动态温度周期内的降温部分与升温部分,依据上述温度变化率进行计算与设定。
表2为基于动态温度方法的各周期内调控时刻的温度设定值(当房间温度低于最适温度时不做调整),由于10min内室外温度的变化量不大,所以对于±10min内的温度取值与最近时刻的温度值相同。
表2动态温度调控方法下各周期内调控时刻的温度设定值(单位:℃)
Figure BDA0002306267010000092
Figure BDA0002306267010000101
由于早晨8:00室外温度在是否进行动态温度调控的判断范围内,故对此时的温度值不做调整。由图3可以看出室内温度设定值在每个调控周期内以升温时间与降温时间为10:7的比例在变化,同时,整个工作日内的室内环境温度设定值变化规律与室外温度变化规律大致相似,8:00-15:00时温度处于逐渐上升的状态,在15:00-19:00时温度开始逐渐下降,因此温度调控周期的升温阶段为8:00-15:00,降温阶段为15:00-19:00。
室外升温时间段为8:00-15:00,共计7个小时,室内映射值上界为29,下界为27.4,则升温阶段的温度表化率η1=(29-27.4)÷7=0.34度/小时;室外降温时间段为15:00-20:00,共计5个小时,室内映射值上界为29,下界为27.8,则降温阶段的温度变化率η2=(29-27.8)÷5=0.24度/小时。
将表2中的温度序列带入,对不同周期的热舒适值与不满意率进行计算检验,计算结果如表3所示。
表3 PMV-PPD模型检验结果
Figure BDA0002306267010000102
Figure BDA0002306267010000111
由表3可知,采用上述动态温度调控方法,室内热环境满足人体的热舒适性要求,即-0.5<PMV<0.5,PPD<10%。

Claims (6)

1.基于集中空调系统的室内温度动态调控方法,其特征在于,该调控方法是将空调系统的工作期间划分为多个温度调控周期进行调控,每个温度调控周期时长大于等于2h;
对于当前调控周期T,对T的开始整点时刻、结束整点时刻以及与中间转换时刻最接近的整点时刻所对应的室外温度t进行判断,并对室内温度进行升温或降温调控,所述中间转换时刻为R*T对应的时刻,R表示升降温时间比;
若t∈[t'min,t'max],则不对当前t所属时刻进行温度调控;其中:t'min≥27,t'max≤29;
若t∈(-∞,t'min)∪(t'max,+∞),则对当前t所属时刻的室内温度进行升温或降温调控,并且,
对于当前T的开始整点时刻,若当前T为首周期,则当前T的开始整点时刻的室内温度设定值等于室内温度映射值g(t),
Figure FDA0002781544440000011
tmax表示历史室外温度最大值,tmin表示历史室外温度最小值;若当前T不是首周期,则当前T的开始整点时刻的室内温度设定值等于上一个T的结束整点时刻室内温度设定值;
当前T的中间转换时刻的室内温度设定值等于与该中间转换时刻最邻近的整点时刻的g(t);
当前T的结束整点时刻的室内温度设定值根据式Ⅰ获得:
Figure FDA0002781544440000012
t1表示待调控日升温阶段中当前T的结束整点时刻的温度设定值,所述升温阶段为t持续增大的阶段,η1表示历史升温阶段的室内温度映射值变化率,t2表示待调控日降温阶段中每个T的结束整点时刻的温度设定值,所述降温阶段为t持续减小的阶段,η2表示历史降温阶段的室内温度映射值变化率;
升/降温阶段:室外温度连续升高的阶段为升温阶段,温度连续降低的阶段为降温阶段;
中间转换时刻:在每个温度调控周期内的温度发生上下波动的时刻,即升温与降温或降温与升温进行转换的时刻;
历史室外温度最大值:待调控日前一天或去年今日采集到的室外温度最大值;
历史室外温度最小值:待调控日前一天或去年今日采集到的室外温度最小值;
历史升温阶段的室内温度映射值变化率:待调控日前一天或去年今日计算得到的在升温阶段的各调控周期内室内温度映射值的变化率;
历史降温阶段的室内温度映射值变化率:待调控日前一天或去年今日计算得到的在降温阶段的各调控周期内室内温度映射值的变化率。
2.如权利要求1所述的基于集中空调系统的室内温度动态调控方法,其特征在于,所述
Figure FDA0002781544440000021
3.如权利要求1所述的基于集中空调系统的室内温度动态调控方法,其特征在于,所述的温度调控周期
Figure FDA0002781544440000022
Δt表示温度变化量且Δt=t'max-t'min,cph表示待调控日的室内温度映射值变化率且cph<1.15。
4.如权利要求3所述的基于集中空调系统的室内温度动态调控方法,其特征在于,所述cph=1。
5.如权利要求3所述的基于集中空调系统的室内温度动态调控方法,其特征在于,所述温度调控周期T=2h。
6.如权利要求1所述的基于集中空调系统的室内温度动态调控方法,其特征在于,所述待调控日的升温阶段和降温阶段对应的整点时刻为8:00至20:00,总时长为12h。
CN201911241162.3A 2019-12-06 2019-12-06 基于集中空调系统的室内温度动态调控方法 Active CN110986315B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911241162.3A CN110986315B (zh) 2019-12-06 2019-12-06 基于集中空调系统的室内温度动态调控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911241162.3A CN110986315B (zh) 2019-12-06 2019-12-06 基于集中空调系统的室内温度动态调控方法

Publications (2)

Publication Number Publication Date
CN110986315A CN110986315A (zh) 2020-04-10
CN110986315B true CN110986315B (zh) 2021-03-16

Family

ID=70090727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911241162.3A Active CN110986315B (zh) 2019-12-06 2019-12-06 基于集中空调系统的室内温度动态调控方法

Country Status (1)

Country Link
CN (1) CN110986315B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113357828B (zh) * 2021-06-17 2022-10-28 青岛海尔空调器有限总公司 温度控制方法、装置和暖风机
CN116661528A (zh) * 2023-06-15 2023-08-29 中原工学院 基于热健康与生理节律的建筑室内环境调节装置与方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07253236A (ja) * 1994-03-17 1995-10-03 Toshiba Corp 空気調和機
CN104748318A (zh) * 2015-04-09 2015-07-01 广东美的制冷设备有限公司 空调器恒温除湿控制方法及装置、空调器
CN104896660A (zh) * 2015-05-20 2015-09-09 中南大学 一种办公建筑室内空调温度优化设定方法
CN105241035A (zh) * 2015-11-26 2016-01-13 湘潭大学 基于动态热舒适的空调器控制系统及其控制方法
CN107120721A (zh) * 2017-05-25 2017-09-01 河北健特建筑安装工程有限公司 一种集中供热动态气候补偿方法
CN107576001A (zh) * 2017-08-18 2018-01-12 广东美的暖通设备有限公司 预估能耗计算方法、系统及空调
CN110425691A (zh) * 2019-07-16 2019-11-08 江苏江荣智能科技有限公司 一种水冷中央空调系统的自适应节能控制系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07253236A (ja) * 1994-03-17 1995-10-03 Toshiba Corp 空気調和機
CN104748318A (zh) * 2015-04-09 2015-07-01 广东美的制冷设备有限公司 空调器恒温除湿控制方法及装置、空调器
CN104896660A (zh) * 2015-05-20 2015-09-09 中南大学 一种办公建筑室内空调温度优化设定方法
CN105241035A (zh) * 2015-11-26 2016-01-13 湘潭大学 基于动态热舒适的空调器控制系统及其控制方法
CN107120721A (zh) * 2017-05-25 2017-09-01 河北健特建筑安装工程有限公司 一种集中供热动态气候补偿方法
CN107120721B (zh) * 2017-05-25 2019-08-09 石家庄华浩能源科技有限公司 一种集中供热动态气候补偿方法
CN107576001A (zh) * 2017-08-18 2018-01-12 广东美的暖通设备有限公司 预估能耗计算方法、系统及空调
CN110425691A (zh) * 2019-07-16 2019-11-08 江苏江荣智能科技有限公司 一种水冷中央空调系统的自适应节能控制系统

Also Published As

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

Similar Documents

Publication Publication Date Title
Humphreys et al. An Adaptive Guideline for UK Office Temperatures: UK Adaptive Temperature Guideline
CN108954741A (zh) 酒店房间空调器舒适度控制方法
CN101949576B (zh) 节能空调和空调风扇联动系统
CN109579232A (zh) 提高空调睡眠模式舒适性的控制方法和系统
KR101162582B1 (ko) 습도 추정 장치 및 습도 추정 방법
CN110986315B (zh) 基于集中空调系统的室内温度动态调控方法
CN109855253B (zh) 用于空调器的控制方法
JP2007285579A (ja) 空調制御装置
CN110631211A (zh) 空调控制方法及装置
CN104566774A (zh) 一种空调器的控制方法
CN105157170A (zh) 变频空调控制方法
Su et al. Evaluation method of natural ventilation system based on thermal comfort in China
CN109523081A (zh) 一种地源热泵系统运行温度调节优化方法
CN106765986B (zh) 空调及其出风控制方法
WO2021169109A1 (zh) 空调温度控制方法、控制装置及空调
CN201859327U (zh) 节能空调和空调风扇联动系统
CN109138721B (zh) 一种基于空调运行的电动窗系统及其控制方法
JP2007120889A (ja) 空調制御装置
CN112880152B (zh) 基于用户状态的个性化空调系统动态温度调控方法及系统
JP2017003205A (ja) エアコンの省エネ冷房方法および冷房制御装置
CN116989379A (zh) 一种基于动态室温设定值和模型预测控制的室温控制方法
CN112923529B (zh) 基于热经历的供暖空调温度调控方法及系统
CN107860076B (zh) 基于人工智能的多用户动态调温中央空调系统及方法
CN116358114A (zh) 一种基于深度强化学习的空调温度控制方法
CN112923530B (zh) 基于人体热感觉的智能温控方法及系统

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