2011数学建模城市表层土壤重金属污染分析.docVIP

城市表层土壤重金属污染分析模型 摘要 本文针对城市土壤表层土壤重金属浓度的分析，运用单因子污染指数法和内梅罗综合污染指数法对城区的污染程度进行综合评价，结合重金属污染物的传播特性建立扩散模型并确定污染源的位置。 对于问题一，采用单因子污染指数法和内梅罗综合污染指数，对城市表层土壤重金属污染进行综合评价。然后根据土地功能区的划分，分别对各区域进行综合评价，得到各个区域的污染程度：生活区、工业区和主干道路区为重度污染，生活区为中度污染，山区为轻度污染。利用atlab作图，画出以坐标x、y及重金属浓度大小的等高线图。可整个城中某一重金属元素的大致分布。工业三废引起的重金属污染汽车尾气的排放 现要求你们通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 模型假设 1.假设不考虑重金属的衰变； 2.假设污染源的重金属浓度不再增加； 3.假设不考虑风向对污染物传播的影响； 4.假设土壤重金属污染程度仅和这8种重金属有关； 5.假设采样点的数据较好的反映了该地区的污染物浓度； 6.假设每种重金属都是可以在外界因素下可以自由传播的，不存在重金属元素在一个地区堆积，只有传播能力的不同 符号说明 采样点i第j种金属的实测浓度(i=1,2…319;j=1.2…8) 第j种重金属的土壤环境质量标准临界值（j=1.2…8) 采样点i的第j种重金属的环境质量指数(i=1,2…319;j=1.2…8) 采样点i的第j种重金属污染指数的平均值(i=1,2…319;j=1.2…8) 采样点i的j种重金属污染指数中的最大值(i=1,2…319;j=1.2…8) 各采样点的重金属污染的综合指数 标准化后的采样点i第j种重金属的实测浓度(i=1,2…319;j=1.2…8) 第j种重金属浓度的平均值(j=1.2…8) 第j种为重金属浓度的标准差(j=1.2…8) D(x,y,z) 扩散系数 C(x,y,z,) 时刻污染物的浓度 C(x,y,z,) 时刻污染物的浓度 问题分析 对于问题一：要求8种主要重金属在该城区的空间分布，但由于采样点是离散的，可以通过插值得到8种重金属含量在xy平面上的等值线分布图。对于土壤的污染程度分析，可以单因子质量指数法得到所有样本点的各种金属污染指数。于要分析的是不同区域的污染程度，可以由内梅罗综合污染指数法计算所有样点的综合污染指数，但要考虑的是5个区域的污染度，根据功能区分成5个区。又可由统计分析各区域的平均值的污染程度和各污染等级占的百分含量，结合污染综合指数的等级标准得到各区域的污染程度。 对于问题二：则要求通过数据分析说明重金属污染的主要原因，即要分析重金属污染的主要来源和成因。可以通过因子分析法得到影响土地重金属污染的主因子，然后根据旋转后的载荷矩阵计算主因子得分，分析各重金属污染的主要原因。由于要说明重金属污染的主要原因，即最后综合考虑5个区域的污染程度，分析城区总体重金属污染的主要原因。 对于问题三：要分析重金属污染物的传播特征，通过查询资料可知重金属污染物的传播主要由高浓度向浓度，由高海拔向地海拔传播，即污染物的传播为扩散模型。而污染物的传播又主要通过大气和液体的渗透，即满足扩散定律。最后根据污染源在重金属高浓度集中的区域，结合各种重金属浓度的插值图进行计算，可得到各种重金属污染源的位置。 对于问题四：尽管问题三建立的扩散模型使得数据处理更加简单，但是没有充分考虑影响地质环境变化的因素。因此为了更好地研究城市地质环境的演变模式，可以收集大气重金属自然沉降、大气风向和污染物的形态的信息量。在有了这些信息量之后，通过对数据的分析拟合，可以构造一个三维曲面，并且结合问题三建立的模型可以更准确的确定污染源的位置。 模型的建立 问题1 由采样点的各种重金属分布的不同浓度，结合采样点的地理位置，对于污染程度的评价模型，可以根据单因子质量指数法得到各采样点的污染指数，单因子污染指数法是以土壤元素背景值为评价标准来评价重金属元素的污染程度，表达式为： （公式一） （设为土壤中重金属i的环境质量指数；为重金属i 的实际浓度, 为土壤重金属j环境质量标准临界值，即背景值中范围的最大临界值，式中1表示污染；≦1表示未污染，且越大，则污染越严重） 由上可知所有采样点8种重金属的环境质量指数，为