A method and apparatus for establishing a substation load model

Abstract

The application provides a method and device for establishing a substation load model, wherein the method is used for a substation in a power supply area including multiple grid-connected energy storage devices, and the method collects grid-connected energy storage system data, static load device data, induction motor data and substation sending end data of the substation power supply area, respectively calculates aggregated equivalent parameters of the grid-connected energy storage system, equivalent parameters of static loads in the power supply area, equivalent parameters of dynamic loads and dynamic reactive power compensation parameters, and then establishes a comprehensive load model including the grid-connected energy storage system, which is composed of an equivalent grid-connected energy storage system model, an equivalent static load model, an equivalent dynamic load model and a reactive power compensation model. The method and device can accurately represent the comprehensive load model after the system including multiple grid-connected energy storage devices is added, and overcome the shortcoming that the traditional dynamic load model cannot describe the influence of the grid-connected energy storage on the power grid characteristics.

Description

Technical Field

[0001] The present invention relates to the technical field of power system safety and stability analysis, and more particularly to a method and device for establishing a substation load model. Background Art

[0002] The low inertia and low damping characteristics of dual-high power systems present significant challenges in terms of frequency and voltage stability. In the grid-following (GFL) control method of the converter, phase information at the point of common coupling (PCC) must be measured and acquired through a phase-locked loop (PLL) to maintain synchronization with the grid. However, in weak grids with low physical inertia, this approach can face stability issues. In these situations, the converter is more suitable for a grid-forming (GFM) control method. Grid-forming converters have frequency regulation and voltage control capabilities, enabling them to provide inertia support similar to synchronous generators. Grid-forming energy storage, with its advantages of energy storage and rapid power response, not only provides energy balancing services for the grid but also offers a wider range of stable support and a longer duration. Consequently, grid-forming energy storage devices have been increasingly used in power systems in recent years and have become a hot topic of research and development. Current research on grid-connected energy storage focuses on modeling grid-connected energy storage devices and analyzing the impact of grid-connected energy storage devices on power system stability. This work provides important technical support for the research and analysis of the grid-connected characteristics of centralized grid-connected energy storage systems and the development of new control methods. However, many 330kV (or 220kV) substations are equipped with multiple grid-connected energy storage devices. If each grid-connected energy storage device is modeled in detail in the simulation calculation of the large power grid of the power system, it will cause a dimensionality disaster and fail to meet the requirements in terms of simulation speed and parameter setting. Therefore, it is necessary to establish a comprehensive load model that takes into account the grid-connected energy storage system for substations equipped with a grid-connected energy storage system composed of multiple grid-connected energy storage devices to overcome the shortcomings of traditional dynamic load models that cannot describe the impact of the grid-connected energy storage system on the characteristics of the power grid. Summary of the Invention

[0003] In order to solve the technical problems in the prior art that, when performing large-scale power grid simulation calculations using a substation load model, modeling each grid-type energy storage device will cause a dimensional disaster, resulting in the simulation speed and parameter settings being unable to meet the simulation requirements, and the traditional dynamic load model is unable to describe the impact of the energy storage device on the grid characteristics when the grid is configured, the present invention provides a method and device for establishing a substation load model.

[0004] According to one aspect of the present invention, the present invention provides a method for establishing a substation load model, the method comprising:

[0005] Collect grid-connected energy storage system data, static load equipment data, induction motor data, and substation sending-end data in the substation power distribution area;

[0006] Calculating aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishing an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system;

[0007] Calculating static load equivalent parameters of the power distribution area according to the static load equipment data, and establishing a static load equivalent model based on the static load equivalent parameters;

[0008] Calculating dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters;

[0009] Calculating reactive compensation parameters of the power distribution area based on the grid-type energy storage system data, induction motor data, static load equipment data, and substation sending-end data, and establishing a reactive compensation model based on the reactive compensation parameters;

[0010] A comprehensive load model including a grid-type energy storage system is constructed based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive power compensation model.

[0011] Optionally, the collecting of grid-connected energy storage system data, induction motor data, static load equipment data and substation sending-end data in the power distribution area of ​​the substation includes:

[0012] Collect the rated capacity, maximum active output, maximum reactive output, actual active output, actual reactive output, connection reactance of the N energy storage devices to the common connection point PCC, line reactance and transformer reactance of the N energy storage devices in the grid-type energy storage system, as well as the transient time constant of the model of the N energy storage devices, the time constant of the equivalent virtual synchronization module, virtual damping, the regulator gain of the equivalent virtual excitation control module, the proportional integral selection factor, the time constant of the first voltage regulator and the second voltage regulator, the filter time constant, the synchronous reactance, and the d-axis transient reactance;

[0013] Collect the rated active power, rated reactive power and rated voltage of several static load devices in the substation, as well as the actual reactive output;

[0014] Collect the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy, rotor winding copper loss and actual absorbed reactive power of several induction motors in the substation;

[0015] Collect reactive power at the sending end of the substation.

[0016] Optionally, calculating the aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishing an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system, includes:

[0017] Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system;

[0018] Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula:

[0019]

[0020] X′ C,n =X C,n +X L,n +X T,n

[0021] Where X′ C,n , X C,n , X L,n and X T,n are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, X C,EQ is the equivalent reactance of the grid-type energy storage system;

[0022] According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is:

[0023]

[0024] Where D n is the virtual damping of the equal-valued virtual synchronization module in the model of the n-th grid-type energy storage device, D EQ is the equivalent virtual damping of the grid-type energy storage system;

[0025] The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is:

[0026]

[0027] Where, T j,EQ is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, P max,n is the maximum active output of the nth grid-connected energy storage device, T j,n is the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device;

[0028] According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is:

[0029]

[0030] Where, X d,n and X′ d,n are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, X d,EQ and X′ d,EQ are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively;

[0031] According to the maximum active output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are:

[0032]

[0033]

[0034] Where Q max,n is the maximum active output of the nth energy storage device, T′ d0,n , k n , T 1,n , T 2,n , T R,n and k V,n are the transient time constant of the model of the nth energy storage device, the proportional integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, T′ d0,EQ , k EQ , T 1,eQ , T2,EQ , T R,EQ and k V,EQ are respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system;

[0035] An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system.

[0036] Optionally, calculating static load equivalent parameters of the power distribution area according to the static load equipment data, and establishing a static load equivalent model based on the static load equivalent parameters includes:

[0037] Based on the ZIP model principle of constant impedance, constant current and constant power, the static load equivalent parameters are determined according to the rated active power, rated reactive power and rated voltage of the substation static load equipment;

[0038] An equivalent static load model is established based on the equivalent static load parameters.

[0039] Optionally, calculating the dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters, includes:

[0040] Based on the principle of keeping the total rated active power and reactive power absorbed by the load nodes in the power distribution area unchanged, the total electromagnetic power, maximum electromagnetic power and total kinetic energy unchanged, and the total rotor winding copper loss unchanged, the dynamic load equivalent parameters are calculated based on the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy and rotor winding copper loss of the substation induction motor;

[0041] An equivalent dynamic load model is established based on the dynamic load equivalent parameters.

[0042] Optionally, calculating reactive compensation parameters of the power distribution area according to the grid-type energy storage system data, induction motor data, static load equipment data, and substation sending-end data, and establishing a reactive compensation model based on the reactive compensation parameters includes:

[0043] Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is:

[0044] -Q SC =Q-(Q IM +Q Z +Q I +Q P -Q EQ )

[0045] Where Q SC is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, Q IM Q is the sum of the reactive power actually absorbed by several induction motors in the substation, Z , Q I and Q P are the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices, Q EQ is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system;

[0046] Based on the reactive compensation parameter Q SC Establish reactive power compensation model.

[0047] According to another aspect of the present invention, there is provided a device for establishing a substation load model, the device comprising:

[0048] The data acquisition module is used to collect grid-connected energy storage system data, static load equipment data, induction motor data, and substation sending-end data in the substation power distribution area;

[0049] A first model module is used to calculate the aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establish an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system;

[0050] A second model module is configured to calculate static load equivalent parameters of the power distribution area according to the static load equipment data, and establish a static load equivalent model based on the static load equivalent parameters;

[0051] a third model module, configured to calculate dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establish an equivalent dynamic load model based on the dynamic load equivalent parameters;

[0052] a fourth model module, configured to calculate reactive compensation parameters of the power distribution area based on the grid-type energy storage system data, the induction motor data, the static load equipment data, and the substation sending-end data, and to establish a reactive compensation model based on the reactive compensation parameters;

[0053] The load model module is used to construct a comprehensive load model including the grid-type energy storage system based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive compensation model.

[0054] Optionally, the data acquisition module collects grid-type energy storage system data, induction motor data, static load equipment data and substation sending-end data in the power distribution area of ​​the substation, including:

[0055] Collect the rated capacity, maximum active output, maximum reactive output, actual active output, actual reactive output, connection reactance of the N energy storage devices to the common connection point PCC, line reactance and transformer reactance of the N energy storage devices in the grid-type energy storage system, as well as the transient time constant of the model of the N energy storage devices, the time constant of the equivalent virtual synchronization module, virtual damping, the regulator gain of the equivalent virtual excitation control module, the proportional integral selection factor, the time constant of the first voltage regulator and the second voltage regulator, the filter time constant, the synchronous reactance, and the d-axis transient reactance;

[0056] Collect the rated active power, rated reactive power and rated voltage of several static load devices in the substation, as well as the actual reactive output;

[0057] Collect the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy, rotor winding copper loss and actual absorbed reactive power of several induction motors in the substation;

[0058] Collect reactive power at the sending end of the substation.

[0059] Optionally, the first model module calculates aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system, including:

[0060] Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system;

[0061] Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula:

[0062]

[0063] X′ C,n =X C,n +X L,n +X T,n

[0064] Where X′ C,n , X C,n , X L,n and X T,n are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, X C,EQ is the equivalent reactance of the grid-type energy storage system;

[0065] According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is:

[0066]

[0067] Where D n is the virtual damping of the equal-valued virtual synchronization module in the model of the n-th grid-type energy storage device, D EQ is the equivalent virtual damping of the grid-type energy storage system;

[0068] The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is:

[0069]

[0070] Where, T j,EQ is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, P max,n is the maximum active output of the nth grid-connected energy storage device, T j,n is the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device;

[0071] According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is:

[0072]

[0073] Where, X d,n and X′ d,n are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, X d,EQ and X′ d,EQ are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively;

[0074] According to the maximum active output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are:

[0075]

[0076] Where Q max,n is the maximum active output of the nth energy storage device, T′ d0,n , k n , T 1,n , T 2,n , T R,n and k V,n are the transient time constant of the model of the nth energy storage device, the proportional integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, T′ d0,EQ , k EQ , T 1,EQ , T 2,EQ , T R,EQ and k V,EQ are respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system;

[0077] An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system.

[0078] Optionally, the second model module calculates static load equivalent parameters of the power distribution area according to the static load equipment data, and establishes a static load equivalent model based on the static load equivalent parameters, including:

[0079] Based on the ZIP model principle of constant impedance, constant current and constant power, the static load equivalent parameters are determined according to the rated active power, rated reactive power and rated voltage of the substation static load equipment;

[0080] An equivalent static load model is established based on the equivalent static load parameters.

[0081] Optionally, the third model module calculates dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishes an equivalent dynamic load model based on the dynamic load equivalent parameters, including:

[0082] Based on the principle of keeping the total rated active power and reactive power absorbed by the load nodes in the power distribution area unchanged, the total electromagnetic power, maximum electromagnetic power and total kinetic energy unchanged, and the total rotor winding copper loss unchanged, the dynamic load equivalent parameters are calculated based on the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy and rotor winding copper loss of the substation induction motor;

[0083] An equivalent dynamic load model is established based on the dynamic load equivalent parameters.

[0084] Optionally, the fourth model module calculates reactive compensation parameters of the power distribution area based on the grid-type energy storage system data, induction motor data, static load equipment data, and substation sending-end data, and establishes a reactive compensation model based on the reactive compensation parameters, including:

[0085] Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is:

[0086] -Q SC =Q-(Q IM +QZ +Q I +Q P -Q EQ )

[0087] Where Q SC is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, Q IM Q is the sum of the reactive power actually absorbed by several induction motors in the substation, Z , Q I and Q P are the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices, Q EQ is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system;

[0088] Based on the reactive compensation parameter Q SC Establish reactive power compensation model.

[0089] In the method and device for establishing a substation load model provided by the present invention, the method collects grid-type energy storage system data, static load equipment data, induction motor data and substation sending end data of the substation's power distribution area for a substation including multiple grid-type energy storage devices; calculates the aggregated equivalent parameters of the grid-type energy storage system based on the grid-type energy storage system data, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system; calculates the static load equivalent parameters of the power distribution area based on the static load equipment data, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system; calculates the static load equivalent parameters of the power distribution area based on the static load equivalent parameters, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system; The invention provides a method for establishing a static load equivalent model based on the value parameters; calculating the dynamic load equivalent parameters of the power distribution area based on the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters; calculating the reactive compensation parameters of the power distribution area based on the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data, and establishing a reactive compensation model based on the reactive compensation parameters; constructing a comprehensive load model including the grid-type energy storage system based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive compensation model. The method and device can accurately characterize the comprehensive load model after the system composed of multiple grid-type energy storage devices is added, overcoming the shortcoming that the traditional dynamic load model cannot describe the impact of grid-type energy storage on the characteristics of the power grid. The results of the invention can be widely applied to load modeling software including grid-type energy storage systems, and the generated load model parameters including grid-type energy storage systems can be fully applied to power grid simulation calculations, planning and design, improving the accuracy of power grid simulation calculations, and providing more accurate simulation analysis results for power grid analysis and decision-making. BRIEF DESCRIPTION OF THE DRAWINGS

[0090] A more complete understanding of exemplary embodiments of the present invention may be obtained by referring to the following drawings:

[0091] Figure 1 Flow chart of a method for establishing a substation load model according to a preferred embodiment of the present invention;

[0092] Figure 2 Schematic diagram of the structure of a model of a single grid-type energy storage device according to a preferred embodiment of the present invention;

[0093] Figure 3 A schematic structural diagram of a comprehensive load model considering a grid-type energy storage system according to a preferred embodiment of the present invention;

[0094] Figure 4 A schematic diagram of a 220 kV substation including multiple grid-type energy storage devices according to another preferred embodiment of the present invention;

[0095] Figure 5 A schematic diagram showing a comparison of grid connection point voltage waveforms according to another preferred embodiment of the present invention;

[0096] Figure 6 A schematic diagram showing comparison of active power at a grid connection point according to another preferred embodiment of the present invention;

[0097] Figure 7 A schematic diagram showing a comparison of reactive power at a grid connection point according to another preferred embodiment of the present invention;

[0098] Figure 8 Schematic diagram of the structure of a device for establishing a substation load model according to a preferred embodiment of the present invention;

[0099] Figure 9 Schematic diagram of the structure of an electronic device according to a preferred embodiment of the present invention. DETAILED DESCRIPTION

[0100] Exemplary embodiments of the present invention will now be described with reference to the accompanying drawings. However, the present invention may be embodied in many different forms and is not limited to the embodiments described herein. These embodiments are provided to provide a thorough and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments shown in the accompanying drawings is not intended to limit the present invention. In the accompanying drawings, identical elements are denoted by the same reference numerals.

[0101] Unless otherwise specified, the terms used herein (including technical terms) have the meanings commonly understood by those skilled in the art. In addition, it is understood that terms defined in commonly used dictionaries should be understood to have the same meanings as those in the context of the relevant fields, and should not be understood as idealized or overly formal meanings.

[0102] Exemplary Methods

[0103] Figure 1 FIG. 1 is a flow chart of a method for establishing a substation load model according to a preferred embodiment of the present invention. Figure 1 As shown, the method for establishing a substation load model according to this preferred embodiment starts from step 101 .

[0104] In step 101, grid-connected energy storage system data, static load equipment data, induction motor data and substation sending-end data in the power distribution area of ​​the substation are collected.

[0105] In practical applications, many 330kV or 220kV substations are equipped with multiple grid-type energy storage devices. If each grid-type energy storage device is modeled in detail during large-scale power grid simulation calculations, it will cause a dimensionality curse, making the simulation speed and parameter settings unsatisfactory. To address this issue, this preferred embodiment designs an equivalent grid-type energy storage system model for an energy storage system composed of multiple grid-type energy storage devices. By adding the equivalent grid-type energy storage system model to the substation load model, the characteristics of the large number of grid-type energy storage devices in the power supply and distribution area of ​​the load site are simulated.

[0106] In step 102 , based on the grid-type energy storage system data, aggregated equivalent parameters of the grid-type energy storage system are calculated, and an equivalent grid-type energy storage system model is established based on the aggregated equivalent parameters of the grid-type energy storage system.

[0107] Figure 2 Schematic diagram of the structure of a single grid-type energy storage device according to a preferred embodiment of the present invention. Figure 2 As shown, this preferred embodiment proposes a grid-type energy storage device model that considers the virtual co-frequency control characteristics. The model of a single grid-type energy storage device mainly includes two control modules, one is a virtual co-frequency control module, and the other is a virtual excitation control module. There are 11 model parameters, which are the connection reactance X of the grid-type energy storage device connected to the PCC point, c , transient time constant T d0 ′、Inertia time constant T of virtual synchronization module j , virtual damping D, regulator gain K of virtual excitation control module, proportional integral selection factor K v, first voltage regulator time constant T1, second voltage regulator time constant T2, filter time constant T R , synchronous reactance X d , d-axis transient reactance X d ′.

[0108] The virtual synchronous control module achieves synchronization with the grid by following the traditional synchronous machine rotor motion equations, thus introducing the risk of power angle instability associated with traditional synchronous machines. The differences lie in the inertia time constant Tj and the damping coefficient D. While these are fixed physical parameters for synchronous machines, they are manually adjustable control parameters for virtual synchronous control, potentially resulting in differences in characteristics from those of real synchronous machines. For a virtual synchronous controller, its transfer function can be transformed into a first-order inertia and integral series link. Since the damping coefficient D represents the total damping torque of the virtual synchronous machine, its value determines the key characteristics of the network-based controller. When Tj / D is extremely small, the virtual synchronous control characteristic becomes a droop control characteristic. When Tj / D is large, the virtual synchronous control characteristic maintains a characteristic close to that of a synchronous machine. Because the support capacity of the grid-type energy storage device for the system frequency is mainly determined by the maximum active power of the grid-type energy storage device, therefore, in this preferred embodiment, when equating the grid-type energy storage system composed of multiple grid-type energy storage devices, the virtual synchronization parameters of the equivalent grid-type energy storage system model should be calculated by weighted average using the proportion of the maximum active power of each grid-type energy storage device relative to the total maximum active power of all grid-type energy storage devices in the entire load site. Similarly, the virtual excitation control module determines the voltage support capacity of the grid-type energy storage. The impact of the grid-type energy storage on voltage stability is mainly reflected by the amplitude of the internal potential. During the transient process, voltage support can be divided into two stages: 1) Maintaining a stable internal potential during the fault period and passively providing voltage support to the system; 2) Under the action of the virtual excitation controller, the internal potential actively rises, the voltage support capacity is enhanced and continues until the system voltage returns to normal. Because the supporting capacity of the grid-type energy storage for the system voltage is mainly determined by the maximum reactive power of the grid-type energy storage, when the grid-type energy storage system composed of multiple grid-type energy storage devices is equal, the virtual excitation parameters of the equalized grid-type energy storage model should be calculated by weighted average using the proportion of the maximum reactive power of each grid-type energy storage relative to the total maximum reactive power of all grid-type energy storages in the entire 220kV load site.

[0109] Preferably, the calculating of aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishing an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system, includes:

[0110] Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system;

[0111] Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula:

[0112]

[0113] X′ C,n =X C,n +X L,n +X T,n

[0114] Where X′ C,n , X C,n , X L,n and X T,n are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, X C,EQ is the equivalent reactance of the grid-type energy storage system;

[0115] According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is:

[0116]

[0117] Where D n is the virtual damping of the equal-valued virtual synchronization module in the model of the n-th grid-type energy storage device, D EQ is the equivalent virtual damping of the grid-type energy storage system;

[0118] The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is:

[0119]

[0120] Where, T j,EQ is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, P max,n is the maximum active output of the nth grid-connected energy storage device, T j,n is the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device;

[0121] According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is:

[0122]

[0123] Where, X d,n and X′ d,n are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, X d,EQ and X′ d,EQ are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively;

[0124] According to the maximum active output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are:

[0125]

[0126] Where Q max,n is the maximum active output of the nth energy storage device, T′ d0,n , k n , T 1,n , T 2,n , T R,n and k V,n are the transient time constant of the model of the nth energy storage device, the proportional integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, T′ d0,EQ , k EQ , T 1,EQ , T 2,EQ , T R,EQ and k V,EQ are respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system;

[0127] An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system.

[0128] In step 103, static load equivalent parameters of the power distribution area are calculated according to the static load equipment data, and a static load equivalent model is established based on the static load equivalent parameters.

[0129] Preferably, calculating the static load equivalent parameters of the power distribution area according to the static load equipment data, and establishing a static load equivalent model based on the static load equivalent parameters, includes:

[0130] Based on the ZIP model principle of constant impedance, constant current and constant power, the static load equivalent parameters are determined according to the rated active power, rated reactive power and rated voltage of the substation static load equipment;

[0131] An equivalent static load model is established based on the equivalent static load parameters.

[0132] Specifically, the equivalent static load model described in this preferred embodiment describes the relationship between load power and voltage as a polynomial load model in the form of a polynomial equation. The general form of the model is as follows:

[0133]

[0134] The above load model is called a ZIP model because it includes constant impedance Z, constant current I, and constant power P. This model is used to describe specific load devices or load elements, that is, load devices or load elements that have any one of the following characteristics: constant impedance Z, constant current I, and constant power P. In the above formula, the polynomial active power load model coefficients are a, b, and c, and the reactive power load model coefficients are α, β, and γ. Vo represents the rated voltage of the load, and Po and Qo represent the rated active power and reactive power of the load at the rated voltage Vo, respectively. If this model is used to describe the comprehensive load of a bus, Vo, Po, and Qo are usually used to represent the values ​​under the system's initial operating conditions.

[0135] The equivalent values ​​for static loads are mainly the equivalent values ​​of coefficients Po, a, b, c and Qo, α, β, γ. The equivalent values ​​for polynomial load models are based on the sensitivity of load power to load terminal voltage, that is,

[0136]

[0137] P1, P2…Pn and Q1, Q2…Qn are the active power and reactive power of each static load, and the corresponding polynomial load model coefficients are Po1…Pon, a1…an, b1…bn, c1…cn and Qo1…Qon, α1…αn, β1…βn, γ1…γn. When V=Vo, we have:

[0138]

[0139] According to the above formula, by substituting the collected data values, the equivalent static load parameters can be solved, and then an equivalent static load model can be established based on the solved equivalent static load parameters.

[0140] In step 104 , dynamic load equivalent parameters of the power distribution area are calculated according to the induction motor data, and an equivalent dynamic load model is established based on the dynamic load equivalent parameters.

[0141] Preferably, calculating the dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters, includes:

[0142] Based on the principle of keeping the total rated active power and reactive power absorbed by the load nodes in the power distribution area unchanged, the total electromagnetic power, maximum electromagnetic power and total kinetic energy unchanged, and the total rotor winding copper loss unchanged, the dynamic load equivalent parameters are calculated based on the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy and rotor winding copper loss of the substation induction motor;

[0143] An equivalent dynamic load model is established based on the dynamic load equivalent parameters.

[0144] Specifically, in this preferred embodiment, the calculation of dynamic load equivalent parameters of the power distribution area based on the collected induction motor data includes:

[0145] The formulas for calculating the total stator winding copper loss ∑Pcu1, the rated slip Sn of the equivalent motor, and the equivalent inertia time constant H are as follows:

[0146]

[0147] Where, Pemn and Pn are the rated electromagnetic power and rated active power of the n-th motor, Pcu2 and Eenergy are the rotor winding copper loss and kinetic energy of the n-th motor, respectively.

[0148] Among them, ∑Pemn-∑Pcu2 is the rated mechanical power output by the equivalent motor, which remains unchanged;

[0149] Calculate the electrical parameters of the equivalent motor model, including stator resistance Rs, stator leakage reactance Xs, rotor resistance Rr, rotor leakage reactance Xr, and magnetizing reactance Xm. Assuming the rated phase voltage is Un, the calculation process of the electrical parameters is as follows:

[0150] Let Pemt_max = ∑Pemn_max; Pent_max is the total maximum electromagnetic power, Pemn_max is the maximum electromagnetic power of the nth motor;

[0151] Assume the total stator phase current is I n , which is calculated as follows:

[0152]

[0153] Where Pn, Qn and Un are the rated active power, rated reactive power and rated phase voltage of the nth motor respectively;

[0154]

[0155] Among them, P cul is the total stator winding copper loss of the equivalent motor, In is the total stator phase current;

[0156] Calculate the equivalent impedance Zdeq of the equivalent motor, the calculation formula is:

[0157]

[0158] Rdeq=real(Zdeq)

[0159] Xdeq=imag(Zdeq)

[0160] Among them, R dep is the equivalent resistance of the equivalent motor, X dep is the equivalent reactance of the equivalent motor, j is the imaginary unit of the complex number;

[0161] Calculate the stator leakage reactance Xs and rotor leakage reactance Xr using the following formula:

[0162]

[0163] Xr=Xs

[0164] In this algorithm, it is always assumed that Xr = Xs, and Xs and Xr calculated according to the above formula are necessarily too small, because the maximum electromagnetic power calculated according to the simplified maximum electromagnetic power formula is larger than the actual maximum electromagnetic power, so Xs and Xr need to be corrected through an iterative method;

[0165] According to the calculated stator resistance Rs, stator leakage reactance Xs, rotor leakage reactance Xr and equivalent impedance Zdeq = Rdeq + jXdeq, the rotor resistance Rr and magnetizing reactance Xm are calculated so that Pem = ∑Pemn holds true, where K r K is the difference between the equivalent resistance of the equivalent motor and the stator resistance. x is the difference between the equivalent reactance and stator reactance of an equivalent motor:

[0166] Kr=Rdeq-Rs

[0167] Kx=Xdeq-Xs

[0168]

[0169] Based on the obtained Rs, Xs, Rt, Xr and Xm, the maximum electromagnetic power is recalculated according to the simplified formula:

[0170]

[0171] Calculate the actual maximum electromagnetic power under the new parameters based on the Thevenin equivalent circuit:

[0172] Thevenin equivalent impedance is:

[0173]

[0174] Rdp=real(Zdp)

[0175] Xdp=imag(Zdp)

[0176] Among them, R dp is the Thevenin equivalent resistance, X dp is the Thevenin equivalent reactance, Z dp is the Thevenin equivalent impedance; the condition for generating maximum electromagnetic power is:

[0177]

[0178] Among them, S m is the critical slip, R pm is the Thevenin equivalent impedance value corresponding to the maximum electromagnetic power; the open-circuit voltage of the Thevenin equivalent circuit is:

[0179]

[0180] Recalculate the actual maximum electromagnetic torque corresponding to the new parameters according to the following formula:

[0181]

[0182] Calculate the ratio of Pemt_maxi to Pem_maxi and correct Pemt_max:

[0183]

[0184] Pemt_max=kmaxiPem_max

[0185] Compare the difference between Pem_maxi and Pem_max:

[0186] ErrPem_max=|Pem_max-[em_maxi|

[0187] If ErrPem_max ≥ 1.0e -5 , then recalculate the stator leakage reactance Xs and rotor leakage reactance Xr, otherwise end the calculation.

[0188] After obtaining the total stator winding copper loss, rated slip, equivalent inertia time constant, and electrical parameters of the above-mentioned equivalent motor, an equivalent dynamic load model is generated according to the above-mentioned parameters.

[0189] In step 105, reactive compensation parameters of the power distribution area are calculated based on the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data, and a reactive compensation model is established based on the reactive compensation parameters.

[0190] Preferably, the reactive compensation parameters of the power distribution area are calculated based on the grid-type energy storage system data, induction motor data, static load equipment data and substation sending-end data, and a reactive compensation model is established based on the reactive compensation parameters, including:

[0191] Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is:

[0192] -Q SC =Q-(Q IM +Q Z +Q I +Q P -Q EQ )

[0193] Where Q SC is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, Q IM Q is the sum of the reactive power actually absorbed by several induction motors in the substation, Z , Q I and Q Pare the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices, Q EQ is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system;

[0194] Based on the reactive compensation parameter Q SC Establish reactive power compensation model.

[0195] In step 106, a comprehensive load model including the grid-type energy storage system is constructed based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive power compensation model.

[0196] Figure 3 Schematic diagram of the load model of the grid-type energy storage system according to the preferred embodiment of the present invention. Figure 3 As shown, the load model in this preferred embodiment is divided into four main parts: an equivalent grid-type energy storage model, an equivalent static load model, an equivalent dynamic load model, and a reactive compensation model. The equivalent static load includes a constant impedance load, a constant current load, and a constant power load.

[0197] Figure 4 FIG. 1 is a schematic diagram of a 220 kV substation including multiple grid-type energy storage devices according to another preferred embodiment of the present invention. Figure 4 As shown, in this preferred embodiment, the 220kV substation obtained through investigation is taken as an example for analysis and explanation. The 220kV substation is equipped with 8 grid-type energy storage devices G1-G8, and their operating characteristic parameters are shown in Table 1 and Table 1.

[0198] Table 1

[0199]

[0200]

[0201] Table 1

[0202] K Kv <![CDATA[T1(s)]]> <![CDATA[T2(s)]]> <![CDATA[T R (s)]]> <![CDATA[X d (pu)]]> <![CDATA[X d ′(pu)]]> Network type G1 7 0.9 2.0 0.5 0.03 2.41 0.259 Network type G2 12 0.5 2.0 0.5 0.01 3.5 0.45 Network type G3 10 0.8 2.0 0.5 0.02 2.5 0.31 Network type G4 5 0.7 2.0 0.5 0.05 2.1 0.2 Network type G5 6 0.6 2.0 0.5 0.04 2.3 0.21 Network type G6 7 0.9 2.0 0.5 0.03 2.41 0.259 Network type G7 12 0.5 2.0 0.5 0.01 3.5 0.45 Network type G8 5 0.7 2.0 0.5 0.05 2.1 0.2

[0203] The aggregated equivalent parameters obtained by the method for constructing an equivalent grid-type energy storage system model described in the present invention are shown in Table 2 and Table 2.

[0204] Table 2

[0205]

[0206] Table 2

[0207] K Kv <![CDATA[T1(s)]]> <![CDATA[T2(s)]]> <![CDATA[T R (s)]]> <![CDATA[X d (pu)]]> <![CDATA[X d ′(pu)]]> 8.5 0.69 2.0 0.5 0.027 0.31 0.033

[0208] After obtaining the static load equivalent parameters, dynamic load equivalent parameters and reactive compensation parameters using the method described in this preferred embodiment, an equivalent model considering the grid-type energy storage system is established to carry out power system simulation, and the results of transient simulation are compared with those of a detailed distribution network model using a model of 8 grid-type energy storage devices. Figure 5 This is a schematic diagram showing a comparison of grid connection point voltage waveforms according to another preferred embodiment of the present invention. Figure 6 This is a schematic diagram of active power comparison of grid connection points according to another preferred embodiment of the present invention. Figure 7 FIG. 1 is a schematic diagram showing a comparison of reactive power at a grid connection point according to another preferred embodiment of the present invention. Figure 5 , Figure 6 and Figure 7 As shown, the voltage waveforms, active power waveforms, and reactive power waveforms of the comprehensive load model and the detailed distribution network model generated by the method for establishing a load model described in the present invention are highly consistent, indicating that the model considering the grid-type energy storage system established based on the method described in the present invention can well reflect the load characteristics of the active distribution network containing multiple grid-type energy storage devices.

[0209] The method for establishing a load model described in this preferred embodiment is for a substation whose power distribution area includes a grid-type energy storage system composed of multiple grid-type energy storage devices, and collects grid-type energy storage system data, static load equipment data, induction motor data and substation sending-end data of the substation power distribution area; calculates the aggregated equivalent parameters of the grid-type energy storage system based on the grid-type energy storage system data, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system; calculates the static load equivalent parameters of the power distribution area based on the static load equipment data, and establishes an equivalent grid-type energy storage system model based on the static load Equivalent parameters are used to establish a static load equivalent model; based on the induction motor data, the dynamic load equivalent parameters of the power distribution area are calculated, and an equivalent dynamic load model is established based on the dynamic load equivalent parameters; based on the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data, the reactive compensation parameters of the power distribution area are calculated, and a reactive compensation model is established based on the reactive compensation parameters; based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive compensation model, a comprehensive load model including the grid-type energy storage system is constructed. The method can accurately characterize the comprehensive load model after the system composed of multiple grid-type energy storage devices is added, overcoming the shortcoming that the traditional dynamic load model cannot describe the impact of grid-type energy storage on the characteristics of the power grid. The results of this preferred embodiment can be widely applied to load modeling software including grid-type energy storage systems, and the generated load model parameters including grid-type energy storage systems can be fully applied to power grid simulation calculations, planning and design, improving the accuracy of power grid simulation calculations, and providing more accurate simulation analysis results for power grid analysis and decision-making.

[0210] Exemplary devices

[0211] Figure 8 FIG. 1 is a schematic diagram of a structure of a device for establishing a substation load model according to a preferred embodiment of the present invention. Figure 8 As shown, the device 800 for establishing a substation load model in this preferred embodiment includes:

[0212] Data acquisition module 801, used to collect grid-type energy storage system data, static load equipment data, induction motor data and substation sending end data in the power distribution area of ​​the substation;

[0213] A first model module 802 is configured to calculate aggregate equivalent parameters of the grid-type energy storage system based on the grid-type energy storage system data, and establish an equivalent grid-type energy storage system model based on the aggregate equivalent parameters of the grid-type energy storage system;

[0214] The second model module 803 is configured to calculate static load equivalent parameters of the power distribution area according to the static load equipment data, and establish a static load equivalent model based on the static load equivalent parameters;

[0215] A third model module 804 is configured to calculate dynamic load equivalent parameters of the power distribution area based on the induction motor data, and establish an equivalent dynamic load model based on the dynamic load equivalent parameters;

[0216] A fourth model module 805 is configured to calculate reactive compensation parameters for the power distribution area based on the grid-type energy storage system data, the induction motor data, the static load equipment data, and the substation sending-end data, and establish a reactive compensation model based on the reactive compensation parameters.

[0217] The load model module 806 is used to construct a comprehensive load model including the grid-type energy storage system based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive compensation model.

[0218] Preferably, the data acquisition module 801 collects grid-type energy storage system data, induction motor data, static load equipment data and substation sending-end data in the power distribution area of ​​the substation, including:

[0219] Collect the rated capacity, maximum active output, maximum reactive output, actual active output, actual reactive output, connection reactance of the N energy storage devices to the common connection point PCC, line reactance and transformer reactance of the N energy storage devices in the grid-type energy storage system, as well as the transient time constant of the model of the N energy storage devices, the time constant of the equivalent virtual synchronization module, virtual damping, the regulator gain of the equivalent virtual excitation control module, the proportional integral selection factor, the time constant of the first voltage regulator and the second voltage regulator, the filter time constant, the synchronous reactance, and the d-axis transient reactance;

[0220] Collect the rated active power, rated reactive power and rated voltage of several static load devices in the substation, as well as the actual reactive output;

[0221] Collect the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy, rotor winding copper loss and actual absorbed reactive power of several induction motors in the substation;

[0222] Collect reactive power at the sending end of the substation.

[0223] Preferably, the first model module 802 calculates the aggregated equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishes an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system, including:

[0224] Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system;

[0225] Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula:

[0226]

[0227] X′ C,n =X C,n +X L,n +X T,n

[0228] Where X′ C,n , X C,n , X L,n and X T,n are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, X C,EQ is the equivalent reactance of the grid-type energy storage system;

[0229] According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is:

[0230]

[0231] Where D n is the virtual damping of the equal-valued virtual synchronization module in the model of the n-th grid-type energy storage device, D EQ is the equivalent virtual damping of the grid-type energy storage system;

[0232] The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is:

[0233]

[0234] Where, T j,EQ is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, P max,n is the maximum active output of the nth grid-connected energy storage device, T j,nis the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device;

[0235] According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is:

[0236]

[0237] Where, X d,n and X′ d,n are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, X d,EQ and X′ d,EQ are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively;

[0238] According to the maximum active output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are:

[0239]

[0240]

[0241] Where Q max,n is the maximum active output of the nth energy storage device, T′ d0,n , k n , T 1,n , T 2,n , T R,n and k V,n are the transient time constant of the model of the nth energy storage device, the proportional integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, T′ d0,EQ , k EQ , T 1,EQ , T 2,EQ , T R,EQ and k V,EQare respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system;

[0242] An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system.

[0243] Preferably, the second model module 803 calculates the static load equivalent parameters of the power distribution area according to the static load equipment data, and establishes a static load equivalent model based on the static load equivalent parameters, including:

[0244] Based on the ZIP model principle of constant impedance, constant current and constant power, the static load equivalent parameters are determined according to the rated active power, rated reactive power and rated voltage of the substation static load equipment;

[0245] An equivalent static load model is established based on the equivalent static load parameters.

[0246] Preferably, the third model module 804 calculates the dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishes an equivalent dynamic load model based on the dynamic load equivalent parameters, including:

[0247] Based on the principle of keeping the total rated active power and reactive power absorbed by the load nodes in the power distribution area unchanged, the total electromagnetic power, maximum electromagnetic power and total kinetic energy unchanged, and the total rotor winding copper loss unchanged, the dynamic load equivalent parameters are calculated based on the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy and rotor winding copper loss of the substation induction motor;

[0248] An equivalent dynamic load model is established based on the dynamic load equivalent parameters.

[0249] Preferably, the fourth model module 805 calculates reactive compensation parameters of the power distribution area based on the grid-type energy storage system data, induction motor data, static load equipment data, and substation sending-end data, and establishes a reactive compensation model based on the reactive compensation parameters, including:

[0250] Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is:

[0251] -Q SC =Q-(Q IM +Q Z +Q I +Q P -Q EQ )

[0252] Where Q SC is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, Q IM Q is the sum of the reactive power actually absorbed by several induction motors in the substation, Z , Q I and Q P are the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices, Q EQ is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system;

[0253] Based on the reactive compensation parameter Q SC Establish reactive power compensation model.

[0254] The device for establishing a substation load model described in this preferred embodiment is aimed at a substation including a grid-type energy storage system composed of multiple grid-type energy storage devices. By collecting grid-type energy storage system data, static load equipment data, induction motor data and substation sending-end data in the substation power distribution area, the aggregation equivalent parameters of the grid-type energy storage system, the static load equivalent parameters, dynamic load equivalent parameters and reactive compensation parameters of the power distribution area are calculated respectively. The steps of establishing a comprehensive load model including an equivalent grid-type energy storage model are the same as the steps taken by the method for establishing a substation load model described in the present invention, and the technical effects achieved are also the same, which will not be repeated here.

[0255] Exemplary electronic devices

[0256] Figure 9 1 is a schematic diagram of the structure of an electronic device according to a preferred embodiment of the present invention. The electronic device can be either or both of the first device and the second device, or a standalone device independent of them, and the standalone device can communicate with the first device and the second device to receive collected input signals from them. Figure 9 FIG2 is a block diagram of an electronic device according to an embodiment of the present disclosure. Figure 9 As shown, the electronic device includes one or more processors 901 and a memory 902 .

[0257] The processor 901 may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions.

[0258] The memory 902 may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, and the processor 901 may run the program instructions to implement the energy consumption anomaly diagnosis method based on the enterprise energy consumption space of the various embodiments disclosed above and / or other desired functions. In one example, the electronic device may further include: an input device 903 and an output device 904, and these components are interconnected via a bus system and / or other forms of connection mechanisms (not shown).

[0259] In addition, the input device 903 may also include, for example, a keyboard, a mouse, and the like.

[0260] The output device 904 can output various information to the outside, and can include, for example, a display, a speaker, a printer, a communication network and a remote output device connected thereto.

[0261] Of course, to simplify, Figure 9 Only some of the components related to the present disclosure in the electronic device are shown, and components such as buses, input / output interfaces, etc. are omitted. In addition, the electronic device may further include any other appropriate components according to specific application scenarios.

[0262] Exemplary computer program products and computer-readable storage media

[0263] In addition to the above-mentioned methods and devices, an embodiment of the present disclosure may also be a computer program product, which includes computer program instructions, which, when executed by a processor, enable the processor to execute the steps in the method of establishing a substation load model according to various embodiments of the present disclosure described in the above-mentioned "Exemplary Method" section of this specification.

[0264] The computer program product may be written in any combination of one or more programming languages ​​to implement the operations of the disclosed embodiments, including object-oriented programming languages ​​such as Java, C++, and conventional procedural programming languages ​​such as C or similar programming languages. The program code may be executed entirely on the user's computing device, partially on the user's computing device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0265] In addition, an embodiment of the present disclosure may also be a computer-readable storage medium having computer program instructions stored thereon, which, when executed by a processor, enable the processor to execute the steps in the method for establishing a substation load model according to various embodiments of the present disclosure described in the above "Exemplary Method" section of this specification.

[0266] The computer-readable storage medium can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can, for example, include but is not limited to a system, device or component of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination thereof. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof.

[0267] The basic principles of the present disclosure have been described above in conjunction with specific embodiments. However, it should be noted that the advantages, strengths, and effects mentioned in this disclosure are merely illustrative and not restrictive, and should not be construed as necessarily possessed by each embodiment of the present disclosure. Furthermore, the specific details disclosed above are provided for illustrative purposes and to facilitate understanding, rather than as limitations. These details do not limit the present disclosure to necessarily being implemented using these specific details.

[0268] Each embodiment in this specification is described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. References to the same or similar parts between the various embodiments are sufficient. For system embodiments, since they largely correspond to method embodiments, their description is relatively simple. For relevant parts, references to the description of the method embodiments are sufficient.

[0269] The block diagrams of the devices, devices, equipment, and systems involved in this disclosure are merely illustrative examples and are not intended to require or imply that they must be connected, arranged, or configured in the manner shown in the block diagrams. As will be appreciated by those skilled in the art, these devices, devices, equipment, and systems can be connected, arranged, or configured in any manner. Words such as "include," "comprise," "have," and the like are open-ended words, meaning "including but not limited to," and can be used interchangeably therewith. The words "or" and "and" used herein refer to the words "and / or" and can be used interchangeably therewith, unless the context clearly indicates otherwise. The word "such as" used herein refers to the phrase "such as but not limited to," and can be used interchangeably therewith.

[0270] The apparatus and method of the present disclosure may be implemented in many ways. For example, the apparatus and method of the present disclosure may be implemented by software, hardware, firmware, or any combination of software, hardware, and firmware. The above order of steps for the method is for illustration only, and the steps of the method of the present disclosure are not limited to the order specifically described above unless otherwise specified. In addition, in some embodiments, the present disclosure may also be implemented as programs recorded in a recording medium, which include machine-readable instructions for implementing the method according to the present disclosure. Therefore, the present disclosure also covers recording media that store programs for executing the method according to the present disclosure.

[0271] It should also be noted that, in the apparatus, equipment and method of the present disclosure, each component or each step can be decomposed and / or recombined. These decompositions and / or recombinations should be regarded as equivalent schemes of the present disclosure. The above description of the disclosed aspects is provided to enable any technician in this field to make or use the present disclosure. Various modifications to these aspects will be very obvious to those skilled in the art, and the general principles defined herein can be applied to other aspects without departing from the scope of the present disclosure. Therefore, the present disclosure is not intended to be limited to the aspects shown here, but to the widest range consistent with the principles and novel features disclosed herein.

[0272] The above description has been provided for the purpose of illustration and description. In addition, this description is not intended to limit the embodiments of the present disclosure to the forms disclosed herein. Although a number of example aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations thereof.

Claims

1. A method for establishing a substation load model, characterized in that: The method comprises: Collect grid-connected energy storage system data, static load equipment data, induction motor data, and substation sending-end data in the substation power distribution area, including: Collect the rated capacity, maximum active output, maximum reactive output, actual active output, actual reactive output, connection reactance of the N energy storage devices to the common connection point PCC, line reactance and transformer reactance of the N energy storage devices in the grid-type energy storage system, as well as the transient time constant of the model of the N energy storage devices, the time constant of the equivalent virtual synchronization module, virtual damping, the regulator gain of the equivalent virtual excitation control module, the proportional integral selection factor, the time constant of the first voltage regulator and the second voltage regulator, the filter time constant, the synchronous reactance, and the d-axis transient reactance; Collect the rated active power, rated reactive power and rated voltage of several static load devices in the substation, as well as the actual reactive output; Collect the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy, rotor winding copper loss and actual absorbed reactive power of several induction motors in the substation; Collect reactive power at the sending end of the substation; Calculating aggregate equivalent parameters of the grid-type energy storage system according to the grid-type energy storage system data, and establishing an equivalent grid-type energy storage system model based on the aggregate equivalent parameters of the grid-type energy storage system, including: Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system; Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula: Where, , , and are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, is the equivalent reactance of the grid-type energy storage system; According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is: Where, is the virtual damping of the equal-valued virtual synchronization module in the model of the nth grid-type energy storage device, is the equivalent virtual damping of the grid-type energy storage system; The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is: Where, is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, is the maximum active output of the nth grid-connected energy storage device, is the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device; According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is: Where, and are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, and are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively; According to the maximum reactive output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are: Where, is the maximum reactive power output of the nth energy storage device, , , , , and are the transient time constant of the model of the nth energy storage device, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, , , , , and are respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system; An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system; Calculating static load equivalent parameters of the power distribution area according to the static load equipment data, and establishing an equivalent static load model based on the static load equivalent parameters; Calculating dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters; Calculating reactive power compensation parameters of the power distribution area based on the grid-type energy storage system data, induction motor data, static load equipment data, and substation sending-end data, and establishing a reactive power compensation model based on the reactive power compensation parameters, including: Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is: Where, is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, is the sum of the reactive powers actually absorbed by several induction motors in the substation. , and They are the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices. is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system; Based on the reactive compensation parameters Establish reactive power compensation model; A comprehensive load model including a grid-type energy storage system is constructed based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive power compensation model.

2. The method according to claim 1, characterized in that Calculating the static load equivalent parameters of the power distribution area according to the static load equipment data, and establishing an equivalent static load model based on the static load equivalent parameters, includes: Based on the ZIP model principle of constant impedance, constant current and constant power, the static load equivalent parameters are determined according to the rated active power, rated reactive power and rated voltage of the substation static load equipment; An equivalent static load model is established based on the equivalent static load parameters.

3. The method according to claim 1, characterized in that Calculating the dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establishing an equivalent dynamic load model based on the dynamic load equivalent parameters, includes: Based on the principle of keeping the total rated active power and reactive power absorbed by the load nodes in the power distribution area unchanged, the total electromagnetic power, maximum electromagnetic power and total kinetic energy unchanged, and the total rotor winding copper loss unchanged, the dynamic load equivalent parameters are calculated based on the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy and rotor winding copper loss of the substation induction motor; An equivalent dynamic load model is established based on the dynamic load equivalent parameters.

4. A device for establishing a substation load model, characterized in that: The device comprises: The data acquisition module is used to collect grid-connected energy storage system data, static load equipment data, induction motor data, and substation sending-end data in the substation power distribution area, including: Collect the rated capacity, maximum active output, maximum reactive output, actual active output, actual reactive output, connection reactance of the N energy storage devices to the common connection point PCC, line reactance and transformer reactance of the N energy storage devices in the grid-type energy storage system, as well as the transient time constant of the model of the N energy storage devices, the time constant of the equivalent virtual synchronization module, virtual damping, the regulator gain of the equivalent virtual excitation control module, the proportional integral selection factor, the time constant of the first voltage regulator and the second voltage regulator, the filter time constant, the synchronous reactance, and the d-axis transient reactance; Collect the rated active power, rated reactive power and rated voltage of several static load devices in the substation, as well as the actual reactive output; Collect the rated active power, rated reactive power, rated phase voltage, rated electromagnetic power, maximum electromagnetic power, kinetic energy, rotor winding copper loss and actual absorbed reactive power of several induction motors in the substation; Collect reactive power at the sending end of the substation; The first model module is configured to calculate, based on the grid-type energy storage system data, aggregated equivalent parameters of the grid-type energy storage system, and establish an equivalent grid-type energy storage system model based on the aggregated equivalent parameters of the grid-type energy storage system, including: Summing the rated capacity, maximum active output, maximum reactive output, actual active output, and actual reactive output of N energy storage devices in the grid-type energy storage system, respectively, to determine the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, and equivalent actual reactive output of the grid-type energy storage system; Based on the connection reactance, line reactance and transformer reactance of N energy storage devices to the common connection point PCC in the grid-type energy storage system, the equivalent reactance is determined by the following calculation formula: Where, , , and are the equivalent reactance, connection reactance, line reactance and transformer reactance of the nth grid-connected energy storage device connected to the PCC, is the equivalent reactance of the grid-type energy storage system; According to the virtual damping of the equivalent virtual synchronization module in the model of N energy storage devices in the grid-type energy storage system, the equivalent virtual damping of the grid-type energy storage system is determined, and the calculation formula is: Where, is the virtual damping of the equal-valued virtual synchronization module in the model of the nth grid-type energy storage device, is the equivalent virtual damping of the grid-type energy storage system; The time constant of the equivalent virtual synchronization module of the grid-type energy storage system is determined based on the maximum active output of the N energy storage devices in the grid-type energy storage system and the time constant of the equivalent virtual synchronization module of the corresponding device model. The calculation formula is: Where, is the time constant of the equivalent virtual synchronization module of the grid-type energy storage system, is the maximum active output of the nth grid-connected energy storage device, is the time constant of the equivalent virtual synchronization module of the model of the nth grid-connected energy storage device; According to the synchronous reactance and d-axis transient reactance of the models of the N energy storage devices in the grid-type energy storage system, the equivalent synchronous reactance and equivalent d-axis transient reactance of the grid-type energy storage system are determined respectively, and the calculation formula is: Where, and are the equivalent synchronous reactance and equivalent d-axis transient reactance of the model of the nth grid-connected energy storage device, and are the equivalent synchronous reactance and the equivalent d-axis transient reactance of the grid-type energy storage system respectively; According to the maximum reactive output of N energy storage devices in the grid-type energy storage system and the transient time constant of the corresponding device model, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, the equivalent transient time constant, the equivalent proportional-integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the regulator gain of the equivalent virtual excitation control module of the grid-type energy storage system are determined, and the calculation formulas are: Where, is the maximum reactive power output of the nth energy storage device, , , , , and are the transient time constant of the model of the nth energy storage device, the proportional-integral selection factor, the time constants of the first and second voltage regulators, the filter time constant and the regulator gain of the equivalent virtual excitation control module, , , , , and are respectively the equivalent transient time constant, the equivalent proportional integral selection factor, the equivalent first voltage regulator time constant, the equivalent second voltage regulator time constant, the equivalent filter time constant and the equivalent regulator gain of the virtual excitation control module of the grid-type energy storage system; An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system; An equivalent grid-type energy storage system model is established based on the equivalent rated capacity, equivalent maximum active output, equivalent maximum reactive output, equivalent actual active output, equivalent actual reactive output, equivalent reactance, equivalent virtual damping, equivalent synchronous reactance and equivalent d-axis transient reactance, equivalent transient time constant, equivalent proportional-integral selection factor, equivalent first voltage regulator time constant, equivalent second voltage regulator time constant, equivalent filter time constant and equivalent virtual excitation control module regulator gain of the grid-type energy storage system; A second model module is configured to calculate static load equivalent parameters of the power distribution area according to the static load equipment data, and establish an equivalent static load model based on the static load equivalent parameters; a third model module, configured to calculate dynamic load equivalent parameters of the power distribution area according to the induction motor data, and establish an equivalent dynamic load model based on the dynamic load equivalent parameters; The fourth model module is configured to calculate reactive compensation parameters for the power distribution area based on the grid-connected energy storage system data, the induction motor data, the static load equipment data, and the substation sending-end data, and establish a reactive compensation model based on the reactive compensation parameters, including: Based on the reactive power balance principle, the reactive compensation parameters of the power distribution area are calculated according to the grid-type energy storage system data, induction motor data, static load equipment data and substation sending end data. The calculation formula of the reactive compensation parameters is: Where, is the reactive power compensation parameter of the power distribution area, Q is the reactive power of the substation sending end, is the sum of the reactive powers actually absorbed by several induction motors in the substation. , and They are the sum of the actual reactive power outputs of all static constant impedance load devices in the substation, the sum of the actual reactive power outputs of all static constant current load devices in the substation, and the actual reactive power outputs of all static constant power load devices. is the sum of the actual reactive power outputs of N grid-type energy storage devices in the grid-type energy storage system; Based on the reactive compensation parameters Establish reactive power compensation model; The load model module is used to construct a comprehensive load model including the grid-type energy storage system based on the equivalent grid-type energy storage system model, the equivalent static load model, the equivalent dynamic load model and the reactive compensation model.

5. A computer-readable storage medium, characterized in that The storage medium stores a computer program, and the computer program is used to execute the method according to any one of claims 1 to 3.

6. An electronic device, characterized in that: The electronic device comprises: processor; a memory for storing instructions executable by the processor; The processor is configured to read the executable instructions from the memory and execute the instructions to implement the method according to any one of claims 1 to 3.

Images