• <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
    欧美激情久久久久久爽电影| 啦啦啦免费观看视频1| 亚洲av成人av| 黄片大片在线免费观看| 级片在线观看| 午夜两性在线视频| 男女做爰动态图高潮gif福利片| 少妇裸体淫交视频免费看高清| 国产美女午夜福利| 超碰av人人做人人爽久久 | 国产精品 欧美亚洲| 日本熟妇午夜| 亚洲美女黄片视频| 淫秽高清视频在线观看| 在线看三级毛片| 国产精品野战在线观看| 免费在线观看影片大全网站| 亚洲av美国av| 男人的好看免费观看在线视频| 黄色丝袜av网址大全| 国产毛片a区久久久久| 村上凉子中文字幕在线| 日韩欧美免费精品| 色视频www国产| 99久久99久久久精品蜜桃| 午夜福利高清视频| 亚洲无线在线观看| 99热6这里只有精品| 国内精品久久久久精免费| 三级男女做爰猛烈吃奶摸视频| 一级a爱片免费观看的视频| 深爱激情五月婷婷| 国产淫片久久久久久久久 | 国产成人a区在线观看| 午夜免费男女啪啪视频观看 | 男女做爰动态图高潮gif福利片| 99在线人妻在线中文字幕| 亚洲av第一区精品v没综合| 又粗又爽又猛毛片免费看| 欧美日韩综合久久久久久 | 亚洲,欧美精品.| 亚洲真实伦在线观看| 国产野战对白在线观看| 九九久久精品国产亚洲av麻豆| 久久香蕉国产精品| 欧美激情在线99| 麻豆国产av国片精品| 亚洲成a人片在线一区二区| 亚洲人成网站在线播| 高潮久久久久久久久久久不卡| 好看av亚洲va欧美ⅴa在| av福利片在线观看| 欧美最新免费一区二区三区 | 午夜视频国产福利| 日本五十路高清| 久久久国产精品麻豆| 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 国产99白浆流出| 国产高清激情床上av| 日本a在线网址| 性欧美人与动物交配| 亚洲成av人片在线播放无| 嫩草影视91久久| 国产精品日韩av在线免费观看| 亚洲人成电影免费在线| 韩国av一区二区三区四区| 神马国产精品三级电影在线观看| 香蕉丝袜av| 欧美丝袜亚洲另类 | 桃色一区二区三区在线观看| 久久久久国产精品人妻aⅴ院| 久久国产乱子伦精品免费另类| 欧美一区二区亚洲| 亚洲avbb在线观看| 色老头精品视频在线观看| 日韩成人在线观看一区二区三区| 亚洲av第一区精品v没综合| 国产精品久久久久久久久免 | 亚洲av五月六月丁香网| 黄色女人牲交| 青草久久国产| 日日夜夜操网爽| 91字幕亚洲| 欧美性猛交黑人性爽| 国产成人aa在线观看| 免费一级毛片在线播放高清视频| 三级国产精品欧美在线观看| 琪琪午夜伦伦电影理论片6080| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 舔av片在线| 老司机午夜十八禁免费视频| 亚洲av五月六月丁香网| 九色成人免费人妻av| 99热6这里只有精品| 国产午夜精品论理片| 日韩欧美精品v在线| 日本免费a在线| 欧美3d第一页| 国产毛片a区久久久久| 亚洲美女视频黄频| 夜夜爽天天搞| 99热这里只有是精品50| 91字幕亚洲| 舔av片在线| 性欧美人与动物交配| 免费看光身美女| 国产一区在线观看成人免费| 久久天躁狠狠躁夜夜2o2o| 高清在线国产一区| 亚洲专区中文字幕在线| 日韩av在线大香蕉| 99热这里只有是精品50| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 99精品在免费线老司机午夜| 18禁黄网站禁片免费观看直播| 亚洲国产精品成人综合色| 99久久成人亚洲精品观看| 国产99白浆流出| 国产亚洲欧美在线一区二区| 老司机福利观看| 国产老妇女一区| 长腿黑丝高跟| 精品日产1卡2卡| 久久国产精品影院| 久久久精品欧美日韩精品| 亚洲精品456在线播放app | 两性午夜刺激爽爽歪歪视频在线观看| 丰满人妻一区二区三区视频av | 女同久久另类99精品国产91| 国产一区二区亚洲精品在线观看| 色av中文字幕| 午夜福利高清视频| 两个人视频免费观看高清| 国产高清videossex| av在线蜜桃| 中文字幕熟女人妻在线| 特大巨黑吊av在线直播| av在线蜜桃| 极品教师在线免费播放| 国产精品 欧美亚洲| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 99久久综合精品五月天人人| 国产免费一级a男人的天堂| 天天添夜夜摸| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品人妻偷拍中文字幕| 老汉色av国产亚洲站长工具| a在线观看视频网站| 欧美乱妇无乱码| 成年女人毛片免费观看观看9| 美女被艹到高潮喷水动态| 欧美不卡视频在线免费观看| 亚洲精品在线观看二区| АⅤ资源中文在线天堂| 中亚洲国语对白在线视频| 蜜桃亚洲精品一区二区三区| 日韩欧美在线乱码| 噜噜噜噜噜久久久久久91| 俺也久久电影网| 在线播放国产精品三级| 欧美不卡视频在线免费观看| 一级a爱片免费观看的视频| 性色avwww在线观看| a级一级毛片免费在线观看| 精品久久久久久,| 国产精品久久视频播放| 欧美日韩国产亚洲二区| 欧美午夜高清在线| 欧美丝袜亚洲另类 | 51国产日韩欧美| 亚洲国产精品合色在线| 麻豆国产97在线/欧美| 最新在线观看一区二区三区| 久久精品91蜜桃| 九九热线精品视视频播放| 成熟少妇高潮喷水视频| 日韩精品青青久久久久久| 日韩av在线大香蕉| 亚洲中文日韩欧美视频| 他把我摸到了高潮在线观看| 非洲黑人性xxxx精品又粗又长| 激情在线观看视频在线高清| 中文字幕熟女人妻在线| 亚洲在线观看片| 欧美av亚洲av综合av国产av| 日韩亚洲欧美综合| 手机成人av网站| 亚洲精品色激情综合| 老司机午夜十八禁免费视频| 动漫黄色视频在线观看| 亚洲国产精品成人综合色| 啦啦啦韩国在线观看视频| 搡老岳熟女国产| 久久午夜亚洲精品久久| 国产乱人伦免费视频| 麻豆国产97在线/欧美| 午夜福利成人在线免费观看| 此物有八面人人有两片| 亚洲五月婷婷丁香| 少妇高潮的动态图| 欧美性感艳星| 丁香六月欧美| 色吧在线观看| 精品久久久久久成人av| 日本五十路高清| 十八禁人妻一区二区| 国产午夜精品论理片| 亚洲,欧美精品.| 女人被狂操c到高潮| 动漫黄色视频在线观看| 天天躁日日操中文字幕| 国产伦在线观看视频一区| 久久中文看片网| 在线看三级毛片| 久久香蕉精品热| 亚洲性夜色夜夜综合| 九色国产91popny在线| 国产真实伦视频高清在线观看 | 亚洲黑人精品在线| 国产亚洲精品综合一区在线观看| 欧美午夜高清在线| 精品午夜福利视频在线观看一区| 国产三级中文精品| 精品欧美国产一区二区三| 99久久99久久久精品蜜桃| 99久久久亚洲精品蜜臀av| 3wmmmm亚洲av在线观看| 成年女人永久免费观看视频| 美女cb高潮喷水在线观看| 久久久国产成人精品二区| 国产精品亚洲美女久久久| 少妇人妻一区二区三区视频| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 亚洲欧美日韩东京热| 国产真人三级小视频在线观看| 两个人视频免费观看高清| 熟女少妇亚洲综合色aaa.| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 免费看十八禁软件| 九色国产91popny在线| 国产免费一级a男人的天堂| 淫秽高清视频在线观看| 欧美一级毛片孕妇| 变态另类成人亚洲欧美熟女| 中文字幕人妻丝袜一区二区| 尤物成人国产欧美一区二区三区| 久久人人精品亚洲av| 禁无遮挡网站| 婷婷精品国产亚洲av| 99热只有精品国产| 成人午夜高清在线视频| 3wmmmm亚洲av在线观看| 制服丝袜大香蕉在线| 热99re8久久精品国产| 中文资源天堂在线| 亚洲精品日韩av片在线观看 | 51国产日韩欧美| 免费看光身美女| 成人三级黄色视频| 久久亚洲真实| 国产综合懂色| 在线观看av片永久免费下载| 国产精品精品国产色婷婷| 久久久国产成人精品二区| 亚洲av二区三区四区| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 欧美乱色亚洲激情| 久久精品国产清高在天天线| 久99久视频精品免费| 久久精品亚洲精品国产色婷小说| 国产精品三级大全| 精品一区二区三区av网在线观看| 精品国产美女av久久久久小说| 久久久国产成人精品二区| 好男人电影高清在线观看| 午夜老司机福利剧场| 日韩国内少妇激情av| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 国产欧美日韩精品亚洲av| 亚洲第一电影网av| 国产高清激情床上av| 欧美中文日本在线观看视频| 男女午夜视频在线观看| 国产中年淑女户外野战色| 动漫黄色视频在线观看| 最近最新免费中文字幕在线| 国产精品99久久99久久久不卡| 国产黄色小视频在线观看| 在线国产一区二区在线| 国产精品爽爽va在线观看网站| 天堂网av新在线| 国产欧美日韩一区二区三| 欧美成人免费av一区二区三区| 亚洲真实伦在线观看| 女警被强在线播放| ponron亚洲| 久久九九热精品免费| 在线a可以看的网站| 狠狠狠狠99中文字幕| 好看av亚洲va欧美ⅴa在| 国产成人福利小说| 日韩大尺度精品在线看网址| 美女高潮喷水抽搐中文字幕| 国产三级在线视频| av女优亚洲男人天堂| 中文亚洲av片在线观看爽| 免费人成视频x8x8入口观看| 天天一区二区日本电影三级| 91在线精品国自产拍蜜月 | 午夜福利在线观看吧| 色老头精品视频在线观看| 在线观看美女被高潮喷水网站 | 小说图片视频综合网站| 性色avwww在线观看| 成人三级黄色视频| 国产精品一区二区三区四区免费观看 | 日韩国内少妇激情av| 久久久久久久久久黄片| aaaaa片日本免费| 小说图片视频综合网站| 综合色av麻豆| 蜜桃久久精品国产亚洲av| 人妻夜夜爽99麻豆av| 欧美性猛交╳xxx乱大交人| 一卡2卡三卡四卡精品乱码亚洲| 老司机在亚洲福利影院| 国产成人福利小说| 99久国产av精品| 成人特级黄色片久久久久久久| 亚洲性夜色夜夜综合| 日韩欧美国产在线观看| 人人妻,人人澡人人爽秒播| 一个人观看的视频www高清免费观看| 精品一区二区三区av网在线观看| 成人三级黄色视频| tocl精华| 国产免费男女视频| 深夜精品福利| www国产在线视频色| 国产精品久久久久久久电影 | 亚洲国产欧美网| 日本 av在线| 很黄的视频免费| 午夜两性在线视频| 欧美日韩乱码在线| 夜夜躁狠狠躁天天躁| 搡老岳熟女国产| 国产97色在线日韩免费| 搡老岳熟女国产| 成年女人毛片免费观看观看9| 久久国产乱子伦精品免费另类| 美女高潮的动态| 国产午夜精品久久久久久一区二区三区 | 成人欧美大片| 99在线人妻在线中文字幕| 人妻丰满熟妇av一区二区三区| 美女黄网站色视频| 一级毛片女人18水好多| 国产成人福利小说| 国产伦精品一区二区三区四那| 亚洲精品粉嫩美女一区| 久久久久久久久中文| 亚洲av电影在线进入| netflix在线观看网站| 亚洲av电影在线进入| 亚洲成av人片在线播放无| 听说在线观看完整版免费高清| 欧美高清成人免费视频www| 欧美日韩中文字幕国产精品一区二区三区| 小蜜桃在线观看免费完整版高清| 深爱激情五月婷婷| 日韩精品中文字幕看吧| 国产精品 国内视频| 麻豆成人av在线观看| 久久天躁狠狠躁夜夜2o2o| 精品不卡国产一区二区三区| 午夜精品在线福利| 91字幕亚洲| 久久久久性生活片| 亚洲欧美日韩无卡精品| 亚洲久久久久久中文字幕| 久久婷婷人人爽人人干人人爱| 丰满人妻熟妇乱又伦精品不卡| 欧美午夜高清在线| 精品久久久久久久久久免费视频| 三级国产精品欧美在线观看| 亚洲av美国av| 精品一区二区三区av网在线观看| 九色成人免费人妻av| a级一级毛片免费在线观看| 国产精品爽爽va在线观看网站| 中文资源天堂在线| 精品免费久久久久久久清纯| 婷婷亚洲欧美| a在线观看视频网站| 国产一级毛片七仙女欲春2| 国产真实乱freesex| 嫩草影视91久久| 国产亚洲欧美98| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆 | 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 中文字幕久久专区| 丰满乱子伦码专区| 免费看光身美女| 啪啪无遮挡十八禁网站| 动漫黄色视频在线观看| 午夜激情欧美在线| 精华霜和精华液先用哪个| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av美国av| 中文字幕高清在线视频| а√天堂www在线а√下载| 国产乱人伦免费视频| 高清毛片免费观看视频网站| 精品电影一区二区在线| 夜夜夜夜夜久久久久| 成人三级黄色视频| 级片在线观看| 久久香蕉国产精品| 久久精品国产亚洲av涩爱 | 中文字幕人妻熟人妻熟丝袜美 | 久久精品国产亚洲av涩爱 | xxx96com| 午夜福利免费观看在线| 在线十欧美十亚洲十日本专区| 天天躁日日操中文字幕| 老汉色∧v一级毛片| 国产中年淑女户外野战色| 国产真实伦视频高清在线观看 | 亚洲精品456在线播放app | 九九热线精品视视频播放| 欧美日韩国产亚洲二区| 国产精品 国内视频| 国产极品精品免费视频能看的| 99精品久久久久人妻精品| 女人高潮潮喷娇喘18禁视频| 法律面前人人平等表现在哪些方面| 免费电影在线观看免费观看| 99热精品在线国产| 观看美女的网站| 波野结衣二区三区在线 | av天堂中文字幕网| 99精品久久久久人妻精品| 午夜影院日韩av| 亚洲国产精品久久男人天堂| 99在线人妻在线中文字幕| 91九色精品人成在线观看| 色播亚洲综合网| 国产精品久久久久久久电影 | avwww免费| 久久人妻av系列| 欧美一区二区国产精品久久精品| 欧美精品啪啪一区二区三区| 国产av一区在线观看免费| 在线免费观看的www视频| 国产激情偷乱视频一区二区| 日本免费a在线| 久久国产精品影院| bbb黄色大片| 亚洲国产精品sss在线观看| 成人18禁在线播放| 少妇的逼好多水| 欧美+日韩+精品| 欧美成人一区二区免费高清观看| 精品人妻偷拍中文字幕| 人妻久久中文字幕网| 丰满人妻一区二区三区视频av | 天堂动漫精品| 91在线观看av| 国模一区二区三区四区视频| 91字幕亚洲| 老司机福利观看| 国产成年人精品一区二区| 亚洲国产日韩欧美精品在线观看 | 人妻夜夜爽99麻豆av| 成年免费大片在线观看| 高清日韩中文字幕在线| 亚洲午夜理论影院| 久久精品91无色码中文字幕| 91字幕亚洲| 亚洲av电影不卡..在线观看| 亚洲av成人av| 99国产精品一区二区三区| 一本久久中文字幕| 国产av在哪里看| www.熟女人妻精品国产| 久久久久久大精品| 欧美精品啪啪一区二区三区| 好男人在线观看高清免费视频| 男女午夜视频在线观看| 国产又黄又爽又无遮挡在线| 黄色视频,在线免费观看| 波野结衣二区三区在线 | 精品国产超薄肉色丝袜足j| 亚洲欧美激情综合另类| 精品国产美女av久久久久小说| 看片在线看免费视频| 欧美zozozo另类| 午夜福利在线观看免费完整高清在 | 欧美极品一区二区三区四区| 十八禁网站免费在线| 内射极品少妇av片p| 性色avwww在线观看| 神马国产精品三级电影在线观看| 亚洲第一电影网av| 精品久久久久久久毛片微露脸| 国产精品98久久久久久宅男小说| 18禁黄网站禁片午夜丰满| 少妇丰满av| 啦啦啦韩国在线观看视频| 最近最新中文字幕大全免费视频| 免费在线观看影片大全网站| 人妻久久中文字幕网| 国产av不卡久久| 欧美激情久久久久久爽电影| 精品熟女少妇八av免费久了| eeuss影院久久| 啦啦啦韩国在线观看视频| av视频在线观看入口| 人妻夜夜爽99麻豆av| 亚洲欧美精品综合久久99| 免费高清视频大片| 亚洲欧美日韩东京热| 久9热在线精品视频| 亚洲欧美日韩无卡精品| 高潮久久久久久久久久久不卡| 18禁美女被吸乳视频| 色哟哟哟哟哟哟| 亚洲第一欧美日韩一区二区三区| 国产国拍精品亚洲av在线观看 | 国产一区二区三区在线臀色熟女| 精品不卡国产一区二区三区| 桃色一区二区三区在线观看| 亚洲av免费在线观看| 老司机深夜福利视频在线观看| 久久午夜亚洲精品久久| 3wmmmm亚洲av在线观看| 99久久精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 老熟妇仑乱视频hdxx| 内地一区二区视频在线| 久久久久久国产a免费观看| 婷婷精品国产亚洲av| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线| 又黄又粗又硬又大视频| 久久久国产成人免费| 午夜a级毛片| 美女被艹到高潮喷水动态| 天天躁日日操中文字幕| 国产av在哪里看| 亚洲av成人精品一区久久| а√天堂www在线а√下载| 久久久久久人人人人人| 蜜桃久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 九九在线视频观看精品| 又爽又黄无遮挡网站| 99在线人妻在线中文字幕| 国语自产精品视频在线第100页| 一级黄片播放器| 欧美黄色片欧美黄色片| 欧美绝顶高潮抽搐喷水| 最近在线观看免费完整版| 中文字幕人成人乱码亚洲影| 国产av麻豆久久久久久久| 一区二区三区免费毛片| 日韩成人在线观看一区二区三区| 脱女人内裤的视频| 天天躁日日操中文字幕| 人人妻人人澡欧美一区二区| 午夜精品一区二区三区免费看| 久久久久久久精品吃奶| 神马国产精品三级电影在线观看| 日本一本二区三区精品| 亚洲aⅴ乱码一区二区在线播放| 国产三级中文精品| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区黑人| 午夜福利高清视频| 午夜久久久久精精品| 精品免费久久久久久久清纯| 亚洲专区中文字幕在线| 国产 一区 欧美 日韩| 岛国在线免费视频观看| 神马国产精品三级电影在线观看| 亚洲av成人精品一区久久| 亚洲黑人精品在线| 男女午夜视频在线观看| 欧美最黄视频在线播放免费| 久久精品91无色码中文字幕| 在线视频色国产色| 麻豆一二三区av精品| 久久久久九九精品影院| 午夜福利欧美成人| 18+在线观看网站| 久久九九热精品免费| 热99在线观看视频| 欧美日本视频| 国产欧美日韩一区二区精品| 丰满的人妻完整版| 天天一区二区日本电影三级| 中文字幕精品亚洲无线码一区| 欧美黑人欧美精品刺激|