• <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)狀與展望
    探秘“大徑流”

    精品久久久久久久久久久久久| 一个人看的www免费观看视频| 看免费成人av毛片| 男女边摸边吃奶| 丰满人妻一区二区三区视频av| 视频中文字幕在线观看| 91精品国产九色| 大又大粗又爽又黄少妇毛片口| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 卡戴珊不雅视频在线播放| 最近2019中文字幕mv第一页| 看黄色毛片网站| 精品一区二区三区人妻视频| 国产 亚洲一区二区三区 | 国产高清国产精品国产三级 | 亚洲欧美成人综合另类久久久| 舔av片在线| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 国产黄色免费在线视频| 日韩强制内射视频| 国产中年淑女户外野战色| 成人综合一区亚洲| 国产一级毛片七仙女欲春2| 国产成人免费观看mmmm| 日本一二三区视频观看| 男女那种视频在线观看| 免费无遮挡裸体视频| 日韩电影二区| 毛片一级片免费看久久久久| 国产亚洲最大av| 久久久精品94久久精品| 午夜老司机福利剧场| 久久久久久久国产电影| 亚洲成色77777| 少妇被粗大猛烈的视频| 免费观看无遮挡的男女| 亚洲精品中文字幕在线视频 | 91aial.com中文字幕在线观看| 天堂√8在线中文| 国产亚洲精品久久久com| 直男gayav资源| 欧美最新免费一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲va在线va天堂va国产| 大陆偷拍与自拍| 日本色播在线视频| 免费大片18禁| 亚洲国产成人一精品久久久| 免费av不卡在线播放| 亚洲欧洲国产日韩| 一本久久精品| 夜夜爽夜夜爽视频| 亚洲精品视频女| 激情五月婷婷亚洲| 色视频www国产| 国产伦理片在线播放av一区| 日本午夜av视频| 啦啦啦中文免费视频观看日本| 男人狂女人下面高潮的视频| 国产精品一区二区在线观看99 | 国产淫片久久久久久久久| 精品一区二区三卡| 神马国产精品三级电影在线观看| 超碰av人人做人人爽久久| 日韩 亚洲 欧美在线| 午夜精品国产一区二区电影 | 国产片特级美女逼逼视频| 国产 一区 欧美 日韩| 精品国产露脸久久av麻豆 | 超碰97精品在线观看| 国产午夜精品一二区理论片| 亚洲美女搞黄在线观看| 国产老妇女一区| 日韩欧美 国产精品| 久久久久久久久久黄片| 午夜精品在线福利| 亚洲人成网站在线观看播放| 日韩精品青青久久久久久| 狂野欧美激情性xxxx在线观看| 亚洲成人久久爱视频| a级一级毛片免费在线观看| 日韩人妻高清精品专区| 一级二级三级毛片免费看| 天堂俺去俺来也www色官网 | 午夜亚洲福利在线播放| 亚洲精品一区蜜桃| 边亲边吃奶的免费视频| 激情 狠狠 欧美| 天堂俺去俺来也www色官网 | av在线播放精品| 午夜福利成人在线免费观看| 高清午夜精品一区二区三区| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 午夜视频国产福利| av免费观看日本| 毛片一级片免费看久久久久| 丝袜喷水一区| 欧美xxⅹ黑人| 超碰av人人做人人爽久久| 18禁动态无遮挡网站| 99久久九九国产精品国产免费| 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 99热这里只有是精品50| 2022亚洲国产成人精品| 毛片一级片免费看久久久久| 成年女人在线观看亚洲视频 | 青青草视频在线视频观看| 又爽又黄无遮挡网站| 久久久久国产网址| 在线免费观看的www视频| 亚洲国产精品成人久久小说| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 美女内射精品一级片tv| 免费黄频网站在线观看国产| 亚洲精品视频女| 天美传媒精品一区二区| 在线观看人妻少妇| 亚洲av日韩在线播放| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 黄色日韩在线| 国国产精品蜜臀av免费| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 色综合站精品国产| 久久久久久久国产电影| 中文天堂在线官网| 免费观看无遮挡的男女| 联通29元200g的流量卡| 欧美日韩精品成人综合77777| 欧美日韩国产mv在线观看视频 | 久久精品国产鲁丝片午夜精品| 欧美三级亚洲精品| 超碰av人人做人人爽久久| 高清毛片免费看| 肉色欧美久久久久久久蜜桃 | 秋霞在线观看毛片| 国产精品一区二区在线观看99 | av黄色大香蕉| 91久久精品电影网| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆| 国产黄色免费在线视频| 性插视频无遮挡在线免费观看| 久久久a久久爽久久v久久| 在线观看av片永久免费下载| 国产免费视频播放在线视频 | 日韩人妻高清精品专区| 国内精品美女久久久久久| 亚洲成人精品中文字幕电影| 亚洲欧洲国产日韩| 18禁在线播放成人免费| 人妻少妇偷人精品九色| 十八禁国产超污无遮挡网站| 亚洲最大成人中文| 日韩精品青青久久久久久| 婷婷色综合大香蕉| 日本免费a在线| 亚洲成人久久爱视频| 亚洲在线观看片| 少妇高潮的动态图| 亚洲成人中文字幕在线播放| 亚洲欧美日韩卡通动漫| 国产欧美日韩精品一区二区| 久久精品久久久久久噜噜老黄| 禁无遮挡网站| 欧美激情国产日韩精品一区| 国产男人的电影天堂91| 国产精品一二三区在线看| 成人亚洲精品一区在线观看 | 欧美精品一区二区大全| 亚洲aⅴ乱码一区二区在线播放| 在线观看人妻少妇| 亚洲真实伦在线观看| 午夜福利视频精品| 免费看不卡的av| 最新中文字幕久久久久| 非洲黑人性xxxx精品又粗又长| 久久鲁丝午夜福利片| 国产黄频视频在线观看| 嫩草影院精品99| 国产精品久久久久久精品电影小说 | 91精品一卡2卡3卡4卡| 亚洲精品乱码久久久久久按摩| 国产av不卡久久| 亚洲精品国产av成人精品| av线在线观看网站| 亚洲精华国产精华液的使用体验| 身体一侧抽搐| 国产精品人妻久久久影院| 欧美xxxx黑人xx丫x性爽| 伊人久久国产一区二区| 在线a可以看的网站| 精品久久久久久久久亚洲| 日日干狠狠操夜夜爽| 尤物成人国产欧美一区二区三区| 色网站视频免费| 亚洲av成人av| 国产老妇伦熟女老妇高清| 亚洲在线观看片| 国产精品无大码| 2021少妇久久久久久久久久久| 亚洲性久久影院| 美女黄网站色视频| 久久久久久久大尺度免费视频| 我的老师免费观看完整版| 免费观看性生交大片5| 黄色欧美视频在线观看| 亚洲美女视频黄频| 麻豆av噜噜一区二区三区| 国产熟女欧美一区二区| 两个人视频免费观看高清| 精品一区二区三卡| 91午夜精品亚洲一区二区三区| 我要看日韩黄色一级片| 天天躁日日操中文字幕| 午夜福利在线在线| 五月玫瑰六月丁香| 一个人观看的视频www高清免费观看| 成人av在线播放网站| 国产在线一区二区三区精| 精品少妇黑人巨大在线播放| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 亚洲国产色片| 免费黄频网站在线观看国产| 国产白丝娇喘喷水9色精品| 成年女人在线观看亚洲视频 | 婷婷色av中文字幕| 免费av观看视频| 中文字幕av成人在线电影| 亚洲国产欧美在线一区| 午夜免费激情av| 伊人久久国产一区二区| 久久久久久久久久人人人人人人| 热99在线观看视频| 美女内射精品一级片tv| 伊人久久国产一区二区| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影| 午夜福利在线观看吧| 久久久久久久久久久丰满| 精品不卡国产一区二区三区| 日本wwww免费看| 插阴视频在线观看视频| 亚洲在久久综合| 精品久久久久久久久av| 天堂网av新在线| 又大又黄又爽视频免费| 久久草成人影院| 边亲边吃奶的免费视频| 国内精品宾馆在线| 夫妻性生交免费视频一级片| 国产免费视频播放在线视频 | 在线免费观看的www视频| 久久久久网色| 国产精品不卡视频一区二区| 亚洲欧美精品专区久久| 亚洲av电影在线观看一区二区三区 | 欧美另类一区| 免费黄色在线免费观看| 国产精品人妻久久久久久| 亚洲欧美精品专区久久| 午夜亚洲福利在线播放| 久热久热在线精品观看| 国产精品一区二区在线观看99 | 熟女人妻精品中文字幕| 中文字幕av在线有码专区| 午夜日本视频在线| 欧美另类一区| 如何舔出高潮| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| 大香蕉久久网| 久久99热这里只频精品6学生| 免费少妇av软件| 欧美精品一区二区大全| 成人国产麻豆网| 在线观看免费高清a一片| 99热这里只有是精品在线观看| 免费av观看视频| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区成人| 免费黄频网站在线观看国产| ponron亚洲| 舔av片在线| 中文字幕亚洲精品专区| 大陆偷拍与自拍| 精品午夜福利在线看| 国产精品一区www在线观看| 久久亚洲国产成人精品v| 欧美日韩精品成人综合77777| 免费观看的影片在线观看| 亚洲人成网站在线观看播放| 国产一区二区在线观看日韩| 2021少妇久久久久久久久久久| 亚洲av中文字字幕乱码综合| 国产黄色免费在线视频| 亚洲av一区综合| 最近的中文字幕免费完整| 在线播放无遮挡| 久久精品国产自在天天线| 国产亚洲精品久久久com| 可以在线观看毛片的网站| 久久99精品国语久久久| 青春草亚洲视频在线观看| 中文字幕制服av| 欧美区成人在线视频| 国产成年人精品一区二区| 水蜜桃什么品种好| 一个人看视频在线观看www免费| 久久久成人免费电影| 丝袜美腿在线中文| 久久精品国产鲁丝片午夜精品| 身体一侧抽搐| 久久久久久久午夜电影| 久久99热6这里只有精品| 国产精品.久久久| 日韩欧美三级三区| 非洲黑人性xxxx精品又粗又长| 欧美激情久久久久久爽电影| 久久久久久久久久黄片| 天堂俺去俺来也www色官网 | 国产一区二区在线观看日韩| 一二三四中文在线观看免费高清| 高清在线视频一区二区三区| 简卡轻食公司| 大又大粗又爽又黄少妇毛片口| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 观看免费一级毛片| 日韩一本色道免费dvd| 69人妻影院| 亚洲在线观看片| 中文欧美无线码| 女人十人毛片免费观看3o分钟| 99热这里只有是精品50| 国产三级在线视频| 久久韩国三级中文字幕| 午夜精品一区二区三区免费看| av在线观看视频网站免费| 一个人观看的视频www高清免费观看| 青春草视频在线免费观看| 亚洲真实伦在线观看| 精品久久久精品久久久| 三级经典国产精品| 精品久久久久久久末码| 久久久久久久久久人人人人人人| 一级毛片我不卡| 大陆偷拍与自拍| 欧美成人一区二区免费高清观看| 国产av码专区亚洲av| 亚洲婷婷狠狠爱综合网| 三级男女做爰猛烈吃奶摸视频| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 亚洲美女视频黄频| 国产淫片久久久久久久久| 国产成人午夜福利电影在线观看| 欧美日韩精品成人综合77777| 99久久精品国产国产毛片| 色5月婷婷丁香| 日日撸夜夜添| 国产视频内射| 精品国产露脸久久av麻豆 | 国产精品国产三级国产专区5o| 精品酒店卫生间| 国产一区二区亚洲精品在线观看| 中文字幕亚洲精品专区| 午夜福利高清视频| 久久韩国三级中文字幕| 99久久精品一区二区三区| 国产伦精品一区二区三区视频9| 国产黄色视频一区二区在线观看| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 成年av动漫网址| 2022亚洲国产成人精品| 国产精品一区www在线观看| 欧美精品一区二区大全| 一级黄片播放器| 只有这里有精品99| 国产单亲对白刺激| 成人国产麻豆网| 亚洲av在线观看美女高潮| 亚洲国产av新网站| 国产精品麻豆人妻色哟哟久久 | 中国美白少妇内射xxxbb| 舔av片在线| 国产成人一区二区在线| 日本av手机在线免费观看| 特级一级黄色大片| 欧美激情在线99| 黄色日韩在线| 午夜精品在线福利| 亚洲精品乱码久久久久久按摩| 秋霞在线观看毛片| 免费观看a级毛片全部| 欧美日韩亚洲高清精品| 草草在线视频免费看| 精品久久久久久久久久久久久| 免费电影在线观看免费观看| 成年女人在线观看亚洲视频 | 韩国av在线不卡| 只有这里有精品99| 国产91av在线免费观看| 免费电影在线观看免费观看| 可以在线观看毛片的网站| 91久久精品国产一区二区三区| 91精品一卡2卡3卡4卡| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 女人十人毛片免费观看3o分钟| 黄色日韩在线| 好男人在线观看高清免费视频| 午夜福利网站1000一区二区三区| 好男人在线观看高清免费视频| 亚洲经典国产精华液单| av在线观看视频网站免费| 国产国拍精品亚洲av在线观看| 免费看日本二区| 全区人妻精品视频| 一级二级三级毛片免费看| 国产淫语在线视频| 亚洲一级一片aⅴ在线观看| 美女高潮的动态| 国产高潮美女av| 久久97久久精品| 午夜免费男女啪啪视频观看| eeuss影院久久| 国产美女午夜福利| 成年女人看的毛片在线观看| 搡老妇女老女人老熟妇| 欧美潮喷喷水| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕人妻熟人妻熟丝袜美| 久久99热6这里只有精品| 寂寞人妻少妇视频99o| 国产成人a∨麻豆精品| 91精品伊人久久大香线蕉| 中文字幕久久专区| 日韩不卡一区二区三区视频在线| 久久精品国产亚洲网站| 三级国产精品欧美在线观看| 亚洲国产最新在线播放| 国产精品99久久久久久久久| 久久久久久国产a免费观看| 国产一区有黄有色的免费视频 | 色哟哟·www| 精品久久久噜噜| 亚洲av男天堂| 亚洲久久久久久中文字幕| 日韩欧美精品v在线| 99热6这里只有精品| 22中文网久久字幕| 欧美成人一区二区免费高清观看| 看免费成人av毛片| 超碰av人人做人人爽久久| 亚洲欧美成人精品一区二区| 国产精品日韩av在线免费观看| 又黄又爽又刺激的免费视频.| 亚洲国产av新网站| 内地一区二区视频在线| 日韩欧美精品v在线| 久久精品国产亚洲av天美| 亚洲精品色激情综合| 欧美高清成人免费视频www| 亚洲欧美日韩卡通动漫| 寂寞人妻少妇视频99o| 久久久久性生活片| 精品国产三级普通话版| 啦啦啦韩国在线观看视频| 亚洲av一区综合| av福利片在线观看| 国产人妻一区二区三区在| 深夜a级毛片| 国产在视频线在精品| 夫妻午夜视频| 欧美激情久久久久久爽电影| 国产精品嫩草影院av在线观看| 国产极品天堂在线| 丝瓜视频免费看黄片| 天堂√8在线中文| 高清视频免费观看一区二区 | 亚洲精品视频女| 亚洲国产成人一精品久久久| 国产男人的电影天堂91| 熟妇人妻久久中文字幕3abv| 丰满少妇做爰视频| 91午夜精品亚洲一区二区三区| 如何舔出高潮| 欧美精品国产亚洲| 免费看光身美女| eeuss影院久久| 免费少妇av软件| 美女xxoo啪啪120秒动态图| 黄色一级大片看看| av女优亚洲男人天堂| 91av网一区二区| 亚洲人成网站高清观看| 久久综合国产亚洲精品| 全区人妻精品视频| 国产亚洲一区二区精品| 色吧在线观看| 国产伦一二天堂av在线观看| 久久6这里有精品| 久久久久国产网址| 亚洲精品影视一区二区三区av| 成人毛片a级毛片在线播放| 成年免费大片在线观看| 久久国产乱子免费精品| 中文字幕免费在线视频6| 99久久中文字幕三级久久日本| 97超碰精品成人国产| 午夜福利在线观看吧| 亚洲精品亚洲一区二区| 国产在线一区二区三区精| 夫妻午夜视频| 人人妻人人澡人人爽人人夜夜 | 亚洲欧美一区二区三区黑人 | 国产一区亚洲一区在线观看| 久久亚洲国产成人精品v| 午夜老司机福利剧场| 国产精品国产三级专区第一集| 日日啪夜夜爽| 久久久精品欧美日韩精品| 亚洲av二区三区四区| 久久精品夜色国产| 久久这里只有精品中国| 国产激情偷乱视频一区二区| 一二三四中文在线观看免费高清| ponron亚洲| 欧美高清性xxxxhd video| 久久精品综合一区二区三区| 波野结衣二区三区在线| 麻豆乱淫一区二区| 美女黄网站色视频| 精品少妇黑人巨大在线播放| 亚洲精品日本国产第一区| 美女主播在线视频| 久久这里有精品视频免费| 91精品伊人久久大香线蕉| 日韩一本色道免费dvd| 亚洲av在线观看美女高潮| 国产精品爽爽va在线观看网站| 国产老妇女一区| 成年免费大片在线观看| 亚洲自偷自拍三级| 国产亚洲精品av在线| 两个人视频免费观看高清| 亚洲成色77777| 久久国内精品自在自线图片| 亚洲精品aⅴ在线观看| 97精品久久久久久久久久精品| 国产乱人视频| 国产69精品久久久久777片| av免费观看日本| 亚洲av在线观看美女高潮| 亚州av有码| 秋霞伦理黄片| 18禁在线播放成人免费| 亚洲在久久综合| 日日摸夜夜添夜夜爱| 大香蕉久久网| 男插女下体视频免费在线播放| 午夜福利在线观看免费完整高清在| 2021少妇久久久久久久久久久| 1000部很黄的大片| 三级毛片av免费| 亚洲自拍偷在线| 蜜桃久久精品国产亚洲av| av免费在线看不卡| 成人欧美大片| 久久综合国产亚洲精品| 岛国毛片在线播放| 日日干狠狠操夜夜爽| 女人久久www免费人成看片| 精品久久久精品久久久| 一区二区三区高清视频在线| 韩国高清视频一区二区三区| 国产不卡一卡二| 国产综合懂色| 色综合站精品国产| 国产精品99久久久久久久久| 国产精品久久视频播放| 身体一侧抽搐| 国产精品一二三区在线看| 日韩精品有码人妻一区| 日本一本二区三区精品| 麻豆国产97在线/欧美| 麻豆成人av视频| 精品久久久久久久末码| 成人毛片60女人毛片免费| 麻豆精品久久久久久蜜桃| 男女下面进入的视频免费午夜| 免费黄频网站在线观看国产| 男女国产视频网站| 成人高潮视频无遮挡免费网站| 婷婷色综合www| 免费高清在线观看视频在线观看| 久久久久久久大尺度免费视频| 美女黄网站色视频| 91精品一卡2卡3卡4卡| 丝袜美腿在线中文| 久久久久久九九精品二区国产| 日韩欧美国产在线观看| 久久精品久久精品一区二区三区| 亚洲精品乱码久久久久久按摩| 亚洲人成网站高清观看| 天天躁日日操中文字幕|