• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Detecting Icing on the Blades of a Wind Turbine Using a Deep Neural Network

    2023-01-24 02:51:20TingshunLiJiaohuiXuZesanLiuDadiWangandWenTan

    Tingshun Li,Jiaohui Xu,★,Zesan Liu,Dadi Wang and Wen Tan

    1School of Control and Computer Engineering,North China Electric Power University,Beijing,102206,China

    2R&D,State Grid Information&Telecommunication Group Co.,Beijing,102211,China

    ABSTRACT The blades of wind turbines located at high latitudes are often covered with ice in late autumn and winter,where this affects their capacity for power generation as well as their safety.Accurately identifying the icing of the blades of wind turbines in remote areas is thus important,and a general model is needed to this end.This paper proposes a universal model based on a Deep Neural Network(DNN)that uses data from the Supervisory Control and Data Acquisition(SCADA)system.Two datasets from SCADA are first preprocessed through undersampling,that is,they are labeled,normalized,and balanced.The features of icing of the blades of a turbine identified in previous studies are then used to extract training data from the training dataset.A middle feature is proposed to show how a given feature is correlated with icing on the blade.Performance indicators for the model,including a reward function,are also designed to assess its predictive accuracy.Finally,the most suitable model is used to predict the testing data,and values of the reward function and the predictive accuracy of the model are calculated.The proposed method can be used to relate continuously transferred features with a binary status of icing of the blades of the turbine by using variables of the middle feature.The results here show that an integrated indicator system is superior to a single indicator of accuracy when evaluating the prediction model.

    KEYWORDS DNN;predicting blade icing;SCADA data;wind power;reward function

    1 Introduction

    China and the United States have recently released the Joint Glasgow Declaration on Enhancing Climate Action at the United Nations (UN) climate conference,known as COP26,in Scotland [1].China had set itself the goals of achieving a carbon peak and carbon neutralization prior to this[2].As it is a major form of renewable energy,advances in the use of wind power are expected to continue towards this end(wind power generation in China had increased by 34.4%as of May 2021[3]).Most of these wind turbines have been installed in high-latitude areas,where the climate is severe.The blades of these turbines are often covered with a thick layer of ice in late autumn and winter,as shown in Fig.1.The icing of the blades is a critical problem that affects the capacity of the turbine for power generation as well as its safety[4–6].Some researchers have claimed that the power output of a wind turbine declines by 27%in the case of icing[7].

    Figure 1:A wind turbine covered with ice

    Many solutions have been proposed to the problem of ice buildup upon the blades of a wind turbine.Rizk et al.[8]proposed using hyperspectral imaging to detect the icing of the blades of a wind turbine,where this requires the construction of a hyperspectral imaging system.Research on the icing of the blades of wind turbines has focused on the analysis of data from SCADA(Supervisory Control and Data Acquisition)without using any additional devices.For example,engineers determine whether the blades of a turbine are covered with ice by monitoring SCADA data [9].Some researchers have designed machine learning models to predict icing by using this data.Xiong et al.[10]built a model to predict the icing of the blades of a wind turbine in the short term based on Bi-LSTM and the SVM.Xu et al.[7]also proposed a prediction model based on the SVM,and demonstrated how to obtain better parameters for it by using particle swarm optimization.Such models can be used to detect the status of wind turbines,but most of them run into the problem of overfitting and are thus unsuitable for practical use.Moreover,the accuracy of their predictions decreases when they are transferred to the operating environments of wind turbines other than the ones for which they were designed.A universal machine learning-based model is thus needed for predicting ice buildup upon the blades of wind turbines.

    In this paper,the authors first select features from SCADA data that are closely correlated with ice buildup upon the blades of a wind turbine.Big data from SCADA are then prepared for use in machine learning,and are labeled and preprocessed.Following this,a model is built based on a deep neural network(DNN).This model features an integrated indicator system that can not only measure the accuracy of the model,but also its ability to be used on datasets of wind turbines other than the specific one on which it is based.The remainder of this paper is organized as follows: Section 2 details the data preprocessing,including data labeling and normalization,the balancing of the data distribution,and the division of the data into groups.Section 3 discusses feature selection based on previous studies in the area,and Section 4 introduces the model based on the DNN and parameters of the model.Section 5 describes the experimental process used to verify the proposed model and its results.It involves building an integrated indicator system that can assess the accuracy of the model,a description of the experimental process,and a comparative analysis of the results.Finally,Section 6 concludes by summarizing the limitations of the proposed method and its implications for future work in the area.

    2 Data Preprocessing

    2.1 Data Source

    A model based on machine learning involves the use of big data.Two datasets of SCADA,corresponding to two wind turbines in a wind farm in North China,were used.Data from one,called turbine#1,was used to train the model.The dataset contained more than 190,000 items.Data from the other,called turbine#2,was used to validate the model.Its dataset contained nearly 170,000 data items.The distributions of data in the two sets are presented in Table 1.

    Table 1: Distribution of datasets of wind turbines

    The sources of data from SCADA have 26 measurements,as shown in Table 2.

    Table 2: SCADA measurements of wind turbines

    Table 2 (continued)

    These SCADA data were automatically generated by using sensors and programs in the empirical environment.They were highly structured and time-serialized,and thus were suitable for analysis and data mining.

    2.2 Balancing the Data

    The distribution of the data shown in Table 1 indicates that the datasets contained a significantly greater amount of data on icing-free blades than those representing their icing.This imbalance is reasonable because most wind turbines operate normally most of the time.However,most machine learning algorithms used for classification are sensitive to the distribution of the training samples,which means that this imbalance would lead to a poor model that is complex and inaccurate.It was thus necessary to balance the distribution of the sampled data.Oversampling,undersampling,and other sampling methods based on both are recommended to balance the data.

    Methods to balance the data have been provided in [7,11].Among the normal data,we first eliminated the icing-free data and then eliminated redundant data.For data representing icing,we first removed data representing heavy icing and then eliminated redundant data.The resulting amount of icing-free data in the sample dataset was approximately equal to the data representing icing.

    2.3 Normalizing the Data

    The magnitudes of different monitoring points are in general different.Some are very large and others are quite small.This numerical size can affect the weights of the model,and can even cause some features that are related to icing to be ignored during model training.This reduces the predictive accuracy of the model.To reduce the impact of this numerical difference,the sample data needs to be normalized[12].

    There are many methods to normalize data—for example,min-max normalization and the Z-score[12].Given that the model will be used in empirical environments,where it is impossible to compute the standard deviation,the min-max standard method is more suitable for use.

    In the above,“min”is the minimum value of variable x while“max”is its maximum value,andx*is its normalized value.

    The values of the characteristic variable x can be mapped to the range[0,1]by using Eq.(1).

    3 Selecting Features Correlated with Icing

    The process of ice buildup upon the blades of a turbine is complex,and is directly related to such environmental factors as the ambient temperature,wind speed,and humidity.It is challenging to establish a mathematical model to predict the ice buildup upon wind turbine blades by using only these factors.Moreover,the aerodynamic performance of the blades of wind turbines changes when they are frozen.For example,the power generated at the same wind speed is reduced.These factors can be used to predict icing,but it is difficult to identify the variables related to the icing of the blades.

    There are 26 measurements in SCADA for wind turbines,as shown in Table 2,and some of them are correlated with the icing of the blades.Only a model formulated by using the features related to icing can accurately predict it when transferred to other turbines.

    3.1 Transferring the Features

    Even if the manner in which a given feature affects icing is known,that in which a group of features impacts icing cannot be determined by simply aggregating the effects of all features due to the influence that the variables have on one another.Transfer learning is a useful solution to the problem of an insufficient number of samples.Pandeya et al.provided a model to detect the needs of pet cats according to their emitted sounds by transferring some features of such sounds from their previous research[13].This suggests that somefeatures of the icing of blades can also be transferred over from similar past studies and utilized as well.

    Wei et al.[9] provided features related to the icing of the blades based on practical experience.Li et al.[11]proposed a method to identify features correlated with the icing of the blades based on logistic regression and experiments.We use these two studies to extract some features related to icing,as shown in Table 3.

    Table 3: Variables correlated with icing of the blades of a turbine

    Some of the feature variables in Table 3 are copied from Table 2 while the others were calculated by using the variables in Table 2.

    3.2 Feature Correlation Analysis

    It was necessary to check each feature in Table 3 to determine the extent to which it was correlated with the icing of the blades.A chart is a visible and direct way to show the relationships between the relevant variables and the status of the icing of blades.However,whether or not icing has occurred is a binary proposition in time rather than a continuous one,such as wind speed,ambient temperature,and power output.We thus seek a“middle feature”that is continuous in time and closely correlated with the icing of the blades,and use it to relate the continuous features in time with the icing of the blades.

    The rotation of the blades of the turbine generates power.The stronger the wind is,the more power is generated.When the blades are covered with ice,the aerodynamic characteristics of the wind turbine change and its power output decreases under the same wind-induced force.Goshima et al.claimed that the power output at a rated wind speed in case of icing is smaller than that without icing,as is the maximum thrust of the rotor,whereas the wind speed at which the turbine’s rated power can be reached in case of icing is higher than that without icing[14].The power output of the wind turbine thus directly reflects whether icing has occurred on it,and is a suitable middle feature.

    By taking power as the middle feature,the relation between a given variable and the power output can reflect the status of the wind turbine.An example is the variable of wind speed,which is one of the most apparent features related to the power output.The relationship between wind speed,power,and the icing of the blades is clearly presented in Fig.2.Certainly,the data in the below figures are normalized to ignore the influence of the measurement scale.

    Figure 2: (Continued)

    Figure 2:Relation between wind speed and power output of a wind turbine.(a)Distribution of wind speed with changes in time(gray points)and power(red points)(the upper horizontal blue lines in the figure represent periods of icing of the wind turbine).(b)Partial magnification of the diagram of the relationship between wind speed and power(the circled part represents no icing while the rectangular part represents icing)

    The relation between power and Wind_Speed_Face in Fig.2a indicates that the power output of the wind turbine is approximately linearly related to wind speed when there is no icing on its blades,but this linear relationship is lost in case of icing.Part of Fig.2a is enlarged to show the relationship between wind speed and power in Fig.2b.In the case of icing,the power output of the turbine departs further in relation to wind speed than in the icing-free case.

    To clearly show the relationship between wind speed and power,a comparison between the cases with and without icing is presented in Fig.3.In the case of no icing,the wind speed is almost directly proportional to power,whereas,in the case of icing,the relation between the power and the wind speed is lost just as shown in the lower picture in Fig.3.

    Figure 3: (Continued)

    Figure 3:Comparison between cases with and without the icing of the blades of the turbine

    Fig.2 shows the relationship between certain variables and the middle feature in time,and Fig.3 presents a direct comparison between cases with and without icing.They can be combined into one single figure.

    Similarly,we can obtain the relationships between power,and the pitch angle(Fig.4),as well as the ambient temperature(Fig.5),and the difference in temperature(Fig.6).

    Figure 4:Relation between the pitch angle and power

    Fig.4 shows that when the pitch angle was greater than 0.5,the likelihood of icing of the blades was very low.

    The icing of the blades mostly occurs at low temperatures,as shown in Fig.5.

    Figure 5:Relation between ambient temperature and power

    Fig.6 shows that the icing of the blades often occurs in a certain range of temperature,leading to reduced power output.

    Figure 6:Relationship between temperature difference and power

    4 Constructing the DNN Model

    The relationship between the icing of the blades and the feature variables in SCADA is non-linear,and the DNN can be used to fit this relationship in theory [15].The model of icing the blades of a turbine based on the DNN is shown in Fig.7.

    Figure 7:Model of icing of the blades based on the DNN

    The model has five hidden layers and uses ReLU as the activation function cross-entropy as the cost function[16].The sigmoid function is used in the output layer for logistic regression.

    To avoid overfitting when training the model on the data samples and enhance its capability of being used on other datasets,the dropout parameter can be added to the hidden layer [17,18].The values of dropout were[0,0.4,0.3,0.2,0.2].

    5 Experimental Process and Results

    5.1 Indicators of Model Performance

    Predictive accuracy is the most important criterion to assess the performance of a model in general,but is insufficient for the problem of binary classification.Therefore,Li et al.[11]proposed a reward function to measure how a model is fitted.It is shown in Eq.(2):

    where P1 is the predictive accuracy of data with label 1 and P0 is that of data with label 0.Each occupied 50%of the total weight.The weights here are changeable according to the given circumstance.

    The reward function can help avoid the influence of unbalanced classification distribution in the sample data.Its value together with the predictive accuracy determines how well the model fits the data.

    5.2 Experimental Process

    The experimental process is presented in Fig.8.First,the sample data was prepared.After extracting data from the dataset for wind turbine #1 by using the transferred features,they were preprocessed as the sample data.Second,we divided the sample data into training and testing data.The training data was used to train the icing model and the testing data to validate it.Third,we identified the icing model that yielded the highest value of the reward function within the range of the maximum number of epochs.Finally,we used it to predict the data on wind turbine#2,and calculated the values of the reward function as well as predictive accuracy.

    Figure 8:Experiment process of the proposed model

    5.3 Experimental Results and Analysis

    5.3.1 Effects of Different Models

    The data on wind turbines is called time-series data,and the status of ice buildup upon the wind turbine blades is binary.So,the Bi-directional long short-term memory (Bi-LSTM) appears to be a reasonable choice to predict wind fan icing.Xiong et al.built a model of icing the blades of a wind turbine based on Bi-LSTM and the SVM,and reported accurate predictions[9].However,the effect when the model is used to predict the data for a different turbine is unknown.Therefore,we tested different models,including the neural network (NN),DNN,and LSTM,by using the above experimental process.The results are shown in Table 4.

    Table 4: Main indicators of the results of predictions of different models using transferred variables

    Table 4 indicates that when only predictive accuracy was considered,LSTM delivered the best performance.But the reward value is lower when LSTM is transferred to a#2 wind fan.It means that there are many blade-icing alarms missing when the model is transferred to another wind fan,which is not acceptable in practical environments.So,DNN is actually the best one.

    5.3.2 Effect of Transferred Features

    Following the above experimental process,we carried out a contrast test of the 26 transferred features,and the results are presented in Table 5.

    Table 5:Main indicators of the results of prediction obtained by using the DNN model trained on the transferred variables

    We can draw the flowing conclusions from the results in Table 5:1)The results obtained when the transferred features were used were better than those when all 26 features were used,both in terms of accuracy and the reward function.2)The accuracy obtained when using wind turbine#1 data was lower than that when using wind turbine #2 data.But Table 4 shows that the value of the reward function for data from wind turbine#1 was higher than that for turbine#2.This is because the target of the test was to obtain a high reward in addition to high accuracy.

    Fig.9 uses a column chart to directly compare the results of wind turbine#2 and the results of the prediction of the DNN model using the transferred features.

    Figure 9:Comparison charts between the results of wind turbine#2 and the predictions of the DNN model trained on transferred features

    5.3.3 Effect of Further Optimization

    The above DNN model still has room for improvement.Some engineers have suggested that the pitch angle must be changed when wind turbines are covered by ice [6].The variables of the pitch angle(pitch1_angle,pitch2_angle,pitch3_angle)should be combined according to this.For example,the relation between pitch1_angle and power are shown in Fig.10.

    Fig.10 shows that the icing of the blades occurred only when pitch1_angle was small,and the power output was small in this case.This means that pitch1_angle was closely correlated with the icing of the blades.

    A better DNN model is obtained after combining the variables of the pitch angle and using the above experimental process,as shown in Fig.8,and the results are shown in Table 6.

    The DNN model could thus be optimized,which shows that the use of suitable feature variables related to the icing of the blades can improve the model.Admittedly,there is room to improve the current model.

    Figure 10:Relation between pitch1_angle and power of the wind turbine

    Table 6: Main indicators of results of prediction obtained by using the optimized DNN model

    5.3.4 Comparison with Recent Work

    In 2020,Xu et al.[7] proposed a model of blade-ice development by using particle swarm optimization and the support vector machine.They obtained reasonably accurate predictions as shown in Table 7.

    Table 7: Predictions of icing of the blades in[7]

    A comparison of Tables 6 and 7 shows that the model of the blades icing proposed in[7]had higher accuracy(nearly 97%)on the validation dataset than the optimized DNN model(91.58%)proposed in this paper.But an interesting result was obtained on the test dataset:the accuracy of icing free in Table 7 (94.5%) is higher than one in Table 6 (90.03%),on the other hand,the accuracy of icing in Table 7(80.8%)is lower than one in Table 6(88.69%).So,it is difficult to decide which model is the better one only by accuracy indicator.However,when the reward indicator was used,we found that its value in Table 6(89.36)was higher than that in Table 7(87.65).This indicates that the proposed model obtained more balanced predictions than the one proposed in[7].Furthermore,it highlights the need for an integrated system of indicators.

    6 Conclusions

    Models obtained by using machine learning can help improve products in a variety of industries.In this study,the authors transferred some feature variables related to the icing of the blades of a wind turbine from previous studies to obtain a model that can accurately predict icing.To establish a correlation between a given feature and the status of icing,power was used as a middle feature.This was used to establish a correlation between a continuous variable and the binary status of icing on the blades.The results of experiments showed that the transferred feature variables can obtain better results than the original features.The conclusions of this study are as follows:

    1.An optimized model of blade icing upon turbines,based on DNN,was developed that delivered accurate predictions when transferred to data on other turbines.Once the value of the reward function was considered,the DNN was more suitable for such transference of data than LSTM,NN,and logistic regression.

    2.It is suitable to build an integrated system of indicators,in addition to simply those regarding accuracy according to the given problem because this can help measure the model along different dimensions.Although the predictive accuracy of the machine learning model is an important indicator of its performance,it is not adequate on its own.The best model maybe not have the highest accuracy.

    3.The feature variables,in addition to the training data,directly affected the predictions of the model.The predictions of the model were improved by transferring feature variables related to the icing of the blades from past studies to it,which means that knowledge of the domain can help quickly obtain a better model.Furthermore,the use of the middle feature shows that it can serve as a bridge between feature variables.Such indirect analysis can provide a solution to the problem considered here.

    4.A better model can be obtained by using features that are more closely related to the icing of the blades than using original ones,and a better indicators system can help to determine which algorithm is more feasible than using only the metric of accuracy.

    Credit Authorship Contribution Statement:Tingshun Li&Jiaohui Xu:Writing–review&editing.Wen Tan:Conceptualization.Zesan Liu&Dadi Wang:Project administration.

    Funding Statement:This work was supported by the National Natural Science Foundation of China under Grant No.61573138.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    日韩精品青青久久久久久| av天堂在线播放| 精品卡一卡二卡四卡免费| 中文字幕人妻熟女乱码| 亚洲精品粉嫩美女一区| 日韩欧美 国产精品| 可以在线观看毛片的网站| 久9热在线精品视频| 国产亚洲欧美在线一区二区| 国产av在哪里看| 欧美精品啪啪一区二区三区| 午夜福利成人在线免费观看| cao死你这个sao货| 婷婷六月久久综合丁香| 亚洲中文字幕日韩| 精品国内亚洲2022精品成人| bbb黄色大片| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片免费观看直播| 老熟妇乱子伦视频在线观看| 最新美女视频免费是黄的| 欧美色视频一区免费| 露出奶头的视频| 国产不卡一卡二| 国产亚洲精品久久久久久毛片| 啦啦啦韩国在线观看视频| 欧美色欧美亚洲另类二区| 亚洲欧美一区二区三区黑人| 国产精品久久久久久精品电影 | 亚洲第一青青草原| 亚洲精品中文字幕在线视频| 在线国产一区二区在线| 老司机午夜福利在线观看视频| 国产精品 国内视频| 久久久水蜜桃国产精品网| 国产野战对白在线观看| 欧美精品亚洲一区二区| 午夜免费鲁丝| 国产区一区二久久| 亚洲精华国产精华精| 免费高清视频大片| 午夜福利18| 黄色成人免费大全| 久久青草综合色| 精品一区二区三区av网在线观看| 国产精品美女特级片免费视频播放器 | 听说在线观看完整版免费高清| 巨乳人妻的诱惑在线观看| 国产亚洲精品久久久久久毛片| 中文字幕高清在线视频| 成人三级做爰电影| 99国产综合亚洲精品| 一a级毛片在线观看| 99久久久亚洲精品蜜臀av| 香蕉久久夜色| 国产又黄又爽又无遮挡在线| 99热6这里只有精品| 精品欧美一区二区三区在线| 国产av又大| 国产伦人伦偷精品视频| 欧美午夜高清在线| 日韩精品青青久久久久久| 国语自产精品视频在线第100页| 成人一区二区视频在线观看| 国产v大片淫在线免费观看| 国产亚洲av嫩草精品影院| 国内毛片毛片毛片毛片毛片| 午夜免费观看网址| 午夜福利18| 国产一区二区三区视频了| 免费在线观看视频国产中文字幕亚洲| 两个人看的免费小视频| 最新在线观看一区二区三区| 午夜免费激情av| 精品国产乱码久久久久久男人| 在线永久观看黄色视频| 久久久久亚洲av毛片大全| 精品少妇一区二区三区视频日本电影| 色播在线永久视频| 欧美大码av| 国产v大片淫在线免费观看| 99在线视频只有这里精品首页| 黄频高清免费视频| 亚洲熟妇中文字幕五十中出| 巨乳人妻的诱惑在线观看| 午夜视频精品福利| 熟女电影av网| 好看av亚洲va欧美ⅴa在| 日日爽夜夜爽网站| 老司机午夜福利在线观看视频| 一本一本综合久久| 久久久久久久精品吃奶| 成人免费观看视频高清| 国产精品香港三级国产av潘金莲| 午夜福利高清视频| 一级黄色大片毛片| 中文资源天堂在线| 国产精品一区二区免费欧美| 久久久久国产一级毛片高清牌| 久久久国产成人免费| 中亚洲国语对白在线视频| 女性生殖器流出的白浆| 日韩欧美一区二区三区在线观看| av天堂在线播放| 国产爱豆传媒在线观看 | 韩国av一区二区三区四区| 黄色a级毛片大全视频| 老司机福利观看| xxxwww97欧美| 久9热在线精品视频| 麻豆av在线久日| 国产亚洲精品久久久久5区| 哪里可以看免费的av片| 香蕉av资源在线| 国产成人系列免费观看| 久久久久久人人人人人| 熟女电影av网| 老熟妇乱子伦视频在线观看| 国产熟女xx| 亚洲精品av麻豆狂野| 亚洲第一电影网av| www国产在线视频色| 免费在线观看黄色视频的| a在线观看视频网站| 最近最新中文字幕大全免费视频| 成人亚洲精品一区在线观看| 亚洲av成人av| 欧美黄色片欧美黄色片| 黄色女人牲交| 啦啦啦免费观看视频1| 精品人妻1区二区| 久久久久久九九精品二区国产 | 激情在线观看视频在线高清| 久久久久免费精品人妻一区二区 | 在线视频色国产色| 日韩欧美国产一区二区入口| 精品国产超薄肉色丝袜足j| 婷婷精品国产亚洲av| 欧美日韩瑟瑟在线播放| 搡老熟女国产l中国老女人| 99国产精品99久久久久| 久久伊人香网站| 最新在线观看一区二区三区| 美女 人体艺术 gogo| 麻豆国产av国片精品| 十八禁人妻一区二区| 一级片免费观看大全| 久久久久久久精品吃奶| 日本a在线网址| 欧美日韩乱码在线| a级毛片在线看网站| 搡老妇女老女人老熟妇| 嫩草影院精品99| 亚洲av电影在线进入| 国产午夜精品久久久久久| 岛国视频午夜一区免费看| 丰满的人妻完整版| 婷婷六月久久综合丁香| 日本 av在线| 午夜福利欧美成人| 国产又爽黄色视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产三级在线视频| 男人舔奶头视频| 亚洲 国产 在线| 成人国产一区最新在线观看| 久久久久久九九精品二区国产 | 久久人妻av系列| 亚洲欧美激情综合另类| 色精品久久人妻99蜜桃| 一个人观看的视频www高清免费观看 | 国产亚洲欧美在线一区二区| 国产午夜精品久久久久久| 久久香蕉国产精品| 啦啦啦 在线观看视频| 国产黄a三级三级三级人| 欧美激情极品国产一区二区三区| 深夜精品福利| 长腿黑丝高跟| 国产av不卡久久| www.999成人在线观看| 日日干狠狠操夜夜爽| 欧美成人性av电影在线观看| 狠狠狠狠99中文字幕| 99国产精品一区二区蜜桃av| 51午夜福利影视在线观看| 男女之事视频高清在线观看| а√天堂www在线а√下载| cao死你这个sao货| 久久天躁狠狠躁夜夜2o2o| 又大又爽又粗| 黄片大片在线免费观看| 在线观看免费日韩欧美大片| 亚洲色图 男人天堂 中文字幕| 一级作爱视频免费观看| 午夜精品久久久久久毛片777| 日韩有码中文字幕| 午夜激情福利司机影院| 亚洲一区高清亚洲精品| 国产97色在线日韩免费| 亚洲欧美精品综合久久99| 熟女少妇亚洲综合色aaa.| 高清在线国产一区| 国产精品一区二区免费欧美| 日韩大尺度精品在线看网址| 一级毛片精品| 亚洲 欧美 日韩 在线 免费| 国产精品影院久久| 99在线人妻在线中文字幕| 极品教师在线免费播放| 91国产中文字幕| 一级作爱视频免费观看| 中文字幕精品亚洲无线码一区 | 啦啦啦观看免费观看视频高清| 黄色毛片三级朝国网站| 91九色精品人成在线观看| 色婷婷久久久亚洲欧美| 在线观看66精品国产| 久久久久久免费高清国产稀缺| ponron亚洲| 男女做爰动态图高潮gif福利片| 午夜a级毛片| 黄片小视频在线播放| 19禁男女啪啪无遮挡网站| 精品久久久久久久毛片微露脸| 黄频高清免费视频| 中亚洲国语对白在线视频| 欧美成狂野欧美在线观看| 午夜精品在线福利| bbb黄色大片| 欧美一级毛片孕妇| 中文字幕最新亚洲高清| 中文字幕av电影在线播放| 一本大道久久a久久精品| 国产精品一区二区精品视频观看| 亚洲 欧美 日韩 在线 免费| 亚洲全国av大片| 亚洲激情在线av| 午夜精品在线福利| 人人妻人人澡人人看| 午夜福利成人在线免费观看| 久久天堂一区二区三区四区| 他把我摸到了高潮在线观看| 亚洲国产欧洲综合997久久, | 黑人欧美特级aaaaaa片| 成人精品一区二区免费| 99国产精品一区二区三区| 久久香蕉精品热| 韩国精品一区二区三区| 久久国产精品男人的天堂亚洲| 18禁国产床啪视频网站| 国产精品爽爽va在线观看网站 | 很黄的视频免费| www.精华液| 可以在线观看毛片的网站| 天天一区二区日本电影三级| 亚洲 欧美 日韩 在线 免费| 国产成人系列免费观看| 91av网站免费观看| 很黄的视频免费| 中文资源天堂在线| 神马国产精品三级电影在线观看 | 亚洲av日韩精品久久久久久密| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看 | 国产精品野战在线观看| 亚洲国产精品合色在线| 亚洲精品在线美女| 人人妻人人澡欧美一区二区| 精品国产美女av久久久久小说| 欧美激情极品国产一区二区三区| 日日干狠狠操夜夜爽| 国产欧美日韩一区二区精品| 免费看日本二区| 国产精品一区二区免费欧美| 国产精品亚洲美女久久久| 欧美三级亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品二区激情视频| 狠狠狠狠99中文字幕| 欧美色视频一区免费| 色av中文字幕| 欧美日本视频| 中文亚洲av片在线观看爽| 欧美又色又爽又黄视频| 国产av一区二区精品久久| or卡值多少钱| 大香蕉久久成人网| 久久久久国产一级毛片高清牌| 中文字幕久久专区| 激情在线观看视频在线高清| 日本免费一区二区三区高清不卡| 久久精品aⅴ一区二区三区四区| 日韩三级视频一区二区三区| 久久人妻福利社区极品人妻图片| 欧美成人一区二区免费高清观看 | 看免费av毛片| 欧美 亚洲 国产 日韩一| 一级黄色大片毛片| 男人舔女人的私密视频| 美女午夜性视频免费| 亚洲av成人av| 精品免费久久久久久久清纯| 亚洲成av人片免费观看| 成年免费大片在线观看| 亚洲欧洲精品一区二区精品久久久| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影| 国产爱豆传媒在线观看 | 国产免费男女视频| 亚洲电影在线观看av| 少妇被粗大的猛进出69影院| 琪琪午夜伦伦电影理论片6080| 亚洲一码二码三码区别大吗| 精品熟女少妇八av免费久了| 日本三级黄在线观看| 久久精品国产综合久久久| 久久人人精品亚洲av| 亚洲久久久国产精品| 亚洲aⅴ乱码一区二区在线播放 | 麻豆成人午夜福利视频| 国内久久婷婷六月综合欲色啪| 香蕉久久夜色| 999久久久国产精品视频| 免费在线观看完整版高清| 国产成人精品无人区| 日韩一卡2卡3卡4卡2021年| 久久久久久久久免费视频了| 成年版毛片免费区| 国产野战对白在线观看| 欧美激情极品国产一区二区三区| 这个男人来自地球电影免费观看| 99久久综合精品五月天人人| 午夜福利成人在线免费观看| 中文资源天堂在线| 国产av在哪里看| 国产精品日韩av在线免费观看| 大香蕉久久成人网| √禁漫天堂资源中文www| 熟妇人妻久久中文字幕3abv| 大型黄色视频在线免费观看| 男女床上黄色一级片免费看| 男女视频在线观看网站免费 | 非洲黑人性xxxx精品又粗又长| 黄色a级毛片大全视频| 久久国产精品影院| 中文字幕人成人乱码亚洲影| 亚洲精品中文字幕一二三四区| 女警被强在线播放| 人人澡人人妻人| 中文亚洲av片在线观看爽| 草草在线视频免费看| 午夜亚洲福利在线播放| 国产日本99.免费观看| 国产精品野战在线观看| 午夜福利高清视频| 日韩三级视频一区二区三区| 桃红色精品国产亚洲av| 亚洲精品国产一区二区精华液| 日本成人三级电影网站| 免费高清视频大片| 欧美成人免费av一区二区三区| 日韩免费av在线播放| 男女视频在线观看网站免费 | 日本熟妇午夜| 变态另类丝袜制服| 夜夜夜夜夜久久久久| 岛国在线观看网站| 亚洲av美国av| 99久久精品国产亚洲精品| 99国产精品99久久久久| 精品国产美女av久久久久小说| 国产亚洲欧美98| 欧美三级亚洲精品| 久久久久久久精品吃奶| 久久 成人 亚洲| 日韩欧美一区视频在线观看| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 极品教师在线免费播放| 国产人伦9x9x在线观看| 黄色成人免费大全| 免费在线观看完整版高清| 国产成人av教育| 成人国产综合亚洲| 国产精品,欧美在线| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 日韩视频一区二区在线观看| 久久久久久久久免费视频了| 日韩精品免费视频一区二区三区| 免费在线观看亚洲国产| 亚洲 欧美 日韩 在线 免费| 日韩av在线大香蕉| 国产伦在线观看视频一区| 天堂动漫精品| 国产又黄又爽又无遮挡在线| 他把我摸到了高潮在线观看| e午夜精品久久久久久久| 婷婷精品国产亚洲av| 免费看美女性在线毛片视频| 精品国内亚洲2022精品成人| 亚洲国产毛片av蜜桃av| 99在线视频只有这里精品首页| 婷婷精品国产亚洲av| 国产一区二区三区在线臀色熟女| 国产精品一区二区免费欧美| 久久亚洲精品不卡| 久久久久久亚洲精品国产蜜桃av| 亚洲国产高清在线一区二区三 | 亚洲人成电影免费在线| 国产精品野战在线观看| 男人舔女人下体高潮全视频| 欧美日韩福利视频一区二区| 成年版毛片免费区| 黄色视频不卡| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 欧美zozozo另类| 99久久99久久久精品蜜桃| 亚洲狠狠婷婷综合久久图片| 少妇熟女aⅴ在线视频| 欧美黑人巨大hd| 亚洲中文字幕日韩| 国产国语露脸激情在线看| 亚洲av成人不卡在线观看播放网| www国产在线视频色| 一区福利在线观看| 久久久久久久久免费视频了| 欧美性猛交黑人性爽| 成人永久免费在线观看视频| a级毛片在线看网站| 日韩中文字幕欧美一区二区| 国产麻豆成人av免费视频| 久久人人精品亚洲av| 身体一侧抽搐| 欧美人与性动交α欧美精品济南到| 一进一出抽搐gif免费好疼| 三级毛片av免费| 波多野结衣av一区二区av| 熟女少妇亚洲综合色aaa.| 亚洲激情在线av| 一级毛片女人18水好多| 亚洲九九香蕉| 国产精品免费视频内射| 好男人在线观看高清免费视频 | 亚洲男人的天堂狠狠| 亚洲成人国产一区在线观看| 亚洲av中文字字幕乱码综合 | 熟女电影av网| 99久久国产精品久久久| 中文字幕人成人乱码亚洲影| 欧美久久黑人一区二区| 一个人免费在线观看的高清视频| 欧美性长视频在线观看| 久久精品成人免费网站| 亚洲全国av大片| 久久精品国产99精品国产亚洲性色| 久久精品夜夜夜夜夜久久蜜豆 | 国产伦在线观看视频一区| 亚洲国产精品成人综合色| 亚洲国产高清在线一区二区三 | 欧美日韩黄片免| 美女 人体艺术 gogo| 亚洲第一av免费看| 精品第一国产精品| 嫩草影院精品99| 美女国产高潮福利片在线看| tocl精华| 午夜福利18| 在线播放国产精品三级| 亚洲精品粉嫩美女一区| 色婷婷久久久亚洲欧美| 午夜精品在线福利| 亚洲美女黄片视频| 国产人伦9x9x在线观看| 午夜福利免费观看在线| 老司机在亚洲福利影院| 久久久国产成人免费| 91成人精品电影| 成人18禁在线播放| 久久久久久国产a免费观看| 超碰成人久久| 日韩大码丰满熟妇| 韩国精品一区二区三区| www.精华液| 国产av不卡久久| 熟妇人妻久久中文字幕3abv| 一区二区三区精品91| 国产精品久久久久久亚洲av鲁大| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉国产精品| 成人国产综合亚洲| 麻豆成人午夜福利视频| 精华霜和精华液先用哪个| 男女那种视频在线观看| 叶爱在线成人免费视频播放| 国产精品亚洲一级av第二区| 精品无人区乱码1区二区| www日本黄色视频网| xxx96com| 成在线人永久免费视频| 在线观看免费视频日本深夜| 无人区码免费观看不卡| 婷婷六月久久综合丁香| 成人精品一区二区免费| 国产精品 国内视频| 欧美一区二区精品小视频在线| 国产v大片淫在线免费观看| 免费在线观看完整版高清| 亚洲国产日韩欧美精品在线观看 | 好男人电影高清在线观看| 精品久久蜜臀av无| 1024香蕉在线观看| 亚洲av成人一区二区三| 久久精品成人免费网站| 久久久久亚洲av毛片大全| 欧洲精品卡2卡3卡4卡5卡区| 国产三级黄色录像| 成年版毛片免费区| 日本一区二区免费在线视频| 男人舔女人下体高潮全视频| 成人欧美大片| 成年女人毛片免费观看观看9| 香蕉av资源在线| 亚洲av成人一区二区三| 一边摸一边抽搐一进一小说| 午夜福利在线在线| 午夜久久久久精精品| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 亚洲成人久久爱视频| 久久精品夜夜夜夜夜久久蜜豆 | 岛国在线观看网站| 别揉我奶头~嗯~啊~动态视频| 美女免费视频网站| 国产99久久九九免费精品| 亚洲精品中文字幕一二三四区| 亚洲一区二区三区色噜噜| 成年女人毛片免费观看观看9| 亚洲久久久国产精品| 免费在线观看视频国产中文字幕亚洲| tocl精华| av在线播放免费不卡| 国产在线精品亚洲第一网站| 日本免费a在线| 亚洲精品在线美女| 精品欧美一区二区三区在线| 久久久水蜜桃国产精品网| 色老头精品视频在线观看| 精品熟女少妇八av免费久了| 久久精品aⅴ一区二区三区四区| 一级毛片高清免费大全| av超薄肉色丝袜交足视频| 神马国产精品三级电影在线观看 | 在线观看www视频免费| 首页视频小说图片口味搜索| 亚洲成av片中文字幕在线观看| 99国产精品一区二区三区| 麻豆成人午夜福利视频| 欧美黑人精品巨大| netflix在线观看网站| 午夜免费激情av| 亚洲国产毛片av蜜桃av| 欧美最黄视频在线播放免费| 男女视频在线观看网站免费 | 亚洲九九香蕉| 国产一区二区三区视频了| 成熟少妇高潮喷水视频| 亚洲无线在线观看| 午夜两性在线视频| 老司机午夜福利在线观看视频| 哪里可以看免费的av片| 国产精品一区二区精品视频观看| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| 欧美日韩瑟瑟在线播放| 热99re8久久精品国产| 国产精品,欧美在线| 国产精品乱码一区二三区的特点| 好男人在线观看高清免费视频 | 可以在线观看毛片的网站| 日本 欧美在线| 国产三级在线视频| 一级毛片女人18水好多| 日本 欧美在线| 午夜影院日韩av| 怎么达到女性高潮| 亚洲国产欧美一区二区综合| 免费看a级黄色片| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 久久欧美精品欧美久久欧美| 国产aⅴ精品一区二区三区波| 亚洲欧美一区二区三区黑人| 亚洲 国产 在线| 观看免费一级毛片| 成人午夜高清在线视频 | 国产久久久一区二区三区| 99国产综合亚洲精品| 老司机靠b影院| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美免费精品| 国产精品综合久久久久久久免费| 正在播放国产对白刺激| 51午夜福利影视在线观看| 久久精品人妻少妇| 中文在线观看免费www的网站 | 亚洲男人的天堂狠狠| 亚洲性夜色夜夜综合| 丝袜在线中文字幕| 国产av一区在线观看免费| 51午夜福利影视在线观看| 韩国精品一区二区三区| 国产成人av激情在线播放|