1、大氣穩(wěn)定度分級(jí)常用的大氣穩(wěn)定度分類方法有帕斯奎爾(Pasquill)法和國(guó)標(biāo)原子能機(jī)構(gòu) IAEA 推薦的方法。這里介紹的是中國(guó)現(xiàn)有法規(guī)中推薦的修訂帕斯奎爾分類法（簡(jiǎn)記 PS） ，分為強(qiáng)不穩(wěn)定、不穩(wěn)定、弱不穩(wěn)定、中性、較穩(wěn)定和穩(wěn)定六級(jí)。它們分別表示為A、B、C、D、E、F。確定等級(jí)時(shí)首先計(jì)算出太陽(yáng)高度角按表 B1查出太陽(yáng)輻射等級(jí)數(shù)，再由太陽(yáng)輻射等級(jí)數(shù)與地面風(fēng)速按表 B2查找穩(wěn)定等級(jí)。云量，1/10總云量/低云量夜間ho15 15ho35 35ho65 ho6544-2-11235-74-1012384-1001155-7000018800000表 B1 太陽(yáng)輻射等級(jí)數(shù)太 陽(yáng) 輻 射 等 級(jí) 數(shù)注

2、：云量（全天空十分制）觀測(cè)規(guī)則與現(xiàn)國(guó)家氣象局編定的地面氣象觀測(cè)規(guī)范相同。3210121.9AA-B BDEF2-2.9A-BBCDEF3-4.9BB-CCDDE5-5.9CC-DDDDD6DDDDDD地 面 風(fēng)速，m/s太 陽(yáng) 輻 射 等 級(jí)表 B2 大氣穩(wěn)定度的等級(jí)注：地面風(fēng)速（m/s）系指距地面 10m 高度處 10min 平均風(fēng)速，如使用氣象臺(tái)（站）資料，其觀測(cè)規(guī)則與國(guó)家氣象局編定的地面氣象觀測(cè)規(guī)范相同。太陽(yáng)高度角ho使用下式計(jì)算： .(30015coscoscossinsinarcsinthoB1)式中： ho-太陽(yáng)高度角，deg ；-當(dāng)?shù)鼐暥龋琩eg.；-當(dāng)?shù)亟?jīng)度；deg ；t-進(jìn)行

3、觀測(cè)時(shí)的北京時(shí)間；-太陽(yáng)傾角，deg，可按下式計(jì)算：/1803sin001480. 03cos002697. 02sin000907. 02cos006758. 0sin070257. 0cos39912. 0006918. 0oooooo式中： o-360dn/365，deg； dn- 一年中日期序數(shù)，0、1、2、364。1.1.1.1 地形對(duì)煙羽的影響此前的擴(kuò)散模式都假設(shè)地面是完全平整的（煙囪底部是一個(gè)無(wú)限大的水平面，其高程為 0） ，因此在擴(kuò)散過程中煙羽的中心線可保持水平不變。但如果在預(yù)測(cè)點(diǎn)（x,y,z）處，地面有一定的高程 hT(0z)，則在對(duì)（x,y,z）式應(yīng)用以上模式時(shí)，應(yīng)對(duì)有效煙

4、羽高度進(jìn)行一些修正。假定煙羽路徑始終與起伏的地形保持平行，或者假設(shè)煙羽軸線保持固定的海拔高度，并與高于煙羽的地形相交，都是不正確的，實(shí)際情況應(yīng)該是介于上述二者之間。具體的修正方法如下。（1）中性和不穩(wěn)定天氣條件令：hT為凸出的地形高度；He為煙軸高度（即有效高度） ；T為煙軸高度修正系數(shù)（或地形系數(shù)） ，修正后的煙囪有效高度應(yīng)該是 THe oT 則應(yīng)按下式取值：時(shí)當(dāng)時(shí)當(dāng)TeeTeTehHHhHThHT,/ )2/(, 2/1（2）穩(wěn)定天氣條件（D-E、E、F）在穩(wěn)定天氣條件下，當(dāng)煙羽逼近孤立山體時(shí)，煙羽以臨界高度 Hc為界分成兩部分，臨界高度以上的煙羽有足夠的動(dòng)能爬越山體，而臨界高度以下的煙羽

5、，只能被迫繞著山體過去。臨界高度 H c 可由下式確定：Hc=Hm-u/(gd/dz)1/2式中 Hm-孤立山體高度，m； Hc-臨界高度，m； -z 高度處大氣位溫，K； d/dz-z 高度處位溫梯度，K/m； u-平均風(fēng)速，m/s； g-重力加速度，m/s2。例如，=300K，d/dz=0.01K/m，u=2m/s，Hm=200m；則Hc=Hm-111=89m，煙囪有效高度大于 89m 時(shí)，煙羽將有足夠的動(dòng)能爬越山體。對(duì)于山體高度 Hm已定的情況，大氣越穩(wěn)定，則 Hc 越小。所以一般只需計(jì)算在 F 穩(wěn)定度下的 Hc，如果煙羽有效高度 HeHc(F)，則可認(rèn)為煙羽能夠爬越山體。1.2點(diǎn)源擴(kuò)散

6、模式1.2.1持續(xù)排放源1.2.1.1 有風(fēng)模式（U101.5m/s）1.自由空間中的連續(xù)點(diǎn)源實(shí)際上絕大多數(shù)污染源都是連續(xù)的，對(duì)于連續(xù)排放源，可理解為在時(shí)間上依次連續(xù)釋放無(wú)窮多個(gè)煙團(tuán)。因此，連續(xù)排放源的擴(kuò)散模式可以通過將瞬時(shí)單煙團(tuán)模式對(duì) to從到 t 積分后求得。以煙團(tuán)初始空間坐標(biāo)為原點(diǎn)，下風(fēng)方為 x 軸，煙羽軸線與 x 軸一直保持重合，都是 x 的函數(shù)，將對(duì) to的積分變zyx,換為對(duì)（xuT）/x的積分，可得最基本的煙羽擴(kuò)散模式：22222exp2exp2,zyzyzyuQzyxC適用條件為：自由空間；風(fēng)速要比較大（u101.5m/s） ；當(dāng)大氣不穩(wěn)定狀態(tài)時(shí)，可能帶來一定的誤差。2.地面反

7、射用像源法，假想地平線為一鏡面，在其下方有一與真實(shí)源完全對(duì)稱的虛源，則這兩個(gè)源疊加后的效果和真實(shí)源考慮到地面反射的結(jié)果是等價(jià)的。以煙囪地面位置的中心點(diǎn)為坐標(biāo)原點(diǎn)，在考慮到地面反射后，污染源下風(fēng)方任一點(diǎn)小于 24 小時(shí)取樣時(shí)間的污染物濃度 C（x，y，z）由下式給出：2222222exp2exp2exp2,zezeyzyHzHzyuQzyxCz=0 時(shí)的地面濃度 C（x,y,0） ，可簡(jiǎn)化為；222222exp0 ,zeyzyHyuQyxC下風(fēng)方 X 軸線上（y=0）的地面濃度 C（x,0,0）為：22212exp0 , 0 ,zeHuQxC對(duì)于較低的排放源（例如 He50m，具體限值由地面粗糙

8、度、混合層高度等因素決定） ，一般可直接應(yīng)以上式子計(jì)算。3.混合層反射對(duì)于高架源，當(dāng)超過一定的下風(fēng)距離時(shí)，需對(duì)煙羽在混合層頂?shù)姆瓷溥M(jìn)行修正。同考慮地面反射類似，用像源法修正后，污染源下風(fēng)方任一點(diǎn)小于 24 小時(shí)取樣時(shí)間的污染物濃度C（x,y,z）可表示為：FyuQzyxcyzy222exp2),(kknzezezHnhzHnhF222222exp22exp式中 h混合層高度； k反射次數(shù)，一、二級(jí)項(xiàng)目取 k=4 已足夠。4.側(cè)面反射詳見狹長(zhǎng)山谷擴(kuò)散模式。1.2.1.2 小風(fēng)靜風(fēng)模式（U101.5m/s）小風(fēng)靜風(fēng)時(shí)，污染物地面濃度 C(x,y,0)可用下式計(jì)算： GQYXcL2022/322,式

9、中 和 G 按下式計(jì)算：，2202201222eHYX seseGU2/32/2201221， dtestsoo2/22101UXS 和分別是橫向和鉛直向擴(kuò)散參數(shù)的回歸系數(shù)(y=z=0102T,z=T),T為擴(kuò)散時(shí)間(s),和的定值見 HJT2.2-199301020102附錄 B3。注意，上式中 He是煙筒有效高度相對(duì)于預(yù)測(cè)點(diǎn)的高度。若預(yù)測(cè)點(diǎn)高度坐標(biāo)為 z，煙筒有效高度處坐標(biāo)為 z0，則 He=z0-z。1.2.2非正常排放源非正常排放是指建設(shè)項(xiàng)目生產(chǎn)運(yùn)行階段的開車、停車、檢修、一般性事故和發(fā)生漏泄等情況時(shí)的污染物的不正常排放。非正常排放常發(fā)生在有限時(shí)間（T）內(nèi)。以瞬時(shí)單煙團(tuán)正態(tài)擴(kuò)散式，對(duì)

10、t0在有限時(shí)間 T 內(nèi)積分，經(jīng)整理后可得非正常排放模式。1.2.2.1 有風(fēng)情況（U101.5m/s）非正常排放條件下的地面濃度ca(mg/m3)建議按下列各式計(jì)算。以排氣筒地面位置為原點(diǎn),有效源高為He，平均風(fēng)向軸為X軸，源強(qiáng)為Q（mg/s） ，開始非正常排放時(shí)的時(shí)間為 t，非正常排放持續(xù)時(shí)間為T，預(yù)測(cè)時(shí)刻的時(shí)間為t。1.有風(fēng)情況（U101.5m/s）t時(shí)刻任一點(diǎn)（x,y,z）的濃度，以持續(xù)排放源模式為基礎(chǔ)，乘上一個(gè)系數(shù) G1，按下式計(jì)算：1222exp2),(GFyuQzyxcyzykknzezezHnhzHnhF222222exp22exp11xxxxxxUtxUTUtxUtG式中 F混

11、合層反射項(xiàng)； G1非正常排放項(xiàng)； h混合層高度； k反射次數(shù)，一、二級(jí)項(xiàng)目取 k=4 已足夠。擴(kuò)散參數(shù)各指數(shù)、系數(shù)的定值見導(dǎo)則,2121XXzyx附錄 B。1.2.2.2 小風(fēng)靜風(fēng)（U10=1.5m/s）的高斯持續(xù)排放點(diǎn)源，并且要求穩(wěn)定度較不穩(wěn)定、混合層反射可忽略等條件，其計(jì)算結(jié)果 Xm必須在擴(kuò)散參數(shù)系數(shù) y1、y2和指數(shù) a1、a2的應(yīng)用范圍之內(nèi)。例如，以 1000m 范圍內(nèi)的擴(kuò)散參數(shù)計(jì)算系數(shù)和指數(shù)計(jì)算得Xm=3000m，則結(jié)果是不可靠的。由于解析式應(yīng)用范圍有限，應(yīng)用條件較苛刻，實(shí)際計(jì)算中常用數(shù)值法。借助計(jì)算機(jī)的快速計(jì)算，可使用分段逼近法求最大落地濃度。該原理是：將最大濃度可能的出現(xiàn)范圍，如

12、下風(fēng)向0，1000000m，分成 10 段共 11 個(gè)點(diǎn)位，每點(diǎn)求一濃度，再以最大濃度點(diǎn)位左、右兩點(diǎn)位作為新的計(jì)算范圍起、止點(diǎn)，再分10 段，如此循環(huán)計(jì)算，直到每段長(zhǎng)度小于要求的精度（如0.01m） 。一般計(jì)算幾十次即可得到結(jié)果。這種方法算法很簡(jiǎn)單，無(wú)須求函數(shù)導(dǎo)數(shù)，同時(shí)靈活性極強(qiáng)，因?yàn)閷?duì)每一個(gè)點(diǎn)求濃度的時(shí)候，可以按調(diào)用任意模式按任意方法計(jì)算。因此，這一方法可用于有風(fēng)、小風(fēng)靜風(fēng)、面源體源的點(diǎn)源修正法和非正常排放等所有單源模式的求最大濃度值。1.2.4對(duì)源強(qiáng)和有效源高的修正1.2.4.1 地形對(duì)煙羽的影響此前的擴(kuò)散模式都假設(shè)地面是完全平整的（煙囪底部是一個(gè)無(wú)限大的水平面，其高程為 0） ，因此在擴(kuò)

13、散過程中煙羽的中心線可保持水平不變。但如果在預(yù)測(cè)點(diǎn)（x,y,z）處，地面有一定的高程 hT(0z)，則在對(duì)（x,y,z）式應(yīng)用以上模式時(shí)，應(yīng)對(duì)有效煙羽高度進(jìn)行一些修正。假定煙羽路徑始終與起伏的地形保持平行，或者假設(shè)煙羽軸線保持固定的海拔高度，并與高于煙羽的地形相交，都是不正確的，實(shí)際情況應(yīng)該是介于上述二者之間。具體的修正方法如下。（1）中性和不穩(wěn)定天氣條件令：hT為凸出的地形高度；He為煙軸高度（即有效高度） ；T為煙軸高度修正系數(shù)（或地形系數(shù)） ，修正后的煙囪有效高度應(yīng)該是 THe oT 則應(yīng)按下式取值：時(shí)當(dāng)時(shí)當(dāng)TeeTeTehHHhHThHT,/ )2/(, 2/1（2）穩(wěn)定天氣條件（D-

14、E、E、F）在穩(wěn)定天氣條件下，當(dāng)煙羽逼近孤立山體時(shí)，煙羽以臨界高度 Hc為界分成兩部分，臨界高度以上的煙羽有足夠的動(dòng)能爬越山體，而臨界高度以下的煙羽，只能被迫繞著山體過去。臨界高度 H c 可由下式確定：Hc=Hmu/(gd/dz)1/2式中 Hm孤立山體高度，m； Hc臨界高度，m； z 高度處大氣位溫，K； d/dzz 高度處位溫梯度，K/m； u平均風(fēng)速，m/s； g重力加速度，m/s2。例如，=300K，d/dz=0.01K/m，u=2m/s，Hm=200m；則Hc=Hm-111=89m，煙囪有效高度大于 89m 時(shí)，煙羽將有足夠的動(dòng)能爬越山體。對(duì)于山體高度 Hm已定的情況，大氣越穩(wěn)定

15、，則 Hc 越小。所以一般只需計(jì)算在 F 穩(wěn)定度下的 Hc，如果煙羽有效高度HeHc(F)，則可認(rèn)為煙羽能夠爬越山體。1.2.4.2 熱浮力對(duì)煙羽的修正1.混合層頂穿透問題高架源熱浮力煙羽對(duì)混合層并非只有“完全穿透”和“完全不穿透”兩種情況，還須考慮到“部分穿透”問題。定義 P 為穿透系數(shù)：P=1.5-(h-H)/H當(dāng) P0 時(shí)，H1=minH,2(h-H)/3,Q1=Q；當(dāng) P 1 時(shí)，H1=h-H,Q1=0；當(dāng) 0P1 時(shí)，H1=(2/3+P/3)(h-H),Q1=Q(1-P)式中，Q1，H1分別為修正后的源強(qiáng)和抬升高度； h混合層高度。2.風(fēng)向切變和熱浮力問題高浮力煙羽在 z 方向受風(fēng)向

16、切變及熱浮力的影響可按下式對(duì)y和 z進(jìn)行修正：2yc=2y+(H/3.5)2+(0.032x2)2zc=2z+(H/3.5)2式中，yc、zc分別為修正后 y、z 方向的擴(kuò)散參數(shù)； (H/3.5)2熱浮力修正項(xiàng)； (0.032x2)風(fēng)向切變修正項(xiàng)。 按下表確定：表 風(fēng)向切變修正參數(shù) /Z（度/m）穩(wěn)定度ABCDEF/Z0.0050.0100.0150.0200.0250.0351.2.4.3 干沉積修正1.干沉積的含義在重力、湍流擴(kuò)散、分子擴(kuò)散、靜電引力以及其他生物學(xué)、化學(xué)和物理學(xué)等因素的作用下，大氣中的顆粒物或某些氣體隨時(shí)會(huì)被地表（土壤、植物、水體）滯留或吸收，使這些物質(zhì)連續(xù)不斷地從大氣向地

17、表作質(zhì)量轉(zhuǎn)移，從而減少其在空氣中的濃度。通常，把這一與降水作用無(wú)關(guān)的質(zhì)量轉(zhuǎn)移過程，稱之為干沉積。在對(duì)中距離大氣污染物輸送的研究中發(fā)現(xiàn)，差不多一半以上的質(zhì)量轉(zhuǎn)移是干沉積過程引起的。為些，對(duì)一些物理、化學(xué)特性和空氣相差較大的大氣污染物（例如顆粒物） ，或者對(duì)于輸送距離比較遠(yuǎn)的高架源的濃度預(yù)測(cè)，都應(yīng)當(dāng)考慮其干沉積的影響。下面將介紹兩種常用的干沉積擴(kuò)散模式。2.源虧損模式源虧損模式主要用于粒徑小于 10m 易產(chǎn)生沉積的顆粒物或氣體。假定因各種機(jī)理造成的大氣污染物向地面的沉積通量Wmg/(sm2)由下式表示：W=VdC.(4.9-79)式中 Vd沉積速度，m/s； C大氣污染物地面濃度，常取自距地面 1

18、m 高度處。大氣污染物自煙囪出口排出后，其初始源強(qiáng) Q（0）因沉積作用將隨下風(fēng)距離逐漸減弱（虧損） 。根據(jù)式（4.9-79）可導(dǎo)出虧損后的源強(qiáng) Q（x） ；因?yàn)閷?W 從到對(duì) y 積分恰是dQ(x)/dx，再對(duì) x 積分則可得 dxHuVQxQxzezd02212exp2exp0由式（4.9-80）可見，問題的關(guān)鍵是給出沉積速度 Vd。長(zhǎng)期以來，許多研究者為求得 Vd，在實(shí)驗(yàn)或模式化方面都進(jìn)行了大量的工作。對(duì)于 SO2，Hanna 曾推薦如下 Vd(cm/s)的野外實(shí)驗(yàn)值：高約 1m 的莊稼地、城市、水面為 0.7；高 0.1m 的草地為 0.5；酸性土壤為 0.4（干）或 0.6（濕） ；石

19、灰質(zhì)土壤為 0.8；干積雪為 0.1。在模式化方面，通常是把沉積速度和電流類比，假定沉積速度與各種阻力的總和成反比。上式可用數(shù)值積分計(jì)算。采用 1 小時(shí)取樣時(shí)間的國(guó)標(biāo)擴(kuò)散參數(shù)計(jì)算結(jié)果表明，D 穩(wěn)定度時(shí)虧損最明顯，不穩(wěn)定或更穩(wěn)定時(shí)虧損依次減少。設(shè) He=100m,Q（0）為 100，Vd為 0.007m/s，U為 3m/s，則在 10000m 處 Q 為下表所示：穩(wěn)定度ABB-CCC-DDD-EEFQ(10km)99.0896.2795.1793.9392.2191.9092.8194.7698.86一般來說，3000m 以內(nèi)的擴(kuò)散源虧損小于 1%，可以忽略不計(jì)。3.部分反射（斜煙羽）部分反射模

20、式主要用于粒徑大于 10m 的氣載顆粒物。因?yàn)檫@種污染物有明顯的沉降作用，其煙羽重心是是逐漸降低的，此外，地面只能反射一部分的污染物。對(duì)某一粒徑的顆粒污染物，對(duì) He 和 Q 進(jìn)行以下調(diào)整后，再按相應(yīng)氣態(tài)物模式計(jì)算：He(x)=He(0)-Vg x/u。Q(x)=Q(0) (1+)/2式中 顆粒物的地面反射系數(shù)（參閱下表） ； Vg顆粒物沉降速度，用 STOCKS 公式計(jì)算：；182gdVg d、分別為顆粒物的直徑和密度； g重力加速度； 空氣粘性系數(shù)。表 地面反射系數(shù) 粒徑范圍/m153031476076100平均粒徑/m22385085反射系數(shù) 0.80.50.301.2.4.4 濕沉積修

21、正濕沉積系指大氣污染物因降水而減少其在空氣中濃度的過程。嚴(yán)格地講，濕沉積分云中和云下兩種清除機(jī)制。在工程應(yīng)用中，常把這兩種機(jī)制合起來考慮。假設(shè)大氣污染物的初始源強(qiáng)Q（0）因降水隨下風(fēng)距離 x 成指數(shù)衰減，則Q（x）=Q(0)exp(x/u)式中 x接受點(diǎn)的下風(fēng)距離； u煙囪出口處的風(fēng)速；清除系數(shù)，S1。令 J 為降水強(qiáng)度（mm/h） ，根據(jù)實(shí)驗(yàn)定測(cè)量結(jié)果，對(duì)于SO2=1.7104J0.61.2.4.5 化學(xué)遷移的修正化學(xué)遷移的修正，在工程應(yīng)用上與濕沉積類似，可令修正后的源強(qiáng) Q（x）等于初始源強(qiáng) Q（0）乘以修正因子fc 。ccuTxf/exp 式中 Tc大氣污染物的時(shí)間常數(shù)；其他符號(hào)同前。如

22、果定義fc=1/2 所對(duì)應(yīng)的時(shí)間為該大氣污染物的半衰期td，將其代入式（44.9-85） ，可得：dcutxf/693. 0exp 對(duì)于城市尺度，SO2的 td4 小時(shí)；長(zhǎng)距離輸送 SO2的 td典型值為 4 天。1.3面源(體源)的點(diǎn)源修正算法把面源的排放當(dāng)作一個(gè)位于其幾何中心的點(diǎn)源的排放，對(duì)擴(kuò)散參數(shù)適當(dāng)修正后，采用點(diǎn)源模式直接計(jì)算，用以近似代表該面源的擴(kuò)散。但這種方法僅適用于有風(fēng)條件（U101.5m/s） 。在小風(fēng)靜風(fēng)的條件下，面源擴(kuò)散如何計(jì)算，尚未有定論，這里暫按小風(fēng)靜風(fēng)的點(diǎn)源方法計(jì)算。如果測(cè)點(diǎn)位于面源（體源）之外，這種近似算法較為準(zhǔn)確，但如果測(cè)點(diǎn)在面源內(nèi)部，則計(jì)算結(jié)果就很不理想。比如，

23、在有風(fēng)時(shí)，位于面源幾何中心上風(fēng)向各點(diǎn)（仍處于面源內(nèi)部） ，按這種近似算法濃度應(yīng)是 0，實(shí)際不然。此外，這種算法不能反映面源的形狀。這些都是點(diǎn)源修正算法的不足之外。有風(fēng)條件下，對(duì)于面源或體源內(nèi)部的預(yù)測(cè)點(diǎn)，可采用下式計(jì)算：,2jsQC 式中，x 為測(cè)點(diǎn)離面源上風(fēng)邊界的距離。這種方axH2222/法使得在有風(fēng)時(shí)，面源內(nèi)任意一點(diǎn)均有濃度。面源。采用直接修正法。如果面源的面積較?。⊿1km2） ，面源外的 Cs可按點(diǎn)源擴(kuò)散模式計(jì)算，只是應(yīng)附加一個(gè)初始擾動(dòng)。這一初始擾動(dòng)使煙羽在 x=0 就有一個(gè)和面源橫向?qū)挾认嗟鹊臋M向尺度，以及和面源高度相等的垂直向尺度。注意到煙羽的半寬度等于 2.15y或 2.15z，

24、則修正后的 y和 z分別為：15. 2/3 . 4/2211Hxaxazyay式中 x自接受點(diǎn)至面源中心點(diǎn)的距離； ay面源在 y 方向的長(zhǎng)度； H面源的平均排放高度。體源當(dāng)無(wú)組織排放源為體源時(shí)，地面濃度可按點(diǎn)源直接修正法計(jì)算；令 y、z分別表示體源在 y 和 z 方向的邊長(zhǎng)，則修正后的 y和 z分別為：3 . 4/3 . 4/2211zzyyxxx自接受點(diǎn)至面源中心點(diǎn)的距離。如果面源或體源為非正常排放源，可按上述修正后，采用相應(yīng)的非正常排放點(diǎn)源模式。但當(dāng)風(fēng)速 U10 1.5m/s 時(shí)，對(duì)于體源，可用在實(shí)際的時(shí)刻 t 中加一個(gè)初始時(shí)間 t0的方法進(jìn)行修正：3/10220103 . 41zyxa

25、aat式中，ax、ay、az分別為體源在平均風(fēng)向(x)、橫向(y)、垂直方向(z)的邊長(zhǎng)。1.4面源(體源)的數(shù)值積分算法將面源（體源）分成一系列微小面源（dLdW）,分別以點(diǎn)源代替這些微源，再用相應(yīng)點(diǎn)源模式計(jì)算每一個(gè)微源對(duì)預(yù)測(cè)點(diǎn)的影響，再疊加所有微源在測(cè)點(diǎn)的濃度。這種算法思維直接明了，算法準(zhǔn)確，但無(wú)法用手工完成，以前計(jì)算機(jī)速度太慢也不現(xiàn)實(shí)，而現(xiàn)在的 PC 機(jī)的計(jì)算速度已經(jīng)允許采用二次積分方法來實(shí)現(xiàn)這一算法，同時(shí)控制較高的精度。設(shè)面源 是一封閉的區(qū)間，坐標(biāo)系是以風(fēng)向?yàn)檎?X 軸，幾何中心為原點(diǎn)，其幾何中心坐標(biāo)為（0,0） 。預(yù)測(cè)點(diǎn)在下風(fēng)向X，橫風(fēng)向?yàn)?Y，測(cè)點(diǎn)坐標(biāo)(X,Y)。對(duì) 內(nèi)的任意一個(gè)點(diǎn)

26、(x,y)，設(shè)其代表一個(gè)面積為 dxdy的微源，則該微源對(duì)預(yù)測(cè)點(diǎn)的濃度貢獻(xiàn)為：C(X,Y)(x,y)=Qdxdy f(X-x,Y-y)式中，Q 是面源單位面積的源強(qiáng)mg/(s.m2)，Qdxdy 為微元的源強(qiáng)(mg/s)。f()為點(diǎn)源計(jì)算公式，與風(fēng)速、排放是否連續(xù)、穩(wěn)定度等因素有關(guān)。因此整個(gè)面源對(duì)（X,Y）的濃度等于 內(nèi)所有點(diǎn)對(duì)（X,Y）濃度的疊加值，用積分式表示為：QdxdyyYxXfYXC),(),(f(X-x,Y-y)代表了面源中點(diǎn) x,y 對(duì)測(cè)點(diǎn)（X,Y）的點(diǎn)源計(jì)算公式。這個(gè)公式可以是有風(fēng)模式，也可以是小風(fēng)靜風(fēng)模式；可以是正常持續(xù)排放模式，也可是非正常排放的間斷源模式。還可以對(duì)源有效高

27、和源強(qiáng)進(jìn)行適當(dāng)修正（如干沉、濕沉、化學(xué)遷移等等） ，總之，可以完全運(yùn)用點(diǎn)源計(jì)算中的公式，詳細(xì)見點(diǎn)源章節(jié)。因 Q 實(shí)際上是與位置有關(guān)的（若要考慮源衰減） ，因此將Q 歸入 f()函數(shù)中，形成一個(gè)新的函數(shù) F()，則 C（X,Y）的濃度為：dxdyyYxXFYXC),(),(本軟件中對(duì)上式采用變步長(zhǎng)辛普生二重積分法進(jìn)行積分，計(jì)算速度約為點(diǎn)源修正方法的百分之一。如果采用 P300 以上處理器，這一速度還是可以接受的。但若計(jì)算長(zhǎng)期平均濃度時(shí)，一般 P計(jì)算機(jī)顯得太慢。為簡(jiǎn)化求積過程和方便定義面源形狀，本軟件中所有面源（體源）均要求是矩形，但角度可作任意旋轉(zhuǎn)。優(yōu)點(diǎn)：可以準(zhǔn)確描述面源與風(fēng)向的關(guān)系，即使面源的

28、邊與風(fēng)向不垂直時(shí)也能較好計(jì)算，面源內(nèi)部的點(diǎn)也能算；可以直接使用點(diǎn)源公式，小風(fēng)靜風(fēng)時(shí)也能計(jì)算。從理論上來說，計(jì)算結(jié)果也比一般方法準(zhǔn)確。缺點(diǎn)：計(jì)算量很大；要求準(zhǔn)確定義面源位置和形狀。1.5多源疊加模式如果需要評(píng)價(jià)的點(diǎn)源多于一個(gè)，計(jì)算濃度時(shí)，應(yīng)將各個(gè)源對(duì)接受點(diǎn)濃度的貢獻(xiàn)進(jìn)行疊加。在評(píng)價(jià)區(qū)內(nèi)選一原點(diǎn)，以平均風(fēng)的上風(fēng)方為正 x 軸，各個(gè)源（坐標(biāo)為 xr,yr,0）對(duì)評(píng)價(jià)區(qū)內(nèi)任一地面點(diǎn)（x,y）的濃度總貢獻(xiàn) Cn可按下式計(jì)算：)529 . 4.(.,0 ,rrrnyyxxCyxC式中 Cr是第 r 個(gè)點(diǎn)源對(duì)（x,y,0）點(diǎn)的濃度貢獻(xiàn)，其計(jì)算公式可根據(jù)不同條件選用本章給出的有關(guān)點(diǎn)源模式，但是注意坐標(biāo)變換，

29、（x,y,0）代以（x-xr，y-yr） 。1.6長(zhǎng)期平均濃度計(jì)算方法1.6.1孤立源長(zhǎng)期平均濃度公式當(dāng)平均時(shí)間超過 1 小時(shí)之后，由于風(fēng)向的擺動(dòng)，任一風(fēng)方位內(nèi)的污染物濃度在橫向都將趨于均勻分布。為此，可將連續(xù)點(diǎn)源模式對(duì) y 積分，并除以接受點(diǎn)所在位置的風(fēng)方位寬度（或弧線長(zhǎng)度） 。對(duì)于孤立排放源，以煙囪地面位置為原點(diǎn)，在某一穩(wěn)定度（序號(hào)為 j）和平均風(fēng)速（序號(hào)為 k）時(shí)，任意風(fēng)向方位 i 的下風(fēng)方 x 處的長(zhǎng)期平均濃度（季、期或年均值）Cijk（x）（mgm3）為： FnxuQCzijk12/3/2式中 n風(fēng)向方位數(shù)，一般取 16（2x/16 即為 x 處22.5圓心角對(duì)應(yīng)的弧長(zhǎng)） ；其他符號(hào)

30、同前。在可能出現(xiàn)的穩(wěn)定度和平均風(fēng)速條件下，任意風(fēng)向位 i 的下風(fēng)方 x 處的長(zhǎng)期平均濃度 Ci(x)(mg m3)為： LijkLijkkijkijkkiifCfCxC式中 fijk有風(fēng)時(shí)風(fēng)向方位、穩(wěn)定度、風(fēng)速聯(lián)合頻率； Cijk對(duì)應(yīng)于該聯(lián)合頻率在下風(fēng)方 x 處有風(fēng)時(shí)的濃度值，由式（4.9-35）給出； fLijk靜風(fēng)或小風(fēng)時(shí)，不同風(fēng)方位和穩(wěn)定度的出現(xiàn)頻率（下標(biāo) k 只含有靜風(fēng)和小風(fēng)兩個(gè)風(fēng)速段） ；CLijk對(duì)應(yīng)于fLijk的靜風(fēng)或小風(fēng)時(shí)的地面濃度。因?yàn)殪o風(fēng)或小風(fēng)時(shí)的風(fēng)脈動(dòng)角本來就比較大，CLijk可直接按小風(fēng)靜風(fēng)模式計(jì)算。式（4.9-36）中的 j 和 k 的加總數(shù)取決于所劃分的穩(wěn)定度和風(fēng)速

31、段數(shù)目，j 的總數(shù)不宜少于 3（穩(wěn)定、中性、不穩(wěn)定） ；有風(fēng)時(shí) k 的總數(shù)一般也宜少于 3。在估算每個(gè)風(fēng)速段的平均風(fēng)速時(shí)，由于平均風(fēng)速出現(xiàn)在公式的分母中，因而平均風(fēng)速應(yīng)等于單次風(fēng)速倒數(shù)的平均值倒數(shù)。其表達(dá)式為：1/1/1iiuNu式中 ui第 i 個(gè)風(fēng)速值； N總個(gè)數(shù)。1.6.2多源長(zhǎng)期平均濃度公式如果評(píng)價(jià)區(qū)的煙囪多于一個(gè)，則任一接受點(diǎn)（x，y）的長(zhǎng)期平均濃度為：LijkrLijkrijkrijkrkjifCfCyxC,式中 Crijk、CLrijk分別是在接受點(diǎn)上風(fēng)方對(duì)應(yīng)于fijk和fLijk聯(lián)合頻率的第 r 個(gè)源對(duì)接收點(diǎn)的濃度貢獻(xiàn)。Crijk、CLrijk的公式形式分別和 Cijk、CLi

32、jk相同（參閱 4.9.6.1） ，但應(yīng)注意坐標(biāo)變換，將坐標(biāo)轉(zhuǎn)換到以接受點(diǎn)為原點(diǎn)，i 風(fēng)方位為正 x 軸的新坐標(biāo)系后，再應(yīng)用 Cijk或 CLijk公式。1點(diǎn)源擴(kuò)散模式1.1持續(xù)排放源1.1.1有風(fēng)模式（U101.5m/s）1.1.2小風(fēng)靜風(fēng)模式（U101.5m/s）1.2非正常排放源1.2.1有風(fēng)情況（U101.5m/s）1.2.2小風(fēng)靜風(fēng)（U101.5m/s）1.3單源擴(kuò)散的地面軸線最大濃度1.4對(duì)源強(qiáng)和有效源高的修正1.4.1地形對(duì)煙羽的影響1.4.2熱浮力對(duì)煙羽的修正1.4.3干沉積修正1.4.4濕沉積修正1.4.5化學(xué)遷移的修正2面源(體源)的點(diǎn)源修正算法3面源(體源)的數(shù)值積分算法

33、4多源疊加模式5長(zhǎng)期平均濃度計(jì)算方法5.1孤立源長(zhǎng)期平均濃度公式5.2多源長(zhǎng)期平均濃度公式2 點(diǎn)源擴(kuò)散模式2.1持續(xù)排放源2.1.1有風(fēng)模式（U101.5m/s）1.自由空間中的連續(xù)點(diǎn)源實(shí)際上絕大多數(shù)污染源都是連續(xù)的，對(duì)于連續(xù)排放源，可理解為在時(shí)間上依次連續(xù)釋放無(wú)窮多個(gè)煙團(tuán)。因此，連續(xù)排放源的擴(kuò)散模式可以通過將瞬時(shí)單煙團(tuán)模式對(duì) to從-到 t 積分后求得。以煙團(tuán)初始空間坐標(biāo)為原點(diǎn)，下風(fēng)方為 x 軸，煙羽軸線與 x 軸一直保持重合，都是 x 的函數(shù)，將對(duì) to的積分變zyx,換為對(duì)（x-uT）/x的積分，可得最基本的煙羽擴(kuò)散模式：22222exp2exp2,zyzyzyuQzyxC適用條件為：自

34、由空間；風(fēng)速要比較大（u101.5m/s） ；當(dāng)大氣不穩(wěn)定狀態(tài)時(shí)，可能帶來一定的誤差。2.地面反射用像源法，假想地平線為一鏡面，在其下方有一與真實(shí)源完全對(duì)稱的虛源，則這兩個(gè)源疊加后的效果和真實(shí)源考慮到地面反射的結(jié)果是等價(jià)的。以煙囪地面位置的中心點(diǎn)為坐標(biāo)原點(diǎn)，在考慮到地面反射后，污染源下風(fēng)方任一點(diǎn)小于 24 小時(shí)取樣時(shí)間的污染物濃度 C（x，y，z）由下式給出：2222222exp2exp2exp2,zezeyzyHzHzyuQzyxCz=0 時(shí)的地面濃度 C（x,y,0） ，可簡(jiǎn)化為；222222exp0 ,zeyzyHyuQyxC下風(fēng)方 X 軸線上（y=0）的地面濃度 C（x,0,0）為：2

35、2212exp0 , 0 ,zeHuQxC對(duì)于較低的排放源（例如 He50m，具體限值由地面粗糙度、混合層高度等因素決定） ，一般可直接應(yīng)以上式子計(jì)算。3.混合層反射對(duì)于高架源，當(dāng)超過一定的下風(fēng)距離時(shí)，需對(duì)煙羽在混合層頂?shù)姆瓷溥M(jìn)行修正。同考慮地面反射類似，用像源法修正后，污染源下風(fēng)方任一點(diǎn)小于 24 小時(shí)取樣時(shí)間的污染物濃度C（x,y,z）可表示為：FyuQzyxcyzy222exp2),(kknzezezHnhzHnhF222222exp22exp式中 h-混合層高度； k-反射次數(shù)，一、二級(jí)項(xiàng)目取 k=4 已足夠。4.側(cè)面反射詳見狹長(zhǎng)山谷擴(kuò)散模式。2.1.2小風(fēng)靜風(fēng)模式（U101.5m/s

36、）小風(fēng)靜風(fēng)時(shí)，污染物地面濃度 C(x,y,0)可用下式計(jì)算： GQYXcL2022/322,式中 和 G 按下式計(jì)算：，2202201222eHYX seseGU2/32/2201221， dtestsoo2/22101UXS 和分別是橫向和鉛直向擴(kuò)散參數(shù)的回歸系數(shù)(y=z=0102T,z=T),T為擴(kuò)散時(shí)間(s),和的定值見 HJT2.2-199301020102附錄 B3。注意，上式中 He是煙筒有效高度相對(duì)于預(yù)測(cè)點(diǎn)的高度。若預(yù)測(cè)點(diǎn)高度坐標(biāo)為 z，煙筒有效高度處坐標(biāo)為 z0，則 He=z0-z。2.2非正常排放源非正常排放是指建設(shè)項(xiàng)目生產(chǎn)運(yùn)行階段的開車、停車、檢修、一般性事故和發(fā)生漏泄等情

37、況時(shí)的污染物的不正常排放。非正常排放常發(fā)生在有限時(shí)間（T）內(nèi)。以瞬時(shí)單煙團(tuán)正態(tài)擴(kuò)散式，對(duì) t0在有限時(shí)間 T 內(nèi)積分，經(jīng)整理后可得非正常排放模式。2.2.1有風(fēng)情況（U101.5m/s）非正常排放條件下的地面濃度ca(mg/m3)建議按下列各式計(jì)算。以排氣筒地面位置為原點(diǎn),有效源高為He，平均風(fēng)向軸為X軸，源強(qiáng)為Q（mg/s） ，開始非正常排放時(shí)的時(shí)間為 t，非正常排放持續(xù)時(shí)間為T，預(yù)測(cè)時(shí)刻的時(shí)間為t。1.有風(fēng)情況（U101.5m/s）t時(shí)刻任一點(diǎn)（x,y,z）的濃度，以持續(xù)排放源模式為基礎(chǔ)，乘上一個(gè)系數(shù) G1，按下式計(jì)算：1222exp2),(GFyuQzyxcyzykknzezezHnhz

38、HnhF222222exp22exp11xxxxxxUtxUTUtxUtG式中 F-混合層反射項(xiàng)； G1-非正常排放項(xiàng)； h-混合層高度； k-反射次數(shù)，一、二級(jí)項(xiàng)目取 k=4 已足夠。擴(kuò)散參數(shù)各指數(shù)、系數(shù)的定值見導(dǎo)則,2121XXzyx附錄 B。2.2.2小風(fēng)靜風(fēng)（U10=1.5m/s）的高斯持續(xù)排放點(diǎn)源，并且要求穩(wěn)定度較不穩(wěn)定、混合層反射可忽略等條件，其計(jì)算結(jié)果 Xm必須在擴(kuò)散參數(shù)系數(shù) y1、y2和指數(shù) a1、a2的應(yīng)用范圍之內(nèi)。例如，以 1000m 范圍內(nèi)的擴(kuò)散參數(shù)計(jì)算系數(shù)和指數(shù)計(jì)算得Xm=3000m，則結(jié)果是不可靠的。由于解析式應(yīng)用范圍有限，應(yīng)用條件較苛刻，實(shí)際計(jì)算中常用數(shù)值法。借助計(jì)

39、算機(jī)的快速計(jì)算，可使用分段逼近法求最大落地濃度。該原理是：將最大濃度可能的出現(xiàn)范圍，如下風(fēng)向0，1000000m，分成 10 段共 11 個(gè)點(diǎn)位，每點(diǎn)求一濃度，再以最大濃度點(diǎn)位左、右兩點(diǎn)位作為新的計(jì)算范圍起、止點(diǎn)，再分10 段，如此循環(huán)計(jì)算，直到每段長(zhǎng)度小于要求的精度（如0.01m） 。一般計(jì)算幾十次即可得到結(jié)果。這種方法算法很簡(jiǎn)單，無(wú)須求函數(shù)導(dǎo)數(shù)，同時(shí)靈活性極強(qiáng)，因?yàn)閷?duì)每一個(gè)點(diǎn)求濃度的時(shí)候，可以按調(diào)用任意模式按任意方法計(jì)算。因此，這一方法可用于有風(fēng)、小風(fēng)靜風(fēng)、面源體源的點(diǎn)源修正法和非正常排放等所有單源模式的求最大濃度值。2.4對(duì)源強(qiáng)和有效源高的修正2.4.1地形對(duì)煙羽的影響此前的擴(kuò)散模式都假

40、設(shè)地面是完全平整的（煙囪底部是一個(gè)無(wú)限大的水平面，其高程為 0） ，因此在擴(kuò)散過程中煙羽的中心線可保持水平不變。但如果在預(yù)測(cè)點(diǎn)（x,y,z）處，地面有一定的高程 hT(0z)，則在對(duì)（x,y,z）式應(yīng)用以上模式時(shí)，應(yīng)對(duì)有效煙羽高度進(jìn)行一些修正。假定煙羽路徑始終與起伏的地形保持平行，或者假設(shè)煙羽軸線保持固定的海拔高度，并與高于煙羽的地形相交，都是不正確的，實(shí)際情況應(yīng)該是介于上述二者之間。具體的修正方法如下。（1）中性和不穩(wěn)定天氣條件令：hT為凸出的地形高度；He為煙軸高度（即有效高度） ；T為煙軸高度修正系數(shù)（或地形系數(shù)） ，修正后的煙囪有效高度應(yīng)該是 THe oT 則應(yīng)按下式取值：時(shí)當(dāng)時(shí)當(dāng)Te

41、eTeTehHHhHThHT,/ )2/(, 2/1（2）穩(wěn)定天氣條件（D-E、E、F）在穩(wěn)定天氣條件下，當(dāng)煙羽逼近孤立山體時(shí)，煙羽以臨界高度 Hc為界分成兩部分，臨界高度以上的煙羽有足夠的動(dòng)能爬越山體，而臨界高度以下的煙羽，只能被迫繞著山體過去。臨界高度 H c 可由下式確定：Hc=Hm-u/(gd/dz)1/2式中 Hm-孤立山體高度，m； Hc-臨界高度，m； -z 高度處大氣位溫，K； d/dz-z 高度處位溫梯度，K/m； u-平均風(fēng)速，m/s； g-重力加速度，m/s2。例如，=300K，d/dz=0.01K/m，u=2m/s，Hm=200m；則Hc=Hm-111=89m，煙囪有效

42、高度大于 89m 時(shí)，煙羽將有足夠的動(dòng)能爬越山體。對(duì)于山體高度 Hm已定的情況，大氣越穩(wěn)定，則 Hc 越小。所以一般只需計(jì)算在 F 穩(wěn)定度下的 Hc，如果煙羽有效高度HeHc(F)，則可認(rèn)為煙羽能夠爬越山體。2.4.2熱浮力對(duì)煙羽的修正1.混合層頂穿透問題高架源熱浮力煙羽對(duì)混合層并非只有“完全穿透”和“完全不穿透”兩種情況，還須考慮到“部分穿透”問題。定義 P 為穿透系數(shù)：P=1.5-(h-H)/H當(dāng) P0 時(shí)，H1=minH,2(h-H)/3,Q1=Q；當(dāng) P 1 時(shí)，H1=h-H,Q1=0；當(dāng) 0P1 時(shí)，H1=(2/3+P/3)(h-H),Q1=Q(1-P)式中，Q1，H1-分別為修正后

43、的源強(qiáng)和抬升高度； h-混合層高度。2.風(fēng)向切變和熱浮力問題高浮力煙羽在 z 方向受風(fēng)向切變及熱浮力的影響可按下式對(duì)y和 z進(jìn)行修正：2yc=2y+(H/3.5)2+(0.032x2)2zc=2z+(H/3.5)2式中，yc、zc-分別為修正后 y、z 方向的擴(kuò)散參數(shù)； (H/3.5)2-熱浮力修正項(xiàng)； (0.032x2)-風(fēng)向切變修正項(xiàng)。 按下表確定：表 風(fēng)向切變修正參數(shù) /Z（度/m）穩(wěn)定度ABCDEF/Z0.0050.0100.0150.0200.0250.0352.4.3干沉積修正1.干沉積的含義在重力、湍流擴(kuò)散、分子擴(kuò)散、靜電引力以及其他生物學(xué)、化學(xué)和物理學(xué)等因素的作用下，大氣中的顆

44、粒物或某些氣體隨時(shí)會(huì)被地表（土壤、植物、水體）滯留或吸收，使這些物質(zhì)連續(xù)不斷地從大氣向地表作質(zhì)量轉(zhuǎn)移，從而減少其在空氣中的濃度。通常，把這一與降水作用無(wú)關(guān)的質(zhì)量轉(zhuǎn)移過程，稱之為干沉積。在對(duì)中距離大氣污染物輸送的研究中發(fā)現(xiàn)，差不多一半以上的質(zhì)量轉(zhuǎn)移是干沉積過程引起的。為些，對(duì)一些物理、化學(xué)特性和空氣相差較大的大氣污染物（例如顆粒物） ，或者對(duì)于輸送距離比較遠(yuǎn)的高架源的濃度預(yù)測(cè)，都應(yīng)當(dāng)考慮其干沉積的影響。下面將介紹兩種常用的干沉積擴(kuò)散模式。2.源虧損模式源虧損模式主要用于粒徑小于 10m 易產(chǎn)生沉積的顆粒物或氣體。假定因各種機(jī)理造成的大氣污染物向地面的沉積通量Wmg/(sm2)由下式表示：W=Vd

45、C.(4.9-79)式中 Vd-沉積速度，m/s； C-大氣污染物地面濃度，常取自距地面 1m 高度處。大氣污染物自煙囪出口排出后，其初始源強(qiáng) Q（0）因沉積作用將隨下風(fēng)距離逐漸減弱（虧損） 。根據(jù)式（4.9-79）可導(dǎo)出虧損后的源強(qiáng) Q（x） ；因?yàn)閷?W 從到對(duì) y 積分恰是-dQ(x)/dx，再對(duì) x 積分則可得 dxHuVQxQxzezd02212exp2exp0由式（4.9-80）可見，問題的關(guān)鍵是給出沉積速度 Vd。長(zhǎng)期以來，許多研究者為求得 Vd，在實(shí)驗(yàn)或模式化方面都進(jìn)行了大量的工作。對(duì)于 SO2，Hanna 曾推薦如下 Vd(cm/s)的野外實(shí)驗(yàn)值：高約 1m 的莊稼地、城市、

46、水面為 0.7；高 0.1m 的草地為 0.5；酸性土壤為 0.4（干）或 0.6（濕） ；石灰質(zhì)土壤為 0.8；干積雪為 0.1。在模式化方面，通常是把沉積速度和電流類比，假定沉積速度與各種阻力的總和成反比。上式可用數(shù)值積分計(jì)算。采用 1 小時(shí)取樣時(shí)間的國(guó)標(biāo)擴(kuò)散參數(shù)計(jì)算結(jié)果表明，D 穩(wěn)定度時(shí)虧損最明顯，不穩(wěn)定或更穩(wěn)定時(shí)虧損依次減少。設(shè) He=100m,Q（0）為 100，Vd為 0.007m/s，U為 3m/s，則在 10000m 處 Q 為下表所示：穩(wěn)定度ABB-CCC-DDD-EEFQ(10km)99.0896.2795.1793.9392.2191.9092.8194.7698.86一

47、般來說，3000m 以內(nèi)的擴(kuò)散源虧損小于 1%，可以忽略不計(jì)。3.部分反射（斜煙羽）部分反射模式主要用于粒徑大于 10m 的氣載顆粒物。因?yàn)檫@種污染物有明顯的沉降作用，其煙羽重心是是逐漸降低的，此外，地面只能反射一部分的污染物。對(duì)某一粒徑的顆粒污染物，對(duì) He 和 Q 進(jìn)行以下調(diào)整后，再按相應(yīng)氣態(tài)物模式計(jì)算：He(x)=He(0)-Vg x/u。Q(x)=Q(0) (1+)/2式中 -顆粒物的地面反射系數(shù)（參閱下表） ； Vg-顆粒物沉降速度，用 STOCKS 公式計(jì)算：；182gdVg d、-分別為顆粒物的直徑和密度； g-重力加速度； -空氣粘性系數(shù)。表 地面反射系數(shù) 粒徑范圍/m1530

48、31476076100平均粒徑/m22385085反射系數(shù) 0.80.50.302.4.4濕沉積修正濕沉積系指大氣污染物因降水而減少其在空氣中濃度的過程。嚴(yán)格地講，濕沉積分云中和云下兩種清除機(jī)制。在工程應(yīng)用中，常把這兩種機(jī)制合起來考慮。假設(shè)大氣污染物的初始源強(qiáng)Q（0）因降水隨下風(fēng)距離 x 成指數(shù)衰減，則Q（x）=Q(0)exp(-x/u)式中 x-接受點(diǎn)的下風(fēng)距離； u-煙囪出口處的風(fēng)速；-清除系數(shù)，S1。令 J 為降水強(qiáng)度（mm/h） ，根據(jù)實(shí)驗(yàn)定測(cè)量結(jié)果，對(duì)于SO2=1.7104J0.62.4.5化學(xué)遷移的修正化學(xué)遷移的修正，在工程應(yīng)用上與濕沉積類似，可令修正后的源強(qiáng) Q（x）等于初始源強(qiáng)

49、 Q（0）乘以修正因子fc 。ccuTxf/exp 式中 Tc-大氣污染物的時(shí)間常數(shù)；其他符號(hào)同前。如果定義fc=1/2 所對(duì)應(yīng)的時(shí)間為該大氣污染物的半衰期td，將其代入式（44.9-85） ，可得：dcutxf/693. 0exp 對(duì)于城市尺度，SO2的 td4 小時(shí)；長(zhǎng)距離輸送 SO2的 td典型值為 4 天。3 面源(體源)的點(diǎn)源修正算法把面源的排放當(dāng)作一個(gè)位于其幾何中心的點(diǎn)源的排放，對(duì)擴(kuò)散參數(shù)適當(dāng)修正后，采用點(diǎn)源模式直接計(jì)算，用以近似代表該面源的擴(kuò)散。但這種方法僅適用于有風(fēng)條件（U101.5m/s） 。在小風(fēng)靜風(fēng)的條件下，面源擴(kuò)散如何計(jì)算，尚未有定論，這里暫按小風(fēng)靜風(fēng)的點(diǎn)源方法計(jì)算。如

50、果測(cè)點(diǎn)位于面源（體源）之外，這種近似算法較為準(zhǔn)確，但如果測(cè)點(diǎn)在面源內(nèi)部，則計(jì)算結(jié)果就很不理想。比如，在有風(fēng)時(shí)，位于面源幾何中心上風(fēng)向各點(diǎn)（仍處于面源內(nèi)部） ，按這種近似算法濃度應(yīng)是 0，實(shí)際不然。此外，這種算法不能反映面源的形狀。這些都是點(diǎn)源修正算法的不足之外。有風(fēng)條件下，對(duì)于面源或體源內(nèi)部的預(yù)測(cè)點(diǎn)，可采用下式計(jì)算：,2jsQC 式中，x 為測(cè)點(diǎn)離面源上風(fēng)邊界的距離。這種方axH2222/法使得在有風(fēng)時(shí)，面源內(nèi)任意一點(diǎn)均有濃度。面源。采用直接修正法。如果面源的面積較?。⊿1km2） ，面源外的 Cs可按點(diǎn)源擴(kuò)散模式計(jì)算，只是應(yīng)附加一個(gè)初始擾動(dòng)。這一初始擾動(dòng)使煙羽在 x=0 就有一個(gè)和面源橫向?qū)?/p>

51、度相等的橫向尺度，以及和面源高度相等的垂直向尺度。注意到煙羽的半寬度等于 2.15y或 2.15z，則修正后的 y和 z分別為：15. 2/3 . 4/2211Hxaxazyay式中 x-自接受點(diǎn)至面源中心點(diǎn)的距離； ay-面源在 y 方向的長(zhǎng)度； H-面源的平均排放高度。體源當(dāng)無(wú)組織排放源為體源時(shí)，地面濃度可按點(diǎn)源直接修正法計(jì)算；令 y、z分別表示體源在 y 和 z 方向的邊長(zhǎng)，則修正后的 y和 z分別為：3 . 4/3 . 4/2211zzyyxxx-自接受點(diǎn)至面源中心點(diǎn)的距離。如果面源或體源為非正常排放源，可按上述修正后，采用相應(yīng)的非正常排放點(diǎn)源模式。但當(dāng)風(fēng)速 U10 1.5m/s 時(shí)，

52、對(duì)于體源，可用在實(shí)際的時(shí)刻 t 中加一個(gè)初始時(shí)間 t0的方法進(jìn)行修正：3/10220103 . 41zyxaaat式中，ax、ay、az分別為體源在平均風(fēng)向(x)、橫向(y)、垂直方向(z)的邊長(zhǎng)。4 面源(體源)的數(shù)值積分算法將面源（體源）分成一系列微小面源（dLdW）,分別以點(diǎn)源代替這些微源，再用相應(yīng)點(diǎn)源模式計(jì)算每一個(gè)微源對(duì)預(yù)測(cè)點(diǎn)的影響，再疊加所有微源在測(cè)點(diǎn)的濃度。這種算法思維直接明了，算法準(zhǔn)確，但無(wú)法用手工完成，以前計(jì)算機(jī)速度太慢也不現(xiàn)實(shí)，而現(xiàn)在的 PC 機(jī)的計(jì)算速度已經(jīng)允許采用二次積分方法來實(shí)現(xiàn)這一算法，同時(shí)控制較高的精度。設(shè)面源 是一封閉的區(qū)間，坐標(biāo)系是以風(fēng)向?yàn)檎?X 軸，幾何中心為

53、原點(diǎn)，其幾何中心坐標(biāo)為（0,0） 。預(yù)測(cè)點(diǎn)在下風(fēng)向X，橫風(fēng)向?yàn)?Y，測(cè)點(diǎn)坐標(biāo)(X,Y)。對(duì) 內(nèi)的任意一個(gè)點(diǎn)(x,y)，設(shè)其代表一個(gè)面積為 dxdy的微源，則該微源對(duì)預(yù)測(cè)點(diǎn)的濃度貢獻(xiàn)為：C(X,Y)(x,y)=Qdxdy f(X-x,Y-y)式中，Q 是面源單位面積的源強(qiáng)mg/(s.m2)，Qdxdy 為微元的源強(qiáng)(mg/s)。f()為點(diǎn)源計(jì)算公式，與風(fēng)速、排放是否連續(xù)、穩(wěn)定度等因素有關(guān)。因此整個(gè)面源對(duì)（X,Y）的濃度等于 內(nèi)所有點(diǎn)對(duì)（X,Y）濃度的疊加值，用積分式表示為：QdxdyyYxXfYXC),(),(f(X-x,Y-y)代表了面源中點(diǎn) x,y 對(duì)測(cè)點(diǎn)（X,Y）的點(diǎn)源計(jì)算公式。這個(gè)公式可

54、以是有風(fēng)模式，也可以是小風(fēng)靜風(fēng)模式；可以是正常持續(xù)排放模式，也可是非正常排放的間斷源模式。還可以對(duì)源有效高和源強(qiáng)進(jìn)行適當(dāng)修正（如干沉、濕沉、化學(xué)遷移等等） ，總之，可以完全運(yùn)用點(diǎn)源計(jì)算中的公式，詳細(xì)見點(diǎn)源章節(jié)。因 Q 實(shí)際上是與位置有關(guān)的（若要考慮源衰減） ，因此將Q 歸入 f()函數(shù)中，形成一個(gè)新的函數(shù) F()，則 C（X,Y）的濃度為：dxdyyYxXFYXC),(),(本軟件中對(duì)上式采用變步長(zhǎng)辛普生二重積分法進(jìn)行積分，計(jì)算速度約為點(diǎn)源修正方法的百分之一。如果采用 P300 以上處理器，這一速度還是可以接受的。但若計(jì)算長(zhǎng)期平均濃度時(shí)，一般 P計(jì)算機(jī)顯得太慢。為簡(jiǎn)化求積過程和方便定義面源形狀

55、，本軟件中所有面源（體源）均要求是矩形，但角度可作任意旋轉(zhuǎn)。優(yōu)點(diǎn)：可以準(zhǔn)確描述面源與風(fēng)向的關(guān)系，即使面源的邊與風(fēng)向不垂直時(shí)也能較好計(jì)算，面源內(nèi)部的點(diǎn)也能算；可以直接使用點(diǎn)源公式，小風(fēng)靜風(fēng)時(shí)也能計(jì)算。從理論上來說，計(jì)算結(jié)果也比一般方法準(zhǔn)確。缺點(diǎn)：計(jì)算量很大；要求準(zhǔn)確定義面源位置和形狀。5 多源疊加模式如果需要評(píng)價(jià)的點(diǎn)源多于一個(gè)，計(jì)算濃度時(shí)，應(yīng)將各個(gè)源對(duì)接受點(diǎn)濃度的貢獻(xiàn)進(jìn)行疊加。在評(píng)價(jià)區(qū)內(nèi)選一原點(diǎn)，以平均風(fēng)的上風(fēng)方為正 x 軸，各個(gè)源（坐標(biāo)為 xr,yr,0）對(duì)評(píng)價(jià)區(qū)內(nèi)任一地面點(diǎn)（x,y）的濃度總貢獻(xiàn) Cn可按下式計(jì)算：)529 . 4.(.,0 ,rrrnyyxxCyxC式中 Cr是第 r 個(gè)點(diǎn)

56、源對(duì)（x,y,0）點(diǎn)的濃度貢獻(xiàn)，其計(jì)算公式可根據(jù)不同條件選用本章給出的有關(guān)點(diǎn)源模式，但是注意坐標(biāo)變換， （x,y,0）代以（x-xr，y-yr） 。6 長(zhǎng)期平均濃度計(jì)算方法6.1孤立源長(zhǎng)期平均濃度公式當(dāng)平均時(shí)間超過 1 小時(shí)之后，由于風(fēng)向的擺動(dòng)，任一風(fēng)方位內(nèi)的污染物濃度在橫向都將趨于均勻分布。為此，可將連續(xù)點(diǎn)源模式對(duì) y 積分，并除以接受點(diǎn)所在位置的風(fēng)方位寬度（或弧線長(zhǎng)度） 。對(duì)于孤立排放源，以煙囪地面位置為原點(diǎn)，在某一穩(wěn)定度（序號(hào)為 j）和平均風(fēng)速（序號(hào)為 k）時(shí)，任意風(fēng)向方位 i 的下風(fēng)方 x 處的長(zhǎng)期平均濃度（季、期或年均值）Cijk（x）（mgm3）為： FnxuQCzijk12/3/

57、2式中 n-風(fēng)向方位數(shù)，一般取 16（2x/16 即為 x 處22.5圓心角對(duì)應(yīng)的弧長(zhǎng)） ；其他符號(hào)同前。在可能出現(xiàn)的穩(wěn)定度和平均風(fēng)速條件下，任意風(fēng)向位 i 的下風(fēng)方 x 處的長(zhǎng)期平均濃度 Ci(x)(mg m3)為： LijkLijkkijkijkkiifCfCxC式中 fijk-有風(fēng)時(shí)風(fēng)向方位、穩(wěn)定度、風(fēng)速聯(lián)合頻率； Cijk-對(duì)應(yīng)于該聯(lián)合頻率在下風(fēng)方 x 處有風(fēng)時(shí)的濃度值，由式（4.9-35）給出； fLijk-靜風(fēng)或小風(fēng)時(shí)，不同風(fēng)方位和穩(wěn)定度的出現(xiàn)頻率（下標(biāo) k 只含有靜風(fēng)和小風(fēng)兩個(gè)風(fēng)速段） ；CLijk-對(duì)應(yīng)于fLijk的靜風(fēng)或小風(fēng)時(shí)的地面濃度。因?yàn)殪o風(fēng)或小風(fēng)時(shí)的風(fēng)脈動(dòng)角本來就比較大

58、，CLijk可直接按小風(fēng)靜風(fēng)模式計(jì)算。式（4.9-36）中的 j 和 k 的加總數(shù)取決于所劃分的穩(wěn)定度和風(fēng)速段數(shù)目，j 的總數(shù)不宜少于 3（穩(wěn)定、中性、不穩(wěn)定） ；有風(fēng)時(shí) k 的總數(shù)一般也宜少于 3。在估算每個(gè)風(fēng)速段的平均風(fēng)速時(shí)，由于平均風(fēng)速出現(xiàn)在公式的分母中，因而平均風(fēng)速應(yīng)等于單次風(fēng)速倒數(shù)的平均值倒數(shù)。其表達(dá)式為：1/1/1iiuNu式中 ui-第 i 個(gè)風(fēng)速值； N-總個(gè)數(shù)。6.2多源長(zhǎng)期平均濃度公式如果評(píng)價(jià)區(qū)的煙囪多于一個(gè)，則任一接受點(diǎn)（x，y）的長(zhǎng)期平均濃度為：LijkrLijkrijkrijkrkjifCfCyxC,式中 Crijk、CLrijk-分別是在接受點(diǎn)上風(fēng)方對(duì)應(yīng)于fijk和

59、fLijk聯(lián)合頻率的第 r 個(gè)源對(duì)接收點(diǎn)的濃度貢獻(xiàn)。Crijk、CLrijk的公式形式分別和 Cijk、CLijk相同（參閱 4.9.6.1） ，但應(yīng)注意坐標(biāo)變換，將坐標(biāo)轉(zhuǎn)換到以接受點(diǎn)為原點(diǎn)，i 風(fēng)方位為正 x 軸的新坐標(biāo)系后，再應(yīng)用 Cijk或 CLijk公式。公路汽車尾氣預(yù)測(cè)模式1.1.1.1概述1.1.1.2污染物排放強(qiáng)度1.1.1.3坐標(biāo)變換與計(jì)算的起、止點(diǎn)1.1.1.4HIWAY-2 模式1.1.1.5CALINE4 模式1.1.1.6小風(fēng)靜風(fēng)模式1.1.1.7垂直風(fēng)和平行風(fēng)時(shí)的解析式6.2.1.1 概述對(duì)汽車尾氣的擴(kuò)散預(yù)測(cè)，以高斯擴(kuò)散模式為基礎(chǔ)，各國(guó)曾推導(dǎo)出多種實(shí)用的氣態(tài)污染物擴(kuò)散

60、模式，如我國(guó)常用的近似式、內(nèi)插式，美國(guó) EPA 的 HIWAY-2、加州運(yùn)輸部的 CALINE4，得克薩斯州的 TXLINE 和英國(guó)的簡(jiǎn)單桌面模式等。經(jīng)過監(jiān)測(cè)、驗(yàn)算和對(duì)比，除內(nèi)插式和桌面式差別較大外，對(duì)平原微丘地區(qū)的直線公路，其它模式的計(jì)算結(jié)果相差并不大。這里采用了 HIWAY-2 積分模式和 CALINE4 模式。HIWAY-2 模式的基本思想是：將公路線源劃分成一系列微元，每一微元可作為點(diǎn)源處理，以高斯點(diǎn)源擴(kuò)散公式計(jì)算每一微元對(duì)預(yù)測(cè)點(diǎn)的貢獻(xiàn)值，然后線性疊加，結(jié)果作為公路線源對(duì)預(yù)測(cè)點(diǎn)的形成濃度的近似值。而 CALINE4 模式的基本思想是：將公路線源劃分成一系列線元，每一線元看作一個(gè)中點(diǎn)在線