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

    環(huán)渤海排污口臨近海域水交換能力研究*

    2014-04-17 07:48:30張宇銘宋朝陽吳克儉張淑芳
    關(guān)鍵詞:小潮大潮排污口

    張宇銘,宋朝陽,吳克儉**,張淑芳

    (1.中國海洋大學(xué)物理海洋實(shí)驗(yàn)室,山東 青島266100;2.國家海洋環(huán)境監(jiān)測中心,遼寧 大連116023)

    渤海作為我國最大的內(nèi)陸海,以有限的海洋資源為沿海地區(qū)創(chuàng)造了巨大的經(jīng)濟(jì)財(cái)富。然而,經(jīng)濟(jì)和社會快速發(fā)展的1個(gè)負(fù)面影響是自然環(huán)境的污染與破壞。前人的研究結(jié)果顯示,渤海沿岸海域污染嚴(yán)重,海洋生態(tài)環(huán)境遭到嚴(yán)重破壞,赤潮等災(zāi)害頻發(fā)[1],而陸源排放,包括排污口排污和入海河流排污等,是渤海污染物的1個(gè)主要來源[2]。在對陸源入海排污口進(jìn)行管理和整治時(shí),如何在不影響經(jīng)濟(jì)發(fā)展的前提下實(shí)現(xiàn)改善海洋環(huán)境的目的成為沿海各級政府部門關(guān)心的重要問題。由于海洋本身有一定的自凈能力,海洋可以通過自身的物理、化學(xué)和生物的作用,使污染物的濃度自然地逐漸降低。因此海洋行政部門可以在排污口臨近海域的海水自凈能力的基礎(chǔ)上探討陸源入海排污口的減排方案,根據(jù)海洋水動(dòng)力條件合理分配排污量和規(guī)劃排污時(shí)間。

    前人對入海排污口臨近海域水交換及物理自凈能力進(jìn)行了一系列的研究。從研究方法的角度大體可分為兩類:一是基于水質(zhì)濃度模型分析在不同的排污量下的污染物擴(kuò)散情況[3-4],計(jì)算環(huán)境容量,以對排污量進(jìn)行合理分配;二是利用質(zhì)點(diǎn)追蹤方法分析入海污染物的輸運(yùn)軌跡[5-6],分析結(jié)果可對排污口選址起到一定的指導(dǎo)作用。然而前人的研究多針對小區(qū)域內(nèi)的特定排污口進(jìn)行分析,如果考慮較大區(qū)域和排污口密集的情況,上述評價(jià)方式將會失效。本文在質(zhì)點(diǎn)追蹤方法的基礎(chǔ)上改進(jìn)了前人的量化評價(jià)指標(biāo),定義了影響距離的概念,在歸并排污口的基礎(chǔ)上分析了環(huán)渤海主要排污口鄰近海域的水交換狀況。若排污口附近水交換狀況良好,則污染物可通過海水的對流擴(kuò)散作用迅速遠(yuǎn)離排污口,與外海水混合,以實(shí)現(xiàn)排污口臨近海域的自凈;反之,污染物將在排污口附近堆積,造成海洋環(huán)境的惡化。因此該研究可為相關(guān)海洋環(huán)境管理部門對各排污口排污量的分配以及未來排污口的選址提供參考。同時(shí)本文探討了季節(jié)和潮汐對排污口臨近海域水交換能力的影響,以期為管理者規(guī)劃設(shè)閘排污口的開閘排污時(shí)間提供一定的參考依據(jù)。

    1 渤海三維水動(dòng)力模型

    本文利用區(qū)域海洋模式ROMS(Regional Oceanic Modeling System)建立渤海三維水動(dòng)力模型。模型的計(jì)算范圍為整個(gè)渤海與部分北黃海,水深數(shù)據(jù)來自渤海海圖,網(wǎng)格分布與地形見圖1。除東邊界為開邊界外,其他均為閉邊界。模型垂向?yàn)镾坐標(biāo)系統(tǒng),均勻分為10層。

    圖1 模型計(jì)算網(wǎng)格與地形Fig.1 Model grid and bathymetry

    計(jì)算區(qū)域的初始條件涉及水位、流速、溫度和鹽度的初始值。由于水位和流速對外界動(dòng)力響應(yīng)較快,初值均取為零;溫度初值和鹽度初值按照《渤海、黃海、東海海洋圖集》[7]的溫鹽分布插值得到。渤海海峽設(shè)置為開邊界,并以水位變化作為模式的驅(qū)動(dòng),開邊界不考慮入流及溫度、鹽度。在本模式中取8個(gè)分潮來計(jì)算水位,包括 M2、S2、K1、O1、K2、N2、P1和 Q1。其中M2、S2、K1、O1分潮的調(diào)和常數(shù)來自《渤海、黃海、東海海洋圖集》,K2、N2、P1和Q1分潮的調(diào)和常數(shù)由華東師范大學(xué)河口海岸國家重點(diǎn)實(shí)驗(yàn)室基于FVCOM建立的東中國海模型的潮汐計(jì)算結(jié)果[8]線性插值得到。模型中采用的海氣熱通量及風(fēng)應(yīng)力驅(qū)動(dòng)數(shù)據(jù)由NCEP再分析資料根據(jù)塊體公式[9]計(jì)算得到。黃河、海河、遼河的總流量占渤海淡水輸入的80%[10],是渤海最主要的3條河流,因此模型中選用這3條河流的流量數(shù)據(jù)作為模型的點(diǎn)源輸入數(shù)據(jù)。河流流量數(shù)據(jù)來自國家水利部發(fā)布的《中國河流泥沙公報(bào)》[11-17],鹽度統(tǒng)一取為1。為了使模式更貼近渤海近幾年的實(shí)際情況,本文采用2001—2010年多年6h平均的海表通量和風(fēng)應(yīng)力數(shù)據(jù)以及2003—2009的多年月平均河流流量數(shù)據(jù)循環(huán)驅(qū)動(dòng)模型,模式積分24個(gè)月達(dá)到穩(wěn)定。垂向采用GLS[18]紊流封閉方程。

    圖2為ROMS模擬得到的M2、S2、K1、O1分潮的等振幅線和同潮時(shí)線,模擬結(jié)果與《渤海、黃海、東海海洋圖集》較為相似。M2分潮振幅在渤海最大可達(dá)1.2m左右,其他分潮振幅均在0.5m以下。

    圖2 渤海 M2(a)、S2(b)、K1(c)、O1(d)分潮等振幅線(紅線;單位:cm)和同潮時(shí)線(藍(lán)線;單位:度)Fig.2 Co-tidal lines(blue line,unit:degree)and iso-amplitude lines(red line,unit:cm)for M2(a),S2(b),K1(c)and O1(d)tidal components in the Bohai Sea

    2010年5~6 月和9~10月分別在金州-普蘭店灣和渤海灣進(jìn)行了海流的半個(gè)月連續(xù)觀測。測流儀器為海床基固定的ADCP,垂向間隔為1m,時(shí)間間隔為1h,站點(diǎn)位置見表1。其中J3、J4分別位于遼寧省附近的普蘭店灣和金州灣,H7位于渤海灣。通過對模擬結(jié)果與大、小潮期間底層監(jiān)測海流數(shù)據(jù)進(jìn)行對比可知(見圖3),該模型可以較好地模擬海流情況。

    2 排污口臨近海域水交換狀況評價(jià)

    2.1 方法介紹

    《2008年渤海海洋環(huán)境質(zhì)量公報(bào)》顯示,本年度國家海洋局共監(jiān)測了96個(gè)環(huán)渤海入海排污口。考慮到環(huán)渤海入海排污口數(shù)量較多,對臨近海域水交換狀況評價(jià)帶來不便,因此首先將環(huán)渤海主要入海排污口按位置歸并至31個(gè)點(diǎn)(見圖4)。各歸并點(diǎn)代表其附近的若干排污口,并將行政邊界作為劃分歸并點(diǎn)的界限,即各歸并點(diǎn)不跨省代表排污口。

    表1 監(jiān)測站點(diǎn)位置Table 1 Observation stations location

    圖3 站位J3、J4、H7底層實(shí)測(線)與模擬(點(diǎn))流速、流向?qū)Ρ菷ig.3 Comparison of observed(line)and calculated(point)current speed and direction of the bottom layer at station J2,J3and H7

    假設(shè)在歸并點(diǎn)附近存在1個(gè)初始高污染區(qū),本文通過Lagrangian質(zhì)點(diǎn)追蹤的方法計(jì)算位于初始高污染區(qū)的質(zhì)點(diǎn)經(jīng)過一段時(shí)間所能到達(dá)的最遠(yuǎn)距離,并將其定義為影響距離。若影響距離較大,證明水交換能力較強(qiáng)。以此來評價(jià)各歸并點(diǎn)所代表的排污口臨近海域的水交換狀況。

    Lagrangian質(zhì)點(diǎn)追蹤方法通過如下方程進(jìn)行質(zhì)點(diǎn)軌跡的求解:

    圖4 歸并點(diǎn)位置Fig.4 The location of the attaching points

    本文為了體現(xiàn)季節(jié)和潮汐對排污口臨近海域水交換狀況的影響,釋放質(zhì)點(diǎn)的時(shí)間定為三維水動(dòng)力模型計(jì)算穩(wěn)定后的夏季(7月)大潮、夏季小潮、冬季(1月)大潮和冬季小潮。由于歸并點(diǎn)分布于渤海沿岸的不同區(qū)域,大、小潮時(shí)間略有差別,所以各歸并點(diǎn)釋放質(zhì)點(diǎn)的時(shí)間起點(diǎn)是根據(jù)數(shù)值模擬得到的該點(diǎn)附近的水位數(shù)據(jù)篩選而定。從選定的時(shí)間起點(diǎn)開始,每2h釋放一次質(zhì)點(diǎn),持續(xù)1d(即13次),這種方式也比較符合多數(shù)排污口連續(xù)排污的特點(diǎn)。由于對排污方式(表層排污還是底層排污)的不確定,本文選取垂向網(wǎng)格的第1、3、5、7、9層釋放質(zhì)點(diǎn)。

    為了定量地評價(jià)排污口臨近海域的水交換能力,本文定義從釋放首個(gè)質(zhì)點(diǎn)的時(shí)刻起30d之內(nèi)全部質(zhì)點(diǎn)所達(dá)的最遠(yuǎn)距離為“30d影響距離”,以符號D30表示。對于1個(gè)指定的排污口,30d影響距離的計(jì)算如下:從釋放首個(gè)質(zhì)點(diǎn)起的30d內(nèi),每天計(jì)算1次從表層到底層的所有質(zhì)點(diǎn)距離排污口的位置,并得到當(dāng)天的質(zhì)點(diǎn)到達(dá)最遠(yuǎn)距離,將這30d的最大距離定義為30d影響距離。

    2.1.1 初始高污染區(qū)范圍的敏感性實(shí)驗(yàn) 在渤海三維水動(dòng)力模型的基礎(chǔ)上,在各歸并點(diǎn)附近選取距離最近的N個(gè)水平網(wǎng)格點(diǎn)代表該歸并點(diǎn)的初始高污染區(qū),通過進(jìn)行敏感實(shí)驗(yàn)選取N值。渤海三維水動(dòng)力模型采用的是正交曲線網(wǎng)格,導(dǎo)致了不同歸并點(diǎn)附近網(wǎng)格密度有所差別,本文采用敏感實(shí)驗(yàn)檢驗(yàn)N值的選取對結(jié)果在多大程度上有影響。隨機(jī)選取了2個(gè)歸并點(diǎn)(8、24號歸并點(diǎn)),在其附近分別選擇距離最近的10、6和3個(gè)網(wǎng)格點(diǎn),并按照上述方式在夏季大潮持續(xù)1d釋放質(zhì)點(diǎn)。結(jié)果顯示(見圖5),初始選取水平網(wǎng)格點(diǎn)的個(gè)數(shù)對于模型計(jì)算30d質(zhì)點(diǎn)所達(dá)到的最大距離影響很小,而模型計(jì)算30d質(zhì)點(diǎn)所達(dá)范圍內(nèi)的質(zhì)點(diǎn)密度和影響范圍略有差別,然而該項(xiàng)指標(biāo)并不影響對排污口臨近海域水交換狀況的評價(jià)。因此綜合考慮計(jì)算資源等問題,本文在水平上選取距離歸并點(diǎn)最近的6個(gè)網(wǎng)格點(diǎn)作為初始高污染區(qū),并在此釋放質(zhì)點(diǎn)。

    圖5 夏季大潮30d影響范圍Fig.5 The 30-day influence distances in the period of spring tide in summer

    2.1.2 質(zhì)點(diǎn)釋放時(shí)間的敏感性實(shí)驗(yàn) 考慮到質(zhì)點(diǎn)追蹤方法存在一定的隨機(jī)性,在同一時(shí)刻、同一位置釋放質(zhì)點(diǎn)可能出現(xiàn)不同的運(yùn)動(dòng)路徑,因此實(shí)驗(yàn)結(jié)果的可重復(fù)性關(guān)系到方法的有效性和結(jié)論的可靠性。本文在模擬夏季大潮期間,進(jìn)行了2次質(zhì)點(diǎn)釋放實(shí)驗(yàn),質(zhì)點(diǎn)釋放方式見上文(即在N=6的初始高污染區(qū)內(nèi)選取垂向5層,并持續(xù)1d釋放質(zhì)點(diǎn))。2次實(shí)驗(yàn)的首個(gè)質(zhì)點(diǎn)釋放時(shí)間均位于大潮最大水位出現(xiàn)附近,兩者相差24h,以評估質(zhì)點(diǎn)釋放時(shí)間對于結(jié)果的影響。經(jīng)過對比,2次實(shí)驗(yàn)中31個(gè)歸并點(diǎn)的30d影響距離幾乎一致。在統(tǒng)計(jì)意義上,2次實(shí)驗(yàn)得到的1~31號歸并點(diǎn)30d影響距離的標(biāo)準(zhǔn)偏差僅為0.46%,說明在1次潮汐過程中,質(zhì)點(diǎn)釋放時(shí)間相差24h不會引起30d影響距離的計(jì)算結(jié)果出現(xiàn)大的偏差,從而證明了該方法對于評估水交換能力的可靠性。

    2.2 模擬結(jié)果分析

    從釋放首個(gè)質(zhì)點(diǎn)的時(shí)刻起30d之內(nèi)全部質(zhì)點(diǎn)所達(dá)的最遠(yuǎn)距離即為30d影響距離D30。對于1個(gè)指定的排污口,30d影響距離的計(jì)算如下:從釋放首個(gè)質(zhì)點(diǎn)起的30d內(nèi),每天計(jì)算一次從表層到底層的所有質(zhì)點(diǎn)距離排污口的位置,并得到當(dāng)天的質(zhì)點(diǎn)到達(dá)最遠(yuǎn)距離,取這30d的最大距離定義為30d影響距離。由于模型計(jì)算區(qū)域有限,對于30d內(nèi)質(zhì)點(diǎn)運(yùn)移已到達(dá)計(jì)算邊界的排污口,其影響距離在表2中以“>”(大于)的形式表示(如29、30、31號歸并點(diǎn))。

    由表2可見,位于遼東灣頂(遼寧省管轄海域)的5號歸并點(diǎn)的平均30d影響距離最小,僅為40.39km,說明該點(diǎn)所代表的排污口臨近海域水交換能力最差,這意味著該排污口的入海污染物極易在排污口臨近海域堆積,造成近岸水環(huán)境的惡化。其他影響距離較小的有位于河北省管轄海域的22號歸并點(diǎn)和位于山東省管轄海域的23號歸并點(diǎn),2個(gè)排污口位于河北與山東交界處。對于這類臨近海域水交換能力差的排污口,應(yīng)采取相應(yīng)的措施限制其排海污染物總量,如提高污水處理技術(shù)等,同時(shí)新建排污口應(yīng)避免在其附近選址。

    平均30d影響距離最大的是位于河北省的9號歸并點(diǎn)。說明該歸并點(diǎn)所代表排污口的污水在排放后可以迅速遠(yuǎn)離排污口,與海水進(jìn)行混合,從而降低排污口附近的污染物濃度。其附近的8、10號歸并點(diǎn)的30d影響距離也較大。然而,影響距離較大并不意味著該排污口可以肆意排污,因?yàn)楹K募{污及自凈能力是有限的。本文僅從水動(dòng)力角度分析水交換能力,但排污量的設(shè)定還需綜合考慮海洋功能區(qū)劃的要求等多方面因素。

    2.2.1 大、小潮對排污口臨近海域水交換能力的影響

    由表2可見,總體上,大潮期所釋放質(zhì)點(diǎn)的影響距離相對于小潮期的影響距離更大,但是夏季有8個(gè)歸并點(diǎn)(1、2、5、14、16、17、25、26)例外,冬季有6個(gè)歸并點(diǎn)(2、5、13、24、25、29)例外。

    從水動(dòng)力角度分析,大潮期釋放的質(zhì)點(diǎn)在最初的幾天時(shí)間運(yùn)移距離較遠(yuǎn)。相對而言,離岸越遠(yuǎn),質(zhì)點(diǎn)本身的運(yùn)動(dòng)能力就越強(qiáng)。因此大潮期釋放的質(zhì)點(diǎn)迅速離岸后獲得更強(qiáng)的運(yùn)動(dòng)能力,從而導(dǎo)致其影響距離相對于小潮期更大。但是由于地形、岸線等其他因素的影響,個(gè)別歸并點(diǎn)也出現(xiàn)了例外的情況。

    表2 30天影響距離Table 2 30-day influence distances /km

    為了表征大、小潮對于30d影響距離的作用,定義了大潮相對于小潮的距離增加率Rtide

    對滿足大潮影響距離大于小潮影響距離的歸并點(diǎn),分別按照上述定義計(jì)算夏季和冬季的大小潮距離增加率Rtide,并對所有歸并點(diǎn)取平均,得到大小潮影響的平均距離增加率為15.8%,意味著大潮期所釋放質(zhì)點(diǎn)的影響距離相對于小潮期所釋放質(zhì)點(diǎn)的影響距離有15.8%的提高。因此,對于滿足該特性的歸并點(diǎn)所代表的排污口可以通過設(shè)閘的方式控制污染物的排放時(shí)間和排放量,以達(dá)到降低污染的目的。

    2.2.2 季節(jié)對排污口臨近海域交換能力的影響 季節(jié)對排污口臨近海域水交換能力的影響主要體現(xiàn)為季風(fēng)對流場的影響。

    將無風(fēng)情況和季風(fēng)作用情況下排污的30d影響范圍進(jìn)行對比可知(見圖6,7),在季風(fēng)和岸線的聯(lián)合作用下,7~24號歸并點(diǎn)均體現(xiàn)為夏季的30d影響范圍位于歸并點(diǎn)偏北側(cè),而冬季的30d影響范圍位于偏南側(cè),這是由于渤海在東亞季風(fēng)的影響下,夏季盛行東南風(fēng),冬季盛行西北風(fēng)。由Ekman理論可知,對于海面邊界層,北半球的風(fēng)引起向右的體積輸送。

    此外,當(dāng)排污口位于海域的上風(fēng)向時(shí),污染物有更大的空間擴(kuò)散,因此影響距離會更大。由圖4可知,位于渤海西北(遼寧及河北沿岸)的6~12號歸并點(diǎn)在冬季處于渤海的上風(fēng)向,而夏季處于渤海的下風(fēng)向,因此6~12號歸并點(diǎn)冬季30天影響距離大于夏季(見表2)。

    圖6 7號(a)、11號(b)、24號(c)和29號(d)歸并點(diǎn)的30d影響范圍Fig.6 The 30-day influence scopes of No.7(a),11(b),24(c)and 29(d)attaching point

    圖7 無風(fēng)情況下7號(a)、11號(b)、24號(c)和29號(d)歸并點(diǎn)的30d影響范圍Fig.7 The 30-day influence scopes of No.7(a),11(b),24(c)and 29(d)attaching point under windless condition

    3 結(jié)語

    本文將環(huán)渤海主要排污口按位置劃分至31個(gè)歸并點(diǎn),在渤海三維水動(dòng)力模型的基礎(chǔ)上,利用質(zhì)點(diǎn)追蹤的方法研究了初始高污染區(qū)的30d影響距離,并以此作為排污口臨近海域水交換能力的衡量標(biāo)準(zhǔn)。若計(jì)算得到的30d影響距離較大,說明該歸并點(diǎn)所代表的排污口的污水在排放后30d運(yùn)移距離較遠(yuǎn),臨近海域水交換能力強(qiáng);反之,則水交換能力弱,大量污染物易聚集,水環(huán)境易惡化。

    通過在夏季大潮、夏季小潮、冬季大潮、冬季小潮4個(gè)時(shí)間段釋放質(zhì)點(diǎn),研究了季節(jié)和潮汐對影響距離的作用??傮w上看,大潮期的30d影響距離大于小潮期,說明大潮期間入海排污口臨近海域的水交換能力更強(qiáng)。管理部門可據(jù)此對排污行為進(jìn)行合理性規(guī)劃,控制排污口的污水排放時(shí)間和排污量,如大潮期間適當(dāng)增加排污量,小潮期間相對減少排污量等,以使污水可以盡快擴(kuò)散,降低污染危害。然而,上述建議只適用于設(shè)置閘門的排污口,而現(xiàn)階段環(huán)渤海入海排污口并非均設(shè)置了閘門。顯然,以設(shè)置閘門的方式對排污行為進(jìn)行控制對海洋環(huán)境保護(hù)有積極意義。因此,首先應(yīng)建議管理部門對入海排污口增設(shè)閘門,進(jìn)而對排污行為合理規(guī)劃。

    環(huán)渤海入海排污口臨近海域水交換能力存在明顯的季節(jié)變化,在東亞季風(fēng)的作用下,各歸并點(diǎn)附近釋放的污染物在擴(kuò)散范圍、方向和距離上于冬季和夏季存在明顯差別。然而本文在季節(jié)分析方面具有一定的局限性,主要原因在于渤海冬季海冰覆蓋率較高,而本文并未考慮海冰的作用,因此冬季的分析結(jié)果可能存在較大偏差。

    [1]于保華,陳瑛.渤海海洋污染與赤潮[J].海洋信息,2000(4):25-26.

    [2]劉學(xué)海.渤海近岸水域環(huán)境污染狀況分析[J].環(huán)境保護(hù)科學(xué),2010,36(1):14-18.

    [3]劉燦生.湄洲灣海域污水排海工程環(huán)境容量評價(jià)的研究[J].哈爾濱建筑大學(xué)學(xué)報(bào),2002,35(4):47-49.

    [4]劉金貴,李瑞杰,曹晶晶,等.蘇北近海水動(dòng)力場及污水稀釋擴(kuò)散的模擬[J].海洋學(xué)研究,2009,27(3):8-13.

    [5]孫英蘭,孫長青,趙可勝.青島東部開發(fā)區(qū)排污口優(yōu)選[J].青島海洋大學(xué)學(xué)報(bào),1994(S1):134-141.

    [6]褚芹芹,李磊,李培良.膠州灣潮流場的示蹤粒子數(shù)值模擬研究[J].中國海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2010,40(11):29-34.

    [7]海洋圖集編委會.渤海、黃海、東海海洋圖集(水文分冊)[M].北京:海洋出版社,1992.

    [8]葛建忠.Multi-scale FVCOM model system for the East China Sea and Changjiang Estuary and its applications[D].上海:華東師范大學(xué),2010.

    [9]Fairall C W,Bradley E F,Rogers D P,et al.Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere Coupled-Ocean Atmosphere Response Experiment [J].J Geophys Res,1996,101:3747-3764.

    [10]Liu H.Fate of three major rivers in the Bohai Sea:A model study[J].Continental Shelf Research,2011,31:1490-1499.

    [11]中華人民共和國水利部.2003年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2003.

    [12]中華人民共和國水利部.2004年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2004.

    [13]中華人民共和國水利部.2005年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2005.

    [14]中華人民共和國水利部.2006年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2006.

    [15]中華人民共和國水利部.2007年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2007.

    [16]中華人民共和國水利部.2008年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2008.

    [17]中華人民共和國水利部.2009年中國河流泥沙公報(bào)[M].北京:中國水利水電出版社,2009.

    [18]Warner J C,Sherwood C R,Arango H G,et al.Performance of four turbulence closure methods implemented using ageneric length scale method[J].Ocean Modelling,2005,8:81-113.

    [19]Hunter J R,Craig P D,Philips H E.On the use of random walk models with spatially variable diffusivity[J].Journal of Computational Physics,1993,106:366-376.

    [20]Visser A W.Using random walk models to simulate the vertical distribution of particles in a turbulent water column[J].Marine Ecology-Progress,1997,158:275-281.

    [21]陳昞睿.質(zhì)點(diǎn)追蹤法及其應(yīng)用[D].上海:華東師范大學(xué),2005.

    [22]Bilgili A,Proehl J A,Lynch D R,et al.Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary:A Lagrangian modeling study[J].Estuarine,Coastal and Shelf Science,2005,65:607-624.

    [23]龍紹橋,婁安剛,譚海濤,等.海上溢油粒子追蹤預(yù)測模型中的兩種數(shù)值方法比較[J].中國海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2006,36(Sup):157-162.

    [24]Pinones A,Hofmann E E,Dinniman M S,et al.Lagrangian sim-ulation of transport pathways and residence times along the western Antarctic Peninsula [J].Deep-Sea Research II,2011,58:1524-1539.

    猜你喜歡
    小潮大潮排污口
    第五課:中國青年投身改革開放的大潮
    直擊入河入海排污口排查:“十四五”持續(xù)改善水環(huán)境質(zhì)量的“牛鼻子”
    5G手機(jī)大潮即來
    海峽姐妹(2019年1期)2019-03-23 02:42:54
    希 望
    希望
    河南省淮河流域入河排污口調(diào)查及布設(shè)研究
    新一季流行色已發(fā)布?快來入手同色系數(shù)碼3C小潮物!
    《神威啟示錄》系列報(bào)道四 大潮退去看“裸泳”
    河南省2014年入河湖排污口數(shù)量
    海河流域入河排污口分布強(qiáng)度分析
    海河水利(2014年3期)2014-02-20 10:02:28
    91九色精品人成在线观看| 一级a爱片免费观看的视频| 久久天躁狠狠躁夜夜2o2o| av有码第一页| 国产在线观看jvid| 亚洲自偷自拍图片 自拍| 桃色一区二区三区在线观看| 18禁黄网站禁片午夜丰满| 女性生殖器流出的白浆| 在线免费观看的www视频| 亚洲美女黄片视频| 一级毛片高清免费大全| 人妻久久中文字幕网| 成人亚洲精品av一区二区| 国产区一区二久久| 久久久国产成人精品二区| 在线观看免费午夜福利视频| 99在线视频只有这里精品首页| 麻豆成人午夜福利视频| 国产熟女午夜一区二区三区| 很黄的视频免费| 夜夜爽天天搞| 日韩大码丰满熟妇| 欧美色视频一区免费| 大香蕉久久成人网| 欧美人与性动交α欧美精品济南到| 国产真实乱freesex| 别揉我奶头~嗯~啊~动态视频| 成年人黄色毛片网站| 国产熟女午夜一区二区三区| 免费在线观看黄色视频的| 免费在线观看黄色视频的| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 精品一区二区三区av网在线观看| 级片在线观看| svipshipincom国产片| 在线播放国产精品三级| 精品日产1卡2卡| 国产v大片淫在线免费观看| 亚洲av电影不卡..在线观看| 日本三级黄在线观看| 国产人伦9x9x在线观看| 亚洲精品av麻豆狂野| 成人欧美大片| 人人妻人人看人人澡| 老司机深夜福利视频在线观看| 国语自产精品视频在线第100页| 免费看日本二区| 黑人欧美特级aaaaaa片| 好男人在线观看高清免费视频 | 操出白浆在线播放| 国产单亲对白刺激| 亚洲第一电影网av| 一区二区三区国产精品乱码| 久久精品亚洲精品国产色婷小说| 国产免费男女视频| 侵犯人妻中文字幕一二三四区| 在线观看www视频免费| 国产精品久久久久久亚洲av鲁大| 啦啦啦免费观看视频1| 天堂√8在线中文| 免费看美女性在线毛片视频| av视频在线观看入口| 国产高清有码在线观看视频 | 1024香蕉在线观看| 久久久久久大精品| 午夜亚洲福利在线播放| 欧美黑人巨大hd| 免费看十八禁软件| 亚洲色图 男人天堂 中文字幕| 国内精品久久久久久久电影| 欧美国产日韩亚洲一区| 亚洲在线自拍视频| 日韩大尺度精品在线看网址| 9191精品国产免费久久| 亚洲国产高清在线一区二区三 | 久久久国产精品麻豆| 满18在线观看网站| 欧美一级毛片孕妇| 国产激情久久老熟女| 精品久久久久久久久久免费视频| 欧美乱码精品一区二区三区| 午夜激情av网站| 一进一出抽搐gif免费好疼| 免费看a级黄色片| 在线免费观看的www视频| 色综合站精品国产| 成人午夜高清在线视频 | 亚洲最大成人中文| 在线观看66精品国产| 成人av一区二区三区在线看| 男人舔奶头视频| 少妇被粗大的猛进出69影院| 啦啦啦 在线观看视频| 亚洲av第一区精品v没综合| 大香蕉久久成人网| aaaaa片日本免费| 国产麻豆成人av免费视频| av片东京热男人的天堂| 色尼玛亚洲综合影院| 窝窝影院91人妻| 亚洲一区中文字幕在线| 久久婷婷成人综合色麻豆| 女人高潮潮喷娇喘18禁视频| 久久婷婷人人爽人人干人人爱| 亚洲精品美女久久久久99蜜臀| or卡值多少钱| 免费看十八禁软件| 一进一出抽搐动态| 国产激情久久老熟女| 看片在线看免费视频| 欧美精品啪啪一区二区三区| 欧美日韩黄片免| 美女午夜性视频免费| АⅤ资源中文在线天堂| 性色av乱码一区二区三区2| 久久九九热精品免费| 久久久久久亚洲精品国产蜜桃av| 精品久久久久久久久久免费视频| 久久精品亚洲精品国产色婷小说| 久久香蕉精品热| 99国产极品粉嫩在线观看| 亚洲激情在线av| av有码第一页| 天天添夜夜摸| 一级片免费观看大全| 午夜精品久久久久久毛片777| 久久久国产成人精品二区| 夜夜躁狠狠躁天天躁| 色婷婷久久久亚洲欧美| 日韩欧美 国产精品| 国产私拍福利视频在线观看| 十八禁网站免费在线| 久久中文字幕一级| 1024视频免费在线观看| 在线看三级毛片| 国产aⅴ精品一区二区三区波| 91九色精品人成在线观看| 动漫黄色视频在线观看| 女人高潮潮喷娇喘18禁视频| 国产免费av片在线观看野外av| 久久久久久大精品| 中国美女看黄片| 亚洲成a人片在线一区二区| 欧美大码av| 国产熟女xx| 少妇粗大呻吟视频| 丰满的人妻完整版| 亚洲精品中文字幕一二三四区| 在线观看免费午夜福利视频| 国产aⅴ精品一区二区三区波| 中文字幕精品免费在线观看视频| 国产一区二区三区视频了| 18禁黄网站禁片免费观看直播| 美女高潮到喷水免费观看| 国产亚洲av嫩草精品影院| 欧美丝袜亚洲另类 | 色哟哟哟哟哟哟| 婷婷六月久久综合丁香| 亚洲av电影在线进入| 99精品久久久久人妻精品| 国产精品久久电影中文字幕| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品第一综合不卡| 久久精品影院6| 满18在线观看网站| 欧美日韩福利视频一区二区| 级片在线观看| 老司机深夜福利视频在线观看| 悠悠久久av| 日韩 欧美 亚洲 中文字幕| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| 欧美激情久久久久久爽电影| 国产精品综合久久久久久久免费| 久久精品夜夜夜夜夜久久蜜豆 | 天天躁狠狠躁夜夜躁狠狠躁| 美女高潮喷水抽搐中文字幕| cao死你这个sao货| 免费在线观看亚洲国产| 国内精品久久久久精免费| 一二三四在线观看免费中文在| 亚洲av成人不卡在线观看播放网| 黄色女人牲交| 99久久久亚洲精品蜜臀av| 精华霜和精华液先用哪个| 国产午夜精品久久久久久| 这个男人来自地球电影免费观看| 久久久久国产一级毛片高清牌| 中文字幕av电影在线播放| 亚洲电影在线观看av| 国产精品久久久久久精品电影 | 一级毛片高清免费大全| 国产99白浆流出| 久久天堂一区二区三区四区| 欧美激情高清一区二区三区| 精品欧美国产一区二区三| 级片在线观看| 激情在线观看视频在线高清| 男女午夜视频在线观看| 欧美黑人欧美精品刺激| 国产精品久久视频播放| 亚洲第一欧美日韩一区二区三区| 人人妻人人看人人澡| 大型黄色视频在线免费观看| 黄色视频不卡| 禁无遮挡网站| 观看免费一级毛片| 国产精品一区二区精品视频观看| 国产av在哪里看| 国产av不卡久久| 久久九九热精品免费| 久久午夜综合久久蜜桃| 老熟妇乱子伦视频在线观看| 免费看a级黄色片| 国产精品九九99| 午夜免费激情av| 欧美国产日韩亚洲一区| 白带黄色成豆腐渣| 国产一区在线观看成人免费| 亚洲av电影在线进入| 亚洲人成电影免费在线| 伦理电影免费视频| 欧美日韩精品网址| 亚洲自拍偷在线| 欧美性猛交╳xxx乱大交人| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| 波多野结衣av一区二区av| 欧美性猛交╳xxx乱大交人| 久久久国产成人免费| 怎么达到女性高潮| 999精品在线视频| 一二三四社区在线视频社区8| 日韩高清综合在线| 特大巨黑吊av在线直播 | 欧美黄色片欧美黄色片| 男人操女人黄网站| 天堂影院成人在线观看| 女性被躁到高潮视频| 国产一区在线观看成人免费| 欧美黑人巨大hd| 国产人伦9x9x在线观看| 中文亚洲av片在线观看爽| 国产私拍福利视频在线观看| 俺也久久电影网| 亚洲欧洲精品一区二区精品久久久| 国产亚洲av高清不卡| 免费在线观看完整版高清| a在线观看视频网站| 精品第一国产精品| 亚洲国产高清在线一区二区三 | 国产精品 国内视频| 亚洲精品国产一区二区精华液| 一本综合久久免费| 人人妻人人澡欧美一区二区| 黄色视频,在线免费观看| x7x7x7水蜜桃| 免费一级毛片在线播放高清视频| 悠悠久久av| 久久香蕉国产精品| 免费在线观看成人毛片| 一区二区三区国产精品乱码| 伊人久久大香线蕉亚洲五| 日韩欧美三级三区| 两性夫妻黄色片| 黄片播放在线免费| 妹子高潮喷水视频| 男女床上黄色一级片免费看| 欧美人与性动交α欧美精品济南到| 国产激情欧美一区二区| 精品电影一区二区在线| 淫秽高清视频在线观看| 美女 人体艺术 gogo| 看黄色毛片网站| 午夜福利18| 国产成人精品久久二区二区91| 国产成人一区二区三区免费视频网站| 亚洲一区高清亚洲精品| 国内毛片毛片毛片毛片毛片| 在线观看66精品国产| 美女 人体艺术 gogo| 天堂动漫精品| 免费在线观看视频国产中文字幕亚洲| 国产精品免费视频内射| 村上凉子中文字幕在线| 美女高潮到喷水免费观看| 人人妻人人澡人人看| 一本综合久久免费| 丁香欧美五月| 欧美黑人巨大hd| 日本一本二区三区精品| 999精品在线视频| a级毛片a级免费在线| 夜夜夜夜夜久久久久| 岛国在线观看网站| 两个人免费观看高清视频| 大香蕉久久成人网| 精品国内亚洲2022精品成人| 亚洲一区二区三区色噜噜| 国产黄a三级三级三级人| 欧美日韩福利视频一区二区| 精品久久久久久久人妻蜜臀av| 亚洲国产毛片av蜜桃av| 欧美在线一区亚洲| 精品久久久久久成人av| 在线观看免费午夜福利视频| 99久久国产精品久久久| 好看av亚洲va欧美ⅴa在| 婷婷六月久久综合丁香| 亚洲五月婷婷丁香| 欧美久久黑人一区二区| 99国产极品粉嫩在线观看| 亚洲黑人精品在线| 91大片在线观看| 最好的美女福利视频网| 欧美成人性av电影在线观看| 亚洲第一电影网av| 不卡一级毛片| 国产午夜精品久久久久久| 免费在线观看亚洲国产| 99国产精品一区二区三区| 草草在线视频免费看| 亚洲五月婷婷丁香| 国产激情欧美一区二区| 黄色 视频免费看| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 亚洲精品一卡2卡三卡4卡5卡| 99精品欧美一区二区三区四区| 欧美日韩福利视频一区二区| 香蕉久久夜色| 久热这里只有精品99| 男人的好看免费观看在线视频 | 国产av一区在线观看免费| 欧美日韩黄片免| 亚洲av电影在线进入| 午夜福利成人在线免费观看| 黄网站色视频无遮挡免费观看| 国产精品一区二区三区四区久久 | 亚洲精品av麻豆狂野| 精品久久久久久,| 搡老妇女老女人老熟妇| 亚洲精品国产一区二区精华液| 日本黄色视频三级网站网址| 亚洲中文av在线| 国产成人欧美在线观看| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 亚洲国产毛片av蜜桃av| www.www免费av| 亚洲 欧美一区二区三区| 女人高潮潮喷娇喘18禁视频| 69av精品久久久久久| 淫妇啪啪啪对白视频| 最好的美女福利视频网| 亚洲黑人精品在线| 午夜免费鲁丝| 午夜福利在线在线| 久久精品成人免费网站| 十分钟在线观看高清视频www| 国产精品国产高清国产av| 国产精品一区二区三区四区久久 | 日韩大尺度精品在线看网址| 97碰自拍视频| а√天堂www在线а√下载| 亚洲欧洲精品一区二区精品久久久| 日本免费一区二区三区高清不卡| 99国产综合亚洲精品| svipshipincom国产片| 亚洲精品在线观看二区| 午夜福利在线观看吧| 亚洲精品av麻豆狂野| 久久国产乱子伦精品免费另类| 一本久久中文字幕| 韩国av一区二区三区四区| 国产精品久久久人人做人人爽| 很黄的视频免费| 侵犯人妻中文字幕一二三四区| 久久久久久久午夜电影| 国产成人啪精品午夜网站| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久末码| 精品卡一卡二卡四卡免费| 草草在线视频免费看| 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区mp4| 久久中文字幕人妻熟女| 亚洲天堂国产精品一区在线| 国产三级在线视频| 99久久无色码亚洲精品果冻| 国内精品久久久久精免费| 丁香欧美五月| 真人一进一出gif抽搐免费| 国产成+人综合+亚洲专区| 变态另类丝袜制服| 欧美中文综合在线视频| 波多野结衣高清无吗| 高清毛片免费观看视频网站| 国产伦一二天堂av在线观看| 天堂√8在线中文| 精品久久久久久成人av| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| av免费在线观看网站| 精品国产乱码久久久久久男人| 91老司机精品| 亚洲成人精品中文字幕电影| 久久久久国内视频| 少妇熟女aⅴ在线视频| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 国产真人三级小视频在线观看| bbb黄色大片| 一二三四在线观看免费中文在| 听说在线观看完整版免费高清| 黄片播放在线免费| 亚洲欧美日韩无卡精品| 中文字幕人妻熟女乱码| 亚洲五月色婷婷综合| 精品免费久久久久久久清纯| 窝窝影院91人妻| 久久天堂一区二区三区四区| 狂野欧美激情性xxxx| 日韩免费av在线播放| 久久天躁狠狠躁夜夜2o2o| 波多野结衣巨乳人妻| 国产私拍福利视频在线观看| 精品久久蜜臀av无| 波多野结衣高清无吗| 后天国语完整版免费观看| 日韩精品中文字幕看吧| 久久人妻福利社区极品人妻图片| 人人澡人人妻人| 国产高清激情床上av| 亚洲免费av在线视频| 欧美亚洲日本最大视频资源| 久久久久国内视频| 亚洲国产精品合色在线| 国产成+人综合+亚洲专区| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 国产亚洲精品av在线| 国产精品九九99| 50天的宝宝边吃奶边哭怎么回事| 亚洲男人天堂网一区| 亚洲第一电影网av| 久久狼人影院| 国产v大片淫在线免费观看| 亚洲专区中文字幕在线| 欧美成人一区二区免费高清观看 | 操出白浆在线播放| 国产一级毛片七仙女欲春2 | 欧美精品啪啪一区二区三区| a级毛片a级免费在线| 50天的宝宝边吃奶边哭怎么回事| 在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 久久香蕉精品热| 亚洲七黄色美女视频| 99久久99久久久精品蜜桃| 美女高潮到喷水免费观看| 观看免费一级毛片| 亚洲精品久久成人aⅴ小说| 老司机午夜福利在线观看视频| 男人舔奶头视频| 麻豆成人午夜福利视频| 国产亚洲精品久久久久久毛片| 久久精品国产亚洲av高清一级| 波多野结衣巨乳人妻| 亚洲久久久国产精品| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| 99久久国产精品久久久| 身体一侧抽搐| 国产成人影院久久av| 国产精品免费视频内射| 亚洲一码二码三码区别大吗| av电影中文网址| 亚洲久久久国产精品| 久久国产精品人妻蜜桃| 欧美黑人精品巨大| 欧美一级毛片孕妇| 色在线成人网| 亚洲人成伊人成综合网2020| 日本免费一区二区三区高清不卡| 一本久久中文字幕| 老司机在亚洲福利影院| 精品免费久久久久久久清纯| 精品卡一卡二卡四卡免费| 国产一区二区激情短视频| 黄色a级毛片大全视频| 在线免费观看的www视频| 国产激情久久老熟女| www日本黄色视频网| 在线av久久热| 欧美+亚洲+日韩+国产| 欧美日韩福利视频一区二区| 国产99白浆流出| 久久精品成人免费网站| 精品少妇一区二区三区视频日本电影| 两个人视频免费观看高清| 男男h啪啪无遮挡| 在线免费观看的www视频| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av香蕉五月| 精品一区二区三区四区五区乱码| 色哟哟哟哟哟哟| 桃红色精品国产亚洲av| 国产亚洲精品久久久久久毛片| 日韩高清综合在线| 亚洲三区欧美一区| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看| 久久久久九九精品影院| 亚洲中文字幕日韩| www.www免费av| 亚洲熟女毛片儿| 高清毛片免费观看视频网站| 午夜两性在线视频| 一本久久中文字幕| 一区二区三区激情视频| 高清毛片免费观看视频网站| 午夜两性在线视频| 午夜免费激情av| 亚洲 欧美 日韩 在线 免费| 99riav亚洲国产免费| 欧美成人午夜精品| 一边摸一边抽搐一进一小说| 在线观看66精品国产| 色精品久久人妻99蜜桃| 国产精品98久久久久久宅男小说| 久久亚洲真实| 欧美精品啪啪一区二区三区| 日本a在线网址| 婷婷亚洲欧美| 99国产精品一区二区蜜桃av| 又黄又爽又免费观看的视频| 夜夜爽天天搞| 国产精品免费视频内射| av超薄肉色丝袜交足视频| 国产精品国产高清国产av| 99久久无色码亚洲精品果冻| 国产精品99久久99久久久不卡| 99久久无色码亚洲精品果冻| 精品不卡国产一区二区三区| 在线观看免费日韩欧美大片| 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 自线自在国产av| 桃色一区二区三区在线观看| 亚洲国产精品999在线| 亚洲第一欧美日韩一区二区三区| 欧美在线黄色| 成人亚洲精品av一区二区| 日本 av在线| 久久午夜综合久久蜜桃| 精品乱码久久久久久99久播| 久久青草综合色| 久久久水蜜桃国产精品网| 一级毛片精品| 美女高潮到喷水免费观看| 久久精品国产综合久久久| 精品国产乱码久久久久久男人| 欧美激情久久久久久爽电影| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 久久九九热精品免费| 久久久久久人人人人人| 国产国语露脸激情在线看| av天堂在线播放| 精品久久蜜臀av无| 亚洲片人在线观看| svipshipincom国产片| 国产av一区二区精品久久| 成人18禁高潮啪啪吃奶动态图| 天堂√8在线中文| 看黄色毛片网站| 国产伦人伦偷精品视频| 激情在线观看视频在线高清| 日韩高清综合在线| 久久久久久九九精品二区国产 | 亚洲成人久久爱视频| 亚洲全国av大片| 国产单亲对白刺激| 草草在线视频免费看| 亚洲人成网站高清观看| 国产av不卡久久| 精品卡一卡二卡四卡免费| 精品免费久久久久久久清纯| 婷婷亚洲欧美| 欧美黑人巨大hd| 日韩中文字幕欧美一区二区| 亚洲成av人片免费观看| 少妇粗大呻吟视频| 黄色片一级片一级黄色片| 波多野结衣av一区二区av| 国产精品久久电影中文字幕| 免费电影在线观看免费观看| 丁香六月欧美| 性色av乱码一区二区三区2| 老熟妇乱子伦视频在线观看| 91国产中文字幕| 好男人电影高清在线观看| 国产视频内射| 久久精品夜夜夜夜夜久久蜜豆 | 99久久国产精品久久久| 精品久久久久久,| 日本在线视频免费播放|