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

    基于DEM的大清河子流域劃分方法

    2020-10-12 14:29:16張曉嬌焦裕飛劉佳李旺林李傳哲
    人民黃河 2020年6期
    關(guān)鍵詞:大清河

    張曉嬌 焦裕飛 劉佳 李旺林 李傳哲

    摘?要:針對(duì)流域水循環(huán)模擬中的子流域劃分問(wèn)題,以大清河流域?yàn)槔紫然贏rcGIS軟件,以30 m分辨率的DEM作為數(shù)據(jù)源,進(jìn)行填洼處理、流向判定、匯流累計(jì)量計(jì)算等;然后采用河網(wǎng)密度法和水系分維法確定最佳集水面積閾值,提取出與實(shí)際最符合的河網(wǎng),以此更準(zhǔn)確地提取流域水系;最后根據(jù)水流流向、匯流累計(jì)以及出水口的位置劃分子流域,并基于MATLAB程序?qū)ψ恿饔蚣捌鋬?nèi)部柵格進(jìn)行編碼。結(jié)果表明:當(dāng)集水面積閾值為270 km2時(shí),提取的大清河流域河網(wǎng)水系與實(shí)際較為符合,且對(duì)子流域的劃分及編碼可以實(shí)現(xiàn)快速準(zhǔn)確定位每一個(gè)柵格,方便計(jì)算機(jī)的識(shí)別及演算。因此,以DEM為基礎(chǔ)實(shí)現(xiàn)大流域尺度下子流域劃分是可行的,對(duì)子流域的劃分及編碼可以為開(kāi)發(fā)構(gòu)建大尺度流域分布式水文模型提供技術(shù)支撐。

    關(guān)鍵詞:DEM;水系提取;子流域劃分;大清河

    中圖分類號(hào):TV213.4文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.06.003

    Study on Method of Sub-Basin Partition of Daqing River Based on DEM

    ZHANG Xiaojiao1,2,JIAO Yufei2,LIU Jia2,LI Wanglin1,LI Chuanzhe2

    (1. School of Water Conservancy and Environment,University of Jinan,Jinan 250022,China;

    2. State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin,China

    Institute of Water Resources and Hydropower Research,Beijing 100038,China)

    Abstract:Aiming at the sub-basin partition in the basin water cycle simulation, taking Daqing River basin as an example, based on ArcGIS software and taking the DEM with the resolution of 30 m as data source for depression filling, determination of flow direction and calculation of cumulative amount of confluence, this paper adopted the river network density method and water system fractal method to determine the optimal catchment area threshold and extract the most practical river network, in order to extract the drainage system of Daqing River basin more accurately. Finally, the molecular watershed was divided according to the flow direction, confluence accumulation and the location of the outlet, and MATLAB program was used for the sub-watershed coding. The results show that when the threshold of catchment area is 270 km2, the extracted river network system of Daqing River basin is more consistent with the actual situation, and the partition and coding of sub-basins can achieve rapid and accurate positioning to each grid, which is convenient for computer identification and calculation. Therefore, it is feasible to realize subcatchment of large watershed scale based on DEM, and the coding of subcatchments can also provide technical support for the development and construction of distributed hydrological model of large watershed scale.

    Key words: DEM;ditch extraction;subcatchments partition;Daqing River

    子流域劃分是構(gòu)建分布式水文模型的基礎(chǔ),分布式水文模型以子流域作為計(jì)算單元,可以清晰地了解單元之間的水文過(guò)程,并可很好地與傳統(tǒng)水文模型相結(jié)合,縮短模型開(kāi)發(fā)時(shí)間[1-2]。Jong-Gun K.等[3]研究認(rèn)為,子流域劃分對(duì)水文模擬結(jié)果的影響大于模型參數(shù)的影響,而且大尺度流域所劃分的子流域的數(shù)量不同使模擬精度存在不確定性。目前,基于數(shù)字高程模型(DEM)的流域地形特征提取方法已經(jīng)成為獲取流域信息的主要手段。OCallaghan等[4]基于水文學(xué)原理確定水流的路徑,并引入集水面積閾值的概念提取連續(xù)河網(wǎng)。傳統(tǒng)方法下作為子流域和河網(wǎng)劃分依據(jù)的集水面積閾值主要靠人為經(jīng)驗(yàn)設(shè)定,劃分結(jié)果與實(shí)際誤差較大。為了提高精度,學(xué)者們進(jìn)行了許多研究,目前確定最佳集水面積閾值的方法主要有河道平均坡降法[5]、河網(wǎng)密度法[6]、流域?qū)挾确植挤╗7]、水系分形法[8]等。 筆者針對(duì)大尺度流域的子流域劃分方法進(jìn)行研究,借助ArcGIS與MATLAB軟件,對(duì)大清河流域進(jìn)行河網(wǎng)提取、子流域劃分及編碼,以期為開(kāi)發(fā)構(gòu)建大清河流域分布式水文模型提供技術(shù)支撐。

    1?子流域劃分方法

    借助ArcGIS中的Hydrology工具箱以及MATLAB程序,對(duì)流域原始DEM進(jìn)行處理,步驟如下。

    (1)地形預(yù)處理。原始的DEM中大多存在洼地或尖峰,會(huì)導(dǎo)致提取的水系出現(xiàn)斷裂或者水流方向錯(cuò)誤的現(xiàn)象[9],因此需要對(duì)原始的DEM進(jìn)行填洼處理。

    (2)提取水流流向。D8算法[4]分別用數(shù)字 1、2、4、8、16、32、64、128 來(lái)表示東、東南、南、西南、西、西北、北、東北8個(gè)方向(見(jiàn)圖1),利用該算法生成流域水流流向。

    (3)計(jì)算匯流累計(jì)量。柵格匯流累計(jì)量的大小代表形成地表徑流的難易程度,數(shù)值越大,越易形成地表徑流。

    (4)提取河網(wǎng)水系。通過(guò)設(shè)定一個(gè)合理的集水面積閾值,根據(jù)水流流向及匯流累計(jì)量提取河網(wǎng)水系[10]。采用河網(wǎng)密度法及水系分維法兩種方法來(lái)確定最佳集水面積閾值:①河網(wǎng)密度法,選取不同的集水面積閾值對(duì)流域水系進(jìn)行提取,計(jì)算不同閾值下的河網(wǎng)密度及河源密度,分別繪制集水面積閾值與河網(wǎng)密度、河源密度的關(guān)系曲線,曲線趨于平穩(wěn)的點(diǎn)即為最佳集水面積閾值;②水系分維法,采用盒維數(shù)法[11],使用邊長(zhǎng)為r的正方形網(wǎng)格來(lái)覆蓋水系,確定覆蓋到水系的網(wǎng)格數(shù)N,分別取邊長(zhǎng)及網(wǎng)格數(shù)的對(duì)數(shù)進(jìn)行線性回歸,擬合直線的斜率即為水系的分維值,由此可以通過(guò)計(jì)算不同集水面積閾值下的水系分維值,與實(shí)際水系的分維值進(jìn)行比較,來(lái)確定最佳的集水面積閾值。

    (5)河網(wǎng)分級(jí)。采用Strahler分級(jí)法對(duì)流域提取的水系進(jìn)行分級(jí)[12],分級(jí)的原理是將河網(wǎng)中所有沒(méi)有支流匯入的河流作為第1級(jí)河流,兩個(gè)1級(jí)河流匯成第2級(jí)河流……依此類推,一直到河流出水口。在Strahler分級(jí)中,有且僅有同級(jí)別的兩條河流匯成一條河流時(shí),河流的級(jí)別才會(huì)增大。

    (6)確定子流域范圍。子流域的劃分是基于河流流向柵格和河流分級(jí)柵格,子流域內(nèi)匯流量最大的柵格為子流域的出水口,將所有能流向同一個(gè)子流域出水口的柵格劃在同一個(gè)子流域范圍內(nèi)[13]。

    (7)子流域編碼。通過(guò)一定的編碼規(guī)則確定柵格之間的拓?fù)潢P(guān)系,以此來(lái)反映柵格的計(jì)算順序,便于計(jì)算機(jī)識(shí)別[14-15]?;跂鸥裰g的水流流向以及匯流關(guān)系,對(duì)每個(gè)柵格賦予一個(gè)5層編碼:第一層編碼A表示子流域的編碼;第二層編碼B表示該柵格在子流域中的匯流層數(shù),子流域出水口柵格匯流層數(shù)為1,流入出水口柵格的柵格層數(shù)為2……依此類推;第三層編碼C表示匯流層內(nèi)的柵格順序,用來(lái)區(qū)分同一層內(nèi)的柵格;第四層編碼D表示流入柵格的匯流層;第五層編碼E用來(lái)區(qū)分河道與非河道柵格,其中河道柵格自出水口處從1開(kāi)始編碼,非河道柵格編碼為0。

    2?應(yīng)用實(shí)例

    2.1?研究區(qū)概況

    大清河流域位于海河流域中部,流域面積4.3萬(wàn)km2,地跨京、津、冀、晉4?。ㄊ校?。大清河發(fā)源于太行山,分為白溝河水系和趙王河水系。該流域?qū)儆诘湫偷臏貛Ъ撅L(fēng)性氣候區(qū),四季分明,降水量年際變化大、年內(nèi)分配不均,冬季和春季降雨少,夏季暴雨頻繁。

    2.2?河網(wǎng)提取

    借助ArcGIS軟件,對(duì)大清河流域原始DEM進(jìn)行填洼處理(見(jiàn)圖2),并利用D8算法提取水流流向(見(jiàn)圖3),計(jì)算匯流累計(jì)量。

    河網(wǎng)提取的重點(diǎn)是通過(guò)設(shè)定一個(gè)最佳的集水面積閾值,使得提取的河網(wǎng)與實(shí)際河網(wǎng)水系相吻合。對(duì)于同一個(gè)流域,選取不同的集水面積閾值將得到不同的河網(wǎng),閾值的選取也有很大的主觀性和隨意性,因此本文采用兩種方法確定最佳集水面積閾值。

    (1)河網(wǎng)密度法。分別選取9、90、225、270、315、360、405、450 km2作為集水面積閾值,對(duì)大清河流域進(jìn)行水系提取,統(tǒng)計(jì)不同閾值下的河流長(zhǎng)度、河源數(shù)量,計(jì)算出相應(yīng)的河網(wǎng)密度與河源密度,結(jié)果見(jiàn)表1。

    將集水面積閾值分別與河網(wǎng)密度、河源密度進(jìn)行趨勢(shì)線擬合(見(jiàn)圖4、圖5)。集水面積閾值與河網(wǎng)密度、河源密度都存在冪函數(shù)關(guān)系,擬合優(yōu)度R2都大于0.99。對(duì)所得到的冪函數(shù)求二階導(dǎo)數(shù)并再次進(jìn)行趨勢(shì)線擬合,將拐點(diǎn)作為最佳集水面積閾值(見(jiàn)圖6、圖7)。當(dāng)集水面積閾值為270 km2時(shí),河網(wǎng)密度與河源密度二階導(dǎo)數(shù)都趨于0且之后不再變化,因此確定大清河流域水系提取的最佳集水面積閾值為270 km2。

    (2)水系分形法。分別用邊長(zhǎng)r為100、500、1 000、2 000 m的網(wǎng)格對(duì)不同集水面積閾值下提取的水系進(jìn)行覆蓋,統(tǒng)計(jì)覆蓋水系的網(wǎng)格數(shù)目N,并用邊長(zhǎng)及網(wǎng)格數(shù)的對(duì)數(shù)之間的線性關(guān)系計(jì)算分維值,結(jié)果見(jiàn)表2。

    分別用邊長(zhǎng)為100、500、1 000、2 000 m的網(wǎng)格在ArcGIS中對(duì)實(shí)際水系矢量圖進(jìn)行覆蓋,得到覆蓋水系的網(wǎng)格數(shù),并用邊長(zhǎng)及網(wǎng)格數(shù)的對(duì)數(shù)之間的線性關(guān)系計(jì)算分維值,結(jié)果見(jiàn)表3。由表3可知,實(shí)際水系分維值為1.032,與集水面積閾值為270 km2時(shí)提取的水系分維值相同,因此確定最佳集水面積閾值為270 km2。

    綜合上述兩種方法,可確定對(duì)大清河流域進(jìn)行水系提取的最佳集水面積閾值為270 km2。

    2.3?河網(wǎng)分級(jí)

    大清河流域河流分級(jí)見(jiàn)圖8。當(dāng)最佳集水面積閾值為270 km2時(shí),大清河流域提取的河流一共有85條,其中:1級(jí)河流7條,占總河流數(shù)量的7.6%;2級(jí)河流10條,占總河流數(shù)量的10.7%;3級(jí)河流24條,占總河流數(shù)量的29.8%;4級(jí)河流44條,占總河流數(shù)量的51.9%。由此可見(jiàn),大清河流域水系密布、發(fā)育較好。

    2.4?子流域劃分

    基于DEM提取的大清河流域面積為43 209 km2,與實(shí)際收集的數(shù)據(jù)基本吻合,可以認(rèn)為該DEM精度能夠滿足提取流域特征的要求。集水面積閾值為270 km2時(shí),可將大清河流域劃分為56個(gè)子流域(見(jiàn)圖9),其中:子流域面積小于500 km2的占51%,在500~750 km2之間的占21%,在750~1 000 km2之間的占17%,大于1 000 km2的占11%。

    2.5?子流域編碼

    對(duì)子流域的編碼主要通過(guò)MATLAB程序來(lái)實(shí)現(xiàn)。圖10為大清河流域某一子流域內(nèi)所有柵格的編碼,以流域中某一柵格(18,3,5,3,150)為例,第一層編碼18表示第18個(gè)子流域,第二層編碼3為該柵格在第18個(gè)子流域中的第3層匯流層,第三層編碼5表示柵格為第3層匯流層中的第5個(gè)柵格,第四層編碼3表示流入該柵格的匯流所在層數(shù)為3,第五層編碼150表示第150條河道編碼。

    3?結(jié)?論

    本文采用河網(wǎng)密度法與水系分維法相互驗(yàn)證,確定大清河流域最佳集水面積閾值為270 km2,該條件下提取的水系與實(shí)際的水系吻合度很高。在此基礎(chǔ)上,應(yīng)用ArcGIS軟件對(duì)流域河網(wǎng)水系進(jìn)行分級(jí)并劃分子流域,然后借助MATLAB程序?qū)ψ恿饔蜻M(jìn)行編碼,可以實(shí)現(xiàn)對(duì)流域內(nèi)每一個(gè)柵格快速、精準(zhǔn)定位,有利于計(jì)算機(jī)的識(shí)別及演算。

    通過(guò)對(duì)大清河流域水文特征進(jìn)行提取與研究,發(fā)現(xiàn)大尺度流域的子流域劃分關(guān)鍵在于確定最佳的集水面積閾值。由DEM數(shù)據(jù)提取的河網(wǎng)與實(shí)際誤差較小,可以滿足水文模擬的需要。通過(guò)基于水流流向與匯流計(jì)算的子流域劃分及編碼,可以充分了解分布式水文模型構(gòu)建過(guò)程中柵格之間水量交換等水文過(guò)程,為開(kāi)發(fā)適合大尺度流域的分布式水文模型奠定基礎(chǔ)。

    參考文獻(xiàn):

    [1]?賈仰文,王浩,嚴(yán)登華.黑河流域水循環(huán)系統(tǒng)的分布式模擬(I):模型開(kāi)發(fā)與驗(yàn)證[J].水利學(xué)報(bào),2006,37(5):534-542.

    [2]?王中根,鄭紅星,劉昌明,等.黃河典型流域分布式水文模型及應(yīng)用研究[J].中國(guó)科學(xué)E輯:技術(shù)科學(xué),2004,34(增刊1):49-59.

    [3]?JONG-GUN K,YOUNSHIK P,DONGSUN Y,et al. Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds[J]. Journal of the American Water Resources Association,2009,45(4):963-972.

    [4]?OCALLAGHAN J F. The Extraction of Drainage Networks from Digital Elevation Data[J].Computer Vision Graphics and Image Processing,1984,28(3):323-344.

    [5]?熊立華,郭生練.基于DEM的數(shù)字河網(wǎng)生成方法的探討[J].長(zhǎng)江科學(xué)院院報(bào),2003,27(4):14-17.

    [6]?孔凡哲,李莉莉.利用DEM提取河網(wǎng)時(shí)集水面積閾值的確定[J].水電能源科學(xué),2005,23(4):65-67,93.

    [7]?楊邦,任立良.集水面積閾值確定方法的比較研究[J].水電能源科學(xué),2009,27(5):11-14.

    [8]?李麗. 分布式水文模型的匯流演算研究[D].南京:河海大學(xué),2007:32-53.

    [9]?徐新良,莊大方,賈紹鳳,等.GIS環(huán)境下基于DEM的中國(guó)流域自動(dòng)提取方法[J].長(zhǎng)江流域資源與環(huán)境,2004,13(4):343-348.

    [10]?宋向陽(yáng),吳發(fā)啟,趙龍山,等.基于DEM的延河流域水文特征提取與分析[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(4):200-206.

    [11]?SKUBALSKA-RAFAJOWICZ E. A New Method of Estimation of the Box-Counting Dimension of Multivariate Objects Using Space-Filling Curves[J]. Nonlinear Analysis: Theory,Methods & Applications,2005,63(5-7):e1281-e1287.

    [12]?MELTON M A. A Derivation of Strahlers Channel-Ordering System[J]. Journal of Geology,1959,67(3):345-346.

    [13]?葉愛(ài)中,夏軍,王綱勝,等.基于數(shù)字高程模型的河網(wǎng)提取及子流域生成[J].水利學(xué)報(bào),2005,36(5):531-537.

    [14]?劉佳嘉,周祖昊,賈仰文,等.分布式水文模型子流域編碼方法對(duì)比分析[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,45(1):22-29

    [15]?王春,江嶺,陳泰生,等.基于Pfafstetter規(guī)則的流域編碼算法并行化方法[J].地球信息科學(xué)學(xué)報(bào),2015,17(5):556-561.

    【責(zé)任編輯?張華興】

    猜你喜歡
    大清河
    大清河流域徑流量年際變化特征研究
    哦,大清河
    大清河流域綜合規(guī)劃獲水利部批復(fù)
    海河水利(2022年1期)2023-01-02 09:19:31
    給大清河一個(gè)美麗的承諾
    滹沱河、大清河(白洋淀)補(bǔ)水提前實(shí)現(xiàn)全線貫通目標(biāo)
    海河水利(2021年3期)2021-01-07 10:49:13
    基于單因子水質(zhì)標(biāo)識(shí)指數(shù)法的大清河流域府河段水質(zhì)評(píng)價(jià)
    綠色科技(2020年2期)2020-06-24 12:50:52
    大清河西海段河道風(fēng)場(chǎng)要素的確定
    大清河水環(huán)境質(zhì)量綜合評(píng)價(jià)
    【走遍山東】清河灘頭做釣翁
    走向世界(2017年3期)2017-06-30 15:38:24
    大清河中下游地下水水源評(píng)價(jià)區(qū)水量安全性評(píng)價(jià)
    欧美在线黄色| 国产精品一区二区三区四区久久| 国产真人三级小视频在线观看| 婷婷六月久久综合丁香| 亚洲18禁久久av| 日本 av在线| 亚洲18禁久久av| 国产亚洲精品久久久久5区| 久久亚洲真实| 国产精品久久视频播放| 亚洲人成77777在线视频| 老熟妇仑乱视频hdxx| 黑人操中国人逼视频| 91大片在线观看| 午夜亚洲福利在线播放| 日韩欧美精品v在线| 亚洲国产欧美网| 在线观看日韩欧美| 国产在线观看jvid| 欧美成人免费av一区二区三区| 听说在线观看完整版免费高清| 在线观看66精品国产| 精品久久久久久久久久久久久| 麻豆av在线久日| 免费av毛片视频| 女同久久另类99精品国产91| 亚洲精品av麻豆狂野| 又大又爽又粗| 午夜a级毛片| 亚洲最大成人中文| 成人精品一区二区免费| 制服人妻中文乱码| 观看免费一级毛片| 九色国产91popny在线| 淫妇啪啪啪对白视频| 99精品在免费线老司机午夜| 成人欧美大片| 国产探花在线观看一区二区| www日本在线高清视频| 神马国产精品三级电影在线观看 | 成熟少妇高潮喷水视频| 又黄又粗又硬又大视频| 午夜福利在线在线| 亚洲国产欧美人成| 好男人在线观看高清免费视频| 精品久久久久久久久久久久久| 激情在线观看视频在线高清| 亚洲av片天天在线观看| 成人av一区二区三区在线看| 麻豆国产97在线/欧美 | 欧美绝顶高潮抽搐喷水| 日本五十路高清| 黄色视频不卡| 久久久国产成人免费| 91大片在线观看| 久热爱精品视频在线9| 欧美在线黄色| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 99国产精品99久久久久| 亚洲成人久久性| 美女黄网站色视频| 桃红色精品国产亚洲av| 中文字幕最新亚洲高清| 又紧又爽又黄一区二区| 在线视频色国产色| 精品国产乱码久久久久久男人| 午夜老司机福利片| 麻豆一二三区av精品| 亚洲成人久久性| 亚洲熟妇中文字幕五十中出| 操出白浆在线播放| 色av中文字幕| 最近最新中文字幕大全电影3| 最新美女视频免费是黄的| 一级毛片女人18水好多| 国产亚洲精品综合一区在线观看 | 淫秽高清视频在线观看| 欧美zozozo另类| 国产精品久久视频播放| 制服诱惑二区| 婷婷丁香在线五月| 亚洲av成人精品一区久久| 丁香六月欧美| 欧美 亚洲 国产 日韩一| 免费无遮挡裸体视频| 亚洲精品一卡2卡三卡4卡5卡| 久久精品91蜜桃| 欧美又色又爽又黄视频| 日韩欧美一区二区三区在线观看| 露出奶头的视频| 亚洲精品中文字幕一二三四区| 这个男人来自地球电影免费观看| 精品久久久久久久久久免费视频| 国产男靠女视频免费网站| 免费看美女性在线毛片视频| 成人亚洲精品av一区二区| 老司机在亚洲福利影院| 黄色女人牲交| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看 | www.熟女人妻精品国产| 小说图片视频综合网站| 日本精品一区二区三区蜜桃| 中文字幕精品亚洲无线码一区| 看黄色毛片网站| 国产三级中文精品| 免费看美女性在线毛片视频| 黄色丝袜av网址大全| 黑人巨大精品欧美一区二区mp4| 国产精品一及| 美女 人体艺术 gogo| 女人爽到高潮嗷嗷叫在线视频| 欧美日本亚洲视频在线播放| 国产黄色小视频在线观看| 精品欧美一区二区三区在线| 亚洲欧美精品综合一区二区三区| 18禁国产床啪视频网站| 国产精华一区二区三区| 日日摸夜夜添夜夜添小说| 一进一出抽搐gif免费好疼| 亚洲精品国产精品久久久不卡| 三级毛片av免费| 色综合婷婷激情| 色在线成人网| 日日干狠狠操夜夜爽| 蜜桃久久精品国产亚洲av| 国产亚洲av高清不卡| www.999成人在线观看| 最近在线观看免费完整版| 久久 成人 亚洲| 97碰自拍视频| cao死你这个sao货| 淫妇啪啪啪对白视频| 午夜福利高清视频| 香蕉丝袜av| 成熟少妇高潮喷水视频| 国产精品久久久久久久电影 | 在线观看一区二区三区| 久久久久九九精品影院| 天堂av国产一区二区熟女人妻 | 别揉我奶头~嗯~啊~动态视频| 国产激情偷乱视频一区二区| 国产精品久久久av美女十八| 久久九九热精品免费| 欧美最黄视频在线播放免费| 18禁黄网站禁片免费观看直播| 一级片免费观看大全| 中文字幕最新亚洲高清| 日韩精品青青久久久久久| 99国产综合亚洲精品| 色av中文字幕| 欧美一区二区精品小视频在线| 俄罗斯特黄特色一大片| 国产精品 欧美亚洲| 全区人妻精品视频| 九色成人免费人妻av| 日本撒尿小便嘘嘘汇集6| 在线十欧美十亚洲十日本专区| 欧美黑人精品巨大| 亚洲人成电影免费在线| 欧美乱色亚洲激情| 免费av毛片视频| 国产欧美日韩一区二区精品| 两性午夜刺激爽爽歪歪视频在线观看 | 美女黄网站色视频| 亚洲欧美激情综合另类| 变态另类成人亚洲欧美熟女| 国产一区二区三区在线臀色熟女| 色尼玛亚洲综合影院| 麻豆久久精品国产亚洲av| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 欧美三级亚洲精品| 黄色视频,在线免费观看| 九色国产91popny在线| 熟女电影av网| 香蕉丝袜av| 亚洲第一电影网av| 男人舔女人的私密视频| 亚洲精品中文字幕一二三四区| 中文字幕最新亚洲高清| av天堂在线播放| 少妇裸体淫交视频免费看高清 | 亚洲美女黄片视频| 色播亚洲综合网| 日韩大尺度精品在线看网址| 人人妻人人澡欧美一区二区| 一本一本综合久久| 国产精品亚洲美女久久久| 国产主播在线观看一区二区| 美女 人体艺术 gogo| 欧美日韩国产亚洲二区| 他把我摸到了高潮在线观看| 日本撒尿小便嘘嘘汇集6| 激情在线观看视频在线高清| a级毛片a级免费在线| 欧美日韩亚洲国产一区二区在线观看| 久久久久久人人人人人| 国产成人精品久久二区二区91| 观看免费一级毛片| 亚洲精品一区av在线观看| 99久久无色码亚洲精品果冻| 久久精品国产综合久久久| 久久精品夜夜夜夜夜久久蜜豆 | 午夜激情av网站| 久久久水蜜桃国产精品网| 欧美丝袜亚洲另类 | 搡老岳熟女国产| 午夜福利视频1000在线观看| 一本大道久久a久久精品| 国产主播在线观看一区二区| 老司机午夜十八禁免费视频| 精品国产超薄肉色丝袜足j| 欧美日韩中文字幕国产精品一区二区三区| 高清在线国产一区| 男人舔奶头视频| 亚洲欧美一区二区三区黑人| 成人一区二区视频在线观看| 国产主播在线观看一区二区| 国产区一区二久久| 亚洲电影在线观看av| 欧美性猛交╳xxx乱大交人| 在线播放国产精品三级| av天堂在线播放| 日韩高清综合在线| 欧美高清成人免费视频www| 男女做爰动态图高潮gif福利片| 麻豆成人午夜福利视频| 激情在线观看视频在线高清| 国产精品精品国产色婷婷| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 午夜福利视频1000在线观看| 国产精品99久久99久久久不卡| 欧美日韩中文字幕国产精品一区二区三区| 久久热在线av| 亚洲美女视频黄频| xxxwww97欧美| 天天躁夜夜躁狠狠躁躁| 午夜免费激情av| 亚洲专区字幕在线| 老汉色av国产亚洲站长工具| 亚洲天堂国产精品一区在线| 精品久久久久久,| 亚洲欧美日韩高清专用| 欧美精品亚洲一区二区| 欧美日韩乱码在线| 一级毛片高清免费大全| 国产私拍福利视频在线观看| www.精华液| 嫩草影视91久久| 麻豆av在线久日| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 成人三级做爰电影| 天天躁夜夜躁狠狠躁躁| 18美女黄网站色大片免费观看| 亚洲国产欧美网| 国产精品免费一区二区三区在线| 熟女少妇亚洲综合色aaa.| 午夜福利高清视频| 久久久久九九精品影院| 白带黄色成豆腐渣| 久久中文字幕人妻熟女| 亚洲精品av麻豆狂野| 婷婷六月久久综合丁香| 91麻豆av在线| 女人高潮潮喷娇喘18禁视频| 欧美日韩一级在线毛片| 777久久人妻少妇嫩草av网站| 欧美性猛交黑人性爽| 最新美女视频免费是黄的| 美女免费视频网站| 桃色一区二区三区在线观看| 亚洲专区中文字幕在线| 在线观看午夜福利视频| 欧美性长视频在线观看| 嫩草影院精品99| 久久人妻av系列| 啦啦啦免费观看视频1| 久久久久九九精品影院| 波多野结衣巨乳人妻| 波多野结衣高清无吗| 久久中文字幕人妻熟女| 色尼玛亚洲综合影院| 亚洲熟女毛片儿| 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 亚洲 欧美一区二区三区| 日韩免费av在线播放| 色综合站精品国产| 精品日产1卡2卡| 一二三四在线观看免费中文在| 丰满的人妻完整版| 精品欧美国产一区二区三| 国产激情欧美一区二区| 午夜福利18| 人成视频在线观看免费观看| 淫妇啪啪啪对白视频| 大型黄色视频在线免费观看| 国产成人精品久久二区二区91| 国产高清videossex| 中文亚洲av片在线观看爽| 久久久水蜜桃国产精品网| 免费在线观看亚洲国产| 亚洲欧美精品综合一区二区三区| 久久精品综合一区二区三区| 色尼玛亚洲综合影院| 悠悠久久av| 免费电影在线观看免费观看| svipshipincom国产片| 婷婷精品国产亚洲av| 麻豆国产97在线/欧美 | 成人欧美大片| 波多野结衣高清无吗| 国产高清有码在线观看视频 | 欧美乱码精品一区二区三区| 成人国产一区最新在线观看| 久久中文字幕一级| 色尼玛亚洲综合影院| www.自偷自拍.com| 久久精品91无色码中文字幕| 国产亚洲精品av在线| 国产一区二区在线观看日韩 | 我的老师免费观看完整版| 欧美性猛交黑人性爽| 国产欧美日韩一区二区精品| 免费无遮挡裸体视频| 国产一区二区激情短视频| 欧美黄色淫秽网站| 成人特级黄色片久久久久久久| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| www.999成人在线观看| 老司机福利观看| 日韩欧美 国产精品| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 亚洲欧美精品综合一区二区三区| 亚洲成人国产一区在线观看| 91国产中文字幕| 欧美另类亚洲清纯唯美| 国产精品亚洲一级av第二区| 欧美黄色片欧美黄色片| 国产熟女午夜一区二区三区| 黄色a级毛片大全视频| 神马国产精品三级电影在线观看 | 国产99久久九九免费精品| 成人国产一区最新在线观看| 亚洲性夜色夜夜综合| 欧美成人免费av一区二区三区| 亚洲一区二区三区不卡视频| 国产蜜桃级精品一区二区三区| 国产成人精品久久二区二区91| 男女那种视频在线观看| 男女午夜视频在线观看| 成人国语在线视频| 国产在线观看jvid| 亚洲av电影不卡..在线观看| 日本 av在线| 三级毛片av免费| 国语自产精品视频在线第100页| 性色av乱码一区二区三区2| 欧美高清成人免费视频www| 久久久国产精品麻豆| 精品久久久久久久人妻蜜臀av| 国产私拍福利视频在线观看| 久久草成人影院| 亚洲中文字幕一区二区三区有码在线看 | 少妇被粗大的猛进出69影院| www.精华液| 久久精品aⅴ一区二区三区四区| 国产激情久久老熟女| 午夜福利欧美成人| 色尼玛亚洲综合影院| 黄色女人牲交| 黑人欧美特级aaaaaa片| 国产免费av片在线观看野外av| 国内久久婷婷六月综合欲色啪| 国产蜜桃级精品一区二区三区| 高清在线国产一区| 久久久久久久久久黄片| 日韩大尺度精品在线看网址| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品av在线| av在线天堂中文字幕| 一本综合久久免费| 久久久久久久精品吃奶| 久久国产精品影院| 俺也久久电影网| 在线观看www视频免费| 精品欧美国产一区二区三| 老司机靠b影院| 1024手机看黄色片| 日韩欧美国产一区二区入口| 桃色一区二区三区在线观看| 亚洲国产看品久久| 国产精品 欧美亚洲| 久久天堂一区二区三区四区| 午夜精品一区二区三区免费看| 久久久久久久久久黄片| 免费一级毛片在线播放高清视频| 久久中文字幕人妻熟女| 欧美激情久久久久久爽电影| 1024视频免费在线观看| 亚洲精品中文字幕在线视频| 午夜老司机福利片| 少妇人妻一区二区三区视频| 欧美极品一区二区三区四区| 国产欧美日韩一区二区三| 成人18禁高潮啪啪吃奶动态图| 日本精品一区二区三区蜜桃| 精品高清国产在线一区| 制服丝袜大香蕉在线| 久久久久国产精品人妻aⅴ院| 亚洲黑人精品在线| 精品少妇一区二区三区视频日本电影| 在线观看免费日韩欧美大片| 精品电影一区二区在线| 欧美丝袜亚洲另类 | 国产伦人伦偷精品视频| 99久久99久久久精品蜜桃| 久久久水蜜桃国产精品网| 国产aⅴ精品一区二区三区波| 国产一区二区在线观看日韩 | 一级黄色大片毛片| 欧美黑人巨大hd| 成人高潮视频无遮挡免费网站| 国产黄a三级三级三级人| 50天的宝宝边吃奶边哭怎么回事| 色综合亚洲欧美另类图片| 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 在线a可以看的网站| 亚洲七黄色美女视频| 人人妻人人看人人澡| 亚洲男人天堂网一区| 夜夜夜夜夜久久久久| 一二三四在线观看免费中文在| 久久中文看片网| 首页视频小说图片口味搜索| www日本在线高清视频| 日本精品一区二区三区蜜桃| 非洲黑人性xxxx精品又粗又长| 亚洲熟女毛片儿| 国模一区二区三区四区视频 | 欧美午夜高清在线| 免费在线观看日本一区| 国产一区二区三区在线臀色熟女| 亚洲国产精品合色在线| www.www免费av| 少妇粗大呻吟视频| 又爽又黄无遮挡网站| 在线观看一区二区三区| 亚洲在线自拍视频| 舔av片在线| 欧美 亚洲 国产 日韩一| 中亚洲国语对白在线视频| 精品国产超薄肉色丝袜足j| 可以在线观看的亚洲视频| 久久久久亚洲av毛片大全| 老司机福利观看| 国产精品免费一区二区三区在线| 波多野结衣巨乳人妻| 香蕉丝袜av| 黄色成人免费大全| 午夜福利视频1000在线观看| 天天躁夜夜躁狠狠躁躁| 国产区一区二久久| 国产精品久久久久久精品电影| 国产在线精品亚洲第一网站| 舔av片在线| 老熟妇乱子伦视频在线观看| 毛片女人毛片| 国产三级在线视频| 成人18禁在线播放| 哪里可以看免费的av片| 99国产精品99久久久久| 18禁美女被吸乳视频| 在线观看免费视频日本深夜| 国产成人av教育| 一进一出抽搐gif免费好疼| 欧美色欧美亚洲另类二区| 欧美黑人精品巨大| 亚洲精品av麻豆狂野| 性色av乱码一区二区三区2| 国产精品1区2区在线观看.| 亚洲成av人片免费观看| 99国产极品粉嫩在线观看| 欧美中文综合在线视频| 亚洲精品美女久久av网站| 国产成人精品久久二区二区免费| 极品教师在线免费播放| 特级一级黄色大片| 琪琪午夜伦伦电影理论片6080| 国产精品野战在线观看| 久久精品综合一区二区三区| 99国产精品99久久久久| 草草在线视频免费看| 国产激情久久老熟女| 日本 av在线| 国产aⅴ精品一区二区三区波| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 国产91精品成人一区二区三区| 最新美女视频免费是黄的| 国产成年人精品一区二区| 亚洲人成网站高清观看| 国产精品,欧美在线| www.精华液| avwww免费| 人妻久久中文字幕网| 午夜福利欧美成人| 午夜福利成人在线免费观看| 少妇人妻一区二区三区视频| 国产精品亚洲av一区麻豆| 午夜影院日韩av| 97超级碰碰碰精品色视频在线观看| av天堂在线播放| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 亚洲人成伊人成综合网2020| 99国产极品粉嫩在线观看| 99久久精品国产亚洲精品| 男女午夜视频在线观看| 国产区一区二久久| 亚洲自拍偷在线| 精品欧美一区二区三区在线| 97人妻精品一区二区三区麻豆| 精品久久蜜臀av无| 天天躁狠狠躁夜夜躁狠狠躁| xxxwww97欧美| 国产精品久久久久久精品电影| 一区二区三区国产精品乱码| 久久精品91无色码中文字幕| 午夜福利在线在线| 天堂影院成人在线观看| 俄罗斯特黄特色一大片| 亚洲电影在线观看av| 国产97色在线日韩免费| 日韩欧美三级三区| 欧美乱色亚洲激情| 制服诱惑二区| 99re在线观看精品视频| 99精品在免费线老司机午夜| 免费人成视频x8x8入口观看| 三级男女做爰猛烈吃奶摸视频| 十八禁人妻一区二区| 亚洲中文字幕日韩| 一级作爱视频免费观看| 国产精品久久久久久久电影 | 三级国产精品欧美在线观看 | 97超级碰碰碰精品色视频在线观看| 欧美日本亚洲视频在线播放| 亚洲第一欧美日韩一区二区三区| www日本黄色视频网| 国产单亲对白刺激| 美女午夜性视频免费| 黄色丝袜av网址大全| 免费在线观看完整版高清| 中文字幕久久专区| 久久精品综合一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 亚洲七黄色美女视频| 国产黄片美女视频| 午夜免费激情av| 日韩欧美 国产精品| x7x7x7水蜜桃| 欧美在线黄色| 国产精品野战在线观看| 国产精品久久视频播放| 国产亚洲精品第一综合不卡| 91大片在线观看| 精品国产乱码久久久久久男人| 不卡一级毛片| 人妻久久中文字幕网| 久久久久久国产a免费观看| 午夜激情av网站| 免费在线观看亚洲国产| 久久久久久久久久黄片| 在线十欧美十亚洲十日本专区| 在线观看免费午夜福利视频| 亚洲精品久久成人aⅴ小说| 人人妻人人看人人澡| 男人的好看免费观看在线视频 | 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 亚洲一区高清亚洲精品| www日本黄色视频网| svipshipincom国产片| 精品乱码久久久久久99久播| 日本一本二区三区精品| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 亚洲精品美女久久av网站| 窝窝影院91人妻| videosex国产| 国产免费男女视频| 免费搜索国产男女视频| 亚洲精品在线美女| 国产真实乱freesex| 亚洲精品在线观看二区| 亚洲欧美精品综合久久99| 在线永久观看黄色视频| 亚洲av成人精品一区久久| 国产一区二区三区在线臀色熟女| 不卡av一区二区三区| 久久精品成人免费网站| 精品国产超薄肉色丝袜足j| 亚洲人成伊人成综合网2020| av国产免费在线观看| 日本一本二区三区精品|