B大氣污染評價模型-夏深寶、李品、胡世鋒

1、 /242015年科技大學(xué)模擬訓(xùn)練二承諾書我們仔細閱讀了數(shù)學(xué)建模選拔賽的規(guī)則.我們完全明白，在做題期間不能以任何方式（包括、電子、網(wǎng)上咨詢等）與隊外的任何人研究、討論與選拔題有關(guān)的問題。我們知道，抄襲別人的成果是違反選拔規(guī)則的,如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們重承諾，嚴格遵守選拔規(guī)則，以保證選拔的公正、公平性。如有違反選拔規(guī)則的行為，我們將受到嚴肅處理。我們選擇的題號是（從A/B/C/D中選擇一項填寫）：B隊員簽名：1.2.3.日期：2015/8/212015年科技大學(xué)數(shù)學(xué)建模競賽選拔編號專用頁評閱

2、編號（評閱前進行編號）：評閱記錄（評閱時使用）： /24關(guān)于大氣污染預(yù)報問題的數(shù)學(xué)建模摘要近年來，伴隨著社會經(jīng)濟的迅速發(fā)展，但經(jīng)濟增長模式相對傳統(tǒng)落后，對生態(tài)平衡和自然環(huán)境造成一定破壞，因此，控制大氣污染，提高環(huán)境空氣污染質(zhì)量就成了迫切的需要。本文正是以京津冀地區(qū)的空氣污染問題為背景，針對單污染源空氣污染擴散和多污染源空氣污染擴散分別建立不同的模型來進行研究的。對于問題1，我們借鑒國家空氣質(zhì)量標(biāo)準，并采用加權(quán)平均的思想引入了綜合空氣質(zhì)量指數(shù)來評價某個地區(qū)的空氣質(zhì)量環(huán)境綜合空氣質(zhì)量指數(shù)數(shù)值越大，則說明綜合污染程度越重。對于問題2，根據(jù)京津冀地區(qū)的污染數(shù)據(jù)統(tǒng)計與分析，最后將大氣污染源分為三大類：工

3、業(yè)源，生活源，交通運輸源，并分別對京津冀三個地區(qū)的主要污染物和他們的產(chǎn)業(yè)結(jié)構(gòu)做研究和分析。對于問題3，采用了經(jīng)典的高斯模式來建立單污染源空氣污染擴散模型，然后用Matlab模擬出工廠分別在早上8點，中午12點，晚上9點空氣污染濃度分布和空氣質(zhì)量等級。對于問題4，先建立了多污染源空氣污染擴散模型，然后以汽車尾氣污染源為例，利用線源模式對汽車尾氣擴散濃度分布進行預(yù)測，最后將我們用模型模擬出來的結(jié)果與實際結(jié)果進行對比，結(jié)果發(fā)現(xiàn)誤差不是很大，說明模型比較符合實際。最后，針對京津冀地區(qū)空氣污染的具體情況和對自己建立模型的分析，給京線源模式津冀地區(qū)環(huán)保部門提出了實現(xiàn)“APEC”的可行性措施和建議。關(guān)鍵詞：

4、綜合空氣質(zhì)量指數(shù)污染源高斯模式Matlab一、問題重述大氣是指包圍在地球外圍的空氣層，是地球自然環(huán)境的重要組成部分之一。人類生活在大氣里，潔凈大氣是人類賴于生存的必要條件。一個人在五個星期不吃飯或5天不喝水，尚能維持生命，但超過5分鐘不呼吸空氣，便會死亡。隨著地球上人口的急劇增加，人類經(jīng)濟增長的急速增大，對生態(tài)平衡和自然環(huán)境造成一定的破壞，地球上的大氣污染日趨嚴重，其影響也日趨深刻，如由于一些有害氣體的大量排放，不僅造成局部地區(qū)大氣的污染，而且影響到全球性的氣候變化?？諝馕廴镜谋撞∪找嫱怀?，特別是日益加重的霧霾天氣已經(jīng)干擾到社會的出行秩序和生活質(zhì)量。因此，加強大氣質(zhì)量的監(jiān)測和預(yù)報是非常必要。目

5、前對大氣質(zhì)量的監(jiān)測主要是監(jiān)測大氣中5、X、CO、；二、懸浮顆粒物（主要為PM10）等的濃度，研究說明，城市空氣質(zhì)量好壞與季節(jié)與氣象條件的關(guān)系十分密切。以城市京、津、冀地區(qū)的污染物含量與氣象參數(shù)的數(shù)據(jù)來解決以下問題：（1）參考現(xiàn)有國標(biāo)和美標(biāo)，建立衡量空氣質(zhì)量優(yōu)劣程度等級的數(shù)學(xué)模型。（2）查找數(shù)據(jù)并列出京津冀地區(qū)主要污染源與其污染參數(shù)，分析影響空氣質(zhì)量的主要污染源的性質(zhì)和種類。（3）建立單污染源空氣污染擴散模型，描述其對周圍空氣污染的動態(tài)影響規(guī)律?，F(xiàn)有境某一工廠廢氣排放煙囪高50m,主要排放物為氮氧化物。早上9點至下午3點期間的排放濃度為406.92mg/m3，排放速度為1200m3/h；晚上10

6、點-凌晨4點期間的排放濃度為1160mg/m3,排放速度為5700m3/h；通過你的擴散模型求解該工廠方圓51公里分別在早上8點、中午12點、晚上9點空氣污染濃度分布和空氣質(zhì)量等級。（4）建立多污染源空氣污染擴散模型，并以汽車尾氣污染源為例求解分析以下問題：在2015年1月15日已經(jīng)連續(xù)三天發(fā)生重污染，假設(shè)從16日開始啟動汽車單雙號限行交通管制措施，求解市二環(huán)、四環(huán)、六環(huán)路在16日早上8點、中午12點、晚上9點時空氣污染濃度梯度變化與空氣質(zhì)量等級。（5）根據(jù)你們的模型和求解結(jié)果，分析總結(jié)影響空氣質(zhì)量的關(guān)鍵參數(shù)，為京津冀地區(qū)環(huán)保部門撰寫一份建議報告，給出實現(xiàn)“APEC”的可行性措施和建議。二、模

7、型假設(shè)1、各組數(shù)據(jù)真實可信，各參數(shù)不隨測量地點的改變而改變，且是在同一地點同一時間采集，不考慮人為因素，具有統(tǒng)計、預(yù)測意義。2、假設(shè)京津冀三個地區(qū)發(fā)展?fàn)顩r相同，各個城市不會出現(xiàn)重大空氣污染事故，或環(huán)境不會變好，環(huán)境不會出現(xiàn)強烈波動，即發(fā)展速度沒有明顯差異。3、各城市的檢測時間具有隨機性，城市空氣質(zhì)量好壞與氣象條件的關(guān)系十分密切，與其他因素關(guān)系不大。4、污染物在擴散過程中，污染物本身是被動、保守的，即污染物和空氣無相對運動。5、擴散時擴散面平坦、開闊、性質(zhì)均勻，平均流場平直、穩(wěn)定，不考慮風(fēng)場對擴散的影響。6、假設(shè)啟動汽車單雙號限行交通管制措施后，車輛數(shù)較少為原來的一半，忽略其他干擾因素。7、通常

8、情況下，車輛在有坡度的道路上行駛，其排出的污染物要比平路上的高假設(shè)制定排放標(biāo)準時不考慮道路坡度的問題。三、符號說明C：污染物項目i的空氣質(zhì)量分指數(shù)iI：污染物項目i的空氣質(zhì)量分指數(shù)iBI：與C相近的污染物濃度限值的高位值HiBI：與C相近的污染物濃度限值的低位值LiHI：與BI對應(yīng)的空氣質(zhì)量分指數(shù)HLI：與BI對應(yīng)的空氣質(zhì)量分指數(shù)Lw:第i個目標(biāo)污染物的權(quán)系數(shù)CAQI：綜合空氣質(zhì)量指數(shù)C（x,y,z）:空間點（x,y,z）的污染物濃度，mg/sQ：源強，單位時間污染物排放量，mg/sU風(fēng)速平均速度m/sV：顆粒的重力沉降速度，單位m/sta顆粒的地面反射系數(shù)gg:煙氣的擴散系數(shù)，與大氣穩(wěn)定度和

9、水平距離X有關(guān)xyd顆粒直徑mpp：顆粒密度Kg/mp卩：空氣粘度Pa*sg:重力加速度m/丁v:煙流出口處速度，m/sgD:煙囪出口徑，mT:煙囪出口處的煙氣溫度，KgT:環(huán)境大氣溫度，KQ:某條線源道路，j類型車p型燃料w種污染物排放源強，g/hijpwL:第i段路長，Kmiq:j類型車p型燃料在j段路上的車流量/hijpn：某條線源道路上劃分的總段數(shù)C:j類型車每輛每公里p型燃料消耗量，Kg/KmjpE:p類型燃料，w種污染物的排放因子，Kg/1fpw四、問題分析4.1問題一的分析查找現(xiàn)有國家標(biāo)準，考慮國情，以國標(biāo)為基準找到國家大氣質(zhì)量標(biāo)準中環(huán)境空氣質(zhì)量的標(biāo)準規(guī)定，進一步分析可建立衡量空

10、氣質(zhì)量優(yōu)劣程度等級的數(shù)學(xué)模型。4.2問題二的分析查找、三個代表城市一年主要污染源與其污染參數(shù)，用Excel處理數(shù)據(jù)并根據(jù)問題立的模型進行分析。4.3問題三的分析對于問題3建立單污染源空氣污染擴散模型，污染濃度的擴散模式是一種以處理大氣污染物在大氣中輸送和擴散問題的數(shù)學(xué)模型，由于影響擴散過程的氣象條件、地形、下墜面狀況與污染本身的復(fù)雜性，到目前為止，還不能找到一個適用于各種條件的大氣擴散模式來描述所有這些復(fù)雜條件下的大氣擴散問題，所以我們必須建立一個比較研究污染物從高處向四周擴散的模型，然后將所研究的工廠廢氣排放煙囪看作是高架點源。利用Matlab編程，將題中數(shù)據(jù)帶入公式即可求得空氣污染濃度分布

11、和空氣質(zhì)量等級。4.4問題四的分析問題4建立多污染源空氣污染擴散模型，可以對問題3建立的單污染源進行疊加。針對汽車尾氣污染源，可以將汽車看成是一個一個的污染源，將其轉(zhuǎn)換成多污染源問題，再利用多污染源對其進行求解。4.5問題五的分析根據(jù)分析京津冀三個地區(qū)主要污染源的性質(zhì)和種類與其擴散模型和求解結(jié)果，可以給出合理化措施和建議。五、模型建立與求解衡量空氣質(zhì)量優(yōu)劣程度等級模型模型建立我國現(xiàn)在通常采用AQI來衡量空氣質(zhì)量，AQI計算方便，根據(jù)對人體健康影響最大的某項污染物來衡量空氣質(zhì)量狀況，有助于人們直觀的了解空氣污染狀況。但當(dāng)空氣中各類污染物質(zhì)量濃度差異較大時，AQI會掩蓋污染物對空氣的總體影響狀況，

12、丟失大量有用信息，這對研究各類污染物較為不利。因此可進一步研究空氣質(zhì)量評價存在的模糊和灰色區(qū)域。根據(jù)環(huán)境空氣質(zhì)量標(biāo)準當(dāng)中的所規(guī)定的污染物排放限值，依據(jù)加權(quán)平均的思想來建立衡量空氣質(zhì)量優(yōu)良等級的評價模型。在本模型中，我們引入了綜合空氣質(zhì)量指數(shù)來評價某個地區(qū)的空氣質(zhì)量，指數(shù)數(shù)值越大說明綜合污染程度越重。首先假設(shè)空氣中只有一種污染物i,根據(jù)空氣質(zhì)量分指數(shù)的計算方法，我們得到污染物i的空氣質(zhì)量分指數(shù)HI-LISIh-BIl實際空氣中含有多種污染物，當(dāng)對多種污染物進行綜合評價時，我們要對各種污染物值分別進行加權(quán)平均，從而確定評價指標(biāo)的權(quán)重值，此權(quán)重是衡量因子集中某一因子對空氣污染影響程度相對大小的量，權(quán)

13、重系數(shù)越大，則該污染對空氣的影響程度越大。本文以相近的污染物濃度限值的平均值來決定權(quán)重，各污染物的監(jiān)測值相對于大氣質(zhì)量標(biāo)準限額的超標(biāo)越大，則對污染的貢獻越大，從而權(quán)重也就越大，因此權(quán)重系數(shù)可以表示為：2)其中P表示污染物i的權(quán)重值,B表示與Cj相近的污染物濃度限值的平均值,得到了權(quán)重值后，為方便進行復(fù)合運算，還應(yīng)該對n個污染物因子的權(quán)重結(jié)果進行歸一化處理，即第i個目標(biāo)污染物的權(quán)系數(shù)可以表示為：j二丁I二亍:二匸CifE；(3)最后，根據(jù)模糊運算關(guān)系，用每個污染物的空氣質(zhì)量分指數(shù)乘以相應(yīng)的權(quán)重系數(shù)就可以得到衡量空氣質(zhì)量的參數(shù)CAQICAQI=vr=1i.5.1.2空氣質(zhì)量指數(shù)級別的劃分類似于國家

14、空氣質(zhì)量指數(shù)級別的劃分方法，根據(jù)中國環(huán)境監(jiān)測總站發(fā)布的2015年3月份京津冀地區(qū)PM2.5、PM10、5：、:、C0-95per和O3-8H-90per六種污染物的月均濃度數(shù)據(jù)，大致算出它們分別的空氣質(zhì)量分指數(shù)，然后將他們累加得到綜合空氣質(zhì)量指數(shù)，再根據(jù)環(huán)境監(jiān)測總站發(fā)布的空氣級別對其進行修正。根據(jù)修正后得到的綜合空氣質(zhì)量指數(shù)也將空氣質(zhì)量劃分為六個級別，依次為優(yōu)、良、輕度污染、中度污染、重度污染和嚴重污染。表1綜合空氣質(zhì)量指數(shù)級別與相關(guān)信息綜合空氣質(zhì)量指數(shù)空氣質(zhì)量指數(shù)級別空氣質(zhì)量指數(shù)類別與表示顏色對健康影響情況建議采取的措施0-40一級優(yōu)綠色無影響各類人群可正?；顒?0-85二級良黃色可承受，對

15、少數(shù)異常敏感人群健康有影響極少數(shù)敏感人群減少戶外活動85T40三級輕度污染橙色易感人群癥狀輕度加劇，健康人群出現(xiàn)刺激癥狀易感人群減少時間、高強度的戶外鍛煉140-195四級中度污染紅色進一步加劇易感人群癥狀，對健康人呼吸系統(tǒng)有影響易感人群應(yīng)避免高強度的戶外鍛煉，一般人群適量減少戶外活動195-270五級重度污染紫色心臟病和肺病患者癥狀顯著加劇，健康人群出現(xiàn)癥狀易感人群停止戶外活動，一般人群減少戶外活動270六級嚴重污染褐紅色健康人群運動耐受力降低，有明顯強烈癥狀易感人群停留在戶，一般人群避免戶外活動5.1.3首要污染物與超標(biāo)污染物的確定方法二：對空氣污染程度最大的污染物i;n：污染物的項目Fj

16、50時，I】；最大的污染物為首要污染物；1：-最大的污染物為兩項或兩項以上時，則他們并列為首要污染物;口140的污染物稱為超標(biāo)污染物。5.2京津冀地區(qū)主要污染源與其污染參數(shù)PM2.5的排放源主要有：揚塵、交通源、燃煤源、建筑塵、生物質(zhì)燃燒和工業(yè)無組織排放等。此外，還存在5門、：二與有機物轉(zhuǎn)化形成顆粒物的二次排放源。揚塵和燃煤污染源貢獻率分別為18.11%和16.14%，是PM2.5的主要來源。4.50%3305.60%比例未丸源有機物燃煤機動車排旅生物質(zhì)燃姥二次硫酸鹽和硝酸鹽圖1京津冀PM2.5主要來源的年均貢獻率數(shù)據(jù)顯示，2012年，京津冀燃煤消費總量38927萬噸，煤炭消費量占其能源消費總

17、量的88.8%，遠遠高于的4.3%和的6.8%，煤炭消費排放出大量的二氧化硫，對大氣環(huán)境造成很大影響。圖2京津冀燃煤消費總量比例冀區(qū)氮氧化物的比重達45%，分別高于28.8個百分點，26.2個百分點。比例26%北京45%天津石家莊29%圖3京津冀機動車氮氧化物排放占比根據(jù)2015年3月與第一季度京津冀空氣質(zhì)量報告，我們統(tǒng)計出京津冀地區(qū)3月份PM2.5、PM10、二氧化氮、二氧化硫、一氧化碳、臭氧等污染物月均濃度分布表。圖4京津冀2015年3月份各污染物月均濃度通過上面的分析，我們把大氣污染源分三類：（1）工業(yè)源：工業(yè)生產(chǎn)是大氣污染的一個重要來源，工業(yè)生產(chǎn)排放到大氣中的污染物種類繁多，有煙塵，硫

18、的氧化物，氮的氧化物，有機化合物，鹵化物，碳化合物等，有的是煙塵，有的是氣體。（2）生活爐灶與采暖鍋爐源：城市量民用生活爐灶和采暖鍋爐需要消耗大量煤炭，煤炭在燃燒的過程中要釋放大量的煙塵，二氧化硫，一氧化碳等有害物質(zhì)污染大氣。特別是在冬季采暖時，造成很大污染。城市和人口密集的居住區(qū)是人類消費活動集中地，是主要的生活污染源。（3）交通運輸源：汽車、火車、飛機、輪船是當(dāng)代主要的運輸工具，它們燒煤或石油產(chǎn)生的廢氣也是重要的污染源，特別是城市中的汽車，量大而集中，尾氣所排放的污染物能直接侵襲人的呼吸器官，對城市的空氣污染很嚴重，成為大城市空氣的主要污染源之一，汽車排放的廢氣主要有一氧化碳，二氧化硫，氮

19、氧化物和碳氫化合物等，前三種物質(zhì)危害很大。5.3單污染源空氣擴散模型污染濃度的擴散模式是一種用以處理大氣污染物在大氣中輸送和擴散問題的物理和數(shù)學(xué)模型。由于影響擴散過程的氣象條件、地形、下墊面狀況與污染本身的復(fù)雜性，我們采用了經(jīng)典的高斯模式來建立單污染源空氣污染擴散模型。高斯擴展模式是在大量的實測資料的分析基礎(chǔ)上，應(yīng)用湍流統(tǒng)計理論得到正態(tài)分布假設(shè)下的擴散模式。一般總是把排放口或高架源在地面上的投影點作為坐標(biāo)原點，x軸正向沿平均風(fēng)向水平延伸，y軸在水平面上垂直于x軸，x軸左側(cè)為正,z軸垂直于水平面，向上為正，形成右手坐標(biāo)系。模型的建立連續(xù)點源一般指排放大量污染物的煙囪，放散管，通風(fēng)口等，本模型就是

20、針對于地面連續(xù)點源或高架點源。針對排放口處于高空位置的高架點源，將點源在地面上的投影點作為坐標(biāo)原點，有效源位于Z軸上某點，Z=H,高架點源有效的高度有兩部分組成，H二h+二其中h為排放口的有效高度，二；是熱煙流的浮升力和煙氣以一定的速度豎直離開排放口的沖力使煙流抬升的一個附加高度。假設(shè)污染物到達地面后被完全吸收，那污染物的濃度為：C&y,z,H)-益旳卜扌看詈”若要計算高架點源的地面濃度公式，則可令z=0,得：C(x,y,0,H)=exp-1gj(7)煙囪、放散管一般為顆粒物，要根據(jù)其粒徑大小來判斷其計算方法，若顆粒的粒徑小于15：;,其地面濃度可按上述氣體擴散模式計算，而對于大于15：;的顆

21、粒物，由于具有明顯的重力沉降作用，將使?jié)舛戎涤兴淖?，可按傾斜煙流模式計算濃度：嗎一exp四其中，匚顆粒的地面反射系數(shù)可查表：表2粒子圍反射系數(shù)粒子圍/pm.15-3031-4748-7576-100平均粒徑/pm.22386085反射系數(shù)肛0.80.50.30式中的表示顆粒的重力沉降速度，單位m/s,計算公式:煙流抬升高度與擴散參數(shù)的確定高斯擴散公式應(yīng)用的效果依賴于公式中各個參數(shù)的準確程度，尤其是煙流抬升高度H與擴散參數(shù)的估算，其中，平均風(fēng)速u取觀測的常規(guī)氣象數(shù)據(jù)，源強Q可以計算或測定，即H與氣象條件和地面狀況密切相關(guān)。煙流抬升高度的計算煙流抬升高度是確定高架源的位置，準確判斷大氣污染擴散與

22、估計地面污染濃度的重要參數(shù)之一。連續(xù)點源的排放大部分是采用煙囪排放的，通常具有一定速度的熱煙氣從煙囪排出口后，會繼續(xù)上升到很高的高度，相當(dāng)于增加了煙囪的幾何高度。圖5煙氣抬升示意圖因此煙囪的有效高度H,應(yīng)該為煙囪的幾何高度h與煙氣抬升高度二.之和。即：H二h+二】(10)對于一個煙囪來說，h是固定的，只要能計算出二】，就可以確定H的數(shù)值了。產(chǎn)生煙氣抬升有兩方面原因：一是由于煙囪的出口煙氣具有一定的初始動量；二是由于煙氣溫度高于周圍氣溫而產(chǎn)生一定的浮力。初始動量的大小決定于煙氣出口流速和煙囪出口徑，而浮力大小則主要決定于煙氣與周圍大氣之間的溫差。此外，平均風(fēng)速，風(fēng)速垂直切變以與大氣穩(wěn)定度等因素，

23、都會對煙氣抬升產(chǎn)生影響。由于煙流抬升受諸多因素的影響，二：的廣泛計算方法用霍蘭德公式：Ah=號(1.5+2.7=扌(1.5VSD+9.6宜103Qh)(11)上式的計算結(jié)果對很強的熱源比較適中甚至偏高，而對中小型熱(QhuOOMW)的估計偏低。當(dāng)大氣處于不穩(wěn)定或穩(wěn)定狀態(tài)時，可在上式計算結(jié)果的基礎(chǔ)上分別增加或減少10%-20%。擴散參數(shù)的估算擴散參數(shù)丁工-是表示擴散圍與速率大小的特征量，即正態(tài)分布函數(shù)的標(biāo)準差，可采用以下經(jīng)驗公式確定：=(12)7.=(13)表3橫向擴散參數(shù)冪函數(shù)表達式系數(shù)值二=匸：-（取樣時間是0.5h）穩(wěn)定度Yi下風(fēng)距離，mA0.9010740.4258090-10000.8

24、509340.6020521000B0.9143700.2818460-10000.8750860.3963531000B-C0.9193250.2295000-10000.8750860.3142381000C0.9242790.1771541-10000.8851570.2321231000C-D0.9268490.1439401-10000.8869400.1893961000D0.9294180.1107261-10000.8887230.1466691000D-E0.9251180.09856311-10000.8927940.1243081000E0.9208180.0864001

25、1-10000.8968640.1019471000F0.9294180.05536340-10000.8887230.07333481000表4垂直擴散參數(shù)冪函數(shù)表達式系數(shù)值二、；：-：（取樣時間是0.5h）穩(wěn)定度下風(fēng)距離，m1.121540.07999040-300A1.513600.00854771300-5002.108810.000211545500B0.9644350.1271900-5001.093560.057025500B-C0.9410150.1146820-5001.007700.075718500C0.9175950.1068030C-D0.8386280.126152

26、0-20000.7564100.2356672000-100000.8155750.136659100000.8262120.1046341-1000D0.6320230.4001671000-100000.555360.810763100000.7768640.1117710-2000D-E0.5723470.52899222000-100000.4991491.03810100000.7883700.09275290-1000E0.5651880.4333841000-100000.4147431.73421100000.7844000.06207650-1000F0.5259690.37

27、00151000-100000.3226592.4069110000擴散參數(shù)選取方法如下：平原地區(qū)農(nóng)村和城市遠郊區(qū)，A、B、C級穩(wěn)定度按表4直接查算，D、E、F級穩(wěn)定度則需要向不穩(wěn)定方向提半級后按表7和表8.工業(yè)區(qū)域或城區(qū)中的點源，A、B不提級，C級提到B級，D、E、F級向不穩(wěn)定方向提一級，再按表4和表5丘陵山區(qū)的農(nóng)村或城市，擴散參數(shù)選取方法同工業(yè)區(qū)。模型求解已知：現(xiàn)有境某一工廠廢氣排放煙囪高50m，主要排放物為氮氧化物。早上9點至下午3點期間的排放濃度為406.92mg/k，排放速度為1200k/h；晚上10點-凌晨4點期間的排放濃度為1160mg/:，排放速度為5700u:/h；求解：該工

28、廠方圓51公里分別在早上8點、中午12點、晚上9點空氣污染濃度分布和空氣質(zhì)量等級。根據(jù)資料，取煙囪出口徑D=1.7m、Ts=370K、Ta=293K，用Matlab模擬出風(fēng)速與煙流抬升高度間的關(guān)系。然后根據(jù)查找的境風(fēng)速方面的資料，假設(shè)取u=6.7m/s，模擬出在正下風(fēng)向0-5000m距離的擴散趨勢圖。圖6正下風(fēng)向氮氧化物濃度的擴散分布趨勢圖圖形顯示出，氮氧化物濃度約在750m處達到最大值，而在0-750m氮氧化物濃度是逐漸遞增的，說明地面反射的作用是很強的。在750m以后，氮氧化物濃度在向遠處擴散過程中逐漸衰減。最后查到了的氣象資料，得到24小時整點溫度分布圖，得到境早上8點，中午12點，晚上

29、9點的溫度分別為18攝氏度，24攝氏度，21攝氏度。9牢01總Q3巻吋點05空06戸QF戌D6點DO*1Q點li口12亡U點i*頁15F?KH空LE戸肝卓2Q點21匣圖7境氣象資料然后將這三個時間點的溫度帶入到模型中求解，可得到該工廠方圓51公里分別在早上8點、中午12點、晚上9點空氣污染濃度分布和空氣質(zhì)量等級。圖8早上8點圖9中午12點圖10晚上9點表5空氣質(zhì)量等級時間距煙囪出口正下方的距離空氣質(zhì)量等級早上8點0-3Km四級（中度污染）3-12Km三級（輕度污染）12-510Km二級（良）中午12點0-4.5Km四級（中度污染）4.5-14Km三級（輕度污染）0-3.5Km二級（良）晚上9點

30、0-3.5Km四級（中度污染）3.5-12.5Km三級（輕度污染）12.5-510Km二級（良）5.4多污染源空氣污染擴散模型若需要評價的點源多于一個，計算濃度時，將各個源對接收點濃度的貢獻進行疊加。在評價區(qū)選一原點，以平均風(fēng)的上風(fēng)方為正X軸，各個污染源對評價區(qū)任一地面店（x,y,0）的濃度總貢獻C，按下式計算：C（x.y,CJi=VQ0（14）式中，Cj是第i個點源對（x,y,0）點的濃度貢獻。5.4.1建立模型如果以汽車尾氣污染源為例，可以利用線源模式進行汽車尾氣擴散濃度分布的預(yù)測，將道路劃分成一系列面積不等的路段單元，每個路段單元等效成與風(fēng)向垂直的有限長線源，道路對承受點的貢獻濃度等效為

31、所有有限長線源對承受點濃度之和。以汽車為移動電源，當(dāng)其在固定路線上運行并達到一定車流量時，可視為固定線源。根據(jù)實際應(yīng)用，任意風(fēng)向時的線源模式采用高斯點源迭加的方法可得到較高的計算精度。1、風(fēng)向與線源垂直時的地面濃度值的計算設(shè)x軸與風(fēng)向一致，坐標(biāo)原點在線源的中點，線源在y軸上的長度為Zy,由它所產(chǎn)生的任一點（x,y,0）的濃度為(15)為誤差函數(shù)。當(dāng)y-時，從（1）式可導(dǎo)出線源排放高度為h時，無限長線源的地面濃公式Qi:(17)2、風(fēng)向與線源平行時的地面濃度值的計算取x軸與線源一致，坐標(biāo)原點和線源中點重合，并設(shè)線源長度為2，近距離(18)可作如下假設(shè)：=b（常數(shù)）3卷（19）其中a,b為常數(shù)，則

32、線源高度為h的地面濃度值的計算公式為:C（x,yr0）=.?Lerffr&erf（_-（20）若為無線長線源，因此只有上風(fēng)向的線源才對計算點的濃度有貢獻，所以Qlv2titFUCFE(r)(21)3-風(fēng)向與線源成任意交角時的地面濃度值的計算設(shè)風(fēng)向與線源間的交角為（三覚度），用簡單插方法估算兩者成任意交角時的地面濃度，其計算公式為：；=:垂直平行：心：（）（22）對無限長線源，C（垂直）按（2）式計算，C（平行）按（4）式計算。5.4.2模型中各參數(shù)的確定風(fēng)速的計算考慮汽車排氣從低位置排出時，其擴散容易受到地表面的影響，因此，線源排放高度的風(fēng)速，是使污染物擴散的風(fēng)速，按下式求出：加lm計算，為高

33、度h（m）時的實測風(fēng)速（m/s）,K為風(fēng)速高度數(shù).擴散參數(shù)的確定擴散參數(shù)可通過查表4和5得到。研究汽車排氣擴散時，必須考慮汽車行駛時大氣的機械混合作用和因汽車排熱而產(chǎn)生的大氣熱紊亂以與浮力效應(yīng)等現(xiàn)象，也必須考慮在建筑密集地因道路兩旁建筑物產(chǎn)生的大氣擾亂現(xiàn)象。當(dāng)風(fēng)向與線源垂直時，水平方向的擴散參數(shù)二對道路周圍的濃度影響不大，但當(dāng)平行風(fēng)時，二的大小對濃度有較大影響。因此，開闊的道路和狹小的道路,其初始擴散參數(shù)值：是不同的。如果考慮接近平行線源的風(fēng)向時隨車道的增大而增大。即一=二-爛（24）式中，二為汽車行駛產(chǎn)生的擾亂而引起的初始擴散參數(shù)值（m）,般設(shè)為3m；三左為把繁忙的公路視為一個線源而引起的初

34、始擴散值（m）,按下式式中；是從道路中心到行駛于第i車道的中心線的距離。綜上所述，在考慮初始擴散參數(shù)后，垂直和水平方向的擴散參數(shù)按下式計算兀=*：1一海卩（26）壬二6J廠&：（27）線源源強的計算移動線源源強通?？梢杂孟旅婀竭M行計算Qijpw=QijpLj,CjpEfpW(28)Qjpwr=社iQijpw(29)式中j車型P燃料類型Qjip,.某條線源道路第i段路上j類型車p型燃料w種污染物排放源強,g/h;Qjpw某條線源道路，j類型車p型燃料w種污染物排放源強，g/h;-：第i段路長，km.寸二：j類型車p型燃料在i段路上的車流量，輛/h.n:某條線源道路上劃分的總段數(shù)。Cip：j類型

35、車每輛每公里p型燃料消耗量，Kg/Km.E.j-:p類型燃料，w種污染物的排放因子,Kg/1.表6研究區(qū)各類機動車每輛單位里程燃燒消耗量項目輕型汽油出租車轎車輕型汽油車中型汽油車重型汽油車輕型柴油車低速貨車中型柴油車重型柴油車摩托車單位0.010.050.040.090.040.210.20.050.170.23里程8994368油耗(Kg/Km)表7研究區(qū)各類機動車HC、CO、NOx的單位里程平均排放因子項目輕型汽油出租車轎車輕型汽油車中型汽油車重型汽油車輕型柴油車與三輪中型柴油車重型柴油車摩托車HC0.350.040.3612.6711.021.331.454.190.16CO3.010.

36、603.0473.6786.622.612.675.631.252.303.452.294.669.562.392.5824.100.215.4.3模型求解已知：在2015年1月15日已經(jīng)連續(xù)三天發(fā)生重污染，假設(shè)從16日開始啟動汽車單雙號限行交通管制措施，求解市二環(huán)、四環(huán)、六環(huán)路在16日早上8點、中午12點、晚上9點時空氣污染濃度梯度變化與空氣質(zhì)量等級。我們根據(jù)得到的二環(huán)路的車流量數(shù)據(jù)來近似代替四環(huán)、六環(huán)路的車流量。表8二環(huán)路7:30-19:30時段車流量統(tǒng)計7:30-8:30234011:30-12:30196015:30-16:3016208:30-9:3098012:30-13:3010

37、8016:30-17:3016309:30-10:30120013:30-14:30143017:30-18:30260010:30-11:30141014:30-15:30187018:30-19:302400然后我們將相關(guān)數(shù)據(jù)代入到上面模型中，模擬的結(jié)果與實際檢測得到的結(jié)果比較相近，說明建立的多污染源空氣污染擴散模型具有一定的可靠性，符合實際結(jié)果。tt：R:輕履污黒：-中園珥黠：竄亙珥卓：嚴眉虧聚圖111月15日24小時空氣質(zhì)量指數(shù)30020050023El204?It：R；輕變污樂：中度厲粢：重虔活秦：嚴重污棗：寅很擲400圖121月16日24小時空氣質(zhì)量指數(shù)趨勢5.5影響空氣質(zhì)量的關(guān)鍵

38、參數(shù)以與對京津冀地區(qū)環(huán)保部門的建議報告5.5.1影響空氣質(zhì)量點哦關(guān)鍵參數(shù)通過對模型的分析，共總結(jié)出了三個影響空氣質(zhì)量的關(guān)鍵參數(shù)，源強Q,風(fēng)速u,高架點源的高度H。源強Q：污染物的源強指污染物的排放速率，污染物的濃度與源強成正比，源強越大，污染越嚴重，若要研究某個地區(qū)污染情況，清楚地掌握該地區(qū)的源強對于研究該地區(qū)污染狀況十分重要，此外，如果我們減少某個地區(qū)的污染強度，可以通過減少源強來達到減少污染的目的，比如減少汽車尾氣的排放，控制煤炭的燃燒都是通過較小源強強度來減少大氣污染。風(fēng)速u：污染物在大氣中排放的濃度與總排放量成正比，而與平均風(fēng)速成反比。若風(fēng)速增加一倍，則下風(fēng)側(cè)有害氣體濃度就減少一半。因

39、為風(fēng)力的加大，使單位時間通過煙波斷面的空氣量增大和湍流擴散增強，起著稀釋煙塵污染的作用。高架點源的高度H：高架點源的有效高度由兩部分組成，即日二二-二其中h為排放口的有效高度，二.是熱煙流的浮升力和煙氣以一定速度豎直離開排放口的沖力使煙流抬升的一個附加高度。高架點源的高度H越大，污染物傳播的圍越大。5.5.2對京津冀地區(qū)環(huán)保部門的建議報告根據(jù)本文對大氣污染的數(shù)學(xué)模型以與參考的相關(guān)數(shù)據(jù)，我們對環(huán)境保護提出以下建議：（1）壓減燃煤，調(diào)整能源結(jié)構(gòu)，推行天然氣，電力等清潔能源，降低原煤消耗所占的比例。特別是在冬季供暖季節(jié)，改造居民采取燒煤取暖的狀況，完善壓減燃煤工作管理制度，強化排放標(biāo)準，嚴格排放監(jiān)管。努力擴大天然氣、煤氣等清潔能源消費量，強化能源節(jié)約。（2）加強工業(yè)污染的防治，以循環(huán)經(jīng)濟模式發(fā)展工業(yè)經(jīng)濟。如：加強大型火電廠的脫硫、除塵以與低氮燃燒等措施。（3