WO2019150805A1 - Wind power generation device and wind power generation system - Google Patents
Wind power generation device and wind power generation system Download PDFInfo
- Publication number
- WO2019150805A1 WO2019150805A1 PCT/JP2018/046459 JP2018046459W WO2019150805A1 WO 2019150805 A1 WO2019150805 A1 WO 2019150805A1 JP 2018046459 W JP2018046459 W JP 2018046459W WO 2019150805 A1 WO2019150805 A1 WO 2019150805A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- wind
- load
- shear
- power generation
- wind power
- 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.)
- Ceased
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
-
- 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
Definitions
- the present invention relates to a wind power generation apparatus and a wind power generation system having a function of estimating a wind condition around the wind power generation apparatus.
- a megawatt-class wind power generator is equipped with a rotor that is mounted radially on a hub that rotates blades, a nacelle that supports the rotor via a main shaft, and a tower that allows the nacelle to support yaw rotation from below. Is frequently used. Wind power generators generate electricity using wind that changes every moment as an energy source. Therefore, when the wind speed or turbulence of the wind that actually flows into the wind turbine generator is severer than the design conditions, the load on the wind turbine generator may increase and damage to the components may be accelerated.
- Patent Literature 1 a table that associates wind speed, power generation amount or pitch angle, which has been created in advance, from wind speed data that can be obtained with a simple device configuration, and power generation amount or blade pitch angle data, and wind condition parameters.
- a device that estimates the life of components of a wind turbine generator and outputs maintenance information using a table that associates wind parameters and fatigue loads applied to the wind turbine generator has been proposed.
- a microwave or radar wave emitting device generally called a Doppler lidar is attached to a nacelle or hub of a wind power generator, and the wind speed distribution in front and rear of the wind power generator is measured, thereby improving the power generation efficiency.
- Methods have been proposed to control the pitch angle of the blades so as to maximize or minimize the load on the wind turbine generator.
- the present invention provides a wind power generation apparatus and a wind power generation system that can accurately estimate wind shear, which is a wind speed distribution, with a simple configuration.
- a wind turbine generator includes at least a rotor, a nacelle, and a wind turbine generator that has a tower that rotatably supports the nacelle, and a control device that controls the wind turbine generator,
- the control device includes a load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and the wind shear, and a window based on the load and the wind shear function.
- a wind condition estimation device having a wind speed distribution calculation unit for calculating shear is provided.
- the wind power generation system includes at least one wind power generation device, a control device that controls the wind power generation device, an electronic terminal that includes a display device, and a communication network that connects these devices so that they can communicate with each other.
- the control device includes a load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and wind shear, and the load and wind shear function And a wind speed estimation device having a wind speed distribution calculation unit for calculating wind shear based on the above.
- another wind power generation system includes at least one wind power generation device, a control device that controls the wind power generation device, an electronic terminal that includes a display device, and communication that connects these to each other so that they can communicate with each other.
- a network including a load measuring unit that measures a load applied to the wind turbine generator, and the electronic terminal storing a wind shear function that defines a relationship between the load and the wind shear
- a wind speed estimation device that calculates a wind shear based on the load input from the load measurement unit via the communication network and a wind shear function stored in the storage unit It is characterized by that.
- FIG. 1 is an overall schematic configuration diagram of a wind power generation system according to an embodiment of the present invention. It is a figure which shows the structure of the wind power generator of Example 1 which concerns on one Example of this invention. It is a figure for demonstrating the wind shear in Example 1.
- FIG. It is a functional block diagram of the control apparatus which comprises the wind power generator of Example 1.
- FIG. It is a figure for demonstrating the principle of the wind condition estimation apparatus which comprises the control apparatus shown in FIG. It is a figure for demonstrating the relationship between a power index and the moment in a hub center. It is a figure for demonstrating the influence of the wind speed in the relationship between the hail index and the moment in a hub center.
- FIG. 1 is an overall schematic configuration diagram of a wind power generation system according to an embodiment of the present invention.
- the wind power generation system 1 includes a wind power generation device 2 and an electronic terminal 4 installed in the operation management center 3 or a server (not shown), which includes a communication network 5 that can communicate with each other. Connected through.
- the communication network 5 may be wired or wireless.
- the wind power generator 2 includes a blade 24 that rotates by receiving wind, a hub 23 that supports the blade 24, a nacelle 22, and a tower 21 that rotatably supports the nacelle 22.
- a main shaft 26 that is connected to the hub 23 and rotates together with the hub 23, a speed increasing device 27 that is connected to the main shaft 26 and increases the rotational speed, and a rotor at a rotational speed increased by the speed increasing device 27.
- a generator 28 that rotates to generate electricity is provided.
- a portion that transmits the rotational energy of the blade 24 to the generator 28 is referred to as a power transmission unit.
- 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.
- a rotor 25 is constituted by the blade 24 and the hub 23.
- a power converter 30 for converting the frequency of power, a switching switch and a transformer for switching current (not shown), and control A device 31 and the like are arranged at the bottom (lower part) of the tower 21, a power converter 30 for converting the frequency of power, a switching switch and a transformer for switching current (not shown), and control A device 31 and the like are arranged.
- the control device 31 for example, a control panel or SCADA (Supervision Control And Data Acquisition) is used.
- FIG. 1 shows an example in which the rotor 25 is constituted by three blades 24 and a hub 23.
- the present invention is not limited to this, and the rotor 25 is connected to the hub 23 and at least one blade 24. May be configured. Embodiments of the present invention will be described below with reference to the drawings.
- FIG. 2 is a diagram showing the configuration of the wind turbine generator of Example 1 according to an example of the present invention.
- the structure of the wind speed distribution 11 of the height direction around the wind power generator 2 which concerns on a present Example is shown.
- FIG. 2 shows a state in which the wind turbine generator 2 is viewed from the side, and it is assumed that the wind is blowing from the left to the right of the page.
- the wind turbine generator 2 includes a rotor 25 radially attached to a hub 23 that rotates blades 24, a nacelle 22 that allows the rotor 25 to rotate and supports the hub 23 from the side, and a nacelle And a tower 21 which rotatably supports 22 with respect to a vertical axis from the lower part.
- the wind power generator 2 is a downwind type wind power generator in which the rotor 25 is located on the leeward side of the tower 21, but the upwind type wind power generator in which the rotor 25 is located on the windward side of the tower 21.
- a device may be used.
- the wind power generator 2 includes a strain sensor 7 attached to the tower 21.
- the strain sensor 7 is not limited to the tower 21 and may be installed on the nacelle 22 or the hub 23. Instead of the strain sensor 7, another load sensor such as an acceleration sensor may be used.
- the wind power generator 2 also includes a pitch angle control mechanism 6 for controlling the pitch angle of the blade 24, an anemometer 8 installed on the top of the nacelle 22, a thermometer 9 installed in the nacelle 22, and a barometer 10. May be provided.
- the anemometer 8, the thermometer 9, and the barometer 10 may be installed at other positions of the wind power generator 2, or may be installed outside the wind power generator 2 as long as they are in the vicinity of the wind power generator 2. good.
- 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.
- the wind speed distribution 11 usually changes in the height direction (Z direction), and generally the wind speed tends to increase toward the sky due to the boundary layer of the atmosphere. If the change in the wind speed in the height direction (Z direction) is called wind shear, and the power index representing the strength is ⁇ WS , the wind speed distribution in the height direction (Z direction) is expressed by the following equation (1). Can be assumed.
- V (z) is a wind speed at a height z from the ground surface
- z ref is a height that defines a reference wind speed
- V (z ref ) is a reference wind speed.
- the change in the wind speed due to the height increases as the power index ⁇ WS increases.
- the wind index measured by V (z ref ) in equation (1) and the height from the ground surface of the anemometer 8 are used as z ref.
- ⁇ WS is known, a wind speed distribution in the height direction can be obtained. That is, when the wind speed distribution is assumed as shown in Equation (1), the wind speed distribution estimation problem results in the power index ⁇ WS estimation problem.
- the wind speed distribution of Expression (1) is assumed and the power index ⁇ WS representing the strength of wind shear is estimated.
- the assumption of the wind speed distribution is not limited to Expression (1), A logarithm, a polynomial, or the like may be used, and a wind speed distribution may be assumed using a plurality of parameters.
- FIG. 4 is a functional block diagram of the control device 31 constituting the wind power generator 2 of the present embodiment.
- the control device 31 includes a wind speed measurement unit 12, a load measurement unit 13, a temperature measurement unit 14, an atmospheric pressure measurement unit 15, a storage unit 16 that stores a wind shear function 33, an atmospheric density calculation unit 17, a wind A shear estimation unit 18, a wind speed distribution calculation unit 19, an input I / F 34, an output I / F 35, and a communication I / F 36 are provided and are connected to each other via an internal bus 37.
- the wind speed measurement unit 12, the load measurement unit 13, the temperature measurement unit 14, the atmospheric pressure measurement unit 15, the storage unit 16 that stores the wind shear function 33, the atmospheric density calculation unit 17, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 are
- the wind condition estimating device 32 is configured not only to estimate the power index ⁇ WS that defines the strength of wind shear but also to determine the actual wind speed distribution. Note that if only the power index ⁇ WS is estimated, the wind condition estimation device 32 is stored in the wind speed measurement unit 12, the load measurement unit 13, the temperature measurement unit 14, the atmospheric pressure measurement unit 15, and the wind shear function 33.
- the air density calculation unit 17 and the wind shear estimation unit 18 may be used.
- the estimation device 32 may be configured.
- the wind speed measuring unit 12, the load measuring unit 13, the temperature measuring unit 14, the atmospheric pressure measuring unit 15, the atmospheric density calculating unit 17, the wind shear estimating unit 18, and the wind speed distribution calculating unit 19 constituting the wind condition estimating device 32 are, for example, , A processor such as a CPU (Central Processing Unit) (not shown), a ROM that stores various programs, a RAM that temporarily stores calculation process data, a storage device such as an external storage device, and a processor such as a CPU Reads out and executes various programs stored in the ROM, and stores the calculation result as the execution result in the RAM or the external storage device.
- CPU Central Processing Unit
- the wind speed measuring unit 12 calculates time series data of wind speed and average wind speed averaged over 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.
- the wind speed measurement 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.
- the average wind speed may be calculated every 10 minutes, or the average value may be output continuously using a moving average.
- 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 measurement value of the strain of the tower 21 is measured at the measurement position. By multiplying the section modulus, it is converted into a bending moment. Further, the bending moment is averaged over a predetermined time and transferred to the wind shear estimation unit 18 via the internal bus 37.
- the strain sensor 7 may be installed at one location, but when the wind direction changes, wind power generation is performed along the direction of the nacelle 22.
- strain sensors 7 It is desirable to install two or more strain sensors 7 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 toppled.
- the position in the height direction to 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 of the top. Therefore, it is not necessary to limit the measurement position to the top.
- the temperature measurement unit 14 inputs the temperature measured by the thermometer 10 via the input I / F 34 and the internal bus 37, and calculates an average temperature for a predetermined time.
- the temperature measurement unit 14 transfers the calculated average temperature to the atmospheric density calculation unit 17 via the internal bus 37.
- the atmospheric pressure measurement unit 15 inputs the atmospheric pressure measured by the barometer 11 via the input I / F 34 and the internal bus 37, and calculates an averaged atmospheric pressure over a predetermined time.
- the atmospheric pressure measurement unit 15 transfers the calculated average 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 wind shear estimation unit 18 via the internal bus 37 every 10 minutes, for example. To do. Note that other methods may be used as long as the atmospheric density can be calculated.
- the wind shear estimating unit 18 receives the bending moment transferred from the load measuring unit 13, the average wind speed transferred from the wind speed measuring unit 12, and the average value of the atmospheric density transferred from the atmospheric density calculating unit 17 as inputs.
- the wind shear function 33 that defines the relationship between the value of and the power index ⁇ WS representing the strength of wind shear
- an estimated value of the power index ⁇ WS is output for each predetermined time used for averaging.
- the average value of the bending moment, the average wind speed, and the atmospheric density is used as the input to the wind shear estimation unit 18.
- the load measurement unit Only the bending moment transferred from 13 may be input, and in addition to the bending moment, one of the average wind speed or the average value of the atmospheric density may be used. Further, instead of the average wind speed, another physical quantity that changes in accordance with the wind speed may be used. For example, the amount of power generation, pitch angle, rotor rotational speed, etc. can be used.
- the wind shear function 33 stored in the storage unit 16 stores a function that can uniquely output the power exponent ⁇ WS with respect to the average value of the bending moment, the average wind speed, and the atmospheric density.
- the function storage method is not particularly limited, and may be a response surface, a neural network, a mathematical expression such as a polynomial, or a data table in which variables are separated by bins.
- the function creation method is not particularly limited, but, for example, a Doppler lidar is temporarily attached to the top of the nacelle 22 constituting the wind power generator 2 and the power exponent ⁇ WS is calculated from the measured wind speed distribution by the least square method or the like.
- the wind speed distribution calculation unit 19 uses the above formula (from the estimated value of the power 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. 1) is used to calculate the time series data of the wind speed distribution. It is assumed that the height z ref from the ground surface of the anemometer 8 is known.
- the output I / F 35 outputs 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 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 by using the time series data of the wind speed calculated by the wind speed distribution calculation unit 19 and transferred via the internal bus 37. Or it transmits to the server which is not illustrated.
- b is the length of the blade 24
- ⁇ is the atmospheric density
- C t is a local thrust coefficient at a distance r from the center of the hub 23 of the blade 24
- V is a position r of a distance r from the center of the hub 23 of the blade 24
- C is a chord length at a distance r from the center of the hub 23 of the blade 24.
- b and c are values inherent to the blade 24 and are constant regardless of the operating conditions.
- Ct is also a value inherent to the blade 24, but varies according to the pitch angle, wind speed, and rotational speed.
- the wind speed has the distribution of the above formula (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 coning angle and rotates in a vertical plane. At this time, in the range where the blade 24 does not stall because the angle of attack of the blades 24 with increasing wind speed increases, C t is also monotonically increasing. Therefore, as can be seen from the above formula (2), when the local wind speed V increases, the thrust force T increases as C t V 2 increases.
- z ref represents the height of the hub 23 in FIG. 3, and the wind speed is increased above the hub 23 as the power index ⁇ 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 power index ⁇ WS increases.
- the moment in the direction in which the wind turbine generator 2 is toppled along the direction of the nacelle 22 at the center of the hub 23 when the blade 24 is at an arbitrary position in the plane of rotation, an increase in the power index ⁇ WS It can be seen that the moment increases monotonously with this.
- the average value of moments at the center of the hub 23 during one rotation of the rotor 25 also increases monotonously as shown in FIG. 6 as the power index ⁇ WS increases. Therefore, when the air density and the wind speed at the center of the hub 23 are constant, the power index ⁇ WS can be estimated from the moment at the center of the hub 23.
- the figure for demonstrating the relationship between the power index shown in FIG. 6 and the moment in a hub center is produced
- the hub 23 is obtained by extrapolating the moment measured at two positions in the height direction of the tower 21 instead of the moment at the tower 21 measured by the strain sensor 7. It is also possible to calculate the moment at the center of. As for the wind speed and the atmospheric density, it is not necessary to use the value at the center of the hub 23, and it may be measured near the wind power generator 2.
- the pitch angle and the rotational speed are constant, since these are generally determined by control according to the wind speed, when estimating the power index ⁇ WS using the wind speed, the pitch angle Even if the rotational speed changes, the influence on the estimation accuracy is small.
- the tilt angle and the coning angle are values inherent to the wind power generator 2 and are constant regardless of the operating conditions, and thus do not affect the estimation accuracy.
- the storage unit 16 that stores the wind shear function 33, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 that constitute the wind condition estimation device 32 are included in the control device 31.
- the present invention is not limited to this.
- the storage unit 16 storing the wind shear function 33, the wind shear estimation unit 18, and the wind speed distribution calculation unit 19 are connected to the electronic terminal 4 installed in the operation management center 3 shown in FIG. It is good also as a structure to mount.
- FIG. 9 is a functional block diagram of the control device 31a constituting the wind turbine generator of the second embodiment according to another embodiment of the present invention.
- the control device includes a reliability evaluation apparatus 40 that evaluates the reliability of the wind turbine generator 2 using the wind shear estimated by the wind shear estimation section 18 and the wind speed distribution calculation section 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 denoted by the same reference numerals, and the description overlapping with that in the first embodiment is omitted below.
- FIG. 9 only the input I / F 34, the wind shear estimation unit 18, the output I / F 35, the communication I / F 36, and the internal bus 37 are the same as in the first embodiment. Show.
- the control device 31a in addition to the input I / F 34, the wind shear estimation unit 18, the output I / F 35, and the communication I / F 36, the control device 31a according to the present embodiment includes an operation 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, which are connected to each other via an internal bus 37.
- a reliability evaluation device 40 is configured by the operating condition acquisition unit 41, the load calculation unit 42, the storage unit 16a that stores the design information 43, the reliability evaluation unit 44, and the information output unit 45. Note that the reliability evaluation device 40 may not include the information output unit 45.
- the wind shear estimation unit 18, the operation condition acquisition unit 41, the load calculation unit 42, the reliability evaluation unit 44, and the information output unit 45 are, for example, a processor such as a CPU (Central Processing Unit) (not shown) and a ROM that stores various programs.
- a processor such as a CPU (Central Processing Unit) (not shown) and a ROM that stores various programs.
- a processor such as a CPU reads out and executes various programs stored in the ROM, and an operation result that is an execution result Are stored in RAM or an external storage device.
- the reliability evaluation device 40 will be described by taking as an example the reliability evaluation of the blades 24 of the wind turbine generator 2. Note that the application destination of the reliability evaluation device 40 is not necessarily limited to the blade 24, and other configurations of the wind power generator 2 such as the nacelle 22 and the tower 21 as long as the wind shear affects the reliability. Parts may be used.
- the operating condition acquisition unit 41 acquires time history data of the operating conditions of the wind turbine generator 2 related to the load acting on the blade 24.
- the operating conditions include, for example, the wind speed and wind direction at the anemometer 8, the pitch angle and azimuth angle of the blade 24, the rotational speed of the rotor 25, the azimuth angle of the nacelle 22, and the amount of power generated by the wind power generator 2.
- a load sensor such as a strain sensor or an acceleration sensor
- the time history data thereof may be included.
- SCADA is used as the control device 31a
- the operating conditions may be acquired from this SCADA.
- the load calculating unit 42 uses the operating conditions transferred from the operating condition acquiring unit 41, the hull index ⁇ WS transferred from the wind shear estimating unit 18, and the design information 43 of the blade 24, and the load acting on the blade 24. Calculate the time history of.
- the wind speed distribution (wind shear) flowing into the rotor 25 is calculated from the wind speed at the anemometer 8 and the power index ⁇ WS using the above formula (1), and the wind speed distribution and operating conditions are calculated. Is used as an input for aeroelastic simulation based on blade element momentum theory or multibody dynamics.
- the design information 43 stored in the storage unit 16a includes, for example, the dimensions and mass distribution, rigidity distribution, aerodynamic coefficient of the blade 24 and the tower 21, the tilt angle and coning angle of the rotor 25, the dimension and mass distribution of the nacelle 22, and the aerodynamic coefficient. And design data such as a control program for the wind turbine generator 2.
- the design information 25 stored in the storage unit 16a includes data relating to the reliability of the component parts of the blade 24 as reliability information, for example, the dimensions, elastic modulus, section coefficient, stress concentration factor, S- N diagram is included.
- the reliability evaluation unit 44 performs reliability evaluation of the components of the blade 24 using the time history of the load acting on the blade 24 output from the load calculation unit 42 and the reliability information.
- the fatigue damage degree, remaining life, failure probability, etc. are calculated for the component parts.
- the degree of fatigue damage from load time history data 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 rain flow method is applied to the stress time history data to convert it into a stress amplitude appearance frequency distribution, and the S—of the component material stored in the obtained stress amplitude appearance frequency distribution and reliability information is stored. From the N diagram, the fatigue damage degree in the time history is calculated using the linear cumulative damage degree rule.
- 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.
- data output from the operating condition acquisition unit 41, the wind shear estimation unit 18, and the load calculation unit 42 may be displayed.
- the wind speed distribution flowing into the rotor 25 according to the above equation (1) using the wind speed at the anemometer 8 output from the operating condition acquisition unit 41 and the power index ⁇ WS output from the wind shear estimation unit 18. May be displayed.
- the reliability evaluation device 40 is provided in the control device 31a.
- the present invention is not limited to this, and the electronic terminal 4 installed in the operation management center 3 shown in FIG. 1 or not shown. It is good also as a structure mounted in a server.
- the reliability evaluation device 40 in addition to the effects of the first embodiment, by using the reliability evaluation device 40, the wind shear time history data estimated from the operating conditions of the wind turbine generator and the load of the wind turbine generator. Based on the above, it becomes possible to evaluate the reliability of the wind turbine generator.
- FIG. 10 is a functional block diagram of the control device 31b constituting the wind turbine generator of the third embodiment according to another embodiment of the present invention.
- the present embodiment is different from the first embodiment in that the wind power generator 2 is controlled 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 denoted by the same reference numerals, and the description overlapping with that in the first embodiment is omitted below.
- FIG. 10 only the input I / F 34, the wind shear estimation unit 18, the output I / F 35, the communication I / F 36, and the internal bus 37 are the same as in the first embodiment. Show.
- the control device 31b in addition to the input I / F 34, the wind shear estimation unit 18, the output I / F 35, and the communication I / F 36, the control device 31b according to the present embodiment includes an operation condition acquisition unit 41, a control amount calculation unit. 51 and a storage unit 16b for storing control information 52, which are connected to each other via an internal bus 37.
- the wind shear estimation unit 18, the operation condition acquisition unit 41, and the control amount calculation unit 51 temporarily store, for example, a processor (not shown) such as a CPU (Central Processing Unit), a ROM that stores various programs, and calculation process data.
- a processor such as a CPU reads out and executes various programs stored in the ROM, and stores an operation result as an execution result in the RAM or the external storage device.
- the operating condition acquisition unit 41 acquires the operating conditions of the wind turbine generator 2 as in the second embodiment.
- the control amount calculation unit 51 includes an operation condition transferred from the operation condition acquisition unit 41 via the internal bus 37, a power index ⁇ WS transferred from the wind shear estimation unit 18 via the internal bus 37, and control information 52.
- the pitch angle of the blade 24 and the rotational speed of the rotor 25 can be determined by, for example, maximizing the power generation amount or minimizing the load fluctuation of the blade 24. Note that the pitch angle may change while the rotor 25 rotates once.
- the wind speed distribution flowing into the rotor 4 is calculated from the wind speed at the anemometer 8 and the power exponent ⁇ WS using the above-described equation (1), and the flow into the blade 24 is obtained for any rotor rotational speed.
- the storage unit 16 b stores constants used for control and aerodynamic characteristics of the blade 24 as the control information 52.
- the above-described 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.
- the wind turbine generator is controlled based on the operating conditions of the wind turbine generator and the wind shear time history data estimated from the load of the wind turbine generator. It is possible to maximize the power generation amount and minimize the blade load fluctuation.
- FIG. 11 is a diagram showing a configuration of a wind turbine generator of Example 4 according to another example of the present invention.
- the present embodiment is different from the first embodiment in that horizontal wind shear is estimated using the wind condition estimating device 32 of the wind turbine generator 2 of the first embodiment described above.
- the same components as those in the first embodiment are denoted by the same reference numerals, and the description overlapping with that in the first embodiment is omitted below.
- FIG. 11 shows the configuration of the wind power generator 2 and the horizontal wind speed distribution 61 around it.
- the state which looked at the wind power generator 2 from upper direction is shown, and the wind shall blow from the paper surface left to the right.
- the wind power generator 2 has the same configuration as that of the first embodiment.
- the horizontal wind speed distribution 61 shown in FIG. 11 changes in the horizontal direction. If the coefficient representing the strength of wind shear in the horizontal direction is ⁇ WS , the horizontal wind speed distribution 61 is expressed by the following equation (3). Can be assumed as follows.
- V (y) is a wind speed at a position y in a direction perpendicular to the wind speed vector on the horizontal plane
- y ref is a position defining a reference wind speed
- V (y ref ) is a reference wind speed.
- the wind speed distribution estimation problem results in the coefficient ⁇ WS estimation problem.
- the coefficient ⁇ WS representing the strength of the wind shear in the horizontal direction is estimated on the assumption of the wind speed distribution of Expression (3), but the assumption of the wind speed distribution is limited to Expression (3). There is no need to do so, and a wind speed distribution may be assumed using a plurality of parameters.
- the configuration of the wind condition estimating device 32 of the wind turbine generator 2 in this embodiment for estimating the coefficient ⁇ WS that defines the strength of wind shear in the horizontal direction is the same as that in the first embodiment.
- the load measuring unit 13 measures the shear strain of the tower 21 using the strain sensor 7 and twists using the sectional secondary pole moment and the radius of the tower 21 at the measurement position. Convert to moment. Further, the twisting moment is averaged over a predetermined time and output. Since the torsional moment acts uniformly in the cross section of the tower 21, the strain sensor 7 may be installed at one place, and the position in the height direction need not be particularly limited.
- a wind shear estimation unit 18 receives the average values of torsional moment, wind speed, and atmospheric density as input, and defines a relationship between these values and a coefficient ⁇ WS representing the strength of wind shear in the horizontal direction.
- an estimated value of the coefficient ⁇ WS is output at every predetermined time used for averaging.
- the torsional moment, wind speed, and atmospheric density are used as the input to the wind shear estimation unit 18.
- the estimated accuracy of the coefficient ⁇ WS can be lowered, only the torsional moment may be input.
- either the wind speed or the atmospheric density may be used.
- the amount of power generation, pitch angle, rotor rotational speed, etc. can be used.
- the torsional moment is also generated due to a deviation between the nacelle 22 called the yaw error and the wind direction. Therefore, by adding the yaw error to the input of the wind shear estimation unit 18, high accuracy can be achieved.
- the wind shear function 33 is defined as a function that can uniquely output the coefficient ⁇ WS with respect to the input values of torsional moment, wind speed, and atmospheric density.
- the function storage method and the creation method are not particularly limited, and the same method as in the first embodiment may be used.
- the wind speed distribution calculating unit 19 uses the estimated value of the coefficient ⁇ WS obtained by the wind shear estimating unit 18 and the time series data of the wind speed obtained by the wind speed measuring unit 12 to calculate the time of the wind speed distribution using the above equation (3). Series data is calculated.
- the principle of horizontal 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, the moment acting on the center of the hub 23 by the horizontal wind shear is generated in the direction in which the nacelle 22 is rotated in the horizontal plane. Therefore, since the torsional moment acts on the tower 21, it is necessary to measure the torsional moment in order to estimate the coefficient ⁇ WS . In addition, by simultaneously measuring the bending moment and the torsional moment of the tower 21, the wind shear in the height direction (Z direction) and the horizontal direction (y direction) in the first embodiment and the fourth embodiment may be estimated simultaneously. The estimated horizontal wind shear may be input to the reliability evaluation device 40 of the second embodiment and the control device 31b of the third embodiment.
- FIG. 12 is a diagram showing a configuration of a wind turbine generator of Example 5 according to another example 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 in the blade 24.
- the same components as those in the first embodiment are denoted by the same reference numerals, and the description overlapping with that in the first embodiment is omitted below.
- FIG. 12 shows the configuration of the wind power generator 2 according to the present embodiment and the wind speed distribution 11 in the height direction around the wind power generator 2.
- the wind power generator 2 includes the strain sensor 7 on the blade 24 in place of the tower 21 in addition to the configuration similar to that of the first embodiment.
- the strain sensor 7 only needs to 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.
- the configuration of the wind condition estimating device 32 of the wind turbine generator 2 according to the present embodiment for estimating the power index ⁇ WS that defines the strength of wind shear is the same as that of the first embodiment.
- the load measuring unit 13 measures the strain of the blade 24 using the strain sensor 7 and multiplies the section coefficient at the measurement position to convert it into a bending moment. Further, the bending moment is statistically processed and a statistical value is output.
- the statistical value for example, a standard deviation of a bending moment at a predetermined time, a difference between a maximum value and a minimum value (maximum amplitude), or the like is used.
- the average value for each bin may be used by dividing the azimuth angle by bins, or the strain sensors 7 are attached to all the blades 24 constituting the rotor 25, and the center of the hub 23 is determined from the bending moment of each blade 24. A value obtained by calculation of the moment may be used.
- the strain sensor 7 may be installed at one place, it is desirable to install two or more strain sensors 7 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.
- the strain sensor 7 is preferably attached near the blade root of the blade 24, but the measurement position need not be limited to the vicinity of the blade root.
- the wind shear estimation unit 18 receives as input the statistical value of the bending moment of the blade 24, the wind speed, and the average value of the atmospheric density, and defines a relationship between these values and the power index ⁇ WS representing the strength of the wind shear. By using 33, an estimated value of the power exponent ⁇ WS is output every predetermined time.
- the bending moment, the wind speed, and the atmospheric density are used as the input to the wind shear estimation unit 18.
- the estimated accuracy of the power index ⁇ WS can be lowered, only the bending moment may be input.
- one of wind speed or atmospheric density may be used.
- the wind speed other physical quantities that change in accordance with the wind speed may be used. For example, the amount of power generation, pitch angle, rotor rotational speed, etc. can be used.
- the principle of wind shear estimation by the wind condition estimating device 32 is the same as that of the first embodiment except for the following points. That is, the thrust force acting on the blade 24 by wind shear changes according to the azimuth angle, as shown in FIG. 5, and therefore the bending moment measured by the strain sensor 7 also changes according to the azimuth angle. Therefore, the power index ⁇ 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 using the standard deviation and the maximum amplitude. Further, the power exponent ⁇ WS can be estimated 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.
- the bending moment and the thrust force when the blade 24 is positioned directly above the hub 23 monotonically increase with the increase of the power index ⁇ WS , and therefore, from these values, the power index ⁇ WS Can be estimated.
- the strain sensors 7 are 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 operated on the same principle as in the first embodiment. Shea can be estimated.
- the horizontal wind shear can be estimated, or the height direction and horizontal wind shear may be estimated simultaneously. Further, the estimated wind shear may be input to the reliability evaluation device 40 of the second embodiment and the control device 31b of the third embodiment.
- this invention is not limited to the above-mentioned Example, Various modifications are included.
- the above-described embodiments have been described in detail for easy understanding of the present invention, and are not necessarily limited to those having all the configurations described.
- a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment.
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)
Abstract
Description
本発明は、風力発電装置の周りの風況を推定する機能を有する風力発電装置及び風力発電システムに関する。 The present invention relates to a wind power generation apparatus and a wind power generation system having a function of estimating a wind condition around the wind power generation apparatus.
再生可能エネルギー活用への関心の高まりから、風力発電装置の世界的な市場拡大が予測されている。メガワット級の風力発電装置としては、ブレードを回転するハブに放射状に取りつけたロータと、主軸を介してロータを支持するナセルと、ナセルを下部からヨー回転を許容し支持するタワーを備えているものが頻繁に用いられる。
風力発電装置では、時々刻々と変化する風をエネルギー源として発電を行う。従って、実際に風力発電装置に流入する風の風速や乱れが設計条件よりも厳しい場合、風力発電装置の負荷が増大し、構成部品の損傷が加速する可能性がある。損傷の加速により、予期せぬ故障が発生した場合、故障部品の交換に要する時間に加えて、交換用部品の手配や工事用機材、作業員の手配に時間を要するため、計画された部品交換を実施する場合よりも、風力発電装置の稼動停止時間が増大し、発電量が減少する懸念がある。
この対策として、エネルギー源となる流入風を推定・計測することで、風力発電装置の構成部品の損傷を推定する方法や、負荷を軽減するように制御する方法がある。例えば、特許文献1では、簡便な装置構成で入手可能な風速と、発電量若しくはブレードのピッチ角のデータから、予め作成しておいた風速、発電量若しくはピッチ角と、風況パラメータを関連付けるテーブルと、風況パラメータと風力発電装置に加わる疲労荷重を関連付けるテーブルを用いて、風力発電装置の構成部品の寿命を推定し、メンテナンス情報を出力する装置が提案されている。
一方、特許文献2では、一般にドップラーライダと称されるマイクロ波またはレーダ波発射装置を風力発電装置のナセルまたはハブに取り付け、風力発電装置前方および後方の風速分布を計測することで、発電効率を最大化または風力発電装置の負荷を最小化するようにブレードのピッチ角を制御する方法が提案されている。
Growing interest in the use of renewable energy is expected to expand the global market for wind turbines. A megawatt-class wind power generator is equipped with a rotor that is mounted radially on a hub that rotates blades, a nacelle that supports the rotor via a main shaft, and a tower that allows the nacelle to support yaw rotation from below. Is frequently used.
Wind power generators generate electricity using wind that changes every moment as an energy source. Therefore, when the wind speed or turbulence of the wind that actually flows into the wind turbine generator is severer than the design conditions, the load on the wind turbine generator may increase and damage to the components may be accelerated. If an unexpected failure occurs due to acceleration of damage, in addition to the time required to replace the failed part, it takes time to arrange replacement parts, construction equipment, and workers. There is a concern that the operation stop time of the wind power generator increases and the amount of power generation decreases compared to the case of implementing the above.
As countermeasures, there are a method for estimating damage to the components of the wind power generator and a method for controlling the load to be reduced by estimating and measuring the inflow wind as an energy source. For example, in
On the other hand, in
しかしながら、特許文献1に開示される構成では、発電量やピッチ角は概ね風速に対して一意の値になるよう制御されているため、風速、発電量、ピッチ角から風速分布であるウィンドシアを高精度に推定することは困難であり、構成部品の寿命推定精度への影響も懸念される。また、特許文献2に開示される構成では、ドップラーライダにより風速分布を正確に計測できるものの、このような計測機器は高価であるため、風力発電所内(ウィンドファーム内)の全ての風力発電装置に設置することはコスト面から現実的ではない。
However, in the configuration disclosed in
そこで、本発明は、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供する。 Therefore, the present invention provides a wind power generation apparatus and a wind power generation system that can accurately estimate wind shear, which is a wind speed distribution, with a simple configuration.
上記課題を解決するため、本発明に係る風力発電装置は、少なくともロータ及びナセル並びにナセルを回動可能に支持するタワーを有する風力発電装置と、風力発電装置を制御する制御装置とを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記荷重及びウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
また、本発明に係る風力発電システムは、少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記荷重及びウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
また、本発明に係る他の風力発電システムは、少なくとも1基の風力発電装置と、風力発電装置を制御する制御装置と、表示装置を有する電子端末と、これらを相互に通信可能に接続する通信ネットワークを備え、前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部を有し、前記電子端末は、荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重及び前記前記記憶部に格納されるウィンドシア関数に基づきウィンドシアを計算する風速分布計算部とを有する風況推定装置を備えることを特徴とする。
In order to solve the above problems, a wind turbine generator according to the present invention includes at least a rotor, a nacelle, and a wind turbine generator that has a tower that rotatably supports the nacelle, and a control device that controls the wind turbine generator, The control device includes a load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and the wind shear, and a window based on the load and the wind shear function. A wind condition estimation device having a wind speed distribution calculation unit for calculating shear is provided.
In addition, the wind power generation system according to the present invention includes at least one wind power generation device, a control device that controls the wind power generation device, an electronic terminal that includes a display device, and a communication network that connects these devices so that they can communicate with each other. The control device includes a load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and wind shear, and the load and wind shear function And a wind speed estimation device having a wind speed distribution calculation unit for calculating wind shear based on the above.
In addition, another wind power generation system according to the present invention includes at least one wind power generation device, a control device that controls the wind power generation device, an electronic terminal that includes a display device, and communication that connects these to each other so that they can communicate with each other. A network including a load measuring unit that measures a load applied to the wind turbine generator, and the electronic terminal storing a wind shear function that defines a relationship between the load and the wind shear And a wind speed estimation device that calculates a wind shear based on the load input from the load measurement unit via the communication network and a wind shear function stored in the storage unit It is characterized by that.
本発明によれば、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供することが可能となる。
上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
ADVANTAGE OF THE INVENTION According to this invention, it becomes possible to provide the wind power generator and wind power generation system which can estimate the wind shear which is a wind speed distribution with high precision by simple structure.
Problems, configurations, and effects other than those described above will be clarified by the following description of embodiments.
図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
また、風力発電装置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)が用いられる。
The
なお、図1に示す風力発電装置2は、3枚のブレード24とハブ23にてロータ25を構成する例を示すが、これに限られず、ロータ25はハブ23と少なくとも1枚のブレード24にて構成しても良い。
以下、図面を用いて本発明の実施例について説明する。
1 shows an example in which the
Embodiments of the present invention will be described below 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 the configuration of the wind turbine generator of Example 1 according to an example of the present invention. In FIG. 2, the structure of the
風力発電装置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
風速分布11は、通常、高さ方向(Z方向)に変化しており、一般的には大気の境界層によって上空ほど風速が大きくなる傾向がある。この高さ方向(Z方向)の風速の変化をウィンドシアと称し、その強弱を表す冪指数をαWSとすると、高さ方向(Z方向)の風速分布を以下の式(1)のように仮定することができる。
The
ここで、V(z)は地表から高さzにおける風速であり、zrefは基準となる風速を定義する高さ、V(zref)は基準となる風速を表している。図3に示すように、冪指数αWSが大きくなるほど、高さによる風速の変化は大きくなる。例えば、基準となる風速を風速計8で計測する場合、式(1)のV(zref)に計測した風速、zrefに風速計8の地表からの高さを用いることで、ある冪指数αWSがわかれば高さ方向の風速分布を得ることができる。すなわち、風速分布を式(1)のように仮定する場合、風速分布の推定問題は、冪指数αWSの推定問題に帰着する。なお、本実施例では、式(1)の風速分布を仮定して、ウィンドシアの強弱を表す冪指数αWSを推定することを考えるが、風速分布の仮定は式(1)に限らず、対数や多項式などを用いても良く、複数のパラメータを用いて風速分布を仮定しても良い。
Here, V (z) is a wind speed at a height z from the ground surface, z ref is a height that defines a reference wind speed, and V (z ref ) is a reference wind speed. As shown in FIG. 3, the change in the wind speed due to the height increases as the power index α WS increases. For example, when the
図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
The wind
風速計測部12は、風速計8にて計測され入力I/F34及び内部バス37介して入力される風速から、風速の時系列データ及び所定時間で平均化処理を施した平均風速を算出する。そして、風速計測部12は、内部バス37を介して平均風速をウィンドシア推定部18へ転送すると共に、風速の時系列データを風速分布計算部19へ転送する。なお、平均化処理では、例えば10分ごとに平均風速を算出しても良く、移動平均を用いて平均値を連続的に出力しても良い。
The wind
荷重計測部13は、ひずみセンサ7にて計測されるタワー21のひずみを、入力I/F34及び内部バス37介して入力し、入力されたタワー21のひずみの計測値に対し、計測位置での断面係数を乗ずることで、曲げモーメントに変換する。更に、曲げモーメントを所定時間での平均化処理を施して、ウィンドシア推定部18へ内部バス37を介して転送する。風力発電装置2の建設位置での風向が概ね一方向である場合には、ひずみセンサ7の設置は一箇所でもかまわないが、風向が変化する場合には、ナセル22の方向に沿って風力発電装置2を転倒させる向きの曲げモーメント成分を算出するために、タワー21の同一高さ地点に二箇所以上ひずみセンサ7を設置することが望ましい。また、ひずみセンサ7を取り付ける高さ方向の位置は、タワー21の頂部付近であることが望ましいが、直接頂部のモーメントを計測しなくとも、風速とロータ25のスラスト係数から補正を施すこともできるため、計測位置を頂部に限定する必要はない。
The
温度計測部14は、温度計10にて計測された温度を、入力I/F34及び内部バス37介して入力し、所定時間での平均化温度を算出する。温度計測部14は、算出した平均化温度を大気密度計算部17へ内部バス37を介して転送する。気圧計測部15は、気圧計11にて計測された気圧を、入力I/F34及び内部バス37介して入力し、所定時間での平均化気圧を算出する。気圧計測部15は、算出した平均化気圧を大気密度計算部17へ内部バス37を介して転送する。
大気密度計算部17は、温度と気圧から気体の状態方程式を用いて、大気密度を計算し、例えば、10分ごとに大気密度の平均値をウィンドシア推定部18へ内部バス37を介して転送する。なお、大気密度を算出できれば、他の方法を用いてもかまわない。
The
The atmospheric
ウィンドシア推定部18は、荷重計測部13より転送される曲げモーメント、風速計測部12より転送される平均風速、及び、大気密度計算部17より転送される大気密度の平均値を入力として、これらの値とウィンドシアの強弱を表す冪指数αWSとの関係を定義するウィンドシア関数33を用いることにより、平均化に用いた所定時間ごとに冪指数αWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、曲げモーメント、平均風速、大気密度の平均値を用いているが、冪指数αWSの推定精度低下を許容できる場合には、荷重計測部13より転送される曲げモーメントのみを入力としても良く、曲げモーメントに加えて平均風速または大気密度の平均値の一方を用いても良い。また、平均風速に代えて、風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。
The wind
記憶部16に格納されるウィンドシア関数33は、曲げモーメント、平均風速、大気密度の平均値に対して冪指数αWSを一意に出力できる関数を保存している。関数の保存方法は特に限定されず、応答曲面やニューラルネットワークのほか、多項式のような数式でも、各変数をビンで区切ったデータテーブルでも構わない。関数の作成方法についても特に限定されないが、例えば、風力発電装置2を構成するナセル22の上部に一時的にドップラーライダを取り付けて、計測した風速分布から冪指数αWSを最小二乗法などで計算し、曲げモーメント、平均風速、大気密度の平均値を設計変数、冪指数αWSを目的関数として応答曲面を作成する方法がある。他の方法としては、数値シミュレーションを用いて、冪指数αWS、平均風速、大気密度の平均値を変化させた場合の曲げモーメントを計算し、計算結果をニューラルネットワークに学習させて、冪指数αWSを推定することもできる。
The
風速分布計算部19は、ウィンドシア推定部18で得られた冪指数αWSの推定値、及び、風速計測部12より内部バス37を介して転送された風速の時系列データから、上記式(1)を用いて風速分布の時系列データを算出する。なお、風速計8の地表からの高さzrefは既知であるとする。
出力I/F35は、風速分布計算部19により算出され、内部バス37を介して転送された風速分布の時系列データを図示しない表示部に出力する。
通信I/F36は、風速分布計算部19により算出され、内部バス37を介して転送された風速分布の時系列データを、通信ネットワーク5を介して運転管理センター3内に設置される電子端末4或いは図示しないサーバへ送信する。
The wind speed
The output I /
The communication I /
次に、風況推定装置32によるウィンドシアの推定原理について説明する。ロータ25に流入する風速がウィンドシアにより、上記式(1)のような分布をとると仮定する。ブレード24に作用する風向方向の空気力(スラスト力)Tは式(2)で表すことができる。
Next, the principle of wind shear estimation by the wind
ここで、bはブレード24の長さ、ρは大気密度、Ctはブレード24のハブ23の中心から距離rの位置における局所スラスト係数、Vはブレード24のハブ23の中心から距離rの位置における局所風速、cはブレード24のハブ23の中心から距離rの位置における翼弦長である。これらのうち、b及びcはブレード24に固有の値であり、運転条件によらず一定である。また、Ctもブレード24に固有の値であるが、ピッチ角、風速、回転数に応じて変化する。ここで簡単のため、風速は高さ方向に上述の式(1)の分布を有するが時間変化せず、ピッチ角、回転数は一定であると仮定する。また、ロータ25はチルト角、コーニング角を持たず鉛直平面内を回転すると仮定する。このとき、ブレード24が失速しない範囲では風速の増加に伴ってブレード24の迎え角が増加するため、Ctも単調増加する。したがって、上記式(2)からわかるように、局所風速Vが増加する場合、CtV2が増加することによってスラスト力Tも増加する。
Here, b is the length of the
今、大気密度及びハブ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
実際には回転するハブ23でのモーメント計測は困難であるため、ひずみセンサ7によって計測したタワー21でのモーメントで代用する、タワー21の高さ方向二箇所で計測したモーメントを外挿してハブ23の中心でのモーメントを算出するなどしても良い。風速と大気密度については、ハブ23の中心の値を用いる必要はなく、風力発電装置2の近傍で計測すれば良い。また、ピッチ角、回転数は一定と仮定したが、これらは風速に応じた制御によって決定されることが一般的であるため、風速を用いて冪指数αWSを推定する場合には、ピッチ角、回転数が変化しても推定精度への影響は小さい。チルト角、コーニング角は、風力発電装置2に固有の値であり、運転条件によらず一定であるため、推定精度への影響はない。
In practice, since it is difficult to measure the moment at the rotating
なお、本実施例では図4に示したように、風況推定装置32を構成する、ウィンドシア関数33を格納する記憶部16、ウィンドシア推定部18、及び風速分布計算部19を制御装置31内に実装する場合を一例として説明したがこれに限られるものではない。例えば、これらウィンドシア関数33を格納する記憶部16、ウィンドシア推定部18、及び風速分布計算部19を、図1に示した運転管理センター3内に設置される電子端末4或いは図示しないサーバに実装する構成としても良い。
In the present embodiment, as shown in FIG. 4, the
以上の通り、本実施例によれば、簡易な構成にて、風速分布であるウィンドシアを高精度に推定し得る風力発電装置及び風力発電システムを提供することが可能となる。 As described above, according to the present embodiment, it is possible to provide a wind power generation apparatus and a wind power generation system that can accurately estimate wind shear, which is a wind speed distribution, 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
図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, in addition to the input I /
以下では、風力発電装置2のブレード24の信頼性評価を一例として、信頼性評価装置40について説明する。なお、信頼性評価装置40の適用先は、ブレード24に限定される必要はなく、ウィンドシアが信頼性に影響を与える部位であれば、ナセル22やタワー21など風力発電装置2の他の構成部品でも良い。
Hereinafter, the
運転条件取得部41は、ブレード24に作用する荷重に関連する風力発電装置2の運転条件の時刻歴データを取得する。運転条件とは、例えば、風速計8での風速や風向、ブレード24のピッチ角やアジマス角、ロータ25の回転速度、ナセル22の方位角、風力発電装置2の発電量などである。風力発電装置2にひずみセンサや加速度センサなどの荷重センサが取り付けられている場合には、それらの時刻歴データを含んでも良い。また、制御装置31aとして、例えばSCADAが用いられる場合には、このSCADAから運転条件を取得しても良い。
The operating
荷重計算部42は、運転条件取得部41から転送される運転条件、ウィンドシア推定部18から転送される冪指数αWSと、ブレード24の設計情報43とを用いて、ブレード24に作用する荷重の時刻歴を計算する。荷重の計算方法としては、例えば、風速計8での風速と冪指数αWSから上述の式(1)を用いてロータ25に流入する風速分布(ウィンドシア)を計算し、風速分布と運転条件の時刻歴データを翼素運動量理論或いはマルチボディダイナミクスなどに基づく空力弾性シミュレーションの入力とする方法がある。
記憶部16aに格納される設計情報43は、例えば、ブレード24とタワー21の寸法や質量分布、剛性分布、空力係数、ロータ25のチルト角やコーニング角、ナセル22の寸法や質量分布、空力係数、風力発電装置2の制御プログラムなどの設計データである。なお、記憶部16aに格納される設計情報25には、信頼性情報としてのブレード24の構成部品の信頼性に関するデータ、例えば、構成部品の寸法、弾性係数、断面係数、応力集中係数、S―N線図などが含まれる。
The
The
信頼性評価部44は、荷重計算部42から出力されるブレード24に作用する荷重の時刻歴と上記信頼性情報とを用いて、ブレード24の構成部品の信頼性評価を行う。信頼性評価としては、構成部品に対する疲労損傷度や余寿命、破壊確率などの計算を行う。例えば、荷重の時刻歴データから疲労損傷度を計算する場合には、以下のような方法がある。まず、ブレード24に作用する荷重から構成部品に作用する応力の時刻歴を計算する。次に、応力の時刻歴データにレインフロー法を適用して応力振幅の出現頻度分布に変換し、得られた応力振幅の出現頻度分布と信頼性情報に保存される構成部品の材料のS―N線図から、線形累積損傷度則を用いて、時刻歴における疲労損傷度を計算する。
The
情報出力部45は、信頼性評価部44から内部バス37を介して転送される評価結果を表やグラフ、コンター図として表示する。なお、運転条件取得部41、ウィンドシア推定部18、荷重計算部42から出力されるデータを表示してもよい。例えば、運転条件取得部41から出力される風速計8での風速と、ウィンドシア推定部18から出力される冪指数αWSを用いて上述の式(1)により、ロータ25に流入する風速分布を表示しても良い。
なお、本実施例では、信頼性評価装置40を、制御装置31a内に設ける構成としたが、これに限られず、図1に示した運転管理センター3内に設置される電子端末4或いは図示しないサーバに実装する構成としても良い。
The
In this embodiment, the
以上の通り本実施例によれば、実施例1の効果に加え、信頼性評価装置40を用いることで、風力発電装置の運転条件と、風力発電装置の荷重から推定したウィンドシアの時刻歴データに基づいて、風力発電装置の信頼性を評価することが可能となる。
As described above, according to the present embodiment, in addition to the effects of the first embodiment, by using the
図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 the
図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, in addition to the input I /
運転条件取得部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
The storage unit 16 b stores constants used for control and aerodynamic characteristics of the
以上の通り本実施例によれは、実施例1の効果に加え、風力発電装置の運転条件と、風力発電装置の荷重から推定したウィンドシアの時刻歴データに基づいて、風力発電装置を制御し、発電量の最大化やブレードの荷重変動の最小化が可能となる。 As described above, according to the present embodiment, in addition to the effects of the first embodiment, the wind turbine generator is controlled based on the operating conditions of the wind turbine generator and the wind shear time history data estimated from the load of the wind turbine generator. It is possible to maximize the power generation amount and minimize the blade load fluctuation.
図11は、本発明の他の実施例に係る実施例4の風力発電装置の構成を示す図である。本実施例では、上述の実施例1の風力発電装置2の風況推定装置32を用いて、水平方向のウィンドシアを推定する点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。
FIG. 11 is a diagram showing a configuration of a wind turbine generator of Example 4 according to another example of the present invention. The present embodiment is different from the first embodiment in that horizontal wind shear is estimated using the wind
図11では、風力発電装置2及びその周りでの水平方向の風速分布61の構成を示している。図11では、風力発電装置2を上方から眺めた状態を示しており、風は紙面左から右に吹いているものとする。なお、風力発電装置2は、上述の実施例1と同様の構成を備えている。
FIG. 11 shows the configuration of the
図11に示す水平方向の風速分布61は、水平方向に変化しており、この水平方向のウィンドシアの強弱を表す係数をβWSとすると、水平方向の風速分布61を以下の式(3)のように仮定することができる。
The horizontal
ここで、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 a wind speed at a position y in a direction perpendicular to the wind speed vector on the horizontal plane, y ref is a position defining a reference wind speed, and V (y ref ) is a reference wind speed. Yes. As can be seen from Equation (3), the larger the coefficient β WS , the greater the change in the wind speed in the horizontal direction. For example, when the
水平方向のウィンドシアの強弱を定義する係数βWSを推定するための、本実施例における風力発電装置2の風況推定装置32の構成は実施例1と同様である。以下では、風況推定装置32の詳細について、実施例1と異なる点について述べる。
本実施例において、荷重計測部13(図示せず)では、ひずみセンサ7を用いてタワー21のせん断ひずみを計測し、計測位置での断面二次極モーメントとタワー21の半径を用いて、捻りモーメントに変換する。更に、捻りモーメントを所定時間での平均化処理を施して出力する。捻りモーメントはタワー21の断面内で均一に作用するので、ひずみセンサ7の設置は一箇所でもよく、高さ方向の位置も特に限定する必要はない。
The configuration of the wind
In the present embodiment, the load measuring unit 13 (not shown) measures the shear strain of the
ウィンドシア推定部18(図示せず)では、捻りモーメント、風速、大気密度の平均値を入力として、これらの値と水平方向のウィンドシアの強弱を表す係数βWSとの関係を定義するウィンドシア関数33を用いることにより、平均化に用いた所定時間ごとに係数βWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、捻りモーメント、風速、大気密度を用いているが、係数βWSの推定精度低下を許容できる場合には、捻りモーメントのみを入力としても良く、捻りモーメントに加えて風速または大気密度の一方を用いても良い。また、風速の代わりに風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。なお、捻りモーメントは、ヨーエラーと呼ばれるナセル22と風向とのずれによっても発生するため、ヨーエラーをウィンドシア推定部18の入力に加えることで、高精度化が可能である。
A wind shear estimation unit 18 (not shown) receives the average values of torsional moment, wind speed, and atmospheric density as input, and defines a relationship between these values and a coefficient β WS representing the strength of wind shear in the horizontal direction. By using the
ウィンドシア関数33は、捻りモーメント、風速、大気密度の入力値に対して係数βWSを一意に出力できる関数として定義されている。関数の保存方法、作成方法は特に限定されず、成就の実施例1と同様の方法を用いても良い。
風速分布計算部19では、ウィンドシア推定部18で得られた係数βWSの推定値と風速計測部12で得られた風速の時系列データから、上記式(3)を用いて風速分布の時系列データを算出する。
The
The wind speed
風況推定装置32による水平方向のウィンドシアの推定原理は、以下の点を除いて実施例1と同様である。すなわち、水平方向のウィンドシアによってハブ23の中心に作用するモーメントは、ナセル22を水平面内で回転させる向きに発生する。したがって、タワー21には、捻りモーメントが作用するため、係数βWSの推定には、捻りモーメントを計測する必要がある。なお、タワー21の曲げモーメントと捻りモーメントを同時に計測することにより、実施例1と実施例4の高さ方向(Z方向)、水平方向(y方向)のウィンドシアを同時に推定しても良く、推定した水平方向のウィンドシアを、上述の実施例2の信頼性評価装置40、実施例3の制御装置31bの入力としても良い。
The principle of horizontal wind shear estimation by the wind
以上の通り本実施例によれば、風力発電装置の荷重を計測することにより、簡易な構成にて、水平方向のウィンドシアを高精度に推定することが可能となる。 As described above, according to the present embodiment, it is possible to estimate the wind shear in the horizontal direction with high accuracy by measuring the load of the wind turbine generator.
図12は、本発明の他の実施例に係る実施例5の風力発電装置の構成を示す図である。本実施例では、上述の実施例1に示したタワー21に設置されるひずみセンサ7をブレード24に設置した点が実施例1と異なる。実施例1と同様の構成要素に同一符号を付し、以下では、実施例1と重複する説明を省略する。
FIG. 12 is a diagram showing a configuration of a wind turbine generator of Example 5 according to another example of the present invention. This embodiment is different from the first embodiment in that the
図12に本実施例に係る風力発電装置2及びその周りでの高さ方向の風速分布11の構成を示す。図12では、風力発電装置2を側方から眺めた状態を示しており、風は紙面左から右に吹いているものとする。風力発電装置2は、実施例1と同様の構成に加えて、タワー21に代えて、ブレード24にひずみセンサ7を備えている。ひずみセンサ7は、ロータ25を構成する少なくとも一つのブレード24に取り付けられていれば良く、ひずみセンサ7に代えて加速度センサのような他の荷重センサを用いても良い。
FIG. 12 shows the configuration of the
ウィンドシアの強弱を定義する冪指数α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
In the present embodiment, the
ウィンドシア推定部18は、ブレード24の曲げモーメントの統計量、風速、大気密度の平均値を入力として、これらの値とウィンドシアの強弱を表す冪指数αWSとの関係を定義するウィンドシア関数33を用いることにより、所定時間ごとに冪指数αWSの推定値を出力する。本実施例では、ウィンドシア推定部18への入力として、曲げモーメント、風速、大気密度を用いているが、冪指数αWSの推定精度低下を許容できる場合には、曲げモーメントのみを入力としても良く、曲げモーメントに加えて風速または大気密度の一方を用いても良い。また、風速に代えて風速に対応して変化する他の物理量を用いても良い。例えば、発電量やピッチ角、ロータ回転数などを用いることができる。
The wind
風況推定装置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
なお、ブレード24が水平となるときの曲げモーメントを用いれば、水平方向のウィンドシアを推定することもでき、高さ方向、水平方向のウィンドシアを同時に推定しても良い。また、推定したウィンドシアを、上述の実施例2の信頼性評価装置40、実施例3の制御装置31bの入力としても良い。
In addition, if the bending moment when the
以上の通り本実施例によっても、風力発電装置の荷重を計測することにより、簡易な構成にて、ウィンドシアを高精度に推定することが可能となる。 As described above, according to this embodiment as well, it is possible to estimate the wind shear with high accuracy with a simple configuration by measuring the load of the wind turbine generator.
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。 In addition, this invention is not limited to the above-mentioned Example, Various modifications are included. For example, the above-described embodiments have been described in detail for easy understanding of the present invention, and are not necessarily limited to those having all the configurations described. Further, a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added 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…水平方向の風速分布
DESCRIPTION OF
35 ... Output I / F
36 ... Communication I / F
37 ...
Claims (15)
前記制御装置は、
風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記荷重及びウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電装置。 A wind turbine generator having at least a rotor and a nacelle and a tower that rotatably supports the nacelle, and a control device that controls the wind turbine generator,
The controller is
A load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and the wind shear, and a wind shear that is calculated based on the load and the wind shear function. A wind power generation apparatus comprising a wind condition estimation device having a wind shear estimation unit.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電装置。 The wind turbine generator according to claim 1,
The wind shear estimation unit calculates a wind shear that is a wind velocity distribution in a height direction based on a wind shear function stored in the storage unit and the load input from the load measurement unit. Power generation device.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電装置。 The wind turbine generator according to claim 1,
The wind shear estimation unit calculates a wind shear, which is a wind velocity distribution in a horizontal direction, based on a wind shear function stored in the storage unit and the load input from the load measurement unit. apparatus.
前記制御装置は、
前記高さ方向の風速分布であるウィンドシアを用いて、前記風力発電装置の信頼性を評価する信頼性評価装置を有することを特徴とする風力発電装置。 The wind turbine generator according to claim 2,
The controller is
A wind turbine generator having a reliability evaluation device that evaluates reliability of the wind turbine generator using wind shear that is a wind velocity distribution in the height direction.
前記荷重計測部は、前記タワーに付加される荷重を計測することを特徴とする風力発電装置。 In the wind turbine generator according to claim 2 or claim 3,
The said load measurement part measures the load added to the said tower, The wind power generator characterized by the above-mentioned.
前記荷重計測部は、前記風力発電装置に付加される荷重として、ひずみまたは加速度を計測することを特徴とする風力発電装置。 In the wind turbine generator according to claim 2 or claim 3,
The said load measurement part measures distortion or acceleration as a load added to the said wind power generator, The wind power generator characterized by the above-mentioned.
前記制御装置は、
風力発電装置に付加される荷重を計測する荷重計測部と、前記荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記荷重及びウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電システム。 At least one wind power generator, a control device that controls the wind power generator, an electronic terminal that includes a display device, and a communication network that connects these to each other so that they can communicate with each other,
The controller is
A load measuring unit that measures a load applied to the wind turbine generator, a storage unit that stores a wind shear function that defines a relationship between the load and the wind shear, and a wind shear that is calculated based on the load and the wind shear function. A wind power generation system comprising a wind condition estimation device having a wind shear estimation unit.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。 The wind power generation system according to claim 7,
The wind shear estimation unit calculates a wind shear that is a wind speed distribution in a height direction based on a wind shear function stored in the storage unit and the load input from the load measurement unit via the communication network. Wind power generation system characterized by that.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。 The wind power generation system according to claim 7,
The wind shear estimation unit calculates a wind shear that is a wind speed distribution in a horizontal direction based on a wind shear function stored in the storage unit and the load input from the load measurement unit via the communication network. Wind power generation system characterized by
前記制御装置は、
前記高さ方向の風速分布であるウィンドシアを用いて、前記風力発電装置の信頼性を評価する信頼性評価装置を有することを特徴とする風力発電システム。 The wind power generation system according to claim 8,
The controller is
A wind power generation system comprising a reliability evaluation device that evaluates the reliability of the wind power generation device using wind shear which is a wind velocity distribution in the height direction.
前記荷重計測部は、前記風力発電装置のタワーに付加される荷重を計測することを特徴とする風力発電システム。 In the wind power generation system according to claim 8 or claim 9,
The said load measurement part measures the load added to the tower of the said wind power generator, The wind power generation system characterized by the above-mentioned.
前記荷重計測部は、前記風力発電装置に付加される荷重として、ひずみまたは加速度を計測することを特徴とする風力発電システム。 In the wind power generation system according to claim 8 or claim 9,
The said load measurement part measures distortion or acceleration as a load added to the said wind power generator, The wind power generation system characterized by the above-mentioned.
前記制御装置は、風力発電装置に付加される荷重を計測する荷重計測部を有し、
前記電子端末は、荷重とウィンドシアとの関係を定義するウィンドシア関数を格納する記憶部と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重及び前記記憶部に格納されるウィンドシア関数に基づきウィンドシアを計算するウィンドシア推定部とを有する風況推定装置を備えることを特徴とする風力発電システム。 At least one wind power generator, a control device that controls the wind power generator, an electronic terminal that includes a display device, and a communication network that connects these to each other so that they can communicate with each other,
The control device has a load measuring unit that measures a load applied to the wind turbine generator,
The electronic terminal includes a storage unit that stores a wind shear function that defines a relationship between a load and wind shear, the load that is input from the load measurement unit via the communication network, and a window that is stored in the storage unit. A wind power generation system comprising a wind condition estimation device having a wind shear estimation unit for calculating wind shear based on a shear function.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき高さ方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。 The wind power generation system according to claim 13,
The wind shear estimation unit calculates a wind shear that is a wind speed distribution in a height direction based on a wind shear function stored in the storage unit and the load input from the load measurement unit via the communication network. Wind power generation system characterized by that.
前記ウィンドシア推定部は、前記記憶部に格納されるウィンドシア関数と、前記通信ネットワークを介して前記荷重計測部から入力される前記荷重に基づき水平方向の風速分布であるウィンドシアを計算することを特徴とする風力発電システム。 The wind power generation system according to claim 13,
The wind shear estimation unit calculates a wind shear that is a wind speed distribution in a horizontal direction based on a wind shear function stored in the storage unit and the load input from the load measurement unit via the communication network. Wind power generation system characterized by
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018-014369 | 2018-01-31 | ||
| JP2018014369A JP7009237B2 (en) | 2018-01-31 | 2018-01-31 | Wind power generation equipment and wind power generation system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2019150805A1 true WO2019150805A1 (en) | 2019-08-08 |
Family
ID=67477994
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2018/046459 Ceased WO2019150805A1 (en) | 2018-01-31 | 2018-12-18 | Wind power generation device and wind power generation system |
Country Status (3)
| Country | Link |
|---|---|
| JP (1) | JP7009237B2 (en) |
| TW (1) | TWI729349B (en) |
| WO (1) | WO2019150805A1 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220145856A1 (en) * | 2020-11-06 | 2022-05-12 | IFP Energies Nouvelles | Method of determining the wind speed in the rotor plane of a wind turbine |
| CN119914461A (en) * | 2024-12-16 | 2025-05-02 | 明阳智慧能源集团股份公司 | A method and system for identifying wind shear conditions of a wind turbine generator set |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111044754B (en) * | 2019-12-09 | 2021-07-16 | 大唐山东烟台电力开发有限公司 | A wind condition estimation device for wind turbines |
| 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)
| 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 |
-
2018
- 2018-01-31 JP JP2018014369A patent/JP7009237B2/en active Active
- 2018-12-18 WO PCT/JP2018/046459 patent/WO2019150805A1/en not_active Ceased
-
2019
- 2019-01-23 TW TW108102524A patent/TWI729349B/en active
Patent Citations (3)
| 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 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220145856A1 (en) * | 2020-11-06 | 2022-05-12 | IFP Energies Nouvelles | Method of determining the wind speed in the rotor plane of a wind turbine |
| US11885305B2 (en) * | 2020-11-06 | 2024-01-30 | IFP Energies Nouvelles | Method of determining the wind speed in the rotor plane of a wind turbine |
| CN119914461A (en) * | 2024-12-16 | 2025-05-02 | 明阳智慧能源集团股份公司 | A method and system for identifying wind shear conditions of a wind turbine generator set |
Also Published As
| Publication number | Publication date |
|---|---|
| JP7009237B2 (en) | 2022-01-25 |
| JP2019132178A (en) | 2019-08-08 |
| TW201934870A (en) | 2019-09-01 |
| TWI729349B (en) | 2021-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3613982B1 (en) | Method for controlling operation of a wind turbine | |
| EP3465359B1 (en) | System and method for controlling a dynamic system, in particular a wind turbine | |
| US10202965B2 (en) | Method and system for estimating the surface condition of a rotating blade | |
| US8096762B2 (en) | Wind turbine with pitch control arranged to reduce life shortening loads on components thereof | |
| US20100014969A1 (en) | Wind turbine with blade pitch control to compensate for wind shear and wind misalignment | |
| TWI708893B (en) | Wind power system | |
| US11022100B2 (en) | System and method for controlling wind turbines | |
| US20130259682A1 (en) | Method of rotor-stall prevention in wind turbines | |
| CN113279907B (en) | Method for operating a wind turbine | |
| WO2019150805A1 (en) | Wind power generation device and wind power generation system | |
| EP3619425B1 (en) | Adjustment factor for aerodynamic performance map | |
| CN110067697B (en) | Twist correction factor for aerodynamic performance maps used in wind turbine control | |
| KR20160036214A (en) | Wind turbine generator and control method for wind turbine generator | |
| EP3259472B1 (en) | Control for a wind turbine | |
| JP7684107B2 (en) | Remaining life assessment method and device for wind power generation equipment | |
| CN105041584A (en) | Method for calculating gradient of wind generation set tower | |
| EP4560138A1 (en) | Methods for determining wake behavior and controllers | |
| EP4667741A1 (en) | Control support system for wind power generation apparatus and control support method for wind power generation apparatus | |
| EP4339452A1 (en) | Method for operating a wind turbine, control system and wind turbine | |
| CN113137334B (en) | System and method for optimizing auxiliary loads based on operational use cases | |
| CN119923518A (en) | Wind turbine blade, wind turbine and method for operating a wind turbine | |
| Koopman | Accelerating Wind Energy | |
| Willmer | Wind tunnel tests on a 3m diameter Musgrove windmill | |
| JP2021038729A (en) | Wind power generation system and diagnostic method of wind power generation system | |
| BR112018073343B1 (en) | METHOD AND SYSTEM FOR PREDICTING A PARK-LEVEL ENERGY OUTPUT OF A WIND FARM HAVING A PLURALITY OF WIND TURBINES |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 18904330 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 18904330 Country of ref document: EP Kind code of ref document: A1 |