污染物擴散模型-深圳數(shù)學建模.doc

賽區(qū)評閱編號（由賽區(qū)組委會填寫）：2015高教社杯全國大學生數(shù)學建模競賽承 諾 書我們仔細閱讀了全國大學生數(shù)學建模競賽章程和全國大學生數(shù)學建模競賽參賽規(guī)則（以下簡稱為“競賽章程和參賽規(guī)則”，可從全國大學生數(shù)學建模競賽網(wǎng)站下載）。我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反競賽章程和參賽規(guī)則的，如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽章程和參賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽章程和參賽規(guī)則的行為，我們將受到嚴肅處理。我們授權(quán)全國大學生數(shù)學建模競賽組委會，可將我們的論文以任何形式進行公開展示（包括進行網(wǎng)上公示，在書籍、期刊和其他媒體進行正式或非正式發(fā)表等）。我們參賽選擇的題號（從A/B/C/D中選擇一項填寫）： C 我們的報名參賽隊號（12位數(shù)字全國統(tǒng)一編號）： 參賽學校（完整的學校全稱，不含院系名）： 溫州醫(yī)科大學 參賽隊員 (打印并簽名) ：1. 章成俊 2. 楊超 3. 謝錦 指導教師或指導教師組負責人 (打印并簽名)： 日期： 年 月 日賽區(qū)評閱編號（由賽區(qū)組委會填寫）：2015高教社杯全國大學生數(shù)學建模競賽編 號 專 用 頁賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）：評閱人備注送全國評閱統(tǒng)一編號（由賽區(qū)組委會填寫）：全國評閱隨機編號（由全國組委會填寫）：對垃圾處理廠污染的動態(tài)監(jiān)控及居民補償摘 要城市垃圾處理問題是一個世界性難題。目前垃圾焚燒正逐步成為中國垃圾處理的主要手段之一。本論文構(gòu)根據(jù)題目設(shè)置的垃圾處理廠規(guī)模，建立了環(huán)境動態(tài)監(jiān)控體系，并根據(jù)潛在污染風險對周圍居民進行了合理經(jīng)濟補償?shù)脑O(shè)計。對于問題（1），為了實現(xiàn)對垃圾焚燒廠煙氣排放及相關(guān)環(huán)境影響狀況的動態(tài)監(jiān)控，本論文在高斯煙羽模型的基礎(chǔ)上進行改進，引入溫度、降雨對污染物擴散的影響，建立了新的污染物擴散模型。本論文創(chuàng)新性的提出了風雨影響指數(shù),用來衡量風向、降雨對顆粒物擴散的影響。本論文將抽象的污染物含量形象化，利用空氣污染指數(shù)API描述具體的污染程度及其給周圍居民帶來的影響。并且從不同角度給出了模型檢驗，驗證了所建模型的準確性。對于問題（1）具體賠償方案的制定，在綜合考慮了不同方位風向頻率、受污染時間、受污染程度的基礎(chǔ)上，本論文使用了層次分析法，并且進行了一致性檢驗，使得賠償方案具有說服力。通過MATLAB編程，計算出當政府和垃圾處理廠共支付風險賠償金為時，得出居住地的每位居民應(yīng)得的賠償金額計算公式。對于監(jiān)測點的設(shè)置，經(jīng)計算共需21個，具體布置情況見后文。對于問題（2），在題目所述的發(fā)生事故的情況下，對污染物的具體含量進行了合理的預(yù)測與假設(shè)。模擬出酸性物質(zhì)與顆粒物的影響范圍，并根據(jù)具體的污染程度設(shè)置不同的污染區(qū)。對每個污染區(qū)的不同情況設(shè)置更改監(jiān)測點的設(shè)置，并且在問題（1）的基礎(chǔ)上對居民的經(jīng)濟補償進行合理修改。關(guān)鍵詞：高斯煙羽模型，層次分析法，空氣污染指數(shù)，煙氣抬升公式一、 問題重述“垃圾圍城”是世界性難題，在今天的中國顯得尤為突出。數(shù)據(jù)顯示，目前全國三分之二以上的城市面臨“垃圾圍城”問題，垃圾堆放累計侵占土地75萬畝。因此，垃圾焚燒正逐步成為中國垃圾處理的主要手段之一。然而，由于政府監(jiān)管不力、投資者目光短淺等多方面的原因，致使前些年各地建設(shè)的垃圾焚燒電廠在運營中出現(xiàn)了環(huán)境污染問題，給垃圾焚燒技術(shù)在我國的推廣造成了很大阻力，許多城市的新建垃圾焚燒廠選址都出現(xiàn)因居民反對而難以落地的局面。在垃圾焚燒廠運行監(jiān)管方面，目前主要是在垃圾焚燒廠內(nèi)進行測量監(jiān)控，缺少從周邊環(huán)境視角出發(fā)的外圍動態(tài)監(jiān)控，因而難以形成為民眾所信服的全方位垃圾焚燒廠環(huán)境監(jiān)控體系。深圳市某地點計劃建立一個中型的垃圾焚燒廠，計劃處理垃圾量1950噸/天（設(shè)置三臺可處理垃圾650噸/天的焚燒爐，排煙口高度80米，每天24小時運轉(zhuǎn)）。從構(gòu)建環(huán)境動態(tài)監(jiān)控體系、并根據(jù)潛在污染風險對周圍居民進行合理經(jīng)濟補償?shù)男枨蟪霭l(fā)，有關(guān)部門希望能綜合考慮垃圾焚燒廠對周圍帶來環(huán)境污染以及其他危害的多種因素（例如，焚燒爐的污染物排放量、居住點離開垃圾焚燒廠的距離、風力和風向及降雨等氣象條件、地形地貌以及建筑物的遮擋程度等等），在進行科學定量分析的基礎(chǔ)上，確立一套可行的垃圾焚燒廠環(huán)境影響動態(tài)監(jiān)控評估方法，并針對潛在環(huán)境風險制定出合理的經(jīng)濟補償方案。請你在收集相關(guān)資料的基礎(chǔ)上考慮以下問題：(1) 假定焚燒爐的排放符合國家新的污染物排放標準（參見附件1），根據(jù)垃圾焚燒廠周邊環(huán)境設(shè)計一種環(huán)境指標監(jiān)測方法，實現(xiàn)對垃圾焚燒廠煙氣排放及相關(guān)環(huán)境影響狀況的動態(tài)監(jiān)控。以你設(shè)計的環(huán)境動態(tài)監(jiān)控體系實際監(jiān)控結(jié)果為依據(jù)，設(shè)計合理的周圍居民風險承擔經(jīng)濟補償方案。(2) 由于各種因素焚燒爐的除塵裝置（如袋式除塵器）損壞或出現(xiàn)其他故障導致污染物的排放增加，致使相關(guān)各項指標將嚴重超標（如：煙塵濃度、二氧化硫、氮氧化物、一氧化碳、二惡英類及重金屬等排放超標，附件2給出了一臺可處理垃圾350噸/天的焚燒爐正常運作時的在線排放監(jiān)測記錄）。請在考慮故障發(fā)生概率的情況下修正你設(shè)計的監(jiān)測方法和補償方案。二、模型假設(shè)（1）.假設(shè)污染物在所建的坐標軸的、風向上分布為正態(tài)分布；（2）.假設(shè)全部高度風速均勻穩(wěn)定；（3）.假設(shè)污染物在擴散中不會相互轉(zhuǎn)化；（4）.假設(shè)源強是連續(xù)均勻穩(wěn)定的；（5）.假設(shè)附件4所給風速距地面10米處所測得風速；三、符號約定任意點的污染物濃度，單位：mg/ m3；源強，單位時間內(nèi)污染物排放量，單位：mg/s；側(cè)向擴散系數(shù)；豎向擴散系數(shù)；排放口的平均風速，單位：m/s；煙囪的有效高度；污染源排放點至下風向上任一點的距離，單位：m；污染物的中心軸在直角水平方向上到任一點的距離，單位：m；從地表到任一點的高度，單位：m；煙囪的幾何高度，單位：m；污染物的抬升高度，單位：m；大氣壓力，單位：；煙氣熱釋放率，單位：KJ/s；實際排煙率，單位： m3/s；排氣筒出口處污染物排出速度，單位：m/s； 煙氣出口溫度，單位：K；煙氣出口溫度與環(huán)境溫度差，單位：K；。：表示第月的風向為（代表風向分別從正北方向沿順時針依次到到西北方向）風雨影響指數(shù)；四、模型建立與求解4.1問題一：4.1.1問題一分析根據(jù)垃圾焚燒廠周邊環(huán)境設(shè)計一種環(huán)境指標監(jiān)測方法，以實現(xiàn)對垃圾焚燒廠煙氣排放及相關(guān)環(huán)境影響狀況的動態(tài)監(jiān)控，這首先需要調(diào)查出各類污染源對這片地區(qū)的污染影響。而在高斯煙羽擴散模型上進行改進，即可模擬出污染源對周邊空氣影響的大致范圍及污染程度。而后，根據(jù)具體的污染程度及污染物類型，選取合適的位置設(shè)置監(jiān)控點，并制定相應(yīng)的周圍居民風險承擔經(jīng)濟補償方案。高斯煙羽模型計算各個方位、顆粒物的濃度計算相應(yīng)的空氣污染指數(shù)計算不同風向頻率確定賠償方法確定檢測方案圖1.對于問題一的分析流程圖4.1.2.模型建立（1）.高斯煙羽模型基本形式的數(shù)學表達式為： 公式（1）公式（2）其中為任意點的污染物濃度，為單位時間內(nèi)污染物排放量，為排放口的平均風速，為煙囪的有效高度，為煙囪的幾何高度，為污染物的抬升高度，與參數(shù)。污染物在排放出來時，具有初始動量，且受空氣浮力，所以污染物會有抬升高度。根據(jù)煙氣抬升公式，煙氣抬升高度為：公式（3）公式（4）公式（5）公式（6）公式（7）以上公式中，為煙氣熱釋放率，為實際排煙率，為排氣筒出口處污染物排出速度，為煙氣出口溫度，為煙氣出口溫度與環(huán)境溫度差，、分別為參數(shù)，具體值見下表：、的選?。▎挝唬篕J/s）地表狀況（平原）2100農(nóng)村或城市遠郊區(qū)1.4271/32/3城市及近郊區(qū)1.3031/32/32100240m時，按240m計算。P為風速高度指數(shù)，依賴于大氣穩(wěn)定度和地面粗糙度。應(yīng)根據(jù)觀測結(jié)果，利用統(tǒng)計學方法求出。根據(jù)具體的觀測數(shù)據(jù)，也可采用風速隨高度變化的對數(shù)律或其它半經(jīng)驗分式。無實測數(shù)據(jù)時，P值可按下表選取。由附件5推測，此垃圾焚燒廠在深圳市的鄉(xiāng)村。地區(qū)各穩(wěn)定等級下的P值A(chǔ)BCDE或F城市0.10.150.20.250.3鄉(xiāng)村0.070.070.10.150.25表2.各穩(wěn)定度等級下的P值在本文中，根據(jù)廣東沿海地區(qū)大氣穩(wěn)定度的分類方法探討，估算出深圳市大氣穩(wěn)定度等級頻率，具體數(shù)值如下表：穩(wěn)定度等級ABCDEF頻率1.833.8624.2748.6315.245.55表3.深圳市大氣穩(wěn)定等級頻率對于參數(shù)與的取值參考下表得出：大氣穩(wěn)定度A0.22/(1+0.0001)0.2B0.16/(1+0.0001)0.12C0.11/(1+0.0001)0.08/(1+0.0002)D0.08/(1+0.0001)0.06/(1+0.0015)E0.06/(1+0.0001)0.03/(1+0.0003)F0.04/(1+0.0001)0.016/(1+0.0003)表4.參數(shù)與的取值當污染熱釋放率大于或等于2100KJ/s，且污染物溫度與環(huán)境溫度的差值大于或等于35K時，污染物抬升高度采用公式（3）、（4）計算；當1700KJ/s300重污染健康人運動耐受力降低，有明顯強烈癥狀，提前出現(xiàn)某些疾病老年人和病人應(yīng)當留在室內(nèi)，避免體力消耗，一般人群應(yīng)避免戶外活動表5.API指數(shù)及其對應(yīng)的空氣質(zhì)量和人類生活的影響利用層次分析法，計算不同空氣污染等級對用的權(quán)重。設(shè)定API數(shù)值越高及其對應(yīng)的賠償重要性越重要。得出正互反矩陣：通過MATLAB計算得出相應(yīng)的賠償權(quán)重如下表：API3002003001002000100賠償權(quán)重0.46730.27720.16010.0954表6.不同API等級對應(yīng)的賠償權(quán)重經(jīng)計算，CI=0.0103，CR=0.0115，RI=0.90。CIRI且CR0.1，通過一致性檢驗。通過對附件4的風向統(tǒng)計，得出不同風向頻率，具體數(shù)值如下表：風向北東北東東南南西南西西北頻率0.0670730.0701219510.0335365850.0426830.0792680.3841460.2408540.082317表7.不同風向?qū)?yīng)的頻率假設(shè)當?shù)卣屠贌龔S共需對當?shù)鼐用褓r償金額，經(jīng)過資料查詢，得知深圳2011年人口密度為5359/平方公里，而經(jīng)過模型計算，污染范圍大約為78.5平方公里，估算受影響人口大致為420890人。由此得出具體居住地的每位居民應(yīng)得的賠償金額為。對于監(jiān)控點的設(shè)置，污染程度高的地方理應(yīng)設(shè)置更多的監(jiān)測點。根據(jù)MATLAB計算，不同方向的最大API值以及其距焚燒爐的距離見下表：方向北東北東東南南西南西西北最大API164.24141.57132.47134.55129.8175.68179.85165.32距焚燒爐距離1350m1250m1200m1200m1200m1450m1450m1350m表8.最高API值及其出現(xiàn)位置考慮到不同風向頻率的影響，最終給出具體監(jiān)測點的數(shù)量見下表：方向北東北東東南南西南西西北監(jiān)控點數(shù)量2個1個1個1個2個7個5個2個表9.各個方向設(shè)置的監(jiān)控點個數(shù)4.1.4模型的檢驗垃圾焚燒爐所排放的顆粒物，往往帶有汞、二惡英等重金屬及有毒物質(zhì)，所以顆粒物是周圍居民生活活動和生命健康的主要威脅污染物之一，并且顆粒物陳降到土壤里，也會對土壤產(chǎn)生嚴重影響。風速、風向及降雨對顆粒物的沖刷作用都會影響顆粒物的沉降。本文定義風雨影響指數(shù)來衡量顆粒物的沉降對土壤的污染。可衡量八個不同方位全年土壤因風速、風向、降雨而受的顆粒物的影響。的計算公式如下：=風向頻率對應(yīng)方位平均風速大小降雨量公式（9）.綜合考慮每個月不同風向的頻率得出矩陣，用來描述每個月不同風向的頻率。每個月不同風向的頻率再乘以對應(yīng)的風速，得到全年不同方向的土壤受風的影響，其中這里的表示一年的十二個月，表示從正北順時針到西北的八個方向。通過訪問深圳氣象局在2011年4月到2012年3月降雨量的資料（見附錄），計算得到一年不同方位的風雨影響指數(shù)。可用來描述八個方向土壤受垃圾焚燒爐顆粒物排放的影響。具體情況見下表：方向北東北東東南南西南西西北234.010 214.2890296.0116167.3187166.3451148.578161.993747.2536 表5.不同方向的風雨影響指數(shù)的值越大，則對應(yīng)方向的土地受顆粒物的污染越嚴重。由上表可看出，正東方向土地污染情況最為嚴重，而正北方向、東北方向土地污染情況也較嚴重。這與上文經(jīng)過高斯煙羽模型計算的污染物濃度分布情況一致，相互證明了所建模型的準確性。4.2問題二：4.2.1.問題分析本題考慮故障發(fā)生概率的情況下修正設(shè)計的監(jiān)測方法和補償方案，通過查閱相關(guān)文獻得到垃圾場對煙塵、二氧化硫、氮氧化物的去除率分別為99%，60%，25%7，發(fā)電廠年故障時間60h，通過調(diào)整源強，調(diào)用高斯擴散模型得到3種污染物的濃度范圍，繪制API等值線圖并對監(jiān)測方法和補償方案進行修正。4.2.2.模型建立與求解1）源強Q的確定煙塵、二氧化硫、氮氧化物的去除率分別為99%，60%，25%，設(shè)定在故障發(fā)生的情況下煙塵、二氧化硫、氮氧化物的源強分別為正常情況下的100倍，5/3倍，4倍。2）故障發(fā)生概率P的確定發(fā)電廠年故障時間3002003001002000100級別一級區(qū)二級區(qū)三級區(qū)四級區(qū)賠償權(quán)重0.46730.27720.16010.0954表6區(qū)域內(nèi)煙塵濃度嚴重超標，根據(jù)嚴重程度進行賠償。危險程度65544331級別一級區(qū)二級區(qū)三級區(qū)四級區(qū)賠償權(quán)重0.46730.27720.16