• <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é)舍己救人的英雄事跡
    精品久久蜜臀av无| 午夜成年电影在线免费观看| 一区福利在线观看| 在线观看一区二区三区| 亚洲精品国产精品久久久不卡| 午夜免费观看网址| 国产99白浆流出| 麻豆成人av在线观看| 亚洲成人久久爱视频| 变态另类成人亚洲欧美熟女| 日本免费一区二区三区高清不卡| 国产精品久久久久久亚洲av鲁大| 国产一区二区在线av高清观看| 99久久无色码亚洲精品果冻| 免费在线观看亚洲国产| 天堂av国产一区二区熟女人妻| 亚洲精品一卡2卡三卡4卡5卡| 日本五十路高清| 男女床上黄色一级片免费看| 色精品久久人妻99蜜桃| 首页视频小说图片口味搜索| 韩国av一区二区三区四区| 国产午夜福利久久久久久| 久久国产乱子伦精品免费另类| 哪里可以看免费的av片| 男人和女人高潮做爰伦理| 亚洲精品乱码久久久v下载方式 | 久久久精品大字幕| 午夜两性在线视频| 1024香蕉在线观看| 十八禁人妻一区二区| 特大巨黑吊av在线直播| 免费一级毛片在线播放高清视频| 亚洲国产高清在线一区二区三| 亚洲国产看品久久| 搡老妇女老女人老熟妇| 国产69精品久久久久777片 | 最近视频中文字幕2019在线8| 国产一区二区激情短视频| 老司机午夜福利在线观看视频| 在线免费观看的www视频| 亚洲 国产 在线| 亚洲欧洲精品一区二区精品久久久| 嫩草影院入口| 国产精品综合久久久久久久免费| 免费无遮挡裸体视频| 人妻丰满熟妇av一区二区三区| 国产久久久一区二区三区| 一级作爱视频免费观看| 中文字幕人成人乱码亚洲影| 不卡av一区二区三区| 一二三四在线观看免费中文在| 亚洲国产高清在线一区二区三| 午夜福利高清视频| 午夜福利欧美成人| av中文乱码字幕在线| 国产精品香港三级国产av潘金莲| 91麻豆精品激情在线观看国产| 欧美又色又爽又黄视频| 亚洲午夜理论影院| 国产精品电影一区二区三区| 久久国产精品人妻蜜桃| 久久热在线av| 免费在线观看日本一区| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 中国美女看黄片| 99热这里只有精品一区 | 婷婷精品国产亚洲av| 亚洲精品在线观看二区| 国产精品女同一区二区软件 | 国产一区在线观看成人免费| 黄频高清免费视频| 99久久综合精品五月天人人| 曰老女人黄片| 欧美午夜高清在线| 成年人黄色毛片网站| 伦理电影免费视频| 日韩中文字幕欧美一区二区| 床上黄色一级片| 国产欧美日韩一区二区精品| 国产精品永久免费网站| 亚洲第一欧美日韩一区二区三区| 欧美+亚洲+日韩+国产| 又黄又粗又硬又大视频| 久久午夜综合久久蜜桃| 黄色丝袜av网址大全| 欧美国产日韩亚洲一区| 国产精品一区二区免费欧美| 国内精品美女久久久久久| 色吧在线观看| 国产精品 国内视频| 欧美3d第一页| 亚洲av五月六月丁香网| 他把我摸到了高潮在线观看| 亚洲成av人片免费观看| 99热精品在线国产| 欧美又色又爽又黄视频| 搡老妇女老女人老熟妇| 丰满人妻熟妇乱又伦精品不卡| 国产毛片a区久久久久| 国产一区在线观看成人免费| 99riav亚洲国产免费| 久久久久精品国产欧美久久久| 婷婷精品国产亚洲av| 草草在线视频免费看| 国产成人精品久久二区二区91| 两人在一起打扑克的视频| netflix在线观看网站| 成人av一区二区三区在线看| www日本黄色视频网| 日韩国内少妇激情av| 日韩欧美免费精品| 极品教师在线免费播放| 又大又爽又粗| 美女 人体艺术 gogo| 国产精品 欧美亚洲| 国产精品野战在线观看| 搡老熟女国产l中国老女人| 久久香蕉精品热| 亚洲av片天天在线观看| 人妻久久中文字幕网| 精品国产乱子伦一区二区三区| 99久久99久久久精品蜜桃| 日本黄色视频三级网站网址| 午夜福利18| 岛国视频午夜一区免费看| 国产精品一及| 99精品欧美一区二区三区四区| 九九热线精品视视频播放| 国产精品久久久久久人妻精品电影| 最近在线观看免费完整版| 亚洲欧美精品综合久久99| 亚洲avbb在线观看| 麻豆久久精品国产亚洲av| 久久精品综合一区二区三区| 亚洲男人的天堂狠狠| 在线a可以看的网站| 免费电影在线观看免费观看| 可以在线观看的亚洲视频| 9191精品国产免费久久| 久久99热这里只有精品18| 亚洲五月婷婷丁香| 国产亚洲精品一区二区www| 91av网一区二区| 男女视频在线观看网站免费| 国产激情欧美一区二区| 五月伊人婷婷丁香| 久久午夜亚洲精品久久| 搞女人的毛片| 波多野结衣高清作品| 女警被强在线播放| 婷婷精品国产亚洲av| 国产成年人精品一区二区| bbb黄色大片| 亚洲自偷自拍图片 自拍| avwww免费| 男人的好看免费观看在线视频| 免费在线观看日本一区| 亚洲成人免费电影在线观看| 久久香蕉精品热| 制服丝袜大香蕉在线| 天天躁日日操中文字幕| 一进一出好大好爽视频| 久久这里只有精品中国| 久久欧美精品欧美久久欧美| 黄色丝袜av网址大全| 法律面前人人平等表现在哪些方面| 亚洲国产日韩欧美精品在线观看 | 久久久久国内视频| 黄色片一级片一级黄色片| 两性夫妻黄色片| 在线免费观看不下载黄p国产 | 丰满人妻一区二区三区视频av | 在线国产一区二区在线| 一个人看的www免费观看视频| 国产精品乱码一区二三区的特点| 女同久久另类99精品国产91| 成年女人永久免费观看视频| 悠悠久久av| 日韩欧美 国产精品| 91在线观看av| 国产成人精品久久二区二区免费| 18禁黄网站禁片午夜丰满| 国产亚洲欧美98| 色av中文字幕| 色综合欧美亚洲国产小说| 成年免费大片在线观看| 国产人伦9x9x在线观看| 黄色丝袜av网址大全| 色综合婷婷激情| 久99久视频精品免费| 日韩成人在线观看一区二区三区| 999精品在线视频| 国产高清视频在线播放一区| 国产乱人视频| 精品人妻1区二区| 久久中文字幕人妻熟女| 老司机午夜十八禁免费视频| 天天一区二区日本电影三级| 天天躁狠狠躁夜夜躁狠狠躁| 男人和女人高潮做爰伦理| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| 亚洲av免费在线观看| 国产精品av久久久久免费| 国产久久久一区二区三区| 亚洲色图av天堂| 国产真人三级小视频在线观看| 日韩欧美免费精品| 性欧美人与动物交配| 动漫黄色视频在线观看| 国产成人福利小说| 一进一出抽搐动态| 一本综合久久免费| 1024香蕉在线观看| 欧美色视频一区免费| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看日本一区| 九九热线精品视视频播放| 亚洲国产精品合色在线| 亚洲午夜精品一区,二区,三区| 又黄又爽又免费观看的视频| 女同久久另类99精品国产91| 亚洲人成网站在线播放欧美日韩| 国内精品美女久久久久久| 免费在线观看视频国产中文字幕亚洲| 午夜福利18| 99riav亚洲国产免费| 亚洲天堂国产精品一区在线| 亚洲国产高清在线一区二区三| 十八禁网站免费在线| 天天添夜夜摸| 在线十欧美十亚洲十日本专区| 久久久国产成人精品二区| 男人的好看免费观看在线视频| 色播亚洲综合网| 男插女下体视频免费在线播放| 1024香蕉在线观看| 精品久久久久久久久久免费视频| 欧美在线黄色| 免费在线观看影片大全网站| 成人鲁丝片一二三区免费| av女优亚洲男人天堂 | 亚洲中文字幕一区二区三区有码在线看 | 欧美在线黄色| cao死你这个sao货| 久久精品人妻少妇| 俄罗斯特黄特色一大片| 成人18禁在线播放| 女人高潮潮喷娇喘18禁视频| 精品一区二区三区视频在线 | 99热6这里只有精品| 国产一级毛片七仙女欲春2| 欧美日韩乱码在线| www.自偷自拍.com| 九九在线视频观看精品| 不卡av一区二区三区| 欧美最黄视频在线播放免费| www国产在线视频色| 中国美女看黄片| 香蕉av资源在线| 欧美精品啪啪一区二区三区| 男插女下体视频免费在线播放| 国产又色又爽无遮挡免费看| 国内精品一区二区在线观看| 国内精品久久久久精免费| 精品午夜福利视频在线观看一区| 中文字幕熟女人妻在线| 男女午夜视频在线观看| 亚洲午夜精品一区,二区,三区| 国产97色在线日韩免费| 深夜精品福利| 久久国产精品影院| 国内精品久久久久精免费| 我要搜黄色片| 国产99白浆流出| 亚洲男人的天堂狠狠| 少妇裸体淫交视频免费看高清| 欧美日韩国产亚洲二区| 欧美乱妇无乱码| 久久久久久九九精品二区国产| 免费在线观看亚洲国产| 舔av片在线| 亚洲欧美日韩无卡精品| 欧美日韩综合久久久久久 | 欧美日韩福利视频一区二区| 香蕉av资源在线| 精品国产亚洲在线| 99国产精品99久久久久| 99热这里只有是精品50| 国产亚洲欧美98| 免费搜索国产男女视频| 国产精品久久久av美女十八| 在线国产一区二区在线| 成人18禁在线播放| 久久草成人影院| 国产美女午夜福利| 黄频高清免费视频| 亚洲精品粉嫩美女一区| 1024香蕉在线观看| 一a级毛片在线观看| 动漫黄色视频在线观看| 亚洲午夜精品一区,二区,三区| svipshipincom国产片| 精品久久蜜臀av无| 亚洲人成电影免费在线| 午夜亚洲福利在线播放| 国产成人精品无人区| 亚洲欧美日韩高清专用| 操出白浆在线播放| 日本在线视频免费播放| 亚洲欧洲精品一区二区精品久久久| 一级毛片高清免费大全| 日日摸夜夜添夜夜添小说| 亚洲在线观看片| 老熟妇乱子伦视频在线观看| 亚洲最大成人中文| 亚洲色图av天堂| 精品国产超薄肉色丝袜足j| 在线视频色国产色| 黄色 视频免费看| 国产乱人伦免费视频| 成人鲁丝片一二三区免费| 亚洲国产欧美一区二区综合| 亚洲精品国产精品久久久不卡| 此物有八面人人有两片| 国产精品国产高清国产av| 精品久久蜜臀av无| 亚洲aⅴ乱码一区二区在线播放| 久久中文看片网| 嫩草影视91久久| 国产精品久久视频播放| 久久久久久久久久黄片| av在线蜜桃| 一级毛片女人18水好多| 97超级碰碰碰精品色视频在线观看| 久久久久九九精品影院| 嫩草影院入口| 亚洲18禁久久av| avwww免费| 中文字幕人成人乱码亚洲影| 国产三级在线视频| 12—13女人毛片做爰片一| 少妇熟女aⅴ在线视频| 18禁裸乳无遮挡免费网站照片| 精品国内亚洲2022精品成人| 嫩草影院入口| 精品一区二区三区视频在线观看免费| 亚洲人成网站高清观看| 特大巨黑吊av在线直播| 偷拍熟女少妇极品色| 国产精品九九99| 国产成年人精品一区二区| 中文字幕人成人乱码亚洲影| 亚洲av免费在线观看| 岛国在线免费视频观看| 很黄的视频免费| 久久久久久久久免费视频了| 亚洲午夜精品一区,二区,三区| 波多野结衣高清无吗| 国产午夜精品久久久久久| 两性夫妻黄色片| 国产探花在线观看一区二区| 少妇人妻一区二区三区视频| 嫩草影院入口| 欧美日本视频| 国产三级在线视频| 亚洲欧美精品综合久久99| 国产视频一区二区在线看| 九九热线精品视视频播放| 校园春色视频在线观看| 国产野战对白在线观看| 禁无遮挡网站| 欧美高清成人免费视频www| 一a级毛片在线观看| 精品一区二区三区视频在线 | 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 大型黄色视频在线免费观看| 日韩高清综合在线| 精品国内亚洲2022精品成人| 欧美一区二区精品小视频在线| 久久亚洲精品不卡| www日本在线高清视频| 婷婷亚洲欧美| 日韩 欧美 亚洲 中文字幕| 黄色视频,在线免费观看| 午夜福利视频1000在线观看| 欧美黄色淫秽网站| 波多野结衣巨乳人妻| 国产99白浆流出| 欧美黑人欧美精品刺激| 黄色视频,在线免费观看| 亚洲精品美女久久久久99蜜臀| 亚洲成人中文字幕在线播放| www日本黄色视频网| 午夜a级毛片| 麻豆国产97在线/欧美| 淫秽高清视频在线观看| 在线观看午夜福利视频| 欧美极品一区二区三区四区| 啦啦啦观看免费观看视频高清| 黑人操中国人逼视频| 亚洲 欧美 日韩 在线 免费| 日本成人三级电影网站| 人人妻人人澡欧美一区二区| 男女床上黄色一级片免费看| 久久精品91无色码中文字幕| 精品无人区乱码1区二区| 老熟妇乱子伦视频在线观看| av在线蜜桃| 国产精品av视频在线免费观看| 在线观看美女被高潮喷水网站 | 国产精品女同一区二区软件 | 一本综合久久免费| 白带黄色成豆腐渣| 亚洲片人在线观看| 露出奶头的视频| 亚洲午夜理论影院| 91字幕亚洲| 在线免费观看的www视频| 在线永久观看黄色视频| 嫩草影院入口| 天堂影院成人在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利欧美成人| 久久中文字幕一级| 天堂网av新在线| 黄色 视频免费看| 日韩三级视频一区二区三区| 观看美女的网站| 最新在线观看一区二区三区| 亚洲成人免费电影在线观看| 91麻豆精品激情在线观看国产| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲18禁久久av| 欧美极品一区二区三区四区| 男人舔女人的私密视频| 91久久精品国产一区二区成人 | 亚洲一区二区三区色噜噜| 国产精品98久久久久久宅男小说| 少妇丰满av| a在线观看视频网站| 免费在线观看日本一区| 我的老师免费观看完整版| 五月玫瑰六月丁香| xxxwww97欧美| 欧美日韩瑟瑟在线播放| 亚洲激情在线av| 午夜福利在线观看吧| 国产高清三级在线| 老汉色∧v一级毛片| 亚洲精品国产精品久久久不卡| 久久久久性生活片| 国内精品久久久久精免费| av天堂在线播放| 日韩av在线大香蕉| 成人永久免费在线观看视频| 亚洲av免费在线观看| 免费看日本二区| 久久久久亚洲av毛片大全| 黄片大片在线免费观看| 成年版毛片免费区| 男女之事视频高清在线观看| 99久久成人亚洲精品观看| 亚洲18禁久久av| 久久久久精品国产欧美久久久| 久久国产乱子伦精品免费另类| 亚洲人成电影免费在线| 熟女电影av网| 99riav亚洲国产免费| 国产精品久久久久久久电影 | 精品国产美女av久久久久小说| 性欧美人与动物交配| 12—13女人毛片做爰片一| 男女视频在线观看网站免费| 香蕉久久夜色| 日韩欧美国产在线观看| 色综合婷婷激情| 亚洲精品美女久久av网站| 久久久久精品国产欧美久久久| 少妇丰满av| 免费在线观看成人毛片| 色噜噜av男人的天堂激情| 日本黄色片子视频| 日韩精品青青久久久久久| 色播亚洲综合网| 国产免费男女视频| 国产成人福利小说| 欧美绝顶高潮抽搐喷水| 99国产综合亚洲精品| 日日摸夜夜添夜夜添小说| 午夜福利欧美成人| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 一区二区三区高清视频在线| 国产精品自产拍在线观看55亚洲| 成人亚洲精品av一区二区| 免费av不卡在线播放| 黄色女人牲交| 狠狠狠狠99中文字幕| 午夜精品一区二区三区免费看| 91字幕亚洲| 狠狠狠狠99中文字幕| 亚洲avbb在线观看| 午夜福利在线观看免费完整高清在 | 日韩欧美一区二区三区在线观看| 亚洲色图 男人天堂 中文字幕| 一卡2卡三卡四卡精品乱码亚洲| h日本视频在线播放| 欧美黄色片欧美黄色片| aaaaa片日本免费| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 久久国产精品人妻蜜桃| 悠悠久久av| 1024手机看黄色片| 老熟妇乱子伦视频在线观看| 久久久久久久精品吃奶| 国产欧美日韩精品一区二区| 亚洲国产欧洲综合997久久,| 久久久国产成人精品二区| 国产精品亚洲一级av第二区| 久久人人精品亚洲av| 中亚洲国语对白在线视频| 成人一区二区视频在线观看| 美女大奶头视频| 老司机午夜福利在线观看视频| 亚洲 欧美 日韩 在线 免费| 欧美日韩乱码在线| 日本熟妇午夜| 国产伦精品一区二区三区视频9 | 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| 国产精品亚洲美女久久久| 久久久久久久久久黄片| 99久久国产精品久久久| 麻豆久久精品国产亚洲av| а√天堂www在线а√下载| 国产av一区在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| 久久午夜亚洲精品久久| 亚洲av五月六月丁香网| 免费人成视频x8x8入口观看| 亚洲av成人一区二区三| 精品久久久久久久毛片微露脸| 国产久久久一区二区三区| 日本 av在线| 非洲黑人性xxxx精品又粗又长| 熟女电影av网| 国产伦人伦偷精品视频| 中文资源天堂在线| 国产精品一及| 中文字幕精品亚洲无线码一区| 男人舔奶头视频| 亚洲国产欧美一区二区综合| 看片在线看免费视频| 少妇人妻一区二区三区视频| 国产精品久久视频播放| 精品国产亚洲在线| 亚洲aⅴ乱码一区二区在线播放| 丁香六月欧美| 亚洲av免费在线观看| 精品一区二区三区四区五区乱码| 老司机在亚洲福利影院| 国产不卡一卡二| 激情在线观看视频在线高清| 男女做爰动态图高潮gif福利片| 在线播放国产精品三级| 一区二区三区国产精品乱码| 国产精品久久久久久精品电影| 天堂影院成人在线观看| 亚洲成人中文字幕在线播放| 嫁个100分男人电影在线观看| av在线蜜桃| 天天躁日日操中文字幕| 怎么达到女性高潮| 亚洲av电影在线进入| 亚洲天堂国产精品一区在线| 国产精品香港三级国产av潘金莲| 男人和女人高潮做爰伦理| 国产精品99久久久久久久久| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 亚洲国产精品999在线| 亚洲色图 男人天堂 中文字幕| 嫩草影视91久久| 天天躁日日操中文字幕| 亚洲国产日韩欧美精品在线观看 | 亚洲精品一区av在线观看| 国产精品女同一区二区软件 | 欧美黑人欧美精品刺激| 法律面前人人平等表现在哪些方面| 欧美激情久久久久久爽电影| 亚洲专区国产一区二区| 夜夜看夜夜爽夜夜摸| 麻豆国产av国片精品| 成人国产综合亚洲| 亚洲精品粉嫩美女一区| 国产高清videossex| 亚洲国产欧洲综合997久久,| 色尼玛亚洲综合影院| 老熟妇仑乱视频hdxx| 亚洲欧美激情综合另类| 亚洲成人久久性| 最近最新免费中文字幕在线| 国产视频一区二区在线看| 免费大片18禁| 99久久精品一区二区三区| 丰满人妻一区二区三区视频av | 一二三四在线观看免费中文在| 一个人免费在线观看的高清视频| 国产av一区在线观看免费|