達(dá)標(biāo)污水離岸排海末端處置技術(shù)研究綜述

彭士濤 ，王心海

(1.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，天津 300456;2.國(guó)家海洋局第一海洋研究所，青島 266061)

根據(jù)2010年全國(guó)環(huán)境統(tǒng)計(jì)公報(bào)［1］，我國(guó)“十一五”污水排放量已達(dá)到617億t，預(yù)計(jì)到“十二五”末污水排放量將超過700億t，這意味著我國(guó)將有大量污水進(jìn)入自然水域，嚴(yán)重影響水資源環(huán)境，污水的出路問題已成為制約國(guó)家經(jīng)濟(jì)發(fā)展的重要因素。達(dá)標(biāo)污水采用離岸排海，通過放流系統(tǒng)排放到離岸較遠(yuǎn)的海域，一方面利用海洋強(qiáng)大而豐富的物理、化學(xué)和生物能力稀釋凈化，彌補(bǔ)了陸域污水處理投資大、占地多、運(yùn)行費(fèi)用高等的不足，大量工程應(yīng)用表明離岸排海比二級(jí)處理廠費(fèi)用降低將近一半;另一方面，通過污水排?？谖恢?、擴(kuò)散器射流參數(shù)及結(jié)構(gòu)的優(yōu)化，污水排海的初始稀釋度可達(dá)到100倍左右，相比二級(jí)處理廠的稀釋倍數(shù)將近增加了10倍，而且排海后污染物被快速帶離海岸，對(duì)周邊海岸的環(huán)境影響較小，具有低成本、高處理率等優(yōu)點(diǎn)。

達(dá)標(biāo)污水離岸排海處置技術(shù)研究最早可追溯到19世紀(jì)末英國(guó)的 Harwich污水排海工程，經(jīng)歷了從污水不進(jìn)行處理到污水進(jìn)行一級(jí)、二級(jí)預(yù)處理后排放，直到現(xiàn)在廣泛采用的從近海延伸到深?？梢詽M足一級(jí)處理深海排放的處置技術(shù)［2］。國(guó)內(nèi)沿海城市也有采用，但相關(guān)研究起步較晚，對(duì)于直接關(guān)系到污水排海后的稀釋擴(kuò)散輸移效果的末端排放口位置選劃、擴(kuò)散器射流參數(shù)及水力結(jié)構(gòu)優(yōu)化的研究報(bào)道的相對(duì)較少，且目前國(guó)內(nèi)還沒有擴(kuò)散器設(shè)計(jì)的規(guī)范可以參考，制約了污水離岸處置工程技術(shù)的發(fā)展與應(yīng)用［3-4］。針對(duì)現(xiàn)狀，文中對(duì)國(guó)內(nèi)外達(dá)標(biāo)污水離岸排海口位置選劃、擴(kuò)散器射流參數(shù)及水力結(jié)構(gòu)優(yōu)化等技術(shù)的研究方法及存在問題進(jìn)行了分析，便于后續(xù)研究的參考。

1 排污口位置選劃

排污口位置選劃是離岸排海末端處置的起步工程，直接決定著排海處置效果、工程安全性、經(jīng)濟(jì)性等。目前，排污口位置選劃采用先定性選劃后定量研究，主要采用建立海域水動(dòng)力模型、污染物輸移擴(kuò)散模型、質(zhì)點(diǎn)追蹤模型、泥沙場(chǎng)模型等數(shù)值模型和物理模型的方法［2-18］。早期主要采用 ROBERTS、BROOKS等經(jīng)驗(yàn)公式預(yù)測(cè)近區(qū)稀釋度和海域輸移擴(kuò)散的方法進(jìn)行排污口位置選劃，由于公式預(yù)測(cè)精度不足，目前已很少使用［2］。Blumberg［5］首次采用遠(yuǎn)區(qū)環(huán)流數(shù)值模型對(duì)污水排海稀釋度進(jìn)行了模擬，該數(shù)值模型較多用于海域流場(chǎng)和污染物輸移擴(kuò)散的研究，用于羽流的研究還需進(jìn)一步論證。Etemadshahidi［6］采用CORMIX2和VISJET模型對(duì)排海羽流的上升高度和稀釋度進(jìn)行了模擬，這些特征參數(shù)的確定可以為排污口位置確定提供依據(jù)。Li等［7-9］通過采用近遠(yuǎn)區(qū)復(fù)合模型模擬了污水排海的稀釋擴(kuò)散，預(yù)測(cè)了污水排海的輸移擴(kuò)散和羽流紊動(dòng)稀釋，這將是污染物稀釋擴(kuò)散研究的重要方向。孫長(zhǎng)青等［10-11］考慮到公式選劃的局限性，首次通過建立水質(zhì)二維擴(kuò)散預(yù)測(cè)模型及拉格朗日余流模型，對(duì)污染物輸移擴(kuò)散、濃度變化及質(zhì)點(diǎn)跟蹤進(jìn)行了研究，預(yù)測(cè)了污染物影響包絡(luò)面積和質(zhì)點(diǎn)遷移，并提出先定性后定量的方法，大大節(jié)約了人力物力投入，對(duì)后續(xù)研究起到重要的指導(dǎo)作用。趙俊杰［12］采用MIKE軟件的潮流模型和水動(dòng)力模型模擬分析入海COD的稀釋擴(kuò)散濃度分布并綜合考慮工程施工和經(jīng)濟(jì)性等進(jìn)行排污口位置選劃。侯?。?3］針對(duì)二維模型難以很好適應(yīng)大水深的表層底層流速相差大的問題，提出采用三維水動(dòng)力和水質(zhì)模型更好的模擬了污水排海不同水深的輸移擴(kuò)散規(guī)律。宋強(qiáng)等［14-15］通過研究過渡區(qū)模擬方法，嘗試將近區(qū)噴口射流模擬模型嵌入遠(yuǎn)區(qū)擴(kuò)散模型，形成全場(chǎng)數(shù)學(xué)模型，更真實(shí)的反應(yīng)污水排海的稀釋變化，但后續(xù)相關(guān)研究較少，嵌入方法還有待深入研究。楊樹森等［16］鑒于排海工程應(yīng)用中經(jīng)常發(fā)生噴口淤埋導(dǎo)致排污不正常運(yùn)行問題，通過建立泥沙數(shù)學(xué)模型選取地形穩(wěn)定、沖刷變化小排污口，為選取穩(wěn)定安全的排污口提供了重要依據(jù)。

前人采用先定性后定量的方法，從單一的遠(yuǎn)區(qū)污染物輸移擴(kuò)散模擬到綜合近遠(yuǎn)區(qū)的全場(chǎng)模擬，甚至到考慮排污口安全穩(wěn)定出發(fā)的泥沙沖淤研究，但是這些研究多是就某個(gè)方面的分析，筆者［17-18］通過對(duì)區(qū)域功能符合性、海域稀釋擴(kuò)散、泥沙沖淤、生態(tài)環(huán)境影響、管道工程技術(shù)經(jīng)濟(jì)等多因素綜合分析，選取了更加符合區(qū)域發(fā)展需求和環(huán)境安全要求的排污口。由此，綜合考慮水動(dòng)力、遷移軌跡、稀釋擴(kuò)散、泥沙沖淤，環(huán)境區(qū)域符合性，生態(tài)影響，工程經(jīng)濟(jì)性等多個(gè)角度多方面的因素分析方法是今后排污口科學(xué)選劃的趨勢(shì)。

但總體而言，目前排污口選劃還存在以下問題需要解決:一是排污口選劃多以海域水動(dòng)力條件和稀釋擴(kuò)散預(yù)測(cè)為主，而實(shí)際還與生態(tài)敏感區(qū)等區(qū)域規(guī)劃、工程沖刷穩(wěn)定和運(yùn)行投資密切相關(guān)，應(yīng)綜合加以考慮;二是選劃多以物理過程模擬為主，而實(shí)際中污水排放經(jīng)過幾個(gè)潮周后已經(jīng)被完全稀釋擴(kuò)散，此后的生物作用顯得更為重要，應(yīng)分析污水對(duì)水域生物造成的影響，并重點(diǎn)分析長(zhǎng)期排污可能產(chǎn)生的生物累積效應(yīng);三是現(xiàn)有排海限制指標(biāo)要求少，只對(duì)初始稀釋度、混合區(qū)面積、水深做了簡(jiǎn)單的要求，與實(shí)際排海工程有較大差別，急需結(jié)合污水的排放量和類別制定相應(yīng)的規(guī)范和標(biāo)準(zhǔn)，指導(dǎo)離岸排海工程的健康發(fā)展。

2 擴(kuò)散器射流參數(shù)優(yōu)化

離岸排海末端擴(kuò)散器射流結(jié)構(gòu)決定著擴(kuò)散器百米區(qū)域的紊動(dòng)稀釋和高濃度場(chǎng)的存在與否，是排海處置工程實(shí)現(xiàn)環(huán)保達(dá)標(biāo)排放的保障。擴(kuò)散器射流參數(shù)主要有擴(kuò)散管長(zhǎng)度、上升管個(gè)數(shù)、間距、噴口角度、噴口型式等，通過優(yōu)化這些參數(shù)可使初始稀釋倍數(shù)達(dá)到幾百倍，對(duì)污水排海工程建設(shè)的成敗起著決定性作用。目前，對(duì)于擴(kuò)散器射流參數(shù)優(yōu)化主要采用經(jīng)驗(yàn)公式、物理模型試驗(yàn)、射流數(shù)值模擬、量綱分析以及基礎(chǔ)理論分析相結(jié)合的研究方法［2-4，19-35］。擴(kuò)散器射流參數(shù)優(yōu)化早期多以經(jīng)驗(yàn)公式為主，采用ROBERTS等稀釋度預(yù)測(cè)經(jīng)驗(yàn)公式，提出了影響擴(kuò)散器長(zhǎng)度的排放量、稀釋度、水深、水動(dòng)力等主要因素，推出擴(kuò)散器長(zhǎng)度和起始稀釋度的經(jīng)驗(yàn)公式，在初期擴(kuò)散器射流結(jié)構(gòu)簡(jiǎn)易設(shè)計(jì)中應(yīng)用較多［2-4］。Duer Michael［19］通過物理模型研究了鴨嘴閥稀釋特性，并發(fā)現(xiàn)擴(kuò)散器通過使用鴨嘴閥能取得更好的稀釋效果，保證了鴨嘴閥在擴(kuò)散器應(yīng)用中的可行性。Wang［20］采用DPIV和 PLIF技術(shù)進(jìn)行了擴(kuò)散器周邊區(qū)域的濃度場(chǎng)測(cè)量分析，得到了浮射流的稀釋度，為浮射流研究提供了重要的研究手段。Tian等［21-23］通過建立擴(kuò)散器物理模型，采用3DLIF裝置獲取了擴(kuò)散器排入密度均勻流、密度分層流海域后的污染物稀釋擴(kuò)散圖像，并通過圖像數(shù)據(jù)處理取得了稀釋度，該種方法成為了近區(qū)擴(kuò)散器射流研究的重要方法。Lai等［24-25］采用水槽試驗(yàn)和VISJET數(shù)值模型對(duì)新型的玫瑰型擴(kuò)散器的射流擴(kuò)散稀釋特性進(jìn)行了研究。Patalano［26］采用 ADV物理測(cè)量技術(shù)和 Delft3D數(shù)值模型研究了擴(kuò)散器羽流的混合和稀釋度特性。徐高田等［27-28］通過物理模型試驗(yàn)結(jié)合量綱分析方法得出半理論半經(jīng)驗(yàn)公式，比純經(jīng)驗(yàn)公式有了進(jìn)步，并初步揭示了環(huán)境橫流和噴口射流相互作用機(jī)理，對(duì)擴(kuò)散器結(jié)構(gòu)優(yōu)化具有一定指導(dǎo)意義。唐洪武等［29］運(yùn)用超聲流速儀對(duì)鴨嘴閥射流特性進(jìn)行了系列試驗(yàn)研究，并發(fā)現(xiàn)鴨嘴閥射流近區(qū)紊動(dòng)強(qiáng)度明顯高于圓射流，有利于達(dá)標(biāo)污水稀釋擴(kuò)散，為鴨嘴閥在離岸排海的應(yīng)用上提供了基礎(chǔ)。周連偉等［30］將遠(yuǎn)區(qū)水動(dòng)力結(jié)果引入U(xiǎn)M模型和JETLAG3模型計(jì)算初始稀釋度，對(duì)擴(kuò)散器參數(shù)進(jìn)行優(yōu)化，具有一定的實(shí)用性，但數(shù)模仍難以真實(shí)模擬環(huán)境射流稀釋狀況。方神光［31］提出了靜水中有阻力的浮力射流，貼壁浮力射流，單孔卷吸設(shè)置，射流預(yù)稀釋裝置等增加擴(kuò)散器稀釋度的若干措施，為擴(kuò)散器改進(jìn)提供了理論基礎(chǔ)，研究發(fā)現(xiàn)通過阻力裝置能成倍提高稀釋度，尤其能改善表層的污水場(chǎng)狀況。高柱［32］采用FLUENT軟件建立的同向流動(dòng)環(huán)境中的橢圓射流數(shù)學(xué)模型和LDV技術(shù)量測(cè)射流速度場(chǎng)的方法，對(duì)擴(kuò)散器排海浮射流進(jìn)行了研究和驗(yàn)證。肖洋［33-34］除了采用FLUENT數(shù)模和LDV技術(shù)，還采用PIV測(cè)速技術(shù)對(duì)射流進(jìn)行了研究，對(duì)橫向流動(dòng)條件下多孔水平動(dòng)量射流的濃度場(chǎng)、速度場(chǎng)、渦量場(chǎng)和壓強(qiáng)場(chǎng)的特性進(jìn)行了較為全面的研究。

由于多孔擴(kuò)散器深海排放稀釋擴(kuò)散非常復(fù)雜，受到排污條件、環(huán)境條件和擴(kuò)散器參數(shù)等多種因素的影響，因此，只依據(jù)經(jīng)驗(yàn)公式或數(shù)學(xué)模擬還很難真正解決，物理模型仍然是最重要的研究方法，筆者［18，35］鑒于擴(kuò)散器長(zhǎng)度設(shè)計(jì)主要依賴經(jīng)驗(yàn)公式的現(xiàn)狀，在大量水槽試驗(yàn)研究的基礎(chǔ)上，提出了數(shù)模物模相結(jié)合確定擴(kuò)散器長(zhǎng)度的方法，較為準(zhǔn)確的提出了擴(kuò)散器的有效長(zhǎng)度，取得了很好的效果，圖1為天津某排海工程通過三維數(shù)值模擬的帶有10，20個(gè)和30個(gè)擴(kuò)散器上升管方案的污染物表層濃度分布，圖中明顯看到10、20、30個(gè)方案的表層最大污染物濃度分別為0.0098、0.0083 mg/L 和0.006 mg/L，10 個(gè)方案明顯存在高于臨界濃度0.0083 mg/L的深色超標(biāo)區(qū)域，不能滿足擴(kuò)散器排污表層稀釋濃度要求，20個(gè)方案介于臨界濃度0.0083 mg/L之中，30個(gè)方案明顯小于臨界濃度0.0083 mg/L，鑒于上升管建設(shè)的高投資問題，通過建立水槽物理模型試驗(yàn)?zāi)M10個(gè)和20個(gè)方案浮射流作用下的擴(kuò)散器排污濃度分布，并用分光光度計(jì)進(jìn)行了近區(qū)濃度測(cè)量，發(fā)現(xiàn)10個(gè)方案仍存在大于0.0083 mg/L的超標(biāo)濃度，20個(gè)方案的超標(biāo)濃度消失，且最大濃度僅為0.0071 mg/L，綜合考慮，認(rèn)為擴(kuò)散器采用20個(gè)上升管方案可以滿足環(huán)保要求。水質(zhì)數(shù)值模擬和水槽物理模型相結(jié)合的方法可以進(jìn)一步優(yōu)化了擴(kuò)散器的射流結(jié)構(gòu)，并確保研究結(jié)論的準(zhǔn)確性。

圖1 三維數(shù)值模擬污染物表層濃度分布Fig.1 the surface concentration field of pollutants by three-dimensional numerical simulation

總之，對(duì)于擴(kuò)散器長(zhǎng)度、上升管間距、噴口角度等已有一些研究成果，但這些成果多針對(duì)某個(gè)工程進(jìn)行，不具有普遍性，沒有形成相應(yīng)的規(guī)范和標(biāo)準(zhǔn)，不利于很好的推廣;此外，擴(kuò)散器經(jīng)歷了由單孔排放到設(shè)置多個(gè)出口，再到設(shè)置上升管及多個(gè)噴口，每一步改進(jìn)都使污水離岸排海處理效果大幅度提高，海域環(huán)境得到更好的保護(hù)利用，但目前對(duì)于工程較多采用的多噴口型式的射流理論研究相對(duì)尚少，需要相關(guān)的研究為工程建設(shè)提供更多的技術(shù)支持。

3 擴(kuò)散器水力結(jié)構(gòu)優(yōu)化

離岸排海擴(kuò)散器水力結(jié)構(gòu)是整個(gè)排海工程正常運(yùn)行的關(guān)鍵。隨著排海工程的大規(guī)模應(yīng)用，排放量和排污口水深的要求也越來越大，擴(kuò)散器出現(xiàn)噴口出流不均勻甚至部分管道沒有出流、泥沙生物入侵管道并在管道內(nèi)沉積嚴(yán)重等現(xiàn)象，都可能導(dǎo)致排海工程運(yùn)行差甚至無法運(yùn)行。擴(kuò)散器水力結(jié)構(gòu)優(yōu)化主要通過對(duì)復(fù)雜管道的管徑設(shè)計(jì)和添加附加設(shè)置，使得污水能在各個(gè)出口均勻出流，同時(shí)避免海水倒灌和泥沙生物入侵現(xiàn)象的發(fā)生，為污水排海提供穩(wěn)定可行的運(yùn)行環(huán)境。主要包括兩方面內(nèi)容:一是擴(kuò)散器流量分配的設(shè)計(jì)，實(shí)現(xiàn)流量按擴(kuò)散器具體工程要求進(jìn)行沿程分配，使污水排出后與受納水體進(jìn)行最大程度的摻混稀釋;二是海水泥沙入侵的防止及清除設(shè)計(jì)，海水密度高于污水密度，極易發(fā)生海水倒灌擴(kuò)散器系統(tǒng)，嚴(yán)重影響擴(kuò)散器排放量，腐蝕排海管道內(nèi)壁，同時(shí)攜帶的泥沙還會(huì)造成擴(kuò)散器內(nèi)泥沙長(zhǎng)期淤積，發(fā)生擴(kuò)散器堵塞。這些問題的解決需要了解擴(kuò)散器的內(nèi)部水流運(yùn)動(dòng)規(guī)律，并對(duì)影響因素有深刻的理解，才能實(shí)現(xiàn)擴(kuò)散器水力結(jié)構(gòu)的針對(duì)性設(shè)計(jì)。

擴(kuò)散器水力結(jié)構(gòu)優(yōu)化中管道流量分配水力參數(shù)已有較多的研究報(bào)道［2-4，36-42］。Bleninger 等［36-39］通過采用CFD軟件模擬了擴(kuò)散器內(nèi)部流場(chǎng)和壓力場(chǎng)，得出了擴(kuò)散器的出口流量分布均勻度、水頭損失以及局部阻力系數(shù)等。賴翼峰等［40-42］采用出流系數(shù)法和動(dòng)壓頭法兩種水力方法進(jìn)行了末端處置擴(kuò)散器的流場(chǎng)計(jì)算，并通過物理模型試驗(yàn)驗(yàn)證，這是擴(kuò)散器流量分配水力設(shè)計(jì)的主要方法。

目前，筆者［18］采用FLUNENT軟件建立擴(kuò)散器管道水力模型，并通過物理模型驗(yàn)證得出很好的模擬結(jié)果，F(xiàn)LUNENT軟件的應(yīng)用也大大提高了流量分配設(shè)計(jì)便捷性和準(zhǔn)確性。圖2為天津某排海工程數(shù)值模擬擴(kuò)散器管內(nèi)流速分布圖，圖中可見最大流速出現(xiàn)在擴(kuò)散器噴口處，其流速為3 m/s左右，主管內(nèi)流速為2 m/s左右，滿足擴(kuò)散器噴口和主管流速的要求，通過數(shù)值模型進(jìn)行擴(kuò)散器水力結(jié)構(gòu)方案比選，能有效減少物理模型的試驗(yàn)方案。在此基礎(chǔ)上，通過建立擴(kuò)散器管道水力物理模型，采用體積流量測(cè)量法和測(cè)壓管法分別得到了噴口的流量，計(jì)算得到了擴(kuò)散器的阻力系數(shù)和管道流速，發(fā)現(xiàn)數(shù)值模擬和物理模擬的誤差為1%—5%，進(jìn)一步驗(yàn)證細(xì)化了數(shù)值模擬結(jié)果，說明數(shù)模和物理模型相結(jié)合的方法可有效的用于擴(kuò)散器結(jié)構(gòu)優(yōu)化。

圖2 數(shù)值模擬的擴(kuò)散器速度流場(chǎng)分布圖Fig.2 Velocity field of the diffuser by FLUENT numerical simulation

曾小清［43］針對(duì)排海工程經(jīng)常在非設(shè)計(jì)工況下運(yùn)行易發(fā)生海水入侵的問題，提出了海水入侵沖洗的臨界水力特點(diǎn)和辨別式，為非設(shè)計(jì)條件的運(yùn)行管理提供依據(jù)。劉成等［44］對(duì)排海工程泥沙入侵問題進(jìn)行了系統(tǒng)研究，提出了改進(jìn)尾部和定期沖洗來降低改善淤積的措施，并提出臨界入侵、沖洗密度數(shù)公式，并提出在控制密度傅汝德數(shù)的同時(shí)增設(shè)文丘里管和安裝鴨嘴閥等改進(jìn)措施，很好的改善了海水泥沙入侵的問題。錢達(dá)仁等［45-46］研究了波浪對(duì)管道出流的影響，發(fā)現(xiàn)小流量時(shí)管道出流明顯受波浪影響，極易發(fā)生海水入侵，應(yīng)引起后續(xù)研究的重視。綜上，海水泥沙入侵?jǐn)U散器的主要因素是密度差大和波浪變化，防止措施主要是通過優(yōu)化噴口參數(shù)增大傅汝徳數(shù)來防止擴(kuò)散器入侵，同時(shí)噴口采用安裝鴨嘴閥也能較好的解決擴(kuò)散器入侵問題。

目前研究中，擴(kuò)散器水力參數(shù)優(yōu)化還存在著海水泥沙入侵機(jī)理及泥沙入侵清除技術(shù)不夠完善，這是對(duì)于擴(kuò)散器大量推廣應(yīng)用的一個(gè)重要制約因素，有待深入研究。

4 結(jié)論

通過闡述國(guó)內(nèi)外離岸排污口位置選劃、擴(kuò)散器射流參數(shù)及水力結(jié)構(gòu)優(yōu)化等技術(shù)的發(fā)展現(xiàn)狀，分析了達(dá)標(biāo)污水離岸排海末端處置技術(shù)中存在的問題及相應(yīng)的解決方法:

(1)離岸排海擴(kuò)散器位置選劃作為排海的起始工程，決定著整個(gè)工程的成敗。綜合考慮水動(dòng)力、遷移軌跡、稀釋擴(kuò)散、泥沙沖淤，環(huán)境區(qū)域符合性，生態(tài)影響，工程經(jīng)濟(jì)性等多個(gè)角度多方面的因素分析方法是今后排污口科學(xué)選劃的趨勢(shì)。而污水對(duì)生物影響的累積效應(yīng)的研究是今后需要不斷豐富和突破的研究?jī)?nèi)容。

(2)擴(kuò)散器射流參數(shù)優(yōu)化對(duì)離岸排海工程建設(shè)至關(guān)重要，需要通過擴(kuò)散器射流結(jié)構(gòu)和射流理論，尤其是噴口浮射流理論的系統(tǒng)研究，進(jìn)一步完善規(guī)范擴(kuò)散器射流參數(shù)的優(yōu)化技術(shù)和型式設(shè)計(jì)。

(3)擴(kuò)散器管道水流理論基礎(chǔ)研究尚少，海水泥沙入侵清除機(jī)理和海泥沙入侵清除技術(shù)還不夠完善，成為制約擴(kuò)散器推廣應(yīng)用的重要因素，有待深入研究。

(4)目前為止還沒有離岸排海擴(kuò)散器設(shè)計(jì)的規(guī)范可以參考，設(shè)計(jì)中隨意性比較大，設(shè)計(jì)規(guī)范化、標(biāo)準(zhǔn)化成為發(fā)展污水離岸排海技術(shù)急待解決的問題。

［1］China Environment Statistics Bulletin. Beijing: Ministry of Environmental Protection of the People's Republic of China，2010.

［2］Wood I R，Bell R G，Willkinson D L.Ocean Disposal of Wastewater.Singapore，Incorporated:World Scientific Publishing Company，1993:0-425.

［3］Zhang Y L，Yan H B.Technical Guidelines of Wastewater Marine Disposal.Beijing:China Environmental Science Press，1995.

［4］Xia Q.Engineering Design Theory and Method of Wastewater Marine Disposal.Beijing: China Environmental Science Press，1996.

［5］Blumberg A F，Ji Z G，Ziegler C K.Modeling outfall plume behavior using far field circulation model.Journal of Hydraulic Engineering，1996，122(11):610-616.

［6］Etemad-Shahidi A，Azimi A.Testing the CORMIX2 and VISJET models to predict the dilution of Sanfrancisco outfall //Proceedings of the 7th IWA International Conference on Diffuse Pollution and Basin Management.Dublin: Diffuse Pollution Conference，2003:129-133.

［7］Li S，Hodgins D O.A dynamically coupled outfall plumecirculation model for effluent dispersion in Burrard Inlet，British Columbia.Journal of Environmental Engineering and Science，2004，3(5):433-449.

［8］Bleninger T，Jirka G H.Near-and far-field model coupling methodology for wastewater discharges//Lee J H W，Lam K M.Environmental Hydraulics and Sustainable Water Management.London:Balkema，2004:447-453.

［9］Morelissen R，van der Kaaij T，van Vossen B，Bleninger T.Dynamic coupling of near-field and far-field models for simulating effluent discharges.Water Science Technology，2013，67(10):2210-2220.

［10］Sun C Q，Zhao K S，Sun Y L.Selection of pollutant discharge outlet and prediction of water quality in the offshore areaof Basuo.Journal of Ocean University of Qingdao，1994，(S1):209-216.

［11］Yin Y，Zhong WN，Chang N H，Li W，Xu GS.Study of sewage draining schemes Qinhuangdao city in Beidai River Area.Coastal Engineering，1997，16(2):6-12.

［12］Zhao J J，Bai J，Kang S H.A modeling study on selecting the effluent discharge points of the Nantong seawage discharging project.Coastal Engineering，2012，31(2):31-38.

［13］Hou P.Three Dimensional Mathematical Model of Coastal Sewage Ocean Discharge in Small Yangkou Port［D］.Nanjing:Hohai University，2008.

［14］Song Q，Wu H.Research on the marine environmental impact from Shenzhen municipal sewage discharge by numerical simulation analysis.ACTA Scientiarum Naturalium Universitatis Sunyatseni，2001，40(S2):126-129.

［15］Wu H，Wang Z L，Huang J，Zhao Y.Determination of pollutant discharge region's area on municipal wastewater release to sea project of Shenzhen.The Administration and Technique of Environmental Monitoring，2002，14(6):37-40.

［16］Yang S S.Simulation on the Flow Field，Suspended Sediment，Sediment Transport of the Sewage Marine Disposal Engineering in Yangkou Port Economic Development Zone of Nantong.Tianjin:Tianjin Research Institute for Water Transport Engineering，2010.

［17］Wang X H，Zhan SF，Peng ST，Yu H.Research on discharge position selection of sewage marine disposal in Yangkou Port.Journal of Waterway and Harbor，2011，32(2):140-143.

［18］Bai JF，Peng S T，Wang X H.Simulation and Analysis on the Diffuser Structure of the Standard Sewage Marine Disposal Engineering in Tianjin Nangang Industrial Zone.Tianjin:Tianjin Research Institute for Water Transport Engineering，2013.

［19］Duer-Michael J.Use of variable orifice “duckbill”valves for hydraulic and dilution optimization of multiport diffusers.Water Science and Technology，1998，38(10):277-284.

［20］Wang H W.Investigation of Buoyant Jet Discharges Using DPIV and PLIF ［D］. Singapore: Nanyang Technological University，2000.

［21］Tian X D，Philip J W.Marine wastewater discharges from multiport diffusers.Ⅰ:Unstratified stationary water.Journal of Hydraulic Engineering，2005，130(12):1137-1146.

［22］Tian X D，Philip J W.Marine wastewater discharges from multiport diffusers.Ⅱ:Unstratified flowing water.Journal of Hydraulic Engineering，2005，130(12):1147-1155.

［23］Daviero G J，Roberts P J.Marine wastewater discharges from multiport diffusers.III:Stratified stationary water.Journal of Hydraulic Engineering，2006，132(4):404-410.

［24］Lai A C H，Yu D，Lee J H W.Initial dilution of rosette buoyantjet group in crossflow //Proceedings of the 6th International Symposium on Stratified Flow. Perth，Australia，2006.

［25］Lai A C H，Yu D，Lee J H W.Near and intermediate field mixing of arosette jet group//Proceedings of the5th International Symposium on Environmental Hydraulics. Tempe，Arizona，2007.

［26］Patalano A，Corral M，Rodriguez A，Garcia M，Bleninger T.Submarine outfall design methodology for Argentinean Coasts//Proceedings of International Symposium on Outfall Systems.Argentinean，2011.

［27］Xu G T. Research on the Dilution and Diffusion Law of Wastewater Ocean Disposal with Multi Nozzle Jet［D］.Shanghai:Tongji University，1999.

［28］Xu G T，Wei H P.Physical and mathematical model study on dilution and dispersion in the near field of Jiaxing wastewater ocean disposal project.Marine Environmental Science，2001，20(2):42-46.

［29］Tang H W，Li X W.Experimental study on characteristics of turbulent jet for duck-bill valve diffuser.Advances in Water Science，2001，12(2):153-159.

［30］Zhou L W. Mathematics Simulation on Three-dimensional Unsteady Flow Field of Circle Floating Jet［D］.Dalian:Dalian University of Technology，2004.

［31］Fang S G.Research on the Behavior of Obstacled Buoyant Jets［D］.Wanhan:Wuhan University，2005.

［32］Gao Z.Study on Mixing Characteristics of Elliptical Jet and Its Application［D］.Nanjing:Hohai University，2005.

［33］Xiao Y.The Study on the Mixing Characteristics of A Tandem Multiple Jet Group in Cross-Flow ［D］. Nanjing: Hohai University，2005.

［34］Xiao Y，Tang H W，AYIDEN.Experiments on the scalar concentration field of a multiple jet group in cross-flow.Journal of Experiments in Fluid Mechanics，2011，25(5):35-39.

［35］ Bai JF，Yu H，Wang X H，Zhao H X，Chen Y H L，Jing L，He Z H，Zhang C Y.Diffuser with Multi-nozzles and Big Jet Angle Used on Sewage Marine Disposal:China，201220492554.4，2013-4-17.

［36］Bleninger T，Avanzini C A，Jirka G H.Hydraulic and technical evaluation of single diameter diffusers with flow rate control through calibrated，replaceable port exits//Proceedinds of the 3rd International Conference on Marine Waste Water Discharges and Marine Environment.Catania，Italy，2004.

［37］Law A K W，Lee C C，Qu Y.CFD modeling of a multi-port diffuser in an oblique current//Proceedings of the 2nd International Conference on Marine Waste Water Discharges.Istanbul，Turkey，2002:1-8.

［38］Bleninger T，Lipari G，Jirka G H.Design and optimization program for internal diffuser hydraulics//Proceedings of the 2nd International Conference on Marine Waste Water Discharges.Istanbul，Turkey，2002:9-15.

［39］Shannon N R，Mackinnon PA，Hamill G A.Evaluation of a CFD model of saline intrusion in marine outfalls//Proceedings of the 2nd International Conference on Marine Waste Water Discharges.Istanbul，2002:16-20.

［40］Lai Y F.Experimental study of hydraulic characteristics on a sewage discharge diffuser. Shanghai Environmental Science，1995，14(5):11-15.

［41］Wang C J，Wang K Q.Hydraulic design and model test of a diffuser for wastewater discharge.Marine Environmental Science，2001，20(2):47-50.

［42］Wu W.Study on Characteristics of Saline Intrusion and Purging in Marine Sewage Outfalls［D］.Nanjing:Hohai University，2006.

［43］Zeng X Q，Guo ZR，Huang Z F.Study on the intrusion and purge of seawater of Shenzhen ocean outfall under the undesigned condition.Marine Environmental Science，1997，16(3):29-33.

［44］Liu C.Sediment Research on Sewage Marine Disposal［D］.Shanghai:Tongji University，1999.

［45］Qian D R，Chen Z J，Wei H P.Study on hydraulic characteristics of the ocean diffusers with equal section.Sichuan Environment，2001，20(1):9-12.

［46］Zhao Y S，Liu WL，F(xiàn)eng Q.The Impact of the waves on the flow of the sewage diffuser.Chinese Quarterly of Mechanics，2009，30(1):154-163.

參考文獻(xiàn):

［1］全國(guó)環(huán)境統(tǒng)計(jì)公報(bào).北京:中華人民共和國(guó)環(huán)境保護(hù)部，2010.

［3］張永良，閻鴻邦.污水海洋處置技術(shù)指南.北京:中國(guó)環(huán)境科學(xué)出版社，1995.

［4］夏青.污水海洋處置工程設(shè)計(jì)理論與方法.北京:中國(guó)環(huán)境科學(xué)出版社，1996.

［10］孫長(zhǎng)青，趙可勝，孫英蘭.八所近岸海域排污口選址及水質(zhì)預(yù)測(cè).青島海洋大學(xué)學(xué)報(bào)，1994，(S1):209-216.

［11］尹毅，仲維妮，常乃環(huán)，李偉，許廣森.秦皇島市北戴河區(qū)污水排海方案的研究.海岸工程，1997，16(2):6-12.

［12］趙俊杰，白靜，康蘇海.南通污水排海工程排放點(diǎn)比選的數(shù)值模擬.海岸工程，2012，31(2):31-38.

［13］侯堋.小洋口近岸污水排海及其影響三維數(shù)學(xué)模型［D］.南京:河海大學(xué)，2008.

［14］宋強(qiáng)，吳航.深圳市政污水排海工程對(duì)海域環(huán)境影響的數(shù)值模擬研究.中山大學(xué)學(xué)報(bào):自然科學(xué)版，2001，40(S2):126-129.

［15］吳航，王澤良，黃劍，趙云.深圳市政污水排海工程排污混合區(qū)范圍的確定.環(huán)境監(jiān)測(cè)管理與技術(shù)，2002，14(6):37-40.

［16］楊樹森.南通市洋口港經(jīng)濟(jì)開發(fā)區(qū)一期污水處理排海工程潮流場(chǎng)、懸沙場(chǎng)數(shù)模及泥沙沖淤分析計(jì)算研究.天津:交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，2010.

［17］王心海，詹水芬，彭士濤，于航.洋口港污水排海工程排污口選劃研究.水道港口，2011，32(2):140-143.

［18］白景峰，彭士濤，王心海.天津南港工業(yè)區(qū)達(dá)標(biāo)尾水深海排放工程擴(kuò)散器模擬分析及結(jié)構(gòu)型式研究.天津:交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，2013.

［27］徐高田.污水海洋處置多噴口浮射流稀釋擴(kuò)散規(guī)律研究［D］.上海:同濟(jì)大學(xué)，1999.

［28］徐高田，韋鶴平.嘉興市污水海洋處置工程近區(qū)稀釋擴(kuò)散物理模型和數(shù)學(xué)模型研究.海洋環(huán)境科學(xué)，2001，20(2):42-46.

［29］唐洪武，李行偉.鴨嘴閥擴(kuò)散器紊動(dòng)射流特性試驗(yàn)研究.水科學(xué)進(jìn)展，2001，12(2):153-159.

［30］周連偉.三維非恒定流場(chǎng)中圓形浮射流的數(shù)值模擬研究［D］.大連:大連理工大學(xué)，2004.

［31］方神光.有障礙浮力射流特性的研究［D］.武漢:武漢大學(xué)，2005.

［32］高柱.橢圓射流摻混特性研究及應(yīng)用［D］.南京:河海大學(xué)，2005.

［33］肖洋.橫向流動(dòng)條件下多孔水平動(dòng)量射流摻混特性研究［D］.南京:河海大學(xué)，2005.

［34］肖洋，唐洪武，阿衣丁別克·居馬拜.橫流中多孔射流的稀釋特性實(shí)驗(yàn)研究.實(shí)驗(yàn)流體力學(xué)，2011，25(5):35-39.

［35］白景峰，于航，王心海，趙宏鑫，陳瑤泓伶，井亮，何澤慧，張春意.用于污水深海排放的多噴口大射流角度擴(kuò)散器:中國(guó)，201220492554.4，2013-4-17.

［40］賴翼峰.某排污擴(kuò)散器水力特性的實(shí)驗(yàn)研究.上海環(huán)境科學(xué)，1995，14(5):11-15.

［41］王春節(jié)，王可欽.污水排海工程的擴(kuò)散器水力設(shè)計(jì)與模型試驗(yàn)研究.海洋環(huán)境科學(xué)，2001，20(2):47-50.

［42］吳瑋.污水海洋擴(kuò)散器海水入侵及清除特性研究［D］.南京:河海大學(xué)，2006.

［43］曾小清，郭振仁，黃章富.深圳媽灣排海工程非設(shè)計(jì)條件下海水入侵與沖洗研究.海洋環(huán)境科學(xué)，1997，16(3):29-33.

［44］劉成.污水海洋處置工程泥沙問題研究［D］.上海:同濟(jì)大學(xué)，1999.

［45］錢達(dá)仁，陳祖軍，韋鶴平.污水海洋處置等截面擴(kuò)散器水力特性的研究.四川環(huán)境，2001，20(1):9-12.

［46］趙毅山，劉維祿，馮奇.波浪對(duì)污水?dāng)U散器出流量影響.力學(xué)季刊，2009，30(1):154-163.