CN105977968B - Optimal configuration method for power quality monitor of annular multi-source power distribution network - Google Patents
Optimal configuration method for power quality monitor of annular multi-source power distribution network Download PDFInfo
- Publication number
- CN105977968B CN105977968B CN201610393772.5A CN201610393772A CN105977968B CN 105977968 B CN105977968 B CN 105977968B CN 201610393772 A CN201610393772 A CN 201610393772A CN 105977968 B CN105977968 B CN 105977968B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- endpoint
- matrix
- busbar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 98
- 238000012544 monitoring process Methods 0.000 claims abstract description 70
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 238000012512 characterization method Methods 0.000 claims description 8
- 238000007596 consolidation process Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 230000000875 corresponding effect Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 238000000265 homogenisation Methods 0.000 claims description 3
- 230000006855 networking Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000010937 topological data analysis Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000012806 monitoring device Methods 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 5
- 230000002457 bidirectional effect Effects 0.000 abstract description 3
- 230000006698 induction Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 3
- 238000013139 quantization Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Power Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The method is suitable for new concept definition of topology analysis of the annular multi-source power grid structure, and comprises a positive direction, a bidirectional equivalent line and an endpoint; constructing an annular multisource power grid structure representation matrix, including correcting an adjacent matrix, a system structure matrix, a tide direction matrix and a monitoring observable matrix; defining weight factors, including general rule induction of power quality monitor arrangement, attribute classification of different endpoints and definition of corresponding weight coefficients; ambiguity index definition, including the system overall ambiguity, partial maximum ambiguity and monitoring system ambiguity; establishing an optimal configuration model of the power quality monitor; and solving the optimization model to determine the optimal configuration scheme of the system power quality monitor.
Description
Technical field
The present invention relates to a kind of power quality monitor Optimal Configuration Method of ring-type multiple source power distribution net, belong to electrical engineering and
Power quality field.
Background technology
In following intelligent distribution network, structure can realize the network of continuous multiple spot, on-line measurement and information fusion at the same time
Change electric energy quality monitoring system (networked power quality monitoring system, NPQMS) becomes development
Trend.Configure yet with fund, the limitation of technology and on a large scale power quality monitor (power quality
Monitor, PQM) data redundancy that brings and the problems such as heavy communication task, how quantity to PQM in NPQMS and installation position
Put and make optimal selection, that is, realize that distributing rationally for PQM is of great significance.
At present, it is less to distribute relevant achievement in research rationally with PQM.A small amount of document is proposed based on Graph Theory Analysis Method, section
The method that point analysis method etc. realizes PQM preferred arrangements, but the radial distribution networks network of single supply is generally only applicable to, and only consider
The ornamental of power quality index, and fail to consider power quality disturbance (power quality disturbance
Source, PQDS) positioning etc. advanced intelligent diagnostics demand.Actually, on the one hand, to improve the security reliability of power supply, more
Carry out more important area loads and use ring connection form, while as the high permeability of more and more distributed generation resources connects
Enter, power distribution network will generally become the complex topology structure of multi-source;On the other hand, NPQMS will be wanted more in following intelligent grid
The power quality high-grade intelligent diagnostic functions such as the realistic automatic accurate positioning of existing PQDS.In terms of patent achievement, although and power quality
Relevant achievement is very much, but the overwhelming majority concentrates on power quality device, and such as 2015107407598,2014105143834;Electricity
Can mass monitoring system, such as 2013100476831,2015109143522;Power Quality Disturbance is analyzed, such as
2015108224005,2014107853645;Electricity quality evaluation, such as 2011102114770,2014106381100;Electric energy matter
Measure disturbed depth, such as 2007100517107,2009102167982;Utility power quality control, such as 2011103534079,
2012105091794;Etc..And with PQM distribute rationally it is relevant seldom, the invention of only Application No. 201610017639X is special
Profit proposes a kind of distribution network electric energy quality monitoring point configuration side based on Modified particle swarm optimization algorithm, meter and distributed generation resource
Method, but it carries the distribution network that algorithm is also only applicable to the radial topological structure of single supply.Patent of the present invention optimizes for PQM
Allocation problem is studied, and proposes that one kind can be suitable for ring-type multi-source topological structure distribution network systems, and takes into account PQDS positioning etc.
The PQM Optimal Configuration Methods of high-grade intelligent diagnostic function demand.
The content of the invention
The problem of present invention will overcome existing PQM Optimal Configuration Methods to be effectively applicable to contain more power supply power network with circle structure,
Based on the definition of positive direction, two-way equivalent circuit and endpoint, and correct the structure realization system of adjacency matrix, direction of tide matrix
The Efficient Characterization for structural information of uniting, defines the importance degree of different each attribute endpoints of weight factor quantization signifying, establishes system
Monitoring point optimization allocation models and by monitoring the system ambiguous degree of considerable matrix check, filters out the PQM of ring-type multi-source network most
Excellent allocation plan.
The present invention is to achieve the above object, it is proposed that a kind of PQM optimizations that can be suitable for ring-type multiple source power distribution network are matched somebody with somebody
Method is put, as shown in Figure 1, its process includes the following steps:
1st, new concept of the definition suitable for ring-type multi-source electric network composition topological analysis:With radial electric network composition feature not
Together, ring-type multi-source power system operating mode and state have changeability, and part circuit current direction is there are a variety of possibility, so as to influence
The monitoring range of each node configuration PQM, it is desirable to account for various possible situations;
Step 101, the regulation of positive direction:For matching PQDS direction determining methods, it is necessary to establish a fixed coordinate system, if
It is positive direction to determine to be directed toward line direction by busbar;
Step 102, the definition of two-way equivalent circuit:By every line equivalent in addition to power supply, load lead-out wire into two
The circuit with fixed current flow direction of parallel bidirectional;
Step 103, the definition of endpoint:Line current inflow end is arranged to the endpoint of configurable PQM to represent the circuit,
I.e. every circuit can equivalently represented be two endpoints;
2nd, the structure of ring-type multi-source system structural characterization matrix:Analyzed based on step 1 ring-network topology, structure is corrected adjacent respectively
Matrix, system structure matrix, direction of tide matrix, the considerable matrix description network structure of monitoring are connect, offer number is distributed rationally for PQM
According to support;
Step 201, the structure of adjacency matrix is corrected:Adjacency matrix AorgIt is usually used in Power Network Topology Analysis, its form isi1,i2=1,2 ..., NBNumbered for bus nodes, NBFor system busbar number of nodes;Matrix element
Assignment principle such as formula (1) shown in;Because load and power supply are all important monitorings pair in networking electric energy quality monitoring system
As, it is necessary to suitably be corrected to adjacency matrix:Represent to load with busbar and the connecting line being equivalent between the earth of busbar 0, use
Busbar represents power supply with the connecting line of its own;Correct adjacency matrixForm be For
Add the system node numbering after the busbar 0 for representing the earth;ElementRepresent the structure of each endpoint in ring-type multi-source network
Numbering, its value represent that the endpoint whether there is, shown in assignment principle such as formula (2);
Step 202, system structure matrix is built:WillAll valid endpoints in matrix, i.e.,Endpoint, press
Row order extracts, and it is sequentially labeled as to new system end point symbol t successivelyj, wherein j=1,2 ... T is successively sequentially
T after markjNumbering, T is valid endpoint total quantity;It can be collectively formed with busbar using it and be fixed between each endpoint of description and busbar
The system structure matrix S of position relationshipSM;SSMIt is a T × NBMatrix, shown in its element assignment principle such as formula (3);
Step 203, direction of tide matrix is built:Although the ring-type multi-source system method of operation and state there are a variety of possibility,
But its most time is in normal operating condition, PQM configurations need to preferentially ensure topological structure under system normal operating condition
The whole network it is considerable;, can be to S by comparing current direction and setting positive directionSMMatrix is modified, and obtains believing comprising system load flow
The direction of tide Matrix C of breathDM;CDMMatrix element cjiAssignment principle such as formula (4) shown in, wherein j is tjNumbering, i is system
Median generatrix is numbered;
Step 204, system list is divided away from block:Consider in system topology, circuit closeness configures monitoring point
Influence;The topological structure for reflecting direction of tide under system normal operating condition is primarily based on, equivalent separation is carried out to system line
Handle, scale of the shortest circuit as cell distance between adjacent busbar in selection topological structure, circuit is equivalent between each busbar
Spacing is integral multiple cell distance, and divides block with cell spacing, is named as single away from block;List is away from unit wires way in block
Amount refers to the number of, lines not comprising power supply and load lead-out wire;
Step 205, considerable matrix structure is monitored:Target grid completes monitor configuration, according under normal operating condition
Trend topological structure, can build the considerable matrix M of monitoringCM, for describing the position relationship between each PQM and each circuit, characterize each prison
Survey the considerable scope of device;MCMMatrix is the matrix of a T × M, and M is monitor quantity, its element mjrAssignment principle such as formula (5)
Shown, wherein r is that PQM is numbered in system;
3rd, the definition of weight factor:Based on PQM arrangement universal criterious, different endpoints have heterogeneity in system, deposit
In different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, enables the property of each endpoint
Embodied in optimal configuration algorithm;
Step 301, PQM preferred arrangements principle is concluded:According to ring-type multi-source network design feature and take PQDS into consideration and determine
Position demand, summarizes 6 PQM arrangement universal criterious;
A1) criterion 1:(N-1) criterion, N are the circuit sum connected on same busbar;
A2) criterion 2:There is the weight of higher on load branch;
A3) criterion 3:Power outlet has highest weight;
A4) criterion 4:Preferentially a PQM is only filled on same circuit;
A5) criterion 5:The endpoint being connected in addition to lead-out wire is loaded with the busbar be less than or equal to two, such endpoint weight compared with
It is low;Correspondingly, the endpoint weight that same circuit is belonged to such endpoint is higher;
A6) criterion 6:Preferentially meet that the whole network is considerable under normal operating condition;
Step 302, valid endpoint carries out computing numbering and attributive classification:To improve efficiency of algorithm, filtered out to all
Valid endpoint carries out computing numbering and attributive classification;The attribute of endpoint is divided into power supply endpoint, load endpoint, opposite end point and common
Endpoint;Opposite end point, which refers to, removes the endpoint that the outer endpoint quantity being connected with busbar of load is less than or equal to 2, and load endpoint includes load
Draw endpoint and belong to the endpoint of same circuit with opposite end point, popular endpoints do not belong to the endpoint of remaining 3 type;Fortune
It is that endpoint presses consideration busbar numbering and S to calculate numberingSMValid endpoint sequence is numbered in matrix;If certain end point is at the same time comprising two
Attribute, then carry out homogenization processing to its property value;
Step 303, weight factor is defined:The explanation of universal criterious based on step 301, different endpoints tools in system
There is heterogeneity, there are different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, makes each
Properties of end vertex can be embodied in optimal configuration algorithm;With reference to endpoint attribute and numbering, α is definedij, βij, γijThree weights
Coefficient quantization characterization PQM configuration universal criterious:βijCharacterize criterion 1;αijCharacterize criterion 2,3,5;γijCharacterize criterion 4,6;αij,
γij, βijMathematic(al) representation respectively as shown in formula (6)~(12);
Wherein,
In formula:I represents busbar numbering;J represents valid endpoint sequence number;aijNumbered for endpoint computing, represent to connect on i busbares
J sequence number valid endpoints;TijRepresent that computing numbering is aijEndpoint corresponding property value when being heterogeneity, such as formula (12) institute
Show;Respectively power supply endpoint, load endpoint, popular endpoints and the corresponding property value of opposite end point, and its
Relation is WithRepresent that busbar i powers on source endpoint, load side respectively
Point, popular endpoints and opposite end point quantity;CiRepresent adding up for the upper all heterogeneity endpoint quantity of busbar i and property value product
With as shown in formula (10);Represent single away from unit line quantity in block k;K represents single away from block number;L represents Dan Ju areas
Number of blocks;DkSingle ratio away from unit line quantity in block k Yu all number of, lines of system is represented, as shown in formula (11);NBij
Represent endpoint aijThe busbar quantity that lower region includes;NiRepresent busbar i upper extreme point total quantitys;AijAll endpoint quantity of expression system
With the ratio of system busbar quantity, as shown in formula (9);BijRepresent endpoint aijThe endpoint quantity that lower region includes;
4th, fuzziness index definition:The monitor quantity arranged in system is inadequate, and part endpoint will be caused effective
Monitoring, endpoint will appear from topology when identifying fuzzy;Define different fuzziness exponential forms and carry out system end fuzzy degree description;
Step 401 is, it is specified that the category that PQM monitorings obscure;
B1 1) is provided:The endpoint for being provided with monitor does not obscure;
B2 2) are provided:The N being connected with same busbariA endpoint, the endpoint quantity for being provided with monitor are more than or wait
In Ni- 1, the then NiA endpoint does not obscure;
B3 3) are provided:Remove and installed when still suffering from two endpoints being connected with busbar after the endpoint of monitor, and if its
Connected endpoint does not obscure, then the two endpoints also do not obscure;
Step 402, ring-type multi-source system monitoring fuzzy indicator definition:Two kinds of fuzziness indexes are defined, system integrally obscures
The maximum fuzziness of degree and part;System entirety fuzziness abbreviation IEA, the maximum fuzziness abbreviation IMA in part;Respectively such as formula (13),
(14) shown in;
Wherein:Represent that k blocks cannot identified circuit end points quantity;bkRepresent the overall fog-level of block k,MkRepresent the monitor quantity being monitored to block k;K is system block sum;[0,1] consolidation is carried out to IMA
Processing, obtains the normalized part maximum fuzziness IMA' as shown in formula (15);
The all consolidations of the value range of IEA and IMA' are to [0,1] section;IEA explanations are when electric energy matter occurs for power grid arbitrary point
During amount disturbance, monitoring system judges the ability of disturbance source position;IMA' illustrates equalization of the monitoring system to each area monitoring ability
Degree;According to the two indexs, configuration N can be calculatedMThe system ambiguous degree of monitoring of a monitoring point, monitors system ambiguous degree referred to as
IA, as shown in formula (16);
Step 403, using the considerable matrix checking fuzziness of monitoring:Each fuzziness index of monitoring system can by monitoring can
See matrix to be obtained and verified, its main flow is:1) the considerable matrix M of system monitoring is builtCM;2) to MCMMatrix is divided
Regional processing;3) each PQM is excluded from each region according to fuzzy regulation 1,2,3;4) b of each block is calculatedk;5) IEA is calculated,
IMA', IA;
The foundation of step 5, PQM Optimal Allocation Models:Based on network characteristic analyze, and eigenmatrix, weight coefficient and
The definition of fuzziness index, establishes the Optimized model of PQM the optimum layouts, object function and constraints respectively such as formula (17),
(18) shown in;
Wherein:F (x) is object function;NMminAnd NMmaxFor the minimum number and maximum quantity of monitoring point;Decision variable xij
Assignment principle such as formula (19) shown in;
6th, the solution of Optimized model:The basic ideas that PQM Optimal Allocation Models solve are to calculate weight factor, base first
Object function is established in it and is optimized under given monitor quantity, determines detector positions;Then, based on known prison
Device position is surveyed, structure monitors considerable matrix, and then calculates IA indexes based on fuzziness testing process;By successive ignition, finally
Determine optimal number and the position of monitor, IA values is minimized under specifying constraint.
Beneficial effects of the present invention are mainly manifested in:1st, for the topological structure feature of ring-type multi-source network, based on pros
Definition to, two-way equivalent circuit and endpoint, and the structure feasible system structure of system structure matrix, direction of tide matrix
The Efficient Characterization of information;2nd, summarize suitable for ring-type multi-source network PQM configure universal criterious, and reasonable definition weight because
The importance degree of sub- all kinds of attribute endpoints of quantization signifying;3rd, target is minimised as with the total weight of system to establish system monitoring point excellent
Change allocation models, and by monitoring the system ambiguous degree of considerable matrix check, filter out optimal monitoring point allocation plan;4th, carried
Method can take into account the high level diagnostics demand of the follow-up progress PQDS positioning of system, and can provide and determine for the stage construction schemes of NPQMS
Plan foundation.
Brief description of the drawings
Fig. 1 is the general frame of the method for the present invention.
Fig. 2 is the definition figure of positive direction, two-way equivalent circuit and endpoint.
Fig. 3 is the computing numbering and attributive classification figure of equivalent network.
Fig. 4 is 14 node system topology diagrams of IEEE.
Fig. 5 is 14 node system trend topology diagrams of IEEE under normal operating condition.
Embodiment
With reference to embodiment and attached drawing, the present invention is described in further detail, but embodiments of the present invention are unlimited
In this.In embodiment a kind of the general frame of the power quality monitor Optimal Configuration Method of meshed distribution network as shown in Figure 1,
Comprise the following steps:
1st, new concept of the definition suitable for ring-type multi-source electric network composition topological analysis:With radial electric network composition feature not
Together, ring-type multi-source power system operating mode and state have changeability, and part circuit current direction is there are a variety of possibility, so as to influence
The monitoring range of each node configuration PQM, it is desirable to account for various possible situations;
Step 101, the regulation of positive direction:For matching PQDS direction determining methods, it is necessary to establish a fixed coordinate system,
Set and line direction is directed toward as positive direction by busbar;
Step 102, the definition of two-way equivalent circuit:By every line equivalent in addition to power supply, load lead-out wire into two
The circuit with fixed current flow direction of parallel bidirectional;
Step 103, the definition of endpoint:Line current inflow end is arranged to the endpoint of configurable PQM to represent the circuit,
I.e. every circuit can equivalently represented be two endpoints;
The topological structure of one typical 6 bus nodes dual power supply power network with circle structure, and its positive direction, two-way equivalent circuit and
The definition of endpoint is as shown in Figure 2;
2nd, the structure of ring-type multi-source system structural characterization matrix:Analyzed based on step 1 ring-network topology, structure is corrected adjacent respectively
Matrix, system structure matrix, direction of tide matrix, the considerable matrix description network structure of monitoring are connect, offer number is distributed rationally for PQM
According to support;
Step 201, the structure of adjacency matrix is corrected:Adjacency matrix AorgIt is usually used in Power Network Topology Analysis, its form isi1,i2=1,2 ..., NBNumbered for bus nodes, NBFor system busbar number of nodes;Matrix element
Assignment principle such as formula (1) shown in;Because load and power supply are all important monitorings pair in networking electric energy quality monitoring system
As, it is necessary to suitably be corrected to adjacency matrix:Represent to load with busbar and the connecting line being equivalent between the earth of busbar 0, use
Busbar represents power supply with the connecting line of its own;Correct adjacency matrixForm be To increase
The system node numbering after the busbar 0 for representing the earth is added;ElementRepresent that the structure of each endpoint in ring-type multi-source network is compiled
Number, its value represents that the endpoint whether there is, shown in assignment principle such as formula (2);The adjoining square of 6 node example systems in step 1
Battle array AorgWith amendment adjacency matrixAs shown in formula (20);
Step 202, system structure matrix is built:WillAll valid endpoints in matrix, i.e.,Endpoint, press
Row order extracts, and it is sequentially labeled as to new system end point symbol t successivelyj, wherein j=1,2 ... T is successively sequentially
T after markjNumbering, T is valid endpoint total quantity;It can be collectively formed with busbar using it and be fixed between each endpoint of description and busbar
The system structure matrix S of position relationshipSM;SSMIt is a T × NBMatrix, shown in its element assignment principle such as formula (3);Step 1
In 6 node example systems system structure matrix SSMAs shown in formula (21);
Step 203, direction of tide matrix is built:Although the ring-type multi-source system method of operation and state there are a variety of possibility,
But its most time is in normal operating condition, PQM configurations need to preferentially ensure topological structure under system normal operating condition
The whole network it is considerable;, can be to S by comparing current direction and setting positive directionSMMatrix is modified, and obtains believing comprising system load flow
The direction of tide Matrix C of breathDM;CDMMatrix element cjiAssignment principle such as formula (4) shown in, wherein j is tjNumbering, i is system
Median generatrix is numbered;The direction of tide Matrix C of 6 node example systems in step 1DMAs shown in formula (22);
Step 204, system list is divided away from block:Consider in system topology, circuit closeness configures monitoring point
Influence;The topological structure for reflecting direction of tide under system normal operating condition is primarily based on, equivalent separation is carried out to system line
Handle, scale of the shortest circuit as cell distance between adjacent busbar in selection topological structure, circuit is equivalent between each busbar
Spacing is integral multiple cell distance, and divides block with cell spacing, is named as single away from block;List is away from unit wires way in block
Amount refers to the number of, lines not comprising power supply and load lead-out wire;
Step 205, considerable matrix structure is monitored:Target grid completes monitor configuration, according under normal operating condition
Trend topological structure, can build the considerable matrix M of monitoringCM, for describing the position relationship between each PQM and each circuit, characterize each prison
Survey the considerable scope of device;MCMMatrix is the matrix of a T × M, and M is monitor quantity, its element mjrAssignment principle such as formula (5)
Shown, wherein r is that PQM is numbered in system;
3rd, the definition of weight factor:Based on PQM arrangement universal criterious, different endpoints have heterogeneity in system, deposit
In different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, enables the property of each endpoint
Embodied in optimal configuration algorithm;
Step 301, PQM preferred arrangements principle is concluded:According to ring-type multi-source network design feature and take PQDS into consideration and determine
Position demand, summarizes 6 PQM arrangement universal criterious;
A1) criterion 1:(N-1) criterion, N are the circuit sum connected on same busbar;
A2) criterion 2:There is the weight of higher on load branch;
A3) criterion 3:Power outlet has highest weight;
A4) criterion 4:Preferentially a PQM is only filled on same circuit;
A5) criterion 5:The endpoint being connected in addition to lead-out wire is loaded with the busbar be less than or equal to two, such endpoint weight compared with
It is low;Correspondingly, the endpoint weight that same circuit is belonged to such endpoint is higher;
A6) criterion 6:Preferentially meet that the whole network is considerable under normal operating condition;
Step 302, valid endpoint carries out computing numbering and attributive classification:To improve efficiency of algorithm, filtered out to all
Valid endpoint carries out computing numbering and attributive classification;The attribute of endpoint is divided into power supply endpoint, load endpoint, opposite end point and common
Endpoint;Opposite end point, which refers to, removes the endpoint that the outer endpoint quantity being connected with busbar of load is less than or equal to 2, and load endpoint includes load
Draw endpoint and belong to the endpoint of same circuit with opposite end point, popular endpoints do not belong to the endpoint of remaining 3 type;Fortune
It is that endpoint presses consideration busbar numbering and S to calculate numberingSMValid endpoint sequence is numbered in matrix;If certain end point is at the same time comprising two
Attribute, then carry out homogenization processing to its property value;The equivalent network computing numbering and category of 6 node example systems in step 1
Property classification it is as shown in Figure 3;
Step 303, weight factor is defined:The explanation of universal criterious based on step 301, different endpoints have in system
Heterogeneity, there are different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, makes each end
Point property can be embodied in optimal configuration algorithm;With reference to endpoint attribute and numbering, α is definedij, βij, γijThree weight systems
Quantification characterization PQM configuration universal criterious:βijCharacterize criterion 1;αijCharacterize criterion 2,3,5;γijCharacterize criterion 4,6;αij,
γij, βijMathematic(al) representation respectively as shown in formula (6)~(12);
The topological structure of one 14 node ring-type multi-source system of IEEE is as shown in Figure 4;Through simulation study, its normal operation
System load flow topological structure under state is as shown in Figure 5;Based on step 1,2 system structure information description and the weight of step 3
Because of sub-definite, each total weighted value δ of endpointij=αij·βij·γijResult of calculation it is as shown in table 1;
The total weighted value δ of 1 each endpoint of tableijResult of calculation
4th, fuzziness index definition:The monitor quantity arranged in system is inadequate, and part endpoint will be caused effective
Monitoring, endpoint will appear from topology when identifying fuzzy;Define different fuzziness exponential forms and carry out system end fuzzy degree description;
Step 401 is, it is specified that the category that PQM monitorings obscure;
B1 1) is provided:The endpoint for being provided with monitor does not obscure;
B2 2) are provided:The N being connected with same busbariA endpoint, the endpoint quantity for being provided with monitor are more than or wait
In Ni- 1, the then NiA endpoint does not obscure;
B3 3) are provided:Remove and installed when still suffering from two endpoints being connected with busbar after the endpoint of monitor, and if its
Connected endpoint does not obscure, then the two endpoints also do not obscure;
Step 402, ring-type multi-source system monitoring fuzzy indicator definition:Two kinds of fuzziness indexes are defined, system integrally obscures
The maximum fuzziness of degree and part;System entirety fuzziness abbreviation IEA, the maximum fuzziness abbreviation IMA in part;Respectively such as formula (13),
(14) shown in;[0,1] consolidation processing is carried out to IMA, is obtained shown in IMA' such as formulas (15);
The all consolidations of the value range of IEA and IMA' are to [0,1] section;IEA explanations are when electric energy matter occurs for power grid arbitrary point
During amount disturbance, monitoring system judges the ability of disturbance source position;IMA' illustrates equalization of the monitoring system to each area monitoring ability
Degree;According to the two indexs, configuration N can be calculatedMThe system ambiguous degree of monitoring of a monitoring point, monitors system ambiguous degree referred to as
IA, as shown in formula (16);
Step 403, using the considerable matrix checking fuzziness of monitoring:Each fuzziness index of monitoring system can by monitoring can
See matrix to be obtained and verified, its main flow is:1) the considerable matrix M of system monitoring is builtCM;2) to MCMMatrix is divided
Regional processing;3) each PQM is excluded from each region according to fuzzy regulation 1,2,3;4) b of each block is calculatedk;5) IEA is calculated,
IMA', IA;
5th, the foundation of PQM Optimal Allocation Models:Analyzed based on network characteristic, and eigenmatrix, weight coefficient and fuzzy
The definition of index is spent, establishes the Optimized model of PQM the optimum layouts, object function and constraints are respectively such as formula (17), (18) institute
Show;Decision variable xijAssignment principle such as formula (19) shown in;
6th, the solution of Optimized model:The basic ideas that PQM Optimal Allocation Models solve are to calculate weight factor, base first
Object function is established in it and is optimized under given monitor quantity, determines detector positions;Then, based on known prison
Device position is surveyed, structure monitors considerable matrix, and then calculates IA indexes based on fuzziness testing process;By successive ignition, finally
Determine optimal number and the position of monitor, IA values is minimized under specifying constraint.
For the node systems of IEEE 14 shown in Fig. 4, the system ambiguous degree index of weight coefficient and step 4 based on step 3
Computational methods, system PQM Optimal Allocation Models are established using step 5, and the PQM optimizations of step 6 are carried out using Matlab softwares
Allocation models solves, and it is { t to obtain the optimal PQM configuration results of system11,t12,t22,t23,t24,t33,t34,t47,t52,t54,
t60,t66,t6.12,t88,t94,t9.10,t9.14,t10.0,t11.0,t12.0,t13.0,t13.6};System monitoring fuzzy indicator IA=at this time
0, it can be achieved that PQDS is pin-pointed on specific circuit;When the PQM quantity that system can be provided be less than 22, will appear from monitoring
It is fuzzy;Equally spacedly choose different monitoring points quantity and optimize configuration, obtained PQM the optimum layouts result and corresponding monitoring system
Fuzziness of uniting IA is as shown in table 2.
IA and the optimum layout scheme during 2 PQM lazy weights of table
Sample calculation analysis shows that institute's extracting method of the present invention can carry out the effective table of topological structure of ring-type multiple source power distribution network
Sign, and the definition of the weight factor based on different attribute endpoint and different system fuzzy targets, structure majorized function realize mesh
Mark system PQM is distributed rationally;Particularly, institute's extracting method can take into account the power quality high-grade intelligent of the follow-up progress PQDS positioning of system
Diagnostic requirements.
As described above, the present invention can be better realized, above-described embodiment is only the exemplary embodiments of the present invention, is not used
To limit the practical range of the present invention, i.e., all equivalent changes and modifications made according to present invention all will for right of the present invention
Scope claimed is asked to be covered.
Claims (1)
1. a kind of power quality monitor Optimal Configuration Method of ring-type multiple source power distribution net, includes the following steps:
Step 1, definition are suitable for the new concept of ring-type multi-source electric network composition topological analysis:With radial electric network composition feature not
Together, ring-type multi-source power system operating mode and state have changeability, and part circuit current direction is there are a variety of possibility, so as to influence
The monitoring range of each node configuration power quality monitor, it is desirable to account for various possible situations;Electric energy quality monitoring
Device abbreviation PQM;
Step 101, the regulation of positive direction:For matching power quality disturbance direction determining method, it is necessary to establish a fixed seat
Mark system, sets and is directed toward line direction as positive direction by busbar;Power quality disturbance abbreviation PQDS;
Step 102, the definition of two-way equivalent circuit:Every line equivalent in addition to power supply, load lead-out wire is parallel into two
The two-way circuit with fixed current flow direction;
Step 103, the definition of endpoint:Line current inflow end is arranged to the endpoint of configurable PQM to represent the circuit, i.e., often
Bar line equivalent is expressed as two endpoints;
The structure of step 2, ring-type multi-source system structural characterization matrix:Analyzed based on step 1 ring-network topology, structure is corrected adjacent respectively
Matrix, system structure matrix, direction of tide matrix, the considerable matrix description network structure of monitoring are connect, offer number is distributed rationally for PQM
According to support;
Step 201, the structure of adjacency matrix is corrected:Adjacency matrix AorgIt is usually used in Power Network Topology Analysis, its form isNumbered for bus nodes, NBFor system busbar number of nodes;Matrix element's
Shown in assignment principle such as formula (1);Because load and power supply are all important monitoring objects in networking electric energy quality monitoring system,
Need suitably to correct adjacency matrix:Represent to load with busbar and the connecting line being equivalent between the earth of busbar 0, use busbar
Power supply is represented with the connecting line of its own;Correct adjacency matrixForm be To add
Represent the system node numbering after the busbar 0 of the earth;ElementRepresent the structure number of each endpoint in ring-type multi-source network,
Its value represents that the endpoint whether there is, shown in assignment principle such as formula (2);
Step 202, system structure matrix is built:WillAll valid endpoints in matrix, i.e.,Endpoint, by row order
Extract, and it is sequentially labeled as to new system end point symbol t successivelyj, wherein j=1,2 ... after T is sequentially marks successively
tjNumbering, T is valid endpoint total quantity;Fixed position relation between each endpoint of description and busbar is collectively formed using itself and busbar
System structure matrix SSM;SSMIt is a T × NBMatrix, shown in its element assignment principle such as formula (3);
Step 203, direction of tide matrix is built:Although the ring-type multi-source system method of operation and state there are a variety of possibility, its
Most times are in normal operating condition, and PQM configurations need to preferentially ensure the complete of topological structure under system normal operating condition
Net considerable;By comparing current direction and setting positive direction, to SSMMatrix is modified, and obtains including the tide of system load flow information
Flow direction matrix CDM;CDMMatrix element cjiAssignment principle such as formula (4) shown in, wherein j is tjNumbering, i is system median generatrix
Numbering;
Step 204, system list is divided away from block:Consider in system topology, the influence that circuit closeness configures monitoring point;
The topological structure for reflecting direction of tide under system normal operating condition is primarily based on, equivalent separation processing is carried out to system line,
Choose in topological structure scale of the shortest circuit as cell distance between adjacent busbar, the equivalent separation of circuit is between each busbar
Integral multiple cell distance, and block is divided with cell spacing, it is named as single away from block;It is single to refer to not away from unit line quantity in block
Number of, lines comprising power supply and load lead-out wire;
Step 205, considerable matrix structure is monitored:Target grid completes monitor configuration, according to the trend under normal operating condition
Topological structure, structure monitor considerable matrix MCM, for describing the position relationship between each PQM and each circuit, characterize each monitor
Considerable scope;MCMMatrix is the matrix of a T × M, and M is monitor quantity, its element mjrShown in assignment principle such as formula (5), its
Middle r is that PQM is numbered in system;
The definition of step 3, weight factor:Based on PQM arrangement universal criterious, different endpoints have heterogeneity in system, deposit
In different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, enables the property of each endpoint
Embodied in optimal configuration algorithm;
Step 301, PQM preferred arrangements principle is concluded:According to ring-type multi-source network design feature and take PQDS positioning need into consideration
Ask, summarize 6 PQM arrangement universal criterious;
A1) criterion 1:(N-1) criterion, N are the circuit sum connected on same busbar;
A2) criterion 2:There is the weight of higher on load branch;
A3) criterion 3:Power outlet has highest weight;
A4) criterion 4:Preferentially a PQM is only filled on same circuit;
A5) criterion 5:The endpoint being connected in addition to lead-out wire is loaded with the busbar is less than or equal to two, and such endpoint weight is relatively low;Phase
Ying Di, the endpoint weight that same circuit is belonged to such endpoint are higher;
A6) criterion 6:Preferentially meet that the whole network is considerable under normal operating condition;
Step 302, valid endpoint carries out computing numbering and attributive classification:For improve efficiency of algorithm, to it is all filter out it is effective
Endpoint carries out computing numbering and attributive classification;The attribute of endpoint is divided into power supply endpoint, load endpoint, opposite end point and common end
Point;Opposite end point, which refers to, removes the endpoint that the outer endpoint quantity being connected with busbar of load is less than or equal to 2, and load endpoint draws including load
Go out endpoint and belong to the endpoint of same circuit with opposite end point, popular endpoints do not belong to the endpoint of remaining 3 type;Computing
Numbering is endpoint by consideration busbar numbering and SSMValid endpoint sequence is numbered in matrix;If certain end point is at the same time comprising two kinds
Attribute, then carry out homogenization processing to its property value;
Step 303, weight factor is defined:The explanation of universal criterious based on step 301, different endpoints have difference in system
Property, there are different significance levels;The significance level of each endpoint of reasonable definition weight factor quantificational description, makes each endpoint
Mass-energy is embodied in optimal configuration algorithm;With reference to endpoint attribute and numbering, α is definedij, βij, γijThree weight coefficient amounts
Change characterization PQM configuration universal criterious:βijCharacterize criterion 1;αijCharacterize criterion 2,3,5;γijCharacterize criterion 4,6;αij, γij,
βijMathematic(al) representation respectively as shown in formula (6)~(12);
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</munderover>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>N</mi>
<mn>1</mn>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>N</mi>
<mn>2</mn>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<msub>
<mi>t</mi>
<mn>3</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>N</mi>
<mn>3</mn>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<msub>
<mi>t</mi>
<mn>4</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>N</mi>
<mn>4</mn>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>N</mi>
<mi>D</mi>
<mi>k</mi>
</msubsup>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<msubsup>
<mi>N</mi>
<mi>D</mi>
<mi>k</mi>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:I represents busbar numbering;J represents valid endpoint sequence number;aijNumbered for endpoint computing, represent the j connected on i busbares
Sequence number valid endpoint;TijRepresent that computing numbering is aijEndpoint corresponding property value when being heterogeneity, as shown in formula (12);Respectively power supply endpoint, load endpoint, popular endpoints and the corresponding property value of opposite end point, and its relation
For WithRepresent that busbar i powers on source endpoint, load endpoint, general respectively
Go side point and opposite end point quantity;CiRepresent the upper all heterogeneity endpoint quantity of busbar i and property value product add up and, such as
Shown in formula (10);Represent single away from unit line quantity in block k;K represents single away from block number;L represents single away from block counts
Amount;DkSingle ratio away from unit line quantity in block k Yu all number of, lines of system is represented, as shown in formula (11);NBijRepresent
Endpoint aijThe busbar quantity that lower region includes;NiRepresent busbar i upper extreme point total quantitys;AijAll endpoint quantity of expression system are with being
The ratio for busbar quantity of uniting, as shown in formula (9);BijRepresent endpoint aijThe endpoint quantity that lower region includes;
Step 4, fuzziness index definition:The monitor quantity arranged in system is inadequate, and part endpoint will be caused effective
Monitoring, endpoint will appear from topology when identifying fuzzy;Define different fuzziness exponential forms and carry out system end fuzzy degree description;
Step 401 is, it is specified that the category that PQM monitorings obscure;
B1 1) is provided:The endpoint for being provided with monitor does not obscure;
B2 2) are provided:The N being connected with same busbariA endpoint, the endpoint quantity for being provided with monitor are greater than or equal to Ni-
1, the then NiA endpoint does not obscure;
B3 3) are provided:Remove and installed when still suffering from two endpoints being connected with busbar after the endpoint of monitor, if coupled
Endpoint do not obscure, then the two endpoints also do not obscure;
Step 402, ring-type multi-source system monitoring fuzzy indicator definition:Define two kinds of fuzziness indexes, system entirety fuzziness and
The maximum fuzziness in part;System entirety fuzziness abbreviation IEA, the maximum fuzziness abbreviation IMA in part;Respectively such as formula (13), (14)
It is shown;
<mrow>
<mi>I</mi>
<mi>E</mi>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>N</mi>
<mrow>
<mi>N</mi>
<mi>D</mi>
</mrow>
<mi>k</mi>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</munderover>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>I</mi>
<mi>M</mi>
<mi>A</mi>
<mo>&equiv;</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mi>k</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>K</mi>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein:Represent that k blocks cannot identified circuit end points quantity;bkRepresent the overall fog-level of block k,MkRepresent the monitor quantity being monitored to block k;K is system block sum;[0,1] consolidation is carried out to IMA
Processing, obtains the normalized part maximum fuzziness IMA' as shown in formula (15);
<mrow>
<msup>
<mi>IMA</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>N</mi>
<mrow>
<mi>N</mi>
<mi>D</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>&times;</mo>
<mi>I</mi>
<mi>M</mi>
<mi>A</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>N</mi>
<mrow>
<mi>N</mi>
<mi>D</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
The all consolidations of the value range of IEA and IMA' are to [0,1] section;IEA explanations are disturbed when power quality occurs for power grid arbitrary point
When dynamic, monitoring system judges the ability of disturbance source position;IMA' illustrates impartial degree of the monitoring system to each area monitoring ability;
According to the two indexs, configuration N is calculatedMThe system ambiguous degree of monitoring of a monitoring point, monitors system ambiguous degree abbreviation IA, such as formula
(16) shown in;
Step 403, using the considerable matrix checking fuzziness of monitoring:Each fuzziness index of monitoring system is by monitoring considerable matrix
Obtained and verified, its main flow is:1) the considerable matrix M of system monitoring is builtCM;2) to MCMMatrix is carried out at subregion
Reason;3) each PQM is excluded from each region according to fuzzy regulation 1,2,3;4) b of each block is calculatedk;5) IEA, IMA', IA are calculated;
The foundation of step 5, PQM Optimal Allocation Models:Analyzed based on network characteristic, and eigenmatrix, weight coefficient and fuzzy
The definition of index is spent, establishes the Optimized model of PQM the optimum layouts, object function and constraints are respectively such as formula (17), (18) institute
Show;
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>min</mi>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</munder>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</munderover>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>NM</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</munderover>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</munderover>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>NM</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein:F (x) is object function;NMminAnd NMmaxFor the minimum number and maximum quantity of monitoring point;Decision variable xijTax
It is worth shown in principle such as formula (19);
The solution of step 6, Optimized model:The basic ideas that PQM Optimal Allocation Models solve are to calculate weight factor, base first
Object function is established in it and is optimized under given monitor quantity, determines detector positions;Then, based on known prison
Device position is surveyed, structure monitors considerable matrix, and then calculates IA indexes based on fuzziness testing process;By successive ignition, finally
Determine optimal number and the position of monitor, IA values is minimized under specifying constraint.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610393772.5A CN105977968B (en) | 2016-06-06 | 2016-06-06 | Optimal configuration method for power quality monitor of annular multi-source power distribution network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610393772.5A CN105977968B (en) | 2016-06-06 | 2016-06-06 | Optimal configuration method for power quality monitor of annular multi-source power distribution network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105977968A CN105977968A (en) | 2016-09-28 |
CN105977968B true CN105977968B (en) | 2018-05-08 |
Family
ID=57010830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610393772.5A Active CN105977968B (en) | 2016-06-06 | 2016-06-06 | Optimal configuration method for power quality monitor of annular multi-source power distribution network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105977968B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107181254B (en) * | 2016-11-23 | 2020-04-14 | 国电南瑞科技股份有限公司 | Matrix type quantization index calculation and alternating current-direct current coordination control method for multiple stability problems |
CN108964282B (en) * | 2018-08-28 | 2020-09-29 | 西门子电力自动化有限公司 | Method and device for deploying power quality monitoring device and computer readable medium |
CN113691433B (en) * | 2021-08-31 | 2023-03-03 | 深圳云天励飞技术股份有限公司 | Data transmission system, method, device, electronic equipment and storage medium |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101232182A (en) * | 2008-01-18 | 2008-07-30 | 清华大学 | Three-dimensional coordinated electric network energy managing system and method for controlling and evaluating electric network |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9336338B2 (en) * | 2012-03-07 | 2016-05-10 | Siemens Aktiengesellschaft | Global solutions of smart building-grid energy management models |
-
2016
- 2016-06-06 CN CN201610393772.5A patent/CN105977968B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101232182A (en) * | 2008-01-18 | 2008-07-30 | 清华大学 | Three-dimensional coordinated electric network energy managing system and method for controlling and evaluating electric network |
Non-Patent Citations (1)
Title |
---|
基于状态估计的电能质量监测点优化配置研究;吴丽丽等;《机电工程》;20160331;第33卷(第3期);第331-335页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105977968A (en) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106505593B (en) | A kind of analysis of distribution transforming three-phase imbalance and the method for load adjustment based on big data | |
CN109274095B (en) | Mutual information-based low-voltage distribution network user topology estimation method and system | |
Jing et al. | Distributed or centralized? Designing district-level urban energy systems by a hierarchical approach considering demand uncertainties | |
CN104835018B (en) | The reliability estimation method and device of medium voltage distribution network rack | |
CN104573947A (en) | Comprehensive evaluation method for low-voltage transformer areas of regional intelligent distribution network | |
CN105977968B (en) | Optimal configuration method for power quality monitor of annular multi-source power distribution network | |
CN105354680B (en) | GIS-based power grid section auxiliary determination analysis method | |
Alam et al. | A SAX-based advanced computational tool for assessment of clustered rooftop solar PV impacts on LV and MV networks in smart grid | |
CN102623987A (en) | Multiple-DC (direct current)-droppoint selection method based on multiple feed-in short circuit ratios | |
CN103996147A (en) | Comprehensive evaluation method for power distribution network | |
Qi et al. | A fuzzy comprehensive evaluation and entropy weight decision-making based method for power network structure assessment | |
CN105633948A (en) | Random fuzzy power flow algorithm for distributed wind power, photovoltaic power generation and other uncertain energy sources connected to power system | |
CN104573848A (en) | Power demand prediction and planning and reliability-based power distribution network construction method | |
WO2023165348A1 (en) | Method and device for determining switching state of reactive compensation of distribution transformer | |
CN112329973A (en) | Space-time load prediction method based on graph neural network and regional gridding | |
CN106786529A (en) | A kind of distributed static security analysis method | |
Ma et al. | Economic operation evaluation of active distribution network based on fuzzy Borda method | |
Wang et al. | A two-step load disaggregation algorithm for quasi-static time-series analysis on actual distribution feeders | |
Miraftabzadeh et al. | K-means and alternative clustering methods in modern power systems | |
CN106022957A (en) | Power grid coordinated development evaluation method for power system | |
CN104867063A (en) | Power grid multiple major failure recognition method based on a linear weighting method | |
CN107732902A (en) | Power distribution network economical operation monitoring and evaluation method | |
CN105224993A (en) | A kind of substation planning method optimized power distribution network and run | |
Hasan et al. | A new smart approach of an efficient energy consumption management by using a machinelearning technique | |
Zhang et al. | Distributed Architecture of Power Grid Asset Management and Future Research Directions |
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 |