CN108090616A - A kind of electric system Active Splitting optimal section searching method - Google Patents

A kind of electric system Active Splitting optimal section searching method Download PDF

Info

Publication number
CN108090616A
CN108090616A CN201711399556.2A CN201711399556A CN108090616A CN 108090616 A CN108090616 A CN 108090616A CN 201711399556 A CN201711399556 A CN 201711399556A CN 108090616 A CN108090616 A CN 108090616A
Authority
CN
China
Prior art keywords
node
nodes
splitting
power
section
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.)
Pending
Application number
CN201711399556.2A
Other languages
Chinese (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.)
Southeast University
Liyang Research Institute of Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201711399556.2A priority Critical patent/CN108090616A/en
Publication of CN108090616A publication Critical patent/CN108090616A/en
Pending legal-status Critical Current

Links

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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply

Abstract

The present invention discloses a kind of electric system Active Splitting optimal section searching method, includes the following steps:Step 1, system weighted-graph is established according to operation states of electric power system information;Step 2, each load bus is calculated to the electrical link of each electric power generator group to ask for the Floyd algorithms of beeline in graph theory;Step 3, load bus is divided by general node and common node according to the electrical link degree difference between each load bus and different electric power generator groups;Step 4, the minimum target of system isolated island imbalance power asks for optimal off-the-line section after in common node region with off-the-line.Such method can rapidly and accurately determine Active Splitting section, while take into account isolated island imbalance power, off-the-line the operation two aspect problems of trend impact caused by system after off-the-line.

Description

Method for searching optimal section of active splitting of power system
Technical Field
The invention belongs to the field of safety and stability analysis of power systems, and particularly relates to a method for searching an optimal section of a power system by active splitting.
Background
The expansion of the scale of the grid frame of the electric power system and the diversification of the operation mode lead the complexity of the operation control of the electric power system to be increased day by day. In recent years, research and analysis on a series of major power failure accidents show that the running state of the power system can be accurately predicted in time, corresponding control measures are taken, and the risk resistance capability of the power system can be effectively improved. With the continuous promotion of intelligent power grid construction, the levels of an online data measurement technology and a high-speed communication technology are continuously improved, and important technical support is provided for constructing response-based wide-area safe and stable control of a power system and realizing real-time decision and real-time control.
The active splitting technology is one of important means for the safe and stable control of the power system, and can effectively avoid large-area power failure accidents caused by the spread of accidents. The research aiming at the active splitting technology mainly comprises two aspects of splitting criterion and splitting scheme. The existing research aiming at the search of the splitting section is mainly based on a method of graph theory, and methods such as an early ordered binary decision tree and the like obtain good results in the research of the search of the splitting section from the perspective of the graph theory, but the problem of higher calculation complexity needs to be solved when the system scale is enlarged. In order to break through the limitation of complexity, related research converts the search of the splitting section into the problem of solving the global minimum cut in the graph theory, and the calculation rate is effectively improved. Further, it is studied to simplify the search space of the splitting strategy by reducing the search space to quickly determine the splitting section, so as to further improve the solving speed, but the optimal solution may be lost when the search space is simplified.
Disclosure of Invention
The invention aims to provide a method for searching an optimal active splitting section of a power system, which can quickly and accurately determine the active splitting section and simultaneously take the two problems of unbalanced power of an island after splitting and tidal current impact on the system caused by splitting operation into consideration.
In order to achieve the above purpose, the solution of the invention is:
an active splitting optimal section searching method for a power system comprises the following steps:
step 1, establishing a system authorized undirected graph according to the running state information of the power system;
step 2, calculating the electrical connection from each load node to each generator group by using a Floyd algorithm for obtaining the shortest distance in the graph theory;
step 3, dividing the load nodes into general nodes and public nodes according to the difference of the electrical contact degree between each load node and different generator groups;
and 4, solving an optimal splitting section in the common node area by taking the minimum unbalanced power of the split system island as a target.
In the step 1, the ratio X of the reactance value of the line to the transmitted active power is taken into consideration of the energy transfer size and distance ij /P ij And measuring the tightness degree of the electrical connection between the nodes, wherein the smaller the value of the tightness degree is, namely the smaller the reactance is, the larger the transmitted active power is, the tighter the electrical connection between the adjacent nodes is, and further establishing a system weighted undirected graph adjacency matrix.
The details of the step 2 are as follows: for the power system weighted undirected graph, when the edge weight is set as the ratio of the line reactance to the transmission power, the average minimum edge weight summation of the load node to each generator in the generator group is obtained through the Floyd algorithm to measure the electrical connection strength of the load node to the generator group, as shown in the following formula:
in the formula, N m Representing the number of generators in the cluster m, W ij And representing the minimum edge weight accumulated sum between the node i and the generator node j in the cluster m.
In step 3, it is assumed that m is present for both clusters 1 And m 2 Let us orderWherein, the first and the second end of the pipe are connected with each other, are respectively a cluster m 1 、m 2 Strength of electrical connection of; taking a threshold value of epsilon, if gamma&Epsilon, the node i belongs to a public node with similar relation degree with the two groups; gamma ray to&And e, i belongs to a common node, i is classified into a cluster which is closely connected with the common node by comparing the positive and negative of gamma, wherein the threshold value e is specifically set according to statistical data and operation experience, and the common node is about 20 percent of the load node.
The value of epsilon is such that the common node is 20% of the load nodes.
In the step 4, the splitting section is obtained in the common node area by taking the minimum unbalanced power as a target: endowing the generator node with a weight P G The load node weight is-P L And sequentially searching common nodes from a certain cluster and a general node set area contained in the cluster, carrying out objective function check on the common nodes, if unbalanced power in the area is reduced, bringing the common nodes into the area until the common nodes meeting the conditions do not exist, and obtaining a final splitting section.
After adopting the scheme, compared with the prior art, the invention has the beneficial effects that:
the search space of the splitting section is reduced based on the electrical connection, so that the searching speed of the splitting section is improved. The method accurately measures the electrical connection between the nodes according to the ratio of the line reactance to the transmission power, does not cause large tidal current impact because the active power is introduced when the nodes are divided according to the electrical connection, takes the two problems of the unbalanced power of the island after splitting and the tidal current impact caused by the splitting operation to the system into consideration, and can better adapt to the requirement of the power system on the splitting scheme.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a diagram of a new England 39 node system architecture used to validate the present invention;
FIG. 3 is a cross section of an optimal solution found in accordance with the present invention.
Detailed Description
The technical solution and the advantages of the present invention will be described in detail with reference to the accompanying drawings.
As shown in fig. 1, the present invention provides a method for searching an optimal active splitting section of a power system, which comprises the following steps:
step 1, establishing a system authorized undirected graph according to the running state information of the power system;
considering the energy transfer size and distance, and taking the ratio X of the line reactance value to the transmitted active power size ij /P ij And measuring the tightness degree of the electrical connection between the nodes, wherein the smaller the value of the tightness degree is, namely the smaller the reactance is, the larger the transmitted active power is, the tighter the electrical connection between the adjacent nodes is, and further establishing a system weighted undirected graph adjacency matrix.
Constructing an adjacent matrix of a weighted undirected graph of the power system, and if any two nodes i and j in the system are adjacent, forming an edge of the undirected graph by a line connecting the two nodes, wherein an edge weight is set as a ratio of a line reactance to active power transmitted between i and j; if the two nodes are not adjacent, the corresponding value in the adjacent matrix is infinite; if two nodes actually represent the same node, the value in the adjacency matrix takes 0.
Step 2, calculating the electrical connection from each load node to each generator group by using a Floyd algorithm for obtaining the shortest distance in the graph theory;
the Floyd algorithm is a classical method for solving the shortest distance between two nodes in graph theory, and the core idea is to select an arc section with the smallest distance from unmarked points, further perform cyclic comparison and finally determine the shortest distance between the two nodes. For the power system weighted undirected graph, when the edge weight is set as the ratio of the line reactance to the transmission power, the minimum edge weight sum between two nodes obtained by the Floyd algorithm reflects the electrical connection strength between the two nodes to a certain extent. And then, the average minimum edge weight summation of the load node to each generator in the generator group can be obtained to measure the electrical connection strength from the load node to the generator group, as shown in formula (2):
in the formula, N m Representing the number of generators in the cluster m, W ij And representing the minimum edge weight accumulated sum between the node i and the generator node j in the cluster m.
Step 3, dividing the load nodes into general nodes and common nodes according to the difference of the electrical contact degree between each load node and different generator groups;
the load nodes in the system after the instability can be divided into two types, one type is a general node which has strong electrical connection with a certain cluster, and the other type is a common node which has similar electrical connection degree with two or more clusters. The majority of common nodes are divided into clusters closely connected with the common nodes, and the searching and splitting section substantially determines the splitting scheme to reasonably divide the common nodes into the corresponding clusters, thereby effectively reducing the searching space.
Suppose that for two clusters m 1 And m 2 Let us orderTaking a threshold value of epsilon, if gamma&If the node i belongs to a public node with similar relation degree with the two computer groups, the node i is a node with similar relation degree with the two computer groups; gamma ray&And e, i belongs to a common node, i is classified into a cluster which is closely connected with the common node by comparing the positive and negative of gamma, wherein the threshold value e is specifically set according to statistical data and operation experience, and the common node is about 20 percent of the load node.
And 4, solving an optimal splitting section in the common node area by taking the minimum unbalanced power of the split system island as a target.
The objective functions of the splitting section search mainly include two types, one is minimum unbalanced power between the generator and the load in the split region, and the other is minimum tidal current impact, and expressions are respectively shown in formulas (3) and (4):
min|P G -P L | (3)
in the formula, P G 、P L Respectively representing generator power and load power in the split region, P ij Representing the power transmission size between node i and node j, wherein i and j belong to different areas V 1 、V 2 . The minimum unbalanced power ensures less generator tripping and load shedding of the system after the splitting, and the minimum tide impact aims to reduce the impact on the system caused by the splitting operation. In the simplified search space, namely, the common node area, the invention obtains the splitting section by taking the minimum unbalanced power as a target: endowing the generator node with a weight P G The load node weight is-P L And sequentially searching common nodes from a certain cluster and a general node set area contained in the cluster, carrying out objective function check on the common nodes, if unbalanced power in the area is reduced, bringing the common nodes into the area until the common nodes meeting the conditions do not exist, and obtaining a final splitting section.
Taking the new england 39-node standard system as an example for testing, the system structure is shown in fig. 2, and the system is a standard system for analyzing the transient safety of the system, and comprises 10 generators and 46 power lines.
A fault condition: the 0s three-phase short circuit at the line 23-24, the 0.4s fault disappears, the generator under the fault is divided into two groups, wherein the generators 33-36 are the advanced group.
Firstly, an adjacency matrix is calculated according to line reactance and power flow, a system weighted undirected graph is established, the average minimum edge weight cumulative sum from a load node to two clusters is calculated respectively by adopting a Floyd algorithm, and the difference between the edge weight cumulative sum from each load node to the two clusters is obtained, as shown in Table 1.
TABLE 1 difference between the cumulative sums of the minimum edge weights
Taking the threshold e =0.035, the common nodes 2, 3, 10, 11, 12, 13, 14, 25 are available, while the clusters 33-36 comprise the common nodes 15, 16, 17, 18, 19, 20, 21, 22, 23, 24. The splitting of the common nodes into respective clusters with the aim of minimum unbalanced power, the resulting splitting pattern is shown in fig. 3.
In the splitting mode, the minimum unbalanced power of the two islands is 0.329p.u and 0.130p.u respectively, and the transmission power of a line is considered when the electrical connection between the nodes is measured, so that the splitting scheme can effectively avoid the line with larger power flow, and the condition that the system is subjected to larger power flow impact due to the splitting operation is avoided, and the power flow impact under the scheme is 0.805p.u. And because the invention reduces the search space of the column-searching section based on the electrical connection, the solving complexity is greatly reduced, the time has obvious advantages, the solving time in the process is about 11ms, and the invention can better adapt to the requirement of the actual system on the response rate.
In summary, the invention discloses a method for searching the optimal splitting section of an active splitting of a power system, so as to quickly and accurately search the optimal splitting section after the system is unstable. Firstly, measuring the electrical connection from a load node to a generator group according to an electrical connection index between nodes fusing line reactance and transmission power, and dividing the load node into a general node and a common node so as to reduce a search space of a splitting section; and then, quickly determining an optimal splitting section by taking the minimum unbalanced power of the split system island as a target. Compared with the existing splitting section searching method, the method provided by the invention has higher efficiency, can take two problems of unbalanced power of an island after splitting and tidal current impact on a system caused by splitting operation into consideration, and can better adapt to the requirement of a power system on a splitting scheme.
The above embodiments are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modifications made on the basis of the technical scheme according to the technical idea of the present invention fall within the protection scope of the present invention.

Claims (6)

1. A method for searching an optimal section for active splitting of a power system is characterized by comprising the following steps:
step 1, establishing a system authorized undirected graph according to the running state information of the power system;
step 2, calculating the electrical connection from each load node to each generator group by using a Floyd algorithm for obtaining the shortest distance in the graph theory;
step 3, dividing the load nodes into general nodes and public nodes according to the difference of the electrical contact degree between each load node and different generator groups;
and 4, solving an optimal splitting section in the common node area by taking the minimum unbalanced power of the split system island as a target.
2. The method for searching for the optimal section for active splitting of the power system as claimed in claim 1, wherein: in the step 1, the ratio X of the reactance value of the line to the transmitted active power is taken into consideration of the energy transfer size and the distance ij /P ij And measuring the tightness degree of the electrical connection between the nodes, wherein the smaller the value of the tightness degree is, namely the smaller the reactance is, the larger the transmitted active power is, the tighter the electrical connection between the adjacent nodes is, and further establishing a system weighted undirected graph adjacency matrix.
3. The method for searching for the optimal section for active splitting of the power system as claimed in claim 1, wherein: the details of the step 2 are as follows: for the power system weighted undirected graph, when the edge weight is set as the ratio of the line reactance to the transmission power, the average minimum edge weight summation of the load node to each generator in the generator group is obtained through the Floyd algorithm to measure the electrical connection strength of the load node to the generator group, as shown in the following formula:
in the formula, N m Representing the number of generators in the cluster m, W ij And representing the minimum edge weight accumulated sum between the node i and the generator node j in the cluster m.
4. The method for searching for the optimal section for active splitting of the power system as claimed in claim 1, wherein: in said step 3, it is assumed that m is for two clusters 1 And m 2 Let us orderWherein, the first and the second end of the pipe are connected with each other, are respectively a cluster m 1 、m 2 Strength of electrical connection of; taking a threshold value of epsilon, if gamma&Epsilon, the node i belongs to a public node with similar relation degree with the two groups; gamma ray&And e, i belongs to a common node, i is classified into a cluster which is closely connected with the common node by comparing the positive and negative of gamma, wherein the threshold value e is specifically set according to statistical data and operation experience, and the common node is about 20 percent of the load node.
5. The method for searching for the optimal section for active splitting of the power system as claimed in claim 4, wherein: the value of epsilon is such that the common node is 20% of the load nodes.
6. The method for searching for the optimal section for active splitting of the power system as claimed in claim 1, wherein: in the step 4, a splitting section is obtained in the common node area by taking the minimum unbalanced power as a target: endowing the generator node with a weight P G The weight of the load node is-P L From a certain machineAnd starting from the cluster and the general node set area contained by the cluster to sequentially search the common nodes, carrying out objective function check on the common nodes, if the unbalanced power in the area is reduced, bringing the common nodes into the area until the common nodes meeting the conditions do not exist, and obtaining a final splitting section.
CN201711399556.2A 2017-12-22 2017-12-22 A kind of electric system Active Splitting optimal section searching method Pending CN108090616A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711399556.2A CN108090616A (en) 2017-12-22 2017-12-22 A kind of electric system Active Splitting optimal section searching method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711399556.2A CN108090616A (en) 2017-12-22 2017-12-22 A kind of electric system Active Splitting optimal section searching method

Publications (1)

Publication Number Publication Date
CN108090616A true CN108090616A (en) 2018-05-29

Family

ID=62178379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711399556.2A Pending CN108090616A (en) 2017-12-22 2017-12-22 A kind of electric system Active Splitting optimal section searching method

Country Status (1)

Country Link
CN (1) CN108090616A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109871587A (en) * 2019-01-21 2019-06-11 南京铭越创信电气有限公司 A kind of method of the controlled off-the-line of electric system under extreme weather conditions
CN111585280A (en) * 2020-03-06 2020-08-25 山东大学 Temporary power grid splitting method and system for blocking multi-direct-current commutation failure linkage

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102244394A (en) * 2011-06-24 2011-11-16 山东大学 Two-stage initiative separation method based on normalized spectral clustering and constrained spectral clustering
CN102280885A (en) * 2011-08-29 2011-12-14 天津大学 Method for islanding electric power system into a plurality of subsystems
CN103259265A (en) * 2013-06-06 2013-08-21 山东大学 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix
CN106506357A (en) * 2016-12-06 2017-03-15 国网河南省电力公司信息通信公司 A kind of double route collocation method of power telecom network and device
CN106712059A (en) * 2017-02-15 2017-05-24 国网江苏省电力公司电力科学研究院 Initiative splitting optimal section searching method for power system based on convex optimization theory
CN107301469A (en) * 2017-05-02 2017-10-27 广东电网有限责任公司惠州供电局 A kind of power network bulk transmission grid recognition methods

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102244394A (en) * 2011-06-24 2011-11-16 山东大学 Two-stage initiative separation method based on normalized spectral clustering and constrained spectral clustering
CN102280885A (en) * 2011-08-29 2011-12-14 天津大学 Method for islanding electric power system into a plurality of subsystems
CN103259265A (en) * 2013-06-06 2013-08-21 山东大学 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix
CN106506357A (en) * 2016-12-06 2017-03-15 国网河南省电力公司信息通信公司 A kind of double route collocation method of power telecom network and device
CN106712059A (en) * 2017-02-15 2017-05-24 国网江苏省电力公司电力科学研究院 Initiative splitting optimal section searching method for power system based on convex optimization theory
CN107301469A (en) * 2017-05-02 2017-10-27 广东电网有限责任公司惠州供电局 A kind of power network bulk transmission grid recognition methods

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋洪磊 等: "考虑同调约束的电力系统主动解列断面的搜索方法", 《电力系统自动化》 *
王乙斐 等: "基于"搜索 优化"的主动解列断面选择方法", 《电力自动化设备》 *
谭健聪 等: "Floyd 算法在电力光缆路由分析选择中的应用研究", 《自动化与仪器仪表》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109871587A (en) * 2019-01-21 2019-06-11 南京铭越创信电气有限公司 A kind of method of the controlled off-the-line of electric system under extreme weather conditions
CN111585280A (en) * 2020-03-06 2020-08-25 山东大学 Temporary power grid splitting method and system for blocking multi-direct-current commutation failure linkage

Similar Documents

Publication Publication Date Title
CN107679716B (en) Interconnected power grid cascading failure risk assessment and alarm method considering communication vulnerability
AU2020103195A4 (en) A Method for Detecting Vulnerability of Large-scale Power Grid Based On Complex Network
CN109687438B (en) Power grid fragile line identification method considering high-speed rail impact load effect
CN105912835B (en) A kind of power circuit fragility identification model method for building up counted and wind power plant accesses
CN107292481B (en) Power grid key node evaluation method based on node importance
CN112288326B (en) Fault scene set reduction method suitable for toughness evaluation of power transmission system
CN110034581A (en) The electrical betweenness vulnerability assessment method in the section of electric system under wind-electricity integration
CN107368966B (en) Method for evaluating vulnerability of power network
CN109919801B (en) Coupling method and device based on node importance of power system
CN105354675A (en) Key transmission section identification-based cascading failure analysis method for alternating current/direct current power network
CN109033507A (en) A kind of Model in Reliability Evaluation of Power Systems method considering the failure of information system function for monitoring
CN107622360A (en) A kind of critical circuits recognition methods for considering subjective and objective factor
CN115545479A (en) Method and device for determining important nodes or important lines of power distribution network
CN111753420A (en) Cascade fault simulation method, system and storage medium for power information physical system
CN107591842B (en) Station domain island detection method
CN108090616A (en) A kind of electric system Active Splitting optimal section searching method
CN111062569A (en) Low-current fault discrimination method based on BP neural network
CN110768260B (en) Power grid cascading failure model building method based on electrical betweenness
CN117236030A (en) Power system toughness evaluation modeling method considering cascading overload fault occurrence under typhoon disaster
CN111191867B (en) Reliability evaluation method for complex network of power system
CN111429043A (en) Voltage sag characteristic quantity random evaluation method based on scene construction
CN103065193A (en) On-line intelligent identification method of provincial level power grid cascading failures
AU2021106383A4 (en) Method for Calculating Reliability of Distribution System Based on Equipment Reliability Cloud Model
CN115577292A (en) Vulnerable line identification method based on abstract dual network and cascading failure super network
CN114386222A (en) Power distribution network cooperative fault recovery method and system based on information physical system

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
TA01 Transfer of patent application right

Effective date of registration: 20201201

Address after: 213300 room 428, building a, 218 Hongkou Road, Kunlun Street, Liyang City, Changzhou City, Jiangsu Province (in Zhongguancun Science and Technology Industrial Park, Jiangsu Province)

Applicant after: Liyang Research Institute of Southeast University

Applicant after: SOUTHEAST University

Address before: Four pailou Nanjing Xuanwu District of Jiangsu Province, No. 2 210096

Applicant before: SOUTHEAST University

TA01 Transfer of patent application right
RJ01 Rejection of invention patent application after publication

Application publication date: 20180529

RJ01 Rejection of invention patent application after publication