JP7009237B2 - Wind power generation equipment and wind power generation system - Google Patents

Wind power generation equipment and wind power generation system Download PDF

Info

Publication number
JP7009237B2
JP7009237B2 JP2018014369A JP2018014369A JP7009237B2 JP 7009237 B2 JP7009237 B2 JP 7009237B2 JP 2018014369 A JP2018014369 A JP 2018014369A JP 2018014369 A JP2018014369 A JP 2018014369A JP 7009237 B2 JP7009237 B2 JP 7009237B2
Authority
JP
Japan
Prior art keywords
wind
power generation
wind power
shear
load
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
Application number
JP2018014369A
Other languages
Japanese (ja)
Other versions
JP2019132178A (en
Inventor
伸夫 苗村
満 佐伯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2018014369A priority Critical patent/JP7009237B2/en
Priority to PCT/JP2018/046459 priority patent/WO2019150805A1/en
Priority to TW108102524A priority patent/TWI729349B/en
Publication of JP2019132178A publication Critical patent/JP2019132178A/en
Application granted granted Critical
Publication of JP7009237B2 publication Critical patent/JP7009237B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors 
    • F03D7/02Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/04Automatic control; Regulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)

Description

本発明は、風力発電装置の周りの風況を推定する機能を有する風力発電装置及び風力発電システムに関する。 The present invention relates to a wind power generation device and a wind power generation system having a function of estimating the wind condition around the wind power generation device.

再生可能エネルギー活用への関心の高まりから、風力発電装置の世界的な市場拡大が予測されている。メガワット級の風力発電装置としては、ブレードを回転するハブに放射状に取りつけたロータと、主軸を介してロータを支持するナセルと、ナセルを下部からヨー回転を許容し支持するタワーを備えているものが頻繁に用いられる。
風力発電装置では、時々刻々と変化する風をエネルギー源として発電を行う。従って、実際に風力発電装置に流入する風の風速や乱れが設計条件よりも厳しい場合、風力発電装置の負荷が増大し、構成部品の損傷が加速する可能性がある。損傷の加速により、予期せぬ故障が発生した場合、故障部品の交換に要する時間に加えて、交換用部品の手配や工事用機材、作業員の手配に時間を要するため、計画された部品交換を実施する場合よりも、風力発電装置の稼動停止時間が増大し、発電量が減少する懸念がある。
この対策として、エネルギー源となる流入風を推定・計測することで、風力発電装置の構成部品の損傷を推定する方法や、負荷を軽減するように制御する方法がある。例えば、特許文献1では、簡便な装置構成で入手可能な風速と、発電量若しくはブレードのピッチ角のデータから、予め作成しておいた風速、発電量若しくはピッチ角と、風況パラメータを関連付けるテーブルと、風況パラメータと風力発電装置に加わる疲労荷重を関連付けるテーブルを用いて、風力発電装置の構成部品の寿命を推定し、メンテナンス情報を出力する装置が提案されている。
一方、特許文献2では、一般にドップラーライダと称されるマイクロ波またはレーダ波発射装置を風力発電装置のナセルまたはハブに取り付け、風力発電装置前方および後方の風速分布を計測することで、発電効率を最大化または風力発電装置の負荷を最小化するようにブレードのピッチ角を制御する方法が提案されている。
The growing interest in renewable energy utilization is expected to expand the global market for wind power plants. Megawatt-class wind power generators are equipped with a rotor that is radially attached to a hub that rotates blades, a nacelle that supports the rotor via a spindle, and a tower that allows and supports yaw rotation from below. Is frequently used.
Wind power generators generate electricity using the ever-changing wind as an energy source. Therefore, if the wind speed or turbulence of the wind that actually flows into the wind power generation device is stricter than the design conditions, the load on the wind power generation device may increase and the damage to the components may be accelerated. If an unexpected failure occurs due to accelerated damage, it will take time to arrange replacement parts, construction equipment, and workers in addition to the time required to replace the failed parts, so planned parts replacement There is a concern that the downtime of the wind power generation equipment will increase and the amount of power generation will decrease compared to the case of implementing.
As a countermeasure, there are a method of estimating damage to the components of the wind power generation device by estimating and measuring the inflow wind as an energy source, and a method of controlling so as to reduce the load. For example, in Patent Document 1, a table that associates the wind speed, the amount of power generation, or the pitch angle created in advance from the data of the wind speed, the amount of power generation, or the pitch angle of the blade, which can be obtained with a simple device configuration, with the wind condition parameter. A device that estimates the life of the components of the wind power generation device and outputs maintenance information using a table that associates the wind condition parameters with the fatigue load applied to the wind power generation device has been proposed.
On the other hand, in Patent Document 2, a microwave or radar wave launcher generally called a Doppler lidar is attached to a nacelle or hub of a wind power generation device, and the wind speed distribution in front of and behind the wind power generation device is measured to improve power generation efficiency. Methods have been proposed to control the pitch angle of the blades to maximize or minimize the load on the wind turbine.

特開2015-117682号公報JP-A-2015-117682 特表2015-519516号公報Japanese Patent Application Laid-Open No. 2015-591516

しかしながら、特許文献1に開示される構成では、発電量やピッチ角は概ね風速に対して一意の値になるよう制御されているため、風速、発電量、ピッチ角から風速分布であるウィンドシアを高精度に推定することは困難であり、構成部品の寿命推定精度への影響も懸念される。また、特許文献2に開示される構成では、ドップラーライダにより風速分布を正確に計測できるものの、このような計測機器は高価であるため、風力発電所内(ウィンドファーム内)の全ての風力発電装置に設置することはコスト面から現実的ではない。 However, in the configuration disclosed in Patent Document 1, since the power generation amount and the pitch angle are controlled so as to be substantially unique values with respect to the wind speed, the wind shear which is the wind speed distribution from the wind speed, the power generation amount, and the pitch angle is obtained. It is difficult to estimate with high accuracy, and there is concern about the influence on the life estimation accuracy of components. Further, in the configuration disclosed in Patent Document 2, although the wind speed distribution can be accurately measured by the Doppler lidar, such a measuring device is expensive, so that it can be used for all wind power generation devices in a wind power plant (in a wind farm). Installation is not realistic in terms of cost.

そこで、本発明は、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供する。 Therefore, the present invention provides a wind power generation device and a wind power generation system capable of estimating wind shear, which is a wind speed distribution, with high accuracy with a simple configuration.

上記課題を解決するため、本発明に係る風力発電装置は、少なくともロータ及びナセル並びにナセルを回動可能に支持するタワーを有する風力発電装置と、風力発電装置を制御する制御装置とを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、前記荷重及び前記大気密度の平均値並びにウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
また、本発明に係る風力発電システムは、少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、前記荷重及び前記大気密度の平均値並びにウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
また、本発明に係る他の風力発電システムは、少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、を有し、前記電子端末は、荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重及び前記大気密度計算部から入力される大気密度の平均値並びに前記前記記憶部に格納されるウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
In order to solve the above problems, the wind power generation device according to the present invention includes at least a wind power generation device having a rotor, a nacelle, and a tower that rotatably supports the nacelle, and a control device for controlling the wind power generation device. The control device has a load measuring unit that measures the load applied to the wind power generator, a storage unit that stores the wind shear function that defines the relationship between the load and the wind shear, and a gas state equation from temperature and pressure. A wind condition having an air density calculation unit that calculates the air density using the wind shear and obtains the average value of the air density, and a wind speed distribution calculation unit that calculates the wind shear based on the load, the average value of the air density, and the wind shear function. It is characterized by being provided with an estimation device.
Further, the wind power generation system according to the present invention includes at least one wind power generation device, a control device for controlling the wind power generation device, an electronic terminal having a display device, and a communication network connecting these to each other so as to be communicable. The control device includes a load measuring unit that measures the load applied to the wind power generator, a storage unit that stores a wind shear function that defines the relationship between the load and the wind shear, and a gas from temperature and pressure. An air density calculation unit that calculates the air density using the state equation and obtains the average value of the air density, and a wind speed distribution calculation unit that calculates the wind shear based on the load, the average value of the air density, and the wind shear function. It is characterized by being provided with a wind condition estimation device.
Further, another wind power generation system according to the present invention is a communication that connects at least one wind power generation device, a control device for controlling the wind power generation device, and an electronic terminal having a display device so as to be communicable with each other. Equipped with a network, the control device has a load measuring unit that measures the load applied to the wind power generator , and the atmosphere that calculates the air density from the temperature and pressure using the gas state equation and obtains the average value of the air density. The electronic terminal has a density calculation unit, a storage unit for storing a wind shear function that defines a relationship between a load and a wind power, and the load input from the load measurement unit via the communication network. A wind condition estimation device including an average value of the air density input from the air density calculation unit and a wind speed distribution calculation unit for calculating wind power based on the wind shear function stored in the storage unit. And.

本発明によれば、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供することが可能となる。
上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
According to the present invention, it is possible to provide a wind power generation device and a wind power generation system capable of estimating wind shear, which is a wind speed distribution, with high accuracy with a simple configuration.
Issues, configurations and effects other than those described above will be clarified by the following description of the embodiments.

本発明の一実施形態に係る風力発電システムの全体概略構成図である。It is an overall schematic block diagram of the wind power generation system which concerns on one Embodiment of this invention. 本発明の一実施例に係る実施例1の風力発電装置の構成を示す図である。It is a figure which shows the structure of the wind power generation apparatus of Example 1 which concerns on one Example of this invention. 実施例1におけるウィンドシアを説明するための図である。It is a figure for demonstrating the wind shear in Example 1. FIG. 実施例1の風力発電装置を構成する制御装置の機能ブロック図である。It is a functional block diagram of the control device which constitutes the wind power generation device of Example 1. FIG. 図4に示す制御装置を構成する風況推定装置の原理を説明するための図である。It is a figure for demonstrating the principle of the wind condition estimation apparatus which constitutes the control apparatus shown in FIG. 冪指数とハブ中心でのモーメントの関係を説明するための図である。It is a figure for demonstrating the relationship between the exponentiation and the moment at the center of a hub. 冪指数とハブ中心でのモーメントの関係における風速の影響を説明するための図である。It is a figure for demonstrating the influence of the wind speed on the relationship between the exponent and the moment at the center of the hub. 冪指数とハブ中心でのモーメントの関係における空気密度の影響を説明するための図である。It is a figure for demonstrating the influence of the air density on the relationship between the exponent and the moment at the center of the hub. 本発明の他の実施例に係る実施例2の風力発電装置を構成する制御装置の機能ブロック図である。It is a functional block diagram of the control device which constitutes the wind power generation device of Example 2 which concerns on another Example of this invention. 本発明の他の実施例に係る実施例3の風力発電装置を構成する制御装置の機能ブロック図である。It is a functional block diagram of the control device which constitutes the wind power generation device of Example 3 which concerns on another Example of this invention. 本発明の他の実施例に係る実施例4の風力発電装置の構成を示す図である。It is a figure which shows the structure of the wind power generation apparatus of Example 4 which concerns on another Example of this invention. 本発明の他の実施例に係る実施例5の風力発電装置の構成を示す図である。It is a figure which shows the structure of the wind power generation apparatus of Example 5 which concerns on another Example of this invention.

図1は、本発明の一実施形態に係る風力発電システムの全体概略構成図である。図1に示すように、風力発電システム1は、風力発電装置2、及び、運転管理センター3内に設置される電子端末4或いは図示しないサーバを備え、これらは相互に通信可能に通信ネットワーク5を介して接続されている。なお、通信ネットワーク5は有線か無線かを問わない。 FIG. 1 is an overall schematic configuration diagram of a wind power generation system according to an embodiment of the present invention. As shown in FIG. 1, the wind power generation system 1 includes a wind power generation device 2 and an electronic terminal 4 or a server (not shown) installed in an operation management center 3, which enable a communication network 5 to communicate with each other. Connected via. The communication network 5 may be wired or wireless.

また、風力発電装置2は、風を受けて回転するブレード24、ブレード24を支持するハブ23、ナセル22、及びナセル22を回動可能に支持するタワー21を備える。ナセル22内に、ハブ23に接続されハブ23と共に回転する主軸26、主軸26に接続され回転速度を増速する増速機27、及び増速機27により増速された回転速度で回転子を回転させて発電運転する発電機28を備えている。ブレード24の回転エネルギーを発電機28に伝達する部位は、動力伝達部と称され、本実施形態では、主軸26及び増速機27が動力伝達部に含まれる。そして、増速機27及び発電機28は、メインフレーム29上に保持されている。また、ブレード24及びハブ23によりロータ25が構成される。図1に示すように、タワー21内の底部(下部)に、電力の周波数を変換する電力変換器30、電流の開閉を行うスイッチング用の開閉器及び変圧器など(図示せず)、及び制御装置31などが配されている。制御装置31として、例えば、制御盤又はSCADA(Supervisory Control And Data Acquisition)が用いられる。 Further, the wind power generation device 2 includes a blade 24 that rotates in response to wind, a hub 23 that supports the blade 24, a nacelle 22, and a tower 21 that rotatably supports the nacelle 22. In the nacelle 22, a spindle 26 connected to the hub 23 and rotating together with the hub 23, a speed increasing machine 27 connected to the spindle 26 to increase the rotation speed, and a rotor at a rotation speed increased by the speed increasing machine 27. It is equipped with a generator 28 that is rotated to generate power. The portion that transmits the rotational energy of the blade 24 to the generator 28 is referred to as a power transmission unit, and in the present embodiment, the main shaft 26 and the speed increaser 27 are included in the power transmission unit. The speed increaser 27 and the generator 28 are held on the main frame 29. Further, the rotor 25 is composed of the blade 24 and the hub 23. As shown in FIG. 1, at the bottom (lower part) in the tower 21, a power converter 30 that converts the frequency of electric power, a switch and a transformer for switching that switches and switches current (not shown), and a control are provided. The device 31 and the like are arranged. As the control device 31, for example, a control panel or SCADA (Supervision Control And Data Acquisition) is used.

なお、図1に示す風力発電装置2は、3枚のブレード24とハブ23にてロータ25を構成する例を示すが、これに限られず、ロータ25はハブ23と少なくとも1枚のブレード24にて構成しても良い。
以下、図面を用いて本発明の実施例について説明する。
The wind power generation device 2 shown in FIG. 1 shows an example in which the rotor 25 is composed of three blades 24 and a hub 23, but the rotor 25 is not limited to this, and the rotor 25 includes a hub 23 and at least one blade 24. May be configured.
Hereinafter, examples of the present invention will be described with reference to the drawings.

図2は、本発明の一実施例に係る実施例1の風力発電装置の構成を示す図である。図2では、本実施例に係る風力発電装置2及びその周りでの高さ方向の風速分布11の構成を示している。また、図2では、風力発電装置2を側方から眺めた状態を示しており、風は紙面左から右に吹いているものとする。図2に示すように、風力発電装置2は、ブレード24を回転するハブ23に放射状に取りつけたロータ25と、ロータ25の回転を許容してハブ23を横方向から支持するナセル22と、ナセル22を下部から垂直軸に対して回転自在に支持するタワー21とを備える。風力発電装置2は、図2ではロータ25がタワー21の風下側に位置するダウンウインド型の風力発電装置を示しているが、ロータ25がタワー21の風上側に位置するアップウィンド型の風力発電装置でも良い。 FIG. 2 is a diagram showing a configuration of a wind power generation device according to a first embodiment of the present invention. FIG. 2 shows the configuration of the wind power generation device 2 according to the present embodiment and the wind speed distribution 11 in the height direction around the wind power generator 2. Further, FIG. 2 shows a state in which the wind power generation device 2 is viewed from the side, and it is assumed that the wind is blowing from the left to the right of the paper. As shown in FIG. 2, the wind power generation device 2 includes a rotor 25 radially attached to a hub 23 that rotates the blade 24, a nacelle 22 that allows the rotor 25 to rotate and supports the hub 23 from the lateral direction, and a nacelle. A tower 21 that rotatably supports the 22 from the lower part with respect to the vertical axis is provided. The wind power generation device 2 shows a downwind type wind power generation device in which the rotor 25 is located on the leeward side of the tower 21 in FIG. 2, but the rotor 25 is an upwind type wind power generation device located on the windward side of the tower 21. It may be a device.

風力発電装置2は、タワー21に取り付けられたひずみセンサ7を備えている。ひずみセンサ7については、タワー21に限らず、ナセル22またはハブ23に設置しても良く、ひずみセンサ7に代えて、例えば加速度センサのような他の荷重センサを用いても良い。また、風力発電装置2は、ブレード24のピッチ角を制御するためのピッチ角制御機構6、ナセル22の上部に設置された風速計8、ナセル22内に設置された温度計9、気圧計10を備えていても良い。なお、風速計8、温度計9、気圧計10は、風力発電装置2の他の位置に設置しても良く、風力発電装置2の近傍であれば風力発電装置2の外部に設置しても良い。また、温度計9、気圧計10については計測せず、気象予測などに用いられる公の観測データを用いても良い。 The wind power generator 2 includes a strain sensor 7 attached to the tower 21. The strain sensor 7 may be installed not only in the tower 21 but also in the nacelle 22 or the hub 23, and instead of the strain sensor 7, another load sensor such as an acceleration sensor may be used. Further, the wind power generator 2 includes a pitch angle control mechanism 6 for controlling the pitch angle of the blade 24, an anemometer 8 installed above the nacelle 22, a thermometer 9 installed in the nacelle 22, and a barometer 10. May be equipped. The anemometer 8, the thermometer 9, and the barometer 10 may be installed at other positions of the wind power generation device 2, or may be installed outside the wind power generation device 2 if they are in the vicinity of the wind power generation device 2. good. Further, the thermometer 9 and the barometer 10 may not be measured, and public observation data used for weather prediction or the like may be used.

風速分布11は、通常、高さ方向(Z方向)に変化しており、一般的には大気の境界層によって上空ほど風速が大きくなる傾向がある。この高さ方向(Z方向)の風速の変化をウィンドシアと称し、その強弱を表す冪指数をαWSとすると、高さ方向(Z方向)の風速分布を以下の式(1)のように仮定することができる。 The wind speed distribution 11 usually changes in the height direction (Z direction), and generally, the wind speed tends to increase toward the sky depending on the boundary layer of the atmosphere. This change in wind speed in the height direction (Z direction) is called wind shear, and if the index representing the strength is α WS , the wind speed distribution in the height direction (Z direction) is as shown in the following equation (1). You can assume.

Figure 0007009237000001
Figure 0007009237000001

ここで、V(z)は地表から高さzにおける風速であり、zrefは基準となる風速を定義する高さ、V(zref)は基準となる風速を表している。図3に示すように、冪指数αWSが大きくなるほど、高さによる風速の変化は大きくなる。例えば、基準となる風速を風速計8で計測する場合、式(1)のV(zref)に計測した風速、zrefに風速計8の地表からの高さを用いることで、ある冪指数αWSがわかれば高さ方向の風速分布を得ることができる。すなわち、風速分布を式(1)のように仮定する場合、風速分布の推定問題は、冪指数αWSの推定問題に帰着する。なお、本実施例では、式(1)の風速分布を仮定して、ウィンドシアの強弱を表す冪指数αWSを推定することを考えるが、風速分布の仮定は式(1)に限らず、対数や多項式などを用いても良く、複数のパラメータを用いて風速分布を仮定しても良い。 Here, V (z) is the wind speed at the height z from the ground surface, z ref is the height that defines the reference wind speed, and V (z ref ) is the reference wind speed. As shown in FIG. 3, the larger the exponent α WS , the larger the change in wind speed depending on the height. For example, when the reference wind speed is measured by the anemometer 8, the measured wind speed is used for V (z ref ) in the equation (1), and the height from the ground surface of the anemometer 8 is used for z ref . If the α WS is known, the wind speed distribution in the height direction can be obtained. That is, when the wind speed distribution is assumed as in the equation (1), the problem of estimating the wind speed distribution results in the problem of estimating the exponent α WS . In this embodiment, it is considered to estimate the wind shear index α WS representing the strength of the wind shear by assuming the wind speed distribution of the equation (1), but the assumption of the wind speed distribution is not limited to the equation (1). A logarithm, a polynomial, or the like may be used, or a plurality of parameters may be used to assume a wind speed distribution.

図4は、本実施例の風力発電装置2を構成する制御装置31の機能ブロック図である。図4に示すように、制御装置31は、風速計測部12、荷重計測部13、温度計測部14、気圧計測部15、ウィンドシア関数33を格納する記憶部16、大気密度計算部17、ウィンドシア推定部18、風速分布計算部19、入力I/F34、出力I/F35、及び、通信I/F36を備え、これらは相互に内部バス37にてアクセス可能に接続されている。風速計測部12、荷重計測部13、温度計測部14、気圧計測部15、ウィンドシア関数33を格納する記憶部16、大気密度計算部17、ウィンドシア推定部18、及び風速分布計算部19は、ウィンドシアの強弱を定義する冪指数αWSの推定のみならず、実際の風速分布を求めるため、風況推定装置32を構成する。なお、冪指数αWSの推定のみであれば、風況推定装置32を、風速計測部12、荷重計測部13、温度計測部14、気圧計測部15、ウィンドシア関数33を格納する記憶部16、大気密度計算部17、及びウィンドシア推定部18にて構成しても良い。更には、推定精度の低下を許容して冪指数αWSのみを推定する場合には、荷重計測部13、ウィンドシア関数33を格納する記憶部16、及びウィンドシア推定部18のみにて風況推定装置32を構成しても良い。
風況推定装置32を構成する、風速計測部12、荷重計測部13、温度計測部14、気圧計測部15、大気密度計算部17、ウィンドシア推定部18、及び風速分布計算部19は、例えば、図示しないCPU(Central Processing Unit)などのプロセッサ、各種プログラムを格納するROM、演算過程のデータを一時的に格納するRAM、外部記憶装置などの記憶装置にて実現されると共に、CPUなどのプロセッサがROMに格納された各種プログラムを読み出し実行し、実行結果である演算結果をRAM又は外部記憶装置に格納する。以下に、制御装置31の各部の詳細につき説明する。
FIG. 4 is a functional block diagram of the control device 31 constituting the wind power generation device 2 of this embodiment. As shown in FIG. 4, the control device 31 includes a wind speed measuring unit 12, a load measuring unit 13, a temperature measuring unit 14, a barometric pressure measuring unit 15, a storage unit 16 for storing a wind shear function 33, an atmospheric density calculation unit 17, and a window. A shear estimation unit 18, a wind speed distribution calculation unit 19, an input I / F34, an output I / F35, and a communication I / F36 are provided, and these are connected to each other so as to be accessible by an internal bus 37. The wind speed measurement unit 12, the load measurement unit 13, the temperature measurement unit 14, the barometric pressure measurement unit 15, the storage unit 16 for storing the wind shear function 33, the air density calculation unit 17, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 , The wind condition estimation device 32 is configured to obtain not only the estimation of the barometric pressure index α WS that defines the strength of the wind shear but also the actual wind speed distribution. If only the estimation of the exponentiation α WS is performed, the wind condition estimation device 32 is a storage unit 16 that stores the wind speed measuring unit 12, the load measuring unit 13, the temperature measuring unit 14, the barometric pressure measuring unit 15, and the wind shear function 33. , The atmospheric density calculation unit 17, and the wind shear estimation unit 18 may be configured. Furthermore, when only the exponent α WS is estimated by allowing the estimation accuracy to decrease, the wind condition is determined only by the load measuring unit 13, the storage unit 16 that stores the wind shear function 33, and the wind shear estimation unit 18. The estimation device 32 may be configured.
The wind speed measurement unit 12, the load measurement unit 13, the temperature measurement unit 14, the pressure measurement unit 15, the air density calculation unit 17, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 constituting the wind condition estimation device 32 are, for example,. , A processor such as a CPU (Central Processing Unit) (not shown), a ROM for storing various programs, a RAM for temporarily storing data of an arithmetic process, a storage device such as an external storage device, and a processor such as a CPU. Reads and executes various programs stored in the ROM, and stores the calculation result, which is the execution result, in the RAM or an external storage device. The details of each part of the control device 31 will be described below.

風速計測部12は、風速計8にて計測され入力I/F34及び内部バス37介して入力される風速から、風速の時系列データ及び所定時間で平均化処理を施した平均風速を算出する。そして、風速計測部12は、内部バス37を介して平均風速をウィンドシア推定部18へ転送すると共に、風速の時系列データを風速分布計算部19へ転送する。なお、平均化処理では、例えば10分ごとに平均風速を算出しても良く、移動平均を用いて平均値を連続的に出力しても良い。 The wind speed measuring unit 12 calculates the time-series data of the wind speed and the average wind speed averaged at a predetermined time from the wind speed measured by the anemometer 8 and input via the input I / F 34 and the internal bus 37. Then, the wind speed measuring unit 12 transfers the average wind speed to the wind shear estimation unit 18 via the internal bus 37, and transfers the time series data of the wind speed to the wind speed distribution calculation unit 19. In the averaging process, for example, the average wind speed may be calculated every 10 minutes, or the average value may be continuously output using a moving average.

荷重計測部13は、ひずみセンサ7にて計測されるタワー21のひずみを、入力I/F34及び内部バス37介して入力し、入力されたタワー21のひずみの計測値に対し、計測位置での断面係数を乗ずることで、曲げモーメントに変換する。更に、曲げモーメントを所定時間での平均化処理を施して、ウィンドシア推定部18へ内部バス37を介して転送する。風力発電装置2の建設位置での風向が概ね一方向である場合には、ひずみセンサ7の設置は一箇所でもかまわないが、風向が変化する場合には、ナセル22の方向に沿って風力発電装置2を転倒させる向きの曲げモーメント成分を算出するために、タワー21の同一高さ地点に二箇所以上ひずみセンサ7を設置することが望ましい。また、ひずみセンサ7を取り付ける高さ方向の位置は、タワー21の頂部付近であることが望ましいが、直接頂部のモーメントを計測しなくとも、風速とロータ25のスラスト係数から補正を施すこともできるため、計測位置を頂部に限定する必要はない。 The load measuring unit 13 inputs the strain of the tower 21 measured by the strain sensor 7 via the input I / F 34 and the internal bus 37, and the input measured value of the strain of the tower 21 is measured at the measurement position. By multiplying by the section modulus, it is converted into a bending moment. Further, the bending moment is averaged in a predetermined time and transferred to the wind shear estimation unit 18 via the internal bus 37. If the wind direction at the construction position of the wind power generation device 2 is generally unidirectional, the strain sensor 7 may be installed at one place, but if the wind direction changes, wind power is generated along the direction of the nacelle 22. It is desirable to install strain sensors 7 at two or more locations at the same height point of the tower 21 in order to calculate the bending moment component in the direction in which the device 2 is overturned. Further, the position in the height direction in which the strain sensor 7 is attached is preferably near the top of the tower 21, but correction can be made from the wind speed and the thrust coefficient of the rotor 25 without directly measuring the moment at the top. Therefore, it is not necessary to limit the measurement position to the top.

温度計測部14は、温度計10にて計測された温度を、入力I/F34及び内部バス37介して入力し、所定時間での平均化温度を算出する。温度計測部14は、算出した平均化温度を大気密度計算部17へ内部バス37を介して転送する。気圧計測部15は、気圧計11にて計測された気圧を、入力I/F34及び内部バス37介して入力し、所定時間での平均化気圧を算出する。気圧計測部15は、算出した平均化気圧を大気密度計算部17へ内部バス37を介して転送する。
大気密度計算部17は、温度と気圧から気体の状態方程式を用いて、大気密度を計算し、例えば、10分ごとに大気密度の平均値をウィンドシア推定部18へ内部バス37を介して転送する。なお、大気密度を算出できれば、他の方法を用いてもかまわない。
The temperature measuring unit 14 inputs the temperature measured by the thermometer 10 via the input I / F 34 and the internal bus 37, and calculates the average temperature in a predetermined time. The temperature measuring unit 14 transfers the calculated averaged temperature to the atmospheric density calculation unit 17 via the internal bus 37. The barometric pressure measuring unit 15 inputs the barometric pressure measured by the barometer 11 via the input I / F 34 and the internal bus 37, and calculates the average barometric pressure at a predetermined time. The atmospheric pressure measuring unit 15 transfers the calculated averaged atmospheric pressure to the atmospheric density calculation unit 17 via the internal bus 37.
The atmospheric density calculation unit 17 calculates the atmospheric density from the temperature and the atmospheric pressure using the gas state equation, and transfers the average value of the atmospheric density to the Windsia estimation unit 18 via the internal bus 37 every 10 minutes, for example. do. If the atmospheric density can be calculated, another method may be used.

ウィンドシア推定部18は、荷重計測部13より転送される曲げモーメント、風速計測部12より転送される平均風速、及び、大気密度計算部17より転送される大気密度の平均値を入力として、これらの値とウィンドシアの強弱を表す冪指数αWSとの関係を定義するウィンドシア関数33を用いることにより、平均化に用いた所定時間ごとに冪指数αWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、曲げモーメント、平均風速、大気密度の平均値を用いているが、冪指数αWSの推定精度低下を許容できる場合には、荷重計測部13より転送される曲げモーメントのみを入力としても良く、曲げモーメントに加えて平均風速または大気密度の平均値の一方を用いても良い。また、平均風速に代えて、風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。 The wind shear estimation unit 18 inputs the bending moment transferred from the load measurement unit 13, the average wind speed transferred from the wind speed measurement unit 12, and the average value of the air density transferred from the air density calculation unit 17 as inputs. By using the wind shear function 33 that defines the relationship between the value of wind shear and the power index α WS representing the strength of the wind shear, the estimated value of the power index α WS is output at predetermined time intervals used for averaging. In this embodiment, the average values of bending moment, average wind speed, and atmospheric density are used as inputs to the wind shear estimation unit 18, but if a decrease in the estimation accuracy of the exponential α WS can be tolerated, the load measurement unit Only the bending moment transferred from 13 may be used as an input, and either the average wind speed or the average value of the atmospheric density may be used in addition to the bending moment. Further, instead of the average wind speed, another physical quantity that changes in response to the wind speed may be used. For example, the amount of power generation, the pitch angle, the rotor rotation speed, and the like can be used.

記憶部16に格納されるウィンドシア関数33は、曲げモーメント、平均風速、大気密度の平均値に対して冪指数αWSを一意に出力できる関数を保存している。関数の保存方法は特に限定されず、応答曲面やニューラルネットワークのほか、多項式のような数式でも、各変数をビンで区切ったデータテーブルでも構わない。関数の作成方法についても特に限定されないが、例えば、風力発電装置2を構成するナセル22の上部に一時的にドップラーライダを取り付けて、計測した風速分布から冪指数αWSを最小二乗法などで計算し、曲げモーメント、平均風速、大気密度の平均値を設計変数、冪指数αWSを目的関数として応答曲面を作成する方法がある。他の方法としては、数値シミュレーションを用いて、冪指数αWS、平均風速、大気密度の平均値を変化させた場合の曲げモーメントを計算し、計算結果をニューラルネットワークに学習させて、冪指数αWSを推定することもできる。 The wind shear function 33 stored in the storage unit 16 stores a function that can uniquely output the index α WS with respect to the average value of the bending moment, the average wind speed, and the atmospheric density. The storage method of the function is not particularly limited, and may be a response surface, a neural network, a mathematical formula such as a polynomial, or a data table in which each variable is separated by a bin. The method of creating the function is also not particularly limited, but for example, a Doppler lidar is temporarily attached to the upper part of the nacelle 22 constituting the wind power generator 2, and the exponentiation α WS is calculated from the measured wind speed distribution by the least squares method or the like. Then, there is a method of creating a response curved surface by using the average values of bending moment, average wind speed, and atmospheric density as design variables and exponentiation α WS as an objective function. Another method is to use a numerical simulation to calculate the bending moment when the average value of the exponent α WS , mean wind velocity, and atmospheric density is changed, and let the neural network learn the calculation result to learn the exponent α. WS can also be estimated.

風速分布計算部19は、ウィンドシア推定部18で得られた冪指数αWSの推定値、及び、風速計測部12より内部バス37を介して転送された風速の時系列データから、上記式(1)を用いて風速分布の時系列データを算出する。なお、風速計8の地表からの高さzrefは既知であるとする。
出力I/F35は、風速分布計算部19により算出され、内部バス37を介して転送された風速分布の時系列データを図示しない表示部に出力する。
通信I/F36は、風速分布計算部19により算出され、内部バス37を介して転送された風速分布の時系列データを、通信ネットワーク5を介して運転管理センター3内に設置される電子端末4或いは図示しないサーバへ送信する。
The wind speed distribution calculation unit 19 uses the above equation (from the estimated value of the wind shear index α WS obtained by the wind shear estimation unit 18 and the time series data of the wind speed transferred from the wind speed measurement unit 12 via the internal bus 37. Calculate the time series data of the wind speed distribution using 1). It is assumed that the height zref of the anemometer 8 from the ground surface is known.
The output I / F 35 is calculated by the wind speed distribution calculation unit 19, and outputs the time series data of the wind speed distribution transferred via the internal bus 37 to a display unit (not shown).
The communication I / F 36 is an electronic terminal 4 installed in the operation management center 3 via the communication network 5 for the time series data of the wind speed distribution calculated by the wind speed distribution calculation unit 19 and transferred via the internal bus 37. Alternatively, it is sent to a server (not shown).

次に、風況推定装置32によるウィンドシアの推定原理について説明する。ロータ25に流入する風速がウィンドシアにより、上記式(1)のような分布をとると仮定する。ブレード24に作用する風向方向の空気力(スラスト力)Tは式(2)で表すことができる。 Next, the principle of wind shear estimation by the wind condition estimation device 32 will be described. It is assumed that the wind speed flowing into the rotor 25 has a distribution as shown in the above equation (1) due to the wind shear. The aerodynamic force (thrust force) T acting on the blade 24 in the wind direction can be expressed by the equation (2).

Figure 0007009237000002
Figure 0007009237000002

ここで、bはブレード24の長さ、ρは大気密度、Cはブレード24のハブ23の中心から距離rの位置における局所スラスト係数、Vはブレード24のハブ23の中心から距離rの位置における局所風速、cはブレード24のハブ23の中心から距離rの位置における翼弦長である。これらのうち、b及びcはブレード24に固有の値であり、運転条件によらず一定である。また、Cもブレード24に固有の値であるが、ピッチ角、風速、回転数に応じて変化する。ここで簡単のため、風速は高さ方向に上述の式(1)の分布を有するが時間変化せず、ピッチ角、回転数は一定であると仮定する。また、ロータ25はチルト角、コーニング角を持たず鉛直平面内を回転すると仮定する。このとき、ブレード24が失速しない範囲では風速の増加に伴ってブレード24の迎え角が増加するため、Cも単調増加する。したがって、上記式(2)からわかるように、局所風速Vが増加する場合、Cが増加することによってスラスト力Tも増加する。 Here, b is the length of the blade 24, ρ is the air density, Ct is the local thrust coefficient at the position of the distance r from the center of the hub 23 of the blade 24, and V is the position of the distance r from the center of the hub 23 of the blade 24. The local wind speed, c, is the chord length at a distance r from the center of the hub 23 of the blade 24. Of these, b and c are values unique to the blade 24 and are constant regardless of operating conditions. Further, C t is also a value peculiar to the blade 24, but changes depending on the pitch angle, the wind speed, and the rotation speed. Here, for the sake of simplicity, it is assumed that the wind speed has the distribution of the above equation (1) in the height direction, but does not change with time, and the pitch angle and the rotation speed are constant. Further, it is assumed that the rotor 25 does not have a tilt angle and a cornering angle and rotates in a vertical plane. At this time, as long as the blade 24 does not stall, the angle of attack of the blade 24 increases as the wind speed increases, so that Ct also increases monotonically. Therefore, as can be seen from the above equation (2), when the local wind speed V increases, the thrust force T also increases as the C t V 2 increases.

今、大気密度及びハブ23の中心での風速が一定であると仮定すると、図3においてzrefがハブ23の高さを表し、冪指数αWSの増加に伴ってハブ23の上方では風速が増加し、ハブ23の下方では風速が低下することがわかる。従って、図5に示すようにスラスト力Tも冪指数αWSの増加に伴って、ハブ23の上方では増加し、ハブ23の下方では低下する。ここで、ハブ23の中心でのナセル22の方向に沿って風力発電装置2を転倒させる向きのモーメントを考えると、ブレード24が回転平面内の任意の位置にあるとき、冪指数αWSの増加に伴ってモーメントが単調増加することがわかる。ゆえに、ロータ25が一回転する間のハブ23の中心でのモーメントの平均値も冪指数αWSの増加に伴って図6に示すように単調増加する。従って、大気密度及びハブ23の中心での風速が一定である場合には、ハブ23の中心でのモーメントから冪指数αWSを推定することができる。なお、図6に示す冪指数とハブ中心でのモーメントの関係を説明するための図は、一般的な風力発電装置のモデルに対しシミュレーションにより生成したものである。実際の風力発電装置2では、ハブ23の中心での風速と大気密度も変化するため、これらの増加によってもモーメントは単調増加する(図7、図8)。よって、冪指数αWSを高精度に推定するためには、モーメントに加えて、ある位置での風速と大気密度を計測する必要がある。風速と大気密度がわかれば、図7及び図8に示すような、複数のモーメント-冪指数αWS関係式の中から、風速と大気密度に合った関係式を選択し、モーメントから冪指数αWSを推定することができる。ただし、図8に示す大気密度の変化によるモーメントの変動は、図7に示す冪指数αWSと風速の変化と比較すると小さいため、大気密度を用いなくとも冪指数αWSを推定することは可能である。また、風速の変動幅が小さい場合や冪指数αWSの推定精度低下を許容できる場合には、風速と大気密度を使用せず、モーメントのみから冪指数αWSを推定しても良い。 Assuming that the air density and the wind speed at the center of the hub 23 are constant, zref represents the height of the hub 23 in FIG. 3, and the wind speed is higher above the hub 23 as the exponent α WS increases. It can be seen that the wind speed increases and the wind speed decreases below the hub 23. Therefore, as shown in FIG. 5, the thrust force T also increases above the hub 23 and decreases below the hub 23 as the exponent α WS increases. Here, considering the moment in the direction of overturning the wind turbine generator 2 along the direction of the nacelle 22 at the center of the hub 23, the exponentiation α WS increases when the blade 24 is at an arbitrary position in the rotation plane. It can be seen that the moment increases monotonically with this. Therefore, the average value of the moments at the center of the hub 23 during one rotation of the rotor 25 also monotonically increases as the exponent α WS increases, as shown in FIG. Therefore, when the atmospheric density and the wind speed at the center of the hub 23 are constant, the exponent α WS can be estimated from the moment at the center of the hub 23. The figure for explaining the relationship between the exponent and the moment at the center of the hub shown in FIG. 6 is generated by simulation with respect to a model of a general wind power generation device. In the actual wind power generation device 2, the wind speed and the atmospheric density at the center of the hub 23 also change, and the moment increases monotonically due to these increases (FIGS. 7 and 8). Therefore, in order to estimate the exponential α WS with high accuracy, it is necessary to measure the wind speed and atmospheric density at a certain position in addition to the moment. Once the wind speed and air density are known, select a relational expression suitable for the wind speed and air density from a plurality of moment-exponentiation α WS relational expressions as shown in FIGS. 7 and 8, and the exponentiation α from the moment. WS can be estimated. However, since the change in moment due to the change in air density shown in FIG. 8 is smaller than the change in power index α WS and wind speed shown in FIG. 7, it is possible to estimate the power index α WS without using the air density. Is. Further, when the fluctuation range of the wind speed is small or when the estimation accuracy of the exponentiation α WS can be reduced, the exponentiation α WS may be estimated only from the moment without using the wind speed and the atmospheric density.

実際には回転するハブ23でのモーメント計測は困難であるため、ひずみセンサ7によって計測したタワー21でのモーメントで代用する、タワー21の高さ方向二箇所で計測したモーメントを外挿してハブ23の中心でのモーメントを算出するなどしても良い。風速と大気密度については、ハブ23の中心の値を用いる必要はなく、風力発電装置2の近傍で計測すれば良い。また、ピッチ角、回転数は一定と仮定したが、これらは風速に応じた制御によって決定されることが一般的であるため、風速を用いて冪指数αWSを推定する場合には、ピッチ角、回転数が変化しても推定精度への影響は小さい。チルト角、コーニング角は、風力発電装置2に固有の値であり、運転条件によらず一定であるため、推定精度への影響はない。 Since it is actually difficult to measure the moment with the rotating hub 23, the moment measured at two points in the height direction of the tower 21 is used as a substitute for the moment measured by the strain sensor 7, and the hub 23 is externally mounted. You may calculate the moment at the center of. It is not necessary to use the values at the center of the hub 23 for the wind speed and the atmospheric density, and the values may be measured in the vicinity of the wind power generation device 2. In addition, although the pitch angle and the number of revolutions are assumed to be constant, these are generally determined by control according to the wind speed, so when estimating the exponent α WS using the wind speed, the pitch angle Even if the rotation speed changes, the effect on the estimation accuracy is small. Since the tilt angle and the cornering angle are values peculiar to the wind power generation device 2 and are constant regardless of the operating conditions, there is no influence on the estimation accuracy.

なお、本実施例では図4に示したように、風況推定装置32を構成する、ウィンドシア関数33を格納する記憶部16、ウィンドシア推定部18、及び風速分布計算部19を制御装置31内に実装する場合を一例として説明したがこれに限られるものではない。例えば、これらウィンドシア関数33を格納する記憶部16、ウィンドシア推定部18、及び風速分布計算部19を、図1に示した運転管理センター3内に設置される電子端末4或いは図示しないサーバに実装する構成としても良い。 In this embodiment, as shown in FIG. 4, the control device 31 is a storage unit 16 for storing the wind shear function 33, a wind shear estimation unit 18, and a wind speed distribution calculation unit 19 that constitute the wind condition estimation device 32. The case of implementing within is described as an example, but the present invention is not limited to this. For example, the storage unit 16 storing the wind shear function 33, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 can be installed in the electronic terminal 4 installed in the operation management center 3 shown in FIG. 1 or a server (not shown). It may be a configuration to be implemented.

以上の通り、本実施例によれば、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供することが可能となる。 As described above, according to the present embodiment, it is possible to provide a wind power generation device and a wind power generation system capable of estimating wind shear, which is a wind speed distribution, with high accuracy with a simple configuration.

図9は、本発明の他の実施例に係る実施例2の風力発電装置を構成する制御装置31aの機能ブロック図である。本実施例では、実施例1で示したウィンドシア推定部18及び風速分布計算部19にて推定されたウィンドシアを用いて風力発電装置2の信頼性を評価する信頼性評価装置40を制御装置内に設けた点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。なお、説明の便宜上、図9に示すように、実施例1と同様の構成要素として、入力I/F34、ウィンドシア推定部18、出力I/F35、通信I/F36、及び内部バス37のみを示す。 FIG. 9 is a functional block diagram of the control device 31a constituting the wind power generation device of the second embodiment according to another embodiment of the present invention. In this embodiment, the control device controls the reliability evaluation device 40 for evaluating the reliability of the wind power generation device 2 using the wind shear estimated by the wind shear estimation unit 18 and the wind speed distribution calculation unit 19 shown in the first embodiment. The point provided inside is different from the first embodiment. The same components as those in the first embodiment are designated by the same reference numerals, and the description overlapping with the first embodiment will be omitted below. For convenience of explanation, as shown in FIG. 9, only the input I / F34, the wind shear estimation unit 18, the output I / F35, the communication I / F36, and the internal bus 37 are used as the same components as in the first embodiment. show.

図9に示すように、本実施例に係る制御装置31aは、入力I/F34、ウィンドシア推定部18、出力I/F35、通信I/F36、に加え、運転条件取得部41、荷重計算部42、設計情報43を格納する記憶部16a、信頼性評価部44、及び情報出力部45を備え、これらは相互に内部バス37にてアクセス可能に接続されている。運転条件取得部41、荷重計算部42、設計情報43を格納する記憶部16a、信頼性評価部44、及び情報出力部45から信頼性評価装置40が構成されている。なお、信頼性評価装置40が情報出力部45を有しない構成としても良い。ウィンドシア推定部18、運転条件取得部41、荷重計算部42、信頼性評価部44、及び情報出力部45は、例えば、図示しないCPU(Central Processing Unit)などのプロセッサ、各種プログラムを格納するROM、演算過程のデータを一時的に格納するRAM、外部記憶装置などの記憶装置にて実現されると共に、CPUなどのプロセッサがROMに格納された各種プログラムを読み出し実行し、実行結果である演算結果をRAM又は外部記憶装置に格納する。 As shown in FIG. 9, the control device 31a according to the present embodiment has an input I / F34, a wind shear estimation unit 18, an output I / F35, a communication I / F36, an operating condition acquisition unit 41, and a load calculation unit. 42, a storage unit 16a for storing design information 43, a reliability evaluation unit 44, and an information output unit 45 are provided, and these are connected to each other so as to be accessible by an internal bus 37. The reliability evaluation device 40 is composed of an operating condition acquisition unit 41, a load calculation unit 42, a storage unit 16a for storing design information 43, a reliability evaluation unit 44, and an information output unit 45. The reliability evaluation device 40 may be configured not to have the information output unit 45. The wind shear estimation unit 18, the operating condition acquisition unit 41, the load calculation unit 42, the reliability evaluation unit 44, and the information output unit 45 are, for example, a ROM for storing a processor such as a CPU (Central Processing Unit) (not shown) and various programs. , A RAM that temporarily stores data in the calculation process, a storage device such as an external storage device, and a processor such as a CPU that reads and executes various programs stored in the ROM, and is the calculation result that is the execution result. Is stored in a RAM or an external storage device.

以下では、風力発電装置2のブレード24の信頼性評価を一例として、信頼性評価装置40について説明する。なお、信頼性評価装置40の適用先は、ブレード24に限定される必要はなく、ウィンドシアが信頼性に影響を与える部位であれば、ナセル22やタワー21など風力発電装置2の他の構成部品でも良い。 Hereinafter, the reliability evaluation device 40 will be described by taking the reliability evaluation of the blade 24 of the wind power generation device 2 as an example. The application destination of the reliability evaluation device 40 does not have to be limited to the blade 24, and if the wind shear affects the reliability, other configurations of the wind power generation device 2 such as the nacelle 22 and the tower 21 It may be a part.

運転条件取得部41は、ブレード24に作用する荷重に関連する風力発電装置2の運転条件の時刻歴データを取得する。運転条件とは、例えば、風速計8での風速や風向、ブレード24のピッチ角やアジマス角、ロータ25の回転速度、ナセル22の方位角、風力発電装置2の発電量などである。風力発電装置2にひずみセンサや加速度センサなどの荷重センサが取り付けられている場合には、それらの時刻歴データを含んでも良い。また、制御装置31aとして、例えばSCADAが用いられる場合には、このSCADAから運転条件を取得しても良い。 The operating condition acquisition unit 41 acquires time history data of the operating conditions of the wind power generation device 2 related to the load acting on the blade 24. The operating conditions include, for example, the wind speed and direction of the anemometer 8, the pitch angle and azimuth angle of the blade 24, the rotation speed of the rotor 25, the azimuth angle of the nacelle 22, and the amount of power generated by the wind power generation device 2. When a load sensor such as a strain sensor or an acceleration sensor is attached to the wind power generator 2, the time history data thereof may be included. Further, when, for example, SCADA is used as the control device 31a, the operating conditions may be acquired from this SCADA.

荷重計算部42は、運転条件取得部41から転送される運転条件、ウィンドシア推定部18から転送される冪指数αWSと、ブレード24の設計情報43とを用いて、ブレード24に作用する荷重の時刻歴を計算する。荷重の計算方法としては、例えば、風速計8での風速と冪指数αWSから上述の式(1)を用いてロータ25に流入する風速分布(ウィンドシア)を計算し、風速分布と運転条件の時刻歴データを翼素運動量理論或いはマルチボディダイナミクスなどに基づく空力弾性シミュレーションの入力とする方法がある。
記憶部16aに格納される設計情報43は、例えば、ブレード24とタワー21の寸法や質量分布、剛性分布、空力係数、ロータ25のチルト角やコーニング角、ナセル22の寸法や質量分布、空力係数、風力発電装置2の制御プログラムなどの設計データである。なお、記憶部16aに格納される設計情報25には、信頼性情報としてのブレード24の構成部品の信頼性に関するデータ、例えば、構成部品の寸法、弾性係数、断面係数、応力集中係数、S―N線図などが含まれる。
The load calculation unit 42 uses the operating conditions transferred from the operating condition acquisition unit 41, the exponentiation α WS transferred from the wind shear estimation unit 18, and the design information 43 of the blade 24 to act on the blade 24. Calculate the time history of. As a method of calculating the load, for example, the wind speed distribution (wind shear) flowing into the rotor 25 is calculated from the wind speed in the anemometer 8 and the aerodynamic index α WS using the above equation (1), and the wind speed distribution and operating conditions are used. There is a method of inputting the time history data of the above into an aeroelastic simulation based on the wind shear momentum theory or multibody dynamics.
The design information 43 stored in the storage unit 16a is, for example, the dimensions and mass distribution of the blade 24 and the tower 21, the rigidity distribution, the aerodynamic coefficient, the tilt angle and cornering angle of the rotor 25, the dimensions and mass distribution of the nacelle 22, and the aerodynamic coefficient. , Design data such as a control program for the wind power generator 2. The design information 25 stored in the storage unit 16a includes data on the reliability of the components of the blade 24 as reliability information, for example, the dimensions of the components, the elastic modulus, the section modulus, the stress concentration coefficient, and S-. N-ray diagram and the like are included.

信頼性評価部44は、荷重計算部42から出力されるブレード24に作用する荷重の時刻歴と上記信頼性情報とを用いて、ブレード24の構成部品の信頼性評価を行う。信頼性評価としては、構成部品に対する疲労損傷度や余寿命、破壊確率などの計算を行う。例えば、荷重の時刻歴データから疲労損傷度を計算する場合には、以下のような方法がある。まず、ブレード24に作用する荷重から構成部品に作用する応力の時刻歴を計算する。次に、応力の時刻歴データにレインフロー法を適用して応力振幅の出現頻度分布に変換し、得られた応力振幅の出現頻度分布と信頼性情報に保存される構成部品の材料のS―N線図から、線形累積損傷度則を用いて、時刻歴における疲労損傷度を計算する。 The reliability evaluation unit 44 evaluates the reliability of the components of the blade 24 by using the time history of the load acting on the blade 24 output from the load calculation unit 42 and the reliability information. For reliability evaluation, the degree of fatigue damage, remaining life, and fracture probability of components are calculated. For example, when calculating the degree of fatigue damage from the time history data of the load, there are the following methods. First, the time history of the stress acting on the component is calculated from the load acting on the blade 24. Next, the rainflow method is applied to the stress time history data to convert it into the appearance frequency distribution of the stress amplitude, and the S- of the component material stored in the obtained stress amplitude appearance frequency distribution and reliability information. From the N diagram, the degree of fatigue damage in the time history is calculated using the linear cumulative damage degree rule.

情報出力部45は、信頼性評価部44から内部バス37を介して転送される評価結果を表やグラフ、コンター図として表示する。なお、運転条件取得部41、ウィンドシア推定部18、荷重計算部42から出力されるデータを表示してもよい。例えば、運転条件取得部41から出力される風速計8での風速と、ウィンドシア推定部18から出力される冪指数αWSを用いて上述の式(1)により、ロータ25に流入する風速分布を表示しても良い。
なお、本実施例では、信頼性評価装置40を、制御装置31a内に設ける構成としたが、これに限られず、図1に示した運転管理センター3内に設置される電子端末4或いは図示しないサーバに実装する構成としても良い。
The information output unit 45 displays the evaluation results transferred from the reliability evaluation unit 44 via the internal bus 37 as a table, a graph, or a contour diagram. The data output from the operating condition acquisition unit 41, the wind shear estimation unit 18, and the load calculation unit 42 may be displayed. For example, the wind speed distribution in the rotor 25 by the above equation (1) using the wind speed in the anemometer 8 output from the operating condition acquisition unit 41 and the exponentiation α WS output from the wind shear estimation unit 18. May be displayed.
In this embodiment, the reliability evaluation device 40 is provided in the control device 31a, but the present invention is not limited to this, and the electronic terminal 4 installed in the operation control center 3 shown in FIG. 1 or not shown. It may be configured to be implemented on the server.

以上の通り本実施例によれば、実施例1の効果に加え、信頼性評価装置40を用いることで、風力発電装置の運転条件と、風力発電装置の荷重から推定したウィンドシアの時刻歴データに基づいて、風力発電装置の信頼性を評価することが可能となる。 As described above, according to the present embodiment, in addition to the effect of the first embodiment, by using the reliability evaluation device 40, the time history data of Windsia estimated from the operating conditions of the wind power generation device and the load of the wind power generation device. Based on the above, it is possible to evaluate the reliability of the wind power generation device.

図10は、本発明の他の実施例に係る実施例3の風力発電装置を構成する制御装置31bの機能ブロック図である。本実施例では、実施例1で示したウィンドシア推定部18及び風速分布計算部19にて推定されたウィンドシアを用いて風力発電装置2を制御する点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。なお、説明の便宜上、図10に示すように、実施例1と同様の構成要素として、入力I/F34、ウィンドシア推定部18、出力I/F35、通信I/F36、及び内部バス37のみを示す。 FIG. 10 is a functional block diagram of a control device 31b constituting the wind power generation device of the third embodiment according to another embodiment of the present invention. This embodiment is different from the first embodiment in that the wind shear generator 2 is controlled by using the wind shear estimated by the wind shear estimation unit 18 and the wind speed distribution calculation unit 19 shown in the first embodiment. The same components as those in the first embodiment are designated by the same reference numerals, and the description overlapping with the first embodiment will be omitted below. For convenience of explanation, as shown in FIG. 10, only the input I / F34, the wind shear estimation unit 18, the output I / F35, the communication I / F36, and the internal bus 37 are used as the same components as in the first embodiment. show.

図10に示すように本実施例に係る制御装置31bは、入力I/F34、ウィンドシア推定部18、出力I/F35、通信I/F36、に加え、運転条件取得部41、制御量計算部51、制御情報52を格納する記憶部16bを備え、これらは相互に内部バス37にてアクセス可能に接続されている。ウィンドシア推定部18、運転条件取得部41、及び制御量計算部51は、例えば、図示しないCPU(Central Processing Unit)などのプロセッサ、各種プログラムを格納するROM、演算過程のデータを一時的に格納するRAM、外部記憶装置などの記憶装置にて実現されると共に、CPUなどのプロセッサがROMに格納された各種プログラムを読み出し実行し、実行結果である演算結果をRAM又は外部記憶装置に格納する。 As shown in FIG. 10, the control device 31b according to the present embodiment has an input I / F34, a wind shear estimation unit 18, an output I / F35, a communication I / F36, an operating condition acquisition unit 41, and a control amount calculation unit. 51, a storage unit 16b for storing control information 52 is provided, and these are connected to each other so as to be accessible by an internal bus 37. The wind shear estimation unit 18, the operating condition acquisition unit 41, and the control amount calculation unit 51 temporarily store, for example, a processor such as a CPU (Central Processing Unit) (not shown), a ROM for storing various programs, and data of an arithmetic process. It is realized by a storage device such as a RAM or an external storage device, and a processor such as a CPU reads and executes various programs stored in the ROM, and stores the calculation result which is the execution result in the RAM or the external storage device.

運転条件取得部41は、上述の実施例2と同様に風力発電装置2の運転条件を取得する。制御量計算部51は、運転条件取得部41から内部バス37を介して転送される運転条件と、ウィンドシア推定部18から内部バス37を介して転送される冪指数αWSと、制御情報52とを用いて、例えばブレード24のピッチ角やロータ25の回転速度などを決定し、これらを実現するように制御量を決定する。ブレード24のピッチ角やロータ25の回転速度は、例えば、発電量を最大化する、ブレード24の荷重変動を最小化する、というようにして決定できる。なお、ピッチ角はロータ25が一回転する間に変化しても良い。このような制御は、風速計8での風速と冪指数αWSから上述の式(1)を用いてロータ4に流入する風速分布を計算し、任意のロータ回転速度に対してブレード24に流入する風速、風向の分布を高精度に推定することで実現できる。
記憶部16bは、制御情報52として、制御に用いる定数やブレード24の空力特性を格納している。制御量計算部51から出力される上述の制御量は、内部バス37及び出力I/F35を介して、発電機28やピッチ角制御機構6などへ出力される。
The operating condition acquisition unit 41 acquires the operating conditions of the wind power generation device 2 in the same manner as in the second embodiment described above. The control amount calculation unit 51 includes an operating condition transferred from the operating condition acquisition unit 41 via the internal bus 37, a wind shear estimation unit 18 transferred via the internal bus 37, and a control information 52. For example, the pitch angle of the blade 24 and the rotation speed of the rotor 25 are determined by using and, and the control amount is determined so as to realize these. The pitch angle of the blade 24 and the rotation speed of the rotor 25 can be determined, for example, by maximizing the amount of power generation and minimizing the load fluctuation of the blade 24. The pitch angle may change during one rotation of the rotor 25. In such control, the wind speed distribution flowing into the rotor 4 is calculated from the wind speed of the anemometer 8 and the index α WS using the above equation (1), and the wind speed flows into the blade 24 for an arbitrary rotor rotation speed. This can be achieved by estimating the distribution of wind speed and direction with high accuracy.
The storage unit 16b stores constants used for control and aerodynamic characteristics of the blade 24 as control information 52. The above-mentioned control amount output from the control amount calculation unit 51 is output to the generator 28, the pitch angle control mechanism 6, and the like via the internal bus 37 and the output I / F 35.

以上の通り本実施例によれは、実施例1の効果に加え、風力発電装置の運転条件と、風力発電装置の荷重から推定したウィンドシアの時刻歴データに基づいて、風力発電装置を制御し、発電量の最大化やブレードの荷重変動の最小化が可能となる。 As described above, according to the present embodiment, in addition to the effect of the first embodiment, the wind power generation device is controlled based on the operating conditions of the wind power generation device and the time history data of Windsia estimated from the load of the wind power generation device. , It is possible to maximize the amount of power generation and minimize the load fluctuation of the blade.

図11は、本発明の他の実施例に係る実施例4の風力発電装置の構成を示す図である。本実施例では、上述の実施例1の風力発電装置2の風況推定装置32を用いて、水平方向のウィンドシアを推定する点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。 FIG. 11 is a diagram showing a configuration of a wind power generation device according to a fourth embodiment according to another embodiment of the present invention. This embodiment differs from the first embodiment in that the wind shear in the horizontal direction is estimated by using the wind condition estimation device 32 of the wind power generation device 2 of the above-mentioned embodiment 1. The same components as those in the first embodiment are designated by the same reference numerals, and the description overlapping with the first embodiment will be omitted below.

図11では、風力発電装置2及びその周りでの水平方向の風速分布61の構成を示している。図11では、風力発電装置2を上方から眺めた状態を示しており、風は紙面左から右に吹いているものとする。なお、風力発電装置2は、上述の実施例1と同様の構成を備えている。 FIG. 11 shows the configuration of the wind power generation device 2 and the horizontal wind speed distribution 61 around it. FIG. 11 shows a state in which the wind power generation device 2 is viewed from above, and it is assumed that the wind is blowing from the left to the right of the paper. The wind power generation device 2 has the same configuration as that of the first embodiment described above.

図11に示す水平方向の風速分布61は、水平方向に変化しており、この水平方向のウィンドシアの強弱を表す係数をβWSとすると、水平方向の風速分布61を以下の式(3)のように仮定することができる。 The horizontal wind velocity distribution 61 shown in FIG. 11 changes in the horizontal direction, and if the coefficient representing the strength of the horizontal wind shear is β WS , the horizontal wind velocity distribution 61 is given by the following equation (3). Can be assumed as follows.

Figure 0007009237000003
Figure 0007009237000003

ここで、V(y)は水平面上で風速ベクトルに直交する向きの位置yにおける風速であり、yrefは基準となる風速を定義する位置、V(yref)は基準となる風速を表している。式(3)からわかるように、係数βWSが大きくなるほど、水平方向の風速の変化は大きくなる。例えば、基準となる風速を風速計8で計測する場合、式(3)のV(yref)を計測した風速、yrefを0(原点)として、yは図11において風速計8からのy軸方向の距離を用いることで、ある係数βWSがわかれば水平方向の風速分布61を得ることができる。すなわち、風速分布を式(3)のように仮定する場合、風速分布の推定問題は、係数βWSの推定問題に帰着する。なお、本実施例では、式(3)の風速分布を仮定して、水平方向のウィンドシアの強弱を表す係数βWSを推定することを考えるが、風速分布の仮定は式(3)に限定する必要はなく、複数のパラメータを用いて風速分布を仮定しても良い。 Here, V (y) is the wind speed at the position y in the direction orthogonal to the wind speed vector on the horizontal plane, y ref is the position that defines the reference wind speed, and V (y ref ) is the reference wind speed. There is. As can be seen from the equation (3), the larger the coefficient β WS , the larger the change in the horizontal wind speed. For example, when the reference wind speed is measured by the anemometer 8, y is the wind speed measured by V (y ref ) in the equation (3), y ref is 0 (origin), and y is y from the anemometer 8 in FIG. By using the axial distance, a horizontal wind speed distribution 61 can be obtained if a certain coefficient β WS is known. That is, when the wind speed distribution is assumed as in the equation (3), the problem of estimating the wind speed distribution results in the problem of estimating the coefficient β WS . In this embodiment, it is considered to estimate the coefficient β WS representing the strength of the wind shear in the horizontal direction by assuming the wind speed distribution in the equation (3), but the assumption of the wind speed distribution is limited to the equation (3). It is not necessary to do so, and the wind speed distribution may be assumed using a plurality of parameters.

水平方向のウィンドシアの強弱を定義する係数βWSを推定するための、本実施例における風力発電装置2の風況推定装置32の構成は実施例1と同様である。以下では、風況推定装置32の詳細について、実施例1と異なる点について述べる。
本実施例において、荷重計測部13(図示せず)では、ひずみセンサ7を用いてタワー21のせん断ひずみを計測し、計測位置での断面二次極モーメントとタワー21の半径を用いて、捻りモーメントに変換する。更に、捻りモーメントを所定時間での平均化処理を施して出力する。捻りモーメントはタワー21の断面内で均一に作用するので、ひずみセンサ7の設置は一箇所でもよく、高さ方向の位置も特に限定する必要はない。
The configuration of the wind condition estimation device 32 of the wind power generation device 2 in this embodiment for estimating the coefficient β WS that defines the strength of the wind shear in the horizontal direction is the same as that of the first embodiment. Hereinafter, the details of the wind condition estimation device 32 will be described as different from those of the first embodiment.
In this embodiment, the load measuring unit 13 (not shown) measures the shear strain of the tower 21 using the strain sensor 7, and twists using the moment of inertia of area at the measurement position and the radius of the tower 21. Convert to moment. Further, the twisting moment is averaged in a predetermined time and output. Since the twisting moment acts uniformly within the cross section of the tower 21, the strain sensor 7 may be installed at only one place, and the position in the height direction does not need to be particularly limited.

ウィンドシア推定部18(図示せず)では、捻りモーメント、風速、大気密度の平均値を入力として、これらの値と水平方向のウィンドシアの強弱を表す係数βWSとの関係を定義するウィンドシア関数33を用いることにより、平均化に用いた所定時間ごとに係数βWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、捻りモーメント、風速、大気密度を用いているが、係数βWSの推定精度低下を許容できる場合には、捻りモーメントのみを入力としても良く、捻りモーメントに加えて風速または大気密度の一方を用いても良い。また、風速の代わりに風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。なお、捻りモーメントは、ヨーエラーと呼ばれるナセル22と風向とのずれによっても発生するため、ヨーエラーをウィンドシア推定部18の入力に加えることで、高精度化が可能である。 The wind shear estimation unit 18 (not shown) takes the average values of the twisting moment, wind speed, and atmospheric density as inputs, and defines the relationship between these values and the coefficient β WS that represents the strength of the horizontal wind shear. By using the function 33, the estimated value of the coefficient β WS is output at predetermined time intervals used for averaging. In this embodiment, the twisting moment, the wind speed, and the atmospheric density are used as the inputs to the wind shear estimation unit 18, but if the estimation accuracy of the coefficient β WS can be reduced, only the twisting moment may be input. , Either wind speed or atmospheric density may be used in addition to the twisting moment. Further, instead of the wind speed, another physical quantity that changes in response to the wind speed may be used. For example, the amount of power generation, the pitch angle, the rotor rotation speed, and the like can be used. Since the twisting moment is also generated by the deviation between the nacelle 22 and the wind direction, which is called a yaw error, high accuracy can be achieved by adding the yaw error to the input of the wind shear estimation unit 18.

ウィンドシア関数33は、捻りモーメント、風速、大気密度の入力値に対して係数βWSを一意に出力できる関数として定義されている。関数の保存方法、作成方法は特に限定されず、成就の実施例1と同様の方法を用いても良い。
風速分布計算部19では、ウィンドシア推定部18で得られた係数βWSの推定値と風速計測部12で得られた風速の時系列データから、上記式(3)を用いて風速分布の時系列データを算出する。
The wind shear function 33 is defined as a function that can uniquely output the coefficient β WS with respect to the input values of the twisting moment, the wind speed, and the atmospheric density. The method of saving and creating the function is not particularly limited, and the same method as in Example 1 of fulfillment may be used.
In the wind speed distribution calculation unit 19, the wind speed distribution is performed using the above equation (3) from the estimated value of the coefficient β WS obtained by the wind shear estimation unit 18 and the time series data of the wind speed obtained by the wind speed measurement unit 12. Calculate series data.

風況推定装置32による水平方向のウィンドシアの推定原理は、以下の点を除いて実施例1と同様である。すなわち、水平方向のウィンドシアによってハブ23の中心に作用するモーメントは、ナセル22を水平面内で回転させる向きに発生する。したがって、タワー21には、捻りモーメントが作用するため、係数βWSの推定には、捻りモーメントを計測する必要がある。なお、タワー21の曲げモーメントと捻りモーメントを同時に計測することにより、実施例1と実施例4の高さ方向(Z方向)、水平方向(y方向)のウィンドシアを同時に推定しても良く、推定した水平方向のウィンドシアを、上述の実施例2の信頼性評価装置40、実施例3の制御装置31bの入力としても良い。 The principle of estimating the horizontal wind shear by the wind condition estimation device 32 is the same as that of the first embodiment except for the following points. That is, the moment acting on the center of the hub 23 by the horizontal wind shear is generated in the direction of rotating the nacelle 22 in the horizontal plane. Therefore, since the twisting moment acts on the tower 21, it is necessary to measure the twisting moment in order to estimate the coefficient β WS . By simultaneously measuring the bending moment and the twisting moment of the tower 21, the wind shears in the height direction (Z direction) and the horizontal direction (y direction) of Examples 1 and 4 may be estimated at the same time. The estimated horizontal wind shear may be used as an input for the reliability evaluation device 40 of Example 2 and the control device 31b of Example 3 described above.

以上の通り本実施例によれば、風力発電装置の荷重を計測することにより、簡易な構成にて、水平方向のウィンドシアを高精度に推定することが可能となる。 As described above, according to the present embodiment, by measuring the load of the wind power generation device, it is possible to estimate the wind shear in the horizontal direction with high accuracy with a simple configuration.

図12は、本発明の他の実施例に係る実施例5の風力発電装置の構成を示す図である。本実施例では、上述の実施例1に示したタワー21に設置されるひずみセンサ7をブレード24に設置した点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。 FIG. 12 is a diagram showing a configuration of a wind power generation device according to a fifth embodiment according to another embodiment of the present invention. This embodiment is different from the first embodiment in that the strain sensor 7 installed in the tower 21 shown in the first embodiment is installed on the blade 24. The same components as those in the first embodiment are designated by the same reference numerals, and the description overlapping with the first embodiment will be omitted below.

図12に本実施例に係る風力発電装置2及びその周りでの高さ方向の風速分布11の構成を示す。図12では、風力発電装置2を側方から眺めた状態を示しており、風は紙面左から右に吹いているものとする。風力発電装置2は、実施例1と同様の構成に加えて、タワー21に代えて、ブレード24にひずみセンサ7を備えている。ひずみセンサ7は、ロータ25を構成する少なくとも一つのブレード24に取り付けられていれば良く、ひずみセンサ7に代えて加速度センサのような他の荷重センサを用いても良い。 FIG. 12 shows the configuration of the wind power generation device 2 according to the present embodiment and the wind speed distribution 11 in the height direction around the wind power generator 2. FIG. 12 shows a state in which the wind power generation device 2 is viewed from the side, and it is assumed that the wind is blowing from the left to the right of the paper. In addition to the same configuration as in the first embodiment, the wind power generation device 2 includes a strain sensor 7 on the blade 24 instead of the tower 21. The strain sensor 7 may be attached to at least one blade 24 constituting the rotor 25, and another load sensor such as an acceleration sensor may be used instead of the strain sensor 7.

ウィンドシアの強弱を定義する冪指数αWSを推定するための、本実施例に係る風力発電装置2の風況推定装置32の構成は実施例1と同様である。以下では、風況推定装置32の詳細について、実施例1と異なる点について述べる。
本実施例において、荷重計測部13では、ひずみセンサ7を用いてブレード24のひずみを計測し、計測位置での断面係数を乗ずることで、曲げモーメントに変換する。更に、曲げモーメントを統計処理して統計値を出力する。統計値としては、例えば、所定時間での曲げモーメントの標準偏差や、最大値と最小値との差(最大振幅)などを用いる。また、アジマス角をビンで区切ってビンごとの平均値を用いてもよいし、ロータ25を構成する全てのブレード24にひずみセンサ7を取り付けて、各ブレード24の曲げモーメントからハブ23の中心でのモーメントを計算によって求めた値を用いても良い。ひずみセンサ7の設置は一箇所でもかまわないが、ブレード24を風向方向に曲げるモーメント成分を算出するために、ブレード24の長手方向の同一位置に二箇所以上ひずみセンサ7を設置することが望ましい。また、ひずみセンサ7はブレード24の翼根付近に取り付けることが望ましいが、計測位置を翼根付近に限定する必要はない。
The configuration of the wind condition estimation device 32 of the wind power generation device 2 according to the present embodiment for estimating the exponential α WS that defines the strength of the wind shear is the same as that of the first embodiment. Hereinafter, the details of the wind condition estimation device 32 will be described as different from those of the first embodiment.
In this embodiment, the load measuring unit 13 measures the strain of the blade 24 using the strain sensor 7 and converts it into a bending moment by multiplying it by the cross-sectional coefficient at the measurement position. Further, the bending moment is statistically processed and the statistical value is output. As the statistical value, for example, the standard deviation of the bending moment at a predetermined time, the difference between the maximum value and the minimum value (maximum amplitude), and the like are used. Further, the azimuth angle may be divided by bins and the average value for each bin may be used, or the strain sensors 7 may be attached to all the blades 24 constituting the rotor 25, and the bending moment of each blade 24 may be used at the center of the hub 23. You may use the value obtained by the calculation of the moment of. The strain sensor 7 may be installed at one place, but it is desirable to install the strain sensor 7 at two or more places at the same position in the longitudinal direction of the blade 24 in order to calculate the moment component that bends the blade 24 in the wind direction. Further, although it is desirable that the strain sensor 7 is attached near the blade root, it is not necessary to limit the measurement position to the vicinity of the blade root.

ウィンドシア推定部18は、ブレード24の曲げモーメントの統計量、風速、大気密度の平均値を入力として、これらの値とウィンドシアの強弱を表す冪指数αWSとの関係を定義するウィンドシア関数33を用いることにより、所定時間ごとに冪指数αWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、曲げモーメント、風速、大気密度を用いているが、冪指数αWSの推定精度低下を許容できる場合には、曲げモーメントのみを入力としても良く、曲げモーメントに加えて風速または大気密度の一方を用いても良い。また、風速に代えて風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。 The wind shear estimation unit 18 takes the statistics of the bending moment of the blade 24, the wind speed, and the average value of the air density as inputs, and defines the relationship between these values and the power index α WS representing the strength of the wind shear. By using 33, the estimated value of the wind shear index α WS is output at predetermined time intervals. In this embodiment, the bending moment, the wind speed, and the atmospheric density are used as the inputs to the wind shear estimation unit 18, but if the estimation accuracy of the exponent α WS can be reduced, only the bending moment may be input. At best, either wind speed or atmospheric density may be used in addition to the bending moment. Further, instead of the wind speed, another physical quantity that changes according to the wind speed may be used. For example, the amount of power generation, the pitch angle, the rotor rotation speed, and the like can be used.

風況推定装置32によるウィンドシアの推定原理は、以下の点を除いて実施例1と同様である。すなわち、ウィンドシアによってブレード24に作用するスラスト力は図5に示すように、アジマス角に応じて変化するため、ひずみセンサ7で計測される曲げモーメントもアジマス角に応じて変化する。従って、ロータ25が一回転する間にブレード24作用する曲げモーメントの変動幅を、標準偏差や最大振幅を用いて表すことで、冪指数αWSを推定できる。また、一回転での曲げモーメントの変動幅を用いなくても、あるアジマス角での曲げモーメントの値の大小から冪指数αWSを推定することもできる。例えば、図5に示すようにブレード24がハブ23の真上に位置するときの曲げモーメントやスラスト力は、冪指数αWSの増加に伴って単調増加するため、これらの値から冪指数αWSを推定できる。なお、ロータ25を構成する全てのブレード24にひずみセンサ7を取り付けて、各ブレード24の曲げモーメントからハブ23の中心でのモーメントを計算によって求める場合には、実施例1と同様の原理でウィンドシアを推定できる。 The principle of wind shear estimation by the wind condition estimation device 32 is the same as that of the first embodiment except for the following points. That is, as shown in FIG. 5, the thrust force acting on the blade 24 by the wind shear changes according to the azimuth angle, so that the bending moment measured by the strain sensor 7 also changes according to the azimuth angle. Therefore, the exponent α WS can be estimated by expressing the fluctuation range of the bending moment acting on the blade 24 during one rotation of the rotor 25 by using the standard deviation and the maximum amplitude. Further, it is also possible to estimate the exponent α WS from the magnitude of the value of the bending moment at a certain azimuth angle without using the fluctuation range of the bending moment in one rotation. For example, as shown in FIG. 5, the bending moment and the thrust force when the blade 24 is located directly above the hub 23 increase monotonically with the increase of the exponent α WS , and therefore, from these values, the exponent α WS Can be estimated. When the strain sensor 7 is attached to all the blades 24 constituting the rotor 25 and the moment at the center of the hub 23 is calculated from the bending moment of each blade 24, the window is based on the same principle as in the first embodiment. Shear can be estimated.

なお、ブレード24が水平となるときの曲げモーメントを用いれば、水平方向のウィンドシアを推定することもでき、高さ方向、水平方向のウィンドシアを同時に推定しても良い。また、推定したウィンドシアを、上述の実施例2の信頼性評価装置40、実施例3の制御装置31bの入力としても良い。 By using the bending moment when the blade 24 becomes horizontal, the wind shear in the horizontal direction can be estimated, and the wind shear in the height direction and the horizontal direction may be estimated at the same time. Further, the estimated wind shear may be used as an input for the reliability evaluation device 40 of the above-mentioned Example 2 and the control device 31b of the third embodiment.

以上の通り本実施例によっても、風力発電装置の荷重を計測することにより、簡易な構成にて、ウィンドシアを高精度に推定することが可能となる。 As described above, also in this embodiment, by measuring the load of the wind power generation device, it is possible to estimate the wind shear with high accuracy with a simple configuration.

なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。 The present invention is not limited to the above-described embodiment, and includes various modifications. For example, the above-described embodiment has been described in detail in order to explain the present invention in an easy-to-understand manner, and is not necessarily limited to the one including all the described configurations. Further, it is possible to replace a part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment.

1…風力発電システム
2…風力発電装置
3…運転管理センター
4…電子端末
5…通信ネットワーク
6…ピッチ角制御機構
7…ひずみセンサ
8…風速計
9…温度計
10…気圧計
11…風速分布
12…風速計測部
13…荷重計測部
14…温度計測部
15…気圧計測部
16,16a,16b…記憶部
17…大気密度計算部
18…ウィンドシア推定部
19…風速分布計算部
21…タワー
22…ナセル
23…ハブ
24…ブレード
25…ロータ
26…主軸
27…増速機
28…発電機
29…メインフレーム
30…電力変換器
31,31a,31b…制御装置
32…風況推定装置
33…ウィンドシア関数
34…入力I/F
35…出力I/F
36…通信I/F
37…内部バス
40…信頼性評価装置
41…運転条件取得部
42…荷重計算部
43…設計情報
44…信頼性評価部
45…情報出力部
51…制御量計算部
52…制御情報
61…水平方向の風速分布
1 ... Wind power generation system 2 ... Wind power generation device 3 ... Operation management center 4 ... Electronic terminal 5 ... Communication network 6 ... Pitch angle control mechanism 7 ... Strain sensor 8 ... Anemometer 9 ... Thermometer 10 ... Pressure gauge 11 ... Wind speed distribution 12 ... Wind speed measurement unit 13 ... Load measurement unit 14 ... Temperature measurement unit 15 ... Pressure measurement unit 16, 16a, 16b ... Storage unit 17 ... Atmospheric density calculation unit 18 ... Wind shear estimation unit 19 ... Wind speed distribution calculation unit 21 ... Tower 22 ... Nacelle 23 ... Hub 24 ... Blade 25 ... Rotor 26 ... Main shaft 27 ... Accelerator 28 ... Generator 29 ... Main frame 30 ... Power converters 31, 31a, 31b ... Control device 32 ... Wind condition estimation device 33 ... Wind shear function 34 ... Input I / F
35 ... Output I / F
36 ... Communication I / F
37 ... Internal bus 40 ... Reliability evaluation device 41 ... Operating condition acquisition unit 42 ... Load calculation unit 43 ... Design information 44 ... Reliability evaluation unit 45 ... Information output unit 51 ... Control amount calculation unit 52 ... Control information 61 ... Horizontal direction Wind speed distribution

Claims (15)

少なくともロータ及びナセル並びにナセルを回動可能に支持するタワーを有する風力発電装置と、風力発電装置を制御する制御装置とを備え、
前記制御装置は、
風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、前記荷重及び前記大気密度の平均値並びにウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電装置。
It comprises at least a wind power generator having a rotor and a nacelle and a tower that rotatably supports the nacelle, and a control device for controlling the wind power generator.
The control device is
A load measuring unit that measures the load applied to the wind power generator, a storage unit that stores the wind shear function that defines the relationship between the load and the wind shear, and an atmospheric density using a gas state equation from temperature and pressure. A wind shear estimation device having an atmospheric density calculation unit for calculating the average value of the atmospheric density and a wind shear estimation unit for calculating the wind shear based on the load, the average value of the atmospheric density, and the wind shear function. A wind shear device that is characterized by that.
請求項1に記載の風力発電装置において、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電装置。
In the wind power generation device according to claim 1,
The wind shear estimation unit is characterized in that a wind shear function stored in the storage unit and wind shear, which is a wind speed distribution in the height direction, are calculated based on the load input from the load measuring unit. Power generator.
請求項1に記載の風力発電装置において、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電装置。
In the wind power generation device according to claim 1,
The wind shear estimation unit calculates wind shear, which is a horizontal wind speed distribution, based on the wind shear function stored in the storage unit and the load input from the load measuring unit. Device.
請求項2に記載の風力発電装置において、
前記制御装置は、
前記高さ方向の風速分布であるウィンドシアを用いて、前記風力発電装置の信頼性を評価する信頼性評価装置を有することを特徴とする風力発電装置。
In the wind power generation device according to claim 2,
The control device is
A wind power generation device comprising a reliability evaluation device for evaluating the reliability of the wind power generation device using wind shear which is a wind speed distribution in the height direction.
請求項2又は請求項3に記載の風力発電装置において、
前記荷重計測部は、前記タワーに付加される荷重を計測することを特徴とする風力発電置。
In the wind power generation device according to claim 2 or claim 3.
The load measuring unit is a wind power generator that measures the load applied to the tower.
請求項2又は請求項3に記載の風力発電装置において、
前記荷重計測部は、前記風力発電装置に付加される荷重として、ひずみまたは加速度を計測することを特徴とする風力発電装置。
In the wind power generation device according to claim 2 or claim 3.
The load measuring unit is a wind power generation device that measures strain or acceleration as a load applied to the wind power generation device.
少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、
前記制御装置は、
風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、前記荷重及び前記大気密度の平均値並びにウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電システム。
It is provided with at least one wind power generation device, a control device for controlling the wind power generation device, an electronic terminal having a display device, and a communication network for connecting these to each other so as to be communicable.
The control device is
A load measuring unit that measures the load applied to the wind power generator, a storage unit that stores the wind shear function that defines the relationship between the load and the wind shear, and an atmospheric density using a gas state equation from temperature and pressure. A wind shear estimation device having an atmospheric density calculation unit for calculating the average value of the atmospheric density and a wind shear estimation unit for calculating the wind shear based on the load, the average value of the atmospheric density, and the wind shear function. A wind shear system characterized by that.
請求項7に記載の風力発電システムにおいて、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。
In the wind power generation system according to claim 7.
The wind shear estimation unit calculates wind shear, which is a wind speed distribution in the height direction, based on the wind shear function stored in the storage unit and the load input from the load measuring unit via the communication network. A wind power generation system characterized by that.
請求項7に記載の風力発電システムにおいて、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。
In the wind power generation system according to claim 7.
The wind shear estimation unit calculates wind shear, which is a horizontal wind speed distribution, based on the wind shear function stored in the storage unit and the load input from the load measuring unit via the communication network. A wind power generation system featuring.
請求項8に記載の風力発電システムにおいて、
前記制御装置は、
前記高さ方向の風速分布であるウィンドシアを用いて、前記風力発電装置の信頼性を評価する信頼性評価装置を有することを特徴とする風力発電システム。
In the wind power generation system according to claim 8,
The control device is
A wind power generation system comprising a reliability evaluation device for evaluating the reliability of the wind power generation device by using the wind shear which is the wind speed distribution in the height direction.
請求項8又は請求項9に記載の風力発電システムにおいて、
前記荷重計測部は、前記風力発電装置のタワーに付加される荷重を計測することを特徴とする風力発電システム。
In the wind power generation system according to claim 8 or 9.
The load measuring unit is a wind power generation system characterized by measuring a load applied to a tower of the wind power generation device.
請求項8又は請求項9に記載の風力発電システムにおいて、
前記荷重計測部は、前記風力発電装置に付加される荷重として、ひずみまたは加速度を計測することを特徴とする風力発電システム。
In the wind power generation system according to claim 8 or 9.
The load measuring unit is a wind power generation system characterized by measuring strain or acceleration as a load applied to the wind power generation device.
少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、
前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、温度と気圧から気体の状態方程式を用いて大気密度を計算し、大気密度の平均値を求める大気密度計算部と、を有し、
前記電子端末は、荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重及び前記大気密度計算部から入力される大気密度の平均値並びに前記記憶部に格納されるウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電システム。
It is provided with at least one wind power generation device, a control device for controlling the wind power generation device, an electronic terminal having a display device, and a communication network for connecting these to each other so as to be communicable.
The control device includes a load measuring unit that measures the load applied to the wind power generation device, and an atmospheric density calculation unit that calculates the atmospheric density from the temperature and pressure using the equation of state of the gas and obtains the average value of the atmospheric density. , Has,
The electronic terminal is input from a storage unit that stores a wind shear function that defines a relationship between a load and a wind shear, and a load input from the load measurement unit and an atmospheric density calculation unit via the communication network. A wind power generation system including a wind shear estimation device having a wind shear estimation unit that calculates wind shear based on an average value of atmospheric densities and a wind shear function stored in the storage unit.
請求項13に記載の風力発電システムにおいて、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。
In the wind power generation system according to claim 13,
The wind shear estimation unit calculates wind shear, which is a wind speed distribution in the height direction, based on the wind shear function stored in the storage unit and the load input from the load measuring unit via the communication network. A wind power generation system characterized by that.
請求項13に記載の風力発電システムにおいて、
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。
In the wind power generation system according to claim 13,
The wind shear estimation unit calculates wind shear, which is a horizontal wind speed distribution, based on the wind shear function stored in the storage unit and the load input from the load measuring unit via the communication network. A wind power generation system featuring.
JP2018014369A 2018-01-31 2018-01-31 Wind power generation equipment and wind power generation system Active JP7009237B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2018014369A JP7009237B2 (en) 2018-01-31 2018-01-31 Wind power generation equipment and wind power generation system
PCT/JP2018/046459 WO2019150805A1 (en) 2018-01-31 2018-12-18 Wind power generation device and wind power generation system
TW108102524A TWI729349B (en) 2018-01-31 2019-01-23 Wind power generation device and wind power generation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018014369A JP7009237B2 (en) 2018-01-31 2018-01-31 Wind power generation equipment and wind power generation system

Publications (2)

Publication Number Publication Date
JP2019132178A JP2019132178A (en) 2019-08-08
JP7009237B2 true JP7009237B2 (en) 2022-01-25

Family

ID=67477994

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018014369A Active JP7009237B2 (en) 2018-01-31 2018-01-31 Wind power generation equipment and wind power generation system

Country Status (3)

Country Link
JP (1) JP7009237B2 (en)
TW (1) TWI729349B (en)
WO (1) WO2019150805A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044754B (en) * 2019-12-09 2021-07-16 大唐山东烟台电力开发有限公司 Wind condition estimation device for wind generating set
FR3116123B1 (en) * 2020-11-06 2022-10-14 Ifp Energies Now Method for determining the wind speed in the plane of the rotor of a wind turbine
EP4067649A1 (en) * 2021-03-31 2022-10-05 Siemens Gamesa Renewable Energy A/S Operating a wind turbine in a severe weather condition

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100092292A1 (en) 2008-10-10 2010-04-15 Jacob Johannes Nies Apparatus and method for continuous pitching of a wind turbine
JP2015117682A (en) 2013-12-20 2015-06-25 三菱重工業株式会社 Wind generator monitoring system and wind generator monitoring method
WO2017163562A1 (en) 2016-03-25 2017-09-28 株式会社日立製作所 Remaining service life evaluating device and method, and wind power electricity generating system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100092292A1 (en) 2008-10-10 2010-04-15 Jacob Johannes Nies Apparatus and method for continuous pitching of a wind turbine
JP2015117682A (en) 2013-12-20 2015-06-25 三菱重工業株式会社 Wind generator monitoring system and wind generator monitoring method
WO2017163562A1 (en) 2016-03-25 2017-09-28 株式会社日立製作所 Remaining service life evaluating device and method, and wind power electricity generating system

Also Published As

Publication number Publication date
TWI729349B (en) 2021-06-01
TW201934870A (en) 2019-09-01
WO2019150805A1 (en) 2019-08-08
JP2019132178A (en) 2019-08-08

Similar Documents

Publication Publication Date Title
CN109219782B (en) System and method for controlling dynamic system
EP3613982B1 (en) Method for controlling operation of a wind turbine
EP3139038B1 (en) A method for estimating the surface condition of a rotating blade
CA2755154C (en) Method and system for adjusting a power parameter of a wind turbine
US9371819B2 (en) Method and system to actively pitch to reduce extreme loads on wind turbine
KR101476986B1 (en) Control device for a wind power plant
EP3511567B1 (en) State monitoring device, system and method for wind power generating device
JP7009237B2 (en) Wind power generation equipment and wind power generation system
US20130259682A1 (en) Method of rotor-stall prevention in wind turbines
JP2020502424A (en) Wind turbine power station level load management control strategy
US10724499B2 (en) Controlling wind turbines according to reliability estimates
TWI708893B (en) Wind power system
WO2012149984A1 (en) System and method for operating a wind turbine using an adaptive speed reference
US20180320660A1 (en) Adjustment Factor for Aerodynamic Performance Map
EP3608538A1 (en) Model-based repowering solutions for wind turbines
EP3259472B1 (en) Control for a wind turbine
JP6811150B2 (en) Wind farm control system and wind farm control method
US11728654B2 (en) Systems and methods for operating power generating assets
Ossmann et al. Multivariable controller design verification for a Liberty wind turbine
Barlas et al. Atmospheric full scale testing of a morphing trailing edge flap system for wind turbine blades
WO2022270025A1 (en) Method and device for diagnosing remaining service life of wind power generator
EP4339452A1 (en) Method for operating a wind turbine, control system and wind turbine

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200915

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210601

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210721

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: 20211221

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220112

R150 Certificate of patent or registration of utility model

Ref document number: 7009237

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150