CN106845757A - A kind of electric network swim shifts evaluating severity method - Google Patents
A kind of electric network swim shifts evaluating severity method Download PDFInfo
- Publication number
- CN106845757A CN106845757A CN201611106863.2A CN201611106863A CN106845757A CN 106845757 A CN106845757 A CN 106845757A CN 201611106863 A CN201611106863 A CN 201611106863A CN 106845757 A CN106845757 A CN 106845757A
- Authority
- CN
- China
- Prior art keywords
- circuit
- power
- trend
- power flow
- flow transfer
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000012546 transfer Methods 0.000 claims abstract description 159
- 230000006378 damage Effects 0.000 claims abstract description 37
- 238000012544 monitoring process Methods 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 208000027418 Wounds and injury Diseases 0.000 claims abstract description 9
- 208000014674 injury Diseases 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 208000015181 infectious disease Diseases 0.000 claims description 55
- 230000007257 malfunction Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 238000010845 search algorithm Methods 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 230000000630 rising effect Effects 0.000 claims 1
- 230000033228 biological regulation Effects 0.000 abstract description 7
- 238000011156 evaluation Methods 0.000 abstract description 6
- 230000008569 process Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 230000001360 synchronised effect Effects 0.000 abstract description 4
- 238000005070 sampling Methods 0.000 abstract description 3
- 238000000205 computational method Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 125000002015 acyclic group Chemical group 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Emergency Protection Circuit Devices (AREA)
Abstract
A kind of electric network swim shifts evaluating severity method.The method is using the synchronous real-time dynamic power delta data of the high sampling rate wide area that wide area measurement system is obtained, based on circuit power swing information in power network day-to-day operation, using statistical technique, the extent of injury of each interline power flow transfer relationship and its safety and stability to power network in power network is calculated, and the circuit that may select from out generation line disconnection and have serious harm to power network safety operation carries out key monitoring.The method is not based on Load flow calculation, and rely only on power grid operation personnel can direct access measured data, it is contemplated that various chain regulation devices and dynamic process are to the resultant effect of power flow transfer, so as to realize the Dynamic Comprehensive Evaluation to the power flow transfer extent of injury.
Description
Technical field
The application belongs to power system safety and stability evaluation areas, and more particularly to a kind of safety caused to power flow transfer is steady
Determine the appraisal procedure of problem.
Background technology
During the failure such as be short-circuited in power network, the stabilization control device of relay protection or part is generally by cut-offfing fault wire
Road protects the local devices or partial electric grid not compromised.However, cut-offfing for circuit, often causes via the line transmission
Power is transferred on the All other routes being connected, if the trend for receiving the circuit of transfer trend exceedes the safety and stability of its permission
Capacity limit value, then can cause further line disconnection, so as to cascading failure occur, cause the expansion of power grid accident.In recent years
The power grid accident on a large scale for occurring both at home and abroad is almost directed to local line and cut-offs the power flow transfer for causing, and causes accident scope to expand
Greatly.Therefore, the extent of injury that the power flow transfer that Pre-Evaluation line disconnection causes is caused to power network is to prevention and the expansion of control accident
It is big significant.
According to existing document, two classes are mainly to the appraisal procedure of electric network influencing degree on power flow transfer, a class is to borrow
Help power flow equation, by calculating variable between Jacobian matrix i.e. sensitivity coefficient calculate other lines that line disconnection causes
The change of power flow on road；Another kind of is to utilize emulation mode, by artificially cut-offfing a certain circuit, observes it and is caused in All other routes
Influence.But this two classes method is respectively provided with following shortcoming：
(1) Jacobian matrix is calculated and simulation modeling is calculated and is required to by electric network model parameter and network real-time topology,
When these data are difficult to obtain, or when differing larger with actual value, it is impossible to carry out corresponding mathematical modeling and calculating, or
Result of calculation is inaccurate.
(2) actual electric network is due to when failure occurs power flow transfer, may relate to the interlocking action of regulation device, and unit
The dynamic regulation of equipment, such as connection of steady control cuts operation, the dynamic regulation of AGC and AVC, power electronic equipment it is idle and active
Quick support etc..These chain dynamic behaviours are known and modeled all to have difficulties comprehensively.
Additionally, wide area measurement system obtains extensively should in Provincial and above power network since two thousand and ten
With.Each wide area measurement system is generally with the speed of 50 frames/second or 25 frames/second from each measurement substation high speed real-time collecting voltage
The metric data such as electric current phasor, frequency, power, and each data all have measurement markers of the high precision up to 1 microsecond, therefore can be with
Realize the real-time synchronization observation to the whole network dynamic process, target when at this moment conventional several seconds 1 measurement points, and data are not measured
SCADA to be had observation advantage.
It should also be mentioned that power flow transfer takes place mostly in the network of loop-net operation, and in order to avoid electromagnetic looped network, only
High pressure main grid structure operationally has a looped network, 220kV be following power network generally using looped network construction, the principle of open loop operation, because
This also all has been equipped with wide area and mutually measures it may be said that the power network with extensive loop-net operation is all provincial and above power network at present
Amount system.
The present invention proposes one on wide area measurement system platform, based in power network day-to-day operation in view of the foregoing
Line power fluctuation information, using statistical technique, calculates each Line Flow transfer case and its peace to power network in power network
The method of the extent of injury of full stabilization, and may select from out generation line disconnection and have serious harm to power network safety operation
Circuit carry out key monitoring.The method relies only on the measured data that power grid operation personnel can directly grasp, and can comprehensively examine
Consider the resultant effect of various chain regulation devices and dynamic process, therefore the present invention can be realized to the power flow transfer extent of injury
Objective, Real-Time Evaluation.
The content of the invention
To solve above mentioned problem of the prior art, the invention discloses a kind of electric network swim transfer evaluating severity side
Method, the real-time dynamic power delta data of the high sampling rate wide area synchronization obtained using wide area measurement system platform, base
The circuit power swing information in power network day-to-day operation, using statistical technique, each interline power flow transfer is closed in calculating power network
The extent of injury of system and its safety and stability to power network, and generation line disconnection is may select from out to electricity net safety stable fortune
The circuit that row has serious harm carries out key monitoring.
The present invention specifically uses following technical scheme：
A kind of electric network swim shifts evaluating severity method, it is characterised in that the described method comprises the following steps：
Step 1：In actual electric network, recognize all with the circuit of power flow transfer relation；
Step 2：Statistics produces circuit and is transferred to active power lifting number of times on circuit with trend with period trend, recognizes trend
Transfer association factor；
Step 3：Circuit is produced according to trend and is transferred to the statistics of line power variable quantity with trend, calculate credible power flow transfer
The factor；
Step 4：Calculate the active power transfer amount caused when actual track cut-offs；
Step 5：According to the influence after line disconnection to All other routes Capacity Margin, the trend turn that line disconnection causes is calculated
Move direct density of infection index；
Step 6：According to whether occurring that chain capacity is out-of-limit, the chain density of infection index of power flow transfer of line disconnection is calculated；
Step 7：According to the statistics to historical failure in power network, the power flow transfer harm risk that calculating line disconnection causes refers to
Mark；
Step 8：Each circuit is carried out by the chain density of infection index of power flow transfer and power flow transfer harm risk indicator respectively
Sequence, and different Control Measures are given to each circuit according to ranking results.
The present invention further includes following preferred scheme：
In step 1, all with the circuit of power flow transfer relation, its content includes for identification：
(1.1) transformer between the circuit and the two voltage class of highest in power network and secondary voltage levels is turned
The line segment in network is turned to, searched network is formed；
(1.2) whole mesh are found out in searched network using any Mesh Search Algorithm；
(1.3) using the circuit on all mesh as the statistics monitoring wire of power flow transfer, each circuit of stream calculation is taken turns
Power flow transfer density of infection, when select wherein 1 circuit as evaluated power flow transfer density of infection circuit when, the circuit is called
Trend produces circuit, and All other routes are transferred to circuit as trend.
In step 2, herein below is specifically included：
(2.1) there is the active power of each circuit of power flow transfer relation, in setting time section Δ t interior lines in monitoring power network
Road active power declines change and exceedes threshold value Δ PthWhen, judge the line power bust, it is that trend produces circuit；Finding out circuit has
The circuit i of work(power bust, i.e. trend produce circuit i；
(2.2) in the setting time section Δ t that above-mentioned trend produces the decline of circuit i active power, accordingly search at this
There is the circuit j of active power rise in time period；
(2.3) the power flow transfer association factor that circuit i line relateds are produced with trend is calculated：In statistical time range, if damp
The active reduction of circuit i of circulating out is that trend produces number of times for NDi, the active power rise times of the circuit j that correspondence is detected
It is NRi_j, calculateThen trend produces power flow transfer association factor Rs of the circuit i to circuit ji-jIt is as follows：
Wherein, the setting time section Δ t is 1 second；
The statistical time range is 1 day, 1 week or January.
Threshold value Δ PthPreset according to voltage class or network load situation.For the power network of 220kV and 500kV, institute
State threshold value Δ PthPreferably 50MW.
In step 3, calculating credible flow transferring relativity factor includes herein below：
(3.1) the active power drop-out value Δ P of circuit i n-ths is produced according to trendi_drop_nWhat is caused on circuit j has
Work(power rise value Δ Pj_rise_n(≠ 0), calculate n-th actual measurement trend produce power flow transfer coefficient from circuit i to circuit j i.e.
Active power transfer ratio ki_j_n, wherein, active power drop-out value Δ Pi_drop_nMore than threshold value Δ Pth, calculate according to this each
Secondary trend produces power flow transfer coefficient k when circuit i power drops and circuit j power risesi_j_n：
(3.2) calculate trend and produce average power flow transfer coefficient ks of the circuit i to circuit ji_j：
Wherein, NRi_jIt is the active power of the circuit j that trend produces circuit i power drops and causes in statistical time range
Rise times；
(3.3) trend is calculated according to the following formula produce credible flow transferring relativity factor Ks of the circuit i to circuit ji_j：
Ki_j=Ri_j·ki_j
Wherein, Ri-jFor trend produces power flow transfer association factors of the circuit i to circuit j.
In step 4, when trend produces circuit i, and line disconnection occurs, according to credible flow transferring relativity factor Ki_jCalculate line
The i.e. active power transfer amount of active power incrementss caused on the j of road is：
ΔPj_rise_F=Ki_j·ΔPi_drop_C；
Now the prediction trend of circuit j is：
Pj_F=Δ Pj_rise_F+Pj_C；
Wherein, Δ Pj_rise_FIt is the active power incrementss caused on circuit j, Δ Pi_drop_CRepresent that circuit i occurs circuit
Active power drop-out value when cut-offfing, Ki_jRepresent that trend produces credible flow transferring relativity factors of the circuit i to circuit j, Pj_CRepresent
Trend produces the trend i.e. active power value of circuit j before circuit i open circuits, Pj_FRepresent that trend produces circuit j after circuit i open circuits
Prediction trend.
In steps of 5, the direct density of infection index of power flow transfer is calculated in accordance with the following methods：
(5.1) it is calculated as follows the nargin coefficient D that the circuit j of power rise in the case of circuit i cut-offs is produced in trendi_j,
Work as Di_j1 is more less than, then the Capacity Margin of circuit j is more；Work as Di_jWhen >=1, circuit j effective power flows meet or exceed capacity
Limit, it is believed that circuit j effective power flows are out-of-limit；
Wherein, Pj_FThe prediction trend of circuit j, P after the i open circuits of expression circuitj_LimitRepresent the capacity limit of circuit j；
(5.2) cause m bar circuit effective power flows out-of-limit if trend produces circuit i and cut-offs, i.e., the D of corresponding m bars circuiti_j≥
1, then trend produce circuit i and cut-off the direct density of infection index D of the power flow transfer for causingi_DFor：
Di_D=m
(5.3) circuit i is produced for cut-offfing the trend for not causing any circuit effective power flow out-of-limit, its power flow transfer is direct
Density of infection index Di_DFor：
Di_D=max { Di_j}。
In step 6, the chain density of infection index of power flow transfer is calculated in accordance with the following methods：
(6.1) if trend is produced after circuit i cut-offs, directly cause the effective power flow of m bar circuits out-of-limit, this m bar circuit meeting
Further cut-off because of overload, and further cause power flow transfer；If in the m bars overload circuit that circuit i causes, thering are M bars to enter one
Step is cut-off All other routes effective power flow can be caused out-of-limit, then trend produces circuit i and cut-offs the chain density of infection of the power flow transfer for causing
Index DT_iFor：
DT_i=m+10M；
(6.2) circuit i is produced for trend, only directly causes part circuit effective power flow out-of-limit if it cut-offs, and these
Even if the out-of-limit circuit of effective power flow cut-offs and also no longer further result in that other circuit effective power flows are out-of-limit, then trend produces circuit
The chain density of infection index of power flow transfer of i is set to：
DT_i=Di_D；
(6.3) circuit i is produced for the trend for not causing any Line Flow out-of-limit, the chain density of infection of its power flow transfer refers to
It is demarcated as：
DT_i=Di_D=max { Di_j＜ 1.
In step 7, trend produces the power flow transfer harm risk indicator that circuit causes and calculates by the following method：
(7.1) according to historical statistics, the probability of malfunction that trend produces circuit is obtained, for the failure that trend produces circuit i
Probability Gfault_iIt is defined as follows：
(7.2) the chain density of infection index of power flow transfer and probability of malfunction when producing line disconnection according to each trend, draw
The power flow transfer harm risk indicator of each circuit, the power flow transfer harm risk indicator F of circuit iT_iIt is expressed as follows：
FT_i=Gfault_i·DT_i
Wherein, DT_iThe chain density of infection index of power flow transfer for causing is cut-off for trend produces circuit i.
In step 8, the chain density of infection index of power flow transfer and power flow transfer harm risk indicator are pressed respectively to each circuit
It is ranked up：
It is more than or equal to 4 to the chain density of infection index of power flow transfer or comes first 3, and its power flow transfer harm risk refers to
Mark comes the circuit of top 10, takes Control Measure；
4 are more than or equal to the chain density of infection index of power flow transfer, and power flow transfer risk indicator comes 10 lines afterwards
Road, takes emphasis to monitor measure；
It is more than 1 for the chain density of infection index of power flow transfer, but circuit less than 4 is classified as time emphasis monitoring wire；
Circuit for the chain density of infection of power flow transfer less than 1 uses common monitoring measure.
Implementation of the invention would help power grid operation personnel be independent of the concrete model of power network, parameter, topology with
And in the case of protection control device acting characteristic, rely only on statistics letter of the wide area measurement system in power network day-to-day operation
Breath, just can carry out quantitative assessment, and find out it to the power flow transfer density of infection and risk that cause after each line disconnection in power network
Middle have most serious harm or the circuit of risk the prevention and control and monitoring of emphasis carried out as critical circuits.Relative to traditional
Power flow transfer computational methods based on Load flow calculation, the Statistics-Based Method in the present invention can contemplate various chain regulation and control dresses
Put and dynamic process is to the resultant effect of power flow transfer, so as to realize the Dynamic Comprehensive Evaluation to the power flow transfer extent of injury.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of power flow transfer evaluating severity method of the present invention.
Specific embodiment
Technical scheme is described in further detail with reference to Figure of description and specific embodiment.
The present invention is using the synchronous real-time dynamic power of the high sampling rate wide area that wide area measurement system platform is obtained
Delta data, based on circuit power swing information in power network day-to-day operation, using statistical technique, calculates each circuit tide in power network
The extent of injury of stream transfer case and its safety and stability to power network, and generation line disconnection is may select from out to power network peace
The circuit that full stable operation has serious harm carries out key monitoring.In actual electric network, opened based on method proposed by the invention
The power flow transfer evaluating severity program of hair will run on the wide area measurement of provincial, branch center level and national grid dispatching center
In the senior application service of system WAMS.WAMS obtains control centre and is managed with the speed of 50 frames/second or 25 frames/second
The synchronous metric data such as voltage phasor, electric current phasor, power, the frequency of each circuit measured by PMU on linchpin mains network, is adjusting
The WAMS main websites at degree center, according to the high-precise synchronization markers that data possess, are synchronized after these information are received to data
Alignment, and it is stored in time series real-time database and history library.Meanwhile, WAMS main website is stored with the real-time network of power network
Topological model, the metric data in database is associated with each line facility, so that each senior application of WAMS main websites can be with
Carry out the electrical network analysis based on network topological information.In this application platform, according to the present invention exploitation program can start into
Row power flow transfer evaluating severity.Due to needing the statistics based on a period of time, tide between relatively reliable circuit can be just obtained
Stream transfer relationship, therefore typically at least to have the historical data of 1 hour, as the basis of statistical analysis.But in order to avoid
In measurement period there is big change in network structure, thus it is excessive to power flow transfer relationship affect, therefore measurement period is also unsuitable
Obtain long, be advisable with being no longer than for 1 week.It is noted that for the elements such as circuit probability of malfunction count when
It is long, can not be constrained by this measurement period.
Based on above-mentioned actual application environment, the power flow transfer density of infection run on the senior application server of WAMS platforms is commented
Estimating program can circulate execution following steps (flow chart is shown in accompanying drawing 1), realize the power flow transfer danger to each circuit of administered power network
Evil degree is estimated, and cycle period can be elected as 5 minutes.
Step 1：Using the identification of mesh search method all with the circuit of power flow transfer relation；
Because power flow transfer is occurred over just on looped network, and in actual electric network, looped network is primarily present in the power network highest
Or on the major network of secondary voltage levels (when highest voltage level is very weak, secondary voltage levels allow looped network), therefore
Need to find out all of looped network path on two voltage class power networks of power network highest, (by the side for counting only on loop grid
Method) identification power flow transfer relation.Circuit on acyclic networking footpath is charge circuit or equivalent generating machine circuit, in loop grid hair
After life is cut-off, its trend is it can also happen that change, but this change is caused due to the change of system side equivalent impedance, not very
Make power flow transfer, therefore do not considered when power flow transfer path is searched.The algorithm of the circuit on searching looped network path is as follows：
(1.1) transformer between the circuit and the two voltage class of highest in power network and secondary voltage levels is turned
Turn to line segment (for example：500kV circuits, 220kV circuits, 500/220 transformer are converted into line segment), form searched network N.
(1.2) whole mesh are found out in searched network N using any Mesh Search Algorithm.Paper can for example be used
《Design of the mesh automatic search algorithm in Hydropower Simulation and realization》Northeastern University's journal (natural science edition) September in 2008
The phase of volume 29 the 9th, find out all mesh.
(1.3) using the circuit on all mesh as power flow transfer statistics monitoring wire.Each line of stream calculation can be taken turns
The power flow transfer density of infection on road.When circuit of wherein 1 circuit of selection as evaluated power flow transfer density of infection, claim the circuit
For trend produces circuit, All other routes are transferred to circuit as trend.
Step 2：Statistics produces circuit and is transferred to active power lifting number of times on circuit with trend with period trend, recognizes trend
Transfer association factor；
(2.1) there is the active power of each circuit of power flow transfer relation, in setting short time Δ t in monitoring power network
(being usually taken to be 1 second) interior circuit power drop change exceedes threshold value Δ Pth(can be changed according to voltage class or network load situation,
50MW is usually taken to be in the power network of 220kV and 500kV) when, judge the line power bust, it is that trend produces circuit；Find out line
The circuit i of road active power bust, i.e. trend produce circuit i；
(2.2) in the time period Δ t that above-mentioned trend produces the decline of circuit i active power, accordingly search in the period
There is the circuit j of active power rise；
(2.3) the power flow transfer association factor that circuit i line relateds are produced with trend is calculated：In the statistical time range of setting
(such as 1 day, 1 week or January, if the active reduction number of times of circuit i is NDi, the rise times of circuit j that correspondence is detected are
NRi_j, calculateIf ri_j>0.8, show that the trend of circuit i is certain to turn to circuit j after hindered, almost
Move, therefore from from the point of view of guarding, take trend and produce power flow transfer association factor Rs of the circuit i to circuit ji_jIt is 1；If ri_j<
0.2, show almost there is no power flow transfer relation between circuit i and circuit j, the increase of the circuit j power for observing is not from
Circuit i, takes trend and produces power flow transfer association factor Rs of the circuit i to circuit ji_j=0；Trend is taken in the case of other and produces circuit
Power flow transfer association factor R from i to circuit ji_j=ri_j, in sum trend produce power flow transfer from circuit i to circuit j associate
Factor Ri-jIt is as follows：
Step 3：Circuit is produced according to trend and is transferred to the statistics of line power variable quantity with trend, calculate credible power flow transfer
The factor；
(3.1) active power for producing circuit i n-ths according to trend declines Δ Pi_drop_nThe power caused on circuit j
Rise Δ Pj_rise_n(≠ 0), the trend for calculating n-th actual measurement produces power flow transfer coefficient i.e. active power of the circuit i to circuit j
Transfer ratio ki_j_n, wherein, active power drop-out value Δ Pi_drop_nMore than threshold value Δ Pth, each secondary line i work(is calculated according to this
Power flow transfer coefficient k when rate jump drops and circuit j power rises toi_j_n：
(3.2) calculate trend and produce the average power flow transfer coefficient k that circuit i thinks circuit ji_j：
(3.3) trend is calculated according to the following formula produce credible flow transferring relativity factor Ks of the circuit i to circuit ji_j：
Ki_j=Ri_j·ki_j
Step 4：Calculate the active power transfer amount caused when actual track cut-offs:
When trend produces circuit i, and line disconnection occurs, according to credible flow transferring relativity factor Ki_jCalculate and cause on circuit j
Active power incrementss be：
ΔPj_rise_F=Ki_j·ΔPi_drop_C
At this moment, the prediction trend of circuit j is：
Pj_F=Δ Pj_rise_F+Pj_C
Wherein, Δ Pi_drop_CRepresent active power drop-out value during circuit i generation line disconnections, Ki_jRepresent that trend is produced
Credible flow transferring relativity factor from circuit i to circuit j, Pj_CThe trend for representing the preceding circuit j of circuit i open circuits is active power value, Pj_F
Represent that trend produces the prediction trend of circuit j after circuit i open circuits.
Step 5：According to the influence after line disconnection to All other routes Capacity Margin, the trend turn that line disconnection causes is calculated
Move direct density of infection index；
(5.1) it is calculated as follows the nargin coefficient D that circuit j impacted in the case of circuit i cut-offs is produced in trendi_j, when
Di_j1 is more less than, then the Capacity Margin of circuit is more；Work as Di_jWhen >=1, circuit effective power flow meets or exceeds capacity limit, it is believed that
Circuit j effective power flows are out-of-limit；
(5.2) if trend produces circuit line disconnection and causes m bar circuit effective power flows out-of-limit, i.e., the D of corresponding m bars circuiti_j
≥1.It is D that then trend produce circuit i to cut-off the direct density of infection index of the power flow transfer for causingi_D：
Di_D=m
(5.3) circuit i is produced for cut-offfing the trend for not causing any circuit effective power flow out-of-limit, its power flow transfer is direct
Density of infection index is less than 1, is taken as the Capacity Margin index of the minimum circuit of Capacity Margin, i.e.,
Di_D=max { Di_j}。
Step 6：According to whether occurring that chain capacity is out-of-limit, the chain density of infection index of power flow transfer of line disconnection is calculated；
(6.1) if trend is produced after circuit i cut-offs, directly cause the effective power flow of m bar circuits out-of-limit, in this m bar circuit
Certain circuit j trend it is out-of-limit and cut-off the trend transition that can further cause All other routes, according to circuit j to circuit l's
Flow transferring relativity factor Kj_l, can obtain the trend after the circuit l transition caused after j cut-offs is：
Pi_j_l_F=Kj_l·Pj_F+Pl_C
Predict it is continuous cut-off circuit i and circuit j after, the trend of circuit l is Pi_j_l_FIf it exceedes circuit l's
Capacity limits Pl_limit, then cut-offfing for circuit l can further be caused.According to the chain method cut-off of above-mentioned prediction circuit, if i causes
M bars overload circuit in, having M bars further to cut-off can cause All other routes effective power flow out-of-limit, then trend produces circuit i and cut-offs
The chain density of infection index circuit of power flow transfer for causing is：
DT_i=m+10M
(6.2) circuit i is produced for trend, only directly causes part circuit effective power flow out-of-limit if it cut-offs, and these
Even if the out-of-limit circuit of effective power flow cut-offs and also no longer further result in that other circuit effective power flows are out-of-limit, then trend produces circuit
The chain density of infection index of power flow transfer of i is set to：：
DT_i=Di_D；
(6.3) circuit i is produced for the trend for not causing any Line Flow out-of-limit, the chain density of infection of its power flow transfer refers to
It is demarcated as：
DT_i=Di_D=max { Di_j＜ 1.
DT_iBigger explanation, the harm of the line disconnection is bigger.
Step 7：According to the statistics to historical failure in power network, the power flow transfer harm risk that calculating line disconnection causes refers to
Mark；
(7.1) according to historical statistics, the probability of malfunction that trend produces circuit is obtained, for the probability of malfunction of circuit i
Gfault_iIt is defined as follows：
(7.2) the chain density of infection index of power flow transfer and probability of malfunction when producing line disconnection according to each trend, draw
The power flow transfer harm risk indicator of each circuit, the power flow transfer harm risk indicator F of circuit iT_iIt is expressed as follows：
FT_i=Gfault_i·DT_i。
Step 8：Each circuit is ranked up by the chain density of infection index of power flow transfer and power flow transfer harm risk indicator,
And different Control Measures are given to each circuit according to ranking results.
According to abovementioned steps, the power flow transfer that can calculate the circuit that looped network or mesh in current operation power network are related to is endangered
Evil degree index and harm risk indicator, then carry out following sequence：
(1) the chain density of infection index of power flow transfer according to each circuit is ranked up to circuit；
(2) the power flow transfer harm risk indicator according to each circuit is ranked up to circuit.
It is more than or equal to 4 to the chain density of infection index of power flow transfer or comes first 3, and its power flow transfer harm risk refers to
Mark comes the circuit of top 10, takes Control Measure；To the chain density of infection index of power flow transfer more than or equal to 4, and trend
The index that shifts risk comes 10 circuits afterwards, takes emphasis to monitor measure；It is big for the chain density of infection index of power flow transfer
In 1, but circuit less than 4 can be classified as time emphasis monitoring wire；For circuit of the chain density of infection of power flow transfer less than 1 using general
Logical monitoring measure.
Relative to traditional power flow transfer computational methods based on Load flow calculation, the Statistics-Based Method in the present invention can
To consider the resultant effect of various chain regulation devices and dynamic process to power flow transfer, so as to realize endangering power flow transfer
The Dynamic Comprehensive Evaluation of degree.
Applicant describes in detail and describes with reference to Figure of description to embodiments of the invention, but this area skill
It should be understood that above example is only the preferred embodiments of the invention, explanation in detail is intended merely to help reader art personnel
More fully understand that the present invention is spiritual, and not limiting the scope of the invention, conversely, any based on invention of the invention essence
Any improvement or modification that god is made should all fall within the scope and spirit of the invention.
Claims (13)
1. a kind of electric network swim shifts evaluating severity method, it is characterised in that：
The real-time dynamic power delta data of the wide area synchronization obtained using wide area measurement system, based on the daily fortune of power network
Circuit power swing information in row, using statistical technique, calculate in power network each interline power flow transfer relationship and its to electricity
The extent of injury of the safety and stability of net, and may select from out generation line disconnection and have serious harm to power network safety operation
Circuit carry out key monitoring.
2. a kind of electric network swim shifts evaluating severity method, it is characterised in that the described method comprises the following steps：
Step 1：In actual electric network, recognize all with the circuit of power flow transfer relation；
Step 2：Statistics produces circuit and is transferred to active power lifting number of times on circuit with trend with period trend, recognizes power flow transfer
Association factor；
Step 3：Circuit and trend are produced according to trend and are transferred to the statistics of line power variable quantity, calculate credible power flow transfer because
Son；
Step 4：Calculate the active power transfer amount caused when actual track cut-offs；
Step 5：According to the influence after line disconnection to All other routes Capacity Margin, the power flow transfer that calculating line disconnection causes is straight
Connect density of infection index；
Step 6：According to whether occurring that chain capacity is out-of-limit, the chain density of infection index of power flow transfer of line disconnection is calculated；
Step 7：According to the statistics to historical failure in power network, the power flow transfer harm risk indicator that line disconnection causes is calculated；
Step 8：Each circuit is ranked up by the chain density of infection index of power flow transfer and power flow transfer harm risk indicator respectively,
And different Control Measures are given to each circuit according to ranking results.
3. electric network swim according to claim 2 shifts evaluating severity method, it is characterised in that：
In step 1, all with the circuit of power flow transfer relation, its content includes for identification：
(1.1) transformer between the circuit and the two voltage class of highest in power network and secondary voltage levels is converted into
Line segment in network, forms searched network；
(1.2) whole mesh are found out in searched network using any Mesh Search Algorithm；
(1.3) using the circuit on all mesh as the statistics monitoring wire of power flow transfer, each tide of circuit of stream calculation is taken turns
Circulation moves density of infection, and when circuit of wherein 1 circuit of selection as evaluated power flow transfer density of infection, the circuit is called trend
Circuit is produced, All other routes are transferred to circuit as trend.
4. electric network swim according to claim 2 shifts evaluating severity method, it is characterised in that：
In step 2, herein below is specifically included：
(2.1) there is the active power of each circuit of power flow transfer relation, circuit has in setting time section Δ t in monitoring power network
The change of work(power drop exceedes threshold value Δ PthWhen, judge the line power bust, it is that trend produces circuit；Find out circuit wattful power
The circuit i of rate bust, i.e. trend produce circuit i；
(2.2) in the setting time section Δ t that above-mentioned trend produces the decline of circuit i active power, accordingly search in the time
There is the circuit j of active power rise in section；
(2.3) the power flow transfer association factor that circuit i line relateds are produced with trend is calculated：In statistical time range, if trend turns
The active reduction i.e. trend for going out circuit i produces number of times for NDi, the active power rise times of circuit j that correspondence is detected are
NRi_j, calculateThen trend produces power flow transfer association factor Rs of the circuit i to circuit ji-jIt is as follows：
5. electric network swim according to claim 4 shifts evaluating severity method, it is characterised in that：
The setting time section Δ t is 1 second；
The statistical time range is 1 day, 1 week or January.
6. the electric network swim according to claim 4 or 5 shifts evaluating severity method, it is characterised in that：
Threshold value Δ PthPreset according to voltage class or network load situation.
7. electric network swim according to claim 6 shifts evaluating severity method, it is characterised in that：
For the power network of 220kV and 500kV, the threshold value Δ PthIt is taken as 50MW.
8. the electric network swim according to claim 2 or 4 shifts evaluating severity method, it is characterised in that：
In step 3, calculating credible flow transferring relativity factor includes herein below：
(3.1) the active power drop-out value Δ P of circuit i n-ths is produced according to trendi_drop_nThe wattful power caused on circuit j
Rate rising value Δ Pj_rise_n(≠ 0), it is i.e. active to the power flow transfer coefficient of circuit j that the trend that calculating n-th is surveyed produces circuit i
Power transfer ratio ki_j_n, wherein, active power drop-out value Δ Pi_drop_nMore than threshold value Δ Pth, tide each time is calculated according to this
Power flow transfer coefficient k when the circuit i power drops that circulate out and circuit j power risesi_j_n：
(3.2) calculate trend and produce average power flow transfer coefficient ks of the circuit i to circuit ji_j：
Wherein, NRi_jIt is that in statistical time range, the active power of the circuit j that trend produces circuit i power drops and causes rises
Number of times；
(3.3) trend is calculated according to the following formula produce credible flow transferring relativity factor Ks of the circuit i to circuit ji_j：
Ki_j=Ri_j·ki_j
Wherein, Ri-jFor trend produces power flow transfer association factors of the circuit i to circuit j.
9. the electric network swim according to claim 4 or 8 shifts evaluating severity method, it is characterised in that：
In step 4, when trend produces circuit i, and line disconnection occurs, according to credible flow transferring relativity factor Ki_jCalculate circuit j
On the i.e. active power transfer amount of active power incrementss that causes be：
ΔPj_rise_F=Ki_j·ΔPi_drop_C；
Now the prediction trend of circuit j is：
Pj_F=Δ Pj_rise_F+Pj_C；
Wherein, Δ Pj_rise_FIt is the active power incrementss caused on circuit j, Δ Pi_drop_CRepresent that circuit i occurs line disconnection
When active power drop-out value, Ki_jRepresent that trend produces credible flow transferring relativity factors of the circuit i to circuit j, Pj_CRepresent trend
Produce the trend i.e. active power value of the preceding circuit j of circuit i open circuits, Pj_FRepresent that trend produces the prediction of circuit j after circuit i open circuits
Trend.
10. electric network swim according to claim 9 shifts evaluating severity method, it is characterised in that：
In steps of 5, the direct density of infection index of power flow transfer is calculated in accordance with the following methods：
(5.1) it is calculated as follows the nargin coefficient D that the circuit j of power rise in the case of circuit i cut-offs is produced in trendi_j, when
Di_j1 is more less than, then the Capacity Margin of circuit j is more；Work as Di_jWhen >=1, circuit j effective power flows meet or exceed capacity limit,
Think that circuit j effective power flows are out-of-limit；
Wherein, Pj_FThe prediction trend of circuit j, P after the i open circuits of expression circuitj_LimitRepresent the capacity limit of circuit j；
(5.2) cause m bar circuit effective power flows out-of-limit if trend produces circuit i and cut-offs, i.e., the D of corresponding m bars circuiti_j>=1, then
Trend produces circuit i and cut-offs the direct density of infection index D of the power flow transfer for causingi_DFor：
Di_D=m
(5.3) circuit i is produced for cut-offfing the trend for not causing any circuit effective power flow out-of-limit, its power flow transfer directly endangers
Degree index Di_DFor：
Di_D=max { Di_j}。
11. electric network swim transfer evaluating severity methods according to claim 10, it is characterised in that：
In step 6, the chain density of infection index of power flow transfer is calculated in accordance with the following methods：
(6.1) if trend is produced after circuit i cut-offs, directly cause the effective power flow of m bar circuits out-of-limit, this m bars circuit can be because of mistake
Carry and further cut-off, and further cause power flow transfer；If in the m bars overload circuit that circuit i causes, thering are M bars further to open
Disconnected that All other routes effective power flow can be caused out-of-limit, then trend produces circuit i and cut-offs the chain density of infection index of the power flow transfer for causing
DT_iFor：
DT_i=m+10M；
(6.2) circuit i is produced for trend, only directly causes part circuit effective power flow out-of-limit if it cut-offs, and these are active
Even if the out-of-limit circuit of trend cut-offs and also no longer further result in that other circuit effective power flows are out-of-limit, then trend produces circuit i's
The chain density of infection index of power flow transfer is set to：
DT_i=Di_D；
(6.3) circuit i is produced for the trend for not causing any Line Flow out-of-limit, the chain density of infection index of its power flow transfer is determined
For：
DT_i=Di_D=max { Di_j＜ 1.
The 12. electric network swim transfer evaluating severity method according to claim 2 or 11, it is characterised in that：
In step 7, trend produces the power flow transfer harm risk indicator that circuit causes and calculates by the following method：
(7.1) according to historical statistics, the probability of malfunction that trend produces circuit is obtained, for the probability of malfunction that trend produces circuit i
Gfault_iIt is defined as follows：
(7.2) the chain density of infection index of power flow transfer and probability of malfunction when producing line disconnection according to each trend, draw each line
The power flow transfer harm risk indicator on road, the power flow transfer harm risk indicator F of circuit iT_iIt is expressed as follows：
FT_i=Gfault_i·DT_i
Wherein, DT_iThe chain density of infection index of power flow transfer for causing is cut-off for trend produces circuit i.
13. electric network swim transfer evaluating severity methods according to claim 2, it is characterised in that：
In step 8, each circuit is carried out by the chain density of infection index of power flow transfer and power flow transfer harm risk indicator respectively
Sequence：
It is more than or equal to 4 to the chain density of infection index of power flow transfer or comes first 3, and its power flow transfer harm risk indicator row
In the circuit of top 10, Control Measure is taken；
4 are more than or equal to the chain density of infection index of power flow transfer, and power flow transfer risk indicator comes 10 circuits afterwards, adopts
Take emphasis monitoring measure；
It is more than 1 for the chain density of infection index of power flow transfer, but circuit less than 4 is classified as time emphasis monitoring wire；
Circuit for the chain density of infection of power flow transfer less than 1 uses common monitoring measure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611106863.2A CN106845757B (en) | 2016-12-06 | 2016-12-06 | Power grid power flow transfer hazard degree evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611106863.2A CN106845757B (en) | 2016-12-06 | 2016-12-06 | Power grid power flow transfer hazard degree evaluation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106845757A true CN106845757A (en) | 2017-06-13 |
| CN106845757B CN106845757B (en) | 2020-09-29 |
Family
ID=59145465
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611106863.2A Active CN106845757B (en) | 2016-12-06 | 2016-12-06 | Power grid power flow transfer hazard degree evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106845757B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107316135A (en) * | 2017-06-16 | 2017-11-03 | 广西电网有限责任公司 | A kind of risk quantification appraisal procedure of switch-oriented and switch operation process |
| CN109687438A (en) * | 2018-12-19 | 2019-04-26 | 西南交通大学 | It is a kind of meter and high-speed rail impact load effect under power grid vulnerable line discrimination method |
| CN110265998A (en) * | 2019-05-31 | 2019-09-20 | 国网辽宁省电力有限公司 | A kind of method and device that stabilized control system strategy definite value is checked |
| CN111402075B (en) * | 2020-02-26 | 2023-04-07 | 国网浙江省电力有限公司经济技术研究院 | Power transmission network operation flexibility evaluation method considering line distribution factors |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101615786A (en) * | 2009-07-24 | 2009-12-30 | 清华大学 | Grid system protection method based on local information of transformer substation |
| CN101895124A (en) * | 2010-05-27 | 2010-11-24 | 中国电力科学研究院 | Coupling assessment method for judging mutual influence of alternating current-direct current parallel/series-parallel systems |
| CN103279808A (en) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | Method for forecasting and analyzing circuit state transition probability of electric power system under severe weather |
| CN103887810A (en) * | 2014-03-12 | 2014-06-25 | 国家电网公司 | Method for judging whether continuous phase changing of direct-current system fails on basis of short-circuit ratio dynamic tracking |
| CN104376506A (en) * | 2014-11-17 | 2015-02-25 | 三峡大学 | Electric power system risk assessment method based on deterministic analysis |
| CN104934971A (en) * | 2015-06-09 | 2015-09-23 | 广东电网有限责任公司电力调度控制中心 | Dynamic section control method based on power flow transfer ratio |
-
2016
- 2016-12-06 CN CN201611106863.2A patent/CN106845757B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101615786A (en) * | 2009-07-24 | 2009-12-30 | 清华大学 | Grid system protection method based on local information of transformer substation |
| CN101895124A (en) * | 2010-05-27 | 2010-11-24 | 中国电力科学研究院 | Coupling assessment method for judging mutual influence of alternating current-direct current parallel/series-parallel systems |
| CN103279808A (en) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | Method for forecasting and analyzing circuit state transition probability of electric power system under severe weather |
| CN103887810A (en) * | 2014-03-12 | 2014-06-25 | 国家电网公司 | Method for judging whether continuous phase changing of direct-current system fails on basis of short-circuit ratio dynamic tracking |
| CN104376506A (en) * | 2014-11-17 | 2015-02-25 | 三峡大学 | Electric power system risk assessment method based on deterministic analysis |
| CN104934971A (en) * | 2015-06-09 | 2015-09-23 | 广东电网有限责任公司电力调度控制中心 | Dynamic section control method based on power flow transfer ratio |
Non-Patent Citations (1)
| Title |
|---|
| 徐慧明等: "基于潮流转移因子的广域后备保护方案", 《电网技术》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107316135A (en) * | 2017-06-16 | 2017-11-03 | 广西电网有限责任公司 | A kind of risk quantification appraisal procedure of switch-oriented and switch operation process |
| CN107316135B (en) * | 2017-06-16 | 2020-09-11 | 广西电网有限责任公司 | Risk quantitative evaluation method for switch and disconnecting link operation process |
| CN109687438A (en) * | 2018-12-19 | 2019-04-26 | 西南交通大学 | It is a kind of meter and high-speed rail impact load effect under power grid vulnerable line discrimination method |
| CN109687438B (en) * | 2018-12-19 | 2022-06-24 | 西南交通大学 | Power grid fragile line identification method considering high-speed rail impact load effect |
| CN110265998A (en) * | 2019-05-31 | 2019-09-20 | 国网辽宁省电力有限公司 | A kind of method and device that stabilized control system strategy definite value is checked |
| CN110265998B (en) * | 2019-05-31 | 2023-05-16 | 国网辽宁省电力有限公司 | Method and device for checking strategy fixed value of stability control system |
| CN111402075B (en) * | 2020-02-26 | 2023-04-07 | 国网浙江省电力有限公司经济技术研究院 | Power transmission network operation flexibility evaluation method considering line distribution factors |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106845757B (en) | 2020-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106845757A (en) | A kind of electric network swim shifts evaluating severity method | |
| CN109118098A (en) | The cascading failure methods of risk assessment and system of high proportion wind-electricity integration | |
| CN103219727B (en) | Partition power grid structure dynamics method of adjustment based on PMU actual measurement | |
| CN101800426A (en) | Safety level evaluation method of power grid | |
| CN106548265B (en) | Power transmission network reliability assessment method based on cascading failure accident chain search | |
| CN104505855B (en) | A kind of wind farm grid-connected characteristic online rating system | |
| CN106529791A (en) | Evaluation method for evaluating branch importance of power system | |
| CN102854415A (en) | Method for assessing lightning flashover risks of regional power grid lines | |
| CN109687438A (en) | It is a kind of meter and high-speed rail impact load effect under power grid vulnerable line discrimination method | |
| CN106408438A (en) | Method for automatically checking safety and stability control strategy based on dynamic simulation of whole process | |
| CN107453354B (en) | A kind of weak link recognition methods of power distribution network | |
| CN107368966A (en) | A kind of appraisal procedure of electric power networks fragility | |
| CN106655292A (en) | Ultra-high-voltage AC/DC hybrid grid stability analysis method | |
| CN107067127A (en) | A kind of power system excessive risk cascading failure fault chains searching method | |
| CN114336608A (en) | Unit blocking optimization method and system considering dynamic capacity increase and reconstruction | |
| CN107239896A (en) | A kind of security of relay protection appraisal procedure | |
| CN104156879A (en) | Electric power system protection and control method based on risk assessment | |
| Thomas et al. | Islanding detection using decision tree approach | |
| CN107067126A (en) | It is a kind of based on power flow transfer than thermally-stabilised key transmission channel recognition method | |
| CN105375453B (en) | A kind of emergency control method based on cascading failure destructiveness index | |
| CN111275312B (en) | Method and system for comprehensively evaluating and sequencing multiple fault severity of power grid | |
| CN111369388B (en) | Anti-vulnerability assessment method and system for urban power grid | |
| Tiwari et al. | Event identification and contingency assessment for voltage stability via PMU | |
| CN103218530B (en) | The reliability estimation method of a kind of DC transmission system stoppage in transit maintenance mode and system | |
| Nanab et al. | A quantitative method for assessing the resilience of infrastructure systems |
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 |