Power distribution network fault recovery method considering island fusion condition

Abstract

The invention discloses a power distribution network fault recovery method considering an island fusion condition, and the method comprises the steps: 1, building a target function of a power distribution network fault recovery model with the maximization of a weighted load recovery amount and the minimization of line loss as targets; 2, comprehensively considering a generator operation constraint, a radial constraint, a power distribution network safe operation constraint and a load shedding constraint, and establishing a power distribution network operation constraint; 3, establishing a first-aid repair team path and scheduling constraints; 4, establishing scheduling and operation constraints of the mobile power supply vehicle; 5, establishing a topology change constraint based on a first-aid repair event; 6, considering the frequency and voltage safety problems during steady-state frequency calculation and island fusion, and establishing island fusion condition constraints; and 7, solving the power distribution network fault recovery model to obtain a power distribution network recovery scheme. According to the method, the safety of voltage and frequency during island fusion is ensured, the mobile power supply vehicle and the emergency repair team are considered cooperatively, the power recovery method of the power distribution network in the fault state is effectively obtained, the load recovery is accelerated, and the power failure loss is reduced.

Description

Technical Field

[0001] The invention relates to the field of distribution network fault recovery, and in particular to a distribution network fault recovery method taking into account island fusion conditions. Background Art

[0002] In recent years, global climate change has led to frequent natural disasters, making the operating environment of the power system more complex. Under disaster conditions, the risks faced by the distribution network have increased significantly, especially when the local power grid is operating in an isolated island, the stability and recovery ability of the system are severely tested. Specifically, the impact of disasters on the distribution network is manifested in many aspects. First, the frequent failures make it difficult for conventional load recovery strategies to cope with complex and changeable failure scenarios. Second, during the island fusion process, the fluctuations of key parameters such as voltage and frequency of the power grid intensify, which further increases the difficulty of recovery. In addition, the scheduling optimization of the repair team and mobile power vehicles also faces many technical challenges, which affects the efficiency of post-disaster recovery. The current recovery method has not fully studied the load recovery strategy under the condition of island fusion, and lacks a comprehensive analysis of the coordinated scheduling of mobile power vehicles and repair teams. This situation makes it difficult to effectively guarantee the recovery efficiency and reliability of the distribution network after a disaster. Therefore, it is necessary to explore a more systematic recovery plan to cope with the ever-changing power system operating environment. Summary of the invention

[0003] In order to overcome the deficiencies in the above-mentioned prior art, the present invention provides a distribution network fault recovery method taking into account the island fusion condition, so as to consider the impact of voltage and frequency safety during island fusion on the load recovery of the distribution network after a fault, and collaboratively consider the impact of emergency repair teams and mobile power vehicles on the distribution network fault recovery, while taking into account the unsafe factors caused by frequent topology changes, so as to obtain a distribution network fault recovery strategy recovery plan taking into account the island fusion condition, thereby accelerating load recovery and reducing power outage losses.

[0004] In order to achieve the above object, the technical solution adopted by the present invention is:

[0005] The method for fault recovery of a distribution network considering island fusion conditions of the present invention is characterized in that it is performed according to the following steps:

[0006] Step 1: With the goal of maximizing the weighted load recovery amount and minimizing the network loss, the objective function F of the distribution network fault recovery model is established using formula (1):

[0007] (1)

[0008] In formula (1): is the node number set of the distribution network; is a set of discrete recovery time interval sequence numbers; is the line number set of the distribution network; is the weight of the load at the i-th node; is the active power restored by the load at the i-th node at the t-th time step; is the line l between the i-th node and the j-th node ij The magnitude of the square of the current at the last t-th time step; For line l ij The size of the resistance;

[0009] Step 2: Establish distribution network operation constraints of distribution network fault recovery model;

[0010] Step 3: Establish the repair team path and dispatch constraints of the distribution network fault recovery model;

[0011] Step 4: Establish the dispatching and operation constraints of the mobile power vehicle for the distribution network fault recovery model;

[0012] Step 5: Establish topology change constraints based on emergency repair events for the distribution network fault recovery model;

[0013] Step 6: Establish island fusion constraints for distribution network fault recovery model;

[0014] Step 7: Solve the distribution network fault recovery model and obtain the distribution network fault recovery plan considering the island fusion condition:

[0015] Step 7.1, initialization parameters;

[0016] Step 7.1.1, obtain the initial data of the distribution network in the load recovery stage, including: the upper and lower limits of the output of the generator set, the climbing data of the generator set, the active and reactive power requirements of each node, the line resistance reactance data, the maximum load shedding rate of each node, the upper limit of the line transmission power, and the upper and lower limits of the distribution network voltage;

[0017] Step 7.1.2, obtain the initial data of the repair team's route and dispatch, including: the repair time of each fault location, the time for the repair team to move between the fault locations, and the number of repair teams in each warehouse;

[0018] Step 7.1.3, obtain the dispatch and initial operation data of the mobile power vehicle, including: the moving time of the mobile power vehicle between the candidate locations, the number of mobile power vehicles that can be connected to each candidate location at the same time, and the maximum active and reactive power output data of each mobile power vehicle;

[0019] Step 7.1.4, obtain the initial data of the island fusion conditions, including: the steady-state frequency reference value of the distribution network, the upper and lower limits of the distribution network frequency, and the unit regulation power value of each generator set;

[0020] Step 7.2: Use the solver to solve the distribution network fault recovery model and obtain the distribution network fault recovery plan that takes into account the island fusion conditions, including: the recovery time of nodes and lines, the dispatching output value of generator sets, the path and dispatching value of the repair team, and the dispatching of mobile power vehicles and their power generation.

[0021] The method for fault recovery of a distribution network considering island fusion conditions described in the present invention is also characterized in that step 2 is performed according to the following steps:

[0022] Step 2.1: Use equations (2) to (4) to establish generator operation constraints:

[0023] (2)

[0024] (3)

[0025] (4)

[0026] In formula (2) to formula (4): is the node number set where the generator set is located; For the A Boolean variable representing the node recovery state of the generator set at the node at the tth time step. =0, indicating the The generator set at the node has not been restored at the tth time step. =1, indicating the The generator set at the node is restored at the tth time step; For the The minimum active output value of the generator set at each node; For the The maximum active output value of the generator set at each node; For the The active power output value of the generator set at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For the The minimum reactive power output value of the generator set at each node; For the The maximum reactive power output value of the generator set at each node; For the The reactive power output value of the generator set at the node at the tth time step; For the The lower limit of the active ramp power of the generator set at each node; For the The upper limit of the active ramp power of the generator set at each node; For the The active power output value of the generator set at the node at the t+1th moment;

[0027] Step 2.2: Use equations (5) to (9) to establish radial constraints:

[0028] (5)

[0029] (6)

[0030] (7)

[0031] (8)

[0032] (9)

[0033] In formula (5) to formula (9): is the sequence number set of candidate nodes where the mobile power vehicle is located; For line l ij Boolean variable indicating the line recovery state at the tth time step. =0, indicating line l ij If it is not recovered at the tth time step, =1, indicating line l ij Recover at the tth time step; For the The load at the node is a Boolean variable representing the node recovery state at the tth time step. If =0, indicating the The load at the node has not recovered in the tth time step. =1, indicating the The load at the node is restored at the tth time step; The generator set and mobile power supply vehicle are located in the A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the generator set and mobile power supply vehicle are located in the A node is not a root node at the tth time step, if =1, indicating that the generator set and mobile power supply vehicle are located in the The node is the root node at the tth time step; For line l ij Virtual transmission power at time step t; is the line l between the jth node and the ith node ji Virtual transmission power at time step t; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; M is an infinite parameter;

[0034] Step 2.3: Use equations (10) to (13) to establish distribution network safety operation constraints:

[0035] (10)

[0036] (11)

[0037] (12)

[0038] (13)

[0039] In formula (10) to formula (13): For line l ij The maximum transmittable active power; For line l ij The active power transmitted at the tth time step; For line l ij The maximum transferable reactive power; For line l ij The reactive power transmitted at the tth time step; For line l ij The maximum square value of the current that can be transmitted; is the square value of the minimum voltage of the node; is the square value of the maximum voltage of the node; is the square value of the voltage of the ith node at the tth time step;

[0040] Step 2.4: Use equations (14) to (17) to establish load shedding constraints:

[0041] (14)

[0042] (15)

[0043] (16)

[0044] (17)

[0045] In formula (14) to formula (17): is the rated active power demand of the load at the i-th node; is the active load shedding power of the load at the i-th node at the t-th time step; is the rated reactive power demand of the load at the i-th node; is the reactive load shedding power of the load at the i-th node at the t-th time step; is the maximum load shedding rate of the load at the i-th node and is a parameter between 0 and 1;

[0046] Step 2.5: Use equations (18) to (22) to establish Distflow constraints:

[0047] (18)

[0048] (19)

[0049] (20)

[0050] (twenty one)

[0051] (twenty two)

[0052] In formula (18)-formula (22): For the The active power output value of the generator set at the node at the tth time step; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For line l ji The active power transmitted at the tth time step; For line l ji The size of the resistance; For the The reactive power output value of the generator set at the node at the tth time step; For line l ji The magnitude of the square of the current at the last t-th time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step; For line l ji The reactive power transmitted at the tth time step; For line l ij The size of the reactance; For line l ji The size of the reactance; is the square value of the voltage at the jth node at the tth time step.

[0053] Further, the step 3 is performed as follows:

[0054] Step 3.1: Use equations (23) to (26) to establish the path constraints for the repair team:

[0055] (twenty three)

[0056] (twenty four)

[0057] (25)

[0058] (26)

[0059] In formula (23) to formula (26): is a set of fault location serial numbers; Gather the repair team and the warehouse serial number where the mobile power vehicle is located; is a Boolean variable indicating that the repair team moves from the sth warehouse to the kth fault location, =0 means that the repair team did not move from the sth warehouse to the kth fault location, =1 means the repair team moves from the sth warehouse to the kth fault location; is the number of emergency repair teams in the s-th warehouse; is a Boolean variable indicating that the repair team moves from the kth fault location to the sth warehouse, =0 means that the repair team did not move from the kth fault location to the sth warehouse. =1 means the repair team moves from the kth fault location to the sth warehouse; For the repair team from A Boolean variable that moves the kth fault location to the kth fault location, =0 means the repair team did not start from The fault location moves to the kth fault location, =1 means the repair team starts from The fault location moves to the kth fault location; For the repair team to move from the kth fault location to the A Boolean variable for the fault location shift, = 0 means that the repair team did not move from the kth fault location to the The fault location moves, =1 means the repair team moves from the kth fault location to the The fault location moves;

[0060] Step 3.2: Use equations (27) to (32) to establish the dispatch constraints for the repair team:

[0061] (27)

[0062] (28)

[0063] (29)

[0064] (30)

[0065] (31)

[0066] (32)

[0067] In formula (27) to formula (32): is the time when the repair team arrives at the kth fault location; For the repair team to arrive The moment of the fault location; Repair the first The repair time required for each fault location; For the repair team from The moving time required for the kth fault location to move to the kth fault location; For the The node and Faulty line between nodes A Boolean variable indicating whether the repair is performed at the tth time step. =0, indicating a faulty line If it is not repaired at the tth time step, =1, indicating a faulty line Repaired at the tth time step; The repair time required for the repair team to repair the kth fault location; is the set of faulty lines corresponding to the kth fault location; is a positive number close to 0.

[0068] Further, step 4 is performed as follows:

[0069] Step 4.1: Use equations (33) to (36) to establish the mobile power vehicle dispatch constraints:

[0070] (33)

[0071] (34)

[0072] (35)

[0073] (36)

[0074] In formula (33)-formula (36): It is the serial number collection of the mobile power vehicle; Is the mth mobile power vehicle connected to the first A Boolean variable of a node, if =0, indicating that the mth mobile power vehicle is not connected to the tth time step. nodes, if =1, indicating that the mth mobile power vehicle is connected to the nodes; For the The maximum number of mobile power vehicles that a node can access at the same time; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For mobile power car from Node to The moving time of each node; is the total time length; The mobile power car is located in A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the mobile power supply vehicle is located in A node is not a root node at the tth time step, if =1, indicating that the mobile power supply vehicle is located in The node is the root node at the tth time step;

[0075] Step 4.2: Use equations (37) to (40) to establish the operation constraints of the mobile power vehicle:

[0076] (37)

[0077] (38)

[0078] (39)

[0079] (40)

[0080] In formula (37)-formula (40): is the active power output value of the mth mobile power vehicle at the tth time step; is the maximum active power output of the mth mobile power vehicle; is the reactive power output value of the mth mobile power vehicle at the tth time step; is the maximum reactive power output of the mth mobile power vehicle; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step.

[0081] Further, the step 5 is performed as follows:

[0082] Step 5.1: Use equations (41) to (43) to establish the extraction constraints of emergency repair event information:

[0083] (41)

[0084] (42)

[0085] (43)

[0086] In formula (41)-formula (43): For fault line A Boolean variable indicating whether power can be restored at the tth time step. =0, indicating a faulty line If power cannot be restored at the tth time step, =1, indicating a faulty line Power can be restored at the tth time step; For fault line Is it in the A Boolean variable that is fixed in time steps. If =0, indicating a faulty line In the If the time step is not repaired, =1, indicating a faulty line In the time steps are repaired; is a Boolean variable indicating whether the island topology can be changed at the tth time step. =0, indicating that the island topology cannot be changed at the tth time step. =1, indicating that the island topology can be changed at the tth time step; For fault line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating a faulty line If it is not recovered at the t+1th time step, =1, indicating a faulty line Recover at the t+1th time step;

[0087] Step 5.2: Use equations (44) to (47) to establish topology change constraints:

[0088] (44)

[0089] (45)

[0090] (46)

[0091] (47)

[0092] In formula (44)-formula (47): is a Boolean variable indicating the node recovery state of the load at the i-th node at the t+1-th time step. If = 0, indicating that the load at the i-th node has not recovered at the t+1th time step. =1, indicating that the load at the i-th node is restored at the t+1th time step; For line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating line If it is not recovered at the t+1th time step, =1, indicating line Restore at the t+1th time step.

[0093] Further, step 6 is performed as follows:

[0094] Step 6.1: Use equations (48) to (53) to establish line node and number constraints:

[0095] (48)

[0096] (49)

[0097] (50)

[0098] (51)

[0099] (52)

[0100] (53)

[0101] In formula (48)-formula (53): Is the i-th node in the t-th time step whether it belongs to the The Boolean variable of the island formed by the nodes. = 0, indicating that the i-th node does not belong to the t-th node at the t-th time step. If the island formed by the nodes =1, indicating that the i-th node belongs to the The islands formed by nodes; For line l ij At the tth time step, does it belong to the The Boolean variable of the island formed by the nodes. =0, indicating line l ij At the tth time step, it does not belong to If the island formed by the nodes =1, indicating line l ij At the tth time step, The islands formed by nodes; Is the jth node in the tth time step whether it belongs to the The Boolean variable of the island formed by the nodes, if = 0, indicating that the jth node does not belong to the tth time step. If the island formed by the nodes =1, indicating that the jth node belongs to the The islands formed by nodes; For the Whether the node belongs to the The Boolean variable of the island formed by the nodes. =0, indicating the The node does not belong to the If the island formed by the nodes =1, indicating the The node belongs to the The islands formed by nodes;

[0102] Step 6.2: Use equations (54) to (58) to establish steady-state frequency calculation constraints:

[0103] (54)

[0104] (55)

[0105] (56)

[0106] (57)

[0107] (58)

[0108] In formula (54)-formula (58): For the The deviation between the frequency of the island formed by the nodes and the reference frequency at the tth time step; is the unit regulation power of the generator set at the i-th node; For the The steady-state frequency of the island formed by the nodes at the tth time step; is the reference frequency of the distribution network; is the lower frequency limit of the distribution network; is the upper frequency limit of the distribution network; is the node steady-state frequency of the i-th node at the t-th time step;

[0109] Step 6.3: Use equations (59) to (62) to establish frequency and voltage safety constraints during island fusion:

[0110] (59)

[0111] (60)

[0112] (61)

[0113] (62)

[0114] In formula (59)-formula (62): is a Boolean variable indicating whether the i-th node belongs to the island formed by the h-th node at the t-th time step. = 0, indicating that the i-th node does not belong to the island formed by the h-th node at the t-th time step. =1, indicating that the i-th node belongs to the island formed by the h-th node at the t-th time step; For line l ij A Boolean variable that indicates whether the node belongs to the island formed by the hth node at the t+1th time step. =0, indicating line l ij At the t+1th time step, it does not belong to the island formed by the hth node. If =1, indicating line l ij The island formed by the hth node at the t+1th time step; For the A Boolean variable indicating whether a node is the root node at the t+1th time step. =0, indicating the A node is not a root node at the t+1th time step, if =1, indicating the The node is the root node at the t+1th time step; For line l ij A Boolean variable indicating whether the line is an island boundary line at the tth time step. =0, indicating line l ij At the tth time step, it is not an island boundary line. If =1, indicating line l ij At the tth time step, it is the island boundary line; is a Boolean variable indicating the node recovery state of the load of the jth node at the tth time step. =0, indicating that the load of the jth node has not recovered at the tth time step. =1, indicating that the load of the jth node is restored at the tth time step; It is the maximum allowable voltage deviation value of the nodes at both ends of the connecting line when island fusion occurs; It is the maximum allowable frequency deviation value of two fused islands when islands are fused; is the node steady-state frequency of the j-th node at the t-th time step.

[0115] An electronic device of the present invention includes a memory and a processor, wherein the memory is used to store a program that supports the processor to execute the distribution network fault recovery method, and the processor is configured to execute the program stored in the memory.

[0116] The present invention provides a computer-readable storage medium, wherein a computer program is stored on the computer-readable storage medium, and the computer program executes the steps of the distribution network fault recovery method when the computer program is executed by a processor.

[0117] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0118] 1. The present invention establishes a model of island fusion conditions by establishing line node and numbering constraints, steady-state frequency calculation constraints, and frequency and voltage safety constraints during island fusion, thereby ensuring the safety of islands during fusion.

[0119] 2. The present invention establishes a collaborative optimization model of the repair team and the mobile power supply vehicle by considering the impact of the repair team and the mobile power supply vehicle on the recovery of the distribution network fault, which speeds up the recovery process of the distribution network and reduces the power outage losses.

[0120] 3. In order to solve the problem of insecurity of distribution network caused by frequent topology changes, the present invention establishes a topology change model based on emergency repair events to make the distribution network safer. BRIEF DESCRIPTION OF THE DRAWINGS

[0121] Figure 1 It is a framework diagram of the method of the present invention. DETAILED DESCRIPTION

[0122] In this embodiment, a distribution network fault recovery method considering island fusion conditions is used in the load recovery stage of the distribution network, and a commercial solver can be called for efficient solution, including: 1. Establishing the objective function of the distribution network fault recovery model with the goal of maximizing the weighted load recovery amount and minimizing the line loss; 2. Comprehensively considering the generator operation constraints, radial constraints, distribution network safety operation constraints, and load shedding constraints, using the Distflow power flow model to establish the distribution network operation constraints; 3. Taking into account the movement path of the repair team, considering the repair team scheduling strategy, and establishing the repair team path and scheduling constraints; 4. Taking into account the impact of the dispatching and operation of mobile power vehicles on the operation of the power system, establish the dispatching and operation constraints of mobile power vehicles; 5. Taking into account the dangers brought by frequent topology updates, establish topology change model constraints based on emergency repair events; 6. Taking into account the line and node numbering constraints, consider the frequency and voltage safety issues during steady-state frequency calculation and island fusion, and establish island fusion condition constraints; 7. Model on AMPL software and call the solver to solve, and obtain the distribution network recovery plan, including: node and line recovery time, generator group dispatch output value, repair team path and dispatch value, mobile power vehicle dispatch and its power generation. Specifically, if Figure 1 As shown, the method is performed in the following steps:

[0123] Step 1: With the goal of maximizing the weighted load recovery amount and minimizing the network loss, the objective function F of the distribution network fault recovery model is established using formula (1):

[0124] (1)

[0125] Formula (1) is the objective function of the model, where the first term is to maximize the weighted load recovery amount and the second term is to minimize the line loss; in formula (1): is the node number set of the distribution network; is a set of discrete recovery time interval sequence numbers; is the line number set of the distribution network; is the weight of the load at the i-th node; is the active power restored by the load at the i-th node at the t-th time step; is the line l between the i-th node and the j-th node ij The magnitude of the square of the current at the last t-th time step; For line l ij The size of the resistance.

[0126] Step 2: Establish the distribution network operation constraints of the distribution network fault recovery model:

[0127] Step 2.1: Use equations (2) to (4) to establish generator operation constraints:

[0128] (2)

[0129] (3)

[0130] (4)

[0131] Formula (2) is the active power output limit of the generator set, where the dispatched active power output of the generator set plus the active power generated by the primary frequency regulation cannot exceed the upper and lower limits of the active power of the generator set; Formula (3) is the reactive power output limit of the generator set, where the reactive power output of the generator set cannot exceed the upper and lower limits of the reactive power of the generator set; Formula (4) is the ramp limit of the active power of the generator set; In Formulas (2) to (4): is the node number set where the generator set is located; For the A Boolean variable representing the node recovery state of the generator set at the node at the tth time step. =0, indicating the The generator set at the node has not been restored at the tth time step. =1, indicating the The generator set at the node is restored at the tth time step; For the The minimum active output value of the generator set at each node; For the The maximum active output value of the generator set at each node; For the The active power output value of the generator set at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For the The minimum reactive power output value of the generator set at each node; For the The maximum reactive power output value of the generator set at each node; For the The reactive power output value of the generator set at the node at the tth time step; For the The lower limit of the active ramp power of the generator set at each node; For the The upper limit of the active ramp power of the generator set at each node; For the The active power output value of the generator set at the node at the t+1th moment.

[0132] Step 2.2: Use equations (5) to (9) to establish radial constraints:

[0133] (5)

[0134] (6)

[0135] (7)

[0136] (8)

[0137] (9)

[0138] Formula (5) ensures the radial shape of the distribution network and requires that the number of restored lines is equal to the number of restored nodes minus the number of root nodes; Formula (6)-Formula (9) ensure the connectivity of the distribution network and is a virtual power flow method, which is achieved by introducing a virtual network; Formula (6) is a virtual power flow balance constraint, the root node sends virtual power, and the restored non-root node absorbs virtual power; Formula (7) means that only the root node can send virtual power; Formula (8) means that the non-root node cannot send virtual power; Formula (9) means that only the restored line can transmit virtual power; In Formula (5)-Formula (9): is the sequence number set of candidate nodes where the mobile power vehicle is located; For line l ij Boolean variable indicating the line recovery state at the tth time step. =0, indicating line l ij If it is not recovered at the tth time step, =1, indicating line l ij Recover at the tth time step; For the The load at the node is a Boolean variable representing the node recovery state at the tth time step. If =0, indicating the The load at the node has not recovered in the tth time step. =1, indicating the The load at the node is restored at the tth time step; The generator set and mobile power supply vehicle are located in the A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the generator set and mobile power supply vehicle are located in the A node is not a root node at the tth time step, if =1, indicating that the generator set and mobile power supply vehicle are located in the The node is the root node at the tth time step; For line l ij Virtual transmission power at time step t; is the line l between the jth node and the ith node ji Virtual transmission power at time step t; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; M is an infinite parameter.

[0139] Step 2.3: Use equations (10) to (13) to establish distribution network safety operation constraints:

[0140] (10)

[0141] (11)

[0142] (12)

[0143] (13)

[0144] Formula (10) means that only the restored line can transmit active power, and it cannot exceed the maximum allowable active transmission power of the line; Formula (11) means that only the restored line can transmit reactive power, and it cannot exceed the maximum allowable reactive transmission power of the line; Formula (12) means that only the restored line can have current passing through, and it cannot exceed the maximum allowable transmission current of the line; Formula (13) is the constraint of the upper and lower limits of the node voltage; In Formula (10)-Formula (13): For line l ij The maximum transmittable active power; For line l ij The active power transmitted at the tth time step; For line l ij The maximum transferable reactive power; For line l ij The reactive power transmitted at the tth time step; For line l ij The maximum square value of the current that can be transmitted; is the square value of the minimum voltage of the node; is the square value of the maximum voltage of the node; is the square value of the voltage of the ith node at the tth time step.

[0145] Step 2.4: Use equations (14) to (17) to establish load shedding constraints:

[0146] (14)

[0147] (15)

[0148] (16)

[0149] (17)

[0150] Formula (14) is used to calculate the actual active power restored at each node; Formula (15) is used to calculate the actual reactive power restored at each node; Formula (16) is the maximum active load shedding limit of each node, and the active load shedding cannot exceed the maximum active load shedding limit of its node; Formula (17) is the calculation formula for reactive load shedding; In Formula (14)-Formula (17): is the rated active power demand of the load at the i-th node; is the active load shedding power of the load at the i-th node at the t-th time step; is the rated reactive power demand of the load at the i-th node; is the reactive load shedding power of the load at the i-th node at the t-th time step; is the maximum load shedding rate of the load at the i-th node and is a parameter between 0 and 1.

[0151] Step 2.5: Use equations (18) to (22) to establish Distflow constraints:

[0152] (18)

[0153] (19)

[0154] (20)

[0155] (twenty one)

[0156] (twenty two)

[0157] Formula (18) and Formula (19) are the power flow balance formulas of active power and reactive power; Formula (20) and Formula (21) are the calculation formulas for the node voltage drop on the restoration line; Formula (22) is the line current expression based on the second-order cone programming technology; In Formula (18)-Formula (22): For the The active power output value of the generator set at the node at the tth time step; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For line l ji The active power transmitted at the tth time step; For line l ji The size of the resistance; For the The reactive power output value of the generator set at the node at the tth time step; For line l ji The magnitude of the square of the current at the last t-th time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step; For line l ji The reactive power transmitted at the tth time step; For line l ij The size of the reactance; For line l ji The size of the reactance; is the square value of the voltage at the jth node at the tth time step.

[0158] Step 3: Establish the repair team path and dispatch constraints of the distribution network fault recovery model:

[0159] Step 3.1: Use equations (23) to (26) to establish the path constraints for the repair team:

[0160] (twenty three)

[0161] (twenty four)

[0162] (25)

[0163] (26)

[0164] Formula (23) and Formula (24) mean that all repair teams start from the warehouse and return to the warehouse, and the number of repair teams dispatched cannot exceed the maximum number of repair teams in the warehouse; Formula (25) means that each repair team can only visit each fault location once at most; Formula (26) means that the repair team must leave a fault location after arriving at it; In Formula (23)-(26): is a set of fault location serial numbers; Gather the repair team and the warehouse serial number where the mobile power vehicle is located; is a Boolean variable indicating that the repair team moves from the sth warehouse to the kth fault location, =0 means that the repair team did not move from the sth warehouse to the kth fault location, =1 means the repair team moves from the sth warehouse to the kth fault location; is the number of emergency repair teams in the s-th warehouse; is a Boolean variable indicating that the repair team moves from the kth fault location to the sth warehouse, =0 means that the repair team did not move from the kth fault location to the sth warehouse. =1 means the repair team moves from the kth fault location to the sth warehouse; For the repair team from A Boolean variable that moves the kth fault location to the kth fault location, =0 means the repair team did not start from The fault location moves to the kth fault location, =1 means the repair team starts from The fault location moves to the kth fault location; For the repair team to move from the kth fault location to the A Boolean variable for the fault location shift, = 0 means that the repair team did not move from the kth fault location to the The fault location moves, =1 means the repair team moves from the kth fault location to the The fault location moves.

[0165] Step 3.2: Use equations (27) to (32) to establish the dispatch constraints for the repair team:

[0166] (27)

[0167] (28)

[0168] (29)

[0169] (30)

[0170] (31)

[0171] (32)

[0172] Formula (27) and Formula (28) are the time for the repair team to arrive at each fault location, which is equal to the time for the repair team to arrive at the previous fault location plus the repair time of a fault location and the moving time from the previous fault location to this fault location; Formula (29) means that each fault location can only be repaired once at most; Formula (30) and Formula (31) are used to calculate the repair completion time of each fault location, and the repair completion time is equal to the time for the repair team to arrive at the fault location plus the repair time of the fault location; Formula (32) means that when no repair team visits a fault location, its arrival time will be set to 0; In Formula (27)-Formula (32): is the time when the repair team arrives at the kth fault location; For the repair team to arrive The moment of the fault location; Repair the first The repair time required for each fault location; For the repair team from The moving time required for the kth fault location to move to the kth fault location; For the The node and Faulty line between nodes A Boolean variable indicating whether the repair is performed at the tth time step. =0, indicating a faulty line If it is not repaired at the tth time step, =1, indicating a faulty line Repaired at the tth time step; The repair time required for the repair team to repair the kth fault location; is the set of faulty lines corresponding to the kth fault location; is a positive number close to 0.

[0173] Step 4: Establish the dispatching and operation constraints of the mobile power vehicle for the distribution network fault recovery model:

[0174] Step 4.1: Use equations (33) to (36) to establish the mobile power vehicle dispatch constraints:

[0175] (33)

[0176] (34)

[0177] (35)

[0178] (36)

[0179] Formula (33) means that a mobile power vehicle can only access one candidate node at the same time step; Formula (34) means that a candidate location can access multiple mobile power vehicles at the same time step; Formula (35) ensures that the transfer of mobile power vehicles between different nodes meets the required transfer time; Formula (36) is the candidate node It can become the root node only when a mobile power vehicle is connected; in equations (33) to (36): It is the serial number collection of the mobile power vehicle; Is the mth mobile power vehicle connected to the first A Boolean variable of a node, if = 0, indicating that the mth mobile power vehicle is not connected to the tth time step. nodes, if =1, indicating that the mth mobile power vehicle is connected to the nodes; For the The maximum number of mobile power vehicles that a node can access at the same time; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For mobile power car from Node to The moving time of each node; is the total time length; The mobile power car is located in A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the mobile power supply vehicle is located in A node is not a root node at the tth time step, if =1, indicating that the mobile power supply vehicle is located in The node is the root node at the tth time step.

[0180] Step 4.2: Use equations (37) to (40) to establish the operation constraints of the mobile power vehicle:

[0181] (37)

[0182] (38)

[0183] (39)

[0184] (40)

[0185] Formula (37) is the active power output limit of the mobile power vehicle, and the output cannot be greater than its maximum active power output value; Formula (38) is the reactive power output limit of the mobile power vehicle, and the output cannot be greater than its maximum reactive power output value; Formula (39) is used to calculate the active power output value of the mobile power vehicle at each node at each time; Formula (40) is used to calculate the reactive power output value of the mobile power vehicle at each node at each time; In Formula (37)-Formula (40): is the active power output value of the mth mobile power vehicle at the tth time step; is the maximum active power output of the mth mobile power vehicle; is the reactive power output value of the mth mobile power vehicle at the tth time step; is the maximum reactive power output of the mth mobile power vehicle; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step.

[0186] Step 5: Establish topology change constraints based on emergency repair events for the distribution network fault recovery model:

[0187] Step 5.1: Use equations (41) to (43) to establish the extraction constraints of emergency repair event information:

[0188] (41)

[0189] (42)

[0190] (43)

[0191] Formula (41) means that the repaired fault line can be energized in all time steps after the repair; Formula (42) means that if the fault is repaired in the tth time step, the island topology can be adjusted in the next time step; Formula (43) means that only when the fault line is repaired can the line be energized; In Formula (41)-Formula (43): For fault line A Boolean variable indicating whether power can be restored at the tth time step. =0, indicating a faulty line If power cannot be restored at the tth time step, =1, indicating a faulty line Power can be restored at the tth time step; For fault line Is it in the A Boolean variable that is fixed in time steps. If =0, indicating a faulty line In the If the time step is not repaired, =1, indicating a faulty line In the time steps are repaired; is a Boolean variable indicating whether the island topology can be changed at the tth time step. =0, indicating that the island topology cannot be changed at the tth time step. =1, indicating that the island topology can be changed at the tth time step; For fault line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating a faulty line If it is not recovered at the t+1th time step, =1, indicating a faulty line Restore at the t+1th time step.

[0192] Step 5.2: Use equations (44) to (47) to establish topology change constraints:

[0193] (44)

[0194] (45)

[0195] (46)

[0196] (47)

[0197] Equations (44) to (47) show that island topology adjustment can only be performed at the next time step after any faulty line is repaired; Equations (44) and (45) show that if =1, the state of the node can be changed at the next moment. = 0, the state of the node cannot be changed; Equations (46) and (47) are =1, the line status can be changed at the next moment. =0, the state of the line cannot be changed; in equations (44) to (47): is a Boolean variable indicating the node recovery state of the load at the i-th node at the t+1-th time step. If = 0, indicating that the load at the i-th node has not recovered at the t+1th time step. =1, indicating that the load at the i-th node is restored at the t+1th time step; For line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating line If it is not recovered at the t+1th time step, =1, indicating line Restore at the t+1th time step.

[0198] Step 6: Establish island fusion constraints for distribution network fault recovery model:

[0199] Step 6.1: Use equations (48) to (53) to establish line node and number constraints:

[0200] (48)

[0201] (49)

[0202] (50)

[0203] (51)

[0204] (52)

[0205] (53)

[0206] Formula (48) means that each restored node of the distribution network belongs to and only belongs to an island; Formula (49) means that each restored line of the distribution network belongs to and only belongs to an island; Formulas (50) and (51) mean that once the node is If the node is in the island formed by the nodes, then its line must also be in this island; Formula (52) ensures that the root node of the island must be in the island to which it belongs, and the location of the root node is the number of the island; Formula (53) is when the node If it is not the root node, then no node belongs to The island formed by the nodes; in formula (48)-formula (53): Is the i-th node in the t-th time step whether it belongs to the The Boolean variable of the island formed by the nodes. = 0, indicating that the i-th node does not belong to the t-th node at the t-th time step. If the island formed by the nodes =1, indicating that the i-th node belongs to the The islands formed by nodes; For line l ij At the tth time step, does it belong to the The Boolean variable of the island formed by the nodes. =0, indicating line l ij At the tth time step, it does not belong to If the island formed by the nodes =1, indicating line l ij At the tth time step, The islands formed by nodes; Is the jth node in the tth time step whether it belongs to the The Boolean variable of the island formed by the nodes. = 0, indicating that the jth node does not belong to the tth time step. If the island formed by the nodes =1, indicating that the jth node belongs to the The islands formed by nodes; For the Whether the node belongs to the The Boolean variable of the island formed by the nodes. =0, indicating the The node does not belong to the If the island formed by the nodes =1, indicating the The node belongs to the The islands formed by nodes;

[0207] Step 6.2: Use equations (54) to (58) to establish steady-state frequency calculation constraints:

[0208] (54)

[0209] (55)

[0210] (56)

[0211] (57)

[0212] (58)

[0213] Formula (54) is the calculation formula of the steady-state frequency deviation of the isolated island; Formula (55) is used to calculate the frequency of each isolated island; Formula (56) is used to limit the frequency of the isolated island to not exceed the set maximum and minimum values; Formula (57) is used to calculate the frequency of each node, and the frequencies of the nodes in each isolated island are the same; Formula (58) is used to calculate the power generated by the primary frequency regulation of each generator set; In Formula (54)-Formula (58): For the The deviation between the frequency of the island formed by the nodes and the reference frequency at the tth time step; is the unit regulation power of the generator set at the i-th node; For the The steady-state frequency of the island formed by the nodes at the tth time step; is the reference frequency of the distribution network; is the lower frequency limit of the distribution network; is the upper frequency limit of the distribution network; is the node steady-state frequency of the i-th node at the t-th time step.

[0214] Step 6.3: Use equations (59) to (62) to establish frequency and voltage safety constraints during island fusion:

[0215] (59)

[0216] (60)

[0217] (61)

[0218] (62)

[0219] Formula (59) is the restriction that there must be a line directly connecting two islands when island fusion is restricted; Formula (60) is the calculation formula for identifying the boundary line of the island; Formula (61) is used to determine whether the boundary line of the two fused islands meets the voltage fusion condition; Formula (62) is used to determine whether the two fused islands meet the frequency fusion condition. If the voltage and frequency fusion conditions are met at the same time, then the two islands can be fused at the next moment; In Formula (59)-Formula (62): is a Boolean variable indicating whether the i-th node belongs to the island formed by the h-th node at the t-th time step. = 0, indicating that the i-th node does not belong to the island formed by the h-th node at the t-th time step. =1, indicating that the i-th node belongs to the island formed by the h-th node at the t-th time step; For line l ij A Boolean variable that indicates whether the node belongs to the island formed by the hth node at the t+1th time step. =0, indicating line l ij At the t+1th time step, it does not belong to the island formed by the hth node. If =1, indicating line l ij The island formed by the hth node at the t+1th time step; For the A Boolean variable indicating whether a node is the root node at the t+1th time step. =0, indicating the A node is not a root node at the t+1th time step, if =1, indicating the The node is the root node at the t+1th time step; For line l ij A Boolean variable indicating whether the line is an island boundary line at the tth time step. =0, indicating line l ij At the tth time step, it is not an island boundary line. If =1, indicating line l ij At the tth time step, it is the island boundary line; is a Boolean variable indicating the node recovery state of the load of the jth node at the tth time step. =0, indicating that the load of the jth node has not recovered at the tth time step. =1, indicating that the load of the jth node is restored at the tth time step; It is the maximum allowable voltage deviation value of the nodes at both ends of the connecting line when island fusion occurs; It is the maximum allowable frequency deviation value of two fused islands when islands are fused; is the node steady-state frequency of the j-th node at the t-th time step.

[0220] Step 7: Solve the distribution network fault recovery model and obtain the distribution network fault recovery plan considering the island fusion condition:

[0221] Step 7.1, initialization parameters;

[0222] Step 7.1.1, obtain the initial data of the distribution network in the load recovery stage, including: the upper and lower limits of the output of the generator set, the climbing data of the generator set, the active and reactive power requirements of each node, the line resistance reactance data, the maximum load shedding rate of each node, the upper limit of the line transmission power, and the upper and lower limits of the distribution network voltage;

[0223] Step 7.1.2, obtain the initial data of the repair team's route and dispatch, including: the repair time of each fault location, the time for the repair team to move between the fault locations, and the number of repair teams in each warehouse;

[0224] Step 7.1.3, obtain the dispatch and initial operation data of the mobile power vehicle, including: the moving time of the mobile power vehicle between the candidate locations, the number of mobile power vehicles that can be connected to each candidate location at the same time, and the maximum active and reactive power output data of each mobile power vehicle;

[0225] Step 7.1.4, obtain the initial data of the island fusion conditions, including: the steady-state frequency reference value of the distribution network, the upper and lower limits of the distribution network frequency, and the unit regulation power value of each generator set;

[0226] Step 7.2: Use the solver to solve the distribution network fault recovery model and obtain the distribution network fault recovery plan that takes into account the island fusion conditions, including: the recovery time of nodes and lines, the dispatching output value of generator sets, the path and dispatching value of the repair team, and the dispatching of mobile power vehicles and their power generation.

[0227] In this embodiment, an electronic device includes a memory and a processor, wherein the memory is used to store a program that supports the processor to execute the above method, and the processor is configured to execute the program stored in the memory.

[0228] In this embodiment, a computer-readable storage medium stores a computer program on the computer-readable storage medium, and the computer program executes the steps of the above method when executed by a processor.

Claims

1. A distribution network fault recovery method considering island fusion conditions, characterized in that: The steps are as follows: Step 1: With the goal of maximizing the weighted load recovery amount and minimizing the network loss, the objective function F of the distribution network fault recovery model is established using formula (1): (1) In formula (1): is the node number set of the distribution network; is a set of discrete recovery time interval sequence numbers; is the line number set of the distribution network; is the weight of the load at the i-th node; is the active power restored by the load at the i-th node at the t-th time step; is the line l between the i-th node and the j-th node ij The magnitude of the square of the current at the last t-th time step; For line l ij The size of the resistance; Step 2: Establish distribution network operation constraints for distribution network fault recovery model; Step 3: Establish the repair team path and dispatch constraints of the distribution network fault recovery model; Step 4: Establish the dispatching and operation constraints of the mobile power vehicle for the distribution network fault recovery model; Step 5: Establish topology change constraints based on emergency repair events for the distribution network fault recovery model; Step 6: Establish island fusion constraints for distribution network fault recovery model; Step 7: Solve the distribution network fault recovery model and obtain the distribution network fault recovery plan considering the island fusion condition: Step 7.1, initialization parameters; Step 7.1.1, obtain the initial data of the distribution network in the load recovery stage, including: the upper and lower limits of the output of the generator set, the climbing data of the generator set, the active and reactive power requirements of each node, the line resistance reactance data, the maximum load shedding rate of each node, the upper limit of the line transmission power, and the upper and lower limits of the distribution network voltage; Step 7.1.2, obtain the initial data of the repair team's route and dispatch, including: the repair time of each fault location, the time for the repair team to move between the fault locations, and the number of repair teams in each warehouse; Step 7.1.3, obtain the dispatch and initial operation data of the mobile power vehicle, including: the moving time of the mobile power vehicle between the candidate locations, the number of mobile power vehicles that can be connected to each candidate location at the same time, and the maximum active and reactive power output data of each mobile power vehicle; Step 7.1.4, obtain the initial data of the island fusion conditions, including: the steady-state frequency reference value of the distribution network, the upper and lower limits of the distribution network frequency, and the unit regulation power value of each generator set; Step 7.2: Use the solver to solve the distribution network fault recovery model and obtain the distribution network fault recovery plan that takes into account the island fusion conditions, including: the recovery time of nodes and lines, the dispatching output value of generator sets, the path and dispatching value of the repair team, and the dispatching of mobile power vehicles and their power generation.

2. A distribution network fault recovery method considering island fusion conditions according to claim 1, characterized in that: Described step 2 is carried out as follows: Step 2.1: Use equations (2) to (4) to establish generator operation constraints: (2) (3) (4) In formula (2) to formula (4): is the node number set where the generator set is located; For the A Boolean variable representing the node recovery state of the generator set at the node at the tth time step. =0, indicating the The generator set at the node has not been restored at the tth time step. =1, indicating the The generator set at the node is restored at the tth time step; For the The minimum active output value of the generator set at each node; For the The maximum active output value of the generator set at each node; For the The active power output value of the generator set at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For the The minimum reactive power output value of the generator set at each node; For the The maximum reactive power output value of the generator set at each node; For the The reactive power output value of the generator set at the node at the tth time step; For the The lower limit of the active ramp power of the generator set at each node; For the The upper limit of the active ramp power of the generator set at each node; For the The active power output value of the generator set at the node at the t+1th moment; Step 2.2: Use equations (5) to (9) to establish radial constraints: (5) (6) (7) (8) (9) In formula (5) to formula (9): is the sequence number set of candidate nodes where the mobile power vehicle is located; For line l ij Boolean variable indicating the line recovery state at the tth time step. =0, indicating line l ij If it is not recovered at the tth time step, =1, indicating line l ij Recover at the tth time step; For the The load at the node is a Boolean variable representing the node recovery state at the tth time step. If =0, indicating the The load at the node has not recovered in the tth time step. =1, indicating the The load at the node is restored at the tth time step; The generator set and mobile power supply vehicle are located in the A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the generator set and mobile power supply vehicle are located in the A node is not a root node at the tth time step, if =1, indicating that the generator set and mobile power supply vehicle are located in the The node is the root node at the tth time step; For line l ij Virtual transmission power at time step t; is the line l between the jth node and the ith node ji Virtual transmission power at time step t; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; For the The virtual power generation of the generator set and mobile power vehicle at the node at the tth time step; M is an infinite parameter; Step 2.3: Use equations (10) to (13) to establish distribution network safety operation constraints: (10) (11) (12) (13) In formula (10) to formula (13): For line l ij The maximum transmittable active power; For line l ij The active power transmitted at the tth time step; For line l ij The maximum transferable reactive power; For line l ij The reactive power transmitted at the tth time step; For line l ij The maximum square value of the current that can be transmitted; is the square value of the minimum voltage of the node; is the square value of the maximum voltage of the node; is the square value of the voltage of the ith node at the tth time step; Step 2.4: Use equations (14) to (17) to establish load shedding constraints: (14) (15) (16) (17) In formula (14) to formula (17): is the rated active power demand of the load at the i-th node; is the active load shedding power of the load at the i-th node at the t-th time step; is the rated reactive power demand of the load at the i-th node; is the reactive load shedding power of the load at the i-th node at the t-th time step; is the maximum load shedding rate of the load at the i-th node and is a parameter between 0 and 1; Step 2.5: Use equations (18) to (22) to establish Distflow constraints: (18) (19) (20) (21) (22) In formula (18)-formula (22): For the The active power output value of the generator set at the node at the tth time step; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The active power generated by the generator set at the node in the frequency regulation at the tth time step; For line l ji The active power transmitted at the tth time step; For line l ji The size of the resistance; For the The reactive power output value of the generator set at the node at the tth time step; For line l ji The magnitude of the square of the current at the last t-th time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step; For line l ji The reactive power transmitted at the tth time step; For line l ij The size of the reactance; For line l ji The size of the reactance; is the square value of the voltage at the jth node at the tth time step.

3. A distribution network fault recovery method considering island fusion conditions according to claim 2, characterized in that: The step 3 is carried out as follows: Step 3.1: Use equations (23) to (26) to establish the path constraints for the repair team: (23) (24) (25) (26) In formula (23) to formula (26): is a set of fault location serial numbers; Gather the repair team and the warehouse serial number where the mobile power vehicle is located; is a Boolean variable indicating that the repair team moves from the sth warehouse to the kth fault location, =0 means that the repair team did not move from the sth warehouse to the kth fault location, =1 means the repair team moves from the sth warehouse to the kth fault location; is the number of emergency repair teams in the s-th warehouse; is a Boolean variable indicating that the repair team moves from the kth fault location to the sth warehouse, =0 means that the repair team did not move from the kth fault location to the sth warehouse. =1 means the repair team moves from the kth fault location to the sth warehouse; For the repair team from A Boolean variable that moves the kth fault location to the kth fault location, =0 means the repair team did not start from The fault location moves to the kth fault location, =1 means the repair team starts from The fault location moves to the kth fault location; For the repair team to move from the kth fault location to the A Boolean variable for the fault location shift, = 0 means that the repair team did not move from the kth fault location to the The fault location moves, =1 means the repair team moves from the kth fault location to the The fault location moves; Step 3.2: Use equations (27) to (32) to establish the dispatch constraints for the repair team: (27) (28) (29) (30) (31) (32) In formula (27) to formula (32): is the time when the repair team arrives at the kth fault location; For the repair team to arrive The moment of the fault location; Repair the first The repair time required for each fault location; For the repair team from The moving time required for the kth fault location to move to the kth fault location; For the The node and Faulty line between nodes A Boolean variable indicating whether the repair is performed at the tth time step. =0, indicating a faulty line If it is not repaired at the tth time step, =1, indicating a faulty line Repaired at the tth time step; The repair time required for the repair team to repair the kth fault location; is the set of faulty lines corresponding to the kth fault location; is a positive number approaching 0.

4. A distribution network fault recovery method considering island fusion conditions according to claim 3, characterized in that: The step 4 is carried out as follows: Step 4.1: Use equations (33) to (36) to establish the mobile power vehicle dispatch constraints: (33) (34) (35) (36) In formula (33)-formula (36): It is the serial number collection of the mobile power vehicle; Is the mth mobile power vehicle connected to the first A Boolean variable of a node, if = 0, indicating that the mth mobile power vehicle is not connected to the tth time step. nodes, if =1, indicating that the mth mobile power vehicle is connected to the nodes; For the The maximum number of mobile power vehicles that a node can access at the same time; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For the mth mobile power car in the Is the time step connected to the A Boolean variable of a node, if =0, indicating that the mth mobile power supply vehicle is in the The time step is not connected to the nodes, if =1, indicating that the mth mobile power supply vehicle is in the Time step access nodes; For mobile power car from Node to The moving time of each node; is the total time length; The mobile power car is located in A Boolean variable indicating whether a node is a root node at the tth time step. =0, indicating that the mobile power supply vehicle is located in A node is not a root node at the tth time step, if =1, indicating that the mobile power supply vehicle is located in The node is the root node at the tth time step; Step 4.2: Use equations (37) to (40) to establish the operation constraints of the mobile power vehicle: (37) (38) (39) (40) In formula (37)-formula (40): is the active power output value of the mth mobile power vehicle at the tth time step; is the maximum active power output of the mth mobile power vehicle; is the reactive power output value of the mth mobile power vehicle at the tth time step; is the maximum reactive power output of the mth mobile power vehicle; For the The active power output value of the mobile power vehicle at the node at the tth time step; For the The reactive power output value of the mobile power vehicle at the node at the tth time step.

5. A distribution network fault recovery method considering island fusion conditions according to claim 4, characterized in that: The step 5 is carried out as follows: Step 5.1: Use equations (41) to (43) to establish the extraction constraints of emergency repair event information: (41) (42) (43) In formula (41)-formula (43): For fault line A Boolean variable indicating whether power can be restored at the tth time step. =0, indicating a faulty line If power cannot be restored at the tth time step, =1, indicating a faulty line Power can be restored at the tth time step; For fault line Is it in the A Boolean variable that is fixed in time steps. If =0, indicating a faulty line In the If the time step is not repaired, =1, indicating a faulty line In the time steps are repaired; is a Boolean variable indicating whether the island topology can be changed at the tth time step. =0, indicating that the island topology cannot be changed at the tth time step. =1, indicating that the island topology can be changed at the tth time step; For fault line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating a faulty line If it is not recovered at the t+1th time step, =1, indicating a faulty line Recover at the t+1th time step; Step 5.2: Use equations (44) to (47) to establish topology change constraints: (44) (45) (46) (47) In formula (44)-formula (47): is a Boolean variable indicating the node recovery state of the load at the i-th node at the t+1-th time step. If = 0, indicating that the load at the i-th node has not recovered at the t+1th time step. =1, indicating that the load at the i-th node is restored at the t+1th time step; For Line Boolean variable of the line recovery state at the t+1th time step, if =0, indicating line If it is not recovered at the t+1th time step, =1, indicating line Restore at the t+1th time step.

6. A distribution network fault recovery method considering island fusion conditions according to claim 5, characterized in that: The step 6 is carried out as follows: Step 6.1: Use equations (48) to (53) to establish line node and number constraints: (48) (49) (50) (51) (52) (53) In formula (48) to formula (53): Is the i-th node in the t-th time step whether it belongs to the The Boolean variable of the island formed by the nodes. = 0, indicating that the i-th node does not belong to the t-th node at the t-th time step. If the island formed by the nodes =1, indicating that the i-th node belongs to the The islands formed by nodes; For line l ij At the tth time step, does it belong to the The Boolean variable of the island formed by the nodes. =0, indicating line l ij At the tth time step, it does not belong to If the island formed by the nodes =1, indicating line l ij At the tth time step, The islands formed by nodes; Is the jth node in the tth time step whether it belongs to the The Boolean variable of the island formed by the nodes. = 0, indicating that the jth node does not belong to the tth time step. If the island formed by the nodes =1, indicating that the jth node belongs to the The islands formed by nodes; For the Whether the node belongs to the The Boolean variable of the island formed by the nodes. =0, indicating the The node does not belong to the If the island formed by the nodes =1, indicating the The node belongs to the The islands formed by nodes; Step 6.2: Use equations (54) to (58) to establish steady-state frequency calculation constraints: (54) (55) (56) (57) (58) In formula (54)-formula (58): For the The deviation between the frequency of the island formed by the nodes and the reference frequency at the tth time step; is the unit regulation power of the generator set at the i-th node; For the The steady-state frequency of the island formed by the nodes at the tth time step; is the reference frequency of the distribution network; is the lower frequency limit of the distribution network; is the upper frequency limit of the distribution network; is the node steady-state frequency of the i-th node at the t-th time step; Step 6.3: Use equations (59) to (62) to establish frequency and voltage safety constraints during island fusion: (59) (60) (61) (62) In formula (59)-formula (62): is a Boolean variable indicating whether the i-th node belongs to the island formed by the h-th node at the t-th time step. = 0, indicating that the i-th node does not belong to the island formed by the h-th node at the t-th time step. =1, indicating that the i-th node belongs to the island formed by the h-th node at the t-th time step; For line l ij A Boolean variable that indicates whether the node belongs to the island formed by the hth node at the t+1th time step. =0, indicating line l ij At the t+1th time step, it does not belong to the island formed by the hth node. If =1, indicating line l ij The island formed by the hth node at the t+1th time step; For the A Boolean variable indicating whether a node is the root node at the t+1th time step. =0, indicating the A node is not a root node at the t+1th time step, if =1, indicating the The node is the root node at the t+1th time step; For line l ij A Boolean variable indicating whether the line is an island boundary line at the tth time step. =0, indicating line l ij At the tth time step, it is not an island boundary line. If =1, indicating line l ij At the tth time step, it is the island boundary line; is a Boolean variable indicating the node recovery state of the load of the jth node at the tth time step. = 0, indicating that the load of the jth node has not recovered at the tth time step. =1, indicating that the load of the jth node is restored at the tth time step; It is the maximum allowable voltage deviation value of the nodes at both ends of the connecting line when island fusion occurs; It is the maximum allowable frequency deviation value of two fused islands when islands are fused; is the node steady-state frequency of the j-th node at the t-th time step.

7. An electronic device, comprising a memory and a processor, characterized in that: The memory is used to store a program that supports the processor to execute the distribution network fault recovery method described in any one of claims 1 to 6, and the processor is configured to execute the program stored in the memory.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the distribution network fault recovery method according to any one of claims 1 to 6 are executed.

Images