CN115800271B - Power distribution system parameter correction method and system based on self-adaptive Kalman filtering - Google Patents

Power distribution system parameter correction method and system based on self-adaptive Kalman filtering Download PDF

Info

Publication number
CN115800271B
CN115800271B CN202310048743.5A CN202310048743A CN115800271B CN 115800271 B CN115800271 B CN 115800271B CN 202310048743 A CN202310048743 A CN 202310048743A CN 115800271 B CN115800271 B CN 115800271B
Authority
CN
China
Prior art keywords
covariance
distribution system
value
calculating
adaptive
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
CN202310048743.5A
Other languages
Chinese (zh)
Other versions
CN115800271A (en
Inventor
祖国强
黄旭
丁琪
李治
晋萃萃
张春晖
魏然
贺春
冯郁竹
徐智
吴俣
张军
邱巧红
魏炜
季文文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd
Chengdong Power Supply Co of State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd
Chengdong Power Supply Co of State Grid Tianjin Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Tianjin Electric Power Co Ltd, Electric Power Research Institute of State Grid Tianjin Electric Power Co Ltd, Chengdong Power Supply Co of State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202310048743.5A priority Critical patent/CN115800271B/en
Publication of CN115800271A publication Critical patent/CN115800271A/en
Application granted granted Critical
Publication of CN115800271B publication Critical patent/CN115800271B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a power distribution system parameter correction method and system based on self-adaptive Kalman filtering, wherein the method comprises the following steps: establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering; calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation; calculating residual errors and covariance thereof according to the predicted values and covariance thereof; calculating a mahalanobis distance according to the residual error and the covariance; carrying out self-adaptive correction on the residual error according to the mahalanobis distance; then, calculating Kalman filtering gain through residual error and covariance thereof; and calculating a posterior state estimation value according to the Kalman filtering gain. The invention improves the observability and the computing capability of the medium-voltage distribution network.

Description

Power distribution system parameter correction method and system based on self-adaptive Kalman filtering
Technical Field
The invention belongs to the field of situation awareness of medium-voltage power distribution systems, and particularly relates to a power distribution system parameter correction method and system based on self-adaptive Kalman filtering.
Background
Along with the development of society and the improvement of the living standard of people, the requirements on the power supply reliability and the power supply quality of a power system are higher and higher. The medium voltage distribution system as an important component of the power system is directly oriented to end users, and the perfection of the medium voltage distribution system is directly related to the electricity reliability and the electricity quality of the vast users. In the conventional calculation of the power flow for a medium-voltage power distribution system, the parameters of the medium-voltage power distribution system are generally approximated unchanged. In fact, the parameters of the medium voltage distribution system change at any time according to seasons and time, and the observability and the computing capacity of the medium voltage distribution network are affected by the traditional approximation process. Therefore, how to improve the observability of the medium-voltage distribution network and ensure the operation safety and the power supply reliability of the power network becomes a problem to be solved urgently.
Disclosure of Invention
In view of the above problems, the present invention provides a power distribution system parameter correction method based on adaptive kalman filtering, including:
establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering;
calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation;
calculating residual errors and covariance thereof according to the predicted values and the covariance and the measured values;
calculating a mahalanobis distance according to the residual error and the covariance thereof, and carrying out self-adaptive correction on the residual error according to the mahalanobis distance, wherein the method comprises the following steps:
determining whether the corresponding measurement value is uncertain noise according to the mahalanobis distance, and if so, updating the residual error and covariance thereof through the self-adaptive factor;
determining whether the corresponding measured value is a wild value according to the Marsh distance, if so, discarding the wild value, and recalculating residual errors and covariance thereof by replacing the measured value with a predicted value;
then, calculating Kalman filtering gain through the latest residual error and covariance thereof;
and calculating a posterior state estimation value according to the Kalman filtering gain.
Further, establishing a state equation and a measurement equation of the system parameters of the medium-voltage distribution system based on the adaptive Kalman filtering comprises:
establishing a state equation and a measurement equation of added noise:
Figure SMS_1
(10)
Figure SMS_6
state variables representing system parameters, +.>
Figure SMS_10
Representing the previous timeI.e. +.>
Figure SMS_14
Status of system parameters at time, +.>
Figure SMS_3
Is a state transition matrix, ">
Figure SMS_9
,/>
Figure SMS_13
Zero mean and covariance +.>
Figure SMS_19
Is white gaussian noise; />
Figure SMS_5
Representing the measured value->
Figure SMS_15
Is a measurement matrix->
Figure SMS_21
,/>
Figure SMS_25
Represents the voltage amplitude of the line current inflow at time k,/->
Figure SMS_7
Represents the current through the line at time k, < >>
Figure SMS_12
Representing the current at the line current inflow at time k-1, < >>
Figure SMS_18
Representing the line current inflow terminal voltage amplitude at time k-1,
Figure SMS_23
representing the voltage amplitude of the line current inflow terminal at the time k-2; />
Figure SMS_20
Zero mean and covariance +.>
Figure SMS_24
Is Gaussian white noise, < >>
Figure SMS_26
Is extra noise->
Figure SMS_27
To measure the wild value, wherein->
Figure SMS_2
Obeying the parameter +.>
Figure SMS_11
Is characterized by a Bernoulli distribution,
Figure SMS_17
,/>
Figure SMS_22
is a unit pulse function>
Figure SMS_4
Obeying the parameter +.>
Figure SMS_8
Bernoulli distribution, ->
Figure SMS_16
Wherein the system parameters include line parameters and/or transformer parameters.
Further, calculating the predicted value and covariance of the line parameter at a certain moment includes: obtaining initial values and covariance of line parameters of medium-voltage distribution system
Figure SMS_28
and />
Figure SMS_29
Calculated by the following formula:
Figure SMS_30
(11)
Figure SMS_31
(12)
knowing the initial value
Figure SMS_32
And covariance + ->
Figure SMS_33
The value and the covariance of the state variable at the next moment are obtained, then the value and the covariance of the state variable at the later moment are obtained through the value and the covariance of the state variable at the next moment, and iteration is continuously carried out to obtain +.>
Figure SMS_34
Predicted value of system parameter of time medium voltage distribution system>
Figure SMS_35
And covariance + ->
Figure SMS_36
Further, calculate the residual
Figure SMS_37
And covariance + ->
Figure SMS_38
The formula of (2) is as follows:
Figure SMS_39
(13)
Figure SMS_40
(14)。
further, calculating the mahalanobis distance from the residual and its covariance includes:
Figure SMS_41
(15)
wherein (r) represents the r-th element, (r, r) represents the r-th row, and the r-th column of elements.
Further, if the mahalanobis distance is within the specified confidence interval, indicating that uncertain noise appears in measurement, updating the residual error and covariance thereof by the self-adaptive factor;
adaptive factor
Figure SMS_42
And its updated residual covariance->
Figure SMS_43
The calculation formula of (2) is as follows: />
Figure SMS_44
(16)
Figure SMS_45
(17)
wherein ,
Figure SMS_47
indicating that the lower boundary of the assigned confidence interval, < ->
Figure SMS_51
Is the residual covariance before the update,
Figure SMS_55
for the updated residual covariance +.>
Figure SMS_48
For vector->
Figure SMS_52
Is>
Figure SMS_56
Element(s)>
Figure SMS_58
Is vector quantity
Figure SMS_49
Is>
Figure SMS_53
Element(s)>
Figure SMS_57
、/>
Figure SMS_59
、/>
Figure SMS_46
For matrix->
Figure SMS_50
Diagonal->
Figure SMS_54
The elements.
Further, calculating the kalman filter gain from the residual and its covariance includes:
kalman filtering gain of medium voltage distribution system parameters without occurrence of uncertain noise
Figure SMS_60
The calculation formula of (2) is as follows:
Figure SMS_61
(18)
if uncertain noise occurs, then:
Figure SMS_62
(19)。
further, posterior state estimation
Figure SMS_63
And covariance + ->
Figure SMS_64
Computing means of (a)The formula is as follows:
Figure SMS_65
(20)
Figure SMS_66
(21) 。
further, if the mahalanobis distance is greater than the upper limit of the specified confidence interval, the corresponding measured value is considered to be a wild value and discarded, and the predicted value is adopted instead.
The invention also provides a distribution system parameter correction system based on the self-adaptive Kalman filtering, which comprises:
a model determination unit configured to: establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering;
a prediction unit for: calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation;
a residual calculation unit for: calculating residual errors and covariance thereof according to the predicted values and covariance thereof;
the residual error correction unit is used for calculating a mahalanobis distance according to the residual error and covariance thereof, and carrying out self-adaptive correction on the residual error according to the mahalanobis distance, and comprises the following steps:
determining whether the corresponding measurement value is uncertain noise according to the mahalanobis distance, and if so, updating the residual error and covariance thereof through the self-adaptive factor;
determining whether the corresponding measured value is a wild value according to the Marsh distance, if so, discarding the wild value, and recalculating residual errors and covariance thereof by replacing the measured value with a predicted value;
gain calculation unit for: calculating a Kalman filtering gain by using the residual error and the covariance determined by the residual error correction unit;
an estimated value calculation unit configured to: and calculating a posterior state estimation value according to the Kalman filtering gain.
According to the power distribution system parameter correction method and system based on the self-adaptive Kalman filtering, the self-adaptive Kalman filtering can utilize a large amount of measured value data to estimate and correct the model and noise statistics characteristics, and a large amount of historical data is used for estimating the parameters of the medium-voltage power distribution system to obtain a changed new parameter value, so that the observability and the computing capacity of the medium-voltage power distribution network are improved.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 shows a schematic diagram of an equivalent circuit structure of a medium voltage distribution system line according to an embodiment of the present invention;
fig. 2 is a schematic diagram showing an equivalent circuit structure of a transformer of a medium-voltage distribution system according to an embodiment of the present invention;
FIG. 3 illustrates a flow chart of a method for on-line correction of parameters of a medium voltage power distribution system in accordance with an embodiment of the present invention;
fig. 4 shows a schematic structural diagram of an on-line parameter correction system for a medium voltage distribution system according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention provides a power distribution system parameter correction method based on self-adaptive Kalman filtering, which can be used for on-line correction of medium-voltage power distribution system parameters. The method corrects parameters by utilizing a large amount of historical measurement data recorded in the power grid, and judges and processes the outlier and the uncertainty noise existing in the measurement information by utilizing a hypothesis testing method. Compared with the existing calculation by using unchanged system parameters, the method can consider the condition that the parameters of the medium-voltage distribution system change at any time along with the difference of seasons and time, and improves the observability and the calculation capability of the medium-voltage distribution network.
In the embodiment of the invention, the equivalent circuit structure of the medium voltage distribution system circuit is shown in figure 1. The equivalent resistance R and the equivalent reactance X of the circuit are connected in series, and the voltages at the two ends of the circuit are respectively
Figure SMS_67
The current flowing into the line is +.>
Figure SMS_68
The current in the line is I. The equivalent capacitance to ground at both ends of the line (ground points GND1 and GND2 respectively) is C.
The equivalent circuit structure of the medium voltage distribution system transformer is shown in fig. 2. The equivalent circuit structure comprises an excitation part and a gate-type equivalent circuit which are connected in series. The transformation ratio of the transformer is T. The exciting section is represented by a parallel branch formed by connecting a conductance G and a capacitance C in parallel, one end of the parallel branch is grounded, and the other end is used as a current inflow end of the circuit structure. Other performance uses of transformers
Figure SMS_81
、/>
Figure SMS_70
、/>
Figure SMS_76
Constitution ofGate (n-type) equivalent circuit representation, wherein->
Figure SMS_73
,/>
Figure SMS_75
Figure SMS_74
, wherein ,/>
Figure SMS_80
One end is grounded, the other end is connected with->
Figure SMS_79
Is connected with the current inflow end of (a)>
Figure SMS_83
One end is grounded, the other end is connected with->
Figure SMS_69
Is connected to the current outlet terminal of the (c). Wherein Z represents impedance and consists of an equivalent resistor R and an equivalent reactor X which are connected in series. The voltages at the two ends of the line are respectively->
Figure SMS_77
The current flowing into the line at the head end is +.>
Figure SMS_72
Flow through->
Figure SMS_78
The current of (2) is I, flow->
Figure SMS_82
Is +.>
Figure SMS_84
The current of the excitation part is +.>
Figure SMS_71
An on-line correction flow chart of the medium voltage distribution system parameters based on the adaptive Kalman filtering is shown in figure 3.
The method comprises the following steps:
step 1: establishing system parameters of medium voltage distribution system
Figure SMS_85
State equation and measurement equation based on adaptive Kalman filtering, and gives (can be directly obtained according to the system parameter performance) the process noise covariance +.>
Figure SMS_86
Measurement noise covariance->
Figure SMS_87
In the embodiment of the invention, the system parameters comprise line parameters and/or transformer parameters, and the processing logic of the two parameters is the same by the power distribution system parameter correction method based on the self-adaptive Kalman filtering, so that the parameters can be corrected in an omnibearing and multi-angle manner, and the accuracy of power grid observation and calculation is improved.
The establishment of state equations and measurement equations of system parameters of the medium-voltage distribution system is divided into the following two parts:
(1) Establishing a state equation and a measurement equation without noise interference
The state equation is established as follows:
Figure SMS_88
(1)
wherein ,
Figure SMS_89
discrete time series>
Figure SMS_90
Representing the state of the system parameter at time k +.>
Figure SMS_91
Representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>
Figure SMS_92
Is a state transition matrix, ">
Figure SMS_93
. For line parameters, state variables ∈ ->
Figure SMS_94
The measurement equation is established below, and equivalent reactance is set for the line parameters
Figure SMS_95
The voltage at two ends is +>
Figure SMS_96
In the continuous time domain, according to kirchhoff's law:
easily-known
Figure SMS_97
(2)
Figure SMS_98
(3)
Figure SMS_99
(4)/>
Figure SMS_100
Representing the voltage amplitude of the line current outflow terminal at the moment t, < >>
Figure SMS_101
Representing the line current inflow terminal voltage amplitude at time t, < >>
Figure SMS_102
Representing the equivalent reactance +.>
Figure SMS_103
Voltage at two ends>
Figure SMS_104
Representing the flow at time tCurrent through the line->
Figure SMS_105
The current at the line current inflow end at time t is shown. Since the sampling points are discrete, it is necessary to cross-domain the continuous time domain to the discrete time domain:
Figure SMS_106
(5)
Figure SMS_107
(6)
Figure SMS_108
(7)
wherein ,
Figure SMS_111
are all continuous time functions, and are sampled (in discrete time domain) as
Figure SMS_114
Abbreviated as +.>
Figure SMS_116
,/>
Figure SMS_110
Is the sampling time. I.e.)>
Figure SMS_113
Represents the current at the line current inflow at time k, < >>
Figure SMS_115
Representing the current at the line current inflow at time k-1,
Figure SMS_117
represents the line current inflow terminal voltage amplitude at time k-1,/->
Figure SMS_109
Representing the line at time kCurrent inflow terminal voltage amplitude, ">
Figure SMS_112
Representing the voltage amplitude of the line current outflow end at the moment k;
substituting (4), (5) and (6) into (1) to obtain:
Figure SMS_118
Figure SMS_119
Figure SMS_120
(8)
the preparation method comprises the following steps of:
Figure SMS_121
(9)
Figure SMS_122
/>
t is transformer transformation ratio, impedance
Figure SMS_123
J represents the imaginary part of the complex number, < ->
Figure SMS_124
Figure SMS_125
Figure SMS_126
(2) And (3) establishing a measurement equation and a state equation with noise influence as a system parameter correction model of the adaptive Kalman filtering.
Adding noise interference to a state equation and a measurement equation respectively to obtain:
Figure SMS_127
(10)
Figure SMS_145
zero mean and covariance +.>
Figure SMS_130
Is a gaussian white noise of (c). Measurement value->
Figure SMS_139
,
Figure SMS_143
The voltage amplitude at the end of the line (line of the system), respectively +.>
Figure SMS_149
Is a measurement matrix of the measurement data,
Figure SMS_144
,/>
Figure SMS_148
representing the line current inflow terminal voltage amplitude at time k-2. />
Figure SMS_132
Zero mean and covariance +.>
Figure SMS_136
Is Gaussian white noise, < >>
Figure SMS_128
Is extra noise->
Figure SMS_134
To measure the wild value, wherein->
Figure SMS_131
Obeying the parameter +.>
Figure SMS_135
Is characterized by a Bernoulli distribution,/>
Figure SMS_141
,/>
Figure SMS_147
for a constant of larger amplitude, illustratively, < +.>
Figure SMS_133
Between 100 and 1000%>
Figure SMS_137
Is a unit pulse function>
Figure SMS_140
Obeying the parameter +.>
Figure SMS_146
Is characterized by a Bernoulli distribution,
Figure SMS_129
the initial value and covariance of the system parameters of the medium voltage distribution system are +.>
Figure SMS_138
and />
Figure SMS_142
Step 2: system parameter correction model (namely formula (10)) based on adaptive Kalman filtering and calculation
Figure SMS_150
Predicted value of system parameter of time medium voltage distribution system>
Figure SMS_151
And covariance + ->
Figure SMS_152
And the initial value of the state variable +.>
Figure SMS_153
And covariance +.>
Figure SMS_154
The predicted value and covariance calculation formula are as follows:
Figure SMS_155
(11)
Figure SMS_156
(12)/>
knowing the initial value
Figure SMS_157
And covariance + ->
Figure SMS_158
The value and the covariance of the state variable at the next moment are obtained, then the value and the covariance of the state variable at the later moment are obtained through the value and the covariance of the state variable at the next moment, and iteration is continuously carried out to obtain +.>
Figure SMS_159
Predicted value of system parameter of time medium voltage distribution system>
Figure SMS_160
And covariance + ->
Figure SMS_161
. Initial value->
Figure SMS_162
And covariance + ->
Figure SMS_163
Is obtained by consulting a standard for the known cable used, i.e. the value obtained by subjecting the cable to a factory test.
Step 3: measuring matrix obtained by using current value and load power obtained by measuring device, distribution network three-phase tide algorithm and state estimation
Figure SMS_164
And calculates the residual +.>
Figure SMS_165
Covariance->
Figure SMS_166
Calculating residual errors
Figure SMS_167
And covariance + ->
Figure SMS_168
The formula of (2) is as follows:
Figure SMS_169
(13)
Figure SMS_170
(14)
step 4: calculation of the mahalanobis distance from the resulting residual
Figure SMS_171
If->
Figure SMS_172
Then go to step 5; otherwise, go to step 6;
utilizing residual errors
Figure SMS_173
And covariance + ->
Figure SMS_174
The calculated mahalanobis distance is calculated as follows:
Figure SMS_175
(15)
wherein (r) represents an r-th element, (r, r) represents an r-th row, and an r-th column of elements;
step 5: comparison of
Figure SMS_176
and />
Figure SMS_177
The size of (1)>
Figure SMS_178
The corresponding measurement value is regarded as uncertain noise, and the adaptive factor is calculated>
Figure SMS_179
And updating the prior residual, calculating the updated residual covariance +.>
Figure SMS_180
Then, go to step 7; otherwise, directly turning to the step 7;
judging and processing modes of the outlier and the uncertainty noise are as follows: calculation of the Marshall distance of the residual
Figure SMS_181
If (if)
Figure SMS_182
Then the measurement outlier is considered to appear, is removed and is replaced by a predicted value (obtained by the formula (11)) to replace the posterior state estimation value (used in the formula (20)); when the outlier appears, the covariance is recalculated by using the predicted value instead of the outlier, the covariance is 0, and the kalman filter gain is 0 according to the subsequent calculation, which corresponds to replacing the posterior state estimation value by the predicted value. If->
Figure SMS_183
Then it is considered that uncertainty noise is present at this time. Adaptive factor->
Figure SMS_184
And its updated residual covariance->
Figure SMS_185
The calculation formula of (2) is as follows: />
Figure SMS_186
(16)
Figure SMS_187
(17)
wherein ,
Figure SMS_189
、/>
Figure SMS_195
the lower and upper boundaries of the confidence interval are represented, illustratively by values of 0 and 1, respectively.
Figure SMS_199
Is the residual covariance before update, +.>
Figure SMS_190
For the updated residual covariance +.>
Figure SMS_194
For vector->
Figure SMS_198
Is>
Figure SMS_202
Element(s)>
Figure SMS_188
For vector->
Figure SMS_192
Is>
Figure SMS_196
Element(s)>
Figure SMS_200
、/>
Figure SMS_191
Figure SMS_193
For matrix->
Figure SMS_197
Diagonal->
Figure SMS_201
The elements.
The invention judges the wild value and the uncertain noise through the mahalanobis distance and carries out corresponding processing, thereby realizing the self-adaptive correction of the parameters.
Step 6: taking the measured value as an outlier, discarding, replacing the posterior state estimation value with the predicted value obtained in the step 2, and turning to the step 7;
step 7: kalman filtering gain for calculating parameters of medium voltage distribution system
Figure SMS_203
Kalman filtering gain for medium voltage distribution system parameters
Figure SMS_204
The calculation formula of (2) is as follows:
Figure SMS_205
(18)
if uncertain noise occurs, then:
Figure SMS_206
(19)。
step 8: by means of
Figure SMS_207
Calculated->
Figure SMS_208
Posterior state estimation value of time-of-day system parameters +.>
Figure SMS_209
And covariance + ->
Figure SMS_210
(matrix) and will->
Figure SMS_211
As a result of the correction of system parameters.
Medium voltage power distribution system
Figure SMS_212
Posterior state estimation value of time-of-day system parameters +.>
Figure SMS_213
And covariance + ->
Figure SMS_214
The calculation formula of (2) is as follows:
Figure SMS_215
(20)
Figure SMS_216
(21)
based on the same inventive concept, the embodiment of the invention also provides a power distribution system parameter correction system based on adaptive Kalman filtering, as shown in fig. 4, comprising:
a model determination unit configured to: establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering;
a prediction unit for: calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation;
a residual calculation unit for: calculating residual errors and covariance thereof according to the predicted values and covariance thereof;
the residual error correction unit is used for calculating a mahalanobis distance according to the residual error and covariance thereof, and carrying out self-adaptive correction on the residual error according to the mahalanobis distance, and comprises the following steps of;
determining whether the corresponding measurement value is uncertain noise according to the mahalanobis distance, and if so, updating the residual error and covariance thereof through the self-adaptive factor;
determining whether the corresponding measured value is a wild value according to the Marsh distance, if so, discarding the wild value, and recalculating residual errors and covariance thereof by replacing the measured value with a predicted value;
gain calculation unit for: calculating a Kalman filtering gain by using the residual error and the covariance determined by the residual error correction unit;
an estimated value calculation unit configured to: and calculating a posterior state estimation value according to the Kalman filtering gain.
The residual error correction unit is used for judging whether the mahalanobis distance is larger than the upper limit value of the appointed confidence interval, if so, the corresponding measured value is considered to be a wild value and is discarded, and the predicted value is adopted to replace the posterior state estimated value in the estimated value calculation unit; and judging whether the Marsh distance is within a specified confidence interval, if so, updating the residual error and covariance thereof through the self-adaptive factor, wherein the uncertainty noise appears in measurement.
Each unit specific implementation of the distribution system parameter correction system based on the adaptive Kalman filter can be obtained according to the distribution system parameter correction method based on the adaptive Kalman filter according to the embodiment of the invention, and the description is omitted.
Although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (10)

1. The power distribution system parameter correction method based on the adaptive Kalman filtering is characterized by comprising the following steps of:
establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering;
calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation;
calculating residual errors and covariance thereof according to the predicted values and the covariance and the measured values;
calculating a mahalanobis distance according to the residual error and the covariance thereof, and carrying out self-adaptive correction on the residual error according to the mahalanobis distance, wherein the method comprises the following steps:
determining whether the corresponding measurement value is uncertain noise according to the mahalanobis distance, and if so, updating the residual error and covariance thereof through the self-adaptive factor;
determining whether the corresponding measured value is a wild value according to the Marsh distance, if so, discarding the wild value, and recalculating residual errors and covariance thereof by replacing the measured value with a predicted value;
then, calculating Kalman filtering gain through the latest residual error and covariance thereof;
calculating a posterior state estimation value according to the Kalman filtering gain;
wherein the system parameters include line parameters;
in the equivalent circuit structure of the medium-voltage distribution system circuit, an equivalent resistor R and an equivalent reactance X are connected in series, and voltages at two ends of the circuit are respectively
Figure QLYQS_1
The current flowing into the line is +.>
Figure QLYQS_2
The current in the circuit is I, and the equivalent capacitance to the ground at the two ends of the circuit is C; the state equation and the measurement equation based on the adaptive Kalman filtering for establishing the system parameters of the medium-voltage distribution system comprise:
establishing a state equation and a measurement equation of added noise:
Figure QLYQS_3
(10)
Figure QLYQS_7
a state variable representing a parameter of the system,
Figure QLYQS_12
,/>
Figure QLYQS_19
representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>
Figure QLYQS_5
Is a state transition matrix, ">
Figure QLYQS_13
,/>
Figure QLYQS_20
Zero mean and covariance +.>
Figure QLYQS_26
Is white gaussian noise; />
Figure QLYQS_9
Representing the measured value->
Figure QLYQS_16
Is a measurement matrix of the measurement data,
Figure QLYQS_24
,/>
Figure QLYQS_29
represents the voltage amplitude of the line current inflow at time k,/->
Figure QLYQS_10
Represents the current at the line current inflow at time k, +.>
Figure QLYQS_17
Representing the current at the line current inflow at time k-1, < >>
Figure QLYQS_23
Representing the line current inflow terminal voltage amplitude at time k-1,
Figure QLYQS_28
representing the voltage amplitude of the line current inflow terminal at the time k-2; />
Figure QLYQS_6
Zero mean and covariance +.>
Figure QLYQS_14
Is Gaussian white noise, < >>
Figure QLYQS_21
Is extra noise->
Figure QLYQS_27
To measure the wild value, wherein->
Figure QLYQS_4
Obeying the parameter +.>
Figure QLYQS_11
Is characterized by a Bernoulli distribution,
Figure QLYQS_18
,/>
Figure QLYQS_25
is a unit pulse function>
Figure QLYQS_8
Obeying the parameter +.>
Figure QLYQS_15
Bernoulli distribution, ->
Figure QLYQS_22
2. The adaptive Kalman filtering-based power distribution system parameter modification method of claim 1, wherein the system parameters further comprise transformer parameters,
for transformer parameters, the equivalent circuit structure of the transformer of the medium-voltage distribution system comprises an excitation part and a gate-type equivalent circuit, the transformation ratio of the transformer is T, the excitation part is represented by a parallel branch formed by connecting a conductance G and a capacitance C in parallel, one end of the parallel branch is grounded, the other end of the parallel branch is used as a current inflow end of the circuit structure, and other performances of the transformer are used
Figure QLYQS_38
、/>
Figure QLYQS_32
、/>
Figure QLYQS_46
A gate-type equivalent circuit representation is constructed, wherein +.>
Figure QLYQS_31
,/>
Figure QLYQS_42
,/>
Figure QLYQS_34
, wherein ,/>
Figure QLYQS_40
One end is grounded, the other end is connected with->
Figure QLYQS_33
Is connected with the current inflow end of (a)>
Figure QLYQS_44
One end is grounded, the other end is connected with->
Figure QLYQS_30
Wherein Z represents impedance and consists of an equivalent resistor R and an equivalent reactor X which are connected in series, and voltages at two ends of the line are respectively +.>
Figure QLYQS_47
The current flowing into the line at the head end is +.>
Figure QLYQS_35
Flow through->
Figure QLYQS_43
The current of (2) is I, flow->
Figure QLYQS_37
Is +.>
Figure QLYQS_41
The current of the excitation part is +.>
Figure QLYQS_39
,/>
Figure QLYQS_45
Is connected with the current inflow end of the equivalent circuit structure of the transformer, +.>
Figure QLYQS_36
The current outflow end of the transformer equivalent circuit structure;
in the state equation and the measurement equation,
Figure QLYQS_48
3. the adaptive kalman filter based power distribution system parameter correction method according to claim 2, wherein calculating the predicted value of the system parameter and the covariance thereof at a certain moment comprises: obtaining initial values and covariance of system parameters of a medium voltage distribution system
Figure QLYQS_49
and />
Figure QLYQS_50
Calculated by the following formula:
Figure QLYQS_51
(11)
Figure QLYQS_52
(12)
knowing the initial value
Figure QLYQS_53
And covariance + ->
Figure QLYQS_54
Continuously iterating to obtain ++>
Figure QLYQS_55
Predicted value of system parameter of time medium voltage distribution system>
Figure QLYQS_56
And covariance + ->
Figure QLYQS_57
,/>
Figure QLYQS_58
and />
Figure QLYQS_59
The posterior state estimates and covariance of the k-1 time are obtained.
4. The adaptive kalman filter based power distribution system parameter correction method according to claim 3, wherein residual error is calculated
Figure QLYQS_60
And covariance + ->
Figure QLYQS_61
The formula of (2) is as follows:
Figure QLYQS_62
(13)
Figure QLYQS_63
(14)。
5. the adaptive kalman filter based power distribution system parameter correction method according to claim 4, wherein calculating the mahalanobis distance based on the residual error and its covariance comprises:
Figure QLYQS_64
(15)
where (r) represents the r-th element, (r, r) represents the r-th row, and the r-th column.
6. The adaptive Kalman filtering-based power distribution system parameter modification method of claim 5, wherein,
if the mahalanobis distance is within the specified confidence interval, indicating that uncertain noise appears in measurement, updating the residual error and covariance thereof through the self-adaptive factor;
adaptive factor
Figure QLYQS_65
And its updated residual covariance->
Figure QLYQS_66
The calculation formula of (2) is as follows:
Figure QLYQS_67
(16)
Figure QLYQS_68
(17)
wherein ,
Figure QLYQS_72
representing the lower boundary of the specified confidence interval, +.>
Figure QLYQS_75
Is the residual covariance before update, +.>
Figure QLYQS_79
For the updated residual covariance +.>
Figure QLYQS_82
For vector->
Figure QLYQS_69
Is>
Figure QLYQS_73
Element(s)>
Figure QLYQS_77
Figure QLYQS_81
、/>
Figure QLYQS_70
Respectively is matrix->
Figure QLYQS_76
and />
Figure QLYQS_80
Diagonal->
Figure QLYQS_83
The elements.
7. The adaptive kalman filter based power distribution system parameter correction method according to claim 6, wherein calculating the kalman filter gain by the latest residual and covariance thereof includes:
kalman filtering gain of medium voltage distribution system parameters without occurrence of uncertain noise
Figure QLYQS_84
The calculation formula of (2) is as follows:
Figure QLYQS_85
(18)
if uncertain noise occurs, then:
Figure QLYQS_86
(19)。
8. the adaptive kalman filter based power distribution system parameter correction method according to claim 7, wherein the posterior state estimation value
Figure QLYQS_87
And covariance + ->
Figure QLYQS_88
The calculation formula of (2) is as follows:
Figure QLYQS_89
(20)
Figure QLYQS_90
(21)。
9. an adaptive kalman filter based power distribution system parameter correction method as in any of the claims 6-8, wherein,
if the mahalanobis distance is greater than the upper limit value of the specified confidence interval, the corresponding measurement value is considered as a wild value and discarded, and the predicted value is adopted for replacement.
10. An adaptive kalman filter-based power distribution system parameter correction system, comprising:
a model determination unit configured to: establishing a state equation and a measurement equation of system parameters of the medium-voltage distribution system based on the self-adaptive Kalman filtering;
a prediction unit for: calculating a predicted value and covariance of a system parameter at a certain moment based on the state equation and the measurement equation;
a residual calculation unit for: calculating residual errors and covariance thereof according to the predicted values and covariance thereof;
the residual error correction unit is used for calculating a mahalanobis distance according to the residual error and covariance thereof, and carrying out self-adaptive correction on the residual error according to the mahalanobis distance, and comprises the following steps:
determining whether the corresponding measurement value is uncertain noise according to the mahalanobis distance, and if so, updating the residual error and covariance thereof through the self-adaptive factor;
determining whether the corresponding measured value is a wild value according to the Marsh distance, if so, discarding the wild value, and recalculating residual errors and covariance thereof by replacing the measured value with a predicted value;
gain calculation unit for: calculating a Kalman filtering gain by using the residual error and the covariance determined by the residual error correction unit;
an estimated value calculation unit configured to: calculating a posterior state estimation value according to the Kalman filtering gain;
wherein the system parameters include line parameters;
in the equivalent circuit structure of the medium-voltage distribution system circuit, an equivalent resistor R and an equivalent reactance X are connected in series, and voltages at two ends of the circuit are respectively
Figure QLYQS_91
The current flowing into the line is +.>
Figure QLYQS_92
The current in the line is IThe equivalent capacitance to ground at the two ends of the circuit is C; the state equation and the measurement equation based on the adaptive Kalman filtering for establishing the system parameters of the medium-voltage distribution system comprise:
establishing a state equation and a measurement equation of added noise:
Figure QLYQS_93
(10)
Figure QLYQS_95
a state variable representing a parameter of the system,
Figure QLYQS_103
,/>
Figure QLYQS_110
representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>
Figure QLYQS_100
Is a state transition matrix, ">
Figure QLYQS_106
,/>
Figure QLYQS_113
Zero mean and covariance +.>
Figure QLYQS_118
Is white gaussian noise; />
Figure QLYQS_97
Representing the measured value->
Figure QLYQS_105
Is a measurement matrix of the measurement data,
Figure QLYQS_112
,/>
Figure QLYQS_117
represents the voltage amplitude of the line current inflow at time k,/->
Figure QLYQS_99
Represents the current at the line current inflow at time k, +.>
Figure QLYQS_107
Representing the current at the line current inflow at time k-1, < >>
Figure QLYQS_114
Representing the line current inflow terminal voltage amplitude at time k-1,
Figure QLYQS_119
representing the voltage amplitude of the line current inflow terminal at the time k-2; />
Figure QLYQS_96
Zero mean and covariance +.>
Figure QLYQS_102
Is Gaussian white noise, < >>
Figure QLYQS_109
Is extra noise->
Figure QLYQS_115
To measure the wild value, wherein->
Figure QLYQS_94
Obeying the parameter +.>
Figure QLYQS_101
Is characterized by a Bernoulli distribution,
Figure QLYQS_108
,/>
Figure QLYQS_116
is a unit pulse function>
Figure QLYQS_98
Obeying the parameter +.>
Figure QLYQS_104
Bernoulli distribution, ->
Figure QLYQS_111
。/>
CN202310048743.5A 2023-02-01 2023-02-01 Power distribution system parameter correction method and system based on self-adaptive Kalman filtering Active CN115800271B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310048743.5A CN115800271B (en) 2023-02-01 2023-02-01 Power distribution system parameter correction method and system based on self-adaptive Kalman filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310048743.5A CN115800271B (en) 2023-02-01 2023-02-01 Power distribution system parameter correction method and system based on self-adaptive Kalman filtering

Publications (2)

Publication Number Publication Date
CN115800271A CN115800271A (en) 2023-03-14
CN115800271B true CN115800271B (en) 2023-05-23

Family

ID=85429410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310048743.5A Active CN115800271B (en) 2023-02-01 2023-02-01 Power distribution system parameter correction method and system based on self-adaptive Kalman filtering

Country Status (1)

Country Link
CN (1) CN115800271B (en)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110530365B (en) * 2019-08-05 2021-05-18 浙江工业大学 Human body attitude estimation method based on adaptive Kalman filtering

Also Published As

Publication number Publication date
CN115800271A (en) 2023-03-14

Similar Documents

Publication Publication Date Title
CN110187282B (en) Battery state of charge estimation method and estimation device
CN111366855A (en) Battery equivalent circuit model disturbance-resistant parameterization method
CN107525526B (en) Device and method for providing a measurement of the uncertainty of a transfer function
EP3132514A1 (en) Transformer parameter estimation using terminal measurements
CN111046327A (en) Prony analysis method suitable for low-frequency oscillation and subsynchronous oscillation identification
Abbas A new fast algorithm to estimate real-time phasors using adaptive signal processing
CN110749835A (en) Power transmission line fault positioning method based on Kalman filter
CN108535597B (en) Line model-based single-phase earth fault section positioning method
CN113625182A (en) Method for on-line estimating battery state
CN115800271B (en) Power distribution system parameter correction method and system based on self-adaptive Kalman filtering
Roj et al. Method of measurement of capacitance and dielectric loss factor using artificial neural networks
CN117310455B (en) Oral cavity scanner circuit board fault detection method based on electric signals
CN112083299B (en) DC system insulation fault prediction method based on Kalman filtering
CN109655775A (en) A kind of multiple dimensioned calibration method of arbitrary waveform generator amplitude frequency sweep and device
CN112213560A (en) High-precision power grid broadband signal measuring method based on Z-ADALINE
CN112009252A (en) Fault diagnosis and fault-tolerant control method for power battery system
CN116794517A (en) Lithium ion battery SOC estimation method and system based on fractional order Kalman filtering
Serra et al. Fast ADC testing by repetitive histogram analysis
CN116340810A (en) Transformer outgoing line part partial discharge diagnosis method and device based on random forest
CN112601968A (en) Charge state online estimation system for 48V light-mixed automobile lithium ion battery
CN113985210B (en) Voltage sag domain calculation method considering voltage sag amplitude and duration
Grimaldi et al. Identification of ADC error model by testing of the chosen code bins
CN112700039B (en) Steady state detection and extraction method for load operation data of thermal power plant
CN115130550A (en) Distribution transformer fault identification method based on gradient lifting decision tree
CN109327026B (en) Low-voltage distribution network interval state estimation method

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