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

    基于形態(tài)指數(shù)的山丘區(qū)中小河流灘地分類方法及演變分析

    2016-09-21 00:40:04伊紫函夏繼紅汪穎俊畢利東余根聽曹偉杰林立懷
    中國水土保持科學 2016年4期
    關鍵詞:龍游縣灘地靈山

    伊紫函,夏繼紅?,汪穎俊,畢利東,余根聽,曹偉杰,林立懷

    (1.河海大學水利水電學院,210098,南京;2.龍游縣水利局,324400,浙江龍游)

    ?

    基于形態(tài)指數(shù)的山丘區(qū)中小河流灘地分類方法及演變分析

    伊紫函1,夏繼紅1?,汪穎俊2,畢利東1,余根聽1,曹偉杰1,林立懷1

    (1.河海大學水利水電學院,210098,南京;2.龍游縣水利局,324400,浙江龍游)

    灘地是河道生態(tài)系統(tǒng)的重要組成部分,具有保護生物多樣性、維持生態(tài)系統(tǒng)的動態(tài)平衡、提高河流的自我修復能力等作用,對河道水土保持與河流健康具有重要意義。為探究中小河流灘地的分類方法,筆者選取浙江省龍游縣靈山港灘地為研究對象,利用ImageJ獲取河道灘地幾何特征參數(shù),在計算圓形率、緊湊度、延伸率、形狀率和平均曲率等基本形態(tài)指數(shù)的基礎上,應用主成分分析法,獲得描述灘地形態(tài)的2個主成分,以其構建平面直角坐標系,根據(jù)形態(tài)指數(shù)數(shù)對所在象限,對灘地的幾何形態(tài)進行定量分析。結果表明:山丘區(qū)中小河流灘地可劃分為短寬規(guī)則型、窄長規(guī)則型、窄長不規(guī)則型與短寬不規(guī)則型等4種類型;分別在2003、2010及2013年,對龍游靈山港各類型灘地的面積與數(shù)量進行統(tǒng)計,發(fā)現(xiàn)靈山港灘地面積減少約17.82%;其中,短寬不規(guī)則型灘地面積減少幅度最大,約達60.76%,且灘地幾何形態(tài)明顯呈窄長化與不規(guī)則化趨勢。本文提出的灘地形態(tài)分類方法,能夠反映灘地的縱橫比例和邊緣規(guī)則程度,并用于分析不同類型灘地的增長與衰退趨勢。

    山丘河道; 灘地; 形態(tài)指數(shù); 分類; 主成分分析

    河道灘地是河道水域與陸域的過渡區(qū)域,具有調蓄洪水、穩(wěn)固河道、凈化水質、生物棲息和資源供給等豐富功能,對河道水土保持生態(tài)建設,以及河流健康保護具有重要意義[1-3]。由于過去不合理的管理、建設、開發(fā)與利用,河道灘地資源萎縮,功能退化,這已引起人們的廣泛關注和重視,一些學者對灘地的類型、發(fā)育、治理與利用等問題開展了研究。例如:李志威等2012年通過對長江、漢江等7條大型河道的653個沙洲衛(wèi)星圖片的研究,將沙洲形態(tài)概化為橢圓形、竹葉形和鐮刀形3種類型[4];高進[5]通過對沙洲形狀的概化研究,提出了沙洲的發(fā)育長度公式;李大鵬等[6]通過對山區(qū)河流灘地的研究,提出了細平整、厚填土碾壓實等治理方案;饒良懿等[7]研究發(fā)現(xiàn)河岸帶種植植被,能有效移除氮、磷、鈣、鉀、硫、鎂等營養(yǎng)物質,以及一些污染物?,F(xiàn)有研究大多針對大江大河中的灘地,以定性研究為主,但對中小河流灘地的定量研究較少;而中小河流已成為我國當前河流建設和管理的重點,尤其是灘地治理和保護更是治理中小河流的難點。針對這一現(xiàn)狀,本文以浙江省龍游縣靈山港灘地為研究對象,討論灘地類型的定量劃分方法,通過形態(tài)指數(shù)分析灘地的變化,以期為灘地治理提供參考。

    1 研究區(qū)概況

    靈山港位于浙江省西部仙霞嶺山系東北緣,介于E119°07′-119°12′,N28°49′-29°03′之間,是錢塘江上游衢江右岸的一條支流,發(fā)源于遂昌縣高坪鄉(xiāng)和尚嶺,流域總面積726.9 km2。靈山港流域屬亞熱帶季風氣候區(qū),雨量豐沛,四季變化明顯。靈山港在龍游縣境內的流域面積為367.6 km2,主流長43.79 km,河道比降為4.24‰。河道自南往北流經(jīng)沐塵、溪口、靈山、官潭、寺后和上圩頭等地,南北貫穿龍游縣城至東北部驛前、湖底葉兩村間匯入衢江。河道上游沐塵村建有沐塵水庫,為大(二)型水庫,流域內年平均降水量為1 815 mm,年平均徑流總量為4.49億m3,年平均徑流深1 130 mm,徑流系數(shù)為0.62。流域內植被較好,森林覆蓋率較高,屬亞熱帶竹林區(qū)和常綠闊葉林區(qū)。該河道在長期的演化過程中,形成了豐富的灘地資源,對河流的健康發(fā)揮了重要作用;但近年來,由于過度的開發(fā)利用,特別是采砂量的增加,導致部分灘地面積減小,破碎化程度增高,功能退化嚴重。

    2 研究方法

    2.1常用平面圖形幾何形態(tài)指標

    平面圖形的形態(tài)通常用圓形率、緊湊度、延伸率、形狀率和平均曲率等指標進行定量計算[8-9]。各指標計算方法如表1所示。

    2.2數(shù)據(jù)獲取方法

    選取2003年靈山港河道實測地形圖、2010及2013年Google Earth影像圖(比例尺為1∶1萬)為數(shù)據(jù)源,應用美國國立衛(wèi)生研究院(National Institutes of Health)開發(fā)的圖像分析工具ImageJ,對相關圖形進行分析與處理[10]。筆者主要應用ImageJ軟件,讀取河道全線灘地的周長、面積、橫徑及縱徑等參數(shù)。首先將相關圖形導入ImageJ,并經(jīng)校準、清噪、灰度化等處理,再利用圈取、抓讀等功能,讀取灘地幾何特征參數(shù),最后存入相應數(shù)據(jù)文件。具體流程如圖1所示。

    表1 常用平面圖形幾何形態(tài)指標的計算公式

    注:A為圖形面積,m2;P為圖形周長,m;A′為該圖形最小外接圓面積,m2;L為圖形最長軸長度,m; L′為圖形最短軸長度,m。Note:A presents the area of the graph,m2; P presents the perimeter of the graph,m; A′ presents the minimum circumscribed circle of the graph,m2; L presents the major axis of the graph,m; L′ presents the minor axis of the graph,m.

    圖1 灘地圖像分析流程Fig.1 Process of floodplain image analysis

    2.3主成分分析法

    主成分分析法是利用降維的思想,將多個指標轉化為幾個綜合指標的數(shù)理統(tǒng)計方法[11-12]。對于n個樣本和p個指標的數(shù)據(jù),可將其構造為矩陣X=(xij)n×p;其中:i=1,2,…,n,j=1,2,…,p。計算X樣本的相關系數(shù)陣R=XTX/(n-1),特征值|R-λIp|=0,可得p 個特征值λ1≥λ2≥…≥λp≥0,各特征值對應的標準化正交特征向量為γ1,γ2,…,γp,應用特征值通過貢獻率=λi/Σλi,計算各個特征值的貢獻率,選取貢獻率85%以上的成分,即為主成分,則第i個主成分為Y=γ1iX1+γ2iX2+…+γpiXp(i=1,2…,p)。本文對圓形率、緊湊度、延伸率、形狀率和平均曲率等指標進行主成分分析,依次算出協(xié)方差矩陣、協(xié)方差矩陣的特征根、特征向量、負荷量、中心化數(shù)據(jù)、各個指標的權重因子及標準化后主成分數(shù)值等,最后確定出描述灘地形態(tài)的主成分。

    3 結果與分析

    3.1灘地形態(tài)指數(shù)

    以2003年靈山港灘地幾何參數(shù)為數(shù)據(jù)源,分別計算出河道沿線各灘地的K1、K2、K3、K4及K5的值。對2003年幾何形態(tài)指標數(shù)據(jù)進行主成分分析,計算出K1、K2、K3、K4及K5構成的協(xié)方差矩陣,由協(xié)方差矩陣計算特征值、特征向量及負荷量(表2和表3)。表2列出了5個指標的特征值、貢獻率和累積貢獻率,成分1和成分2的累積貢獻率達94.01%。其中:成分1的貢獻率為82.20%;成分2的貢獻率為11.81%。選取成分1與成分2作為灘地形態(tài)綜合判斷的2個主成分,分別記為P1和P2。

    表2 協(xié)方差矩陣的特征值及其貢獻率與累積貢獻率

    表3 各指標主成分的特征向量及負荷量

    P1、P2的特征根對應的特征向量,以及每個指標對特征向量的負荷量見表3。可以看出:對P1貢獻最大的是K1、K2、K3和K4,負荷量分別為0.720、0.917、0.940和0.917,則P1基本代表了K1、K2、K3和K4;對P2貢獻最大的是K5,負荷量為0.943,則P2基本代表了K5。因此,P1與P2可用于描述灘地形態(tài),將其定義為灘地形態(tài)指數(shù)向量(P1,P2)。其中:P1主要描述灘地形態(tài)的緊湊性與帶狀程度,反映灘地形態(tài)的“趨圓性”,其值越大灘地越短寬,值越小灘地越窄長;P2主要描述灘地邊界輪廓的不規(guī)則程度,其值越大灘地邊界越規(guī)則,值越小灘地邊界越不規(guī)則。P1與P2值是由K1、K2、K3、K4和K5標準化成ZK1、ZK2、ZK3、ZK4和ZK5后,將各特征向量數(shù)據(jù)中心化復合而得。計算公式如下:

    P1=0.447ZK1+0.485ZK2+0.448ZK3+

    0.485ZK4+0.360ZK5;

    (1)

    P2=0.137ZK1-0.193ZK2-0.410ZK3-

    0.193ZK4+0.861ZK5。

    (2)

    3.2灘地的主要類型

    以P1為橫軸,P2為縱軸構建二維坐標系,應用式(1)、(2)分別計算2003年靈山港沿線灘地P1、P2值,點(P1,P2)分布如圖2??梢钥闯觯械?P1,P2)點分布于P1-P2坐標系的4個象限內,隨著P1值的減小,灘地的形狀變得窄長,隨著P2值的減小,灘地的邊界變得不規(guī)則。據(jù)此本文將分布于第1象限內的灘地歸為短寬規(guī)則型灘地,分布于第2象限內的灘地歸為窄長規(guī)則型灘地,分布于第3象限內的灘地歸為窄長不規(guī)則型灘地,分布于第4象限內的灘地歸為短寬不規(guī)則型灘地,各類型灘地概化圖如圖例所示。

    圖2 灘地類型劃分Fig.2 Floodplain classification analysis

    3.3灘地類型變化分析

    應用上述分類方法,對2010和2013年灘地進行分類,各類型灘地面積和個數(shù)變化見圖3和圖4。由圖3可見,灘地總面積顯著減少,其中,2003年灘地總面積為137.8萬m2,2010年灘地總面積為123.9萬m2,2013年灘地總面積為113.2萬m2。2013與2003年比灘地總面積減少17.82%。其中:短寬不規(guī)則型灘地面積減少最多,約減少60.76%;窄長不規(guī)則型灘地面積約減少5.75%。

    圖3 不同年份各類型灘地面積變化Fig.3 Interannual area proportion change of floodplain

    由圖4看出:2003至2013年,短寬型灘地逐漸減少,窄長型灘地逐漸增多。至2013年時,窄長不規(guī)則型灘地在4類灘地中所占比例最大。2013與2003年相比,窄長規(guī)則型灘地與窄長不規(guī)則灘地分別增加約3.64%與4.40%,短寬規(guī)則型灘地與短寬不規(guī)則型灘地分別減少約0.60%與7.36%。表明靈山港灘地變得狹長,邊界不規(guī)則程度也增加。例如上游沐塵村灘地2003—2013年的變化如圖5所示。2003年時,沐塵村灘地為短寬規(guī)則型灘地,而到2010年、2013年時,該灘地變?yōu)槎虒挷灰?guī)則型灘地??梢?,2003—2013年,該灘地邊緣不規(guī)則程度顯著增加。這主要是由于該灘地位于沐塵水庫下游,且距灘地頭部500 m建有堰壩,這些建筑物使該灘地邊緣水流紊動性增強,造成灘地邊緣不規(guī)則性增強。另外,梅村灘地2003—2013年的形態(tài)變化如圖6所示,梅村灘地由2003年的整塊灘地變?yōu)?010年和2013年的相互分離的小灘地,面積減少約60%。2003年以來,梅村灘地邊界越來越不規(guī)則,2003年時梅村灘地為窄長規(guī)則型,而到2010年灘地離散成短寬不規(guī)則型(①號灘地)、窄長不規(guī)則型(②號灘地)和窄長規(guī)則型(③號灘地),到2013年灘地演變?yōu)檎L不規(guī)則型(④號灘地)、短寬不規(guī)則型(⑤號灘地)和窄長規(guī)則型(⑥號灘地)。這主要是由于在2003—2010年間,該河段內受到采砂等人為因素的擾動,使灘地局部沖刷深度加大,甚至出現(xiàn)深坑,最終導致灘地支離破碎。

    圖4 不同年份各類型灘地數(shù)量比例變化Fig.4 Interannual quantity proportion change of floodplain

    圖5 不同年份沐塵村灘地演變Fig.5 Interannual change of floodplain in Muchen village

    圖6 不同年份梅村灘地演變Fig.6 Interannual change of floodplain in Mei village

    4 結論

    1)以浙江省龍游縣靈山港為研究對象,應用主成分分析法,通過復合描述平面形態(tài)的指標,構建了描述灘地形態(tài)的形態(tài)指數(shù)向量,并建立了由形態(tài)指數(shù)向量為坐標軸的坐標系,根據(jù)形態(tài)指數(shù)在坐標系中的分布象限情況,將灘地分為短寬規(guī)則型、窄長規(guī)則型、窄長不規(guī)則型與短寬不規(guī)則型4種類型。與大型河道的灘地形狀概化分類方式相比,本文提出的山丘區(qū)中小河流灘地的分類方法,既能反映灘地縱橫比例,也能反映灘地邊緣規(guī)則程度,適用于形態(tài)多樣的中小河流灘地。該分類方法為山丘中小河道規(guī)劃、設計和管理提供了一定依據(jù)。

    2)在灘地類型的動態(tài)變化研究中,各類型灘地的數(shù)量和面積變化,能夠準確的反映不同類型灘地的增長與衰退。分析表明:靈山港灘地總面積呈下降趨勢,尤其是短寬型灘地面積急劇下降,灘地形態(tài)總體呈現(xiàn)窄長化與不規(guī)則化,導致靈山港河勢不太穩(wěn)定;因此,需根據(jù)河道來水來沙和人為擾動情況,采取有效生態(tài)治理措施,制定科學的水庫調度運行制度,才能有效遏制灘地破碎化和功能退化的趨勢。

    [1]夏繼紅,林俊強,陳永明,等.國外河流潛流層研究的發(fā)展過程及研究方法[J].水利水電科技進展,2013,33(4):73.

    Xia Jihong,Lin Junqiang,Chen Yongming,et al.Advancing processes and methods of abroad research on hyporheic zone [J].Advances in Science and Technology of Water Resources,2013,33(4):73.(in Chinese)

    [2]夏繼紅,林俊強,姚莉,等.河岸帶的邊緣結構特征與邊緣效應[J].河海大學學報(自然科學版),2010,38(2):215.

    Xia Jihong,Lin Junqiang,Yao Li,et al.Edge structure and edge effect of riparian zones[J].Journal of Hohai University (Natural Sciences),2010,38(2):215.(in Chinese)

    [3]夏繼紅,陳永明,王為木,等.河岸帶潛流層動態(tài)過程與生態(tài)修復[J].水科學進展,2013,24(4):589.

    Xia Jihong,Chen Yongming,Wang Weimu,et al.Dynamic processes and ecological restoration of hyporheic layer in riparian zone [J].Advances in Water Science,2013,24(4):589.(in Chinese)

    [4]李志威,王兆印,張康.典型沙洲形態(tài)與河道的關系[J].泥沙研究,2012,(1):68.

    Li Zhiwei,Wang Zhaoyi,Zhang Kang.Relationship between morphology of typical sand bars and river channels [J].Journal of Sediment Research,2012,(1):68.(in Chinese)

    [5]高進.河流沙洲發(fā)育的理論分析[J].水利學報,1999(6):66.

    Gao Jin.On the formation of shoal in river [J].Journal of Hydraulic Engineering,1999(6):66.(in Chinese)

    [6]李大鵬,王靜.淺談山區(qū)河流及灘涂綜合治理方法[J].農業(yè)與技術,2007,27(2):107.

    Li Dapeng,Wang Jing.Introduction to mountain rivers and beaches comprehensive treatment method [J].Agriculture & Technology,2007,27(2):107.(in Chinese)

    [7]饒良懿,崔建國.河岸植被緩沖帶生態(tài)水文功能研究進展[J].中國水土保持科學,2008,6(4):121.

    Rao Liangyi,Cui Jianguo.Research advances on the eco-hydrological functions of riparian buffer [J].Science of Soil and Water Conservation,2008,6(4):121.(in Chinese)

    [8]潘竟虎,韓文超.近20a中國省會及以上城市空間形態(tài)演變[J].自然資源學報,2013,28(3):470.

    Pan Jinghu,Han Chaowen.Spatial-temporal changes of urban morphology of provincial capital cities or above in China [J].Journal of Natural Researches,2013,28(3):470.(in Chinese)

    [9]張鳳太,王臘春,冷輝,等.典型天然與人工湖泊形態(tài)特征比較分析[J].中國農村水利水電,2012(7):38.

    Zhang Fengtai,Wang Lachun,Leng Hui,et al.A comparative analysis of morphological characteristics of typical natural and artificial lakes [J].China Rural Water and Hydropower,2012(7):38.(in Chinese)

    [10] 畢利東,張斌,潘繼花.運用ImageJ軟件分析土壤結構特征[J].土壤,2009,41(4):654.

    Bi Lidong,Zhang Bin,Pan Jihua.Analysis of soil structural properties by using Image-J software [J].Soils,2009,41(4):654.(in Chinese)

    [11] Sebok E,Duque C,Engesgaard P,et al.Spatial variability in streambed hydraulic conductivity of contrasting stream morphologies:channel bend and straight channel [J].Hydrological Processes,2015,29(3):458.

    [12] 趙宏亮,陳凱,張強,等.應用主成分分析和聚類分析的水稻源庫特性研究[J].沈陽農業(yè)大學學報,2015(2):135.

    Zhao Hongliang,Chen Kai,Zhang Qiang,et al.Application of principal component analysis and cluster analysis of source-sink characteristics of rice [J].Journal of Shenyang Agricultural University,2015(2):135.(in Chinese)

    Morphological indexes-based quantitative classification of floodplains in hilly areas

    Yi Zihan1,Xia Jihong1,Wang Yingjun2,Bi Lidong1,Yu Genting1,Cao Weijie1,Lin Lihuai1

    (1.College of Water Conservancy and Hydropower Engineering,Hohai University,210098,Nanjing,China;2.Longyou Water Conservancy Bureau,324400,Longyou,Zhejiang,China)

    [Background] The floodplain which is an important part of the river ecosystem,can protect the biodiversity,decrease the over-exploitation of sources,keep the homeostasis of the ecological system,and improve the abilities of the self-adjustment and self-healing of the system for the improvement of the human living environment,then has an important significance for the river water and soil conservation and the river health protection.In order to prove the classification method of medium or small river floodplains in hilly areas,this paper takes the geometric feature parameters of the Lingshan River floodplains in Longyou County of Zhejiang Province as research objects.[Methods] The geometric feature parameters of the Lingshan River floodplains were obtained through ImageJ software.And based on the basic morphology indexes,including circularity ratio,compact ration,elongation,form ratio,and mean curvature,the principal component analysis was used to get two principal components,P1 and P2,the two principal components can represent the information of all the basic morphology indexes.Using the two principal components as the plane coordinate axes and according to their values,the quantitative analysis of floodplain geometrical morphology was carried out.[Results] The research result indicated that the floodplains were classified into four types:the short wide regular floodplains,the narrow long regular floodplains,the narrow long irregular floodplains,and the short wide irregular floodplains.Analyzing the area and amount of floodplains in 2003,2010 and 2013,it could be concluded that the area of floodplains decreased by about 17.82% in Lingshan River from 2003 to 2013,and the area of short wide irregular floodplains declined fastest by about 60.76%.The geometrical morphology of floodplain obviously presented the narrow long and irregular trend.Then,taking the Muchen village floodplain and Mei village floodplain as examples,the Muchen village floodplain was short wide and regular in 2003.However,in 2010 and 2013,it became a short wide irregular floodplain.This is mainly because that the Muchen village floodplain is located in the lower reaches of Muchen reservoir.The dam enhanced the turbulent fluctuation of the stream of the floodplain fringe and the irregularity of the floodplain fringe.The whole Mei village floodplain in 2003 was separated from some unrelatedly small floodplains in 2010 and 2013,and the whole area reduced by about 60%.Since 2003,due to the effects of human activities and excessive utilization,the Mei village floodplain fringe has been becoming more and more irregular.[Conclusions] Therefore,the classification method proposed in this paper can reflect the length and width rate and the degree of edge rules of the floodplains.In addition,it can be applied to analyze the morphology change tendency of floodplains in medium or small rivers.

    hilly stream; floodplain; morphological index; classification; principal component analysis

    2016-04-11

    2016-06-27

    項目名稱:國家自然科學基金“蜿蜒型河岸帶潛流層水動力學機制及溶質運移規(guī)律研究”(41471069);浙江省水利科技項目“龍游縣中小河流灘地時空演化機理及生態(tài)修復技術研究”(RC1527)

    伊紫函(1992—),女,碩士研究生。主要研究方向:河流生態(tài)建設。E-mail:guodonghan143@sohu.com

    簡介:夏繼紅(1970—),男,博士,教授。主要研究方向:河流健康理論與生態(tài)修復。E-mail:syjhxia@hhu.edu.cn

    TV147

    A

    1672-3007(2016)04-0128-06

    10.16843/j.sswc.2016.04.016

    猜你喜歡
    龍游縣灘地靈山
    浙江省典型河流灘地時空演變分析*
    彎道灘地植被作用下漫灘洪水水沙因子橫向分布研究
    洞庭湖區(qū)灘地及稻田洼地資源分布研究
    靈山勝景
    寶藏(2022年1期)2022-08-01 02:12:42
    龍游縣冬季霧霾污染過程中碳質氣溶膠特征及區(qū)域傳輸分析
    能源工程(2021年3期)2021-08-05 07:26:12
    與教材“對話” 讓復習 “精準”——從回歸教材視角談高三數(shù)學復習“五部曲”
    山丘區(qū)中小河流灘地時空演化與擾動因子分析
    叫響靈山灣
    商周刊(2017年11期)2017-06-13 07:32:30
    靈山道中即景(外四首)
    岷峨詩稿(2017年4期)2017-04-20 06:26:32
    龍游縣種植紅肉蜜柚的氣候可行性分析
    一级毛片女人18水好多| 久久久久亚洲av毛片大全| 亚洲国产精品sss在线观看| 中文字幕人成人乱码亚洲影| 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 久久久久久久久免费视频了| 亚洲精品在线美女| 99在线人妻在线中文字幕| 久久久久九九精品影院| 国产av一区在线观看免费| 久久精品亚洲熟妇少妇任你| 在线视频色国产色| 最近最新中文字幕大全免费视频| 成人国语在线视频| 中文字幕精品免费在线观看视频| 好男人电影高清在线观看| 免费av毛片视频| 夜夜爽天天搞| 女人精品久久久久毛片| tocl精华| 老司机在亚洲福利影院| 亚洲av片天天在线观看| 亚洲精品一卡2卡三卡4卡5卡| 91老司机精品| 亚洲av片天天在线观看| 在线观看舔阴道视频| 精品卡一卡二卡四卡免费| 制服人妻中文乱码| 欧美日韩亚洲国产一区二区在线观看| 国内精品久久久久久久电影| 国产极品粉嫩免费观看在线| 亚洲欧美精品综合久久99| 欧美久久黑人一区二区| 亚洲 欧美一区二区三区| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美精品综合一区二区三区| 精品电影一区二区在线| 精品国产一区二区三区四区第35| 国产精品久久视频播放| 色尼玛亚洲综合影院| 美女国产高潮福利片在线看| 神马国产精品三级电影在线观看 | 国语自产精品视频在线第100页| 亚洲第一电影网av| 欧美国产日韩亚洲一区| 老司机午夜福利在线观看视频| 国产精品永久免费网站| 黄片播放在线免费| 热re99久久国产66热| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久人人人人人| 午夜福利18| 黄片小视频在线播放| 一个人免费在线观看的高清视频| 色播亚洲综合网| 欧美色视频一区免费| 99久久综合精品五月天人人| 中文字幕久久专区| 亚洲午夜精品一区,二区,三区| 亚洲中文字幕日韩| 久久欧美精品欧美久久欧美| 日本欧美视频一区| 日日摸夜夜添夜夜添小说| 成人三级做爰电影| 91大片在线观看| 欧美不卡视频在线免费观看 | 一级作爱视频免费观看| a在线观看视频网站| 亚洲人成电影免费在线| 亚洲国产精品久久男人天堂| 在线国产一区二区在线| 一边摸一边抽搐一进一出视频| 日韩免费av在线播放| 国产蜜桃级精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 久久久久久久久久久久大奶| 18禁黄网站禁片午夜丰满| 欧美另类亚洲清纯唯美| 日韩大尺度精品在线看网址 | 18禁美女被吸乳视频| 别揉我奶头~嗯~啊~动态视频| 成熟少妇高潮喷水视频| 九色国产91popny在线| 熟妇人妻久久中文字幕3abv| 老汉色∧v一级毛片| 国产精品久久久人人做人人爽| 久久精品91无色码中文字幕| 欧美最黄视频在线播放免费| 亚洲精品中文字幕在线视频| 91精品三级在线观看| 真人一进一出gif抽搐免费| 一级毛片精品| 免费看a级黄色片| 久久人妻福利社区极品人妻图片| 校园春色视频在线观看| 一级a爱视频在线免费观看| 国产主播在线观看一区二区| 午夜免费激情av| 国产在线精品亚洲第一网站| www.自偷自拍.com| 极品教师在线免费播放| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 天堂√8在线中文| 在线观看免费午夜福利视频| 在线国产一区二区在线| 婷婷六月久久综合丁香| 一级作爱视频免费观看| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 丝袜人妻中文字幕| 日日摸夜夜添夜夜添小说| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 在线播放国产精品三级| 男人舔女人的私密视频| 国产伦人伦偷精品视频| 亚洲国产精品久久男人天堂| 在线观看www视频免费| 19禁男女啪啪无遮挡网站| 很黄的视频免费| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 很黄的视频免费| 在线十欧美十亚洲十日本专区| 两个人免费观看高清视频| 高清毛片免费观看视频网站| 精品卡一卡二卡四卡免费| 美女免费视频网站| 麻豆av在线久日| 黄色a级毛片大全视频| 纯流量卡能插随身wifi吗| 叶爱在线成人免费视频播放| 制服人妻中文乱码| 欧美日本亚洲视频在线播放| 桃红色精品国产亚洲av| 国产亚洲精品av在线| 国产又爽黄色视频| 国产熟女午夜一区二区三区| 国产精品日韩av在线免费观看 | 亚洲久久久国产精品| 精品国产超薄肉色丝袜足j| 久久精品91蜜桃| 免费av毛片视频| а√天堂www在线а√下载| 亚洲成av人片免费观看| 久久这里只有精品19| 中文字幕最新亚洲高清| 不卡av一区二区三区| 伊人久久大香线蕉亚洲五| 18禁国产床啪视频网站| 久久精品国产亚洲av高清一级| 久久久精品国产亚洲av高清涩受| 色av中文字幕| 在线观看舔阴道视频| 午夜视频精品福利| 一二三四社区在线视频社区8| 久久精品成人免费网站| 中出人妻视频一区二区| 亚洲av片天天在线观看| 99久久久亚洲精品蜜臀av| 久久精品人人爽人人爽视色| www.精华液| 丝袜在线中文字幕| 51午夜福利影视在线观看| 正在播放国产对白刺激| 无遮挡黄片免费观看| a级毛片在线看网站| 中文字幕久久专区| 久久亚洲精品不卡| 很黄的视频免费| 欧美精品亚洲一区二区| 淫秽高清视频在线观看| 天堂影院成人在线观看| 亚洲少妇的诱惑av| 给我免费播放毛片高清在线观看| 757午夜福利合集在线观看| 88av欧美| 婷婷丁香在线五月| 18禁观看日本| 亚洲精品久久国产高清桃花| 久久久久久久久中文| 丁香欧美五月| bbb黄色大片| 少妇熟女aⅴ在线视频| 欧美激情 高清一区二区三区| 三级毛片av免费| 欧美日韩瑟瑟在线播放| videosex国产| 淫妇啪啪啪对白视频| av天堂在线播放| 18美女黄网站色大片免费观看| 日日摸夜夜添夜夜添小说| 精品国产一区二区三区四区第35| 国产成人一区二区三区免费视频网站| 国产亚洲精品久久久久久毛片| 亚洲无线在线观看| 久久人妻av系列| av天堂在线播放| 久久久久国内视频| 欧美性长视频在线观看| 亚洲精品国产精品久久久不卡| 久久久久国产一级毛片高清牌| 日韩中文字幕欧美一区二区| 高清黄色对白视频在线免费看| 亚洲在线自拍视频| 男女做爰动态图高潮gif福利片 | 一级黄色大片毛片| 免费女性裸体啪啪无遮挡网站| 久久午夜综合久久蜜桃| 巨乳人妻的诱惑在线观看| 欧美乱妇无乱码| 国产精品野战在线观看| 1024香蕉在线观看| 国产单亲对白刺激| 老司机靠b影院| 少妇裸体淫交视频免费看高清 | 久久精品91无色码中文字幕| 国产伦一二天堂av在线观看| 国产精品二区激情视频| 色av中文字幕| 国产亚洲精品第一综合不卡| 亚洲五月婷婷丁香| 九色亚洲精品在线播放| 免费高清在线观看日韩| 成人欧美大片| 亚洲第一青青草原| www.精华液| 亚洲精品在线观看二区| av福利片在线| 动漫黄色视频在线观看| 国产高清有码在线观看视频 | 亚洲人成伊人成综合网2020| 亚洲人成网站在线播放欧美日韩| 高清黄色对白视频在线免费看| 黄色毛片三级朝国网站| 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 神马国产精品三级电影在线观看 | 国产精品久久久久久人妻精品电影| 国产成人免费无遮挡视频| 精品人妻在线不人妻| 老司机午夜十八禁免费视频| 一区二区三区激情视频| 国产精品久久久久久人妻精品电影| 久久国产亚洲av麻豆专区| 他把我摸到了高潮在线观看| 少妇的丰满在线观看| 不卡一级毛片| 久久久久国内视频| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 黄色a级毛片大全视频| 一a级毛片在线观看| 亚洲黑人精品在线| 国产高清videossex| 国产精品亚洲美女久久久| 天堂影院成人在线观看| 中文字幕久久专区| 日韩免费av在线播放| 国产亚洲欧美在线一区二区| 国产精品亚洲av一区麻豆| 在线十欧美十亚洲十日本专区| 久久婷婷成人综合色麻豆| 长腿黑丝高跟| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 午夜福利18| 啪啪无遮挡十八禁网站| 亚洲第一电影网av| 午夜福利高清视频| 少妇裸体淫交视频免费看高清 | 亚洲成人久久性| 黑丝袜美女国产一区| 国产精品 欧美亚洲| 激情视频va一区二区三区| 亚洲中文av在线| 首页视频小说图片口味搜索| 一二三四在线观看免费中文在| 99国产精品免费福利视频| 成人三级黄色视频| 老鸭窝网址在线观看| 久久久久久久久免费视频了| 97碰自拍视频| 免费看美女性在线毛片视频| 黄色丝袜av网址大全| 日本黄色视频三级网站网址| 一二三四社区在线视频社区8| x7x7x7水蜜桃| 757午夜福利合集在线观看| 日韩欧美一区视频在线观看| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 久久精品成人免费网站| 少妇 在线观看| 少妇裸体淫交视频免费看高清 | 精品无人区乱码1区二区| 麻豆久久精品国产亚洲av| 久久久久国内视频| 黑人巨大精品欧美一区二区mp4| 久久人妻熟女aⅴ| 亚洲国产欧美网| 一二三四社区在线视频社区8| 级片在线观看| 国产私拍福利视频在线观看| 久久久久精品国产欧美久久久| 精品久久久久久,| 很黄的视频免费| 午夜福利高清视频| bbb黄色大片| 国产成人欧美| tocl精华| 老司机靠b影院| 大香蕉久久成人网| 亚洲三区欧美一区| 在线观看午夜福利视频| 一边摸一边做爽爽视频免费| 亚洲人成电影观看| 久久香蕉精品热| 久久亚洲真实| 欧美日韩精品网址| 黄色视频不卡| 成人三级做爰电影| 精品久久久久久久人妻蜜臀av | 91九色精品人成在线观看| 亚洲国产毛片av蜜桃av| 精品国产乱码久久久久久男人| 免费久久久久久久精品成人欧美视频| 中文亚洲av片在线观看爽| 91麻豆精品激情在线观看国产| 男人舔女人的私密视频| videosex国产| 最新美女视频免费是黄的| 动漫黄色视频在线观看| 88av欧美| 久久性视频一级片| 窝窝影院91人妻| 91成人精品电影| 久久久久久久久久久久大奶| 高潮久久久久久久久久久不卡| 久久国产乱子伦精品免费另类| 18美女黄网站色大片免费观看| 热re99久久国产66热| 国产高清激情床上av| 女警被强在线播放| 欧美日韩精品网址| 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产亚洲在线| 深夜精品福利| 久久午夜亚洲精品久久| 亚洲av成人一区二区三| 成人av一区二区三区在线看| 极品教师在线免费播放| 成人18禁高潮啪啪吃奶动态图| 91成人精品电影| 给我免费播放毛片高清在线观看| 欧美乱色亚洲激情| 亚洲第一青青草原| 老司机午夜十八禁免费视频| 国产精品1区2区在线观看.| 99国产精品一区二区三区| 久久影院123| 欧美日韩黄片免| 亚洲国产精品成人综合色| 日日夜夜操网爽| 久久精品aⅴ一区二区三区四区| 真人一进一出gif抽搐免费| 韩国av一区二区三区四区| 国产精品免费一区二区三区在线| 国产成人av教育| 又紧又爽又黄一区二区| 国产精品日韩av在线免费观看 | 免费无遮挡裸体视频| 人人妻人人爽人人添夜夜欢视频| 成年女人毛片免费观看观看9| 丝袜在线中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 在线观看午夜福利视频| 在线观看舔阴道视频| av视频免费观看在线观看| 美女 人体艺术 gogo| 一级黄色大片毛片| 极品教师在线免费播放| 黄色a级毛片大全视频| 国产不卡一卡二| 久99久视频精品免费| 国产亚洲av高清不卡| а√天堂www在线а√下载| 一级毛片精品| 国产私拍福利视频在线观看| 婷婷丁香在线五月| 神马国产精品三级电影在线观看 | 亚洲熟女毛片儿| 大陆偷拍与自拍| 欧美一级a爱片免费观看看 | 久久草成人影院| 黄色 视频免费看| 黄色视频,在线免费观看| 精品欧美一区二区三区在线| 精品不卡国产一区二区三区| 琪琪午夜伦伦电影理论片6080| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区中文字幕在线| 精品免费久久久久久久清纯| 国产日韩一区二区三区精品不卡| 女人被躁到高潮嗷嗷叫费观| 黄色视频,在线免费观看| 伦理电影免费视频| 久久久久久人人人人人| 日本免费a在线| 成人18禁高潮啪啪吃奶动态图| 丁香欧美五月| 国产亚洲欧美98| 久久久久久人人人人人| av网站免费在线观看视频| 亚洲专区字幕在线| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 热99re8久久精品国产| 久久国产精品影院| 免费高清在线观看日韩| 精品国内亚洲2022精品成人| 精品乱码久久久久久99久播| 精品久久久久久成人av| 纯流量卡能插随身wifi吗| 老鸭窝网址在线观看| 在线观看舔阴道视频| 久久久久精品国产欧美久久久| 日韩欧美国产一区二区入口| tocl精华| 黄网站色视频无遮挡免费观看| 国产成人精品无人区| 99久久99久久久精品蜜桃| 变态另类丝袜制服| 国产日韩一区二区三区精品不卡| 久久狼人影院| 亚洲一区高清亚洲精品| 韩国精品一区二区三区| 香蕉国产在线看| 在线观看免费视频日本深夜| cao死你这个sao货| 一本大道久久a久久精品| 美女高潮到喷水免费观看| 亚洲在线自拍视频| 日韩欧美一区视频在线观看| 法律面前人人平等表现在哪些方面| 一本综合久久免费| 久久人人97超碰香蕉20202| 18禁裸乳无遮挡免费网站照片 | 亚洲avbb在线观看| 麻豆久久精品国产亚洲av| 一边摸一边抽搐一进一小说| 成人国产综合亚洲| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人一区二区三| 久久香蕉国产精品| 一级毛片精品| АⅤ资源中文在线天堂| 久久午夜综合久久蜜桃| av视频免费观看在线观看| 国产亚洲欧美精品永久| 亚洲美女黄片视频| 亚洲欧美激情综合另类| 国产精品久久久久久亚洲av鲁大| 一进一出抽搐动态| 一进一出抽搐gif免费好疼| av网站免费在线观看视频| 757午夜福利合集在线观看| 露出奶头的视频| 亚洲情色 制服丝袜| 久久天堂一区二区三区四区| 国产单亲对白刺激| 妹子高潮喷水视频| 18禁裸乳无遮挡免费网站照片 | 啪啪无遮挡十八禁网站| 亚洲一区中文字幕在线| av视频在线观看入口| 国产伦一二天堂av在线观看| 亚洲性夜色夜夜综合| 男男h啪啪无遮挡| 亚洲精品粉嫩美女一区| 中文字幕久久专区| 长腿黑丝高跟| 亚洲欧美日韩高清在线视频| 精品乱码久久久久久99久播| 夜夜躁狠狠躁天天躁| 变态另类丝袜制服| 一区二区三区国产精品乱码| 久久久国产精品麻豆| 国产熟女xx| 欧美激情高清一区二区三区| 在线十欧美十亚洲十日本专区| cao死你这个sao货| 一级黄色大片毛片| 在线观看免费日韩欧美大片| 亚洲人成网站在线播放欧美日韩| 久久午夜亚洲精品久久| 女人被躁到高潮嗷嗷叫费观| 国产成人一区二区三区免费视频网站| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 人人妻,人人澡人人爽秒播| x7x7x7水蜜桃| 色在线成人网| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费视频网站a站| 亚洲男人的天堂狠狠| 他把我摸到了高潮在线观看| 亚洲天堂国产精品一区在线| 岛国视频午夜一区免费看| 久久青草综合色| 成人亚洲精品一区在线观看| 国产亚洲精品av在线| 亚洲欧美一区二区三区黑人| 777久久人妻少妇嫩草av网站| 一区二区三区精品91| 老司机靠b影院| 国产成人av教育| 亚洲自拍偷在线| 少妇的丰满在线观看| 国产又爽黄色视频| 高清在线国产一区| 亚洲熟妇中文字幕五十中出| 黄色毛片三级朝国网站| 天天添夜夜摸| 可以在线观看毛片的网站| 长腿黑丝高跟| 久久国产亚洲av麻豆专区| 久久精品人人爽人人爽视色| 99久久99久久久精品蜜桃| 国产精品日韩av在线免费观看 | 国产亚洲欧美在线一区二区| 男男h啪啪无遮挡| 国产成人精品在线电影| 女人爽到高潮嗷嗷叫在线视频| 男女之事视频高清在线观看| 激情在线观看视频在线高清| 午夜福利在线观看吧| 一区二区三区激情视频| 久久久国产成人免费| 成人国产一区最新在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲国产毛片av蜜桃av| 熟女少妇亚洲综合色aaa.| 精品人妻在线不人妻| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久久久99蜜臀| 免费在线观看日本一区| 成人免费观看视频高清| 好看av亚洲va欧美ⅴa在| 搡老妇女老女人老熟妇| 精品福利观看| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 一区在线观看完整版| 99国产极品粉嫩在线观看| 精品午夜福利视频在线观看一区| 久久婷婷人人爽人人干人人爱 | 免费在线观看完整版高清| 淫妇啪啪啪对白视频| 99在线人妻在线中文字幕| 亚洲av成人一区二区三| 国产高清有码在线观看视频 | 18禁黄网站禁片午夜丰满| 久久婷婷成人综合色麻豆| 亚洲免费av在线视频| 97人妻精品一区二区三区麻豆 | 欧美丝袜亚洲另类 | 免费高清视频大片| 丝袜美足系列| 男女之事视频高清在线观看| 亚洲成人免费电影在线观看| 无人区码免费观看不卡| 欧美在线黄色| 国产日韩一区二区三区精品不卡| 琪琪午夜伦伦电影理论片6080| 欧美 亚洲 国产 日韩一| 欧美成人一区二区免费高清观看 | 老司机靠b影院| 日韩一卡2卡3卡4卡2021年| 日韩欧美一区视频在线观看| 我的亚洲天堂| 天天添夜夜摸| 国产精品98久久久久久宅男小说| 精品欧美国产一区二区三| 国产精品影院久久| 午夜久久久在线观看| 亚洲,欧美精品.| 啦啦啦韩国在线观看视频| 国产黄a三级三级三级人| 99久久99久久久精品蜜桃| 高清黄色对白视频在线免费看| 在线观看免费视频日本深夜| 99精品久久久久人妻精品| 国产真人三级小视频在线观看| 伦理电影免费视频| 亚洲中文av在线| 欧美性长视频在线观看| 高清毛片免费观看视频网站| 午夜精品国产一区二区电影| 久久久国产成人精品二区| 国产精品香港三级国产av潘金莲| 伦理电影免费视频| 变态另类丝袜制服| 久久精品国产综合久久久| 大型av网站在线播放| 国产视频一区二区在线看| 别揉我奶头~嗯~啊~动态视频| 悠悠久久av| 欧美丝袜亚洲另类 | 国产精品精品国产色婷婷| 麻豆成人av在线观看|