• <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)脈易損斑塊的療效
    更正聲明
    十八禁人妻一区二区| 亚洲成人久久性| 午夜精品久久久久久毛片777| 人人妻,人人澡人人爽秒播| 韩国av一区二区三区四区| 在线播放国产精品三级| 午夜福利在线观看吧| 久久国产乱子伦精品免费另类| 成人欧美大片| 三级毛片av免费| 少妇裸体淫交视频免费看高清 | 无遮挡黄片免费观看| 欧美精品亚洲一区二区| 亚洲情色 制服丝袜| 亚洲av熟女| 宅男免费午夜| 少妇粗大呻吟视频| 国产成人精品久久二区二区免费| 老司机午夜十八禁免费视频| 一个人免费在线观看的高清视频| 一a级毛片在线观看| 999久久久精品免费观看国产| 久久久久国产精品人妻aⅴ院| 亚洲精品国产色婷婷电影| 日本免费a在线| 亚洲av美国av| 国产精品久久久久久亚洲av鲁大| 免费无遮挡裸体视频| 精品欧美一区二区三区在线| 69av精品久久久久久| 天堂√8在线中文| 美女午夜性视频免费| 12—13女人毛片做爰片一| 搞女人的毛片| 国产在线精品亚洲第一网站| 午夜老司机福利片| 久久久精品欧美日韩精品| 国产三级在线视频| 久久久久久免费高清国产稀缺| 国产亚洲精品av在线| √禁漫天堂资源中文www| 日韩av在线大香蕉| 国产一区在线观看成人免费| 麻豆久久精品国产亚洲av| 老司机深夜福利视频在线观看| 一级作爱视频免费观看| 久久人妻熟女aⅴ| 国产91精品成人一区二区三区| 纯流量卡能插随身wifi吗| 看黄色毛片网站| 成人18禁在线播放| 久久久久久久久久久久大奶| 女生性感内裤真人,穿戴方法视频| 亚洲第一av免费看| 亚洲精品一卡2卡三卡4卡5卡| 成人18禁高潮啪啪吃奶动态图| 亚洲 国产 在线| 欧美日本视频| 久久热在线av| 亚洲国产精品久久男人天堂| 国产麻豆69| 免费人成视频x8x8入口观看| 欧美激情极品国产一区二区三区| 国产精品香港三级国产av潘金莲| 女人被躁到高潮嗷嗷叫费观| 黄网站色视频无遮挡免费观看| 女性被躁到高潮视频| 波多野结衣巨乳人妻| 国产激情欧美一区二区| 制服诱惑二区| 日韩欧美免费精品| 久久久国产精品麻豆| 美女午夜性视频免费| 999精品在线视频| 欧美性长视频在线观看| 少妇的丰满在线观看| 黄片大片在线免费观看| 一二三四社区在线视频社区8| 精品人妻在线不人妻| 伊人久久大香线蕉亚洲五| 999久久久精品免费观看国产| 老司机深夜福利视频在线观看| 真人做人爱边吃奶动态| 可以免费在线观看a视频的电影网站| 欧美成人一区二区免费高清观看 | 国产99久久九九免费精品| 99久久综合精品五月天人人| 国产欧美日韩综合在线一区二区| 大香蕉久久成人网| av在线天堂中文字幕| 国产又色又爽无遮挡免费看| 午夜福利欧美成人| 亚洲成人久久性| 999久久久精品免费观看国产| 可以在线观看毛片的网站| 可以免费在线观看a视频的电影网站| 欧美日韩福利视频一区二区| 精品国产乱码久久久久久男人| 他把我摸到了高潮在线观看| 日本免费一区二区三区高清不卡 | 在线观看一区二区三区| 在线观看免费视频网站a站| 嫩草影视91久久| 欧美日韩黄片免| 波多野结衣高清无吗| 中文亚洲av片在线观看爽| 国产精品国产高清国产av| 女人爽到高潮嗷嗷叫在线视频| 亚洲avbb在线观看| 国产97色在线日韩免费| 国产高清有码在线观看视频 | 久久人妻熟女aⅴ| 18禁美女被吸乳视频| 久久久精品欧美日韩精品| 久久久久久大精品| 午夜福利,免费看| 亚洲精品一区av在线观看| 欧美日韩黄片免| 动漫黄色视频在线观看| 欧美精品亚洲一区二区| 无限看片的www在线观看| 亚洲最大成人中文| 国产成人av激情在线播放| 色av中文字幕| x7x7x7水蜜桃| x7x7x7水蜜桃| 日本一区二区免费在线视频| 欧美日韩一级在线毛片| 好男人在线观看高清免费视频 | 亚洲久久久国产精品| 巨乳人妻的诱惑在线观看| 99国产极品粉嫩在线观看| 亚洲国产高清在线一区二区三 | 制服诱惑二区| 一区二区三区精品91| 人成视频在线观看免费观看| 麻豆av在线久日| 青草久久国产| 久久久久九九精品影院| 久久精品成人免费网站| 亚洲国产精品成人综合色| 热99re8久久精品国产| 午夜免费激情av| 女警被强在线播放| 午夜福利在线观看吧| 一个人免费在线观看的高清视频| 亚洲精品美女久久久久99蜜臀| 一区福利在线观看| 国产高清视频在线播放一区| 国产亚洲av高清不卡| 97人妻天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 欧美性长视频在线观看| 国产精品国产高清国产av| 女人爽到高潮嗷嗷叫在线视频| 国产激情欧美一区二区| 日日摸夜夜添夜夜添小说| 黑人操中国人逼视频| 日日爽夜夜爽网站| 精品国内亚洲2022精品成人| 亚洲色图综合在线观看| 日本在线视频免费播放| bbb黄色大片| 丝袜人妻中文字幕| 可以在线观看的亚洲视频| 久久久精品国产亚洲av高清涩受| 国产极品粉嫩免费观看在线| 久久久国产成人免费| 啪啪无遮挡十八禁网站| 亚洲精品久久国产高清桃花| 18美女黄网站色大片免费观看| 丝袜美腿诱惑在线| 亚洲九九香蕉| 亚洲性夜色夜夜综合| 国产一区二区三区综合在线观看| 一进一出好大好爽视频| 国产一区二区激情短视频| 欧美黑人精品巨大| 19禁男女啪啪无遮挡网站| 久久久久久国产a免费观看| av有码第一页| 日本a在线网址| 18禁裸乳无遮挡免费网站照片 | 亚洲少妇的诱惑av| 亚洲片人在线观看| 黄片小视频在线播放| 欧美黄色片欧美黄色片| 久久久久久久午夜电影| av在线天堂中文字幕| 午夜免费激情av| 性少妇av在线| 午夜成年电影在线免费观看| 啦啦啦观看免费观看视频高清 | 国产男靠女视频免费网站| 国产不卡一卡二| 一二三四社区在线视频社区8| 久久久久久久午夜电影| 男女之事视频高清在线观看| 91成年电影在线观看| 丝袜美腿诱惑在线| 制服诱惑二区| 免费在线观看黄色视频的| 99久久国产精品久久久| 久久久精品国产亚洲av高清涩受| 国产欧美日韩综合在线一区二区| x7x7x7水蜜桃| 亚洲男人的天堂狠狠| 亚洲天堂国产精品一区在线| 午夜免费鲁丝| netflix在线观看网站| 天天一区二区日本电影三级 | 国产主播在线观看一区二区| 精品乱码久久久久久99久播| 亚洲久久久国产精品| or卡值多少钱| 香蕉丝袜av| 欧美日韩亚洲综合一区二区三区_| 亚洲性夜色夜夜综合| 国产av精品麻豆| 精品久久久久久久毛片微露脸| 国产伦人伦偷精品视频| 满18在线观看网站| 中亚洲国语对白在线视频| 久久性视频一级片| 一级作爱视频免费观看| 日韩一卡2卡3卡4卡2021年| 色av中文字幕| 午夜成年电影在线免费观看| 午夜久久久在线观看| 成人18禁高潮啪啪吃奶动态图| 黄色丝袜av网址大全| 欧美国产日韩亚洲一区| 99久久99久久久精品蜜桃| 精品福利观看| 欧美日韩黄片免| 免费在线观看影片大全网站| 波多野结衣av一区二区av| 欧美激情 高清一区二区三区| 18禁国产床啪视频网站| 天堂动漫精品| 淫妇啪啪啪对白视频| 天堂影院成人在线观看| 韩国av一区二区三区四区| 老熟妇乱子伦视频在线观看| 99久久久亚洲精品蜜臀av| 国产欧美日韩一区二区三| 搡老妇女老女人老熟妇| 麻豆av在线久日| 国产成人欧美在线观看| 国产在线精品亚洲第一网站| 黄色a级毛片大全视频| 欧美 亚洲 国产 日韩一| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品久久久久久毛片| 欧美av亚洲av综合av国产av| 99国产精品免费福利视频| а√天堂www在线а√下载| 久久久国产成人免费| 老司机福利观看| 亚洲国产欧美网| 欧美成人午夜精品| 色在线成人网| 美女高潮到喷水免费观看| 日韩高清综合在线| 亚洲一区高清亚洲精品| 国产精品一区二区免费欧美| 女性被躁到高潮视频| 神马国产精品三级电影在线观看 | 午夜精品久久久久久毛片777| 色尼玛亚洲综合影院| 亚洲av片天天在线观看| 午夜两性在线视频| 亚洲欧美日韩无卡精品| 久久香蕉国产精品| 久久久久九九精品影院| 色综合欧美亚洲国产小说| 香蕉久久夜色| av在线播放免费不卡| 成人手机av| 最新美女视频免费是黄的| 国产精品久久久久久亚洲av鲁大| av网站免费在线观看视频| 国产日韩一区二区三区精品不卡| 变态另类成人亚洲欧美熟女 | 1024视频免费在线观看| 91麻豆av在线| 日韩国内少妇激情av| 久久欧美精品欧美久久欧美| 日本免费a在线| 日本在线视频免费播放| 久久香蕉精品热| av视频免费观看在线观看| 桃红色精品国产亚洲av| 看免费av毛片| 国产乱人伦免费视频| 国产xxxxx性猛交| 99国产综合亚洲精品| 99久久国产精品久久久| 欧美最黄视频在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 黄色女人牲交| 91av网站免费观看| 午夜精品在线福利| 国产成人啪精品午夜网站| 搡老熟女国产l中国老女人| 亚洲五月色婷婷综合| 成人三级黄色视频| 琪琪午夜伦伦电影理论片6080| 午夜精品在线福利| 黄色丝袜av网址大全| 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 国产欧美日韩精品亚洲av| 久久天躁狠狠躁夜夜2o2o| 精品午夜福利视频在线观看一区| 国产熟女xx| 自线自在国产av| 国产欧美日韩精品亚洲av| a在线观看视频网站| www.精华液| 久久久久国产精品人妻aⅴ院| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 后天国语完整版免费观看| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品一区二区三区| 国产三级在线视频| or卡值多少钱| 午夜老司机福利片| 一级毛片精品| 老熟妇仑乱视频hdxx| 天堂动漫精品| 国产精华一区二区三区| 欧美大码av| 久热爱精品视频在线9| 亚洲国产看品久久| 男女之事视频高清在线观看| 99久久精品国产亚洲精品| 99国产精品99久久久久| 悠悠久久av| 国产一区二区三区综合在线观看| 成人手机av| 级片在线观看| 国产精品99久久99久久久不卡| 在线视频色国产色| 国产av精品麻豆| 女警被强在线播放| 成人欧美大片| 午夜精品国产一区二区电影| 88av欧美| 香蕉国产在线看| 日日爽夜夜爽网站| 纯流量卡能插随身wifi吗| 19禁男女啪啪无遮挡网站| av福利片在线| 香蕉国产在线看| 两人在一起打扑克的视频| 午夜免费鲁丝| 国产成人欧美| 日韩欧美在线二视频| 国产亚洲精品综合一区在线观看 | 久久久久久国产a免费观看| 别揉我奶头~嗯~啊~动态视频| 国内精品久久久久精免费| 亚洲欧美激情在线| 国产精品,欧美在线| 久久精品91无色码中文字幕| 亚洲国产精品久久男人天堂| 国产乱人伦免费视频| 精品一品国产午夜福利视频| 人成视频在线观看免费观看| 亚洲精品国产一区二区精华液| 欧美绝顶高潮抽搐喷水| 在线视频色国产色| 欧美成人性av电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美不卡视频在线免费观看 | e午夜精品久久久久久久| 国产成人精品久久二区二区免费| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 精品无人区乱码1区二区| 看片在线看免费视频| 欧美一区二区精品小视频在线| 黄色视频不卡| 一进一出抽搐动态| netflix在线观看网站| 人人妻人人澡欧美一区二区 | 精品一区二区三区av网在线观看| 亚洲电影在线观看av| 一二三四社区在线视频社区8| 波多野结衣一区麻豆| 免费在线观看日本一区| ponron亚洲| 久久久久久大精品| 夜夜夜夜夜久久久久| 一本综合久久免费| 丁香欧美五月| 啪啪无遮挡十八禁网站| 国产精品永久免费网站| 一区二区三区激情视频| www.999成人在线观看| 日韩有码中文字幕| 国产免费男女视频| 精品人妻1区二区| 国产一卡二卡三卡精品| 久久精品成人免费网站| 亚洲欧洲精品一区二区精品久久久| 亚洲国产高清在线一区二区三 | 又黄又粗又硬又大视频| 男女下面进入的视频免费午夜 | 亚洲少妇的诱惑av| 国内毛片毛片毛片毛片毛片| x7x7x7水蜜桃| 搞女人的毛片| 久热这里只有精品99| 大香蕉久久成人网| 中文字幕av电影在线播放| 亚洲av成人不卡在线观看播放网| 91av网站免费观看| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区激情视频| 美女高潮到喷水免费观看| 十分钟在线观看高清视频www| 久久久久国产精品人妻aⅴ院| 国产主播在线观看一区二区| 91国产中文字幕| 看片在线看免费视频| 日韩成人在线观看一区二区三区| 高清黄色对白视频在线免费看| www.熟女人妻精品国产| 国产一级毛片七仙女欲春2 | 黄色女人牲交| 国产欧美日韩一区二区精品| 亚洲av片天天在线观看| 久久狼人影院| 如日韩欧美国产精品一区二区三区| 精品卡一卡二卡四卡免费| 男女下面进入的视频免费午夜 | 亚洲国产看品久久| 久久久久久久久久久久大奶| 级片在线观看| 亚洲国产高清在线一区二区三 | 动漫黄色视频在线观看| 成人亚洲精品一区在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 夜夜看夜夜爽夜夜摸| 亚洲人成77777在线视频| 成熟少妇高潮喷水视频| 久久精品91蜜桃| or卡值多少钱| 日韩欧美在线二视频| 免费在线观看影片大全网站| 天天添夜夜摸| 国产精品久久电影中文字幕| 国产成人欧美| 精品免费久久久久久久清纯| 国产精品爽爽va在线观看网站 | 丁香欧美五月| www.www免费av| 成人精品一区二区免费| 日本欧美视频一区| 两性夫妻黄色片| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 日日干狠狠操夜夜爽| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 亚洲精华国产精华精| 国产1区2区3区精品| 亚洲男人的天堂狠狠| 亚洲国产高清在线一区二区三 | 韩国av一区二区三区四区| 欧美在线黄色| 桃红色精品国产亚洲av| 91九色精品人成在线观看| 久久天堂一区二区三区四区| 国内精品久久久久久久电影| 精品少妇一区二区三区视频日本电影| 久久香蕉国产精品| 国产成人啪精品午夜网站| 麻豆久久精品国产亚洲av| 大陆偷拍与自拍| 欧洲精品卡2卡3卡4卡5卡区| 亚洲电影在线观看av| 久久久国产精品麻豆| 久久久久精品国产欧美久久久| 国产精品久久久久久亚洲av鲁大| 中亚洲国语对白在线视频| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片 | 亚洲成国产人片在线观看| 黑人巨大精品欧美一区二区mp4| 国产精品爽爽va在线观看网站 | 免费观看人在逋| 淫妇啪啪啪对白视频| 国产一区在线观看成人免费| 一进一出抽搐动态| 午夜激情av网站| 极品教师在线免费播放| 我的亚洲天堂| 精品久久久久久,| 麻豆国产av国片精品| 手机成人av网站| 淫秽高清视频在线观看| 日韩欧美国产在线观看| 婷婷精品国产亚洲av在线| 18禁裸乳无遮挡免费网站照片 | 好男人在线观看高清免费视频 | 成人精品一区二区免费| 免费av毛片视频| 99久久精品国产亚洲精品| 国产一区二区三区视频了| 免费搜索国产男女视频| 精品卡一卡二卡四卡免费| 亚洲久久久国产精品| 欧美中文综合在线视频| 久久人人精品亚洲av| 热99re8久久精品国产| 久久婷婷成人综合色麻豆| 久99久视频精品免费| 久久久久亚洲av毛片大全| 国产区一区二久久| 搡老岳熟女国产| 成人三级做爰电影| 国产男靠女视频免费网站| 99久久精品国产亚洲精品| 欧美亚洲日本最大视频资源| 最近最新免费中文字幕在线| 欧美乱色亚洲激情| 视频区欧美日本亚洲| 50天的宝宝边吃奶边哭怎么回事| 精品国产美女av久久久久小说| 国产激情久久老熟女| 亚洲国产精品999在线| 免费不卡黄色视频| 亚洲专区中文字幕在线| 如日韩欧美国产精品一区二区三区| 久久国产精品影院| 人成视频在线观看免费观看| 久久精品国产综合久久久| 久久久久久国产a免费观看| 一级毛片女人18水好多| 亚洲激情在线av| svipshipincom国产片| 亚洲精品久久成人aⅴ小说| 欧美成人免费av一区二区三区| 国产三级在线视频| 国产麻豆成人av免费视频| 国产一区二区三区视频了| 又大又爽又粗| 九色亚洲精品在线播放| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www| 国产午夜精品久久久久久| 夜夜躁狠狠躁天天躁| 国产精品免费一区二区三区在线| svipshipincom国产片| 国产不卡一卡二| 国产单亲对白刺激| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码| 少妇粗大呻吟视频| 国产av又大| 动漫黄色视频在线观看| 日韩高清综合在线| 亚洲欧美精品综合一区二区三区| 电影成人av| 亚洲午夜精品一区,二区,三区| 日本a在线网址| 搞女人的毛片| 国产亚洲精品久久久久5区| 亚洲国产日韩欧美精品在线观看 | 两个人看的免费小视频| 欧美黑人精品巨大| 日韩欧美在线二视频| 一区福利在线观看| 免费久久久久久久精品成人欧美视频| 在线视频色国产色| 欧美一级毛片孕妇| 日韩成人在线观看一区二区三区| 性少妇av在线| 国产男靠女视频免费网站| x7x7x7水蜜桃| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| 日韩免费av在线播放| 国产麻豆成人av免费视频| 18禁观看日本| av福利片在线| 好男人电影高清在线观看| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 别揉我奶头~嗯~啊~动态视频| 欧美一区二区精品小视频在线| 中文字幕人妻熟女乱码| 999久久久国产精品视频| 一进一出抽搐动态| 韩国精品一区二区三区| 久久久精品国产亚洲av高清涩受| 久久欧美精品欧美久久欧美| 老司机午夜十八禁免费视频| 国产亚洲精品第一综合不卡| 成人国语在线视频| 一级a爱片免费观看的视频| 色综合欧美亚洲国产小说| 日韩 欧美 亚洲 中文字幕| 欧美中文综合在线视频| 国产精品久久久久久精品电影 | 成人国语在线视频| 日本欧美视频一区| 国产精品永久免费网站| 88av欧美| 99国产精品99久久久久| 黑丝袜美女国产一区|