• <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à)
    国内毛片毛片毛片毛片毛片| 最新的欧美精品一区二区| 精品国产一区二区三区久久久樱花| 欧美亚洲 丝袜 人妻 在线| 波多野结衣一区麻豆| 午夜激情久久久久久久| av电影中文网址| 男人舔女人的私密视频| 亚洲精品国产精品久久久不卡| 交换朋友夫妻互换小说| 少妇的丰满在线观看| 黄色毛片三级朝国网站| 91字幕亚洲| 高清视频免费观看一区二区| a在线观看视频网站| 免费不卡黄色视频| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 国产精品二区激情视频| 中文欧美无线码| 免费av中文字幕在线| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区三区在线| 女人久久www免费人成看片| 亚洲精华国产精华精| 后天国语完整版免费观看| 国产男靠女视频免费网站| 国产99久久九九免费精品| 老司机在亚洲福利影院| 伊人久久大香线蕉亚洲五| 天天躁夜夜躁狠狠躁躁| 老汉色∧v一级毛片| 久久国产精品男人的天堂亚洲| 乱人伦中国视频| 亚洲欧美色中文字幕在线| 午夜福利,免费看| 成年人免费黄色播放视频| 母亲3免费完整高清在线观看| 美女主播在线视频| 美女扒开内裤让男人捅视频| 18禁裸乳无遮挡动漫免费视频| 亚洲avbb在线观看| 啦啦啦在线免费观看视频4| 国产精品电影一区二区三区 | 中国美女看黄片| 色94色欧美一区二区| 亚洲熟妇熟女久久| 午夜福利,免费看| 久久精品国产亚洲av香蕉五月 | 男人舔女人的私密视频| 午夜福利影视在线免费观看| 热re99久久国产66热| 日韩免费高清中文字幕av| 一区二区三区国产精品乱码| 国产精品影院久久| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 五月天丁香电影| 丁香欧美五月| 狠狠婷婷综合久久久久久88av| 日本a在线网址| 少妇粗大呻吟视频| 9色porny在线观看| 在线观看免费日韩欧美大片| 777米奇影视久久| 精品久久久久久电影网| 熟女少妇亚洲综合色aaa.| 看免费av毛片| 丰满少妇做爰视频| 美国免费a级毛片| 精品第一国产精品| 日韩一卡2卡3卡4卡2021年| 一区二区av电影网| 日本一区二区免费在线视频| 黄色片一级片一级黄色片| 在线观看免费日韩欧美大片| 欧美乱妇无乱码| 国产精品偷伦视频观看了| 亚洲人成电影免费在线| 男人操女人黄网站| 亚洲国产av新网站| 国产麻豆69| 1024视频免费在线观看| 男人操女人黄网站| 人人澡人人妻人| 咕卡用的链子| 777久久人妻少妇嫩草av网站| 国产成人啪精品午夜网站| 久久性视频一级片| 日韩人妻精品一区2区三区| 国产极品粉嫩免费观看在线| 成人黄色视频免费在线看| 欧美大码av| 91精品国产国语对白视频| 久久国产精品影院| 久久国产精品男人的天堂亚洲| 久久久国产精品麻豆| 免费观看人在逋| 久久精品成人免费网站| 久久久国产精品麻豆| 性高湖久久久久久久久免费观看| 俄罗斯特黄特色一大片| 亚洲久久久国产精品| 精品一区二区三卡| 久久亚洲精品不卡| avwww免费| 亚洲 国产 在线| 亚洲美女黄片视频| 搡老岳熟女国产| 国产淫语在线视频| 久久久久网色| 大码成人一级视频| 美女扒开内裤让男人捅视频| 国产一区二区在线观看av| 国产av又大| 亚洲天堂av无毛| 亚洲avbb在线观看| 久久久久国产一级毛片高清牌| 国产成人一区二区三区免费视频网站| 久9热在线精品视频| 亚洲人成电影免费在线| 桃花免费在线播放| 少妇的丰满在线观看| 国产亚洲欧美精品永久| 久久ye,这里只有精品| 中文字幕最新亚洲高清| 国产精品 欧美亚洲| 欧美国产精品一级二级三级| 免费一级毛片在线播放高清视频 | 日本av手机在线免费观看| 757午夜福利合集在线观看| 精品一区二区三卡| 久久精品亚洲av国产电影网| 精品高清国产在线一区| 国产精品久久久人人做人人爽| 成人手机av| 免费一级毛片在线播放高清视频 | 中文字幕最新亚洲高清| 婷婷成人精品国产| 99精国产麻豆久久婷婷| 久久午夜综合久久蜜桃| 成人国产一区最新在线观看| 男女下面插进去视频免费观看| 一级毛片电影观看| 亚洲少妇的诱惑av| 又大又爽又粗| 亚洲欧美一区二区三区久久| 18禁国产床啪视频网站| 老司机午夜十八禁免费视频| 在线观看免费视频日本深夜| 日韩视频一区二区在线观看| 欧美久久黑人一区二区| 亚洲国产毛片av蜜桃av| 久9热在线精品视频| 精品久久久精品久久久| 亚洲男人天堂网一区| 69精品国产乱码久久久| 色老头精品视频在线观看| 超碰97精品在线观看| 91成年电影在线观看| 国产97色在线日韩免费| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 伦理电影免费视频| 三上悠亚av全集在线观看| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 日韩制服丝袜自拍偷拍| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区久久| 99精品欧美一区二区三区四区| 日本wwww免费看| 激情在线观看视频在线高清 | 欧美成狂野欧美在线观看| 热re99久久国产66热| av在线播放免费不卡| 两性夫妻黄色片| 久久精品国产综合久久久| 男女无遮挡免费网站观看| 中文欧美无线码| 咕卡用的链子| 亚洲久久久国产精品| 国产区一区二久久| 精品国产超薄肉色丝袜足j| 亚洲av日韩精品久久久久久密| 老司机午夜福利在线观看视频 | 久久99一区二区三区| 久久午夜综合久久蜜桃| 久久国产精品男人的天堂亚洲| 国产精品偷伦视频观看了| 午夜日韩欧美国产| 日韩精品免费视频一区二区三区| 中文字幕人妻熟女乱码| 在线 av 中文字幕| 国产精品一区二区在线观看99| 一级片'在线观看视频| 老汉色∧v一级毛片| 少妇裸体淫交视频免费看高清 | 久久精品成人免费网站| 亚洲专区字幕在线| 日本精品一区二区三区蜜桃| av网站在线播放免费| 动漫黄色视频在线观看| 国产成人精品无人区| 91老司机精品| 亚洲国产看品久久| 国产成人一区二区三区免费视频网站| 动漫黄色视频在线观看| 亚洲七黄色美女视频| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| a在线观看视频网站| 咕卡用的链子| 国产一区二区激情短视频| 久久久久国产一级毛片高清牌| 精品国产超薄肉色丝袜足j| 黄频高清免费视频| 91精品国产国语对白视频| 久久精品aⅴ一区二区三区四区| 国产精品美女特级片免费视频播放器 | 亚洲欧美激情在线| 男女免费视频国产| 丝袜美足系列| 久久毛片免费看一区二区三区| 精品国产一区二区久久| 动漫黄色视频在线观看| 一级,二级,三级黄色视频| 国产成人免费无遮挡视频| 亚洲国产看品久久| 亚洲免费av在线视频| 女同久久另类99精品国产91| 午夜福利免费观看在线| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 免费观看人在逋| 免费不卡黄色视频| 欧美黑人欧美精品刺激| 亚洲成人国产一区在线观看| 精品久久久久久久毛片微露脸| 狠狠婷婷综合久久久久久88av| 精品一区二区三区视频在线观看免费 | av天堂在线播放| 久久中文看片网| 成人手机av| 亚洲av片天天在线观看| 亚洲精品一卡2卡三卡4卡5卡| 91成人精品电影| 免费在线观看完整版高清| 在线观看免费日韩欧美大片| 国产激情久久老熟女| 午夜两性在线视频| 91成年电影在线观看| 亚洲第一欧美日韩一区二区三区 | 热99久久久久精品小说推荐| 日韩大码丰满熟妇| 免费不卡黄色视频| 精品国产国语对白av| 亚洲伊人色综图| 女性生殖器流出的白浆| 欧美日韩视频精品一区| 欧美精品av麻豆av| 天堂俺去俺来也www色官网| 日韩 欧美 亚洲 中文字幕| 成人亚洲精品一区在线观看| 精品熟女少妇八av免费久了| 国产高清videossex| 一级a爱视频在线免费观看| 日韩 欧美 亚洲 中文字幕| 久久人妻熟女aⅴ| 成年女人毛片免费观看观看9 | 国产在线精品亚洲第一网站| 丁香欧美五月| 首页视频小说图片口味搜索| 极品少妇高潮喷水抽搐| 一区二区av电影网| 大型黄色视频在线免费观看| 波多野结衣一区麻豆| 亚洲专区国产一区二区| 高清黄色对白视频在线免费看| 国产成人精品久久二区二区免费| 国产欧美亚洲国产| 久久青草综合色| 欧美午夜高清在线| 精品欧美一区二区三区在线| 亚洲伊人色综图| 一区在线观看完整版| 老司机在亚洲福利影院| 天天操日日干夜夜撸| 香蕉久久夜色| 国产一区二区三区在线臀色熟女 | 国产男女内射视频| 激情视频va一区二区三区| 日本a在线网址| 亚洲av片天天在线观看| 色94色欧美一区二区| 高潮久久久久久久久久久不卡| av欧美777| av片东京热男人的天堂| 自线自在国产av| 高清毛片免费观看视频网站 | 亚洲成a人片在线一区二区| 欧美日韩黄片免| 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 国产亚洲精品久久久久5区| 精品高清国产在线一区| 久久久久精品人妻al黑| 看免费av毛片| 啦啦啦 在线观看视频| 最近最新中文字幕大全电影3 | 精品少妇黑人巨大在线播放| 国产精品一区二区免费欧美| 午夜福利,免费看| 高清在线国产一区| 久久久精品免费免费高清| 精品欧美一区二区三区在线| 国产又爽黄色视频| 91麻豆av在线| 亚洲欧美色中文字幕在线| 日本撒尿小便嘘嘘汇集6| 大型黄色视频在线免费观看| 无人区码免费观看不卡 | av又黄又爽大尺度在线免费看| 99国产极品粉嫩在线观看| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 九色亚洲精品在线播放| 麻豆av在线久日| 欧美日韩av久久| 久久精品国产亚洲av高清一级| 搡老熟女国产l中国老女人| 岛国毛片在线播放| 亚洲色图综合在线观看| 国产三级黄色录像| 精品国内亚洲2022精品成人 | 日韩欧美免费精品| 在线观看舔阴道视频| 大香蕉久久网| av一本久久久久| 一区二区日韩欧美中文字幕| 免费久久久久久久精品成人欧美视频| 日韩欧美一区视频在线观看| av一本久久久久| 久久精品国产亚洲av香蕉五月 | 亚洲专区国产一区二区| 又大又爽又粗| 丰满饥渴人妻一区二区三| 又大又爽又粗| 精品一区二区三区视频在线观看免费 | 999精品在线视频| 久久性视频一级片| 亚洲国产欧美日韩在线播放| 女警被强在线播放| 韩国精品一区二区三区| 三级毛片av免费| 成人av一区二区三区在线看| 一进一出抽搐动态| 三上悠亚av全集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜视频精品福利| 热99re8久久精品国产| 天天躁夜夜躁狠狠躁躁| 99久久精品国产亚洲精品| 好男人电影高清在线观看| 欧美亚洲 丝袜 人妻 在线| av免费在线观看网站| 啪啪无遮挡十八禁网站| 亚洲欧美精品综合一区二区三区| 最新的欧美精品一区二区| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 亚洲精品在线美女| 亚洲天堂av无毛| 飞空精品影院首页| 日韩有码中文字幕| 亚洲少妇的诱惑av| 777米奇影视久久| 日韩 欧美 亚洲 中文字幕| 一级片'在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 中国美女看黄片| 亚洲欧洲精品一区二区精品久久久| 国产成人影院久久av| 18禁裸乳无遮挡动漫免费视频| 亚洲午夜精品一区,二区,三区| 大香蕉久久网| 国产亚洲午夜精品一区二区久久| 成人av一区二区三区在线看| 在线十欧美十亚洲十日本专区| 九色亚洲精品在线播放| 免费在线观看完整版高清| 国产精品秋霞免费鲁丝片| 麻豆乱淫一区二区| 性色av乱码一区二区三区2| 久久久久久久久免费视频了| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 亚洲中文字幕日韩| 99国产综合亚洲精品| 日本av手机在线免费观看| 两个人免费观看高清视频| 免费不卡黄色视频| 飞空精品影院首页| 一边摸一边抽搐一进一小说 | 国产av国产精品国产| 国产精品麻豆人妻色哟哟久久| 少妇的丰满在线观看| 少妇精品久久久久久久| 国产精品欧美亚洲77777| 久热这里只有精品99| 国产精品一区二区在线不卡| 亚洲精品美女久久av网站| 真人做人爱边吃奶动态| 肉色欧美久久久久久久蜜桃| 色老头精品视频在线观看| 最近最新免费中文字幕在线| 超碰97精品在线观看| av福利片在线| 波多野结衣av一区二区av| 日韩中文字幕视频在线看片| 国产精品久久久人人做人人爽| 精品第一国产精品| 成年人免费黄色播放视频| 搡老乐熟女国产| 国产xxxxx性猛交| 亚洲一码二码三码区别大吗| 90打野战视频偷拍视频| 国产黄频视频在线观看| 中文字幕最新亚洲高清| 亚洲国产成人一精品久久久| 99国产精品一区二区蜜桃av | 久久ye,这里只有精品| 999精品在线视频| av网站免费在线观看视频| 脱女人内裤的视频| 国产一区有黄有色的免费视频| 国产欧美日韩一区二区精品| 侵犯人妻中文字幕一二三四区| 男女之事视频高清在线观看| 亚洲国产欧美网| 国产精品国产av在线观看| 成人国产av品久久久| 日本黄色日本黄色录像| 男女高潮啪啪啪动态图| 在线观看一区二区三区激情| 老熟妇乱子伦视频在线观看| 久久中文看片网| 国产精品免费大片| 国产日韩欧美亚洲二区| 精品久久久久久电影网| 99久久人妻综合| 国产麻豆69| 韩国精品一区二区三区| 国产精品秋霞免费鲁丝片| 天天躁夜夜躁狠狠躁躁| 99久久国产精品久久久| 国产无遮挡羞羞视频在线观看| 18禁美女被吸乳视频| 一本大道久久a久久精品| 国产欧美日韩精品亚洲av| 国产精品 欧美亚洲| 国产成人系列免费观看| 精品久久久久久电影网| 久久久久久久大尺度免费视频| 国产精品成人在线| 又紧又爽又黄一区二区| 满18在线观看网站| 美国免费a级毛片| 国产欧美亚洲国产| 成人18禁高潮啪啪吃奶动态图| 99久久人妻综合| 69精品国产乱码久久久| 少妇 在线观看| 欧美精品一区二区大全| 久久精品亚洲精品国产色婷小说| 人人妻人人澡人人爽人人夜夜| 91老司机精品| 国产精品免费一区二区三区在线 | 免费人妻精品一区二区三区视频| 麻豆成人av在线观看| 最近最新免费中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 久久天躁狠狠躁夜夜2o2o| 成年人午夜在线观看视频| 精品福利永久在线观看| 老司机在亚洲福利影院| 亚洲av欧美aⅴ国产| 成年人午夜在线观看视频| 天堂中文最新版在线下载| 在线观看免费高清a一片| 国产午夜精品久久久久久| 欧美久久黑人一区二区| 人妻 亚洲 视频| 国产精品欧美亚洲77777| 国产高清激情床上av| 自拍欧美九色日韩亚洲蝌蚪91| 欧美av亚洲av综合av国产av| 国产成人欧美在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 操美女的视频在线观看| 亚洲专区中文字幕在线| 两性夫妻黄色片| 激情视频va一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 纯流量卡能插随身wifi吗| 别揉我奶头~嗯~啊~动态视频| 建设人人有责人人尽责人人享有的| a级毛片在线看网站| 久久久久网色| 久久久久久久精品吃奶| 成年人午夜在线观看视频| 男女下面插进去视频免费观看| 久久久久网色| 亚洲av欧美aⅴ国产| 97人妻天天添夜夜摸| 啪啪无遮挡十八禁网站| 国产精品久久久久久精品电影小说| 大型av网站在线播放| 国产精品久久久久久精品电影小说| 国产精品久久久久久人妻精品电影 | 精品亚洲乱码少妇综合久久| 色婷婷av一区二区三区视频| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 欧美+亚洲+日韩+国产| 最黄视频免费看| 久久久久久久大尺度免费视频| 后天国语完整版免费观看| 国产高清视频在线播放一区| 国产欧美日韩一区二区精品| 免费在线观看影片大全网站| 黄色片一级片一级黄色片| 亚洲精品在线观看二区| 日韩免费av在线播放| 亚洲精品中文字幕一二三四区 | 午夜福利乱码中文字幕| 精品欧美一区二区三区在线| 日韩视频在线欧美| 国产激情久久老熟女| 日韩欧美一区视频在线观看| 高清在线国产一区| 精品乱码久久久久久99久播| 50天的宝宝边吃奶边哭怎么回事| 亚洲成av片中文字幕在线观看| xxxhd国产人妻xxx| av有码第一页| 黄色成人免费大全| 日韩一卡2卡3卡4卡2021年| 国产欧美日韩一区二区三| 午夜91福利影院| 久久久久国内视频| 亚洲成人免费av在线播放| 亚洲情色 制服丝袜| 国产伦人伦偷精品视频| 热re99久久精品国产66热6| 日韩中文字幕视频在线看片| 欧美日本中文国产一区发布| 亚洲午夜精品一区,二区,三区| 日韩有码中文字幕| 欧美变态另类bdsm刘玥| 91九色精品人成在线观看| 久久性视频一级片| 天天影视国产精品| 免费一级毛片在线播放高清视频 | 久久久久久久大尺度免费视频| 免费久久久久久久精品成人欧美视频| 女人久久www免费人成看片| 一级a爱视频在线免费观看| 亚洲人成伊人成综合网2020| 亚洲第一青青草原| 99国产精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利乱码中文字幕| netflix在线观看网站| 亚洲五月婷婷丁香| 在线观看一区二区三区激情| 啦啦啦在线免费观看视频4| 国产精品 国内视频| 老熟女久久久| 制服人妻中文乱码| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美日韩在线播放| 多毛熟女@视频| 亚洲av电影在线进入| 国产亚洲一区二区精品| 久久精品国产99精品国产亚洲性色 | 男女边摸边吃奶| 老熟女久久久| 91老司机精品| 欧美精品一区二区大全| av天堂在线播放| 黄色a级毛片大全视频| 亚洲七黄色美女视频| 亚洲 国产 在线| 久久人妻av系列| 18禁美女被吸乳视频| 男人舔女人的私密视频| 国产aⅴ精品一区二区三区波| 狠狠婷婷综合久久久久久88av| 欧美日韩亚洲高清精品| 亚洲欧洲日产国产| 俄罗斯特黄特色一大片| 老司机靠b影院| 黄色丝袜av网址大全| 男女之事视频高清在线观看| 18禁观看日本| 精品第一国产精品| 欧美精品一区二区大全| 免费看a级黄色片| 亚洲 欧美一区二区三区| 一级片免费观看大全| 丝袜喷水一区| 少妇猛男粗大的猛烈进出视频| 亚洲五月婷婷丁香|