Water power-grid power flow coupling simulation method and device based on multi-time scale and HLA architecture and computer equipment

Abstract

The application relates to a water power-grid power flow coupling simulation method and device based on a multi-time scale and HLA architecture, computer equipment, a computer readable storage medium and a computer program product. The method comprises the following steps: receiving a water power parameter matrix for a region to be analyzed at a current simulation time forwarded by an RTI; determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed according to the water power parameter matrix, determining the equipment state of the power grid equipment based on the failure probability, updating the power grid topology based on the equipment state, and obtaining the updated power grid topology; obtaining the power grid power flow calculation result of the power grid equipment in the updated power grid topology, performing protection checking based on the power grid power flow calculation result, obtaining the protection checking result, and updating the updated power grid topology again based on the protection checking result until there is no new equipment state change at the current simulation time. By using the method, the accuracy of the simulation evaluation of the flood risk of the power grid can be improved.

Description

Technical Field

[0001] This application relates to the field of power grid technology, and in particular to a hydrodynamic-power grid flow coupling simulation method, device, computer equipment, computer-readable storage medium, and computer program product based on a multi-timescale and HLA architecture. Background Technology

[0002] With the frequent occurrence of extreme rainstorms and floods, a large number of power distribution equipment such as substations, distribution cabinets, and poles in urban areas are susceptible to damage from water immersion and water flow, which can lead to equipment failure, line tripping, and even large-scale cascading power outages. Therefore, it is necessary to carry out hydrodynamic and power flow coupling simulation to dynamically assess the risks of power grid operation under flood scenarios.

[0003] In traditional techniques, offline serial simulation is mainly used when conducting hydrodynamic-power flow coupling simulation. However, offline serial simulation can only calculate the power grid operating state once at a single hydrodynamic simulation moment, and cannot simulate the continuous evolution of power grid cascading faults at a small step scale between two hydrodynamic moments, resulting in poor simulation integrity. Moreover, the data structures, time scales, and spatial scales of the hydrodynamic model and the power flow model are inconsistent. Simple offline import cannot accurately represent the temporal correlation between the evolution of flood disasters and the triggering of power grid faults, ultimately causing the simulation results of power grid fault evolution under flood scenarios to be distorted, which significantly reduces the accuracy of power grid flood risk simulation assessment. Summary of the Invention

[0004] Therefore, it is necessary to provide a hydrodynamic-power flow coupling simulation method, device, computer equipment, computer-readable storage medium, and computer program product based on a multi-timescale and HLA architecture that can improve the accuracy of power grid flood risk simulation assessment, addressing the aforementioned technical problems.

[0005] Firstly, this application provides a hydrodynamic-grid power flow coupling simulation method based on a multi-timescale and HLA architecture, applied to power federation members in the HLA architecture; the HLA architecture also includes hydrodynamic federation members and RTIs, wherein the power federation members are registered by a grid power flow calculation program, and the hydrodynamic federation members are registered by a hydrodynamic simulation program; including:

[0006] The system receives the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time, forwarded by the RTI. The hydrodynamic federation member is used to perform simulation calculations on the flood inundation process of the region to be analyzed according to a preset hydrodynamic calculation step size, to obtain the hydrodynamic parameter matrix carrying the corresponding simulation time, and publishes the hydrodynamic parameter matrix to the RTI.

[0007] Based on the hydrodynamic parameter matrix, the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed is determined. Based on the failure probability, the equipment status of the power grid equipment is determined. Based on the equipment status, the power grid topology is updated to obtain the updated power grid topology.

[0008] Obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain the protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time;

[0009] The sub-update simulation time of the current simulation time is determined based on the preset power flow calculation step size, and then the process jumps to the step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0010] The next simulation moment is taken as the new current simulation moment, and the process jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation moment forwarded by RTI, until the simulation end condition is met.

[0011] In one embodiment, determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix includes:

[0012] The hydrodynamic parameter matrix is ​​analyzed to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed;

[0013] According to the reference coordinate system corresponding to the power grid topology, the coordinate information of each hydrodynamic calculation grid is transformed to obtain the transformed coordinate information of each hydrodynamic calculation grid.

[0014] Based on the transformed coordinate information of each hydrodynamic calculation grid, the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology is determined;

[0015] Based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, the hydrodynamic parameters of the location of the power grid equipment are determined.

[0016] The failure probability of the power grid equipment is determined based on the hydrodynamic parameters of its location.

[0017] In one embodiment, determining the hydrodynamic parameters of the location of the power grid equipment based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology includes:

[0018] When the number of target hydrodynamic calculation grids corresponding to the power grid equipment is one, the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment are used as the hydrodynamic parameters of the location of the power grid equipment.

[0019] When there are multiple target hydrodynamic calculation grids corresponding to the power grid equipment, the hydrodynamic parameters of each target hydrodynamic calculation grid corresponding to the power grid equipment are fused according to the superposition area of ​​the power grid equipment and each target hydrodynamic calculation grid to obtain the hydrodynamic parameters of the location of the power grid equipment.

[0020] In one embodiment, determining the failure probability of the power grid equipment based on the hydrodynamic parameters of the location of the power grid equipment includes:

[0021] Based on the equipment type of the power grid equipment, the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment are determined; the vulnerability parameter is used to characterize the influence weight of the hydrodynamic parameters on the failure probability.

[0022] The failure probability of the power grid equipment is obtained by inputting the hydrodynamic parameters, the hydrodynamic parameter thresholds, and the vulnerability parameters into a pre-built failure probability statistical model.

[0023] In one embodiment, determining the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment based on the equipment type of the power grid equipment includes:

[0024] Based on the equipment type of the power grid equipment, the correspondence between the equipment type and the hydrodynamic parameter threshold is queried to obtain the hydrodynamic parameter threshold corresponding to the equipment type, which is then used as the hydrodynamic parameter threshold corresponding to the location of the power grid equipment.

[0025] Based on the equipment type of the power grid equipment, the correspondence between the equipment type and the vulnerability parameter is queried to obtain the vulnerability parameter corresponding to the equipment type, which is then used as the vulnerability parameter corresponding to the hydrodynamic parameter of the location of the power grid equipment.

[0026] In one embodiment, determining the equipment status of the power grid equipment based on the failure probability includes:

[0027] If the failure probability is greater than a preset probability, the power grid equipment is determined to be in a failure state.

[0028] If the failure probability is less than or equal to the preset probability, the power grid equipment is determined to be in a normal state.

[0029] Secondly, this application also provides a hydrodynamic-grid power flow coupling simulation device based on a multi-timescale and HLA architecture, applied to power federation members in the HLA architecture; the HLA architecture also includes hydrodynamic federation members and RTIs, wherein the power federation members are registered by a grid power flow calculation program, and the hydrodynamic federation members are registered by a hydrodynamic simulation program; including:

[0030] Force parameter matrix; The hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the area to be analyzed according to a preset hydrodynamic calculation step size, to obtain a hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI;

[0031] The topology update module is used to determine the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, and update the power grid topology based on the equipment status to obtain the updated power grid topology.

[0032] The first processing module is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0033] The second processing module is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0034] The third processing module is used to take the next simulation time as the new current simulation time and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

[0035] Thirdly, this application also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0036] Force parameter matrix; The hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the area to be analyzed according to a preset hydrodynamic calculation step size, to obtain a hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI;

[0037] The topology update module is used to determine the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, and update the power grid topology based on the equipment status to obtain the updated power grid topology.

[0038] The first processing module is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0039] The second processing module is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0040] The third processing module is used to take the next simulation time as the new current simulation time and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

[0041] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, performs the following steps:

[0042] Force parameter matrix; The hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the area to be analyzed according to a preset hydrodynamic calculation step size, to obtain a hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI;

[0043] The topology update module is used to determine the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, and update the power grid topology based on the equipment status to obtain the updated power grid topology.

[0044] The first processing module is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0045] The second processing module is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0046] The third processing module is used to take the next simulation time as the new current simulation time and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

[0047] Fifthly, this application also provides a computer program product, including a computer program that, when executed by a processor, performs the following steps:

[0048] Force parameter matrix; The hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the area to be analyzed according to a preset hydrodynamic calculation step size, to obtain a hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI;

[0049] The topology update module is used to determine the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, and update the power grid topology based on the equipment status to obtain the updated power grid topology.

[0050] The first processing module is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0051] The second processing module is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0052] The third processing module is used to take the next simulation time as the new current simulation time and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

[0053] The aforementioned hydrodynamic-power grid flow coupling simulation method, apparatus, computer equipment, computer-readable storage medium, and computer program product based on multi-timescale and HLA architecture first receive the hydrodynamic parameter matrix for the region to be analyzed at the current simulation moment, forwarded by the RTI. Hydrodynamic federation members then perform simulation calculations on the flood inundation process of the region to be analyzed according to a preset hydrodynamic calculation step size, obtaining a hydrodynamic parameter matrix carrying the corresponding simulation moment, and publish the hydrodynamic parameter matrix to the RTI. Next, based on the hydrodynamic parameter matrix, the failure probability of power grid equipment in the corresponding power grid topology of the region to be analyzed is determined. Based on the failure probability, the equipment status of the power grid equipment is determined, and the power grid topology is updated based on the equipment status to obtain the updated power grid topology. Then, more... The power flow calculation results of the power grid equipment in the new power grid topology are used to perform protection verification based on the power flow calculation results. The updated power grid topology is then updated again based on the protection verification results until no new equipment status changes at the current simulation time. Next, the sub-update simulation time of the current simulation time is determined based on the preset power flow calculation step size, and the process jumps to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time is reached. Finally, the next simulation time is taken as the new current simulation time, and the process jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time forwarded by RTI, until the simulation termination condition is met. In this way, under the HLA architecture, the hydrodynamic simulation program and the power flow calculation program are coupled and simulated based on RTI. The hydrodynamic parameter matrix for the flood inundation process of the area to be analyzed is updated by pre-setting the hydrodynamic calculation step size. At the same time, the sub-update simulation time is obtained by using the pre-set power flow calculation step size within the adjacent simulation time for hydrodynamics. The entire process of power grid equipment failure, topology update and protection verification is continuously iterated and deduced. It can completely reproduce the dynamic evolution details of the power grid cascading faults between two adjacent simulation time points, which makes up for the shortcomings of traditional offline serial simulation, which can only calculate the power grid state at a single point and thus has poor simulation integrity. At the same time, the transmission and reception of the hydrodynamic parameter matrix is ​​realized through RTI, which solves the problem of time series correlation distortion caused by the mismatch of multi-model data structure and spatiotemporal scale. It accurately matches the hydrodynamic parameters corresponding to the power grid equipment and quantifies the equipment failure probability, effectively avoiding the problem of distortion of traditional simulation results. This makes the simulation results of the power grid fault evolution process under flood scenarios more in line with the actual working conditions, thereby improving the accuracy of power grid flood risk simulation assessment. Attached Figure Description

[0054] To more clearly illustrate the technical solutions in the embodiments of this application or related technologies, the drawings used in the description of the embodiments of this application or related technologies will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0055] Figure 1 This is an application environment diagram of a hydrodynamic-power grid coupling simulation method based on multi-timescale and HLA architecture in one embodiment;

[0056] Figure 2 This is a flowchart illustrating a hydrodynamic-power grid flow coupling simulation method based on a multi-timescale and HLA architecture in one embodiment.

[0057] Figure 3 This is a flowchart illustrating the steps for determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed in one embodiment.

[0058] Figure 4 This is a structural block diagram of a hydrodynamic-electric power flow coupling simulation device based on a multi-timescale and HLA architecture in one embodiment.

[0059] Figure 5 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0060] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0061] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0062] Urban flooding and river overflows can alter the operating environment of power distribution network equipment. Substations, ring main units, switch stations, low-voltage distribution boxes, and some tower foundations may experience insulation degradation, control circuit failures, protection malfunctions, or manual shutdowns under conditions of water accumulation, scouring, and prolonged immersion. After equipment is decommissioned, the original load will be redistributed along available lines, and some feeders may experience overload or voltage exceeding limits, which may trigger further tripping.

[0063] Therefore, traditional technologies mainly suffer from the following technical problems: 1. Difficulty in exchanging data between heterogeneous models in chronological order: Hydrodynamic models often use two-dimensional grids or irregular meshes to represent water depth and flow velocity, with calculation steps typically on the order of minutes; Power flow models take buses, branches, and switches as objects, and state changes and protection checks often require steps on the order of seconds or even smaller; the data structures, time scales, and spatial scales of the two are inconsistent, and simple offline import is insufficient to express the sequential relationships during a disaster. 2. Overly simplistic equipment failure criteria: Some schemes only use water depth thresholds to determine whether equipment has failed, for example, determining shutdown when water exceeds a certain height, without considering flow velocity impact, duration of flooding, and differences between different equipment types; such processing can estimate the number of flooded devices, but it is difficult to describe the gradual changes in equipment status during the flood process. 3. The scope of power outage impact is often treated as a static result: In actual operation, the exit of a node may cause a power flow shift, which may in turn cause other lines or equipment to exceed limits; if the simulation process does not place equipment failure, topology update, power flow recalculation, and protection actions in the same loop, it is difficult to analyze the successive tripping and load transfer processes.

[0064] Therefore, a method is needed to organize hydrodynamic processes, grid topology states, and power flow calculations along a unified simulation timeline, enabling the hydrodynamic calculation results to be used by the grid model in a timely manner and to continue power flow and protection verification after equipment status changes. Based on this, this application provides a hydrodynamic-grid power flow coupled simulation method based on multiple timescales and an HLA (High Level Architecture) architecture. This method can describe the mutual influence between flood evolution, equipment failure, and power flow transfer on the same simulation timeline, and is suitable for analyzing the power outage range of distribution networks and comparing emergency dispatch schemes under flood disasters.

[0065] The hydrodynamic-power grid coupling simulation method based on multi-timescale and HLA architecture provided in this application can be applied to, for example... Figure 1 The application environment shown is illustrated. The HLA architecture includes RTI (Run-Time Infrastructure) 101, a hydrodynamic federation member 102, and an electric federation member 103. RTI 101 is connected to both hydrodynamic federation member 102 and electric federation member 103. For details, refer to... Figure 1Hydrodynamic Federation member 102 performs simulation calculations of the flood inundation process in the area to be analyzed according to the preset hydrodynamic calculation step size, obtains the hydrodynamic parameter matrix carrying the corresponding simulation time, and publishes the hydrodynamic parameter matrix to RTI101. RTI101 sends the hydrodynamic parameter matrix for the area to be analyzed at the current simulation time to the power federation member 103. Based on the hydrodynamic parameter matrix, member 103 of the power federation determines the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed. Based on the failure probability, it determines the equipment status of the power grid equipment and updates the power grid topology based on the equipment status, obtaining the updated power grid topology. It obtains the power flow calculation results of the power grid equipment in the updated power grid topology, performs protection verification based on the power flow calculation results, obtains the protection verification results, and updates the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time. Based on the preset power flow calculation step size, it determines the sub-update simulation time of the current simulation time and jumps to the step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time is reached. The next simulation time is taken as the new current simulation time, and it jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time forwarded by RTI, until the simulation termination condition is met.

[0066] The HLA architecture refers to a standardized distributed simulation framework that separates a hydrodynamic simulation program and a power flow calculation program—which are completely different in principle, step size, and software—into two independent federation members. These two members do not require modification of the underlying code or customization of interface; they only communicate and synchronize data across multiple time scales through the middleware RTI. RTI101 is the core middleware in the HLA architecture, acting as the central scheduling platform for the entire architecture. It is primarily responsible for data interaction and timing synchronization between hydrodynamic federation member 102 and power federation member 103, as well as configuring their publish / subscribe relationships. It should be noted that hydrodynamic federation member 102 and power federation member 103 cannot communicate directly. Hydrodynamic federation member 102 is a federation member registered by the hydrodynamic simulation program. The hydrodynamic simulation program is an independent simulation software / program module specifically designed to solve two-dimensional shallow water dynamic equations and simulate the evolution of regional flooding and waterlogging. In the HLA architecture, it runs as an independent federation member and can be deployed on a terminal or server. Power Federation Member 103 refers to a federation member registered by the power grid flow calculation program. The power grid flow calculation program is an independent simulation software / program module that specifically receives the hydrodynamic parameter matrix forwarded by RTI and completes the entire process of equipment failure judgment, topology update, flow calculation, and protection verification. It runs as another independent federation member in the HLA architecture and can be deployed on a terminal or server.

[0067] In one exemplary embodiment, such as Figure 2 As shown, a hydrodynamic-grid power flow coupling simulation method based on multi-timescale and HLA architecture is presented, which is then applied to... Figure 1 The following steps are used as an example of the power federation members in the HLA architecture, including steps S201 to S205. Wherein:

[0068] Step S201: Receive the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time forwarded by RTI; Hydrodynamic Federation members are used to perform simulation calculations of the flood inundation process of the region to be analyzed according to the preset hydrodynamic calculation step size, obtain the hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to RTI.

[0069] The HLA architecture also includes hydrodynamic federation members and RTIs. The power federation members are registered by the power grid flow calculation program, and the hydrodynamic federation members are registered by the hydrodynamic simulation program.

[0070] The preset hydrodynamic calculation step size refers to the pre-set hydrodynamic calculation step size, specifically the interval between two adjacent simulation moments for hydrodynamics, usually on the order of minutes, such as 5 minutes. The preset hydrodynamic calculation step size is larger than the preset tidal current calculation step size.

[0071] The area to be analyzed refers to the area where power grid flood risk simulation is required, such as a city, an administrative region, or a power distribution network area.

[0072] The flooding process in the area to be analyzed can be simulated using two-dimensional shallow water dynamic equations.

[0073] Each hydrodynamic parameter matrix carries a corresponding simulation time, such as the 5th minute, the 10th minute, the 15th minute, etc. The hydrodynamic parameter matrix includes the hydrodynamic parameters corresponding to each hydrodynamic calculation grid of the region to be analyzed; the hydrodynamic parameters include water depth, flow velocity, and duration of inundation (i.e., inundation duration).

[0074] For example, the hydrodynamic federation member performs simulation calculations of the flood inundation process in the area to be analyzed according to a preset hydrodynamic calculation step size and based on two-dimensional shallow water dynamics equations. This yields the hydrodynamic parameters of each hydrodynamic calculation grid in the area to be analyzed at the corresponding simulation time. Then, the hydrodynamic parameters of each hydrodynamic calculation grid in the area to be analyzed at the corresponding simulation time are summarized to obtain the grid numerical matrix of the area to be analyzed at the corresponding simulation time. This matrix serves as the hydrodynamic parameter matrix carrying the corresponding simulation time, such as at minute 5 and minute 10, and is published to the RTI. Upon receiving the hydrodynamic parameter matrix carrying the corresponding simulation time, the RTI sends the hydrodynamic parameter matrix for the area to be analyzed at the current simulation time (e.g., minute 5) to the power federation member.

[0075] Step S202: Based on the hydrodynamic parameter matrix, determine the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed, determine the equipment status of the power grid equipment based on the failure probability, update the power grid topology based on the equipment status, and obtain the updated power grid topology.

[0076] Among them, the failure probability of power grid equipment is used to measure whether the power grid equipment has failed. For example, if the failure probability of power grid equipment is greater than the preset probability, the power grid equipment is confirmed to have failed.

[0077] The power grid topology includes multiple power grid devices, such as substations, ring main units, switch stations, and low-voltage distribution boxes.

[0078] Among them, based on the hydrodynamic parameter matrix, the hydrodynamic parameters of the location of the power grid equipment in the power grid topology corresponding to the region to be analyzed can be determined, and the failure probability of the power grid equipment can be determined based on the hydrodynamic parameters of the location of the power grid equipment.

[0079] The equipment status of power grid equipment includes either a failure status or a normal status.

[0080] Updating the power grid topology refers to setting power grid equipment in a failed state to be out of operation and disconnecting the corresponding electrical connections of the power grid equipment.

[0081] For example, members of the power federation parse the hydrodynamic parameter matrix according to the matrix parsing instructions to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed. Then, according to the reference coordinate system corresponding to the power grid topology, the coordinate information of each hydrodynamic calculation grid is transformed to obtain the transformed coordinate information of each hydrodynamic calculation grid. Next, based on the transformed coordinate information of each hydrodynamic calculation grid, the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology is determined from each hydrodynamic calculation grid. Then, based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, the hydrodynamic parameters of the location of the power grid equipment are determined. Finally, based on the equipment type of the power grid equipment, the hydrodynamic parameter threshold and vulnerability parameter corresponding to that equipment type are determined as the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment. The hydrodynamic parameters, hydrodynamic parameter thresholds, and vulnerability parameters are then input into a pre-constructed failure probability statistical model to obtain the failure probability of the power grid equipment.

[0082] Furthermore, based on the relationship between the failure probability and the preset probability, the members of the power federation determine the equipment status of the power grid equipment and update the power grid topology based on the equipment status. For example, they set the power grid equipment with the failure status in the power grid topology to be out of operation and disconnect the electrical connection relationship corresponding to the power grid equipment to obtain the updated power grid topology.

[0083] Step S203: Obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain the protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0084] Among them, the power flow calculation results refer to the electrical operating parameters obtained after solving the power flow equations based on the updated power grid topology, including the voltage amplitude of each distribution network node, the active and reactive power of each line / branch, and branch current, etc. Based on the results, overcurrent, voltage over-limit protection verification can be carried out to determine whether the power grid topology needs to be updated again.

[0085] Among them, protection verification refers to verifying whether the distribution network lines and equipment have experienced faults or abnormalities such as overcurrent or voltage exceeding limits based on the power grid flow calculation results, and simulating the tripping judgment logic of relay protection.

[0086] Among them, the protection verification result is the judgment conclusion after verification; if the protection verification result is that there is an electrical over-limit fault, the relay protection will disconnect the corresponding electrical branch, and the power grid topology needs to be updated again and the power flow recalculated; if the protection verification result is that there is no electrical abnormality, the power grid operating state at the current sub-update simulation time remains unchanged.

[0087] The statement "until the current simulation time no new device status changes" means that no new device status changes within the current power flow calculation step. For example, if the current simulation time is 5 minutes and 10 seconds and no new device status changes, then the sub-update simulation time is determined to be 5 minutes and 20 seconds. The power flow calculation and protection verification process is repeated until the next simulation time is reached, such as the 10th minute.

[0088] For example, power federation members calculate the power flow calculation results of the grid devices in the updated grid topology based on the grid flow calculation equations or grid flow calculation instructions. Then, based on the grid flow calculation results, they perform protection checks on the grid devices in the grid topology to determine whether there are operating conditions such as branch overload, node voltage exceeding limits, or islanding, and obtain the corresponding protection check results. Next, based on the protection check results, if it is confirmed that the grid devices in the updated grid topology meet the protection action or operation constraints, the states of the corresponding switches, lines, or nodes in the grid topology are updated to update the updated grid topology again, and the above power flow calculation and protection check steps are executed cyclically until no new device states change at the current simulation time.

[0089] Step S204: Determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0090] The preset power flow calculation step size refers to the pre-set power flow calculation step size, which specifically refers to the interval between two adjacent sub-update simulation times for power flow calculation. It is usually on the order of seconds, such as 10 seconds.

[0091] Wherein, the sub-update simulation time = current simulation time + preset power flow calculation step size; for example, if the current simulation time is 5 minutes and the preset power flow calculation step size is 10 seconds, then the sub-update simulation time is 5 minutes and 10 seconds.

[0092] It should be noted that within each preset power flow calculation step between the current simulation time and the next simulation time, the power federation members cyclically perform power flow calculation and protection verification until no new equipment status changes within the current power flow calculation step. Then, they jump to the next power flow calculation step and continue with equipment failure judgment, topology update, power flow recalculation, and protection actions. This process is repeated continuously until the next simulation time is reached. When the next simulation time is reached, if the hydrodynamic parameter matrix has been updated, step S201 is executed again until the simulation termination condition is met, at which point the simulation stops.

[0093] For example, the power federation member adds the preset power flow calculation step size to the current simulation time to obtain the updated simulation time, which is used as the sub-updated simulation time of the current simulation time. Then, it jumps to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, that is, it jumps to step S202 and continues to repeat the above steps until the next simulation time of the current simulation time is reached.

[0094] Step S205: Take the next simulation time as the new current simulation time, and jump to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

[0095] The simulation termination condition can refer to reaching a preset simulation duration, such as 6 hours.

[0096] For example, the power federation member takes the next simulation moment as the new current simulation moment and jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation moment forwarded by RTI, that is, jumps to step S201 and continues to repeat the above steps until the simulation end condition is met, such as reaching the preset simulation duration, then the simulation stops.

[0097] For example, this application also provides another hydrodynamic-power flow coupling simulation method based on multi-timescale and HLA architecture, specifically including the following steps: Step 1: Establish a collaborative simulation federation conforming to the High-Level Architecture (HLA), load the Federation Object Model (FOM), register the hydrodynamic simulation program as a member of the hydrodynamic federation, register the power flow calculation program as a member of the power federation, and configure the publish / subscribe relationship between the two in the Runtime Structure (RTI); Step 2: Set up a multi-timescale time advancement mechanism, setting the hydrodynamic calculation step size for the hydrodynamic federation member and the power flow calculation step size for the power federation member respectively, and authorizing the time advancement by the RTI based on timestamp interaction messages, time advancement requests, and preset lookahead values; Step 3: The hydrodynamic federation member calculates the flood inundation process, generates a timestamped hydrodynamic parameter matrix, and publishes it to the RTI; the power federation member subscribes to the hydrodynamic parameter matrix and, based on the power grid GIS (Geographic Information System),... The Geographic Information System (GIS) topology maps water depth, flow velocity, and inundation duration to the corresponding power grid equipment nodes, and calculates the failure probability and equipment status of the equipment nodes at the current simulation time. Step 4: When the status of equipment nodes or lines changes, the power federation members update the power grid topology and perform AC power flow calculations to determine branch overload, voltage overrun, and islanding. If the protection action or operation constraints are met, the status of the corresponding switches, lines, or nodes is updated, and the power flow check is repeated until there are no new status changes in the current time step.

[0098] The federated object model includes hydrodynamic frame objects, power grid equipment objects, and topology event objects. Hydrodynamic frame objects include at least the fields of simulation time, calculation grid number, water depth, x-direction flow velocity, y-direction flow velocity, and inundation duration. Power grid equipment objects include at least the fields of equipment number, equipment type, bus number, geographic coordinates, and equipment operating status. Topology event objects include at least the fields of event time, event type, associated equipment number, and post-action status.

[0099] The multi-timescale time advancement mechanism includes a conservative synchronization strategy. Under the conservative synchronization strategy, RTI determines the safe time boundary for advancement based on the next request time and lookahead value submitted by each federation member, and returns time advancement authorization to the corresponding federation member to avoid power federation members performing power flow calculations in advance before receiving the corresponding hydrodynamic frame.

[0100] Among them, members of the hydrodynamic federation calculated the flood inundation process based on two-dimensional shallow water dynamic equations, which include a continuity equation and momentum equations in two directions, as detailed below:

[0101] The continuity equation is:

[0102] ;

[0103] The momentum equation in the x-direction is:

[0104] ;

[0105] The momentum equation in the y-direction is:

[0106] ;

[0107] Where h is the water depth, u and v are the average vertical flow velocities in the x and y directions, respectively, q is the source-sink term, g is the gravitational acceleration, Z0 is the bottom elevation, and n is the Manning roughness coefficient.

[0108] The process by which members of the power federation map the hydrodynamic parameter matrix to the power grid equipment nodes includes: converting the hydrodynamic calculation grid to the power grid GIS topology layer according to the coordinate reference system; performing grid sampling or surface intersection calculation on the location points of substations, distribution cabinets, switch stations and towers; and when equipment spans multiple hydrodynamic calculation grids, calculating the corresponding water depth, flow velocity and inundation duration according to the equipment's occupied area, equipment center point or preset key points.

[0109] The AC power flow calculation adopts the Newton-Raphson method. When the equipment status changes and causes a topology change, the node admittance matrix is ​​reconstructed, the voltage amplitude and phase angle of each node are solved iteratively, and the line and switch status are updated according to the branch thermal stability limit, the upper and lower limits of the node voltage and the relay protection action delay.

[0110] In the aforementioned hydrodynamic-grid power flow coupling simulation method based on multi-timescale and HLA architecture, the following steps are first taken: The hydrodynamic parameter matrix for the region to be analyzed at the current simulation time is received from the RTI. Hydrodynamic federation members then perform simulation calculations of the flood inundation process in the region to be analyzed according to a preset hydrodynamic calculation step size, obtaining a hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI. Next, based on the hydrodynamic parameter matrix, the failure probability of grid equipment in the corresponding grid topology of the region to be analyzed is determined. Based on the failure probability, the equipment status of the grid equipment is determined, and the grid topology is updated based on the equipment status to obtain the updated grid topology. Finally, the electrical parameters of the grid equipment in the updated grid topology are obtained. The power flow calculation results are used to perform protection verification, resulting in protection verification results. Based on these results, the updated power grid topology is updated again until no new equipment states change at the current simulation time. Next, the sub-update simulation time is determined based on the preset power flow calculation step size, and the process jumps to the step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time is reached. Finally, the next simulation time is taken as the new current simulation time, and the process jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time forwarded by RTI, until the simulation termination condition is met. In this way, under the HLA architecture, the hydrodynamic simulation program and the power flow calculation program are coupled and simulated based on RTI. The hydrodynamic parameter matrix for the flood inundation process of the area to be analyzed is updated by pre-setting the hydrodynamic calculation step size. At the same time, the sub-update simulation time is obtained by using the pre-set power flow calculation step size within the adjacent simulation time for hydrodynamics. The entire process of power grid equipment failure, topology update and protection verification is continuously iterated and deduced. It can completely reproduce the dynamic evolution details of the power grid cascading faults between two adjacent simulation time points, which makes up for the shortcomings of traditional offline serial simulation, which can only calculate the power grid state at a single point and thus has poor simulation integrity. At the same time, the transmission and reception of the hydrodynamic parameter matrix is ​​realized through RTI, which solves the problem of time series correlation distortion caused by the mismatch of multi-model data structure and spatiotemporal scale. It accurately matches the hydrodynamic parameters corresponding to the power grid equipment and quantifies the equipment failure probability, effectively avoiding the problem of distortion of traditional simulation results. This makes the simulation results of the power grid fault evolution process under flood scenarios more in line with the actual working conditions, thereby improving the accuracy of power grid flood risk simulation assessment.

[0111] In one exemplary embodiment, such as Figure 3 As shown, step S202 above, which determines the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, specifically includes the following steps S301 to S305. Wherein:

[0112] Step S301: Analyze the hydrodynamic parameter matrix to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed.

[0113] Step S302: According to the reference coordinate system corresponding to the power grid topology, the coordinate information of each hydrodynamic calculation grid is transformed to obtain the transformed coordinate information of each hydrodynamic calculation grid.

[0114] Step S303: Based on the transformed coordinate information of each hydrodynamic calculation grid, determine the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology.

[0115] Step S304: Based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, determine the hydrodynamic parameters of the location of the power grid equipment.

[0116] Step S305: Determine the failure probability of the power grid equipment based on the hydrodynamic parameters of the location of the power grid equipment.

[0117] The target hydrodynamic calculation grid refers to the hydrodynamic calculation grid corresponding to the location of the power grid equipment. The hydrodynamic parameters at the location of the power grid equipment refer to the water depth, flow velocity, and duration of submersion at that location.

[0118] For example, members of the power federation obtain matrix parsing instructions and parse the hydrodynamic parameter matrix according to the instructions to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed. Next, they obtain the reference coordinate system corresponding to the power grid topology and perform coordinate transformation processing on the coordinate information of each hydrodynamic calculation grid according to the reference coordinate system corresponding to the power grid topology to obtain the transformed coordinate information of each hydrodynamic calculation grid. Next, based on the transformed coordinate information of each hydrodynamic calculation grid, they determine the hydrodynamic calculation grid corresponding to the location of the power grid equipment in the power grid topology from each hydrodynamic calculation grid, as the target hydrodynamic calculation grid corresponding to the power grid equipment. Next, based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, and the number of target hydrodynamic calculation grids, they determine the hydrodynamic parameters of the location of the power grid equipment. Finally, they input the hydrodynamic parameters of the location of the power grid equipment into a pre-built failure probability statistical model to obtain the failure probability of the power grid equipment.

[0119] In this embodiment, when determining the failure probability of power grid equipment, the coordinate information of each hydrodynamic calculation grid is first transformed, and the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology are comprehensively considered, which helps to improve the accuracy of failure probability determination.

[0120] In an exemplary embodiment, step S304, which determines the hydrodynamic parameters of the location of the power grid equipment based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, specifically includes the following: when there is only one target hydrodynamic calculation grid corresponding to the power grid equipment, the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment are used as the hydrodynamic parameters of the location of the power grid equipment; when there are multiple target hydrodynamic calculation grids corresponding to the power grid equipment, the hydrodynamic parameters of each target hydrodynamic calculation grid corresponding to the power grid equipment are fused according to the superposition area of ​​the power grid equipment and each target hydrodynamic calculation grid to obtain the hydrodynamic parameters of the location of the power grid equipment.

[0121] If the power grid equipment is a point-like device, the corresponding target hydrodynamic calculation grid is 1. If the power grid equipment is a surface-like device, the corresponding target hydrodynamic calculation grid is multiple.

[0122] It should be noted that when performing fusion processing on hydrodynamic parameters, the fusion processing is performed separately for each hydrodynamic parameter. For example, the water depth of each target hydrodynamic calculation grid corresponding to the power grid equipment is fused, and the flow velocity of each target hydrodynamic calculation grid corresponding to the power grid equipment is fused.

[0123] For example, the power federation members first determine whether the power grid equipment is a point-like or area-like device. If the power grid equipment is a point-like device, the number of target hydrodynamic calculation grids corresponding to the power grid equipment is confirmed to be one, and the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment are directly used as the hydrodynamic parameters of the location of the power grid equipment. If the power grid equipment is an area-like device, the number of target hydrodynamic calculation grids corresponding to the power grid equipment is confirmed to be multiple, and the superposition area of ​​the power grid equipment and each target hydrodynamic calculation grid is obtained. Based on the superposition area of ​​the power grid equipment and each target hydrodynamic calculation grid, the weights corresponding to the hydrodynamic parameters of each target hydrodynamic calculation grid are confirmed. According to the weights corresponding to the hydrodynamic parameters of each target hydrodynamic calculation grid, the hydrodynamic parameters of each target hydrodynamic calculation grid are fused to obtain the hydrodynamic parameters of the location of the power grid equipment. Of course, in the case of the power grid equipment being an area-like device, the power federation members can also directly use the hydrodynamic parameters of the target hydrodynamic calculation grids corresponding to the center point or preset key point of the power grid equipment as the hydrodynamic parameters of the location of the power grid equipment.

[0124] For example, suppose the area to be analyzed is divided into a 10m×10m hydrodynamic calculation grid, with three adjacent grids A, B, and C. If the power grid equipment is a point device (single grid), such as a tower whose coordinates fall within grid A, then grid sampling is performed, and the hydrodynamic parameters of grid A are directly used as the hydrodynamic parameters corresponding to the location of the power grid equipment. If the power grid equipment is a planar device (spanning multiple grids), such as a medium-sized substation covering both grids A and B, which is considered a grid-spanning device, there are three ways to determine the hydrodynamic parameters corresponding to the location of the power grid equipment: Method 1 (Area): 70% of the substation is in grid A and 30% is in grid B. The weighted average of these values ​​is used to obtain the comprehensive hydrodynamic parameters, which are then used as the hydrodynamic parameters corresponding to the location of the power grid equipment. Method 2 (Center Point): The center point of the substation falls in grid B. The hydrodynamic parameters of grid B are directly used as the hydrodynamic parameters corresponding to the location of the power grid equipment. Method 3 (Key Point): The low-lying distribution room of the substation is a preset key point. This point is in grid A. The hydrodynamic parameters of grid A are used as the hydrodynamic parameters corresponding to the location of the power grid equipment.

[0125] In this embodiment, when determining the hydrodynamic parameters of the location of the power grid equipment, the number of target hydrodynamic calculation grids corresponding to the power grid equipment and the hydrodynamic parameters of the target hydrodynamic calculation grids corresponding to the power grid equipment are taken into account, which helps to improve the accuracy of determining the hydrodynamic parameters of the location of the power grid equipment.

[0126] In an exemplary embodiment, step S305 above, which determines the failure probability of the power grid equipment based on the hydrodynamic parameters of the location of the power grid equipment, specifically includes the following: determining the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment based on the equipment type of the power grid equipment; the vulnerability parameter is used to characterize the influence weight of the hydrodynamic parameters on the failure probability; and inputting the hydrodynamic parameters, hydrodynamic parameter threshold, and vulnerability parameter into a pre-built failure probability statistical model to obtain the failure probability of the power grid equipment.

[0127] Different equipment types correspond to different hydrodynamic parameter thresholds and vulnerability parameters. Hydrodynamic parameter thresholds can refer to water depth thresholds, flow velocity thresholds, and duration of inundation thresholds. Vulnerability parameters can refer to vulnerability parameters corresponding to water depth (such as the first vulnerability parameter), vulnerability parameters corresponding to flow velocity (such as the second vulnerability parameter), and vulnerability parameters corresponding to duration of inundation (such as the third vulnerability parameter).

[0128] Among them, the vulnerability parameter refers to the parameter used to measure the weight of the influence of hydrodynamic parameters on the failure probability.

[0129] Among them, the failure probability statistical model refers to the statistical model constructed based on the time-varying vulnerability function.

[0130] The failure probability of the power grid equipment is 1 / {1+exp[-(first vulnerability parameter (water depth - water depth threshold) + second vulnerability parameter (flow velocity - flow velocity threshold) + third vulnerability parameter (continuous flooding time - continuous flooding time threshold)]}.

[0131] For example, based on the equipment type of the power grid equipment, members of the power federation query the correspondence between equipment type and hydrodynamic parameter thresholds, as well as the correspondence between equipment type and vulnerability parameters, to obtain the hydrodynamic parameter thresholds and vulnerability parameters corresponding to the location of the power grid equipment. Then, the hydrodynamic parameters, hydrodynamic parameter thresholds, and vulnerability parameters are input into a pre-built failure probability statistical model. Through the failure probability statistical model, based on the hydrodynamic parameters, hydrodynamic parameter thresholds, and vulnerability parameters, failure probability statistical processing is performed to obtain the failure probability of the power grid equipment.

[0132] For example, the failure probability is calculated using a time-varying vulnerability function with different parameters according to the equipment type. The time-varying vulnerability function is as follows:

[0133] ;

[0134] Where Pf(i,t) is the failure probability of device i at simulation time t, h_i is the water depth at the device location, v_i is the synthetic flow velocity at the device location, τ_i is the duration of submersion, h0_i, v0_i and τ0_i are the water depth threshold, flow velocity threshold and duration of submersion threshold for the corresponding type of device i, respectively, and a_i, b_i and c_i are pre-calibrated vulnerability parameters.

[0135] In this embodiment, when obtaining the failure probability of power grid equipment, hydrodynamic parameters, hydrodynamic parameter thresholds, and vulnerability parameters are comprehensively considered, which makes the calculated failure probability more accurate and further improves the accuracy of failure probability determination.

[0136] In an exemplary embodiment, based on the equipment type of the power grid equipment, the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters at the location of the power grid equipment are determined. Specifically, this includes: querying the correspondence between the equipment type and the hydrodynamic parameter threshold based on the equipment type of the power grid equipment to obtain the hydrodynamic parameter threshold corresponding to the equipment type, which is used as the hydrodynamic parameter threshold corresponding to the location of the power grid equipment; and querying the correspondence between the equipment type and the vulnerability parameter based on the equipment type of the power grid equipment to obtain the vulnerability parameter corresponding to the equipment type, which is used as the vulnerability parameter corresponding to the hydrodynamic parameters at the location of the power grid equipment.

[0137] The correspondence between equipment type and hydrodynamic parameter thresholds is used to indicate that different equipment types correspond to different hydrodynamic parameter thresholds, such as water depth threshold, flow velocity threshold, and duration of submersion threshold.

[0138] The correspondence between device type and vulnerability parameter is used to indicate that different device types correspond to different vulnerability parameters, such as the first vulnerability parameter, the second vulnerability parameter, and the third vulnerability parameter.

[0139] For example, the power federation members first obtain the correspondence between equipment type and hydrodynamic parameter thresholds. Then, based on the equipment type of the power grid equipment, they query the correspondence between equipment type and hydrodynamic parameter thresholds to obtain the hydrodynamic parameter threshold corresponding to the equipment type, and use it as the hydrodynamic parameter threshold corresponding to the location of the power grid equipment. Finally, they obtain the correspondence between equipment type and vulnerability parameters, and based on the equipment type of the power grid equipment, they query the correspondence between equipment type and vulnerability parameters to obtain the vulnerability parameter corresponding to the equipment type, and use it as the vulnerability parameter corresponding to the hydrodynamic parameter corresponding to the location of the power grid equipment.

[0140] In this embodiment, by querying the correspondence between equipment type and hydrodynamic parameter threshold, as well as the correspondence between equipment type and vulnerability parameter, it is beneficial to quickly determine the hydrodynamic parameter threshold corresponding to the location of the power grid equipment, and the vulnerability parameter corresponding to the hydrodynamic parameter at the location of the power grid equipment, thereby improving the efficiency of determining the hydrodynamic parameter threshold and vulnerability parameter.

[0141] In an exemplary embodiment, step S202 above, which determines the equipment status of the power grid equipment based on the failure probability, specifically includes the following: if the failure probability is greater than a preset probability, the equipment status of the power grid equipment is determined to be a failed state; if the failure probability is less than or equal to the preset probability, the equipment status of the power grid equipment is determined to be a normal state.

[0142] Among them, the equipment status of power grid equipment is the failure status, which is used to indicate that the power grid equipment needs to be taken out of operation.

[0143] For example, members of the power federation compare the failure probability with a preset probability. If the failure probability is greater than the preset probability, the power grid equipment is determined to be in a failed state and needs to be taken out of operation. If the failure probability is less than or equal to the preset probability, the power grid equipment is determined to be in a normal state.

[0144] In this embodiment, the comparison between the failure probability and the preset probability helps to accurately determine the equipment status of the power grid equipment, thereby improving the accuracy of equipment status determination.

[0145] In one exemplary embodiment, this application also provides a hydrodynamic-power grid flow coupling simulation system based on multi-timescale and HLA architecture, specifically including:

[0146] The Federation Management module is used to load the Federation object model, initialize RTI connections, and register hydrodynamic and power federation members.

[0147] The time management module is used to configure the hydrodynamic calculation step size, tidal current calculation step size, and lookahead value, and to handle time advance requests, time advance authorizations, and timestamped interactive messages.

[0148] The hydrodynamic calculation module is used to solve the two-dimensional shallow water dynamic equations and generate a hydrodynamic parameter matrix that includes water depth, flow velocity, and inundation duration.

[0149] The hydrodynamic-grid mapping module is used to spatially match the hydrodynamic parameter matrix with the power grid GIS topology to form the flooding exposure status of the power grid equipment nodes.

[0150] The equipment status assessment module is used to calculate the probability of equipment failure and update the equipment status based on equipment type, exposure status, and time-varying vulnerability function.

[0151] The power flow and protection simulation module is used to perform AC power flow calculations after a change in the power grid topology, and to update the power grid operating status based on overload, voltage limit exceedance, islanding, and protection action rules.

[0152] The aforementioned hydrodynamic-power flow coupling simulation system based on multi-timescale and HLA architecture can describe the mutual influence between flood evolution, equipment failure and power flow transfer on the same simulation time axis. It is suitable for analyzing the power outage range of distribution networks and comparing emergency dispatch schemes under flood disasters.

[0153] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages in other steps. It is understood that the steps in different embodiments can be freely combined as needed, and all non-contradictory solutions formed by such combinations are within the scope of protection of this application.

[0154] Based on the same inventive concept, this application also provides a hydrodynamic-power grid flow coupling simulation device based on a multi-timescale and HLA architecture for implementing the above-mentioned hydrodynamic-power grid flow coupling simulation method based on a multi-timescale and HLA architecture. The solution provided by this device is similar to the solution described in the above method. Therefore, the specific limitations of one or more embodiments of the hydrodynamic-power grid flow coupling simulation device based on a multi-timescale and HLA architecture provided below can be found in the limitations of the hydrodynamic-power grid flow coupling simulation method based on a multi-timescale and HLA architecture described above, and will not be repeated here.

[0155] In one exemplary embodiment, such as Figure 4 As shown, a hydrodynamic-grid power flow coupling simulation device 400 based on a multi-timescale and HLA architecture is provided, applied to the power federation members in the HLA architecture; the HLA architecture also includes hydrodynamic federation members and RTIs, the power federation members are registered by the grid power flow calculation program, and the hydrodynamic federation members are registered by the hydrodynamic simulation program; it includes: a matrix receiving module 401, a topology update module 402, a first processing module 403, a second processing module 404, and a third processing module 405, wherein:

[0156] The matrix receiving module 401 is used to receive the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time, which is forwarded by the RTI; the hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the region to be analyzed according to the preset hydrodynamic calculation step size, obtain the hydrodynamic parameter matrix carrying the corresponding simulation time, and publish the hydrodynamic parameter matrix to the RTI.

[0157] The topology update module 402 is used to determine the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, update the power grid topology based on the equipment status, and obtain the updated power grid topology.

[0158] The first processing module 403 is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain the protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time.

[0159] The second processing module 404 is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached.

[0160] The third processing module 405 is used to take the next simulation moment as the new current simulation moment and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation moment forwarded by RTI, until the simulation end condition is met.

[0161] In an exemplary embodiment, the topology update module 402 is further configured to parse the hydrodynamic parameter matrix to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed; transform the coordinate information of each hydrodynamic calculation grid according to the reference coordinate system corresponding to the power grid topology to obtain the transformed coordinate information of each hydrodynamic calculation grid; determine the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology based on the transformed coordinate information of each hydrodynamic calculation grid; determine the hydrodynamic parameters of the location of the power grid equipment based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology; and determine the failure probability of the power grid equipment based on the hydrodynamic parameters of the location of the power grid equipment.

[0162] In an exemplary embodiment, the topology update module 402 is further configured to, when there is only one target hydrodynamic calculation grid corresponding to the power grid device, use the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid device as the hydrodynamic parameters of the location of the power grid device; when there are multiple target hydrodynamic calculation grids corresponding to the power grid device, perform fusion processing on the hydrodynamic parameters of each target hydrodynamic calculation grid corresponding to the power grid device according to the superposition area of ​​the power grid device and each target hydrodynamic calculation grid to obtain the hydrodynamic parameters of the location of the power grid device.

[0163] In an exemplary embodiment, the topology update module 402 is further configured to determine the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment based on the equipment type of the power grid equipment; the vulnerability parameter is used to characterize the influence weight of the hydrodynamic parameters on the failure probability; and the hydrodynamic parameters, hydrodynamic parameter threshold and vulnerability parameter are input into a pre-built failure probability statistical model to obtain the failure probability of the power grid equipment.

[0164] In an exemplary embodiment, the topology update module 402 is further configured to query the correspondence between the equipment type and the hydrodynamic parameter threshold based on the equipment type of the power grid equipment, and obtain the hydrodynamic parameter threshold corresponding to the equipment type as the hydrodynamic parameter threshold corresponding to the location of the power grid equipment; and query the correspondence between the equipment type and the vulnerability parameter based on the equipment type of the power grid equipment, and obtain the vulnerability parameter corresponding to the equipment type as the vulnerability parameter corresponding to the location of the power grid equipment.

[0165] In an exemplary embodiment, the topology update module 402 is further configured to determine the equipment status of the power grid equipment as a failed state when the failure probability is greater than a preset probability; and to determine the equipment status of the power grid equipment as a normal state when the failure probability is less than or equal to the preset probability.

[0166] The modules in the aforementioned hydrodynamic-power grid flow coupling simulation device based on multi-timescale and HLA architecture can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the computer device's memory as software, so that the processor can call and execute the corresponding operations of each module.

[0167] In one exemplary embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 5 As shown, this computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operating system and computer programs stored in the non-volatile storage media. The database stores data such as hydrodynamic parameter matrices. The I / O interfaces are used for exchanging information between the processor and external devices. The communication interface is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements a hydrodynamic-power grid flow coupling simulation method based on a multi-timescale and HLA architecture.

[0168] Those skilled in the art will understand that Figure 5 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0169] In one exemplary embodiment, a computer device is also provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0170] In one exemplary embodiment, a computer-readable storage medium is provided having a computer program stored thereon that, when executed by a processor, implements the steps in the above-described method embodiments.

[0171] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0172] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, artificial intelligence (AI) processors, etc., and are not limited to these.

[0173] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0174] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A hydrodynamic-power grid flow coupling simulation method based on multi-timescale and HLA architecture, characterized in that, The method applies to power federation members in an HLA architecture; the HLA architecture also includes hydrodynamic federation members and RTIs, wherein the power federation members are registered by a power grid flow calculation program, and the hydrodynamic federation members are registered by a hydrodynamic simulation program; the method includes: The system receives the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time, forwarded by the RTI. The hydrodynamic federation member is used to perform simulation calculations on the flood inundation process of the region to be analyzed according to a preset hydrodynamic calculation step size, to obtain the hydrodynamic parameter matrix carrying the corresponding simulation time, and publishes the hydrodynamic parameter matrix to the RTI. Based on the hydrodynamic parameter matrix, the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed is determined. Based on the failure probability, the equipment status of the power grid equipment is determined. Based on the equipment status, the power grid topology is updated to obtain the updated power grid topology. Obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain the protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time; The sub-update simulation time of the current simulation time is determined based on the preset power flow calculation step size, and then the process jumps to the step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached. The next simulation moment is taken as the new current simulation moment, and the process jumps to the step of receiving the hydrodynamic parameter matrix for the region to be analyzed at the current simulation moment forwarded by RTI, until the simulation end condition is met.

2. The method according to claim 1, characterized in that, The step of determining the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix includes: The hydrodynamic parameter matrix is ​​analyzed to obtain the coordinate information and hydrodynamic parameters of each hydrodynamic calculation grid corresponding to the region to be analyzed; According to the reference coordinate system corresponding to the power grid topology, the coordinate information of each hydrodynamic calculation grid is transformed to obtain the transformed coordinate information of each hydrodynamic calculation grid. Based on the transformed coordinate information of each hydrodynamic calculation grid, the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology is determined; Based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology, the hydrodynamic parameters of the location of the power grid equipment are determined. The failure probability of the power grid equipment is determined based on the hydrodynamic parameters of its location.

3. The method according to claim 2, characterized in that, The process of determining the hydrodynamic parameters of the location of the power grid equipment based on the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment in the power grid topology includes: When the number of target hydrodynamic calculation grids corresponding to the power grid equipment is one, the hydrodynamic parameters of the target hydrodynamic calculation grid corresponding to the power grid equipment are used as the hydrodynamic parameters of the location of the power grid equipment. When there are multiple target hydrodynamic calculation grids corresponding to the power grid equipment, the hydrodynamic parameters of each target hydrodynamic calculation grid corresponding to the power grid equipment are fused according to the superposition area of ​​the power grid equipment and each target hydrodynamic calculation grid to obtain the hydrodynamic parameters of the location of the power grid equipment.

4. The method according to claim 2, characterized in that, Determining the failure probability of the power grid equipment based on the hydrodynamic parameters of its location includes: Based on the equipment type of the power grid equipment, the hydrodynamic parameter threshold and vulnerability parameter corresponding to the hydrodynamic parameters of the location of the power grid equipment are determined; the vulnerability parameter is used to characterize the influence weight of the hydrodynamic parameters on the failure probability. The failure probability of the power grid equipment is obtained by inputting the hydrodynamic parameters, the hydrodynamic parameter thresholds, and the vulnerability parameters into a pre-built failure probability statistical model.

5. The method according to claim 4, characterized in that, The step of determining the hydrodynamic parameter thresholds and vulnerability parameters corresponding to the hydrodynamic parameters at the location of the power grid equipment based on the equipment type of the power grid equipment includes: Based on the equipment type of the power grid equipment, the correspondence between the equipment type and the hydrodynamic parameter threshold is queried to obtain the hydrodynamic parameter threshold corresponding to the equipment type, which is then used as the hydrodynamic parameter threshold corresponding to the location of the power grid equipment. Based on the equipment type of the power grid equipment, the correspondence between the equipment type and the vulnerability parameter is queried to obtain the vulnerability parameter corresponding to the equipment type, which is then used as the vulnerability parameter corresponding to the hydrodynamic parameter of the location of the power grid equipment.

6. The method according to any one of claims 1 to 5, characterized in that, Determining the equipment status of the power grid equipment based on the failure probability includes: If the failure probability is greater than a preset probability, the power grid equipment is determined to be in a failure state. If the failure probability is less than or equal to the preset probability, the power grid equipment is determined to be in a normal state.

7. A hydrodynamic-electric power flow coupling simulation device based on multi-timescale and HLA architecture, characterized in that, The device includes: power federation members applied in an HLA architecture; the HLA architecture also includes hydrodynamic federation members and RTIs, wherein the power federation members are registered by a power grid flow calculation program, and the hydrodynamic federation members are registered by a hydrodynamic simulation program; the device includes: The matrix receiving module is used to receive the hydrodynamic parameter matrix for the region to be analyzed at the current simulation time, which is forwarded by the RTI. The hydrodynamic federation members are used to perform simulation calculations on the flood inundation process of the region to be analyzed according to a preset hydrodynamic calculation step size, to obtain the hydrodynamic parameter matrix carrying the corresponding simulation time, and to publish the hydrodynamic parameter matrix to the RTI. The topology update module is used to determine the failure probability of power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, determine the equipment status of the power grid equipment based on the failure probability, and update the power grid topology based on the equipment status to obtain the updated power grid topology. The first processing module is used to obtain the power flow calculation results of the power grid equipment in the updated power grid topology, perform protection verification based on the power flow calculation results, obtain protection verification results, and update the updated power grid topology again based on the protection verification results until no new equipment status changes at the current simulation time. The second processing module is used to determine the sub-update simulation time of the current simulation time based on the preset power flow calculation step size, and jump to the step of determining the failure probability of the power grid equipment in the power grid topology corresponding to the region to be analyzed based on the hydrodynamic parameter matrix, until the next simulation time of the current simulation time is reached. The third processing module is used to take the next simulation time as the new current simulation time and jump to the step of receiving the hydrodynamic parameter matrix of the region to be analyzed at the current simulation time forwarded by RTI, until the simulation end condition is met.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

Images