Power grid operation state simulation method, device, equipment, medium and program product

Abstract

The application discloses a power grid operation state simulation method, device, equipment, medium and program product. The method obtains preset operation state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path; obtains target first time sequence data corresponding to the preset operation state information, and performs phase-mode conversion on the target first time sequence data to obtain first modulus data; based on the unit transmission parameter information, second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated; and the second modulus data is subjected to mode-phase conversion to obtain target second time sequence data. According to the embodiment of the application, the simulation result can be refined, more accurate power distribution network state simulation is realized, the power distribution network electromagnetic state simulation calculation in a long time period or under a plurality of state superpositions can be applied, and strong support can be provided for power distribution network state monitoring and fault diagnosis.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of electronic information electromagnetic calculation, and particularly relates to a simulation method, device, equipment, medium and program product of power grid operation state. BACKGROUND

[0002] The operation state monitoring of the distribution network is an important means to improve the reliability of power supply of users, and can effectively avoid the power loss and the increase of operation and maintenance costs caused by long-time fault finding.

[0003] The existing simulation of the operation state of the power system depends on industrial software, and the calculation of the existing industrial software is mainly for the whole global expansion, cannot customize the distribution parameters of a single transmission grid, the simulation result is insufficient in detail information, and cannot meet the demand of in-depth analysis of the operation state of the distribution network; and cannot be calculated in a unified framework, so the electromagnetic state simulation calculation of the distribution network in a long period of time and in combination of multiple states cannot be realized, and the electromagnetic simulation framework research for the normal-non-normal operation state of the power distribution network is lacking.

[0004] Therefore, when simulating the operation state of the distribution network, the existing technology has the problems of insufficient detail information of the simulation result and inability to realize the electromagnetic state simulation in a long time and in combination of multiple states. SUMMARY

[0005] The embodiments of the application provide a simulation method, device, equipment, medium and program product of the operation state of the power grid, which can refine the simulation result, realize more accurate simulation of the state of the distribution network, and is also applicable to the electromagnetic state simulation calculation of the distribution network in a long period of time or in combination of multiple states, and can provide strong support for the state monitoring and fault diagnosis of the distribution network.

[0006] In a first aspect, the embodiments of the application provide a simulation method of the operation state of the power grid, comprising:

[0007] obtaining preset operation state information of a distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path;

[0008] In the correspondence between the preset operation state and the first time sequence data, target first time sequence data corresponding to the preset operation state information is obtained, and the target first time sequence data is subjected to a phase-mode transformation to obtain first modulus data, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the distribution network;

[0009] Based on the unit transmission parameter information, second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated;

[0010] The second modulus data is subjected to a phase-modulus conversion to obtain target second time sequence data, and the target second time sequence data includes second time sequence data of each phase voltage signal.

[0011] In a possible implementation of the first aspect, the preset operating state information of the power distribution network and the unit transmission parameter information of each unit transmission path in the preset total transmission path are obtained, including:

[0012] The first page is displayed.

[0013] In response to a first input in the first page, the preset operating state information of the power distribution network and preset operating environment configuration information in the first input are obtained.

[0014] According to the preset operating environment configuration information, the unit transmission parameter information of each unit transmission path in the preset total transmission path is calculated.

[0015] In a possible implementation of the first aspect, the preset operating environment configuration information includes terrain information, line distribution information and line material information in the preset total transmission path; and the unit transmission parameter information includes resistance values of distribution resistances, inductance values of distribution inductances, capacitance values of distribution capacitances and conductance values of distribution conductances in the unit transmission path.

[0016] In a possible implementation of the first aspect, the first input further includes first configuration information of the preset total transmission path and second configuration information of the unit transmission path; and the method further includes:

[0017] The first configuration information of the preset total transmission path and the second configuration information of the unit transmission path in the first input are obtained.

[0018] According to the first configuration information, the preset total transmission path is determined; and according to the second configuration information, the unit transmission path is determined.

[0019] In a possible implementation of the first aspect, the target first time sequence data is subjected to a phase-modulus conversion to obtain first modulus data, including:

[0020] The target first time sequence data is multiplied by a preset conversion matrix to obtain the first modulus data.

[0021]

[0022] Q is the conversion matrix.

[0023] In a possible implementation of the first aspect, the second modulus data is subjected to a phase-modulus conversion to obtain target second time sequence data, including:

[0024] The second modulus data is multiplied by a preset inverse conversion matrix to obtain the target second time sequence data.

[0025]

[0026] wherein Q -1 is the inverse transform matrix.

[0027] In a possible implementation of the first aspect, the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated based on the unit transmission parameter information, and includes:

[0028] Obtaining information of a preset terminal reflection coefficient;

[0029] Calculating a terminal matching impedance value according to the information of the terminal reflection coefficient;

[0030] The second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated based on the unit transmission parameter information and the terminal matching impedance value.

[0031] In a possible implementation of the first aspect, the terminal matching impedance value is calculated according to the information of the terminal reflection coefficient, and includes:

[0032] Calculating an equivalent impedance value of the last unit transmission path according to the unit transmission parameter information of the last unit transmission path in the preset total transmission path;

[0033] In a case where the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is taken as the terminal matching impedance value;

[0034]

[0035] wherein Z L is the terminal matching impedance value, Zc is the equivalent impedance value of the last unit transmission path, L, R, C and G are the unit transmission parameter information of the last unit transmission path.

[0036] In a possible implementation of the first aspect, the method further includes:

[0037] In a case where the terminal reflection coefficient is not 0, the terminal matching impedance value is calculated according to the terminal reflection coefficient;

[0038]

[0039] wherein Γ is the terminal reflection coefficient.

[0040] In a possible implementation of the first aspect, the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated based on the unit transmission parameter information and the terminal matching impedance value, and includes:

[0041] Obtaining information of a source internal resistance;

[0042] calculate, based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path, first sub-modulus data formed after the first modulus data is transmitted via the first preset unit transmission path;

[0043] calculate, based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data, N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths, where N is an integer greater than or equal to 1;

[0044] In the case where the N+1th unit transmission path is the last unit transmission path, calculate, based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth sub-modulus data and the terminal matching impedance value, second modulus data formed after the first modulus data is transmitted via each unit transmission path in the preset total transmission path.

[0045] In a possible embodiment of the first aspect, the calculation of the first sub-modulus data formed after the first modulus data is transmitted via the first preset unit transmission path based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path comprises:

[0046]

[0047] where Δx is preset length information of the unit transmission path, Δt is preset time interval information between time points of the first time sequence data, Rs is the information of the source internal resistance, C is the unit transmission parameter information of the first unit transmission path, N represents the serial number of the unit transmission path, k represents the time point, I represents the current, and V represents the voltage.

[0048] In a possible embodiment of the first aspect, the calculation of the N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data comprises:

[0049]

[0050] where L0, R0, C0 and G0 are the unit transmission parameter information of the N+1th unit transmission path, I N and V N are the Nth sub-modulus data, I N+1 and V N+1 are the N+1th sub-modulus data.

[0051] In a possible implementation of the first aspect, in a case where the N+1th unit transmission path is the last unit transmission path, the second modulus data formed after the first modulus data is transmitted via each unit transmission path in the preset total transmission path is calculated based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth sub-modulus data, and the terminal matching impedance value, and includes:

[0052]

[0053] wherein C is the unit transmission parameter information of the last unit transmission path, V N+1 is the second modulus data.

[0054] In a possible implementation of the first aspect, before the phase-to-modulus transformation is performed on the target first time-series data to obtain the first modulus data, the method further includes:

[0055] In a case where the target first time-series data includes voltage data of a mutation type, a mutation time length range in which a time point corresponding to the voltage data of the mutation type is located is determined.

[0056] The voltage data corresponding to the second type of time point in the mutation time length range is obtained by interpolation calculation on the voltage data corresponding to the first type of time point in the mutation time length range as known data.

[0057] The voltage data corresponding to the second type of time point in the target first time-series data is replaced by the corresponding interpolation voltage data to obtain updated target first time-series data.

[0058] The phase-to-modulus transformation is performed on the updated target first time-series data to obtain the first modulus data.

[0059] The phase-to-modulus transformation is performed on the updated target first time-series data to obtain the first modulus data.

[0060] In a possible implementation of the first aspect, the method further includes:

[0061] In the second page, the waveform information corresponding to the target second time-series data is displayed.

[0062] Based on the same inventive concept, in the second aspect, the embodiments of the present application further provide a simulation device for power grid operation state, including:

[0063] The acquisition module is configured to acquire preset operation state information of the power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path.

[0064] The conversion module is configured to, in a preset operating state and a corresponding relationship of the first time sequence data, obtain target first time sequence data corresponding to preset operating state information, and perform phase-to-magnitude conversion on the target first time sequence data to obtain first magnitude data, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network.

[0065] The calculation module is configured to calculate second magnitude data formed after the first magnitude data is transmitted via each unit transmission path in the preset total transmission path based on each unit transmission parameter information.

[0066] The conversion module is further configured to perform magnitude-to-phase conversion on the second magnitude data to obtain target second time sequence data, and the target second time sequence data includes second time sequence data of each phase voltage signal.

[0067] Based on the same inventive concept, in a third aspect, the embodiments of the present application further provide a power grid operating state simulation device, which comprises a processor and a memory storing computer program instructions; the processor implements the power grid operating state simulation method in the first aspect or any of the embodiments of the first aspect when executing the computer program instructions.

[0068] Based on the same inventive concept, in a fourth aspect, the embodiments of the present application further provide a computer storage medium, and the computer readable storage medium stores computer program instructions, and the computer program instructions are executed by a processor to implement the power grid operating state simulation method in the first aspect or any of the embodiments of the first aspect.

[0069] Based on the same inventive concept, in a fifth aspect, the embodiments of the present application further provide a computer program product, and instructions in the computer program product are executed by a processor of a device to enable the device to execute the power grid operating state simulation method in the first aspect or any of the embodiments of the first aspect.

[0070] The simulation method, device, equipment, medium and program product of the power grid operation state provided by the embodiments of the present application can more finely simulate the actual operation state of the power distribution network, which is beneficial to subsequent in-depth analysis of the operation state of the power distribution network, and can also integrate multiple states into a system for comprehensive calculation. Then, in the correspondence between the preset operation state and the first time sequence data, the target first time sequence data corresponding to the preset operation state information is obtained, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network, and the first time sequence data of each phase voltage signal has an associated relationship. Next, the target first time sequence data is subjected to phase-to-magnitude transformation to obtain first magnitude data, wherein each module component in the first magnitude data is independent of each other, so that the decoupling of the target first time sequence data can be realized through phase-to-magnitude transformation, and since each module component is independent of each other, subsequent calculation and analysis are facilitated. Next, based on each unit transmission parameter information, second magnitude data formed after the first magnitude data is transmitted through each unit transmission path in the preset total transmission path is calculated. Next, the second magnitude data is subjected to magnitude-to-phase transformation, which can realize coupling reconstruction to obtain target second time sequence data, and the target second time sequence data includes second time sequence data of each phase voltage signal, and the second time sequence data of each phase voltage signal has an associated relationship, which restores the associated relationship between each phase voltage signal, so that the simulation result is more in line with the actual situation. Through the acquisition of the preset operation state information of the power distribution network and the unit transmission parameter information of each unit transmission path, the embodiments of the present application can refine the simulation result, realize more accurate simulation of the power distribution network state; at the same time, by using the phase-to-magnitude transformation and magnitude-to-phase transformation technology, the decoupling and coupling reconstruction of the data are realized, the complexity of the calculation is reduced, the space overhead in the calculation process is reduced, a complete simulation process is constructed, a unified simulation framework is formed, which can be applied to the electromagnetic state simulation calculation of the power distribution network in a long time period or multiple state superpositions, can realize the electromagnetic simulation of the normal-non-normal operation state of the power system power distribution network, and can provide strong support for the state monitoring and fault diagnosis of the power distribution network. BRIEF DESCRIPTION OF DRAWINGS

[0071] Other features, objects, and advantages of the present application will become more apparent from the following detailed description of non-limiting embodiments with reference to the drawings, in which like or similar elements are designated with like or similar reference numerals, and in which the drawings have not necessarily been drawn to scale.

[0072] Figure 1 is a flowchart of the simulation method of the power grid operation state provided by the embodiments of the present application;

[0073] Figure 2is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0074] Figure 3 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0075] Figure 4 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0076] Figure 5 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0077] Figure 6 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0078] Figure 7 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0079] Figure 8 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0080] Figure 9 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0081] Figure 10-A is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0082] Figure 10-B is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0083] Figure 10-C is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0084] Figure 11 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0085] Figure 12 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application;

[0086] Figure 13is a normal- abnormal state modulus comparison chart in phase-mode conversion in the power grid operation state simulation method provided by the embodiment of the present application;

[0087] Figure 14 is a result chart of first modulus data obtained by phase-mode conversion of target first time sequence data in the power grid operation state simulation method provided by the embodiment of the present application;

[0088] Figure 15-A is a waveform schematic diagram of second time sequence data received in a normal state in the power grid operation state simulation method provided by the embodiment of the present application;

[0089] Figure 15-B is a waveform schematic diagram of second time sequence data received when a single-phase grounding fault occurs in the power grid operation state simulation method provided by the embodiment of the present application;

[0090] Figure 15-C is a waveform schematic diagram of second time sequence data received when an overvoltage fault occurs in the power grid operation state simulation method provided by the embodiment of the present application;

[0091] Figure 16 is a structural schematic diagram of the power grid operation state simulation device provided by the embodiment of the present application;

[0092] Figure 17 is a structural schematic diagram of the power grid operation state simulation device provided by the embodiment of the present application. DETAILED DESCRIPTION

[0093] The features and exemplary embodiments of various aspects of the present application will be described in detail below, in order to make the purposes, technical solutions and advantages of the present application more clear and apparent, the present application will be further described in detail below in combination with the drawings and specific embodiments. It should be understood that the specific embodiments described herein are only configured to explain the present application, and are not configured to limit the present application. For those skilled in the art, the present application can be implemented without some of these specific details. The following description of the embodiments is only to provide a better understanding of the present application by showing examples of the present application.

[0094] It should be noted that the relative terms, such as first and second, and the like, are used herein only to distinguish one entity or operation from another, and do not necessarily require or imply any actual such relationship or order between or among these entities or operations. Moreover, the terms "comprises", "comprising", or any other variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by an "includes" statement does not exclude the existence of additional elements of the same nature as those included in the process, method, article, or apparatus.

[0095] It should be understood that the term "and / or" used herein only describes an association relationship of associated objects, which means that there can be three relationships, for example, A and / or B can represent three cases of A alone, A and B together, and B alone. In addition, the character " / " herein generally represents an "or" relationship between the front and rear associated objects.

[0096] Various modifications and changes can be made to the present application without departing from the spirit or scope of the application. It is therefore intended that the present application cover all such modifications and changes as fall within the scope of the corresponding claims (technical solutions claimed to be protected) and their equivalents. It should be noted that the embodiments provided by the present application can be combined with each other without contradiction, if possible.

[0097] Before describing the technical solutions provided by the embodiments of the present application, the problems existing in the related art will be described in detail to facilitate the understanding of the embodiments of the present application:

[0098] Power distribution network operation state monitoring is an important means to improve user power supply reliability, which can effectively avoid power loss and increase operation and maintenance costs caused by long-time fault finding.

[0099] In the related art, simulation of the operation state of the power system depends on industrial software, and the calculation of the related industrial software is mainly for global expansion, which cannot customize the distribution parameters of a single transmission grid, the simulation result lacks detailed information, and cannot meet the demand for in-depth analysis of the operation state of the power distribution network; and it cannot be calculated in a unified framework, so it cannot realize electromagnetic state simulation calculation of the power distribution network in a long period of time and in combination with multiple states, and lacks research on electromagnetic simulation framework for normal-non-normal operation states of the power distribution network of the power system. Therefore, when simulating the operation state of the power distribution network, the related art has the problems of insufficient detailed information of the simulation result and inability to realize long-time, multi-state superimposed electromagnetic state simulation.

[0100] Based on this, the embodiment of the application provides a power grid operation state simulation method, device, equipment, medium and program product, which can refine the simulation result, realize more accurate power distribution network state simulation, and is also applicable to power distribution network electromagnetic state simulation calculation in a long time period or a plurality of state superpositions, and can provide strong support for power distribution network state monitoring and fault diagnosis.

[0101] The power grid operation state simulation method provided by the embodiment of the application will be described in detail below with reference to the accompanying drawings.

[0102] Figure 1 is a flowchart of a power grid operation state simulation method provided by the embodiment of the application, as Figure 1 shown, the method can include steps S110-S140.

[0103] S110, obtaining preset operation state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path.

[0104] The power distribution network can distribute electric energy to each user, including residential, commercial and industrial users.

[0105] The preset operation state information refers to the operation state information of the power system set according to actual demand or assumed conditions in the simulation process. For example, the operation state of the power distribution network can include normal operation state, single-phase grounding state, operating overvoltage state, etc.

[0106] The preset total transmission path is a self-defined total transmission path. The total transmission path refers to the entire power transmission network path through which electric energy is transmitted from the power distribution network to the user terminal.

[0107] The unit transmission path is a basic unit constituting the total transmission path.

[0108] The unit transmission parameter information is the transmission parameter information in the unit transmission path, such as the resistance value of the distributed resistance, the inductance value of the distributed inductance, the capacitance value of the distributed capacitance, etc. in the unit transmission path. The unit transmission parameter information of each unit transmission path is different and can be self-defined, so as to refine the simulation result and realize more accurate power distribution network state simulation.

[0109] Specifically, first, the self-defined preset operation state information about the power distribution network needs to be obtained, for example, the normal operation state, single-phase grounding state, operating overvoltage state, etc. can be self-defined, and at the same time, the specific transmission parameters of each subdivided unit transmission path in the preset total transmission path, such as resistance, inductance, capacitance, etc. also need to be obtained. These parameters are crucial for accurately describing the transmission characteristics of electric energy in the power distribution network.

[0110] S120, in the preset operating state and the corresponding relationship of the first time sequence data, obtain the target first time sequence data corresponding to the preset operating state information, and perform phase-mode transformation on the target first time sequence data to obtain first modulus data, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network.

[0111] The preset operating state and the corresponding relationship of the first time sequence data is a mapping relationship between a series of preset operating states and time sequence data (i.e. first time sequence data) corresponding to these states. The first time sequence data is data of each phase voltage signal changing with time, which is used to reflect the dynamic behavior of the power system in different states.

[0112] The target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network, and the first time sequence data of each phase voltage signal has a correlation relationship.

[0113] Phase-mode transformation is a method of converting electrical quantities (such as voltage, current) of a three-phase system from phase components (A, B, C phases) to modulus components (such as positive, negative, zero sequence components), which can realize decoupling of the three-phase system and help simplify subsequent analysis and calculation, especially when dealing with asymmetric faults and unbalanced loads.

[0114] The first modulus data is modulus data corresponding to the original target first time sequence data after phase-mode transformation processing, and each modulus component in the first modulus data is independent of each other. The first modulus data can include independent positive, negative, and zero sequence component data.

[0115] Specifically, according to the preset operating state and the corresponding relationship of the first time sequence data, the target first time sequence data matching the specific preset operating state can be found. The target first time sequence data mainly records the change of each phase voltage signal with time in the power distribution network. The first time sequence data of each phase voltage signal in the target first time sequence data is correlated. Subsequently, the target first time sequence data can be decoupled, i.e. phase-mode transformation processing, to obtain first modulus data, which includes independent positive, negative, and zero sequence component data, facilitating power system analysis under conditions such as asymmetric faults and unbalanced loads.

[0116] S130, based on each unit transmission parameter information, calculate second modulus data formed by the first modulus data transmitted through each unit transmission path in the preset total transmission path.

[0117] Each unit transmission path refers to each transmission unit subdivided from the preset total transmission path. Each unit transmission path has its specific transmission parameter.

[0118] The second modulus data is new modulus data calculated according to the unit transmission parameter information after the first modulus data sequentially passes through each unit transmission path in the preset total transmission path, and reflects the electrical quantity situation of the electric energy after passing through the transmission path.

[0119] Specifically, the specific parameter information of each unit transmission path is used to calculate the new modulus data, i.e., the second modulus data, of the initial first modulus data after sequentially passing through each unit path on the preset total transmission path, which can refine the simulation results and more accurately reflect the changes and losses in the electric energy transmission process.

[0120] S140, performing modulus-phase conversion on the second modulus data to obtain target second time sequence data, the target second time sequence data including second time sequence data of each phase voltage signal.

[0121] The modulus-phase conversion is a method of coupled reconstruction, which is used to convert modulus data (such as positive sequence, negative sequence, and zero sequence components) back to phase components (such as A, B, and C phase voltages or currents), and can convert the results obtained from modulus analysis back to phase component form which is more intuitive and easier to understand.

[0122] The target second time sequence data includes second time sequence data of each phase voltage signal, and the second time sequence data of each phase voltage signal has a correlation relationship. In other words, the second time sequence data of each phase voltage signal is data of the voltage signal of each phase (A, B, and C phase) changing with time after the first modulus data is transmitted and modulus-phase converted.

[0123] Specifically, by performing modulus-phase conversion on the second modulus data (including positive sequence, negative sequence, and zero sequence components), it can be converted back to the time sequence data form of each phase voltage signal. These converted data can be referred to as target second time sequence data, and the target second time sequence data records in detail the change of each phase voltage signal with time (second time sequence data).

[0124] The preset state information of the power distribution network and the unit transmission parameter information of each unit transmission path in the preset total transmission path are acquired, so that the actual operation state of the power distribution network can be simulated in more detail, which is beneficial to subsequent in-depth analysis of the operation state of the power distribution network, and various states can be integrated into a system for comprehensive calculation. Then, in the correspondence between the preset operation state and the first time sequence data, target first time sequence data corresponding to the preset operation state information is acquired, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network, and the first time sequence data of each phase voltage signal has a correlation relationship. Next, phase-mode transformation is performed on the target first time sequence data to obtain first modulus data, wherein each modulus component in the first modulus data is independent of each other, so that decoupling of the target first time sequence data can be realized through phase-mode transformation, and subsequent calculation and analysis are facilitated due to the independence of each modulus component. Next, based on the unit transmission parameter information of each unit transmission path, second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated. Next, the second modulus data is subjected to mode-phase transformation to realize coupling reconstruction and obtain target second time sequence data, wherein the target second time sequence data includes second time sequence data of each phase voltage signal, and the second time sequence data of each phase voltage signal has a correlation relationship, so that the correlation relationship between each phase voltage signal is restored, and the simulation result is more in line with the actual situation. Through acquisition of the preset operation state information of the power distribution network and the unit transmission parameter information of each unit transmission path, the simulation result can be refined, and more accurate power distribution network state simulation can be realized. At the same time, through the use of phase-mode transformation and mode-phase transformation technology, decoupling and coupling reconstruction of data are realized, the complexity of calculation is reduced, the space overhead in the calculation process is reduced, a complete simulation process is constructed, a unified simulation framework is formed, and the simulation framework can be applied to electromagnetic state simulation calculation of the power distribution network in a long time period or under superposition of various states, electromagnetic simulation of the normal-non-normal operation state of the power distribution network can be realized, and strong support can be provided for state monitoring and fault diagnosis of the power distribution network.

[0125] It should be noted that, according to the preset state, corresponding phase voltage signal time sequence data is acquired, phase-mode transformation technology is applied to convert the phase voltage signal time sequence data into independent modulus component data, new modulus data is obtained based on the unit transmission parameter and the change of the modulus component in the transmission path, and finally, mode-phase transformation is used to restore the correlation relationship between the phase voltage signals, so as to form a complete simulation process. Not only is the simulation result refined, but also the calculation complexity is reduced through the use of phase-mode and mode-phase transformation technology, more accurate power distribution network state simulation is realized, and strong support is provided for state monitoring and fault diagnosis.

[0126] Figure 2is another flowchart of the simulation method of the power grid operation state provided in the embodiments of the present application.

[0127] In some embodiments, as shown in Figure 2 The step S110 of acquiring the preset operation state information of the power distribution network and the unit transmission parameter information of each unit transmission path in the preset total transmission path can include steps S111-S113.

[0128] S111, displaying a first page.

[0129] Specifically, the first page can be displayed to the user, and the user can input the customized information in the first page, and the customized information is taken as the first input in the first page.

[0130] S112, in response to the first input in the first page, acquiring the preset operation state information of the power distribution network and the preset operation environment configuration information in the first input.

[0131] Specifically, in the case that the user inputs the customized information in the first page, the specific content of the customized information, such as the preset operation state information of the power distribution network and the preset operation environment configuration information, can be acquired.

[0132] S113, according to the preset operation environment configuration information, calculating the unit transmission parameter information of each unit transmission path in the preset total transmission path.

[0133] Specifically, according to the customized preset operation environment configuration information, the specific unit transmission parameter information of each unit transmission path in the preset total transmission path can be calculated.

[0134] The embodiments of the present application can flexibly define the preset operation state information and the preset operation environment configuration information of the power distribution network to be simulated by displaying the first page to the user to input the customized information, and the specific unit transmission parameter information of each unit transmission path in the preset total transmission path can be accurately calculated accordingly, thereby improving the accuracy and flexibility of the power grid operation state simulation.

[0135] In some embodiments, the preset operation environment configuration information includes the terrain information, the line distribution information and the line material information in the preset total transmission path. The unit transmission parameter information includes the resistance value of the distributed resistance, the inductance value of the distributed inductance, the capacitance value of the distributed capacitance and the conductance value of the distributed conductance in the unit transmission path.

[0136] Specifically, according to the terrain information, the line distribution information and the line material information in the preset total transmission path in the preset operation environment configuration information, the unit transmission parameter information of each unit transmission path in the preset total transmission path, such as the resistance value of the distributed resistance, the inductance value of the distributed inductance, the capacitance value of the distributed capacitance and the conductance value of the distributed conductance, can be calculated.

[0137] The preset operation environment configuration information including the terrain information, the line distribution information and the line material information can be used for accurately calculating unit transmission path information such as distribution resistance, distribution inductance, distribution capacitance and distribution conductance, refining the simulation result, more accurately simulating the actual operation state of the power distribution network, and facilitating the accurate analysis of the operation state of the power distribution network.

[0138] Figure 3 FIG. 6 is another flow diagram of the simulation method of the operation state of the power grid provided in the embodiments of the present application.

[0139] In some embodiments, as shown in FIG. 6, the first input further includes first configuration information of a preset total transmission path and second configuration information of a unit transmission path. The simulation method of the operation state of the power grid can further include steps S151 and S152. Figure 3

[0140] S151, obtaining the first configuration information of the preset total transmission path and the second configuration information of the unit transmission path in the first input.

[0141] The configuration information of the preset total transmission path can be referred to as the first configuration information, for example, the start position and the end position of the total transmission path. In other words, the spatial limit of the total transmission path can be customized.

[0142] The configuration information of the unit transmission path can be referred to as the second configuration information, for example, the length of the unit transmission path. In other words, the spatial precision of the total transmission path can be customized.

[0143] Specifically, the user can input the preset operation state information and the preset operation environment configuration information of the power distribution network in the first page, and can also input the first configuration information of the preset total transmission path and the second configuration information of the unit transmission path, to provide more comprehensive data support for subsequent simulation calculation.

[0144] S152, determining the preset total transmission path according to the first configuration information, and determining the unit transmission path according to the second configuration information.

[0145] Specifically, the preset total transmission path in the entire power grid simulation is determined by analyzing the first configuration information input by the user, and meanwhile, the specific unit transmission path under the total path can be determined according to the second configuration information, to provide a basis for subsequent simulation of the operation state of the power grid.

[0146] The embodiments of the present application can more accurately define the power grid structure by allowing the user to input the configuration information including the preset total transmission path and the unit transmission path, to provide more comprehensive data support for subsequent simulation calculation. ​

[0147] Figure 4 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application.

[0148] In some embodiments, as shown in Figure 4 , the phase-to-magnitude transformation on the target first time-series data in step S120 to obtain the first magnitude data can include:

[0149] multiplying the target first time-series data by a preset transformation matrix to obtain the first magnitude data.

[0150]

[0151] wherein Q is the transformation matrix.

[0152] Specifically, by multiplying the target first time-series data by the preset transformation matrix Q, the phase-to-magnitude transformation on the target first time-series data is realized, thereby obtaining the first magnitude data.

[0153] The embodiment of the present application utilizes the preset transformation matrix Q to perform the phase-to-magnitude transformation on the target first time-series data, which can efficiently calculate the first magnitude data, realizes the mutual decoupling of three-phase voltages, simplifies the calculation process, and improves the efficiency of the power grid operation state simulation.

[0154] In some embodiments, referring back to Figure 4 , the magnitude-to-phase transformation on the second magnitude data in step S140 to obtain the target second time-series data can include:

[0155] multiplying the second magnitude data by a preset inverse transformation matrix to obtain the target second time-series data.

[0156]

[0157] wherein Q -1 is the inverse transformation matrix.

[0158] Specifically, by utilizing the preset inverse transformation matrix Q -1 (that is, the inverse matrix of Q) to perform the magnitude-to-phase transformation on the second magnitude data, the coupling reconstruction can be accurately realized.

[0159] The embodiment of the present application utilizes the preset inverse transformation matrix Q -1 to perform the magnitude-to-phase transformation on the second magnitude data, which can accurately realize the coupling reconstruction of the phase voltages, can convert the magnitude data (such as positive sequence, negative sequence, and zero sequence components) back to the phase components (such as A, B, and C phase voltages or currents), and thereby convert the results obtained from the magnitude analysis back to the phase component form which is more intuitive and easier to understand.

[0160] Figure 5is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application.

[0161] In the first embodiment, as shown in Figure 5 S130, the step of calculating the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information can include steps S131-S133.

[0162] S131, obtaining information of the preset terminal reflection coefficient.

[0163] The terminal reflection coefficient is a parameter for measuring the degree of signal reflection at the end of the transmission line.

[0164] Specifically, the information of the preset terminal reflection coefficient, i.e., the parameter information for describing the degree of signal reflection at the end of the transmission line, can be set according to the terminal condition.

[0165] S132, calculating the terminal matching impedance value according to the information of the terminal reflection coefficient.

[0166] Specifically, according to the information of the terminal reflection coefficient, the matching impedance value of the terminal can be calculated.

[0167] S133, calculating the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information and the terminal matching impedance value.

[0168] Specifically, using the specific unit transmission parameter information of each unit transmission path and the impedance value of the terminal, the new modulus data, i.e., the second modulus data, of the first modulus data after being transmitted through each unit transmission path in the preset total transmission path can be calculated.

[0169] By obtaining the preset terminal reflection coefficient information, calculating the terminal matching impedance value, and calculating the second modulus data based on the information, the embodiment of the present application can more accurately simulate the signal transmission process in the power grid and improve the accuracy and reliability of the simulation result.

[0170] Figure 6 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application.

[0171] In some embodiments, as shown in Figure 6 S132, the step of calculating the terminal matching impedance value according to the information of the terminal reflection coefficient can include steps S1321 and S1322.

[0172] S1321, calculate the equivalent impedance value of the last unit transmission path according to the unit transmission parameter information of the last unit transmission path in the preset total transmission path.

[0173] Specifically, the equivalent impedance value of the last unit transmission path can be calculated according to the unit transmission parameter information of the last unit transmission path in the preset total transmission path, which describes the impedance characteristics of the signal in the path transmission.

[0174] S1322, in the case of terminal reflection coefficient being 0, the equivalent impedance value of the last unit transmission path is taken as the terminal matching impedance value. Please refer to formula 3.

[0175]

[0176] Wherein, Z L is the terminal matching impedance value, Zc is the equivalent impedance value of the last unit transmission path, L, R, C and G are the unit transmission parameter information of the last unit transmission path.

[0177] Wherein, when the terminal reflection coefficient is 0, it means that the signal is not reflected at the end of the transmission line, that is, the signal is completely absorbed or the terminal is matched.

[0178] Specifically, when the signal is not reflected at the end of the transmission line, that is, the terminal reflection coefficient is 0, the equivalent impedance value Zc of the last unit transmission path in the preset total transmission path can be taken as the terminal matching impedance value Z L . Here, Zc is calculated based on the unit transmission parameter information (including inductance L, resistance R, capacitance C and conductance G) of the last unit transmission path, which describes the impedance characteristics of the signal in the path transmission.

[0179] The embodiment of the application can more accurately simulate the transmission of the signal in the power grid by calculating the equivalent impedance value of the last unit transmission path and directly taking it as the terminal matching impedance value when the terminal reflection coefficient is 0, especially in the special condition that the signal is not reflected at the end of the transmission line (i.e. terminal matching), which improves the accuracy and practicability of the simulation results.

[0180] In some embodiments, the power grid operation state simulation method can further include step S1323. Figure 6

[0181] S1323, in the case of terminal reflection coefficient not being 0, calculate the terminal matching impedance value according to the terminal reflection coefficient. Please refer to formula 4.

[0182]

[0183] ​Wherein, Γ is a terminal reflection coefficient.

[0184] Wherein, when the terminal reflection coefficient is not 0, it means that the signal will be reflected at the end of the transmission line, and the degree of reflection is determined by the value of the terminal reflection coefficient. The inventor found that in the case where the terminal reflection coefficient is not 0, if the equivalent impedance value of the last unit transmission path is simply taken as the terminal matching impedance value, the influence of the reflected signal will be ignored.

[0185] Specifically, in the case where the terminal reflection coefficient is not 0, the terminal reflection coefficient Γ and the equivalent impedance value Zc based on the last unit transmission path can be substituted into formula 4 to calculate the terminal matching impedance value Z L .

[0186] In the embodiment of the present application, when the terminal reflection coefficient is not 0, the signal will be reflected at the end of the transmission line, so the terminal reflection coefficient Γ and the equivalent impedance value Zc based on the last unit transmission path can be substituted into formula 4 to calculate the terminal matching impedance value Z L , which can obtain a more accurate terminal matching impedance value, and improve the accuracy and practicability of the simulation result.

[0187] Figure 7 is another flowchart of the simulation method of the power grid operation state provided by the embodiment of the present application.

[0188] In some embodiments, as shown in Figure 7 , step S133 calculates the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on each unit transmission parameter information and the terminal matching impedance value, which can include steps S1331-S1334.

[0189] S1331, obtain the information of the source internal resistance.

[0190] Specifically, the information of the internal resistance of the power source of the power distribution network can be obtained, which represents the hindering effect of the power source on the electric signal, and the source internal resistance can be denoted as Rs.

[0191] S1332, based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path, calculate the first sub-modulus data formed after the first modulus data is transmitted through the first preset unit transmission path.

[0192] Specifically, the first sub-modulus data formed after the first modulus data is transmitted through the first unit transmission path can be calculated by using the source internal resistance information and the unit transmission parameter information of the first unit transmission path in the preset total transmission path.

[0193] S1333, based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth submodulus data, calculating N+1th submodulus data formed by the first modulus data after transmission via the first N+1 preset unit transmission paths, wherein N is an integer greater than or equal to 1.

[0194] Specifically, the N+1th submodulus data formed by the first modulus data after transmission via the first N+1 unit transmission paths can be calculated according to the unit transmission parameter information of the N+1th unit transmission path on the preset total transmission path and the Nth submodulus data calculated before, wherein N represents an integer sequence number greater than or equal to 1.

[0195] For example, the second submodulus data formed by the first modulus data after transmission via the first two unit transmission paths can be calculated according to the unit transmission parameter information of the second unit transmission path on the preset total transmission path and the first submodulus data calculated before. Then, the third submodulus data formed by the first modulus data after transmission via the first three unit transmission paths can be calculated according to the unit transmission parameter information of the third unit transmission path on the preset total transmission path and the second submodulus data calculated before. Similarly, the calculation can be carried out.

[0196] S1334, in the case that the N+1th unit transmission path is the last unit transmission path, based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth submodulus data and the terminal matching impedance value, calculating the second modulus data formed by the first modulus data after transmission via each unit transmission path in the preset total transmission path.

[0197] Specifically, when the N+1th unit transmission path is the last one in the preset total transmission path, the N+1th submodulus data after the Nth submodulus data via the last unit transmission path, that is, the final second modulus data formed by the first modulus data after transmission via the entire preset total transmission path (including all unit transmission paths), can be calculated by using the unit transmission parameter information of the path, the Nth submodulus data calculated before and the terminal matching impedance value.

[0198] The embodiments of the present application can accurately simulate the transmission process of electrical signals in the power grid by gradually calculating the transmission results of the first modulus data in each unit transmission path (including the first submodulus data, the second submodulus data,..., to the final second modulus data), and particularly considering the influence of the source internal resistance information and the terminal matching impedance value, thereby improving the accuracy of simulation.

[0199] In some embodiments, the step S1332 of calculating the first sub-modulus data formed after the first modulus data is transmitted via the first preset unit transmission path based on the information of the source internal resistance, and the unit transmission parameter information of the first unit transmission path in the preset total transmission path can include:

[0200]

[0201] wherein Δx is preset length information of the unit transmission path, Δt is preset time interval information between time points of the first time sequence data, Rs is information of the source internal resistance, C is unit transmission parameter information of the first unit transmission path, N represents the serial number of the unit transmission path, k represents a time point corresponding to the time sequence data, I represents the current, and V represents the voltage.

[0202] Specifically, the first sub-modulus data can be calculated according to the formula 5, by using the information of the source internal resistance Rs, the unit transmission parameter information of the first unit transmission path in the preset total transmission path, and the preset length information Δx of the unit transmission path and the preset time interval information Δt between the time points of the first time sequence data, to calculate the change of the first modulus data after being transmitted through the first preset unit transmission path, i.e., the first sub-modulus data formed.

[0203] The embodiments of the present application can calculate the change of the first modulus data after being transmitted through the first preset unit transmission path, i.e., the first sub-modulus data formed, by the formula 5, comprehensively considering multiple key factors, including the information of the source internal resistance Rs, the unit transmission parameter information of the first unit transmission path in the preset total transmission path, the preset length information Δx of the unit transmission path, and the preset time interval information Δt between the time points of the first time sequence data, to accurately predict the change of the first modulus data after being transmitted through the first unit transmission path.

[0204] In some embodiments, the step S1333 of calculating the N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data can include:

[0205]

[0206] wherein L0, R0, C0, and G0 are unit transmission parameter information of the N+1th unit transmission path, I N and V N are the Nth sub-modulus data, I N+1 and V N+1 are the N+1th sub-modulus data.

[0207] Specifically, the change of the first modulus data after passing through the first N+1 preset unit transmission paths, i.e., the N+1th sub-modulus data (including current I N and voltage V N ) is accurately calculated based on the unit transmission parameter information (including inductance L0, resistance R0, capacitance C0 and conductance G0) of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data (including current I N+1 and voltage V N+1 ) calculated before. The specific operation mode can be referred to formulas (6) and (7), which describe the transmission characteristics of the electrical signal in the unit transmission path and calculate the state of the electrical signal after passing through the specific path.

[0208] The embodiments of the present application can systematically calculate the change of the first modulus data after passing through the continuous unit transmission paths, i.e., the new sub-modulus data, by using the unit transmission parameter information (inductance, resistance, capacitance and conductance) of each unit transmission path in the preset total transmission path and the sub-modulus data (current and voltage) calculated in the previous step. This process can be realized by a specific mathematical formula (such as formulas 6 and 7), thereby ensuring the accuracy and reliability of the calculation.

[0209] In some embodiments, when the N+1th unit transmission path is the last unit transmission path, the step S1334 can calculate the second modulus data formed after the first modulus data passing through each unit transmission path in the preset total transmission path based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth sub-modulus data and the terminal matching impedance value, which can include:

[0210]

[0211] Wherein, C is the unit transmission parameter information of the last unit transmission path, V N+1 is the second modulus data.

[0212] Specifically, when the N+1th unit transmission path is the last unit transmission path in the preset total transmission path, the final state of the first modulus data after passing through the entire preset total transmission path (i.e., all unit transmission paths), i.e., the second modulus data, can be calculated based on the unit transmission parameter information of the last unit transmission path, the Nth sub-modulus data (including current I N and voltage V N ) calculated before and the terminal matching impedance value Z L by formula 8.

[0213] The embodiment of the present application can comprehensively utilize the unit transmission parameter information of the last unit transmission path, the sub-modulus data (including current I N and voltage V N calculated in the previous step, and the terminal matching impedance value Z L when the electric signal reaches the last unit transmission path of the preset total transmission path, to accurately calculate the final state of the signal after passing through the entire transmission path, i.e., the second modulus data, thereby ensuring the simulation integrity and accuracy of signal transmission.

[0214] Figure 8 is another flowchart of the simulation method of the power grid operating state provided by the embodiment of the present application.

[0215] In some embodiments, as shown in Figure 8 , before the phase-mode transformation of the target first time-series data to obtain the first modulus data, the simulation method of the power grid operating state further includes steps S161-S163.

[0216] S161, in the case that the target first time-series data includes voltage data of the mutation type, determining a mutation duration range in which a time point corresponding to the voltage data of the mutation type is located.

[0217] Specifically, when there is voltage data of the mutation type in the target first time-series data, a mutation duration range can be determined according to the time point at which the mutation voltage data appears, for example, 100 time points can be taken around the time point at which the mutation voltage data appears as the mutation duration range.

[0218] S162, performing interpolation calculation on voltage data corresponding to the first type of time point in the mutation duration range as known data to obtain interpolation voltage data corresponding to the second type of time point in the mutation duration range.

[0219] The first type of time point can be a preset number of time points in the mutation duration range. The time points in the mutation duration range other than the first type of time point can be referred to as the second type of time point. The distribution form of the first type of time point can be set according to requirements, which is not limited by the embodiment of the present application. For example, the first type of time point can be a preset number of time points uniformly distributed in the mutation duration range.

[0220] Specifically, based on the voltage data corresponding to the first type of time point in the known mutation duration range, the voltage data corresponding to the second type of time point in the range can be calculated by using the interpolation method, which can be referred to as interpolation voltage data.

[0221] S163, replacing the voltage data corresponding to the second type of time point in the target first time-series data with the corresponding interpolation voltage data to obtain the updated target first time-series data.

[0222] Specifically, the voltage data originally corresponding to the second type of time point in the target first time sequence data can be replaced by the corresponding interpolation voltage data calculated by interpolation, so as to obtain the updated target first time sequence data.

[0223] The phase-mode transformation of the target first time sequence data in step S120 to obtain the first modulus data can include:

[0224] The phase-mode transformation of the updated target first time sequence data to obtain the first modulus data.

[0225] Specifically, after the voltage data originally corresponding to the second type of time point in the target first time sequence data is replaced by the corresponding interpolation voltage data calculated by interpolation to obtain the updated target first time sequence data, the phase-mode transformation can be performed on the updated target first time sequence data to obtain the first modulus data.

[0226] The embodiments of the present application add the step of processing the sudden voltage data before phase-mode transformation, that is, determining the sudden duration range, performing interpolation calculation to obtain interpolation voltage data, and replacing the voltage data originally corresponding to the second type of time point in the target first time sequence data with the corresponding interpolation voltage data calculated by interpolation, so that the updated target first time sequence data realizes gradual transition, which can be more consistent with the situation that instantaneous voltage and current cannot be suddenly changed in actual situation, so as to obtain more accurate target first time sequence data.

[0227] Figure 9 is another flowchart of the simulation method of the power grid operation state provided by the embodiments of the present application.

[0228] In some embodiments, as shown in Figure 9 The simulation method of the power grid operation state can further include step S170.

[0229] S170, displaying the waveform information corresponding to the target second time sequence data on the second page.

[0230] Specifically, the waveform information corresponding to the target second time sequence data calculated can be displayed on the second page for the user to view or further analyze.

[0231] The embodiments of the present application can make the user more intuitively obtain the simulation result of the power grid operation state by displaying the waveform information corresponding to the target second time sequence data calculated on the second page, including the waveform information of each phase voltage, so as to facilitate the user to view and analyze, and further improve the efficiency and accuracy of power grid management and decision-making.

[0232] In one embodiment, the simulation method of the power grid operation state provided by the embodiments of the present application can include:

[0233] 1) Obtain preset operation state information of the power distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path.

[0234] In one example, custom information can be input in the first page, which can include simulation time, spatial limits (start and end positions of the total transmission path), spatial accuracy (length of the unit transmission path, such as Δx), time accuracy (Δt), and other parameter information. Among them, the time accuracy and spatial accuracy need to meet the stability condition.

[0235] In the case of a preset operation state being a normal operation state, the preset operation state information can include the initial value and initial phase of the three-phase voltage and current in the normal working state. The waveform diagram of the three-phase voltage signal of the excitation source in the normal working state can be seen in Figure 10-A .

[0236] In the case of a preset operation state being a single-phase grounding fault, the preset operation state information can include the fault occurrence time and the phase where the fault occurs. The waveform diagram of the three-phase voltage signal of the excitation source in the single-phase grounding fault can be seen in Figure 10-B .

[0237] In the case of a preset operation state being an overvoltage fault, the preset operation state information can include a fault fitting waveform, including a Gaussian pulse, a single exponential decay pulse, a double exponential decay pulse, and a decaying sinusoidal pulse. The waveform diagram of the three-phase voltage signal of the excitation source in the overvoltage fault can be seen in Figure 10-C , including a Gaussian pulse, a single exponential decay pulse, a double exponential decay pulse, and a decaying sinusoidal pulse.

[0238] It should be noted that for different preset operation state information, such as normal working state, single-phase grounding state, and operating overvoltage state, corresponding functions of line voltage and time can be generated, that is, first time sequence data.

[0239] It should also be noted that considering the uneven distribution of parameters on the power transmission line in the actual scenario, the power distribution network is equivalent using a transmission line model, which can be seen in Figure 11 , Figure 11 is a schematic diagram of the transmission line model used by the unit transmission path provided by the embodiment of the present application. The unit transmission parameter information (resistance R0, inductance L0, capacitance C0, conductance G0) of different unit transmission paths is different.

[0240] 2) In the preset operation state and the corresponding relationship of the first time sequence data, obtain the target first time sequence data corresponding to the preset operation state information, wherein the target first time sequence data includes the first time sequence data of each phase voltage signal of the power distribution network.

[0241] 3) Interpolation replacement.

[0242] In the case where the target first time series data includes voltage data of the mutation type, interpolation voltage data in the mutation duration range is calculated, and the voltage data corresponding to the second type time point in the target first time series data is replaced with the corresponding interpolation voltage data to obtain updated target first time series data.

[0243] In one example, Hermite interpolation is used for gradual transition of initial values. The comparison chart before and after interpolation replacement can be seen in Figure 12 It can be more consistent with the actual situation that instantaneous voltage and current cannot mutate.

[0244] For example, for a single-phase ground fault signal waveform, a range of 100 time steps before or after the voltage mutation point is selected as the mutation duration range, representing the transition state of fault-non-fault. Then, ten points are extracted before or after the voltage mutation point as known points, i.e., as the first type time points in the mutation duration range. Then, using the extracted known points, the voltage data corresponding to the second type time points in the mutation step range is processed by Hermite interpolation, so that the first-order continuity of the result can be guaranteed. Then, the signal generated by the interpolation processing is used to replace the corresponding part of the original signal, and finally the complete fault source signal is synthesized. In other words, the voltage data corresponding to the second type time point in the target first time series data is replaced with the corresponding interpolation voltage data to obtain updated target first time series data.

[0245] 4) The updated target first time series data is subjected to phase-mode transformation to obtain first modulus data, which facilitates further independent calculation of the transmission of each component. The modulus comparison chart of normal-non-normal state during phase-mode transformation can be seen in Figure 13 .

[0246] The result chart of phase-mode transformation of the target first time series data to obtain the first modulus data can be seen in Figure 14 .

[0247] 5) Based on the unit transmission parameter information, the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated. For example, the finite difference time domain algorithm is used to update the voltage and current alternately, and the update formula can be seen in the above formula 6 and formula 7. In other words, since the distribution parameters of the power transmission line in the actual scene are not balanced, the power distribution network can be equivalent using the transmission line model. Based on the transmission line distribution parameter model, time domain difference iteration calculation is performed.

[0248] It should be noted that, in order to avoid additional reflections at the source and terminal side, additional processing is required for the source and terminal point, and a matching resistance is selected at the source and load end, and the impedance value is set to be equal to the characteristic impedance of the line, which can be seen from the above formula 3. Alternatively, if it is desired to add a reflection simulation at the terminal to simulate the connection of the distribution network, the reflection coefficient at the connection can be calculated using the actual topology structure, and the matching resistance value of the terminal can be calculated according to the above formula 4.

[0249] It should also be noted that, according to the impedance values of the source and terminal, the voltage and current values at the source side and the terminal side are updated using the above formula 5 and formula 8, and the voltage parameters at the boundary are set.

[0250] 6) The second modulus data is subjected to modulus-phase transformation to obtain the final actual voltage and current results on the line, i.e., target second time series data, and the target second time series data includes second time series data of each phase voltage signal.

[0251] It should be noted that the calculation result by the finite difference time domain algorithm is modulus, and in order to convert it to the actual voltage and current on the line, phase modulus inverse transformation is required, and the inverse transformation matrix can be seen from formula 2.

[0252] 7) The inverse transformation result is visualized as the final calculation result output. The waveform information of the received second time series data under normal conditions can be seen from Figure 15-A . The waveform information of the received second time series data when a single-phase ground fault occurs can be seen from Figure 15-B . The waveform information of the received second time series data when an overvoltage fault occurs can be seen from Figure 15-C .

[0253] The simulation method of the power grid operating state provided by the embodiment of the application is a normal-non-normal working state electromagnetic simulation method of the distribution network based on the transmission line model combined with the finite difference time domain method. The preset operating state information and unit transmission parameter information of the distribution network are obtained, the target first time series data is subjected to interpolation replacement to process the sudden voltage data, and then phase modulus transformation and transmission calculation are performed. Finally, the actual voltage and current results on the line (target second time series data) are obtained through modulus-phase transformation, and visual output is performed. The power grid under different operating states can be accurately simulated, sudden voltage can be processed, simulation accuracy can be improved, strong support can be provided for the safe and stable operation of the power grid, visual output is convenient for users to intuitively understand and analyze the power grid state. On the one hand, the embodiment of the application decouples the three-phase voltage and current by using phase modulus transformation, reduces the calculation complexity, reduces the space overhead in the calculation process, and unifies the calculation under different conditions into one architecture. On the other hand, the distribution network system is approximated by using the transmission line model, the influence of the distributed parameters is considered in the calculation process, the calculation result is closer to the actual situation, has the characteristics of small space overhead, wide use scene and high calculation accuracy.

[0254] In one example, a simulation of power grid operation was implemented using MATLAB on a host computer with 128GB of RAM, a 3.4GHz CPU, and an NVIDIA RTX 4070Ti graphics card. The experiment tested normal operation, single-phase ground fault, and overvoltage conditions. The simulation scenario used a three-phase AC power system with a rated voltage of 220V. The three phases were 120° out of phase, and the power was transmitted via overhead conductors. Specifically:

[0255] 1) To verify the ability of this invention to perform electromagnetic simulation of normal and abnormal operating states of distribution networks based on the transmission line model and the finite-difference time-domain method, the simulation calculation of the normal operating state is first performed.

[0256] The normal operating voltage is set to three-phase power, 220V, with a phase difference of 120° between the three phases. The excitation source waveform is as follows: Figure 10-A As shown. The transmission line model used in the simulation process is as follows. Figure 11 As shown, the total simulation length is 3000m, with a unit length of 3m. The transmission line parameters are represented as distributed capacitance, distributed inductance, distributed resistance, and distributed transconductance, respectively. The simulation time is one period of 0.02s, with a time step Δt of 0.5e. -8 It satisfies the computational stability condition.

[0257] It should be noted that since the power distribution network transmission system is a three-conductor system, there is mutual coupling during transmission. Related technologies need to consider the mutual inductance and transconductance between the three phases, increasing the computational load and space overhead. In this embodiment, we choose to first decouple the three-phase voltages by performing phase-mode transformation, and then perform calculations to determine the electromagnetic characteristics of the power distribution network under normal and abnormal operating conditions. The phase-mode transformation results of the three-phase electricity at a certain moment are as follows: Figure 13 As shown. Figure 13 The first row represents the three-phase components in normal operation, where the three phases are symmetrically distributed. After phase mode transformation, only the positive-sequence component exists; zero-sequence and reverse-sequence components are absent. Figure 14 The first row shows the result of the phase-mode transformation of the entire time series, after which only the positive-sequence component remains. Since the three phases have been decoupled during the phase-mode transformation, subsequent calculations do not need to consider the coupling between the three phases, greatly reducing the computational load.

[0258] The receiving location is set at a distance of 2000m. A matching resistor is connected to the terminal. Terminal reflection is ignored. After calculation, a phase-mode inverse transformation is performed to obtain the final phasor result as follows: Figure 15-A As shown, the present invention can simulate and calculate the normal working state of the power distribution network. The calculation is stable, and the waveform will not be distorted during transmission when the power system is working stably.

[0259] 2) To verify the ability of this invention to handle fault voltages in power system distribution networks, simulation calculations were performed on single-phase grounding fault voltages.

[0260] The fault voltage is set to a ground fault occurring in phase A during operation. After a brief period, the power supply to the faulted area is cut off by the power system, and the external signal is that the power system has restored stable power transmission. The fault waveform is as follows: Figure 10-B As shown, due to the coupling between the three phases, when a ground fault occurs in phase A, the voltages of phases B and C rise, and the phase difference between the three phases no longer meets the 120-degree requirement.

[0261] Figure 13 The second line shows the phase-mode transformation result at the moment of the fault. Due to the ground fault, the three phases are no longer rotationally symmetrical, and the zero-sequence component in the transformation result is no longer zero. The system simultaneously contains both positive-sequence and zero-sequence components. Figure 14 The second line is a phase-mode transformation of the fault voltage time series. When a fault occurs, the transformed waveform contains both positive-sequence and zero-sequence components.

[0262] It should be noted that, since voltage and current cannot change abruptly in real-world scenarios, and such changes would introduce additional high-frequency components in the calculations, a time step range of 100 steps is selected before or after the voltage abrupt change point as the transition state between fault and non-fault conditions for the single-phase ground fault signal waveform. First, ten points are extracted before or after the voltage abrupt change point as known points. Then, using these ten points, Hermite interpolation is performed on other signals within the 100 time step range to ensure first-order continuity of the results. Finally, the interpolated signal replaces the corresponding part of the original signal, ultimately synthesizing the complete fault source signal. In other words, the voltage data corresponding to the second type of time point in the target first time series data is replaced with the corresponding interpolated voltage data to obtain the updated target first time series data. The comparison between the interpolated result (the target first time series data after interpolation replacement) and the original result (the target first time series data before interpolation replacement) is shown below. Figure 12 As shown.

[0263] The receiving location is set at a distance of 2000m. A matching resistor is connected to the terminal. Terminal reflection is ignored. After calculation, a phase-mode inverse transformation is performed to obtain the final phasor result as follows: Figure 15-B As shown, the observation results reveal that when a ground fault occurs, a voltage surge will excite a high-frequency component on the line. The transmission of this high-frequency component is affected by distributed parameters. Compared with traditional industrial software such as MATLAB's Simulink simulation, the distributed parameters can be considered in the calculation process of this invention, and the simulation results contain more detailed information.

[0264] 3) In order to verify the generalization ability of the electromagnetic simulation of the normal and abnormal working state of the power distribution network of the power system in the application, the overvoltage fault is simulated.

[0265] The overvoltage fault excitation signal waveform is as shown in Figure 10-C The fault signal respectively considers Gaussian pulse, single exponential decay signal, double exponential decay signal and decaying sinusoidal signal, the simulation conditions are consistent with the above normal / abnormal state, and the calculation result is as shown in Figure 15-C The results show that the overvoltage fault can be processed under the framework of the embodiment of the application, and the normal / abnormal working state of the power distribution network of the power system can be simulated, which expands the application scenarios and application scope.

[0266] The power grid operation state simulation method realized by the embodiment of the application successfully simulates the electromagnetic characteristics of the three-phase alternating current of the power system under the normal, single-phase grounding fault and overvoltage state by combining the transmission line model and the finite difference time domain algorithm. The method not only reduces the calculation complexity caused by the coupling between the three phases, but also optimizes the mutation processing of the fault voltage through Hermite interpolation processing, and improves the simulation accuracy. Compared with the traditional simulation tool, the method can more carefully consider the influence of the distributed parameters, and the simulation result contains more detailed information, which verifies the stability and generalization ability of the method in the electromagnetic simulation of the normal and abnormal working state of the power distribution network of the power system.

[0267] It should be noted that the solution scheme based on the industrial software such as Simulink in MATLAB in the related art can only solve the global voltage and current, and cannot define and set the distributed parameter information. In actual scenarios, due to the different materials of the power distribution network lines, the different heights of overhead lines and the different ground features, the distributed parameters in the entire line are not balanced, so the inability to define the distributed parameters will cause some information to be missing in the calculation result. The embodiment of the application can define the distributed parameters on each calculation grid based on the transmission line model using the finite difference time domain algorithm, and the calculation result contains more information and has higher accuracy.

[0268] It should be further noted that the related calculation method does not include various faults in a unified framework for solving, which makes it difficult to simulate the electromagnetic state of the power distribution network with long time and multiple states stacked. The embodiment of the application uses phase-mode transformation to unify the calculation under different conditions into one architecture, and the use scenario is wider.

[0269] It should be further explained that in actual transmission process, since the power distribution network is a three-conductor system, coupling exists among the three phases, and the related technology needs to consider the coupling matrix in the calculation process, which has large calculation amount. The embodiment of the application decouples the three-phase voltage and current through phase-to-magnitude transformation, and only needs to consider the corresponding magnitude in the calculation, without considering the coupling among the three phases, thereby effectively reducing the calculation complexity and reducing the space overhead in the calculation process.

[0270] Based on the same inventive concept, the embodiment of the application further provides a simulation device for power grid operation state, as shown in Figure 16 The device 1600 can include an acquisition module 1610, a conversion module 1620 and a calculation module 1630.

[0271] The acquisition module 1610 is configured to acquire preset operation state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path.

[0272] The conversion module 1620 is configured to acquire target first time sequence data corresponding to the preset operation state information in the preset correspondence between the preset operation state and the first time sequence data, and perform phase-to-magnitude transformation on the target first time sequence data to obtain first magnitude data, wherein the target first time sequence data includes first time sequence data of each phase voltage signal of the power distribution network.

[0273] The calculation module 1630 is configured to calculate second magnitude data formed by the first magnitude data transmitted through each unit transmission path in the preset total transmission path based on each unit transmission parameter information.

[0274] The conversion module 1620 is further configured to perform magnitude-to-phase transformation on the second magnitude data to obtain target second time sequence data, and the target second time sequence data includes second time sequence data of each phase voltage signal.

[0275] In some embodiments, the acquisition module is configured to acquire preset operation state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path, and can be specifically used for:

[0276] Displaying a first page;

[0277] In response to a first input in the first page, acquiring preset operation state information of a power distribution network and preset running environment configuration information in the first input;

[0278] According to the preset running environment configuration information, calculating the unit transmission parameter information of each unit transmission path in the preset total transmission path.

[0279] In some embodiments, the preset running environment configuration information includes terrain information, line distribution information and line material information in a preset total transmission path; and the unit transmission parameter information includes resistance value of distributed resistance, inductance value of distributed inductance, capacitance value of distributed capacitance and conductance value of distributed conductance in a unit transmission path.

[0280] In some embodiments, the first input further includes first configuration information of the preset total transmission path and second configuration information of the unit transmission path; and the apparatus further includes a determination module:

[0281] The acquisition module is further configured to acquire the first configuration information of the preset total transmission path and the second configuration information of the unit transmission path in the first input;

[0282] The determination module is configured to determine the preset total transmission path according to the first configuration information, and determine the unit transmission path according to the second configuration information.

[0283] In some embodiments, the conversion module is configured to perform phase-mode conversion on the target first time-series data to obtain first modulus data, and specifically can be configured to:

[0284] multiply the target first time-series data by a preset conversion matrix to obtain the first modulus data;

[0285]

[0286] wherein Q is the conversion matrix.

[0287] In some embodiments, the conversion module is configured to perform mode-phase conversion on the second modulus data to obtain target second time-series data, and specifically can be configured to:

[0288] multiply the second modulus data by a preset inverse conversion matrix to obtain the target second time-series data;

[0289]

[0290] wherein Q -1 is the inverse conversion matrix.

[0291] In some embodiments, the calculation module is configured to calculate second modulus data formed by the first modulus data after being transmitted through each unit transmission path in the preset total transmission path based on each unit transmission parameter information, and specifically can be configured to:

[0292] acquire information of a preset terminal reflection coefficient;

[0293] calculate a terminal matching impedance value according to the information of the terminal reflection coefficient;

[0294] The second modulus data is formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information and the terminal matching impedance value.

[0295] In some embodiments, the computing module is configured to calculate the terminal matching impedance value according to the information of the terminal reflection coefficient, and specifically can be configured to:

[0296] According to the unit transmission parameter information of the last unit transmission path in the preset total transmission path, the equivalent impedance value of the last unit transmission path is calculated.

[0297] In the case that the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is taken as the terminal matching impedance value.

[0298]

[0299] Wherein, Z L is the terminal matching impedance value, Zc is the equivalent impedance value of the last unit transmission path, L, R, C and G are the unit transmission parameter information of the last unit transmission path.

[0300] In some embodiments, the apparatus further comprises:

[0301] The computing module is further configured to calculate the terminal matching impedance value according to the terminal reflection coefficient in the case that the terminal reflection coefficient is not 0.

[0302]

[0303] Wherein, Γ is the terminal reflection coefficient.

[0304] In some embodiments, the computing module is configured to calculate the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information and the terminal matching impedance value, and specifically can be configured to:

[0305] Obtain the information of the source internal resistance;

[0306] According to the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path, the first sub-modulus data formed after the first modulus data is transmitted through the first preset unit transmission path is calculated.

[0307] According to the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data, the N+1th sub-modulus data formed after the first modulus data is transmitted through the first N+1 preset unit transmission paths is calculated, wherein N is an integer greater than or equal to 1.

[0308] In the case that the N+1th unit transmission path is the last unit transmission path, the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path is calculated based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth modulus data and the terminal matching impedance value.

[0309] In some embodiments, the calculation module is configured to calculate the first modulus data formed after the first modulus data is transmitted through the first preset unit transmission path based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path, and specifically can be used for:

[0310]

[0311] Where Δx is preset length information of the unit transmission path, Δt is preset time interval information between time points of the first time sequence data, Rs is information of the source internal resistance, C is unit transmission parameter information of the first unit transmission path, N represents the serial number of the unit transmission path, k represents the time point, I represents the current, and V represents the voltage.

[0312] In some embodiments, the calculation module is configured to calculate the N+1th modulus data formed after the first modulus data is transmitted through the first N+1 preset unit transmission paths based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth modulus data, and specifically can be used for:

[0313]

[0314] Where L0, R0, C0 and G0 are unit transmission parameter information of the N+1th unit transmission path, I N and V N are the Nth modulus data, I N+1 and V N+1 are the N+1th modulus data.

[0315] In some embodiments, the calculation module is configured to calculate the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth modulus data and the terminal matching impedance value in the case that the N+1th unit transmission path is the last unit transmission path, and specifically can be used for:

[0316]

[0317] Where C is unit transmission parameter information of the last unit transmission path, V N+1 is the second modulus data.

[0318] In some embodiments, before the conversion module is configured to perform a phase-mode transformation on the target first time-series data to obtain first mode data, the apparatus further comprises an interpolation module and a replacement module:

[0319] The determination module is further configured to, in a case where the target first time-series data comprises mutation type voltage data, determine a mutation duration range in which a time point corresponding to the mutation type voltage data is located;

[0320] The interpolation module is configured to perform interpolation calculation on voltage data corresponding to a first type of time point in the mutation duration range as known data to obtain interpolation voltage data corresponding to a second type of time point in the mutation duration range;

[0321] The replacement module is configured to replace voltage data corresponding to the second type of time point in the target first time-series data with corresponding interpolation voltage data to obtain updated target first time-series data;

[0322] The conversion module is configured to perform a phase-mode transformation on the target first time-series data to obtain first mode data, and can be specifically configured to:

[0323] perform a phase-mode transformation on the updated target first time-series data to obtain first mode data.

[0324] In some embodiments, the apparatus further comprises a display module:

[0325] The display module is configured to display waveform information corresponding to the target second time-series data on a second page.

[0326] The various modules in the power grid operation state simulation apparatus provided by the embodiments of the present application can achieve Figures 1-9 the functions of the various steps of the power grid operation state simulation method provided by the embodiments of the present application and achieve the corresponding technical effects. For brevity, the functions of the various steps of the power grid operation state simulation method provided by the embodiments of the present application will not be described here.

[0327] Figure 17 A hardware structure schematic diagram of a power grid operation state simulation device provided by the embodiments of the present application is shown.

[0328] The power grid operation state simulation device can comprise a processor 1701 and a memory 1702 having stored computer program instructions.

[0329] Specifically, the processor 1701 can comprise a central processing unit (CPU), or an application specific integrated circuit (ASIC), or can be configured as one or more integrated circuits implementing the embodiments of the present application.

[0330] The memory 1702 can include mass storage for data or instructions. As an example and not by way of limitation, the memory 1702 can include a hard disk drive (HDD), a floppy disk drive, flash memory, an optical disc (e.g., a compact disc (CD) or a digital versatile disc (DVD)), a solid-state drive (SSD), a USB drive, or a combination of two or more of these. Where appropriate, the memory 1702 can include removable or non-removable (or fixed) media, where appropriate. The memory 1702 can be internal or external to the emulated power grid operating state device, as appropriate. In particular embodiments, the memory 1702 is non-volatile, solid-state memory.

[0331] The memory can include read-only memory (ROM), random access memory (RAM), magnetic disk storage mediums, optical storage mediums, flash memory devices, electrical, optical, or other physical / tangible memory storage devices. Thus, in general, the memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., a memory device) encoded with software that, when executed (by one or more processors), is operable to

[0332] The processor 1701 implements any one of the emulated power grid operating state methods in the above embodiments by reading and executing computer program instructions stored in the memory 1702.

[0333] In one example, the emulated power grid operating state device further includes a communication interface 1703 and a bus 1704. As shown, the processor 1701, the memory 1702, and the communication interface 1703 are connected through the bus 1704 and complete communication among each other. Figure 17

[0334] The communication interface 1703 is mainly used to realize the communication between various modules, devices, units and / or equipment in the embodiments of the present application.

[0335] ​Bus 1704 includes a hardware, software, or both, that couples components of the grid operating state simulation device to each other. As an example and not by way of limitation, the bus can include an Accelerated Graphics Port (AGP) or other graphics bus, an Extended Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a Hyper Transport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an InfiniBand interconnect, a Low Pin Count (LPC) bus, a memory bus, a Micro Channel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association local (VESA) bus, or another suitable bus or a combination of two or more of these. Where suitable, bus 1704 can include one or more buses. Although this embodiment describes and shows a particular bus, this embodiment contemplates any suitable bus or interconnect.

[0336] The device can perform the grid operating state simulation method in this embodiment based on each unit / component in the grid operating state simulation device, thereby realizing the grid operating state simulation method in this embodiment Figures 1-9 The grid operating state simulation method described.

[0337] In addition, in combination with the grid operating state simulation method in the above embodiments, this embodiment can provide a computer storage medium to realize. The computer storage medium has computer program instructions stored thereon; the computer program instructions are executed by a processor to realize any one of the grid operating state simulation methods in the above embodiments.

[0338] This embodiment also provides a computer program product, instructions in the computer program product are executed by a processor of an electronic device, so that the electronic device performs each process of realizing the grid operating state simulation method embodiments of any one of the above.

[0339] It is to be understood that the application is not limited to the particular configurations and processes described hereinabove and shown in the figures. For the sake of brevity, detailed descriptions of known methods and apparatuses will not be repeated here. In the above embodiments, several specific steps are described and shown as examples. However, the method processes of the present application are not limited to the specific steps described and shown, and various changes, modifications and additions can be made thereto by those skilled in the art without departing from the spirit of the present application, or changing the order of the steps.

[0340] The functional blocks shown in the above described block diagrams can be implemented as hardware, software, firmware, or a combination thereof. When implemented in hardware, they can be, for example, electronic circuits, application specific integrated circuits (ASICs), appropriate firmware, plug-ins, functional cards, etc. When implemented in software, the elements of the present application are program or code segments that are used to perform the required tasks. The program or code segments can be stored in a machine-readable medium or transmitted through a data signal carried in a carrier wave over a transmission medium or communication link. A "machine-readable medium" includes any medium that can store or transport information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, read-only memory (ROM), flash memory, erasable read-only memory (EROM), floppy disks, compact discs (CD-ROM), optical disks, hard disks, fiber-optic media, radio frequency (RF) links, etc. The code segments can be downloaded via computer networks such as the Internet, intranet, etc.

[0341] It is also to be understood that the exemplary embodiments mentioned in the present application describe some methods or systems based on a series of steps or devices. However, the present application is not limited to the order of the above steps, that is, the steps can be performed in the order mentioned in the embodiments, or in an order different from the embodiments, or several steps can be performed simultaneously.

[0342] The computer program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other processing device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other processing device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0343] The above is merely specific implementation of the present application, and those skilled in the art can clearly understand the specific working processes of the system, module and unit described above for the convenience and brevity of description, which can refer to the corresponding processes in the foregoing method embodiments, and will not be described herein. It should be understood that the protection scope of the present application is not limited in this way, and any person skilled in the art can easily think of various equivalent modifications or replacements within the technical range disclosed in the present application, and these modifications or replacements should be covered within the protection scope of the present application.

Claims

1. A method of simulating an operating state of an electrical network, characterized in that The method comprises: obtaining preset operating state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path; in the preset operating state and the corresponding relationship of the first time sequence data, obtaining the target first time sequence data corresponding to the preset operating state information, and performing phase-to-magnitude transformation on the target first time sequence data to obtain first magnitude data, wherein the target first time sequence data comprises first time sequence data of each phase voltage signal of the power distribution network; based on each unit transmission parameter information, calculating second magnitude data formed after the first magnitude data is transmitted through each unit transmission path in the preset total transmission path; performing magnitude-to-phase transformation on the second magnitude data to obtain target second time sequence data, wherein the target second time sequence data comprises second time sequence data of each phase voltage signal.

2. The method of claim 1, wherein, The method comprises: displaying a first page; in response to a first input in the first page, obtaining preset operating state information of a power distribution network and preset operating environment configuration information in the first input; based on the preset operating environment configuration information, calculating unit transmission parameter information of each unit transmission path in a preset total transmission path.

3. The method of claim 2, wherein, The preset operating environment configuration information comprises terrain information, line distribution information and line material information in the preset total transmission path; and the unit transmission parameter information comprises resistance value of a distributed resistance, inductance value of a distributed inductor, capacitance value of a distributed capacitor and conductance value of a distributed conductance in the unit transmission path.

4. The method of claim 2, wherein, The first input further comprises first configuration information of the preset total transmission path and second configuration information of the unit transmission path; and the method further comprises: obtaining the first configuration information of the preset total transmission path and the second configuration information of the unit transmission path in the first input; based on the first configuration information, determining the preset total transmission path; and based on the second configuration information, determining the unit transmission path.

5. The method of claim 1, wherein, The method comprises: multiplying the target first time sequence data by a preset transformation matrix to obtain first magnitude data; wherein Q is the transformation matrix.

6. The method of claim 5, wherein, The method comprises: multiplying the second magnitude data by a preset inverse transformation matrix to obtain target second time sequence data; wherein Q -1 is the inverse transformation matrix.

7. The method of claim 1, wherein, The method comprises: obtaining information of a preset terminal reflection coefficient; based on the information of the terminal reflection coefficient, calculating a terminal matching impedance value; based on each unit transmission parameter information and the terminal matching impedance value, calculating second magnitude data formed after the first magnitude data is transmitted through each unit transmission path in the preset total transmission path.

8. The method of claim 7, wherein, The method comprises: based on the information of the terminal reflection coefficient, calculating a terminal matching impedance value; calculating an equivalent impedance value of the last unit transmission path according to the unit transmission parameter information of the last unit transmission path in the preset total transmission path; in the case that the terminal reflection coefficient is 0, taking the equivalent impedance value of the last unit transmission path as the terminal matching impedance value; Wherein, the Z L is the matching impedance value of the terminal, the Zc is the equivalent impedance value of the last unit transmission path, the L, the R, the C and the G are unit transmission parameter information of the last unit transmission path.

9. The method of claim 8, wherein, further comprising: in the case that the terminal reflection coefficient is not 0, calculating the terminal matching impedance value according to the terminal reflection coefficient; wherein the Γ is the terminal reflection coefficient.

10. The method of claim 7, wherein, the first modulus data after being transmitted through each unit transmission path in the preset total transmission path, comprises: obtaining information of a source internal resistance; calculating first sub-modulus data of the first modulus data after being transmitted through the first preset unit transmission path based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path; calculating N+1th sub-modulus data of the first modulus data after being transmitted through the first N+1 preset unit transmission paths based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data, wherein the N is an integer greater than or equal to 1; in the case that the N+1th unit transmission path is the last unit transmission path, calculating second modulus data of the first modulus data after being transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth sub-modulus data and the terminal matching impedance value.

11. The method of claim 10, wherein, the first modulus data after being transmitted through the first preset unit transmission path based on the information of the source internal resistance and the unit transmission parameter information of the first unit transmission path in the preset total transmission path, comprises: wherein the Δx is preset length information of the unit transmission path, the Δt is preset time interval information between time points of the first time sequence data, the Rs is the information of the source internal resistance, the C is the unit transmission parameter information of the first unit transmission path, the N represents the serial number of the unit transmission path, the k represents the time point, I represents current, and V represents voltage.

12. The method of claim 11, wherein, the N+1th sub-modulus data of the first modulus data after being transmitted through the first N+1 preset unit transmission paths based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path and the Nth sub-modulus data, comprises: Wherein, the L0, R0, C0, G0 are unit transmission parameter information of the N+1th unit transmission path, the I N and the V N are the Nth submodulus data, the I N+1 and the V N+1 are the N+1th submodulus data.

13. The method of claim 12, wherein, the second modulus data of the first modulus data after being transmitted through each unit transmission path in the preset total transmission path based on the unit transmission parameter information of the N+1th unit transmission path in the preset total transmission path, the Nth sub-modulus data and the terminal matching impedance value in the case that the N+1th unit transmission path is the last unit transmission path, comprises: Wherein, the C is the unit transmission parameter information of the last unit transmission path, the V N+1 is the second modulus data.

14. The method of claim 1, wherein, Before the phase-to-magnitude transformation on the target first time-series data to obtain first magnitude data, the method further comprises: In the case that the target first time-series data comprises voltage data of mutation type, determining a mutation duration range in which a time point corresponding to the voltage data of mutation type is located; Performing interpolation calculation on voltage data corresponding to first type time points in the mutation duration range as known data to obtain interpolation voltage data corresponding to second type time points in the mutation duration range; Replacing voltage data corresponding to second type time points in the target first time-series data with corresponding interpolation voltage data to obtain updated target first time-series data; The phase-to-magnitude transformation on the target first time-series data to obtain first magnitude data comprises: Performing phase-to-magnitude transformation on the updated target first time-series data to obtain first magnitude data.

15. The method according to any one of claims 1 to 14, characterized in that, Further comprising: Displaying waveform information corresponding to the target second time-series data on the second page.

16. An apparatus for simulating an operating state of an electric power grid, characterized by Comprise: An acquisition module configured to acquire preset operation state information of a power distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path; A conversion module configured to, in a preset correspondence between operation state and first time-series data, acquire target first time-series data corresponding to the preset operation state information, and perform phase-to-magnitude transformation on the target first time-series data to obtain first magnitude data, wherein the target first time-series data comprises first time-series data of each phase voltage signal of the power distribution network; A calculation module configured to calculate, based on each unit transmission parameter information, second magnitude data formed by the first magnitude data after being transmitted via each unit transmission path in the preset total transmission path; The conversion module is further configured to perform magnitude-to-phase transformation on the second magnitude data to obtain target second time-series data, wherein the target second time-series data comprises second time-series data of each phase voltage signal.

17. A simulation device for power grid operating states, characterized by The device comprises a processor and a memory having computer program instructions stored therein; the processor implements the power grid operation state simulation method according to any one of claims 1 to 15 when executing the computer program instructions.

18. A computer-readable storage medium, characterized in that, The computer readable storage medium has computer program instructions stored thereon, and the computer program instructions are executed by the processor to implement the power grid operation state simulation method according to any one of claims 1 to 15.

19. A computer program product, characterised in that, The instructions in the computer program product are executed by the processor of the device, so that the device can execute the power grid operation state simulation method according to any one of claims 1 to 15.

Images