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

    基于GIS的潰壩洪水模擬

    2021-10-26 04:52:48高瞇瞇張永彬
    關(guān)鍵詞:潰壩元胞洪水

    高瞇瞇,張永彬

    (華北理工大學(xué) 礦業(yè)工程學(xué)院,河北 唐山 063210)

    引言

    為了便于耕種灌溉、漁民養(yǎng)殖以及滿足日常用電需求等,人類在天然湖泊上修建了許多水庫(kù)大壩。地震、沙塵暴、洪澇等自然災(zāi)害的高頻發(fā)生,使得各類極端自然災(zāi)害事件大幅增加,在遍及全球的各種自然災(zāi)害中,危害最大的當(dāng)屬洪澇災(zāi)害,屬于一種最為尋常、破壞力極大的災(zāi)害。這種災(zāi)害發(fā)生的頻率極高,影響范圍極廣,歷年來(lái)導(dǎo)致了不可估量的經(jīng)濟(jì)損失,且一直有上升的趨勢(shì)[1]。我國(guó)大壩總量位居世界第一,潰壩率亦是世界第一。因此,國(guó)家將防洪工作作為關(guān)乎民生安全保障的一項(xiàng)重要事業(yè),在水庫(kù)大壩的防洪建設(shè)方面提出了多種措施方案,在防洪減災(zāi)工作中具有巨大的影響力,為國(guó)家?guī)?lái)了巨大的經(jīng)濟(jì)效益和社會(huì)效益,與此同時(shí)也促進(jìn)了人民生活水平的改善和提升。人類為了更好地保證自身安全,使災(zāi)害損失達(dá)到最小化,應(yīng)該充分了解洪水活動(dòng)變化所有過(guò)程,深入研究洪水的運(yùn)動(dòng)規(guī)律,找出洪水的變化的特征,此類研究對(duì)人民生活和國(guó)家經(jīng)濟(jì)建設(shè)與發(fā)展具有重要意義[2,3]。

    現(xiàn)代技術(shù)在洪災(zāi)中的運(yùn)用越來(lái)越普及,因而產(chǎn)生了治理洪災(zāi)更多的可能性,愈來(lái)愈多學(xué)者關(guān)注并研究潰壩洪水淹沒(méi)問(wèn)題,并分別在研究的內(nèi)容、方法、范圍等方面獲得一定成就[4-6]。但該項(xiàng)研究仍有諸多不足之處,例如,一部分研究?jī)H展示了洪水的二維、三維可視化表達(dá)的效果,并沒(méi)有進(jìn)行詳細(xì)地分析;另一部分研究由于一些客觀條件的限制,比如:計(jì)算方法、實(shí)驗(yàn)條件等,僅僅是針對(duì)洪水淹沒(méi)范圍的一些數(shù)據(jù)分析,沒(méi)有將淹沒(méi)的范圍給出直觀的可視化展示,或者展示的效果并不理想,導(dǎo)致預(yù)警預(yù)測(cè)洪水淹沒(méi)的過(guò)程難以進(jìn)行[7-10]。在大數(shù)據(jù)時(shí)代,為了更加合理、有效地利用潰壩洪水信息資源,為水文學(xué)研究人員,防洪減災(zāi)相關(guān)業(yè)務(wù)部門的科研人員直觀了解洪水淹沒(méi)的過(guò)程以及相關(guān)信息提供技術(shù)支撐。該項(xiàng)研究開(kāi)發(fā)了集數(shù)據(jù)處理、分析、計(jì)算、模擬與預(yù)警于一體的潰壩洪水模擬系統(tǒng)。

    1模型與算法描述

    1.1 SCS模型

    SCS模型根據(jù)采樣區(qū)的實(shí)際滲水量(F)與真正的徑流量(Q)之比等于本區(qū)域該場(chǎng)降雨前的潛在入滲量(S)與潛在徑流之比的假定基礎(chǔ)上建立的[11],即:

    (1)

    假定潛在徑流量為降雨量(P)與徑流產(chǎn)生前植物截流、初滲水和填洼水構(gòu)成了集水區(qū)初損量Ia的差值,即:

    Qm=P-Ia

    (2)

    實(shí)際入滲量為降雨量減去初損和徑流量,即:

    F=P-Ia-Q

    (3)

    由式(1)~式(3)可得出:

    (4)

    假定集水區(qū)降雨的初損為潛在入滲量的20%,即:

    Ia=0.2S

    (5)

    則式(4)可寫為:

    (6)

    1.2 元胞自動(dòng)機(jī)模型

    該模型包含以下4個(gè)參數(shù)CA=。式中:L為將采樣區(qū)離散化的規(guī)則元胞聚集,該模型采用規(guī)則格網(wǎng)元胞,邊長(zhǎng)d為5 m。N為馮諾依曼鄰域關(guān)系,該元胞下一時(shí)間段的形態(tài)是由自身和相鄰元胞的形態(tài)共同確定的,在所有時(shí)間點(diǎn)該元胞均能夠與周圍的元胞存在水流的傳遞,中心元胞下一秒的形態(tài)是由其相連接的8個(gè)方向上的元胞以及自身共同決定的。P是用來(lái)描述轉(zhuǎn)換功能和元胞形態(tài)的元胞自動(dòng)機(jī)模型參數(shù)的限定集,其中包括:全局參數(shù),如元胞邊長(zhǎng)d、歷時(shí)步長(zhǎng)t;靜態(tài)參數(shù),如DEM高程和曼寧糙率系數(shù);動(dòng)態(tài)參數(shù),如降雨量、水深、水位高程等。R表示決定一個(gè)元胞到另一個(gè)元胞流向的轉(zhuǎn)換規(guī)則[12]。在單位時(shí)間步長(zhǎng)t內(nèi),可通過(guò)曼寧公式進(jìn)行計(jì)算分配的水量:

    (7)

    (8)

    式中:V為水流速度(m/s);h為水深(m);j為坡度(°);n為曼寧粗糙系數(shù);E為交換水量(m3);S為水?dāng)嗝婷娣e(m2);d為元胞邊長(zhǎng)(m)。

    1.3 曼寧公式

    曼寧公式指的是明渠道流量或速度經(jīng)驗(yàn)公式,常用于水利建設(shè)等活動(dòng)中。該項(xiàng)研究中采用其計(jì)算水流速度[13],公式如下:

    (9)

    式中:V為水流速度(m/s);k為轉(zhuǎn)換常數(shù),國(guó)際單位制中取值為1;n為糙率,綜合反映表面粗糙情況對(duì)水流的影響;Rh為水力半徑(m);S為坡度(°),表示2個(gè)格網(wǎng)之間的坡度[14]。

    1.4 D8單向流法

    該項(xiàng)研究中采用D8單向流算法決定流向,假設(shè)中心單元中的水只能匯入與之相鄰的8個(gè)不同方向的單元格中。利用最大水位高程落差法來(lái)確定水流的方向,即在3×3的DEM 單元格網(wǎng)中,計(jì)算中心單元格與各相鄰單元格的水位高程之間的關(guān)系,取與中心單元差值最大的鄰域單元格作為下游單元格流入,如圖1所示。

    圖1 D8單向流

    2案例研究

    2.1 研究區(qū)概況

    實(shí)驗(yàn)區(qū)溪源宮水庫(kù)地處福州市,常年氣溫適宜,溫暖濕潤(rùn),雨量充沛。福州市地處閩江下游,屬于沿海城市,全市海域面積高達(dá)11.09×104km2。溪源宮水庫(kù)屬中型水庫(kù),位于福州市閩侯上街溪源宮,其特殊的地理位置,時(shí)常出現(xiàn)持續(xù)較強(qiáng)降雨,連續(xù)的降雨導(dǎo)致水位上漲,降雨、融雪、潰壩崩塌等原因引發(fā)水庫(kù)沿海水位上升,水庫(kù)壩體存在著決堤的可能性。

    2.2 數(shù)據(jù)來(lái)源及預(yù)處理

    圖2所示為10 m分辨率的原始遙感影像。該項(xiàng)研究選用源于哨兵二號(hào)官網(wǎng)下載的遙感影像European Space Agency,Sentinel-2 衛(wèi)星數(shù)據(jù)(Band 2、Band3、Band4組合)以及在地理空間信息云上下載的福州市dem數(shù)據(jù)。福州市1:100 萬(wàn)矢量地圖。降水量、水深、水位高程信息均來(lái)自于福建水利信息網(wǎng)。

    圖2 10 m分辨率原始遙感影像

    圖3所示為監(jiān)督分類后的遙感影像。利用ENVI軟件對(duì)10 m分辨率的原始遙感影像進(jìn)行監(jiān)督分類,將影像中的地物大致分為六大類,定義紅色為道路,綠色為耕地,藍(lán)色為河道,黃色為居民地,鮮青色為建設(shè)用地,紫色為林地。

    圖3 監(jiān)督分類后的圖

    2.3 系統(tǒng)功能設(shè)計(jì)

    系統(tǒng)設(shè)計(jì)數(shù)據(jù)處理、數(shù)據(jù)管理、洪水模擬、數(shù)據(jù)分析、查詢、動(dòng)態(tài)推演與預(yù)警六大模塊,實(shí)現(xiàn)了數(shù)據(jù)的預(yù)處理、地圖的加載、圖層的管理、圖像的放大縮小平移、暴雨的模擬、潰壩的模擬、淹沒(méi)面積的計(jì)算、淹沒(méi)體積的計(jì)算、淹沒(méi)水深的查詢、降雨量的查詢、動(dòng)態(tài)推演與預(yù)警等重要功能。圖4所示為系統(tǒng)技術(shù)路線圖。

    圖4 系統(tǒng)技術(shù)路線圖

    系統(tǒng)設(shè)計(jì)2種不同形式的模擬即:暴雨洪水淹沒(méi)模擬與潰壩洪水淹沒(méi)模擬。圖5所表示的是洪水淹沒(méi)模擬推演的一個(gè)過(guò)程。整個(gè)系統(tǒng)是基于元胞自動(dòng)機(jī)的方法來(lái)實(shí)現(xiàn)的,為了模擬洪水徑流的過(guò)程,選取了SCS水文模型,來(lái)模擬產(chǎn)生匯流的過(guò)程,在該研究區(qū)存在幾個(gè)雨量站,經(jīng)過(guò)IDW雨量插值可以得到對(duì)應(yīng)點(diǎn)的降雨量,通過(guò)空間插值算法得到研究區(qū)面雨量。獲取得到面雨量之后,在SCS模型中進(jìn)一步進(jìn)行數(shù)據(jù)處理,因?yàn)榻德涞挠晁⒉皇侨柯涞降厣闲纬煞e水或者河流,其中一部分會(huì)留在樹上或滲透到土壤里,這部分被稱作是降雨初損。因此就要進(jìn)行凈雨量的計(jì)算,降雨初損使用20%最大下滲量來(lái)計(jì)算。在洪水淹沒(méi)模擬元胞自動(dòng)機(jī)模型中,降雨從中心柵格向相鄰柵格流入的分布形成了匯流。在這個(gè)過(guò)程中,采用D8單向流算法推算雨水的流動(dòng)方向,該方向是由中心單元與相鄰單元水位落差之間的關(guān)系進(jìn)行確定的,基本規(guī)則如下:若中心元胞水位高程低于或者等于所有鄰域元胞水位高程,則不進(jìn)行水量分配;反之選擇高差最大的鄰域元胞作為下游。雨水流動(dòng)速度采用曼寧公式確定,每一個(gè)單元看作是一個(gè)元胞,整個(gè)流域分布分成很多個(gè)單元,通過(guò)計(jì)算每個(gè)單元的水量進(jìn)而計(jì)算出整個(gè)流域的水量。

    圖5 洪水淹沒(méi)模擬推演數(shù)據(jù)流圖

    2.4 系統(tǒng)實(shí)現(xiàn)

    2.4.1洪水淹沒(méi)模擬功能的實(shí)現(xiàn)

    暴雨洪水是由較大強(qiáng)度的降雨而形成的大區(qū)域洪水,引發(fā)水庫(kù)水量迅速增加并伴隨水位急劇上升的現(xiàn)象[15];潰壩洪水是小區(qū)域大水量的輸入,壩體或者其他擋水建筑物瞬間潰決發(fā)生水體突泄所形成的洪水現(xiàn)象。基于元胞自動(dòng)機(jī)的模擬中,利用DEM高差采用D8單向流法計(jì)算坡度和流向,并利用曼寧值來(lái)計(jì)算流速,流速乘以時(shí)間就是此單元向下一單元此時(shí)間段內(nèi)流動(dòng)的水量,將當(dāng)前的淹沒(méi)情況結(jié)果進(jìn)行累加模擬。對(duì)洪水淹沒(méi)事件做出模擬響應(yīng),找到即將發(fā)生潰壩決堤的位置,確定影響范圍等。模擬結(jié)果如圖6、圖7所示。

    圖6 暴雨洪水淹沒(méi)模擬結(jié)果 圖7 潰壩洪水淹沒(méi)模擬結(jié)果

    2.4.2查詢功能的實(shí)現(xiàn)

    查詢功能具體包含:對(duì)研究區(qū)斑塊的水深信息進(jìn)行查詢、對(duì)研究區(qū)降雨量進(jìn)行查詢、對(duì)洪水淹沒(méi)模擬動(dòng)態(tài)推演過(guò)程中斑塊水深信息的查詢、對(duì)動(dòng)態(tài)推演過(guò)程中任意時(shí)刻的洪水淹沒(méi)情況進(jìn)行查詢等。提供點(diǎn)擊查詢及按照給定條件查詢,查詢結(jié)果可直觀數(shù)字顯示或者直接在地圖上進(jìn)行顯示標(biāo)記。選擇地圖上的水域,即可獲取降雨量的信息;可以通過(guò)識(shí)別功能點(diǎn)選某一區(qū)域的水深量、面降雨量;點(diǎn)擊下拉框進(jìn)行篩選不同模式下同一地段的雨水量。以米為單位,效果圖如8所示。

    圖8 任意時(shí)刻的淹沒(méi)情況查詢

    2.4.3數(shù)據(jù)分析功能的實(shí)現(xiàn)

    經(jīng)過(guò)不同模式的洪水淹沒(méi)模擬,將模擬的結(jié)果柵格圖加載進(jìn)行分析,根據(jù)用戶自定義淹沒(méi)判斷的數(shù)值,可分別計(jì)算出該區(qū)域不同模式下符合淹沒(méi)判定的淹沒(méi)區(qū)域的體積、面積、水深情況,將符合淹沒(méi)判定的淹沒(méi)區(qū)域直觀地顯示出來(lái),利于使用系統(tǒng)的研究人員能夠即刻地分析出研究區(qū)的洪水淹沒(méi)的情況,以便于做出相應(yīng)的應(yīng)對(duì)措施。效果如圖9所示。

    圖9 淹沒(méi)面積體積的計(jì)算

    2.4.4動(dòng)態(tài)推演與預(yù)警功能的實(shí)現(xiàn)

    根據(jù)暴雨洪水淹沒(méi)模擬和潰壩洪水淹沒(méi)模擬的結(jié)果,將其給予模擬的起始時(shí)間點(diǎn)發(fā)布wms 服務(wù),按照模擬時(shí)候確定的模式,模擬時(shí)間步長(zhǎng)、模擬總時(shí)長(zhǎng)規(guī)定按照多久出一景結(jié)果為前提進(jìn)行動(dòng)態(tài)推演。根據(jù)輸入的模擬時(shí)長(zhǎng)進(jìn)行動(dòng)態(tài)推演,可以通過(guò)前進(jìn)后退逐時(shí)間段播放洪水演進(jìn)的情況,隨時(shí)點(diǎn)擊暫停查看該時(shí)刻的淹沒(méi)情況,分別展示了第40 min、第80 min、第120 min、第200 min的洪水淹沒(méi)情況,效果如圖10所示。

    圖10 洪水淹沒(méi)動(dòng)態(tài)推演效果圖

    當(dāng)水量將要逼近水庫(kù)的最大臨界點(diǎn)時(shí)將會(huì)發(fā)生潰壩現(xiàn)象,系統(tǒng)會(huì)立刻給予預(yù)警提示并標(biāo)記出將要決堤的地方以及決堤后的漫流方向。系統(tǒng)提前預(yù)警預(yù)測(cè)以高亮紅色顯出出來(lái)外加提示窗彈出,以使相關(guān)部門及時(shí)做出撤離行動(dòng)及相應(yīng)的解決方案,如圖11所示。

    圖11 洪水淹沒(méi)預(yù)警

    3結(jié)論

    根據(jù)洪水災(zāi)害的特點(diǎn),選擇實(shí)驗(yàn)區(qū)福州市溪源宮水庫(kù)進(jìn)行實(shí)驗(yàn),基于C/S +B/S架構(gòu)構(gòu)建基于GIS 的潰壩洪水模擬系統(tǒng)。實(shí)現(xiàn)了洪水淹沒(méi)模擬的動(dòng)態(tài)推演與預(yù)警,該系統(tǒng)能進(jìn)行流態(tài)模擬、負(fù)荷分析、流向狀況分析,輔助決策,對(duì)突發(fā)潰壩事件做出快速響應(yīng)與預(yù)警,快速找到潰壩臨界點(diǎn),確定淹沒(méi)方向、影響范圍等。系統(tǒng)的建立為洪水淹沒(méi)分析與預(yù)警提供了技術(shù)支撐,為水庫(kù)下游人員撤離提供參考依據(jù)。

    猜你喜歡
    潰壩元胞洪水
    洪水時(shí)遇到電線低垂或折斷該怎么辦
    基于元胞自動(dòng)機(jī)下的交通事故路段仿真
    又見(jiàn)洪水(外二首)
    徐家河尾礦庫(kù)潰壩分析
    潰壩涌浪及其對(duì)重力壩影響的數(shù)值模擬
    潰壩波對(duì)單橋墩作用水力特性研究
    基于改進(jìn)控制方程的土石壩潰壩洪水演進(jìn)數(shù)值模擬
    洪水來(lái)了
    幼兒畫刊(2016年8期)2016-02-28 21:00:52
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機(jī)制
    論設(shè)計(jì)洪水計(jì)算
    国产精品久久久久成人av| 久久中文字幕一级| 热99久久久久精品小说推荐| 国产高清激情床上av| 97在线人人人人妻| 97在线人人人人妻| 老汉色av国产亚洲站长工具| 亚洲中文av在线| 一边摸一边抽搐一进一小说 | 超色免费av| av又黄又爽大尺度在线免费看| 亚洲欧美日韩高清在线视频 | 欧美日韩福利视频一区二区| 久久午夜综合久久蜜桃| 国产精品亚洲一级av第二区| 亚洲精品一二三| 天天操日日干夜夜撸| 两人在一起打扑克的视频| 日本av手机在线免费观看| 女同久久另类99精品国产91| 亚洲精品一二三| 国产一区二区三区视频了| 免费少妇av软件| 国产在线视频一区二区| 成年人黄色毛片网站| 最新美女视频免费是黄的| 9热在线视频观看99| 99久久人妻综合| 国产精品麻豆人妻色哟哟久久| 午夜精品国产一区二区电影| 国产成人一区二区三区免费视频网站| 午夜激情久久久久久久| 久久久国产一区二区| 黄色丝袜av网址大全| 黄色a级毛片大全视频| 色老头精品视频在线观看| 午夜91福利影院| 在线观看66精品国产| 精品少妇久久久久久888优播| 亚洲情色 制服丝袜| 亚洲欧美日韩高清在线视频 | 亚洲欧美色中文字幕在线| 亚洲 国产 在线| 老熟女久久久| 不卡一级毛片| 久久久久久久大尺度免费视频| 丰满饥渴人妻一区二区三| 亚洲熟女毛片儿| 亚洲三区欧美一区| 女同久久另类99精品国产91| 另类精品久久| av线在线观看网站| 91字幕亚洲| 中文字幕制服av| 久久久久久免费高清国产稀缺| 久久久精品区二区三区| 亚洲三区欧美一区| 精品国产亚洲在线| 久久婷婷成人综合色麻豆| 国产精品98久久久久久宅男小说| 色婷婷久久久亚洲欧美| 美女福利国产在线| 99精品在免费线老司机午夜| 777米奇影视久久| 国产极品粉嫩免费观看在线| 丝袜人妻中文字幕| 久久99一区二区三区| 狠狠狠狠99中文字幕| 国产免费现黄频在线看| 两个人免费观看高清视频| 久久精品国产a三级三级三级| 每晚都被弄得嗷嗷叫到高潮| 日韩中文字幕欧美一区二区| 欧美激情久久久久久爽电影 | 岛国在线观看网站| 不卡av一区二区三区| 亚洲精品久久午夜乱码| 精品国产乱子伦一区二区三区| 国产亚洲欧美精品永久| 桃花免费在线播放| 免费观看a级毛片全部| 满18在线观看网站| 久久久欧美国产精品| 国产日韩欧美亚洲二区| 伊人久久大香线蕉亚洲五| 亚洲欧美精品综合一区二区三区| 色综合婷婷激情| 久久人人爽av亚洲精品天堂| 国产激情久久老熟女| 成人国产av品久久久| 精品一区二区三区av网在线观看 | 又紧又爽又黄一区二区| 成人永久免费在线观看视频 | 少妇精品久久久久久久| 999久久久精品免费观看国产| 91老司机精品| 一区二区av电影网| www.999成人在线观看| 亚洲美女黄片视频| 国产成+人综合+亚洲专区| 高清在线国产一区| 一进一出抽搐动态| 国产欧美日韩综合在线一区二区| 精品国产乱码久久久久久男人| 国产一区二区三区综合在线观看| h视频一区二区三区| 精品亚洲乱码少妇综合久久| 国产无遮挡羞羞视频在线观看| 一本久久精品| 天堂8中文在线网| 国产片内射在线| 亚洲免费av在线视频| 在线观看免费高清a一片| 日日摸夜夜添夜夜添小说| 国产高清videossex| 欧美乱妇无乱码| 满18在线观看网站| 日韩视频一区二区在线观看| 19禁男女啪啪无遮挡网站| 欧美变态另类bdsm刘玥| 日韩 欧美 亚洲 中文字幕| 超碰97精品在线观看| 岛国在线观看网站| 9热在线视频观看99| 王馨瑶露胸无遮挡在线观看| 我的亚洲天堂| 夫妻午夜视频| 国产精品久久电影中文字幕 | 日韩欧美三级三区| 少妇精品久久久久久久| 新久久久久国产一级毛片| 下体分泌物呈黄色| 97人妻天天添夜夜摸| 国产精品免费视频内射| 悠悠久久av| 国产成人系列免费观看| 久久毛片免费看一区二区三区| 成年人免费黄色播放视频| 18禁美女被吸乳视频| 久久久久久久久久久久大奶| 国产成人精品久久二区二区91| 丝袜人妻中文字幕| 日韩一卡2卡3卡4卡2021年| 久久久精品国产亚洲av高清涩受| 大陆偷拍与自拍| 精品久久久久久久毛片微露脸| cao死你这个sao货| 精品少妇黑人巨大在线播放| 国产精品九九99| 中文字幕制服av| 啦啦啦免费观看视频1| 国产在线视频一区二区| 精品少妇黑人巨大在线播放| 99在线人妻在线中文字幕 | 国产精品麻豆人妻色哟哟久久| 精品一区二区三卡| 免费在线观看完整版高清| 久久亚洲精品不卡| 亚洲精品在线观看二区| 国产成人系列免费观看| 国产精品香港三级国产av潘金莲| 热99久久久久精品小说推荐| 一本一本久久a久久精品综合妖精| 99re6热这里在线精品视频| 欧美黑人欧美精品刺激| 91麻豆精品激情在线观看国产 | 老司机福利观看| 自拍欧美九色日韩亚洲蝌蚪91| 91成人精品电影| 国产一区二区三区视频了| 丝袜在线中文字幕| 极品少妇高潮喷水抽搐| 在线观看人妻少妇| 亚洲国产毛片av蜜桃av| 国产免费av片在线观看野外av| 日韩欧美三级三区| 老汉色av国产亚洲站长工具| 人妻一区二区av| 美女高潮到喷水免费观看| 人人妻人人澡人人爽人人夜夜| 亚洲av第一区精品v没综合| 桃红色精品国产亚洲av| 一区在线观看完整版| 18禁黄网站禁片午夜丰满| 女警被强在线播放| 成年人黄色毛片网站| 黄片大片在线免费观看| 国产精品成人在线| 国产av又大| aaaaa片日本免费| 成人国语在线视频| 黑人猛操日本美女一级片| 看免费av毛片| 黑人操中国人逼视频| 国产精品久久久久成人av| 亚洲午夜精品一区,二区,三区| 少妇 在线观看| 自线自在国产av| 在线观看一区二区三区激情| 欧美变态另类bdsm刘玥| 大香蕉久久网| 少妇裸体淫交视频免费看高清 | 欧美成人午夜精品| 久久久精品区二区三区| 韩国精品一区二区三区| 国产不卡一卡二| 久久天堂一区二区三区四区| 色94色欧美一区二区| 性色av乱码一区二区三区2| 亚洲天堂av无毛| 国产成人av教育| 亚洲欧洲日产国产| 757午夜福利合集在线观看| 在线观看人妻少妇| 亚洲精品粉嫩美女一区| 日韩欧美一区视频在线观看| 精品久久久久久久毛片微露脸| 日本撒尿小便嘘嘘汇集6| 我的亚洲天堂| 18禁观看日本| 国产男女超爽视频在线观看| 午夜老司机福利片| 咕卡用的链子| 久久狼人影院| 欧美精品一区二区大全| av又黄又爽大尺度在线免费看| 欧美 亚洲 国产 日韩一| 久久久久久久大尺度免费视频| 啦啦啦免费观看视频1| 日日摸夜夜添夜夜添小说| 亚洲欧美一区二区三区久久| 国产精品国产高清国产av | 脱女人内裤的视频| 亚洲欧美一区二区三区久久| 国产精品一区二区精品视频观看| 每晚都被弄得嗷嗷叫到高潮| av欧美777| 天天添夜夜摸| 亚洲熟女精品中文字幕| 欧美中文综合在线视频| 啦啦啦在线免费观看视频4| www日本在线高清视频| 免费在线观看影片大全网站| 天堂动漫精品| 亚洲色图 男人天堂 中文字幕| 亚洲情色 制服丝袜| 一本久久精品| www.自偷自拍.com| 一区二区日韩欧美中文字幕| 国产精品自产拍在线观看55亚洲 | 啦啦啦 在线观看视频| 欧美精品人与动牲交sv欧美| 欧美黄色淫秽网站| 天堂8中文在线网| 在线观看66精品国产| 亚洲成人免费av在线播放| 亚洲熟妇熟女久久| 91av网站免费观看| 中文欧美无线码| 亚洲国产中文字幕在线视频| 91精品国产国语对白视频| 老司机在亚洲福利影院| 亚洲欧美激情在线| 一个人免费看片子| 五月开心婷婷网| 亚洲av国产av综合av卡| e午夜精品久久久久久久| 久久人人爽av亚洲精品天堂| 国产精品久久久久成人av| 黄色 视频免费看| 久久亚洲精品不卡| 国内毛片毛片毛片毛片毛片| 久久免费观看电影| 精品久久久久久电影网| 777米奇影视久久| 亚洲午夜理论影院| 久久婷婷成人综合色麻豆| 啦啦啦免费观看视频1| 国产精品av久久久久免费| 69精品国产乱码久久久| 高清av免费在线| 国产成人精品久久二区二区免费| 老熟妇乱子伦视频在线观看| 国产无遮挡羞羞视频在线观看| 亚洲少妇的诱惑av| 欧美黄色片欧美黄色片| 亚洲av欧美aⅴ国产| 欧美黄色淫秽网站| 怎么达到女性高潮| 九色亚洲精品在线播放| 超碰97精品在线观看| 美女午夜性视频免费| 狠狠婷婷综合久久久久久88av| 男人操女人黄网站| 精品午夜福利视频在线观看一区 | 搡老乐熟女国产| 一边摸一边抽搐一进一出视频| 国产日韩欧美视频二区| 蜜桃国产av成人99| 国产高清视频在线播放一区| 嫩草影视91久久| 美女高潮喷水抽搐中文字幕| 久久天躁狠狠躁夜夜2o2o| 成人三级做爰电影| 黄色毛片三级朝国网站| 在线观看www视频免费| 久久精品熟女亚洲av麻豆精品| 中文字幕人妻丝袜制服| 69精品国产乱码久久久| 久久人人爽av亚洲精品天堂| 变态另类成人亚洲欧美熟女 | 巨乳人妻的诱惑在线观看| 成人免费观看视频高清| 欧美乱妇无乱码| 久久久精品94久久精品| 女性被躁到高潮视频| 免费少妇av软件| 91成年电影在线观看| 制服诱惑二区| 久久久久久久久免费视频了| 亚洲欧美精品综合一区二区三区| avwww免费| 久久久久久久精品吃奶| 一本—道久久a久久精品蜜桃钙片| 中文字幕高清在线视频| 黄色 视频免费看| 男女床上黄色一级片免费看| 成人手机av| www.熟女人妻精品国产| 在线十欧美十亚洲十日本专区| 色在线成人网| 90打野战视频偷拍视频| bbb黄色大片| 久久久久久久久免费视频了| 成人特级黄色片久久久久久久 | 夜夜爽天天搞| 久久久久久亚洲精品国产蜜桃av| 极品人妻少妇av视频| av片东京热男人的天堂| 麻豆国产av国片精品| 视频区欧美日本亚洲| 久久久久网色| 一本久久精品| 欧美性长视频在线观看| 国产高清激情床上av| 国产男女超爽视频在线观看| 波多野结衣一区麻豆| 日韩大片免费观看网站| 少妇猛男粗大的猛烈进出视频| 成人特级黄色片久久久久久久 | 五月开心婷婷网| 欧美日韩一级在线毛片| 国产精品电影一区二区三区 | 蜜桃国产av成人99| 成年人免费黄色播放视频| 亚洲精品自拍成人| 国产成人啪精品午夜网站| 老司机亚洲免费影院| 91成人精品电影| 视频区欧美日本亚洲| 99香蕉大伊视频| 亚洲精品国产色婷婷电影| 欧美激情久久久久久爽电影 | 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 欧美乱妇无乱码| 精品一区二区三区av网在线观看 | 欧美精品高潮呻吟av久久| 90打野战视频偷拍视频| 午夜老司机福利片| 色综合婷婷激情| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 亚洲色图 男人天堂 中文字幕| 91老司机精品| 飞空精品影院首页| 99国产精品免费福利视频| 久久香蕉激情| 啦啦啦视频在线资源免费观看| 搡老岳熟女国产| a级毛片在线看网站| 亚洲国产中文字幕在线视频| 麻豆成人av在线观看| 人妻 亚洲 视频| 亚洲 欧美一区二区三区| 老汉色∧v一级毛片| 99精品欧美一区二区三区四区| 一区在线观看完整版| 真人做人爱边吃奶动态| 97在线人人人人妻| 18禁美女被吸乳视频| 午夜成年电影在线免费观看| 欧美精品亚洲一区二区| 日本欧美视频一区| 啪啪无遮挡十八禁网站| 人妻久久中文字幕网| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 高清毛片免费观看视频网站 | 亚洲精品国产色婷婷电影| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 国产熟女午夜一区二区三区| 久久国产精品人妻蜜桃| 老熟妇仑乱视频hdxx| 丝袜美腿诱惑在线| 婷婷丁香在线五月| www日本在线高清视频| 另类亚洲欧美激情| 丝袜人妻中文字幕| av国产精品久久久久影院| 丝瓜视频免费看黄片| 欧美日韩亚洲综合一区二区三区_| 高清av免费在线| 国产日韩欧美视频二区| 国产在线免费精品| 国产成人欧美在线观看 | 欧美人与性动交α欧美软件| 精品一区二区三卡| 国产黄频视频在线观看| 久久精品国产综合久久久| 久久国产精品男人的天堂亚洲| 超色免费av| 少妇猛男粗大的猛烈进出视频| 国产又爽黄色视频| 肉色欧美久久久久久久蜜桃| 夫妻午夜视频| 少妇粗大呻吟视频| 亚洲av日韩精品久久久久久密| 免费不卡黄色视频| 国产老妇伦熟女老妇高清| 国产一卡二卡三卡精品| 亚洲色图 男人天堂 中文字幕| 美女福利国产在线| tube8黄色片| 1024视频免费在线观看| 久久 成人 亚洲| 亚洲少妇的诱惑av| 亚洲国产欧美一区二区综合| 黄片播放在线免费| 三上悠亚av全集在线观看| 精品免费久久久久久久清纯 | 欧美国产精品一级二级三级| 黑人巨大精品欧美一区二区mp4| 国产在线视频一区二区| 热99re8久久精品国产| a级毛片黄视频| 亚洲精品美女久久久久99蜜臀| 日本vs欧美在线观看视频| 一边摸一边抽搐一进一出视频| 亚洲,欧美精品.| 国产成人一区二区三区免费视频网站| 高清欧美精品videossex| 国产精品国产av在线观看| 午夜老司机福利片| 欧美精品亚洲一区二区| 女性生殖器流出的白浆| 一本—道久久a久久精品蜜桃钙片| 久久国产精品人妻蜜桃| 亚洲精品成人av观看孕妇| 啦啦啦 在线观看视频| 欧美一级毛片孕妇| 久久婷婷成人综合色麻豆| 精品免费久久久久久久清纯 | 伦理电影免费视频| 天天影视国产精品| 99久久精品国产亚洲精品| 成年人午夜在线观看视频| 亚洲熟妇熟女久久| 亚洲综合色网址| 欧美日韩亚洲高清精品| 人成视频在线观看免费观看| 中文字幕av电影在线播放| tube8黄色片| 男人操女人黄网站| 久久久久网色| 欧美激情极品国产一区二区三区| 国产精品欧美亚洲77777| 成人国产一区最新在线观看| a级毛片在线看网站| 丰满少妇做爰视频| 欧美黑人欧美精品刺激| 国产单亲对白刺激| 亚洲精品在线美女| 黄频高清免费视频| 脱女人内裤的视频| 久久精品aⅴ一区二区三区四区| aaaaa片日本免费| 在线观看www视频免费| 男女免费视频国产| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 日韩中文字幕视频在线看片| 嫩草影视91久久| 视频在线观看一区二区三区| 夜夜骑夜夜射夜夜干| 日本精品一区二区三区蜜桃| 欧美中文综合在线视频| 日韩欧美国产一区二区入口| tocl精华| 一二三四在线观看免费中文在| 90打野战视频偷拍视频| 人人妻人人添人人爽欧美一区卜| h视频一区二区三区| av片东京热男人的天堂| 飞空精品影院首页| 18在线观看网站| 美女国产高潮福利片在线看| 一个人免费在线观看的高清视频| 丝瓜视频免费看黄片| 欧美在线一区亚洲| 亚洲av片天天在线观看| 精品人妻熟女毛片av久久网站| 成人18禁在线播放| 新久久久久国产一级毛片| 在线观看免费午夜福利视频| 精品少妇一区二区三区视频日本电影| 一夜夜www| 久久精品国产综合久久久| 国产精品九九99| 欧美大码av| 欧美黑人精品巨大| 少妇的丰满在线观看| av线在线观看网站| 别揉我奶头~嗯~啊~动态视频| 黄网站色视频无遮挡免费观看| cao死你这个sao货| 国产成+人综合+亚洲专区| av有码第一页| 精品一区二区三区av网在线观看 | 精品久久久久久久毛片微露脸| 色视频在线一区二区三区| 久久久欧美国产精品| 欧美日韩亚洲国产一区二区在线观看 | 精品乱码久久久久久99久播| 成人国语在线视频| 欧美日韩亚洲高清精品| 欧美精品一区二区免费开放| 变态另类成人亚洲欧美熟女 | a级片在线免费高清观看视频| 午夜激情久久久久久久| 新久久久久国产一级毛片| 国产成人av教育| 中文字幕制服av| 如日韩欧美国产精品一区二区三区| 国产精品98久久久久久宅男小说| 满18在线观看网站| 一本大道久久a久久精品| 日韩视频一区二区在线观看| 国产成人免费无遮挡视频| 国产不卡av网站在线观看| 夫妻午夜视频| 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频 | 在线天堂中文资源库| 91大片在线观看| 久热这里只有精品99| 夫妻午夜视频| 久久精品熟女亚洲av麻豆精品| 可以免费在线观看a视频的电影网站| 97在线人人人人妻| 国产成人一区二区三区免费视频网站| 国产成人欧美| 大香蕉久久成人网| 中文亚洲av片在线观看爽 | 男女下面插进去视频免费观看| 欧美 亚洲 国产 日韩一| 少妇裸体淫交视频免费看高清 | 国产成人啪精品午夜网站| 国产又爽黄色视频| 18禁裸乳无遮挡动漫免费视频| 国产日韩欧美亚洲二区| 18禁黄网站禁片午夜丰满| 日韩免费av在线播放| 美女主播在线视频| 飞空精品影院首页| 99精品久久久久人妻精品| 飞空精品影院首页| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 久久性视频一级片| 天堂俺去俺来也www色官网| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 一级黄色大片毛片| 视频在线观看一区二区三区| 麻豆乱淫一区二区| 国产精品偷伦视频观看了| 一区二区三区激情视频| 欧美亚洲 丝袜 人妻 在线| 成人特级黄色片久久久久久久 | www.精华液| 欧美黑人精品巨大| 在线播放国产精品三级| 考比视频在线观看| 91大片在线观看| 国产在线视频一区二区| 窝窝影院91人妻| 90打野战视频偷拍视频| 一本综合久久免费| 国产成人精品久久二区二区91| 国产精品98久久久久久宅男小说| 亚洲精品自拍成人| 久久久久久免费高清国产稀缺| 美女主播在线视频| av天堂在线播放| 亚洲第一欧美日韩一区二区三区 | 狠狠婷婷综合久久久久久88av| 国产老妇伦熟女老妇高清| 一级,二级,三级黄色视频| 久久国产精品影院| 香蕉国产在线看| 天天躁夜夜躁狠狠躁躁| av片东京热男人的天堂|