2022年数学建模A题国家一等奖 .pdf
《2022年数学建模A题国家一等奖 .pdf》由会员分享,可在线阅读,更多相关《2022年数学建模A题国家一等奖 .pdf(28页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料) ,必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。我们参赛选择的题号是(从A/B/C/D 中选择一项填写):A 我们的参赛报名号为(如果赛区设置报名
2、号的话):所属学校(请填写完整的全名) :参赛队员(打印并签名 ) :1. 2. 3. 指导教师或指导教师组负责人(打印并签名 ):日期: 2011 年 09 月 12 日赛区评阅编号(由赛区组委会评阅前进行编号):名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页,共 28 页 - - - - - - - - - 2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由
3、赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页,共 28 页 - - - - - - - - - 1 城市表层土壤重金属污染分析摘要本文针对城市表层土壤受重金属污染问题进行综合分析。首先运用Matlab软件求解8 种主要重金属元素在城区空间分布,土壤受污染主要原因,传播特征以及为今后如何更好研究地质演变问题分别建立了相应的数学模型,并对其求解结果作出了分析。针对问题 1,根据各种污染物浓度在不同区域内分布
4、的随机性,利用空间内插值法,以城市位置为平面“横纵向” ,污染物浓度为“竖向”建立三维空间模型,得到这8 种重金属元素的空间分布图,直观反映出污染物浓度的等值线,得到污染物浓度的分布规律和各种重金属元素的污染浓度范围区块。然后通过建立污染负荷指数法模型算出各区的)(PLI值,生活区)(PLI为 1.8336,工业区)(PLI为 2.1573,山区)(PLI为 1.0602,交通区)(PLI为 1.9209,公园绿地区)(PLI为 1.5780;结果表明工业区的污染程度最高为 2 级强污染,其他区为1 级中等污染。针对问题 2,根据问题1 中不同功能区重金属的污染程度与附件2 的数据,利用Pea
5、rson相关性分析和主成分分析法,建立了相关性模型与主成分模型,对应找出各区重金属污染的主要因子,得到该城区的重金属污染主要元素为:ugCdCHPb、,污染主要原因为冶炼、化物生产等工业废水、污泥,汽车尾气排放。针对问题 3,根据重金属污染物的圆型传播特征,利用物理热力学第二定理,分别建立了一维土壤迁移模型、沉降模型和最优解模型,确定了生活区有3 个污染源:分别为样本点 20、306、259;工业区的污染源为样本点261 ;山区的污染源为样本点62 ;交通区有 2 个污染源:分别为样本点245、292 ;公园绿地区的污染源为样本点315 。其结果与重金属元素在该城区的空间分布图作对比,吻合得很
6、好。针对问题 4,通过对问题 1,2,3所建模型优缺点的分析, 为更好地研究城市地质环境的演变模式,还应收集城区每年生活、工业等重要污染源的垃圾排放放量、每年的生物降解量、排污企业个数的增减值、PH值、河流所经区域等信息。根据这些信息建立时间序列模型,用此模型预测未来时间段该城区的土壤中重金属含量变化情况,结合问题1,2,3 所建模型可以更好地为今后如何使用专业降浓剂保证土壤良好性的用剂量提供参考依据。关键词: 重金属空间分布污染负荷指数主成份分析最优解名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - -
7、 - - - 第 3 页,共 28 页 - - - - - - - - - 2 一、问题重述随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据资料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关注的焦点。按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等,分别记为1 类区、 2 类区、5 类区,不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此,将所考察的城区划分为间距1 公里左右的网格子区域,按照每平方公里1 个采样点
8、对表层土( 010 厘米深度)进行取样、编号,并用GPS 记录采样点的位置。应用专门仪器测试分析,获得了每个样本所含的多种化学元素的浓度数据。另一方面,按照2 公里的间距在那些远离人群及工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值。附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息,附件2 列出了 8 种主要重金属元素在采样点处的浓度,附件3 列出了 8 种主要重金属元素的背景值。通过数学建模来完成以下任务:(1) 给出 8 种主要重金属元素在该城区的空间分布,并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析,说明重金属污染的主要原因。(3) 分析重金属污染物
9、的传播特征,由此建立模型,确定污染源的位置。(4) 分析你所建立模型的优缺点,为更好地研究城市地质环境的演变模式,还应收集什么信息?有了这些信息,如何建立模型解决问题?二、问题分析2.1 问题一的分析通过对附件 1 的数据进行处理求出该城区的城区状貌,并对该城区利用三维插值法求出 8 种主要重金属元素在该城区的空间分布。并对 8 种金属分别进行查找出各自范围区块。利用最高污染系数)(CF分别求出 8 种重金属的各自污染系数,然后利用污染负荷指数法)(PLI求出 319 个地点的重金属污染程度。根据这些地点的污染程度对5 个不同区域进去综合评价得出不同区域的污染程度。2.2 问题二的分析研究重金
10、属污染的主要原因既为研究重金属污染的主要来源,而同一活动可以产生许多污染物,同一种污染物的来源也不尽相同。因此利用Pearson相关分析和主成份分析方法解析重金属污染的元素和分析土壤重金属污染物的来源和类别。通过对附件 2、3的数据进行分析综合运用这些统计分析方法来分析5个不同区域的 8种重金属污染物的来源,并对 5 个不同区域查找相关性比较强的区域。最后综合5 个不同区域的来源分析出该城区的重金属污染的主要原因。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页,共 28 页
11、- - - - - - - - - 3 2.3 问题三的分析对于重金属的传播特征形式, 利用生物学基础对各种金属在该城区的传播形式进行假设分析出三种传播特征:圆型传播、条型传播、散点传播。结合前两问求出的结果和8 种主要金属元素浓度含量的空间分布图分别求出属于圆型传播、条型传播、散点传播的金属元素,由此可对这8 种重金属元素综合归类为那种传播特征。由于污染源在不同区域的位置对各个区域的影响有一个直接影响值,在污染源数目较多的场合下,要确定最优解的模型污染源是很困难的,而且在处理数据时还有考虑该地区的地势情况中的地下金属传播途径的迁移形式和大气中大气降沉形式,则综合考虑该城区的污染源就更困难。因
12、此通过对模型的研究,建立迁移模型、沉降模型对土壤金属流动进行分析,再结合最优解模型找出该城区5个不同区域的污染源。2.4 问题四的分析对于前三问建立的模型和结果进行分析查找出所建立的模型对该城区和不同区域的优缺点。由于该城区的城市布局和降水量,PH值等不知道,因此需要去完善模型并预测地质的演变,就得收集城区每年生活工业等重要污染源的垃圾排放放量以及每年的生物降解量,降雨量,PH值,河流所经区域等。通过这些数据建立起各重金属元素与时间有关的时间预测模型,由此可预测出未来时间段内的土壤重金属浓度,即知道该城区地质环境的演变模式。且根据不同功能区的不同安全浓度指标,以及考虑大气压降和地势海拔等因素,
13、制定出未来时间段内如何使用专业降浓剂保证土壤良好性的用剂量,这样,既达保证了经济效益又保护了土壤环境。三、问题假设1.该城区城市布局合理;2.该城区土壤重金属元素来源不受其他城区影响,数据采样都是可靠准确;3.假设每个样本点都能很好的反映该平方公里的实际情况;4.污染物的排放瞬时完成,且排放速率恒定;5.气体的传播服从扩散定律,非稳态;6.污染物的沉降速率恒定;四、符号说明iCF:元素的最高污染系数;iC:元素 i 的实测含量;iC0:元素 i 的评价标准,既背景值。P L I:污染负荷指数 ; c:土壤重金属在液相的浓度;s:土壤重金属在固相中的浓度;ijM:在区域点, i j重金属浓度;M
14、 :总重金属的浓度;tQ:土壤重某金属在 t 年后的含量;t :预测年限;n:评价元素的个数;名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页,共 28 页 - - - - - - - - - 4 z o n eP L I为区域污染负荷指数;m为评价点的个数(采样点的个数) 。五、问题一的解答5.1 空间分布对附件 1的 319个随机地点中的)海拔( mmymx),(),(做一个三维曲线拟合得该城区的城区状貌。如图1. 0123x 10400.511.52x 104-200-1
15、000100200300400图 1 该城区状貌由图 1 可以得到该城区的城区状貌, 然后利用附件 2 的 8 种金属元素浓度代替海拔()m就可以得出 8 种金属元素在该地区的空间分布图如图2 所示。00.511.522.5x 104020004000600080001000012000140001600018000-15-10-5051015202530As00.511.522.5x 104020004000600080001000012000140001600018000-1500-1000-500050010001500Cd名师资料总结 - - -精品资料欢迎下载 - - - - - -
16、 - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页,共 28 页 - - - - - - - - - 5 00.511.522.5x 104020004000600080001000012000140001600018000-1000-800-600-400-2000200400600800Cr00.511.522.5x 104020004000600080001000012000140001600018000-1000-50005001000150020002500Cu00.511.522.5x 1040200040006000800010
17、00012000140001600018000-5000050001000015000Hg00.511.522.5x 104020004000600080001000012000140001600018000-100-50050100Ni00.511.522.5x 104020004000600080001000012000140001600018000-1500-1000-5000500100015002000250030003500Zn00.511.522.5x 1040200040006000800010000120001400016000180000100200300400500Pb图
18、 2 8 种主要金属元素浓度的空间分布由图 2 的 8 种主要金属元素浓度空间分布图可对比出每种元素的浓度高涵盖点区,再结合附件 1、2 对数据进行处理,则可以得出该种元素的污染浓度范围区域,其具体情况如表 1 所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页,共 28 页 - - - - - - - - - 6 表 1 8 种元素的污染浓度范围元素区域数污染浓度范围As3 交通区 95 号;交通区 22 号和生活区 1 号地点组合;工业区29 号Cd4 交通区 95 号
19、;工业区 223号;交通区 22 号和生活区 1 号组合;工业区 6、8 号和交通区的 9 号组合Cr1 生活区 20 号和交通区 22 号组合Cu1 工业区 8 号和交通区 22 号组合Hg3 工业区 8 号和交通区 9 号组合;交通区 182 号;交通区 257 号Ni2 交通区 22 号;山区 135 号Pb2 生活区 16、20 号和公园绿地区 143 号组合;工业区 6、8 号组合Zn2 生活区 36 号;交通区 61 号5.2 污染程度污染负荷指数法模型污染负荷指数法是Tomlinson等人从事重金属污染水平的分级研究中提出来的一种评价方法。该指数又评价区域所包含的多种重金属成分共
20、同构成,并使用了求积的统计法,通过这种方法能对整个区域各个点位各种重金属进行定量评价,并对各点的污染程度进行分级,能直观反映对环境污染最严重的元素和各种元素对环境污染的贡献程度,以及重金属在时间、空间上的变化趋势,应用比较方便 1。能避免污染指数加和关系造成的对评价结果歪曲的现象,并能对任意给定的区域进行定量的判断。他也采用研究区土壤背景值为评价标准,因此它能很好的判断土壤的综合人为污染情况。首先根据某一点的测量金属含量,进行最高污染系数(CF)的计算2: iiiCCCF0(1)某一点的污染负荷指数(PLI)为:nnCFCFCFCFPLI321(2)某一区域的污染负荷指数(zonePLI)为:
21、mmz o n eP L IP L IP L IP L IP L I321(3)表 2 污染负荷指数与污染程度之间的关系PLI值121323污染等级0 1 2 3 污染程度无污染中等污染强污染极强污染名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页,共 28 页 - - - - - - - - - 7 通过Matlab程序可以计算出各个地点的最高污染系数和污染负荷指数(程序见附录1) , 并结合这些点的数据综合评价得出如表3所示的五个区域的污染负荷指数 (zonePLI) 。表
22、 3 五个区域的污染负荷指数(zonePLI)生活区工业区山区交通区公园绿地区PLI值1.8336 2.1573 1.0602 1.9209 1.5780 污染级数1 2 1 1 1 污染程度中等污染强污染中等污染中等污染中等污染由上表可知该城区内不同区域重金属的污染程度为中等污染的区域有生活区、山区、交通区、公园绿地区。强污染的区域有工业区。六、问题二的解答研究土壤重金属污染的主要原因需要调查污染土壤中各种污染物的影响因素。土壤重金属污染影响因素主要指产生重金属污染物的设备、位置、场所等污染源 3。这些污染源产生的重金属污染有的直接进入土壤,有的则是通过大气和水体再沉降进入土壤,造成城区土壤
23、的重金属污染。按重金属污染物产生的部门,土壤重金属污染源主要分为工业污染源、交通污染源、生活污染源、农药和肥料等。具体来源参见表4 中国土壤重金属的主要来源4。表 4 中国土壤重金属的主要来源重金属来源As硫酸、农药、医药、化肥、玻璃等工业废水、废气、农药Cd冶炼、电镀、燃料等工业废水、污泥和废气、肥料杂质Cr冶炼、电镀、制革、印染等工业废水和污泥Cu冶炼、铜制品生产等废水、废渣和污泥,含铜农药Hg制烧碱、 Hg 化物生产等工业废水和污泥、含Hg 农药、 Hg 蒸Ni冶炼、电镀、炼油、染料等工业废水和污泥Pb颜料、冶炼等工业废水,汽车防爆燃料尾气,农药Zn冶炼、镀锌、纺织等工业废水和污泥、废渣
24、,含锌农药同一区域可以产生许多污染物,同一种污染物的来源也不尽相同。相关分析和主成分分析等常规统计方法常用于污染物源解析中。为了更好地分析土壤重金属污染物的来源和类别,综合运用这些统计分析方法来分析5个不同区域的 8种重金属污染物的来源。最后综合 5 个不同区域的来源分析出该城区的重金属污染的主要因素。6.1 相关分析模型相关分析就是研究两个或两个以上变量之间相互关系的统计分析方法。在环境研究中常用来定性的分析研究各环境变量之间的相互关系,而变量之间的相互关系用相关系数和相关系数的显著水平来表示6。 相关分析法有多种,环境研究中最常用的是Pearson相关分析法。名师资料总结 - - -精品资
25、料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页,共 28 页 - - - - - - - - - 8 相关系数:考察两个事物(在数据里我们称之为变量)之间的相关程度。如果有两个变量: X、Y,最终计算出的相关系数的含义可以有如下理解:(1) 、当相关系数为 0 时,X和 Y两变量无关系。(2) 、当 X的值增大(减小), Y值增大(减小),两个变量为正相关,相关系数在0.00与 1.00 之间。(3) 、当 X的值增大(减小) ,Y值减小(增大),两个变量为负相关,相关系数在-1.00与 0.00 之间。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年数学建模A题国家一等奖 2022 数学 建模 国家 一等奖
限制150内