CN114970072A - Partition scheme evaluation method and system applied to power distribution network simulation - Google Patents
Partition scheme evaluation method and system applied to power distribution network simulation Download PDFInfo
- Publication number
- CN114970072A CN114970072A CN202110219443.XA CN202110219443A CN114970072A CN 114970072 A CN114970072 A CN 114970072A CN 202110219443 A CN202110219443 A CN 202110219443A CN 114970072 A CN114970072 A CN 114970072A
- Authority
- CN
- China
- Prior art keywords
- simulation
- distribution network
- parallel
- power distribution
- partition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 331
- 238000005192 partition Methods 0.000 title claims abstract description 173
- 238000009826 distribution Methods 0.000 title claims abstract description 160
- 238000011156 evaluation Methods 0.000 title claims abstract description 125
- 238000004364 calculation method Methods 0.000 claims abstract description 150
- 238000000034 method Methods 0.000 claims abstract description 72
- 238000000638 solvent extraction Methods 0.000 claims abstract description 46
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims description 36
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000001133 acceleration Effects 0.000 claims description 27
- 238000004891 communication Methods 0.000 claims description 24
- 238000003860 storage Methods 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000013215 result calculation Methods 0.000 description 4
- 101000912561 Bos taurus Fibrinogen gamma-B chain Proteins 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 235000008446 instant noodles Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method and a system for evaluating a partition scheme applied to power distribution network simulation, wherein the method comprises the following steps: acquiring structural parameters of a power distribution network simulation partitioning scheme, a parallel simulation result after partitioning of the power distribution network and a serial simulation result of the power distribution network; calculating each evaluation index value for evaluating the power distribution network partition scheme by using the structural parameters, the serial simulation result and the parallel simulation result; calculating scores of the power distribution network partition schemes based on the evaluation index values and weights corresponding to the evaluation indexes predetermined by an analytic hierarchy process, and obtaining the quality grades of the power distribution network partition schemes based on the scores and preset partition standards; wherein the evaluation index comprises: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision. The invention is beneficial to the selection of the optimal partition scheme of the power distribution network simulation.
Description
Technical Field
The invention relates to the field of power distribution network simulation partitioning, in particular to an evaluation method and system of a partitioning scheme applied to power distribution network simulation.
Background
The power distribution network has numerous devices, complex structure, large scale and complex analysis and calculation. Especially, with the development of the smart power grid, a large number of distributed power sources and different types of loads are connected to the power distribution network, so that the simulation analysis and calculation of the power distribution network become more complex. The simulation calculation scale can be greatly reduced through network partition, and a large amount of calculation time is saved. However, different partition methods consume different computer resources, running time and calculation accuracy when performing subsequent calculations. The better partitioning mode is to make the calculation scale of each partitioned area similar, otherwise, computer resources are wasted due to mutual waiting, so that the parallel simulation partitioning standard in the prior art is difficult to accurately obtain the optimal partitioning scheme, and the accuracy of the power distribution network simulation partitioning and the precision and efficiency of the power distribution network in simulation analysis calculation cannot be further improved.
Disclosure of Invention
Aiming at the problem that the optimal partition scheme of the power distribution network is difficult to accurately obtain in the prior art, the invention provides an evaluation method of the partition scheme applied to power distribution network simulation, which comprises the following steps:
acquiring structural parameters of a power distribution network simulation partitioning scheme, a parallel simulation result after partitioning of the power distribution network and a serial simulation result of the power distribution network;
calculating each evaluation index value for evaluating the power distribution network partition scheme by using the structural parameters, the serial simulation result and the parallel simulation result;
calculating scores of the power distribution network partition schemes based on the evaluation index values and weights corresponding to the evaluation indexes predetermined by an analytic hierarchy process, and obtaining quality grades of the power distribution network partition schemes based on the scores and preset partition standards;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
Preferably, the calculating an evaluation result of the power distribution network simulation partitioning scheme based on the evaluation index values and the weights corresponding to the evaluation indexes acquired in advance includes:
summing the products of the evaluation index values obtained by calculation and the weights corresponding to the evaluation index values to obtain the evaluation scores of the power distribution network simulation partition scheme;
and classifying according to preset grades based on the evaluation scores of the power distribution network simulation partition scheme to obtain the evaluation result of the power distribution network simulation partition scheme.
Preferably, the evaluation score of the power distribution network simulation partitioning scheme is calculated according to the following formula:
Y＝ω 1 y 1 +ω 2 y 2 +ω 3 y 3 +ω 4 y 4 +ω 5 y 5 +ω 6 y 6 +ω 7 y 7
in the formula, Y is the evaluation score of the power distribution network simulation partition scheme, omega 1 Simulating weight values, y, of the computation load balance for the nodes in each partition 1 Simulating a computation of a metric balance, ω, for nodes in each partition 2 Dividing the weight value of degree for the feeder node, y 2 Degree of division, omega, for feeder nodes 3 Weight value of degree of division for substation feeder node, y 3 Degree of division, omega, for substation feeder nodes 4 Computing weight values, y, of acceleration ratios for parallel simulations 4 Computing acceleration ratio, omega, for parallel simulation 5 Weight value of parallel simulation computing efficiency, y 5 For parallel simulation computational efficiency, omega 6 Weight of expansibility of power distribution network in parallel simulationValue y 6 For the expansibility, omega, of distribution networks in parallel simulation 7 Computing weight values of precision, y, for parallel simulations 7 The accuracy is calculated for the parallel simulation.
Preferably, the structural parameters of the power distribution network simulation partition scheme include:
the number of PQ nodes, the number of PQ node loads, the number of PV nodes, the number of loads of PV node equivalent PQ nodes, the number of ZIP nodes, the number of loads of ZIP node equivalent PQ nodes, the number of partitions, the node numbers of all feeders and the node numbers of all substations in the power distribution network simulation partitioning scheme;
the serial simulation result comprises: analyzing and calculating the required time and the voltage result of each node under serial simulation by serial simulation;
the parallel simulation result comprises the following steps: the method comprises the steps of time required by each partition for performing parallel simulation analysis and calculation, overhead time such as communication and storage required by parallel simulation calculation and analysis, cost for communication and synchronization in parallel calculation, total time required by parallel calculation, and voltage results of each node under parallel simulation.
Preferably, the balance degree y of the simulation calculated amount of all the partition nodes 1 Calculated as follows:
in the formula, n cal,i The equivalent PQ node load number of the ith partition;
wherein the equivalent PQ node load number n of the ith partition cal,i Determined as follows:
n cal,i ＝N PQ,i ×n cal,PQ +N PV,i ×n cal,PV +N ZIP,i ×n cal,ZIP
in the formula, N PQ,i The number of PQ nodes in the ith partition, n cal,PQ Is the PQ node load number, N PV,i The number of PV nodes in the ith partition, n cal,PV Equivalent PQ node load for PV nodeNumber, N ZIP,i Is the number of ZIP nodes in the ith partition, n cal,ZIP The load number of the PQ node is equivalent to that of the ZIP node.
Preferably, the parallel simulation calculates an acceleration ratio y 4 Calculated as follows:
in the formula, T ∞ Time required for serial simulation analysis calculation, T i Time required for parallel simulation analysis calculation for ith partition, T t Overhead time such as communication and storage required for parallel simulation calculation and analysis, wherein i is the number of partitions;
preferably, the parallel simulation has a computational efficiency y 5 Determined as follows:
preferably, the expansibility y of the power distribution network in parallel simulation 6 Calculated as follows:
in the formula, T 0 Charges for communication and synchronization in parallel computing, T 0 ' total time required for parallel computation;
wherein the total time T required for the parallel computation 0 ', calculated as:
preferably, the parallel simulation computing precision y 7 Calculated as follows:
y 7 ＝1-η
in the formula, eta is the average percentage of voltage errors between the serial simulation and the parallel simulation;
wherein the average percentage η of voltage errors between the serial simulation and the parallel simulation is calculated as follows:
in the formula, n is the number of nodes, U par,m 、U ser,m The voltage results of the mth node under parallel and serial simulations, respectively.
Preferably, the determining of the weight corresponding to each evaluation index includes:
s1: evaluating the importance of each evaluation index by adopting a reciprocity scaling method in an analytic hierarchy process, and forming a judgment matrix;
s2: carrying out consistency check on the judgment matrix, and entering S3 if the consistency check requirement is met, or entering S1 to adjust the value of an element in the judgment matrix;
s3: and calculating the index weight of each evaluation index based on the judgment matrix passing the consistency test.
Preferably, the index weight calculation formula under the evaluation index system is as follows:
in the formula, ω a The index weight of the a-th index, v a A characteristic vector of the a index; v. of k Is the feature vector of the k index, and m is the evaluation index.
Based on the same invention concept, the invention also provides an evaluation system of the partition scheme applied to the power distribution network simulation, which comprises the following steps:
the acquisition module is used for acquiring the structural parameters of the power distribution network simulation partitioning scheme, the parallel simulation result after the power distribution network is partitioned and the serial simulation result of the power distribution network;
the calculation module is used for calculating each evaluation index value for evaluating the power distribution network partition scheme by utilizing the structural parameters, the serial simulation result and the parallel simulation result;
the evaluation module is used for calculating the scores of the power distribution network partition schemes based on the evaluation index values and the weights corresponding to the evaluation indexes predetermined by utilizing an analytic hierarchy process, and obtaining the quality grades of the power distribution network partition schemes based on the scores and the preset partition standards;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention provides a method and a system for evaluating a partition scheme applied to power distribution network simulation, wherein the method comprises the following steps: acquiring structural parameters of a power distribution network simulation partitioning scheme, a parallel simulation result after partitioning of the power distribution network and a serial simulation result of the power distribution network; calculating each evaluation index value for evaluating the power distribution network partition scheme by using the structural parameters, the serial simulation result and the parallel simulation result; calculating scores of the power distribution network partition schemes based on the evaluation index values and weights corresponding to the evaluation indexes predetermined by an analytic hierarchy process, and obtaining quality grades of the power distribution network partition schemes based on the scores and preset partition standards; wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision. The method is beneficial to guiding the selection of the optimal scheme of the simulation partition of the power distribution network.
2. In the partition scheme for comprehensively evaluating the power distribution network simulation, whether each partitioned partition can make the calculation scale (time scale) of each iterative operation as small as possible in the parallel simulation calculation or not is judged, namely, a larger parallel calculation acceleration ratio is obtained; whether the precision of the calculation result after partitioning is approximately the same as that before partitioning is evaluated so as to ensure the accuracy of the calculation result; meanwhile, considering the topology of a communication network, evaluating whether the feeders of the shared communication channel are distributed in the same subarea as much as possible; and whether the feeder lines of the same transformer substation are distributed in the same partition as much as possible or not is judged, so that the calculation speed and the calculation precision of the parallel simulation calculation analysis of the power distribution network are improved.
Drawings
FIG. 1 is a schematic diagram of an evaluation method of a partitioning scheme applied to power distribution network simulation according to the present invention;
FIG. 2 is a schematic diagram of a power distribution network simulation partition structure according to the present invention;
fig. 3 is a schematic diagram of an evaluation system of a partitioning scheme applied to power distribution network simulation according to the present invention.
Detailed Description
Aiming at the problems in the prior art, the invention provides an evaluation method of a partition scheme applied to power distribution network simulation, as shown in fig. 1, the evaluation method comprises the following steps:
step 2, calculating each evaluation index value for evaluating the power distribution network partition scheme by using the structural parameters, the serial simulation result and the parallel simulation result;
step 3, calculating scores of the power distribution network partition schemes based on the evaluation index values and weights corresponding to the evaluation indexes predetermined by an analytic hierarchy process, and combining the scores with preset partition standards to obtain the quality grades of the power distribution network partition schemes;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
In this embodiment, the distribution network includes two or more distribution areas and their subsidiary lines; the parallel simulation analysis refers to dividing the power distribution network into a plurality of areas and simultaneously performing simulation analysis calculation on the plurality of areas when the power distribution network is subjected to simulation analysis, as shown in fig. 2.
In step 1, obtaining the structural parameters of the power distribution network simulation partitioning scheme includes: the number of PQ nodes, the number of PQ node loads, the number of PV nodes, the number of loads of PV node equivalent PQ nodes, the number of ZIP nodes, the number of loads of ZIP node equivalent PQ nodes, the number of partitions, the node numbers of all feeders and the node numbers of all substations in the power distribution network simulation partitioning scheme;
acquiring the serial simulation result comprises: analyzing and calculating the required time and the voltage result of each node under serial simulation by serial simulation;
the obtaining of the parallel simulation result comprises: the method comprises the steps of time required by each partition for performing parallel simulation analysis and calculation, overhead time such as communication and storage required by parallel simulation calculation and analysis, cost for communication and synchronization in parallel calculation, total time required by parallel calculation, and voltage results of each node under parallel simulation.
In the step 2, based on the influence of the power distribution network simulation partitioning scheme on the parallel analysis and calculation of the power distribution network, the influence factors are used as partitioning evaluation indexes;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation quantity of all partition nodes, dividing degree of feeder nodes, dividing degree of substation nodes, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
The simulation calculated amount balance degree of the parallel simulation partition nodes is as follows: the method is used for evaluating whether the calculated amount of each node in the power distribution network simulation partitioning scheme is balanced or not;
because different partitions can contain the PQ node, the PV node and the ZIP node, the simulation calculation amount of the three types of load nodes is different, and the simulation calculation amount of each partition is the sum of the simulation calculation amount of all the nodes in the partition;
in the embodiment, simulation calculation amount of each PQ node is defined as 1 according to simulation time through simulink software modeling simulation;
defining simulation calculated amount of each PV node according to simulation time, converting the simulation calculated amount into equivalent PQ node load number, and determining according to the following formula:
n cal,PV ＝α PV ×n cal,PQ (1)
in the formula, n cal,PV Converting PV node into equivalent PQ node load number, alpha PV As PV-PQ conversion factor, n cal,PQ The PQ node load number;
defining simulation calculated quantity of each ZIP node according to simulation time, converting the simulation calculated quantity into equivalent PQ node load quantity, and determining according to the following formula:
n cal,ZIP ＝α ZIP ×n cal,PQ (2)
in the formula, n cal,ZIP Conversion to an equivalent PQ node load number, α, for ZIP nodes ZIP Is the ZIP-PQ conversion coefficient;
the ZIP load model based on the voltage static characteristics is determined according to the following formula:
wherein P (V), Q (V) are the active and reactive power of the load model with ZIP as a whole, V, V 0 Respectively an actual voltage and a reference voltage, P 0 、Q 0 Active power and reactive power when the voltage is reference voltage respectively, and alpha (alpha '), beta (beta ') and gamma (gamma ') are active (reactive) Z, I, P load proportionality coefficients respectively;
wherein α + β + γ ═ 1, α ' + β ' + γ ' ═ 1;
the equivalent PQ node load number calculated by the ith partition can be obtained from the formulas (1), (2) and (3), and is determined according to the following formula:
n cal,i ＝N PQ,i ×n cal,PQ +N PV,i ×n cal,PV +N ZIP,i ×n cal,ZIP (4)
in the formula, n cal,i The equivalent PQ node load number, N, calculated by the ith partition PQ,i The number of PQ nodes in the ith partition, N PV,i The number of PV nodes in the ith partition, N ZIP,i The number of ZIP nodes in the ith partition.
Therefore, the simulation calculation amount balance degree of the parallel simulation partition nodes is determined according to the following formula:
in the formula, y 1 Simulating the balance of calculated quantity for the parallel simulation partition nodes, wherein i is the number of partitions;
when the simulation calculation amount of the parallel simulation partition node calculated in the formula (5) is closer to 1, the balance degree of the simulation calculation amount of the parallel simulation partition node is better.
Degree of division of feeder nodes: the method is used for judging whether nodes on the same feeder are divided into the same partition as much as possible so as to evaluate the influence of the division of the feeder nodes on the parallel computing analysis in the partition scheme;
the division degree of the feeder line nodes is determined according to the following formula:
y 2 ＝L i,end -L i+1,st (6)
in the formula, y 2 For degree of division of feeder node, L i,end Is the feeder number, L, to which the last node of the ith partition calculation belongs i+1,st The feeder number of the first node calculated by the (i + 1) th partition is the feeder number of the first node;
when the division degree of the feeder node calculated in the formula (6) is 1, the division degree of the feeder is the best, and the division degree also indicates that the nodes on the same feeder in the power distribution network simulation partition scheme are divided into the same partition as much as possible.
The division degree of the transformer substation nodes is as follows: the method is used for judging whether nodes of the same substation in the power distribution network simulation partition scheme are partitioned in the same partition as much as possible so as to evaluate the influence of the node partitioning of the substation on parallel computing analysis in the partition scheme;
the division degree of the transformer substation nodes is determined according to the following formula:
y 3 ＝M i,end -M i+1,st (7)
in the formula, y 3 Degree of division for substation nodes, M i,end Is the substation number to which the last node of the ith zone calculation belongs, M i+1,st The number of the substation to which the first node of the (i + 1) th zone calculation belongs;
when the division degree result of the substation nodes obtained by the formula (7) is 1, the best result is obtained, which indicates that the nodes in the same substation in the power distribution network simulation zoning scheme are divided into the same zones as much as possible.
Parallel computing acceleration ratio: a measure of the fraction of time that the server is effectively utilized, defined as the ratio of the serial system runtime to the parallel system runtime, is determined as follows:
in the formula, T s For serial system running time, T p Running time for the parallel system;
however, since the parallel time obtained after the partition modeling operation is now independent for each partition (for example, 10 partitions are divided to obtain 10 partition calculation times), and the synchronous communication time of the partitions is not included, the right-side fractional numerator is used to estimate the operation time of the partition task in a multi-core processor, and when calculating the weight, the optimal value of each index is expected to be 1, so the ratio of the acceleration ratio to the number of partitions is defined, therefore, the acceleration ratio is calculated by the parallel simulation according to the following formula:
in the formula, y 4 For parallel calculation of acceleration ratio, T ∞ Time required for serial simulation analysis calculation, T i Parallel emulation for ith partitionTime required for true analysis calculation, T t The overhead time for communication, storage and the like required during parallel computing and analysis, wherein i is the number of partitions;
wherein, when the parallel computation acceleration ratio obtained by the equation (9) is equal to 1, the acceleration of the parallel computation is the highest.
Parallel simulation calculation efficiency: the measure of the reduction degree of the operation time of the server after parallel optimization is originally defined as the ratio of the acceleration ratio to the number of calculation cores, and is determined according to the following formula
In the formula, e is the original parallel simulation calculation efficiency, s is the original acceleration ratio, and n is the original calculation kernel number;
but due to the uncertainty of the parallel time, the computational efficiency of the parallel simulation is determined according to the following formula:
in the formula, y 5 Calculating efficiency for parallel simulation;
wherein, the parallel overhead is caused by storage, communication and other factors, and in the formula (11), the time T required by the serial simulation analysis calculation ∞ Determined as follows:
T ∞ ＜i×(max(T 1 ,T 2 ...T i-1 ,T i )+T t )
when the result of the parallel simulation calculation efficiency calculated by the equation (11) is close to 1, the efficiency of the parallel simulation calculation is optimal.
The expansibility of the power distribution network in parallel simulation: representing a distribution network in a parallel simulation such a distribution network is called a scalable system if it can be scaled (scaled down) to a larger (smaller) distribution network and linearly increase (decrease) its performance (cost). According to the definition, the expansion of the system to improve the system performance or the reduction of the system to improve the system cost performance belongs to the expandability category;
the expansibility of the power distribution network in parallel simulation is determined according to the following formula:
in the formula, y 6 For the expansibility of a distribution network in parallel simulation, T 0 For communication and synchronisation costs in parallel operation of the system, T 0 ' total time required for parallel operation;
where the total time T required for the parallel operation 0 ' maximum time max (T) by parallel calculation after partitioning 1 ,T 2 ...T i-1 ,T i ) Time T for communication and synchronization in parallel operation with system t Is expressed by the sum of the two, so that the total time T required for parallel operation 0 ', determined as follows:
the scalability of the power distribution network reflected in equation (12) is determined by the parallel overhead and the relative time of the parallel computation. When the parallel overhead is smaller than the computing time, the expandability of the parallel computing power distribution network is stronger. The expansibility measurement index can more intuitively reflect the expandability of the power distribution network, and the expansibility index is close to 1 in a system with good expansibility.
Parallel computing precision: the device is used for judging the calculation precision of the evaluated subarea under the parallel simulation calculation;
in this embodiment, in order to compare the simulation accuracy of parallel and serial computations, the same simulation platform is used for verification, and both the simulation model and the solving method have no difference. In the serial simulation mode, the calculation accuracy is highest, so that the actual value of the calculation result can be used as the real value of the system operation, the simulation result difference value of the parallel calculation and the serial calculation is defined as the simulation error, the ratio of the simulation error average value to the measured measurement average value is defined as the simulation error percentage, and the voltage error average percentage is determined according to the following formula:
in the formula, eta is the average percentage of voltage errors between the serial simulation and the parallel simulation, n is the number of nodes, and U par,m 、U ser,m Respectively carrying out parallel and serial simulation voltage results on the mth node;
the parallel computing precision index can be obtained from the formula (13) and is determined according to the following formula:
y 7 ＝1-η (14)。
in the formula, y 7 The accuracy is calculated in parallel.
In the step 3, based on each evaluation index value of the power distribution network partition scheme calculated in the step 2, determining index weight corresponding to each evaluation index by using an analytic hierarchy process;
the Analytic Hierarchy Process (AHP) is a practical multi-factor decision method proposed in 1971. The method decomposes some complex fuzzy multi-target decision problems into a plurality of targets or criteria, further decomposes the problems into a plurality of layers of multi-index, then carries out fuzzy quantization on qualitative indexes, calculates the single-layer ordering (weight) and the total ordering of the layers and obtains the final evaluation result. The analytic hierarchy process is suitable for the problems which are difficult to completely quantify and analyze, and the method is flexible and easy to understand. The method is a qualitative and quantitative combined decision-making method, and is suitable for the decision-making problem that the expert qualitatively judges the important function of the method and the decision-making result is difficult to directly and accurately measure.
Firstly, a judgment matrix for determining the index weight is constructed, namely, each index of the weight to be determined is simply compared, judged and calculated, so that the corresponding weight of each index is obtained, and a decision basis is provided for the evaluation of the decomposition scheme of the power distribution network. Firstly, the importance degrees of the evaluation indexes of the power distribution network decomposition scheme are compared pairwise, and a questionnaire is filled in by experts, as shown in table 1.
TABLE 1
| | Index | 1 | Index 2 | … | Index i | … | Index m |
| Index 1 | b 11 | b 12 | … | b 1j | … | b 1m | |
| Index 2 | b 21 | b 22 | … | b 2j | … | b 2m | |
| … | … | … | … | … | … | … | |
| Index i | b i1 | b i2 | … | b ij | … | b im | |
| … | … | … | … | … | … | … | |
| Index m | b m1 | b m2 | … | b mj | … | b mm |
In Table 1, bij reflects the degree of importance of index i relative to index j;
in this embodiment, a reciprocal 1-9 scaling method is used, and the two-by-two comparison results are written into numerical values to form a judgment matrix, and the scaling meaning is shown in table 2.
TABLE 2
A decision matrix is established according to table 1 and determined as follows:
wherein B is a judgment matrix, B mm The importance degree of the mth evaluation index to the mth evaluation index is set;
in this embodiment, when the expert scores the importance of the evaluation index, the score is difficult to be completely consistent due to strong subjectivity. In order to ensure the reasonability and the correctness of the comment result, the consistency check is carried out on the judgment matrix, and the matrix which does not pass the consistency check is corrected and then checked.
The random consistency ratio is determined as follows:
in the formula, CR is a random consistency ratio, CI is a consistency index, and RI is an average random consistency index;
m evaluation indexes for comparison are set, and for any judgment matrix: maximum feature root λ when the matrices are completely consistent max M; when the matrix does not remain exactly the same, λ max More than m; then, there is a consistency index CI, which is determined as follows:
in the formula, m is the order of the judgment matrix;
wherein, the average random consistency index RI is shown in Table 3.
TABLE 3
And when CR is less than 0.1, the judgment matrix is considered to meet the consistency requirement, otherwise, the value of the element in the judgment matrix needs to be adjusted, and the consistency check is carried out again.
Calculating the weight of the index
Firstly, solving the maximum eigenvalue of the judgment matrix B and the corresponding eigenvector thereof, and calculating according to the following formula:
Bv＝λ max v (18)
in the formula, λ max V is the maximum eigenvalue of the judgment matrix B, and v is the maximum eigenvalue eigenvector of the judgment matrix B;
then, the weight of each evaluation index is obtained according to the feature vector v, and the calculation is performed according to the following formula
In the formula, ω a The index weight of the a-th index, v a A characteristic vector of the a index; v. of k Is the feature vector of the k index, and m is the evaluation index.
After the weight of each evaluation index is obtained, determining an evaluation result calculation formula of the power distribution network simulation partitioning scheme, carrying out weighted summation on the extracted index values, and evaluating the power distribution network partitioning scheme according to the obtained result;
the evaluation result calculation formula of the power distribution network parallel simulation partition scheme is determined according to the following formula:
Y＝ω 1 y 1 +ω 2 y 2 +ω 3 y 3 +ω 4 y 4 +ω 5 y 5 +ω 6 y 6 +ω 7 y 7 (20)
in the formula, ω 1 Simulating weight values, omega, of the calculated quantity balance for the nodes in each partition 2 Dividing the weight value of degree for the feeder node, omega 3 Dividing the weight value of degree, omega, for the transformer substation feeder node 4 Computing weight values, omega, of acceleration ratios for parallel simulations 5 Weight value of parallel simulation calculation efficiency, omega 6 Power distribution network expansibility weight value omega in parallel simulation 7 Calculating a weight value of the precision for the parallel simulation;
for different network partition schemes, the evaluation scores corresponding to the simulation partition scheme can be obtained through the formula (20), the partition scheme is divided into four grades of high, good, medium and poor according to the advantages and the disadvantages according to different evaluation scores, and the specific division standard is shown in table 4.
TABLE 4
| Rating of evaluation | Youyou (an instant noodle) | Good wine | In | Difference (D) |
| Fractional range | 0.7～1 | 0.55～0.7 | 0.4～0.55 | 0～0.4 |
The invention provides an evaluation method applied to a power distribution network simulation partition scheme, wherein power distribution network simulation partition evaluation indexes are determined based on influence factors of power distribution network simulation partitions on parallel simulation calculation, and the evaluation indexes comprise the balance degree of simulation calculation quantity of different partition nodes, the division degree of feeder nodes, the division degree of substation nodes, the acceleration ratio of parallel simulation calculation, the parallel simulation calculation efficiency, the expansibility of a power distribution network under parallel simulation and the parallel simulation calculation precision, the advantages and the disadvantages of the parallel calculation efficiency after partitioning and the height of the calculation precision can be reflected more accurately through the evaluation indexes, and the reliability of an evaluation result is higher; and the evaluation of the power distribution network simulation partition scheme is carried out by adopting an analytic hierarchy process, and the input mode of pairwise comparison and judgment is applied, so that the factors which are not effectively quantified can be easily processed, and the method is flexible and easy to understand.
Example 2
Based on the same inventive concept, the present invention further provides an evaluation system for a partitioning scheme applied to power distribution network simulation, as shown in fig. 3, including:
the acquisition module is used for acquiring the structural parameters of the power distribution network simulation partitioning scheme, the parallel simulation result after the power distribution network is partitioned and the serial simulation result of the power distribution network;
the calculation module is used for calculating each evaluation index value for evaluating the power distribution network partition scheme by utilizing the structural parameters, the serial simulation result and the parallel simulation result;
the evaluation module is used for calculating the scores of the power distribution network partition schemes based on the evaluation index values and the weights corresponding to the evaluation indexes predetermined by utilizing an analytic hierarchy process, and obtaining the quality grades of the power distribution network partition schemes based on the scores and the preset partition standards;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
In this embodiment, the distribution network includes two or more distribution areas and their subsidiary lines; the parallel simulation analysis refers to that when the power distribution network is subjected to simulation analysis, the power distribution network is divided into a plurality of areas, and meanwhile simulation analysis calculation is carried out on the plurality of areas.
The obtaining module is configured to obtain a structural parameter of the power distribution network simulation partitioning scheme, and includes: the number of PQ nodes, the number of PQ node loads, the number of PV nodes, the number of loads of PV node equivalent PQ nodes, the number of ZIP nodes, the number of loads of ZIP node equivalent PQ nodes, the number of partitions, the node numbers of all feeders and the node numbers of all substations in the power distribution network simulation partitioning scheme;
acquiring the serial simulation result comprises: analyzing and calculating the required time and the voltage result of each node under serial simulation by serial simulation;
the obtaining of the parallel simulation result comprises: the method comprises the steps of calculating time required by each partition for performing parallel simulation analysis, overhead time such as communication and storage required by parallel simulation calculation analysis, cost for communication and synchronization in parallel calculation, total time required by parallel calculation and voltage results of each node under parallel simulation.
The calculation module is used for calculating the balance degree of the node simulation calculation amount in all the subareas, the division degree of the nodes of the feeder line and the transformer substation, the parallel simulation calculation acceleration ratio, the parallel simulation calculation efficiency, the expansibility of the power distribution network in parallel simulation and the parallel simulation calculation precision according to the parameters acquired by the acquisition module;
the simulation calculation amount balance degree of the parallel simulation partition nodes is as follows: the method is used for evaluating whether the calculated amount of each node in the power distribution network simulation partitioning scheme is balanced or not;
because different partitions can contain the PQ node, the PV node and the ZIP node, the simulation calculation amount of the three types of load nodes is different, and the simulation calculation amount of each partition is the sum of the simulation calculation amount of all the nodes in the partition;
in this embodimentModeling simulation through simulink software, and defining simulation calculated quantity n of each PQ node according to simulation time cal,PQ Is 1;
defining simulation calculated quantity n of each PV node according to simulation time cal,PV And converting the equivalent PQ node load number into an equivalent PQ node load number, and determining the equivalent PQ node load number according to the following formula:
n cal,PV ＝α PV ×n cal,PQ (1)
in the formula, n cal,PV Converting PV node into equivalent PQ node load number, alpha PV For the PV-PQ conversion coefficient, n cal,PQ The PQ node load number;
defining simulation calculated quantity n of each ZIP node according to simulation time cal,ZIP And converting the equivalent PQ node load number into an equivalent PQ node load number, and determining the equivalent PQ node load number according to the following formula:
n cal,ZIP ＝α ZIP ×n cal,PQ (2)
in the formula, n cal,ZIP Conversion to an equivalent PQ node load number, α, for ZIP nodes ZIP Is the ZIP-PQ conversion coefficient;
the ZIP load model based on the voltage static characteristics is determined according to the following formula:
wherein P (V), Q (V) are the active and reactive power of the load model with ZIP as a whole, V, V 0 Respectively an actual voltage and a reference voltage, P 0 、Q 0 Respectively active power and reactive power when the voltage is reference voltage, and respectively alpha (alpha '), beta (beta '), and gamma (gamma ') are active (reactive) Z, I, P load proportionality coefficients;
wherein α + β + γ ═ 1, α ' + β ' + γ ' ═ 1;
the equivalent PQ node load number calculated by the ith partition can be obtained from the formulas (1), (2) and (3), and is determined according to the following formula:
n cal,i ＝N PQ,i ×n cal,PQ +N PV,i ×n cal,PV +N ZIP,i ×n cal,ZIP (4)
in the formula, n cal,i The equivalent PQ node load number, N, calculated by the ith partition PQ,i The number of PQ nodes in the ith partition, N PV,i The number of PV nodes in the ith partition, N ZIP,i The number of ZIP nodes in the ith partition.
Therefore, the simulation calculation amount balance degree of the parallel simulation partition nodes is determined according to the following formula:
in the formula, y 1 Simulating the balance of calculated quantity for the parallel simulation partition nodes, wherein i is the number of partitions;
when the simulation calculation amount of the parallel simulation partition node calculated in the formula (5) is closer to 1, the balance degree of the simulation calculation amount of the parallel simulation partition node is better.
The division degree of the feeder nodes is as follows: the method is used for judging whether nodes on the same feeder are divided into the same partition as much as possible so as to evaluate the influence of the division of the feeder nodes in the partition on the parallel computing analysis;
the division degree of the feeder line nodes is determined according to the following formula:
y 2 ＝L i,end -L i+1,st (6)
in the formula, y 2 Degree of division for feeder nodes, L i,end Is the feeder number, L, to which the last node of the ith partition calculation belongs i+1,st The feeder number of the first node calculated by the (i + 1) th partition is the feeder number of the first node;
when the division degree of the feeder node calculated in the formula (6) is 1, the division degree of the feeder is the best, and the nodes on the same feeder are also shown to be divided into the same partition as much as possible.
Dividing degree of the transformer substation nodes: the method is used for judging whether the nodes of the same transformer substation are divided into the same subarea as much as possible or not so as to evaluate the influence of the division of the nodes of the transformer substation in the subareas on the parallel computing analysis;
the division degree of the transformer substation nodes is determined according to the following formula:
y 3 ＝M i,end -M i+1,st (7)
in the formula, y 3 Degree of division for substation nodes, M i,end Is the substation number to which the last node of the ith partition calculation belongs, M i+1,st The number of the substation to which the first node calculated by the (i + 1) th partition belongs;
when the division degree result of the substation nodes obtained by the formula (7) is 1, the node division degree result is the best, and the nodes in the same substation are divided into the same subarea as much as possible.
Parallel computing acceleration ratio: a measure of the fraction of time that the server is effectively utilized, defined as the ratio of the serial system runtime to the parallel system runtime, is determined as follows:
in the formula, T s For serial system running time, T p Running time for the parallel system;
however, since the parallel time obtained after the partition modeling operation is now independent for each partition (for example, 10 partitions are divided to obtain 10 partition calculation times), and does not include the synchronous communication time of the partitions, the right-side fractional numerator is used to estimate the operation time of the partition task in a multi-core processor, and when calculating the weight, the optimal value of each index is expected to be 1, so the ratio of the acceleration ratio to the number of partitions is defined, and therefore, the acceleration ratio is calculated by the parallel simulation according to the following formula:
in the formula, y 4 For parallel calculation of acceleration ratio, T ∞ Time required for serial simulation analysis calculation, T i Time required for parallel simulation analysis calculation for ith partition, T t Overhead time for communication, storage and the like required during parallel computing and analysis;
wherein, when the parallel computation acceleration ratio obtained by the equation (9) is equal to 1, the acceleration of the parallel computation is the highest.
Parallel simulation calculation efficiency: the method is characterized in that the measure of the reduction degree of the operation time of the server after parallel optimization is originally defined as the ratio of the acceleration ratio to the number of calculation cores, and is determined according to the following formula
In the formula, e is the original parallel simulation calculation efficiency, s is the original acceleration ratio, and n is the original calculation kernel number;
but due to the uncertainty of the parallel time, the computational efficiency of the parallel simulation is determined according to the following formula:
in the formula, y 5 Calculating efficiency for parallel simulation;
wherein, the storage, communication and other factors may cause the parallel overhead, and in the formula (11), the T ∞ Determined as follows:
T ∞ ＜i×(max(T 1 ,T 2 ...T i-1 ,T i )+T t )
when the result of the parallel simulation calculation efficiency calculated by the equation (11) is close to 1, the efficiency of the parallel simulation calculation is optimal.
The expansibility of the power distribution network in parallel simulation: a distribution network is referred to as a scalable system if it can be scaled (scaled down) to a larger (smaller) distribution network and linearly increase (decrease) its performance (cost) in parallel simulations. According to the definition, the expansion of the system to improve the system performance or the reduction of the system to improve the system cost performance belongs to the expandability category;
the expansibility of the power distribution network in parallel simulation is determined according to the following formula:
in the formula, y 6 For the expansibility of a distribution network in parallel simulation, T 0 For communication and synchronisation costs in parallel operation of the system, T 0 ' total time required for parallel operation;
the total time T required for the parallel operation 0 ' maximum time max (T) by parallel calculation after partitioning 1 ,T 2 ...T i-1 ,T i ) Time T for communication and synchronization in parallel operation with system t Is expressed by the sum of the two, so that the total time T required for parallel operation 0 ', determined as follows:
the scalability of the power distribution network reflected in equation (12) is determined by the parallel overhead and the relative time of the parallel computation. When the parallel overhead is smaller than the computing time, the expandability of the parallel computing power distribution network is stronger. The expansibility measurement index can more intuitively reflect the expandability of the power distribution network, and the expansibility index is close to 1 in a system with good expansibility.
Parallel computing precision: the device is used for judging the calculation precision of the evaluated subarea under the parallel simulation calculation;
in this embodiment, in order to compare the simulation accuracy of parallel and serial computations, the same simulation platform is used for verification, and both the simulation model and the solving method have no difference. In the serial simulation mode, the calculation accuracy is highest, so that the actual value of the calculation result can be used as the real value of the system operation, the simulation result difference value of the parallel calculation and the serial calculation is defined as the simulation error, the ratio of the simulation error average value to the measured measurement average value is defined as the simulation error percentage, and the voltage error average percentage is determined according to the following formula:
in the formula, eta is the average percentage of voltage errors between the serial simulation and the parallel simulation, n is the number of nodes, and U par,m 、U ser,m Respectively carrying out parallel and serial simulation voltage results on the mth node;
the parallel computing precision index can be obtained from the formula (13) and is determined according to the following formula:
y 7 ＝1-η (14)
in the formula, y 7 The accuracy is calculated in parallel.
The evaluation module comprises: a weight determination submodule and an evaluation result determination submodule;
the weight determination submodule is used for determining index weights corresponding to the evaluation indexes by utilizing an analytic hierarchy process based on the evaluation index values of the power distribution network partition scheme calculated in the step 2;
the Analytic Hierarchy Process (AHP) is a practical multi-factor decision method proposed in 1971. The method decomposes some complex fuzzy multi-target decision problems into a plurality of targets or criteria, further decomposes the problems into a plurality of layers of multi-index, then carries out fuzzy quantization on qualitative indexes, calculates the single-layer ordering (weight) and the total ordering of the layers and obtains the final evaluation result. The analytic hierarchy process is suitable for the problems which are difficult to completely quantify and is flexible and easy to understand. The method is a qualitative and quantitative combined decision-making method, and is suitable for the decision-making problem that experts judge the important function qualitatively and are difficult to measure the decision-making result directly and accurately.
Firstly, a judgment matrix for determining the index weight is constructed, namely, each index of the weight to be determined is simply compared, judged and calculated, so that the corresponding weight of each index is obtained, and a decision basis is provided for the evaluation of the decomposition scheme of the power distribution network. Firstly, the importance degrees of the evaluation indexes of the power distribution network decomposition scheme are compared pairwise, and a questionnaire is filled in by experts, as shown in table 1.
TABLE 1
| | Index | 1 | Index 2 | … | Index i | … | Index m |
| Index 1 | b 11 | b 12 | … | b 1j | … | b 1m | |
| Index 2 | b 21 | b 22 | … | b 2j | … | b 2m | |
| … | … | … | … | … | … | … | |
| Index i | b i1 | b i2 | … | b ij | … | b im | |
| … | … | … | … | … | … | … | |
| Index m | b m1 | b m2 | … | b mj | … | b mm |
In Table 1, bij reflects the degree of importance of index i relative to index j;
in this embodiment, a reciprocal 1-9 scaling method is used, and the two-by-two comparison results are written into numerical values to form a judgment matrix, and the scaling meaning is shown in table 2.
TABLE 2
A decision matrix B is established according to table 1, determined as follows:
where B is a judgment matrix, B mm The importance degree of the mth evaluation index to the mth evaluation index is set;
in this embodiment, when the expert scores the importance of the evaluation index, the score is difficult to be completely consistent due to strong subjectivity. In order to ensure the reasonability and the correctness of the comment result, the consistency check is carried out on the judgment matrix, and the matrix which does not pass the consistency check is corrected and then checked.
The random consistency ratio is determined as follows:
in the formula, CR is a random consistency ratio, CI is a consistency index, and RI is an average random consistency index;
m evaluation indexes for comparison are set, and for any judgment matrix: maximum feature root λ when the matrix remains completely uniform max M; when the matrix does not remain exactly the same, λ max More than m; then, there is a consistency index CI, which is determined as follows:
in the formula, m is the order of the judgment matrix;
wherein the average random consistency index RI is shown in table 3.
TABLE 3
And when CR is less than 0.1, the judgment matrix is considered to meet the consistency requirement, otherwise, the value of the element in the judgment matrix needs to be adjusted, and the consistency check is carried out again.
When calculating the weight index, firstly, the maximum eigenvalue of the judgment matrix B and the corresponding eigenvector are solved, and the calculation is carried out according to the following formula:
Bv＝λ max v (18)
in the formula, λ max V is the maximum eigenvalue eigenvector of the judgment matrix B;
secondly, the weight of each evaluation index is obtained according to the characteristic vector v, and the weight is calculated according to the following formula:
in the formula, ω a The index weight of the a-th index, v a A characteristic vector of the a index; v. of k Is the feature vector of the k index, and m is the evaluation index.
The evaluation result determination submodule is used for determining an evaluation result calculation formula of the power distribution network simulation partitioning scheme after the weight of each evaluation index is utilized, carrying out weighted summation on the provided index values, and determining the evaluation result of the power distribution network partitioning scheme according to the obtained scores;
the evaluation result calculation formula of the power distribution network parallel simulation partition scheme is determined according to the following formula:
Y＝ω 1 y 1 +ω 2 y 2 +ω 3 y 3 +ω 4 y 4 +ω 5 y 5 +ω 6 y 6 +ω 7 y 7 (20)
in the formula，ω 1 Weight value omega of the balance of the simulation calculated quantity of the nodes in each subarea 2 Dividing the weight value of degree for the feeder node, omega 3 Weight value of degree, omega, for substation feeder node division 4 Computing weight values, omega, of acceleration ratios for parallel simulations 5 Weight value of parallel simulation computing efficiency, omega 6 Power distribution network expansibility weight value omega in parallel simulation 7 Calculating a weight value of the precision for the parallel simulation;
for different network partition schemes, the evaluation scores corresponding to the simulation partition scheme can be obtained through the formula (20), the partition scheme is divided into four grades of high, good, medium and poor according to the advantages and the disadvantages according to different evaluation scores, and the specific division standard is shown in table 4.
TABLE 4
| Rating of evaluation | Superior food | Good quality | In | Difference (D) |
| Fractional range | 0.7～1 | 0.55～0.7 | 0.4～0.55 | 0～0.4 |
It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The present invention is not limited to the above embodiments, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention are included in the scope of the claims of the present invention which are filed as the application.
Claims (11)
1. A partition scheme evaluation method applied to power distribution network simulation is characterized by comprising the following steps:
acquiring structural parameters of a power distribution network simulation partitioning scheme, a parallel simulation result after partitioning of the power distribution network and a serial simulation result of the power distribution network;
calculating each evaluation index value for evaluating the power distribution network partition scheme by using the structural parameters, the serial simulation result and the parallel simulation result;
calculating scores of the power distribution network partition schemes based on the evaluation index values and weights corresponding to the evaluation indexes predetermined by an analytic hierarchy process, and obtaining quality grades of the power distribution network partition schemes based on the scores and preset partition standards;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
2. The method according to claim 1, wherein the calculating an evaluation result of the power distribution network simulation partitioning scheme based on the evaluation index values and the weights corresponding to the evaluation indexes acquired in advance comprises:
summing the products of the evaluation index values obtained by calculation and the weights corresponding to the evaluation index values to obtain the evaluation scores of the power distribution network simulation partition scheme;
and classifying according to preset grades based on the evaluation scores of the power distribution network simulation partition scheme to obtain the evaluation result of the power distribution network simulation partition scheme.
3. The method of claim 2, wherein the evaluation score of the power distribution network simulation partition scheme is calculated according to the following formula:
Y＝ω 1 y 1 +ω 2 y 2 +ω 3 y 3 +ω 4 y 4 +ω 5 y 5 +ω 6 y 6 +ω 7 y 7
in the formula, Y is the evaluation score of the power distribution network simulation partition scheme, omega 1 Simulating weight values, y, of the computation load balance for the nodes in each partition 1 Simulating the computation of the metric balance, ω, for the nodes in each partition 2 Dividing the weight value of degree for the feeder node, y 2 Degree of division, omega, for feeder nodes 3 Weight value y for division degree of substation feeder node 3 Degree of division, omega, for substation feeder nodes 4 Computing weight values, y, of acceleration ratios for parallel simulations 4 Computing acceleration ratio, omega, for parallel simulation 5 Weight value of parallel simulation computing efficiency, y 5 For parallel simulation computational efficiency, omega 6 Weight value of expansibility, y, of power distribution network in parallel simulation 6 For the expansibility, omega, of distribution networks in parallel simulation 7 Computing weight values of precision, y, for parallel simulations 7 The accuracy is calculated for the parallel simulation.
4. The method of claim 3, wherein the configuration parameters of the power distribution network simulation partitioning scheme comprise:
the number of PQ nodes, the number of PQ node loads, the number of PV nodes, the number of loads of PV node equivalent PQ nodes, the number of ZIP nodes, the number of loads of ZIP node equivalent PQ nodes, the number of partitions, the node numbers of all feeders and the node numbers of all substations in the power distribution network simulation partitioning scheme;
the serial simulation result comprises the following steps: analyzing and calculating the required time and the voltage result of each node under serial simulation by serial simulation;
the parallel simulation result comprises: the method comprises the steps of calculating time required by each partition for performing parallel simulation analysis, overhead time such as communication and storage required by parallel simulation calculation analysis, cost for communication and synchronization in parallel calculation, total time required by parallel calculation and voltage results of each node under parallel simulation.
5. The method of claim 4, wherein the balance y of the simulation computation of all the partition nodes is 1 Calculated as follows:
in the formula, n cal,i The equivalent PQ node load number of the ith partition;
wherein the equivalent PQ node load number n of the ith partition cal,i Determined as follows:
n cal,i ＝N PQ,i ×n cal,PQ +N PV,i ×n cal,PV +N ZIP,i ×n cal,ZIP
in the formula, N PQ,i The number of PQ nodes in the ith partition, n cal,PQ Is the PQ node load number, N PV,i The number of PV nodes in the ith partition, n cal,PV Number of loads for PV node equivalent PQ node, N ZIP,i Is the number of ZIP nodes in the ith partition, n cal,ZIP The number of loads of the PQ node equivalent to the ZIP node is represented by i, and the number of partitions is represented by i.
6. The method of claim 5, wherein the parallel simulation calculates an acceleration ratio y 4 Calculated as follows:
in the formula, T ∞ Time required for serial simulation analysis calculation, T i Time required for parallel simulation analysis calculation for ith partition, T t Overhead time such as communication and storage required for parallel simulation calculation and analysis;
preferably, the parallel simulation has a computational efficiency y 5 Determined as follows:
7. the method of claim 6, wherein the power distribution network has a scalability y in parallel simulation 6 Calculated as follows:
in the formula, T 0 Cost for communication and synchronization in parallel computing, T 0 ' total time required for parallel computation;
wherein the total time T required for the parallel computation 0 ', calculated as:
8. the method of claim 4, wherein the parallel simulation has a computational accuracy y 7 Calculated as follows:
y 7 ＝1-η
in the formula, eta is the average percentage of voltage errors between the serial simulation and the parallel simulation;
wherein the average percentage η of voltage errors between the serial simulation and the parallel simulation is calculated as follows:
in the formula, n is the number of nodes, U par,m 、U ser,m The voltage results of the mth node under parallel and serial simulations, respectively.
9. The method according to claim 2, wherein the determining of the weight corresponding to each evaluation index includes:
s1: evaluating the importance of each evaluation index by adopting a reciprocity scaling method in an analytic hierarchy process, and forming a judgment matrix;
s2: carrying out consistency check on the judgment matrix, and entering S3 if the consistency check requirement is met, or entering S1 to adjust the value of an element in the judgment matrix;
s3: and calculating the index weight of each evaluation index based on the judgment matrix passing the consistency test.
10. The method according to claim 9, wherein the index weight under the evaluation index system is calculated as follows:
in the formula, ω a The index weight of the a-th index, v a A characteristic vector of the a index; v. of k Is the feature vector of the k index, and m is the evaluation index.
11. The evaluation system of the partition scheme applied to the power distribution network simulation is characterized by comprising the following steps of:
the acquisition module is used for acquiring the structural parameters of the power distribution network simulation partitioning scheme, the parallel simulation result after the power distribution network is partitioned and the serial simulation result of the power distribution network;
the calculation module is used for calculating each evaluation index value for evaluating the power distribution network partition scheme by utilizing the structural parameters, the serial simulation result and the parallel simulation result;
the evaluation module is used for calculating the scores of the power distribution network partition schemes based on the evaluation index values and the weights corresponding to the evaluation indexes predetermined by utilizing an analytic hierarchy process, and obtaining the quality grades of the power distribution network partition schemes based on the scores and the preset partition standards;
wherein the evaluation index includes: the method comprises the steps of balancing simulation calculation amount of nodes in all the partitions, dividing degrees of the nodes to which the feeder lines and the transformer substations belong, accelerating ratio of parallel simulation calculation, parallel simulation calculation efficiency, expansibility of a power distribution network in parallel simulation and parallel simulation calculation precision.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110219443.XA CN114970072A (en) | 2021-02-26 | 2021-02-26 | Partition scheme evaluation method and system applied to power distribution network simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110219443.XA CN114970072A (en) | 2021-02-26 | 2021-02-26 | Partition scheme evaluation method and system applied to power distribution network simulation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114970072A true CN114970072A (en) | 2022-08-30 |
Family
ID=82973558
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110219443.XA Pending CN114970072A (en) | 2021-02-26 | 2021-02-26 | Partition scheme evaluation method and system applied to power distribution network simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114970072A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115827227A (en) * | 2022-11-28 | 2023-03-21 | 中国华能集团清洁能源技术研究院有限公司 | Serial-parallel resource optimal allocation method and system for simulation system |
-
2021
- 2021-02-26 CN CN202110219443.XA patent/CN114970072A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115827227A (en) * | 2022-11-28 | 2023-03-21 | 中国华能集团清洁能源技术研究院有限公司 | Serial-parallel resource optimal allocation method and system for simulation system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112566196B (en) | Heterogeneous network access selection method based on smart grid and related equipment | |
| CN110417011B (en) | Online dynamic security assessment method based on mutual information and iterative random forest | |
| US20170300546A1 (en) | Method and Apparatus for Data Processing in Data Modeling | |
| JP6784780B2 (en) | How to build a probabilistic model for large-scale renewable energy data | |
| CN108052387B (en) | Resource allocation prediction method and system in mobile cloud computing | |
| CN107360026B (en) | Distributed message middleware performance prediction and modeling method | |
| Li et al. | Production task queue optimization based on multi-attribute evaluation for complex product assembly workshop | |
| CN111882198A (en) | Project performance evaluation method and system | |
| CN108090615B (en) | Minimum frequency prediction method after power system fault based on cross entropy integrated learning | |
| CN107069696A (en) | A kind of parallel calculating method of Power system state estimation | |
| CN108336739A (en) | A kind of Probabilistic Load Flow on-line calculation method based on RBF neural | |
| CN111966574A (en) | Architecture scheme evaluation method of avionics general computing platform | |
| CN111967183A (en) | Method and system for calculating line loss of distribution network area | |
| CN109981749A (en) | A kind of cloud workflow task running time prediction method promoted based on limit gradient | |
| CN106612511A (en) | Wireless network throughput evaluation method and device based on support vector machine | |
| CN113505534A (en) | Load prediction method considering demand response | |
| CN114970072A (en) | Partition scheme evaluation method and system applied to power distribution network simulation | |
| CN109540089B (en) | Bridge deck elevation fitting method based on Bayes-Kriging model | |
| CN113328467B (en) | Probability voltage stability evaluation method, system, terminal device and medium | |
| CN114492742A (en) | Neural network structure searching method, model issuing method, electronic device, and storage medium | |
| CN117034851A (en) | Region decomposition method based on pruning and graph segmentation optimization | |
| CN108268982B (en) | Large-scale active power distribution network decomposition strategy evaluation method and device | |
| CN116523640A (en) | Financial information management system based on scheduling feedback algorithm | |
| CN113627655B (en) | Method and device for simulating and predicting pre-disaster fault scene of power distribution network | |
| CN114549233A (en) | Floating population prediction method based on combination of LGB algorithm and ARIMA algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |