CN111864730B - Tracking and identifying method for equivalent parameters of power grid - Google Patents

Tracking and identifying method for equivalent parameters of power grid Download PDF

Info

Publication number
CN111864730B
CN111864730B CN201910363344.1A CN201910363344A CN111864730B CN 111864730 B CN111864730 B CN 111864730B CN 201910363344 A CN201910363344 A CN 201910363344A CN 111864730 B CN111864730 B CN 111864730B
Authority
CN
China
Prior art keywords
voltage
current
equivalent
order differential
differential equation
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
CN201910363344.1A
Other languages
Chinese (zh)
Other versions
CN111864730A (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.)
Chongqing University
State Grid Corp of China SGCC
State Grid Chongqing Electric Power Co Ltd
Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd
Chongqing University of Technology
Original Assignee
Chongqing University
State Grid Corp of China SGCC
State Grid Chongqing Electric Power Co Ltd
Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd
Chongqing University of Technology
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 Chongqing University, State Grid Corp of China SGCC, State Grid Chongqing Electric Power Co Ltd, Dongguan Power Supply Bureau of Guangdong Power Grid Co Ltd, Chongqing University of Technology filed Critical Chongqing University
Priority to CN201910363344.1A priority Critical patent/CN111864730B/en
Publication of CN111864730A publication Critical patent/CN111864730A/en
Application granted granted Critical
Publication of CN111864730B publication Critical patent/CN111864730B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a tracking and identifying method for equivalent parameters of a power grid, which comprises the following steps of: s1: extracting a voltage and current signal at a first moment after disturbance and a voltage and current signal at a second moment after disturbance from a port of the power network to be equalized, calculating voltage and current transient state electric quantity data, establishing a voltage and current high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current high-order differential equation of the whole power network to be equalized, and solving coefficients of the high-order differential equations; s2: and obtaining and simplifying the high-order differential equations of the voltage and current steady-state electric quantity at two sides of the port of the power network to be equivalent from each established high-order differential equation, enabling the high-order differential equations to be consistent with the forms of the power network equivalent simple network differential equations at two sides of the port, calculating power network equivalent simple network parameters according to the coefficients of the high-order differential equations, and realizing the tracking identification of the power network equivalent parameters. The method realizes tracking identification of equivalent parameters of the power grid by using the voltage, current and electric quantity in the running state of the power grid, can track the parameter characteristics of the running state of the power grid, and can provide an analysis basis for static stability and dynamic stability judgment of the power grid.

Description

Tracking and identifying method for equivalent parameters of power grid
Technical Field
The invention relates to the field of power grid analysis and operation, in particular to a power grid equivalent parameter tracking identification method.
Background
Real-time tracking and online analysis of the operation state of the power grid are necessary in aspects of safety and stability analysis of the power grid, safety weak link discovery and the like. With the continuous expansion of the interconnection scale of the power grid, the structure of the power grid becomes more and more complex, the real-time tracking and online analysis of the state of the power grid become more and more difficult, the online analysis and calculation of the existing complex power grid are usually completed on the basis of the simplification of the power grid, and a power grid equivalence method is a main means for the simplification of the power grid. The accuracy of the power grid equivalence depends on the effectiveness of a parameter identification method in the power grid equivalence model.
The existing power grid equivalent model parameter identification method generally identifies power grid equivalent model parameters by using steady-state electric quantity reflected by measurement data of a power grid in different running states, and the identification method has the contradiction between parameter identifiability and parameter change, namely, the parameters of the power grid equivalent model are changed along with factors such as power supply or power grid load change and power grid topological structure change, so that the change of the power grid equivalent model parameters needs to be tracked in real time in the parameter identification process. However, in a normal state, the change of the power grid operation state is smooth, so that the difference of the measured steady-state electric quantity data is not obvious, and the identification of model parameters cannot be realized; after the change of the operation condition of the power grid is enough to cause the obvious difference of the steady-state electric quantity data, the parameter to be identified of the model is changed, and the parameter identification cannot be realized by using the measured data on the premise of assuming that the parameter is not changed.
The method for realizing the tracking and identification of the equivalent parameters of the power grid based on the transient electric quantity generated by the disturbance of the power grid can solve the problems existing in the parameter identification based on the steady-state electric quantity, realizes the real-time tracking of the running state of the power grid, and is an important research direction of the power grid research. However, the existing power grid equivalent parameter identification based on the transient power has to collect power measurement data before and after the disturbance occurs, and the identification is established on the premise that the power grid parameters before and after the disturbance are not changed, so that the existing method is narrow in application range, and the transient power data is relatively complex to acquire.
Disclosure of Invention
The invention aims to provide a tracking and identifying method for equivalent parameters of a power grid, which realizes tracking and identifying of the equivalent parameters of the power grid by utilizing voltage, current and electric quantity in a power grid running state, and can perform multiple simple network equivalence on two sides of a port of a power network to be equivalent, so that the parameter characteristics of the power grid running state can be tracked, and an analysis basis can be provided for static stability and dynamic stability judgment of the power grid.
In order to solve the technical problem, the invention provides a power grid equivalent parameter tracking and identifying method, which comprises the following steps:
s1: extracting a voltage signal and a current signal at a first moment after disturbance and a voltage signal and a current signal at a second moment after disturbance from a port of a power network to be equalized, calculating voltage and current transient state electric quantity data according to the voltage signal and the current signal, establishing a voltage and current high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current high-order differential equation of the whole power network to be equalized, stripping a voltage and current transient state electric quantity part from each established high-order differential equation, obtaining a voltage and current transient state electric quantity high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current transient state electric quantity high-order differential equation of the whole power network to be equalized, and solving coefficients of each high-order differential equation according to the voltage and current transient state electric quantity data;
s2: stripping voltage and current steady state electric quantity parts from the established high-order differential equations to obtain voltage and current steady state electric quantity high-order differential equations at two sides of the port of the electric power network to be equalized, simplifying the voltage and current steady state electric quantity high-order differential equations at two sides of the port of the electric power network to be equalized to enable the voltage and current steady state electric quantity high-order differential equations to be consistent with the form of the electric network equivalent simple network differential equations at two sides of the port, associating the simplified steady state electric quantity high-order differential equations with corresponding coefficients of the electric network equivalent simple network differential equations, calculating electric network equivalent simple network parameters according to the coefficients of the high-order differential equations, and realizing tracking identification of the electric network equivalent parameters.
Further, the step S1 specifically includes:
s11: extracting a voltage signal u (t1) and a current signal i (t1) at the first time t1 and a voltage signal u (t2) and a current signal i (t2) at the second time t2 after a disturbance occurs from the grid; when the t2 is different from the t1 by one cycle, subtracting the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs from the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs; or when the difference between t2 and t1 is half a cycle, adding the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs to the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs to obtain a transient voltage signal Δ u (t) and a transient current signal Δ i (t);
s12: establishing a voltage-current high-order differential equation at two sides of a port of the power network to be equivalent:
Figure GDA0002159535640000031
Figure GDA0002159535640000032
and a voltage and current high-order differential equation of the whole power network to be equivalent:
Figure GDA0002159535640000033
wherein K is the order of the high order differential equation, u (t)(k)And i (t)(k)Is k-order derivative of port voltage signal and current signal, e (t) is equivalent power source of power network to be equivalent at right end of port,
Figure GDA0002159535640000034
an equivalent power supply of the power network to be equivalent at the left end of the port, akAnd bkThe high-order differential equation coefficient to be solved of the power network to be equivalent at the right end,
Figure GDA0002159535640000035
and
Figure GDA0002159535640000036
a high-order differential equation coefficient to be solved of the power network to be equivalent at the left end is represented by t, wherein K is 0, … and K;
here, the voltage and current signals of the port are composed of voltage and current steady-state signals and voltage and current transient signals, i.e.
Figure GDA0002159535640000037
Wherein u iss(t)、is(t) is a steady state voltage signal and a steady state current signal;
substituting the expression (4) into the expressions (1), (2) and (3) to obtain a voltage current steady-state electric quantity high-order differential equation:
Figure GDA0002159535640000038
Figure GDA0002159535640000039
Figure GDA00021595356400000310
and a voltage current transient state electric quantity high-order differential equation:
Figure GDA0002159535640000041
Figure GDA0002159535640000042
Figure GDA0002159535640000043
s13: according to the voltage and current transient state electric quantity data of S11, solving the coefficient a of the transient state electric quantity high-order differential equation of the formulas (8), (9) and (10) by adopting a least square methodk、bkAnd
Figure GDA0002159535640000044
further, step S2 specifically includes:
s21: will be provided with
Figure GDA0002159535640000045
Substituting the voltage current steady state electric quantity high-order differential equation at the right side of the port of the equation (5) to obtain:
Figure GDA0002159535640000046
the two sides of the formula (11) are simplified and combined to obtain
Figure GDA0002159535640000047
When K is an even number, let N be K/2, N be 0, …, N, and obtain:
Figure GDA0002159535640000048
Figure GDA0002159535640000051
when K is an odd number, let N be (K-1)/2, N be 0, …, N, and then:
Figure GDA0002159535640000052
Figure GDA0002159535640000053
s22: the power grid equivalent simple network differential equation of the power network to be equivalent on the right side of the port is set as follows:
Figure GDA0002159535640000054
wherein K 'is the order of the power grid equivalent simple network differential equation, obviously K'<K,us(t)(k)And is(t)(k)Is k order derivative of the port voltage steady state signal and the current steady state signal, eeq(t) is an equivalent power source a 'of the equivalent simple network of the right-end power grid'kAnd b'kThe coefficient of the equivalent simple network differential equation of the power grid is shown, t is time, and K is 0, … and K'; will be provided with
Figure GDA0002159535640000055
Substituting the equation (17) into the port right side power grid equivalent simple network differential equation to obtain:
Figure GDA0002159535640000056
by simply combining the two sides of formula (18), the product can be obtained
Figure GDA0002159535640000057
When K ' is an even number, let N ═ K '/2, N ═ 0, …, N ', we can find:
Figure GDA0002159535640000061
Figure GDA0002159535640000062
when K 'is an odd number, let N ═ 1)/2, N ═ 0, …, N', we can find:
Figure GDA0002159535640000063
Figure GDA0002159535640000064
s23: and (4) correlating the corresponding coefficients in the formula (12) in the step S21 and the formula (19) in the step S22, and calculating and solving equivalent simple network parameters according to the coefficients of the high-order differential equation to realize the tracking identification of the equivalent parameters of the power grid on the right side of the port.
S24: will be provided with
Figure GDA0002159535640000065
And substituting the voltage and current steady-state electric quantity high-order differential equation into the port left side voltage and current steady-state electric quantity high-order differential equation in the formula (6), and adopting the same methods of S21, S22 and S23 to realize the tracking identification of the equivalent parameters of the port left side power grid.
The invention has the beneficial effects that: the invention realizes the tracking identification of the equivalent parameters of the power grid by utilizing the voltage, current and electric quantity in the running state of the power grid, and can carry out multiple simple network equivalence on two sides of the port of the power network to be equivalent, thereby tracking the parameter characteristics of the running state of the power grid and providing an analysis basis for the static stability and dynamic stability judgment of the power grid. The method is completely based on the voltage, current and electric quantity of a power grid in the running state, so that data measurement is simplified, and the additional condition that the equivalent parameters are not changed before and after disturbance does not need to be observed.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. In the drawings:
FIG. 1 is a diagram of an electric power network with waiting equivalence
FIG. 2 is a circuit diagram of a Thevenin equivalent model according to an embodiment of the present invention;
FIG. 3 is a circuit diagram of a Norton equivalent model according to an embodiment of the present invention;
FIG. 4 is a circuit diagram of a generic equivalent model according to an embodiment of the present invention.
Detailed Description
A power grid equivalent parameter tracking and identifying method comprises the following steps:
s1: extracting a voltage signal and a current signal at a first moment after disturbance and a voltage signal and a current signal at a second moment after disturbance from a port of the power network to be equalized, which are shown in FIG. 1, calculating voltage and current transient state electric quantity data according to the voltage signal and the current signal, establishing a voltage and current high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current high-order differential equation of the whole power network to be equalized, stripping a voltage and current transient state electric quantity part from each established high-order differential equation, obtaining the voltage and current transient state electric quantity high-order differential equation at two sides of the port of the power network to be equalized and the voltage and current transient state electric quantity high-order differential equation of the whole power network to be equalized, and solving coefficients of each high-order differential equation according to the voltage and current transient state electric quantity data;
s2: stripping voltage and current steady state electric quantity parts from the established high-order differential equations to obtain voltage and current steady state electric quantity high-order differential equations at two sides of the port of the electric power network to be equalized, simplifying the voltage and current steady state electric quantity high-order differential equations at two sides of the port of the electric power network to be equalized to enable the voltage and current steady state electric quantity high-order differential equations to be consistent with the form of the electric network equivalent simple network differential equations at two sides of the port, associating the simplified steady state electric quantity high-order differential equations with corresponding coefficients of the electric network equivalent simple network differential equations, calculating electric network equivalent simple network parameters according to the coefficients of the high-order differential equations, and realizing tracking identification of the electric network equivalent parameters.
Wherein, the step S1 specifically includes:
s11: for the power network to be equalized shown in fig. 1, the power network to be equalized includes power networks to be equalized on the left and right sides of a port; extracting from the port the voltage signal u (t1) and the current signal i (t1) at said first instant t1 and the voltage signal u (t2) and the current signal i (t2) at said second instant t2 after the occurrence of the disturbance; when the t2 is different from the t1 by one cycle, subtracting the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs from the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs; or when the difference between t2 and t1 is half a cycle, adding the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs to the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs to obtain a transient voltage signal Δ u (t) and a transient current signal Δ i (t);
s12: establishing a voltage-current high-order differential equation at two sides of a port of the power network to be equivalent:
Figure GDA0002159535640000081
Figure GDA0002159535640000082
and a voltage and current high-order differential equation of the whole power network to be equivalent:
Figure GDA0002159535640000083
wherein K is the order of the high order differential equation, u (t)(k)And i (t)(k)Is k-order derivative of port voltage signal and current signal, e (t) is equivalent power source of right-end equivalent power network,
Figure GDA0002159535640000084
for the equivalent power supply of the left-end standby equivalent power network, akAnd bkThe high-order differential equation coefficient to be solved of the power network to be equivalent at the right end,
Figure GDA0002159535640000085
and
Figure GDA0002159535640000086
a high-order differential equation coefficient to be solved of the power network to be equivalent at the left end is represented by t, wherein K is 0, … and K;
here, the voltage and current signals of the port are composed of voltage and current steady-state signals and voltage and current transient signals, i.e.
Figure GDA0002159535640000087
Wherein u iss(t)、is(t) is stableA steady state voltage signal and a steady state current signal.
Substituting the expression (4) into the expressions (1), (2) and (3) to obtain a voltage current steady-state electric quantity high-order differential equation:
Figure GDA0002159535640000091
Figure GDA0002159535640000092
Figure GDA0002159535640000093
and a voltage current transient state electric quantity high-order differential equation:
Figure GDA0002159535640000094
Figure GDA0002159535640000095
Figure GDA0002159535640000096
s13: according to the voltage and current transient state electric quantity data of S11, solving the coefficient a of the transient state electric quantity high-order differential equation of the formulas (8), (9) and (10) by adopting a least square methodk、bkAnd
Figure GDA0002159535640000097
the step S2 specifically includes:
s21: will be provided with
Figure GDA0002159535640000098
Substituting the voltage and current steady state electric quantity high-order differential equation on the right side of the port in the formula (5) to obtain:
Figure GDA0002159535640000099
The two sides of the formula (11) are simplified and combined to obtain
Figure GDA00021595356400000910
When K is an even number, let N be K/2, N be 0, …, N, and obtain:
Figure GDA0002159535640000101
Figure GDA0002159535640000102
when K is an odd number, let N be (K-1)/2, N be 0, …, N, and then:
Figure GDA0002159535640000103
Figure GDA0002159535640000104
s22: the equivalent simple network differential equation of the power grid on the right side of the port is set as follows:
Figure GDA0002159535640000105
wherein K 'is the order of the power grid equivalent simple network differential equation, obviously K'<K,us(t)(k)And is(t)(k)Is k order derivative of the port voltage steady state signal and the current steady state signal, eeq(t) is an equivalent power source a 'of the equivalent simple network of the right-end power grid'kAnd b'kThe coefficient of the equivalent simple network differential equation of the power grid is shown, t is time, and K is 0, … and K'; will be provided with
Figure GDA0002159535640000106
Substituting the equation (17) into the port right side power grid equivalent simple network differential equation to obtain:
Figure GDA0002159535640000107
by simply combining the two sides of formula (18), the product can be obtained
Figure GDA0002159535640000111
When K ' is an even number, let N ═ K '/2, N ═ 0, …, N ', we can find:
Figure GDA0002159535640000112
Figure GDA0002159535640000113
when K 'is an odd number, let N ═ 1)/2, N ═ 0, …, N', we can find:
Figure GDA0002159535640000114
Figure GDA0002159535640000115
s23: correlating the corresponding coefficients in the formula (12) in the step S21 and the formula (19) in the step S22, calculating and solving equivalent simple network parameters according to the coefficients of the high-order differential equation, and realizing the tracking identification of the equivalent parameters of the power grid on the right side of the port; if the number of the power grid equivalent simple network parameters is less than 4, the power grid equivalent simple network parameters can be uniquely solved, and if the number of the power grid equivalent simple network parameters is more than 4, the power grid equivalent simple network parameters have multiple groups of solutions.
S24: will be provided with
Figure GDA0002159535640000116
Substituting the voltage and current steady-state electric quantity high-order differential equation of the left side of the port in the formula (6), and adopting the same methods of S21, S22 and S23 to realize the tracking identification of equivalent parameters of the left side of the port; if the number of the power grid equivalent simple network parameters is less than 4, the power grid equivalent simple network parameters can be uniquely solved, and if the number of the power grid equivalent simple network parameters is more than 4, the power grid equivalent simple network parameters have multiple groups of solutions.
Specific applications of the above-described method are specifically exemplified below.
Example 1
The power network shown in fig. 2 is used as a power network to be equalized, a voltage-current transient state electric quantity high-order differential equation of the power network to be equalized is obtained through the method step S1, and each order derivative Δ u (t) of the transient state electric quantity is obtained through a numerical method(k)、Δi(t)(k)Solving the coefficient a of the higher order differential equation by using a least square methodk、bk(ii) a The high-order differential equation of the voltage and current steady-state electric quantity at two sides of the power network port to be equalized is obtained through the method step S2, the high-order differential equation of the voltage and current steady-state electric quantity at two sides of the power network port to be equalized is simplified, and four parameters A under the steady-state condition can be obtained through the formula (12) in S210、A1、B0、B1
From thevenin's theorem, for any port power network, the voltage-current characteristics of the power network port can be characterized in the form of a power source series impedance (as shown in fig. 2). The Thevenin equivalent network differential equation under the steady state is obtained from FIG. 2:
Figure GDA0002159535640000121
will be provided with
Figure GDA0002159535640000122
Substituting formula (24) to obtain:
Figure GDA0002159535640000123
in comparison with formula (19), it is found that: a'0=1、A′1=0、B′0=Req、B′1=ωLeq. The formula (12) is converted into a form consistent with the equivalent simple network differential equation (25) of the power grid by adopting mathematical transformation to obtain
Figure GDA0002159535640000124
And then, correlating the simplified high-order differential equation (26) with a corresponding coefficient of the power grid equivalent simple network differential equation (25), and solving equivalent simple network parameters to obtain Thevenin equivalent parameters:
Figure GDA0002159535640000131
according to the voltage and current values measured by the power network port, the Thevenin equivalent power supply can be obtained by using the formula (25)
Figure GDA0002159535640000132
And completing the tracking identification of the Thevenin equivalent network parameters of the power grid.
Example 2
As shown in fig. 3, the power network to be equalized is a power network to be equalized, and as known from norton's theorem, for any port power network, it can be equivalently modeled as a current source parallel reactance (as shown in fig. 3), and a norton's equivalent network differential equation in a steady state is obtained from fig. 3:
Figure GDA0002159535640000133
will be provided with
Figure GDA0002159535640000134
Substituting formula (24) to obtain:
Figure GDA0002159535640000135
in comparison with formula (19), it is found that: a'0=1/Req、A′1=ωCeq、B′0=1、B′10. The formula (12) is converted into a form consistent with the equivalent simple network differential equation (29) of the power grid by adopting mathematical transformation, and the result is obtained
Figure GDA0002159535640000136
Then, the simplified high-order differential equation (30) is associated with the corresponding coefficient of the equivalent simple network differential equation (29) of the power grid, and the equivalent simple network parameters are solved to obtain
Figure GDA0002159535640000141
According to the voltage and current values measured by the power network port, the Noton equivalent power supply can be obtained by using the formula (29)
Figure GDA0002159535640000142
And completing the tracking identification of the network parameters of the power grid Noton equivalent.
Example 3
As shown in fig. 4, the power network is used as a power network to be equivalent, in an actual power grid equivalent method, a complex system is further equivalent to a model (as shown in fig. 4) in which a resistance series inductance is added with a capacitance to ground, and an equivalent network differential equation in a steady state of the embodiment 3 is obtained from fig. 4:
Figure GDA0002159535640000143
will be provided with
Figure GDA0002159535640000144
Substituting formula (32) to obtain:
Figure GDA0002159535640000145
correlating the simplified differential equation (12) with the corresponding coefficient of the equivalent simple network differential equation (33) of the power grid, and solving the equivalent simple network parameters to obtain
Figure GDA0002159535640000146
The equivalent power source of example 3 can be obtained from the voltage and current values measured at the power network port by using equation (33)
Figure GDA0002159535640000147
And completing the tracking identification of the equivalent network parameters in example 3.
The method can theoretically solve the equivalent parameters of different equivalent models under the condition that the number of the to-be-equivalent power network parameters is less than or equal to 4. The method can track and identify parameters aiming at various equivalent power network models.
Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all of them should be covered in the claims of the present invention.

Claims (1)

1. A power grid equivalent parameter tracking and identifying method is characterized by comprising the following steps:
s1: extracting a voltage signal and a current signal at a first moment after disturbance and a voltage signal and a current signal at a second moment after disturbance from a port of a power network to be equalized, calculating voltage and current transient state electric quantity data according to the voltage signal and the current signal, establishing a voltage and current high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current high-order differential equation of the whole power network to be equalized, stripping a voltage and current transient state electric quantity part from each established high-order differential equation, obtaining a voltage and current transient state electric quantity high-order differential equation at two sides of the port of the power network to be equalized and a voltage and current transient state electric quantity high-order differential equation of the whole power network to be equalized, and solving coefficients of each high-order differential equation according to the voltage and current transient state electric quantity data;
s2: stripping a voltage and current steady state electric quantity part from each established high-order differential equation to obtain voltage and current steady state electric quantity high-order differential equations at two sides of a port of the electric power network to be equivalent, simplifying the voltage and current steady state electric quantity high-order differential equations at two sides of the port of the electric power network to be equivalent so as to be consistent with the form of the electric network equivalent simple network differential equations at two sides of the port, associating the simplified steady state electric quantity high-order differential equations with corresponding coefficients of the electric network equivalent simple network differential equations, calculating electric network equivalent simple network parameters according to the coefficients of the high-order differential equations, and realizing tracking identification of the electric network equivalent parameters;
wherein, the step S1 specifically includes:
s11: extracting a voltage signal u (t1) and a current signal i (t1) at the first time t1 and a voltage signal u (t2) and a current signal i (t2) at the second time t2 after a disturbance occurs from the grid; when the t2 is different from the t1 by one cycle, subtracting the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs from the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs; or when the difference between t2 and t1 is half a cycle, adding the voltage signal u (t2) and the current signal i (t2) after the disturbance occurs to the voltage signal u (t1) and the current signal i (t1) after the disturbance occurs to obtain a transient voltage signal Δ u (t) and a transient current signal Δ i (t);
s12: establishing a voltage-current high-order differential equation at two sides of a port of the power network to be equivalent:
Figure FDA0003549127750000011
Figure FDA0003549127750000021
and a voltage and current high-order differential equation of the whole power network to be equivalent:
Figure FDA0003549127750000022
wherein K is the order of the high order differential equation, u (t)(k)And i (t)(k)Is k-order derivative of port voltage signal and current signal, e (t) is equivalent power source of power network to be equivalent at right end of port,
Figure FDA0003549127750000023
an equivalent power supply of the power network to be equivalent at the left end of the port, akAnd bkThe high-order differential equation coefficient to be solved of the power network to be equivalent at the right end,
Figure FDA0003549127750000024
and
Figure FDA0003549127750000025
a high-order differential equation coefficient to be solved of the power network to be equivalent at the left end is represented by t, wherein K is 0, … and K;
here, the voltage and current signals of the port are composed of voltage and current steady-state signals and voltage and current transient signals, i.e.
Figure FDA0003549127750000026
Wherein u iss(t)、is(t) is a steady state voltage signal and a steady state current signal;
substituting the expression (4) into the expressions (1), (2) and (3) to obtain a voltage current steady-state electric quantity high-order differential equation:
Figure FDA0003549127750000027
Figure FDA0003549127750000028
Figure FDA0003549127750000029
and a voltage current transient state electric quantity high-order differential equation:
Figure FDA00035491277500000210
Figure FDA00035491277500000211
Figure FDA00035491277500000212
s13: according to the voltage and current transient state electric quantity data of S11, solving the coefficient a of the transient state electric quantity high-order differential equation of the formulas (8), (9) and (10) by adopting a least square methodk、bkAnd
Figure FDA0003549127750000031
further, step S2 specifically includes:
s21: will be provided with
Figure FDA0003549127750000032
Substituting into the port of formula (5)And (3) obtaining a high-order differential equation of the steady-state electric quantity of the voltage and the current on the right side:
Figure FDA0003549127750000033
the two sides of the formula (11) are simplified and combined to obtain
Figure FDA0003549127750000034
When K is an even number, let N be K/2, N be 0, …, N, and obtain:
Figure FDA0003549127750000035
Figure FDA0003549127750000036
when K is an odd number, let N be (K-1)/2, N be 0, …, N, and then:
Figure FDA0003549127750000037
Figure FDA0003549127750000041
s22: the power grid equivalent simple network differential equation of the power network to be equivalent on the right side of the port is set as follows:
Figure FDA0003549127750000042
wherein K 'is the order of the power grid equivalent simple network differential equation, obviously K'<K,us(t)(k)And is(t)(k)Is k order derivative of the port voltage steady state signal and the current steady state signal, eeq(t) is an equivalent power source a 'of the equivalent simple network of the right-end power grid'kAnd b'kThe coefficient of the equivalent simple network differential equation of the power grid is shown, t is time, and K is 0, … and K'; will be provided with
Figure FDA0003549127750000043
Substituting the equation (17) into the port right side power grid equivalent simple network differential equation to obtain:
Figure FDA0003549127750000044
by simply combining the two sides of formula (18), the product can be obtained
Figure FDA0003549127750000045
When K ' is an even number, let N ═ K '/2, N ═ 0, …, N ', we can find:
Figure FDA0003549127750000046
Figure FDA0003549127750000047
when K 'is an odd number, let N ═ 1)/2, N ═ 0, …, N', we can find:
Figure FDA0003549127750000051
Figure FDA0003549127750000052
s23: correlating the corresponding coefficients in the formula (12) in the step S21 and the formula (19) in the step S22, calculating and solving equivalent simple network parameters according to the coefficients of the high-order differential equation, and realizing the tracking identification of the equivalent parameters of the power grid on the right side of the port;
s24: will be provided with
Figure FDA0003549127750000053
And substituting the voltage and current steady-state electric quantity high-order differential equation into the port left side voltage and current steady-state electric quantity high-order differential equation in the formula (6), and adopting the same methods of S21, S22 and S23 to realize the tracking identification of the equivalent parameters of the port left side power grid.
CN201910363344.1A 2019-04-30 2019-04-30 Tracking and identifying method for equivalent parameters of power grid Active CN111864730B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910363344.1A CN111864730B (en) 2019-04-30 2019-04-30 Tracking and identifying method for equivalent parameters of power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910363344.1A CN111864730B (en) 2019-04-30 2019-04-30 Tracking and identifying method for equivalent parameters of power grid

Publications (2)

Publication Number Publication Date
CN111864730A CN111864730A (en) 2020-10-30
CN111864730B true CN111864730B (en) 2022-05-03

Family

ID=72966685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910363344.1A Active CN111864730B (en) 2019-04-30 2019-04-30 Tracking and identifying method for equivalent parameters of power grid

Country Status (1)

Country Link
CN (1) CN111864730B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104391198B (en) * 2014-11-26 2017-06-23 国网四川省电力公司自贡供电公司 A kind of low pressure power network safety monitoring method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104391198B (en) * 2014-11-26 2017-06-23 国网四川省电力公司自贡供电公司 A kind of low pressure power network safety monitoring method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
低压配电网三相不平衡负载参数的辨识方法研究;张鹏等;《电力***保护与控制》;20151216;第43卷(第24期);全文 *
基于暂态电压电流故障分量的电网等值模型参数辨识方法研究;王学斌;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20141215;第1-84页 *

Also Published As

Publication number Publication date
CN111864730A (en) 2020-10-30

Similar Documents

Publication Publication Date Title
CN108667048B (en) Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system
CN102567603B (en) Method for automatically generating BPA calculation file based on actual measurement topology and measured data
JP2013208051A (en) Method and system for determining state of bus of power system
CN106443246B (en) The on-line identification method of small interference stability parameter based on PMU metric data
CN107577870A (en) The distribution network voltage power sensitivity robust estimation method measured based on synchronized phasor
CN105375484A (en) PMU-based electric power system distributed dynamic-state estimation method
CN103928925A (en) Power distribution network load flow calculation method based on forward-backward sweep
CN111141995B (en) Line double-end steady-state distance measuring method and system based on amplitude comparison principle
Vignesh et al. An experimental study on the load modelling using PMU measurements
Shahraeini et al. A survey on topological observability of power systems
CN103793556B (en) Comprehensively go out transformer station&#39;s load modeling method of specific electric load
CN109358268A (en) A kind of accurate fault distance-finding method of double-circuit line on same pole independent of line parameter circuit value
CN109494726A (en) Stability of power system online evaluation method based on DLRNN neural network
CN108614194A (en) A kind of double line down distance measuring method and device based on distributed parameter model
CN105160058B (en) A kind of network boundary equivalence method calculated for electric network electromagnet transient
CN111864730B (en) Tracking and identifying method for equivalent parameters of power grid
CN103823998A (en) Weak cross section determination method taking influence of network topology changes on power transmission capacity into consideration
CN109270405A (en) A kind of Zero sequence parameter calculation method and system based on double-circuit line characteristic equation
CN109858061A (en) Distribution network equivalence and simplifying method towards voltage power sensitivity estimation
CN117310369A (en) Fault section identification method and system based on main feeder node voltage difference matrix
CN104240036B (en) A kind of Transient Voltage Stability quantitative estimation method based on critical system equiva lent impedance
CN104991202B (en) A kind of more battery combination system estimation SOC method and relevant apparatus
CN107947197B (en) Low-frequency oscillation mode and mode identification method for power system
CN104062501A (en) Double-transformer substation harmonic wave state estimation method
CN106920047A (en) Power network statistics based on level storage model derives from computational methods and system

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