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

    基于反地形坡長(zhǎng)地表匯流模擬逼近算法

    2023-10-09 08:43:10桑云云羅明良白雷超姜佳艷劉維明
    安徽農(nóng)學(xué)通報(bào) 2023年17期
    關(guān)鍵詞:樣區(qū)坡長(zhǎng)延川

    桑云云 羅明良 楊 慧 白雷超 姜佳艷 劉維明

    (1西華師范大學(xué)生命科學(xué)學(xué)院,四川南充 637009;2西華師范大學(xué)地理科學(xué)學(xué)院,四川南充 637009;3西華師范大學(xué)四川省干旱河谷土壤侵蝕檢測(cè)與控制工程實(shí)驗(yàn)室,四川南充 637009;4中國(guó)科學(xué)院水利部、成都山地災(zāi)害與環(huán)境研究所,四川成都 610000)

    我國(guó)土壤侵蝕分布范圍較廣,是土壤侵蝕較為嚴(yán)重的國(guó)家之一[1-2]。坡長(zhǎng)是土壤侵蝕領(lǐng)域研究中的重要因子之一[3],當(dāng)其他條件相同時(shí),土壤侵蝕的強(qiáng)度由坡面的長(zhǎng)度來決定,坡面越長(zhǎng),匯聚的流量越大,侵蝕力和沖刷力就越強(qiáng)[4-5]。目前,關(guān)于坡度因子的計(jì)算已經(jīng)相對(duì)成熟,而坡長(zhǎng)因子的精確、快速提取是目前研究的重點(diǎn)之一[6-7]。

    對(duì)于坡長(zhǎng)的計(jì)算,學(xué)者提出了多種方法,例如:Hickey 等[8]在水文分析的基礎(chǔ)上提出了非累積流量的直接算法;Desmet 等[9]提出了基于累積流量的單位匯水面積算法;楊昕等[10]基于Arcview空間分析提出坡長(zhǎng)快速近似算法。隨著計(jì)算機(jī)和GIS 研究方法的不斷創(chuàng)新,在DEM 數(shù)據(jù)和GIS 空間分析方法的支持下,坡長(zhǎng)提取的方法得以不斷改進(jìn)。其中,非累積流量的直接算法可以通過按照不同的流向進(jìn)行賦值來體現(xiàn)坡長(zhǎng)的準(zhǔn)確特征[11],但此方法需要大量的計(jì)算和復(fù)雜地形因子的提取;快速近似算法雖然計(jì)算過程簡(jiǎn)便,但由于地形類型復(fù)雜,局部高程點(diǎn)無法準(zhǔn)確識(shí)別,導(dǎo)致計(jì)算結(jié)果誤差較大[8]。基于網(wǎng)格DEM 的坡長(zhǎng)計(jì)算可用來模擬地表的水流路徑,但要求水流在DEM 模擬的地形面上自由流動(dòng),因此地表的自然凹陷對(duì)坡長(zhǎng)的計(jì)算結(jié)果有很大影響,需要在確定水流方向之前對(duì)DEM凹陷處進(jìn)行填充和平整。此外,坡長(zhǎng)的計(jì)算還需要考慮很多因素,不僅包括地形特征,還需要考慮具體地塊的土地利用類型[12-13]。坡長(zhǎng)作為土壤侵蝕計(jì)算中不可忽視的因子,迫切需要一個(gè)高度精確的計(jì)算方法。

    為此,本研究以黃土高原甘泉和延川小流域?yàn)檠芯繀^(qū),基于5 m分辨率數(shù)據(jù),提出在反地形條件下不斷逼近真實(shí)坡長(zhǎng)結(jié)果的算法。該研究不僅降低了原有坡長(zhǎng)計(jì)算方法中由于填洼造成的誤差,也有效減輕了計(jì)算過程中的工作量,為坡長(zhǎng)計(jì)算提供了一定的思路和借鑒。

    1 研究區(qū)概況

    本研究選取黃土高原侵蝕劇烈、地形起伏大的甘泉和延川2個(gè)典型的丘陵溝壑地貌區(qū)作為研究對(duì)象[14-15]。研究區(qū)海拔為600~1 800 m,甘泉位于陜西省延安市中部,屬于陜西省北部黃土高原黃土丘陵,以黃土梁狀丘陵溝壑地貌為主,溝壑密度約5.6 km/100 hm2。延川則以梁峁?fàn)钋鹆甑孛矠橹鳎横蛊鸱?,殘塬、梁、峁相間分布,溝壑縱橫、河谷深切,溝壑密度在所選樣區(qū)中最大。研究區(qū)內(nèi)以大陸性季風(fēng)氣候?yàn)橹?,大陸性氣候特征明顯,年均降水量約447 mm,多短時(shí)強(qiáng)降雨,土壤侵蝕強(qiáng)烈,如圖1所示。

    圖1 研究區(qū)地形

    2 研究方法

    2.1 研究原理

    坡長(zhǎng)一般定義為從坡面徑流的起點(diǎn)到徑流被攔截點(diǎn)或流路中斷點(diǎn)的水平距離[8,16]。如圖2 中A′B′是徑流線在水平面上的投影,即坡長(zhǎng)。

    圖2 坡長(zhǎng)示意

    本研究在反地形條件下計(jì)算坡長(zhǎng)時(shí),為確定水流的起點(diǎn),將原始地形分解為若干個(gè)集水盆地。當(dāng)反轉(zhuǎn)地形時(shí),原始正向地貌的山地反轉(zhuǎn)成為負(fù)向流域及集水區(qū)系統(tǒng),相應(yīng)的山頂點(diǎn)也演化為各集水區(qū)中的洼地點(diǎn),原有正向地形上的洼地成為向上凸起的小丘。反地形水流的起點(diǎn)在原始地形中集水盆地的凹陷點(diǎn)上。該方法雖省略了填洼過程,但并不會(huì)對(duì)原始地形流路和相關(guān)數(shù)據(jù)的提取產(chǎn)生影響。這種計(jì)算過程不僅簡(jiǎn)單易懂,而且可以減小多種數(shù)據(jù)融合后計(jì)算過程中產(chǎn)生的誤差。

    根據(jù)坡長(zhǎng)的定義,如果將坡上的水流路徑(M)設(shè)為直角三角形的斜邊,則坡長(zhǎng)可定義為坡上流線的長(zhǎng)度乘以cos(α),其中α為流向與水平方向的夾角[17-19]。為計(jì)算α,坡長(zhǎng)逼近算法(ASL)首先計(jì)算從水流流動(dòng)起點(diǎn)到流路上某一點(diǎn)的高差,定義為δh,然后用反正弦函數(shù)求解α,即α=arcsin(δh/M)。按照這種思路,坡長(zhǎng)逼近算法ASL的設(shè)計(jì)步驟如下。

    反地形和相關(guān)水文分析(圖3),如流向、匯流累積量和分水嶺等[20]。其中,默認(rèn)的流向算法為D8流向算法,即單元梯度的最佳代表值是以其為中心的3×3窗口中周圍8個(gè)方向的最大梯度,水流的流動(dòng)方向與最大梯度相同。

    圖3 水文分析

    分區(qū)統(tǒng)計(jì)后利用柵格計(jì)算器計(jì)算前一步分區(qū)統(tǒng)計(jì)結(jié)果與原始DEM的差值,利用重分類工具將上一步柵格計(jì)算的結(jié)果進(jìn)行賦值顯示后,復(fù)制為1 的區(qū)域即為區(qū)域最高點(diǎn)。接著使用水文分析中的水流長(zhǎng)度工具基于所得到的流向數(shù)據(jù)計(jì)算出反地形下的水流長(zhǎng)度(上坡長(zhǎng)度:計(jì)算沿水流路徑從每個(gè)像元到分水嶺頂部的最長(zhǎng)上坡距離),即為地形坡面長(zhǎng)度。

    利用網(wǎng)格計(jì)算器構(gòu)建地圖表達(dá)式α=arcsin(δh/M),并計(jì)算坡面與水平面之間的包含角α;使用三角函數(shù)cos構(gòu)建函數(shù)表達(dá)式L=M×cos(α),計(jì)算結(jié)果L即坡長(zhǎng)(圖4)。

    圖4 坡長(zhǎng)計(jì)算示意

    2.2 坡長(zhǎng)逼近算法驗(yàn)證

    為驗(yàn)證ASL計(jì)算坡長(zhǎng)的精度,以相對(duì)差值系數(shù)和線性函數(shù)為描述指標(biāo),對(duì)3種坡長(zhǎng)計(jì)算結(jié)果進(jìn)行分析。

    2.2.1 LS_TOOL(LST)LS_TOOL 算法是基于GIS 技術(shù)和多方位地形系數(shù)的估算方法,用于獲取區(qū)域尺度的地形坡長(zhǎng)坡度參數(shù)等數(shù)據(jù)。以LS_TOOL 算法計(jì)算的坡長(zhǎng)(L)作為實(shí)際坡長(zhǎng),驗(yàn)證本文提出的近似算法的可靠性[18]。

    2.2.2 SAGA(SGL)SAGA軟件中的坡長(zhǎng)計(jì)算工具可以計(jì)算特殊情況下的坡長(zhǎng)。首先用Zevenbergen &Thorne 方法計(jì)算每個(gè)單元的斜率,然后使用D8單流向算法和閾值標(biāo)準(zhǔn),將坡長(zhǎng)值向下累積。

    2.2.3 相關(guān)性分析 利用線性回歸方法擬合出算法間的一元線性函數(shù),當(dāng)2 種算法之間的坡長(zhǎng)值完全一致時(shí),則函數(shù)關(guān)系為X=Y,XY 散點(diǎn)分布圖是一條斜率為1的直線;當(dāng)2種算法之間的坡長(zhǎng)值存在差異時(shí),另一種函數(shù)關(guān)系為Y=a+bx,表明2種算法之間存在系統(tǒng)差異。接著計(jì)算相對(duì)差異系數(shù),相對(duì)差異系數(shù)通過逐個(gè)柵格單元進(jìn)行對(duì)比,可分區(qū)域分析算法之間的相似性。當(dāng)2 種算法進(jìn)行比較時(shí),選擇其中一個(gè)作為基準(zhǔn)算法,另外一個(gè)則作為比較算法,它們之間的差值θ如公式 (1)所示,可以通過減去相應(yīng)位置的結(jié)果而得到(是基準(zhǔn)算法在網(wǎng)格單元i的坡長(zhǎng)為整個(gè)DEM的平均坡長(zhǎng),是比較算法在相應(yīng)位置的坡長(zhǎng)。)

    為能夠更清楚地比較坡長(zhǎng)算法之間的相關(guān)關(guān)系,將其轉(zhuǎn)換為相對(duì)差異系數(shù)α,表示比較算法和基準(zhǔn)算法之間的偏差,其公式如下:

    當(dāng)α=1 時(shí),2 種計(jì)算方法的結(jié)果完全一致,沒有差異;α值越小,則算法之間的差異越大;當(dāng)α<0時(shí),則2種算法基本沒有可比性。

    3 結(jié)果與分析

    3.1 不同區(qū)域3種坡長(zhǎng)算法

    ASL、LST和SGL 3種計(jì)算方法得到的甘泉和延川樣區(qū)的坡長(zhǎng)結(jié)果如圖5~6 所示。從圖5~6 可以看出,同一樣區(qū)3種算法的坡長(zhǎng)變化大致相同,但坡長(zhǎng)結(jié)果均有差異,主要表現(xiàn)在3 種坡長(zhǎng)算法的最大值上,ASL算法的最大值分別為564和468,LST算法的最大值為493 和353,SGL 算法的最大值為399 和389,ASL 算法相對(duì)于其他2 種算法較大的坡長(zhǎng)值主要集中在山脊處。

    圖6 延川樣區(qū)坡長(zhǎng)(aASL、bLST、cSGL)

    經(jīng)過抽稀提取2 個(gè)研究區(qū)20%的坡長(zhǎng)結(jié)果,并將結(jié)果劃分為3個(gè)部分繪制累積坡長(zhǎng)(圖7)。從圖7可以看出,3種算法得到的坡長(zhǎng)在0~300 m,占提取坡長(zhǎng)總數(shù)的98%左右;而大于300 m的坡長(zhǎng)較少,占2%左右。ASL受流向和網(wǎng)格分辨率的影響,在甘泉樣區(qū)有超過500 m的離群值,但其累計(jì)數(shù)量不大,僅占坡長(zhǎng)總數(shù)的0.001 7%;而延川樣區(qū)超過400 m 的離群值占總數(shù)的0.002 4%,這些離群值均出現(xiàn)在山脊處。

    圖7 坡長(zhǎng)累積頻率分布(a甘泉坡長(zhǎng)結(jié)果頻率分布,b延川坡長(zhǎng)結(jié)果頻率分布)

    3.2 相關(guān)性分析

    為檢驗(yàn)ASL算法的可靠性,利用ArcGIS 10.2中的抽稀工具提取了3 種算法坡長(zhǎng)結(jié)果20%的數(shù)據(jù),并分別計(jì)算ASL 與LST、ASL 與SGL、SGL 與LST 算法之間的相關(guān)關(guān)系。

    3.2.1 線性回歸分析 每2種坡長(zhǎng)算法間的一元線性函數(shù)以及R2值見表1。

    表1 3種算法的相關(guān)性分析結(jié)果

    算法之間的線性函數(shù)的R2值均在0.900 以上,其中最大的R2值出現(xiàn)在延川樣區(qū)的ASL 與LST 之間,為0.993;最小R2值出現(xiàn)在延川樣區(qū)的SGL 與LST之間,為0.937,表明3種算法的坡長(zhǎng)結(jié)果之間存在強(qiáng)相關(guān)性。

    3.2.2 相對(duì)差異系數(shù) 利用抽稀得到的2個(gè)樣區(qū)20%坡長(zhǎng)數(shù)據(jù)計(jì)算其相對(duì)差異系數(shù)(a,b),結(jié)果見表2。

    表2 3種算法的相對(duì)差異系數(shù)

    從表2 可以看出,算法間的相對(duì)差異系數(shù)均在0.80 以上,說明每2 個(gè)樣區(qū)的坡長(zhǎng)結(jié)果高度相關(guān)。其中,甘泉區(qū)ASL 與LST 的相對(duì)差異系數(shù)最高,為0.96;延川樣區(qū)ASL與SGL的相對(duì)差異系數(shù)最低,為0.80。說明整體上這2 個(gè)樣區(qū)的坡長(zhǎng)結(jié)果相似度很高,但不同樣區(qū)算法間的相對(duì)差異系數(shù)不完全相同,表示算法之間存在一定的差異。

    4 討論和結(jié)論

    4.1 討論

    影響坡長(zhǎng)計(jì)算結(jié)果并產(chǎn)生差異的原因,首先是地形表達(dá)的有限分辨率限制了地形高度的表達(dá)精度[21-22]。受網(wǎng)格大小的限制,一些山峰或山脊表達(dá)為平坦地形,導(dǎo)致集水區(qū)的劃定結(jié)果出現(xiàn)一定誤差。由圖8 中的等高線可以看出,地形變化比較明顯,頂部有大面積的平坦區(qū)域。結(jié)合谷歌圖片,可以看出坡長(zhǎng)出現(xiàn)較大值時(shí)存在大量的梯田、黃土脊?fàn)钌角鸬绕教沟孛差愋停匦蜗鄬?duì)破碎,這對(duì)DEM反映真實(shí)地形帶來了較大挑戰(zhàn),進(jìn)而影響了坡長(zhǎng)的計(jì)算。如樣本區(qū)有較大的平坦區(qū)域,則會(huì)影響反地形下的分割精度,從而導(dǎo)致坡長(zhǎng)計(jì)算結(jié)果誤差較大[23]。本研究通過反地形計(jì)算坡長(zhǎng),去除了在原始地形中填充洼地帶來的誤差,同時(shí)利用函數(shù)方法計(jì)算坡長(zhǎng)提高了坡長(zhǎng)計(jì)算的準(zhǔn)確度,降低了試驗(yàn)過程中可能帶來的誤差,同時(shí)也減輕了計(jì)算過程中的工作量。

    圖8 研究區(qū)坡度和等高線(a甘泉,b延川)

    其次,D8流向算法的局限性可能會(huì)影響計(jì)算結(jié)果。流向是8 個(gè)可能的網(wǎng)格方向之一,它們之間的間隔為45°,流向是確定的且不連續(xù)的。因此,D8算法具有高度敏感性,即較小的高度誤差也會(huì)導(dǎo)致水流方向的改變,進(jìn)而導(dǎo)致坡長(zhǎng)計(jì)算結(jié)果出現(xiàn)偏差。其他2種坡長(zhǎng)算法同樣采用D8流向算法,因此不同流向算法對(duì)坡長(zhǎng)的具體影響,還需要在后續(xù)的研究中做進(jìn)一步的分析。

    4.2 結(jié)論

    本研究基于5 m 分辨率高程數(shù)據(jù)、利用ASL、LST、SGL 3種坡長(zhǎng)計(jì)算方法對(duì)黃土高原丘陵溝壑區(qū)2 個(gè)樣區(qū)的坡長(zhǎng)進(jìn)行計(jì)算,并利用線性回歸和相對(duì)差異系數(shù)對(duì)計(jì)算結(jié)果進(jìn)行了相似性分析,得出以下結(jié)論。

    相同精度的DEM數(shù)據(jù)得到2個(gè)樣區(qū)內(nèi)不同計(jì)算方法結(jié)果具有差異。本研究提出的逼近算法其坡長(zhǎng)結(jié)果98%的數(shù)據(jù)主要集中在0~300 m,并且線性回歸以及相對(duì)差異系數(shù)分析結(jié)果表明3種算法計(jì)算得到的坡長(zhǎng)結(jié)果非常相似。3個(gè)坡長(zhǎng)計(jì)算方法通過回歸函數(shù)得到的R2值均在0.9 以上,同時(shí)2 個(gè)樣區(qū)的3個(gè)算法間的相對(duì)差異系數(shù)也主要集中在0.9附近,表明3 種坡長(zhǎng)計(jì)算結(jié)果擬合效果均較好。因此,為獲取更高精度的坡長(zhǎng)結(jié)果,可將3 種算法的坡長(zhǎng)平均值作為最逼近真實(shí)的坡長(zhǎng)用于分析地區(qū)的土壤侵蝕狀況。

    坡長(zhǎng)因子作為估計(jì)土壤侵蝕狀況的重要參數(shù),尋求一個(gè)不斷逼近真實(shí)坡長(zhǎng)結(jié)果的算法對(duì)提高土壤侵蝕治理水平十分重要。本研究提出的坡長(zhǎng)逼近算法為坡長(zhǎng)計(jì)算提供了更多的思路。

    猜你喜歡
    樣區(qū)坡長(zhǎng)延川
    促進(jìn)大果沙棗扦插育苗生長(zhǎng)的最佳施肥措施
    桂林市銀杏綠化調(diào)查與分析
    地邊截水地物對(duì)黑土區(qū)小流域坡長(zhǎng)因子計(jì)算的影響
    野生植物對(duì)陜北黃土丘陵區(qū)土壤石油污染影響研究
    桂北油茶早實(shí)豐產(chǎn)林營(yíng)建現(xiàn)狀調(diào)查
    坡長(zhǎng)對(duì)貴州喀斯特區(qū)黃壤坡耕地土壤侵蝕的影響
    美麗的延川
    “井工廠”鉆井技術(shù)在延川南煤層氣開發(fā)中的應(yīng)用
    延川南區(qū)塊煤層氣鉆井防漏堵漏技術(shù)優(yōu)化
    延安至延川高速公路即將建成
    日本欧美视频一区| 一级作爱视频免费观看| 啦啦啦免费观看视频1| 成人三级黄色视频| 日日夜夜操网爽| 色综合欧美亚洲国产小说| 一区在线观看完整版| 母亲3免费完整高清在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品精品国产色婷婷| 女性生殖器流出的白浆| 亚洲精品粉嫩美女一区| 18禁国产床啪视频网站| 亚洲激情在线av| 亚洲,欧美精品.| 久久狼人影院| 国产亚洲av高清不卡| 久久久久久免费高清国产稀缺| 麻豆久久精品国产亚洲av| 国产欧美日韩一区二区三区在线| av中文乱码字幕在线| 亚洲片人在线观看| 亚洲中文av在线| 搡老妇女老女人老熟妇| 亚洲人成电影观看| 精品熟女少妇八av免费久了| 欧美黄色淫秽网站| 操出白浆在线播放| 久久国产乱子伦精品免费另类| 国产精品亚洲一级av第二区| 国产精品一区二区精品视频观看| 日韩中文字幕欧美一区二区| 韩国精品一区二区三区| 亚洲片人在线观看| 欧美激情极品国产一区二区三区| 一夜夜www| 欧美成狂野欧美在线观看| 日本vs欧美在线观看视频| 在线播放国产精品三级| 成人精品一区二区免费| 国产精品免费视频内射| 黑人巨大精品欧美一区二区蜜桃| 国产av一区二区精品久久| 女生性感内裤真人,穿戴方法视频| 大型av网站在线播放| 制服诱惑二区| 久久久久国产一级毛片高清牌| 中文字幕久久专区| 免费在线观看视频国产中文字幕亚洲| 欧美成人午夜精品| 一二三四在线观看免费中文在| АⅤ资源中文在线天堂| 国产精品爽爽va在线观看网站 | 日韩精品青青久久久久久| av视频免费观看在线观看| 黑丝袜美女国产一区| 日韩大码丰满熟妇| 国产亚洲精品av在线| 九色亚洲精品在线播放| 亚洲国产欧美日韩在线播放| 亚洲成人久久性| 久久国产乱子伦精品免费另类| 99热只有精品国产| 亚洲国产精品成人综合色| 美女 人体艺术 gogo| 久久伊人香网站| 午夜日韩欧美国产| 国产99白浆流出| 日韩 欧美 亚洲 中文字幕| 日本免费a在线| 久久久国产成人免费| 女人被狂操c到高潮| 亚洲国产中文字幕在线视频| 亚洲色图综合在线观看| 久久国产亚洲av麻豆专区| 国产精品一区二区精品视频观看| 欧美亚洲日本最大视频资源| 一级a爱片免费观看的视频| 亚洲第一青青草原| 黄频高清免费视频| 精品熟女少妇八av免费久了| 国产激情久久老熟女| 精品久久久久久久人妻蜜臀av | 男人舔女人下体高潮全视频| 女性生殖器流出的白浆| 满18在线观看网站| 99香蕉大伊视频| 亚洲av电影不卡..在线观看| av电影中文网址| 女人被躁到高潮嗷嗷叫费观| 亚洲视频免费观看视频| 午夜亚洲福利在线播放| 狠狠狠狠99中文字幕| 极品人妻少妇av视频| 久久精品国产综合久久久| 欧美日本亚洲视频在线播放| 亚洲精华国产精华精| 国产高清videossex| 一边摸一边做爽爽视频免费| 亚洲熟妇熟女久久| 亚洲成国产人片在线观看| 成人国产综合亚洲| 18禁国产床啪视频网站| 亚洲欧美激情综合另类| 曰老女人黄片| 国产黄a三级三级三级人| 女警被强在线播放| 女性被躁到高潮视频| 亚洲精品久久成人aⅴ小说| 久久久久国产一级毛片高清牌| 欧美日韩瑟瑟在线播放| 欧美黑人欧美精品刺激| 18禁黄网站禁片午夜丰满| 亚洲av成人一区二区三| 男女做爰动态图高潮gif福利片 | 中文字幕精品免费在线观看视频| 午夜免费观看网址| 免费不卡黄色视频| 变态另类成人亚洲欧美熟女 | 在线观看舔阴道视频| 丝袜人妻中文字幕| 看黄色毛片网站| 涩涩av久久男人的天堂| 久久精品国产亚洲av高清一级| 亚洲男人的天堂狠狠| 欧美日本视频| 久热这里只有精品99| 这个男人来自地球电影免费观看| 真人做人爱边吃奶动态| 一区福利在线观看| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| 国产精品,欧美在线| 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 国产av又大| 男男h啪啪无遮挡| 涩涩av久久男人的天堂| 久久久久久久久久久久大奶| 真人做人爱边吃奶动态| 一区福利在线观看| 欧美人与性动交α欧美精品济南到| 国产亚洲精品av在线| 久久狼人影院| 99香蕉大伊视频| 国产精品亚洲美女久久久| av欧美777| 黄色a级毛片大全视频| 韩国精品一区二区三区| 久久精品91无色码中文字幕| 午夜日韩欧美国产| 1024香蕉在线观看| 精品一区二区三区四区五区乱码| 少妇熟女aⅴ在线视频| 国产一卡二卡三卡精品| bbb黄色大片| 757午夜福利合集在线观看| 黄色女人牲交| 非洲黑人性xxxx精品又粗又长| 超碰成人久久| 一二三四社区在线视频社区8| 久久这里只有精品19| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 男人操女人黄网站| 色婷婷久久久亚洲欧美| 女性生殖器流出的白浆| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 国产片内射在线| 久久久久久久久久久久大奶| 日韩 欧美 亚洲 中文字幕| 日韩大码丰满熟妇| 欧美日本视频| 91精品三级在线观看| e午夜精品久久久久久久| 九色亚洲精品在线播放| 久久中文字幕人妻熟女| 精品第一国产精品| 欧美在线黄色| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 免费在线观看完整版高清| 极品人妻少妇av视频| 黄色成人免费大全| 丝袜美腿诱惑在线| 黄片播放在线免费| 国产精品久久久av美女十八| 久久香蕉激情| 久久国产精品影院| 免费女性裸体啪啪无遮挡网站| av天堂久久9| 亚洲精品中文字幕在线视频| 黄片小视频在线播放| 午夜老司机福利片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品粉嫩美女一区| av天堂久久9| 精品人妻1区二区| 国产精品久久视频播放| av电影中文网址| 亚洲av美国av| 熟女少妇亚洲综合色aaa.| 欧美日韩黄片免| 日本免费a在线| 日韩欧美国产在线观看| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 国产精品1区2区在线观看.| 国产午夜精品久久久久久| 国产亚洲欧美精品永久| 成人国产一区最新在线观看| 久久性视频一级片| 97人妻天天添夜夜摸| 国产亚洲欧美在线一区二区| 久久人妻av系列| 91成人精品电影| 久久久久九九精品影院| 女人爽到高潮嗷嗷叫在线视频| 亚洲,欧美精品.| 亚洲成国产人片在线观看| 欧美日韩黄片免| 最近最新中文字幕大全免费视频| 亚洲午夜理论影院| 久久久水蜜桃国产精品网| 欧美成人免费av一区二区三区| 身体一侧抽搐| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 视频区欧美日本亚洲| 亚洲第一电影网av| 亚洲国产精品999在线| 香蕉国产在线看| 日韩精品中文字幕看吧| 嫩草影院精品99| 最近最新中文字幕大全电影3 | 国产精品久久久人人做人人爽| 黄色片一级片一级黄色片| 精品少妇一区二区三区视频日本电影| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 88av欧美| svipshipincom国产片| 亚洲少妇的诱惑av| 国产极品粉嫩免费观看在线| 身体一侧抽搐| 一区二区三区精品91| 日韩国内少妇激情av| 国产欧美日韩精品亚洲av| 校园春色视频在线观看| 神马国产精品三级电影在线观看 | 色播亚洲综合网| 亚洲av片天天在线观看| 亚洲久久久国产精品| 亚洲国产中文字幕在线视频| 久久久国产欧美日韩av| 国内精品久久久久久久电影| 日韩免费av在线播放| 亚洲五月天丁香| 午夜福利影视在线免费观看| 精品久久久久久久毛片微露脸| 九色国产91popny在线| 国产一区在线观看成人免费| 欧美中文日本在线观看视频| 精品人妻1区二区| 后天国语完整版免费观看| 免费无遮挡裸体视频| 亚洲人成77777在线视频| 亚洲情色 制服丝袜| 国产又爽黄色视频| 欧美黄色淫秽网站| 亚洲 国产 在线| 亚洲五月婷婷丁香| 午夜两性在线视频| 69av精品久久久久久| 两人在一起打扑克的视频| 久久香蕉精品热| 99国产精品一区二区蜜桃av| 老司机深夜福利视频在线观看| 天堂√8在线中文| 国产高清videossex| 亚洲第一av免费看| www国产在线视频色| 自线自在国产av| 国产亚洲av嫩草精品影院| 97人妻天天添夜夜摸| 免费观看人在逋| 91精品国产国语对白视频| 国产精品,欧美在线| 亚洲avbb在线观看| 精品人妻1区二区| 久9热在线精品视频| 又紧又爽又黄一区二区| 国产精品久久视频播放| 黄色成人免费大全| 成年女人毛片免费观看观看9| 18禁观看日本| 久99久视频精品免费| 精品国产亚洲在线| a级毛片在线看网站| 一级毛片精品| 国产精品99久久99久久久不卡| 人妻丰满熟妇av一区二区三区| 手机成人av网站| 宅男免费午夜| 麻豆国产av国片精品| 国产成人精品久久二区二区免费| 国产三级在线视频| 中亚洲国语对白在线视频| 国产激情欧美一区二区| 久久久国产欧美日韩av| 1024视频免费在线观看| 日日摸夜夜添夜夜添小说| 啦啦啦免费观看视频1| 男人的好看免费观看在线视频 | 极品教师在线免费播放| 色尼玛亚洲综合影院| 一进一出抽搐gif免费好疼| 精品国产乱码久久久久久男人| 丰满人妻熟妇乱又伦精品不卡| 午夜精品国产一区二区电影| 51午夜福利影视在线观看| 国产人伦9x9x在线观看| 精品人妻在线不人妻| 精品日产1卡2卡| 日韩成人在线观看一区二区三区| www国产在线视频色| 涩涩av久久男人的天堂| 国产精品亚洲av一区麻豆| 此物有八面人人有两片| 日本撒尿小便嘘嘘汇集6| 老熟妇仑乱视频hdxx| 香蕉国产在线看| 一区二区三区激情视频| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 大型av网站在线播放| 欧美性长视频在线观看| 最近最新中文字幕大全免费视频| 国产精品国产高清国产av| 国产亚洲精品久久久久5区| 9191精品国产免费久久| 90打野战视频偷拍视频| 国产熟女xx| 亚洲人成77777在线视频| 一区二区三区国产精品乱码| 精品久久久久久久毛片微露脸| 国产在线精品亚洲第一网站| 丁香欧美五月| 十八禁人妻一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩精品亚洲av| √禁漫天堂资源中文www| 一级a爱片免费观看的视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 精品卡一卡二卡四卡免费| 久久精品国产清高在天天线| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 91国产中文字幕| 少妇 在线观看| 日日干狠狠操夜夜爽| 看免费av毛片| 在线观看一区二区三区| 免费在线观看亚洲国产| 老司机靠b影院| 久久亚洲真实| 精品高清国产在线一区| 国产高清激情床上av| 18禁黄网站禁片午夜丰满| 国语自产精品视频在线第100页| 天堂√8在线中文| 老司机福利观看| 日本免费一区二区三区高清不卡 | 色婷婷久久久亚洲欧美| 国产黄a三级三级三级人| 午夜久久久在线观看| 亚洲专区中文字幕在线| 99国产精品一区二区蜜桃av| 真人做人爱边吃奶动态| 丰满的人妻完整版| 免费不卡黄色视频| 不卡一级毛片| 国产精品久久久久久人妻精品电影| 久久国产亚洲av麻豆专区| 午夜免费观看网址| 首页视频小说图片口味搜索| 男人舔女人的私密视频| 桃红色精品国产亚洲av| 色婷婷久久久亚洲欧美| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 一区二区日韩欧美中文字幕| 两人在一起打扑克的视频| 午夜影院日韩av| 成年女人毛片免费观看观看9| 精品国产一区二区三区四区第35| 在线播放国产精品三级| 亚洲成av片中文字幕在线观看| 国内精品久久久久精免费| 变态另类丝袜制服| 国产熟女xx| 日本a在线网址| 亚洲午夜精品一区,二区,三区| 日本 欧美在线| 精品不卡国产一区二区三区| 如日韩欧美国产精品一区二区三区| 亚洲一区二区三区色噜噜| 日韩有码中文字幕| 国产亚洲精品第一综合不卡| 国产亚洲欧美精品永久| 丰满的人妻完整版| 国产99久久九九免费精品| 国产一区二区激情短视频| 99久久99久久久精品蜜桃| 国产亚洲精品久久久久5区| 色老头精品视频在线观看| 久久精品91无色码中文字幕| 国产成人欧美| 日韩一卡2卡3卡4卡2021年| 午夜影院日韩av| 成人三级做爰电影| 黄片小视频在线播放| 欧美最黄视频在线播放免费| 国产精品av久久久久免费| 99国产精品一区二区三区| 啪啪无遮挡十八禁网站| 亚洲中文字幕一区二区三区有码在线看 | 久久精品亚洲精品国产色婷小说| 51午夜福利影视在线观看| 精品福利观看| 好男人在线观看高清免费视频 | 久久久国产精品麻豆| 在线观看舔阴道视频| 两个人视频免费观看高清| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利视频1000在线观看 | 色综合欧美亚洲国产小说| 色综合亚洲欧美另类图片| 两个人看的免费小视频| 国产精品爽爽va在线观看网站 | 国产一级毛片七仙女欲春2 | 国产伦一二天堂av在线观看| 国产xxxxx性猛交| 老司机午夜福利在线观看视频| 久久精品国产亚洲av香蕉五月| 午夜福利一区二区在线看| 色尼玛亚洲综合影院| 日本精品一区二区三区蜜桃| 亚洲av成人一区二区三| 色播亚洲综合网| 精品国产亚洲在线| 极品人妻少妇av视频| 中文字幕人成人乱码亚洲影| 多毛熟女@视频| 亚洲精品久久成人aⅴ小说| 久久久久国产一级毛片高清牌| 亚洲精品一区av在线观看| 日韩有码中文字幕| 国产一区二区三区视频了| av超薄肉色丝袜交足视频| 亚洲一区二区三区不卡视频| 亚洲av日韩精品久久久久久密| 久久人妻av系列| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 亚洲少妇的诱惑av| 精品国产一区二区久久| 国产欧美日韩一区二区三区在线| 制服诱惑二区| 在线观看www视频免费| 免费看十八禁软件| 国产精品一区二区精品视频观看| 日韩av在线大香蕉| 免费在线观看视频国产中文字幕亚洲| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 久久久久久久午夜电影| 窝窝影院91人妻| 亚洲精品美女久久av网站| 精品国产乱子伦一区二区三区| 国产精品二区激情视频| 国产精品一区二区免费欧美| 国产精品免费一区二区三区在线| 久久午夜综合久久蜜桃| 久久精品国产清高在天天线| svipshipincom国产片| 久热这里只有精品99| 高清在线国产一区| 国产av一区二区精品久久| av天堂久久9| 亚洲无线在线观看| 成年版毛片免费区| 欧美日韩一级在线毛片| 中文字幕人成人乱码亚洲影| 人妻丰满熟妇av一区二区三区| 叶爱在线成人免费视频播放| 欧美成人午夜精品| 性少妇av在线| 欧美午夜高清在线| 男男h啪啪无遮挡| 免费人成视频x8x8入口观看| 老汉色∧v一级毛片| 亚洲成国产人片在线观看| 亚洲成人国产一区在线观看| 亚洲国产日韩欧美精品在线观看 | 国内久久婷婷六月综合欲色啪| 黄片小视频在线播放| 国产99久久九九免费精品| 成人免费观看视频高清| 久久精品亚洲精品国产色婷小说| 亚洲五月婷婷丁香| 午夜福利18| 伦理电影免费视频| 亚洲无线在线观看| 欧美精品亚洲一区二区| 无人区码免费观看不卡| 国产av一区二区精品久久| 久久国产亚洲av麻豆专区| 免费人成视频x8x8入口观看| 精品电影一区二区在线| 熟女少妇亚洲综合色aaa.| 久久久水蜜桃国产精品网| 少妇熟女aⅴ在线视频| 后天国语完整版免费观看| 国产成人系列免费观看| 欧美久久黑人一区二区| 久久中文字幕人妻熟女| 一区二区三区激情视频| 欧美中文日本在线观看视频| 亚洲久久久国产精品| 欧美性长视频在线观看| 久久久国产成人免费| 精品国产一区二区三区四区第35| 国产激情欧美一区二区| 老鸭窝网址在线观看| 日韩免费av在线播放| 国产成人啪精品午夜网站| 美女高潮到喷水免费观看| 亚洲av五月六月丁香网| 人人妻人人澡人人看| 国产免费av片在线观看野外av| 亚洲av片天天在线观看| 51午夜福利影视在线观看| 免费无遮挡裸体视频| 一二三四社区在线视频社区8| 无人区码免费观看不卡| 亚洲美女黄片视频| 午夜久久久在线观看| 亚洲,欧美精品.| 亚洲 欧美一区二区三区| 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 久久久精品欧美日韩精品| 久久久久国产精品人妻aⅴ院| 男人的好看免费观看在线视频 | 欧美成人免费av一区二区三区| 久久精品aⅴ一区二区三区四区| 久久久久国产一级毛片高清牌| 色在线成人网| 女人精品久久久久毛片| 人人妻,人人澡人人爽秒播| 亚洲成av人片免费观看| 岛国在线观看网站| 国内毛片毛片毛片毛片毛片| 夜夜躁狠狠躁天天躁| 午夜影院日韩av| 婷婷丁香在线五月| 国产精品免费一区二区三区在线| 男男h啪啪无遮挡| 国产精品久久久久久人妻精品电影| 国产一区二区三区视频了| 老司机午夜十八禁免费视频| 一区福利在线观看| 国产99久久九九免费精品| 欧美日韩乱码在线| 侵犯人妻中文字幕一二三四区| 一级作爱视频免费观看| 男女午夜视频在线观看| 精品久久蜜臀av无| 日韩欧美一区视频在线观看| 亚洲情色 制服丝袜| 视频区欧美日本亚洲| 老司机深夜福利视频在线观看| 日本 欧美在线| 亚洲一区二区三区色噜噜| 日本vs欧美在线观看视频| 男女下面进入的视频免费午夜 | 色尼玛亚洲综合影院| 天堂√8在线中文| 一本大道久久a久久精品| 久久天躁狠狠躁夜夜2o2o| 国产高清视频在线播放一区| 狂野欧美激情性xxxx| 日韩欧美三级三区| 欧美乱码精品一区二区三区| 午夜福利在线观看吧| 18禁观看日本| 亚洲avbb在线观看| 亚洲av成人av| 国产熟女xx| www.www免费av| 此物有八面人人有两片| 国产激情久久老熟女| 变态另类丝袜制服| 亚洲国产精品合色在线| 禁无遮挡网站| 无遮挡黄片免费观看| 久久午夜综合久久蜜桃| 一个人免费在线观看的高清视频| 午夜激情av网站| 可以在线观看的亚洲视频| 男人舔女人下体高潮全视频|