数学建模国赛A题Word文档下载推荐.docx

1、1.假设采样数据合理，能真实的反应当地的情况。2.假设各重金属在传播过程中不与其他物质反应，也不会有衰减。3.假设污染严重的局部地区仅存在唯一的污染源。4.重金属的传播互不影响，且污染源的排放量不变。四符号说明 ： 第i项重金属的污染指数 第i项重金属浓度 第i项金属的背景值浓度 : 第n个观测点 污染源中心 污染源排污总量: 第t年污染源总排污量 观测点位置为x的污染物五模型建立与求解5.1问题一模型的建立与求解我们首先用Surfer软件画出该城区简单的空间分布图，x，y表示采样点坐标位置，z轴表示该采样点对应的海拔高度，所得图形如下：图1 空间分布图 图中，向上为北，向右为东，且等高线越密

2、集表示地势越高。由图可看出，东北方向等高线密集，因此可知该方向为山区；而西面地势较平坦，为平原。使用surfer分别绘制重金属在该城区的分布图如下： （1） As分布图 （2） Cd分布图 （3） Cr分布图 （4） Cu分布图 （5） Hg分布图 （6） Ni分布图 （7） Pb分布图 （8） Zn分布图图2 各重金属元素的空间分布图由图可知，As元素主要分布在三块区域：横坐标为4000-5000,纵坐标为6500-7500的区域；横坐标为12500-13500，纵坐标为2500-3500的区域；横坐标为17500-18500，纵坐标为9500-10500的区域，其他地域相对分布较少。Cd元

3、素主要分布在该城区的西南部、中部，其他地方分布很少。Cr元素总体分布较少，主要集中在以坐标（4500,5000）为圆心的1000米范围内。Cu元素的分布主要集中在西南面。Hg元素主要集中在三个区域，分别位于横坐标，但大部分地区都不含此元素。Ni元素主要分布在横坐标4000-5000，纵坐标为5000-6000的区域，其他区域分布较少。Pb元素主要分布在西南部，其他地区相对较少。Zn元素主要分布在中部及西南部，其中西南部较集中。因此，总的来说，西南面污染较为严重。单项污染指数公式： （1）分析该城区不同区域重金属的污染指数，如下表：表1 各区域重金属元素单项污染指数As Cd Cr Cu Hg

4、Ni Pb 生活区1.7912.3782.2453.4883.0301.5022.265工业区2.0143.0241.7239.66218.351.6113.001山区1.1231.1721.2571.3121.1701.2561.179交通区1.5862.7421.8614.67612.7281.4302.033公园绿地区1.7402.1581.4082.2873.2851.2431.958建立内梅罗环境污染指数模型：（2）其中，代入数据可得五个功能区的内梅罗环境污染指数如下表：表2梅罗污染指数功能区综合污染指数1.8026683.4481871.1187911.9096121.631696

5、表3土壤内染指数评价标准等级一级二级三级四级五级内梅罗污染指数等级p清洁警戒线轻度污染中度污染重污染结合内梅罗指数评价标准可知，生活区、山区、交通区及公园绿地区受到了轻度污染，工业区受到了重度污染。5.2 问题二模型的建立与求解为使污染指数更明显，将表1转化为条形图，如下：图3各区域重金属元素单项污染指数条形图由图可知各种重金属对各区域的污染程度都不尽相同，有着较大的差异。重金属污染主要来源于工业污染，其次是交通污染和生活污染重金属污染大多通过工业三废，即废水、废气、废渣进入环境，污染土地；交通污染主要由于汽油中含有重金属元素，这些元素通过汽车尾气进入空气中，对环境造成污染；而生活污染主要是生

6、活垃圾：废电池、废弃金属制品、废弃化妆品、塑料制品等。根据以上分析我们可以得出：人类活动越密集的地方，重金属污染越严重，因此可以初步确定是人类的各项活动造成了重金属污染。表4 各行业主要污染物类型工业污染源采矿冶金等行业三废排放：Hg、Cr、Cd、As、Pb、Zn交通污染源汽车尾气排放：Pb、Zn、Cd、Cr、Cu生活污染源生活废弃物：Cr、Cu、Zn、Pb对照表1知，As,Cd,Cr,Ni含量较为接近背景值，而Cu,Hg,Pb,Zn含量在人类活动频繁的区域明显增加，其中Hg,Cu含量增加明显，工业区Hg平均超值18.4倍，Cu超标9.7倍，说明此城市存在Hg,Cu排放严重超标的企业，并且其排

7、放的污染物已经对该城市造成了严重的影响，可推测该城市可能存在大型化工厂或电池制造场等高等高排放量企业或工厂。各土壤重金属之间存在某种依存关系，而这种依存关系能帮助揭示重金属元素的污染来源，我们选用受人类活动影响较大的Cu,Hg,Zn,Pb为例，对它们进行相关性分析。由于功能区三和功能区五重金属元素含量超标不多，因此我们选用功能区一、二、四进行分析。使用SPSS分别对功能区一，功能区二和功能区四的Hg,Cu,Pb,Zn做相关性分析，得到下表：表5 Hg,Cu,Pb,Zn相关性相关系数RHg-CuHg-PbHg-ZnCu-PbCu-ZnPb-Zn功能区一0.9910.9900.9830.9960.

8、9840.982功能区二06120.5900.6700.6220.739功能区四0. 0330.2670.1190.5070.4330.480分析上表可知：第一：Cu，Zn，Pb间的相关性较强，其中Zn,Pb自相关性较强，有资料表明受交通污染源控制时，交通道路断面土壤重金属污染程度取决于其距离交通道路的远近。离公路越远，受交通污染源的影响越小，土壤Pb、Zn污染越轻，递减趋势明显。土壤中Zn和Pb含量随道路远近的变化趋势具有较强的一致性，说明土壤Pb和Zn源于相同的污染源为交通污染源。第二：不同功能区土壤重金属间依存关系各不相同，反映各地区重金属来源的途径存在差异。在工业区，部分土壤重金属间依

9、附关系不明显，而部分金属间的依附关系显著，这反映该功能区受到复合型污染源（即同时排放不同重金属元素的污染源）影响，而公园绿地区、山区重金属污染相对小。因此，该城市工业区，交通区，生活区污染均较为严重，其中重金属元素Hg含量最为突出，可推测该地区存在大型或数量较多的Hg严重超标的企业；Cu，Zn，Pb等重金属元素浓度也严重超标，由于Cu ,Zn ,Pb在各地区含量均存在很强的相关性，故推测Cu ,Zn ,Pb等重金属元素来源于超标排放的企业或工厂。不同功能区，反映出污染源的性质不同，如工业区反映重金属污染源的多样性，生活区和交通区反映污染源单一，分别为生活废弃物和汽车尾气排放。5.3 问题三模型

10、的建立与求解考虑到重金属在水体中溶解态、悬浮态和底泥态的物理运动物质交换，吸附和解吸物质交换，以及沉降和再悬浮物质交换过程复杂，而且对城市的河流，悬浮物，底泥无相关信息。故从总体角度将重金属输送迁移基本方程简化为平面扩散模型。因为重金属的传播以污染源为中心，向四周散开，因此浓度最大的地方就是污染源，将问题转化为求局部最大值。假设在局部足够小的区域内对于某类重金属存在唯一的重金属污染源Q，此污染源对周围地区造成了重金属污染，在该地区寻找n个观测点，其重金属污染物含量为5.3.1 一维污染源模型对于这一特定的重金属污染源Q,假设污染源的总排放量为对污染物的直线扩散方式进行分析，记坐标轴为x的观测点

11、的污染物含量为Q（x），若受体距离污染源越远，则污染源对其作用效果越不明显，即传播系数为随距离递减的曲线，则在的观测点，其污染物的一维分布形式满足：（3）为重金属的传播系数，为重金属污染物的传播范围，因为重金属污染在传播过程中，性质稳定，且具有富集的特点,在空间传播速度较慢，故一般认为为一常数，由重金属污染物的性质决定。当时有：（4）即：（5）分离变量得：（6）因此模型可以简化为:（7）发现函数模型图形呈倒扣钟形，且近似高斯分布曲线，其中当时（为污染源的位置），函数取得最大值，5.3.2 二维污染源传播模型重金属污染物分布在平面上，依据上述分析，对一维污染物模型进行推广，忽略空间的各向异性，可

12、知重金属在空间上的分布为二维高斯分布形式，其形状为一倒扣钟形，其扩散模型为：（8）表示污染物的空间分布状态取n个样本点,对函数进行拟合。因为局部地区的污染源位置大致位于极大值的所处周围区域，故选取极大值附近的点作为拟合函数的样本点具有实际意义。以As为例进行说明：图4 As污染分布图如图所示，存在三个等浓度线密集的区域，表明空间分布明显存在三个污染源地区，分别位于,；和的三个区域，故在该范围内寻找极值点。选取周围10个点作为样本点，算法实现见附录，在VC6.0中运行程序，得出以下结果：表6 As的污染源样本点 （1）污染源一样本点 （2）污染源二样本点编号As元素含量X轴Y轴8430.1318

13、13410046836.351781410707871.96183939183853.891719898108210.991873810921762.72170441069186291171449081816.141900711488803.6918413117211995.83179048287As元素浓度17823.721269630241774.791273440151796,471240020601766.051164935151824.171369423571836.26138553345641242243291755.62114152585653.11130934339515.31

14、117024480（3） 污染源三样本点对模型两边取对数，得：（9）对该模型进行拟合（代码见附录），结果如下：污染源1:（10）污染源2：（11）污染源3：（12）得到三个污染源分别为：污染源一：（15598,9345）；污染源二：（14766,8021）；污染源三：（1163.4,770.1）。经查表分析知污染源附近均存在污染特别严重的区域，故模型较为可信。同理，对其他重金属进行分析，所得污染源如下表所示：表7 其他几种重金属污染源坐标污染物污染源Cd（2723,2315.4），（21442,11381）Cr（3381.2 ,6001.8）Cu（2459.2, 3632.7）Hg（15247

15、,9111），（2459.2 ,3632.7） （-20777 , 93031）Ni（424.38,869.41），（22193,12185）Pb（4778.7 ,4891.5）Zn（8182.7,4493.0），（13800 , 9623）由表可知，部分金属的污染源坐标为负值，说明该金属的污染源不在该城区内，而应在该城的城区外，城内土壤重金属含量超标是因为受到了城区外污染源的影响。5.4 问题四模型的建立与求解问题3中的重金属污染源传播模型以及污染源位置的确定模型：优点：通过污染物的传播特点，建立二维污染物空间分布模型，准确反映了污染物的平面分布状态，考虑到污染源的位于污染严重的附近区域，仅

16、在污染严重的区域搜索样本点，极大地简化算法，且能准确反映污染源的位置。采用多元非线性拟合二维高斯分布函数，确定局部极值所对应的坐标轴，确定污染源所在地。缺点：二维的污染物扩散模型仅反映了空间的污染物积累对于周围地区污染扩散分布的影响，却不能准确反映重金属污染传播的时空分布关系，而污染物扩散模型是对污染物扩散的模拟，是对土壤重金属污染物随时间和空间迁移转化规律的描述，即用于描述土壤中重金属污染物质与时间、空间定量关系的数学方程。5.4.1 一维时空扩散模型由重金属污染源的空间分布，确定污染源的大致区域，然后收集该地区工业，交通业发展的状况，引进时间因素t，对二维的污染物扩散模型进行修正。污染物排

17、放的扩散与在空间的分布形式与污染物排放的时间存在相互制约的关系，当污染物排放时间较短，污染物集中分布在u位置及其周围地区；当污染物排放时间较长时，污染物相对均匀分布在u位置及其周围地区。二维扩散模型是对污染物在空间的富集状态的描述，对于某一年份污染物排量在空间上的分布也服从二维扩散模型，为了体现不同年份的影响效果，由高斯分布模型知道钟形的形状即污染物的空间分布状况由参数决定。引入时间的影响，即第t年份污染物排放总量对空间分布的影响模型如下：（13）已知随t的增大而增大，且当时有,可考虑线性模型或指数模型：采用上述搜索算法，在局部区域内搜索n个观测点，因为重金属性质稳定，不易发生化学变化等改变物

18、质属性的，且传播过程易出现富集，人工治理难度大，所以各观测点的污染物的含量为污染源多年排放量的线性累积和，假设该地区污染企业存在m年，故得到：（14）5.4.2 二维扩散模型假设局部地区的污染源坐标为,对上述一维模型推广应用得到重金属的二维空间分布模型为:（15）模型综合考虑空间和时间对于污染物扩散的影响，更加精确描述重金属污染物的传播方式，为研究城市地质环境的演变模式提供理论基础。六模型评价1.考虑到污染物扩散的累积效果，使结果更加精确，更加贴近实际，得出的结果对重金属污染防治具有现实指导意义。2.问题一采用Surfer软件绘图，能清晰的反映出不同地区的浓度。3.问题二从相关性分析污染源性质

19、，能清晰的反映出主要污染及其污染源源头。4.假设局部地区存在单一污染源很大程度上简化了问题。1.重金属的传播模型没有体现基于功能区土壤，气候，水流等对金属污染物传播方式的影响。2. 模型忽视了污染源分布的随机性。3. 重金属污染的非确定性因素很多，难以考虑周全。七模型的改进与推广在现有模型的基础上搜集该城市各功能区的土壤，气候，水流等信息，分析传播媒介的主要性质，确定传播系数等数值，建立基于土壤，气候，水流的信息分析的对流扩散模型。以此作为重金属污染的传播模型，能更加准确研究城市地质环境的演变模式。改进后的模型不但对重金属的污染源寻址提供较大的帮助，且在对其他污染物的扩散与治理的相关研究也具有

20、巨大的参考意义。八. 参考文献【1】石东伟、陈娜，扩散模型在确定污染源位置中的应用，Journal of Henan Institute of Science and Technology,2012【2】刘衍君、汤庆新，基于地质积累与内梅罗指数的耕地重金属污染研究，中国农学报，2009【3】中国环保信息网，城市土壤的重金属污染及生物治理，【4】李亚松、张兆吉、费宇红、王昭，内梅罗指数评价法的修正及其应用，中国地质科学院水文地质环境地质研究所，2009年附录一：VC6.0 找污染源所在，及十个样本点/程序计算各功能区各类金属的单因子指数以为各功能区的内梅罗综合系数/搜索样本点/程序在VC6.0运行通过 #include#includefstreamstdlib.hiomanipcmathusing namespace std;struct Area /各污染物种类见命名 int ID;