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

    變化環(huán)境下岷江黑水河流域徑流演變特征分析

    2022-04-20 04:16:00王文圣余思怡
    人民珠江 2022年4期
    關(guān)鍵詞:岷江年際徑流

    王 寧,王文圣,余思怡

    (四川大學(xué)水利水電學(xué)院,四川 成都 610065)

    全球及區(qū)域氣候變化顯著,受人類社會影響的流域下墊面變化劇烈,導(dǎo)致河川徑流在時空上都呈現(xiàn)出明顯的變異特征[1-4]。河川徑流是中國主要的水資源,關(guān)系著區(qū)域生態(tài)環(huán)境、農(nóng)業(yè)灌溉以及社會經(jīng)濟發(fā)展和穩(wěn)定。因此,研究和明晰變化環(huán)境下的徑流演變特征,對水資源的合理開發(fā)、優(yōu)化配置及促進區(qū)域社會經(jīng)濟發(fā)展具有重要的意義[5]。

    水文學(xué)者對徑流演變特征開展了大量的研究。張建云等[6]研究表明,1956—2018年間中國主要江河徑流量均不同程度減少。王道席等[7]歸納了黃河源區(qū)徑流演變研究進展,黃河源區(qū)徑流總體呈下降趨勢,年徑流存在3個周期且在1990年左右發(fā)生突變。Chen等[8]研究了珠江流域徑流多尺度變化特征,指出年最小流量呈增加趨勢,年最大流量相比年最小流量的標(biāo)準(zhǔn)偏差較大,這不利于珠江三角洲的洪水減緩,極端水事件如干旱或洪水,發(fā)生的可能性較高。王文圣等[9]探討了長江上游主要河流年徑流的變化特性,嘉陵江、岷江、沱江和涪江的年徑流呈下降趨勢,金沙江年徑流無明顯的趨勢變化,主要河流年徑流存在長度不等的周期震蕩特征,多數(shù)測站間的年徑流在空間上具有豐枯同步變化特點。梁淑琪等[10]研究表明岷江上游徑流年內(nèi)分配逐年由均勻向不均勻再向均勻轉(zhuǎn)變,年徑流總體呈顯著下降趨勢,且下降趨勢在未來具有持續(xù)性。趙軍凱等[11]分析了宜昌站徑流量年際變化、年內(nèi)變化及三峽工程運行對宜昌徑流變化的影響,宜昌汛期徑流呈顯著減少趨勢,枯季則呈不明顯增加趨勢,徑流年內(nèi)分配趨于均勻,三峽工程運行加劇了長江上游徑流汛期減少與枯季增加的趨勢,使年內(nèi)分配差異減小。近年來,環(huán)境變化特別是降水和水電開發(fā)的影響使得徑流的影響發(fā)生明顯變化[12-13],進一步影響區(qū)域水資源的管理。

    岷江上游是長江流域的重要生態(tài)屏障,也是成都平原的重要水源,岷江上游的水電開發(fā)進度在1990年之后逐漸加快,隨著岷江上游水電站數(shù)量的增加,該流域天然狀態(tài)下的徑流變化更加顯著[14-15]。黑水河為岷江上游典型流域,水能資源相對豐富,流域氣候變化敏感,其徑流的演變趨勢、變異特征、周期變化、環(huán)境變化特別是降水和水電開發(fā)影響下徑流突變前后年際、年內(nèi)變化特征研究,事關(guān)成都平原水資源開發(fā)利用與管理。研究成果為提高對岷江上游徑流變化規(guī)律的認(rèn)識水平,以期為區(qū)域水資源利用及水電開發(fā)決策提供理論依據(jù)。

    1 研究流域概況與數(shù)據(jù)

    黑水河是岷江上游右岸的一級支流,主要流經(jīng)阿壩州黑水縣和茂縣(31°42′~32°54′N,102°35′~103°41′E,見圖1)。黑水河上游有西、北兩源,西源發(fā)源于黑水縣西部的羊拱山麓,北源發(fā)源于毛爾蓋草原,干流長122 km,流域面積為7 240 km2,落差達(dá)1 048 m,平均比降8.6‰,河口多年平均流量為140 m3/s[16]。黑水河具有高山峽谷區(qū)的主要地貌,氣候變化敏感,徑流補給來源以降水為主。

    圖1 黑水河流域水系

    沙壩水文站為黑水河流域控制出口站,控制流域面積為7 218 km2,有1960—2015年實測徑流資料供研究使用。

    2 研究方法

    2.1 突變點診斷法

    采用改進的啟發(fā)式分割算法對年徑流序列的突變點進行診斷,該方法基于雙重統(tǒng)計檢驗進行人工干預(yù),能夠快速且準(zhǔn)確找到序列中的跳躍點。在啟發(fā)式分割算法的基礎(chǔ)上,引入Z檢驗和SNHT檢驗(標(biāo)準(zhǔn)正態(tài)均一性檢驗)[17],提出了改進的啟發(fā)式分割算法,并將其用于年徑流時間序列的突變點檢測。改進算法如下。

    步驟一將序列x1,x2,…,xn標(biāo)準(zhǔn)化為Ri(i=2,3,…,n)。

    步驟二將Ri分割為前后2個子序列,分割點記為i(3≤i≤n-3),分割點前后序列記為Rt(i)、Rr(i),對應(yīng)的均值分別為μt(i)、μr(i),方差分別為st(i)、sr(i)。按照Z檢驗法構(gòu)建統(tǒng)計量Z0為

    (1)

    其中,nt、nr為前后兩子序列的長度。

    步驟三給定顯著性水平α,若|Z0|>Zα/2,則表明序列x1,x2,…,xn可能在分割點i處發(fā)生突變。利用SNHT檢驗法對突變點進行顯著性檢驗,構(gòu)建統(tǒng)計量T0=nt×μt(i)2+nr×μr(i)2,根據(jù)顯著性水平α和序列長度,在文獻(xiàn)[18]中查找臨界值,若T0大于臨界值,則i點為序列的顯著突變點。

    步驟四i=i+1,返回執(zhí)行步驟二。

    2.2 趨勢性識別方法

    利用線性傾向估計法、Mann-Kendall趨勢檢驗法、R/S分析法對徑流序列進行趨勢分析,線性傾向估計法[19]原理不在此贅述。

    2.2.1Mann-Kendall趨勢檢驗法

    Mann-Kendall趨勢檢驗法假定序列x1,x2,…,xn是獨立的,定義中間變量S并構(gòu)造統(tǒng)計量Z如下[20]:

    (2)

    (3)

    式中,sgn為符號函數(shù),同上。給定顯著性水平α,若|Z|≥Zα/2,則判定序列變化趨勢顯著,反之不顯著;Z>0,表明序列呈上升趨勢,反之呈下降趨勢。

    2.2.2R/S分析法

    R/S分析法是利用Hurst指數(shù)來分析水文序列的持續(xù)性和反持續(xù)性,從而揭示將來的變化趨勢[19-20]。

    將序列x1,x2,…,xn等分為長度為τ(3≤τ≤n)的M個(M為n/τ的整數(shù)部分)連續(xù)子序列,記Ik(k=1,2,…,M),子序列的元素為xk,m(m=1,2,…,τ)。計算Ik序列的均值Ek、累計離差Ck,a、極差RIa及其標(biāo)準(zhǔn)差SIa。

    (4)

    (5)

    RIa=max(Ck,a)-min(Ck,a)

    (6)

    (7)

    計算Ik序列的重標(biāo)極差RSk,并計算M個連續(xù)子序列重標(biāo)極差的平均值RS。

    (8)

    (9)

    重復(fù)以上計算過程,直至τ=n,得到(RS)n序列,繪制以ln(n)為自變量、ln[(RS)n]為因變量的散點圖,以最小二乘法擬合直線,該直線的斜率即為Hurst指數(shù)值H。

    若0.5

    2.3 周期性識別方法

    小波分析法是一種通過調(diào)整時頻窗來識別序列周期成分的方法[21]。對于年徑流序列,其連續(xù)小波變換定義為:

    (10)

    小波方差Var(a)可以識別序列的周期成分,計算公式為:

    (11)

    小波方差隨尺度a變化的過程稱為小波方差圖,圖上波峰對應(yīng)的尺度為時間序列的主周期。

    3 結(jié)果及其分析

    3.1 徑流突變性分析

    應(yīng)用改進的啟發(fā)式分割算法對黑水河年徑流突變年份進行診斷,根據(jù)顯著性水平(α=5%)和序列長度(n=56)查找臨界值為8.566,當(dāng)分割點i=1995年時,統(tǒng)計量T0=8.762,統(tǒng)計量計算結(jié)果大于臨界值,突變點通過顯著性檢驗,說明黑水河年徑流在1995年發(fā)生顯著突變。

    3.2 突變前后徑流趨勢性分析

    圖2繪制了黑水河年徑流變化曲線,并給出了年徑流5 a滑動平均過程,表1列出了年徑流在1960—2015年、突變點前(1960—1995年)、突變點后(1996—2015年)的趨勢分析結(jié)果。從圖2和表1可知,年徑流在1960—2015年間以-0.32 m3/(s·a)的速率呈顯著下降趨勢,在突變點前以0.12 m3/(s·a)的速率呈不顯著下降趨勢,這與整體趨勢一致,但變化速率變?。煌蛔凕c后的徑流變化趨勢與整體趨勢相反,年徑流以0.80 m3/(s·a)的速率呈不顯著上升趨勢。

    圖2 黑水河年徑流變化曲線

    表1 1960—2015年黑水河年徑流趨勢分析結(jié)果

    R/S分析法計算得到1960—2015年間年徑流Hurst指數(shù)H為0.64,大于0.5,表明黑水河年徑流在未來一段時間內(nèi)變化趨勢可能與該時期的變化趨勢一致,即呈下降的變化(圖3)。

    圖3 黑水河年徑流R/S分析法計算結(jié)果

    3.3 突變前后徑流年際與年內(nèi)變化特征分析

    黑水河流域重要的人類活動之一是水電能源開發(fā),黑水河干流馬橋至白溪河段,按照 “二庫五級”規(guī)劃方案開發(fā),其中馬橋水庫和毛爾蓋水庫的調(diào)節(jié)庫容為2.8億m3[22],自上而下梯級電站分別為馬橋、竹格多、毛爾蓋、色爾古和柳坪電站,各梯級電站基本情況見表2。黑水河徑流補給以降水為主,已有學(xué)者研究發(fā)現(xiàn)黑水河年降水突變年份為1995年[23],與徑流突變一致,由此這里探討降水和水電開發(fā)影響下黑水河徑流在突變點前后的年際、年內(nèi)變化情況。表3給出了突變點前后及整個研究時段年徑流和年降水的均值和Cv值,圖4給出了突變點前后月平均流量分配過程及其變幅。表3和圖4表明:①年徑流突變點前的均值為141 m3/s,變點后的均值為128 m3/s,前后相差較大,減少9.22%,突變點前Cv值為0.10,變點后Cv值為0.13,突變后年徑流年際變化較突變前大;②年降水量的年際變化與年徑流量變化一致,突變點后均值較突變點前小,減少4.63%,突變點后的Cv值較突變前大,突變后年降水年際變化較突變前大;③突變前后徑流的年內(nèi)分配形式略有不同,突變前呈“雙峰型”,突變后呈“單峰型”;④突變前后徑流的年內(nèi)變幅呈現(xiàn)明顯的差異,突變后汛期和非汛期的最大月平均流量、最小月平均流量、月平均流量較突變前均不同程度減少,其中汛期的徑流量變幅較非汛期大。

    表2 黑水河干流二庫五級方案各梯級水電站基本情況

    表3 不同時期黑水河徑流和降水年際變化指標(biāo)計算結(jié)果

    圖4 突變點前和突變點后黑水河年內(nèi)月徑流分配過程

    以上研究說明,受降水、水電開發(fā)的影響,黑水河徑流在1995年前后的演變特征不同,這種不同是顯著的。黑水河汛期徑流量的減少有效減輕了該時段的防洪壓力,但非汛期的最小月平均流量的減少不利于枯季水資源利用。

    3.4 突變前后徑流周期性分析

    利用小波分析法分別計算1960—2015、1960—1995、1996—2015年的年徑流序列的多時間尺度變化和小波方差變化過程(圖5)??梢钥闯?,黑水河年徑流演變過程中呈現(xiàn)多時間尺度變化特征,不同時期時間尺度的強弱分布各有不同。年徑流在1960—2015年間存在著6、16、28 a的主周期,其中6 a左右的周期震蕩最強,為第一主周期;年徑流在1960—1995年間存在7、16 a的主周期,在1996—2015年間存在2、5 a的主周期。

    a)1960—2015年年徑流小波系數(shù)時頻分布

    d)1960—1995年年徑流小波方差

    4 結(jié)論

    基于黑水河沙壩水文站1960—2015年實測徑流資料,采用改進的啟發(fā)式分割算法、Mann-Kendall趨勢檢驗法、R/S分析法和小波分析法等對年徑流序列的演變特征進行研究,主要結(jié)論如下。

    a)黑水河年徑流在1995年發(fā)生突變。

    b)年徑流在1960—2015年以0.32 m3/s/a的速率呈顯著下降趨勢,且在未來一段時間內(nèi)仍可能保持下降趨勢,在1995年前后兩個時期的變化趨勢不顯著。

    c)受環(huán)境變化特別是降水和水電開發(fā)的影響,黑水河徑流突變前均值較突變后大,突變后的年際變化較突變前大;年內(nèi)分配形式在突變點前表現(xiàn)為雙峰,在變點后為單峰;突變前后徑流的年內(nèi)變幅呈現(xiàn)明顯的差異,突變后的汛期徑流量顯著減少,可以有效減輕該流域的防洪壓力,但非汛期最小月平均流量的減少不利于枯季水資源利用。

    d)年徑流呈現(xiàn)多時間尺度變化特征,在1960—2015年間存在6、16、28 a的主周期,在1960—1995年間存在7、16 a的主周期,在1996—2015年間存在2、5 a的主周期。

    猜你喜歡
    岷江年際徑流
    岷江
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    岷江行(外一首)
    岷峨詩稿(2019年4期)2019-04-20 09:02:06
    瘋狂的蟲草,瘋狂的松茸和瘋狂的岷江柏*——專訪作家阿來
    阿來研究(2018年1期)2018-07-10 09:39:08
    岷江同大渡河相會樂山
    廉政瞭望(2016年13期)2016-08-11 11:22:02
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    探秘“大徑流”
    攻克“大徑流”
    亚洲精品久久久久久婷婷小说| 亚洲国产精品成人久久小说| 日韩欧美国产在线观看| 成人亚洲精品一区在线观看 | 成人一区二区视频在线观看| 永久免费av网站大全| 国内揄拍国产精品人妻在线| 国产av国产精品国产| 麻豆av噜噜一区二区三区| 天堂av国产一区二区熟女人妻| 美女脱内裤让男人舔精品视频| 亚洲真实伦在线观看| 一级a做视频免费观看| 18禁在线无遮挡免费观看视频| 晚上一个人看的免费电影| 日日摸夜夜添夜夜添av毛片| 午夜日本视频在线| 舔av片在线| 国产麻豆成人av免费视频| 日本一本二区三区精品| av.在线天堂| 免费av毛片视频| 观看美女的网站| 建设人人有责人人尽责人人享有的 | 男女下面进入的视频免费午夜| 偷拍熟女少妇极品色| 久久久久性生活片| 欧美日本视频| 亚洲精品成人久久久久久| 亚洲av一区综合| a级毛色黄片| 免费看不卡的av| 91午夜精品亚洲一区二区三区| 九色成人免费人妻av| 在现免费观看毛片| 日韩欧美一区视频在线观看 | 一级av片app| 亚洲精品一二三| 亚洲av中文av极速乱| 久久精品综合一区二区三区| 一级黄片播放器| 80岁老熟妇乱子伦牲交| 只有这里有精品99| 嫩草影院新地址| 国产精品一二三区在线看| 中文在线观看免费www的网站| 七月丁香在线播放| 最近最新中文字幕免费大全7| 国国产精品蜜臀av免费| 国产成人a∨麻豆精品| 久久久色成人| 特大巨黑吊av在线直播| 午夜久久久久精精品| 小蜜桃在线观看免费完整版高清| 国产伦在线观看视频一区| 免费大片黄手机在线观看| 波野结衣二区三区在线| 欧美xxⅹ黑人| 91午夜精品亚洲一区二区三区| 欧美日韩在线观看h| 菩萨蛮人人尽说江南好唐韦庄| 亚洲熟妇中文字幕五十中出| 久久精品综合一区二区三区| 国精品久久久久久国模美| 春色校园在线视频观看| 国产v大片淫在线免费观看| 精品久久久噜噜| 亚洲欧美日韩东京热| 国产免费一级a男人的天堂| av在线蜜桃| 国产真实伦视频高清在线观看| 一级a做视频免费观看| 国产精品久久久久久精品电影小说 | 久久久久久久亚洲中文字幕| 亚洲无线观看免费| 永久网站在线| 国产免费福利视频在线观看| 免费av毛片视频| 一二三四中文在线观看免费高清| 18禁裸乳无遮挡免费网站照片| 午夜激情欧美在线| 亚洲av电影不卡..在线观看| 国产成人a区在线观看| 大香蕉久久网| 久久热精品热| av又黄又爽大尺度在线免费看| 成人无遮挡网站| 国产精品熟女久久久久浪| 亚洲欧美清纯卡通| 亚洲美女搞黄在线观看| 国产淫语在线视频| 天堂中文最新版在线下载 | 特级一级黄色大片| 黄色一级大片看看| 全区人妻精品视频| 中文天堂在线官网| 热99在线观看视频| 一级毛片电影观看| videossex国产| 最近2019中文字幕mv第一页| 内地一区二区视频在线| 国产欧美日韩精品一区二区| 99久久九九国产精品国产免费| 最近中文字幕2019免费版| 男女边摸边吃奶| 菩萨蛮人人尽说江南好唐韦庄| 青春草亚洲视频在线观看| 婷婷色综合www| 日本三级黄在线观看| 免费少妇av软件| 中文天堂在线官网| 日日摸夜夜添夜夜爱| av线在线观看网站| 精品人妻偷拍中文字幕| 免费看美女性在线毛片视频| 亚洲av电影不卡..在线观看| 永久免费av网站大全| 国产高清国产精品国产三级 | 亚洲精品一二三| 国产精品福利在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 日本三级黄在线观看| 日本-黄色视频高清免费观看| 床上黄色一级片| 国产乱来视频区| 天天一区二区日本电影三级| 欧美精品国产亚洲| 亚洲精品456在线播放app| 女人久久www免费人成看片| 久久久久久久久大av| 男人和女人高潮做爰伦理| a级一级毛片免费在线观看| 人人妻人人澡欧美一区二区| 亚洲精品影视一区二区三区av| 久久久a久久爽久久v久久| 插逼视频在线观看| 一级二级三级毛片免费看| 欧美日韩亚洲高清精品| 精品一区二区三区人妻视频| 97在线视频观看| 国产视频内射| 亚洲电影在线观看av| 床上黄色一级片| 欧美丝袜亚洲另类| 欧美区成人在线视频| 国产亚洲91精品色在线| 久热久热在线精品观看| 国产免费一级a男人的天堂| 青春草亚洲视频在线观看| 成人亚洲精品av一区二区| .国产精品久久| 天堂√8在线中文| 99热全是精品| 国产老妇女一区| 国产国拍精品亚洲av在线观看| 韩国av在线不卡| 欧美不卡视频在线免费观看| 色综合站精品国产| 亚洲图色成人| freevideosex欧美| 禁无遮挡网站| 18+在线观看网站| 看非洲黑人一级黄片| 亚洲国产最新在线播放| 日日摸夜夜添夜夜添av毛片| 久久久久久久久久黄片| 国产单亲对白刺激| av在线蜜桃| 超碰av人人做人人爽久久| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看| 久久97久久精品| 亚洲天堂国产精品一区在线| 亚洲乱码一区二区免费版| 一级毛片久久久久久久久女| 最近手机中文字幕大全| av国产免费在线观看| 国产 亚洲一区二区三区 | 久久99热这里只有精品18| 深爱激情五月婷婷| 亚洲欧美中文字幕日韩二区| 99久久人妻综合| 少妇熟女aⅴ在线视频| 天天躁夜夜躁狠狠久久av| 亚洲自偷自拍三级| 午夜福利在线观看吧| 国产91av在线免费观看| 久久精品国产亚洲av涩爱| 亚洲精品,欧美精品| 国产精品日韩av在线免费观看| 在线免费十八禁| 免费少妇av软件| 午夜福利网站1000一区二区三区| 日本av手机在线免费观看| 午夜久久久久精精品| 亚洲av国产av综合av卡| 嫩草影院入口| 91狼人影院| 亚洲一级一片aⅴ在线观看| 国产中年淑女户外野战色| 三级经典国产精品| 欧美性猛交╳xxx乱大交人| 国产乱人偷精品视频| 日韩av不卡免费在线播放| 亚洲国产精品成人综合色| 成人一区二区视频在线观看| 久久99热这里只有精品18| 免费黄色在线免费观看| 午夜福利在线观看免费完整高清在| 日本免费a在线| 99热这里只有是精品50| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 可以在线观看毛片的网站| 黄色一级大片看看| 一级爰片在线观看| 欧美变态另类bdsm刘玥| 国产乱人视频| 日本猛色少妇xxxxx猛交久久| 能在线免费看毛片的网站| 国产v大片淫在线免费观看| 99re6热这里在线精品视频| 午夜福利在线在线| 久久久久久久午夜电影| 麻豆成人午夜福利视频| 干丝袜人妻中文字幕| 国产综合懂色| 伦理电影大哥的女人| 熟女电影av网| 免费电影在线观看免费观看| 亚洲高清免费不卡视频| 国产久久久一区二区三区| 三级经典国产精品| av免费在线看不卡| 韩国高清视频一区二区三区| 欧美精品一区二区大全| 好男人视频免费观看在线| 日韩成人伦理影院| 男女边吃奶边做爰视频| 日本一本二区三区精品| 99久久精品热视频| 嘟嘟电影网在线观看| 午夜福利高清视频| 毛片女人毛片| av天堂中文字幕网| 亚洲aⅴ乱码一区二区在线播放| 国产在线男女| 精品午夜福利在线看| 久久久久久久大尺度免费视频| 91久久精品电影网| 春色校园在线视频观看| 免费av不卡在线播放| 听说在线观看完整版免费高清| 丝袜喷水一区| 亚洲国产av新网站| 又爽又黄无遮挡网站| ponron亚洲| 午夜福利在线观看吧| 欧美人与善性xxx| 麻豆国产97在线/欧美| 丝瓜视频免费看黄片| 蜜臀久久99精品久久宅男| 日本一二三区视频观看| 欧美激情国产日韩精品一区| av一本久久久久| 久久久久久久国产电影| 国产 一区精品| 熟女人妻精品中文字幕| 永久网站在线| 美女cb高潮喷水在线观看| 国产毛片a区久久久久| 激情五月婷婷亚洲| 啦啦啦韩国在线观看视频| 热99在线观看视频| 亚洲精品国产av蜜桃| 汤姆久久久久久久影院中文字幕 | 男女那种视频在线观看| 久久久久久久久大av| 蜜桃亚洲精品一区二区三区| 国产亚洲5aaaaa淫片| 毛片女人毛片| 欧美三级亚洲精品| 一级毛片我不卡| 日本熟妇午夜| 舔av片在线| 一级毛片aaaaaa免费看小| 99久久精品热视频| 国产久久久一区二区三区| 深夜a级毛片| 免费黄色在线免费观看| 久久精品人妻少妇| 久99久视频精品免费| 夜夜看夜夜爽夜夜摸| 卡戴珊不雅视频在线播放| 亚洲无线观看免费| 偷拍熟女少妇极品色| 亚洲av中文字字幕乱码综合| 精品久久久久久久久av| 成年女人在线观看亚洲视频 | 久久久久九九精品影院| av天堂中文字幕网| 好男人视频免费观看在线| 亚洲最大成人中文| 啦啦啦中文免费视频观看日本| 97超碰精品成人国产| 色哟哟·www| 观看免费一级毛片| 国产精品一及| 国产在视频线在精品| 少妇丰满av| 美女高潮的动态| 午夜日本视频在线| 少妇人妻精品综合一区二区| 亚洲av中文av极速乱| 亚洲欧美中文字幕日韩二区| 成人毛片a级毛片在线播放| 午夜免费激情av| 十八禁网站网址无遮挡 | av福利片在线观看| 好男人视频免费观看在线| 岛国毛片在线播放| 午夜亚洲福利在线播放| 精品少妇黑人巨大在线播放| 久久99热这里只频精品6学生| 18+在线观看网站| 国产一区二区亚洲精品在线观看| 国产在线男女| 色综合站精品国产| 亚洲av成人精品一区久久| 亚洲av国产av综合av卡| 精品久久久久久久久av| 亚洲精品日韩av片在线观看| 日日干狠狠操夜夜爽| 亚洲不卡免费看| 两个人的视频大全免费| 在线观看美女被高潮喷水网站| 免费看日本二区| 在线天堂最新版资源| 国产精品日韩av在线免费观看| 中文字幕亚洲精品专区| 精品国产露脸久久av麻豆 | 嫩草影院入口| 亚洲人成网站在线观看播放| 国内揄拍国产精品人妻在线| 中文乱码字字幕精品一区二区三区 | 国产高潮美女av| 嫩草影院新地址| 黄色欧美视频在线观看| 黄片wwwwww| 简卡轻食公司| 内射极品少妇av片p| 日韩av免费高清视频| 国产成人a区在线观看| 日韩三级伦理在线观看| 天堂网av新在线| 亚洲一级一片aⅴ在线观看| 中文在线观看免费www的网站| 精品一区二区免费观看| 国产综合懂色| 国产亚洲91精品色在线| 午夜免费观看性视频| 高清av免费在线| 免费人成在线观看视频色| 国产一区亚洲一区在线观看| 成人一区二区视频在线观看| 精品一区在线观看国产| 国产亚洲精品av在线| 中国国产av一级| 亚洲精品一二三| 成人漫画全彩无遮挡| 欧美高清成人免费视频www| 91久久精品国产一区二区三区| 亚洲人成网站高清观看| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 狂野欧美白嫩少妇大欣赏| 久久精品国产亚洲av涩爱| 一级a做视频免费观看| 久久久久久久亚洲中文字幕| 网址你懂的国产日韩在线| 午夜激情久久久久久久| 欧美极品一区二区三区四区| 久久久亚洲精品成人影院| 能在线免费看毛片的网站| 十八禁网站网址无遮挡 | 亚洲av一区综合| 免费观看精品视频网站| 亚洲国产精品专区欧美| 麻豆av噜噜一区二区三区| 91久久精品电影网| 免费观看无遮挡的男女| 中文资源天堂在线| 国产精品一区二区三区四区久久| 精品亚洲乱码少妇综合久久| 国产一区二区在线观看日韩| 国产中年淑女户外野战色| 身体一侧抽搐| 精品久久久久久久末码| www.av在线官网国产| 夫妻午夜视频| 真实男女啪啪啪动态图| 最新中文字幕久久久久| 久久人人爽人人片av| 欧美三级亚洲精品| 性色avwww在线观看| 国产精品一区二区在线观看99 | 日韩亚洲欧美综合| 亚洲精品一二三| 精品人妻视频免费看| 精品人妻熟女av久视频| 大话2 男鬼变身卡| 毛片女人毛片| 精品不卡国产一区二区三区| 一夜夜www| 少妇的逼水好多| 亚洲国产欧美在线一区| 别揉我奶头 嗯啊视频| 亚洲欧美日韩无卡精品| 亚洲欧美中文字幕日韩二区| 欧美日韩一区二区视频在线观看视频在线 | 欧美日韩综合久久久久久| 国产黄色小视频在线观看| 精品久久久久久久人妻蜜臀av| 国产成年人精品一区二区| 夜夜看夜夜爽夜夜摸| 亚洲最大成人av| 毛片女人毛片| av线在线观看网站| 大话2 男鬼变身卡| 亚洲av电影在线观看一区二区三区 | 国产精品久久视频播放| 欧美潮喷喷水| 美女国产视频在线观看| 超碰97精品在线观看| 精品不卡国产一区二区三区| 视频中文字幕在线观看| 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 久久精品国产亚洲网站| 少妇熟女aⅴ在线视频| 久久人人爽人人爽人人片va| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 久久国内精品自在自线图片| 一个人看视频在线观看www免费| 国产成人福利小说| 伊人久久国产一区二区| 中文欧美无线码| 亚洲国产精品成人久久小说| 欧美日韩综合久久久久久| 国产免费视频播放在线视频 | 嫩草影院精品99| 中文字幕人妻熟人妻熟丝袜美| 久久国产乱子免费精品| 在线免费十八禁| 三级国产精品欧美在线观看| 性插视频无遮挡在线免费观看| 麻豆成人av视频| 久久久精品94久久精品| 国产综合懂色| 亚洲av成人精品一二三区| 嘟嘟电影网在线观看| 汤姆久久久久久久影院中文字幕 | 91精品一卡2卡3卡4卡| 亚洲av福利一区| 极品少妇高潮喷水抽搐| 日本色播在线视频| 麻豆成人午夜福利视频| 成人亚洲精品一区在线观看 | av线在线观看网站| 亚洲熟女精品中文字幕| 少妇熟女欧美另类| 欧美精品一区二区大全| 午夜免费激情av| 色综合色国产| 亚洲欧美日韩卡通动漫| 国产在线男女| 老女人水多毛片| 国产高清不卡午夜福利| 一区二区三区四区激情视频| 欧美日本视频| 丰满乱子伦码专区| 亚洲人成网站高清观看| 亚洲国产精品专区欧美| 久久久久久久久久久丰满| 男女下面进入的视频免费午夜| 啦啦啦啦在线视频资源| 蜜臀久久99精品久久宅男| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| av女优亚洲男人天堂| 九九在线视频观看精品| 91精品伊人久久大香线蕉| 大香蕉久久网| 搡老妇女老女人老熟妇| 久久久精品免费免费高清| 亚洲18禁久久av| 亚洲国产精品国产精品| 欧美日韩精品成人综合77777| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 成人性生交大片免费视频hd| 国产女主播在线喷水免费视频网站 | 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 日本黄色片子视频| freevideosex欧美| 久久久久性生活片| 汤姆久久久久久久影院中文字幕 | 久久久精品免费免费高清| 亚洲国产欧美在线一区| 又黄又爽又刺激的免费视频.| 国产成人免费观看mmmm| 久久久久久久国产电影| 久久精品夜夜夜夜夜久久蜜豆| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 亚洲一区高清亚洲精品| freevideosex欧美| 久久99热6这里只有精品| 亚洲欧美精品专区久久| 水蜜桃什么品种好| 亚洲成人av在线免费| 能在线免费看毛片的网站| av天堂中文字幕网| 18禁在线无遮挡免费观看视频| 欧美激情久久久久久爽电影| 免费观看性生交大片5| kizo精华| 欧美性感艳星| 久久久久国产网址| av在线播放精品| 精品少妇黑人巨大在线播放| 日韩成人av中文字幕在线观看| 欧美日本视频| 美女xxoo啪啪120秒动态图| 色5月婷婷丁香| 国产成人91sexporn| 99热这里只有精品一区| 黄色日韩在线| 亚洲精品国产成人久久av| 蜜桃久久精品国产亚洲av| 日本色播在线视频| 午夜精品在线福利| 一级毛片黄色毛片免费观看视频| 禁无遮挡网站| 肉色欧美久久久久久久蜜桃 | 97热精品久久久久久| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 大陆偷拍与自拍| 日本一二三区视频观看| 天天一区二区日本电影三级| 美女内射精品一级片tv| av天堂中文字幕网| 国产av国产精品国产| 欧美激情久久久久久爽电影| 免费观看性生交大片5| 18+在线观看网站| 最近视频中文字幕2019在线8| 大又大粗又爽又黄少妇毛片口| 国产三级在线视频| 免费在线观看成人毛片| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 亚洲精品,欧美精品| 亚洲国产成人一精品久久久| 亚洲怡红院男人天堂| 国产精品麻豆人妻色哟哟久久 | 大又大粗又爽又黄少妇毛片口| or卡值多少钱| 久久鲁丝午夜福利片| 欧美zozozo另类| h日本视频在线播放| 性色avwww在线观看| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 日产精品乱码卡一卡2卡三| 极品教师在线视频| 99久久精品热视频| 免费少妇av软件| 丰满乱子伦码专区| 尾随美女入室| 男插女下体视频免费在线播放| 夜夜看夜夜爽夜夜摸| 最近2019中文字幕mv第一页| 只有这里有精品99| 久久精品久久久久久久性| 最新中文字幕久久久久| 色尼玛亚洲综合影院| 久久这里只有精品中国| 成年av动漫网址| 街头女战士在线观看网站| 久久99热6这里只有精品| 只有这里有精品99| 麻豆av噜噜一区二区三区| 国产三级在线视频| 日本黄色片子视频| 国产av在哪里看| 日本黄色片子视频| 欧美日韩在线观看h| 亚洲精品中文字幕在线视频 | 亚洲最大成人中文| 一个人观看的视频www高清免费观看| 黄色日韩在线| 亚洲欧美清纯卡通| 我的老师免费观看完整版| 国产精品久久久久久精品电影小说 | 免费电影在线观看免费观看| 天堂网av新在线| 淫秽高清视频在线观看| 国产精品不卡视频一区二区| 综合色丁香网| 久久鲁丝午夜福利片| 97热精品久久久久久|