CN114912767A - Three-phase active power distribution network power flow analysis method and device suitable for mass scenes - Google Patents
Three-phase active power distribution network power flow analysis method and device suitable for mass scenes Download PDFInfo
- Publication number
- CN114912767A CN114912767A CN202210415820.1A CN202210415820A CN114912767A CN 114912767 A CN114912767 A CN 114912767A CN 202210415820 A CN202210415820 A CN 202210415820A CN 114912767 A CN114912767 A CN 114912767A
- Authority
- CN
- China
- Prior art keywords
- power
- power flow
- distribution network
- determining
- multivariable
- 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
- 238000000034 method Methods 0.000 title claims abstract description 104
- 238000009826 distribution Methods 0.000 title claims abstract description 85
- 238000005206 flow analysis Methods 0.000 title claims abstract description 69
- 230000006870 function Effects 0.000 claims description 59
- 238000003860 storage Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000009795 derivation Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 41
- 239000000243 solution Substances 0.000 description 22
- 238000010586 diagram Methods 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 11
- 238000006467 substitution reaction Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000027311 M phase Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- 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
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Strategic Management (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The application provides a three-phase active power distribution network power flow analysis method and device suitable for massive scenes, and the method comprises the following steps: acquiring operation state data of a mass scene to be analyzed, wherein the operation state data comprises power proportion control factors corresponding to the mass scene to be analyzed; and inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed. The scheme can be used for rapidly carrying out load flow calculation on massive scenes.
Description
Technical Field
The invention belongs to the technical field of power system scheduling, and particularly relates to a three-phase active power distribution network power flow analysis method and device suitable for massive scenes.
Background
The generation of electricity by replacing the traditional fossil fuel with renewable energy is an important way to solve the problem of carbon emission. In order to ensure the cleanness and sustainable power supply of the power distribution network, a large number of distributed power supplies such as photovoltaic power, wind power and the like are connected into the power distribution network, and the traditional radial passive power distribution network is gradually changed into an active power distribution network. However, the high permeability and uncertainty of the distributed power supply increase the complexity of the power distribution network, and also lead to a drastic increase in the number of operating scenarios of the power distribution network, which brings great challenges to the analysis of the power distribution network. Whether power distribution network planning, operation or reliability analysis, massive operation scenes are comprehensively analyzed. Load flow calculations, as a fundamental analysis tool for modern power distribution network analysis and distribution management systems, should be done quickly and reliably.
Based on the characteristics of radial and weak mesh structures, high R/X ratio, multi-phase and unbalanced operation, more buses and branches and the like of the power distribution network, the currently proposed three-phase power distribution network power flow algorithm comprises a forward-push back substitution method, an improved Newton method, a direct method and an implicit Z method Bus The gaussian method, the linear power flow method, the current injection method, and the like. However, these methods are generally based on iterative methods, and the resulting trend solution is a numerical solution. And if a massive operation scene is to be analyzed, load flow calculation must be carried out repeatedly. In addition, another class of all-pure embedding methods based on a recursive method can represent a power flow solution as an all-pure function of a power proportional control factor, so as to obtain an analytic solution of the power flow. However, in a physical sense, since the method only embeds one power proportional control factor as a variable, the analyzed operation scene is required to be changed only along the same direction, and the spatial correlation and the difference among a large number of distributed power supplies are difficult to be described.
Therefore, if a fast power flow analysis is to be performed on a large number of arbitrary operation scenes, the above-mentioned problems of repeated calculation by an iterative method and a single direction of variable change by a recursive method must be solved.
Disclosure of Invention
The embodiment of the specification aims to provide a three-phase active power distribution network power flow analysis method and device suitable for massive scenes.
In order to solve the above technical problem, the embodiments of the present application are implemented as follows:
in a first aspect, the present application provides a three-phase active power distribution network power flow analysis method adapted to a massive scene, including:
acquiring operation state data of a mass scene to be analyzed, wherein the operation state data comprises power proportion control factors corresponding to the mass scene to be analyzed;
and inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
In one embodiment, the step of determining the offline power flow analysis expression includes:
dividing the system into at least two areas according to the type of a bus of the three-phase active power distribution network and the geographic position of the distributed power supply or the load, wherein each area corresponds to at least one independent power proportional control factor variable, and the power proportional control factor variables are respectively used for representing the power change characteristics of the distributed power supply or the load in the corresponding area;
based on a multidimensional pure embedding method, an offline power flow analytical expression with a plurality of power proportional control factor variables as multivariable is determined.
In one embodiment, the determining an offline power flow analytic expression with a plurality of power proportional control factor variables as multivariable based on a multidimensional full-pure embedding method includes:
according to a plurality of power proportional control factor variables, representing the voltage of each bus in the three-phase active power distribution network as a multivariable pure function with unknown coefficients based on a multidimensional pure embedding method;
determining an initial power balance equation of each bus according to the connection mode of the buses;
determining a fully-pure embedding equation according to a multivariable fully-pure function and an initial power balance equation;
determining unknown coefficients of the multivariable holomorphic functions according to the multivariable holomorphic functions and the holomorphic embedding equations;
and determining an offline power flow analytical expression according to the unknown coefficient of the multivariable holomorphic function.
In one embodiment, determining the unknown coefficients of the multivariate holohedral function according to the multivariate holohedral function and the holohedral embedding equation comprises:
substituting a multivariable pureness function into a pureness embedding equation to obtain a group of equation sets;
after the derivation transformation of the equation set, obtaining a recursion equation through the equality of coefficients with the same power on both sides of the equation;
and solving the unknown coefficient of the multivariable holo-pure function according to a recurrence equation.
In one embodiment, determining an offline power flow analytic expression according to unknown coefficients of a multivariate holohedral function includes:
determining a fractional expression by a multivariable pure function according to a generalized rational approximation technology;
determining the coefficient of the fractional expression according to the unknown coefficient of the multivariable holomorphic function;
and determining an offline power flow analysis expression according to the coefficient of the fractional expression.
In one embodiment, the bus bars are connected in a manner that includes a Y-connection or a delta-connection.
In a second aspect, the present application provides a three-phase active power distribution network power flow analysis device adapted to a massive scene, the device including:
the acquisition module is used for acquiring the running state data of the mass scene to be analyzed, and the running state data comprises a power proportion control factor corresponding to the mass scene to be analyzed;
and the processing module is used for inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
In a third aspect, the present application provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the method for analyzing a power flow of a three-phase active power distribution network adapted to a massive scene in the first aspect is implemented.
In a fourth aspect, the present application provides a readable storage medium, on which a computer program is stored, where the program, when executed by a processor, implements the method for power flow analysis of a three-phase active power distribution network adapted to mass scenes according to the first aspect.
As can be seen from the technical solutions provided in the embodiments of the present specification, the solution: the method can solve the problems of repeated calculation of the iterative load flow calculation method and single change direction of the power proportion control factor of the recursive load flow calculation method in the analysis of the existing mass scenes, and realizes the rapid load flow analysis of mass and arbitrary operation scenes.
The scheme can be suitable for massive scene power flow analysis of the large-scale three-phase active power distribution network, and has high calculation efficiency, so that decision support is provided for planning, probability power flow analysis, reliability analysis and the like of the three-phase active power distribution network.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort.
Fig. 1 is a schematic flow diagram of a three-phase active power distribution network power flow analysis method suitable for a mass scene provided by the present application;
FIG. 2 is a schematic diagram of a modified IEEE-123 bus three-phase unbalanced distribution network provided by an embodiment of the application;
fig. 3 is a variation characteristic diagram of a power proportional control factor variable in the three-phase active power distribution network power flow analysis method adapted to mass scenes provided by the application;
fig. 4 is a schematic view of a power flow result of the three-phase active power distribution network power flow analysis method adapted to mass scenes provided by the present application;
fig. 5 is a schematic structural diagram of a three-phase active power distribution network power flow analysis device adapted to a mass scene provided by the present application;
fig. 6 is a schematic structural diagram of an electronic device provided in the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments of the present specification, and not all of the embodiments. All other embodiments obtained by a person skilled in the art based on the embodiments in the present specification without any inventive step should fall within the scope of protection of the present specification.
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
It will be apparent to those skilled in the art that various modifications and variations can be made in the specific embodiments described herein without departing from the scope or spirit of the application. Other embodiments will be apparent to the skilled person from the description of the present application. The specification and examples are exemplary only.
As used herein, the terms "comprising," "including," "having," "containing," and the like are open-ended terms that mean including, but not limited to.
In the present application, "parts" are in parts by mass unless otherwise specified.
In the related technology, the iterative load flow calculation method must repeatedly perform load flow calculation during massive scene analysis, and the power proportional control factor of the recursive load flow calculation method has a single change direction.
Based on the defects, the embodiment of the application provides the three-phase active power distribution network power flow analysis method suitable for the mass scenes, and power flow calculation of the mass scenes can be rapidly carried out.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Referring to fig. 1, a flow diagram of a three-phase active power distribution network power flow analysis method suitable for mass scenes provided in the embodiment of the present application is shown.
As shown in fig. 1, a three-phase active power distribution network power flow analysis method adapted to a massive scene may include:
s110, obtaining operation state data of the mass scene to be analyzed, wherein the operation state data comprises power proportion control factors corresponding to the mass scene to be analyzed.
Specifically, the mass scene is the operation state of the system, and different mass scenes correspond to different operation state data.
The power proportional control factors are parameters for controlling the power change characteristics of the corresponding distributed power sources or loads, and when the mass scene to be analyzed is determined, the corresponding group of power proportional control factors are known values. And the system is divided into several areas, and correspondingly, there are several power proportion control factors.
And S120, inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
According to the power flow analysis method for the three-phase active power distribution network adapting to the mass scenes, the power flow solution of the mass scenes to be analyzed can be obtained by inputting the running state data of the mass scenes to be analyzed into the offline power flow analysis expression, the power flow solution can be quickly obtained without repeated calculation, the calculation efficiency is higher, and therefore decision support is quickly provided for planning, probability power flow analysis, reliability analysis and the like of the three-phase active power distribution network. The three-phase active power distribution network power flow analysis method suitable for the mass scenes provided by the embodiment can be used for rapidly analyzing the power flow of the mass and any operation scenes.
Specifically, the offline power flow analysis expression is predetermined.
The operation state data is input parameters of an offline power flow analysis expression, and power flow of a mass scene to be analyzed is resolved into output values.
In one embodiment, the step of determining the offline power flow analysis expression comprises:
dividing the system into at least two areas according to the type of a bus of the three-phase active power distribution network and the geographic position of the distributed power supply or the load, wherein each area corresponds to at least one independent power proportional control factor variable, and the power proportional control factor variables are respectively used for representing the power change characteristics of the distributed power supply or the load in the corresponding area;
based on a multidimensional pure embedding method, an offline power flow analytical expression with a plurality of power proportional control factor variables as multivariable is determined.
In particular, distributed power sources or loads in the same area in a three-phase active power distribution network tend to have strong spatial correlation, and their characteristics in different areas are very different. Therefore, the whole system is divided into a plurality of areas in consideration of factors such as bus type, geographic position and the like, and a plurality of independent power scale factor variables are adopted to represent power change characteristics of distributed power supplies or loads in different areas. The division of the areas can be divided according to the actual geographic requirements of the power distribution network system. For example, the power distribution network system is a county power distribution network, and may be divided into a production type area, a business type area, a general type area, and the like according to the type and function of a power utilization building in the county.
The multidimensional pure embedding method is a load flow calculation method based on complex analysis, a recursion theory and analytic extension. The original implicit power balance equation is converted into an explicit tidal flow expression through the steps of multivariate embedding, recursive calculation, generalized approximation and the like.
In one embodiment, the method for determining the offline power flow analytic expression with a plurality of power proportional control factor variables as multivariable based on the multidimensional full-pure embedding method comprises the following steps:
according to a plurality of power proportion control factor variables, representing each bus voltage in the three-phase active power distribution network as a multivariable pure function with unknown coefficients based on a multidimensional pure embedding method;
determining an initial power balance equation of each bus according to the connection mode of the buses;
determining a fully-pure embedding equation according to a multivariable fully-pure function and an initial power balance equation;
and determining an unknown coefficient of the offline power flow analytic expression according to the multivariable holomorphic function and the holomorphic embedding equation. The bus bar connection mode may include Y-connection, delta-connection, or the like.
Determining an unknown coefficient of an offline power flow analytic expression according to a multivariable holohedral function and a holohedral embedding equation may include:
substituting a multivariable pureness function into a pureness embedding equation to obtain a group of equation sets;
after the derivation transformation of the equation set, obtaining a recursion equation through the equality of coefficients with the same power on both sides of the equation;
and solving the unknown coefficient of the offline power flow analytic expression according to a recursion equation.
Determining an offline power flow analytic expression according to an unknown coefficient of a multivariable holohedral function may include:
determining a fractional expression of a multivariable pure function according to a generalized rational approximation technology;
determining the coefficient of the fractional expression according to the unknown coefficient of the multivariable holo-pure function;
and determining an offline power flow analysis expression according to the coefficient of the fractional expression.
Specifically, assume that a three-phase unbalanced active power distribution network has N buses (one slack bus and N-1 PQ buses) and D zones, where both the distributed power supply and the load are PQ buses, and a Y-type connection or a delta-type connection is used. The power balance equation (i.e. the initial power balance equation of the bus) of the three-phase unbalanced active power distribution network is as follows:
wherein
Of m-phase voltage at node k and of node i, respectively
The phase voltages are set to be equal to each other,
is the conjugate of the phase-to-phase voltages of the z-phase and x-phase of node i,
being node i
Admittance between the phase and the m-phase of node k,
respectively injected into node i
The active power and the reactive power of the phases,
the active power and the reactive power of the x phase and the y phase of the node i are respectively,
respectively being loose bus-bar
The voltage amplitude of the phase is the conjugate of the complex number with the phase angle, (·).
Because the three-phase active power distribution network is divided into D areas, two complex variables can be embedded into a power balance equation of each area to respectively control the changes of active power and reactive power of the three-phase active power distribution network. Therefore, based on the multidimensional all-pure embedding method, the unknown phase voltage
And phase-to-phase voltage
The fully pure function to be expressed as multivariate on 2D variables is as follows:
wherein s is 1 To s D And s D+1 To s 2D Respectively representing power proportional control factor variables corresponding to active power and reactive power; n is e Is s e To the power of (c);
is an expression
In
The coefficient of the term. Based on this, the unknown coefficients are solved
Becoming the current primary task.
All-pure function of multiple variables
And
introducing equations (1) and (2), the all-pure embedding equation of the initial power balance equation can be written as:
wherein x (i) represents the zone index to which the bus i belongs; s x(i) And s D+x(i) Respectively representing the power proportional control factor variables of the active power and the reactive power at the bus i;
and
respectively given initial operating conditions
And
the value of (c).
Are used separately
And
replacement of
And
the reciprocal of (a) is then:
wherein
Are respectively an expression
In
The coefficient of the term.
Substituting (8) and (9) into (6) and (7) to obtain:
since both sides of the equations (10) - (11) are related to(s) 1 ,s 2 ,...,s 2D ) The recursion equation can be obtained by making the coefficients of the same term on both sides of the equation equal as follows:
equations (12) and (13) are expressions
And
the recurrence relation between the coefficients of (c).
It can be found that the coefficients of the order to be solved
And
they can always be represented by the found low order coefficients. In addition, the method can be used for producing a composite material
And
and
the relationship (2) can be further obtained by (8) to (9). Let equation
The coefficients of the same terms on both sides are equal, and the following can be obtained:
wherein
And
is that
And
the abbreviations of (a), which are constant terms of their corresponding power series, are readily available under initial operating conditions; m represents
The overall order of the term, Conv (-) for M ≧ 2, can be specifically expressed as:
K＝k 1 +k 2 +…+k 2D ,0≤k s ≤n s ,s＝1～2D (18)
in conclusion, the (12) to (18) jointly form a recurrence equation of the rapid power flow analysis technology of the three-phase active power distribution network. The coefficients of the solution of any order (left side of the equation) can be represented linearly by coefficients of a lower order. Once the physical initial solution of the power flow is determined, the unknown coefficients of the multivariate all-pure function can be computed step by step from low order to high order. And then substituting the unknown coefficients of the multivariable pure function into the multivariable pure function to obtain an off-line power flow analytical expression.
However, the multivariate perfectly pure function is essentially a multivariate function expression with infinite coefficients, and is a completely accurate expression, but in practical calculation, the infinite coefficients are difficult to obtain completely, and the above processes (12) - (18) can only obtain the multivariate function expression with finite coefficients, which is called a truncated series, and the solution is performed by the truncated series instead of the multivariate perfectly pure function in the present application. Thus, cutting offThe number of stages is not precise. In order to obtain a more accurate function expression than the truncated series, it is necessary to apply a multivariable function expression having finite term coefficients (i.e., the truncated series, in equation (19))
) Using generalized rational approximation technique to convert into fractional expression (in formula (19))
) Then, the unknown coefficient of the fractional expression is solved according to the coefficient of the truncation series
And
coefficient of fractional expression
And
the coefficient of the truncated series obtained by the solution is substituted into equation (19), and the coefficients of the same terms on both sides of equation (19) are made equal. Coefficient of the obtained fractional expression
And
substituting into a fractional expression, namely an off-line tide analysis expression, wherein the expression is more accurate than a truncation series.
Experimental verification
(1) IEEE-123 bus distribution network verification
Three-phase active power distribution adapting to mass scenesThe power grid power flow analysis method is verified on the modified IEEE-123 bus power distribution network, and the topology and the configuration of the distributed power supply are shown in figure 2. The numerical result may be in the form of
Core TM i5-11300H CPU @3.1GHz was obtained on a computer using MATLAB.
Dividing the system into 6 areas according to the spatial correlation of the distributed power supply and the load, and embedding a power proportional control factor variable s 1 ～s 6 To independently specify regional power variations. Based on a three-phase active power distribution network power flow analysis method adaptive to massive scenes, a multivariable holohedral function expression containing 1716 variables and an offline power flow analysis expression (namely a fractional expression) with the maximum power of 6 are obtained through calculation.
Table 1 shows a comparison between the three-phase active power distribution network power flow analysis method and the forward-backward substitution method, which are suitable for mass scenes. The maximum error in table 1 is a difference between a calculation result (including a voltage amplitude V, a voltage phase angle θ, an active power P, and a reactive power Q) of the three-phase active power distribution network power flow analysis method adapted to the massive scenes and a calculation result of a forward-backward substitution method, and as can be seen from table 1, the maximum errors are very small and are approximately close to 0, so that the three-phase active power distribution network power flow analysis method adapted to the massive scenes has high calculation accuracy.
The total calculation time of the three-phase active power distribution network power flow analysis method suitable for massive scenes can be divided into the solving time (corresponding to step 1 in table 1) of a physical initial solution (namely, a power flow solution in an initial specified running state, including voltage amplitude, voltage phase angle, active power and reactive power), the recursive calculation (namely, formula (12) -formula (18)) time (corresponding to step 2 in table 1), the rational approach (namely, formula (19)) time (corresponding to step 3 in table 1) in the embodiment of the application) and the calculation time (corresponding to step 4 in table 1) of the power flow solution in the massive scenes to be analyzed. In single scene analysis, the total calculation time of the three-phase active power distribution network power flow analysis method suitable for massive scenes is longer than that of a forward-backward substitution method. However, the three-phase active power distribution network power flow analysis method suitable for massive scenes can execute the step 1 to the step 3 in advance in an off-line mode, and an off-line power flow expression is obtained in advance. The three steps, step 1 to step 3, are fixed and only calculated once, once given the network topology. Therefore, the power flow calculation time in a single scenario is only the time for substituting the numerical value (i.e., the operation state data) of the distributed power supply or the load into the offline power flow expression (i.e., the online calculation time of 0.0006 seconds in step 4), which is shorter than the calculation time of 0.0035 seconds in the forward-backward substitution method.
TABLE 1 comparison of the calculated Performance of the present application with the Forward pushback substitution method
In addition, power flow analysis under a massive scene is performed, and as shown in fig. 3, a power proportional control factor variable s in the power flow analysis method of the three-phase active power distribution network applicable to the massive scene is shown 1 ～s 6 The change profile of (2). Fig. 4 is a schematic view of a power flow result of the three-phase active power distribution network power flow analysis method adapted to mass scenes.
For 1X 10 5 In each operation scene, the online calculation time of the three-phase active power distribution network power flow analysis method suitable for the massive scenes is only 2.52 seconds, while the online calculation time of the forward-backward substitution method is 340.83 seconds, which fully shows that the three-phase active power distribution network power flow analysis method suitable for the massive scenes has higher efficiency in power flow analysis of the massive operation scenes.
(2) Mass scene load flow calculation performance analysis of power distribution networks of different scales
Table 2 compares the performance of the three-phase active power distribution network power flow analysis method adapted to the massive scene in the calculation of four unbalanced three-phase active power distribution networks (including modified 37 bus, 123 bus, 906 bus and 1811 bus power distribution networks) with 5 other power flow calculation methods, so as to compare the accuracy and calculation time of power flow calculation in the massive operation scene.
Table 2 comparison of the calculation performance of the present application and other five load flow calculation methods in four power distribution networks of different scales
Note: 'BFS' -forward-backward substitution, 'MF _ BFS' -forward-backward substitution in a matrix form, 'DA' -direct method, 'LPF' -linearized power flow method, 'CIM' -injection current method.
As can be seen from table 2, in the power flow calculation, besides the linearized power flow method, the three-phase active power distribution network power flow analysis method applicable to the massive scenes has the same high calculation accuracy as other methods, and the calculation efficiency of the three-phase active power distribution network power flow analysis method applicable to the massive scenes is much higher than that of other methods. Particularly, in a large-scale 1811 bus three-phase active power distribution network, the three-phase active power distribution network power flow analysis method suitable for massive scenes calculates 1 x 10 5 The time (35.86s) spent in each operation scene is still far lower than that of other methods (more than or equal to 1115s), and the three-phase active power distribution network power flow analysis method suitable for the mass scenes is fully suitable for rapid power flow calculation of the mass scenes of the large-scale system.
Referring to fig. 5, a schematic structural diagram of a three-phase active power distribution network power flow analysis device adapting to a massive scene is shown according to an embodiment of the present application.
As shown in fig. 5, a three-phase active power distribution network power flow analysis device 500 adapted to a massive scene may include:
an obtaining module 510, configured to obtain operation state data of a mass scene to be analyzed, where the operation state data includes a power proportional control factor corresponding to the mass scene to be analyzed;
the processing module 520 is configured to input the operation state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
Optionally, the three-phase active power distribution network power flow analysis device 500 adapted to the massive scenes further includes a determining module, configured to determine an offline power flow analysis expression, including:
dividing the system into at least two areas according to the type of a bus of the three-phase active power distribution network and the geographic position of the distributed power supply or the load, wherein each area corresponds to at least one independent power proportional control factor variable, and the power proportional control factor variables are respectively used for representing the power change characteristics of the distributed power supply or the load in the corresponding area;
based on a multidimensional pure embedding method, an offline power flow analytical expression with a plurality of power proportional control factor variables as multivariable is determined.
Optionally, the determining module is further configured to:
according to a plurality of power proportional control factor variables, representing the voltage of each bus in the three-phase active power distribution network as a multivariable pure function with unknown coefficients based on a multidimensional pure embedding method;
determining an initial power balance equation of each bus according to the connection mode of the buses;
determining a fully-pure embedding equation according to a multivariable fully-pure function and an initial power balance equation;
determining unknown coefficients of the multivariable holomorphic functions according to the multivariable holomorphic functions and the holomorphic embedding equations;
and determining an offline power flow analytical expression according to the unknown coefficient of the multivariable holomorphic function.
Optionally, the determining module is further configured to:
substituting a multivariable pureness function into a pureness embedding equation to obtain a group of equation sets;
after the derivation transformation of the equation set, obtaining a recursion equation through the equality of coefficients with the same power on both sides of the equation;
and solving the unknown coefficient of the multivariable holo-pure function according to a recurrence equation.
Optionally, the determining module is further configured to:
determining a fractional expression of a multivariable pure function according to a generalized rational approximation technology;
determining the coefficient of the fractional expression according to the unknown coefficient of the multivariable holomorphic function;
and determining an offline power flow analysis expression according to the coefficient of the fractional expression.
Optionally, the bus bar connection mode includes Y-type connection or delta-type connection.
The three-phase active power distribution network power flow analysis device suitable for massive scenes provided by the embodiment can execute the embodiment of the method, the implementation principle and the technical effect are similar, and the details are not repeated herein.
Fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the present invention. As shown in fig. 6, a schematic structural diagram of an electronic device 300 suitable for implementing the embodiments of the present application is shown.
As shown in fig. 6, the electronic apparatus 300 includes a Central Processing Unit (CPU)301 that can perform various appropriate actions and processes in accordance with a program stored in a Read Only Memory (ROM)302 or a program loaded from a storage section 308 into a Random Access Memory (RAM) 303. In the RAM 303, various programs and data necessary for the operation of the apparatus 300 are also stored. The CPU 301, ROM 302, and RAM 303 are connected to each other via a bus 304. An input/output (I/O) interface 305 is also connected to bus 304.
The following components are connected to the I/O interface 305: an input portion 306 including a keyboard, a mouse, and the like; an output section 307 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage section 308 including a hard disk and the like; and a communication section 309 including a network interface card such as a LAN card, a modem, or the like. The communication section 309 performs communication processing via a network such as the internet. A drive 310 is also connected to the I/O interface 306 as needed. A removable medium 311 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 310 as necessary, so that a computer program read out therefrom is mounted into the storage section 308 as necessary.
In particular, the process described above with reference to fig. 1 may be implemented as a computer software program, according to an embodiment of the present disclosure. For example, embodiments of the present disclosure include a computer program product comprising a computer program tangibly embodied on a machine-readable medium, the computer program comprising program code for performing the above-described three-phase active power distribution network power flow resolution method that accommodates massive scenarios. In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 309, and/or installed from the removable medium 311.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units or modules described in the embodiments of the present application may be implemented by software or hardware. The described units or modules may also be provided in a processor. The names of these units or modules do not in some cases constitute a limitation of the unit or module itself.
The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. One typical implementation device is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a mobile phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
As another aspect, the present application also provides a storage medium, which may be the storage medium contained in the foregoing device in the above embodiment; or may be a storage medium that exists separately and is not assembled into the device. The storage medium stores one or more programs, and the programs are used by one or more processors to execute the three-phase active power distribution network power flow analysis method adaptive to the massive scenes.
Storage media, including permanent and non-permanent, removable and non-removable media, may implement the information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It is to be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in the process, method, article, or apparatus that comprises the element.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, as for the system embodiment, since it is substantially similar to the method embodiment, the description is relatively simple, and reference may be made to the partial description of the method embodiment for relevant points.
Claims (9)
1. A three-phase active power distribution network power flow analysis method adaptive to massive scenes is characterized by comprising the following steps:
acquiring running state data of a mass scene to be analyzed, wherein the running state data comprises a power proportion control factor corresponding to the mass scene to be analyzed;
and inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
2. The method of claim 1, wherein the step of determining the offline power flow analysis expression comprises:
dividing the system into at least two areas according to the type of a bus of the three-phase active power distribution network and the geographic position of the distributed power supply or the load, wherein each area corresponds to at least one independent power proportional control factor variable, and the power proportional control factor variables are respectively used for representing the power change characteristics of the distributed power supply or the load in the corresponding area;
and determining the offline power flow analytical expression with the power proportional control factor variables as multivariable based on a multidimensional pure embedding method.
3. The method of claim 2, wherein the determining the offline power flow analytic expression with the power proportional control factor variables as multivariate based on the multidimensional all-pure embedding method comprises:
according to the power proportional control factor variables, representing the voltage of each bus in the three-phase active power distribution network as a multivariable pure function with unknown coefficients based on a multidimensional pure embedding method;
determining an initial power balance equation of each bus according to the connection mode of the buses;
determining a fully-pure embedding equation according to the multivariable fully-pure function and the initial power balance equation;
determining unknown coefficients of the multivariable holo-pure function according to the multivariable holo-pure function and the holo-pure embedding equation;
and determining the offline power flow analytical expression according to the unknown coefficient of the multivariable holomorphic function.
4. The method of claim 3, wherein said determining unknown coefficients of said multivariate holopure function from said multivariate holopure function and said holopure embedding equation comprises:
substituting the multivariable pureness function into the pureness embedding equation to obtain a group of equation sets;
after the derivation transformation of the equation set, obtaining a recursion equation through the equality of coefficients with the same power on both sides of the equation;
and solving the unknown coefficient of the multivariable holomorphic function according to the recursion equation.
5. The method of claim 3, wherein determining the offline power flow analytic expression according to the unknown coefficients of the multivariate holograte function comprises:
determining a fractional expression by the multivariable holo-pure function according to a generalized rational approximation technology;
determining the coefficient of the fractional expression according to the unknown coefficient of the multivariable holomorphic function;
and determining the offline power flow analysis expression according to the coefficient of the fractional expression.
6. The method of claim 3, wherein the bus bars are connected in a manner comprising a Y-connection or a delta-connection.
7. The utility model provides a three-phase active distribution network trend analytical equipment of adaptation magnanimity scene which characterized in that, the device includes:
the system comprises an acquisition module, a processing module and a processing module, wherein the acquisition module is used for acquiring operation state data of a mass scene to be analyzed, and the operation state data comprises a power proportion control factor corresponding to the mass scene to be analyzed;
and the processing module is used for inputting the running state data into an offline power flow analysis expression to obtain a power flow solution of the mass scene to be analyzed.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the method for flow analysis of a three-phase active power distribution network according to any of claims 1 to 6 adapted to mass scenarios.
9. A readable storage medium, on which a computer program is stored, wherein the program, when executed by a processor, implements the method for power flow analysis of a three-phase active power distribution network according to any one of claims 1 to 6, adapted to mass scenarios.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210415820.1A CN114912767A (en) | 2022-04-20 | 2022-04-20 | Three-phase active power distribution network power flow analysis method and device suitable for mass scenes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210415820.1A CN114912767A (en) | 2022-04-20 | 2022-04-20 | Three-phase active power distribution network power flow analysis method and device suitable for mass scenes |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114912767A true CN114912767A (en) | 2022-08-16 |
Family
ID=82765318
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210415820.1A Pending CN114912767A (en) | 2022-04-20 | 2022-04-20 | Three-phase active power distribution network power flow analysis method and device suitable for mass scenes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114912767A (en) |
-
2022
- 2022-04-20 CN CN202210415820.1A patent/CN114912767A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zakzouk et al. | Improved performance low‐cost incremental conductance PV MPPT technique | |
| Tang et al. | Dimension-adaptive sparse grid interpolation for uncertainty quantification in modern power systems: Probabilistic power flow | |
| CN104156542B (en) | It is a kind of based on the active distribution system Simulation of stability method implicitly projected | |
| Zhang et al. | Multiobjective optimization of a fractional‐order pid controller for pumped turbine governing system using an improved nsga‐Iii algorithm under multiworking conditions | |
| Yuan et al. | Research on dynamic load modelling based on power quality monitoring system | |
| Xia et al. | Computationally efficient multi‐step direct predictive torque control for surface‐mounted permanent magnet synchronous motor | |
| Mareček et al. | A low-rank coordinate-descent algorithm for semidefinite programming relaxations of optimal power flow | |
| Sang et al. | Design and implementation of perturbation observer‐based robust passivity‐based control for VSC‐MTDC systems considering offshore wind power integration | |
| Abdel‐Akher et al. | Initialised load‐flow analysis based on Lagrange polynomial approximation for efficient quasi‐static time‐series simulation | |
| CN111327048A (en) | Robust operation optimization method for power distribution network containing three-terminal SNOP | |
| Zhang et al. | Short‐circuit current constrained unit commitment and transmission switching model for improving renewable integration: An MILP formulation | |
| Luo et al. | Direct prediction compensation strategy of unified power quality conditioner based on FCS–MPC | |
| Meng et al. | Stability analysis and control optimisation based on particle swarm algorithm of modular multilevel matrix converter in fractional frequency transmission system | |
| Iweh et al. | Assessment of the optimum location and hosting capacity of distributed solar PV in the southern interconnected grid (SIG) of Cameroon | |
| Sarker et al. | Modified harmonic minimisation technique for doubly fed induction generators with solar‐wind hybrid system using biogeography‐based optimisation | |
| CN114912767A (en) | Three-phase active power distribution network power flow analysis method and device suitable for mass scenes | |
| Jin et al. | Model predictive control strategy for NPC grid‐connected inverters in unbalanced grids | |
| CN116865318A (en) | Power transmission network and energy storage joint planning method and system based on two-stage random optimization | |
| Bhujade et al. | Economic analysis of annual load loss due to voltage sags in industrial distribution networks with distributed PVs | |
| Wang et al. | Unbalanced responsibility division considering renewable energy integration | |
| Hong et al. | Optimal design of IT2‐FCS‐based STATCOM controller applied to power system with wind farms using Taguchi method | |
| Hussain et al. | Real‐time implementation of three‐phase single‐stage SPV grid‐tied system using TL‐VSC | |
| Chen et al. | Dynamic optimal power flow model incorporating interval uncertainty applied to distribution network | |
| CN110707689B (en) | Stability analysis method and device suitable for full-clean energy power generation power grid | |
| Jabari et al. | Probabilistic power flow analysis of distribution systems using Monte Carlo simulations |
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 |