• <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)化
    延安至延川高速公路即將建成
    国产综合精华液| 亚洲成人一二三区av| 久久久久精品性色| 欧美 日韩 精品 国产| 亚洲天堂国产精品一区在线| 看黄色毛片网站| 欧美日本视频| 国产高潮美女av| 免费看a级黄色片| 亚洲精品乱码久久久久久按摩| 国产极品天堂在线| 99热国产这里只有精品6| 一本色道久久久久久精品综合| 听说在线观看完整版免费高清| 久久韩国三级中文字幕| 国产中年淑女户外野战色| 成年av动漫网址| 日韩电影二区| 高清午夜精品一区二区三区| 精华霜和精华液先用哪个| 国产成年人精品一区二区| 亚洲综合色惰| 97精品久久久久久久久久精品| 丝袜美腿在线中文| 亚洲伊人久久精品综合| 日韩欧美 国产精品| 伦精品一区二区三区| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在| 成人毛片a级毛片在线播放| 亚洲自拍偷在线| 成人国产av品久久久| 少妇人妻 视频| 97超视频在线观看视频| 国产黄片视频在线免费观看| 97精品久久久久久久久久精品| 91精品国产九色| 亚洲国产色片| 热99国产精品久久久久久7| 亚洲欧美中文字幕日韩二区| 久久久久久久国产电影| 精品久久久久久久人妻蜜臀av| 国产成人精品福利久久| 男女啪啪激烈高潮av片| 亚洲精品成人久久久久久| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 日本黄色片子视频| 美女被艹到高潮喷水动态| 人妻制服诱惑在线中文字幕| 九九在线视频观看精品| 黄片wwwwww| 大片电影免费在线观看免费| 亚洲精品影视一区二区三区av| 亚洲欧美中文字幕日韩二区| 国产在线一区二区三区精| 小蜜桃在线观看免费完整版高清| 国产精品一区二区性色av| 亚洲精品成人久久久久久| 亚洲国产欧美在线一区| 国语对白做爰xxxⅹ性视频网站| 少妇熟女欧美另类| 亚洲三级黄色毛片| 在线天堂最新版资源| 国产精品成人在线| 五月玫瑰六月丁香| 女人十人毛片免费观看3o分钟| 亚洲最大成人手机在线| 特大巨黑吊av在线直播| 91狼人影院| 国产女主播在线喷水免费视频网站| 国内揄拍国产精品人妻在线| 欧美日韩视频高清一区二区三区二| 亚洲最大成人手机在线| 涩涩av久久男人的天堂| 久久久久久久久久久免费av| 中文字幕av成人在线电影| 日韩伦理黄色片| 国产国拍精品亚洲av在线观看| 激情五月婷婷亚洲| 3wmmmm亚洲av在线观看| 51国产日韩欧美| 国产亚洲av嫩草精品影院| 插阴视频在线观看视频| 日本一二三区视频观看| 欧美zozozo另类| 搡老乐熟女国产| 日韩大片免费观看网站| 国产综合精华液| 国产精品一区www在线观看| 免费黄网站久久成人精品| 国产精品国产av在线观看| 国产老妇女一区| 日本与韩国留学比较| 国产欧美另类精品又又久久亚洲欧美| 美女内射精品一级片tv| 在线亚洲精品国产二区图片欧美 | 亚洲欧美中文字幕日韩二区| 成人亚洲精品av一区二区| av专区在线播放| 99视频精品全部免费 在线| 亚洲熟女精品中文字幕| a级毛片免费高清观看在线播放| 大片免费播放器 马上看| 久久人人爽人人片av| 国产亚洲av片在线观看秒播厂| 欧美高清成人免费视频www| 制服丝袜香蕉在线| 国产成人91sexporn| 在线观看美女被高潮喷水网站| 在线免费十八禁| 又爽又黄a免费视频| 午夜视频国产福利| 高清午夜精品一区二区三区| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 国产v大片淫在线免费观看| 亚洲天堂国产精品一区在线| 国产精品久久久久久久电影| 亚洲av男天堂| 最近手机中文字幕大全| 国产色爽女视频免费观看| 人妻制服诱惑在线中文字幕| 国产 一区精品| 夫妻午夜视频| 亚洲成人久久爱视频| 观看免费一级毛片| 亚洲电影在线观看av| 又大又黄又爽视频免费| 亚洲国产精品成人综合色| 97精品久久久久久久久久精品| 一级爰片在线观看| 成人亚洲精品一区在线观看 | 涩涩av久久男人的天堂| 成人高潮视频无遮挡免费网站| 久久久亚洲精品成人影院| 一级毛片 在线播放| 青春草亚洲视频在线观看| 各种免费的搞黄视频| 在线观看一区二区三区激情| 18禁动态无遮挡网站| 国产大屁股一区二区在线视频| 国产91av在线免费观看| 99热这里只有精品一区| 嘟嘟电影网在线观看| 麻豆精品久久久久久蜜桃| 国产中年淑女户外野战色| 久久久久精品性色| 国产黄色视频一区二区在线观看| 丰满少妇做爰视频| 欧美成人一区二区免费高清观看| 欧美日韩国产mv在线观看视频 | 深夜a级毛片| 自拍偷自拍亚洲精品老妇| 街头女战士在线观看网站| 亚洲精品久久久久久婷婷小说| 97人妻精品一区二区三区麻豆| 激情 狠狠 欧美| 亚洲最大成人手机在线| 伊人久久国产一区二区| 亚洲色图av天堂| 久久国产乱子免费精品| 久久综合国产亚洲精品| 国产精品国产三级国产专区5o| 国产日韩欧美在线精品| 国产一区二区三区综合在线观看 | 狂野欧美激情性bbbbbb| 色5月婷婷丁香| 蜜臀久久99精品久久宅男| 日本爱情动作片www.在线观看| 欧美丝袜亚洲另类| 尾随美女入室| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩无卡精品| 黄片无遮挡物在线观看| 中文精品一卡2卡3卡4更新| 精品久久久久久久末码| 亚洲人成网站高清观看| 在线天堂最新版资源| 综合色av麻豆| 亚洲精品成人av观看孕妇| 成人欧美大片| 亚洲精品,欧美精品| 在线观看国产h片| 国产中年淑女户外野战色| 欧美高清成人免费视频www| 婷婷色麻豆天堂久久| 成人特级av手机在线观看| 日韩 亚洲 欧美在线| 免费看不卡的av| 男人舔奶头视频| 国产视频首页在线观看| 久久久久久久大尺度免费视频| 国产精品99久久久久久久久| 日本wwww免费看| 国产亚洲91精品色在线| 欧美日本视频| 中文字幕制服av| 色视频www国产| 性色av一级| 亚洲久久久久久中文字幕| 国产成人a区在线观看| 特级一级黄色大片| 国产老妇女一区| 久久久亚洲精品成人影院| 视频区图区小说| 久久久久久久久久成人| 高清午夜精品一区二区三区| av网站免费在线观看视频| 男女边摸边吃奶| 一区二区av电影网| 精品久久久久久电影网| 亚洲av.av天堂| 国产黄a三级三级三级人| 男人爽女人下面视频在线观看| 午夜免费观看性视频| 亚洲一区二区三区欧美精品 | 超碰97精品在线观看| a级一级毛片免费在线观看| 男的添女的下面高潮视频| 亚洲精品国产成人久久av| 人妻系列 视频| 麻豆成人av视频| 婷婷色综合大香蕉| 插阴视频在线观看视频| 久久久欧美国产精品| 免费观看a级毛片全部| av在线播放精品| 内地一区二区视频在线| 岛国毛片在线播放| 国产成人福利小说| 国产爽快片一区二区三区| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 日本爱情动作片www.在线观看| 久久久久久久大尺度免费视频| 白带黄色成豆腐渣| 国产一区亚洲一区在线观看| 老司机影院成人| 国产探花在线观看一区二区| 日韩欧美精品v在线| 久久韩国三级中文字幕| 国产探花在线观看一区二区| 亚洲成人精品中文字幕电影| 水蜜桃什么品种好| 日本黄大片高清| 国产高清有码在线观看视频| 在线播放无遮挡| 日韩三级伦理在线观看| 久久精品国产鲁丝片午夜精品| 自拍偷自拍亚洲精品老妇| 热99国产精品久久久久久7| 岛国毛片在线播放| 天美传媒精品一区二区| 欧美高清成人免费视频www| 99久久中文字幕三级久久日本| 热re99久久精品国产66热6| 国产成人精品久久久久久| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 欧美精品国产亚洲| av线在线观看网站| 3wmmmm亚洲av在线观看| 一级av片app| 久久久久久九九精品二区国产| 国产免费一区二区三区四区乱码| 国产欧美亚洲国产| 久久久久久久精品精品| 国产毛片在线视频| 日韩强制内射视频| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱久久久久久| 性色av一级| 欧美激情久久久久久爽电影| 美女国产视频在线观看| 日韩欧美精品免费久久| 午夜视频国产福利| 少妇被粗大猛烈的视频| 国语对白做爰xxxⅹ性视频网站| 成年女人在线观看亚洲视频 | 国产在视频线精品| 人人妻人人看人人澡| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 精品酒店卫生间| 大码成人一级视频| 如何舔出高潮| 久久久亚洲精品成人影院| 69人妻影院| 精品久久久噜噜| 国产精品福利在线免费观看| a级毛色黄片| av福利片在线观看| 99九九线精品视频在线观看视频| 国产亚洲一区二区精品| 亚洲av日韩在线播放| 天美传媒精品一区二区| 国产一区有黄有色的免费视频| 午夜福利在线在线| 国语对白做爰xxxⅹ性视频网站| 午夜视频国产福利| 少妇人妻 视频| 热re99久久精品国产66热6| 高清午夜精品一区二区三区| 在线观看国产h片| 一级毛片aaaaaa免费看小| 亚洲欧美日韩无卡精品| 麻豆成人午夜福利视频| 日本黄大片高清| 91在线精品国自产拍蜜月| 一级爰片在线观看| av国产免费在线观看| 男人和女人高潮做爰伦理| 亚洲国产精品国产精品| 大话2 男鬼变身卡| 亚洲色图综合在线观看| 欧美潮喷喷水| 国产精品久久久久久精品电影| 国产精品熟女久久久久浪| 日韩在线高清观看一区二区三区| 97超视频在线观看视频| 美女视频免费永久观看网站| 成人毛片60女人毛片免费| 欧美性感艳星| 在线观看一区二区三区激情| 丰满人妻一区二区三区视频av| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 视频中文字幕在线观看| 亚洲久久久久久中文字幕| 国产高潮美女av| 欧美zozozo另类| 日本三级黄在线观看| 91久久精品电影网| 性插视频无遮挡在线免费观看| av专区在线播放| 日韩不卡一区二区三区视频在线| 亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 视频中文字幕在线观看| 成人综合一区亚洲| 久久99热这里只有精品18| 麻豆成人午夜福利视频| 亚洲经典国产精华液单| av在线天堂中文字幕| 亚洲不卡免费看| 日本-黄色视频高清免费观看| 丰满人妻一区二区三区视频av| 六月丁香七月| 成年免费大片在线观看| 中国美白少妇内射xxxbb| 日韩av在线免费看完整版不卡| 亚州av有码| 99热这里只有是精品在线观看| 日本免费在线观看一区| 久久综合国产亚洲精品| 欧美高清性xxxxhd video| 日本与韩国留学比较| 国产av国产精品国产| 国产欧美日韩精品一区二区| 视频区图区小说| 秋霞在线观看毛片| 国产精品久久久久久久电影| 久久97久久精品| a级一级毛片免费在线观看| 国产伦精品一区二区三区视频9| 久久久久性生活片| 在线看a的网站| 偷拍熟女少妇极品色| 国产一区有黄有色的免费视频| 欧美日韩在线观看h| 亚洲av免费高清在线观看| 日本熟妇午夜| 久久久精品欧美日韩精品| 久久精品国产鲁丝片午夜精品| 自拍偷自拍亚洲精品老妇| 欧美xxxx黑人xx丫x性爽| kizo精华| 精品久久久久久电影网| 欧美成人精品欧美一级黄| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 午夜精品一区二区三区免费看| 2022亚洲国产成人精品| 国产亚洲av嫩草精品影院| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 最近最新中文字幕大全电影3| 国产精品偷伦视频观看了| 免费观看的影片在线观看| 可以在线观看毛片的网站| 国内揄拍国产精品人妻在线| 人妻系列 视频| av在线蜜桃| 久久6这里有精品| 免费av观看视频| 亚洲自偷自拍三级| 亚洲精品乱码久久久v下载方式| 日韩av不卡免费在线播放| 18禁在线无遮挡免费观看视频| 看黄色毛片网站| 又粗又硬又长又爽又黄的视频| 精品久久久噜噜| 精品久久久久久久久av| 熟女av电影| 亚洲精品久久午夜乱码| .国产精品久久| 日韩免费高清中文字幕av| 免费观看在线日韩| 简卡轻食公司| av播播在线观看一区| 国产精品久久久久久久久免| 亚洲怡红院男人天堂| 久久精品国产鲁丝片午夜精品| 免费电影在线观看免费观看| 联通29元200g的流量卡| 欧美性猛交╳xxx乱大交人| 国产成人freesex在线| av在线播放精品| 波野结衣二区三区在线| 日韩亚洲欧美综合| 免费高清在线观看视频在线观看| 美女视频免费永久观看网站| 亚洲伊人久久精品综合| 别揉我奶头 嗯啊视频| 国产午夜精品久久久久久一区二区三区| 老司机影院成人| 亚洲欧美日韩东京热| 欧美少妇被猛烈插入视频| 国产大屁股一区二区在线视频| 黄片wwwwww| 极品教师在线视频| 国产精品秋霞免费鲁丝片| 丝袜喷水一区| 成人黄色视频免费在线看| 麻豆成人午夜福利视频| 大陆偷拍与自拍| 一个人观看的视频www高清免费观看| 国产色爽女视频免费观看| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 中文字幕av成人在线电影| 丝袜美腿在线中文| 2021少妇久久久久久久久久久| 国产有黄有色有爽视频| 国产成人aa在线观看| 日产精品乱码卡一卡2卡三| 成人二区视频| 久久久久久久久大av| 亚洲三级黄色毛片| 偷拍熟女少妇极品色| 99久久九九国产精品国产免费| 日韩欧美精品免费久久| 日韩人妻高清精品专区| 女人被狂操c到高潮| 亚洲欧洲国产日韩| 日韩视频在线欧美| 舔av片在线| 看免费成人av毛片| 一区二区三区乱码不卡18| 日韩欧美 国产精品| 国产亚洲一区二区精品| 日韩制服骚丝袜av| 一级毛片aaaaaa免费看小| 18禁在线播放成人免费| 三级男女做爰猛烈吃奶摸视频| 在线 av 中文字幕| 亚洲成人精品中文字幕电影| 新久久久久国产一级毛片| 2018国产大陆天天弄谢| 内射极品少妇av片p| 六月丁香七月| 国产爱豆传媒在线观看| 超碰97精品在线观看| 精品一区二区免费观看| 精品视频人人做人人爽| 综合色丁香网| 九九在线视频观看精品| 最近的中文字幕免费完整| av免费观看日本| 国产 一区 欧美 日韩| 久久久精品94久久精品| 成年版毛片免费区| 五月天丁香电影| av播播在线观看一区| 久久精品国产a三级三级三级| 大码成人一级视频| 久久久久久久久久成人| 婷婷色综合www| 神马国产精品三级电影在线观看| 高清毛片免费看| 777米奇影视久久| 91狼人影院| 一区二区三区免费毛片| 91午夜精品亚洲一区二区三区| 久久久久久久大尺度免费视频| 一本色道久久久久久精品综合| 18+在线观看网站| 91精品伊人久久大香线蕉| 欧美xxⅹ黑人| 美女高潮的动态| 国产精品久久久久久精品电影| 成人午夜精彩视频在线观看| 少妇猛男粗大的猛烈进出视频 | 中文乱码字字幕精品一区二区三区| 亚洲欧美精品专区久久| 丝袜喷水一区| 久久综合国产亚洲精品| 国产熟女欧美一区二区| 青春草国产在线视频| 精品国产三级普通话版| 一级毛片aaaaaa免费看小| 91在线精品国自产拍蜜月| 成年女人看的毛片在线观看| 免费看av在线观看网站| 亚洲av成人精品一二三区| 美女被艹到高潮喷水动态| av国产精品久久久久影院| 国产免费福利视频在线观看| 国产精品伦人一区二区| 五月玫瑰六月丁香| 不卡视频在线观看欧美| 老司机影院成人| 国产精品偷伦视频观看了| 久久精品国产a三级三级三级| 少妇的逼好多水| 国产成人一区二区在线| 亚洲最大成人手机在线| 成人黄色视频免费在线看| 一级爰片在线观看| 校园人妻丝袜中文字幕| 天堂中文最新版在线下载 | 在线观看av片永久免费下载| 熟妇人妻不卡中文字幕| 成人特级av手机在线观看| 亚洲成色77777| 51国产日韩欧美| 99久久精品热视频| 亚洲精品一二三| 97热精品久久久久久| 亚洲精品日韩在线中文字幕| av黄色大香蕉| 国产精品女同一区二区软件| 插逼视频在线观看| 亚洲成人一二三区av| 五月玫瑰六月丁香| 人妻制服诱惑在线中文字幕| 国产精品一区www在线观看| 视频区图区小说| 国产午夜精品一二区理论片| 欧美成人a在线观看| 亚洲精品亚洲一区二区| 精品人妻一区二区三区麻豆| 久久久精品欧美日韩精品| 内地一区二区视频在线| 久久久精品94久久精品| 大香蕉久久网| 午夜激情福利司机影院| 我的女老师完整版在线观看| 欧美激情国产日韩精品一区| 在现免费观看毛片| 久久精品久久久久久噜噜老黄| 成人欧美大片| 欧美xxxx黑人xx丫x性爽| 只有这里有精品99| 一级爰片在线观看| 亚洲,一卡二卡三卡| 亚洲av欧美aⅴ国产| 蜜桃亚洲精品一区二区三区| 国产精品秋霞免费鲁丝片| 真实男女啪啪啪动态图| av免费在线看不卡| 老女人水多毛片| 国产成年人精品一区二区| 小蜜桃在线观看免费完整版高清| 亚洲色图综合在线观看| 精品熟女少妇av免费看| 亚洲精品成人av观看孕妇| 亚洲成人中文字幕在线播放| 亚洲美女搞黄在线观看| 在线观看三级黄色| 欧美一级a爱片免费观看看| 久久久精品94久久精品| 一个人看视频在线观看www免费| 毛片一级片免费看久久久久| 老司机影院毛片| 97热精品久久久久久| 三级国产精品欧美在线观看| 又爽又黄无遮挡网站| 国产成人精品婷婷| 人人妻人人澡人人爽人人夜夜| 99热这里只有是精品在线观看| xxx大片免费视频| 三级经典国产精品| 美女内射精品一级片tv| 人体艺术视频欧美日本| 一本久久精品| 久久久久久九九精品二区国产| 欧美3d第一页| 在线观看三级黄色| 99久久九九国产精品国产免费| 老司机影院成人| 菩萨蛮人人尽说江南好唐韦庄| 日韩一区二区视频免费看| 亚洲,一卡二卡三卡| 婷婷色麻豆天堂久久| 国产大屁股一区二区在线视频| 国产高清有码在线观看视频| 亚洲av中文av极速乱| 在线观看一区二区三区激情| 午夜福利视频精品| 亚洲精品,欧美精品| 成年女人在线观看亚洲视频 | 日本黄大片高清|