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 PDFInfo
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 43
- 238000012937 correction Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims description 27
- 230000003044 adaptive effect Effects 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000002715 modification method Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 101100489713 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) GND1 gene Proteins 0.000 description 1
- 101100489717 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) GND2 gene Proteins 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
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:
state variables representing system parameters, +.>Representing the previous timeI.e. +.>Status of system parameters at time, +.>Is a state transition matrix, ">,/>Zero mean and covariance +.>Is white gaussian noise; />Representing the measured value->Is a measurement matrix->,/>Represents the voltage amplitude of the line current inflow at time k,/->Represents the current through the line at time k, < >>Representing the current at the line current inflow at time k-1, < >>Representing the line current inflow terminal voltage amplitude at time k-1,representing the voltage amplitude of the line current inflow terminal at the time k-2; />Zero mean and covariance +.>Is Gaussian white noise, < >>Is extra noise->To measure the wild value, wherein->Obeying the parameter +.>Is characterized by a Bernoulli distribution,,/>is a unit pulse function>Obeying the parameter +.>Bernoulli distribution, ->;
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 and />,
Calculated by the following formula:
knowing the initial valueAnd covariance + ->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 +.>Predicted value of system parameter of time medium voltage distribution system>And covariance + ->。
further, calculating the mahalanobis distance from the residual and its covariance includes:
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 factorAnd its updated residual covariance->The calculation formula of (2) is as follows: />
wherein ,indicating that the lower boundary of the assigned confidence interval, < ->Is the residual covariance before the update,for the updated residual covariance +.>For vector->Is>Element(s)>Is vector quantityIs>Element(s)>、/>、/>For matrix->Diagonal->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 noiseThe calculation formula of (2) is as follows:
if uncertain noise occurs, then:
further, posterior state estimationAnd covariance + ->Computing means of (a)The formula is as follows:
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 respectivelyThe current flowing into the line is +.>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、/>、/>Constitution ofGate (n-type) equivalent circuit representation, wherein->,/>,, wherein ,/>One end is grounded, the other end is connected with->Is connected with the current inflow end of (a)>One end is grounded, the other end is connected with->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->The current flowing into the line at the head end is +.>Flow through->The current of (2) is I, flow->Is +.>The current of the excitation part is +.>。
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 systemState 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 +.>Measurement noise covariance->;
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:
wherein ,discrete time series>Representing the state of the system parameter at time k +.>Representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>Is a state transition matrix, ">. For line parameters, state variables ∈ ->。
The measurement equation is established below, and equivalent reactance is set for the line parametersThe voltage at two ends is +>In the continuous time domain, according to kirchhoff's law:
Representing the voltage amplitude of the line current outflow terminal at the moment t, < >>Representing the line current inflow terminal voltage amplitude at time t, < >>Representing the equivalent reactance +.>Voltage at two ends>Representing the flow at time tCurrent through the line->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:
wherein ,are all continuous time functions, and are sampled (in discrete time domain) asAbbreviated as +.>,/>Is the sampling time. I.e.)>Represents the current at the line current inflow at time k, < >>Representing the current at the line current inflow at time k-1,represents the line current inflow terminal voltage amplitude at time k-1,/->Representing the line at time kCurrent inflow terminal voltage amplitude, ">Representing the voltage amplitude of the line current outflow end at the moment k;
substituting (4), (5) and (6) into (1) to obtain:
the preparation method comprises the following steps of:
t is transformer transformation ratio, impedanceJ represents the imaginary part of the complex number, < ->;
(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:
zero mean and covariance +.>Is a gaussian white noise of (c). Measurement value->,The voltage amplitude at the end of the line (line of the system), respectively +.>Is a measurement matrix of the measurement data,,/>representing the line current inflow terminal voltage amplitude at time k-2. />Zero mean and covariance +.>Is Gaussian white noise, < >>Is extra noise->To measure the wild value, wherein->Obeying the parameter +.>Is characterized by a Bernoulli distribution,/>,/>for a constant of larger amplitude, illustratively, < +.>Between 100 and 1000%>Is a unit pulse function>Obeying the parameter +.>Is characterized by a Bernoulli distribution,the initial value and covariance of the system parameters of the medium voltage distribution system are +.> and />。
Step 2: system parameter correction model (namely formula (10)) based on adaptive Kalman filtering and calculationPredicted value of system parameter of time medium voltage distribution system>And covariance + ->And the initial value of the state variable +.>And covariance +.>;
The predicted value and covariance calculation formula are as follows:
knowing the initial valueAnd covariance + ->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 +.>Predicted value of system parameter of time medium voltage distribution system>And covariance + ->. Initial value->And covariance + ->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 estimationAnd calculates the residual +.>Covariance->;
step 4: calculation of the mahalanobis distance from the resulting residualIf->Then go to step 5; otherwise, go to step 6;
utilizing residual errorsAnd covariance + ->The calculated mahalanobis distance is calculated as follows:
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 and />The size of (1)>The corresponding measurement value is regarded as uncertain noise, and the adaptive factor is calculated>And updating the prior residual, calculating the updated residual covariance +.>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 residualIf (if)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->Then it is considered that uncertainty noise is present at this time. Adaptive factor->And its updated residual covariance->The calculation formula of (2) is as follows: />
wherein ,、/>the lower and upper boundaries of the confidence interval are represented, illustratively by values of 0 and 1, respectively.Is the residual covariance before update, +.>For the updated residual covariance +.>For vector->Is>Element(s)>For vector->Is>Element(s)>、/>、For matrix->Diagonal->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;
Kalman filtering gain for medium voltage distribution system parametersThe calculation formula of (2) is as follows:
if uncertain noise occurs, then:
step 8: by means ofCalculated->Posterior state estimation value of time-of-day system parameters +.>And covariance + ->(matrix) and will->As a result of the correction of system parameters.
Medium voltage power distribution systemPosterior state estimation value of time-of-day system parameters +.>And covariance + ->The calculation formula of (2) is as follows:
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 respectivelyThe current flowing into the line is +.>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:
a state variable representing a parameter of the system,,/>representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>Is a state transition matrix, ">,/>Zero mean and covariance +.>Is white gaussian noise; />Representing the measured value->Is a measurement matrix of the measurement data,,/>represents the voltage amplitude of the line current inflow at time k,/->Represents the current at the line current inflow at time k, +.>Representing the current at the line current inflow at time k-1, < >>Representing the line current inflow terminal voltage amplitude at time k-1,representing the voltage amplitude of the line current inflow terminal at the time k-2; />Zero mean and covariance +.>Is Gaussian white noise, < >>Is extra noise->To measure the wild value, wherein->Obeying the parameter +.>Is characterized by a Bernoulli distribution,,/>is a unit pulse function>Obeying the parameter +.>Bernoulli distribution, ->。
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、/>、/>A gate-type equivalent circuit representation is constructed, wherein +.>,/>,/>, wherein ,/>One end is grounded, the other end is connected with->Is connected with the current inflow end of (a)>One end is grounded, the other end is connected with->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 +.>The current flowing into the line at the head end is +.>Flow through->The current of (2) is I, flow->Is +.>The current of the excitation part is +.>,/>Is connected with the current inflow end of the equivalent circuit structure of the transformer, +.>The current outflow end of the transformer equivalent circuit structure;
in the state equation and the measurement equation,
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 and />,
Calculated by the following formula:
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:
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;
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 noiseThe calculation formula of (2) is as follows:
if uncertain noise occurs, then:
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 respectivelyThe current flowing into the line is +.>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:
a state variable representing a parameter of the system,,/>representing the state of the system parameter at the previous instant, i.e. at time k-1 +.>Is a state transition matrix, ">,/>Zero mean and covariance +.>Is white gaussian noise; />Representing the measured value->Is a measurement matrix of the measurement data,,/>represents the voltage amplitude of the line current inflow at time k,/->Represents the current at the line current inflow at time k, +.>Representing the current at the line current inflow at time k-1, < >>Representing the line current inflow terminal voltage amplitude at time k-1,representing the voltage amplitude of the line current inflow terminal at the time k-2; />Zero mean and covariance +.>Is Gaussian white noise, < >>Is extra noise->To measure the wild value, wherein->Obeying the parameter +.>Is characterized by a Bernoulli distribution,,/>is a unit pulse function>Obeying the parameter +.>Bernoulli distribution, ->。/>
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)
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 |
-
2023
- 2023-02-01 CN CN202310048743.5A patent/CN115800271B/en active Active
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 |