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

    基于GF衛(wèi)星解譯巴丹吉林沙漠湖泊水量變化

    2020-10-12 14:30曹樂聶振龍姜高磊
    人民黃河 2020年7期

    曹樂 聶振龍 姜高磊

    摘 要:遙感技術(shù)作為科學(xué)、快速、大面積調(diào)查監(jiān)測手段,在湖泊演化、動態(tài)變化研究中應(yīng)用廣泛。通過解譯國產(chǎn)GF-1、GF-2衛(wèi)星遙感影像信息,結(jié)合實地勘測,解譯了巴丹吉林沙漠110個湖泊地質(zhì)歷史時期與現(xiàn)代的湖面高程和邊界,總結(jié)了湖泊萎縮規(guī)律,估算了兩期水量變化。結(jié)果顯示:湖面高程平均降低9.76 m,面積共減少61.052 km2(占古湖總面積的75.49%),湖泊群水量共減少4.9億m3,說明湖泊萎縮程度高,沙漠干旱化趨勢明顯;研究區(qū)古、今地下水等水位線均表現(xiàn)出東南高、西北低的宏觀特征,反映了古、今沙漠區(qū)地下水相同的補給來源與徑流條件;基底凹陷區(qū)湖泊萎縮程度較低(如蘇木吉林湖),基底隆起區(qū)湖泊萎縮程度較高(如雅布賴山前湖泊),萎縮程度的差異性受區(qū)域構(gòu)造基底的控制與影響。

    關(guān)鍵詞:遙感解譯;古湖高程;湖泊萎縮;水量變化;巴丹吉林沙漠

    中圖分類號:TV213.3 ? 文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.07.009

    Abstract:Remote sensing technology is widely used in the study of lake evolution and dynamic changes as a scientific, rapid and large-scale investigation and monitoring method. This study obtained the elevations and boundaries of 110 lakes in ancient and modern in the Badain Jaran Desert, based on the GF-1 and GF-2 satellite remote sensing images and regional survey. The results show that the average lakes elevation is decreased by 9.76 m and the area of lakes reduced totally 61.052 km2, which accounts for 75.49% of the ancient lakes area. The total water volume of lakes has been reduced 490 million m3, showing high degree of shrinkage of lakes and obvious drought trend. The counter lines of groundwater have the characters of higher in southeast and lower in northwest generally whether in ancient or at present and those demonstrate the source of groundwater recharge and the groundwater runoff condition remains unchanged. The lakes shrinking are obviously in the basement uplift areas (e.g. Yabulai piedmont), while lightly in the basement depression areas (e.g. Lake Sumu Jaran). The differences of lakes shrinkage are controlled by the regional tectonic basement.

    Key words: remote sensing interpretation; ancient lakes elevation; lakes shrinking; status change; Badain Jaran Desert

    巴丹吉林沙漠位于阿拉善高原西部,面積4.9萬km2,沙漠中有高大沙山以及100多個湖泊,為世界罕見[1-3]。關(guān)于湖泊成因、演變、補給來源等問題一直是第四紀(jì)地質(zhì)與水文地質(zhì)界研究的熱點,但尚無定論[4-13]。遙感技術(shù)能夠全方位、大區(qū)域地進(jìn)行湖泊季節(jié)性、年際動態(tài)變化監(jiān)測,是湖泊研究的重要技術(shù)手段[14]。朱金峰等[15]利用ETM+遙感影像,提取了湖泊年內(nèi)季節(jié)變化信息,認(rèn)為湖泊面積和數(shù)量在年內(nèi)隨季節(jié)依次減少,到第二年春季又恢復(fù)到前一年春季狀態(tài)。張振瑜等[16-17]利用Landsat影像分析了1973—2010年湖泊的時空變化特征,認(rèn)為地下水補給源水量變化和補給方式的不同控制著湖泊面積的時空變化。Zhang等[18]利用Landsat遙感影像分析了1967—2010年50個湖泊的面積變化特征,認(rèn)為目前當(dāng)?shù)亟邓畬囱a給有限,地下水補給源來自于古湖殘余水體。Jiao等[19]利用衛(wèi)星數(shù)據(jù)與重力恢復(fù)試驗證明了沙漠湖泊水位與地下水水量在緩慢減少。

    根據(jù)湖水化學(xué)特征推測,該地區(qū)湖泊演化了800~7 000 a[1, 20],已有的高于現(xiàn)代湖水位的湖相沉積測年結(jié)果表明湖泊區(qū)在全新世中期處于泛湖期[21]。借助遙感技術(shù),對巴丹吉林沙漠湖泊的季節(jié)性變化和近幾十年的年際變化特征已有了初步的認(rèn)識[15-18],但對于地質(zhì)歷史時期湖泊面積、水量變化還未開展研究。筆者利用國產(chǎn)GF-1、GF-2衛(wèi)星PMS傳感器影像資料,結(jié)合Google Earth影像和實地調(diào)查對巴丹吉林沙漠110個湖泊進(jìn)行了演化解譯,確定了最大古湖的高程和邊界,以期為研究該區(qū)古地下水流場、湖泊萎縮特征等提供參考。

    1 研究區(qū)概況

    巴丹吉林沙漠位于雅布賴山與北大山以北、宗乃山以西、古日乃與弱水以東、拐子湖以南(見圖1)。

    沙漠區(qū)整體地勢東南高、西北低,海拔900~1 600 m,地形起伏較大,東南部高大沙山居多,相對高差一般為200~300 m[5, 22]。高大沙山間分布有100多個湖泊,呈串珠狀排列,多為咸水湖,面積一般小于1 km2[21]。沙漠區(qū)降水稀少,多年平均降水量為30~120 mm,由東南向西北逐漸減少[23-24],年平均水面蒸發(fā)量為1 450 mm[25]。沙漠區(qū)主要發(fā)育北東東向、北西向斷裂,控制著區(qū)內(nèi)基底結(jié)構(gòu)、隆起格架及中生代盆地的展布,使盆地形成凹凸相間的結(jié)構(gòu)特征,由南向北主要有雅布賴山隆起、蘇亥圖凹陷、宗乃山隆起、陶勒特凹陷、特羅西灘低隆起等構(gòu)造單元[26-27]。

    2 研究方法

    2.1 遙感解譯數(shù)據(jù)選取

    湖泊在演化過程中處于各穩(wěn)定期次時,在地表殘留的湖相沉積物、各類植被及其微地貌形態(tài)所共同構(gòu)成的規(guī)則閉合邊界只有在高分辨率遙感影像中才能準(zhǔn)確解譯。為保證解譯準(zhǔn)確性,將用于湖泊演化期次解譯的國產(chǎn)GF-1衛(wèi)星16 m分辨率、GF-1衛(wèi)星2 m分辨率(融合后)、GF-2衛(wèi)星1 m分辨率(融合后)數(shù)據(jù)以及Google Earth提取的0.6 m分辨率數(shù)據(jù)進(jìn)行比對。選取音德爾圖作為比對湖泊,以此確定湖泊演化期次提取所需最優(yōu)空間分辨率的區(qū)間范圍(見圖2)。結(jié)果顯示,演化期次邊界在亞米級影像上可以準(zhǔn)確識別,最優(yōu)空間分辨率為0.6 ~1.0 m (分辨率高于0.6 m對解譯結(jié)果準(zhǔn)確性的提升影響不大);而分辨率低于5 m的影像,難以識別湖泊演化期次邊界。故GF-1衛(wèi)星16 m分辨率影像不能滿足要求,GF-2衛(wèi)星及Google Earth所獲取的遙感數(shù)據(jù)較為適用。

    2.2 湖泊邊界的確定

    結(jié)合遙感影像與實地調(diào)查結(jié)果,將湖泊演化期次邊界分為古湖邊界、中間穩(wěn)定期次邊界、今湖邊界3類,各演化期次邊界及區(qū)域分布關(guān)系如圖3所示,古湖邊界以影像上可見面積最大的湖岸線為準(zhǔn),今湖邊界以2015年3月影像中湖泊水體面積的極大值為準(zhǔn)。

    各類演化區(qū)域邊界鑒定特征見圖4:①網(wǎng)格狀沙丘區(qū),由微型曲弧狀沙梁及下凹沙窩構(gòu)成,網(wǎng)格及蜂窩狀形態(tài)明顯,沙窩內(nèi)發(fā)育灌叢及草本植物,整體色調(diào)較淺,多為土黃色,與外側(cè)緩起伏沙地及沙山區(qū)別明顯,其邊緣指示古湖邊界;②草灌叢沙堆區(qū)、植被覆蓋區(qū)、鹽漬化區(qū)、潛水蒸發(fā)區(qū)影像上形態(tài)起伏明顯,表面平滑與粗糙相異,色彩深淺不同,兩區(qū)交界邊緣均指示某一湖泊穩(wěn)定期邊界;③季節(jié)性湖泊區(qū),影像表面光滑,水體多呈不規(guī)則形態(tài),水體間裸露少量土體,可見斑狀白色鹽堿,整體色調(diào)呈深藍(lán)色,水體季節(jié)性變化明顯,故其邊界年內(nèi)變化較大,指示今湖邊界;④常年性湖泊區(qū),影像表面光滑,水體形態(tài)規(guī)則完整,色調(diào)為深藍(lán)色,水體面積季節(jié)性變化小,邊界年內(nèi)變化小,指示今湖邊界。

    2.3 湖泊高程提取與驗證

    對于研究古湖及今湖高程特征而言,其相對高程比絕對高程更為重要,為此從不同來源的數(shù)據(jù)獲取固定點對的高程差,判斷相對高程的誤差特征。比對野外移動GPS獲取的高程數(shù)據(jù)與室內(nèi)提取的DEM高程數(shù)據(jù),采用RTK(Real-Time Kinematic)方法所采集的高精度高程值進(jìn)行驗證(控制點相對高程誤差控制在5 cm內(nèi))。表1顯示,以RTK數(shù)據(jù)為基準(zhǔn),DEM數(shù)據(jù)所得點對相對高程誤差小于2 m,無正負(fù)號偏差,在區(qū)域研究尺度上是可以接受的;而移動GPS的相對高程出現(xiàn)正負(fù)號偏差,且數(shù)值與DEM高程相差可達(dá)10 m以上。因此,本研究統(tǒng)一采用DEM高程數(shù)據(jù)。

    古湖/今湖邊界的高程提取方法是利用已有的古湖/今湖矢量邊界,將其等分為多條(如100條)矢量線段后,獲取各矢量線段中點,分別提取各中點對應(yīng)的DEM高程數(shù)據(jù)并求其平均值,最終獲得的邊界高程均值為其高程值。高程提取主要步驟如圖5所示。

    2.4 野外驗證

    選擇代表性湖泊進(jìn)行演化期次邊界野外驗證,驗證古湖67個,占湖泊總數(shù)的61%;通過野外典型標(biāo)志的確定(見圖6),證明了解譯的準(zhǔn)確性。以昂斯格湖與巴丹湖區(qū)驗證為例,昂斯格湖影像(圖6(b))中箭頭①②為古湖邊界,實地調(diào)查時發(fā)現(xiàn)①處的古湖邊界處于沙山半坡上,明顯可見古湖邊界近水平出露于沙山表面(圖6(a));②處指示湖泊東部古湖邊界,實地可見古湖沉積物出露地表,其高程與對岸古湖邊界一致(圖6(a));③處指示網(wǎng)格狀沙丘區(qū)(圖6(a))。巴丹湖區(qū)東北部出露的古湖沉積物與遙感劃定古湖邊界位置相符(圖6(c)),據(jù)野外調(diào)查,巴丹湖及附近湖泊可能曾經(jīng)為同一個大湖,在后期演化萎縮過程中,古湖區(qū)被沙體掩埋覆蓋,最終形成現(xiàn)在的多個湖泊。

    3 結(jié)果分析

    3.1 湖泊高程與面積變化

    圖7為所有古湖與今湖的高程與面積,可以看出,所有今湖湖面高程下降,平均降幅為9.76 m,面積萎縮。湖面高程下降0~10 m的有61個湖泊,下降10~20 m的有43個湖泊,下降20~25 m的有6個湖泊。從解譯的古湖面積來看,最大的為昂斯格,面積3.412 6 km2,最小的為葫蘆斯臺(東湖),面積僅0.007 2 km2;從解譯的今湖面積來看,最大的為東諾爾圖,面積1.417 5 km2,最小的為葫蘆斯臺(中湖),面積僅為0.000 4 km2。作為地表水體,地質(zhì)歷史時期最大古湖面積總計80.877 km2,今湖現(xiàn)存水體面積總計19.825 km2,減小61.052 km2,占古湖總面積的75.49%。表2顯示,110個湖泊的萎縮程度((1-今湖面積/古湖面積)×100%)均大于25%,大于90%的湖泊有48個,大于95%的有33個。不同湖泊的萎縮程度各不相同,萎縮程度最大的湖泊是達(dá)布蘇圖巴潤敖格欽,今湖水體面積僅為古湖面積的0.12%;萎縮程度最小的湖泊是扎拉特,今湖水體面積占古湖面積的74.96%。

    3.2 湖泊水量平衡

    圖8為110個湖泊的水量平衡計算結(jié)果(湖泊編號按湖泊面積從大到小排列),可以看出,所有湖泊體積減小,面積越大的古湖水量減少越多。計算得到110個湖泊水量總共減少4.9億m3,其中:蘇木吉林與巴潤蘇木吉林作為一個古湖萎縮分裂成的兩個湖泊,減少的水量最大,為0.32億m3;未分裂的單個湖泊中,伊和吉格德減小量最大,為0.27億m3。

    3.3 古、今湖泊等水位線特征

    地表水位的分布特征一定程度上受地下水系統(tǒng)影響,因此可利用湖面高程對地下水系統(tǒng)特征進(jìn)行研究。圖9顯示,利用古、今湖邊界平均高程數(shù)據(jù)得出的等水位線分布特征基本一致,二者均顯示出地下水位東南高、西北低的宏觀特征;東南部地下水水力梯度較大,西北部相對較小,說明地下水流向未有明顯變化。在圖9中選取3點P1、P2、P3,分段計算水力梯度,結(jié)果見表3,顯示今湖的地下水水力梯度均大于古湖的,P1—P2段的水力梯度較大,為0.707%;從東南部水力梯度較大區(qū)到西北部,水力梯度逐漸減小,P2—P3段今湖的平均水力梯度為0.070%,P1—P3段今湖的地下水平均水力梯度為0.219%。

    0.070%~0.707%的水力梯度變幅與張競等[13, 29]的研究結(jié)果0.08%~0.79%很相近,說明古湖與今湖地下水由東南向西北的主流向一致,反映地下水補給來源與補給模式未發(fā)生改變。地下水是歷史時期東南部雅布賴山區(qū)降水徑流及山前古河道滲漏補給而來,馬金珠等[30]、黃天明等[31]利用環(huán)境示蹤劑研究認(rèn)為是更新世晚期至全新世早期周邊山區(qū)降水徑流補給的古水。

    3.4 湖泊萎縮特征

    圖10顯示了湖泊萎縮程度區(qū)的分布特征,為湖泊萎縮程度差分模擬結(jié)果。區(qū)內(nèi)東南部、東北部、中北部及西部湖泊面積萎縮程度較大,而中部、中東部、南部及北部湖泊面積萎縮程度較小。結(jié)合區(qū)域地質(zhì)構(gòu)造特征(見圖1(a)),凹陷區(qū)基底較深,地下水匯集于此,水量豐富,當(dāng)氣候干旱時,有地下水補給的湖泊萎縮較慢;相反,隆起區(qū)湖泊對氣候干旱、水量減少反應(yīng)更為敏感。圖10中Ⅰ區(qū)、Ⅱ區(qū)、Ⅲ區(qū)分別位于雅布賴山隆起、因格井南凸起及宗乃山隆起帶,整體萎縮程度較大。Ⅳ區(qū)、Ⅴ區(qū)、Ⅵ區(qū)屬于蘇亥圖凹陷區(qū),其中:Ⅴ區(qū)為沙漠大湖分布區(qū),核心區(qū)為蘇木吉林湖,位于蘇亥圖凹陷區(qū)中部,萎縮程度最為緩慢;Ⅳ區(qū)、Ⅵ區(qū)萎縮程度略大于Ⅴ區(qū),表現(xiàn)出凹陷區(qū)中心至邊界基底深度的漸變(Ⅵ區(qū))與局部變化(Ⅳ區(qū)),造成湖泊萎縮程度有差異。Ⅶ區(qū)位于因格井凹陷南端,相對較深的基底保證了更厚的含水層,湖泊萎縮程度相對較小,地下水接受其東側(cè)的宗乃山山區(qū)降水徑流補給。

    4 結(jié) 論

    (1)巴丹吉林沙漠現(xiàn)代湖泊湖面高程平均低于古湖面9.76 m,湖泊總面積萎縮減小61.052 km2,占古湖總面積的75.49%,湖泊水量共減少4.9億m3,說明湖泊萎縮程度較高,沙漠干旱化趨勢明顯。

    (2)湖泊區(qū)古、今等水位線均表現(xiàn)出東南高、西北低的宏觀特征,古、今水力梯度變化不明顯,反映了長時間尺度地下水有相同的補給來源,可能是周邊山區(qū)特別是東南部雅布賴山區(qū)古降水徑流及山前古河道下滲補給。

    (3)湖泊萎縮的差異性受區(qū)域內(nèi)構(gòu)造基底的控制與影響,凹陷區(qū)萎縮程度相對較低,是大湖的主要分布區(qū)。

    參考文獻(xiàn):

    [1] YANG X P, WILLIAMS M A J. The Ion Chemistry of Lakes and Late Holocene Desiccation in the Badain Jaran Desert, Inner Mongolia, China[J]. Catena, 2003,51(1):45-60.

    [2] ANDREOTTI B, FOURRIERE A, OULD-KADDOU F, et al. Giant Aeolian Dune Size Determined by the Average Depth of the Atmospheric Boundary Layer[J]. Nature, 2009,457:1120-1123.

    [3] DONG Z B, QIAN G Q, LYU P, et al. Investigation of the Sand Sea with the Tallest Dunes on Earth: Chinas Badain Jaran Sand Sea[J]. Earth-Science Reviews, 2013,120:20-39.

    [4] CHEN J S, LI L, WANG J Y, et al. Groundwater Maintains Dune Landscape[J]. Nature,2004,432:459-460.

    [5] 陳建生, 趙霞, 盛雪芬, 等.巴丹吉林沙漠湖泊群與沙山形成機理研究[J].科學(xué)通報,2006,51(23):2789-2796.

    [6] YANG X P, SCUDERI L A. Hydrological and Climatic Changes in Deserts of China Since the Late Pleistocene[J]. Quaternary Research,2010,73(1):1-9.

    [7] GATES J B, EDMUNDS W M, DARLING W G, et al. Conceptual Model of Recharge to Southeastern Badain Jaran Desert Groundwater and Lakes from Environmental Tracers[J]. Applied Geochemistry, 2008,23(12):3519-3534.

    [8] GATES J B, EDMUNDS W M, MA J Z, et al. Estimating Groundwater Recharge in a Cold Desert Environment in Northern China Using Chloride[J]. Hydrogeology Journal, 2008,16(5):893-910.

    [9] MA J Z, EDMUNDS W M. Groundwater and Lake Evolution in the Badain Jaran Desert Ecosystem, Inner Mongolia[J]. Hydrogeology Journal,2006,14(7):1231-1243.

    [10] MA J Z, DING Z, GATES J B, et al. Chloride and the Environmental Isotopes as the Indicators of the Groundwater Recharge in the Gobi Desert, Northwest China[J]. Environmental Geology,2008,55(7):1407-1419.

    [11] WANG Z T, CHEN T Y, LIU S W, et al. Aeolian Origin of Inter Dune Lakes in the Badain Jaran Desert, China[J]. Arabian Journal of Geosciences,2016,9(3):190.

    [12] 丁宏偉,郭瑞,藍(lán)永超,等.再論巴丹吉林沙漠湖泊水的補給來源、補給模式與高大沙山的形成機理[J].冰川凍土,2015,37(3):783-792.

    [13] 張競,王旭升,胡曉農(nóng),等.巴丹吉林沙漠地下水流場的宏觀特征[J].中國沙漠,2015,35(3):774-782.

    [14] 王海波,馬明國.基于遙感的湖泊水域動態(tài)變化監(jiān)測研究進(jìn)展[J].遙感技術(shù)與應(yīng)用,2009,24(5):674-684.

    [15] 朱金峰,王乃昂,李卓侖,等.巴丹吉林沙漠湖泊季節(jié)變化的遙感監(jiān)測[J].湖泊科學(xué),2011,23(4):657-664.

    [16] 張振瑜,王乃昂,馬寧,等.近40 a巴丹吉林沙漠腹地湖泊面積變化及其影響因素[J].中國沙漠,2012,32(6):1743-1750.

    [17] 張振瑜,王乃昂,吳月,等.1973—2010年巴丹吉林沙漠腹地湖泊面積空間變化的遙感分析[J].湖泊科學(xué),2013,25(4):514-520.

    [18] ZHANG Z, DONG Z B, YAN C Z, et al. Change of Lake Area in the Southeastern Part of Chinas Badain Jaran Sand Sea and Its Implications for Recharge Sources[J]. Journal of Arid Land,2015,7(1):1-9.

    [19] JIAO J J, ZHANG X T, WANG X S. Satellite-Based Estimates of Groundwater Depletion in the Badain Jaran Desert, China[J]. Scientific Reports,2015,5(1):8960.

    [20] GONG Y P, WANG X S, HU B X, et al. Groundwater Contributions in Water-Salt Balances of the Lakes in the Badain Jaran Desert, China[J]. Journal of Arid Land,2016,8(5):694-706.

    [21] WANG N A, NING K, LI Z L, et al. Holocene High Lake-Levels and Pan-Lake Period on Badain Jaran Desert[J]. Science China:Earth Sciences,2016,59(8):1633-1641.

    [22] 朱金峰,王乃昂,陳紅寶,等. 基于遙感的巴丹吉林沙漠范圍與面積分析[J].地理科學(xué)進(jìn)展,2010,29(9):1087-1094.

    [23] 王乃昂,馬寧,陳紅寶,等.巴丹吉林沙漠腹地降水特征的初步分析[J].水科學(xué)進(jìn)展,2013,24(2):153-160.

    [24] 張克存,姚正毅,安志山,等.巴丹吉林沙漠及其毗鄰地區(qū)降水特征及風(fēng)沙環(huán)境分析[J].中國沙漠,2012,32(6):1507-1511.

    [25] 胡文峰,王乃昂,趙力強,等.巴丹吉林沙漠典型湖泊湖氣界面水-熱交換特征[J].地理科學(xué)進(jìn)展,2015,34(8):1061-1071.

    [26] 劉建利,申安斌,陳小龍.大地電磁測深方法在內(nèi)蒙古西部銀根-額濟納旗盆地石炭系—二疊系油氣地質(zhì)調(diào)查中的應(yīng)用[J].地質(zhì)通報,2011,30(6):993-1000.

    [27] 嚴(yán)云奎,袁炳強,楊高印,等.內(nèi)蒙古西部銀根-額濟納旗盆地重力場與斷裂構(gòu)造的特征[J].地質(zhì)通報,2011,30(12):1962-1968.

    [28] 張國慶,XIE Hongjie,姚檀棟,等.基于ICESat和Landsat的中國十大湖泊水量平衡估算[J].科學(xué)通報,2013,58(26):2664-2678.

    [29] 張競,王旭升,賈鳳超,等.對內(nèi)蒙古阿拉善西部地下水流向問題的新認(rèn)識[J].現(xiàn)代地質(zhì),2015,29(1):213-219.

    [30] 馬金珠,黃天明,丁貞玉,等.同位素指示的巴丹吉林沙漠南緣地下水補給來源[J].地球科學(xué)進(jìn)展,2007,22(9):922-930.

    [31] 黃天明,龐忠和.應(yīng)用環(huán)境示蹤劑探討巴丹吉林沙漠及古日乃綠洲地下水補給[J].現(xiàn)代地質(zhì),2007,21(4):624-631.

    【責(zé)任編輯 張華興】

    最近最新中文字幕大全免费视频| 国产精品一区二区三区四区久久| 可以在线观看毛片的网站| 国产亚洲欧美98| 日本 av在线| 国产在线观看jvid| 国产精品亚洲一级av第二区| 日韩欧美国产在线观看| 欧美3d第一页| 亚洲18禁久久av| 国产高清视频在线观看网站| 亚洲av成人一区二区三| 日韩欧美一区二区三区在线观看| 视频区欧美日本亚洲| 18禁黄网站禁片午夜丰满| 国内少妇人妻偷人精品xxx网站 | 亚洲精华国产精华精| 国产精品日韩av在线免费观看| 一二三四社区在线视频社区8| 成人国产综合亚洲| 搞女人的毛片| 亚洲午夜精品一区,二区,三区| 久久久国产成人精品二区| 97超级碰碰碰精品色视频在线观看| 热99re8久久精品国产| 欧美乱妇无乱码| 男人的好看免费观看在线视频 | 国产熟女xx| 欧美乱色亚洲激情| 搡老熟女国产l中国老女人| 一区二区三区国产精品乱码| 国产1区2区3区精品| 成人欧美大片| 免费一级毛片在线播放高清视频| 午夜精品一区二区三区免费看| 黄色丝袜av网址大全| 亚洲欧美精品综合久久99| 国产区一区二久久| 久久久精品国产亚洲av高清涩受| 国产成人aa在线观看| 国产成人一区二区三区免费视频网站| tocl精华| 亚洲五月婷婷丁香| 给我免费播放毛片高清在线观看| 成人午夜高清在线视频| 国产高清有码在线观看视频 | 精品熟女少妇八av免费久了| 制服诱惑二区| tocl精华| 久久久精品国产亚洲av高清涩受| 国产真实乱freesex| 一进一出抽搐动态| 国内精品久久久久久久电影| 亚洲欧美激情综合另类| 久久久精品国产亚洲av高清涩受| 女警被强在线播放| e午夜精品久久久久久久| 婷婷丁香在线五月| 精品一区二区三区四区五区乱码| 欧美日韩一级在线毛片| 免费在线观看黄色视频的| 国产麻豆成人av免费视频| 久9热在线精品视频| 午夜激情福利司机影院| 亚洲成a人片在线一区二区| 最新在线观看一区二区三区| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 高潮久久久久久久久久久不卡| 亚洲av中文字字幕乱码综合| 真人做人爱边吃奶动态| 日韩欧美在线乱码| 国产成人精品久久二区二区免费| 午夜激情福利司机影院| 久久久久九九精品影院| 国产欧美日韩一区二区精品| 在线观看日韩欧美| 97超级碰碰碰精品色视频在线观看| 91麻豆av在线| 成人18禁在线播放| 国产在线精品亚洲第一网站| 亚洲av电影不卡..在线观看| 日韩欧美在线二视频| www.熟女人妻精品国产| 一区二区三区国产精品乱码| 国产熟女xx| 日本熟妇午夜| 在线a可以看的网站| 大型黄色视频在线免费观看| 国产一区二区在线观看日韩 | 免费人成视频x8x8入口观看| 国产精品久久久久久人妻精品电影| 禁无遮挡网站| 国产熟女午夜一区二区三区| 精品久久蜜臀av无| 国产人伦9x9x在线观看| 日韩大尺度精品在线看网址| 国产成+人综合+亚洲专区| 欧美日韩精品网址| 两个人的视频大全免费| 久久人人精品亚洲av| 国产91精品成人一区二区三区| 19禁男女啪啪无遮挡网站| 99国产综合亚洲精品| 中亚洲国语对白在线视频| 桃红色精品国产亚洲av| 亚洲成a人片在线一区二区| 九色国产91popny在线| 久久精品国产综合久久久| 久热爱精品视频在线9| 两个人的视频大全免费| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| 亚洲av日韩精品久久久久久密| 亚洲国产精品成人综合色| 成年人黄色毛片网站| 99久久综合精品五月天人人| 淫妇啪啪啪对白视频| 欧美日韩国产亚洲二区| 亚洲 欧美 日韩 在线 免费| 国产精品久久久av美女十八| 欧美黄色片欧美黄色片| 久久精品综合一区二区三区| 黄色 视频免费看| 欧美 亚洲 国产 日韩一| 一级黄色大片毛片| 两个人的视频大全免费| 啦啦啦免费观看视频1| 婷婷六月久久综合丁香| 色综合站精品国产| 亚洲国产欧美一区二区综合| 欧美日韩福利视频一区二区| 国产v大片淫在线免费观看| 成人欧美大片| 小说图片视频综合网站| 国产精品久久久久久亚洲av鲁大| 久久精品国产清高在天天线| 国产黄色小视频在线观看| 国产精华一区二区三区| 在线看三级毛片| 中文在线观看免费www的网站 | 男女床上黄色一级片免费看| 日韩欧美免费精品| 色精品久久人妻99蜜桃| 国内精品久久久久久久电影| 欧美大码av| 欧美日韩国产亚洲二区| 在线观看www视频免费| 午夜福利18| 国产欧美日韩精品亚洲av| 搡老妇女老女人老熟妇| 麻豆成人午夜福利视频| 一卡2卡三卡四卡精品乱码亚洲| 脱女人内裤的视频| 国产精品1区2区在线观看.| 久久久久亚洲av毛片大全| 老汉色av国产亚洲站长工具| 欧美zozozo另类| 亚洲人成伊人成综合网2020| 91麻豆av在线| 黄色丝袜av网址大全| 亚洲午夜理论影院| 又紧又爽又黄一区二区| 久久久国产欧美日韩av| 一进一出好大好爽视频| 亚洲人成77777在线视频| 午夜福利视频1000在线观看| 日韩欧美在线二视频| 亚洲国产日韩欧美精品在线观看 | 国产精品美女特级片免费视频播放器 | 国产高清视频在线播放一区| 香蕉av资源在线| 男女之事视频高清在线观看| 国产成人一区二区三区免费视频网站| 精品久久久久久久久久久久久| 亚洲av片天天在线观看| 亚洲精品在线美女| 黄色女人牲交| 亚洲精品中文字幕一二三四区| 亚洲人成电影免费在线| 伦理电影免费视频| 国产伦一二天堂av在线观看| 老汉色∧v一级毛片| 一区福利在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本一二三区视频观看| 日韩国内少妇激情av| 一区二区三区高清视频在线| 亚洲成人国产一区在线观看| a级毛片a级免费在线| 在线观看免费日韩欧美大片| 国产单亲对白刺激| 啦啦啦韩国在线观看视频| 变态另类成人亚洲欧美熟女| 成人三级黄色视频| 国产精品日韩av在线免费观看| 日韩大码丰满熟妇| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看 | 亚洲五月婷婷丁香| 日韩欧美在线二视频| 女人爽到高潮嗷嗷叫在线视频| 国内少妇人妻偷人精品xxx网站 | 久久这里只有精品中国| 一个人免费在线观看电影 | 成年免费大片在线观看| 久久久久久久久免费视频了| 国产黄片美女视频| 美女免费视频网站| 久久久久久国产a免费观看| 一本久久中文字幕| 伦理电影免费视频| 亚洲av五月六月丁香网| 18禁黄网站禁片免费观看直播| 久久中文字幕一级| 女警被强在线播放| 日日爽夜夜爽网站| 麻豆久久精品国产亚洲av| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 午夜福利成人在线免费观看| 91九色精品人成在线观看| 熟女电影av网| 欧美黑人精品巨大| 亚洲精华国产精华精| 精品日产1卡2卡| 成人高潮视频无遮挡免费网站| 无限看片的www在线观看| 天天一区二区日本电影三级| 国产av在哪里看| 极品教师在线免费播放| 久久久精品国产亚洲av高清涩受| 国产精品一区二区免费欧美| 白带黄色成豆腐渣| 国产一区二区在线观看日韩 | 最近最新中文字幕大全电影3| www日本在线高清视频| 在线视频色国产色| 91av网站免费观看| 欧美日韩福利视频一区二区| 午夜成年电影在线免费观看| 国内精品一区二区在线观看| 少妇人妻一区二区三区视频| 久久中文字幕一级| 中出人妻视频一区二区| 香蕉久久夜色| av天堂在线播放| or卡值多少钱| 国产区一区二久久| 亚洲一区二区三区不卡视频| 欧美又色又爽又黄视频| 色综合站精品国产| 亚洲人成电影免费在线| 一级毛片女人18水好多| 99久久久亚洲精品蜜臀av| 麻豆国产97在线/欧美 | 成人欧美大片| 午夜视频精品福利| 欧美黄色片欧美黄色片| 午夜免费激情av| 亚洲在线自拍视频| 一区二区三区高清视频在线| 久久久久久国产a免费观看| 在线观看免费午夜福利视频| 久久 成人 亚洲| 精品电影一区二区在线| 久久久久久久午夜电影| 国内精品久久久久久久电影| 看黄色毛片网站| 欧美日本亚洲视频在线播放| 欧美日韩亚洲综合一区二区三区_| 国产精品美女特级片免费视频播放器 | 一边摸一边抽搐一进一小说| 精品久久蜜臀av无| 久久久久免费精品人妻一区二区| 婷婷精品国产亚洲av在线| 听说在线观看完整版免费高清| 两人在一起打扑克的视频| 人人妻人人看人人澡| 此物有八面人人有两片| 亚洲人成77777在线视频| 久久久久性生活片| 国产亚洲精品一区二区www| 午夜福利18| 成人午夜高清在线视频| 黄色丝袜av网址大全| 一二三四在线观看免费中文在| 1024香蕉在线观看| 制服诱惑二区| 亚洲 欧美一区二区三区| 一区二区三区高清视频在线| 亚洲av日韩精品久久久久久密| 国产精品美女特级片免费视频播放器 | a在线观看视频网站| 国产又黄又爽又无遮挡在线| 少妇粗大呻吟视频| 麻豆国产97在线/欧美 | 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看| 国产成人精品久久二区二区91| 婷婷丁香在线五月| 国产成+人综合+亚洲专区| 午夜免费成人在线视频| 久久中文字幕人妻熟女| 又爽又黄无遮挡网站| 婷婷亚洲欧美| 亚洲成人精品中文字幕电影| 亚洲精品中文字幕一二三四区| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 香蕉av资源在线| 高清毛片免费观看视频网站| 久久国产精品人妻蜜桃| 成人高潮视频无遮挡免费网站| 国产一区在线观看成人免费| 欧美黑人巨大hd| 男女那种视频在线观看| 午夜视频精品福利| 免费观看人在逋| 精品久久久久久成人av| 三级男女做爰猛烈吃奶摸视频| 亚洲成人国产一区在线观看| 久久久久国内视频| 黑人欧美特级aaaaaa片| 两个人的视频大全免费| 亚洲九九香蕉| 亚洲av电影不卡..在线观看| 久热爱精品视频在线9| 国产午夜精品久久久久久| 久久久久性生活片| 曰老女人黄片| 欧美黄色片欧美黄色片| 女警被强在线播放| 久久久久免费精品人妻一区二区| 日韩成人在线观看一区二区三区| 精品国产乱码久久久久久男人| 老司机福利观看| 国产黄片美女视频| 日韩欧美三级三区| 亚洲人与动物交配视频| 在线观看舔阴道视频| 99re在线观看精品视频| 桃红色精品国产亚洲av| 久久 成人 亚洲| 亚洲自拍偷在线| 亚洲精品久久国产高清桃花| 50天的宝宝边吃奶边哭怎么回事| 男人舔女人的私密视频| 露出奶头的视频| 国产成人精品久久二区二区91| 久久午夜亚洲精品久久| 真人做人爱边吃奶动态| 精品久久久久久,| 丝袜人妻中文字幕| 国产精品野战在线观看| 国产真人三级小视频在线观看| 欧美黑人巨大hd| 国产激情欧美一区二区| 亚洲国产精品sss在线观看| 男插女下体视频免费在线播放| 国产成人精品久久二区二区免费| 丁香欧美五月| 禁无遮挡网站| 黄色毛片三级朝国网站| 亚洲国产看品久久| 香蕉国产在线看| 欧美乱码精品一区二区三区| 男女下面进入的视频免费午夜| 亚洲成av人片免费观看| 免费在线观看日本一区| 久久久国产欧美日韩av| 夜夜爽天天搞| 啦啦啦韩国在线观看视频| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 啦啦啦韩国在线观看视频| 五月伊人婷婷丁香| 国产精华一区二区三区| 亚洲一区二区三区不卡视频| 在线永久观看黄色视频| 女警被强在线播放| 日韩国内少妇激情av| 久久天堂一区二区三区四区| 国产精品美女特级片免费视频播放器 | 久久香蕉激情| 国产伦在线观看视频一区| 国产91精品成人一区二区三区| 久久这里只有精品19| 曰老女人黄片| 床上黄色一级片| 亚洲avbb在线观看| 妹子高潮喷水视频| 天天添夜夜摸| 国产人伦9x9x在线观看| 国产亚洲精品第一综合不卡| 日本一二三区视频观看| 亚洲国产日韩欧美精品在线观看 | 久久人妻av系列| 亚洲 国产 在线| 一级片免费观看大全| 人人妻人人看人人澡| 久久精品国产清高在天天线| 国产区一区二久久| 欧美激情久久久久久爽电影| 亚洲精品美女久久av网站| 国产成人av教育| 手机成人av网站| 琪琪午夜伦伦电影理论片6080| 国内久久婷婷六月综合欲色啪| www.自偷自拍.com| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 久久中文看片网| 久久久国产成人精品二区| 人妻夜夜爽99麻豆av| 成人午夜高清在线视频| 看免费av毛片| 在线观看一区二区三区| 精品熟女少妇八av免费久了| 白带黄色成豆腐渣| 亚洲精品粉嫩美女一区| 国产精品美女特级片免费视频播放器 | 少妇的丰满在线观看| 三级国产精品欧美在线观看 | 啪啪无遮挡十八禁网站| 午夜日韩欧美国产| 国产私拍福利视频在线观看| 免费看日本二区| 大型黄色视频在线免费观看| 麻豆av在线久日| 在线观看舔阴道视频| 天堂√8在线中文| 国产av一区在线观看免费| 无遮挡黄片免费观看| 在线观看日韩欧美| 一二三四在线观看免费中文在| av在线天堂中文字幕| 看免费av毛片| 国产主播在线观看一区二区| 免费在线观看成人毛片| 香蕉久久夜色| 国产三级黄色录像| 天堂av国产一区二区熟女人妻 | 看片在线看免费视频| 脱女人内裤的视频| 又粗又爽又猛毛片免费看| 丝袜美腿诱惑在线| 看免费av毛片| 久久天堂一区二区三区四区| 中出人妻视频一区二区| 亚洲熟女毛片儿| 亚洲av熟女| 色综合婷婷激情| 五月玫瑰六月丁香| 国产激情久久老熟女| 午夜福利在线在线| 亚洲欧美日韩高清专用| 亚洲中文日韩欧美视频| 日本五十路高清| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷亚洲欧美| 露出奶头的视频| 中文字幕最新亚洲高清| 女同久久另类99精品国产91| 亚洲性夜色夜夜综合| 国产精品 欧美亚洲| 国产精品一区二区精品视频观看| 天天添夜夜摸| 怎么达到女性高潮| 日韩欧美一区二区三区在线观看| 少妇的丰满在线观看| av中文乱码字幕在线| 久久午夜亚洲精品久久| 在线十欧美十亚洲十日本专区| 亚洲性夜色夜夜综合| 丁香六月欧美| 亚洲男人天堂网一区| 亚洲国产中文字幕在线视频| 日韩成人在线观看一区二区三区| 在线观看免费午夜福利视频| 国产精品亚洲一级av第二区| 人成视频在线观看免费观看| 特级一级黄色大片| 午夜福利免费观看在线| 高潮久久久久久久久久久不卡| 丰满的人妻完整版| 成熟少妇高潮喷水视频| 国产三级在线视频| 热99re8久久精品国产| 每晚都被弄得嗷嗷叫到高潮| 午夜福利在线观看吧| 少妇熟女aⅴ在线视频| 欧美日韩亚洲综合一区二区三区_| 久久 成人 亚洲| 99精品欧美一区二区三区四区| 国产久久久一区二区三区| 亚洲国产精品sss在线观看| 精品无人区乱码1区二区| 黄色 视频免费看| 成人手机av| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 亚洲国产精品合色在线| 男女之事视频高清在线观看| 欧美乱码精品一区二区三区| 国产精品九九99| 一级毛片精品| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频| 午夜免费观看网址| 国产精品 欧美亚洲| www日本黄色视频网| 久久久久久亚洲精品国产蜜桃av| 欧美又色又爽又黄视频| 中文字幕精品亚洲无线码一区| 三级国产精品欧美在线观看 | 这个男人来自地球电影免费观看| 国产乱人伦免费视频| 最近最新中文字幕大全免费视频| 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀| 久久国产精品人妻蜜桃| 男女做爰动态图高潮gif福利片| 欧美中文综合在线视频| 在线观看一区二区三区| e午夜精品久久久久久久| 国产亚洲精品第一综合不卡| 成人国语在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 在线看三级毛片| 日本免费一区二区三区高清不卡| 国产精品美女特级片免费视频播放器 | 精品国内亚洲2022精品成人| 亚洲av成人av| 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| 丁香六月欧美| 午夜福利免费观看在线| 日韩欧美在线乱码| 国产精品精品国产色婷婷| 国产成人精品久久二区二区免费| 视频区欧美日本亚洲| 黑人操中国人逼视频| 丝袜美腿诱惑在线| 老司机深夜福利视频在线观看| 日日干狠狠操夜夜爽| 精品福利观看| 俄罗斯特黄特色一大片| 国产亚洲欧美在线一区二区| 一级毛片高清免费大全| 九九热线精品视视频播放| 9191精品国产免费久久| 久久久精品欧美日韩精品| 成人永久免费在线观看视频| 搡老妇女老女人老熟妇| 欧美性长视频在线观看| 国产成人精品无人区| 亚洲精品一卡2卡三卡4卡5卡| 国模一区二区三区四区视频 | 999久久久精品免费观看国产| 亚洲色图 男人天堂 中文字幕| 精品国产亚洲在线| 午夜精品在线福利| 我的老师免费观看完整版| 成人av在线播放网站| 嫩草影视91久久| 久久久久久亚洲精品国产蜜桃av| 18禁黄网站禁片免费观看直播| www.精华液| 国产精品久久久久久精品电影| 欧美久久黑人一区二区| 亚洲18禁久久av| 麻豆一二三区av精品| 听说在线观看完整版免费高清| 99国产精品一区二区三区| 欧美黑人巨大hd| 熟女电影av网| 1024手机看黄色片| 久久精品人妻少妇| 欧美色欧美亚洲另类二区| 成人三级黄色视频| 久久伊人香网站| 少妇被粗大的猛进出69影院| 国产人伦9x9x在线观看| 99在线人妻在线中文字幕| 看黄色毛片网站| 久久人妻福利社区极品人妻图片| 亚洲一区中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 不卡一级毛片| 狂野欧美白嫩少妇大欣赏| 亚洲美女黄片视频| 国产精品久久电影中文字幕| 久久中文看片网| 国产精品野战在线观看| 18禁黄网站禁片免费观看直播| 成人18禁在线播放| 亚洲精品色激情综合| 成人精品一区二区免费| 亚洲av第一区精品v没综合| 女同久久另类99精品国产91| 久久久国产成人免费| 午夜福利18| 日韩欧美在线二视频| 99热只有精品国产| 亚洲av中文字字幕乱码综合| 亚洲av日韩精品久久久久久密| 两性午夜刺激爽爽歪歪视频在线观看 | 国产视频内射| av欧美777| 91老司机精品| 精品久久久久久久毛片微露脸| 在线十欧美十亚洲十日本专区|