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

    基于地貌視角的喀斯特流域水系分維估算方法適應(yīng)性分析

    2019-12-02 01:13:54賀中華趙翠薇焦樹林楊朝暉
    自然資源遙感 2019年4期
    關(guān)鍵詞:分維漁網(wǎng)數(shù)據(jù)源

    安 全,賀中華,趙翠薇,梁 虹,焦樹林,楊朝暉

    (1.貴州師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,貴陽 550001;2.貴州師范大學(xué)國(guó)家喀斯特石漠化防治工程技術(shù)研究中心,貴陽 550001;3.貴州省山地資源與環(huán)境遙感應(yīng)用重點(diǎn)實(shí)驗(yàn)室,貴陽 550001;4.貴州省應(yīng)急管理廳,貴陽 550001)

    0 引言

    自然界許多系統(tǒng)、要素的空間組合結(jié)構(gòu)具有某種意義上的分形特點(diǎn)。諸多河流發(fā)育的支流與干流在空間結(jié)構(gòu)上呈現(xiàn)出一種典型的分形特征。通過傳統(tǒng)計(jì)算水系分維方法工作量大、費(fèi)時(shí)耗力,利用地理信息系統(tǒng)(geographic information system,GIS)技術(shù)并結(jié)合1∶5萬地形圖水系與30 m空間分辨率ASTER-GDEM數(shù)字高程數(shù)據(jù)估算流域水系分維,探討流域地貌發(fā)育與水文水系特征的關(guān)系是當(dāng)前研究的熱點(diǎn)問題。分形是自然界諸多客體表現(xiàn)局部特征與整體的相似性客觀幾何規(guī)律,在現(xiàn)代數(shù)學(xué)的發(fā)展中逐步演化成一支新的具有獨(dú)特理論體系與認(rèn)識(shí)的世界觀與方法論。在地學(xué)領(lǐng)域,地貌發(fā)育與河流水系的演化規(guī)律中本質(zhì)上都潛藏著分形的自相似性特點(diǎn)。分形理論在地學(xué)界的快速發(fā)展始于20世紀(jì)90年代,早在1977年Mandelbrot首先將分形理論引入水文學(xué),成為水系分維理論的創(chuàng)始人。進(jìn)入20世紀(jì)90年代后,國(guó)內(nèi)外對(duì)水系分維的研究迅速展開。何隆華等[1]利用計(jì)盒法對(duì)全國(guó)14個(gè)大流域與67個(gè)小流域進(jìn)行了計(jì)算,并得出水系分維值1.6是流域地貌發(fā)育階段的臨界值。進(jìn)入21世紀(jì)后,基于數(shù)字高程模型(digital elevation model,DEM)提取水系估算水系分維的研究得到快速發(fā)展[2-5]。往后隨著科學(xué)技術(shù)的快速發(fā)展,特別是GIS與DEM的發(fā)展與建立使得基于ArcGIS系列技術(shù)與DEM模型的流域水系分維研究進(jìn)入大發(fā)展時(shí)期[6-15]。如王倩[16]等利用GIS技術(shù),對(duì)秦淮河流域水系分維展開了研究,王林等[17]基于ArcGIS8.3和ArcVIEW3.2進(jìn)行了DEM模型的流域水系分維估算分析。

    目前,水系分維的研究主要是區(qū)域流域水系分維與地貌發(fā)育關(guān)系,人工渠系與自然水系分維對(duì)比,水系分維與洪澇、徑流過程,水系分維與巖性識(shí)別,水系分維與居名點(diǎn)分布特征,水系分維的生態(tài)意義等方面的研究[18-20]。而針對(duì)中國(guó)西南喀斯特筑壩流域地區(qū)與地貌發(fā)育對(duì)水文特征影響方面的研究報(bào)道還較少。鑒于此,本文以黔中筑壩工程區(qū)龍場(chǎng)橋流域?yàn)檠芯繉?duì)象,以1∶5萬地形圖水系、30 m空間分辨率的ASTER-GDEM和1∶10萬貴州省綜合地貌圖為數(shù)據(jù)源,利用基于ArcGIS10.2的Horton-Strahler理論、水系柵格法、漁網(wǎng)法估算黔中筑壩工程區(qū)龍場(chǎng)橋流域水系分維,探討筑壩區(qū)地貌發(fā)育對(duì)流域水文特征的影響,以期為喀斯特地區(qū)的生態(tài)建設(shè),經(jīng)濟(jì)的可持續(xù)發(fā)展,筑壩帶的洪澇風(fēng)險(xiǎn)評(píng)估以及流域綜合管理等方面提供一定參考。

    1 研究區(qū)概況及數(shù)據(jù)源

    1.1 研究區(qū)概況

    研究區(qū)是以龍場(chǎng)橋水文控制斷面提取的流域,位于黔中水利工程區(qū)上游地區(qū)(圖1)。黔中水利樞紐工程是貴州省首個(gè)大型跨地區(qū)、跨流域長(zhǎng)距離水利調(diào)水工程,也是黔中地區(qū)生存和發(fā)展的生命線工程[23]。工程區(qū)內(nèi)是長(zhǎng)江流域和珠江流域兩大流域交錯(cuò)地帶,研究區(qū)流域面積4 074.15 km2,地形西高東低,海拔在1 138~2 800 m之間,流域的西部地貌類型為典型的喀斯特巖溶高原[27],中東部峽谷、峰叢洼地相間分布,氣候?qū)儆趤啛釒Ъ撅L(fēng)氣候,夏季高溫多雨,冬季溫和濕潤(rùn),四季分明,多年平均降雨量在1 100~1 400 mm之間,最熱月平均氣溫22 ℃以上,最冷月0~15 ℃左右。植被為亞熱帶常綠闊葉林、灌木林、草地和裸巖相間分布。土壤以黃壤為主,兼有非地帶性石灰土分布。流域內(nèi)發(fā)育的烏江水系南源三岔河干流全長(zhǎng)325.6 km,流域集水面積7 264 km2,是黔中地區(qū)大型水利工程集中地。

    (a)龍場(chǎng)橋流域位置 (b)龍場(chǎng)橋流域DEM

    圖1 研究區(qū)概況圖

    Fig.1Researchareaoverview

    1.2 數(shù)據(jù)源及其預(yù)處理

    本文數(shù)據(jù)源包括1∶5萬地形圖水系(以下簡(jiǎn)稱地形圖水系)和30 m空間分辨率的ASTER-GDEM等。在對(duì)數(shù)據(jù)源的坐標(biāo)進(jìn)行一致性處理的基礎(chǔ)上,使用30 m空間分辨率 ASTER-GDEM數(shù)據(jù)進(jìn)行水系提取,利用ArcGIS10.2軟件的Hydrology水文分析工具集對(duì)從中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心的地理空間數(shù)據(jù)云(http://www.gscloud.cn)獲取的ASTER-GDEM數(shù)字高程數(shù)據(jù)進(jìn)行填洼、流向、匯流累積量計(jì)算處理,利用ArcGIS10.2軟件的Con命令設(shè)定匯流閾值1 500提取研究區(qū)河網(wǎng)柵格水系(以下簡(jiǎn)稱ASTER-GDEM水系)。

    2 研究方法

    利用ASTER-GDEM水系、結(jié)合0.6 m空間分辨率Google Earth影像、30 m空間分辨率Landsat 8 OLI遙感影像,對(duì)提取的筑壩水庫(kù)區(qū)平行狀水系依據(jù)Strahler理論進(jìn)行了分級(jí)。依次將提取水系柵格的邊長(zhǎng)設(shè)定為500~10 000 m,步長(zhǎng)為500 m;然后利用ArcGIS軟件的SQL查詢功能分別統(tǒng)計(jì)不同邊長(zhǎng)水系柵格所對(duì)應(yīng)的柵格數(shù)目,再對(duì)統(tǒng)計(jì)的柵格邊長(zhǎng)、數(shù)目與Strahler分級(jí)的矢量水系導(dǎo)入相關(guān)統(tǒng)計(jì)軟件,分別求柵格法與Horton定理的水系分維[28-29]。

    2.1 Horton定理

    Horton-Strahler計(jì)算原理為利用ArcGIS 10.2 提取1∶5萬地形圖水系,應(yīng)用Strahler分級(jí)原理對(duì)提取的水系進(jìn)行河流分級(jí);同理,以ASTER-GDEM為數(shù)據(jù)源,利用ArcGIS 10.2的con函數(shù)對(duì)30 m空間分辨率DEM分別設(shè)定匯流閾值100~10 000自動(dòng)提取水系,并與1∶5萬地形圖尺度下提取的水系不斷做疊加試驗(yàn),通過匯流閾值與河網(wǎng)密度的擬合曲線變點(diǎn)分析,確定用以提取研究區(qū)河網(wǎng)水系的最佳匯流閾值為1 500,再根據(jù)Strahler河網(wǎng)分級(jí)理論對(duì)矢量河網(wǎng)進(jìn)行分級(jí);最后統(tǒng)計(jì)地形圖水系和ASTER-GDEM通過最佳閾值1 500提取水系的每一級(jí)河流分級(jí)的河道數(shù)及其河道長(zhǎng)度,計(jì)算2個(gè)數(shù)據(jù)源提取水系的的分叉比與河長(zhǎng)比,對(duì)分叉比與河長(zhǎng)比取對(duì)數(shù),比值即為所求分維。計(jì)算公式如下:

    Rb=Ni-1/Ni,i=2,3,…,n,

    (1)

    (2)

    (3)

    式中:Rb為分叉比;RL為河長(zhǎng)比;Ni為各級(jí)河道數(shù);i為河道等級(jí);DH為Horton-Strahler水系分維值。

    2.2 柵格法

    柵格法又稱覆蓋法、網(wǎng)格法、計(jì)盒維數(shù)法,其原理是:取邊長(zhǎng)為r的正方形網(wǎng)格與水系圖求取交集,得到水系覆蓋的網(wǎng)格數(shù)目為N(r),當(dāng)r不斷變化時(shí),與之相對(duì)應(yīng)會(huì)得到一系列的N(r)值,二者關(guān)系為

    N(r)∝r-D。

    (4)

    左右兩邊求取對(duì)數(shù),以點(diǎn)(lgr,lgN(r))為坐標(biāo)作雙對(duì)數(shù)圖,采用最小二乘法可擬合出一條直線,即

    lgN(r)=-Dlgr+b,

    (5)

    式中:r為正方形網(wǎng)格的邊長(zhǎng);N(r)為對(duì)應(yīng)邊長(zhǎng)網(wǎng)格與水系圖求取交集所覆蓋的網(wǎng)格數(shù)目;b為待定系數(shù);D為雙對(duì)數(shù)曲線的斜率值,即所求的水系分維[25]。

    3 結(jié)果與分析

    3.1 黔中筑壩流域水系分維提取分析

    利用ArcGIS10.2軟件對(duì)提取水系進(jìn)行Strahler二次分級(jí),通過Horton-Strahler法求得地形圖水系與研究區(qū)ASTER-GDEM水系的水系分維,再通過水系柵格法、漁網(wǎng)法對(duì)提取地形圖水系與ASTER-GDEM水系,分別依次將柵格邊長(zhǎng)設(shè)定為500~10 000 m,步長(zhǎng)為500 m,分別統(tǒng)計(jì)不同柵格邊長(zhǎng)的柵格個(gè)數(shù),并求取柵格邊長(zhǎng)與不同邊長(zhǎng)柵格數(shù)的對(duì)數(shù),利用Origing9.1統(tǒng)計(jì)軟件作柵格邊長(zhǎng)與不同邊長(zhǎng)柵格數(shù)雙對(duì)數(shù)擬合線,求其斜率即為分維值,統(tǒng)計(jì)結(jié)果如表1—4所示。

    表1 基于地形圖水系的Horton-Strahler法參數(shù)Tab.1 Horton-Strahler method parameters based on water system of topographic map

    表2 基于ASTER-GDEM水系的Horton-Strahler法參數(shù)Tab.2 Horton-Strahler method parameters based on water system of ASTER-GDEM

    表1和表2的統(tǒng)計(jì)結(jié)果表明,河流數(shù)與河流平均長(zhǎng)度隨著河流級(jí)數(shù)的增大而減少。河流的分叉比隨著河流級(jí)數(shù)的增大而增大,與河流數(shù)、河流平均長(zhǎng)度成反比;河長(zhǎng)比隨河流級(jí)數(shù)越大而減少??傮w而言,河流分叉比、河長(zhǎng)比除了受到河流數(shù)及其河流長(zhǎng)度的影響,與河流的分級(jí)存在一定的相關(guān)性。

    表3和表4是地形圖水系與ASTER-GDEM水系柵格參數(shù)和漁網(wǎng)法參數(shù)。統(tǒng)計(jì)結(jié)果表明,柵格水系法與漁網(wǎng)法統(tǒng)計(jì)的柵格數(shù)目隨著試驗(yàn)設(shè)定的柵格邊長(zhǎng)的增加而減少,柵格邊長(zhǎng)與柵格數(shù)目之間存在明顯的線性關(guān)系。

    表3 地形圖水系與ASTER-GDEM水系柵格參數(shù)Tab.3 Water system of topographic map and ASTER-GDEM grid parameters

    表4 地形圖水系與ASTER-GDEM漁網(wǎng)法參數(shù)Tab.4 Water system of topographic map and ASTER-GDEM fishing net method parameters

    表5是結(jié)合表1至表4統(tǒng)計(jì)數(shù)據(jù),利用分維計(jì)算公式根據(jù)表1和表2分別求出1∶5萬地形圖提取水系和ASTER-GDEM提取水系的Horton-Strahler分維;利用Origing9.1統(tǒng)計(jì)軟件根據(jù)表3和表4分別作1∶5萬水系和ASTER-GDEM的柵格水系法及漁網(wǎng)水系法試驗(yàn)得到柵格邊長(zhǎng)與柵格數(shù)的雙對(duì)數(shù)擬合曲線。結(jié)果見表5。

    表5 不同方法提取地形圖水系與ASTER-GDEM分維值對(duì)比表Tab.5 Comparison of fractal values from water system of topographic map and ASTER- GDEM by different methods

    從表5中可以看出,同一地區(qū)不同方法、不同數(shù)據(jù)源估算的水系分維相差較大。其中Horton-Strahler法估算地形圖水系差別顯著,根據(jù)Horton-Strahler法估算地形圖水系與ASTER-GDEM水系的分維值分別為1.69和0.66,二者相差1.03。根據(jù)河數(shù)定律與河長(zhǎng)定律的分叉比、河長(zhǎng)比理論取值范圍可知,針對(duì)二元三維地表地下結(jié)構(gòu),溶洞、裂隙發(fā)育的典型喀斯特地貌而言,由于諸多干流沿途遇到溶洞,大型溶蝕管道、裂隙時(shí)會(huì)轉(zhuǎn)入地下流一段距離等特殊現(xiàn)象,使得在地形圖上提取水系時(shí)無法判斷地下的水流狀況導(dǎo)致提取水系誤差增大,根據(jù)地形圖水系計(jì)算的河流的分叉比、河長(zhǎng)比的比值范圍未在理論值范圍,以至于所計(jì)算出的水系分維值與水系柵格法、漁網(wǎng)法估算的分維值大相徑庭,不能穩(wěn)定的判斷研究區(qū)地貌發(fā)育狀況。與此同時(shí),表1和表2根據(jù)Horton-Strahler法提取的地形圖水系與ASTER-GDEM水系的河段數(shù)差別尤為顯著,總體上相差1 174段,河流分級(jí)的河段數(shù)差別主要體現(xiàn)在1~3級(jí)。這種顯著差異可能是由于復(fù)雜的喀斯特地貌組合結(jié)構(gòu)所致。

    然而,柵格水系法與漁網(wǎng)法估算的水系分維無論是地形圖水系還是ASTER-GDEM水系,估算的分維值均比較接近。就表5數(shù)據(jù)表明,通過ASTER-GDEM水系利用2種方法估算的分維略比地形圖水系大。另外,從水系提取的R2值來看,1∶5萬地形圖尺度下柵格水系法與漁網(wǎng)法提取水系分維的R2值分別為0.996 6和0.996 4,而基于ASTER-GDEM采用柵格水系法與漁網(wǎng)法提取水系分維的R2值分別為0.994 1和0.993 4,這表明前者提取水系的無標(biāo)度區(qū)間在0.996 4~0.996 6之間,而后者為0.993 4~0.994 1之間。

    3.2 黔中筑壩流域水系分維與地貌發(fā)育關(guān)系分析

    圖2為研究區(qū)根據(jù)1∶10萬貴州省綜合地貌圖矢量化的研究區(qū)地貌發(fā)育現(xiàn)狀圖和30 m空間分辨率ASTER-GDEM數(shù)字高程數(shù)據(jù)利用ArcGIS10.2制作的研究區(qū)山體陰影圖。表6為研究區(qū)各類型地貌類型及其面積比例。

    (a)研究區(qū)地貌類型 (b)研究區(qū)山體陰影

    圖2 研究區(qū)各流域地貌類型與山體陰影圖

    Fig.2Landformtypesandmountainshadowmapsofvariouswatershedsinthestudyarea

    表6 不同流域地貌發(fā)育參數(shù)及其面積比例Tab.6 Geomorphological development parameters and area ratio of different watersheds

    ①字母F,K1,K2分別表示地貌發(fā)育的侵蝕-剝蝕類型、溶蝕為主型、溶蝕-侵蝕類型。

    從圖2(a)中可以看出,整個(gè)研究區(qū)地貌發(fā)育類型以中山谷地、峰叢谷地、峰叢洼地為主,其中龍場(chǎng)橋流域、陽長(zhǎng)流域少部分地區(qū)出現(xiàn)深切中山型地貌。從圖2(b)山體陰影圖可以看出,研究區(qū)地貌發(fā)育破碎程度大,河谷多呈“V”型,越往研究區(qū)的東部即陽長(zhǎng)、龍場(chǎng)橋一帶,地勢(shì)起伏越大。地貌發(fā)育組合以深切中山、中山谷地、峰叢洼地復(fù)型組合為主。從表6中可以看出,向陽流域地貌發(fā)育組合以峰叢谷地、峰林溶原(盆地)、淺切中山為主,所占流域面積比例分別為23.31 %,20.83%和38.04%,研究區(qū)比例分別為4.76%,4.25%和7.77%;陽長(zhǎng)流域地貌發(fā)育組合以K化中山谷地、峰叢谷地、峰叢洼地為主,所占流域面積比例分別為26.75%,47.43%和11.62%,研究區(qū)比例分別為10.48%,18.59%和4.56%;龍場(chǎng)橋流域地貌發(fā)育以K化中山谷地、峰叢谷地、峰叢洼地為主,所占流域面積比例分別為43.23%,11.95%和41.14%,研究區(qū)比例分別為17.46%,4.83%和16.61%。

    3.3 不同流域地貌發(fā)育成因分析

    從表6中可以看出,向陽流域、陽長(zhǎng)流域和龍場(chǎng)橋流域3個(gè)流域的地貌發(fā)育成因主要為侵蝕-剝蝕類型,溶蝕為主類型、溶蝕-侵蝕類型,向陽流域由于侵蝕-剝蝕成因下發(fā)育的地貌類型深切中山面積占整個(gè)流域的38.04%,這說明整個(gè)流域地貌的發(fā)育主要以侵蝕-剝蝕為主,其次為溶蝕、溶蝕-侵蝕型。陽長(zhǎng)流域由于溶蝕、溶蝕-侵蝕成因下發(fā)育的地貌類型峰叢谷地面積占整個(gè)流域的47.43%,其地貌發(fā)育主要以溶蝕、溶蝕-侵蝕為主。龍場(chǎng)橋流域由于溶蝕、溶蝕-侵蝕成因下發(fā)育的地貌類型K化中山谷地、峰叢洼地面積占整個(gè)流域的34.07%,故其地貌發(fā)育主要以溶蝕、溶蝕-侵蝕為主。從地貌發(fā)育的角度看皆屬于地貌發(fā)育幼年期的中晚期階段。

    綜合水系估算分維值來看,根據(jù)根據(jù)何隆華的計(jì)算的水系分維與地貌發(fā)育階段的關(guān)系可知,當(dāng)D或DH≤1.6時(shí),流域地貌發(fā)育屬于侵蝕發(fā)育階段的幼年期;當(dāng)D>1.6或DH≤1.89時(shí),流域地貌發(fā)育屬于侵蝕發(fā)育階段的壯年期;當(dāng)D或DH>1.89時(shí),流域地貌發(fā)育屬于侵蝕發(fā)育階段的老年期。Horton-Strahler法、水系柵格法、漁網(wǎng)法提取除了Horton-Strahler法估算水系分維與研究區(qū)實(shí)際出入較大外,水系柵格法、漁網(wǎng)法估算的水系分維與研究區(qū)地貌發(fā)育較為吻合。結(jié)合表5和表6數(shù)據(jù)來看,綜合分析典型喀斯特復(fù)雜的地貌結(jié)構(gòu)與水系分維關(guān)系,漁網(wǎng)法估算的水系分維與研究區(qū)實(shí)際地貌現(xiàn)狀最為吻合。根據(jù)漁網(wǎng)估算的水系分維可知,研究區(qū)利用漁網(wǎng)法估算的1∶5萬地形圖提取水系分維值為1.54,通過ASTER-GDEM提取水系估算的分維值約為1.60,這說明研究區(qū)正處于地貌發(fā)育階段的幼年晚期、壯年期早期。此結(jié)果與研究區(qū)實(shí)際發(fā)育地貌吻合。

    4 結(jié)論與展望

    本文以黔中筑壩工程區(qū)龍場(chǎng)橋流域?yàn)檠芯繉?duì)象,以1∶5萬地形圖水系、30 m空間分辨率的ASTER-GDEM和1∶10萬貴州省綜合地貌圖為數(shù)據(jù)源,利用基于ArcGIS10.2的Horton-Strahler理論、水系柵格法、漁網(wǎng)法估算流域水系分維,探討筑壩區(qū)地貌發(fā)育對(duì)流域水文特征的影響,結(jié)果表明:

    1)喀斯特地區(qū)復(fù)雜地貌組合結(jié)構(gòu)下不同方法、不同數(shù)據(jù)源估算的水系分維相差較大。Horton-Strahler法、水系柵格法、漁網(wǎng)法估算地形圖水系分維值分別為1.69,1,53和1.54;估算ASTER-GDEM水系的分維值分別為0.66,1.59和1.60。其中Horton-Strahler法估算分維值差別顯著,差值達(dá)到1.03。

    2)綜合分析Horton-Strahler理論、水系柵格法、漁網(wǎng)法估算喀斯特筑壩區(qū)不同數(shù)據(jù)源水系分維與實(shí)際地貌發(fā)育的關(guān)系可知,漁網(wǎng)法估算的水系分維與研究區(qū)實(shí)際地貌現(xiàn)狀最為吻合。根據(jù)漁網(wǎng)估算的水系分維可知,研究區(qū)利用漁網(wǎng)法估算的地形圖水系分維值為1.54,ASTER-GDEM水系估算的分維值約為1.60,這說明研究區(qū)正處于地貌發(fā)育階段的幼年晚期、壯年期早期。此結(jié)果與研究區(qū)實(shí)際發(fā)育地貌吻合。此外,3種方法估算喀斯特筑壩流域的水系分維精度排序?yàn)椋簼O網(wǎng)法>水系柵格法>Horton-Strahler法。

    猜你喜歡
    分維漁網(wǎng)數(shù)據(jù)源
    木質(zhì)材料視覺感觀語義分維量化評(píng)價(jià)研究
    包裝工程(2023年16期)2023-08-25 11:36:32
    爸爸的漁網(wǎng)
    改進(jìn)的投影覆蓋方法對(duì)遼河河道粗糙床面分維量化研究
    透視漁網(wǎng)
    Web 大數(shù)據(jù)系統(tǒng)數(shù)據(jù)源選擇*
    “漁網(wǎng)”來了,快跑!
    基于不同網(wǎng)絡(luò)數(shù)據(jù)源的期刊評(píng)價(jià)研究
    獨(dú)創(chuàng)的漁網(wǎng)
    基于分形滲流模型的導(dǎo)電瀝青混凝土的分維計(jì)算
    基于真值發(fā)現(xiàn)的沖突數(shù)據(jù)源質(zhì)量評(píng)價(jià)算法
    av国产免费在线观看| 久久天堂一区二区三区四区| 99re在线观看精品视频| av在线播放免费不卡| 亚洲美女视频黄频| 天堂√8在线中文| 亚洲av电影在线进入| 搡老妇女老女人老熟妇| 免费在线观看黄色视频的| 九色国产91popny在线| 成年版毛片免费区| 在线永久观看黄色视频| 日韩大码丰满熟妇| 九色国产91popny在线| 在线a可以看的网站| 99re在线观看精品视频| 亚洲五月天丁香| 国产v大片淫在线免费观看| 在线观看舔阴道视频| 好男人电影高清在线观看| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 99久久综合精品五月天人人| 亚洲色图av天堂| 18禁黄网站禁片免费观看直播| 天天添夜夜摸| 亚洲欧美激情综合另类| 国产精品亚洲美女久久久| 国产精品久久久久久亚洲av鲁大| av在线天堂中文字幕| 中文字幕精品亚洲无线码一区| 成人一区二区视频在线观看| 久久久久久九九精品二区国产 | 国产探花在线观看一区二区| ponron亚洲| 亚洲熟妇中文字幕五十中出| 亚洲成人久久爱视频| 中文字幕久久专区| 此物有八面人人有两片| 欧美日韩中文字幕国产精品一区二区三区| 性色av乱码一区二区三区2| 日韩大码丰满熟妇| 黄色片一级片一级黄色片| 久久热在线av| 老熟妇仑乱视频hdxx| www.999成人在线观看| 日本一二三区视频观看| 制服人妻中文乱码| 悠悠久久av| 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清专用| 久久久久久久久中文| 97碰自拍视频| 91成年电影在线观看| 成年人黄色毛片网站| 国产v大片淫在线免费观看| 91麻豆av在线| 欧美色欧美亚洲另类二区| 黄色片一级片一级黄色片| 少妇的丰满在线观看| 最新美女视频免费是黄的| 亚洲一区中文字幕在线| 国产精品av视频在线免费观看| 国产精品久久久久久人妻精品电影| 久久久久国产一级毛片高清牌| 欧美 亚洲 国产 日韩一| 成人三级黄色视频| 免费人成视频x8x8入口观看| 国产精品久久电影中文字幕| 在线观看免费视频日本深夜| 久久精品影院6| 欧美成狂野欧美在线观看| 最近最新中文字幕大全电影3| 国产亚洲av高清不卡| 一区福利在线观看| 日韩高清综合在线| 久久久久久久久免费视频了| 久久精品人妻少妇| 一个人免费在线观看电影 | 欧美高清成人免费视频www| 国产三级中文精品| 亚洲七黄色美女视频| 精品日产1卡2卡| 老熟妇乱子伦视频在线观看| 精品国产美女av久久久久小说| 亚洲精华国产精华精| 国内久久婷婷六月综合欲色啪| 国产精品免费视频内射| 宅男免费午夜| 国产欧美日韩一区二区精品| 丰满人妻熟妇乱又伦精品不卡| 麻豆一二三区av精品| 日本熟妇午夜| 亚洲 国产 在线| 欧美日韩乱码在线| 午夜免费激情av| 亚洲欧美日韩高清在线视频| 女人高潮潮喷娇喘18禁视频| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 一级a爱片免费观看的视频| 成人午夜高清在线视频| 欧美一级a爱片免费观看看 | 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 国产精品久久久人人做人人爽| 高清毛片免费观看视频网站| 啦啦啦观看免费观看视频高清| 亚洲精品一卡2卡三卡4卡5卡| 老司机靠b影院| 每晚都被弄得嗷嗷叫到高潮| 欧美又色又爽又黄视频| 免费电影在线观看免费观看| 国产成人av教育| 欧美乱码精品一区二区三区| 国产精品久久久人人做人人爽| 亚洲成av人片免费观看| 国产熟女xx| 一级a爱片免费观看的视频| a级毛片在线看网站| 九色国产91popny在线| 久久久久久九九精品二区国产 | 欧美精品亚洲一区二区| 1024香蕉在线观看| 99在线视频只有这里精品首页| 午夜日韩欧美国产| 日韩 欧美 亚洲 中文字幕| 丝袜人妻中文字幕| 91九色精品人成在线观看| 99re在线观看精品视频| 成人手机av| 黄色毛片三级朝国网站| avwww免费| 国产成人影院久久av| 啦啦啦韩国在线观看视频| 久久欧美精品欧美久久欧美| 久久精品91无色码中文字幕| 中文字幕av在线有码专区| 亚洲欧美日韩无卡精品| 久久天躁狠狠躁夜夜2o2o| 两个人看的免费小视频| 久久国产乱子伦精品免费另类| 母亲3免费完整高清在线观看| 亚洲人成网站高清观看| 人人妻人人澡欧美一区二区| 黄片大片在线免费观看| 99国产精品99久久久久| 国产精品av视频在线免费观看| 久久午夜亚洲精品久久| 男男h啪啪无遮挡| 美女大奶头视频| 男女视频在线观看网站免费 | 免费观看人在逋| 特大巨黑吊av在线直播| 亚洲人成网站高清观看| 1024手机看黄色片| 亚洲人成电影免费在线| 一卡2卡三卡四卡精品乱码亚洲| 久久久精品大字幕| 亚洲国产精品合色在线| 午夜免费成人在线视频| 中文字幕人妻丝袜一区二区| 日本一二三区视频观看| 99热这里只有是精品50| 在线观看66精品国产| 一a级毛片在线观看| 亚洲av成人av| 18美女黄网站色大片免费观看| 亚洲欧美精品综合久久99| 久久久久免费精品人妻一区二区| 午夜免费成人在线视频| 亚洲国产欧美人成| 久久精品91无色码中文字幕| 88av欧美| 99国产精品一区二区蜜桃av| 99久久精品国产亚洲精品| 99精品在免费线老司机午夜| 国产精品一区二区精品视频观看| 天天躁狠狠躁夜夜躁狠狠躁| 99国产综合亚洲精品| 亚洲成av人片免费观看| 最近视频中文字幕2019在线8| 一个人免费在线观看的高清视频| 亚洲精品久久国产高清桃花| www.熟女人妻精品国产| 变态另类成人亚洲欧美熟女| 国产精品久久久久久精品电影| 欧美乱妇无乱码| 国产精品电影一区二区三区| 亚洲 欧美 日韩 在线 免费| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 国产精品一区二区免费欧美| 妹子高潮喷水视频| 欧美黑人欧美精品刺激| 99在线视频只有这里精品首页| 久久精品国产综合久久久| 午夜福利18| 我要搜黄色片| 桃色一区二区三区在线观看| 老司机深夜福利视频在线观看| 日韩 欧美 亚洲 中文字幕| 午夜福利免费观看在线| 最新在线观看一区二区三区| 91九色精品人成在线观看| 美女午夜性视频免费| а√天堂www在线а√下载| av片东京热男人的天堂| 亚洲一区二区三区不卡视频| 亚洲精品中文字幕在线视频| 九色成人免费人妻av| 欧美日韩福利视频一区二区| netflix在线观看网站| 久久婷婷成人综合色麻豆| 久久婷婷人人爽人人干人人爱| 午夜福利在线观看吧| 亚洲精品国产一区二区精华液| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 国产欧美日韩精品亚洲av| 欧美成人一区二区免费高清观看 | 男人的好看免费观看在线视频 | 婷婷六月久久综合丁香| 男女视频在线观看网站免费 | 久久精品国产亚洲av香蕉五月| 男女视频在线观看网站免费 | 亚洲精品久久成人aⅴ小说| 少妇裸体淫交视频免费看高清 | 国产精品一区二区精品视频观看| 亚洲激情在线av| 制服诱惑二区| 日日摸夜夜添夜夜添小说| 天天躁夜夜躁狠狠躁躁| 婷婷亚洲欧美| 人妻丰满熟妇av一区二区三区| 精品人妻1区二区| 88av欧美| 国产欧美日韩一区二区精品| 午夜成年电影在线免费观看| 不卡一级毛片| av有码第一页| 欧美成人免费av一区二区三区| 在线a可以看的网站| 久久精品影院6| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| 国产av在哪里看| 我要搜黄色片| 校园春色视频在线观看| 一二三四在线观看免费中文在| 国产精品免费一区二区三区在线| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 此物有八面人人有两片| 香蕉丝袜av| 熟女少妇亚洲综合色aaa.| 三级男女做爰猛烈吃奶摸视频| 国产91精品成人一区二区三区| 香蕉国产在线看| 亚洲av中文字字幕乱码综合| 精品电影一区二区在线| 给我免费播放毛片高清在线观看| 欧美精品啪啪一区二区三区| 欧美日韩乱码在线| 韩国av一区二区三区四区| 五月玫瑰六月丁香| 国产精品永久免费网站| 97人妻精品一区二区三区麻豆| 别揉我奶头~嗯~啊~动态视频| 三级男女做爰猛烈吃奶摸视频| 91av网站免费观看| 精品久久久久久久人妻蜜臀av| 精品国产美女av久久久久小说| 久久草成人影院| 久久久久久九九精品二区国产 | 欧美日韩精品网址| 久久精品国产亚洲av高清一级| 免费一级毛片在线播放高清视频| 男女床上黄色一级片免费看| 最新美女视频免费是黄的| 美女午夜性视频免费| 免费在线观看成人毛片| 97碰自拍视频| 国产成人一区二区三区免费视频网站| 悠悠久久av| 99riav亚洲国产免费| 欧美性长视频在线观看| 欧美性猛交黑人性爽| 两个人看的免费小视频| 在线免费观看的www视频| 91麻豆av在线| 亚洲成人中文字幕在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲美女黄片视频| 欧美丝袜亚洲另类 | 九色成人免费人妻av| 中文在线观看免费www的网站 | 久久久精品大字幕| 国产精品免费一区二区三区在线| 中文在线观看免费www的网站 | 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 最新美女视频免费是黄的| 日本 av在线| 少妇的丰满在线观看| 国产三级中文精品| 精品无人区乱码1区二区| 悠悠久久av| 在线观看免费日韩欧美大片| 精品国产亚洲在线| 香蕉国产在线看| 日日干狠狠操夜夜爽| 99国产精品99久久久久| 亚洲无线在线观看| 丁香六月欧美| 亚洲av成人av| 两个人看的免费小视频| 深夜精品福利| 淫妇啪啪啪对白视频| 国产欧美日韩一区二区三| 精品午夜福利视频在线观看一区| 亚洲成人中文字幕在线播放| 精品日产1卡2卡| 日韩欧美精品v在线| 午夜激情福利司机影院| 国产精品综合久久久久久久免费| 日本一二三区视频观看| 小说图片视频综合网站| 国模一区二区三区四区视频 | 久久婷婷人人爽人人干人人爱| 亚洲欧美精品综合久久99| 精品久久久久久,| 国产视频内射| 少妇人妻一区二区三区视频| 日本撒尿小便嘘嘘汇集6| 一进一出抽搐gif免费好疼| 一本一本综合久久| 老熟妇仑乱视频hdxx| 丝袜人妻中文字幕| 国产亚洲av高清不卡| 色精品久久人妻99蜜桃| 色老头精品视频在线观看| 免费观看人在逋| 国产午夜精品论理片| 国产精品精品国产色婷婷| 亚洲av中文字字幕乱码综合| 天天添夜夜摸| 国产免费男女视频| 日日干狠狠操夜夜爽| 天天一区二区日本电影三级| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美最黄视频在线播放免费| 亚洲精品粉嫩美女一区| 国产探花在线观看一区二区| 亚洲中文av在线| 国产99久久九九免费精品| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| 免费高清视频大片| 免费一级毛片在线播放高清视频| 精品久久久久久久毛片微露脸| 麻豆一二三区av精品| 亚洲国产欧美人成| 韩国av一区二区三区四区| 舔av片在线| 亚洲成av人片在线播放无| 最近视频中文字幕2019在线8| 9191精品国产免费久久| 欧美日韩亚洲综合一区二区三区_| 搡老熟女国产l中国老女人| 丝袜美腿诱惑在线| 在线十欧美十亚洲十日本专区| 日本a在线网址| 免费在线观看视频国产中文字幕亚洲| 国产成人精品无人区| 亚洲av五月六月丁香网| 中国美女看黄片| 国产亚洲精品久久久久久毛片| 久久久久国内视频| 亚洲中文字幕一区二区三区有码在线看 | 90打野战视频偷拍视频| 国产午夜精品论理片| 黄频高清免费视频| 国产精品一区二区三区四区免费观看 | 久久久国产欧美日韩av| 婷婷丁香在线五月| 久久精品国产亚洲av香蕉五月| 夜夜夜夜夜久久久久| 999精品在线视频| 欧美日韩乱码在线| 观看免费一级毛片| 男女午夜视频在线观看| 久久久久国内视频| 国产伦在线观看视频一区| 88av欧美| 欧美极品一区二区三区四区| 男人舔奶头视频| 91国产中文字幕| 欧美激情久久久久久爽电影| 日韩三级视频一区二区三区| 成人18禁在线播放| 亚洲国产精品999在线| 一进一出抽搐动态| 三级男女做爰猛烈吃奶摸视频| 日本免费a在线| 欧美一级a爱片免费观看看 | 国产久久久一区二区三区| 人人妻人人澡欧美一区二区| 精品熟女少妇八av免费久了| 成人特级黄色片久久久久久久| 国产成人精品久久二区二区91| 好看av亚洲va欧美ⅴa在| 99精品欧美一区二区三区四区| 久久久精品大字幕| 露出奶头的视频| 国产午夜精品论理片| 国产高清有码在线观看视频 | 18禁美女被吸乳视频| 免费无遮挡裸体视频| 久久精品成人免费网站| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 女人高潮潮喷娇喘18禁视频| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 成人精品一区二区免费| 欧美日韩乱码在线| 三级男女做爰猛烈吃奶摸视频| 欧美成狂野欧美在线观看| 男人舔女人的私密视频| 久久久久国产一级毛片高清牌| 草草在线视频免费看| 最新美女视频免费是黄的| 床上黄色一级片| 91九色精品人成在线观看| 成人午夜高清在线视频| 日本黄色视频三级网站网址| 中文字幕精品亚洲无线码一区| 国产在线精品亚洲第一网站| 国产av不卡久久| 91在线观看av| 老鸭窝网址在线观看| 国产视频一区二区在线看| 欧美人与性动交α欧美精品济南到| 日韩欧美国产一区二区入口| 精品午夜福利视频在线观看一区| 国产精品免费一区二区三区在线| 男男h啪啪无遮挡| 男女之事视频高清在线观看| 麻豆成人av在线观看| 青草久久国产| 国产av在哪里看| 亚洲精品国产一区二区精华液| 久久中文看片网| 成人欧美大片| 国产亚洲精品久久久久5区| 老熟妇乱子伦视频在线观看| 午夜精品在线福利| 男女午夜视频在线观看| 高清在线国产一区| 久久香蕉国产精品| 欧美精品啪啪一区二区三区| 毛片女人毛片| 国产亚洲精品第一综合不卡| 伦理电影免费视频| 亚洲av日韩精品久久久久久密| 日韩国内少妇激情av| www.www免费av| 久久久久久九九精品二区国产 | 国产精品,欧美在线| 亚洲 欧美一区二区三区| 欧美日韩亚洲综合一区二区三区_| 精品福利观看| 禁无遮挡网站| 欧美日本亚洲视频在线播放| www.熟女人妻精品国产| 琪琪午夜伦伦电影理论片6080| 亚洲va日本ⅴa欧美va伊人久久| www.www免费av| 伊人久久大香线蕉亚洲五| 特大巨黑吊av在线直播| 国产高清视频在线观看网站| 搡老熟女国产l中国老女人| av超薄肉色丝袜交足视频| 每晚都被弄得嗷嗷叫到高潮| 欧美一级毛片孕妇| 美女高潮喷水抽搐中文字幕| 女人爽到高潮嗷嗷叫在线视频| 免费搜索国产男女视频| 午夜视频精品福利| cao死你这个sao货| 身体一侧抽搐| 成年女人毛片免费观看观看9| 1024香蕉在线观看| 91在线观看av| www.www免费av| 波多野结衣高清无吗| 欧美另类亚洲清纯唯美| 变态另类丝袜制服| 亚洲一区中文字幕在线| 午夜精品久久久久久毛片777| 欧美av亚洲av综合av国产av| 观看免费一级毛片| 久久中文字幕一级| 亚洲国产高清在线一区二区三| 在线看三级毛片| 99久久精品热视频| 午夜日韩欧美国产| 在线观看午夜福利视频| 午夜免费成人在线视频| 天堂影院成人在线观看| 亚洲av电影在线进入| 国产精品99久久99久久久不卡| 精品熟女少妇八av免费久了| avwww免费| 丁香六月欧美| 久久亚洲真实| 日韩免费av在线播放| 色av中文字幕| 国产av不卡久久| 757午夜福利合集在线观看| 制服诱惑二区| 国产真人三级小视频在线观看| 国产精品一区二区免费欧美| 日本免费a在线| 亚洲精品在线美女| 久久精品国产99精品国产亚洲性色| 免费在线观看完整版高清| 欧美日韩国产亚洲二区| 毛片女人毛片| 巨乳人妻的诱惑在线观看| 国内精品久久久久精免费| 亚洲欧美一区二区三区黑人| 亚洲中文字幕一区二区三区有码在线看 | 亚洲avbb在线观看| 国产一区二区激情短视频| 亚洲国产精品合色在线| 久99久视频精品免费| 久久精品夜夜夜夜夜久久蜜豆 | 黄片小视频在线播放| 亚洲一区二区三区不卡视频| 亚洲成av人片免费观看| 欧美在线黄色| 看片在线看免费视频| 中亚洲国语对白在线视频| 神马国产精品三级电影在线观看 | 九九热线精品视视频播放| 成年版毛片免费区| 久久这里只有精品19| 草草在线视频免费看| 亚洲av片天天在线观看| 久久久久九九精品影院| 亚洲国产精品sss在线观看| 人妻丰满熟妇av一区二区三区| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 久久伊人香网站| 天天添夜夜摸| 亚洲中文字幕一区二区三区有码在线看 | 国产人伦9x9x在线观看| 国产精品av久久久久免费| 精品国产亚洲在线| 91麻豆av在线| 琪琪午夜伦伦电影理论片6080| 亚洲专区国产一区二区| 亚洲免费av在线视频| 波多野结衣高清作品| 亚洲午夜精品一区,二区,三区| 国产成人av教育| 国产私拍福利视频在线观看| 欧美绝顶高潮抽搐喷水| 精品欧美国产一区二区三| 亚洲国产看品久久| 人人妻人人看人人澡| 黄片大片在线免费观看| 国产真人三级小视频在线观看| 国产激情偷乱视频一区二区| 久久久久国内视频| 少妇人妻一区二区三区视频| 成人欧美大片| 久久精品成人免费网站| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片精品| 欧美成人一区二区免费高清观看 | 两个人视频免费观看高清| 91字幕亚洲| 国产高清videossex| 丝袜美腿诱惑在线| 99久久久亚洲精品蜜臀av| 丰满人妻熟妇乱又伦精品不卡| 亚洲狠狠婷婷综合久久图片| 999久久久精品免费观看国产| 亚洲精品国产一区二区精华液| 亚洲狠狠婷婷综合久久图片| 国产一区二区三区视频了| 欧美中文综合在线视频| 亚洲人成77777在线视频| 一区二区三区激情视频| 成人亚洲精品av一区二区| 天天添夜夜摸| 每晚都被弄得嗷嗷叫到高潮| 伦理电影免费视频| 亚洲七黄色美女视频| 国产三级中文精品| 日本三级黄在线观看| АⅤ资源中文在线天堂| 中文资源天堂在线| 久久久久久久精品吃奶| 黑人巨大精品欧美一区二区mp4| 毛片女人毛片| 日韩欧美精品v在线| xxx96com| 国内毛片毛片毛片毛片毛片|