2011数学建模A题---城市表层土壤重金属污染分析

32011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：重庆交通大学参赛队员(打印并签名)：1.陈训教2.范雷3.陈芮指导教师或指导教师组负责人(打印并签名)：胡小虎日期：2011年9月12日赛区评阅编号（由赛区组委会评阅前进行编号）：42011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：5城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染做出了详细的分析，对于本题中所提出的问题一，我们利用MATLAB软件对所给的数值进行空间作图，然后分别作出了八种重金属元素的空间分布特征，然后，我们利用综合指数（内梅罗指数）评价的方法，对五个区域进行了综合评价，得出结果令人满意。对于问题二，我们根据第一问和题目所给的数据进行综合分析，得出了重金属污染的主要原因来自于交通区含铅为主的大量排放，和工业区污水的大量排放等等。对于问题三，我们通过对问题一中的八张重金属元素空间分布的图可以看出，发现大多数金属都呈中心发散性传播，同时经过分析，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析和拟合检验，发现在参数变化在10%左右，模型的稳定性良好。最后我们全面分析了模型的优缺点，，最后可以用MATLAB软件得出相应的结果。为更好地研究城市地质环境的演变模式，测定污染源范围还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，地下水流动方向以及每年的生物降解量，降雨量对重金属元素扩散的影响。一但有污染证据，我们可以在该污染源附近沿地下水流动方向设定更多采样点，由此，我们可以构造一个三维公式来计算污染物质浓度的浮动就可以模拟三维空间内的重金属分布影响。关键字：表层土壤重金属污染MATLAB内梅罗指数偏微分方程稳定性检验灵敏性分析地质演变生物降解量6一、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点，对城区重金属污染分析以及传播越来越有其必要性。对于本题中所提出的问题一，我们利用MATLAB软件对所给的数值进行空间作图，分别作出了八种重金属元素的空间分布特征图，然后，我们运用综合指数（内梅罗指数）评价的方法，对五个区域进行了综合评价，得出结果令人满意。对于问题二，我们根据第一问和题目所给的数据进行综合分析，得出了重金属污染的主要原因来自于交通区大量排放的含铅为主污染物，和工业区污水的大量排放等等。对问题三，通过分析第一问中所给定各种元素空间分布规律，和查阅大量资料，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析，发现在参数变化在10%左右，模型的稳定性良好。为更好地研究城市地质环境的演变模式，测定污染源范围还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，地下水流动方向以及每年的生物降解量，降雨量对重金属元素扩散的影响。一旦有污染证据，我们可以在该污染源附近沿地下水流动方向设定更多采样点，由此，我们可以构造一个三维公式来计算污染物质浓度的浮动就可以模拟三维空间内的重金属分布影响。7二、问题分析针对本题所提出的问题一，我们利用题给的数据运用MATLAB软件作出相应的八中重金属元素空间分布图，同时考虑到分析和评价城区累不同金属的污染度，我们根据数据应用内梅罗指数综合评价城区内不同区域的重金属污染度。得到的结果令人满意。针对问题二，通过对第一问和原始数据的分析，我们得出了重金属污染的主要原因还是交通区汽车含铅气体的大量排放，工业区不达标的污水大量排放等等。对问题三，通过分析第一问中所给定各种元素空间分布规律，和查阅大量资料，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析，发现在参数变化在10%左右，模型的稳定性良好。证明了我们模型的正确性。对于问题四，三、模型假设假设一：问题中附件给出的原始数据真实，有效。假设二：（0～10cm）的地表是由土和沙砾所构成。假设三：所有的土和蓄水层的性质在浸透到水层和未浸透水层都是均匀的。假设四：稳定、均匀内的水流只发生在整个未浸透水层的垂直方向，以及在浸透水层水平（纵向）平面中。沿地表水流动方向。假设五：物理过程起着重要作用，在此模型中只考虑物理过程（扩散），不考虑在其过程中的化学反应。假设六：所有的污染源均为点源。8四、符号说明iP：单项内梅罗指标数iaveP：iP的平均值maxiP：iP的最大值P：内梅罗综合指数IC：污染物实测值aX：污染金属背景值的最小值cX：污染金属背景值的平均值pX：污染金属背景值最大值la：水平扩散系数（m）ta：垂直扩散系数（m）C：污染物质的浓度（litermg）bC：背景浓度（litermg）0C:污染物源处的浓度（litermg）D：渗透系数（sm2）K：水力传导性（2ftdaygal）I：水力梯度m：污染物质的排放率（daymg）9n：有效多孔性q:污染物质的排放率（daymg）dR：阻滞因子S:复合参数0t：污染开始的时间（yr）dV：地下水的流速（dayft）W：hantush函数),(00yx：污染源的坐标10五、模型的建立与求解5.1.1问题一的模型建立从附件-1中得到的数据，运用MATLAB软件，绘制了下面的三维地表图形，图形中的黑点就是取样的地点。（z代表的是海拔）。1112131415161718下面建立该城区内不同区域重金属的污染程度，为了更好地评价该城区不同区域的重金属污染程度，我们引入目前比较广泛应用的评价指标——内梅罗指数，进行综合评定，单项污染指数法能够比较直观地反映环境中各项污染指标的情况；内梅罗综合指数法不仅考虑到了所有评价因子单项污染程度的平均水平，而且还考虑到了最大污染指数，因此能够更为科学、综合的反映评价区域内总体环境质量状况。评价公式如下：单项指标数iP：19picppiipiccpaiiciaacaiiaiiiiXCXXXCPXCXXXXCPXCXXXXCPXCXCP)/()(3)/()(2)/()(1/（1）其中：iC为污染物实测值，pcaXXX,,是依据附件给出的平均值，偏差，以及范围所定出的界限值。综合污染指数（P）采用内梅罗综合指数法：222maxiaveiPPP（2）其中：maxiP为单项指标数iP的最大值，iaveP为单项指标数iP的平均值。给定污染综合指数等级的划分评定表，见表一：表一：土壤污染评价分级标准等级划分P污染等级污染水平1P0.7安全清洁217.0P警戒线尚清洁321P轻度污染土壤轻度污染432P中度污染土壤中度污染53P重度污染土壤重度污染5.1.2问题一的模型求解依据原始数据和公式（1）（2），利用MATLAB编程求得八种重金属元素的不同区域的单项内梅罗指数，见表二：表二：八种重金属元素的不同区域的单项内梅罗指数As(μg/g)Cd(ng/g)Cr(μg/g)Cu(μg/g)Hg(ng/g)Ni(μg/g)Pb(μg/g)Zn(μg/g)生活区0.4180.9330670.7668891.970.6202670.45851.1586512.123367工业区0.4833332.1330140.5934442.2215432.28530.495250.2699532.259767山区0.2693331.0154670.4328890.5773330.2730670.386251.0072560.862235交通区0.3806672.0857290.6452.0348862.122350.44051.1326982.142833公园绿地区0.4173331.8702670.4848891.00960.76660.382251.119581.602087为了能更好评出等级，我们求出各区域的内梅罗综合指数，也即求每个区域的八种重金属元素的P值的平均值，最终评价出污染等级，其果见表二，最后得出五20大区域的污染等级见（表三）表三：五大地区的污染等级城区不同区域区域内梅罗综合指数污染等级污染水平生活区1.718093轻度污染工业区1.9205133轻度污染山区0.8350922警戒线交通区1.79966263轻度污染公园绿地区1.4854143轻度污染5.2分析说明重金属污染的主要原因5.2．1该市表层土壤重金属含量基本状况分析将实际测得的该市五大区域土壤样品重金属含量与自然区土壤重金属含量背景值进行比较发现，五大区8大重金属含量平均值均高于背景值（表4）并且8大金属含量大部分样点超过自然区表层土壤均值，占总样点数的百分比最小为67.71%，最大为88.09%。这说明表土层8大重金属均有外源物质的进入，并有了一定的积累。表四5.2．2表层土壤重金属污染总体评价根据五大区域8种重金属内梅罗指数（表二）求出各元素的单项内梅罗平均指数（见表五），再对照土壤的等级评价标准（表一）对表层土壤污染进行评价，评价得出Cd、Cu、Hg、Zn这四种重金属指数达到了污染等级指数，污染级污染样点分别占79.94%、88.08%、67.71%、78.68%，均属于大面积轻度污染，其中Pb临近警戒线。这就可以看出，Cd、Cu、Hg、Zn是构成污染的的主要原因，需要加强控制并采取相应措施进行降排，而Pb需要预防性控制，其余的As、Cr、Ni均较为安全。重金属种类AsCdCrCuHgNiPbZn实测重金属含量均值5.68302.453.5155.06299.7117.2661.74201.21重金属含量背景值3.61303113.23512.33169超过均值样点数246255257281216240260251占样点总数百分比（%）0.77120.79940.805640.8808780.677120.75240.8150470.7868321表五：各元素的单项内梅罗平均指数元素种类AsCdCrCuHgNiPbZn平均指数0.39371.60750.584621.5626721.213520.43260.9376281.798065.2．3地区对四种轻度污染元素的影响经过对样点金属含量的详细分析，可以看出各种重金属在不同的区域的分不存在着很大的差异性。于是针对达到污染等级的四种重