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

    基于隨機(jī)數(shù)值模擬的溢油對(duì)敏感資源風(fēng)險(xiǎn)評(píng)價(jià)研究

    2017-02-27 11:43:48
    關(guān)鍵詞:污染區(qū)域環(huán)境

    鄧 健 孫 浩 謝 澄

    (武漢理工大學(xué)航運(yùn)學(xué)院1) 武漢 430063) (內(nèi)河航運(yùn)技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室2) 武漢 430063) (交通運(yùn)輸部長(zhǎng)江航運(yùn)技術(shù)行業(yè)研發(fā)中心3) 武漢 430000)

    基于隨機(jī)數(shù)值模擬的溢油對(duì)敏感資源風(fēng)險(xiǎn)評(píng)價(jià)研究

    鄧 健1,2,3)孫 浩1)謝 澄1)

    (武漢理工大學(xué)航運(yùn)學(xué)院1)武漢 430063) (內(nèi)河航運(yùn)技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室2)武漢 430063) (交通運(yùn)輸部長(zhǎng)江航運(yùn)技術(shù)行業(yè)研發(fā)中心3)武漢 430000)

    針對(duì)船舶溢油事故對(duì)沿海經(jīng)濟(jì)體系及海洋生態(tài)環(huán)境所造成的嚴(yán)重?fù)p害,以防城港8萬(wàn)t級(jí)化工碼頭為研究對(duì)象,基于潮流擴(kuò)散、漂移等模型,采用隨機(jī)統(tǒng)計(jì)模擬法,在變化潮流、風(fēng)力等影響溢油運(yùn)行軌跡的相關(guān)因素情況下,對(duì)港口區(qū)域溢油污染的概率及溢油污染的最短時(shí)間2個(gè)風(fēng)險(xiǎn)因素進(jìn)行統(tǒng)計(jì)與分析,確定評(píng)價(jià)碼頭發(fā)生溢油事故后對(duì)環(huán)境敏感資源的風(fēng)險(xiǎn)大小.

    安全工程;隨機(jī)模擬;溢油風(fēng)險(xiǎn);風(fēng)險(xiǎn)評(píng)價(jià)

    0 引 言

    隨著航運(yùn)業(yè)的不斷發(fā)展,溢油事故時(shí)有發(fā)生,其對(duì)軌跡覆蓋區(qū)域內(nèi)養(yǎng)殖區(qū)、水域環(huán)境等造成重大影響,所造成的直接或間接損失不可估量.溢油事故一般為船舶碰撞或擱淺的次生事故,具有突發(fā)性,事故發(fā)生后,溢油將迅速擴(kuò)散,對(duì)岸線及相關(guān)水域造成嚴(yán)重污染,因此,對(duì)溢油行為的預(yù)測(cè)研究一直是各學(xué)者研究的熱點(diǎn).李穎等[1]利用數(shù)值建模結(jié)合地理信息系統(tǒng)(GIS)分析了在墨西哥灣油品泄漏對(duì)岸線及生態(tài)的影響預(yù)測(cè).Mohammad[2]結(jié)合環(huán)境因素、船舶交通事故等因素建立了一個(gè)動(dòng)態(tài)的風(fēng)險(xiǎn)模型,來(lái)確定區(qū)域風(fēng)險(xiǎn)狀況,并在實(shí)時(shí)海況下對(duì)船舶溢油事故進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià).目前許多學(xué)者的研究主要集中于單次油粒子的擴(kuò)散模型,忽略風(fēng)場(chǎng)流場(chǎng)等因素的隨機(jī)性,將復(fù)雜多變的隨機(jī)輸入、輸出關(guān)系以確定性的輸入和輸出求解數(shù)學(xué)模型,得到確定性的輸出,對(duì)輸出的確定性結(jié)果也不作概率統(tǒng)計(jì)方面的分析,這種模擬方法與真實(shí)結(jié)果存在較大誤差.

    隨機(jī)數(shù)值模擬是一種通過(guò)設(shè)定隨機(jī)過(guò)程,反復(fù)生成序列,研究其分布的方法,主要用于易變指標(biāo)(如風(fēng)、流等多變條件)的風(fēng)險(xiǎn)預(yù)測(cè).Bardossy等[3]使用隨機(jī)數(shù)值模擬和ARMA過(guò)程的參數(shù),將德國(guó)西南32個(gè)降水站的40年歷史觀察與100個(gè)匹配時(shí)間進(jìn)行對(duì)比,使用大量不同的診斷從而對(duì)德國(guó)降水量進(jìn)行估計(jì).隨機(jī)模擬在水上交通也早有相關(guān)應(yīng)用.Moonjin[4]采用油粒子模型,基于隨機(jī)模擬的方法,對(duì)溢油到達(dá)指定區(qū)域的概率和時(shí)間進(jìn)行了分析,得出了韓國(guó)大山港的溢油評(píng)價(jià)總體結(jié)論.從國(guó)內(nèi)外研究可知,在溢油的預(yù)測(cè)方面,研究目標(biāo)以溢油的到達(dá)指定區(qū)域的概率和最短時(shí)間為主.

    隨機(jī)數(shù)值模擬計(jì)算結(jié)果只依賴于隨機(jī)過(guò)程中序列的分布,受人為和特定自然環(huán)境因素影響較小,能得到較為客觀、全面的預(yù)測(cè)結(jié)果.文中選取國(guó)內(nèi)典型港口作為樣本,分析其敏感資源分布情況,以300次隨機(jī)溢油擴(kuò)散場(chǎng)景為序列基礎(chǔ),得到溢油事故污染影響的概率分布、到達(dá)敏感區(qū)域擴(kuò)散時(shí)間分布,以期為溢油風(fēng)險(xiǎn)快速預(yù)測(cè)提供理論支撐.

    1 評(píng)價(jià)模型的建立

    1.1 潮流模型

    潮流是近岸淺水海域的基本流動(dòng),潮流數(shù)值計(jì)算是模擬入海溢油隨潮漂移和歸宿的基礎(chǔ).文中采用三維DELFT潮流模型來(lái)模擬評(píng)價(jià)碼頭附近水域的潮流場(chǎng).三維潮流模型用于描述正壓海洋的深度平均運(yùn)動(dòng),方程組由連續(xù)性方程和動(dòng)量方程組成[5].

    沿水深積分的連續(xù)方程:

    (1)

    ξ和η方向的動(dòng)量方程:

    ξ方向:

    (2)

    η方向:

    (3)

    式中:ζ為水位;d為水深;Gζζ,Gη為直角坐標(biāo)系與正交曲線坐標(biāo)系轉(zhuǎn)換系數(shù);U,V為ζ和η方向上的水深平均速度分量;Q為單位面積上取排水、蒸發(fā)和降水對(duì)流量的貢獻(xiàn);f為柯氏力系數(shù);Fζ,F(xiàn)η為ζ和η方向的紊動(dòng)動(dòng)量通量;Pζ,Pη為ζ和η方向的水壓力梯度;Mζ,Mη為ζ和η方向的動(dòng)量的源或匯;Vv為垂向渦動(dòng)系數(shù).

    1.2 溢油事故后果預(yù)測(cè)模型

    事故后果采用油粒子模型對(duì)溢油的漂移、風(fēng)化、擴(kuò)散等一系列過(guò)程進(jìn)行模擬.其采用的數(shù)值模型主要包括:漂移模型和擴(kuò)散模型.

    假設(shè)溢油可分成許多獨(dú)立溢油粒,將每個(gè)溢油??醋骼窭嗜樟W?溢油粒在t時(shí)刻的位置分量用Xt表示.

    (4)

    式中:Δt為時(shí)間步長(zhǎng);Xt-1為t-1時(shí)刻溢油粒子位置;Uoil為溢油粒漂移速度.

    溢油粒子漂移速度Uoil可定義為

    (5)

    式中:Ut為水流產(chǎn)生的速度分量;Uw為風(fēng)產(chǎn)生的速度分量;Udisp為分散過(guò)程產(chǎn)生的速度分量.

    風(fēng)產(chǎn)生的速度分量Uw可定義為

    (6)

    式中:風(fēng)漂移因子C1是一個(gè)常量,在1%~4.5%變化,由實(shí)測(cè)確定,3%~3.5%對(duì)開(kāi)闊海域較合適.對(duì)封閉或半封閉海域,C1值較小;Uwindspeed為風(fēng)速.

    東方向擴(kuò)散速度分量可定義為

    (7)

    北方向擴(kuò)散速度分量可定義為

    (8)

    式中:Dx為東西方向水平擴(kuò)散系數(shù);Dy為南北方向水平擴(kuò)散系數(shù);Δt表示時(shí)間步長(zhǎng);r為隨機(jī)系數(shù)(-1~1).水平擴(kuò)散系數(shù)Dx,Dy通常相等.

    1.3 風(fēng)險(xiǎn)評(píng)價(jià)方法

    依據(jù)國(guó)內(nèi)外學(xué)者的研究情況,溢油污染風(fēng)險(xiǎn)可以根據(jù)水面遭受溢油污染的概率和溢油污染的最短時(shí)間兩方面進(jìn)行評(píng)價(jià).

    溢油污染的概率是評(píng)定敏感區(qū)域是否受到溢油污染的重要指標(biāo),該信息可以為有效的分配溢油應(yīng)急資源提供支持.溢油影響的概率可定義為溢油粒子到達(dá)網(wǎng)格區(qū)域的次數(shù)與模擬總次數(shù)的比值,可以用公式

    (9)

    式中:P為溢油影響的概率;n為溢油達(dá)到網(wǎng)格區(qū)域的次數(shù);N為模擬總次數(shù).

    溢油污染的最短時(shí)間可以提供優(yōu)先響應(yīng)信息,為提前應(yīng)對(duì)溢油污染做準(zhǔn)備.溢油達(dá)到的最短時(shí)間可定義為溢油粒子從事故地點(diǎn)到達(dá)到每個(gè)網(wǎng)格區(qū)域所需要的時(shí)間,可以用公式

    t=T-T0

    (10)

    式中:t為溢油污染的最短時(shí)間;T為溢油達(dá)到網(wǎng)格的時(shí)間;T0為發(fā)生事故時(shí)間.由于在不同條件下進(jìn)行每次隨機(jī)模擬實(shí)驗(yàn)時(shí)溢油達(dá)到同一區(qū)域的時(shí)間均不相同,因此取其最小值.

    2 典型碼頭實(shí)例分析

    2.1 評(píng)價(jià)區(qū)域概述

    防城港背靠大西南、面向東南亞、東臨粵港澳,是我國(guó)沿海主要港口之一[6].港區(qū)內(nèi)有危險(xiǎn)品碼頭、水產(chǎn)養(yǎng)殖區(qū)、環(huán)境保護(hù)區(qū)等敏感資源.一旦發(fā)生溢油事故,將會(huì)對(duì)海洋生態(tài)環(huán)境和港口經(jīng)濟(jì)造成難以挽回的巨大損失.港區(qū)內(nèi)化工碼頭眾多,因其貨種的高風(fēng)險(xiǎn)性,油品泄漏后可能遇熱源產(chǎn)生爆炸等連鎖事故.風(fēng)險(xiǎn)性相較于其他散貨碼頭更高,因此文中選取防城港東灣8萬(wàn)t級(jí)化工碼頭作為事故發(fā)生的模擬位置.

    碼頭主要為周邊地區(qū)的企業(yè)提供成品油、化學(xué)品等原材料和產(chǎn)品運(yùn)輸服務(wù).據(jù)統(tǒng)計(jì),碼頭年吞吐量147萬(wàn)t,主力船型為5萬(wàn)t級(jí)原油船.據(jù)海事部門統(tǒng)計(jì),1994—2013年防城港及其附近海域共發(fā)生船舶交通事故33起,船舶交通流量共計(jì)66 667艘次,港口貨物吞吐量共計(jì)43 478萬(wàn)t.在此時(shí)間內(nèi),防城港港海域發(fā)生操作性船舶溢油事故1起,海難性船舶溢油事故2起.可計(jì)算出船舶污染次數(shù)發(fā)生的概率,見(jiàn)表1.

    表1 船舶污染事故發(fā)生頻率

    根據(jù)《船舶污染海洋環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)規(guī)范》給出的預(yù)測(cè)方法,操作性事故最可能發(fā)生溢油量以3分鐘泄漏量為準(zhǔn),海難性最可能發(fā)生泄漏量以主力船型單個(gè)艙室油品泄漏量為準(zhǔn),結(jié)合本碼頭的運(yùn)輸船舶的主要船型、噸位和實(shí)載率,確定本碼頭最可能發(fā)生的操作性溢油量約為25 t,最可能發(fā)生事故溢油量原油360 t.

    2.2 環(huán)境敏感資源概述

    根據(jù)《防城港市海洋功能區(qū)劃》《廣西壯族自治區(qū)海洋功能區(qū)劃》對(duì)各類海洋資源的劃分[7],可得防城港區(qū)環(huán)境敏感目標(biāo)如表2

    表2 環(huán)境敏感目標(biāo)

    環(huán)境敏感指數(shù)圖與溢油事故應(yīng)急響應(yīng)息息相關(guān),是選擇溢油污染風(fēng)險(xiǎn)評(píng)價(jià)區(qū)域的重要因素.根據(jù)《船舶污染海洋環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)規(guī)范(試行)》中環(huán)境資源敏感性分類表,將評(píng)價(jià)區(qū)域環(huán)境敏感度分為5檔(見(jiàn)圖1),其中環(huán)境敏感系數(shù)0表示環(huán)境敏感程度很低,1表示環(huán)境敏感程度低,5表示環(huán)境敏感程度中等,20表示環(huán)境敏感程度高,50表示環(huán)境敏感程度非常高.環(huán)境敏感程度非常高區(qū)域主要為紅樹(shù)林保護(hù)區(qū),環(huán)境敏感程度高的區(qū)域主要為海洋保護(hù)區(qū)和國(guó)家級(jí)的養(yǎng)殖區(qū),環(huán)境敏感程度中等區(qū)域主要為一般的風(fēng)景旅游區(qū)和潮間帶等,環(huán)境敏感程度低的區(qū)域主要為敏感度的岸線,環(huán)境敏感程度很低的區(qū)域主要為開(kāi)闊的海域.

    圖1 海域環(huán)境保護(hù)目標(biāo)敏感系數(shù)及網(wǎng)格分布圖

    2.3 隨機(jī)數(shù)值模擬方法

    圖2為評(píng)價(jià)區(qū)域漲潮表層流場(chǎng)圖,由圖2可知,對(duì)評(píng)價(jià)碼頭及其周邊評(píng)價(jià)區(qū)域自西向東,自南向北劃分了51×53個(gè)網(wǎng)格.確定各環(huán)境敏感資源所在網(wǎng)格,根據(jù)風(fēng)、潮流、溢油泄漏數(shù)量等不同條件下,依據(jù)1.3方法進(jìn)行評(píng)價(jià)區(qū)域內(nèi)的溢油模擬.

    2.3.1 隨機(jī)模擬實(shí)驗(yàn)參數(shù)選擇

    通過(guò)查詢國(guó)內(nèi)外相關(guān)研究可知,對(duì)于船舶溢油污染的影響為流場(chǎng)、風(fēng)場(chǎng)、時(shí)間、地點(diǎn)、數(shù)量的集合[8-9],為此,在隨機(jī)模擬參數(shù)的選取可按如下函數(shù)表示.

    通過(guò)美國(guó)地球物理數(shù)據(jù)中心(NGDC)提供的高精度海岸線邊界數(shù)據(jù),基于1.1潮流模型,采用曲線正交網(wǎng)格剖分的方法模擬評(píng)價(jià)區(qū)域范圍內(nèi)的潮流情況,見(jiàn)圖2~3.

    圖2 評(píng)價(jià)區(qū)域漲潮表層流場(chǎng)圖

    圖3 評(píng)價(jià)區(qū)域落潮表層流場(chǎng)圖

    由模擬結(jié)果可見(jiàn)評(píng)價(jià)區(qū)域內(nèi)大潮時(shí)落潮流速大于漲潮流速,防城港灣外海域漲潮以偏NE向?yàn)橹鳎涑币云玈W向?yàn)橹?,外海潮流流速比?nèi)灣大,灣頸處達(dá)最大.模擬結(jié)果與2013年11月大潮期觀測(cè)結(jié)果吻合較好.

    圖4為防城港市的風(fēng)玫瑰圖,由圖4可知,全年常風(fēng)向?yàn)镹NE,頻率為30.5%,次常風(fēng)向?yàn)镾SW,其頻率為8.4%;強(qiáng)風(fēng)向?yàn)镋,頻率為4.9%,其最大風(fēng)速為36 m/s,次強(qiáng)風(fēng)向?yàn)镹NE,其最大風(fēng)速為27 m/s.多年平均風(fēng)速約為5 m/s.

    圖4 防城港風(fēng)玫瑰圖

    結(jié)合時(shí)間、地點(diǎn)等影響因素,給出如表3所示隨機(jī)模擬參數(shù).

    表3 隨機(jī)統(tǒng)計(jì)模擬計(jì)算參數(shù)表

    2.3.2 隨機(jī)模擬個(gè)例分析

    隨機(jī)選取夏季SW風(fēng)向落潮,風(fēng)速2.9 m/s,原油15 t,等相關(guān)條件進(jìn)行模擬,模擬結(jié)果見(jiàn)圖5.

    圖5 溢油軌跡

    2.3.3 溢油風(fēng)險(xiǎn)評(píng)價(jià)

    根據(jù)2.3.2模擬方法進(jìn)行隨機(jī)條件下300次的隨機(jī)模擬實(shí)驗(yàn),對(duì)結(jié)果進(jìn)行統(tǒng)計(jì),得出水面受到溢油污染的概率及最短時(shí)間分布,見(jiàn)圖6~7.

    由圖可見(jiàn),事故發(fā)生后,72 h內(nèi)可能污染的面積范圍較大,其中污染概率較大的為本工程岸線所在水域,概率超過(guò)70%以上.溢油對(duì)環(huán)境敏感資源影響的概率及時(shí)間具體見(jiàn)表4.

    圖6 水面溢油污染概率分布圖

    圖7 水面受污染最短時(shí)間分布圖

    環(huán)境敏感資源溢油污染概率/%溢油污染最短時(shí)間/h工程所在岸線902東灣海洋保護(hù)區(qū)307江山半島東岸旅游休閑娛樂(lè)區(qū)208江山半島南部農(nóng)漁業(yè)區(qū)209企沙半島南部農(nóng)漁業(yè)區(qū)1020東灣紅樹(shù)林528

    2.3.4 整體風(fēng)險(xiǎn)評(píng)價(jià)

    風(fēng)險(xiǎn)為概率與風(fēng)險(xiǎn)后果的乘積:

    風(fēng)險(xiǎn)R=概率P×后果C

    其風(fēng)險(xiǎn)概率由溢油污染概率可知,而風(fēng)險(xiǎn)后果除與遭受溢油污染的最短時(shí)間有關(guān)外,敏感資源自身的敏感系數(shù)也起著相當(dāng)重要的作用.因此我們可定義評(píng)價(jià)碼頭的風(fēng)險(xiǎn)為

    R=Pn/T

    式中:R為風(fēng)險(xiǎn);P為溢油影響的概率;n為敏感系數(shù);T為溢油影響的時(shí)間.

    得出各敏感區(qū)域的溢油風(fēng)險(xiǎn),見(jiàn)表5.

    通過(guò)表5所計(jì)算出風(fēng)險(xiǎn)大小,確定該工程所在海域的高風(fēng)險(xiǎn)區(qū)域?yàn)楣こ趟诎毒€、東灣海洋保護(hù)區(qū)、江山半島東岸旅游休閑娛樂(lè)區(qū)、江山半島南部農(nóng)漁業(yè)區(qū).

    表5 敏感區(qū)域溢油風(fēng)險(xiǎn)

    3 結(jié) 論

    1) 溢油污染風(fēng)險(xiǎn)受到實(shí)時(shí)海況下的潮流、風(fēng)等因素的影響.利用潮流模型和預(yù)測(cè)模型對(duì)漂移、分散、擴(kuò)展等一系列過(guò)程進(jìn)行模擬,分析溢油泄漏后的運(yùn)行軌跡.由于海洋環(huán)境變化性,為提高模擬系統(tǒng)的精確性,引入隨機(jī)模擬理論,在不同條件下對(duì)模擬事故發(fā)生地點(diǎn)進(jìn)行300次的隨機(jī)模擬實(shí)驗(yàn),分析水面遭受溢油污染的概率和溢油污染的最短時(shí)間,很好的反應(yīng)出評(píng)價(jià)區(qū)域里的敏感資源的風(fēng)險(xiǎn)情況

    2) 溢油風(fēng)險(xiǎn)研究雖已十分成熟,但因其受特定環(huán)境(特定風(fēng)速、流場(chǎng)等)因素較大,不具備普遍的適用性和預(yù)測(cè)性.隨機(jī)模擬在評(píng)價(jià)過(guò)程中只考慮所有序列的分布情況,其溢油影響預(yù)測(cè)可借鑒程度高,為制定港口溢油污染對(duì)策及建立溢油應(yīng)急響應(yīng)機(jī)制提供有力的支持.

    [1]李穎,蘭國(guó)新,劉丙新,等.基于遙感和GIS技術(shù)的墨西哥灣溢油動(dòng)態(tài)分析[J].海洋環(huán)境科學(xué),2012(1):79-82.

    [2]MOHAMMAD T T. Experimental investigation of oil leakage from damaged ships due to collision and grounding[J]. Ocean Engineering,2011,38(17):1894-1907.

    [3]BARDOSSY A,PEGRAM G S.Copula based multisite model for daily precipitation simulation[J]. Hydrologic Earth System Science,2009(1):55-58.

    [4]MOONJIN L.Pollution risk assessment of oil spill accidents in Garorim bay of Korea[J]. Marine Pollution Bulletin,2015(2):297-303.

    [5]應(yīng)超.二維潮流數(shù)學(xué)模型建設(shè)及其可視化研究[D].天津:天津大學(xué),2012.

    [6]祁昆源.防城港轄區(qū)船舶防污染管理研究[D].大連:大連海事大學(xué),2012.

    [7]楊紅,韓眾,李曰嵩,等.東海西湖凹陷區(qū)溢油行為數(shù)值模擬[J].海洋環(huán)境科學(xué),2012(2):88-93.

    [8]CHENG Y C, LI X F, XU Q. SAR observation and modeltracking of an oil spill event in coastal waters[J]. Marine Pollution Bulletin,2011,62(2):350-363.

    [9]SEBASTIAO P, GUEDES S C. Uncertainty in predictions of oil spill trajectories in open sea[J]. Ocean Engineering,2007,34(3):576-584.

    Study on the Risk Assessment of Oil Spill Based on Stochastic Numerical Simulation for the Pollution Risk of Sensitive Resources

    DENG Jian1,2,3)SUN Hao1)XIE Cheng1)

    (Schoolofnavigation,WuhanUniversityofTechnology,Wuhan430063,China)1)(HubeiInlandShippingTechnologyKeyLaboratory,Wuhan430063,China)2)(TheMinistryoftransportandshippingoftheYangtseRiverTechnologyIndustryResearchandDevelopmentCenter,Wuhan430000,China)3)

    Aiming at the serious damage of the coastal economic system and marine ecological environment caused by the oil spill accidents of ships and taking the 80 000-ton chemical port in Fang Cheng-gang as an example, two risk factors including the probability of oil spill pollution in the port area and the shortest time of oil spill pollution are analyzed on the condition that the related factors of the running track of oil spill are affected by the changed flow and wind. Meanwhile, the analysis is based on the model of tidal current diffusion and drift, and the stochastic simulation method is used. Thereby, the risk of environmental sensitive resources after the oil spill accidents can be evaluated and determined.

    safety engineering; stochastic simulation; oil spill risk; risk assessment

    2016-10-25

    U698.7

    10.3963/j.issn.2095-3844.2017.01.007

    鄧健(1977—):男,博士,副教授,主要研究領(lǐng)域?yàn)樗辖煌ò踩c防污染

    猜你喜歡
    污染區(qū)域環(huán)境
    長(zhǎng)期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    孕期遠(yuǎn)離容易致畸的環(huán)境
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    環(huán)境
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    對(duì)抗塵污染,遠(yuǎn)離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    97人妻精品一区二区三区麻豆| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91久久精品国产一区二区成人| 我要搜黄色片| 亚洲欧美清纯卡通| 中文字幕制服av| 99热这里只有是精品50| 最近2019中文字幕mv第一页| 热99re8久久精品国产| 亚洲国产色片| 51国产日韩欧美| 性色avwww在线观看| 亚洲最大成人手机在线| 男人的好看免费观看在线视频| 亚洲国产精品成人综合色| 两个人视频免费观看高清| a级毛片免费高清观看在线播放| 久久精品影院6| 天堂√8在线中文| 日韩中字成人| 免费av不卡在线播放| 国产精品一区二区三区四区免费观看| 97超碰精品成人国产| 免费av观看视频| 久久这里有精品视频免费| 噜噜噜噜噜久久久久久91| 久久这里只有精品中国| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品一区二区三区视频在线| 美女内射精品一级片tv| 我要搜黄色片| 青青草视频在线视频观看| 亚洲七黄色美女视频| 97热精品久久久久久| 午夜精品在线福利| 国产成人福利小说| 欧美性感艳星| 99九九线精品视频在线观看视频| 偷拍熟女少妇极品色| a级毛片a级免费在线| 亚洲欧洲日产国产| 国产av不卡久久| 女人十人毛片免费观看3o分钟| 国产淫片久久久久久久久| 性欧美人与动物交配| 国产精品不卡视频一区二区| 国内精品美女久久久久久| 黄片wwwwww| 长腿黑丝高跟| 一级黄片播放器| 国产精品麻豆人妻色哟哟久久 | 欧美成人精品欧美一级黄| 欧美潮喷喷水| 99国产极品粉嫩在线观看| a级毛片免费高清观看在线播放| 一本久久中文字幕| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 国产午夜福利久久久久久| 3wmmmm亚洲av在线观看| 91麻豆精品激情在线观看国产| 卡戴珊不雅视频在线播放| 国产中年淑女户外野战色| 国产 一区 欧美 日韩| 午夜免费男女啪啪视频观看| 91aial.com中文字幕在线观看| 变态另类丝袜制服| 最近手机中文字幕大全| 一本精品99久久精品77| 婷婷亚洲欧美| 一区福利在线观看| 日韩av在线大香蕉| 成熟少妇高潮喷水视频| 99视频精品全部免费 在线| 免费观看a级毛片全部| 亚洲av.av天堂| 伦理电影大哥的女人| 91久久精品国产一区二区成人| 精品久久久久久久久av| 九色成人免费人妻av| 亚洲av免费在线观看| 国产亚洲91精品色在线| 国产精品精品国产色婷婷| 综合色丁香网| 99国产精品一区二区蜜桃av| 女人十人毛片免费观看3o分钟| 国产成人a∨麻豆精品| 性插视频无遮挡在线免费观看| 午夜老司机福利剧场| 国产黄片视频在线免费观看| 少妇猛男粗大的猛烈进出视频 | 黑人高潮一二区| 国产成人a∨麻豆精品| 晚上一个人看的免费电影| 69av精品久久久久久| 欧美色视频一区免费| 国产精品野战在线观看| av国产免费在线观看| 亚洲欧洲日产国产| 亚洲av不卡在线观看| 日韩,欧美,国产一区二区三区 | 校园人妻丝袜中文字幕| 夜夜夜夜夜久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品国产亚洲av香蕉五月| 亚洲婷婷狠狠爱综合网| 日韩大尺度精品在线看网址| 亚洲精品自拍成人| 嫩草影院精品99| 青春草国产在线视频 | 哪里可以看免费的av片| 秋霞在线观看毛片| 国产成人一区二区在线| 久久人人精品亚洲av| 久久久久久国产a免费观看| 亚洲欧美清纯卡通| 国产亚洲av片在线观看秒播厂 | 亚洲国产色片| 国产成人精品一,二区 | 秋霞在线观看毛片| 在线播放无遮挡| 久久久久久九九精品二区国产| 久久99热这里只有精品18| 国产高清有码在线观看视频| 国产精品蜜桃在线观看 | 一边亲一边摸免费视频| 亚洲五月天丁香| 女的被弄到高潮叫床怎么办| 免费看光身美女| 最新中文字幕久久久久| 91在线精品国自产拍蜜月| 国产精品一二三区在线看| 欧美激情国产日韩精品一区| 国产三级在线视频| 日日摸夜夜添夜夜爱| 三级毛片av免费| 免费人成在线观看视频色| 免费无遮挡裸体视频| 热99在线观看视频| 18禁在线无遮挡免费观看视频| 欧美另类亚洲清纯唯美| 久久这里有精品视频免费| 你懂的网址亚洲精品在线观看 | 深夜精品福利| 女人被狂操c到高潮| 婷婷精品国产亚洲av| 日韩中字成人| 偷拍熟女少妇极品色| 国产精品久久电影中文字幕| 夫妻性生交免费视频一级片| 在线免费观看的www视频| 人妻少妇偷人精品九色| 日韩一区二区视频免费看| 精品人妻偷拍中文字幕| 久久人人爽人人爽人人片va| 亚洲国产色片| 国产免费男女视频| 日日啪夜夜撸| 成人特级av手机在线观看| 欧美3d第一页| 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 午夜免费男女啪啪视频观看| 日韩中字成人| 国产精华一区二区三区| 日韩国内少妇激情av| 日韩中字成人| 国产成人aa在线观看| 一个人看视频在线观看www免费| 美女大奶头视频| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| .国产精品久久| 高清日韩中文字幕在线| 桃色一区二区三区在线观看| 日韩,欧美,国产一区二区三区 | 五月玫瑰六月丁香| 人妻少妇偷人精品九色| 91狼人影院| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 精品国内亚洲2022精品成人| 国产精品综合久久久久久久免费| 69人妻影院| 亚洲久久久久久中文字幕| 国产激情偷乱视频一区二区| 成人永久免费在线观看视频| 欧美日韩国产亚洲二区| 免费一级毛片在线播放高清视频| 内射极品少妇av片p| 精品熟女少妇av免费看| 欧美又色又爽又黄视频| 中文精品一卡2卡3卡4更新| 免费观看人在逋| 给我免费播放毛片高清在线观看| 国产精品蜜桃在线观看 | 亚洲在久久综合| 欧美日本亚洲视频在线播放| 国产成人精品婷婷| 夜夜爽天天搞| 免费观看的影片在线观看| 国产69精品久久久久777片| 成人永久免费在线观看视频| 久久婷婷人人爽人人干人人爱| 国产成人影院久久av| 国产在视频线在精品| 91午夜精品亚洲一区二区三区| 久久综合国产亚洲精品| 国内精品久久久久精免费| 久久99蜜桃精品久久| 一本久久中文字幕| 少妇高潮的动态图| 成年女人永久免费观看视频| 直男gayav资源| 精品人妻偷拍中文字幕| 美女cb高潮喷水在线观看| 欧美高清成人免费视频www| 国产高清激情床上av| 能在线免费观看的黄片| 美女脱内裤让男人舔精品视频 | 别揉我奶头 嗯啊视频| 麻豆久久精品国产亚洲av| 少妇裸体淫交视频免费看高清| 中文字幕人妻熟人妻熟丝袜美| 亚洲激情五月婷婷啪啪| 亚洲精品国产av成人精品| 97超碰精品成人国产| 国产蜜桃级精品一区二区三区| 最近的中文字幕免费完整| 伦理电影大哥的女人| 国国产精品蜜臀av免费| 尤物成人国产欧美一区二区三区| 日韩欧美精品v在线| 久久国内精品自在自线图片| 日本熟妇午夜| 丰满乱子伦码专区| 人妻久久中文字幕网| 亚洲第一区二区三区不卡| 偷拍熟女少妇极品色| 又爽又黄a免费视频| 又粗又爽又猛毛片免费看| 99久久无色码亚洲精品果冻| 亚洲国产日韩欧美精品在线观看| 欧美精品国产亚洲| 亚洲国产精品国产精品| 国产一区二区三区av在线 | 又粗又爽又猛毛片免费看| 日韩av在线大香蕉| 91狼人影院| 亚洲av成人av| 级片在线观看| av免费观看日本| 在线免费十八禁| 久久鲁丝午夜福利片| 最新中文字幕久久久久| 日本黄色片子视频| 亚洲第一区二区三区不卡| 免费观看人在逋| 亚洲欧美日韩高清在线视频| 99国产极品粉嫩在线观看| 性色avwww在线观看| 欧美潮喷喷水| 婷婷六月久久综合丁香| 色视频www国产| 亚洲中文字幕日韩| 日韩强制内射视频| 如何舔出高潮| 亚洲中文字幕一区二区三区有码在线看| 亚洲av一区综合| 久久久久久久亚洲中文字幕| 精品一区二区免费观看| 男插女下体视频免费在线播放| 国产伦一二天堂av在线观看| 欧美精品国产亚洲| 亚洲美女搞黄在线观看| 成人三级黄色视频| 亚洲精品久久国产高清桃花| 大又大粗又爽又黄少妇毛片口| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看| 国产综合懂色| 免费av不卡在线播放| av国产免费在线观看| 毛片一级片免费看久久久久| 久久久国产成人精品二区| 欧美变态另类bdsm刘玥| 边亲边吃奶的免费视频| 亚洲av成人av| 午夜福利视频1000在线观看| 直男gayav资源| 黄色一级大片看看| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 国产成人午夜福利电影在线观看| 精华霜和精华液先用哪个| videossex国产| 国产av麻豆久久久久久久| 热99在线观看视频| 青春草国产在线视频 | 日本av手机在线免费观看| 亚洲av中文字字幕乱码综合| 高清日韩中文字幕在线| 干丝袜人妻中文字幕| 精品一区二区三区人妻视频| 久久精品国产清高在天天线| 一进一出抽搐动态| 麻豆一二三区av精品| 国内精品美女久久久久久| 中国国产av一级| 国产精品一二三区在线看| 在线播放国产精品三级| 不卡视频在线观看欧美| 久久久精品欧美日韩精品| 亚洲久久久久久中文字幕| 成年女人看的毛片在线观看| 夜夜爽天天搞| 亚洲性久久影院| 国产成人福利小说| 天堂影院成人在线观看| 久久精品国产清高在天天线| 国产精品蜜桃在线观看 | 超碰av人人做人人爽久久| 亚洲人成网站在线播| 亚洲乱码一区二区免费版| 久久精品夜色国产| 两个人的视频大全免费| 日本免费a在线| 一级毛片aaaaaa免费看小| 日韩欧美精品v在线| 国产黄色视频一区二区在线观看 | 简卡轻食公司| 国产av麻豆久久久久久久| 国产高清不卡午夜福利| 爱豆传媒免费全集在线观看| 嫩草影院新地址| 国产精品美女特级片免费视频播放器| 国产大屁股一区二区在线视频| 岛国毛片在线播放| av.在线天堂| 内射极品少妇av片p| 最近2019中文字幕mv第一页| 亚洲三级黄色毛片| 2021天堂中文幕一二区在线观| а√天堂www在线а√下载| 老女人水多毛片| av天堂在线播放| 五月玫瑰六月丁香| 久久精品国产99精品国产亚洲性色| 岛国毛片在线播放| 日韩人妻高清精品专区| 成年女人永久免费观看视频| 欧美另类亚洲清纯唯美| 亚洲精品亚洲一区二区| 一区二区三区四区激情视频 | 亚洲第一区二区三区不卡| 亚洲内射少妇av| 蜜桃亚洲精品一区二区三区| 三级男女做爰猛烈吃奶摸视频| 女人十人毛片免费观看3o分钟| 亚洲人成网站在线观看播放| 欧美性猛交黑人性爽| 成人一区二区视频在线观看| 国国产精品蜜臀av免费| 日韩在线高清观看一区二区三区| 久久久久久大精品| 国产精品电影一区二区三区| 国产白丝娇喘喷水9色精品| а√天堂www在线а√下载| 国产爱豆传媒在线观看| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 成人毛片a级毛片在线播放| 欧美性猛交╳xxx乱大交人| 别揉我奶头 嗯啊视频| 日韩av不卡免费在线播放| 男女视频在线观看网站免费| 成人综合一区亚洲| 亚洲自拍偷在线| 午夜精品在线福利| 欧美人与善性xxx| 国产男人的电影天堂91| 国产精品久久视频播放| 日韩成人伦理影院| 国产高清视频在线观看网站| 国产精品蜜桃在线观看 | 国产伦一二天堂av在线观看| 热99在线观看视频| 午夜精品在线福利| 特级一级黄色大片| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 亚洲成人久久性| 夜夜爽天天搞| 国产激情偷乱视频一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 天天躁日日操中文字幕| 麻豆一二三区av精品| 国产精品久久久久久精品电影| 三级毛片av免费| 在线观看免费视频日本深夜| 噜噜噜噜噜久久久久久91| 最近最新中文字幕大全电影3| 日本五十路高清| 国内精品宾馆在线| 老熟妇乱子伦视频在线观看| 久久久久久大精品| 精品久久久久久久久久免费视频| 欧美又色又爽又黄视频| 一区二区三区免费毛片| 日韩在线高清观看一区二区三区| 久久精品国产清高在天天线| 日本撒尿小便嘘嘘汇集6| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 亚洲在线自拍视频| 一区福利在线观看| 午夜激情福利司机影院| 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一电影网av| 18禁黄网站禁片免费观看直播| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看 | 国产精品.久久久| 亚洲精品国产成人久久av| 亚洲在线观看片| 日韩中字成人| 精品人妻熟女av久视频| 久久久久久久久久久免费av| 综合色av麻豆| 天堂网av新在线| 欧美日韩精品成人综合77777| 熟女人妻精品中文字幕| 夫妻性生交免费视频一级片| 观看美女的网站| 中文字幕av在线有码专区| www.色视频.com| 性欧美人与动物交配| 激情 狠狠 欧美| 男女那种视频在线观看| 亚洲18禁久久av| 男人狂女人下面高潮的视频| 亚洲精品成人久久久久久| 可以在线观看的亚洲视频| 亚洲激情五月婷婷啪啪| 亚洲成人久久爱视频| 国产在线精品亚洲第一网站| 免费av毛片视频| 国产伦精品一区二区三区视频9| 乱码一卡2卡4卡精品| 久久久久久伊人网av| 亚洲人与动物交配视频| 大香蕉久久网| 国产精品综合久久久久久久免费| 少妇人妻一区二区三区视频| 少妇猛男粗大的猛烈进出视频 | 色视频www国产| 少妇高潮的动态图| 美女国产视频在线观看| 简卡轻食公司| 中文字幕av在线有码专区| 国产精品一区二区三区四区久久| 熟女人妻精品中文字幕| 免费看光身美女| 欧美bdsm另类| 久久九九热精品免费| 国产伦精品一区二区三区四那| 国产精品.久久久| 国产私拍福利视频在线观看| 国产麻豆成人av免费视频| 亚洲在久久综合| 97超视频在线观看视频| 亚洲人成网站在线播| 哪个播放器可以免费观看大片| 日韩亚洲欧美综合| 国产成人freesex在线| АⅤ资源中文在线天堂| 中出人妻视频一区二区| 一级黄色大片毛片| 黑人高潮一二区| 国产视频内射| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 亚洲第一电影网av| avwww免费| 久久久国产成人免费| 精品国内亚洲2022精品成人| 亚洲四区av| 尾随美女入室| 美女黄网站色视频| 国产老妇伦熟女老妇高清| 成人高潮视频无遮挡免费网站| 97超视频在线观看视频| 国产精品一二三区在线看| .国产精品久久| 老司机影院成人| 亚洲电影在线观看av| 麻豆乱淫一区二区| 亚洲成人av在线免费| 在线观看午夜福利视频| 欧美精品一区二区大全| 男人狂女人下面高潮的视频| 亚洲精品色激情综合| 91精品一卡2卡3卡4卡| 男女视频在线观看网站免费| 高清午夜精品一区二区三区 | 国产在线精品亚洲第一网站| 黄色视频,在线免费观看| 精品日产1卡2卡| av.在线天堂| 国产真实乱freesex| 99精品在免费线老司机午夜| 免费观看人在逋| 欧美不卡视频在线免费观看| 日本av手机在线免费观看| 久久精品91蜜桃| 69人妻影院| 精品一区二区三区人妻视频| 高清毛片免费观看视频网站| 精品不卡国产一区二区三区| 一区二区三区免费毛片| 国产熟女欧美一区二区| 夜夜夜夜夜久久久久| 91av网一区二区| 国产淫片久久久久久久久| 亚洲成av人片在线播放无| 一级毛片aaaaaa免费看小| 国产爱豆传媒在线观看| 色播亚洲综合网| 一级毛片aaaaaa免费看小| 国产成人91sexporn| 可以在线观看的亚洲视频| 国产极品天堂在线| 三级毛片av免费| 日韩一本色道免费dvd| 国产亚洲av嫩草精品影院| 国产免费一级a男人的天堂| 日日干狠狠操夜夜爽| 在线播放无遮挡| 亚洲精品色激情综合| 99久久精品国产国产毛片| 欧美日韩综合久久久久久| 搡老妇女老女人老熟妇| 欧美日本视频| 18禁黄网站禁片免费观看直播| videossex国产| 可以在线观看的亚洲视频| .国产精品久久| 久久久精品欧美日韩精品| 级片在线观看| 五月玫瑰六月丁香| 国产 一区 欧美 日韩| 国产黄a三级三级三级人| 国产一级毛片在线| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 国产精品1区2区在线观看.| 亚洲国产欧美人成| 亚洲国产精品国产精品| 少妇熟女aⅴ在线视频| 日韩三级伦理在线观看| 又爽又黄无遮挡网站| 免费无遮挡裸体视频| 久久综合国产亚洲精品| 联通29元200g的流量卡| 18禁黄网站禁片免费观看直播| 日本一二三区视频观看| 白带黄色成豆腐渣| 午夜亚洲福利在线播放| 国产精品麻豆人妻色哟哟久久 | 网址你懂的国产日韩在线| 久久综合国产亚洲精品| 久久这里有精品视频免费| 国产不卡一卡二| 99在线视频只有这里精品首页| 小说图片视频综合网站| 日韩高清综合在线| 精品国内亚洲2022精品成人| 美女国产视频在线观看| 中文在线观看免费www的网站| 我的女老师完整版在线观看| 午夜精品一区二区三区免费看| 91av网一区二区| 国产色婷婷99| 国产不卡一卡二| kizo精华| 国产成年人精品一区二区| 国产在线男女| 午夜亚洲福利在线播放| 99视频精品全部免费 在线| 国产麻豆成人av免费视频| 精品久久久久久久久亚洲| 最近最新中文字幕大全电影3| 久久久久久国产a免费观看| 国产亚洲av片在线观看秒播厂 | 国产精品,欧美在线| 青春草视频在线免费观看| 97超视频在线观看视频| 天堂av国产一区二区熟女人妻| 一个人观看的视频www高清免费观看| 亚洲第一区二区三区不卡| 精品日产1卡2卡| 亚洲一级一片aⅴ在线观看| 亚洲性久久影院| or卡值多少钱| 亚洲最大成人中文| 在线观看美女被高潮喷水网站| 成人欧美大片| 亚洲国产欧洲综合997久久,| 亚洲av不卡在线观看| 久久久午夜欧美精品| 五月玫瑰六月丁香| 好男人视频免费观看在线| 搞女人的毛片|