(57)【要約】【課題】  風車に衝突する羽数を定量的に評価する。【解決手段】  風力発電設備の設置域で、風車回転高度幅の容積を算出する（Ｓ２）。風車のブレード回転面を中心とする球体を風車接触可能域と想定し、全風車の接触可能域の容積を算出する（Ｓ３）。後者の容積を前者の容積で除算することで、鳥類に依存しない空間占有率を算出する（Ｓ４）。個々の鳥類について（Ｓ５）、風車接触可能域を通過する通過頻度を算出し（Ｓ６）、その内でブレード回転面を通過する回数を算出し（Ｓ７）、回避率及び接触率を加味して、衝突羽数を算出する（Ｓ８）。【選択図】  図２
2014037795.tif 000002
【特許請求の範囲】
【請求項１】
風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価装置であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出手段と、
注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出し、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出手段
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価装置。
【請求項２】
当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１に記載の鳥類衝突数評価装置。
【請求項３】
当該衝突羽数算出手段は、
当該注目する鳥類が当該風車接触可能域を通過する頻度を算出する手段と、
当該頻度と、当該ブレード回転面への進入を回避する回避率とから当該ブレード回転面を通過する回数に相当する値を算出する手段と、
当該ブレード回転面で当該ブレードに接触する接触率を算出する手段
とを具備することを特徴とする請求項１又は２に記載の鳥類衝突数評価装置。
【請求項４】
当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項３に記載の鳥類衝突数評価装置。
【請求項５】
当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項３に記載の鳥類衝突数評価装置。
【請求項６】
更に、当該鳥類とは別の鳥類について、当該衝突羽数算出手段に当該衝突羽数を算出させる制御手段を具備することを特徴とする請求項１乃至５の何れか１項に記載の鳥類衝突数評価装置。
【請求項７】
下記式に従い、風力発電設備の風車に衝突する鳥類の羽数を定量評価する鳥類衝突数評価装置。
衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（１−回避率）
ただし、
Ｃｎ＝（Ｖｓ／Ｖａ）×Ｍｄ／（ｄａｖｇ×２）
Ｖａ＝Ｓａ×Ｈｃ
Ｖｓ＝Ｍ×（４π／３）Ｒａ３
Ｓａ：風力発電設備の設置域の面積
Ｈｃ：風車回転高度幅
Ｍ：風車台数
Ｒａ：風車のブレード半径又は風車軸からブレード端までの距離
Ｍｄ：注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
ｄａｖｇ：風車接触可能域を通過する平均距離
接触率：ブレード回転面に進入した個体がブレードに接触する確率
回避率：風車設置域、風車接触可能域及びブレード回転面への進入を回避する合計の割合
【請求項８】
当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項７に記載の鳥類衝突数評価装置。
【請求項９】
当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項７に記載の鳥類衝突数評価装置。
【請求項１０】
風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出ステップと、
注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出し、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出ステップ
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価方法。
【請求項１１】
当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１０に記載の鳥類衝突数評価方法。
【請求項１２】
当該衝突羽数算出ステップは、
当該注目する鳥類が当該風車接触可能域を通過する頻度を算出するステップと、
当該頻度と、当該ブレード回転面への進入を回避する回避率とから当該ブレード回転面を通過する回数に相当する値を算出するステップと、
当該ブレード回転面で当該ブレードに接触する接触率を算出するステップ
とを具備することを特徴とする請求項１０又は１１に記載の鳥類衝突数評価方法。
【請求項１３】
当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項１２に記載の鳥類衝突数評価方法。
【請求項１４】
当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項１２に記載の鳥類衝突数評価方法。
【請求項１５】
更に、当該鳥類とは別の鳥類について、当該衝突羽数算出手段に当該衝突羽数を算出させる制御ステップを具備することを特徴とする請求項１０乃至１４の何れか１項に記載の鳥類衝突数評価方法。
【請求項１６】
風力発電設備の風車への鳥類衝突数を定量評価するコンピュータプログラムであって、コンピュータに、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出させる空間占有率算出機能と、
注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出させ、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出させる衝突羽数算出機能
とを実現させ、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価プログラム。
【請求項１７】
当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１６に記載の鳥類衝突数評価プログラム。
【請求項１８】
当該衝突羽数算出機能は、
当該コンピュータに当該注目する鳥類が当該風車接触可能域を通過する頻度を算出させる機能と、
当該コンピュータに、当該頻度と、当該ブレード回転面への進入を回避する回避率とから当該ブレード回転面を通過する回数に相当する値を算出させる機能と、
当該コンピュータに、当該ブレード回転面で当該ブレードに接触する接触率を算出させる機能
とを具備することを特徴とする請求項１６又は１７に記載の鳥類衝突数評価プログラム。
【請求項１９】
当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項１８に記載の鳥類衝突数評価プログラム。
【請求項２０】
当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項１８に記載の鳥類衝突数評価プログラム。
【請求項２１】
当該鳥類とは別の鳥類について、当該衝突羽数算出機能に当該衝突羽数を算出させる制御機能を当該コンピュータに実現させることを特徴とする請求項１６乃至２０の何れか１項に記載の鳥類衝突数評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、風力発電の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法、システム及びプログラムに関する。
【背景技術】
【０００２】
風力発電では、しばしば周囲を飛翔する鳥類が風車に衝突する、いわゆるバードストライクを引き起こす。これは、風車を傷めて稼働率を低下させるばかりでなく、希少鳥類が周囲に生息する場合に、その生存を脅かすことになる。
【０００３】
そこで、風力発電事業のアセスメントとして、風車に衝突する鳥類衝突数を定量評価することが求められる。環境省は、ホームページ上に、平成２３年１月７日付けで「鳥類等に関する風力発電施設立地適正化のための手引き」を公開している（非特許文献１）。
【０００４】
同手引きには、鳥類衝突数定量評価モデルとして、鳥類の飛翔高度と風車のブレードの回転域の高さとの関係から衝突数を定量評価する手法が記載されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】ｈｔｔｐ：／／ｗｗｗ．ｅｎｖ．ｇｏ．ｊｐ／ｐｒｅｓｓ．ｐｈｐ？ｓｅｒｉａｌ＝１３３３１
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従来の鳥類衝突数定量評価モデルでは、風車ブレードの回転面の高度を飛翔し、回転面に進入する鳥類数を計数している。
【０００７】
しかし、実際には、回転面に対して斜めに進入する鳥類もあり、これらは無視できない。
【０００８】
本発明は、このような斜め進入も考慮して鳥類衝突数を定量評価する鳥類衝突数評価装置、方法及びプログラムを提示することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る鳥類衝突数評価装置は、風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価装置であって、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出手段と、注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出し、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出手段とを具備し、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。
【００１０】
本発明に係る鳥類衝突数評価装置は、下記式に従い、風力発電設備の風車に衝突する鳥類の羽数を定量評価する鳥類衝突数評価装置であって、
衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（１−回避率）
である。ただし、
Ｃｎ＝（Ｖｓ／Ｖａ）×Ｍｄ／（ｄａｖｇ×２）
Ｖａ＝Ｓａ×Ｈｃ
Ｖｓ＝Ｍ×（４π／３）Ｒａ３
Ｓａ：風力発電設備の設置域の面積
Ｈｃ：風車回転高度幅
Ｍ：風車台数
Ｒａ：風車のブレード半径又は風車軸からブレード端までの距離
Ｍｄ：注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
ｄａｖｇ：風車接触可能域を通過する平均距離
接触率：ブレード回転面に進入した個体がブレードに接触する確率
回避率：風車設置域、風車接触可能域及びブレード回転面への進入を回避する合計の割合。
【００１１】
本発明に係る鳥類衝突数評価方法は、風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法であって、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出ステップと、注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出し、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出ステップとを具備し、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。
【００１２】
本発明に係る鳥類衝突数評価プログラムは、風力発電設備の風車への鳥類衝突数を定量評価するコンピュータプログラムであって、コンピュータに、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出させる空間占有率算出機能と、注目する鳥類について、当該風車接触可能域に進入して当該風車のブレード回転面でブレードに接触する接触数を算出させ、当該接触数に当該空間占有率及び当該風車の稼働率を乗算して衝突羽数を算出させる衝突羽数算出機能とを実現させ、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、個々の風車に対して球体モデルによる風車接触可能域を設定したことで、全球方向から風車に向かう個体の衝突確率を適切に定量評価することが可能になる。これにより、衝突羽数をより精確に見込むことが可能になる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施例の概略構成ブロック図である。
【図２】本実施例の動作フローチャートである。
【図３】風力発電設備情報テーブルの構成例である。
【図４】鳥類情報テーブルの構成例である。
【図５】風車のパラメータを示す正面図である。
【図６】風車接触可能域への進入確率を説明する、風車を上方向から見た平面図である。
【図７】接触確率を説明する模式図である。
【図８】進入角度対接触確率のグラフ例である。
【発明を実施するための形態】
【００１５】
以下、図面を参照して、本発明の実施例を詳細に説明する。
【実施例１】
【００１６】
図１は、本発明の一実施例の概略構成ブロック図を示し、図２は、本実施例による鳥類衝突数評価の動作フローチャートを示す。本実施例では、コンピュータ上のソフトウエアにより実現されうるが、もちろん、そのようなソフトウエアで実現される機能を組み込んだ専用装置によっても実現可能である。
【００１７】
本実施例では、風力発電設備の風車に対して全球方向から鳥類が進入する可能性が一般的に存在することを考慮し、個々の風車に対して衝突可能性のある鳥類の進入エリアとして風車ブレードの回転に外接する球体を風車接触可能域として想定する。
【００１８】
コンピュータ１０は、ＣＰＵ１２、主記憶としてのメモリ１４、補助記憶としてのＨＤＤ１６、表示装置１８、キーボード２０及び、ポインティングデバイスとしてのマウス２２を具備し、鳥類衝突数評価装置の一実施例として機能する。図２に示すフローチャートに相当する鳥類衝突数評価のコンピュータプログラムがＨＤＤ１６に格納されており、ＣＰＵ１２は、その鳥類衝突数評価プログラムをＨＤＤ１６からメモリ１４に読み出し、実行する。これにより、ＣＰＵ１２は、図２に示すフローチャートで説明される各機能を実現する。
【００１９】
コンピュータ１０には、予め、設置しようとする風力発電設備情報と、風力発電設備を設置しようとする地域に生息し、衝突可能性のある鳥類の情報が格納される。ＨＤＤ１６の風力発電設備情報テーブル１６ａには風力発電設備情報が格納され、鳥類情報テーブル１６ｂには鳥類情報が格納される。
【００２０】
図３は、風力発電設備情報テーブル１６ａの構造例を示す。設置面積（Ｓａ）は、風力発電設備情報を設置しようとする地域の面積であって、例えば、設置する風車に衝突する可能性のある鳥類の生息域をカバーする範囲で設定される。明確さの点で、設置面積（Ｓａ）を、設置する複数の風力発電設備を外接する多角形その他の線図で定義しても良い。稼働率は１年で風車が回転する時間の割合を示す。風車軸高度（Ｈａｘｉｓ）は、風車の回転軸の、地上からの高度を示す。ブレード半径（Ｒａ）は、風車を構成する羽（ブレード）の半径、即ち、風車軸からブレードの端までの距離を示す。平均回転速度は、想定される風速に対する平均的な回転速度を示す。平均回転速度は、風車を通過する鳥類がブレードに接触する接触確率を計算するのに使用される。図５は、風車を正面から見た状態の模式図であって、風車軸高度（Ｈａｘｉｓ）、ブレード半径（Ｒａ）及び風車回転高度幅（２Ｒａ）の関係を示す。
【００２１】
図４は、鳥類情報テーブル１６ｂのデータ構造例を示す。鳥類情報テーブル１６ｂは、衝突可能性のある鳥類毎に、サイズ（体長Ｉ，翼開長Ｗ）、飛翔高度範囲（Ｈａ〜Ｈｂ）、飛翔速度（ｖ）、飛翔距離（Ｍｄ）及び回避率（ＰＬ）の各データからなり、専ら実測により決定される。サイズは、羽を拡げた飛翔時の飛翔方向の大きさ（体長：Ｉ）と羽方向の大きさ（翼開長Ｗ）を示す。飛翔高度範囲（Ｈａ〜Ｈｂ）は、通常、飛翔する高度の範囲を示す。飛翔速度（ｖ）は、風車のブレードが回転する範囲の高度での飛翔速度を示す。飛翔距離（Ｍｄ）は、単位期間（例えば、１年）内に風車回転高度範囲（Ｈａ〜Ｈｂ）を飛翔する場合の設置地域内での全飛翔個体の合計距離を示す。回避率（ＰＬ）は、回転する風車に向かって飛びつつも風車を回避する確率を示し、より具体的には、風車設置域、風車接触可能域及びブレード回転面への進入を回避する合計の割合である。例えば、オジロワシの場合、回避率は９５％程度と言われている。
【００２２】
図２を参照して、本実施例による鳥類衝突数の定量評価動作を説明する。先ず、オペレータは、キーボード２０及びマウス２２を使い、設置しようとする風力発電設備情報とその地域の鳥類情報をＨＤＤ１６に格納する（Ｓ１）。これにより、ＨＤＤ１６に風力発電設備情報テーブル１６ａと鳥類情報テーブル１６ｂが格納され、準備される。勿論、遠隔地のサーバにある同様のテーブルを参照してもよいが、ここでは理解を容易にするために、ローカルにこれらテーブル１６ａ，１６ｂが存在するものとしている。
【００２３】
ＣＰＵ１２は、風力発電設備情報テーブル１６ａを参照し、風車回転高度容積（Ｖａ）を算出する（Ｓ２）。風車回転高度容積（Ｖａ）は、テーブル１６ａの設置面積（Ｓａ）に風車回転高度幅（Ｈｃ）を乗算して得られる。すなわち、
Ｖａ＝Ｓａ×Ｈｃ
とする。風車回転高度幅（Ｈｃ）は、図５に示すように、回転する風車の上端から下端までの距離であり、一般的には、ブレード半径（Ｒａ）の２倍に相当する。風車回転高度容積（Ｖａ）は、高度で見た空間的な存在可能性を示す指数となる。
【００２４】
風車回転高度幅の範囲に鳥類が飛翔していたとしても、ある程度に風車に接近しない限り、風車に衝突することはあり得ない。また、鳥類が風車に衝突するのは、鳥類が風車の正面から風車に向かって進入する場合に限らず、斜めに進入する場合もありうる。却って、風車回転面に直角に進入する場合に比べ、斜めに入射する方が、風車のブレードに接触する時間的空間的確率が高いので、風車に衝突する確率が高くなる。そこで、本実施例では、個々の風車に対して、風車のブレード回転面を中心とする球体を風車接触可能域として想定する球体モデルを採用する。
【００２５】
ＣＰＵ１２は、風車接触可能域として想定される全球体の容積を風車接触可能容積（Ｖｓ）として、
Ｖｓ＝Ｍ×（４π／３）Ｒａ３
を計算する（Ｓ３）。但し、Ｍは風車台数を示し、Ｒａは、風車のブレード半径（又は風車軸からブレード端までの距離）を示す。
【００２６】
風車回転高度容積（Ｖａ）に対する風車接触可能容積（Ｖｓ）の比は、空間的な広さだけから見た場合に風車に衝突する可能性のある風車接触可能域に鳥類が存在する確率を示す。この確率を空間占有率と定義し、ＣＰＵ１２は、下記式に従い、設置地域における空間的な空間占有率（Ｐｖ）を算出する（Ｓ４）。すなわち、
Ｐｖ＝Ｖｓ／Ｖａ
とする。
【００２７】
空間占有率（Ｐｖ）は、鳥類の種別に依存しない指数であり、本実施例は、この空間占有率（Ｐｖ）に鳥類毎の事情を加味して、個々の鳥類の衝突羽数を推計する。
【００２８】
先ず、ＣＰＵ１２は、鳥類情報テーブル１６ｂを参照し、鳥類情報テーブル１６ｂから評価対象の鳥類を選択し、選択した鳥類の緒データ（サイズ（Ｉ，Ｗ）、飛翔高度範囲（Ｈａ〜Ｈｂ）、飛翔速度（ｖ）、飛翔距離（Ｍｄ）及び回避率（ＰＬ））を読み込む（Ｓ５）。
【００２９】
ＣＰＵ１２は、選択された鳥類が風車接触可能域を通過する通過頻度（Ｆｑ）を算出する（Ｓ６）。この通過頻度（Ｆｑ）は、選択された鳥類が所定の期間内に風車接触可能域（風車接触可能容積（Ｖｓ））を通過する回数であり、選択された鳥類について、風車回転高度幅Ｈａ〜Ｈｂでの全飛翔距離（Ｍｄ）を、１台の風車の風車接触可能域を通過する平均距離（ｄａｖｇ）で除算することで得られる。すなわち、ＣＰＵ１２は、
Ｆｑ＝Ｐｖ×Ｍｄ／ｄａｖｇ
で頻度（Ｆｑ）を算出する。頻度（Ｆｑ）の単位は、鳥類種類別の所定期間内の回数となる。
【００３０】
球体をランダムに通過するとした場合の平均距離（ｄａｖｇ）は、球体にランダムに１本の線を通過させたときの、始点（入口）から終点（出口）までの距離の平均であり、通過軌跡の総和を軌跡数で除算することで得られる。従って、平均距離（ｄａｖｇ）は、球の体積、すなわち、風車接触可能域を構成する球体の体積を、その中心を通る円の面積で除算した値となる。
【００３１】
次に、ＣＰＵ１２は、何れかの風車の風車接触可能域に進入した鳥類が更に、そのブレード回転面を通過する回数（Ｃｎ）を算出する（Ｓ７）。ステップＳ６で計算した通過頻度（Ｆｑ）は、風車のブレード回転面を中心とする球体に進入する頻度を示すが、球体に進入した鳥がそのまま直進したとしても、全てがブレード回転面に進入するわけではなく、以下に説明するように、ブレード回転面に進入する確率は１／２となる。従って、ブレード回転面通過回数（Ｃｎ）は、
Ｃｎ＝Ｆｑ／２＝Ｐｖ×Ｍｄ／（ｄａｖｇ×２）
となる。
【００３２】
図６を参照して、ブレード回転面に進入する確率が１／２であることを説明する。図６は、風車をその上空から見た平面図を示す。本実施例の球体モデルでは、ブレード回転面５０に対して、外接する球体５２を想定する。この球体５２の内部が風車接触可能域に対応する。この球体５２に外部から任意の角度で鳥が線５４で示すように進入するとする。入射点とブレード回転面の両端を結ぶ線のなす角は周知の通り９０度となり、球体５２の球面への進入点での接線で見ると、余角、すなわちブレード面方向に進まない角度も９０度である。これは、進入する個体がブレード面を見込む立体角と、ブレード面以外を見込む立体角の比で見ても、同様である。従って、先に説明したように、風車接触可能域を構成する球体５２に進入した鳥類がブレード面に突入する確率は常に１／２となる。
【００３３】
ここまでは、鳥類が直進することを前提としてきた。しかし、実際には、ブレード回転面に向かって飛翔する鳥類が全てブレード回転面に進入するとは限らない。回避行動があり、回避しない個体も、ブレード回転面を無事に通過できる場合がある。回避できず、ブレード回転面を無事に通過できなかった個体が、ブレードに衝突することになる。
【００３４】
ブレード回転面に接近した鳥類も、その飛翔能力によりブレード回転面への進入を回避できる。このような回避率は、経験的又は実測により定量評価可能であり、本実施例でも、既存の方法で計測又は決定された知られた数値を採用する。例えば、オジロワシの場合で、風車設置域、風車接触可能域及びブレード回転面への進入を回避する合計の割合である回避率は、９５％程度と言われている。
【００３５】
回避せずにブレード回転面に進入しても、一部は、ブレード回転面を無事に通過できる。無事に通過できるかどうかは、鳥類のサイズ（Ｉ，Ｗ）と、飛翔速度（ｖ）、ブレード面への入射角度（θ）、風車の回転速度及びブレード数に依存する。図７は、ブレード回転面を上方から見た場合の２種類の進入角度例であって、図７（ａ）はブレード回転面への進入角度が浅い（θが小さい）場合を示し、同（ｂ）は進入角度が大きい（θが大きい）場合を示す。
【００３６】
鳥類の頭がブレード回転面に入ってから尾が抜けるまでの間にブレードが鳥類に接触する場合に、鳥類がブレードに衝突したことになる。進入角度θが小さいほどブレード回転面に位置する時間が長くなるので、接触確率が増大する。進入角度θが大きくなると、短時間でブレード回転面を通過することが可能になり、接触確率が低下する。進入角度θを０度から９０度の範囲で接触確率を算定すると、平均接触率Ｔを算定できる。すなわち、進入角度θの接触確率は、｛（ブレード回転面通過距離）／（飛翔速度ｖ）｝×｛（ブレード羽数）／（１回転の時間）｝で計算できる。平均衝突率Ｔは、各進入角度の衝突確率を進入角度で積分して平均化することで得られる。平均の代わりに、代表的な進入角度における接触確率を代表値として採用しても良い。
【００３７】
ブレード回転面通過距離は、以下のように算出できる。進入角度θで野鳥がブレード面に進入する場合で、体長Ｉの１／２の箇所で翼開長Ｗを十字クロスさせると、
通過距離＝Ｗ／ｔａｎθ
となる。但し、オジロワシ、オオワシ、イヌワシ及びチュウヒの４種で推定した通過距離は、進入角度θが約６８度以上で各種類の体長（Ｉ）と同じ値になる。ブレード厚みｂを考慮する場合、
通過距離＝Ｗ／ｔａｎθ＋ｂ／ｓｉｎθ
となる。
【００３８】
オジロワシを例に説明する。オジロワシのオスメス平均の体長及び翼開長から求めたブレード面を通過するのに要する距離、すなわちブレード面通過距離は、平均で２．３ｍ程であり、飛翔速度は１０．６ｍ／ｓである。このとき、ブレード回転面を垂直に通過する時間ｔは、
ｔ＝（ブレード面通過距離）／飛翔速度
であるから、０．２１７秒（＝２．３／１０．６）となる。
【００３９】
風車のブレードが２．５秒に１回転すると仮定すると、ブレード１枚がｔ秒間に移動する間にブレード回転面を掃引する面積率は０．２１７／２．５となる。風車のブレードが３枚である場合、合計の掃引面積率は、
３×０．２１７／２．５＝０．２６
となる。つまり、合計の掃引面積率は、
ｔ×（ブレード数）／（１回転の秒数）
で得られる。この合計掃引面積率が９０度進入の場合の接触確率となる。進入角度θを変更して同様の計算を行い、０〜９０度の角度別の接触確率を平均化することで、平均接触率が得られる。図８は、進入角度対接触確率の一例を示す。
【００４０】
以上を踏まえて、ＣＰＵ１２は、所定期間内に風車に衝突する羽数、すなわち期間内衝突羽数を、下記式に従って算定する（Ｓ８）。すなわち、
期間内衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（回転面進入率）
ただし、
回転面進入率＝１−回避率
である。
【００４１】
別種の野鳥についても衝突羽数を算定すべき場合（Ｓ９）、ＣＰＵ１２は、別の鳥類を選択して、緒データを鳥類情報テーブル１６ｂから読み込み（Ｓ１０）、当該鳥類について、ステップＳ６〜Ｓ８を実行する。
【００４２】
全ての野鳥について衝突羽数を算出したら（Ｓ９）、ＣＰＵ１２は、所定書式のレポートを生成し、ＨＤＤ１６に保存して（Ｓ１１）、終了する。
【００４３】
オペレータは、算出結果に疑問がある場合には、テーブル１６ａ，１６ｂのデータを修正して図２に示すフローを再実行すればよい。
【００４４】
上記実施例では、ステップＳ６で通過頻度の計算にＳ４で得られた空間占有率を乗算しているが、ステップＳ６、Ｓ７では割合のみを計算し、最終的にステップＳ８において空間占有率を乗算するようにしてもよいことは、明らかである。
【００４５】
このように、本実施例によれば、個々の風車に対して球体モデルによる風車接触可能域を設定したことで、全球方向から風車に向かう個体の衝突確率を適切に定量評価することが可能になり、従って、衝突羽数をより精確に見込むことが可能になる。
【００４６】
特定の説明用の実施例を参照して本発明を説明したが、特許請求の範囲に規定される本発明の技術的範囲を逸脱しないで、上述の実施例に種々の変更・修整を施しうることは、本発明の属する分野の技術者にとって自明であり、このような変更・修整も本発明の技術的範囲に含まれる。
【符号の説明】
【００４７】
１０：コンピュータ
１２：ＣＰＵ
１４：メモリ
１６：ＨＤＤ１
１６ａ：風力発電設備情報テーブル
１６ｂ：鳥類情報テーブル
１８：表示装置
２０：キーボード
２２：マウス
５０：ブレード回転面
５２：風車接触可能域（球体）
５４：進入方向
【図１】
2014037795.tif 000003
【図２】
2014037795.tif 000004
【図３】
2014037795.tif 000005
【図４】
2014037795.tif 000006
【図５】
2014037795.tif 000007
【図６】
2014037795.tif 000008
【図７】
2014037795.tif 000009
【図８】
2014037795.tif 000010