CN111950868B - Comprehensive energy system load scene generation method based on generation countermeasure network - Google Patents

Comprehensive energy system load scene generation method based on generation countermeasure network Download PDF

Info

Publication number
CN111950868B
CN111950868B CN202010735071.1A CN202010735071A CN111950868B CN 111950868 B CN111950868 B CN 111950868B CN 202010735071 A CN202010735071 A CN 202010735071A CN 111950868 B CN111950868 B CN 111950868B
Authority
CN
China
Prior art keywords
scene
load
model
network
energy system
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
CN202010735071.1A
Other languages
Chinese (zh)
Other versions
CN111950868A (en
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.)
State Grid Shandong Electric Power Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
State Grid Electric Power Research Institute
Original Assignee
State Grid Shandong Electric Power Co Ltd
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
State Grid Electric Power 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 State Grid Shandong Electric Power Co Ltd, Nari Technology Co Ltd, NARI Nanjing Control System Co Ltd, State Grid Electric Power Research Institute filed Critical State Grid Shandong Electric Power Co Ltd
Priority to CN202010735071.1A priority Critical patent/CN111950868B/en
Publication of CN111950868A publication Critical patent/CN111950868A/en
Application granted granted Critical
Publication of CN111950868B publication Critical patent/CN111950868B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Economics (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Feedback Control In General (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a comprehensive energy system load scene generation method based on a generation countermeasure network, which comprises the following steps: s1: acquiring load data of a comprehensive energy system, and establishing a cold, heat and power load sample data set; s2: constructing and generating a confrontation network model; s3: alternately training a generator network and a discriminator network D; s4: and generating a comprehensive energy system load scene. The method utilizes the implicit relation of the generated confrontation network learning sample data to capture the internal complex deep nonlinear structural features. And when the real data is lacked or the diversity of the real data is insufficient, generating a comprehensive energy system multi-load scene with the statistical characteristics similar to the real scene. The scene generated by the method of the invention can provide decision support for load prediction, abnormal detection and operation scheduling, and the flexibility and reliability of the system are improved.

Description

Comprehensive energy system load scene generation method based on generation countermeasure network
Technical Field
The invention relates to a comprehensive energy system load scene generation method based on a generation countermeasure network, and belongs to the technical field of operation and analysis of comprehensive energy systems.
Background
In recent years, the degree of informatization and digitization has been increasing, and mining technology based on large data has become mature, and the value of data has been recognized and utilized. The comprehensive energy data with multi-source heterogeneity, complex relevance and real-time interactivity has huge information value, value information is extracted from the energy big data, a large number of scenes are calculated and analyzed, the change characteristics of loads are reflected, and decision support can be provided for key problems such as fault detection, load prediction, safety assessment and energy management. However, such data is often not readily available due to limitations in data privacy, data security, acquisition costs, etc., which prevents large-scale availability of data, posing challenges to further deployment of data-driven technologies in the energy field.
Therefore, it is necessary to provide a method for generating a comprehensive energy system scene, which finds and learns the potential association and the complex coupling in the training data when the real data is absent or the diversity of the real data is insufficient, mines the characteristics of the data, directly learns the probability distribution of the load time series data, and automatically generates a comprehensive energy system multi-load scene with the statistical characteristics similar to that of the real scene without explicitly modeling the probability distribution.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a comprehensive energy system load scene generation method based on a generation countermeasure network, which is used for generating a comprehensive energy system multi-element load scene with similar statistical characteristics to a real scene when real data is lacked or the diversity of the real data is insufficient.
In order to achieve the above object, the present invention provides a method for generating a comprehensive energy system load scenario based on a generation countermeasure network, comprising the following steps:
s1: acquiring load data of a comprehensive energy system, and establishing a cold, heat and power load sample data set;
s2: constructing and generating a confrontation network model;
s3: alternately training a generator network and a discriminator network D;
s4: and generating a comprehensive energy system load scene.
Preferably, the step S1 comprises the steps of:
s11: acquiring cold, heat and electricity load data of each day, acquiring N load values in one day, and acquiring the load values at the same interval;
s12: preprocessing the acquired load value, and completing lost data or correcting abnormal data;
s13: the load values are normalized according to the following equation, scaling the value of each datum to a value between 0 and 1:
Figure BDA0002604697750000021
in the formula: x denotes a sample data set, x = [ x = [ x ] 1 ,x 2 ,…,x n ]N is the total number of samples; x is the number of min Representing a minimum value in the sample data set; x is the number of max Representing a maximum value in the sample data set; x is the number of i Is the ith data in the sample data set, and the value range of i is [1, n ]]。
Preferably, the step S2 comprises the steps of:
s21: defining a generator model G consisting of a full-connection layer and four deconvolution layers, performing batch normalization processing after each calculation layer, and finally outputting by using a Sigmoid activation function;
wherein, the Sigmoid activation function is:
Figure BDA0002604697750000022
s22: defining a discriminator model D consisting of three convolution layers and a full-connection layer, wherein except the convolution layer of the first layer and the convolution layer of the third layer, the convolution layer of the second layer and the full-connection layer are subjected to batch normalization processing, and finally LeakyRelu is used as an activation function to be output;
wherein, leakyRelu activation function is:
Figure BDA0002604697750000023
preferably, the input to the generator model G is a set of random arrays of size 96 × 1.
Preferably, the filter sizes of the first, second and third convolutional layers are all 5 × 5.
Preferably, the step S3 comprises the steps of:
s31: the noise generating unit N generates a random noise sequence z of a certain length i
S32: random noise sequence z i Inputting a generator model G, and mapping the generator model G into a scene G (z) through a network of the generator model G;
s33: scene G (z) generated by generator model G and sampled training sample scene x i Inputting the values into a discriminator network D, wherein the discriminator network D respectively provides values D (x) and D (G (z)) representing the scene as true and false;
s34: respectively calculating loss functions of the discriminator model D and the generator model G according to the following formula;
L G =-E Z [D(G(z))],
L D =-E X [D(x)+E Z [D(G(z))]],
wherein L is G Representing the loss function of the generator, L D Representing the loss function of the arbiter, E Z Expected value, E, representing the corresponding random noise distribution X An expected value representing a corresponding true sample distribution;
s35: and reversely transmitting the result of the S34 to the generator model G and the discriminator model, and updating and adjusting the network parameter theta according to the following formula g And a network parameter theta d
Figure BDA0002604697750000031
Figure BDA0002604697750000032
Wherein, theta g,t-1 And theta d,t-1 The parameters of the generator network and the parameters of the discriminator network D are generated before the t-th round of training begins;
s36: repeating steps S31-S35 until the network parameter theta g And a network parameter theta d Until convergence.
Preferably, in step S4, the generator model G and the discriminator model D after training are fixed, and the generator model G is used for the random noise sequence z input from the input i And simulating to generate and output phase cooling heat and power load data.
The invention achieves the following beneficial effects:
the invention discloses a comprehensive energy system load scene generation method based on generation of a countermeasure network, which utilizes the implicit relation of the generated countermeasure network learning sample data to capture the internal complex deep nonlinear structure characteristics and generate new data with the statistical distribution characteristics similar to the real data. When the real data are lacked or the diversity of the real data is insufficient, the comprehensive energy system multi-load scene with the statistical characteristics similar to the real scene is automatically generated, explicit modeling on probability distribution is not needed, decision support is provided for abnormal detection and operation scheduling, the flexibility and the reliability of the system are improved, and the problem that the diversity of the real data is lacked or the diversity of the real data is insufficient is solved.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a graph of loss function variation for the generator and arbiter training process of the present invention;
FIG. 3 is a comparison of generated scene and real scene sample data for the present invention;
fig. 4 is a comparison graph of probability distribution characteristics of generated scene and real scene data of the present invention.
Detailed Description
The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
A comprehensive energy system load scene generation method based on generation of a countermeasure network comprises the following steps:
s1: acquiring load data of a comprehensive energy system, and establishing a cold, heat and power load sample data set;
s2: constructing and generating a confrontation network model;
s3: alternately training a generator network and a discriminator network D;
s4: and generating a comprehensive energy system load scene.
Further, the step S1 includes the steps of:
s11: acquiring cold, heat and electricity load data of each day, acquiring N load values in one day, and acquiring the load values at the same interval;
s12: preprocessing the acquired load value, and completing lost data or correcting abnormal data;
s13: the load values are normalized, scaling the value of each datum to a value between 0 and 1, according to the following equation:
Figure BDA0002604697750000041
in the formula: x represents a sample data set, x = [ x ] 1 ,x 2 ,…,x n ]N is the total number of samples; x is the number of min Representing a minimum value in the sample data set; x is the number of max Representing a maximum value in the sample data set; x is a radical of a fluorine atom i Is the ith data in the sample data set, and the value range of i is [1, n ]]。
Further, the step S2 includes the steps of:
s21: defining a generator model G consisting of a full-connection layer and four deconvolution layers, performing batch normalization processing after each calculation layer, and finally outputting by using a Sigmoid activation function;
wherein, the Sigmoid activation function is:
Figure BDA0002604697750000042
s22: defining a discriminator model D consisting of three convolution layers and a full-link layer, wherein except the first convolution layer and the third convolution layer, the second convolution layer and the full-link layer are subjected to batch normalization processing, and finally LeakyRelu is used as an activation function to be output;
wherein, the LeakyRelu activation function is:
Figure BDA0002604697750000043
further, the input to the generator model G is a set of random arrays of size 96 × 1.
Furthermore, the filter sizes of the first layer of convolutional layer, the second layer of convolutional layer and the third layer of convolutional layer are all 5 × 5.
Further, the step S3 includes the steps of:
s31: the noise generation unit N generates a random noise sequence z of a certain length i
S32: random noise sequence z i Inputting a generator model G and mapping the generator model G into a scene G (z) through a network of the generator model G;
s33: scene G (z) generated by generator model G and sampled training sample scene x i Inputting the values into a discriminator network D, wherein the discriminator network D respectively provides values D (x) and D (G (z)) representing the scene as true and false;
s34: respectively calculating loss functions of the discriminator model G and the generator model G according to the following formula;
L G =-E Z [D(G(z))],
L D =-E X [D(x)+E Z [D(G(z))]],
wherein L is G Representing the loss function of the generator, L D Representing the loss function of the arbiter, E Z Expected value, E, representing the corresponding random noise distribution X An expected value representing a corresponding true sample distribution;
s35: the result of S34 is reversely propagated to the generator model G and the discriminator model, and the network parameter theta is updated and adjusted according to the following formula g And a network parameter theta d
Figure BDA0002604697750000051
Figure BDA0002604697750000052
Wherein, theta g,t-1 And theta d,t-1 The parameters of the generator network and the parameters of the discriminator network D are generated before the t round of training begins;
s36: repeating steps S31-S35 until the network parameter theta g And a network parameter theta d Until convergence.
Further, in step S4, the generator model G and the discriminator model D after training are fixed, and the generator model G is used for the random noise sequence z input from the input i And simulating to generate and output phase cooling heat and power load data.
Example (b):
the embodiment uses the load data of residents in a certain area from 2013 to 2017 to generate a typical scene.
The generator and arbiter training process loss function variation graph of this embodiment is shown in fig. 2. After 2000 times of iterative training, the loss functions of the generator and the discriminator are basically converged, which shows that the generated countermeasure network reaches the balance at the moment. The loss function of the generator steadily fluctuates around the value 0 after continuously decreasing, which indicates that the artificial sample distribution synthesized by the generator is very close to the real sample distribution.
And selecting partial curves from the test set to verify the generated scene. A comparison graph of sample data of a generated scene and a real scene is shown in fig. 3, a curve generated by the model has an output trend consistent with that of a real sample, and is in accordance with an actual load characteristic but not completely consistent with an original sample, and a consistent fluctuation rule is maintained between the two. This means that the model does not over-fit the data in the training set, but learns the inherent relevance of the sample data, and then generates completely new data, which indicates that the model has strong generalization capability.
And analyzing the probability distribution characteristics of the generated scene by adopting a Cumulative Distribution Function (CDF). The comparison graph of the probability distribution characteristics of the generated scene data and the real scene data of the invention is shown in FIG. 4. The CDF curve of the generated load scene set is highly fitted, the probability distribution characteristics of the generated scene are very similar to those of the historical scene, and the GAN model is proved to have the learning capability of learning the potential probability distribution of given data.
According to the comprehensive energy system load scene generation method, the comprehensive energy system load scene is generated by adopting the generation countermeasure network, the fluctuation and probability distribution characteristics of the load under different conditions can be effectively represented, and the load scene can be directly generated, so that the model has strong generalization capability. The samples are made to exhibit diversity of mode changes while ensuring that the generated data has a high degree of similarity to the actual scene.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.

Claims (4)

1. A comprehensive energy system load scene generation method based on generation countermeasure network is characterized by comprising the following steps:
s1: acquiring load data of a comprehensive energy system, and establishing a cold, heat and power load sample data set;
s2: constructing and generating a confrontation network model;
s3: alternately training a generator model and a discriminator model D;
s4: generating a comprehensive energy system load scene;
the step S1 includes the steps of:
s11: acquiring cold, heat and electricity load data of each day, acquiring N load values in one day, and acquiring the load values at the same interval;
s12: preprocessing the acquired load value, and completing lost data or correcting abnormal data;
s13: the load values are normalized according to the following equation, scaling the value of each datum to a value between 0 and 1:
Figure FDA0003767206580000011
in the formula: x represents a sample data set, x = [ x ] 1 ,x 2 ,…,x i ,…,x n ]N is the total number of samples; x is the number of min Representing a minimum value in the sample data set; x is the number of max Representing a maximum value in the sample data set; x is a radical of a fluorine atom i Is the ith data in the sample data set, and the value range of i is [1, n ]];
The step S2 includes the steps of:
s21: defining a generator model G consisting of a full-connection layer and four deconvolution layers, performing batch normalization processing after each calculation layer, and finally outputting by using a Sigmoid activation function;
wherein, the Sigmoid activation function is:
Figure FDA0003767206580000012
s22: defining a discriminator model D consisting of three convolution layers and a full-link layer, wherein except the first convolution layer and the third convolution layer, the second convolution layer and the full-link layer are subjected to batch normalization processing, and finally LeakyRelu is used as an activation function to be output;
wherein, the LeakyRelu activation function is:
Figure FDA0003767206580000013
the step S3 includes the steps of:
s31: the noise generation unit N generates a random noise sequence z of a certain length i
S32: random noise sequence z i Inputting a generator model G and mapping the generator model G into a scene G (z) through a network of the generator model G;
s33: scene G (z) generated by generator model G and sampled training sample scene x i Inputting the values into a discriminator model D, wherein the discriminator model D respectively gives values D (x) and D (G (z)) representing the scene as true and false;
s34: respectively calculating loss functions of the discriminator model D and the generator model G according to the following formula;
L G =-E Z [D(G(z))],
L D =-E X [D(x)+E Z [D(G(z))]],
wherein L is G Representing the loss function of the generator, L D Loss function representing the discriminator, E Z Expected value, E, representing the corresponding random noise distribution X An expected value representing a corresponding true sample distribution;
s35: the result of S34 is reversely propagated to the generator model G and the discriminator model, and the network parameter theta is updated and adjusted according to the following formula g,t And a network parameter theta d,t
θ g,t =θ g,t-1 -▽ln(1-D(G(z))),
θ d,t =θ d,t-1 -▽ln(D(x))+ln(1-D(G(z))),
Wherein, theta g,t-1 And theta d,t-1 The parameters of the generator network and the parameters of the discriminator network D are generated before the t-th round of training begins;
s36: repeating steps S31-S35 untilTo network parameter θ g,t And a network parameter θ d,t Until convergence.
2. The method as claimed in claim 1, wherein the generator model G is a set of random arrays with a size of 96 x 1.
3. The method for generating the integrated energy system load scene based on the generation countermeasure network of claim 1, wherein the filter sizes of the first layer of convolutional layer, the second layer of convolutional layer and the third layer of convolutional layer are all 5 x 5.
4. The method as claimed in claim 1, wherein in step S4, the generator model G and the discriminator model D are fixed, and the generator model G is used for the random noise sequence z input from the input i And (4) generating and outputting phase cooling, heating and power load data in a simulation mode.
CN202010735071.1A 2020-07-28 2020-07-28 Comprehensive energy system load scene generation method based on generation countermeasure network Active CN111950868B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010735071.1A CN111950868B (en) 2020-07-28 2020-07-28 Comprehensive energy system load scene generation method based on generation countermeasure network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010735071.1A CN111950868B (en) 2020-07-28 2020-07-28 Comprehensive energy system load scene generation method based on generation countermeasure network

Publications (2)

Publication Number Publication Date
CN111950868A CN111950868A (en) 2020-11-17
CN111950868B true CN111950868B (en) 2022-11-15

Family

ID=73339673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010735071.1A Active CN111950868B (en) 2020-07-28 2020-07-28 Comprehensive energy system load scene generation method based on generation countermeasure network

Country Status (1)

Country Link
CN (1) CN111950868B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112598188A (en) * 2020-12-29 2021-04-02 沃太能源南通有限公司 Neural network generation method, power prediction method, device and storage medium
CN113256017A (en) * 2021-06-07 2021-08-13 国网山东省电力公司电力科学研究院 Short-term load prediction method and system
CN113269314B (en) * 2021-06-16 2022-10-28 福州大学 New energy power generation scene data migration method based on generation countermeasure network
CN113673159A (en) * 2021-08-19 2021-11-19 东北电力大学 Renewable energy space-time scene generation method based on federal deep-generation learning
CN114862123A (en) * 2022-04-12 2022-08-05 国网江苏省电力有限公司电力科学研究院 Comprehensive energy system scene generation method and device

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110490659A (en) * 2019-08-21 2019-11-22 广东电网有限责任公司 A kind of customer charge curve generation method based on GAN

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110490659A (en) * 2019-08-21 2019-11-22 广东电网有限责任公司 A kind of customer charge curve generation method based on GAN

Also Published As

Publication number Publication date
CN111950868A (en) 2020-11-17

Similar Documents

Publication Publication Date Title
CN111950868B (en) Comprehensive energy system load scene generation method based on generation countermeasure network
DE102017128693A1 (en) Feature and limit setting for threat detection in an industrial plant control system
CN104869126A (en) Network intrusion anomaly detection method
Moayyed et al. A Cyber-Secure generalized supermodel for wind power forecasting based on deep federated learning and image processing
CN113486078A (en) Distributed power distribution network operation monitoring method and system
CN113191429A (en) Power transformer bushing fault diagnosis method and device
CN111310722A (en) Power equipment image fault identification method based on improved neural network
CN112232488A (en) Deep learning and data driving-based new energy output scene generation method
Costilla-Enriquez et al. Attack power system state estimation by implicitly learning the underlying models
CN111192158A (en) Transformer substation daily load curve similarity matching method based on deep learning
Cao et al. Deep learning-based online small signal stability assessment of power systems with renewable generation
Li et al. A LSTM-based method for comprehension and evaluation of network security situation
CN117218425A (en) Power generation loss analysis method and system for photovoltaic power station
CN113988189B (en) Migration fault diagnosis method for cross-wind turbine generator
CN116523001A (en) Method, device and computer equipment for constructing weak line identification model of power grid
Xiang et al. An Improved Multiple Imputation Method Based on Chained Equations for Distributed Photovoltaic Systems
CN114298300A (en) Uncertainty prediction method, uncertainty prediction device, electronic device, and storage medium
CN113536508A (en) Method and system for classifying manufacturing network nodes
Li et al. Deep-learning-assisted power loss management for active distribution networks
CN113642784A (en) Wind power ultra-short term prediction method considering fan state
CN115392102A (en) Method and device for establishing energy consumption prediction model and method and system for predicting energy consumption
Chen et al. Electric power forecasting based on logistic regression and gradient boost regression tree
CN111008584A (en) Electric energy quality measurement deficiency repairing method of fuzzy self-organizing neural network
Xia et al. Smart substation network fault classification based on a hybrid optimization algorithm
Moradpour et al. A New Method Based on Symbolic Regression to Detect The Probability of False Data Injection Attacks on PV Generation

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