(57)【要約】【課題】風予測情報を取得し、予測発電量を出力できる風力発電量予測システムを提供する。【解決手段】風力発電量予測システムについて、地点別に風向および風速を判別可能な風観測情報を取得する地点別風観測情報取得手段と、地点別に風向および風速を判別可能な風予測情報を取得する地点別風予測情報取得手段と、風力発電装置のタービン出力特性を取得するタービン出力特性取得手段と、前記地点別風観測情報取得手段、前記地点別風予測情報取得手段、および前記タービン出力特性取得手段により取得した情報を親ノードとして条件付確率表を用いて処理し予測発電量を子ノードとして出力するベイジアンネットワーク処理手段とを備え、前記ベイジアンネットワーク処理手段は、季節毎の風速安定度をノードとして用いる構成とした。【選択図】図８
2017102760.tif 000002
【特許請求の範囲】
【請求項１】
地点別に風向および風速を判別可能な風観測情報を取得する地点別風観測情報取得手段と、
地点別に風向および風速を判別可能な風予測情報を取得する地点別風予測情報取得手段と、
風力発電装置のタービン出力特性を取得するタービン出力特性取得手段と、
前記地点別風観測情報取得手段、前記地点別風予測情報取得手段、および前記タービン出力特性取得手段により取得した情報を親ノードとして条件付確率表を用いて処理し予測発電量を子ノードとして出力するベイジアンネットワーク処理手段とを備え、
前記ベイジアンネットワーク処理手段は、季節毎の風速安定度をノードとして用いる構成である
風力発電量予測システム。
【請求項２】
前記ベイジアンネットワーク処理手段は、極端気象の発生確率をノードとして用いる構成である
請求項１記載の風力発電量予測システム。
【請求項３】
各地点の風速値とバリアンス値を求め、風速値が大きくバリアンス値が小さい地点を前記風力発電装置の風車の設置に最適な風車設置好適地点として出力する風車設置好適地点出力手段を備えた
請求項１または２記載の風力発電量予測システム。
【請求項４】
前記ベイジアンネットワーク処理手段は、前記風力発電装置による予測発電量を子ノードとして用い、
予測発電量を出力する予測発電量出力手段を備えた
請求項１、２、または３記載の風力発電量予測システム。
【請求項５】
コンピュータを、
地点別に風向および風速を判別可能な風観測情報を取得する地点別風観測情報取得手段と、
地点別に風向および風速を判別可能な風予測情報を取得する地点別風予測情報取得手段と、
風力発電装置のタービン出力特性を取得するタービン出力特性取得手段と、
前記地点別風観測情報取得手段、前記地点別風予測情報取得手段、および前記タービン出力特性取得手段により取得した情報を親ノードとして条件付確率表を用いて処理し予測発電量を子ノードとして出力するベイジアンネットワーク処理手段として機能させ、
前記ベイジアンネットワーク処理手段を、季節毎の風速安定度をノードとして用いる構成である
風力発電量予測プログラム。
【請求項６】
地点別風観測情報取得手段により地点別に風向および風速を判別可能な風観測情報を取得し、
地点別風予測情報取得手段により地点別に風向および風速を判別可能な風予測情報を取得し、
タービン出力特性取得手段により風力発電装置のタービン出力特性を取得し、
前記地点別風観測情報取得手段、前記地点別風予測情報取得手段、および前記タービン出力特性取得手段により取得した情報を親ノードとして条件付確率表を用いて処理し、ベイジアンネットワーク処理手段が予測発電量を子ノードとして出力する
風力発電量予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えば、風力発電の発電量を予測するような風力発電量予測システム、風力発電量予測プログラム、および風力発電量予測方法に関する。
【背景技術】
【０００２】
従来、風力によって発電を行う様々な風力発電装置が提供されている。この風力発電装置は、海洋に設置されるもの、山林に設置されるもの、都市部に設置されるものなど、様々な種類が提供されている。
【０００３】
このような風力発電装置は、得られる風力によって発電量が変化する。このため、風力発電装置を設置する場所は、発電量を増大する上で重要である。
【０００４】
このような風力発電装置を有効に活用するために、風の向きや風の強さの変化を予測し、姿勢やプロペラ角度を最適化する風力発電システムが提案されている（特許文献１参照）。
【０００５】
この風力発電システムは、複数の風力発電装置を用い、各風力発電装置において風向、風力、温度、湿度、発電量などを測定して相互にやりとりし、自身の配置場所における風の変化タイミングを予測して姿勢、プロペラ角度等を最適化するものである。これにより、発電効率を常時最大に維持できるとされている。
【０００６】
しかし、この風力発電システムは、風力発電装置が複数設置されなければ予測できないものであった。このため、例えば、１台だけ風力発電装置を設置する際に、最も発電量が多くなる場所を特定するといったことはできなかった。また、発電した電力を効率良く利生するために発電量を事前に予測するといったことに利用できるものでもなかった。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０００−８７８４１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
この発明は、上述の問題に鑑みて、風予測情報を取得し、予測発電量を出力できる風力発電量予測システム、風力発電量予測プログラム、および風力発電量予測方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
この発明は、地点別に風向および風速を判別可能な風観測情報を取得する地点別風観測情報取得手段と、地点別に風向および風速を判別可能な風予測情報を取得する地点別風予測情報取得手段と、風力発電装置のタービン出力特性を取得するタービン出力特性取得手段と、前記地点別風観測情報取得手段、前記地点別風予測情報取得手段、および前記タービン出力特性取得手段により取得した情報を親ノードとして条件付確率表を用いて処理し予測発電量を子ノードとして出力するベイジアンネットワーク処理手段とを備え、前記ベイジアンネットワーク処理手段は、季節毎の風速安定度をノードとして用いる構成である風力発電量予測システムであることを特徴とする。
【発明の効果】
【００１０】
この発明により、風予測情報を取得し、予測発電量を出力できる風力発電量予測システム、風力発電量予測プログラム、および風力発電量予測方法を提供できる。
【図面の簡単な説明】
【００１１】
【図１】予測システムのシステム構成を示すブロック図。
【図２】風力発電量予測システムのベイジアンネットワークのブロック図。
【図３】消費電力予測システムのベイジアンネットワークのブロック図。
【図４】ベイジアンネットワークのブロック図。
【図５】条件付確率表（ＣＰＴ）の構成を示す説明図。
【図６】ベイズフィルタ推定の説明図。
【図７】風車設置場所選定処理を実行する制御部の動作を示すフローチャート。
【図８】ウインドファームの設置・運用を行うシステムのシステム構成図。
【図９】風力発電量予測と消費電力量予測のインバランスを最小に抑えるシステムのシステム構成図。
【図１０】最適な需要供給バランスを維持するシステムのシステム構成図。
【発明を実施するための形態】
【００１２】
以下、この発明の一実施形態を図面と共に説明する。
【実施例１】
【００１３】
図１は、予測システム１（風力発電量予測システム）のシステム構成を示すブロック図である。  
予測システム１は、風力発電予測システム１Ａ、数値予測モデルによる予測システム３、気象観測システム４、風力発電装置５、スマートメータ６、ユーザ端末８、およびインターネット９により構成されている。
【００１４】
また、風力発電予測システム１Ａは、ドップラーライダー２（地点別風観測情報取得手段，地点別風予測情報取得手段）および風力発電予測サーバ１０（風力発電予測装置）により構成されている。
【００１５】
ドップラーライダー２は、大気中にレーザー光を発射し、エアロゾルからの散乱光（反射光）のドップラーシフトを検出することで、エアロゾルの動きを基準に風況（風向き、風速）を観測する。このドップラーライダー２は、降雨の有無に関係なく風況を観測することができる。また、ドップラーライダー２は、レーザー光を用いてエアロゾルの観測を行うため、空間分解能及び時間分解能に優れた観測を行うことができる。
【００１６】
このドップラーライダー２は、複数台（例えば２台以上が好ましく、３台以上がより好ましい）が互いに所定の距離離間した略同一高さにそれぞれ設置されている。各ドップラーライダー２は、観測された視線風速に基づいて風況（風向及び風速）の水平分布情報を算出する。
【００１７】
各ドップラーライダー２は、その視線方向（レーザー光の発信方向）が水平方向に３６０°回転するように構成されており、約１分間に１周（３６０°）の周期で定速回転しつつ、例えば、半径約２０ｋｍの範囲の視線風速を連続的（１分毎）に観測する。
【００１８】
各ドップラーライダー２で得られた視線風速のデータは、ドップラーライダー２が備える送信器を介して予測サーバ１０に送信される。
【００１９】
予測サーバ１０の制御部１１（ベイジアンネットワーク処理手段、風車設置好適地点出力手段）は、各ドップラーライダー２の観測結果を合成し、観測領域における風況情報として、例えば、東西方向の風速画像（東西風速画像）及び南北方向の風速画像（南北風速画像）を取得する。
【００２０】
ここで、風況分布情報を生成する方法としては、種々の手法を採用することができ、例えば、３台以上のドップラーライダー２で観測されたデータに基づいて分布情報を生成する場合には最尤推定法により算出すればよく、また、２台のドップラーライダーで観測されたデータに基づいて分布情報を生成する場合には、線形重み付き最小二乗法により算出すればよい。このようにして風況分布情報を算出し、風況パターン情報を作製した場合には、１分間毎に風況パターン情報を記憶部１２に記憶すると良い。
【００２１】
上記風況分布情報は、適宜の水平分解能（例えば約１００ｍ程度など）とすることができる。また、このドップラーライダー２は、部品の汎用性の高さからよく使われている波長１．５μｍ帯の赤外線レーザーを用いたドップラーライダーなど、適宜のドップラーライダーを使用することができる。
【００２２】
また、各ドップラーライダー２は、設置高度の差が０〜１０ｍ程度になるように設置されることが好ましい。この程度の高度差であれば測定高度が異なっていても、風況に大きな差がないからである。
【００２３】
２台のドップラーライダー２で観測領域内を観測した場合、観測領域内の風向きによっては測定結果に誤差が生じる場合があるのに対し、３台以上のドップラーライダーを用いて観測領域内を観測することにより、より少ない誤差で観測領域の風況を計測することができるからである。
【００２４】
数値予測モデルによる予測システム３は、確率論的予測により大きなスケールでの気象予報を行うシステムである。
【００２５】
気象観測システム４は、天気、降雨量、風速、風向等の気象情報を観測する。
【００２６】
風力発電装置５は、風力によって風車が回転し、この回転により発生した電力を供給する発電装置として機能する。
【００２７】
スマートメータ６は、電力使用量を所定タイミング（例えば毎日、毎月、毎年等）等で送信する。
【００２８】
ユーザ端末８は、パーソナルコンピュータ、タブレット端末、あるいはスマートフォン等の適宜のコンピュータにより構成されており、通信を行う通信部と、入力を受け付ける入力部と、表示を行う表示部と、情報を記憶する記憶する記憶部と、各種制御動作を実行する制御部とを有している。
【００２９】
インターネット９は、公衆電気通信回線であり、各種装置を相互に通信可能に接続している。
【００３０】
風力発電予測サーバ１０は、制御部１１、記憶部１２、通信部１３（タービン出力特性取得手段、予測発電量出力手段）、入力部１４、及び表示部１５を備えている。
【００３１】
制御部１１は、ＣＰＵ、ＲＯＭ、ＲＡＭにより構成され、記憶部１２に記憶されている各種プログラムを起動して動作し、記憶部１２に記憶されている各種データおよび通信部１３から取得した各種データを用いて各種の演算と各部の制御動作を実行する。
【００３２】
記憶部１２は、風力発電量予測プログラム２１、消費電力量予測プログラム２２、風車設置場所選定プログラム２３、ベイズノード２４、および条件付確率分布表２５を記憶している。
【００３３】
通信部１３は、制御部１１の制御に従ってインターネット９に接続された各種装置とデータを送受信する。
【００３４】
入力部１４は、キーボードとマウス、あるいはタッチパネル等の入力装置によって構成され、ユーザによる操作入力を受け付けて入力信号を制御部１１へ送信する。
【００３５】
表示部１５は、ＣＲＴモニタあるいは液晶モニタ等の各種モニタ装置によって構成されており、制御部１１から出力される文字、図形、写真等の画像を表す画像信号に基づいて画像を表示する。
【００３６】
風力発電量予測プログラム２１は、風の予測と、風力発電装置５による発電量の予測を実行する。
【００３７】
消費電力量予測プログラム２２は、風力発電装置５により発電された電力を消費する側の消費電力量の予測を実行する。
【００３８】
風車設置場所選定プログラム２３は、風力発電装置５の風車をどの場所へ設置すれば多くの発電量を得られるかを選定する演算を実行する。
【００３９】
ベイズノード２４は、多数のベイズノード群を記憶している。  
条件付確率分布表２５は、ベイジアンネットワークに必要な条件付確率分布表を複数記憶している。
【００４０】
図２は、風力発電量を予測する風力発電量予測システムとして機能する際に必要なベイジアンネットワーク３０Ａの構成を示すブロック図であり、図３は、消費電力を予測する消費電力予測システムとして機能する際に必要なベイジアンネットワーク３０Ｂの構成を示すブロック図である。
【００４１】
風力発電量予測システムとして機能するベイジアンネットワーク３０Ａは、子がない親ノードとなる入力ノード群（３１〜４２）と、子ノードとなる出力ノード群（４６〜５１）と、自ノードとなるベイズノード群（４３〜４５）とを備えている。
【００４２】
親ノードとなる入力ノード群（３１〜４２）は、ドップラーライダー観測データ３１（風観測情報）に含まれている風向分布３２および風速分布３３等、数値予報モデルの予測データ３４（風予測情報）に含まれている風向・風速３５および気温・湿度３６等、タービン出力特性データ３７に含まれている最大出力値３８および平均出力値３９等、実際のタービン出力値データ４０に含まれている最大出力値４１および平均出力値４２等を有している。
【００４３】
自ノードとなるベイズノード群（４３〜４５）は、それぞれのタービンの発電量実況と予測４３と、極端気象の発生確率４４と、季節ごとの風速安定度４５を備えている。
【００４４】
ここで、極端気象の発生確率４４は、ゲリラ豪雨や突風や竜巻や台風など、極端な気象の発生確率で構成されている。これにより、風力発電が適切に行えない状態（例えば突風によって風力発電ができない状態など）を予測することができる。
【００４５】
また、季節ごとの風速安定度４５は、予め定めた単位時間あたりの風速変化の最大値が小さいほど風速が安定しているとする等、適宜の方法によって算出することができる。
【００４６】
子ノードとなる出力ノード群（４６〜５１）は、ウインドファーム全体の発電量のナウキャスト（リアルタイム値）４６と、２日前に予測される風車１基毎とウインドファーム全体の発電量４７と、１日前に予測される風車１基毎とウインドファーム全体の発電量４８と、１２時間前に予測される風車１基毎とウインドファーム全体の発電量４９と、３０分前に予測される風車１基毎とウインドファーム全体の発電量５０と、ウインドファームが複数ある場合は各ファーム毎で得られると予測される風力発電量５１とを有している。
【００４７】
消費電力予測システムとして機能するベイジアンネットワーク３０Ｂは、子がない親ノードとなる入力ノード群（６０〜６７）と、子ノードとなる出力ノード群（７３〜７７）と、自ノードとなるベイズノード群（６８〜７２）とを備えている。
【００４８】
親ノードとなる入力ノード群（６０〜６７）は、都市の気象状態データ６０に含まれている気温・湿度６１と雨・雪６２等、都市の気象モデル予測データ６３に含まれている気温・湿度６４と雨・雪６５等、スマートメータ値データに含まれているリアルタイム値６７等を有している。
【００４９】
自ノードとなるベイズノード群（６８〜７２）は、家庭での消費電気量の推定（モデル１）：大邸宅６８と、家庭での消費電気量の推定（モデル２）：中間層６９と、家庭での消費電気量の推定（モデル３）：貧困層７０と、工場７１と、ビル７２等を有している。
【００５０】
子ノードとなる出力ノード群（７３〜７７）は、住宅地・工場・ビル等で消費されている電力のナウキャスト（リアルタイム値）７３と、２日前に予測される住宅地・工場・ビル等で消費されている電力量７４と、１日前に予測される住宅地・工場・ビル等で消費されている電力量７５と、１２時間前に予測される住宅地・工場・ビル等で消費されている電力量７６と、３０分前に予測される住宅地・工場・ビル等で消費されている電力量７７とを有している。
【００５１】
ベイジアンネットワーク３０（３０Ａ，３０Ｂ）を構成するこれらの各ノード（３１〜７７）は、各ノードがとり得る状３態に応じた事前確率表又は条件付確率表（ＣＰＴ：Ｃｏｎｄｉｔｉｏｎａｌ ｐｒｏｂａｂｉｌｉｔｙ ｔａｂｌｅ）を保持しており、この確率表は、ノードがとる状態が増えるごとに（例えば、観測結果が蓄積される毎に）に更新される。従って、データ蓄積によってベイジアンネットワーク３０が成長する。
【００５２】
図４は、ベイジアンネットワーク５４の構成をデータのインプット及びアウトプットと演算式の観点から表記したブロック図である。  
データ入力部５３は、ライダー観測による風向風速データＡＦ、数値予測モデルによる予報ＡＬ、タービン出力特性ＡＲ、実際のタービン出力ＡＰ、及び都市の気象状態ＡＧ等を取得し、ベイジアンネットワーク５４に入力する。このデータ入力部５３により取得するデータは、観測値の全集合をＡ（親の無いノード）としている。
【００５３】
ここで、ライダー観測による風向風速データＡＦとして入力されるデータには、風向および風速のデータに加えて、上昇気流発生の原因となる風の収束に関する情報が含まれていることが好ましい。この風の収束に関する情報を加えておくことで、ゲリラ豪雨や突風等の突発気象を予測できる。
【００５４】
また、各データは、地点別のデータであり、ＸＹＺの関数としてあらわされる、このため、例えばライダー観測による風向風速データＡＦであれば、その１つの要素ａＦ１は、ａ（ＸＦ１，ＹＦ１，ＺＦ１）とＸ，Ｙ，Ｚの関数で表される。なお、地点についてＸＹＺの３次元座標としているが、高さを固定する場合であれば、Ｚ座標を固定あるいは省略してＺＹの２次元座標として扱ってもよい。
【００５５】
発生確率出力部５５は、ベイジアンネットワーク５４から算出される各ウインドファームの出力ナウキャストＢＨ、発電量予測ＢＢ、および各地区毎の使用電力量ＢＷ等の事象の発生確率を出力する。この発生確率出力部５５により出力するデータは、推定すべき量の全集合をＢ（子の無いノード）としている。
【００５６】
なお、風力発電装置５（図１参照）を海洋上に配置し、都市部の発電量を予測する場合、データ入力部５３の各要素については海洋上の地点を示す座標とし、発生確率出力部５５の各要素については都市部での地点を示す座標とすることが好ましい。これにより、発電場所と電力使用場所が異なっていても、適切に発電量および使用電力量を予測することができる。
【００５７】
ベイジアンネットワーク５４は、データ入力部５３から入力された各種データ（親のないノード）から、発生確率出力部５５で発生確率を出力する対象である事象（子のないノード）に状態移行する確率を示すベイズノード群で構成されており、事後確率Ｐ（Ｙ｜Ｘ）を推定する。  
このベイジアンネットワーク５４は、次に説明する条件付確率表（ＣＰＴ）を用いる。
【００５８】
図５は、条件付確率表（ＣＰＴ）の構成を示す説明図である。  
各ノードがとり得る状態はノード毎に適宜選択される要素であり、親ノード群ＲがＲｍ個の状態を取り得るとき、親ノードの先にある自ノードはＮｃ個の離散的状態をとる。従って、自ノードの条件的確率表は、図５に図示するように、Ｎｃ行、ＮＲ１＊ＮＲ２・・・＊ＮＲｍ列の行列で表現される（＊は乗算を示す）。  
この条件付確率表では、親ノード状態が列に配置されるため、例えばライダー観測による風向風速データが観測地点毎に左から右へ並んでいき、次に数値予測モデルによる予報が予報地点毎に左から右へ並んでいくというように、入力データ５３が列方向に並んでいく。そして、行としては、自ノードが並んでいく。この自ノードは、自ノードから別の自ノードへ状態遷移することもあれば、自ノードから出力ノードへ状態遷移することもある。
この条件的確率表のそれぞれの行列要素には、条件的確率Ｐ（Ｓ＝ｓｉ｜Ｒ＝ｒｊ）が入る。  
このように構成された条件付確率表は、ベイズノードが増加する毎に行方向に要素が増加していく。
【００５９】
図６は、条件付確立表（ＣＰＴ）を用いたベイズフィルタ推定を説明する説明図である。図６（Ａ）は、ベイズフィルタ推定の流れを示すフローチャートであり、図６（Ｂ）は、各ノードの状態推定法におけるベクトルを説明する説明図である。
【００６０】
図６（Ａ）に示すように、風力発電予測システム１Ａは、観測値を取得し、制御部１１（図１参照）により風向、風速、および風の安定性を算出するなど推定すべき量の結果データベースを作成しておく（ステップＳ４０１）。
【００６１】
風力発電予測システム１Ａは、制御部１１（図１参照）により、ステップＳ４０１で新しく得られた観測値集合（Ｘ）をノードにセットする（ステップＳ４０２）。ここで、親ノードも観測値も持たないノードがでたときには、事前確率を代入しておく。また、下端のノードにも主観確率（主観によって予め定めた確率）を入れておく。
【００６２】
風力発電予測システム１Ａは、制御部１１（図１参照）により、各ノードの状態変数を順に得ていく（ステップＳ４０３）。  
風力発電予測システム１Ａは、制御部１１（図１参照）により、知りたい対象Ｙの各要素の事後確率Ｐ（ｙ＿ｉ｜Ｘ）を導出する（ステップＳ４０４）。  
風力発電予測システム１Ａは、制御部１１（図１参照）により、条件付確立表（ＣＰＴ）をアップデートして学習させる。ここでノード構造を変更する場合、全データベースを用いて条件付確立表（ＣＰＴ）を全て作り直す（ステップＳ４０５）。
【００６３】
また、風力発電予測システム１Ａは、制御部１１（図１参照）により、Ｙについての事前分布から、よく一致する確率分布（ディレクトリ分布など）を得ておく（ステップＳ４０６）。
【００６４】
風力発電予測システム１Ａは、制御部１１（図１参照）により、前記ステップＳ４０４で得られたＰ（ｙｉ｜Ｘ）を取得した後、すべてのｙｉの値について確率を導出し、確率分布関数に最小二乗法にてフィッティングさせる（Ｓ４０７）。
【００６５】
風力発電予測システム１Ａは、制御部１１（図１参照）により、再尤値を求めるべき確率として出力する（ステップＳ４０８）。
【００６６】
以上のフローチャートにより、条件付確立表（ＣＰＴ）を常時更新することができるとともに、再尤値を求める確率として出力できる。
【００６７】
なお、データ量が多くなるに従って各ノードの状態数を増やすことが可能である。また、観測が欠測であっても統計的に得られている事前確率分布を使って推定することができる。
【００６８】
図６（Ｂ）は、各ノードの状態推定法におけるベクトルを説明する説明図である。
【００６９】
上流側の観測によって得たＫ個のデータ（風向分布３２、風速分布３３、風向・風速３５、気温・湿度３６、最大出力値３８、平均出力値３９、最大出力値４１、および平均出力値４２等）により、上流観測ベクトルは、Ｘｕ＝［ｘｑ１，ｘｑ２，…，ｘｑｋ］として与えられる。  
上流状態ベクトルは、Ｒ＝［ｒ１，ｒ２，…，ｒＭ］として与えられる。この上流状態ベクトルは、データ入力部５３（図４参照）の各要素に対応している。  
Ｓの状態数は、Ｎとする。このＳは、ベイズノード群７８（図４参照）に対応している。  
下流状態ベクトルは、Ｔ＝［ｔ１，ｔ２，…，ｔＬ］として与えられる。この下流状態ベクトルは、発生確率出力部５５（図４参照）の各要素に対応している。  
下流側の観測によって得られるＬ個のデータ（住宅地・工場・ビル等で消費されている電力のナウキャスト（リアルタイム値）７３、２日前に予測される住宅地・工場・ビル等で消費されている電力量７４、１日前に予測される住宅地・工場・ビル等で消費されている電力量７５、１２時間前に予測される住宅地・工場・ビル等で消費されている電力量７６、３０分前に予測される住宅地・工場・ビル等で消費されている電力量７７）として、下流観測ベクトルは、Ｘｄ＝［ｘ１，ｘ２，…，ｘｕＬ］として与えられる。
【００７０】
ここで、上流のＸｕと下流のＸｄは条件的独立であるから、次の式（数Ａ）が成立する。  
［数Ａ］
Ｐ（Ｓ＝Ｓｎ）  
＝Ａ・Ｐ（Ｘｕ｜Ｓ＝Ｓｎ）Ｐ（Ｓ＝Ｓｎ｜Ｘｄ）  
＝Ａ・ωｕ（Ｓ＝Ｓｎ）ωｄ（Ｓ＝Ｓｎ）  
※Ａは尤度を確率値に正規化する定数を示す。  
ωｕ（Ｓｎ）、ωｄ（Ｓｎ）は上流と下流の観測のＳへの寄与分を示す。
【００７１】
そして、ωｕ（Ｓ＝Ｓｎ）、ωｄ（Ｓ＝Ｓｎ）は、次の式（数１）,式（数２）によって計算することができる。
【００７２】
【数１】
2017102760.tif 000003
【００７３】
【数２】
2017102760.tif 000004
この解き方として、観測値が与えられるノードはその値を導入する。最上流で観測値が与えられないノードには、事前確率を導入する。最下流で観測値が与えられないノードには、無情報として一様分布を導入する。さらに、上流や下流のノードがある場合は、再帰的に上式を適用してゆくことでネットワークのすべてのノードの状態確率が得られる。
【００７４】
なお、矢印の向きを考えない場合にループがないネットワークは、上記の方法のみベイジングネットワークを解くことができるが、今回の実装では矢印の方向を考慮しない場合にループが発生する。従って、例えばＬｏｏｐｙＢＰ法等の手法を導入し、複結合ネットワークに対応する。
【００７５】
図７は、風力発電装置５の風車を設置するのに適切な設置場所を選定する風車設置場所選定処理を実行する制御部１１（図１参照）の動作を示すフローチャートである。
【００７６】
制御部１１は、ドップラーライダー２（図１参照）で観測できる範囲において数１０〜１００ｍ程の分解能で風速の値を得る（ステップＳ５０１）。このとき得る値は、例えば格子点毎の格子点値等とすることができる。
【００７７】
制御部１１は、ドップラーライダー２で観測できる範囲において数１０〜１００ｍ程の分解能で風速バリアンス値を得る（ステップＳ５０２）。このとき得る値は、例えば格子点毎の格子点値等とすることができる。
【００７８】
制御部１１は、取得した風速値とバリアンス値のデータセットを全点探索する（ステップＳ５０３）。
【００７９】
制御部１１は、探索結果を風速値とバリアンス値のグラフにする（ステップＳ５０４）。このとき、例えば縦軸を風速値、横軸をバリアンス値とする二次元グラフを作成し、この座標上に各格子点をプロットすることでグラフ化することができる。
【００８０】
制御部１１は、風速値が大きくバリアンス値の小さい点があるか否か判定するＵステップＳ５０５）。この判定は、例えば風速値が予め定めた所定値（風速基準値）より大きく、かつ、バリアンス値が予め定めた所定値（バリアンス基準値）より小さいものが存在するか否かとするなど、適宜の判定方法とすることができる。
【００８１】
制御部１１は、風速値が大きくバリアンス値の小さい点があれば（ステップＳ５０５：Ｙｅｓ）、この点を風車設置場所として選定する（ステップＳ５０６）。この選定の結果は、風車設置場所を示す座標（緯度、経度、高度）で出力するなど、場所を特定する情報として出力する。なお、風速値が大きくバリアンス値の小さい点が複数ある場合は、判定の基準値に対する風速値の差とバリアンス値の差を加算した値が最も大きい点を採用するなど、適宜の条件によって最も効率の良い点に定めると良い。
【００８２】
制御部１１は、風速値が大きくバリアンス値の小さい点がなければ（ステップＳ５０５：Ｎｏ）、設置場所候補なしとして（ステップＳ５０７）、処理を終了する。
【００８３】
以上の構成および動作により、気象条件の変化によって発電量が変化する風力発電において、風力発電装置の発電量を精度良く予測することができる。
【００８４】
また、極端気象を予測して風力発電の予測に役立てるため、風力発電量をより精度よく予測することができる。すなわち、強い突風が断続的に吹いた場合、従来の方法で風速の合計から発電量を算出すると、実際に発電される発電量よりも多く算出される。これは、突風等では風力発電装置が適切に作動せず、発電量が低くなること等が原因である。このような不適切な予測を極力排除し、精度を向上させることができる。
【００８５】
また、地点別に風力や風向を検出し、さらに地点別に風力や風向を予測できるため、風力発電装置の風車をどこに設置すれば最も効率よく発電できるか算出することができる。
【００８６】
詳述すると、地表上における風は、様々な要素が複雑にからみあって発生している。このため、一見すれば広い海洋上であっても、位置が少し異なれば、年間を通じて発生する風速や風向が全く異なってくる。そうすると、全く同じ風力発電装置であっても、設置する位置が少し異なるだけで得られる発電力が全く異なるということにもなる。
【００８７】
また、ベイジアンネットワーク５４は、ベイズノード群７８を適宜増加させて強化学習していくことができる。これにより、ベイズノードを増加させて、さらに精度を高めることができる。
【００８８】
このようにして、図８のシステム構成図に示すように、都市の気象状態データ６０、都市の気象モデル予測データ６３、スマートメータ値６６、実際のタービンの出力値４０等のデータをベイジアンネットワーク５４に入力し、全体の電力量や、風車毎の電力量、時間別の発電量の予想等を得ることができる。
【００８９】
また、図９のシステム構成図に示すように、発電会社の発電量、電子卸市場の売買情報、顧客の希望情報等を新電力会社にて算出し、このバランスをとって電力卸売市場から不足する電力を購入し、あるいは余剰の電力を販売するといったことができる。また、再生可能エネルギーの発電量を精度よく予測でき、需要と供給のバランスが常に求められている新電力会社に発電量・需要量予測情報を提供することができて、３０分市場（即時市場）に対応することができる。また、前もって再生可能エネルギーの情報を提供することで、インバランスを最小に抑えることができる。
【００９０】
また、図１０のシステム構成図に示すように、電力を卸売りする卸売市場や、発電量の実況や予測を行う電力会社は、いつ電力を買うのが良く、いつ電力を売るのが良いか悩むことがある。これに対して、発電量の通知、購入した電力量等を知らせ、電力を買うタイミングの検討等に時間をとられることを防止でき、エラーバーを３％以内に収めるといったことができる。
【００９１】
なお、この発明は、上述した実施形態に限られず、他の様々な実施形態とすることができる。
【産業上の利用可能性】
【００９２】
この発明は、風の状況を予測して利用するような産業に利用することができる。
【符号の説明】
【００９３】
１…予測システム
２…ドップラーライダー
１１…制御部
１３…通信部
３１…ドップラーライダー観測データ
３４…予測データ
【図１】
2017102760.tif 000005
【図２】
2017102760.tif 000006
【図３】
2017102760.tif 000007
【図４】
2017102760.tif 000008
【図５】
2017102760.tif 000009
【図６】
2017102760.tif 000010
【図７】
2017102760.tif 000011
【図８】
2017102760.tif 000012
【図９】
2017102760.tif 000013
【図１０】
2017102760.tif 000014