• <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ì)算
    欧美日本中文国产一区发布| 天天影视国产精品| 欧美 亚洲 国产 日韩一| 久久久久网色| 一区二区三区激情视频| 纯流量卡能插随身wifi吗| 亚洲美女搞黄在线观看| avwww免费| 亚洲精品久久成人aⅴ小说| 国产深夜福利视频在线观看| 嫩草影视91久久| 久久久久久久大尺度免费视频| 欧美精品亚洲一区二区| 不卡av一区二区三区| 精品国产国语对白av| 人妻人人澡人人爽人人| 丰满少妇做爰视频| 国产精品蜜桃在线观看| 91精品三级在线观看| 亚洲av欧美aⅴ国产| 在线观看免费午夜福利视频| 亚洲欧美一区二区三区黑人| 国产深夜福利视频在线观看| 国产精品麻豆人妻色哟哟久久| 人妻一区二区av| 制服丝袜香蕉在线| 亚洲欧美中文字幕日韩二区| av国产精品久久久久影院| av天堂久久9| 国产亚洲最大av| 久久人人爽av亚洲精品天堂| 中文字幕av电影在线播放| 色婷婷久久久亚洲欧美| 亚洲国产欧美日韩在线播放| 午夜精品国产一区二区电影| 国产精品久久久久久久久免| 九色亚洲精品在线播放| 在线 av 中文字幕| 韩国av在线不卡| 九九爱精品视频在线观看| 女人高潮潮喷娇喘18禁视频| 丝袜在线中文字幕| 久久狼人影院| 99国产精品免费福利视频| 女性生殖器流出的白浆| 美女高潮到喷水免费观看| 天天添夜夜摸| 国产在线视频一区二区| 99精国产麻豆久久婷婷| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 观看av在线不卡| 亚洲国产av新网站| 777米奇影视久久| 亚洲欧洲精品一区二区精品久久久 | 赤兔流量卡办理| 国产毛片在线视频| 国产一区二区三区综合在线观看| 十八禁网站网址无遮挡| 1024香蕉在线观看| 考比视频在线观看| 蜜桃在线观看..| 秋霞伦理黄片| 亚洲成av片中文字幕在线观看| 国产高清不卡午夜福利| 亚洲人成77777在线视频| 天天操日日干夜夜撸| 天天操日日干夜夜撸| av免费观看日本| 日韩不卡一区二区三区视频在线| 一本色道久久久久久精品综合| 在线观看三级黄色| 久久综合国产亚洲精品| 国产av国产精品国产| 黄色毛片三级朝国网站| 日韩免费高清中文字幕av| 欧美精品亚洲一区二区| 超碰97精品在线观看| 国产片特级美女逼逼视频| 狂野欧美激情性bbbbbb| 欧美亚洲日本最大视频资源| 欧美在线黄色| videosex国产| 亚洲激情五月婷婷啪啪| 免费观看人在逋| 精品酒店卫生间| 日本黄色日本黄色录像| 久久精品国产亚洲av高清一级| 免费观看a级毛片全部| 成人三级做爰电影| 校园人妻丝袜中文字幕| 欧美精品av麻豆av| 久久久久人妻精品一区果冻| 亚洲精品av麻豆狂野| 美国免费a级毛片| 成人影院久久| 国产欧美日韩综合在线一区二区| 国产人伦9x9x在线观看| 中文字幕av电影在线播放| 无限看片的www在线观看| 国产精品99久久99久久久不卡 | 国产欧美日韩综合在线一区二区| 中文字幕最新亚洲高清| 国产成人a∨麻豆精品| 午夜91福利影院| 亚洲av福利一区| 久久人人爽av亚洲精品天堂| 99香蕉大伊视频| 午夜福利网站1000一区二区三区| 伊人亚洲综合成人网| 亚洲国产av影院在线观看| 欧美国产精品一级二级三级| 最近中文字幕高清免费大全6| 亚洲精品日本国产第一区| 亚洲美女搞黄在线观看| av线在线观看网站| 国产一卡二卡三卡精品 | 国产成人啪精品午夜网站| 看非洲黑人一级黄片| 国产成人av激情在线播放| 久久青草综合色| 一本久久精品| 飞空精品影院首页| 美女中出高潮动态图| 亚洲中文av在线| 超色免费av| 美女主播在线视频| 男女高潮啪啪啪动态图| 人人澡人人妻人| 丰满少妇做爰视频| 国产成人a∨麻豆精品| 一级黄片播放器| 精品亚洲乱码少妇综合久久| 欧美日韩亚洲高清精品| 搡老岳熟女国产| 亚洲国产精品999| 大话2 男鬼变身卡| 亚洲一区二区三区欧美精品| 国产成人精品无人区| 国产熟女午夜一区二区三区| 成人国语在线视频| 国产欧美日韩综合在线一区二区| 国产黄色免费在线视频| 久久99一区二区三区| 国产国语露脸激情在线看| 国产精品免费大片| 美女国产高潮福利片在线看| 国产亚洲av片在线观看秒播厂| 国产精品二区激情视频| 少妇精品久久久久久久| 十八禁网站网址无遮挡| 麻豆av在线久日| 亚洲精品乱久久久久久| 久久综合国产亚洲精品| 国产老妇伦熟女老妇高清| 中文字幕制服av| 精品少妇黑人巨大在线播放| 中文字幕色久视频| 国产成人欧美| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 一级a爱视频在线免费观看| 精品国产一区二区三区四区第35| 日韩,欧美,国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 日韩人妻精品一区2区三区| 又大又爽又粗| 一区在线观看完整版| 少妇人妻 视频| 国产一区二区 视频在线| av不卡在线播放| 中文字幕制服av| 黄色视频不卡| 伦理电影免费视频| 在线观看人妻少妇| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 久久久国产欧美日韩av| 欧美精品一区二区免费开放| 亚洲一区中文字幕在线| 欧美激情极品国产一区二区三区| 一级片'在线观看视频| 妹子高潮喷水视频| 男女免费视频国产| 国产成人a∨麻豆精品| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| 国产免费一区二区三区四区乱码| kizo精华| 爱豆传媒免费全集在线观看| 中文字幕制服av| 国产成人精品无人区| 午夜免费鲁丝| 在线观看免费视频网站a站| 免费观看a级毛片全部| 久久精品亚洲熟妇少妇任你| 亚洲中文av在线| 2018国产大陆天天弄谢| 人人妻,人人澡人人爽秒播 | 国产精品一国产av| 国产精品一区二区精品视频观看| 欧美最新免费一区二区三区| 黄色怎么调成土黄色| 日本91视频免费播放| 日韩大码丰满熟妇| 亚洲四区av| 久久久久精品人妻al黑| 精品一区在线观看国产| 晚上一个人看的免费电影| 亚洲国产精品999| 免费观看av网站的网址| 国产成人精品无人区| 男女边摸边吃奶| 日韩av免费高清视频| 国产精品久久久久久精品电影小说| 国产精品麻豆人妻色哟哟久久| 男女免费视频国产| 欧美中文综合在线视频| videos熟女内射| 亚洲人成电影观看| 中文字幕精品免费在线观看视频| 91精品三级在线观看| 亚洲综合色网址| 国产一级毛片在线| 最新的欧美精品一区二区| 日韩欧美精品免费久久| 欧美日本中文国产一区发布| 久久影院123| 一个人免费看片子| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩另类电影网站| 51午夜福利影视在线观看| 亚洲精品久久午夜乱码| 国产成人系列免费观看| 亚洲美女视频黄频| 久久精品久久久久久久性| 欧美激情高清一区二区三区 | 夫妻性生交免费视频一级片| 啦啦啦中文免费视频观看日本| 在线观看免费视频网站a站| 亚洲美女视频黄频| 国产日韩欧美视频二区| 欧美精品一区二区大全| 精品少妇内射三级| 男女床上黄色一级片免费看| 一级片'在线观看视频| 亚洲专区中文字幕在线 | 中文字幕制服av| 久久天躁狠狠躁夜夜2o2o | 少妇精品久久久久久久| 精品一区二区免费观看| 在线亚洲精品国产二区图片欧美| 男女边摸边吃奶| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 久久ye,这里只有精品| 欧美av亚洲av综合av国产av | 免费少妇av软件| 99香蕉大伊视频| 亚洲精品第二区| 久久人妻熟女aⅴ| 黄色毛片三级朝国网站| 精品人妻熟女毛片av久久网站| 国产亚洲午夜精品一区二区久久| 亚洲少妇的诱惑av| 日本一区二区免费在线视频| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 一本久久精品| 国产av精品麻豆| 国产成人午夜福利电影在线观看| 国产激情久久老熟女| 99久久人妻综合| 亚洲精品,欧美精品| a级毛片黄视频| 色94色欧美一区二区| 人人澡人人妻人| 欧美黑人精品巨大| 国产国语露脸激情在线看| 色视频在线一区二区三区| 久久人妻熟女aⅴ| 最新的欧美精品一区二区| 亚洲情色 制服丝袜| 悠悠久久av| 精品国产露脸久久av麻豆| 777米奇影视久久| 老司机在亚洲福利影院| 国产精品国产av在线观看| 国产女主播在线喷水免费视频网站| 国产精品亚洲av一区麻豆 | 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区| 精品国产乱码久久久久久小说| 丁香六月欧美| 国产野战对白在线观看| 精品福利永久在线观看| 国产人伦9x9x在线观看| 欧美亚洲日本最大视频资源| 丝袜美腿诱惑在线| 大片电影免费在线观看免费| 丰满饥渴人妻一区二区三| 欧美激情极品国产一区二区三区| 十分钟在线观看高清视频www| 亚洲精品日本国产第一区| a级片在线免费高清观看视频| 别揉我奶头~嗯~啊~动态视频 | 汤姆久久久久久久影院中文字幕| www日本在线高清视频| 精品酒店卫生间| 人人澡人人妻人| 国产又色又爽无遮挡免| 国产亚洲午夜精品一区二区久久| 日韩一区二区视频免费看| 亚洲成国产人片在线观看| 男人舔女人的私密视频| 亚洲av成人精品一二三区| 999精品在线视频| 亚洲人成电影观看| 亚洲欧美激情在线| 亚洲精品自拍成人| 男人爽女人下面视频在线观看| 成人亚洲欧美一区二区av| 免费观看性生交大片5| 一本一本久久a久久精品综合妖精| 精品一区在线观看国产| 亚洲视频免费观看视频| 中国国产av一级| 国产精品国产av在线观看| 亚洲国产中文字幕在线视频| 精品视频人人做人人爽| 欧美中文综合在线视频| 国产av一区二区精品久久| 一级片'在线观看视频| 国产一区有黄有色的免费视频| 亚洲av男天堂| 国产日韩欧美亚洲二区| 一本色道久久久久久精品综合| 女人高潮潮喷娇喘18禁视频| 成人国产av品久久久| 少妇被粗大猛烈的视频| 男人舔女人的私密视频| 99香蕉大伊视频| 悠悠久久av| 亚洲色图 男人天堂 中文字幕| 亚洲三区欧美一区| 香蕉国产在线看| 性色av一级| 精品福利永久在线观看| 97精品久久久久久久久久精品| 伊人亚洲综合成人网| 亚洲国产精品成人久久小说| 人成视频在线观看免费观看| 亚洲欧美中文字幕日韩二区| 亚洲一区二区三区欧美精品| 黄网站色视频无遮挡免费观看| 日韩大码丰满熟妇| 亚洲av综合色区一区| 欧美精品一区二区免费开放| 国产精品蜜桃在线观看| 一区二区av电影网| 国产在线一区二区三区精| 日韩精品有码人妻一区| 99久久人妻综合| 精品久久久久久电影网| 日本91视频免费播放| 一本—道久久a久久精品蜜桃钙片| 97在线人人人人妻| 久久久久视频综合| 美女中出高潮动态图| 国产亚洲欧美精品永久| 人妻 亚洲 视频| 男女高潮啪啪啪动态图| 国产毛片在线视频| 啦啦啦在线免费观看视频4| 啦啦啦 在线观看视频| 最近中文字幕高清免费大全6| 亚洲精品美女久久av网站| 看免费av毛片| 女性生殖器流出的白浆| 王馨瑶露胸无遮挡在线观看| 男女床上黄色一级片免费看| 欧美日韩成人在线一区二区| 欧美精品一区二区大全| 国产一区亚洲一区在线观看| 激情视频va一区二区三区| 午夜福利在线免费观看网站| 久久久精品94久久精品| 成年av动漫网址| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 亚洲成人国产一区在线观看 | 日本vs欧美在线观看视频| 丰满饥渴人妻一区二区三| 国产成人精品久久久久久| 高清不卡的av网站| 免费看不卡的av| 各种免费的搞黄视频| 女人被躁到高潮嗷嗷叫费观| 久久久亚洲精品成人影院| 天天躁日日躁夜夜躁夜夜| 黄片小视频在线播放| 777米奇影视久久| 亚洲图色成人| 只有这里有精品99| 飞空精品影院首页| 日本vs欧美在线观看视频| 免费黄频网站在线观看国产| 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| 9色porny在线观看| 成人免费观看视频高清| 美国免费a级毛片| 久久精品国产亚洲av高清一级| 丝袜脚勾引网站| 18在线观看网站| 婷婷成人精品国产| 亚洲美女视频黄频| 巨乳人妻的诱惑在线观看| 男女边吃奶边做爰视频| 国产毛片在线视频| 91精品伊人久久大香线蕉| 叶爱在线成人免费视频播放| 国产片特级美女逼逼视频| 亚洲精品自拍成人| 欧美成人精品欧美一级黄| 亚洲国产av新网站| 欧美97在线视频| 国产精品 欧美亚洲| 亚洲精华国产精华液的使用体验| 91老司机精品| 人人妻人人澡人人爽人人夜夜| 亚洲综合精品二区| 母亲3免费完整高清在线观看| 亚洲精品第二区| 亚洲色图 男人天堂 中文字幕| 国产爽快片一区二区三区| 国产麻豆69| 免费黄频网站在线观看国产| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区| 大片免费播放器 马上看| 久久人妻熟女aⅴ| 最近2019中文字幕mv第一页| 国产 一区精品| av有码第一页| 国精品久久久久久国模美| 老司机影院毛片| 五月天丁香电影| 国产黄色免费在线视频| 18禁观看日本| 国产男人的电影天堂91| 免费在线观看完整版高清| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 日韩伦理黄色片| 国产精品无大码| 曰老女人黄片| 亚洲欧美一区二区三区黑人| 美国免费a级毛片| 国产男女超爽视频在线观看| 国产99久久九九免费精品| 人妻一区二区av| 日韩av免费高清视频| 成年av动漫网址| 婷婷色综合大香蕉| 精品少妇一区二区三区视频日本电影 | 在线观看一区二区三区激情| 亚洲精品美女久久av网站| 欧美中文综合在线视频| 欧美在线黄色| 在线观看三级黄色| 成人午夜精彩视频在线观看| 国精品久久久久久国模美| 中文字幕人妻丝袜制服| 亚洲精品第二区| 欧美日韩福利视频一区二区| av在线老鸭窝| 亚洲国产精品一区三区| 天天添夜夜摸| 侵犯人妻中文字幕一二三四区| 亚洲色图 男人天堂 中文字幕| 97在线人人人人妻| 男女高潮啪啪啪动态图| 成人免费观看视频高清| 国产精品久久久久久精品电影小说| 91国产中文字幕| 精品亚洲乱码少妇综合久久| 美女扒开内裤让男人捅视频| 久久久久精品人妻al黑| 2021少妇久久久久久久久久久| 伊人亚洲综合成人网| 亚洲欧美一区二区三区国产| 亚洲一码二码三码区别大吗| 欧美国产精品一级二级三级| 最近手机中文字幕大全| 人体艺术视频欧美日本| 亚洲自偷自拍图片 自拍| 大话2 男鬼变身卡| 国产精品久久久久成人av| 美女脱内裤让男人舔精品视频| 欧美 亚洲 国产 日韩一| 国产免费福利视频在线观看| 久久久久网色| 欧美精品人与动牲交sv欧美| 狂野欧美激情性bbbbbb| 最新的欧美精品一区二区| 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 欧美日韩精品网址| 亚洲av欧美aⅴ国产| 老司机在亚洲福利影院| 久久久久久久久久久久大奶| 国产精品久久久久久精品电影小说| 最近最新中文字幕大全免费视频 | 久久精品国产综合久久久| 色网站视频免费| 一本久久精品| 亚洲欧美精品自产自拍| 2021少妇久久久久久久久久久| 在线看a的网站| 老司机在亚洲福利影院| 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 亚洲精品日韩在线中文字幕| 国产精品免费视频内射| 日韩大码丰满熟妇| 热re99久久精品国产66热6| 成人手机av| 晚上一个人看的免费电影| 中文精品一卡2卡3卡4更新| 十八禁网站网址无遮挡| 成年人午夜在线观看视频| 婷婷色综合大香蕉| avwww免费| 黄色毛片三级朝国网站| 久久免费观看电影| 亚洲伊人色综图| 中文字幕最新亚洲高清| 美女脱内裤让男人舔精品视频| 亚洲三区欧美一区| 亚洲精品乱久久久久久| 一本一本久久a久久精品综合妖精| 无限看片的www在线观看| 亚洲美女黄色视频免费看| 一级爰片在线观看| av在线播放精品| 亚洲av欧美aⅴ国产| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久 | 日韩精品有码人妻一区| 国产精品久久久人人做人人爽| 精品卡一卡二卡四卡免费| 久久天躁狠狠躁夜夜2o2o | av.在线天堂| 韩国高清视频一区二区三区| 国产精品人妻久久久影院| 久久久久精品人妻al黑| 久久久久精品性色| 久久这里只有精品19| av视频免费观看在线观看| 亚洲国产欧美网| 一区二区三区精品91| 在线观看一区二区三区激情| 国产av国产精品国产| 久久免费观看电影| 涩涩av久久男人的天堂| 天堂8中文在线网| 国产成人精品久久二区二区91 | 国产精品一国产av| 男女无遮挡免费网站观看| 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| 亚洲伊人色综图| 9色porny在线观看| 满18在线观看网站| 亚洲专区中文字幕在线 | 国产精品一区二区在线观看99| 麻豆乱淫一区二区| 精品久久蜜臀av无| 伦理电影大哥的女人| 99热全是精品| 国产精品人妻久久久影院| 国产在线免费精品| 一级黄片播放器| 精品久久蜜臀av无| 美女高潮到喷水免费观看| 亚洲精品第二区| 少妇的丰满在线观看| 亚洲精品美女久久av网站| 国产精品99久久99久久久不卡 | 免费看不卡的av| 超碰成人久久| 又大又黄又爽视频免费| 久久天堂一区二区三区四区| 欧美 亚洲 国产 日韩一| 国产欧美日韩综合在线一区二区| 国产不卡av网站在线观看| av不卡在线播放| 国产国语露脸激情在线看| 欧美激情 高清一区二区三区| 五月天丁香电影| 亚洲国产av影院在线观看| 在线精品无人区一区二区三| 亚洲欧美中文字幕日韩二区| 国产成人免费观看mmmm| 欧美日本中文国产一区发布| 97精品久久久久久久久久精品| 亚洲成国产人片在线观看| av有码第一页| 天堂俺去俺来也www色官网| 2018国产大陆天天弄谢| 99国产综合亚洲精品| 久久久久久久大尺度免费视频| 成人影院久久| 国产成人精品无人区|