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

Abstract

The invention discloses a power grid operation state simulation method, device and equipment, a medium and a program product. The method comprises the following steps: acquiring 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; obtaining target first time sequence data corresponding to the preset operation state information, and performing phase-mode transformation on the target first time sequence data to obtain first modulus data; calculating second modulus data formed after the first modulus data is transmitted through each unit transmission path in a preset total transmission path based on each unit transmission parameter information; and carrying out modulus phase transformation on the second modulus data to obtain target second time sequence data. According to the embodiment of the invention, the simulation result can be refined, more accurate state simulation of the power distribution network can be realized, the method can also be suitable for electromagnetic state simulation calculation of the power distribution network in a long time period or with multiple superposed states, and powerful support can be provided for state monitoring and fault diagnosis of the power distribution network.

Description

Technical Field

[0001] The present application belongs to the field of electronic information electromagnetic computing technology, and in particular, relates to a simulation method, device, equipment, medium and program product for power grid operation status. Background Art

[0002] Monitoring the operating status of distribution networks is an important means to improve the reliability of power supply to users. It can effectively avoid power outage losses and increased operation and maintenance costs caused by long-term fault detection.

[0003] The existing simulation of the operating status of the power system relies on industrial software. The calculations of the existing industrial software are mainly carried out globally. It is impossible to customize the distribution parameters of a single transmission grid. The simulation results lack detailed information and cannot meet the needs of in-depth analysis of the operating status of the distribution network. Moreover, it is impossible to incorporate it into a unified framework for calculation. Therefore, it is impossible to realize the simulation calculation of the electromagnetic state of the distribution network with a long period of time and multiple states. There is a lack of research on the electromagnetic simulation framework for the normal-abnormal operating state of the power system distribution network.

[0004] Therefore, when simulating the operating status of the distribution network, the existing technology has the problem of insufficient detailed information of the simulation results and inability to realize the electromagnetic state simulation with long time and multi-state superposition. Summary of the invention

[0005] The embodiments of the present application provide a method, device, equipment, medium and program product for simulating the operating status of a power grid, which can refine the simulation results and achieve more accurate simulation of the distribution network status. It can also be applied to the simulation calculation of the electromagnetic status of the distribution network over a long period of time or with multiple states superimposed, and can provide strong support for the status monitoring and fault diagnosis of the distribution network.

[0006] In a first aspect, an embodiment of the present application provides a method for simulating a power grid operation state, comprising:

[0007] Obtaining preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path;

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

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

[0010] The second modulus data is subjected to an analog-to-phase transformation to obtain target second time series data, wherein the target second time series data includes second time series data of each phase voltage signal.

[0011] In a possible embodiment of the first aspect, obtaining preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path includes:

[0012] Display the first page;

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

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

[0015] In a possible embodiment of the first aspect, the preset operating environment configuration information includes terrain information, line distribution information and line material information in a preset total transmission path; the unit transmission parameter information includes the resistance value of the distributed resistor, the inductance value of the distributed inductance, the capacitance value of the distributed capacitor and the conductance value of the distributed conductance in the unit transmission path.

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

[0017] Acquire first configuration information of a preset total transmission path and second configuration information of a unit transmission path in the first input;

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

[0019] In a possible embodiment of the first aspect, performing phase modulus transformation on the target first time series data to obtain first modulus data includes:

[0020] Multiplying the target first time series data by a preset transformation matrix to obtain first modulus data;

[0021]

[0022] Among them, Q is the transformation matrix.

[0023] In a possible embodiment of the first aspect, performing a modulus-to-phase transformation on the second modulus data to obtain target second time series data includes:

[0024] Multiplying the second modulus data by a preset inverse transformation matrix to obtain target second time series data;

[0025]

[0026] Among them, Q -1 is the inverse transformation matrix.

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

[0028] Obtaining information of preset terminal reflection coefficient;

[0029] According to the information of the terminal reflection coefficient, the terminal matching impedance value is calculated;

[0030] Based on each unit transmission parameter information 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 is calculated.

[0031] In a possible embodiment of the first aspect, calculating the terminal matching impedance value according to information of the terminal reflection coefficient 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] When the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is used as the terminal matching impedance value;

[0034]

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

[0036] In a possible embodiment of the first aspect, the present invention further includes:

[0037] When the terminal reflection coefficient is not 0, the terminal matching impedance value is calculated according to the terminal reflection coefficient;

[0038]

[0039] Where Γ is the terminal reflection coefficient.

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

[0041] Get information about the source internal resistance;

[0042] 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 the first sub-modulus data formed after the first modulus data is transmitted via the first preset unit transmission path;

[0043] 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, calculating the 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] When 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 sub-modulus data and the terminal matching impedance value, 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.

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

[0046]

[0047] Among them, Δx is the preset length information of the unit transmission path, Δt is the preset time interval information between the time points of the first timing data, Rs is the source internal resistance information, C is the unit transmission parameter information of the first unit transmission path, N represents the sequence 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, 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, calculating the N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths includes:

[0049]

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

[0051] In a possible embodiment of the first aspect, when 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 sub-modulus data, and the terminal matching impedance value, calculating the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path, including:

[0052]

[0053] Among them, 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 embodiment of the first aspect, before performing phase modulus transformation on the target first time series data to obtain the first modulus data, the method further includes:

[0055] In the case where the target first time series data includes voltage data of a mutation type, determining a mutation duration range where a time point corresponding to the voltage data of the mutation type is located;

[0056] The voltage data corresponding to the first type of time point within the mutation duration range is used as known data for interpolation calculation to obtain the interpolation voltage data corresponding to the second type of time point within the mutation duration range;

[0057] Replacing the voltage data corresponding to the second type of time point in the target first time series data with the corresponding interpolated voltage data to obtain updated target first time series data;

[0058] Performing phase modulus transformation on the target first time series data to obtain first modulus data includes:

[0059] Perform phase modulus transformation on the updated target first time series data to obtain first modulus data.

[0060] In a possible embodiment of the first aspect, the present invention further includes:

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

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

[0063] An acquisition module, used to acquire preset operation status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path;

[0064] A conversion module, used for obtaining target first time series data corresponding to preset operation state information in the corresponding relationship between the preset operation state and the first time series data, and performing phase-mode conversion on the target first time series data to obtain first mode data, wherein the target first time series data includes first time series data of each phase voltage signal of the distribution network;

[0065] A calculation module, used for calculating, based on each unit transmission parameter information, second modulus data formed after the first modulus data is transmitted via each unit transmission path in a preset total transmission path;

[0066] The conversion module is also used to perform analog-to-phase conversion on the second analog data to obtain target second time series data, and the target second time series data includes the second time series data of each phase voltage signal.

[0067] Based on the same inventive concept, in the third aspect, an embodiment of the present application also provides a simulation device for the operation status of a power grid, the device comprising a processor and a memory storing computer program instructions; when the processor executes the computer program instructions, the simulation method for the operation status of the power grid in the first aspect or any embodiment of the first aspect is implemented.

[0068] Based on the same inventive concept, in the fourth aspect, an embodiment of the present application further provides a computer storage medium, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the simulation method of the power grid operation status in the first aspect or any embodiment of the first aspect is implemented.

[0069] Based on the same inventive concept, in the fifth aspect, an embodiment of the present application also provides a computer program product. When the instructions in the computer program product are executed by the processor of the device, the device can execute the simulation method of the power grid operation status in the first aspect or any embodiment of the first aspect.

[0070] The simulation method, device, equipment, medium and program product of the power grid operation state of the embodiment of the present application, by obtaining the preset operation state information of the distribution network and the unit transmission parameter information of each unit transmission path in the preset total transmission path, can simulate the actual operation state of the distribution network more carefully, which is conducive to the subsequent in-depth analysis of the operation state of the distribution network, and can also incorporate multiple states into a system for comprehensive calculation. Then, in the corresponding relationship between the preset operation state and the first time series data, the target first time series data corresponding to the preset operation state information is obtained, wherein the target first time series data includes the first time series data of each phase voltage signal of the distribution network, and the first time series data of each phase voltage signal have an associated relationship. Then, the target first time series data is subjected to phase-mode transformation to obtain the first modulus data, wherein each modulus component in the first modulus data is independent of each other, thereby the decoupling of the target first time series data can be achieved through phase-mode transformation, and since each modulus component is independent of each other, it is convenient for subsequent calculation and analysis. Next, based on each 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. Next, the second modulus data is subjected to a phase-modulus transformation, which can realize coupling reconstruction and obtain the target second time series data, which includes the second time series data of each phase voltage signal, and the second time series data of each phase voltage signal have a correlation relationship, and the correlation relationship between each phase voltage signal is restored, so that the simulation result is more in line with the actual situation. The embodiment of the present application can refine the simulation results and realize more accurate distribution network state simulation by obtaining the preset operation state information of the distribution network and the unit transmission parameter information of each unit transmission path; at the same time, the phase-modulus transformation and the phase-modulus transformation technology are used to realize the decoupling and coupling reconstruction of the data, reduce the complexity of the calculation, reduce the space overhead in the calculation process, construct a complete set of simulation processes, and form a unified simulation framework, which can be applied to the simulation calculation of the electromagnetic state of the distribution network in a long period of time or with multiple states superimposed, and can realize the electromagnetic simulation of the normal-abnormal operation state of the power system distribution network, and can provide strong support for the state monitoring and fault diagnosis of the distribution network. BRIEF DESCRIPTION OF THE DRAWINGS

[0071] Other features, objects and advantages of the present application will become more apparent by reading the following detailed description of non-limiting embodiments with reference to the accompanying drawings, in which the same or similar reference numerals represent the same or similar features and the accompanying drawings are not drawn to scale.

[0072] Figure 1 It is a flow chart of a method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0073] Figure 2It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0074] Figure 3 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0075] Figure 4 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0076] Figure 5 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0077] Figure 6 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0078] Figure 7 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0079] Figure 8 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0080] Fig. 9 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0081] Figure 10-A It is a waveform diagram of a three-phase voltage signal of an excitation source under a normal working state in a simulation method for a power grid operating state provided in an embodiment of the present application;

[0082] Figure 10-B It is a waveform diagram of a three-phase voltage signal of an excitation source when a single-phase grounding fault occurs in the simulation method of the power grid operation state provided in an embodiment of the present application;

[0083] Figure 10-C It is a waveform diagram of a three-phase voltage signal of an excitation source when an overvoltage fault occurs in the simulation method of the power grid operation state provided in an embodiment of the present application;

[0084] Fig.11 It is a schematic diagram of a transmission line model used in a unit transmission path in a simulation method for a power grid operation state provided in an embodiment of the present application;

[0085] Fig.12 It is a comparison diagram before and after interpolation replacement in the simulation method of the power grid operation state provided in the embodiment of the present application;

[0086] Fig.13It is a comparison diagram of the modulus of normal and abnormal states during phase mode transformation in the simulation method of the power grid operation state provided in the embodiment of the present application;

[0087] Fig.14 It is a result diagram of obtaining first modulus data by performing phase modulus transformation on target first time series data in the simulation method of power grid operation status provided by an embodiment of the present application;

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

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

[0090] Figure 15-C is a waveform diagram of second time series data received when an overvoltage fault occurs in the method for simulating the operation status of a power grid provided in an embodiment of the present application;

[0091] Fig.16 It is a structural schematic diagram of a simulation device for power grid operation status provided in an embodiment of the present application;

[0092] Fig.17 It is a structural diagram of a simulation device for the operation status of a power grid provided in an 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 purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying 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 the need for some of these specific details. The following description of the embodiments is only to provide a better understanding of the present application by illustrating the examples of the present application.

[0094] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the statement "include..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0095] It should be understood that the term "and / or" used in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0096] It is obvious to those skilled in the art that various modifications and changes can be made in the present application without departing from the spirit or scope of the present application. Therefore, the present application is intended to cover modifications and changes of the present application that fall within the scope of the corresponding claims (technical solutions for protection) and their equivalents. It should be noted that the implementation methods provided in the embodiments of the present application can be combined with each other without contradiction.

[0097] Before describing the technical solutions provided by the embodiments of the present application, in order to facilitate the understanding of the embodiments of the present application, the present application first specifically describes the problems existing in the related art:

[0098] Monitoring the operating status of distribution networks is an important means to improve the reliability of power supply to users. It can effectively avoid power outage losses and increased operation and maintenance costs caused by long-term fault detection.

[0099] The simulation of the operation status of the power system in the related technology relies on industrial software, and the calculation of the related industrial software is mainly carried out globally. It is impossible to customize the distribution parameters of a single transmission grid, and the simulation results are insufficient in detail information, which cannot meet the needs of in-depth analysis of the operation status of the distribution network; and it is impossible to incorporate it into a unified framework for calculation, so it is impossible to achieve a long period of time, a combination of multiple states of the distribution network electromagnetic state simulation calculation, lack of research on the electromagnetic simulation framework for the normal-abnormal operation state of the power system distribution network. Therefore, when the related technology simulates the operation status of the distribution network, there is a problem that the simulation results are insufficient in detail information, and it is impossible to achieve a long-term, multi-state superposition electromagnetic state simulation.

[0100] Based on this, the embodiments of the present application provide a simulation method, device, equipment, medium and program product for the operation status of a power grid, which can refine the simulation results and achieve more accurate simulation of the distribution network status. It can also be applied to the simulation calculation of the electromagnetic status of the distribution network over a long period of time or with multiple states superimposed, and can provide strong support for the status monitoring and fault diagnosis of the distribution network.

[0101] The simulation method of the power grid operation status provided in the embodiment of the present application is described in detail below with reference to the accompanying drawings.

[0102] Figure 1 is a flow chart of a simulation method for power grid operation status provided in an embodiment of the present application, such as Figure 1 As shown, the method may include steps S110 to S140.

[0103] S110, obtaining preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path.

[0104] Among them, the distribution network can distribute electrical energy to various users, including residential, commercial and industrial users.

[0105] The preset operation status information refers to the operation status information of the power system set according to actual needs or assumed conditions during the simulation process. For example, the operation status of the distribution network may include normal operation status, single-phase grounding status, operation overvoltage status, etc.

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

[0107] The unit transmission path is the basic unit that constitutes 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 resistor, the inductance value of the distributed inductor, the capacitance value of the distributed capacitor, etc. The unit transmission parameter information of each unit transmission path is different and customizable, so the simulation results can be refined to achieve more accurate distribution network state simulation.

[0109] Specifically, it is first necessary to obtain customized preset operating status information about the distribution network, such as customizable normal operating status, single-phase grounding status, operational overvoltage status, etc. At the same time, it is also necessary to obtain the specific transmission parameters of each subdivided unit transmission path within the preset total transmission path, such as resistance, inductance, capacitance, etc. These parameters are crucial for accurately describing the transmission characteristics of electric energy in the distribution network.

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

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

[0112] The target first time series data includes the first time series data of each phase voltage signal of the power distribution network, and the first time series data of each phase voltage signal are associated with each other.

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

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

[0115] Specifically, according to the correspondence between the preset operating state and the first time series data, the target first time series data matching the specific preset operating state can be found. The target first time series data mainly records the changes of each phase voltage signal in the distribution network over time. The first time series data of each phase voltage signal in the target first time series data are interrelated. Subsequently, the target first time series data can be decoupled, that is, phase module transformation processing can be performed to obtain the first module data. The first module data includes mutually independent positive sequence, negative sequence, and zero sequence component data, which is convenient for power system analysis under conditions such as asymmetric faults and unbalanced loads.

[0116] S130, based on each unit transmission parameter information, calculating second modulus data formed after the first modulus data is transmitted via each unit transmission path in the preset total transmission path.

[0117] The unit transmission paths refer to each transmission unit subdivided from the preset total transmission path, and each unit transmission path has its own specific transmission parameters.

[0118] The second modulus data is new modulus data calculated based on unit transmission parameter information after the first modulus data is transmitted in sequence through each unit transmission path in the preset total transmission path, reflecting the electrical quantity 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 obtained after the initial first modulus data is transmitted through each unit path in turn on the preset total transmission path, that is, the second modulus data. This can refine the simulation results and more accurately reflect the changes and losses in the power transmission process.

[0120] S140, performing an analog-to-phase transformation on the second analog data to obtain target second time series data, where the target second time series data includes second time series data of each phase voltage signal.

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

[0122] The target second time series data includes the second time series data of each phase voltage signal, and the second time series data of each phase voltage signal are associated with each other. In other words, the second time series data of each phase voltage signal is the data of the voltage signal of each phase (A, B, C phase) changing with time after the first modulus data is transmitted and modulus-to-phase converted.

[0123] Specifically, by performing analog-to-phase transformation on the second analog data (including positive-sequence, negative-sequence, and zero-sequence components), it can be converted back into the timing data form of each phase voltage signal. These converted data can be called target second timing data, which records in detail the changes of each phase voltage signal over time (second timing data).

[0124] The embodiment of the present application can simulate the actual operation state of the distribution network more meticulously by obtaining the preset operation state information of the distribution network and the unit transmission parameter information of each unit transmission path in the preset total transmission path, which is conducive to the subsequent in-depth analysis of the operation state of the distribution network, and can also incorporate multiple states into a system for comprehensive calculation. Then, in the corresponding relationship between the preset operation state and the first time series data, the target first time series data corresponding to the preset operation state information is obtained, wherein the target first time series data includes the first time series data of each phase voltage signal of the distribution network, and the first time series data of each phase voltage signal have a correlation relationship. Then, the target first time series data is subjected to phase-mode transformation to obtain the first modulus data, wherein each modulus component in the first modulus data is independent of each other, thereby the decoupling of the target first time series data can be achieved through phase-mode transformation, and since each modulus component is independent of each other, it is convenient for subsequent calculation and analysis. Next, based on the information of each unit transmission parameter, 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. Next, the second modulus data is subjected to a phase-modulus transformation, which can realize coupling reconstruction and obtain the target second time series data, which includes the second time series data of each phase voltage signal, and the second time series data of each phase voltage signal have a correlation relationship, and the correlation relationship between each phase voltage signal is restored, so that the simulation result is more in line with the actual situation. The embodiment of the present application can refine the simulation results and realize more accurate distribution network state simulation by obtaining the preset operation state information of the distribution network and the unit transmission parameter information of each unit transmission path; at the same time, the phase-modulus transformation and the phase-modulus transformation technology are used to realize the decoupling and coupling reconstruction of the data, reduce the complexity of the calculation, reduce the space overhead in the calculation process, construct a complete set of simulation processes, and form a unified simulation framework, which can be applied to the simulation calculation of the electromagnetic state of the distribution network in a long period of time or with multiple states superimposed, and can realize the electromagnetic simulation of the normal-abnormal operation state of the power system distribution network, and can provide strong support for the state monitoring and fault diagnosis of the distribution network.

[0125] It should be noted that the embodiment of the present application obtains the corresponding phase voltage signal time series data according to the preset state, and converts it into independent modulus component data by applying the phase-mode transformation technology, and then calculates the change of the modulus component in the transmission path based on the unit transmission parameter to obtain the new modulus data, and finally restores the correlation between the phase voltage signals through the modulus-phase transformation to form a complete set of simulation processes. Not only does it refine the simulation results, but it also uses the phase-mode and modulus-phase transformation technology to reduce the computational complexity, and achieves a more accurate distribution network state simulation, which is suitable for electromagnetic simulation of long time periods or superposition of multiple states, and provides strong support for state monitoring and fault diagnosis.

[0126] Figure 2It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0127] In some embodiments, Figure 2 As shown, step S110 obtains preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path, and may include steps S111 to S113.

[0128] S111, display the first page.

[0129] Specifically, a first page may be displayed to the user, and the user may input custom information in the first page, and the custom information may be used as the first input in the first page.

[0130] S112: In response to a first input in the first page, obtain preset operating state information and preset operating environment configuration information of the power distribution network in the first input.

[0131] Specifically, when the user enters custom information in the first page, the specific content of the custom information, such as preset operating status information and preset operating environment configuration information of the power distribution network, can be obtained.

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

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

[0134] The embodiment of the present application displays a first page to allow the user to input custom information, thereby flexibly defining the preset operating status information and preset operating environment configuration information of the distribution network to be simulated, and accurately calculating the specific unit transmission parameter information of each unit transmission path in the preset total transmission path, thereby improving the accuracy and flexibility of the power grid operating status simulation.

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

[0136] Specifically, based on the terrain information, line distribution information and line material information in the preset total transmission path in the preset operating environment configuration information, the unit transmission parameter information of each unit transmission path in the preset total transmission path can be calculated, such as the resistance value of the distributed resistor, the inductance value of the distributed inductor, the capacitance value of the distributed capacitor and the conductance value of the distributed conductance.

[0137] In the embodiment of the present application, by including preset operating environment configuration information of terrain information, line distribution information and line material information, the unit transmission parameter information such as distributed resistance, distributed inductance, distributed capacitance and distributed conductance in each unit transmission path can be accurately calculated, the simulation results can be refined, and the actual operating status of the distribution network can be simulated more carefully, which is conducive to the subsequent accurate analysis of the operating status of the distribution network.

[0138] Figure 3 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0139] In some embodiments, Figure 3 As shown, the first input also includes first configuration information of a preset total transmission path and second configuration information of a unit transmission path. The method for simulating the operation state of a power grid may further include steps S151 and S152.

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

[0141] The configuration information of the preset total transmission path may be referred to as first configuration information, for example, the starting position and the ending position of the total transmission path. In other words, the space limit of the total transmission path may be customized.

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

[0143] Specifically, the user can enter the preset operating status information and preset operating environment configuration information of the distribution network in the first page, and can also enter 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 calculations.

[0144] S152, determining a preset total transmission path according to the first configuration information; and determining a 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 parsing the first configuration information input by the user, and at the same time, the specific unit transmission paths under the total path can be clarified based on the second configuration information, providing a basis for subsequent power grid operation status simulation.

[0146] The embodiment of the present application allows the user to input configuration information including a preset total transmission path and a unit transmission path, so as to more accurately define the power grid structure and provide more comprehensive data support for subsequent simulation calculations.

[0147] Figure 4 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0148] In some embodiments, Figure 4 As shown, in step S120, performing phase modulus transformation on the target first time series data to obtain first modulus data may include:

[0149] The target first time series data is multiplied by a preset transformation matrix to obtain first modulus data.

[0150]

[0151] Among them, Q is the transformation matrix.

[0152] Specifically, by multiplying the target first time series data with the preset transformation matrix Q, the phase modulus transformation of the target first time series data is achieved, thereby obtaining the first modulus data.

[0153] The embodiment of the present application utilizes a preset transformation matrix Q to perform phase-modulus transformation on the target first time series data, which can efficiently calculate the first modulus data, realize the mutual decoupling of the three-phase voltages, simplify the calculation process, and improve the efficiency of the power grid operation status simulation.

[0154] In some embodiments, you can continue to see Figure 4 Step S140 performs a modulus-phase transformation on the second modulus data to obtain target second time series data, which may include:

[0155] The second modulus data is multiplied by a preset inverse transformation matrix to obtain target second time series data.

[0156]

[0157] Among them, Q -1 is the inverse transformation matrix.

[0158] Specifically, by using the preset inverse transformation matrix Q -1 (i.e., the inverse matrix of Q) performs a modulus-phase transformation on the second modulus data to accurately achieve coupled reconstruction.

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

[0160] Figure 5It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0161] In some embodiments, Figure 5 As shown, step S130 calculates the second modulus data formed after the first modulus data is transmitted via each unit transmission path in the preset total transmission path based on each unit transmission parameter information, and may include steps S131 to S133.

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

[0163] Among them, the terminal reflection coefficient is a parameter that measures the extent to which the signal is reflected back at the end of the transmission line.

[0164] Specifically, the information of the preset terminal reflection coefficient is obtained, that is, the parameter information preset to describe the degree to which the signal is reflected back at the terminal of the transmission line is obtained, which can be set according to the terminal situation.

[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 obtained by calculation.

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

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

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

[0170] Figure 6 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

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

[0172] S1321, 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.

[0173] Specifically, 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 can be calculated, and this value describes the blocking characteristics of the signal when it is transmitted in the path.

[0174] S1322: When the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is used as the terminal matching impedance value. See Formula 3.

[0175]

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

[0177] When the terminal reflection coefficient is 0, it means that the signal is not reflected at all 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 transmission line terminal, 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 used as the terminal matching impedance value Zc. 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, and describes the obstruction characteristics of the signal when it is transmitted along the path.

[0179] The embodiment of the present application calculates the equivalent impedance value of the last unit transmission path and directly uses it as the terminal matching impedance value when the terminal reflection coefficient is 0, so as to more accurately simulate the transmission of signals in the power grid, especially under the special condition that the signal is not reflected at all at the end of the transmission line (i.e., terminal matching), thereby improving the accuracy and practicality of the simulation results.

[0180] In some embodiments, you can continue to see Figure 6 The simulation method of the power grid operation status may further include step S1323.

[0181] S1323: When the terminal reflection coefficient is not 0, the terminal matching impedance value is calculated according to the terminal reflection coefficient. Please refer to Formula 4.

[0182]

[0183] Where Γ is the terminal reflection coefficient.

[0184] Among them, 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 inventors have found that when the terminal reflection coefficient is not 0, if the equivalent impedance value of the last unit transmission path is simply used as the terminal matching impedance value, the influence of the reflected signal will be ignored.

[0185] Specifically, when 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, since the signal will be reflected at the end of the transmission line, 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 more accurate terminal matching impedance value and improve the accuracy and practicality of simulation results.

[0187] Figure 7 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0188] In some embodiments, Figure 7 As shown, 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, and may include steps S1331 to S1334.

[0189] S1331, obtaining information of source internal resistance.

[0190] Specifically, the information of the internal resistance of the power supply of the distribution network can be obtained, which represents the obstruction of the power supply to the electrical signal. The internal resistance of the source can be recorded 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 via the first preset unit transmission path.

[0192] Specifically, 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, the first sub-modulus data formed after the first modulus data is transmitted through the first unit transmission path can be calculated.

[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 sub-modulus data, calculate the 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.

[0194] Specifically, the N+1th sub-modulus data formed after the first modulus data is transmitted through the first N+1 unit transmission paths can be inferred based on the unit transmission parameter information of the N+1th unit transmission path on the preset total transmission path and the Nth sub-modulus data calculated previously, where N represents an integer serial number greater than or equal to 1.

[0195] For example, the second sub-modulus data formed by the first modulus data after being transmitted through the first two unit transmission paths can be deduced based on the unit transmission parameter information of the second unit transmission path on the preset total transmission path and the first sub-modulus data calculated previously. Then, the third sub-modulus data formed by the first modulus data after being transmitted through the first three unit transmission paths can be deduced based on the unit transmission parameter information of the third unit transmission path on the preset total transmission path and the second sub-modulus data calculated previously. And so on.

[0196] S1334, when 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 sub-modulus data and the terminal matching impedance value, calculate the second modulus data formed after the first modulus data is transmitted through 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 unit transmission parameter information of the path, the Nth sub-modulus data calculated previously and the terminal matching impedance value can be used, and the N+1th sub-modulus data after the Nth sub-modulus data passes through the last unit transmission path, that is, the final second modulus data formed by the first modulus data after passing through the entire preset total transmission path (including all unit transmission paths) is calculated.

[0198] The embodiment of the present application gradually calculates the transmission results of the first modulus data in each unit transmission path (including the first sub-modulus data, the second sub-modulus data, ..., to the final second modulus data), and specially considers the influence of the source internal resistance information and the terminal matching impedance value, so as to accurately simulate the transmission process of the electrical signal in the power grid and improve the accuracy of the simulation.

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

[0200]

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

[0202] Specifically, the first sub-modulus data can be calculated according to Formula 5, by utilizing the source internal resistance information 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 timing data, to calculate the change of the first modulus data after being transmitted through the first preset unit transmission path, that is, the first sub-modulus data is formed.

[0203] The embodiment of the present application can calculate the change of the first modulus data after being transmitted through the first preset unit transmission path, that is, the first sub-modulus data formed, through Formula 5, which comprehensively considers multiple key factors, including the source internal resistance information 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 timing data, and can accurately predict the change of the first modulus data after passing through the first unit transmission path.

[0204] In some embodiments, step S1333 calculates 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, and may include:

[0205]

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

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

[0208] The embodiment 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 formed, 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 implemented through a specific mathematical formula (as shown in 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, step S1334 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 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 may include:

[0210]

[0211] Among them, 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 Nth sub-module data (including the current I N and voltage V N ) and terminal matching impedance value Z L , to calculate the final state of the first modulus data after passing through the entire preset total transmission path (ie, all unit transmission paths), that is, the formed second modulus data.

[0213] The embodiment of the present application can comprehensively utilize the unit transmission parameter information of the last unit transmission path and the sub-module data (including the current I N and voltage V N ) and terminal matching impedance value Z L , accurately calculate the final state of the signal after passing through the entire transmission path, that is, the second modulus data, ensuring the simulation integrity and accuracy of signal transmission.

[0214] Figure 8 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0215] In some embodiments, Figure 8 As shown, before performing phase-modulus transformation on the target first time series data to obtain the first modulus data, the method for simulating the operation status of the power grid further includes steps S161 to S163.

[0216] S161, when the target first time series data includes voltage data of a 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 mutation type voltage data in the target first time series data, a mutation duration range can be determined according to the time point when the mutation voltage data appears. For example, 100 time points around the time point when the mutation voltage data appears can be taken as the mutation duration range.

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

[0219] Among them, the first type of time point can be a preset number of time points within the mutation duration range. The time points other than the first type of time points within the mutation duration range can be called second type of time points. The distribution form of the first type of time points can be set according to demand, and the embodiment of the present application is not limited to this. For example, the first type of time point can be a preset number of time points evenly distributed within the mutation duration range.

[0220] Specifically, based on the voltage data corresponding to the first type of time point within the known mutation duration range, the voltage data corresponding to the second type of time point within the range can be calculated using an interpolation method, and the voltage data can be called interpolated 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 interpolated voltage data to obtain 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 series data may be replaced with the corresponding interpolated voltage data obtained by interpolation calculation, thereby obtaining updated target first time series data.

[0223] In step S120, performing phase-modulus transformation on the target first time series data to obtain first modulus data may include:

[0224] Perform phase modulus transformation on the updated target first time series data to obtain first modulus data.

[0225] Specifically, the voltage data originally corresponding to the second type of time point in the target first time series data is replaced with the corresponding interpolated voltage data obtained by interpolation calculation. After obtaining the updated target first time series data, the updated target first time series data can be subjected to phase modulus transformation to obtain the first modulus data.

[0226] The embodiment of the present application adds a step of processing sudden voltage data before phase mode transformation, that is, determining the sudden change time 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 series data with the corresponding interpolation voltage data obtained by interpolation calculation, so that the updated target first time series data can achieve a gradual transition, which can be more in line with the actual situation where the instantaneous voltage and current cannot change suddenly, thereby obtaining more accurate target first time series data.

[0227] Fig. 9 It is another flow chart of the method for simulating the operation status of a power grid provided in an embodiment of the present application.

[0228] In some embodiments, Fig. 9 As shown, the method for simulating the power grid operation status may further include step S170.

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

[0230] Specifically, the waveform information corresponding to the target second time series data obtained through calculation can be displayed on the second page, so that the user can view or conduct further analysis.

[0231] The embodiment of the present application enables users to more intuitively obtain the simulation results of the power grid operation status by displaying the waveform information corresponding to the calculated target second time series data on the second page, including the waveform information of each phase voltage, which is convenient for users to view and analyze, thereby improving 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 in the embodiment of the present application may include:

[0233] 1) Obtaining preset operating status information of the 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 entered in the first page, including simulation time, space limit (the starting position and the ending position of the total transmission path), space accuracy (the length of the unit transmission path, such as Δx), time accuracy (Δt), and other parameter information. Among them, the time accuracy and space accuracy need to meet the stability condition.

[0235] When the preset operation state is the normal operation state, the preset operation state information may include the initial values ​​and initial phases of the three-phase voltage and current in the normal operation state. The waveform diagram of the three-phase voltage signal of the excitation source in the normal operation state can be found in Figure 10-A .

[0236] When the preset operating state is a single-phase grounding fault, the preset operating state information may include the time of the fault occurrence and the phase where the fault is located. The waveform diagram of the three-phase voltage signal of the excitation source when a single-phase grounding fault occurs can be found in Figure 10-B .

[0237] When the preset operating state is an overvoltage fault, the preset operating state information may include a fault fitting waveform, including a Gaussian pulse, a single exponential decay pulse, a double exponential decay pulse, an attenuated sinusoidal pulse, etc. The waveform diagram of the three-phase voltage signal of the excitation source when an overvoltage fault occurs can be found in Figure 10-C , including Gaussian pulse, single exponential decay pulse, double exponential decay pulse, decayed sine pulse and other waveforms.

[0238] It should be noted that for different preset operating status information, such as normal working state, single-phase grounding state, operation overvoltage state, etc., a corresponding function of online voltage and time, that is, the first time series data, can be generated.

[0239] It should also be noted that considering the imbalance of distribution parameters on the transmission lines of the distribution network in actual scenarios, the transmission line model is used to equate the distribution network, which can be seen in Fig.11 , Fig.11 Schematic diagram of a transmission line model used in a unit transmission path provided in an embodiment of the present application. Different unit transmission paths have different unit transmission parameter information (resistance R0, inductance L0, capacitance C0, and conductance G0).

[0240] 2) In the correspondence between the preset operating state and the first time series data, the target first time series data corresponding to the preset operating state information is obtained, wherein the target first time series data includes the first time series data of each phase voltage signal of the distribution network.

[0241] 3) Interpolation replacement.

[0242] When the target first time series data includes voltage data of a mutation type, interpolated voltage data within the mutation duration range is calculated, and 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 updated target first time series data.

[0243] In an example, Hermite interpolation is used to gradually transition the initial value. The comparison chart before and after interpolation replacement can be seen in Fig.12 , which can better conform to the actual situation that the instantaneous voltage and current cannot change suddenly.

[0244] For example, for the single-phase ground fault signal waveform, a range of 100 time steps is selected before or after the voltage mutation point as the mutation time range to represent the fault-non-fault transition state. Then, ten points are extracted before or after the voltage mutation point as known points, that is, as the first type of time points within the mutation time range. Subsequently, the voltage data corresponding to the second type of time points within the mutation step range are subjected to Hermite interpolation processing using the extracted known points, thereby ensuring the first-order continuity of the result. Then, the signal generated by the interpolation processing replaces the corresponding part of the signal in the original signal, and finally the complete fault source signal is synthesized. In other words, the voltage data corresponding to the second type of time points in the target first time series data is replaced with the corresponding interpolated voltage data to obtain the updated target first time series data.

[0245] 4) Perform phase-mode transformation on the updated target first time series data to obtain first modulus data, so as to facilitate further independent calculation of the transmission of each component. The modulus comparison diagram of normal and abnormal states during phase-mode transformation can be found in Fig.13 .

[0246] The result diagram of the first modulus data obtained by performing phase modulus transformation on the target first time series data can be seen in Fig.14 .

[0247] 5) Based on the transmission parameter information of each unit, calculate the second modulus data formed after the first modulus data is transmitted through each unit transmission path in the preset total transmission path. For example, the voltage and current are updated alternately using the time-domain finite difference algorithm, and the update formula can refer to the above formula 6 and formula 7. In other words, due to the uneven distribution parameters on the transmission lines of the distribution network in actual scenarios, the transmission line model can be used to make the distribution network equivalent, and the time-domain difference iterative calculation is performed based on the transmission line distributed parameter model.

[0248] It should be noted that in order to avoid additional reflections on the source and terminal sides, additional processing is required for the source and terminal points. Matching resistors are added to the source and load ends, and their impedance values ​​are set equal to the line characteristic impedance. See the above formula 3. Alternatively, if you want to add reflections to the terminals to simulate the connection of the distribution network, you can use the actual topology to calculate the reflection coefficient of the connection, and calculate the matching resistance value of the terminal according to the above formula 4.

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

[0250] 6) Performing analog-to-phase transformation on the second analog data to obtain the final actual voltage and current results on the line, that is, the target second time series data, which includes the second time series data of each phase voltage signal.

[0251] It should be noted that the result calculated by the finite difference time domain algorithm is a modulus. In order to convert it into the voltage and current on the actual line, an inverse phase-to-modulus transformation is required. The inverse transformation matrix can be found in Formula 2.

[0252] 7) Visualize the inverse transformation result and output it as the final calculation result. The waveform information of the second time series data received under normal conditions can be found in Figure 15-A The waveform information of the second time series data received when a single-phase grounding fault occurs can be found in Figure 15-B The waveform information of the second time series data received when an overvoltage fault occurs can be found in Figure 15-C .

[0253] The simulation method of the power grid operation state provided in the embodiment of the present application is an electromagnetic simulation method for the normal-abnormal working state of the distribution network based on the transmission line model combined with the time-domain finite difference method. By obtaining the preset operation state information and unit transmission parameter information of the distribution network, the target first time series data is interpolated and replaced to process the sudden voltage data, and then the phase model transformation and transmission calculation are performed. Finally, the actual voltage and current results (target second time series data) on the line are obtained through the model-phase transformation, and the output is visualized. It can accurately simulate the power grid conditions under different operating states, process sudden voltages, improve the simulation accuracy, and provide strong support for the safe and stable operation of the power grid. At the same time, the visual output facilitates users to intuitively understand and analyze the power grid state. On the one hand, the embodiment of the present application uses phase model transformation to decouple the three-phase voltage and current, reduce the calculation complexity, reduce the space overhead in the calculation process, and unify the calculations in different situations into one architecture. On the other hand, the distribution network system is approximated by using the transmission line model, and the influence of distributed parameters can be considered in the calculation process to make the calculation results closer to the actual situation, with small space overhead, wide application scenarios and high calculation accuracy.

[0254] In one example, MATLAB was used on a host with 128GB memory, 3.4GHz CPU, and NVIDIA RTX 4070Ti graphics card to implement a simulation method for the operation status of the power grid. In the experiment, the normal working state, single-phase ground fault state, and overvoltage state were tested respectively. The simulation scenario was the three-phase AC power used in the power system with a rated voltage of 220V. The three phases differed by 120°, and the power transmission method was to use overhead wires. Specifically:

[0255] 1) In order to verify the capability of the present invention to simulate the electromagnetic state of the distribution network under normal and abnormal working conditions based on the transmission line model combined with the finite difference time domain method, the normal working condition is firstly simulated and calculated.

[0256] The normal working voltage is set to three-phase electricity, voltage 220V, the phase difference between the three phases is 120°, the excitation source waveform is as follows Figure 10-A The transmission line model used in the simulation process is shown in Fig.11 As shown in the figure, the total length of the simulation is 3000m and the unit length is 3m. The transmission line parameters are equivalent to distributed capacitance, distributed inductance, distributed resistance and distributed transconductance, the simulation time is one cycle 0.02s, and the time step Δt is 0.5e -8 , which 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 the transmission process. The relevant technology needs to consider the mutual inductance and inter-phase transconductance between the three phases, which increases the amount of calculation and space overhead. In the embodiment of the present application, the three-phase voltage is first decoupled by phase mode transformation, and then the electromagnetic characteristics of the normal-abnormal working state of the distribution network are calculated by calculation. The phase mode transformation result of the three-phase electricity at a certain moment is as follows: Fig.13 shown. Fig.13 The first line is the three-phase component in the normal state, in which the three phases are symmetrically distributed. After the phase mode transformation, only the positive sequence component exists, and there is no zero sequence and reverse sequence component. Fig.14 The first line is the result of phase mode transformation of the entire time series. After the transformation, only the positive sequence component exists. Since the three phases have been decoupled during the phase mode transformation, the subsequent calculation does not need to consider the coupling between the three phases, which greatly reduces the amount of calculation.

[0258] The receiving position is set to be 2000m away, the terminal is connected to a matching resistor, and the terminal reflection is not considered. After calculation, the phase inverse transformation is performed to obtain the final phase result as follows Figure 15-A As shown, it can be known that the present invention can simulate and calculate the normal working state of the distribution network, and the calculation is stable. When the power system works stably, the waveform will not be distorted during the transmission process.

[0259] 2) In order to verify the ability of the present invention in processing the fault voltage of the power system distribution network, a simulation calculation is performed on the single-phase grounding fault voltage.

[0260] The fault voltage is set as a ground fault in phase A during operation. After a short period of time, the power system cuts off the power supply interrupted by the fault, and the power system resumes 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 120 degrees.

[0261] Fig.13 The second line is the result of phase transformation at the fault time. 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 has both positive-sequence and zero-sequence components. Fig.14 The second line is the phase mode transformation of the fault voltage time series. When a fault occurs, the transformed waveform has both positive-sequence and zero-sequence components.

[0262] It should be noted that, since the voltage and current cannot change suddenly in the actual scenario, and the sudden change of voltage and current will cause additional high-frequency components in the calculation, for the single-phase grounding fault signal waveform, a range of 100 time steps is selected before or after the voltage mutation point as the fault-non-fault transition state; first, ten points are extracted before or after the voltage mutation point as known points; then, the ten extracted points are used to perform Hermite interpolation processing on other signals within the range of 100 time steps to ensure the first-order continuity of the result; finally, the signal generated by the interpolation processing replaces the corresponding part of the signal in the original signal, and finally a complete fault source signal is synthesized. In other words, the voltage data corresponding to the second type of time point in the target first time series data is replaced by the corresponding interpolated voltage data to obtain the updated target first time series data. The comparison between the interpolation result (the target first time series data after interpolation replacement) and the original result (the target first time series data before interpolation replacement) is as follows. Fig.12 shown.

[0263] The receiving position is set to be 2000m away, the terminal is connected to a matching resistor, and the terminal reflection is not considered. After calculation, the phase inverse transformation is performed to obtain the final phase result as follows Figure 15-B As shown, the observation results show that when a ground fault occurs, due to the voltage mutation, a high-frequency component will be excited on the line. The transmission of the high-frequency component is affected by the distributed parameters. Compared with traditional industrial software such as MATLAB's Simulink simulation, the distributed parameters can be considered in the calculation process of the present invention, and the simulation results have more detailed information.

[0264] 3) In order to verify the generalization capability of the electromagnetic simulation of the power system distribution network under normal and abnormal working conditions of the present invention, an overvoltage fault is simulated.

[0265] The overvoltage fault excitation signal waveform is as follows: Figure 10-C As shown, the fault signals consider Gaussian pulse, single exponential decay signal, double exponential decay signal and decayed sine signal respectively. The simulation conditions are consistent with the above normal / abnormal state. The calculation results are as follows Figure 15-C As shown, the results show that under the framework of the embodiment of the present application, overvoltage faults can be handled, and the normal / abnormal working conditions of the power system distribution network can be collaboratively simulated, which expands the application scenarios and scope of application.

[0266] The embodiment of the present application uses a grid operation state simulation method implemented by MATLAB, which successfully simulates the electromagnetic characteristics of the three-phase AC power in the power system under normal, single-phase ground fault and overvoltage conditions by combining the transmission line model with the time-domain finite difference algorithm. This method not only reduces the computational complexity caused by the coupling between the three phases, but also optimizes the sudden change processing of the fault voltage through Hermite interpolation processing, thereby improving the simulation accuracy. Compared with traditional simulation tools, this method can consider the influence of distributed parameters more carefully, and the simulation results contain more detailed information, which verifies its stability and generalization ability in the electromagnetic simulation of the normal and abnormal working conditions of the power system distribution network.

[0267] It should be noted that the solution based on industrial software such as Simulink in MATLAB in the related technology can only solve the global voltage and current, and cannot customize the distributed parameter information. In actual scenarios, due to the different line materials of the distribution network, the different heights of the overhead lines, and the different ground characteristics, the distributed parameters are not balanced in the entire line. Therefore, the inability to customize the distributed parameters will cause some information to be missing in the calculation results. The embodiment of the present application uses the time-domain finite difference algorithm based on the transmission line model to customize the distributed parameters on each calculation grid, and the calculation results contain more information and higher accuracy.

[0268] It should also be noted that the relevant calculation method does not incorporate various types of faults into a unified framework for solution, which makes it difficult to simulate the electromagnetic state of the distribution network with long-term multi-state superposition. The embodiment of the present application uses phase mode transformation to unify the calculations in different situations into one architecture, which has a wider range of usage scenarios.

[0269] It should also be noted that in the actual transmission process, since the distribution network is a three-conductor system, there is coupling between the three phases. The relevant technology needs to consider the coupling matrix in the calculation process, which requires a large amount of calculation. The embodiment of the present application decouples the three-phase voltage and current through phase modulus transformation. In the calculation, only the corresponding modulus needs to be considered separately, and there is no need to consider the coupling between the three phases, which effectively reduces the calculation complexity and reduces the space overhead in the calculation process.

[0270] Based on the same inventive concept, the embodiment of the present application also provides a simulation device for the operation status of a power grid, such as Fig.16 As shown, the device 1600 may include an acquisition module 1610 , a conversion module 1620 and a calculation module 1630 .

[0271] An acquisition module 1610 is used to acquire preset operation status information of the distribution network and unit transmission parameter information of each unit transmission path in a preset total transmission path;

[0272] The conversion module 1620 is used to obtain the target first time series data corresponding to the preset operating state information in the corresponding relationship between the preset operating state and the first time series data, and perform phase-mode transformation on the target first time series data to obtain the first mode data, wherein the target first time series data includes the first time series data of each phase voltage signal of the distribution network;

[0273] The calculation module 1630 is used to calculate the second modulus data formed after the first modulus data is transmitted via each unit transmission path in the preset total transmission path based on each unit transmission parameter information;

[0274] The conversion module 1620 is further used to perform analog-to-phase conversion on the second analog data to obtain target second time series data, where the target second time series data includes the second time series data of each phase voltage signal.

[0275] In some embodiments, the acquisition module is used to acquire preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path, which can be specifically used to:

[0276] Display the first page;

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

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

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

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

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

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

[0283] In some embodiments, the conversion module is used to perform phase-modulus transformation on the target first time series data to obtain first modulus data, which can be specifically used for:

[0284] Multiplying the target first time series data by a preset transformation matrix to obtain first modulus data;

[0285]

[0286] Among them, Q is the transformation matrix.

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

[0288] Multiplying the second modulus data by a preset inverse transformation matrix to obtain target second time series data;

[0289]

[0290] Among them, Q -1 is the inverse transformation matrix.

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

[0292] Obtaining information of preset terminal reflection coefficient;

[0293] According to the information of the terminal reflection coefficient, the terminal matching impedance value is calculated;

[0294] Based on each unit transmission parameter information 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 is calculated.

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

[0296] 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;

[0297] When the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is used as the terminal matching impedance value;

[0298]

[0299] Among them, Z L is the terminal matching impedance value, Zc is the equivalent impedance value of the last unit transmission path, and 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 calculation module is further used to calculate the terminal matching impedance value according to the terminal reflection coefficient when the terminal reflection coefficient is not 0;

[0302]

[0303] Where Γ is the terminal reflection coefficient.

[0304] In some embodiments, the calculation module is used 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 each unit transmission parameter information and the terminal matching impedance value, which can be specifically used to:

[0305] Get information about the source internal resistance;

[0306] 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 the first sub-modulus data formed after the first modulus data is transmitted via the first preset unit transmission path;

[0307] 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, calculating the 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;

[0308] When 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 sub-modulus data and the terminal matching impedance value, 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.

[0309] In some embodiments, the calculation module is used to calculate 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, which can be specifically used to:

[0310]

[0311] Among them, Δx is the preset length information of the unit transmission path, Δt is the preset time interval information between the time points of the first timing data, Rs is the source internal resistance information, C is the unit transmission parameter information of the first unit transmission path, N represents the sequence 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 used to calculate 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, which can be specifically used to:

[0313]

[0314] Among them, L0, R0, C0, G0 are the unit transmission parameter information of the N+1th unit transmission path, I N and V N is the Nth submodulus data, I N+1 and V N+1 It is the N+1th sub-modulus data.

[0315] In some embodiments, the calculation module is used 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 sub-modulus data and the terminal matching impedance value when the N+1th unit transmission path is the last unit transmission path, which can be specifically used to:

[0316]

[0317] Among them, C is the 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 used to perform phase-modulus transformation on the target first time series data to obtain the first modulus data, the device further includes an interpolation module and a replacement module:

[0319] The determination module is further used to determine the mutation duration range where the time point corresponding to the voltage data of the mutation type is located when the target first time series data includes voltage data of the mutation type;

[0320] An interpolation module is used to perform interpolation calculation using the voltage data corresponding to the first type of time point within the mutation duration range as known data to obtain interpolation voltage data corresponding to the second type of time point within the mutation duration range;

[0321] A replacement module, used to replace 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 updated target first time series data;

[0322] The conversion module is used to perform phase-modulus transformation on the target first time series data to obtain first modulus data, which can be specifically used for:

[0323] Perform phase modulus transformation on the updated target first time series data to obtain first modulus data.

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

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

[0326] Each module in the simulation device for power grid operation status provided in the embodiment of the present application can realize Figures 1 to 9 The functions of each step of the provided method for simulating the operation status of a power grid and the corresponding technical effects achieved are described briefly and will not be elaborated on here.

[0327] Fig.17 A schematic diagram of the hardware structure of a simulation device for the operation status of a power grid provided in an embodiment of the present application is shown.

[0328] The simulation device in the grid operation state may include a processor 1701 and a memory 1702 storing computer program instructions.

[0329] Specifically, the processor 1701 may include a central processing unit (CPU), or an application specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application.

[0330] Memory 1702 may include a large capacity memory for data or instructions. By way of example and not limitation, memory 1702 may include a hard disk drive (HDD), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a tape, or a universal serial bus (USB) drive or a combination of two or more of these. Where appropriate, memory 1702 may include a removable or non-removable (or fixed) medium. Where appropriate, memory 1702 may be inside or outside a simulation device for the operation state of a power grid. In a particular embodiment, memory 1702 is a non-volatile solid-state memory.

[0331] The memory may include a read-only memory (ROM), a random access memory (RAM), a magnetic disk storage medium device, an optical storage medium device, a flash memory device, an electrical, optical or other physical / tangible memory storage device. Therefore, generally, the memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., memory devices) encoded with software including computer-executable instructions, and when the software is executed (e.g., by one or more processors), it is operable to perform the operations described with reference to the method according to one aspect of the present disclosure.

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

[0333] In one example, the simulation device for the power grid operation state may further include a communication interface 1703 and a bus 1704. Fig.17 As shown, the processor 1701, the memory 1702, and the communication interface 1703 are connected via a bus 1704 and communicate with each other.

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

[0335] The bus 1704 includes hardware, software or both, coupling the components of the simulation device of the power grid operation state to each other. For example and not limitation, the bus may 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 Linear Predictive Coding (LPC) bus, a memory bus, a MicroChannel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (Peripheral Component Interconnect-X, PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local (VESA Local Bus, VLB) bus or other suitable buses or a combination of two or more of these. Where appropriate, the bus 1704 may include one or more buses. Although embodiments of the present application describe and illustrate a particular bus, the present application contemplates any suitable bus or interconnect.

[0336] The device can execute the simulation method of the power grid operation state in the embodiment of the present application based on each unit / component in the simulation device of the power grid operation state, thereby realizing the combination of Figures 1 to 9 A simulation method for describing the operation status of a power grid.

[0337] In addition, in combination with the simulation method of the power grid operation state in the above embodiment, the embodiment of the present application can provide a computer storage medium for implementation. The computer storage medium stores computer program instructions; when the computer program instructions are executed by the processor, any one of the simulation methods of the power grid operation state in the above embodiment is implemented.

[0338] The present application also provides a computer program product. When the instructions in the computer program product are executed by a processor of an electronic device, the electronic device executes each process of an embodiment of a simulation method for implementing any of the above-mentioned power grid operation states.

[0339] It should be clear that the present application is not limited to the specific configuration and processing described above and shown in the figures. For the sake of simplicity, a detailed description of the known method is omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present application is not limited to the specific steps described and shown, and those skilled in the art can make various changes, modifications and additions, or change the order between the steps after understanding the spirit of the present application.

[0340] The functional blocks shown in the above-described block diagram can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, it can be, for example, an electronic circuit, an application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), appropriate firmware, plug-in, function card, etc. When implemented in software, the elements of the present application are programs or code segments used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or communication link by a data signal carried in a carrier. "Machine-readable medium" may include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, read-only memory (Read-Only Memory, ROM), flash memory, erasable read-only memory (Erasable Read Only Memory, EROM), floppy disks, compact disc read-only memory (Compact Disc Read-Only Memory, CD-ROM), optical discs, hard disks, optical fiber media, radio frequency (Radio Frequency, RF) links, etc. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

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

[0342] Aspects of the present disclosure are described above with reference to the flowchart and / or block diagram of the method, device (system) and computer program product according to the embodiment of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device to produce a machine so that these instructions executed by the processor of the computer or other programmable data processing device enable the implementation of the function / action specified in one or more boxes of the flowchart and / or block diagram. Such a processor can be, but is not limited to, a general-purpose processor, a special-purpose processor, a special application processor, or a field programmable logic circuit. It can also be understood that each box in the block diagram and / or flowchart and the combination of boxes in the block diagram and / or flowchart can also be implemented by dedicated hardware that performs a specified function or action, or can be implemented by a combination of dedicated hardware and computer instructions.

[0343] The above are only specific implementation methods of the present application. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, modules and units described above can refer to the corresponding processes in the aforementioned method embodiments, and will not be repeated here. It should be understood that the protection scope of the present application is not limited to this. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed in this application, and these modifications or replacements should be included in the protection scope of this application.

Claims

1. A method for simulating the operation status of a power grid, characterized in that: include: Obtaining preset operating status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path; In the correspondence between the preset operating state and the first time series data, the target first time series data corresponding to the preset operating state information is obtained, and the target first time series data is subjected to phase-modulus transformation to obtain the first modulus data, wherein the target first time series data includes the first time series data of each phase voltage signal of the power distribution network; Based on each of the unit transmission parameter information, calculating second modulus data formed after the first modulus data is transmitted via each of the unit transmission paths in the preset total transmission path; The second modulus data is subjected to an analog-to-phase transformation to obtain target second time series data, wherein the target second time series data includes second time series data of each phase voltage signal.

2. The method according to claim 1, characterized in that The obtaining of preset operation status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path includes: Display the first page; In response to a first input in the first page, obtaining preset operating state information and preset operating environment configuration information of the power distribution network in the first input; According to the preset operating environment configuration information, unit transmission parameter information of each unit transmission path in the preset total transmission path is calculated.

3. The method according to claim 2, characterized in that The preset operating environment configuration information includes terrain information, line distribution information and line material information in the preset total transmission path; the unit transmission parameter information includes the resistance value of the distributed resistor, the inductance value of the distributed inductor, the capacitance value of the distributed capacitor and the conductance value of the distributed conductance in the unit transmission path.

4. The method according to claim 2, characterized in that: The first input also includes first configuration information of the preset total transmission path and second configuration information of the unit transmission path; the method further includes: Acquire first configuration information of the preset total transmission path and second configuration information of the unit transmission path in the first input; The preset total transmission path is determined according to the first configuration information; and the unit transmission path is determined according to the second configuration information.

5. The method according to claim 1, characterized in that The performing phase modulus transformation on the target first time series data to obtain first modulus data includes: Multiplying the target first time series data by a preset transformation matrix to obtain first modulus data; Wherein, Q is the transformation matrix.

6. The method according to claim 5, characterized in that The performing a modulus-to-phase transformation on the second modulus data to obtain target second time series data includes: Multiplying the second modulus data by a preset inverse transformation matrix to obtain target second time series data; Among them, the Q -1 is the inverse transformation matrix.

7. The method according to claim 1, characterized in that The calculating, based on each of the unit transmission parameter information, second modulus data formed after the first modulus data is transmitted via each of the unit transmission paths in the preset total transmission path comprises: Obtaining information of preset terminal reflection coefficient; Calculating the terminal matching impedance value according to the information of the terminal reflection coefficient; Based on each of the unit transmission parameter information and the terminal matching impedance value, second modulus data formed after the first modulus data is transmitted via each of the unit transmission paths in the preset total transmission path is calculated.

8. The method according to claim 7, characterized in that The calculating the terminal matching impedance value according to the information of the terminal reflection coefficient includes: 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; When the terminal reflection coefficient is 0, the equivalent impedance value of the last unit transmission path is used as the terminal matching impedance value; Among them, the Z L is the terminal matching impedance value, Zc is the equivalent impedance value of the last unit transmission path, and L, R, C and G are unit transmission parameter information of the last unit transmission path.

9. The method according to claim 8, characterized in that Also includes: When the terminal reflection coefficient is not 0, calculating the terminal matching impedance value according to the terminal reflection coefficient; Wherein, Γ is the terminal reflection coefficient.

10. The method according to claim 7, characterized in that The calculating, based on each of the unit transmission parameter information and the terminal matching impedance value, second modulus data formed after the first modulus data is transmitted via each of the unit transmission paths in the preset total transmission path comprises: Get information about the source internal resistance; Calculate, based on the source internal resistance information 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 via the first preset unit transmission path; 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, calculating the N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths, wherein 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, 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, the second modulus data formed after the first modulus data is transmitted through each of the unit transmission paths in the preset total transmission path is calculated.

11. The method according to claim 10, characterized in that The calculating, based on the source internal resistance information 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 via the first preset unit transmission path comprises: Among them, the Δx is the preset length information of the unit transmission path, the Δt is the preset time interval information between the time points of the first timing 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 the current, and V represents the voltage.

12. The method according to claim 11, characterized in that The calculating, 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, the N+1th sub-modulus data formed after the first modulus data is transmitted via the first N+1 preset unit transmission paths comprises: Wherein, L0, R0, C0, G0 are unit transmission parameter information of the N+1th unit transmission path, and I N and the V N is the Nth sub-module data, the I N+1 and the V N+1 is the N+1th sub-modulus data.

13. The method according to claim 12, characterized in that In the case where 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 sub-modulus data and the terminal matching impedance value, calculating the second modulus data formed after the first modulus data is transmitted through each of the unit transmission paths in the preset total transmission path, including: Wherein, C is the unit transmission parameter information of the last unit transmission path, and V N+1 is the second modulus data.

14. The method according to claim 1, characterized in that Before performing phase-modulus transformation on the target first time series data to obtain first modulus data, the method further includes: In a case where the target first time series data includes voltage data of a mutation type, determining a mutation duration range where a time point corresponding to the voltage data of the mutation type is located; Performing interpolation calculation using the voltage data corresponding to the first type of time point within the mutation duration range as known data to obtain interpolation voltage data corresponding to the second type of time point within the mutation duration range; Replacing the voltage data corresponding to the second type of time point in the target first time series data with the corresponding interpolated voltage data to obtain updated target first time series data; The performing phase modulus transformation on the target first time series data to obtain first modulus data includes: Perform phase modulus transformation on the updated target first time series data to obtain first modulus data.

15. The method according to any one of claims 1 to 14, characterized in that Also includes: On the second page, waveform information corresponding to the target second time series data is displayed.

16. A simulation device for power grid operation status, characterized in that: include: An acquisition module, used to acquire preset operation status information of the distribution network and unit transmission parameter information of each unit transmission path in the preset total transmission path; A conversion module, used for obtaining target first time series data corresponding to the preset operating state information in the corresponding relationship between the preset operating state and the first time series data, and performing phase-mode conversion on the target first time series data to obtain first mode data, wherein the target first time series data includes first time series data of each phase voltage signal of the power distribution network; A calculation module, configured to calculate, based on each of the unit transmission parameter information, second modulus data formed after the first modulus data is transmitted via each of the unit transmission paths in the preset total transmission path; The conversion module is further used to perform analog-to-phase conversion on the second analog data to obtain target second time series data, wherein the target second time series data includes second time series data of each phase voltage signal.

17. A simulation device for power grid operation status, characterized in that: The device comprises: a processor and a memory storing computer program instructions; when the processor executes the computer program instructions, the method for simulating the operation status of a power grid as claimed in any one of claims 1 to 15 is implemented.

18. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer program instructions, and when the computer program instructions are executed by a processor, the method for simulating the operation status of a power grid according to any one of claims 1 to 15 is implemented.

19. A computer program product, characterized in that When the instructions in the computer program product are executed by a processor of a device, the device is enabled to execute the method for simulating the operation status of a power grid as described in any one of claims 1 to 15.

Images