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

    基于SWAT 模型和Budyko 假設的石羊河流域徑流變化歸因?qū)Ρ确治?/h1>
    2023-08-27 15:41:53李瑤董增川張城邵逸卿倪效寬孟金玉徐亦凡魏園
    人民黃河 2023年7期
    關鍵詞:人類活動石羊河徑流

    李瑤 董增川 張城 邵逸卿 倪效寬 孟金玉 徐亦凡 魏園

    摘 要:探究流域徑流變化及其影響因素對流域水資源的可持續(xù)開發(fā)利用具有重要意義。利用1960—2019 年石羊河流域水文氣象資料,通過Mann-Kendall 突變檢驗和Pettitt 檢驗確定蔡旗站徑流序列的突變年份,采用SWAT 模型和基于Budyko 假設的彈性系數(shù)法進行徑流變化歸因?qū)Ρ确治?,定量辨識氣候變化和人類活動對徑流變化的影響。結(jié)果表明:近60 a 蔡旗站年徑流量呈現(xiàn)先減小后增大的趨勢,并在1972 年和2002 年發(fā)生突變;SWAT 模型與彈性系數(shù)法對徑流變化歸因分析結(jié)果基本一致,SWAT 模型分析結(jié)果表明,變化期1973—2002 年、2003—2019 年氣候變化對徑流變化的貢獻率分別為32.18%、46.09%,人類活動對徑流變化的貢獻率分別為-67.82%、-53.91%,而彈性系數(shù)法分析結(jié)果表明,兩個時期氣候變化對徑流變化的貢獻率分別為30.80%、41.44%,人類活動對徑流變化的貢獻率分別為-69.20%、-58.56%,相對于1960—1972 年,1973—2019 年石羊河流域徑流變化受人類活動的影響大于氣候變化的影響。人類活動是引起石羊河流域徑流變化的主要因素。

    關鍵詞:氣候變化;人類活動;徑流;SWAT 模型;Budyko;石羊河

    中圖分類號: TV121 文獻標志碼:A doi:10.3969/ j.issn.1000-1379.2023.07.006

    引用格式:李瑤,董增川,張城,等.基于SWAT 模型和Budyko 假設的石羊河流域徑流變化歸因?qū)Ρ确治觯郏剩荩嗣顸S河,2023,45(7):30-35.

    河川徑流不僅是水文循環(huán)的重要組成部分,而且是影響生態(tài)環(huán)境變化和社會經(jīng)濟發(fā)展的重要因素。在氣候變化與人類活動影響下,流域水循環(huán)過程和水資源時空分布規(guī)律發(fā)生了明顯改變,社會經(jīng)濟用水格局逐步演化,加劇了河流水文特性的復雜性[1] 。石羊河流域是西北干旱區(qū)重要的內(nèi)流河流域,《石羊河流域重點治理規(guī)劃》的實施,在緩解用水緊張、提高用水保證程度、改善流域生態(tài)環(huán)境方面取得了顯著成效,但受氣候變化和人類活動影響,仍面臨著水資源短缺、河湖生態(tài)功能萎縮、局部地下水位持續(xù)下降等問題。因此,研究變化環(huán)境下石羊河流域徑流演變規(guī)律及其驅(qū)動機制,對認識流域水文情勢以及探討流域水資源合理開發(fā)利用等具有重要意義。

    目前,徑流演變歸因分析的方法主要有數(shù)理統(tǒng)計法[2] 、水文模擬法[3] 和彈性系數(shù)法[4] 。其中:數(shù)理統(tǒng)計法主要利用回歸分析、神經(jīng)網(wǎng)絡等數(shù)學方法或統(tǒng)計模型量化主要驅(qū)動因子對徑流變化的貢獻率,該方法計算簡單,要求長序列的水文數(shù)據(jù);水文模擬法是基于流域水文模型或陸面過程模型定量分析各驅(qū)動因子對徑流變化的影響,該方法物理機制強,可在不同土地利用條件、多種氣候情景模式下進行模擬,分析精度較高;彈性系數(shù)法是依據(jù)Budyko 假設和彈性系數(shù)概念,計算氣候變化和人類活動作用下徑流的變化量,該方法物理意義明確,參數(shù)易獲取且計算簡單。

    近年來,許多學者對石羊河流域徑流變化進行了研究,研究多基于石羊河長期實測徑流資料,對徑流演變特征進行分析。在徑流變化歸因方面,薛東香[5] 采用改進的雙累積曲線法進行研究,發(fā)現(xiàn)石羊河上游產(chǎn)流區(qū)徑流變化主要受氣候變化影響,中下游徑流變化主要受人類活動影響;賈飛飛等[6] 采用累積量斜率變化率比較法進行研究,發(fā)現(xiàn)人類活動是黃羊河流域徑流變化的主導因素;陳忠升[7] 采用基于Budyko 假設的彈性系數(shù)法進行研究,發(fā)現(xiàn)石羊河中下游徑流減少的主要因素是人類活動。以往研究多從定性角度評價氣候變化和人類活動對徑流變化的影響,雖然也有一部分從定量角度進行分析,但研究方法主要集中于數(shù)理統(tǒng)計法和彈性系數(shù)法,運用基于物理機理的水文模擬法量化徑流變化貢獻率的研究還較少。因此,本文采用SWAT 模型和基于Budyko 假設的彈性系數(shù)法進行徑流變化歸因?qū)Ρ确治?,定量辨識氣候變化和人類活動對石羊河流域徑流變化的影響,以期為石羊河流域生態(tài)保護和治理規(guī)劃提供科學依據(jù)。

    1 研究區(qū)概況、數(shù)據(jù)和突變點分析

    1.1 研究區(qū)概況

    石羊河流域地處東經(jīng)101° 41′—104° 16′、北緯36°29′—39°27′之間,是甘肅省三大內(nèi)陸河流域之一,流域涉及武威、金昌、張掖和白銀四市的9 個縣(區(qū))。流域內(nèi)地勢南高北低,自西南向東北傾斜。石羊河發(fā)源于祁連山東沿線的冷龍嶺大雪山,河長300 km,流域面積4.16 萬km2,其中產(chǎn)流面積1.11 萬km2,多年平均徑流量15.60 億m3。石羊河水系由西大河、東大河、西營河、金塔河、雜木河、黃羊河、古浪河與大靖河等8 條支流和眾多小溝河組成。

    1.2 數(shù)據(jù)來源

    本研究所用數(shù)據(jù)有數(shù)字高程模型(DEM)數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)、氣象數(shù)據(jù)和水文數(shù)據(jù)。其中:DEM 數(shù)據(jù)為SRTMGL1.003 30 m 分辨率數(shù)字高程數(shù)據(jù),來源于美國地質(zhì)調(diào)查局(USGS);土地利用數(shù)據(jù)為1980 年甘肅省土地利用遙感監(jiān)測數(shù)據(jù),分辨率為1 km,來源于中國科學院資源環(huán)境科學與數(shù)據(jù)中心;土壤數(shù)據(jù)采用世界土壤數(shù)據(jù)庫(HWSD V1.1),分辨率為1 km,來源于國家冰川凍土沙漠科學數(shù)據(jù)中心;氣象數(shù)據(jù)為研究區(qū)及周邊地區(qū)11 個氣象站1960—2019年逐日降水、最高氣溫、最低氣溫、平均氣溫、日照時數(shù)、相對濕度和平均風速等,來源于中國氣象數(shù)據(jù)網(wǎng);水文數(shù)據(jù)為蔡旗水文站1960—2019 年日徑流數(shù)據(jù),來源于內(nèi)陸河水文年鑒。水文站與氣象站分布如圖1所示。

    1.3 突變點識別

    蔡旗站年徑流量變化過程如圖2 所示。由圖2 可知,年徑流量在2002 年前后發(fā)生明顯變化,總體呈先減小后增大趨勢,多年平均徑流量為2.73 億m3,最大值出現(xiàn)在1967 年,為5.96 億m3,最小值出現(xiàn)在2002 年,為0.85 億m3。由5 a 滑動平均曲線可知,蔡旗站徑流量呈減少(1964—1966 年)—增加(1967—1971 年)—減少(1972—2002 年)—增加(2003—2019 年)的變化趨勢。通過時間序列離散度計算,蔡旗站年徑流量的變差系數(shù)CV值為0.47,徑流年際間波動較為劇烈。

    采用Mann-Kendall 檢驗(M-K 檢驗)[8] 和Pettitt突變檢驗[9] 結(jié)合的方法對徑流序列進行突變性分析,確定突變點。徑流突變檢驗結(jié)果如圖3 所示,蔡旗站M- K 統(tǒng)計量序列曲線交點發(fā)生在1972 年, 另有1974—2019 年UF 曲線超過置信區(qū)間,是可能發(fā)生突變的區(qū)域;經(jīng)過Pettitt 突變點識別,突變點可能發(fā)生的年份為1972 年,統(tǒng)計量為558,置信水平p = 1.356×10-5<0.05,通過了置信水平為0.05 的顯著性檢驗。結(jié)合圖2 可知2002 年是較明確的突變點,且在M-K 檢驗可能發(fā)生突變的區(qū)域,因此可以認為蔡旗站年徑流量的突變點為1972 年和2002 年。

    2 徑流變化歸因分析

    2.1 SWAT 模型歸因分析

    1)SWAT 建模與評價。SWAT 模型是美國農(nóng)業(yè)部農(nóng)業(yè)研究中心(USDA-ARS)開發(fā)的基于長時段、適用于較大流域尺度的分布式水文模型[10] ,可以同時考慮地形、土地利用、土壤、氣象等對水文過程的影響,被廣泛應用于流域水文過程模擬以及流域在氣候和土地利用變化下的水文響應等。

    基于石羊河流域DEM 數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)、氣象數(shù)據(jù)等,通過土壤水文特性軟件、SWAT?Weather 等構建模型的土地利用、土壤、氣象數(shù)據(jù)庫,劃分子流域和水文響應單元,建立流域SWAT 水文模型,最終將流域劃分為47 個子流域和327 個水文響應單元。建模結(jié)果如圖4 所示。

    采用SWAT 模型對流域1960—1972 年徑流量進行月尺度模擬。設置1960 年為模型預熱期,1961—1967 年為率定期,1968—1972 年為驗證期。采用SWAT-CUP 軟件的SUFI2 算法對模型進行率定和驗證,通過t-Stat 和P-Value 兩個指標判別參數(shù)的敏感性,敏感參數(shù)及率定結(jié)果見表1。徑流主要受河道深層地下水損失率、地下水蒸發(fā)系數(shù)、地下水滯后系數(shù)、降雪氣溫、側(cè)向流移動時間、土壤有效含水率等參數(shù)影響。

    采用決定系數(shù)(R2 ) 和Nash-Sutcliffe 效率系數(shù)(NSE)對計算結(jié)果進行精度分析,綜合評價SWAT 模型的模擬效果,一般認為,當R2 與NSE 均大于0.5 時模型模擬結(jié)果可接受[11] 。由表2 可知,率定期和驗證期的R2和NSE 均大于0.8,表明SWAT 模型能夠較好地模擬石羊河流域的徑流過程。

    率定期和驗證期的月均流量模擬值與實測值對比及其相關關系如圖5 和圖6 所示,可以直觀看出SWAT 模型在率定期和驗證期月徑流模擬結(jié)果與實測徑流過程吻合較好。

    2)貢獻率分析。運用SWAT 模型[12] 定量評估氣候變化和人類活動對徑流變化的貢獻率。具體步驟:利用基準期(1960—1972 年)的氣象、水文和土地利用等數(shù)據(jù)率定水文模型參數(shù),將變化期(1973—2002 年、2003—2019 年)的氣象數(shù)據(jù)輸入模型中還原出天然徑流量Rsim,并計算氣候變化導致的徑流變化量ΔRC 和人類活動導致的徑流變化量ΔRω,利用突變前后徑流總變化量ΔR,分別計算二者的貢獻率ηC、ηω。具體計算公式如下:

    以徑流量突變前的1960—1972 年為基準期,突變后的1973—2002 年、2003—2019 年為變化期。由表4可知,變化期1973—2002 年、2003—2019 年年均降水量相對基準期分別增加55.21、95.03 mm,年均徑流深相對基準期分別減少23.05、18.27 mm,年均潛在蒸散發(fā)量相對基準期分別增加9.17、55.22 mm。變化期的徑流系數(shù)(R / P)、干旱指數(shù)(ET0 / P)均小于基準期的,而下墊面參數(shù)絕對值較基準期明顯增大。

    根據(jù)式(12) ~式(14)計算得到徑流的降水量、潛在蒸散發(fā)量、下墊面參數(shù)彈性系數(shù),徑流對降水量的彈性系數(shù)為正值,而對潛在蒸散發(fā)量和下墊面參數(shù)的彈性系數(shù)為負值,表示降水量增加將導致徑流量增大,潛在蒸散發(fā)量、下墊面參數(shù)的增大則會導致徑流量減小。變化期相對基準期,各項關于徑流變化彈性系數(shù)絕對值均增大,徑流變化對各因子敏感性:下墊面參數(shù)>降水量>潛在蒸散發(fā)量。

    由式(16)計算徑流變化貢獻率,結(jié)果見表5。降水量和潛在蒸散發(fā)量對徑流變化的貢獻率之和是氣候變化對徑流變化的貢獻率,人類活動對徑流變化的貢獻率主要由下墊面參數(shù)來反映。降水對徑流變化的貢獻率為正值,潛在蒸散發(fā)量、下墊面參數(shù)的貢獻率為負值。下墊面參數(shù)對徑流影響最大。相對于1960—1972 年,1973—2019 年人類活動對石羊河流域徑流變化影響大于氣候變化的影響,其中1973—2002 年、2003—2019 年2 個時期氣候變化對徑流變化的貢獻率分別為30.80%、41.44%,人類活動對徑流變化的貢獻率分別為-69.20%、-58.56%。

    2.3 SWAT 模型與Budyko 假設歸因?qū)Ρ确治?/p>

    SWAT 模型徑流變化歸因分析結(jié)果表明,變化期1973—2002 年、2003—2019 年氣候變化對徑流變化的貢獻率分別為32.18%、46.09%,人類活動對徑流變化的貢獻率分別為-67.82%、-53.91%,而彈性系數(shù)法結(jié)果表明,兩個時期氣候變化對徑流變化的貢獻率分別為30.80%、41.44%,人類活動對徑流變化的貢獻率分別為-69.20%、-58.56%。對比發(fā)現(xiàn),兩種方法所得徑流變化歸因分析結(jié)果基本一致,相對于1960—1972年,1973—2019 年石羊河流域徑流變化受人類活動的影響大于氣候變化的影響。與彈性系數(shù)法相比,SWAT 模型物理機制強,分析精度更高,能夠?qū)λ淖兞考捌溲葑円?guī)律進行準確的描述,可進一步采用SWAT 模型對石羊河流域未來徑流進行預測。

    3 結(jié)論

    1)石羊河流域年徑流量呈現(xiàn)先減小后增大的趨勢,且在1972 年和2002 年發(fā)生突變。

    2)SWAT 模型與基于Budyko 假設的彈性系數(shù)法徑流變化歸因分析結(jié)果基本一致,SWAT 模型分析結(jié)果表明,變化期1973—2002 年、2003—2019 年氣候變化對徑流變化的貢獻率分別為32.18%、46.09%,人類活動對徑流變化的貢獻率分別為-67.82%、-53.91%,而彈性系數(shù)法結(jié)果表明,兩個時期氣候變化對徑流變化的貢獻率分別為30.80%、41.44%,人類活動對徑流變化的貢獻率分別為-69.20%、-58.56%。相對于1960—1972 年,1973—2019 年石羊河流域徑流變化受人類活動的影響大于氣候變化的影響,人類活動是影響石羊河流域徑流變化的主要因素。

    3)人類活動是導致石羊河流域徑流減少的主要原因,科學合理地規(guī)劃、開發(fā)、利用水資源對緩解石羊河流域水資源短缺、改善流域生態(tài)環(huán)境具有重要意義。

    參考文獻:

    [1] ZHANG H L,MENG C C,WANG Y J,et al.ComprehensiveEvaluation of the Effects of Climate Change and Land Useand Land Cover Change Variables on Runoff and SedimentDischarge [ J]. Science of the Total Environment, 2020,702:134401.

    [2] OCHOA?TOCACHI B F,BUYTAERT W,BIEVRE B D.Re?gionalization of Land?Use Impacts on Streamflow Using a Net?work of Paired Catchments[J].Water Resources Research,2016,52(9):6710-6729.

    [3] 徐宗學,姜瑤.變化環(huán)境下的徑流演變與影響研究:回顧與展望[J].水利水運工程學報,2022(1):9-18.

    [4] 慕星,趙勇,劉歡,等.氣候變化和人類活動影響下徑流演變研究進展[J].人民黃河,2021,43(5):35-41.

    [5] 薛東香.石羊河流域徑流變化及歸因分析[D].蘭州:西北師范大學,2021:53-67.

    [6] 賈飛飛,李鑫,孫翠洋,等.1951—2016 年黃羊河流域徑流變化特征及影響因素分析[J].遼寧師范大學學報(自然科學版),2019,42(2):262-267.

    [7] 陳忠升.中國西北干旱區(qū)河川徑流變化及歸因定量辨識[D].上海:華東師范大學,2016:161-202.

    [8] MANN H B.Nonparametric Tests Against Trend[J].Econometri?ca:Journal of the Econometric Society,1945,13:245-259.

    [9] 鄒磊,夏軍,張印,等.海河流域降水時空演變特征及其驅(qū)動力分析[J].水資源保護,2021,37(1):53-60.

    [10] 芮孝芳.水文學研究進展[M].南京:河海大學出版社,2007:10-15.

    [11] MORIASI D,ARNOLD J,LIEW M W V,et al.Model Evalu?ation Guidelines for Systematic Quantification of Accuracyin Watershed Simulations[J].Transactions of the ASABE,2007,50(3):885-900.

    [12] DOSDOGRU F,KALIN L,WANG R Y,et al.Potential Im?pacts of Land Use/ Cover and Climate Changes on Ecologi?cally Relevant Flows [ J]. Journal of Hydrology, 2020,584:124654.

    [13] 劉昌明,魏豪杉,張永強,等.黃河干流徑流變化歸因分析與有關問題商榷[J].人民黃河,2021,43(10):1-6,16.

    [14] 傅抱璞.論陸面蒸發(fā)的計算[J].大氣科學,1981,5(1):23-31.

    [15] XU X Y,YANG D W,YANG H B,et al.Attribution AnalysisBased on the Budyko Hypothesis for Detecting the DominantCause of Runoff Decline in Haihe Basin[J].Journal of Hy?drology,2014,510:530-540.

    【責任編輯 張 帥】

    猜你喜歡
    人類活動石羊河徑流
    基于不同旱情指數(shù)的石羊河流域春旱監(jiān)測研究
    石羊河流域永昌縣地下水及水資源供需平衡分析
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    氣候變化和人類活動對漢江上游徑流變化影響的定量研究
    Topmodel在布哈河流域徑流模擬中的應用
    加強投入 推進石羊河流域源頭治理
    石羊河中下游流域水環(huán)境化學特征研究
    地下水(2015年1期)2015-12-14 09:22:44
    人為因素主導下海岸帶生態(tài)系統(tǒng)變遷探討
    綠色科技(2015年6期)2015-08-05 18:01:22
    山地表生地球化學研究現(xiàn)狀與展望
    探秘“大徑流”

    一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片| 在线天堂最新版资源| 亚洲精品日韩av片在线观看| 日本免费a在线| 欧美色视频一区免费| 国产精品永久免费网站| 色综合色国产| 国产中年淑女户外野战色| 婷婷色av中文字幕| 免费av观看视频| 我的老师免费观看完整版| 九九久久精品国产亚洲av麻豆| 美女xxoo啪啪120秒动态图| 国内精品宾馆在线| 99热这里只有是精品50| 国产精品一区二区性色av| 国产精品av视频在线免费观看| 少妇裸体淫交视频免费看高清| 国产黄色视频一区二区在线观看 | a级一级毛片免费在线观看| 欧美日韩一区二区视频在线观看视频在线 | 床上黄色一级片| 亚洲精品影视一区二区三区av| 中出人妻视频一区二区| 久久久久久久久久成人| 亚洲精品久久久久久婷婷小说 | 一进一出抽搐gif免费好疼| 免费观看a级毛片全部| 九九久久精品国产亚洲av麻豆| av视频在线观看入口| 成人二区视频| 国产精品一区二区三区四区免费观看| 欧美日韩精品成人综合77777| 在线观看美女被高潮喷水网站| 久久人人爽人人片av| 国产伦精品一区二区三区四那| 成人亚洲欧美一区二区av| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 久久草成人影院| 亚洲性久久影院| av天堂在线播放| 成人特级av手机在线观看| 亚洲欧美中文字幕日韩二区| 看片在线看免费视频| 国产高潮美女av| 精品午夜福利在线看| 最近手机中文字幕大全| 99精品在免费线老司机午夜| 国产精品一区二区性色av| 国产亚洲91精品色在线| 狂野欧美白嫩少妇大欣赏| 成人一区二区视频在线观看| 欧美+亚洲+日韩+国产| 国产一区二区在线av高清观看| 国产片特级美女逼逼视频| 成人美女网站在线观看视频| 男人舔奶头视频| 少妇的逼好多水| 久久久久网色| 国产午夜精品一二区理论片| 国产三级在线视频| 日日干狠狠操夜夜爽| 嘟嘟电影网在线观看| 内射极品少妇av片p| 国产精品无大码| 青青草视频在线视频观看| 久久精品久久久久久久性| 久久久久免费精品人妻一区二区| 亚洲真实伦在线观看| 最近视频中文字幕2019在线8| 午夜激情欧美在线| 亚洲性久久影院| 91av网一区二区| av视频在线观看入口| 蜜桃久久精品国产亚洲av| 国产毛片a区久久久久| 国产一区二区在线观看日韩| 最后的刺客免费高清国语| 中文在线观看免费www的网站| 小说图片视频综合网站| 亚洲内射少妇av| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线观看播放| 久久99蜜桃精品久久| 欧美性感艳星| 性欧美人与动物交配| 内射极品少妇av片p| 色噜噜av男人的天堂激情| 色播亚洲综合网| 亚洲图色成人| 永久网站在线| a级毛片a级免费在线| 久久人人爽人人片av| 亚洲性久久影院| 精品久久久久久久久亚洲| 免费搜索国产男女视频| 国产精品麻豆人妻色哟哟久久 | 99久久九九国产精品国产免费| 亚洲欧美成人综合另类久久久 | 一本精品99久久精品77| 91久久精品电影网| 激情 狠狠 欧美| 免费看日本二区| 九草在线视频观看| 国产私拍福利视频在线观看| 午夜爱爱视频在线播放| 成人永久免费在线观看视频| 18禁黄网站禁片免费观看直播| 欧美激情在线99| 日本三级黄在线观看| 国产熟女欧美一区二区| 97在线视频观看| 免费在线观看成人毛片| 岛国毛片在线播放| 精品午夜福利在线看| 国产精品美女特级片免费视频播放器| 国产午夜精品论理片| 欧美又色又爽又黄视频| 亚洲国产欧美在线一区| 国产亚洲91精品色在线| 青春草亚洲视频在线观看| 国国产精品蜜臀av免费| 中文字幕熟女人妻在线| 日韩高清综合在线| 晚上一个人看的免费电影| 亚洲一区高清亚洲精品| 99热这里只有是精品50| 成人欧美大片| 国产精品1区2区在线观看.| 国产一级毛片七仙女欲春2| 日韩 亚洲 欧美在线| 又爽又黄a免费视频| 日韩大尺度精品在线看网址| 欧美人与善性xxx| 日本黄大片高清| 国产精品福利在线免费观看| 欧美+亚洲+日韩+国产| 黄片无遮挡物在线观看| 亚州av有码| 丰满乱子伦码专区| 好男人视频免费观看在线| 一级二级三级毛片免费看| 男女边吃奶边做爰视频| 亚洲av免费在线观看| 美女国产视频在线观看| 国产一区二区激情短视频| 男女视频在线观看网站免费| 能在线免费观看的黄片| 女人被狂操c到高潮| 天堂影院成人在线观看| 国产精品久久久久久久久免| 成年免费大片在线观看| 国产片特级美女逼逼视频| 国产片特级美女逼逼视频| 麻豆精品久久久久久蜜桃| 男人狂女人下面高潮的视频| 国产又黄又爽又无遮挡在线| 国产色爽女视频免费观看| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| 寂寞人妻少妇视频99o| 亚洲国产高清在线一区二区三| 成人三级黄色视频| 亚洲美女视频黄频| 九九热线精品视视频播放| 在线观看午夜福利视频| 嘟嘟电影网在线观看| 97超碰精品成人国产| 国产熟女欧美一区二区| 男女下面进入的视频免费午夜| 国产成人福利小说| 卡戴珊不雅视频在线播放| 亚洲av不卡在线观看| 亚洲,欧美,日韩| 国产成人精品久久久久久| 全区人妻精品视频| 小说图片视频综合网站| 你懂的网址亚洲精品在线观看 | 可以在线观看毛片的网站| 人人妻人人看人人澡| 欧美日本视频| 变态另类成人亚洲欧美熟女| 国产亚洲5aaaaa淫片| 亚洲自拍偷在线| 国产一区二区在线观看日韩| 亚洲精品成人久久久久久| 日韩欧美一区二区三区在线观看| 日产精品乱码卡一卡2卡三| 99国产极品粉嫩在线观看| 少妇人妻一区二区三区视频| 亚洲真实伦在线观看| 日韩制服骚丝袜av| 国产精品免费一区二区三区在线| 亚洲一级一片aⅴ在线观看| 日日撸夜夜添| 五月玫瑰六月丁香| 一个人免费在线观看电影| 亚洲精品影视一区二区三区av| 欧美色欧美亚洲另类二区| 国产综合懂色| .国产精品久久| 午夜a级毛片| 好男人视频免费观看在线| 高清毛片免费看| 女同久久另类99精品国产91| www.色视频.com| 久久6这里有精品| 极品教师在线视频| 一个人看视频在线观看www免费| 亚洲真实伦在线观看| 免费看美女性在线毛片视频| 成年免费大片在线观看| 伦精品一区二区三区| 亚洲乱码一区二区免费版| 亚洲最大成人av| 麻豆精品久久久久久蜜桃| 亚洲欧美日韩东京热| 国产一区二区在线观看日韩| 男女视频在线观看网站免费| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 国产高清有码在线观看视频| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 久久精品夜色国产| 国产成人精品久久久久久| 亚洲七黄色美女视频| 亚洲丝袜综合中文字幕| 久久久久久久久久黄片| 亚洲久久久久久中文字幕| 国产爱豆传媒在线观看| 97热精品久久久久久| a级一级毛片免费在线观看| 久久99热这里只有精品18| 国产中年淑女户外野战色| 国产老妇女一区| 我要看日韩黄色一级片| 天堂网av新在线| 自拍偷自拍亚洲精品老妇| 久久久久久大精品| 嫩草影院入口| 中文欧美无线码| 国产私拍福利视频在线观看| 精品免费久久久久久久清纯| 久久久久性生活片| 色5月婷婷丁香| 3wmmmm亚洲av在线观看| 性插视频无遮挡在线免费观看| 亚洲美女搞黄在线观看| 国产日韩欧美在线精品| 精品熟女少妇av免费看| 亚洲国产欧洲综合997久久,| 亚洲国产精品成人久久小说 | 色综合亚洲欧美另类图片| 久久久色成人| 久久久欧美国产精品| 免费看a级黄色片| 爱豆传媒免费全集在线观看| 成人特级黄色片久久久久久久| 人人妻人人澡欧美一区二区| 精品一区二区免费观看| a级毛片a级免费在线| 亚洲aⅴ乱码一区二区在线播放| 在线观看午夜福利视频| 乱码一卡2卡4卡精品| 国产伦一二天堂av在线观看| 亚洲欧美精品专区久久| 少妇裸体淫交视频免费看高清| 狂野欧美白嫩少妇大欣赏| 观看免费一级毛片| 在线天堂最新版资源| 色吧在线观看| 性插视频无遮挡在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一区二区三区四区免费观看| 国产蜜桃级精品一区二区三区| 一本久久中文字幕| 99精品在免费线老司机午夜| 自拍偷自拍亚洲精品老妇| 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产亚洲av涩爱 | 国产成年人精品一区二区| 亚洲av免费高清在线观看| 亚洲av电影不卡..在线观看| 国产精品乱码一区二三区的特点| 国内精品美女久久久久久| 国产蜜桃级精品一区二区三区| 成人无遮挡网站| 亚洲三级黄色毛片| 神马国产精品三级电影在线观看| a级毛片免费高清观看在线播放| 欧美xxxx性猛交bbbb| av黄色大香蕉| 日本-黄色视频高清免费观看| 婷婷色综合大香蕉| 欧美一区二区国产精品久久精品| 国产老妇伦熟女老妇高清| 成人欧美大片| av在线蜜桃| 国产精品免费一区二区三区在线| 国内精品宾馆在线| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 毛片女人毛片| 人人妻人人看人人澡| 国产一区二区在线观看日韩| 禁无遮挡网站| 国产在线男女| 久久亚洲精品不卡| 国产亚洲精品久久久com| 国产成人午夜福利电影在线观看| 美女大奶头视频| 最近手机中文字幕大全| 国产伦在线观看视频一区| 少妇丰满av| av女优亚洲男人天堂| 亚洲欧美日韩东京热| 国产精品一区二区三区四区久久| 啦啦啦啦在线视频资源| 国产 一区 欧美 日韩| 国内久久婷婷六月综合欲色啪| 两个人的视频大全免费| 久久韩国三级中文字幕| 看黄色毛片网站| 国产精品一及| 一级av片app| 熟女电影av网| 欧美一区二区亚洲| 99久久成人亚洲精品观看| 日韩一本色道免费dvd| 日日摸夜夜添夜夜爱| 波多野结衣巨乳人妻| avwww免费| 欧美色视频一区免费| 国产成人freesex在线| 狂野欧美激情性xxxx在线观看| 美女被艹到高潮喷水动态| 国产男人的电影天堂91| 国产日韩欧美在线精品| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 一区福利在线观看| 在线播放国产精品三级| 免费av毛片视频| 日本在线视频免费播放| 国产私拍福利视频在线观看| 国产在线精品亚洲第一网站| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 亚洲18禁久久av| 久久人妻av系列| 美女黄网站色视频| 国产色婷婷99| 成年女人永久免费观看视频| 少妇的逼水好多| 色综合站精品国产| 久久久久久久久中文| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| 成年女人永久免费观看视频| 亚洲欧美成人精品一区二区| 青青草视频在线视频观看| 国产午夜精品论理片| 国产精品女同一区二区软件| 日韩,欧美,国产一区二区三区 | 好男人视频免费观看在线| 国产精品一区二区在线观看99 | 91麻豆精品激情在线观看国产| 久久久久久久久久久丰满| 久久久成人免费电影| 色尼玛亚洲综合影院| 精品一区二区免费观看| av在线老鸭窝| av又黄又爽大尺度在线免费看 | 久久久久久久久久成人| 成人午夜高清在线视频| 亚洲av二区三区四区| 搡老妇女老女人老熟妇| av专区在线播放| 久久婷婷人人爽人人干人人爱| 日本与韩国留学比较| 长腿黑丝高跟| 日本欧美国产在线视频| 狂野欧美白嫩少妇大欣赏| 成年版毛片免费区| 九色成人免费人妻av| 婷婷色av中文字幕| 国内精品美女久久久久久| 69av精品久久久久久| 国产激情偷乱视频一区二区| 日韩制服骚丝袜av| 国产成人aa在线观看| 波野结衣二区三区在线| 欧美性猛交╳xxx乱大交人| 久久久精品大字幕| 国产高清激情床上av| 麻豆乱淫一区二区| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 久久久欧美国产精品| 99久久中文字幕三级久久日本| 午夜福利成人在线免费观看| 精品一区二区免费观看| 亚洲在线观看片| 国产av不卡久久| 国产午夜精品一二区理论片| 校园人妻丝袜中文字幕| 三级男女做爰猛烈吃奶摸视频| 高清毛片免费观看视频网站| av在线蜜桃| 欧美xxxx性猛交bbbb| 2021天堂中文幕一二区在线观| 免费看av在线观看网站| 久久婷婷人人爽人人干人人爱| 国产又黄又爽又无遮挡在线| 男女边吃奶边做爰视频| 亚洲,欧美,日韩| 69av精品久久久久久| 在线观看66精品国产| 一本久久精品| 日韩av不卡免费在线播放| 国产视频内射| 黄片wwwwww| 在线播放国产精品三级| 一边亲一边摸免费视频| 成人永久免费在线观看视频| 亚洲高清免费不卡视频| 1000部很黄的大片| 成年免费大片在线观看| 在线观看一区二区三区| 久久久久久久亚洲中文字幕| 狂野欧美激情性xxxx在线观看| 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 女同久久另类99精品国产91| 久99久视频精品免费| 波多野结衣高清无吗| 国产一级毛片在线| 国产私拍福利视频在线观看| 成人午夜精彩视频在线观看| 69人妻影院| 一个人看的www免费观看视频| 日本一二三区视频观看| 亚洲精品自拍成人| 你懂的网址亚洲精品在线观看 | 婷婷色av中文字幕| www.色视频.com| 成人综合一区亚洲| 国产精品女同一区二区软件| 亚洲av二区三区四区| 给我免费播放毛片高清在线观看| 久久久午夜欧美精品| 精品日产1卡2卡| 我的女老师完整版在线观看| 天堂av国产一区二区熟女人妻| 欧美性猛交╳xxx乱大交人| 少妇高潮的动态图| 亚洲最大成人手机在线| 国产精品久久久久久亚洲av鲁大| 欧美丝袜亚洲另类| 国产一区二区在线观看日韩| 亚洲人成网站高清观看| 欧美性猛交╳xxx乱大交人| 大香蕉久久网| 欧美色视频一区免费| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 深夜精品福利| 91狼人影院| 亚洲国产精品成人综合色| 国产极品精品免费视频能看的| 乱系列少妇在线播放| 村上凉子中文字幕在线| 一本精品99久久精品77| 日韩强制内射视频| 熟妇人妻久久中文字幕3abv| 亚洲av一区综合| 久久久久久久亚洲中文字幕| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 欧美成人精品欧美一级黄| 亚洲一区二区三区色噜噜| av视频在线观看入口| 尾随美女入室| 亚洲内射少妇av| 久久精品国产鲁丝片午夜精品| 亚洲真实伦在线观看| 精品人妻熟女av久视频| 日本-黄色视频高清免费观看| 丝袜美腿在线中文| 97超碰精品成人国产| 国产精品免费一区二区三区在线| 午夜福利成人在线免费观看| 2021天堂中文幕一二区在线观| 1024手机看黄色片| 女同久久另类99精品国产91| 亚洲av熟女| 不卡视频在线观看欧美| or卡值多少钱| 黄色欧美视频在线观看| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 亚洲精品日韩在线中文字幕 | 尾随美女入室| 欧美一区二区国产精品久久精品| 国产成人精品一,二区 | 99久久人妻综合| 3wmmmm亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片| 天天躁夜夜躁狠狠久久av| 亚洲欧美中文字幕日韩二区| 成人毛片a级毛片在线播放| 一本一本综合久久| 你懂的网址亚洲精品在线观看 | 国产成人aa在线观看| 精品人妻偷拍中文字幕| 亚洲国产色片| 国产免费一级a男人的天堂| 在现免费观看毛片| ponron亚洲| 最近的中文字幕免费完整| 国产av在哪里看| 久久精品夜色国产| 又粗又爽又猛毛片免费看| 久久久国产成人免费| 国产精品一区二区性色av| 春色校园在线视频观看| 欧美激情久久久久久爽电影| 丰满乱子伦码专区| 久久这里有精品视频免费| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| 菩萨蛮人人尽说江南好唐韦庄 | 国产黄a三级三级三级人| 日韩大尺度精品在线看网址| 国产精品一区二区三区四区久久| 国产日韩欧美在线精品| 亚州av有码| 成人欧美大片| av又黄又爽大尺度在线免费看 | 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放 | 国产亚洲91精品色在线| 久久人妻av系列| 亚洲综合色惰| 色尼玛亚洲综合影院| 精品久久久久久久久亚洲| 性欧美人与动物交配| 国产三级在线视频| 日韩一区二区三区影片| 久久人人精品亚洲av| 国产日本99.免费观看| 国产亚洲5aaaaa淫片| 99热6这里只有精品| 国产三级中文精品| 美女 人体艺术 gogo| 国产精品久久久久久亚洲av鲁大| 欧美色欧美亚洲另类二区| 日产精品乱码卡一卡2卡三| а√天堂www在线а√下载| 久久久久免费精品人妻一区二区| 日本免费a在线| 一个人看的www免费观看视频| 非洲黑人性xxxx精品又粗又长| 九九久久精品国产亚洲av麻豆| 亚洲欧美成人精品一区二区| 午夜福利视频1000在线观看| 亚洲成人中文字幕在线播放| 国产欧美日韩精品一区二区| 亚洲高清免费不卡视频| 人人妻人人澡欧美一区二区| 亚洲人成网站高清观看| 亚洲av熟女| 日韩精品青青久久久久久| 午夜久久久久精精品| 欧美日韩乱码在线| 成人毛片a级毛片在线播放| 日本三级黄在线观看| av在线天堂中文字幕| 寂寞人妻少妇视频99o| 亚洲人成网站在线播| 激情 狠狠 欧美| 永久网站在线| 人妻久久中文字幕网| 亚洲精品国产成人久久av| 高清午夜精品一区二区三区 | 久久婷婷人人爽人人干人人爱| 欧美日本视频| 在现免费观看毛片| 国内揄拍国产精品人妻在线| 日韩成人伦理影院| 国产片特级美女逼逼视频| 国产蜜桃级精品一区二区三区| 国产伦一二天堂av在线观看| 成人av在线播放网站| 99热网站在线观看| 亚洲欧美日韩东京热| 美女 人体艺术 gogo| 国产一区亚洲一区在线观看| 看十八女毛片水多多多| 国产精品一二三区在线看| 国产国拍精品亚洲av在线观看| 看十八女毛片水多多多| 国产69精品久久久久777片| 国产成人福利小说| 亚洲人成网站在线播放欧美日韩| 晚上一个人看的免费电影| 欧美一区二区国产精品久久精品| 你懂的网址亚洲精品在线观看 | 黄色日韩在线|