METHOD FOR MAKING SPATIO-TEMPORAL COMBINED OPTIMAL SCHEDULING STRATEGY OF MOBILE ENERGY STORAGE (MES) SYSTEM

Abstract

A method for making a spatio-temporal combined optimal scheduling strategy of a mobile energy storage (MES) system includes: inputting data of a power system, a traffic system, and an MES system; setting a time interval, and initializing a time interval counter; inputting real-time fault, traffic, and MES data; and performing rolling optimization and solving, and delivering regulation decision instructions of the MES system, till a fault is removed. The core of the present disclosure is to propose a spatio-temporal combined optimal model of the MES system to describe spatio-temporal coupling statuses of an energy storage vehicle, a traffic network, and a power distribution network. The present disclosure provides guidance for an optimal scheduling decision of the MES system by properly regulating a traveling path and charging and discharging power of the MES system, thereby supporting high-reliability operation of the power distribution network.

Claims

1. A method for making a spatio-temporal combined optimal scheduling strategy of a mobile energy storage (MES) system, comprising the following steps: step 1: inputting data of a power system, a traffic system, and an MES system, specifically comprising: a line parameter, load data, generator data, and a topological structure of the power system, a road parameter and a topological structure of the traffic system, a quantity of MES vehicles, a capacity of an MES battery, an upper limit of output active power, an upper limit of output reactive power, upper and lower capacity limits of the battery, charging and discharging efficiency of the battery, and an average traveling speed of an MES vehicle; step 2: setting a time interval Δt, and initializing a time interval counter N to 0 and initial time t to 0; step 3: inputting real-time fault, traffic, and MES data, specifically comprising: a position of a power line fault, expected repair time, a road congestion coefficient of the traffic system, and a real-time position and battery capacity of the MES vehicle; step 4: optimizing and solving an optimal regulation model of the MES system, and delivering regulation decision instructions of the MES system, comprising a path instruction and a power instruction; and step 5: if the fault still exists after NΔt, setting N=N+1, and returning to step 3; otherwise, terminating a process of making a scheduling strategy.

2. The method for making a spatio-temporal combined optimal scheduling strategy of an MES system according to claim 1, wherein step 4 comprises: establishing the optimal regulation model of the MES system, wherein an objective function of the model is set to a minimum load reduction in a time set [NΔt, T], and is specifically expressed as: $\begin{matrix} \min f = \underset{t \in T}{.Math.} \underset{i \in E_{n}}{.Math.} P_{i, t}^{L S} Δ t + λ \underset{t \in T}{.Math.} \underset{i j \in E_{b}}{.Math.} r_{i j} I_{ij, t} Δ t & (1) \end{matrix}$ wherein in the formula (1), E.sub.n and E.sub.b represent a node set and a branch set in a power distribution network respectively, PLS i,t represents a load reduction of a node i at a time point t, r.sub.ij and I.sub.ij,t represent resistance of a branch ij and a square of a current of the branch ij at the time point t respectively, Δt represents the time interval, $\underset{t \in T}{.Math.} \underset{i j \in E_{b}}{.Math.} P_{i, t}^{LS} Δ t$ represents a load reduction of the power distribution network, $λ \underset{t \in T}{.Math.} \underset{i j \in E_{b}}{.Math.} r_{i j} I_{i, j, t} Δ t$ is used to ensure accuracy of solving the model after second-order cone relaxation optimization, and λ represents a real number greater than 10; and setting model constraints, comprising a path traffic constraint, an energy storage battery constraint, and a power flow constraint of the power distribution network, wherein the path traffic constraint is that the MES vehicle needs to complete a traveling process according to the delivered path instruction; p.sub.i,s,t is a binary decision variable, indicating whether an MES system s is located on the node i at the time point t; and when p.sub.i,s,t=1, it indicates that the MES system s is located on the node i at the time point t; $\begin{matrix} \underset{i \in W_{n}}{.Math.} p_{i, s, t} \leq 1 \forall s \in Ω_{s}, t \in T & (2) \end{matrix}$ wherein in the formula (2), W.sub.n represents a set of nodes in a traffic network, Ω.sub.s represents a set of MES vehicles, and a traveling track of the MES vehicle also needs to meet a constraint of the traffic network to ensure that the MES vehicle can reach an expected position within specified time, as shown below: $\begin{matrix} L_{i, j, t} = round (\underset{k \in W_{ij}}{.Math.} \frac{d_{k}}{γ_{k, t} .Math. v}) \forall i, j \in W_{n}, t \in T & (3) \end{matrix}$
p.sub.i,s,t+p.sub.j,s,t+τ≤1∀i∈W.sub.n,j∈W.sub.n,τ∈[1, . . . ,min{L.sub.i,j,t,T−t}],s∈Ω.sub.s (4) wherein in the formulas (3) and (4), L.sub.i,j,t represents time required by the MES vehicle to travel from the node i to a node j at the time point t; W.sub.ij represents a set of roads within a least time-consuming path between the node i and the node j, and is obtained based on a floyd algorithm; γ.sub.k,t represents a congestion coefficient of a road k at the time point t, and is used to measure an actual traffic capacity of the road based on traffic congestion, road construction, and other factors; d.sub.k and v represent a length of the road k and the average traveling speed of the MES vehicle respectively; and P.sub.j,s,t+τ represents whether the mobile storage energy system s is located on the node j at the time point t; the energy storage battery constraint is related to the power instruction, the power instruction indicates active discharging power PM i,s,t and reactive discharging power custom-character M i,s,t of the energy storage battery after the MES vehicle reaches a destination and is connected to a power grid, and the energy storage battery constraint comprises capacity constraints represented by formulas (5) and (6) and power constraints represented by formulas (7) to (9), as shown below: $\begin{matrix} E_{s, t} = E_{s, t - 1} - \underset{i \in E_{n}}{.Math.} \frac{P_{i, s, t}^{M}}{η_{d i s}} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (5) \\ σ_{\min} E_{s} \leq E_{s, t} \leq σ_{\max} E_{s} \forall s \in Ω_{s}, t \in T & (6) \\ - M .Math. p_{sta, s, t} .Math. P_{s}^{\max} \leq P_{s t a, s, t}^{M} \leq 0 \forall s \in Ω_{s}, t \in T & (7) \\ 0 \leq \frac{P_{i, s, t}^{M}}{η_{dis}} \leq p_{j, s, t} .Math. P_{s}^{\max} \forall i \in E_{n}, j \in W_{n}, s \in Ω_{s}, t \in T & (8) \\ - p_{j, s, t} Q_{s}^{\max} \leq Q_{i, s, t}^{M} \leq p_{j, s, t} Q_{s}^{\max} \forall i \in E_{n}, j \in W_{n}, s \in Ω_{s}, t \in T & (9) \end{matrix}$ wherein in the formulas (5) to (9), E.sub.s,t represents a battery capacity of the MES system s at the time point t, PM i,s,t and custom-character M i,s,t represent active and reactive charging and discharging power of the MES system s on the node i at the time point t respectively, η.sub.dis represents discharging efficiency, E.sub.s represents the capacity of the MES battery, σ.sub.max and σ.sub.min represent the upper and lower capacity limits of the MES battery respectively, Pmax s and custom-character max s represent upper limits of active and reactive discharging power of the MES system s respectively, M represents a real number greater than 100, and sta represents a position of a central MES station; and the power flow constraint of the power distribution network is that a power distribution system needs to meet power balance constraints represented by formulas (10) to (13) and operation safety constraints represented by formulas (14) and (15) after the MES system is orderly connected to the power distribution network in a plug and play manner, as shown below: $\begin{matrix} \underset{j i \in E_{b}}{.Math.} (P_{j i, t} - r_{i j} I_{j i, t}) + P_{i, t}^{L} - P_{i, t}^{L S} - P_{i, t}^{G} - \underset{s \in Ω_{s}}{.Math.} P_{i, s, t}^{M} = \underset{ik \in E_{b}}{.Math.} P_{i k, t} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (10) \\ \underset{j i \in E_{b}}{.Math.} (Q_{j i, t} - x_{i j} I_{j i, t}) + Q_{i, t}^{L} - Q_{i, t}^{L S} - Q_{i, t}^{G} - \underset{s \in Ω_{s}}{.Math.} Q_{i, s, t}^{M} = \underset{i k \in E_{b}}{.Math.} Q_{i k, t} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (11) \\ U_{i, t} - U_{j, t} + (r_{i j}^{2} + x_{i j}^{2}) I_{i j, t} - 2 (r_{i j} P_{i j, t} + x_{i j} Q_{i j, t}) = 0 \forall i \in E_{n}, j \in E_{n}, ij \in E_{b}, t \in T & (12) \\ {.Math. {[\begin{matrix} 2 P_{i j, t} & 2 Q_{i j, t} & I_{i j, t} & - U_{i, t} \end{matrix}]}^{T} .Math.}_{2} \leq I_{i j, t} + U_{i, t} \forall i \in E_{n}, ij \in E_{b}, t \in T & (13) \end{matrix}$ wherein in the formulas (10) to (13), x.sub.ij represents reactance of the branch ij, P.sub.ij,t and custom-character .sub.ij,t represent active power and reactive power flowing from the node i to the node j at the time point t respectively, U.sub.i,t represents a square of a voltage of the node i at the time point t, PL i,t and L i,t represent active power and reactive power of load on the node i at the time point t respectively, PLS i,t and custom-character LS i,t respectively represent active power and reactive power that are reduced for the load on the node i at the time point t, and PG i,t and G i,t respectively represent active power and reactive power that are injected by a generator on the node i at the time point t; and
U.sub.i.sup.min≤U.sub.i,t≤U.sub.i.sup.max∀i∈E.sub.n,t∈T (14)
I.sub.ij,t≤I.sub.ij.sup.max∀ij∈E.sub.b,t∈T (15) wherein in the formulas (14) and (15), Umax i and Umin i represent upper and lower limits of a square of a voltage of the node i respectively, and Imax ij represents a square of a maximum allowable current of the branch ij.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] FIG. 1 is a schematic diagram of operation of an MES system;

[0031] FIG. 2 is a framework diagram of rolling optimization of a spatio-temporal combined scheduling strategy of an MES system;

[0032] FIG. 3 is a flowchart of rolling optimization of a spatio-temporal combined scheduling strategy of an MES system;

[0033] FIG. 4 is a schematic diagram of a topology of a test system;

[0034] FIG. 5 shows an optimal regulation strategy of an MES system in scenario 1;

[0035] FIG. 6 shows an optimal regulation strategy of an MES system in scenario 2;

[0036] FIG. 7 shows an optimal regulation strategy of an MES system in scenario 3; and

[0037] FIG. 8 shows an optimal regulation strategy of an MES system in scenario 4.

DETAILED DESCRIPTION

[0038] A method for making a spatio-temporal combined optimal scheduling strategy of an MES system in the present disclosure is described in more detail below with reference to embodiments and accompanying drawings.

[0039] A method for making a spatio-temporal combined optimal scheduling strategy of an MES system is provided. An operation mode of the MES system is shown in FIG. 1, and a framework and a process of rolling optimization of the spatio-temporal combined scheduling strategy of the MES system are shown in FIG. 2 and FIG. 3. The method includes the following steps:

[0040] Step 1: Input data of a power system, a traffic system, and the MES system. The input data specifically includes: a line parameter, load data, generator data, and a topological structure of the power system, a road parameter and a topological structure of the traffic system, a quantity of MES vehicles, a capacity of an MES battery, an upper limit of output active power, an upper limit of output reactive power, upper and lower capacity limits of the battery, charging and discharging efficiency of the battery, and an average traveling speed of an MES vehicle.

[0041] Step 2: Set a time interval Δt, and initialize a time interval counter N to 0 and initial time t to 0.

[0042] Step 3: Input real-time fault, traffic, and MES data. The input data specifically includes: a position of a power line fault, expected repair time, a road congestion coefficient of the traffic system, and a real-time position and battery capacity of the MES vehicle.

[0043] Step 4: Optimize and solve an optimal regulation model of the MES system, and deliver regulation decision instructions of the MES system, including a path instruction and a power instruction. A best regulation solution of the MES vehicle is obtained from a perspective of reliability under a constraint on coupling between a power distribution network and a traffic network, where corresponding optimal regulation instructions include a path instruction and a power instruction.

[0044] The optimal regulation model of the MES system is established, where an objective function of the model is set to a minimum load reduction in a time set [NΔt, T], and is specifically expressed as:

[00008] $\begin{matrix} \min f = \underset{t \in T}{.Math.} \underset{i \in E_{n}}{.Math.} P_{i, t}^{L S} Δ t + λ \underset{t \in T}{.Math.} \underset{i j \in E_{b}}{.Math.} r_{i j} I_{ij, t} Δ t & (1) \end{matrix}$

[0045] where in the above formula, E.sub.n and E.sub.b represent a node set and a branch set in the power distribution network respectively, PLS i,t represents a load reduction of a node i at a time point t, r.sub.ij and I.sub.ij,t represent resistance of a branch ij and a square of a current of the branch ij at the time point t respectively, Δt represents the time interval, the first item on the right represents a load reduction of the power distribution network, the second item on the right is used to ensure accuracy of solving the model after second-order cone relaxation optimization, and λ, represents a large positive number, and may be a real number greater than 10.

[0046] Model constraints are set, including a path traffic constraint, an energy storage battery constraint, and a power flow constraint of the power distribution network.

[0047] The path traffic constraint is that the MES vehicle needs to complete a traveling process according to the delivered path instruction. p.sub.i,s,t is a binary decision variable, indicating whether an MES system s is located on the node i at the time point t; and when p.sub.i,s,t=1, it indicates that the MES system s is located on the node i at the time point t, as shown below:

[00009] $\begin{matrix} \underset{i \in W_{n}}{.Math.} p_{i, s, t} \leq 1 \forall s \in Ω_{s}, t \in T & (2) \end{matrix}$

[0048] where in the above formula, W.sub.n represents a set of nodes in the traffic network, Ω.sub.s represents a set of MES vehicles, and a traveling track of the MES vehicle also needs to meet a constraint of the traffic network to ensure that the MES vehicle can reach an expected position within specified time, as shown below:

[00010] $\begin{matrix} L_{i, j, t} = round (\underset{k \in W_{ij}}{.Math.} \frac{d_{k}}{γ_{k, t} .Math. v}) \forall i, j \in W_{n}, t \in T & (3) \\ p_{i, s, t} + p_{j, s, t + τ} \leq 1 \forall i \in W_{n}, j \in W_{n}, τ \in [1, .Math., \min {L_{i, j, t}, T - t}], s \in Ω_{s} & (4) \end{matrix}$

[0049] where in the above formulas, L.sub.i,j,t represents time required by the MES vehicle to travel from the node i to a node j at the time point t; W.sub.ij represents a set of roads within a least time-consuming path between the node i and the node j, and is obtained based on a floyd algorithm; γ.sub.k,t represents a congestion coefficient of a road k at the time point t, and is used to measure an actual traffic capacity of the road based on traffic congestion, road construction, and other factors; and d.sub.k and v represent a length of the road k and the average traveling speed of the MES vehicle respectively. P.sub.j,s,t+τ represents whether the MES system s is located on the node j at the time point t.

[0050] The energy storage battery constraint is related to the power instruction, the power instruction indicates active discharging power PM i,s,t and reactive discharging power custom-character M i,s,t of the energy storage battery after the MES vehicle reaches a destination and is connected to a power grid, and the energy storage battery constraint includes capacity constraints represented by formulas (5) and (6) and power constraints represented by formulas (7) to (9), as shown below:

[00011] $\begin{matrix} E_{s, t} = E_{s, t - 1} - \underset{i \in E_{n}}{.Math.} \frac{P_{i, s, t}^{M}}{η_{d i s}} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (5) \\ σ_{\min} E_{s} \leq E_{s, t} \leq σ_{\max} E_{s} \forall s \in Ω_{s}, t \in T & (6) \\ - M .Math. p_{sta, s, t} .Math. P_{s}^{\max} \leq P_{s t a, s, t}^{M} \leq 0 \forall s \in Ω_{s}, t \in T & (7) \\ 0 \leq \frac{P_{i, s, t}^{M}}{η_{dis}} \leq p_{j, s, t} .Math. P_{s}^{\max} \forall i \in E_{n}, j \in W_{n}, s \in Ω_{s}, t \in T & (8) \\ - p_{j, s, t} Q_{s}^{\max} \leq Q_{i, s, t}^{M} \leq p_{j, s, t} Q_{s}^{\max} \forall i \in E_{n}, j \in W_{n}, s \in Ω_{s}, t \in T & (9) \end{matrix}$

[0051] where in the above formulas, E.sub.s,t represents a battery capacity of the MES system s at the time point t, PM i,s,t and custom-character M i,s,t represent active and reactive charging and discharging power of the MES system s on the node i at the time point t respectively, η.sub.dis represents discharging efficiency, E.sub.s represents the capacity of the MES battery, σ.sub.max and σ.sub.min represent the upper and lower capacity limits of the MES battery respectively, Pmax s and custom-character max s represent upper limits of active and reactive discharging power of the MES system s respectively, M represents a large positive number, and may be a real number greater than 100, and sta represents a position of a central MES station.

[0052] The power flow constraint of the power distribution network is that a power distribution system needs to meet power balance constraints represented by formulas (10) to (13) and operation safety constraints represented by formulas (14) and (15) after the MES system is orderly connected to the power distribution network in a plug and play manner, as shown below:

[00012] $\begin{matrix} \underset{j i \in E_{b}}{.Math.} (P_{j i, t} - r_{i j} I_{j i, t}) + P_{i, t}^{L} - P_{i, t}^{L S} - P_{i, t}^{G} - \underset{s \in Ω_{s}}{.Math.} P_{i, s, t}^{M} = \underset{ik \in E_{b}}{.Math.} P_{i k, t} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (10) \\ \underset{j i \in E_{b}}{.Math.} (Q_{j i, t} - x_{i j} I_{j i, t}) + Q_{i, t}^{L} - Q_{i, t}^{L S} - Q_{i, t}^{G} - \underset{s \in Ω_{s}}{.Math.} Q_{i, s, t}^{M} = \underset{i k \in E_{b}}{.Math.} Q_{i k, t} \forall i \in E_{n}, s \in Ω_{s}, t \in T & (11) \\ U_{i, t} - U_{j, t} + (r_{i j}^{2} + x_{i j}^{2}) I_{i j, t} - 2 (r_{i j} P_{i j, t} + x_{i j} Q_{i j, t}) = 0 \forall i \in E_{n}, j \in E_{n}, ij \in E_{b}, t \in T & (12) \\ {.Math. {[\begin{matrix} 2 P_{i j, t} & 2 Q_{i j, t} & I_{i j, t} & - U_{i, t} \end{matrix}]}^{T} .Math.}_{2} \leq I_{i j, t} + U_{i, t} \forall i \in E_{n}, ij \in E_{b}, t \in T & (13) \end{matrix}$

[0053] where in the above formulas, x.sub.i, represents reactance of the branch if, P.sub.ij,t and custom-character .sub.ij,t represent active power and reactive power flowing from the node i to the node j at the time point t respectively, U.sub.i,t represents a square of a voltage of the node i at the time point t, PL i,t and L i,t represent active power and reactive power of load on the node i at the time point t respectively, PLS i,t and custom-character LS i,t respectively represent active power and reactive power that are reduced for the load on the node i at the time point t, and PG i,t and G i,t respectively represent active power and reactive power that are injected by a generator on the node i at the time point t; and

U.sub.i.sup.min≤U.sub.i,t≤U.sub.i.sup.max∀i∈E.sub.n,t∈T (14)

I.sub.ij,t≤I.sub.ij.sup.max∀ij∈E.sub.b,t∈T (15)

[0054] where in the above formulas, Umax i and Umin i represent upper and lower limits of a square of a voltage of the node i respectively, and Imax ij represents a square of a maximum allowable current of the branch ij.

[0055] Step 5: If the fault still exists after NΔt, set N=N+1, and return to step 3; otherwise, terminate a process of making a scheduling strategy.

[0056] The method for making a spatio-temporal combined optimal scheduling strategy of an MES system in the present disclosure realizes a rolling optimal scheduling decision of each MES vehicle when a fault occurs, fully excavates a support potentiality of flexibility of the MES system for system reliability, and provides guidance for making a scheduling strategy of the MES system.

[0057] In the embodiment of the present disclosure, some trunk roads in an urban area of a city are selected as the traffic network including 29 nodes and 49 roads, and a standard power distribution network composed of IEEE33 nodes is used as a test power distribution system. The power distribution network includes 33 nodes, one generator, and 37 branches. A topological relationship between the traffic network and the power distribution network is shown in FIG. 4. For example, three MES vehicles (MES 1, MES 2, and MES 3) are disposed, and each of them adopts a container-type lithium iron phosphate battery. For related parameters, refer to Table 1. The central MES station is located on node 1. Scheduling solutions of the MES system in four scenarios are analyzed and compared. FIG. 5 to FIG. 8 show results of optimal regulation strategies of the MES system in these scenarios.

[0058] In scenario 1, power cuts of lines (30, 31) and (17, 18) are scheduled 3 h in advance, and estimated fault repair time is 5 h.

[0059] In scenario 2, lines (30, 31) and (17, 18) are in a power-cut state, and estimated fault repair time is 5 h.

[0060] In scenario 3, at a first time point in scenario 2, driving on a road (4, 7) is prohibited, and a congestion coefficient of a road (7, 8) is 30%.

[0061] In scenario 4, a fault is the same as the fault in scenario 2, and the time interval is 15 min.

[0062] Impact of a line fault type is analyzed as follows:

[0063] In scenario 1, destinations of the MES vehicles are nodes 18, 31, and 32 respectively. A shortest moving path is 1.fwdarw.19.fwdarw.32.fwdarw.18 for MES 1, 1.fwdarw.19.fwdarw.20.fwdarw.31 for MES 2, and 1.fwdarw.19.fwdarw.32.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.32 for MES 3. The three MES vehicles are charged in advance and connected to the power grid in a full-capacity manner when the fault occurs. MES 3 returns to the central station at a fifth time point to use a new full-capacity battery. According to the above regulation strategy, the system has no load reduction. It can be learned from FIG. 5 that under this strategy, the MES system directly supplies power to the nodes (18, 31, 32) with heaviest load, to reduce a line loss. For related line power in a whole process, refer to Table 2.

[0064] In scenario 2, a moving path is 1.fwdarw.19.fwdarw.20.fwdarw.31 for MES 1, 1.fwdarw.19.fwdarw.20.fwdarw.31.fwdarw.19.fwdarw.20.fwdarw.1.fwdarw.19.fwdarw.32 for MES 2, and 1.fwdarw.19.fwdarw.32.fwdarw.18.fwdarw.33.fwdarw.18.fwdarw.32.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.32.fwdarw.18 for MES 3, and load of the system is reduced by 0.255 MWh. After the fault occurs, nodes 18, 31, 32, and 33 form an isolated network, and a load of 0.51 MW cannot be supplied. Because maximum power supply power of a single MES vehicle is 0.3 MW, at least two MES vehicles need to be connected to the power grid at the same time to meet the electrical load. As a result, MES 2 and MES 3 go to the destination immediately, with a primary goal of minimizing a load reduction, and then return to the central station after 1.5 h and 2.5 h respectively to replace their own batteries. Referring to the regulation strategy in scenario 1, when the battery capacity is sufficient, an optimal access point of the MES vehicle is still a main load point due to a line transmission loss. Table 3 shows related line power in scenario 2.

[0065] It can be learned by comparing FIG. 5 and FIG. 6 that in case of a pre-scheduled power cut, regulation space of the MES system is more abundant, and more emergency support can be provided for the system. In addition, the battery capacity and fault duration are key factors affecting the path instruction of the MES system. If an MES vehicle returns to the central station to replace the battery, load of a corresponding node is supplied by an MES vehicle of another node, which may cause a line transmission loss. Therefore, a fewer quantity of times that the MES vehicle replaces the battery leads to a less power transmission loss. It can be learned that the proposed regulation strategy of the MES system can ensure a minimum load reduction and a minimum transmission loss at the same time.

[0066] Impact of a traffic network fault is analyzed as follows:

[0067] In scenario 3, due to congestion in the traffic network, a moving path is 31.fwdarw.22.fwdarw.33.fwdarw.18 for MES 1, 31.fwdarw.20.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.20.fwdarw.31 for MES 2, and 33.fwdarw.22.fwdarw.31.fwdarw.20.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.20.fwdarw.31.fwdarw.32 for MES 3. To ensure that there is no load reduction in the system, both MES 1 and MES 2 are connected to node 31, because a path from the central station to node 31 is least time-consuming. Then, to reduce a transmission loss, the MES vehicles are connected to nodes 18, 31, and 32 respectively. FIG. 6 and FIG. 7 show that the proposed rolling optimization model can fully consider uncertain factors such as the traffic network fault and modify the regulation strategy in real time, thereby realizing efficient access of the MES system.

[0068] Impact of scheduling time is analyzed as follows:

[0069] In scenario 4, a moving path is 1.fwdarw.19.fwdarw.32.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.32 for MES 1, 1.fwdarw.19.fwdarw.20.fwdarw.31 for MES 2, and 1.fwdarw.19.fwdarw.32.fwdarw.19.fwdarw.1.fwdarw.19.fwdarw.32.fwdarw.18 for MES 3. The regulation time is reduced to more accurately predict operation time of the path, and a path from the central station to node 32 is the shortest, so that the MES vehicles are first connected to node 32, thereby ensuring a minimum load reduction of the system. Correspondingly, to reduce a network loss, the MES vehicles are finally connected to nodes 18, 31, and 32 respectively, which is consistent with the strategies in scenarios 1 to 3. To sum up, proper optimal regulation time needs to be set with reference to an actual site situation, while accuracy of the regulation instruction is improved.

[0070] It should be emphasized that the embodiment in the present disclosure is illustrative rather than restrictive. Therefore, the present disclosure includes but is not limited to the embodiment in the detailed description. All other implementations obtained by those skilled in the art according to the technical solutions of the present disclosure also fall within the protection scope of the present disclosure.

TABLE-US-00001 TABLE 1 Parameter setting Parameter name Value Capacity E.sub.s of an energy storage battery (MWh) 1 Upper capacity limit σ.sub.max of the battery (%) 90 Lower capacity limit σ.sub.min of the battery (%) 10 Discharging efficiency η.sub.dis of the battery (%) 95 Maximum active discharging power Pmax s (MW) 0.3 Maximum reactive discharging power Qmax s (Mvar) 0.3 Time interval Δt (h) 0.5 Average traveling speed v of an MES vehicle (km/h) 20

TABLE-US-00002 TABLE 2 Line power (MW) in scenario 1 Time (h) Line 0 1 2 3 4 5 31-32 0 −0.031 −0.031 0.145 0.135 −0.031 32-33 0 0.036 0.036 −0.075 −0.075 0.036 18-33 0 0.024 0.024 0.135 0.135 0.024 (interconnection switch)

TABLE-US-00003 TABLE 3 Line power (MW) in scenario 2 Time (h) Line 0 1 2 3 4 5 31-32 0 0.122 0.099 0.075 −0.007 −0.007 32-33 0 −0.088 −0.111 0.15 0.02 0.020 18-33 0 −0.09 −0.09 −0.09 0.04 0.040 (interconnection switch)

TABLE-US-00004 TABLE 4 Node parameters of an IEEE 33 node system Upper Lower Active Reactive Reference voltage voltage Node Node load load voltage limit limit No. type (MW) (MVar) (kV) (p.u.) (p.u.) 1 3 0 0 12.66 1.05 0.95 2 1 0.1 0.06 12.66 1.05 0.95 3 1 0.09 0.04 12.66 1.05 0.95 4 1 0.12 0.08 12.66 1.05 0.95 5 1 0.06 0.03 12.66 1.05 0.95 6 1 0.06 0.02 12.66 1.05 0.95 7 1 0.2 0.1 12.66 1.05 0.95 8 1 0.2 0.1 12.66 1.05 0.95 9 1 0.06 0.02 12.66 1.05 0.95 10 1 0.06 0.02 12.66 1.05 0.95 11 1 0.045 0.03 12.66 1.05 0.95 12 1 0.06 0.035 12.66 1.05 0.95 13 1 0.06 0.035 12.66 1.05 0.95 14 1 0.12 0.08 12.66 1.05 0.95 15 1 0.06 0.01 12.66 1.05 0.95 16 1 0.06 0.02 12.66 1.05 0.95 17 1 0.06 0.02 12.66 1.05 0.95 18 1 0.09 0.04 12.66 1.05 0.95 19 1 0.09 0.04 12.66 1.05 0.95 20 1 0.09 0.04 12.66 1.05 0.95 21 1 0.09 0.04 12.66 1.05 0.95 22 1 0.09 0.04 12.66 1.05 0.95 23 1 0.09 0.05 12.66 1.05 0.95 24 1 0.42 0.2 12.66 1.05 0.95 25 1 0.42 0.2 12.66 1.05 0.95 26 1 0.06 0.025 12.66 1.05 0.95 27 1 0.06 0.025 12.66 1.05 0.95 28 1 0.06 0.02 12.66 1.05 0.95 29 1 0.12 0.07 12.66 1.05 0.95 30 1 0.2 0.6 12.66 1.05 0.95 31 1 0.15 0.07 12.66 1.05 0.95 32 1 0.21 0.1 12.66 1.05 0.95 33 1 0.06 0.04 12.66 1.05 0.95

[0071] Values 1, 2, and 3 of the node type represent a PQ node, a PV node, and a slack bus respectively.

TABLE-US-00005 TABLE 5 Generator parameters of an IEEE 33 node system Upper limit Lower limit Upper limit Lower limit of an of the of a of the active active reactive reactive power power power power Generator output output output output position (MW) (MW) (MVar) (MVar) 1 10 0 10 −10

TABLE-US-00006 TABLE 6 Power line parameters of an IEEE 33 node system Line Start End Resistance Reactance No. node node (p.u.) (p.u.) 1 1 2 0.005753 0.002932 2 2 3 0.03076 0.015667 3 3 4 0.022836 0.01163 4 4 5 0.023778 0.01211 5 5 6 0.051099 0.044112 6 6 7 0.01168 0.038608 7 7 8 0.044386 0.014668 8 8 9 0.064264 0.04617 9 9 10 0.065138 0.04617 10 10 11 0.012266 0.004056 11 11 12 0.02336 0.007724 12 12 13 0.091592 0.072063 13 13 14 0.033792 0.04448 14 14 15 0.036874 0.032818 15 15 16 0.046564 0.034004 16 16 17 0.080424 0.107378 17 17 18 0.02321 0.035813 18 2 19 0.010232 0.009764 19 19 20 0.093851 0.084567 20 20 21 0.02555 0.029849 21 21 22 0.04423 0.058481 22 3 23 0.028152 0.019236 23 23 24 0.056028 0.044243 24 24 25 0.055904 0.043743 25 6 26 0.012666 0.006451 26 26 27 0.017732 0.009028 27 27 28 0.066074 0.058256 28 28 29 0.050176 0.043712 29 29 30 0.031664 0.016128 30 30 31 0.060795 0.060084 31 31 32 0.019373 0.02258 32 32 33 0.021276 0.033455 33 8 21 0.124785 0.124785 34 9 15 0.124785 0.124785 35 12 22 0.124785 0.124785 36 18 33 0.031196 0.031196 37 25 29 0.031196 0.031196

TABLE-US-00007 TABLE 7 Corresponding table of node numbers of a power distribution network and a traffic network Node Node in Node in in the Node the power Node in the power Node in power in the distribution the traffic distribution the traffic distribution traffic network network network network network network 1 1 12 26 23 6 2 2 13 27 24 10 3 3 14 28 25 13 4 11 15 17 26 — 5 — 16 18 27 — 6 12 17 19 28 — 7 — 18 21 29 14 8 15 19 4 30 24 9 16 20 5 31 8 10 29 21 9 32 7 11 25 22 23 33 22 Note: “—” indicates that the node in the power distribution network is not coupled to the node in the traffic network.

TABLE-US-00008 TABLE 8 Road parameters of a traffic system Start End Road Start End Road node node length node node length 1 2 2 14 20 2.2 1 4 1.7 14 24 3 2 3 2.2 15 16 5.8 2 5 1.8 16 17 9.2 3 6 1.8 16 29 3.4 3 11 3.3 17 18 2.6 4 5 2 17 28 5.3 4 7 2.3 18 19 2.4 5 6 2.1 18 27 4.5 5 8 1.7 19 21 2.6 6 9 2.7 20 25 2.8 6 10 2.6 20 29 4.4 7 8 2.8 21 22 2.7 7 21 3 22 23 2.3 8 9 2.8 22 27 2.6 8 23 2.6 23 24 3.6 9 10 2.2 23 26 2.2 9 24 2.1 24 25 3.7 10 11 1.9 25 26 2.4 10 13 4 25 29 1.4 11 12 4 26 27 2.2 12 13 1.5 26 28 1.8 12 15 6.6 27 28 1.8 13 14 3.6 28 29 2.4 14 15 2.6

METHOD FOR MAKING SPATIO-TEMPORAL COMBINED OPTIMAL SCHEDULING STRATEGY OF MOBILE ENERGY STORAGE (MES) SYSTEM

Inventors

Cpc classification

Classification Explorer

H02J3/0073

ELECTRICITY

Classification Explorer

G06Q50/06

PHYSICS

Classification Explorer

H02J3/08

ELECTRICITY

Classification Explorer

H02J3/001

ELECTRICITY

Classification Explorer

H02J2203/10

ELECTRICITY

Classification Explorer

H02J3/003

ELECTRICITY

International classification

Classification Explorer

H02J3/00

ELECTRICITY

Classification Explorer

H02J3/08

ELECTRICITY

Abstract

Claims

Description