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

    基于BP神經(jīng)網(wǎng)絡的船舶主機能效狀態(tài)評估

    2018-08-14 15:07:12陳偉南黃連忠張勇路通
    中國艦船研究 2018年4期
    關鍵詞:油耗航行船舶

    陳偉南,黃連忠,張勇,路通

    大連海事大學輪機工程學院,遼寧大連116026

    0 引 言

    對船舶主機而言,其能效狀態(tài)由主機的設定轉(zhuǎn)速、外界的海洋氣象狀況、船舶載重、船舶姿態(tài)等多個因素共同決定。在傳統(tǒng)的船舶管理模式中,船舶操縱人員大多依靠自身經(jīng)驗對主機狀態(tài)進行評估,并依據(jù)評估結(jié)果進行操作。這種評估具有極大的不確定性,并降低了船舶的能源利用效率[1]。

    劉伯運等[2]和孫宜權(quán)等[3]針對車用柴油機的狀態(tài)評估進行了分析。在船用柴油機方面,周根明等[4]運用層次分析法對船舶主機的整體健康狀態(tài)建立了評估體系;孫峰等[5]基于數(shù)據(jù)挖掘技術(shù)提出了一種評估船用柴油機性能的方法。這些方法雖然對船舶主機的能效狀態(tài)評估進行了研究,但并不能對航行中的船舶主機狀態(tài)進行實時、迅速的評估。葉睿等[6]利用神經(jīng)網(wǎng)絡構(gòu)建了客滾船的主機油耗模型,但其輸入?yún)?shù)較為復雜,沒有進行簡化處理,同時其輸入?yún)?shù)僅為船舶航行時的船體姿態(tài)數(shù)據(jù),并未包含主機的運行參數(shù),這也影響了模型的可靠度與準確性。

    基于此,本文將以狀態(tài)良好的船舶運行記錄作為樣本,結(jié)合主成分分析法和反向傳播(Back-Propagation,BP)神經(jīng)網(wǎng)絡智能算法,構(gòu)建船舶航行狀態(tài)識別模型和主機油耗模型。在構(gòu)建油耗模型時,將利用主成分分析法對模型輸入?yún)?shù)進行簡化,從而降低數(shù)據(jù)噪聲和網(wǎng)絡復雜性。此外,本文油耗模型的輸入?yún)?shù)還將包含主機的運行參數(shù)和船舶的航行姿態(tài)參數(shù),用以更準確地反映主機的能效狀態(tài)。本文構(gòu)建的航行狀態(tài)識別模型可用于識別船舶當前的航行狀態(tài),而油耗模型可用于判斷船舶主機當前的能耗狀態(tài)是否正常,從而對主機能效進行智能、準確的評估,可為智能船舶的能效評估研究提供一定的參考。

    1 構(gòu)建模型

    采用船舶的航行狀態(tài)模型和主機油耗模型對主機狀態(tài)進行評估的流程如圖1所示。首先,采集船舶的實時運行數(shù)據(jù),將其輸入到基于神經(jīng)網(wǎng)絡構(gòu)建的主機航行狀態(tài)識別模型中,運用該模型識別出主機不同運行數(shù)據(jù)對應的航行狀態(tài);然后,選取正常航行狀態(tài)數(shù)據(jù),運用主成分分析法對數(shù)據(jù)進行預處理,用以降低后續(xù)計算的復雜度,加快神經(jīng)網(wǎng)絡的收斂速度并減小誤差;接著,將處理后的航行狀態(tài)數(shù)據(jù)輸入到基于神經(jīng)網(wǎng)絡構(gòu)建的主機油耗模型中,得出主機的功率、轉(zhuǎn)速、油耗等參數(shù);最后,將主機的實際油耗值x2與理論油耗值x1進行實時對比,得到油耗殘差值δ,即可判斷主機當前的運行狀態(tài)是否正常。

    1.1 目標船及數(shù)據(jù)

    本文的船舶實時運行數(shù)據(jù)均來源于某目標船——30萬噸級遠洋礦砂運輸船。該船總長327 m,型寬55 m,型深29 m,設計吃水21 m,航速14.5 kn,載重 298 000 t,主機型號為 MAN B&W 6S80MC-C,額定功率為22 360 kW,額定轉(zhuǎn)速為73 r/min。

    通過船載傳感器進行數(shù)據(jù)采集工作,目標船的狀態(tài)參數(shù)如表1所示。其中,風速分解為沿船舶航行方向(X方向)和垂直于航行方向(Y方向)2個風速,并分別以船艏來風和右舷來風為正。

    表1 目標船狀態(tài)參數(shù)Table 1 State parameters of target ship

    傳感器每隔30 s采集一次數(shù)據(jù),為保證模型的穩(wěn)定性,以5 min時長為間隔來處理所采集的數(shù)據(jù)。其中,除主機油耗量數(shù)據(jù)為5 min內(nèi)的累計值之外,其他數(shù)據(jù)均為5 min內(nèi)的均值。

    X方向和Y方向的風速分別為

    式中:N為5 min內(nèi)傳感器采集的總次數(shù);Vwc為風速儀第c次采集的相對風速,其中c=1,2,…,N;θc為船舶航向與風向之間的夾角。

    主機缸套水的進口、出口溫差平均值ΔT為

    式中:Tin為主機缸套水的進口溫度;Toutg為1~6號氣缸缸套水的出口溫度,其中g(shù)=1,2,…,6。

    1~6號氣缸排氣溫度的平均值Tˉ為

    式中,Tg為1~6號氣缸的排氣溫度。

    1.2 BP神經(jīng)網(wǎng)絡

    對于輸入?yún)?shù)繁多的船舶航行狀態(tài)識別模型和油耗模型而言,若構(gòu)建一般的函數(shù),則難以準確描述各物理量之間的因果關系,且普通函數(shù)對各型船舶的通用性較差[6]。本文將采用BP神經(jīng)網(wǎng)絡來構(gòu)建模型。BP神經(jīng)網(wǎng)絡是一種非線性智能算法,具有運算速度快、容錯能力高、自學能力強等特點,是目前應用最廣泛的神經(jīng)網(wǎng)絡模型之一。BP神經(jīng)網(wǎng)絡可以學習和存儲大量的輸入—輸出模式映射關系,而無需事前揭示描述這種映射關系的數(shù)學方程[7-9]。

    BP神經(jīng)網(wǎng)絡是一種多層前饋型神經(jīng)網(wǎng)絡,其拓撲結(jié)構(gòu)如圖2所示。圖中:P1,P2,…,Pn為BP神經(jīng)網(wǎng)絡的輸入?yún)?shù);b1,b2,…,bm和θ1,θ2,…,θk為網(wǎng)絡的偏置;a1,a2,…,ak為網(wǎng)絡的預測值;ωij和ωjh為BP神經(jīng)網(wǎng)絡的權(quán)值。其中i=1,2,…,n;j=1,2,…,m;h=1,2,…,k。

    圖2中的網(wǎng)絡包括n個輸入節(jié)點和k個輸出節(jié)點,表達了從n個自變量到k個因變量的函數(shù)映射關系。網(wǎng)絡中隱含層第j個節(jié)點的輸出結(jié)果yj和輸出層第h個節(jié)點的輸出結(jié)果ah分別為

    式中,f1和f2分別為輸入關系和輸出關系的傳遞函數(shù)。

    本文將采用Matlab軟件來構(gòu)建神經(jīng)網(wǎng)絡,選取f1為雙曲線正切函數(shù)tansig,f2為線性變換函數(shù)purelin。BP神經(jīng)網(wǎng)絡的學習分為2個階段:第1個階段是根據(jù)輸入的已知樣本,從第1層向后計算各神經(jīng)元的輸出;第2個階段是從最后一層向前計算各權(quán)值和閾值對總誤差的影響,進而修改各權(quán)值和閾值。這2個階段反復交替進行,直至函數(shù)收斂。為了提高網(wǎng)絡的收斂速度,本文選用列文伯格—馬夸爾特(Levenberg-Marquardt,LM)算法作為第2個階段的計算方法,并以均方誤差函數(shù)作為性能函數(shù)。

    1.3 主成分分析法

    航行狀態(tài)識別模型的輸出結(jié)果分為5種,即停泊、機動航行、正常航行、沖車航行和大風浪航行,其設計輸入為相對風速Xˉ和Yˉ、船舶航速Vs、主機轉(zhuǎn)速ne、主機曲軸扭矩M、主機曲軸功率Pe及主機燃油消耗量F,運用BP神經(jīng)網(wǎng)絡算法可以實現(xiàn)較高準確度的航行狀態(tài)識別。

    然而,對于主機油耗模型而言,其輸出結(jié)果是具體的燃油消耗量,其輸入?yún)?shù)多達10項(表1中前10項)。在此情況下,變量之間有一定的關聯(lián)性,故變量提供給BP神經(jīng)網(wǎng)絡的信息會存在重疊,從而增加了計算的復雜性,減緩了網(wǎng)絡的收斂速度,降低了模型的準確性并提高了BP神經(jīng)網(wǎng)絡過擬合的可能性[10-11]。如果神經(jīng)網(wǎng)絡過于復雜,會使信號和噪聲同時進行擬合,從而導致神經(jīng)網(wǎng)絡的預測值與實際值嚴重不符,而采用大量的數(shù)據(jù)對網(wǎng)絡進行訓練并降低訓練數(shù)據(jù)的噪聲則是避免出現(xiàn)過擬合問題的有效方法[12]。因此,本文除了采用大量數(shù)據(jù)對模型進行訓練之外,還將應用主成分分析法(Principal Component Analysis)對油耗模型的輸入數(shù)據(jù)進行預處理以降低數(shù)據(jù)噪聲。

    主成分分析法的實質(zhì)是在不改變樣本數(shù)據(jù)結(jié)構(gòu)的情況下,通過旋轉(zhuǎn)多維空間的坐標將原變量轉(zhuǎn)換成兩兩不相關的主成分,同時盡可能反映原變量所包含的信息,從而簡化計算[13-14]。本文將采用主成分分析法與BP神經(jīng)網(wǎng)絡算法相結(jié)合的方法,假設樣本矩陣X為

    式中:樣本矩陣X中第v個參數(shù)類型中的第u個樣本數(shù)據(jù)為xuv,其中u=1,2,…,n且v=1,2,…,p;xu=[xu1xu2…xup]T,為p維列向量。

    樣本矩陣的均值xˉ為

    協(xié)方差矩陣S的表達式為

    對于n行p列的樣本矩陣X而言,經(jīng)過矩陣變換后即得到樣本協(xié)方差矩陣S,其為p行p列的方陣,故式(8)之后的u=1,2,…,p且v=1,2,…,p。

    協(xié)方差矩陣S中的任一元素Suv為

    式中:xu′u和xu′v分別為樣本矩陣X中的某個元素,其中u′=1,2,…,p;

    設Suv的p個特征值分別為且其單位正交特征向量t?為

    則第u個樣本的主成分yu為

    2 航行狀態(tài)識別模型

    對正常航行的船舶而言,在實時評估其主機能耗狀態(tài)之前,應先識別當前船舶的航行狀態(tài)。在停泊狀態(tài)下,船舶航速、主機轉(zhuǎn)速、耗油量基本為0;在機動航行狀態(tài)下,各項參數(shù)基本在停泊狀態(tài)和正常航行狀態(tài)的數(shù)值特征范圍內(nèi)波動;在正常航行狀態(tài)下,各項參數(shù)基本穩(wěn)定在一定范圍內(nèi);在沖車航行狀態(tài)下,各項參數(shù)將出現(xiàn)波動且有所增加;在大風浪航行狀態(tài)下,相對風速將明顯增加,而船舶航速將有所降低。根據(jù)不同航行狀態(tài)的參數(shù)特征,將現(xiàn)有的歷史數(shù)據(jù)進行統(tǒng)計分類,即可得到各個航行狀態(tài)的樣本數(shù)據(jù)。

    為了實現(xiàn)較高準確度的航行狀態(tài)識別,本文將以船舶航速、相對風速、主機轉(zhuǎn)速、曲軸扭矩、曲軸功率和主機油耗量作為輸入,船舶航行狀態(tài)作為輸出。同時,還需要設定神經(jīng)網(wǎng)絡計算矩陣的目標值,即停泊為[1 0 0 0 0]T、機動航行為[0 1 0 0 0]T、正常航行為[0 0 1 0 0]T、沖車航行為[0 0 0 1 0]T、大風浪航行為[0 0 0 0 1]T,然后進行訓練學習。在Matlab軟件中建立并訓練BP神經(jīng)網(wǎng)絡時,設定訓練函數(shù)為trainlm、學習函數(shù)為learngdm、性能函數(shù)為mse、隱含層和輸出層的傳遞函數(shù)分別為tansig和purelin,其中2層隱含層的神經(jīng)元數(shù)量均為10。訓練函數(shù)trainlm采用LM算法,該方法結(jié)合了梯度下降法和牛頓法的優(yōu)點,可以有效避免陷入局部最優(yōu)解并保證收斂速度[15]。

    本文將選取目標船在新加坡—圣路易斯航段(2015年03月20日至2015年04月06日)作為訓練樣本,即該時間段內(nèi)船舶和主機均處于良好狀態(tài)下的航行數(shù)據(jù),其中停泊記錄22條、機動航行記錄33條、正常航行記錄153條、沖車航行記錄14條、大風浪航行記錄114條。

    為了驗證模型的準確性,本文選擇了4個月后的航行記錄模型作為檢驗樣本,即目標船在馬跡山—新加坡航段(2015年07月19日至2015年07月22日)內(nèi)連續(xù)77 h的樣本數(shù)據(jù),仿真結(jié)果如圖3所示。圖中,各分圖的橫坐標均為傳感器數(shù)據(jù)采集點的序號。

    圖3(g)所示的評估航行狀態(tài)圖即模型識別所得的船舶航行狀態(tài),圖3(h)所示的實際航行狀態(tài)圖即通過查詢該目標船輪機日志和航海日志后統(tǒng)計得出的船舶實際航行狀態(tài),其縱坐標數(shù)值1,2,3,4,5分別代表停泊、機動航行、正常航行、沖車航行、大風浪航行這5個狀態(tài)。由圖3可知,在第0~120采集點之間,船舶航速、主機轉(zhuǎn)速、扭矩、功率、油耗等均為0,識別模型判斷為停泊狀態(tài),判斷準確;在第120~135采集點之間,船舶和主機的各項參數(shù)均出現(xiàn)了波動,模型將其識別為機動航行狀態(tài),判斷準確;在第128點左右,由于船舶的短暫提速,模型將其識別為正常航行狀態(tài),判斷錯誤;在第135~486點之間,模型判斷船舶為正常航行狀態(tài),經(jīng)查詢航海日志,該目標船確實為正常航行狀態(tài),判斷準確;在第486~719點之間,風速提高、船速降低,模型判斷為大風浪航行狀態(tài),經(jīng)查詢航海日志,判斷準確;在第885~892點之間,各項參數(shù)均出現(xiàn)了波動且有所增加,模型判斷為沖車航行狀態(tài),經(jīng)查詢輪機日志,該目標船在此時間段確實進行了沖車操作,故模型判斷準確。由此可知,模型判斷整個樣本的識別正確率為98.05%,從而為下一步建立主機油耗模型提供了基礎。

    3 主機油耗模型

    3.1 主成分分析處理

    應用BP神經(jīng)網(wǎng)絡構(gòu)建主機油耗模型時,需要先對訓練樣本進行主成分分析。本文選取了目標船在2015年3~5月的船舶和主機運行數(shù)據(jù)作為訓練樣本,該時間段內(nèi)目標船和主機均處于良好運行狀態(tài)。選取表1的前10個參數(shù)作為主成分分析的原始變量,設定主成分Z1~Z10,計算所得的成分矩陣如表2所示,各主成分對應的特征值、方差貢獻率及累計方差貢獻率如表3所示。

    由表3可知,Z1~Z5的累計貢獻率為93.87%,而Z6~Z10的貢獻率則很小。因此,本文將選擇Z1~Z5作為表述樣本內(nèi)容的主要數(shù)據(jù)來計算各項主成分的系數(shù),結(jié)果如表4所示。根據(jù)各主成分的系數(shù),將原訓練樣本的10種原始變量轉(zhuǎn)換為5項主成分,即可作為BP神經(jīng)網(wǎng)絡的輸入數(shù)據(jù)。

    3.2 建立油耗模型

    在matlab中設定BP網(wǎng)絡油耗模型的訓練函數(shù)為trainlm、學習函數(shù)為learngdm、性能函數(shù)為mse、隱含層和輸出層的傳遞函數(shù)分別為tansig和purelin。對于隱含層的層數(shù)和神經(jīng)元節(jié)點數(shù)目,目前缺乏設定最優(yōu)數(shù)目的明確指導方法,本文將設定神經(jīng)元節(jié)點數(shù)目的取值范圍為0~20,然后采用循環(huán)程序遍歷2層隱含層中所有神經(jīng)元數(shù)目組合的模型訓練誤差,最后根據(jù)計算所得的訓練誤差值確定最優(yōu)的神經(jīng)元數(shù)目設定。如圖4所示,將3.1節(jié)中已處理的主成分數(shù)據(jù)信息作為BP網(wǎng)絡的設計輸入,通過循環(huán)程序即可得到不同神經(jīng)元數(shù)目組合的訓練誤差結(jié)果。

    表2 成分矩陣Table 2 Component matrix

    表3 總方差解釋Table 3 Total variance interpretation

    表4 各項主成分的系數(shù)Table 4 The coefficients of each principal component

    經(jīng)計算,當?shù)?層隱含層神經(jīng)元數(shù)目為16且第2層為11時,數(shù)據(jù)整體的均方誤差最小,其值為6.325 9×10-5。因此,本文設定BP神經(jīng)網(wǎng)絡中2層隱含層的神經(jīng)元數(shù)目分別為16和11,然后,即可建立并訓練網(wǎng)絡,最終得出正常航行狀態(tài)下的主機油耗模型。

    需要指出的是,由于不同的船舶和主機的匹配關系與性能均有所差別,所以針對其他船舶構(gòu)建航行狀態(tài)識別模型和主機油耗模型時,還需要重新通過計算機程序進行神經(jīng)網(wǎng)絡訓練,所需的時間成本也應在可接受范圍內(nèi)。本文已經(jīng)確定了神經(jīng)網(wǎng)絡訓練的模式和函數(shù),故可以推廣應用至其他船舶和主機。

    3.3 模型驗證與分析

    選取目標船在圣路易斯—馬跡山—新加坡航段(2015年05月19日至2015年05月24日)的部分航程作為驗證樣本,經(jīng)航行狀態(tài)識別模型提取正常航行狀態(tài)數(shù)據(jù),并根據(jù)各主成分系數(shù)計算得出5項主成分參數(shù),輸入BP神經(jīng)網(wǎng)絡即可得出對應的理論油耗值,結(jié)果如圖5所示。

    由于主機燃油流量計數(shù)器的輸出數(shù)據(jù)精度為0.01 m3,故圖5中的實際油耗值曲線存在幅值為0.01左右的波動,但總體趨勢與計算油耗值曲線一致。在第150個數(shù)據(jù)點附近,計算油耗值出現(xiàn)了小幅降低,船舶的曲軸扭矩和功率分別在1 620 kN·m和9 600 kW左右波動,略低于整體平均水平1 673 kN·m和10 022 kW,而實際油耗值也出現(xiàn)了降低現(xiàn)象,可見計算結(jié)果與實際運行情況相符。在第500個數(shù)據(jù)點附近,實際油耗值出現(xiàn)了異常高峰點,經(jīng)查為傳感器的暫時故障所致,而其他數(shù)據(jù)點的傳感器工作正常。由圖5可知,主機油耗模型的總體平均相對誤差為3.47%,且未產(chǎn)生明顯的過擬合現(xiàn)象。

    將實際油耗值和計算油耗值進行對比,以殘差形式體現(xiàn)主機的能效狀態(tài),如圖6所示。

    由圖6可知,在該段航程內(nèi)的各項工況下,主機的性能特征與良好能效狀態(tài)相吻合,即主機能效狀態(tài)正常。同時,由于主機燃油流量計數(shù)器的輸出數(shù)據(jù)精度為0.01 m3,故大部分樣本數(shù)據(jù)點的耗油量殘差都在-0.01~0.01范圍內(nèi)波動。

    4 結(jié) 語

    本文結(jié)合實船航行數(shù)據(jù),運用BP神經(jīng)網(wǎng)絡算法建立了船舶航行狀態(tài)識別模型,可以有效識別船舶的5種航行狀態(tài),其準確率達98.05%。同時,依據(jù)船舶和主機的10個運行參數(shù)與主機燃油消耗量之間的關系,采用主成分分析法降低了各個變量之間的關聯(lián)性和復雜性,并運用BP神經(jīng)網(wǎng)絡算法構(gòu)建了船舶主機正常值油耗模型,其總體平均相對誤差為3.47%。通過綜合運用這2種模型,可以對正常航行狀態(tài)下的主機實時狀態(tài)進行評估,可為智能船舶的能效管理提供參考。

    猜你喜歡
    油耗航行船舶
    不談油耗 只講運動 試駕第十一代思域e:HEV
    車主之友(2022年5期)2022-11-23 07:22:20
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    到慧骃國的航行
    船舶!請加速
    BOG壓縮機在小型LNG船舶上的應用
    降低內(nèi)燃裝卸機械油耗措施的探討
    雙管齊下 YarisL致享綜合油耗測試
    車迷(2017年12期)2018-01-18 02:16:10
    小舟在河上航行
    中學生英語(2017年6期)2017-07-31 21:28:55
    船舶壓載水管理系統(tǒng)
    中國船檢(2017年3期)2017-05-18 11:33:09
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    人妻久久中文字幕网| 亚洲国产精品合色在线| 国产深夜福利视频在线观看| 国产黄a三级三级三级人| 一级片免费观看大全| aaaaa片日本免费| 精品国产乱子伦一区二区三区| 777久久人妻少妇嫩草av网站| 999久久久国产精品视频| 国产精品免费一区二区三区在线| 一级片免费观看大全| 亚洲成国产人片在线观看| 看片在线看免费视频| 一级黄色大片毛片| 1024香蕉在线观看| 午夜福利在线观看吧| 大码成人一级视频| 一本综合久久免费| 新久久久久国产一级毛片| 一二三四社区在线视频社区8| 又大又爽又粗| www.精华液| 国产免费男女视频| 女人高潮潮喷娇喘18禁视频| 午夜福利免费观看在线| 国产精品秋霞免费鲁丝片| xxxhd国产人妻xxx| 妹子高潮喷水视频| 天堂√8在线中文| 在线观看免费日韩欧美大片| 夜夜躁狠狠躁天天躁| 中文亚洲av片在线观看爽| 日韩欧美一区视频在线观看| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| 日韩中文字幕欧美一区二区| 精品国产美女av久久久久小说| 大型av网站在线播放| 日本vs欧美在线观看视频| 国产精品电影一区二区三区| 好男人电影高清在线观看| 国产精品二区激情视频| 久久午夜亚洲精品久久| 久久午夜综合久久蜜桃| 国产成人欧美| 女生性感内裤真人,穿戴方法视频| 午夜免费激情av| √禁漫天堂资源中文www| 色婷婷久久久亚洲欧美| 国产精品香港三级国产av潘金莲| 亚洲国产欧美一区二区综合| 欧美大码av| 国产精品98久久久久久宅男小说| 日韩欧美一区视频在线观看| 嫁个100分男人电影在线观看| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 精品国产一区二区三区四区第35| 欧美成人性av电影在线观看| 桃色一区二区三区在线观看| 无限看片的www在线观看| 久久久久亚洲av毛片大全| 在线天堂中文资源库| 久久性视频一级片| 欧美久久黑人一区二区| 一进一出抽搐gif免费好疼 | 老司机深夜福利视频在线观看| 男女午夜视频在线观看| 久久中文字幕一级| 亚洲欧美精品综合久久99| 欧美激情久久久久久爽电影 | 99国产精品一区二区蜜桃av| 亚洲avbb在线观看| 丰满饥渴人妻一区二区三| 人人妻人人爽人人添夜夜欢视频| 无人区码免费观看不卡| 欧美乱色亚洲激情| 在线观看一区二区三区| 一级,二级,三级黄色视频| 久久久久精品国产欧美久久久| 变态另类成人亚洲欧美熟女 | 人成视频在线观看免费观看| 欧美黑人精品巨大| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 999久久久国产精品视频| 色综合站精品国产| 国产精品亚洲av一区麻豆| 在线观看舔阴道视频| 91成人精品电影| 久久人妻熟女aⅴ| 首页视频小说图片口味搜索| 又紧又爽又黄一区二区| 国产一区二区激情短视频| 精品国产一区二区久久| 久久久久久久久中文| 纯流量卡能插随身wifi吗| www.999成人在线观看| 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 亚洲国产欧美日韩在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 香蕉丝袜av| 国产单亲对白刺激| 亚洲人成网站在线播放欧美日韩| videosex国产| 成熟少妇高潮喷水视频| 激情视频va一区二区三区| 国产91精品成人一区二区三区| 制服诱惑二区| 国产成+人综合+亚洲专区| 国产在线观看jvid| 大香蕉久久成人网| 可以在线观看毛片的网站| 伊人久久大香线蕉亚洲五| 黄色视频,在线免费观看| 首页视频小说图片口味搜索| 日韩中文字幕欧美一区二区| 久久中文字幕人妻熟女| 久久久久国产一级毛片高清牌| 亚洲一区二区三区色噜噜 | 欧美日本亚洲视频在线播放| 午夜免费鲁丝| 免费一级毛片在线播放高清视频 | 亚洲精品av麻豆狂野| 亚洲av片天天在线观看| 久久久久国产精品人妻aⅴ院| 黄色 视频免费看| 乱人伦中国视频| 日韩欧美一区视频在线观看| 美女高潮到喷水免费观看| cao死你这个sao货| 国产成人av教育| 99久久综合精品五月天人人| 免费在线观看日本一区| 一级a爱视频在线免费观看| 嫩草影视91久久| 成年人免费黄色播放视频| 亚洲欧美激情综合另类| 亚洲成人免费电影在线观看| 黄色视频,在线免费观看| 久久精品国产99精品国产亚洲性色 | 亚洲男人的天堂狠狠| 又紧又爽又黄一区二区| 99国产综合亚洲精品| 老司机午夜福利在线观看视频| 亚洲国产精品合色在线| 亚洲第一青青草原| 女人精品久久久久毛片| 久久精品91蜜桃| 免费高清视频大片| 成人亚洲精品av一区二区 | 欧美日韩福利视频一区二区| 一本综合久久免费| 欧美午夜高清在线| 亚洲成a人片在线一区二区| 色老头精品视频在线观看| 国产91精品成人一区二区三区| 免费在线观看黄色视频的| 欧美乱码精品一区二区三区| 一进一出好大好爽视频| 中文字幕人妻熟女乱码| 欧美不卡视频在线免费观看 | 久久精品国产综合久久久| 日本欧美视频一区| 久9热在线精品视频| 精品国产国语对白av| 国产xxxxx性猛交| 高清毛片免费观看视频网站 | 精品乱码久久久久久99久播| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| 国产精品爽爽va在线观看网站 | 欧美在线一区亚洲| 久久国产亚洲av麻豆专区| 丁香欧美五月| 人人妻人人澡人人看| 老司机深夜福利视频在线观看| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 亚洲自偷自拍图片 自拍| 久久中文看片网| bbb黄色大片| 热99国产精品久久久久久7| 他把我摸到了高潮在线观看| 极品人妻少妇av视频| 一边摸一边抽搐一进一出视频| 亚洲男人天堂网一区| 日韩一卡2卡3卡4卡2021年| 色老头精品视频在线观看| 一二三四社区在线视频社区8| 亚洲成国产人片在线观看| 亚洲中文字幕日韩| 90打野战视频偷拍视频| 久久久久久久久免费视频了| 日韩精品免费视频一区二区三区| 日本wwww免费看| 777久久人妻少妇嫩草av网站| 成年人免费黄色播放视频| 亚洲欧美精品综合久久99| 不卡一级毛片| 久久亚洲真实| 久久精品国产综合久久久| 在线永久观看黄色视频| 热99re8久久精品国产| 免费日韩欧美在线观看| 午夜福利在线免费观看网站| 高清在线国产一区| 美女大奶头视频| 久久久久九九精品影院| 无限看片的www在线观看| 美女高潮喷水抽搐中文字幕| 久久久国产一区二区| 激情视频va一区二区三区| 亚洲欧美精品综合久久99| 999精品在线视频| 亚洲七黄色美女视频| 最近最新中文字幕大全电影3 | 我的亚洲天堂| 精品人妻在线不人妻| 午夜老司机福利片| 国产成人欧美| 亚洲人成电影观看| 女性被躁到高潮视频| www.999成人在线观看| 精品午夜福利视频在线观看一区| 国产高清视频在线播放一区| 最近最新免费中文字幕在线| 久久精品亚洲熟妇少妇任你| 国产精品日韩av在线免费观看 | 日日摸夜夜添夜夜添小说| av天堂在线播放| 露出奶头的视频| 99久久国产精品久久久| 看免费av毛片| 精品一区二区三区视频在线观看免费 | 欧美一区二区精品小视频在线| 国产av一区二区精品久久| 欧美成人性av电影在线观看| 亚洲男人的天堂狠狠| 欧美日韩瑟瑟在线播放| 亚洲自拍偷在线| 啦啦啦 在线观看视频| 日本wwww免费看| 午夜福利影视在线免费观看| 亚洲自拍偷在线| 一二三四社区在线视频社区8| 欧美日韩福利视频一区二区| 亚洲中文日韩欧美视频| 久久久久久人人人人人| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 一区二区三区国产精品乱码| av中文乱码字幕在线| 国产成人精品在线电影| 国产区一区二久久| 不卡一级毛片| 亚洲人成电影观看| 在线观看66精品国产| 久久精品人人爽人人爽视色| av国产精品久久久久影院| 久久精品aⅴ一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区视频在线观看免费 | 欧美黄色淫秽网站| 国产精品影院久久| 国产成人av教育| 精品国产超薄肉色丝袜足j| 亚洲欧美一区二区三区黑人| 国产成+人综合+亚洲专区| 黄色视频不卡| 高潮久久久久久久久久久不卡| 在线观看免费视频网站a站| 久久这里只有精品19| 精品国产亚洲在线| 18禁裸乳无遮挡免费网站照片 | 成人三级做爰电影| 成人三级黄色视频| xxx96com| 日韩欧美一区视频在线观看| 成人三级做爰电影| 精品人妻1区二区| 午夜精品久久久久久毛片777| 麻豆成人av在线观看| 日本免费a在线| 免费一级毛片在线播放高清视频 | 两性夫妻黄色片| 欧美日韩av久久| 丝袜美腿诱惑在线| 99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| 国产片内射在线| 午夜福利欧美成人| 精品欧美一区二区三区在线| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 999久久久国产精品视频| 国产精品成人在线| 黄色成人免费大全| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 91av网站免费观看| 又黄又粗又硬又大视频| 国产精品爽爽va在线观看网站 | 日韩有码中文字幕| 国产亚洲精品久久久久久毛片| 欧美在线一区亚洲| 欧美不卡视频在线免费观看 | 亚洲少妇的诱惑av| 国产精品久久久久久人妻精品电影| 免费久久久久久久精品成人欧美视频| 两性午夜刺激爽爽歪歪视频在线观看 | 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久av美女十八| 69av精品久久久久久| 男女高潮啪啪啪动态图| 老汉色∧v一级毛片| 久久久久久大精品| 色综合欧美亚洲国产小说| 亚洲欧美激情在线| 国产精品九九99| 欧美中文日本在线观看视频| 久热爱精品视频在线9| 99热只有精品国产| 老司机午夜十八禁免费视频| 伊人久久大香线蕉亚洲五| 波多野结衣一区麻豆| 黄网站色视频无遮挡免费观看| 日韩av在线大香蕉| 韩国av一区二区三区四区| 美女高潮到喷水免费观看| 黄色怎么调成土黄色| 日韩欧美三级三区| 久久热在线av| 精品乱码久久久久久99久播| 91老司机精品| 黄片小视频在线播放| svipshipincom国产片| 久久国产精品男人的天堂亚洲| 女人精品久久久久毛片| 久久久国产欧美日韩av| 亚洲第一欧美日韩一区二区三区| 国产xxxxx性猛交| www国产在线视频色| av中文乱码字幕在线| 日韩有码中文字幕| 午夜视频精品福利| 真人一进一出gif抽搐免费| 一边摸一边做爽爽视频免费| a在线观看视频网站| 麻豆成人av在线观看| 亚洲五月天丁香| 久久久久国产精品人妻aⅴ院| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 亚洲人成77777在线视频| 丁香六月欧美| 久久久久久久午夜电影 | 精品第一国产精品| 精品国产美女av久久久久小说| 国产精品 欧美亚洲| 午夜福利欧美成人| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 午夜老司机福利片| 乱人伦中国视频| 极品人妻少妇av视频| 神马国产精品三级电影在线观看 | 日韩成人在线观看一区二区三区| 大码成人一级视频| 久久亚洲精品不卡| 黄色丝袜av网址大全| 亚洲激情在线av| 国产人伦9x9x在线观看| 黄片小视频在线播放| 亚洲精品国产色婷婷电影| 国产av一区二区精品久久| 午夜久久久在线观看| 午夜a级毛片| 精品久久久久久成人av| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级| 亚洲精品美女久久久久99蜜臀| 久久青草综合色| 亚洲av美国av| 看免费av毛片| 亚洲人成77777在线视频| 亚洲欧美一区二区三区久久| 男女之事视频高清在线观看| 又黄又爽又免费观看的视频| 啪啪无遮挡十八禁网站| 午夜激情av网站| 麻豆av在线久日| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩一区二区精品| 久久欧美精品欧美久久欧美| 亚洲免费av在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 法律面前人人平等表现在哪些方面| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 中文字幕人妻丝袜一区二区| 亚洲中文日韩欧美视频| 欧美日韩亚洲综合一区二区三区_| 中文字幕精品免费在线观看视频| 国产成人精品无人区| 手机成人av网站| 免费搜索国产男女视频| 国产精品亚洲一级av第二区| 高清毛片免费观看视频网站 | 亚洲黑人精品在线| 少妇 在线观看| 亚洲五月婷婷丁香| 男人舔女人的私密视频| 性欧美人与动物交配| 国产野战对白在线观看| 久久久久久亚洲精品国产蜜桃av| 天天影视国产精品| 人妻久久中文字幕网| 久热这里只有精品99| 国产高清videossex| 国产欧美日韩一区二区三| 国产av一区在线观看免费| 岛国视频午夜一区免费看| 黄片播放在线免费| 亚洲av电影在线进入| 亚洲自偷自拍图片 自拍| 国产精品久久久人人做人人爽| 亚洲欧美精品综合久久99| 在线观看免费视频网站a站| 国产人伦9x9x在线观看| 国产黄a三级三级三级人| 色播在线永久视频| 高清黄色对白视频在线免费看| 亚洲一区中文字幕在线| 亚洲人成电影观看| 欧美一区二区精品小视频在线| 精品一区二区三区四区五区乱码| 久久久国产欧美日韩av| www.熟女人妻精品国产| 啦啦啦 在线观看视频| 91在线观看av| 69av精品久久久久久| 亚洲欧美激情在线| bbb黄色大片| 日韩成人在线观看一区二区三区| 日韩欧美在线二视频| 久久中文看片网| 中文字幕av电影在线播放| 久久中文字幕人妻熟女| 亚洲精品一卡2卡三卡4卡5卡| 久久性视频一级片| 久久婷婷成人综合色麻豆| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 精品午夜福利视频在线观看一区| a级片在线免费高清观看视频| 久久青草综合色| av在线天堂中文字幕 | 日韩欧美三级三区| 国产成人精品久久二区二区免费| 天堂中文最新版在线下载| 欧美不卡视频在线免费观看 | 精品久久久久久,| 午夜成年电影在线免费观看| 中文字幕人妻丝袜一区二区| 欧美日韩国产mv在线观看视频| 男女下面插进去视频免费观看| 日韩中文字幕欧美一区二区| 午夜老司机福利片| 久久狼人影院| 免费在线观看日本一区| 午夜亚洲福利在线播放| 少妇 在线观看| 亚洲人成电影免费在线| 国产又爽黄色视频| 91老司机精品| 欧美激情高清一区二区三区| 亚洲人成网站在线播放欧美日韩| 三级毛片av免费| 香蕉丝袜av| 国产成+人综合+亚洲专区| 午夜福利在线观看吧| 老熟妇仑乱视频hdxx| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 日韩一卡2卡3卡4卡2021年| 嫩草影视91久久| 午夜精品在线福利| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜制服| 热re99久久国产66热| 美女高潮到喷水免费观看| 怎么达到女性高潮| 亚洲成av片中文字幕在线观看| 制服人妻中文乱码| avwww免费| 在线av久久热| 婷婷精品国产亚洲av在线| x7x7x7水蜜桃| 国产高清视频在线播放一区| 一区在线观看完整版| 午夜两性在线视频| 久久久久九九精品影院| 女同久久另类99精品国产91| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 国产精品 国内视频| 九色亚洲精品在线播放| 一级片免费观看大全| 777久久人妻少妇嫩草av网站| 母亲3免费完整高清在线观看| 老熟妇乱子伦视频在线观看| 母亲3免费完整高清在线观看| 九色亚洲精品在线播放| 国产熟女午夜一区二区三区| 午夜影院日韩av| 制服人妻中文乱码| 99国产综合亚洲精品| 国产精品亚洲一级av第二区| 日本撒尿小便嘘嘘汇集6| 久久狼人影院| 国产成人免费无遮挡视频| 欧美在线一区亚洲| 午夜精品国产一区二区电影| 成人亚洲精品av一区二区 | 亚洲av片天天在线观看| 在线观看午夜福利视频| 国产av在哪里看| 黄色丝袜av网址大全| 露出奶头的视频| 久久国产精品影院| 日本vs欧美在线观看视频| 午夜免费观看网址| 国产亚洲精品一区二区www| 午夜免费激情av| 久久人人精品亚洲av| 成人18禁高潮啪啪吃奶动态图| 久久久精品欧美日韩精品| 国产精品1区2区在线观看.| 三级毛片av免费| 国产免费av片在线观看野外av| 男男h啪啪无遮挡| 18美女黄网站色大片免费观看| 国产精品影院久久| 最新美女视频免费是黄的| av在线天堂中文字幕 | 国产精品 国内视频| 成人国产一区最新在线观看| 欧洲精品卡2卡3卡4卡5卡区| 露出奶头的视频| 女同久久另类99精品国产91| x7x7x7水蜜桃| 精品乱码久久久久久99久播| 又黄又粗又硬又大视频| 久9热在线精品视频| 一a级毛片在线观看| av国产精品久久久久影院| 亚洲一码二码三码区别大吗| 身体一侧抽搐| 国产成人欧美在线观看| 757午夜福利合集在线观看| 亚洲色图 男人天堂 中文字幕| 巨乳人妻的诱惑在线观看| 日日爽夜夜爽网站| 18美女黄网站色大片免费观看| 国产深夜福利视频在线观看| 高清在线国产一区| 国产精品av久久久久免费| 丰满饥渴人妻一区二区三| 久久久久国产一级毛片高清牌| 日韩三级视频一区二区三区| 国产精品一区二区精品视频观看| 亚洲欧美激情综合另类| 美女高潮喷水抽搐中文字幕| 色婷婷久久久亚洲欧美| 一进一出抽搐动态| 在线观看免费视频网站a站| 亚洲国产欧美网| 成人手机av| 99re在线观看精品视频| 一区二区三区国产精品乱码| 99riav亚洲国产免费| 亚洲国产毛片av蜜桃av| 窝窝影院91人妻| 免费看十八禁软件| 欧美av亚洲av综合av国产av| 日日夜夜操网爽| 1024视频免费在线观看| 亚洲成人免费av在线播放| 免费观看精品视频网站| 精品免费久久久久久久清纯| 欧美一区二区精品小视频在线| 高清在线国产一区| 欧美激情高清一区二区三区| 中文亚洲av片在线观看爽| 99久久综合精品五月天人人| 亚洲av成人av| 变态另类成人亚洲欧美熟女 | 老司机亚洲免费影院| 国产成人av教育| 亚洲av美国av| 久久久久国产精品人妻aⅴ院| 国产极品粉嫩免费观看在线| 宅男免费午夜| 亚洲一区中文字幕在线| 极品人妻少妇av视频| 亚洲熟妇中文字幕五十中出 | 久久九九热精品免费| avwww免费| 欧美乱码精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩欧美一区二区三区在线观看| 无限看片的www在线观看| 免费在线观看黄色视频的|