CN112670987A - Power grid three-phase harmonic flow phasor matrix calculation method - Google Patents

Power grid three-phase harmonic flow phasor matrix calculation method Download PDF

Info

Publication number
CN112670987A
CN112670987A CN202011601327.6A CN202011601327A CN112670987A CN 112670987 A CN112670987 A CN 112670987A CN 202011601327 A CN202011601327 A CN 202011601327A CN 112670987 A CN112670987 A CN 112670987A
Authority
CN
China
Prior art keywords
harmonic
phase
power flow
matrix
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011601327.6A
Other languages
Chinese (zh)
Other versions
CN112670987B (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.)
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
State Grid Fujian Electric Power Co Ltd
Longyan Power Supply Co of State Grid Fujian Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
State Grid Fujian Electric Power Co Ltd
Longyan Power Supply Co of State Grid Fujian 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 Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd, State Grid Fujian Electric Power Co Ltd, Longyan Power Supply Co of State Grid Fujian Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Priority to CN202011601327.6A priority Critical patent/CN112670987B/en
Publication of CN112670987A publication Critical patent/CN112670987A/en
Application granted granted Critical
Publication of CN112670987B publication Critical patent/CN112670987B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a method for calculating a three-phase harmonic flow phasor matrix of a power grid, which comprises the following steps: setting harmonic times, and setting an initial value of a harmonic source statistical model and an algorithm iteration error limit; collecting parameters including transformer parameters, line parameters, node voltage, generator parameters and load parameters, and establishing a fundamental wave and harmonic wave three-phase admittance matrix; calculating fundamental wave three-phase power flow; detecting nodes with harmonic sources, predicting harmonic data by a harmonic source statistical model by using a fundamental wave three-phase power flow result, and calculating harmonic three-phase power flow; updating the fundamental wave three-phase load power by utilizing the harmonic wave three-phase network loss power; and judging whether the injection power error of the front node and the back node is smaller than a preset iteration error limit, if so, calculating convergence, outputting a harmonic three-phase power flow calculation result, and otherwise, returning to recalculate the fundamental three-phase power flow. The method carries out program architecture based on the thought of the phasor matrix, and realizes the rapid development and the efficient calculation of the harmonic power flow simulation of the large-scale three-phase power grid.

Description

Power grid three-phase harmonic flow phasor matrix calculation method
Technical Field
The invention relates to the technical field of power system safety, in particular to a power grid three-phase harmonic power flow phasor matrix calculation method.
Background
The harmonic problem of the power system has attracted people's attention in the 30 s of the 20 th century, and then due to the development of modern industrial technologies, the nonlinear load in the power grid has increased greatly, such as electric arc and contact welding equipment, submerged arc furnaces, ferrosilicon furnaces, high-frequency furnaces and the like which are widely used in industry, and meanwhile, due to the rapid development of power electronic technologies, thyristor rectification and current conversion technologies have been widely used, which leads to the substantial increase of the harmonic in the power grid. The harmonic waves have serious influence and harm on the power system, and are mainly manifested by reducing the power factor of the power supply system, overheating electrical equipment, insulating aging, generating vibration and noise, shortening the service life of the electrical equipment, causing parallel and series resonance of the power system to damage a capacitor and the like. In order to solve the harmonic hazard problem, the harmonic problem is researched by various input forces of industrial countries since the 70 s, and as an important content of the harmonic problem research, harmonic trend calculation is correspondingly developed.
The three-phase harmonic current calculation is an important content for researching harmonic problems, the three-phase harmonic current distribution of the electric network can be described through the three-phase harmonic current calculation, harmonic indexes of each node of the three phases of the electric network are obtained, the harmonic indexes are important basis for evaluating the safe operation of the electric power system, the reason for generating the harmonic waves can be analyzed according to the three-phase harmonic current calculation result, and harmonic treatment measures are further researched.
The current three-phase harmonic power flow calculation is not suitable for being applied to a large-scale three-phase network, and the harmonic source is not conveniently expanded, so that the efficiency of the program is not high, and the method is not suitable for the actual situation of the project.
Disclosure of Invention
In view of the above, the invention aims to provide a power grid three-phase harmonic power flow phasor matrix calculation method, which performs a program architecture based on the thought of a phasor matrix and realizes the rapid development and efficient calculation of large-scale three-phase power grid harmonic power flow simulation.
The invention is realized by adopting the following scheme: a power grid three-phase harmonic power flow phasor matrix calculation method specifically comprises the following steps:
setting harmonic times, and setting an initial value of a harmonic source statistical model and an algorithm iteration error limit;
collecting parameters including transformer parameters, line parameters, node voltage, generator parameters and load parameters, and establishing a fundamental wave and harmonic wave three-phase admittance matrix;
calculating fundamental wave three-phase power flow;
detecting nodes with harmonic sources, predicting harmonic data by a harmonic source statistical model by using a fundamental wave three-phase power flow result, and calculating harmonic three-phase power flow;
updating the fundamental wave three-phase load power by utilizing the harmonic wave three-phase network loss power;
and judging whether the injection power error of the front node and the back node is smaller than a preset iteration error limit, if so, calculating convergence, outputting a harmonic three-phase power flow calculation result, and otherwise, returning to recalculate the fundamental three-phase power flow.
Further, the setting of the harmonic times, the setting of the initial value of the harmonic source statistical model, and the algorithm iteration error limit specifically include: setting the harmonic times to be calculated, inputting the initial values of statistical models of different harmonic sources, and iterative convergence error values of fundamental wave and harmonic algorithm.
Further, the established fundamental wave three-phase admittance matrix
Figure BDA0002871496610000021
Comprises the following steps:
Figure BDA0002871496610000022
in the formula,
Figure BDA0002871496610000031
is an admittance matrix of the fundamental three-phase line,
Figure BDA0002871496610000032
the admittance matrix is a fundamental wave three-phase transformer;
wherein, the admittance matrix of the fundamental wave three-phase line
Figure BDA0002871496610000033
Comprises the following steps:
Figure BDA0002871496610000034
wherein nL represents the total number of lines, [ C ]Lf]And [ CLt]In order to be a correlation matrix, the correlation matrix,
Figure BDA0002871496610000035
Figure BDA0002871496610000036
the three-phase susceptance matrixes are respectively a t end and an f end of the line, and the upper mark i represents the ith line;
wherein,
Figure BDA0002871496610000037
Figure BDA0002871496610000038
Figure BDA0002871496610000039
in the formula, rL1,xL1Respectively representing the positive sequence resistance and the reactance of the line; r isL2,xL2Respectively representing the negative sequence resistance and the reactance of the line; r isL0,xL0Respectively representing zero sequence resistance and reactance of the circuit;
wherein, the admittance matrix of the fundamental wave three-phase transformer
Figure BDA00028714966100000310
Comprises the following steps:
Figure BDA0002871496610000041
wherein nT is the total number of transformers, [ CTf]And [ CTt]Respectively are incidence matrixes;
Figure BDA0002871496610000042
wherein f denotes the head end, t denotes the tail end, kf,ktRespectively representing the transformation ratio of an f-end transformer and the transformation ratio of a t-end transformer, and the superscript q represents the zero sequence resistance and the reactance of a q-th transformer;
Figure BDA0002871496610000043
wherein,
Figure BDA0002871496610000044
Figure BDA0002871496610000045
respectively the positive sequence resistance and reactance of the q-th transformer,
Figure BDA0002871496610000046
negative sequence resistance and reactance of q-th transformer respectively；
Figure BDA0002871496610000047
The zero sequence resistance and the reactance of the q-th transformer are respectively.
Further, the harmonic three-phase admittance matrix is established
Figure BDA0002871496610000048
Comprises the following steps:
Figure BDA0002871496610000049
in the formula,
Figure BDA00028714966100000410
is a harmonic admittance matrix of the transformer,
Figure BDA00028714966100000411
is a harmonic line admittance matrix, and is,
Figure BDA00028714966100000412
is a load harmonic admittance matrix, and the load harmonic admittance matrix,
Figure BDA00028714966100000413
is a generator harmonic admittance matrix.
Further, the calculating the fundamental wave three-phase power flow specifically includes: modifying a basic three-phase power flow equation into:
Figure BDA0002871496610000051
in the formula of UabcxIs the real part of three-phase voltage, UabcyIs the imaginary part of the three-phase voltage,
Figure BDA0002871496610000052
given a voltage, Δ P, for the PV nodeabcIs the three-phase active unbalance, Delta QabcThe three-phase reactive unbalance is obtained; both pvbus and pqbus are vector index matrixes from beginningInput node naming sequence determination; let the error correction amount be: [ Delta xabc]＝[J]\[F](ii) a J represents a Jacobian matrix, and F represents a power unbalance matrix;
the correction voltage is as follows:
Figure BDA0002871496610000053
in the formula i1Index vector i corresponding to real voltage part of pqbus2Index vector corresponding to imaginary voltage part of pqbus i3Index vector i corresponding to real voltage part of pvbus4The index vector corresponding to the voltage imaginary part of the pvbus is determined by the naming sequence of the initial input nodes; Δ xabcIs a voltage error correction matrix;
judging whether F is larger than a preset fundamental wave three-phase iteration error value, if so, circularly executing the calculation; if not, the three-phase fundamental wave power flow calculation is finished, and the three-phase fundamental wave power flow result is stored.
Further, the detecting the node with the harmonic source, predicting harmonic data by using a harmonic source statistical model according to a fundamental wave three-phase power flow result, and performing harmonic three-phase power flow calculation specifically comprises:
detecting nodes containing harmonic source marks by using the parameter marks, and judging whether a plurality of harmonic sources are in the same node; if so, considering the mutual influence among harmonic sources; if not, the mutual influence among harmonic sources is not considered;
each harmonic current:
Figure BDA0002871496610000068
if a plurality of harmonic sources are considered to be influenced by each other when the same node has the harmonic sources, the formula for considering the influence of the harmonic sources is obtained from the national standard as follows:
Figure BDA0002871496610000061
in the formula, harmonichIs the h-th harmonic current content, Iabcx、IabcyThe power is obtained by the fundamental current,
Figure BDA0002871496610000062
harmonic three-phase currents of the harmonic source 1 and the harmonic source 2 respectively; khIs a given coefficient;
by
Figure BDA0002871496610000063
Obtaining each harmonic three-phase voltage of each node, obtaining harmonic three-phase loss according to the harmonic three-phase voltage, and updating the three-phase power of the node by utilizing the harmonic loss of the harmonic source node; wherein,
Figure BDA0002871496610000064
is a three-phase harmonic voltage matrix and is provided with a plurality of harmonic voltage matrixes,
Figure BDA0002871496610000065
is a network harmonic admittance matrix, and the harmonic admittance matrix,
Figure BDA0002871496610000066
three-phase harmonic power.
Further, the step of judging whether the error of the injected power of the front node and the back node is smaller than a preset iteration error limit, if so, calculating convergence and outputting a harmonic three-phase power flow calculation result, otherwise, returning to recalculate the fundamental three-phase power flow specifically comprises the following steps:
Figure BDA0002871496610000067
in the formula,
Figure BDA0002871496610000071
for the three-phase power value matrix obtained for the k-th cycle,
Figure BDA0002871496610000072
obtaining a three-phase power value matrix for the k-1 th cycle; judgment of [ Delta Sabc(k)]Whether or not less than a preset harmonic iteration errorAnd if not, returning to recalculate the fundamental wave three-phase load flow, and if so, ending and outputting the harmonic wave three-phase load flow calculation result.
The invention also provides a power grid three-phase harmonic power flow phasor matrix calculation system, which comprises a memory, a processor and computer program instructions stored on the memory and capable of being executed on the processor, wherein when the processor executes the computer program instructions, the method steps as described above can be realized.
Compared with the prior art, the invention has the following beneficial effects: the method carries out program architecture based on the thought of the phasor matrix, and realizes the rapid development and the efficient calculation of the harmonic power flow simulation of the large-scale three-phase power grid.
Detailed Description
The present invention will be further described with reference to the following examples.
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
The embodiment provides a method for calculating a three-phase harmonic power flow phasor matrix of a power grid, which specifically comprises the following steps:
setting harmonic times, and setting an initial value of a harmonic source statistical model and an algorithm iteration error limit;
collecting parameters including transformer parameters, line parameters, node voltage, generator parameters and load parameters, and establishing a fundamental wave and harmonic wave three-phase admittance matrix;
calculating fundamental wave three-phase power flow;
detecting nodes with harmonic sources, predicting harmonic data by a harmonic source statistical model by using a fundamental wave three-phase power flow result, and calculating harmonic three-phase power flow;
updating the fundamental wave three-phase load power by utilizing the harmonic wave three-phase network loss power;
and judging whether the injection power error of the front node and the back node is smaller than a preset iteration error limit, if so, calculating convergence, outputting a harmonic three-phase power flow calculation result, and otherwise, returning to recalculate the fundamental three-phase power flow.
In this embodiment, the setting of the harmonic frequency, the setting of the initial value of the harmonic source statistical model, and the algorithm iteration error limit specifically include: setting the harmonic times to be calculated, inputting the initial values of statistical models of different harmonic sources, and iterative convergence error values of fundamental wave and harmonic algorithm.
In the present embodiment, a fundamental three-phase admittance matrix is established
Figure BDA0002871496610000081
Comprises the following steps:
Figure BDA0002871496610000082
in the formula,
Figure BDA0002871496610000083
is an admittance matrix of the fundamental three-phase line,
Figure BDA0002871496610000084
the admittance matrix is a fundamental wave three-phase transformer;
wherein, the admittance matrix of the fundamental wave three-phase line
Figure BDA0002871496610000085
Comprises the following steps:
Figure BDA0002871496610000086
wherein nL represents the total number of lines, [ C ]Lf]And [ CLt]In order to be a correlation matrix, the correlation matrix,
Figure BDA0002871496610000087
Figure BDA0002871496610000091
the three-phase susceptance matrixes are respectively a t end and an f end of the line, and the upper mark i represents the ith line;
wherein,
Figure BDA0002871496610000092
Figure BDA0002871496610000093
Figure BDA0002871496610000094
in the formula, rL1,xL1Respectively representing the positive sequence resistance and the reactance of the line; r isL2,xL2Respectively representing the negative sequence resistance and the reactance of the line; r isL0,xL0Respectively representing zero sequence resistance and reactance of the circuit;
wherein, the admittance matrix of the fundamental wave three-phase transformer
Figure BDA0002871496610000095
For (because the admittance model of the transformer is related to the wiring of the transformer, the wiring of Ynyn12 is only used as an example here):
Figure BDA0002871496610000096
wherein nT is the total number of transformers, [ CTf]And [ CTt]Respectively are incidence matrixes;
Figure BDA0002871496610000097
wherein f denotes the head end, t denotes the tail end, kf,ktRespectively representing the transformation ratio of an f-end transformer and the transformation ratio of a t-end transformer, and the superscript q represents the zero sequence resistance and the reactance of a q-th transformer;
Figure BDA0002871496610000101
wherein,
Figure BDA0002871496610000102
Figure BDA0002871496610000103
respectively the positive sequence resistance and reactance of the q-th transformer,
Figure BDA0002871496610000104
negative sequence resistance and reactance of the q-th transformer respectively;
Figure BDA0002871496610000105
the zero sequence resistance and the reactance of the q-th transformer are respectively.
Wherein,
Figure BDA0002871496610000106
in the formula,
Figure BDA0002871496610000107
is the zero-sequence impedance of the ith transformer,
Figure BDA0002871496610000108
is the negative sequence impedance of the ith transformer,
Figure BDA0002871496610000109
is the positive sequence impedance of the ith transformer,
Figure BDA00028714966100001010
for the zero sequence susceptance of the ith transformer,
Figure BDA00028714966100001011
is the negative sequence susceptance of the ith transformer,
Figure BDA00028714966100001012
the positive sequence susceptance is the ith transformer;
the incidence matrix is:
Figure BDA00028714966100001013
Figure BDA0002871496610000111
in the formula, Lf is line head end arrangement, nB is node number, Lt is line tail end arrangement, Tf is transformer head end arrangement, and Tt is transformer tail end arrangement;
in the present embodiment, a harmonic three-phase admittance matrix is established
Figure BDA0002871496610000112
Comprises the following steps:
Figure BDA0002871496610000113
in the formula,
Figure BDA0002871496610000114
is a harmonic admittance matrix of the transformer,
Figure BDA0002871496610000115
is a harmonic line admittance matrix, and is,
Figure BDA0002871496610000116
is a load harmonic admittance matrix, and the load harmonic admittance matrix,
Figure BDA0002871496610000117
is a generator harmonic admittance matrix.
Specifically, the generator harmonic impedance (neglecting the generator resistance) is:
Figure BDA0002871496610000118
e represents the e-th generator;
Figure BDA0002871496610000119
is the zero sequence reactance of the e-th generator,
Figure BDA00028714966100001110
is the negative-sequence reactance of the e-th generator,
Figure BDA00028714966100001111
is the positive sequence reactance of the e-th generator;
Figure BDA00028714966100001112
Figure BDA00028714966100001113
Figure BDA00028714966100001114
nB is the number of nodes of the system, and when a generator exists in a node, the h-th harmonic admittance corresponding to the node is
Figure BDA0002871496610000121
If the generator does not exist, the corresponding harmonic admittance is 0;
Figure BDA0002871496610000122
nh is the total number of harmonics to be calculated.
Wherein the load harmonic impedance is:
Figure BDA0002871496610000123
Figure BDA0002871496610000124
p represents the p-th load; sa,Sb,ScThe power of a phase, b phase and c phase is respectively given to the node, Pa,Pb,PcGiving the node active power of a phase, b phase and c phase, Qa,Qb,QcThe reactive power of a phase, b phase and c phase is given to the node,
Figure BDA0002871496610000125
the voltage of a phase of a p load node;
Figure BDA0002871496610000126
Figure BDA0002871496610000127
Figure BDA0002871496610000128
when the node has load, the h-th harmonic admittance corresponding to the node is
Figure BDA0002871496610000129
If no load exists, the corresponding harmonic admittance is 0;
Figure BDA00028714966100001210
wherein, the harmonic circuit is:
Figure BDA0002871496610000131
Figure BDA0002871496610000132
rL1,xL1respectively a line positive sequence resistor and a reactance; r isL2,xL2Respectively a line negative sequence resistance and a reactance; r isL0,xL0I represents the ith line for zero sequence resistance and reactance of the line respectively
Figure BDA0002871496610000133
Figure BDA0002871496610000134
Figure BDA0002871496610000135
Association matrix
Figure BDA0002871496610000136
Figure BDA0002871496610000137
Figure BDA0002871496610000138
Figure BDA0002871496610000139
Wherein, the transformer harmonic matrix is:
Figure BDA0002871496610000141
Figure BDA0002871496610000142
Figure BDA0002871496610000143
Figure BDA0002871496610000144
and (3) correlation matrix:
Figure BDA0002871496610000145
Figure BDA0002871496610000146
Figure BDA0002871496610000147
Figure BDA0002871496610000148
in summary, the harmonic total admittance matrix is:
Figure BDA0002871496610000149
in this embodiment, the calculating the fundamental three-phase power flow specifically includes:
voltage current cell phasor
Figure BDA0002871496610000151
Figure BDA0002871496610000152
The basic three-phase power flow equation is expressed as
Figure BDA0002871496610000153
Figure BDA0002871496610000154
Figure BDA0002871496610000155
Figure BDA0002871496610000156
Figure BDA0002871496610000157
[Iabcx]＝[Gabc][Uabcx]-[Babc][Uabcy]
Iabcy＝[Gabc][Uabcy]+[Babc][Uabcx]
Figure BDA0002871496610000158
Figure BDA0002871496610000159
Figure BDA0002871496610000161
Figure BDA0002871496610000162
Figure BDA0002871496610000163
Figure BDA0002871496610000164
[Iax],[Ibx],[Icx]The real part matrixes of the currents of the phase a, the phase b and the phase c are respectively [ Iay],[Iby],[Icy]The imaginary part matrixes of the currents of the phase a, the phase b and the phase c are respectively [ U ]ax],[Ubx],[Ucx]A phase a, a phase b and a phase c voltage real part matrix, [ U ]ay],[Uby],[Ucy]Respectively an imaginary part matrix of a phase, b phase and c phase voltage,
Figure BDA0002871496610000165
an active power matrix is injected for the node,
Figure BDA0002871496610000166
reactive power matrix is injected for the node, [ B ]abc]Is a susceptance matrix, [ G ]abc]In the form of a conductance matrix,
Figure BDA0002871496610000167
is a three-phase power matrix, [ U ]abcx]Is a matrix of real parts of three-phase voltages, [ I ]abcx]Is a matrix of real parts of three-phase current, [ U ]abcy]Is a three-phase voltage imaginary part matrix, [ I ]abcy]Is a three-phase current imaginary part matrix,
Figure BDA0002871496610000168
is a three-phase fundamental wave admittance matrix;
the voltage of the generator bus is controlled in general in the operation of the power grid, and the reactive power at the voltage-controlled bus is not controllable, so that the basic equation is revised as follows.
Modifying a basic three-phase power flow equation into:
Figure BDA0002871496610000171
in the formula of UabcxIs the real part of three-phase voltage, UabcyIs the imaginary part of the three-phase voltage,
Figure BDA0002871496610000172
given a voltage, Δ P, for the PV nodeabcIs the three-phase active unbalance, Delta QabcThe three-phase reactive unbalance is obtained; both pvbus and pqbus are vector index matrixes and are determined by the naming sequence of initial input nodes; let the error correction amount be: [ Delta xabc]＝[J]\[F](ii) a J represents a three-phase Jacobian matrix, and F represents a three-phase power unbalance matrix;
the correction voltage is as follows:
Figure BDA0002871496610000173
in the formula i1Index vector i corresponding to real voltage part of pqbus2Index vector corresponding to imaginary voltage part of pqbus i3Index vector i corresponding to real voltage part of pvbus4The index vector corresponding to the voltage imaginary part of the pvbus is determined by the naming sequence of the initial input nodes; Δ xabcThe error correction quantity of the three-phase voltage is obtained;
judging whether F is larger than a preset fundamental wave three-phase iteration error value, if so, circularly executing the calculation; if not, the three-phase fundamental wave power flow calculation is finished, and the three-phase fundamental wave power flow result is stored.
In this embodiment, the detecting a node where a harmonic source exists, predicting harmonic data by using a harmonic source statistical model according to a fundamental three-phase power flow result, and performing harmonic three-phase power flow calculation specifically includes:
detecting nodes containing harmonic source marks by using the parameter marks, and judging whether a plurality of harmonic sources are in the same node; if so, considering the mutual influence among harmonic sources; if not, the mutual influence among harmonic sources is not considered;
each harmonic current:
Figure BDA0002871496610000181
if a plurality of harmonic sources are considered to be influenced by each other when the same node has the harmonic sources, the formula for considering the influence of the harmonic sources is obtained from the national standard as follows:
Figure BDA0002871496610000182
in the formula, harmonichIs the h-th harmonic current content, Iabcx、IabcyThe power is obtained by the fundamental current,
Figure BDA0002871496610000183
harmonic three-phase currents of the harmonic source 1 and the harmonic source 2 respectively; khIs a given coefficient;
by
Figure BDA0002871496610000184
Obtaining each harmonic three-phase voltage of each node, obtaining harmonic three-phase loss according to the harmonic three-phase voltage, and updating the three-phase power of the node by utilizing the harmonic loss of the harmonic source node; wherein,
Figure BDA0002871496610000185
is a three-phase harmonic voltage matrix and is provided with a plurality of harmonic voltage matrixes,
Figure BDA0002871496610000186
is a three-phase harmonic admittance matrix,
Figure BDA0002871496610000187
is a three-phase harmonic power matrix.
In this embodiment, the determining whether the error of the injected power of the front and rear nodes is smaller than a preset iteration error limit, if yes, calculating convergence, and outputting a harmonic three-phase power flow calculation result, otherwise, returning to recalculate the fundamental three-phase power flow specifically includes:
Figure BDA0002871496610000188
in the formula,
Figure BDA0002871496610000189
for the three-phase power value matrix obtained for the k-th cycle,
Figure BDA00028714966100001810
obtaining a three-phase power value matrix for the k-1 th cycle; judgment of [ Delta Sabc(k)]And if not, returning to recalculate the fundamental wave three-phase power flow, and if so, ending and outputting a harmonic wave three-phase power flow calculation result.
The embodiment also provides a power grid three-phase harmonic power flow phasor matrix calculation system, which includes a memory, a processor and computer program instructions stored on the memory and capable of being executed on the processor, and when the computer program instructions are executed by the processor, the method steps as described above can be implemented.
The foregoing is directed to preferred embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow. However, any simple modification, equivalent change and modification of the above embodiments according to the technical essence of the present invention are within the protection scope of the technical solution of the present invention.

Claims (8)

1. A power grid three-phase harmonic power flow phasor matrix calculation method is characterized by comprising the following steps:
setting harmonic times, and setting an initial value of a harmonic source statistical model and an algorithm iteration error limit;
collecting parameters including transformer parameters, line parameters, node voltage, generator parameters and load parameters, and establishing a fundamental wave and harmonic wave three-phase admittance matrix;
calculating fundamental wave three-phase power flow;
detecting nodes with harmonic sources, predicting harmonic data by a harmonic source statistical model by using a fundamental wave three-phase power flow result, and calculating harmonic three-phase power flow;
updating the fundamental wave three-phase load power by utilizing the harmonic wave three-phase network loss power;
and judging whether the injection power error of the front node and the back node is smaller than a preset iteration error limit, if so, calculating convergence, outputting a harmonic three-phase power flow calculation result, and otherwise, returning to recalculate the fundamental three-phase power flow.
2. The method for calculating the phasor matrix of the three-phase harmonic power flow of the power grid according to claim 1, wherein the setting of the harmonic times, the setting of the initial value of the harmonic source statistical model and the algorithm iteration error limit are specifically as follows: setting the harmonic times to be calculated, inputting the initial values of statistical models of different harmonic sources, and iterative convergence error values of fundamental wave and harmonic algorithm.
3. The method for calculating the phasor matrix of the three-phase harmonic power flow of the power grid according to claim 1, wherein the established fundamental three-phase admittance matrix
Figure FDA0002871496600000011
Comprises the following steps:
Figure FDA0002871496600000012
in the formula,
Figure FDA0002871496600000013
is an admittance matrix of the fundamental three-phase line,
Figure FDA0002871496600000014
the admittance matrix is a fundamental wave three-phase transformer;
wherein, the admittance matrix of the fundamental wave three-phase line
Figure FDA0002871496600000015
Comprises the following steps:
Figure FDA0002871496600000021
wherein nL represents the total number of lines, [ C ]Lf]And [ CLt]In order to be a correlation matrix, the correlation matrix,
Figure FDA0002871496600000022
Figure FDA0002871496600000023
and
Figure FDA0002871496600000024
the three-phase susceptance matrixes are respectively a t end and an f end of the line, and the upper mark i represents the ith line;
wherein,
Figure FDA0002871496600000025
Figure FDA0002871496600000026
Figure FDA0002871496600000027
in the formula, rL1,xL1Respectively representing the positive sequence resistance and the reactance of the line; r isL2,xL2Respectively representing the negative sequence resistance and the reactance of the line; r isL0,xL0Respectively representing zero sequence resistance and reactance of the circuit;
wherein, the admittance matrix of the fundamental wave three-phase transformer
Figure FDA0002871496600000028
Comprises the following steps:
Figure FDA0002871496600000029
wherein nT is the total number of transformers, [ CTf]And [ CTt]Respectively are incidence matrixes;
Figure FDA0002871496600000031
wherein f denotes the head end, t denotes the tail end, kf,ktRespectively representing the transformation ratio of an f-end transformer and the transformation ratio of a t-end transformer, and the superscript q represents the zero sequence resistance and the reactance of a q-th transformer;
Figure FDA0002871496600000032
wherein,
Figure FDA0002871496600000033
Figure FDA0002871496600000034
respectively the positive sequence resistance and reactance of the q-th transformer,
Figure FDA0002871496600000035
negative sequence resistance and reactance of the q-th transformer respectively;
Figure FDA0002871496600000036
the zero sequence resistance and the reactance of the q-th transformer are respectively.
4. The method for calculating the phasor matrix of the three-phase harmonic power flow of the power grid as claimed in claim 1, wherein the established harmonic three-phase admittance matrix
Figure FDA0002871496600000037
Comprises the following steps:
Figure FDA0002871496600000038
in the formula,
Figure FDA0002871496600000039
is a harmonic admittance matrix of the transformer,
Figure FDA00028714966000000310
is a harmonic line admittance matrix, and is,
Figure FDA00028714966000000311
is a load harmonic admittance matrix, and the load harmonic admittance matrix,
Figure FDA00028714966000000312
is a generator harmonic admittance matrix.
5. The method for calculating the phasor matrix of the three-phase harmonic power flow of the power grid according to claim 1, wherein the calculating of the fundamental three-phase power flow specifically comprises: modifying a basic three-phase power flow equation into:
Figure FDA0002871496600000041
in the formula of UabcxIs the real part of three-phase voltage, UabcyIs the imaginary part of the three-phase voltage,
Figure FDA0002871496600000042
voltage set for PV node, Δ PabcIs the abc three-phase active unbalance quantity, delta QabcThe balance of abc three-phase reactive unbalance; both pvbus and pqbus are vector index matrixes and are determined by the naming sequence of initial input nodes; let the error correction amount be: [ Delta xabc]＝[J]\[F](ii) a J represents a Jacobian matrix, and F represents a power unbalance matrix;
the correction voltage is as follows:
Figure FDA0002871496600000043
in the formula i1Index vector i corresponding to real voltage part of pqbus2Index vector corresponding to imaginary voltage part of pqbus i3Index vector i corresponding to real voltage part of pvbus4The index vector corresponding to the voltage imaginary part of the pvbus is determined by the naming sequence of the initial input nodes; Δ xabcIs a voltage error correction matrix;
judging whether F is larger than a preset fundamental wave three-phase iteration error value, if so, circularly executing the calculation; if not, the three-phase fundamental wave power flow calculation is finished, and the three-phase fundamental wave power flow result is stored.
6. The method for calculating the phasor matrix of the three-phase harmonic power flow of the power grid according to claim 1, wherein the step of detecting the node with the harmonic source is implemented by predicting harmonic data by using a harmonic source statistical model by using a fundamental three-phase power flow result, and the step of performing harmonic three-phase power flow calculation specifically comprises the following steps:
detecting nodes containing harmonic source marks by using the parameter marks, and judging whether a plurality of harmonic sources are in the same node; if so, considering the mutual influence among harmonic sources; if not, the mutual influence among harmonic sources is not considered;
each harmonic current:
Figure FDA0002871496600000051
if a plurality of harmonic sources are considered to be influenced by each other when the same node has the harmonic sources, the formula for considering the influence of the harmonic sources is obtained from the national standard as follows:
Figure FDA0002871496600000052
in the formula, harmonichIs the h-th harmonic current content, Iabcx、IabcyThe real and imaginary parts of the current resulting from the fundamental power flow,
Figure FDA0002871496600000053
harmonic three-phase currents of the harmonic source 1 and the harmonic source 2 respectively; khIs a given coefficient;
by
Figure FDA0002871496600000054
Obtaining each harmonic three-phase voltage of each node, obtaining harmonic three-phase loss according to the harmonic three-phase voltage, and updating the three-phase power of the node by utilizing the harmonic loss of the harmonic source node; wherein,
Figure FDA0002871496600000055
is a three-phase harmonic voltage, and is,
Figure FDA0002871496600000056
is a network harmonic admittance matrix, and the harmonic admittance matrix,
Figure FDA0002871496600000057
three-phase harmonic power.
7. The method according to claim 1, wherein the determining whether the error of the injected power of the front and rear nodes is smaller than a preset iteration error limit is performed, if yes, convergence is performed, a harmonic three-phase power flow calculation result is output, and otherwise, the recalculating of the fundamental three-phase power flow is performed specifically:
Figure FDA0002871496600000061
in the formula,
Figure FDA0002871496600000062
for the three-phase power value matrix obtained for the k-th cycle,
Figure FDA0002871496600000063
obtaining a three-phase power value matrix for the k-1 th cycle; judgment of [ Delta Sabc(k)]And if not, returning to recalculate the fundamental wave three-phase power flow, and if so, ending and outputting a harmonic wave three-phase power flow calculation result.
8. A grid three-phase harmonic power flow phasor matrix calculation system comprising a memory, a processor and computer program instructions stored on the memory and executable on the processor, the computer program instructions when executed by the processor being operable to implement the method steps of any of claims 1 to 7.
CN202011601327.6A 2020-12-30 2020-12-30 Power grid three-phase harmonic current phasor matrix calculation method Active CN112670987B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011601327.6A CN112670987B (en) 2020-12-30 2020-12-30 Power grid three-phase harmonic current phasor matrix calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011601327.6A CN112670987B (en) 2020-12-30 2020-12-30 Power grid three-phase harmonic current phasor matrix calculation method

Publications (2)

Publication Number Publication Date
CN112670987A true CN112670987A (en) 2021-04-16
CN112670987B CN112670987B (en) 2023-06-23

Family

ID=75410613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011601327.6A Active CN112670987B (en) 2020-12-30 2020-12-30 Power grid three-phase harmonic current phasor matrix calculation method

Country Status (1)

Country Link
CN (1) CN112670987B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113690893A (en) * 2021-08-26 2021-11-23 广东电网有限责任公司广州供电局 Load flow calculation method and system of microgrid with voltage regulator
CN115036925A (en) * 2022-04-02 2022-09-09 深圳市宝安任达电器实业有限公司 Harmonic analysis admittance modeling method for nonlinear load of converter
CN113783211B (en) * 2021-10-12 2023-10-31 国网福建省电力有限公司龙岩供电公司 Power grid three-phase voltage unbalance analysis method based on three-phase tide technology

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102118026A (en) * 2011-03-21 2011-07-06 江苏省电力试验研究院有限公司 ETAP-based harmonic and negative sequence power flow method for large-scale power grid
CN106786608A (en) * 2017-03-22 2017-05-31 国网福建省电力有限公司 A kind of uncertain harmonic flow calculation method accessed suitable for distributed power source
CN109842118A (en) * 2019-01-30 2019-06-04 南方电网科学研究院有限责任公司 Power grid probability harmonic load flow calculation method considering correlation
CN110601198A (en) * 2019-10-30 2019-12-20 国网浙江省电力有限公司宁波供电公司 Hybrid micro-grid optimized operation method considering harmonic and voltage unbalance constraints

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102118026A (en) * 2011-03-21 2011-07-06 江苏省电力试验研究院有限公司 ETAP-based harmonic and negative sequence power flow method for large-scale power grid
CN106786608A (en) * 2017-03-22 2017-05-31 国网福建省电力有限公司 A kind of uncertain harmonic flow calculation method accessed suitable for distributed power source
CN109842118A (en) * 2019-01-30 2019-06-04 南方电网科学研究院有限责任公司 Power grid probability harmonic load flow calculation method considering correlation
CN110601198A (en) * 2019-10-30 2019-12-20 国网浙江省电力有限公司宁波供电公司 Hybrid micro-grid optimized operation method considering harmonic and voltage unbalance constraints

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
严;王金全;黄克峰;徐晔;: "计及谐波功率修正的独立微电网三相潮流计算" *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113690893A (en) * 2021-08-26 2021-11-23 广东电网有限责任公司广州供电局 Load flow calculation method and system of microgrid with voltage regulator
CN113690893B (en) * 2021-08-26 2022-08-02 广东电网有限责任公司广州供电局 Load flow calculation method and system of microgrid with voltage regulator
CN113783211B (en) * 2021-10-12 2023-10-31 国网福建省电力有限公司龙岩供电公司 Power grid three-phase voltage unbalance analysis method based on three-phase tide technology
CN115036925A (en) * 2022-04-02 2022-09-09 深圳市宝安任达电器实业有限公司 Harmonic analysis admittance modeling method for nonlinear load of converter

Also Published As

Publication number Publication date
CN112670987B (en) 2023-06-23

Similar Documents

Publication Publication Date Title
CN112670987B (en) Power grid three-phase harmonic current phasor matrix calculation method
Herman et al. A proportional-resonant current controller for selective harmonic compensation in a hybrid active power filter
Lei et al. An improved virtual resistance damping method for grid-connected inverters with LCL filters
CN103018534B (en) Determine the method and system of harmonic voltage
CN110943456B (en) Small signal harmonic state space modeling method and device for LCC-HVDC system
Yang et al. An overview of DC cable modelling for fault analysis of VSC-HVDC transmission systems
Luo et al. Analytical calculation of transient short-circuit currents for MMC-based MTDC grids
Lin et al. Multiple-frequency three-phase load flow for harmonic analysis
CN112667962B (en) Power grid harmonic flow phasor matrix calculation method
CN106532703B (en) Energy storage current transformation system, current transformer, control method of bidirectional current transformer and controller
CN109921422A (en) Active Power Filter-APF non-singular terminal sliding-mode control based on single Feedback Neural Network
CN110221168B (en) Method for positioning leading harmonic source and tracking harmonic pollution propagation path
CN111769592A (en) Virtual harmonic resistance control method of grid-connected inverter based on parabola method
Lenwari Optimized design of modified proportional-resonant controller for current control of active filters
Ranjbar et al. Comparison of compensation strategies for shunt active power filter control in unbalanced tree-phase four-wire systems
CN113642179A (en) Alternating-current impedance modeling method of modular multilevel converter
CN109462326B (en) The overcurrent protection Poewr control method of gird-connected inverter under a kind of unbalanced power grid
Selmi et al. A simple method for the steady state performances of self-excited induction generators
CN116148589B (en) Simplified analysis method and system for fault current of low-frequency power transmission system
Zhu et al. Virtual arm impedance emulation and stability improvement in modular multilevel converters
CN111769562A (en) Virtual impedance-based power quality frequency division control method and system
CN113783211B (en) Power grid three-phase voltage unbalance analysis method based on three-phase tide technology
Frater et al. Linearisation of nonlinear loads by phase dependent frequency coupling admittance matrices
Mayordomo et al. Decoupled Newton algorithms in the harmonic domain for the harmonic interaction of line commutated converters with AC systems
Godbole et al. An Efficient Harmonic Load Flow Solution with Reduced Jacobean Matrix

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