CN107122562B - 基于多fpga的有源配电网实时仿真器串行通讯方法 - Google Patents

基于多fpga的有源配电网实时仿真器串行通讯方法 Download PDF

Info

Publication number
CN107122562B
CN107122562B CN201710325399.4A CN201710325399A CN107122562B CN 107122562 B CN107122562 B CN 107122562B CN 201710325399 A CN201710325399 A CN 201710325399A CN 107122562 B CN107122562 B CN 107122562B
Authority
CN
China
Prior art keywords
simulation
time
fpga
real
power distribution
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
CN201710325399.4A
Other languages
English (en)
Other versions
CN107122562A (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.)
Tianjin University
State Grid Economic and Technological Research Institute
Original Assignee
Tianjin University
State Grid Economic and Technological Research Institute
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 Tianjin University, State Grid Economic and Technological Research Institute filed Critical Tianjin University
Priority to CN201710325399.4A priority Critical patent/CN107122562B/zh
Publication of CN107122562A publication Critical patent/CN107122562A/zh
Application granted granted Critical
Publication of CN107122562B publication Critical patent/CN107122562B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Design And Manufacture Of Integrated Circuits (AREA)

Abstract

一种基于多FPGA的有源配电网实时仿真器串行通讯方法，包括：1)将各子系统信息下载到对应FPGA中；2)设置仿真时刻t＝0，启动仿真；3)仿真时间向前推进一个步长，t＝t+Δt；4)各FPGA完成仿真计算；5)各FPGA发送计算得到的仿真接口数据；6)各FPGA接收并将仿真接口数据写入数据存储器中；7)对实时仿真器进行通讯结束校验；8)判断物理时间是否达到仿真时间t，如达到则进入下一步，否则实时仿真器待机至t后，进入下一步；9)判断仿真时间t是否达到设定的仿真终了时刻，如达到则仿真结束，否则返回步骤3)。本发明采用串行通讯的方法，在满足实时仿真器仿真速度要求的前提下，保证了仿真精度，为实现基于多FPGA的大规模有源配电网实时仿真奠定了基础。

Description

基于多FPGA的有源配电网实时仿真器串行通讯方法
技术领域
本发明涉及一种有源配电网实时仿真器通讯方法。特别是涉及一种基于多FPGA的有源配电网实时仿真器串行通讯方法。
背景技术
随着分布式电源、储能装置、微电网等配电侧各种资源的大量接入，有源配电网的组织结构和运行特性发生了广泛而深刻的变化。有源配电网的上述特点使其在规划设计、运行优化、保护控制、仿真分析等方面与传统配电系统相比存在较大差异。在仿真分析层面，有源配电系统中广泛接入的各种分布式电源、储能、电力电子装置等新型设备使得其自身的动态特性更加复杂，针对传统配电网的稳态仿真分析已不能满足需求，需要借助精细的暂态仿真来深入了解有源配电网的运行机理与动态特征。
有源配电网详细动态特性的分析与研究仅依靠离线暂态仿真是不够的，尤其是对各种控制器、保护装置、智能终端、新型能量管理系统等软硬件装置与系统的试验、测试均需要在硬件在环(hardware-in-the-loop，HIL)环境中完成，必须借助实时仿真器实现。目前，国外开发的商业实时仿真器有RTDS、ARENE、HYPERSIM、NETOMAC、RT-LAB等，这些实时仿真器全部采用DSP(Digital Signal Processor)、CPU(Central Processing Unit)、PowerPC等串行处理器作为底层硬件计算资源，通过多个处理器的并行计算，从而达到实时仿真的计算能力。
有源配电网复杂的网络结构和庞大的系统规模对实时仿真器的仿真精度、仿真速度、硬件资源等提出了新的挑战。在有源配电网中，电力电子开关具有高频动作特性，对该类元件的仿真需要较小的仿真步长；分布式电源及储能元件自身的控制、电力电子变流器的控制增加了系统的仿真规模，给硬件计算资源带来了较大的负担。基于串行处理器的实时仿真器囿于信号处理速度、物理结构的限制，实时仿真计算能力较为有限，同时，多个处理器之间数据的传输延时限制了仿真步长的选择。并行处理器现场可编程门阵列(FieldProgrammable Gate Array，FPGA)为高性能有源配电网实时仿真器的开发提供了一种新的思路。
FPGA是一款具有固有并行性质的可编程逻辑器件，可实现大量线程的并行处理，通过采用流水线的操作方式，极大地提高了数字信号的处理速度。同时，FPGA的I/O资源丰富，包括全双工LVDS通道、用户自定义I/O接口、高速收发器等，可实现大量数据的板级交互，使多FPGA联合实时仿真成为可能。
针对含大量电力电子以及非线性元件的有源配电网，要实现仿真步长在几个到十几个微秒之间的实时仿真，采用多个FPGA并行处理的方式提高底层硬件的计算能力，可以满足详细建模的有源配电网实时仿真的需求。将有源配电网划分为若干个子系统，分别在不同的FPGA上仿真，采用子系统之间的自然延迟时间补偿FPGA之间的通讯时间是合适的。考虑到子系统之间的自然延迟时间将会影响实时仿真器仿真精度，因此需要设计一种合适的数据通讯方式，以降低数据通讯时间对仿真器仿真精度的影响。
发明内容
本发明所要解决的技术问题是，提供一种能够满足高精度实时仿真的基于多FPGA的有源配电网实时仿真器串行通讯方法。
本发明所采用的技术方案是：一种基于多FPGA的有源配电网实时仿真器串行通讯方法，包括如下步骤：
1)在由N个FPGA构成的有源配电网实时仿真器的上位机中，将待仿真的有源配电系统依据拓扑连接关系及FPGA的计算资源划分为N个子系统，其中，N>1，读取各子系统元件的基本参数，形成各子系统电气部分的节点电导矩阵和控制部分的计算矩阵，设定实时仿真步长，将各子系统的相关信息分别下载到对应的FPGA中，根据子系统之间的连接关系及数据接口，设置第i个FPGA发送到与所述第i个FPGA直接相连的第j个FPGA的仿真接口数据的个数为Mi,j，其中i＝1,2,…,N，j＝1,2,…,N，仿真接口数据的传输延迟时间为L个时钟周期，各FPGA的仿真计算时间为Ki个时钟周期；
2)初始化实时仿真器，并设置仿真时刻t＝0，启动仿真；
3)仿真时间向前推进一个步长，t＝t+Δt；
4)每一个FPGA都从所述FPGA的串行通讯数据存储器中读出仿真所需的仿真接口数据，经过Ki个时钟周期完成步骤3)所述步长的仿真计算；
5)将每一个FPGA计算得到的Mi,j个仿真接口数据，发送到与所述FPGA直接相连的FPGA中；
6)每一个FPGA在接收到与所述FPGA直接相连的FPGA发送的仿真接口数据后，将仿真接口数据写入所述FPGA的串行通讯数据存储器中；
7)对实时仿真器进行仿真接口数据通讯结束校验，如果仿真接口数据通讯结束，则进入下一步，否则等待直至所有FPGA仿真接口数据通讯结束；
8)判断物理时间是否达到仿真时间t，如达到仿真时间t，则进入下一步，否则实时仿真器待机至仿真时间t后，进入下一步；
9)判断仿真时间t是否达到设定的仿真终了时刻T，如达到设定的仿真终了时刻T，则仿真结束，否则返回步骤3)。
步骤4)所述的串行通讯数据存储器，是由随机存取存储器RAMi构成，随机存取存储器RAMi的读使能信号ena_rdi在每一仿真时步仿真开始时刻为高电平，并持续
Figure BDA0001291039680000021
个时钟周期，随机存取存储器RAMi的读地址addr_rdi为连续整数
Figure BDA0001291039680000022
随机存取存储器RAMi的写使能信号ena_wri在每一仿真时步仿真计算完成后开始接收仿真接口数据时刻为高电平，并持续
Figure BDA0001291039680000023
个时钟周期，随机存取存储器RAMi的写地址addr_wri设置为连续整数
Figure BDA0001291039680000024
随机存取存储器RAMi的写使能信号ena_wri比随机存取存储器RAMi的读使能信号ena_rdi延迟L+Ki个时钟周期。
步骤7)所述的仿真接口数据通讯结束校验，是各FPGA将仿真接口数据写入串行通讯数据存储器结束后，分别生成通讯结束信号end_comm_sigi，通讯结束信号end_comm_sigi为高电平有效，当所有FPGA的通讯结束信号end_comm_sigi均为高电平时，仿真接口数据通讯结束，否则等待，直至仿真接口数据通讯结束。
本发明的基于多FPGA的有源配电网实时仿真器串行通讯方法，充分考虑了FPGA的自身硬件特性和有源配电网的结构特点，采用串行通讯的方法，在满足基于多FPGA的有源配电网实时仿真器仿真速度要求的前提下，保证了仿真器仿真精度，为实现基于多FPGA的大规模有源配电网实时仿真奠定了基础。
附图说明
图1是本发明基于多FPGA的有源配电网实时仿真器串行通讯方法的流程图；
图2是基于多FPGA的有源配电网实时仿真平台示意图；
图3是含光伏的有源配电网测试算例图；
图4是单极光伏发电单元详细结构图；
图5是光伏单元并网点A相电压Va仿真结果图；
图6是光伏单元并网点A相电流Ia仿真结果图；
图7是光伏单元有功功率Pinv仿真结果图。
具体实施方式
下面结合实施例和附图对本发明的基于多FPGA的有源配电网实时仿真器串行通讯方法做出详细说明。
如图1所示，本发明的基于多FPGA的有源配电网实时仿真器串行通讯方法，包括如下步骤：
1)在由N个FPGA构成的有源配电网实时仿真器的上位机中，将待仿真的有源配电系统依据拓扑连接关系及FPGA的计算资源划分为N个子系统，其中，N>1，读取各子系统元件的基本参数，形成各子系统电气部分的节点电导矩阵和控制部分的计算矩阵，设定实时仿真步长，将各子系统的相关信息分别下载到对应的FPGA中，根据子系统之间的连接关系及数据接口，设置第i个FPGA发送到与所述第i个FPGA直接相连的第j个FPGA的仿真接口数据的个数为Mi,j，其中i＝1,2,…,N，j＝1,2,…,N，仿真接口数据的传输延迟时间为L个时钟周期，各FPGA的仿真计算时间为Ki个时钟周期；
2)初始化实时仿真器，并设置仿真时刻t＝0，启动仿真；
3)仿真时间向前推进一个步长，t＝t+Δt；
4)每一个FPGA都从所述FPGA的串行通讯数据存储器中读出仿真所需的仿真接口数据，经过Ki个时钟周期完成步骤3)所述步长的仿真计算；
所述的串行通讯数据存储器，是由随机存取存储器RAMi构成，随机存取存储器RAMi的读使能信号ena_rdi在每一仿真时步仿真开始时刻为高电平，并持续
Figure BDA0001291039680000031
个时钟周期，随机存取存储器RAMi的读地址addr_rdi为连续整数
Figure BDA0001291039680000032
随机存取存储器RAMi的写使能信号ena_wri在每一仿真时步仿真计算完成后开始接收仿真接口数据时刻为高电平，并持续
Figure BDA0001291039680000033
个时钟周期，随机存取存储器RAMi的写地址addr_wri设置为连续整数
Figure BDA0001291039680000034
随机存取存储器RAMi的写使能信号ena_wri比随机存取存储器RAMi的读使能信号ena_rdi延迟L+Ki个时钟周期。
5)将每一个FPGA计算得到的Mi,j个仿真接口数据，发送到与所述FPGA直接相连的FPGA中；
6)每一个FPGA在接收到与所述FPGA直接相连的FPGA发送的仿真接口数据后，将仿真接口数据写入所述FPGA的串行通讯数据存储器中；
7)对实时仿真器进行仿真接口数据通讯结束校验，如果仿真接口数据通讯结束，则进入下一步，否则等待直至所有FPGA仿真接口数据通讯结束；
所述的仿真接口数据通讯结束校验，是各FPGA将仿真接口数据写入串行通讯数据存储器结束后，分别生成通讯结束信号end_comm_sigi，通讯结束信号end_comm_sigi为高电平有效，当所有FPGA的通讯结束信号end_comm_sigi均为高电平时，仿真接口数据通讯结束，否则等待，直至仿真接口数据通讯结束。
8)判断物理时间是否达到仿真时间t，如达到仿真时间t，则进入下一步，否则实时仿真器待机至仿真时间t后，进入下一步；
9)判断仿真时间t是否达到设定的仿真终了时刻T，如达到设定的仿真终了时刻T，则仿真结束，否则返回步骤3)。
下面给出具体实例：
本发明实施例中基于多FPGA的实时仿真器采用四块Altera公司的Stratix V系列FPGA5SGSMD5K2F40C2N及其配套官方开发板完成含光伏的有源配电网实时仿真。仿真平台如图2所示，各FPGA开发板之间采用光纤实现信号传输。整个实时仿真器通过125MHz的时钟驱动，发送通道并行时钟为125MHz，发送通道串行时钟为2500MHz，接收通道并行时钟为125MHz，接收通道串行时钟为2500MHz，FPGA之间单通道数据传输速率为2500Mbps。
测试算例为含光伏的有源配电网，如图3所示，在节点12处接入单级式光伏发电单元，光伏单元的详细结构如图4所示。光伏单元中光伏电池采用单二极管等效电路模拟，逆变器采用Vdc-Q控制，温度设置为298K，光伏电压参考值设为700V，无功功率参考值设为0var。仿真场景设置为1.2s时时光照强度从500W/m2增加到1000W/m2。
整个算例在多FPGA实时仿真器上进行仿真，其中网络部分占用FPGA1，光伏单元占用FPGA2，光伏单元与网络之间采用线路的贝瑞隆模型进行网络分割。FPGA1的仿真计算时间为2.208μs，通讯时间为0.488μs，FPGA2的仿真计算时间为1.792μs，通讯时间为0.488μs，整个仿真器仿真步长设置为3μs。
基于多FPGA的实时仿真器与商业软件PSCAD/EMTDC的仿真结果对比如图5～图7所示，图5是光伏单元并网点A相电压Va的仿真结果，图6是光伏单元并网点A相电流Ia的仿真结果，图7是光伏单元有功功率Pinv的仿真结果，PSCAD/EMTDC的仿真步长为3μs。从图中可以看出，两个仿真系统给出的结果基本一致，从而验证了本专利提出的基于多FPGA的有源配电网实时仿真器串行通讯方法的正确性。

Claims (2)

1.一种基于多FPGA的有源配电网实时仿真器串行通讯方法，其特征在于，包括如下步骤：
1)在由N个FPGA构成的有源配电网实时仿真器的上位机中，将待仿真的有源配电系统依据拓扑连接关系及FPGA的计算资源划分为N个子系统，其中，N>1，读取各子系统元件的基本参数，形成各子系统电气部分的节点电导矩阵和控制部分的计算矩阵，设定实时仿真步长，将各子系统的相关信息分别下载到对应的FPGA中，根据子系统之间的连接关系及数据接口，设置第i个FPGA发送到与所述第i个FPGA直接相连的第j个FPGA的仿真接口数据的个数为Mi,j，其中i＝1,2,…,N，j＝1,2,…,N，仿真接口数据的传输延迟时间为L个时钟周期，各FPGA的仿真计算时间为Ki个时钟周期；
2)初始化实时仿真器，并设置仿真时刻t＝0，启动仿真；
3)仿真时间向前推进一个步长，t＝t+Δt；
4)每一个FPGA都从所述FPGA的串行通讯数据存储器中读出仿真所需的仿真接口数据，经过Ki个时钟周期完成步骤3)所述步长的仿真计算；
所述的串行通讯数据存储器，是由随机存取存储器RAMi构成，随机存取存储器RAMi的读使能信号ena_rdi在每一仿真时步仿真开始时刻为高电平，并持续
Figure FDA0002324220160000011
个时钟周期，随机存取存储器RAMi的读地址addr_rdi为连续整数
Figure FDA0002324220160000012
随机存取存储器RAMi的写使能信号ena_wri在每一仿真时步仿真计算完成后开始接收仿真接口数据时刻为高电平，并持续
Figure FDA0002324220160000013
个时钟周期，随机存取存储器RAMi的写地址addr_wri设置为连续整数
Figure FDA0002324220160000014
随机存取存储器RAMi的写使能信号ena_wri比随机存取存储器RAMi的读使能信号ena_rdi延迟L+Ki个时钟周期；
5)将每一个FPGA计算得到的Mi,j个仿真接口数据，发送到与所述FPGA直接相连的FPGA中；
6)每一个FPGA在接收到与所述FPGA直接相连的FPGA发送的仿真接口数据后，将仿真接口数据写入所述FPGA的串行通讯数据存储器中；
7)对实时仿真器进行仿真接口数据通讯结束校验，如果仿真接口数据通讯结束，则进入下一步，否则等待直至所有FPGA仿真接口数据通讯结束；
8)判断物理时间是否达到仿真时间t，如达到仿真时间t，则进入下一步，否则实时仿真器待机至仿真时间t后，进入下一步；
9)判断仿真时间t是否达到设定的仿真终了时刻T，如达到设定的仿真终了时刻T，则仿真结束，否则返回步骤3)。
2.根据权利要求1所述的基于多FPGA的有源配电网实时仿真器串行通讯方法，其特征在于，步骤7)所述的仿真接口数据通讯结束校验，是各FPGA将仿真接口数据写入串行通讯数据存储器结束后，分别生成通讯结束信号end_comm_sigi，通讯结束信号end_comm_sigi为高电平有效，当所有FPGA的通讯结束信号end_comm_sigi均为高电平时，仿真接口数据通讯结束，否则等待，直至仿真接口数据通讯结束。
CN201710325399.4A 2017-05-10 2017-05-10 基于多fpga的有源配电网实时仿真器串行通讯方法 Expired - Fee Related CN107122562B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710325399.4A CN107122562B (zh) 2017-05-10 2017-05-10 基于多fpga的有源配电网实时仿真器串行通讯方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710325399.4A CN107122562B (zh) 2017-05-10 2017-05-10 基于多fpga的有源配电网实时仿真器串行通讯方法

Publications (2)

Publication Number Publication Date
CN107122562A CN107122562A (zh) 2017-09-01
CN107122562B true CN107122562B (zh) 2020-07-31

Family

ID=59726828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710325399.4A Expired - Fee Related CN107122562B (zh) 2017-05-10 2017-05-10 基于多fpga的有源配电网实时仿真器串行通讯方法

Country Status (1)

Country Link
CN (1) CN107122562B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108021741A (zh) * 2017-11-22 2018-05-11 天津大学 基于多fpga有源配电网多速率实时仿真器同步机制设计方法
CN109508479B (zh) * 2018-10-19 2022-12-16 天津大学 基于fpga的有源配电网实时仿真器参数配置通用化方法
CN110232227B (zh) * 2019-05-28 2022-11-04 天津大学 基于fpga的风电场实时仿真器超标量流水线设计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793562A (zh) * 2014-01-05 2014-05-14 天津大学 基于fpga的有源配电网暂态实时仿真系统设计方法
CN103942372A (zh) * 2014-04-04 2014-07-23 天津大学 基于fpga的有源配电网暂态实时仿真多速率接口方法
CN104715103A (zh) * 2015-01-12 2015-06-17 国家电网公司 基于fpga的光伏电池实时仿真模型设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8640070B2 (en) * 2010-11-08 2014-01-28 International Business Machines Corporation Method and infrastructure for cycle-reproducible simulation on large scale digital circuits on a coordinated set of field-programmable gate arrays (FPGAs)

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793562A (zh) * 2014-01-05 2014-05-14 天津大学 基于fpga的有源配电网暂态实时仿真系统设计方法
CN103942372A (zh) * 2014-04-04 2014-07-23 天津大学 基于fpga的有源配电网暂态实时仿真多速率接口方法
CN104715103A (zh) * 2015-01-12 2015-06-17 国家电网公司 基于fpga的光伏电池实时仿真模型设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于FPGA 的配电网暂态实时仿真研究(一)：功能模块实现;王成山 等;《中国电机工程学报》;20140105;第34卷(第1期);全文 *
基于FPGA 的配电网暂态实时仿真研究(二)：系统架构与算例验证;王成山 等;《中国电机工程学报》;20140205;第34卷(第4期);全文 *

Also Published As

Publication number Publication date
CN107122562A (zh) 2017-09-01

Similar Documents

Publication Publication Date Title
CN107423476B (zh) 基于多fpga的有源配电网实时仿真器并行通讯方法
CN103942372B (zh) 基于fpga的有源配电网暂态实时仿真多速率接口方法
CN111709107B (zh) 一种含大规模电力电子设备的电力网络仿真方法及系统
US11476667B2 (en) Hybrid electromagnetic transient simulation method for microgrid real-time simulation
CN106844256B (zh) 一种基于多fpga的有源配电网实时仿真器内部接口设计方法
CN102436534B (zh) 一种电力系统仿真方法
CN107122562B (zh) 基于多fpga的有源配电网实时仿真器串行通讯方法
Cao et al. Real‐time cyber− physical system co‐simulation testbed for microgrids control
CN108448631A (zh) 基于模型降阶的含分布式电源接入配电网的动态相量建模方法
CN106909732B (zh) 一种基于多fpga的有源配电网实时仿真器同步机制设计方法
CN107784158B (zh) 一种基于fpga的有源配电网实时仿真解算器的设计方法
Li et al. Synchronisation mechanism and interfaces design of multi‐FPGA‐based real‐time simulator for microgrids
CN113448890A (zh) 分布式仿真系统和控制方法
CN109583013B (zh) 基于fpga的有源配电网实时仿真器输出模块设计方法
CN101404040B (zh) 基于子图同构的对电力系统实时仿真的计算资源划分方法
CN108520105B (zh) 一种基于fpga的有源配电网多速率实时仿真方法
CN106845041A (zh) 基于mmc的实时仿真系统、仿真方法及mmc阀仿真器
CN110210106A (zh) 基于fpga的风电场实时仿真器模块级流水线设计方法
CN106844900B (zh) 电磁暂态仿真系统的搭设方法
CN109508479B (zh) 基于fpga的有源配电网实时仿真器参数配置通用化方法
CN110569558B (zh) 适用于微电网实时仿真的混合电磁暂态仿真方法
CN204856067U (zh) 功率在环型数字与物理混合实时仿真的试验系统
CN108563588B (zh) 基于fpga的有源配电网实时仿真器多速率接口设计方法
CN110489798A (zh) 面向有源配电网实时仿真的解算器细粒度高效设计方法
CN203551692U (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
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: 20200731