CN110690716B - Method and system for positioning active splitting section of power grid based on voltage trajectory information - Google Patents
Method and system for positioning active splitting section of power grid based on voltage trajectory information Download PDFInfo
- Publication number
- CN110690716B CN110690716B CN201910760005.7A CN201910760005A CN110690716B CN 110690716 B CN110690716 B CN 110690716B CN 201910760005 A CN201910760005 A CN 201910760005A CN 110690716 B CN110690716 B CN 110690716B
- Authority
- CN
- China
- Prior art keywords
- node
- track
- node voltage
- voltage
- time sequence
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
Abstract
The invention discloses a method and a system for positioning an active splitting section of a power grid based on voltage trajectory information, wherein the method comprises the following steps: constructing a vector offset feature space of a node voltage phase track; acquiring the time sequence evolution characteristic of the node voltage when the power system is in step-out oscillation by utilizing the vector offset characteristic space of the node voltage trajectory; and according to the time sequence evolution characteristics of the node voltage, evaluating the similarity of the node voltage change based on a track clustering algorithm, and positioning one or more splitting sections by tracking the one or more splitting sections. According to the technical scheme, a time sequence evolution rule of node voltage change is extracted by constructing an offset characteristic space, the rationality of the rule is explained by a two-machine equivalent system, the similarity evaluation of node voltage tracks is realized by a track distance-based adaptive clustering algorithm, and then the splitting section is accurately positioned on line through cluster expansion and power self-balancing constraint.
Description
Technical Field
The invention relates to the technical field of transient stability control of a large power grid, in particular to a method and a system for positioning an active splitting section of the power grid based on voltage trajectory information.
Background
With the rapid development of ultrahigh voltage alternating current and direct current power grids in China, the grid pattern and the power supply structure are greatly changed, the operating characteristics of the power grids are deeply changed, the traditional power grid online safety defense concept and stability control technology taking modeling simulation and expected faults as the core are difficult to adapt to the power grid development requirements, and three defense systems for safety and stability of the current power grids face severe challenges. The out-of-step separation is used as the last defense line for ensuring the safe and stable operation of the power grid, and is one of important prevention and control measures for restraining further propagation of the power failure accident of the interconnected power grid. However, with the continuous expansion of the system scale and the formation of the cross-regional interconnected power grid, the passive splitting section positioning method based on the oscillation center positioning is difficult to adapt to complex fault forms, and is very easy to cause serious consequences.
In recent years, with the popularization and application of Wide Area Measurement Systems (WAMS), high-precision real-time monitoring of the operation state of a power grid becomes possible. Active splitting methods based on online decision-making become current research hotspots, which are mainly classified into the following three categories:
(1) the method is mainly characterized in that the method comprises a network simplification method and a rapid network division method. The network simplification method sacrifices the integrity of the system in order to reduce the search scale and accelerate the solving speed, which may cause the loss of a feasible solution; the rapid network partition method can effectively avoid traversal search of feasible solutions, but the final solution often does not meet the connectivity constraint; (2) the method comprises the steps of searching a splitting section based on a slow coherent theory, extracting a dynamic mode of a power system and analyzing weak connection among clusters by constructing a singular perturbation model with double time scales, and identifying coherent clusters and weak links. However, the method needs to calculate characteristic values and characteristic vectors, the calculation amount is too large, and the solving speed of the optimal splitting section is forced to be reduced; (3) and (4) searching a splitting section based on an intelligent algorithm. The active splitting of the power grid has higher requirements on timeliness, and the optimal splitting section needs to be calculated within a limited time. The optimal splitting section solving problem based on the optimization theory is substantially an NP-hard problem, and researches prove that the optimal splitting section solving problem does not have an accurate solution in linear time complexity, so that the searching speed and the accuracy can be improved by solving an approximate solution or converting an objective function, and the stability of an isolated island after splitting is ensured.
In the actual large power grid splitting process, the splitting section candidate space is a set of all branch circuit breaking combinations, so that the splitting scheme is increased in a geometric exponential mode along with scale improvement. For the slow coherent theory and the intelligent optimization algorithm, apart from inherent defects of the slow coherent theory and the intelligent optimization algorithm, the slow coherent theory and the intelligent optimization algorithm are limited by system scale in different degrees and face huge calculation pressure, so that the second-level solving speed of practical power grid splitting control is difficult to meet, and the slow coherent theory and the intelligent optimization algorithm can be applied to on-line operation only by further research and improvement.
Therefore, a technique is needed to realize the positioning of the active splitting section of the power grid based on the voltage trajectory information.
Disclosure of Invention
The technical scheme of the invention provides a method and a system for positioning an active splitting section of a power grid based on voltage trajectory information, so as to solve the problem of how to position the active splitting section of the power grid based on the voltage trajectory information.
In order to solve the above problem, the present invention provides a method for positioning an active splitting section of a power grid based on voltage trajectory information, the method including:
constructing a vector offset feature space of a node voltage phase track;
acquiring the time sequence evolution characteristic of the node voltage when the power system is in step-out oscillation by utilizing the vector offset characteristic space of the node voltage trajectory;
and according to the time sequence evolution characteristics of the node voltage, evaluating the similarity of the node voltage change based on a track clustering algorithm, and positioning one or more splitting sections by tracking the one or more splitting sections.
Preferably, the constructing a vector offset feature space of the nodal voltage phase trajectories comprises:
with d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset feature space psi of a node voltage phase track d Let a node start point V mi The coordinate is (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) And the vector offset characteristic track of the vector between the starting point and the end point forming node i in a certain time interval is as follows:
constructed ofAs a node at d l 、d θ 、d p Judging basis of movement direction change and speed change in a certain time interval in space;component values in the x, y, and z axes characterize the rate, direction, and position magnitudes of node offsets, respectively.
Preferably, the obtaining, by using a vector offset feature space of the node voltage trajectory, a time sequence evolution feature of the node voltage during step-out oscillation of the power system includes:
collecting the time sequence information of the node voltage in the power grid through a wide area measurement system to obtain the time sequence of the node voltage in a certain time intervalInformation forming a group UR of voltage phasor trajectories for all nodes in the voltage complex space within the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; and n is the number of all nodes of the power grid.
Preferably, the estimating, according to the time-series evolution characteristic of the node voltage, the similarity of the node voltage changes based on a trajectory clustering algorithm, and locating one or more splitting sections by tracking the one or more splitting sections, further includes:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree between time sequence tracks and measuring the difference degree between different tracks; the distance between two discrete points at a certain time in the space is known as follows:
the distance D can be used to generalize to the inter-track distance D within a certain time interval Traj Is calculated as follows:
D Traj =D 1 +D 2 +···+D n
in the formula D n The distance between two points at the nth moment is expressed, and n is the window length of the time sliding window;
the distance between any two tracks can be obtained by sliding the window length by one unit.
Preferably, the method further comprises the following steps:
performing cluster expansion on the neighborhood track group, and offsetting tracks of each vector in the neighborhoodPerforming secondary traversal and label assignment by applying to the neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
According to another aspect of the present invention, there is provided a system for positioning an active splitting section of a power grid based on voltage trajectory information, the system comprising:
the building unit is used for building a vector offset characteristic space of the node voltage phase track;
the acquisition unit is used for acquiring the time sequence evolution characteristics of the node voltage during the step-out oscillation of the power system by utilizing the vector offset characteristic space of the node voltage trajectory;
and the positioning unit is used for evaluating the similarity of the node voltage change based on a track clustering algorithm according to the time sequence evolution characteristics of the node voltage and positioning one or more splitting sections by tracking the one or more splitting sections.
Preferably, the construction unit is configured to construct a vector offset eigenspace of the node voltage phase trajectories, and further configured to:
with d l 、dθ、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset feature space psi of a node voltage phase track d Let a node start point V mi The coordinates are (d) lmi ,dθ mi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,dθ (m+1)i ,d p(m+1)i ) The vector between the starting point and the end point forms a node i in a certain time intervalThe vector offset feature trajectory of (a) is:
constructed ofAs a node at d l 、dθ、d p Judging basis of motion direction change and speed change in a certain time interval in space;component values in the x, y, and z axes characterize the rate, direction, and position magnitudes of node offsets, respectively.
Preferably, the obtaining unit is configured to obtain a time sequence evolution characteristic of the node voltage when the power system is in step-out oscillation by using a vector offset feature space of the node voltage trajectory, and is further configured to:
collecting the time sequence information of the node voltage in the power grid through a wide area measurement system, acquiring the time sequence information of the node voltage in a certain time interval, and forming a voltage phasor track group UR of all nodes in a voltage complex space in the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; and n is the number of all nodes of the power grid.
Preferably, the positioning unit is configured to evaluate similarity of changes of the node voltages based on a trajectory clustering algorithm according to a time-series evolution characteristic of the node voltages, and position one or more splitting sections by tracking the one or more splitting sections, and is further configured to:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree between time sequence tracks and measuring the difference degree between different tracks; the distance between two discrete points at a certain time in the known space is:
the distance D can be used to generalize to the inter-track distance D within a certain time interval Traj Is calculated as follows:
D Traj =D 1 +D 2 +···+D n
in the formula D n Expressed as the distance between two points at the nth moment, n is the window length of the time sliding window;
the distance between any two tracks can be obtained by sliding the window length by one unit.
Preferably, the positioning unit is further configured to:
performing cluster expansion on the neighborhood track group, and offsetting tracks of each vector in the neighborhoodPerforming secondary traversal and label endowing by aiming at neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
The technical scheme of the invention provides a method and a system for positioning an active splitting section of a power grid based on voltage track information, wherein the method comprises the following steps: constructing a vector offset feature space of a node voltage phase track; acquiring a time sequence evolution characteristic of the node voltage when the power system is in step-out oscillation by utilizing a vector offset characteristic space of a node voltage track; and evaluating the change of the node voltage based on a track clustering algorithm according to the time sequence evolution characteristics of the node voltage, and positioning one or more splitting sections by tracking the one or more splitting sections. The technical scheme of the invention is switched in from a brand-new visual angle of a voltage track, provides a quick positioning method of an active splitting section and aims to save the survival capability of a power grid main body. According to the technical scheme, a time sequence evolution rule of node voltage change is extracted by constructing an offset characteristic space, the rationality of the rule is explained by a two-machine equivalent system, the evaluation of the node voltage track similarity is realized by a track distance-based adaptive clustering algorithm, and then the splitting section is accurately positioned on line through cluster expansion and power self-balancing constraint. The technical scheme of the invention is independent of a mathematical model, is not limited by a system operation mode and a fault form, does not need complex calculation, has strong section connectivity, and has important significance for realizing transient stability active defense of an information-driven large power grid.
Drawings
A more complete understanding of exemplary embodiments of the present invention may be had by reference to the following drawings in which:
FIG. 1 is a flow chart of a method for locating an active splitting section of a power grid based on voltage trajectory information in accordance with a preferred embodiment of the present invention;
FIG. 2 is a schematic diagram of feature space construction according to a preferred embodiment of the present invention;
FIG. 3 is a system diagram of an IEEE-9 node in accordance with a preferred embodiment of the present invention;
FIG. 4 is a graph of a characteristic trace for a system steady after a fault in accordance with a preferred embodiment of the present invention;
FIG. 5 is a characteristic trace diagram for a post-fault system instability condition in accordance with a preferred embodiment of the present invention;
FIG. 6 is a schematic diagram of a two-machine model of a power system in accordance with a preferred embodiment of the present invention;
FIG. 7 is a graph of node voltage traces in accordance with a preferred embodiment of the present invention;
FIG. 8 is a graph of node voltage phase angles in accordance with a preferred embodiment of the present invention;
FIG. 9 is a flow chart of a clustering algorithm according to a preferred embodiment of the present invention;
FIG. 10 is a diagram illustrating a cluster expansion procedure in accordance with a preferred embodiment of the present invention;
FIG. 11 is a graph showing simulation results according to a preferred embodiment of the present invention;
FIG. 12 is a graph showing simulation results according to a preferred embodiment of the present invention;
FIG. 13 is a graph showing simulation results according to a preferred embodiment of the present invention;
FIG. 14 is a cross-sectional transition diagram illustrating active splitting according to a preferred embodiment of the present invention;
FIG. 15 is a graph of simulation results for an example of an actual grid in accordance with a preferred embodiment of the present invention; and
fig. 16 is a system configuration diagram for locating an active grid splitting section based on voltage trajectory information according to a preferred embodiment of the present invention.
Detailed Description
Example embodiments of the present invention will now be described with reference to the accompanying drawings, however, the invention may be embodied in many different forms and should not be construed as limited to the embodiments set forth herein, which are provided for a complete and complete disclosure of the invention and to fully convey the scope of the invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. In addition, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their context in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
Fig. 1 is a flowchart of a method for locating an active splitting section of a power grid based on voltage trajectory information according to a preferred embodiment of the present invention. In recent years, the grid pattern and the power supply structure are greatly changed, the operating characteristics of the power grid are deeply changed, the traditional power grid online safety defense concept and the stability control technology are difficult to adapt to the power grid development requirements, and the current three-defense system for the safety and the stability of the power grid faces a severe challenge. The implementation mode of the application provides a quick positioning method of an active splitting section by cutting in with a brand-new visual angle of a voltage track, and aims to save the survival capacity of a power grid main body. The method comprises the steps of extracting a time sequence evolution rule of node voltage change by constructing an offset characteristic space, evaluating the similarity of node voltage tracks by a track distance-based adaptive clustering algorithm, and accurately positioning a splitting section on line by cluster expansion and power self-balancing constraint. The effectiveness of the method is verified by the IEEE-39 calculation and the practical regional interconnected network calculation, complex calculation is not needed, time consumption is short, section connectivity is strong, and the method has a certain engineering practice value. The invention provides a novel method for quickly positioning an active splitting section of a power grid, which achieves the purpose of the invention by constructing an offset characteristic space, extracting a time sequence evolution rule of node voltage change, clustering and positioning the splitting section based on a characteristic track and adaptively adjusting an algorithm. As shown in fig. 1, a method for positioning an active splitting section of a power grid based on voltage trajectory information includes:
preferably, in step 101: and constructing a vector offset feature space of the node voltage phase tracks. Preferably, constructing a vector offset eigenspace of the node voltage phase trajectories comprises:
with d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset feature space psi of a node voltage phase track d Let a node start point V mi The coordinate is (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) The vectors between the starting and ending points constituting the node i in a certain time intervalThe vector offset feature trajectory is:
constructed ofAs a node at d l 、d θ 、d p Judging basis of motion direction change and speed change in a certain time interval in space;component values in the x, y, and z axes characterize the rate, direction, and position magnitudes of node offsets, respectively.
The voltage phase locus is a directed locus of bus voltage phasor in a complex number space in the space-time motion behavior of the power grid. The track has the characteristic of complete vector and can be regarded as continuous directional splicing of infinite direction vectors on a time scale. The method is based on the existing vector distance definition and adjacent time intervals T in a voltage complex space m 、T m+1 Fig. 2(a) shows a 2-segment trajectory from which the node i is obtained. The vector offset of the 2 nd segment with respect to the 1 st segment track can be represented by d l(m+1)i 、dθ (m+1)i And d p(m+1)i And (5) characterizing. Similarly, the vector offset of the 1 st segment of the track relative to the track in the previous interval can be represented by d lmi 、d θmi And d pmi And (5) characterizing.
In the present application d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset feature space psi of a node voltage phase track d As shown in fig. 2 (b). And set a starting point V mi The coordinates are (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) Section of vector construction between two pointsThe Vector offset feature Trajectory (VVMCT) of point i in a time interval is:
constructed by the formula (1)Can be regarded as a node at d l -d θ -d p And judging the change of the moving direction and the change of the speed in a certain time interval in the space.Component values in the x, y, and z axes characterize the rate, direction, and position magnitudes of node offsets, respectively.
Preferably, at step 102: the method for acquiring the time sequence evolution characteristics of the node voltage during the out-of-step oscillation of the power system by utilizing the vector offset characteristic space of the node voltage trajectory comprises the following steps:
collecting time sequence information of node voltage in the power grid through a wide area measurement system, acquiring the time sequence information of the node voltage in a certain time interval, and forming a voltage phasor track group UR of all nodes in a voltage complex space in the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; and n is the number of all nodes of the power grid.
The time sequence information of the node voltage in the power grid can be acquired through the wide area measurement system. By obtaining the node voltage time sequence information in a certain time interval, the voltage complex space full-scale in the time interval can be formedVoltage phasor trajectory group UR of partial nodes Tm :
In the formula (I), the compound is shown in the specification,represents the voltage vector of the node i in the mth time interval Tm; and n is the number of all nodes of the power grid.
In order to illustrate the necessity of constructing the offset feature space and verify the validity of the method, the present application takes an I EEE 9 node system shown in fig. 3 as an example, and analyzes the node voltage motion conditions of the system after the fault in two forms, namely, stable state and unstable state.
1) Three-phase short circuit faults occur at 50% of the positions of the connecting lines between the buses 2 and 4, the fault lines are cut off after 0.1s, and the system is kept stable. Based on the above method, the motion process of the node voltage in the transient process 1s is projected to the feature space, and the trajectories in the feature space and the three feature tangential planes are shown in fig. 4. The motion tracks of all nodes in the system form 1 characteristic track cluster, no obvious layering classification phenomenon occurs, and the track of the cluster has motion convergence, which shows that all nodes generally keep the consistency of motion directions from the moment of fault removal, so that the system starts to transit to the next stable operation point; if the observation is carried out by cutting in different characteristic planes, the difference between the tracks in the clusters explains the relative swing in the movement process of the nodes, so that the difference of the direction, the speed and the offset on the position exists between the nodes.
2) Three-phase short circuit faults occur at 50% of the positions on the connecting lines between the buses 2 and 4, after 0.16s, the fault lines are cut off, and the generator G1 and the system are out of step. The trajectory in the feature space and the three feature tangential planes is shown in fig. 5. The motion tracks of all nodes in the system form 2 track clusters, and an obvious classification phenomenon occurs; if the two clusters of tracks are cut into and observed from different characteristic planes, a certain distance is always kept between the three attributes, and the difference of the two clusters of tracks in the motion process indicates the difference of the node motion in the system, so that the system instability is a specific expression.
The aggregation effect shown in fig. 5 is regressed to the network topology of the system, and it can be easily found that 4 tracks in the sibling track cluster (ii) are respectively G1, Bus1, Bus4 and Bus5, a cut set section can be formed exactly in space, and G1 in the cut set is a destabilizing unit. The reason for this is that the instability unit G1 pulls its surrounding nodes to move in a homodromous manner during the desynchronization process, and deviates from the motion direction of the synchronous cluster and its surrounding nodes, so that the offset characteristic trajectories of the nodes are clustered into two clusters, and the motion offset difference of the two clusters of nodes is embodied as cut set formation in the topological space.
The trajectory in the offset feature space substantially represents the offset attribute of continuous motion of the node in a certain period, and is a concrete embodiment of continuous accumulation of motion offset. Because the similarity of the offset attributes essentially represents the motion similarity, it is necessary to extract the time sequence evolution rule of the node voltage variation by using the offset feature space as a medium and observe the aggregation condition of the tracks.
When a serious fault occurs in an actual complex power grid, more than one area threatening the safety of the main grid is possible, and the stability of the main grid can be effectively ensured only by splitting a plurality of power transmission sections urgently. Therefore, the accurate positioning of the multiple splitting sections plays an extremely important role in improving the survival capability of the complex power grid in the face of serious faults. The invention realizes the similarity evaluation of node motion by analyzing the motion conditions of node voltages of a simple system under two typical stable forms and based on a characteristic space, and the aggregation condition of the track cluster can provide important reference for disconnection control.
Taking an equivalent two-machine system as an example, the rationality of the evolution rule of the extracted node voltage change time sequence is analyzed. For convenience, it is assumed that the two potentials of the system shown in fig. 6 are equal in magnitude and the oscillation center falls at the center c, i.e., the center of symmetry. a. b is a node far away from the oscillation center and close to the sending end, d and e are nodes close to the receiving end, and X is system connection impedance.
Knowing the node voltage amplitudeThe voltage phase angle θ is arctan (Im/Re), where Re and Im are real and imaginary parts of the voltage, respectively. Bringing it into dV/dt and d θ/dt, respectively, there are:
the power angle difference of the two motors changes periodically within the range of 0-360 degrees, the voltage phase tracks of the nodes on the contact section are circles with different radiuses, as shown in fig. 7, and the phase angle change of each node in the oscillation process is as shown in fig. 8. Considering the symmetry of the oscillation process, the present application addresses the case where δ is at [0 °,180 ° ]]Detailed analysis of the oscillation process of (1), S in the figure 1 The phase angle in the corresponding oscillation process is maximum.
As can be seen from the voltage phase locus and vector relationship of fig. 7, the voltage variation amplitude at the oscillation center c is the largest. If and only if Re Im is 0, the voltage amplitude at c is zero, so the phase locus circle of the oscillation center crosses the origin, and the imaginary axis and the phase locus circle are tangent to the origin, where the phase angle at node c is 90 °.
The phase track circle of the node M is the outer boundary of the circle cluster, the voltage amplitude is kept unchanged in the oscillation process, and the relation between the electrical compaction and imaginary parts of the node and the angular speed direction can be obtained by substituting the condition into a formula (3) and then connecting the formula (4):
for nodes d and e far away from the oscillation center and close to the receiving end, the phase locus circle is in the locus circle of the oscillation center c, and the closer to the c, the larger the circle radius is, and the larger the voltage change is. The phase angles of the nodes have maximum values in the oscillation process, and the formula (4) is zero to derive the formula (6) to show the extreme value point S 1 The tangent line of the point (A) must pass through the origin, and the phase angle begins to decrease after passing through the point (A), and finallyAnd finally to zero.
For nodes a and b far away from the oscillation center and close to the sending end, the phase locus circle is between locus circles of M and c, and the closer to c, the smaller the radius of the circle is, and the larger the voltage change is. The phase angle of the class node monotonically increases to 180 ° during oscillation.
In combination with the above analysis: the phase locus circle of the oscillation center c is a boundary line of two types of node motions, and finally a and b which move along with the node M intersect with the virtual axis positive half shaft in the phase angle increasing process, the voltage locus of the type of node crosses the first quadrant and the second quadrant, the locus mode is relatively long, and the motion characteristics of the locus can be captured by the length distance and the position distance in the deviation characteristics. D and e which finally move along with the node N, the phase angle of the d and e has a maximum value, and an extreme point S 1 The tangent line of the position crosses the origin, the voltage phase track is shown to move only in the first quadrant, the included angle of two vectors in adjacent time intervals reaches a zero value in advance, and the motion characteristic of the track can be captured by the position distance and the angle distance.
The motion difference of the two types of nodes is expressed as the difference of the geometric characteristics of the phase tracks, and the motion similarity of the nodes of the same type is expressed as the similarity of the geometric characteristics of the phase tracks. By mining the geometric characteristics of the phase trajectory, extracting the time sequence evolution rule of the node state, analyzing the time-space linkage relation of the node state, further identifying the splitting section, and having feasibility.
Preferably, in step 103: according to the time sequence evolution characteristics of the node voltage, the similarity of the change of the node voltage is evaluated based on a track clustering algorithm, and one or more splitting sections are positioned by tracking the one or more splitting sections, and the method also comprises the following steps:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree between time sequence tracks and measuring the difference degree between different tracks; known hollowThe distance between two discrete points at a certain time in the interval is as follows:
the distance D can be used to generalize to the inter-track distance D within a certain time interval Traj Is calculated as follows:
D Traj =D 1 +D 2 +···+D n
in the formula D n The distance between two points at the nth moment is expressed, and n is the window length of the time sliding window;
the distance between any two tracks can be obtained every time the window length slides by one unit.
Preferably, cluster expansion is performed on the neighborhood track group, and the track is shifted for each vector in the neighborhoodPerforming secondary traversal and label endowing by aiming at neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
According to the method, the aggregation condition of the node voltage in the feature space under different stable forms is considered, the time sequence track in the feature space is extracted through a time sliding window, the node motion similarity is evaluated on line based on a track clustering algorithm, the migration condition of a multi-solution section is tracked in real time, and decision support is provided for the solution control of the system.
Selecting a time sequence track in a certain time interval compared with characteristic scatter clustering of all nodes under a certain time sectionClustering is better at capturing the time-evolution characteristics of node motion. Therefore, a concept of "track distance" is proposed, which is defined as the sum of the distances between points of different tracks at the same time scale and is recorded as D Traj The method is used for representing the similarity degree between the time sequence tracks and measuring the difference degree between different tracks. The distance between two discrete points at a certain time in the space is known as follows:
the distance D can be used to generalize to the inter-track distance D within a certain time interval Traj Is calculated as follows:
D Traj =D 1 +D 2 +···+D n (8)
in the formula D n Is expressed as the distance between two points at the nth time, and n is the window length of the time sliding window.
From the above analysis, it can be known that the distance between any two tracks can be obtained every time the time window slides by one unit, and the VVMCT clustering algorithm is proposed based on the distance, and the steps are shown in fig. 9. Of unstable power generation nodesIs likely to be in contact with the partThe large track distance exists because the moving direction of the off-line unit is completely deviated from the main network, the average track distance between the off-line unit and the nodes around the off-line unit linked by the off-line unit is large because the moving speed is too high, the off-line unit is higher than a clustering threshold value and cannot be clustered, and the unstable power generation nodes are isolated, which is a normal result of clustering. Therefore, the isolated track of the clustering result is optimized, and the optimization constraint purpose of the power generation nodes is set, namely, the isolated track close to the unstable power generation nodes is clustered again in space.
If the method of the application only evaluates the strong association relationship between the motion trails according to the result of one traversal, the result may be too many to be one-sided. Therefore, it is necessary to perform cluster expansion on the neighborhood trajectory groupFor each strip in the neighborhoodThe specific steps of performing the second traversal and the label assignment are shown in fig. 10. By making a neighborhoodInner trackPerforming secondary traversal, and searching and matching in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag. The cluster expansion ensures the accuracy and global optimality of cluster division, and is a key step of the clustering algorithm.
The clustering result of the algorithm is basically based on the classification evaluation of the node voltage trajectory similarity, and one or more sections are positioned through boundary nodes in a cluster where a destabilizing machine is located. In a large power grid with complex topology, some intermediate and load nodes with certain electrical distance from each power generation node are far away from the unstable power generation node and close to the oscillation center, track characteristics are not obvious, and clusters can be formed by self during clustering. When the splitting control is carried out, the cluster nodes can completely interrupt power supply due to lack of support of a power supply source, and if important loads which can not be powered off are contained in the cluster, the social production life is seriously influenced by power failure, so that huge economic loss is brought. Based on the above analysis, the present invention proposes three constraint principles for this type of node, aiming to ensure uninterrupted power supply of important loads, the principles are as follows: principle (I): the load nodes are preferentially divided into clusters where the cluster centers closest to the load nodes are located; principle II: after the load nodes are subdivided into adjacent regions, the unbalanced power in the regions is required to be ensured not to exceed 10%; principle (c): the intermediate nodes are allowed to stand alone if necessary.
(4) Adaptive adjustment of the algorithm:
the track distance threshold epsilon and the track cluster threshold lambda are two keys in the clustering algorithmThe reasonability of the values of the parameters directly influences the accuracy of the result. The invention provides a parameter self-adaptive adjusting method, firstly for D Traj Carrying out weighted average on the sample space, and if the number of all tracks in the VVMCT group is n, then D Traj The number of trajectory distances contained in the sample space Ω is:
in calculating D Traj The samples with overlarge values need to be excluded during the mean value, so the samples in omega are screened, and all D are detected Traj After the samples are sorted from small to large, the first 25 percent of D is extracted Traj Samples, make up the optimization space:
in the formula, N 1 N/4, i.e. the rounded number of optimized spatial samples, di being D in the optimized space Traj And (4) sampling. The optimization space largely preserves the threshold information because the sample spacing within each class is relatively small. Arranging samples in the optimized space from small to large and dividing the samples into four groups of subspaces, arranging threshold information contained in the four groups of subspaces from large to small, and calculating a weighted average value after giving different weights to the subspaces, wherein the weighted average value comprises the following steps:
in the formula (I), the compound is shown in the specification,for the weight value of each subspace, the weight a is taken as the neighborhood of the value 1 in principle, and as the threshold information density of each subspace is sequentially decreased, the weight selection should also be sequentially decreased, which can be determined according to the operation experience of the regulation and control personnel.
If the control person wishes to split the system as small as possible, i.e. the threshold values ε are relativeSmaller, at this time, increase a 1 Values are taken to enhance the effect of the 1 st subspace on the optimization space, but it should be noted that too large a value will result in isolated trajectories. To sum up, this application gets: a is a 1 ,a 2 ,a 3 ,a 4 1.6, 1.2, 0.8, 0.4, wherein N 2 =N 1 And/4, the number of samples in the subspace after the whole treatment is taken.
The adaptive adjustment method of the epsilon value dynamically adjusts the target threshold according to the sample space, and realizes the adaptive adjustment of the parameters of different VVMCT groups. And if the power grid topologies are the same, the node numbers are the same, and the lambda value is unchanged. In order to adapt to different numbers of nodes in different power grid topologies, the method selects the following steps:
λ=n/10 (12)
in the formula, lambda also needs rounding treatment. The track clustering number matched with the power grid node number can be rapidly calculated through the formula (19), so that the reasonability of the cluster threshold value is improved.
The application provides a new method for rapidly positioning an information-driven active splitting section of a power grid, which has the following characteristics:
d constructed in this application l -d θ -d p Characteristic space psi d The deviation attribute of the node motion track can be effectively represented, the time-space dynamics characteristics of a power grid and the time sequence evolution rule of node voltage change are reflected, and the rationality of the rule is explained through a two-machine equivalent system;
the self-adaptive clustering algorithm based on the track distance is provided, the data-driven node motion similarity classification evaluation is substantial, the method is independent of a model, is not limited by a system operation mode and a fault form, and is high in applicability;
according to the method and the device, the migration condition of the splitting section is tracked in real time through the sliding time window, the calculation is simple, the evaluation timeliness is good, the survival capability of a main body of a power grid is effectively preserved, and meanwhile, the online requirement of active splitting of the large power grid is met.
Fig. 16 is a system configuration diagram for locating an active grid splitting section based on voltage trajectory information according to a preferred embodiment of the present invention. The application provides a system for positioning power grid initiative splitting section based on voltage track information, the system includes:
and the constructing unit 601 is used for constructing a vector offset feature space of the node voltage phase trajectory. Preferably, the construction unit 601 is configured to construct a vector offset eigenspace of the node voltage phase trajectories, and further configured to:
with d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset feature space psi of a node voltage phase track d Let a node start point V mi The coordinate is (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) The vector offset characteristic track of the node i in a certain time interval formed by the vectors between the starting point and the end point is as follows:
constructed ofAs a node at d l 、d θ 、d p Judging basis of movement direction change and speed change in a certain time interval in space;component values in the x, y, and z axes characterize the rate, direction, and position magnitudes of node offsets, respectively.
The obtaining unit 602 is configured to obtain a time sequence evolution characteristic of the node voltage during step-out oscillation of the power system by using a vector offset characteristic space of a node voltage trajectory. Preferably, the obtaining unit 602 is configured to obtain a time-sequence evolution characteristic of the node voltage when the power system is out-of-step oscillated, by using a vector offset feature space of the node voltage trajectory, and further configured to:
collecting time sequence information of node voltage in the power grid through a wide area measurement system to obtain the node voltage in a certain time intervalForming a voltage phasor trajectory group UR of all nodes in the voltage complex space within the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; and n is the number of all nodes of the power grid.
And the positioning unit 603 is configured to evaluate similarity of node voltage changes based on a trajectory clustering algorithm according to a time sequence evolution characteristic of the node voltage, and position one or more splitting sections by tracking the one or more splitting sections.
Preferably, the positioning unit 603 is configured to evaluate a change of the node voltage based on a trajectory clustering algorithm according to a time-series evolution characteristic of the node voltage, and position the one or more splitting sections by tracking the one or more splitting sections, and further configured to:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree between time sequence tracks and measuring the difference degree between different tracks; the distance between two discrete points at a certain time in the space is known as follows:
the distance D can be used to generalize to the inter-track distance D within a certain time interval Traj Is calculated as follows:
D Traj =D 1 +D 2 +···+D n
in the formula D n The distance between two points at the nth moment is expressed, and n is the window length of the time sliding window;
the distance between any two tracks can be obtained every time the window length slides by one unit.
Preferably, the positioning unit 603 is further configured to:
performing cluster expansion on the neighborhood track group, and offsetting tracks of each vector in the neighborhoodPerforming secondary traversal and label assignment by applying to the neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
The system 600 for positioning an active splitting section of a power grid based on voltage trajectory information in the preferred embodiment of the present invention corresponds to the method 100 for positioning an active splitting section of a power grid based on voltage trajectory information in the preferred embodiment of the present invention, and is not described herein again.
The invention has been described with reference to a few embodiments. However, other embodiments of the invention than the ones disclosed above are equally possible within the scope of these appended patent claims, as these are known to those skilled in the art.
Generally, all terms used in the claims are to be interpreted according to their ordinary meaning in the technical field, unless explicitly defined otherwise herein. All references to "a/an/the [ device, component, etc ]" are to be interpreted openly as referring to at least one instance of said device, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless explicitly stated.
Claims (6)
1. A method for locating an active splitting section of a power grid based on voltage phase trajectory information, the method comprising:
constructing a vector offset feature space of the node voltage phase trajectory, comprising:
with d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset eigenspace psi of a node voltage phase track d Let a node start point V mi The coordinates are (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) And the vector offset characteristic track of the vector between the starting point and the end point forming node i in a certain time interval is as follows:
wherein the content of the first and second substances,as a node at d l 、d θ 、d p Judging basis of movement direction change and speed change in a certain time interval in space;the component values on the x, y and z axes respectively represent the rate change amplitude, the direction change amplitude and the position change amplitude of the node offset;
acquiring the time sequence evolution characteristic of the node voltage during the out-of-step oscillation of the power system by utilizing the vector offset characteristic space of the node voltage trajectory, wherein the time sequence evolution characteristic comprises the following steps:
collecting the time sequence information of the node voltage in the power grid through a wide area measurement system, acquiring the time sequence information of the node voltage in a certain time interval, and forming a voltage phasor track group UR of all nodes in a voltage complex space in the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; n is the number of all nodes of the power grid;
and according to the time sequence evolution characteristics of the node voltage, evaluating the similarity of the node voltage change based on a track clustering algorithm, and positioning one or more splitting sections by tracking the one or more splitting sections.
2. The method of claim 1, wherein the estimating the similarity of the node voltage changes based on a trajectory clustering algorithm according to the time-series evolution characteristics of the node voltages, and the locating of one or more splitting sections through tracking of the one or more splitting sections further comprises:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree among time sequence tracks and measuring the difference degree among different tracks; the distance between two discrete points at a certain time in the space is known as follows:
then generalizing to the inter-track distance within a certain time interval using distance D, D Traj Is calculated as follows:
D Traj =D 1 +D 2 +…+D t
in the formula D t The distance between two discrete points at the tth moment is represented, and t is the window length of the time sliding window;
the distance between any two tracks can be obtained by sliding the window length by one unit.
3. The method of claim 2, further comprising:
performing cluster expansion on the neighborhood track group, and offsetting tracks of each vector in the neighborhoodPerforming secondary traversal and label assignment by applying to the neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
4. A system for locating an active splitting section of a power grid based on voltage trajectory information, the system comprising:
a building unit, configured to build a vector offset eigenspace of the node voltage phase trajectory, further configured to:
with d l 、d θ 、d p Respectively serving as x, y and z coordinate axes of a three-dimensional space, and constructing a vector offset eigenspace psi of a node voltage phase track d Let a node start point V mi The coordinate is (d) lmi ,d θmi ,d pmi ) End point V (m+1)i The coordinate is (d) l(m+1)i ,d θ(m+1)i ,d p(m+1)i ) And the vector offset characteristic track of the vector between the starting point and the end point forming node i in a certain time interval is as follows:
wherein the content of the first and second substances,as a node at d l 、d θ 、d p Judging basis of motion direction change and speed change in a certain time interval in space;the component values on the x, y and z axes respectively represent the rate change amplitude, the direction change amplitude and the position change amplitude of the node offset;
the acquisition unit is used for acquiring the time sequence evolution characteristics of the node voltage during the out-of-step oscillation of the power system by utilizing the vector offset characteristic space of the node voltage trajectory, and is also used for:
collecting the time sequence information of the node voltage in the power grid through a wide area measurement system, acquiring the time sequence information of the node voltage in a certain time interval, and forming a voltage phasor track group UR of all nodes in a voltage complex space in the time interval Tm :
In the formula (I), the compound is shown in the specification,represents the m-th time interval T m Inner, voltage vector of node i; n is the number of all nodes of the power grid;
and the positioning unit is used for evaluating the similarity of the node voltage change based on a track clustering algorithm according to the time sequence evolution characteristics of the node voltage and positioning one or more splitting sections by tracking the one or more splitting sections.
5. The system of claim 4, wherein the positioning unit is configured to evaluate similarity of the node voltage changes based on a trajectory clustering algorithm according to a time-series evolution characteristic of the node voltage, and to position one or more separation sections by tracking the one or more separation sections, and further configured to:
selecting time sequence tracks in a preset time interval for clustering, defining the time sequence tracks as the sum of distances between points of different tracks under the same time scale, and recording the sum as D Traj The method is used for representing the similarity degree among time sequence tracks and measuring the difference degree among different tracks; the distance between two discrete points at a certain time in the known space is:
then generalizing to the inter-track distance within a certain time interval using distance D, D Traj Is calculated as follows:
D Traj =D 1 +D 2 +…+D t
in the formula D t The distance between two discrete points at the tth moment is represented, and t is the window length of the time sliding window;
the distance between any two tracks can be obtained by sliding the window length by one unit.
6. The system of claim 5, the positioning unit further to:
performing cluster expansion on the neighborhood track group, and shifting tracks of each vector in the neighborhoodPerforming secondary traversal and label endowing by aiming at neighborhoodInner track vector offset trackPerforming a second traversal to find an AND in the original clusterIndirectly similar and not yet assigned to the trajectory of the tag.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910760005.7A CN110690716B (en) | 2019-08-16 | 2019-08-16 | Method and system for positioning active splitting section of power grid based on voltage trajectory information |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910760005.7A CN110690716B (en) | 2019-08-16 | 2019-08-16 | Method and system for positioning active splitting section of power grid based on voltage trajectory information |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110690716A CN110690716A (en) | 2020-01-14 |
CN110690716B true CN110690716B (en) | 2022-09-27 |
Family
ID=69108339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910760005.7A Active CN110690716B (en) | 2019-08-16 | 2019-08-16 | Method and system for positioning active splitting section of power grid based on voltage trajectory information |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110690716B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111766462B (en) * | 2020-05-14 | 2023-03-10 | 中国计量大学 | Non-invasive load identification method based on V-I track |
CN111896869B (en) * | 2020-07-11 | 2023-03-24 | 西北工业大学 | Five-phase permanent magnet synchronous motor open-circuit fault diagnosis method |
CN116243097B (en) * | 2023-05-11 | 2023-08-15 | 新风光电子科技股份有限公司 | Electric energy quality detection method based on big data |
CN117559449B (en) * | 2024-01-12 | 2024-03-26 | 武汉华飞智能电气科技有限公司 | Power grid power transmission stability control method, system and storage medium |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7427815B1 (en) * | 2003-11-14 | 2008-09-23 | General Electric Company | Method, memory media and apparatus for detection of grid disconnect |
CN105226637B (en) * | 2015-09-01 | 2017-02-22 | 三峡大学 | Oscillation grouping identification method based low-frequency oscillation dispatch control method |
CN107482621B (en) * | 2017-08-02 | 2019-09-27 | 清华大学 | A kind of Transient Voltage Stability in Electric Power System appraisal procedure based on voltage sequential track |
CN107611996B (en) * | 2017-09-15 | 2020-04-21 | 华北电力大学 | Multi-frequency oscillation out-of-step oscillation center positioning and migration tracking system and method thereof |
-
2019
- 2019-08-16 CN CN201910760005.7A patent/CN110690716B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110690716A (en) | 2020-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110690716B (en) | Method and system for positioning active splitting section of power grid based on voltage trajectory information | |
CN110311376B (en) | Dynamic safety assessment comprehensive model and space-time visualization method for power system | |
CN101984533B (en) | Method for assessing power distribution reliability of large-scale medium-voltage distribution network based on modes | |
CN111628494B (en) | Low-voltage distribution network topology identification method and system based on logistic regression method | |
Zhang et al. | Optimal PMU placement evaluation for power system dynamic state estimation | |
CN104578053A (en) | Power system transient stability prediction method based on disturbance voltage trajectory cluster features | |
Bai et al. | Measurement-based frequency dynamic response estimation using geometric template matching and recurrent artificial neural network | |
Leng et al. | Short-term forecasting model of traffic flow based on GRNN | |
Romero et al. | First results in leak localization in water distribution networks using graph-based clustering and deep learning | |
CN103311960A (en) | Method for dividing coherent stability regions of strong-connection power grids | |
Kamaludin et al. | Cyber Physical System in the industry | |
CN112234627B (en) | Active splitting control method integrating off-line simulation and real-time information monitoring | |
CN104505827B (en) | Complicated electric power system closed loop control method based on response message | |
CN104901328A (en) | Multi-terminal flexible DC control mode automatic identification method based on complex control network | |
CN110675276A (en) | Method and system for inversion droop control of direct current power transmission system | |
Li et al. | Node localisation for wireless sensor networks in smart distribution automation | |
Hu et al. | Power grid's Intelligent Stability Analysis based on big data technology | |
CN107340454B (en) | Power system fault positioning analysis method based on RuLSIF variable point detection technology | |
CN108988319B (en) | Rapid emergency control method based on deep feedforward neural network and numerical integration sensitivity | |
Han et al. | Online Transfer Learning-based Method for Predicting Remaining Useful Life of Aero-engines | |
Wang et al. | Analysis and calculation of line loss data based on hybrid clustering | |
Hang et al. | Discussion on multi-sensor information fusion in greenhouse detection system | |
Wang et al. | Abnormal Diagnosis Method of Line Loss Rate in Low-voltage Station Area Based on Data-driven | |
Bai et al. | Abnormal Detection Scheme of Substation Equipment based on Intelligent Fusion Terminal | |
Rahman et al. | Distributed dynamic state estimation for smart grid transmission 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |