Frequency characteristic evaluation method and device suitable for high-proportion new energy system

Abstract

The invention relates to the technical field of new energy power generation operation control, and particularly provides a frequency characteristic evaluation method and device suitable for a high-proportion new energy system, and the method comprises the steps: obtaining the real-time frequency variation of a generator in the high-proportion new energy system; substituting the real-time frequency variation of the generator in the high-proportion new energy system into a pre-constructed frequency variation function of a bus node and solving to obtain the frequency variation of the bus node; wherein the pre-constructed frequency variation function of the bus node is obtained based on the voltage phasor of the bus node. The technical scheme provided by the invention can quickly and effectively evaluate the frequency characteristic of each node of the high-proportion new energy power generation system, can be used for quickly evaluating the frequency characteristic of a power transmission network, and is also suitable for accurately acquiring the frequency characteristic of a power distribution network; and the method has certain significance for a power grid dispatching department to quickly and accurately position the power grid frequency support weak points and effectively take countermeasures.

Description

Technical Field

[0001] The present invention relates to the field of renewable energy power generation operation control technology, and in particular to a frequency characteristic evaluation method and device suitable for a high-proportion renewable energy system. Background Art

[0002] After a power disturbance occurs in the power system, the frequency dynamics of different nodes in the system are actually different. The root cause is that the frequency dynamics are accompanied by the angle of attack swing process. Before the system reaches synchronization, the relative deviation of the speed of different generator nodes will cause differences in the frequencies of different nodes, causing uneven distribution of frequencies in space. This uneven distribution in space will lead to the asynchronous response of the frequency modulation resources of each node in time, which in turn causes differences in the spatial and temporal distribution of system frequency.

[0003] In traditional power systems, the differentiated characteristics of frequency temporal and spatial distribution objectively exist, but their impact on the system is relatively small. However, in dual-high power systems, the difference in the spatial distribution of grid frequency after power disturbance is more significant. On the one hand, the decrease in system inertia makes the frequency fluctuation after disturbance more intense, and the uneven distribution of unbalanced power during the power angle swing may lead to significant differences in the frequency change rate of each node at the beginning of the power disturbance; on the other hand, the frequency modulation characteristics of new energy are largely determined by its control method when it is connected to the grid through power electronic equipment, and the physical structure of power electronic equipment is different from that of synchronous machines, which makes the heterogeneous characteristics of the frequency modulation characteristics of each node obvious during the frequency dynamic process. Summary of the invention

[0004] In order to overcome the above-mentioned defects, the present invention proposes a frequency characteristic evaluation method and device suitable for a high-proportion new energy system.

[0005] In a first aspect, a frequency characteristic evaluation method suitable for a high-proportion new energy system is provided, the frequency characteristic evaluation method suitable for a high-proportion new energy system comprising:

[0006] Obtain the real-time frequency change of the generator in the high-proportion new energy system;

[0007] Substituting the real-time frequency change of the generator in the high-proportion new energy system into a pre-constructed frequency change function of the bus node and solving it to obtain the frequency change of the bus node;

[0008] The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

[0009] Preferably, the process of obtaining the pre-constructed frequency variation function of the bus node includes:

[0010] Determine the calculation formula of the bus node voltage phasor based on the internal potential phasor of each generator in the power grid and the complex coefficient matrix corresponding to the bus node;

[0011] The pre-constructed frequency variation function of the bus node is determined based on the calculation formula of the bus node voltage phasor.

[0012] Furthermore, the calculation formula of the bus node voltage is as follows:

[0013]

[0014] In the above formula, is the voltage phasor of bus node n, is the element corresponding to the mth column in the complex coefficient matrix corresponding to busbar node n, is the internal potential phasor of generator m, and M is the total number of generators.

[0015] Further, the determining of the pre-constructed frequency variation function of the bus node based on the calculation formula of the bus node voltage phasor includes:

[0016] Obtain the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator, and take the derivative thereof to obtain the first function;

[0017] Deriving the calculation formula of the bus node voltage and inserting it into the first function to obtain a second function;

[0018] The second function is developed using the Euler formula, and the imaginary part is taken to obtain the frequency variation function of the pre-constructed bus node.

[0019] Furthermore, the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator is as follows:

[0020]

[0021] In the above formula, E m is the internal potential amplitude of generator m, e is a natural constant, δ m is the internal potential phase angle of generator m, f N is the rated frequency of the system, △f G is the frequency change of the generator, δ 0 is the initial phase angle of the potential inside the generator, U Bn is the voltage amplitude of bus node n, θ n is the voltage phase angle of bus node n, △f B is the frequency change of the bus node, θ 0 is the initial phase angle of the voltage at the bus node, and j is an imaginary unit.

[0022] Furthermore, the first function is as follows:

[0023]

[0024] Furthermore, the second function is as follows:

[0025]

[0026] In the above formula, △f Gm is the frequency change of generator m, △f Bn is the frequency change of bus node n, D nm is the element amplitude corresponding to the mth column in the complex coefficient matrix corresponding to the busbar node n, α nm is the phase angle of the element corresponding to the mth column in the complex coefficient matrix corresponding to the bus node n.

[0027] Furthermore, the frequency variation function of the pre-built bus node is as follows:

[0028]

[0029] In the above formula, △f Bn is the frequency change of bus node n.

[0030] In a second aspect, a frequency characteristic evaluation device suitable for a high-proportion new energy system is provided, wherein the frequency characteristic evaluation device suitable for a high-proportion new energy system comprises:

[0031] An acquisition module is used to acquire the real-time frequency change of the generator in the high-proportion new energy system;

[0032] An analysis module, used for substituting the real-time frequency variation of the generator in the high-proportion new energy system into a pre-built frequency variation function of the bus node and solving it to obtain the frequency variation of the bus node;

[0033] The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

[0034] In a third aspect, a computer device is provided, comprising: one or more processors;

[0035] The processor is used to store one or more programs;

[0036] When the one or more programs are executed by the one or more processors, the frequency characteristic evaluation method suitable for a high-proportion new energy system is implemented.

[0037] In a fourth aspect, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed, the frequency characteristic evaluation method suitable for a high-proportion new energy system is implemented.

[0038] The above one or more technical solutions of the present invention have at least one or more of the following beneficial effects:

[0039] The present invention provides a method and device for evaluating the frequency characteristics of a high-proportion renewable energy system, including: obtaining the real-time frequency variation of the generator in the high-proportion renewable energy system; substituting the real-time frequency variation of the generator in the high-proportion renewable energy system into a pre-constructed frequency variation function of a bus node and solving it to obtain the frequency variation of the bus node; wherein the pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor. The technical solution provided by the present invention can quickly and effectively evaluate the frequency characteristics of each node in a high-proportion renewable energy power generation system, can be used for the rapid evaluation of the frequency characteristics of the transmission network, and is also suitable for the accurate acquisition of the frequency characteristics of the distribution network, which is of certain significance for the power grid dispatching department to quickly and accurately locate the weak points of the power grid frequency support and effectively take countermeasures. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 It is a schematic flow chart of the main steps of the frequency characteristic evaluation method suitable for a high-proportion new energy system according to an embodiment of the present invention;

[0041] Figure 2 is a structural diagram of the IEEE39 improved model for verifying the effectiveness of the proposed method according to an embodiment of the present invention;

[0042] Figure 3 is a curve chart comparing the time domain simulation and fast calculation of scenario 1 of an embodiment of the present invention;

[0043] Figure 4 is a curve chart comparing the time domain simulation and fast calculation of scenario 2 of an embodiment of the present invention;

[0044] Figure 5 It is a main structural block diagram of a frequency characteristic evaluation device suitable for a high-proportion new energy system according to an embodiment of the present invention. DETAILED DESCRIPTION

[0045] The specific implementation modes of the present invention will be further described in detail below in conjunction with the accompanying drawings.

[0046] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0047] As disclosed in the background technology, after a power disturbance occurs in the power system, the frequency dynamic process of different nodes in the system is actually different. The root cause is that the frequency dynamic process is accompanied by the angle of attack swing process. Before the system reaches synchronization, the relative deviation of the rotation speed of different generator nodes will cause differences in the frequencies of different nodes, causing uneven distribution of frequencies in space. This uneven distribution in space will lead to the asynchronous response of the frequency modulation resources of each node in time, which in turn causes differences in the spatial and temporal distribution of system frequency.

[0048] In traditional power systems, the differentiated characteristics of frequency temporal and spatial distribution objectively exist, but their impact on the system is relatively small. However, in dual-high power systems, the difference in the spatial distribution of grid frequency after power disturbance is more significant. On the one hand, the decrease in system inertia makes the frequency fluctuation after disturbance more intense, and the uneven distribution of unbalanced power during the power angle swing may lead to significant differences in the frequency change rate of each node at the beginning of the power disturbance; on the other hand, the frequency modulation characteristics of new energy are largely determined by its control method when it is connected to the grid through power electronic equipment, and the physical structure of power electronic equipment is different from that of synchronous machines, which makes the heterogeneous characteristics of the frequency modulation characteristics of each node obvious during the frequency dynamic process.

[0049] In order to improve the above problems, the present invention provides a frequency characteristic evaluation method and device suitable for a high-proportion renewable energy system, including: obtaining the real-time frequency change of the generator in the high-proportion renewable energy system; substituting the real-time frequency change of the generator in the high-proportion renewable energy system into a pre-constructed frequency change function of the bus node and solving it to obtain the frequency change of the bus node; wherein the pre-constructed frequency change function of the bus node is obtained based on the bus node voltage phasor. The technical solution provided by the present invention can quickly and effectively evaluate the frequency characteristics of each node in a high-proportion renewable energy power generation system, can be used for the rapid evaluation of the frequency characteristics of the transmission network, and is also suitable for the accurate acquisition of the frequency characteristics of the distribution network, which is of certain significance for the power grid dispatching department to quickly and accurately locate the weak points of the power grid frequency support and effectively take countermeasures.

[0050] The above scheme is described in detail below.

[0051] Example 1

[0052] See attached Figure 1 , Figure 1 FIG. 1 is a flow chart showing the main steps of a frequency characteristic evaluation method for a high-proportion new energy system according to an embodiment of the present invention. Figure 1 As shown, the frequency characteristic evaluation method suitable for a high-proportion new energy system in the embodiment of the present invention mainly includes the following steps:

[0053] Step S101: obtaining the real-time frequency change of the generator in the high-proportion new energy system;

[0054] Step S102: Substituting the real-time frequency variation of the generator in the high-proportion new energy system into a pre-constructed frequency variation function of the bus node and solving it to obtain the frequency variation of the bus node;

[0055] The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

[0056] In this embodiment, the process of obtaining the pre-constructed frequency variation function of the bus node includes:

[0057] Determine the calculation formula of the bus node voltage phasor based on the internal potential phasor of each generator in the power grid and the complex coefficient matrix corresponding to the bus node;

[0058] The pre-constructed frequency variation function of the bus node is determined based on the calculation formula of the bus node voltage phasor.

[0059] In one embodiment, the bus node voltage is calculated as follows:

[0060]

[0061] In the above formula, is the voltage phasor of bus node n, is the element corresponding to the mth column in the complex coefficient matrix corresponding to busbar node n, is the internal potential phasor of generator m, and M is the total number of generators.

[0062] In one embodiment, the determining the pre-constructed frequency variation function of the bus node based on the calculation formula of the bus node voltage phasor includes:

[0063] Obtain the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator, and take the derivative thereof to obtain the first function;

[0064] Deriving the calculation formula of the bus node voltage and inserting it into the first function to obtain a second function;

[0065] The second function is developed using the Euler formula, and the imaginary part is taken to obtain the frequency variation function of the pre-constructed bus node.

[0066] The relationship between the voltage phasor of the bus node and the internal potential phasor of the generator is as follows:

[0067]

[0068] In the above formula, E m is the internal potential amplitude of generator m, e is a natural constant, δ mis the internal potential phase angle of generator m, f N is the rated frequency of the system, △f G is the frequency change of the generator, δ 0 is the initial phase angle of the potential inside the generator, U Bn is the voltage amplitude of bus node n, θ n is the voltage phase angle of bus node n, △f B is the frequency change of the bus node, θ 0 is the initial phase angle of the voltage at the bus node, and j is an imaginary unit.

[0069] Wherein, the first function is as follows:

[0070]

[0071] Wherein, the second function is as follows:

[0072]

[0073] In the above formula, △f Gm is the frequency change of generator m, △f Bn is the frequency change of bus node n, D nm is the element amplitude corresponding to the mth column in the complex coefficient matrix corresponding to the busbar node n, α nm is the phase angle of the element corresponding to the mth column in the complex coefficient matrix corresponding to the bus node n.

[0074] The frequency variation function of the pre-built bus node is as follows:

[0075]

[0076] In the above formula, △f Bn is the frequency change of bus node n.

[0077] Based on the frequency operation data of each grid node (including generator nodes and bus nodes), the lowest frequency point values ​​of different system nodes can be easily obtained, and the weak nodes that affect the frequency stability of the grid can be quickly evaluated. On this basis, auxiliary frequency regulation measures can be used for the weak frequency nodes to enhance the node's ability to support the stable operation of the grid frequency and improve the grid frequency stability performance.

[0078] In order to further illustrate the effectiveness of the proposed method, the effectiveness of the proposed method is verified based on the improved IEEE39-node system. In order to make the frequency characteristics of different nodes in the system more obvious after the power disturbance, the transmission lines 17-18 and 25-26 in the original system are disconnected during the simulation verification. Figure 2 As shown in the figure, two different scenarios are set up to verify the scenario universality of the proposed method:

[0079] Scenario 1: No new energy is connected to the system, and the key frequency regulation parameters of conventional units in the power grid are set to: primary frequency regulation dead zone ±0.033Hz, primary frequency regulation gain 20;

[0080] Scenario 2: Nodes 33 and 34 are connected to a grid-following new energy power generation system, and nodes 35 and 36 are connected to a grid-forming new energy power generation system. The primary frequency regulation dead zone is set to ±0.05Hz, and the primary frequency regulation gain is 20. In this scenario, the grid penetration rate of new energy is 33.72%.

[0081] Obtain the system grid structure and parameters, the internal impedance parameters of the synchronous generator, and the short-circuit impedance parameters of the grid-forming new energy in the two scenarios, and calculate the related admittance matrix D between each bus node and the generator (including synchronous machine and grid-forming new energy) node of the system according to formula (4). Figure 2 The model system shown performs power flow calculations for different operating conditions to obtain the voltage amplitude and phase angle of each node in the network, and the active and reactive power flow values ​​of the synchronous machine.

[0082] Set the simulation event: At t = 2s, the Bus3 node suddenly increases the load disturbance by 200MW, obtain the simulation data of the rotation speed of each generator in the system and the frequency of each bus node, and then obtain the dynamic frequency data of each bus node in the example system according to the pre-built frequency change function of the bus node, and compare it with the actual simulation data of the PowerFactory simulation platform, and select the frequency curves of nodes 1 and 22 in the system that are closer to the generator in electrical distance, and the central node 16 of the grid structure for comparison, such as Figure 3 , Figure 4 shown.

[0083] according to Figure 3 , Figure 4 From the comparison curves, it can be seen that no matter the system is a system with all synchronous power supplies or a system with asynchronous power supplies, the frequency curve obtained by time domain simulation and the frequency curve estimated according to the frequency space correlation coefficient matrix G have a high degree of overlap, indicating that the frequencies of each node in the system have spatial similarity, which also verifies the accuracy of the proposed method.

[0084] Example 2

[0085] Based on the same inventive concept, the present invention also provides a frequency characteristic evaluation device suitable for a high-proportion new energy system, such as Figure 5 As shown,

[0086] The frequency characteristic evaluation device suitable for a high-proportion new energy system comprises:

[0087] An acquisition module is used to acquire the real-time frequency change of the generator in the high-proportion new energy system;

[0088] An analysis module, used for substituting the real-time frequency variation of the generator in the high-proportion new energy system into a pre-built frequency variation function of the bus node and solving it to obtain the frequency variation of the bus node;

[0089] The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

[0090] Preferably, the process of obtaining the pre-constructed frequency variation function of the bus node includes:

[0091] Determine the calculation formula of the bus node voltage phasor based on the internal potential phasor of each generator in the power grid and the complex coefficient matrix corresponding to the bus node;

[0092] The pre-constructed frequency variation function of the bus node is determined based on the calculation formula of the bus node voltage phasor.

[0093] Furthermore, the calculation formula of the bus node voltage is as follows:

[0094]

[0095] In the above formula, is the voltage phasor of bus node n, is the element corresponding to the mth column in the complex coefficient matrix corresponding to busbar node n, is the internal potential phasor of generator m, and M is the total number of generators.

[0096] Further, the determining of the pre-constructed frequency variation function of the bus node based on the calculation formula of the bus node voltage phasor includes:

[0097] Obtain the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator, and take the derivative thereof to obtain the first function;

[0098] Deriving the calculation formula of the bus node voltage and inserting it into the first function to obtain a second function;

[0099] The second function is developed using the Euler formula, and the imaginary part is taken to obtain the frequency variation function of the pre-constructed bus node.

[0100] Furthermore, the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator is as follows:

[0101]

[0102] In the above formula, E mis the internal potential amplitude of generator m, e is a natural constant, δ m is the internal potential phase angle of generator m, f N is the rated frequency of the system, △f G is the frequency change of the generator, δ 0 is the initial phase angle of the potential inside the generator, U Bn is the voltage amplitude of bus node n, θ n is the voltage phase angle of bus node n, △f B is the frequency change of the bus node, θ 0 is the initial phase angle of the voltage at the bus node, and j is an imaginary unit.

[0103] Furthermore, the first function is as follows:

[0104]

[0105] Furthermore, the second function is as follows:

[0106]

[0107] In the above formula, △f Gm is the frequency change of generator m, △f Bn is the frequency change of bus node n, D nm is the element amplitude corresponding to the mth column in the complex coefficient matrix corresponding to the busbar node n, α nm is the phase angle of the element corresponding to the mth column in the complex coefficient matrix corresponding to the bus node n.

[0108] Furthermore, the frequency variation function of the pre-built bus node is as follows:

[0109]

[0110] In the above formula, △f Bn is the frequency change of bus node n.

[0111] Example 3

[0112] Based on the same inventive concept, the present invention also provides a computer device, which includes a processor and a memory, the memory is used to store a computer program, the computer program includes program instructions, and the processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, which is suitable for implementing one or more instructions, and is specifically suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function, so as to implement the steps of a frequency characteristic evaluation method suitable for a high-proportion new energy system in the above embodiment.

[0113] Example 4

[0114] Based on the same inventive concept, the present invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device for storing programs and data. It is understandable that the computer-readable storage medium here can include both built-in storage media in a computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space that stores the operating system of the terminal. In addition, one or more instructions suitable for being loaded and executed by a processor are also stored in the storage space, and these instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the steps of a frequency characteristic evaluation method suitable for a high-proportion new energy system in the above embodiment.

[0115] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0116] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0117] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0118] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0119] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A frequency characteristic evaluation method suitable for a high-proportion new energy system, characterized in that: The method comprises: Obtain the real-time frequency change of the generator in the high-proportion new energy system; Substituting the real-time frequency change of the generator in the high-proportion new energy system into a pre-constructed frequency change function of the bus node and solving it to obtain the frequency change of the bus node; The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

2. The method according to claim 1, characterized in that The process of obtaining the frequency variation function of the pre-built bus node includes: Determine the calculation formula of the bus node voltage phasor based on the internal potential phasor of each generator in the power grid and the complex coefficient matrix corresponding to the bus node; The pre-constructed frequency variation function of the bus node is determined based on the calculation formula of the bus node voltage phasor.

3. The method according to claim 2, characterized in that The calculation formula of the bus node voltage is as follows: In the above formula, is the voltage phasor of bus node n, is the element corresponding to the mth column in the complex coefficient matrix corresponding to busbar node n, is the internal potential phasor of generator m, and M is the total number of generators.

4. The method according to claim 3, characterized in that The determining the pre-constructed frequency variation function of the bus node based on the calculation formula of the bus node voltage phasor comprises: Obtain the relationship between the voltage phasor of the bus node and the internal potential phasor of the generator, and take the derivative thereof to obtain the first function; Deriving the calculation formula of the bus node voltage and inserting it into the first function to obtain a second function; The second function is developed using the Euler formula, and the imaginary part is taken to obtain the frequency variation function of the pre-constructed bus node.

5. The method according to claim 4, characterized in that The relationship between the voltage phasor of the bus node and the internal potential phasor of the generator is as follows: In the above formula, E m is the internal potential amplitude of generator m, e is a natural constant, δ m is the internal potential phase angle of generator m, f N is the rated frequency of the system, △f G is the frequency change of the generator, δ0 is the initial phase angle of the potential inside the generator, U Bn is the voltage amplitude of bus node n, θ n is the voltage phase angle of bus node n, △f B is the frequency change of the bus node, θ0 is the initial phase angle of the voltage at the bus node, and j is an imaginary unit.

6. The method according to claim 5, characterized in that The first function is as follows:

7. The method according to claim 6, characterized in that The second function is as follows: In the above formula, △f Gm is the frequency change of generator m, △f Bn is the frequency change of bus node n, D nm is the element amplitude corresponding to the mth column in the complex coefficient matrix corresponding to the busbar node n, α nm is the phase angle of the element corresponding to the mth column in the complex coefficient matrix corresponding to the bus node n.

8. The method according to claim 7, characterized in that The frequency variation function of the pre-built bus node is as follows:

9. A frequency characteristic evaluation device suitable for a high-proportion new energy system, characterized in that: The device comprises: An acquisition module is used to acquire the real-time frequency change of the generator in the high-proportion new energy system; An analysis module, used for substituting the real-time frequency variation of the generator in the high-proportion new energy system into a pre-built frequency variation function of the bus node and solving it to obtain the frequency variation of the bus node; The pre-constructed frequency variation function of the bus node is obtained based on the bus node voltage phasor.

10. A computer device, characterized in that: include: one or more processors; The processor is used to store one or more programs; When the one or more programs are executed by the one or more processors, the frequency characteristic evaluation method suitable for a high-proportion new energy system as claimed in any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and when the computer program is executed, the frequency characteristic evaluation method suitable for a high-proportion new energy system as claimed in any one of claims 1 to 8 is implemented.

Images