CN117748501B - A wind power prediction method and system for energy storage assisted black start - Google Patents

A wind power prediction method and system for energy storage assisted black start Download PDF

Info

Publication number
CN117748501B
CN117748501B CN202410182403.6A CN202410182403A CN117748501B CN 117748501 B CN117748501 B CN 117748501B CN 202410182403 A CN202410182403 A CN 202410182403A CN 117748501 B CN117748501 B CN 117748501B
Authority
CN
China
Prior art keywords
wind speed
modal
component
components
obtaining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410182403.6A
Other languages
Chinese (zh)
Other versions
CN117748501A (en
Inventor
李志鹏
兀鹏越
赵俊博
寇水潮
郝博瑜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Thermal Power Research Institute Co Ltd
Original Assignee
Xian Thermal Power Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Thermal Power Research Institute Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202410182403.6A priority Critical patent/CN117748501B/en
Publication of CN117748501A publication Critical patent/CN117748501A/en
Application granted granted Critical
Publication of CN117748501B publication Critical patent/CN117748501B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请属于风电技术领域,具体提出一种储能辅助黑启动的风功率预测方法和系统,该方法包括获取风速序列;利用CEEMD算法对风速序列进行分解获得多个模态分量;基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;基于各误差系数的正负性对所有模态分量进行重构以得到多个互补模态分量;基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值进而利用风速预测值得到风功率预测值以在停电故障时利用风功率预测值参与储能辅助黑启动。利用本申请的方法能够提高风功率的预测精度。

The present application belongs to the field of wind power technology, and specifically proposes a wind power prediction method and system for energy storage-assisted black start, the method comprising obtaining a wind speed sequence; decomposing the wind speed sequence using the CEEMD algorithm to obtain multiple modal components; obtaining the corresponding wind speed component prediction value based on each modal component using the first recurrent neural network model, and obtaining the error coefficient corresponding to each modal component based on each wind speed component prediction value and each modal component; reconstructing all modal components based on the positive and negative properties of each error coefficient to obtain multiple complementary modal components; obtaining the corresponding target wind speed component prediction value based on each complementary modal component using the second recurrent neural network model; obtaining the wind speed prediction value of the next sampling point of the current sampling point based on each target wind speed component prediction value, and then using the wind speed prediction value to obtain the wind power prediction value so as to use the wind power prediction value to participate in energy storage-assisted black start in the event of a power outage. The method of the present application can improve the prediction accuracy of wind power.

Description

一种储能辅助黑启动的风功率预测方法和系统A wind power prediction method and system for energy storage assisted black start

技术领域Technical Field

本申请涉及风电技术领域,尤其涉及一种储能辅助黑启动的风功率预测方法和系统。The present application relates to the field of wind power technology, and in particular to a method and system for predicting wind power during energy storage-assisted black start.

背景技术Background technique

随着电网规模增大,大停电事故后果越来越严重;受自然资源约束,风光充裕地区的传统黑启动电源不足。将风光储系统作为黑启动电源,可以提高区域电网的黑启动能力。风光储发电系统作为黑启动电源时,考虑储能装置的充/放电功率约束和电量约束,在黑启动过程中,风电场和光伏电站出力不足或波动剧烈时,可能会出现储能过充过放的情况,导致储能无法继续利用,使黑启动失败。因此为了更好地在停电事故后利用风光储发电系统进行黑启动,需要根据历史风速数据,对风电场的持续有效出力概率风功率进行评估。然而传统的风功率预测存在预测精度差的问题。As the scale of the power grid increases, the consequences of large-scale power outages are becoming more and more serious; due to natural resource constraints, traditional black start power sources in areas with abundant wind and solar power are insufficient. Using wind, solar and storage systems as black start power sources can improve the black start capability of regional power grids. When the wind, solar and storage power generation system is used as a black start power source, the charge/discharge power constraints and power constraints of the energy storage device are considered. During the black start process, when the output of wind farms and photovoltaic power stations is insufficient or fluctuates violently, the energy storage may be overcharged or over-discharged, resulting in the inability to continue to use the energy storage, causing the black start to fail. Therefore, in order to better use the wind, solar and storage power generation system for black start after a power outage, it is necessary to evaluate the wind power of the wind farm's continuous effective output probability based on historical wind speed data. However, traditional wind power forecasting has the problem of poor prediction accuracy.

发明内容Summary of the invention

本申请旨在至少在一定程度上解决相关技术中的技术问题之一。The present application aims to solve one of the technical problems in the related art at least to some extent.

为此,本申请的第一个目的在于提出一种储能辅助黑启动的风功率预测方法,以提高风功率的预测精度。To this end, the first objective of the present application is to propose a wind power prediction method for energy storage-assisted black start to improve the prediction accuracy of wind power.

本申请的第二个目的在于提出一种储能辅助黑启动的风功率预测系统。The second objective of the present application is to propose a wind power prediction system for energy storage-assisted black start.

本申请的第三个目的在于提出一种电子设备。The third objective of the present application is to provide an electronic device.

本申请的第四个目的在于提出一种计算机可读存储介质。A fourth objective of the present application is to provide a computer-readable storage medium.

为达上述目的,本申请第一方面实施例提出了一种储能辅助黑启动的风功率预测方法,包括:To achieve the above-mentioned purpose, the first embodiment of the present application proposes a wind power prediction method for energy storage-assisted black start, comprising:

获取风速序列,所述风速序列包括当前采样点和多个历史采样点的风速测量值;Acquire a wind speed sequence, wherein the wind speed sequence includes wind speed measurement values of a current sampling point and multiple historical sampling points;

利用CEEMD算法对所述风速序列进行分解获得多个模态分量;Decomposing the wind speed sequence by using CEEMD algorithm to obtain multiple modal components;

基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;Based on each modal component, a corresponding wind speed component prediction value is obtained using a first recurrent neural network model, and based on each wind speed component prediction value and each modal component, an error coefficient corresponding to each modal component is obtained;

基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量;Based on the positive and negative properties of each error coefficient, all modal components are reconstructed to obtain multiple complementary modal components;

基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;Based on each complementary modal component, a second recurrent neural network model is used to obtain a corresponding target wind speed component prediction value;

基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用所述风功率预测值参与储能辅助黑启动。The wind speed prediction value of the next sampling point of the current sampling point is obtained based on the prediction value of each target wind speed component, and the wind power prediction value is then used to obtain the wind speed prediction value, so as to participate in energy storage-assisted black start when a power outage occurs.

在本申请的第一方面的方法中,每个模态分量由所述风速序列的所有采样点的风速测量值分解得到的风速分量组成,风速分量数量等于所述风速序列的采样点数量;设置第一循环神经网络模型的参数,以使第一循环神经网络模型输出预设数量的风速分量预测值,模型输出的预设数量的风速分量预测值对应的实际值为输入的模态分量中当前采样点及之前的相同数量的风速分量;所述基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数,包括:针对任一模态分量,将该模态分量输入第一循环神经网络模型获得对应的预设数量的风速分量预测值;基于预设数量的风速分量预测值和对应的实际值获得该模态分量的误差系数,进而得到各模态分量的误差系数。In the method of the first aspect of the present application, each modal component is composed of wind speed components obtained by decomposing the wind speed measurement values of all sampling points of the wind speed sequence, and the number of wind speed components is equal to the number of sampling points of the wind speed sequence; the parameters of the first recurrent neural network model are set so that the first recurrent neural network model outputs a preset number of wind speed component prediction values, and the actual values corresponding to the preset number of wind speed component prediction values output by the model are the same number of wind speed components at the current sampling point and before in the input modal component; the method of obtaining the corresponding wind speed component prediction value based on each modal component using the first recurrent neural network model, and obtaining the error coefficient corresponding to each modal component based on each wind speed component prediction value and each modal component, includes: for any modal component, inputting the modal component into the first recurrent neural network model to obtain the corresponding preset number of wind speed component prediction values; obtaining the error coefficient of the modal component based on the preset number of wind speed component prediction values and the corresponding actual values, and then obtaining the error coefficient of each modal component.

在本申请的第一方面的方法中,所述基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量,包括:基于各模态分量的误差系数的正负性,将所有模态分量划分为非负数组和负数组;按误差系数,将非负数组中模态分量从大到小进行排列得到目标非负数组,将负数组中模态分量从小到大进行排列得到目标负数组;将目标非负数组与目标负数组中对应位置的模态分量相加从而得到多个互补模态分量。In the method of the first aspect of the present application, all modal components are reconstructed based on the positive or negative nature of each error coefficient to obtain a plurality of complementary modal components, including: based on the positive or negative nature of the error coefficient of each modal component, all modal components are divided into a non-negative array and a negative array; according to the error coefficient, the modal components in the non-negative array are arranged from large to small to obtain a target non-negative array, and the modal components in the negative array are arranged from small to large to obtain a target negative array; and the modal components at corresponding positions in the target non-negative array and the target negative array are added to obtain a plurality of complementary modal components.

在本申请的第一方面的方法中,设置第二循环神经网络模型的参数,以使第二循环神经网络模型输出当前采样点的下一采样点的目标风速分量预测值;所述基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值,包括:针对任一互补模态分量,将该互补模态分量输入第二循环神经网络模型获得对应的当前采样点的下一采样点的目标风速分量预测值,进而得到各互补模态分量对应的目标风速分量预测值。In the method of the first aspect of the present application, the parameters of the second recurrent neural network model are set so that the second recurrent neural network model outputs a target wind speed component prediction value for the next sampling point of the current sampling point; the method of obtaining the corresponding target wind speed component prediction value based on each complementary modal component using the second recurrent neural network model includes: for any complementary modal component, inputting the complementary modal component into the second recurrent neural network model to obtain the corresponding target wind speed component prediction value for the next sampling point of the current sampling point, and then obtaining the target wind speed component prediction value corresponding to each complementary modal component.

在本申请的第一方面的方法中,所述基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,包括:将各目标风速分量预测值进行求和从而得到当前采样点的下一采样点的风速预测值。In the method of the first aspect of the present application, obtaining the wind speed prediction value of the next sampling point of the current sampling point based on each target wind speed component prediction value includes: summing the target wind speed component prediction values to obtain the wind speed prediction value of the next sampling point of the current sampling point.

在本申请的第一方面的方法中,所述第一循环神经网络模型和所述第二循环神经网络模型分别采用GRU模型。In the method of the first aspect of the present application, the first recurrent neural network model and the second recurrent neural network model respectively adopt GRU models.

在本申请的第一方面的方法中,各模态分量的误差系数满足:,其中/>为第k个模态分量的误差系数,N为第k个模态分量对应的模型输出的预设数量,y k,n为第k个模态分量对应的模型输出的第n个风速分量预测值,s k,ny k,n对应的实际值。In the method of the first aspect of the present application, the error coefficient of each modal component satisfies: , where/> is the error coefficient of the kth modal component, N is the preset number of model outputs corresponding to the kth modal component, yk ,n is the predicted value of the nth wind speed component output by the model corresponding to the kth modal component, and sk,n is the actual value corresponding to yk ,n .

为达上述目的,本申请第二方面实施例提出了一种储能辅助黑启动的风功率预测系统,包括:To achieve the above-mentioned purpose, the second embodiment of the present application proposes a wind power prediction system for energy storage-assisted black start, comprising:

获取模块,用于获取风速序列,所述风速序列包括当前采样点和多个历史采样点的风速测量值;An acquisition module, used to acquire a wind speed sequence, wherein the wind speed sequence includes wind speed measurement values of a current sampling point and multiple historical sampling points;

分解模块,用于利用CEEMD算法对所述风速序列进行分解获得多个模态分量;A decomposition module, used for decomposing the wind speed sequence by using a CEEMD algorithm to obtain multiple modal components;

误差计算模块,用于基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;An error calculation module is used to obtain a corresponding wind speed component prediction value based on each modal component using the first recurrent neural network model, and obtain an error coefficient corresponding to each modal component based on each wind speed component prediction value and each modal component;

重构模块,用于基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量;A reconstruction module, used for reconstructing all modal components to obtain multiple complementary modal components based on the positive and negative properties of each error coefficient;

预测模块,用于基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;A prediction module, used for obtaining a corresponding target wind speed component prediction value using a second recurrent neural network model based on each complementary modal component;

控制模块,用于基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用所述风功率预测值参与储能辅助黑启动。The control module is used to obtain the wind speed prediction value of the next sampling point of the current sampling point based on the prediction value of each target wind speed component, and then use the wind speed prediction value to obtain the wind power prediction value, so as to use the wind power prediction value to participate in energy storage assisted black start when a power outage occurs.

为达上述目的,本申请第三方面实施例提出了一种电子设备,包括:处理器,以及与所述处理器通信连接的存储器;所述存储器存储计算机执行指令;所述处理器执行所述存储器存储的计算机执行指令,以实现本申请第一方面提出的方法。To achieve the above-mentioned purpose, the third aspect embodiment of the present application proposes an electronic device, comprising: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method proposed in the first aspect of the present application.

为达上述目的,本申请第四方面实施例提出了一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,所述计算机执行指令被处理器执行时用于实现本申请第一方面提出的方法。To achieve the above-mentioned purpose, the fourth aspect embodiment of the present application proposes a computer-readable storage medium, in which computer execution instructions are stored. When the computer execution instructions are executed by a processor, they are used to implement the method proposed in the first aspect of the present application.

本申请提供的储能辅助黑启动的风功率预测方法、系统、电子设备及存储介质,通过获取风速序列,风速序列包括当前采样点和多个历史采样点的风速测量值;利用CEEMD算法对风速序列进行分解获得多个模态分量;基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量;基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用风功率预测值参与储能辅助黑启动。在这种情况下,在利用CEEMD算法获得多个模态分量后,还利用循环神经网络模型得到风速分量预测值进而获得各模态分量对应的误差系数,通过各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量,在利用互补模态分量得到当前采样点的下一采样点的风速预测值,避免像现有技术那样仅利用传统CEEMD算法分解得到分量来进行预测,降低了各分量的误差,提高了风速预测的精度,进而提高了风功率的预测精度。The wind power prediction method, system, electronic device and storage medium for energy storage assisted black start provided in the present application obtain a wind speed sequence, which includes wind speed measurement values of a current sampling point and multiple historical sampling points; decompose the wind speed sequence using a CEEMD algorithm to obtain multiple modal components; obtain a corresponding wind speed component prediction value based on each modal component using a first recurrent neural network model, and obtain an error coefficient corresponding to each modal component based on each wind speed component prediction value and each modal component; reconstruct all modal components based on the positive and negative properties of each error coefficient to obtain multiple complementary modal components; obtain a corresponding target wind speed component prediction value based on each complementary modal component using a second recurrent neural network model; obtain a wind speed prediction value of the next sampling point of the current sampling point based on each target wind speed component prediction value, and then use the wind speed prediction value to obtain a wind power prediction value, so that the wind power prediction value can be used to participate in the energy storage assisted black start when a power outage occurs. In this case, after using the CEEMD algorithm to obtain multiple modal components, the recurrent neural network model is also used to obtain the predicted value of the wind speed component and then obtain the error coefficient corresponding to each modal component. Through the positivity and negativity of each error coefficient, all modal components are reconstructed to obtain multiple complementary modal components. The complementary modal components are used to obtain the predicted wind speed value of the next sampling point of the current sampling point, avoiding the use of only the traditional CEEMD algorithm to decompose and obtain components for prediction as in the prior art, reducing the errors of each component, improving the accuracy of wind speed prediction, and thus improving the prediction accuracy of wind power.

本申请附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

本申请上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

图1为本申请实施例所提供的一种储能辅助黑启动的风功率预测方法的流程示意图;FIG1 is a schematic flow chart of a method for predicting wind power during energy storage-assisted black start provided in an embodiment of the present application;

图2为本申请实施例所提供的风速预测过程的具体流程示意图;FIG2 is a schematic diagram of a specific flow chart of a wind speed prediction process provided in an embodiment of the present application;

图3为本申请实施例所提供的一种储能辅助黑启动的风功率预测系统的框图。FIG3 is a block diagram of a wind power prediction system for energy storage-assisted black start provided in an embodiment of the present application.

具体实施方式Detailed ways

下面详细描述本申请的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。The embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

下面参考附图描述本申请实施例的储能辅助黑启动的风功率预测方法和系统。The following describes the wind power prediction method and system for energy storage-assisted black start in embodiments of the present application with reference to the accompanying drawings.

本申请实施例提供了一种储能辅助黑启动的风功率预测方法,以提高风功率的预测精度。本申请的实施例先预测风速,再利用预测风速得到所需的风功率预测值。具体过程如下。The embodiment of the present application provides a wind power prediction method for energy storage-assisted black start to improve the prediction accuracy of wind power. The embodiment of the present application first predicts the wind speed, and then uses the predicted wind speed to obtain the required wind power prediction value. The specific process is as follows.

图1为本申请实施例所提供的一种储能辅助黑启动的风功率预测方法的流程示意图。图2为本申请实施例所提供的风速预测过程的具体流程示意图。Fig. 1 is a schematic flow chart of a method for predicting wind power during energy storage-assisted black start provided in an embodiment of the present application. Fig. 2 is a schematic flow chart of a specific wind speed prediction process provided in an embodiment of the present application.

如图1所示,该储能辅助黑启动的风功率预测方法包括以下步骤:As shown in FIG1 , the energy storage-assisted black start wind power prediction method includes the following steps:

步骤S101,获取风速序列,风速序列包括当前采样点和多个历史采样点的风速测量值。Step S101, obtaining a wind speed sequence, where the wind speed sequence includes wind speed measurement values of a current sampling point and multiple historical sampling points.

具体地,在步骤S101中,设置采样点数量为M,M个采样点包括当前采样点和M-1个历史采样点,获取M个采样点处采集到的风速测量值得到所需的风速序列(也称原始风速序列)。换言之风速序列中风速测量值数量等于M。Specifically, in step S101, the number of sampling points is set to M, including the current sampling point and M-1 historical sampling points, and the wind speed measurement values collected at the M sampling points are obtained to obtain the required wind speed sequence (also called the original wind speed sequence). In other words, the number of wind speed measurement values in the wind speed sequence is equal to M.

步骤S102,利用CEEMD算法对风速序列进行分解获得多个模态分量。Step S102: Decompose the wind speed sequence using the CEEMD algorithm to obtain multiple modal components.

易于理解地,CEEMD(Complementary Ensemble Empirical Mode Decomposition,互补集合经验模态分解)算法是一种自适应信号处理方法。其利用1组符号相反的辅助白噪声和原始信号,得到1组新的正负混合序列,再对每个序列进行EMD(Empirical ModeDecomposition,经验模态分解)分解,通过多次分解最终将原始信号(如风速序列)分解为有限个数量的模态分量(Intrinsic Mode Function,IMF)。在步骤S102中,利用CEEMD算法对风速序列进行分解得到的模态分量数量为K,其中第k个模态分量表示为IMFk,k取1~K。It is easy to understand that the CEEMD (Complementary Ensemble Empirical Mode Decomposition) algorithm is an adaptive signal processing method. It uses a group of auxiliary white noises with opposite signs and the original signal to obtain a group of new positive and negative mixed sequences, and then performs EMD (Empirical Mode Decomposition) decomposition on each sequence. Through multiple decompositions, the original signal (such as the wind speed sequence) is finally decomposed into a finite number of modal components (Intrinsic Mode Function, IMF). In step S102, the number of modal components obtained by decomposing the wind speed sequence using the CEEMD algorithm is K, where the kth modal component is represented as IMF k , and k is 1~K.

由于风速序列包括M个风速测量值,故步骤S102中分解得到的每个模态分量可以看作是由风速序列的所有采样点的风速测量值分解得到的风速分量组成。风速分量数量等于风速序列的采样点数量。即每个模态分量的序列长度等于M,模态分量中的每个子序列为一个风速分量,每个模态分量对应的采样点与获取的风速序列一致,每个风速分量由对应采样点处的风速测量值分解得到。Since the wind speed sequence includes M wind speed measurement values, each modal component decomposed in step S102 can be regarded as a wind speed component decomposed from the wind speed measurement values of all sampling points of the wind speed sequence. The number of wind speed components is equal to the number of sampling points of the wind speed sequence. That is, the sequence length of each modal component is equal to M, each subsequence in the modal component is a wind speed component, the sampling point corresponding to each modal component is consistent with the acquired wind speed sequence, and each wind speed component is decomposed from the wind speed measurement value at the corresponding sampling point.

步骤S103,基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数。Step S103, based on each modal component, using the first recurrent neural network model to obtain the corresponding wind speed component prediction value, and based on each wind speed component prediction value and each modal component, obtaining the error coefficient corresponding to each modal component.

具体地,在步骤S103中,设置第一循环神经网络模型的参数,以使第一循环神经网络模型输出预设数量的风速分量预测值,模型输出的预设数量的风速分量预测值对应的实际值为输入的模态分量中当前采样点及之前的相同数量的风速分量;基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数,包括:针对任一模态分量,将该模态分量输入第一循环神经网络模型获得对应的预设数量的风速分量预测值;基于预设数量的风速分量预测值和对应的实际值获得该模态分量的误差系数,进而得到各模态分量的误差系数。Specifically, in step S103, the parameters of the first recurrent neural network model are set so that the first recurrent neural network model outputs a preset number of wind speed component prediction values, and the actual values corresponding to the preset number of wind speed component prediction values output by the model are the same number of wind speed components in the input modal components at the current sampling point and before; based on each modal component, the corresponding wind speed component prediction value is obtained using the first recurrent neural network model, and based on each wind speed component prediction value and each modal component, the error coefficient corresponding to each modal component is obtained, including: for any modal component, the modal component is input into the first recurrent neural network model to obtain the corresponding preset number of wind speed component prediction values; based on the preset number of wind speed component prediction values and the corresponding actual values, the error coefficient of the modal component is obtained, and then the error coefficient of each modal component is obtained.

以第一个模态分量IMF1、预设数量N,模态分量的序列长度等于M,其中M>N,为例,模型输入为第一个模态分量IMF1,第一个模态分量IMF1包括当前采样点和M-1个历史采样点下的M个风速分量,将该模态分量输入第一循环神经网络模型获得N个风速分量预测值,N个风速分量预测值对应的实际值为第一个模态分量IMF1中当前采样点及之前的N-1个历史采样点下的风速分量,基于N个风速分量预测值和对应的N个风速分量获得第一个模态分量IMF1的误差系数Taking the first modal component IMF 1 , the preset number N, and the sequence length of the modal component equal to M, where M>N, as an example, the model input is the first modal component IMF 1 , and the first modal component IMF 1 includes M wind speed components at the current sampling point and M-1 historical sampling points. The modal component is input into the first recurrent neural network model to obtain N wind speed component prediction values. The actual values corresponding to the N wind speed component prediction values are the wind speed components at the current sampling point and the previous N-1 historical sampling points in the first modal component IMF 1. Based on the N wind speed component prediction values and the corresponding N wind speed components, the error coefficient of the first modal component IMF 1 is obtained. .

在步骤S103中,各模态分量的误差系数满足:,其中/>为第k个模态分量的误差系数,N为第k个模态分量对应的模型输出的预设数量,y k,n为第k个模态分量对应的模型输出的第n个风速分量预测值,n取1~N。s k,ny k,n对应的实际值。In step S103, the error coefficient of each modal component satisfies: , where/> is the error coefficient of the kth modal component, N is the preset number of model outputs corresponding to the kth modal component, y k,n is the predicted value of the nth wind speed component output by the model corresponding to the kth modal component, and n ranges from 1 to N. s k,n is the actual value corresponding to y k,n .

在步骤S103中第一循环神经网络模型采用门控循环单元网络(gated recurrentunit,GRU)模型。易于理解地,GRU网络是LSTM(Long Short-Term Memory,长短期记忆网络)的改进模型,通过将遗忘门和输入门集成为更新门,一定程度上减少了网络的训练参数,同时又能保证对有效信息的记忆。更新门和重置门状态分别为z tr t,输入向量为x t,隐层状态为h t,数学表达式为:In step S103, the first recurrent neural network model adopts a gated recurrent unit (GRU) model. It is easy to understand that the GRU network is an improved model of the LSTM (Long Short-Term Memory). By integrating the forget gate and the input gate into an update gate, the training parameters of the network are reduced to a certain extent, while ensuring the memory of effective information. The update gate and reset gate states are z t and r t respectively, the input vector is x t , the hidden layer state is h t , and the mathematical expression is:

式中,t时刻输入状态与上一时刻隐层状态/>的过程向量;/>表示sigmoid函数;/>、/>分别表示输入向量对更新门和重置门的权重矩阵,/>、/>分别表示隐层状态对更新门和重置门的权重矩阵;/>表示隐层状态对过程向量的权重矩阵;W表示输入向量对过程向量的权重矩阵;/>表示输出权重矩阵,b为偏置;/>为t时刻的输出。In the formula, is the input state at time t and the hidden state at the previous time/> The process vector of Represents the sigmoid function; /> 、/> Represent the weight matrices of the input vector to the update gate and the reset gate, respectively./> 、/> Respectively represent the weight matrices of the hidden layer state to the update gate and the reset gate;/> represents the weight matrix of hidden layer state to process vector; W represents the weight matrix of input vector to process vector;/> represents the output weight matrix, b is the bias; /> is the output at time t.

如图2所示,在互补序列的建立阶段,原始风速序列经CEEMD分解,得到K个模态分量IMF(IMF1,...,IMFK),将K个IMF分量分别放入GRU模型进行预测,利用得到的风速分量预测值和对应的实际值获得各模态分量的误差系数(即,...,/>)。As shown in Figure 2, in the stage of establishing the complementary sequence, the original wind speed sequence is decomposed by CEEMD to obtain K modal components IMF (IMF 1 ,...,IMF K ). The K IMF components are respectively put into the GRU model for prediction, and the error coefficient of each modal component is obtained by using the predicted value of the wind speed component and the corresponding actual value (i.e. ,...,/> ).

步骤S104,基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量。Step S104: reconstruct all modal components based on the positive and negative properties of each error coefficient to obtain a plurality of complementary modal components.

具体地,在步骤S104中,基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量,包括:基于各模态分量的误差系数的正负性,将所有模态分量划分为非负数组和负数组;按误差系数,将非负数组中模态分量从大到小进行排列得到目标非负数组,将负数组中模态分量从小到大进行排列得到目标负数组;将目标非负数组与目标负数组中对应位置的模态分量相加从而得到多个互补模态分量。Specifically, in step S104, based on the positive or negative property of each error coefficient, all modal components are reconstructed to obtain a plurality of complementary modal components, including: based on the positive or negative property of the error coefficient of each modal component, all modal components are divided into a non-negative array and a negative array; according to the error coefficient, the modal components in the non-negative array are arranged from large to small to obtain a target non-negative array, and the modal components in the negative array are arranged from small to large to obtain a target negative array; and the modal components at corresponding positions in the target non-negative array and the target negative array are added to obtain a plurality of complementary modal components.

以非负数组中模态分量的数量为L,负数组中模态分量的数量为K-L为例,将K个误差系数按照正负号分成两组,其中L个误差系数大于等于0,K-L个误差系数小于0,然后将大于等于0的误差系数对应的模态分量,按照误差系数从大到小进行排列,得到目标非负数组[IMFmax-L+1,…,IMFmax-2,IMFmax-1,IMFmax],将小于0的误差系数对应的模态分量,按照误差系数从小到大进行排列,得到目标负数组[IMFmin-K+L+1,…,IMFmin-2,IMFmin-1,IMFmin],然后将这两组对应位置的IMF两两相加,直到序列长度较小的一组加完为止,形成多个互补模态分量[H1,H2,…,Hmax(L,K-L)](参见图2)。其中,两两相加时,从各组起始位置开始相加(即从IMFmax-L+1+IMFmin-K+L+1开始相加),直到序列长度较小的一组加完为止,序列长度较大的一组中剩余的模态分量直接保留。Taking the case where the number of modal components in the non-negative array is L and the number of modal components in the negative array is KL as an example, the K error coefficients are divided into two groups according to their positive and negative signs, where L error coefficients are greater than or equal to 0 and KL error coefficients are less than 0. Then, the modal components corresponding to the error coefficients greater than or equal to 0 are arranged from large to small according to the error coefficients to obtain the target non-negative array [IMF max-L+1 ,…,IMF max-2 ,IMF max-1 ,IMF max ]. The modal components corresponding to the error coefficients less than 0 are arranged from small to large according to the error coefficients to obtain the target negative array [IMF min-K+L+1 ,…,IMF min-2 ,IMF min-1 ,IMF min ]. Then, the IMFs at corresponding positions in these two groups are added in pairs until the group with the smaller sequence length is added, so as to form multiple complementary modal components [ H1 , H2 ,…, Hmax(L,KL) ] (see Figure 2). When adding two by two, the addition starts from the starting position of each group (that is, starting from IMF max-L+1 +IMF min-K+L+1 ) until the group with a smaller sequence length is added, and the remaining modal components in the group with a larger sequence length are directly retained.

步骤S105,基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值。Step S105, obtaining the corresponding target wind speed component prediction value using the second recurrent neural network model based on each complementary modal component.

具体地,在步骤S105中,设置第二循环神经网络模型的参数,以使第二循环神经网络模型输出当前采样点的下一采样点的目标风速分量预测值;基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值,包括:针对任一互补模态分量,将该互补模态分量输入第二循环神经网络模型获得对应的当前采样点的下一采样点的目标风速分量预测值,进而得到各互补模态分量对应的目标风速分量预测值。Specifically, in step S105, the parameters of the second recurrent neural network model are set so that the second recurrent neural network model outputs a target wind speed component prediction value of the next sampling point of the current sampling point; based on each complementary modal component, the second recurrent neural network model is used to obtain the corresponding target wind speed component prediction value of the next sampling point of the current sampling point, including: for any complementary modal component, the complementary modal component is input into the second recurrent neural network model to obtain the corresponding target wind speed component prediction value of the next sampling point of the current sampling point, and then the target wind speed component prediction value corresponding to each complementary modal component is obtained.

在步骤S105中,第二循环神经网络模型可以采用GRU模型。如图2所示,在预测阶段将步骤S104得到的多个互补模态分量[H1,H2,…,Hmax(L,K-L)]分别输入GRU模型分别得到对应互补模态分量的目标风速分量预测值(即Y1 , Y2 ,…, Ymax(L,K-L) )。In step S105, the second recurrent neural network model may adopt a GRU model. As shown in FIG2, in the prediction stage, the multiple complementary modal components [ H1 , H2 , ..., Hmax (L, KL) ] obtained in step S104 are respectively input into the GRU model to obtain the target wind speed component prediction values (i.e., Y1 ' , Y2 ' , ..., Ymax (L, KL) ' ) of the corresponding complementary modal components.

步骤S106,基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用风功率预测值参与储能辅助黑启动。Step S106, based on the predicted values of each target wind speed component, the predicted wind speed value of the next sampling point of the current sampling point is obtained, and then the predicted wind speed value is used to obtain the predicted wind power value, so as to use the predicted wind power value to participate in energy storage assisted black start when a power outage occurs.

具体地,在步骤S106中,基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,包括:将各目标风速分量预测值进行求和从而得到当前采样点的下一采样点的风速预测值。如图2所示,将所有目标风速分量预测值(即Y1 , Y2 ,…, Ymax(L,K-L) )求和从而得到风速预测最终结果,该最终结果即为当前采样点的下一采样点的风速预测值。Specifically, in step S106, the wind speed prediction value of the next sampling point of the current sampling point is obtained based on each target wind speed component prediction value, including: summing each target wind speed component prediction value to obtain the wind speed prediction value of the next sampling point of the current sampling point. As shown in FIG2, all target wind speed component prediction values (i.e., Y 1 ' , Y 2 ' , ..., Y max(L,KL) ' ) are summed to obtain the final result of wind speed prediction, which is the wind speed prediction value of the next sampling point of the current sampling point.

在步骤S106中,利用风速预测值得到风功率预测值包括:基于风速预测值、风机额定输出功率、额定风速、切入风速和切出风速获得风功率预测值。其中,风速与风功率(即风电机组的发电功率)的关系式满足:In step S106, using the wind speed prediction value to obtain the wind power prediction value includes: obtaining the wind power prediction value based on the wind speed prediction value, the rated output power of the wind turbine, the rated wind speed, the cut-in wind speed and the cut-out wind speed. The relationship between wind speed and wind power (i.e., the power generation power of the wind turbine) satisfies:

其中,P r为风机额定输出功率(kW);v为风机轮毂高度处的风速(m/s),v r为额定风速(m/s);v ci为切入风速(m/s),v co为切出风速(m/s)。Wherein, P r is the rated output power of the wind turbine (kW); v is the wind speed at the height of the wind turbine hub (m/s), v r is the rated wind speed (m/s); v ci is the cut-in wind speed (m/s), and v co is the cut-out wind speed (m/s).

在步骤S106中,将风速预测值看做风机轮毂高度处的风速v,按照风速与风功率的关系式获得风功率预测值P wIn step S106, the predicted wind speed value is regarded as the wind speed v at the height of the wind turbine hub, and the predicted wind power value P w is obtained according to the relationship between wind speed and wind power.

为了验证本申请的方法的效果,进行验证。其中,获取的实验数据(即原始风速序列)来自具有储能辅助黑启动能力的风电场。从2022年9月1日至2016年9月14日,每隔15分钟获取一次采样风速数据。对于每组案例,风速序列分为训练集和测试集。因此,选择7天的数据,并提供总共以15分钟为间隔共672个数据样本来训练预测模型;接下来的96个数据(对应于1天的数据)用于测试所提出模型的性能。In order to verify the effect of the method of the present application, a verification is performed. Among them, the experimental data obtained (i.e., the original wind speed sequence) comes from a wind farm with energy storage assisted black start capability. From September 1, 2022 to September 14, 2016, sampled wind speed data is obtained every 15 minutes. For each group of cases, the wind speed sequence is divided into a training set and a test set. Therefore, 7 days of data are selected, and a total of 672 data samples at intervals of 15 minutes are provided to train the prediction model; the next 96 data (corresponding to 1 day of data) are used to test the performance of the proposed model.

为了验证本申请所提出模型的性能,通过实验来进一步评估所提出的模型。其中GRU,CEEMD-GRU称为基准模型,用来与本申请所提出模型进行比较。在本申请的实验中,选择平均绝对百分比误差(Mean absolute percentage error, MAPE)及均方根误差(Rootmean square error, RMSE)为各模型的评价标准。本申请的实验中,实验结果如表1。In order to verify the performance of the model proposed in this application, experiments are conducted to further evaluate the proposed model. Among them, GRU and CEEMD-GRU are called benchmark models, which are used to compare with the model proposed in this application. In the experiments of this application, mean absolute percentage error (MAPE) and root mean square error (RMSE) are selected as the evaluation criteria for each model. In the experiments of this application, the experimental results are shown in Table 1.

表1 各模型的评价指标表Table 1 Evaluation index table of each model

如表1所示,当与其他所有的基准模型对比时,本申请提出的模型具有最小的误差系数。本申请的模型与CEEMD-GRU模型对比,2个评价指标的误差系数均大幅下降,证明了本申请所提出的方法能够有效的提高预测性能,将本申请的模型与GRU相比,2个评价指标的误差系数均大幅下降,证明数据预处理能够提高预测精度。As shown in Table 1, when compared with all other benchmark models, the model proposed in this application has the smallest error coefficient. Compared with the CEEMD-GRU model, the error coefficients of the two evaluation indicators of the model of this application have been greatly reduced, proving that the method proposed in this application can effectively improve the prediction performance. Compared with the GRU model of this application, the error coefficients of the two evaluation indicators have been greatly reduced, proving that data preprocessing can improve the prediction accuracy.

为了实现上述实施例,本申请还提出一种储能辅助黑启动的风功率预测系统。In order to implement the above embodiments, the present application also proposes a wind power prediction system for energy storage-assisted black start.

图3为本申请实施例所提供的一种储能辅助黑启动的风功率预测系统的框图。FIG3 is a block diagram of a wind power prediction system for energy storage-assisted black start provided in an embodiment of the present application.

如图3所示,该储能辅助黑启动的风功率预测系统包括获取模块11、分解模块12、误差计算模块13、重构模块14、预测模块15和控制模块16,其中:As shown in FIG3 , the energy storage assisted black start wind power prediction system includes an acquisition module 11, a decomposition module 12, an error calculation module 13, a reconstruction module 14, a prediction module 15 and a control module 16, wherein:

获取模块11,用于获取风速序列,风速序列包括当前采样点和多个历史采样点的风速测量值;An acquisition module 11 is used to acquire a wind speed sequence, where the wind speed sequence includes wind speed measurement values of a current sampling point and multiple historical sampling points;

分解模块12,用于利用CEEMD算法对风速序列进行分解获得多个模态分量;A decomposition module 12 is used to decompose the wind speed sequence using a CEEMD algorithm to obtain multiple modal components;

误差计算模块13,用于基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;The error calculation module 13 is used to obtain the corresponding wind speed component prediction value based on each modal component using the first recurrent neural network model, and obtain the error coefficient corresponding to each modal component based on each wind speed component prediction value and each modal component;

重构模块14,用于基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量;A reconstruction module 14, configured to reconstruct all modal components based on the positive and negative properties of each error coefficient to obtain a plurality of complementary modal components;

预测模块15,用于基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;A prediction module 15, configured to obtain a corresponding target wind speed component prediction value based on each complementary modal component using a second recurrent neural network model;

控制模块16,用于基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用风功率预测值参与储能辅助黑启动。The control module 16 is used to obtain the wind speed prediction value of the next sampling point of the current sampling point based on the predicted values of each target wind speed component, and then use the wind speed prediction value to obtain the wind power prediction value, so as to use the wind power prediction value to participate in energy storage assisted black start when a power outage occurs.

进一步地,在本申请实施例的一种可能的实现方式中,每个模态分量由风速序列的所有采样点的风速测量值分解得到的风速分量组成,风速分量数量等于风速序列的采样点数量,设置第一循环神经网络模型的参数,以使第一循环神经网络模型输出预设数量的风速分量预测值,模型输出的预设数量的风速分量预测值对应的实际值为输入的模态分量中当前采样点及之前的相同数量的风速分量;误差计算模块13,具体用于:针对任一模态分量,将该模态分量输入第一循环神经网络模型获得对应的预设数量的风速分量预测值;基于预设数量的风速分量预测值和对应的实际值获得该模态分量的误差系数,进而得到各模态分量的误差系数。Furthermore, in a possible implementation of an embodiment of the present application, each modal component is composed of wind speed components obtained by decomposing the wind speed measurement values of all sampling points of the wind speed sequence, the number of wind speed components is equal to the number of sampling points of the wind speed sequence, and the parameters of the first recurrent neural network model are set so that the first recurrent neural network model outputs a preset number of wind speed component prediction values, and the actual values corresponding to the preset number of wind speed component prediction values output by the model are the same number of wind speed components in the input modal components at the current sampling point and before; the error calculation module 13 is specifically used to: for any modal component, input the modal component into the first recurrent neural network model to obtain the corresponding preset number of wind speed component prediction values; obtain the error coefficient of the modal component based on the preset number of wind speed component prediction values and the corresponding actual values, and then obtain the error coefficient of each modal component.

进一步地,在本申请实施例的一种可能的实现方式中,第一循环神经网络模型和第二循环神经网络模型分别采用GRU模型。Furthermore, in a possible implementation of an embodiment of the present application, the first recurrent neural network model and the second recurrent neural network model respectively adopt GRU models.

进一步地,在本申请实施例的一种可能的实现方式中,误差计算模块13中,各模态分量的误差系数满足:,其中/>为第k个模态分量的误差系数,N为第k个模态分量对应的模型输出的预设数量,y k,n为第k个模态分量对应的模型输出的第n个风速分量预测值,s k,ny k,n对应的实际值。Further, in a possible implementation of the embodiment of the present application, in the error calculation module 13, the error coefficient of each modal component satisfies: , where/> is the error coefficient of the kth modal component, N is the preset number of model outputs corresponding to the kth modal component, yk ,n is the predicted value of the nth wind speed component output by the model corresponding to the kth modal component, and sk,n is the actual value corresponding to yk ,n .

进一步地,在本申请实施例的一种可能的实现方式中,重构模块14,具体用于:基于各模态分量的误差系数的正负性,将所有模态分量划分为非负数组和负数组;按误差系数,将非负数组中模态分量从大到小进行排列得到目标非负数组,将负数组中模态分量从小到大进行排列得到目标负数组;将目标非负数组与目标负数组中对应位置的模态分量相加从而得到多个互补模态分量。Furthermore, in a possible implementation of an embodiment of the present application, the reconstruction module 14 is specifically used to: divide all modal components into a non-negative array and a negative array based on the positive or negative nature of the error coefficient of each modal component; arrange the modal components in the non-negative array from large to small according to the error coefficient to obtain a target non-negative array, and arrange the modal components in the negative array from small to large to obtain a target negative array; add the modal components at corresponding positions in the target non-negative array and the target negative array to obtain multiple complementary modal components.

进一步地,在本申请实施例的一种可能的实现方式中,设置第二循环神经网络模型的参数,以使第二循环神经网络模型输出当前采样点的下一采样点的目标风速分量预测值;预测模块15,具体用于:针对任一互补模态分量,将该互补模态分量输入第二循环神经网络模型获得对应的当前采样点的下一采样点的目标风速分量预测值,进而得到各互补模态分量对应的目标风速分量预测值。Furthermore, in a possible implementation of an embodiment of the present application, the parameters of the second recurrent neural network model are set so that the second recurrent neural network model outputs a predicted value of the target wind speed component of the next sampling point of the current sampling point; the prediction module 15 is specifically used to: for any complementary modal component, input the complementary modal component into the second recurrent neural network model to obtain the corresponding predicted value of the target wind speed component of the next sampling point of the current sampling point, and then obtain the predicted value of the target wind speed component corresponding to each complementary modal component.

进一步地,在本申请实施例的一种可能的实现方式中,控制模块16,具体用于:将各目标风速分量预测值进行求和从而得到当前采样点的下一采样点的风速预测值。Furthermore, in a possible implementation of the embodiment of the present application, the control module 16 is specifically configured to: sum the predicted values of each target wind speed component to obtain a predicted wind speed value of a sampling point next to the current sampling point.

需要说明的是,前述对储能辅助黑启动的风功率预测方法实施例的解释说明也适用于该实施例的储能辅助黑启动的风功率预测系统,此处不再赘述。It should be noted that the aforementioned explanation of the embodiment of the wind power prediction method for energy storage-assisted black start is also applicable to the wind power prediction system for energy storage-assisted black start of this embodiment, and will not be repeated here.

本申请实施例中,通过获取风速序列,风速序列包括当前采样点和多个历史采样点的风速测量值;利用CEEMD算法对风速序列进行分解获得多个模态分量;基于各模态分量利用第一循环神经网络模型获得对应的风速分量预测值,基于各风速分量预测值和各模态分量获得各模态分量对应的误差系数;基于各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量;基于各互补模态分量利用第二循环神经网络模型获得对应的目标风速分量预测值;基于各目标风速分量预测值获得当前采样点的下一采样点的风速预测值,进而利用风速预测值得到风功率预测值,以在发生停电故障时利用风功率预测值参与储能辅助黑启动。在这种情况下,在利用CEEMD算法获得多个模态分量后,还利用循环神经网络模型得到风速分量预测值进而获得各模态分量对应的误差系数,通过各误差系数的正负性,对所有模态分量进行重构以得到多个互补模态分量,在利用互补模态分量得到当前采样点的下一采样点的风速预测值,避免像现有技术那样仅利用传统CEEMD算法分解得到分量来进行预测,降低了各分量的误差,提高了风速预测的精度,进而提高了风功率的预测精度。本申请的方法通过正负误差抵消法解决了传统数据分解技术(CEEMD)分解的子序列存在较大的误差,会对预测产生影响,而采用组合模型预测又会增加运行时间的问题,提高了风功率的预测精度。In an embodiment of the present application, a wind speed sequence is obtained, and the wind speed sequence includes wind speed measurement values of a current sampling point and multiple historical sampling points; the wind speed sequence is decomposed using a CEEMD algorithm to obtain multiple modal components; a corresponding wind speed component prediction value is obtained based on each modal component using a first recurrent neural network model, and an error coefficient corresponding to each modal component is obtained based on each wind speed component prediction value and each modal component; based on the positive and negative properties of each error coefficient, all modal components are reconstructed to obtain multiple complementary modal components; based on each complementary modal component, a corresponding target wind speed component prediction value is obtained using a second recurrent neural network model; based on each target wind speed component prediction value, a wind speed prediction value of the next sampling point of the current sampling point is obtained, and then a wind power prediction value is obtained using the wind speed prediction value, so that the wind power prediction value can be used to participate in energy storage-assisted black start when a power outage occurs. In this case, after using the CEEMD algorithm to obtain multiple modal components, the recurrent neural network model is also used to obtain the wind speed component prediction value and then obtain the error coefficient corresponding to each modal component. Through the positive and negative properties of each error coefficient, all modal components are reconstructed to obtain multiple complementary modal components. The wind speed prediction value of the next sampling point of the current sampling point is obtained by using the complementary modal components, avoiding the use of only the traditional CEEMD algorithm to decompose and obtain components for prediction as in the prior art, reducing the error of each component, improving the accuracy of wind speed prediction, and thus improving the prediction accuracy of wind power. The method of the present application solves the problem that the subsequence decomposed by the traditional data decomposition technology (CEEMD) has large errors, which will affect the prediction, and the use of a combined model prediction will increase the running time, thereby improving the prediction accuracy of wind power.

为了实现上述实施例,本申请还提出一种电子设备,包括:处理器,以及与处理器通信连接的存储器;存储器存储计算机执行指令;处理器执行存储器存储的计算机执行指令,以实现执行前述实施例所提供的方法。In order to implement the above embodiments, the present application also proposes an electronic device, comprising: a processor, and a memory communicatively connected to the processor; the memory stores computer-executable instructions; the processor executes the computer-executable instructions stored in the memory to implement the method provided by the above embodiments.

为了实现上述实施例,本申请还提出一种计算机可读存储介质,计算机可读存储介质中存储有计算机执行指令,计算机执行指令被处理器执行时用于实现前述实施例所提供的方法。In order to implement the above embodiments, the present application also proposes a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the methods provided by the above embodiments.

为了实现上述实施例,本申请还提出一种计算机程序产品,包括计算机程序,该计算机程序被处理器执行时实现前述实施例所提供的方法。In order to implement the above embodiments, the present application also proposes a computer program product, including a computer program, which implements the methods provided by the above embodiments when executed by a processor.

在前述各实施例描述中,参考术语“一个实施例”、“一些实施例”、 “示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。In the description of the aforementioned embodiments, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of this application, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本申请的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本申请的实施例所属技术领域的技术人员所理解。Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。The logic and/or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples (non-exhaustive list) of computer-readable media include the following: an electrical connection with one or more wires (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

应当理解,本申请的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如,如果用硬件来实现和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned 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, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

此外,在本申请各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

上述提到的存储介质可以是只读存储器,磁盘或光盘等。尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。The storage medium mentioned above may be a read-only memory, a disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims (9)

1. The wind power prediction method for energy storage auxiliary black start is characterized by comprising the following steps of:
Acquiring a wind speed sequence, wherein the wind speed sequence comprises wind speed measurement values of a current sampling point and a plurality of historical sampling points;
decomposing the wind speed sequence by using a CEEMD algorithm to obtain a plurality of modal components;
Obtaining a corresponding wind speed component predicted value by using a first cyclic neural network model based on each modal component, and obtaining an error coefficient corresponding to each modal component based on each wind speed component predicted value and each modal component;
Reconstructing all modal components based on the positive and negative of each error coefficient to obtain a plurality of complementary modal components;
Obtaining a corresponding target wind speed component predicted value by using a second cyclic neural network model based on each complementary modal component;
Obtaining a wind speed predicted value of a next sampling point of the current sampling point based on each target wind speed component predicted value, and further obtaining a wind power predicted value by utilizing the wind speed predicted value so as to participate in energy storage auxiliary black start when a power failure occurs;
wherein reconstructing all modal components based on the positive and negative of each error coefficient to obtain a plurality of complementary modal components comprises:
dividing all modal components into a non-negative array and a negative array based on the positive and negative of error coefficients of each modal component;
According to the error coefficient, arranging the modal components in the non-negative array from large to small to obtain a target non-negative array, and arranging the modal components in the negative array from small to large to obtain a target negative array;
And adding the modal components of the corresponding positions in the target non-negative array and the target negative array to obtain a plurality of complementary modal components.
2. The energy storage assisted black start wind power prediction method according to claim 1, wherein each modal component consists of wind speed components obtained by decomposing wind speed measurement values of all sampling points of the wind speed sequence, and the number of wind speed components is equal to the number of sampling points of the wind speed sequence;
Setting parameters of a first cyclic neural network model, so that the first cyclic neural network model outputs a preset number of wind speed component predicted values, and the actual values corresponding to the preset number of wind speed component predicted values output by the model are the current sampling points in the input modal components and the same number of wind speed components before the current sampling points;
The method for obtaining the corresponding wind speed component predicted value based on each modal component by using the first cyclic neural network model, obtaining the error coefficient corresponding to each modal component based on each wind speed component predicted value and each modal component comprises the following steps:
inputting any modal component into a first cyclic neural network model to obtain a corresponding preset number of wind speed component predicted values; and obtaining error coefficients of the modal components based on the predicted values of the wind speed components of the preset quantity and the corresponding actual values, and further obtaining the error coefficients of the modal components.
3. The energy storage-assisted black-start wind power prediction method according to claim 1, wherein parameters of the second cyclic neural network model are set so that the second cyclic neural network model outputs a target wind speed component predicted value of a sampling point next to the current sampling point;
The obtaining a corresponding target wind speed component predicted value based on each complementary modal component by using a second cyclic neural network model comprises the following steps:
And inputting the complementary modal component into a second cyclic neural network model for any complementary modal component to obtain a target wind speed component predicted value of a next sampling point of the corresponding current sampling point, and further obtaining the target wind speed component predicted value corresponding to each complementary modal component.
4. A method of energy storage assisted black start wind power prediction according to claim 3, wherein the obtaining a wind speed prediction value for a next sampling point to a current sampling point based on each target wind speed component prediction value comprises:
and summing the predicted values of the target wind speed components to obtain the predicted value of the wind speed of the next sampling point of the current sampling point.
5. The energy storage assisted black start wind power prediction method of claim 1, wherein the first and second recurrent neural network models each employ a GRU model.
6. The energy storage assisted black start wind power prediction method according to claim 2, wherein error coefficients of each modal component satisfy: wherein/> For the error coefficient of the kth modal component, N is the preset number of model outputs corresponding to the kth modal component, y k,n is the nth wind speed component predicted value of the model outputs corresponding to the kth modal component, and s k,n is the actual value corresponding to y k,n.
7. An energy storage assisted black start wind power prediction system, comprising:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring a wind speed sequence, and the wind speed sequence comprises wind speed measurement values of a current sampling point and a plurality of historical sampling points;
The decomposition module is used for decomposing the wind speed sequence by using a CEEMD algorithm to obtain a plurality of modal components;
the error calculation module is used for obtaining a corresponding wind speed component predicted value by utilizing the first cyclic neural network model based on each modal component and obtaining an error coefficient corresponding to each modal component based on each wind speed component predicted value and each modal component;
The reconstruction module is used for reconstructing all the modal components based on the positive and negative of each error coefficient to obtain a plurality of complementary modal components;
The prediction module is used for obtaining a corresponding target wind speed component predicted value by using a second cyclic neural network model based on each complementary modal component;
the control module is used for obtaining a wind speed predicted value of a next sampling point of the current sampling point based on each target wind speed component predicted value, and further obtaining a wind power predicted value by utilizing the wind speed predicted value so as to participate in energy storage auxiliary black start when a power failure occurs;
wherein reconstructing all modal components based on the positive and negative of each error coefficient to obtain a plurality of complementary modal components comprises:
dividing all modal components into a non-negative array and a negative array based on the positive and negative of error coefficients of each modal component;
According to the error coefficient, arranging the modal components in the non-negative array from large to small to obtain a target non-negative array, and arranging the modal components in the negative array from small to large to obtain a target negative array;
And adding the modal components of the corresponding positions in the target non-negative array and the target negative array to obtain a plurality of complementary modal components.
8. An electronic device, comprising: a processor, and a memory communicatively coupled to the processor;
The memory stores computer-executable instructions;
the processor executes computer-executable instructions stored in the memory to implement the method of any one of claims 1-6.
9. A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to carry out the method of any one of claims 1-6.
CN202410182403.6A 2024-02-19 2024-02-19 A wind power prediction method and system for energy storage assisted black start Active CN117748501B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410182403.6A CN117748501B (en) 2024-02-19 2024-02-19 A wind power prediction method and system for energy storage assisted black start

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410182403.6A CN117748501B (en) 2024-02-19 2024-02-19 A wind power prediction method and system for energy storage assisted black start

Publications (2)

Publication Number Publication Date
CN117748501A CN117748501A (en) 2024-03-22
CN117748501B true CN117748501B (en) 2024-05-07

Family

ID=90279703

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410182403.6A Active CN117748501B (en) 2024-02-19 2024-02-19 A wind power prediction method and system for energy storage assisted black start

Country Status (1)

Country Link
CN (1) CN117748501B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118311692B (en) * 2024-03-28 2024-08-20 西安热工研究院有限公司 Wind speed prediction method for energy storage assisted black start based on two-dimensional error compensation method
CN117972363B (en) * 2024-03-29 2024-06-07 西安热工研究院有限公司 Wind speed prediction method and system for energy storage auxiliary black start based on stability
CN117977583B (en) * 2024-04-01 2024-07-12 西安热工研究院有限公司 Energy storage auxiliary black start wind speed prediction method based on double-order error compensation
CN118523318B (en) * 2024-07-19 2024-10-18 西安热工研究院有限公司 A black start wind power prediction method based on super capacity energy storage

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016127643A (en) * 2014-12-26 2016-07-11 株式会社日立パワーソリューションズ Wind power generation amount prediction method
CN106203723A (en) * 2016-07-19 2016-12-07 河海大学 Wind power short-term interval prediction method based on RT reconstruct EEMD RVM built-up pattern
CN109118000A (en) * 2018-08-07 2019-01-01 广东工业大学 A kind of short-term wind speed forecasting method based on CEEMD-VMD-GA-ORELM model
CN112669168A (en) * 2020-12-15 2021-04-16 国网辽宁省电力有限公司阜新供电公司 Short-term wind power prediction method
WO2021164112A1 (en) * 2020-02-18 2021-08-26 山东大学 Frequency control method and system during using wind farm as black-start power source by means of optimal configuration of energy storage
CN113688567A (en) * 2021-08-10 2021-11-23 华北电力大学(保定) Two-stage optimal scheduling method of virtual power plant considering impact load
CN113783207A (en) * 2021-08-20 2021-12-10 清华大学 Hybrid energy storage method, device, electronic device and storage medium for suppressing wind power fluctuations
CN113988481A (en) * 2021-12-23 2022-01-28 南京鼐威欣信息技术有限公司 Wind power prediction method based on dynamic matrix prediction control
CN115549202A (en) * 2022-09-27 2022-12-30 国网辽宁省电力有限公司沈阳供电公司 Real-time electricity price based virtual power plant multi-energy coordination low-carbon optimization scheduling method and system and storage medium
CN116885711A (en) * 2023-07-20 2023-10-13 贵州电网有限责任公司 Wind power prediction method, device, equipment and readable storage medium
CN116960941A (en) * 2023-06-26 2023-10-27 华电电力科学研究院有限公司 Power climbing event prediction method, device, equipment and storage medium
CN117039895A (en) * 2023-10-09 2023-11-10 西安热工研究院有限公司 Wind power prediction method and system for energy storage auxiliary black start
CN117150705A (en) * 2023-01-31 2023-12-01 浙大城市学院 Marine wind farm short-term wind speed prediction method based on characteristic deconstruction
CN117498352A (en) * 2023-12-29 2024-02-02 西安热工研究院有限公司 Wind speed prediction method and device based on energy storage auxiliary black start capability

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016127643A (en) * 2014-12-26 2016-07-11 株式会社日立パワーソリューションズ Wind power generation amount prediction method
CN106203723A (en) * 2016-07-19 2016-12-07 河海大学 Wind power short-term interval prediction method based on RT reconstruct EEMD RVM built-up pattern
CN109118000A (en) * 2018-08-07 2019-01-01 广东工业大学 A kind of short-term wind speed forecasting method based on CEEMD-VMD-GA-ORELM model
WO2021164112A1 (en) * 2020-02-18 2021-08-26 山东大学 Frequency control method and system during using wind farm as black-start power source by means of optimal configuration of energy storage
CN112669168A (en) * 2020-12-15 2021-04-16 国网辽宁省电力有限公司阜新供电公司 Short-term wind power prediction method
CN113688567A (en) * 2021-08-10 2021-11-23 华北电力大学(保定) Two-stage optimal scheduling method of virtual power plant considering impact load
CN113783207A (en) * 2021-08-20 2021-12-10 清华大学 Hybrid energy storage method, device, electronic device and storage medium for suppressing wind power fluctuations
CN113988481A (en) * 2021-12-23 2022-01-28 南京鼐威欣信息技术有限公司 Wind power prediction method based on dynamic matrix prediction control
CN115549202A (en) * 2022-09-27 2022-12-30 国网辽宁省电力有限公司沈阳供电公司 Real-time electricity price based virtual power plant multi-energy coordination low-carbon optimization scheduling method and system and storage medium
CN117150705A (en) * 2023-01-31 2023-12-01 浙大城市学院 Marine wind farm short-term wind speed prediction method based on characteristic deconstruction
CN116960941A (en) * 2023-06-26 2023-10-27 华电电力科学研究院有限公司 Power climbing event prediction method, device, equipment and storage medium
CN116885711A (en) * 2023-07-20 2023-10-13 贵州电网有限责任公司 Wind power prediction method, device, equipment and readable storage medium
CN117039895A (en) * 2023-10-09 2023-11-10 西安热工研究院有限公司 Wind power prediction method and system for energy storage auxiliary black start
CN117498352A (en) * 2023-12-29 2024-02-02 西安热工研究院有限公司 Wind speed prediction method and device based on energy storage auxiliary black start capability

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Wind Speed Data Repairing Method Based on Bidirectional Prediction;Xincheng Shen et al.;IEEE;20211231;全文 *
基于特征变权的超短期风电功率预测;王晓东 等;太阳能学报;20230228;全文 *

Also Published As

Publication number Publication date
CN117748501A (en) 2024-03-22

Similar Documents

Publication Publication Date Title
CN117748501B (en) A wind power prediction method and system for energy storage assisted black start
CN117829202B (en) A wind speed prediction method and system for energy storage-assisted black start based on double errors
CN117744893B (en) A wind speed prediction method and system for energy storage assisted black start
CN116306798A (en) A method and system for ultra-short-term wind speed prediction
CN112906757B (en) Model training method, optical power prediction method, device, equipment and storage medium
CN111709112B (en) Offshore wind power operation simulation method, device and storage medium
CN118013235B (en) Energy storage auxiliary black start wind speed prediction method based on cyclic integrated error compensation
CN112669168B (en) Short-term prediction method for wind power
CN109872248B (en) A kind of wind farm cluster output calculation method and system
CN117708547B (en) Method and system for processing vibration signals of steam turbine unit
CN117787509B (en) A wind speed prediction method and system for energy storage assisted black start
CN117390967B (en) Offshore wind power short-term power prediction method, device, equipment and storage medium
CN118569424A (en) A short-term wind power integrated prediction method and device
CN117972363B (en) Wind speed prediction method and system for energy storage auxiliary black start based on stability
CN118013236B (en) Wind speed prediction method and system for energy storage assisted black start based on error compensation method
CN117977583B (en) Energy storage auxiliary black start wind speed prediction method based on double-order error compensation
CN118311692B (en) Wind speed prediction method for energy storage assisted black start based on two-dimensional error compensation method
CN116484743A (en) Method and device for evaluating the impact of wind power grid-connected on power system reliability
CN115409291B (en) Wind power prediction method and system for correcting wind speed
CN116611011A (en) A method, system, terminal and medium for predicting transient overvoltage amplitude
CN117728504B (en) A black start system and method for a diesel-generator combined combustion engine
Xu et al. A new approach for fast reliability evaluation of composite power system considering wind farm
CN114399402B (en) Method, system, equipment and medium for evaluating the power generation increase rate before and after unit transformation
Li et al. Pre-trained Multivariate Time Series Graph Neural Networks for Wind Power Forecasting
CN113408193B (en) New energy consumption capacity evaluation method and device, terminal device and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant