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

    基于Copula 函數(shù)風(fēng)險(xiǎn)組合的西峽縣山洪災(zāi)害預(yù)警分析

    2022-11-28 06:38:12張玉昶榮富強(qiáng)王明閃王建國(guó)
    中國(guó)農(nóng)村水利水電 2022年11期
    關(guān)鍵詞:山洪模數(shù)雨量

    張玉昶,榮富強(qiáng),王明閃,王建國(guó),萬(wàn) 芳

    (1.河南省地質(zhì)礦產(chǎn)勘查開(kāi)發(fā)局第二地質(zhì)環(huán)境調(diào)查院,河南鄭州 450000;2.華北水利水電大學(xué)水資源學(xué)院,河南鄭州 450011)

    近年來(lái),極端氣候頻發(fā),導(dǎo)致山洪災(zāi)害越來(lái)越頻繁發(fā)生,山洪災(zāi)害已經(jīng)成為世界各國(guó)自然災(zāi)害中的主要災(zāi)種之一,造成了大量的人員傷亡和社會(huì)經(jīng)濟(jì)損失[1],為了防御山洪災(zāi)害,減輕災(zāi)害損失,國(guó)內(nèi)外許多專(zhuān)家學(xué)者對(duì)山洪災(zāi)害預(yù)警進(jìn)行了研究。王協(xié)康等[2]采用TRI-GRS 模型分析了小流域的滑坡易發(fā)性,將暴雨山洪洪水劃分為預(yù)警區(qū)、中、高風(fēng)險(xiǎn)區(qū)。我國(guó)下墊面等條件復(fù)雜多樣,山洪災(zāi)害成災(zāi)機(jī)理復(fù)雜,而臨界雨量是山洪災(zāi)害的主要影響因素。目前,臨界雨量的計(jì)算方法主要有數(shù)據(jù)驅(qū)動(dòng)法及水文水力學(xué)法[3,4]。但是,臨界雨量受前期影響雨量影響較大[5],江錦紅等[6]采用雙曲函數(shù)確定的暴雨臨界雨量最為山洪災(zāi)害的預(yù)警標(biāo)準(zhǔn);劉媛媛等[7]針對(duì)小流域,利用臨界流量及水位計(jì)算相應(yīng)的臨界頻率;張珊珊等[8]基于HEC-HMS 模型,研究選用雨量預(yù)警指標(biāo),綜合考慮土壤含水量、匯流時(shí)間等,反推臨界雨量。臨界雨量的研究成果豐碩,但是未考慮到幾種風(fēng)險(xiǎn)組合情況的臨界雨量及相應(yīng)的預(yù)警分析。水文科學(xué)領(lǐng)域的種隨機(jī)變量之間往往呈現(xiàn)出各種復(fù)雜的線(xiàn)性、非線(xiàn)性的相關(guān)關(guān)系,本文引入Copula 函數(shù)來(lái)描述降雨量與峰值雨強(qiáng)變量間的相關(guān)性,Copula 函數(shù)可將聯(lián)合分布的函數(shù)和它各自的邊緣分布函數(shù)聯(lián)系在一起,基于變量之間的非線(xiàn)性相關(guān)關(guān)系而建立的,可以描述變量間非線(xiàn)性、非對(duì)稱(chēng)和對(duì)稱(chēng)的相關(guān)關(guān)系,目前在水文科學(xué)中得到廣泛的應(yīng)用[9-11]。

    1 研究區(qū)概況

    西峽縣位于河南省西南邊陲,與陜西省接界,地處秦嶺山系東南余脈伏牛山南。東南與本市的內(nèi)鄉(xiāng)、淅川相連,北西與欒川、盧氏和陜西省的商南縣毗鄰。地處北亞熱帶向暖溫帶的過(guò)度地區(qū),屬于亞熱帶季風(fēng)型大陸性氣候,季風(fēng)影響明顯,四季分明,春季溫度升降劇烈,雨水分布不均,夏季熱雨集中,暴雨頻發(fā),地質(zhì)地貌復(fù)雜,加上受人類(lèi)活動(dòng)的影響,導(dǎo)致山洪災(zāi)害發(fā)生頻繁。山洪災(zāi)害在縣境內(nèi)分布廣,具有來(lái)勢(shì)猛、流速快、破壞力大、突發(fā)性強(qiáng)的特點(diǎn),預(yù)報(bào)、預(yù)測(cè)、預(yù)防難度較大,不僅對(duì)山區(qū)的基礎(chǔ)設(shè)施造成毀滅性破壞,而且對(duì)人民群眾的生命安全構(gòu)成極大的威脅,已經(jīng)成為當(dāng)前防災(zāi)減災(zāi)工作中的突出問(wèn)題,是山區(qū)經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的重要制約因素之一。

    西峽縣山洪災(zāi)害致災(zāi)因素具有自然地理和經(jīng)濟(jì)社會(huì)的雙重屬性,降雨因素是誘發(fā)山洪災(zāi)害的直接因素和激發(fā)條件。降雨量大,多數(shù)情況下意味著雨強(qiáng)高、激發(fā)力強(qiáng),在一定的下墊面條件下,易產(chǎn)生溪河洪水災(zāi)害、泥石流災(zāi)害和滑坡災(zāi)害。根據(jù)西峽縣近年來(lái)的統(tǒng)計(jì)資料表明,每次山洪災(zāi)害的發(fā)生,都是在降雨量較大這一基本前提條件下形成的;高強(qiáng)度的降雨是引起山洪災(zāi)害的最主要原因之一。

    在相同的條件下,降雨歷時(shí)越長(zhǎng),降雨量越大,產(chǎn)生的徑流量越大,山洪災(zāi)害損失也越嚴(yán)重。在溪河洪水、泥石流、滑坡三種山洪災(zāi)害類(lèi)型中,溪河洪水災(zāi)害和泥石流災(zāi)害特點(diǎn)相似,在下墊面條件滿(mǎn)足的情況下,只要有足夠大的降雨強(qiáng)度和降雨量,溪河洪水和泥石流災(zāi)害就可發(fā)生,降雨歷時(shí)越長(zhǎng),所產(chǎn)生的洪量越大,災(zāi)害損失越大。滑坡災(zāi)害與降雨歷時(shí)的關(guān)系更為明顯,一般而言,滑坡和降雨并不是同時(shí)發(fā)生,滑坡一般滯后于降雨,由溪河洪水誘發(fā)的滑坡災(zāi)害會(huì)隨著降雨歷時(shí)的增長(zhǎng),洪流對(duì)坡腳的掏蝕,滑坡也會(huì)越來(lái)越多。

    2 混合Copula 函數(shù)的風(fēng)險(xiǎn)組合概率的臨界雨量確定

    2.1 Copula理論及Copula函數(shù)分布模型

    1959年Copula函數(shù)理論提出,Sklar提出可把一個(gè)聯(lián)合分布函數(shù)分解為一個(gè)Copula 函數(shù)和K個(gè)邊緣分布函數(shù),而分解出來(lái)的Copula 函數(shù)可用來(lái)描述變量之間相關(guān)性,也就是可將聯(lián)合分布的函數(shù)和它各自的邊緣分布函數(shù)聯(lián)系在一起的函數(shù)。N元Copula函數(shù)C(u1,u2,…,uN)具有以下性質(zhì):

    常用的Copula 數(shù)有很多種,本文主要介紹常用的多元正態(tài)Copula 函數(shù)、阿基米德Copula 函數(shù)Archimedean 的三類(lèi)常用函數(shù)(Clayton 函數(shù)、Gumbel Copula 函數(shù)、Frank Copula 函數(shù))和極值Copula函數(shù)。

    (1)多元正態(tài)Copula函數(shù)。其中N元正態(tài)Copula函數(shù)的分布函數(shù)與密度函數(shù)可表達(dá)為:

    式中:ρ為對(duì)角線(xiàn)上的元素為1的對(duì)稱(chēng)正定矩陣;|ρ|表示與矩陣ρ相對(duì)應(yīng)的行列式的值;Φρ(?)表示相關(guān)系數(shù)矩陣為ρ的標(biāo)準(zhǔn)正態(tài)分布函數(shù);Φ-1(?)為標(biāo)準(zhǔn)多元正態(tài)函數(shù)Φ(?)的逆函數(shù);?=(?1,?2,…,?N)′,?n=Φ-1(un),n=1,2,…,N,I為單位矩陣。

    (2)阿基米德Copula函數(shù)。其中阿基米德Copula函數(shù)可表示為:

    Gumbel Copula 函數(shù)、Clayton Copula 函數(shù)、Frank Copula 函數(shù)是常用的二元阿基米德Copula 函數(shù),分別由它們擴(kuò)展到N元阿基米德Copula函數(shù)可表示為:

    ①N元Gumbel Copula函數(shù)的表達(dá)式為:

    ②N元Clayton Copula函數(shù)的表達(dá)式為:

    ③N元Frank Copula函數(shù)的表達(dá)式為:

    其中,α,θ,λ為相關(guān)參數(shù)。

    (3)極值Copula 函數(shù)。極值Copula 函數(shù)(Extreme Value Copula,EVC)可表達(dá)為:

    2.2 混合Copula函數(shù)的構(gòu)造與相關(guān)性分析

    (1)不同類(lèi)型的Copula 函數(shù)比較分析。一般情況下,多元正態(tài)Copula 函數(shù)通常描述變量之間的相關(guān)關(guān)系,但是由于此函數(shù)有對(duì)稱(chēng)性的特點(diǎn),所以對(duì)于變量的非對(duì)稱(chēng)相關(guān)關(guān)系很難擬合。

    而Gumbel Copula 函數(shù)和Clayton Copula 函數(shù)都具有非對(duì)稱(chēng)性的特點(diǎn),其中Gumbel Copula函數(shù)為“J”字形分布,下尾低而上尾高,對(duì)水文變量的下尾部變化不太敏感,而對(duì)上尾部的分布變化比較敏感,所以很難描述下尾部的相關(guān)變化情況,即當(dāng)一個(gè)水文變量出現(xiàn)極大值時(shí),另外兩個(gè)水文變量也出現(xiàn)極大值的概率增大;

    Clayton Copula 函數(shù)為“L”字形分布,下尾高而上尾低,對(duì)水文變量的下尾部變化比較敏感,而對(duì)上尾部的分布變化不太敏感,所以很難描述上尾部的相關(guān)變化情況,即當(dāng)一個(gè)水文變量出現(xiàn)極小值時(shí),其他4個(gè)水文變量也出現(xiàn)極小值的概率增大。

    Frank Copula 函數(shù)為“U”字形分布,它有對(duì)稱(chēng)性的特點(diǎn),所以很難描述水文變量之間的非對(duì)稱(chēng)關(guān)系,F(xiàn)rank Copula 函數(shù)只適合描述具有對(duì)稱(chēng)相關(guān)結(jié)構(gòu)的變量之間的相關(guān)關(guān)系,即各個(gè)水文變量間極大值相關(guān)性與極小值相關(guān)性是對(duì)稱(chēng)增長(zhǎng)的,但是其尾部的分布變量是比較獨(dú)立的,所以Frank Copula 函數(shù)無(wú)論在描述上尾部還是下尾部的相關(guān)性中都是不敏感的,故無(wú)法描述尾部變化的相關(guān)性。

    (2)混合Copula 函數(shù)的構(gòu)造與相關(guān)性分析。Gumbel Copu‐la 函數(shù)、Clayton Copula 函數(shù)、Frank Copula 函數(shù)三類(lèi)常用的阿基米德函數(shù)能夠捕捉尾部相關(guān)的情形:上尾部的相關(guān)、下尾部的相關(guān)、上尾部和下尾部的對(duì)稱(chēng)相關(guān)。這些Copula 函數(shù)具有描述各種水文變量模式之間的關(guān)系,特別是尾部相關(guān)關(guān)系,水文系統(tǒng)中各變量之間的關(guān)系是復(fù)雜多變的,不是拘泥于某種特定關(guān)系,它們只能反映水文變量間相關(guān)性的某個(gè)側(cè)面,因此很難用一個(gè)簡(jiǎn)單的Copula 函數(shù)來(lái)全面的刻畫(huà)水文系統(tǒng)中各變量之間的相關(guān)模式,所以一個(gè)更加靈活的Copula 函數(shù)需要構(gòu)造,來(lái)描述各種水文變量模式之間的關(guān)系。應(yīng)用不同Copula 函數(shù)的優(yōu)點(diǎn),本文選用Gumbel Copula 函數(shù)、Clayton Copula 函數(shù)、Frank Copula 函數(shù)的線(xiàn)性組合來(lái)構(gòu)造混合的Copula 函數(shù),可以更加靈活的刻畫(huà)水文系統(tǒng)中各變量之間的相關(guān)關(guān)系。

    混合Copula函數(shù)M-Copula表達(dá)式為:

    式中:CG,CF,Ccl分別為Gumbel Copula、Frank Copula、Clayton Copula 函數(shù);wG,wF,wcl為相應(yīng)Copula 函數(shù)的權(quán)重系數(shù)。由式(8)~(10)可知,MC3中包括了6 個(gè)參數(shù),參數(shù)向量(α,λ,θ)用來(lái)描述變量間相關(guān)的程度,變量間相關(guān)的模式由線(xiàn)性權(quán)重參數(shù)向量(wG,wF,wcl)表示。

    2.3 Copula模型的參數(shù)估計(jì)、檢驗(yàn)與評(píng)價(jià)

    (1)Copula 模型的參數(shù)估計(jì)。Copula 函數(shù)模型中的參數(shù)估計(jì)有很多種,一般采用極大似然和矩估計(jì),而其中極大似然估計(jì)是較常用的Copula 模型參數(shù)估計(jì)方法,其聯(lián)合分布函數(shù)的密度函數(shù)為:

    式中:θc為Copula 函數(shù)的1 ×mc維參數(shù)向量;Fn(xn;θn)為邊緣分布函數(shù);θn為邊緣分布函數(shù)Fn(xn;θn)的1 ×mn維參數(shù)向量;θ=(θ1,θ2,…,θN,θc)′;n=1,2,…,N。

    因此,可以得到樣本(x1t,x2t,…,xNt),t=1,2,…,T的對(duì)數(shù)似然函數(shù)為:

    使似然函數(shù)取的最大值的θ即是最大似然估計(jì)值。

    (2)Copula 模型的檢驗(yàn)與評(píng)價(jià)。應(yīng)用的Copula 分布函數(shù)是否能夠很好的擬合變量之間的相關(guān)結(jié)構(gòu)及分布,所以Copula 函數(shù)的檢驗(yàn)與擬合優(yōu)度評(píng)價(jià)需要建立。K-S用于檢驗(yàn)樣本是否服從同一分布,故應(yīng)用其對(duì)Copula 分布函數(shù)進(jìn)行檢驗(yàn);Copula 函數(shù)的擬合度評(píng)價(jià)采用均方根誤差RSME最小準(zhǔn)則來(lái)計(jì)算,其定義如式(14):

    式中:N是樣本容量;i為樣本序號(hào);pc為模型計(jì)算的理論頻率;p0為聯(lián)合分布的經(jīng)驗(yàn)頻率。

    其中K-S檢驗(yàn)的統(tǒng)計(jì)量D計(jì)算公式如(15)式所示:

    式中:N為樣本容量;Ck為樣本xk=(x1k,x2k,x3k)的Copula 值;uk為樣本中滿(mǎn)足條件x≤xk的個(gè)數(shù),即滿(mǎn)足:x1≤x1k,x2≤x2k,x3≤x3k。

    2.4 風(fēng)險(xiǎn)組合概率的臨界雨量確定

    設(shè)隨機(jī)變量X,Y分別服從邊緣分布FX(x)、FY(y),其聯(lián)合分布函數(shù)為F(x,y),其中X表示總雨量,Y表示峰值雨強(qiáng),則存在一個(gè)Copula函數(shù):

    臨界雨量中總雨量和峰值雨強(qiáng)存在多種耦合情況,可以運(yùn)用上述的聯(lián)合概率分布模型計(jì)算出它們的聯(lián)合分布概率,本文主要考慮超過(guò)概率和條件概率,即當(dāng)總雨量X超過(guò)某特定值的條件下,峰值雨強(qiáng)超過(guò)某特定值時(shí)該雨型發(fā)生的條件概率,如式(17)所示。

    式中:XA表示概率為AA的總雨量值;y表示某特定峰值雨強(qiáng)值。

    由此,雨型風(fēng)險(xiǎn)組合概率(A,B)的定義為:

    在山洪災(zāi)害預(yù)警預(yù)報(bào)中,通常利用重現(xiàn)期對(duì)暴雨量級(jí)進(jìn)行描述,對(duì)于超過(guò)概率,重現(xiàn)期計(jì)算公式如下式所示。

    3 西峽縣基于臨界雨量的山洪災(zāi)害預(yù)警分析

    3.1 Copula函數(shù)的選擇

    分別采用Gumbel、Clayton、Frank 及混合Copula(M-Copula)函數(shù)構(gòu)建12 h 場(chǎng)次降雨量和峰值雨強(qiáng)聯(lián)合概率分布模型,K-S檢驗(yàn)統(tǒng)計(jì)量及擬合優(yōu)度評(píng)價(jià)指標(biāo)計(jì)算結(jié)果如表1所示。

    表1 計(jì) 算 表 明,Gumbel-Copula、Clayton-Copula、Frank-Copula、M-Copula 函數(shù)均通過(guò)K-S 檢驗(yàn),能較好的擬合臨界雨強(qiáng)的邊緣分布,根據(jù)均方根誤差RSME最小準(zhǔn)則,選取M-Copu‐la函數(shù)為連接函數(shù)來(lái)擬合降雨量與峰值雨強(qiáng)聯(lián)合分布情況。

    表1 Copula概率分布函數(shù)的計(jì)算、檢驗(yàn)與評(píng)價(jià)結(jié)果Tab.1 Calculation results,inspection and evaluation of Copula cumulative distribution function

    3.2 組合概率的臨界雨量的風(fēng)險(xiǎn)預(yù)警計(jì)算

    綜合考慮防災(zāi)對(duì)象所處河段河谷形態(tài)、洪水上漲速率、預(yù)警響應(yīng)時(shí)間和站點(diǎn)位置等因素,在臨界雨量的基礎(chǔ)上綜合確定準(zhǔn)備轉(zhuǎn)移和立即轉(zhuǎn)移的預(yù)警指標(biāo);并利用該預(yù)警指標(biāo)進(jìn)行暴雨洪水復(fù)核校正,以避免與成災(zāi)水位及相應(yīng)的暴雨洪水頻率差異過(guò)大。

    通常情況下,由于臨界雨量是從成災(zāi)水位對(duì)應(yīng)流量的洪水推算得到,故在數(shù)值上認(rèn)為臨界雨量即立即轉(zhuǎn)移指標(biāo);對(duì)于準(zhǔn)備轉(zhuǎn)移指標(biāo),是在臨界雨量基礎(chǔ)上根據(jù)準(zhǔn)備轉(zhuǎn)移時(shí)間及洪水過(guò)程線(xiàn)綜合進(jìn)行“折減”處理。

    以軍馬河鄉(xiāng)白果村回龍灣組為例,在不同的土壤含水量的情況下,利用風(fēng)險(xiǎn)概率雨型分別計(jì)算1、3、6、24 h 臨界雨量,其臨界雨量預(yù)警閾值見(jiàn)表2 所示;水位預(yù)警閾值見(jiàn)圖1 所示;準(zhǔn)備轉(zhuǎn)移預(yù)警和立即轉(zhuǎn)移預(yù)警雨量臨界線(xiàn)圖見(jiàn)圖2所示。

    圖1 水位預(yù)警閾值圖Fig.1 Threshold map of water level warning

    圖2 準(zhǔn)備轉(zhuǎn)移預(yù)警和立即轉(zhuǎn)移預(yù)警雨量臨界線(xiàn)圖Fig.2 Threshold maps for the preparation and immediate diversion of early warning rainfall

    表2 臨界雨量預(yù)警閾值表Tab.2 Table of critical rainfall warning thresholds

    3.3 結(jié)果合理性分析

    對(duì)于山丘區(qū)河流洪峰模數(shù)具有一定的規(guī)律,流域面積越大,洪峰模數(shù)越小,面積越小,洪峰模數(shù)越大。由此可以根據(jù)河南省山丘區(qū)水文斷面的洪峰模數(shù)統(tǒng)計(jì)資料,間接評(píng)估設(shè)計(jì)洪峰流量的合理性。根據(jù)以往經(jīng)驗(yàn)資料,50年一遇的中小河流斷面的設(shè)計(jì)洪峰流量對(duì)應(yīng)的洪峰模數(shù)一般在10.0 左右,并且10.0 以下占多數(shù)?;佚垶辰M計(jì)算得到洪峰模數(shù)為8.86,因此設(shè)計(jì)流量計(jì)算成果是合理的。

    表3 回龍灣組計(jì)算成果合理性分析Tab.3 Analysis on rationality of calculation results of Huilongwan

    圖3 各級(jí)水位對(duì)比圖Fig.3 Comparison chart of water level at different stages

    4 結(jié)論

    通過(guò)對(duì)西峽縣的山洪災(zāi)害預(yù)警分析可知:混合Copula 函數(shù)可以更加靈活的刻畫(huà)水文系統(tǒng)中各變量之間的相關(guān)關(guān)系,Copula 分布函數(shù)能夠很好的擬合總雨量和峰值雨強(qiáng)的二維聯(lián)合分布函數(shù)變量之間的相關(guān)結(jié)構(gòu)及分布,根據(jù)均方根誤差RSME最小準(zhǔn)則,選取M-Copula 函數(shù)為連接函數(shù)來(lái)擬合臨界雨量聯(lián)合分布情況;基于對(duì)西峽縣軍馬河鄉(xiāng)白果村回龍灣組為例進(jìn)行的計(jì)算及分析,計(jì)算組合概率的臨界雨量,并對(duì)山洪災(zāi)害的水位、準(zhǔn)備轉(zhuǎn)移預(yù)警及立即轉(zhuǎn)移預(yù)警雨量臨界進(jìn)行預(yù)警及分析。結(jié)果表明:對(duì)于山丘區(qū)河流洪峰模數(shù)具有一定的規(guī)律,流域面積越大,洪峰模數(shù)越小,面積越小,洪峰模數(shù)越大,設(shè)計(jì)流量為639 m3∕s 的計(jì)算成果是合理的。因此,對(duì)山洪災(zāi)害預(yù)警分析要最大限度的保證人民群眾的安全,分析評(píng)價(jià)中以各沿河村落漲洪歷時(shí)為基礎(chǔ),分析沿河村落暴雨洪水預(yù)警響應(yīng)時(shí)間及其分布情況,這對(duì)山洪災(zāi)害防治預(yù)案編制及預(yù)警等均具有重要作用。

    猜你喜歡
    山洪模數(shù)雨量
    寧夏紅柳溝流域水沙變化及產(chǎn)沙分析
    基于單片機(jī)和模數(shù)化設(shè)計(jì)的低壓側(cè)電壓監(jiān)視與保護(hù)裝置
    能源工程(2021年2期)2021-07-21 08:40:02
    模數(shù)化設(shè)計(jì)方法在景觀鋪裝設(shè)計(jì)中的應(yīng)用
    綠色科技(2020年11期)2020-08-01 02:23:58
    優(yōu)雅地表達(dá)
    基于小波去噪的稱(chēng)重雨量數(shù)據(jù)分析
    基于LID模式的城區(qū)排澇模數(shù)探析
    一種新型的RSA密碼體制模數(shù)分解算法
    SL—1 型雨量傳感器故障分析排除和維護(hù)
    西藏科技(2016年5期)2016-09-26 12:16:40
    遭遇暴雨山洪如何避險(xiǎn)自救
    湖北省山洪溝治理思路淺析
    国产无遮挡羞羞视频在线观看| 99国产精品免费福利视频| 国产又色又爽无遮挡免| 丝袜美腿诱惑在线| 免费看十八禁软件| 午夜福利免费观看在线| 久久中文字幕一级| 中文字幕人妻熟女乱码| 欧美日韩亚洲国产一区二区在线观看 | 高清视频免费观看一区二区| 日本a在线网址| 一区二区三区精品91| e午夜精品久久久久久久| 熟女少妇亚洲综合色aaa.| 国产成人av激情在线播放| 一区二区三区乱码不卡18| 亚洲一区中文字幕在线| 91成人精品电影| 久久天堂一区二区三区四区| 国产有黄有色有爽视频| 亚洲天堂av无毛| 中文精品一卡2卡3卡4更新| 欧美精品av麻豆av| 在线 av 中文字幕| av在线app专区| kizo精华| 天天操日日干夜夜撸| 亚洲人成电影免费在线| 亚洲国产精品一区三区| 天天躁日日躁夜夜躁夜夜| 婷婷色综合www| 热re99久久精品国产66热6| 亚洲欧美色中文字幕在线| 亚洲成av片中文字幕在线观看| 免费观看人在逋| 国产精品成人在线| 老鸭窝网址在线观看| 手机成人av网站| 久9热在线精品视频| 少妇精品久久久久久久| 五月天丁香电影| 大码成人一级视频| 免费高清在线观看视频在线观看| 久久久国产精品麻豆| 美女高潮到喷水免费观看| 欧美激情极品国产一区二区三区| 日韩av不卡免费在线播放| 少妇人妻 视频| 亚洲免费av在线视频| 色视频在线一区二区三区| 欧美 亚洲 国产 日韩一| 超碰97精品在线观看| 国产又爽黄色视频| 成在线人永久免费视频| 别揉我奶头~嗯~啊~动态视频 | 熟女少妇亚洲综合色aaa.| 色播在线永久视频| 国产爽快片一区二区三区| 久久毛片免费看一区二区三区| 国产精品av久久久久免费| 日本av免费视频播放| av网站在线播放免费| 国产免费一区二区三区四区乱码| 午夜视频精品福利| 国产精品国产av在线观看| 你懂的网址亚洲精品在线观看| 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 国产黄色视频一区二区在线观看| 美女午夜性视频免费| 国产精品.久久久| 成人手机av| 亚洲综合色网址| 男男h啪啪无遮挡| 亚洲三区欧美一区| 国产主播在线观看一区二区 | 久久国产精品人妻蜜桃| 国产老妇伦熟女老妇高清| 免费在线观看影片大全网站 | 免费高清在线观看视频在线观看| 久久久精品94久久精品| 水蜜桃什么品种好| 青春草视频在线免费观看| 90打野战视频偷拍视频| bbb黄色大片| 精品欧美一区二区三区在线| 免费在线观看视频国产中文字幕亚洲 | 成人影院久久| 亚洲专区国产一区二区| 日韩精品免费视频一区二区三区| 久久国产精品大桥未久av| 国产精品一国产av| 国产精品秋霞免费鲁丝片| xxxhd国产人妻xxx| 欧美亚洲日本最大视频资源| 久久亚洲国产成人精品v| 色精品久久人妻99蜜桃| www.精华液| 免费日韩欧美在线观看| 巨乳人妻的诱惑在线观看| 色94色欧美一区二区| 欧美精品一区二区免费开放| 一级毛片电影观看| 最近手机中文字幕大全| 日本91视频免费播放| 国产三级黄色录像| 欧美在线一区亚洲| 亚洲国产中文字幕在线视频| 亚洲一区中文字幕在线| 在线观看一区二区三区激情| 男女免费视频国产| e午夜精品久久久久久久| 精品一区在线观看国产| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲中文字幕日韩| 青春草视频在线免费观看| 成年美女黄网站色视频大全免费| 大片电影免费在线观看免费| 欧美日韩av久久| 午夜福利视频在线观看免费| cao死你这个sao货| 欧美中文综合在线视频| 日韩 亚洲 欧美在线| 美女中出高潮动态图| 中文字幕人妻丝袜制服| 成年女人毛片免费观看观看9 | 久久精品久久久久久久性| 精品国产超薄肉色丝袜足j| 久久久欧美国产精品| 久久99一区二区三区| 十分钟在线观看高清视频www| 亚洲精品久久午夜乱码| 国产爽快片一区二区三区| 九草在线视频观看| 国产熟女午夜一区二区三区| 国产免费视频播放在线视频| 久久精品国产a三级三级三级| 亚洲国产日韩一区二区| 国产欧美日韩综合在线一区二区| 久久 成人 亚洲| 亚洲精品中文字幕在线视频| 亚洲欧美精品自产自拍| 50天的宝宝边吃奶边哭怎么回事| 丁香六月天网| 这个男人来自地球电影免费观看| 一区二区三区乱码不卡18| 久久久久久久大尺度免费视频| 国产淫语在线视频| 人人妻人人澡人人爽人人夜夜| 香蕉国产在线看| 国产视频首页在线观看| 国产精品一二三区在线看| 国产又色又爽无遮挡免| 中文字幕人妻丝袜制服| 成人黄色视频免费在线看| 99热国产这里只有精品6| 老司机午夜十八禁免费视频| 日韩大片免费观看网站| 国产高清videossex| 免费在线观看日本一区| √禁漫天堂资源中文www| 黄色怎么调成土黄色| 女性被躁到高潮视频| 亚洲欧美色中文字幕在线| 51午夜福利影视在线观看| 亚洲国产av影院在线观看| 久久国产亚洲av麻豆专区| 男人舔女人的私密视频| 国产成人精品在线电影| 日韩中文字幕视频在线看片| 国产一区亚洲一区在线观看| 欧美精品一区二区免费开放| 啦啦啦在线免费观看视频4| 国产成人免费无遮挡视频| 999久久久国产精品视频| 国产在线免费精品| 精品一区二区三区四区五区乱码 | 18禁裸乳无遮挡动漫免费视频| 看免费av毛片| 999精品在线视频| 亚洲,欧美,日韩| 色精品久久人妻99蜜桃| 色视频在线一区二区三区| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| 在线观看免费午夜福利视频| 欧美 亚洲 国产 日韩一| 两个人免费观看高清视频| 免费看十八禁软件| 黄色a级毛片大全视频| 少妇粗大呻吟视频| 大香蕉久久网| 免费观看a级毛片全部| 国产成人免费无遮挡视频| www.自偷自拍.com| 性色av乱码一区二区三区2| 国产男人的电影天堂91| 无限看片的www在线观看| 欧美精品高潮呻吟av久久| 欧美97在线视频| 人妻一区二区av| 波野结衣二区三区在线| bbb黄色大片| 亚洲国产最新在线播放| 久久久国产欧美日韩av| 精品免费久久久久久久清纯 | 伦理电影免费视频| 黄色 视频免费看| 国产精品九九99| 国产一区二区在线观看av| 午夜福利影视在线免费观看| 国产精品久久久av美女十八| 大香蕉久久成人网| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 欧美激情高清一区二区三区| 一级毛片电影观看| 少妇粗大呻吟视频| 国产成人精品久久二区二区91| 亚洲精品国产色婷婷电影| 性色av一级| 国产欧美亚洲国产| 成年女人毛片免费观看观看9 | 一本综合久久免费| avwww免费| 日韩大码丰满熟妇| 好男人视频免费观看在线| 欧美97在线视频| 超色免费av| 一区二区三区四区激情视频| 黄片小视频在线播放| 欧美黄色淫秽网站| 一级黄片播放器| 91精品三级在线观看| 日韩伦理黄色片| 中文字幕精品免费在线观看视频| 久久精品熟女亚洲av麻豆精品| 在线观看国产h片| www日本在线高清视频| 欧美变态另类bdsm刘玥| 久9热在线精品视频| 国语对白做爰xxxⅹ性视频网站| 免费少妇av软件| 一级片免费观看大全| 久久久久久久久免费视频了| 老熟女久久久| 1024香蕉在线观看| 少妇裸体淫交视频免费看高清 | 精品高清国产在线一区| a级片在线免费高清观看视频| 十分钟在线观看高清视频www| 蜜桃国产av成人99| 国产黄色免费在线视频| 在现免费观看毛片| 亚洲av美国av| 成年女人毛片免费观看观看9 | 国产精品欧美亚洲77777| 欧美日韩黄片免| 国产精品一国产av| 一边摸一边抽搐一进一出视频| 午夜福利,免费看| 亚洲av男天堂| 日韩 亚洲 欧美在线| 欧美激情 高清一区二区三区| 日韩伦理黄色片| 国产亚洲av高清不卡| 十八禁网站网址无遮挡| 一级毛片 在线播放| 超碰成人久久| 欧美黑人欧美精品刺激| 我要看黄色一级片免费的| 新久久久久国产一级毛片| 日韩人妻精品一区2区三区| 日韩av免费高清视频| 999精品在线视频| 久久久亚洲精品成人影院| netflix在线观看网站| 亚洲精品美女久久久久99蜜臀 | 80岁老熟妇乱子伦牲交| 婷婷丁香在线五月| 99久久人妻综合| 国产一区二区三区av在线| 亚洲午夜精品一区,二区,三区| 天天添夜夜摸| av在线播放精品| 国产91精品成人一区二区三区 | 亚洲av男天堂| 精品人妻1区二区| 欧美日韩视频精品一区| 午夜两性在线视频| 日韩大码丰满熟妇| 一边摸一边做爽爽视频免费| 母亲3免费完整高清在线观看| 午夜福利影视在线免费观看| 人体艺术视频欧美日本| 亚洲成人手机| 亚洲成色77777| 青青草视频在线视频观看| 久久精品熟女亚洲av麻豆精品| 成人国产av品久久久| 大片免费播放器 马上看| 91精品三级在线观看| 亚洲欧洲精品一区二区精品久久久| 久久久国产欧美日韩av| 女人精品久久久久毛片| 精品人妻1区二区| 操出白浆在线播放| 亚洲国产成人一精品久久久| 男人操女人黄网站| 悠悠久久av| 最新的欧美精品一区二区| 一边摸一边抽搐一进一出视频| 国产熟女欧美一区二区| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 午夜免费鲁丝| 亚洲精品美女久久av网站| 国产精品一区二区免费欧美 | 亚洲精品自拍成人| 亚洲五月色婷婷综合| 最近最新中文字幕大全免费视频 | 18在线观看网站| 一边摸一边抽搐一进一出视频| 肉色欧美久久久久久久蜜桃| 美国免费a级毛片| 久久精品人人爽人人爽视色| 精品高清国产在线一区| 亚洲av电影在线进入| xxxhd国产人妻xxx| 欧美黑人欧美精品刺激| xxxhd国产人妻xxx| 大香蕉久久网| 99久久综合免费| 性高湖久久久久久久久免费观看| 国产免费视频播放在线视频| 午夜福利免费观看在线| 亚洲国产欧美网| 国产欧美亚洲国产| 美女福利国产在线| 亚洲av日韩精品久久久久久密 | 国产不卡av网站在线观看| 熟女少妇亚洲综合色aaa.| 在线观看一区二区三区激情| 亚洲精品乱久久久久久| 国产亚洲av片在线观看秒播厂| 国产成人一区二区三区免费视频网站 | 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 日本五十路高清| 日韩制服骚丝袜av| 精品人妻1区二区| 国产麻豆69| 亚洲精品国产av蜜桃| 欧美日韩一级在线毛片| 好男人电影高清在线观看| 免费在线观看黄色视频的| 丝袜在线中文字幕| 女警被强在线播放| 日本五十路高清| 国产成人91sexporn| 91字幕亚洲| 日韩伦理黄色片| 国产爽快片一区二区三区| 99热国产这里只有精品6| 午夜精品国产一区二区电影| 精品一区二区三区四区五区乱码 | 日本色播在线视频| 男人舔女人的私密视频| 老司机午夜十八禁免费视频| av在线老鸭窝| 亚洲国产中文字幕在线视频| 人妻 亚洲 视频| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| xxxhd国产人妻xxx| 久久精品人人爽人人爽视色| 国产免费一区二区三区四区乱码| 精品久久久久久久毛片微露脸 | 又紧又爽又黄一区二区| 丝袜人妻中文字幕| 丰满饥渴人妻一区二区三| 国产成人影院久久av| 亚洲成人免费电影在线观看 | 中文字幕亚洲精品专区| 国产一区二区在线观看av| 精品人妻1区二区| 一级毛片 在线播放| tube8黄色片| 国产一区亚洲一区在线观看| 男人添女人高潮全过程视频| √禁漫天堂资源中文www| 欧美变态另类bdsm刘玥| 精品一区二区三区av网在线观看 | 国产成人欧美| 99热国产这里只有精品6| 丝瓜视频免费看黄片| 亚洲国产av新网站| 久久久国产欧美日韩av| 日本猛色少妇xxxxx猛交久久| 考比视频在线观看| 午夜免费鲁丝| 午夜福利免费观看在线| 2021少妇久久久久久久久久久| 天堂中文最新版在线下载| 亚洲精品日本国产第一区| 久久久久久久久久久久大奶| 激情视频va一区二区三区| 麻豆国产av国片精品| 久久久久久久大尺度免费视频| 国产成人av教育| 乱人伦中国视频| 久久久久久亚洲精品国产蜜桃av| 欧美97在线视频| 久久久精品94久久精品| kizo精华| 免费一级毛片在线播放高清视频 | 国产一区二区三区综合在线观看| 国产免费视频播放在线视频| 成在线人永久免费视频| 亚洲午夜精品一区,二区,三区| 中国美女看黄片| 亚洲欧美激情在线| 国产熟女午夜一区二区三区| av天堂在线播放| 首页视频小说图片口味搜索 | 看免费av毛片| 精品国产一区二区三区久久久樱花| 日韩免费高清中文字幕av| 国产女主播在线喷水免费视频网站| 老司机在亚洲福利影院| 免费人妻精品一区二区三区视频| 天天躁夜夜躁狠狠躁躁| 国产一区亚洲一区在线观看| 波多野结衣av一区二区av| 可以免费在线观看a视频的电影网站| 日本欧美视频一区| 欧美另类一区| 啦啦啦 在线观看视频| 精品国产国语对白av| 国产成人精品在线电影| 国产亚洲午夜精品一区二区久久| 黄网站色视频无遮挡免费观看| 精品国产一区二区三区久久久樱花| 久久 成人 亚洲| kizo精华| 欧美大码av| 91麻豆精品激情在线观看国产 | 国产精品一二三区在线看| 别揉我奶头~嗯~啊~动态视频 | 免费女性裸体啪啪无遮挡网站| 下体分泌物呈黄色| 女人被躁到高潮嗷嗷叫费观| 最黄视频免费看| 亚洲伊人色综图| xxx大片免费视频| 人妻 亚洲 视频| 色视频在线一区二区三区| 人人妻人人澡人人爽人人夜夜| 免费在线观看视频国产中文字幕亚洲 | 久久天躁狠狠躁夜夜2o2o | 一本—道久久a久久精品蜜桃钙片| 日本欧美视频一区| 亚洲av电影在线观看一区二区三区| 国产欧美日韩精品亚洲av| 日日爽夜夜爽网站| 大陆偷拍与自拍| 男人爽女人下面视频在线观看| 在线观看免费日韩欧美大片| 自线自在国产av| 国产男人的电影天堂91| 亚洲,欧美精品.| 久久精品久久精品一区二区三区| 黑丝袜美女国产一区| 亚洲午夜精品一区,二区,三区| 大香蕉久久成人网| 男人舔女人的私密视频| xxx大片免费视频| 欧美激情 高清一区二区三区| 真人做人爱边吃奶动态| 狠狠婷婷综合久久久久久88av| 一级a爱视频在线免费观看| 免费观看a级毛片全部| 国产野战对白在线观看| 免费av中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 日韩一本色道免费dvd| 老司机靠b影院| 国产97色在线日韩免费| 18禁国产床啪视频网站| 99久久综合免费| 国产精品久久久久久人妻精品电影 | 欧美大码av| 久久久国产精品麻豆| 久久久久网色| 欧美人与性动交α欧美精品济南到| 91成人精品电影| 男女边吃奶边做爰视频| 美女高潮到喷水免费观看| 久久精品亚洲熟妇少妇任你| 欧美成狂野欧美在线观看| 韩国精品一区二区三区| 多毛熟女@视频| www.av在线官网国产| 色婷婷久久久亚洲欧美| 波多野结衣一区麻豆| 少妇精品久久久久久久| 青青草视频在线视频观看| 人体艺术视频欧美日本| 一边摸一边做爽爽视频免费| 日本vs欧美在线观看视频| 久久毛片免费看一区二区三区| 国产免费视频播放在线视频| av网站免费在线观看视频| 999精品在线视频| 日韩中文字幕视频在线看片| 高清视频免费观看一区二区| 我要看黄色一级片免费的| 免费人妻精品一区二区三区视频| 涩涩av久久男人的天堂| 丰满迷人的少妇在线观看| 99精国产麻豆久久婷婷| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 免费看不卡的av| 国产免费视频播放在线视频| 母亲3免费完整高清在线观看| 中文字幕制服av| 2018国产大陆天天弄谢| 国产一区二区三区综合在线观看| 五月天丁香电影| 观看av在线不卡| 美女午夜性视频免费| 国产99久久九九免费精品| 亚洲一码二码三码区别大吗| 久久久久久久大尺度免费视频| 午夜激情久久久久久久| 久久久久精品人妻al黑| 美女国产高潮福利片在线看| 最黄视频免费看| 午夜福利乱码中文字幕| 亚洲成人手机| 国产精品国产三级专区第一集| av在线app专区| 久久热在线av| 日本a在线网址| 亚洲av综合色区一区| 亚洲国产欧美在线一区| 久久精品国产亚洲av涩爱| 老司机亚洲免费影院| 午夜av观看不卡| 国产淫语在线视频| 国产成人a∨麻豆精品| av网站免费在线观看视频| 伊人久久大香线蕉亚洲五| 欧美少妇被猛烈插入视频| 国产男女超爽视频在线观看| 又紧又爽又黄一区二区| 一级毛片 在线播放| 男女床上黄色一级片免费看| 五月天丁香电影| 99热国产这里只有精品6| 午夜福利一区二区在线看| av国产久精品久网站免费入址| 久久99精品国语久久久| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 永久免费av网站大全| 在线亚洲精品国产二区图片欧美| 熟女少妇亚洲综合色aaa.| 久久精品人人爽人人爽视色| 国产日韩欧美视频二区| 女人精品久久久久毛片| 十八禁网站网址无遮挡| 777久久人妻少妇嫩草av网站| 日韩制服骚丝袜av| 亚洲男人天堂网一区| 老司机深夜福利视频在线观看 | 国精品久久久久久国模美| 天天操日日干夜夜撸| 久久精品国产a三级三级三级| 国产一区二区 视频在线| 亚洲七黄色美女视频| 久久久久久久大尺度免费视频| 黑人猛操日本美女一级片| 大陆偷拍与自拍| 婷婷色av中文字幕| 亚洲欧洲日产国产| 叶爱在线成人免费视频播放| 久久精品久久久久久久性| 大片电影免费在线观看免费| 国产精品一区二区在线观看99| 韩国高清视频一区二区三区| 一级片'在线观看视频| 国产一区二区激情短视频 | 亚洲,欧美,日韩| 欧美日韩av久久| 两人在一起打扑克的视频| 久久精品国产亚洲av涩爱| 777久久人妻少妇嫩草av网站| 免费不卡黄色视频| 亚洲精品美女久久久久99蜜臀 | 黄色片一级片一级黄色片| 国产欧美日韩一区二区三 | 激情视频va一区二区三区| 大片电影免费在线观看免费| 老司机影院毛片| 激情视频va一区二区三区| 国产精品免费大片| 丝袜美腿诱惑在线| 中文字幕制服av| 亚洲国产精品一区三区| 欧美在线一区亚洲| 精品人妻一区二区三区麻豆| 脱女人内裤的视频| 婷婷色麻豆天堂久久|