CN114188945B - Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply - Google Patents

Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply Download PDF

Info

Publication number
CN114188945B
CN114188945B CN202210144000.3A CN202210144000A CN114188945B CN 114188945 B CN114188945 B CN 114188945B CN 202210144000 A CN202210144000 A CN 202210144000A CN 114188945 B CN114188945 B CN 114188945B
Authority
CN
China
Prior art keywords
node
photovoltaic
current
voltage
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210144000.3A
Other languages
Chinese (zh)
Other versions
CN114188945A (en
Inventor
潘本仁
张妍
邹进
谢国强
万勇
周宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202210144000.3A priority Critical patent/CN114188945B/en
Publication of CN114188945A publication Critical patent/CN114188945A/en
Application granted granted Critical
Publication of CN114188945B publication Critical patent/CN114188945B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • 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
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/22—The renewable source being solar energy
    • H02J2300/24—The renewable source being solar energy of photovoltaic origin
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method and a device for calculating short-circuit current of a power distribution network containing a photovoltaic power supply, wherein the method comprises the steps of calculating to obtain a normal operation value of voltage and current of a photovoltaic node according to load flow before a fault; respectively carrying out equivalent admittance processing suitable for Gaussian iterative short circuit calculation on the load nodes and the fault nodes to form an admittance array containing the photovoltaic power distribution network; calculating and correcting the photovoltaic node injection current in the injection current vector according to the photovoltaic node voltage after each Gaussian iteration based on the voltage-controlled current source characteristic of the photovoltaic power supply; and after the short circuit calculation iterative convergence, the voltage value of each node after the fault is obtained, and the short circuit current of the whole network is obtained according to the voltage value. The method can better embed the Gaussian iteration algorithm into the short circuit calculation process, and improve the calculation precision of the short circuit current of the power distribution network containing the photovoltaic power supply.

Description

Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply
Technical Field
The invention belongs to the technical field of electric power system analysis, and particularly relates to a method and a device for calculating short-circuit current of a power distribution network with a photovoltaic power supply.
Background
At present, in order to solve the increasingly severe energy crisis and environmental problems, the new energy power generation technology represented by photovoltaic and wind power is rapidly developed, and by the end of 2020, the installed capacity of new energy in China is up to 43%, wherein the installed capacity of grid-connected photovoltaic is up to 25343 ten thousand kilowatts and is 12% of the total installed capacity. After the photovoltaic power supply is connected to a power grid in a high-density mode, the topological structure and the tide distribution of the power grid are changed, and the safe operation level of the system under the disturbance of the short-circuit fault of the power grid is influenced. However, in the current process of calculating the short-circuit current of the power distribution network and setting the protection configuration, new energy power supplies such as wind power and photovoltaic power supplies are generally simply viewed as constant current sources or loads, which is inconsistent with the actual operation condition, and the accuracy of calculating the fault current after the short circuit and the action performance of the relay protection device are seriously influenced.
With the rapid increase of the power access proportion of photovoltaic power supplies and the like, because the output current of the photovoltaic power supply is related to the grid-connected point voltage, the short-circuit characteristic of the photovoltaic power supply cannot be considered by the conventional short-circuit current calculation method, and the calculation accuracy of the short-circuit current cannot meet the requirement of actual production operation, an iterative calculation method capable of improving the calculation accuracy of the short-circuit current is urgently needed to be researched.
Disclosure of Invention
The invention provides a method and a device for calculating short-circuit current of a power distribution network containing a photovoltaic power supply, which are used for solving at least one of the technical problems.
In a first aspect, the invention provides a method for calculating a short-circuit current of a power distribution network containing a photovoltaic power supply, which comprises the following steps: step A: the photovoltaic node is regarded as a PQ node, and the power factor is1, obtaining node voltage and injection current of each node including photovoltaic nodes in normal operation by using load flow calculation, wherein the voltage vector of the ith node is
Figure 542229DEST_PATH_IMAGE001
The injection current vector of the i-th node is
Figure 578449DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 112199DEST_PATH_IMAGE003
(ii) a And B: after the system generates three-phase symmetrical grounding short circuit, the fault point is regarded as a newly added node in the network, and the load flow calculation node admittance array is modified into the photovoltaic short circuit current calculation admittance array by adopting an impedance simulation method suitable for the Gaussian iterative photovoltaic short circuit calculation
Figure 166742DEST_PATH_IMAGE004
Constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation; and C: obtaining a photovoltaic node voltage vector according to the k-1 iteration
Figure 127745DEST_PATH_IMAGE005
According to the voltage-controlled current source characteristic of the photovoltaic power supply, the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage is obtained through calculation
Figure 149797DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 221658DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 396287DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 13344DEST_PATH_IMAGE009
The short-circuit current of each branch can be obtained.
In a second aspect, the present invention provides a device for calculating short-circuit current of a power distribution network including a photovoltaic power supply, including: the calculation module is configured to regard the photovoltaic node as a PQ node, the power factor is 1, and the node voltage and the injection current of each node including the photovoltaic node during normal operation are obtained by utilizing load flow calculation, wherein the voltage vector of the ith node is PQ node
Figure 273424DEST_PATH_IMAGE001
The injection current vector of the i-th node is
Figure 883397DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 912533DEST_PATH_IMAGE003
(ii) a And the modification module is configured to modify the load flow calculation node admittance array into the photovoltaic short-circuit current calculation admittance array by adopting an impedance simulation method suitable for Gaussian iterative photovoltaic short-circuit calculation after the three-phase symmetrical grounding short circuit of the system occurs and the fault point is regarded as a newly-added node in the network
Figure 464606DEST_PATH_IMAGE004
Constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation; an iteration module configured to obtain a photovoltaic node voltage vector according to the k-1 th iteration
Figure 211982DEST_PATH_IMAGE005
Calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage according to the voltage-controlled current source characteristic of the photovoltaic power supply
Figure 625646DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 509288DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 733728DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 702821DEST_PATH_IMAGE009
The short-circuit current of each branch can be obtained.
In a third aspect, an electronic device is provided, comprising: the photovoltaic power distribution network short-circuit current calculation method comprises at least one processor and a memory which is in communication connection with the at least one processor, wherein the memory stores instructions which can be executed by the at least one processor, and the instructions are executed by the at least one processor, so that the at least one processor can execute the steps of the photovoltaic power distribution network short-circuit current calculation method according to any embodiment of the invention.
In a fourth aspect, the present invention also provides a computer-readable storage medium, on which a computer program is stored, the computer program comprising program instructions which, when executed by a computer, cause the computer to perform the steps of the method for calculating a short-circuit current of a distribution network including a photovoltaic power supply according to any of the embodiments of the present invention.
According to the method and the device for calculating the short-circuit current of the power distribution network containing the photovoltaic power supply, a Gaussian iteration algorithm can be well embedded into the short-circuit calculation process, and the calculation precision of the short-circuit current of the power distribution network containing the new energy power supply is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a flowchart of a method for calculating a short-circuit current of a power distribution network including a photovoltaic power supply according to an embodiment of the present invention;
FIG. 2 is a flow chart of a Gaussian iterative solution provided by an embodiment of the present invention;
fig. 3 is a block diagram illustrating a short-circuit current calculation apparatus for a power distribution network including a photovoltaic power supply according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. 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.
Referring to fig. 1, a flowchart of a method for calculating a short-circuit current of a power distribution network including a photovoltaic power supply according to the present application is shown.
As shown in fig. 1, the method for calculating the short-circuit current of the distribution network including the photovoltaic power supply specifically includes the following steps:
step A, regarding the photovoltaic node as a PQ node, wherein the power factor is 1, and obtaining the node voltage and the injection current of each node including the photovoltaic node during normal operation by utilizing load flow calculation, wherein the voltage vector of the ith node is PQ node
Figure 654596DEST_PATH_IMAGE001
The injection current vector of the i-th node is
Figure 173171DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 817779DEST_PATH_IMAGE003
。
Step B, after the three-phase symmetrical grounding short circuit occurs in the system, the fault point is regarded as a newly added node in the network, and the method is suitable for Gaussian iterationThe impedance simulation method for photovoltaic short circuit calculation modifies the load flow calculation node admittance array into the photovoltaic short circuit current calculation admittance array
Figure 274168DEST_PATH_IMAGE004
And constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation.
In the present embodiment, step B1: the node load simulation admittance obtained by using the node injection power and the node voltage under the normal network state is expressed as follows:
Figure 29634DEST_PATH_IMAGE010
，
in the formula (I), the compound is shown in the specification,
Figure 373022DEST_PATH_IMAGE011
the admittance is simulated for the load at the ith node,
Figure 188532DEST_PATH_IMAGE012
being the conjugate of the injected power at the ith node,
Figure 397796DEST_PATH_IMAGE013
is the voltage vector of the i-th node,
Figure 956953DEST_PATH_IMAGE014
is the conjugate of the voltage of the ith node;
after equivalent load admittance, the node injects a current vector
Figure 653383DEST_PATH_IMAGE015
；
Step B2: setting a fault point three-phase grounding admittance as
Figure 374214DEST_PATH_IMAGE016
；
Step B3: after the node analog impedance and the three-phase grounding admittance are determined, summing photovoltaic node self-admittance and photovoltaic node load analog admittance in the admittance array, and summing short-circuit fault node self-admittance and short-circuit fault node grounding admittance to obtain a photovoltaic short-circuit current admittance array expression:
Figure 805195DEST_PATH_IMAGE017
，
in the formula (I), the compound is shown in the specification,
Figure 902464DEST_PATH_IMAGE018
for the self-admittance of the node 1 in the load flow calculation,
Figure 751603DEST_PATH_IMAGE019
for the self-admittance of node n in the load flow calculation,
Figure 377756DEST_PATH_IMAGE020
for the self-admittance of the photovoltaic node 1 in the power flow calculation,
Figure 296034DEST_PATH_IMAGE021
for the mutual admittance of the photovoltaic node m and the photovoltaic node n in the load flow calculation,
Figure 196994DEST_PATH_IMAGE022
simulating admittance for the load of node n;
step B4: the flow calculation before the fault obtains the vector of the current injected into the node 1 as
Figure 891453DEST_PATH_IMAGE023
Node 1 voltage vector of
Figure 219666DEST_PATH_IMAGE024
Wherein the node 1 is a balance node, and calculating the potential vector of the generator
Figure 359660DEST_PATH_IMAGE025
The expression of (a) is:
Figure 64311DEST_PATH_IMAGE026
，
in the formula，
Figure 356883DEST_PATH_IMAGE027
Is the equivalent reactance of the generator,
Figure 855998DEST_PATH_IMAGE028
is mathematical imaginary symbols;
generator potential vector before and after supposing short-circuit fault
Figure 483288DEST_PATH_IMAGE025
Constant, equivalent reactance
Figure 726051DEST_PATH_IMAGE027
Without change, the fault-time node 1 injection current vector can be expressed as:
Figure 371664DEST_PATH_IMAGE029
，
in the formula (I), the compound is shown in the specification,
Figure 41680DEST_PATH_IMAGE030
in order to be a function of the sign,
Figure 156267DEST_PATH_IMAGE031
is the voltage of node 1;
a function of a photovoltaic node injection current vector and a photovoltaic node voltage vector, expressed as:
Figure 937141DEST_PATH_IMAGE032
，
Figure 204305DEST_PATH_IMAGE033
a voltage vector of a photovoltaic node;
step B5: after the injection currents of the balance node and the photovoltaic node are determined, a node injection current vector is obtained
Figure 45222DEST_PATH_IMAGE034
Expression:
Figure 647105DEST_PATH_IMAGE035
，
step B6: the expression of the node impedance network equation is obtained as follows:
Figure 231670DEST_PATH_IMAGE036
wherein, in the step (A),
Figure 851876DEST_PATH_IMAGE037
,
in the formula (I), the compound is shown in the specification,
Figure 863695DEST_PATH_IMAGE038
is the injected current vector of the photovoltaic node 1,
Figure 687294DEST_PATH_IMAGE039
is the vector of the injected current at photovoltaic node m,
Figure 75550DEST_PATH_IMAGE040
is the voltage vector of the photovoltaic node 1,
Figure 51727DEST_PATH_IMAGE041
is the voltage vector of the photovoltaic node m.
In the method of this embodiment, in order to effectively calculate the short-circuit current by using the Zbus gaussian iteration algorithm, the impedance simulation method suitable for the photovoltaic short-circuit calculation of the gaussian iteration is adopted to modify the load flow calculation node admittance array into the short-circuit current calculation admittance array
Figure 234447DEST_PATH_IMAGE042
And writing node impedance network equations suitable for Gaussian iterative photovoltaic short circuit calculation in parallel, and then iteratively solving to obtain each node voltage.
Step C, obtaining a photovoltaic node voltage vector according to the k-1 iteration
Figure 810922DEST_PATH_IMAGE043
According to the voltage-controlled current source characteristic of the photovoltaic power supply, the voltage-controlled current source characteristic is calculated and obtained under the directional control of the d-axis voltageOutput current vector of photovoltaic inverter
Figure 737290DEST_PATH_IMAGE044
Make the node after the k-1 iteration of replacement inject the current vector
Figure 551662DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 420130DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 218321DEST_PATH_IMAGE009
Thus, the short-circuit current of each branch can be obtained.
In this embodiment, step C1: assuming that the output of the photovoltaic power supply is unchanged, setting the short-circuit voltage value of each node obtained by solving after the fault as the voltage vector iteration initial value of the Gaussian iteration short-circuit current node admittance equation
Figure 948380DEST_PATH_IMAGE045
Wherein:
Figure 882838DEST_PATH_IMAGE046
，
in the formula (I), the compound is shown in the specification,
Figure 423672DEST_PATH_IMAGE047
iteration initial values of the n voltage of the node are obtained;
step C2: the voltage iteration formula is as follows:
Figure 709160DEST_PATH_IMAGE048
，
in the formula (I), the compound is shown in the specification,
Figure 977330DEST_PATH_IMAGE049
the voltage at node i after the (k + 1) th iteration,
Figure 31874DEST_PATH_IMAGE050
for the voltage at node i after the kth iteration,
Figure 242144DEST_PATH_IMAGE051
in order to inject the current, the current is injected,
Figure 749349DEST_PATH_IMAGE052
is a photovoltaic short-circuit current admittance array;
step C3: before each iteration, according to the photovoltaic node voltage vector obtained after the k-1 iteration
Figure 821210DEST_PATH_IMAGE053
And calculating the voltage value drop degree of the photovoltaic node
Figure 995839DEST_PATH_IMAGE054
，
Figure 612897DEST_PATH_IMAGE055
Is the voltage vector of the photovoltaic node m after the (k-1) th iteration,
Figure 607397DEST_PATH_IMAGE056
for the voltage vector of the photovoltaic node m, according to the voltage-controlled source characteristic of the photovoltaic output current, when the voltage drop degree is more than or equal to 0.9, the photovoltaic power supply outputs active current with a normal operation value and does not provide reactive current, when the voltage drop degree is less than 0.9 and more than or equal to 0.2, the voltage drops by 1% every time, the photovoltaic power supply provides 1.5% of reactive current, before the current reaches the amplitude limit, the active current output is the normal operation value, after the current reaches the amplitude limit, the active current output is reduced, when the voltage drop degree is less than 0.2, the photovoltaic power supply does not provide active current, only provides reactive current with the amplitude of 1.05 times of the amplitude of the normal operation value, calculates the photovoltaic output current in the current iteration, replaces the node injection current vector after the k-1 iteration, and
Figure 482949DEST_PATH_IMAGE057
photovoltaic node m injected current vector
Figure 246506DEST_PATH_IMAGE058
Obtaining the node injection current vector after the kth iteration
Figure 533000DEST_PATH_IMAGE059
Wherein the photovoltaic node m after the kth iteration injects a current vector
Figure 280376DEST_PATH_IMAGE060
The mode length and the phase angle are respectively:
Figure 694040DEST_PATH_IMAGE061
，
the phase angle of the photovoltaic output current is:
Figure 577682DEST_PATH_IMAGE062
，
in the formula (I), the compound is shown in the specification,
Figure 802121DEST_PATH_IMAGE063
is the voltage mode length of the photovoltaic node m,
Figure 771214DEST_PATH_IMAGE064
is the phase angle of the voltage at the photovoltaic node m after the (k-1) th iteration,
Figure 722990DEST_PATH_IMAGE065
for the voltage mode length of photovoltaic node m after the (k-1) th iteration,
Figure 726718DEST_PATH_IMAGE066
is the current mode length of the photovoltaic node m,
Figure 355014DEST_PATH_IMAGE067
for the mode length of the photovoltaic node m output current after the kth iteration,
Figure 811403DEST_PATH_IMAGE068
for photovoltaic node mInputting a current vector;
step C4: successive iteration calculation is carried out until the voltage difference value of each node is iterated for two times
Figure 566870DEST_PATH_IMAGE069
Is less than
Figure 425104DEST_PATH_IMAGE070
To obtain a warp
Figure 991346DEST_PATH_IMAGE071
Node i voltage vector after sub-iteration
Figure 935031DEST_PATH_IMAGE072
Can be based on
Figure 228609DEST_PATH_IMAGE073
Obtaining short-circuit current vectors of each branch circuit, wherein
Figure 941350DEST_PATH_IMAGE074
Is composed of
Figure 177029DEST_PATH_IMAGE071
The voltage vector at node j after the second iteration,
Figure 608010DEST_PATH_IMAGE075
is the impedance between node i and node j,
Figure 705279DEST_PATH_IMAGE076
is the vector of the injected current between node i and node j (as shown in fig. 2).
According to the method, the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage is obtained through calculation according to the voltage-controlled current source characteristic of the photovoltaic power supply
Figure 272526DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 180571DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 98848DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 734229DEST_PATH_IMAGE009
The short-circuit current of each branch can be obtained.
In summary, the method obtains the normal operation value of the voltage and the current of the photovoltaic node according to the load flow calculation before the fault; after the three-phase symmetrical grounding short circuit occurs to the power distribution network, respectively carrying out equivalent admittance processing suitable for Gaussian iterative short circuit calculation on the load nodes and the fault nodes to form a power distribution network admittance array containing a photovoltaic power supply; calculating and correcting the photovoltaic node injection current in the injection current vector according to the photovoltaic node voltage after each Gaussian iteration based on the voltage-controlled current source characteristic of the photovoltaic power supply; and after the short circuit calculation iterative convergence, the voltage value of each node after the fault is obtained, and the short circuit current of the whole network is obtained according to the voltage value. Therefore, the calculation scale of the short-circuit current calculation cannot be increased, original load flow calculation software can be well embedded, the method has good adaptability, the characteristics of the voltage-controlled current source of the new energy power supply are fully considered, the calculation accuracy can be improved when the short-circuit current calculation is carried out on the network containing the new energy power supply, the method has good reference value for accurately describing the fault characteristics of the power system containing the new energy power supply, and the setting and the verification of a relay protection principle are facilitated.
Referring to fig. 3, a block diagram of a short-circuit current calculation apparatus for a power distribution network including a photovoltaic power supply according to the present application is shown.
As shown in fig. 3, the distribution network short-circuit current calculation apparatus 200 includes a calculation module 210, a modification module 220, and an iteration module 230.
The calculation module 210 is configured to regard the photovoltaic node as a PQ node, have a power factor of 1, and obtain a node voltage and an injection current of each node including the photovoltaic node during normal operation by using load flow calculation, where a voltage vector of an ith node is
Figure 155983DEST_PATH_IMAGE001
The injection current vector of the i-th node is
Figure 484196DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 873458DEST_PATH_IMAGE003
(ii) a A modification module 220 configured to modify the load flow calculation node admittance array into a photovoltaic short-circuit current calculation admittance array by adopting an impedance simulation method suitable for Gaussian iterative photovoltaic short-circuit calculation after the three-phase symmetrical grounding short circuit of the system occurs and the fault point is regarded as a newly added node in the network
Figure 578109DEST_PATH_IMAGE004
Constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation; an iteration module 230 configured to obtain a photovoltaic node voltage vector according to the k-1 th iteration
Figure 119949DEST_PATH_IMAGE005
Calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage according to the voltage-controlled current source characteristic of the photovoltaic power supply
Figure 619063DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 997086DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 974269DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 636195DEST_PATH_IMAGE009
Thus, the short-circuit current of each branch can be obtained.
It should be understood that the modules depicted in fig. 3 correspond to various steps in the method described with reference to fig. 1. Thus, the operations and features described above for the method and the corresponding technical effects are also applicable to the modules in fig. 3, and are not described again here.
In other embodiments, an embodiment of the present invention further provides a computer-readable storage medium, where computer-executable instructions are stored, and the computer-executable instructions may execute the method for calculating a short-circuit current of a power distribution network including a photovoltaic power supply in any of the method embodiments described above;
as one embodiment, the computer-readable storage medium of the present invention stores computer-executable instructions configured to:
the photovoltaic node is regarded as a PQ node, the power factor is 1, the node voltage and the injection current of each node including the photovoltaic node during normal operation are obtained by load flow calculation, wherein the voltage vector of the ith node is PQ
Figure 306210DEST_PATH_IMAGE001
The injection current vector of the ith node is
Figure 670065DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 450939DEST_PATH_IMAGE003
；
After the system generates three-phase symmetrical grounding short circuit, the fault point is regarded as a newly added node in the network, and the load flow calculation node admittance array is modified into the photovoltaic short circuit current calculation admittance array by adopting an impedance simulation method suitable for the Gaussian iterative photovoltaic short circuit calculation
Figure 967371DEST_PATH_IMAGE004
Constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation;
obtaining a photovoltaic node voltage vector according to the k-1 iteration
Figure 808288DEST_PATH_IMAGE005
Calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage according to the voltage-controlled current source characteristic of the photovoltaic power supply
Figure 160903DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 745468DEST_PATH_IMAGE007
Photovoltaic node m of (1) inject current vector
Figure 116406DEST_PATH_IMAGE008
Successive iteration is carried out until the voltage variation of each node of the two previous iterations and the two subsequent iterations is less than
Figure 128225DEST_PATH_IMAGE009
The short-circuit current of each branch can be obtained.
The computer-readable storage medium may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created from use of a distribution network short-circuit current calculation device including a photovoltaic power supply, and the like. Further, the computer-readable storage medium may include high speed random access memory, and may also include memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid state storage device. In some embodiments, the computer readable storage medium optionally includes memory remotely located from the processor, and the remote memory may be connected to a distribution grid short circuit current calculation device containing a photovoltaic power source over a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 4, the electronic device includes: a processor 310 and a memory 320. The electronic device may further include: an input device 330 and an output device 340. The processor 310, the memory 320, the input device 330, and the output device 340 may be connected by a bus or other means, such as the bus connection in fig. 4. The memory 320 is the computer-readable storage medium described above. The processor 310 executes various functional applications and data processing of the server by running the nonvolatile software programs, instructions and modules stored in the memory 320, namely, the method for calculating the short-circuit current of the power distribution network including the photovoltaic power supply of the above method embodiment is realized. The input device 330 may receive input numeric or character information and generate key signal inputs related to user settings and function control of the distribution grid short circuit current calculation device including the photovoltaic power source. The output device 340 may include a display device such as a display screen.
The electronic device can execute the method provided by the embodiment of the invention, and has the corresponding functional modules and beneficial effects of the execution method. For technical details that are not described in detail in this embodiment, reference may be made to the method provided by the embodiment of the present invention.
As an embodiment, the electronic device is applied to a power distribution network short-circuit current calculation device including a photovoltaic power supply, and is used for a client, and the electronic device includes: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to:
the photovoltaic node is regarded as a PQ node, the power factor is 1, the node voltage and the injection current of each node including the photovoltaic node during normal operation are obtained by load flow calculation, wherein the voltage vector of the ith node is PQ node
Figure 201092DEST_PATH_IMAGE001
The injection current vector of the ith node is
Figure 589348DEST_PATH_IMAGE002
The injection power vector of the ith node is
Figure 814793DEST_PATH_IMAGE003
；
After the system generates three-phase symmetrical grounding short circuit, the fault point is regarded as a newly added node in the network, and the load flow calculation node admittance array is modified into the photovoltaic short circuit current calculation admittance array by adopting an impedance simulation method suitable for the Gaussian iterative photovoltaic short circuit calculation
Figure 997512DEST_PATH_IMAGE004
Constructing a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation;
obtaining a photovoltaic node voltage vector according to the k-1 iteration
Figure 324720DEST_PATH_IMAGE005
Calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage according to the voltage-controlled current source characteristic of the photovoltaic power supply
Figure 251087DEST_PATH_IMAGE006
Make the node after the k-1 iteration of replacement inject the current vector
Figure 331039DEST_PATH_IMAGE007
Photovoltaic node m injected current vector
Figure 950239DEST_PATH_IMAGE008
Successively iterating until the voltage variation of each node of the two iterations is less than
Figure 997698DEST_PATH_IMAGE009
The short-circuit current of each branch can be obtained.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods of the various embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, and not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (6)

1. A method for calculating short-circuit current of a power distribution network containing a photovoltaic power supply is characterized by comprising the following steps:
step A: the photovoltaic node is regarded as a PQ node, the power factor is 1, the node voltage and the injection current of each node including the photovoltaic node during normal operation are obtained by load flow calculation, wherein the voltage vector of the ith node is PQ
Figure FDA0003579914780000011
The injection current vector of the i-th node is
Figure FDA0003579914780000012
The injection power vector of the ith node is
Figure FDA0003579914780000013
And B: after a three-phase symmetrical grounding short circuit occurs in a system, a fault point is regarded as a node newly added in the network, a load flow calculation node admittance array is modified into a photovoltaic short circuit current calculation admittance array Y by adopting an impedance simulation method suitable for Gaussian iterative photovoltaic short circuit calculation, and a node impedance network equation suitable for Gaussian iterative photovoltaic short circuit calculation is constructed, wherein the construction of the photovoltaic short circuit calculation node impedance network equation suitable for Gaussian iteration comprises the following steps:
step B1: the node load simulation admittance obtained by using the node injection power and the node voltage under the normal network state is expressed as follows:
Figure FDA0003579914780000014
in the formula, YLoadiThe admittance is simulated for the load at the ith node,
Figure FDA0003579914780000015
being the conjugate of the injected power at the ith node,
Figure FDA0003579914780000016
is the voltage vector of the i-th node,
Figure FDA0003579914780000017
is the conjugate of the voltage of the ith node;
after equivalent load admittance, the node injects a current vector
Figure FDA0003579914780000018
Step B2: setting the three-phase grounding admittance at the fault point as yf；
Step B3: after the node analog impedance and the three-phase grounding admittance are determined, summing photovoltaic node self-admittance and photovoltaic node load analog admittance in the admittance array, and summing short-circuit fault node self-admittance and short-circuit fault node grounding admittance to obtain a photovoltaic short-circuit current admittance array expression:
Figure FDA0003579914780000021
in the formula, y11For node 1 self-admittance in load flow calculation, ynnFor the self-admittance of node n in the load flow calculation,
Figure FDA0003579914780000022
for the self-admittance of the photovoltaic node 1 in the power flow calculation,
Figure FDA0003579914780000023
for mutual admittance, y, of photovoltaic node m and node n in load flow calculationloadnFor load simulation admittance, y, of node ngIs the generator equivalent admittance;
step B4: the flow calculation before the fault obtains the vector of the current injected into the node 1 as
Figure FDA0003579914780000024
Node 1 voltage vector is
Figure FDA0003579914780000025
Wherein, the node 1 is a balance node, and the potential vector of the generator is calculated
Figure FDA0003579914780000026
The expression of (a) is:
Figure FDA0003579914780000027
in the formula, XgIs the equivalent reactance of the generator, and j is a mathematical imaginary number symbol;
generator potential vector before and after supposing short-circuit fault
Figure FDA0003579914780000028
Constant, equivalent reactance XgWithout change, the fault-time node 1 injection current vector can be expressed as:
Figure FDA0003579914780000029
wherein f (—) is a function symbol,
Figure FDA00035799147800000210
is the voltage vector of node 1;
a function of a photovoltaic node injection current vector and a photovoltaic node voltage vector, expressed as:
Figure FDA00035799147800000211
is a voltage vector of the photovoltaic node;
step B5: after the injection currents of the balance node and the photovoltaic node are determined, a node injection current vector is obtained
Figure FDA00035799147800000212
Expression:
Figure FDA0003579914780000031
step B6: the expression of the node impedance network equation is obtained as follows:
Figure FDA0003579914780000032
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003579914780000033
in the formula (I), the compound is shown in the specification,
Figure FDA0003579914780000034
is the injected current vector of the photovoltaic node 1,
Figure FDA0003579914780000035
is the vector of the injected current at photovoltaic node m,
Figure FDA0003579914780000036
is the voltage vector of the photovoltaic node 1,
Figure FDA0003579914780000037
is the voltage vector of the photovoltaic node m;
step C: obtaining a photovoltaic node voltage vector according to the k-1 iteration
Figure FDA0003579914780000038
According to the voltage-controlled current source characteristic of the photovoltaic power supply, calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage
Figure FDA0003579914780000041
Injecting current vector into node after k-1 iteration of replacement
Figure FDA0003579914780000042
Photovoltaic node m injected current vector
Figure FDA0003579914780000043
Successively iterating until the voltage variation of each node of the two iterations is less than epsilon 10-6The short-circuit current of each branch can be obtained.
2. The method as claimed in claim 1, wherein in step C, the output current vector of the pv inverter under d-axis voltage-oriented control is calculated according to the voltage-controlled current source characteristics of the pv power source
Figure FDA0003579914780000044
Injecting the current vector into the replacement node after the k-1 iteration
Figure FDA0003579914780000045
Photovoltaic node m injected current vector
Figure FDA0003579914780000046
The method comprises the following steps:
step C1: assuming photovoltaic power outputThe short-circuit voltage value of each node obtained by solving after the fault is set as the voltage vector iteration initial value of the Gaussian iteration short-circuit current node admittance equation
Figure FDA0003579914780000047
Wherein:
Figure FDA0003579914780000048
in the formula (I), the compound is shown in the specification,
Figure FDA0003579914780000049
iteration initial values are the n voltage vectors of the nodes;
step C2: the voltage iteration formula is as follows:
Figure FDA00035799147800000410
in the formula (I), the compound is shown in the specification,
Figure FDA00035799147800000411
the voltage at node i after the (k + 1) th iteration,
Figure FDA00035799147800000412
the voltage of a node I after the kth iteration is obtained, wherein I is an injection current, and Y is a photovoltaic short-circuit current admittance array;
step C3: before each iteration, according to the photovoltaic node voltage vector obtained after the k-1 iteration
Figure FDA00035799147800000413
Calculating the voltage value drop degree of the photovoltaic node
Figure FDA00035799147800000414
Figure FDA00035799147800000415
Is the voltage vector of the photovoltaic node m after the (k-1) th iteration,
Figure FDA00035799147800000416
for the voltage vector of the photovoltaic node m, according to the voltage-controlled source characteristic of the photovoltaic output current, when the voltage drop degree is more than or equal to 0.9, the photovoltaic power supply outputs active current with a normal operation value and does not provide reactive current, when the voltage drop degree is less than 0.9 and more than or equal to 0.2, the voltage drops by 1% every time, the photovoltaic power supply provides 1.5% of reactive current, before the current reaches the amplitude limit, the active current output is the normal operation value, after the current reaches the amplitude limit, the active current output is reduced, when the voltage drop degree is less than 0.2, the photovoltaic power supply does not provide active current, only provides reactive current with the amplitude of 1.05 times of the amplitude of the normal operation value, calculates the photovoltaic output current in the current iteration, replaces the node injection current vector after the k-1 iteration, and
Figure FDA00035799147800000417
photovoltaic node m injected current vector
Figure FDA00035799147800000418
Obtaining a node injection current vector after the kth iteration
Figure FDA00035799147800000419
Wherein the photovoltaic node m after the kth iteration injects a current vector
Figure FDA00035799147800000420
The mode length and the phase angle of (d) are respectively:
Figure FDA0003579914780000051
the phase angle of the photovoltaic output current is:
Figure FDA0003579914780000052
in the formula (I), the compound is shown in the specification,
Figure FDA0003579914780000053
is the voltage mode length of the photovoltaic node m,
Figure FDA0003579914780000054
is the phase angle of the voltage at the photovoltaic node m after the (k-1) th iteration,
Figure FDA0003579914780000055
the voltage modulus length of the photovoltaic node m after the (k-1) th iteration,
Figure FDA0003579914780000056
is the current mode length of the photovoltaic node m,
Figure FDA0003579914780000057
is the modular length of the photovoltaic node m output current after the kth iteration,
Figure FDA0003579914780000058
is the injection current vector of the photovoltaic node m.
3. The method for calculating the short-circuit current of the power distribution network with the photovoltaic power supply according to claim 1, wherein in the step C, the successive iteration is performed until the voltage variation of each node is less than 10 ∈ ═ in two iterations before and after the successive iteration-6The method for determining the short-circuit current of each branch circuit comprises the following steps:
step C4: successive iteration calculation is carried out until the voltage difference value of each node is iterated for two times
Figure FDA0003579914780000059
Less than 10 ∈ ═ 10-6Obtaining the voltage vector of the node i after k iterations
Figure FDA00035799147800000510
Can be based on
Figure FDA00035799147800000511
Obtaining the short-circuit current vector of each branch circuit, wherein
Figure FDA00035799147800000512
For the node j voltage vector, Z, after k iterationsijIs the impedance between node i and node j,
Figure FDA00035799147800000513
is the injected current vector between node i and node j.
4. The utility model provides a distribution network short-circuit current accounting device who contains photovoltaic power supply which characterized in that includes:
the calculation module is configured to regard the photovoltaic node as a PQ node, the power factor is 1, and the node voltage and the injection current of each node including the photovoltaic node during normal operation are obtained by utilizing load flow calculation, wherein the voltage vector of the ith node is PQ node
Figure FDA00035799147800000514
The injection current vector of the i-th node is
Figure FDA00035799147800000515
The injection power vector of the ith node is
Figure FDA0003579914780000061
The modification module is configured to modify a tidal current calculation node admittance array into a photovoltaic short-circuit current calculation admittance array Y by adopting an impedance simulation method suitable for Gaussian iterative photovoltaic short-circuit calculation after a three-phase symmetrical grounding short-circuit occurs in a system, and construct a node impedance network equation suitable for Gaussian iterative photovoltaic short-circuit calculation, wherein the construction of the photovoltaic short-circuit calculation node impedance network equation suitable for Gaussian iteration comprises the following steps:
the node load simulation admittance obtained by using the node injection power and the node voltage under the normal network state is expressed as follows:
Figure FDA0003579914780000062
in the formula, YLoadiThe admittance is simulated for the load at the ith node,
Figure FDA0003579914780000063
being the conjugate of the injected power at the ith node,
Figure FDA0003579914780000064
is the voltage vector of the i-th node,
Figure FDA0003579914780000065
is the conjugate of the voltage of the ith node;
after equivalent load admittance, the node injects current vector
Figure FDA0003579914780000066
Setting the three-phase grounding admittance at the fault point as yf；
After the node analog impedance and the three-phase grounding admittance are determined, summing photovoltaic node self-admittance and photovoltaic node load analog admittance in the admittance array, and summing short-circuit fault node self-admittance and short-circuit fault node grounding admittance to obtain a photovoltaic short-circuit current admittance array expression:
Figure FDA0003579914780000067
in the formula, y11For node 1 self-admittance in load flow calculation, ynnFor the self-admittance of node n in the load flow calculation,
Figure FDA0003579914780000068
for the self-admittance of the photovoltaic node 1 in the power flow calculation,
Figure FDA0003579914780000069
for mutual admittance, y, of photovoltaic node m and node n in load flow calculationloadnFor load simulation admittance, y, of node ngEquivalent admittance of the generator;
the flow calculation before the fault obtains the vector of the current injected into the node 1 as
Figure FDA0003579914780000071
Node 1 voltage vector is
Figure FDA0003579914780000072
Wherein, the node 1 is a balance node, and the potential vector of the generator is calculated
Figure FDA0003579914780000073
The expression of (a) is:
Figure FDA0003579914780000074
in the formula, XgJ is the equivalent reactance of the generator and is a mathematical imaginary number symbol;
generator potential vector before and after supposing short-circuit fault
Figure FDA0003579914780000075
Constant, equivalent reactance XgWithout change, the fault-time node 1 injection current vector can be expressed as:
Figure FDA0003579914780000076
wherein f (—) is a function symbol,
Figure FDA0003579914780000077
is the voltage vector of node 1;
a function of a photovoltaic node injection current vector and a photovoltaic node voltage vector, expressed as:
Figure FDA0003579914780000078
Figure FDA0003579914780000079
is a voltage vector of the photovoltaic node;
after the injection currents of the balance node and the photovoltaic node are determined, a node injection current vector is obtained
Figure FDA00035799147800000710
Expression:
Figure FDA00035799147800000711
the expression of the node impedance network equation is obtained as follows:
Figure FDA00035799147800000712
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003579914780000081
in the formula (I), the compound is shown in the specification,
Figure FDA0003579914780000082
is the injected current vector of the photovoltaic node 1,
Figure FDA0003579914780000083
is the vector of the injected current at photovoltaic node m,
Figure FDA0003579914780000084
is the voltage vector of the photovoltaic node 1,
Figure FDA0003579914780000085
is the voltage vector of the photovoltaic node m;
an iteration module configured to obtain a photovoltaic node voltage vector according to the k-1 th iteration
Figure FDA0003579914780000086
According to the voltage-controlled current source characteristic of the photovoltaic power supply, calculating to obtain the output current vector of the photovoltaic inverter under the directional control of the d-axis voltage
Figure FDA0003579914780000087
Injecting current vector into node after k-1 iteration of replacement
Figure FDA0003579914780000088
Photovoltaic node m of (1) inject current vector
Figure FDA0003579914780000089
The iteration is carried out successively until the voltage variation of each node of the two iterations is less than epsilon to 10-6The short-circuit current of each branch can be obtained.
5. An electronic device, comprising: at least one processor, and a memory communicatively coupled to the at least one processor, wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of any of claims 1-3.
6. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method of any one of claims 1 to 3.
CN202210144000.3A 2022-02-17 2022-02-17 Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply Active CN114188945B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210144000.3A CN114188945B (en) 2022-02-17 2022-02-17 Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210144000.3A CN114188945B (en) 2022-02-17 2022-02-17 Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply

Publications (2)

Publication Number Publication Date
CN114188945A CN114188945A (en) 2022-03-15
CN114188945B true CN114188945B (en) 2022-06-14

Family

ID=80546092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210144000.3A Active CN114188945B (en) 2022-02-17 2022-02-17 Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply

Country Status (1)

Country Link
CN (1) CN114188945B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115207920B (en) * 2022-09-13 2023-03-14 国网江西省电力有限公司电力科学研究院 Photovoltaic power supply-containing power distribution network short-circuit current calculation method and system based on factor table

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110199707A1 (en) * 2010-02-16 2011-08-18 Greenvolts, Inc Photovoltaic array ground fault detection method for utility-scale grounded solar electric power generating systems
CN106066918A (en) * 2016-06-06 2016-11-02 国网山东省电力公司经济技术研究院 Based on homotopy containing distributed power source and the short-circuit current calculation method of nonlinear-load
CN110336327A (en) * 2019-07-29 2019-10-15 国网上海市电力公司 Consider the power distribution network short circuit current acquisition methods of distributed photovoltaic low voltage crossing
CN111625914A (en) * 2020-05-25 2020-09-04 广东电网有限责任公司 Short-circuit current calculation method and system
CN113742907A (en) * 2021-08-24 2021-12-03 国网河南省电力公司电力科学研究院 Photovoltaic power station short-circuit current unified calculation method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110199707A1 (en) * 2010-02-16 2011-08-18 Greenvolts, Inc Photovoltaic array ground fault detection method for utility-scale grounded solar electric power generating systems
CN106066918A (en) * 2016-06-06 2016-11-02 国网山东省电力公司经济技术研究院 Based on homotopy containing distributed power source and the short-circuit current calculation method of nonlinear-load
CN110336327A (en) * 2019-07-29 2019-10-15 国网上海市电力公司 Consider the power distribution network short circuit current acquisition methods of distributed photovoltaic low voltage crossing
CN111625914A (en) * 2020-05-25 2020-09-04 广东电网有限责任公司 Short-circuit current calculation method and system
CN113742907A (en) * 2021-08-24 2021-12-03 国网河南省电力公司电力科学研究院 Photovoltaic power station short-circuit current unified calculation method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A_new_energy_efficiency_management_and_control_strategy_of_grid-friendly-based_intelligent_electricity;Wei Xie;《2014 International Conference on Power System Technology》;20141222;全文 *
Admittance_Extraction_of_Loaded_Grid_Power_Distribution_Network_with_Arbitrary_Shapes;SHI YAO;《2019 IEEE Asia-Pacific Microwave Conference (APMC)》;20200319;全文 *
含分布式电源的配电网潮流算法研究;刘冰;《中国优秀硕士学位论文全文数据库》;20190115;全文 *
适应逆变型分布式电源接入的配电网保护方法;潘本仁;《南方电网技术》;20181031;第12卷(第10期);全文 *

Also Published As

Publication number Publication date
CN114188945A (en) 2022-03-15

Similar Documents

Publication Publication Date Title
Ju et al. Loop‐analysis‐based continuation power flow algorithm for distribution networks
CVS et al. Enhanced modelling of doubly fed induction generator in load flow analysis of distribution systems
CN111797510A (en) Method and system for calculating short circuit ratio of new energy station
CN114204564B (en) Short-circuit current calculation method and device for power grid containing inversion type new energy
CN114520511A (en) Photovoltaic-considered distribution network three-phase asymmetric short-circuit current calculation method and device
CN114188945B (en) Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply
CN103956735A (en) Harmonic power flow analysis method of distributed power generation system
Chang et al. Data-driven estimation of voltage-to-power sensitivities considering their mutual dependency in medium voltage distribution networks
CN117498433A (en) Transient stability power limit quantization method for hybrid parallel system
Li et al. A loop-analysis theory based power flow method and its linear formulation for low-voltage DC grid
Fu et al. Exponential integration algorithm for large-scale wind farm simulation with Krylov subspace acceleration
CN112234598B (en) Electromagnetic transient simulation initialization method
CN114597902A (en) Power flow calculation method of power system based on convex relaxation
Bhimarasetti et al. A new contribution to distribution load flow analysis for radial and mesh distribution systems
Husain et al. Load flow analysis of radial and mesh distribution system using ZIP model
CN107425519B (en) Method for calculating maximum power supply capacity of three-phase power distribution network containing distributed power supply
Xu et al. Artificial neural network model of photovoltaic generator for power flow analysis in PSS® SINCAL
CN115207920B (en) Photovoltaic power supply-containing power distribution network short-circuit current calculation method and system based on factor table
CN107086603A (en) A kind of Random-fuzzy Continuation power flow of power system containing DFIG
CN114188944A (en) Method and device for calculating short-circuit current of power grid with distributed power supply
CN114091623A (en) Transient voltage partitioning method and device based on multi-binary-table-AP clustering
CN110046450B (en) Initial value selection method suitable for Newton method load flow calculation of superconducting cable-containing power grid
Babu et al. An efficient power flow method for distribution system studies under various load models
Guo et al. Fast optimal power allocation algorithm for multi‐terminal AC/DC hybrid grids with wind power integration
CN110224408B (en) Power system reduced interval Krawczyk iterative power flow method containing new energy power generation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant