Optimization Method and System for Post-Earthquake Repair Sequence of Substations

Abstract

This invention discloses a method and system for optimizing the post-earthquake repair sequence of a substation. The method includes: acquiring data information of the target substation; constructing a topology model of the target substation and determining the repair time and post-earthquake recovery process of each electrical device in the target substation; acquiring data information of damaged equipment in the target substation after the main shock of the earthquake; performing post-earthquake repair on the equipment in the target substation and updating the equipment status in real time during aftershocks; obtaining the functional recovery curve of the target substation and calculating the toughness loss index of the target substation; and solving for the post-earthquake repair sequence of the target substation based on a genetic algorithm and a variable neighborhood search scheme, thus completing the optimization of the post-earthquake repair sequence of the target substation. This invention not only optimizes the post-earthquake repair sequence of the equipment in the target substation but also achieves higher reliability and better accuracy.

Description

Technical Field

[0001] This invention belongs to the field of digital signal processing, specifically relating to a method and system for optimizing post-earthquake repair sequences in substations. Background Technology

[0002] With economic and technological development and the improvement of people's living standards, electricity has become an indispensable secondary energy source in people's production and daily life, bringing endless convenience. Therefore, ensuring a stable and reliable supply of electricity has become one of the most important tasks of the power system.

[0003] Substations are a critical component of the power system, and their rapid recovery capability after an earthquake is of paramount importance. During an earthquake, multiple aftershocks often follow the main shock, forming a mainshock-aftershock sequence. This sequence not only exacerbates damage to equipment in substations but may also interrupt the repair process, leading to a high degree of uncertainty in post-earthquake repair work.

[0004] Currently, traditional post-earthquake repair sequence optimization schemes for substations often only consider the impact of the mainshock while ignoring the influence of aftershocks; or they rely on empirical rules or static analysis based on a single mainshock, failing to fully consider the randomness and sequence characteristics of seismic action. Therefore, existing post-earthquake repair sequence optimization schemes for substations often suffer from low reliability and poor accuracy. Summary of the Invention

[0005] One of the objectives of this invention is to provide a highly reliable and accurate method for optimizing the post-earthquake repair sequence of substations.

[0006] The second objective of this invention is to provide a system for implementing the post-earthquake repair sequence optimization method for substations.

[0007] The post-earthquake repair sequence optimization method for substations provided by this invention includes the following steps:

[0008] S1. Obtain data information from the target substation;

[0009] S2. Based on the data obtained in step S1, construct the topology model of the target substation and determine the repair time and post-earthquake recovery process of each electrical device in the target substation.

[0010] S3. After the main shock of the earthquake occurs, acquire data on the damaged equipment of the target substation;

[0011] S4. Based on the data obtained in steps S2 and S3, perform post-earthquake repairs on the target substation equipment and update the equipment status in real time when aftershocks occur.

[0012] S5. Based on the data obtained in step S4, obtain the functional recovery curve of the target substation and calculate the resilience loss index of the target substation.

[0013] S6. Based on the resilience loss index obtained in step S5, construct a fitness function, and solve for the post-earthquake repair sequence of the target substation using a genetic algorithm and a variable neighborhood search scheme, thus completing the optimization of the post-earthquake repair sequence of the target substation.

[0014] Step S1 specifically includes the following steps:

[0015] Obtain data information about the target substation; the data information includes the transmission capacity of the target substation, the electrical nodes of the target substation, and the electrical equipment of the target substation.

[0016] Step S2 includes the following steps:

[0017] Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation;

[0018] The relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameters was established.

[0019] Set the repair time for each electrical device in the target substation;

[0020] Based on the time stepping method, the state recovery process of each electrical device in the target substation is set.

[0021] Step S2 specifically includes the following steps:

[0022] Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation. for , where nodes Represents electrical nodes and edges Represents electrical equipment, the j-th edge Corresponding dynamic capacity for , For the j-th edge The rated capacity of the corresponding electrical equipment, For the j-th edge The equipment status of the corresponding electrical equipment. The rules for determining the value are as follows: ;

[0023] Peak ground acceleration (PGA) is used as the seismic intensity parameter; a functional relationship is established between the failure probability of each electrical device in the target substation and the seismic intensity parameter, following the median value. The logarithmic standard deviation is The log-normal cumulative distribution;

[0024] Set the repair time for each electrical device in the target substation, following the median value. The logarithmic standard deviation is The log-normal cumulative distribution;

[0025] Based on the time-stepping method, the state recovery function of each electrical device in the target substation is set as follows:

[0026] In the formula For the i-th electrical device in the target substation at time... The state value, if but ,like but ; The length of the time step; Let be the repair rate of the i-th electrical device in the target substation; Let i be the maintenance instruction function for the i-th electrical device in the target substation. The rules for determining the value are as follows: .

[0027] Step S3 specifically includes the following steps:

[0028] After the main shock of the earthquake occurs, manual on-site investigation is conducted to obtain information on the status of damaged equipment in the target substation.

[0029] or

[0030] After the mainshock of the earthquake occurs, the peak ground acceleration corresponding to the mainshock is obtained. Based on the functional relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameter, the failure probability of each electrical device in the target substation is obtained. The obtained failure probability is compared with a generated random number: if the failure probability is less than or equal to the generated random number, the corresponding electrical device is determined to be in a normal state; if the failure probability is greater than the generated random number, the corresponding electrical device is determined to be in a damaged state. The generated random number is randomly generated and falls within the range of... A uniformly distributed random number.

[0031] Step S4 specifically includes the following steps:

[0032] After the main shock occurs, repairs will be carried out on all electrical equipment in the target substation.

[0033] After an aftershock occurs, the status of each electrical device in the target substation is updated using the following formula:

[0034] In the formula For the i-th electrical device in the target substation at time... The state after; the moment The time when the aftershocks occurred; For the i-th electrical device in the target substation at time... The previous state; Let i be the aftershock vulnerability of the i-th electrical equipment in the target substation. , Let be the cumulative probability density function. For the peak ground acceleration, Let this be the median of the number of damaged electrical devices. Let be the logarithmic standard deviation of the failure of the i-th electrical device.

[0035] Step S5 specifically includes the following steps:

[0036] Let n be the total number of devices to be repaired; repair sequence Represented as , This refers to the nth device that was repaired. To repair equipment Required repair time; For equipment The restored function values; This is the moment when the fault occurs and repairs begin;

[0037] Calculate the cumulative time to complete each repair action: before repair is complete The moment after the device It means that you thought , To repair equipment Required repair time;

[0038] Calculate the cumulative restored function value after each repair action is completed: before repair is complete The cumulative restored function value after each device for , For equipment The restored function values;

[0039] Then, a set of discrete function recovery points is obtained, denoted as ;

[0040] Plot time on the x-axis, with a time range from... to The function recovery curve of the target substation is obtained by plotting the cumulative restored function value on the vertical axis.

[0041] The toughness loss index of the target substation is calculated using the following formula. :

[0042] In the formula The time when the target substation's functions are restored to pre-earthquake levels; The initial power supply capacity of the target substation; Let be the power supply capacity of the target substation at time t.

[0043] Step S6 includes the following steps:

[0044] A fitness function is constructed based on the toughness loss index obtained in step S5.

[0045] A population is constructed based on the repair order of the damaged electrical equipment in the target substation. During the optimization process, a variable neighborhood search scheme is used for individual selection, and a genetic algorithm is used to optimize the repair order of the damaged electrical equipment in the target substation.

[0046] Finally, the post-earthquake repair sequence of the target substation was obtained, and the post-earthquake repair sequence optimization of the target substation was completed.

[0047] Step S6 specifically includes the following steps:

[0048] A. Initialize the population and use it as the parent population P: Number the damaged electrical devices and randomly sort them to generate the initial population; generate a total of Each population has a different initialization population; each population corresponds to a repair order.

[0049] B. Perform genetic operations on the parent population P to generate the offspring population O: calculate the toughness loss index corresponding to each parent population, and select several populations with the best toughness loss index. Through crossover and mutation, randomly change the repair order of several devices to generate new populations.

[0050] C. In the newly generated population from step B, calculate the resilience loss index for each population and retain the one with the smallest resilience loss index. Each population serves as the parent population for the next generation;

[0051] D. For the new generation of the population, the fitness value of each individual in the population is calculated using the following formula:

[0052] In the formula Let be the fitness value of the i1th individual; This is the resilience loss index for the i1th individual; This is a set minimum value to prevent the denominator from being 0;

[0053] For individuals whose fitness values ​​are better than the average fitness values ​​of all individuals in the population, a variable neighborhood search scheme is used. By switching between different neighborhood structures, optimization is performed within a set local range to improve the quality of the solution and avoid premature convergence.

[0054] The aforementioned variable neighborhood search scheme specifically includes the following steps:

[0055] D1. Set the selected individual to be represented as Set the current neighborhood number. =1;

[0056] D2. In the current neighborhood Randomly generate new solutions ;

[0057] D3. In of In the neighborhood, a local search scheme is used to search for a better solution, and the better solution obtained is taken as the new solution. The local search schemes mentioned include greedy search and 2-opt search schemes.

[0058] D4. If Then The value is updated to And keep the value of k unchanged; otherwise, keep The value of remains unchanged, and the value of k is increased by 1;

[0059] D5. Determine the value of k: if k is greater than 4, terminate the selection; otherwise, return to step D2 for iterative search.

[0060] The neighborhood includes four categories, defined as follows:

[0061] Exchange Neighborhood Randomly swap the positions of two devices;

[0062] Insert Neighborhood : Randomly select a device and insert it into a randomly set position, and the randomly set position is different from the original position of the selected device;

[0063] Inverted Neighborhood Select location in the current repair sequence. and , n1 is the length of the repair sequence, and The repair order for some parts is reversed; and The selection rule is random selection;

[0064] Block swap neighborhood : Set the block length to L, , This is a rounding function; in the current repair order, a first subsequence of length L and a second subsequence of length L are randomly selected, and the first and second subsequences are swapped;

[0065] E. Repeat steps B to D until the set iteration termination condition is met; output the final optimal population to obtain the corresponding post-earthquake repair sequence of the target substation, thus completing the post-earthquake repair sequence optimization of the target substation.

[0066] This invention also provides a system for implementing the post-earthquake repair sequence optimization method for substations, comprising a data acquisition module, a model building module, a damage confirmation module, a status update module, an index calculation module, and a sequence optimization module; the data acquisition module, model building module, damage confirmation module, status update module, index calculation module, and sequence optimization module are connected in series; the data acquisition module is used to acquire data information of the target substation and upload the data information to the model building module; the model building module is used to construct a topology model of the target substation based on the received data information and the acquired data information, and determine the repair time and post-earthquake status recovery process of each electrical device in the target substation, and upload the data information to the damage confirmation module; the damage confirmation module is used to determine the repair time and post-earthquake status recovery process of each electrical device in the target substation based on the received data information after the main shock of the earthquake. The system acquires data on damaged equipment in the target substation and uploads this data to the status update module. The status update module performs post-earthquake repairs on the equipment based on the received data, updates the equipment status in real-time during aftershocks, and uploads this data to the index calculation module. The index calculation module calculates the functional recovery curve of the target substation based on the received data and the obtained data, calculates the resilience loss index of the target substation, and uploads this data to the sequence optimization module. The sequence optimization module constructs a fitness function based on the obtained resilience loss index, and uses a genetic algorithm and a variable neighborhood search scheme to solve for the post-earthquake repair sequence of the target substation, thus optimizing the post-earthquake repair sequence of the target substation.

[0067] The post-earthquake repair sequence optimization method and system provided by this invention models the repair process of the target substation and equipment, and considers the effects of the mainshock and aftershocks of the earthquake process. By combining the functional recovery target of the target substation, genetic algorithm and variable neighborhood search scheme, it not only optimizes the repair sequence of the equipment of the target substation after the earthquake, but also achieves higher reliability and better accuracy. Attached Figure Description

[0068] Figure 1 This is a schematic diagram of the method flow of the present invention.

[0069] Figure 2 This is a schematic diagram of a substation model according to an embodiment of the method of the present invention.

[0070] Figure 3 This is a schematic diagram illustrating the prediction comparison results of an embodiment of the method of the present invention, wherein... Figure 3 (a) is a diagram comparing the repair times of various solutions. Figure 3 (b) is a comparative diagram of the toughness loss index of various schemes.

[0071] Figure 4 This is a schematic diagram illustrating the prediction comparison results considering parameter uncertainties in an embodiment of the method of the present invention, wherein... Figure 4 (a) A comparative diagram of the toughness loss indices of various schemes considering parameter uncertainties. Figure 4 (b) A comparative diagram of the repair times of various schemes considering parameter uncertainties.

[0072] Figure 5 This is a schematic diagram of the functional modules of the system of the present invention. Detailed Implementation

[0073] like Figure 1 The diagram shown is a flowchart of the method of the present invention: The post-earthquake repair sequence optimization method for substations disclosed in this invention includes the following steps:

[0074] S1. Obtain data information from the target substation; specifically including the following steps:

[0075] Acquire data information about the target substation; the data information includes the transmission capacity of the target substation, the electrical nodes of the target substation, and the electrical equipment of the target substation;

[0076] S2. Based on the data obtained in step S1, construct a topology model of the target substation and determine the repair time and post-earthquake recovery process of each electrical device in the target substation; including the following steps:

[0077] Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation;

[0078] The relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameters was established.

[0079] Set the repair time for each electrical device in the target substation;

[0080] Based on the time-stepping method, the state recovery process of each electrical device in the target substation is defined;

[0081] In practice, the following steps can be taken:

[0082] Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation. for , where nodes Represents electrical nodes (including busbars, connection points, etc.), edges This represents the j-th edge of an electrical equipment (including circuit breakers, transformers, and overhead contact lines). Corresponding dynamic capacity for , For the j-th edge The rated capacity of the corresponding electrical equipment, For the j-th edge The equipment status of the corresponding electrical equipment. The rules for determining the value are as follows: ;

[0083] Peak ground acceleration (PGA) is used as the seismic intensity parameter; a functional relationship is established between the failure probability of each electrical device in the target substation and the seismic intensity parameter, following the median value. The logarithmic standard deviation is The log-normal cumulative distribution;

[0084] Set the repair time for each electrical device in the target substation, following the median value. The logarithmic standard deviation is The log-normal cumulative distribution;

[0085] Based on the time-stepping method, the state recovery function of each electrical device in the target substation is set as follows:

[0086] In the formula For the i-th electrical device in the target substation at time... The state value, if but ,like but ; The length of the time step; Let be the repair rate of the i-th electrical device in the target substation; Let i be the maintenance instruction function for the i-th electrical device in the target substation. The rules for determining the value are as follows: ;

[0087] In practical implementation, the state value A value of 1 indicates that the equipment is under maintenance; status value. A value of 0 indicates that the equipment has not been repaired;

[0088] S3. After the main shock of the earthquake, acquire data on the damaged equipment of the target substation; specifically including the following steps:

[0089] After the main shock of the earthquake occurs, manual on-site investigation is conducted to obtain information on the status of damaged equipment in the target substation.

[0090] or

[0091] After the mainshock of the earthquake occurs, the peak ground acceleration corresponding to the mainshock is obtained. Based on the functional relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameter, the failure probability of each electrical device in the target substation is obtained. The obtained failure probability is compared with a generated random number: if the failure probability is less than or equal to the generated random number, the corresponding electrical device is determined to be in a normal state; if the failure probability is greater than the generated random number, the corresponding electrical device is determined to be in a damaged state. The generated random number is randomly generated and falls within the range of... A uniformly distributed random number;

[0092] S4. Based on the data obtained in steps S2 and S3, perform post-earthquake repairs on the target substation equipment and update the equipment status in real time when aftershocks occur; specifically including the following steps:

[0093] After the main shock occurs, repairs will be carried out on all electrical equipment in the target substation.

[0094] After an aftershock occurs, the status of each electrical device in the target substation is updated using the following formula:

[0095] In the formula For the i-th electrical device in the target substation at time... The state after; the moment The time when the aftershocks occurred; For the i-th electrical device in the target substation at time... The previous state; Let i be the aftershock vulnerability of the i-th electrical equipment in the target substation. , Let be the cumulative probability density function. For the peak ground acceleration, Let this be the median of the number of damaged electrical devices. Let be the logarithmic standard deviation of the failure of the i-th electrical device;

[0096] S5. Based on the data obtained in step S4, obtain the functional recovery curve of the target substation and calculate the resilience loss index of the target substation; specifically including the following steps:

[0097] Let n be the total number of devices to be repaired; repair sequence Represented as , This refers to the nth device that was repaired. To repair equipment Required repair time; For equipment The restored function values; This is the moment when the fault occurs and repairs begin;

[0098] Calculate the cumulative time to complete each repair action: before repair is complete The moment after the device It means that you thought , To repair equipment Required repair time;

[0099] Calculate the cumulative restored function value after each repair action is completed: before repair is complete The cumulative restored function value after each device for , For equipment The restored function values;

[0100] Then, a set of discrete function recovery points is obtained, denoted as ;

[0101] Plot time on the x-axis, with a time range from... to The function recovery curve of the target substation is obtained by plotting the cumulative restored function value on the vertical axis.

[0102] The toughness loss index of the target substation is calculated using the following formula. :

[0103] In the formula The time when the target substation's functions are restored to pre-earthquake levels; The initial power supply capacity of the target substation; Let be the power supply capacity of the target substation at time t;

[0104] S6. Based on the resilience loss index obtained in step S5, construct a fitness function, and solve for the post-earthquake repair sequence of the target substation using a genetic algorithm and a variable neighborhood search scheme, thus completing the optimization of the post-earthquake repair sequence of the target substation; including the following steps:

[0105] A fitness function is constructed based on the toughness loss index obtained in step S5.

[0106] A population is constructed based on the repair order of the damaged electrical equipment in the target substation. During the optimization process, a variable neighborhood search scheme is used for individual selection, and a genetic algorithm is used to optimize the repair order of the damaged electrical equipment in the target substation.

[0107] Finally, the post-earthquake repair sequence of the target substation was obtained, and the post-earthquake repair sequence optimization of the target substation was completed;

[0108] In practice, the following steps can be taken:

[0109] A. Initialize the population and use it as the parent population P: Number the damaged electrical devices and randomly sort them to generate the initial population; generate a total of Each population has a different initialization population; each population corresponds to a repair order.

[0110] B. Perform genetic operations on the parent population P to generate the offspring population O: calculate the toughness loss index corresponding to each parent population, and select several populations with the best toughness loss index. Through crossover and mutation, randomly change the repair order of several devices to generate new populations.

[0111] C. In the newly generated population from step B, calculate the resilience loss index for each population and retain the one with the smallest resilience loss index. Each population serves as the parent population for the next generation;

[0112] D. For the new generation of the population, the fitness value of each individual in the population is calculated using the following formula:

[0113] In the formula Let be the fitness value of the i1th individual; This is the resilience loss index for the i1th individual; This is a set minimum value to prevent the denominator from being 0;

[0114] For individuals whose fitness values ​​are better than the average fitness values ​​of all individuals in the population, a variable neighborhood search scheme is used. By switching between different neighborhood structures, optimization is performed within a set local range to improve the quality of the solution and avoid premature convergence.

[0115] The aforementioned variable neighborhood search scheme specifically includes the following steps:

[0116] D1. Set the selected individual to be represented as Set the current neighborhood number. =1;

[0117] D2. In the current neighborhood Randomly generate new solutions ;

[0118] D3. In of In the neighborhood, a local search scheme is used to search for a better solution, and the better solution obtained is taken as the new solution. The local search schemes mentioned include greedy search and 2-opt search schemes.

[0119] D4. If Then The value is updated to And keep the value of k unchanged; otherwise, keep The value of remains unchanged, and the value of k is increased by 1;

[0120] D5. Determine the value of k: if k is greater than 4, terminate the selection; otherwise, return to step D2 for iterative search.

[0121] The neighborhood includes four categories, defined as follows:

[0122] Exchange Neighborhood Randomly swap the positions of two devices;

[0123] Insert Neighborhood : Randomly select a device and insert it into a randomly set position, and the randomly set position is different from the original position of the selected device;

[0124] Inverted Neighborhood Select location in the current repair sequence. and , n1 is the length of the repair sequence, and The repair order for some parts is reversed; and The selection rule is random selection;

[0125] Block swap neighborhood : Set the block length to L, , This is a rounding function; in the current repair order, a first subsequence of length L and a second subsequence of length L are randomly selected, and the first and second subsequences are swapped;

[0126] E. Repeat steps B to D until the set iteration termination condition is met; output the final optimal population to obtain the corresponding post-earthquake repair sequence of the target substation, thus completing the post-earthquake repair sequence optimization of the target substation.

[0127] This invention establishes a directed graph logic model of the substation system, integrating the seismic vulnerability of equipment, probabilistic repair time, and the dynamic coupling effects of the main shock and aftershocks, thus achieving a refined simulation of the post-earthquake system functional recovery trajectory. Furthermore, this invention constructs a hybrid intelligent optimization algorithm that combines the global search capability of a genetic algorithm with the local optimization capability of a variable neighborhood search. This algorithm can efficiently solve for the optimal equipment repair order that minimizes system resilience loss and maximizes recovery from a massive number of possible sequences. Compared with traditional empirical strategies or other single algorithms, this invention significantly improves post-earthquake recovery efficiency and power supply reliability, and exhibits excellent robustness and stability under parameter uncertainty conditions, providing scientific and reliable quantitative decision support for post-earthquake emergency resource scheduling in power systems.

[0128] The method of the present invention will be further described below with reference to an embodiment:

[0129] by Figure 2 Taking a typical substation as an example, it includes two parallel "220 kV incoming line - main transformer - 27.5 kV outgoing line" units, one as the main power supply and the other as a backup power supply. Each unit mainly includes the following equipment: circuit breaker (QF), disconnector (QS), current transformer (TA), voltage transformer (TV), main transformer (T), and surge arrester (BL). The No. 1 and No. 2 incoming lines of this substation operate in a mutual backup mode, and are equipped with backup transformer #2T and mutual backup transformer #4T. A directed graph logic model of the substation system is established to illustrate the layout and logical relationships of the electrical equipment in the substation. Since the surge arrester does not affect current transmission in the equipment, it is not included in the logic model.

[0130] The vulnerability parameters of the substation are shown in Table 1.

[0131] The substation repair time parameters are shown in Table 2.

[0132] After the main shock occurs, the damaged equipment in the substation is located and collected. In this embodiment, the equipment seismic vulnerability curve is used to obtain the failure probability P (0.4g) of the equipment under the current seismic input intensity. The relative magnitude of this probability P is then compared with a randomly generated number uniformly distributed between 0 and 1 to determine the working status of the equipment. The damaged equipment is shown in Table 3.

[0133] Substation equipment repair began. After an aftershock occurred (the peak ground acceleration was 0.3g, and the aftershock occurred 25 days after the main shock), the equipment damage status was updated in real time, and the newly damaged equipment caused by the aftershock was shown in Table 4.

[0134] To verify the advantages of the method of this invention, the scheme of this invention (labeled "GA-VNS") is compared with three commonly used strategies: stochastic recovery method (RS), greedy algorithm recovery method (GD), and simulated annealing algorithm recovery method (SA). Based on the simulated system function recovery curve (function-time curve), the resilience loss index is calculated; the recovery time of the substation system under different algorithms is statistically obtained, and the results are as follows: Figure 3 As shown, the proposed solution (GA-VNS) is significantly superior to other algorithms, and the repair sequences of each algorithm are shown in Table 5.

[0135] Furthermore, to verify the stability of the present invention (GA-VNS) under uncertain conditions, this embodiment considers the uncertainties in the quantity of spare parts for each piece of equipment and the efficiency of the construction team, and performs 2000 random simulations. The results are as follows: Figure 4 As shown. (Through) Figure 4 As can be seen, the proposed solution (GA-VNS) is significantly superior to other algorithms in terms of both substation resilience loss and repair time.

[0136] like Figure 5The diagram shows the functional modules of the system of this invention: The system disclosed in this invention for implementing the post-earthquake repair sequence optimization method for substations includes a data acquisition module, a model building module, a damage confirmation module, a status update module, an index calculation module, and a sequence optimization module; these modules are connected in series. The data acquisition module acquires data information of the target substation and uploads it to the model building module. The model building module constructs a topology model of the target substation based on the received and acquired data information, determines the repair time and post-earthquake status recovery process of each electrical device in the target substation, and uploads the data information to the damage confirmation module. The damage confirmation module, based on the received data information, when... After the main shock of the earthquake, data on the damaged equipment of the target substation is acquired and uploaded to the status update module. The status update module performs post-earthquake repairs on the equipment based on the received data and updates the equipment status in real time during aftershocks, uploading the data to the index calculation module. The index calculation module calculates the functional recovery curve of the target substation based on the received data and the obtained data, and calculates the resilience loss index of the target substation, uploading the data to the sequence optimization module. The sequence optimization module constructs a fitness function based on the obtained resilience loss index, and solves for the post-earthquake repair sequence of the target substation using a genetic algorithm and a variable neighborhood search scheme, thus completing the post-earthquake repair sequence optimization of the target substation.

Claims

1. A method for optimizing the post-earthquake repair sequence of a substation, characterized in that... Includes the following steps: S1. Obtain data information from the target substation; S2. Based on the data obtained in step S1, construct a topology model of the target substation and determine the repair time and post-earthquake recovery process of each electrical device in the target substation; including the following steps: Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation; The relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameters was established. Set the repair time for each electrical device in the target substation; Based on the time-stepping method, the state recovery process of each electrical device in the target substation is defined; S3. After the main shock of the earthquake occurs, acquire data on the damaged equipment of the target substation; S4. Based on the data obtained in steps S2 and S3, perform post-earthquake repairs on the target substation equipment and update the equipment status in real time when aftershocks occur. S5. Based on the data obtained in step S4, obtain the functional recovery curve of the target substation and calculate the resilience loss index of the target substation. S6. Based on the resilience loss index obtained in step S5, construct a fitness function, and solve for the post-earthquake repair sequence of the target substation using a genetic algorithm and a variable neighborhood search scheme, thus completing the optimization of the post-earthquake repair sequence of the target substation.

2. The substation post-earthquake repair sequence optimization method according to claim 1, characterized in that... Step S1 specifically includes the following steps: Obtain data information about the target substation; the data information includes the transmission capacity of the target substation, the electrical nodes of the target substation, and the electrical equipment of the target substation.

3. The substation post-earthquake repair sequence optimization method according to claim 2, characterized in that... Step S2 specifically includes the following steps: Based on the layout and logical relationships of the electrical equipment in the target substation, construct a directed graph model of the target substation. for , where nodes Represents electrical nodes and edges. Represents electrical equipment, the j-th edge Corresponding dynamic capacity for , For the j-th edge The rated capacity of the corresponding electrical equipment, For the j-th edge The equipment status of the corresponding electrical equipment. The rules for determining the value are as follows: ; Peak ground acceleration (PGA) is used as the seismic intensity parameter; a functional relationship is established between the failure probability of each electrical device in the target substation and the seismic intensity parameter, following a median value. The logarithmic standard deviation is The log-normal cumulative distribution; Set the repair time for each electrical device in the target substation, following the median value. The logarithmic standard deviation is The log-normal cumulative distribution; Based on the time-stepping method, the state recovery function of each electrical device in the target substation is set as follows: In the formula For the i-th electrical device in the target substation at time... The state value, if but ,like but ; The length of the time step; Let be the repair rate of the i-th electrical device in the target substation; Let i be the maintenance instruction function for the i-th electrical device in the target substation. The rules for determining the value are as follows: .

4. The substation post-earthquake repair sequence optimization method according to claim 3, characterized in that... Step S3 specifically includes the following steps: After the main shock of an earthquake occurs, a manual on-site investigation is conducted to obtain information on the status of damaged equipment in the target substation.

5. The substation post-earthquake repair sequence optimization method according to claim 3, characterized in that... Step S3 specifically includes the following steps: After the mainshock of the earthquake occurs, the peak ground acceleration corresponding to the mainshock is obtained. Based on the functional relationship between the failure probability of each electrical device in the target substation and the seismic intensity parameter, the failure probability of each electrical device in the target substation is obtained. The obtained failure probability is compared with a generated random number: if the failure probability is less than or equal to the generated random number, the corresponding electrical device is determined to be in a normal state; if the failure probability is greater than the generated random number, the corresponding electrical device is determined to be in a damaged state. The generated random number is randomly generated and falls within the range of... A uniformly distributed random number.

6. The substation post-earthquake repair sequence optimization method according to claim 5, characterized in that... Step S4 specifically includes the following steps: After the main shock occurs, repairs will be carried out on all electrical equipment in the target substation. After an aftershock occurs, the status of each electrical device in the target substation is updated using the following formula: In the formula For the i-th electrical device in the target substation at time... The state after; the moment The time when the aftershocks occurred; For the i-th electrical device in the target substation at time... The previous state; Let i be the aftershock vulnerability of the i-th electrical equipment in the target substation. , Let be the cumulative probability density function. For the peak ground acceleration, Let this be the median of the number of damaged electrical devices. Let be the logarithmic standard deviation of the failure of the i-th electrical device.

7. The substation post-earthquake repair sequence optimization method according to claim 6, characterized in that... Step S5 specifically includes the following steps: Let n be the total number of devices to be repaired; repair sequence Represented as , This refers to the nth device that was repaired. To repair equipment Required repair time; For equipment The restored function values; This is the moment when the fault occurs and repairs begin; Calculate the cumulative time to complete each repair action: before repair is complete The moment after the device It means that you thought , To repair equipment Required repair time; Calculate the cumulative restored function value after each repair action is completed: before repair is complete The cumulative restored function value after each device for , For equipment The restored function values; Then, a set of discrete function recovery points is obtained, denoted as ; Plot time on the x-axis, with a time range from... to The function recovery curve of the target substation is obtained by plotting the cumulative restored function value on the vertical axis. The toughness loss index of the target substation is calculated using the following formula. : In the formula The time when the target substation's functions are restored to pre-earthquake levels; The initial power supply capacity of the target substation; Let be the power supply capacity of the target substation at time t.

8. The substation post-earthquake repair sequence optimization method according to claim 7, characterized in that... Step S6 includes the following steps: A fitness function is constructed based on the toughness loss index obtained in step S5. A population is constructed based on the repair order of the damaged electrical equipment in the target substation. During the optimization process, a variable neighborhood search scheme is used for individual selection, and a genetic algorithm is used to optimize the repair order of the damaged electrical equipment in the target substation. Finally, the post-earthquake repair sequence of the target substation was obtained, and the post-earthquake repair sequence optimization of the target substation was completed.

9. The method for optimizing the post-earthquake repair sequence of a substation according to claim 8, characterized in that... Step S6 specifically includes the following steps: A. Initialize the population and use it as the parent population P: Number the damaged electrical devices and randomly sort them to generate the initial population; generate a total of Each population has a different initialization population; each population corresponds to a repair order. B. Perform genetic operations on the parent population P to generate the offspring population O: calculate the toughness loss index corresponding to each parent population, and select several populations with the best toughness loss index. Through crossover and mutation, randomly change the repair order of several devices to generate new populations. C. In the newly generated population from step B, calculate the resilience loss index for each population and retain the one with the smallest resilience loss index. Each population serves as the parent population for the next generation; D. For the new generation of the population, the fitness value of each individual in the population is calculated using the following formula: In the formula Let be the fitness value of the i1th individual; This is the resilience loss index for the i1th individual; This is a set minimum value to prevent the denominator from being 0; For individuals whose fitness values ​​are better than the average fitness values ​​of all individuals in the population, a variable neighborhood search scheme is used. By switching between different neighborhood structures, optimization is performed within a set local range to improve the quality of the solution and avoid premature convergence. The aforementioned variable neighborhood search scheme specifically includes the following steps: D1. Set the selected individual to be represented as Set the current neighborhood number. =1; D2. In the current neighborhood Randomly generate new solutions ; D3. In of In the neighborhood, a local search scheme is used to search for a better solution, and the better solution obtained is taken as the new solution. ; D4. If Then The value is updated to And keep the value of k unchanged; otherwise, keep The value of remains unchanged, and the value of k is increased by 1; D5. Determine the value of k: if k is greater than 4, terminate the selection; otherwise, return to step D2 for iterative search. The neighborhood includes four categories, defined as follows: Exchange Neighborhood Randomly swap the positions of two devices; Insert Neighborhood : Randomly select a device and insert it into a randomly set position, and the randomly set position is different from the original position of the selected device; Inverted Neighborhood Select location in the current repair sequence. and , n1 is the length of the repair sequence, and The repair order for some parts is reversed; and The selection rule is random selection; Block swap neighborhood Set the block length to L. , This is a rounding function; in the current repair order, a first subsequence of length L and a second subsequence of length L are randomly selected, and the first and second subsequences are swapped; E. Repeat steps B to D until the set iteration termination condition is met; output the final optimal population to obtain the corresponding post-earthquake repair sequence of the target substation, thus completing the post-earthquake repair sequence optimization of the target substation.

10. A system for implementing the substation post-earthquake repair sequence optimization method according to any one of claims 1 to 9, characterized in that... It includes a data acquisition module, a model building module, a damage confirmation module, a status update module, an index calculation module, and a sequence optimization module; these modules are connected in series. The data acquisition module acquires data information from the target substation and uploads it to the model building module. The model building module constructs a topology model of the target substation based on the received and acquired data information, determines the repair time and post-earthquake recovery process of each electrical device in the target substation, and uploads the data to the damage confirmation module. The damage confirmation module is used to obtain data information on the damaged equipment of the target substation after the main shock of the earthquake, based on the received data information, and upload the data information to the status update module. The status update module is used to perform post-earthquake repair of the target substation equipment based on the received data information, update the equipment status in real time when aftershocks occur, and upload the data information to the index calculation module. The index calculation module is used to obtain the functional recovery curve of the target substation based on the received data information, calculate the resilience loss index of the target substation, and upload the data information to the sequence optimization module. The sequence optimization module is used to construct a fitness function based on the received data information and the obtained resilience loss index. Based on the genetic algorithm and the variable neighborhood search scheme, it solves the post-earthquake repair sequence of the target substation and completes the post-earthquake repair sequence optimization of the target substation.

Images