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

    考慮非平穩(wěn)性特征的雅礱江流域洪水量化分析

    2020-12-21 01:55:34潔,楊
    水力發(fā)電 2020年9期
    關(guān)鍵詞:洪峰流量時變水文站

    高 潔,楊 龍

    (1.水電水利規(guī)劃設(shè)計總院,北京 100120;2.南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023)

    0 引 言

    設(shè)計洪水是確定水電工程規(guī)模和制定防洪標(biāo)準(zhǔn)的重要基礎(chǔ)。基于頻率分析的設(shè)計洪水計算在分布線型、參數(shù)估計等方面積累了豐富的研究成果[1]。氣候變化和人類活動對水文時間序列的平穩(wěn)性特征產(chǎn)生影響,傳統(tǒng)水文設(shè)計方法面臨變化環(huán)境帶來的挑戰(zhàn)[2]。郭生練等[1]將非一致性洪水頻率分析、設(shè)計洪水不確定性研究、基于水文物理機(jī)制的洪水頻率分析等內(nèi)容列為當(dāng)前設(shè)計洪水研究前沿與熱點。魯帆等[3]在水文氣象領(lǐng)域非平穩(wěn)時間序列極值模型研究綜述中介紹,傳統(tǒng)水文頻率分析通常采用皮爾遜III型(P-III),對數(shù)皮爾遜III型,廣義極值分布(GEV),廣義帕累托分布(GPD)等理論線型。我國工程界設(shè)計洪水標(biāo)準(zhǔn)主要采用P-III分布。以年最大為樣本的GEV分布也在水文極值統(tǒng)計中逐漸被采用。關(guān)于GEV分布參數(shù),GEV-CDN模型[4]認(rèn)為形狀參數(shù)決定了分布的尾部形態(tài),具有較大的不確定性,可將形狀參數(shù)取值設(shè)定在[-0.5,0.5]。GEV模型應(yīng)用于洞庭湖區(qū)、岷江流域年最大日降水系列[5- 6],發(fā)現(xiàn)各地的形狀參數(shù)差異較大,形狀參數(shù)為0的極值I型Gumbel分布、形狀參數(shù)<0的極值III型Weibull分布,形狀參數(shù)>0的極值II型Frechet分布均有適用。通過對全國降雨資料統(tǒng)計,Yang等[7]發(fā)現(xiàn)極端降雨在長江流域呈現(xiàn)形狀參數(shù)<0的薄尾分布;在環(huán)渤海地區(qū)形狀參數(shù)>0,呈現(xiàn)厚尾分布。長江中下游控制站大通水文站[8]1950年至2012年歷年最大和非汛期最大流量序列符合形狀參數(shù)<0的Weibull分布,新疆天山開都河[9]1957年至2011年歷年、春季和夏季最大日流量序列服從形狀參數(shù)≥0的Gumbel分布和Frechet分布。通過對珠江流域北江水系武江上游坪石水文站和下游犁市(二)水文站實測年最大徑流系列分別采用薄尾P-III分布、Weibull分布、厚尾廣義邏輯分布(GLO)、三參數(shù)混合尾GEV、Gumbel、兩參數(shù)對數(shù)正態(tài)分布(LN2)[10]擬合發(fā)現(xiàn),厚尾GLO和三參數(shù)混合尾GEV分布在擬合時可使高水端上翹,提高了高水端擬合優(yōu)度。

    非平穩(wěn)時間序列的概率分布均值、方差或協(xié)方差隨時間變化,重現(xiàn)期和風(fēng)險度有別于傳統(tǒng)水文頻率計算結(jié)果,需要對變化序列進(jìn)行還原/還現(xiàn),或針對時變序列進(jìn)行頻率計算。謝平等[11]按照“分離”“合成”的思路對水文序列非一致性進(jìn)行診斷??紤]時間因素對傳統(tǒng)頻率分布的影響,劉德地等[12]推導(dǎo)了非平穩(wěn)條件下洪水重現(xiàn)期的計算公式。顧西輝等[13]選取珠江流域28個水文站1951年至2010年歷年最大流量數(shù)據(jù),在突變點分析的基礎(chǔ)上采用GEV、GPD、Gumbel、P-III、GLO、廣義對數(shù)正態(tài)分布(GNO)6種分布進(jìn)行擬合。研究發(fā)現(xiàn),最優(yōu)分布以GEV居多。珠江流域內(nèi)西江、北江頻率設(shè)計值增大,且隨著重現(xiàn)期增加而增大,基于整體樣本的水利水電工程規(guī)劃設(shè)計,可能低估大洪水影響,導(dǎo)致防洪風(fēng)險。顧西輝等[14]又將GEV-CDN模型[4]應(yīng)用于珠江流域6個水文站近45年歷年最大日流量序列。研究發(fā)現(xiàn),50年一遇以下洪水,平穩(wěn)和非平穩(wěn)模型成果差異不大;但50年一遇以上洪水,不同模型頻率設(shè)計值差異較大,存在潛在的防洪風(fēng)險。梁忠民等[15]針對非平穩(wěn)環(huán)境下工程水文設(shè)計值發(fā)生變化的問題,以等可靠度建立平穩(wěn)序列與非平穩(wěn)序列的聯(lián)系,提出了待建工程、已建工程推求和協(xié)調(diào)頻率設(shè)計值的思路。

    考慮位置參數(shù)、尺度參數(shù)、形狀參數(shù)的廣義可變模型GAMLSS為變化環(huán)境水文頻率計算提供了直接可靠的流程和方法。顧西輝等[16-17]將GAMLSS模型應(yīng)用于全國839個氣象站極端降水發(fā)生率及新疆塔里木河洪水發(fā)生率的非平穩(wěn)性研究。高潔[18-19]將GAMLSS模型應(yīng)用于雅礱江流域。研究發(fā)現(xiàn),全流域降水特性無顯著變化,但是流域中下游極端降水量增加,并伴隨著最長連續(xù)干期變長、連續(xù)濕期變短。流域中上游極端降水平穩(wěn)性較好,下游不同站點極端降水非平穩(wěn)性特征差異較大。小流域及支流洪水的平穩(wěn)性特征易受局地極端降水影響。雖然有建壩修庫等人類活動,但隨著支流入?yún)R,流域干流水文站洪峰流量回復(fù)平穩(wěn)狀態(tài)。

    關(guān)于水文時間序列,基于平穩(wěn)假設(shè)和考慮非平穩(wěn)性,均開展了大量的研究工作。根據(jù)傳統(tǒng)規(guī)范要求,我國水電水利工程行業(yè),水文設(shè)計選擇P-III分布作為頻率計算的基礎(chǔ)??紤]極值理論和變化環(huán)境,傳統(tǒng)頻率計算成果與GEV分布、GAMLSS模型成果的差異尚缺乏定量分析。本文以雅礱江流域為研究對象,通過GEV統(tǒng)計參數(shù)和GAMLSS最優(yōu)模擬,量化分析GEV分布、GAMLSS模型成果與傳統(tǒng)P-III適線成果差異。雅礱江干流河道天然落差約3 870 m,水能資源豐富,是我國規(guī)劃建設(shè)的十四大水電基地之一。根據(jù)河流水電規(guī)劃,雅礱江流域水電技術(shù)可開發(fā)量約2 900萬kW,主要集中在流域中下游河段,分布著眾多水電水利工程。因此,充分認(rèn)識流域洪水時變特征、洪水頻率計算成果偏差,對于流域內(nèi)眾多涉水工程設(shè)計運行具有重要意義。

    1 數(shù)據(jù)和方法

    1.1 研究對象

    雅礱江流域(96°52′E~102°48′E,26°32′N~33°58′N)是長江上游金沙江最大的一級支流,地處青藏高原東南部,位于金沙江和大渡河之間的狹長地帶。流域面積約13萬km2,谷嶺高差大[20],地形地貌以高山峽谷為主,基于90 m空間分辨率的DEM數(shù)據(jù),海拔高程位于966~5 888 m[20],地勢西北高東南低。本研究收集了流域內(nèi)8個水文站:甘孜、雅江、瀘寧、道孚、朱巴、濯桑、莊房(二)、孫水關(guān)長系列年最大洪峰流量系列。

    1.2 研究方法

    1.2.1洪水系列趨勢檢驗和突變分析

    針對洪水系列的單調(diào)變化和突變特征,分別進(jìn)行趨勢檢驗和突變點分析。

    MK檢驗作為廣泛應(yīng)用于水文氣象時間序列趨勢分析的非參數(shù)化方法,得到了國際氣象組織的推薦[21],被廣泛應(yīng)用。在突變點分析方面,分別采用Pettitt、AMOC和ECP三種檢驗方法。Pettitt法通過檢驗時間序列均值變化,確定序列躍變發(fā)生時間[22]。基于R語言平臺,changepoint程序包[23]的AMOC檢驗法是基于似然函數(shù)框架的突變點檢測,時間序列無需滿足正態(tài)分布假設(shè)。ECP程序包[24]是基于凝聚和分裂算法的聚類分析法,可用于時間序列的多突變點檢驗。

    1.2.2洪水系列非平穩(wěn)性分析

    Stasinopoulos和Rigby于2005年提出(半)參數(shù)回歸模型GAMLSS(基于位置參數(shù)、尺度參數(shù)和形狀參數(shù)的廣義可加模型),該模型[25]可模擬統(tǒng)計參數(shù)隨時間的變化模式。本文為了簡化模型結(jié)構(gòu),采用兩參數(shù)分布函數(shù),基于對數(shù)正態(tài)(LOGNO)、Gamma(GA)、Gumbel(GU)和Weibull(WEI)分布,擬定位置參數(shù)、尺度參數(shù)與時間的關(guān)系,通過AIC準(zhǔn)則選擇最優(yōu)分布函數(shù),計算水文序列的時變頻率設(shè)計值。

    1.2.3洪水系列適線

    我國傳統(tǒng)設(shè)計洪水頻率適線通常采用P-III型曲線,極值統(tǒng)計近年來也得到了廣泛關(guān)注[26]。GEV作為極值統(tǒng)計的一個重要分類,包含了Gumbel、Frechet和Weibull三種分布類型。

    H(x;μ,σ,ξ)=exp{-[1+ξ(x-μ)/σ]-1/ξ},[1+ξ(x-μ)]/σ>0

    (1)

    式中,μ為位置參數(shù);σ為尺度參數(shù);ξ為形狀參數(shù)。μ和σ反映洪水的量級和變幅。ξ決定了分布函數(shù)的尾部形態(tài)[7]。ξ>0的極值II型(Frechet)呈現(xiàn)厚尾分布,尾部較長且無上邊界。ξ<0極值的III型分布(Weibull)具有上端點μ-σ/ξ,呈現(xiàn)薄尾分布。當(dāng)統(tǒng)計變量x趨近于無窮,ξ=0的極值I型分布(Gumbel)為無上下邊界的薄尾分布。

    表1 各水文站趨勢檢驗和突變點分析

    2 成果分析

    2.1 時間序列特性

    根據(jù)MK趨勢檢驗,各水文站的年最大洪峰流量序列均無顯著變化趨勢(顯著度水平5%)。由于收集的洪峰流量系列長度僅在37年至65年,本研究檢測一個突變點,并根據(jù)變點位置將序列劃分為前后兩部分。經(jīng)Pettitt、AMOC(mean)和ECP突變檢驗,甘孜、濯桑、莊房(二)、孫水關(guān)4站通過顯著水平5%突變檢驗,雅江站滿足顯著水平10%突變檢驗。變化趨勢和突變時間點如表1所示。

    2.2 平穩(wěn)性分析

    對各水文站年最大洪峰流量系列進(jìn)行GAMLSS非平穩(wěn)性分析,在LOGNO、GA、GU和WEI四種分布中,以LOGNO分布占優(yōu)。雖然各水文站年最大洪峰流量的MK趨勢檢驗無顯著變化,但是與甘孜、濯桑、莊房、孫水關(guān)、雅江5站通過突變檢驗相協(xié)調(diào)。GAMLSS模型顯示,其中的甘孜、濯桑和孫水關(guān)3站的年最大洪峰流量呈現(xiàn)非平穩(wěn)特征,如圖1所示。

    (1)甘孜站位置參數(shù)表征的年最大洪峰流量均值隨時間單調(diào)減小,MK檢驗呈不顯著減小趨勢,1986年以后均值突變減小。

    (2)濯桑站尺度參數(shù)表征的不同重現(xiàn)期年最大流量離散程度逐漸減小,體現(xiàn)為常遇頻率設(shè)計值增加、稀遇頻率設(shè)計值減小,并逐漸穩(wěn)定。MK檢驗常遇頻率設(shè)計值呈不顯著增加。

    (3)孫水關(guān)站位置參數(shù)表征的年最大洪峰流量均值在1980年~1990年出現(xiàn)拐點,先減小后增大。MK檢驗常遇頻率設(shè)計值呈不顯著增加,在1985年以后突變增加。

    受限于觀測資料多為常遇頻率,MK趨勢檢驗和突變分析主要反映的是常遇頻率下洪水事件。GAMLSS成果可以體現(xiàn)不同重現(xiàn)期表征的常遇頻率和稀遇頻率設(shè)計值差異化的時變特征。年最大洪峰流量的非平穩(wěn)性分析(見表2),基本映射出漸變性和突變性檢驗中尚不顯著的發(fā)展趨勢,具有一定的前瞻性。

    時間序列趨勢檢驗無顯著變化,采用傳統(tǒng)的P-III適線和穩(wěn)態(tài)GEV分布擬合仍然合適。基于GAMLSS非平穩(wěn)性分析和考慮突變點的GEV分布,可反映變化環(huán)境下洪峰流量頻率設(shè)計值隨時間的變化。

    表2 各水文站GAMLSS非平穩(wěn)性分析成果

    2.3 頻率適線成果對比

    P-III分布目估適線,通常微調(diào)Cv和Cs兩個參數(shù),控制適線的傾斜程度和彎曲程度,使高水端上翹,傾向于頭部特大值,使稀遇頻率設(shè)計偏安全,得到比較保守的頻率設(shè)計值。GAMLSS兩參數(shù)分布由于參數(shù)較少,靈活性低,不能人為提高高水端的擬合優(yōu)度,可能出現(xiàn)低頻稀遇設(shè)計洪水洪峰流量小于實測頭部特大值,該方法調(diào)整裕度較少。GEV分布通過擬合位置參數(shù)、尺度參數(shù)和形狀參數(shù)計算頻率設(shè)計值。位置參數(shù)和尺度參數(shù)均與集水面積具有相關(guān)性,關(guān)于形狀參數(shù)的影響因素,學(xué)術(shù)界尚無定論。但是,形狀參數(shù)又對于GEV尾部形態(tài)具有關(guān)鍵作用,以正負(fù)為界,區(qū)別厚尾和薄尾特征,體現(xiàn)了分布極值的偏離程度和密集程度,對于研究極端事件具有重要意義。

    2.3.1平穩(wěn)序列適線成果對比

    年最大洪峰流量序列相對平穩(wěn)的雅江、瀘寧、朱巴、道孚、莊房5站經(jīng)三參數(shù)P-III目估適線、GEV分布擬合及兩參數(shù)GAMLSS模型優(yōu)選,除了道孚站在低頻稀遇洪水時GEV分布推算頻率設(shè)計值偏大,其他4站均為GEV分布和GAMLSS模型優(yōu)選成果較為接近,普遍小于P-III適線成果,且部分站點(如瀘寧、莊房)隨著重現(xiàn)期增大,稀遇洪水設(shè)計值的偏差增加。三種方法成果對比見圖2。

    圖2 平穩(wěn)序列GAMLSS、P-III、GEV 3種方法計算成果比較

    根據(jù)平穩(wěn)序列5個水文站的擬合成果,三種方法對于20年一遇以下洪水洪峰流量的擬合偏差低于10%,50年一遇至200年一遇洪水洪峰流量的擬合偏差≤15%,1 000年一遇以下洪水洪峰流量的擬合誤差基本不超過20%。對于5 000年、10 000年等稀遇特大洪水的擬合偏差在25%至30%之間(見圖2f)。

    2.3.2非平穩(wěn)序列適線成果對比

    通過GAMLSS模型推算年最大洪峰流量時間序列非平穩(wěn)的甘孜、濯桑和孫水關(guān)3站頻率設(shè)計值的時變特性。隨著重現(xiàn)期增大,頻率設(shè)計值變幅增大,體現(xiàn)了采用有限資料推算稀遇事件的外延性導(dǎo)致不確定性增加(見圖3a、3b、3c),可通過頻率設(shè)計值時變范圍與設(shè)計值的比例來反映。

    甘孜站1956年~2016年洪峰流量頻率設(shè)計值時變范圍與均值比例,對于不同重現(xiàn)期基本在14.8%~15.5%。濯桑站1966年~2016年洪峰流量頻率設(shè)計值隨著重新期增大,時變范圍比例從10年一遇洪水23%增至10 000年一遇洪水的57%。孫水關(guān)站時變范圍比例從10年一遇洪水45%降至10 000年一遇洪水的39%。

    濯桑站和孫水關(guān)站由于集水面積較小,實測年最大洪峰流量量級相對較小(均值分別為287 m3/s和610 m3/s);因此,頻率設(shè)計值時變范圍比例較大(23%~57%)。GAMLSS頻率設(shè)計值的時變特征是否具有尺度效應(yīng),小流域非平穩(wěn)性普遍較大,有待更多站點驗證。

    隨著重現(xiàn)期增大,甘孜站和孫水關(guān)站洪峰流量時變范圍與均值的比例基本穩(wěn)定。濯桑站洪峰流量時變范圍比例明顯增大,主要由于低頻稀遇洪水時變減小趨勢(見圖1b、3b),體現(xiàn)為10年一遇以上洪水隨著重現(xiàn)期增大增幅減小,稀遇洪水致災(zāi)風(fēng)險隨時間降低。

    P-III目估適線成果均在GAMSS測算區(qū)間內(nèi)。但是,GEV分布在低頻大值的估算上與P-III適線和GAMLSS模型成果(見圖3a和3c)產(chǎn)生偏差。GEV分布對于非平穩(wěn)系列的擬合較為敏感,需要考慮參數(shù)時變特性和尾部特點。

    圖3 非平穩(wěn)序列GAMLSS、P-III、GEV三種方法計算成果比較

    表3 各水文站突變前后及完整系列GEV參數(shù)

    2.4 考慮突變的GEV變化分析

    2.4.1GEV參數(shù)分析

    Yang[7]發(fā)現(xiàn)中國極端降水GEV分布的位置參數(shù)、尺度參數(shù)與集水面積相關(guān),長江流域極端降水的GEV形狀參數(shù)以負(fù)值居多。與該結(jié)論一致,流域內(nèi)8個水文站原始系列、突變前后系列的位置參數(shù)和尺度參數(shù)均與集水面積(見表3)呈現(xiàn)雙對數(shù)相關(guān)(見圖4)。本流域位于長江流域上游支流,年最大洪峰流量GEV分布也以形狀參數(shù)<0的薄尾分布為主。僅道孚站和孫水關(guān)站的形狀參數(shù)>0,呈現(xiàn)厚尾分布,由此造成兩站GEV擬合頻率設(shè)計值比P-III分布和GAMLSS上限值均偏大(見圖2c、3c)。GEV對于極大值的過分高估,是否由于形狀參數(shù)>0所致,值得研究。

    圖4 GEV統(tǒng)計參數(shù)與集水面積雙對數(shù)線性相關(guān)

    2.4.2突變前后GEV分析

    選取形狀參數(shù)<0的甘孜、雅江、濯桑、莊房4個站點,基于突變點分析,對比完整系列、突變前、突變后3種系列的GEV位置參數(shù)、尺度參數(shù)和形狀參數(shù),推算頻率設(shè)計值。

    (1)甘孜站和濯桑站分別在1986年突變后和1989年突變前,系列的GEV形狀參數(shù)>0,由此導(dǎo)致GEV擬合值偏離較大,分析無效。

    (2)雅江站在突變點1998年以后,各頻率設(shè)計洪水均有增加(見圖5a)。因此1998年以前設(shè)計的工程有必要復(fù)核現(xiàn)狀防洪安全。

    (3)莊房站在2006年以前,各頻率設(shè)計洪水略大于現(xiàn)狀水平,2006年以前的設(shè)計仍可滿足現(xiàn)狀相應(yīng)的防洪標(biāo)準(zhǔn)(見圖5b)。

    圖5 基于完整系列、突變前后系列推算設(shè)計洪水頻率值

    3 結(jié) 論

    基于MK趨勢檢驗,Pettitt、AMOC和ECP突變分析,研究雅礱江流域8個水文站年最大洪峰流量漸變和突變特征,基于GAMLSS模型分析流域洪水的平穩(wěn)性特征。按照平穩(wěn)序列和非平穩(wěn)序列,對各水文站年最大洪峰流量系列進(jìn)行P-III目估適線、GEV分布擬合和GAMLSS模型成果對比分析,探索GEV分布位置參數(shù)、尺度參數(shù)和形狀參數(shù)特點,結(jié)合突變點分析,研究流域防洪安全。研究發(fā)現(xiàn):

    (1)GAMLSS擬合以對數(shù)正態(tài)最優(yōu)分布居多。其中,甘孜站、濯桑站和孫水關(guān)站年最大洪峰流量時間序列非平穩(wěn)。GAMLSS模型成果與實測系列漸變、突變分析統(tǒng)籌考慮,體現(xiàn)了不同重現(xiàn)期設(shè)計值時變特征,揭示了實測點難以反映的低頻洪水特性。

    (2)比較了3參數(shù)P-III目估適線、GEV分布和兩參數(shù)GAMLSS模型成果,針對平穩(wěn)序列5個站點,傾向于頭部特大值的P-III目估適線成果普遍大于以穿過點據(jù)中央為擬合原則的GEV分布和GAMLSS模擬成果。3種途徑推算20年一遇頻率洪水偏差<10%,200年一遇頻率洪水偏差<15%,1 000年一遇頻率洪水偏差基本不超過20%。

    (3)非平穩(wěn)序列的3個站點,P-III適線成果基本在GAMLSS模型擬合時變范圍內(nèi)。GEV擬合非平穩(wěn)序列建議考慮參數(shù)的時變特性和尾部特點。

    (4)雅礱江流域年最大洪峰流量序列GEV分布形狀參數(shù)以<0的薄尾分布為主。呈現(xiàn)厚尾分布的道孚站和孫水關(guān)站GEV頻率設(shè)計值遠(yuǎn)大于其他模型成果,存在極大值被高估的現(xiàn)象。雅江站1998年以前設(shè)計的工程有必要復(fù)核現(xiàn)狀防洪安全。

    猜你喜歡
    洪峰流量時變水文站
    SL流量計在特殊河段的應(yīng)用——以河源水文站為例
    退耕還林工程對渭河洪峰流量的影響
    基于時變Copula的股票市場相關(guān)性分析
    智富時代(2017年4期)2017-04-27 17:08:47
    郭家屯水文站單斷沙關(guān)系分析
    佛岡縣潖江流域年洪峰流量P-Ⅲ分布參數(shù)估算
    草壩水文站兩種蒸發(fā)器對比分析
    煙氣輪機(jī)復(fù)合故障時變退化特征提取
    大南川流域設(shè)計洪峰流量計算分析
    某特小流域設(shè)計洪峰流量計算分析
    韓府灣水文站報汛方案
    亚洲第一青青草原| 久久国产亚洲av麻豆专区| 色综合婷婷激情| 国产视频一区二区在线看| tube8黄色片| 欧美激情极品国产一区二区三区| 久久精品国产亚洲av高清一级| 国产成人av教育| 黄片播放在线免费| 欧美亚洲日本最大视频资源| 人妻 亚洲 视频| 满18在线观看网站| 成人三级做爰电影| 国产精品99久久99久久久不卡| 丁香欧美五月| 欧美日韩亚洲综合一区二区三区_| 一进一出抽搐动态| 国产1区2区3区精品| 国产精品久久久久久精品古装| 精品人妻在线不人妻| 日韩欧美免费精品| 色综合欧美亚洲国产小说| 午夜福利影视在线免费观看| 无遮挡黄片免费观看| 精品久久久久久久毛片微露脸| 成人黄色视频免费在线看| 丰满迷人的少妇在线观看| 又紧又爽又黄一区二区| 日韩欧美三级三区| 亚洲精品自拍成人| 在线播放国产精品三级| 一级a爱视频在线免费观看| 午夜福利欧美成人| 人成视频在线观看免费观看| 日韩欧美在线二视频 | 亚洲精品美女久久久久99蜜臀| 国产黄色免费在线视频| 黄片播放在线免费| 18禁观看日本| 日韩熟女老妇一区二区性免费视频| 欧美日韩视频精品一区| 少妇猛男粗大的猛烈进出视频| 久久婷婷成人综合色麻豆| 99精品欧美一区二区三区四区| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 三上悠亚av全集在线观看| 久久久久精品人妻al黑| bbb黄色大片| 黄色毛片三级朝国网站| 伊人久久大香线蕉亚洲五| 天堂动漫精品| 日韩有码中文字幕| 在线观看舔阴道视频| 国产高清视频在线播放一区| 天天影视国产精品| 91字幕亚洲| 91字幕亚洲| 国产不卡av网站在线观看| 国产成人精品在线电影| 国产精品亚洲av一区麻豆| 午夜福利影视在线免费观看| 国产又色又爽无遮挡免费看| 亚洲国产欧美一区二区综合| 校园春色视频在线观看| 曰老女人黄片| 国产欧美日韩一区二区三| 久久精品成人免费网站| 人妻丰满熟妇av一区二区三区 | 免费观看人在逋| 在线av久久热| 热99国产精品久久久久久7| 老熟妇仑乱视频hdxx| 啦啦啦 在线观看视频| 久久香蕉激情| 亚洲午夜理论影院| 亚洲成人免费av在线播放| 国产精品av久久久久免费| e午夜精品久久久久久久| 欧美日韩国产mv在线观看视频| 99久久人妻综合| 老司机影院毛片| 国产区一区二久久| 欧美精品亚洲一区二区| 久久精品国产亚洲av香蕉五月 | 999精品在线视频| 一二三四社区在线视频社区8| 一区在线观看完整版| 亚洲 欧美一区二区三区| 欧美日韩亚洲高清精品| 99精国产麻豆久久婷婷| 女警被强在线播放| 免费黄频网站在线观看国产| 在线免费观看的www视频| av有码第一页| 老熟妇乱子伦视频在线观看| 欧美 日韩 精品 国产| 美女高潮喷水抽搐中文字幕| 中亚洲国语对白在线视频| 欧美丝袜亚洲另类 | 人妻丰满熟妇av一区二区三区 | 深夜精品福利| 免费观看人在逋| 最新在线观看一区二区三区| 欧美最黄视频在线播放免费 | 脱女人内裤的视频| 999久久久精品免费观看国产| 黄片小视频在线播放| 高清黄色对白视频在线免费看| 一级作爱视频免费观看| 午夜久久久在线观看| 又黄又爽又免费观看的视频| 国产亚洲精品第一综合不卡| 久久香蕉国产精品| 久久精品91无色码中文字幕| av福利片在线| 久久亚洲真实| 午夜老司机福利片| 亚洲自偷自拍图片 自拍| 啦啦啦免费观看视频1| 交换朋友夫妻互换小说| xxx96com| 免费在线观看影片大全网站| 久久99一区二区三区| 亚洲色图综合在线观看| 波多野结衣av一区二区av| 高清黄色对白视频在线免费看| 伦理电影免费视频| 99久久人妻综合| 午夜两性在线视频| 欧美+亚洲+日韩+国产| 极品教师在线免费播放| 亚洲综合色网址| 黄色成人免费大全| 成人亚洲精品一区在线观看| 久久久国产一区二区| 国产视频一区二区在线看| 日本黄色日本黄色录像| 高清av免费在线| 一区二区三区激情视频| 丝袜美腿诱惑在线| 黄色丝袜av网址大全| 99在线人妻在线中文字幕 | 黑人巨大精品欧美一区二区mp4| 一本一本久久a久久精品综合妖精| 国产精品美女特级片免费视频播放器 | 国产极品粉嫩免费观看在线| 美国免费a级毛片| 午夜福利乱码中文字幕| x7x7x7水蜜桃| a级毛片黄视频| 18禁国产床啪视频网站| 丝袜美足系列| 啦啦啦免费观看视频1| 嫁个100分男人电影在线观看| av国产精品久久久久影院| 久久婷婷成人综合色麻豆| 天天添夜夜摸| 欧美另类亚洲清纯唯美| 欧美亚洲日本最大视频资源| 国产精品久久久久久人妻精品电影| 久久人人爽av亚洲精品天堂| 天天影视国产精品| 久久性视频一级片| 757午夜福利合集在线观看| 亚洲,欧美精品.| 久久这里只有精品19| 欧美性长视频在线观看| 亚洲精品美女久久久久99蜜臀| 免费不卡黄色视频| 成人特级黄色片久久久久久久| 中文字幕最新亚洲高清| 中文字幕人妻熟女乱码| 久久草成人影院| 黄色怎么调成土黄色| 色婷婷久久久亚洲欧美| 十八禁高潮呻吟视频| 亚洲精品粉嫩美女一区| 嫁个100分男人电影在线观看| 免费av中文字幕在线| 久久久久国产一级毛片高清牌| 黄片播放在线免费| 亚洲专区国产一区二区| 精品视频人人做人人爽| 免费在线观看亚洲国产| 黄网站色视频无遮挡免费观看| 岛国毛片在线播放| 夫妻午夜视频| 亚洲伊人色综图| 精品福利永久在线观看| 精品卡一卡二卡四卡免费| 久久亚洲精品不卡| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 免费女性裸体啪啪无遮挡网站| 成人18禁在线播放| 人人妻人人澡人人爽人人夜夜| 亚洲人成电影观看| 亚洲自偷自拍图片 自拍| 亚洲男人天堂网一区| 精品无人区乱码1区二区| 欧美不卡视频在线免费观看 | 成在线人永久免费视频| av天堂久久9| 法律面前人人平等表现在哪些方面| 在线观看66精品国产| 久久久水蜜桃国产精品网| 精品欧美一区二区三区在线| 自线自在国产av| 国产成人欧美在线观看 | 美女 人体艺术 gogo| 欧美不卡视频在线免费观看 | 老司机午夜十八禁免费视频| 热re99久久精品国产66热6| 俄罗斯特黄特色一大片| 91大片在线观看| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看| 国产精品久久久久成人av| 精品乱码久久久久久99久播| 高清av免费在线| 天天躁夜夜躁狠狠躁躁| 大香蕉久久成人网| 亚洲一区中文字幕在线| 黑丝袜美女国产一区| 久久国产精品人妻蜜桃| 日韩欧美一区视频在线观看| 男男h啪啪无遮挡| 自线自在国产av| 国产黄色免费在线视频| 色老头精品视频在线观看| 亚洲成人国产一区在线观看| 精品久久蜜臀av无| 韩国精品一区二区三区| 久9热在线精品视频| 成年版毛片免费区| 国产有黄有色有爽视频| 国产成人精品久久二区二区免费| 一二三四社区在线视频社区8| 久久久国产一区二区| 自线自在国产av| 一级片免费观看大全| 精品午夜福利视频在线观看一区| 黑人巨大精品欧美一区二区蜜桃| 成人手机av| 在线十欧美十亚洲十日本专区| 亚洲熟妇熟女久久| 丰满迷人的少妇在线观看| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 777久久人妻少妇嫩草av网站| av有码第一页| 脱女人内裤的视频| 深夜精品福利| 国产aⅴ精品一区二区三区波| 亚洲国产精品sss在线观看 | √禁漫天堂资源中文www| 99国产极品粉嫩在线观看| 99riav亚洲国产免费| 日韩免费高清中文字幕av| 亚洲久久久国产精品| 一级片'在线观看视频| 18禁观看日本| 久久香蕉激情| 国产一区二区三区视频了| 久久99一区二区三区| 看黄色毛片网站| 深夜精品福利| 超碰97精品在线观看| 色94色欧美一区二区| 制服人妻中文乱码| 国产又爽黄色视频| 狂野欧美激情性xxxx| 18禁观看日本| 在线观看日韩欧美| 激情视频va一区二区三区| 男女免费视频国产| 香蕉久久夜色| 女性生殖器流出的白浆| 丁香六月欧美| 日本黄色日本黄色录像| 我的亚洲天堂| 亚洲,欧美精品.| 国产精品免费大片| 搡老岳熟女国产| 老司机亚洲免费影院| 日韩熟女老妇一区二区性免费视频| 欧美久久黑人一区二区| 日韩中文字幕欧美一区二区| 在线观看免费视频网站a站| 很黄的视频免费| 亚洲精品国产区一区二| 十八禁人妻一区二区| 亚洲va日本ⅴa欧美va伊人久久| 色94色欧美一区二区| 精品少妇久久久久久888优播| 久热这里只有精品99| 首页视频小说图片口味搜索| 国产在线一区二区三区精| 如日韩欧美国产精品一区二区三区| 欧美日韩精品网址| 久久中文字幕一级| 韩国av一区二区三区四区| 久久久水蜜桃国产精品网| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 精品卡一卡二卡四卡免费| 国产精品二区激情视频| 黑人欧美特级aaaaaa片| 一本一本久久a久久精品综合妖精| 久久久久视频综合| 久久香蕉精品热| 视频区欧美日本亚洲| 三级毛片av免费| 香蕉久久夜色| 国产精品久久久人人做人人爽| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 亚洲人成电影观看| 黄网站色视频无遮挡免费观看| 日韩人妻精品一区2区三区| 欧美 亚洲 国产 日韩一| 国产精品一区二区在线观看99| 男人舔女人的私密视频| 国产成人影院久久av| 激情在线观看视频在线高清 | 水蜜桃什么品种好| 天天操日日干夜夜撸| 国产成+人综合+亚洲专区| 午夜福利免费观看在线| 男男h啪啪无遮挡| 19禁男女啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 国产精品免费大片| 国产av一区二区精品久久| 成熟少妇高潮喷水视频| avwww免费| 中文字幕色久视频| 99re在线观看精品视频| 1024视频免费在线观看| 亚洲avbb在线观看| avwww免费| 成熟少妇高潮喷水视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧洲精品一区二区精品久久久| 国产主播在线观看一区二区| 丰满人妻熟妇乱又伦精品不卡| av欧美777| 欧美日韩亚洲国产一区二区在线观看 | 久久久久精品人妻al黑| 看免费av毛片| 久久久水蜜桃国产精品网| 高潮久久久久久久久久久不卡| 久久精品国产清高在天天线| 制服人妻中文乱码| 99re6热这里在线精品视频| 在线看a的网站| 他把我摸到了高潮在线观看| 免费久久久久久久精品成人欧美视频| 国产av一区二区精品久久| 黑人猛操日本美女一级片| 午夜福利欧美成人| 国产精品欧美亚洲77777| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 日本撒尿小便嘘嘘汇集6| av天堂在线播放| 日本vs欧美在线观看视频| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 亚洲精品在线观看二区| 国产亚洲欧美98| 极品少妇高潮喷水抽搐| e午夜精品久久久久久久| 女人高潮潮喷娇喘18禁视频| 中出人妻视频一区二区| avwww免费| 多毛熟女@视频| 国产蜜桃级精品一区二区三区 | 亚洲久久久国产精品| 丝袜在线中文字幕| 在线看a的网站| 免费日韩欧美在线观看| 99国产精品一区二区蜜桃av | 欧美精品一区二区免费开放| 国产成人欧美| av福利片在线| 极品少妇高潮喷水抽搐| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利欧美成人| 国产片内射在线| 欧美老熟妇乱子伦牲交| 夜夜躁狠狠躁天天躁| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| 一边摸一边抽搐一进一出视频| 精品免费久久久久久久清纯 | a在线观看视频网站| 丰满饥渴人妻一区二区三| 国产精品偷伦视频观看了| 亚洲精品国产区一区二| 美女 人体艺术 gogo| 免费看a级黄色片| 成人国语在线视频| a在线观看视频网站| 亚洲欧美精品综合一区二区三区| 又大又爽又粗| 欧美 日韩 精品 国产| 99热国产这里只有精品6| 亚洲伊人色综图| 女人被狂操c到高潮| 在线观看一区二区三区激情| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一二三| 黄色a级毛片大全视频| 亚洲熟妇中文字幕五十中出 | 韩国精品一区二区三区| 久久草成人影院| 午夜福利乱码中文字幕| 精品久久久久久久久久免费视频 | 黄频高清免费视频| 国产精品99久久99久久久不卡| 18禁裸乳无遮挡免费网站照片 | 视频在线观看一区二区三区| 婷婷成人精品国产| 久久久精品国产亚洲av高清涩受| 天天影视国产精品| 亚洲av第一区精品v没综合| 香蕉久久夜色| 亚洲五月天丁香| 午夜两性在线视频| 日本一区二区免费在线视频| 在线视频色国产色| 国产精品欧美亚洲77777| 亚洲精品一二三| 精品国内亚洲2022精品成人 | 制服人妻中文乱码| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 日韩人妻精品一区2区三区| 亚洲中文字幕日韩| 无人区码免费观看不卡| 精品久久久久久,| 成熟少妇高潮喷水视频| 午夜福利,免费看| 日日夜夜操网爽| 国产精品av久久久久免费| 国产精品国产高清国产av | 男人舔女人的私密视频| 精品欧美一区二区三区在线| 亚洲第一青青草原| 热re99久久国产66热| 51午夜福利影视在线观看| 在线观看66精品国产| 成人黄色视频免费在线看| 国产在线一区二区三区精| 久久精品国产99精品国产亚洲性色 | 欧美精品亚洲一区二区| 国产av又大| 99国产极品粉嫩在线观看| av网站在线播放免费| 一进一出好大好爽视频| 真人做人爱边吃奶动态| 午夜免费鲁丝| 久久久久久久午夜电影 | 巨乳人妻的诱惑在线观看| 在线播放国产精品三级| 久久香蕉国产精品| 欧美国产精品va在线观看不卡| 天堂中文最新版在线下载| 精品免费久久久久久久清纯 | 国产精品.久久久| 久久热在线av| 国产精品国产av在线观看| 欧美精品av麻豆av| 欧美日韩中文字幕国产精品一区二区三区 | а√天堂www在线а√下载 | 欧美黄色淫秽网站| 国产精品久久久人人做人人爽| 脱女人内裤的视频| 天堂中文最新版在线下载| 国产视频一区二区在线看| 18禁黄网站禁片午夜丰满| 啦啦啦免费观看视频1| 最新美女视频免费是黄的| videosex国产| 又黄又爽又免费观看的视频| 亚洲色图 男人天堂 中文字幕| 国产在线一区二区三区精| 精品午夜福利视频在线观看一区| 99精品欧美一区二区三区四区| 在线观看舔阴道视频| 一级,二级,三级黄色视频| 国产欧美日韩一区二区三区在线| 成人三级做爰电影| 捣出白浆h1v1| 亚洲国产精品sss在线观看 | 怎么达到女性高潮| 国产亚洲精品久久久久5区| av免费在线观看网站| 一级作爱视频免费观看| 巨乳人妻的诱惑在线观看| 久久久久久久久免费视频了| 黄色视频,在线免费观看| 水蜜桃什么品种好| bbb黄色大片| av欧美777| 国产激情欧美一区二区| 欧美在线一区亚洲| 一进一出抽搐gif免费好疼 | 欧美日韩精品网址| 国产欧美日韩一区二区精品| 巨乳人妻的诱惑在线观看| 我的亚洲天堂| 亚洲av第一区精品v没综合| 激情视频va一区二区三区| 国产成人免费观看mmmm| 黄片小视频在线播放| 国产成人精品久久二区二区免费| 色精品久久人妻99蜜桃| 久久久久久久国产电影| 国产一区二区激情短视频| 黑人猛操日本美女一级片| 欧美日本中文国产一区发布| 成年人午夜在线观看视频| 欧美成狂野欧美在线观看| 男女午夜视频在线观看| 精品少妇一区二区三区视频日本电影| 久久性视频一级片| 亚洲国产精品sss在线观看 | 午夜91福利影院| 国产精品 国内视频| 999久久久国产精品视频| 美女扒开内裤让男人捅视频| 女人久久www免费人成看片| 久久久久国产精品人妻aⅴ院 | 天天影视国产精品| 精品国内亚洲2022精品成人 | 搡老熟女国产l中国老女人| a在线观看视频网站| 男女午夜视频在线观看| 久久久久久久久免费视频了| 天堂中文最新版在线下载| 999精品在线视频| 老司机福利观看| 亚洲精品国产精品久久久不卡| 成年版毛片免费区| 日本wwww免费看| 两个人免费观看高清视频| 啦啦啦视频在线资源免费观看| 操出白浆在线播放| 一级a爱片免费观看的视频| 人人澡人人妻人| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| 久久久国产一区二区| 国产精品.久久久| 久久国产精品影院| 99热网站在线观看| 一本一本久久a久久精品综合妖精| 国产免费男女视频| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 亚洲国产欧美网| 在线观看午夜福利视频| 久久久国产欧美日韩av| 人人澡人人妻人| 亚洲一码二码三码区别大吗| 精品电影一区二区在线| 国产高清videossex| 香蕉丝袜av| 12—13女人毛片做爰片一| 夫妻午夜视频| 丝袜美腿诱惑在线| 美女视频免费永久观看网站| www.精华液| 亚洲性夜色夜夜综合| 午夜两性在线视频| 777米奇影视久久| 亚洲精品粉嫩美女一区| 变态另类成人亚洲欧美熟女 | 黄色片一级片一级黄色片| 久久精品国产亚洲av香蕉五月 | 在线观看午夜福利视频| 在线视频色国产色| 日日摸夜夜添夜夜添小说| 最近最新免费中文字幕在线| 亚洲五月色婷婷综合| 成人精品一区二区免费| 中国美女看黄片| 精品久久久久久,| 窝窝影院91人妻| 久久精品成人免费网站| 欧美精品高潮呻吟av久久| 免费观看精品视频网站| av超薄肉色丝袜交足视频| 91成人精品电影| 99精品在免费线老司机午夜| 亚洲精华国产精华精| videos熟女内射| 亚洲色图综合在线观看| 又大又爽又粗| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区免费| 久久九九热精品免费| 日韩 欧美 亚洲 中文字幕| 亚洲av成人一区二区三| tocl精华| 国产精品久久久人人做人人爽| 亚洲自偷自拍图片 自拍| 日本一区二区免费在线视频| 日韩欧美三级三区| 成年版毛片免费区|