999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于情境感知注意機制的出租車需求深度學習預測模型

2023-08-10 07:01:22樊云閣呂玉輝韓紅旗
計算機應用與軟件 2023年7期
關鍵詞:區域模型

樊云閣 呂玉輝 韓紅旗

1(河南農業職業學院信息工程學院 河南 鄭州 451450) 2(洛陽理工學院計算機學院 河南 洛陽 471023) 3(中國科學技術信息研究所 北京 100038)

0 引 言

城市出租車作為一種綠色公共交通工具,為廣大乘客提供舒適、便捷的服務[1]。然而,出租汽車服務中存在的一個主要問題是搭乘服務供給與乘客需求之間的時空平衡。為了保持城市道路交通供給與乘客需求的平衡,準確預測城市出租車需求具有重要意義[2]。

出租車需求與包括道路網、興趣點、氣象及假日資料等數據之間的相關性已經得到了很好的研究,但是,要做出準確的預測還存在一些挑戰[3]。首先,出租車需求具有高度的動態性,即隨著時間的推移變化迅速且有一定的隨機性;第二,城市不同功能區的出租車需求變化趨勢存在顯著差異;第三,出租車需求受一些外部因素的影響較大,如景點、天氣狀況和節假日等[4]。

為了預測出租車需求,學者們提出了一系列方法。例如,利用自回歸綜合移動平均(ARIMA)及其改進方法,將出租車需求預測問題建模為時間序列預測[5]。文獻[6]提出一個具有高維特征的統一線性回歸模型來預測各地區的出租車需求。由于缺乏非線性建模能力,這些方法通常預測精度較低。近年來,深層神經網絡(DNN)被引入到出租車需求預測中,它提取了復雜的時空相關性,并利用外部因素輔助預測。例如,文獻[7]提出了一種基于長短時記憶(LSTM)網絡的序列學習框架,利用瞬時空間相關性,將卷積運算與LSTM相結合。文獻[8]同時提取空間相關性和時間相關性,利用外部因素進一步提高預測精度。Tang等[9]利用模型訓練過程中學習到的固定參數矩陣,嘗試將時空相關性和外部因素結合起來。以上現有研究成果都局限于將不同的時空特征和外部因素結合在一起,這是因為在沒有考慮到情景信息的影響下對參數設置了固定的值,僅限于通過理想化的出租車需求多視角特征以及單純外部因素的情境感知情形下,才可以簡單地確定網絡中的可學習參數,所以并未充分挖掘數據中蘊含的信息,并且適用范圍也有限。

為解決上述問題,本文提出一種基于改進遞歸神經網絡的情境感知注意預測算法,用于預測城市出租車需求。通過紐約市的多個實際數據集實驗驗證了本文方法能夠有效提升預測精度,方法更具一般性。

1 問題描述

定義1(道路網絡) 市區的路網由一組路段組成,每個路段都有兩個終點,即十字路口的交叉點相關聯,并通過共享相同的終點與其他路段相連,所有路段均以圖形格式構成道路網絡。

為了在空間和時間維度上更完善地描述細粒度的出租車需求,將時間離散化為一組等間隔的時間段,用τ={t1,t2,…,tτ,…}表示,其中tτ表示當前時間段。還根據其道路網絡將整個市區劃分為不相連的區域。每個區域都是被多個路段包圍的不規則多邊形。區域集由R={r1,r2,…,rN}表示,其中N代表第幾個區域。根據時間段和區域的定義,進一步提出了出租車的細粒度需求的正式定義如下。

定義2(出租車需求) 使用Xn,τ表示在時間段為tτ∈τ、區域為rn∈R的出租車乘客人數。然后,將在時間段tτ的出租車需求定義為Xτ=[X1,τ,X2,τ,…,XN,τ]。

定義3(出租車行程記錄) 用{tr}表示一組出租車歷史行程記錄。每個記錄tr包含上下車的位置和具體時間,可以用元組tr=(tr.pl,tr.pt,tr.dl,tr.dt)表示。在此元組中,用緯度和經度表示上車地點tr.pl和下車地點tr.dl,而用年月日、小時和分鐘的形式表示具體的上車時間tr.pt和下車時間tr.dt。

為了預測未來的出租車需求,首先從定義3中的歷史出租車行程記錄中提取出歷史時間段的出租車需求。給定一個市區出租車行程記錄的數據集,可以根據某區域rn和某時間點tτ的歷史出租車需求,近似估計得到出租車的乘客數量,公式如下:

Xn,τ=|{tr|tr.pl∈rn∧tr.pt∈tτ}|

(1)

式中:tr.pl∈rn和tr.pt∈tτ分別代表記錄tr在區域rn與時間段tτ的具體乘車地點和乘車時間;函數|·|表示集合的基數。

定義4(興趣點) 興趣點(POI)即位于市區的場所,如購物中心。每個興趣點pi與一個位置pi.l和一個范疇pi.c∈C相關聯,C表示類別集。

(2)

式中:pi.l∈rn表示POI的pi位于區域rn。此外,還為每個區域rn定義了一組鄰域集合rn.neig,該區域由相鄰區域組成。

2 數據描述與分析

2.1 道路網絡數據

本文采用了2016年調查的紐約市道路網絡的數據集[10],該數據集包含87 898個十字路口和91 649個路段。每個交叉點都是矢量,按其緯度和經度定位。利用文獻[11]中提出的基于形態學的地圖分割方法,根據路網將紐約市的整個區域劃分為972個不相交的區域。

2.2 出租車行程數據

紐約市的出租車和轎車委員會已經開放了紐約市出租車出行記錄的數據集[12],其中包含出租車的詳細駕駛信息,例如每輛出租車的上下車地點和具體時間。2016年1月1日至2016年6月30日(130個工作日和52個周末)的出租車行程數據記錄了87 866 988次旅行,覆蓋紐約市的972個地區。

圖1分別描繪了三月和四月各地區工作日和周末的出租車平均需求量。一個地區的顏色越深,該地區的出租車需求就越多。可以發現,不同地區的出租車需求之間存在嚴重的不平衡。此外,工作日和周末的出租車需求有顯著差異。

(a) 工作日 (b) 周末 圖1 紐約市平均出租車需求量

圖2和圖3分別顯示了短期(一天)和長期(一周)出租車需求的周期性,其中分別繪制了格林威治村和中央公園兩個地區的出租車需求。

(a) 格林威治村

(a) 格林威治村

尤其是圖2顯示了三月中不同的四天中(包含兩個工作日和兩個周末)每小時出租車需求變化。不難發現,周三和周四兩個地區的出租車需求具有相似的變化趨勢,而周六和周日也是如此。這表明每個地區的出租車需求在短期內具有較強的周期性。此外,由于人們在工作日和周末的出行方式不同,各地區在工作日和周末的周期變化趨勢顯然也不一樣。

同樣,從圖3中可以看出,由于兩個地區每周三出租車需求的變化趨勢相似,因此出租車需求在長期內具有較強的周期性。短期和長期周期性的存在可以幫助預測未來每個地區的出租車需求。

2.3 POI數據

POI中包含的信息,如一個區域中不同類別的POI數量,反映了該區域的社會功能,這也與出租車需求相關。在紐約使用了一個有670 916個POI的POI數據集,分為16個不同的類別,如表1所示。

表1 興趣點的類別

圖4(a)和圖4(b)分別描繪了藝術與娛樂類POI的地理分布以及周六中午12點至下午6點的出租車需求。可以發現,在某些地區(虛線圓圈標出部分),某些類別的POI數量與某些時期的出租車需求有很強的相關性。

(a) POI:藝術娛樂

為了進一步探索出租車需求與POI之間的關系,以兩個標志性建筑為例:哥倫比亞大學(位于347號區域)和紐約城市大學(位于382號區域)。將POI的比例分為六個類別,即不同企業間對比、專業人員和服務、家鄉和家庭、食品與餐飲、教育、購物,跟蹤出租車需求在一天內的變化如圖5所示。可以發現這兩個地標具有相似的POI分布和相似的出租車需求變化趨勢。通過這種觀察,可以利用POI數據來輔助預測出租車需求。

(a) POI比例

2.4 氣象和假期數據

除POI之外,出租車需求還受到其他外部因素的影響,例如天氣狀況和節假日。還在工作中使用了2016年1月1日至2016年6月30日的氣象記錄數據集,其中包含每六小時記錄的天氣狀況,例如,晴天和陰雨、溫度、風速和濕度等信息。

圖6描繪了紐約市在不同天氣下一天的總體出租車需求,這表明天氣狀況對出租車需求有很大的影響。特別地,從圖6(a)中可以發現,暴風雪后出租車需求顯著減少。相反,圖6(b)顯示了當早上7點到9點下雨時,更多的人會選擇乘坐出租車。

3 模型設計

3.1 CACRNN模型概述

圖7提供了前文提出的深度神經網絡模型的概述,該模型包含四個模塊。

圖7 深度神經網絡模型結構

(1) 實時時空模塊。該模塊由一系列局部卷積(LC)層和門控循環單元(GRU)組成,這兩個層用于在封閉時段內提取出租車需求的時空相關性。特別地,它以o個連續時間段中的出租車需求為輸入,由Yi=[Xτ+1-0,Xτ+2-0,…,Xτ]∈R0×N表示,并預測出下一個時間段fl∈Ro×N中出租車需求。

(2) 短期周期模塊。該模塊根據短期周期的出租車需求進行預測。通過GRU來學習短期周期性,以p為周期性時段中的出租車需求序列、間隔Δs作為輸入,即Ys=[Xτ+1-pΔs,Xτ+1-(p-1)Δs,…,Xτ+1-Δs]∈Rp×N。此外,該模塊還考慮了一些外部因素,如天氣和假日對周期性的影響。將tτ中的外部因素的特征表示為向量uτ∈Rl×ω,并將其與出租車需求連接起來作為GRU的輸入。在此基礎上,預測出了出租車需求fs∈Rl×N。

(3) 長期周期模塊。該模塊繪制了長期(如幾周)周期性模式的出租車需求。與最后一個模塊類似,將出租車需求序列Yl=[Xτ+1-qΔl,Xτ+1-(q-1)Δl,…,Xτ+1-Δl]∈Rq×N與外部因素的特征相結合輸入GRU,并預測出時刻tτ+1處的出租車需求,用fl∈Rl×N表示。

(4) 情景感知注意模塊。利用此模塊將上述模塊的輸出合并到最終的出租車需求預測中,這是一種新的嘗試。特別地,此模塊可以根據天氣信息、時刻索引和區域功能等上下文信息,為每個區域的預測分配不同的權重。

3.2 實時時空模塊

該模塊的結構基于局部卷積層和GRU,用于提取區域內的潛在空間相關性和特定周期內的時間相關性。

(3)

然后,基于rn和rm的函數向量,計算了它們之間的函數相似性。特別使用了兩個向量的加權余弦值來度量相似度,即:

(4)

式中:A(·)表示給定區域的面積。式(4)表明,在出租車需求預測中,POI密度較高的區域對其功能相似區域的影響更大。

(5)

3.3 短期/長期周期性模塊

這兩個模塊分別對出租車需求的短、長期周期性進行預測。如圖8所示,它們共享一個基于GRU的結構,以出租車需求和外部因素特征的組合作為輸入。

圖8 感知機網絡結構

時間維度上的外部因素,如天氣狀況和節假日、周末,對出租車需求的周期性有很大影響。例如,通過調查發現,周末或節假日出租車需求的每小時變化與工作日有顯著差異。此外,雨天和晴天出租車需求高峰時段也不同。為了捕捉外部因素的特征,采用了嵌入方法將這些因素在每個時刻的值轉換為外部特征向量[13],用ut表示。這種嵌入方法被廣泛用于將分類值映射到低維向量。

然后,將出租車需求序列與相應時間內的外部特征向量串聯起來作為GRU的輸入,并通過FC層將輸出轉化為下一時間段內出租車需求的預測。短期周期模塊預測結果的公式定義如下:

fs=FC(GRU(Xτ+1-pΔs⊕uτ+1-pΔs,

Xτ+1-(p-1)Δs⊕uτ+1-(p-1)Δs,…,Xτ+1-Δs⊕uτ+1-Δs))

(6)

fl=FC(GRU(Xτ+1-qΔl⊕uτ+1-qΔl

Xτ+1-(q-1)Δl⊕uτ+1-(q-1)Δl,…,Xτ+1-Δl⊕uτ+1-Δl))

(7)

3.4 情境感知注意模塊

前面的模塊已經預測出了三個結果,即fi、fs、fl。在本節中,利用注意力機制,通過考慮情景信息來整合這三種預測結果,這是第一次嘗試。

在時刻tτ+1處為每個區域rn構造了情景特征向量,用gn表示。在這里,考慮了三個主要的環境因素,包括tτ+1時刻的天氣情況、tτ+1時刻的索引、rn區域的功能,這三個環境因素對出租車需求產生不同影響。外部因素的特征提取中專門采用了同樣的方法,將時間索引嵌入到一個低維向量中,并與天氣狀況向量和區域函數向量相連接。

(8)

注意,權重表明預測結果的程度大小。

由于出租車需求預測是一個回歸問題,采用均方誤差作為損失函數,并通過將預測結果與真實情況的誤差最小化來訓練網絡模型,即:

(9)

式中:Ω是網絡模型中所有可學習參數的集合。

4 實驗與結果分析

4.1 比較方法

將本文模型與以下方法進行比較。

(1) 歷史平均值(HA):通過同一時段的平均歷史出租車需求,預測每個地區下一個時段的出租車需求。

(2) 自回歸綜合移動平均值(ARIMA)[5]:是一種廣泛應用于時間序列預測問題的方法,它通過結合移動平均值和自回歸分量對時間相關性進行建模。

(3) 長短期記憶(LSTM)[7]:是循環神經網絡的一種變體,可以從序列數據中有效地學習長期和短期的潛在相關性。

(4) 擴散卷積循環神經網絡(DCRNN)[14]:將圖形卷積集成到門控循環單元中,以預測道路網絡上的交通流量。在該模型中,采用雙向圖隨機游走操作,提取交通流量的空間動態,并用RNN捕捉時間動態。

(5) 時空圖卷積網絡(STGCN)[15]:由幾個時空模塊組成,這些模塊用完全卷積層構建,用于處理車流量預測任務。具體地說,每個塊由圖像卷積和門控時間卷積組成,共同處理圖結構的時間序列。

注意,為了通過DCRNN和STGCN預測出租車需求,為這兩種方法構造了一個圖,其中每個節點代表一個區域,每條邊代表兩個區域中心之間的距離。

并對模型中不同模塊的性能進行分析,研究其在出租汽車需求預測中的有效性。

? 實時時空模塊(ISTM):只使用實時時空模塊,此模塊包括LC層和GRU。

? 短/長周期模塊和情境感知注意模塊(PM+CAAM):只使用兩個周期性模塊,同時考慮短周期和長周期。兩個模塊的輸出由情景感知注意模塊進行整合。

? 無LC的實時模塊、短/長周期模塊和情景感知注意模塊(IM+PM+CAAM):分別使用無LC層的實時時空模塊和周期性模塊預測出租車需求。三個模塊的輸出由情景感知注意模塊進行整合。

? 實時時空模塊和短/長周期模塊(ISTM+PM):只使用實時時空模塊和周期模塊,它們的輸出由網絡訓練過程中學習的權值張量進行整合。

? 基于注意力機制的卷積循環網絡(ACRN):只考慮實時時空模塊中相鄰區域之間的空間相關性,忽略了功能相似區域。網絡的其他部分與CACRNN相同。

4.2 實驗設置

實驗中參數的默認值設置如下。將時間間隔設置為15分鐘,在實時時空模塊中,使用了6個連續的時間片段,即o=6。此外,將參數默認設置為k=16、M=15、L=3和K=5。在短/長期周期模塊中,分別采用4個和2個周期時間片段,時間間隔Δs=96(一天)、Δl=96×7(一周)。將天氣狀況、假期狀況和周末分別用1×3向量表示。三個模塊中GRU中的隱藏節點數均設置為512個。在情景感知注意模塊中,將時間索引用1×5向量表示。

使用紐約市2016年6月的歷史記錄作為測試數據,其余記錄作為訓練數據。此外,采用中國大型城市成都的實際數據集進行了實驗,以進一步驗證本文方法與上文提及的基準相比的性能。成都所使用的數據集如下。

道路網絡數據集成都市道路網絡由756條道路組成。在此實驗中,采用了成都市部分區域(由以下坐標包圍: [北緯30.727 818°,東經104.043 333°],[北緯30.726 490°,東經104.129 076°],[北緯30.655 191°,東經104.129 591°],[北緯30.652 828°,東經104.042 102°]),按路段劃分為623個區域。

出租車行程數據集滴滴網發布了成都市2016年11月1日至2016年11月30日出租車出行記錄數據集,共4 141 504次。

POI數據集使用了成都地區146 587個POI數據,將其分為13個不同的類別。

氣象和節假日數據集還使用了2016年11月1日至2016年11月30日的氣象記錄數據,包含22個工作日和8個周末。

通過均方根誤差(RMSE)和平均絕對誤差(MAE)評價每種方法的性能。此外,對于所有基于DNN的模型,均采用Adam優化算法來訓練參數。Adam的學習率設置為10-4,訓練的批量大小為64。還將訓練提前停止,其中回合數和訓練次數分別設置為6和100。所有的實驗都是在NVIDIA RTX20700顯卡上進行的,實驗結果是在相同參數下使用不同的隨機種子進行5次運行計算,然后取這5次計算結果的平均值。

4.3 實驗結果

4.3.1各模型在基準下的比較

表2和表3分別顯示了本文模型和默認參數設置下的基準所實現的性能。很容易發現,與所有基準相比,本文模型在不同的數據集下達到了最低的RMSE和MAE。具體來說,HA和ARIMA表現最差。結果表明,深層神經網絡(如LSTM)可以有效地應用于城市數據預測。此外,由于LSTM只對出租車需求的時間相關性進行建模,因此LSTM的性能比本文模型差。在基準中,STGCN和DCRNN都獲得了良好的性能,它們同時捕獲了空間和時間的相關性。與DCRNN相比,模型在紐約市數據集上的RMSE(MAE)降低了9.4%(13.6%),在成都市數據集上的RMSE(MAE)降低了13.0%(8.5%)。與STGCN相比,模型在紐約市數據集上的RMSE(MAE)降低了8.9%(12.9%),在成都市數據集上的RMSE(MAE)降低了14.4%(23.6%)。

表2 性能比較——紐約市

表3 性能比較——成都市

4.3.2模塊評估

評估本文模型中不同模塊的有效性,結果分別如表4和表5所示。很容易發現,在本文模型中,每一個模塊都有更好的預測性能。特別地,通過比較ISTM和ISTM+PM的結果,確認了周期性模塊的工作有效性。由于PM+CAAC的性能比本文模型差,可以得知實時時空模塊是起到一定作用的。通過IM+PM+CAAM與本文模型的比較,驗證了LC層提取不同地區出租車需求空間相關性的有效性。很明顯,本文模型取得了比ISTM+PM更好的性能,這證實了情景感知注意模塊的有用性。此外,CACRNN考慮了功能相似區域的空間相關性,提高了ACRN的性能。

表4 模塊評估——紐約市

表5 模塊評估——成都市

4.3.3時間效率評估

分別評估了本文模型和較先進的基于DNN(即STGCN和DCRNN)所獲得的時間效率,以及訓練和測試過程的時間消耗。如表6所示,盡管DCRNN的可訓練參數數量最少,但其訓練和測試所花費的時間遠遠大于STGCN和CACRNN。此外,與STGCN相比,本文提出的CACRNN在訓練和測試過程中的時間消耗更高。這是因為CACRNN的可訓練參數幾乎是STGCN的18倍。眾所周知,在訓練過程中計算反向傳播過程中的參數梯度是非常耗時的。在實際應用中,由于DNN模型可以預先訓練,所以需要更注重測試的時間效率。幸運的是,可以觀察到CACRNN在測試過程中只比STGCN多花費17 ms,這是可以接受的。

表6 時間消耗

4.3.4參數的影響

還評估了CACRNN不同模塊中一些關鍵參數的影響,即L個LC層數、實時模塊o中的時間片段數、短周期模塊p中的時間片段數、長周期模塊q中的時間片段數。即保持其他參數不變,改變研究對象參數,通過改變參數大小觀察參數對預測精度指標的影響。在圖9中,可以發現當L從0到3變化時,MAE和RMSE減小,從而驗證了相鄰區域之間的深層空間相關性的提取。另一方面,當L從3變為5時,RMSE和MAE略有增加,因為隨著L的增加,更多不相關的區域會進行卷積運算,并對預測結果產生負面影響。 如圖10-圖12所示,當o從0變為8、p設為8、q設為4時,MAE和RMSE先減小然后增大。這表明在特定的時間段(例如90分鐘)內,出租車需求之間存在時空相關性,這有助于細粒度出租車需求的預測。此外,出租車需求具有一定的日周期性和周周期性,其過去4天和兩周的出租車需求包含了最有用的信息。

圖9 LC層數的影響

圖10 o長度的影響

圖11 p長度的影響

圖12 q長度的影響

4.3.5時空相關因素的影響

進一步評估天氣、節假日、局部區域等因素對預測精度的影響,分別將考慮天氣、節假日、局部區域的數據與不考慮上述因素的數據作為訓練數據,最后通過測試數據進行預測,相應的對比結果如表7所示。

表7 時空相關因素對預測結果的影響

可以看出,考慮了天氣、節假日與局部區域因素的預測精度明顯高于不考慮這些因素的預測精度。這進一步證明了充分考慮各個因素對預測精度的提升有較大的作用,也更符合實際情況。

5 結 語

針對出租車需求預測方法中由于數據挖掘深度不夠導致的預測精度較低問題,提出一種基于情境感知注意機制的循環卷積神經網絡出租車需求預測方法。通過紐約市和成都市的實際數據集綜合實驗結果可得如下結論:

(1) 情景感知注意模塊、功能相似區域、時空相關性等因素的引入,有利于出租車需求預測精度以及預測速度的提高。

(2) 引入短期/長期周期性模塊,使得方法對于長期短期出租車需求預測均具備較好的預測性能。

(3) 與現有的方法相比,本文方法在預測精度以及預測速度上具有較大的優勢,進一步證明了本文算法的有效性。

猜你喜歡
區域模型
一半模型
永久基本農田集中區域“禁廢”
今日農業(2021年9期)2021-11-26 07:41:24
分割區域
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
3D打印中的模型分割與打包
關于四色猜想
分區域
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
基于嚴重區域的多PCC點暫降頻次估計
電測與儀表(2015年5期)2015-04-09 11:30:52
主站蜘蛛池模板: 狠狠五月天中文字幕| 日韩国产另类| 亚洲精选高清无码| 中文字幕首页系列人妻| 久无码久无码av无码| 中文字幕久久精品波多野结| 欧美日韩免费在线视频| 亚洲精品欧美重口| 久热中文字幕在线观看| 高清无码不卡视频| 免费无码在线观看| 欧美亚洲一区二区三区导航| A级毛片高清免费视频就| 成人国产免费| 91视频首页| 欧美中文字幕无线码视频| 女人18一级毛片免费观看 | 国产精品任我爽爆在线播放6080 | 老熟妇喷水一区二区三区| 性色一区| 欧美在线三级| 国产JIZzJIzz视频全部免费| 欧美成人怡春院在线激情| 亚洲国产日韩欧美在线| 免费精品一区二区h| 这里只有精品在线| 亚洲国产一区在线观看| 欧美日韩动态图| 免费无码又爽又黄又刺激网站| 国产精品久久久久久搜索 | 亚洲swag精品自拍一区| 99视频在线免费观看| 国产精品一老牛影视频| 久热中文字幕在线观看| 亚洲经典在线中文字幕| 亚洲成a人片77777在线播放| 国产情精品嫩草影院88av| 免费看美女自慰的网站| 伊人中文网| 91精品网站| 亚洲制服丝袜第一页| 情侣午夜国产在线一区无码| 午夜福利视频一区| 亚洲精品视频免费| 欧美天天干| 在线国产毛片手机小视频| 中国成人在线视频| 重口调教一区二区视频| 国产在线拍偷自揄观看视频网站| 无码一区中文字幕| 伊人久久青草青青综合| 国产精品久久久精品三级| 国产成人高清在线精品| 青青草一区| 一级毛片不卡片免费观看| 国产无人区一区二区三区| 亚洲激情区| 91极品美女高潮叫床在线观看| 亚洲国产中文精品va在线播放| 国产在线自揄拍揄视频网站| 永久免费精品视频| 啪啪永久免费av| 国产黄色免费看| 色网站在线免费观看| 全色黄大色大片免费久久老太| 亚洲男人天堂2020| 亚洲性视频网站| 国产精品大白天新婚身材| 亚洲人在线| 少妇露出福利视频| 精品人妻一区二区三区蜜桃AⅤ| 精品无码国产一区二区三区AV| 黄色网在线| 91精品国产无线乱码在线 | 2020精品极品国产色在线观看| 高清视频一区| 亚洲综合亚洲国产尤物| 91蝌蚪视频在线观看| 国产成人综合在线观看| 99re在线视频观看| 亚洲美女AV免费一区| 国产一级α片|