CN104281919A - Method for evaluating control performance of grid system - Google Patents

Method for evaluating control performance of grid system Download PDF

Info

Publication number
CN104281919A
CN104281919A CN201410504980.9A CN201410504980A CN104281919A CN 104281919 A CN104281919 A CN 104281919A CN 201410504980 A CN201410504980 A CN 201410504980A CN 104281919 A CN104281919 A CN 104281919A
Authority
CN
China
Prior art keywords
control
control zone
frequency
deviation
network system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410504980.9A
Other languages
Chinese (zh)
Other versions
CN104281919B (en
Inventor
李卫东
程凯
刘娆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201410504980.9A priority Critical patent/CN104281919B/en
Publication of CN104281919A publication Critical patent/CN104281919A/en
Application granted granted Critical
Publication of CN104281919B publication Critical patent/CN104281919B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (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)

Abstract

The invention discloses a method for evaluating the control performance of a grid system. The method for evaluating the control performance of the grid system comprises the following steps that the grid system is divided into multiple interconnected control areas, and the frequency deviation factors of all the control areas are obtained; whether the absolute value of the difference value betweenthe frequency deviation factors of any two control areas is larger than a preset threshold or not is judged, and all the control areas are evaluated according to a corrected CPS1 standard if the absolute value of the difference value between the frequency deviation factors of any two control areas is larger than the preset threshold, wherein the corrected CPS1 standard is obtained through an expression (please see the expression in the specifications); otherwise, all the control areas are evaluated according to a CPS1 standard, wherein the CPS1 standard is obtained through an expression (please see the expression in the specifications). According to the method for evaluating the control performance of the grid system, the problem that the situation is likely to happen that a frequency control target cannot be distributed fairly when control performance evaluation is carried out on control areas of different load scales can be solved well.

Description

A kind of network system control performance evaluation method
Technical field
The present invention relates to a kind of network system control performance evaluation method.
Background technology
In interconnected electric power system, the frequency control effect of each regional power grid is evaluated and specification by control performance assessment criteria, control performance assessment criteria demand fulfillment system operation reliability and the objectivity requirement evaluated regional power grid controlling behavior; At present, the CPS standard of what the control performance evaluation of each large regions electrical network of China adopted is North America electric reliability association (NERC), relative to the A standard previously adopted, CPS standard is the effect of behavior to interconnected electric power system in order to evaluate each control zone more objectively, guarantee that exchange power between control zone is within planned value scope, makes the frequency departure of interconnected electric power system within the limits prescribed; But, along with CPS standard economizes the enforcement of netting in part, find to there is particular law between CPS index and the scale of control area: the CPS index of larger control zone is generally better than the CPS index of the less control zone of scale, namely adopt CPS standard of the prior art to carry out easy frequency of occurrences control objectives ε when control performance is evaluated to interconnected electric power system 2can not fair allocat, and the adaptability problem of CPS standard evaluation scale great disparity Region control performance, as pointed out the regularity of distribution of CPS index between the region of different scales in existing document " CPS is for the applicability analysis of different load scale ", point out the average fluctuation due to district control deviation ACE and the non-linear proportional relation of frequency response coefficient simultaneously, cause the adaptability problem utilizing CPS standard evaluation scale great disparity Region control performance.
The formulation of CPS standard maintains system frequency by the district control deviation ACE of restriction interconnected electric power system to stablize, the specific requirement of described CPS standard in formula, Δ f is network system actual frequency f and scheduled frequency f 0deviation, be the one-minute average value of Δ f, ε be in given 1 year period, root-mean-square value; Suppose that an interconnected electric power system contains M control zone, there is identical frequency individual features each control zone, namely each control zone is under same service condition, its free-running frequency characteristic coefficient β is identical, and under the pattern that controls in order wire frequency deviation of interconnected network system, because the β of control zone is identical, therefore there is identical frequency bias coefficient B any moment control zone.If the root mean square of each control zone frequency departure when independent operating is the same, and the controlling behavior of control zone remains unchanged, so control zone interconnected after have the responsibility of equality to go to maintain common interconnected network frequency departure, according to district control deviation ACE and between relation, formula (1) is converted, draws:
RMS 2 { Δf ‾ } = AVG ( Δf ‾ 2 ) = AVG { [ Σ i ( ACE i ‾ k ) ] × Δf ‾ } = AVG { [ . . . + ( ACE i ‾ k ) + . . . ] × Δf ‾ } = . . . + AVG { ( ACE i ‾ k ) × Δf ‾ } + . . . , ACE in formula ibe the district control deviation of i-th control zone, district control deviation ACE ione-minute average value, k be the frequency response coefficient of system; In order to meet frequency control target, get k=-10MB, B is the frequency bias coefficient of region i, then because each control zone has identical frequency response characteristic, the responsibility of each control zone to maintenance interconnected network frequency is identical, as long as therefore meet (2) responsibility, so maintaining system frequency is fair to each control zone; In summary it can be seen that said process is based on all identical situation of the frequency bias coefficient of control zone in interconnected network system, if when contain the different control zone of frequency bias coefficient in interacted system, set a large control zone and comprise L control area, namely this large control zone is made up of L standard control area linear superposition, its load scale is L times of standard control zone, therefore frequency bias coefficient B c=LB, B cfor the frequency bias coefficient of large control zone, if the ACE now during each control zone independent operating meets the requirement of formula (2), L control zone requires it is the same to the control performance of frequency, therefore when they integrally run, as long as this large control zone meets AVG { ( ACE 1 + . . . + ACE L ‾ - 10 B ) × Δf ‾ ≤ L × ϵ 2 } - - - ( 3 ) , This large control zone is to safeguarding that the frequency stabilization responsibility of interconnected network and standard control zone are to safeguarding that the responsibility of interconnected network frequency stabilization is identical, the district control deviation ACE of this large control zone cbe the district control deviation sum of 1 to region, region L, substituted into formula (3) and obtain large control zone frequency responsibility born as a whole can find out thus, for the control zone that interconnected network load scale is different, as long as control this region to meet formula (4), it is fair that each control zone meets to interconnected network frequency the responsibility that frequency control target bears.
CPS1 standard focuses on the long-term assessment to system frequency quality, primary goal is the stability maintaining system frequency, namely require that the root-mean-square value of system frequency is less than system frequency control objectives value, because each control zone cannot be implemented to frequency requirement, by the relation of frequency and district control deviation, analysis conversion is carried out to formula (1), and then CPS1 index is formed to the ACE restriction of control zone; Lower mask body is divided into load scale identical and different two kinds of situations in each control zone to analyze:
When each control zone load scale is all identical, in order to the physical significance of every variable in formula is better described, first by formula the control zone number M of interconnected network is multiplied by both sides simultaneously, obtains m ε in formula (5) 2be defined as frequency adjustment responsibility, primary and foremost purpose frequency adjustment responsibility being assigned to M control zone is the allocation scheme confirming responsibility, frequency adjustment responsibility has two kinds of allocation scheme: distribute according to each control zone load scale, and distributes according to each control zone load fluctuation size; In employing load scale allocation scheme situation, the load scale of control zone is equal to the size of control zone to frequency regulation capability, and therefore load scale in control zone is larger, and the frequency regulating duty born is also larger; Under this allocation model, the frequency bias coefficient linear proportional example relation of the frequency regulating duty that control zone is born in interconnected network and this control zone.Under employing control zone load fluctuation size distribution condition, load fluctuation is larger, and the difficulty of control zone frequency adjustment is larger, to get frequency control responsibility larger; Under this allocation model, the frequency regulating duty that control zone is born in interconnected network should distribute according to the long-time statistical characteristic of this control zone load fluctuation.When carrying out assigned frequency control responsibility to formula (5) by load scale, because the load scale of each control zone is all identical, therefore the frequency control responsibility born of each control zone is also identical, therefore the frequency adjustment responsibility that each control zone is distributed is ε 2; When assigned frequency adjustment responsibility is come according to load fluctuation size to formula (5), because the service condition of control zone is all similar, can think the long-time statistical characteristic close of each control zone load fluctuation, the frequency adjustment responsibility that therefore control zone is born under this allocation model is all ε 2; Show that, when interconnected electric power system is when the load scale of each control zone is all identical, carry out assigned frequency regulating and controlling responsibility from load scale and load fluctuation two kinds of angles, formula (4) all can reach the requirement of fair allocat thus.
When in interconnected network containing load scale not identical control zone time, when assigned frequency adjustment responsibility is come according to load scale to formula (5), because the load scale of this large control zone is L times of standard control zone, therefore the frequency responsibility that this control zone is distributed is L ε 2, again due to the frequency bias coefficient B of this control zone c=LB, so under this allocation model, think that each control zone of interconnected network meets formula (4) and has general adaptability to meeting system to frequency quality requirement and control zone responsibility distributional equity, to formula (5) when carrying out assigned frequency according to load fluctuation size and regulating responsibility, along with the increase of load scale, load fluctuation is not increase according to the proportional linearity of load scale, because there is the smoothing effect of load fluctuation in interconnected network, namely along with load scale in control zone increases, because in region, the load fluctuation of diverse location is cancelled out each other, load overall undulatory property in control zone is weakened, such load scale is larger, load smoothing effect is further remarkable, and now the unit capacity of control zone is that the proportional linearity of same load scale increases, the control zone that therefore load scale is larger, it is more prone to the control of frequency, and formula (4) is not suitable for according to load fluctuation allocation model.
Summary of the invention
The present invention is directed to the proposition of above problem, and develop a kind of network system control performance evaluation method.
Technological means of the present invention is as follows:
A kind of network system control performance evaluation method, comprises the steps:
Step 1: described network system is divided into several intarconnected cotrol districts, draws the frequency bias coefficient of each control zone;
Step 2: judge whether the absolute difference of the frequency bias coefficient of any two control zones is greater than predetermined threshold value, is perform step 3, otherwise perform step 4;
Step 3: evaluate each control zone respectively according to revised CPS1 standard, described revised CPS1 standard is AVG { ( VCE i ‾ - 10 B s ) × Δf ‾ } ≤ - 10 B i Σ i = 1 n - 10 B i × ϵ 2 , In formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
Step 4: evaluate each control zone respectively according to CPS1 standard, described CPS1 standard is in formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
Further, described step 3 comprises the steps:
1. Real-time Collection network system actual frequency f, obtains network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
2. the Tie line Power P of Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
3. according to formula calculate in the revised CPS1 standard of control zone i and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B sfor the frequency bias coefficient that network system is total, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
What 4. 3. calculate when step meets factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What 3. calculate when step meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What 3. calculate when step meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor the frequency bias coefficient of control zone i, B sfor the frequency bias coefficient that network system is total, ε 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period;
Further, described step 4 comprises the steps:
I, Real-time Collection network system actual frequency f, obtain network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
The Tie line Power P of II, Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
III, according to formula calculate in the iCPS1 standard of control zone and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
IV, calculate when step III meet factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What calculate when step III meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What calculate when step III meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor frequency bias coefficient, the B of control zone i sfor frequency bias coefficient, ε that network system is total 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period;
Further, described revised CPS1 standard draws as follows:
A, detect corresponding time point t kcontrol zone i load fluctuation value Δ L i(t k), described time point forms a set T n={ t k, k=1,2 ..., K, have time interval W between adjacent time point, the setting of W makes K meet central limit theorem;
B, calculated load undulating quantity Δ L i(t k) expectation value and standard deviation ;
C, structure control zone i load fluctuation model t in formula kfor measuring the time point of load fluctuation, Δ L i(t k) equal t for time point ktime load fluctuation value;
D, calculated load volatility model expectation value and standard deviation wherein &sigma; &Delta;L &OverBar; n = 1 K &sigma; &Delta; L I , Draw load fluctuation model submit to normal distribution;
E, calculate the expected value and standard deviation of the large control zone be made up of M control zone, the expectation value of described large control zone is the standard deviation of described large control zone is draw the increase along with load scale, the expectation value of load fluctuation increases according to increasing the consistent multiple of multiple with load scale, and the standard deviation of load fluctuation increases according to increasing the consistent multiple of multiple with the load scale after evolution;
F, by the frequency control target ε of interconnected electric power system 2increase multiple according to the load scale after evolution to distribute, draw revised CPS1 standard AVG { ( VCE i &OverBar; - 10 B s ) &times; &Delta;f &OverBar; } &le; - 10 B i &Sigma; i = 1 n - 10 B i &times; &epsiv; 2 .
Owing to have employed technique scheme, a kind of network system control performance evaluation method provided by the invention, predetermined threshold value whether is greater than as decision condition using the absolute difference of the frequency bias coefficient of any two control zones, and then each control zone with different load scale otherness can be distinguished, thus adopt CPS1 standard or revised CPS1 standard to evaluate to each control zone correspondence with different load scale otherness, solve the adaptability problem adopting CPS standard evaluation scale of the prior art great disparity Region control performance, the CPS1 index of control zone load scale being differed greatly by the method for the invention is similar in the allocation proportion in each interval.Formula medium frequency after adjustment controls responsibility allocation scheme and considers the impact of load smoothing effect on the long-time statistical characteristic of load fluctuation, easy frequency of occurrences control objectives to be carried out when control performance is evaluated in the control zone that can solve preferably between different load scale can not the problem of fair allocat, the revised CPS1 standard of the present invention considers the smoothing effect of load fluctuation to the impact of the fluctuation mean value of district control deviation, the Region control performance of different load scale can be evaluated more objectively, and the frequency control target of interconnected electric power system can be met.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of the method for the invention.
Embodiment
A kind of network system control performance evaluation method as shown in Figure 1, comprises the steps:
A kind of network system control performance evaluation method, comprises the steps:
Step 1: described network system is divided into several intarconnected cotrol districts, draws the frequency bias coefficient of each control zone;
Step 2: judge whether the absolute difference of the frequency bias coefficient of any two control zones is greater than predetermined threshold value, is perform step 3, otherwise perform step 4;
Step 3: evaluate each control zone respectively according to revised CPS1 standard, described revised CPS1 standard is AVG { ( VCE i &OverBar; - 10 B s ) &times; &Delta;f &OverBar; } &le; - 10 B i &Sigma; i = 1 n - 10 B i &times; &epsiv; 2 , In formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n, n be the number of the control zone that network system marks off;
Step 4: evaluate each control zone respectively according to CPS1 standard, described CPS1 standard is in formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
Further, described step 3 comprises the steps:
1. Real-time Collection network system actual frequency f, obtains network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
2. the Tie line Power P of Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
3. according to formula calculate in the revised CPS1 standard of control zone i and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B sfor the frequency bias coefficient that network system is total, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
What 4. 3. calculate when step meets factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What 3. calculate when step meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What 3. calculate when step meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor the frequency bias coefficient of control zone i, B sfor the frequency bias coefficient that network system is total, ε 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period;
Further, described step 4 comprises the steps:
I, Real-time Collection network system actual frequency f, obtain network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
The Tie line Power P of II, Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
III, according to formula calculate in the iCPS1 standard of control zone and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
IV, calculate when step III meet factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What calculate when step III meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What calculate when step III meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor frequency bias coefficient, the B of control zone i sfor frequency bias coefficient, ε that network system is total 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period;
Further, described revised CPS1 standard draws as follows:
A, detect corresponding time point t kcontrol zone i load fluctuation value Δ L i(t k), described time point forms a set T n={ t k, k=1,2 ..., K, have time interval W between adjacent time point, the setting of W makes K meet central limit theorem;
B, calculated load undulating quantity Δ L i(t k) expectation value and standard deviation ; Be be converted to another steady state (SS) from current steady state (SS) under network system load disturbance at any time, load fluctuation be not independent of each other in the same time as stochastic variable, set { the Δ L of load fluctuation i(t k): k=1,2 ..., K} is one group of independently random distribution data;
C, structure control zone i load fluctuation model t in formula kfor measuring the time point of load fluctuation, Δ L i(t k) for time point be t ktime load fluctuation value;
D, calculated load volatility model expectation value and standard deviation wherein &sigma; &Delta;L &OverBar; n = 1 K &sigma; &Delta; L I , Draw load fluctuation model submit to normal distribution;
E, calculate the expected value and standard deviation of the large control zone be made up of M control zone, the expectation value of described large control zone is the standard deviation of described large control zone is draw the increase along with load scale, the expectation value of load fluctuation increases according to increasing the consistent multiple of multiple with load scale, and the standard deviation of load fluctuation increases according to increasing the consistent multiple of multiple with the load scale after evolution; A large control zone is made up of several control zones linear superposition, and then the load fluctuation of this large control zone is also be made up of the load linear superposition of each control zone, and the load fluctuation of each control zone is separate, according to Probability Theory and Math Statistics relevant knowledge, the load fluctuation of this large control zone is also Normal Distribution;
F, by the frequency control target ε of interconnected electric power system 2increase multiple according to the load scale after evolution to distribute, draw revised CPS1 standard AVG { ( VCE i &OverBar; - 10 B s ) &times; &Delta;f &OverBar; } &le; - 10 B i &Sigma; i = 1 n - 10 B i &times; &epsiv; 2 .
The present invention adopts the standard deviation of load fluctuation to characterize load fluctuation characteristic, standard deviation has better intuitive, and proximity test value itself more, the fluctuation size of a sample can be weighed, the dispersion degree of reflected sample exactly, and then characterize load fluctuation characteristic fully and effectively.
The present invention gives the allocation scheme of two kinds of frequency adjustment responsibilities according to the difference of each control zone load scale otherness, lower mask body passes judgment on the responsibility distributional equity sex chromosome mosaicism of each control zone under two kinds of allocation scheme, namely under these two kinds of allocation models, whether the Region control effect of each control zone is the same, the index of evaluation region control effects is the statistical property from seeing the frequency departure of this control zone within longer time period, represents with the root-mean-square value of frequency departure; Under the pattern of distributing according to load scale, due to be with load scale for distribute foundation, Region control performance has nothing to do with the load fluctuation size of control zone, therefore when each control zone independent operating, keep and interconnected operation time consistent, to the requirement of the root-mean-square value of frequency departure all the same (frequency control target is the same) during each control zone independent operating, be then fair when control zone is interconnected to the distribution of frequency adjustment responsibility; Under the pattern of distributing according to load fluctuation, due to be with load fluctuation be distribute foundation, therefore Region control performance is relevant with load fluctuation, the service condition when service condition when each control area is in isolated operation and their interconnected operations is consistent, and the root-mean-square value of the frequency departure of each control zone independent operating time all the same, each control zone just arrives fair principle to the distribution of frequency adjustment responsibility.
When control zone c independent operating (each control zone interconnected after frequency departure be consistent), according to control zone control deviation ACE ccomputing formula, can obtain the frequency deviation f during c independent operating of control zone cone-minute average value, therefore can change into being the one-minute average value of network system frequency departure, is the coefficient result in each control zone, frequency departure during reflection control zone c independent operating and the relation between the frequency departure of interconnected network; And then AVG { ( ACE c - 10 B c ) &OverBar; &times; &Delta;f &OverBar; } &le; &epsiv; 2 Be converted into AVG { &Delta;f &OverBar; c &times; &Delta;f &OverBar; } &epsiv; 2 - - - ( 6 ) , From formula (6), after the frequency adjustment responsibility distribution of interconnected network, frequency departure when limiting control zone independent operating and the relation between the frequency departure of system, frequency departure due to system is identical, therefore after distributing according to load scale pattern, each control zone is when service condition is constant, and Region control effect during independent operating all needs to be not more than identical frequency control target ε 0, namely under this allocation model, meet the principle of frequency fair allocat, therefore CPS1 standard is applicable to load scale as distributing the pattern of foundation, according under load fluctuation allocation model, according to CPS1 standard, each control zone all demand fulfillment formula (7) during independent operating after frequency adjustment responsibility is distributed, because the existence of load fluctuation smoothing effect in interconnected network, make the control zone that load scale is larger, load fluctuation is relatively more level and smooth, AGC (control of automatic generation amount) performance is better, and then frequency departure during the independent operating of control zone is less, if or according to load scale, frequency control general objective is distributed, then each control zone demand fulfillment formula (7) requirement, but the control zone that scale is larger, the one-minute average value of the frequency departure of independent operating is less, control zone also less, therefore the requirement of formula (7) more easily meets, and this is obviously unfair to the control zone of small scale.
Existing CPS1 standard does not consider the smoothing effect of load, thinks the long-time statistical characteristic of load fluctuation and the frequency bias coefficient linear proportional relation of this control zone, thus distributes the overhead control target ε of interconnected electric power system according to frequency bias coefficient size 2, a kind of network system control performance evaluation method provided by the invention, predetermined threshold value whether is greater than as decision condition using the absolute difference of the frequency bias coefficient of any two control zones, and then each control zone with different load scale otherness can be distinguished, thus adopt CPS1 standard or revised CPS1 standard to evaluate to each control zone correspondence with different load scale otherness, solve the adaptability problem adopting CPS standard evaluation scale of the prior art great disparity Region control performance, the CPS1 index of control zone load scale being differed greatly by the method for the invention is similar in the allocation proportion in each interval.Formula medium frequency after adjustment controls responsibility allocation scheme and considers the impact of load smoothing effect on the long-time statistical characteristic of load fluctuation, easy frequency of occurrences control objectives to be carried out when control performance is evaluated in the control zone that can solve preferably between different load scale can not the problem of fair allocat, the revised CPS1 standard of the present invention considers the smoothing effect of load fluctuation to the impact of the fluctuation mean value of district control deviation, the Region control performance of different load scale can be evaluated more objectively, and the frequency control target of interconnected electric power system can be met.
MATLAB/Simulink is utilized to carry out simulation analysis to the network system being divided into two intarconnected cotrol districts below, two intarconnected cotrol districts are control zone A and B, wherein the frequency bias coefficient of control zone A is 42.66, the frequency bias coefficient of control zone B is 21.33, and namely the frequency bias coefficient of control zone A is 2 times of control zone B; The load fluctuation model of control zone B submits to normal distribution, and the load fluctuation model of control zone A is made up of the load fluctuation model linear superposition of 2 control zone B.Table 1 is each interval statistics value (the evaluation period is 1min) of CPS1 standard, as known from Table 1, the CPS1 index entirety of control zone A is better than the index of control zone B, and the ratio being less than 100% interval that is distributed in of control zone A is starkly lower than the ratio that control zone B distributes in this interval; And the ratio that control zone A is distributed in 100% ~ 200% interval is significantly higher than the ratio that control zone B distributes in this interval, control zone A is being greater than the ratio in 200% interval a little less than the ratio of control zone B in this interval.
Table 1. liang intarconnected cotrol district CPS1 index distribution.
Between CPS1 Index areas <100% 100%~200% >200%
Control zone A 0.0945 0.7120 0.1934
Control zone B 0.1689 0.6121 0.2189
Table 2 carries out revised simulation result to CPS1 standard, now the ratio that CPS1 index is being less than 100% interval of control zone A is increased to 11.34% by 9.45%, distribution in 100% ~ 200% interval drops to 69.32% by 71.2%, does not become in the distribution being greater than 200% interval; Be distributed in the ratio being less than 100% interval in the B index set of control zone and drop to 11.26% by 16.89%, be increased to 66.84% at the interval distribution proportion of 100% ~ 200% by 61.21%, do not change in the interval being greater than 200%; After A and B each interval Distribution Indexes ratio in control zone changes, control zone A is substantially identical with the distribution proportion that control zone B is being less than 100% interval, and the distribution proportion in 100% ~ 200% interval is also close; Therefore revised CPS1 standard can solve the malapportioned problem in control zone between different load scale preferably.
The revised CPS1 index distribution in table 2. liang intarconnected cotrol district.
Between CPS1 Index areas <100% 100%~200% >200%
Control zone A 0.1134 0.6932 0.1934
Control zone B 0.1126 0.6684 0.2189
The above; be only the present invention's preferably embodiment; but protection scope of the present invention is not limited thereto; anyly be familiar with those skilled in the art in the technical scope that the present invention discloses; be equal to according to technical scheme of the present invention and inventive concept thereof and replace or change, all should be encompassed within protection scope of the present invention.

Claims (4)

1. a network system control performance evaluation method, is characterized in that comprising the steps:
Step 1: described network system is divided into several intarconnected cotrol districts, draws the frequency bias coefficient of each control zone;
Step 2: judge whether the absolute difference of the frequency bias coefficient of any two control zones is greater than predetermined threshold value, is perform step 3, otherwise perform step 4;
Step 3: evaluate each control zone respectively according to revised CPS1 standard, described revised CPS1 standard is AVG { ( VCE i &OverBar; - 10 B s ) &times; &Delta;f &OverBar; } &le; - 10 B i &Sigma; i = 1 n - 10 B i &times; &epsiv; 2 , In formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
Step 4: evaluate each control zone respectively according to CPS1 standard, described CPS1 standard is in formula for the control deviation ACE of control zone i ione-minute average value, B sfor the frequency bias coefficient that network system is total, for the one-minute average value of each control zone frequency deviation f, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n.
2. a kind of power grid control method of evaluating performance according to claim 1, is characterized in that described step 3 comprises the steps:
1. Real-time Collection network system actual frequency f, obtains network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
2. the Tie line Power P of Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
3. according to formula calculate in the revised CPS1 standard of control zone i and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B sfor the frequency bias coefficient that network system is total, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
What 4. 3. calculate when step meets factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What 3. calculate when step meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What 3. calculate when step meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor the frequency bias coefficient of control zone i, B sfor the frequency bias coefficient that network system is total, ε 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period.
3. a kind of power grid control method of evaluating performance according to claim 1, is characterized in that described step 4 comprises the steps:
I, Real-time Collection network system actual frequency f, obtain network system actual frequency f and scheduled frequency f 0difference DELTA f=f-f 0, then the frequency departure of each control zone is Δ f, and calculates the one-minute average value of each control zone frequency deviation f wherein f is interconnected network system frequency, the f of Real-time Collection 0for network system scheduled frequency;
The Tie line Power P of II, Real-time Collection connection control district i t, obtain the Tie line Power P of control zone i twith scheduled net interchange P 0difference DELTA P t=P t-P 0, according to ACE i=-10B iΔ f+ Δ P tobtain the control deviation ACE of control zone i i, wherein ACE ifor the control deviation of control zone i, Δ f is the frequency departure of control zone i, Δ P tfor Tie line Power P twith scheduled net interchange P 0difference, B ifor the frequency bias coefficient of control zone i, i=1,2 ..., n;
III, according to formula calculate in the iCPS1 standard of control zone and meet factor CF 1, wherein for the control deviation ACE of control zone i ione-minute average value, for the one-minute average value of control zone i frequency deviation f, B ifor the frequency bias coefficient of control zone i, ε is the root-mean-square value of the one-minute average value of Δ f in given 1 year period, i=1,2 ..., n;
IV, calculate when step III meet factor CF 1when≤0, automatic generation amount control mode is the control deviation ACE of retentive control district i iconstant;
What calculate when step III meets the factor 0 < CF 1when≤1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ithe direction change of trend zero;
What calculate when step III meets factor CF 1during > 1, automatic generation amount control mode is the active power adjusting corresponding automatic-generation-control unit, makes the control deviation ACE of control zone i ibe less than wherein b ifor frequency bias coefficient, the B of control zone i sfor frequency bias coefficient, ε that network system is total 10for the root-mean-square value of ten minutes mean value of network system actual frequency and scheduled frequency deviation in given 1 year period.
4. a kind of power grid control method of evaluating performance according to claim 1, is characterized in that described revised CPS1 standard draws as follows:
A, detect corresponding time point t kcontrol zone i load fluctuation value Δ L i(t k), described time point forms a set T n={ t k, k=1,2 ..., K, have time interval W between adjacent time point, the setting of W makes K meet central limit theorem;
B, calculated load undulating quantity Δ L i(t k) expectation value and standard deviation ;
C, structure control zone i load fluctuation model t in formula kfor measuring the time point of load fluctuation, Δ L i(t k) equal t for time point ktime load fluctuation value;
D, calculated load volatility model expectation value and standard deviation wherein &sigma; &Delta;L &OverBar; n = 1 K &sigma; &Delta; L I , Draw load fluctuation model submit to normal distribution;
E, calculate the expected value and standard deviation of the large control zone be made up of M control zone, the expectation value of described large control zone is the standard deviation of described large control zone is draw the increase along with load scale, the expectation value of load fluctuation increases according to increasing the consistent multiple of multiple with load scale, and the standard deviation of load fluctuation increases according to increasing the consistent multiple of multiple with the load scale after evolution;
F, by the frequency control target ε of interconnected electric power system 2increase multiple according to the load scale after evolution to distribute, draw revised CPS1 standard AVG { ( VCE i &OverBar; - 10 B s ) &times; &Delta;f &OverBar; } &le; - 10 B i &Sigma; i = 1 n - 10 B i &times; &epsiv; 2 .
CN201410504980.9A 2014-09-26 2014-09-26 A kind of network system control performance evaluation method Expired - Fee Related CN104281919B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410504980.9A CN104281919B (en) 2014-09-26 2014-09-26 A kind of network system control performance evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410504980.9A CN104281919B (en) 2014-09-26 2014-09-26 A kind of network system control performance evaluation method

Publications (2)

Publication Number Publication Date
CN104281919A true CN104281919A (en) 2015-01-14
CN104281919B CN104281919B (en) 2017-10-24

Family

ID=52256777

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410504980.9A Expired - Fee Related CN104281919B (en) 2014-09-26 2014-09-26 A kind of network system control performance evaluation method

Country Status (1)

Country Link
CN (1) CN104281919B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104898644A (en) * 2015-04-29 2015-09-09 中国电力科学研究院 Evaluation method of frequency control performance of control areas in interconnected power grid
CN105808889A (en) * 2016-04-25 2016-07-27 中国电力科学研究院 Frequency deviation factor simulation collocation method
CN105958524A (en) * 2016-07-11 2016-09-21 大连理工大学 DC frequency modulation participated area control performance evaluation method
CN106936132A (en) * 2017-04-17 2017-07-07 中国石油大学(华东) A kind of multi-area Interconnected Power System STOCHASTIC CONTROL device for evaluating performance and method based on entropy
CN109301845A (en) * 2018-12-11 2019-02-01 国网江苏省电力有限公司扬州供电分公司 Method is stabilized in the active fluctuation of extra-high voltage interconnection based on master-slave mode energy storage coordinated control
CN112711237A (en) * 2020-12-29 2021-04-27 华润电力技术研究院有限公司 Automatic control quality online evaluation method and system for thermal power generating unit

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010193605A (en) * 2009-02-18 2010-09-02 Tokyo Electric Power Co Inc:The Load estimating method of power distribution section and power distribution system control method
CN102324885A (en) * 2011-09-20 2012-01-18 大连理工大学 Method and system for controlling CPS (Control Performance Standard) oriented automatic generation control unit
CN102497147A (en) * 2011-11-25 2012-06-13 辽宁省电力有限公司 Control method of AGC facing control performance standard (CPS) and control system thereof
CN103166221A (en) * 2013-03-30 2013-06-19 大连理工大学 Ultra-high voltage grid operating control mode for controlling in-series execution of performance evaluation criterion
CN103744291A (en) * 2014-01-13 2014-04-23 武汉大学 AGC (Automatic Generation Control) logic control system and method based on CPS (Control Performance Standard)

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010193605A (en) * 2009-02-18 2010-09-02 Tokyo Electric Power Co Inc:The Load estimating method of power distribution section and power distribution system control method
CN102324885A (en) * 2011-09-20 2012-01-18 大连理工大学 Method and system for controlling CPS (Control Performance Standard) oriented automatic generation control unit
CN102497147A (en) * 2011-11-25 2012-06-13 辽宁省电力有限公司 Control method of AGC facing control performance standard (CPS) and control system thereof
CN103166221A (en) * 2013-03-30 2013-06-19 大连理工大学 Ultra-high voltage grid operating control mode for controlling in-series execution of performance evaluation criterion
CN103744291A (en) * 2014-01-13 2014-04-23 武汉大学 AGC (Automatic Generation Control) logic control system and method based on CPS (Control Performance Standard)

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
熊卫斌 等: "CPS对于不同负荷规模的适用性分析", 《电力系统自动化》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104898644A (en) * 2015-04-29 2015-09-09 中国电力科学研究院 Evaluation method of frequency control performance of control areas in interconnected power grid
CN104898644B (en) * 2015-04-29 2018-01-19 中国电力科学研究院 The FREQUENCY CONTROL method of evaluating performance of each control zone in a kind of interconnected network
CN105808889A (en) * 2016-04-25 2016-07-27 中国电力科学研究院 Frequency deviation factor simulation collocation method
CN105808889B (en) * 2016-04-25 2023-05-23 中国电力科学研究院 Frequency deviation coefficient simulation configuration method
CN105958524A (en) * 2016-07-11 2016-09-21 大连理工大学 DC frequency modulation participated area control performance evaluation method
CN106936132A (en) * 2017-04-17 2017-07-07 中国石油大学(华东) A kind of multi-area Interconnected Power System STOCHASTIC CONTROL device for evaluating performance and method based on entropy
CN106936132B (en) * 2017-04-17 2019-08-02 中国石油大学(华东) A kind of multi-area Interconnected Power System STOCHASTIC CONTROL device for evaluating performance and method based on entropy
CN109301845A (en) * 2018-12-11 2019-02-01 国网江苏省电力有限公司扬州供电分公司 Method is stabilized in the active fluctuation of extra-high voltage interconnection based on master-slave mode energy storage coordinated control
CN109301845B (en) * 2018-12-11 2020-11-03 国网江苏省电力有限公司扬州供电分公司 Active fluctuation stabilizing method of extra-high voltage tie line based on master-slave type energy storage coordination control
CN112711237A (en) * 2020-12-29 2021-04-27 华润电力技术研究院有限公司 Automatic control quality online evaluation method and system for thermal power generating unit

Also Published As

Publication number Publication date
CN104281919B (en) 2017-10-24

Similar Documents

Publication Publication Date Title
CN104281919A (en) Method for evaluating control performance of grid system
CN103559655B (en) The Forecasting Methodology of the novel feeder line load of microgrid based on data mining
Banerjee et al. Reliability based optimum location of distributed generation
US9906030B2 (en) Optimal control method for reactive voltage of wind power and photovoltaic power centralized grid connection
US20160253768A1 (en) Method for determining integrated network loss rate in uhv ac cross-regional electricity trading
CN106066423A (en) A kind of analysis method of opposing electricity-stealing based on Loss allocation suspicion analysis
Péan et al. Environmental and economic impact of demand response strategies for energy flexible buildings
KR20140105506A (en) Adaptation of a power generation capacity and determining of an energy storage unit size
CN103544537B (en) Based on the cumulative short-term load forecasting method of comprehensive subnet of reliability forecasting assessment
US20150194812A1 (en) Metering optimal sampling
CN110011358B (en) Distribution network load state adjustment controller
CN105225021A (en) The optimum choice method of power distribution network project yet to be built
Wang et al. Load aggregation methods for quasi-static power flow analysis on high PV penetration feeders
Triebs et al. From heat demand to heat supply: How to obtain more accurate feed-in time series for district heating systems
CN116644920A (en) Intelligent power failure management method and system
Konishi et al. Optimal allocation of photovoltaic systems and energy storage systems considering constraints of both transmission and distribution systems
Dong et al. A research of real-time pricing mechanism and its characteristics
CN105322543A (en) Incidence matrix based load distribution system and method
CN108334990A (en) A kind of bulk power grid reactive-load compensation addressing and capacity optimization method and system
CN114899942A (en) Monitoring system of active power distribution network
Yang et al. The regulating characteristic analysis for distributed electric heating load in Northern China
CN113657032A (en) Low-frequency load shedding method and system for pre-centralized coordination and real-time distributed control
US20190081478A1 (en) Method for monitoring and controlling an electrical network
Liu et al. Aggregation model of air conditioning load considering temperature sensor accuracy
Zhang et al. A robust interval economic dispatch model accommodating large-scale wind power generation with consideration of price-based demand response

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171024

Termination date: 20200926