JP6016211B2 - 鳥類衝突数評価装置、方法及びプログラム - Google Patents
鳥類衝突数評価装置、方法及びプログラム Download PDFInfo
- Publication number
- JP6016211B2 JP6016211B2 JP2012180105A JP2012180105A JP6016211B2 JP 6016211 B2 JP6016211 B2 JP 6016211B2 JP 2012180105 A JP2012180105 A JP 2012180105A JP 2012180105 A JP2012180105 A JP 2012180105A JP 6016211 B2 JP6016211 B2 JP 6016211B2
- Authority
- JP
- Japan
- Prior art keywords
- windmill
- bird
- rate
- collision
- contact
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Wind Motors (AREA)
Description
衝突羽数=Cn×(接触率)×(稼働率)×(1−回避率)
である。ただし、
Cn=(Vs/Va)×Md/(davg×2)
Va=Sa×Hc
Vs=M×(4π/3)Ra3
Sa:風力発電設備の設置域の面積
Hc:風車回転高度幅
M:風車台数
Ra:風車のブレード半径又は風車軸からブレード端までの距離
Md:注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
davg:風車接触可能域を通過する平均距離
接触率:ブレード回転面に進入した1個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率
回避率:当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該ブレード回転面への進入を回避する個体数の割合
Va=Sa×Hc
とする。風車回転高度幅(Hc)は、図5に示すように、回転する風車の上端から下端までの距離であり、一般的には、ブレード半径(Ra)の2倍に相当する。風車回転高度容積(Va)は、高度で見た空間的な存在可能性を示す指数となる。
Vs=M×(4π/3)Ra3
を計算する(S3)。但し、Mは風車台数を示し、Raは、風車のブレード半径(又は風車軸からブレード端までの距離)を示す。
Pv=Vs/Va
とする。
Fq×d avg =Pv×Md
であり、これを変形して、
Fq=Pv×Md/d avg
となる。選択された鳥類について、CPU12は、上式により通過頻度(Fq)を算出する。通過頻度(Fq)の単位は、鳥類種類別の所定期間内の回数となる。
Cn=Fq/2=Pv×Md/(davg×2)
となる。
回避率=(回避個体数)/Tx
である。なお、(Tx?回避個体数)がブレード面に進入する個体数に相当するので、(Tx?回避個体数)/Txが回転面進入率に相当する。本実施例でも、既存の方法で計測又は決定された知られた数値を採用する。例えば、オジロワシの場合で、回避率は、95%程度と言われている。
通過距離=W/tanθ
となる。但し、オジロワシ、オオワシ、イヌワシ及びチュウヒの4種で推定した通過距離は、進入角度θが約68度以上で各種類の体長(I)と同じ値になる。ブレード厚みbを考慮する場合、
通過距離=W/tanθ+b/sinθ
となる。
t=(ブレード面通過距離)/飛翔速度
であるから、0.217秒(=2.3/10.6)となる。
3×0.217/2.5=0.26
となる。つまり、合計の掃引面積率は、
t×(ブレード数)/(1回転の秒数)
で得られる。この合計掃引面積率が90度進入の場合の接触確率となる。進入角度θを変更して同様の計算を行い、0〜90度の角度別の接触確率を平均化することで、平均接触率が得られる。図8は、進入角度対接触確率の一例を示す。
期間内衝突羽数=Cn×(接触率)×(稼働率)×(回転面進入率)
ただし、
回転面進入率=1−回避率
である。
12:CPU
14:メモリ
16:HDD1
16a:風力発電設備情報テーブル
16b:鳥類情報テーブル
18:表示装置
20:キーボード
22:マウス
50:ブレード回転面
52:風車接触可能域(球体)
54:進入方向
Claims (21)
- 風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価装置であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出手段と、
注目する1以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力手段と、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか1台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出手段と、
当該注目する鳥類について、当該ブレード回転面に進入した1個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出手段と、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出手段
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価装置。 - 当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項1に記載の鳥類衝突数評価装置。
- さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する手段と、
当該回避率を1から減算した値を当該回転面進入率とする演算手段
とを具備することを特徴とする請求項1または2に記載の鳥類衝突数評価装置。 - 当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項1から3の何れか1項に記載の鳥類衝突数評価装置。
- 当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項1から3の何れか1項に記載の鳥類衝突数評価装置。
- 更に、当該鳥類とは別の鳥類について、当該通過頻度算出手段、当該接触率算出手段及び当該衝突羽数算出手段を使って当該衝突羽数を算出させる制御手段を具備することを特徴とする請求項1乃至5の何れか1項に記載の鳥類衝突数評価装置。
- 下記式に従い、風力発電設備の風車に衝突する鳥類の羽数を定量評価する鳥類衝突数評価装置。
衝突羽数=Cn×(接触率)×(稼働率)×(1−回避率)
ただし、
Cn=(Vs/Va)×Md/(davg×2)
Va=Sa×Hc
Vs=M×(4π/3)Ra3
Sa:風力発電設備の設置域の面積
Hc:風車回転高度幅
M:風車台数
Ra:風車のブレード半径又は風車軸からブレード端までの距離
Md:注目する鳥類について風車回転高度幅での所定期間の全飛翔距離
davg:風車接触可能域を通過する平均距離
接触率:ブレード回転面に進入した1個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率
回避率:当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該ブレード回転面への進入を回避する個体数の割合 - 当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項7に記載の鳥類衝突数評価装置。
- 当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項7に記載の鳥類衝突数評価装置。
- 風力発電設備の風車への鳥類衝突数を定量評価する鳥類衝突数評価方法であって、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出ステップと、
注目する1以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力ステップと、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか1台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出ステップと、
当該注目する鳥類について、当該ブレード回転面に進入した1個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出ステップと、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出ステップ
とを具備し、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価方法。 - 当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項10に記載の鳥類衝突数評価方法。
- さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する入力ステップと、
当該回避率を1から減算した値を当該回転面進入率とする演算ステップ
とを有することを特徴とする請求項10または11に記載の鳥類衝突数評価方法。 - 当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項10から12の何れか1項に記載の鳥類衝突数評価方法。
- 当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項10から12の何れか1項に記載の鳥類衝突数評価方法。
- 更に、当該鳥類とは別の鳥類について、当該通過頻度算出ステップ、当該接触率算出ステップ及び当該衝突羽数算出ステップに当該衝突羽数を算出させる制御ステップを具備することを特徴とする請求項10乃至14の何れか1項に記載の鳥類衝突数評価方法。
- 風力発電設備の風車への鳥類衝突数を定量評価するコンピュータプログラムであって、コンピュータに、
当該風力発電設備の設置域における当該風車の風車接触可能域の空間占有率を算出する空間占有率算出機能と、
注目する1以上の鳥類のそれぞれについて、サイズ、当該風車の風車回転高度範囲での飛翔速度、及び単位期間内に当該風車回転高度範囲を飛翔する場合の当該設置域内での全飛翔個体の合計距離を示す全飛翔距離を入力する入力機能と、
注目する鳥類について、当該全飛翔距離及び当該空間占有率から、当該設置域を飛翔する個体が当該設置域におけるいずれか1台の風車のブレード回転面を通過する通過頻度を算出する通過頻度算出機能と、
当該注目する鳥類について、当該ブレード回転面に進入した1個体が当該ブレード回転面を通過しきれずに当該ブレード回転面を構成するブレードに接触する確率を示す接触率を、当該鳥類のサイズ及び飛翔速度、当該ブレード回転面への進入角度、並びに当該風車の回転速度及びブレード数に基づき算出する接触率算出機能と、
当該通過頻度に当該接触率、当該ブレード回転面を回避せずに当該ブレード回転面に進入してしまう経験的または実測される確率を示す回転面進入率、及び当該風車の稼働率を乗算して衝突羽数を算出する衝突羽数算出機能
とを実現させ、
当該風車接触可能域が当該風車のブレード回転面を中心とする球体で規定される
ことを特徴とする鳥類衝突数評価プログラム。 - 当該空間占有率が、当該風車接触可能域の総容積を、当該設置域の当該風車の風車回転高度容積で除算した結果であることを特徴とする請求項16に記載の鳥類衝突数評価プログラム。
- さらに、
当該設置域で当該風車回転高度範囲に侵入する個体数のうち、当該設置域に入らなくなる個体数と、当該設置域に入っても当該風車接触可能域に入らない個体数と、当該風車接触可能域に入っても当該ブレード回転面に進入しない個体数の合計の割合を経験的または実測により示す回避率を入力する入力機能と、
当該回避率を1から減算した値を当該回転面進入率とする演算機能
とを当該コンピュータに実現させることを特徴とする請求項16又は17に記載の鳥類衝突数評価プログラム。 - 当該接触率が、当該ブレード回転面に所定の進入角度範囲で計算される各進入角度の接触確率の平均値であることを特徴とする請求項16から18の何れか1項に記載の鳥類衝突数評価プログラム。
- 当該接触率が、当該ブレード回転面の代表的な進入角度で計算される接触確率であることを特徴とする請求項16から18の何れか1項に記載の鳥類衝突数評価プログラム。
- 当該鳥類とは別の鳥類について、当該通過頻度算出機能、当該接触率算出機能及び当該衝突羽数算出機能に当該衝突羽数を算出させる制御機能を当該コンピュータに実現させることを特徴とする請求項16乃至20の何れか1項に記載の鳥類衝突数評価プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012180105A JP6016211B2 (ja) | 2012-08-15 | 2012-08-15 | 鳥類衝突数評価装置、方法及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012180105A JP6016211B2 (ja) | 2012-08-15 | 2012-08-15 | 鳥類衝突数評価装置、方法及びプログラム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2014037795A JP2014037795A (ja) | 2014-02-27 |
JP2014037795A5 JP2014037795A5 (ja) | 2015-05-28 |
JP6016211B2 true JP6016211B2 (ja) | 2016-10-26 |
Family
ID=50286106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012180105A Active JP6016211B2 (ja) | 2012-08-15 | 2012-08-15 | 鳥類衝突数評価装置、方法及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6016211B2 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6316638B2 (ja) * | 2014-04-04 | 2018-04-25 | アジア航測株式会社 | 監視装置、監視方法および監視プログラム |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005046860A1 (de) * | 2005-09-29 | 2007-04-05 | Daubner & Stommel GbR Bau-Werk-Planung (vertretungsberechtigter Gesellschafter: Matthias Stommel, 27777 Ganderkesee) | Verfahren zur Regelung einer Windenergieanlage |
JP2010071100A (ja) * | 2008-09-16 | 2010-04-02 | Tokyo Electric Power Co Inc:The | 高速衝突判定処理システムおよびコンピュータプログラム |
FR2939902A1 (fr) * | 2008-12-16 | 2010-06-18 | Henri Pierre Roche | Systeme de detection d'oiseaux et d'arret automatise d'eolienne industrielle |
-
2012
- 2012-08-15 JP JP2012180105A patent/JP6016211B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2014037795A (ja) | 2014-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eichhorn et al. | Model-based estimation of collision risks of predatory birds with wind turbines | |
EP2667022A2 (en) | Method and system for wind velocity field measurements on a wind farm | |
JP4626265B2 (ja) | 風力発電装置、風力発電装置の制御方法およびコンピュータプログラム | |
Bevelhimer et al. | Hydroacoustic assessment of behavioral responses by fish passing near an operating tidal turbine in the east river, New York | |
WO2012097076A2 (en) | Methods and apparatus for monitoring complex flow fields for wind turbine applications | |
Subramanian et al. | Aerodynamics of wind turbine wakes in flat and complex terrains | |
Carneiro et al. | Investigation of possible societal risk associated with wind power generation systems | |
MX2010008512A (es) | Sistemas y metodos para mitigar los efectos de las turbinas eolicas en un radar. | |
US20180171972A1 (en) | System and method for preventing collisions between wind turbine blades and flying objects | |
Machefaux et al. | Investigation of wake interaction using full‐scale lidar measurements and large eddy simulation | |
JP7061874B2 (ja) | 風力事象の到着を予測するためのシステムおよび方法 | |
Smales et al. | A description of the Biosis model to assess risk of bird collisions with wind turbines | |
JP2010071100A (ja) | 高速衝突判定処理システムおよびコンピュータプログラム | |
Li et al. | Investigation of the dynamic radar signatures of a vertical-axis wind turbine | |
JP6016211B2 (ja) | 鳥類衝突数評価装置、方法及びプログラム | |
Kaikwong et al. | Aircraft wake fortex observations in Hong Kong | |
Holmstrom et al. | Assessing avian-wind turbine collision risk: an approach angle dependent model | |
Santos et al. | Factors influencing wind turbine avoidance behaviour of a migrating soaring bird | |
Nguyen et al. | A Monte Carlo simulation study of wind turbine loads in thunderstorm downbursts | |
EP3009670A1 (en) | Method to control the operation of a wind energy installation, including generating a 3d wind speed model, and wind energy installation | |
Rogers et al. | A method for defining wind turbine setback standards | |
Gallacher et al. | Lidar measurements and visualisation of turbulence and wake decay length | |
Wang et al. | Impact of wind turbines on ATC radars and mitigation results | |
Prakash et al. | Development and testing of a three‐dimensional ballistics model for bat strikes on wind turbines | |
JP3236051U (ja) | バードストライク回避・検証システム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150409 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150409 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160203 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160301 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160418 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160920 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160920 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6016211 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |