上傳人：建***

認證信息

認證類型：個人認證

認證主體：曹**（實名認證）

IP屬地：上海

IP屬地：上海 上傳時間：2022-04-13 格式：DOCX 頁數(shù)：23 大?。?78.62KB 積分：20 舉報 版權(quán)申訴

1、*大學(xué)（第七屆）數(shù)學(xué)建模競賽承 諾 書本次競賽的章程與全國大學(xué)生數(shù)學(xué)建模競賽章程（以下簡稱為“競賽章程和參賽規(guī)則”，可從全國大學(xué)生數(shù)學(xué)建模競賽網(wǎng)站下載）相同，我們已認真閱讀。我們明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反參賽規(guī)則的。如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽章程和參賽規(guī)則，以保證競賽的公正、公平性。如有違反參賽規(guī)則的行為，我們的比賽成績將被取消。我們參賽選擇的題號是（從

2、A/B中選擇一項填寫）： B 參賽隊員 (打印并簽名) ：1 2. 3. （論文紙質(zhì)版與電子版中的以上信息必須一致，只是電子版中無需簽名。以上內(nèi)容請仔細核對，提交后將不再允許做任何修改。如填寫錯誤，論文可能被取消評獎資格。）為匯總參賽信息，請認真填寫以下電子表格（不是手寫），如表格空間不夠，可縮小字體：姓名學(xué)院專業(yè)學(xué)號手機郵箱張健胡靜嫻刁夢雯B題摘要本文主要分析津京冀地區(qū)主要污染源擴散（分為點源擴散模型和線源擴散模型）及其對大氣空氣質(zhì)量的影響，在現(xiàn)有的國家最新空氣污染物監(jiān)測標準（HJ633-2012環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定）的基礎(chǔ)上利用層次分析原理對現(xiàn)有空氣質(zhì)量計算模型進行改進。利用

3、改進后的評價模型對空氣質(zhì)量進行評估，并根據(jù)污染源的相關(guān)特征提出可行性建議與措施。對于問題一，聯(lián)系實際說明現(xiàn)有模型的局限性，然后引入層次分析法根據(jù)污染物的危害性對不同的污染物賦予相應(yīng)的權(quán)重得到優(yōu)化后的空氣質(zhì)量評價模型。對于問題二，通過相關(guān)網(wǎng)站查的津京冀地區(qū)的主要污染源以及污染源的主要成分，分析污染源的種類、成因、性質(zhì)、危害。對于問題三，以高斯煙羽擴散模型作為工廠污染氣體排放的基本模型，并作相應(yīng)修改使之適用于高架點源擴散，核心問題在于解決無污染氣體排放階段到有污染氣體排放階段中間過渡段污染氣體濃度的變化，最終得出了不同時刻工廠周圍污染源的濃度值。對于問題四，以高斯線源擴散為基本模型，綜合考慮二環(huán)、

4、四環(huán)、六環(huán)所處位置周圍地形地貌的差別，風(fēng)向的影響以及城區(qū)與郊區(qū)不同擴散因素，車流量情況等求出濃度分布與梯度分布，并以圖表的方式形象的展示了污染物濃度隨時間地點的梯度變化，從理論上證明了政府限車令對于緩解環(huán)境壓力的作用。在前四問的基礎(chǔ)上，我們綜合統(tǒng)籌模型的優(yōu)缺點，結(jié)合津京冀地區(qū)發(fā)展實際，提出了一系列有助于藍天計劃的詳實措施?！娟P(guān)鍵詞】AQI 高斯點源擴散 高斯線源擴散 污染物 津京冀 層次分析法1問題重述近年來，我國GDP持續(xù)快速增長，對生態(tài)平衡和自然環(huán)境造成一定的破壞，空氣污染的加重，霧霾天氣嚴重影響人們的生活和健康。一般認為影響空氣質(zhì)量的主要因素有PM2.5、PM10、二氧化氮、二氧化硫、一

5、氧化碳、臭氧、硫化氫、碳氫化合物和煙塵等，以京津冀地區(qū)為研究對象解決以下問題：（1）參考現(xiàn)有國標和美標，建立衡量空氣質(zhì)量優(yōu)劣程度等級的數(shù)學(xué)模型。（2）分析京津冀地區(qū)主要污染源及其污染參數(shù)，分析影響空氣質(zhì)量的主要污染源的性質(zhì)和種類。（3）建立單污染源空氣污染擴散模型，描述其對周圍空氣污染的動態(tài)影響規(guī)律?，F(xiàn)有河北境內(nèi)某一工廠廢氣排放煙囪高50m，主要排放物為氮氧化物。早上9點至下午3點期間的排放濃度為406.92mg/m3，排放速度為1200m3/h；晚上10點-凌晨4點期間的排放濃度為1160mg/m3，排放速度為5700m3/h；通過你的擴散模型求解該工廠方圓51公里分別在早上8點、中午12點

6、、晚上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”藍天的可行性措施和建議。2基本假設(shè)風(fēng)的平均流場穩(wěn)定，風(fēng)速均勻，風(fēng)向平直；污染物的濃度在y、z軸方向符合正態(tài)分布；污染物在輸送擴散中質(zhì)量守恒；污染源的源強均

7、勻、連續(xù)；假設(shè)地面對放射性氣體起全反射作用3符號說明I空氣質(zhì)量指數(shù)，即AQI，輸出值z垂直方向的標準差，即z方向的擴散參數(shù)C污染物濃度，輸入值u平均風(fēng)速，m/sClow小于或等于C的濃度限值，常量h是熱煙流的浮升力和煙氣以一定速度豎直離開排放口的沖力使煙流抬升的一個附加高度Chigh大于或等于C的濃度限值，常量q源強，即單位時間內(nèi)排放的污染物，g/sIlow對應(yīng)于Clow的指數(shù)限值，常量Ihigh對應(yīng)于Chigh的指數(shù)限值，常量y分別為水平方向的標準差，即y方向的擴散參數(shù)4問題分析問題一的分析我國現(xiàn)有空氣質(zhì)量計算模型主要依據(jù)HJ633-2012環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定里一維插值分析法

8、。一維差值公式，可計算出每一種污染物對應(yīng)的AQI:注：中美兩國濃度限值見附錄一我國的空氣質(zhì)量評價表與美國的空氣質(zhì)量評價表有一點區(qū)別，主要在于評價空氣質(zhì)量標準中各污染物參考上限值不一樣，美國的空氣質(zhì)量優(yōu)良的標準偏高，中國偏低，這是由我國的基本國情決定的。針對津京冀地區(qū)而言，我們應(yīng)以國標為重，美標為輔建立空氣質(zhì)量評價標準。根據(jù)美標和國標AQI算法，觀測點取各種污染物所對應(yīng)AQI的最大值來評價整體的空氣質(zhì)量。公式如下：AQI=maxIAQI1, IAQI2, IAQI3, IAQI4. IAQIn 這就出現(xiàn)了一個問題：以津京冀地區(qū)北京某一時刻的各主要污染氣體的IAQ為例（表3），其中主要污染氣體PM

9、2.5、PM10在不同的時刻均保持較高值，且O3的AQI分指數(shù)也相對高，如果僅僅按照最大值來評估，該最大值只是其中的一種污染物，其他污染物對空氣造成的影響沒有考慮在內(nèi)。試想想，如果有一個地區(qū)各項AQI分指標為（75，1，1，2，3，10），那么最終公布的AQI=75。而另一個地區(qū)的各項AQI分指數(shù)為（75，75，75，75，75，75）那么按國標算法公布的AQI也應(yīng)為75，很顯然，這兩個地區(qū)的環(huán)境質(zhì)量差距是可觀的，按照現(xiàn)有算法很可能會給群眾帶來誤解，認為這兩個地區(qū)環(huán)境質(zhì)量相當(dāng)。為了解決這個問題，我們利用了層次分析法，賦予每一個污染項一定的權(quán)值，最后以加權(quán)的形式給出最終的AQI值。表1：北京某一

10、時刻各主要污染氣體的AQIi污染物分指數(shù)（AQI）7:009:0012:0015:0020:00PM2.57328222275SO201101NO2766614O32527303340CO33336PM1015172504573問題二的分析通過網(wǎng)絡(luò)、圖書館等途徑查找我們所需要的具體數(shù)據(jù)，針對京津冀地區(qū)的地理位置和經(jīng)濟狀況分析其主要污染源，找到具體有力的相關(guān)參數(shù)，分析出影響空氣質(zhì)量的污染源的性質(zhì)及其種類。同時也為后面問題的解答提供基礎(chǔ)知識，從而更好地驗證第三第四小題解答的正確性。問題三的分析對于問題3，為了處理單污染源擴散模型，我們建立了高斯煙羽擴散模型，考慮到題目中要解決的是工廠煙囪排放廢氣問

11、題，有效源高（煙囪的架高與煙云抬升高度之和）將在問題中起到一定影響，于是我們將高斯模型進行修改，綜合考慮有效源高和地面反射的作用，最終確定了高架點源高斯擴散模型來對工廠污染物持續(xù)釋放時間內(nèi)周圍污染物濃度進行求解。至于題目中要求對不同時間點濃度值的計算，此時持續(xù)排放的時刻與未進行排放時刻，工廠周圍污染物濃度必然會有所不同，另外，不同的時刻大氣穩(wěn)定度也不一樣，這會導(dǎo)致擴散參數(shù)的不同，為了區(qū)分這些不同，建立了高斯擴散下排放源強為0mg/s時的污染物擴散模型。問題四的分析對于建立多污染源空氣污染擴散模型，我們只需在問題三的基礎(chǔ)上，將高斯點源擴散模型轉(zhuǎn)化為高斯線源擴散模型即可，線源模型下的污染物濃度C即

12、為高斯點源模型下污染物濃度沿線源方向的定積分或廣義積分。為了解決題目中政策實施后，不同地點，不同時間污染物濃度C的變化，對于不同環(huán)路早上8點，中午12點，晚上9點這三個特殊時間點我們綜合考慮了氣象（大氣穩(wěn)定度）和社會方面（上下班高峰期等）的因素，從而使算出的數(shù)據(jù)更具實際意義。最后，通過對空氣質(zhì)量等級的評估，來闡述政府調(diào)控在保護環(huán)境當(dāng)中的重要作用。問題五的分析在問題一中我們通過層次分析法建立了新的空氣質(zhì)量評價標準，根據(jù)表3中的數(shù)據(jù)，我們?nèi)庋劭梢灾庇^的看出PM2.5、PM10、O3等污染氣體濃度較高。問題二中對津京冀地區(qū)分別進行分析得出：表2對津京冀地區(qū)主要污染源的統(tǒng)計地區(qū)主要污染源北京機動車尾氣

13、排放天津工業(yè)廢氣污染河北燃煤、工業(yè)廢氣污染津京冀地區(qū)的主要污染原因不盡相同，問題三、四針對北京和河北舉出典型的兩個例子分析濃度分布變化，我們可以直觀的看出不同時間段濃度變化的程度并可以在本題中提出針對各個時間段的改善措施。模型的建立與求解問題一的建模和求解1 國標與美標濃度限度的比較見附錄一2 在研究各種指標集對某物影響的過程中，不僅指標集中的最大值具有最重要的作用，次大值等的作用也不容忽視，甚至具有與最大值類似的影響。因此，為了在評價體系中引入次大值等的影響。我們采用層次分析加權(quán)法得到各個污染物IAQI所占的權(quán)重。ai即為各個污染物IAQI的權(quán)值。具體方法如下：1構(gòu)造判斷矩陣 通過對指標之間

14、兩兩重要程度進行比較和分析判斷，構(gòu)造判斷矩陣。層次分析法在對指標的相對重要程度進行測量時，我們引入了九分位的相對重要的比例標度。令A(yù)為判斷矩陣，用以表示同一層次各個指標的相對重要性的判斷值。根據(jù)心理學(xué)家提出的“人區(qū)分信息等級能力為7+2”的研究結(jié)論，有如下評分規(guī)則：表3權(quán)重的評分規(guī)則甲指標與乙指標比較極端重要強烈重要明顯重要比較重要重要較不重要不重要很不重要極不重要甲指標評分9753113151719注：取8,6,4,2,1/2,1/4,1/6,1/8為上述評價值的中間值。在綜合考慮了京津冀地區(qū)的污染情況后，我們建立了比較矩陣A如下：表4空氣質(zhì)量指數(shù)空氣評價判斷矩陣AO3SO2NO2COPM1

15、0PM2.5O314331/21/2SO21/413/43/41/81/8NO21/34/3111/61/6CO1/34/3111/61/6PM10286611PM2.52866113對判斷矩陣進行一致性檢驗：3.1用matlab求出判斷矩陣A的特征向量和最大特征值，最大特征值對應(yīng)的特征向量歸一化后即為各分指標IAQI的權(quán)重。3.2計算判斷矩陣的一致性指標。具體公式為：CI=(-n)/( -1)CI為0時，A一致；CI越大，A的不一致程度就越高。表5隨機一致性指標RI：N1234567891011RI000.580.901.121.241.321.411.451.491.51CR=CI/RI,

16、當(dāng)CR<0.1時，A的不一致性程度在允許范圍內(nèi)，此時可用A的特征向量作為權(quán)向量。表6各指標權(quán)重和一次性指標值O3SO2NO2COPM10PM2.5CICR0.16900.04230.05630.05630.33800.3380-3.5527e-16-2.8651e-166然后把各污染物的IAQI與其相應(yīng)的權(quán)重相乘得到綜合指數(shù)AQI。注：中國國家AQI評價標準見附錄二表7各污染物AQI分指標SO2NO2O3COPM2.5PM10地區(qū)一5830185073地區(qū)二32523073地區(qū)一與地區(qū)二按國標算法兩者AQI都為，但從表中顯然可以看出地區(qū)二的環(huán)境優(yōu)于地區(qū)一，因此國標算法在這個問題上存在缺陷

17、。而依據(jù)我們改進的算法，地區(qū)一和地區(qū)二的AQI分別為48.32和36.00。改進后的算法很好的辨別出了兩地的環(huán)境差異。當(dāng)然基于此算法環(huán)境優(yōu)良等級也應(yīng)做相應(yīng)調(diào)整。問題二的建模和求解：有關(guān)數(shù)據(jù)顯示，京津冀地區(qū)主要污染源分別為機動車、燃煤、揚塵、流動源以及工業(yè)生產(chǎn)，其中燃煤和機動車占主要因素，見以下表格及餅圖 。近幾年，北京主要污染源為機動車，見下圖及表格。2012年，京津冀機動車氮氧化物排放量68.2萬噸，占氮氧化物排放總量的30%，其中北京機動車氮氧化物排放量占本地區(qū)氮氧化物的比重達45%，分別高于天津28.8個和河北13.9個百分點。天津工業(yè)污染影響最大,數(shù)據(jù)顯示，2012年，京津冀工業(yè)二氧化

18、硫排放量占二氧化硫排放總量的91.2%；工業(yè)氮氧化物排放量占氮氧化物排放總量的68.4%；工業(yè)煙（粉）塵排放量占煙（粉）塵排放總量的82.6%。分省市看，天津工業(yè)污染影響最大。河北燃煤污染尤其嚴重。河北二氧化硫排放占京津冀八成,數(shù)據(jù)顯示，2012年，京津冀燃煤消費總量38927萬噸，河北煤炭消費量占其能源消費總量的88.8%，遠遠高于北京的25.4%和天津的59.6%。煤炭消費排放出大量二氧化硫，對大氣環(huán)境造成很大影響，2012年河北二氧化硫排放量占京津冀的80.8%。表8地點主要污染源所含主要污染成分危害北京汽車尾氣固體懸浮微粒、一氧化碳、碳氫化合物、氮氧化合物、鉛及硫氧化合物等造成霧霾，

19、顆粒物對人體危害大天津工業(yè)生產(chǎn)SO2 、粉塵、NOx對工農(nóng)業(yè)建設(shè)和人民健康危害極大河北燃煤SO2、CO、粉塵、NOx酸雨， 破壞生態(tài)環(huán)境表9京津冀區(qū)域2.5來源所占比例燃煤34%機動車16%工業(yè)、外來輸送、揚塵、餐飲和其他50%34%16%50%表10北京霧霾天氣的成因汽車尾氣燃煤污染建筑揚塵工業(yè)排放22.20%16.70%16.30%15.70%圖1京津冀區(qū)域2.5來源所占比例 圖2北京霧霾成因由以上查找數(shù)據(jù)可以總結(jié)：主要污染源機動車尾氣的所含污染成分為CO和氮氧化物，燃煤所含污染物為二氧化硫和一氧化碳，工業(yè)生產(chǎn)含有的污染成分為二氧化硫和氮氧化物及粉塵。問題三的建模和求解：對于單污染源的擴散

20、，我們以高斯煙羽擴散為基本模型：本題為高架點源的高斯擴散模式。有效源位于坐標原點o處，平均風(fēng)向與x軸平行，并與x軸正向同向。假設(shè)點源在沒有任何障礙物的自由空間擴散，不考慮下墊面的存在。大氣中的擴散是具有y與z兩個坐標方向的二維正態(tài)分布，當(dāng)兩坐標方向的隨機變量獨立時，分布密度為每個坐標方向的一維正態(tài)分布密度函數(shù)的乘積。由正態(tài)分布的假設(shè)條件，參照正態(tài)分布函數(shù)的基本形式，取0，則在點源下風(fēng)向任一點的濃度分布函數(shù)為： （3-1）式中 C空間點（x，y，z）的污染物的濃度，mg/m3； A（x）待定函數(shù)； y、z分別為水平、垂直方向的標準差，即y、x方向的擴散參數(shù)，m。 由守恒和連續(xù)假設(shè)條件和，在任一垂

21、直于x軸的煙流截面上有： （3-2）式中 q源強，即單位時間內(nèi)排放的污染物，g/s； u平均風(fēng)速，m/s。 將式（3-1）代入式（3-2）, 由風(fēng)速穩(wěn)定假設(shè)條件，A與y、z無關(guān)，考慮到和，積分可得待定函數(shù)A（x）： （3-3）將式（3-3）代入式（3-1），得大空間連續(xù)點源的高斯擴散模式 （3-4）式中，擴散系數(shù)y、z與大氣穩(wěn)定度和水平距離x有關(guān)，并隨x的增大而增加。為了解決題中煙囪排放污染物問題，我們對高斯點擴散模型做了一些修改，得到適用于高架源的高斯高架源擴散模型見下圖3。點源在地面上的投影點o作為坐標原點，有效源位于z軸上某點， zH。高架有效源的高度由兩部分組成，即Hhh，其中h為排放

22、口的有效高度，h是熱煙流的浮升力和煙氣以一定速度豎直離開排放口的沖力使煙流抬升的一個附加高度，如圖3所示。當(dāng)污染物到達地面后被全部反射時，可以按照全反射 圖3 地面全反射高架點源模型原理，用“像源法”來求解空間某點k的濃度。圖3中k點的濃度顯然比大空間點源擴散公式（3-4）計算值大，它是位于(0，0，H)的實源在k點擴散的濃度和反射回來的濃度的疊加。反射濃度可視為由一與實源對稱的位于(0，0，H)的像源（假想源）擴散到k點的濃度。由圖可見，k點在以實源為原點的坐標系中的垂直坐標為(z-H)，則實源在k點擴散的濃度為式（3-4）的坐標沿z軸向下平移距離H： （3-5）k點在以像源為原點的坐標系中

23、的垂直坐標為(zH)，則像源在k點擴散的濃度為式（3-4）的坐標沿z軸向上平移距離H： （3-6）由此，實源Cs與像源Cx之和即為k點的實際污染物濃度：（3-7） 若污染物到達地面后被完全吸收，則Cx0，污染物濃度C（x，y，z，H）Cs，即式（3-5）。 （2）地面全部反射時的地面濃度。實際中，高架點源擴散問題中最關(guān)心的是地面濃度的分布狀況，尤其是地面最大濃度值和它離源頭的距離。在式（3-7）中，令z0，可得高架點源的地面濃度公式： （3-8） 上式中進一步令y0則可得到沿x軸線上的濃度分布： （3-9）考慮到題目中對于不同時間點污染物濃度的計算，我們對高斯擴散模型進行了一些修改，得到了與時

24、刻t有關(guān)的高斯擴散模型（此模型更加適用于源強為0時的高斯擴散）:（3-10）題中問題求解：1. 參數(shù)的確定。高斯擴散公式的應(yīng)用效果依賴于公式中的各個參數(shù)的準確程度，尤其是擴散參數(shù)y、z及煙流抬升高度h的估算。其中，平均風(fēng)速u取多年觀測的常規(guī)氣象數(shù)據(jù)；源強q可以計算或測定，而y、z及h與氣象條件和地面狀況密切相關(guān)。h的確定：h可用如下公式近似計算：在查閱了河北省歷年平均風(fēng)速以及污染物排放指數(shù)后，大致計算出h=12.8m則H=h+h=50+12.8=62.8m擴散參數(shù)y、z的確定：y、z的確定依賴于當(dāng)?shù)氐娘L(fēng)速，太陽照射角，云量，對于它們的確定一般都是在確定了穩(wěn)定度（依據(jù)大氣穩(wěn)定度的等級表）的情況下

25、，根據(jù)擴散參數(shù)函數(shù)表（添加表名稱）來確定y、z關(guān)于X的具體表達式。表11為Pasquill穩(wěn)定度分級方法不論何種天空狀況，日落前一個小時和日出后一個小時穩(wěn)定度為中性，且仲夏晴天中午為強日照，寒冬晴天中午為弱日照。表12橫向擴散參數(shù) 冪函數(shù)表達式數(shù)據(jù)表13 Briggs擴散參數(shù)根據(jù)河北省氣象站1971-2009年統(tǒng)計數(shù)據(jù)，我們大致確定出早上8點大氣穩(wěn)定度等級為D，中午12點大氣穩(wěn)定度等級為B，晚上9點大氣穩(wěn)定度等級為F，則不同時刻不同地點y、z可以相繼算出。確定參數(shù)后，我們用matlab對題中問題進行了簡答并對不同時刻工廠周圍的污染物濃度作出圖形，以便形象的觀測。 圖4早上8點污染物濃度分布 圖

26、5中午12點污染物濃度分布 圖6晚上9點污染物濃度分布從圖中我們大致可以得到結(jié)論：污染物濃度整體變化趨勢是隨著距工廠的距離的增加遞減的，另外，從中午點與早上點同一地點的污染物濃度比較來看，早上點污染物的濃度比起中午點要低，這說明正在排污時污染物對周圍環(huán)境的破壞是最大的，因此從保護環(huán)境的角度來看，限制污染源的排放是一個很有效的措施。對于3個點質(zhì)量等級的判斷，我們不妨取3個時刻污染物濃度最大的地點進行計算，因為污染物濃度最大的地點是最容易形成污染物聚集而對環(huán)境造成破壞的。根據(jù)matlab所計算數(shù)據(jù)，IAQI計算公式以及第一問建立的AQI模型可以對題中3個時間點進行環(huán)境質(zhì)量等級評估。早上點濃度最大值

27、為：0.0082/m3其對周圍氮氧化物IAQI的貢獻值為4.1；中午點濃度最大值為：0.035mg/m3其對周圍氮氧化物IAQI的貢獻值為17.5；晚上點濃度最大值為：0.0028mg/m3其對周圍氮氧化物IAQI的貢獻值為1.4。根據(jù)河北環(huán)境監(jiān)測局網(wǎng)站數(shù)據(jù)河北市郊區(qū)各個污染源平均IAQI分布如下：表14 河北省郊區(qū)各污染物AQI分指標PM2.5SO2NO2O3COPM107317253072根據(jù)matlab計算后可知污染物排放后各個污染源AQI分布如下：表15河北省郊區(qū)各污染物AQI分指標PM2.5SO2NO2O3COPM108:0073111.125307212:0073124.52530

28、7221:007318.4253072表16 根據(jù)題一所建模型，各時刻AQI值如下AQI排污前54.96638:0055.22312:0055.977521:0055.071注：對應(yīng)代碼見附錄三問題四的建模求解在問題三的基礎(chǔ)上，我們將高斯點源模型擴展為高斯線源模型，并對不同情況下線源模型的適用進行分類說明（為了簡化計算，下列模型均建立在風(fēng)向與線源垂直的基礎(chǔ)上）。1.有風(fēng)時(u>0.4m/s)采用高斯煙流模式設(shè)x軸與風(fēng)向一致，坐標原點在線源的中點，線源在y軸上的長度為2y0。由它所產(chǎn)生的在任一點（x,y,0）的濃度為：C(x,y,0；h)= *exp(-)* (4-1)式中；u為平均風(fēng)速（

29、m/s）；y、z為水平和垂直方向的擴散參數(shù)；H為有效源高；C為污染物濃度；q為線源強度。當(dāng)為無限長線源時，地面濃度公式為： (4-2)2. 準靜風(fēng)時(u < 0.4 m / s ) 采用高斯煙團模式通常平均風(fēng)速小于1m/s時，適用于連續(xù)點源的準定常煙流模式不再適用。此時用高斯煙團模式代替煙流模式。取y軸與線源一致，坐標遠點與線源重點重合時，對于線源排放高度為h的地面濃度計算公式為： (4-3)當(dāng)為無限長線源時，準靜止風(fēng)條件下的地面濃度計算公式為： (4-4)參數(shù)的確定：1. 有風(fēng)時，參數(shù)y、z的確定與問題三中確定方法相同，在這不作重復(fù)說明。2.準靜風(fēng)時（u<0.4m/s）擴散參數(shù)的

30、確定可參照Turner分類法進行的Pasquill穩(wěn)定度分類表。2. 線源源強的計算q=車輛平均排放因子*車流量排放因子：指車輛行駛單位距離平均排出污染物的量車流量：每小時通過的車輛數(shù)目實例解答：在查閱了北京歷年氣象資料后，發(fā)現(xiàn)北京的情況適合用有風(fēng)時的高斯線源模型來求解。通過查閱相關(guān)資料，我們得出了北京二，四，六環(huán)在題中3個時間點的車流量如下表所示：表17北京政策實施前車流量 單位（輛/小時） 環(huán)數(shù) 時間8:0012:0021:00二環(huán)1200060008600四環(huán)600030004000六環(huán)300020002000表18北京政策實施后車流量 單位（輛/小時） 環(huán)數(shù) 時間8:0012:0021

31、:00二環(huán)800040005600四環(huán)400020003000六環(huán)200015001500表19北京車輛平均排放因子為：汽車尾氣中不同污染物排放因子 單位（毫克/米.輛）污染物CONOx排放因子40.052.03二環(huán)：綜合考慮地理位置和時間的影響，二環(huán)早上8點大氣穩(wěn)定度為C，中午12點為B，晚上9點為E根據(jù)不同時刻大氣穩(wěn)定度可確定出不同地點y、z。用matlab將各數(shù)據(jù)帶入高斯線源擴散模型即可求出各個時刻的汽車尾氣濃度變化。圖7二環(huán)CO濃度曲線圖8二環(huán)NOx濃度曲線四環(huán)：綜合考慮地理位置和時間的影響，四環(huán)早上8點大氣穩(wěn)定度為C，中午12點為B，晚上9點為D根據(jù)不同時刻大氣穩(wěn)定度可確定出不同地點

32、y、z。用matlab將各數(shù)據(jù)帶入高斯線源擴散模型即可求出各個時刻的汽車尾氣污染物濃度分布情況。 圖9四環(huán)CO濃度曲線 圖10四環(huán)NOx濃度曲線六環(huán)：綜合考慮地理位置和時間的影響，六環(huán)平均風(fēng)速2.1m/s，六環(huán)早上8點大氣穩(wěn)定度為D，中午12點為C，晚上9點為D根據(jù)不同時刻大氣穩(wěn)定度可確定出不同地點y、z。用matlab將各數(shù)據(jù)帶入高斯線源擴散模型即可求出各個時刻的汽車尾氣污染物濃度分布情況。 圖11六環(huán)CO濃度曲線 圖12六環(huán)NOx濃度曲線分析：經(jīng)過上面幾幅圖的對比，我們發(fā)現(xiàn)，就同一地點而言，比如二環(huán)，早晨8點時刻比其它兩個時刻CO和NOx的濃度都要高，這正符合實際情況，因為早晨8點處于上班

33、高峰期，車流量的增多必然會引起污染物濃度相較其他時間偏高。而從不同地點即不同環(huán)度來比較，同一時刻污染物濃度也不盡相同，這也正說明了氣象條件（風(fēng)速）對污染物擴散的影響。但從整體來看，車流量是影響污染物擴散的主要因素，因此北京單雙號限行措施會在很大程度上緩解北京的空氣環(huán)境壓力。題中問題解答：污染物在同一地點不同時間的梯度變化可以上面圖形中看出。至于空氣質(zhì)量等級的評估，我們更愿意去考慮線源對于周圍環(huán)境影響的大小。故列表如下：表20距污染源距離X政策實施前CO濃度mg/s單項API值政策實施后CO濃度mg/s單項API值1016.28112.5610.85101.7503.5435.42.3523.5

34、1001.7717.71.1811.82000.898.90.595.9從圖表中可以看出，車流量的減少在一定程度上釋放了周圍空氣環(huán)境的壓力，同一時間點，政策實施前后AQI值有所改變，這在一定程度說明了限行政策的可行性。另外從空間來看，當(dāng)離線源距離超過200時，線源對周圍環(huán)境的影響便不再明顯，而影響較大的區(qū)域為50M以內(nèi)，這也就指導(dǎo)人們房屋建筑應(yīng)與道路拉開一定距離，不能靠的太近。盡管車流數(shù)據(jù)不大準確，但我們還是可以看出，政府政策調(diào)控對于環(huán)境保護的有力作用。注：對應(yīng)代碼見附錄三問題五的求解由第二問的分析以及第三問的高斯煙羽擴散模型和第四問的高斯線源模型所獲得結(jié)果和數(shù)據(jù)可以直觀的對比出影響津京冀地區(qū)

35、空氣質(zhì)量的主要參數(shù)。表21津京冀地區(qū)主要污染源成分參數(shù)分析地區(qū)主要污染源主要成分北京機動車尾氣排放PM2.5天津工業(yè)廢氣污染SO2 、粉塵、NOx 河北燃煤、工業(yè)廢氣污染SO2、CO、粉塵、NOx對津京冀地區(qū)環(huán)保部分的建議：地區(qū)建議北京根據(jù)第四問可知在晚上九點鐘的時候解除限行指令；推廣使用天然氣；改用優(yōu)質(zhì)低硫低灰煤天津發(fā)電廠燃煤機組的清潔化利用、開發(fā)新型能源、嚴格執(zhí)行相關(guān)法律法規(guī)；定期清洗馬路控制粉塵河北大力開發(fā)使用新型能源，例如風(fēng)能、太陽能、利用燃煤企業(yè)余熱發(fā)電等，嚴格監(jiān)督工廠做到合格排污，根據(jù)第三問分析可知，中午要注重調(diào)控工廠廢氣的釋放對于APEC藍的可行性措施和建議：車輛限行，提高車輛質(zhì)

36、量，加大對電力汽車的研發(fā)和使用力度。在華北地區(qū)建立起聯(lián)動控制工期質(zhì)量的工作機制，加大環(huán)境執(zhí)法的執(zhí)行力度。模型的優(yōu)缺點及推廣優(yōu)點：在建模過程中，我們針對不同情況，對模型進行了相應(yīng)修改，達到了基礎(chǔ)模型的通用性和特定模型的針對性的結(jié)合，并且由模型求出的解在一定程度上與實際吻合較好，對實際生活有一定的指導(dǎo)意義。缺點：以上模型大多為了求解方便都是建立在理想情況下的，未能細節(jié)深入到實際細小的方面，比如污染源周圍建筑物，污染物的沉降效應(yīng)等都未作考慮，這就使得模型在一些特定情況下不適用。推廣：以上模型還可以推廣到重金屬的污染擴散求解等問題，另若將模型在做一維擴展，可以用來求解面源擴散問題。參考文獻1 

37、;姜啟源 數(shù)學(xué)模型（第三版），高等教育出版社，2003年8月2謝紹東等 汽車排放污染物濃度的預(yù)測，中國環(huán)境科學(xué)，1995年10月3 卓金武，魏永生，秦建，李必文Matlab在數(shù)學(xué)建模中的應(yīng)用 M 北京：北京航空航天大學(xué)出版社，20114 張斌才，趙軍大氣污染擴散的高斯煙羽模型及其GIS集成研究J 環(huán)境監(jiān)測管理與技術(shù)，2008，20（5）：17-195 環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定（試行）HJ633-2012 Z.北京：中國環(huán)境科學(xué)出版社 6 環(huán)境空氣質(zhì)量標準GB3095-2012 Z.北京：中國環(huán)境科學(xué)出版社附錄附錄一附錄二中國空氣質(zhì)量等級分類空氣質(zhì)量指數(shù)值中國級別中國分類0-

38、50一級優(yōu)51-100二級良101-150三級輕度污染151-200四級中度污染201-250五級重度污染251-300301-350六級嚴重污染351-400401-450451-500500+附錄三問題一對應(yīng)MATLAB代碼%利用層次分析加權(quán)法計算各污染物的權(quán)重A=1,6,2,5,1/2,1/21/6,1,1/3,5/6,1/12,1/121/2,3,1,2/5,1/4,1/41/5,6/5,2/5,1,1/10,1/102,12,4,10,1,12,12,4,10,1,1;D,X=eig(A);ZD=max(max(X);m,m=size(D);p=max(X);for i=1:mif

39、p(i)=ZDq=iendendz=sum(D);for i=1:mfor j=1:mD(i,j)=D(i,j)/z(j);endendB=(D(:,q)'RI=0,0,0.58,0.9,1.12,1.24,1.32,1.41,1.45,1.49,1.51;CI=(ZD-m)/(m-1)CR=CI/RI(m)Q=zeros(11,1);for i=1:11Q(i,1)=0;for j=1:6Q(i,1)=Q(i,1)+sum(i,j)*B(1,j);endEnd問題三對應(yīng)matlab碼%中午12點圖形x=0:100:20000;oy=0.433384.*(x).0.565188;oz=0.0570251.*(x).1.09356;H=60.8;u=2.5;q=135.64;C=q./(2*pi*u*oy.*oz).*exp(-H2./(oz.2);plot(x,C,'bd',x,C,'b-')%晚上9點圖形x=0:100:20000;oy=0.433384.*(x).0.565188;oz=0.0570251.*(x).1.09356;H=60.8;u=2.5;q=1836.7;C_=q*3600*6./(4*pi*1