A Fault Restoration Method for Distribution Networks Considering Island Integration Conditions

Abstract

The present invention discloses a method for distribution network fault recovery considering island integration conditions, including: 1. Establishing an objective function of a distribution network fault recovery model with the goal of maximizing the weighted load recovery amount and minimizing line losses; 2. Establishing distribution network operation constraints by comprehensively considering generator operation constraints, radial constraints, distribution network safe operation constraints, and load shedding constraints; 3. Establishing repair team path and scheduling constraints; 4. Establishing dispatching and operation constraints of mobile power supply vehicles; 5. Establishing topology change constraints based on repair events; 6. Considering steady-state frequency calculation and frequency and voltage safety issues during island integration, establishing island integration condition constraints; 7. Solving the distribution network fault recovery model to obtain a distribution network recovery plan. The present invention ensures the safety of voltage and frequency during island integration, collaboratively considers mobile power supply vehicles and repair teams, effectively obtains a power recovery method for the distribution network in a fault state, helps to accelerate load recovery, and reduces power outage losses.

Description

Technical Field

[0001] The present invention relates to the field of distribution network fault recovery, and specifically to a distribution network fault recovery method considering island integration conditions. Background Art

[0002] In recent years, frequent natural disasters caused by global climate change have made 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 operates in island mode, 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, frequent faults make it difficult for conventional load recovery strategies to cope with complex and changeable fault situations. Second, during the island integration process, key parameters such as the voltage and frequency of the power grid fluctuate more severely, which further increases the difficulty of recovery. In addition, there are also many technical challenges in the scheduling optimization of repair teams and mobile power supply vehicles, which affect the post-disaster recovery efficiency. The current recovery methods have not fully studied the load recovery strategy under island integration conditions and lack a comprehensive analysis of the coordinated scheduling of mobile power supply vehicles and repair teams. This situation makes it difficult to effectively guarantee the recovery efficiency and reliability of the distribution network after disasters. Therefore, it is necessary to explore a more systematic recovery plan to cope with the changing operating environment of the power system. Summary of the Invention

[0003] In order to overcome the above deficiencies in the prior art, the present invention provides a distribution network fault recovery method considering island integration conditions, aiming to consider the impact of the safety of voltage and frequency during island integration on the load recovery of the distribution network after a fault, and coordinately consider the impact of repair teams and mobile power supply vehicles on the distribution network fault recovery. At the same time, considering the unsafe factors brought by frequent topology changes, a distribution network fault recovery strategy recovery plan considering island integration conditions is obtained, so as to accelerate load recovery and reduce power outage losses.

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

[0005] The distribution network fault recovery method considering island integration conditions of the present invention is characterized in that it is carried out according to the following steps:

[0006] Step 1: Taking the maximization of the weighted load recovery amount and the minimization of the network loss as the objectives, use Equation (1) to establish the objective function F of the distribution network fault recovery model:

[0007] (1)

[0008] In Equation (1): is the set of node numbers of the distribution network; is the set of discrete recovery time interval numbers; is the set of line numbers 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 is the magnitude of the square of the current at the t-th time step on line l; is line l ij is the resistance of line l;

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

[0010] Step 3. Establish the path and scheduling constraints of the repair team for the distribution network fault recovery model;

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

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

[0013] Step 6. Establish the island integration constraints for the distribution network fault recovery model;

[0014] Step 7. Solve the distribution network fault recovery model to obtain the distribution network fault recovery plan considering the island integration conditions:

[0015] Step 7.1. Initialize the 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 sets, the ramp data of the generator sets, the active and reactive power demands of each node, the line resistance and 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 path and scheduling, including: the repair time at each fault location, the time for the repair team to move between each fault location, and the number of repair teams available in each warehouse;

[0018] Step 7.1.3. Obtain the initial data of the scheduling and operation of the mobile power vehicle, including: the moving time of the mobile power vehicle between each candidate location, the data on the number of mobile power vehicles that can be connected simultaneously at each candidate location, 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 integration 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 restoration model to obtain a distribution network fault restoration plan considering island integration conditions, including: the restoration time of nodes and lines, the dispatching output value of the generator set, the path and dispatching value of the emergency repair team, and the dispatching and power generation of the mobile power vehicle.

[0021] The feature of a distribution network fault restoration method considering island integration conditions according to the present invention also lies in that the step 2 is carried out according to the following steps:

[0022] Step 2.1: Establish the generator operation constraints by using equations (2)-(4):

[0023] (2)

[0024] (3)

[0025] (4)

[0026] In equations (2)-(4): is the set of node numbers where the generator sets are located; is the th Boolean variable of the node restoration state of the generator set at the th node at the tth time step. If =0, it means that the generator set at the th node is not restored at the tth time step. If =1, it means that the generator set at the th node is restored at the tth time step; is the minimum active power output value of the generator set at the th node; is the maximum active power output value of the generator set at the th node; is the active power output value of the generator set at the th node at the tth time step; is the additional active power generated by the primary frequency regulation of the generator set at the th node at the tth time step; is the minimum reactive power output value of the generator set at the th node; is the maximum reactive power output value of the generator set at the th node; is the reactive power output value of the generator set at the th node at the tth time step; is the lower limit of the active power ramp rate of the generator set at the th node; The upper limit of the active power ramp rate of the generator set at a node; is the active power output value of the generator set at the th node at the (t + 1)-th moment;

[0027] Step 2.2. Establish the radial constraint by using Equations (5) - (9):

[0028] (5)

[0029] (6)

[0030] (7)

[0031] (8)

[0032] (9)

[0033] In Equations (5) - (9): is the set of serial numbers of the candidate nodes where the mobile power vehicle is located; is the line l ij at the t-th time step, the Boolean variable of the line restoration state. If = 0, it means that the line l ij is not restored at the t-th time step. If = 1, it means that the line l ij is restored at the t-th time step; is the Boolean variable of the node restoration state of the load at the th node at the t-th time step. If = 0, it means that the load at the th node is not restored at the t-th time step. If = 1, it means that the load at the th node is restored at the t-th time step; is the Boolean variable indicating whether the th node where the generator set and the mobile power vehicle are located is the root node at the t-th time step. If = 0, it means that the th node where the generator set and the mobile power vehicle are located is not the root node at the t-th time step. If = 1, it means that the th node where the generator set and the mobile power vehicle are located is the root node at the t-th time step; is the virtual transmission power of the line l ij at the t-th time step; is the virtual transmission power of the line l ji between the j-th node and the i-th node at the t-th time step; is the The virtual power generation of the generator set and mobile power supply vehicle at the th node at the t-th time step; The virtual power generation of the generator set and mobile power supply vehicle at the

[0034] Step 2.3. Establish the safe operation constraints of the distribution network using equations (10)-(13):

[0035] (10)

[0036] (11)

[0037] (12)

[0038] (13)

[0039] In equations (10)-(13): is the maximum active power that can be transmitted by line l ij ; is the active power transmitted by line l ij at the t-th time step; is the maximum reactive power that can be transmitted by line l ij ; is the reactive power transmitted by line l ij at the t-th time step; is the maximum value of the square of the current that can be transmitted by line l ij ; 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 i-th node at the t-th time step;

[0040] Step 2.4. Establish the load shedding constraints using equations (14)-(17):

[0041] (14)

[0042] (15)

[0043] (16)

[0044] (17)

[0045] In equations (14)-(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 in 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 in 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. Establish the Distflow power flow constraint using Equations (18)-(22):

[0047] (18)

[0048] (19)

[0049] (20)

[0050] (21)

[0051] (22)

[0052] In Equations (18)-(22): is the active power output value of the generator set at the -th node in the t-th time step; is the active power output value of the mobile power supply vehicle at the -th node in the t-th time step; is the additional active power generated by the primary frequency regulation of the generator set at the -th node in the t-th time step; is the active power transmitted by line l ji in the t-th time step; is the resistance value of line l ji ; is the reactive power output value of the generator set at the -th node in the t-th time step; is the magnitude of the square of the current on line l ji at the t-th time step; is the reactive power output value of the mobile power supply vehicle at the -th node in the t-th time step; is the reactive power transmitted by line l ji in the t-th time step; is the reactance value of line l ij ; is the reactance value of line l ji ; is the squared voltage value of the j-th node at the t-th time step.

[0053] Furthermore, step 3 is carried out as follows:

[0054] Step 3.1: Establish the repair team path constraint using equations (23) - (26):

[0055] (23)

[0056] (24)

[0057] (25)

[0058] (26)

[0059] In equations (23) - (26): is the set of fault location numbers; is the set of warehouse numbers where the repair teams and mobile power supply vehicles are located; is a Boolean variable indicating that the repair team moves from the s-th warehouse to the k-th fault location, = 0 means the repair team does not move from the s-th warehouse to the k-th fault location, = 1 means the repair team moves from the s-th warehouse to the k-th fault location; is the number of repair teams in the s-th warehouse; is a Boolean variable indicating that the repair team moves from the k-th fault location to the s-th warehouse, = 0 means the repair team does not move from the k-th fault location to the s-th warehouse, = 1 means the repair team moves from the k-th fault location to the s-th warehouse; is a Boolean variable indicating that the repair team moves from the -th fault location to the k-th fault location, = 0 means the repair team does not move from the -th fault location to the k-th fault location, = 1 means the repair team moves from the -th fault location to the k-th fault location; is a Boolean variable indicating that the repair team moves from the k-th fault location to the -th fault location, = 0 means the repair team does not move from the k-th fault location to the -th fault location, = 1 means the repair team moves from the k-th fault location to the -th fault location;

[0060] Step 3.2. Establish the dispatching constraints for the emergency repair teams by using Equations (27) - (32):

[0061] (27)

[0062] (28)

[0063] (29)

[0064] (30)

[0065] (31)

[0066] (32)

[0067] In Equations (27) - (32): is the time when the emergency repair team arrives at the k-th fault location; is the time when the emergency repair team arrives at the -th fault location; is the repair time required for the emergency repair team to repair the -th fault location; is the moving time required for the emergency repair team to move from the -th fault location to the k-th fault location; is the -th node and the -th node between the fault line at the t-th time step whether the repaired Boolean variable, if = 0, indicating that the fault line is not repaired at the t-th time step, if = 1, indicating that the fault line is repaired at the t-th time step; is the repair time required for the emergency repair team to repair the k-th fault location; is the set of fault lines corresponding to the k-th fault location; is a positive number approaching 0.

[0068] Furthermore, the said Step 4 is carried out according to the following steps:

[0069] Step 4.1. Establish the dispatching constraints for the mobile power supply vehicles by using Equations (33) - (36):

[0070] (33)

[0071] (34)

[0072] (35)

[0073] (36)

[0074] In equations (33) - (36): is the set of serial numbers of mobile power supply vehicles; is a Boolean variable indicating whether the m-th mobile power supply vehicle is connected to the -th node at the t-th time step. If = 0, it means that the m-th mobile power supply vehicle is not connected to the -th node at the t-th time step. If = 1, it means that the m-th mobile power supply vehicle is connected to the -th node at the t-th time step; is the maximum number of mobile power supply vehicles that can be connected to the -th node at the same moment; is a Boolean variable indicating whether the m-th mobile power supply vehicle is connected to the -th node at the -th time step. If = 0, it means that the m-th mobile power supply vehicle is not connected to the -th node at the -th time step. If = 1, it means that the m-th mobile power supply vehicle is connected to the -th node at the -th time step; is a Boolean variable indicating whether the m-th mobile power supply vehicle is connected to the -th node at the -th time step. If = 0, it means that the m-th mobile power supply vehicle is not connected to the -th node at the -th time step. If = 1, it means that the m-th mobile power supply vehicle is connected to the -th node at the -th time step; is the moving time of the mobile power supply vehicle from the -th node to the -th node; is the total time length; is a Boolean variable indicating whether the -th node where the mobile power supply vehicle is located is the root node at the t-th time step. If = 0, it means that the -th node where the mobile power supply vehicle is located is not the root node at the t-th time step. If = 1, it means that the -th node where the mobile power supply vehicle is located is the root node at the t-th time step;

[0075] Step 4.2: Establish the operation constraints of the mobile power supply vehicle by using Equations (37) - (40):

[0076] (37)

[0077] (38)

[0078] (39)

[0079] (40)

[0080] In Equations (37) - (40): is the active power output value of the m-th mobile power supply vehicle at the t-th time step; is the maximum active power output of the m-th mobile power supply vehicle; is the reactive power output value of the m-th mobile power supply vehicle at the t-th time step; is the maximum reactive power output of the m-th mobile power supply vehicle; is for the active power output value of the mobile power supply vehicle at the n-th node at the t-th time step; is for the reactive power output value of the mobile power supply vehicle at the n-th node at the t-th time step.

[0081] Furthermore, Step 5 is carried out as follows:

[0082] Step 5.1: Establish the extraction constraints of the emergency repair event information by using Equations (41) - (43):

[0083] (41)

[0084] (42)

[0085] (43)

[0086] In Equations (41) - (43): is the Boolean variable indicating whether the faulty line can be restored to power at the t-th time step. If = 0, it means that the faulty line cannot be restored to power at the t-th time step. If = 1, it means that the faulty line can be restored to power at the t-th time step; is the Boolean variable indicating whether the faulty line is repaired at the -th time step. If = 0 indicates the faulty line At the th time step, it is not repaired. If = 1, it indicates that the faulty line is repaired at the th time step; is a Boolean variable indicating whether the island topology structure can be changed at the t-th time step. If = 0, it means that the island topology structure cannot be changed at the t-th time step. If = 1, it means that the island topology structure can be changed at the t-th time step; is a Boolean variable representing the line restoration status of the faulty line at the (t + 1)-th time step. If = 0, it indicates that the faulty line is not restored at the (t + 1)-th time step. If = 1, it indicates that the faulty line is restored at the (t + 1)-th time step;

[0087] Step 5.2: Establish topology change constraints using equations (44) - (47):

[0088] (44)

[0089] (45)

[0090] (46)

[0091] (47)

[0092] In equations (44) - (47): is a Boolean variable representing the node restoration status of the load at the i-th node at the (t + 1)-th time step. If = 0, it indicates that the load at the i-th node is not restored at the (t + 1)-th time step. If = 1, it indicates that the load at the i-th node is restored at the (t + 1)-th time step; is a Boolean variable representing the line restoration status of the line at the (t + 1)-th time step. If = 0, it indicates that the line is not restored at the (t + 1)-th time step. If = 1, it indicates that the line is restored at the (t + 1)-th time step.

[0093] Furthermore, the said Step 6 is carried out as follows:

[0094] Step 6.1: Establish line node and number constraints using equations (48) - (53):

[0095] (48)

[0096] (49)

[0097] (50)

[0098] (51)

[0099] (52)

[0100] (53)

[0101] In equations (48) to (53): is a Boolean variable indicating whether the i-th node belongs to the island formed by the -th node at the t-th time step. If = 0, it means that the i-th node does not belong to the island formed by the -th node at the t-th time step. If = 1, it means that the i-th node belongs to the island formed by the -th node at the t-th time step; is a Boolean variable indicating whether line l ij belongs to the island formed by the -th node at the t-th time step. If = 0, it means that line l ij does not belong to the island formed by the -th node at the t-th time step. If = 1, it means that line l ij belongs to the island formed by the -th node at the t-th time step; is a Boolean variable indicating whether the j-th node belongs to the island formed by the -th node at the t-th time step. If = 0, it means that the j-th node does not belong to the island formed by the -th node at the t-th time step. If = 1, it means that the j-th node belongs to the island formed by the -th node at the t-th time step; is a Boolean variable indicating whether the -th node belongs to the island formed by the -th node at the t-th time step. If = 0, it means that the -th node does not belong to the island formed by the -th node at the t-th time step. If = 1, it means that the The node at the t-th time step belongs to the island formed by the

[0102] Step 6.2. Establish the steady-state frequency calculation constraint by using Equations (54) - (58):

[0103] (54)

[0104] (55)

[0105] (56)

[0106] (57)

[0107] (58)

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

[0109] Step 6.3. Establish the frequency and voltage safety constraints during island integration by using Equations (59) - (62):

[0110] (59)

[0111] (60)

[0112] (61)

[0113] (62)

[0114] In Equations (59) - (62): is a Boolean variable indicating whether the i-th node at the t-th time step belongs to the island formed by the h-th node. If = 0, it means that the i-th node at the t-th time step does not belong to the island formed by the h-th node. If = 1 indicates that the i-th node belongs to the island formed by the h-th node at the t-th time step; is for line l ij a Boolean variable indicating whether it belongs to the island formed by the h-th node at the (t + 1)-th time step. If = 0, it means that line l ij does not belong to the island formed by the h-th node at the (t + 1)-th time step. If = 1, it means that line l ij belongs to the island formed by the h-th node at the (t + 1)-th time step; is a Boolean variable indicating whether the -th node is the root node at the (t + 1)-th time step. If = 0, it means that the -th node is not the root node at the (t + 1)-th time step. If = 1, it means that the -th node is the root node at the (t + 1)-th time step; is a Boolean variable indicating whether line l ij is an island boundary line at the t-th time step. If = 0, it means that line l ij is not an island boundary line at the t-th time step. If = 1, it means that line l ij is an island boundary line at the t-th time step; is a Boolean variable indicating the node restoration status of the load of the j-th node at the t-th time step. If = 0, it means that the load of the j-th node is not restored at the t-th time step. If = 1, it means that the load of the j-th node is restored at the t-th time step; is the maximum allowable voltage deviation value for connecting the nodes at both ends of the line during island integration; is the maximum allowable frequency deviation value for two integrating islands during island integration; is the steady-state frequency of the j-th node at the t-th time step.

[0115] An electronic device according to the present invention, comprising a memory and a processor, is characterized in that the memory is used to store a program for supporting the processor to execute the distribution network fault recovery method, and the processor is configured to execute the program stored in the memory.

[0116] A computer-readable storage medium according to the present invention, on which a computer program is stored, is characterized in that the computer program executes the steps of the distribution network fault recovery method when run by a processor.

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

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

[0119] 2. By considering the impact of repair teams and mobile power supply vehicles on the fault recovery of the distribution network, the present invention establishes a collaborative optimization model for repair teams and mobile power supply vehicles, accelerating the recovery process of the distribution network and reducing power outage losses.

[0120] 3. Aiming at the problem of distribution network insecurity caused by frequent topology changes, the present invention establishes a topology change model based on repair events, making 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 OF THE EMBODIMENTS

[0122] In this embodiment, a method for fault recovery of a distribution network considering island integration conditions is used in the load recovery stage of the distribution network, and a commercial solver can be called for efficient solution, including: 1. Establish an 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 consider generator operation constraints, radial constraints, distribution network safe operation constraints, load shedding constraints, and use the Distflow power flow model to establish distribution network operation constraints; 3. Consider the movement path of the repair team and the repair team scheduling strategy to establish the path and scheduling constraints of the repair team; 4. Consider the impact of the scheduling and operation of the mobile power supply vehicle on the operation of the power system to establish the scheduling and operation constraints of the mobile power supply vehicle; 5. Consider the danger caused by frequent topology updates to establish topology change model constraints based on repair events; 6. Consider line and node numbering constraints, and consider the frequency and voltage safety issues during steady-state frequency calculation and island integration to establish island integration condition constraints; 7. Model on the AMPL software and call the solver for solution to obtain the distribution network recovery plan, including: the recovery time of nodes and lines, the scheduling output value of the generator set, the path and scheduling value of the repair team, and the scheduling and power generation amount of the mobile power supply vehicle. Specifically, as Figure 1 shown, the method is carried out in the following steps:

[0123] Step 1: With the goal of maximizing the weighted load recovery amount and minimizing the network loss, use Equation (1) to establish the objective function F of the distribution network fault recovery model:

[0124] (1)

[0125] Equation (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 Equation (1): is the set of node numbers of the distribution network; is the set of discrete restoration time interval numbers; is the set of line numbers of the distribution network; is the weight value 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 t-th time step; is the line l ij the resistance value.

[0126] Step 2. Establish the operation constraints of the distribution network for the distribution network fault restoration model:

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

[0128] (2)

[0129] (3)

[0130] (4)

[0131] Equation (2) is the active power output limit of the generator set, where the scheduled active power output of the generator set plus the additional active power from primary frequency regulation cannot exceed the upper and lower limits of the active power of the generator set; Equation (3) is the reactive power output limit of the generator set, and the reactive power output of the generator set cannot exceed the upper and lower limits of the reactive power of the generator set; Equation (4) is the ramp limit of the active power of the generator set; in equations (2) - (4): is the set of node numbers where the generator sets are located; is the Boolean variable of the node restoration state of the generator set at the -th node at the t-th time step. If = 0, it means that the generator set at the -th node is not restored at the t-th time step. If = 1, it means that the generator set at the -th node is restored at the t-th time step; is the minimum active power output value of the generator set at the -th node; is the maximum active power output value of the generator set at the -th node; is the active power output value of the generator set at the -th node at the t-th time step; The active power increase of the generator set at the at the Minimum reactive power output value of the generator set at the is the Maximum reactive power output value of the generator set at the is the Reactive power output value of the generator set at the is the Lower limit of the active power ramp rate of the generator set at the is the Upper limit of the active power ramp rate of the generator set at the is the Active power output value of the generator set at the (t + 1)-th moment at the

[0132] Step 2.2. Establish the radial constraint using Equations (5) - (9):

[0133] (5)

[0134] (6)

[0135] (7)

[0136] (8)

[0137] (9)

[0138] Equation (5) ensures the radiality of the distribution network, requiring the number of restored lines to be equal to the number of restored nodes minus the number of root nodes; Equations (6) - (9) ensure the connectivity of the distribution network, which is a virtual power flow method implemented by introducing a virtual network; Equation (6) is the virtual power flow balance constraint, where the root node emits virtual power and the restored non-root nodes absorb virtual power; Equation (7) indicates that only the root node can emit virtual power; Equation (8) indicates that non-root nodes cannot emit virtual power; Equation (9) indicates that only the restored lines can transmit virtual power flow; in Equations (5) - (9): is the set of serial numbers of the candidate nodes where the mobile power vehicle is located; is for line l ij Boolean variable of the line restoration status at the = 0, indicating that line l ij is not restored at the = 1, indicating that line l ij is restored at the is the A Boolean variable indicating the restoration status of the node load at the t-th time step. If = 0, it means that the load at the -th node is not restored at the t-th time step. If = 1, it means that the load at the -th node is restored at the t-th time step; is a Boolean variable indicating whether the -th node where the generator set and mobile power vehicle are located is the root node at the t-th time step. If = 0, it means that the -th node where the generator set and mobile power vehicle are located is not the root node at the t-th time step. If = 1, it means that the -th node where the generator set and mobile power vehicle are located is the root node at the t-th time step; is the virtual transmission power of line l ij at the t-th time step; is the virtual transmission power of line l ji between the j-th node and the i-th node at the t-th time step; is the virtual power generation power of the generator set and mobile power vehicle at the -th node at the t-th time step; is the virtual power generation power of the generator set and mobile power vehicle at the -th node at the t-th time step; is the virtual power generation power of the generator set and mobile power vehicle at the -th node at the t-th time step; M is a parameter of infinity.

[0139] Step 2.3: Establish the safe operation constraints of the distribution network using Equations (10)-(13):

[0140] (10)

[0141] (11)

[0142] (12)

[0143] (13)

[0144] Equation (10) means that only restored lines can transmit active power and cannot exceed the maximum allowable active power transmission of the line; Equation (11) means that only restored lines can transmit reactive power and cannot exceed the maximum allowable reactive power transmission of the line; Equation (12) means that only restored lines can have current passing through and cannot exceed the maximum allowable transmission current of the line; Equation (13) is the constraint of the upper and lower limits of the node voltage; in Equations (10)-(13): is the maximum active power that can be transmitted by line l ij ; is the active power transmitted by line l ij at the t-th time step; is the maximum reactive power that can be transmitted by line l ij ; is the reactive power transmitted by line l ij at the t-th time step; is the maximum value of the square of the current that can be transmitted by line l ij ; 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 i-th node at the t-th time step.

[0145] Step 2.4. Establish the load shedding constraint using equations (14) - (17):

[0146] (14)

[0147] (15)

[0148] (16)

[0149] (17)

[0150] Equation (14) is used to calculate the actual active power restored at each node; Equation (15) is used to calculate the actual reactive power restored at each node; Equation (16) is the limit of the maximum active load shedding at each node, and the active load shedding cannot exceed the limit of the maximum active load shedding at its node; Equation (17) is the calculation formula for the reactive load shedding amount; in equations (14) - (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 it is a parameter between 0 and 1.

[0151] Step 2.5. Establish the Distflow power flow constraint using equations (18) - (22):

[0152] (18)

[0153] (19)

[0154] (20)

[0155] (21)

[0156] (22)

[0157] Equations (18) and (19) are the power flow balance equations for active power and reactive power; Equations (20) and (21) are the calculation formulas for the node voltage drop on the restored line; Equation (22) is the line current expression based on the second-order cone programming technique; in Equations (18)-(22): is the active power output value of the generator set at the -th node at the t-th time step; is the active power output value of the mobile power vehicle at the -th node at the t-th time step; is the additional active power generated by the primary frequency regulation of the generator set at the -th node at the t-th time step; is the active power transmitted by line l ji at the t-th time step; is the resistance value of line l ji ; is the reactive power output value of the generator set at the -th node at the t-th time step; is the magnitude of the square of the current on line l ji at the t-th time step; is the reactive power output value of the mobile power vehicle at the -th node at the t-th time step; is the reactive power transmitted by line l ji at the t-th time step; is the reactance value of line l ij ; is the reactance value of line l ji ; is the square value of the voltage at the j-th node at the t-th time step.

[0158] Step 3. Establish the repair team path and scheduling constraints for the distribution network fault restoration model:

[0159] Step 3.1. Use Equations (23)-(26) to establish the repair team path constraints:

[0160] (23)

[0161] (24)

[0162] (25)

[0163] (26)

[0164] Equations (23) and (24) are for all emergency repair teams to start from the warehouse and return to the warehouse, and the number of dispatched emergency repair teams cannot exceed the maximum number of emergency repair teams owned by the warehouse; Equation (25) is that each emergency repair team can visit a fault location at most once; Equation (26) is that the emergency repair team must leave a fault location after arriving at it; in Equations (23)-(26): is the set of fault location sequence numbers; is the set of the sequence numbers of the warehouses where the emergency repair teams and mobile power supply vehicles are located; is a Boolean variable for the emergency repair team to move from the sth warehouse to the kth fault location, = 0 indicates that the emergency repair team does not move from the sth warehouse to the kth fault location, = 1 indicates that the emergency repair team moves from the sth warehouse to the kth fault location; is the number of emergency repair teams in the sth warehouse; is a Boolean variable for the emergency repair team to move from the kth fault location to the sth warehouse, = 0 indicates that the emergency repair team does not move from the kth fault location to the sth warehouse, = 1 indicates that the emergency repair team moves from the kth fault location to the sth warehouse; is for the emergency repair team to move from the th fault location to the kth fault location, = 0 indicates that the emergency repair team does not move from the th fault location to the kth fault location, = 1 indicates that the emergency repair team moves from the th fault location to the kth fault location; is for the emergency repair team to move from the kth fault location to the th fault location, = 0 indicates that the emergency repair team does not move from the kth fault location to the th fault location, = 1 indicates that the emergency repair team moves from the kth fault location to the th fault location.

[0165] Step 3.2. Establish the dispatching constraints of the emergency repair teams by using Equations (27)-(32):

[0166] (27)

[0167] (28)

[0168] (29)

[0169] (30)

[0170] (31)

[0171] (32)

[0172] Equations (27) and (28) are the times for the emergency repair team to reach each fault location, which is equal to the time for the emergency repair team to reach the previous fault location plus the repair time of one fault location and the moving time from the previous fault location to this fault location; Equation (29) is that each fault location can only be repaired once at most; Equations (30) and (31) are to calculate the emergency repair completion time of each fault location, and the emergency repair completion time is equal to the time for the emergency repair team to reach this fault location plus the repair time of this fault location; Equation (32) is that when no emergency repair team visits a certain fault location, then its arrival time will be set to 0; in Equations (27)-(32): is the moment when the emergency repair team reaches the k-th fault location; is the moment when the emergency repair team reaches the -th fault location; is the repair time required for the emergency repair team to repair the -th fault location; is the moving time required for the emergency repair team to move from the -th fault location to the k-th fault location; is the -th node and the -th node between the fault line at the t-th time step whether the repaired Boolean variable, if = 0, indicating that the fault line at the t-th time step is not repaired, if = 1, indicating that the fault line at the t-th time step is repaired; is the repair time required for the emergency repair team to repair the k-th fault location; is the set of fault lines corresponding to the k-th fault location; is a positive number approaching 0.

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

[0174] Step 4.1: Use Equations (33)-(36) to establish the mobile power vehicle scheduling constraints:

[0175] (33)

[0176] (34)

[0177] (35)

[0178] (36)

[0179] Equation (33) means that a mobile power vehicle can only access one candidate node at the same time step; Equation (34) means that multiple mobile power vehicles can access a candidate location at the same time step; Equation (35) ensures that the transfer of the mobile power vehicle between different nodes meets the required transfer time; Equation (36) is for the candidate node It can only become the root node when a mobile power vehicle accesses it; in Equations (33)-(36): is the set of serial numbers of mobile power vehicles; is a Boolean variable indicating whether the m-th mobile power vehicle accesses the -th node at the t-th time step. If = 0, it means that the m-th mobile power vehicle does not access the -th node at the t-th time step. If = 1, it means that the m-th mobile power vehicle accesses the -th node at the t-th time step; is the maximum number of mobile power vehicles that can access the -th node at the same moment; is a Boolean variable indicating whether the m-th mobile power vehicle accesses the -th node at the -th time step. If = 0, it means that the m-th mobile power vehicle does not access the -th node at the -th time step. If = 1, it means that the m-th mobile power vehicle accesses the -th node at the -th time step; is a Boolean variable indicating whether the m-th mobile power vehicle accesses the -th node at the -th time step. If = 0, it means that the m-th mobile power vehicle does not access the -th node at the -th time step. If = 1, it means that the m-th mobile power vehicle accesses the -th node at the -th time step; is the time for the mobile power vehicle to transfer from the The moving time of a node to the th node; is the total time length; is the th node where the mobile power vehicle is located. It is a Boolean variable indicating whether it is the root node at the t-th time step. If = 0, it means that the th node where the mobile power vehicle is located is not the root node at the t-th time step. If = 1, it means that the th node where the mobile power vehicle is located is the root node at the t-th time step.

[0180] Step 4.2: Establish the operation constraints of the mobile power vehicle by using Equations (37) - (40):

[0181] (37)

[0182] (38)

[0183] (39)

[0184] (40)

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

[0186] Step Five: Establish the topology change constraint based on the repair event of the distribution network fault recovery model:

[0187] Step 5.1: Establish the extraction constraint of the repair event information by using Equations (41) - (43):

[0188] (41)

[0189] (42)

[0190] (43)

[0191] Equation (41) means that the repaired faulty line can be energized at all time steps after repair; Equation (42) means that if a fault is repaired at the t-th time step, the island topology can be adjusted at the next time step; Equation (43) means that only when the faulty line is repaired can this line be energized; In Equations (41)-(43): is the faulty line is a Boolean variable indicating whether the faulty line can be restored to power at the t-th time step. If = 0, it means that the faulty line cannot be restored to power at the t-th time step. If = 1, it means that the faulty line can be restored to power at the t-th time step; is the faulty line is a Boolean variable indicating whether the faulty line is repaired at the -th time step. If = 0, it means that the faulty line is not repaired at the -th time step. If = 1, it means that the faulty line is repaired at the -th time step; is a Boolean variable indicating whether the island topology can be changed at the t-th time step. If = 0, it means that the island topology cannot be changed at the t-th time step. If = 1, it means that the island topology can be changed at the t-th time step; is the faulty line is a Boolean variable indicating the line restoration status of the faulty line at the (t + 1)-th time step. If = 0, it means that the faulty line is not restored at the (t + 1)-th time step. If = 1, it means that the faulty line is restored at the (t + 1)-th time step.

[0192] Step 5.2. Establish topology change constraints using Equations (44)-(47):

[0193] (44)

[0194] (45)

[0195] (46)

[0196] (47)

[0197] Equations (44) - (47) mean that the island topology adjustment can only be carried out at the next time step after any failed line is repaired; equations (44) and (45) mean that if = 1, the state of the node can be changed at the next moment, if = 0, the state of the node cannot be changed; equations (46) and (47) mean that if = 1, the state of the line can be changed at the next moment, if = 0, the state of the line cannot be changed; in equations (44) - (47): is the Boolean variable of the node recovery state of the load at the i-th node at the (t + 1)-th time step. If = 0, it means that the load at the i-th node is not recovered at the (t + 1)-th time step. If = 1, it means that the load at the i-th node is recovered at the (t + 1)-th time step; is the Boolean variable of the line recovery state of line at the (t + 1)-th time step. If = 0, it means that line is not recovered at the (t + 1)-th time step. If = 1, it means that line is recovered at the (t + 1)-th time step.

[0198] Step 6. Establish the island integration constraint of the distribution network fault recovery model:

[0199] Step 6.1. Use equations (48) - (53) to establish the line node and number constraint:

[0200] (48)

[0201] (49)

[0202] (50)

[0203] (51)

[0204] (52)

[0205] (53)

[0206] Equation (48) means that each restored node in the distribution network belongs to and only belongs to one island; Equation (49) means that each restored line in the distribution network belongs to and only belongs to one island; Equations (50) and (51) mean that once a node is in the island formed by the th node, its line must also be in this island; Equation (52) ensures that the root node of the island must be within the island to which it belongs, and the position of the root node is the number of the island; Equation (53) means that when the node is not the root node, no node belongs to the island formed by the th node; In Equations (48)-(53): is a Boolean variable indicating whether the i-th node belongs to the island formed by the th node at the t-th time step. If = 0, it means that the i-th node does not belong to the island formed by the th node at the t-th time step. If = 1, it means that the i-th node belongs to the island formed by the th node at the t-th time step; is a Boolean variable indicating whether the line l ij belongs to the island formed by the th node at the t-th time step. If = 0, it means that the line l ij does not belong to the island formed by the th node at the t-th time step. If = 1, it means that the line l ij belongs to the island formed by the th node at the t-th time step; is a Boolean variable indicating whether the j-th node belongs to the island formed by the th node at the t-th time step. If = 0, it means that the j-th node does not belong to the island formed by the th node at the t-th time step. If = 1, it means that the j-th node belongs to the island formed by the th node at the t-th time step; is a Boolean variable indicating whether the th node belongs to the island formed by the th node at the t-th time step. If = 0, it means that the th node does not belong to the island formed by the th node at the t-th time step. If = 1, it means that the th node belongs to the island formed by the th node at the t-th time step;

[0207] Step 6.2. Establish the steady-state frequency calculation constraints using Equations (54) - (58):

[0208] (54)

[0209] (55)

[0210] (56)

[0211] (57)

[0212] (58)

[0213] Equation (54) is the calculation formula for the steady-state frequency deviation of the island; Equation (55) is used to calculate the frequency of each island; Equation (56) is to limit the island frequency not to exceed the set maximum and minimum values; Equation (57) calculates the frequency of each node, and the frequencies of the nodes within each island are the same; Equation (58) is to calculate the additional power generated by the primary frequency regulation of each generator set. In Equations (54) - (58): is the deviation between the frequency of the island formed by the th node and the reference frequency at the t-th time step; is the unit regulation power of the generator set at the i-th node; is the steady-state island frequency of the island formed by the th node at the t-th 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 steady-state node frequency of the i-th node at the t-th time step.

[0214] Step 6.3. Establish the frequency and voltage safety constraints during island integration using Equations (59) - (62):

[0215] (59)

[0216] (60)

[0217] (61)

[0218] (62)

[0219] Equation (59) means that when restricting island integration, there must be a line directly connecting two islands; Equation (60) is the calculation formula for identifying the boundary line of the island; Equation (61) is to judge whether the boundary line of the two integrated islands meets the voltage integration condition; Equation (62) is to judge whether the two integrated islands meet the frequency integration condition. If the voltage and frequency integration conditions are met simultaneously, then the two islands will be able to integrate at the next moment. In Equations (59) - (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. If = 0, it means that the i-th node does not belong to the island formed by the h-th node at the t-th time step. If = 1, it means that the i-th node belongs to the island formed by the h-th node at the t-th time step; is the line l ij is a Boolean variable indicating whether it belongs to the island formed by the h-th node at the (t + 1)-th time step. If = 0, it means that the line l ij does not belong to the island formed by the h-th node at the (t + 1)-th time step. If = 1, it means that the line l ij belongs to the island formed by the h-th node at the (t + 1)-th time step; is the -th node's Boolean variable indicating whether it is the root node at the (t + 1)-th time step. If = 0, it means that the -th node is not the root node at the (t + 1)-th time step. If = 1, it means that the -th node is the root node at the (t + 1)-th time step; is the line l ij is a Boolean variable indicating whether it is the island boundary line at the t-th time step. If = 0, it means that the line l ij is not the island boundary line at the t-th time step. If = 1, it means that the line l ij is the island boundary line at the t-th time step; is the Boolean variable of the node restoration status of the load of the j-th node at the t-th time step. If = 0, it means that the load of the j-th node is not restored at the t-th time step. If = 1, it means that the load of the j-th node is restored at the t-th time step; is the maximum allowable voltage deviation value for connecting the two end nodes of the line during island integration; is the maximum allowable frequency deviation value for the two integrated islands during island integration; is the steady-state frequency of the j-th node at the t-th time step.

[0220] Step 7. Solve the distribution network fault restoration model to obtain a distribution network fault restoration plan considering island integration conditions:

[0221] Step 7.1. Initialize parameters;

[0222] Step 7.1.1. Obtain the initial data of the distribution network during the load restoration stage, including: the upper and lower limits of the output of the generator sets, the ramp data of the generator sets, the active and reactive power demands of each node, the line resistance and 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 path and dispatching, including: the repair time of each fault location, the time for the repair team to move between each fault location, and the number of repair teams in each warehouse;

[0224] Step 7.1.3. Obtain the initial data of the dispatching and operation of the mobile power vehicle, including: the moving time of the mobile power vehicle between each candidate location, the data on the number of mobile power vehicles that can be connected simultaneously at each candidate location, 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 integration 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 restoration model to obtain a distribution network fault restoration plan considering island integration conditions, including: the restoration time of nodes and lines, the dispatching output value of the generator sets, the path and dispatching value of the repair team, and the dispatching and power generation amount of the mobile power vehicle.

[0227] In this embodiment, an electronic device includes a memory and a processor. 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, and when the computer program is run by a processor, it executes the steps of the above method.

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 of the 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 the 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 plan 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