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

    黃土高原坡度因子的多重分形譜特征

    2021-08-09 12:54:42夏積德叢佩娟李俊鋒
    水土保持通報 2021年3期
    關鍵詞:黃土塬分形坡度

    夏積德, 叢佩娟, 李俊鋒

    (1.楊凌職業(yè)技術學院, 陜西 楊凌 712100; 2.水利部 水土保持監(jiān)測中心, 北京 100053; 3.陜西測繪地理信息局, 陜西 西安 710054)

    地形因子是水土流失、水文過程模擬的參數(shù),科學量化地形因子對于水土流失這一地理界面過程的研究具有極為重要的意義。大量研究發(fā)現(xiàn),自然界中大部分區(qū)域的地表形態(tài)復雜多變,呈現(xiàn)出明顯的不均勻性和分形特征,即具有以非整數(shù)維形式充填空間的形態(tài)特征。若僅采用普通分形維數(shù)的方法并不能詳細且全面地反映不同區(qū)域條件、不同尺度下的地形因子特征。在多重分形理論被提出后,由于其采用廣義維數(shù)和多重分形譜來描述分形客體,考慮了對象屬性在幾何支集的空間奇異性分布這一特點[1],即能考慮到復雜分形體內(nèi)部不同區(qū)域呈現(xiàn)出不同分形特征,因而逐漸被用于包括地形地貌學在內(nèi)的各種地理學相關研究。沈中原等[2]借助多重分形理論構建了海拔信息空間分布概率計算模型;張建興等[3]探討了黃土高原若干流域的河網(wǎng)分形特征并基于此建立起產(chǎn)沙模數(shù)模型;王民等[4-5]運用多重分形理論及其實現(xiàn)方法分析了大理河流域地表形態(tài)的多重分形特征。這些研究多運用分形理論表征區(qū)域地貌形態(tài)的綜合特征[6],但對于黃土高原地區(qū)地形因子(如坡度)的多重分形研究還很少涉及。因此,本文運用多重分形理論來研究黃土高原地區(qū)地形因子的特點,旨在為黃土高原地區(qū)開展地形因子方面的研究提供參考。本研究對區(qū)域水文過程、水土流失預報模型的建立具有積極意義。

    1 研究區(qū)概況

    黃土高原地處中國中北部地區(qū),轄內(nèi)蒙古、山西、陜西、寧夏、甘肅、青海、河南7個省區(qū),東西跨越1 000 km多,南北跨越近750 km,總面積約為6.40×104km2,海拔高度介于1 000~3 000 m,是中國4大高原之一。黃土高原整體表現(xiàn)出西北高東南低的特點,寧夏西南方的六盤山和山西西部的呂梁山將整個黃土高原劃分為3個部分,在此劃分基礎上進一步按照地形地貌特點將黃土高原劃分為山地區(qū)、黃土丘陵區(qū)、黃土塬區(qū)以及河谷平原區(qū)4種類型。山地區(qū)主要包括馬銜山地區(qū)、拉脊山地區(qū)、六盤山地區(qū)、子午嶺地區(qū)、太行山地區(qū)和呂梁山地區(qū),海拔跨度大、地形地貌差異較大,由于長期受重力、流水以及風力等侵蝕因素的影響,山地區(qū)的地貌侵蝕現(xiàn)象嚴重,溝壑明顯,坡度較大;丘陵區(qū)作為黃土高原地區(qū)較為典型的地形地貌區(qū)域有著十分廣闊的分布區(qū)域,其占據(jù)了超過黃土高原地區(qū)一半的面積。丘陵區(qū)在黃土高原分布范圍遼闊,不同區(qū)域的丘陵區(qū)之間還存在不小的差異;黃土塬區(qū)中分布著大量塬面,塬面頂部十分平坦,地形復雜程度較低。與之形成強烈對比的是塬面周圍陡峭的溝壑,這些溝壑分布在一個又一個塬面之間,使得黃土塬區(qū)內(nèi)部的地形地貌形成了十分大的反差,地形復雜度對比強烈;河谷平原區(qū)以關中平原最具代表性,關中平原位于陜西省中部,呈東西帶狀分布,地形較為單一,多為平原,在長期的沖擊下形成了平坦開闊的地貌景觀。

    2 研究方法

    在黃土高原不同地貌類型區(qū)內(nèi)選擇合適的研究區(qū),獲取各子研究區(qū)的DEM(來源于地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/,空間分辨率為30 m×30 m,并在GIS中提取得到坡度因子的柵格數(shù)據(jù);隨后,采用盒維數(shù)的計算方法,計算得到坡度因子的多重分形譜并對其主要參數(shù)及圖像特征進行對比分析。

    2.1 多重分形譜的計算

    盒維數(shù)方法是一種基于統(tǒng)計學原理的計算多重分形的方法,具有簡便且易于理解的特點,經(jīng)常被用于地形地貌多重分形的計算[7]。盒維數(shù)是用一系列緊密相連的正方形盒子去覆蓋分形體,即將研究對象分為很多個(假設有i行j列總計N個)尺度大小為r的正方形網(wǎng)格單元[8],然后計算每個網(wǎng)格單元內(nèi)的概率測度。對于坡度因子來說,概率測度的大小可以用坡度概率函數(shù)表示,每個小網(wǎng)格單元的坡度概率函數(shù)的計算公式為:

    (1)

    式中:Pij(r) 表示位于第i行第j列的網(wǎng)格的坡度概率函數(shù);Air(r)表示位于第i行第j列的網(wǎng)格內(nèi)的坡度統(tǒng)計值;N表示網(wǎng)格單元的總個數(shù),其在數(shù)值上滿足:N=i×j。

    根據(jù)公式(1)并結合多重分形理論[9-10]可知,坡度概率函數(shù)可以粗略衡量一個網(wǎng)格單元區(qū)域的地形特點,數(shù)值越小表示該網(wǎng)格區(qū)域內(nèi)坡度越小,具有越平坦的地形特點,如塬面、平原等;反之,坡度概率函數(shù)的數(shù)值越大表示網(wǎng)格區(qū)域內(nèi)的坡度越大,地形起伏程度越大,如峁、坡面等。

    在盒維數(shù)方法中,檢驗一個研究對象是否具有多重分形特征,有一個配分函數(shù)的概念[11],配分函數(shù)的計算公式如下:

    (2)

    式中:Xq(r)表示分形體的配分函數(shù);Pij(r)表示位于第i行第j列的網(wǎng)格的坡度概率函數(shù);q表示階矩,可以取任意實數(shù),配分函數(shù)中q的取值直接影響著不同的概率測度在整個配分函數(shù)上所占比重的大小。

    對于給定的一對r和q,都可以通過公式(1)和公式(2)計算得到一個配分函數(shù)。在盒維數(shù)方法中,如果研究對象具有多重分形特征,那么配分函數(shù)Xq(r)和r之間應該滿足以下的關系[12]:

    Xq(r)∝rτ(q)

    (3)

    式中:τ(q)是關于q的一個表達式,在探討研究對象是否具有多重分形特征的過程中無關緊要。

    根據(jù)上述關系可知,對兩邊同時取對數(shù),lnXq(r)和lnr滿足線性關系。因此繪制以lnr為橫軸,lnXq(r)為縱軸的散點分布圖,對其擬合,如果lnXq(r)—lnr圖像為一系列斜率不等的直線,則表明研究對象具有多重分形特征[13],可以進一步計算其多重分形譜及其特征參數(shù)。

    本研究利用Chhabra等[14]提出的一種建立在盒維數(shù)方法基礎上的通過擬合參數(shù)計算多重分形譜的方法來計算多重分形譜。這種方法需要設立一個經(jīng)過歸一化處理的測度集μij(q,r),其計算公式如下:

    (4)

    在這個測度集的基礎上結合多重分形的概念,可以通過下面的公式來計算多重分形中的奇異指數(shù)α和分形維數(shù)f(α):

    (5)

    (6)

    奇異指數(shù)α是一個標度指數(shù),它是反映分形體內(nèi)各個小區(qū)間的奇異程度的一個量,α的數(shù)值與其所在的位置有關。將分形上具有相同α值的小區(qū)間數(shù)記為f(α),基于此可以得到α—f(α)的函數(shù)圖像,稱為奇異譜。橫坐標的范圍是多重分形譜奇異指數(shù)α的分布范圍,α的大小反映該部分的奇異程度。

    在計算時,常常按照公式(5)和公式(6)的數(shù)學幾何意義借助擬合的方法來得到α和f(α)的值。將公式(5)處理可得:

    (7)

    公式(7)根據(jù)一次函數(shù)相關概念及其圖像性質可知α為∑μij(q,r)·lnPij(r)—lnr圖像的斜率,類似地,f(α)為∑〔μij(q,r)·lnμij(q,r)〕—lnr圖像的斜率。

    首先確定一個q值,保持階矩q不變而改變尺度r,每選取一個r進行計算,按照公式(1),(2),(4)就可以得到一個對應的∑μij(q,r)·lnPij(r)—lnr和∑〔μij(q,r)·lnμij(q,r)〕,通過在二維坐標系中進行擬合便可得到每個q對應的奇異指數(shù)α和分形維數(shù)f(α)。改變階矩q的值,便能得到一系列的α和f(α),而多重分形譜就是以α為橫坐標,f(α)為縱坐標的圖像。

    2.2 特征參數(shù)的取值

    (1) 格網(wǎng)大小(尺度)r的選取。根據(jù)多重分形理論,格網(wǎng)的大小直接影響到格網(wǎng)的總數(shù)量和每個格網(wǎng)的概率測度,這兩個因素決定著多重分形計算過程中的效率和準確性,因此從簡便、可實施的角度出發(fā),需選取大小為30 m整數(shù)倍的格網(wǎng)來覆蓋研究區(qū)。結合所選4個研究區(qū)坡度數(shù)據(jù)的實際情況,規(guī)定一個柵格大小即30 m為一個單位,通過反復試驗最終選取格網(wǎng)尺度大小為20,30, 40,…,190,200個單位,亦即600,900,1 200,…,5 700,6 000 m,共19種。

    (2) 階距q取值范圍的確定。公式(3)配分函數(shù)中q的取值直接影響著不同的概率測度在整個配分函數(shù)上所占比重的大小[15],當q取正值且取值越大,配分函數(shù)越能反映具有較大坡度概率函數(shù)的那些網(wǎng)格單元的特征和性質;當q取負值且取值越小,則配分函數(shù)越能反映具有較小坡度概率函數(shù)的那些網(wǎng)格單元的特征和性質。因而q的取值對配分函數(shù)的計算結果影響很大,不同區(qū)域坡度的位置分布和大小范圍不同,各研究區(qū)q的取值范圍也不同,研究表明,在地形地貌多重分形的計算中,q值分布在0左右且通常介于±30之間[4]。

    本研究對每一個子研究區(qū)在-50~50之間等間距選取不同的q值,通過盒維數(shù)計算多重分形譜的方法驗證q是否在合適的區(qū)間。在初步確定q的合適區(qū)間后繼續(xù)加密q的取值,最終得到每個研究區(qū)階矩q的最佳取值范圍見表1。

    表1 研究區(qū)q最佳取值范圍

    3 結果與分析

    3.1 坡度概率分布特點

    為更詳細、系統(tǒng)地了解各研究區(qū)坡度大小的分布情況,對4個研究區(qū)的坡度進行統(tǒng)計(圖1)。可以看出,山地區(qū)、黃土丘陵區(qū)、黃土塬區(qū)、河谷平原區(qū)的坡度介于0°~10°之間的柵格占比分別為35.80%,35.50%,60.75%,>99.99%;介于10°~20°之間的柵格占比分別為24.23%,40.59%,16.94%,<0.01%;介于20°~30°之間的柵格占比分別為20.48%,20.92%,13.74%,0%;介于30°~40°之間的柵格占比分別為15.76%,2.89%,7.54%,0%;介于40°~50°之間的柵格占比分別為3.63%,0.10%,1.02%,0%;介于50°~60°之間的柵格占比分別為0.1%,0%,0.01%,0%。

    圖1 研究區(qū)不同地形坡度統(tǒng)計

    由統(tǒng)計結果得知:①山地區(qū)的坡度分布區(qū)間最大,大部分柵格處于0°~20°這一中小坡度區(qū)間,坡度大于40°的區(qū)間在4個研究區(qū)中占比最高;②黃土丘陵區(qū)坡度介于0°~10°區(qū)間內(nèi)的占比最少,坡度大部分位于10°~20°這一坡度范圍內(nèi)。同時,大于30°的坡度占比較小,因此在黃土丘陵研究區(qū)內(nèi)大坡度的地形較少,地形特點較為均一;③黃土塬區(qū)有大面積的平坦塬面,坡度大部分介于0°~10°之間,因此坡度值往往很小,由于所選研究區(qū)面積固定為30 km×30 km,在這樣一個較大的范圍內(nèi)不可避免地包括了塬面周圍的切溝、坡面等非塬面地形,因此也存在少量坡度大于20°的非塬面特征柵格;④河谷平原區(qū)的坡度統(tǒng)計圖特征十分明顯,幾乎所有柵格的坡度均介于0°~10°之間,在4個研究區(qū)中地形最為平坦。

    3.2 多重分形特征檢驗結果

    在計算坡度信息多重分形譜前需要先驗證研究對象是否具有多重分形特征。為此,在平面坐標系中,以lnr為橫坐標,lnXq(r)為縱坐標,繪制出每個研究區(qū)的lnXq(r)曲線,4個研究區(qū)坡度地形因子繪制的lnXq(r)—lnr曲線見圖2。

    從圖2中可以看出,每個研究區(qū)在各自合適的q范圍內(nèi),當尺度對數(shù)滿足3.0≤lnr≤5.3時,配分函數(shù)對數(shù)與尺度對數(shù)之間存在很好的線性關系,表明黃土高原區(qū)域上述4個研究區(qū)的坡度信息具有明顯的多重分形特征,可以進行多重分形譜的計算。

    圖2 研究區(qū)不同地形lnXq(r)-lnr特征

    3.3 坡度信息多重分形譜及其參數(shù)

    以坡度數(shù)據(jù)為基礎,利用多重分形理論原理和計算方法得到研究區(qū)坡度信息多重分形譜散點分布圖(圖3)及其參數(shù)結果(表2)。由圖3的多重分形譜可以看出,不同研究區(qū)的多重分形譜均為凸函數(shù)。在一個研究區(qū)內(nèi),多重分形譜函數(shù)值隨著坡度奇異指數(shù)的增加均呈現(xiàn)出先遞增后遞減的變化特征。而在各研究區(qū)之間,坡度多重分形譜的對稱性、圖像分布范圍均有不同,這是由于各研究區(qū)多重分形譜參數(shù)的差異導致的。通過對多重分形譜圖的二維坐標系進行擬合,得到各研究區(qū)詳細的多重分形譜參數(shù)信息見表2。

    圖3 研究區(qū)不同地形多重分形譜特征

    表2 研究區(qū)多重分形特征參數(shù)

    從多重分形譜參數(shù)統(tǒng)計表中可以看出,最小坡度奇異指數(shù)的最大值為1.94,出現(xiàn)在黃土丘陵區(qū),最小值為1.82,出現(xiàn)在黃土塬區(qū),說明在4個研究區(qū)中,坡度最大的地貌單元出現(xiàn)在黃土塬區(qū)(奇異指數(shù)最小值對應的是坡度最大地貌單元的性質);最小坡度奇異指數(shù)相應的多重分形譜值的最大值為1.78,出現(xiàn)在黃土塬區(qū),最小值為1.45,出現(xiàn)在黃土丘陵區(qū),這說明最小坡度奇異指數(shù)對應的地貌類型單元在黃土塬區(qū)所占比例最大,這與魯克新等[16]的研究結果相一致,在黃土丘陵區(qū)所占比例最小,即坡度最大的子區(qū)域在黃土塬區(qū)占比最高,而在黃土丘陵區(qū)占比最低;最大坡度奇異指數(shù)的最大值為2.90,出現(xiàn)在黃土塬區(qū),最小值為2.07,出現(xiàn)在黃土丘陵區(qū),說明地形特征最簡單的地貌單元出現(xiàn)在黃土塬區(qū);最大坡度奇異指數(shù)所對應的多重分形譜值最大值為1.63,出現(xiàn)在河谷平原區(qū),最小值為1.56,出現(xiàn)在黃土塬區(qū),說明最大坡度奇異指數(shù)對應的地貌類型單元在河谷平原區(qū)所占比例最大,在黃土塬區(qū)所占比例最小。黃土塬區(qū)對應的坡度奇異指數(shù)的極值最大,為1.08,山地區(qū)次之,為0.50,然后是河谷平原區(qū)的0.35,而黃土丘陵區(qū)為0.13。

    在多重分形譜左端點值與右端點值之差這個參數(shù)上,從絕對值的大小來看,黃土塬區(qū)最大,數(shù)值為0.22,說明該研究區(qū)最小坡度奇異指數(shù)對應的地貌單元所占比例與最大坡度奇異指數(shù)對應的地貌單元所占比例的差異性最大,黃土丘陵區(qū)次之,差異性最小的是河谷平原區(qū)。從正負的角度來看,黃土丘陵區(qū)和河谷平原區(qū)為負值,山地區(qū)、黃土塬區(qū)為正值,說明前兩個研究區(qū)的大坡度地貌單元所占比重較后兩個研究區(qū)要小。

    4 討論與結論

    目前,已有大量研究者采用簡單分形維數(shù)、多重分形譜和聯(lián)合多重分形等方法來表征復雜地貌的形態(tài)及其特征[17],但對于坡度因子的多重分形研究還很少涉及。本研究在DEM數(shù)據(jù)的基礎上,將黃土高原按照不同地形地貌劃分為山地區(qū)、黃土丘陵區(qū)、黃土塬區(qū)、河谷平原區(qū),并對坡度因子的多重分形譜圖像和多重分形特征參數(shù)進行了對比分析,總結出各自區(qū)域的多重分形譜特點,為后續(xù)開展黃土高原地形因子方面的研究打下了基礎,對區(qū)域水文過程、水土流失預報模型的建立具有積極意義。

    像元尺度的選取是地形因子分形維數(shù)計算的關鍵,像元尺度在范圍之外時,地形因子的分形計算沒有意義。本研究對研究區(qū)的多重分形特征進行驗證,結果表明尺度對數(shù)滿足3.0≤lnr≤5.3時,黃土高原地區(qū)坡度地形因子的配分函數(shù)對數(shù)與尺度對數(shù)之間的線性關系良好,表明黃土高原區(qū)域上述4個研究區(qū)的坡度信息具有明顯的多重分形特征,可以進行多重分形譜的計算。

    不同地形的多重分形譜圖以及特征參數(shù)都有著明顯的差異,其橫坐標奇異指數(shù)α的大小反映該區(qū)域的奇異程度,這一參數(shù)的數(shù)值越大,表明分形體坡度特征的分異現(xiàn)象越大,其內(nèi)部越不均勻、越復雜,兩極化趨勢越明顯,α最小值、最大值所對應的多重分形譜值記為f(αmin),f(αmax),描述了小坡度單元區(qū)域和大坡度單元區(qū)域在整個研究區(qū)域內(nèi)所占的比例高低。通過多重分形譜的研究可以發(fā)現(xiàn),最小坡度奇異指數(shù)的最小值為1.82,相應的多重分形譜值的最大值為1.78,最大坡度奇異指數(shù)在4個研究區(qū)內(nèi)最大,所對應的多重分形譜在4個研究區(qū)最小為1.57,說明坡度最小地貌單元在黃土塬區(qū)所占比例最小,坡度最大的地貌單元出現(xiàn)在黃土塬區(qū)且占比最高。黃土丘陵區(qū)、河谷平原區(qū)則表現(xiàn)出較小的奇異指數(shù)變化范圍,黃土丘陵區(qū)最小坡度奇異指數(shù)的最小值為1.94,最大值為2.07,河谷平原區(qū)最小坡度奇異指數(shù)的最小值為1.85,最大值為2.19,相比之下地形較為均一。

    通過多重分形譜與統(tǒng)計直方圖對照研究可以發(fā)現(xiàn),在4個研究區(qū)中,黃土塬區(qū)有不少坡度達到了40°,且坡度最大的地貌單元出現(xiàn)在黃土塬區(qū)且占比最高,在地勢平坦的塬面出現(xiàn)這種情況是由于所選研究區(qū)面積固定為30 km×30 km,而單個塬面通常面積有限,在這樣一個較大且規(guī)則的范圍內(nèi)不可避免地包括了塬面周圍的切溝、坡面等非塬面地形,因此其坡度統(tǒng)計中也存在少量坡度大于20°的非塬面特征柵格,這也是黃土塬區(qū)的多重分形譜中奇異指數(shù)最小、對應的多重分形譜值在4個研究區(qū)中最大的原因,同時也反映了研究區(qū)地形復雜的特點。結合多重分形理論的相關知識發(fā)現(xiàn),本研究結果與實際情況有著較好的對應,表明基于多重分形理論與方法對黃土高原地區(qū)的地形因子進行量化研究是可行的,為黃土高原地區(qū)開展地形因子方面的研究提供了參考,對區(qū)域水文過程模擬、水土流失預報模型的建立具有積極意義??蔀辄S土高原的地貌學研究、土壤侵蝕研究和水土保持規(guī)劃工作的開展提供技術支持和幫助。

    猜你喜歡
    黃土塬分形坡度
    感受分形
    全數(shù)字高密度三維地震勘探技術在黃土塬區(qū)的應用研究
    科技視界(2019年16期)2019-07-23 01:50:52
    關于公路超高漸變段合成坡度解析與應用
    分形之美
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應用
    黃土塬上唱黃土謠
    青年歌聲(2017年8期)2017-02-09 01:35:14
    節(jié)水工程在黃土塬地質災害治理中的應用
    基于圖像處理的定位器坡度計算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    壓制黃土塬區(qū)復雜地表條件下折射多次波的組合激發(fā)技術
    坡度在巖石風化層解譯中的應用
    河北遙感(2015年2期)2015-07-18 11:11:14
    午夜影院日韩av| 国产成人a区在线观看| 午夜福利成人在线免费观看| 又黄又爽又刺激的免费视频.| 99热精品在线国产| 久久精品国产亚洲av涩爱 | 乱人视频在线观看| 丝袜美腿在线中文| 国产美女午夜福利| 成熟少妇高潮喷水视频| 天堂av国产一区二区熟女人妻| 老司机深夜福利视频在线观看| 又紧又爽又黄一区二区| .国产精品久久| 亚洲人成网站在线播放欧美日韩| 色在线成人网| 亚洲国产欧美人成| 免费大片18禁| 春色校园在线视频观看| 国产单亲对白刺激| 亚洲va日本ⅴa欧美va伊人久久| 人妻夜夜爽99麻豆av| 91av网一区二区| 国产成年人精品一区二区| 偷拍熟女少妇极品色| 国产亚洲精品久久久com| 国产精品伦人一区二区| 看片在线看免费视频| 成人美女网站在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 大型黄色视频在线免费观看| 欧美黑人欧美精品刺激| 亚洲国产精品成人综合色| 99精品久久久久人妻精品| 久久久国产成人精品二区| 十八禁国产超污无遮挡网站| 免费无遮挡裸体视频| 久久热精品热| 88av欧美| 久久亚洲真实| 三级毛片av免费| 国产精品久久久久久久久免| 中文字幕av在线有码专区| 欧美绝顶高潮抽搐喷水| 国产女主播在线喷水免费视频网站 | 午夜免费激情av| 久99久视频精品免费| 国产久久久一区二区三区| 真人做人爱边吃奶动态| 夜夜看夜夜爽夜夜摸| 色在线成人网| 有码 亚洲区| 成人午夜高清在线视频| 国产视频内射| 免费在线观看日本一区| 欧美三级亚洲精品| 日本色播在线视频| 亚洲性夜色夜夜综合| 99riav亚洲国产免费| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 不卡视频在线观看欧美| 日本熟妇午夜| 色5月婷婷丁香| 在线观看舔阴道视频| 亚洲人成网站高清观看| 久久久久性生活片| 亚洲欧美日韩高清专用| 九九爱精品视频在线观看| 99精品在免费线老司机午夜| 欧美精品啪啪一区二区三区| 中文字幕精品亚洲无线码一区| 久久精品人妻少妇| 黄色日韩在线| 久久久久久久亚洲中文字幕| 亚洲不卡免费看| 老司机福利观看| 国产精品伦人一区二区| 国产精品久久久久久久久免| 蜜桃亚洲精品一区二区三区| 欧美3d第一页| 两人在一起打扑克的视频| 久久精品国产亚洲网站| 99热6这里只有精品| 亚洲五月天丁香| 麻豆av噜噜一区二区三区| 国产成人a区在线观看| 国产黄片美女视频| 婷婷精品国产亚洲av在线| 99热只有精品国产| 夜夜爽天天搞| 国产精品人妻久久久影院| 伦精品一区二区三区| 国产午夜精品论理片| 国产欧美日韩精品一区二区| 久久久久国产精品人妻aⅴ院| 九九在线视频观看精品| 国产亚洲精品av在线| 国产精品日韩av在线免费观看| 3wmmmm亚洲av在线观看| 国产不卡一卡二| 看十八女毛片水多多多| 成人亚洲精品av一区二区| 91久久精品国产一区二区成人| 国产精品久久久久久av不卡| 中亚洲国语对白在线视频| 毛片一级片免费看久久久久 | 亚洲av日韩精品久久久久久密| 国产一区二区在线av高清观看| 成年女人永久免费观看视频| 国产精品日韩av在线免费观看| 我要看日韩黄色一级片| 91久久精品国产一区二区三区| 国产一区二区激情短视频| www.www免费av| 欧美成人免费av一区二区三区| 人人妻人人看人人澡| 国产极品精品免费视频能看的| 成人国产麻豆网| 黄色视频,在线免费观看| 人妻少妇偷人精品九色| 免费看日本二区| 欧美成人免费av一区二区三区| 一本一本综合久久| 免费看av在线观看网站| 亚洲最大成人av| 日本五十路高清| 人妻丰满熟妇av一区二区三区| 乱人视频在线观看| 国产精品,欧美在线| 亚洲人成网站在线播放欧美日韩| 日韩欧美精品免费久久| av天堂在线播放| 老女人水多毛片| 久久午夜福利片| 最近视频中文字幕2019在线8| 91av网一区二区| 全区人妻精品视频| 别揉我奶头 嗯啊视频| 最近视频中文字幕2019在线8| 亚洲av.av天堂| 日本爱情动作片www.在线观看 | 欧美色视频一区免费| 中文字幕免费在线视频6| 国产成人一区二区在线| 日韩欧美国产一区二区入口| 狂野欧美激情性xxxx在线观看| 简卡轻食公司| 国产精品野战在线观看| 22中文网久久字幕| 久久久久久国产a免费观看| 禁无遮挡网站| 中文字幕久久专区| 久久热精品热| 深爱激情五月婷婷| 久久久久性生活片| 精品久久久久久久久久免费视频| 成人毛片a级毛片在线播放| 日日撸夜夜添| 久久午夜福利片| 国产精品永久免费网站| 亚洲av二区三区四区| 97碰自拍视频| 亚洲成人久久性| 国产美女午夜福利| 简卡轻食公司| 亚洲最大成人手机在线| 久久久精品欧美日韩精品| 哪里可以看免费的av片| 久久99热这里只有精品18| 天堂网av新在线| 色视频www国产| 桃色一区二区三区在线观看| 99热这里只有是精品50| 日本黄色视频三级网站网址| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx在线观看| 99在线视频只有这里精品首页| 高清在线国产一区| 国产成人一区二区在线| 色在线成人网| 日本免费a在线| 午夜爱爱视频在线播放| 国产在线男女| 成人av一区二区三区在线看| 久久久久久国产a免费观看| 国产精华一区二区三区| 尤物成人国产欧美一区二区三区| 免费无遮挡裸体视频| 欧美日韩精品成人综合77777| 国产男人的电影天堂91| 高清毛片免费观看视频网站| 日韩欧美国产在线观看| 欧美日韩瑟瑟在线播放| 国产精品av视频在线免费观看| 69av精品久久久久久| 精品日产1卡2卡| 中文资源天堂在线| 神马国产精品三级电影在线观看| 国产精品女同一区二区软件 | 最近最新中文字幕大全电影3| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| 成人三级黄色视频| 国产一区二区三区在线臀色熟女| bbb黄色大片| 国产亚洲精品av在线| 亚洲精品一区av在线观看| 国产一区二区亚洲精品在线观看| 乱系列少妇在线播放| 老女人水多毛片| 成人特级黄色片久久久久久久| 97超视频在线观看视频| 俺也久久电影网| av女优亚洲男人天堂| 亚洲在线自拍视频| 国产三级中文精品| 一个人免费在线观看电影| 亚洲成a人片在线一区二区| 亚洲成av人片在线播放无| 亚洲欧美日韩无卡精品| 久久亚洲真实| 免费观看精品视频网站| 五月伊人婷婷丁香| 国产视频内射| 精品人妻偷拍中文字幕| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 国产精品人妻久久久影院| 精品午夜福利视频在线观看一区| 又粗又爽又猛毛片免费看| 一边摸一边抽搐一进一小说| 成人国产综合亚洲| 黄色视频,在线免费观看| 国产色爽女视频免费观看| 最好的美女福利视频网| 午夜福利视频1000在线观看| 天堂av国产一区二区熟女人妻| 亚洲精品成人久久久久久| 69av精品久久久久久| 成熟少妇高潮喷水视频| 国产av麻豆久久久久久久| 日日啪夜夜撸| 乱人视频在线观看| 精品一区二区三区视频在线| 国产探花在线观看一区二区| 自拍偷自拍亚洲精品老妇| 精品免费久久久久久久清纯| 在线播放国产精品三级| 精品一区二区三区视频在线观看免费| or卡值多少钱| 欧洲精品卡2卡3卡4卡5卡区| 成年女人看的毛片在线观看| 淫秽高清视频在线观看| 国产中年淑女户外野战色| 真人做人爱边吃奶动态| avwww免费| 国产精品久久久久久av不卡| 亚洲在线自拍视频| 久久人人精品亚洲av| 亚洲熟妇熟女久久| 日本撒尿小便嘘嘘汇集6| 中文资源天堂在线| 久久久久久伊人网av| 好男人在线观看高清免费视频| 看黄色毛片网站| av在线老鸭窝| 国产69精品久久久久777片| a级毛片a级免费在线| 波多野结衣高清作品| 中文字幕熟女人妻在线| 欧美黑人巨大hd| 欧美色欧美亚洲另类二区| 男插女下体视频免费在线播放| 久久久久九九精品影院| 毛片女人毛片| 成人av在线播放网站| 成人鲁丝片一二三区免费| 精品久久久久久成人av| 18禁黄网站禁片午夜丰满| 干丝袜人妻中文字幕| 成人一区二区视频在线观看| 99久久九九国产精品国产免费| 色哟哟哟哟哟哟| 精品一区二区三区视频在线观看免费| 99久久久亚洲精品蜜臀av| 久久久久久久午夜电影| 国产单亲对白刺激| 一区二区三区激情视频| 十八禁网站免费在线| 美女 人体艺术 gogo| 天堂av国产一区二区熟女人妻| 午夜福利18| 免费人成视频x8x8入口观看| 日韩大尺度精品在线看网址| 人人妻,人人澡人人爽秒播| 综合色av麻豆| 免费观看的影片在线观看| 在线免费十八禁| 国产高清激情床上av| 日本精品一区二区三区蜜桃| 女同久久另类99精品国产91| 欧美高清性xxxxhd video| 国产综合懂色| 国产伦一二天堂av在线观看| 欧美激情久久久久久爽电影| 久久香蕉精品热| 无遮挡黄片免费观看| 搡老妇女老女人老熟妇| 亚洲在线观看片| 日韩人妻高清精品专区| 亚洲av一区综合| 欧美另类亚洲清纯唯美| 不卡一级毛片| 日韩欧美免费精品| 黄色丝袜av网址大全| 精品一区二区三区视频在线观看免费| 听说在线观看完整版免费高清| 欧美高清性xxxxhd video| 国产免费一级a男人的天堂| 舔av片在线| 欧美日韩亚洲国产一区二区在线观看| 我的老师免费观看完整版| 精品久久久久久久久av| 久久久久久大精品| 高清在线国产一区| 成人国产一区最新在线观看| 久久久久国内视频| 小说图片视频综合网站| 国产精品电影一区二区三区| 亚洲专区中文字幕在线| 99精品久久久久人妻精品| 欧美日韩综合久久久久久 | 国产精品伦人一区二区| av在线观看视频网站免费| 国产人妻一区二区三区在| 不卡一级毛片| 99精品久久久久人妻精品| 我的老师免费观看完整版| 成年人黄色毛片网站| 小说图片视频综合网站| 日本一二三区视频观看| 日本黄色视频三级网站网址| 男人和女人高潮做爰伦理| 免费观看在线日韩| 精品欧美国产一区二区三| 欧美区成人在线视频| 中国美女看黄片| av福利片在线观看| 欧美成人一区二区免费高清观看| 亚洲成av人片在线播放无| 亚洲美女视频黄频| 日本三级黄在线观看| 嫩草影视91久久| 免费人成视频x8x8入口观看| 欧美不卡视频在线免费观看| 久久午夜福利片| 国产黄片美女视频| 性色avwww在线观看| 国产精品98久久久久久宅男小说| 中亚洲国语对白在线视频| 婷婷精品国产亚洲av| 夜夜爽天天搞| 久久香蕉精品热| 久久久久九九精品影院| 欧美三级亚洲精品| av女优亚洲男人天堂| 我的老师免费观看完整版| 日本欧美国产在线视频| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频| 成年女人看的毛片在线观看| 在现免费观看毛片| 国产高清视频在线播放一区| 欧美xxxx性猛交bbbb| 欧美+亚洲+日韩+国产| 在线免费十八禁| 我的老师免费观看完整版| 久9热在线精品视频| 亚洲av五月六月丁香网| 亚洲无线在线观看| 女的被弄到高潮叫床怎么办 | 三级毛片av免费| 欧美日韩综合久久久久久 | 亚洲国产色片| 麻豆成人av在线观看| 亚洲黑人精品在线| 成人永久免费在线观看视频| 一区福利在线观看| 欧美黑人巨大hd| netflix在线观看网站| 亚州av有码| 成人亚洲精品av一区二区| 日本黄大片高清| 亚洲av第一区精品v没综合| 91在线观看av| 欧美一区二区国产精品久久精品| 国产aⅴ精品一区二区三区波| 亚洲精品成人久久久久久| 国产人妻一区二区三区在| 黄色日韩在线| 欧美激情在线99| 亚洲中文日韩欧美视频| 深夜精品福利| 尤物成人国产欧美一区二区三区| 亚洲 国产 在线| 午夜精品在线福利| 欧美不卡视频在线免费观看| 色哟哟·www| 久久人妻av系列| 欧美日韩精品成人综合77777| 中亚洲国语对白在线视频| 真实男女啪啪啪动态图| 美女免费视频网站| 在线看三级毛片| 国内精品美女久久久久久| 成人午夜高清在线视频| 亚洲av中文字字幕乱码综合| 免费电影在线观看免费观看| 99热这里只有精品一区| 欧美三级亚洲精品| 国产精品久久久久久av不卡| 伦精品一区二区三区| 亚洲国产欧洲综合997久久,| 夜夜夜夜夜久久久久| 一级黄色大片毛片| av视频在线观看入口| 少妇的逼水好多| 色吧在线观看| 成年免费大片在线观看| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| 2021天堂中文幕一二区在线观| 国产淫片久久久久久久久| 大型黄色视频在线免费观看| 国产精品久久久久久久久免| 久久精品国产清高在天天线| 国产精品野战在线观看| 国产黄色小视频在线观看| 亚洲av二区三区四区| 禁无遮挡网站| 亚洲av中文av极速乱 | 精品久久久久久久久亚洲 | 女生性感内裤真人,穿戴方法视频| 欧美高清成人免费视频www| 日本爱情动作片www.在线观看 | 亚洲av.av天堂| 啦啦啦观看免费观看视频高清| a级毛片免费高清观看在线播放| 伦精品一区二区三区| 99久久无色码亚洲精品果冻| 日本黄色视频三级网站网址| 午夜激情福利司机影院| 午夜影院日韩av| 免费在线观看成人毛片| 国内精品一区二区在线观看| 亚洲第一区二区三区不卡| av视频在线观看入口| 国产主播在线观看一区二区| 日韩欧美在线二视频| 亚洲自偷自拍三级| 赤兔流量卡办理| 亚洲av第一区精品v没综合| 麻豆成人av在线观看| 午夜精品久久久久久毛片777| 此物有八面人人有两片| 99热这里只有精品一区| 国产亚洲91精品色在线| 国产伦在线观看视频一区| 国产精品久久久久久av不卡| 两个人的视频大全免费| 97超级碰碰碰精品色视频在线观看| 久久国产乱子免费精品| 久久人妻av系列| 久久99热这里只有精品18| 免费在线观看日本一区| 校园人妻丝袜中文字幕| 久久中文看片网| 窝窝影院91人妻| 一个人免费在线观看电影| 亚洲经典国产精华液单| 真人一进一出gif抽搐免费| 日韩欧美国产一区二区入口| 日韩欧美在线二视频| 国产av不卡久久| 在线观看舔阴道视频| 亚洲精华国产精华精| 日韩欧美 国产精品| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 亚洲三级黄色毛片| 一夜夜www| av.在线天堂| 欧美xxxx性猛交bbbb| 国产一区二区在线观看日韩| 黄色欧美视频在线观看| 男女之事视频高清在线观看| 国产视频一区二区在线看| 51国产日韩欧美| 精华霜和精华液先用哪个| 男女那种视频在线观看| 国内毛片毛片毛片毛片毛片| 在线观看66精品国产| 国产淫片久久久久久久久| 亚洲国产精品sss在线观看| 最近最新中文字幕大全电影3| 露出奶头的视频| 国产精品1区2区在线观看.| 精品一区二区三区视频在线| av国产免费在线观看| 精品久久久久久久久久免费视频| 我要搜黄色片| 久久精品91蜜桃| 久久精品国产亚洲网站| 草草在线视频免费看| 免费一级毛片在线播放高清视频| 狠狠狠狠99中文字幕| av.在线天堂| 最新在线观看一区二区三区| 亚洲综合色惰| 免费av观看视频| 男女视频在线观看网站免费| 不卡一级毛片| 久久精品国产99精品国产亚洲性色| 香蕉av资源在线| 国产精品不卡视频一区二区| 欧美xxxx黑人xx丫x性爽| 国产三级中文精品| 欧美高清性xxxxhd video| 亚洲国产欧美人成| 婷婷亚洲欧美| 国产麻豆成人av免费视频| 国产高清激情床上av| 可以在线观看毛片的网站| 亚洲性夜色夜夜综合| 在线天堂最新版资源| 韩国av在线不卡| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 一本久久中文字幕| 最好的美女福利视频网| 老师上课跳d突然被开到最大视频| a级毛片a级免费在线| 亚洲av二区三区四区| 在线免费观看不下载黄p国产 | 日本一二三区视频观看| 深夜a级毛片| 黄色欧美视频在线观看| 女人十人毛片免费观看3o分钟| 日韩欧美 国产精品| 亚洲男人的天堂狠狠| 久久精品综合一区二区三区| 看十八女毛片水多多多| 午夜激情福利司机影院| 国产私拍福利视频在线观看| 国产精品99久久久久久久久| 国产欧美日韩一区二区精品| 日韩精品有码人妻一区| 午夜视频国产福利| 亚洲精品久久国产高清桃花| 日本五十路高清| 久久久久久久亚洲中文字幕| а√天堂www在线а√下载| 国产高清视频在线播放一区| 别揉我奶头 嗯啊视频| 嫩草影院精品99| 动漫黄色视频在线观看| 欧美成人免费av一区二区三区| 国产精品三级大全| 国产高清不卡午夜福利| 啦啦啦啦在线视频资源| 长腿黑丝高跟| 黄色日韩在线| 久久99热6这里只有精品| 精品午夜福利在线看| 天天一区二区日本电影三级| 乱系列少妇在线播放| 久久久成人免费电影| 国产伦一二天堂av在线观看| 麻豆国产av国片精品| 国产主播在线观看一区二区| 中文资源天堂在线| 日韩亚洲欧美综合| 亚洲精品在线观看二区| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 日韩欧美三级三区| 大又大粗又爽又黄少妇毛片口| or卡值多少钱| 国内精品一区二区在线观看| 男女视频在线观看网站免费| av福利片在线观看| 午夜免费激情av| 精品一区二区三区视频在线观看免费| 久久午夜亚洲精品久久| 日本免费一区二区三区高清不卡| 久久久久性生活片| 99久久中文字幕三级久久日本| 免费看光身美女| 三级毛片av免费| 亚洲精品色激情综合| 白带黄色成豆腐渣| 成年免费大片在线观看| 亚洲内射少妇av| .国产精品久久| 成人二区视频| 国产精品日韩av在线免费观看| 熟女电影av网| 一卡2卡三卡四卡精品乱码亚洲| 亚洲狠狠婷婷综合久久图片| av在线蜜桃| 又爽又黄a免费视频| 日日干狠狠操夜夜爽| 欧美又色又爽又黄视频| av在线老鸭窝|