搜索
2012。[1]319、8、、、1～5。8。。。。。：2012-03-21：2011。（，210044）：Surfer、。：：X51：A：1008-5327（2012）02-0064-05PollutantsPropagationModelBasedonDiffusionEquationHUJieGAOMing-jieJINYi-shuai(NanjingUniversityofInformationEngineering,Nanjing210044,China)Abstract:Thispaperdrawsupurbanheavymetaldistributionmaps,accordingtosamplepointdata,usingSurfersoftwareandKrigingInterpolationMethodtooptimizeinterpolation.ItconstructstheassessmentsystemofurbansoilpollutionbyadoptingCarmeloCompositeIndexAssessment,takingtheimpactofdifferentheavymetalelementsonhumanhealthintoconsideration,andrevisestheweightsofassessment.Itexplainsthemainreasonofheavymetalpollutionthroughtheuseofclusteranalysismodelandanalysisofvariancemodel.ItestablishesPartialDifferentialEquationModelsforpollutantdensitydiffusion.Thelocationofpollutionsourcescanbefoundthroughzonedivision,equationsimplificationandLeastSquaresMethod.Theauthorproposestheideatopromoteandexpanditthroughtheanalysisoftheadvantagesanddisadvantagesofpollu－tantdispersionequation.Keywords:urbansoils;heavymetal;pollution;diffusionequation;nonlinearfittingVol．26No．2Jun.201226220126!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!JOURNALOFNANTONGVOCATIONALCOLLEGEdoi:10.3969/j.issn.1008-5327.2012.02.014642。11.1SPSSPerson1。8α=0.01HgAs、Cr、Ni。。1.2SurferSurfer3198。Zx0=Ni=1ΣμiZ（xi）1Zx0Z（xi）Nμii[2]。12CrCu6。AsCdCrCuHgNiPbZn。1.3。、、1994。Pz=P2max+P2have2姨Pmax=maxPiPhave=1nni=1ΣPi姨姨姨姨姨姨姨姨姨2PznPi=Ci/SiiCiimg/kgSimg/kgCdCrCuHgNiPbZnAs0.260.190.160.060.320.290.25Cd0.350.400.270.330.660.43Cr0.530.100.720.380.42Cu0.420.500.520.39Hg0.100.300.20Ni0.310.44Pb0.5011Cr2Cu960900840780720660600540480420360300240180120600151050km510152025km28002400200012001000800600400200140120100806040200151050km510152025kmμg/gμg/g652012[1]。1.42。。[3]。GB15618－1995Bii=12…nBi-1iτi=Bi-1B1-1+B2-1+…+Bn-13832。P′have=ni=1ΣτiPi。PmaxPiPsP′max=0.5Pmax+Ps。P=P′2max+P′2have2姨41.51“”12]“”23]“”3“”。3194P3。352.2%、77.8%、10.6%、60.7%、40%＞＞＞＞。5P3.21、10.18、1.31、7.01、2.62。。2。2.1。X1X2…XnY1Y2…Yn[4]ni=1ΣXi-Yi5。SPSS8AsNiCrPbCuCdZnHg。5As。CdZn2/mg·kg-1As300.0062Cd0.30.6168Cr250.0009Cu2000.0009Hg0.50.3701Ni500.0037Pb3000.0006Zn2500.00073444171013360811176623367013864829553521977662CdZn。Hg、HgHgHg。2.2Xijiji=12…5j=12…niXij=μ+ai+ξij。ai=μi-μμ=1n5i=1Σniμin=5i=1Σniξij~N0σ2。。SAS84。“2∧3”23。4HgNiCdCd。33.1。13。。。3.2、、。、、、、[5]坠Cij坠t+v（坠Cij坠x+坠Cij坠y）=Ej△Cij+fijxyt-λjCij6CijijvEjjfijxyttijλjj。v=0、λj=06Ej△Cij+fijxy=0。Ej-△Cij=fijxyCij=Cijxy。。-△Cij=Nm=1Σβmδr-rmrrmβm。4As2∧12∧32∧43∧43∧5Cd1∧21∧32∧32∧53∧43∧54∧5Cr1∧3Cu1∧22∧32∧42∧5HgNi2∧3Pb1∧21∧32∧32∧42∧53∧43∧5Zn1∧32∧33∧46720125/km/kmAs2.446.84As12.92.84As18.410.1Cd21.611.6Cr4.484.58Cr3.096.12Cu3.275.98Hg3.123.91Hg14.12.88Hg15.09.79Ni3.135.86Pb2.373.66Pb3.046.94Zn12.86.84Zn13.89.62Gr=-12πrijCij=-Ω乙GrNm=1Σβmδr-rm=-Nm=1ΣβmGr-rmCij=12πNm=1Σβm（x-xm）2+（y-ym）2姨劬F（x，y）73.3Cij3Nxyminsk=1Σ[Cij（xk，yk）-F（xk，yk）]28Cij（xk，yk）（xk，yk）s。N=1231、23。MatlabNlinfit[6]。155。5As、Cr、Cu、Ni、Pb、、。44.1、。15。。4.2Ejλj。坠Cij坠t+v（坠Cij坠x+坠Cij坠y+坠Cij坠z）=Ejxyz+fijxyzt-λjxyz+kuuCij-aρ9ukuρa。。LQldlxydlC=L乙βrQldlS（74）682012QsdσxyC=βrQsdσ：[1]．2011CUMCMA[EB/OL]．[2011-09-09].．[2]．ArcGIS[J]．20081236-38.[3]．[J]．2009648-50.[4]．[M]．2003.[5]．[J]．20104574-575.[6]．Matlab[M]．2010.乙乙S。4.31μσ2。23。3“”。：[1]．2011CUMCMA[EB/OL]．[2011-09-09].．[2]．[D]．2008159-163.[3]．[D]．200873-77[4]．[D]．200742-69[5]．SPSS[M]．2003．[6]．[EB/OL]．2011-08-12[2011-09-09].=like/．[7]．[D]．200930-34．（68）乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙乙74