Energy storage frequency modulation instruction prediction method and system

Abstract

The invention relates to the technical field of energy storage frequency modulation, in particular to an energy storage frequency modulation instruction prediction method and system, and the method comprises the steps: decomposing a frequency modulation instruction sequence into a plurality of subsequences; determining the sub-sequence with the maximum frequency change in all the sub-sequences as a target sub-sequence, and generating a target virtual sequence to replace the target sub-sequence; predicting the target virtual sequence and other subsequences by adopting a GRU network to obtain a plurality of corresponding prediction results; and superposing all the prediction results to obtain a target frequency modulation instruction prediction result. The method has the beneficial effects that the frequency modulation instruction sequence is decomposed into a plurality of subsequences by using the variational mode decomposition algorithm, and then the subsequences are processed and predicted through the frequency method and the virtual sequence method, so that the frequency modulation instruction is predicted in advance, an energy storage system can respond more quickly, the response time difference is reduced, and the frequency modulation performance index is improved; and higher economic benefits are created for a power plant.

Description

Technical Field

[0001] The present invention relates to the technical field of energy storage frequency modulation, and in particular to a method and system for predicting energy storage frequency modulation commands. Background Art

[0002] With the continuous increase in the penetration rate of renewable energy in the power grid, the frequency modulation demand of the power grid has also increased accordingly. Energy storage systems, especially hybrid energy storage systems (such as the combination of supercapacitors and lithium-ion batteries), play an important role in assisting thermal power units in frequency modulation. However, there are some technical limitations in existing hybrid energy storage assisted frequency modulation methods.

[0003] In traditional frequency modulation methods, the difference between the frequency modulation command and the output of the thermal power unit is transmitted to the hybrid energy storage system, where the low-frequency part is borne by the battery and the high-frequency part is borne by the supercapacitor. A major problem with this method is that both signal transmission and the response of the energy storage system require time, resulting in a response time delay, which in turn affects the benefit value of frequency modulation and the overall economic efficiency of the power plant.

[0004] In addition, existing prediction methods, such as using variational mode decomposition (VMD) to decompose the original sequence into multiple subsequences, then predicting each subsequence separately, and finally superimposing the prediction results, have some defects. In particular, some of the subsequences after VMD decomposition have extremely high frequencies, and these high-frequency subsequences have strong randomness and nonlinearity, resulting in difficult prediction and low prediction accuracy. Summary of the Invention

[0005] To solve the above technical problems, the present invention provides the following technical solutions: In a first aspect, the present invention provides a method for predicting energy storage frequency modulation commands, including decomposing the frequency modulation command sequence into multiple subsequences; Determining the subsequence with the largest frequency change among all subsequences as the target subsequence, and generating a target virtual sequence to replace the target subsequence; Using a GRU network to predict the target virtual sequence and the remaining subsequences to obtain corresponding multiple prediction results; Superimposing all the prediction results to obtain the target frequency modulation command prediction result; Among them, determining the subsequence with the largest frequency change among all subsequences as the target subsequence includes: Calculating the frequency change function corresponding to each subsequence; Calculating the corresponding integral value according to each frequency change function; Determining the target subsequence according to the magnitudes of all integral values.

[0006] As a preferred embodiment of the method for predicting energy storage frequency modulation commands of the present invention, where: decomposing the frequency modulation command sequence into multiple subsequences includes: Decompose the frequency modulation command sequence using the variational mode decomposition algorithm to obtain multiple intrinsic mode function subsequences.

[0007] As a preferred embodiment of the energy storage frequency modulation command prediction method of the present invention, wherein: determining the target subsequence according to the magnitudes of all integral values includes: Calculate the absolute value corresponding to each integral value; Determine the subsequence corresponding to the maximum value among all absolute values as the target subsequence.

[0008] As a preferred embodiment of the energy storage frequency modulation command prediction method of the present invention, wherein: calculating the frequency change function corresponding to each subsequence includes: Perform a difference operation on each intrinsic mode function subsequence to generate a difference sequence; Perform an exponential calculation on each value of the difference sequence to generate a difference exponential sequence.

[0009] As a preferred embodiment of the energy storage frequency modulation command prediction method of the present invention, wherein: generating a target virtual sequence to replace the target subsequence includes: In the coordinate system, perform a symmetry transformation on the target subsequence based on the horizontal line of Y = value M to generate a target virtual sequence; Value M is the value of M corresponding to the minimum sum of the target function; The target function is: , wherein, X i represents the i-th data point in the target subsequence, N is the total number of data points in the target subsequence, and M is an optimal value used to generate the target virtual sequence.

[0010] As a preferred embodiment of the energy storage frequency modulation command prediction method of the present invention, wherein: using a GRU network to predict the target virtual sequence and the remaining subsequences to obtain corresponding multiple prediction results further includes: correcting the prediction result of the target virtual sequence.

[0011] In a second aspect, the present invention provides an energy storage frequency modulation command prediction system, including a decomposition module: used for decomposing the frequency modulation command sequence into multiple subsequences; A generation module: used for determining the subsequence with the largest frequency change among all subsequences as the target subsequence, and generating a target virtual sequence to replace the target subsequence; A prediction module: used for using a GRU network to predict the target virtual sequence and the remaining subsequences to obtain corresponding multiple prediction results; An overlay module: used for overlaying all prediction results to obtain the target frequency modulation command prediction result.

[0012] In a third aspect, the present invention provides a computer device, including a memory and a processor. The memory stores a computer program, and is characterized in that the processor implements the steps of an energy storage frequency modulation instruction prediction method when executing the computer program.

[0013] In a fourth aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, and is characterized in that the computer program implements the steps of an energy storage frequency modulation instruction prediction method when being executed by a processor.

[0014] Advantages of the present invention: Effectively solve the problems of delayed frequency modulation response time and low prediction accuracy in the prior art for hybrid energy storage assisting thermal power units in frequency modulation. By using the variational mode decomposition (VMD) algorithm to decompose the frequency modulation instruction sequence into multiple subsequences, and then processing and predicting the subsequences through the frequency method and the virtual sequence method, the frequency modulation instructions can be predicted in advance, enabling the energy storage system to respond more quickly, reducing the response time difference, improving the frequency modulation performance index, and creating higher economic benefits for the power plant. In addition, by optimizing the algorithm to find the optimal value to generate a symmetric virtual sequence, the volatility of the high-frequency sequence is further reduced, significantly improving the prediction accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following will briefly introduce the drawings required for description in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0016] Figure 1 It is a schematic flowchart of the energy storage frequency modulation instruction prediction method.

[0017] Figure 2 It is a schematic flowchart of the traditional prediction method.

[0018] Figure 3 It is a schematic diagram before the change of the target subsequence.

[0019] Figure 4 It is a schematic diagram after the change of the target subsequence. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0020] To make the above objects, features, and advantages of the present invention more obvious and understandable, the following will give a detailed description of the specific embodiments of the present invention in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the scope of protection of the present invention.

[0021] Example 1, referring to Figures 1 to 4 , which is the first embodiment of the present invention. This embodiment provides a method for predicting energy storage frequency modulation commands, including: S1. Decompose the frequency modulation command sequence into multiple subsequences.

[0022] Furthermore, decomposing the frequency modulation command sequence into multiple subsequences includes: Use the variational mode decomposition algorithm to decompose the frequency modulation command sequence to obtain multiple intrinsic mode function subsequences.

[0023] It should be noted that assuming the frequency modulation command is Pt, after VMD decomposition, it becomes IMF1, IMF2, IMF3,..., IMF k (where: IMF k is the k-th subsequence), and further assuming that a subsequence IMF i = [X i1 , X i2 , X i3 ,..., X iN (where: X iN is the value of the i-th IMF subsequence at the N-th time point).

[0024] S2. Determine the subsequence with the largest frequency change among all subsequences as the target subsequence, and generate a target virtual sequence to replace the target subsequence.

[0025] Furthermore, determining the subsequence with the largest frequency change among all subsequences as the target subsequence includes: Calculate the frequency change function corresponding to each subsequence.

[0026] Furthermore, calculating the frequency change function corresponding to each subsequence includes: Perform differential processing on each intrinsic mode function subsequence to generate a difference sequence, where the difference sequence is obtained by calculating the differences between pairwise elements; Perform exponential calculation on each value of the difference sequence to generate a difference index sequence, where generating the difference numerical sequence is to quantify the relative importance of the differences and enhance the influence of larger differences through an exponential function; Perform fitting or simulation on the difference index sequence to generate a frequency change function.

[0027] It should be noted that the difference sequence obtained after differential processing of the i-th subsequence is: IMF ci = [X i1 - X i2 , X i1 - X i3 , X i1 - Xi4 ,...,X i1 -X iN ,X i2 -X i3， X i2 -X i4 ,...,X i(N-1) -X iN (There are N(N - 1) / 2 values) (where: X iN is the value of the current subsequence at the Nth time point, and X i(N-1) is the value of the current subsequence at the (N - 1)th time point) Furthermore, multiply each value in the above difference sequence by the exponential function e () The resulting difference exponential sequence is: =[(X i1 -X i2 )(X i1 -X i3 ),(X i1 -X i4 ),...,(X i1 -X iN ),(X i2 -X i3 ),(X i2 -X i4 ),...,(X i(N-1) -X iN )] Furthermore, construct a function f(t) through fitting or simulation software i , and this function takes different values at different time points t, and these values are determined by the difference exponential sequence. Specifically: When t = 1, f(t) i =(X i1 -X i2 ); When t = 2, f(t) i =(X i1 -X i3 ); When t = 3, f(t) i =(X i1 -X i4 ); …… When t = N(N - 1) / 2, f(t) i =(X i(N-1) -X iN )。

[0028] Furthermore, the frequency change function is: , Among them, F(t) is the frequency variation function; A is the augmentation coefficient, specifically the natural logarithm of the ratio of the maximum value to the average value of the value difference index sequence; is the exponential decay function, used to calculate the integral value of the frequency variation function F(t).

[0029] Calculate the corresponding integral value according to each frequency variation function.

[0030] It should be noted that the integral value represents the area enclosed by f(t) and the X-axis, and is used to quantify the frequency variation of the subsequence. The specific calculation formula is: , where T is the total length of the time series; A i represents the f(t) i area enclosed by the i-th subsequence and the X-axis; f(t) i represents the frequency variation function of the i-th subsequence.

[0031] Determine the target subsequence according to the magnitudes of all integral values.

[0032] Furthermore, determining the target subsequence according to the magnitudes of all integral values includes: Calculate the absolute value corresponding to each integral value; Determine the subsequence corresponding to the maximum value among all absolute values as the target subsequence.

[0033] Furthermore, generating a target virtual sequence to replace the target subsequence includes: In the coordinate system, perform a symmetry transformation on the target subsequence based on the horizontal line of Y = value M to generate a target virtual sequence; The value M is the value of M corresponding to the minimum of the target function; The target function is: , where X i represents the i-th data point in the target subsequence, N is the total number of data points in the target subsequence, M is an optimal value used to generate the target virtual sequence, sin(X i ) is the sine function, and its input is the data point X i .

[0034] It should be noted that the solution method of the value M can use iterative loops or optimization algorithms such as (WOA, PSO), etc. Assume the target subsequence is IMF i = [X i1 , X i2, X i3 , ……, X iN (where X iN is the value of the current subsequence at the Nth time point), establish a coordinate system, the abscissa of the coordinate system represents the subscript of the data points in the target subsequence, and the ordinate represents the value of the specific data points in the target subsequence, attached Figure 4 is the illustration after symmetric change through the horizontal line of Y = value M, and the target virtual sequence after the change is used to replace the target subsequence for prediction. Figure 3 After passing through the horizontal line of Y = value M for symmetric change, the illustration is shown. The target virtual sequence after the change is used to replace the target subsequence for prediction.

[0035] S3. Use the GRU network to predict the target virtual sequence and the remaining subsequences to obtain corresponding multiple prediction results.

[0036] Furthermore, using the GRU network to predict the target virtual sequence and the remaining subsequences to obtain corresponding multiple prediction results also includes: correcting the prediction results of the target virtual sequence.

[0037] It should be noted that let the prediction result be [X i(N+1) , X i(N+2) , ……, X i(N+0.1N) , where N is the length of the sequence, i is the prediction step, N + 0.1N represents adding 10% of the sequence length to the original length N of the sequence as the prediction time point, and the correction formula is: , In the formula, sigmoid(N + i) is an activation function used to control the correction amplitude; Q is a constant used to adjust the output of the activation function to ensure that the corrected predicted value is within a reasonable range; X i(N+1) : represents the predicted value of the ith subsequence at the N + 1 time point; X i(N+1) (the second one on the right side of the equation): represents the predicted value after correction.

[0038] S4. Superimpose all the prediction results to obtain the target frequency modulation command prediction result.

[0039] It should be noted that superimpose all the prediction results predicted separately (using the above prediction method) (this is the prior art and will not be elaborated here) to obtain the final target frequency modulation command prediction result.

[0040] In summary, a method for predicting energy storage frequency modulation commands of the present invention effectively solves the problems of delayed frequency modulation response time and low prediction accuracy of hybrid energy storage assisting thermal power units in the prior art. By using the variational mode decomposition (VMD) algorithm to decompose the frequency modulation command sequence into multiple subsequences, and then processing and predicting the subsequences through the frequency method and the virtual sequence method, the frequency modulation commands can be predicted in advance, enabling the energy storage system to respond more quickly, reducing the response time difference, improving the frequency modulation performance index, and creating higher economic benefits for the power plant. In addition, by using an optimization algorithm to find the optimal value to generate a symmetric virtual sequence, the volatility of the high-frequency sequence is further reduced, significantly improving the prediction accuracy.

[0041] Example 2 is the second example of the present invention. To verify the advantages of the present invention, the method of the present invention and the prediction method of VMD-GRU (a hybrid model combining variational mode decomposition and gated recurrent unit) are used to predict the frequency modulation sequence, and the results are shown in Table 1: Table 1 Comparison table of prediction results , It should be noted that the calculation formulas and definitions of the above four evaluation indexes are shown in Table 2: Table 2 Operation formulas and definitions of each index , Among them, N represents the sample size, and and respectively represent the actual value and the predicted value at time n.

[0042] In summary, according to the experimental data, the errors of the prediction using this method are all smaller than the results of the VMD-GRU prediction, fully demonstrating that the prediction using the present invention has better effects.

[0043] Example 3 is the third example of the present invention, which is different from the previous two examples in that: If the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0044] The logic and / or steps represented in the flowchart or otherwise described herein can be considered, for example, as a definitional sequence of executable instructions for implementing a logical function, which can be specifically implemented in any computer-readable medium for use by an instruction execution system, apparatus, or device (such as a computer-based system, a system including a processor, or other systems that can fetch and execute instructions from the instruction execution system, apparatus, or device), or used in conjunction with these instruction execution systems, apparatus, or devices. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0045] More specific examples (non-exhaustive list) of computer-readable media include the following: an electrical connection portion with one or more wirings (electronic device), a portable computer disk cartridge (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, a computer-readable medium can even be paper or other suitable media on which a program can be printed, because the program can be obtained electronically, for example, by optically scanning the paper or other media, followed by editing, interpretation, or other appropriate processing as necessary, and then stored in a computer memory.

[0046] It should be understood that various parts of the present invention can be implemented by hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, any one or a combination of the following techniques well known in the art can be used: discrete logic circuits with logic gate circuits for implementing logical functions on data signals, application-specific integrated circuits with appropriate combinational logic gate circuits, programmable gate arrays (PGAs), field programmable gate arrays (FPGAs), etc.

[0047] Embodiment 4 is the fourth embodiment of the present invention. This embodiment provides an energy storage frequency modulation instruction prediction system, including a decomposition module: used to decompose a frequency modulation instruction sequence into multiple subsequences; A generation module: used to determine the subsequence with the largest frequency change among all subsequences as the target subsequence, and generate a target virtual sequence to replace the target subsequence; A prediction module: used to predict the target virtual sequence and the remaining subsequences using a GRU network to obtain corresponding multiple prediction results; An overlay module: used to overlay all the prediction results to obtain the target frequency modulation instruction prediction result.

[0048] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit them. Although the present invention has been described in detail with reference to the preferred embodiments, those of ordinary skill in the art should understand that the technical solutions of the present invention can be modified or equivalently replaced without departing from the spirit and scope of the technical solutions of the present invention, and they should all be covered within the scope of the claims of the present invention.

Claims

1. A method for predicting energy storage frequency modulation instructions, characterized in that: include, Decomposing the frequency modulation instruction sequence into multiple subsequences; Determine the subsequence with the largest frequency change among all subsequences as the target subsequence, and generate a target virtual sequence to replace the target subsequence; The GRU network is used to predict the target virtual sequence and other subsequences to obtain corresponding multiple prediction results; Superimposing all the prediction results to obtain the target frequency modulation instruction prediction result; The step of determining the subsequence with the largest frequency change among all subsequences as the target subsequence includes: Calculate the frequency change function corresponding to each subsequence; Calculate the corresponding integral value according to each frequency variation function; The target subsequence is determined according to the size of all integral values.

2. A method for predicting energy storage frequency modulation instructions according to claim 1, characterized in that: Decomposing the frequency modulation instruction sequence into a plurality of subsequences comprises: The frequency modulation instruction sequence is decomposed using the variational mode decomposition algorithm to obtain multiple intrinsic mode function subsequences.

3. A method for predicting energy storage frequency modulation instructions as claimed in claim 2, characterized in that: Determining the target subsequence according to the magnitude of all integral values ​​includes: Calculate the absolute value corresponding to each integral value; The subsequence corresponding to the maximum value among all absolute values ​​is determined as the target subsequence.

4. A method for predicting energy storage frequency modulation instructions as claimed in claim 3, characterized in that: The calculating of the frequency change function corresponding to each subsequence comprises: Performing difference processing on each intrinsic mode function subsequence to generate a difference sequence; Perform exponential calculation on each value of the difference sequence to generate a difference exponential sequence; Fit or simulate the difference exponential series to generate a frequency change function.

5. A method for predicting energy storage frequency modulation instructions as claimed in claim 4, characterized in that: Generating a target virtual sequence to replace the target subsequence comprises: In the coordinate system, the target subsequence is symmetrically transformed based on the horizontal line of Y=value M to generate a target virtual sequence; The value M is the value of M corresponding to the minimum objective function sum; The objective function is: ， Among them, X i represents the i-th data point in the target subsequence, N is the total number of data points in the target subsequence, and M is an optimal value used to generate the target virtual sequence.

6. A method for predicting energy storage frequency modulation instructions as claimed in claim 5, characterized in that: The method of using the GRU network to predict the target virtual sequence and the remaining subsequences to obtain the corresponding multiple prediction results also includes: correcting the prediction results of the target virtual sequence.

7. A system using the energy storage frequency modulation instruction prediction method according to any one of claims 1 to 6, characterized in that: include, Decomposition module: used to decompose the frequency modulation instruction sequence into multiple sub-sequences; A generation module is used to determine the subsequence with the largest frequency change among all subsequences as a target subsequence, and generate a target virtual sequence to replace the target subsequence; Prediction module: used to use the GRU network to predict the target virtual sequence and other subsequences to obtain corresponding multiple prediction results; Superposition module: used to superimpose all the prediction results to obtain the target frequency modulation instruction prediction result.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of a method for predicting energy storage frequency modulation instructions as described in any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of a method for predicting energy storage frequency modulation instructions as described in any one of claims 1 to 6 are implemented.

Images