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

    基于洪水演進(jìn)模擬的潰堤洪災(zāi)損失評(píng)估

    2023-08-02 14:52:56李文歡支歡樂(lè)蔣水華
    人民珠江 2023年7期
    關(guān)鍵詞:模型

    李文歡,支歡樂(lè),蔣水華*,雷 聲,黃 河

    (1.南昌大學(xué)工程建設(shè)學(xué)院,江西 南昌 330031;2.中鐵水利水電規(guī)劃設(shè)計(jì)集團(tuán)有限公司,江西 南昌 330029;3.江西省水利科學(xué)院,江西 南昌 330029)

    洪水作為世界上最常見(jiàn)且最具威脅的自然災(zāi)害之一,嚴(yán)重影響了人們生產(chǎn)生活,制約了中國(guó)經(jīng)濟(jì)社會(huì)的發(fā)展。中國(guó)水旱災(zāi)害公報(bào)資料表明[1],2019年,中國(guó)因洪澇災(zāi)害共造成4 766.6萬(wàn)人受災(zāi),658人死亡失蹤,直接經(jīng)濟(jì)損失1 922.7億元。因此,開(kāi)展洪水災(zāi)害風(fēng)險(xiǎn)評(píng)估,進(jìn)而建立洪水災(zāi)害預(yù)警系統(tǒng),制定有效的防災(zāi)減災(zāi)措施具有重要的工程應(yīng)用價(jià)值。

    目前,眾多學(xué)者在洪水演進(jìn)分析方面開(kāi)展了大量有益的研究[2-4],并發(fā)現(xiàn)相較于Delft3D[2]、HEC-RAS[3]和Delft3D-Flow[4]等軟件,MIKE 21軟件不僅可以考慮風(fēng)速和入滲等多種因素的影響,而且模擬結(jié)果更準(zhǔn)確,計(jì)算效率更高。例如,劉衛(wèi)林等[5]針對(duì)水文資料缺乏的中小河流通過(guò)小流域水文計(jì)算和MIKE 11水動(dòng)力模型進(jìn)行研究,提出了臨界水位預(yù)警指標(biāo)確定方法。在水動(dòng)力學(xué)模型方面,詹明強(qiáng)等[6]研究了潰壩對(duì)下游淹沒(méi)范圍和程度的影響,建立了一維非恒定流水動(dòng)力學(xué)模型。

    洪澇災(zāi)害損失評(píng)估主要分為生命損失、經(jīng)濟(jì)損失和生態(tài)環(huán)境損失這3類(lèi)[7]。對(duì)于生命損失評(píng)估,目前主要參考了國(guó)外經(jīng)驗(yàn)并結(jié)合國(guó)內(nèi)情況估算洪災(zāi)人口死亡率和生命損失[8]。郭磊等[9]通過(guò)夜間燈光的遙感數(shù)據(jù)和土地利用類(lèi)型推求人口空間展布數(shù)據(jù),再疊加歷史洪水模擬得到受災(zāi)人口數(shù)量。經(jīng)濟(jì)損失評(píng)估常分為直接和間接損失評(píng)估,如Pariartha等[10]采用MIKE FLOOD計(jì)算洪水淹沒(méi)水深和淹沒(méi)面積,疊加分析住宅地理信息,再根據(jù)住房損失曲線(xiàn)和平均房?jī)r(jià)來(lái)評(píng)估洪水造成的經(jīng)濟(jì)損失。目前關(guān)于洪水造成的生態(tài)環(huán)境損失研究較為欠缺,這方面的成果不多,國(guó)外則主要以定性方法評(píng)估洪水對(duì)于生態(tài)環(huán)境的影響。近年來(lái)發(fā)展迅猛的地理信息技術(shù)受到了學(xué)者們的青睞,與洪水災(zāi)害損失評(píng)估結(jié)合進(jìn)行洪水淹沒(méi)空間分布數(shù)據(jù)分析和承災(zāi)體空間信息災(zāi)害評(píng)估[8-10]。

    采用MIKE 21軟件模擬潰堤洪水演進(jìn)過(guò)程,根據(jù)土地利用類(lèi)型不同合理確定下墊面糙率值,并基于Python對(duì)ArcGIS進(jìn)行二次開(kāi)發(fā),在此基礎(chǔ)上建立具有明確物理意義的損失定量評(píng)估模型,進(jìn)而量化經(jīng)濟(jì)財(cái)產(chǎn)價(jià)值和生態(tài)系統(tǒng)服務(wù)價(jià)值,結(jié)合潰堤洪水時(shí)空分布數(shù)據(jù),建立洪水災(zāi)害損失評(píng)估模型。最后,以鄱陽(yáng)湖區(qū)某重點(diǎn)實(shí)際工程為例驗(yàn)證方法的有效性,并評(píng)估潰堤洪水造成的生命、經(jīng)濟(jì)和生態(tài)環(huán)境損失。

    1 潰堤洪水演進(jìn)模擬方法

    水動(dòng)力學(xué)模型常被用于模擬洪水過(guò)程演進(jìn)過(guò)程。MIKE 21軟件是目前國(guó)際上流行的水動(dòng)力學(xué)模擬軟件[11],常用于模擬潰堤/潰壩洪水演進(jìn)過(guò)程。Symonds等[2]分析對(duì)比Delft3D和MIKE 21的模擬效果和計(jì)算效率,表明兩者的模擬結(jié)果都與實(shí)測(cè)數(shù)據(jù)吻合良好,但是MIKE 21在多核運(yùn)行時(shí)使用非結(jié)構(gòu)化模型,計(jì)算效率更高。Shrestha等[3]分別采用HEC-RAS和MIKE 21對(duì)比分析了布倫特伍德地區(qū)Deer Creek河流的洪水演進(jìn)過(guò)程,表明兩者都可以精確模擬洪水演進(jìn)過(guò)程,但是相對(duì)于HEC-RAS,MIKE 21還可以考慮風(fēng)速、入滲等因素,更真實(shí)地模擬實(shí)際洪水過(guò)程。Parsapour等[4]分別采用Delft3D-Flow和MIKE 21進(jìn)行水動(dòng)力建模,在同等網(wǎng)格分辨率的情況下,MIKE 21會(huì)獲得更準(zhǔn)確的模擬結(jié)果。MIKE 21采用的二維水流模型基本方程組包括二維流水流連續(xù)方程與水流動(dòng)量方程,其中二維流水流連續(xù)方程見(jiàn)式(1):

    (1)

    式中ξ——水面到基準(zhǔn)面的高度,m;t——計(jì)算時(shí)間,s;u——流速在x方向上的分量,m/s;v——流速在y方向上的分量,m/s。

    x和y方向上的二維水流動(dòng)量方程見(jiàn)式(2)、(3):

    (2)

    (3)

    式中h——靜止水深,m;g——重力加速度,m2/s;V——風(fēng)速;Vx——x方向的風(fēng)速矢量,m/s;Vy——y方向的風(fēng)速矢量,m/s;C——阻力系數(shù);f——風(fēng)摩擦因素,f=γa2ρa(bǔ);γa——風(fēng)應(yīng)力系數(shù);ρa(bǔ)——空氣密度;Ω——科氏力系數(shù),與緯度相關(guān)。

    2 潰堤洪水災(zāi)害損失評(píng)估方法

    基于研究區(qū)域洪災(zāi)承災(zāi)體類(lèi)別、價(jià)值和空間分布情況,構(gòu)建空間人群分布、經(jīng)濟(jì)價(jià)值量化、生態(tài)系統(tǒng)服務(wù)價(jià)值分布數(shù)據(jù)空間化拓?fù)鋽?shù)據(jù)庫(kù),發(fā)展可融合洪水時(shí)空分布數(shù)據(jù)的災(zāi)害損失評(píng)估方法,并結(jié)合洪水淹沒(méi)的時(shí)空分布特征,動(dòng)態(tài)評(píng)估潰堤洪水所造成的生命、經(jīng)濟(jì)和生態(tài)環(huán)境損失。對(duì)于生命損失,一方面需要考慮人群的空間分布數(shù)據(jù),另一方面還需要考慮洪災(zāi)的生命損失率,采用的生命損失LOL計(jì)算見(jiàn)式(4)[12]:

    (4)

    式中 PAR——風(fēng)險(xiǎn)人口數(shù)量,即研究區(qū)域內(nèi)居民受洪水影響的總數(shù)量;a——人群分布數(shù)據(jù)的像元網(wǎng)格總數(shù)量;IRi——第i個(gè)像元網(wǎng)格上的風(fēng)險(xiǎn)個(gè)體生命損失率,計(jì)算見(jiàn)式(5)、(6)[12]:

    IR=f0αβ

    (5)

    α=qm1+(1-q)m2

    (6)

    式中f0——中國(guó)風(fēng)險(xiǎn)人口死亡率;α——洪災(zāi)嚴(yán)重程度系數(shù);β——修正系數(shù),一般取值為1.4;m1——災(zāi)難程度主要影響因子;m2——災(zāi)難程度次要影響因子;q——權(quán)重系數(shù),取值為0.8。

    其中,m1可視為各個(gè)主要影響因素對(duì)風(fēng)險(xiǎn)人口死亡率的影響程度與其權(quán)重系數(shù)的乘積,見(jiàn)式(7)[13]:

    (7)

    式中si——第i個(gè)主要影響因素對(duì)風(fēng)險(xiǎn)人口死亡率的影響程度,主要包含風(fēng)險(xiǎn)人口數(shù)量、潰堤洪水嚴(yán)重性SD、預(yù)警時(shí)間WT以及對(duì)洪水嚴(yán)重性的理解程度UD;θi——權(quán)重系數(shù),分別取值為θ1=θ2=0.2,θ3=θ4=0.3。

    m2為各個(gè)次要影響因素對(duì)風(fēng)險(xiǎn)人口死亡率的影響程度與其權(quán)重系數(shù)的乘積,見(jiàn)式(8)[13]:

    (8)

    式中ti——第i個(gè)次要影響因素對(duì)風(fēng)險(xiǎn)人口死亡率的影響程度;ni——ti對(duì)應(yīng)的權(quán)重系數(shù)值。

    經(jīng)濟(jì)損失評(píng)估主要包括直接經(jīng)濟(jì)損失和間接經(jīng)濟(jì)損失兩方面。其中,直接經(jīng)濟(jì)損失采取損失率方法計(jì)算,通過(guò)建立各類(lèi)財(cái)產(chǎn)洪水損失率,結(jié)合洪水淹沒(méi)特征、經(jīng)濟(jì)損失統(tǒng)計(jì)指標(biāo)以及淹沒(méi)水深等級(jí),進(jìn)行疊加計(jì)算見(jiàn)式(9)[14]:

    (9)

    式中D——總直接經(jīng)濟(jì)損失;Sp——第p類(lèi)財(cái)產(chǎn)的直接經(jīng)濟(jì)損失;e——研究區(qū)域財(cái)產(chǎn)分布在二維空間上的像元網(wǎng)格數(shù)量;h——第j個(gè)像元網(wǎng)格上的財(cái)產(chǎn)種類(lèi);l——不同淹沒(méi)水深等級(jí);βijk——第i個(gè)像元網(wǎng)格上第j類(lèi)財(cái)產(chǎn)在第k種淹沒(méi)水深對(duì)應(yīng)的損失率;Wijk——第i個(gè)像元網(wǎng)格上第j類(lèi)財(cái)產(chǎn)對(duì)應(yīng)的財(cái)產(chǎn)價(jià)值。

    目前間接經(jīng)濟(jì)損失R主要通過(guò)洪災(zāi)間接系數(shù)進(jìn)行確定,根據(jù)直接經(jīng)濟(jì)損失與間接經(jīng)濟(jì)損失的關(guān)系計(jì)算見(jiàn)式(10)[14]:

    (10)

    式中Ki——第p類(lèi)財(cái)產(chǎn)對(duì)應(yīng)的間接系數(shù),農(nóng)業(yè)取值為15%~28%,住宅區(qū)取值為15%,公路取值為25%。

    對(duì)于生態(tài)環(huán)境損失計(jì)算,改進(jìn)了文獻(xiàn)[15]方法,見(jiàn)式(11):

    E=c×ESV

    (11)

    式中c——經(jīng)驗(yàn)損失率,取c為0.1;ESV——被洪水淹沒(méi)區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值;E——按經(jīng)驗(yàn)系數(shù)法計(jì)算的生態(tài)環(huán)境損失。

    3 洪災(zāi)損失評(píng)估系統(tǒng)實(shí)現(xiàn)

    為定量評(píng)估潰堤洪水災(zāi)害損失,借助面向?qū)ο蟮挠?jì)算機(jī)程序語(yǔ)言——Python,調(diào)用Arcpy站點(diǎn)包實(shí)現(xiàn)內(nèi)嵌于ArcGIS的內(nèi)置插件條,基于堤防保護(hù)區(qū)域人群空間分布數(shù)據(jù)、經(jīng)濟(jì)財(cái)產(chǎn)價(jià)值量化空間分布數(shù)據(jù)和生態(tài)系統(tǒng)服務(wù)價(jià)值空間分布數(shù)據(jù),有機(jī)結(jié)合洪水時(shí)空分布特征,建立洪水災(zāi)害損失動(dòng)態(tài)評(píng)估模型,實(shí)現(xiàn)對(duì)生命、經(jīng)濟(jì)以及生態(tài)環(huán)境損失的動(dòng)態(tài)定量評(píng)估。

    在損失評(píng)估模型的實(shí)現(xiàn)過(guò)程中,洪水分布數(shù)據(jù)和承載體空間分布數(shù)據(jù)的類(lèi)型及精度是不一致的。在后續(xù)的損失計(jì)算過(guò)程中,要求數(shù)據(jù)可以相互疊加計(jì)算,需要保證洪水?dāng)?shù)據(jù)和承載體數(shù)據(jù)的投影坐標(biāo)系、柵格像元大小以及數(shù)據(jù)類(lèi)型等一致,所以需要提前進(jìn)行數(shù)據(jù)柵格標(biāo)準(zhǔn)化處理。具體步驟如下:①根據(jù)潰堤洪水演進(jìn)模型,提取未被標(biāo)準(zhǔn)化的洪水?dāng)?shù)據(jù);②將洪水?dāng)?shù)據(jù)進(jìn)行面轉(zhuǎn)柵格處理,然后將柵格化的洪水?dāng)?shù)據(jù)與承載體數(shù)據(jù)空間分布數(shù)據(jù)進(jìn)一步進(jìn)行標(biāo)準(zhǔn)化處理。

    對(duì)于生命損失評(píng)估模塊,在數(shù)據(jù)標(biāo)準(zhǔn)化處理模塊的基礎(chǔ)上,實(shí)現(xiàn)動(dòng)態(tài)生命損失評(píng)估的具體步驟如下:①輸入洪水的時(shí)間尺度,本文輸入時(shí)間是潰堤洪水演進(jìn)過(guò)程中的模型起始時(shí)間和末尾時(shí)間,設(shè)定的間隔時(shí)間為2 h,生命損失模塊可根據(jù)間隔時(shí)間循環(huán)讀取淹沒(méi)水深數(shù)據(jù)和流速數(shù)據(jù);②設(shè)定模型參數(shù),計(jì)算個(gè)體風(fēng)險(xiǎn)生命損失率IR分布情況,累加計(jì)算所有像元網(wǎng)格的個(gè)體風(fēng)險(xiǎn)生命損失率總值,疊加計(jì)算人群數(shù)據(jù),將獲得的生命損失值統(tǒng)計(jì)情況存入表格中。

    對(duì)于經(jīng)濟(jì)損失評(píng)估模塊,洪水淹沒(méi)是1個(gè)動(dòng)態(tài)變化的過(guò)程,某一瞬時(shí)的洪水?dāng)?shù)據(jù)僅能分析這個(gè)時(shí)間點(diǎn)的經(jīng)濟(jì)損失空間分布情況。經(jīng)濟(jì)損失評(píng)估具體步驟如下:①為實(shí)現(xiàn)經(jīng)濟(jì)損失動(dòng)態(tài)計(jì)算,輸入潰堤洪水演進(jìn)過(guò)程中起始、末尾時(shí)間以及間隔時(shí)間,表示每過(guò)一次間隔時(shí)間計(jì)算一次經(jīng)濟(jì)損失;②輸入經(jīng)濟(jì)分布數(shù)據(jù)和淹沒(méi)水深數(shù)據(jù),而財(cái)產(chǎn)種類(lèi)數(shù)據(jù)和損失率關(guān)系可根據(jù)蓄滯洪區(qū)或者防洪保護(hù)區(qū)的實(shí)際社會(huì)情況進(jìn)行修改;③輸出各類(lèi)財(cái)產(chǎn)損失的統(tǒng)計(jì)表以及經(jīng)濟(jì)損失的分布數(shù)據(jù)。

    對(duì)于生態(tài)環(huán)境損失評(píng)估模塊,根據(jù)式(11)可知,需要輸入生態(tài)系統(tǒng)服務(wù)價(jià)值空間分布數(shù)據(jù)。生態(tài)環(huán)境損失評(píng)估具體步驟如下:①輸入始末以及間隔時(shí)間;②通過(guò)疊加洪水淹沒(méi)的面積篩選受影響的生態(tài)系統(tǒng)服務(wù)價(jià)值空間分布數(shù)據(jù),設(shè)定經(jīng)驗(yàn)損失率,計(jì)算生態(tài)損失;③輸出生態(tài)損失統(tǒng)計(jì)表。

    本文建立的洪災(zāi)損失動(dòng)態(tài)評(píng)估模型有效利用了潰堤洪水演進(jìn)模擬獲取的洪水空間分布數(shù)據(jù),其主要優(yōu)勢(shì)如下:①基于ArcGIS平臺(tái)構(gòu)建空間化拓?fù)鋽?shù)據(jù)庫(kù),并且與國(guó)家數(shù)據(jù)統(tǒng)計(jì)年鑒內(nèi)容緊密結(jié)合,保證了數(shù)據(jù)有效性,同時(shí)解決了數(shù)據(jù)的及時(shí)更新問(wèn)題;②考慮了洪水和承載體數(shù)據(jù)的空間分布差異特征,有效提升了洪水損失評(píng)估結(jié)果準(zhǔn)確性;③通過(guò)有限的輸入和輸出以及模型參數(shù)修改反映了洪災(zāi)損失動(dòng)態(tài)變化情況;④可以根據(jù)研究對(duì)象的洪災(zāi)損失典型資料,調(diào)查分析和調(diào)整洪災(zāi)損失率關(guān)系,兼顧了模型的普適性;⑤實(shí)現(xiàn)了潰堤洪水演進(jìn)模擬、空間化拓?fù)鋽?shù)據(jù)庫(kù)動(dòng)態(tài)連接以及洪災(zāi)損失的實(shí)時(shí)評(píng)估。

    4 工程應(yīng)用

    鄱陽(yáng)湖區(qū)某重點(diǎn)堤防地處江西省九江市永修縣境內(nèi),圩堤起于永修縣城郊的觀音巖,止于楊公腦,全長(zhǎng)33.57 km,保護(hù)面積56.28 km2,保護(hù)耕地5.03×104畝(1畝約等于666.67 m2,下同),常住人口2.34萬(wàn)人。據(jù)2005年現(xiàn)場(chǎng)勘察,三角聯(lián)圩各類(lèi)險(xiǎn)工險(xiǎn)段多達(dá)6處,累計(jì)堤線(xiàn)長(zhǎng)30多km,大多堤段的斷面較為單薄且內(nèi)外坡均較陡,汛期水流逐年沖刷,極易導(dǎo)致堤身滲透破壞,存在部分堤段填筑壓實(shí)不均勻且高程不足,脫坡險(xiǎn)情突出。2020年7月12日19時(shí)40分,潰口寬度達(dá)20余m,潰堤后,圩內(nèi)平均水位在0.5 h內(nèi)上漲0.5 m多,由于洪水沖刷,潰口迅速擴(kuò)大,截至7月13日10點(diǎn),三角聯(lián)圩潰口寬度已經(jīng)擴(kuò)大到200余m,并于7月14日14時(shí)啟動(dòng)封堵作業(yè),7月16日21時(shí)43分,三角聯(lián)圩完成合攏。淹沒(méi)耕地面積5.54×104畝,水產(chǎn)業(yè)面積1.05×104畝,淹沒(méi)5 100余棟房屋,其中倒塌了40余棟,潰口經(jīng)過(guò)55 h的搶修后合攏。

    4.1 模型計(jì)算范圍確定

    基于中國(guó)大地坐標(biāo)系,建立堤防保護(hù)區(qū)域數(shù)字高程模型見(jiàn)圖1。采用三角形網(wǎng)格進(jìn)行模型剖分見(jiàn)圖2,將東南部堤防潰口(紅色邊界)設(shè)置為入流邊界,北部三角聯(lián)圩和西部陸地(綠色邊界)設(shè)置為閉邊界。

    圖1 堤防保護(hù)區(qū)域數(shù)字高程模型

    圖2 堤防水動(dòng)力學(xué)模型網(wǎng)格剖分

    本數(shù)值模型共存在4個(gè)片區(qū),分別為堤防南部片區(qū)、北部片區(qū)、東部片區(qū)和西部片區(qū),在每個(gè)片區(qū)內(nèi)分別設(shè)置4個(gè)特征點(diǎn)(即自然村)。假設(shè)決口發(fā)生在三角聯(lián)圩南部片區(qū)下方見(jiàn)圖3,其具體地理位置為九江市永修縣三角鄉(xiāng)、永豐墾殖場(chǎng)和南昌市新建區(qū)大塘坪鄉(xiāng)交界,樁號(hào)27k+880~28k+010 m處。

    圖3 三角堤防區(qū)域特征點(diǎn)和潰口位置

    4.2 模型參數(shù)設(shè)置及有效性

    提出的洪水演進(jìn)模擬方法涉及的物理參數(shù)(包括邊界條件、網(wǎng)格形式、初始條件、渦黏系數(shù))見(jiàn)表1。

    表1 洪水演進(jìn)數(shù)值模型主要物理參數(shù)取值[16]

    本模型中土地利用類(lèi)型主要包括旱地、水田、草地、灘地、居民區(qū)和水面6種,其相應(yīng)糙率參數(shù)按照下墊面分布情況,并參照相關(guān)經(jīng)驗(yàn)進(jìn)行率定選取[19],見(jiàn)表2。為了驗(yàn)證模型的有效性,圖4比較了采用所建模型模擬的潰口淹沒(méi)水位過(guò)程線(xiàn)和實(shí)測(cè)的潰口區(qū)域內(nèi)淹沒(méi)水位過(guò)程線(xiàn)。由圖4可知,在水位基本不變后,水位模擬值為22.65 m,水位實(shí)測(cè)值為22.58 m,兩條水位過(guò)程線(xiàn)較為吻合,表明基于MIKE 21軟件能夠合理模擬潰堤洪水在堤防保護(hù)區(qū)域內(nèi)的演進(jìn)過(guò)程。此外,所建模型模擬的淹沒(méi)面積為56.15 km2,與堤防潰口合攏后的實(shí)測(cè)淹沒(méi)面積(56.28 km2)相差僅0.13 km2,進(jìn)一步說(shuō)明了本模型的有效性,可以滿(mǎn)足計(jì)算精度要求。

    表2 糙率參數(shù)參考取值

    圖4 潰口區(qū)域內(nèi)模擬水位與實(shí)測(cè)水位的比較

    4.3 潰堤洪水演進(jìn)過(guò)程分析

    圖5給出了三角聯(lián)圩潰口流量隨潰決時(shí)間的變化情況。2020年7月12日19時(shí)為圖5中的起始時(shí)刻,在7月12日19時(shí)40分堤防潰決,隨后堤防潰口寬度逐漸擴(kuò)大至最大值200 m,潰堤后的21 h內(nèi)流量從0急劇增加到最大值1 724 m3/s,隨后流量開(kāi)始下降,于7月14日14時(shí)啟動(dòng)封堵作業(yè),7月16日21時(shí)43分,三角聯(lián)圩潰口完成合攏,流量接近于零,整個(gè)堤防保護(hù)區(qū)域進(jìn)洪過(guò)程在98 h內(nèi)完成。

    圖5 三角堤防潰口流量過(guò)程線(xiàn)

    為了進(jìn)一步展示潰堤洪水演進(jìn)的淹沒(méi)過(guò)程,圖6給出了6個(gè)典型進(jìn)洪時(shí)段的淹沒(méi)水深分布。洪水在決堤后10 h內(nèi)到達(dá)堤防南部片區(qū),雖然紅旗村處于潰口附近,但是因紅旗村的地勢(shì)較高,在潰決40 h內(nèi)受災(zāi)情況較輕。由于堤防北部片區(qū)地勢(shì)較低,洪水在20 h內(nèi)向周家湖、三角圩等村莊挺進(jìn),潰堤40 h后,堤防北部片區(qū)基本被淹沒(méi),洪水逐漸向東西部蔓延,堤防保護(hù)區(qū)域內(nèi)水深高達(dá)5.2 m。潰堤60 h后,堤防北部片區(qū)淹沒(méi)水深最高達(dá)6.0 m;潰堤80 h后,因堤防西部片區(qū)比東部片區(qū)地勢(shì)較高,故東部區(qū)域淹沒(méi)深度最深達(dá)到6.5 m,而西部區(qū)域仍舊有部分地區(qū)未被淹沒(méi)。進(jìn)洪100 h后,堤防保護(hù)區(qū)域絕大部分被淹沒(méi),此時(shí)需要盡快對(duì)圩內(nèi)居民實(shí)施避險(xiǎn)轉(zhuǎn)移搶險(xiǎn)救援等舉措。

    a)10 h

    4.4 潰堤洪水損失評(píng)估

    根據(jù)收集的洪水影響對(duì)象統(tǒng)計(jì)數(shù)據(jù),對(duì)人群數(shù)據(jù)進(jìn)行空間展布,并量化經(jīng)濟(jì)財(cái)產(chǎn)和生態(tài)系統(tǒng)服務(wù)價(jià)值。然后,搭建空間化拓?fù)鋽?shù)據(jù)庫(kù)來(lái)反映生命、經(jīng)濟(jì)、生態(tài)環(huán)境在空間上的分布差異,接著在潰堤洪水演進(jìn)模型中提取洪水淹沒(méi)時(shí)空分布數(shù)據(jù),并將數(shù)據(jù)柵格進(jìn)行標(biāo)準(zhǔn)化處理。最后,進(jìn)行包括生命、經(jīng)濟(jì)和生態(tài)環(huán)境損失在內(nèi)的綜合損失動(dòng)態(tài)評(píng)估。圖7給出了潰堤造成的生命損失隨潰口歷時(shí)的變化曲線(xiàn)。由圖7可知,隨著預(yù)警時(shí)間的延長(zhǎng),生命損失不斷減少。同時(shí),居民對(duì)洪水嚴(yán)重性的理解程度也對(duì)生命損失具有重要的影響因素。例如,在預(yù)警時(shí)間不超過(guò)0.25 h且居民對(duì)洪水嚴(yán)重性理解程度為模糊的情況下,生命損失高達(dá)454人,見(jiàn)圖7a;相較之下,如果預(yù)警時(shí)間足夠長(zhǎng),且居民對(duì)洪水嚴(yán)重性的理解程度為明確的情況下,生命損失驟降至2人,見(jiàn)圖7b。

    a)居民模糊理解洪水嚴(yán)重性情況下生命損失曲線(xiàn)

    圖8給出了潰堤造成的經(jīng)濟(jì)損失隨潰口歷時(shí)的變化關(guān)系曲線(xiàn)。由圖8可知,直接經(jīng)濟(jì)損失隨著潰口歷時(shí)的增加而增大,其中居民住房直接經(jīng)濟(jì)損失最大,高達(dá)2.86億元,其次為家庭財(cái)產(chǎn)直接經(jīng)濟(jì)損失。同時(shí),圖9給出了堤防保護(hù)區(qū)域內(nèi)總經(jīng)濟(jì)損失分布情況,不難看出居民用地?fù)p失情況最為嚴(yán)重。

    a)耕地-家庭財(cái)產(chǎn)-住房直接經(jīng)濟(jì)損失

    圖9 潰堤造成的總經(jīng)濟(jì)損失分布

    5 結(jié)論

    以江西省鄱陽(yáng)湖區(qū)某重點(diǎn)堤防為例,建立了基于MIKE 21軟件的洪水演進(jìn)數(shù)值模型,獲得了潰堤洪水下淹沒(méi)水深和流速等洪水空間分布數(shù)據(jù),豐富了潰堤洪水演進(jìn)模擬方法。同時(shí),考慮生命、經(jīng)濟(jì)、生態(tài)環(huán)境的空間分布特性,基于Python平臺(tái)建立了一種具有明確物理意義的,可考慮生命、經(jīng)濟(jì)和環(huán)境損失的定量動(dòng)態(tài)評(píng)估模型。主要結(jié)論如下。

    a)根據(jù)研究區(qū)下墊面土地利用情況設(shè)置不同的糙率值,以及與實(shí)測(cè)的潰口處淹沒(méi)水位過(guò)程線(xiàn)進(jìn)行對(duì)比,驗(yàn)證了洪水演進(jìn)數(shù)值模型的有效性,表明所建模型具有較高的計(jì)算精度,可有效用于獲取淹沒(méi)水深、流速等洪水空間分布數(shù)據(jù)。

    b)發(fā)展了潰堤洪水災(zāi)害損失評(píng)估方法,并采用Python語(yǔ)言編寫(xiě)了基于ArcGIS的嵌入式耦合插件工具,建立了洪災(zāi)損失動(dòng)態(tài)評(píng)估模型,可以有效評(píng)估不同工況下潰堤洪水造成的生命、經(jīng)濟(jì)以及生態(tài)環(huán)境損失,為洪水災(zāi)害防控提供了有效途徑。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線(xiàn)三等角』
    重尾非線(xiàn)性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    99热精品在线国产| 国产熟女xx| 色综合欧美亚洲国产小说| 给我免费播放毛片高清在线观看| 欧美成人性av电影在线观看| 国产成人福利小说| 国产精品av视频在线免费观看| 精品久久久久久,| 国产99白浆流出| 婷婷亚洲欧美| 又黄又粗又硬又大视频| 啦啦啦观看免费观看视频高清| 中文在线观看免费www的网站| 99热这里只有精品一区| 熟女少妇亚洲综合色aaa.| 狂野欧美激情性xxxx| 狂野欧美白嫩少妇大欣赏| 午夜精品在线福利| 午夜福利视频1000在线观看| 国产成人影院久久av| 一个人观看的视频www高清免费观看| 亚洲国产精品成人综合色| 超碰av人人做人人爽久久 | 波多野结衣高清无吗| ponron亚洲| 色哟哟哟哟哟哟| 免费在线观看影片大全网站| 色哟哟哟哟哟哟| 亚洲电影在线观看av| 欧美一区二区精品小视频在线| 欧美最新免费一区二区三区 | 欧美3d第一页| 一本久久中文字幕| 一个人观看的视频www高清免费观看| av在线天堂中文字幕| 国产精品日韩av在线免费观看| 97超视频在线观看视频| 精品国产三级普通话版| 3wmmmm亚洲av在线观看| eeuss影院久久| 国产不卡一卡二| 男女午夜视频在线观看| 欧美一区二区精品小视频在线| 国产免费一级a男人的天堂| 1000部很黄的大片| 婷婷亚洲欧美| 欧美日韩瑟瑟在线播放| 9191精品国产免费久久| 少妇的逼水好多| 中文字幕av在线有码专区| 日韩欧美在线乱码| 韩国av一区二区三区四区| 黄色成人免费大全| 12—13女人毛片做爰片一| 午夜福利在线在线| 久久亚洲精品不卡| 久久久久久久久大av| 草草在线视频免费看| 少妇丰满av| 九九久久精品国产亚洲av麻豆| 欧美在线黄色| 亚洲片人在线观看| 婷婷六月久久综合丁香| 欧美在线黄色| 18禁黄网站禁片免费观看直播| 欧美又色又爽又黄视频| 精品人妻1区二区| 美女大奶头视频| 午夜福利18| 日韩成人在线观看一区二区三区| 欧美av亚洲av综合av国产av| 在线视频色国产色| 男女午夜视频在线观看| 国产男靠女视频免费网站| 人人妻,人人澡人人爽秒播| 成人国产一区最新在线观看| 男人舔奶头视频| 国产亚洲精品av在线| 午夜免费男女啪啪视频观看 | 成人精品一区二区免费| 精品国内亚洲2022精品成人| 欧美一级a爱片免费观看看| 国产av一区在线观看免费| 亚洲男人的天堂狠狠| 99久久成人亚洲精品观看| 欧美中文综合在线视频| 高清在线国产一区| 最新美女视频免费是黄的| 日韩大尺度精品在线看网址| 亚洲中文日韩欧美视频| 脱女人内裤的视频| 色综合婷婷激情| 亚洲va日本ⅴa欧美va伊人久久| 免费无遮挡裸体视频| 精品不卡国产一区二区三区| 国产精品亚洲一级av第二区| 熟女电影av网| av在线蜜桃| 日韩欧美 国产精品| 91字幕亚洲| 五月伊人婷婷丁香| 亚洲av免费在线观看| 男女床上黄色一级片免费看| 久久久久久久午夜电影| 欧美bdsm另类| 他把我摸到了高潮在线观看| 天堂√8在线中文| 淫秽高清视频在线观看| 叶爱在线成人免费视频播放| 午夜亚洲福利在线播放| 99久久无色码亚洲精品果冻| 欧美成人一区二区免费高清观看| a在线观看视频网站| 国产真实乱freesex| 好男人在线观看高清免费视频| 午夜精品久久久久久毛片777| 亚洲精品美女久久久久99蜜臀| 亚洲黑人精品在线| 亚洲人成网站在线播| 在线a可以看的网站| 真人一进一出gif抽搐免费| 最近在线观看免费完整版| 18禁裸乳无遮挡免费网站照片| 国产免费av片在线观看野外av| 极品教师在线免费播放| 精品久久久久久久末码| 欧美色欧美亚洲另类二区| 日韩欧美国产在线观看| 俄罗斯特黄特色一大片| 老汉色av国产亚洲站长工具| 99国产精品一区二区蜜桃av| tocl精华| 国产精品久久久久久久久免 | 亚洲国产日韩欧美精品在线观看 | 观看美女的网站| 成年女人永久免费观看视频| 午夜福利在线观看吧| 成人性生交大片免费视频hd| 欧美另类亚洲清纯唯美| 中文字幕久久专区| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品sss在线观看| 国产一区二区亚洲精品在线观看| 毛片女人毛片| 在线播放国产精品三级| 欧美+日韩+精品| av中文乱码字幕在线| 变态另类丝袜制服| 在线观看午夜福利视频| 18禁黄网站禁片午夜丰满| 国产在线精品亚洲第一网站| 少妇熟女aⅴ在线视频| 1000部很黄的大片| 亚洲激情在线av| 欧美激情在线99| 99热这里只有精品一区| 日本熟妇午夜| 精华霜和精华液先用哪个| 久久久久久国产a免费观看| 国产精品嫩草影院av在线观看 | 欧美av亚洲av综合av国产av| 九色成人免费人妻av| www.www免费av| 成人鲁丝片一二三区免费| 成人一区二区视频在线观看| 国产亚洲精品久久久久久毛片| 3wmmmm亚洲av在线观看| 国产一区二区亚洲精品在线观看| 国产精华一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 在线国产一区二区在线| 国产三级中文精品| 亚洲熟妇熟女久久| 在线观看日韩欧美| 日本一本二区三区精品| 动漫黄色视频在线观看| 99久久精品热视频| 一夜夜www| 又爽又黄无遮挡网站| 欧美+日韩+精品| 此物有八面人人有两片| 老汉色av国产亚洲站长工具| 成人无遮挡网站| 亚洲国产精品久久男人天堂| 中文资源天堂在线| 1000部很黄的大片| 久久婷婷人人爽人人干人人爱| 麻豆国产97在线/欧美| 中亚洲国语对白在线视频| 亚洲不卡免费看| 国产成+人综合+亚洲专区| 精品人妻1区二区| 久久久国产精品麻豆| 国产亚洲精品久久久com| 99久久成人亚洲精品观看| 国产精品乱码一区二三区的特点| 精品久久久久久,| 天堂√8在线中文| 久久精品国产亚洲av香蕉五月| 男女做爰动态图高潮gif福利片| 久久久精品欧美日韩精品| 性色avwww在线观看| 久久久久久九九精品二区国产| h日本视频在线播放| av福利片在线观看| 不卡一级毛片| 久久这里只有精品中国| 国产精品久久久久久精品电影| 日本精品一区二区三区蜜桃| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久久久电影 | 嫩草影院精品99| 1024手机看黄色片| 久久久精品大字幕| 国产视频一区二区在线看| 热99在线观看视频| 成人午夜高清在线视频| www国产在线视频色| 欧美日韩国产亚洲二区| 国产精品影院久久| 亚洲av日韩精品久久久久久密| 久9热在线精品视频| 日本一本二区三区精品| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 亚洲精华国产精华精| 欧美另类亚洲清纯唯美| 亚洲最大成人中文| 亚洲自拍偷在线| 国产一区二区三区在线臀色熟女| 99在线视频只有这里精品首页| 97人妻精品一区二区三区麻豆| 欧美大码av| 午夜福利免费观看在线| 夜夜爽天天搞| 一级黄片播放器| 淫妇啪啪啪对白视频| 成人欧美大片| 久久国产精品影院| 色噜噜av男人的天堂激情| 18禁美女被吸乳视频| 桃色一区二区三区在线观看| av专区在线播放| 大型黄色视频在线免费观看| 亚洲精品国产精品久久久不卡| 欧美成人一区二区免费高清观看| 欧美日韩亚洲国产一区二区在线观看| 免费人成在线观看视频色| 不卡一级毛片| 岛国视频午夜一区免费看| 日本一二三区视频观看| 亚洲avbb在线观看| 日本熟妇午夜| 欧美在线黄色| 欧美成人a在线观看| 成熟少妇高潮喷水视频| 一个人看视频在线观看www免费 | 啪啪无遮挡十八禁网站| 亚洲av不卡在线观看| 亚洲欧美日韩东京热| 国产三级在线视频| 噜噜噜噜噜久久久久久91| 狂野欧美白嫩少妇大欣赏| 中文字幕高清在线视频| 亚洲精品亚洲一区二区| or卡值多少钱| 久久久久性生活片| 亚洲成av人片在线播放无| av中文乱码字幕在线| 黄色视频,在线免费观看| 久久久久性生活片| 国产熟女xx| 一级a爱片免费观看的视频| 免费看美女性在线毛片视频| 两人在一起打扑克的视频| 欧美高清成人免费视频www| 免费看a级黄色片| 免费大片18禁| 成人亚洲精品av一区二区| h日本视频在线播放| 欧美午夜高清在线| 国产成人影院久久av| 观看免费一级毛片| 国产99白浆流出| 有码 亚洲区| 欧美日韩福利视频一区二区| 午夜久久久久精精品| 一级黄色大片毛片| 一进一出好大好爽视频| 男人的好看免费观看在线视频| 国产精品亚洲av一区麻豆| 在线十欧美十亚洲十日本专区| 日韩欧美精品v在线| 免费大片18禁| 国产精品久久电影中文字幕| 黄色日韩在线| 黄片小视频在线播放| 最近在线观看免费完整版| 观看免费一级毛片| 亚洲专区中文字幕在线| 中文字幕熟女人妻在线| 久久草成人影院| 国产黄片美女视频| 日韩欧美一区二区三区在线观看| 久久久久国内视频| 看免费av毛片| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 91麻豆av在线| 婷婷精品国产亚洲av| 欧美黑人欧美精品刺激| 成人性生交大片免费视频hd| 国产伦精品一区二区三区四那| 男女午夜视频在线观看| 51午夜福利影视在线观看| 欧洲精品卡2卡3卡4卡5卡区| 少妇高潮的动态图| 国产精品久久久久久精品电影| 美女cb高潮喷水在线观看| 人妻丰满熟妇av一区二区三区| 人妻夜夜爽99麻豆av| 欧美性猛交黑人性爽| 狂野欧美激情性xxxx| 亚洲精品一卡2卡三卡4卡5卡| 国产伦精品一区二区三区视频9 | 国产真实乱freesex| a级一级毛片免费在线观看| 午夜免费成人在线视频| 欧美一区二区亚洲| 国产精品久久电影中文字幕| 国产视频一区二区在线看| 欧美黄色片欧美黄色片| 久久午夜亚洲精品久久| a级一级毛片免费在线观看| 欧美日本亚洲视频在线播放| а√天堂www在线а√下载| 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 最后的刺客免费高清国语| 特大巨黑吊av在线直播| 欧美不卡视频在线免费观看| 久久精品国产清高在天天线| 日韩亚洲欧美综合| 精品国内亚洲2022精品成人| 在线观看av片永久免费下载| av在线天堂中文字幕| 美女高潮的动态| 国产精品一区二区三区四区免费观看 | av天堂中文字幕网| 日本与韩国留学比较| 久久6这里有精品| 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 69人妻影院| 精品久久久久久久久久免费视频| 两个人的视频大全免费| 日本 av在线| 亚洲av免费高清在线观看| 国产成+人综合+亚洲专区| 久久久久久久久大av| 狂野欧美白嫩少妇大欣赏| 内射极品少妇av片p| 国产免费av片在线观看野外av| 久久久国产成人免费| 9191精品国产免费久久| 在线视频色国产色| 99国产精品一区二区蜜桃av| 亚洲在线自拍视频| 午夜福利在线在线| 亚洲黑人精品在线| 日韩大尺度精品在线看网址| 日韩欧美 国产精品| 观看美女的网站| 最新美女视频免费是黄的| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩卡通动漫| 最近最新中文字幕大全电影3| 亚洲 国产 在线| 亚洲欧美一区二区三区黑人| xxx96com| 97人妻精品一区二区三区麻豆| 亚洲国产色片| 免费无遮挡裸体视频| 亚洲国产欧洲综合997久久,| 成年女人毛片免费观看观看9| 1000部很黄的大片| 免费观看人在逋| 最近最新中文字幕大全电影3| 国产乱人视频| 亚洲黑人精品在线| 国产激情欧美一区二区| 亚洲成人中文字幕在线播放| www.www免费av| 午夜免费男女啪啪视频观看 | 国产午夜精品久久久久久一区二区三区 | 午夜两性在线视频| 国产精品,欧美在线| 国产97色在线日韩免费| 亚洲专区国产一区二区| 国产av不卡久久| 国产97色在线日韩免费| 国产精品亚洲一级av第二区| 久久精品综合一区二区三区| 久久精品国产99精品国产亚洲性色| 久久国产精品影院| 亚洲av日韩精品久久久久久密| 亚洲男人的天堂狠狠| 精品一区二区三区视频在线观看免费| 九九热线精品视视频播放| 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 久久久精品欧美日韩精品| 国产在视频线在精品| 欧美日本亚洲视频在线播放| 51午夜福利影视在线观看| 国产 一区 欧美 日韩| 国产三级中文精品| 国产精品久久视频播放| 在线国产一区二区在线| 欧美午夜高清在线| 热99re8久久精品国产| 国产精品99久久99久久久不卡| 国产av不卡久久| 亚洲avbb在线观看| 欧美一区二区精品小视频在线| 免费大片18禁| 人妻丰满熟妇av一区二区三区| 男女下面进入的视频免费午夜| 青草久久国产| 特级一级黄色大片| 亚洲av免费在线观看| 免费看日本二区| 日韩欧美 国产精品| 精品久久久久久久末码| 99久久精品一区二区三区| 少妇人妻一区二区三区视频| a在线观看视频网站| 啦啦啦韩国在线观看视频| 亚洲精品亚洲一区二区| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 波多野结衣高清作品| 99热这里只有精品一区| 久久久久亚洲av毛片大全| 午夜免费激情av| 人人妻人人看人人澡| 日本与韩国留学比较| 人人妻人人澡欧美一区二区| 国产精品久久视频播放| 国产毛片a区久久久久| 国产精品三级大全| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 亚洲国产色片| 精品电影一区二区在线| 麻豆国产av国片精品| 老汉色∧v一级毛片| www.www免费av| h日本视频在线播放| 色播亚洲综合网| 国产欧美日韩精品一区二区| 久久久久亚洲av毛片大全| 国产av不卡久久| 91麻豆av在线| 男人舔奶头视频| 男人和女人高潮做爰伦理| 老司机在亚洲福利影院| 欧美性感艳星| 国产激情欧美一区二区| 宅男免费午夜| 欧美成人免费av一区二区三区| 久久久久久九九精品二区国产| 十八禁人妻一区二区| 香蕉久久夜色| 国产一区二区在线av高清观看| 尤物成人国产欧美一区二区三区| 啦啦啦免费观看视频1| 午夜激情欧美在线| 日本黄色视频三级网站网址| 99国产精品一区二区蜜桃av| 久久精品国产综合久久久| 日韩欧美三级三区| 欧美国产日韩亚洲一区| 日本与韩国留学比较| 久9热在线精品视频| 国产精品亚洲一级av第二区| 热99re8久久精品国产| 老熟妇乱子伦视频在线观看| 噜噜噜噜噜久久久久久91| 国产成人a区在线观看| 久久婷婷人人爽人人干人人爱| 99久国产av精品| 午夜福利18| 操出白浆在线播放| 999久久久精品免费观看国产| 夜夜爽天天搞| 国产精品一区二区三区四区免费观看 | 少妇熟女aⅴ在线视频| 国产av不卡久久| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 黄色丝袜av网址大全| 欧美在线黄色| 欧美黄色淫秽网站| 此物有八面人人有两片| 久99久视频精品免费| 国内精品美女久久久久久| 日韩高清综合在线| 亚洲av五月六月丁香网| 欧美一区二区亚洲| 久久精品国产综合久久久| 综合色av麻豆| 亚洲中文日韩欧美视频| 床上黄色一级片| 免费观看精品视频网站| 亚洲av免费在线观看| 国产 一区 欧美 日韩| 欧美午夜高清在线| svipshipincom国产片| 中文字幕人成人乱码亚洲影| 国产97色在线日韩免费| 久久欧美精品欧美久久欧美| 麻豆成人午夜福利视频| 精品一区二区三区视频在线观看免费| 一个人观看的视频www高清免费观看| 欧美日本视频| 国产精品爽爽va在线观看网站| 99久久精品一区二区三区| 长腿黑丝高跟| 两个人看的免费小视频| 黄色女人牲交| 亚洲国产精品久久男人天堂| 精品国产三级普通话版| 熟女少妇亚洲综合色aaa.| 亚洲国产精品sss在线观看| 97超级碰碰碰精品色视频在线观看| 国产乱人视频| 美女 人体艺术 gogo| 婷婷精品国产亚洲av在线| 可以在线观看的亚洲视频| 日本免费一区二区三区高清不卡| 99精品久久久久人妻精品| 国产爱豆传媒在线观看| 国产真实乱freesex| 99久久久亚洲精品蜜臀av| 欧美成狂野欧美在线观看| 日本成人三级电影网站| 久久久久国内视频| 国产老妇女一区| 最近最新免费中文字幕在线| 三级男女做爰猛烈吃奶摸视频| 免费看十八禁软件| 一级黄色大片毛片| a在线观看视频网站| 国产极品精品免费视频能看的| 免费观看的影片在线观看| 全区人妻精品视频| 男人和女人高潮做爰伦理| 身体一侧抽搐| 亚洲最大成人手机在线| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 亚洲中文字幕一区二区三区有码在线看| 给我免费播放毛片高清在线观看| 精品国产超薄肉色丝袜足j| 男人的好看免费观看在线视频| 国产亚洲精品久久久久久毛片| 女人十人毛片免费观看3o分钟| 成人国产综合亚洲| 男女视频在线观看网站免费| 精品福利观看| 欧美xxxx黑人xx丫x性爽| 在线观看免费视频日本深夜| 午夜免费激情av| 精品一区二区三区视频在线观看免费| 久久午夜亚洲精品久久| 1000部很黄的大片| 国产美女午夜福利| 成人国产一区最新在线观看| 狂野欧美白嫩少妇大欣赏| 久久久久亚洲av毛片大全| 国产69精品久久久久777片| 亚洲成人精品中文字幕电影| 老司机午夜福利在线观看视频| 少妇人妻一区二区三区视频| 狂野欧美激情性xxxx| 最新中文字幕久久久久| 亚洲久久久久久中文字幕| 黄色女人牲交| xxx96com| 亚洲国产欧洲综合997久久,| 午夜精品一区二区三区免费看| 九色成人免费人妻av| 亚洲午夜理论影院| 香蕉丝袜av| 激情在线观看视频在线高清| 午夜福利在线观看吧| 全区人妻精品视频| 好男人在线观看高清免费视频| 精品电影一区二区在线| 国产日本99.免费观看| 午夜福利18| 少妇的丰满在线观看| 日本 av在线| 亚洲av成人精品一区久久| 天堂√8在线中文| 中文字幕高清在线视频| 在线观看美女被高潮喷水网站 | bbb黄色大片| 亚洲欧美日韩无卡精品| 国产久久久一区二区三区| 亚洲18禁久久av| 窝窝影院91人妻| 欧美性猛交╳xxx乱大交人| 少妇熟女aⅴ在线视频|