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

    基于GRACE 的黃河流域陸地水儲(chǔ)量時(shí)空變化研究

    2022-05-08 05:51:12張成業(yè)宋子恒郭添玉
    煤田地質(zhì)與勘探 2022年4期
    關(guān)鍵詞:黃河流域儲(chǔ)量振幅

    楊 飛,張成業(yè),李 軍,宋子恒,郭添玉

    (1.中國(guó)礦業(yè)大學(xué)(北京) 煤炭資源與安全開采國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100083;2.中國(guó)礦業(yè)大學(xué)(北京) 地球科學(xué)與測(cè)繪工程學(xué)院,北京 100083;3.煤炭開采水資源保護(hù)與利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102211)

    黃河流域是我國(guó)重要的經(jīng)濟(jì)地帶,構(gòu)成了我國(guó)重要的生態(tài)屏障。推進(jìn)黃河流域水資源節(jié)約集約利用是黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展的主要目標(biāo)任務(wù)之一。彭蘇萍等[1]指出,水資源是支撐黃河流域煤炭開發(fā)、區(qū)域經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境協(xié)調(diào)的基礎(chǔ)資源,而水資源短缺也是黃河流域中上游生態(tài)環(huán)境治理的瓶頸。因此,了解并掌握黃河流域水資源及其變化顯得尤為重要,是當(dāng)前黃河流域及其礦區(qū)生態(tài)環(huán)境研究的熱點(diǎn)方向。傳統(tǒng)監(jiān)測(cè)手段大多利用地面水資源動(dòng)態(tài)監(jiān)測(cè)網(wǎng)點(diǎn)及相應(yīng)的監(jiān)測(cè)系統(tǒng)來(lái)實(shí)現(xiàn),會(huì)受到監(jiān)測(cè)站點(diǎn)分布和數(shù)目的影響。

    GRACE(Gravity Recovery and Climate Experiment)重力衛(wèi)星計(jì)劃由美國(guó)國(guó)家宇航局和德國(guó)空間飛行中心聯(lián)合實(shí)施,衛(wèi)星于2002 年發(fā)射并于2004 年公開其重力數(shù)據(jù)[2]。GRACE 衛(wèi)星采用極地近圓軌道低衛(wèi)星跟蹤技術(shù),通過獲取星間距離的變化及其殘差來(lái)恢復(fù)時(shí)變重力場(chǎng)信息[3-4]。在消除天體引潮力對(duì)海洋、大氣和相關(guān)地球動(dòng)力過程引起的時(shí)變重力場(chǎng)影響后,GRACE 時(shí)變重力場(chǎng)反映的主要是兩極冰蓋、山岳冰川、陸地水儲(chǔ)量以及海平面變化信息[5-6]。因此,利用GRACE 時(shí)變重力場(chǎng)可以反演陸地水儲(chǔ)量變化信息,成為當(dāng)前進(jìn)行中長(zhǎng)尺度陸地水資源時(shí)空變化分析的重要手段[7-9]。J.Wahr 等[10-11]推導(dǎo)了GRACE 數(shù)據(jù)計(jì)算陸地水儲(chǔ)量變化的基礎(chǔ)理論公式,首次利用GRACE衛(wèi)星11 個(gè)月的數(shù)據(jù)反演密西西比河流域、亞馬孫河流域和孟加拉灣地區(qū)的儲(chǔ)水量變化,證明GRACE 具有在月尺度上探測(cè)到這些流域儲(chǔ)水量變化的潛力。隨后,眾多學(xué)者提出多種基于GRACE 的陸地水儲(chǔ)量反演的改進(jìn)方法[12-17],并成功應(yīng)用于世界范圍內(nèi)不同區(qū)域陸地水儲(chǔ)量的變化研究中[18-21]。

    當(dāng)前,利用GRACE 衛(wèi)星數(shù)據(jù)計(jì)算陸地水儲(chǔ)量變化的方法主要包括球諧系數(shù)方法和Mascons 方法[22-23]。鄒賢才等[24]詳細(xì)對(duì)比了上述2 種方法并發(fā)現(xiàn)Mascons 方法存在巨大優(yōu)勢(shì);其他學(xué)者的研究結(jié)果同樣證明Mascons 方法在一些流域的表現(xiàn)要優(yōu)于球諧系數(shù)方法[25-26]。2015 年以來(lái),美國(guó)得克薩斯州州立大學(xué)空間研究中心(CSR)和美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(JPL)發(fā)布用Mascons 方法計(jì)算陸地水儲(chǔ)量的產(chǎn)品,分別稱為CSRM 和JPL-M。2 種產(chǎn)品在不同區(qū)域陸地水儲(chǔ)量變化的比較和研究中表明,CSR-M 產(chǎn)品可以更好地反映整個(gè)區(qū)域地表陸地水遷移的實(shí)際情況[27-28]。

    基于此,利用CSR 提供的RL06 Mascons 水儲(chǔ)量變化產(chǎn)品,開展黃河流域上、中、下游近20 a 陸地水儲(chǔ)量變化的研究。通過緯圈長(zhǎng)度加權(quán)平均計(jì)算每個(gè)子流域的水儲(chǔ)量變化均值,采用時(shí)間序列分解方法得到黃河流域水儲(chǔ)量變化的趨勢(shì)、年周期等整體規(guī)律,利用空間分析比較不同區(qū)域水儲(chǔ)量變化的特征,以期為更好地了解和掌握黃河流域水資源及其變化規(guī)律提供依據(jù),為流域礦區(qū)的生態(tài)保護(hù)與可持續(xù)發(fā)展提供參考。

    1 黃河流域概況

    黃河發(fā)源于青海省巴顏喀拉山脈,橫跨青藏高原、內(nèi)蒙古高原、黃土高原和黃淮海平原,先后流經(jīng)青海、四川、甘肅、寧夏、內(nèi)蒙古、陜西、山西、河南和山東9 個(gè)省/自治區(qū)。其東西、南北跨度分別長(zhǎng)約1 900 km和1 100 km,經(jīng)緯度介于96°E-119°E 和32°N-42°N,流域總面積79.5 萬(wàn)km2。受大氣環(huán)流和季風(fēng)環(huán)流的影響作用比較復(fù)雜,黃河流域內(nèi)不同地區(qū)的氣候差異顯著,主要屬于南溫帶、中溫帶和高原氣候區(qū)。根據(jù)黃河水利委員會(huì)的劃分方案[29],黃河上游指的是從源頭到內(nèi)蒙古河口鎮(zhèn),河長(zhǎng)3 472 km,流域面積38.5 萬(wàn)km2;黃河中游指的是從內(nèi)蒙古河口鎮(zhèn)到河南鄭州桃花峪,河長(zhǎng)1 206 km,流域面積34.3 萬(wàn)km2;黃河下游指的是從鄭州桃花峪到渤海,河長(zhǎng)786 km,流域面積2.3 萬(wàn)km2。黃河流域中上游分布有晉北、晉東、晉中、黃隴、陜北、神東和寧東7 個(gè)國(guó)家大型煤炭基地,是我國(guó)目前主要煤炭生產(chǎn)區(qū)和調(diào)出區(qū);而黃河下流分布有河南、魯西等煤炭基地。黃河流域及其分界線如圖1 所示。

    圖1 黃流流域分布Fig.1 Map of the Yellow River Basin

    2 數(shù)據(jù)及分析方法

    本文采用美國(guó)德克薩斯大學(xué)空間研究中心CSR提供的2002 年4 月到2017 年6 月的最新版本RL06 Mascons 產(chǎn)品。較RL05 版本,RL06 使用了改進(jìn)的背景模型、新的地球系統(tǒng)參考框架和新的數(shù)據(jù)處理策略,顯著提升了數(shù)據(jù)質(zhì)量。CSR RL06 Mascons 產(chǎn)品是對(duì)GRACE 衛(wèi)星的高精度星間觀測(cè)數(shù)據(jù)Level1B 進(jìn)行直接解算,C20 項(xiàng)采用衛(wèi)星激光測(cè)月(SLR)數(shù)據(jù)的C20項(xiàng)替換,同時(shí)考慮地心一階項(xiàng)和冰后回彈改正,采用Tikhonov 正則化進(jìn)行方程組的求解。最終,CSR RL06 Mascons 數(shù)據(jù)采用NetCDF 格式,以0.25°×0.25°分辨率格網(wǎng)的形式發(fā)布。

    為分析2002-2017 年黃河流域上中下游各個(gè)區(qū)域的水儲(chǔ)量變化,需要計(jì)算每個(gè)區(qū)域的平均水儲(chǔ)量變化值。采用對(duì)每個(gè)格網(wǎng)點(diǎn)的水儲(chǔ)量變化數(shù)據(jù)直接取平均值的傳統(tǒng)方法,相當(dāng)于認(rèn)為每個(gè)格網(wǎng)點(diǎn)等權(quán)。然而,CSR RL06 Mascons 產(chǎn)品為0.25°×0.25°格網(wǎng)點(diǎn),數(shù)據(jù)在地球不同緯度之間代表的實(shí)際面積差異較大,需要考慮地球曲面的影響。本文采用緯圈長(zhǎng)度加權(quán)方法,根據(jù)地理緯度為每個(gè)格網(wǎng)點(diǎn)賦予權(quán)重來(lái)獲取不同區(qū)域的水儲(chǔ)量變化均值。計(jì)算公式如下:

    式中:Wr為某個(gè)區(qū)域的水儲(chǔ)量變化均值;Δh為某一個(gè)格網(wǎng)點(diǎn)的水儲(chǔ)量變化值;λi和φi分別為格網(wǎng)點(diǎn)經(jīng)度和緯度;n為該區(qū)域總格網(wǎng)點(diǎn)數(shù)量。

    本文采用時(shí)間序列分解方法來(lái)分析黃河流域上中下游不同區(qū)域水儲(chǔ)量變化的時(shí)間序列特征。將各個(gè)區(qū)域水儲(chǔ)量變化均值時(shí)間序列W(t)分解為趨勢(shì)項(xiàng)、年周期項(xiàng)和半年周期項(xiàng)。利用最小二乘法進(jìn)行參數(shù)擬合,擬合形式如下:

    式中:t為時(shí)間;a為常數(shù)項(xiàng);b為趨勢(shì)項(xiàng);c1、c2和d1、d2分別為年周期變化和半年周期變化。

    3 結(jié)果分析

    通過緯圈長(zhǎng)度加權(quán)方法得到每個(gè)月黃河流域上中下游的水儲(chǔ)量變化均值,采用等效水高表示[3-4,6-7],其時(shí)間序列及線性擬合趨勢(shì)如圖2 左側(cè)所示,可以看到,2002-2017 年,黃河上游水儲(chǔ)量變化均值的波動(dòng)范圍較小,最大和最小值分別為32.8、-56.9 mm,線性擬合趨勢(shì)項(xiàng)為-2.3 mm/a,水儲(chǔ)量變化遞減的趨勢(shì)不明顯。而黃河流域中游和下游的水儲(chǔ)量變化均值表現(xiàn)出明顯的遞減趨勢(shì),其趨勢(shì)項(xiàng)分別達(dá)到-10.3、-17.5 mm/a。在黃河流域中游,水儲(chǔ)量變化均值的最大和最小值分別為108.9、-168.6 mm;在黃河流域下游,它們分別達(dá)到了143.8、-301.9 mm。上述結(jié)果表明近20 a 黃河流域3 個(gè)流域段水儲(chǔ)量變化的差異性,上游水儲(chǔ)量變化穩(wěn)定,月均值差異較小,趨勢(shì)性弱;中游的水儲(chǔ)量變化呈遞減趨勢(shì);下游的水儲(chǔ)量變化范圍最大且遞減趨勢(shì)尤為明顯。

    去趨勢(shì)后,黃河流域上中下游水儲(chǔ)量變化均值的時(shí)間序列如圖2 右側(cè)所示,3 個(gè)流段都表現(xiàn)出周年性及季節(jié)性變化。值得注意的是,黃河中下游水儲(chǔ)量變化在2003 年秋季有顯著的提升。這是因?yàn)樵?003 年上半年,黃河流域來(lái)水量較常年減少40%~70%,為實(shí)測(cè)記錄以來(lái)最小值;進(jìn)入汛期以后,黃河流域大部分地區(qū)降水量增加,較常年增多10%~70%,洪水時(shí)水庫(kù)得以蓄水;所以在2003 年秋季,黃河流域水儲(chǔ)量變化迅速?gòu)呢?fù)值變化到正值。

    圖2 黃河流域上中下游水儲(chǔ)量變化均值的時(shí)間序列及趨勢(shì)和去趨勢(shì)后的時(shí)間序列Fig.2 Time series and trend of the mean water storage anomalies(left),and the time series after detrending(right) in the upper,middle and lower Yellow River Basin

    為進(jìn)一步分析黃河流域不同流段水儲(chǔ)量變化的月度規(guī)律以及數(shù)據(jù)的變化幅度,本文統(tǒng)計(jì)了2002 年4 月到2017 年6 月3 段流域各個(gè)月份的水儲(chǔ)量變化均值并繪制箱形圖,如圖3 所示??梢园l(fā)現(xiàn),黃河上中下游的水儲(chǔ)量變化月度范圍有顯著的區(qū)域差異,上游整體處于-40~40 mm,而中游和下游分別處于-100~100 mm和-200~200 mm。黃河上中下游每個(gè)月的數(shù)據(jù)都在一個(gè)相對(duì)穩(wěn)定的范圍內(nèi)波動(dòng),春冬季節(jié)的波動(dòng)明顯大于夏秋季節(jié),這與該流域“冬干春旱,夏秋多雨”的氣候特點(diǎn)密切相關(guān)。另外,有少數(shù)波動(dòng)表現(xiàn)為異常值,這些異常值集中于夏季,受到夏季降水和高溫蒸發(fā)等因素影響。3 個(gè)流域段多年平均水儲(chǔ)量變化在9 月份達(dá)到最大值,且在整個(gè)夏季和秋季維持較大值。

    圖3 黃河流域上中下游水儲(chǔ)量變化月數(shù)據(jù)箱形圖Fig.3 The boxplots of monthly mean data of the terrestrial water storage anomalies in the upper,middle and lower Yellow River Basin

    通過對(duì)黃河流域上中下游水儲(chǔ)量變化均值的時(shí)間序列進(jìn)行分解,得到其趨勢(shì)、年周期和半年周期等信息,見表1??梢钥吹?,時(shí)間序列分解得到的趨勢(shì)項(xiàng)與線性擬合的趨勢(shì)差別很小。黃河上游的趨勢(shì)項(xiàng)只有-2.17 mm/a,年振幅和半年振幅都較小,說(shuō)明該區(qū)域水儲(chǔ)量變化在近20 a 一直處于比較穩(wěn)定的狀態(tài)。相較于黃河中游,黃河下游的趨勢(shì)項(xiàng)、年振幅及半年振幅的增幅分別達(dá)到73%、38%和98%。上述結(jié)果表明,黃河流域上中下游水儲(chǔ)量變化存在不同的趨勢(shì)項(xiàng)和周期性。確定系數(shù)的范圍為0~1,表征最小二乘法參數(shù)擬合的可靠程度的統(tǒng)計(jì)指標(biāo)。上中下游的確定系數(shù)均大于0.5,表明時(shí)間序列分解得到較好的模型擬合參數(shù)。其中,中游確定系數(shù)為0.70,參數(shù)擬合結(jié)果優(yōu)于上游和下游。

    表1 黃河流域上中下游水儲(chǔ)量變化均值的時(shí)間序列分解統(tǒng)計(jì)Table 1 Statistics of the water storage anomalies in the upper,middle and lower Yellow River Basin through time series decomposition

    為表現(xiàn)黃河流域水儲(chǔ)量變化的空間分布情況,本文對(duì)該區(qū)域0.25°×0.25°格網(wǎng)點(diǎn)數(shù)據(jù)進(jìn)行時(shí)間序列分解,獲取每個(gè)格網(wǎng)點(diǎn)2002 年6 月到2017 年6 月的時(shí)間序列擬合參數(shù)?;诖?,開展黃河流域水儲(chǔ)量變化空間分析。圖4 表示水儲(chǔ)量變化趨勢(shì)(圖4a)和確定系數(shù)(圖4b)的空間分布??梢钥吹剑S河流域水儲(chǔ)量變化為增加趨勢(shì)的地區(qū)主要集中在黃河上游源頭附近區(qū)域,其年增加量在5 mm 左右。黃河下游及中游靠近黃河的區(qū)域表現(xiàn)出明顯的遞減趨勢(shì),年減少量基本達(dá)到15 mm 以上。其他一些區(qū)域趨勢(shì)項(xiàng)大多處于-10 mm/a 以下。黃河流域水儲(chǔ)量變化存在隨經(jīng)度變化的規(guī)律,由西向東表現(xiàn)出越來(lái)越明顯的遞減趨勢(shì)。從確定系數(shù)的空間分布可以看到黃河流域擬合結(jié)果的南北差異,流域北側(cè)擬合可靠度整體優(yōu)于流域南側(cè)。整個(gè)流域確定系數(shù)大于0.5 的格網(wǎng)點(diǎn)占64%,大于0.3的格網(wǎng)點(diǎn)達(dá)到85%。

    圖4 黃河流域水儲(chǔ)量參數(shù)的空間分布Fig.4 Spatial distribution of water storage factors in the Yellow River Basin

    圖5 表示黃河流域水儲(chǔ)量變化的年周期振幅和半年周期振幅的空間分布??梢钥吹?,整個(gè)流域年周期振幅都明顯大于其半年周期振幅。黃河上游源頭附近區(qū)域、中游桃花峪及黃河下游區(qū)域年周期振幅較大,達(dá)到30 mm 及以上。在上述區(qū)域,高山融雪、降水量等季節(jié)性差異是導(dǎo)致其年周期振幅相對(duì)較大的原因。黃河流域中上游其他區(qū)域的年周期振幅較小,處于10~20 mm。黃河流域水儲(chǔ)量變化的半年周期振幅表現(xiàn)出由西北向東南增大的趨勢(shì),在黃河上游小于10 mm,黃河中游處于10~20 mm,黃河下游達(dá)到了30 mm。這與不同區(qū)域所處的氣候環(huán)境等因素有直接關(guān)系。

    圖5 黃河流域水儲(chǔ)量不同變化周期的空間分布Fig.5 Spatial distribution of water storage anomalies in different period in the Yellow River Basin

    4 結(jié) 論

    a.黃河上游水儲(chǔ)量變化的趨勢(shì)性不明顯,只有-2.3 mm/a;中下游表現(xiàn)出顯著的遞減趨勢(shì),分別達(dá)到-10.3 mm/a 和-17.5 mm/a。GRACE 水儲(chǔ)量變化數(shù)據(jù)很好地反映出2003 年黃河流域來(lái)水量、降水量的特殊性。

    b.GRACE 水儲(chǔ)量變化數(shù)據(jù)表明黃河流域水儲(chǔ)量變化在9 月達(dá)到最大值,而春冬季節(jié)波動(dòng)較大,與該流域“冬干春旱,夏秋多雨”的氣候特點(diǎn)密切相關(guān)。水儲(chǔ)量變化的異常值集中于夏季,與夏季降水和高溫蒸發(fā)等因素相關(guān)。

    c.黃河上游源頭附近區(qū)域的水儲(chǔ)量變化呈增長(zhǎng)趨勢(shì),增長(zhǎng)率為5 mm/a 左右。黃河流域水儲(chǔ)量變化存在隨經(jīng)度變化的趨勢(shì),由西向東表現(xiàn)出越來(lái)越明顯的遞減趨勢(shì),最大遞減率超過20 mm/a。

    d.黃河流域年周期振幅明顯大于其半年周期振幅。受高山融雪、降水量等季節(jié)性差異的影響,年周期振幅在上游源頭附近區(qū)域、中游桃花峪及下游區(qū)域達(dá)到較大值。半年周期振幅表現(xiàn)出由西北向東南增大的趨勢(shì),可能與區(qū)域的氣候環(huán)境等因素相關(guān)。

    猜你喜歡
    黃河流域儲(chǔ)量振幅
    生態(tài)環(huán)境部啟動(dòng)新一年度黃河流域“清廢行動(dòng)”
    《礦產(chǎn)資源儲(chǔ)量技術(shù)標(biāo)準(zhǔn)》修訂對(duì)資源儲(chǔ)量報(bào)告編寫的影響
    基于三維軟件資源儲(chǔ)量估算對(duì)比研究
    黃河流域燦爛的齊家文化——陶器
    收藏界(2019年2期)2019-10-12 08:26:10
    增強(qiáng)大局意識(shí) 提升黃河流域生態(tài)保護(hù)發(fā)展水平
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    三十六計(jì)之順手牽羊
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    欧美激情久久久久久爽电影 | 色综合站精品国产| 久久久久亚洲av毛片大全| av免费在线观看网站| 亚洲情色 制服丝袜| 两人在一起打扑克的视频| 久久久久九九精品影院| 成在线人永久免费视频| 精品乱码久久久久久99久播| 他把我摸到了高潮在线观看| 国产一区二区三区视频了| 久久 成人 亚洲| 久久久久久亚洲精品国产蜜桃av| 一级作爱视频免费观看| 视频区图区小说| 久久人妻熟女aⅴ| 午夜福利在线免费观看网站| 亚洲男人天堂网一区| 亚洲片人在线观看| 乱人伦中国视频| 97碰自拍视频| 亚洲精品一卡2卡三卡4卡5卡| 久久久久国产精品人妻aⅴ院| 国产熟女xx| 久久影院123| 国产精品亚洲av一区麻豆| 久久精品成人免费网站| 亚洲av第一区精品v没综合| 久久精品国产清高在天天线| 国产亚洲精品第一综合不卡| 午夜视频精品福利| 男人舔女人的私密视频| 伊人久久大香线蕉亚洲五| 黑人巨大精品欧美一区二区mp4| 国产99久久九九免费精品| svipshipincom国产片| 热99re8久久精品国产| 成年人黄色毛片网站| 欧美日韩乱码在线| 99精国产麻豆久久婷婷| 在线观看一区二区三区激情| 欧美av亚洲av综合av国产av| 欧美成人午夜精品| 不卡av一区二区三区| 一级毛片女人18水好多| 夫妻午夜视频| 亚洲人成电影免费在线| 一级作爱视频免费观看| 夜夜爽天天搞| 国产亚洲欧美在线一区二区| 久久欧美精品欧美久久欧美| 三上悠亚av全集在线观看| 热re99久久精品国产66热6| 又黄又爽又免费观看的视频| 亚洲欧美精品综合一区二区三区| 亚洲国产中文字幕在线视频| 在线播放国产精品三级| 成熟少妇高潮喷水视频| 少妇的丰满在线观看| 亚洲成国产人片在线观看| 精品熟女少妇八av免费久了| 成人手机av| 中国美女看黄片| 亚洲性夜色夜夜综合| 成人影院久久| 十八禁网站免费在线| 99国产精品免费福利视频| 在线观看免费视频网站a站| 999久久久国产精品视频| 大型黄色视频在线免费观看| 岛国在线观看网站| 高清欧美精品videossex| 在线观看舔阴道视频| 黄色视频,在线免费观看| 日韩欧美免费精品| 丁香欧美五月| 纯流量卡能插随身wifi吗| 少妇粗大呻吟视频| 丝袜在线中文字幕| 丝袜在线中文字幕| 9191精品国产免费久久| 国产成+人综合+亚洲专区| 日韩国内少妇激情av| 咕卡用的链子| 黄色成人免费大全| 中文字幕精品免费在线观看视频| 在线观看66精品国产| 日韩精品中文字幕看吧| 手机成人av网站| 欧美精品啪啪一区二区三区| 国产熟女xx| 在线观看免费高清a一片| 国产在线精品亚洲第一网站| 制服人妻中文乱码| 亚洲在线自拍视频| 成人影院久久| 真人做人爱边吃奶动态| 757午夜福利合集在线观看| 国产深夜福利视频在线观看| 日韩欧美免费精品| 免费不卡黄色视频| 制服人妻中文乱码| 人人澡人人妻人| 亚洲 欧美一区二区三区| 成人影院久久| 91老司机精品| 久久人人爽av亚洲精品天堂| xxx96com| 亚洲色图综合在线观看| 看片在线看免费视频| 黄片播放在线免费| 脱女人内裤的视频| 看片在线看免费视频| 一边摸一边抽搐一进一出视频| 又黄又爽又免费观看的视频| 黄色 视频免费看| 多毛熟女@视频| 欧美大码av| 亚洲精品av麻豆狂野| 国产成人精品久久二区二区免费| 亚洲国产欧美网| 精品卡一卡二卡四卡免费| 91九色精品人成在线观看| 日韩国内少妇激情av| 欧美日韩黄片免| 免费在线观看影片大全网站| www.自偷自拍.com| 色在线成人网| 亚洲第一av免费看| 亚洲视频免费观看视频| 亚洲精品久久成人aⅴ小说| 欧美老熟妇乱子伦牲交| 国产精品综合久久久久久久免费 | 99香蕉大伊视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品美女特级片免费视频播放器 | 中国美女看黄片| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av| cao死你这个sao货| 91国产中文字幕| 操出白浆在线播放| 日本撒尿小便嘘嘘汇集6| 制服人妻中文乱码| 国产1区2区3区精品| 国产激情欧美一区二区| 99久久99久久久精品蜜桃| 91精品国产国语对白视频| 久久精品亚洲精品国产色婷小说| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕最新亚洲高清| 一级片'在线观看视频| 国产精品久久电影中文字幕| 午夜成年电影在线免费观看| 国产在线精品亚洲第一网站| 日韩大尺度精品在线看网址 | 91大片在线观看| 久久 成人 亚洲| 日韩免费高清中文字幕av| 中文字幕精品免费在线观看视频| 国产亚洲精品一区二区www| 老司机深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 国产免费男女视频| 亚洲av五月六月丁香网| 亚洲三区欧美一区| 女人被躁到高潮嗷嗷叫费观| 免费少妇av软件| 黄色女人牲交| 一个人免费在线观看的高清视频| 在线观看午夜福利视频| e午夜精品久久久久久久| 黑人猛操日本美女一级片| 人人澡人人妻人| 国产精品一区二区免费欧美| 午夜91福利影院| 在线观看免费高清a一片| 一边摸一边抽搐一进一出视频| 欧美在线一区亚洲| 欧美人与性动交α欧美精品济南到| 国产97色在线日韩免费| x7x7x7水蜜桃| 国产aⅴ精品一区二区三区波| 日日爽夜夜爽网站| 9热在线视频观看99| 天堂俺去俺来也www色官网| 不卡一级毛片| 妹子高潮喷水视频| 精品人妻1区二区| 黄片小视频在线播放| 女性生殖器流出的白浆| 亚洲国产欧美一区二区综合| a在线观看视频网站| 久久这里只有精品19| 国产成人啪精品午夜网站| 欧美久久黑人一区二区| 国产精品乱码一区二三区的特点 | 国产激情久久老熟女| 色尼玛亚洲综合影院| 麻豆一二三区av精品| 国产免费男女视频| 亚洲三区欧美一区| 国产精品影院久久| 国产蜜桃级精品一区二区三区| 久久久水蜜桃国产精品网| 美女高潮到喷水免费观看| 一夜夜www| 久久性视频一级片| 可以免费在线观看a视频的电影网站| 一级片'在线观看视频| 久久热在线av| 成人永久免费在线观看视频| 12—13女人毛片做爰片一| 欧美在线黄色| 亚洲熟妇中文字幕五十中出 | 久久热在线av| 亚洲一区中文字幕在线| 精品少妇一区二区三区视频日本电影| 国产av在哪里看| 国产精品国产av在线观看| 久久久国产精品麻豆| 欧美日本中文国产一区发布| 男女之事视频高清在线观看| www.www免费av| 两人在一起打扑克的视频| 国产精品亚洲av一区麻豆| 脱女人内裤的视频| 777久久人妻少妇嫩草av网站| 老熟妇仑乱视频hdxx| 在线观看一区二区三区激情| 人人妻人人爽人人添夜夜欢视频| 成人免费观看视频高清| 国产黄色免费在线视频| 日韩一卡2卡3卡4卡2021年| 最近最新中文字幕大全电影3 | 热99国产精品久久久久久7| 侵犯人妻中文字幕一二三四区| 神马国产精品三级电影在线观看 | 午夜福利免费观看在线| 色精品久久人妻99蜜桃| av片东京热男人的天堂| 亚洲精品久久午夜乱码| 欧美不卡视频在线免费观看 | 在线免费观看的www视频| 丝袜美腿诱惑在线| 日韩高清综合在线| 亚洲国产中文字幕在线视频| 首页视频小说图片口味搜索| 国产精品电影一区二区三区| 曰老女人黄片| 国产高清视频在线播放一区| 亚洲第一青青草原| 99国产综合亚洲精品| 最近最新免费中文字幕在线| 99在线视频只有这里精品首页| 夫妻午夜视频| 国产免费男女视频| 亚洲人成网站在线播放欧美日韩| 国产一区在线观看成人免费| 男女之事视频高清在线观看| 电影成人av| 亚洲一卡2卡3卡4卡5卡精品中文| 校园春色视频在线观看| 9191精品国产免费久久| 老熟妇乱子伦视频在线观看| 嫩草影视91久久| 欧美午夜高清在线| 久久中文字幕一级| 欧美乱妇无乱码| 法律面前人人平等表现在哪些方面| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片大片在线免费观看| 色综合欧美亚洲国产小说| 亚洲 欧美一区二区三区| 一边摸一边做爽爽视频免费| 亚洲少妇的诱惑av| 天堂√8在线中文| 欧美成狂野欧美在线观看| 欧美另类亚洲清纯唯美| 亚洲av成人一区二区三| 视频区欧美日本亚洲| 国内毛片毛片毛片毛片毛片| 性少妇av在线| 国产精品电影一区二区三区| 欧美精品亚洲一区二区| 亚洲国产精品sss在线观看 | 香蕉丝袜av| 国产极品粉嫩免费观看在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人一区二区三| 不卡一级毛片| 在线观看免费视频网站a站| 大陆偷拍与自拍| 久久天躁狠狠躁夜夜2o2o| 欧美中文日本在线观看视频| 黄色a级毛片大全视频| 欧美日韩精品网址| 乱人伦中国视频| 美女高潮到喷水免费观看| 久久午夜亚洲精品久久| 女人被狂操c到高潮| 日日夜夜操网爽| 亚洲精品中文字幕在线视频| 亚洲成人精品中文字幕电影 | 啪啪无遮挡十八禁网站| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 成人永久免费在线观看视频| 亚洲精品一区av在线观看| 女性生殖器流出的白浆| 中文字幕人妻熟女乱码| 国产高清videossex| 成年版毛片免费区| 国产不卡一卡二| 久久久久久大精品| 日本wwww免费看| 热99国产精品久久久久久7| 热re99久久精品国产66热6| 国产精品免费一区二区三区在线| 麻豆av在线久日| 精品久久久久久成人av| 女性被躁到高潮视频| 国产精品香港三级国产av潘金莲| 日韩人妻精品一区2区三区| 国产日韩一区二区三区精品不卡| 日韩国内少妇激情av| 美女午夜性视频免费| 高清毛片免费观看视频网站 | 看免费av毛片| 精品国产亚洲在线| 18禁国产床啪视频网站| 免费看十八禁软件| xxx96com| 老汉色∧v一级毛片| 美女高潮喷水抽搐中文字幕| 欧美乱妇无乱码| 久久人人爽av亚洲精品天堂| 精品电影一区二区在线| 90打野战视频偷拍视频| 18禁黄网站禁片午夜丰满| 午夜日韩欧美国产| 亚洲国产毛片av蜜桃av| 51午夜福利影视在线观看| 18禁黄网站禁片午夜丰满| 在线国产一区二区在线| 香蕉久久夜色| 午夜精品国产一区二区电影| 国产精品九九99| 亚洲一卡2卡3卡4卡5卡精品中文| 最近最新中文字幕大全电影3 | 日韩精品青青久久久久久| 久久香蕉国产精品| 国产成人欧美在线观看| 老司机亚洲免费影院| 国产高清激情床上av| 免费女性裸体啪啪无遮挡网站| 精品熟女少妇八av免费久了| 无限看片的www在线观看| 色在线成人网| 亚洲一码二码三码区别大吗| 国产成人啪精品午夜网站| www日本在线高清视频| 欧美黑人精品巨大| 亚洲一区二区三区欧美精品| 国产高清国产精品国产三级| 老司机靠b影院| 大型av网站在线播放| 久久久国产一区二区| 亚洲色图综合在线观看| 国产精品综合久久久久久久免费 | 亚洲欧美激情在线| 午夜福利,免费看| 国产精品美女特级片免费视频播放器 | 免费在线观看影片大全网站| 国产精品九九99| 18禁黄网站禁片午夜丰满| 久久精品国产综合久久久| 午夜a级毛片| 欧美成狂野欧美在线观看| 亚洲情色 制服丝袜| 久久久久亚洲av毛片大全| 久久久国产精品麻豆| 欧美中文日本在线观看视频| 国产精品99久久99久久久不卡| 亚洲,欧美精品.| bbb黄色大片| 欧美乱色亚洲激情| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 女人被狂操c到高潮| www日本在线高清视频| 日韩精品免费视频一区二区三区| 99在线视频只有这里精品首页| 黄色毛片三级朝国网站| 韩国精品一区二区三区| 午夜福利免费观看在线| 一边摸一边抽搐一进一小说| 国产成人av激情在线播放| 成人影院久久| 手机成人av网站| 黄片小视频在线播放| 首页视频小说图片口味搜索| 国产精品国产高清国产av| 91麻豆精品激情在线观看国产 | 看黄色毛片网站| 亚洲国产看品久久| 成年人黄色毛片网站| 91老司机精品| 日韩欧美国产一区二区入口| 久久人妻熟女aⅴ| 久久国产精品人妻蜜桃| 99久久人妻综合| 少妇裸体淫交视频免费看高清 | 国产精品九九99| 欧美日韩亚洲国产一区二区在线观看| 纯流量卡能插随身wifi吗| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 国产精品国产av在线观看| 三级毛片av免费| 97人妻天天添夜夜摸| 91成年电影在线观看| 日韩欧美在线二视频| 少妇粗大呻吟视频| 免费女性裸体啪啪无遮挡网站| 欧美大码av| 好男人电影高清在线观看| 男女下面插进去视频免费观看| 我的亚洲天堂| 丁香欧美五月| 亚洲人成电影观看| 又黄又粗又硬又大视频| 日韩欧美一区视频在线观看| 99热国产这里只有精品6| 国产免费男女视频| 99国产精品免费福利视频| 中文字幕色久视频| 不卡一级毛片| 最近最新免费中文字幕在线| 欧美乱色亚洲激情| 亚洲性夜色夜夜综合| 欧美最黄视频在线播放免费 | 欧美日韩中文字幕国产精品一区二区三区 | 又紧又爽又黄一区二区| 久久久久久大精品| 91老司机精品| 亚洲国产精品999在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品一区二区精品视频观看| 精品国产美女av久久久久小说| 一进一出抽搐动态| 在线观看免费视频日本深夜| 国产精品av久久久久免费| 国产日韩一区二区三区精品不卡| 搡老岳熟女国产| 69精品国产乱码久久久| 欧美日韩精品网址| 极品教师在线免费播放| 欧美久久黑人一区二区| 在线视频色国产色| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 欧美激情极品国产一区二区三区| 最新在线观看一区二区三区| 国产成人精品久久二区二区免费| 国产精品国产av在线观看| 最新在线观看一区二区三区| 久久久久久免费高清国产稀缺| 免费人成视频x8x8入口观看| 在线国产一区二区在线| 亚洲男人天堂网一区| 欧美中文日本在线观看视频| 国产99白浆流出| 一边摸一边抽搐一进一小说| 免费在线观看视频国产中文字幕亚洲| 国产亚洲欧美在线一区二区| 高潮久久久久久久久久久不卡| 久久精品国产清高在天天线| 久久国产亚洲av麻豆专区| 亚洲人成77777在线视频| 精品高清国产在线一区| 成年人免费黄色播放视频| 欧美久久黑人一区二区| 成人三级黄色视频| tocl精华| 一二三四在线观看免费中文在| 亚洲av第一区精品v没综合| 精品高清国产在线一区| 深夜精品福利| 成人18禁在线播放| 欧美av亚洲av综合av国产av| 久9热在线精品视频| 黄色成人免费大全| 国产有黄有色有爽视频| av网站在线播放免费| 乱人伦中国视频| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品久久久久久毛片| 亚洲av熟女| 美国免费a级毛片| 99热只有精品国产| 色综合欧美亚洲国产小说| 久久人人精品亚洲av| 成人手机av| 国产精品免费视频内射| 天堂中文最新版在线下载| 午夜福利影视在线免费观看| 亚洲精品美女久久av网站| 亚洲欧洲精品一区二区精品久久久| 一级a爱片免费观看的视频| 国产成人av教育| 97人妻天天添夜夜摸| 国产精品爽爽va在线观看网站 | 夜夜看夜夜爽夜夜摸 | 热99re8久久精品国产| 国产成人欧美| 桃红色精品国产亚洲av| 亚洲五月色婷婷综合| 老鸭窝网址在线观看| 国产成+人综合+亚洲专区| 性色av乱码一区二区三区2| 日韩免费高清中文字幕av| 成在线人永久免费视频| 成人精品一区二区免费| 国产一区二区在线av高清观看| 国产精品一区二区在线不卡| 巨乳人妻的诱惑在线观看| 老司机午夜十八禁免费视频| 神马国产精品三级电影在线观看 | 人人妻人人澡人人看| 亚洲成人久久性| 欧美日韩视频精品一区| 国产精品久久久av美女十八| 亚洲精华国产精华精| 好看av亚洲va欧美ⅴa在| 精品一区二区三区视频在线观看免费 | 国产精品免费一区二区三区在线| 精品熟女少妇八av免费久了| 精品久久久久久成人av| 亚洲自拍偷在线| 真人一进一出gif抽搐免费| 757午夜福利合集在线观看| 日本免费a在线| 精品久久久精品久久久| 精品免费久久久久久久清纯| 岛国在线观看网站| 日本a在线网址| 757午夜福利合集在线观看| 久久午夜亚洲精品久久| 男女午夜视频在线观看| 法律面前人人平等表现在哪些方面| 亚洲精品粉嫩美女一区| 亚洲一区中文字幕在线| 成人黄色视频免费在线看| 日韩成人在线观看一区二区三区| 搡老乐熟女国产| 黄色视频,在线免费观看| 午夜福利影视在线免费观看| 母亲3免费完整高清在线观看| 一边摸一边做爽爽视频免费| 日本撒尿小便嘘嘘汇集6| 精品国产美女av久久久久小说| 中文字幕高清在线视频| 国产亚洲精品久久久久5区| 久热爱精品视频在线9| 极品人妻少妇av视频| 另类亚洲欧美激情| 9191精品国产免费久久| 桃红色精品国产亚洲av| 亚洲精品一区av在线观看| 亚洲精品国产色婷婷电影| 露出奶头的视频| 国产精品98久久久久久宅男小说| 精品熟女少妇八av免费久了| 国产蜜桃级精品一区二区三区| 在线国产一区二区在线| 国产黄色免费在线视频| 久久人人精品亚洲av| 国产在线精品亚洲第一网站| 久久人妻熟女aⅴ| 久久婷婷成人综合色麻豆| 成人亚洲精品av一区二区 | 一本综合久久免费| 亚洲五月天丁香| 久久久久亚洲av毛片大全| av天堂久久9| 午夜亚洲福利在线播放| 久久婷婷成人综合色麻豆| 多毛熟女@视频| 午夜福利免费观看在线| 国产精品98久久久久久宅男小说| 免费不卡黄色视频| 亚洲激情在线av| 高清欧美精品videossex| 国产熟女午夜一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美网| 在线观看免费日韩欧美大片| 一级a爱片免费观看的视频| 一级毛片精品| 视频在线观看一区二区三区| 午夜影院日韩av| 一级毛片女人18水好多| 国产精品久久视频播放| 村上凉子中文字幕在线| 午夜免费鲁丝| 国产精品久久久久久人妻精品电影| 久久人人97超碰香蕉20202| 成人影院久久| 51午夜福利影视在线观看| 久久久久久人人人人人| 国产精品一区二区在线不卡| 精品免费久久久久久久清纯| 嫩草影视91久久| 亚洲情色 制服丝袜|