CN106097146A - A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status - Google Patents
A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status Download PDFInfo
- Publication number
- CN106097146A CN106097146A CN201610443200.3A CN201610443200A CN106097146A CN 106097146 A CN106097146 A CN 106097146A CN 201610443200 A CN201610443200 A CN 201610443200A CN 106097146 A CN106097146 A CN 106097146A
- Authority
- CN
- China
- Prior art keywords
- probability
- parameter
- wind speed
- wind turbines
- wind
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012544 monitoring process Methods 0.000 claims abstract description 11
- 238000012423 maintenance Methods 0.000 claims abstract description 8
- 210000005036 nerve Anatomy 0.000 claims abstract description 4
- 238000002834 transmittance Methods 0.000 claims abstract description 4
- 230000002159 abnormal effect Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 7
- 238000000205 computational method Methods 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000003062 neural network model Methods 0.000 claims description 2
- 210000004218 nerve net Anatomy 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010205 computational analysis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011058 failure modes and effects analysis Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000002633 protecting effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Wind Motors (AREA)
Abstract
The present invention relates to a kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status, comprise the following steps: obtain Wind turbines state parameter by status monitoring with data collecting system, status monitoring parameter is divided into two classes;For device temperature parameter, set up state parameter forecast model based on reverse transmittance nerve network, calculate protection act probability based on prediction residual distribution character, for remaining parameter, calculate protection act probability according to meter and threshold crossing time;By wind energy turbine set operation maintenance data and SCADA data, stoppage in transit number of times and the moment wind speed of stopping transport of Wind turbines are added up, set up the Wind turbines statistics outage model that wind speed is interdependent;Get over limit information in conjunction with stoppage in transit statistical information and state parameter, calculate Wind turbines short-term stoppage in transit probability.This method can evaluate Wind turbines stoppage in transit risk in a short time exactly, increases substantially the accuracy of Wind turbines short-term outage model, and the short term reliability for whole wind energy turbine set is assessed and safety and economic operation provides Technical Reference.
Description
Technical field
The invention belongs to new forms of energy power equipment safety assessment technology field, relate to a kind of meter and the wind turbine of running status
Group short term reliability Forecasting Methodology.
Background technology
Wind turbines is the key equipment of wind energy turbine set, and its reliability has closely connection with the safe and economical operation of wind energy turbine set
System.Owing to long-term work is in severe natural environment, the outage rate of Wind turbines is far above tradition power transmission and transforming equipment.Effectively
Prediction Wind turbines reliability within the following short time, exerts oneself prediction to wind energy turbine set and control has important value.
Utilizing statistical method to obtain equipment dependability parameter is the primary hand of current Wind turbines reliability assessment and prediction
Section.The reliability data of Wind turbines is added up by external existing multiple mechanisms, utilizes statistical data to establish based on can
By property block-scheme method, failure mode and effect analysis method, FTA, appraisal procedure based on markoff process, pattra leaves
The Wind turbines long-term reliability model of this network technique and Monte Carlo method etc..But this class model is difficult to apply and Wind turbines
Short term reliability prediction, main cause includes: 1) reliability of some mechanical and electric information assembly depends not only on product matter
Amount and maintenance condition, also have substantial connection, existing statistical data not to consider these factors with natural environment and operating condition;
2) due to natural environment and operating condition real-time change, static constant reliability statistics index is difficult to reflect that Wind turbines is current
Stoppage in transit risk.
Summary of the invention
In view of this, it is an object of the invention to provide the Wind turbines short term reliability prediction of a kind of meter and running status
Method, the Wind turbines state parameter that the method is obtained with data acquisition (SCADA) system by status monitoring, according to natural ring
Status monitoring parameter is divided into two classes by the dependency in border;For device temperature parameter, set up based on reverse transmittance nerve network
(BPNN) state parameter forecast model, proposes protection act method for calculating probability based on prediction residual distribution character;For
Remaining parameter, proposes meter and the protection act method for calculating probability of threshold crossing time;By wind energy turbine set operation maintenance data and SCADA
Stoppage in transit number of times and the moment wind speed of stopping transport of Wind turbines are added up by data, and the Wind turbines statistics proposing wind speed interdependent is stopped transport
Model.Finally, get over limit information in conjunction with stoppage in transit statistical information and state parameter, propose the calculating side of Wind turbines short-term stoppage in transit probability
Method, improves the short-term stoppage in transit pre-alerting ability of Wind turbines.
For reaching above-mentioned purpose, the present invention provides following technical scheme:
A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status, the method comprises the following steps:
S1: obtain Wind turbines state parameter with data acquisition (SCADA) system by status monitoring, according to natural environment
Dependency status monitoring parameter is divided into two classes: device temperature parameter and remaining state parameter;
S2: for device temperature parameter, sets up state parameter forecast model based on reverse transmittance nerve network (BPNN),
Protection act probability is calculated based on prediction residual distribution character;For remaining parameter, calculate protection according to meter and threshold crossing time dynamic
Make probability;
S3: by wind energy turbine set operation maintenance data and SCADA data to the stoppage in transit number of times of Wind turbines and moment wind of stopping transport
Speed is added up, and sets up the Wind turbines statistics outage model that wind speed is interdependent;
S4: combining stoppage in transit statistical information and state parameter gets over limit information, calculate Wind turbines short-term stoppage in transit probability, carrying out can
By property prediction.
Further, in the method, unit short-term stoppage in transit probability P is obtained based on state parameter out-of-limit protection act probabilityz
For:
In formula, N is the sum of state parameter, PziOut-of-limit protection act probability for state parameter i;
The mode of operation that input is Wind turbines of statistics outage model and wind speed probabilistic forecasting information, when state parameter is got over
When limit protection act probability is less, stoppage in transit probability based on statistical data can reflect to a certain extent owing to other reasons causes
The probability of compressor emergency shutdown, and when parameter out-of-limit protection act probability is bigger, the reference price of stoppage in transit probability based on statistical data
Being worth less, therefore the computational methods of unit short-term stoppage in transit probability are:
P=max (Pt,Pz)
Wherein, PtFor unit short-term stoppage in transit probability based on statistical data, PzStop transport for the short-term out-of-limit based on state parameter
Probability.
Further, in step s 2, for the two class state parameters out-of-limit protection act probability of calculating:
Described device temperature parameter out-of-limit protection act method for calculating probability specifically includes: 1) set up wind speed short-term forecast mould
Type, it was predicted that the wind speed in following 15min;2) temperature parameter forecast model is set up, it was predicted that temperature parameter during each forecasting wind speed value;
Due to wind speed too low time, the impact of output-power fluctuation is negligible by compressor emergency shutdown, therefore do not consider low wind speed and
The shutdown situation that low temperature causes;3) residual distribution characteristic based on forecasting wind speed model and temperature prediction model, calculates temperature
The out-of-limit probability of parameter, owing to predicted time yardstick is relatively big, out-of-limit protection act probability can be approximately considered out-of-limit equal to parameter generally
Rate;
Remaining state parameter described, owing to such state parameter is affected more weak by environmental factors, the out-of-limit guarantor of this type of parameter
Protecting action probability approximation uses threshold crossing time linearly to express, and is shown below:
In formula, t is parameter threshold crossing time, tlimFor setting time;Owing to setting time is shorter, such parameter out-of-limit generally
Rate evaluation time yardstick is less, generally within 1min.
Further, in the out-of-limit protection act probability calculation of device temperature parameter:
Described 1) wind speed Short-term Forecasting Model is set up, it was predicted that the wind speed in following 15min includes: use BP neutral net mould
Type prediction of wind speed, obtains Parameters of Normal Distribution according to the error statistics data estimation of prediction;Out-of-limit for simplifying calculating temperature parameter
Probability, carries out sliding-model control to forecast error, obtains the probability of each error amount according to Parameters of Normal Distribution;
Described 2) temperature parameter forecast model is set up, it was predicted that temperature parameter during each forecasting wind speed value includes: based on BP god
The forecast model of each temperature parameter is set up through network model;In conjunction with the operation characteristic of Wind turbines, use wind speed, ambient temperature with
And the temperature parameter in upper moment is as input parameter;Consider the seasonality of weather conditions, the individual variation of equipment and fault feelings
Condition, choosing state parameter data when unit is properly functioning under each season is that the parametric prediction model of each unit is entered by training sample
Row training modeling;
Described 3) residual distribution characteristic based on forecasting wind speed model and temperature prediction model, calculates the out-of-limit of temperature parameter
Probability includes: when Wind turbines is properly functioning, and the prediction residual symbol of its temperature parameter is the normal distribution of zero from average, and in short-term
The change of interior (in 15min) ambient temperature is inconspicuous, and therefore the out-of-limit probability of the short-term of temperature may largely be determined by wind speed
The temperature parameter value in predictive value and a upper moment;Normal distribution characteristic according to temperature prediction residual error, can get in prediction
The out-of-limit probability of temperature parameter under wind speed:
Oi(vj)=P (Te> Tlim-Tvj)=1-FN(Tlim-Tvj)
In formula, TvjFor prediction of wind speed vjThe predictive value of lower temperature parameter i, TlimFor the higher limit of temperature parameter, FN() is
The distribution function of temperature prediction error;
When prediction residual falls outside normal distribution 99% confidence interval, it is calculated as once predicting inefficacy;If at one hour
The most at least occur that 3 predictions were lost efficacy, then it is assumed that temperature parameter is abnormal;Data conduct during owing to using Wind turbines properly functioning
The state parameter that the forecast model that sample training obtains is difficult to being in abnormal condition carries out Accurate Prediction, therefore this method root
Predict the outcome according to the intensity of anomaly correction of parameter, the average of former residual error normal distribution is set to the forecast error at upper a moment, and then
Obtain unit when being in abnormal conditions, it was predicted that the out-of-limit probability of the temperature parameter under wind speed:
Oi(vj)=1-FN(Tlim-Tvj-εt-1)
In formula, εt-1Forecast error for upper a moment;
The out-of-limit probabilistic forecasting value of temperature parameter is as follows:
In formula, N is that the centrifugal pump of forecasting wind speed error is counted, P (vj) represent that wind speed is v in a short timejProbability;Prediction wind
Speed value vjWith P (vj) can be obtained by forecasting wind speed result and forecast error.
Further, in step s3, by wind energy turbine set SCADA data and the stoppage in transit of operation maintenance data statistics Wind turbines
Number of times and the wind speed in moment of stopping transport, carry out subregion to wind speed with the interval of 1m/s, and under each wind speed interval, Wind turbines outage rate can
It is calculated by following formula:
In formula, NviFor wind speed viIn the case of the number of times of always stopping transport of all Wind turbines, T in wind energy turbine setviFor all units
Cumulative operation time;
The Wind turbines short-term stoppage in transit probability of meter and wind speed is shown below:
Pt(t, v)=1-eλ(v)t。
The beneficial effects of the present invention is: Wind turbines short term reliability Forecasting Methodology of the present invention is easily programmed reality
Existing, it is possible to the fault stoppage in transit accident evaluating Wind turbines stoppage in transit risk in a short time, the most only Wind turbines exactly carries
For accurate early warning information, increase substantially the accuracy of Wind turbines short-term outage model, be also that the short-term of whole wind energy turbine set can
Technical Reference is provided by property assessment and safety and economic operation.
Accompanying drawing explanation
In order to make the purpose of the present invention, technical scheme and beneficial effect clearer, the present invention provides drawings described below to carry out
Illustrate:
The Wind turbines short term reliability estimation flow that Fig. 1 provides for the present invention;
Fig. 2 is the out-of-limit method for calculating probability of device temperature parameter;
Fig. 3 is the outage rate under each wind speed;
Fig. 4 is the unit outage probability in embodiment.
Detailed description of the invention
Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described in detail.
Fig. 1 is Wind turbines short term reliability estimation flow figure.Wind turbines state parameter abnormal information, unit Working mould
Formula and the input parameter that wind speed probabilistic forecasting information is Wind turbines short term reliability assessment models.Owing to state parameter occurs more
Limit situation may directly result in compressor emergency shutdown, it is also possible to provides alarm signal then to be determined whether to shut down by staff,
State parameter generation is out-of-limit will result directly in unit outage.Obviously any state parameter is out-of-limit, and unit all will be caused to stop
Fortune, therefore, obtains unit short-term stoppage in transit probability P based on state parameter out-of-limit protection act modelzFor:
In formula, N is the sum of state parameter, PziOut-of-limit protection act probability for state parameter i.
The mode of operation that input is Wind turbines of statistics outage model and wind speed probabilistic forecasting information.When state parameter is got over
When limit protection act probability is less, stoppage in transit probability based on statistical data can reflect to a certain extent owing to other reasons causes
The probability of compressor emergency shutdown, and when parameter out-of-limit protection act probability is bigger, the reference price of stoppage in transit probability based on statistical data
It is worth less.Therefore the computational methods of unit short-term stoppage in transit probability are:
P=max (Pt,Pz)
Wherein, PtFor unit short-term stoppage in transit probability based on statistical data, PzStop transport for the short-term out-of-limit based on state parameter
Probability.
1, state parameter out-of-limit protection act probability
Wind turbines SCADA state parameter is as shown in table 1, can be according to the dependency of each state parameter and natural environment by this
A little parameters are divided into two classes, and first kind status monitoring parameter is notable by wind speed and ambient temperature effect, predominantly device temperature parameter;
Equations of The Second Kind status monitoring parameter is more weak with the dependency relation of natural environment.Next will be respectively directed to this two classes state parameter calculate
Out-of-limit protection act probability.
Table 1 status monitoring parameter
1) first kind parameter out-of-limit protection act method for calculating probability:
The computational methods of first kind state parameter out-of-limit protection act probability are as in figure 2 it is shown, calculating process can be divided into three steps:
(1) wind speed in the following 15min of wind speed short-term forecast prediction is set up.(2) temperature parameter forecast model is set up, it was predicted that each wind speed is pre-
Temperature parameter during measured value.Due to wind speed too low time, the impact of output-power fluctuation is negligible, therefore by compressor emergency shutdown
Do not consider low wind speed and shutdown situation that low temperature causes.(3) based on forecasting wind speed model and the residual error of temperature prediction model
Distribution character, calculates the out-of-limit probability of temperature parameter, and owing to predicted time yardstick is relatively big, out-of-limit protection act probability can approximate to be recognized
For equal to the out-of-limit probability of parameter.
Short-term wind speed probabilistic forecasting:
Use BP Neural Network model predictive wind speed, obtain normal distribution ginseng according to the error statistics data estimation of prediction
Number.Calculate the out-of-limit probability of temperature parameter for simplifying, forecast error is carried out sliding-model control, obtains respectively according to Parameters of Normal Distribution
The probability of error amount.Table 2 show the forecast error centrifugal pump of Single-step Prediction in 15min and corresponding probability.
The table 2 forecasting wind speed probability of error
Error/m/s | Probability |
0 | 0.234 |
±0.5 | 0.197 |
±1 | 0.118 |
±1.5 | 0.05 |
±2 | 0.018 |
Temperature parameter probabilistic forecasting:
Forecast model based on each temperature parameter of BP Establishment of Neural Model.In conjunction with the operation characteristic of Wind turbines, adopt
With the temperature parameter in wind speed, ambient temperature and a upper moment as input parameter.Consider the seasonality of weather conditions, equipment
Individual variation and failure condition, choosing state parameter data when unit is properly functioning under each season is that training sample is to each unit
Parametric prediction model be trained modeling.
The out-of-limit probability of temperature parameter:
When Wind turbines is properly functioning, the prediction residual of its temperature parameter generally accord with from average be the normal distribution of zero.And
In short time, the change of (in 15min) ambient temperature is inconspicuous, and therefore the out-of-limit probability of the short-term of temperature may largely be determined by
The predictive value of wind speed and the temperature parameter value in a upper moment.Normal distribution characteristic according to temperature prediction residual error, can get in advance
The out-of-limit probability of temperature parameter under the wind speed surveyed:
Oi(vj)=P (Te> Tlim-Tvj)=1-FN(Tlim-Tvj)
In formula, TvjFor prediction of wind speed vjThe predictive value of lower temperature parameter i, TlimFor the higher limit of temperature parameter, FN() is
The distribution function of temperature prediction error.
When prediction residual falls outside normal distribution 99% confidence interval, it is calculated as once predicting inefficacy.If at one hour
The most at least occur that 3 predictions were lost efficacy, then it is assumed that temperature parameter is abnormal.Data conduct during owing to using Wind turbines properly functioning
The state parameter that the forecast model that sample training obtains is difficult to being in abnormal condition carries out Accurate Prediction, therefore this patent root
Predict the outcome according to the intensity of anomaly correction of parameter, the average of former residual error normal distribution is set to the forecast error at upper a moment, and then
Obtain unit when being in abnormal conditions, it was predicted that the out-of-limit probability of the temperature parameter under wind speed:
Oi(vj)=1-FN(Tlim-Tvj-εt-1)
In formula, εt-1Forecast error for upper a moment.
The out-of-limit probabilistic forecasting value of temperature parameter is as follows:
In formula, N is that the centrifugal pump of forecasting wind speed error is counted, P (vj) represent that wind speed is v in a short timejProbability.Prediction wind
Speed value vjWith P (vj) can be obtained by forecasting wind speed result and forecast error.
2) Equations of The Second Kind parameter out-of-limit protection act method for calculating probability:
Owing to Equations of The Second Kind state parameter is affected more weak by environmental factors, this type of parameter out-of-limit protection act probability can approximate
Use threshold crossing time linearly to express, be shown below:
In formula, t is parameter threshold crossing time, tlimFor setting time.Owing to setting time is shorter, Equations of The Second Kind parameter out-of-limit
Probability assessment time scale is less, generally within 1min.
2, statistics outage model
Fig. 3 is the outage rate under each wind speed, by wind energy turbine set SCADA data and operation maintenance data statistics Wind turbines
Stoppage in transit number of times and the wind speed in moment of stopping transport, carry out subregion to wind speed with the interval of 1m/s, and under each wind speed interval, Wind turbines is stopped transport
Rate can be calculated by following formula::
In formula, NviFor wind speed viIn the case of the number of times of always stopping transport of all Wind turbines, T in wind energy turbine setviFor all units
Cumulative operation time.The Wind turbines short-term stoppage in transit probability of meter and wind speed is shown below:
Pt(t, v)=1-eλ(v)t
Embodiment:
The present embodiment causes compressor emergency shutdown accident with one of certain wind energy turbine set domestic because first kind state parameter is out-of-limit
As a example by Wind turbines, the short-term stoppage in transit probability of unit is studied, by reliable to meter and the short-term of operating states of the units
Property Forecasting Methodology and the relative analysis of short term reliability Forecasting Methodology result based on statistical data, demonstrate the carried side of the present invention
The effectiveness of method and accuracy.
This unit on June 17th, 2014 about 17:00 stop transport because of dynamo bearing B temperature over-range.Stopping transport
Front in several hours, mean wind speed fluctuates near rated value 11.5m/s, and ambient temperature is in the state being stepped up.Machine
Group dynamo bearing B temperature starts occur extremely in early June, and actual value is mostly higher than predictive value.Therefore, unit exposes for a long time
In high temperature, high wind speed natural environment under, the lasting rising of bearing B temperature has ultimately resulted in the generation of compressor emergency shutdown accident.
The probabilistic forecasting result of table 3 wind speed
Out-of-limit probability under each wind speed of table 4
With the data of the moment 17:00 that stops transport, compressor emergency shutdown probability is carried out computational analysis.The mean wind speed prediction of this period
Really definite value is 11.2m/s.According to the forecasting wind speed probability of error, the probabilistic forecasting value of wind speed can be obtained, as shown in table 3.Each wind speed
Under temperature prediction determine that value and corresponding out-of-limit probability are as shown in table 4.Respectively by meter proposed by the invention and running status
Stoppage in transit probabilistic model and be based only upon the outage model of statistical data and calculate the stoppage in transit probability in each moment of this unit, such as Fig. 4
Shown in, visible by relative analysis, unit outage model based on running status can the stoppage in transit in a short time of significant reaction unit
Risk, although and the outage model based on the statistical data high wind speed period before stoppage in transit has shown that of a relatively high stoppage in transit is general
Rate, but still far below this method.Therefore, the outage model that the present invention proposes can effectively reflect unit stoppage in transit wind in a short time
Danger.
Finally illustrate, preferred embodiment above only in order to technical scheme to be described and unrestricted, although logical
Cross above preferred embodiment the present invention to be described in detail, it is to be understood by those skilled in the art that can be
In form and it is made various change, without departing from claims of the present invention limited range in details.
Claims (5)
1. a meter and the Wind turbines short term reliability Forecasting Methodology of running status, it is characterised in that: the method includes following
Step:
S1: obtain Wind turbines state parameter with data acquisition (SCADA) system by status monitoring, according to the phase of natural environment
Status monitoring parameter is divided into two classes by closing property: device temperature parameter and remaining state parameter;
S2: for device temperature parameter, sets up state parameter forecast model based on reverse transmittance nerve network (BPNN), based on
Prediction residual distribution character calculates protection act probability;For remaining parameter, calculate protection act according to meter and threshold crossing time general
Rate;
S3: stoppage in transit number of times and the moment wind speed of stopping transport of Wind turbines are entered by wind energy turbine set operation maintenance data and SCADA data
Row statistics, sets up the Wind turbines statistics outage model that wind speed is interdependent;
S4: combine stoppage in transit statistical information and state parameter gets over limit information, calculate Wind turbines short-term stoppage in transit probability, carry out reliability
Prediction.
A kind of meter the most according to claim 1 and the Wind turbines short term reliability Forecasting Methodology of running status, its feature
It is: in the method, obtains unit short-term stoppage in transit probability P based on state parameter out-of-limit protection act probabilityzFor:
In formula, N is the sum of state parameter, PziOut-of-limit protection act probability for state parameter i;
The mode of operation that input is Wind turbines of statistics outage model and wind speed probabilistic forecasting information, as the out-of-limit guarantor of state parameter
Protect action probability less time, stoppage in transit probability based on statistical data can reflect to a certain extent owing to other reasons causes unit
The probability shut down, and when parameter out-of-limit protection act probability is bigger, the reference value of stoppage in transit probability based on statistical data is relatively
Little, therefore the computational methods of unit short-term stoppage in transit probability are:
P=max (Pt,Pz)
Wherein, PtFor unit short-term stoppage in transit probability based on statistical data, PzStop transport for the short-term out-of-limit based on state parameter general
Rate.
A kind of meter the most according to claim 2 and the Wind turbines short term reliability Forecasting Methodology of running status, its feature
It is: in step s 2, for the two class state parameters out-of-limit protection act probability of calculating:
Described device temperature parameter out-of-limit protection act method for calculating probability specifically includes: 1) set up wind speed Short-term Forecasting Model,
Wind speed in the following 15min of prediction;2) temperature parameter forecast model is set up, it was predicted that temperature parameter during each forecasting wind speed value;By
When wind speed is too low, the impact of output-power fluctuation is negligible by compressor emergency shutdown, does not therefore consider low wind speed and low
The shutdown situation that temperature causes;3) residual distribution characteristic based on forecasting wind speed model and temperature prediction model, calculates temperature ginseng
The out-of-limit probability of number, owing to predicted time yardstick is relatively big, out-of-limit protection act probability can be approximately considered equal to the out-of-limit probability of parameter;
Remaining state parameter described, owing to such state parameter is affected more weak by environmental factors, the out-of-limit protection of this type of parameter is moved
Making probability approximation uses threshold crossing time linearly to express, and is shown below:
In formula, t is parameter threshold crossing time, tlimFor setting time;Owing to setting time is shorter, the out-of-limit probability of such parameter is commented
Estimate time scale less, generally within 1min.
A kind of meter the most according to claim 3 and the Wind turbines short term reliability Forecasting Methodology of running status, its feature
It is: in the out-of-limit protection act probability calculation of device temperature parameter:
Described 1) wind speed Short-term Forecasting Model is set up, it was predicted that the wind speed in following 15min includes: use BP neural network model pre-
Survey wind speed, obtain Parameters of Normal Distribution according to the error statistics data estimation of prediction;The out-of-limit probability of temperature parameter is calculated for simplifying,
Forecast error is carried out sliding-model control, obtains the probability of each error amount according to Parameters of Normal Distribution;
Described 2) temperature parameter forecast model is set up, it was predicted that temperature parameter during each forecasting wind speed value includes: based on BP nerve net
The forecast model of each temperature parameter set up by network model;In conjunction with the operation characteristic of Wind turbines, use wind speed, ambient temperature and on
The temperature parameter in one moment is as input parameter;Consider the seasonality of weather conditions, the individual variation of equipment and failure condition, choosing
Taking state parameter data when unit is properly functioning under each season is that the parametric prediction model of each unit is instructed by training sample
Practice modeling;
Described 3) residual distribution characteristic based on forecasting wind speed model and temperature prediction model, calculates the out-of-limit probability of temperature parameter
Including: when Wind turbines is properly functioning, the prediction residual symbol of its temperature parameter is the normal distribution of zero from average, and in the short time
The change of (in 15min) ambient temperature is inconspicuous, and therefore the out-of-limit probability of the short-term of temperature may largely be determined by the pre-of wind speed
Measured value and the temperature parameter value in a upper moment;Normal distribution characteristic according to temperature prediction residual error, can get the wind speed in prediction
The out-of-limit probability of lower temperature parameter:
Oi(vj)=P (Te> Tlim-Tvj)=1-FN(Tlim-Tvj)
In formula, TvjFor prediction of wind speed vjThe predictive value of lower temperature parameter i, TlimFor the higher limit of temperature parameter, FN() is temperature
The distribution function of forecast error;
When prediction residual falls outside normal distribution 99% confidence interval, it is calculated as once predicting inefficacy;If in one hour extremely
Occur that 3 predictions were lost efficacy less, then it is assumed that temperature parameter is abnormal;Data during owing to using Wind turbines properly functioning are as sample
The state parameter that the forecast model that training obtains is difficult to being in abnormal condition carries out Accurate Prediction, and therefore this method is according to ginseng
The intensity of anomaly correction of number predicts the outcome, and the average of former residual error normal distribution is set to the forecast error at upper a moment, and then obtains
When unit is in abnormal conditions, it was predicted that the out-of-limit probability of the temperature parameter under wind speed:
Oi(vj)=1-FN(Tlim-Tvj-εt-1)
In formula, εt-1Forecast error for upper a moment;
The out-of-limit probabilistic forecasting value of temperature parameter is as follows:
In formula, N is that the centrifugal pump of forecasting wind speed error is counted, P (vj) represent that wind speed is v in a short timejProbability;Prediction of wind speed value
vjWith P (vj) can be obtained by forecasting wind speed result and forecast error.
A kind of meter the most according to claim 4 and the Wind turbines short term reliability Forecasting Methodology of running status, its feature
It is: in step s3, by wind energy turbine set SCADA data and the stoppage in transit number of times of operation maintenance data statistics Wind turbines and stoppage in transit
The wind speed in moment, carries out subregion to wind speed with the interval of 1m/s, and under each wind speed interval, Wind turbines outage rate can be calculated by following formula
Obtain:
In formula, NviFor wind speed viIn the case of the number of times of always stopping transport of all Wind turbines, T in wind energy turbine setviAdding up for all units
Working time;
The Wind turbines short-term stoppage in transit probability of meter and wind speed is shown below:
Pt(t, v)=1-eλ(v)t。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443200.3A CN106097146A (en) | 2016-06-20 | 2016-06-20 | A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443200.3A CN106097146A (en) | 2016-06-20 | 2016-06-20 | A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106097146A true CN106097146A (en) | 2016-11-09 |
Family
ID=57237349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610443200.3A Pending CN106097146A (en) | 2016-06-20 | 2016-06-20 | A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106097146A (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106651035A (en) * | 2016-12-25 | 2017-05-10 | 杭州博烁晟斐智能科技有限公司 | Method and system for prediction of iron tower verticality fault probability caused by strong breeze weather |
CN107944706A (en) * | 2017-11-27 | 2018-04-20 | 国网内蒙古东部电力有限公司检修分公司 | A kind of safe method for commissioning of ultra-high-tension power transmission line and system |
CN109376877A (en) * | 2018-10-11 | 2019-02-22 | 华自科技股份有限公司 | Equipment O&M method for early warning, device, computer equipment and storage medium |
CN109685371A (en) * | 2018-12-25 | 2019-04-26 | 华能陕西定边电力有限公司 | Dynamic based on Bayesian network generally weighs running of wind generating set state comprehensive estimation method |
CN110023850A (en) * | 2016-12-06 | 2019-07-16 | 西门子股份公司 | Method and control device for control technology system |
CN110570122A (en) * | 2019-09-10 | 2019-12-13 | 重庆大学 | Offshore wind power plant reliability assessment method considering wind speed seasonal characteristics and current collection system element faults |
CN110739679A (en) * | 2019-08-27 | 2020-01-31 | 华北电力大学 | Flexible multi-state switch reliability analysis method based on Bayesian network method |
CN111080039A (en) * | 2020-03-17 | 2020-04-28 | 浙江上风高科专风实业有限公司 | Fan cluster fault prediction method and system |
CN114308358A (en) * | 2022-03-17 | 2022-04-12 | 山东金有粮脱皮制粉设备有限公司 | Safe operation monitoring system of corncob grinding device |
CN115596620A (en) * | 2022-12-15 | 2023-01-13 | 深圳鹏锐信息技术股份有限公司(Cn) | Wind generating set fault intelligent analysis method and system based on artificial intelligence |
US11635060B2 (en) | 2021-01-20 | 2023-04-25 | General Electric Company | System for operating a wind turbine using cumulative load histograms based on actual operation thereof |
US11661919B2 (en) | 2021-01-20 | 2023-05-30 | General Electric Company | Odometer-based control of a wind turbine power system |
CN116523349A (en) * | 2023-05-19 | 2023-08-01 | 北京协合运维风电技术有限公司 | Wind power station reliability analysis method and system |
US11728654B2 (en) | 2021-03-19 | 2023-08-15 | General Electric Renovables Espana, S.L. | Systems and methods for operating power generating assets |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427005A (en) * | 2016-01-18 | 2016-03-23 | 重庆大学 | Operation risk assessment method of wind power station |
-
2016
- 2016-06-20 CN CN201610443200.3A patent/CN106097146A/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427005A (en) * | 2016-01-18 | 2016-03-23 | 重庆大学 | Operation risk assessment method of wind power station |
Non-Patent Citations (1)
Title |
---|
李剑等: "一种计及运行状态的风电机组短期停运模型", 《中国电机工程学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110023850A (en) * | 2016-12-06 | 2019-07-16 | 西门子股份公司 | Method and control device for control technology system |
US11340564B2 (en) | 2016-12-06 | 2022-05-24 | Siemens Aktiengesellschaft | Method and control device for controlling a technical system |
CN106651035A (en) * | 2016-12-25 | 2017-05-10 | 杭州博烁晟斐智能科技有限公司 | Method and system for prediction of iron tower verticality fault probability caused by strong breeze weather |
CN107944706A (en) * | 2017-11-27 | 2018-04-20 | 国网内蒙古东部电力有限公司检修分公司 | A kind of safe method for commissioning of ultra-high-tension power transmission line and system |
CN109376877A (en) * | 2018-10-11 | 2019-02-22 | 华自科技股份有限公司 | Equipment O&M method for early warning, device, computer equipment and storage medium |
CN109376877B (en) * | 2018-10-11 | 2021-01-26 | 华自科技股份有限公司 | Equipment operation and maintenance early warning method and device, computer equipment and storage medium |
CN109685371B (en) * | 2018-12-25 | 2021-07-13 | 华能陕西定边电力有限公司 | Bayesian network-based dynamic probability wind turbine generator running state comprehensive evaluation method |
CN109685371A (en) * | 2018-12-25 | 2019-04-26 | 华能陕西定边电力有限公司 | Dynamic based on Bayesian network generally weighs running of wind generating set state comprehensive estimation method |
CN110739679A (en) * | 2019-08-27 | 2020-01-31 | 华北电力大学 | Flexible multi-state switch reliability analysis method based on Bayesian network method |
CN110570122B (en) * | 2019-09-10 | 2023-03-14 | 重庆大学 | Offshore wind power plant reliability assessment method considering wind speed seasonal characteristics and current collection system element faults |
CN110570122A (en) * | 2019-09-10 | 2019-12-13 | 重庆大学 | Offshore wind power plant reliability assessment method considering wind speed seasonal characteristics and current collection system element faults |
CN111080039A (en) * | 2020-03-17 | 2020-04-28 | 浙江上风高科专风实业有限公司 | Fan cluster fault prediction method and system |
CN111080039B (en) * | 2020-03-17 | 2023-10-20 | 浙江上风高科专风实业有限公司 | Fan cluster fault prediction method and system |
US11635060B2 (en) | 2021-01-20 | 2023-04-25 | General Electric Company | System for operating a wind turbine using cumulative load histograms based on actual operation thereof |
US11661919B2 (en) | 2021-01-20 | 2023-05-30 | General Electric Company | Odometer-based control of a wind turbine power system |
US11728654B2 (en) | 2021-03-19 | 2023-08-15 | General Electric Renovables Espana, S.L. | Systems and methods for operating power generating assets |
CN114308358A (en) * | 2022-03-17 | 2022-04-12 | 山东金有粮脱皮制粉设备有限公司 | Safe operation monitoring system of corncob grinding device |
CN114308358B (en) * | 2022-03-17 | 2022-05-27 | 山东金有粮脱皮制粉设备有限公司 | Safe operation monitoring system of corncob grinding device |
CN115596620A (en) * | 2022-12-15 | 2023-01-13 | 深圳鹏锐信息技术股份有限公司(Cn) | Wind generating set fault intelligent analysis method and system based on artificial intelligence |
CN115596620B (en) * | 2022-12-15 | 2023-03-07 | 深圳鹏锐信息技术股份有限公司 | Wind generating set fault intelligent analysis method and system based on artificial intelligence |
CN116523349A (en) * | 2023-05-19 | 2023-08-01 | 北京协合运维风电技术有限公司 | Wind power station reliability analysis method and system |
CN116523349B (en) * | 2023-05-19 | 2024-01-23 | 北京协合运维风电技术有限公司 | Wind power station reliability analysis method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106097146A (en) | A kind of meter and the Wind turbines short term reliability Forecasting Methodology of running status | |
CN102063651B (en) | Urban power grid risk evaluation system based on on-line data acquisition | |
CN103440400B (en) | Power system short-term risk determination method taking disaster factors into account | |
EP2537222B1 (en) | Energy generating system and control thereof | |
Guo et al. | Wind turbine generator condition-monitoring using temperature trend analysis | |
CN102022264B (en) | System and method for wind turbine health management | |
KR101096793B1 (en) | data collecting method for detection and on-time warning system of industrial process | |
CN105894133A (en) | Wind turbine component maintenance and spare part demand forecasting method | |
CN103246939B (en) | Safe operation of electric network risk case on-line identification method based on security margin | |
CN106251045A (en) | Distribution network reliability appraisal procedure based on multiple leading factor | |
CN103218754B (en) | A kind of risk checking method of Forming Electrical Dispatching Command Tickets and device | |
CN106203800A (en) | A kind of power distribution network operational reliability index extraction method | |
CN101923685A (en) | System and method for deciding power shedding load based on line breaking fault rate prediction | |
CN104599023A (en) | Typhoon weather transmission line time-variant reliability calculation method and risk evaluation system | |
CN102915515A (en) | Cascade outage recognition and risk evaluation method of power grid | |
KR102496917B1 (en) | Failure prediction system based on failure prediction model | |
CN102663522A (en) | On-line risk evaluation method of power grid | |
CN104268693A (en) | Power grid dispatching determining method and system | |
CN117728585B (en) | Distributed photovoltaic power station state prediction method and system | |
Wu et al. | Preventive control strategy for an island power system that considers system security and economics | |
Wan | Primer on wind power for utility applications | |
CN103605910A (en) | Reliability evaluation method for single protection device based on consideration of covert fault | |
CN106022594A (en) | Power system online operation safety risk assessment method based on extreme value theory | |
Zeng et al. | Impact of intermittent wind distributed generation on the reliability of distribution system | |
CN107886220A (en) | A kind of distribution line risk probability appraisal procedure based on historical factor analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161109 |
|
RJ01 | Rejection of invention patent application after publication |