2011國賽A題優(yōu)秀論文

1、2011高教社杯全國大學(xué)生數(shù)學(xué)建模競賽承諾書我們仔細閱讀了中國大學(xué)生數(shù)學(xué)建模競賽的競賽規(guī)則.我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導(dǎo)教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反競賽規(guī)則的，如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽規(guī)則的行為，我們將受到嚴肅處理。我們參賽選擇的題號是（從A/B/C/D中選擇一項填寫）：我們的參賽報名號為（如果賽區(qū)設(shè)置報名號的話）

2、：所屬學(xué)校（請?zhí)顚懲暾娜簠①愱爢T（打印并簽名）：1.指導(dǎo)教師或指導(dǎo)教師組負責(zé)人（打印并簽名）：日期：2011年9月J2_日賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：2011高教社杯全國大學(xué)生數(shù)學(xué)建模競賽編號專用頁賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）：評閱人一評分備注OOOO全國統(tǒng)一編號（由賽區(qū)組委會送交全國前編號）：全國評閱編號（由全國組委會評閱前進行編號）：城市表層土壤重金屬污染分析摘要本文以城市表層土壤重金屬污染作為問題的背景，以所給出的樣值點重金屬污染物濃度為基礎(chǔ)，探求重金屬污染物空間分布、傳播特征、污染源位置等問題，我們通過對數(shù)據(jù)分析

3、，利用軟件處理數(shù)據(jù)，建立一系列模型,逐步解決題目提出的四個問題。針對問題一：我們首先對數(shù)據(jù)進行有效處理，利用matlab和GSplus專業(yè)繪圖軟件，畫出八種重金屬的空間分布情況。而后建立內(nèi)梅羅模型對所給數(shù)據(jù)進行量化處理，并與國家標準作比較，劃以等級，并用百分比去分析各個功能區(qū)內(nèi)不同等級的點的比例情況，最后用一個功能區(qū)內(nèi)各個點的期望去表征一個功能區(qū)的污染強度，可知工業(yè)區(qū)污染最嚴重，次之是交通區(qū)，山區(qū)污染程度最低。針對問題二：功能分區(qū)的不同使各個區(qū)內(nèi)污染有著各自的特征，以當(dāng)塊土地的背景值為基礎(chǔ)，加上標準差，建立單因素分析模型，同樣對各個點進行量化，并與背景值作比較，確定污染程度。并對8種重金屬在各

4、個區(qū)的分布情況進行橫向?qū)Ρ?，找出八種元素濃度值最高的區(qū)域，并結(jié)合當(dāng)?shù)毓δ軈^(qū)的特征進行詳細的介紹，具體說明了污染原因。針對第三問，我們根據(jù)圖形判斷其高斯分布特征，建立微分方程，并通過各向同性和各項異性建立方程，去確定污染源的位置，并給出了污染后的傳播特征方程。一共在五個區(qū)找到22個污染源，Cr元素污染源最少，同時給出各個污染源的坐標，具體見模型三。針對問題四，我們分外因和內(nèi)因去尋找新數(shù)據(jù)，去建立城市演變的灰色關(guān)聯(lián)模型，該模型更好的說明了污染物傳播和哪些因素有關(guān)，同時又給出了各個因素之間的關(guān)系，有著很好的前瞻預(yù)測性。最后，通過對四個問題的回答,我們找出了關(guān)于重金屬污染物傳播特性,同時對城市地質(zhì)結(jié)構(gòu)

5、有了一定的了解。根據(jù)所得污染物的情況，合理建議城市主管部門對城市合理分區(qū)，并有效治理污染！關(guān)鍵字：內(nèi)梅羅綜合指數(shù)GS+咼斯模型模糊分析1.問題重述隨著城市經(jīng)濟的快速發(fā)展和城市人口的不斷增加，人類活動對城市環(huán)境質(zhì)量的影響日顯突出。對城市土壤地質(zhì)環(huán)境異常的查證，以及如何應(yīng)用查證獲得的海量數(shù)據(jù)資料開展城市環(huán)境質(zhì)量評價，研究人類活動影響下城市地質(zhì)環(huán)境的演變模式，日益成為人們關(guān)注的焦點。按照功能劃分，城區(qū)一般可分為生活區(qū)、工業(yè)區(qū)、山區(qū)、主干道路區(qū)及公園綠地區(qū)等，分別記為1類區(qū)、2類區(qū)、5類區(qū)，不同的區(qū)域環(huán)境受人類活動影響的程度不同?，F(xiàn)對某城市城區(qū)土壤地質(zhì)環(huán)境進行調(diào)查。為此，將所考察的城區(qū)劃分為間距1公里

6、左右的網(wǎng)格子區(qū)域，按照每平方公里1個采樣點對表層土(010厘米深度)進行取樣、編號，并用GPS記錄采樣點的位置。應(yīng)用專門儀器測試分析，獲得了每個樣本所含的多種化學(xué)元素的濃度數(shù)據(jù)。另一方面，按照2公里的間距在那些遠離人群及工業(yè)活動的自然區(qū)取樣，將其作為該城區(qū)表層土壤中元素的背景值。附件1列出了采樣點的位置、海拔高度及其所屬功能區(qū)等信息，附件2列出了8種主要重金屬元素在采樣點處的濃度，附件3列出了8種主要重金屬元素的背景值?，F(xiàn)要求你們通過數(shù)學(xué)建模來完成以下任務(wù)：給出8種主要重金屬元素在該城區(qū)的空間分布，并分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。通過數(shù)據(jù)分析，說明重金屬污染的主要原因。分析重金屬污染物

7、的傳播特征，由此建立模型，確定污染源的位置。分析你所建立模型的優(yōu)缺點，為更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集什么信息？有了這些信息，如何建立模型解決問題？2問題分析隨著經(jīng)濟的發(fā)展和城市現(xiàn)代化的建設(shè)，人類的生產(chǎn)生活對環(huán)境的影響越來越大，工廠排放的廢棄物品，交通車輛的大量尾氣排放，生活污水的不合理處理，以及對自然的不加限制的開發(fā)，導(dǎo)致大自然自我抵抗能力下降，人類生活環(huán)境不斷惡化。本題實際上是根據(jù)此背景就城市土壤重金屬的污染進行探究性分析，以確定所給題目中城市環(huán)境重金屬的污染情況，并找出污染源，進行有效的治理。對于問題一，用金屬元素濃度代替空間三維分布中的海拔高度，用8張二維半圖可表示出8種元

8、素的空間分布，通過Excel分析每種金屬在各個區(qū)的平均濃度，標準偏差，范圍等參數(shù)，與我國土壤環(huán)境質(zhì)量評價標準作比較，運用尼梅羅綜合指數(shù)法可以對每個采樣點進行綜合污染指數(shù)評級，再分類匯總，評定各個區(qū)域的土壤被污染程度。對于第二問，分析各個區(qū)的每種金屬污染的嚴重程度，首先需要利用背景值求得各元素單因子污染指數(shù)評價分級標準，再對采集點污染等級進行分類匯總，求出各個等級在該功能區(qū)所占比例，比較最大的若干個功能區(qū)域，分析它們與重金屬污染有何聯(lián)系，即為污染主要原因。對于第三問，首先找出金屬濃度的最大值點，一般認為重金屬元素從濃度高處向濃度低處擴散，建立擴散模型模擬驗證。對于第四問，對于城市地質(zhì)環(huán)境的演變模

9、式，考慮的是它的影響因素，PH值，地形，歷史原因，發(fā)展規(guī)劃等，綜合考慮，預(yù)測演變趨勢。3.模型假設(shè)假設(shè)：不考慮pH值對各個重金屬污染分布的影響，即認為pH值對各個地區(qū)各種重金屬元素污染分布的影響均衡。在對樣點進行采集的過程中，該城市沒有自然災(zāi)害，土壤表層的特征保持不變。3在對樣點進行采集的過程中，該城市沒有大型工廠的遷入遷出。4.題目的采集點具有隨機的一般性，能夠充分反映城市的一些特征。符號說明字母代表含義P重金屬污染程度S.1重金屬污染評價標準C.11污染物的實測值P采樣點的綜合污染指數(shù)綜P采樣點所有單項污染指數(shù)的平均值平均P最大采樣點所有單項污染指數(shù)的最大值模型的建立與求解5.1模型一：內(nèi)

10、梅羅綜合污染指數(shù)模型5.1.1模型的建立分析對于第一問，要求給出八種重金屬的分布特征，同時分析不同城市區(qū)域內(nèi)污染程度。首先利用GSPlus地理作圖軟件對該城市八種金屬的分布，用濃度代替海拔高度畫出二維半圖，清晰表征八種元素分布特征。對于污染程度，我們選取內(nèi)梅羅法對該城市重金屬污染程度進行較全面的分析，并污染水平清潔土壤輕度污染土壤中度污染土壤污染相當(dāng)嚴重I通過與國家標準值進行比較，合理劃分等級，做出污染程度合理的評價。5.1.2模型的建立與求解首先根據(jù)國家標準去確立單向污染指數(shù):c,CcCI污染起始值污染起始值C-C+I污染起始值,CyCC中度污染起始值TOC o 1-5 h zC-C污染起始

11、值I中度污染起始值污染起始值C-C+I中度污染起始值,CC-c中度污染起始值重度污染起始值中度污染起始值C-C+I重度污染起始值,CC-c重度污染起始值重度污染起始值中度污染起始值p為采樣點i污染物的污染指數(shù)，c為i污染物的實測值,污染起始值，中度污染起始值，重度污染起始值的評價標準，此問中選用國家環(huán)境土壤質(zhì)量標準。根據(jù)此求出每一個采樣點的每一種重金屬污染后的P值，進一步去求P和P/平均最大：fPP1平均=8P最大二max(P)IP=（P2+P2）/2綜平均最大國家土壤環(huán)境質(zhì)量評價標準（mgkg-1）級別CDAsCuPBCrZnNI0.200.1515353590100400.300.5030

12、10030020025050三1.001.5404005003005001.0國家土壤重金屬污染分級標準污染等級一級二級三級四級PIP1IYP2IYP3P0.85綜0.85YP1.71綜1.71YP2.56綜PA2.56綜各個生活區(qū)的污染情況和比例:功能區(qū)1生活區(qū)2工業(yè)區(qū)3山區(qū)4交通區(qū)5公園綠地區(qū)采樣點個數(shù)44366613835各污染一級2511.1171.2118.8437.14等級百二級505027.2746.3842.86分比三級18.1811.111.5227.5417.14（%）四級6.8227.7807.252.86P綜期望1.5801911.992350.8456041.8470

13、991.2851925.1.3結(jié)果分析平均綜合污染指標：山區(qū)公園綠地區(qū)生活區(qū)交通區(qū)工業(yè)區(qū)，綜合污染程度可以由此判斷山區(qū)的土地重金屬污染程度最低，約四分之三的采樣點為清潔等級。約99%的采樣點為一、二級土壤等級，四級嚴重污染等級所占采樣點百分率為零，土壤環(huán)境質(zhì)量在霧中功能區(qū)中最高公園綠化區(qū)主要是一、二級輕度污染，總量占采樣點的五分之四，但比之山區(qū)，三級中度污染土壤比重明顯偏多生活區(qū)主要是二級輕度污染，占50%，其次是一級清潔等級土壤，占%25，嚴重污染等級占生活區(qū)采樣點的6.82%，比山區(qū)，公園綠地區(qū)的污染明顯嚴重交通區(qū)和工業(yè)區(qū)的污染是最嚴重的，交通區(qū)的二級，三級土壤等級比較多，約占總數(shù)的四分之

14、三，工業(yè)區(qū)則是嚴重污染的比重大，已經(jīng)超過了四分之一。綜上所述，工業(yè)區(qū)的四級嚴重重金屬污染土壤比重最大；其次是交通區(qū)中度污染比重最大；生活區(qū)土壤主要是輕度污染；工業(yè)綠化區(qū)的一、二級土壤等級在五個區(qū)中排第二位；山區(qū)的表層的土壤污染指數(shù)最低，是目前受污染最小。同時通過數(shù)據(jù)的量值加之地理特征，生活區(qū)作為居民生活的地方，其重金屬的污染程度較高是由于靠近工業(yè)區(qū)，并與交通區(qū)混連，受著來著工業(yè)區(qū)和交通區(qū)的雙重污染，故其較大的值是有地理的原因的。海拔三維圖5.1.4各種元素濃度分布圖18449.0.9551.19103.F212299.6150.28654.18449.12299.6150.AsCd18449.

15、12299.6150.9551.19103.28654.18449.0.12299.o6150.0.9551.19103.28654.F2CrCu還有四張圖見附件。5.1.5圖像分析上面為由GS做出的八種元素在全市區(qū)內(nèi)的一個濃度分布情況，可以看出，對于八種元素，基本上濃度最高的地方都是在圖的左下方，也就是工業(yè)區(qū)的地方。由分析可得出，工業(yè)區(qū)由于其本身的特性，各種重金屬的污染都相對較為嚴重，遠離工業(yè)區(qū)的山區(qū)是受污染最輕的地方。5.2模型二：單因子污染指數(shù)法5.2.1模型建立分析522模型的符號說明PiC.1S.1SX土壤中i污染物的污染指數(shù)土壤中1污染物的實測濃度1污染物的評價標準該地區(qū)背景值標準

16、偏差該地區(qū)背景值平均值5.2.3模型建立求解污染指數(shù)的確定：PC/SiiiS二x+si土壤污染分級：由于各元素在土壤中含量的差異性，此處采用能反映各元素在土壤中含量變化幅度的標準差的變化進行土壤污染分級：當(dāng)P-1時，該點未受i污染物污染，為清潔級當(dāng)1Pi(x+2S)/Si時，土壤中i污染物處于背景值范圍內(nèi)，為尚清潔級當(dāng)(無+2S)/SiPi(x+3S)/Si時，土壤中i污染物含量超出背景值范圍，為輕污染級。當(dāng)(無+3S)/SiPi(元+4S)/Si時，土壤中i已構(gòu)成重污染，為重污染8種主要重金屬元素的背景值得下表表格5-1：表格5-18種主要重金屬元素的背景值平均值標準偏差范圍評價標準兀素(x

17、)(S)(x土2S)(S.)As(pg/g)3.60.91.85.44.5Cd(ng/g)1303070190160Cr(pg/g)319134940Cu(pg/g)13.23.66.020.416.8Hg(ng/g)358195143Ni(pg/g)12.33.84.719.916.1Pb(pg/g)316194339Zn(pg/g)6914419783根據(jù)模型處理數(shù)據(jù)，將污染程度分為5個等級，可得：表格5-2土壤中各元素單因子污染指數(shù)評價分級標準表元糸IIIIIIVAs1.6Cd1.5625Cr1.675Cu1.643Hg1.558Ni1.708Pb1.486Zn1.506污染程度清潔尚清

18、潔輕污染中污染重污染模型求解：利用上述模型對各個采集點8種重金屬元素的實測值進行計算，將計算結(jié)果與表4中的P標準值進行比較，確定該點的污染程度，對每個地區(qū)的各種元素進行百分比統(tǒng)計i表格5-3各地區(qū)土壤中重金屬元素污染狀況統(tǒng)計表兀素清潔(%)尚清潔(%)輕污染(%)中污染(%)重污染(%)As生活區(qū)20.4511.362511.3631.82工業(yè)區(qū)2513.898.3311.1141.67山區(qū)71.2112.126.061.529.09交通區(qū)33.3316.6721.0113.0415.94公園綠地區(qū)205.7125.7122.8625.71Cd生活區(qū)22.734.549.0913.6450工業(yè)

19、區(qū)5.568.338.335.5672.22山區(qū)66.677.5610.614.5510.61交通區(qū)17.395.805.0711.5960.14公園綠地區(qū)37.148.57202.8631.43Cr生活區(qū)27.2727.2725020.45工業(yè)區(qū)41.6716.6713.898.3319.44山區(qū)69.707.586.067.589.09交通區(qū)35.5128.2615.225.8015.22公園綠地區(qū)45.7128.5717.1408.57Cu生活區(qū)11.3611.3611.3611.3654.55工業(yè)區(qū)2.782.785.5616.6772.22山區(qū)57.5816.6710.616.069

20、.09交通區(qū)9.427.977.9710.1464.49公園綠地區(qū)17.1414.2911.438.5748.57Hg生活區(qū)38.644.546.826.8243.18工業(yè)區(qū)22.228.332.785.5661.11山區(qū)68.187.587.589.097.58交通區(qū)39.137.256.524.3542.75公園綠地區(qū)42.865.712.865.7142.86Ni生活區(qū)34.0938.6411.369.096.82工業(yè)區(qū)44.4413.8911.118.3322.22山區(qū)68.1812.126.063.0310.61交通區(qū)46.3832.6112.327.251.45公園綠地區(qū)68.57

21、17.145.715.712.86Pb生活區(qū)31.8213.646.829.0938.64工業(yè)區(qū)11.118.338.332.7869.44山區(qū)75.769.093.036.066.06交通區(qū)24.6413.045.075.0752.17公園綠地區(qū)42.8617.148.572.8628.57Zn生活區(qū)27.2711.366.826.8247.73工業(yè)區(qū)5.5616.678.335.5663.89山區(qū)74.2412.127.581.524.55交通區(qū)16.6711.5910.142.9058.70公園綠地區(qū)37.1414.295.7114.2828.57524模型的結(jié)果原因分析As污染主要是工

22、業(yè)區(qū)嚴重污染土壤中As主要來自于化肥工業(yè)廢物磷石膏，因為其可以改良酸性土壤，從而被大量使用，造成土壤中As的含量增加，它也會造成Cr、Pb的含量增加。Cd污染主要來自于工業(yè)區(qū)和交通區(qū)工業(yè)區(qū)主要原因可能是化工冶金工業(yè)與土壤中Cd超標密切相關(guān)交通區(qū)主要是汽車潤滑油在高溫條件下與空氣中的氧氣發(fā)生氧化反應(yīng)生成有機酸、醇、醛及其他有機化合物，這些物質(zhì)能夠腐蝕與其接觸的含Zn、Cu、Cd等合金的金屬部件以及含Cu元素等的油泵，導(dǎo)致其向環(huán)境中釋放相應(yīng)的重金屬元素，導(dǎo)致交通區(qū)Cd超標比例較大。Cr主要來自于工業(yè)區(qū)和生活區(qū)工業(yè)區(qū)主要為化工冶金工業(yè)，Hg、Cr、Cd超標與之密切相關(guān)，且有一些化肥工業(yè)廢物含有Cr元

23、素，如磷石膏，施入土壤后回事Cr含量增加生活區(qū)主要為城市灰塵，其中含Cr元素，且認為受活動干擾嚴重Cu主要來自于工業(yè)區(qū)、交通區(qū)和生活區(qū)工業(yè)區(qū)主要是一些冶銅工業(yè)以及其他的銅相關(guān)產(chǎn)業(yè)。交通區(qū)主要是汽車輪胎及排放的廢氣中含有Pb、Zn、Cu等多種重金屬元素，Pb、Zn和Cu等被看作機動車輛污染源的標志元素，所以Pb、Zn和Cu元素被認為主要來自于交通運輸過程，Pb、Zn和Cu被稱為交通因子。生活區(qū)城市灰塵中受人為活動干擾嚴重的重金屬為Pb、Cu、Cr、ZnHg主要來自于工業(yè)區(qū)，生活區(qū)，公園綠地區(qū)和交通區(qū)汞在工業(yè)上應(yīng)用很廣，造成的污染較嚴重，汞工業(yè)區(qū)主要來自燃煤電廠和金礦開采，有色金屬生產(chǎn)，鋼鐵生產(chǎn)等

24、生活區(qū)、公園綠地區(qū)、交通區(qū)，主要是汞產(chǎn)品在人們的生活中隨處可見。比如一個小的電開關(guān)、電池、熒光燈，這些含汞器械、設(shè)備與人們?nèi)粘Ｉ罱佑|很多，又容易被各界所忽略，就成為一個值得密切關(guān)注的汞污染源。另外，醫(yī)療和城市廢物被燃燒后也會造成一定的汞污染Ni主要來自工業(yè)區(qū)工業(yè)區(qū)主要是鎳礦的開采和冶煉；合金鋼的生產(chǎn)和加工過程；煤、石油燃燒時排放煙塵中；電鍍、鍍鎳等生產(chǎn)過程Pb主要來自工業(yè)區(qū)和交通區(qū)工業(yè)區(qū)中，工廠對周邊環(huán)境的Mn、Zn、V污染有很大的貢獻。汽車輪胎及排放的廢氣中含有Pb、Zn、Cu等多種重金屬元素，進入周圍的土壤環(huán)境,容易造成土壤重金屬污染。生活區(qū)土壤的Pb含量嚴重污染等級占比例也很大，主要是

25、城市灰塵中重金屬Pb，它受人為活動干擾嚴重Zn主要來自交通區(qū)、工業(yè)區(qū)、生活區(qū)在交通區(qū)，汽車輪胎添加劑中的Zn是城市土壤中Zn的重要來源由于人類的生產(chǎn)、生活等活動，Cu應(yīng)用廣泛，導(dǎo)致居民區(qū)和商業(yè)區(qū)主要以Cu和Zn的積累為特征。上述橫向的從元素的角度進行了分析，由特定的重金屬元素去分析功能區(qū)的特點，并簡單的猜測該功能區(qū)可能有哪些原因造成，可以說相信的闡述了重金屬污染主要原因的問題。5.3高斯模型5.3.1模型分析通過對各個區(qū)域不同金屬的濃度三維擬合，與高斯模型1具有很大的一致性，所以此題根據(jù)高斯模型建立濃度擴散方程找到污染源，我們首先建立無差異的傳播模型，即各個方向均勻傳播。q=-6*gradCt

26、JJq*nddt5.3.2模型建立Q=J1tQ_JJJC(X,y,z,t+At)一C(x,y,z,t)dV2Q于JJJJPdVdt00又根據(jù)曲面積分的Gauss公式:JJq-ndb=JJJdtJJqndb_JJJdivqdVAt.JJJC(x，y,21+At)-C(x，y，z，t)AtdV+J-0,-sx,y,zQtiiQt2yQt2zQt2這是無界區(qū)域的拋物線型偏微分方程，根據(jù)假設(shè)1,初始條件為作用在坐標原點的電源函數(shù)，記作C(x,y,z,0)_Q06(x,y,z)最后得到二維高斯分布原始公式：G二Aec（x_a）2+c2（y-b）2+c3（xa）（y-b），由于假設(shè)無差異傳播，（xa）2丄

27、（y_b）2所以x,y不相關(guān)，。3=0,cl=c2=l/c建立咼斯模型：G=Aec，其中A,a，b，c均為待求參數(shù)，G為（x,y）所在點的濃度,A為污染源處的濃度值該模型從坐標（a，b）處向四周均勻擴散，沒有差異性，則最終求得的（a,b）點即為極值點坐標。5.3.3模型求解：分別帶入極值點周圍數(shù)點，用計算機擬合可得高斯模型方程參數(shù)最優(yōu)解，解如下:As元素有三個污染源污染源abcA118205100895842030.13212695302513602434750729827332023Cd元素有四個污染源污染源abcA12145011375314801640212370208039382890

28、3333660691484245211404271523002286281140Cr元素有一個污染源污染源abcA13305602480543630925Cu元素有兩個污染源污染源abcA1330560239651613702237936882153722534Hg元素有三個污染源污染源abcA11524891603464138002137012360216599150103213038902937216000Ni元素有兩個污染源污染源abcA1221961219058777623305602410898144Pb有三個污染源污染源abcA147704905240804802237537001

29、579940031655274070223460Zn有四個污染源:污染源abcA11380396272311237702127103035127741660393154321339513000432906022292901880上述的啊，a,b代表著污染源的坐標，具體的圖如下:21.81.61.410.00.60.40.20.511.522.534x10X10WOO1.15350030002500250020001500200030001000800600400200040003000250020001500100050004000600020005000400021.81.61.41：210

30、.00.60.40.2000.511.522.53.x-10412001000-800-600-4001000200030003000400050005.4基于城市演變的灰色關(guān)聯(lián)模型5.4.1模型的建立分析縱觀上文，我們在根據(jù)題目的要求下，逐步分析該城市重金屬的空間分布，污染程度和污染原因。從宏觀入手建立了重金屬污染物的傳播特征的模型，但是不同金屬元素在土壤中不同形態(tài)的分配情況，以及自然原因等對金屬富集產(chǎn)生影響。同時，有上面建立的一系列模型可知，不同功能區(qū)的重金屬樣值有著不同的特征，且相差很多。說明城市土壤重金屬含量除了與自然環(huán)境有關(guān)，還與城市功能有關(guān)密切的關(guān)系。因此我們建立灰色關(guān)聯(lián)度模型，對

31、影響城市土壤重金屬分布的內(nèi)因和外因與重金屬含量進行灰色關(guān)聯(lián)度分析。5.4.2模型建立進一步收集影響城市土壤重金屬分布的內(nèi)因外果內(nèi)因：土壤有機質(zhì)，土壤粘性顆粒含量，土壤PH值；外因：生產(chǎn)生活的廢棄物，工業(yè)區(qū)遠近，交通車流量。對這些因子按110作等級評價，形成土壤有機質(zhì)、黏貼顆粒、土壤PH值、生產(chǎn)生活的廢棄物、工業(yè)區(qū)遠近、交通車流量的六大城市土壤重金屬含量影響因子。X1X2X3X4X5X6設(shè)六大因子分別為XX2X3X4X5,X6土壤有機質(zhì)含量(g/Kg)土壤顆粒含量()PH值(生產(chǎn)生活廢棄物占土壤總量的百分比)(相對距離)(相對流量，無量綱)生產(chǎn)生活廢物比例離工業(yè)區(qū)遠近相對交通車流量4根據(jù)所選擇的

32、樣本點，進行主要因子分析剖面號X112表格5-4主要因子分析表X2X3X4X5X6通過建立這樣的一張圖，直觀反映各因素與重金屬的關(guān)系，并通過他們與重金屬做關(guān)聯(lián)度分析，搞清楚各因子間的主次關(guān)系，更有利于城市土壤的生態(tài)改良和修復(fù)。5.灰色關(guān)聯(lián)系數(shù)的確定根據(jù)內(nèi)梅羅計算所得的P，X=P/P，同樣將X,X,X,X,X作比較數(shù)列，得TOC o 1-5 h zij0ij1123456IX(0)-X(k)l數(shù)據(jù)列表，對于X在k剖面的關(guān)聯(lián)系數(shù)8(k)有：0108二(MinMinIX(0)-X(k)I+九MaxMaxIX(0)-X(k)I)/(IX(k)-X(k)l)+九MaxMaxIX(k)-X(k)Iiik01ik0101ik01式中8(k)為k剖面污染指數(shù)與i因子的灰色關(guān)聯(lián)系數(shù)，九為分辨系數(shù)，X(k)為k剖面1ii因子的影響值。6確定R評價標準R=2(k)/nii=1用R值去表征關(guān)聯(lián)度，R值越大，其關(guān)聯(lián)度越大。5.4.3模型分析本模型是基于內(nèi)梅羅模型進行的擴展，選取了壤有機質(zhì)、黏貼顆粒、土壤PH值、生產(chǎn)生活的廢棄物、工業(yè)區(qū)遠近、交通車流量六大因