高斯煙羽模型

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

2、方向，整個擴散過程中風(fēng)速的大小、方向保持不變，不隨地點、時間變化而變化；(7)地面對泄漏氣體起全反射作用，不發(fā)生吸收或吸附作用；(8)整個過程中，泄漏氣體不發(fā)生沉降、分解，不發(fā)生任何化學(xué)反應(yīng)等。3、模型公式推導(dǎo) 由正態(tài)分布假設(shè)可以導(dǎo)出下風(fēng)向任意一點x（x,y,z）處泄漏氣體濃度的函數(shù)為： （1） 由概率統(tǒng)計理論可以寫出方差的表達式為： （2）由假設(shè)可以寫出源強的積分公式： （3）式中：、為泄漏氣體在y、z方向分布的標(biāo)準(zhǔn)差，單位為 m；x（x,y,z）為任一點處泄漏氣體的濃度，單位為 kg/m；u為平均風(fēng)速，單位為 m/s；q為源強（即泄漏速度），單位為 kg/s；將（1）式代入（2）式，積分可

3、得： （4）將（1）式和（4）式代入（3）式，積分可得： （5）再將（4）式和（5）式代入（1）式，可得： （6）上式為無界空間連續(xù)點源擴散的高斯模型公式，然而在實際中，由于地面的存在，煙羽的擴散是有界的。根據(jù)假設(shè)可以把地面看做一鏡面，對泄漏氣體起全反射作用，并采用像源法處理，原理如圖2.2所示?？梢园讶我稽cp處的濃度看做兩部分的貢獻之和：一部分是不存在地面時所造成的泄漏物濃度；一部分是由于地面反射作用增加的泄漏物濃度。該處的泄漏物濃度即相當(dāng)于不存在地面時由位于（0,0,h）的實源和位于（0,0,-h）的像源在p點處所造成的泄漏物濃度之和。其中，實源的貢獻為： （7）其中，像源的貢獻為： （8

4、）則該處的實際濃度為： （9）由以上條件公式可得到高架連續(xù)點源擴散的高斯煙羽模型公式為： （10）其中：x（x,y,z）為下風(fēng)向x米、橫向y米、地面上方z米處的擴散的氣體濃度，單位為kg/m；q為源強（即源釋放速率），單位為kg/s；u為平均風(fēng)速，單位為m/s；為水平擴散參數(shù)，單位為m；為垂直擴散參數(shù)，單位為 m；t為泄漏后是時間，單位為s；h為泄漏源有效高度，單位為m；y為橫向距離，單位為m；z為垂直方向距離，單位為m。如：式中，令z=0，即可得到地面氣體濃度計算公式： (11)令y=z=0，即可得到地面軸線氣體濃度計算公式： （12）其中，x（x,y,0）為下風(fēng)向x米、橫向y米處的地面擴散

5、氣體濃度，單位為kg/m；若令y=0，則可以得到下風(fēng)向中心線上的濃度分布。4、泄漏源有效高度（煙云抬升高度的計算）：以上式中的泄漏源有效高度是指泄漏氣體形成的氣云基本上變成水平狀的時候氣云中心的離地高度。實際上，泄漏源有效高度就等于泄漏源幾何高度加泄漏煙云抬升高度。影響煙云抬升高度的因素有很多，主要包括：泄漏氣體的初始速度和方向、初始溫度、泄漏口直徑、環(huán)境風(fēng)速及風(fēng)速歲高度的變化率、環(huán)境溫度及大氣穩(wěn)定度。有效源高： 其中，為泄漏源幾何架高，為煙云抬升高度。實驗表明，泄漏源抬升高度可以用下面公式近似計算： （13）其中，是氣云釋放速度，單位為m/s；是泄漏出口直徑，單位為 m；v為環(huán)境風(fēng)速，單位為

6、m/s；（13）式是20世紀(jì)70年代末、80年代初，wilson根據(jù)管道破裂泄漏實驗所得的實驗公式。實驗時氣體的噴射方向與風(fēng)速垂直并且氣體噴射路徑上無障礙物。實驗表明，當(dāng)氣體噴射方向垂直向上時，預(yù)測值與實際值之比在洗漱2以內(nèi)。計算出泄漏煙云抬升高度以后，將泄漏源抬升高度與泄漏源實際幾何高度相加就得到了泄漏源有效高度。 5、擴散系數(shù)的選?。簲U散系數(shù)、的大小與大氣湍流結(jié)構(gòu)、離地高度、地面粗糙度、泄漏持續(xù)時間、抽樣時間間隔、風(fēng)速以及離開泄漏源的距離等因素有關(guān)。大氣的湍流結(jié)構(gòu)和風(fēng)速在大氣穩(wěn)定度中考慮。大氣穩(wěn)定度由10米高度以上的風(fēng)速、白天的太陽輻射或夜間的云量等參數(shù)決定。按照pasquill的分類方法

7、，隨著氣象條件穩(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類表示氣象條件中等程度穩(wěn)定。大氣穩(wěn)定度具體分類方法見下表1、表2。表1 pasquill大氣穩(wěn)定度的確定地面風(fēng)速（m/s）白天日照夜間條件強中等弱陰天且云層薄，或低空云量為4/8天空云量為3/8aa-bbcca-bbb-cc-ddbccddedddfedd表2 日照強度的確定天空云層的情況日照角 60日照角 35日照角 15且 35天空云量為4/8，或高空有薄云強中等弱天空云量為5/8-7/8，云層高度為2134-4877m中等弱弱天空云量為5/8-7/8，云層高度2134m弱弱弱 表1和表2中的云量是指當(dāng)?shù)靥炜諏痈采w率。例如，云量為3/8是指當(dāng)?shù)?/8的天空有云層覆蓋。日照角是指當(dāng)?shù)靥?/p>