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

    可可西里鹽湖時空演化規(guī)律及漫溢成因分析

    2024-10-25 00:00:00苗學(xué)云米維軍唐占峰程佳沈朝政李勝威
    人民長江 2024年10期

    摘要:可可西里鹽湖漫溢水流給下游工程設(shè)施帶來了較大的安全隱患。通過對鹽湖遙感影像、所在區(qū)域氣溫、降雨、地溫,凍土上限變化,漫溢河流水位和流速的監(jiān)測數(shù)據(jù)分析,研究了鹽湖時空演化規(guī)律、漫溢成因及對下游潛在的安全風(fēng)險。研究表明:在空間上可可西里鹽湖呈“K-C-O”三階段變化特征,2022年鹽湖面積較2011年增加了3.73倍;時間上表現(xiàn)為明顯的“快-慢”兩階段增長期,快速期增長率為慢速期增長率的2.78倍,鹽湖向下游工程走廊漫溢速率為0.31 km/a;隨著季節(jié)性潮漲潮落及流速“峰谷”值變化,給下游橋梁樁基造成了一定的沖刷和熱侵蝕作用;鹽湖漫溢與多年凍土區(qū)氣溫、降雨、地溫等因素呈正相關(guān),近30 a氣溫升溫速率為0.052 ℃/a,年均降雨量線性增速為3.09 mm/a,地溫線性升溫速率為0.036 ℃/a,氣候暖濕化特征是造成鹽湖漫溢的主要成因。研究成果可為可可西里鹽湖下游工程防護(hù)措施制定提供理論與技術(shù)指導(dǎo)。

    關(guān) 鍵 詞:鹽湖漫溢;鹽湖面積;氣候變化;地溫;可可西里鹽湖;青藏高原

    中圖法分類號:P343.3

    文獻(xiàn)標(biāo)志碼:ADOI:10.16232/j.cnki.1001-4179.2024.10.017

    0 引 言

    青藏高原地域廣闊,分布著多種多樣的鹽湖。這些鹽湖大多深居高原腹地或高原盆地之中,以內(nèi)陸湖為主[1,這些湖泊處于動態(tài)變化中,受降雨增加影響,近年來,青海南部和西藏東北部地區(qū)的鹽湖處于穩(wěn)定擴(kuò)張狀態(tài)2。鹽湖的擴(kuò)張會導(dǎo)致高礦化度水流的增加,而高礦化度水流兼具腐蝕性和沖刷作用,對流經(jīng)區(qū)域的公路、鐵路、通訊線路、油氣管線、輸電線路等各類國家重大基礎(chǔ)設(shè)施的安全均會構(gòu)成威脅,也會對該地區(qū)生態(tài)環(huán)境、野生動植物產(chǎn)生不利影響[2-4。針對鹽湖面積持續(xù)性的增長,研究人員從不同角度進(jìn)行了探討:杜玉娥等5應(yīng)用多源衛(wèi)星資料對鹽湖進(jìn)行長期持續(xù)觀測,分析了鹽湖56 a來的變化趨勢;Zhang[6、Baban[7等應(yīng)用衛(wèi)星遙感影像提取了青藏高原湖泊精確的水位數(shù)據(jù),對鹽湖生態(tài)信息進(jìn)行了分類;袁康等[8利用2010~2018年青海省可可西里腹地鹽湖的CryoSat-2衛(wèi)星測高數(shù)據(jù)和Landsat遙感影像數(shù)據(jù),計算并構(gòu)建庫容關(guān)系曲線,結(jié)合氣候變化特征進(jìn)行了驅(qū)動力分析;Lu等[9利用衛(wèi)星遙感資料,建立了湖泊變化遙感綜合分析方法;Li等[10利用衛(wèi)星遙感數(shù)據(jù),分析了30多年來內(nèi)陸封閉流域內(nèi)湖泊變化的時空特征;謝昌衛(wèi)等[11根據(jù)野外考察資料對卓乃湖的潰決過程和潰決后的影響進(jìn)行了研究,并對下游鹽湖可能的潰決方式進(jìn)行了分析;李琳等12采用多源衛(wèi)星遙感技術(shù)對湖泊面積等水文要素進(jìn)行監(jiān)測,分析了湖泊面積與氣象要素、湖泊水量之間的相關(guān)性;劉文惠等13分析了鹽湖面積增大的趨勢,利用水庫潰壩預(yù)測分析模型,計算了鹽湖溢水潰壩洪峰流量,并提出了防治措施。

    國內(nèi)外學(xué)者們通過遙感衛(wèi)星影像分析,研究了鹽湖面積的變化趨勢,并對鹽湖決堤因素進(jìn)行了分析,取得了一系列重要的研究成果,但已有研究成果對可可西里鹽湖宏觀時空演化規(guī)律研究不夠充分,針對鹽湖漫溢成因也未能進(jìn)行相關(guān)研究工作。漫溢引起的高礦化度水流具有腐蝕性,且由于持續(xù)增加的水流勢必對下游道路工程樁基產(chǎn)生沖刷作用,具有較大的安全隱患,亟需在宏觀維度持續(xù)監(jiān)測鹽湖的變化過程。由于可可西里鹽湖地處青藏高原多年凍土區(qū)腹地,深受外界環(huán)境變遷的影響,近年來青藏高原地區(qū)升溫現(xiàn)象明顯,相關(guān)數(shù)據(jù)顯示,近45 a來升溫幅度達(dá)到了2.9 ℃,年均升溫速率達(dá)到了0.052 ℃/a[14。因多年凍土具有溫度敏感性,環(huán)境溫濕度的加劇,給周邊鹽湖堤壩穩(wěn)定性造成了較大的安全隱患,可可西里鹽湖漫溢與環(huán)境要素變化存在必然聯(lián)系,因此需要從環(huán)境變化角度進(jìn)行分析。

    可可西里鹽湖下游為國家進(jìn)藏大動脈,是中國多年凍土區(qū)重要的交通樞紐地帶,開展鹽湖漫溢時空演化規(guī)律研究,以及深入探討鹽湖漫溢的成因,對于保障工程廊道、能源通道的安全性具有十分重要的研究意義?;诖?,本文收集青藏鐵路沿線風(fēng)火山觀測站、格爾木市氣象局實測資料,從環(huán)境氣溫、降水、凍土上限變化等因素,揭示鹽湖漫溢的根本成因,以期為下游工程防護(hù)措施設(shè)計提供參考。

    1 可可西里鹽湖基本概況

    可可西里鹽湖(以下簡稱鹽湖)位于青海省玉樹藏族自治州治多縣西北五道梁鄉(xiāng)西部,位于東經(jīng)93°42′、北緯35°52′,海拔4 500 m,屬于高寒凍土帶寒冷氣候,湖區(qū)附近是荒漠沼澤草原地帶。根據(jù)2022年衛(wèi)星遙感影像分析可知,可可西里鹽湖東西長19.40 km,寬15.06 km,面積191.26 km2,湖盆受大氣降水和地表水補(bǔ)給,湖水化學(xué)成分主要為硫酸鹽型、硫酸鎂亞型,湖水含鹽量偏低[15-16。

    鹽湖下游與工程廊道最近距離為10.63 km,區(qū)域內(nèi)平坦開闊,坡降不足1‰,如圖1所示,受2011年上游卓乃湖決堤影響,湖面面積持續(xù)增大,鹽湖水位通過沖溝向下游清水河開始漫溢,水流順著清水河河道將會對下游道路工程產(chǎn)生較大安全隱患。

    2 數(shù)據(jù)來源及監(jiān)測方案

    2.1 鹽湖面積數(shù)據(jù)處理

    本次研究選取了2011~2022年鹽湖面積衛(wèi)星遙感影像,圖示比例為1∶300,插入AutoCAD軟件,將測量圖像組合形成面域,用area命令法進(jìn)行計算,可求得鹽湖面積;同時應(yīng)用比例系數(shù)法可求得鹽湖至工程走廊的距離。

    2.2 氣象數(shù)據(jù)來源與處理

    中鐵西北科學(xué)研究院風(fēng)火山觀測站地處青藏高原多年凍土區(qū)腹地,海拔4 760 m,青海省五道梁氣象站海拔4 600 m,是青海省地區(qū)重要的高原寒區(qū)氣象站,可可西里鹽湖位置介于五道梁與風(fēng)火山地區(qū)之間,鹽湖的氣象要素采用風(fēng)火山與五道梁觀測站數(shù)據(jù)的均值處理,氣象數(shù)據(jù)采用DZZ4型全自動氣象站。地溫監(jiān)測數(shù)據(jù)來源于青藏線多年凍土區(qū)監(jiān)測數(shù)據(jù),現(xiàn)場布設(shè)深孔地溫監(jiān)測孔,孔內(nèi)布設(shè)熱敏電阻串,地溫孔布設(shè)示意如圖2所示。

    2.3 漫溢河流監(jiān)測方案

    為監(jiān)測漫溢河流的變化情況,分別在河流橋墩及上下游處布設(shè)河流水位監(jiān)測儀,流速儀采用全自動雷達(dá)流速儀,共布設(shè)3套監(jiān)測設(shè)備,分別布置于橋墩的上游、下游及墩臺處,編號S1、S2、S3。具體布置見圖3。

    3 鹽湖漫溢時空變化特征分析

    3.1 空間變化分析

    3.1.1 上游湖泊面積變化分析

    圖4為2011年卓乃湖決堤前后可可西里鹽湖及上游3個湖泊(卓乃湖、庫賽湖、海丁諾爾湖)的面積變化情況。2011年決堤前,最上游卓乃湖面積為270.91 km2,決堤發(fā)生后,卓乃湖在水力梯度作用下向下游漫溢,卓乃湖面積逐漸減小,2017年12月面積變化為154.09 km2,而鹽湖面積增大為150.21 km2,卓乃湖下游的庫賽湖面積也逐步擴(kuò)大,由288.89 km2增大至327.46 km2,到2022年12月鹽湖面積達(dá)到180.66 km2。由于在4個湖泊中,鹽湖海拔最低,在水力梯度作用下,決堤后鹽湖成為上游水量最主要的承載體,隨著鹽湖面積的逐步擴(kuò)大,其溢流勢必會向下游工程廊道漫溢。

    3.1.2 鹽湖面積變化情況

    截選鹽湖2011年11月、2017年及2022年12月的衛(wèi)星遙感影像進(jìn)行疊加,如圖5所示。從圖5可以清晰地看到3個階段鹽湖面積的變化情況,隨著上游湖泊伴隨水力梯度往下游漫溢,鹽湖的面積逐步增大,向四周逐漸擴(kuò)展,由于東南方向地勢較低,水流向低洼地區(qū)流動,導(dǎo)致東南方向成為了漫溢的主要方向,漫溢的河流已經(jīng)與下游的清水河產(chǎn)生了通道,勢必會增加清水河流量,從而對下游工程a58ln/8QxD3WVBFLTpUrTf+OBAQNLlIzXc8DYozEU6M=廊道橋梁樁基工程造成了沖刷和腐蝕作用,給橋梁帶來一定的安全風(fēng)險。

    由圖6鹽湖面積變化分析可知,鹽湖在空間上具有明顯的向下游(東南方向)漫溢變化的過程,上游卓乃湖決堤前鹽湖所示面積近似“K”形,決堤后鹽湖在水力梯度作用下逐漸擴(kuò)大,至2012年12月鹽湖產(chǎn)生了明顯的外溢,向4個方位均有擴(kuò)大的趨勢,鹽湖初步形成“C”形,隨著鹽湖上游水分的持續(xù)補(bǔ)給,鹽湖的面積進(jìn)一步擴(kuò)大,至2017年12月,鹽湖形成了封閉的“O”形,鹽湖擴(kuò)大的趨勢開始向下游漫延,至2022年12月鹽湖的面積相較于2017年存在進(jìn)一步擴(kuò)張的趨勢,并與下游清水河形成了一條水力通道,鹽湖的水流對下游橋梁樁基產(chǎn)生了一定的沖刷作用,且鹽湖腐蝕性水質(zhì)也會對樁基產(chǎn)生鹽蝕作用,對下游橋梁帶來了一定的安全風(fēng)險。

    3.2 時間變化分析

    圖7為鹽湖面積歷時變化,鹽湖是卓乃湖、庫賽湖、海丁諾爾湖的下游湖,地勢較低,由圖6鹽湖面積變化來看,1989~2022年呈明顯的增加趨勢,由上游卓乃湖決堤前的48.45 km2,增加至2022年的180.66 km2,2011~2022年間面積增大速率為12.02 km2/a,面積累計增加了132.21 km2,增大率為272.88%,面積增大為決堤前的3.73倍。

    在卓乃湖決堤前后鹽湖面積產(chǎn)生了較大的變化,2011~2017年為快速增長期,2011年上游湖決堤前為48.45 km2,2017年面積增大為150.21 km2,鹽湖年平均增長率為16.96%;2017~2022年為慢速增長期,由2017年的150.21 km2,增大為2022年的180.66 km2,鹽湖年均面積增長率為6.09%,快速增長期是慢速增長期的2.78倍。上述兩階段明顯變化是因為在上游湖泊水力梯度作用下,漫溢的水量較大,且鹽湖四周地勢較為平坦,決堤后水位迅速向四周漫溢的趨勢明顯,表現(xiàn)為鹽湖面積較快的增長;當(dāng)上游湖泊水位逐漸達(dá)到了新的水頭平穩(wěn),且鹽湖漫溢河流通過水力通道向下游漫溢,鹽湖的面積增長趨勢有所減緩,除此之外,通過人工應(yīng)急工程干預(yù),增設(shè)了漫溢河流通道的導(dǎo)流墻,一定程度上起到了減緩面積進(jìn)一步擴(kuò)展的作用。

    3.3 邊界漫溢速率

    受上游湖泊決堤影響,可可西里鹽湖的面積逐步擴(kuò)大,呈現(xiàn)明顯的漫溢效應(yīng),可可西里鹽湖漫溢隨著時間的推移,逐步向下游清水河大橋趨近發(fā)展,由圖8可知,2011年可可西里鹽湖距離清水河大橋最小距離為13.74 km,隨著可可西里鹽湖面積的增大,其最小間距逐步縮小,至2022年距離為10.38 km,距離減小了3.36 km,表明可可西里鹽湖向清水河大橋漫溢了3.36 km,年均推進(jìn)0.31 km,最小距離相較于2011年,減小了24.46%。

    4 漫溢河流水文變化分析

    4.1 河流水位

    水位變化情況見圖9,根據(jù)水位變化情況可以看出,該處河流流水為常流水,即一年中均處于有水流狀態(tài)。3處水位監(jiān)測點變化較為一致,2021年7~12月,水位較為穩(wěn)定,水位值在0.3~0.4 m區(qū)間,2021年12月至2022年3月,水位呈抬升趨勢發(fā)展,由0.3 m上升至0.8 m,最高值為0.83 m,2~5月水位處于回落狀態(tài),5月以后則又開始緩慢抬升,3處監(jiān)測點的常年水位平均值為0.33~0.39 m。常年有水會對地基多年凍土造成較大的熱干擾,同時水流的季節(jié)性漲落對橋梁墩臺有一定的沖刷作用。

    4.2 河流流速

    圖10為河水流速變化,從圖中可以看到水流基本呈上升-下降兩階段趨勢發(fā)展,自2022年7月監(jiān)測以來,水流流速逐漸增大,S1監(jiān)測點由7月的0.18 m/s逐漸上升至10月的0.52 m/s,流速增大了0.34 m/s,S2監(jiān)測點由0.09 m/s增大至0.34 m/s,增大了0.25 m/s,此后水流流速呈下降趨勢發(fā)展,至2022年11月流速分別為0.50 m/s和0.30 m/s,河流流速的變化對兩側(cè)多年凍土層路基造成了較大影響和熱干擾,且對清水河大橋樁基產(chǎn)生一定的侵蝕和沖刷作用。

    5 鹽湖漫溢因素相關(guān)性分析

    可可西里鹽湖漫溢的直接原因是上游卓乃湖決堤,潰堤后水流迅速向下游聚集,由于可可西里鹽湖處于較低地勢,成為上游湖泊溢流最終接受體,鹽湖面積呈持續(xù)增加趨勢。而上游湖泊的決堤與青藏高原環(huán)境變遷具有密不可分的關(guān)系,可可西里地處多年凍土腹地,是青藏高原環(huán)境氣候變化的“敏感器”。受氣候暖濕化影響,多年凍土退化,活動層厚度持續(xù)加深[17-20,地層融沉病害加劇,原有湖泊湖堤穩(wěn)定性受到擾動,處于不穩(wěn)定狀態(tài),加之降雨量增加,湖面水位持續(xù)升高,成為決堤的關(guān)鍵因素。本節(jié)就青藏高原風(fēng)火山和五道梁氣象站的溫度、降雨、地溫與鹽湖漫溢相關(guān)性進(jìn)行分析。

    5.1 鹽湖漫溢與氣溫變化相關(guān)性分析

    圖11為鹽湖面積與氣溫均值相關(guān)性曲線圖,氣溫數(shù)據(jù)為風(fēng)火山與五道梁氣象站監(jiān)測數(shù)據(jù)平均處理得到。從圖中明顯可以看到,鹽湖面積與氣溫具有正相關(guān)性,相關(guān)性系數(shù)為0.53。由變化趨勢看,氣溫呈明顯的波動式上升過程,近30 a來,氣溫由-5.90 ℃上升至-4.25 ℃,升溫速率為0.052 ℃/a;在2011年前,鹽湖面積基本處于較低值,均值為39.8 km2,上游湖潰決后面積迅速上升,達(dá)到了180.66 km2;鹽湖在潰決前,年均氣溫處于高值,約為-4.20 ℃,比1990~2023年平均值高出0.80 ℃,成為鹽湖漫溢的直接原因之一。隨著2022年搶險應(yīng)急工程的開展,漫溢河道兩旁設(shè)置了導(dǎo)流堤和河床鋪砌,將漫溢的水流引入到清水河中,減緩了鹽湖進(jìn)一步擴(kuò)張的趨勢,鹽湖面積有所減小,至2023年12月顯示,鹽湖面積穩(wěn)定在133.56 km2。

    隨著環(huán)境溫度的升高,勢必會對地層溫度造成直接影響,從而引起地溫冷儲量逐漸散失,活動層厚度持續(xù)增加,地層融沉病害直接影響了湖堤的穩(wěn)定性??煽晌骼锏貐^(qū)其他湖泊還存在著不同程度的潰壩風(fēng)險,應(yīng)積極排查上游鹽湖堤壩防護(hù)狀態(tài),進(jìn)一步提高沿岸防護(hù)措施等級。

    5.2 鹽湖漫溢與降雨變化相關(guān)性分析

    圖12為可可西里鹽湖面積變化與降雨量相關(guān)性曲線圖,由圖可知,降雨量與鹽湖漫溢具有正相關(guān)性,相關(guān)性系數(shù)為0.47,介于0.3~0.5,屬于中等相關(guān)性。

    從雨量分布來看,除1989年降雨量較高外,1990~2022年,降雨量基本處于明顯的上升狀態(tài),年均降雨量線性增速為3.09 mm/a,以30 mm/10a的幅度增加,其中2010~2014年降雨量處于高位,平均降雨量為943.9 mm,比1990~2010年平均降雨量(874.4 mm)多出69.5 mm,降雨量持續(xù)增加成為可可西里地區(qū)湖泊的重要水源補(bǔ)給,雨量的增加也對多年凍土造成了熱侵蝕。

    5.3 鹽湖漫溢與地溫相關(guān)性分析

    從圖13鹽湖面積與地溫相關(guān)性可知,自1990年以來,地溫的升高與鹽湖漫溢具有較強(qiáng)的正相關(guān)性,相關(guān)系數(shù)達(dá)到了0.9,表現(xiàn)為強(qiáng)相關(guān)性。從地溫變化來看,自2005年以后,升溫速率明顯加快,年均地溫達(dá)到了-2.65 ℃,較近30多年平均地溫高出0.38 ℃,線性升溫速率為0.036 ℃/a。同時期,鹽湖也進(jìn)入一個快速增長期,兩者之間具有相同的變化規(guī)律,鹽湖面積擴(kuò)張了3.73倍。

    隨著地溫的不斷升高,多年凍土區(qū)退化態(tài)勢由南北兩界逐步向可可5jWEEvVhUA6yUzbqI7hVUQ==西里地區(qū)發(fā)展[21,可可西里地區(qū)升溫現(xiàn)象更為顯著。地溫的升高會加劇湖岸凍土穩(wěn)定性的減弱,地表徑流及地下水成為鹽湖重要的補(bǔ)給通道,一方面增加了鹽湖水量,另一方面對堤壩工程造成一定的影響,給上游湖泊安全防護(hù)帶來較大的安全隱患。

    6 討 論

    由鹽湖面積的擴(kuò)張過程來看,自2005年前后開始有較大的升溫過程,氣溫升高直接導(dǎo)致周邊冰川消融、積雪融化、凍土圈逐步退化等一系列連鎖反應(yīng),整個青藏高原地區(qū)地表水系處于較為活躍時期,地表徑流和降雨可能成為鹽湖重要的水量補(bǔ)給。隨著氣溫的升高,凍土層逐步融化,活動層厚度增加,融化土層一定程度影響了湖泊沿岸堤壩工程的穩(wěn)定性,成為鹽湖漫溢因素之一。2010年前后,降雨量達(dá)到高值,為429~576 mm,高于統(tǒng)計以來的平均值336~382 mm,這段時期的降雨量增加與鹽湖漫溢具有一致性,雨量的增加一方面補(bǔ)給了鹽湖水位,另一方面雨水侵蝕下給地層凍土造成了熱侵蝕,進(jìn)一步增加了堤壩失穩(wěn)的風(fēng)險。隨著氣候暖濕化加劇,地溫呈明顯的升溫狀態(tài),凍土層凍結(jié)水的能力持續(xù)減弱,也是造成鹽湖擴(kuò)張的一個關(guān)鍵因素。相關(guān)性系數(shù)計算結(jié)果表明,地溫升高是最為關(guān)鍵的因素,其次是氣溫升高,再次為降雨量的增加。顯然鹽湖擴(kuò)張是多方面因素綜合作用的結(jié)果,今后需要進(jìn)一步加強(qiáng)當(dāng)?shù)貧庀?、水文、地溫的監(jiān)測。

    鹽湖漫溢的下泄將會嚴(yán)重影響下游工程廊道的系統(tǒng)性風(fēng)險。鹽湖水體含鹽量高、腐蝕性強(qiáng),一旦可可西里鹽湖大量下泄,將會對地基多年凍土和橋梁造成影響,如由于水量增加,水流的長期沖刷及熱侵蝕會引起河道下切,進(jìn)而引起地基凍土融化和水土流失。根據(jù)現(xiàn)場調(diào)查,由于水流沖刷,兩岸河堤塌陷嚴(yán)重,最高處已超過2 m。同時,水流流經(jīng)橋下時對橋墩形成沖刷,下泄高礦化度湖水不僅會進(jìn)一步加劇凍土融化,同時對橋墩及基礎(chǔ)混凝土產(chǎn)生腐蝕作用,這些作用會導(dǎo)致橋墩不均勻沉降、傾斜、混凝土剝蝕等病害。

    隨著氣候進(jìn)一步暖濕化,鹽湖再次漫溢風(fēng)險呈持續(xù)增加的態(tài)勢,應(yīng)加強(qiáng)當(dāng)?shù)貧庀?、地溫、河道水量的監(jiān)測工作,同時對下游工程樁基、梁體增加長期監(jiān)測系統(tǒng),持續(xù)監(jiān)測漫溢河流變化、橋梁樁基位移變形以及樁基基礎(chǔ)地溫變化,為多年凍土區(qū)工程廊道安全運營提供技術(shù)保障。

    7 結(jié) 論

    本文應(yīng)用衛(wèi)星遙感影像分析,揭示了可可西里鹽湖漫溢的時空演化規(guī)律,并根據(jù)鹽湖附近風(fēng)火山和五道梁氣象站氣溫、降雨、地溫實測數(shù)據(jù),綜合分析了可可西里鹽湖漫溢的成因,并提出了可可西里鹽湖漫溢造成的下游潛在危害。

    (1)可可西里鹽湖在空間上具有明顯的向東南方向漫溢變化的過程,可可西里鹽湖面積分為“K-C-O”形三階段變化過程,可可西里鹽湖的面積由44.83 km2增大至180.66 km2,面積增大為3.73倍。在時間上具有兩階段增長期過程,2011~2017年為快速增長期,可可西里鹽湖年平均增長率為16.69%;2017~2022年為慢增長期,可可西里鹽湖年均面積增長率為6.09%。

    (2)可可西里鹽湖隨著時間的推移,逐步向下游工程走廊漫溢。2011年可可西里鹽湖距離清水河大橋最小距離為13.74 km,隨著可可西里鹽湖面積的增大,其最小間距逐步縮小,至2022年距離為10.38 km,距離減小了3.36 km,表明可可西里鹽湖向清水河大橋漫溢了3.36 km,年均推進(jìn)0.31 km,最小距離相較于2011年,減小了24.46%。

    (3)地溫、氣溫和降雨是鹽湖擴(kuò)張的重要因素,青藏高原近年來溫度呈現(xiàn)逐漸升高的趨勢,環(huán)境溫度的升高,加劇了多年凍土區(qū)凍土圈的融化,上游下泄的水量一旦超越可可西里鹽湖承載水量能力,湖泊進(jìn)一步漫溢的風(fēng)險將呈增加的趨勢。

    (4)鹽湖下泄河流隨著季節(jié)性漲落變化,不僅會進(jìn)一步加劇凍土融化,同時會對橋墩及基礎(chǔ)混凝土產(chǎn)生腐蝕作用。這些作用將導(dǎo)致橋墩不均勻沉降、傾斜、混凝土剝蝕等病害。

    參考文獻(xiàn):

    [1] 孔維剛,卜令忠,鄭綿平,等.青藏高原鹽湖觀測站概述[J].科技導(dǎo)報,2017,35(6):103-107.

    [2] 閆立娟,鄭綿平,齊路晶.青藏高原湖泊湖面變遷及影響因素[J].科技導(dǎo)報,2017,35(6):83-88.

    [3] 劉智琦,潘保柱,韓谞,等.青藏高原湖泊水環(huán)境特征及水質(zhì)評價[J].環(huán)境科學(xué),2022,43(11):5073-5083.

    [4] 黃建蓉.鹽度對青藏高原湖泊微生物群落結(jié)構(gòu)與功能穩(wěn)定性的影響[D].武漢:中國地質(zhì)大學(xué),2021.

    [5] 杜玉娥,劉寶康,賀衛(wèi)國,等.1976~2017年青藏高原可可西里鹽湖面積動態(tài)變化及成因分析[J].冰川凍土,2018,40(1):47-54.

    [6] ZHANG G Q,XIE H J,KANG S C,et al.Monitoring lake level chan-ges on the Tibetan Plateau using ICESat altimetry data(2003~2009)[J].Remote Sensing of Environment,2011,115(7):1733-1742.

    [7] BABAN S M J.Use of remote sensing and geographical information systems in developing lake management strategies[J].Hydrobiologia,1999,395:211-226.

    [8] 袁康,譚德寶,趙靜,等 近十年可可西里鹽湖水量變化及其影響因素分析[J].人民長江,2022,53(5):111-117.

    [9] LU A X,WANG L H,YAO T D.Study on remote sensing method of modern change of Qinghai-Tibet Plateau lakes[J].Remote Sensing Technology and Application,2006,21(3):174-177.

    [10]LI J L,SHENG Y W,LUO J C,et al.Remote sensing mapping of inland lakes change in the Tibetan Plateau[J].Journal of Lake Sciences,2011,23(3):311-320.

    [11]謝昌衛(wèi),張鈺鑫,劉文惠,等.可可西里卓乃湖潰決后湖區(qū)環(huán)境變化及鹽湖可能的潰決方式[J].冰川凍土,2020,42(4):1344-1352.

    [12]李琳,譚德寶,文雄飛,等.氣候變化對可可西里鹽湖流域湖泊水量變化的影響分析[J].長江科學(xué)院院報,2022,39(10):16-23.

    [13]劉文惠,謝昌衛(wèi),王武,等.青藏高原可可西里鹽湖水位上漲趨勢及潰決風(fēng)險分析[J].冰川凍土,2019,41(6):1467-1474.

    [14]苗學(xué)云,李永強(qiáng),米維軍,等.青藏線高原凍土風(fēng)火山綜合觀測研究[R].蘭州:中鐵西北科學(xué)研究院有限公司,2023.

    [15]YAO X J,SUN M P,GONG P,et al.Overflow probability of the salt lake in Hoh Xil Region[J].Journal of Geographical Sciences,2018,28(5):647-655.

    [16]劉勇平,周敬,韓鳳清,等.青??煽晌骼飽|部鹽湖水化學(xué)及沉積特征初步研究[J].鹽湖研究,2009,17(3):10-16.

    [17]趙海鵬.青藏高原多年凍土區(qū)土壤水熱變化及其影響因素研究[D].蘭州:蘭州大學(xué),2021.

    [18]王銀學(xué),趙林,李韌,等.影響多年凍土上限變化的因素探討[J].冰川凍土,2011,33(5):1064-1067.

    [19]韓龍武,楊印海,戴路生,等.青藏鐵路沿線填土路基下多年凍土的演化規(guī)律[J].中國安全科學(xué)學(xué)報,2018,28(增2):6-10.

    [20]楊永鵬,孟進(jìn)寶,韓龍武,等.青藏鐵路路基下多年凍土演化特征及規(guī)律研究[J].鐵道建筑,2018,58(5):81-86.

    [21]孟超,韓龍武,趙相卿,等.氣溫持續(xù)升高對青藏鐵路運輸安全的影響研究[J].中國安全科學(xué)學(xué)報,2018,28(增2):1-5.

    (編輯:謝玲嫻)

    Analysis on spatiotemporal evolution law and overflow causes of Salt Lake in Hoh XilMIAO Xueyun1,2,MI Weijun1,2,TANG Zhanfeng3,CHENG Jia1,2,SHEN Chaozheng2,4,LI Shengwei2,4

    (1.Northwest Research Institute Co.,Ltd.of CREC,Lanzhou 730000,China; 2.Qinghai Province Key Laboratory of Frozen Soil and Environment Engineering,Geermu 816000,China; 3.Engineering Management Office of China Railway Lanzhou Bureau Group Co.,Ltd.,Lanzhou 730002,China; 4.College of Transportation Engineering of Nanjing Tech,Nanjing Tech University,Nanjing 210009,China)

    Abstract: The overflowing water of the Salt Lake in Hoh Xil has brought great potential safety hazards to the downstream rail transit facilities.We analyzed monitoring data of remote sensing images over the Salt Lake in Hoh Xil,temperature,rainfall,ground temperature and the upper limit changes of permafrost in the area,as well as the water level and flow velocity of the overflowing river.Based on the above data,the spatiotemporal evolution law of the Salt Lake in Hoh Xil,the causes of overflow,and the potential safety risks to the downstream were studied.The research showed that spatially,the Salt Lake in Hoh Xil presented a three-stage change characteristic of "K-C-O".In 2022,the area of the Salt Lake in Hoh Xil was 3.73 times larger than that in 2011.Temporally,it showed an obvious two-stage growth period of "fast-slow".The growth rate in the fast period was 2.78 times that in the slow period.The overflow rate of the Salt Lake in Hoh Xil to the downstream engineering corridor was 0.31 km/a.With the seasonal ebb and flow and the change of the "peak-valley" value of flow velocity,it has caused certain scouring and thermal erosion effects on the bridge pile foundation downstream.The overflow of the Salt Lake in Hoh Xil was positively correlated with factors such as temperature,rainfall,and ground temperature in the permafrost region.In the past 30 years,the temperature increase rate was 0.052 °C/a,the linear growth rate of annual average rainfall was 3.09 mm/a,and the linear temperature increase rate of ground temperature was 0.036 °C/a.The characteristics of warm and humid climate were the main causes of Salt Lake in Hoh Xil overflow.The research results can provide theoretical and technical guidance for downstream facility protection.

    Key words: overflow of salt lake;area of salt lake;climate change;ground temperature;Salt Lake in Hoh Xil;Qinghai-Tibet Plateau

    小说图片视频综合网站| 卡戴珊不雅视频在线播放| 春色校园在线视频观看| 免费看a级黄色片| 成人毛片a级毛片在线播放| 又粗又爽又猛毛片免费看| 欧美又色又爽又黄视频| 久久99蜜桃精品久久| 五月伊人婷婷丁香| 成人国产麻豆网| 高清日韩中文字幕在线| 成年女人永久免费观看视频| 欧美一区二区亚洲| 免费人成视频x8x8入口观看| 亚洲国产精品成人久久小说 | 国产爱豆传媒在线观看| 久久久久久久午夜电影| 中文欧美无线码| 悠悠久久av| 色吧在线观看| 免费观看的影片在线观看| 久久久久久久久久久丰满| 国产成人freesex在线| 欧美日本亚洲视频在线播放| 18禁在线播放成人免费| 秋霞在线观看毛片| 久久久久免费精品人妻一区二区| 国产 一区精品| 中国国产av一级| 国产麻豆成人av免费视频| 午夜精品一区二区三区免费看| 日本成人三级电影网站| 亚州av有码| 欧美xxxx黑人xx丫x性爽| 校园春色视频在线观看| .国产精品久久| 亚洲自拍偷在线| 午夜福利在线在线| 亚洲成人精品中文字幕电影| 亚洲中文字幕一区二区三区有码在线看| ponron亚洲| 免费观看a级毛片全部| 免费电影在线观看免费观看| 男女边吃奶边做爰视频| 搡老妇女老女人老熟妇| 高清毛片免费观看视频网站| 麻豆精品久久久久久蜜桃| 亚洲图色成人| 精品国内亚洲2022精品成人| 一个人看视频在线观看www免费| 三级经典国产精品| 精品一区二区免费观看| videossex国产| 亚洲图色成人| 成人毛片60女人毛片免费| 亚洲18禁久久av| 国产亚洲精品久久久com| 国产大屁股一区二区在线视频| 国产高清不卡午夜福利| 国产精品一区二区三区四区久久| 欧美xxxx性猛交bbbb| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 校园人妻丝袜中文字幕| 成人亚洲精品av一区二区| 日韩,欧美,国产一区二区三区 | 日韩三级伦理在线观看| 狂野欧美白嫩少妇大欣赏| 国产av麻豆久久久久久久| 91久久精品电影网| 一级二级三级毛片免费看| 色视频www国产| 中文资源天堂在线| 成人三级黄色视频| 午夜老司机福利剧场| 亚洲美女搞黄在线观看| 身体一侧抽搐| 久久久久久久午夜电影| 桃色一区二区三区在线观看| 国产精品国产高清国产av| 少妇丰满av| 嫩草影院入口| 啦啦啦韩国在线观看视频| 午夜福利高清视频| 青青草视频在线视频观看| 一级av片app| 少妇熟女aⅴ在线视频| 男女下面进入的视频免费午夜| 性插视频无遮挡在线免费观看| 99久久人妻综合| 熟女电影av网| 久久久久久国产a免费观看| а√天堂www在线а√下载| 麻豆精品久久久久久蜜桃| 国产成人精品一,二区 | 亚洲欧美日韩高清专用| 亚洲无线观看免费| 色综合亚洲欧美另类图片| 午夜福利视频1000在线观看| 久久久欧美国产精品| 久久99蜜桃精品久久| 日日干狠狠操夜夜爽| 国产精品久久电影中文字幕| 成人毛片60女人毛片免费| 丰满乱子伦码专区| 国产日本99.免费观看| 色5月婷婷丁香| av国产免费在线观看| 欧美激情在线99| 99久久成人亚洲精品观看| 免费观看在线日韩| 成人漫画全彩无遮挡| 99久国产av精品| 少妇熟女aⅴ在线视频| 99久久精品热视频| 人妻少妇偷人精品九色| 久久鲁丝午夜福利片| 国产高清激情床上av| 免费观看的影片在线观看| 国产av麻豆久久久久久久| 亚洲真实伦在线观看| 天美传媒精品一区二区| 尤物成人国产欧美一区二区三区| 亚洲精品乱码久久久久久按摩| 欧美性感艳星| 99热这里只有精品一区| 日韩欧美 国产精品| a级毛片a级免费在线| 免费av不卡在线播放| 我的女老师完整版在线观看| 亚洲成人久久爱视频| 中国美白少妇内射xxxbb| 深夜a级毛片| 国产黄色小视频在线观看| 亚洲av不卡在线观看| 亚洲人成网站在线播放欧美日韩| 有码 亚洲区| 久久99热6这里只有精品| 看免费成人av毛片| 欧美精品国产亚洲| 亚洲精品456在线播放app| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 亚洲人与动物交配视频| 中文资源天堂在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合色国产| 男女边吃奶边做爰视频| 亚洲不卡免费看| 国产女主播在线喷水免费视频网站 | a级一级毛片免费在线观看| 中文字幕av成人在线电影| 少妇裸体淫交视频免费看高清| 国产精品嫩草影院av在线观看| 真实男女啪啪啪动态图| 99久久久亚洲精品蜜臀av| 亚洲欧洲国产日韩| 日韩 亚洲 欧美在线| 不卡视频在线观看欧美| 国产av不卡久久| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 成人二区视频| 国产免费一级a男人的天堂| 日韩强制内射视频| 级片在线观看| 亚洲自拍偷在线| 亚洲精品久久国产高清桃花| av在线亚洲专区| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 中文资源天堂在线| 99热网站在线观看| 国产黄色视频一区二区在线观看 | eeuss影院久久| 日韩高清综合在线| 不卡视频在线观看欧美| 有码 亚洲区| 国产v大片淫在线免费观看| 在线天堂最新版资源| 久久久久久大精品| 中文字幕人妻熟人妻熟丝袜美| 男女啪啪激烈高潮av片| 日韩高清综合在线| 亚洲精华国产精华液的使用体验 | 99热精品在线国产| 九九久久精品国产亚洲av麻豆| 99在线视频只有这里精品首页| 精品少妇黑人巨大在线播放 | 给我免费播放毛片高清在线观看| 国产精品不卡视频一区二区| 免费不卡的大黄色大毛片视频在线观看 | 九色成人免费人妻av| 女人被狂操c到高潮| 国产亚洲精品久久久久久毛片| 国产一区二区三区在线臀色熟女| 国产高清三级在线| eeuss影院久久| 国产真实乱freesex| 一个人看视频在线观看www免费| 国产精品,欧美在线| 欧美xxxx性猛交bbbb| 网址你懂的国产日韩在线| 变态另类丝袜制服| 国产精品一二三区在线看| 小说图片视频综合网站| 狠狠狠狠99中文字幕| 国内精品久久久久精免费| 哪个播放器可以免费观看大片| 国内精品宾馆在线| 成人高潮视频无遮挡免费网站| 深夜精品福利| 中文字幕免费在线视频6| 两个人的视频大全免费| 国产精品久久视频播放| 国产探花极品一区二区| 九色成人免费人妻av| 亚洲人成网站在线播| 欧美三级亚洲精品| 久久精品国产清高在天天线| 毛片一级片免费看久久久久| 欧美极品一区二区三区四区| 此物有八面人人有两片| 好男人在线观看高清免费视频| 不卡视频在线观看欧美| 嫩草影院精品99| 又爽又黄a免费视频| 只有这里有精品99| 欧美色视频一区免费| 九草在线视频观看| 黄色视频,在线免费观看| 校园春色视频在线观看| 亚洲欧美日韩无卡精品| 非洲黑人性xxxx精品又粗又长| 一级黄色大片毛片| 免费看美女性在线毛片视频| 看片在线看免费视频| 丰满的人妻完整版| 青春草国产在线视频 | 欧美zozozo另类| 99久久精品一区二区三区| 成人国产麻豆网| 麻豆乱淫一区二区| 边亲边吃奶的免费视频| 亚洲精品国产av成人精品| 国产成人一区二区在线| 淫秽高清视频在线观看| 老女人水多毛片| 91aial.com中文字幕在线观看| 在现免费观看毛片| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 久久精品国产亚洲av天美| 日日啪夜夜撸| 国产亚洲av嫩草精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 日韩视频在线欧美| av免费在线看不卡| 免费黄网站久久成人精品| 2022亚洲国产成人精品| 久99久视频精品免费| 国产高清不卡午夜福利| 男女视频在线观看网站免费| 欧美激情久久久久久爽电影| 激情 狠狠 欧美| 久久久a久久爽久久v久久| www日本黄色视频网| 亚洲国产色片| 国产欧美日韩精品一区二区| 国产色爽女视频免费观看| 直男gayav资源| 日本熟妇午夜| 97人妻精品一区二区三区麻豆| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 热99在线观看视频| 最后的刺客免费高清国语| 国产高清三级在线| 国产精品久久久久久亚洲av鲁大| 狠狠狠狠99中文字幕| 国产大屁股一区二区在线视频| 男人舔奶头视频| 狂野欧美激情性xxxx在线观看| 少妇熟女欧美另类| 伦精品一区二区三区| 久久久久久久久大av| 久久久久性生活片| 亚洲经典国产精华液单| 哪里可以看免费的av片| 日韩一区二区三区影片| 亚洲成人av在线免费| 一边亲一边摸免费视频| 欧美色欧美亚洲另类二区| 精品少妇黑人巨大在线播放 | 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 高清毛片免费看| 悠悠久久av| 偷拍熟女少妇极品色| 蜜桃亚洲精品一区二区三区| www.av在线官网国产| 亚洲美女视频黄频| 99riav亚洲国产免费| 韩国av在线不卡| h日本视频在线播放| 长腿黑丝高跟| 国产单亲对白刺激| 久久久久国产网址| 有码 亚洲区| 99久国产av精品| 白带黄色成豆腐渣| 又粗又硬又长又爽又黄的视频 | 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 亚洲精品粉嫩美女一区| 啦啦啦啦在线视频资源| h日本视频在线播放| 亚州av有码| 日韩一区二区三区影片| 夫妻性生交免费视频一级片| 国产不卡一卡二| 一进一出抽搐gif免费好疼| 国产精品一区二区三区四区久久| 国产精品久久视频播放| 亚洲欧美日韩高清在线视频| 精品人妻视频免费看| 春色校园在线视频观看| 免费人成在线观看视频色| 少妇人妻一区二区三区视频| 男人的好看免费观看在线视频| 亚洲成av人片在线播放无| 好男人在线观看高清免费视频| 亚洲av第一区精品v没综合| 老女人水多毛片| 人人妻人人澡欧美一区二区| av视频在线观看入口| 男女那种视频在线观看| 国产三级在线视频| 亚洲婷婷狠狠爱综合网| 两性午夜刺激爽爽歪歪视频在线观看| 69人妻影院| 国产亚洲精品av在线| 综合色av麻豆| 中国美白少妇内射xxxbb| 成年女人永久免费观看视频| 欧美高清成人免费视频www| 国产一区二区亚洲精品在线观看| 国产激情偷乱视频一区二区| 精品久久国产蜜桃| 亚洲av成人精品一区久久| 蜜臀久久99精品久久宅男| 18禁在线播放成人免费| 在线天堂最新版资源| 免费电影在线观看免费观看| 国产单亲对白刺激| 日本免费一区二区三区高清不卡| 高清毛片免费观看视频网站| 网址你懂的国产日韩在线| 最近中文字幕高清免费大全6| 亚洲成av人片在线播放无| 精品久久久久久久久av| 真实男女啪啪啪动态图| 久久婷婷人人爽人人干人人爱| 如何舔出高潮| 亚洲自拍偷在线| 秋霞在线观看毛片| 欧美一区二区精品小视频在线| 最近手机中文字幕大全| 亚洲经典国产精华液单| 成人午夜高清在线视频| 内射极品少妇av片p| 青春草亚洲视频在线观看| 成人欧美大片| 亚洲av免费在线观看| 99久久成人亚洲精品观看| 黄片wwwwww| 18禁裸乳无遮挡免费网站照片| 成人美女网站在线观看视频| 亚洲精品亚洲一区二区| 好男人视频免费观看在线| 啦啦啦韩国在线观看视频| 午夜a级毛片| 最近的中文字幕免费完整| 小说图片视频综合网站| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| 亚洲乱码一区二区免费版| 久久九九热精品免费| 欧美不卡视频在线免费观看| 成人特级av手机在线观看| 乱码一卡2卡4卡精品| av免费观看日本| 国产精品美女特级片免费视频播放器| 午夜福利在线在线| 国产v大片淫在线免费观看| 婷婷色综合大香蕉| 99久久成人亚洲精品观看| 国产精品,欧美在线| 亚洲av成人精品一区久久| 亚洲av成人av| 91精品一卡2卡3卡4卡| 黄色视频,在线免费观看| 特级一级黄色大片| 亚洲国产精品成人综合色| 免费看光身美女| 丝袜美腿在线中文| 国产精品久久久久久久电影| 国产精品一区www在线观看| 欧美+亚洲+日韩+国产| 亚洲人与动物交配视频| 国产真实伦视频高清在线观看| 国产精品伦人一区二区| 热99在线观看视频| 久久久久久久亚洲中文字幕| 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 综合色丁香网| 亚洲av二区三区四区| 最近手机中文字幕大全| 国产乱人视频| 蜜桃亚洲精品一区二区三区| 99热网站在线观看| av专区在线播放| 黄色视频,在线免费观看| 国产成人freesex在线| 男女视频在线观看网站免费| 国产成人精品久久久久久| 成人性生交大片免费视频hd| 黄片无遮挡物在线观看| 国内久久婷婷六月综合欲色啪| 爱豆传媒免费全集在线观看| 中文字幕av在线有码专区| 国产一区二区三区在线臀色熟女| 欧美激情国产日韩精品一区| 一区二区三区免费毛片| 亚洲欧美精品自产自拍| 中文字幕精品亚洲无线码一区| 蜜臀久久99精品久久宅男| 国产精品蜜桃在线观看 | 日本免费a在线| www.av在线官网国产| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看| 99久国产av精品国产电影| 在线天堂最新版资源| 亚洲欧美中文字幕日韩二区| 日韩欧美精品v在线| 高清毛片免费看| 亚洲自偷自拍三级| 欧美日本视频| 亚洲国产精品成人综合色| 你懂的网址亚洲精品在线观看 | 久久久欧美国产精品| 变态另类丝袜制服| 看黄色毛片网站| 我要搜黄色片| 美女脱内裤让男人舔精品视频 | 99视频精品全部免费 在线| 国产淫片久久久久久久久| 2022亚洲国产成人精品| 国产亚洲欧美98| 神马国产精品三级电影在线观看| 午夜精品在线福利| 亚洲经典国产精华液单| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 成人午夜高清在线视频| 尾随美女入室| 国产精品人妻久久久影院| 99热这里只有精品一区| 亚洲自拍偷在线| 久久久久性生活片| 亚洲激情五月婷婷啪啪| 1024手机看黄色片| 久久精品91蜜桃| 国产成人影院久久av| 国产黄片美女视频| 69av精品久久久久久| 久久鲁丝午夜福利片| 亚洲av第一区精品v没综合| av在线亚洲专区| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 爱豆传媒免费全集在线观看| 久久久a久久爽久久v久久| 99热网站在线观看| 久久精品国产鲁丝片午夜精品| 一进一出抽搐gif免费好疼| 91aial.com中文字幕在线观看| 在线免费十八禁| av.在线天堂| 亚洲精品久久国产高清桃花| 亚洲色图av天堂| 一区二区三区四区激情视频 | 日韩欧美三级三区| 观看免费一级毛片| 免费电影在线观看免费观看| 日韩欧美国产在线观看| 免费av毛片视频| 国产乱人偷精品视频| 亚洲美女视频黄频| 亚洲国产精品国产精品| www.色视频.com| 一本精品99久久精品77| 99久久中文字幕三级久久日本| 老女人水多毛片| 亚洲精品色激情综合| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久免费视频| 免费不卡的大黄色大毛片视频在线观看 | 中出人妻视频一区二区| 精品午夜福利在线看| 欧美激情国产日韩精品一区| 毛片女人毛片| 亚洲欧美精品专区久久| 乱码一卡2卡4卡精品| 嫩草影院精品99| 日本成人三级电影网站| 日本撒尿小便嘘嘘汇集6| 久久久成人免费电影| 美女被艹到高潮喷水动态| 99国产精品一区二区蜜桃av| 午夜福利高清视频| 国产91av在线免费观看| 亚洲欧洲日产国产| 一个人看的www免费观看视频| 成人性生交大片免费视频hd| 国内揄拍国产精品人妻在线| 麻豆一二三区av精品| 性插视频无遮挡在线免费观看| 日韩一本色道免费dvd| 国产精品久久电影中文字幕| 欧美一区二区国产精品久久精品| 国产精品.久久久| av女优亚洲男人天堂| 色噜噜av男人的天堂激情| 欧美最黄视频在线播放免费| 丝袜美腿在线中文| 欧美+亚洲+日韩+国产| 亚洲国产精品国产精品| 黄色视频,在线免费观看| 亚洲人成网站在线播放欧美日韩| 麻豆av噜噜一区二区三区| 黑人高潮一二区| 日本熟妇午夜| 午夜爱爱视频在线播放| 久久婷婷人人爽人人干人人爱| 久久久欧美国产精品| 最后的刺客免费高清国语| 国产三级在线视频| 国产成人精品一,二区 | 国产精品一区二区性色av| 欧美成人免费av一区二区三区| 日韩三级伦理在线观看| 久久久久国产网址| 国产麻豆成人av免费视频| www.色视频.com| 尾随美女入室| 国产精品嫩草影院av在线观看| 3wmmmm亚洲av在线观看| 亚洲丝袜综合中文字幕| 精品久久久久久久末码| 亚洲真实伦在线观看| 亚洲人成网站高清观看| 女人十人毛片免费观看3o分钟| 大型黄色视频在线免费观看| 午夜久久久久精精品| 偷拍熟女少妇极品色| 精品一区二区三区人妻视频| 午夜精品国产一区二区电影 | 日韩视频在线欧美| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区精品小视频在线| 色尼玛亚洲综合影院| 国产精品女同一区二区软件| 永久网站在线| 成熟少妇高潮喷水视频| 色综合站精品国产| 在线a可以看的网站| 桃色一区二区三区在线观看| 又粗又硬又长又爽又黄的视频 | 精品免费久久久久久久清纯| 亚洲成人久久性| 国产精品久久电影中文字幕| 最近视频中文字幕2019在线8| 亚洲中文字幕日韩| 国产精品久久久久久久久免| 亚洲精品粉嫩美女一区| 欧美日韩综合久久久久久| 久久精品国产99精品国产亚洲性色| 床上黄色一级片| 国产老妇女一区| 免费看日本二区| 一级毛片电影观看 | 欧美日本亚洲视频在线播放| 成人综合一区亚洲| 18禁在线无遮挡免费观看视频| 自拍偷自拍亚洲精品老妇| 中文字幕久久专区| 毛片一级片免费看久久久久| 亚洲精品久久久久久婷婷小说 | 一个人看视频在线观看www免费| 人妻少妇偷人精品九色| 亚洲av男天堂| 亚洲最大成人中文| 欧美xxxx性猛交bbbb| 深夜精品福利| 亚洲经典国产精华液单| 女人被狂操c到高潮| 精品久久国产蜜桃|