高斯煙羽模型

1、模型假設(shè)：1、坐標(biāo)系高斯模型的坐標(biāo)系如圖2.1所示，原點(diǎn)為排放點(diǎn)若為高架源，原點(diǎn)為排放 點(diǎn)在地面的投影，x軸正向?yàn)轱L(fēng)速方向，y軸在水平面上垂直于x軸，正向在x 軸的左側(cè)，z軸垂直于水平面xoy，向上為正向。在此坐標(biāo)系下煙流中心線或煙 流中心線在xoy面的投影與x軸重合。圖2,1高斯模式的坐標(biāo)系2、模型假設(shè)污染物的濃度在y、z軸上的分布是高斯分布正態(tài)分布的；污染源的源強(qiáng)是連續(xù)且均勻的，初始時(shí)刻云團(tuán)內(nèi)部的濃度、溫度呈均勻 分布；擴(kuò)散過程中不考慮云團(tuán)內(nèi)部溫度的變化，忽略熱傳遞、熱對(duì)流與熱輻射;泄漏氣體是理想氣體，遵守理想氣體狀態(tài)方程；在水平方向，大氣擴(kuò)散系數(shù)呈各向同性；取x軸為平均風(fēng)速方向，整個(gè)擴(kuò)散

2、過程中風(fēng)速的大小、方向保持不變， 不隨地點(diǎn)、時(shí)間變化而變化；地面對(duì)泄漏氣體起全反射作用，不發(fā)生吸收或吸附作用；整個(gè)過程中，泄漏氣體不發(fā)生沉降、分解，不發(fā)生任何化學(xué)反應(yīng)等。3、模型公式推導(dǎo)由正態(tài)分布假設(shè)可以導(dǎo)出下風(fēng)向任意一點(diǎn)Xx,y,z處泄漏氣體濃度的函 數(shù)為：1X(x, y, z) = A(x)e ay 2 e bz21由概率統(tǒng)計(jì)理論可以寫出方差的表達(dá)式為:js y 2 Xdyj Xdyz 2 Xdzjs Xdz由假設(shè)可以寫出源強(qiáng)的積分公式:Q = j-s式中：b y、b z為泄漏氣體在y、z方向分布的標(biāo)準(zhǔn)差，單位為m；Xx,y,z為任一點(diǎn)處泄漏氣體的濃度，單位為kg/m3； u為平均風(fēng)速，單

3、位為m/s； Q為源強(qiáng)即泄漏速度，單位為kg/s；將1式代入2式，1a =將1式代入2式，1a =2b 2yb =上 2b 2z將1式和4式代入3式，積分可得:人（x）= Q2兀ub b,5式代入1式，可得：再將4式和積分可得:4X(X(x, y, z)=exp2兀ub b上式為無界空間連續(xù)點(diǎn)源擴(kuò)散的高斯模型公式，然而在實(shí)際中，由于地面的 存在，煙羽的擴(kuò)散是有界的。根據(jù)假設(shè)可以把地面看做一鏡面，對(duì)泄漏氣體起全 反射作用，并采用像源法處理，原理如圖2.2所示。可以把任一點(diǎn)p處的濃度看 做兩部分的貢獻(xiàn)之和：一部分是不存在地面時(shí)所造成的泄漏物濃度；一部分是由 于地面反射作用增加的泄漏物濃度。該處的泄

4、漏物濃度即相當(dāng)于不存在地面時(shí)由 位于0,0,H的實(shí)源和位于0,0,-H的像源在P點(diǎn)處所造成的泄漏物濃度之 和。圖2.2像源法原理示意圖其中，實(shí)源的貢獻(xiàn)為：Q / 1 y 21 (z - H )2Xy，z) =exp(-2旅)exp(-2-)y zyz其中，像源的貢獻(xiàn)為：/、 Q / 1 y21(z + h)2X*y，z)= 2Sexp(-2ge碩-2-) y zyz則該處的實(shí)際濃度為：X(x, y, z)= X(x, y, z)+ X x, y, z)由以上條件公式可得到高架連續(xù)點(diǎn)源擴(kuò)散的高斯煙羽模型公式為：Q z 1 y2exp( - _ 二)x TOC o 1-5 h z 2兀 u b b

5、2 b 2( 1 (z - H)2: 1 (z + H)2：1exp( -) + exp(-) HYPERLINK l bookmark5 o Current Document 2 b 22 b 2其中：Xx,y,z為下風(fēng)向x米、橫向y米、地面上方z米處的擴(kuò)散的氣體 濃度，單位為kg/m3； Q為源強(qiáng)即源釋放速率，單位為kg/s； u為平均風(fēng)速， 單位為峪；。為水平擴(kuò)散參數(shù)，單位加；。為垂直擴(kuò)散參數(shù)，單位為m；七為泄漏后是時(shí)間，單位為s； H為泄漏源有效高度，單位為m； y為橫向距離，單 位為m； z為垂直方向距離，單位為m。如：式中，令z=0，即可得到地面氣體濃度計(jì)算公式：X(x, X(x,

6、 y,0, H)= Q 兀ub blz y、exp -()22 b1-y(H 2exp 2b 2z(11)令y=z=0，即可得到地面軸線氣體濃度計(jì)算公式:QH 2X (x,0,0, H) = -_2- exp(- -)12y zz其中，Xx,y,0為下風(fēng)向x米、橫向y米處的地面擴(kuò)散氣體濃度，單位為 kg/m3；若令y=0，則可以得到下風(fēng)向中心線上的濃度分布。4、泄漏源有效高度煙云抬升高度的計(jì)算：以上式中的泄漏源有效高度是指泄漏氣體形成的氣云基本上變成水平狀的 時(shí)候氣云中心的離地高度。實(shí)際上，泄漏源有效高度就等于泄漏源幾何高度加泄 漏煙云抬升高度。影響煙云抬升高度的因素有很多，主要包括：泄漏氣體

7、的初始速度和方向、 初始溫度、泄漏口直徑、環(huán)境風(fēng)速與風(fēng)速歲高度的變化率、環(huán)境溫度與大氣穩(wěn)定 度。有效源高：H = HS +AH其中，H s為泄漏源幾何架高，AH為煙云抬升高度。 實(shí)驗(yàn)表明，泄漏源抬升高度可以用下面公式近似計(jì)算：AH = 2.4V d / V13其中，,是氣云釋放速度，單位為m/s； d是泄漏出口直徑，單位為m； V 為環(huán)境風(fēng)速，單位為m/s；13式是20世紀(jì)70年代末、80年代初，Wilson根據(jù)管道破裂泄漏實(shí)驗(yàn) 所得的實(shí)驗(yàn)公式。實(shí)驗(yàn)時(shí)氣體的噴射方向與風(fēng)速垂直并且氣體噴射路徑上無障礙 物。實(shí)驗(yàn)表明，當(dāng)氣體噴射方向垂直向上時(shí)，預(yù)測(cè)值與實(shí)際值之比在洗漱2以內(nèi)。計(jì)算出泄漏煙云抬升高度

8、以后，將泄漏源抬升高度與泄漏源實(shí)際幾何高度相 加就得到了泄漏源有效高度。5、擴(kuò)散系數(shù)的選?。簲U(kuò)散系數(shù)b x、b y、b z擴(kuò)散系數(shù)b x、按照Pasquill的分類方法，隨著氣象條件穩(wěn)定性的增加，大氣穩(wěn)定度可以分為A、B、C、D、E、F六類。其中A、B、C三類表示氣象條件不穩(wěn)定，E、F兩 類表示氣象條件穩(wěn)定，D類表示中性氣象條件，也就是說氣象條件的穩(wěn)定性在穩(wěn) 定和不穩(wěn)定之間。A、B、C三種類型的穩(wěn)定度中，A類表示氣象條件極其不穩(wěn)定， B類表示氣象條件中等程度不穩(wěn)定，C類表示氣象條件弱不穩(wěn)定。E和F兩種類 型的穩(wěn)定度中，E類表示氣象條件弱穩(wěn)定，F(xiàn)類表示氣象條件中等程度穩(wěn)定。大 氣穩(wěn)定度具體分類方

9、法見下表1、表2。表1 Pasquill大氣穩(wěn)定度的確定地面風(fēng)速m/s白天日照夜間條件強(qiáng)中等弱陰天且云 層薄，或 低空云量 為4/8天空云量 為3/8 6CDDDD表2日照強(qiáng)度的確定天空云層的情況日照角 60日照角35日照角15。且35天空云量為4/8，或高空有薄云強(qiáng)中等弱天空云量為 5/8-7/8，云層高 度為 2134-4877m中等弱弱天空云量為 5/8-7/8，云層高 度2134m弱弱弱表1和表2中的云量是指當(dāng)?shù)靥炜諏痈采w率。例如，云量為3/8是指當(dāng)?shù)?3/8的天空有云層覆蓋。日照角是指當(dāng)?shù)靥柟饩€與地平線之間的夾度。例如， 陽光垂直照射地賣弄時(shí)的日照角為90。一般來說，隨著大氣穩(wěn)定度的增加，擴(kuò)散系數(shù)減小。根據(jù)Hanna和Drivas 的建議，化學(xué)危險(xiǎn)品事故泄漏擴(kuò)散系數(shù)與大氣穩(wěn)定度類型和下風(fēng)向的關(guān)系如下表 3。表3擴(kuò)散系數(shù)的計(jì)算方法太氣穩(wěn)定度A0.22a/(I 十叫OOLr)妲OJj：B0.1女/（1+。頌吹）”D.2aCOHj