2011年全国大学生数学建模竞赛A题优秀论文范文模板参考资料-城市表层土壤重金属污染分析

1、城市表层土壤重金属污染分析摘要重金属污染是当今世界关心的又一大问题，自从2021年“血铅超标事件发生后，国家越来越重视重金属对我们生活的影响。因此，本文对某地区的八种重金属的含量分五区进行分析。得到了该地区的重金属扩散模型和参数识别模型。用matlab插值法建立三维模型。得到重金属元素的浓度分布图，算出不同区域内各种重金属元素的污染程度。分析重金属的传播特征。从高浓度区向低浓度区扩散。我们建立扩散模型，求出函数的极值，从而确定污染源的位置。数据以后建立因子分析法，回归分析，曲线拟合等模型解决问题。关键字： MATLAB软件 数据拟合 曲线拟合三维模型内梅尔指数 插值法、均值法、扩散模型、因子分

2、析、回归分析。 一 模型分析问题背景河流、湖泊是人类生活用水的主要来源。2021年全国的环境状况来看，重金属污染形势也比拟严峻。因此通过对某地区的八种重金属的含量进行分析得出该地区的污染程度。问题分析对于问题一，用附件一中给出的数据，用matlab插值法建立三维模型，总共有9个图，一个是取样地点的地形图，另外八个是八种重金属元素的浓度分布图，用matlab插值法建立三维模型然后通过均值法，算出不同区域内各种重金属元素的污染程度。对于问题二，通过对问题一结论的分析得出，生活区和工业区是污染比拟厉害的地区。目前我国由于在重金属的开采、冶炼、加工过程中，造成不少重金属如铅、汞、镉、钴等进入大气、水、

3、土壤引起严重的环境污染。人类生活中各种用品都含有不同量的重金属元素，比方说废旧电池，含有较多的汞、铬、锰、铅、镍、锌等重金属。它们通过自然和生物降解，随着雨水进入到土壤和河流当中。对于问题三，根据前两问的结论分析重金属的传播特征，主要有从高海拔到低海拔，从高浓度区向低浓度区扩散。我们建立扩散模型，求出函数的极值，从而确定污染源的位置。对于问题四，我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式，还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，以及每年的生物降解量，降雨量对重金属元素扩散的影响，空气污染也应该考虑进去。有了这些数据以后建立因子分析法，回归分析，曲线拟合等模型

4、解决问题。二 模型假设 1污染源的重金属浓度不在增加。(2) 取样点的数据较好的反映了该地区的污染物浓度。3测量的个别数据对整体没有影响。4元素的扩散只与元素的含量的上下有关 三 模型参数假设i=1,2,3,8,分别对应8种金属元素，x表示距离参照点的横坐标y表示距离参照点的纵坐标表示i种金属的x,y处的含量j=1,2,3,4,5,表示种元素在j区的平均含量该地区i种元素的背景值k为考虑各差异可能会引起背景值的变动而取的系数这里取1.5。为土壤中污染中污染元素的实测值；为土壤中污染元素的背景值。分别代表土壤污染积累起始值、中度污染起始值和重度污染起始值D ( x, y , z ) 为扩散系数S

5、 为曲面质量 四 模型建立对土壤环境质量进行评价时所采取的评价方法多种多样,主要有污染指数法、污染程度法, T 值分级法、基准分级法、密切值法和模糊数学综合评判法、灰色聚类法等。由于地域差异等各种因素的影响,污染指数法没有统一的评价标准。通常进行土壤环境质量评价时，参照中国土壤环境质量标准，本文采用具体的单因子污染指数法如下9： 1与的计算关系如公式2所示。 2根据式2计算得到污染指数值来对重金属的污染程度进行评价，其中为非污染状态；为轻污染状态；为中度污染状态；为重污染状态。为了全面、综合地反映多种污染物的整体污染水平，因此,需要一种同时考虑多种污染物综合污染水平的多因子评价方法,即将单因子

6、污染指数按一定方法综合。常用的方法有内梅罗指数法,它兼顾了单因子污染指数的平均值和最高值,能较全面地反映环境质量,而且可以突出污染较重的污染物的作用,其计算公式如下： 3式中：为土壤污染元素综合污染指数；为土壤中污染元素的实测值；为土壤中污染元素的评价标准。同单因子污染评价程度一样，为非污染状态；为轻污染状态；为中度污染状态；为重污染状态。下表为比照值地积累指数lego分级污染程度5<lego<=106极严重污染4<lego5强-极严重污染3<lego4强污染2<lego3中等-强污染1<lego2中等污染0<lego1轻度-中等污染lego0无污染

7、使用内梅罗指数评价各样点土壤污染程度的分布情况评价范围污染程度未污染轻污染中污染重污染样本数0081表1 土壤污染评价标准的建议范围与数值元素污染积累起始值中度污染起始值重度污染起始值HgCr90250400CdCu3550400Pb35250500As153060五 模型求解运用matlab插值法建立三维模型得到九幅图问题一：该地区的空间模拟图得出五区的对应位置为下表功能区坐标值*m生活区1xy工业区0x10y山区x3y2交通区1x1y公园绿地区x30y该地区As的含量图该地区Cd含量图该地区Cr含量图该地区Cu含量图该地区Hg含量图该地区Ni含量图该地区Pb含量图该地区Zn含量图通过观察知

8、道五区的污染程度得出八种元素的污染源分别为元素AsCdCrCu污染区工业区、生活区工业区工业区工业区元素HgNiPbZn污染区工业区、生活区工业区工业区、生活区工业区、生活区问题二 用与的计算关系如下公式用与的计算关系如下公式求出五区八种元素的平均含量得出五区的金属含量和土壤污染评价标准的建议范围与数值对比照。得出污染程度污染程度表功能区HgNiPbZnlego分级lego分级lego分级lego分级工业区极严重污染极严重污染极严重污染极严重污染商业区中等-强污染轻度-中等污染轻度-中等污染中等-强污染居民区中等-强污染强污染强-极严重污染强-极严重污染农业区轻度-中等污染轻度-中等污染轻度-

9、中等污染轻度-中等污染旅游区无污染无污染无污染无污染功能区AsCdCrCulego分级lego分级lego分级lego分级工业区极严重污染极严重污染极严重污染极严重污染商业区中等-强污染轻度-中等污染轻度-中等污染中等-强污染居民区轻度-中等污染轻度-中等污染强-极严重污染轻度-中等污染农业区轻度-中等污染中等-强污染轻度-中等污染中等-强污染旅游区轻度-中等污染无污染轻度-中等污染无污染问题三：此题和研究气体的扩散,液体的渗透,半导体材料中的杂质扩散等问题所满足的微分方程相同 扩散物质在单位时间内沿法线方向 n 流过单位面积的曲面的质量与物质浓度 C ( x, y, z , t ) 沿法线方

10、向 n 的方向导数 C 成正比. n 由扩散定律,扩散物质在时段 dt 内沿法线方向 n 流过面积为 dS的曲面的质量 dm 为: dm = D ( x, y, z ) C dS dt n ,出现负号是由于物质总是由浓度高的一侧向浓度低的一侧渗 透. 任取一封闭曲面 ,它所围区域记为 ,那么从时刻 t1 到时刻 t2 进入此闭曲面的物质质 量为由高斯公式 同时,物质渗透到区域 内,使得内部的浓度发生变化,在时间间隔 t1 , t1 内,浓度 由 C ( x, y, z , t1 ) 变化为 C ( x, y, z , t2 ) ,增加的物质质量为由质量守恒即有于是得到扩散方程 假设扩散系数 D

11、 ( x, y , z ) 为常数,那么扩散方程为得出污染物是从高浓度传向低浓度的所以我们知道了各种重金属污染物的污染源位置为 总的来说：工业化程度越高的地区污染越严重，市区高于远郊和农村，地表高于地下，污染区污染时间越长重金属积累就越多，以大气传播媒介土壤重金属污染土壤的具有很强的叠加性,熟化程度越高重金属含量越高。元素AsCdCrCu污染区工业区、生活区工业区工业区工业区元素HgNiPbZn污染区工业区、生活区工业区工业区工业区、生活区问题四：对于问题四，我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式，还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，以及每年的生物降

12、解量，降雨量对重金属元素扩散的影响，污水排放，沙漠化的程度，空气污染也应该考虑进去。有了这些数据以后建立因子分析法，回归分析，曲线拟合等模型解决问题。六 模型的评价 1本文把所解决的问题归结为优化问题，建立的数学模型清晰合理。 2运用MATLAB和LINGO软件处理数据和进行运算，降低运算量，简单易行，有很大的可操作性。且所得数据较为合理可靠。 3运用因子分析模型解题，全面可靠。4但在实际运用本方案中还应考虑后来污染因素对模型的影响，在应用的过程中根据实际情况进行灵活改变。七参考文献 1 肖华勇.?基于MATLAB和LINGO的数学实验?M.西安：西安工业大学出版

13、社，2021(3).1142. 张志涌.?精通MATLAB 6.5版?M.北京：北京航天航空大学出版社，2003(3).3133. 邢文训，谢金星. ?现代优化算法?M. 北京: 清华大学出版社, 1999.4 .姜启源 谢金星 . 叶俊 ?数学模型?M，北京，2003.5. 濮定国、田蔚文主编，数学模型，东南大学出版社，1999. 6. 欧阳亮，系统科学中数学模型，山东大学出版社，2003. 7. 陈义华，数学模型，重庆大学出版社，1999. 8. 朱思铭，李尚廉，数学模型，中山大学出版社，1995. 9. 蔡常丰，数学模型建模分析，科学出版社，1995. 10. 徐全智，杨晋浩，数学建模入

14、门，电子科技大学出版社，1996. 11. 沈继红、施久玉、高振滨、张晓威，数学建模，哈尔滨工程大学出版社，1996. 附录：. x=74，1373，1321，0，1049，1647，2883，2383，2708，2933，4233，4043，2427，3526，5062，4777，5868，6534，5481，4592，2486，3299，3573，4741，5375，5635，5394，5291，4742，4948，5567，7004，7304，7048，8180，9328，9090，8049，8077，8017，6869，7056，7747，8457，9460，9062，9319，106

15、31，10685，10643，11702，11730，11482，10700，10630，11678，11902，13244，12746，12855，13797，14325，15467，12442，13093，13920，14844，16569，16387，16061，15658，14298，14177，15092，12778，17044，17087，17075，17962，18413，19007，18738，17814，18134，17198，17144，18393，19767，21006，21091，22846，23664，22304，21418，21439，20554，20211，210

16、72，20215，18993，19968，21766，22674，22535，25221，26453，26416，27816，25361，24065，25998，27177，26424，26073，24631，24702，25461，24813，26086，26015，27700，27696，27346，26591，27823，27232，24580，24153，22965，23198，24685，28654，24003，21684，22193，17079，15255，15007，3518，3469，3762，3927，4153，3267，4684，5495，5664，5541，5451，40

17、20，4026，5101，5438，5382，5314，5503，5636，6605，7093，7100，6837，7906，8045，8394，8403，8079，9663，9469，9178，9095，10225，10210，10340，11557，11415，11649，12734，12696，12400，12591，13765，13694，13855，14862，14896，15387，15810，16032，15801，15087，16872，17734，16823，17008，17203，17005，16947，16301，17904，18303，18438，18556，18954

18、，18012，19072，20282，21475，21450，20261，19569，19411，19501，20582，19909，21018，22176，23359，23238，22624，21703，5006，5734，6395，7405，8446，7612，7912，8866，9296，9475，9212，8629，7776，8622，9237，8307，7106，6423，7458，8904，10547，10398，10395，11529，11563，11646，12641，14000，14207，14065，12734，12727，14173，15467，15140，15198，1

19、5248，16428，16289，16267，16440，16440，15412，14269，13277，13175，12153，11958，10800，10022，9333，9277，11121，10856，12644，12625，9036，10599，12632，14405，14074，14262，14624，16629，18470，20591，20983，20217，19041，18906，18467，17414，15748，15517，16607，15952，22605，23146，22046，23785，25981，27380，25021，23325，26852，17981，1448

20、2，14318，10352，9095，10510，13954，10142，17765，6924，4678，6182，5985，7653;y=781，731，1791，1787，2127，2728，3617，3692，2295，1767，895，1895，3971，4357，4339，4897，4904，5641，6004，4603，5999，6018，6213，6434，8643，7965，8631，7349，7293，7293，6782，6226，5230，4600，4496，4311，5365，5439，6401，7210，7286，8348，8260，8991，8311，7639，679

21、9，6472，5528，4472，4480，5532，6354，8184，8774，8618，7709，7056，8450，8945，9621，8666，8658，4329，4339，5354，5519，6055，6609，7352，7594，7418，6684，6936，5799，10691，11933，12924，12823，11721，11488，10921，10707，10046，9810，9081，9183，8810，8819，9482，9149，9790，10527，10721，11383，11228，10774，10404，9951，12371，12961，12348，12173

22、，11293，5795，5577，6508，5581，6423，7353，7032，7771，8639，8807，9422，9522，9834，10799，11094，12078，11609，11621，13331，13715，14737，14482，13319，12450，13535，13523，14278，8755，15286，13101，12185，5894，5110，5535，2571，2308，2170，2110，2299，793，1364，1205，1653，2093，2757，2990，3913，4080，3994，3012，2060，1127，133，374，1381，2449

23、，3490，3978，3052，2035，1075，0，1288，2286，3299，3975，3821，2789，1764，1581，2585，3515，4015，3024，2060，1063，1353，2357，3345，2524，1603，729，2307，3061，3966，3512，2798，3629，4207，4775，6218，7212，7487，8299，8287，7385，6539，5588，4874，4414，8519，8590，8540，7555，7586，7348，6934，6091，6548，5300，5764，5492，5325，6502，4818，6591，884

24、6，9659，10443，10981，11200，11938，12840，13143，13102，12000，11305，12086，10613，10638，9872，9726，9467，8831，8920，8868，9591，10360，11203，11243，10298，9381，9560，8970，9980，10987，10344，7691，11941，12080，11101，10100，9106，9069，10072，11058，12068，13232，12982，12877，13204，12238，12336，13313，13282，12204，14631，16148，16432，1

25、4727，14943，16259，17538，17980，17949，18032，16516，15129，14004，14481，14411，13549，15862，17642，15769，16346，17001，15476，15728，17034，17365，18397，14301，15382，17634，17643，18051，18202，16290，16701，16114，18449，12692，13569，17133，16414，15314，5615，1662，3561，5696，3765，2005，2567，1952;， 5.93 ，4.90 ，6.56 ，6.35 ，14.08 ，

26、8.94 ，9.62 ，7.41 ，8.72 ，5.93 ，9.17 ，5.72 ，4.49 ，5.51 ，11.45 ，6.14 ，7.84 ，7.41 ，8.50 ，5.51 ，9.84 ，9.39 ，3.30 ，4.09 ，6.14 ，5.31 ，3.69 ，21.87 ，18.38 ，10.53 ，3.50 ，6.35 ，5.51 ，4.49 ，3.50 ，5.51 ，4.29 ，4.29 ，6.56 ，16.58 ，7.41 ，5.93 ，4.69 ，4.90 ，5.31 ，4.29 ，5.51 ，4.69 ，7.20 ，5.31 ，4.90 ，4.90 ，3.89 ，3.69 ，3

27、.11 ，3.89 ，3.89 ，2.91 ，3.30 ，4.90 ，4.09 ，2.91 ，2.72 ，3.11 ，3.30 ，3.30 ，6.14 ，3.69 ，4.49 ，3.11 ，8.06 ，3.69 ，3.69 ，3.50 ，2.72 ，8.50 ，1.77 ，2.53 ，3.69 ，6.14 ，10.99 ，6.35 ，30.13 ，3.89 ，2.91 ，1.96 ，6.98 ，2.91 ，5.93 ，5.93 ，7.41 ，4.29 ，4.69 ，5.72 ，6.77 ，4.49 ，3.11 ，2.91 ，4.29 ，7.63 ，5.93 ，2.34 ，2.91 ，5.72

28、，2.34 ，6.56 ，4.69 ，6.35 ，5.10 ，4.69 ，3.50 ，4.69 ，4.49 ，3.30 ，2.91 ，4.09 ，5.72 ，4.90 ，4.90 ，4.09 ，5.93 ，2.91 ，2.72 ，2.34 ，2.53 ，3.89 ，2.34 ，2.34 ，1.96 ，2.72 ，3.11 ，3.50 ，2.72 ，1.77 ，3.50 ，2.53 ，9.17 ，7.41 ，5.72 ，8.06 ，5.72 ，8.72 ，6.77 ，6.56 ，6.77 ，6.98 ，6.14 ，6.69 ，8.00 ，6.69 ，8.23 ，9.35 ，8.90 ，3.77

29、，5.41 ，7.78 ，5.62 ，5.41 ，4.58 ，6.91 ，5.00 ，5.62 ，6.91 7.78 ，6.26 ，7.56 ，4.79 ，2.77 ，6.26 ，7.34 ，4.17 ，5.00 ，7.56 ，5.62 ，6.05 ，4.79 ，23.72 ，6.47 ，4.17 ，5.00 ，4.17 ，6.26 ，5.00 ，1.61 ，4.58 ，5.41 ，3.57 ，2.00 ，3.37 ，6.47 ，3.37 ，6.91 ，2.38 ，7.56 ，8.67 ，4.79 ，5.41 ，5.83 ，5.62 ，4.38 ，2.77 ，5.41 ，5.62 ，6.69

30、，7.56 ，2.77 ，1.61 ，5.62 ，2.97 ，4.58 ，5.20 ，3.17 ，5.00 ，7.56 ，5.20 ，6.05 ，2.77 ，2.38 ，5.20 ，8.67 ，5.41 ，7.12 ，4.58 ，8.23 ，4.58 ，6.05 ，8.23 ，4.17 ，3.97 ，10.74 ，3.77 ，10.27 ，5.20 ，6.47 ，6.47 ，8.90 ，3.37 ，6.69 ，4.17 ，8.23 ，4.38 ，5.20 ，5.20 ，5.41 ，5.00 ，9.58 ，3.17 ，3.77 ，9.13 ，7.34 ，5.62 ，4.79 ，5.20 ，3.3

31、7 ，7.34 ，3.57 ，4.38 ，5.41 ，4.38 ，5.62 ，6.05 ，6.26 ，5.20 ，6.91 ，4.58 ，8.67 ，6.47 ，7.12 ，3.77 ，5.41 ，10.97 ，9.81 ，8.23 ，5.41 ，2.77 ，7.78 ，6.47 ，5.62 ，3.77 ，3.57 ，4.38 ，6.91 ，7.12 ，4.38 ，3.97 ，8.00 ，3.57 ，3.37 ，3.17 ，2.97 ，7.12 ，5.62 ，5.83 ，3.97 ，3.17 ，2.77 ，3.17 ，2.97 ，3.17 ，1.80 ，3.17 ，3.77 ，4.79 ，3.

32、57 ，6.47 ，7.34 ，8.23 ，10.74 ，11.68 ，7.34 ，6.05 ，5.41 ，6.26 ，6.47 ，6.47 ，4.79 ，7.56 ，9.35 ;cx=0:100:30000;cy=0:100:30000;cAs=griddata(x,y,As,cx,cy','cubic');meshz(cx,cy,cAs)x=74，1373，1321，0，1049，1647，2883，2383，2708，2933，4233，4043，2427，3526，5062，4777，5868，6534，5481，4592，2486，3299，3573，4741

33、，5375，5635，5394，5291，4742，4948，5567，7004，7304，7048，8180，9328，9090，8049，8077，8017，6869，7056，7747，8457，9460，9062，9319，10631，10685，10643，11702，11730，11482，10700，10630，11678，11902，13244，12746，12855，13797，14325，15467，12442，13093，13920，14844，16569，16387，16061，15658，14298，14177，15092，12778，17044，17087，1707

34、5，17962，18413，19007，18738，17814，18134，17198，17144，18393，19767，21006，21091，22846，23664，22304，21418，21439，20554，20211，21072，20215，18993，19968，21766，22674，22535，25221，26453，26416，27816，25361，24065，25998，27177，26424，26073，24631，24702，25461，24813，26086，26015，27700，27696，27346，26591，27823，27232，24580，2415

35、3，22965，23198，24685，28654，24003，21684，22193，17079，15255，15007，3518，3469，3762，3927，4153，3267，4684，5495，5664，5541，5451，4020，4026，5101，5438，5382，5314，5503，5636，6605，7093，7100，6837，7906，8045，8394，8403，8079，9663，9469，9178，9095，10225，10210，10340，11557，11415，11649，12734，12696，12400，12591，13765，13694，13855，

36、14862，14896，15387，15810，16032，15801，15087，16872，17734，16823，17008，17203，17005，16947，16301，17904，18303，18438，18556，18954，18012，19072，20282，21475，21450，20261，19569，19411，19501，20582，19909，21018，22176，23359，23238，22624，21703，5006，5734，6395，7405，8446，7612，7912，8866，9296，9475，9212，8629，7776，8622，9237，830

37、7，7106，6423，7458，8904，10547，10398，10395，11529，11563，11646，12641，14000，14207，14065，12734，12727，14173，15467，15140，15198，15248，16428，16289，16267，16440，16440，15412，14269，13277，13175，12153，11958，10800，10022，9333，9277，11121，10856，12644，12625，9036，10599，12632，14405，14074，14262，14624，16629，18470，20591，20983

38、，20217，19041，18906，18467，17414，15748，15517，16607，15952，22605，23146，22046，23785，25981，27380，25021，23325，26852，17981，14482，14318，10352，9095，10510，13954，10142，17765，6924，4678，6182，5985，7653;y=781，731，1791，1787，2127，2728，3617，3692，2295，1767，895，1895，3971，4357，4339，4897，4904，5641，6004，4603，5999，6018，6213

39、，6434，8643，7965，8631，7349，7293，7293，6782，6226，5230，4600，4496，4311，5365，5439，6401，7210，7286，8348，8260，8991，8311，7639，6799，6472，5528，4472，4480，5532，6354，8184，8774，8618，7709，7056，8450，8945，9621，8666，8658，4329，4339，5354，5519，6055，6609，7352，7594，7418，6684，6936，5799，10691，11933，12924，12823，11721，11488，109

40、21，10707，10046，9810，9081，9183，8810，8819，9482，9149，9790，10527，10721，11383，11228，10774，10404，9951，12371，12961，12348，12173，11293，5795，5577，6508，5581，6423，7353，7032，7771，8639，8807，9422，9522，9834，10799，11094，12078，11609，11621，13331，13715，14737，14482，13319，12450，13535，13523，14278，8755，15286，13101，12185，58

41、94，5110，5535，2571，2308，2170，2110，2299，793，1364，1205，1653，2093，2757，2990，3913，4080，3994，3012，2060，1127，133，374，1381，2449，3490，3978，3052，2035，1075，0，1288，2286，3299，3975，3821，2789，1764，1581，2585，3515，4015，3024，2060，1063，1353，2357，3345，2524，1603，729，2307，3061，3966，3512，2798，3629，4207，4775，6218，7212，7487

42、，8299，8287，7385，6539，5588，4874，4414，8519，8590，8540，7555，7586，7348，6934，6091，6548，5300，5764，5492，5325，6502，4818，6591，8846，9659，10443，10981，11200，11938，12840，13143，13102，12000，11305，12086，10613，10638，9872，9726，9467，8831，8920，8868，9591，10360，11203，11243，10298，9381，9560，8970，9980，10987，10344，7691，11941，

43、12080，11101，10100，9106，9069，10072，11058，12068，13232，12982，12877，13204，12238，12336，13313，13282，12204，14631，16148，16432，14727，14943，16259，17538，17980，17949，18032，16516，15129，14004，14481，14411，13549，15862，17642，15769，16346，17001，15476，15728，17034，17365，18397，14301，15382，17634，17643，18051，18202，16290，16

44、701，16114，18449，12692，13569，17133，16414，15314，5615，1662，3561，5696，3765，2005，2567，1952;Cd=153.80 ，146.20 ，439.20 ，223.90 ，525.20 ，1092.90 ，269.80 ，1066.20 ，1123.90 ，267.10 ，201.40 ，287.00 ，193.70 ，359.50 ，516.40 ，1044.50 ，445.40 ，347.90 ，345.70 ，614.00 ，257.20 ，1213.50 ，325.80 ，212.10 ，90.50 ，583.40

45、，366.40 ，323.90 ，424.50 ，630.00 ，635.30 ，463.40 ，532.00 ，778.70 ，754.80 ，396.30 ，687.80 ，526.00 ，449.10 852.70 ，459.00 ，337.30 ，568.10 ，599.00 ，635.50 ，600.70 ，567.60 ，228.50 ，568.60 ，214.70 ，151.90 ，343.30 ，293.90 ，312.90 ，315.90 ，416.30 ，374.00 ，344.30 ，252.90 ，503.40 ，303.80 ，127.00 ，265.00 ，278.

46、90 ，751.20 ，361.30 ，488.00 ，227.00 ，347.40 ，136.00 ，327.10 ，113.10 ，270.50 ，160.30 ，305.50 ，70.90 ，261.30 ，119.80 ，468.80 ，150.70 ，100.30 ，109.80 ，91.80 ，743.90 ，416.80 ，369.80 ，194.00 ，50.10 ，198.80 ，886.60 ，128.90 ，114.30 ，232.90 ，132.80 ，1619.80 ，282.50 ，180.60 ，386.60 ，345.00 ，95.60 ，87.10 ，203.

47、00 ，353.00 ，233.20 ，174.30 ，87.60 ，245.60 ，167.90 ，111.40 ，94.60 ，111.60 ，85.50 ，169.50 ，138.00 ，131.40 ，41.00 ，129.70 ，148.00 ，108.10 ，132.50 ，162.90 ，88.80 ，206.70 ，121.80 ，231.10 ，128.60 ，126.80 ，79.70 ，294.10 ，238.30 ，109.30 ，69.40 ，155.70 ，109.30 ，205.50 ，53.20 ，40.00 ，82.80 ，222.40 ，164.80 ，19

48、4.00 ，310.10 ，1024.90 ，492.80 ，202.30 ，97.20 ，208.60 ，143.40 ，158.10 ，332.60 ，401.80 ，756.40 ，407.50 ，307.30 ，242.10 ，178.90 ，315.50 ，134.60 ，235.60 ，203.80 ，568.50 ，506.50 ，880.00 ，250.30 ，161.90 ，249.40 ，248.20 ，156.30 ，120.30 ，429.20 ，205.10 ，113.80 ，221.30 ，580.10 ，171.00 ，365.00 ，289.70 ，452.10

49、 ，847.60 ，310.70 ，317.50 ，357.50 ，387.20 ，196.50 ，295.80 ，129.00 ，204.90 ，174.80 ，197.00 ，106.70 ，738.00 ，788.70 ，133.10 ，149.50 ，206.60 ，78.90 ，178.60 ，250.10 ，89.20 ，453.50 ，92.40 ，198.70 ，260.10 ，101.10 ，146.50 ，66.50 ，119.00 ，434.80 ，188.00 ，77.50 ，130.80 ，134.50 ，101.80 ，232.00 ，121.10 ，238.70

50、，151.20 ，202.60 ，351.60 ，173.80 ，423.30 ，645.30 ，967.70 ，445.70 ，420.60 ，257.40 ，248.70 ，350.30 ，109.90 ，306.20 ，306.50 ，515.60 ，175.90 ，417.80 ，368.20 ，393.00 ，117.50 ，379.10 ，320.30 ，613.40 ，220.50 ，229.30 ，334.30 ，284.30 ，230.10 ，551.10 ，201.00 ，203.80 ，355.20 ，133.70 ，95.60 ，352.60 ，556.20 ，113.

51、70 ，240.50 ，120.30 ，514.10 ，246.30 ，158.00 ，211.30 ，236.50 ，193.00 ，169.70 ，320.10 ，180.20 ，351.80 ，245.70 ，86.80 ，367.80 ，467.10 ，364.70 ，248.50 ，171.80 ，409.90 ，302.50 ，236.20 ，114.50 ，165.20 ，380.40 ，398.40 ，268.60 ，126.50 ，290.30 ，228.40 ，305.50 ，407.60 ，96.60 ，185.50 ，288.70 ，90.40 ，285.40 ，100

52、.00 ，306.90 ，319.50 ，100.10 ，218.60 ，239.80 ，156.80 ，281.10 ，142.50 ，195.50 ，153.90 ，104.20 ，72.10 ，190.80 ，282.90 ，149.00 ，121.30 ，479.20 ，870.50 ，279.00 ，162.00 ，907.00 ，132.90 ，197.00 ，100.70 ，119.10 ，63.50 ，156.00 ;cx=0:100:30000;cy=0:100:30000;cCd=griddata(x,y,Cd,cx,cy','cubic');mes

53、hz(cx,cy,cCd)x=74，1373，1321，0，1049，1647，2883，2383，2708，2933，4233，4043，2427，3526，5062，4777，5868，6534，5481，4592，2486，3299，3573，4741，5375，5635，5394，5291，4742，4948，5567，7004，7304，7048，8180，9328，9090，8049，8077，8017，6869，7056，7747，8457，9460，9062，9319，10631，10685，10643，11702，11730，11482，10700，10630，11678，1

54、1902，13244，12746，12855，13797，14325，15467，12442，13093，13920，14844，16569，16387，16061，15658，14298，14177，15092，12778，17044，17087，17075，17962，18413，19007，18738，17814，18134，17198，17144，18393，19767，21006，21091，22846，23664，22304，21418，21439，20554，20211，21072，20215，18993，19968，21766，22674，22535，25221，26453，2

55、6416，27816，25361，24065，25998，27177，26424，26073，24631，24702，25461，24813，26086，26015，27700，27696，27346，26591，27823，27232，24580，24153，22965，23198，24685，28654，24003，21684，22193，17079，15255，15007，3518，3469，3762，3927，4153，3267，4684，5495，5664，5541，5451，4020，4026，5101，5438，5382，5314，5503，5636，6605，7093，7100

56、，6837，7906，8045，8394，8403，8079，9663，9469，9178，9095，10225，10210，10340，11557，11415，11649，12734，12696，12400，12591，13765，13694，13855，14862，14896，15387，15810，16032，15801，15087，16872，17734，16823，17008，17203，17005，16947，16301，17904，18303，18438，18556，18954，18012，19072，20282，21475，21450，20261，19569，19411，195

57、01，20582，19909，21018，22176，23359，23238，22624，21703，5006，5734，6395，7405，8446，7612，7912，8866，9296，9475，9212，8629，7776，8622，9237，8307，7106，6423，7458，8904，10547，10398，10395，11529，11563，11646，12641，14000，14207，14065，12734，12727，14173，15467，15140，15198，15248，16428，16289，16267，16440，16440，15412，14269，13277，13175，12153，11958，10800，10022，9333，9277，11121，10856，12644，12625，9036，10599，12632，14405，14074，14262，146