TW201635224A - Method of short-term wind power generation forecasting - Google Patents

Method of short-term wind power generation forecasting Download PDF

Info

Publication number
TW201635224A
TW201635224A TW104108689A TW104108689A TW201635224A TW 201635224 A TW201635224 A TW 201635224A TW 104108689 A TW104108689 A TW 104108689A TW 104108689 A TW104108689 A TW 104108689A TW 201635224 A TW201635224 A TW 201635224A
Authority
TW
Taiwan
Prior art keywords
power generation
wind power
layer
node
value
Prior art date
Application number
TW104108689A
Other languages
Chinese (zh)
Other versions
TWI540533B (en
Inventor
張文宇
Original Assignee
聖約翰科技大學
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 聖約翰科技大學 filed Critical 聖約翰科技大學
Priority to TW104108689A priority Critical patent/TWI540533B/en
Application granted granted Critical
Publication of TWI540533B publication Critical patent/TWI540533B/en
Publication of TW201635224A publication Critical patent/TW201635224A/en

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

An intelligent method of forecasting short-term wind power generation includes providing a persistence algorithm to forecast wind speed and wind power generation in the future according to the wind speed and wind power generation in the present; providing a radial basis function neural network algorithm to forecast wind power generation in the future; providing a back propagation neural network algorithm to forecast wind power generation in the future; providing a smart mathematical forecasting model; and providing an artificial bee colony algorithm to solve weighting coefficients of the smart mathematical forecasting model to obtain the optimal wind power generation forecasting value.

Description

智慧型短期電力發電量預測方法Intelligent short-term power generation forecasting method

本發明係有關一種短期電力發電量預測方法,尤指一種以人工蜂群演算法求解最佳化複合預測法的權重係數之智慧型短期電力發電量預測方法。The invention relates to a short-term power generation quantity prediction method, in particular to a smart short-term power generation quantity prediction method which uses an artificial bee colony algorithm to solve the weight coefficient of the optimized compound prediction method.

由於石化能源日漸減少,石油價格不斷提高,能源供給的日漸短絀已是不爭的事實。此外,由於環保議題逐漸受到重視,在2005年「京都議定書」實施後,全國二氧化碳的總排放量已受到管制,因而限制了石化燃料的使用。在傳統熱力發電機組增設不易的情形下,新的替代能源開發就是必須執行的方案。因此,各先進國家莫不以再生能源的應用技術開發為努力目標,而風力發電就是最具商業運轉價值的再生能源之一,尤其在用電量日增的台灣,風力發電的技術開發更是刻不容緩的。近年來風力發電機組已大量的設置在台灣的西部海岸,主要原因是風力發電具有高效率、高功率密度和低污染的特性。As petrochemical energy is declining and oil prices are rising, it is an indisputable fact that energy supply is becoming shorter. In addition, as environmental issues have gradually gained attention, after the implementation of the Kyoto Protocol in 2005, the total emissions of carbon dioxide in the country have been regulated, thus limiting the use of fossil fuels. In the case that the traditional heat generating units are not easy to add, new alternative energy development is a necessary solution. Therefore, all advanced countries do not aim at the development of renewable energy application technology, and wind power generation is one of the most valuable renewable energy sources. Especially in Taiwan, where electricity consumption is increasing, the development of wind power technology is urgent. of. In recent years, wind turbines have been installed on the western coast of Taiwan. The main reason is that wind power has high efficiency, high power density and low pollution.

雖然風力發電具有上述諸多優點,但在應用上仍有一項亟待克服的技術難題,就是風力發電量預測。因為風力發電機組的發電量會隨著氣象狀態和風速而變化,風力發電機組發電量的不可預期變化將導致電力系統備轉容量提高,而使運轉成本增加。電力系統調度人員必須預測風力發電機組發電量的變化以調度備轉容量與管理系統運轉。準確的風力發電量預測有助於提高自由化市場機制的設計、電網的即時運轉控制、併網標準的制定以及提升電力品質。由於風力發電已成為再生能源開發的重要途徑,因此風力發電量預測的相關技術開發也成為電力工程中極為重要的一環。Although wind power has many of the above advantages, there is still a technical problem to be overcome in application, which is the prediction of wind power generation. Because the wind turbine's power generation will change with meteorological conditions and wind speed, unpredictable changes in wind turbine power generation will lead to an increase in power system backup capacity and an increase in operating costs. Power system dispatchers must predict changes in wind turbine generation to schedule the reserve capacity and manage the system. Accurate wind power forecasting helps to improve the design of liberalized market mechanisms, the immediate operation control of the grid, the development of grid-connected standards, and the improvement of power quality. Since wind power has become an important way to develop renewable energy, the development of related technologies for wind power forecasting has become an extremely important part of power engineering.

因此,如何設計出一種智慧型短期電力發電量預測方法,提供具有極大彈性、能適用於各種類型的風力發電機組發電量預測,並且以人工蜂群演算法求解每種個別預測方法的權重係數,有效提高智慧型預測法的預測準確性,進而有效與正確的預測風力發電系統在不同季節狀態下的發電量,提供電力系統調度人員進行相關電力系統運轉控制時的重要工具,以達到提高電力系統穩定度、降低系統運轉成本目的,乃為本案發明人所欲行克服並加以解決的一大課題。Therefore, how to design a smart short-term power generation forecasting method, provide a large elasticity, can be applied to various types of wind turbine generation generation forecast, and use artificial bee colony algorithm to solve the weight coefficient of each individual prediction method, Effectively improve the prediction accuracy of the intelligent forecasting method, and effectively and correctly predict the power generation of the wind power generation system in different seasons, and provide an important tool for the power system dispatcher to control the operation of the relevant power system to improve the power system. The purpose of stability and lowering the operating cost of the system is a major issue that the inventors of the present invention have tried to overcome and solve.

本發明之一目的在於提供一種智慧型短期電力發電量預測方法,以克服習知技術的問題。該預測方法係包含:(a) 提供一持續法,根據目前時間點測得的風速與風力發電量,預測未來的時間點的預測風速與風力發電量:;其中,為目前時間點、()為未來時間點、為目前風速、為未來風速、為目前風力發電量以及為未來風力發電量;(b) 提供一徑向基底函數類神經網路演算法,利用包括一輸入層、一隱藏層以及一輸出層之網路架構,預測未來風力發電量;(c) 提供一倒傳遞類神經網路演算法,利用包括一輸入層、一隱藏層以及一輸出層之網路架構,預測未來風力發電量;(d) 根據該持續法、該徑向基底函數類神經網路演算法以及該倒傳遞類神經網路演算法,提供一智慧型數學預測模型;其中,;以及(t = 1, 2, …,L )為實際的發電量時間序列數據,M為個別預測方法數量,L為樣本數,(i = 1, 2, …,M ,t = 1, 2, …,L )係為第i種預測方法的預測值,=-為預測誤差,為第i種預測方法的權重係數,的估測值,而為智慧型預測法的預測值;及(e) 提供一人工蜂群演算法,求解該智慧型數學預測模型之權重係數,以獲得風力發電量最佳值;其中,It is an object of the present invention to provide a smart short-term power generation amount prediction method that overcomes the problems of the prior art. The prediction method includes: (a) providing a continuous method to predict the predicted wind speed and wind power generation at a future time point based on the wind speed and wind power generation measured at the current time point: , ;among them, For the current time, ( ) for the future time, For the current wind speed, For future wind speed, For current wind power generation as well (b) provide a radial basis function-like neural network algorithm that predicts future wind power generation using a network architecture that includes an input layer, a hidden layer, and an output layer; (c) provides a An inverse transfer-like neural network algorithm predicts future wind power generation using a network architecture including an input layer, a hidden layer, and an output layer; (d) according to the persistence method, the radial basis function-like neural network algorithm And the inverse transfer neural network algorithm provides a smart mathematical prediction model ;among them, , ;as well as ( t = 1, 2, ..., L ) is the actual power generation time series data, M is the number of individual prediction methods, and L is the number of samples. ( i = 1, 2, ..., M , t = 1, 2, ..., L ) is the predicted value of the i-th prediction method, = - For prediction errors, The weighting factor for the i-th prediction method, for Estimated value, and a predictive value for the intelligent predictive method; and (e) providing an artificial bee colony algorithm to solve the weighting coefficient of the intelligent mathematical predictive model To obtain the best value of wind power generation; .

為了能更進一步瞭解本發明為達成預定目的所採取之技術、手段及功效,請參閱以下有關本發明之詳細說明與附圖,相信本發明之目的、特徵與特點,當可由此得一深入且具體之瞭解,然而所附圖式僅提供參考與說明用,並非用來對本發明加以限制者。In order to further understand the technology, the means and the effect of the present invention in order to achieve the intended purpose, refer to the following detailed description of the invention and the accompanying drawings. The detailed description is to be understood as illustrative and not restrictive.

茲有關本發明之技術內容及詳細說明,配合圖式說明如下:The technical content and detailed description of the present invention are as follows:

請參閱圖1係為本發明智慧型短期風力發電系統發電量預測系統之系統架構圖。該發電量預測系統的建構首先由實際風力發電系統的發電資料建立訓練資料庫,而後經由訓練程序建立一套短期風力發電量預測系統。本發明智慧型短期風力發電系統發電量預測系統組合三種短期風力發電系統發電量預測系統:(1)持續法、(2)徑向基底函數類神經網路演算法、以及(3)倒傳遞類神經網路演算法,輸入信號先以三種短期風力發電系統發電量預測系統進行預測,再以複合式預測法之方程式計算出發電量預測值。其中,該短期風力發電量預測系統係以接收前20分鐘發電量、前10分鐘發電量、目前發電量以及風速預測值為輸入資料,並且利用人工蜂群演算法求解該風力發電機組發電量,以獲得10分鐘後發電量的預測值。Please refer to FIG. 1 , which is a system architecture diagram of a power generation quantity prediction system for a smart short-term wind power generation system of the present invention. The construction of the power generation quantity prediction system first establishes a training database from the power generation data of the actual wind power generation system, and then establishes a short-term wind power generation quantity prediction system through the training program. The intelligent short-term wind power generation system power generation quantity prediction system combines three short-term wind power generation system power generation quantity prediction systems: (1) continuous method, (2) radial basis function-like neural network algorithm, and (3) reverse-transfer-like nerve In the network algorithm, the input signal is first predicted by three short-term wind power generation system prediction systems, and the predicted value of the starting electricity is calculated by the equation of the composite prediction method. The short-term wind power generation forecasting system inputs the data of the first 20 minutes of power generation, the first 10 minutes of power generation, the current power generation amount, and the wind speed prediction value, and uses an artificial bee colony algorithm to solve the wind power generation amount. Obtain a predicted value of power generation after 10 minutes.

智慧型預測法的組合理論基礎是以數個不同的個別預測方法共同解決特定的預測問題,每一個別預測方法具有特定權重係數,智慧型預測的預測值則是將這些個別預測方法的權重和(weighting sum)。智慧型預測法的數學模型可以表示如下:The combined theoretical basis of intelligent predictive method is to solve specific prediction problems by several different individual prediction methods. Each individual prediction method has a specific weight coefficient, and the intelligent prediction prediction value is the weight of these individual prediction methods. (weighting sum). The mathematical model of the intelligent prediction method can be expressed as follows:

(1) (1)

(2) (2)

其中,(t = 1, 2, …,L )為實際的發電量時間序列數據,M為個別預測方法數量,L為樣本數,(i = 1, 2, …,M ,t = 1, 2, …,L )係為第i種預測方法的預測值,=-為預測誤差,為第i種預測方法的權重係數,的估測值,而為智慧型預測法的預測值。among them, ( t = 1, 2, ..., L ) is the actual power generation time series data, M is the number of individual prediction methods, and L is the number of samples. ( i = 1, 2, ..., M , t = 1, 2, ..., L ) is the predicted value of the i-th prediction method, = - For prediction errors, The weighting factor for the i-th prediction method, for Estimated value, and The predicted value for the intelligent forecasting method.

在智慧型預測法的結構中對於每種個別預測方法中權重係數的測定是最重要的步驟,權重係數可以經由求解複合預測法的絕對誤差極小值的最佳化問題來實現。最佳化問題的數學模型可以表示如下:In the structure of the intelligent prediction method, the determination of the weight coefficient in each individual prediction method is the most important step, and the weight coefficient can be realized by solving the optimization problem of the absolute error minimum value of the composite prediction method. The mathematical model of the optimization problem can be expressed as follows:

(3) (3)

以下,將對本發明所採用之演算法加以詳細說明。Hereinafter, the algorithm employed in the present invention will be described in detail.

1、人工蜂群演算法(artificial bee colony, ABC)1. Artificial bee colony (ABC)

本發明將以人工蜂群演算法求解(3)式,人工蜂群演算法是依據蜂群覓食之行為模式,達成最佳化演算法之研發,這種生物群體智慧演算法之演算流程兼具開發程序及探索機制特性,不僅避免求解過程陷入局部解,並可提升求解最佳化問題之運算效能。本發明將以人工蜂群演算法決定三種個別預測方法的權重係數組合,再以(2)式進行短期風力發電系統發電量預測。The invention will solve the (3) formula by the artificial bee colony algorithm, and the artificial bee colony algorithm is based on the bee group foraging behavior mode, and the optimization algorithm is developed, and the calculation process of the bio-group wisdom algorithm is The development program and the exploration mechanism feature not only avoid the solution process falling into the local solution, but also improve the computational efficiency of solving the optimization problem. The invention will determine the weight coefficient combination of the three individual prediction methods by the artificial bee colony algorithm, and then predict the power generation of the short-term wind power generation system by the formula (2).

依據蜂群覓食之過程,工蜂可分為採集蜂、待命蜂與偵察蜂等三類角色。在覓食初期,因工蜂群尚無任何食物源資訊,工蜂便以偵察蜂的角色隨機搜尋,一旦找到食物源後,即以採集蜂的角色前往採集該食物源之花蜜,直至該食物源之花蜜採集完畢,再轉化為待命蜂的角色或偵察蜂的角色,以便後續採蜜工作之進行。當工蜂放棄該食物源後,若轉化為待命蜂的角色時,便飛回蜂巢休息及等待下次偵察蜂分享食物源訊息;反之,若轉化為偵察蜂的角色時,則持續隨機搜尋新的食物源位置。According to the process of bee colony feeding, worker bees can be divided into three types of characters: collecting bees, standby bees and scout bees. In the early stage of foraging, because the worker bee group does not have any food source information, the worker bees randomly search for the role of the scout bee. Once the food source is found, the bee is collected to collect the nectar of the food source until the food source After the nectar is collected, it is converted into the role of the standby bee or the role of the scout bee, so that the subsequent honey collecting work can be carried out. When the worker bees give up the food source, if they are converted into the role of the standby bee, they fly back to the hive to rest and wait for the next scout to share the food source message; otherwise, if they are converted into the role of the scout bee, they continue to randomly search for new ones. Food source location.

本發明將人工蜂群演算法用於求解個別預測方法的權重係數組合問題,食物源位置就是個別預測方法的權重係數組合,食物源所具有之花蜜量,可視為最佳化問題之目標函數值亦即複合預測法的絕對誤差。於人工蜂群演算法開始前,需先輸入個別預測方法相關資訊,於人工蜂群演算法的計算過程中,參考實際工蜂搜尋花蜜之行為模式,可將其概分為開發程序及探索程序,且若將此兩種程序合併運行,不僅可執行局部搜尋工作,同時有助於全域解的搜尋。The artificial bee colony algorithm is used to solve the weight coefficient combination problem of the individual prediction methods. The food source position is the weight coefficient combination of the individual prediction methods, and the nectar quantity of the food source can be regarded as the target function value of the optimization problem. That is, the absolute error of the composite prediction method. Before the start of the artificial bee colony algorithm, it is necessary to input the relevant information of the individual prediction methods. In the calculation process of the artificial bee colony algorithm, the actual worker bee is searched for the behavior pattern of the nectar, which can be divided into a development program and an exploration program. And if this program is combined and run, not only can the local search work be performed, but also the search for the global solution.

人工蜂群演算法先隨機產生N 個初始食物源位置與N 隻採集蜂,且每個食物源均對映至一隻採集蜂,食物源的位置就是求解最佳化問題的自變數,在本發明中即為個別預測方法的權重係數組合,而各食物源所蘊藏之花蜜量,則是求解最佳化問題的目標函數值,在本發明中即為複合預測法的絕對誤差。食物源位置初始化完成後,即可執行開發程序與探索程序。The artificial bee colony algorithm randomly generates N initial food source positions and N collecting bees, and each food source is mapped to a collecting bee. The position of the food source is the self-changing parameter for solving the optimization problem. In the invention, the weight coefficient combination of the individual prediction methods, and the amount of nectar contained in each food source is the target function value for solving the optimization problem, which is the absolute error of the composite prediction method in the present invention. Once the food source location is initialized, the development and exploration programs can be executed.

在實際蜂群覓食過程中,當食物源花蜜量不足或已採集完時,採集蜂將轉變為偵察蜂,尋找新的食物源位置。本發明於人工蜂群演算法中,先設定採集蜂採集某食物源花蜜之次數,以評估該食物源之花蜜量是否足夠,亦即目標函數值(複合預測法的絕對誤差),經由設定程序及演算,若目標函數值無法獲得改善時,即放棄該組解,隨機產生另一組新解,以避免演算過程過早收斂或陷入局部解。人工蜂群演算法之重點乃在於採集蜂探索富含花蜜量之食物源,並將此食物資訊分享至其他工蜂,以期在該食物源附近找尋花蜜蘊藏量更豐富之食物源,進而增強區域搜尋能力,若該食物源之花蜜量不足 (目標函數值無法降低),則放棄該食物源並由偵查蜂隨機搜尋新食物源,以增強演算開發新可行解。During the actual bee colony foraging, when the amount of nectar in the food source is insufficient or has been collected, the collecting bee will be turned into a scout bee to find a new food source location. In the artificial bee colony algorithm, the first step is to set the number of nectar of the food source to collect whether the nectar of the food source is sufficient, that is, the target function value (absolute error of the composite prediction method), through the setting procedure. And calculus, if the target function 无法 value can not be improved, the group solution is abandoned, and another set of new solutions is randomly generated to avoid premature convergence or partial solution to the calculation process. The focus of the artificial bee colony algorithm is to collect bees to explore food sources rich in nectar and to share this food information with other worker bees, in order to find a more abundant source of nectar in the vicinity of the food source, thereby enhancing regional search. Ability, if the amount of nectar in the food source is insufficient (the target function cannot be reduced), the food source is abandoned and the scout is randomly searched for new food sources to enhance the calculation of new feasible solutions.

人工蜂群演算法執行開發程序時,待命蜂開發新食物源,即將初始食物源位置更新,以開發其餘富含花蜜之食物源,同時避免原食物源之花蜜蘊藏量耗竭。若待命蜂尋得新食物源後,比較新食物源與原食物源位置之目標函數值,若新食物源較原位置之目標函數值佳時,則以新食物源取代舊食物源位置。執行探索程序時,將召集偵查蜂前往各處探訪,並隨機搜索可能的食物源(新可行解),規劃偵查蜂探索新食物源之行為,並以偵查蜂數目作為決策依據,淘汰目標函數值較差之食物源,並以新食物源取代。直到疊代次數到達程式所設定之疊代次數上限時,即停止人工蜂群演算法演算並輸出最佳解,在本發明中最佳解即為個別預測方法的最佳權重係數組合。When the artificial bee colony algorithm performs the development process, the standby bee develops a new food source, and the initial food source location is updated to develop the remaining nectar-rich food source while avoiding the exhaustion of the original food source nectar. If the standby bee finds a new source of food, compare the target value of the new food source with the original food source position. If the new food source has a better target value than the original position, replace the old food source with the new food source. When performing the exploration process, the investigation bee will be called to visit various places, and random search for possible food sources (new feasible solutions), planning the behavior of the investigation bee to explore new food sources, and using the detection bee as the basis for decision-making, eliminating the target function value Poor food sources and replaced with new food sources. Until the iterative generation reaches the upper limit of the iteration set by the program, the artificial bee colony algorithm calculus is stopped and the optimal solution is output. In the present invention, the optimal solution is the optimal weight coefficient combination of the individual prediction methods.

本發明應用人工蜂群演算法的群體尋優特性,求解個別預測方法的權重係數組合,其基本概念為將個別預測方法的權重係數組合設為人工蜂群演算法每個食物源,透過蜂群在搜尋空間的尋優過程找出最佳解,以確定個別預測方法各權重係數組合的最佳值(如圖2步驟S108、S110所示)。其執行步驟詳述於後。The invention applies the group optimization characteristic of the artificial bee colony algorithm to solve the weight coefficient combination of the individual prediction methods, and the basic concept is that the weight coefficient combination of the individual prediction methods is set as the artificial bee colony algorithm for each food source, through the bee colony. The optimal solution is found in the optimization process of the search space to determine the optimal value of each weight coefficient combination of the individual prediction methods (as shown in steps S108 and S110 of FIG. 2). The execution steps are detailed later.

步驟1:設定初始資料。在此步驟中,設定食物源數量、變數維度、設計變數上下限值、疊代次數上限、食物源限制條件。以(4)式產生N組變數組合,以(5)式計算變數之適應值,並紀錄最佳解。Step 1: Set the initial data. In this step, set the number of food sources, variable dimensions, upper and lower limits of design variables, upper limit of iterations, and food source restrictions. Generate N sets of variable combinations by (4) and calculate variables with (5) Adapt to the value and record the best solution.

(4) (4)

其中,為第j 個變數組合的第i 個變數i =1, 2, …,NN 為食物源數量,也稱為群體數,即為變數組合之數量;j =1, 2, …,DD 為維度;𝛼為一個介於0~1之間的隨機亂數;為變數組合之設計變數的上下限。among them, J-th combination for the i-th variable variable i = 1, 2, ..., N, N is the number of food sources, also referred to as the number of groups, i.e. the number of combinations of variables; j = 1, 2, ... , D, D is a dimension; 𝛼 is a random random number between 0 and 1; versus The upper and lower limits of the design variables for the variable combination.

上述適應值之計算如(5)式所示:The above adaptation values are calculated as shown in equation (5):

(5) (5)

其中,的成本值亦即個別預測方法絕對誤差之權重和。among them, for The cost value is also the sum of the weights of the absolute errors of the individual prediction methods.

步驟2:以鄰域搜尋產生新變數。在此步驟中,以(6)式產生新的設計變數,並限制不可超過上下限值。Step 2: Generate new variables by neighborhood search. In this step, generate new design variables with (6) And the limit cannot exceed the upper and lower limits.

(6) (6)

其中,為第j 個變數組合的第i 個新設計變數,為第j 個變數組合的第i 個舊設計變數,為一個隨機產生的數值介於-1~1之間,為群體中隨機選擇的一個變數組合,among them, The ith new design variable for the jth variable combination, The i-th old design variable for the j- th variable combination, For a randomly generated value between -1 and 1, a combination of variables randomly selected for the population, .

步驟3:決定是否更新變數組合。計算變數組合之適應值,若之適應值較佳則更新變數組合。Step 3: Decide whether to update the variable combination. Calculated variable combination Adaptation value, if Preferably, the adaptation value is updated to change the combination of variables.

步驟4:決定是否進入觀察蜂階段:Step 4: Decide whether to enter the observation bee stage:

變數組合經隨機產生一個0~1之間的亂數x ,並計算變數組合的篩選機率,若x 小於篩選機率則變數組合進入觀察蜂階段。Variable combination Randomly generate a random number x between 0 and 1, and calculate the combination of variables Screening probability, if x is less than the screening probability, the variable combination Enter the observation bee stage.

上述篩選機率之計算如(7)式所示:The above screening probability is calculated as shown in (7):

(7) (7)

其中,為變數組合的篩選機率,的適應值,N 為食物源數量。among them, Combination of variables Screening probability, for The fitness value, N is the number of food sources.

步驟5:觀察蜂階段。在此步驟中,以(6)式產生新的設計變數,並限制不可超過上下限值。計算變數組合之適應值,挑選適應值較好的變數組合。Step 5: Observe the bee stage. In this step, generate new design variables with (6) And the limit cannot exceed the upper and lower limits. Calculated variable combination For the fitness value, select a combination of variables with better fitness values.

步驟6:決定是否進入偵查蜂階段。在此步驟中,判斷是否有食物源達限制條件:如果達限制條件則放棄該食物源進入偵查蜂階段,否則跳過偵查蜂階段。Step 6: Decide whether to enter the scout bee stage. In this step, it is judged whether there is a food source reaching the restriction condition: if the restriction condition is reached, the food source is abandoned to enter the detection bee stage, otherwise the detection bee stage is skipped.

步驟7:偵查蜂階段。針對第j 組變數組合以(4)式產生新的變數組合,計算適應值。Step 7: Detect the bee stage. Generating a new combination of variables with the formula (4) for the j-th set of variables , calculate the fitness value.

步驟8:檢查是否符合結束條件。若搜尋疊代數達到設定之疊代次數上限值,則結束疊代並輸出最佳權重係數組合;否則回到步驟2繼續進行疊代。Step 8: Check if the end condition is met. If the search iteration number reaches the set top generation limit value, the iteration is ended and the optimal weight coefficient combination is output; otherwise, return to step 2 to continue the iteration.

2、持續法(persistence method)2, persistence method (persistence method)

持續法不僅原理簡單,也是最經濟的風速度或功率預測的方法,各國的電力公司也常以持續法作為超短期風力發電量預測的基準方法。持續法的原理是基於一個簡單的假設:當預測的時間區間極短時,在目前量測時間點與未來預測時間點的風速或風力發電量將是相同的。在某些特定狀況下,持續法用於超短期風力發電量預測時比其他風力發電量預測法更準確(如圖2步驟S102所示)。The continuous method is not only a simple principle, but also the most economical method of wind speed or power prediction. Power companies in various countries often use the continuous method as the benchmark method for predicting ultra-short-term wind power generation. The principle of the continuation method is based on a simple assumption: when the predicted time interval is extremely short, the wind speed or wind power generation at the current measurement time point and the future prediction time point will be the same. Under certain conditions, the continuous method is more accurate for the prediction of ultra-short-term wind power than other wind power generation prediction methods (as shown in step S102 of Fig. 2).

持續法的原理係基於一個簡單的假設:如果目前時間點測得的風速與風力發電量為,則未來的時間點()的預測風速與風力發電量為:The principle of the continuation method is based on a simple assumption: if the current time point The measured wind speed and wind power generation are versus , then the future time point ( The predicted wind speed and wind power generation are:

(8) (8)

(9) (9)

3、徑向基底函數(radial basis function, RBF)類神經網路演算法3. Radial basis function (RCF)-like neural network algorithm

類神經網路常被用來解決參數估算的問題,在各類類神經網路中RBF類神經網路具有強健性高與近似能力強的優點,因此特別適合應用於非線性估測,本發明的短期風力發電量預測系統也將採用RBF類神經網路為三種個別預測系統之一(如圖2步驟S104所示)。RBF類神經網路的基本原理來自function approximation,RBF類神經網路的基本架構如圖3所示,係為輸入層、隱藏層與輸出層所構成的三層網路架構。Neural networks are often used to solve the problem of parameter estimation. RBF-like neural networks have the advantages of high robustness and strong approximation in various neural networks. Therefore, they are particularly suitable for nonlinear estimation. The short-term wind power prediction system will also use the RBF-like neural network as one of three individual prediction systems (as shown in step S104 of FIG. 2). The basic principle of RBF neural network comes from function approximation. The basic architecture of RBF neural network is shown in Figure 3. It is a three-layer network architecture composed of input layer, hidden layer and output layer.

若一個RBF類神經網路的輸入層維度為d ,隱藏層維度為q ,輸出層維度為m ,則RBF類神經網路為d 維至m 維的映射關係,如下式所示:If the input layer dimension of an RBF-like neural network is d , the hidden layer dimension is q , and the output layer dimension is m , the RBF-like neural network is a d -dimensional to m -dimensional mapping relationship, as shown in the following equation:

(10) (10)

其中,輸入向量為={,fori = 1, 2, …, d},輸出向量為={,fori = 1, 2, …, m}。Where the input vector is ={ , for i = 1, 2, ..., d}, the output vector is ={ , for i = 1, 2, ..., m}.

RBF類神經網路的隱藏層各節點的徑向基底函數(radial basis function)通常係為高斯函數(Gaussian function) ,如下式所示:Radial basis function of each node of the hidden layer of the RBF-like neural network Usually it is a Gaussian function, as shown in the following equation:

, forj = 1, 2, …, q                             (11) , for j = 1, 2, ..., q (11)

其中,為高斯函數的寬,為高斯函數的中心。among them, For the width of the Gaussian function, Is the center of the Gaussian function.

RBF類神經網路的輸出層各節點的輸出值,如下式所示:Output value of each node of the output layer of the RBF-like neural network , as shown below:

, fork = 1, 2, …, m                                   (12) , for k = 1, 2, ..., m (12)

其中,為隱藏層第j 個節點至輸出層第k 個節點之間的權重,為隱藏層第j 個節點的輸出值。among them, To hide the weight between the jth node of the layer and the kth node of the output layer, The output value of the jth node of the hidden layer.

由RBF類神經網路架構可知,當建立訓練資料對後,則輸入、輸出神經元的數目即已確定。至於隱藏層神經元的數目,通常視為輸入資料的一子集合(subset),亦即以輸入資料的機率密度函數(probability density function)決定其神經元數目。當隱藏層神經元數目過大時,網路學習效果較佳,但訓練時間勢必增加,因而一合適的神經元數目將有助於提升網路學習效果,本發明中應用正交最小平方理論選出最佳的隱藏層節點數目。According to the RBF-like neural network architecture, when the training data pair is established, the number of input and output neurons is determined. As for the number of hidden layer neurons, it is usually regarded as a subset of the input data, that is, the number of neurons is determined by the probability density function of the input data. When the number of hidden layer neurons is too large, the network learning effect is better, but the training time is bound to increase, so a suitable number of neurons will help to improve the network learning effect. In the present invention, the orthogonal least squares theory is used to select the most. The number of good hidden layer nodes.

4、倒傳遞(back propagation, BP)類神經網路演算法4. Back propagation (BP) neural network algorithm

在各類類神經網路中BP類神經網路具有易收斂與對應映射性強的優點,因此特別適合應用於預測,本發明的短期風力發電量預測系統也將採用BP類神經網路為三種短期預測系統之一(如圖2步驟S106所示)。BP類神經網路的基本架構如圖4所示,係由輸入層、隱藏層與輸出層所構成的三層網路架構。In all kinds of neural networks, BP neural networks have the advantages of easy convergence and corresponding mapping, so they are particularly suitable for prediction. The short-term wind power prediction system of the present invention will also adopt BP neural networks for three types. One of the short-term prediction systems (as shown in step S106 of Fig. 2). The basic architecture of the BP-like neural network is shown in Figure 4. It is a three-layer network architecture consisting of an input layer, a hidden layer, and an output layer.

若一個BP類神經網路的輸入層維度為d ,隱藏層維度為q ,輸出層維度為m ,則BP類神經網路為d 維至m 維的映射關係,如(13)式所示:If the input layer dimension of a BP neural network is d , the hidden layer dimension is q , and the output layer dimension is m , the BP neural network is a d -dimensional to m -dimensional mapping relationship, as shown in (13):

(13) (13)

其中,輸入向量為={,fori = 1, 2, …, d},輸出向量為={,fori = 1, 2, …, m}。Where the input vector is ={ , for i = 1, 2, ..., d}, the output vector is ={ , for i = 1, 2, ..., m}.

BP類神經網路的隱藏層第j 個節點的輸入值,如(14)式所示:Input value of the jth node of the hidden layer of the BP-like neural network As shown in (14):

(14) (14)

其中,為輸入層第i 個節點至隱藏層第j 個節點之間的權重,為隱藏層第j 個節點的偏權值。among them, The i-th input layer nodes to the right between the j-th hidden layer node weight, To hide the partial weight of the jth node of the layer.

BP類神經網路的隱藏層第j 個節點的輸出值,如(15)式所示:The output value of the jth node of the hidden layer of the BP-like neural network , as shown in (15):

(15) (15)

其中,為隱藏層第j 個節點的輸入值。among them, To hide the input value of the jth node of the layer.

BP類神經網路的輸出層第j 個節點的輸入值,如(16)式所示:Input value of the jth node of the output layer of the BP-like neural network As shown in (16):

(16) (16)

其中,為隱藏層第i 個節點至輸出層第j 個節點之間的權重,為輸出層第j 個節點的偏權值。among them, To hide the weight between the i- th node of the layer and the j- th node of the output layer, The offset value of the jth node of the output layer.

BP類神經網路的輸出層第j 個節點的輸出值,如(17)式所示:The output value of the jth node of the output layer of the BP-like neural network As shown in (17):

(17) (17)

其中,為輸出層第j 個節點的輸入值。among them, The input value of the jth node of the output layer.

由BP類神經網路架構可知,當建立訓練資料對後,則輸入、輸出神經元的數目即已確定。至於隱藏層神經元的數目,通常視為輸入資料的一子集合(subset),亦即以輸入資料的機率密度函數(probability density function)決定其神經元數目。當隱藏層神經元數目過大時,網路學習效果較佳,但訓練時間勢必增加,因而一合適的神經元數目將有助於提升網路學習效果,本發明中應用正交最小平方理論選出最佳的隱藏層節點數目。According to the BP-like neural network architecture, when the training data pair is established, the number of input and output neurons is determined. As for the number of hidden layer neurons, it is usually regarded as a subset of the input data, that is, the number of neurons is determined by the probability density function of the input data. When the number of hidden layer neurons is too large, the network learning effect is better, but the training time is bound to increase, so a suitable number of neurons will help to improve the network learning effect. In the present invention, the orthogonal least squares theory is used to select the most. The number of good hidden layer nodes.

本發明採用之BP類神經網路架構如圖4所示,共分為三層,其中第一層輸入層共有4個節點,分別風力發電機組的目前實際發電量、前10分鐘的實際發電量、前20分鐘的實際發電量與風速預測值等4種輸入信號;第二層為隱藏層,節點數將依照正交最小平方理論選出之最佳的隱藏層節點數目而調整;第三層為輸出層,只有1個節點即為10分鐘後風力發電量的預測值。The BP-like neural network architecture adopted by the present invention is divided into three layers as shown in FIG. 4, wherein the first input layer has four nodes, respectively, the actual actual power generation of the wind turbine and the actual power generation in the first 10 minutes. 4 kinds of input signals such as the actual power generation and wind speed prediction in the first 20 minutes; the second layer is the hidden layer, and the number of nodes will be adjusted according to the optimal number of hidden layer nodes selected by the orthogonal least squares theory; the third layer is In the output layer, only one node is the predicted value of wind power generation after 10 minutes.

此外,本發明中所採用之正交最小平方法(orthogonal least squares, OLS)決定類神經網路的最佳隱藏層節點數目說明加下。正交最小平方法主要採用Gram-Schmidt正交理論將任意的向量矩陣分解成一組正交基底向量,並使其涵蓋的空間與原始向量矩陣涵蓋的空間相同,因此若將此理論應用於選取RBF類神經網路與BP類神經網路之隱藏層節點數目時,則足以涵蓋原始訓練資料所涵蓋的空間。假設Y 為訓練資料之期望輸出向量,則:In addition, the orthogonal least squares (OLS) method used in the present invention determines the number of optimal hidden layer nodes of the neural network to be added. The orthogonal least squares method mainly uses Gram-Schmidt orthogonal theory to decompose an arbitrary vector matrix into a set of orthogonal basis vectors, and makes the space covered by the same as the space covered by the original vector matrix, so if this theory is applied to the selection of RBF The number of hidden layer nodes of the neural network and the BP neural network is sufficient to cover the space covered by the original training data. Assuming Y is the expected output vector of the training data, then:

(18) (18)

或寫為矩陣型式Or written as a matrix type

(19) (19)

其中,Y ÎÂ m ´ 1B ÎÂ m ´ q E ÎÂ m ´ 1 ,WÎÂ q ´ 1m 為輸出層神經元數目。Wherein, 1, m is the number of neurons in the output layer, Y Î Â m '1, B Î Â m' q, E Î Â m '1, WÎ Â q'.

依據Gram-Schmidt正交理論,B 矩陣分解成一組正交基底向量如下式:According to the Gram-Schmidt orthogonal theory, the B matrix is decomposed into a set of orthogonal basis vectors as follows:

(20) (20)

其中,A ÎÂ q ´ q 為一上三角矩陣,D ÎÂ 1 ´ q 為一正交基底向量矩陣。且Wherein, A Î Â q 'q is an upper triangular matrix, D Î Â 1' q matrix is an orthogonal basis vectors. And

(21) (twenty one)

其中,為正對角矩陣(positive diagonal matrix)H 的元素。整合(7)至(9)式可得:among them, An element of the positive diagonal matrix H. Integration (7) to (9) can be obtained:

(22) (twenty two)

其中,G =AWWhere G = AW .

由於D 為一正交基底向量矩陣,所以期望輸出值Y 的平方和可表示為:Since D is an orthogonal basis vector matrix, the sum of squares of the expected output values Y can be expressed as:

(23) (twenty three)

因此,針對第k 個節點時,其誤差下降率為:Therefore, for the k-th node, the error rate is lowered:

(24) (twenty four)

而節點的篩選流程則是從訓練資料中逐漸挑出誤差下降率最大值的點當作新加入的節點,並將累加,直到滿足所設定的正確率為止。另外,在選取中心節點個數的過程中,本發明將採用伸展度相同的高斯函數,亦即其標準偏差為固定值,如下式所示:The node's screening process is to gradually pick out the maximum error reduction rate from the training data. Point as a newly joined node and Accumulate until the set correct rate is met. In addition, in the process of selecting the number of central nodes, the present invention will adopt a Gaussian function with the same extension, that is, the standard deviation is a fixed value, as shown in the following formula:

(25) (25)

其中,為所有節點中心間最大的距離值,q 為節點數目。among them, The maximum distance value between all node centers, q is the number of nodes.

此外,為進行智慧型風力發電量預測系統的訓練與測試程序,本發明將採用風力發電機組的一年實際發電資料進行測試,風力發電機組的發電資料是每十分鐘量測一次,因此,一年的實際發電資料共有52,560筆資料,本發明將前述正規化之發電資料建立資料庫。由於台灣西部海岸一年四季的風力變化不同,本發明將分別針對四季設計風力發電量預測系統,每季的實際發電資料皆建成一個資料庫,每個資料庫有13,140個資料檔,由各資料庫選出每季前兩個月的資料檔建檔為訓練資料庫,每季後一個月的資料檔則建檔為測試資料庫,這些參數資料庫將提供智慧型風力發電量預測系統的訓練與測試等使用。In addition, in order to carry out the training and testing procedures of the intelligent wind power generation quantity prediction system, the present invention will use the one-year actual power generation data of the wind power generation unit to be tested, and the power generation data of the wind power generation unit is measured every ten minutes, therefore, one The actual power generation data of the year has a total of 52,560 data. The present invention establishes a database of the aforementioned conventional power generation data. Due to the different wind changes in the western coast of Taiwan throughout the year, the present invention will separately design a wind power generation forecasting system for the four seasons. Each quarter of the actual power generation data is built into a database, each database has 13,140 data files, each data The database selects the data files for the first two months of each season as the training database. The data files for each month after each season are filed as test data bases. These parameter databases will provide training and training for intelligent wind power generation forecasting systems. Test and other use.

綜上所述,本發明係具有以下之特色與優點:In summary, the present invention has the following features and advantages:

1、本發明之研究目的為針對風力發電量預測提出一套以人工蜂群演算法結合持續法、徑向基底函數(RBF)類神經網路演算法、倒傳遞(BP)類神經網路演算法的智慧型短期風力發電量預測方法,本方法的特色是具有極大彈性、能適用於各種類型的風力發電機組發電量預測;1. The purpose of the research of the present invention is to propose a set of artificial bee colony algorithm combined with continuous method, radial basis function (RBF) neural network algorithm and reverse transfer (BP) neural network algorithm for wind power generation prediction. Intelligent short-term wind power generation forecasting method, the method is characterized by great flexibility and can be applied to various types of wind turbine generating power generation prediction;

2、以複合方程式計算出發電量預測量預測值,且複合方程式每種個別預測方法的權重係數將以人工蜂群演算法求解,得到最佳化智慧型預測方程式;不僅避免求解過程陷入局部解,並可提升求解最佳化問題之運算效能,可以有效提高智慧型預測法的預測準確性;2. Calculate the predicted value of the starting electricity quantity by the compound equation, and the weighting coefficient of each individual prediction method of the compound equation will be solved by the artificial bee colony algorithm, and the optimized intelligent prediction equation is obtained; not only the solution process is prevented from falling into the local solution. The computational efficiency of solving the optimization problem can be improved, and the prediction accuracy of the intelligent prediction method can be effectively improved;

3、該智慧型預測法以PB類神經網路與RBF類神經網路為基礎,類神經網路常用於處理非線性且複雜的分類或預測問題,通過訓練過程類神經網路可以擷取出變量間的非線性,且預測的準確性高於傳統的統計法;3. The intelligent prediction method is based on the PB neural network and the RBF-like neural network. The neural network is often used to deal with nonlinear and complex classification or prediction problems. The training-like neural network can extract variables. Non-linearity, and the accuracy of prediction is higher than traditional statistical methods;

4、該智慧型預測法中RBF類神經網路與PB類神經網路有關隱藏層節點的數目為決定網路預測效果的重要參數,本發明將應用正交最小平方(Orthogonal Least-Squares)理論選出最佳的隱藏層節點數目,可以提高RBF類神經網路與PB類神經網路的預測準確性; 5、本發明採用風力發電機組的實際歷史發電資料進行RBF類神經網路與PB類神經網路訓練,各種類型的風力發電機組只需將實際歷史發電資料建檔,即可訓練類神經網路,建構專屬該機組的短期風力發電量預測系統,因此,本發明所提出的智慧型風力發電量預測法具有極大的彈性、能適用於各種類型的風力發電機組短期發電量預測;及 6、短期風力發電系統發電量預測系統能有效與正確的預測風力發電系統在不同季節狀態下的發電量,提供電力系統調度人員進行相關電力系統運轉控制時的重要工具,以達到提高電力系統穩定度、降低系統運轉成本目的。4. In the intelligent prediction method, the number of hidden layer nodes related to the RBF neural network and the PB neural network is an important parameter for determining the network prediction effect, and the present invention applies the orthogonal least square (Orthogonal Least-Squares) theory. Selecting the optimal number of hidden layer nodes can improve the prediction accuracy of RBF-like neural networks and PB-like neural networks. 5. The present invention uses the actual historical power generation data of wind turbines to perform RBF-like neural networks and PB-like nerves. Network training, all types of wind turbines only need to archive the actual historical power generation data, and then train the neural network to construct a short-term wind power prediction system dedicated to the unit. Therefore, the intelligent wind power proposed by the present invention The power generation forecasting method has great flexibility and can be applied to the short-term power generation forecast of various types of wind turbines; and 6. The short-term wind power generation system power generation forecasting system can effectively and correctly predict the wind power generation system in different seasons. Quantitatively, providing power system dispatchers with important tools for controlling the operation of relevant power systems to achieve increased power The stability of the force system and the purpose of reducing the operating cost of the system.

惟,以上所述,僅為本發明較佳具體實施例之詳細說明與圖式,惟本發明之特徵並不侷限於此,並非用以限制本發明,本發明之所有範圍應以下述之申請專利範圍為準,凡合於本發明申請專利範圍之精神與其類似變化之實施例,皆應包含於本發明之範疇中,任何熟悉該項技藝者在本發明之領域內,可輕易思及之變化或修飾皆可涵蓋在以下本案之專利範圍。However, the above description is only for the detailed description and the drawings of the preferred embodiments of the present invention, and the present invention is not limited thereto, and is not intended to limit the present invention. The scope of the patent application is intended to be included in the scope of the present invention, and any one skilled in the art can readily appreciate it in the field of the present invention. Variations or modifications may be covered by the patents in this case below.

﹝本發明﹞﹝this invention﹞

S102~S110‧‧‧步驟S102~S110‧‧‧Steps

x1~xd‧‧‧徑向基底函數類神經網路演算法之輸入層節點Input layer node of x1~xd‧‧‧radial basis function neural network algorithm

b1~bq‧‧‧徑向基底函數類神經網路演算法之隱藏層節點B1~bq‧‧‧ hidden layer nodes of radial basis function neural network algorithm

y1~ym‧‧‧徑向基底函數類神經網路演算法之輸出層節點Y1~ym‧‧‧Output layer nodes of radial basis function neural network algorithm

x1~xd‧‧‧倒傳遞類神經網路演算法之輸入層節點Input layer node of x1~xd‧‧‧ inverse transfer neural network algorithm

b1~bq‧‧‧倒傳遞類神經網路演算法之隱藏層節點B1~bq‧‧‧ hidden layer nodes of inverse transfer neural network algorithm

第一圖係為本發明智慧型短期風力發電系統發電量預測系統之系統架構圖;The first figure is a system architecture diagram of the power generation quantity prediction system of the smart short-term wind power generation system of the present invention;

第二圖係為本發明智慧型短期電力發電量預測方法之流程圖;The second figure is a flow chart of the method for predicting the smart short-term power generation amount of the present invention;

第三圖係為本發明徑向基底函數類神經網路演算法之架構圖;及The third figure is an architectural diagram of the radial basis function neural network algorithm of the present invention; and

第四圖係為本發明倒傳遞類神經網路演算法之架構圖。The fourth figure is an architectural diagram of the inverse transfer neural network algorithm of the present invention.

S102~S110‧‧‧步驟 S102~S110‧‧‧Steps

Claims (10)

一種智慧型短期電力發電量預測方法,係包含: (a) 提供一持續法,根據目前時間點測得的風速與風力發電量,預測未來的時間點的預測風速與風力發電量:;其中,為目前時間點、()為未來時間點、為目前風速、為未來風速、為目前風力發電量以及為未來風力發電量; (b) 提供一徑向基底函數類神經網路演算法,利用包括一輸入層、一隱藏層以及一輸出層之網路架構,預測未來風力發電量; (c) 提供一倒傳遞類神經網路演算法,利用包括一輸入層、一隱藏層以及一輸出層之網路架構,預測未來風力發電量; (d) 根據該持續法、該徑向基底函數類神經網路演算法以及該倒傳遞類神經網路演算法,提供一智慧型數學預測模型;其中,;以及(t = 1, 2, …,L )為實際的發電量時間序列數據,M為個別預測方法數量,L為樣本數,(i = 1, 2, …,M ,t = 1, 2, …,L )係為第i種預測方法的預測值,=-為預測誤差,為第i種預測方法的權重係數,的估測值,而為智慧型預測法的預測值;及 (e) 提供一人工蜂群演算法,求解該智慧型數學預測模型之權重係數,以獲得風力發電量最佳值;其中,A smart short-term power generation forecasting method includes: (a) providing a continuous method for predicting predicted wind speed and wind power generation at a future time point based on wind speed and wind power generation measured at current time points: , ;among them, For the current time, ( ) for the future time, For the current wind speed, For future wind speed, For current wind power generation as well (b) provide a radial basis function-like neural network algorithm that predicts future wind power generation using a network architecture that includes an input layer, a hidden layer, and an output layer; (c) provides one An inverse transfer-like neural network algorithm predicts future wind power generation using a network architecture including an input layer, a hidden layer, and an output layer; (d) according to the persistence method, the radial basis function-like neural network algorithm And the inverse transfer neural network algorithm provides a smart mathematical prediction model ;among them, , ;as well as ( t = 1, 2, ..., L ) is the actual power generation time series data, M is the number of individual prediction methods, and L is the number of samples. ( i = 1, 2, ..., M , t = 1, 2, ..., L ) is the predicted value of the i-th prediction method, = - For prediction errors, The weighting factor for the i-th prediction method, for Estimated value, and a predictive value for the intelligent predictive method; and (e) providing an artificial bee colony algorithm to solve the weighting coefficient of the intelligent mathematical predictive model To obtain the best value of wind power generation; . 如申請專利範圍第1項該智慧型短期電力發電量預測方法,其中在步驟(b)中,該徑向基底函數類神經網路演算法(RBF Neural Network)為:, 其中,d為該輸入層維度、q為該隱藏層維度、m為該輸出層維度;以及輸入向量為={,fori = 1, 2, …, d},輸出向量為={,fori = 1, 2, …, m}。For example, the intelligent short-term power generation prediction method according to the first application of the patent scope, wherein in step (b), the RBF Neural Network is: Where d is the input layer dimension, q is the hidden layer dimension, m is the output layer dimension; and the input vector is ={ , for i = 1, 2, ..., d}, the output vector is ={ , for i = 1, 2, ..., m}. 如申請專利範圍第2項該智慧型短期電力發電量預測方法,其中該隱藏層各節點之徑向基底函數為,係為一高斯函數:, forj = 1, 2, …, q, 其中,為高斯函數的寬,為高斯函數的中心; 其中,該隱藏層之節點數目係透過一正交最小平方理論決定。For example, the smart short-term power generation prediction method according to item 2 of the patent application scope, wherein the radial basis function of each node of the hidden layer is , is a Gaussian function: , for j = 1, 2, ..., q, where, For the width of the Gaussian function, It is the center of the Gaussian function; wherein the number of nodes of the hidden layer is determined by an orthogonal least squares theory. 如申請專利範圍第2項該智慧型短期電力發電量預測方法,其中該輸出層之各節點輸出值為, fork = 1, 2, …, m, 其中,為隱藏層第j 個節點至輸出層第k 個節點之間的權重,為隱藏層第j 個節點的輸出值。For example, the smart short-term power generation prediction method according to item 2 of the patent application scope, wherein the output values of the nodes of the output layer are : , for k = 1, 2, ..., m, where, To hide the weight between the jth node of the layer and the kth node of the output layer, The output value of the jth node of the hidden layer. 如申請專利範圍第1項該智慧型短期電力發電量預測方法,其中在步驟(c)中,該倒傳遞類神經網路演算法(BP Neural Network)為:, 其中,d為該輸入層維度、q為該隱藏層維度、m為該輸出層維度;以及輸入向量為={,fori = 1, 2, …, d},輸出向量為={,fori = 1, 2, …, m}。For example, the smart short-term power generation prediction method according to the first application of the patent scope, wherein in step (c), the BP Neural Network algorithm is: Where d is the input layer dimension, q is the hidden layer dimension, m is the output layer dimension; and the input vector is ={ , for i = 1, 2, ..., d}, the output vector is ={ , for i = 1, 2, ..., m}. 如申請專利範圍第5項該智慧型短期電力發電量預測方法,其中該隱藏層各節點之輸入值為:, 其中,為輸入層第i 個節點至隱藏層第j 個節點之間的權重,為隱藏層第j 個節點的偏權值; 其中,該隱藏層第j 個節點的輸出值為:, 其中,係為該隱藏層第j 個節點之輸入值; 其中,該隱藏層之節點數目係透過一正交最小平方理論決定。For example, the intelligent short-term power generation amount prediction method in the fifth application patent scope, wherein the input values of the nodes of the hidden layer are for: , among them, The i-th input layer nodes to the right between the j-th hidden layer node weight, Partial weight of the j th hidden layer node; wherein an output value of the j-th hidden layer node for: , among them, The input value of the jth node of the hidden layer; wherein the number of nodes of the hidden layer is determined by an orthogonal least squares theory. 如申請專利範圍第5項該智慧型短期電力發電量預測方法,其中該輸出層第j 個節點的輸入值為:, 其中,係該隱藏層第i 個節點至輸出層第j 個節點之間的權重,係該輸出層第j 個節點的偏權值; 其中,該輸出層第j 個節點的輸出值為:, 其中,為輸出層第j 個節點的輸入值。For example, the intelligent short-term power generation amount prediction method in the fifth application patent scope, wherein the input value of the j- th node of the output layer for: , among them, The system hidden layer node to the i-th output layer weights between the weight of the j-th node, Partial weight based on the j-th output layer node; wherein the output value of the output layer nodes j for: , among them, The input value of the jth node of the output layer. 如申請專利範圍第1項該智慧型短期電力發電量預測方法,其中在步驟(e)中,更包含: (e1) 設定初始資料,;其中,為第j 個變數組合的第i 個變數i =1, 2, …,NN 為食物源數量;j =1, 2, …,DD 為維度;𝛼為一個介於0~1之間的隨機亂數;為變數組合之設計變數的上下限; (e2) 以鄰域搜尋產生新變數,;其中,為第j 個變數組合的第i 個新設計變數,為第j 個變數組合的第i 個舊設計變數,為一個隨機產生的數值介於-1~1之間,為群體中隨機選擇的一個變數組合,; (e3) 決定是否更新變數組合; (e4) 決定是否進入觀察蜂階段,;其中,為變數組合的篩選機率,的適應值,N 為食物源數量; (e5) 觀察蜂階段; (e6) 決定是否進入偵查蜂階段; (e7) 偵查蜂階段; (e8) 檢查是否符合結束條件,若搜尋疊代數達到設定之疊代次數上限值,則結束疊代並輸出最佳權重係數組合,否則回到步驟(e2)繼續進行疊代。For example, the method for predicting the smart short-term power generation amount in the first application of the patent scope includes, in step (e), further including: (e1) setting initial data, ;among them, The i- th variable i = 1, 2, ..., N , N for the j- th variable is the number of food sources; j = 1, 2, ..., D , D are dimensions; 𝛼 is a Random random number between ~1; versus The upper and lower limits of the design variables for the combination of variables; (e2) the generation of new variables by neighborhood search, ;among them, The ith new design variable for the jth variable combination, The i-th old design variable for the j- th variable combination, For a randomly generated value between -1 and 1, a combination of variables randomly selected for the population, (e3) decide whether to update the variable combination; (e4) decide whether to enter the observation bee stage, ;among them, Combination of variables Screening probability, for Adaptation value, N is the number of food sources; (e5) observation bee stage; (e6) decide whether to enter the detection bee stage; (e7) detection bee stage; (e8) check whether the end condition is met, if the search iteration reaches the set If the upper limit of the iteration is exceeded, the iteration is ended and the optimal weight coefficient combination is output, otherwise it returns to step (e2) to continue the iteration. 如申請專利範圍第1項該智慧型短期電力發電量預測方法,其中在步驟(a)之前,更包含: (a01) 風力發電量資料庫建檔; (a02) 風力發電量資料庫正規化;及 (a03) 訓練資料庫建檔。For example, the smart short-term power generation forecasting method in the first application scope of the patent scope, wherein before step (a), further comprises: (a01) wind power generation database archive; (a02) wind power generation database normalization; And (a03) training database file. 如申請專利範圍第1項該智慧型短期電力發電量預測方法,其中在步驟(e)之後,更包含: (f) 風力發電量預測反正規化。For example, the intelligent short-term power generation quantity prediction method in the first application scope of the patent scope, wherein after step (e), further comprises: (f) anti-normalization of wind power generation prediction.
TW104108689A 2015-03-18 2015-03-18 Method of short-term wind power generation forecasting TWI540533B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW104108689A TWI540533B (en) 2015-03-18 2015-03-18 Method of short-term wind power generation forecasting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW104108689A TWI540533B (en) 2015-03-18 2015-03-18 Method of short-term wind power generation forecasting

Publications (2)

Publication Number Publication Date
TWI540533B TWI540533B (en) 2016-07-01
TW201635224A true TW201635224A (en) 2016-10-01

Family

ID=56997027

Family Applications (1)

Application Number Title Priority Date Filing Date
TW104108689A TWI540533B (en) 2015-03-18 2015-03-18 Method of short-term wind power generation forecasting

Country Status (1)

Country Link
TW (1) TWI540533B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI623890B (en) * 2017-07-04 2018-05-11 台灣電力股份有限公司 System for predicating power generation by utilizing multiple neural networks and method thereof
TWI662423B (en) * 2017-02-06 2019-06-11 台灣電力股份有限公司 Display system and method for wind power prediction
CN110009037A (en) * 2019-04-03 2019-07-12 中南大学 A kind of engineering wind speed Forecasting Approach for Short-term and system based on physical message coupling
TWI678677B (en) * 2017-05-25 2019-12-01 日商日立製作所股份有限公司 Adaptive power generation management
TWI702558B (en) * 2017-02-10 2020-08-21 日商東芝股份有限公司 Operation plan preparation device, operation plan preparation method and storage medium
US11196380B2 (en) 2017-07-04 2021-12-07 Taiwan Power Company Power generation prediction system and method thereof

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108259498B (en) * 2018-01-24 2020-06-23 湖南科技学院 Intrusion detection method and system based on BP algorithm of artificial bee colony optimization

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI662423B (en) * 2017-02-06 2019-06-11 台灣電力股份有限公司 Display system and method for wind power prediction
TWI702558B (en) * 2017-02-10 2020-08-21 日商東芝股份有限公司 Operation plan preparation device, operation plan preparation method and storage medium
TWI678677B (en) * 2017-05-25 2019-12-01 日商日立製作所股份有限公司 Adaptive power generation management
TWI623890B (en) * 2017-07-04 2018-05-11 台灣電力股份有限公司 System for predicating power generation by utilizing multiple neural networks and method thereof
US11196380B2 (en) 2017-07-04 2021-12-07 Taiwan Power Company Power generation prediction system and method thereof
CN110009037A (en) * 2019-04-03 2019-07-12 中南大学 A kind of engineering wind speed Forecasting Approach for Short-term and system based on physical message coupling
CN110009037B (en) * 2019-04-03 2020-10-27 中南大学 Short-term engineering wind speed prediction method and system based on physical information coupling

Also Published As

Publication number Publication date
TWI540533B (en) 2016-07-01

Similar Documents

Publication Publication Date Title
TWI540533B (en) Method of short-term wind power generation forecasting
US11581740B2 (en) Method, system and storage medium for load dispatch optimization for residential microgrid
Dong et al. Wind power day-ahead prediction with cluster analysis of NWP
CN106295899B (en) Wind power probability density Forecasting Methodology based on genetic algorithm Yu supporting vector quantile estimate
CN103065202B (en) Wind power plant ultrashort term wind speed prediction method based on combination kernel function
CN114742294A (en) Neural network algorithm for carbon emission prediction
CN110837915B (en) Low-voltage load point prediction and probability prediction method for power system based on hybrid integrated deep learning
CN111612244B (en) QRA-LSTM-based method for predicting nonparametric probability of photovoltaic power before day
CN104463356A (en) Photovoltaic power generation power prediction method based on multi-dimension information artificial neural network algorithm
CN110866633A (en) Micro-grid ultra-short term load prediction method based on SVR support vector regression
Xu et al. Intelligent forecasting model for regional power grid with distributed generation
CN108694475B (en) Short-time-scale photovoltaic cell power generation capacity prediction method based on hybrid model
CN116544934B (en) Power scheduling method and system based on power load prediction
CN108009684A (en) A kind of micro-grid connection state energy management method comprising short-term load forecasting
Yong et al. Short-term building load forecasting based on similar day selection and LSTM network
CN114611757A (en) Electric power system short-term load prediction method based on genetic algorithm and improved depth residual error network
CN114155044A (en) Power price prediction method and system for power spot market node
CN110991083B (en) Photovoltaic power station model determining method, device, equipment and storage medium
Srinivasan et al. Correlation analysis of solar power and electric demand
CN113283638A (en) Load extreme curve prediction method and system based on fusion model
Chen et al. Short-term wind power forecasting based on spatial correlation and artificial neural network
Chemetova et al. Load peak forecasting in different load patterns situations
Liu et al. Ultra-short-term wind power forecasting based on stacking model
He et al. Fault Causes Identification of Transmission Lines Based on Weighted Naive Bayes Classification Algorithm Combined with Complex Algorithm
CN112308327B (en) Smart city power load estimation method based on self-adaptive feature weight

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees