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

    Traffic flow prediction based on BILSTM model and data denoising scheme

    2022-04-12 03:45:36ZhongYuLi李中昱HongXiaGe葛紅霞andRongJunCheng程榮軍
    Chinese Physics B 2022年4期
    關(guān)鍵詞:榮軍紅霞

    Zhong-Yu Li(李中昱) Hong-Xia Ge(葛紅霞) and Rong-Jun Cheng(程榮軍)

    1Faculty of Maritime and Transportation,Ningbo University,Ningbo 315211,China

    2Jiangsu Provincial Collaborative Innovation Center for Modern Urban Traffic Technologies,Nanjing 210096,China

    3National Traffic Management Engineering and Technology Research Center Ningbo University Subcenter,Ningbo 315211,China

    Keywords: traffic flow prediction,bidirectional long short-term memory network,data denoising

    1. Introduction

    Traffic flow prediction is an important component of intelligent transportation system, which can reduce the number of accidents,improve traffic efficiency,and reduce traffic pollution. The traffic prediction is based on a large number of historical data to predict the future traffic flow.[1]Previous studies show that linear model, nonlinear model, and hybrid model are three typical traffic flow prediction techniques.[2]In the linear models the mathematical methods are used to complete the task of the traffic flow prediction. The autoregressive integrated moving average model (ARIMA) is the most commonly used linear model.[3]Chenet al.[4]proposed an ARIMA model. The ARIMA model was specified by different training data, which represent the different traffic states,according to the different periods in different days. The model from different periods would be used on different data when testing. This would make the model more special and more accurate. The data-driven approach using the ARIMA model in most of studies required sound database for building a model. The seasonal ARIMA method was proposed to predict the short-term traffic flow with limited data.[5,6]Although the ARIMA had good prediction accuracy,there were many nonlinear factors in traffic flow prediction. Therefore,the nonlinear models introduced relevant machine learning methods(including support vector regression(SVR),artificial neural networks(ANN),etc.)to achieve favorable traffic flow prediction performance.The SVR had strict theoretical and mathematical basis. Based on the principle of structural risk minimization,it has strong generalization capability and global optimality. It was a theory for small sample statistics.[7-10]Neural network had strong non-linear fitting capability and could learn rules from large samples. Moreet al.[11]predicted the future traffic volume through the ANN,and used the predicted results to control congestion. Kumaret al.[12]used the ANN to make short-term prediction of traffic volume by using past traffic data (taking the traffic volume, speed and density as inputs).The ANN can produce good results in this study. In order to better capture the spatiotemporal information of traffic flow,convolutional neural network has been proposed.[13-15]In order to better learn the historical information of data,recurrent neural network(RNN)was proposed to be able to transfer the output of the previous time to the next step as input, thereby improving the accuracy of prediction.[16]The RNN can handle certain short-term dependence, but it cannot handle the longterm dependence. To solve this problem,the LSTM was proposed.The LSTM introduced cell state,and the input gate,forgetting gate and output gate were used to maintain and control the information,which can learn long-term information.[17]It is noted that the nonlinear models have achieved great success in dealing with many traffic flow prediction tasks. However,the uncertainty of traffic flow data might have a negative effect on the nonlinear prediction model, which would reduce the performance of the model. The hybrid models were proposed to overcome the disadvantages.[18,19]For example,Tanget al.[20]combined fuzzy c-means and genetic algorithm to predict missing traffic volume data. Zhenget al.[21]used a hybrid model based on a convolutional neural network(CNN)and the LSTM.The experimental results show that the prediction accuracy of the hybrid model is higher than that of the single model.

    Traffic flow data can be obtained by many technical means, such as remote microwave, loop deductive detectors,global positioning system(GPS).[22]Historical data with high fidelity play a critical role in predicting results, but raw data always contain noise. Therefore, one introduced many noise reduction methods,[23]such as the wavelet Kalman filter model[24,25]and wavelet transform.[26,27]Yanget al.[28]selected the appropriate wavelet basis to decompose and reconstruct the traffic flow data. The results showed that prediction effect is better than the results predicted just by neural network in prediction precision and network convergence. Because the wavelet transform needed to select the appropriate wavelet base, sometimes it is difficult to find the appropriate wavelet base. One proposed the EMD,which is an adaptive decomposition method based on the original data.[29]The EMD extracts the intrinsic mode function(IMF)sets from the input data samples, which can be separated into high-frequency (HF) part and low-frequency (LF) part. The EMD is affected by mode mixing,so the EEMD was proposed.[30,31]Chenet al.[32]proposed an integrated framework based on the integrated EEMD and the ANN can predict the traffic flow at different time intervals ahead. The experimental results show that the noise reduction algorithm greatly improves the accuracy of traffic flow prediction.

    The above method combines a certain noise reduction method with the prediction method to achieve better prediction accuracy, but does not use the same set of data to verify the performance of different noise reduction algorithms. The present study intends to combine data denoising algorithms with deep learning model to predict the traffic flow and compare among prediction accuracies of various methods. The contributions of this paper are as follows:

    (i)We use denoising algorithms to process the raw data,which include WL,EMD,and EEMD,and compare the capabilities of these algorithms to suppress the outliers.

    (ii) We build a BILSTM prediction model to predict the data already processed.

    (iii)We consider the predicted effects of the model under different road scenarios(i.e.,mainline,on ramp,off ramp)and on different time scales(i.e.,5 min,10 min,and 15 min).

    The rest of this paper is organized as follows. In Section 2,the denoising methods and BILSTM are introduced. In Section 3,the source of data,the model structure,and prediction results are presented. In Sections 4 some conclusions are drawn from the present study.

    2. Method

    The methods to be used in this article are data denoising and deep learning neural networks to improve prediction accuracy. The frame of this model is shown in Fig.1.

    Fig.1. Framework of prediction model.

    2.1. Denoising methods

    2.1.1. Wavelet

    Wavelet is used to eliminate noise interference in traffic flow data.[33]The purpose of building a noise reduction model is to reduce the influence of interference factors,such as gaussian noise,on the prediction accuracy of the model in a complex,nonlinear traffic data set. One-dimensional noise model can be expressed as

    whereFtis the raw signal with noise;ytis the signal without noise;εtis Gaussian noise;nis the signal length. The principle of wavelet denoising is to suppress theεtpart ofFtand restore theytpart ofFt. Wavelet denoising model preserves the original signal characteristics after denoising. The flow chart of wavelet is shown in Fig.2. The wavelet algorithm is implemented in three steps.

    Step 1 Decomposition

    Choose a wavelet and determine the layerNfor a wavelet decomposition,then compute theN-layer wavelet decomposition for the signal.

    Step 2 Setting a threshold for high-frequency coefficients

    For the high-frequency factor of each layer from layer 1 to layerN,a threshold is selected for thresholding.

    For most of signals,the low-frequency component is very important,and it usually contains the characteristics of the signal, while the high-frequency component contains the details or differences of the signal. The low-frequency information can approximately describe the change trend of the original signal, while high-frequency information represents the detailed information of the signal. Therefore,the original signal can produce two signals (high- and low-frequencies) through two mutual filters.

    Step 3 Data reconstruction

    According to the approximate coefficients and the modified wavelet coefficients,the inverse wavelet transform is used to reconstruct the traffic data after denoising.

    Fig.2. Flowchart of wavelet.

    2.1.2. The empirical mode decomposition (EMD)

    In this paper, we use the EMD, which was first introduced by Huanget al.[34]Different applications in medical signal analysis have shown the effectiveness of this method.The principle is to adaptively decompose the given signal into frequency components,which called intrinsic mode functions(IMF).These components are obtained from the signal by an algorithm. The algorithm extracts the highest frequency oscillation of each mode from the original signal. The more detailed explanations to the EMD model are described as follows.

    (i)Decompose raw traffic flow dataR(t)into IMFs. Find out potential extreme points fromR(t) data series, then employ the cubic spline interpolation method to connect the maximum and minimum points to form the upper envelopeU(t)and lower envelopeL(t).

    (ii) Calculate the mean of the upper envelopeU(t) and lower envelopeL(t) to obtain the mean envelopeM(t) (see Eq.(2)).

    (iii) Calculate the difference betweenR(t) andM(t) to obtainI(t)(see Eq.(3)). IfI(t)satisfies the IMF’s conditions,I(t) is one of the IMFs, otherwise the above three sub-steps need repeating untilI(t) satisfies IMF’s condition. The following IMF’s conditions are satisfied:i)the number difference between extrema and the zero-crossing points is not larger than one;ii)at any point in the IMF,the mean value of the envelope defined by the local maximum and the mean value of the envelope defined by the local minimum should be equal to zero.

    (iv)Obtain the residual partr(t)by calculating the difference betweenR(t)andI(t),repeat the above three sub-steps to obtain new IMFs until the residual data distribution is monotonic or has one extreme point. Finally,R(t) is decomposed into a series of IMFs and a residual partr(t)as shown by the following equations:

    2.1.3. The ensemble empirical mode decomposition(EEMD)

    For some data with too few extreme points, EMD fails to work, so the EEMD method is proposed. The white noise is added into the EEMD to make the data meet the requirements for the EMD. Owing to the zero mean white noise,the zero mean white noise will cancel each other after many times of average calculation, so the calculation result of integrated mean can be directly regarded as the final result. The EEMD can effectively suppress the mode aliasing of EMD.The EEMD algorithm steps are as follows.

    Step 1 The white noisewk(t) of normal distribution is added into the raw traffic dataR(t)(the noise-aided traffic flow data are denoted asRk(t))(see Eq.(5)).

    Step 2 TheRk(t) are decomposed into IMFs and residual components. The data decomposition program will be executed in thek-th round until the conditionk >Enis met(Enis the ensemble number). TheRk(t)data can be decomposed into IMFsek,d(t)and corresponding residual partrk,d(t)(see Eq.(6)).

    Step 3 Repeat steps 1 and 2,add a new normal distribution white noise sequence each time.

    Step 4 Select the noise-free IMFs and residuals to reconstruct the EEMD smoothed traffic flow dataR(t)EEMDwhich are given below:

    2.2. Deep learning

    2.2.1. Long short-term memory (LSTM) network

    The LSTM is an improved neural network based on recurrent neural network (RNN). The LSTM was proposed by Hochreiter and Schmidhuber in 1997.[35]In the RNN the output of the previous moment is taken as the input of the next moment,and the RNN is often used to predict the data of time series. However,the RNN tends to disappear and explode gradients when dealing with long sequences,leading the model to be nonconvergent.To solve this problem,the input gate,forget gate, output gate, and memory cell (C(t)) are added into the LSTM.The LSTM structure is shown in Fig.3.

    Fig.3. LSTM structure.

    The three gates are connected by the Sigmoid activation function. The Sigmoid activation function is a hyperbolic tangent function, and the value field is controlled in an interval of[0,1]. Candidate memory cells(?Ct)are connected by tanh activation function,and the range was controlled in an interval of[-1,1]. The calculation steps of the model are as follows:

    whereI(t),F(t),O(t),and ?C(t)represent the input gate,forget gate, output gate, candidate memory cells, respectively;Wxi,Whi,Wx f,Whf,Wxo,Who,Wxc,andWhcare the weight matrices.bi,bf,bo, andbcare bias vectors.σis a Sigmoid activation function.

    Memory cell (Ct) is obtained fromCtand ?Ct. Then we acquireHtby combining ?Ctwith the tanh activation function as follows:

    2.2.2. Bi-directional long short-term memory network(BILSTM)

    The Long short-term memory is a forward training model,which can only capture the long-term characteristics of forward historical traffic flow data. However,the traffic flow data not only depend on past information, but also correlates with future road conditions. The BILSTM is an optimized LSTM model,which combines forward LSTM and reverse LSTM to obtain the information. Therefore, the BILSTM can process and integrate data forward and backward,which can make the prediction more accurate. The BILSTM structure is illustrated in Fig.4. The forward layer and backward layer are connected to the output layer,which contains six shared weights.

    Fig.4. Structure of BILSTM.

    In the forward layer,the forward calculation is performed from timelto timet, and the output of each forward hidden layer is obtained and saved. In the backward layer, the output of the backward hidden layer at each time is obtained and saved by backward calculation from timetto timel. Finally,the final output is obtained by combining the output results of forward layer and backward layer at each time,the calculation steps of the model are as follows:

    In Fig.4,xtrepresents the input of LSTM,Ofdenotes the output of forward LSTM,andObrefers to the output of reverse LSTM.

    3. Experiments

    3.1. Data source

    The traffic flow data in this paper are obtained from the performance measurement system (PeMS). The PeMS is an intelligent traffic management tool. It collects and stores realtime traffic data from sensors and provides online use. The PeMS data collection covers major metropolitan areas in California. We collected the traffic flow data from three detectors installed in Interstate 15 (I15) freeway segment located in San Diego, California (see in Fig. 5). The three detectors are located in mainline(denoted as detector A),on-ramp(represented as detector B), and off-ramp (marked as detector C)(see in Fig. 6) respectively. Owing to the different road conditions of these three sections, the predictions of these three sections can better illustrate the effectiveness of the model.We selected 2880 data in total from weekdays 1 March 2021 through 12 March 2021, and took the flow every 5 min. The details of the road sections are shown in Table 1.

    Fig.5. Interstate 15(I15)freeway segment.

    Fig.6. Three detectors.

    Table 1. Details of road sections.

    3.2. Experiment setups

    Firstly, the data are subjected to noise reduction. The wavelet denoising is a very common denoising method. It has many wavelet bases. Four commonly used wavelet bases are selected in this paper: WL (db4), WL (coif), WL (sym), and WL (haar). When using wavelet denoising, we set the number of decomposition layers to be 3 and the wavelet bases are assumed to be db4,coif2,haar,and sym2 respectively. In using EMD, we set the number of IMF decompositions to be 5, and the interpolation used is pchip (piecewise cubic Hermite interpolating polynomial method). In using EEMD,there are two parameters (the standard deviation ratio between the added Gaussian white noise and the amplitude of the input signal,An, and the ensemble number,En)need to be set. We setAn=0.2 andEn=1000. We use Adam gradient descent algorithm in BILSTM.The data are used four times to make a prediction once. We divide 70%of the raw traffic data as the training dataset and 30%as the test dataset,use Adam gradient descent algorithm,then train 500 times to obtain the prediction results.

    3.3. Evaluation of predictive performance of model

    The root mean square error (RMSE), the mean absolute error(MAE),and the mean absolute percentage error(MAPE)are used to evaluate the prediction performance of the model.The RMSE,MAE,and MAPE are calculated as follows:

    whereXprediction,tis the prediction results,Xreal,tis the raw traffic data,andNis the data-sampling number.

    3.4. Results

    3.4.1. Traffic flow data denoising for detector A

    We employ the wavelet model, EMD, EEMD to remove outliers from original data. In order to see clearly, the noise reduction results of the first 1000 data are shown in Fig. 7.The reliable data are provided for subsequent prediction,blue lines on the plots represent denoised data, and red lines refer to the raw data. The data of a day is 288,and we can see that these data are periodic. It can be observed that the traffic flow in a day can be divided into peak period and flat peak period.The fluctuation of the data in the peak period is significantly higher than that in the flat peak period. This is because the road occupancy rate is high in the peak period and the driving behaviors of drivers will have a great influence on the road conditions, so the noise of the traffic flow in the peak period will be much higher than that in the flat peak period. It can be seen that the noise interference in the flat hump period is less,and the noise reduction results are almost consistent with the raw data. It can be seen that the curve of data after noise reduction are significantly smoothened than that of the original data,which smoothens the fluctuations from the original data.As shown in Fig.7, WL(coif2)denoised data (see Fig. 7(a))and WL(db4)denoised data(see Fig.7(b))are better than WL(haar)denoised data(see Fig.7(c))and WL(sym2)denoised data(see Fig.7(d))in data denoising. The curve of EEMD denoised data is smoother than that of the EMD denoised data.

    Fig.7. Denoising results for detector A.

    3.4.2. Predictions for detector A

    Now, we come to acquire the smooth data and use the BILSTM model to make predictions for detector A.The traffic flow prediction schemes are denoted as follows: BILSTM,BILSTM+WL (coif2), BILSTM+WL (db4), BILSTM+WL(haar), BILSTM+WL (sym2), BILSTM+EMD, and BILSTM+EEMD. The predicted results are shown in Fig. 8.The black line represents the raw data and the red line denotes the predicted data. The results of RMSE, MAE, and MAPE are shown in Table 2.The RMSE,MAE,and MAPE measured directly from BILSTM are 30.68, 22.70, 8.76%, respectively(see Table 2). The predicted RMSE decreases by 20% after noise reduction. It can be seen that the data after noise reduction can be predicted to better learn the change trend of traffic flow and improve the predict accuracy.

    Fig.8. Results of prediction(detector A).

    Table 2. Results of different prediction models for detector A.

    Furthermore, we can see that the EEMD has the best effect. The RMSE,MAE,and MAPE of BILSTM+EEMD(see Table 2) are 3.98, 3.5, and 2.04% respectively. The performances of EEMD and EMD combined with the BILSTM are better than that of wavelet combined with the BILSTM, indicating the excellent performance of EMD with respect to EEMD in eliminating data anomalies.

    3.4.3. Traffic flow data prediction and denoising for ramp detectors B, C

    To further validate the performance of the model,we will take the data from ramp detectors(detector B,detector C)for prediction.

    Fig.9. Denoising results for detector B.

    Fig.10. Denoising results for detector C.

    Since the road situation of the ramp is different from that of the main line, we verify the predictive capability of the model in different road situations. The data denoising results for the two detectors are shown in Figs.9 and 10,respectively.The noise reduction procedure effectively suppresses outliers from the raw traffic data.

    The prediction performance for detectors B and C are shown in Tables 3 and 4. The RMSE, MAE, and MAPE for BILSTM of off ramp (see in Table 3) are 9.43, 6.76,and 24.98%, respectively. The RMSE, MAE and MAPE for BILSTM+WL(coif2)of off ramp(see Table 3)are 1.96,1.39,and 5.50%,respectively. The accuracy of noise reduction prediction from the wavelet method is improved by about 50%.Table 3 indicates that the BILSTM+EMD scheme is better than BILSTM combined with wavelet. In Table 4,we observe that the prediction distributions are the same as in Table 2.Moreover,the prediction accuracy of BILSTM+EEMD model for detector B and detector C are better than those of other hybrid models.The RMSE,MAE,and MAPE are 0.74,0.62,and 1.75 for the traffic flow data collected from detector B,and the statistical values are 0.68, 0.49, and 2.42 collected from detector C.Through the above analysis,it can be concluded that the BILSTM+EEMD model shows the predicted performance consistent with the traffic flow data collected from mainline,on-ramp, and off-ramp areas. This means that the hybrid model is suitable for different road conditions.

    Table 3. Results of different prediction models for detector B.

    Table 4. Results of different prediction models for detector C.

    3.4.4. Comparison among traffic flow data prediction for detectors A, B, and C

    In order to study the influence of noise on different roads,we will add the prediction effect of noise reduction under three road conditions: mainline, on-ramp, and off-ramp. The predicted results are shown in Tables 2-4, so we select the results of EEMD for comparisons. The comparison is shown in Table 5. It can be seen that the accuracy of the model is greatly improved after noise reduction. The prediction results of off ramp and on ramp are similar. After noise reduction,the MAPE of the ramp decreases greatly in comparison with that of the mainline. It can be seen that the noise has a greater influence on the ramp. According to the analysis of traffic flow state in reality,the ramp is only a one-way line,and there are more uncertain factors in the driving of vehicles. The noise in the data is not only increased with the increase of traffic flow,but also related to many factors, such as the weather and the working state of loop detector. The traffic flow on the ramp is much less than that on the mainline. Slight noise will have a great influence on the ramp flow prediction. However, owing to the large vehicle base on the main road,the influence of noise on the main road is not so serious as that on the ramp.

    Table 5. Comparison of traffic flow data prediction among detectors(A,B,C).

    3.4.5. Traffic flow data prediction on different time scales for detector A

    In order to verify the effect of the model on the long-term traffic flow prediction, we use the model to predict the traffic flows in different time spans,i.e., 5 min, 10 min, 15 min.we use the mainline data (detector A) to predict the traffic flow. The results are shown in Table 6. It can be observed that with the increase of time span,RMES,MAE,and MAPE all become larger. This is because with the increase of time span, the data sampling frequency turns smaller, resulting in greater noise in the sampling data,and some characteristic information of traffic flow may be confirmed. Therefore, the error of prediction results will increase. The results show that EEMD+BILSTM has the best prediction effect,the RMSEs of EEMD+BILSTM for 10 min and 15 min are 8.95 and 15.36, respectively. The prediction trend of wavelet denoising and EEMD de-noising are also consistent with that of EEMD+BILSTM. The results show that the model can also achieve good long-term forecasting results,and the model can be used for long-term prediction.

    Table 6. Statistical performance on different time scales for detector A.

    4. Conclusions

    In this study a novel traffic flow prediction framework is introduced by integrating a data denoising scheme with a deep learning model. This study adopt the traffic volume data collected from three detectors in San Diego, California, which are cited from the PeMS. To clean the data, we employ various popular and effective noise reduction models to suppress outliers from the raw traffic flow data. We chose three denoising schemes: EMD, EEMD, and Wavelet with different basis (haar, db4, sym2, coif2). Then combine with BILSTM to make predictions. Finally, we compare the prediction accuracies of different models by analyzing the RMSE, MAE,and MAPE indicators. Several interesting conclusions can be summarized as follows.

    (I) The model combines the data denoising with the model, which is better than the prediction process of without denoising strategy.

    (II) The performances of WL (db4), WL (coif2), WL(sym2) are better than that of WL (haar). In addition, the performance of WL (haar) in long-term prediction is not so satisfactory as other three wavelet methods.

    (III)The BILSTM+EEMD obtains the best performance in comparison with all the WL models and BILSTM+EMD.The RMSE,MAE,and MAPE indicators of BILSTM+EEMD are all at least 25%lower than those of other methods.

    (IV) These hybrid forecasting methods can be used in long-term prediction. The model proposed in this paper can make the prediction more accurate. It can be used to optimize traffic organization strategy,improve traffic efficiency,and reduce energy consumption.

    (V)In future,there will be several directions for more indepth study. Firstly, this paper has no predictive influence of environmental factors on the forecast,which can be taken into account later. Secondly, we can take some traffic indicators into account in traffic flow prediction,such as speed and road occupancy. Finally, we can combine deep learning with traffic flow models (e.g., car-following models) to predict traffic flow.

    Acknowledgements

    Project supported by the Program of Humanities and Social Science of the Education Ministry of China (Grant No. 20YJA630008), the Natural Science Foundation of Zhejiang Province, China (Grant No. LY20G010004), and the K C Wong Magna Fund in Ningbo University,China.

    猜你喜歡
    榮軍紅霞
    Bifurcation analysis of visual angle model with anticipated time and stabilizing driving behavior
    請(qǐng)你幫個(gè)忙
    An extended smart driver model considering electronic throttle angle changes with memory
    Stabilization strategy of a car-following model with multiple time delays of the drivers?
    《烏鴉喝水》中的“想”
    A Study of Combination of English Language Teaching and Context
    大東方(2018年1期)2018-05-30 01:27:23
    讓動(dòng)作“活”起來(lái)
    “大愛(ài)有聲 義工行動(dòng) 關(guān)愛(ài)榮譽(yù)軍人”慰問(wèn)活動(dòng) 走進(jìn)廣東省第一榮軍醫(yī)院
    大社會(huì)(2016年6期)2016-05-04 03:42:20
    “光的直線傳播”“光的反射”練習(xí)
    接力生命花季無(wú)聲
    ——記王臺(tái)中學(xué)初一158班郭榮軍同學(xué)舍己救人的英雄事跡
    国产 一区精品| 高清日韩中文字幕在线| 一级a做视频免费观看| 久久99热这里只频精品6学生| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜爱| 大香蕉97超碰在线| 亚洲欧美日韩另类电影网站 | 精品久久久久久电影网| 久久久精品欧美日韩精品| 精品久久久久久久末码| 男人狂女人下面高潮的视频| 久久精品国产亚洲av天美| 成年女人在线观看亚洲视频 | 草草在线视频免费看| 麻豆国产97在线/欧美| 日本爱情动作片www.在线观看| 欧美xxxx黑人xx丫x性爽| 国产免费视频播放在线视频| 成人鲁丝片一二三区免费| 在线观看一区二区三区激情| 一个人看的www免费观看视频| 成人毛片a级毛片在线播放| 国产一区二区在线观看日韩| 国产高清国产精品国产三级 | 国产av国产精品国产| 99久久精品国产国产毛片| 又爽又黄无遮挡网站| 亚洲在线观看片| 亚洲国产日韩一区二区| 亚洲真实伦在线观看| 亚洲综合色惰| 国产综合精华液| 亚洲国产色片| 午夜精品国产一区二区电影 | 亚洲天堂国产精品一区在线| 丝袜喷水一区| 国产成人freesex在线| 人妻制服诱惑在线中文字幕| 自拍偷自拍亚洲精品老妇| 亚州av有码| 久久精品久久精品一区二区三区| 在线免费十八禁| 中文字幕免费在线视频6| 日韩欧美一区视频在线观看 | 久久久久久久午夜电影| 国产成人免费无遮挡视频| 午夜爱爱视频在线播放| 午夜精品一区二区三区免费看| 国产欧美亚洲国产| 汤姆久久久久久久影院中文字幕| 少妇人妻久久综合中文| 69人妻影院| 精品人妻熟女av久视频| 精品一区二区三卡| 天堂中文最新版在线下载 | 国产视频首页在线观看| 在线 av 中文字幕| 老女人水多毛片| 日日摸夜夜添夜夜添av毛片| 久久精品人妻少妇| 在线观看美女被高潮喷水网站| a级毛色黄片| 国产精品不卡视频一区二区| 在线观看一区二区三区| 久久午夜福利片| 亚洲精品乱码久久久v下载方式| 亚洲色图av天堂| 亚洲怡红院男人天堂| 男男h啪啪无遮挡| 国产91av在线免费观看| 在线天堂最新版资源| 最近最新中文字幕免费大全7| 国语对白做爰xxxⅹ性视频网站| 涩涩av久久男人的天堂| 午夜福利高清视频| 日本爱情动作片www.在线观看| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 精品一区二区免费观看| 成年人午夜在线观看视频| 欧美一级a爱片免费观看看| 国内精品美女久久久久久| 午夜激情久久久久久久| 99热全是精品| 亚洲性久久影院| 少妇人妻一区二区三区视频| 国产高清有码在线观看视频| 亚洲不卡免费看| 老司机影院毛片| 久久久久久伊人网av| 一本久久精品| 岛国毛片在线播放| 最近最新中文字幕免费大全7| 99久久精品国产国产毛片| 国产探花在线观看一区二区| 国产真实伦视频高清在线观看| 观看美女的网站| 美女视频免费永久观看网站| 欧美最新免费一区二区三区| 2018国产大陆天天弄谢| 国产亚洲一区二区精品| .国产精品久久| 亚洲国产精品专区欧美| 久久久久网色| 美女被艹到高潮喷水动态| 中文字幕av成人在线电影| 91久久精品国产一区二区三区| 青春草国产在线视频| 搞女人的毛片| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频 | 精品午夜福利在线看| 好男人在线观看高清免费视频| 亚洲丝袜综合中文字幕| 国产真实伦视频高清在线观看| 亚洲国产精品成人综合色| 日本三级黄在线观看| 国产av码专区亚洲av| 国产欧美另类精品又又久久亚洲欧美| 午夜激情福利司机影院| 黄色一级大片看看| 亚洲精品乱久久久久久| 亚洲国产色片| 欧美日本视频| 精品久久久久久久久亚洲| 大片免费播放器 马上看| 丝袜喷水一区| 在线观看av片永久免费下载| 久久97久久精品| 直男gayav资源| 国产免费福利视频在线观看| 午夜亚洲福利在线播放| 真实男女啪啪啪动态图| 亚洲真实伦在线观看| av国产精品久久久久影院| 直男gayav资源| 寂寞人妻少妇视频99o| 亚洲国产欧美人成| 又爽又黄无遮挡网站| 亚洲成色77777| 日韩三级伦理在线观看| 免费看av在线观看网站| av天堂中文字幕网| 激情五月婷婷亚洲| 中文字幕制服av| 99视频精品全部免费 在线| 亚洲av福利一区| 青春草国产在线视频| 国产成人精品福利久久| 各种免费的搞黄视频| 国产黄色视频一区二区在线观看| 日韩成人伦理影院| 成人国产av品久久久| 成人毛片a级毛片在线播放| 成人漫画全彩无遮挡| 97精品久久久久久久久久精品| 亚洲精品中文字幕在线视频 | 狂野欧美白嫩少妇大欣赏| 免费黄色在线免费观看| tube8黄色片| 国产欧美亚洲国产| 久久久午夜欧美精品| 97超视频在线观看视频| 天美传媒精品一区二区| 人妻少妇偷人精品九色| av黄色大香蕉| 少妇高潮的动态图| 欧美三级亚洲精品| 成人国产麻豆网| 成人漫画全彩无遮挡| 久久久久精品性色| av黄色大香蕉| 免费大片黄手机在线观看| 久久久午夜欧美精品| 精品久久国产蜜桃| 人妻少妇偷人精品九色| 成人一区二区视频在线观看| 免费观看在线日韩| 久久久久久久午夜电影| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 日本爱情动作片www.在线观看| 免费黄频网站在线观看国产| 成人黄色视频免费在线看| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 成人国产麻豆网| 久久久久精品性色| 看免费成人av毛片| 欧美国产精品一级二级三级 | 少妇熟女欧美另类| 插逼视频在线观看| 久久综合国产亚洲精品| 成人综合一区亚洲| 国产精品av视频在线免费观看| 亚洲四区av| 国产伦在线观看视频一区| 国产精品久久久久久精品电影| 国产精品一区www在线观看| 在线精品无人区一区二区三 | 亚洲最大成人手机在线| 国产精品99久久久久久久久| 黑人高潮一二区| 成年免费大片在线观看| 三级国产精品欧美在线观看| 国产精品久久久久久av不卡| 日韩欧美一区视频在线观看 | 国产 一区精品| 久久国内精品自在自线图片| 日韩,欧美,国产一区二区三区| 日本色播在线视频| 欧美xxxx性猛交bbbb| av在线亚洲专区| 亚洲精品乱久久久久久| 久久久久久久国产电影| 亚洲欧美精品专区久久| 亚洲精品国产色婷婷电影| 99久久人妻综合| 日韩av在线免费看完整版不卡| 最近中文字幕高清免费大全6| 在线观看人妻少妇| 国产午夜精品久久久久久一区二区三区| 最近手机中文字幕大全| 亚洲在久久综合| 色播亚洲综合网| 免费少妇av软件| 国产精品成人在线| 日本一二三区视频观看| 成人亚洲精品一区在线观看 | 好男人视频免费观看在线| 国产伦精品一区二区三区四那| 免费观看a级毛片全部| 2021少妇久久久久久久久久久| 亚洲国产精品999| 七月丁香在线播放| 国产精品蜜桃在线观看| 国内少妇人妻偷人精品xxx网站| 毛片一级片免费看久久久久| 国产精品不卡视频一区二区| 日韩av免费高清视频| 亚洲美女搞黄在线观看| 亚洲精品自拍成人| 色婷婷久久久亚洲欧美| 九九在线视频观看精品| 亚洲,一卡二卡三卡| 亚洲av中文字字幕乱码综合| 一本色道久久久久久精品综合| 三级国产精品片| 免费观看av网站的网址| 国产综合精华液| 国产 一区 欧美 日韩| 永久网站在线| 自拍欧美九色日韩亚洲蝌蚪91 | 97超视频在线观看视频| 高清欧美精品videossex| 国产午夜精品一二区理论片| 国产永久视频网站| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 日日撸夜夜添| 女人被狂操c到高潮| 国内精品美女久久久久久| 精品少妇久久久久久888优播| 午夜精品一区二区三区免费看| 久久久久精品久久久久真实原创| 日韩 亚洲 欧美在线| 大陆偷拍与自拍| 久久99热这里只频精品6学生| 99久久精品一区二区三区| 男女下面进入的视频免费午夜| 在线免费观看不下载黄p国产| 老司机影院毛片| 精品久久久久久久久av| 亚洲精品456在线播放app| 日本免费在线观看一区| 少妇猛男粗大的猛烈进出视频 | 久久久精品免费免费高清| 成人二区视频| 亚洲经典国产精华液单| 嫩草影院新地址| 国产在线男女| 亚洲精品乱久久久久久| 乱码一卡2卡4卡精品| 国产真实伦视频高清在线观看| 久久国内精品自在自线图片| 国产有黄有色有爽视频| 王馨瑶露胸无遮挡在线观看| 伊人久久精品亚洲午夜| 国产精品三级大全| 日本一本二区三区精品| 国产精品国产三级国产av玫瑰| 如何舔出高潮| 91久久精品电影网| 国内精品宾馆在线| 免费观看在线日韩| 久久国内精品自在自线图片| 成人亚洲精品一区在线观看 | 校园人妻丝袜中文字幕| 日韩成人伦理影院| a级毛片免费高清观看在线播放| 97在线人人人人妻| 亚洲最大成人手机在线| 欧美日韩一区二区视频在线观看视频在线 | 国产免费视频播放在线视频| 岛国毛片在线播放| 亚洲色图综合在线观看| 亚洲精华国产精华液的使用体验| 综合色av麻豆| 男女无遮挡免费网站观看| 亚洲国产精品专区欧美| 在线亚洲精品国产二区图片欧美 | 久久精品国产鲁丝片午夜精品| 亚洲成人中文字幕在线播放| 国产 精品1| 久久久久久久久久人人人人人人| 国产成人免费无遮挡视频| 欧美+日韩+精品| 大码成人一级视频| 黄色欧美视频在线观看| 国产成人精品婷婷| 久久久久久久亚洲中文字幕| 免费不卡的大黄色大毛片视频在线观看| 天堂俺去俺来也www色官网| 永久免费av网站大全| 最后的刺客免费高清国语| 看十八女毛片水多多多| 国产精品久久久久久av不卡| 人妻制服诱惑在线中文字幕| 久久久精品欧美日韩精品| 午夜免费观看性视频| 黑人高潮一二区| videos熟女内射| 最新中文字幕久久久久| 中文字幕免费在线视频6| 欧美日韩视频精品一区| 国产亚洲一区二区精品| 亚洲欧美日韩无卡精品| 成人鲁丝片一二三区免费| 久久久午夜欧美精品| 国产美女午夜福利| 亚洲精品日韩在线中文字幕| 麻豆精品久久久久久蜜桃| 久久久亚洲精品成人影院| 五月天丁香电影| 亚洲国产高清在线一区二区三| 欧美日韩国产mv在线观看视频 | 亚洲av男天堂| 99热这里只有精品一区| 搡女人真爽免费视频火全软件| 美女视频免费永久观看网站| 欧美bdsm另类| 精品99又大又爽又粗少妇毛片| 日韩中字成人| 极品教师在线视频| 日韩精品有码人妻一区| 亚洲国产精品国产精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 蜜臀久久99精品久久宅男| 日韩 亚洲 欧美在线| 久久热精品热| 午夜爱爱视频在线播放| 亚洲av国产av综合av卡| 99热这里只有是精品在线观看| 午夜福利在线在线| 青春草国产在线视频| 日产精品乱码卡一卡2卡三| 午夜激情福利司机影院| 亚洲欧美精品专区久久| 如何舔出高潮| 精品国产乱码久久久久久小说| 成人鲁丝片一二三区免费| 91aial.com中文字幕在线观看| 国产亚洲一区二区精品| 3wmmmm亚洲av在线观看| 久久久久久久久久久丰满| 免费播放大片免费观看视频在线观看| 久久精品国产鲁丝片午夜精品| 少妇人妻久久综合中文| 国产精品麻豆人妻色哟哟久久| 亚洲熟女精品中文字幕| 午夜精品国产一区二区电影 | 性色av一级| 亚洲真实伦在线观看| 久久99热6这里只有精品| 国产精品成人在线| 麻豆成人午夜福利视频| 精品酒店卫生间| 欧美另类一区| 国产69精品久久久久777片| 国产黄片美女视频| 一边亲一边摸免费视频| 亚洲人成网站在线播| 国产永久视频网站| 嫩草影院入口| 另类亚洲欧美激情| 日韩不卡一区二区三区视频在线| 亚洲人与动物交配视频| 黄色日韩在线| 18禁在线播放成人免费| 乱码一卡2卡4卡精品| 精品人妻视频免费看| 男人爽女人下面视频在线观看| 色综合色国产| 丰满乱子伦码专区| 精品人妻一区二区三区麻豆| 深夜a级毛片| 黄色视频在线播放观看不卡| 国产精品一区二区三区四区免费观看| 久久久久久久午夜电影| 免费黄频网站在线观看国产| 日本色播在线视频| 午夜激情福利司机影院| 亚洲成色77777| 午夜免费鲁丝| 亚洲人成网站在线观看播放| 国产一区亚洲一区在线观看| 特大巨黑吊av在线直播| 免费高清在线观看视频在线观看| 国产老妇伦熟女老妇高清| h日本视频在线播放| 久久ye,这里只有精品| 精品一区二区三区视频在线| 一级av片app| 中文天堂在线官网| 97热精品久久久久久| 免费大片黄手机在线观看| 高清av免费在线| 最近最新中文字幕大全电影3| 又大又黄又爽视频免费| 国产精品伦人一区二区| 国内少妇人妻偷人精品xxx网站| 一区二区三区免费毛片| 日日摸夜夜添夜夜添av毛片| 色哟哟·www| 99热这里只有是精品50| 亚洲怡红院男人天堂| 欧美日韩亚洲高清精品| 两个人的视频大全免费| 校园人妻丝袜中文字幕| 亚洲性久久影院| 噜噜噜噜噜久久久久久91| 晚上一个人看的免费电影| 亚洲精品成人久久久久久| 丝袜脚勾引网站| 天堂中文最新版在线下载 | 欧美xxxx性猛交bbbb| 你懂的网址亚洲精品在线观看| 女人十人毛片免费观看3o分钟| 内地一区二区视频在线| 免费av观看视频| 麻豆国产97在线/欧美| 中文字幕免费在线视频6| 男人舔奶头视频| 卡戴珊不雅视频在线播放| www.av在线官网国产| 一本色道久久久久久精品综合| 韩国高清视频一区二区三区| 男人舔奶头视频| 精品人妻熟女av久视频| 九九在线视频观看精品| av卡一久久| 三级经典国产精品| 偷拍熟女少妇极品色| 欧美97在线视频| 一级片'在线观看视频| 欧美一级a爱片免费观看看| 亚洲av福利一区| 免费人成在线观看视频色| 五月玫瑰六月丁香| av在线老鸭窝| 日韩欧美精品v在线| 亚洲成人一二三区av| 久久精品人妻少妇| 毛片女人毛片| 内射极品少妇av片p| 日本爱情动作片www.在线观看| 日本熟妇午夜| 国产黄a三级三级三级人| 人妻 亚洲 视频| 大片电影免费在线观看免费| 欧美一级a爱片免费观看看| 日韩制服骚丝袜av| 国产久久久一区二区三区| 99久久精品热视频| 色视频www国产| 寂寞人妻少妇视频99o| 国产人妻一区二区三区在| 亚洲av欧美aⅴ国产| 老司机影院成人| 永久网站在线| 99久国产av精品国产电影| 国产熟女欧美一区二区| 大香蕉97超碰在线| 久久女婷五月综合色啪小说 | 欧美一级a爱片免费观看看| av又黄又爽大尺度在线免费看| 亚洲图色成人| 少妇人妻一区二区三区视频| 日本黄色片子视频| 日韩欧美精品v在线| 51国产日韩欧美| 91午夜精品亚洲一区二区三区| 成年人午夜在线观看视频| 1000部很黄的大片| 男插女下体视频免费在线播放| 欧美日本视频| 在线观看美女被高潮喷水网站| 麻豆成人av视频| 久久久久久久精品精品| 一个人看视频在线观看www免费| 2018国产大陆天天弄谢| 人人妻人人看人人澡| 成人漫画全彩无遮挡| 91狼人影院| 天美传媒精品一区二区| 网址你懂的国产日韩在线| 亚洲国产精品成人综合色| 伦理电影大哥的女人| 王馨瑶露胸无遮挡在线观看| 久久鲁丝午夜福利片| 天美传媒精品一区二区| 人人妻人人看人人澡| 久久影院123| 精品少妇黑人巨大在线播放| 嫩草影院入口| 国产成人精品福利久久| 一二三四中文在线观看免费高清| 国产日韩欧美亚洲二区| 亚洲精品成人久久久久久| 成人鲁丝片一二三区免费| 欧美精品一区二区大全| 亚洲,一卡二卡三卡| 小蜜桃在线观看免费完整版高清| 久久久久久国产a免费观看| 国产伦理片在线播放av一区| 久久ye,这里只有精品| 欧美97在线视频| 精品久久国产蜜桃| 日日摸夜夜添夜夜添av毛片| 亚洲精品久久午夜乱码| 免费大片黄手机在线观看| av国产久精品久网站免费入址| 国产黄片美女视频| 午夜福利视频1000在线观看| 交换朋友夫妻互换小说| 国产乱来视频区| 99热国产这里只有精品6| 中文字幕制服av| 精品国产三级普通话版| 日本一本二区三区精品| 国国产精品蜜臀av免费| 午夜视频国产福利| 国国产精品蜜臀av免费| 一个人看视频在线观看www免费| 免费观看无遮挡的男女| 成人高潮视频无遮挡免费网站| videossex国产| 国内精品美女久久久久久| 国产精品熟女久久久久浪| 永久免费av网站大全| 精品人妻熟女av久视频| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 精品人妻视频免费看| 成年av动漫网址| 日日撸夜夜添| 欧美潮喷喷水| 成人漫画全彩无遮挡| 久久精品国产亚洲av涩爱| 一个人观看的视频www高清免费观看| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 国产成人精品福利久久| 精品久久久久久久末码| 身体一侧抽搐| 久久人人爽人人片av| 日产精品乱码卡一卡2卡三| 嘟嘟电影网在线观看| 精品人妻视频免费看| 久久久久久久久久久丰满| xxx大片免费视频| 尾随美女入室| 天美传媒精品一区二区| 人体艺术视频欧美日本| 简卡轻食公司| 一区二区三区精品91| 国产精品一二三区在线看| 看非洲黑人一级黄片| 亚洲国产精品国产精品| 啦啦啦在线观看免费高清www| 联通29元200g的流量卡| 一级a做视频免费观看| 欧美老熟妇乱子伦牲交| 亚洲av一区综合| 亚洲最大成人中文| 91精品伊人久久大香线蕉| 美女脱内裤让男人舔精品视频| 最近2019中文字幕mv第一页| 自拍偷自拍亚洲精品老妇| 欧美高清性xxxxhd video| 久久99精品国语久久久| 尤物成人国产欧美一区二区三区| av免费观看日本| 亚洲av在线观看美女高潮| 97超视频在线观看视频| 国产精品国产av在线观看| 亚洲色图综合在线观看| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区黑人 | 免费av不卡在线播放| 亚洲国产精品成人久久小说| 日韩欧美精品v在线| 中文字幕人妻熟人妻熟丝袜美|