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

    三峽船閘船舶吃水研究

    2017-11-07 09:21:14譚志榮張球林
    中國航海 2017年1期
    關(guān)鍵詞:過閘船閘樣條

    譚志榮, 張球林, 劉 釗

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院, 武漢 430063; 2.內(nèi)河航運(yùn)技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430063)

    三峽船閘船舶吃水研究

    譚志榮1,2, 張球林1, 劉 釗1

    (1.武漢理工大學(xué) 航運(yùn)學(xué)院, 武漢 430063; 2.內(nèi)河航運(yùn)技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430063)

    將三峽過閘船舶的吃水和載貨量視作一組離散的數(shù)據(jù),對船舶最大吃水、定額噸和實(shí)載噸等3個變量之間的關(guān)系進(jìn)行擬合分析,提出運(yùn)用插值方法確定實(shí)際過閘船舶的最優(yōu)吃水。對連續(xù)觀測的三峽船閘2014年的22 065艘次下行過閘船舶的數(shù)據(jù)進(jìn)行預(yù)處理,構(gòu)建船閘門檻水深與允許過閘船舶最大吃水的計(jì)算模型。運(yùn)用數(shù)理統(tǒng)計(jì)對8個區(qū)間的變量值進(jìn)行聚類分布。運(yùn)用MATLAB中分段線性、邊界條件為Lagrange和二階導(dǎo)數(shù)的三次樣條插值運(yùn)算方法對樣本數(shù)據(jù)進(jìn)行擬合。結(jié)果表明:以二階導(dǎo)數(shù)為邊界條件的三次樣條插值可解決該問題,當(dāng)允許過閘船舶的吃水控制在3.92~4.38 m時,過閘船舶裝載率為最佳。

    三峽船閘;吃水;載貨量;插值計(jì)算法

    Abstract: Fitting analysis is performed with the discrete data set of vessel maximum draft, rated tonnage and real load tonnage from actual statistics of the Three Gorges Lock. An interpolation method is brought up to determine the optimized draft. For the processing, the model which relates the vessel maximum draft to the water depth on ship-lock sill is built on the basis of 22 065 down-going vessel data acquired in 2014. The variables are clustered as eight sections. A piecewise linear three-order interpolation is carried out to find the relationship among variables with the Lagrange and second derivatives related boundary conditions. The fitted curve by cubic spline interpolation with the boundary conditions of second derivatives is shown valid. It concludes that vessels at drafts among 3.92 and 4.38 meters will have the optimal loading rate.

    Keywords: Three Gorges Lock; draft; loading rate; interpolation calculation method

    隨著長江干線航運(yùn)經(jīng)濟(jì)不斷發(fā)展和船型標(biāo)準(zhǔn)化不斷推進(jìn),三峽船閘過閘貨運(yùn)量保持增長態(tài)勢。[1]三峽船閘吃水控制標(biāo)準(zhǔn)關(guān)鍵技術(shù)研究專題成果[2]表明:過閘船舶吃水控制應(yīng)從3.3 m提升到3.9~4.3 m。在統(tǒng)計(jì)年報的基礎(chǔ)上,三峽通航管理部門自2012年開始加強(qiáng)三峽過閘船舶相關(guān)數(shù)據(jù)的統(tǒng)計(jì)工作,逐步完善對過閘船舶的過閘方式、貨種、到錨時間、開閘時間、最大吃水、定額噸及實(shí)載噸等數(shù)據(jù)的統(tǒng)計(jì)。通過對統(tǒng)計(jì)數(shù)據(jù)進(jìn)行挖掘發(fā)現(xiàn):盡管三峽船閘的貨運(yùn)量保持增長態(tài)勢,但從單艘船舶的載貨量上看,過閘船舶實(shí)際的載貨量沒有達(dá)到其額定值,部分船舶甚至是空載過閘,這必然是對三峽船閘過閘時空資源的一種無形中的損耗。該統(tǒng)計(jì)數(shù)據(jù)可為量化分析三峽船閘的船舶載貨量與吃水的關(guān)系提供計(jì)算基礎(chǔ);精細(xì)化控制過閘船舶的吃水有助于提升船閘的貨運(yùn)通過量,并為船型標(biāo)準(zhǔn)化提供技術(shù)支撐。

    目前,針對三峽船閘船舶吃水的研究主要集中在船閘規(guī)范設(shè)計(jì)和船舶實(shí)時吃水測量技術(shù)等方面。張義軍等[3]針對三峽船閘不同水位特征建立多級船閘檻上水深計(jì)算模型,該模型可為過閘船舶吃水控制標(biāo)準(zhǔn)的制訂提供理論根據(jù);鐘輝[4]利用網(wǎng)絡(luò)編程實(shí)現(xiàn)船舶實(shí)時吃水測量;孫坤凱[5]利用超聲波傳感器實(shí)現(xiàn)實(shí)時吃水檢測。已有研究多基于設(shè)計(jì)規(guī)范和試驗(yàn)?zāi)M的數(shù)值計(jì)算進(jìn)行,缺乏實(shí)際過閘船舶吃水與載貨量的實(shí)證分析。

    1 插值法

    在工程和數(shù)學(xué)應(yīng)用中,經(jīng)常有這樣一類數(shù)據(jù)處理問題,即在平面上給定一組離散點(diǎn)列,要求用一條曲線將這些點(diǎn)按次序連接起來,稱之為插值。[6]插值分析的適用范圍很廣:袁輝等[7]形式化描述實(shí)體參數(shù)的插值計(jì)算;高德偉等[8]運(yùn)用插值法進(jìn)行流量插值計(jì)算;姚傳宇等[9]運(yùn)用插值法和變步長滯環(huán)比較法改進(jìn)MPPT算法。

    鑒于高次插值既不收斂又不穩(wěn)定,低次插值既具有收斂性又具有穩(wěn)定性,認(rèn)為低次插值更具有實(shí)用價值。[10]然而,低次插值的光滑性較差,譬如分段線性插值多項(xiàng)式在插值區(qū)間內(nèi)僅具有連續(xù)性,在插值節(jié)點(diǎn)處有棱角,一階導(dǎo)數(shù)不存在。因此,這里在數(shù)據(jù)處理上運(yùn)用三次樣條插值進(jìn)行數(shù)據(jù)處理,既可解決分段線性插值光滑性較差的問題,又具有收斂性和穩(wěn)定性。

    2 數(shù)據(jù)預(yù)處理

    目前三峽水庫水位在145~175 m內(nèi)變化,三峽船閘門檻水深與船舶允許最大吃水密切相關(guān)。

    2.1門檻水深與船舶允許最大吃水的關(guān)系

    《內(nèi)河通航標(biāo)準(zhǔn)》規(guī)定水深、吃水及安全富余量的關(guān)系式為

    H=T+δ+ΔH

    (1)

    式(1)中:H為船閘門檻水深,m;T為允許過閘船舶的吃水深度,m;δ為船舶航行下沉量,m;ΔH為觸底安全富余量,m。

    由此,根據(jù)δ的確定方法,張義軍[11]通過實(shí)船試驗(yàn)得到相關(guān)的函數(shù)關(guān)系為

    δ=(-15.026K2+3.923K+0.026 7)H

    (2)

    (3)

    式(2)和(3)中:δ為船舶航行下沉量,m;H為船閘門檻水深,m;v為船舶對水的相對航速,m/s;F為閘室水下橫斷面面積,m2;f為船舶水下橫斷面面積,m2。

    結(jié)合H的變化情況,代表水深依次為5.00 m,5.20 m,5.25 m,5.50 m,5.8 m,6.00 m和7.00 m;根據(jù)式(1)~式(3),編程計(jì)算得到對應(yīng)的過閘船舶最大吃水T依次為3.92 m,4.11 m,4.15 m,4.38 m,4.65 m,4.84 m和5.76 m。在2014年的統(tǒng)計(jì)樣本中,將船舶滿載吃水的數(shù)值作為最大吃水的樣本值,即最大吃水的區(qū)間為[3.92,5.76]m。

    2.2過閘船舶統(tǒng)計(jì)數(shù)據(jù)預(yù)處理

    在2012—2014年的數(shù)據(jù)中,選取2014年三峽船閘連續(xù)的過閘下行船舶數(shù)據(jù)作為樣本(見表1)。

    表1 2014年三峽船閘下行船舶情況

    經(jīng)過統(tǒng)計(jì)發(fā)現(xiàn):大量樣本船舶的吃水相同,但對應(yīng)的載貨量不同,因此若對每組數(shù)據(jù)都在MATLAB中進(jìn)行相關(guān)的插值計(jì)算,較為混亂和復(fù)雜,且會導(dǎo)致結(jié)果失真。這里對統(tǒng)計(jì)數(shù)據(jù)進(jìn)行預(yù)處理,即:根據(jù)求取的7個吃水特征值,分為8個區(qū)間提取船舶數(shù)據(jù);對每個區(qū)間中所有船舶最大吃水下的船舶定額噸和實(shí)載噸求平均值;對每個區(qū)間中的船舶最大吃水求平均值。具體數(shù)據(jù)預(yù)處理結(jié)果見表2。

    表2 統(tǒng)計(jì)數(shù)據(jù)預(yù)處理結(jié)果

    3 擬合分析

    3.1最大吃水與定額噸

    對處理后的數(shù)據(jù)進(jìn)行擬合計(jì)算,橫坐標(biāo)為最大吃水,縱坐標(biāo)為定額噸。運(yùn)用分段線性(piecewise linear)三次樣條插值(spline)法對其進(jìn)行擬合計(jì)算。

    圖1為最大吃水與定額噸的關(guān)系,其中:第1條曲線(左側(cè))表示線性擬合計(jì)算;第2條曲線(中間)表示Lagrange邊界條件的樣條插值;第3條曲線(右側(cè))為二階導(dǎo)數(shù)邊界條件的樣條插值。

    圖1 最大吃水與定額噸的關(guān)系

    由圖1可知:利用線性插值計(jì)算得到的圖像光滑性較差;而利用三次樣條插值擬合計(jì)算得到的圖像光滑性很好,其中運(yùn)用二階導(dǎo)數(shù)邊界條件(第3條曲線)擬合結(jié)果的精度比Lagrange邊界條件(第2條曲線)擬合結(jié)果的精度高。此外,還可看到在第一段區(qū)間內(nèi)出現(xiàn)一個峰值。文獻(xiàn)[8]中也指出,當(dāng)子區(qū)間較小時(插值節(jié)點(diǎn)較多),樣條插值的結(jié)果精度較高。這說明設(shè)置的第一個區(qū)間跨度較大,因此后面需進(jìn)行單獨(dú)分析。

    3.2最大吃水與實(shí)載噸

    設(shè)定橫坐標(biāo)為最大吃水、縱坐標(biāo)為船舶實(shí)載噸進(jìn)行插值計(jì)算,得到最大吃水與船舶實(shí)載噸的關(guān)系見圖2。

    圖2 最大吃水與船舶實(shí)載噸的關(guān)系

    圖2中運(yùn)用三次樣條插值擬合得到的曲線出現(xiàn)2個峰值,其中:第1個峰值與圖1中的情況相同;對于第2個峰值,隨著船舶最大吃水量增加,船舶實(shí)載貨物量呈先增后減的趨勢。對第2個峰值附近的觀測值進(jìn)行驗(yàn)證,結(jié)果見圖3。

    圖3中第2個峰值附近出現(xiàn)的觀測值與圖2中擬合的基本一致,這表明樣條曲線擬合是合理的,當(dāng)船舶最大吃水達(dá)到5 m后,曲線表示的實(shí)載噸均呈下降趨勢。這表明:通過三峽船閘的大型船舶裝載率偏低,且大型船舶必須通過減載才能滿足三峽船閘過閘的吃水控制要求。

    3.3擬合比較

    為進(jìn)一步尋找最大吃水、船舶定額噸及船舶實(shí)載噸三者之間的關(guān)系,圖4中設(shè)定橫坐標(biāo)為最大吃水、縱坐標(biāo)為載貨量(其中,上部分曲線代表定額噸,下部分曲線代表實(shí)載噸)進(jìn)行曲線對比。

    圖4 定額噸與實(shí)載噸比較

    由圖4可知:

    1)在分段進(jìn)行插值計(jì)算時,三次樣條插值計(jì)算開始出現(xiàn)第1個峰值,也就是最大吃水T<4.11 m時(圖1和圖2中出現(xiàn)峰值的部分),數(shù)據(jù)關(guān)系相差較大,有必要對該分段作進(jìn)一步的分析。

    2)當(dāng)吃水在4.15~4.38 m時,船閘實(shí)際通過的貨物量與船舶滿載時的定額噸相差最小,這表明吃水在該區(qū)間內(nèi),船舶裝載率和船閘過貨效率最高。

    3)根據(jù)現(xiàn)有的三峽船閘吃水控制標(biāo)準(zhǔn),當(dāng)船舶吃水達(dá)到一定值時,隨著船舶最大吃水增加,過閘船舶的實(shí)際裝載率反而下降。

    3.4討論

    當(dāng)船舶最大吃水T<4.11 m時,擬合曲線的差異較大。根據(jù)2014年的實(shí)際數(shù)據(jù)和區(qū)間劃定,對最大吃水分為(3.92,4.11] m和(0.40,3.92] m區(qū)段進(jìn)行再擬合分析。

    1)當(dāng)船舶最大吃水在(3.92,4.11]m區(qū)段時,得到定額噸與實(shí)載噸關(guān)系見圖5。

    由圖5可知:當(dāng)最大吃水在該區(qū)段時,船舶定額噸曲線與實(shí)載噸曲線基本重疊,說明數(shù)值相差很小,船舶裝載率較高。

    2)當(dāng)船舶最大吃水在[0.40,3.92]m時,對此類數(shù)據(jù)進(jìn)行篩選,發(fā)現(xiàn)很多船舶出現(xiàn)空載現(xiàn)象,該區(qū)段對插值擬合關(guān)系的計(jì)算結(jié)論不準(zhǔn)確。因此,對于該區(qū)段的樣本數(shù)據(jù),采用散點(diǎn)分布圖(見圖6)描述定額噸與實(shí)載噸的關(guān)系。

    圖5 定額噸與實(shí)載噸關(guān)系

    圖6 最大吃水在(0.40,3.92]m區(qū)段時實(shí)載噸與定額噸關(guān)系

    由圖6可知:船舶最大吃水為[0.40,2.00]m時,船舶裝載率非常低。這進(jìn)一步佐證了在載貨效率方面,船舶小型化已不適應(yīng)過閘需求。

    由于在線性分段時最大吃水4.11 m和4.15 m為相鄰區(qū)間插值點(diǎn),圖4和圖5的擬合特征基本一致,因此可得到當(dāng)船舶最大吃水在[3.92,4.38]m時船舶實(shí)際裝載率最大。

    4 結(jié)束語

    在三峽船閘實(shí)際運(yùn)行中,過閘船舶的實(shí)際裝載率不高已成為制約過閘貨運(yùn)量提升的瓶頸。運(yùn)用插值方法對連續(xù)觀測的過閘船舶數(shù)據(jù)進(jìn)行分析,得出船舶最大吃水、定額噸和實(shí)載噸的擬合關(guān)系。結(jié)果表明:當(dāng)過閘船舶的吃水在(3.92,4.38]m時,過閘船舶實(shí)際載貨效率為最優(yōu);當(dāng)船舶吃水在

    (0.40,2.00]m時,船閘過貨效率非常低。2014年樣本的統(tǒng)計(jì)結(jié)論的特征與近3 a年報統(tǒng)計(jì)結(jié)論的特征一致,佐證了主管部門2015年提出的過閘船舶吃水控制在4.30 m的經(jīng)驗(yàn)管理具有科學(xué)性。下一階段可在MATLAB模擬仿真中,通過設(shè)置更多的數(shù)據(jù)梯度來得到更加精準(zhǔn)的擬合圖形。此外,選擇的是統(tǒng)計(jì)數(shù)據(jù)的樣本實(shí)際過閘船舶最大吃水,下一階段需對上下行過閘船舶的實(shí)際吃水進(jìn)行分類統(tǒng)計(jì)。該方法僅從過閘船舶實(shí)際裝載率優(yōu)化的角度研究船舶吃水問題,下一階段可結(jié)合船閘服務(wù)質(zhì)量等內(nèi)容進(jìn)行綜合研究。

    [1] 李灼.擴(kuò)大三峽船閘通航能力研究[D].重慶:重慶交通大學(xué),2013.

    [2] 張義軍.過閘船舶吃水控制標(biāo)準(zhǔn)對三峽船閘通過量的影響[J].水運(yùn)管理,2012,34(3):12-14.

    [3] 張義軍,杜嘉立.多級船閘檻上水深計(jì)算模型[J].大連海事大學(xué)學(xué)報,2012,38(4):23-25.

    [4] 鐘輝.三峽通航船舶吃水遠(yuǎn)程監(jiān)測與預(yù)警系統(tǒng)的研究[D].大連:大連海事大學(xué),2011.

    [5] 孫坤凱.三峽船舶吃水量檢測門同步升降系統(tǒng)的研究[D].大連:大連海事大學(xué),2010.

    [6] 司守奎,孫兆亮,孫璽菁.數(shù)學(xué)建模算法與應(yīng)用[M].北京:國防工業(yè)出版社,2011.

    [7] 袁輝,王鳳山,任睿.動荷段實(shí)體參數(shù)插值計(jì)算模型組件[J].解放軍理工大學(xué)學(xué)報(自然科學(xué)版),2016,17(4):386-390.

    [8] 高德偉,黃俊,魏邦記.利用樣條插值法進(jìn)行流量插值計(jì)算[J].東北水利水電,2000,18(7):49-50.

    [9] 姚傳宇,房俊龍,王泉宇,等.基于插值計(jì)算法和變步長滯環(huán)比較法結(jié)合的MPPT算法研究[J].電測與儀表,2015,52(23):52-57.

    [10] 楊永川,楊軻.空間插值法在熱環(huán)境流動觀測中的應(yīng)用[J].中南大學(xué)學(xué)報(自然科學(xué)版),2012,43(9):3741-3748.

    [11] 張義軍.三峽船閘允許過閘船舶最大吃水的確定[J].水運(yùn)管理,2014,36(2):22-24.

    StudyonVesselDraftinThreeGorgesLockTANZhirong

    1,2,ZHANGQiulin1,LIUZhao1

    (1. School of Navigation, Wuhan University of Technology, Wuhan 430063, China; 2.Hubei Key Laboratory of Inland Shipping Technology, Wuhan 430063, China)

    U641.7

    A

    2016-11-21

    國家自然科學(xué)基金(51479156);湖北省科技支撐項(xiàng)目(對外合作類)(2015BHE004)

    譚志榮(1978—), 男, 湖北武漢人,副教授,碩士生導(dǎo)師,博士,從事載運(yùn)工具運(yùn)用工程與交通系統(tǒng)安全性評價研究。 E-mail:tanzhi@whut.edu.cn

    1000-4653(2017)01-0069-04

    猜你喜歡
    過閘船閘樣條
    一元五次B樣條擬插值研究
    江蘇內(nèi)河航道船舶優(yōu)先過閘現(xiàn)狀與經(jīng)驗(yàn)
    三峽船舶過閘計(jì)劃的智能編制與優(yōu)化
    抗疫,在三峽兩壩船閘水域
    水上消防(2020年4期)2021-01-04 01:51:30
    船閘
    五級運(yùn)行閘次排序?qū)θ龒{船閘上行過閘效率的影響
    三次參數(shù)樣條在機(jī)床高速高精加工中的應(yīng)用
    三次樣條和二次刪除相輔助的WASD神經(jīng)網(wǎng)絡(luò)與日本人口預(yù)測
    軟件(2017年6期)2017-09-23 20:56:27
    基于樣條函數(shù)的高精度電子秤設(shè)計(jì)
    用于船閘撞擊防護(hù)微孔塑料的計(jì)算與分析
    久久精品夜夜夜夜夜久久蜜豆| 禁无遮挡网站| 乱系列少妇在线播放| 蜜桃久久精品国产亚洲av| 亚洲成av人片在线播放无| 日本-黄色视频高清免费观看| 欧美日韩精品成人综合77777| 国产精品爽爽va在线观看网站| 男人和女人高潮做爰伦理| 亚洲不卡免费看| 两性午夜刺激爽爽歪歪视频在线观看| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 好男人视频免费观看在线| 如何舔出高潮| 网址你懂的国产日韩在线| 国产 一区精品| 男女国产视频网站| 韩国高清视频一区二区三区| 在线播放无遮挡| 国内少妇人妻偷人精品xxx网站| 最近最新中文字幕免费大全7| 色播亚洲综合网| 又爽又黄无遮挡网站| 欧美日本视频| 国产69精品久久久久777片| 精品久久久久久久人妻蜜臀av| 亚洲第一区二区三区不卡| 色综合亚洲欧美另类图片| 一区二区三区四区激情视频| 午夜精品国产一区二区电影 | 免费黄色在线免费观看| 黄色欧美视频在线观看| av福利片在线观看| 亚州av有码| 亚洲va在线va天堂va国产| 九九爱精品视频在线观看| 成人高潮视频无遮挡免费网站| 亚洲真实伦在线观看| 亚洲18禁久久av| 日本一二三区视频观看| 一级毛片aaaaaa免费看小| 亚洲欧美精品自产自拍| 欧美精品国产亚洲| 国产淫片久久久久久久久| 亚洲av成人精品一区久久| 欧美不卡视频在线免费观看| 久久亚洲精品不卡| 国产在视频线精品| 啦啦啦观看免费观看视频高清| 国产午夜福利久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精华霜和精华液先用哪个| 搡老妇女老女人老熟妇| 国产在线一区二区三区精 | 超碰av人人做人人爽久久| 国产亚洲一区二区精品| 大香蕉久久网| 草草在线视频免费看| 91精品伊人久久大香线蕉| 国产91av在线免费观看| 少妇高潮的动态图| 男女下面进入的视频免费午夜| 岛国在线免费视频观看| 成年av动漫网址| 精品酒店卫生间| 3wmmmm亚洲av在线观看| 一边摸一边抽搐一进一小说| 欧美97在线视频| 成人毛片60女人毛片免费| 国产亚洲91精品色在线| 亚洲欧美一区二区三区国产| 高清视频免费观看一区二区 | 国产黄色小视频在线观看| 美女高潮的动态| 啦啦啦观看免费观看视频高清| 国产精品人妻久久久久久| 亚洲国产精品成人久久小说| 国产片特级美女逼逼视频| 久久国内精品自在自线图片| 日本黄色片子视频| 午夜精品一区二区三区免费看| 秋霞伦理黄片| 精品人妻熟女av久视频| 欧美极品一区二区三区四区| a级一级毛片免费在线观看| 1000部很黄的大片| 亚洲av.av天堂| 日日摸夜夜添夜夜爱| 久久精品夜色国产| 黄色欧美视频在线观看| 看十八女毛片水多多多| 91aial.com中文字幕在线观看| 在线免费观看的www视频| 免费观看性生交大片5| 精品一区二区三区视频在线| 高清在线视频一区二区三区 | 亚洲欧美日韩东京热| 你懂的网址亚洲精品在线观看 | 国产精品一区www在线观看| 秋霞在线观看毛片| 国产精品久久电影中文字幕| 丝袜喷水一区| 国产伦精品一区二区三区四那| 伦精品一区二区三区| 黄色一级大片看看| 一级黄色大片毛片| 午夜福利在线观看免费完整高清在| 永久免费av网站大全| 久久99热这里只有精品18| 国产在视频线精品| 日韩成人伦理影院| 美女国产视频在线观看| 在线免费十八禁| 波多野结衣巨乳人妻| 中文欧美无线码| 两个人视频免费观看高清| 欧美成人午夜免费资源| 国产极品精品免费视频能看的| 变态另类丝袜制服| 久久草成人影院| 国产精品一区二区三区四区免费观看| 国产精品一区二区三区四区久久| 日本黄色片子视频| 日本黄色视频三级网站网址| 一级二级三级毛片免费看| 欧美成人精品欧美一级黄| av福利片在线观看| 婷婷色综合大香蕉| 美女xxoo啪啪120秒动态图| 亚洲va在线va天堂va国产| 中文字幕亚洲精品专区| 中文字幕av在线有码专区| 你懂的网址亚洲精品在线观看 | 日韩欧美精品免费久久| 久久久国产成人免费| 熟女人妻精品中文字幕| 女人被狂操c到高潮| 国产又黄又爽又无遮挡在线| 成人午夜精彩视频在线观看| 免费看光身美女| 永久网站在线| 午夜福利网站1000一区二区三区| 久久精品国产自在天天线| 日韩大片免费观看网站 | 日韩,欧美,国产一区二区三区 | 麻豆一二三区av精品| 国产午夜精品久久久久久一区二区三区| 国产黄色小视频在线观看| 久久久久久久久久久丰满| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 晚上一个人看的免费电影| av福利片在线观看| 久久综合国产亚洲精品| 久久久a久久爽久久v久久| 91午夜精品亚洲一区二区三区| 69人妻影院| 午夜久久久久精精品| 91久久精品电影网| 亚洲五月天丁香| 青春草视频在线免费观看| 自拍偷自拍亚洲精品老妇| 免费搜索国产男女视频| 日韩视频在线欧美| 国产不卡一卡二| 免费观看的影片在线观看| 小蜜桃在线观看免费完整版高清| 中文欧美无线码| 亚洲国产最新在线播放| 色5月婷婷丁香| 精品一区二区免费观看| 久久精品国产亚洲网站| 免费av不卡在线播放| 久久精品影院6| 美女高潮的动态| 精品国产一区二区三区久久久樱花 | 我要搜黄色片| 日韩视频在线欧美| 国产 一区精品| 美女cb高潮喷水在线观看| 久久久欧美国产精品| 久久精品国产自在天天线| 精品欧美国产一区二区三| 在线播放无遮挡| 国产精品久久久久久精品电影| 成人性生交大片免费视频hd| 欧美最新免费一区二区三区| 听说在线观看完整版免费高清| 国产精品永久免费网站| 日本与韩国留学比较| 毛片一级片免费看久久久久| 色5月婷婷丁香| 国产精品,欧美在线| 三级毛片av免费| 国产白丝娇喘喷水9色精品| 九色成人免费人妻av| 如何舔出高潮| 国产精品国产三级专区第一集| 卡戴珊不雅视频在线播放| 国产成人精品一,二区| 日本与韩国留学比较| 男人的好看免费观看在线视频| 国产伦精品一区二区三区四那| 在线观看一区二区三区| 国产亚洲一区二区精品| 两个人的视频大全免费| 97热精品久久久久久| 精品久久久久久久人妻蜜臀av| 亚洲欧美成人综合另类久久久 | 久久久久精品久久久久真实原创| 免费在线观看成人毛片| 亚洲中文字幕一区二区三区有码在线看| 99九九线精品视频在线观看视频| 久久精品久久久久久久性| 国产精品久久电影中文字幕| 一个人观看的视频www高清免费观看| 麻豆精品久久久久久蜜桃| 男女那种视频在线观看| 伊人久久精品亚洲午夜| 日本黄大片高清| 黄片无遮挡物在线观看| 九九久久精品国产亚洲av麻豆| 午夜亚洲福利在线播放| 最近的中文字幕免费完整| 成人综合一区亚洲| 亚洲中文字幕日韩| 国产一区二区在线观看日韩| 欧美一区二区国产精品久久精品| 亚洲av免费高清在线观看| 国产在线一区二区三区精 | 精品人妻视频免费看| 免费观看的影片在线观看| 国产免费男女视频| 亚洲国产欧美人成| 男女那种视频在线观看| 久久久久精品久久久久真实原创| 国产美女午夜福利| av天堂中文字幕网| 亚洲av熟女| 国产欧美日韩精品一区二区| 精品久久久久久电影网 | av国产免费在线观看| 国产真实乱freesex| 女人久久www免费人成看片 | av福利片在线观看| 欧美zozozo另类| 久久久欧美国产精品| 日韩大片免费观看网站 | 蜜臀久久99精品久久宅男| 午夜福利视频1000在线观看| h日本视频在线播放| 欧美性猛交╳xxx乱大交人| 国产一区有黄有色的免费视频 | 高清av免费在线| 欧美另类亚洲清纯唯美| 黄片无遮挡物在线观看| 人妻系列 视频| 日本免费在线观看一区| 亚洲五月天丁香| 国产午夜精品论理片| 高清av免费在线| 99久久成人亚洲精品观看| 中文欧美无线码| 99久国产av精品国产电影| 成人欧美大片| 精品少妇黑人巨大在线播放 | 纵有疾风起免费观看全集完整版 | 午夜福利网站1000一区二区三区| 国产亚洲91精品色在线| 美女cb高潮喷水在线观看| 国产精品麻豆人妻色哟哟久久 | 亚洲精品日韩在线中文字幕| a级毛片免费高清观看在线播放| 好男人在线观看高清免费视频| 直男gayav资源| 看十八女毛片水多多多| 国产乱人视频| 欧美激情在线99| 久久久欧美国产精品| 亚洲自偷自拍三级| 老师上课跳d突然被开到最大视频| 精品人妻视频免费看| 国产成人aa在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产毛片a区久久久久| 嫩草影院精品99| 大又大粗又爽又黄少妇毛片口| 精品一区二区免费观看| 别揉我奶头 嗯啊视频| 蜜臀久久99精品久久宅男| 99在线视频只有这里精品首页| 免费av不卡在线播放| 国产黄片视频在线免费观看| 熟妇人妻久久中文字幕3abv| 狂野欧美激情性xxxx在线观看| 亚洲三级黄色毛片| 亚洲国产欧美在线一区| 久久久精品大字幕| 亚洲国产欧美在线一区| 天堂网av新在线| 七月丁香在线播放| 蜜桃久久精品国产亚洲av| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲第一区二区三区不卡| 国产精品久久视频播放| 亚洲成色77777| 成人午夜高清在线视频| 国产成人freesex在线| 嫩草影院精品99| 免费一级毛片在线播放高清视频| 国产国拍精品亚洲av在线观看| 国产大屁股一区二区在线视频| 黄片wwwwww| 免费看av在线观看网站| 亚洲自拍偷在线| 少妇熟女欧美另类| av在线观看视频网站免费| 午夜爱爱视频在线播放| 国产精品综合久久久久久久免费| av女优亚洲男人天堂| 国产午夜福利久久久久久| 日韩高清综合在线| 国产老妇女一区| 嫩草影院入口| 老司机影院毛片| 国产免费视频播放在线视频 | 日本色播在线视频| 免费看日本二区| 久久99精品国语久久久| 人妻系列 视频| 男人舔奶头视频| 男人和女人高潮做爰伦理| 你懂的网址亚洲精品在线观看 | 国产乱人视频| 非洲黑人性xxxx精品又粗又长| 国内少妇人妻偷人精品xxx网站| 精品欧美国产一区二区三| 日本黄色视频三级网站网址| 性色avwww在线观看| 国产精品1区2区在线观看.| 国产高清有码在线观看视频| www.av在线官网国产| 一个人免费在线观看电影| 国产精品野战在线观看| 久久久精品94久久精品| av专区在线播放| 日本猛色少妇xxxxx猛交久久| 久久久精品大字幕| 床上黄色一级片| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看| 一级毛片aaaaaa免费看小| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 亚洲av福利一区| 亚洲av成人精品一区久久| 在线免费观看的www视频| 最近最新中文字幕大全电影3| 男人舔女人下体高潮全视频| 日韩亚洲欧美综合| 天堂av国产一区二区熟女人妻| 国产成人a区在线观看| 亚洲第一区二区三区不卡| av在线天堂中文字幕| 午夜久久久久精精品| 亚洲欧美成人综合另类久久久 | 精品久久久久久久久亚洲| 91精品国产九色| 日本午夜av视频| 女的被弄到高潮叫床怎么办| 色噜噜av男人的天堂激情| 久久久久久久国产电影| 国产一区二区在线观看日韩| 蜜桃久久精品国产亚洲av| 国产精品av视频在线免费观看| 久久久久久久国产电影| 午夜福利网站1000一区二区三区| 天天躁日日操中文字幕| 伦精品一区二区三区| 成人特级av手机在线观看| 午夜视频国产福利| 免费大片18禁| 黄色配什么色好看| 精品久久久久久久末码| 久久精品91蜜桃| 老司机影院成人| 成人毛片60女人毛片免费| 亚洲欧美中文字幕日韩二区| 神马国产精品三级电影在线观看| 一级av片app| av国产免费在线观看| 99热这里只有精品一区| 亚洲av男天堂| 精品久久久噜噜| 精品久久国产蜜桃| 亚洲无线观看免费| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花 | 中文字幕免费在线视频6| 色噜噜av男人的天堂激情| 九九爱精品视频在线观看| 久久久成人免费电影| 丝袜喷水一区| 亚洲国产欧美人成| 国产精品99久久久久久久久| 嫩草影院入口| 日本黄大片高清| 国产美女午夜福利| 精品人妻视频免费看| 中文字幕免费在线视频6| 久久99热这里只有精品18| 久久久a久久爽久久v久久| 亚洲av福利一区| 看片在线看免费视频| 久久久久久大精品| 亚洲欧美日韩东京热| 国产精品久久久久久av不卡| av天堂中文字幕网| 99在线视频只有这里精品首页| 中文字幕av成人在线电影| 久久99热6这里只有精品| 亚洲成人中文字幕在线播放| 日韩欧美国产在线观看| 国产午夜精品久久久久久一区二区三区| 午夜福利网站1000一区二区三区| 国产精品女同一区二区软件| 真实男女啪啪啪动态图| 村上凉子中文字幕在线| 久久精品夜色国产| 99视频精品全部免费 在线| 一夜夜www| 国产免费一级a男人的天堂| 亚洲欧美日韩高清专用| av在线观看视频网站免费| 99久国产av精品| 亚洲国产精品专区欧美| 国产免费又黄又爽又色| 99九九线精品视频在线观看视频| 国产精品.久久久| 最后的刺客免费高清国语| 精品一区二区三区人妻视频| 99热这里只有是精品在线观看| www日本黄色视频网| 精品人妻一区二区三区麻豆| 精品久久国产蜜桃| 日韩国内少妇激情av| 免费搜索国产男女视频| 精品一区二区三区人妻视频| 99热这里只有是精品在线观看| 精品久久久久久久久亚洲| 尤物成人国产欧美一区二区三区| 国产av不卡久久| 日韩在线高清观看一区二区三区| 亚洲av熟女| 亚洲中文字幕日韩| 精品无人区乱码1区二区| 日本色播在线视频| 久久久a久久爽久久v久久| 国产精品日韩av在线免费观看| 熟妇人妻久久中文字幕3abv| 日韩高清综合在线| 亚洲精品亚洲一区二区| 久热久热在线精品观看| 最近最新中文字幕免费大全7| 成人欧美大片| 我要搜黄色片| 人体艺术视频欧美日本| 久久精品夜色国产| 久久久久久九九精品二区国产| 国产片特级美女逼逼视频| 老师上课跳d突然被开到最大视频| 精品熟女少妇av免费看| 久久午夜福利片| 国语自产精品视频在线第100页| 日本爱情动作片www.在线观看| 国产伦理片在线播放av一区| 中文字幕熟女人妻在线| 99热精品在线国产| 国产成年人精品一区二区| 国产色爽女视频免费观看| 久久这里只有精品中国| 99久国产av精品| 少妇人妻一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 寂寞人妻少妇视频99o| 大香蕉97超碰在线| 亚洲第一区二区三区不卡| 毛片女人毛片| 午夜福利网站1000一区二区三区| 久久久国产成人免费| 久久久久久九九精品二区国产| 一本一本综合久久| 亚洲中文字幕一区二区三区有码在线看| 午夜精品在线福利| 欧美区成人在线视频| 听说在线观看完整版免费高清| 国产在视频线在精品| 精品久久久久久电影网 | 夫妻性生交免费视频一级片| 女人十人毛片免费观看3o分钟| 99久久精品一区二区三区| 国产亚洲午夜精品一区二区久久 | 午夜激情福利司机影院| 国产午夜精品久久久久久一区二区三区| 成年免费大片在线观看| 欧美又色又爽又黄视频| 非洲黑人性xxxx精品又粗又长| 汤姆久久久久久久影院中文字幕 | 69av精品久久久久久| 亚洲欧美日韩东京热| av视频在线观看入口| h日本视频在线播放| 一级毛片我不卡| 人妻夜夜爽99麻豆av| av女优亚洲男人天堂| 午夜福利高清视频| 麻豆av噜噜一区二区三区| 亚洲成人中文字幕在线播放| 精品久久久久久成人av| 日本午夜av视频| 国产69精品久久久久777片| 亚洲自拍偷在线| 美女高潮的动态| 欧美又色又爽又黄视频| 国产精品一区www在线观看| 免费观看在线日韩| 国产在视频线精品| 国产老妇女一区| 午夜老司机福利剧场| 国产精品一区二区在线观看99 | 欧美日韩在线观看h| 色视频www国产| 美女大奶头视频| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| 午夜爱爱视频在线播放| 人妻夜夜爽99麻豆av| 男人的好看免费观看在线视频| 国产精品1区2区在线观看.| 丝袜美腿在线中文| 不卡视频在线观看欧美| av在线天堂中文字幕| 久久久欧美国产精品| 国产一区二区亚洲精品在线观看| 男女那种视频在线观看| 日日撸夜夜添| 国产av一区在线观看免费| 久久亚洲精品不卡| 99热这里只有精品一区| 久久亚洲精品不卡| 蜜桃久久精品国产亚洲av| 3wmmmm亚洲av在线观看| 18禁动态无遮挡网站| 日本色播在线视频| 又爽又黄无遮挡网站| 日本与韩国留学比较| 国产69精品久久久久777片| 只有这里有精品99| 最近的中文字幕免费完整| 日韩人妻高清精品专区| 亚洲欧美中文字幕日韩二区| 亚洲精品aⅴ在线观看| 1024手机看黄色片| 日韩欧美 国产精品| 国产精品久久久久久精品电影小说 | 尤物成人国产欧美一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 老司机福利观看| 日本午夜av视频| 中文字幕制服av| 日本wwww免费看| 国产在线一区二区三区精 | 久久精品久久精品一区二区三区| 国产成人一区二区在线| 亚洲av电影在线观看一区二区三区 | 麻豆久久精品国产亚洲av| 亚洲四区av| 三级国产精品欧美在线观看| 日韩亚洲欧美综合| 国产免费一级a男人的天堂| 亚洲国产精品成人久久小说| 国内少妇人妻偷人精品xxx网站| 亚洲欧美中文字幕日韩二区| 久久热精品热| 韩国av在线不卡| 亚洲欧美精品专区久久| 超碰97精品在线观看| 精品久久久久久久久亚洲| 午夜精品国产一区二区电影 | 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 精品人妻熟女av久视频| 婷婷色麻豆天堂久久 | 国产精品综合久久久久久久免费| 国产人妻一区二区三区在| 日韩av在线免费看完整版不卡| 非洲黑人性xxxx精品又粗又长| 亚洲精品久久久久久婷婷小说 | a级一级毛片免费在线观看| 亚洲最大成人av| 欧美激情国产日韩精品一区| 99久久无色码亚洲精品果冻| 成人亚洲欧美一区二区av| 国产亚洲精品久久久com| 日韩欧美 国产精品| 免费黄色在线免费观看| 国产三级在线视频| 国产精品蜜桃在线观看| 亚洲美女视频黄频| 99久久成人亚洲精品观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 白带黄色成豆腐渣| 国模一区二区三区四区视频|