朱 好，鄭平輝，王曉亮，鄭 偉，白曉平 (中國核電工程有限公司，北京 100840)

現(xiàn)階段我國有一些核設(shè)施采用天然蒸發(fā)池來處理低放射性水平的液態(tài)流出物。通常經(jīng)凈化處理后的液態(tài)流出物采取向環(huán)境水體排放的方式，但對于受納水體稀釋擴散能力較小或者無受納水體的內(nèi)陸廠址，液態(tài)流出物可采用大氣擴散稀釋方式排放，采用天然蒸發(fā)池進行處理是一種可選的技術(shù)方案，該方案不僅在技術(shù)及經(jīng)濟上有一定優(yōu)勢，不增加能耗，還可以盡量減少二次廢物的產(chǎn)生［1］。

在評價天然蒸發(fā)池對周圍環(huán)境及公眾的影響時，需要計算蒸發(fā)池周圍80 km半徑范圍內(nèi)的年均大氣彌散因子(即相對濃度)，通常將蒸發(fā)池作為地面面源考慮，采用基于一維高斯煙羽擴散模式的虛點源法［2］。該方法對于某些低風(fēng)速頻率較高的內(nèi)陸廠址，應(yīng)用中有一些局限性，如煙羽模式中濃度與風(fēng)速成反比，在靜風(fēng)和接近靜風(fēng)條件下會導(dǎo)致計算濃度的不合理增大;煙羽模式是靜態(tài)模式，不考慮擴散的時間過程，無法反應(yīng)低風(fēng)速條件下污染物的累積作用;煙羽模式不考慮氣象場的空間變化，無法很好地反映地形的作用。此外，上述方法對于有風(fēng)條件下靠近面源的受體也是一種較差的近似方法。美國環(huán)保署推薦的CALPUFF模式采用三維拉格朗日高斯煙團結(jié)合數(shù)值積分的方法來處理非均一、非平穩(wěn)氣象條件下面源的擴散，克服了上述方法模擬的局限性［3］，且CALPUFF模式在低風(fēng)速頻率較高的內(nèi)陸廠址通過了現(xiàn)場示蹤物試驗的定量校驗，驗證了小靜風(fēng)下模式模擬結(jié)果的可靠性［4－5］。

美國核管理委員會的技術(shù)報告［6］指出，將模式模擬值與實測值比較是進行模擬結(jié)果驗證的最好方法，但無實測值時，可將簡單模式與高級模式的模擬結(jié)果對比，以判斷簡單模式的適宜性。筆者以某低風(fēng)速頻率較高的內(nèi)陸廠址為例，針對虛擬的天然蒸發(fā)池，分別采用一維虛點源法和三維拉格朗日高斯煙團模式CALPUFF，模擬計算廠址80 km半徑范圍內(nèi)的年均大氣彌散因子，并對結(jié)果進行了對比分析，以期為位于類似廠址的蒸發(fā)池等地面面源的年均大氣彌散計算和后續(xù)環(huán)境影響評價提供參考。

1 計算模型和輸入數(shù)據(jù)

1．1 虛源法 虛源法將面源近似簡化為位于一定上風(fēng)向位置的虛擬點源，面源在下風(fēng)向造成的濃度分布可認(rèn)為是該虛點源的貢獻［7］。橫風(fēng)向源寬可近似看作虛點源擴散在源位置產(chǎn)生的煙羽寬度，因而可對虛點源采用高斯?jié)舛裙?。后退到源上風(fēng)向的虛點源離面源的距離X'對應(yīng)值σy的計算公式為:

式中，d為橫風(fēng)向面源寬度，m。

地面釋放源長期彌散因子計算公式:

對于蒸發(fā)池這類地面面源，采用廠址10 m高度的風(fēng)向、風(fēng)速、大氣穩(wěn)定度三維聯(lián)合頻率以及廠址實際擴散參數(shù)來計算年均大氣彌散因子。由圖1可知，廠址低風(fēng)速頻率較高，靜風(fēng)頻率為14．4%，0．5 ～2 m/s的小風(fēng)頻率高達 64．2%;最多風(fēng)向為 WSW，風(fēng)頻為19．8%，次多風(fēng)向為 SW，風(fēng)頻為12．1%。大氣穩(wěn)定度以中性的D類為主，占60%，不穩(wěn)定類(A、B、C)和穩(wěn)定類(E、F)分別占22%和18%。圖2為根據(jù)廠址風(fēng)洞試驗和數(shù)值模擬綜合得到的大氣擴散參數(shù)。受廠址復(fù)雜地形和建筑物擾動的影響，廠址大氣擴散參數(shù)較對應(yīng)穩(wěn)定度級別的P－G擴散參數(shù)明顯偏大。

1．2 CALPUFF模式 CALPUFF高斯煙團模式為美國國家環(huán)保署(USEPA)推薦的適用于長距離輸送和涉及復(fù)雜流動(如地形和土地利用變化導(dǎo)致的氣象場空間變化、小靜風(fēng)、熏煙、環(huán)流等)近場應(yīng)用的導(dǎo)則模式［3］。CALPUFF模式系統(tǒng)包括3個主要模塊:氣象模式CALMET可利用常規(guī)地面和探空氣象觀測資料，針對關(guān)心區(qū)域的地形、地表狀態(tài)等，利用質(zhì)量守衡方程對風(fēng)場進行診斷;輸送和擴散模式CALPUFF，利用CALMET輸出的風(fēng)場和微氣象學(xué)參數(shù)，通過跟蹤從排放源釋放的離散煙團的運行來模擬煙團移動路徑上的擴散、轉(zhuǎn)化和清除;后處理模塊CALPOST用于計算平均濃度和沉積通量，輸出選定區(qū)域和時段的平均濃度和最大濃度等。CALPUFF采用各向同性分布的高斯煙團和沿風(fēng)向拉伸的煙團兩種積分采樣方案來計算煙團的平流、輸送和擴散等。積分煙團模式適合中等尺度范圍，拉伸煙團模式適合局地尺度?？臻g某點的污染物濃度由源不斷釋放的煙團的擴散貢獻疊加得到。描述某一煙團在一個受體點的濃度貢獻的基本方程為:

式中，C為地面濃度;Q為一個煙團中污染物質(zhì)量;σx、σy、σz分別為縱向、橫向和垂直方向的擴散參數(shù)，da、dc分別為縱向和橫向煙團中心距受體的距離;g為高斯方程中的垂直項;He為煙團中心距離地面的有效高度;h為混合層高度。公式(4)中的求和項代表了混合層頂和地面之間的多次反射。

CALPUFF模式的主體框架具有自然適應(yīng)低風(fēng)速條件擴散過程的能力，通過對常規(guī)算法做若干調(diào)整，如改變拉伸煙團的釋放方式、煙流抬升方式、近場效應(yīng)的模擬方案和擴散參數(shù)的增長方式等，使其適應(yīng)小靜風(fēng)條件。CALPUFF采用二維數(shù)值積分算法來模擬面源下風(fēng)向的濃度分布，其將多邊形面源的各個頂點轉(zhuǎn)化為順風(fēng)向和橫風(fēng)向坐標(biāo)，采用兩個誤差函數(shù)的差值對橫風(fēng)向進行積分，對于順風(fēng)向，采用了更加有效和快速收斂的Romberg數(shù)值積分算法［3］。

針對研究區(qū)域，設(shè)置兩套模擬網(wǎng)格，均以虛擬蒸發(fā)池面源中心為投影中心，小尺度模擬區(qū)域為10 km×10 km，水平氣象網(wǎng)格距為100 m，采樣網(wǎng)格距為50 m，大尺度模擬區(qū)域為160 km×160 km，水平氣象網(wǎng)格距與采樣網(wǎng)格距均設(shè)為1 000 m。兩套模擬網(wǎng)格在垂直方向均劃分為10層，各層高度分別為20、40、80、160、300、600、1 000、1 500、2 200 和3 000 m。模式采用美國地質(zhì)勘探局(USGS)提供的約90 m分辨率的地形高程資料和900 m分辨率的全球土地覆蓋數(shù)據(jù)，氣象資料包括投影中心200 km范圍內(nèi)12個格點的NCEP地面氣象資料和每日4次的垂直廓線資料、廠址周邊13個常規(guī)地面站以及現(xiàn)場氣象鐵塔全年觀測資料。煙團擴散采用CALMET輸出的時空變化的微氣象學(xué)參數(shù)，根據(jù)相似性理論進行計算。

2 結(jié)果與分析

假定虛擬天然蒸發(fā)池的尺寸為300 m×200 m，簡化起見，假定蒸發(fā)池有效蒸發(fā)面積與實際面積相同。CALPUFF利用廠址區(qū)域全年逐時氣象資料結(jié)合面源排放情況，模擬得到全年逐時濃度場，將全年8 760 h的濃度場進行平均，并采用面積平均的方法將各扇區(qū)內(nèi)的網(wǎng)格濃度轉(zhuǎn)換為相應(yīng)的扇區(qū)平均值，以便與虛點源法的計算結(jié)果進行比較。

CALPUFF計算的廠址區(qū)域年均大氣彌散因子分布見圖3。高濃度區(qū)出現(xiàn)在地面面源周圍。受廠址區(qū)域較高的小靜風(fēng)頻率的影響，地面面源在各方位的大氣彌散因子差別并不顯著，主要由于低風(fēng)速條件下風(fēng)向多變，煙羽的水平慢擺作用導(dǎo)致其在水平方向側(cè)向擴散很寬，近源處甚至能出現(xiàn)上風(fēng)向擴散［8］，因而年均大氣彌散因子整體上并未像大風(fēng)或中等風(fēng)速條件表現(xiàn)出明顯的主導(dǎo)風(fēng)向下風(fēng)向的擴散優(yōu)勢。此外，低風(fēng)速條件下，風(fēng)場易受局部地形等因素的影響，空間變化大，因而地形也會影響煙羽的擴散形態(tài)。由圖3可知，在丘陵山區(qū)地形有明顯起伏的地區(qū)，年均大氣彌散因子等值線與地形等高線對應(yīng)較好，表明CALPUFF模式能夠較好地反映廠址區(qū)域復(fù)雜地形對大氣彌散的影響。

對比虛源法與CALPUFF計算的方位最大年均大氣彌散因子隨下風(fēng)距離的變化(圖4)，可見，10 km以內(nèi)，CALPUFF計算結(jié)果較虛源法更為保守。在0．5～1 km，CALPUFF計算結(jié)果是虛源法的2倍左右，隨著距離增加，二者差異逐漸減小;10～30 km，兩種方法計算結(jié)果相當(dāng);30～80 km，虛源法計算結(jié)果大于CALPUFF計算值，隨著距離增加，二者差異增大;30～40 km，虛源法計算結(jié)果較CALPUFF偏大0．5倍;而70～80 km，兩種方法計算結(jié)果相差達3倍左右。此外，虛源法基于高斯直線煙流模式，采用廠址單點氣象觀測資料來代表整個評價區(qū)域的氣象條件，廠址80 km半徑范圍內(nèi)，其計算的方位最大年均大氣彌散因子均出現(xiàn)在主導(dǎo)風(fēng)向下風(fēng)向(即ENE方位)?！逗税踩珜?dǎo)則》HAD 101/02《核電廠廠址選擇的大氣彌散問題》［9］指出，高斯煙羽模式的適用條件是氣象參數(shù)在時間和空間上大致相同，適用范圍大約為10 km，超出該范圍，其適用性就降低。隨著距離的增加，模式中通常采用的均勻性假設(shè)愈來愈不適用。因而，對于廠址區(qū)域地形復(fù)雜、低風(fēng)速頻率較高、風(fēng)場空間變化較大的情形，采用虛源法進行80 km范圍的大氣擴散計算本質(zhì)上來說是不適宜的。CALPUFF這類三維拉格朗日軌跡類模式則可以滿足區(qū)域尺度的應(yīng)用需求，反映煙羽擴散的時空變化特征。上述對比表明，對于環(huán)境影響評價中關(guān)注的方位最大彌散因子及其出現(xiàn)位置，基于高斯煙羽模式的虛源法由于不考慮低風(fēng)速條件下污染物在近區(qū)的累積效應(yīng)以及氣象條件的空間變化，因而其在小靜風(fēng)頻率較高廠址的計算結(jié)果存在不保守情形，對于下墊面非平坦均一的情形，其預(yù)測的年均最大彌散因子出現(xiàn)位置與實際情況有所偏離。

虛源法與CALPUFF計算的各對應(yīng)扇區(qū)年均大氣彌散因子之比隨下風(fēng)距離和方位的變化分布見圖5。從中值分布來看，在45 km以內(nèi)，虛源法計算結(jié)果整體上小于CALPUFF計算值，0．5～1 km扇區(qū)，中值偏低達1個量級左右，55～75 km，二者計算結(jié)果相當(dāng);在NE和ENE方位(即主導(dǎo)風(fēng)向下風(fēng)向)，兩種方法計算值相當(dāng)，其他方位CALPUFF計算值更偏保守。主要由于CALPUFF煙團模式能更好地模擬小靜風(fēng)條件下多方位的彌散，而虛源法這類高斯煙羽模式通常將污染物按照集中在主導(dǎo)風(fēng)下風(fēng)向輸送來處理，因而在其他方位會導(dǎo)致計算結(jié)果的偏低。整體來看，兩種方法計算值之比在各距離和各方位的離散達1個量級左右。

3 結(jié)論

目前，在核設(shè)施天然蒸發(fā)池等地面面源的環(huán)境影響評價中多采用基于高斯煙羽模式的虛源法。對于某些低風(fēng)速頻率較高的內(nèi)陸廠址，該方法在應(yīng)用中有一些局限性。以我國某一低風(fēng)速頻率較高的內(nèi)陸廠址為例，分別采用基于一維高斯煙羽模式的虛源法和三維拉格朗日高斯煙團模式CALPUFF結(jié)合數(shù)值積分法，計算了某虛擬蒸發(fā)池的地面年均大氣彌散因子。結(jié)果表明，虛源法由于不考慮小靜風(fēng)條件下污染物的近區(qū)累積效應(yīng)和多方位彌散，因而在低風(fēng)速頻率較高的廠址，對方位最大彌散因子和各扇區(qū)年均大氣彌散因子的計算結(jié)果存在不保守情形。對于復(fù)雜下墊面，虛源法預(yù)測的年均最大彌散因子出現(xiàn)位置與實際情況有所偏離。對于這類廠址，在評價蒸發(fā)池等面源的環(huán)境影響時，建議選擇其他適宜的大氣擴散模式和方法。

［1］王文靜，李杰，廉冰，等．核設(shè)施蒸發(fā)池對周圍環(huán)境的影響研究［C］//中國核學(xué)會輻射防護分會．輻射防護分會2012年學(xué)術(shù)年會論文集．中國核學(xué)會輻射防護分會，2012:196－199．

［2］劉克強，鄒長貴．天然蒸發(fā)池核素大氣遷移研究(英文)［J］．中國核科技報告，1997(00):991－1003．

［3］SCIRE J S，STRIMAITIS D G，YAMARTINO R J．A user’s guide for the CALPUFF dispersion model(version 5)［R］．Concord，MA:Earth Tech Inc，2000．

［4］朱好，蔡旭暉，張宏升，等．內(nèi)陸丘陵河谷地區(qū)小風(fēng)條件下的大氣擴散模擬研究［J］．環(huán)境科學(xué)學(xué)報，2011，31(3):613 －623．

［5］朱好，張宏升，蔡旭暉，等．CALPUFF在復(fù)雜地形條件下的近場大氣擴散模擬研究［J］．北京大學(xué)學(xué)報:自然科學(xué)版，2013，49(3):452 －462．

［6］MOLENKAMP C R，BIXLER N E，MORROW C W，et al．Comparison of average transport and dispersion among a Gaussian，a two-dimensional，and a three-dimensional model，NUREG/CR － 6853 ［R］．Washington，DC，2004．

［7］蔣維楣．空氣污染氣象學(xué)［M］．南京:南京大學(xué)出版社，2003:95－96．

［8］ARYA S P．Modeling and parameterization of near-source diffusion in weak winds［J］．Journal of Applied Meteorology，1995，34(5):1112 －1122．

［9］國家核安全局．HAD101－02核電廠廠址選擇的大氣彌散問題［S］．［出版地不詳］:［出版社不詳］，1987．