風力発電量予測システム、そのプログラム
【課題】一般的な気象予報データと、データ量が少ない補正用データとに基づいて、各風力発電設備毎の風力発電量を的確に予測できる。
【解決手段】気象予報入力機能部１１は、外部から一般的な気象予報情報としての風向・風速予報値を取得する。補正テーブル記憶部１６には予め風速予報値を補正するための補正用データが記憶される。パワーカーブ記憶部１３には予め風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶される。発電量計算機能部１４は、上記風向・風速予報値と補正用データとパワーカーブとを用いて、加重平均またはファジィ推論によって、風力発電量予測対象の任意の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【解決手段】気象予報入力機能部１１は、外部から一般的な気象予報情報としての風向・風速予報値を取得する。補正テーブル記憶部１６には予め風速予報値を補正するための補正用データが記憶される。パワーカーブ記憶部１３には予め風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶される。発電量計算機能部１４は、上記風向・風速予報値と補正用データとパワーカーブとを用いて、加重平均またはファジィ推論によって、風力発電量予測対象の任意の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、風速予測から風力発電量予測を行う技術に関する。
【背景技術】
【０００２】
風力発電量予測値は、予め作成される、パワーカーブ等と呼ばれる風速と発電量との関係データを用いて、風速予報等による風向・風速予測に基づいて求められる。このパワーカーブは、一般的に風力発電装置の設計値により決定される。
【０００３】
風力発電量に影響を与える主要因は、風速であるため、風力発電予測では実際には風向・風速予測が重要であり、風速等から発電量を求める処理自体は、それほど重要ではない（上記の通り、設計値により決定されたパワーカーブを用いて計算すれば済む）。つまり、風車が受ける風速等を正確に予測できれば、風力発電量予測値が的確に計算できる。
【０００４】
ここで、風速等を予測する従来技術としては、例えば特許文献１，２，３，４に記載の従来技術が知られている。
特許文献１には、風速計および風向計によって風速や風向を観測する装置である観測ポールと、３次元風速計と、予測用コンピュータ等を有する風予測システムが開示されている。特許文献１では、３次元風速計によって参照点の風速および風向を測定するのは、各季節毎に１週間ずつの計４週間で済むようにする手法を提案しており、これによって高価な３次元風速計の利用効率を高め、コストを削減できる。
【０００５】
特許文献２の風力発電出力予測装置は、複数の物理モデルと複数の統計モデルを組み合わせることにより、風力発電出力をより正確に予測できるものである。
特許文献３では、風速予報値を風車高さ相当に補正し、これをパワーカーブを用いて発電量に変換している。
【０００６】
特許文献４では、レーダにより風速を詳細に観測し、ＧＰＶ気象予報を補正するデータベースを持ち、風速を補正している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００５−１００８２号公報
【特許文献２】特開２００７−２３３６３９号公報
【特許文献３】特開２００９−１６７８４８号公報
【特許文献４】特開２００９−１３８５２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
パワーカーブは、実験室や工場などのような理想的な環境で計測され、風力発電設備に関する風速と発電量との関係を正確に反映させたものである。しかし、屋外に設置された実際の風力発電設備では、しばしばパワーカーブどおりの発電量を示さない。その理由は主に以下２つであり、何れの場合も、地形や障害物などの影響で、実際に風車が受ける風速が（たとえ予報自体は的中していても）予報通りにならないことで、発電量予測が外れることになる。
【０００９】
理由１：地形が複雑な場合（近くに建物や山がある場合など）、風向によっては風速が強くなる場合や、弱くなる場合がある。たとえば、いくら風速予報値が大きくても山などの障害物側から吹いている場合は、風車にあたる実際の風速は弱くなるため、風力発電設備の発電量を過大に予測することになる。
【００１０】
理由２：ウィンドファームの場合、風上側風車が稼働していた場合、風下側の風速を弱めるため、風下側風車が受ける実際の風速が弱まる。つまり、風速予測が正確であったとしても風下側風力発電設備の発電量を過大に予測することになる。その一方で、風上側風車がメンテナンス等で停止している場合には、風下側の風速は殆ど弱まることはない。
【００１１】
従来技術として、特許文献１の発明は、観測ポールと３次元風速計とによる実測が必要となる。
特許文献２の発明は、詳細な物理モデルを使用するため、上述理由１、理由２による問題の解決を期待できる有望な方法であるが、詳細な地形情報が必要であり、かつ膨大な計算を必要する欠点がある。
【００１２】
特許文献３の発明は、気象予報で得られる風速を、風車高さ相当に補正するものであり、上記理由１、理由２に係わる問題を解決できるものではない。
特許文献４の発明は、風向ごとにデータベースを持つため、上述理由１、理由２による問題の解決を期待できる方法である。しかし、特殊なレーダによる風速観測が必要である点や、詳細な予測を行う場合には細かな風向ごとにデータベースを持つ必要があり、膨大な観測が必要になる等の欠点がある。つまり、エンジニアリング・コストが膨大になる欠点がある。
【００１３】
本発明の課題は、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減でき、更に、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる風力発電量予測システム等を提供することである。
【課題を解決するための手段】
【００１４】
本発明の風力発電量予測システムは、任意の風力発電設備の発電量を予測する風力発電量予測システムであって、外部から風向・風速予報値を取得する風向・風速予報値取得手段と、予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段とを有する。
【００１５】
上記風力発電量予測システムにおいて、例えば、前記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルであり、前記風力発電量算出手段は、前記風速補正テーブルを参照して、前記風向予報が示す風向に係わる１以上の前記所定の風向毎に、前記風速予報値に対応する前記補正データを用いて風速補正値を求めると共に、前記パワーカーブを用いて、前記求めた１以上の風速補正値毎に対応する発電量を求め、複数の発電量が求められた場合には、該複数の発電量を前記風向予報に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求める。
【００１６】
あるいは、上記風力発電量予測システムにおいて、例えば、前記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、前記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数であり、前記風力発電量算出手段は、前記風速補正テーブルを参照して、前記各風向範囲毎に、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求めると共に、前記メンバシップ関数を参照して、前記各風向範囲毎に、前記風向予報値に応じた前記加重比を求め、前記求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、前記風速予報値を補正した風速補正値を算出し、該算出した風速補正値と前記パワーカーブとによって、前記風力発電設備が受ける風に応じた風力発電量予測値を求める。
【００１７】
上記風力発電量予測システムでは、一般的な風向・風速予報値等と補正用データ等に基づいて、各風力発電設備がそれぞれ実際に受ける風速に応じた風力発電量予測値を求めることが可能となり、地形や障害物等による影響に対応した的確に風力発電量予測を行うことができる。また、加重平均またはファジィ推論を用いることで、補正用データのデータ量を削減できる。例えば全ての風向ではなく上記複数の所定の風向に関してのみ補正データを登録しておけば済むようになる。
【００１８】
また、更に、前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備による影響も考慮して、風力発電量予測値を求めることで、更に的確に風力発電量予測値を求めることができるようにしてもよい。
【発明の効果】
【００１９】
本発明の風力発電量予測システムによれば、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減できる。更に、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる。
【図面の簡単な説明】
【００２０】
【図１】本例の風力発電量予測システムの構成図である。
【図２】パワーカーブの一例を示す図である。
【図３】（ａ）〜（ｃ）は実施例１に係るパワーカーブや風速補正テーブルの一例である。
【図４】（ａ）〜（ｄ）は実施例２に係る風速補正テーブルの一例である。
【図５】（ａ），（ｂ）は、実施例２に係るメンバシップ関数の一例である。
【図６】（ａ），（ｂ）は、実施例３に係るパワーカーブ／風速補正テーブルの一例である。
【図７】（ａ）、（ｂ）は実施例４に係るパワーカーブや発電量補正テーブルの一例である。
【図８】実施例５の効果を示す図である。
【発明を実施するための形態】
【００２１】
以下、図面を参照して本発明の実施の形態について説明する。
図１は、本例の風力発電量予測システムの構成図である。
図示の例の風力発電量予測システム１０は、気象予報入力機能部１１、風力発電設備運転情報取得部１２、パワーカーブ記憶部１３、発電量計算機能部１４、予測値格納部１５、補正テーブル記憶部１６の各種機能部を有する。
【００２２】
尚、風力発電量予測システム１０は、例えば一般的な汎用コンピュータ（パソコン、サーバ装置等）で実現される。つまり、風力発電量予測システム１０は、ハードウェア的には一般的な汎用コンピュータの構成であってよく、特に図示はしないが、例えばＣＰＵ等の演算プロセッサ、メモリやハードディスク等の記憶装置、ネットワークを介して外部のコンピュータ装置と通信可能とする通信部、液晶ディスプレイ等の表示部、キーボードやマウス等の入力部などを有している。
【００２３】
また、上記不図示の記憶装置には予め所定のアプリケーションプログラムが記憶されている。上記不図示の演算プロセッサがこのアプリケーションプログラムを読出し・実行することにより、上記各種機能部１１〜１６の処理機能が実現される。
【００２４】
気象予報入力機能部１１は、外部で作成された気象予報情報（特に、風向・風速予測値など）を、この気象予報情報を記憶している外部のコンピュータ装置等からネットワークを介して入力する入力機能部である（あるいは、上記気象予報情報を見たユーザ等が、上記キーボード等を操作して気象予報情報を入力するものであってもよい）。
【００２５】
入力する気象予報情報は、一般的な気象予報データであり、例えば気象庁やその他の気象業者の発表した気象予報情報そのものであってもよいし、何らかの別の方法で作成された気象予報情報でもよい。何れにしても、一般的な気象予報データは、各風力発電設備（その風車）が実際に受ける風の風速予測値等を示すようなものではなく、ある程度広いエリア毎の風向・風速予測値を示すものである。気象予報入力機能部１１は、例えば風力発電設備が設置されているエリアに関する風向・風速予測値を入力する。
【００２６】
尚、風向は、例えば南、南南西、南西等の１６方位で表すものであってもよいし、真北を基準（０°）にして、東が９０°、南が１８０°、西が２７０°等のように時計周りに表す３６０方位であってもよい。
【００２７】
風力発電設備運転情報取得部１２は、風力発電設備の運転スケジュール情報を取得する機能部である。この運転スケジュール情報は、各風力発電設備の現在や未来の運転状態（運転（ＯＮ）／停止（ＯＦＦ））を示す運転計画データ等である。尚、ここでは、運転停止（ＯＦＦ）は、メンテナンス等に伴って起こるものであり、故障等による運転停止は関係ない。
【００２８】
上記運転スケジュール情報は、本システム１０内に記憶されていてよいし、キーボード等で直接的に入力されてもよいし、不図示の外部装置等に記憶されていてもよい。外部装置に記憶されている場合には、当然、風力発電設備運転情報取得部１２は、不図示のネットワーク等を介して、外部装置から上記運転スケジュール情報等を取得することになる。
【００２９】
パワーカーブ記憶部１３は、予めパワーカーブを記憶しているデータベースである。パワーカーブは、各風力発電設備について風速と発電量の関係を示すデータであり、上記のように従来から存在するものであり、予め任意に作成されて記憶されている。
【００３０】
パワーカーブは、その形式は、テーブル形式であってもよいし、ｎ次関数の形で持っていてもよい。また、１つのパワーカーブであってもよいし、風向毎にパワーカーブをもっていてもよい。風向計を備えており風向きに応じてその向きを変えるタイプの風力発電設備であれば、風向きに関係なく風が吹いてくる方向を向くので（常にほぼ正面から風を受けるので）パワーカーブは１つであってよい。本例では、パワーカーブ自体は１つであるものとする。
【００３１】
図２に、パワーカーブの一例を示す。図示のパワーカーブは、任意の風力発電設備に係る特性データであって、横軸が風速、縦軸が発電量となっており、風速と発電量との関係が示されている。尚、パワーカーブ自体は既存のものであり、これ以上詳細には説明しないものとする。
【００３２】
発電量計算機能部１４は、本システム１０の主要機能であり、その処理内容については後に詳しく説明するが、基本的には上記各機能部１１，１２から得られる入力情報等や上記パワーカーブや後述する補正テーブル等に基づいて、現在から（あるいは未来ｎ時間先から）未来ｍ時間先までの風力発電設備の発電量の予測値を算出・出力する。これは、例えば、所定時間単位の（例えば１時間単位の）各予測対象時間帯毎に、風力発電設備の発電量の予測値を算出・出力する。
【００３３】
風力発電設備が複数ある場合には、各風力発電設備毎にそれぞれ、現在から未来ｍ時間先までの発電量の予測値を算出・出力する。
発電量計算機能部１４によって算出・出力される上記「ｍ時間先までの風力発電設備の発電量の予測値」は、当該予測値を格納するデータベースである予測値格納部１５に格納される。
【００３４】
補正テーブル記憶部１６には、予め、風速予報値を補正するための補正テーブルが記憶されている。これについては、後に詳しく説明する。
上記風力発電設備の運転スケジュール情報は、発電所において複数の風車（風力発電設備）がある場合に、特に必要になるものである。詳しくは後述する。
【００３５】
以下、上記風力発電量予測システム１０の処理機能について、各種実施例を提案・説明する。尚、以下に説明する実施例１〜５における風力発電量予測システム１０運用時の各種処理は、風向・風速予報値や運転スケジュール等を取得する処理以外は、基本的には全て発電量計算機能部１４が実行するものと考えてよい。
【００３６】
まず、実施例１について説明する。
気象予報入力機能部１１からは、例えば１時間先〜ｍ時間先までの１時間毎の、発電所のある代表地点、もしくは近辺の風向・風速予報値が得られる。換言すれば、上記風力発電設備が設置されているエリアの風向・風速予報値が得られるものである。
【００３７】
風力発電設備運転情報取得部１２からは、各風力発電設備の運転スケジュールが得られる。つまり、風力発電設備に関しては、定期点検などのメンテナンス等が理由で運転を停止する期間（よって、事前に決めておくことが可能）があり、この様な運転停止期間（ＯＦＦ期間）と、通常の運転期間（ＯＮ期間）とを示す運転スケジュールが、予め作成されて記憶されているので、これを取得する。これは、勿論、現在の運転状態（ＯＮ／ＯＦＦ）だけでなく、未来の運転予定（ＯＮ／ＯＦＦ）も得られる。例えば、１時間先〜ｍ時間先までの運転スケジュールを得る。尚、運転がＯＦＦ（停止）の期間は、発電量は‘０’である。
【００３８】
基本的には、発電所内にある全ての風力発電設備について、それぞれ、上記運転スケジュールが作成されて記憶されているが、風力発電設備運転情報取得部１２が取得するのは、予測処理対象の風力発電設備の運転スケジュールや、その風上側風力発電設備（風上側風車）の運転スケジュールのみであってよい。また、必ずしも運転スケジュールを取得する必要はない。
【００３９】
尚、例えば、各風力発電設備には予め識別用ＩＤ等が付与されており、特に図示しないが、各風力発電設備毎に、その風上側風力発電設備の識別用ＩＤが登録されているものとする。また、各識別用ＩＤに対応付けてその風力発電設備の運転スケジュールが管理されている。これによって、例えばユーザが任意の風力発電設備を予測処理対象に指定すると、その風上側風力発電設備を判別できると共に、これらの運転スケジュールを取得できる。
【００４０】
また、尚、現在時点の運転状態はＯＮ／ＯＦＦの２値の何れかの状態しかないが、未来の運転予定は、例えば仮に１時間中に６分間停止する予定であれば「稼働率が９０％」等のように、アナログ量として取得してもよい（尚、これは、１時間単位で取得する場合であり、また運転スケジュールは分単位で管理されているものとするが、この例に限らない）。
【００４１】
尚、本例では、上記のように入力・取得するデータは、１時間先〜ｍ時間先まで１時間単位で扱うと共に、これに合わせて発電量予測値も１時間先〜ｍ時間先まで１時間単位で扱うものとする。
【００４２】
ここで、上記運転状態におけるＯＦＦ状態は、あくまでもメンテナンス等で発電できない状態であり、風速が‘０’で発電できない状態とは異なる。
尚、上記風上側風力発電設備（風上側風車）は、予測対象の風力発電設備（対象風車）の風上側にある他の風力発電設備であるが、風上側や風上側風車や風上側風力発電設備とは、ある特定の風向きのときに対象風車の風上となることを意味し、風向きに関係なく常に風上側となるようなことを意味するものではない（換言すれば、風上側風車とは、風上側となったときに影響を及ぼす（対象風車が受ける風を弱める）可能性がある、他の風車を意味するものであり、例えば対象風車の近傍に設けられる他の風車等である）。
【００４３】
また、パワーカーブは、その一例を図２に示したが、例えば図３（ａ）に示すようにテーブル形式であってもよい。図示の例では、１ｍ単位の各風速（ｍ）毎に、対応する発電量（ｋｗ）が格納されているが、この例に限らない。尚、既に従来技術で述べたように、パワーカーブの各値は、例えばメーカが工場等で計測した値である。
【００４４】
また、図３（ｂ）、（ｃ）には、補正テーブル記憶部１６に予め記憶されている風速補正テーブルの具体例を示す。
風速補正テーブルは、一般的な気象情報としての風速予報値を、その風車が受ける風速へと補正するためのテーブルである。尚、以下、風速補正テーブルを簡略化して単に「補正テーブル」と記す場合もあるものとする。
【００４５】
補正テーブルは、例えば、予め、例えば上記気象庁等による一般的な気象予報データから、例えば「発電所のある代表地点もしくは近辺の風速・風向の予報値または現在値」を取得すると共に、各風力発電設備の設置場所で風速計等を用いて風速等を実測する等して、これらに基づいて人間が作成するものである。
【００４６】
すなわち、例えば気象庁等が発表する風速・風向の現在値（観測値）を取得すると共に、各風力発電設備に設けられた風速計で実測された風速データを取得して、両者の差分（実測データ−風速の現在値）を求めて、この差分を上記風速・風向の現在値と共に記録する。または、実測データを上記風速・風向の現在値に対応付けて記憶するようにしてもよい。
【００４７】
あるいは、例えば気象庁等が発表するｍ時間先の風向・風速予報値を取得し、その後、ｍ時間後に各風力発電設備に設けられた風速計等で風速等を実測する。そして、風速予報値と実測データとの差分（実測データ−風速予報値）を求めて、この差分を上記風向・風速予報値に対応付けて記憶するようにしてもよい。または、実測データを上記風向・風速予報値に対応付けて記憶するようにしてもよい。
【００４８】
尚、上記風向は、予報値等の代わりに、各風力発電設備に設けられた風向計等で実測してもよい。
以上の作業を何度も繰り返すことで補正テーブルを作成できる情報が収集できたら、これに基づいて補正テーブルを作成する。尚、風向きは制御できないので、本例では少なくとも東西南北の４方向の風向きについてデータ収集できるまで、上記データ収集を続けることになる。そして、例えば、東西南北の各方向毎に、上記差分の平均値を求める。例えば、東を例にすると、上記収集したデータから風向が東であるデータを全て抽出して、これらを更に上記風速予報値毎（または風速の現在値毎）に分類して、各風速予報値毎（または風速の現在値毎）に、上記差分値を合計してデータ数で除算することで、東に関する各風速毎の上記差分の平均値が求められる（これが後述する図３（ａ）等に示す“東（補正値）”として風速補正テーブルに格納されることになる）。
【００４９】
ここで、本説明では、上記風上側風力発電設備（風上側風車）を考慮する処理について示すものとし（よって、上記運転スケジュールを取得する必要がある）、この場合には風速補正テーブルは、図３（ｂ）、（ｃ）に示す２種類の風速補正テーブルを用いることになる。
【００５０】
図３（ｂ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
図３（ｃ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
【００５１】
上記風力発電設備運転情報取得部１２で得られる各風力発電設備の運転状態のうち、上記予測対象の風力発電設備の近傍等にある（風速等に影響を及ぼし得る）他の風力発電設備（上記風上側風力発電設備）の運転状態（予測対象時間帯の運転状態）に応じて、図３（ｂ）、（ｃ）の補正テーブルの何れか一方を用いることになる。当然、停止（ＯＦＦ）状態の場合は図３（ｃ）を用い、運転（ＯＮ）状態の場合には図３（ｂ）を用いることになる。
【００５２】
図３（ｂ）、（ｃ）の風速補正テーブルは、データ内容が異なるだけであり、テーブル形式は同一である。よって、以下、特に両者を区別する必要がない場合には、図３（ｂ）等の風速補正テーブルなどと記すものとし、特に区別することなく説明する。
【００５３】
図３（ｂ）等の風速補正テーブルは、１ｍ単位の各風速予報（ｍ）に対応して、各風向（東西南北）毎の補正値（風速予報に対する補正値）が格納されている。尚、図では、「風速予報」と記してあるが、これは「風速」に置き換えても構わない。
【００５４】
例えば、図示の「北（補正値）」は、風向が北に対応する補正値であり、例えば風向予報が“北”で風速予報が２ｍのときには、「北（補正値）」＝‘−１’を使用することになるので、この補正値を用いた補正後の風速は、１ｍ（＝２−１）になる。同様に、例えば、風向予報が“西”で風速が１０ｍの場合には、「西（補正値）」＝‘−２’を使用することになるので、この補正値を用いた補正後の風速は、８ｍ（＝１０−２）になる。
【００５５】
但し、この様な例に限らず、後に説明するように、上記“補正値”の代わりに上記“補正後の風速”が、風速補正テーブルに格納されていてもよい。この場合、上記“補正後の風速”を算出する処理は必要なくなる。
【００５６】
従来であれば、上記の例では、それぞれ、「風速２ｍ」、「風速１０ｍ」によりパワーカーブで発電量を求めることになり、発電量予測が実際よりも過大となる場合があったが、本手法では、それぞれ、「風速１ｍ」、「風速８ｍ」に補正してから（予測対象の風車が実際に受けるものと予測される風速に補正してから）パワーカーブで発電量を求めることになり、的確な発電量予測が行える。
【００５７】
尚、ここでは、補正テーブルの補正値が‘０’であり従って実質的には補正が行われなかった場合でも、風速補正値として扱うものとする。
尚、風向予報が東西南北の４方向以外（例えば北東、北北東等）である場合には、後述する加重平均によって上記発電量予測を行う。逆に言えば、加重平均を用いた演算を行うようにすることで、風速補正テーブルには東西南北の４方向のみ登録すれば済み、データベースのデータ量を削減できる。詳しくは後述する。
【００５８】
ここで、図３（ｂ）、（ｃ）に示す例の風速補正テーブルは、北に風上側風車、西に山などの障害物があることを想定している。そのため、風上側風車の運転状態（ＯＮかＯＦＦか）によってテーブルの内容（数値）が変化するのは、図示の通り北のみである。一方図示の通り、東、南、西に関しては、その数値は、図３（ｂ）と図３（ｃ）とで同じである。つまり、障害物が山やビル等である場合には、その運転状態等は関係ないので（そもそも山には運転状態というものはない）、図３（ｂ）と図３（ｃ）とで同じとある。
【００５９】
一方、障害物が風力発電設備である場合には、運転停止（ＯＦＦ）状態のときには殆ど影響しない（風が弱まらない）が、運転状態（ＯＮ）のときには影響を与えることになる（風が弱まる）。これに応じて、上記図３（ｂ）、（ｃ）のような２種類の風速補正テーブルが、予め設定されている。勿論、これは、予め上記風速補正テーブル作成の為のデータ収集の際に、風上側風車が運転時と停止時とに分けてデータ収集することで、２種類の風速補正テーブルが作成されることになる。
【００６０】
発電量計算機能部１４では、上記風向・風速予報値、風上側風車の運転状態、風速補正テーブル、パワーカーブ等に基づいて、風力発電量予測値を求める。この処理の詳細は、後に具体例を示しながら説明する。
【００６１】
ここで既に述べたように、発電量計算機能部１４は、風向予報が東西南北の４方向の何れかの方向であれば、該当する上記「補正後の風速値」を用いてパワーカーブから該当する発電量を取得すれば済むが、これら以外の風向きであれば、パワーカーブから得た複数の発電量に基づいて所定の算出処理を行う。
【００６２】
後に示す（１）式、（２）式は、風向が“北東”の場合を例にした算出処理例を示している。基本的には、北東であれば北と東、北北東も北と東、南西であれば南と西というように、東西南北の４方向以外の風向きは、東西南北の４方向のうちの２方向に関係するものとなる。これより、特に図示しないが、各風向毎に（但し東西南北は除く）、その風向に関係する東西南北のうちの２方向が、予め記憶されていてもよい。
【００６３】
これより、予報された風向に関係する２方向に係る発電量をパワーカーブから取得し、これら２つの発電量を用いて、後に示す（１）式、（２）式等の算出式による算出処理を行う。
【００６４】
また、発電量計算機能部１４は、予測対象の風車の稼働率に応じた、または／及び、その風上側風車の稼働率に応じた、発電量予測値の算出処理を行うものであってもよい。
発電量計算機能部１４の処理機能について、以下、具体例を用いて説明する。
【００６５】
発電量計算機能部１４は、上記風力発電設備運転情報取得部１２を介して得られる風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、上記パワーカーブと上記風速補正テーブルとに基づいて、風力発電量予測値を求める。
【００６６】
たとえば、パワーカーブと風速補正テーブルが、図３（ａ）、（ｂ）、（ｃ）に示す内容であり、風向予報が北東で、風速予報が１０ｍの場合は、風上側風車の運転状態（ＯＮ／ＯＦＦ）に応じて、それぞれ、下記のようにして求める。
【００６７】
尚、図３（ａ）のパワーテーブルでは示されていないが、ここでは仮に風速８（ｍ）のときの発電量は８０（ｋＷ）であるものとする。
（ａ）風上側風車が運転している場合
上記具体例の場合、上記図３（ｂ）の風速補正テーブルを用いて得られる風速補正値、この風速補正値に応じて図３（ａ）のパワーテーブルから得られる発電量は、所定の風向（本例では東西南北）それぞれについては、下記の通りなる。
【００６８】
北の場合： 風速補正値９ｍ ， 発電量 ９０ｋW
東の場合： 風速補正値１０ｍ ， 発電量１００ｋW
南の場合： 風速補正値１０ｍ ， 発電量１００ｋW
西の場合： 風速補正値８ｍ ， 発電量８０ｋW
但し、これらは参考の為に示しており、実際には風向予報に応じたもの（上記２方向）だけを取得すればよい。すなわち、本例では風向予報は北東であるので、北の発電量（＝９０kW）、東の発電量（＝１００kW）を取得すればよい。そして、これらを用いて後述する加重平均によって風力発電量予測値を求める。但し、風向予報が東西南北の４方向の何れかである場合には、その方向に応じた発電量をそのまま風力発電量予測値とすればよい。
【００６９】
尚、上記風速補正値とは、上記「補正後の風速値」に相当する。ここで、図３(ｂ)等の例では、補正値として上記“差分”を格納しており、従って、「風速補正値（補正後の風速値）＝風速予報値＋補正値」を算出する。例えば、図３(ｂ)の例では風速予報１０ｍに対応する「北（補正値）」は‘−１’であるので、風速補正値（補正後の風速値）＝風速予報＋補正値＝１０−１＝９となり、上記の通り風向が北に対しては風速補正値９ｍとなっている。他も同様である。
【００７０】
尚、この例に限らず、例えば上記“差分”の代わりに、上記実測データが格納されていても構わない。つまり、この場合には、図３（ｂ）等の風速補正テーブルには、各風速予報値に対応付けて上記“風速補正値”が格納されることになる。つまり、例えば風速予報１０ｍに対応する「北（補正値）」は‘−９’となるので、これをそのまま取得すればよく、上記「風速補正値（補正後の風速値）＝風速予報＋補正値」の演算を行う必要はない。
【００７１】
上記の通り、一例では、北の発電量（＝９０kW）、東の発電量（＝１００kW）を取得することになり、また予報された風向は北東であり、北と東の丁度中間であることから、北の発電量（＝９０kW）と東の発電量（＝１００kW）との和を２で除算する（両者の平均を求める）ことで、予報された風向（北東）と風速（１０（ｍ））に基づく、予測対象の風力発電設備が受けると予測される風速に応じた発電量予測値が得られることになる。上記の具体例では下記の（１）式の通りとなる。
【００７２】
発電量予測値；（１００ｋW＋９０ｋW）／２＝９５ｋW・・・（１）式
尚、上記（１）式は、東西南北のうち風向予報に係る２方向に関する発電量を用いて、この２つの発電量を加重平均して求める処理の一例である。
【００７３】
上記のように、北東は北と東の丁度中間であることから、５０：５０の加重比で平均化することになり、
発電量予測値；（１００ｋW×５０＋９０ｋW×５０）／１００
＝（１００ｋW＋９０ｋW）／２
となって、上記（１）式となる。
【００７４】
これより、例えば、北北東であれば、２５：７５の加重比で平均化することになり、従って、
発電量予測値；（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４＝９２．５ｋW・・・（１）’式
となる。
【００７５】
尚、予め例えば１６方向の各風向毎に、その風向に関する加重比が記憶されていてもよい。
上述したように、例えば図３（ｂ）等に示す風速補正テーブルには、予め任意に決められた複数の所定の風向毎に（本例では東西南北毎に）、各風速毎に対応する補正データが格納されている。この風速補正テーブルを参照して、予報された風向に係わる１または複数の風速補正値を求める。つまり、上記所定の風向（東西南北）のうち予報された風向に係わる全ての所定の風向（“関係する所定風向”と呼ぶものとする）について、風速予報値に応じた風速補正値を求める。これは、例えば図３（ｂ）等に示すように上記“差分”である補正値が格納されている例では、上記“関係する所定風向”それぞれについて風速予報値に応じた補正値を取得し更に風速予報値を補正値によって補正することで上記風速補正値を求める。尚、上記の通り、風速補正値そのものが格納されている場合には、これをそのまま取得すればよく、算出する処理は必要ない。
【００７６】
上記“関係する所定風向”は、“予報された風向”が上記一例（北東や北北東）のように上記複数の所定の風向（本例では東西南北）以外である場合には、基本的には上記所定の風向のうちの２方向となる。これは、“予報された風向”が例えば北東や北北東である場合には（その他、東北東の場合でも）北と東となる。従って、この例では、北と東それぞれについて上記風速補正値が求められ、以ってパワーテーブルから対応する発電量が得られることになる。
【００７７】
この様に、風向・風速予報値に応じて複数の発電量が求められる場合には、これら複数の発電量を、上記“予報された風向”に応じた加重比で平均化することで、発電量予測値が算出される（一例が上記（１）式である）。上記の例では、例えば“予報された風向”が北北東である場合には、北が７５、東が２５の加重比が用いられることになる（上記（１）’式となる）。これは、後述する図３（ｃ）等の場合も同様である。
【００７８】
（ｂ）風上側風車が運転していない場合
この場合も、補正テーブルを図３（ｂ）の代わりに図３（ｃ）を用いる点以外は、上記（ａ）風上側風車が運転している場合と同様であり、以下、簡略化して説明する。
【００７９】
この場合、所定の風向（本例では東西南北）それぞれについて発電量は、下記の通りなる。
北の場合： 風速補正値１０ｍ ， 発電量 １００ｋW
東の場合： 風速補正値１０ｍ ， 発電量１００ｋW
南の場合： 風速補正値１０ｍ ， 発電量１００ｋW
西の場合： 風速補正値８ｍ ， 発電量８０ｋW
従って、この例では、予報される風向が北東である場合には、発電量予測値は下記の（２）式の通りとなる。
【００８０】
発電量予測値；（１００ｋW＋１００ｋW）／２＝１００ｋW・・・（２）式
となる。
風上側風車の運転状態を、ＯＮ／ＯＦＦで扱う場合には、例えば予報される風向が北東である場合には、上記（１）式、（２）式の何れかによって、風力発電量予測値を求めることになる。
【００８１】
但し、この例に限らず、例えば風上側風車が、予測対象時間帯のなかで９０％の時間は運転し、残りの１０％の時間は停止している（この例を、稼働率＝９０％とする）場合もあり得る（尚、予測対象時間帯とは、例えば１時間単位で予測する場合には、例えば１〜２時間後、２〜３時間後、・・・、ｍ−１〜ｍ時間後等となる。）
この様な場合、（１）式、（２）式の両方の式による算出を行い、これら２つの算出結果に基づいて、風上側風車の稼働率に応じた加重比で平均化することで、風力発電量予測値を求めることになる。
【００８２】
例えば、上記の例において任意の予測対象時間帯における風上側風車の稼働率が９０％の場合には、上記（１）式、（２）式によって２つの発電量（９５ｋWと１００ｋW）が得られたら、これら２つの発電量を、風上側風車の稼働率に応じた加重比（９０：１０の加重比）で平均化する。勿論、この場合、風上側風車が運転している場合の発電量（＝９０kW）に対しては９０％、風上側風車が運転していない場合の発電量（＝１００kW）に対しては１０％が割当てられる。すなわち、下記の通りとなる。
【００８３】
発電量予測値；（９５ｋW×９０＋１００ｋW×１０）／１００＝９５．５ｋW ・・・（４）式
また、更に、予測対象の風車の稼働率を用いて修正を施して、最終的な予測値とするようにしてもよい。
【００８４】
例えば、仮に、予測対象の風車の稼働率が８０％とした場合には、上述したようにして求めた発電量（（１）式では９５ｋW、（２）式では１００ｋW、あるいは上記９５．５ｋW）に対して、０．８を乗じることで、最終的な発電量予測値とする。
【００８５】
尚、上述した処理と予測値は、未来の任意の時間帯（例えば翌日の２時〜３時の１時間分など）に関する発電量予測処理とその結果である。もし、気象予報が１時間ごとｍ時間先まである場合は、上記処理をｍ回行うことで、ｍ時間先まで求めることができる。
【００８６】
以上説明した実施例１によれば、風速補正テーブルには所定の風向毎の補正値が格納される。この所定の風向の数が少なくても（例えば東西南北の４方位のみ）、他の風向きに関しては加重平均にて発電量予測値を求めることができる。したがって、風向補正テーブルのデータ量が少なくても、様々な風向（例えば１６方位）に対応して発電量予測値を求めることができる。このように、予め記憶しておくデータのデータ量が少なくて済み、エンジニアリング・コストを削減できる。
【００８７】
また、風向補正テーブルの補正値を用いた補正を行うことで、一般的な気象予報情報（風向・風速予報値）を用いるにも係らず、予測対象風車が実際に受けると予測される風速に応じた精度の高い発電量予測値を得ることができる。
【００８８】
更に、風上側風車の存在と特徴（運転時と停止時とで、風速等に与える影響が異なる）を考慮して、２種類の風速補正テーブルを設定しておき適宜使用することで、より的確な精度の高い発電量予測値を得ることができる。
【００８９】
次に、以下、実施例２について説明する。
実施例２では、ファジィ推論を用いて発電量を求める。
実施例２に関しても、図１に示す構成図を用いて説明する。
【００９０】
実施例２に関する図１の各処理機能部は、まず、気象予報入力機能部１１、風力発電設備運転情報取得部１２、及びパワーカーブ記憶部１３は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【００９１】
補正テーブル記憶部１６に記憶される補正テーブル、及び発電量計算機能部１４の処理内容が、実施例１とは異なる。上記の通り、実施例２ではファジィ推論を用いて発電量を求める。以下、実施例２における補正テーブル記憶部１６、及び発電量計算機能部１４について説明する。
【００９２】
まず、図４（ａ）、（ｂ）は、実施例２において補正テーブル記憶部１６に記憶されている補正テーブルの一例を示す図である。
図４（ａ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
【００９３】
図４（ｂ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
実施例２においても、実施例１と同様に、風力発電設備運転情報取得部１２が取得する風上側風車の運転スケジュール情報（ＯＮ／ＯＦＦ等）に従って、図４（ａ）、（ｂ）の２つの風速補正テーブルの何れか一方または両方を用いることになる。
【００９４】
ここで、まず、風向きに関しては、例えば、北、北北東、北東等の１６方向で区分するものがよく知られているが、気象庁等から取得できる、風向きに関する特定のデータ（東西成分と南北成分の２種類から成る）に基づいて、３６０°で示される風向に換算することが、既存技術として行われている。尚、ここでは、北風を基点（０°）として時計まわりで示すものとする。よって、東が９０°、南が１８０°、西が２７０°となる。
【００９５】
上記実施例１に係る風速補正テーブル（図３（ｂ）、（ｃ）等）においては、所定の風向（東西南北）毎に、それぞれ、各風速予報値（ｍ）（例えば０（ｍ）、１（ｍ）、２（ｍ）、・・・、１０（ｍ）、・・・のような１ｍ単位）毎に対応する補正値が格納されていた。
【００９６】
実施例２に係る補正テーブル（図４（ａ）、（ｂ）等）においては、上記所定の風向として、上記東西南北の代わりに、上記３６０°表記に係わる特定の角度範囲（風向範囲と呼ぶ場合もあるものとする）を用いている。例えば図４（ａ）、（ｂ）等に示す例では、“２１２°±１０°”の角度範囲と“それ以外”の角度範囲とを示している。尚、“それ以外”とは「“２１２°±１０°”の角度範囲以外の全て」を意味している。
【００９７】
上記の例では、“２１２°±１０°”の角度範囲の方向に、風上側風車が存在しているものとする。この為、“２１２°±１０°”に係る補正値は、基本的には図４（ａ）と図４（ｂ）とで異なる場合があることになる。
【００９８】
実施例２における発電量計算機能部１４も、基本的には実施例１と略同様にして、上記風力発電設備運転情報取得部１２を介して得られる風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、パワーカーブと、風速予報値を補正するための補正用データとに基づいて、風力発電量予測値を求める。但し、実施例２では、風速予報値を補正するための補正用データが、実施例１のものとは異なり、上記図４（ａ）、（ｂ）等に示す風速補正テーブルと、後述する図５に示すメンバシップ関数等となる。そして、これに応じて、当然、風力発電量予測値の算出処理内容も、実施例１とは異なることになる。
【００９９】
実施例２における発電量計算機能部１４は、まず、気象予報入力機能部１１を介して得られる風速予報値に基づいて、図４（ａ）、（ｂ）等に一例を示す風速補正テーブルを参照することで、各風向範囲毎に、風速予報値に対応する補正値を取得すると共にこの補正値によって風速予報値を補正して成る風速補正値を求める。
【０１００】
尚、風速補正テーブルには、補正値ではなく、風速補正値が登録されていてもよく、この場合にはテーブルから該当する風速補正値を取得すればよい。例えば図４（ａ）の例において“２１２°±１０°”に係る風速予報１０（ｍ）に対応する補正値は‘−１’であるので、風速補正値＝１０−１＝９（ｍ）が算出されるが、補正値‘−１’の代わりに風速補正値‘９’が登録されていてもよい。
【０１０１】
本例の場合、風上側風車の運転状況がＯＮ／ＯＦＦのどちらであるかに応じて、図４（ａ）、（ｂ）の２つの風速補正テーブルの何れか一方を選択・使用することになるが、どちらを選択・使用する場合であっても、風速予報値に応じた上記風速補正値として、風向予報に関係なく、各風向範囲毎に風速補正値を求める。よって、図示の例では、上記“２１２°±１０°”に応じた風速補正値ＰＣ２と、“それ以外”に応じた風速補正値ＰＣ１とが求められることになる。尚、風向予報は、後にメンバシップ関数を参照する際に用いられる。
【０１０２】
仮に、風上側風車の運転状況がＯＮであった場合には、図４（ａ）に示す風速補正テーブルが選択・使用されることになり、更に仮に風速予報値が１０（ｍ）であった場合には風速補正値ＰＣ２＝９（ｍ）（＝風速予報値−該当する補正値＝１０−１）、風速補正値ＰＣ１＝１０（ｍ）（＝風速予報値−該当する補正値＝１０−０）が得られることになる。
【０１０３】
発電量計算機能部１４は、更に、例えば図５（ａ）に示すメンバシップ関数を、予め保持している。本例のメンバシップ関数は、上記各風向範囲それぞれに応じた“風速補正の加重比”より成る。図示の例では、上記“２１２°±１０°”に応じた“風速補正の加重比ｐｃ２”と、上記“それ以外”に応じた“風速補正の加重比ｐｃ１”とから成る。尚、これらｐｃ１、ｐｃ２等は、加重比の関数（０°〜３６０°の風向に応じた、加重比の関数）等と呼んでもよい。
【０１０４】
図示のように、これら加重比ｐｃ１、ｐｃ２は、０°〜３６０°の風向に応じて、０．０〜１．０の間の値をとる。基本的には、全ての風向に関して「ｐｃ１＋ｐｃ２」＝‘１’となるように、加重比ｐｃ１、ｐｃ２が規定されている。
【０１０５】
図示の例では、加重比ｐｃ２は、風向２１２°の場合に‘１．０’となり、風向０°〜１９０°および２３４°〜３６０°の範囲内では‘０．０’となり、これら以外の風向では、０．０〜１．０の間の値（但し、０．０と１．０は含まれない）となる。一方、加重比ｐｃ１は、風向２１２°の場合に‘０．０’となり、風向０°〜１９０°および２３４°〜３６０°の範囲内では‘１．０’となり、これら以外の風向では、０．０〜１．０の間の値（但し、０．０と１．０は含まれない）となる。
【０１０６】
つまり、例えば、風向が２１２°の場合は、風速補正値ＰＣ２に対する加重比ｐｃ２を１．０とし、風速補正値ＰＣ１に対する加重比ｐｃ１を０．０とする。一方、風向が１９０°以下や、２３４°以上の場合は、風速補正値ＰＣ１に対する加重比ｐｃ１を１．０とし、風速補正値ＰＣ２に対する加重比ｐｃ２を０．０とする。これらは上記「ｐｃ１＋ｐｃ２」＝‘１’となることは明らかである。
【０１０７】
また、風向が１９０°〜２３４°の範囲内（但し、２１２°は除外する）の場合は、加重比ｐｃ１、ｐｃ２は、それぞれ例えば図示のような０より大きく１未満の範囲内の該当する値（風向に対応する値）となる。これらについては、特に具体的な数値等は示さないが、上記の通り、基本的には「ｐｃ１＋ｐｃ２」＝‘１’の条件を満たすように規定されている。
【０１０８】
尚、２１２°は上記の通り風上側風車が設けられている方向であるが、１０°や１９０°や２３４°は、例えば開発者等が適宜決定した値である。
発電量計算機能部１４は、気象予報入力機能部１１を介して得られる上記風向予報値に基づいて、図５（ａ）に示すメンバシップ関数を参照することで、この風向予報値に対応する加重比ｐｃ１、加重比ｐｃ２を求める。例えば、風向きが２１２°であった場合には、加重比ｐｃ１＝‘０．０’、加重比ｐｃ２＝‘１．０’となる。
【０１０９】
以上、風速予報に応じた風速補正値ＰＣ１、ＰＣ２と、風向予報に応じた、これら風速補正値ＰＣ１、ＰＣ２に対する加重比ｐｃ１、ｐｃ２とを取得したら、これらに基づいて下記の（３）式によって最終的な風速補正値ＰＣを求める。
【０１１０】
風速補正値ＰＣ＝（風速補正値PC1 ×加重比pc1 ＋風速補正値PC2 ×加重比pc2 ）／（加重比pc1＋加重比pc2）・・・（３）式
その後は、例えば図３（ａ）等に示すパワーカーブを用いて、上記（３）式で求めた風速補正値ＰＣに応じた発電量を求める。
【０１１１】
なお、以上の実施例２の説明では、説明を簡単にするため、風上側風車の稼働率などを考慮していないが、実施例１と略同様に、風上側風車の稼働率などを考慮してもよい。仮に、風上側風車が９０％稼働しているとした場合は、発電量予測値は以下のようになる。尚、以下の式は、実質的には上記（４）式と略同様であり、２つの発電量を、風上側風車の稼働率に応じた加重比（９０：１０の加重比）で平均化するものである(これより、(４)’式と記す)。
【０１１２】
発電量予測値＝発電量１×０．９＋発電量２×０．１ ・・・（４）’式
｛ここで、
発電量１：風上側風車が稼働しているときの発電量（式（３）と図４（ａ）の補正テーブルから計算した風速補正値を、図３（ａ）のパワーカーブで発電量に変換して求めたもの）。
【０１１３】
発電量２：風上側風車が停止しているときの発電量（式（３）と図４（ｂ）の補正テーブルから計算した風速補正値を、図３（ａ）のパワーカーブで発電量に変換して求めたもの）。｝
また、図４（ａ），（ｂ）の例では、簡単のため、風速補正に係る風向区分を上記“２１２°±１０°”と“それ以外”の２種類だけとしたが、勿論、この様な例に限らず、３種類以上であっても構わない。例えば図３では西方向に何らかの障害物があるものとしており、この例に準じて、“２７０°±１０°”等もあってもよい。この様な他の例を図４（ｃ）、（ｄ）に示す。
【０１１４】
図４（ｃ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
図４（ｄ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
【０１１５】
図示の通り、これら他の例の風速補正テーブルでは、“２１２°±１０°”と“それ以外”と“２７０°±１０°”の３種類の風向区分（風向範囲）毎に、それぞれ、各風速予報値に応じた補正値が格納されている。勿論、補正値の代わりに上記風速補正値が格納されていてもよい。尚、“２７０°±１０°”に関しては、風上側風車の運転状態は関係ないので、図示の通り図４（ｃ）と図４（ｄ）とで同じ内容となっている。
【０１１６】
尚、本例の場合、“それ以外”は、「“２１２°±１０°”及び“２７０°±１０°”」を除く全ての方向を意味する。
また、この場合には、当然、メンバシップ関数も、３種類の風向区分（風向範囲）に応じたものとなる。図５（ｂ）に、上記図４（ｃ）、（ｄ）の例に応じたメンバシップ関数の一例を示す。図５（ｂ）に示すメンバシップ関数は、上記“２１２°±１０°”に応じた“風速補正の加重比ｐｃ２”と、上記“２７０°±１０°”に応じた“風速補正の加重比ｐｃ３”と、上記“それ以外”に応じた“風速補正の加重比ｐｃ１”とから成る。
【０１１７】
この場合も、上記図５（ａ）に示すメンバシップ関数と略同様に、基本的には、全ての風向に関して「ｐｃ１＋ｐｃ２＋ｐｃ３」＝‘１’となるように、加重比ｐｃ１、ｐｃ２、ｐｃ３が規定されている。
【０１１８】
そして、図４（ａ），（ｂ）及び図５（ａ）の場合と略同様にして、ｍ時間先の風向・風速予報値に基づいて、図４（ｃ）または図４（ｄ）の風速補正テーブルから風速予報値に応じた各風向範囲毎の風速補正値ＰＣ１，ＰＣ２，ＰＣ３を取得すると共に、図５（ｂ）のメンバシップ関数から風向予報値に対応する各風向範囲毎の加重比ｐｃ１、ｐｃ２、ｐｃ３を取得する。
【０１１９】
尚、当然、風上側風車の運転スケジュールに応じて、図４（ｃ）と図４（ｄ）のどちらの風速補正テーブルを使用するのか（あるいは両方とも使用するのか）が決まることになる。
【０１２０】
そして、これら取得したデータに基づいて、上記（３）式と略同様の算出式によって最終的な風速補正値ＰＣを求める。つまり、本例の場合には、下記の（３）’式によって最終的な風速補正値ＰＣを求めることになる。
【０１２１】
風速補正値ＰＣ＝（風速補正値PC1×加重比pc1＋風速補正値PC2×加重比pc2＋風速補正値PC3×加重比pc3）／（加重比pc1＋加重比pc2＋加重比pc3）・・・（３）’式
本実施例では、わずか数種類の風向区分に応じた補正テーブルやメンバシップ関数を予め作成しておくだけで、風向ごとに異なる風速の変化状態を連続的に表現することができる。
【０１２２】
次に、以下、実施例３について説明する。
上記実施例１、２では、風上側風車の存在とその稼働状況を考慮した実施例を示したが、本発明は、この様な例に限らない。例えば、本実施例３のようにしてもよい。
【０１２３】
図６（ａ）は、実施例３におけるパワーカーブの一例であり、これは図３（ａ）に示す実施例１のパワーカーブと、実質的に同じであってよく、ここでは特に説明しない。
図６（ｂ）は、実施例３における風速補正テーブルの一例である。実施例１では、風上側風車の運転／停止に応じた２種類の風速補正テーブルが設けられていたが、実施例３では図示の通り１種類のみである。尚、図示の例では、北方向に風車以外の何らかの障害物があるものとしている。
【０１２４】
実施例１では、上述したように、特定の障害物である風上側風車の運転／停止（ＯＮ／ＯＦＦ）に応じて、２種類の風速補正テーブルを設定して選択的に使用したが、実施例３ではこの様な選択処理は行われない。また、実施例１では風上側風車の稼働率に応じた補正処理等も行ってよいものとしたが、実施例３ではこの様な処理は行われない。また、当然、実施例３では風力発電設備運転情報取得部１２の処理機能は必要ない。これらの点を除けば、実施例３は、実施例１と略同様と見做してもよく、以下、簡単に説明する。
【０１２５】
すなわち、気象予報入力機能部１１を介してｍ時間先の気象予報（特に風向・風速予報値）が入力されると、発電量計算機能部１４は、これらに基づいて、図６（ｂ）に示す風速補正テーブルを参照することで、該当する補正値を得て、以って風速補正値を求める。気象予報値の一例として、例えば風向予報が北で、風速予報値が１０ｍの場合は、図６（ｂ）の風速補正テーブルより、該当する補正値（北で１０ｍの場合、図示の例では補正値＝‘−１’）を得ることで、上記風速補正値を以下のように求める。尚、風速補正値は、風力発電所における処理対象の風力発電設備（その風車）が実際に受けると考えられる風の風速である。
【０１２６】
風速補正値＝９ｍ（＝風速予報値＋補正値＝１０ｍ−１ｍ）
この場合、風力発電量予測値は、図６（ａ）のパワーテーブルより、風速９ｍに対応する９５ｋWとなる。
【０１２７】
なお、上記の説明では、簡単のため、加重平均やファジィ推論を用いなかったが、実施例３において実施例１、２のように加重平均やファジィ推論を用いるようにしてもよい。
例えば、風向予報が北北東で風速予報が１０ｍの場合には、図６（ａ）、（ｂ）に示すパワーテーブルと風速補正テーブルを参照することで、風速補正テーブルに登録されている東西南北の４方向のうち上記風向予報（北北東）に関連する方向である北と東に関して、以下の風速補正値及び発電量が得られる。
【０１２８】
北： 風速補正値９ｍ ， 発電量 ９０ｋW
東： 風速補正値１０ｍ ， 発電量１００ｋW
そして、これら２つの発電量を、風向予報に応じた加重比で平均化することで、風向・風速予報と風速補正テーブルに基づく発電量予測値が算出される。本例では、風向予報が北北東であることから、２５：７５の加重比で平均化することになり、従って、
発電量予測値；（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４ ＝ ９２．５ｋW ・・・（１）’’式
となる。
【０１２９】
尚、上述した例では、実施例１の（１）式、（１）’式も、上記実施例３の（１）’’式も、２つの発電量を求めた後に加重平均算出処理を行ったが、この例に限らない。例えば、上記の例では、北と東の風速補正値（９ｍと１０ｍ）に対して、上記風向予報（北北東）に応じて２５：７５の加重比で平均化することで、風向予報に応じた風速補正値を得るようにしてもよい。
【０１３０】
この例の場合、風向予報（北北東）に応じた風速補正値は、９．２５（ｍ）となる。そして、パワーテーブルを参照することで、風向予報（北北東）に応じた風速補正値（９．２５（ｍ））に対応する発電量を求め、これを上記発電量予測値とするようにしてもよい。但し、図３（ａ）や図６（ａ）に示す例のパワーカーブでは、風速が１ｍ刻みであるので、９．２５（ｍ）に対応する発電量は登録されていない。よって、パワーカーブがこの様な例の場合には、上記（１）式、（１）’式や（１）’’式によって求める手法の方が望ましい。
【０１３１】
この様に、風速補正テーブルには少数の所定方向（本例では東西南北の４方向）に関する補正データを登録しておくだけでよく、これら４方向以外の風向きに関しては加重平均することにより求めることができる。つまり、全ての風向に対してデータベースを持つ必要がないため、事前のデータベース設定が容易であり、エンジニアリング・コストを削減できる。
【０１３２】
あるいは、図６（ｂ）に示す補正テーブルの代わりに、図４に示すような、風向を０°〜３６０°の角度で扱って所定の角度範囲毎に対応する補正値を登録した風速補正テーブルを用いるようにしてもよい。勿論、この場合には、更に、図５に示すようなメンバシップ関数も予め作成・記憶して必要がある。
【０１３３】
尚、実施例３では風上側風車は想定しないので、風速補正テーブルは、図４（ａ）、（ｂ）のような２種類ではなく、１種類のみである。また、それ故、当然、風力発電設備運転情報取得部１２の処理機能は必要ないし、風上側風車のＯＮ／ＯＦＦに応じて２種類の風速補正テーブルの何れか一方（または両方）を選択する処理も必要ない。これらの点を除けば、本例の処理は、実施例２の処理と略同様と見做しても構わない。
【０１３４】
すなわち、実施例２の場合と略同様に、気象予報入力機能部１１を介してｍ時間先の気象予報（風向・風速予報値）を入力すると、風向・風速予報値と不図示の（図４と略同様の）風速補正テーブルに基づいて、上記風速補正値ＰＣ１，ＰＣ２等に相当する風速補正値を求める。更に、風向予報値と不図示の（図５と略同様の）メンバシップ関数を用いて、加重比ｐｃ１、ｐｃ２等に相当する加重比を求める。
【０１３５】
そして、これら求めた風速補正値と加重比とに基づいて、上記（３）式や（３）’式と略同様の算出式によって、最終的な風速補正値ＰＣを求める。その後は、実施例２と略同様に、例えば図３（ａ）等に示すパワーカーブを用いて、求めた風速補正値ＰＣに応じた発電量予測値を求める。
【０１３６】
この様に、実施例３において、ファジィ推論を用いて風速補正値を求め、以って風速発電量予測値を求めるようにしてもよい。
上述した実施例１，２，３に係る図１の風力発電量予測システム１０の各機能部の処理機能は、例えば、以下の様になると言える。
【０１３７】
気象予報入力機能部１１は、外部から風向・風速予報値を取得する。これは、例えば気象庁等による一般的な気象予報情報としての風向・風速予報値を取得するものである。
補正テーブル記憶部１６には、予め、上記風速予報値を補正するための補正用データが記憶される。パワーカーブ記憶部１３には、予め、各風力発電設備毎に、その風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶される。
【０１３８】
発電量計算機能部１４は、上記風向・風速予報値と補正用データとパワーカーブとを用いて、加重平均（実施例１，３）またはファジィ推論（実施例２）によって、風力発電量予測対象の任意の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１３９】
そして、例えば、実施例１，３の場合には、上記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルである。そして、発電量計算機能部１４は、この風速補正テーブルを参照して、風向予報が示す風向に係わる１以上の上記所定の風向毎に、風速予報値に対応する補正データを用いて風速補正値を求めると共に、パワーカーブを用いて、求めた１以上の風速補正値毎に対応する発電量を求める。更に、発電量計算機能部１４は、複数の発電量が求められた場合には、該複数の発電量を風向予報に応じた加重比で平均化することで、風力発電量予測対象の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１４０】
尚、各風向に応じた加重比は、例えば予め人間等が決めて設定・記憶させておけばよい。上述した一例では、複数の所定の風向は東西南北の４方向であり、風向予報がこれら４方向の何れかである場合以外は、複数の発電量が求められることになり、上記加重平均による風力発電量予測値の算出が行われることになる。
【０１４１】
尚、上記風速補正テーブルには、補正データの代わりに風速補正値が格納されていてもよい。この場合、発電量計算機能部１４は、上述した風速予報値に対応する補正データを用いて風速補正値を求める処理は行う必要なく、その代わりに、風速予報値に対応する風速補正値を風速補正テーブルから取得すればよい。
【０１４２】
また、例えば実施例２の場合には、上記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、上記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数である。そして、発電量計算機能部１４は、風速補正テーブルを参照して、各風向範囲毎に、風速予報値に応じた補正データを用いて風向範囲毎風速補正値を求めると共に、メンバシップ関数を参照して、各風向範囲毎に、風向予報値に応じた加重比を求める。更に、求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、風速予報値を補正した風速補正値を算出する。そして、該算出した風速補正値と上記パワーカーブとによって、風力発電量予測対象の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１４３】
尚、上記所定の各風向範囲とは、例えば図４に示す“２１２°±１０°”、“２７０°±１０°”や“それ以外”等であり、この例では風向は３６０°で表されるものとなっている。
【０１４４】
尚、上記風速補正テーブルには、補正データの代わりに風向範囲毎風速補正値が格納されていてもよい。この場合、発電量計算機能部１４は、上述した風速予報値に対応する補正データを用いて風向範囲毎風速補正値を求める処理は行う必要なく、その代わりに、風速予報値に対応する風向範囲毎風速補正値を風速補正テーブルから取得すればよい。
【０１４５】
更に、風上側風車を考慮する場合には、風力発電設備運転情報取得部１２によって、各風力発電設備（少なくとも風上側風車）の運転予定情報を取得する。尚、既に述べたように、風上側風車とは、上記風力発電量予測対象の風力発電設備が受ける風（特にその風速）に影響を及ぼし得る他の風力発電設備（その風車）を意味する。
【０１４６】
この場合、上記風上側風車（風上側風力発電設備）が存在する場合、上記風速補正テーブルは、当該風上側風力発電設備の運転時、停止時それぞれに応じた２種類の風速補正テーブルとなる。
【０１４７】
そして、発電量計算機能部１４は、風上側風力発電設備の上記運転予定情報に基づいて、例えば、該２種類の風速補正テーブルの何れか一方を用いて、上記風速補正値または上記風向範囲毎風速補正値を求める。
【０１４８】
あるいは、発電量計算機能部１４は、上記２種類の風速補正テーブルの両方をそれぞれを用いることで、上記風力発電量予測値を２つ求め、該求めた２つの風力発電量予測値を、上記風上側風力発電設備の稼働率に応じた加重比で平均化することで、最終的な風力発電量予測値とする。
【０１４９】
続いて、以下、実施例４、実施例５について説明する。
まず、実施例４について説明する。
上述した実施例１，３では、パワーカーブと風速補正テーブルを予め作成・記憶しておき、発電量予測処理の際には、まず風速補正テーブルを参照して該当する風速補正値を求めて、次に、パワーカーブを参照することでこの風速補正値に応じた発電量を求めた。
【０１５０】
これに対して、実施例４では、後述するパワーカーブ補正テーブルを予め作成・記憶しておき、発電量予測処理の際には、パワーカーブ補正テーブルを参照することで、一般的な風向・風速予報値に基づき、各風力発電設備毎に的確な発電量予測値を求める。
【０１５１】
実施例４に関しても、図１に示す構成図を用いて説明する。
実施例４に関する図１の各処理機能部は、まず、気象予報入力機能部１１及び風力発電設備運転情報取得部１２は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【０１５２】
パワーカーブ記憶部１３と補正テーブル記憶部１６は削除され、これらの代わりに不図示のパワーカーブ補正テーブル記憶部が設けられる。このパワーカーブ補正テーブル記憶部には、上記パワーカーブ補正テーブル（後に図７（ａ）、（ｂ）に具体例を示す）が、予め記憶されている。
【０１５３】
当然、発電量計算機能部１４の処理内容は、実施例１とは異なる。発電量計算機能部１４は、気象予報入力機能部１１及び風力発電設備運転情報取得部１２によって取得される情報と、上記パワーカーブ補正テーブルとに基づいて、風速発電量予測値を求める。以下、具体例を挙げて説明する。
【０１５４】
図７（ａ）、（ｂ）は、上記パワーカーブ補正テーブルの具体例である。
図７（ａ）は風上側風車が運転（ＯＮ）している場合に対応するパワーカーブ補正テーブルである。
【０１５５】
図７（ｂ）は風上側風車が停止（ＯＦＦ）している場合に対応するパワーカーブ補正テーブルである。
パワーカーブ補正テーブルは、例えば、予め、上記「発電所のある代表地点、もしくは近辺の風速・風向の予報値または現在値」を取得すると共に、各風力発電設備の設置場所で風速等を実測する等し、更にパワーカーブを風車メーカから取得する（あるいは工場等において、風速に応じた発電量を実測する）等して、これらに基づいて人間が作成するものである。
【０１５６】
図７（ａ）、（ｂ）に示すように、パワーカーブ補正テーブルは、１ｍ単位の各風速（ｍ）毎に対応して、設計値（ｋｗ）と、少数の所定の風向（本例では東西南北の４方向）毎の発電量補正値（ｋｗ）（北（補正値）、東（補正値）、南（補正値）、西（補正値））が格納されている。尚、図示の「風速」は、図３等における「風速予報」と異なるようなものではなく、例えば気象予報入力機能部１１が取得する風速予報値が２（ｍ）で風向予報が‘北’である場合には、パワーカーブ補正テーブルにおける「風速」が２（ｍ）の北（補正値）である‘−１’が取得されることになる。よって、図示の「風速」は、「風速予報」に置き換えても良い。
【０１５７】
また、パワーカーブ補正テーブルにおける発電量（設計値）は、風力発電設備のメーカが工場等で計測した値である。つまり、図示の例では、風車が実際に受ける風速が２ｍであったならば、発電量は１０（ｋＷ）になるはずである。
【０１５８】
また、例えば上記北（補正値）とは、風向が北の場合に使用する、発電量に係る補正値であり、例えば図７（ａ）に示す例においては、風速２ｍのときには、本来は１０ｋW発電するが９ｋW（＝１０−１）になることを示している。同様に、風速１０ｍの場合には本来は１００ｋW発電するはずが９０ｋｗ（＝１００−１０）になることを示している。これは、他の風向の補正値（東（補正値）、南（補正値）、西（補正値））についても略同様であり、ここでは特に説明しない。尚、上記“本来は”とは、他の風車や風車以外の障害物の影響で、風速が変化しなければという意味である。
【０１５９】
ここで、図７（ａ）、（ｂ）のパワーカーブ補正テーブルは、上記実施例１の図３（ｂ）、（ｃ）に示す例の風速補正テーブルと略同様に、北に風上側風車、西に山などの障害物があることを想定している。そのため、風上側風車の運転状態（ＯＮかＯＦＦか）によってテーブルの内容（数値）が変化するのは、図示の通り北のみである（図示の通り、東、南、西に関しては、その数値は、図７（ａ）と図７（ｂ）とで同じである）。
【０１６０】
実施例４における発電量計算機能部１４は、上記風力発電設備運転情報取得部１２を介して得られる各風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、上記パワーカーブ補正テーブルとに基づいて、風力発電量予測値を求める。
【０１６１】
例えば、上記風向・風速予報値と各風力発電設備の運転スケジュールに基づいて、図７（ａ）、（ｂ）のパワーカーブ補正テーブルの少なくとも一方は用いて、風力発電量予測値を求める。
【０１６２】
たとえば、風向予報が北東で、風速予報が１０ｍ、かつ予測時間帯における風上側風車が運転予定が“運転（ＯＮ）”の場合は、図７（ａ）に示す（風上側風車が運転している場合の）パワーカーブ補正テーブルを参照することで、各所定風向（東西南北）毎の発電量補正値（ｋｗ）は、下記の通りとなる（尚、ここでは東西南北全てを示すが、風向予報に関係するものだけで構わない。上記一例のように風向予報が北東の例であれば、北と東だけでも構わない）。
【０１６３】
北の場合： ９０ｋW
東の場合：１００ｋW
南の場合：１００ｋW
西の場合：１００ｋW
尚、これらは、図７（ａ）、（ｂ）のパワーカーブ補正テーブルに予め登録されている設計値と補正値とに基づいて算出するものである。但し、この例に限らず、例えば、パワーカーブ補正テーブルには、設計値も補正値も登録されず、代わりに上記発電量補正値（ｋｗ）そのものが予め登録されていてもよい。
【０１６４】
風向予報が東西南北の何れかであれば、その風向きに応じた上記発電量補正値（ｋｗ）をそのまま演算結果（風力発電量予測値）として用いればよい。
風向予報が東西南北以外であれば、上記東西南北の発電量補正値（ｋｗ）のうち風向予報に関係する方向の２つの発電量補正値（ｋｗ）を用いて、これら２つの発電量補正値を風向予報に応じた加重比を用いて加重平均することで、風力発電量予測値を求める。例えば、風向予報が北北東であれば、上記北の場合の発電量補正値（＝９０ｋW）と上記東の場合の発電量補正値（＝１００ｋW）の２つの発電量補正値を用いることになる。
【０１６５】
そして、実施例１等と略同様の加重平均を用いる手法により、これら２つの発電量補正値から風力発電量予測値を求める。
例えば風向予報が北東である場合には、北東は北と東の丁度中間であることから、５０：５０の加重比で平均化することになり、下記の（６）式で風力発電量予測値を求めることになる。
【０１６６】
風力発電量予測値＝（１００ｋW×５０＋９０ｋW×５０）／１００
＝（１００ｋW＋９０ｋW）／２＝９５ｋW ・・・（６）式
となる。
【０１６７】
あるいは、上記の例において、仮に風向予報が北北東であったならば、２５：７５の加重比で平均化することになり、従って、
風力発電量予測値＝（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４＝９２．５ｋW
となる。
【０１６８】
また、仮に予測時間帯における風上側風車の運転予定が停止（ＯＦＦ）であるならば、上記図７（ａ）のテーブルを用いる代わりに図７（ｂ）のテーブルを用いることになる。
また、例えば、風上側の風車の稼働率を考慮する場合には、図７（ａ）、（ｂ）のパワーカーブ補正テーブルそれぞれを用いて、上記のように風力発電量予測値を求めて、これら求めた２つの風力発電量予測値を、風上側風車の稼働率に応じた加重比で平均化することで、最終的な風力発電量予測値を求めることになる。仮に、風上側風車の稼働率が９０％とした場合には、上記２つの風力発電量予測値を、９０：１０の加重比で平均化することで、最終的な風力発電量予測値を求めることになる。
【０１６９】
また、更に、予測対象の風車の稼働率を用いて修正を施して、最終的な予測値とするようにしてもよい。仮に、予測対象の風車の稼働率が８０％の場合は、最終的な風力発電量予測値に対して０．８を乗じることで、最終的な予測値とする。
【０１７０】
尚、上述した処理と予測値は、未来の任意の時間帯（例えば翌日の２時〜３時の１時間分など）に関する発電量予測処理とその結果である。もし、気象予報が１時間ごとｍ時間先まである場合は、上記処理をｍ回行うことで、ｍ時間先まで求めることができる。
【０１７１】
本手法の利点として、風向に関するパワーカーブが少なくても、その間の風速は加重平均にて求めることができるので、精度の高い予測値を得ることができること、および少ないパワーカーブのテーブルだけで予測できるので、エンジニアリング・コストを削減できる点があげられる。
【０１７２】
次に、以下、実施例５について説明する。
本実施例５は、上記実施例２と同様にファジィ推論にて発電量を求めるものであるが、実施例２における図４に示すような風速補正テーブルの代わりに、不図示のパワーカーブ補正テーブルが予め記憶されている。
【０１７３】
不図示のパワーカーブ補正テーブルは、例えば図４のように風上側風車が運転ＯＮの場合とＯＦＦの場合とに応じて２つ用意されると共に、風向きを０°〜３６０°で扱って任意の角度範囲毎（例えば上記“２１２°±１０°等）に対応付けて、補正値が格納されているものである。但し、図４では補正値は風速に係る補正値であったが、パワーカーブ補正テーブルに登録される補正値は、図７と同様に発電量に係る補正値となる。更に、図７と同様に各風速毎の発電量（設計値）も登録される。
【０１７４】
但し、この例に限らず、実施例４と略同様に、補正値と設計値の代わりに、発電量補正値（補正後の発電量）そのものが登録されていてもよい。当然、この場合には、補正値と設計値とに基づいて発電量補正値（補正後の発電量）を算出する処理は必要ない。
【０１７５】
また、例えば図５に示す例のようなメンバシップ関数も、予め記憶されている。
ここで、実施例５に関しても、図１に示す構成図を用いて説明する。
実施例５に関する図１の各処理機能部は、まず、気象予報入力機能部１１及び風力発電設備運転情報取得部１２は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【０１７６】
パワーカーブ記憶部１３と補正テーブル記憶部１６は削除され、これらの代わりに上述した不図示のパワーカーブ補正テーブルが予め記憶されるパワーカーブ補正テーブル記憶部が設けられる。
【０１７７】
パワーカーブ補正テーブルの内容は、例えば以下の２種類である。
通常時（風向２１２°以外の場合）（PC１）： 発電量PC1＝ｆPC１（風速）
風向２１２°の場合（PC２）： 発電量PC2＝ｆPC２（風速）
※ｆPC１、ｆPC２は、特に具体例は示さないが、任意の風速値に応じた発電量を算出するための関数等である。
【０１７８】
※風向は北風が０°で、時計まわり
尚、ここでは、２１２°の方向に、風上側風車が存在する場合を想定している。
発電量計算機能部１４は、まず、上記２つの発電量ＰＣ１、ＰＣ２を求める。これは、気象予報入力機能部１１を介して得られた風速予報値を上記（風速）に代入することで算出される。
【０１７９】
次に、図５（ａ）に示すメンバシップ関数を参照して、気象予報入力機能部１１を介して得られた風向予報に応じた加重比を求める。
つまり、図５（ａ）の例のメンバシップ関数を参照する場合には、予報される風向が２１２°の場合は、上記発電量ＰＣ２に対する加重比ｐｃ２を１．０とし、上記発電量ＰＣ１に対する加重比ｐｃ１を０．０とする。一方、予報される風向が１９０°以下や、２３４°以上の場合は、上記発電量ＰＣ１に対する加重比ｐｃ１を１．０とし、上記発電量ＰＣ２に対する加重比ｐｃ２を０．０とする。風向が１９０°から２３４°の範囲内の任意の値である場合は、図５（ａ）のメンバシップ関数から読み取った加重比ｐｃ１、ｐｃ２を用いる
最終的な発電量予測値ＰＣは次式で求める
発電量予測値ＰＣ＝（発電量ＰＣ１×加重比ｐｃ１＋発電量ＰＣ２×加重比ｐｃ２）／（加重比ｐｃ１＋加重比ｐｃ２）
図８は、本実施例５の効果を示す為の実験結果である。
【０１８０】
図８のグラフは、横軸が時間、縦軸が任意の風力発電設備の発電出力（ｋＷ）である。但し、グラフには風向もあり、この場合、縦軸は右側に示すように０°〜３６０°で表される風向となる。
【０１８１】
図８には、発電出力の実績値を実線で示し、発電出力の予測値を点線または一点鎖線で示す。点線は本実施例５によるファジイ推論を適用した場合であり、一点鎖線はファジイ推論を適用しなかった場合（従来技術）である。
【０１８２】
これより、図８に示すように、ファジィ推論未使用時の発電量は実績値よりも大きいのに対し、ファジィ推論適用時は風向が２１２°付近のときに発電出力が低下し実績値に近い値になっている。
【０１８３】
なお、本実施例５では、説明を簡単にするため単純なファジィ推論により加重比を求めたが、複数のパワーカーブを用いてもよいし、風向だけでなく風上側風車の運転状態も入力条件にしてもよい。
【０１８４】
本実施例５の上記一例では、パワーカーブ補正テーブルにわずか２種類の発電量補正PC1、PC2を用意しただけで、風向ごとに異なる発電量の低下状態を精度よく表現することができる。
【０１８５】
以上説明した本手法の特に実施例１，２，３によれば、以下の効果が得られる。
風力発電予測に関して、風速の予測と、発電量への換算の２段階で行われる。たとえ風速の予測が正確であっても発電量の換算が正確でないと適切な発電量予測値が得られない。その逆に、発電量の換算が正確であっても風速の予測が正確でないと適切な発電量予測値が得られない。
【０１８６】
ここで、例えば日本のような山が多い地形にある風力発電所では、実際に風車が受ける風は、気象予報どおりの風速が得られない場合があった。
これに対して、これまでの風力発電量予測では、特殊な計測装置により実測する必要があったり、風向毎に細かな補正を行うためには多くのテーブルを用意する必要があり、エンジニアリング・コストが掛かった。また、従来の物理モデルを用いる方法は、高精度であるが、膨大な地形情報や計算コストが必要であり、やはり高コストであった。
【０１８７】
本手法では、既設の風力発電所の計測データと一般的な気象予報／気象観測データだけでデータベースを構築でき、特に所定方向（東西南北など）に関する補正情報だけを登録すればよく、更に運用時には一般的な気象予報等に基づいて、データベースを参照した補正を行うことで、各風力発電設備毎に受ける風の風速やそれに応じた発電量を、的確に予測することが可能であり、低コストで高精度の風力発電量予測を可能にする。
【０１８８】
本手法では、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減できる。特に、風向毎に補正するためのデータベースは有するが、全ての風向に対してデータベースを持つ必要がないため、事前のデータベース設定が容易であり、エンジニアリング・コストを削減でき、以って低コストで高精度な予測値を得ることができる。
【０１８９】
本手法を用いることにより、容易・かつ正確に発電量への換算が可能になり、風力発電量予測の精度向上に寄与できる。
更に、特にウィンドファームなどの複数の風車がある場合は、風上側風車の運転状態によって風下側風車が受ける風速が変化するため、風下側風車に関して適切な発電量予測が難しかった。
【０１９０】
これに対して、本手法では、上記特徴・効果に加えて更に、上記データベースとして風上側風車の稼働状態に応じた複数種類のテーブルを設け、風上側風車の運転スケジュールに基づいて該当するテーブルを用いて風速補正を行うことで、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる。
【０１９１】
風上側風車に関しては、それが運転中（ＯＮ）の場合とメンテナンス等で停止中（ＯＦＦ）の場合とで、予測対象の風車が受ける風速に与える影響が異なる。この点を想定して、風上側風車がＯＮの場合とＯＦＦの場合とに応じた２種類の補正テーブルを用意することで、より的確な風力発電量予測を実現できる。
【０１９２】
この様に、本手法は、山や建造物等による影響を考慮して、一般的な気象予報による風速予報値を各風力発電設備毎に補正することで、各風力発電設備（その風車）が実際に受ける風速を的確に予測し、以って発電量を的確に予測することができる。
【０１９３】
また、風速／発電量を補正するためのデータベースは、全ての風向に対する補正用データを登録しなくてもよいので（例えば東西南北の４方向のみ登録）、事前のデータベース設定が容易であり、エンジニアリング・コストを削減することが可能である。
【０１９４】
このように、本発明は、既存の風力発電所設備で得られる観測情報と一般的な気象予報だけで的確な風力発電量予測を実現でき、低コストで高精度な予測値を得ることができる。
【符号の説明】
【０１９５】
１０ 風力発電量予測システム
１１ 気象予報入力機能部
１２ 風力発電設備運転情報取得部
１３ パワーカーブ記憶部
１４ 発電量計算機能部
１５ 予測値格納部
１６ 補正テーブル記憶部
【技術分野】
【０００１】
本発明は、風速予測から風力発電量予測を行う技術に関する。
【背景技術】
【０００２】
風力発電量予測値は、予め作成される、パワーカーブ等と呼ばれる風速と発電量との関係データを用いて、風速予報等による風向・風速予測に基づいて求められる。このパワーカーブは、一般的に風力発電装置の設計値により決定される。
【０００３】
風力発電量に影響を与える主要因は、風速であるため、風力発電予測では実際には風向・風速予測が重要であり、風速等から発電量を求める処理自体は、それほど重要ではない（上記の通り、設計値により決定されたパワーカーブを用いて計算すれば済む）。つまり、風車が受ける風速等を正確に予測できれば、風力発電量予測値が的確に計算できる。
【０００４】
ここで、風速等を予測する従来技術としては、例えば特許文献１，２，３，４に記載の従来技術が知られている。
特許文献１には、風速計および風向計によって風速や風向を観測する装置である観測ポールと、３次元風速計と、予測用コンピュータ等を有する風予測システムが開示されている。特許文献１では、３次元風速計によって参照点の風速および風向を測定するのは、各季節毎に１週間ずつの計４週間で済むようにする手法を提案しており、これによって高価な３次元風速計の利用効率を高め、コストを削減できる。
【０００５】
特許文献２の風力発電出力予測装置は、複数の物理モデルと複数の統計モデルを組み合わせることにより、風力発電出力をより正確に予測できるものである。
特許文献３では、風速予報値を風車高さ相当に補正し、これをパワーカーブを用いて発電量に変換している。
【０００６】
特許文献４では、レーダにより風速を詳細に観測し、ＧＰＶ気象予報を補正するデータベースを持ち、風速を補正している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００５−１００８２号公報
【特許文献２】特開２００７−２３３６３９号公報
【特許文献３】特開２００９−１６７８４８号公報
【特許文献４】特開２００９−１３８５２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
パワーカーブは、実験室や工場などのような理想的な環境で計測され、風力発電設備に関する風速と発電量との関係を正確に反映させたものである。しかし、屋外に設置された実際の風力発電設備では、しばしばパワーカーブどおりの発電量を示さない。その理由は主に以下２つであり、何れの場合も、地形や障害物などの影響で、実際に風車が受ける風速が（たとえ予報自体は的中していても）予報通りにならないことで、発電量予測が外れることになる。
【０００９】
理由１：地形が複雑な場合（近くに建物や山がある場合など）、風向によっては風速が強くなる場合や、弱くなる場合がある。たとえば、いくら風速予報値が大きくても山などの障害物側から吹いている場合は、風車にあたる実際の風速は弱くなるため、風力発電設備の発電量を過大に予測することになる。
【００１０】
理由２：ウィンドファームの場合、風上側風車が稼働していた場合、風下側の風速を弱めるため、風下側風車が受ける実際の風速が弱まる。つまり、風速予測が正確であったとしても風下側風力発電設備の発電量を過大に予測することになる。その一方で、風上側風車がメンテナンス等で停止している場合には、風下側の風速は殆ど弱まることはない。
【００１１】
従来技術として、特許文献１の発明は、観測ポールと３次元風速計とによる実測が必要となる。
特許文献２の発明は、詳細な物理モデルを使用するため、上述理由１、理由２による問題の解決を期待できる有望な方法であるが、詳細な地形情報が必要であり、かつ膨大な計算を必要する欠点がある。
【００１２】
特許文献３の発明は、気象予報で得られる風速を、風車高さ相当に補正するものであり、上記理由１、理由２に係わる問題を解決できるものではない。
特許文献４の発明は、風向ごとにデータベースを持つため、上述理由１、理由２による問題の解決を期待できる方法である。しかし、特殊なレーダによる風速観測が必要である点や、詳細な予測を行う場合には細かな風向ごとにデータベースを持つ必要があり、膨大な観測が必要になる等の欠点がある。つまり、エンジニアリング・コストが膨大になる欠点がある。
【００１３】
本発明の課題は、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減でき、更に、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる風力発電量予測システム等を提供することである。
【課題を解決するための手段】
【００１４】
本発明の風力発電量予測システムは、任意の風力発電設備の発電量を予測する風力発電量予測システムであって、外部から風向・風速予報値を取得する風向・風速予報値取得手段と、予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段とを有する。
【００１５】
上記風力発電量予測システムにおいて、例えば、前記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルであり、前記風力発電量算出手段は、前記風速補正テーブルを参照して、前記風向予報が示す風向に係わる１以上の前記所定の風向毎に、前記風速予報値に対応する前記補正データを用いて風速補正値を求めると共に、前記パワーカーブを用いて、前記求めた１以上の風速補正値毎に対応する発電量を求め、複数の発電量が求められた場合には、該複数の発電量を前記風向予報に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求める。
【００１６】
あるいは、上記風力発電量予測システムにおいて、例えば、前記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、前記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数であり、前記風力発電量算出手段は、前記風速補正テーブルを参照して、前記各風向範囲毎に、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求めると共に、前記メンバシップ関数を参照して、前記各風向範囲毎に、前記風向予報値に応じた前記加重比を求め、前記求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、前記風速予報値を補正した風速補正値を算出し、該算出した風速補正値と前記パワーカーブとによって、前記風力発電設備が受ける風に応じた風力発電量予測値を求める。
【００１７】
上記風力発電量予測システムでは、一般的な風向・風速予報値等と補正用データ等に基づいて、各風力発電設備がそれぞれ実際に受ける風速に応じた風力発電量予測値を求めることが可能となり、地形や障害物等による影響に対応した的確に風力発電量予測を行うことができる。また、加重平均またはファジィ推論を用いることで、補正用データのデータ量を削減できる。例えば全ての風向ではなく上記複数の所定の風向に関してのみ補正データを登録しておけば済むようになる。
【００１８】
また、更に、前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備による影響も考慮して、風力発電量予測値を求めることで、更に的確に風力発電量予測値を求めることができるようにしてもよい。
【発明の効果】
【００１９】
本発明の風力発電量予測システムによれば、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減できる。更に、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる。
【図面の簡単な説明】
【００２０】
【図１】本例の風力発電量予測システムの構成図である。
【図２】パワーカーブの一例を示す図である。
【図３】（ａ）〜（ｃ）は実施例１に係るパワーカーブや風速補正テーブルの一例である。
【図４】（ａ）〜（ｄ）は実施例２に係る風速補正テーブルの一例である。
【図５】（ａ），（ｂ）は、実施例２に係るメンバシップ関数の一例である。
【図６】（ａ），（ｂ）は、実施例３に係るパワーカーブ／風速補正テーブルの一例である。
【図７】（ａ）、（ｂ）は実施例４に係るパワーカーブや発電量補正テーブルの一例である。
【図８】実施例５の効果を示す図である。
【発明を実施するための形態】
【００２１】
以下、図面を参照して本発明の実施の形態について説明する。
図１は、本例の風力発電量予測システムの構成図である。
図示の例の風力発電量予測システム１０は、気象予報入力機能部１１、風力発電設備運転情報取得部１２、パワーカーブ記憶部１３、発電量計算機能部１４、予測値格納部１５、補正テーブル記憶部１６の各種機能部を有する。
【００２２】
尚、風力発電量予測システム１０は、例えば一般的な汎用コンピュータ（パソコン、サーバ装置等）で実現される。つまり、風力発電量予測システム１０は、ハードウェア的には一般的な汎用コンピュータの構成であってよく、特に図示はしないが、例えばＣＰＵ等の演算プロセッサ、メモリやハードディスク等の記憶装置、ネットワークを介して外部のコンピュータ装置と通信可能とする通信部、液晶ディスプレイ等の表示部、キーボードやマウス等の入力部などを有している。
【００２３】
また、上記不図示の記憶装置には予め所定のアプリケーションプログラムが記憶されている。上記不図示の演算プロセッサがこのアプリケーションプログラムを読出し・実行することにより、上記各種機能部１１〜１６の処理機能が実現される。
【００２４】
気象予報入力機能部１１は、外部で作成された気象予報情報（特に、風向・風速予測値など）を、この気象予報情報を記憶している外部のコンピュータ装置等からネットワークを介して入力する入力機能部である（あるいは、上記気象予報情報を見たユーザ等が、上記キーボード等を操作して気象予報情報を入力するものであってもよい）。
【００２５】
入力する気象予報情報は、一般的な気象予報データであり、例えば気象庁やその他の気象業者の発表した気象予報情報そのものであってもよいし、何らかの別の方法で作成された気象予報情報でもよい。何れにしても、一般的な気象予報データは、各風力発電設備（その風車）が実際に受ける風の風速予測値等を示すようなものではなく、ある程度広いエリア毎の風向・風速予測値を示すものである。気象予報入力機能部１１は、例えば風力発電設備が設置されているエリアに関する風向・風速予測値を入力する。
【００２６】
尚、風向は、例えば南、南南西、南西等の１６方位で表すものであってもよいし、真北を基準（０°）にして、東が９０°、南が１８０°、西が２７０°等のように時計周りに表す３６０方位であってもよい。
【００２７】
風力発電設備運転情報取得部１２は、風力発電設備の運転スケジュール情報を取得する機能部である。この運転スケジュール情報は、各風力発電設備の現在や未来の運転状態（運転（ＯＮ）／停止（ＯＦＦ））を示す運転計画データ等である。尚、ここでは、運転停止（ＯＦＦ）は、メンテナンス等に伴って起こるものであり、故障等による運転停止は関係ない。
【００２８】
上記運転スケジュール情報は、本システム１０内に記憶されていてよいし、キーボード等で直接的に入力されてもよいし、不図示の外部装置等に記憶されていてもよい。外部装置に記憶されている場合には、当然、風力発電設備運転情報取得部１２は、不図示のネットワーク等を介して、外部装置から上記運転スケジュール情報等を取得することになる。
【００２９】
パワーカーブ記憶部１３は、予めパワーカーブを記憶しているデータベースである。パワーカーブは、各風力発電設備について風速と発電量の関係を示すデータであり、上記のように従来から存在するものであり、予め任意に作成されて記憶されている。
【００３０】
パワーカーブは、その形式は、テーブル形式であってもよいし、ｎ次関数の形で持っていてもよい。また、１つのパワーカーブであってもよいし、風向毎にパワーカーブをもっていてもよい。風向計を備えており風向きに応じてその向きを変えるタイプの風力発電設備であれば、風向きに関係なく風が吹いてくる方向を向くので（常にほぼ正面から風を受けるので）パワーカーブは１つであってよい。本例では、パワーカーブ自体は１つであるものとする。
【００３１】
図２に、パワーカーブの一例を示す。図示のパワーカーブは、任意の風力発電設備に係る特性データであって、横軸が風速、縦軸が発電量となっており、風速と発電量との関係が示されている。尚、パワーカーブ自体は既存のものであり、これ以上詳細には説明しないものとする。
【００３２】
発電量計算機能部１４は、本システム１０の主要機能であり、その処理内容については後に詳しく説明するが、基本的には上記各機能部１１，１２から得られる入力情報等や上記パワーカーブや後述する補正テーブル等に基づいて、現在から（あるいは未来ｎ時間先から）未来ｍ時間先までの風力発電設備の発電量の予測値を算出・出力する。これは、例えば、所定時間単位の（例えば１時間単位の）各予測対象時間帯毎に、風力発電設備の発電量の予測値を算出・出力する。
【００３３】
風力発電設備が複数ある場合には、各風力発電設備毎にそれぞれ、現在から未来ｍ時間先までの発電量の予測値を算出・出力する。
発電量計算機能部１４によって算出・出力される上記「ｍ時間先までの風力発電設備の発電量の予測値」は、当該予測値を格納するデータベースである予測値格納部１５に格納される。
【００３４】
補正テーブル記憶部１６には、予め、風速予報値を補正するための補正テーブルが記憶されている。これについては、後に詳しく説明する。
上記風力発電設備の運転スケジュール情報は、発電所において複数の風車（風力発電設備）がある場合に、特に必要になるものである。詳しくは後述する。
【００３５】
以下、上記風力発電量予測システム１０の処理機能について、各種実施例を提案・説明する。尚、以下に説明する実施例１〜５における風力発電量予測システム１０運用時の各種処理は、風向・風速予報値や運転スケジュール等を取得する処理以外は、基本的には全て発電量計算機能部１４が実行するものと考えてよい。
【００３６】
まず、実施例１について説明する。
気象予報入力機能部１１からは、例えば１時間先〜ｍ時間先までの１時間毎の、発電所のある代表地点、もしくは近辺の風向・風速予報値が得られる。換言すれば、上記風力発電設備が設置されているエリアの風向・風速予報値が得られるものである。
【００３７】
風力発電設備運転情報取得部１２からは、各風力発電設備の運転スケジュールが得られる。つまり、風力発電設備に関しては、定期点検などのメンテナンス等が理由で運転を停止する期間（よって、事前に決めておくことが可能）があり、この様な運転停止期間（ＯＦＦ期間）と、通常の運転期間（ＯＮ期間）とを示す運転スケジュールが、予め作成されて記憶されているので、これを取得する。これは、勿論、現在の運転状態（ＯＮ／ＯＦＦ）だけでなく、未来の運転予定（ＯＮ／ＯＦＦ）も得られる。例えば、１時間先〜ｍ時間先までの運転スケジュールを得る。尚、運転がＯＦＦ（停止）の期間は、発電量は‘０’である。
【００３８】
基本的には、発電所内にある全ての風力発電設備について、それぞれ、上記運転スケジュールが作成されて記憶されているが、風力発電設備運転情報取得部１２が取得するのは、予測処理対象の風力発電設備の運転スケジュールや、その風上側風力発電設備（風上側風車）の運転スケジュールのみであってよい。また、必ずしも運転スケジュールを取得する必要はない。
【００３９】
尚、例えば、各風力発電設備には予め識別用ＩＤ等が付与されており、特に図示しないが、各風力発電設備毎に、その風上側風力発電設備の識別用ＩＤが登録されているものとする。また、各識別用ＩＤに対応付けてその風力発電設備の運転スケジュールが管理されている。これによって、例えばユーザが任意の風力発電設備を予測処理対象に指定すると、その風上側風力発電設備を判別できると共に、これらの運転スケジュールを取得できる。
【００４０】
また、尚、現在時点の運転状態はＯＮ／ＯＦＦの２値の何れかの状態しかないが、未来の運転予定は、例えば仮に１時間中に６分間停止する予定であれば「稼働率が９０％」等のように、アナログ量として取得してもよい（尚、これは、１時間単位で取得する場合であり、また運転スケジュールは分単位で管理されているものとするが、この例に限らない）。
【００４１】
尚、本例では、上記のように入力・取得するデータは、１時間先〜ｍ時間先まで１時間単位で扱うと共に、これに合わせて発電量予測値も１時間先〜ｍ時間先まで１時間単位で扱うものとする。
【００４２】
ここで、上記運転状態におけるＯＦＦ状態は、あくまでもメンテナンス等で発電できない状態であり、風速が‘０’で発電できない状態とは異なる。
尚、上記風上側風力発電設備（風上側風車）は、予測対象の風力発電設備（対象風車）の風上側にある他の風力発電設備であるが、風上側や風上側風車や風上側風力発電設備とは、ある特定の風向きのときに対象風車の風上となることを意味し、風向きに関係なく常に風上側となるようなことを意味するものではない（換言すれば、風上側風車とは、風上側となったときに影響を及ぼす（対象風車が受ける風を弱める）可能性がある、他の風車を意味するものであり、例えば対象風車の近傍に設けられる他の風車等である）。
【００４３】
また、パワーカーブは、その一例を図２に示したが、例えば図３（ａ）に示すようにテーブル形式であってもよい。図示の例では、１ｍ単位の各風速（ｍ）毎に、対応する発電量（ｋｗ）が格納されているが、この例に限らない。尚、既に従来技術で述べたように、パワーカーブの各値は、例えばメーカが工場等で計測した値である。
【００４４】
また、図３（ｂ）、（ｃ）には、補正テーブル記憶部１６に予め記憶されている風速補正テーブルの具体例を示す。
風速補正テーブルは、一般的な気象情報としての風速予報値を、その風車が受ける風速へと補正するためのテーブルである。尚、以下、風速補正テーブルを簡略化して単に「補正テーブル」と記す場合もあるものとする。
【００４５】
補正テーブルは、例えば、予め、例えば上記気象庁等による一般的な気象予報データから、例えば「発電所のある代表地点もしくは近辺の風速・風向の予報値または現在値」を取得すると共に、各風力発電設備の設置場所で風速計等を用いて風速等を実測する等して、これらに基づいて人間が作成するものである。
【００４６】
すなわち、例えば気象庁等が発表する風速・風向の現在値（観測値）を取得すると共に、各風力発電設備に設けられた風速計で実測された風速データを取得して、両者の差分（実測データ−風速の現在値）を求めて、この差分を上記風速・風向の現在値と共に記録する。または、実測データを上記風速・風向の現在値に対応付けて記憶するようにしてもよい。
【００４７】
あるいは、例えば気象庁等が発表するｍ時間先の風向・風速予報値を取得し、その後、ｍ時間後に各風力発電設備に設けられた風速計等で風速等を実測する。そして、風速予報値と実測データとの差分（実測データ−風速予報値）を求めて、この差分を上記風向・風速予報値に対応付けて記憶するようにしてもよい。または、実測データを上記風向・風速予報値に対応付けて記憶するようにしてもよい。
【００４８】
尚、上記風向は、予報値等の代わりに、各風力発電設備に設けられた風向計等で実測してもよい。
以上の作業を何度も繰り返すことで補正テーブルを作成できる情報が収集できたら、これに基づいて補正テーブルを作成する。尚、風向きは制御できないので、本例では少なくとも東西南北の４方向の風向きについてデータ収集できるまで、上記データ収集を続けることになる。そして、例えば、東西南北の各方向毎に、上記差分の平均値を求める。例えば、東を例にすると、上記収集したデータから風向が東であるデータを全て抽出して、これらを更に上記風速予報値毎（または風速の現在値毎）に分類して、各風速予報値毎（または風速の現在値毎）に、上記差分値を合計してデータ数で除算することで、東に関する各風速毎の上記差分の平均値が求められる（これが後述する図３（ａ）等に示す“東（補正値）”として風速補正テーブルに格納されることになる）。
【００４９】
ここで、本説明では、上記風上側風力発電設備（風上側風車）を考慮する処理について示すものとし（よって、上記運転スケジュールを取得する必要がある）、この場合には風速補正テーブルは、図３（ｂ）、（ｃ）に示す２種類の風速補正テーブルを用いることになる。
【００５０】
図３（ｂ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
図３（ｃ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
【００５１】
上記風力発電設備運転情報取得部１２で得られる各風力発電設備の運転状態のうち、上記予測対象の風力発電設備の近傍等にある（風速等に影響を及ぼし得る）他の風力発電設備（上記風上側風力発電設備）の運転状態（予測対象時間帯の運転状態）に応じて、図３（ｂ）、（ｃ）の補正テーブルの何れか一方を用いることになる。当然、停止（ＯＦＦ）状態の場合は図３（ｃ）を用い、運転（ＯＮ）状態の場合には図３（ｂ）を用いることになる。
【００５２】
図３（ｂ）、（ｃ）の風速補正テーブルは、データ内容が異なるだけであり、テーブル形式は同一である。よって、以下、特に両者を区別する必要がない場合には、図３（ｂ）等の風速補正テーブルなどと記すものとし、特に区別することなく説明する。
【００５３】
図３（ｂ）等の風速補正テーブルは、１ｍ単位の各風速予報（ｍ）に対応して、各風向（東西南北）毎の補正値（風速予報に対する補正値）が格納されている。尚、図では、「風速予報」と記してあるが、これは「風速」に置き換えても構わない。
【００５４】
例えば、図示の「北（補正値）」は、風向が北に対応する補正値であり、例えば風向予報が“北”で風速予報が２ｍのときには、「北（補正値）」＝‘−１’を使用することになるので、この補正値を用いた補正後の風速は、１ｍ（＝２−１）になる。同様に、例えば、風向予報が“西”で風速が１０ｍの場合には、「西（補正値）」＝‘−２’を使用することになるので、この補正値を用いた補正後の風速は、８ｍ（＝１０−２）になる。
【００５５】
但し、この様な例に限らず、後に説明するように、上記“補正値”の代わりに上記“補正後の風速”が、風速補正テーブルに格納されていてもよい。この場合、上記“補正後の風速”を算出する処理は必要なくなる。
【００５６】
従来であれば、上記の例では、それぞれ、「風速２ｍ」、「風速１０ｍ」によりパワーカーブで発電量を求めることになり、発電量予測が実際よりも過大となる場合があったが、本手法では、それぞれ、「風速１ｍ」、「風速８ｍ」に補正してから（予測対象の風車が実際に受けるものと予測される風速に補正してから）パワーカーブで発電量を求めることになり、的確な発電量予測が行える。
【００５７】
尚、ここでは、補正テーブルの補正値が‘０’であり従って実質的には補正が行われなかった場合でも、風速補正値として扱うものとする。
尚、風向予報が東西南北の４方向以外（例えば北東、北北東等）である場合には、後述する加重平均によって上記発電量予測を行う。逆に言えば、加重平均を用いた演算を行うようにすることで、風速補正テーブルには東西南北の４方向のみ登録すれば済み、データベースのデータ量を削減できる。詳しくは後述する。
【００５８】
ここで、図３（ｂ）、（ｃ）に示す例の風速補正テーブルは、北に風上側風車、西に山などの障害物があることを想定している。そのため、風上側風車の運転状態（ＯＮかＯＦＦか）によってテーブルの内容（数値）が変化するのは、図示の通り北のみである。一方図示の通り、東、南、西に関しては、その数値は、図３（ｂ）と図３（ｃ）とで同じである。つまり、障害物が山やビル等である場合には、その運転状態等は関係ないので（そもそも山には運転状態というものはない）、図３（ｂ）と図３（ｃ）とで同じとある。
【００５９】
一方、障害物が風力発電設備である場合には、運転停止（ＯＦＦ）状態のときには殆ど影響しない（風が弱まらない）が、運転状態（ＯＮ）のときには影響を与えることになる（風が弱まる）。これに応じて、上記図３（ｂ）、（ｃ）のような２種類の風速補正テーブルが、予め設定されている。勿論、これは、予め上記風速補正テーブル作成の為のデータ収集の際に、風上側風車が運転時と停止時とに分けてデータ収集することで、２種類の風速補正テーブルが作成されることになる。
【００６０】
発電量計算機能部１４では、上記風向・風速予報値、風上側風車の運転状態、風速補正テーブル、パワーカーブ等に基づいて、風力発電量予測値を求める。この処理の詳細は、後に具体例を示しながら説明する。
【００６１】
ここで既に述べたように、発電量計算機能部１４は、風向予報が東西南北の４方向の何れかの方向であれば、該当する上記「補正後の風速値」を用いてパワーカーブから該当する発電量を取得すれば済むが、これら以外の風向きであれば、パワーカーブから得た複数の発電量に基づいて所定の算出処理を行う。
【００６２】
後に示す（１）式、（２）式は、風向が“北東”の場合を例にした算出処理例を示している。基本的には、北東であれば北と東、北北東も北と東、南西であれば南と西というように、東西南北の４方向以外の風向きは、東西南北の４方向のうちの２方向に関係するものとなる。これより、特に図示しないが、各風向毎に（但し東西南北は除く）、その風向に関係する東西南北のうちの２方向が、予め記憶されていてもよい。
【００６３】
これより、予報された風向に関係する２方向に係る発電量をパワーカーブから取得し、これら２つの発電量を用いて、後に示す（１）式、（２）式等の算出式による算出処理を行う。
【００６４】
また、発電量計算機能部１４は、予測対象の風車の稼働率に応じた、または／及び、その風上側風車の稼働率に応じた、発電量予測値の算出処理を行うものであってもよい。
発電量計算機能部１４の処理機能について、以下、具体例を用いて説明する。
【００６５】
発電量計算機能部１４は、上記風力発電設備運転情報取得部１２を介して得られる風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、上記パワーカーブと上記風速補正テーブルとに基づいて、風力発電量予測値を求める。
【００６６】
たとえば、パワーカーブと風速補正テーブルが、図３（ａ）、（ｂ）、（ｃ）に示す内容であり、風向予報が北東で、風速予報が１０ｍの場合は、風上側風車の運転状態（ＯＮ／ＯＦＦ）に応じて、それぞれ、下記のようにして求める。
【００６７】
尚、図３（ａ）のパワーテーブルでは示されていないが、ここでは仮に風速８（ｍ）のときの発電量は８０（ｋＷ）であるものとする。
（ａ）風上側風車が運転している場合
上記具体例の場合、上記図３（ｂ）の風速補正テーブルを用いて得られる風速補正値、この風速補正値に応じて図３（ａ）のパワーテーブルから得られる発電量は、所定の風向（本例では東西南北）それぞれについては、下記の通りなる。
【００６８】
北の場合： 風速補正値９ｍ ， 発電量 ９０ｋW
東の場合： 風速補正値１０ｍ ， 発電量１００ｋW
南の場合： 風速補正値１０ｍ ， 発電量１００ｋW
西の場合： 風速補正値８ｍ ， 発電量８０ｋW
但し、これらは参考の為に示しており、実際には風向予報に応じたもの（上記２方向）だけを取得すればよい。すなわち、本例では風向予報は北東であるので、北の発電量（＝９０kW）、東の発電量（＝１００kW）を取得すればよい。そして、これらを用いて後述する加重平均によって風力発電量予測値を求める。但し、風向予報が東西南北の４方向の何れかである場合には、その方向に応じた発電量をそのまま風力発電量予測値とすればよい。
【００６９】
尚、上記風速補正値とは、上記「補正後の風速値」に相当する。ここで、図３(ｂ)等の例では、補正値として上記“差分”を格納しており、従って、「風速補正値（補正後の風速値）＝風速予報値＋補正値」を算出する。例えば、図３(ｂ)の例では風速予報１０ｍに対応する「北（補正値）」は‘−１’であるので、風速補正値（補正後の風速値）＝風速予報＋補正値＝１０−１＝９となり、上記の通り風向が北に対しては風速補正値９ｍとなっている。他も同様である。
【００７０】
尚、この例に限らず、例えば上記“差分”の代わりに、上記実測データが格納されていても構わない。つまり、この場合には、図３（ｂ）等の風速補正テーブルには、各風速予報値に対応付けて上記“風速補正値”が格納されることになる。つまり、例えば風速予報１０ｍに対応する「北（補正値）」は‘−９’となるので、これをそのまま取得すればよく、上記「風速補正値（補正後の風速値）＝風速予報＋補正値」の演算を行う必要はない。
【００７１】
上記の通り、一例では、北の発電量（＝９０kW）、東の発電量（＝１００kW）を取得することになり、また予報された風向は北東であり、北と東の丁度中間であることから、北の発電量（＝９０kW）と東の発電量（＝１００kW）との和を２で除算する（両者の平均を求める）ことで、予報された風向（北東）と風速（１０（ｍ））に基づく、予測対象の風力発電設備が受けると予測される風速に応じた発電量予測値が得られることになる。上記の具体例では下記の（１）式の通りとなる。
【００７２】
発電量予測値；（１００ｋW＋９０ｋW）／２＝９５ｋW・・・（１）式
尚、上記（１）式は、東西南北のうち風向予報に係る２方向に関する発電量を用いて、この２つの発電量を加重平均して求める処理の一例である。
【００７３】
上記のように、北東は北と東の丁度中間であることから、５０：５０の加重比で平均化することになり、
発電量予測値；（１００ｋW×５０＋９０ｋW×５０）／１００
＝（１００ｋW＋９０ｋW）／２
となって、上記（１）式となる。
【００７４】
これより、例えば、北北東であれば、２５：７５の加重比で平均化することになり、従って、
発電量予測値；（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４＝９２．５ｋW・・・（１）’式
となる。
【００７５】
尚、予め例えば１６方向の各風向毎に、その風向に関する加重比が記憶されていてもよい。
上述したように、例えば図３（ｂ）等に示す風速補正テーブルには、予め任意に決められた複数の所定の風向毎に（本例では東西南北毎に）、各風速毎に対応する補正データが格納されている。この風速補正テーブルを参照して、予報された風向に係わる１または複数の風速補正値を求める。つまり、上記所定の風向（東西南北）のうち予報された風向に係わる全ての所定の風向（“関係する所定風向”と呼ぶものとする）について、風速予報値に応じた風速補正値を求める。これは、例えば図３（ｂ）等に示すように上記“差分”である補正値が格納されている例では、上記“関係する所定風向”それぞれについて風速予報値に応じた補正値を取得し更に風速予報値を補正値によって補正することで上記風速補正値を求める。尚、上記の通り、風速補正値そのものが格納されている場合には、これをそのまま取得すればよく、算出する処理は必要ない。
【００７６】
上記“関係する所定風向”は、“予報された風向”が上記一例（北東や北北東）のように上記複数の所定の風向（本例では東西南北）以外である場合には、基本的には上記所定の風向のうちの２方向となる。これは、“予報された風向”が例えば北東や北北東である場合には（その他、東北東の場合でも）北と東となる。従って、この例では、北と東それぞれについて上記風速補正値が求められ、以ってパワーテーブルから対応する発電量が得られることになる。
【００７７】
この様に、風向・風速予報値に応じて複数の発電量が求められる場合には、これら複数の発電量を、上記“予報された風向”に応じた加重比で平均化することで、発電量予測値が算出される（一例が上記（１）式である）。上記の例では、例えば“予報された風向”が北北東である場合には、北が７５、東が２５の加重比が用いられることになる（上記（１）’式となる）。これは、後述する図３（ｃ）等の場合も同様である。
【００７８】
（ｂ）風上側風車が運転していない場合
この場合も、補正テーブルを図３（ｂ）の代わりに図３（ｃ）を用いる点以外は、上記（ａ）風上側風車が運転している場合と同様であり、以下、簡略化して説明する。
【００７９】
この場合、所定の風向（本例では東西南北）それぞれについて発電量は、下記の通りなる。
北の場合： 風速補正値１０ｍ ， 発電量 １００ｋW
東の場合： 風速補正値１０ｍ ， 発電量１００ｋW
南の場合： 風速補正値１０ｍ ， 発電量１００ｋW
西の場合： 風速補正値８ｍ ， 発電量８０ｋW
従って、この例では、予報される風向が北東である場合には、発電量予測値は下記の（２）式の通りとなる。
【００８０】
発電量予測値；（１００ｋW＋１００ｋW）／２＝１００ｋW・・・（２）式
となる。
風上側風車の運転状態を、ＯＮ／ＯＦＦで扱う場合には、例えば予報される風向が北東である場合には、上記（１）式、（２）式の何れかによって、風力発電量予測値を求めることになる。
【００８１】
但し、この例に限らず、例えば風上側風車が、予測対象時間帯のなかで９０％の時間は運転し、残りの１０％の時間は停止している（この例を、稼働率＝９０％とする）場合もあり得る（尚、予測対象時間帯とは、例えば１時間単位で予測する場合には、例えば１〜２時間後、２〜３時間後、・・・、ｍ−１〜ｍ時間後等となる。）
この様な場合、（１）式、（２）式の両方の式による算出を行い、これら２つの算出結果に基づいて、風上側風車の稼働率に応じた加重比で平均化することで、風力発電量予測値を求めることになる。
【００８２】
例えば、上記の例において任意の予測対象時間帯における風上側風車の稼働率が９０％の場合には、上記（１）式、（２）式によって２つの発電量（９５ｋWと１００ｋW）が得られたら、これら２つの発電量を、風上側風車の稼働率に応じた加重比（９０：１０の加重比）で平均化する。勿論、この場合、風上側風車が運転している場合の発電量（＝９０kW）に対しては９０％、風上側風車が運転していない場合の発電量（＝１００kW）に対しては１０％が割当てられる。すなわち、下記の通りとなる。
【００８３】
発電量予測値；（９５ｋW×９０＋１００ｋW×１０）／１００＝９５．５ｋW ・・・（４）式
また、更に、予測対象の風車の稼働率を用いて修正を施して、最終的な予測値とするようにしてもよい。
【００８４】
例えば、仮に、予測対象の風車の稼働率が８０％とした場合には、上述したようにして求めた発電量（（１）式では９５ｋW、（２）式では１００ｋW、あるいは上記９５．５ｋW）に対して、０．８を乗じることで、最終的な発電量予測値とする。
【００８５】
尚、上述した処理と予測値は、未来の任意の時間帯（例えば翌日の２時〜３時の１時間分など）に関する発電量予測処理とその結果である。もし、気象予報が１時間ごとｍ時間先まである場合は、上記処理をｍ回行うことで、ｍ時間先まで求めることができる。
【００８６】
以上説明した実施例１によれば、風速補正テーブルには所定の風向毎の補正値が格納される。この所定の風向の数が少なくても（例えば東西南北の４方位のみ）、他の風向きに関しては加重平均にて発電量予測値を求めることができる。したがって、風向補正テーブルのデータ量が少なくても、様々な風向（例えば１６方位）に対応して発電量予測値を求めることができる。このように、予め記憶しておくデータのデータ量が少なくて済み、エンジニアリング・コストを削減できる。
【００８７】
また、風向補正テーブルの補正値を用いた補正を行うことで、一般的な気象予報情報（風向・風速予報値）を用いるにも係らず、予測対象風車が実際に受けると予測される風速に応じた精度の高い発電量予測値を得ることができる。
【００８８】
更に、風上側風車の存在と特徴（運転時と停止時とで、風速等に与える影響が異なる）を考慮して、２種類の風速補正テーブルを設定しておき適宜使用することで、より的確な精度の高い発電量予測値を得ることができる。
【００８９】
次に、以下、実施例２について説明する。
実施例２では、ファジィ推論を用いて発電量を求める。
実施例２に関しても、図１に示す構成図を用いて説明する。
【００９０】
実施例２に関する図１の各処理機能部は、まず、気象予報入力機能部１１、風力発電設備運転情報取得部１２、及びパワーカーブ記憶部１３は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【００９１】
補正テーブル記憶部１６に記憶される補正テーブル、及び発電量計算機能部１４の処理内容が、実施例１とは異なる。上記の通り、実施例２ではファジィ推論を用いて発電量を求める。以下、実施例２における補正テーブル記憶部１６、及び発電量計算機能部１４について説明する。
【００９２】
まず、図４（ａ）、（ｂ）は、実施例２において補正テーブル記憶部１６に記憶されている補正テーブルの一例を示す図である。
図４（ａ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
【００９３】
図４（ｂ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
実施例２においても、実施例１と同様に、風力発電設備運転情報取得部１２が取得する風上側風車の運転スケジュール情報（ＯＮ／ＯＦＦ等）に従って、図４（ａ）、（ｂ）の２つの風速補正テーブルの何れか一方または両方を用いることになる。
【００９４】
ここで、まず、風向きに関しては、例えば、北、北北東、北東等の１６方向で区分するものがよく知られているが、気象庁等から取得できる、風向きに関する特定のデータ（東西成分と南北成分の２種類から成る）に基づいて、３６０°で示される風向に換算することが、既存技術として行われている。尚、ここでは、北風を基点（０°）として時計まわりで示すものとする。よって、東が９０°、南が１８０°、西が２７０°となる。
【００９５】
上記実施例１に係る風速補正テーブル（図３（ｂ）、（ｃ）等）においては、所定の風向（東西南北）毎に、それぞれ、各風速予報値（ｍ）（例えば０（ｍ）、１（ｍ）、２（ｍ）、・・・、１０（ｍ）、・・・のような１ｍ単位）毎に対応する補正値が格納されていた。
【００９６】
実施例２に係る補正テーブル（図４（ａ）、（ｂ）等）においては、上記所定の風向として、上記東西南北の代わりに、上記３６０°表記に係わる特定の角度範囲（風向範囲と呼ぶ場合もあるものとする）を用いている。例えば図４（ａ）、（ｂ）等に示す例では、“２１２°±１０°”の角度範囲と“それ以外”の角度範囲とを示している。尚、“それ以外”とは「“２１２°±１０°”の角度範囲以外の全て」を意味している。
【００９７】
上記の例では、“２１２°±１０°”の角度範囲の方向に、風上側風車が存在しているものとする。この為、“２１２°±１０°”に係る補正値は、基本的には図４（ａ）と図４（ｂ）とで異なる場合があることになる。
【００９８】
実施例２における発電量計算機能部１４も、基本的には実施例１と略同様にして、上記風力発電設備運転情報取得部１２を介して得られる風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、パワーカーブと、風速予報値を補正するための補正用データとに基づいて、風力発電量予測値を求める。但し、実施例２では、風速予報値を補正するための補正用データが、実施例１のものとは異なり、上記図４（ａ）、（ｂ）等に示す風速補正テーブルと、後述する図５に示すメンバシップ関数等となる。そして、これに応じて、当然、風力発電量予測値の算出処理内容も、実施例１とは異なることになる。
【００９９】
実施例２における発電量計算機能部１４は、まず、気象予報入力機能部１１を介して得られる風速予報値に基づいて、図４（ａ）、（ｂ）等に一例を示す風速補正テーブルを参照することで、各風向範囲毎に、風速予報値に対応する補正値を取得すると共にこの補正値によって風速予報値を補正して成る風速補正値を求める。
【０１００】
尚、風速補正テーブルには、補正値ではなく、風速補正値が登録されていてもよく、この場合にはテーブルから該当する風速補正値を取得すればよい。例えば図４（ａ）の例において“２１２°±１０°”に係る風速予報１０（ｍ）に対応する補正値は‘−１’であるので、風速補正値＝１０−１＝９（ｍ）が算出されるが、補正値‘−１’の代わりに風速補正値‘９’が登録されていてもよい。
【０１０１】
本例の場合、風上側風車の運転状況がＯＮ／ＯＦＦのどちらであるかに応じて、図４（ａ）、（ｂ）の２つの風速補正テーブルの何れか一方を選択・使用することになるが、どちらを選択・使用する場合であっても、風速予報値に応じた上記風速補正値として、風向予報に関係なく、各風向範囲毎に風速補正値を求める。よって、図示の例では、上記“２１２°±１０°”に応じた風速補正値ＰＣ２と、“それ以外”に応じた風速補正値ＰＣ１とが求められることになる。尚、風向予報は、後にメンバシップ関数を参照する際に用いられる。
【０１０２】
仮に、風上側風車の運転状況がＯＮであった場合には、図４（ａ）に示す風速補正テーブルが選択・使用されることになり、更に仮に風速予報値が１０（ｍ）であった場合には風速補正値ＰＣ２＝９（ｍ）（＝風速予報値−該当する補正値＝１０−１）、風速補正値ＰＣ１＝１０（ｍ）（＝風速予報値−該当する補正値＝１０−０）が得られることになる。
【０１０３】
発電量計算機能部１４は、更に、例えば図５（ａ）に示すメンバシップ関数を、予め保持している。本例のメンバシップ関数は、上記各風向範囲それぞれに応じた“風速補正の加重比”より成る。図示の例では、上記“２１２°±１０°”に応じた“風速補正の加重比ｐｃ２”と、上記“それ以外”に応じた“風速補正の加重比ｐｃ１”とから成る。尚、これらｐｃ１、ｐｃ２等は、加重比の関数（０°〜３６０°の風向に応じた、加重比の関数）等と呼んでもよい。
【０１０４】
図示のように、これら加重比ｐｃ１、ｐｃ２は、０°〜３６０°の風向に応じて、０．０〜１．０の間の値をとる。基本的には、全ての風向に関して「ｐｃ１＋ｐｃ２」＝‘１’となるように、加重比ｐｃ１、ｐｃ２が規定されている。
【０１０５】
図示の例では、加重比ｐｃ２は、風向２１２°の場合に‘１．０’となり、風向０°〜１９０°および２３４°〜３６０°の範囲内では‘０．０’となり、これら以外の風向では、０．０〜１．０の間の値（但し、０．０と１．０は含まれない）となる。一方、加重比ｐｃ１は、風向２１２°の場合に‘０．０’となり、風向０°〜１９０°および２３４°〜３６０°の範囲内では‘１．０’となり、これら以外の風向では、０．０〜１．０の間の値（但し、０．０と１．０は含まれない）となる。
【０１０６】
つまり、例えば、風向が２１２°の場合は、風速補正値ＰＣ２に対する加重比ｐｃ２を１．０とし、風速補正値ＰＣ１に対する加重比ｐｃ１を０．０とする。一方、風向が１９０°以下や、２３４°以上の場合は、風速補正値ＰＣ１に対する加重比ｐｃ１を１．０とし、風速補正値ＰＣ２に対する加重比ｐｃ２を０．０とする。これらは上記「ｐｃ１＋ｐｃ２」＝‘１’となることは明らかである。
【０１０７】
また、風向が１９０°〜２３４°の範囲内（但し、２１２°は除外する）の場合は、加重比ｐｃ１、ｐｃ２は、それぞれ例えば図示のような０より大きく１未満の範囲内の該当する値（風向に対応する値）となる。これらについては、特に具体的な数値等は示さないが、上記の通り、基本的には「ｐｃ１＋ｐｃ２」＝‘１’の条件を満たすように規定されている。
【０１０８】
尚、２１２°は上記の通り風上側風車が設けられている方向であるが、１０°や１９０°や２３４°は、例えば開発者等が適宜決定した値である。
発電量計算機能部１４は、気象予報入力機能部１１を介して得られる上記風向予報値に基づいて、図５（ａ）に示すメンバシップ関数を参照することで、この風向予報値に対応する加重比ｐｃ１、加重比ｐｃ２を求める。例えば、風向きが２１２°であった場合には、加重比ｐｃ１＝‘０．０’、加重比ｐｃ２＝‘１．０’となる。
【０１０９】
以上、風速予報に応じた風速補正値ＰＣ１、ＰＣ２と、風向予報に応じた、これら風速補正値ＰＣ１、ＰＣ２に対する加重比ｐｃ１、ｐｃ２とを取得したら、これらに基づいて下記の（３）式によって最終的な風速補正値ＰＣを求める。
【０１１０】
風速補正値ＰＣ＝（風速補正値PC1 ×加重比pc1 ＋風速補正値PC2 ×加重比pc2 ）／（加重比pc1＋加重比pc2）・・・（３）式
その後は、例えば図３（ａ）等に示すパワーカーブを用いて、上記（３）式で求めた風速補正値ＰＣに応じた発電量を求める。
【０１１１】
なお、以上の実施例２の説明では、説明を簡単にするため、風上側風車の稼働率などを考慮していないが、実施例１と略同様に、風上側風車の稼働率などを考慮してもよい。仮に、風上側風車が９０％稼働しているとした場合は、発電量予測値は以下のようになる。尚、以下の式は、実質的には上記（４）式と略同様であり、２つの発電量を、風上側風車の稼働率に応じた加重比（９０：１０の加重比）で平均化するものである(これより、(４)’式と記す)。
【０１１２】
発電量予測値＝発電量１×０．９＋発電量２×０．１ ・・・（４）’式
｛ここで、
発電量１：風上側風車が稼働しているときの発電量（式（３）と図４（ａ）の補正テーブルから計算した風速補正値を、図３（ａ）のパワーカーブで発電量に変換して求めたもの）。
【０１１３】
発電量２：風上側風車が停止しているときの発電量（式（３）と図４（ｂ）の補正テーブルから計算した風速補正値を、図３（ａ）のパワーカーブで発電量に変換して求めたもの）。｝
また、図４（ａ），（ｂ）の例では、簡単のため、風速補正に係る風向区分を上記“２１２°±１０°”と“それ以外”の２種類だけとしたが、勿論、この様な例に限らず、３種類以上であっても構わない。例えば図３では西方向に何らかの障害物があるものとしており、この例に準じて、“２７０°±１０°”等もあってもよい。この様な他の例を図４（ｃ）、（ｄ）に示す。
【０１１４】
図４（ｃ）は風上側風車が運転（ＯＮ）している場合に対応する風速補正テーブルである。
図４（ｄ）は風上側風車が停止（ＯＦＦ）している場合に対応する風速補正テーブルである。
【０１１５】
図示の通り、これら他の例の風速補正テーブルでは、“２１２°±１０°”と“それ以外”と“２７０°±１０°”の３種類の風向区分（風向範囲）毎に、それぞれ、各風速予報値に応じた補正値が格納されている。勿論、補正値の代わりに上記風速補正値が格納されていてもよい。尚、“２７０°±１０°”に関しては、風上側風車の運転状態は関係ないので、図示の通り図４（ｃ）と図４（ｄ）とで同じ内容となっている。
【０１１６】
尚、本例の場合、“それ以外”は、「“２１２°±１０°”及び“２７０°±１０°”」を除く全ての方向を意味する。
また、この場合には、当然、メンバシップ関数も、３種類の風向区分（風向範囲）に応じたものとなる。図５（ｂ）に、上記図４（ｃ）、（ｄ）の例に応じたメンバシップ関数の一例を示す。図５（ｂ）に示すメンバシップ関数は、上記“２１２°±１０°”に応じた“風速補正の加重比ｐｃ２”と、上記“２７０°±１０°”に応じた“風速補正の加重比ｐｃ３”と、上記“それ以外”に応じた“風速補正の加重比ｐｃ１”とから成る。
【０１１７】
この場合も、上記図５（ａ）に示すメンバシップ関数と略同様に、基本的には、全ての風向に関して「ｐｃ１＋ｐｃ２＋ｐｃ３」＝‘１’となるように、加重比ｐｃ１、ｐｃ２、ｐｃ３が規定されている。
【０１１８】
そして、図４（ａ），（ｂ）及び図５（ａ）の場合と略同様にして、ｍ時間先の風向・風速予報値に基づいて、図４（ｃ）または図４（ｄ）の風速補正テーブルから風速予報値に応じた各風向範囲毎の風速補正値ＰＣ１，ＰＣ２，ＰＣ３を取得すると共に、図５（ｂ）のメンバシップ関数から風向予報値に対応する各風向範囲毎の加重比ｐｃ１、ｐｃ２、ｐｃ３を取得する。
【０１１９】
尚、当然、風上側風車の運転スケジュールに応じて、図４（ｃ）と図４（ｄ）のどちらの風速補正テーブルを使用するのか（あるいは両方とも使用するのか）が決まることになる。
【０１２０】
そして、これら取得したデータに基づいて、上記（３）式と略同様の算出式によって最終的な風速補正値ＰＣを求める。つまり、本例の場合には、下記の（３）’式によって最終的な風速補正値ＰＣを求めることになる。
【０１２１】
風速補正値ＰＣ＝（風速補正値PC1×加重比pc1＋風速補正値PC2×加重比pc2＋風速補正値PC3×加重比pc3）／（加重比pc1＋加重比pc2＋加重比pc3）・・・（３）’式
本実施例では、わずか数種類の風向区分に応じた補正テーブルやメンバシップ関数を予め作成しておくだけで、風向ごとに異なる風速の変化状態を連続的に表現することができる。
【０１２２】
次に、以下、実施例３について説明する。
上記実施例１、２では、風上側風車の存在とその稼働状況を考慮した実施例を示したが、本発明は、この様な例に限らない。例えば、本実施例３のようにしてもよい。
【０１２３】
図６（ａ）は、実施例３におけるパワーカーブの一例であり、これは図３（ａ）に示す実施例１のパワーカーブと、実質的に同じであってよく、ここでは特に説明しない。
図６（ｂ）は、実施例３における風速補正テーブルの一例である。実施例１では、風上側風車の運転／停止に応じた２種類の風速補正テーブルが設けられていたが、実施例３では図示の通り１種類のみである。尚、図示の例では、北方向に風車以外の何らかの障害物があるものとしている。
【０１２４】
実施例１では、上述したように、特定の障害物である風上側風車の運転／停止（ＯＮ／ＯＦＦ）に応じて、２種類の風速補正テーブルを設定して選択的に使用したが、実施例３ではこの様な選択処理は行われない。また、実施例１では風上側風車の稼働率に応じた補正処理等も行ってよいものとしたが、実施例３ではこの様な処理は行われない。また、当然、実施例３では風力発電設備運転情報取得部１２の処理機能は必要ない。これらの点を除けば、実施例３は、実施例１と略同様と見做してもよく、以下、簡単に説明する。
【０１２５】
すなわち、気象予報入力機能部１１を介してｍ時間先の気象予報（特に風向・風速予報値）が入力されると、発電量計算機能部１４は、これらに基づいて、図６（ｂ）に示す風速補正テーブルを参照することで、該当する補正値を得て、以って風速補正値を求める。気象予報値の一例として、例えば風向予報が北で、風速予報値が１０ｍの場合は、図６（ｂ）の風速補正テーブルより、該当する補正値（北で１０ｍの場合、図示の例では補正値＝‘−１’）を得ることで、上記風速補正値を以下のように求める。尚、風速補正値は、風力発電所における処理対象の風力発電設備（その風車）が実際に受けると考えられる風の風速である。
【０１２６】
風速補正値＝９ｍ（＝風速予報値＋補正値＝１０ｍ−１ｍ）
この場合、風力発電量予測値は、図６（ａ）のパワーテーブルより、風速９ｍに対応する９５ｋWとなる。
【０１２７】
なお、上記の説明では、簡単のため、加重平均やファジィ推論を用いなかったが、実施例３において実施例１、２のように加重平均やファジィ推論を用いるようにしてもよい。
例えば、風向予報が北北東で風速予報が１０ｍの場合には、図６（ａ）、（ｂ）に示すパワーテーブルと風速補正テーブルを参照することで、風速補正テーブルに登録されている東西南北の４方向のうち上記風向予報（北北東）に関連する方向である北と東に関して、以下の風速補正値及び発電量が得られる。
【０１２８】
北： 風速補正値９ｍ ， 発電量 ９０ｋW
東： 風速補正値１０ｍ ， 発電量１００ｋW
そして、これら２つの発電量を、風向予報に応じた加重比で平均化することで、風向・風速予報と風速補正テーブルに基づく発電量予測値が算出される。本例では、風向予報が北北東であることから、２５：７５の加重比で平均化することになり、従って、
発電量予測値；（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４ ＝ ９２．５ｋW ・・・（１）’’式
となる。
【０１２９】
尚、上述した例では、実施例１の（１）式、（１）’式も、上記実施例３の（１）’’式も、２つの発電量を求めた後に加重平均算出処理を行ったが、この例に限らない。例えば、上記の例では、北と東の風速補正値（９ｍと１０ｍ）に対して、上記風向予報（北北東）に応じて２５：７５の加重比で平均化することで、風向予報に応じた風速補正値を得るようにしてもよい。
【０１３０】
この例の場合、風向予報（北北東）に応じた風速補正値は、９．２５（ｍ）となる。そして、パワーテーブルを参照することで、風向予報（北北東）に応じた風速補正値（９．２５（ｍ））に対応する発電量を求め、これを上記発電量予測値とするようにしてもよい。但し、図３（ａ）や図６（ａ）に示す例のパワーカーブでは、風速が１ｍ刻みであるので、９．２５（ｍ）に対応する発電量は登録されていない。よって、パワーカーブがこの様な例の場合には、上記（１）式、（１）’式や（１）’’式によって求める手法の方が望ましい。
【０１３１】
この様に、風速補正テーブルには少数の所定方向（本例では東西南北の４方向）に関する補正データを登録しておくだけでよく、これら４方向以外の風向きに関しては加重平均することにより求めることができる。つまり、全ての風向に対してデータベースを持つ必要がないため、事前のデータベース設定が容易であり、エンジニアリング・コストを削減できる。
【０１３２】
あるいは、図６（ｂ）に示す補正テーブルの代わりに、図４に示すような、風向を０°〜３６０°の角度で扱って所定の角度範囲毎に対応する補正値を登録した風速補正テーブルを用いるようにしてもよい。勿論、この場合には、更に、図５に示すようなメンバシップ関数も予め作成・記憶して必要がある。
【０１３３】
尚、実施例３では風上側風車は想定しないので、風速補正テーブルは、図４（ａ）、（ｂ）のような２種類ではなく、１種類のみである。また、それ故、当然、風力発電設備運転情報取得部１２の処理機能は必要ないし、風上側風車のＯＮ／ＯＦＦに応じて２種類の風速補正テーブルの何れか一方（または両方）を選択する処理も必要ない。これらの点を除けば、本例の処理は、実施例２の処理と略同様と見做しても構わない。
【０１３４】
すなわち、実施例２の場合と略同様に、気象予報入力機能部１１を介してｍ時間先の気象予報（風向・風速予報値）を入力すると、風向・風速予報値と不図示の（図４と略同様の）風速補正テーブルに基づいて、上記風速補正値ＰＣ１，ＰＣ２等に相当する風速補正値を求める。更に、風向予報値と不図示の（図５と略同様の）メンバシップ関数を用いて、加重比ｐｃ１、ｐｃ２等に相当する加重比を求める。
【０１３５】
そして、これら求めた風速補正値と加重比とに基づいて、上記（３）式や（３）’式と略同様の算出式によって、最終的な風速補正値ＰＣを求める。その後は、実施例２と略同様に、例えば図３（ａ）等に示すパワーカーブを用いて、求めた風速補正値ＰＣに応じた発電量予測値を求める。
【０１３６】
この様に、実施例３において、ファジィ推論を用いて風速補正値を求め、以って風速発電量予測値を求めるようにしてもよい。
上述した実施例１，２，３に係る図１の風力発電量予測システム１０の各機能部の処理機能は、例えば、以下の様になると言える。
【０１３７】
気象予報入力機能部１１は、外部から風向・風速予報値を取得する。これは、例えば気象庁等による一般的な気象予報情報としての風向・風速予報値を取得するものである。
補正テーブル記憶部１６には、予め、上記風速予報値を補正するための補正用データが記憶される。パワーカーブ記憶部１３には、予め、各風力発電設備毎に、その風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶される。
【０１３８】
発電量計算機能部１４は、上記風向・風速予報値と補正用データとパワーカーブとを用いて、加重平均（実施例１，３）またはファジィ推論（実施例２）によって、風力発電量予測対象の任意の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１３９】
そして、例えば、実施例１，３の場合には、上記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルである。そして、発電量計算機能部１４は、この風速補正テーブルを参照して、風向予報が示す風向に係わる１以上の上記所定の風向毎に、風速予報値に対応する補正データを用いて風速補正値を求めると共に、パワーカーブを用いて、求めた１以上の風速補正値毎に対応する発電量を求める。更に、発電量計算機能部１４は、複数の発電量が求められた場合には、該複数の発電量を風向予報に応じた加重比で平均化することで、風力発電量予測対象の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１４０】
尚、各風向に応じた加重比は、例えば予め人間等が決めて設定・記憶させておけばよい。上述した一例では、複数の所定の風向は東西南北の４方向であり、風向予報がこれら４方向の何れかである場合以外は、複数の発電量が求められることになり、上記加重平均による風力発電量予測値の算出が行われることになる。
【０１４１】
尚、上記風速補正テーブルには、補正データの代わりに風速補正値が格納されていてもよい。この場合、発電量計算機能部１４は、上述した風速予報値に対応する補正データを用いて風速補正値を求める処理は行う必要なく、その代わりに、風速予報値に対応する風速補正値を風速補正テーブルから取得すればよい。
【０１４２】
また、例えば実施例２の場合には、上記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、上記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数である。そして、発電量計算機能部１４は、風速補正テーブルを参照して、各風向範囲毎に、風速予報値に応じた補正データを用いて風向範囲毎風速補正値を求めると共に、メンバシップ関数を参照して、各風向範囲毎に、風向予報値に応じた加重比を求める。更に、求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、風速予報値を補正した風速補正値を算出する。そして、該算出した風速補正値と上記パワーカーブとによって、風力発電量予測対象の風力発電設備が受ける風に応じた風力発電量予測値を求める。
【０１４３】
尚、上記所定の各風向範囲とは、例えば図４に示す“２１２°±１０°”、“２７０°±１０°”や“それ以外”等であり、この例では風向は３６０°で表されるものとなっている。
【０１４４】
尚、上記風速補正テーブルには、補正データの代わりに風向範囲毎風速補正値が格納されていてもよい。この場合、発電量計算機能部１４は、上述した風速予報値に対応する補正データを用いて風向範囲毎風速補正値を求める処理は行う必要なく、その代わりに、風速予報値に対応する風向範囲毎風速補正値を風速補正テーブルから取得すればよい。
【０１４５】
更に、風上側風車を考慮する場合には、風力発電設備運転情報取得部１２によって、各風力発電設備（少なくとも風上側風車）の運転予定情報を取得する。尚、既に述べたように、風上側風車とは、上記風力発電量予測対象の風力発電設備が受ける風（特にその風速）に影響を及ぼし得る他の風力発電設備（その風車）を意味する。
【０１４６】
この場合、上記風上側風車（風上側風力発電設備）が存在する場合、上記風速補正テーブルは、当該風上側風力発電設備の運転時、停止時それぞれに応じた２種類の風速補正テーブルとなる。
【０１４７】
そして、発電量計算機能部１４は、風上側風力発電設備の上記運転予定情報に基づいて、例えば、該２種類の風速補正テーブルの何れか一方を用いて、上記風速補正値または上記風向範囲毎風速補正値を求める。
【０１４８】
あるいは、発電量計算機能部１４は、上記２種類の風速補正テーブルの両方をそれぞれを用いることで、上記風力発電量予測値を２つ求め、該求めた２つの風力発電量予測値を、上記風上側風力発電設備の稼働率に応じた加重比で平均化することで、最終的な風力発電量予測値とする。
【０１４９】
続いて、以下、実施例４、実施例５について説明する。
まず、実施例４について説明する。
上述した実施例１，３では、パワーカーブと風速補正テーブルを予め作成・記憶しておき、発電量予測処理の際には、まず風速補正テーブルを参照して該当する風速補正値を求めて、次に、パワーカーブを参照することでこの風速補正値に応じた発電量を求めた。
【０１５０】
これに対して、実施例４では、後述するパワーカーブ補正テーブルを予め作成・記憶しておき、発電量予測処理の際には、パワーカーブ補正テーブルを参照することで、一般的な風向・風速予報値に基づき、各風力発電設備毎に的確な発電量予測値を求める。
【０１５１】
実施例４に関しても、図１に示す構成図を用いて説明する。
実施例４に関する図１の各処理機能部は、まず、気象予報入力機能部１１及び風力発電設備運転情報取得部１２は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【０１５２】
パワーカーブ記憶部１３と補正テーブル記憶部１６は削除され、これらの代わりに不図示のパワーカーブ補正テーブル記憶部が設けられる。このパワーカーブ補正テーブル記憶部には、上記パワーカーブ補正テーブル（後に図７（ａ）、（ｂ）に具体例を示す）が、予め記憶されている。
【０１５３】
当然、発電量計算機能部１４の処理内容は、実施例１とは異なる。発電量計算機能部１４は、気象予報入力機能部１１及び風力発電設備運転情報取得部１２によって取得される情報と、上記パワーカーブ補正テーブルとに基づいて、風速発電量予測値を求める。以下、具体例を挙げて説明する。
【０１５４】
図７（ａ）、（ｂ）は、上記パワーカーブ補正テーブルの具体例である。
図７（ａ）は風上側風車が運転（ＯＮ）している場合に対応するパワーカーブ補正テーブルである。
【０１５５】
図７（ｂ）は風上側風車が停止（ＯＦＦ）している場合に対応するパワーカーブ補正テーブルである。
パワーカーブ補正テーブルは、例えば、予め、上記「発電所のある代表地点、もしくは近辺の風速・風向の予報値または現在値」を取得すると共に、各風力発電設備の設置場所で風速等を実測する等し、更にパワーカーブを風車メーカから取得する（あるいは工場等において、風速に応じた発電量を実測する）等して、これらに基づいて人間が作成するものである。
【０１５６】
図７（ａ）、（ｂ）に示すように、パワーカーブ補正テーブルは、１ｍ単位の各風速（ｍ）毎に対応して、設計値（ｋｗ）と、少数の所定の風向（本例では東西南北の４方向）毎の発電量補正値（ｋｗ）（北（補正値）、東（補正値）、南（補正値）、西（補正値））が格納されている。尚、図示の「風速」は、図３等における「風速予報」と異なるようなものではなく、例えば気象予報入力機能部１１が取得する風速予報値が２（ｍ）で風向予報が‘北’である場合には、パワーカーブ補正テーブルにおける「風速」が２（ｍ）の北（補正値）である‘−１’が取得されることになる。よって、図示の「風速」は、「風速予報」に置き換えても良い。
【０１５７】
また、パワーカーブ補正テーブルにおける発電量（設計値）は、風力発電設備のメーカが工場等で計測した値である。つまり、図示の例では、風車が実際に受ける風速が２ｍであったならば、発電量は１０（ｋＷ）になるはずである。
【０１５８】
また、例えば上記北（補正値）とは、風向が北の場合に使用する、発電量に係る補正値であり、例えば図７（ａ）に示す例においては、風速２ｍのときには、本来は１０ｋW発電するが９ｋW（＝１０−１）になることを示している。同様に、風速１０ｍの場合には本来は１００ｋW発電するはずが９０ｋｗ（＝１００−１０）になることを示している。これは、他の風向の補正値（東（補正値）、南（補正値）、西（補正値））についても略同様であり、ここでは特に説明しない。尚、上記“本来は”とは、他の風車や風車以外の障害物の影響で、風速が変化しなければという意味である。
【０１５９】
ここで、図７（ａ）、（ｂ）のパワーカーブ補正テーブルは、上記実施例１の図３（ｂ）、（ｃ）に示す例の風速補正テーブルと略同様に、北に風上側風車、西に山などの障害物があることを想定している。そのため、風上側風車の運転状態（ＯＮかＯＦＦか）によってテーブルの内容（数値）が変化するのは、図示の通り北のみである（図示の通り、東、南、西に関しては、その数値は、図７（ａ）と図７（ｂ）とで同じである）。
【０１６０】
実施例４における発電量計算機能部１４は、上記風力発電設備運転情報取得部１２を介して得られる各風力発電設備の運転スケジュール（ＯＮ／ＯＦＦや稼働率等）と、上記気象予報入力機能部１１を介して得られる上記風向・風速予報値と、上記パワーカーブ補正テーブルとに基づいて、風力発電量予測値を求める。
【０１６１】
例えば、上記風向・風速予報値と各風力発電設備の運転スケジュールに基づいて、図７（ａ）、（ｂ）のパワーカーブ補正テーブルの少なくとも一方は用いて、風力発電量予測値を求める。
【０１６２】
たとえば、風向予報が北東で、風速予報が１０ｍ、かつ予測時間帯における風上側風車が運転予定が“運転（ＯＮ）”の場合は、図７（ａ）に示す（風上側風車が運転している場合の）パワーカーブ補正テーブルを参照することで、各所定風向（東西南北）毎の発電量補正値（ｋｗ）は、下記の通りとなる（尚、ここでは東西南北全てを示すが、風向予報に関係するものだけで構わない。上記一例のように風向予報が北東の例であれば、北と東だけでも構わない）。
【０１６３】
北の場合： ９０ｋW
東の場合：１００ｋW
南の場合：１００ｋW
西の場合：１００ｋW
尚、これらは、図７（ａ）、（ｂ）のパワーカーブ補正テーブルに予め登録されている設計値と補正値とに基づいて算出するものである。但し、この例に限らず、例えば、パワーカーブ補正テーブルには、設計値も補正値も登録されず、代わりに上記発電量補正値（ｋｗ）そのものが予め登録されていてもよい。
【０１６４】
風向予報が東西南北の何れかであれば、その風向きに応じた上記発電量補正値（ｋｗ）をそのまま演算結果（風力発電量予測値）として用いればよい。
風向予報が東西南北以外であれば、上記東西南北の発電量補正値（ｋｗ）のうち風向予報に関係する方向の２つの発電量補正値（ｋｗ）を用いて、これら２つの発電量補正値を風向予報に応じた加重比を用いて加重平均することで、風力発電量予測値を求める。例えば、風向予報が北北東であれば、上記北の場合の発電量補正値（＝９０ｋW）と上記東の場合の発電量補正値（＝１００ｋW）の２つの発電量補正値を用いることになる。
【０１６５】
そして、実施例１等と略同様の加重平均を用いる手法により、これら２つの発電量補正値から風力発電量予測値を求める。
例えば風向予報が北東である場合には、北東は北と東の丁度中間であることから、５０：５０の加重比で平均化することになり、下記の（６）式で風力発電量予測値を求めることになる。
【０１６６】
風力発電量予測値＝（１００ｋW×５０＋９０ｋW×５０）／１００
＝（１００ｋW＋９０ｋW）／２＝９５ｋW ・・・（６）式
となる。
【０１６７】
あるいは、上記の例において、仮に風向予報が北北東であったならば、２５：７５の加重比で平均化することになり、従って、
風力発電量予測値＝（１００ｋW×２５＋９０ｋW×７５）／１００
＝（１００ｋW＋９０ｋW×３）／４＝９２．５ｋW
となる。
【０１６８】
また、仮に予測時間帯における風上側風車の運転予定が停止（ＯＦＦ）であるならば、上記図７（ａ）のテーブルを用いる代わりに図７（ｂ）のテーブルを用いることになる。
また、例えば、風上側の風車の稼働率を考慮する場合には、図７（ａ）、（ｂ）のパワーカーブ補正テーブルそれぞれを用いて、上記のように風力発電量予測値を求めて、これら求めた２つの風力発電量予測値を、風上側風車の稼働率に応じた加重比で平均化することで、最終的な風力発電量予測値を求めることになる。仮に、風上側風車の稼働率が９０％とした場合には、上記２つの風力発電量予測値を、９０：１０の加重比で平均化することで、最終的な風力発電量予測値を求めることになる。
【０１６９】
また、更に、予測対象の風車の稼働率を用いて修正を施して、最終的な予測値とするようにしてもよい。仮に、予測対象の風車の稼働率が８０％の場合は、最終的な風力発電量予測値に対して０．８を乗じることで、最終的な予測値とする。
【０１７０】
尚、上述した処理と予測値は、未来の任意の時間帯（例えば翌日の２時〜３時の１時間分など）に関する発電量予測処理とその結果である。もし、気象予報が１時間ごとｍ時間先まである場合は、上記処理をｍ回行うことで、ｍ時間先まで求めることができる。
【０１７１】
本手法の利点として、風向に関するパワーカーブが少なくても、その間の風速は加重平均にて求めることができるので、精度の高い予測値を得ることができること、および少ないパワーカーブのテーブルだけで予測できるので、エンジニアリング・コストを削減できる点があげられる。
【０１７２】
次に、以下、実施例５について説明する。
本実施例５は、上記実施例２と同様にファジィ推論にて発電量を求めるものであるが、実施例２における図４に示すような風速補正テーブルの代わりに、不図示のパワーカーブ補正テーブルが予め記憶されている。
【０１７３】
不図示のパワーカーブ補正テーブルは、例えば図４のように風上側風車が運転ＯＮの場合とＯＦＦの場合とに応じて２つ用意されると共に、風向きを０°〜３６０°で扱って任意の角度範囲毎（例えば上記“２１２°±１０°等）に対応付けて、補正値が格納されているものである。但し、図４では補正値は風速に係る補正値であったが、パワーカーブ補正テーブルに登録される補正値は、図７と同様に発電量に係る補正値となる。更に、図７と同様に各風速毎の発電量（設計値）も登録される。
【０１７４】
但し、この例に限らず、実施例４と略同様に、補正値と設計値の代わりに、発電量補正値（補正後の発電量）そのものが登録されていてもよい。当然、この場合には、補正値と設計値とに基づいて発電量補正値（補正後の発電量）を算出する処理は必要ない。
【０１７５】
また、例えば図５に示す例のようなメンバシップ関数も、予め記憶されている。
ここで、実施例５に関しても、図１に示す構成図を用いて説明する。
実施例５に関する図１の各処理機能部は、まず、気象予報入力機能部１１及び風力発電設備運転情報取得部１２は、実施例１と略同様であり、ここでは特に説明しない。また、予測値格納部１５も、発電量計算機能部１４によって求められた発電量予測値が格納されるという点では、実施例１と略同様である。
【０１７６】
パワーカーブ記憶部１３と補正テーブル記憶部１６は削除され、これらの代わりに上述した不図示のパワーカーブ補正テーブルが予め記憶されるパワーカーブ補正テーブル記憶部が設けられる。
【０１７７】
パワーカーブ補正テーブルの内容は、例えば以下の２種類である。
通常時（風向２１２°以外の場合）（PC１）： 発電量PC1＝ｆPC１（風速）
風向２１２°の場合（PC２）： 発電量PC2＝ｆPC２（風速）
※ｆPC１、ｆPC２は、特に具体例は示さないが、任意の風速値に応じた発電量を算出するための関数等である。
【０１７８】
※風向は北風が０°で、時計まわり
尚、ここでは、２１２°の方向に、風上側風車が存在する場合を想定している。
発電量計算機能部１４は、まず、上記２つの発電量ＰＣ１、ＰＣ２を求める。これは、気象予報入力機能部１１を介して得られた風速予報値を上記（風速）に代入することで算出される。
【０１７９】
次に、図５（ａ）に示すメンバシップ関数を参照して、気象予報入力機能部１１を介して得られた風向予報に応じた加重比を求める。
つまり、図５（ａ）の例のメンバシップ関数を参照する場合には、予報される風向が２１２°の場合は、上記発電量ＰＣ２に対する加重比ｐｃ２を１．０とし、上記発電量ＰＣ１に対する加重比ｐｃ１を０．０とする。一方、予報される風向が１９０°以下や、２３４°以上の場合は、上記発電量ＰＣ１に対する加重比ｐｃ１を１．０とし、上記発電量ＰＣ２に対する加重比ｐｃ２を０．０とする。風向が１９０°から２３４°の範囲内の任意の値である場合は、図５（ａ）のメンバシップ関数から読み取った加重比ｐｃ１、ｐｃ２を用いる
最終的な発電量予測値ＰＣは次式で求める
発電量予測値ＰＣ＝（発電量ＰＣ１×加重比ｐｃ１＋発電量ＰＣ２×加重比ｐｃ２）／（加重比ｐｃ１＋加重比ｐｃ２）
図８は、本実施例５の効果を示す為の実験結果である。
【０１８０】
図８のグラフは、横軸が時間、縦軸が任意の風力発電設備の発電出力（ｋＷ）である。但し、グラフには風向もあり、この場合、縦軸は右側に示すように０°〜３６０°で表される風向となる。
【０１８１】
図８には、発電出力の実績値を実線で示し、発電出力の予測値を点線または一点鎖線で示す。点線は本実施例５によるファジイ推論を適用した場合であり、一点鎖線はファジイ推論を適用しなかった場合（従来技術）である。
【０１８２】
これより、図８に示すように、ファジィ推論未使用時の発電量は実績値よりも大きいのに対し、ファジィ推論適用時は風向が２１２°付近のときに発電出力が低下し実績値に近い値になっている。
【０１８３】
なお、本実施例５では、説明を簡単にするため単純なファジィ推論により加重比を求めたが、複数のパワーカーブを用いてもよいし、風向だけでなく風上側風車の運転状態も入力条件にしてもよい。
【０１８４】
本実施例５の上記一例では、パワーカーブ補正テーブルにわずか２種類の発電量補正PC1、PC2を用意しただけで、風向ごとに異なる発電量の低下状態を精度よく表現することができる。
【０１８５】
以上説明した本手法の特に実施例１，２，３によれば、以下の効果が得られる。
風力発電予測に関して、風速の予測と、発電量への換算の２段階で行われる。たとえ風速の予測が正確であっても発電量の換算が正確でないと適切な発電量予測値が得られない。その逆に、発電量の換算が正確であっても風速の予測が正確でないと適切な発電量予測値が得られない。
【０１８６】
ここで、例えば日本のような山が多い地形にある風力発電所では、実際に風車が受ける風は、気象予報どおりの風速が得られない場合があった。
これに対して、これまでの風力発電量予測では、特殊な計測装置により実測する必要があったり、風向毎に細かな補正を行うためには多くのテーブルを用意する必要があり、エンジニアリング・コストが掛かった。また、従来の物理モデルを用いる方法は、高精度であるが、膨大な地形情報や計算コストが必要であり、やはり高コストであった。
【０１８７】
本手法では、既設の風力発電所の計測データと一般的な気象予報／気象観測データだけでデータベースを構築でき、特に所定方向（東西南北など）に関する補正情報だけを登録すればよく、更に運用時には一般的な気象予報等に基づいて、データベースを参照した補正を行うことで、各風力発電設備毎に受ける風の風速やそれに応じた発電量を、的確に予測することが可能であり、低コストで高精度の風力発電量予測を可能にする。
【０１８８】
本手法では、一般的な気象予報データと補正用データだけで各風力発電設備毎の風力発電量を的確に予測でき、特に加重平均またはファジィ推論を用いることで補正用データのデータ量削減を図り、以ってエンジニアリング・コストを低減できる。特に、風向毎に補正するためのデータベースは有するが、全ての風向に対してデータベースを持つ必要がないため、事前のデータベース設定が容易であり、エンジニアリング・コストを削減でき、以って低コストで高精度な予測値を得ることができる。
【０１８９】
本手法を用いることにより、容易・かつ正確に発電量への換算が可能になり、風力発電量予測の精度向上に寄与できる。
更に、特にウィンドファームなどの複数の風車がある場合は、風上側風車の運転状態によって風下側風車が受ける風速が変化するため、風下側風車に関して適切な発電量予測が難しかった。
【０１９０】
これに対して、本手法では、上記特徴・効果に加えて更に、上記データベースとして風上側風車の稼働状態に応じた複数種類のテーブルを設け、風上側風車の運転スケジュールに基づいて該当するテーブルを用いて風速補正を行うことで、風上側風車の稼働状態まで考慮して風力発電量を的確に予測することができる。
【０１９１】
風上側風車に関しては、それが運転中（ＯＮ）の場合とメンテナンス等で停止中（ＯＦＦ）の場合とで、予測対象の風車が受ける風速に与える影響が異なる。この点を想定して、風上側風車がＯＮの場合とＯＦＦの場合とに応じた２種類の補正テーブルを用意することで、より的確な風力発電量予測を実現できる。
【０１９２】
この様に、本手法は、山や建造物等による影響を考慮して、一般的な気象予報による風速予報値を各風力発電設備毎に補正することで、各風力発電設備（その風車）が実際に受ける風速を的確に予測し、以って発電量を的確に予測することができる。
【０１９３】
また、風速／発電量を補正するためのデータベースは、全ての風向に対する補正用データを登録しなくてもよいので（例えば東西南北の４方向のみ登録）、事前のデータベース設定が容易であり、エンジニアリング・コストを削減することが可能である。
【０１９４】
このように、本発明は、既存の風力発電所設備で得られる観測情報と一般的な気象予報だけで的確な風力発電量予測を実現でき、低コストで高精度な予測値を得ることができる。
【符号の説明】
【０１９５】
１０ 風力発電量予測システム
１１ 気象予報入力機能部
１２ 風力発電設備運転情報取得部
１３ パワーカーブ記憶部
１４ 発電量計算機能部
１５ 予測値格納部
１６ 補正テーブル記憶部
【特許請求の範囲】
【請求項１】
任意の風力発電設備の発電量を予測する風力発電量予測システムであって、
外部から風向・風速予報値を取得する風向・風速予報値取得手段と、
予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、
予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、
前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段と、
を有することを特徴とする風力発電量予測システム。
【請求項２】
前記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルであり、
前記風力発電量算出手段は、
前記風速補正テーブルを参照して、前記風向予報が示す風向に係わる１以上の前記所定の風向毎に、前記風速予報値に対応する前記補正データを用いて風速補正値を求めると共に、前記パワーカーブを用いて、前記求めた１以上の風速補正値毎に対応する発電量を求め、
複数の発電量が求められた場合には、該複数の発電量を前記風向予報に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項１記載の風力発電量予測システム。
【請求項３】
前記風速補正テーブルには、前記補正データの代わりに前記風速補正値が格納されており、
前記風力発電量算出手段は、前記風速予報値に対応する前記補正データを用いて風速補正値を求める代わりに、前記風速予報値に対応する前記風速補正値を前記風速補正テーブルから取得することを特徴とする請求項２記載の風力発電量予測システム。
【請求項４】
前記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、前記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数であり、
前記風力発電量算出手段は、
前記風速補正テーブルを参照して、前記各風向範囲毎に、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求めると共に、前記メンバシップ関数を参照して、前記各風向範囲毎に、前記風向予報値に応じた前記加重比を求め、
前記求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、前記風速予報値を補正した風速補正値を算出し、
該算出した風速補正値と前記パワーカーブとによって、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項１記載の風力発電量予測システム。
【請求項５】
前記風速補正テーブルには、前記補正データの代わりに前記風向範囲毎風速補正値が格納されており、
前記風力発電量算出手段は、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求める代わりに、前記風速予報値に対応する前記風向範囲毎風速補正値を前記風速補正テーブルから取得することを特徴とする請求項４記載の風力発電量予測システム。
【請求項６】
前記風向は、３６０°で表されることを特徴とする請求項４または５記載の風力発電量予測システム。
【請求項７】
前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備が存在する場合、前記補正用データ記憶手段に記憶される前記風速補正テーブルは、当該他の風力発電設備の運転時、停止時それぞれに応じた２つの風速補正テーブルであり、
前記風力発電量算出手段は、前記他の風力発電設備の運転予定情報に基づいて、該２つの風速補正テーブルの何れか一方を用いて、前記風速補正値または前記風向範囲毎風速補正値を求めることを特徴とする請求項２または４記載の風力発電量予測システム。
【請求項８】
前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備が存在する場合、前記補正用データ記憶手段に記憶される前記風速補正テーブルは、当該他の風力発電設備の運転時、停止時それぞれに応じた２つの風速補正テーブルであり、
前記風力発電量算出手段は、該２つの風速補正テーブルそれぞれを用いることで、前記風力発電量予測値を２つ求め、該求めた２つの風力発電量予測値を、前記他の風力発電設備の稼働率に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項２または４記載の風力発電量予測システム。
【請求項９】
任意の風力発電設備の発電量を予測する風力発電量予測システムのコンピュータを、
外部から風向・風速予報値を取得する風向・風速予報値取得手段と、
予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、
予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、
前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段、
として機能させるためのプログラム。
【請求項１】
任意の風力発電設備の発電量を予測する風力発電量予測システムであって、
外部から風向・風速予報値を取得する風向・風速予報値取得手段と、
予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、
予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、
前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段と、
を有することを特徴とする風力発電量予測システム。
【請求項２】
前記補正用データは、複数の所定の風向毎に、各風速に対応する補正データが格納された風速補正テーブルであり、
前記風力発電量算出手段は、
前記風速補正テーブルを参照して、前記風向予報が示す風向に係わる１以上の前記所定の風向毎に、前記風速予報値に対応する前記補正データを用いて風速補正値を求めると共に、前記パワーカーブを用いて、前記求めた１以上の風速補正値毎に対応する発電量を求め、
複数の発電量が求められた場合には、該複数の発電量を前記風向予報に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項１記載の風力発電量予測システム。
【請求項３】
前記風速補正テーブルには、前記補正データの代わりに前記風速補正値が格納されており、
前記風力発電量算出手段は、前記風速予報値に対応する前記補正データを用いて風速補正値を求める代わりに、前記風速予報値に対応する前記風速補正値を前記風速補正テーブルから取得することを特徴とする請求項２記載の風力発電量予測システム。
【請求項４】
前記補正用データは、所定の各風向範囲毎に、各風速に対応する補正データが格納された風速補正テーブルと、前記各風向範囲毎に、各風向に対応する加重比が規定されたメンバシップ関数であり、
前記風力発電量算出手段は、
前記風速補正テーブルを参照して、前記各風向範囲毎に、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求めると共に、前記メンバシップ関数を参照して、前記各風向範囲毎に、前記風向予報値に応じた前記加重比を求め、
前記求めた複数の風向範囲毎風速補正値および複数の加重比を用いて、前記風速予報値を補正した風速補正値を算出し、
該算出した風速補正値と前記パワーカーブとによって、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項１記載の風力発電量予測システム。
【請求項５】
前記風速補正テーブルには、前記補正データの代わりに前記風向範囲毎風速補正値が格納されており、
前記風力発電量算出手段は、前記風速予報値に応じた前記補正データを用いて風向範囲毎風速補正値を求める代わりに、前記風速予報値に対応する前記風向範囲毎風速補正値を前記風速補正テーブルから取得することを特徴とする請求項４記載の風力発電量予測システム。
【請求項６】
前記風向は、３６０°で表されることを特徴とする請求項４または５記載の風力発電量予測システム。
【請求項７】
前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備が存在する場合、前記補正用データ記憶手段に記憶される前記風速補正テーブルは、当該他の風力発電設備の運転時、停止時それぞれに応じた２つの風速補正テーブルであり、
前記風力発電量算出手段は、前記他の風力発電設備の運転予定情報に基づいて、該２つの風速補正テーブルの何れか一方を用いて、前記風速補正値または前記風向範囲毎風速補正値を求めることを特徴とする請求項２または４記載の風力発電量予測システム。
【請求項８】
前記任意の風力発電設備が受ける風に影響を及ぼし得る他の風力発電設備が存在する場合、前記補正用データ記憶手段に記憶される前記風速補正テーブルは、当該他の風力発電設備の運転時、停止時それぞれに応じた２つの風速補正テーブルであり、
前記風力発電量算出手段は、該２つの風速補正テーブルそれぞれを用いることで、前記風力発電量予測値を２つ求め、該求めた２つの風力発電量予測値を、前記他の風力発電設備の稼働率に応じた加重比で平均化することで、前記風力発電設備が受ける風に応じた風力発電量予測値を求めることを特徴とする請求項２または４記載の風力発電量予測システム。
【請求項９】
任意の風力発電設備の発電量を予測する風力発電量予測システムのコンピュータを、
外部から風向・風速予報値を取得する風向・風速予報値取得手段と、
予め前記風速予報値を補正するための補正用データが記憶される補正用データ記憶手段と、
予め前記風力発電設備に係わる風速と発電量との関係を示すパワーカーブが記憶されるパワーカーブ記憶手段と、
前記風向・風速予報値と前記補正用データと前記パワーカーブとを用いて、加重平均またはファジィ推論によって、前記任意の風力発電設備が受ける風に応じた風力発電量予測値を求める風力発電量算出手段、
として機能させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２０１３−１０８４６２（Ｐ２０１３−１０８４６２Ａ）
【公開日】平成２５年６月６日（２０１３．６．６）
【国際特許分類】
【出願番号】特願２０１１−２５５５６９（Ｐ２０１１−２５５５６９）
【出願日】平成２３年１１月２２日（２０１１．１１．２２）
【出願人】（０００００５２３４）富士電機株式会社 (3,146)
【Ｆターム（参考）】
【公開日】平成２５年６月６日（２０１３．６．６）
【国際特許分類】
【出願日】平成２３年１１月２２日（２０１１．１１．２２）
【出願人】（０００００５２３４）富士電機株式会社 (3,146)
【Ｆターム（参考）】
[ Back to top ]