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

    GIS支持下的小流域山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃分方法研究

    2018-11-01 06:17:44駱建雄
    中國農(nóng)村水利水電 2018年10期
    關(guān)鍵詞:樂安易損易損性

    葛 星,駱建雄

    (河海大學(xué)地理信息科學(xué)與工程研究所,南京 211100)

    1 研究背景

    風(fēng)險(xiǎn)(Risk)概念最早在19世紀(jì)末的西方經(jīng)濟(jì)學(xué)領(lǐng)域提出,現(xiàn)在已經(jīng)廣泛應(yīng)用于經(jīng)濟(jì)學(xué)、社會(huì)學(xué)、自然災(zāi)害、環(huán)境科學(xué)等領(lǐng)域[1]。聯(lián)合國人道主義事務(wù)部(Department of Humainitarian UN)在1991年對(duì)風(fēng)險(xiǎn)作出定義:風(fēng)險(xiǎn)是指在一定區(qū)域和給定的時(shí)間內(nèi),由于特定的自然災(zāi)害引起的人們生命財(cái)產(chǎn)和經(jīng)濟(jì)活動(dòng)的期望損失值,可用公式R=H×V來表示,其中R代表風(fēng)險(xiǎn)度;H代表危險(xiǎn)度;V代表易損度[2]。

    國外有關(guān)洪水風(fēng)險(xiǎn)評(píng)估的研究主要成果有,Joy Sanyal等(2005)結(jié)合GIS與RS技術(shù)利用Landsat及ERS衛(wèi)星數(shù)據(jù)對(duì)山洪災(zāi)害的易損性進(jìn)行了研究[3, 4];Gunter BloschI等(2008)采用分布式水文模型研究了該模型在山洪災(zāi)害預(yù)測中的應(yīng)用[5];Mohammed E1 Bastawesy等(2009)利用GIS與RS技術(shù),通過計(jì)算地表徑流及降雨量模擬出該地區(qū)山洪的影響范圍及受災(zāi)程度[6]。我國學(xué)者在山洪災(zāi)害風(fēng)險(xiǎn)劃分方面最高開始于20世紀(jì)90年代,1996年趙士鵬對(duì)我國山洪災(zāi)害的危險(xiǎn)區(qū)劃做了定性的分析,將全國劃分為西北區(qū)、內(nèi)蒙區(qū)、青藏區(qū)、中部區(qū)、東部平原和東南區(qū)等6個(gè)區(qū)域[7];同年史培軍對(duì)災(zāi)害的致災(zāi)因子、孕災(zāi)環(huán)境和承災(zāi)體作出評(píng)述,形成了災(zāi)害研究的理論體系[8];2004年孫厚才等人對(duì)我國山洪災(zāi)害的特征、危害、成因、風(fēng)險(xiǎn)評(píng)估等進(jìn)行了概括、分析和總結(jié),為我國山洪災(zāi)害防治規(guī)劃的制定工作提供了參考[9];2007年張平倉等對(duì)山洪災(zāi)害的特征和防治對(duì)策做了全面的分析,為我國各地區(qū)山洪災(zāi)害防治對(duì)策提供了參考[10]。

    從我國山洪災(zāi)害風(fēng)險(xiǎn)評(píng)估的研究現(xiàn)狀來看,大部分學(xué)者都是對(duì)大尺度的區(qū)域進(jìn)行山洪災(zāi)害風(fēng)險(xiǎn)劃分,很少有人在縣/區(qū)級(jí)層面上進(jìn)行研究。本研究選擇江西省婺源縣樂安河上游流域?yàn)檠芯繀^(qū),對(duì)其進(jìn)行危險(xiǎn)性劃分和易損性劃分模型構(gòu)建,探索適用于小尺度的流域山洪災(zāi)害風(fēng)險(xiǎn)劃分方法。

    2 研究區(qū)及數(shù)據(jù)準(zhǔn)備

    2.1 研究區(qū)概況

    本研究選取江西省婺源縣樂安河上游流域作為研究區(qū),該區(qū)域?yàn)樯胶槎喟l(fā)區(qū)域,已經(jīng)作為江西省智慧流域建設(shè)的試點(diǎn)之一,數(shù)據(jù)較為齊全。

    樂安河(又稱玉溪河),為饒河支流,正源段莘水發(fā)源于婺源縣北部大庾山、五龍山南麓,南流經(jīng)段莘水庫至武口與古坦水匯合始名樂安河,整個(gè)樂安河流域全長279 km,流域面積8 989 km2。本研究選取的是樂安河上游流域,位于江西省婺源縣境內(nèi)(如圖1),流域面積約為1 416 km2,為亞熱帶季風(fēng)性濕潤氣候,雨量充沛、四季分明;歷年平均氣溫為16.8 ℃,多年平均降雨量為1 810 mm,流域境內(nèi)雨量分布不均,總體趨勢是東北多,西南少。

    圖1 樂安河上游地理位置圖Fig.1 Le'an River upstream location map

    2.2 數(shù)據(jù)準(zhǔn)備

    2.2.1 流域基礎(chǔ)地理數(shù)據(jù)

    為了更好了解樂安河流域下墊面情況,本研究收集了樂安河流域的遙感影像、土地利用數(shù)據(jù)、土壤類型和質(zhì)地?cái)?shù)據(jù)和數(shù)字高程模型。具體的數(shù)據(jù)情況見表1。

    表1 流域基礎(chǔ)數(shù)據(jù)情況表Tab.1 Basic data situation of the basin

    2.2.2 流域社會(huì)經(jīng)濟(jì)數(shù)據(jù)

    流域的社會(huì)經(jīng)濟(jì)數(shù)據(jù),主要包括樂安河上游流域的行政區(qū)劃、道路、人口和GDP等數(shù)據(jù),本研究收集到的數(shù)據(jù)如表2所示。

    表2 流域社會(huì)經(jīng)濟(jì)數(shù)據(jù)Tab.2 Social-economic data of the basin

    2.2.3 流域水文站網(wǎng)數(shù)據(jù)

    樂安河上游流域內(nèi)主要有水文站有2個(gè)(汪口、三都),水位站有6個(gè),雨量站有19個(gè),其地理位置分布見圖2。本研究收集了流域內(nèi)各類水文數(shù)據(jù)包括流域出口點(diǎn)、流域內(nèi)水系、流域最長匯流路徑、流域面、流域邊界、監(jiān)測站點(diǎn)(包括水文站、水位站和雨量站)和降雨數(shù)據(jù),如表3所示。

    表3 流域水文站網(wǎng)數(shù)據(jù)Tab.3 Hydrological station network data of the basin

    圖2 樂安河上游流域水系站網(wǎng)分布Fig.2 Distribution of Water System Network

    在所收集的水文數(shù)據(jù)中,樂安河上游流域被細(xì)化分為103個(gè)小流域(見圖3),每個(gè)小流域的都帶有已經(jīng)處理好的流域?qū)傩?見表4),包括流域的基本屬性如:流域面積、周長、形狀系數(shù)、中心點(diǎn)經(jīng)緯度;流域水文信息如:流域最長匯流路徑及比降、流域糙率和下滲率、洪峰模數(shù)等。

    表4 小流域?qū)傩詳?shù)據(jù)Tab.4 The attribution of the small watershed

    圖3 樂安河上游流域小流域劃分Fig.3 Classification of small watershed

    2.2.4 歷史山洪調(diào)查數(shù)據(jù)

    樂安河上游流域位于江西省4個(gè)暴雨中心之一,地形以中低山、丘陵為主,山區(qū)性河流坡降大,匯流速度快,非常容易受到山洪災(zāi)害的影響,根據(jù)江西省2015完成的山洪災(zāi)害調(diào)查評(píng)價(jià)項(xiàng)目中,得到婺源縣樂安河上游流域歷史山洪災(zāi)害點(diǎn)(見圖4)和山洪災(zāi)害調(diào)查評(píng)價(jià)項(xiàng)目中劃分的危險(xiǎn)區(qū)(見圖5)。

    表5 歷史山洪災(zāi)害發(fā)生情況Fig.5 The situation of the historical floods

    圖4 歷史山洪災(zāi)害點(diǎn)分布圖Fig.4 The distribution of the historical flood disaster

    圖5 山洪災(zāi)害防治單元(危險(xiǎn)區(qū))分布Fig.5 The distribution of the risk flood area

    3 研究方法

    3.1 危險(xiǎn)性劃分

    本文從風(fēng)險(xiǎn)的定義公式R=H×V出發(fā),其中R代表風(fēng)險(xiǎn)度,H代表危險(xiǎn)度,V代表易損度來展開研究。

    3.1.1 危險(xiǎn)性因子選擇

    根據(jù)查閱資料及相關(guān)規(guī)定,自然災(zāi)害的危險(xiǎn)性是指可能導(dǎo)致災(zāi)害的自然過程發(fā)生的規(guī)模頻率集總,是描述潛在的致災(zāi)過程易發(fā)程度的指標(biāo)。促使山洪災(zāi)害發(fā)生的原因主要有暴雨、流域特性、地形地貌、土地利用、土壤質(zhì)地特性等自然因素有關(guān),其中暴雨是誘發(fā)山洪災(zāi)害最重要的外力因素,地形地貌、土地利用、土壤等是山洪災(zāi)害形成的最基本的下墊面因子?;趯?duì)山洪災(zāi)害發(fā)生危險(xiǎn)性的理解和文獻(xiàn)的總結(jié),本研究選取以下8個(gè)指標(biāo)作為樂安河上游山洪災(zāi)害危險(xiǎn)性評(píng)估指標(biāo)(見圖6)。

    表6 危險(xiǎn)性因子指標(biāo)量化Tab.6 The quantification of the risk factor indicator

    圖6 危險(xiǎn)性因子指標(biāo)選取Fig.6 The choice of the risk factor indicator

    3.1.2 邏輯回歸模型構(gòu)建

    本文研究樂安河上游發(fā)生山洪災(zāi)害和未發(fā)生山洪災(zāi)害2個(gè)變量,因此采用二值化回歸模型(Binary Logistics)。在本研究的邏輯回歸分析中,因變量Y是一個(gè)二分類變量,其值取0或1,Y=1代表發(fā)生山洪災(zāi)害;Y=0代表未發(fā)生山洪災(zāi)害。影響Y取值的n個(gè)自變量分別為X1,X2,X3,…,Xn,在n個(gè)自變量的作用下山洪災(zāi)害發(fā)生的條件概率為P=P(Y=1|X1,X2,…,Xn),邏輯回歸模型可表示為:

    zi=a0+a1Xi1+a2Xi2+anXin

    (2)

    式中:zi代表中間變量參數(shù);a0代表回歸常數(shù);ai代表第j個(gè)變量的回歸系數(shù)(i,j=1,2,…,n);Xij代表第i個(gè)單元中第j個(gè)變量的取值,發(fā)生山洪災(zāi)害取1,否則取0;Pi代表第i號(hào)單元內(nèi)發(fā)生山洪災(zāi)害概率的回歸預(yù)測(i=1,2,…,n)。

    以已知的歷史山洪災(zāi)害發(fā)生點(diǎn)的位置及隨機(jī)生成的72個(gè)樣本點(diǎn)共144個(gè)樣本格網(wǎng)作為模型輸入數(shù)據(jù),整個(gè)模型的構(gòu)建和求解都是基于SPSS分析軟件中實(shí)現(xiàn),最終在8個(gè)危險(xiǎn)性指標(biāo)中選擇了3個(gè)(坡度、高程、和河網(wǎng)緩沖指標(biāo))具有顯著性指標(biāo)為最佳模型組合,具體結(jié)果見表7,可以得出3個(gè)指標(biāo)不同權(quán)重的系數(shù),其中每一類指標(biāo)分級(jí)的最后一級(jí)作為參考級(jí),系數(shù)均為0。

    表7 邏輯回歸危險(xiǎn)性指標(biāo)各級(jí)系數(shù)Tab.7 The coefficients of logistic regression model

    根據(jù)邏輯回歸結(jié)果的各級(jí)指標(biāo)系數(shù),求解可得整個(gè)樂安河上游流域山洪災(zāi)害危險(xiǎn)性劃分結(jié)果(如圖7所示),將山洪災(zāi)害發(fā)生的概率P(0~1)的值在0~0.2之間定為山洪災(zāi)害危險(xiǎn)性很低區(qū)域;在0.2~0.4之間定義為山洪災(zāi)害危險(xiǎn)性低區(qū)域;在0.4~0.6之間定義為山洪災(zāi)害危險(xiǎn)性中等區(qū)域;在0.6~0.8之間定義為山洪災(zāi)害危險(xiǎn)性高區(qū)域;在0.8~1之間定義為山洪災(zāi)害危險(xiǎn)性很高區(qū)域。

    圖7 山洪災(zāi)害危險(xiǎn)性邏輯回歸結(jié)果Fig.7 The results of the logistic regression model

    3.2 易損性劃分

    山洪災(zāi)害易損性是指發(fā)生山洪災(zāi)害后,暴露于危險(xiǎn)的人口財(cái)產(chǎn)以及人口財(cái)產(chǎn)對(duì)危險(xiǎn)的抵御能力。本研究根據(jù)文獻(xiàn)查閱及數(shù)據(jù)收集的情況,引用劉希林教授提出的災(zāi)害易損性指標(biāo),并對(duì)其模型進(jìn)行了簡化。由于研究區(qū)為婺源縣樂安河上游流域,在小尺度的研究中對(duì)于地區(qū)的固定資產(chǎn)投資可以忽略,所有本研究中去除了物質(zhì)易損性指標(biāo),保留了樂安河上游三類易損度指標(biāo),分別為人口易損度指標(biāo)(V1)、經(jīng)濟(jì)易損度指標(biāo)(V2)和環(huán)境易損度指標(biāo)(V3)。具體的模型公式如下:

    V=[(FV1+FV2+FV3)/3]0.5

    (3)

    (4)

    (5)

    (6)

    (7)

    式中:V代表山洪災(zāi)害易損程度;V1代表人口易損度指標(biāo);V2代表經(jīng)濟(jì)易損度指標(biāo);V3代表環(huán)境易損度指標(biāo);Bi代表各類土地資源基價(jià),元/m2;Ai代表各類土地利用類型面積,km2;FV1,F(xiàn)V2,F(xiàn)V3代表各類指標(biāo)定量賦值;m,n,l代表指定的常數(shù)項(xiàng)。

    在最終的山洪災(zāi)害易損度劃分結(jié)果中,易損度的取值是0~1之間的值,為了使最終的結(jié)果能與危險(xiǎn)性結(jié)果相匹配,本研究將山洪災(zāi)害易損性分為五個(gè)等級(jí):易損性值在0~0.2之間定為山洪災(zāi)害易損性很低區(qū)域;在0.2~0.4之間定義為山洪災(zāi)害易損性低區(qū)域;在0.4~0.6之間定義為山洪災(zāi)害易損性中等區(qū)域;在0.6~0.8之間定義為山洪災(zāi)害易損性高區(qū)域;在0.8~1之間定義為山洪災(zāi)害易損性很高區(qū)域(見圖9)。

    4 山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃分

    4.1 風(fēng)險(xiǎn)區(qū)劃分

    通過前面對(duì)危險(xiǎn)性和易損性模型的構(gòu)建,已經(jīng)分別得到樂安河上游流域山洪災(zāi)害危險(xiǎn)性和易損性劃分結(jié)果,根據(jù)聯(lián)合國人道主義提出的災(zāi)害風(fēng)險(xiǎn)的定義,本研究對(duì)樂安河流域上游山洪災(zāi)害風(fēng)險(xiǎn)定義為山洪災(zāi)害危險(xiǎn)性與易損性的乘積。由于之前對(duì)危險(xiǎn)性和易損性結(jié)果都分為5級(jí),分級(jí)范圍都是0~0.2;0.2~0.4;0.4~0.6;0.6~0.8;0.8~1,將危險(xiǎn)性與易損性相乘之后,可得到山洪災(zāi)害風(fēng)險(xiǎn)性的5個(gè)等級(jí)劃分(見圖10),分別為:0~0.04表示山洪災(zāi)害風(fēng)險(xiǎn)很低;0.04~0.16表示山洪災(zāi)害風(fēng)險(xiǎn)低;0.16~0.36表示山洪災(zāi)害風(fēng)險(xiǎn)中等;0.36~0.64表示山洪災(zāi)害風(fēng)險(xiǎn)高;0.64~1表示山洪災(zāi)害風(fēng)險(xiǎn)很高。

    圖8 易損性因子指標(biāo)選取Fig.8 The choice of Vulnerability factor indicator

    圖9 山洪災(zāi)害易損性轉(zhuǎn)換賦值函數(shù)結(jié)果Fig.9 The result of the transformation assignment function

    圖10 山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果Fig.10 Flash flood risk classification results

    4.2 精度驗(yàn)證

    在2015年完成的江西省山洪災(zāi)害調(diào)查評(píng)價(jià)項(xiàng)目中,可以提取婺源縣樂安河上游流域范圍內(nèi)居民點(diǎn)的防治類型,其中居民點(diǎn)防治類型分為低、中、高三個(gè)等級(jí)。本研究采用防治類型中和高兩個(gè)等級(jí)的防治點(diǎn)(共134個(gè))對(duì)上述風(fēng)險(xiǎn)劃分的結(jié)果進(jìn)行檢驗(yàn)(見圖11),從表8的統(tǒng)計(jì)結(jié)果來看,在流域面積33.18%的中等、高和很高的山洪災(zāi)害風(fēng)險(xiǎn)區(qū)域中,有79.11%的已劃分的山洪災(zāi)害防治類型點(diǎn)(中和高)被驗(yàn)證,說明樂安河山洪災(zāi)害劃分的結(jié)果具有合理性,本研究的方法適合較小尺度的山丘區(qū)流域山洪災(zāi)害風(fēng)險(xiǎn)劃分。

    表8 山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果驗(yàn)證表Tab.8 The verification table of classification result

    圖11 山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果驗(yàn)證Fig.11 The verification of classification result

    5 結(jié) 論

    本研究根據(jù)聯(lián)合國人道主義事務(wù)部對(duì)風(fēng)險(xiǎn)的定義,將樂安河上游流域山洪災(zāi)害風(fēng)險(xiǎn)定為流域內(nèi)山洪災(zāi)害的危險(xiǎn)性與易損性的乘積,主要研究了邏輯回歸模型對(duì)山洪災(zāi)害危險(xiǎn)性的劃分,轉(zhuǎn)換賦值函數(shù)對(duì)山洪災(zāi)害易損性的劃分,并對(duì)最終的山洪災(zāi)害風(fēng)險(xiǎn)劃分進(jìn)行了合理性驗(yàn)證。具體研究結(jié)果如下:

    (1)對(duì)山洪災(zāi)害的危險(xiǎn)性采用二分類的邏輯回歸模型進(jìn)行模擬求解,以歷史山洪災(zāi)害發(fā)生的實(shí)際位置點(diǎn)結(jié)合流域內(nèi)隨機(jī)生成的等量隨機(jī)點(diǎn)作為模型樣本,將山洪災(zāi)害發(fā)的概率作為因變量,選取流域8個(gè)危險(xiǎn)性指標(biāo)作為自變量,經(jīng)模型回歸后最終選取流域坡度、高程和河網(wǎng)緩沖作為最優(yōu)模型的自變量,得出流域內(nèi)山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果。

    (2)對(duì)山洪災(zāi)害易損性采用轉(zhuǎn)換賦值函數(shù)法,將流域山洪災(zāi)害易損性定義為人口易損度指標(biāo)、經(jīng)濟(jì)易損度指標(biāo)和環(huán)境易損度指標(biāo)的綜合評(píng)價(jià)結(jié)果,得出流域山洪災(zāi)害易損性結(jié)果。

    (3)通過山洪災(zāi)害危險(xiǎn)性和易損性的乘積得到最終樂安河上游流域山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果,并利用山洪災(zāi)害調(diào)查評(píng)價(jià)項(xiàng)目中的防治點(diǎn)對(duì)最終結(jié)果進(jìn)行驗(yàn)證,結(jié)果表明樂安河上游小流域山洪災(zāi)害風(fēng)險(xiǎn)劃分結(jié)果較為合理,證明了該方法對(duì)小流域山洪災(zāi)害風(fēng)險(xiǎn)區(qū)劃分是可行的。

    猜你喜歡
    樂安易損易損性
    基于IDA的預(yù)應(yīng)力混凝土連續(xù)梁橋易損性分析
    江西省樂安老年大學(xué)校歌
    姜黃素誘協(xié)同阿托伐他汀穩(wěn)定易損斑塊的作用與機(jī)制
    詩二首
    中華魂(2017年3期)2017-11-22 02:01:55
    民族英雄趙一曼之歌
    中華魂(2017年5期)2017-11-21 01:53:23
    來自銀時(shí)空的紙條
    更正聲明
    基于PSDM和IDA法的深水隔震橋梁地震易損性分析比較
    超聲造影聯(lián)合常規(guī)超聲觀察阿托伐他汀治療頸動(dòng)脈易損斑塊的療效
    更正聲明
    精品久久久久久久久av| 人人妻人人澡人人爽人人夜夜| 99re6热这里在线精品视频| a级毛色黄片| 最近中文字幕高清免费大全6| 99久久精品热视频| 成人黄色视频免费在线看| 国产成人一区二区在线| 97精品久久久久久久久久精品| 欧美一区二区亚洲| 亚洲av二区三区四区| 国产一区有黄有色的免费视频| 国产精品一二三区在线看| 18禁在线无遮挡免费观看视频| 少妇人妻 视频| 国产探花极品一区二区| 国产一级毛片在线| 午夜日本视频在线| 久久久久久久久久久丰满| 水蜜桃什么品种好| 男女下面进入的视频免费午夜| 国产熟女欧美一区二区| 成人综合一区亚洲| 国产免费一区二区三区四区乱码| 国精品久久久久久国模美| 日韩三级伦理在线观看| xxx大片免费视频| 国产精品成人在线| 国产高清不卡午夜福利| 国产黄片美女视频| 夜夜看夜夜爽夜夜摸| 在线播放无遮挡| 伦精品一区二区三区| 99久久精品一区二区三区| 人妻系列 视频| 国模一区二区三区四区视频| 国产成人a∨麻豆精品| av国产精品久久久久影院| 国产精品免费大片| 国产黄色免费在线视频| 国产成人一区二区在线| 日韩电影二区| 熟妇人妻不卡中文字幕| 极品教师在线视频| 一区二区三区免费毛片| 夫妻午夜视频| 一边亲一边摸免费视频| 亚洲精品国产色婷婷电影| 噜噜噜噜噜久久久久久91| 久久精品国产自在天天线| 久久99精品国语久久久| 国产成人91sexporn| 免费av中文字幕在线| 一边亲一边摸免费视频| 黑人高潮一二区| 亚洲精品色激情综合| 99热网站在线观看| 爱豆传媒免费全集在线观看| 亚洲av男天堂| 高清在线视频一区二区三区| 久久ye,这里只有精品| 国产真实伦视频高清在线观看| 免费观看在线日韩| 国产精品国产av在线观看| 99热全是精品| 99九九线精品视频在线观看视频| 蜜臀久久99精品久久宅男| 色5月婷婷丁香| 欧美丝袜亚洲另类| 欧美日韩国产mv在线观看视频 | 男女无遮挡免费网站观看| 91在线精品国自产拍蜜月| 国产成人精品久久久久久| 国产精品不卡视频一区二区| 国产精品熟女久久久久浪| 在线观看美女被高潮喷水网站| 亚洲天堂av无毛| 精品熟女少妇av免费看| 人妻系列 视频| 免费久久久久久久精品成人欧美视频 | 日韩在线高清观看一区二区三区| 三级国产精品片| 麻豆国产97在线/欧美| 日韩一本色道免费dvd| 如何舔出高潮| 一个人看视频在线观看www免费| 丰满迷人的少妇在线观看| 99热网站在线观看| 在线观看av片永久免费下载| 午夜视频国产福利| 久久久久久久大尺度免费视频| 中文乱码字字幕精品一区二区三区| 精品少妇久久久久久888优播| 日韩精品有码人妻一区| 亚洲真实伦在线观看| 性高湖久久久久久久久免费观看| 日韩在线高清观看一区二区三区| 99久久精品一区二区三区| 99热这里只有精品一区| 免费看日本二区| 又爽又黄a免费视频| 内地一区二区视频在线| 成年免费大片在线观看| 成年人午夜在线观看视频| 久久ye,这里只有精品| 久久久久久伊人网av| 亚洲美女黄色视频免费看| 男女边摸边吃奶| 天天躁夜夜躁狠狠久久av| 高清黄色对白视频在线免费看 | 18禁动态无遮挡网站| 我的老师免费观看完整版| 日韩国内少妇激情av| 国产乱人偷精品视频| 国产免费视频播放在线视频| 国产成人免费无遮挡视频| 这个男人来自地球电影免费观看 | 丰满少妇做爰视频| 亚州av有码| 一级毛片aaaaaa免费看小| 香蕉精品网在线| 久久久久久久国产电影| 国产高清国产精品国产三级 | 男男h啪啪无遮挡| 少妇丰满av| 亚洲电影在线观看av| 多毛熟女@视频| 免费人妻精品一区二区三区视频| 女性生殖器流出的白浆| 99热这里只有是精品在线观看| 天美传媒精品一区二区| 多毛熟女@视频| 欧美+日韩+精品| 美女视频免费永久观看网站| 最近手机中文字幕大全| 97热精品久久久久久| 久久久久久久精品精品| 日韩中文字幕视频在线看片 | 亚洲高清免费不卡视频| 十分钟在线观看高清视频www | 精品一区二区三区视频在线| 久久毛片免费看一区二区三区| 成年av动漫网址| 熟妇人妻不卡中文字幕| 超碰97精品在线观看| 亚洲综合精品二区| 亚洲欧美日韩卡通动漫| 日日撸夜夜添| 高清午夜精品一区二区三区| 欧美高清成人免费视频www| a级毛色黄片| 一区二区三区乱码不卡18| 丰满迷人的少妇在线观看| 26uuu在线亚洲综合色| 久久精品国产自在天天线| 日本午夜av视频| 天天躁夜夜躁狠狠久久av| 欧美xxⅹ黑人| 黄色配什么色好看| 欧美激情国产日韩精品一区| 国产精品麻豆人妻色哟哟久久| 特大巨黑吊av在线直播| 少妇裸体淫交视频免费看高清| 久久久欧美国产精品| 国产成人a区在线观看| 丝袜喷水一区| 国产 一区 欧美 日韩| h视频一区二区三区| 日韩电影二区| 内地一区二区视频在线| 18禁在线无遮挡免费观看视频| 黄片wwwwww| 街头女战士在线观看网站| 亚洲久久久国产精品| 成年女人在线观看亚洲视频| 秋霞伦理黄片| 亚洲精品国产色婷婷电影| 人人妻人人添人人爽欧美一区卜 | 国产午夜精品久久久久久一区二区三区| 高清在线视频一区二区三区| 亚洲精品一区蜜桃| 精品人妻熟女av久视频| 夜夜骑夜夜射夜夜干| 久久国产亚洲av麻豆专区| 亚州av有码| 22中文网久久字幕| 少妇裸体淫交视频免费看高清| 一边亲一边摸免费视频| 最近中文字幕2019免费版| 久久精品国产自在天天线| 大香蕉久久网| 国产毛片在线视频| 大陆偷拍与自拍| 午夜日本视频在线| 九色成人免费人妻av| 卡戴珊不雅视频在线播放| 最后的刺客免费高清国语| av在线app专区| 一级毛片黄色毛片免费观看视频| 视频区图区小说| 美女福利国产在线 | 一级毛片aaaaaa免费看小| 黄色一级大片看看| 欧美精品一区二区大全| 日韩伦理黄色片| 成人特级av手机在线观看| 欧美成人一区二区免费高清观看| 国精品久久久久久国模美| 一本色道久久久久久精品综合| 日韩大片免费观看网站| 大片电影免费在线观看免费| 99久久精品热视频| 观看av在线不卡| 久久热精品热| 在线天堂最新版资源| 九草在线视频观看| 日本-黄色视频高清免费观看| 亚洲第一区二区三区不卡| 亚洲av中文av极速乱| 久久精品国产亚洲av涩爱| 久久鲁丝午夜福利片| 女性生殖器流出的白浆| 熟女av电影| 国产成人午夜福利电影在线观看| 少妇高潮的动态图| 永久免费av网站大全| 国产爱豆传媒在线观看| 日本午夜av视频| av又黄又爽大尺度在线免费看| 欧美国产精品一级二级三级 | 亚洲成人中文字幕在线播放| 久久精品国产自在天天线| 亚洲av成人精品一二三区| 精品一区二区免费观看| 少妇的逼好多水| 人人妻人人添人人爽欧美一区卜 | 久久久色成人| 免费观看的影片在线观看| 国产精品一区www在线观看| 久久精品国产亚洲网站| 一个人免费看片子| 尤物成人国产欧美一区二区三区| 欧美97在线视频| av专区在线播放| 黄色一级大片看看| 永久网站在线| 免费观看av网站的网址| 交换朋友夫妻互换小说| 亚洲性久久影院| 18禁裸乳无遮挡动漫免费视频| 最近最新中文字幕免费大全7| 国产成人精品福利久久| 亚洲av.av天堂| 一本—道久久a久久精品蜜桃钙片| 91久久精品国产一区二区三区| 国产精品一及| 久久久久久伊人网av| 亚洲国产欧美人成| 久久午夜福利片| 大香蕉97超碰在线| 国模一区二区三区四区视频| 91久久精品国产一区二区成人| 国产精品一区二区在线不卡| 美女xxoo啪啪120秒动态图| 最新中文字幕久久久久| 精品午夜福利在线看| 国产精品爽爽va在线观看网站| 久久青草综合色| a 毛片基地| 免费看av在线观看网站| 内地一区二区视频在线| 另类亚洲欧美激情| 内射极品少妇av片p| av网站免费在线观看视频| 国产精品av视频在线免费观看| 最近最新中文字幕免费大全7| 日本爱情动作片www.在线观看| 亚洲国产欧美在线一区| 国产成人精品婷婷| 色视频在线一区二区三区| 国产精品人妻久久久影院| 日韩强制内射视频| 最近中文字幕高清免费大全6| 久久6这里有精品| 美女内射精品一级片tv| 女性被躁到高潮视频| 精品99又大又爽又粗少妇毛片| 国产精品蜜桃在线观看| 亚洲四区av| 久久午夜福利片| 久久久色成人| 欧美成人精品欧美一级黄| 联通29元200g的流量卡| 最近最新中文字幕大全电影3| 成年av动漫网址| 亚洲精品国产av蜜桃| 国产深夜福利视频在线观看| 亚洲第一av免费看| 精品视频人人做人人爽| 欧美人与善性xxx| 免费av不卡在线播放| 亚洲精品久久久久久婷婷小说| 国产一区二区在线观看日韩| 我的女老师完整版在线观看| 青青草视频在线视频观看| 亚洲国产色片| 乱码一卡2卡4卡精品| 伦理电影大哥的女人| 美女脱内裤让男人舔精品视频| 欧美xxxx性猛交bbbb| av国产久精品久网站免费入址| 亚洲成人一二三区av| 亚洲精华国产精华液的使用体验| 五月开心婷婷网| 精品国产乱码久久久久久小说| 精品国产三级普通话版| 综合色丁香网| 亚洲国产色片| 久久国产精品男人的天堂亚洲 | 午夜福利影视在线免费观看| 我要看黄色一级片免费的| 我要看黄色一级片免费的| 国产无遮挡羞羞视频在线观看| 国产高清三级在线| 肉色欧美久久久久久久蜜桃| 久久精品久久久久久噜噜老黄| 一二三四中文在线观看免费高清| 观看美女的网站| 国产精品欧美亚洲77777| 91精品伊人久久大香线蕉| 韩国av在线不卡| 国产精品99久久99久久久不卡 | a 毛片基地| 成人美女网站在线观看视频| 一个人看视频在线观看www免费| 亚洲av中文av极速乱| 日韩av免费高清视频| 日本一二三区视频观看| 尤物成人国产欧美一区二区三区| 大香蕉久久网| 久久精品久久久久久噜噜老黄| 欧美日韩视频高清一区二区三区二| 麻豆成人午夜福利视频| 一级爰片在线观看| 肉色欧美久久久久久久蜜桃| 亚洲天堂av无毛| 人人妻人人添人人爽欧美一区卜 | 老司机影院成人| 日本欧美国产在线视频| 亚洲国产成人一精品久久久| 亚洲欧美成人精品一区二区| 高清黄色对白视频在线免费看 | www.av在线官网国产| 人妻夜夜爽99麻豆av| 国产男女内射视频| 在线观看免费高清a一片| 99久久精品国产国产毛片| 中文字幕免费在线视频6| 在线播放无遮挡| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 一级黄片播放器| 久久精品熟女亚洲av麻豆精品| 夫妻性生交免费视频一级片| 在线观看免费日韩欧美大片 | 人妻少妇偷人精品九色| 女性被躁到高潮视频| 亚洲精品一区蜜桃| 国产黄片美女视频| 精品人妻偷拍中文字幕| 九九久久精品国产亚洲av麻豆| 国产亚洲最大av| 欧美精品亚洲一区二区| 亚洲人成网站高清观看| 亚洲,欧美,日韩| 久久99热这里只有精品18| 亚洲精品国产成人久久av| 久久精品久久精品一区二区三区| 免费看av在线观看网站| 少妇的逼好多水| 亚洲成人中文字幕在线播放| 国产精品秋霞免费鲁丝片| 国产欧美另类精品又又久久亚洲欧美| 国产免费一区二区三区四区乱码| 精品亚洲成a人片在线观看 | 高清日韩中文字幕在线| 青青草视频在线视频观看| 99热这里只有是精品在线观看| 亚洲,欧美,日韩| 成人国产麻豆网| 毛片一级片免费看久久久久| 秋霞在线观看毛片| 亚洲av中文字字幕乱码综合| 国产爱豆传媒在线观看| 七月丁香在线播放| 欧美zozozo另类| 少妇熟女欧美另类| 亚洲性久久影院| 少妇被粗大猛烈的视频| 亚洲婷婷狠狠爱综合网| 极品教师在线视频| 欧美人与善性xxx| 日日啪夜夜撸| 色哟哟·www| 狂野欧美白嫩少妇大欣赏| 美女中出高潮动态图| 亚洲精华国产精华液的使用体验| 精品人妻一区二区三区麻豆| www.av在线官网国产| 日韩不卡一区二区三区视频在线| 一级毛片我不卡| 日韩电影二区| av播播在线观看一区| 亚洲欧美日韩卡通动漫| 国产欧美日韩一区二区三区在线 | 精品久久国产蜜桃| 在线亚洲精品国产二区图片欧美 | 国产精品人妻久久久影院| 亚洲国产最新在线播放| 新久久久久国产一级毛片| 天堂中文最新版在线下载| 大码成人一级视频| 国产中年淑女户外野战色| 色5月婷婷丁香| 亚洲欧美日韩无卡精品| 成人漫画全彩无遮挡| 成年av动漫网址| 欧美xxⅹ黑人| 在线观看国产h片| 久久精品久久久久久久性| 黄色日韩在线| 777米奇影视久久| 干丝袜人妻中文字幕| 边亲边吃奶的免费视频| 97热精品久久久久久| 亚洲av成人精品一区久久| 大香蕉久久网| 精品人妻视频免费看| 91精品伊人久久大香线蕉| 熟妇人妻不卡中文字幕| 亚洲,欧美,日韩| 我的女老师完整版在线观看| 大香蕉97超碰在线| 久久久欧美国产精品| av播播在线观看一区| 久久久久人妻精品一区果冻| 99热这里只有是精品在线观看| 国产探花极品一区二区| 精品国产乱码久久久久久小说| 婷婷色综合www| 小蜜桃在线观看免费完整版高清| 久久99精品国语久久久| 天天躁日日操中文字幕| 日韩电影二区| 深夜a级毛片| 舔av片在线| 国精品久久久久久国模美| 天堂俺去俺来也www色官网| 亚洲欧美日韩东京热| 亚洲综合色惰| 乱系列少妇在线播放| 内射极品少妇av片p| 成人国产麻豆网| 女人十人毛片免费观看3o分钟| 男人舔奶头视频| 大话2 男鬼变身卡| 男女国产视频网站| 极品少妇高潮喷水抽搐| 亚洲精品日韩在线中文字幕| av国产精品久久久久影院| 国产 一区 欧美 日韩| 亚洲精品久久久久久婷婷小说| 国产免费一区二区三区四区乱码| 18+在线观看网站| 视频区图区小说| 五月天丁香电影| 久久久久久久久久成人| 看非洲黑人一级黄片| 精品久久久精品久久久| a级一级毛片免费在线观看| 亚洲成人中文字幕在线播放| 亚洲四区av| 内射极品少妇av片p| 观看av在线不卡| 黄色视频在线播放观看不卡| 免费大片黄手机在线观看| 免费久久久久久久精品成人欧美视频 | 国产精品一二三区在线看| 色网站视频免费| 国产精品一区二区在线不卡| 精品亚洲成a人片在线观看 | 亚洲av免费高清在线观看| 成人高潮视频无遮挡免费网站| 精品久久久久久久久av| 亚洲美女搞黄在线观看| 日韩人妻高清精品专区| 久久午夜福利片| 亚洲精华国产精华液的使用体验| 亚洲av男天堂| 国产无遮挡羞羞视频在线观看| 日本-黄色视频高清免费观看| 人人妻人人添人人爽欧美一区卜 | 亚洲最大成人中文| 国产成人精品婷婷| 国产爽快片一区二区三区| 22中文网久久字幕| 美女内射精品一级片tv| 亚洲精华国产精华液的使用体验| 国产成人免费观看mmmm| 热re99久久精品国产66热6| 三级国产精品欧美在线观看| 亚洲色图综合在线观看| 午夜免费观看性视频| 99久久中文字幕三级久久日本| 日本黄色片子视频| 亚洲av欧美aⅴ国产| 18禁裸乳无遮挡免费网站照片| 欧美激情国产日韩精品一区| 国产亚洲91精品色在线| 狠狠精品人妻久久久久久综合| 婷婷色麻豆天堂久久| 国产在视频线精品| 亚洲最大成人中文| 男女国产视频网站| 99热这里只有是精品50| 男女下面进入的视频免费午夜| 在线免费观看不下载黄p国产| 一级毛片黄色毛片免费观看视频| 人人妻人人添人人爽欧美一区卜 | 精品酒店卫生间| 中国国产av一级| 青春草视频在线免费观看| 又大又黄又爽视频免费| 美女高潮的动态| 久久久久久九九精品二区国产| 18禁裸乳无遮挡免费网站照片| 亚洲最大成人中文| 激情 狠狠 欧美| 欧美97在线视频| 亚洲综合色惰| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 亚洲国产色片| 美女脱内裤让男人舔精品视频| 日韩成人av中文字幕在线观看| 国产av一区二区精品久久 | videossex国产| 99视频精品全部免费 在线| 亚洲国产最新在线播放| 亚洲经典国产精华液单| 久久青草综合色| 国产精品国产av在线观看| 大片电影免费在线观看免费| 中文字幕人妻熟人妻熟丝袜美| 日韩精品有码人妻一区| 天堂俺去俺来也www色官网| av在线蜜桃| 午夜福利视频精品| 亚洲精华国产精华液的使用体验| 亚洲内射少妇av| 日韩电影二区| 欧美精品国产亚洲| 亚洲成人一二三区av| 中文字幕av成人在线电影| 99久国产av精品国产电影| 色网站视频免费| 卡戴珊不雅视频在线播放| 亚洲国产精品999| 国产精品一区二区在线不卡| h视频一区二区三区| 七月丁香在线播放| 直男gayav资源| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 99视频精品全部免费 在线| 久久午夜福利片| 日韩国内少妇激情av| 男女国产视频网站| 久热这里只有精品99| 又粗又硬又长又爽又黄的视频| av黄色大香蕉| 精品久久久久久电影网| 如何舔出高潮| 亚洲欧美精品专区久久| 三级国产精品片| 国产精品一区二区性色av| 日韩免费高清中文字幕av| 成人亚洲精品一区在线观看 | 日韩欧美精品免费久久| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 我要看日韩黄色一级片| 成人午夜精彩视频在线观看| 七月丁香在线播放| 涩涩av久久男人的天堂| 搡女人真爽免费视频火全软件| 国产精品伦人一区二区| 男人舔奶头视频| 一本一本综合久久| 精品人妻熟女av久视频| 永久网站在线| 我要看黄色一级片免费的| 久久久久久久久久人人人人人人| 在线观看免费高清a一片| 欧美激情极品国产一区二区三区 | 亚洲av男天堂| av线在线观看网站| 亚洲伊人久久精品综合| 又黄又爽又刺激的免费视频.| 啦啦啦在线观看免费高清www| 久久婷婷青草| 国产黄片视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美中文字幕日韩二区|