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

    珠江口磨刀門水道鹽度變化與潮汐過程的相關(guān)性分析*

    2013-04-24 00:54:56劉丙軍陳曉宏辛彥博嚴(yán)淑蘭
    關(guān)鍵詞:咸潮潮差河口區(qū)

    章 文,劉丙軍,陳曉宏,辛彥博,嚴(yán)淑蘭

    (1.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院水資源與環(huán)境系∥華南地區(qū)水循環(huán)與水安全廣東省普通高校重點(diǎn)實(shí)驗(yàn)室,廣東 廣州 510275;2.珠江水利委員會水文局,廣東 廣州 510611)

    珠江三角洲地區(qū)具有三江匯流、八門入海的復(fù)雜河網(wǎng),受上游淡水徑流、河口區(qū)潮汐動力以及河口地形、風(fēng)力風(fēng)向、海平面變化、典型人類活動(河道采砂、水利工程調(diào)度)等多重因素影響,河口區(qū)水動力條件存在多種不確定性,鹽度變化過程異常復(fù)雜。同一水道或同一斷面的水體鹽度,會因季節(jié)(上游來水多寡)、潮別(漲或落)、時(shí)程(潮汐漲退過程)、河道與口門地形變化、海平面上升及風(fēng)向風(fēng)力等因素影響而改變[1]。

    近年來,不少學(xué)者對咸潮上溯期河口鹽度變化過程進(jìn)行了大量研究,并取得了諸多創(chuàng)新性成果[2-4]。在珠江河口區(qū),鹽度變化研究進(jìn)展主要集中在鹽度變化驅(qū)動機(jī)理與預(yù)報(bào)上。鹽度變化驅(qū)動機(jī)理研究方面,多以單個(gè)測站分析為主,主要從上游徑流、潮汐動力、河床演變、風(fēng)浪或流速等方面研究各要素與咸潮上溯期鹽度變化的關(guān)系。聞平等[5]運(yùn)用徑流-咸潮響應(yīng)模型,以馬口與三水流量為參照提出了磨刀門水道最小壓咸流量;包蕓等[6]認(rèn)為造成磨刀門水道鹽水強(qiáng)烈上溯的動力機(jī)制與枯季小潮期間特殊的潮汐動力狀態(tài)有關(guān);鹽度變化預(yù)測研究主要以數(shù)學(xué)模擬為主,多采用經(jīng)驗(yàn)?zāi)J?。陳水森等[7]應(yīng)用Brockway一維定常模式建立了磨刀門水道咸潮入侵的概化模型。

    河口區(qū)鹽度變化具有多時(shí)間尺度變化特征。陳榮力等[8]研究表明,珠江河口區(qū)鹽度變化與潮位變化過程密切相關(guān);薛建強(qiáng)[9]指出磨刀門水道咸界運(yùn)動落后于潮位變化,且不同周期的滯后情況有所不同。但目前研究尚未明確河口區(qū)鹽度變化與潮汐過程在時(shí)間上的相位關(guān)系?;诖耍疚睦媒徊嫘〔ǚ治龇椒ㄔ诜治鰰r(shí)間序列共振周期及位相關(guān)系上的優(yōu)勢,利用近十年枯水期磨刀門水道廣昌站逐日鹽度序列與三灶站同期潮位過程資料,研究了磨刀門水道鹽度變化和潮汐過程的相關(guān)性。

    1 研究區(qū)背景與數(shù)據(jù)

    磨刀門水道咸潮上溯期鹽度變化的影響因子眾多,潮汐過程是鹽度變化的主要驅(qū)動因子。本文在前人研究基礎(chǔ)上[10-11],以咸潮上溯較為嚴(yán)重的年份(2003-2004、2006-2007、2007-2008、2008-2009、2009-2010、2010-2011年)作為分析對象,選用磨刀門水道廣昌站枯水期鹽度序列及同期三灶站潮差序列,研究磨刀門水道鹽度變化與潮汐過程的相關(guān)性。

    鹽度觀測站廣昌泵站位于磨刀門水道下游河段,距離河口較近,能較好反映河口區(qū)潮汐過程對鹽度變化的影響[12];潮位觀測站選用三灶站,三灶站位于北緯22°02′,東經(jīng)113°24′,距離磨刀門16 km左右,是珠江河口區(qū)最重要驗(yàn)潮站之一,其潮位變化特征對研究珠江口潮位時(shí)間演變規(guī)律具有較好代表性[5]。磨刀門水道及廣昌站、三灶站位置示意圖見圖1。

    圖1 磨刀門水道示意圖Fig.1 Location of Modaomen waterway

    2 交叉小波分析方法

    傳統(tǒng)的相關(guān)分析如Pearson 相關(guān)分析、Spearman秩相關(guān)分析、Kendall相關(guān)分析等,在描述變量間的相互關(guān)系時(shí)對數(shù)據(jù)有一定要求,如變量需符合雙正態(tài)分布、有序分類等,且難以在時(shí)頻兩域表征不同信號間的相互關(guān)系及局部化特征。交叉小波變換可以克服常規(guī)的相關(guān)分析方法在線性相關(guān)條件下才適用的限制,從多時(shí)間尺度的角度研究揭示兩相關(guān)變量在不同周期上的相關(guān)關(guān)系,以及要素之間同步性問題[13-15]??梢姡徊嫘〔ǚ治龇椒ㄟm宜于鹽度變化與潮汐過程的時(shí)延相關(guān)特征和時(shí)頻位相關(guān)系特征分析。

    2.1 交叉小波變換(XWT)

    (1)

    (2)

    式(2)中,σX和σY分別為兩個(gè)序列的標(biāo)準(zhǔn)差,ZV(p)是關(guān)于p的置信度,v為Morlet小波變換的自由度,如在顯著水平為a=0.05時(shí),ZV(0.95)=3.99。若通過顯著性0.05的紅噪音檢驗(yàn),即認(rèn)為兩序列有較好的相關(guān)性,否則相關(guān)性不顯著。

    2.2 交叉小波位相

    (3)

    通過位相差箭頭方向可以判斷2個(gè)時(shí)間序列各尺度成分間的時(shí)滯相關(guān)性。

    2.3 小波相干譜(WTC)

    兩時(shí)間序列的小波相干譜計(jì)算公式為:

    (4)

    式(4)中,S(W)=Sscale(Stime(Wn(s)))是平滑器,其中,Sscale表示沿著小波伸縮尺度軸平滑,Stime表示沿著小波時(shí)間平移軸平滑。Morlet小波的平滑器表達(dá)式如下:

    Stime(W)S= (Wn(s)×c1-t2/(2s2))s

    (5)

    Sscale(W)n=Wn(s)×c2∏(0.6s)N

    (6)

    式(6)中,∏為一矩陣函數(shù),c1和c2為標(biāo)準(zhǔn)化常數(shù),0.6為經(jīng)驗(yàn)尺度。

    交叉小波變換可以揭示兩個(gè)時(shí)間序列共同的高能量區(qū)以及位相關(guān)系,小波相干譜則可以用來度量兩個(gè)時(shí)間序列在時(shí)頻空間上的局部相關(guān)程度,即使對應(yīng)交叉小波功率譜中的低能量值區(qū),在小波相干譜中也有可能很顯著。

    3 結(jié)果分析

    選用廣昌站近年來枯水期鹽度序列及同期三灶站潮差序列進(jìn)行交叉小波變換,研究廣昌站鹽度變化和與三灶站潮汐過程的相關(guān)性特征,結(jié)果見圖2(交叉小波功率譜)和圖3(小波相干譜圖)。

    圖2 廣昌站鹽度與三灶站潮差的交叉小波功率譜(XWT)Fig.2 The cross-wavelet power for the salinity sequence and the tidal range (XWT)

    1) 磨刀門水道鹽度變化與潮汐過程具有較好的時(shí)延相關(guān)性。各年份枯水期鹽度序列與潮差序列的交叉小波功率譜(XWT)(圖2和表1)在時(shí)頻域上顯示出相似的特征,即在15±2 d的時(shí)間尺度上,廣昌站鹽度過程與三灶站潮差過程出現(xiàn)了顯著的高能量區(qū),且通過顯著水平為0.05的紅噪音檢驗(yàn),表明兩者在15 d左右的周期尺度上顯著相關(guān)。鹽度與潮差的相關(guān)關(guān)系較為穩(wěn)定,一致表現(xiàn)為鹽度變化超前于潮差變化,超前位相約為60°~120°不等,說明鹽度變化超前于潮差變化1/6~1/3個(gè)周期(2.5~5 d)。這與聞平等[5]的研究結(jié)果,即在咸潮入侵的半月周期潮相變化中, 磨刀門水道含氯度日最大值并不出現(xiàn)在潮差最大值日, 而是提前 3~5 d 左右的結(jié)論一致。總體上看,15 d周期在時(shí)域的能量分布較均勻,在大潮期周圍較為顯著,在小潮期附近能量有所減弱。說明潮汐作用強(qiáng)時(shí)潮汐過程與鹽度變化相關(guān)性更為顯著。

    圖3 廣昌站鹽度與三灶站潮差的小波相干譜(WTC)Fig.3 The wavelet coherency spectrum between the salinity sequence and the tidal range (WTC) 圖中錐形細(xì)黑實(shí)線所包圍區(qū)域?yàn)椴皇苓吔缧?yīng)影響;粗黑實(shí)線內(nèi)為通過顯著水平為0.05的紅噪音檢驗(yàn);灰度表示相關(guān)程度的高低

    年份日期上游徑流量/(m3·s-1)潮差/m鹽度/(mg·L-1)位相關(guān)系/(°)2003-20041215-011624141755364900113-021222201726648602006-20071101-0101345017534571500102-0205280418343191002007-20081104-010129831775098900102-0214355117644121202008-20090102-021434491743438602009-20101101-010126341775382450102-021839651732913602010-20111101-12293374172210590

    鹽度變化之所以提前于潮差變化與河口地區(qū)的鹽度輸運(yùn)環(huán)流有關(guān)。在小潮期向大潮期轉(zhuǎn)換期間,潮差每日都在明顯增大,此時(shí)磨刀門口門外海水鹽度也比其他時(shí)期要高。也就是說,在小潮期轉(zhuǎn)大潮期時(shí),水平鹽度往往會較大,使得鹽度向上輸運(yùn)更明顯,導(dǎo)致鹽度變化提前于潮差變化。以2009-2010年為例,對比水平鹽度梯度與潮差變化(圖4)可知,三灶站鹽度最大值出現(xiàn)在11月15日及12月1日(5 055和6 616 mg/L),相應(yīng)的,水平鹽度梯度的最大值也出現(xiàn)在這兩天(797和966 mg/L),而三灶站潮差最大值分別出現(xiàn)在11月19日及12月4日(223和267 cm),即鹽度變化提前于潮差變化3 d左右。

    分析各年份交叉小波功率譜,鹽度變化與潮差過程的時(shí)滯關(guān)系各不相同,具體情況見表1。2003-2004、2010-2011年枯水期,鹽度變化和潮汐過程的交叉小波功率譜在15 d的周期上表現(xiàn)出高能量,兩者位相角為60°~90°; 2006-2007年枯水期,兩者位相關(guān)系為100°~150°;2007-2008年枯水期,兩者位相關(guān)系為90°~120°;2008-2009、2009-2010年枯水期,兩者位相關(guān)系為60°。可見,咸潮上溯過程中,對應(yīng)潮差越大的年份,對應(yīng)的鹽度與潮差之間的位相角越大,即鹽度變化提前潮差變化越多。如2006-2007及2007-2008年的咸潮上溯問題比較嚴(yán)重(潮差分別為1.75和1.75 m),相應(yīng)其位相角也越大,超過了2003-2004、2008-2009及2010-2011年(潮差分別為1.71、1.67和1.69 m)的位相角。對比上游同期徑流過程(三水與馬口合計(jì)流量),位相角的變化與上游同期徑流量的變化亦存在一定對應(yīng)關(guān)系,即當(dāng)上游徑流量由大到小變化時(shí),位相角也表現(xiàn)出由大到小變化的規(guī)律。如2003-2004及2006-2007年枯水期,上游同期徑流量變化情況均為由大到小,與之相對應(yīng),交叉譜的位相角也由大到小變化;2007-2008、2009-2010年枯水期的情況則相反,上游同期徑流量與位相角變化規(guī)律均為由小到大??梢?,徑流對鹽度變化存在逆向驅(qū)動作用,即上游徑流過程變大時(shí),河口區(qū)站點(diǎn)鹽度變小,鹽度變化提前于潮汐過程的時(shí)間越長;反之,上游徑流過程減小,河口區(qū)站點(diǎn)鹽度增大,鹽度變化提前于潮汐過程的時(shí)間則越短。這與劉杰斌等[3]研究成果一致,即上游徑流的大小會影響鹽水上溯速度和距離。

    圖4 水平鹽度梯度與三灶站潮差的關(guān)系Fig.4 Relationship between the horizontal gradient of salinity and tidal range in the MW

    2)磨刀門水道鹽度變化與潮汐過程在15 d的周期尺度上時(shí)延相關(guān)性較強(qiáng)。分析廣昌站鹽度與三灶站潮差的小波相干譜(WTC)(圖3),與交叉小波功率譜對應(yīng),以15 d為周期變化的高能量區(qū)表現(xiàn)明顯,表明在15 d左右的周期尺度上,鹽度序列與潮差序列有較強(qiáng)的時(shí)延相關(guān)性。為進(jìn)一步說明鹽度變化具有15d的周期特征,采用Morlet小波作為基函數(shù)進(jìn)行小波變換,研究磨刀門水道河口區(qū)的鹽度變化周期特征,結(jié)果見圖5??梢钥闯觯サ堕T水道鹽度變化具有與潮汐過程相似的周期變化特征,即鹽度變化存在多時(shí)間周期特征,其中主要周期成分為14.8 d(對于Morlet小波,周期T約為尺度因子a的1.22倍),這與聞平等[5]的研究成果一致,即磨刀門水道氯化物含量與潮汐半月周期有關(guān),具有與潮汐過程相同的半月周期變化特征。

    圖5 枯水期廣昌站鹽度序列小波系數(shù)圖Fig.5 The isograms of the wavelet transform coefficients of the salinity sequence in Guangchang Station

    4 結(jié) 論

    受上游淡水徑流、河口區(qū)潮汐動力、河口地形、風(fēng)力風(fēng)向、海平面變化、河道采砂、水利工程調(diào)度等多重因素影響,河口區(qū)鹽度變化過程異常復(fù)雜。論文選用近10年來珠江河口區(qū)磨刀門水道鹽度與三灶站枯水期同期潮差數(shù)據(jù),運(yùn)用交叉小波分析方法,研究了咸潮上溯期磨刀門水道鹽度周期特征及其與潮汐過程的相關(guān)性,得出如下結(jié)論:

    1) 磨刀門水道鹽度變化具有與潮汐過程相同的多時(shí)間尺度周期變化特征,主要周期成分為14.8 d;

    2) 鹽度變化與潮汐過程密切相關(guān),兩者在15 d左右的周期尺度上顯著相關(guān),鹽度變化超前于潮差變化2.5~3.75 d,且潮汐作用越強(qiáng),其相關(guān)性越顯著。

    3) 徑流對鹽度變化存在逆向驅(qū)動作用,河口區(qū)鹽度變化規(guī)律與上游同期徑流量的變化存在一定正相關(guān)關(guān)系,即徑流量由大到小變化,河口區(qū)鹽度由小變大,鹽度變化提前于潮汐過程的時(shí)間也表現(xiàn)為由大到小變化的規(guī)律。

    河口區(qū)鹽度變化具有較強(qiáng)不確定性與復(fù)雜性,本文僅分析了珠江河口區(qū)磨刀門水道鹽度變化與潮汐過程的相關(guān)性,其他影響因素如上游徑流、風(fēng)向、河口地形等對河口區(qū)水道鹽度變化的影響仍待進(jìn)一步研究和完善。

    參考文獻(xiàn):

    [1] 莫思平,李越,盧素蘭.廣州水道咸潮影響因素分析[J].水利水運(yùn)工程學(xué)報(bào),2007(4):36-42.

    [2] 劉杰斌,包蕓,黃宇銘.豐、枯水年磨刀門水道鹽水上溯運(yùn)動規(guī)律對比[J].力學(xué)學(xué)報(bào),2010,42(6):1098-1103.

    [3] BROCKWAY R, BOWERS D, HOGUANE A, et al. A note on salt intrusion in funnel-shaped estuaries: Application to the Incomati estuary, Mozambique[J]. Estuarine, Coastal and Shelf Science 2006, 66:1-5.

    [4] XU Hongzhou, LIN Jing, WANG Dongxiao. Numerical study on salinity stratification in the Pamlico River Estuary [J]. Estuarine, Coastal and Shelf Science 2008,80:74-84.

    [5] 聞平,陳曉宏,劉斌,等.磨刀門水道咸潮入侵及其變異分析[J].水文, 2007, 27(3): 65-67.

    [6] 包蕓,劉杰斌,任杰等.磨刀門水道鹽水強(qiáng)烈上溯規(guī)律和動力機(jī)制研究[J].中國科學(xué)G輯:物理學(xué) 力學(xué) 天文學(xué),2009,39(10):1527-1534.

    [7] 陳水森,方立剛,李宏麗等.珠江口咸潮入侵分析與經(jīng)驗(yàn)?zāi)P汀阅サ堕T水道為例[J].水科學(xué)進(jìn)展,2007,18(5): 751-755.

    [8] 陳榮力,劉誠,高時(shí)友.磨刀門水道枯季咸潮上溯規(guī)律分析[J].水動力學(xué)研究與進(jìn)展,2011,26(3):312-317.

    [9] 薛建強(qiáng).磨刀門水道咸潮運(yùn)動規(guī)律初步分析[J].廣西水利水電,2012(2):26-27.

    [10] 孔蘭,陳曉宏,聞平等.2009/2010年枯水期珠江口磨刀門水道強(qiáng)咸潮分析[J].自然資源學(xué)報(bào),2011,26(11):1858-1864.

    [11] 尹小玲,張紅武,方紅衛(wèi).枯季磨刀門水道咸潮活動與壓咸控制分析[J].水動力學(xué)研究與進(jìn)展A輯,2008,23(5):554-559.

    [12] 胡溪,毛獻(xiàn)忠.珠江口磨刀門水道咸潮入侵規(guī)律研究[J].水利學(xué)報(bào),2012,43(5):529-536

    [13] 王文圣,熊華康,丁晶.日流量預(yù)測的小波網(wǎng)絡(luò)模型初步探討[J].水科學(xué)進(jìn)展,2004,15(4):382-386.

    [14] 余丹丹,張韌,洪梅,等.基于交叉小波相干的西太平洋副高與東亞夏季風(fēng)系統(tǒng)的關(guān)聯(lián)性分析[J].南京氣象學(xué)院學(xué)報(bào),2007,30(6),755-769.

    [15] TORRENCE, CHRISTOPHER, GILBERT P. COMPO. A Practical guide to wavelet analysis[J]. Bulletin of the American Meteorological Society, 1998, 79: 61-78.

    猜你喜歡
    咸潮潮差河口區(qū)
    天文大潮+寒潮大風(fēng)為何會加劇咸潮
    漢江河口區(qū)四明灘段四明上護(hù)岸坍岸修復(fù)技術(shù)
    海平面變化背景下三大河口咸潮入侵特征及變化淺析
    長江口潮差中長期變化對河口生態(tài)環(huán)境的影響
    基于Copula函數(shù)的甬江流域設(shè)計(jì)潮位過程研究
    “傳承—信仰—未來”教育主題的設(shè)計(jì)者
    ——記山東省東營市河口區(qū)河安小學(xué)校長劉波
    教書育人(2017年23期)2017-09-26 06:47:50
    閩江下游及河口區(qū)表層沉積物多種微量元素特征
    珠江口咸潮影響因素分析
    論地質(zhì)鉆探技術(shù)在海水域工程的方法措施
    地球(2015年8期)2015-10-21 19:30:08
    珠江口咸潮影響因素分析
    a级毛片在线看网站| 免费看a级黄色片| 精品亚洲成国产av| 久久亚洲精品不卡| 亚洲av美国av| 狠狠婷婷综合久久久久久88av| 亚洲欧美一区二区三区黑人| 欧美成狂野欧美在线观看| 别揉我奶头~嗯~啊~动态视频| 桃花免费在线播放| 欧美日韩亚洲国产一区二区在线观看 | 又紧又爽又黄一区二区| 十八禁高潮呻吟视频| 巨乳人妻的诱惑在线观看| 飞空精品影院首页| 亚洲午夜精品一区,二区,三区| 欧美国产精品va在线观看不卡| 成年版毛片免费区| 美女国产高潮福利片在线看| 在线观看www视频免费| 久久久久久人人人人人| 欧美国产精品va在线观看不卡| av一本久久久久| 亚洲人成伊人成综合网2020| 亚洲国产欧美一区二区综合| 免费高清在线观看日韩| 国产免费av片在线观看野外av| 日韩欧美三级三区| 天堂动漫精品| 久久精品人人爽人人爽视色| 国产精品九九99| 国产日韩欧美在线精品| 丝瓜视频免费看黄片| 亚洲国产看品久久| 久久av网站| 国产淫语在线视频| 亚洲欧美色中文字幕在线| 久久狼人影院| 久久天堂一区二区三区四区| 日韩欧美三级三区| 午夜福利免费观看在线| 精品少妇黑人巨大在线播放| 亚洲中文字幕日韩| 成人国语在线视频| 国产成人精品无人区| 亚洲国产毛片av蜜桃av| 一区福利在线观看| 国产不卡av网站在线观看| 亚洲熟妇熟女久久| 欧美激情 高清一区二区三区| 丝袜喷水一区| 黑人巨大精品欧美一区二区mp4| 久久精品国产综合久久久| 国产老妇伦熟女老妇高清| 激情视频va一区二区三区| 女同久久另类99精品国产91| 建设人人有责人人尽责人人享有的| 女人被躁到高潮嗷嗷叫费观| 国产淫语在线视频| 日韩免费高清中文字幕av| 免费一级毛片在线播放高清视频 | 一区二区av电影网| 日韩一卡2卡3卡4卡2021年| 人人妻人人澡人人看| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机午夜十八禁免费视频| av国产精品久久久久影院| 欧美日韩亚洲国产一区二区在线观看 | 下体分泌物呈黄色| 每晚都被弄得嗷嗷叫到高潮| 欧美人与性动交α欧美软件| 免费少妇av软件| 午夜福利一区二区在线看| 高清视频免费观看一区二区| 五月开心婷婷网| 国产成人精品久久二区二区91| 99国产精品一区二区蜜桃av | 香蕉久久夜色| 男人舔女人的私密视频| 亚洲精品美女久久av网站| 黄片大片在线免费观看| 91老司机精品| 国产欧美日韩一区二区三区在线| 19禁男女啪啪无遮挡网站| 日韩大码丰满熟妇| 最新美女视频免费是黄的| 叶爱在线成人免费视频播放| 成人av一区二区三区在线看| 日日夜夜操网爽| 久久影院123| 亚洲avbb在线观看| 不卡一级毛片| 久久久久久亚洲精品国产蜜桃av| 狠狠狠狠99中文字幕| 日本精品一区二区三区蜜桃| 久热爱精品视频在线9| 亚洲欧洲日产国产| bbb黄色大片| 久久天躁狠狠躁夜夜2o2o| 超碰97精品在线观看| 黄色a级毛片大全视频| 日本vs欧美在线观看视频| 亚洲国产欧美一区二区综合| 欧美日本中文国产一区发布| 丰满少妇做爰视频| aaaaa片日本免费| 亚洲欧美一区二区三区黑人| 国产伦理片在线播放av一区| 无遮挡黄片免费观看| 欧美日韩精品网址| 久久 成人 亚洲| 老司机午夜十八禁免费视频| 69av精品久久久久久 | 丝袜人妻中文字幕| 欧美变态另类bdsm刘玥| 色婷婷久久久亚洲欧美| 免费久久久久久久精品成人欧美视频| 亚洲天堂av无毛| 日本精品一区二区三区蜜桃| 亚洲精品国产区一区二| 成年女人毛片免费观看观看9 | 亚洲中文字幕日韩| 亚洲国产欧美网| 成人影院久久| 欧美亚洲日本最大视频资源| 黄色视频,在线免费观看| 国产成人免费观看mmmm| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 国产区一区二久久| 777久久人妻少妇嫩草av网站| 欧美黑人精品巨大| 香蕉国产在线看| 高潮久久久久久久久久久不卡| 国产成人欧美| 国产精品av久久久久免费| 一本色道久久久久久精品综合| av不卡在线播放| 精品一区二区三卡| 国产黄频视频在线观看| 亚洲成人手机| 在线观看一区二区三区激情| a在线观看视频网站| 国产精品一区二区免费欧美| 国产精品影院久久| 国产精品 国内视频| 丰满迷人的少妇在线观看| 精品国产亚洲在线| 黄色视频,在线免费观看| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 一本一本久久a久久精品综合妖精| 久久香蕉激情| av不卡在线播放| 精品人妻在线不人妻| 18在线观看网站| 免费在线观看日本一区| 久久久国产欧美日韩av| a在线观看视频网站| 免费在线观看日本一区| 国产日韩欧美在线精品| 如日韩欧美国产精品一区二区三区| 国产精品电影一区二区三区 | 啦啦啦 在线观看视频| 三级毛片av免费| 国产精品久久久久成人av| 一区福利在线观看| 久9热在线精品视频| 亚洲成人手机| 亚洲国产看品久久| 亚洲色图 男人天堂 中文字幕| 中文字幕av电影在线播放| 久久人人97超碰香蕉20202| 日韩视频一区二区在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线 | 色在线成人网| 女性生殖器流出的白浆| 欧美日韩亚洲综合一区二区三区_| 久久午夜亚洲精品久久| 日本精品一区二区三区蜜桃| 亚洲中文av在线| 欧美成狂野欧美在线观看| 亚洲综合色网址| 少妇猛男粗大的猛烈进出视频| 亚洲欧美精品综合一区二区三区| 大型黄色视频在线免费观看| 另类亚洲欧美激情| 国产成人影院久久av| 亚洲国产欧美在线一区| 黑人操中国人逼视频| 久久久久久久国产电影| 亚洲精品美女久久久久99蜜臀| 日韩大码丰满熟妇| aaaaa片日本免费| 亚洲情色 制服丝袜| 99九九在线精品视频| 黄色视频在线播放观看不卡| 国产成人欧美| 日韩欧美一区视频在线观看| 一二三四社区在线视频社区8| 免费高清在线观看日韩| 99精品久久久久人妻精品| 亚洲精品国产色婷婷电影| 老熟妇仑乱视频hdxx| 亚洲av日韩在线播放| 精品一区二区三区av网在线观看 | 欧美亚洲 丝袜 人妻 在线| 色综合婷婷激情| 精品人妻在线不人妻| 无人区码免费观看不卡 | 亚洲天堂av无毛| 超碰成人久久| 免费高清在线观看日韩| xxxhd国产人妻xxx| 80岁老熟妇乱子伦牲交| 精品久久蜜臀av无| 亚洲 国产 在线| 精品福利永久在线观看| 咕卡用的链子| 欧美日本中文国产一区发布| 欧美大码av| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 咕卡用的链子| 国产亚洲欧美在线一区二区| 欧美久久黑人一区二区| 久久亚洲精品不卡| 日韩大片免费观看网站| 久久久久久久大尺度免费视频| 高清黄色对白视频在线免费看| 精品国产乱子伦一区二区三区| 国产一区有黄有色的免费视频| 国产av又大| a级毛片黄视频| 美女高潮到喷水免费观看| 国产激情久久老熟女| 三上悠亚av全集在线观看| 人人妻人人爽人人添夜夜欢视频| 18禁裸乳无遮挡动漫免费视频| 亚洲精品中文字幕在线视频| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 一本一本久久a久久精品综合妖精| 亚洲第一欧美日韩一区二区三区 | 日韩人妻精品一区2区三区| 男男h啪啪无遮挡| 国产精品国产av在线观看| 亚洲免费av在线视频| 久久久久久免费高清国产稀缺| 国产伦人伦偷精品视频| 午夜日韩欧美国产| 日韩三级视频一区二区三区| 高清在线国产一区| 国产又色又爽无遮挡免费看| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美一区二区三区黑人| 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 午夜福利视频在线观看免费| 日韩视频一区二区在线观看| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 在线观看人妻少妇| 国产午夜精品久久久久久| 亚洲免费av在线视频| 中文字幕精品免费在线观看视频| 波多野结衣av一区二区av| 精品免费久久久久久久清纯 | 国产成人精品久久二区二区91| 十分钟在线观看高清视频www| 色老头精品视频在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲视频免费观看视频| 久久这里只有精品19| 欧美国产精品va在线观看不卡| 两个人免费观看高清视频| 交换朋友夫妻互换小说| 夜夜夜夜夜久久久久| 亚洲av成人一区二区三| 亚洲少妇的诱惑av| 男女之事视频高清在线观看| 久久久久国产一级毛片高清牌| 99国产精品一区二区三区| 精品国产一区二区三区四区第35| 美女主播在线视频| 久久婷婷成人综合色麻豆| 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 国产有黄有色有爽视频| 黑人操中国人逼视频| 国产一区二区三区综合在线观看| 久久国产精品人妻蜜桃| 91麻豆精品激情在线观看国产 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产区一区二| 香蕉国产在线看| 在线观看66精品国产| 一夜夜www| 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀| 国产片内射在线| 在线观看一区二区三区激情| 可以免费在线观看a视频的电影网站| 自线自在国产av| 如日韩欧美国产精品一区二区三区| 国产精品99久久99久久久不卡| 中文字幕人妻丝袜制服| 无限看片的www在线观看| 90打野战视频偷拍视频| 精品高清国产在线一区| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩高清在线视频 | 热re99久久国产66热| 十八禁网站网址无遮挡| 一级毛片女人18水好多| 亚洲久久久国产精品| 高清av免费在线| tube8黄色片| 免费在线观看黄色视频的| 老司机亚洲免费影院| www.999成人在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品久久二区二区91| 免费久久久久久久精品成人欧美视频| 亚洲熟女毛片儿| 宅男免费午夜| 男女边摸边吃奶| 两个人看的免费小视频| 欧美日韩亚洲高清精品| av有码第一页| 日韩有码中文字幕| 国产老妇伦熟女老妇高清| 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩高清在线视频 | 男女高潮啪啪啪动态图| 在线永久观看黄色视频| 久久影院123| 另类亚洲欧美激情| 黑人猛操日本美女一级片| 亚洲成人免费av在线播放| 99热网站在线观看| 欧美av亚洲av综合av国产av| 欧美成人免费av一区二区三区 | 欧美激情 高清一区二区三区| 在线播放国产精品三级| 亚洲精品一二三| 在线播放国产精品三级| 亚洲第一青青草原| 黑人欧美特级aaaaaa片| 国产在线视频一区二区| 久久国产亚洲av麻豆专区| 国产精品九九99| 亚洲色图综合在线观看| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 美女高潮到喷水免费观看| 一级毛片精品| 国产在线免费精品| 亚洲成人手机| 国产有黄有色有爽视频| 岛国在线观看网站| 一级,二级,三级黄色视频| 欧美日韩精品网址| 十八禁人妻一区二区| 亚洲精品成人av观看孕妇| 久久免费观看电影| 亚洲欧美激情在线| 国产av又大| 悠悠久久av| 波多野结衣一区麻豆| 老汉色∧v一级毛片| 国产熟女午夜一区二区三区| 日韩免费av在线播放| 一进一出抽搐动态| 亚洲精品在线观看二区| 老司机午夜福利在线观看视频 | 日本vs欧美在线观看视频| 日韩制服丝袜自拍偷拍| 99热网站在线观看| 啦啦啦免费观看视频1| 欧美日韩亚洲国产一区二区在线观看 | 成年女人毛片免费观看观看9 | 亚洲男人天堂网一区| 国产一区二区 视频在线| 精品乱码久久久久久99久播| 看免费av毛片| 老司机影院毛片| 一区二区三区乱码不卡18| 欧美日韩亚洲综合一区二区三区_| 中文字幕精品免费在线观看视频| 精品福利永久在线观看| 亚洲三区欧美一区| 麻豆国产av国片精品| 女人被躁到高潮嗷嗷叫费观| aaaaa片日本免费| 久久久国产欧美日韩av| 国产成人啪精品午夜网站| 80岁老熟妇乱子伦牲交| 淫妇啪啪啪对白视频| 亚洲av美国av| 777米奇影视久久| 狠狠精品人妻久久久久久综合| 久久亚洲真实| 99久久99久久久精品蜜桃| 亚洲精品一二三| e午夜精品久久久久久久| 色老头精品视频在线观看| 亚洲精品粉嫩美女一区| 欧美国产精品一级二级三级| 久久久久久久国产电影| 久久天躁狠狠躁夜夜2o2o| 搡老熟女国产l中国老女人| 99精品欧美一区二区三区四区| 久久久久国产一级毛片高清牌| 欧美日韩亚洲高清精品| 色精品久久人妻99蜜桃| 天天躁日日躁夜夜躁夜夜| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 久久久久久久国产电影| 久久精品成人免费网站| 国产在视频线精品| 99国产精品免费福利视频| av天堂久久9| 成年动漫av网址| 亚洲第一欧美日韩一区二区三区 | 国产黄色免费在线视频| 老熟妇乱子伦视频在线观看| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 久久久水蜜桃国产精品网| 老司机午夜十八禁免费视频| svipshipincom国产片| 欧美日韩亚洲高清精品| 免费观看av网站的网址| 国产亚洲精品一区二区www | 少妇被粗大的猛进出69影院| 精品一区二区三区四区五区乱码| 一区二区三区激情视频| 天天操日日干夜夜撸| 丁香六月欧美| 91麻豆精品激情在线观看国产 | 一本一本久久a久久精品综合妖精| 日本wwww免费看| 一级a爱视频在线免费观看| 日韩欧美免费精品| 色视频在线一区二区三区| 午夜精品国产一区二区电影| 亚洲av日韩在线播放| 在线观看免费高清a一片| 色在线成人网| 欧美成狂野欧美在线观看| 麻豆乱淫一区二区| 纵有疾风起免费观看全集完整版| 天天躁日日躁夜夜躁夜夜| 91成年电影在线观看| 午夜福利视频精品| 51午夜福利影视在线观看| 侵犯人妻中文字幕一二三四区| 香蕉国产在线看| 亚洲伊人色综图| 婷婷丁香在线五月| 亚洲av美国av| 在线观看舔阴道视频| 国产精品久久久久久人妻精品电影 | 亚洲中文av在线| 高清欧美精品videossex| 中亚洲国语对白在线视频| 婷婷成人精品国产| 国产区一区二久久| av欧美777| 国产一区二区 视频在线| 性高湖久久久久久久久免费观看| 免费在线观看视频国产中文字幕亚洲| 十八禁人妻一区二区| 999久久久国产精品视频| 精品久久蜜臀av无| 人人妻人人爽人人添夜夜欢视频| 91九色精品人成在线观看| 欧美在线一区亚洲| 老司机靠b影院| 免费不卡黄色视频| 欧美变态另类bdsm刘玥| 99国产精品99久久久久| 亚洲黑人精品在线| 悠悠久久av| 桃花免费在线播放| 久久久欧美国产精品| 中文字幕av电影在线播放| 亚洲av第一区精品v没综合| 久久国产精品影院| 国产精品秋霞免费鲁丝片| 国产又色又爽无遮挡免费看| 青青草视频在线视频观看| 极品教师在线免费播放| 国产精品亚洲一级av第二区| 在线观看免费午夜福利视频| 日韩视频在线欧美| 久久青草综合色| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 免费少妇av软件| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美精品综合一区二区三区| 国产成人av教育| 精品一区二区三卡| 亚洲色图综合在线观看| 大型av网站在线播放| 亚洲精品中文字幕在线视频| 可以免费在线观看a视频的电影网站| 久久中文看片网| 亚洲成人手机| 中国美女看黄片| 纯流量卡能插随身wifi吗| 一级片'在线观看视频| 欧美激情 高清一区二区三区| 99久久99久久久精品蜜桃| 午夜福利在线观看吧| 91字幕亚洲| 2018国产大陆天天弄谢| 欧美激情极品国产一区二区三区| 欧美激情久久久久久爽电影 | 男女午夜视频在线观看| 亚洲精品国产色婷婷电影| 91成年电影在线观看| 老司机在亚洲福利影院| 欧美精品一区二区免费开放| 色在线成人网| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩亚洲综合一区二区三区_| 在线观看66精品国产| 一进一出抽搐动态| 考比视频在线观看| netflix在线观看网站| 啪啪无遮挡十八禁网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品国产色婷婷电影| 黑丝袜美女国产一区| 亚洲熟女毛片儿| 国产成人一区二区三区免费视频网站| 日韩 欧美 亚洲 中文字幕| 日日夜夜操网爽| 亚洲一码二码三码区别大吗| 老司机靠b影院| 久久久久国产一级毛片高清牌| 国产欧美日韩精品亚洲av| 亚洲一区中文字幕在线| 视频区图区小说| 成年动漫av网址| 中文字幕av电影在线播放| 老鸭窝网址在线观看| 日韩制服丝袜自拍偷拍| 欧美大码av| 性色av乱码一区二区三区2| 亚洲精品国产区一区二| 高清av免费在线| 日韩熟女老妇一区二区性免费视频| 手机成人av网站| 视频区图区小说| 久久天堂一区二区三区四区| 搡老熟女国产l中国老女人| 青草久久国产| 日本欧美视频一区| 中文字幕av电影在线播放| av超薄肉色丝袜交足视频| 99久久人妻综合| 精品第一国产精品| 久久精品国产亚洲av高清一级| 亚洲国产成人一精品久久久| 久久久精品94久久精品| 国内毛片毛片毛片毛片毛片| 国产精品久久久久久精品古装| 精品午夜福利视频在线观看一区 | 久久人妻av系列| 肉色欧美久久久久久久蜜桃| 国产成人欧美| 免费av中文字幕在线| 亚洲人成77777在线视频| 动漫黄色视频在线观看| 老司机福利观看| 91成年电影在线观看| 丰满饥渴人妻一区二区三| 欧美一级毛片孕妇| 久久久国产精品麻豆| 国产精品亚洲一级av第二区| 久久久久视频综合| 亚洲中文av在线| 日韩一区二区三区影片| 777米奇影视久久| 90打野战视频偷拍视频| 在线观看www视频免费| 亚洲avbb在线观看| 国产精品自产拍在线观看55亚洲 | 成年人免费黄色播放视频| 一级a爱视频在线免费观看| 国产区一区二久久| 19禁男女啪啪无遮挡网站| 女同久久另类99精品国产91| 亚洲第一欧美日韩一区二区三区 | 高潮久久久久久久久久久不卡| 欧美精品一区二区大全| 国产在线一区二区三区精| 国产成人av激情在线播放| 亚洲黑人精品在线| 国产亚洲精品第一综合不卡| 怎么达到女性高潮| 老鸭窝网址在线观看| 午夜福利视频在线观看免费| 国产成人免费无遮挡视频|