特許 6016211 鳥類衝突数評価装置、方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6016211

(24)【登録日】2016年10月7日

(45)【発行日】2016年10月26日

(54)【発明の名称】鳥類衝突数評価装置、方法及びプログラム

(51)【国際特許分類】

   F03D 80/00 20160101AFI20161013BHJP

   F03D 1/06 20060101ALI20161013BHJP

【ＦＩ】

   F03D80/00

   F03D1/06 A

【請求項の数】21

【全頁数】14

(21)【出願番号】特願2012-180105(P2012-180105)

(22)【出願日】2012年8月15日

(65)【公開番号】特開2014-37795(P2014-37795A)

(43)【公開日】2014年2月27日

【審査請求日】2015年4月9日

(73)【特許権者】

【識別番号】397039919

【氏名又は名称】一般財団法人日本気象協会

(73)【特許権者】

【識別番号】512212807

【氏名又は名称】由井 正敏

(74)【代理人】

【識別番号】100090284

【弁理士】

【氏名又は名称】田中 常雄

(72)【発明者】

【氏名】由井 正敏

(72)【発明者】

【氏名】島田 泰夫

【審査官】 山本 崇昭

(56)【参考文献】

【文献】 特開２０１０−０７１１００（ＪＰ，Ａ）

【文献】 特表２０１２−５１２３５２（ＪＰ，Ａ）

【文献】 米国特許出願公開第２００８／０２６０５３１（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｆ０３Ｄ １／００−８０／８０

(57)【特許請求の範囲】

【請求項1】

風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価装置であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出手段と、
注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力手段と、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出手段と、
当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出手段と、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出手段
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価装置。

【請求項2】

当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１に記載の鳥類衝突数評価装置。

【請求項3】

さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する手段と、
当該回避率を１から減算した値を当該回転面進入率とする演算手段
とを具備することを特徴とする請求項１または２に記載の鳥類衝突数評価装置。

【請求項4】

当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項１から３の何れか１項に記載の鳥類衝突数評価装置。

【請求項5】

当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項１から３の何れか１項に記載の鳥類衝突数評価装置。

【請求項6】

更に、当該鳥類とは別の鳥類について、当該通過頻度算出手段、当該接触率算出手段及び当該衝突羽数算出手段を使って当該衝突羽数を算出させる制御手段を具備することを特徴とする請求項１乃至５の何れか１項に記載の鳥類衝突数評価装置。

【請求項7】

下記式に従い、風力発電設備の風車に衝突する鳥類の羽数を定量評価する鳥類衝突数評価装置。
衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（１−回避率）
ただし、
Ｃｎ＝（Ｖｓ／Ｖａ）×Ｍｄ／（ｄ_ａｖｇ×２）
Ｖａ＝Ｓａ×Ｈｃ
Ｖｓ＝Ｍ×（４π／３）Ｒａ^３
Ｓａ：風力発電設備の設置域の面積
Ｈｃ：風車回転高度幅
Ｍ：風車台数
Ｒａ：風車のブレード半径又は風車軸からブレード端までの距離
Ｍｄ：注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
ｄ_ａｖｇ：風車接触可能域を通過する平均距離
接触率：ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率
回避率：当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該ブレード回転面への進入を回避する個体数の割合

【請求項8】

当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項７に記載の鳥類衝突数評価装置。

【請求項9】

当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項７に記載の鳥類衝突数評価装置。

【請求項10】

風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出ステップと、
注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力ステップと、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出ステップと、
当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出ステップと、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出ステップ
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価方法。

【請求項11】

当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１０に記載の鳥類衝突数評価方法。

【請求項12】

さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する入力ステップと、
当該回避率を１から減算した値を当該回転面進入率とする演算ステップ
とを有することを特徴とする請求項１０または１１に記載の鳥類衝突数評価方法。

【請求項13】

当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項１０から１２の何れか１項に記載の鳥類衝突数評価方法。

【請求項14】

当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項１０から１２の何れか１項に記載の鳥類衝突数評価方法。

【請求項15】

更に、当該鳥類とは別の鳥類について、当該通過頻度算出ステップ、当該接触率算出ステップ及び当該衝突羽数算出ステップに当該衝突羽数を算出させる制御ステップを具備することを特徴とする請求項１０乃至１４の何れか１項に記載の鳥類衝突数評価方法。

【請求項16】

風力発電設備の風車への鳥類衝突数を定量評価するコンピュータプログラムであって、コンピュータに、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出機能と、
注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力機能と、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出機能と、
当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出機能と、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出機能
とを実現させ、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価プログラム。

【請求項17】

当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項１６に記載の鳥類衝突数評価プログラム。

【請求項18】

さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する入力機能と、
当該回避率を１から減算した値を当該回転面進入率とする演算機能
とを当該コンピュータに実現させることを特徴とする請求項１６又は１７に記載の鳥類衝突数評価プログラム。

【請求項19】

当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項１６から１８の何れか１項に記載の鳥類衝突数評価プログラム。

【請求項20】

当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項１６から１８の何れか１項に記載の鳥類衝突数評価プログラム。

【請求項21】

当該鳥類とは別の鳥類について、当該通過頻度算出機能、当該接触率算出機能及び当該衝突羽数算出機能に当該衝突羽数を算出させる制御機能を当該コンピュータに実現させることを特徴とする請求項１６乃至２０の何れか１項に記載の鳥類衝突数評価プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、風力発電の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法、システム及びプログラムに関する。

【背景技術】

【0002】

風力発電では、しばしば周囲を飛翔する鳥類が風車に衝突する、いわゆるバードストライクを引き起こす。これは、風車を傷めて稼働率を低下させるばかりでなく、希少鳥類が周囲に生息する場合に、その生存を脅かすことになる。

【0003】

そこで、風力発電事業のアセスメントとして、風車に衝突する鳥類衝突数を定量評価することが求められる。環境省は、ホームページ上に、平成２３年１月７日付けで「鳥類等に関する風力発電施設立地適正化のための手引き」を公開している（非特許文献１）。

【0004】

同手引きには、鳥類衝突数定量評価モデルとして、鳥類の飛翔高度と風車のブレードの回転域の高さとの関係から衝突数を定量評価する手法が記載されている。

【先行技術文献】

【非特許文献】

【0005】

【非特許文献1】ｈｔｔｐ：／／ｗｗｗ．ｅｎｖ．ｇｏ．ｊｐ／ｐｒｅｓｓ．ｐｈｐ？ｓｅｒｉａｌ＝１３３３１

【発明の概要】

【発明が解決しようとする課題】

【0006】

従来の鳥類衝突数定量評価モデルでは、風車ブレードの回転面の高度を飛翔し、回転面に進入する鳥類数を計数している。

【0007】

しかし、実際には、回転面に対して斜めに進入する鳥類もあり、これらは無視できない。

【0008】

本発明は、このような斜め進入も考慮して鳥類衝突数を定量評価する鳥類衝突数評価装置、方法及びプログラムを提示することを目的とする。

【課題を解決するための手段】

【0009】

本発明に係る鳥類衝突数評価装置は、風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価装置であって、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出手段と、注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力手段と、注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出手段と、当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出手段と、当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出手段とを具備し、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。

【0010】

本発明に係る鳥類衝突数評価装置は、下記式に従い、風力発電設備の風車に衝突する鳥類の羽数を定量評価する鳥類衝突数評価装置であって、
衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（１−回避率）
である。ただし、
Ｃｎ＝（Ｖｓ／Ｖａ）×Ｍｄ／（ｄ_ａｖｇ×２）
Ｖａ＝Ｓａ×Ｈｃ
Ｖｓ＝Ｍ×（４π／３）Ｒａ^３
Ｓａ：風力発電設備の設置域の面積
Ｈｃ：風車回転高度幅
Ｍ：風車台数
Ｒａ：風車のブレード半径又は風車軸からブレード端までの距離
Ｍｄ：注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
ｄ_ａｖｇ：風車接触可能域を通過する平均距離
接触率：ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率
回避率：当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該ブレード回転面への進入を回避する個体数の割合

【0011】

本発明に係る鳥類衝突数評価方法は、風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法であって、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出ステップと、注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力ステップと、注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出ステップと、当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出ステップと、当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出ステップとを具備し、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。

【0012】

本発明に係る鳥類衝突数評価プログラムは、風力発電設備の風車への鳥類衝突数を定量評価するコンピュータプログラムであって、コンピュータに、当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出機能と、注目する１以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力機能と、注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか１台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出機能と、当該注目する鳥類について、当該ブレード回転面に進入した１個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出機能と、当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出機能とを実現させ、当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定されることを特徴とする。

【発明の効果】

【0013】

本発明によれば、個々の風車に対して球体モデルによる風車接触可能域を設定したことで、全球方向から風車に向かう個体の衝突確率を適切に定量評価することが可能になる。これにより、衝突羽数をより精確に見込むことが可能になる。

【図面の簡単な説明】

【0014】

【図1】本発明の一実施例の概略構成ブロック図である。

【図2】本実施例の動作フローチャートである。

【図3】風力発電設備情報テーブルの構成例である。

【図4】鳥類情報テーブルの構成例である。

【図5】風車のパラメータを示す正面図である。

【図6】風車接触可能域への進入確率を説明する、風車を上方向から見た平面図である。

【図7】接触確率を説明する模式図である。

【図8】進入角度対接触確率のグラフ例である。

【発明を実施するための形態】

【0015】

以下、図面を参照して、本発明の実施例を詳細に説明する。

【実施例1】

【0016】

図１は、本発明の一実施例の概略構成ブロック図を示し、図２は、本実施例による鳥類衝突数評価の動作フローチャートを示す。本実施例では、コンピュータ上のソフトウエアにより実現されうるが、もちろん、そのようなソフトウエアで実現される機能を組み込んだ専用装置によっても実現可能である。

【0017】

本実施例では、風力発電設備の風車に対して全球方向から鳥類が進入する可能性が一般的に存在することを考慮し、個々の風車に対して衝突可能性のある鳥類の進入エリアとして風車ブレードの回転に外接する球体を風車接触可能域として想定する。

【0018】

コンピュータ１０は、ＣＰＵ１２、主記憶としてのメモリ１４、補助記憶としてのＨＤＤ１６、表示装置１８、キーボード２０及び、ポインティングデバイスとしてのマウス２２を具備し、鳥類衝突数評価装置の一実施例として機能する。図２に示すフローチャートに相当する鳥類衝突数評価のコンピュータプログラムがＨＤＤ１６に格納されており、ＣＰＵ１２は、その鳥類衝突数評価プログラムをＨＤＤ１６からメモリ１４に読み出し、実行する。これにより、ＣＰＵ１２は、図２に示すフローチャートで説明される各機能を実現する。

【0019】

コンピュータ１０には、予め、設置しようとする風力発電設備情報と、風力発電設備を設置しようとする地域に生息し、衝突可能性のある鳥類の情報が格納される。ＨＤＤ１６の風力発電設備情報テーブル１６ａには風力発電設備情報が格納され、鳥類情報テーブル１６ｂには鳥類情報が格納される。

【0020】

図３は、風力発電設備情報テーブル１６ａの構造例を示す。設置面積（Ｓａ）は、風力発電設備情報を設置しようとする地域の面積であって、例えば、設置する風車に衝突する可能性のある鳥類の生息域をカバーする範囲で設定される。明確さの点で、設置面積（Ｓａ）を、設置する複数の風力発電設備を外接する多角形その他の線図で定義しても良い。稼働率は１年で風車が回転する時間の割合を示す。風車軸高度（Ｈ_ａｘｉｓ）は、風車の回転軸の、地上からの高度を示す。ブレード半径（Ｒａ）は、風車を構成する羽（ブレード）の半径、即ち、風車軸からブレードの端までの距離を示す。平均回転速度は、想定される風速に対する平均的な回転速度を示す。平均回転速度は、風車を通過する鳥類がブレードに接触する接触確率を計算するのに使用される。図５は、風車を正面から見た状態の模式図であって、風車軸高度（Ｈ_ａｘｉｓ）、ブレード半径（Ｒａ）及び風車回転高度幅（２Ｒａ）の関係を示す。

【0021】

図４は、鳥類情報テーブル１６ｂのデータ構造例を示す。鳥類情報テーブル１６ｂは、衝突可能性のある鳥類毎に、サイズ（体長Ｉ，翼開長Ｗ）、飛翔高度範囲（Ｈａ〜Ｈｂ）、飛翔速度（ｖ）、全飛翔距離（Ｍｄ）及び回避率（Ｐ_Ｌ）の各データからなり、専ら実測により決定される。サイズは、羽を拡げた飛翔時の飛翔方向の大きさ（体長：Ｉ）と羽方向の大きさ（翼開長Ｗ）を示す。飛翔高度範囲（Ｈａ〜Ｈｂ）は、通常、飛翔する高度の範囲を示す。飛翔速度（ｖ）は、風車のブレードが回転する範囲の高度での飛翔速度を示す。全飛翔距離（Ｍｄ）は、単位期間（例えば、１年）内に風車回転高度範囲（Ｈａ〜Ｈｂ）を飛翔する場合の設置地域内での全飛翔個体の合計距離を示す。回避率（Ｐ_Ｌ）は、風車設置前に風車回転高度範囲（Ｈａ〜Ｈｂ）に侵入することが観察された個体数のうち、回転する風車に向かって飛びつつも風車を回避する割合を示す。風車を回避する個体数は、より具体的には、後述するように、設置された風車を忌避して設置地域に入らなくなる個体数と、設置地域に入っても風車接触可能域に入らない個体数と、風車接触可能域に入ってもブレード回転面に進入しない個体数の合計に相当する。例えば、オジロワシの場合、回避率は９５％程度と言われている。

【0022】

図２を参照して、本実施例による鳥類衝突数の定量評価動作を説明する。先ず、オペレータは、キーボード２０及びマウス２２を使い、設置しようとする風力発電設備情報とその地域の鳥類情報をＨＤＤ１６に格納する（Ｓ１）。これにより、ＨＤＤ１６に風力発電設備情報テーブル１６ａと鳥類情報テーブル１６ｂが格納され、準備される。勿論、遠隔地のサーバにある同様のテーブルを参照してもよいが、ここでは理解を容易にするために、ローカルにこれらテーブル１６ａ，１６ｂが存在するものとしている。

【0023】

ＣＰＵ１２は、風力発電設備情報テーブル１６ａを参照し、風車回転高度容積（Ｖａ）を算出する（Ｓ２）。風車回転高度容積（Ｖａ）は、テーブル１６ａの設置面積（Ｓａ）に風車回転高度幅（Ｈｃ）を乗算して得られる。すなわち、
Ｖａ＝Ｓａ×Ｈｃ
とする。風車回転高度幅（Ｈｃ）は、図５に示すように、回転する風車の上端から下端までの距離であり、一般的には、ブレード半径（Ｒａ）の２倍に相当する。風車回転高度容積（Ｖａ）は、高度で見た空間的な存在可能性を示す指数となる。

【0024】

風車回転高度幅の範囲に鳥類が飛翔していたとしても、ある程度に風車に接近しない限り、風車に衝突することはあり得ない。また、鳥類が風車に衝突するのは、鳥類が風車の正面から風車に向かって進入する場合に限らず、斜めに進入する場合もありうる。却って、風車回転面に直角に進入する場合に比べ、斜めに入射する方が、風車のブレードに接触する時間的空間的確率が高いので、風車に衝突する確率が高くなる。そこで、本実施例では、個々の風車に対して、風車のブレード回転面を中心とする球体を風車接触可能域として想定する球体モデルを採用する。

【0025】

ＣＰＵ１２は、風車接触可能域として想定される全球体の容積を風車接触可能容積（Ｖｓ）として、
Ｖｓ＝Ｍ×（４π／３）Ｒａ^３
を計算する（Ｓ３）。但し、Ｍは風車台数を示し、Ｒａは、風車のブレード半径（又は風車軸からブレード端までの距離）を示す。

【0026】

風車回転高度容積（Ｖａ）に対する風車接触可能容積（Ｖｓ）の比は、空間的な広さだけから見た場合に風車に衝突する可能性のある風車接触可能域に鳥類が存在する確率を示す。この確率を空間占有率と定義し、ＣＰＵ１２は、下記式に従い、設置地域における空間占有率（Ｐｖ）を算出する（Ｓ４）。すなわち、
Ｐｖ＝Ｖｓ／Ｖａ
とする。

【0027】

空間占有率（Ｐｖ）は、鳥類の種別に依存しない指数であり、本実施例は、この空間占有率（Ｐｖ）に鳥類毎の事情を加味して、個々の鳥類の衝突羽数を推計する。

【0028】

先ず、ＣＰＵ１２は、鳥類情報テーブル１６ｂを参照し、鳥類情報テーブル１６ｂから評価対象の鳥類を選択し、選択した鳥類の緒データ（サイズ（Ｉ，Ｗ）、飛翔高度範囲（Ｈａ〜Ｈｂ）、全飛翔速度（ｖ）、飛翔距離（Ｍｄ）及び回避率（Ｐ_Ｌ））を読み込む（Ｓ５）。

【0029】

ＣＰＵ１２は、選択された鳥類が風車接触可能域を通過する通過頻度（Ｆｑ）を算出する（Ｓ６）。通過頻度（Ｆｑ）は、選択された鳥類の全個体が所定の期間内に風車接触可能域（風車接触可能容積（Ｖｓ））を通過する回数である。選択された鳥類について、風車回転高度幅Ｈａ〜Ｈｂを飛翔する全個体の全飛翔距離（Ｍｄ）に空間占有率（Ｐｖ）を乗算した結果は、風車接触可能域を通過するのべ距離であるから、１台の風車の風車接触可能域を通過する平均距離（ｄ_ａｖｇ）に風車接触可能域の通過頻度（Ｆｑ）を乗算した結果に等しい。すなわち、
Ｆｑ×ｄ_ａｖｇ＝Ｐｖ×Ｍｄ
であり、これを変形して、
Ｆｑ＝Ｐｖ×Ｍｄ／ｄ_ａｖｇ
となる。選択された鳥類について、ＣＰＵ１２は、上式により通過頻度（Ｆｑ）を算出する。通過頻度（Ｆｑ）の単位は、鳥類種類別の所定期間内の回数となる。

【0030】

球体をランダムに通過するとした場合の平均距離（ｄ_ａｖｇ）は、球体にランダムに１本の線を通過させたときの、始点（入口）から終点（出口）までの距離の平均であり、通過軌跡の総和を軌跡数で除算することで得られる。従って、平均距離（ｄ_ａｖｇ）は、球の体積、すなわち、風車接触可能域を構成する球体の体積を、その中心を通る円の面積で除算した値となる。

【0031】

次に、ＣＰＵ１２は、何れかの風車の風車接触可能域に進入した鳥類が更に、そのブレード回転面を通過する回数（Ｃｎ）を算出する（Ｓ７）。ステップＳ６で計算した通過頻度（Ｆｑ）は、風車のブレード回転面を中心とする球体に進入する頻度を示すが、球体に進入した鳥がそのまま直進したとしても、全てがブレード回転面に進入するわけではなく、以下に説明するように、ブレード回転面に進入する確率は１／２となる。従って、ブレード回転面通過回数（Ｃｎ）は、
Ｃｎ＝Ｆｑ／２＝Ｐｖ×Ｍｄ／（ｄ_ａｖｇ×２）
となる。

【0032】

図６を参照して、ブレード回転面に進入する確率が１／２であることを説明する。図６は、風車をその上空から見た平面図を示す。本実施例の球体モデルでは、ブレード回転面５０に対して、外接する球体５２を想定する。この球体５２の内部が風車接触可能域に対応する。この球体５２に外部から任意の角度で鳥が線５４で示すように進入するとする。入射点とブレード回転面の両端を結ぶ線のなす角は周知の通り９０度となり、球体５２の球面への進入点での接線で見ると、余角、すなわちブレード面方向に進まない角度も９０度である。これは、進入する個体がブレード面を見込む立体角と、ブレード面以外を見込む立体角の比で見ても、同様である。従って、先に説明したように、風車接触可能域を構成する球体５２に進入した鳥類がブレード面に突入する確率は常に１／２となる。

【0033】

ここまでは、鳥類が直進することを前提としてきた。しかし、実際には、ブレード回転面に向かって飛翔する鳥類が全てブレード回転面に進入するとは限らない。回避行動があり、回避しない個体も、ブレード回転面を無事に通過できる場合がある。回避できず、ブレード回転面を無事に通過できなかった個体が、ブレードに衝突することになる。

【0034】

ブレード回転面に接近した鳥類も、その飛翔能力によりブレード回転面への進入を回避できる。このような回避率は、経験的又は実測により定量評価可能であるが、風車設置前の調査で設置域（設置前であるので、いわば「設置予定域」である。）に風車回転高度範囲に侵入することが観察された個体数（Ｔｘ）のうち、風車ブレード面への進入を回避する個体数（回避個体数）の割合である。回避個体数は、設置された風車を忌避して設置域に入らなくなる個体数と、設置域に入っても風車接触可能域に入らない個体数と、風車接触可能域に入ってもブレード回転面に進入しない個体数の合計として算出または計測できる。すなわち、
回避率＝（回避個体数）／Ｔｘ
である。なお、（Ｔｘ?回避個体数）がブレード面に進入する個体数に相当するので、（Ｔｘ?回避個体数）／Ｔｘが回転面進入率に相当する。本実施例でも、既存の方法で計測又は決定された知られた数値を採用する。例えば、オジロワシの場合で、回避率は、９５％程度と言われている。

【0035】

回避せずにブレード回転面に進入しても、一部は、ブレード回転面を無事に通過できる。無事に通過できるかどうかは、鳥類のサイズ（Ｉ，Ｗ）と、飛翔速度（ｖ）、ブレード面への入射角度（θ）、風車の回転速度及びブレード数に依存する。図７は、ブレード回転面を上方から見た場合の２種類の進入角度例であって、図７（ａ）はブレード回転面への進入角度が浅い（θが小さい）場合を示し、同（ｂ）は進入角度が大きい（θが大きい）場合を示す。

【0036】

鳥類の頭がブレード回転面に入ってから尾が抜けるまでの間にブレードが鳥類に接触する場合に、鳥類がブレードに衝突したことになる。進入角度θが小さいほどブレード回転面に位置する時間が長くなるので、接触確率が増大する。進入角度θが大きくなると、短時間でブレード回転面を通過することが可能になり、接触確率が低下する。進入角度θを０度から９０度の範囲で接触確率を算定すると、平均接触率Ｔを算定できる。すなわち、進入角度θの接触確率は、｛（ブレード回転面通過距離）／（飛翔速度ｖ）｝×｛（ブレード羽数）／（１回転の時間）｝で計算できる。平均衝突率Ｔは、各進入角度の衝突確率を進入角度で積分して平均化することで得られる。平均の代わりに、代表的な進入角度における接触確率を代表値として採用しても良い。

【0037】

ブレード回転面通過距離は、以下のように算出できる。進入角度θで野鳥がブレード面に進入する場合で、体長Ｉの１／２の箇所で翼開長Ｗを十字クロスさせると、
通過距離＝Ｗ／ｔａｎθ
となる。但し、オジロワシ、オオワシ、イヌワシ及びチュウヒの４種で推定した通過距離は、進入角度θが約６８度以上で各種類の体長（Ｉ）と同じ値になる。ブレード厚みｂを考慮する場合、
通過距離＝Ｗ／ｔａｎθ＋ｂ／ｓｉｎθ
となる。

【0038】

オジロワシを例に説明する。オジロワシのオスメス平均の体長及び翼開長から求めたブレード面を通過するのに要する距離、すなわちブレード面通過距離は、平均で２．３ｍ程であり、飛翔速度は１０．６ｍ／ｓである。このとき、ブレード回転面を垂直に通過する時間ｔは、
ｔ＝（ブレード面通過距離）／飛翔速度
であるから、０．２１７秒（＝２．３／１０．６）となる。

【0039】

風車のブレードが２．５秒に１回転すると仮定すると、ブレード１枚がｔ秒間に移動する間にブレード回転面を掃引する面積率は０．２１７／２．５となる。風車のブレードが３枚である場合、合計の掃引面積率は、
３×０．２１７／２．５＝０．２６
となる。つまり、合計の掃引面積率は、
ｔ×（ブレード数）／（１回転の秒数）
で得られる。この合計掃引面積率が９０度進入の場合の接触確率となる。進入角度θを変更して同様の計算を行い、０〜９０度の角度別の接触確率を平均化することで、平均接触率が得られる。図８は、進入角度対接触確率の一例を示す。

【0040】

以上を踏まえて、ＣＰＵ１２は、所定期間内に風車に衝突する羽数、すなわち期間内衝突羽数を、下記式に従って算定する（Ｓ８）。すなわち、
期間内衝突羽数＝Ｃｎ×（接触率）×（稼働率）×（回転面進入率）
ただし、
回転面進入率＝１−回避率
である。

【0041】

別種の野鳥についても衝突羽数を算定すべき場合（Ｓ９）、ＣＰＵ１２は、別の鳥類を選択して、緒データを鳥類情報テーブル１６ｂから読み込み（Ｓ１０）、当該鳥類について、ステップＳ６〜Ｓ８を実行する。

【0042】

全ての野鳥について衝突羽数を算出したら（Ｓ９）、ＣＰＵ１２は、所定書式のレポートを生成し、ＨＤＤ１６に保存して（Ｓ１１）、終了する。

【0043】

オペレータは、算出結果に疑問がある場合には、テーブル１６ａ，１６ｂのデータを修正して図２に示すフローを再実行すればよい。

【0044】

上記実施例では、ステップＳ６で通過頻度の計算にＳ４で得られた空間占有率を乗算しているが、ステップＳ６、Ｓ７では割合のみを計算し、最終的にステップＳ８において空間占有率を乗算するようにしてもよいことは、明らかである。

【0045】

このように、本実施例によれば、個々の風車に対して球体モデルによる風車接触可能域を設定したことで、全球方向から風車に向かう個体の衝突確率を適切に定量評価することが可能になり、従って、衝突羽数をより精確に見込むことが可能になる。

【0046】

特定の説明用の実施例を参照して本発明を説明したが、特許請求の範囲に規定される本発明の技術的範囲を逸脱しないで、上述の実施例に種々の変更・修整を施しうることは、本発明の属する分野の技術者にとって自明であり、このような変更・修整も本発明の技術的範囲に含まれる。

【符号の説明】

【0047】

１０：コンピュータ
１２：ＣＰＵ
１４：メモリ
１６：ＨＤＤ１
１６ａ：風力発電設備情報テーブル
１６ｂ：鳥類情報テーブル
１８：表示装置
２０：キーボード
２２：マウス
５０：ブレード回転面
５２：風車接触可能域（球体）
５４：進入方向

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】