• <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
    龍游縣種植紅肉蜜柚的氣候可行性分析
    久久精品国产a三级三级三级| 成人鲁丝片一二三区免费| av在线观看视频网站免费| 午夜视频国产福利| 国产综合精华液| 夜夜看夜夜爽夜夜摸| 国产免费一级a男人的天堂| 黄色配什么色好看| 国产成人freesex在线| 中文资源天堂在线| 91狼人影院| 国产久久久一区二区三区| 51国产日韩欧美| 午夜爱爱视频在线播放| 人妻少妇偷人精品九色| 亚洲国产色片| av女优亚洲男人天堂| 最近2019中文字幕mv第一页| 99热这里只有是精品在线观看| 国产精品99久久99久久久不卡 | 三级国产精品片| 嫩草影院新地址| 久久久久久伊人网av| 女人久久www免费人成看片| 亚洲av不卡在线观看| 97超视频在线观看视频| 亚洲欧美一区二区三区黑人 | 在线看a的网站| 国产在线一区二区三区精| 亚洲,一卡二卡三卡| 亚洲精品日韩av片在线观看| 国产成人freesex在线| 国产精品熟女久久久久浪| 免费av观看视频| av卡一久久| 国产女主播在线喷水免费视频网站| 国产高清三级在线| 中文字幕亚洲精品专区| 一本一本综合久久| 午夜日本视频在线| 国产白丝娇喘喷水9色精品| 亚洲经典国产精华液单| 熟女电影av网| 大香蕉久久网| 欧美变态另类bdsm刘玥| 一区二区三区四区激情视频| 18禁在线播放成人免费| 色哟哟·www| 精品久久久久久久久av| 亚洲av福利一区| 亚洲伊人久久精品综合| 久久99蜜桃精品久久| 王馨瑶露胸无遮挡在线观看| 国产一区二区在线观看日韩| 久久久久久久国产电影| 丝袜脚勾引网站| 国产精品久久久久久精品电影小说 | 国产精品精品国产色婷婷| 日韩三级伦理在线观看| 一边亲一边摸免费视频| 亚洲自偷自拍三级| 黑人高潮一二区| 国产成人91sexporn| 午夜激情久久久久久久| 亚洲自偷自拍三级| 国产亚洲午夜精品一区二区久久 | 亚洲精品乱码久久久v下载方式| 狠狠精品人妻久久久久久综合| 热99国产精品久久久久久7| 国产精品无大码| 欧美精品一区二区大全| 亚洲人成网站高清观看| 韩国高清视频一区二区三区| 亚洲综合精品二区| 欧美性感艳星| 国产有黄有色有爽视频| 日韩 亚洲 欧美在线| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| av在线老鸭窝| a级毛色黄片| 国产乱来视频区| 国产精品爽爽va在线观看网站| a级一级毛片免费在线观看| 三级国产精品片| 精华霜和精华液先用哪个| 777米奇影视久久| 日本熟妇午夜| 亚洲美女搞黄在线观看| 国产高清有码在线观看视频| www.色视频.com| 亚洲欧美一区二区三区黑人 | 国产女主播在线喷水免费视频网站| 人妻 亚洲 视频| 性色av一级| 欧美成人a在线观看| 久久久久久伊人网av| 亚洲av男天堂| av又黄又爽大尺度在线免费看| 九色成人免费人妻av| 听说在线观看完整版免费高清| 69人妻影院| 日韩欧美精品免费久久| 麻豆成人午夜福利视频| 青春草国产在线视频| 亚洲不卡免费看| 亚洲人成网站高清观看| 久久久欧美国产精品| 成人高潮视频无遮挡免费网站| 日本-黄色视频高清免费观看| 男女啪啪激烈高潮av片| 欧美日韩视频精品一区| 亚洲av.av天堂| 久久久精品欧美日韩精品| 69av精品久久久久久| videossex国产| 春色校园在线视频观看| 99久久精品国产国产毛片| 99热这里只有是精品在线观看| 极品少妇高潮喷水抽搐| 赤兔流量卡办理| 成人亚洲欧美一区二区av| 人妻制服诱惑在线中文字幕| 成年人午夜在线观看视频| 亚洲va在线va天堂va国产| 中文乱码字字幕精品一区二区三区| 亚洲av免费在线观看| 日本-黄色视频高清免费观看| 亚洲一区二区三区欧美精品 | 在线观看美女被高潮喷水网站| 国产午夜福利久久久久久| 日韩不卡一区二区三区视频在线| 麻豆国产97在线/欧美| 插阴视频在线观看视频| 欧美极品一区二区三区四区| 欧美日本视频| 九九在线视频观看精品| 国产男女超爽视频在线观看| 中文字幕免费在线视频6| 青春草国产在线视频| 小蜜桃在线观看免费完整版高清| 日韩欧美 国产精品| 亚洲av福利一区| 欧美xxxx黑人xx丫x性爽| 嘟嘟电影网在线观看| 中文资源天堂在线| 亚洲av国产av综合av卡| 美女主播在线视频| 日本与韩国留学比较| 18+在线观看网站| 精品酒店卫生间| 纵有疾风起免费观看全集完整版| 国产精品av视频在线免费观看| 国产亚洲最大av| 国模一区二区三区四区视频| 在线 av 中文字幕| 亚洲国产精品成人久久小说| 自拍偷自拍亚洲精品老妇| 日韩欧美精品v在线| 国产精品秋霞免费鲁丝片| av又黄又爽大尺度在线免费看| 亚洲国产欧美在线一区| 久久久精品欧美日韩精品| 99久国产av精品国产电影| 久久精品国产亚洲网站| 国产日韩欧美亚洲二区| 人妻 亚洲 视频| 免费av不卡在线播放| 国产精品成人在线| 亚洲精品乱久久久久久| av国产精品久久久久影院| 亚洲国产色片| 久久久久久久国产电影| a级一级毛片免费在线观看| 亚洲av二区三区四区| 国产亚洲5aaaaa淫片| 日日撸夜夜添| 欧美性感艳星| 在线免费十八禁| 涩涩av久久男人的天堂| 欧美高清成人免费视频www| 三级男女做爰猛烈吃奶摸视频| 日本午夜av视频| 亚洲一级一片aⅴ在线观看| 国产一区二区三区综合在线观看 | 纵有疾风起免费观看全集完整版| 国产日韩欧美亚洲二区| 青春草视频在线免费观看| 国产毛片a区久久久久| 寂寞人妻少妇视频99o| 免费av毛片视频| 观看美女的网站| 男人和女人高潮做爰伦理| 国产精品久久久久久久久免| 联通29元200g的流量卡| 最近中文字幕2019免费版| 自拍偷自拍亚洲精品老妇| 欧美日韩精品成人综合77777| 欧美国产精品一级二级三级 | 日韩欧美精品v在线| 日韩三级伦理在线观看| 久久韩国三级中文字幕| 精品亚洲乱码少妇综合久久| 久久女婷五月综合色啪小说 | 91久久精品国产一区二区三区| 日韩免费高清中文字幕av| 高清毛片免费看| 高清在线视频一区二区三区| 欧美成人a在线观看| 国国产精品蜜臀av免费| 成人无遮挡网站| av又黄又爽大尺度在线免费看| 少妇丰满av| 久久久久久国产a免费观看| 国产成人精品婷婷| 亚洲国产欧美在线一区| 欧美日韩综合久久久久久| 国产精品福利在线免费观看| 日本三级黄在线观看| 国产欧美日韩精品一区二区| av在线蜜桃| 高清日韩中文字幕在线| 国产亚洲午夜精品一区二区久久 | 国产片特级美女逼逼视频| 一区二区三区四区激情视频| 一级毛片 在线播放| 国产成人福利小说| 激情 狠狠 欧美| 搡老乐熟女国产| 欧美日韩视频高清一区二区三区二| 午夜激情久久久久久久| 亚洲欧美成人精品一区二区| 久久精品国产亚洲av天美| 精品一区二区三卡| 联通29元200g的流量卡| 久久久精品94久久精品| 视频中文字幕在线观看| 国精品久久久久久国模美| 男女边吃奶边做爰视频| 久久久精品免费免费高清| 久久精品国产亚洲av涩爱| 日韩av在线免费看完整版不卡| 久久久国产一区二区| 欧美成人一区二区免费高清观看| 91精品国产九色| 十八禁网站网址无遮挡 | 欧美成人a在线观看| 亚洲人成网站在线播| 看十八女毛片水多多多| 欧美日韩在线观看h| 国产av不卡久久| 午夜福利视频1000在线观看| 香蕉精品网在线| 一个人看视频在线观看www免费| 人妻夜夜爽99麻豆av| 亚洲精品色激情综合| 国产亚洲精品久久久com| 天堂网av新在线| 国产午夜福利久久久久久| 97超视频在线观看视频| 婷婷色综合www| 精品人妻一区二区三区麻豆| 夜夜看夜夜爽夜夜摸| 大片电影免费在线观看免费| 日本欧美国产在线视频| 亚洲人成网站高清观看| 国产免费福利视频在线观看| xxx大片免费视频| 在线 av 中文字幕| 精品视频人人做人人爽| 我的老师免费观看完整版| av又黄又爽大尺度在线免费看| 亚洲精品日本国产第一区| 男人添女人高潮全过程视频| 寂寞人妻少妇视频99o| 国产精品人妻久久久影院| 性色avwww在线观看| 一级爰片在线观看| 国产精品一区二区性色av| 国产av不卡久久| 免费在线观看成人毛片| 搞女人的毛片| 日韩伦理黄色片| 亚洲,一卡二卡三卡| 亚洲国产精品国产精品| 99视频精品全部免费 在线| 一级片'在线观看视频| 在线免费十八禁| 欧美日本视频| 看非洲黑人一级黄片| 少妇猛男粗大的猛烈进出视频 | 91aial.com中文字幕在线观看| 国产黄a三级三级三级人| 亚洲成人久久爱视频| 日韩欧美一区视频在线观看 | 成人无遮挡网站| 亚洲av在线观看美女高潮| 亚洲国产日韩一区二区| 久久热精品热| 日韩强制内射视频| 亚洲成人一二三区av| 亚洲最大成人中文| 午夜福利在线在线| 亚洲精品国产成人久久av| 久久久久九九精品影院| 国产精品国产三级国产专区5o| 亚洲无线观看免费| 国产精品福利在线免费观看| 国产国拍精品亚洲av在线观看| 嫩草影院精品99| 午夜爱爱视频在线播放| 亚洲欧美日韩另类电影网站 | 性色avwww在线观看| 新久久久久国产一级毛片| 国产白丝娇喘喷水9色精品| 国产成人91sexporn| 日本欧美国产在线视频| 观看免费一级毛片| 最近中文字幕2019免费版| 婷婷色综合大香蕉| 久久久久久久精品精品| 蜜桃亚洲精品一区二区三区| 人妻一区二区av| 精华霜和精华液先用哪个| 亚洲aⅴ乱码一区二区在线播放| 成年av动漫网址| 大陆偷拍与自拍| 在线观看三级黄色| 亚洲aⅴ乱码一区二区在线播放| 国产爽快片一区二区三区| 亚洲成色77777| 一级毛片 在线播放| 日韩精品有码人妻一区| 嫩草影院入口| 婷婷色综合www| 人妻 亚洲 视频| 亚洲精品成人久久久久久| 午夜免费鲁丝| 插阴视频在线观看视频| 国产伦精品一区二区三区四那| 边亲边吃奶的免费视频| 久久久国产一区二区| 日韩精品有码人妻一区| 嫩草影院新地址| 内地一区二区视频在线| 你懂的网址亚洲精品在线观看| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 别揉我奶头 嗯啊视频| 一级毛片黄色毛片免费观看视频| 一级片'在线观看视频| 少妇熟女欧美另类| 久久久久国产网址| 乱系列少妇在线播放| 秋霞伦理黄片| 亚洲欧美日韩东京热| 老女人水多毛片| 欧美三级亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品久久久久久精品电影小说 | 国产v大片淫在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 欧美高清成人免费视频www| 亚洲经典国产精华液单| 久久女婷五月综合色啪小说 | 国产精品久久久久久精品电影小说 | 特级一级黄色大片| 国产久久久一区二区三区| 久久久久国产网址| 亚洲精品乱久久久久久| 最近的中文字幕免费完整| 免费看光身美女| 国产伦精品一区二区三区视频9| 国产免费一区二区三区四区乱码| 日韩一区二区视频免费看| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| videos熟女内射| 精品熟女少妇av免费看| 激情五月婷婷亚洲| 亚洲久久久久久中文字幕| 麻豆久久精品国产亚洲av| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 日韩人妻高清精品专区| 嫩草影院新地址| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区av在线| 春色校园在线视频观看| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 好男人视频免费观看在线| 成人午夜精彩视频在线观看| 国产精品无大码| 国产在视频线精品| 亚洲,欧美,日韩| 免费av毛片视频| 欧美丝袜亚洲另类| 成人免费观看视频高清| 777米奇影视久久| 久久精品久久久久久久性| 免费av不卡在线播放| 国产午夜精品久久久久久一区二区三区| 国产男女内射视频| 免费观看的影片在线观看| 久久久久久久大尺度免费视频| 国产精品人妻久久久久久| 久久热精品热| 18禁裸乳无遮挡动漫免费视频 | 亚洲欧美中文字幕日韩二区| 狠狠精品人妻久久久久久综合| 联通29元200g的流量卡| 色综合色国产| 深夜a级毛片| 国产亚洲最大av| 久久久久国产精品人妻一区二区| 亚洲av成人精品一区久久| 国产黄片美女视频| 日本色播在线视频| 久久综合国产亚洲精品| 狂野欧美白嫩少妇大欣赏| 大又大粗又爽又黄少妇毛片口| 国产精品无大码| 97热精品久久久久久| 成年女人在线观看亚洲视频 | 99热这里只有是精品在线观看| 99精国产麻豆久久婷婷| 欧美性感艳星| 老女人水多毛片| 亚洲欧美精品专区久久| 亚洲最大成人手机在线| 日韩国内少妇激情av| 亚洲内射少妇av| 99久久人妻综合| 久久热精品热| videossex国产| 国产毛片在线视频| 国产成人免费观看mmmm| 日韩一区二区三区影片| 在线天堂最新版资源| 精品少妇久久久久久888优播| 精品一区二区三卡| 精品久久久久久电影网| 久久久久网色| 欧美激情久久久久久爽电影| 九色成人免费人妻av| 在线亚洲精品国产二区图片欧美 | 亚洲精品亚洲一区二区| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 精品久久久久久电影网| 99热这里只有是精品在线观看| 日本黄大片高清| 免费观看无遮挡的男女| 青春草视频在线免费观看| 亚洲伊人久久精品综合| 成年免费大片在线观看| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 免费av观看视频| 精品久久久久久久久av| 国产 一区 欧美 日韩| 国产精品成人在线| 国产av国产精品国产| 亚洲性久久影院| 亚洲欧洲日产国产| videos熟女内射| 视频区图区小说| 国产老妇伦熟女老妇高清| 久久女婷五月综合色啪小说 | 免费电影在线观看免费观看| 如何舔出高潮| 97精品久久久久久久久久精品| 看免费成人av毛片| 久久精品综合一区二区三区| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 日韩亚洲欧美综合| 大香蕉久久网| 一级毛片电影观看| 国产免费又黄又爽又色| 少妇人妻久久综合中文| 免费看日本二区| 你懂的网址亚洲精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产高清三级在线| 久久久精品免费免费高清| 亚洲欧美成人综合另类久久久| 一级av片app| 精品人妻偷拍中文字幕| 国产成人a∨麻豆精品| 日韩国内少妇激情av| 久久97久久精品| 久久久久国产网址| 天美传媒精品一区二区| 亚洲性久久影院| 欧美性猛交╳xxx乱大交人| 国产精品人妻久久久久久| 欧美精品国产亚洲| 如何舔出高潮| 少妇人妻一区二区三区视频| 超碰97精品在线观看| 久久99蜜桃精品久久| 国产中年淑女户外野战色| 麻豆久久精品国产亚洲av| 下体分泌物呈黄色| 我的老师免费观看完整版| 久久久久网色| av天堂中文字幕网| 一区二区三区免费毛片| 亚洲av不卡在线观看| 岛国毛片在线播放| 在线观看国产h片| 亚洲性久久影院| 国产高清国产精品国产三级 | 国产成人a∨麻豆精品| 色综合色国产| 91久久精品国产一区二区成人| 欧美日本视频| 亚洲内射少妇av| 一区二区三区免费毛片| 精品国产露脸久久av麻豆| 日韩av免费高清视频| 国产又色又爽无遮挡免| 深爱激情五月婷婷| 国产毛片a区久久久久| 成人无遮挡网站| 成年女人看的毛片在线观看| 亚洲欧美清纯卡通| av在线蜜桃| 国产 精品1| 免费看a级黄色片| 国产一级毛片在线| 在线观看一区二区三区| 熟妇人妻不卡中文字幕| 一级毛片 在线播放| 亚洲三级黄色毛片| 精品久久国产蜜桃| 一级毛片aaaaaa免费看小| 国产精品人妻久久久影院| 成人欧美大片| 久久综合国产亚洲精品| 岛国毛片在线播放| 一个人观看的视频www高清免费观看| 三级经典国产精品| 久久午夜福利片| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频 | xxx大片免费视频| 国产伦精品一区二区三区四那| 我要看日韩黄色一级片| 国产日韩欧美在线精品| 97在线视频观看| 18+在线观看网站| 青青草视频在线视频观看| 国产伦理片在线播放av一区| 免费大片18禁| 欧美激情久久久久久爽电影| 男男h啪啪无遮挡| 精品国产一区二区三区久久久樱花 | 天堂中文最新版在线下载 | 午夜免费男女啪啪视频观看| 亚洲av电影在线观看一区二区三区 | 国产精品偷伦视频观看了| 又爽又黄无遮挡网站| 少妇被粗大猛烈的视频| 熟妇人妻不卡中文字幕| 高清日韩中文字幕在线| 欧美97在线视频| 国产精品三级大全| 91精品一卡2卡3卡4卡| 日本av手机在线免费观看| 白带黄色成豆腐渣| 王馨瑶露胸无遮挡在线观看| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影小说 | 国产成人freesex在线| 国产成人a∨麻豆精品| 一级毛片黄色毛片免费观看视频| 好男人在线观看高清免费视频| 久久热精品热| 欧美3d第一页| 日本色播在线视频| 亚洲精品一区蜜桃| 97超视频在线观看视频| 青春草视频在线免费观看| 精品一区在线观看国产| 99re6热这里在线精品视频| 国产精品久久久久久av不卡| 日本免费在线观看一区| 在线 av 中文字幕| 男人舔奶头视频| 国产精品久久久久久久久免| 亚洲美女视频黄频| 成人美女网站在线观看视频| 丰满乱子伦码专区| 777米奇影视久久| 久久99热6这里只有精品| 我的老师免费观看完整版| 亚洲人与动物交配视频| 午夜福利在线在线| 成人亚洲精品av一区二区| 亚洲欧洲国产日韩| 久久久久久久久久久丰满| 久久人人爽人人爽人人片va| 秋霞伦理黄片| 欧美性感艳星| 黄色一级大片看看| 能在线免费看毛片的网站| 亚洲欧美精品自产自拍| 99热6这里只有精品| 亚洲自拍偷在线| 直男gayav资源| 国产成人精品婷婷| 亚洲性久久影院|