CN115241901A - Intelligent energy storage soft switch data driving voltage control method considering data quality - Google Patents

Intelligent energy storage soft switch data driving voltage control method considering data quality Download PDF

Info

Publication number
CN115241901A
CN115241901A CN202210893877.2A CN202210893877A CN115241901A CN 115241901 A CN115241901 A CN 115241901A CN 202210893877 A CN202210893877 A CN 202210893877A CN 115241901 A CN115241901 A CN 115241901A
Authority
CN
China
Prior art keywords
energy storage
data
time
voltage
soft switch
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210893877.2A
Other languages
Chinese (zh)
Other versions
CN115241901B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN202210893877.2A priority Critical patent/CN115241901B/en
Publication of CN115241901A publication Critical patent/CN115241901A/en
Application granted granted Critical
Publication of CN115241901B publication Critical patent/CN115241901B/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
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • 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
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to an intelligent energy storage soft switch data driving voltage control method considering data quality, which comprehensively considers the agnostic property of power distribution network line parameters, the position of a distributed power supply and the uncertainty of output conditions, establishes a data driving model through measured data, and adopts a density-based local outlier factor method to identify bad data and introduce an attenuation factor to inhibit the measured disturbance in the data driving and adjusting process because the measured bad data and the measured data disturbance in uploaded data can influence the voltage adjusting effect of a data driving power distribution network. The invention provides a method for controlling the voltage of a power distribution network, which can overcome the adverse effects on data driving voltage control caused by measuring bad data and measuring disturbance in measured data under the condition of poor data quality, regulate the voltage of the power distribution network comprising a multi-terminal intelligent energy storage soft switch, ensure the effectiveness and stability of the data driving voltage control, effectively improve the operation controllability and flexibility of the power distribution network, and have important significance for ensuring the safe and economic operation of the power distribution network.

Description

Intelligent energy storage soft switch data driving voltage control method considering data quality
Technical Field
The invention relates to a voltage regulation method for a power distribution network. In particular to an intelligent energy storage soft switch data driving voltage control method considering data quality.
Background
Due to the high permeability access of the distributed power supply, the operation control of the power distribution network faces huge challenges including voltage out-of-limit, tidal current reverse transmission and the like. The multi-end intelligent energy storage soft switch is a novel flexible power electronic device replacing a traditional feeder line connection switch, can quickly and accurately control tide distribution to achieve flexible exchange of power between feeder lines in a larger range, combines the quick response characteristic of an energy storage link, and provides a better solution for dealing with the problem of safe and economical operation of a power grid of an intermittent power supply under the situation that the permeability of a distributed power supply is improved day by day. The application of the multi-terminal intelligent energy storage soft switch in the power distribution network provides quick response capacity for the power distribution network, helps to stabilize voltage fluctuation and adjust power distribution, and effectively improves operation controllability and flexibility of the power distribution network. Therefore, the reasonable regulation and control of the multi-terminal intelligent energy storage soft switch has important significance for guaranteeing the safe and economic operation of the power distribution network.
However, the multi-terminal intelligent energy storage soft switching optimization strategy based on the physical model cannot adapt to frequent changes of the running state of the power distribution network. And under the actual complicated operating environment, the accurate parameter of distribution network is difficult to obtain. The rapid development of the intelligent measurement terminal and the communication network promotes the high informatization of the power distribution system, and creates favorable conditions for the implementation of the data-driven optimization method. The data driving does not depend on detailed mathematical model information of a controlled system, the input and output relation of a complex link is described statistically by using measured data, the simulation construction of unknown characteristics of the complex link is realized, and the aim of controlling the voltage of the power distribution network is further realized. However, due to the limitations of the accuracy of the measurement device itself, measurement disturbance exists in the measurement data, and due to the appearance of external physical factors and other data sources, bad data also exists in the measurement data, so that the measurement upload data in the actual power distribution network has a data quality problem, while the data driving method completely depends on the measurement data, and the quality of the measurement data inevitably affects the voltage regulation effect of the data-driven flexible power distribution network.
Therefore, the bad data identification link is introduced to carry out online identification on the measured data, and attenuation factors are introduced aiming at the measurement disturbance, so that a data driving voltage control method considering the quality of the measured data is designed, the influence of the bad data and the measurement disturbance on the control effect is weakened, and the voltage regulation effect of the power distribution network is effectively improved.
Disclosure of Invention
The invention aims to solve the technical problem of providing a data quality-considered intelligent energy storage soft switch data driving voltage control method which can regulate the voltage of a power distribution network comprising a multi-terminal intelligent energy storage soft switch under the condition of poor data quality.
The technical scheme adopted by the invention is as follows: an intelligent energy storage soft switching data driving voltage control method considering data quality comprises the following steps:
1) Inputting node voltage, injection power historical data and parameter information of the active power distribution network at each moment of the system per day according to the selected active power distribution network;
2) Obtaining bad data and measurement disturbance w in a period from t-delta t' to t according to the active power distribution network selected in the step 1) t The voltage measurement data of each node is used for establishing a voltage measurement matrix at the time t
Figure BDA0003768626710000021
For is to
Figure BDA0003768626710000022
Identifying bad data from the voltage measurement data, and calculating bad numberAccording to the processed average measured data of the voltage at the time t
Figure BDA0003768626710000023
3) According to the average measured data of the voltage at the time t
Figure BDA0003768626710000024
Combined node voltage reference interval
Figure BDA0003768626710000025
Judging whether voltage out-of-limit occurs or not, if not, turning to the step 7), and if so, calculating the voltage deviation at the time t
Figure BDA0003768626710000026
And are aligned with
Figure BDA0003768626710000027
Carrying out measurement disturbance suppression processing to obtain voltage deviation disturbance suppression data at the time t
Figure BDA0003768626710000028
4) Suppressing data according to the voltage deviation disturbance at the time t
Figure BDA0003768626710000029
Establishing a pseudo Jacobian matrix of the r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the moment t
Figure BDA00037686267100000210
Combining voltage deviation disturbance suppression data at time t
Figure BDA00037686267100000211
And calculating a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at t moment by adopting a multi-layer hierarchical prediction algorithm with data of corresponding moments in active power distribution network node voltage and injection power historical data at m days
Figure BDA00037686267100000212
5) Pseudo Jacobian matrix of the r port of the multi-terminal intelligent energy storage soft switch SOP link according to the time t
Figure BDA00037686267100000213
And t moment multi-terminal intelligent energy storage soft switch energy storage link pseudo Jacobian estimation matrix
Figure BDA00037686267100000214
The data drive multi-terminal intelligent energy storage soft switch data drive voltage control model of considering measured data quality is established, includes: the method comprises the steps of taking the minimum voltage deviation of each node of an active power distribution network as a target function, restraining capacity and loss of converters at each port of a multi-terminal intelligent energy storage soft switch SOP link, restraining upper and lower limits of output of active power and reactive power, restraining charge and discharge efficiency, restraining upper and lower limits of charge and discharge power and restraining a prediction domain T of a multi-terminal intelligent energy storage soft switch energy storage link p Internal charge state constraint and charge-discharge frequency limit constraint;
6) Adopting a mathematical solver to solve a data-driven multi-terminal intelligent energy storage soft switch data driving voltage control model considering the quality of measured data, issuing and executing a solution result, wherein the solution result comprises: charging and discharging power of multi-terminal intelligent energy storage soft switch energy storage link at time t
Figure BDA00037686267100000215
Active and reactive power output of the r-th port of multi-terminal intelligent energy storage soft switch SOP link at time t
Figure BDA00037686267100000216
Figure BDA00037686267100000217
7) Update control time t = t + Δ t, and determine t-t 0 ≥βΔT c If yes, let T p =T p -ΔT c β = β +1, if no, step 8) is performed, where β is the control parameter, Δ t is the control time interval, t 0 Is a starting time, T p To predict the field, Δ T c Is a control domain;
8) According to the control time t in the step 7), judging t-t 0 And (4) whether the time is longer than the optimization time T or not, if not, turning to the step 2), and if so, ending.
According to the intelligent energy storage soft switch data driving voltage control method considering the data quality, the unknown parameters of the power distribution network line, the distributed power supply position and the uncertainty of the output condition are comprehensively considered, the data driving model is established through the measured data, and the measured bad data and the measured data disturbance existing in the uploaded data influence the voltage regulation effect of the data driving power distribution network, so that the bad data are identified by adopting a density-based local outlier factor method, and the measured disturbance is restrained by introducing an attenuation factor in the data driving regulation process. According to the method, under the conditions of actual complex operation environment, frequent change of operation state of the power distribution network and difficulty in obtaining accurate parameters of the power distribution network, adverse effects on data driving voltage control caused by bad data measurement and disturbance in measurement data can be overcome, voltage regulation is performed on the power distribution network with the multi-terminal intelligent energy storage soft switch, effectiveness and stability of data driving voltage control are guaranteed, operation controllability and flexibility of the power distribution network are effectively improved, and the method has important significance for guaranteeing safe and economical operation of the power distribution network.
Drawings
FIG. 1 is a flow chart of the intelligent energy storage soft switch data driving voltage control method considering data quality according to the present invention;
FIG. 2 is a schematic diagram of an example of a flexible interconnection power distribution network with four-terminal intelligent energy storage soft switches in an improved Tianjin demonstration area;
FIG. 3 is a 24 hour distributed power output curve;
FIG. 4 is a 24 hour load curve;
FIG. 5 is a comparison graph of the global voltage maximum results of each scene in 24 hours;
FIG. 6 is a comparison graph of the global voltage minimum results for each scene at 24 hours;
FIG. 7 is a graph of the maximum node voltage distribution for each region in 24 hours;
FIG. 8 is a 24-hour multi-terminal intelligent energy storage soft switch power output diagram of each port;
fig. 9 is a reactive force diagram of each port of the 24-hour multi-port intelligent energy storage soft switch;
fig. 10 is a state of charge change diagram of the 24-hour energy storage link.
Detailed Description
The following describes the intelligent energy storage soft switching data driving voltage control method considering data quality in detail with reference to the embodiments and the accompanying drawings.
As shown in fig. 1, the intelligent energy storage soft switching data driving voltage control method considering data quality of the present invention includes the following steps:
1) Inputting node voltage, injection power historical data and parameter information of the active power distribution network at each time of day of the system according to the selected active power distribution network; the parameter information comprises: capacity, access position, initial charge state of energy storage link, charge-discharge power limit value and node voltage reference interval of multi-terminal intelligent energy storage soft switch
Figure BDA0003768626710000031
Prediction domain T p Control field Δ T c Control time interval Δ T, optimization time T, measured data acquisition time window Δ T', initialization control parameter β =1, start time T 0 Initialization control time t = t 0
2) Obtaining bad data and measurement disturbance w in a period from t-delta t' to t according to the active power distribution network selected in the step 1) t The voltage measurement data of each node is used for establishing a voltage measurement matrix at the time t
Figure BDA0003768626710000032
To pair
Figure BDA0003768626710000033
Identifying bad data by the voltage measurement data in the process, and calculating the average measurement data of the voltage at the time t after the bad data is processed
Figure BDA0003768626710000034
The method specifically comprises the following steps:
(2.1) establishing a voltage measurement matrix at the time t
Figure BDA0003768626710000035
Obtaining bad data and measurement disturbance w of each node in the period from t-delta t' to t t The voltage measurement data of the time point t, and a voltage measurement matrix at the time point t are established
Figure BDA0003768626710000036
The following were used:
Figure BDA0003768626710000037
U Δt′,ξ =[U Δt′,ξ (1),…,U Δt′,ξ (i),…,U Δt′,ξ (I)] T (2)
wherein z represents the voltage measurement data group number uploaded by each node in the period from t-delta t 'to t, I represents the node number of the active power distribution network, and U represents the voltage measurement data group number uploaded by each node in the period from t-delta t' to t Δt′,1 、U Δt′,ξ 、U Δt′,z Respectively representing the voltage measurement data of the 1 st, xi and z groups, U, uploaded by each node in the period from t-delta t' to t Δt′,ξ (1)、U Δt′,ξ (i)、U Δt′,ξ (I) Respectively representing data uploaded by 1 st, I th and I nodes of the active power distribution network in a xi group in a period from t-delta t' to t;
defining a voltage measurement matrix u t The two-dimensional coordinates of the medium voltage data are as follows:
[U Δt′,ξ (i),U Δt′,ξ (i)/σ norm (U Δt′,ξ )] (3)
in the formula, σ norm (U Δt′,ξ ) Indicating voltage measuring matrix
Figure BDA0003768626710000038
Standard deviation of group ξ data.
(2.2) determining the kth distance neighborhood of each voltage measurement data
Voltage measurement data U uploaded by xth node of active power distribution network in xi group in time period from t-delta t' to t Δt′,ξ (x) A k-distance neighborhood centered at N k [U Δt′,ξ (x)]Expressed as follows:
Figure BDA0003768626710000041
in the formula, o is a node.
(2.3) calculating the voltage measurement data U uploaded by the xth node Δt′,ξ (x) Distance d between the voltage measurement data uploaded by other nodes k [U Δt′,ξ (x),U Δt′,ξ (y)]Expressed as follows:
Figure BDA0003768626710000042
in the formula, d [ U ] Δt′,ξ (x),U Δt′,ξ (y)]Voltage measurement data U uploaded by the xth node of the active power distribution network in the xi group in the period from t-delta t' to t Δt′,ξ (x) The voltage measurement data U uploaded to the y node Δt′,ξ Euclidean distance of (y), d k [U Δt′,ξ (x)]Is represented in N k [U Δt′,ξ (x)]Internal and voltage measurement data U Δt′,ξ (x) The euclidean distance from the farthest data point, t, is any node in the k distances from node x.
(2.4) calculating local outlier factor of each voltage measurement data
Voltage measurement data U Δt′,ξ (x) Local outlier factor F of k [U Δt′,ξ (x)]Expressed as follows:
Figure BDA0003768626710000043
in the formula, ρ k [U Δt′,ξ (x)]、ρ k [U Δt′,ξ (y)]Respectively representVoltage measurement data U Δt′,ξ (x) And U Δt′,ξ (y) local density, expressed as follows:
Figure BDA0003768626710000044
(2.5) calculating a local outlier factor threshold
Figure BDA0003768626710000045
If F k [U Δt′,ξ (i)]>F k,max Then voltage measurement data U Δt′,ξ (i) For bad data, let U Δt′,ξ (i)=0
In the formula, F k,max Representing a local outlier factor threshold, Q, set using four-quadrant analysis 1 、Q 3 Representing the local outlier set F k Lower quartile and upper quartile, F k =[F k [U Δt′,ξ (1)],…,F k [U Δt′,ξ (i)],…,F k [U Δt′,ξ (I)]],F k [U Δt′,ξ (1)]、[U Δt′,ξ (i)]、F k [U Δt′,ξ (I)]Respectively represent voltage measurement data U Δt′,ξ (1)、U Δt′,ξ (i)、U Δt′,ξ (I) The local outlier factor of (a) is,
Figure BDA0003768626710000046
representing a local outlier set F k Standard deviation of [, ] 1 、η 2 And represents the standard deviation discrimination coefficient.
(2.6) calculating average measured data of voltage at t moment after bad data processing
Figure BDA0003768626710000047
The new voltage measurement matrix after the identification processing of the measured bad data is
Figure BDA0003768626710000048
Figure BDA0003768626710000049
Respectively representing voltage measurement matrix
Figure BDA00037686267100000410
Group
1, xi and z data;
according to a new voltage measurement matrix of
Figure BDA00037686267100000411
Calculating average measurement data of voltage at t moment
Figure BDA00037686267100000412
The following:
Figure BDA0003768626710000051
3) According to the average measured data of the voltage at the time t
Figure BDA0003768626710000052
Combined node voltage reference interval
Figure BDA0003768626710000053
Judging whether voltage out-of-limit occurs or not, if not, turning to the step 7), and if so, calculating the voltage deviation at the time t
Figure BDA0003768626710000054
And to
Figure BDA0003768626710000055
Carrying out measurement disturbance suppression processing to obtain voltage deviation disturbance suppression data at the time t
Figure BDA0003768626710000056
Said pair
Figure BDA0003768626710000057
Carrying out measurement disturbance suppression processing to obtain voltage deviation disturbance suppression data at the time t
Figure BDA0003768626710000058
Is represented as follows:
Figure BDA0003768626710000059
in the formula, e w,t Indicating a set node voltage reference value U ref Voltage average measurement data with time t
Figure BDA00037686267100000510
The vector of the difference values of (a),
Figure BDA00037686267100000511
mean measurement data of voltage at t
Figure BDA00037686267100000512
Voltage average measurement data at time t-delta t
Figure BDA00037686267100000513
Zeta denotes the attenuation factor, N k =ΔT c At, denotes the control field Δ T c Regulation and control times of internal multi-terminal intelligent energy storage soft switch, e s E represents a unit column vector as a deviation threshold parameter; n is a radical of w Is expressed at DeltaT c Within, from | e w,t |≥e s E phase switch to | E w,t |<e s The number of times that the multi-terminal intelligent energy storage soft switch has completed regulation and control in the stage E; β represents the attenuation factor amplification factor.
4) Suppressing data according to the voltage deviation disturbance at the time t
Figure BDA00037686267100000514
Establishing a pseudo Jacobian matrix of an r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the time t
Figure BDA00037686267100000515
Combining voltage deviation disturbance suppression data at time t
Figure BDA00037686267100000516
And calculating a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the t moment by adopting a multi-layer hierarchical prediction algorithm with data of corresponding moments in historical data of node voltage and injection power of the active power distribution network at each moment of m days
Figure BDA00037686267100000517
Wherein,
the suppression data according to the voltage deviation disturbance at the time t
Figure BDA00037686267100000518
Establishing a pseudo Jacobian matrix of the r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the moment t
Figure BDA00037686267100000519
Is represented as follows:
Figure BDA00037686267100000520
if it is
Figure BDA00037686267100000521
Or
Figure BDA00037686267100000522
Or
Figure BDA00037686267100000523
Then the
Figure BDA00037686267100000524
In the formula,
Figure BDA00037686267100000525
and
Figure BDA00037686267100000526
respectively representing pseudo Jacobian matrixes of an r-th port of an intelligent energy storage soft switch SOP link at the time t and the time t-delta t,
Figure BDA00037686267100000527
is composed of
Figure BDA00037686267100000528
An initial value of (1);
Figure BDA00037686267100000529
Figure BDA00037686267100000530
the variable quantity vector of the active control quantity of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time of t-delta t is shown,
Figure BDA00037686267100000531
Figure BDA00037686267100000532
Figure BDA00037686267100000533
respectively represents the output strategies of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time t-delta t and the time t-2 delta t,
Figure BDA00037686267100000534
respectively representing active power difference vectors and reactive power difference vectors of the r-th port of the multi-port intelligent energy storage soft switch SOP link at the time of t-delta t,
Figure BDA00037686267100000535
η SOP 、μ SoP representing the weight coefficients.
The combined t-time voltage deviation disturbance suppression data
Figure BDA00037686267100000536
Data corresponding to the node voltage and injection power of the active power distribution network at each moment of m days in historical dataCalculating a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the time t by adopting a multi-layer hierarchical prediction algorithm
Figure BDA0003768626710000061
Is represented as follows:
Figure BDA0003768626710000062
in the formula,
Figure BDA0003768626710000063
and (3) representing a pseudo Jacobian matrix of the energy storage link at the time t, wherein the iterative computation solving process comprises the following steps:
Figure BDA0003768626710000064
if it is
Figure BDA0003768626710000065
Or
Figure BDA0003768626710000066
Or
Figure BDA0003768626710000067
Then the
Figure BDA0003768626710000068
In the formula,
Figure BDA0003768626710000069
Figure BDA00037686267100000610
representing the difference vector of the charging and discharging power of the multi-terminal intelligent energy storage soft switch energy storage link at the time of t-delta t,
Figure BDA00037686267100000611
respectively represents the energy storage link at the t-delta t moment and the t-2 delta t momentThe active power output strategy of (a) is,
Figure BDA00037686267100000612
is composed of
Figure BDA00037686267100000613
Of initial value of ES 、μ ES Represents a weight coefficient, and epsilon represents a threshold coefficient;
in the formula (12), the reaction mixture is,
Figure BDA00037686267100000614
represents T + n Δ T c The pseudo jacobian prediction matrix at a time, N ∈ {1, \8230;, N }, N = T p /ΔT c Representing the prediction domain T p And control field Δ T c The calculation formula is as follows:
Figure BDA00037686267100000615
in the formula, theta α,t Denotes the prediction coefficient at time T, where α =1, \8230;, l, T d Representing an estimated time duration in days;
defining a regression coefficient vector θ t =(θ 1,t ,…,θ α,t ,…,θ l,t ) T ,θ 1,t 、θ α,t 、θ l,t Respectively represent time t theta t 1, α, l element, θ t The iterative solution formula is as follows:
Figure BDA00037686267100000616
in the formula, a pseudo-elegant predictive matrix set of t-delta t
Figure BDA00037686267100000617
Figure BDA00037686267100000618
And
Figure BDA00037686267100000619
respectively, using T + n Δ T c -T d Time sum T + n Δ T c -lT d A pseudo Jacobian matrix is obtained by calculating historical measurement data of the active power distribution network at any moment, delta represents a weight coefficient, and theta t-Δt Is the regression coefficient vector at the time t-delta t.
5) Pseudo Jacobian matrix of the r port of the multi-terminal intelligent energy storage soft switch SOP link according to the time t
Figure BDA00037686267100000620
And t moment multi-terminal intelligent energy storage soft switch energy storage link pseudo Jacobian estimation matrix
Figure BDA00037686267100000621
The data drive multi-terminal intelligent energy storage soft switch data drive voltage control model of the quality of the measured data is considered to be established, and the data drive multi-terminal intelligent energy storage soft switch data drive voltage control model comprises the following steps: the method comprises the steps of taking the minimum voltage deviation of each node of an active power distribution network as a target function, restraining capacity and loss of converters at each port of a multi-terminal intelligent energy storage soft switch SOP link, restraining upper and lower limits of output of active power and reactive power, restraining charge and discharge efficiency, restraining upper and lower limits of charge and discharge power and restraining a prediction domain T of a multi-terminal intelligent energy storage soft switch energy storage link p Internal charge state constraint and charge-discharge frequency limit constraint; wherein,
the objective function is expressed as follows:
Figure BDA00037686267100000622
Figure BDA0003768626710000071
Figure BDA0003768626710000072
Figure BDA0003768626710000073
wherein J represents an objective function, J 1 、J′ 1 The adjusted objective function is represented as a function of the target,
Figure BDA0003768626710000074
representing the voltage reference at time t + deltat,
Figure BDA0003768626710000075
represents the average measured data of the voltage at the time t,
Figure BDA0003768626710000076
respectively an upper limit and a lower limit of a voltage reference interval,
Figure BDA0003768626710000077
Figure BDA0003768626710000078
r represents the port number of the multi-terminal intelligent energy storage soft switch,
Figure BDA0003768626710000079
the difference vector of the charging and discharging power of the multi-terminal intelligent energy storage soft switch energy storage link at the moment t is shown,
Figure BDA00037686267100000710
respectively represents the charging and discharging power of the energy storage link at the time t and the time t-delta t,
Figure BDA00037686267100000711
the variable quantity vector of the active and reactive control quantity of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time t is shown,
Figure BDA00037686267100000712
respectively representing active power difference value vector and reactive power difference value vector of the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the moment t, lambda ES 、λ SOP Represents a weight coefficient, ζ represents an attenuation factor;
Figure BDA00037686267100000713
the estimated value of the voltage of each node at the time t + Δ t is expressed as follows:
Figure BDA00037686267100000714
in the formula,
Figure BDA00037686267100000715
representing the vector of the voltage regulation and control variation quantity of the SOP link of the multi-terminal intelligent energy storage soft switch at the time t,
Figure BDA00037686267100000716
a pseudo Jacobian matrix of an r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the time t is represented;
Figure BDA00037686267100000717
denotes from time T to T + N Δ T c The multi-terminal intelligent energy storage soft switch energy storage link voltage regulation variation vector at a moment,
Figure BDA00037686267100000718
respectively represents time T, T + (N-1) DeltaT c Constantly multi-terminal intelligent energy storage soft switch energy storage link voltage regulation and control variable vectors;
Figure BDA00037686267100000719
expressing a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the time t;
Figure BDA00037686267100000720
denotes the time from T to T + (N-1) Δ T c The energy storage link charges and discharges the power variable quantity vector at all times,
Figure BDA00037686267100000721
Figure BDA00037686267100000722
respectively represent time t and t +(N-1)ΔT c And (3) charging and discharging power difference vectors of the moment energy storage link.
The capacity of each port converter, the loss constraint of each port converter and the active and reactive power output upper and lower limit constraints of each port converter in the multi-terminal intelligent energy storage soft switch SOP link are expressed as follows:
Figure BDA00037686267100000723
Figure BDA0003768626710000081
Figure BDA0003768626710000082
in the formula,
Figure BDA0003768626710000083
the actual charging and discharging power of the energy storage link is shown,
Figure BDA0003768626710000084
respectively representing active and reactive power output of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time t,
Figure BDA0003768626710000085
the loss of the current converter at the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the time t is shown,
Figure BDA0003768626710000086
the loss coefficient of the current converter at the r-th port of the SOP link of the multi-terminal intelligent energy storage soft switch is shown,
Figure BDA0003768626710000087
respectively represents the upper limit and the lower limit of active power of a multi-terminal intelligent energy storage soft switch SOP link,
Figure BDA0003768626710000088
the actual injection active power of the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the moment t is shown,
Figure BDA0003768626710000089
respectively representing the upper and lower limits of reactive power of a multi-terminal intelligent energy storage soft switch SOP link at the time t,
Figure BDA00037686267100000810
and the capacity of the r-th port converter of the multi-terminal intelligent energy storage soft switch SOP link is shown.
The charge-discharge efficiency constraint, the charge-discharge power upper and lower limit constraint and the prediction domain T of the multi-terminal intelligent energy storage soft switch energy storage link are adopted p The internal state of charge constraint and the charge-discharge number limit constraint are expressed as follows:
Figure BDA00037686267100000811
if it is
Figure BDA00037686267100000812
And is
Figure BDA00037686267100000813
Figure BDA00037686267100000814
In the formula,
Figure BDA00037686267100000815
the actual charging and discharging power of the energy storage link is shown,
Figure BDA00037686267100000816
respectively represents the charge and discharge power of the energy storage link at the time t and the time t-delta t,
Figure BDA00037686267100000817
respectively representing energy-storage ringsThe charge-discharge efficiency of the node is improved,
Figure BDA00037686267100000818
respectively represents the upper limit and the lower limit of the calculated value of the charging and discharging power of the energy storage link,
Figure BDA00037686267100000819
respectively represent T + DeltaT c Time sum T + (N-1) Δ T c Charging and discharging power difference of moment energy storage link, S ES The capacity of the energy storage link is shown,
Figure BDA00037686267100000820
representing the upper and lower limits of the charge state of the energy storage link,
Figure BDA00037686267100000821
respectively representing the initial state and the prediction domain T of the energy storage system p The state of charge after the end of the cycle,
Figure BDA00037686267100000822
which is indicative of a threshold coefficient of the value,
Figure BDA00037686267100000823
respectively representing the charge states of the multi-terminal intelligent energy storage soft switch energy storage link at T and T-delta T, wherein N = T p /ΔT c Representing the prediction domain T p And control field Δ T c The ratio of (a) to (b),
Figure BDA00037686267100000824
and the charging and discharging power difference of the multi-terminal intelligent energy storage soft switch energy storage link at the moment t is represented.
6) Adopting a mathematical solver to solve a data-driven multi-terminal intelligent energy storage soft switch data driving voltage control model considering the quality of measured data, issuing and executing a solution result, wherein the solution result comprises: charging and discharging power of multi-terminal intelligent energy storage soft switch energy storage link at time t
Figure BDA0003768626710000091
T moment multi-terminal intelligent energy storage soft switch SOP linkActive and reactive power output of the r-th port
Figure BDA0003768626710000092
Figure BDA0003768626710000093
7) Updating control time t = t + Δ t, and judging t-t 0 ≥βΔT c If true, let T p =T p -ΔT c β = β +1, if no, step 8) is performed, where β is the control parameter, Δ t is the control time interval, t is 0 Is a starting time, T p To predict the field, Δ T c Is a control domain;
8) According to the control time t in the step 7), judging t-t 0 And (3) whether the time length is greater than the optimization time length T, if not, turning to the step 2), and if so, ending.
Examples are given below:
the embodiment of the invention is an improved Tianjin demonstration area containing four-terminal intelligent energy storage soft switch flexible interconnection power distribution network example, and the topology is shown in figure 2. The detailed parameters are shown in tables 1-2.
TABLE 1 Tianjin demonstration area distribution network example load access position and power
Figure BDA0003768626710000094
TABLE 2 example line parameters for distribution network in Tianjin demonstration area
Figure BDA0003768626710000101
In the demonstration area, two 110kV transformer substations serve as centers to form a double-loop network structure comprising four feeders, four distribution network areas are flexibly interconnected through four-end intelligent energy storage soft switches, the voltage grades are set to be 10.5kV, and the total active power demand and the total reactive power demand of the load are respectively 9.9880MW and 7.3350Mvar.
In the embodiment, in order to fully consider the influence of the access of a high-permeability distributed power supply on the operation loss and the voltage fluctuation of a flexible power distribution network, 3 groups of photovoltaic systems and 3 wind turbine generators are respectively accessed in an example of a flexible interconnection power distribution network with four-terminal intelligent energy storage soft switches in a demonstration area, and the permeability of the distributed power supply reaches 90.11%. The access positions of the photovoltaic system are nodes 8, 37 and 41, and the capacity is 1MVA; the access positions of the wind turbine generator are nodes 12, 21 and 44, the capacity is 2MVA, and the output reduction of the distributed power supply is not allowed. The daily active power output curve of the distributed power supply is shown in fig. 3. In addition, in the calculation example, the capacity of a current converter of each port of the four-terminal intelligent energy storage soft switch is 3MVA, the capacity of an energy storage link is 2MWh, the limit of charge and discharge power is 0.5MW, and the charge state range is 20% -80%.
Setting each weight factor as lambda according to the convergence and control requirements of the model-free adaptive control method ES =20、μ Es =0.5、λ SOP =0.1、μ Es =10,η Es 、η SOP 、ρ SOP Delta is set to be 1, and the charge and discharge efficiency of the energy storage link is improved
Figure BDA0003768626710000111
Are all set to 99 percent, the loss coefficient of the converter
Figure BDA0003768626710000112
Threshold coefficient e =0.001, optimization time T set to 24 hours, prediction domain T p Set to 24 hours, control field Δ T c Set to 1 hour, the control time interval Δ t is set to 5 minutes, and the metrology data acquisition time window Δ t' is set to 0.5 minutes.
In order to fully verify the advancement of the intelligent energy storage soft switching data driving voltage control method considering the quality of measured data, in the embodiment, the following four scenes are adopted for comparative analysis:
scene 1: the multi-terminal intelligent energy storage soft switch is not controlled, and the initial running state of the power distribution network is obtained;
scene 2: the multi-terminal intelligent energy storage soft switch is subjected to data drive control, and links of bad data processing and measurement disturbance suppression are not included;
scene 3: the method is adopted to carry out data drive control on the multi-terminal intelligent energy storage soft switch, and comprises links of bad data processing and measurement disturbance suppression;
scene 4: and controlling the multi-terminal intelligent energy storage soft switch by adopting a centralized optimization method.
The computer hardware environment for executing the optimized calculation is Intel (R) Core (TM) i7-10700, the dominant frequency is 2.90GHz, and the internal memory is 24GB; the software environment is a Windows10 operating system.
The 24 hour global voltage maximum result pairs for scenarios 1, 2, 3, 4 are shown in fig. 5. The 24 hour global voltage minimum result pairs for scenarios 1, 2, 3, 4 are shown in fig. 6. Fig. 7 shows the maximum node voltage distribution in each region for 24 hours in scenarios 1, 2, 3, and 4. The active output of each port of the 24-hour multi-port intelligent energy storage soft switch in the scene 3 is shown in fig. 8. Reactive output of each port of the 24-hour multi-port intelligent energy storage soft switch in the scene 3 is shown in fig. 9. The state of charge change of the 24-hour energy storage link in the scenario 3 is shown in fig. 10. Table 3 shows voltage index ratio pairs in scenarios 1, 2, 3, and 4.
TABLE 3 comparison of Voltage indicators for various scenes
Figure BDA0003768626710000113
As can be seen from fig. 4 to 6, the scenario 3 can effectively adjust the voltage level of the power distribution network of the present embodiment. In combination with table 3, it can be seen by comparing the global voltage distributions of scenes 2 and 3 in fig. 4 and 5 that, after the bad data processing link and the measurement disturbance suppression link are added in the scene 3, the voltage level of the power distribution network is significantly increased compared with the scene 2, and it can be seen by comparing the global voltage distributions of the scenes 3 and 4 in fig. 4 and 5 that the voltage control effect of the data-driven voltage control scene considering the measured data quality is close to global optimum. It can be seen from fig. 4 to 9 that the intelligent energy storage soft switching data driving voltage control method considering the measured data quality can effectively solve the problem of optimal adjustment of the power distribution network voltage under the condition of poor data quality.

Claims (9)

1. An intelligent energy storage soft switching data driving voltage control method considering data quality is characterized by comprising the following steps:
1) Inputting node voltage, injection power historical data and parameter information of the active power distribution network at each time of day of the system according to the selected active power distribution network;
2) Obtaining bad data and measurement disturbance w in a period from t-delta t' to t according to the active power distribution network selected in the step 1) t The voltage measurement data of each node is used for establishing a voltage measurement matrix at the time t
Figure FDA0003768626700000011
To pair
Figure FDA0003768626700000012
Identifying bad data by the voltage measurement data in the process, and calculating the average measurement data of the voltage at the time t after the bad data is processed
Figure FDA0003768626700000013
3) According to the average measured data of the voltage at the time t
Figure FDA0003768626700000014
Combined node voltage reference interval
Figure FDA0003768626700000015
Judging whether voltage out-of-limit occurs or not, if not, turning to the step 7), and if so, calculating the voltage deviation at the time t
Figure FDA0003768626700000016
And to
Figure FDA0003768626700000017
Carrying out measurement disturbance suppression processing to obtain voltage deviation disturbance suppression data at the time t
Figure FDA0003768626700000018
4) Suppressing data according to the voltage deviation disturbance at the time t
Figure FDA0003768626700000019
Establishing a pseudo Jacobian matrix of the r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the moment t
Figure FDA00037686267000000110
Combined with t time voltage deviation disturbance suppression data
Figure FDA00037686267000000111
And calculating a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the t moment by adopting a multi-layer hierarchical prediction algorithm with data of corresponding moments in historical data of node voltage and injection power of the active power distribution network at each moment of m days
Figure FDA00037686267000000112
5) Pseudo Jacobian matrix of the r port of the multi-terminal intelligent energy storage soft switch SOP link according to the time t
Figure FDA00037686267000000113
And t moment multi-terminal intelligent energy storage soft switch energy storage link pseudo Jacobian estimation matrix
Figure FDA00037686267000000114
The data drive multi-terminal intelligent energy storage soft switch data drive voltage control model of the quality of the measured data is considered to be established, and the data drive multi-terminal intelligent energy storage soft switch data drive voltage control model comprises the following steps: the method comprises the steps of taking the minimum voltage deviation of each node of an active power distribution network as a target function, restraining capacity and loss of converters at each port of a multi-terminal intelligent energy storage soft switch SOP link, restraining upper and lower limits of output of active power and reactive power, restraining charge and discharge efficiency, restraining upper and lower limits of charge and discharge power and restraining a prediction domain T of a multi-terminal intelligent energy storage soft switch energy storage link p Internal charge state constraint and charge-discharge frequency limit constraint;
6) Adopting a mathematical solver to solve a data-driven multi-terminal intelligent energy storage soft switch data driving voltage control model considering the quality of measured data, issuing and executing a solution result, wherein the solution result comprises: charging and discharging power of multi-terminal intelligent energy storage soft switch energy storage link at time t
Figure FDA00037686267000000115
Active and reactive power output of the r-th port of multi-terminal intelligent energy storage soft switch SOP link at time t
Figure FDA00037686267000000116
7) Update control time t = t + Δ t, and determine t-t 0 ≥βΔT c If true, let T p =T p -ΔT c β = β +1, if no, step 8) is performed, where β is the control parameter, Δ t is the control time interval, t is 0 Is a starting time, T p To predict the field, Δ T c Is a control domain;
8) According to the control time t in the step 7), judging t-t 0 And (4) whether the time is longer than the optimization time T or not, if not, turning to the step 2), and if so, ending.
2. The method for controlling the data driving voltage of the intelligent energy storage soft switch considering the data quality as claimed in claim 1, wherein the parameter information in step 1) comprises: capacity, access position, initial charge state of energy storage link, charge-discharge power limit value and node voltage reference interval of multi-terminal intelligent energy storage soft switch
Figure FDA00037686267000000117
Prediction domain T p Control field Δ T c Control time interval Δ T, optimization time T, measured data acquisition time window Δ T', initialization control parameter β =1, start time T 0 Initialization control time t = t 0
3. The intelligent energy storage soft switching data driving voltage control method considering data quality according to claim 1, wherein the step 2) specifically comprises:
(2.1) establishing a voltage measurement matrix at the time t
Figure FDA0003768626700000021
Obtaining bad data and measurement disturbance w of each node in the period from t-delta t' to t t The voltage measurement data of the time point t, and a voltage measurement matrix at the time point t are established
Figure FDA0003768626700000022
The following were used:
Figure FDA0003768626700000023
U Δt′,ξ =[U Δt′,ξ (1),…,U Δt′,ξ (i),…,U Δt′,ξ (I)] T (2)
wherein z represents the voltage measurement data group number uploaded by each node in the period from t-delta t 'to t, I represents the node number of the active power distribution network, and U represents the voltage measurement data group number uploaded by each node in the period from t-delta t' to t Δt′,1 、U Δt′,ξ 、U Δt′,z Respectively representing the voltage measurement data of the 1 st, xi and z groups, U, uploaded by each node in the period from t-delta t' to t Δt′,ξ (1)、U Δt′,ξ (i)、U Δt′,ξ (I) Respectively representing data uploaded by 1 st, I th and I nodes of the active power distribution network in the xi group in a period from t-delta t' to t;
defined voltage measurement matrix
Figure FDA0003768626700000024
The two-dimensional coordinates of the medium voltage data are as follows:
[U Δt′,ξ (i),U Δt′,ξ (i)/σ norm (U Δt′ , ξ )] (3)
in the formula, σ norm (U Δt′,ξ ) Indicating voltage measurement matrix
Figure FDA0003768626700000025
Standard deviation of group xi data;
(2.2) determining the k-th distance neighborhood of each voltage measurement data
Voltage measurement data U uploaded by xth node of active power distribution network in the xi group in the time period from t-delta t' to t Δt′,ξ (x) A k-distance neighborhood centered at N k [U Δt′,ξ (x)]Expressed as follows:
Figure FDA0003768626700000026
in the formula, o is a node;
(2.3) calculating the voltage measurement data U uploaded by the xth node Δt′,ξ (x) Distance d between the voltage measurement data uploaded by other nodes k [U Δt′,ξ (x),U Δt′,ξ (y)]Expressed as follows:
Figure FDA0003768626700000027
in the formula, d [ U ] Δt′,ξ (x),U Δt′,ξ (y)]Voltage measurement data U uploaded by the xth node of the active power distribution network in the xi group in the period from t-delta t' to t Δt′,ξ (x) The voltage measurement data U uploaded to the y node Δt′,ξ Euclidean distance of (y), d k [U Δt′,ξ (x)]Is represented in N k [U Δt′,ξ (x)]Internal and voltage measurement data U Δt′,ξ (x) The Euclidean distance from the farthest data point, y is any node in the k distances from node x;
(2.4) calculating local outlier factors of each voltage measurement data
Voltage measurement data U Δt′,ξ (x) Local outlier factor F of k [U Δt′,ξ (x)]Expressed as follows:
Figure FDA0003768626700000031
in the formula, ρ k [U Δt′,ξ (x)]、ρ k [U Δt′,ξ (y)]Respectively represent voltage measurement data U Δt′,ξ (x) And U Δt′,ξ (y) local density, expressed as follows:
Figure FDA0003768626700000032
(2.5) calculating local outlier threshold
Figure FDA0003768626700000033
If F k [U Δt′,ξ (i)]>F k,max Then voltage measurement data U Δt′,ξ (i) For bad data, order U Δt′,ξ (i)=0
In the formula, F k,max Representing a local outlier factor threshold, Q, set using four-quadrant analysis 1 、Q 3 Representing a local outlier set F k Lower quartile and upper quartile, F k =[F k [U Δt′,ξ (1)],…,F k [U Δt′,ξ (i)],…,F k [U Δt′,ξ (I)]],F k [U Δt′,ξ (1)]、[U Δt′,ξ (i)]、F k [U Δt′,ξ (I)]Respectively represent voltage measurement data U Δt′,ξ (1)、U Δt′,ξ (i)、U Δt′,ξ (I) The local outlier factor of (a) is,
Figure FDA0003768626700000034
representing the local outlier set F k Standard deviation of [, ] 1 、η 2 Representing a standard deviation discrimination coefficient;
(2.6) calculating average measured data of voltage at t moment after bad data processing
Figure FDA0003768626700000035
The new voltage measurement matrix after the identification processing of the measured bad data is
Figure FDA0003768626700000036
Figure FDA0003768626700000037
Respectively represent voltage measurement matrix
Figure FDA0003768626700000038
Group 1, xi and z data;
according to a new voltage measurement matrix of
Figure FDA0003768626700000039
Calculating average measurement data of voltage at t moment
Figure FDA00037686267000000310
The following were used:
Figure FDA00037686267000000311
4. the intelligent energy storage soft-switching data driving voltage control method considering data quality as claimed in claim 1, wherein the pair in step 3)
Figure FDA00037686267000000312
Carrying out measurement disturbance suppression processing to obtain voltage deviation disturbance suppression data at the time t
Figure FDA00037686267000000313
Is represented as follows:
Figure FDA00037686267000000314
in the formula, e w,t Indicating a set node voltage reference value U ref Voltage average measurement data with time t
Figure FDA00037686267000000315
The vector of the difference values of (a),
Figure FDA00037686267000000316
mean measurement data representing voltage at time t
Figure FDA00037686267000000317
Voltage average measurement data at time t-delta t
Figure FDA00037686267000000318
Zeta denotes the attenuation factor, N k =ΔT c At, denotes the control field Δ T c Regulation and control times of internal multi-terminal intelligent energy storage soft switch, e s E represents a unit column vector as a deviation threshold parameter; n is a radical of w Is expressed at DeltaT c Within, from | e w,t |≥e s E phase switch to | E w,t |<e s The number of times that the multi-terminal intelligent energy storage soft switch has completed regulation and control in the stage E; β represents an attenuation factor amplification factor.
5. The method for controlling the data driving voltage of the intelligent energy storage soft switch considering the data quality as claimed in claim 1, wherein the step 4) of suppressing the data according to the voltage deviation disturbance at the time t
Figure FDA0003768626700000041
Establishing a pseudo Jacobian matrix of the r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the moment t
Figure FDA0003768626700000042
Is represented as follows:
Figure FDA0003768626700000043
if it is
Figure FDA0003768626700000044
Or
Figure FDA0003768626700000045
Or
Figure FDA0003768626700000046
Then
Figure FDA0003768626700000047
In the formula,
Figure FDA0003768626700000048
and
Figure FDA0003768626700000049
respectively representing pseudo Jacobian matrixes of an r-th port of an intelligent energy storage soft switch SOP link at the time t and the time t-delta t,
Figure FDA00037686267000000410
is composed of
Figure FDA00037686267000000411
An initial value of (1);
Figure FDA00037686267000000412
Figure FDA00037686267000000413
the variable quantity vector of the active control quantity of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time of t-delta t is represented,
Figure FDA00037686267000000414
Figure FDA00037686267000000415
Figure FDA00037686267000000416
respectively represents the output strategies of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time t-delta t and the time t-2 delta t,
Figure FDA00037686267000000417
respectively representing active power difference vectors and reactive power difference vectors of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time of t-delta t,
Figure FDA00037686267000000418
Figure FDA00037686267000000419
η SOP 、μ SOP representing the weight coefficients.
6. The method for controlling voltage driving according to claim 1, wherein the step 4) is combined with the suppression data of voltage deviation disturbance at time t
Figure FDA00037686267000000420
And calculating a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the t moment by adopting a multi-layer hierarchical prediction algorithm with data of corresponding moments in historical data of node voltage and injection power of the active power distribution network at each moment of m days
Figure FDA00037686267000000421
Is represented as follows:
Figure FDA00037686267000000422
in the formula,
Figure FDA00037686267000000423
and (3) representing a pseudo Jacobian matrix of the energy storage link at the time t, wherein the iterative computation solving process comprises the following steps:
Figure FDA00037686267000000424
if it is
Figure FDA00037686267000000425
Or
Figure FDA00037686267000000426
Or
Figure FDA00037686267000000427
Then
Figure FDA00037686267000000428
In the formula,
Figure FDA00037686267000000429
Figure FDA00037686267000000430
representing the difference vector of the charging and discharging power of the multi-terminal intelligent energy storage soft switch energy storage link at the time of t-delta t,
Figure FDA00037686267000000431
respectively representing the active power output strategies of the energy storage link at the t-delta t moment and the t-2 delta t moment,
Figure FDA00037686267000000432
is composed of
Figure FDA00037686267000000433
Of initial value of ES 、μ ES Represents a weight coefficient, and epsilon represents a threshold coefficient;
in the formula (12), the reaction mixture is,
Figure FDA00037686267000000434
denotes T + n.DELTA.T c The pseudo jacobian prediction matrix at a time, N ∈ {1, \8230;, N }, N = T p /ΔT c Representing the prediction domain T p And control field Δ T c The calculation formula is as follows:
Figure FDA0003768626700000051
in the formula, theta α,t Denotes the prediction coefficient at time T, where α =1, \8230;, l, T d Representing an estimated time duration in days;
defining a regression coefficient vector θ t =(θ 1,t ,…,θ α,t ,…,θ l,t ) T ,θ 1,t 、θ α,t 、θ l,t Respectively representing time t theta t Middle 1, alpha, l element, theta t The iterative solution formula is as follows:
Figure FDA0003768626700000052
in the formula, a pseudo-elegant predictive matrix set of t-delta t
Figure FDA0003768626700000053
Figure FDA0003768626700000054
And
Figure FDA0003768626700000055
respectively, using T + n Δ T c -T d Time sum T + n Δ T c -lT d A pseudo Jacobian matrix is obtained by calculating historical measurement data of the active power distribution network at any moment, delta represents a weight coefficient, and theta t-Δt Is the regression coefficient vector at the time t-delta t.
7. The intelligent energy storage soft switching data driving voltage control method considering data quality according to claim 1, wherein the objective function in step 5) is expressed as follows:
Figure FDA0003768626700000056
wherein J represents an objective function, J 1 、J′ 1 The adjusted objective function is represented as a function of the target,
Figure FDA0003768626700000057
representing the voltage reference at time t + deltat,
Figure FDA0003768626700000058
the average measured voltage data at time t is shown,
Figure FDA0003768626700000059
respectively are the upper limit and the lower limit of a voltage reference interval,
Figure FDA00037686267000000510
Figure FDA00037686267000000511
r represents the port number of the multi-terminal intelligent energy storage soft switch,
Figure FDA00037686267000000512
representing the difference vector of the charging and discharging power of the multi-terminal intelligent energy storage soft switch energy storage link at the time t,
Figure FDA00037686267000000513
respectively represents the charging and discharging power of the energy storage link at the time t and the time t-delta t,
Figure FDA00037686267000000514
the variable quantity vector of the active and reactive control quantity of the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the time t is shown,
Figure FDA00037686267000000515
respectively representing active power difference value vector and reactive power difference value vector of the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the moment t, lambda ES 、λ SOP Represents a weight coefficient, ζ represents an attenuation factor;
Figure FDA00037686267000000516
the estimated voltage value of each node at the time t + Δ t is expressed as follows:
Figure FDA00037686267000000517
Figure FDA0003768626700000061
in the formula,
Figure FDA0003768626700000062
representing the vector of the voltage regulation and control variation quantity of the SOP link of the multi-terminal intelligent energy storage soft switch at the time t,
Figure FDA0003768626700000063
a pseudo Jacobian matrix of an r-th port of a multi-terminal intelligent energy storage soft switch SOP link at the time t is represented;
Figure FDA0003768626700000064
denotes from time T to T + N Δ T c The voltage regulation and control variable vector of the multi-terminal intelligent energy storage soft switch energy storage link at a moment,
Figure FDA0003768626700000065
respectively, time T, T + (N-1) Δ T c Constantly multi-terminal intelligent energy storage soft switch energy storage link voltage regulation and control variable vectors;
Figure FDA0003768626700000066
expressing a pseudo Jacobian estimation matrix of a multi-terminal intelligent energy storage soft switch energy storage link at the time t;
Figure FDA0003768626700000067
represents the time from T to T + (N-1) Δ T c The charge and discharge power variable quantity vector of the energy storage link at any moment,
Figure FDA0003768626700000068
Figure FDA0003768626700000069
respectively, time T, T + (N-1) Δ T c And (3) charging and discharging power difference vectors of the moment energy storage link.
8. The method for controlling driving voltage of intelligent energy storage soft switch data considering data quality as claimed in claim 1, wherein the capacity of each port converter, the loss constraint of each port converter and the active and reactive power output upper and lower limits constraint of each port converter in the multi-terminal intelligent energy storage soft switch SOP link in step 5) are expressed as follows:
Figure FDA00037686267000000610
in the formula,
Figure FDA00037686267000000611
the actual charging and discharging power of the energy storage link is shown,
Figure FDA00037686267000000612
respectively showing the active power of the r port of a multi-terminal intelligent energy storage soft switch SOP link at the moment tAnd the reactive power output,
Figure FDA00037686267000000613
the loss of the current converter at the r-th port of the multi-terminal intelligent energy storage soft switch SOP link at the time t is shown,
Figure FDA00037686267000000614
the loss coefficient of the current converter at the r-th port of the SOP link of the multi-terminal intelligent energy storage soft switch is shown,
Figure FDA00037686267000000615
respectively represents the upper limit and the lower limit of active power of a multi-terminal intelligent energy storage soft switch SOP link,
Figure FDA00037686267000000616
the actual injection active power of the nth port of the multi-terminal intelligent energy storage soft switch SOP link at the time t is shown,
Figure FDA00037686267000000617
respectively representing the upper and lower limits of reactive power of a multi-terminal intelligent energy storage soft switch SOP link at the time t,
Figure FDA00037686267000000618
and the capacity of the r-th port converter of the multi-terminal intelligent energy storage soft switch SOP link is shown.
9. The method for controlling the data driving voltage of the intelligent energy storage soft switch considering the data quality as claimed in claim 1, wherein the charging and discharging efficiency constraint, the charging and discharging power upper and lower limit constraint and the prediction domain T of the multi-terminal intelligent energy storage soft switch energy storage link in the step 5) are implemented by the following steps p The internal state of charge constraint and the charge-discharge number limit constraint are expressed as follows:
Figure FDA00037686267000000619
Figure FDA0003768626700000071
Figure FDA0003768626700000072
Figure FDA0003768626700000073
Figure FDA0003768626700000074
if it is
Figure FDA0003768626700000075
And is
Figure FDA0003768626700000076
Figure FDA0003768626700000077
In the formula,
Figure FDA0003768626700000078
the actual charging and discharging power of the energy storage link is shown,
Figure FDA0003768626700000079
respectively represents the charge and discharge power of the energy storage link at the time t and the time t-delta t,
Figure FDA00037686267000000710
respectively shows the charge and discharge efficiency of the energy storage link,
Figure FDA00037686267000000711
respectively represents the upper limit and the lower limit of the calculated value of the charging and discharging power of the energy storage link,
Figure FDA00037686267000000712
respectively represent T + DeltaT c Time and T + (N-1) Δ T c Charging and discharging power difference S of the moment energy storage link ES The capacity of the energy storage link is shown,
Figure FDA00037686267000000713
the upper and lower limits of the charge state of the energy storage link are shown,
Figure FDA00037686267000000714
respectively representing the initial state and the prediction domain T of the energy storage system p The state of charge after the end of the cycle,
Figure FDA00037686267000000715
the value of the threshold coefficient is represented by,
Figure FDA00037686267000000716
respectively representing the charge states of the multi-terminal intelligent energy storage soft switch energy storage link at T and T-delta T, wherein N = T p /ΔT c Representing the prediction domain T p And control Domain Δ T c The ratio of (a) to (b),
Figure FDA00037686267000000717
and the charging and discharging power difference of the multi-terminal intelligent energy storage soft switch energy storage link at the time t is represented.
CN202210893877.2A 2022-07-27 2022-07-27 Intelligent energy storage soft switch data driving voltage control method considering data quality Active CN115241901B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210893877.2A CN115241901B (en) 2022-07-27 2022-07-27 Intelligent energy storage soft switch data driving voltage control method considering data quality

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210893877.2A CN115241901B (en) 2022-07-27 2022-07-27 Intelligent energy storage soft switch data driving voltage control method considering data quality

Publications (2)

Publication Number Publication Date
CN115241901A true CN115241901A (en) 2022-10-25
CN115241901B CN115241901B (en) 2023-06-20

Family

ID=83676697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210893877.2A Active CN115241901B (en) 2022-07-27 2022-07-27 Intelligent energy storage soft switch data driving voltage control method considering data quality

Country Status (1)

Country Link
CN (1) CN115241901B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249101A (en) * 2016-06-30 2016-12-21 湖南大学 A kind of intelligent distribution network fault identification method
CN110263391A (en) * 2019-05-31 2019-09-20 广东电网有限责任公司 A kind of active power distribution network intelligence Sofe Switch planing method based on chance constraint
CN111443259A (en) * 2020-03-30 2020-07-24 国网山东省电力公司德州供电公司 Active power distribution network fault diagnosis method and system based on local abnormal factor detection
CN111682594A (en) * 2020-06-15 2020-09-18 天津大学 Data-driven model-free adaptive voltage control method for flexible substation of power distribution network
WO2020223883A1 (en) * 2019-05-06 2020-11-12 Abb Schweiz Ag Power converter and control method therefor
CN114142498A (en) * 2021-11-30 2022-03-04 天津大学 Data-driven distributed energy storage self-adaptive prediction control voltage regulation method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249101A (en) * 2016-06-30 2016-12-21 湖南大学 A kind of intelligent distribution network fault identification method
WO2020223883A1 (en) * 2019-05-06 2020-11-12 Abb Schweiz Ag Power converter and control method therefor
CN110263391A (en) * 2019-05-31 2019-09-20 广东电网有限责任公司 A kind of active power distribution network intelligence Sofe Switch planing method based on chance constraint
CN111443259A (en) * 2020-03-30 2020-07-24 国网山东省电力公司德州供电公司 Active power distribution network fault diagnosis method and system based on local abnormal factor detection
CN111682594A (en) * 2020-06-15 2020-09-18 天津大学 Data-driven model-free adaptive voltage control method for flexible substation of power distribution network
CN114142498A (en) * 2021-11-30 2022-03-04 天津大学 Data-driven distributed energy storage self-adaptive prediction control voltage regulation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YANDA HUO等: "Data-Driven Adaptive Operation of Soft Open Points in Active Distribution Networks", IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS, vol. 17, no. 12, pages 8230 - 8242, XP011874023, DOI: 10.1109/TII.2021.3064370 *
刘灏 等: "基于局部离群因子的PMU连续坏数据检测方法", 电力系统自动化, vol. 46, no. 1, pages 25 - 32 *
马世乾 等: "基于机会约束的有源配电网智能软开关电压控制方法", 电力系统及其自动化学报, vol. 31, no. 9, pages 104 - 115 *

Also Published As

Publication number Publication date
CN115241901B (en) 2023-06-20

Similar Documents

Publication Publication Date Title
CN109861202B (en) Dynamic optimization scheduling method and system for flexible interconnected power distribution network
CN107688879B (en) Active power distribution network distributed power supply planning method considering source-load matching degree
CN105846461B (en) Control method and system for large-scale energy storage power station self-adaptive dynamic planning
CN107301472B (en) Distributed photovoltaic planning method based on scene analysis method and voltage regulation strategy
CN114362196B (en) Multi-time-scale active power distribution network voltage control method
CN109361242B (en) Automatic voltage control method for photovoltaic power generation
CN109873447B (en) Multi-time-level active-reactive power regulation and control method for multi-source cooperative active power distribution network
CN109765787B (en) Power distribution network source load rapid tracking method based on intraday-real-time rolling control
CN114362267B (en) Distributed coordination optimization method for AC/DC hybrid power distribution network considering multi-objective optimization
CN114597969B (en) Power distribution network double-layer optimization method considering intelligent soft switch and virtual power plant technology
CN116365506A (en) Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic
CN115000994B (en) Grouping consistency power distribution method for multiple energy storage units
CN112418496B (en) Power distribution station energy storage configuration method based on deep learning
CN106229995B (en) Based on the backup power source shunt reactor parameter optimization method under the Anti-Typhoon operational mode of wind power plant
CN111614110B (en) Receiving-end power grid energy storage optimization configuration method based on improved multi-target particle swarm optimization
CN115588998A (en) Graph reinforcement learning-based power distribution network voltage reactive power optimization method
CN115313403A (en) Real-time voltage regulation and control method based on deep reinforcement learning algorithm
CN111245032B (en) Voltage prediction control method considering loss reduction optimization of wind power plant collector line
CN113890017A (en) Power distribution network voltage self-adaptive control method based on key measurement
CN113890016B (en) Data-driven multi-time scale voltage coordination control method for power distribution network
CN114548597A (en) Optimization method for alternating current-direct current hybrid optical storage and distribution power grid
CN105207220B (en) A kind of tapping voltage regulation and control method based on progressive learning
CN115912372B (en) Voltage control method and system for high-proportion distributed photovoltaic access distribution network
CN115693687B (en) Autonomous regulation and control method for distribution network of transformer area in consideration of comprehensive regulation cost
CN115276067B (en) Distributed energy storage voltage adjusting method suitable for dynamic topological change of power distribution network

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