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

    上海沿海有效波高和最大波高關(guān)系分析

    2017-01-09 06:44:29朱智慧沈其艷常亞楠
    海洋預(yù)報(bào) 2016年6期
    關(guān)鍵詞:波高浮標(biāo)海浪

    朱智慧,沈其艷,常亞楠

    上海沿海有效波高和最大波高關(guān)系分析

    朱智慧,沈其艷,常亞楠

    (上海海洋氣象臺(tái),上海201306)

    運(yùn)用線性回歸方法,對(duì)上海沿海5個(gè)浮標(biāo)站逐時(shí)的有效波高和最大波高觀測(cè)資料進(jìn)行分析,建立了這兩種波高的回歸方程,獲得了估計(jì)最大波高的客觀方法。對(duì)5個(gè)浮標(biāo)站,一元一次回歸方程都能很好地反映最大波高和有效波高的關(guān)系。最大波高通用回歸方程對(duì)長江口船標(biāo)資料的試報(bào)結(jié)果較好,平均絕對(duì)誤差僅為0.1 m。WW3海浪模式對(duì)5個(gè)浮標(biāo)站的有效波高預(yù)報(bào)與實(shí)況比較接近,利用該模式的有效波高進(jìn)行最大波高的預(yù)報(bào),預(yù)報(bào)最大波高與實(shí)測(cè)最大波高具有較好的相關(guān)性,5個(gè)浮標(biāo)站平均誤差最大為-0.17 m,最小為-0.05 m,平均絕對(duì)誤差都為0.4 m,結(jié)果比較理想。

    有效波高;最大波高;回歸分析

    1 引言

    隨著海洋經(jīng)濟(jì)的發(fā)展,及時(shí)提供海上的海洋環(huán)境條件,成為海洋氣象預(yù)報(bào)的迫切任務(wù)。波高對(duì)于海上運(yùn)輸、海水養(yǎng)殖、海事救撈和海洋工程等生產(chǎn)活動(dòng),是必須考慮的重要因素。此外,在信息化時(shí)代,氣象水文條件是海戰(zhàn)環(huán)境的重要組成部分,對(duì)艦艇編隊(duì)活動(dòng)、軍事演習(xí)、登陸作戰(zhàn)、海洋維權(quán)等軍事行動(dòng)有著非常實(shí)用的價(jià)值。鄭崇偉等[1-3]以釣魚島、黃巖島和南海島礁等重要島礁海域?yàn)檠芯繉?duì)象,利用模擬的海浪場(chǎng)數(shù)據(jù),對(duì)研究海域的波候特征進(jìn)行了系統(tǒng)研究,為海洋資源開發(fā)以及戰(zhàn)場(chǎng)環(huán)境建設(shè)提供了科學(xué)指導(dǎo)。如今,數(shù)值預(yù)報(bào)正成為海洋氣象預(yù)報(bào)的最主要手段[4-7],數(shù)值模式輸出的波高為有效波高,但是在實(shí)際中,最大波高是一個(gè)重要的參數(shù),而在目前的海洋氣象預(yù)報(bào)業(yè)務(wù)中,估計(jì)最大波高主要依靠預(yù)報(bào)員的主觀判斷,缺少客觀定量的預(yù)報(bào)方法。文圣常等[8]從理論上分析了平均波高與其他波高的關(guān)系,由瑞利分布導(dǎo)出了不同波高的比值,根據(jù)他的計(jì)算,最大波高約為有效波高的1.67倍,目前,對(duì)中國近海的最大波高和有效波高的基本特征已有一些研究,刁學(xué)賢[9]根據(jù)文圣常的理論結(jié)果對(duì)青島近岸開闊海區(qū)的波高變化進(jìn)行了統(tǒng)計(jì)分析,李合[10]針對(duì)北部灣海區(qū)的實(shí)際波浪狀況,利用測(cè)波資料計(jì)算了譜寬參數(shù)的分布規(guī)律,然后在寬譜的條件下確定了最大波高與有效波高的關(guān)系隨譜寬度參量與波個(gè)數(shù)的變化規(guī)律。國外的學(xué)者Goda[11]提出最大波高可以通過將有效波高乘以1.8倍來獲得,但是Kim等[12]研究指出,這種方法在沒有波破碎的情況下比較適用,在有波破碎的條件下,兩種波高的關(guān)系是高度非線性的,為此,他們引入了一個(gè)隨機(jī)變量來描述兩種波高的關(guān)系。針對(duì)海洋工程中關(guān)心的多年一遇最大設(shè)計(jì)波高和有效波高,國內(nèi)外也有很多研究[13-15],但是,這些最大設(shè)計(jì)波高的計(jì)算方法并不適用于日常海洋氣象業(yè)務(wù)中的最大波高預(yù)報(bào)。因此,分析上海沿海有效波高和最大波高的關(guān)系,建立計(jì)算最大波高的客觀方法,將能夠?yàn)檠睾8鞣N海洋經(jīng)濟(jì)和軍事活動(dòng)提供更加準(zhǔn)確的預(yù)報(bào)。

    2 資料和方法

    2.1 資料

    本文所使用的資料為:

    (1)上海沿海5個(gè)浮標(biāo)站的實(shí)測(cè)資料,分別為長江口船標(biāo)(122.53°E,31.1°N)、海礁浮標(biāo)(123.13°E,30.72°N)、黃澤洋船標(biāo)(122.53°E,30.5°N)、口外浮標(biāo)(122.81°E,30.99°N)和南漕船標(biāo)(122.52°E,30.98°N)。上海沿海的水深和浮標(biāo)位置見圖1。觀測(cè)資料時(shí)間段為2011年1月1日—2015年11月20日,時(shí)間間隔為1 h;

    (2)上海臺(tái)風(fēng)研究所WAVEWATCH III(WW3)海浪模式預(yù)報(bào)有效波高。數(shù)據(jù)時(shí)間段為2014年9月24日—2015年11月20日。模式的區(qū)域?yàn)?°—45°N,105°—145°E,空間分辨率為0.1°×0.1°,空間每點(diǎn)離散化波浪譜的方向分辨率為15°,即24個(gè)方向;時(shí)間分辨率為6 h,每天進(jìn)行兩次預(yù)報(bào),分別為08時(shí)和20時(shí)(北京時(shí),下同),預(yù)報(bào)時(shí)效168 h;所用風(fēng)場(chǎng)是來自美國國家環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction,NCEP)的AVN(Aviation)全球大氣環(huán)流業(yè)務(wù)模式的分析場(chǎng)和預(yù)報(bào)場(chǎng),在輸入模式前還對(duì)模式風(fēng)場(chǎng)的系統(tǒng)性誤差做了訂正,以提高預(yù)報(bào)質(zhì)量。在計(jì)算站點(diǎn)預(yù)報(bào)時(shí),使用雙線性插值方法將格點(diǎn)預(yù)報(bào)插值到站點(diǎn)。

    2.2 回歸分析方法

    設(shè)y為最大波高,x為有效波高,x、y的關(guān)系由一元n次回歸方程確定,即:

    圖1 上海沿海水深及浮標(biāo)位置分布圖

    為了分析不同回歸方程的效果,本文分別建立了x和y的一元一次和一元二次回歸方程。

    2.3 檢驗(yàn)統(tǒng)計(jì)量

    本文使用的檢驗(yàn)統(tǒng)計(jì)參數(shù)包括:相關(guān)系數(shù)、平均誤差、平均絕對(duì)誤差、最大絕對(duì)誤差和均方根誤差。

    (1)相關(guān)系數(shù):

    式(2)—(6)中:Pi為預(yù)報(bào)值,Pˉ為預(yù)報(bào)平均值,Oi為觀測(cè)值,Oˉ為觀測(cè)平均值,N為樣本個(gè)數(shù)。

    3 回歸方程的建立和試報(bào)分析

    表1 各站x和y的相關(guān)系數(shù)

    3.1 單站回歸方程建立

    各單站資料實(shí)測(cè)有效波高x和最大波高y的相關(guān)系數(shù)見表1。從表1中可以看到,5個(gè)站點(diǎn)的最大波高和有效波高的相關(guān)系數(shù)均在0.98以上,通過了信度0.01的相關(guān)性檢驗(yàn)(t0.01=0.24),這說明兩種波高具有明顯的線性相關(guān)性。

    以海礁浮標(biāo)為例,對(duì)2011年1月1日—2013年12月31日的有效波高和最大波高觀測(cè)資料進(jìn)行回歸分析,得到一元一次和一元二次回歸方程為:

    這與文圣常等[8]的理論結(jié)果接近。利用回歸方程對(duì)海礁浮標(biāo)2014年資料進(jìn)行試報(bào)檢驗(yàn),結(jié)果見表2,表中正偏差代表預(yù)報(bào)值比觀測(cè)值大,負(fù)偏差代表預(yù)報(bào)值比觀測(cè)值小。

    從表2中可以看到,對(duì)海礁浮標(biāo),建立的一元一次和一元二次回歸方程得到的最大波高與實(shí)測(cè)最大波高的相關(guān)系數(shù)都為0.988 0,通過了信度0.01的相關(guān)性檢驗(yàn),平均絕對(duì)誤差為0.1 m,說明兩個(gè)回歸方程的擬合效果都非常理想,同時(shí),也有一些預(yù)報(bào)偏差較大的情況,最大絕對(duì)誤差為2.0 m。此外,兩個(gè)回歸方程擬合最大波高的正偏差樣本都多于負(fù)偏差樣本,說明預(yù)報(bào)偏大的情況更多。在實(shí)際業(yè)務(wù)工作中,僅使用一元一次回歸方程就可以取得比較好的預(yù)報(bào)效果。其他4個(gè)浮標(biāo)站的一元一次回歸方程分別為:

    它們的擬合最大波高與實(shí)測(cè)最大波高的相關(guān)系數(shù)分別為0.986 6、0.985 7、0.989 3和0.987 4。對(duì)這4個(gè)浮標(biāo)站,一元一次回歸方程也很好地反映了最大波高和有效波高的關(guān)系。

    3.2 通用回歸方程建立

    利用口外浮標(biāo)、海礁浮標(biāo)、黃澤洋船標(biāo)、南漕船標(biāo)4個(gè)站2011—2013年的有效波高和最大波高資料建立上海沿海最大波高的通用一元一次回歸方程,結(jié)果如下:

    利用建立的最大波高通用回歸方程對(duì)長江口船標(biāo)站進(jìn)行試報(bào),檢驗(yàn)結(jié)果見表3。

    可以看到,利用4個(gè)站建立的通用一元一次回歸方程,其系數(shù)a0和a1與單站一元一次回歸方程的系數(shù)數(shù)值十分接近。對(duì)長江口船標(biāo)進(jìn)行試報(bào),結(jié)果表明,利用通用回歸方程擬合的最大波高與實(shí)測(cè)最大波高相關(guān)系數(shù)為0.988 2,通過了信度0.01的相關(guān)性檢驗(yàn),平均絕對(duì)誤差為0.1 m。這說明,在實(shí)際預(yù)報(bào)業(yè)務(wù)中,對(duì)上海沿海最大波高進(jìn)行估計(jì)時(shí),使用通用回歸方程即可達(dá)到較好的效果。

    表2 海礁浮標(biāo)實(shí)測(cè)最大波高y和回歸分析最大波高Y的有關(guān)統(tǒng)計(jì)量(樣本數(shù)N=7 595)

    表3 最大波高通用一元一次回歸方程對(duì)長江口船標(biāo)資料的試報(bào)檢驗(yàn)結(jié)果

    4 數(shù)值模式預(yù)報(bào)應(yīng)用

    通過以上的分析,已經(jīng)得到估計(jì)最大波高的客觀方程,在實(shí)際業(yè)務(wù)中,只需要使用數(shù)值預(yù)報(bào)的有效波高即可進(jìn)行最大波高的預(yù)報(bào)。數(shù)值模式預(yù)報(bào)與觀測(cè)相比,會(huì)存在誤差,在使用數(shù)值預(yù)報(bào)的結(jié)果前,需要進(jìn)行檢驗(yàn)。利用5個(gè)浮標(biāo)站2014年9月24日—2015年11月20日的觀測(cè)數(shù)據(jù)對(duì)WW3海浪模式的24 h有效波高預(yù)報(bào)進(jìn)行檢驗(yàn),結(jié)果見表4。

    從表4中可以看到,WW3海浪模式對(duì)5個(gè)浮標(biāo)站的有效波高預(yù)報(bào)都比較準(zhǔn)確,相關(guān)系數(shù)為0.788 6—0.915 4,均通過了信度0.01的相關(guān)性檢驗(yàn)。平均誤差最小為-0.04 m,最大也僅為-0.12 m,而平均絕對(duì)誤差數(shù)值也較小,為0.2—0.3 m。從表4中也可以看到,除黃澤洋船標(biāo)外,其他4個(gè)浮標(biāo)站的預(yù)報(bào)以正偏差為主,這可能與黃澤洋船標(biāo)附近島嶼較多、地形更復(fù)雜有關(guān)。此外,有4個(gè)浮標(biāo)站的預(yù)報(bào)以正偏差為主,而5個(gè)浮標(biāo)站的平均誤差均為負(fù)值,說明平均而言,盡管正偏差樣本占多數(shù),負(fù)偏差的絕對(duì)值仍要大于正偏差的絕對(duì)值。

    從表5中可以看到,利用WW3模式的有效波高進(jìn)行最大波高預(yù)報(bào),預(yù)報(bào)最大波高與實(shí)測(cè)最大波高具有較好的相關(guān)性,相關(guān)系數(shù)在0.778 2—0.911 2,均通過了信度0.01的相關(guān)性檢驗(yàn)。平均誤差最大為-0.17 m,最小為-0.05 m。5個(gè)浮標(biāo)站平均絕對(duì)誤差都為0.4 m,對(duì)于一些較大的最大波高值,預(yù)報(bào)結(jié)果偏差較大,最大絕對(duì)誤差出現(xiàn)在海礁浮標(biāo),為3.3 m。與有效波高的預(yù)報(bào)結(jié)果一致,除黃澤洋船標(biāo)外,其他4個(gè)浮標(biāo)站的預(yù)報(bào)以正偏差為主。此外,平均而言,最大波高預(yù)報(bào)的負(fù)偏差的絕對(duì)值也要大于正偏差的絕對(duì)值。

    表4 WW3模式24 h有效波高預(yù)報(bào)檢驗(yàn)

    表5 實(shí)測(cè)最大波高與預(yù)報(bào)最大波高誤差分析

    5 結(jié)論

    通過對(duì)上海沿海有效波高和最大波高關(guān)系進(jìn)行統(tǒng)計(jì)分析,并利用WW3海浪模式進(jìn)行試報(bào),本文主要得出以下幾點(diǎn)結(jié)論:

    (1)對(duì)海礁浮標(biāo),建立的一元一次和一元二次回歸方程擬合的最大波高與實(shí)測(cè)最大波高非常接近,平均絕對(duì)誤差為0.1 m,兩個(gè)回歸方程的擬合效果都非常理想。在實(shí)際業(yè)務(wù)工作中,僅使用一元一次回歸方程就可以取得比較好的預(yù)報(bào)效果。對(duì)其他4個(gè)浮標(biāo)站,一元一次回歸方程也很好地反映了最大波高和有效波高的關(guān)系;

    (2)利用口外浮標(biāo)、海礁浮標(biāo)、黃澤洋船標(biāo)、南漕船標(biāo)4個(gè)站建立的通用一元一次回歸方程,其系數(shù)與單站一元一次回歸方程的系數(shù)數(shù)值十分接近。利用通用方程對(duì)長江口船標(biāo)進(jìn)行試報(bào),擬合的最大波高與實(shí)測(cè)最大波高相關(guān)系數(shù)為0.988 2,平均絕對(duì)誤差為0.1 m。在實(shí)際預(yù)報(bào)業(yè)務(wù)中,對(duì)上海沿海最大波高進(jìn)行估計(jì)時(shí),使用通用回歸方程即可達(dá)到較好的效果;

    (3)WW3海浪模式對(duì)5個(gè)浮標(biāo)站的有效波高預(yù)報(bào)都比較準(zhǔn)確,利用WW3海浪模式的有效波高進(jìn)行最大波高預(yù)報(bào),預(yù)報(bào)最大波高與實(shí)測(cè)最大波高具有較好的相關(guān)性,平均誤差最大為-0.17 m,最小為-0.05 m,平均絕對(duì)誤差都為0.4 m,預(yù)報(bào)效果較好。

    [1]鄭崇偉,游小寶,潘靜,等.釣魚島、黃巖島海域風(fēng)能及波浪能開發(fā)環(huán)境分析[J].海洋預(yù)報(bào),2014,31(1):49-57.

    [2]鄭崇偉,李崇銀.中國南海島礁建設(shè):重點(diǎn)島礁的風(fēng)候、波候特征分析[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,45(9):1-6.

    [3]鄭崇偉,李崇銀.中國南海島礁建設(shè):風(fēng)力發(fā)電、海浪發(fā)電[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,46(9):7-14.

    [4]聞斌,劉俊.海浪數(shù)值模式研究回顧與進(jìn)展[J].海洋預(yù)報(bào),2006, 23(4):76-81.

    [5]鄭沛楠,宋軍,張芳苒,等.常用海洋數(shù)值模式簡(jiǎn)介[J].海洋預(yù)報(bào), 2008,25(4):108-120.

    [6]李本霞,吳淑萍,邢闖,等.近海近岸高精度海浪業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)[J].海洋預(yù)報(bào),2010,27(5):1-6.

    [7]孫明華.一個(gè)集合海浪預(yù)報(bào)系統(tǒng)及其初步試驗(yàn)[J].應(yīng)用氣象學(xué)報(bào),2011,22(6):673-680.

    [8]文圣常,宇宙文.海浪理論和計(jì)算原理[M].北京:科學(xué)出版社, 1985:177-184.

    [9]刁學(xué)賢.青島近岸波高的統(tǒng)計(jì)分析[J].海洋通報(bào),1998,17(4): 19-28.

    [10]李合.北部灣海區(qū)有效波高與最大波高的關(guān)系及最大設(shè)計(jì)波高[J].海洋預(yù)報(bào),1989,6(1):46-49.

    [11]Goda Y.A new Method of Wave Pressure Calculation for the Design of Composite Breakwater[R].PARI Report 012-03-02, 1973.

    [12]Kim D H,Kim T.Probabilistic Distribution of the Maximum Wave Height[J].China Ocean Engineering,2014,28(5):579-586.

    [13]候瑞科.利用正態(tài)分布計(jì)算年最大波高[J].海洋預(yù)報(bào),1991,8 (2):70-75.

    [14]尉華,沈正.海浪日最大波高序列的一種標(biāo)度性質(zhì)[J].海洋湖沼通報(bào),1994(4):1-8.

    [15]曹兵,徐常三,高清清,等.連云港海域極值波高統(tǒng)計(jì)分析方法研究[J].海洋通報(bào),2014,33(1):27-35.

    Analysis of the relationship between significant wave height and maximum wave height in Shanghai coastal area

    ZHU Zhi-hui,SHEN Qi-yan,CHANG Ya-Nan
    (Shanghai Marine Meteorological Center,Shanghai 201306 China)

    Using linear regression method,the relationship between significant wave height and maximum wave height in Shanghai costal area is analyzed based on five buoys data.It is found that the two kinds of wave height has a good linear relationship.Through the establishment of regression equation of the two kinds of wave,the objective method of estimating maximum wave height is obtained.For the 5 buoys,the one-dimensional linear regression equation can well reflect the relationship between significant wave height and maximum wave height. The general forecasting equation of maximum wave height has a good performance for Changjiangkou buoy and the average absolute error is only 0.1m.The forecasted significant wave height byWW3 wave model is close to the observation.Using the significant wave height from this model to forecast maximum wave height,the forecast value has a good correlation with the observation.The averaged error at the five buoys is-0.17m,the minimum is-0.05m,the averaged absolute error is 0.4m,and the result is reasonable.

    significant wave height;maximum wave height;regression analysis

    P731.22

    :A

    :1003-0239(2016)06-0069-05

    10.11737/j.issn.1003-0239.2016.06.008

    2016-03-28

    上海市氣象局面上項(xiàng)目(MS201409,MS201603)。

    朱智慧(1984-),男,高級(jí)工程師,碩士,從事海洋氣象預(yù)報(bào)工作。E-mail:zzh830830@163.com

    猜你喜歡
    波高浮標(biāo)海浪
    基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    受了委屈的浮標(biāo)君
    家教世界(2023年7期)2023-03-22 12:11:24
    受了委屈的浮標(biāo)君
    家教世界(2023年4期)2023-03-04 07:31:28
    丫丫和小海浪
    幼兒園(2021年13期)2021-12-02 05:13:54
    海浪
    小讀者(2021年2期)2021-11-23 07:17:34
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    非平整港池的多向不規(guī)則波試驗(yàn)研究
    樊應(yīng)舉
    書香兩岸(2020年3期)2020-06-29 12:33:45
    淫秽高清视频在线观看| 中国美女看黄片| 亚洲av第一区精品v没综合| 欧美色视频一区免费| 久久久国产成人精品二区| 免费一级毛片在线播放高清视频| 悠悠久久av| 国产91精品成人一区二区三区| 日韩av在线大香蕉| 成熟少妇高潮喷水视频| 亚洲va日本ⅴa欧美va伊人久久| 久久精品人妻少妇| 宅男免费午夜| 男人的好看免费观看在线视频| 午夜福利成人在线免费观看| 国产精品一区二区三区四区免费观看 | 美女大奶头视频| 露出奶头的视频| 久99久视频精品免费| 亚洲在线自拍视频| av在线蜜桃| 亚洲,欧美精品.| 精品久久久久久久末码| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 黄色成人免费大全| 少妇丰满av| 欧美日韩乱码在线| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 美女大奶头视频| 免费观看人在逋| 97超级碰碰碰精品色视频在线观看| av专区在线播放| 亚洲中文字幕一区二区三区有码在线看| 国产精品久久久久久精品电影| 91久久精品国产一区二区成人 | 亚洲中文字幕一区二区三区有码在线看| 免费电影在线观看免费观看| 精品欧美国产一区二区三| 真实男女啪啪啪动态图| 国产一区二区三区视频了| 久久99热这里只有精品18| 欧美zozozo另类| 国产美女午夜福利| 成人一区二区视频在线观看| 免费在线观看影片大全网站| 有码 亚洲区| av女优亚洲男人天堂| 久久国产精品人妻蜜桃| 国产免费一级a男人的天堂| 亚洲精品色激情综合| 午夜老司机福利剧场| 法律面前人人平等表现在哪些方面| 国产亚洲精品一区二区www| 国产精品影院久久| 18禁美女被吸乳视频| 成年免费大片在线观看| 色老头精品视频在线观看| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 在线免费观看不下载黄p国产 | xxxwww97欧美| 日本熟妇午夜| 国产精品久久久人人做人人爽| 在线观看美女被高潮喷水网站 | 最新美女视频免费是黄的| 午夜两性在线视频| 变态另类丝袜制服| 午夜激情欧美在线| 久久精品国产亚洲av涩爱 | 日韩av在线大香蕉| 国产高清videossex| 久久国产精品影院| 91麻豆精品激情在线观看国产| 高清毛片免费观看视频网站| 窝窝影院91人妻| 91久久精品电影网| 欧美一区二区亚洲| 色尼玛亚洲综合影院| 日本成人三级电影网站| 免费高清视频大片| 国产亚洲欧美98| 内射极品少妇av片p| 亚洲国产精品999在线| 51午夜福利影视在线观看| 中出人妻视频一区二区| 久久国产精品影院| 亚洲精品粉嫩美女一区| 日本黄色片子视频| 日本成人三级电影网站| 国产私拍福利视频在线观看| 国产成人aa在线观看| 在线国产一区二区在线| 3wmmmm亚洲av在线观看| 亚洲黑人精品在线| 中文字幕久久专区| 久久国产乱子伦精品免费另类| 搡老妇女老女人老熟妇| 99久久成人亚洲精品观看| 亚洲人成网站在线播放欧美日韩| 国产一区二区亚洲精品在线观看| 午夜激情福利司机影院| 一个人看视频在线观看www免费 | 国产色婷婷99| 亚洲自拍偷在线| а√天堂www在线а√下载| 成人亚洲精品av一区二区| 国产精品99久久99久久久不卡| 日韩欧美一区二区三区在线观看| 亚洲人成电影免费在线| 国产午夜精品论理片| 18禁黄网站禁片免费观看直播| 亚洲国产精品成人综合色| e午夜精品久久久久久久| 国产亚洲精品一区二区www| 午夜福利欧美成人| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品免费一区二区三区在线| x7x7x7水蜜桃| 国产精品1区2区在线观看.| 高清日韩中文字幕在线| 免费在线观看亚洲国产| 婷婷精品国产亚洲av| av欧美777| 成年女人毛片免费观看观看9| 欧美激情久久久久久爽电影| 99热只有精品国产| 桃色一区二区三区在线观看| 国产毛片a区久久久久| 一进一出好大好爽视频| 在线看三级毛片| 亚洲欧美日韩东京热| 99久久精品热视频| 国内精品久久久久久久电影| 国产精品98久久久久久宅男小说| 国产亚洲欧美在线一区二区| aaaaa片日本免费| 国产免费男女视频| 国产蜜桃级精品一区二区三区| 欧美乱妇无乱码| 午夜视频国产福利| 搡女人真爽免费视频火全软件 | 在线观看免费视频日本深夜| 国产精品爽爽va在线观看网站| 久久久国产成人精品二区| 老熟妇乱子伦视频在线观看| 天堂√8在线中文| x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 成年女人永久免费观看视频| 特大巨黑吊av在线直播| 午夜两性在线视频| 国产黄色小视频在线观看| 国产精品永久免费网站| 国产精品综合久久久久久久免费| 99国产精品一区二区蜜桃av| 免费观看精品视频网站| 精品国产三级普通话版| 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 丰满的人妻完整版| 美女高潮的动态| 亚洲最大成人中文| АⅤ资源中文在线天堂| 欧美日韩亚洲国产一区二区在线观看| 夜夜看夜夜爽夜夜摸| 国产亚洲欧美在线一区二区| 久久伊人香网站| av在线蜜桃| 俄罗斯特黄特色一大片| 亚洲人成网站高清观看| 精品人妻1区二区| 色播亚洲综合网| 看免费av毛片| 欧美在线黄色| 日韩有码中文字幕| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 99在线视频只有这里精品首页| 在线播放无遮挡| avwww免费| 男人和女人高潮做爰伦理| 欧美乱码精品一区二区三区| 黄色女人牲交| 最新在线观看一区二区三区| 小说图片视频综合网站| 一个人免费在线观看电影| 中文在线观看免费www的网站| 日韩免费av在线播放| 久久6这里有精品| 99热这里只有是精品50| 久久精品国产自在天天线| 久久久久国内视频| 内射极品少妇av片p| 国产精品一及| 成熟少妇高潮喷水视频| 欧美+亚洲+日韩+国产| 免费搜索国产男女视频| 日韩欧美在线乱码| 丰满的人妻完整版| 久久草成人影院| 动漫黄色视频在线观看| 国产三级黄色录像| 中文在线观看免费www的网站| 久久欧美精品欧美久久欧美| 91在线精品国自产拍蜜月 | 露出奶头的视频| 国内久久婷婷六月综合欲色啪| 12—13女人毛片做爰片一| 两个人视频免费观看高清| 国产探花在线观看一区二区| 精品久久久久久,| 国产亚洲av嫩草精品影院| 91在线观看av| 男女床上黄色一级片免费看| 成人亚洲精品av一区二区| 国产亚洲精品av在线| 久久婷婷人人爽人人干人人爱| 在线观看日韩欧美| 三级男女做爰猛烈吃奶摸视频| 国产精品三级大全| 成人18禁在线播放| 村上凉子中文字幕在线| 精品欧美国产一区二区三| 亚洲欧美日韩东京热| 男女床上黄色一级片免费看| 丰满人妻熟妇乱又伦精品不卡| 亚洲 国产 在线| 免费搜索国产男女视频| 国产免费男女视频| 中文字幕熟女人妻在线| 国内精品久久久久精免费| 国产一区在线观看成人免费| 日本撒尿小便嘘嘘汇集6| 小说图片视频综合网站| 亚洲aⅴ乱码一区二区在线播放| 国产av麻豆久久久久久久| 男女做爰动态图高潮gif福利片| 51国产日韩欧美| 亚洲精品一卡2卡三卡4卡5卡| 九色成人免费人妻av| 校园春色视频在线观看| 村上凉子中文字幕在线| av女优亚洲男人天堂| 老司机在亚洲福利影院| 久久久精品欧美日韩精品| 黑人欧美特级aaaaaa片| 九九久久精品国产亚洲av麻豆| 少妇人妻精品综合一区二区 | 午夜福利在线在线| 一本久久中文字幕| 久久精品国产亚洲av涩爱 | aaaaa片日本免费| 岛国在线观看网站| 18禁在线播放成人免费| 日日摸夜夜添夜夜添小说| 色综合亚洲欧美另类图片| 欧美黄色淫秽网站| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 在线十欧美十亚洲十日本专区| 亚洲精品色激情综合| 非洲黑人性xxxx精品又粗又长| 欧美精品啪啪一区二区三区| 美女被艹到高潮喷水动态| 在线观看舔阴道视频| 日韩欧美 国产精品| 国产精品av视频在线免费观看| 丝袜美腿在线中文| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 热99在线观看视频| 少妇人妻精品综合一区二区 | 波多野结衣巨乳人妻| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女| 成人av在线播放网站| 亚洲国产精品合色在线| 国产成年人精品一区二区| 少妇的丰满在线观看| 亚洲成人免费电影在线观看| 黄色视频,在线免费观看| 精品日产1卡2卡| 欧美不卡视频在线免费观看| 日本 av在线| 欧美成人性av电影在线观看| 亚洲人成网站在线播| 老汉色∧v一级毛片| 51午夜福利影视在线观看| 99热这里只有精品一区| 中文亚洲av片在线观看爽| 国产成人福利小说| 国产精品亚洲av一区麻豆| 国产精品99久久99久久久不卡| 国产精品久久久久久精品电影| 国模一区二区三区四区视频| 中出人妻视频一区二区| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 51午夜福利影视在线观看| 成人亚洲精品av一区二区| 日韩精品中文字幕看吧| 3wmmmm亚洲av在线观看| 久久亚洲真实| 欧洲精品卡2卡3卡4卡5卡区| 伊人久久精品亚洲午夜| 色综合亚洲欧美另类图片| 国产真实乱freesex| 中文资源天堂在线| 女人高潮潮喷娇喘18禁视频| 亚洲七黄色美女视频| 国产av一区在线观看免费| 精品一区二区三区视频在线观看免费| 精品久久久久久久毛片微露脸| 亚洲欧美精品综合久久99| 人人妻,人人澡人人爽秒播| 麻豆国产av国片精品| 熟女电影av网| 免费观看的影片在线观看| h日本视频在线播放| 国产精品爽爽va在线观看网站| 日日夜夜操网爽| 美女黄网站色视频| 久久人妻av系列| 精品国产三级普通话版| 校园春色视频在线观看| 日本 av在线| 欧美日韩亚洲国产一区二区在线观看| av国产免费在线观看| 国产日本99.免费观看| 中文字幕av成人在线电影| 一本一本综合久久| 动漫黄色视频在线观看| 毛片女人毛片| 国产99白浆流出| 久久久久国内视频| h日本视频在线播放| 69av精品久久久久久| 国产精品乱码一区二三区的特点| 男女之事视频高清在线观看| 亚洲最大成人中文| 免费搜索国产男女视频| av福利片在线观看| 嫁个100分男人电影在线观看| 国产黄片美女视频| 亚洲精华国产精华精| 日韩精品中文字幕看吧| 老鸭窝网址在线观看| 成人国产综合亚洲| av视频在线观看入口| 亚洲最大成人中文| 母亲3免费完整高清在线观看| 又黄又爽又免费观看的视频| 欧美黑人巨大hd| 手机成人av网站| 丝袜美腿在线中文| 久久国产精品影院| 亚洲最大成人手机在线| 999久久久精品免费观看国产| 丝袜美腿在线中文| 久久国产精品影院| 青草久久国产| 91久久精品电影网| 免费在线观看亚洲国产| 亚洲欧美日韩无卡精品| 亚洲avbb在线观看| 国内揄拍国产精品人妻在线| 欧美另类亚洲清纯唯美| 一区福利在线观看| 亚洲国产中文字幕在线视频| 精品国产三级普通话版| 午夜久久久久精精品| 搡老岳熟女国产| 亚洲av不卡在线观看| 人人妻人人澡欧美一区二区| 99在线视频只有这里精品首页| 熟女电影av网| 午夜精品久久久久久毛片777| av在线蜜桃| 国产真实乱freesex| 午夜福利高清视频| 国产麻豆成人av免费视频| 此物有八面人人有两片| 国产69精品久久久久777片| 久久人人精品亚洲av| 少妇的丰满在线观看| 亚洲成人中文字幕在线播放| 日韩欧美在线二视频| 久久精品国产亚洲av香蕉五月| 久久伊人香网站| 一进一出好大好爽视频| 日韩大尺度精品在线看网址| 欧美日韩中文字幕国产精品一区二区三区| 91九色精品人成在线观看| 免费观看人在逋| 亚洲最大成人中文| 国模一区二区三区四区视频| 国内精品美女久久久久久| 1024手机看黄色片| 伊人久久大香线蕉亚洲五| 很黄的视频免费| 国产av在哪里看| 最新在线观看一区二区三区| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 久久99热这里只有精品18| 国产又黄又爽又无遮挡在线| 国产成人福利小说| 亚洲一区二区三区不卡视频| 69av精品久久久久久| 变态另类成人亚洲欧美熟女| 国产高清有码在线观看视频| 婷婷精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 麻豆成人av在线观看| www.www免费av| 99久久九九国产精品国产免费| 亚洲av熟女| 日本熟妇午夜| 亚洲精品久久国产高清桃花| 国产极品精品免费视频能看的| 亚洲国产色片| 国产精品久久久久久亚洲av鲁大| 夜夜夜夜夜久久久久| 精品久久久久久,| 婷婷丁香在线五月| 久久人人精品亚洲av| 三级毛片av免费| 国产午夜精品久久久久久一区二区三区 | 国产精品野战在线观看| 国产精品自产拍在线观看55亚洲| 少妇高潮的动态图| 嫁个100分男人电影在线观看| 人妻夜夜爽99麻豆av| 午夜激情福利司机影院| 欧美在线黄色| 男女那种视频在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲午夜理论影院| av视频在线观看入口| 日本成人三级电影网站| 日韩有码中文字幕| 搡老岳熟女国产| www.熟女人妻精品国产| 一二三四社区在线视频社区8| 午夜免费激情av| 99国产精品一区二区三区| 嫩草影院入口| 三级国产精品欧美在线观看| 少妇人妻精品综合一区二区 | 国产91精品成人一区二区三区| 久久人妻av系列| 成年女人永久免费观看视频| 伊人久久精品亚洲午夜| 午夜免费激情av| 国产爱豆传媒在线观看| 国产v大片淫在线免费观看| 一进一出好大好爽视频| a级一级毛片免费在线观看| 亚洲av一区综合| 国产精品98久久久久久宅男小说| 久久午夜亚洲精品久久| 男女午夜视频在线观看| 欧美中文综合在线视频| 长腿黑丝高跟| 一个人看视频在线观看www免费 | 午夜福利在线在线| 欧洲精品卡2卡3卡4卡5卡区| 99在线人妻在线中文字幕| 好看av亚洲va欧美ⅴa在| 99久久无色码亚洲精品果冻| 日本黄大片高清| x7x7x7水蜜桃| 久久这里只有精品中国| 成人一区二区视频在线观看| 一夜夜www| 麻豆成人av在线观看| 日韩有码中文字幕| 国产一区二区在线观看日韩 | 国产伦精品一区二区三区视频9 | 在线免费观看的www视频| 日日摸夜夜添夜夜添小说| 国产黄色小视频在线观看| 国产精品亚洲一级av第二区| 99热6这里只有精品| 中文字幕人成人乱码亚洲影| 精品一区二区三区av网在线观看| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 美女 人体艺术 gogo| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 听说在线观看完整版免费高清| 欧美日韩精品网址| av女优亚洲男人天堂| 亚洲国产精品999在线| 国产aⅴ精品一区二区三区波| 久久久国产精品麻豆| 99久久九九国产精品国产免费| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av在线| 欧美日韩亚洲国产一区二区在线观看| 十八禁网站免费在线| 国产精品99久久久久久久久| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 亚洲久久久久久中文字幕| 亚洲国产欧洲综合997久久,| 两性午夜刺激爽爽歪歪视频在线观看| 给我免费播放毛片高清在线观看| 午夜影院日韩av| 欧美大码av| 乱人视频在线观看| 九九在线视频观看精品| 国产日本99.免费观看| 国产精品久久久人人做人人爽| 真人做人爱边吃奶动态| 一个人观看的视频www高清免费观看| 每晚都被弄得嗷嗷叫到高潮| 精品日产1卡2卡| 国产av不卡久久| 在线十欧美十亚洲十日本专区| 成年人黄色毛片网站| 国产真实伦视频高清在线观看 | 欧美中文综合在线视频| 欧美一区二区国产精品久久精品| 国产精品98久久久久久宅男小说| 国产成人av教育| avwww免费| 国产精品野战在线观看| 天天躁日日操中文字幕| 亚洲 国产 在线| 国产精品 欧美亚洲| 欧美性猛交黑人性爽| 精品福利观看| 女警被强在线播放| 久久久成人免费电影| 人人妻人人看人人澡| 欧美成人a在线观看| 97碰自拍视频| 久久久久精品国产欧美久久久| 琪琪午夜伦伦电影理论片6080| 欧美日韩精品网址| 中文字幕精品亚洲无线码一区| www.色视频.com| 亚洲午夜理论影院| 日本一二三区视频观看| 亚洲国产精品999在线| 亚洲欧美日韩高清在线视频| 女人十人毛片免费观看3o分钟| 国产爱豆传媒在线观看| 国产精品永久免费网站| 国产精品精品国产色婷婷| 在线观看免费视频日本深夜| 最新在线观看一区二区三区| 亚洲18禁久久av| 亚洲国产精品合色在线| 可以在线观看的亚洲视频| 最新中文字幕久久久久| 中文字幕av在线有码专区| www.www免费av| 无遮挡黄片免费观看| 亚洲欧美日韩高清在线视频| 亚洲aⅴ乱码一区二区在线播放| 精品电影一区二区在线| 精品久久久久久久久久久久久| 极品教师在线免费播放| 国产一区二区三区视频了| 亚洲av电影不卡..在线观看| 丰满乱子伦码专区| 午夜两性在线视频| 国产精品久久久久久久久免 | 欧美在线一区亚洲| 波多野结衣高清无吗| 天堂影院成人在线观看| 一个人免费在线观看电影| 成人av一区二区三区在线看| 久久伊人香网站| 欧美性感艳星| 国产一级毛片七仙女欲春2| 成人亚洲精品av一区二区| 久久久国产精品麻豆| 午夜两性在线视频| 91麻豆av在线| 欧美性猛交黑人性爽| 成人特级黄色片久久久久久久| 狂野欧美激情性xxxx| 亚洲国产欧美人成| 亚洲成人免费电影在线观看| 99久久综合精品五月天人人| 色综合站精品国产| 午夜两性在线视频| 精品日产1卡2卡| or卡值多少钱| 国内精品美女久久久久久| 内地一区二区视频在线| 国产精品自产拍在线观看55亚洲| 国产一区二区三区视频了| 亚洲精品乱码久久久v下载方式 | 一级黄色大片毛片| 女人被狂操c到高潮| 成年女人永久免费观看视频| 国产亚洲精品久久久久久毛片| 国产精品久久久久久久久免 | 亚洲五月天丁香| 色在线成人网| av专区在线播放| 午夜免费激情av| 男人舔女人下体高潮全视频| 91在线观看av| 久久久国产成人免费| 免费av毛片视频| 欧美绝顶高潮抽搐喷水| 热99re8久久精品国产|