JP2023075072A - ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステム - Google Patents

ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステム Download PDF

Info

Publication number
JP2023075072A
JP2023075072A JP2022185214A JP2022185214A JP2023075072A JP 2023075072 A JP2023075072 A JP 2023075072A JP 2022185214 A JP2022185214 A JP 2022185214A JP 2022185214 A JP2022185214 A JP 2022185214A JP 2023075072 A JP2023075072 A JP 2023075072A
Authority
JP
Japan
Prior art keywords
wind speed
time series
support vector
rolling
garch
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.)
Granted
Application number
JP2022185214A
Other languages
English (en)
Other versions
JP7441553B2 (ja
Inventor
▲ジュ▼孝国
Xiaoguo Chu
楊政厚
Zhenghou Yang
韓健
Jian Han
周峰
Feng Zhou
岳紅軒
Hongxuan Yue
陳志文
Zhiwen Chen
陳兆聖
Zhaosheng Chen
陳卓
Zhuo Chen
王爽
Shuang Wang
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.)
Beijing Huaneng Xinrui Control Technology Co Ltd
Original Assignee
Beijing Huaneng Xinrui Control Technology 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 Beijing Huaneng Xinrui Control Technology Co Ltd filed Critical Beijing Huaneng Xinrui Control Technology Co Ltd
Publication of JP2023075072A publication Critical patent/JP2023075072A/ja
Application granted granted Critical
Publication of JP7441553B2 publication Critical patent/JP7441553B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/044—Recurrent networks, e.g. Hopfield networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • General Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Biophysics (AREA)
  • Development Economics (AREA)
  • Biomedical Technology (AREA)
  • Game Theory and Decision Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

【課題】風速予測の精度と信頼性を高めるローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステムを提供する。【解決手段】方法は、風速時系列を入力するステップと、位相空間を再構成するステップと、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するステップと、最小のサポートベクトルマシンを適用するステップと、風速の予測値を出力するステップと、を含む。【選択図】図１

Description

本出願は、風力発電の技術分野に属し、特にローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステムに関する。
近年、風力発電は電力システムに大規模に接続されるようになっている。風力発電は風の運動エネルギーを電気エネルギーに変換するが、自然界の風は間欠性と波動性の特徴を有するため、風力発電にも不安定性と波動性が存在する。風力発電の急速な発展に伴い、風力発電が電力ネットにおける発電に占める割合は増加しているが、これも電力ネットの電圧や周波数などの多くの問題と挑戦をもたらしている。風力発電機は、風力発電所の超短期風速予測を正確に行うことで、突然の切り出しによる電圧と周波数の変働を低減することができ、風力発電所のネットワーク化操作と風力発電機の制御に重要な役割を果たし、風力発電機制御の観点から風力発電ネットワーク化の能力を高めている。
従来の風速予測方法には、持続予測法、時系列法、カルマンフィルタ法、人工ニューラルネットワーク法、空間相関法などがある。持続予測法では、最近の何点かのデータを次点の風速予測値として重み付けるが、このような方法では予測誤差が大きく、他の予測方法の精度を評価する基準比較方法として用いられることが多い。ＡＲＭＡモデルとも称される時系列法では、過去のデータを用いて、モデル識別、パラメータ推定及びモデル検証によってモデル化し、このモデルを用いて予測モデルを導き出して、研究した時系列を予報する目的を達成する。時系列法では、必要とするモデル化情報が少なく、演算が簡単で便利、迅速であるが、下位モデルでは予測精度が低い場合が多く、上位モデルではパラメータ推定が困難であるため、時系列法は多くの改善が必要である。カルマンフィルタ法では、線形最小分散推定を用いて、風速の動的システムの状態系列に対して、必要な系列を抽出する方法であり、オンライン風速予測に適したこの方法は、騒音の統計的特性の推定に問題がある。人工ニューラルネットワーク法は、人の脳ニューラルネットワーク構造とその機能を模倣し、大量の簡単な処理ノードの大規模な相互接続から構成され、複雑な非線形問題を解決するために用いられ、その不足点は、過学習と汎化能力が強くないことに表れている。空間相関法では、風力発電所及び周辺の複数地点における複数組の実測の風速データを考慮し、各サイト間の風速の相関性を分析する必要があり、この方法の研究は絶えずに改善されている。
中国特許出願公開号ＣＮ１０５１８４３９Ａには、従来技術である、ウェーブレット分解とサポートベクトルマシンに基づく、風力発電所の風速と電力予測方法が開示されている。ウェーブレットパケット分解技術を利用して、過去の風速時系列をウェーブレットパケット分解し、過去の風速時系列の低周波数帯、中周波数帯及び高周波数帯の成分を得る。灰色サポートベクトルマシンの予測モデルを用いて、過去の風速時系列の各成分を予測し、その後、ウェーブレットパケットの再構成を用いて短期風速の予測データを得る。過去の風力発電電力データと数値天気予報の風速データを訓練セットとして利用し、灰色のサポートベクトルマシンモデルを構築し、風力発電電力の一次予測を行う。得られた風速の予測データ、風力発電電力の予測データに対して、ＲＢＦニューラルネットワークを通じて予測を行い、風力発電電力の最終予測値を得る。中国発明特許出願公開号ＣＮ１０７４４５０７Ａには、先行技術である、サポートベクトルマシンの組合せに基づく、風力発電所における短期風速の予測方法が開示されている。元の風速系列のカオス特性を分析し、Ｃ－Ｃアルゴリズムを用いて、風速系列の埋め込み次元ｍと遅延時間τを計算し、入出力データセットを構築し、訓練セット及び検証セットに分ける。訓練セットの上に、異なる核関数に基づくサポートベクトルマシンの単一予測モデルを構築し、総合学習戦略粒子群の最適化アルゴリズムを用いてキーパラメータを決定する。検証セットで、誘導秩序重み付け調和平均演算子に基づく変重み係数の組合せ予測法を用いて、各単一モデルの重み係数を確定し、それぞれ予測を行い、予測結果に重み付け加算を行うことによって、１歩の風速予測結果を得る。中国発明特許出願公開号ＣＮ１１３４４９８４７４７Ａには、先行技術である、秒単位の時系列の風速変化を考慮した、洋上風力発電電力のローリング予測方法が開示されている。データの前処理を利用して、異常データを処理し、洋上風速と風力発電電力の正規化を完成し、さらに差分平滑パワー系列を通じて、秒単位の時間スケールでの風速予測モデルを構築し、最後にローリングＬＳＴＭ記憶ネットワークを構築し、秒単位の時間系列データの予測を実現する。
しかし、上記の従来技術では風速予測の精度と安定性が低く、ＣＮ１０５１８４３９１Ａでは、風速系列のカオス特性を解釈することができず、ＣＮ１０７４４５０７Ａでは、風速系列のローリング予測を行っておらず、精度の向上に不利であり、ＣＮ１１３４４９８４７４７Ａでは、風速系列のカオス特性に対して分析を展開していない。
風速系列のカオス特性に対して、従来技術の不足を克服し、より正確で安定した風速予測方法を提供することが、当技術分野の早急な解決すべき課題となっている。
上記の従来技術の欠点を克服するために、本発明は、具体的に以下の技術的解決策を採用した、ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステムを提供する。
ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法は、
風速時系列を入力するステップＳ１、
位相空間を再構成するステップＳ２、
ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するステップＳ３、
最小のサポートベクトルマシンを適用するステップＳ４、及び
風速の予測値を出力するステップＳ５を含む。
また、
風速時系列を入力する前記ステップＳ１は、
元の風速時系列Ｐｂを入力するステップを含み、
位相空間を再構成する前記ステップＳ２は、
前記元の風速時系列Ｐｂに対して、以下の式（１）に従って位相空間を再構成し、空間位相を再構成する中間値系列を得るステップを含み、

Figure 2023075072000002
得られた前記空間位相を再構成する中間値系列に対して前処理を行うことによって、データの処理を容易にする。
Figure 2023075072000003
また、時系列モデルＡＲＭＡ（ｐ、ｑ）は、定常系列を予測するために用いられる以下の式（２）であり、

Figure 2023075072000004
前記時系列モデルＡＲＭＡ（ｐ、ｑ）は、モデル残差の平均値が０であり、残差の分散が定数であるという仮定を前提とする。
また、ＧＡＲＣＨモデルは以下の式（３）であり、
Figure 2023075072000005
また、ＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルは以下の式（４）であり、
Figure 2023075072000006
前記時系列モデルＡＲＭＡ（ｐ、ｑ）と前記ＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルを連立させ、以下の式（５）のようなＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデル
Figure 2023075072000007
また、最小のサポートベクトルマシンを適用する前記ステップＳ４は、
学習サンプルを整理して、正規化処理を行い、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速を、６つの入力ベクトルとして選び出すステップを含み、
ａ及びｂは、推定関数を用いて計算し、前記推定関数は以下の式（６）であり、

Figure 2023075072000008
Figure 2023075072000009

Figure 2023075072000010
Figure 2023075072000011
すでに風速がある時系列を訓練に代入することにより、これまでの風速からｆｘが得られているのに対し、逆押し時刻に対応する前のＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速それぞれの対応変数も既知であるため、推定関数を用いてａ及びｂを計算して得、
パラメータａ及びｂを式（６）に代入し、ＬＳ－ＳＶＭモデルの具体的な形式を確定し、
前記６つの入力ベクトルを前記ＬＳ－ＳＶＭモデルに入力し、ｔ時刻の予測風速を得る。
また、時間変数ｔが遷移するに従って、前記ステップＳ４を連続的に繰り返し、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ部分の予測の風速時系列を得る。
本発明は、さらに前記の方法を実現するために用いられる、ローリング時系列及びサポートベクトルマシンに基づく短期風速予測システムに関し、前記システムは、
風速時系列を入力する入力モジュール、
位相空間の再構成に用いられる再構成モジュール、
ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するために用いられるローリング予測モジュール、
最小のサポートベクトルマシンを適用するために用いられる適用モジュール、及び
風速の予測値を出力する出力モジュールを含む。
本発明は、さらにコンピュータプログラムが記憶されているコンピュータ可読記憶媒体に関し、プロセッサによってこのプログラムが実行される場合、前記のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法が実現される。
本発明は、ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法であり、空間位相の再構成法を用いて、風速のカオス特性を再現した後に、時系列モデルとサポートベクトルマシンの反応風速系列の線形と非線形特性を組合わせて、風速予測の精度と信頼性を高める有益な効果を有する。
本発明の方法を示すフローチャートである。 本発明のＡＲＭＡＸ－ＡＲ－ＧＡＲＣＨローリング予測とサポートベクトルマシンモデルの組合せを示す模式図である。 本発明のＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測モデルを示す模式図である。 本発明のＬＳ－ＳＶＭの流れを示す模式図である。
以下、添付図面を参照して本発明をさらに説明する。以下の実施形態は、本発明の技術的解決策をより明確に説明するためにのみ使用され、これによって本発明の保護範囲を制限することはできない。以下の詳細な説明は、すべて例示的なものであり、本願にさらなる説明を提供することが意図されていることに留意されたい。
特に指定がない限り、本明細書で使用されるすべての技術用語及び科学用語は、本願に属する技術分野の一般技術者が通常に理解しているものと同じ意味を有する。なお、本明細書で使用される用語は、本願による例示的な実施形態を限定することを意図するのではなく、特定の実施形態を説明するためにのみ使用される。ここで使用されるように、文脈上特に明示されない限り、単数形は複数形も含むことを意図している。また、本明細書で用語「含む」及び／または「からなる」を使用する場合、特徴、ステップ、操作、デバイス、コンポーネント及び／またはこれらの組み合わせが存在することを示すことも理解されたい。
ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法は、まず、風速系列のカオス特性に対して、元のデータの空間位相を再構成し、得られた変数入力は、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ時系列モデルを用いて、風速の線形中間値を予測し、風速系列の線形関係を明らかにする。同時に、ローリングの予測戦略を用いて、モデルのリアルタイム性を保証する。風速線形及び空間位相を再構成して得られた変数を、最小二乗のサポートベクトルマシンの入力とし、最後の予測風速を得て、風速系列の非線形関係を明らかにする。
図１を参照すると、本発明のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法は、
風速時系列を入力するステップＳ１、
位相空間を再構成するステップＳ２、
ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するステップＳ３、
最小のサポートベクトルマシンを適用するステップＳ４、及び
風速の予測値を出力するステップＳ５を含む。
風速時系列を入力する前記ステップＳ１は、元の風速時系列Ｐｂを入力するステップを含む。
本発明は、空間位相の再構成アルゴリズムを用いて、風速データにカオス特性処理を行い、空間位相を再構成する中間値系列を得、
空間位相を再構成する中間値系列を前処理し、
ＡＩＣを基準探索とした方法を用いて、ＡＲＭＡモデルの次数項（ｐ、ｑ）の探索範囲をそれぞれ（２－２１）、（２－１２）に定め、ＡＩＣに最小を取る次数を選択し、次数項パラメータｐ、ｑの解を求める。
ＡＩＣが最小値を取るように探索するＡＲ部分次数を採用する。
ＧＡＲＣＨの次数選択は（１、１）を選択し、
ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ時系列モデルを確認するため、全てのパラメータを構築する。
図２に示すように、前処理後の風速系列を入力として、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ時系列モデルの出力を求め、出力変数を予測風速にまとめる。
Figure 2023075072000012
図３に示すように、時間変数ｔの遷移に従って、連続的に繰り返し、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ部分の予測の風速時系列を得る。
風速時系列には、非線形とランダム性の特徴があり、元の系列Ｐｂを直接利用して訓練と予測を行うと、モデルの予測精度に影響を与えるため、位相空間を再構成する方法を導入する。
位相空間を再構成する前記ステップＳ２は、
前記元の風速時系列Ｐｂに対して、以下の式（１）に従って位相空間を再構成し、空間位相を再構成する中間値系列を得るステップを含み、

Figure 2023075072000013
得られた前記空間位相を再構成する中間値系列に対して前処理を行うことによって、データの処理を容易にする。
また、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用する前記ステップＳ３は、
Figure 2023075072000014
また、時系列モデルＡＲＭＡ（ｐ、ｑ）は、定常系列を予測するために用いられる以下の式（２）であり、

Figure 2023075072000015
前記時系列モデルＡＲＭＡ（ｐ、ｑ）は、モデル残差の平均値が０であり、残差の分散が定数であるという仮定を前提とする。
また、ＧＡＲＣＨモデルは以下の式（３）であり、
Figure 2023075072000016
ＧＡＲＣＨモデルの増加は、異分散関数のｍ次相関性を考慮したものであり、風速系列の異分散性を効果的に捉まえることができる。
Figure 2023075072000017
このように構成されたＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルは以下の式（４）であり、
Figure 2023075072000018
前記時系列モデルＡＲＭＡ（ｐ、ｑ）と前記ＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルを連立させ、以下の式（５）のようなＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデル
Figure 2023075072000019
ＡＩＣを基準探索とした方法を採用し、ＡＲＭＡモデルの３つの次数項（ｐ、ｑ）の探索範囲をそれぞれ（２－２１）、（２－１２）に定め、ＡＩＣが最小値を取る次数を選択する。
この部分ＡＲの次数決定に対しては、ＡＩＣが最小値を取るように探索する次数を採用する。ＧＡＲＣＨの次数選択は（１、１）を選択する。
最小二乗のサポートベクトルマシン（ＬＳ－ＳＶＭ）は、サポートベクトルマシンの汎化能力が強く、グローバル最適という利点があり、非線形関係の予測に用いることに適している。ラジアル基関数の性能が優れているため、核関数としてラジアル基関数
Figure 2023075072000020
最小二乗のサポートベクトルマシンのモデル化ステップを、図４のように示す。
最小のサポートベクトルマシンを適用する前記ステップＳ４は、
学習サンプルを整理して、正規化処理を行い、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速を、６つの入力ベクトルとして選び出すステップを含む。出力変数は、即ちｔ時刻に対応する予測の風速である。
訓練の前に、交差検証法を通じて、核パラメータ及び正規化パラメータの２つのパラメータを連続的に調整することによって、交差検証誤差が最小のパラメータの組合せを保留する。

Figure 2023075072000021
ＫＫＴ条件に基づいて、
Figure 2023075072000022
すでに風速がある時系列を訓練に代入することにより、これまでの風速からｆｘが得られているのに対し、逆押し時刻に対応する前のＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速それぞれの対応変数も既知であるため、推定関数を用いてａ及びｂを計算して求める。
式（６）にパラメータａ及びｂを代入し、ＬＳ－ＳＶＭモデルの具体的な形式を確定し、
ｔ時刻の風速を予測する場合、前記６つの入力ベクトルを前記ＬＳ－ＳＶＭモデルに入力し、ｔ時刻の予測風速を得る。
時間変数ｔが遷移するに従って、前記ステップＳ４を連続的に繰り返し、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ部分の予測の風速時系列を得る。
本発明は、さらに前記の方法の実現に用いられる、ローリング時系列及びサポートベクトルマシンに基づく短期風速予測システムに関し、前記システムは、
風速時系列を入力する入力モジュール、
位相空間の再構成に用いられる再構成モジュール、
ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するために用いられるローリング予測モジュール、
最小のサポートベクトルマシンを適用するために用いられる適用モジュール、及び
風速の予測値を出力する出力モジュールを含む。
本発明は、さらにコンピュータプログラムが記憶されているコンピュータ可読記憶媒体に関し、プロセッサによってこのプログラムが実行される場合、前記のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法が実現できる。
本発明は、まず風速系列のカオス特性に対して、元のデータの空間位相を再構成し、得られた変数入力は、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ時系列モデルを用いて、風速線形の中間値を予測し、風速系列の線形関係を明らかにする。同時に、ローリングの予測戦略を用いて、モデルのリアルタイム性を保証する。風速線形及び空間位相を再構成して得られた変数を、最小二乗のサポートベクトルマシンの入力とし、最後の予測風速を得、風速系列の非線形関係を明らかにする。
以上は、本発明の好ましい実施形態に過ぎず、当業者にとっては、本発明の技術原理を逸脱することなく、いくつかの改善及び変形を行うことができ、これらの改善及び変形も本発明の保護範囲と見なすべきであることを指摘すべきである。

Claims (10)

  1. ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法であって、
    風速時系列を入力するステップＳ１、
    位相空間を再構成するステップＳ２、
    ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するステップＳ３、
    最小のサポートベクトルマシンを適用するステップＳ４、及び
    風速の予測値を出力するステップＳ５を含むことを特徴とする、ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  2. 風速時系列を入力する前記ステップＳ１は、
    元の風速時系列Ｐｂを入力するステップを含み、
    位相空間を再構成する前記ステップＳ２は、
    前記元の風速時系列Ｐｂに対して、以下の式（１）に従って位相空間を再構成し、空間位相を再構成する中間値系列を得るステップを含み、

    Figure 2023075072000023
    得られた前記空間位相を再構成する中間値系列に対して前処理を行うことによって、データの処理を容易にすることを特徴とする、ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  3. ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用する前記ステップＳ３は、
    Figure 2023075072000024
  4. 時系列モデルＡＲＭＡ（ｐ、ｑ）は、定常系列を予測するために用いられる以下の式（２）であり、
    Figure 2023075072000025
    前記時系列モデルＡＲＭＡ（ｐ、ｑ）は、モデル残差の平均値が０であり、残差の分散が定数であるという仮定を前提とすることを特徴とする、請求項３に記載のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  5. ＧＡＲＣＨモデルは以下の式（３）であり、
    Figure 2023075072000026
  6. ＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルは、以下の式（４）であり、
    Figure 2023075072000027
    前記時系列モデルＡＲＭＡ（ｐ、ｑ）と前記ＡＲ（ｓ）－ＧＡＲＣＨ（ｍ、ｎ）モデルを連立させ、以下の式（５）のようなＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデル
    Figure 2023075072000028
    （５）が得られることを特徴とする、請求項５に記載のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  7. 最小のサポートベクトルマシンを適用する前記ステップＳ４は、
    学習サンプルを整理して、正規化処理を行い、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速を、６つの入力ベクトルとして選び出すステップを含み、
    ａ及びｂは、推定関数を用いて計算し、前記推定関数は以下の式（６）であり、
    Figure 2023075072000029
    核関数としてラジアル基関数
    Figure 2023075072000030
    （１２）を用い、その推定関数は上記式（６）であり、
    すでに風速がある時系列を訓練に代入することにより、これまでの風速からｆｘが得られているのに対し、逆押し時刻に対応する前のＡＲＭＡ－ＡＲ－ＧＡＲＣＨモデルによって予測したｔ－３、ｔ－２、ｔ－１、ｔ時刻の風速、及びｔ－１、ｔ－２時刻の実際の風速それぞれの対応変数も既知であるため、推定関数を用いてａ及びｂを計算して求め、
    パラメータａ及びｂを式（６）に代入し、ＬＳ－ＳＶＭモデルの具体的な形式を確定し、
    前記６つの入力ベクトルを前記ＬＳ－ＳＶＭモデルに入力し、ｔ時刻の予測風速を得ることを特徴とする、請求項１に記載のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  8. 時間変数ｔが遷移するに従って、前記ステップＳ４を連続的に繰り返し、ＡＲＭＡ－ＡＲ－ＧＡＲＣＨ部分の予測風速時系列を得ることを特徴とする、請求項７に記載のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法。
  9. 請求項１～８に記載方法の実現に用いられるローリング時系列及びサポートベクトルマシンに基づく短期風速予測システムであって、前記システムは、
    風速時系列を入力する入力モジュール、
    位相空間の再構成に用いられる再構成モジュール、
    ＡＲＭＡ－ＡＲ－ＧＡＲＣＨローリング予測を適用するために用いられるローリング予測モジュール、
    最小のサポートベクトルマシンを適用するために用いられる適用モジュール、及び
    風速の予測値を出力する出力モジュールを含むことを特徴とする、ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測システム。
  10. コンピュータプログラムが記憶されているコンピュータ可読記憶媒体であって、
    プロセッサによってこのプログラムが実行される場合、請求項１～８に記載のローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法が実現されることを特徴とする、コンピュータ可読記憶媒体。
JP2022185214A 2021-11-18 2022-11-18 ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステム Active JP7441553B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111370513.8 2021-11-18
CN202111370513.8A CN114595858A (zh) 2021-11-18 2021-11-18 基于滚动时间序列和支持向量机的短期风速预测方法和系统

Publications (2)

Publication Number Publication Date
JP2023075072A true JP2023075072A (ja) 2023-05-30
JP7441553B2 JP7441553B2 (ja) 2024-03-01

Family

ID=81803146

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022185214A Active JP7441553B2 (ja) 2021-11-18 2022-11-18 ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステム

Country Status (2)

Country Link
JP (1) JP7441553B2 (ja)
CN (1) CN114595858A (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116662766A (zh) * 2023-08-01 2023-08-29 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备
CN117056692A (zh) * 2023-10-09 2023-11-14 山东芯赛思电子科技有限公司 一种SiC基电机驱动器件的老化预测方法
CN117056735A (zh) * 2023-10-13 2023-11-14 云南电投绿能科技有限公司 一种风电场的短期风速预测方法、装置、设备及存储介质
CN117933363A (zh) * 2024-03-08 2024-04-26 广东工业大学 基于样本迁移的新建风电场短期风电功率预测方法和系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102102626A (zh) * 2011-01-30 2011-06-22 华北电力大学 风电场短期功率预测方法
US20170016430A1 (en) * 2015-07-17 2017-01-19 General Electric Company Systems and methods for improved wind power generation
CN107292446A (zh) * 2017-07-03 2017-10-24 西南交通大学 一种基于考虑分量关联性小波分解的混合风速预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102102626A (zh) * 2011-01-30 2011-06-22 华北电力大学 风电场短期功率预测方法
US20170016430A1 (en) * 2015-07-17 2017-01-19 General Electric Company Systems and methods for improved wind power generation
CN107292446A (zh) * 2017-07-03 2017-10-24 西南交通大学 一种基于考虑分量关联性小波分解的混合风速预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIANG, Y ET AL.: ""Short-term wind speed prediction: Hybrid of ensemble empirical mode decomposition, feature selectio", ENERGY CONVERSION AND MANAGEMENT [ONLINE], vol. 144, JPN6023017625, 2017, pages 340 - 350, XP085029828, ISSN: 0005054710, DOI: 10.1016/j.enconman.2017.04.064 *
飯坂 達也 ほか: ""風力発電の出力予測の現状"", 電気学会研究会資料, vol. MES-11-003, JPN6020025221, 2011, pages 13 - 18, ISSN: 0005054711 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116662766A (zh) * 2023-08-01 2023-08-29 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备
CN116662766B (zh) * 2023-08-01 2023-10-03 浙江大学 基于数据二维重构的风速预测方法及装置、电子设备
CN117056692A (zh) * 2023-10-09 2023-11-14 山东芯赛思电子科技有限公司 一种SiC基电机驱动器件的老化预测方法
CN117056692B (zh) * 2023-10-09 2024-01-23 山东芯赛思电子科技有限公司 一种SiC基电机驱动器件的老化预测方法
CN117056735A (zh) * 2023-10-13 2023-11-14 云南电投绿能科技有限公司 一种风电场的短期风速预测方法、装置、设备及存储介质
CN117933363A (zh) * 2024-03-08 2024-04-26 广东工业大学 基于样本迁移的新建风电场短期风电功率预测方法和系统
CN117933363B (zh) * 2024-03-08 2024-06-11 广东工业大学 基于样本迁移的新建风电场短期风电功率预测方法和系统

Also Published As

Publication number Publication date
CN114595858A (zh) 2022-06-07
JP7441553B2 (ja) 2024-03-01

Similar Documents

Publication Publication Date Title
JP2023075072A (ja) ローリング時系列及びサポートベクトルマシンに基づく短期風速の予測方法及びシステム
Song et al. Hourly heat load prediction model based on temporal convolutional neural network
Mi et al. Wind speed prediction model using singular spectrum analysis, empirical mode decomposition and convolutional support vector machine
Wang et al. Wind speed forecasting based on the hybrid ensemble empirical mode decomposition and GA-BP neural network method
Nguyen et al. Hourly day ahead wind speed forecasting based on a hybrid model of EEMD, CNN-Bi-LSTM embedded with GA optimization
Wang et al. Short-term wind speed forecasting based on information of neighboring wind farms
Suo et al. Wind speed prediction by a swarm intelligence based deep learning model via signal decomposition and parameter optimization using improved chimp optimization algorithm
Liu et al. Short-term load forecasting of multi-energy in integrated energy system based on multivariate phase space reconstruction and support vector regression mode
Wang et al. An innovative hybrid approach for multi-step ahead wind speed prediction
Liu et al. Super resolution perception for smart meter data
CN112651290A (zh) 一种水电机组振动趋势预测方法和系统
Zhang et al. A hybrid prediction model for forecasting wind energy resources
Li et al. Multi-step ahead wind speed forecasting approach coupling maximal overlap discrete wavelet transform, improved grey wolf optimization algorithm and long short-term memory
Wang et al. Short-term wind speed prediction using support vector regression
Abedinia et al. A multi-level model for hybrid short term wind forecasting based on SVM, wavelet transform and feature selection
Zhang et al. Interval prediction of ultra-short-term photovoltaic power based on a hybrid model
Wang et al. A novel wind power prediction model improved with feature enhancement and autoregressive error compensation
CN104239689B (zh) 基于优化相关向量机的短期风速预测方法
Wu et al. Interpretable wind speed forecasting with meteorological feature exploring and two-stage decomposition
Wang et al. Research on ultra-short term forecasting technology of wind power output based on various meteorological factors
Zhu et al. Wind Speed Short-Term Prediction Based on Empirical Wavelet Transform, Recurrent Neural Network and Error Correction
Zeng et al. Short-term wind power prediction based on the combination of numerical weather forecast and time series
Chen et al. Chaotic wind speed series forecasting based on wavelet packet decomposition and support vector regression
Mahin et al. Advanced Forecasting of Wind Energy Generation Through Integration of AE-CLSTM
Liu et al. Wind power prediction based on the stacking model of xgboost and random forest

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221118

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20221202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230509

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20230808

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20231010

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231107

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240206

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240209

R150 Certificate of patent or registration of utility model

Ref document number: 7441553

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150