城市表层土壤重金属污染分析(修改中级版)

1、.PAGE ：.；城市表层土壤重金属污染分析修正版摘要：本文主要是经过对标题给定的种重金属对城市表层土壤污染的一些数据和信息分析建立了分散传播模型、自然沉降模型、对比模型和转换模型处理标题中的问题。针对重金属的空间分布问题，先拟出该城区的地势图，根据所给的数据绘出该地域的三维地势及采样点在其上的空间分布图。再利用matlab散乱插值法得到种重金属元素的空间分布。其次，经过单因子污染指数法和内梅罗综合指数法，得出城区内不同区域重金属的污染程度：工业区交通区生活区公园绿地域山区。即得出结论，该城区的重金属污染主要缘由在于工业释放和交通排放。在分析重金属污染传播的特征时，他们结合大气分散模型、土壤分

2、散模型、水体分散模型三个传播模型，自然沉降模型结合酸雨对重金属元素的影响，从而经过偏微分方程建立重金属元素的对流分散模型。对重金属在土壤中的浓度进展时间和空间上的拟合。当存在n个污染源时，记为第j个污染源坐标，是瞬时点源的质量，c(x,y,z,t)是t时辰x,y,z处的重金属污染物的浓度，得到与重金属污染源位置相关的偏微分方程：其中为相关系数。根据取样点重金属浓度变化确定污染源的位置，发现城市表层土壤重金属污染主要是由工业区、交通区和生活区的污染呵斥。为了对重金属浓度和客观要素优化推行，对城市地质的演化做更准确的分析和预测，需求搜集不然历史时段土壤污染程度、土壤浓度的PH值、地下水重金属的浓度

3、、植物中重金属浓度等多方面信息。经过指数法评价金属元素在土壤、PH、地下水、植物中的污染程度，采用层次分析法综合评价土壤污染程度，阐明地质环境的演化方式，再根据不同历史时间段土壤污染程度，阐明城市的地质环境变化方式。关键词：等高线图 散乱插值法 重金属分散方程 单因子污染指数法 内梅罗综合指数法目 录 TOC o - h z u HYPERLINK l _Toc 城市表层土壤重金属污染分析修正版 PAGEREF _Toc h HYPERLINK l _Toc 一、问题重述 PAGEREF _Toc h HYPERLINK l _Toc 二、问题分析 PAGEREF _Toc h HYPERLI

4、NK l _Toc . 重金属元素的空间分布及不同区域重金属污染程度 PAGEREF _Toc h HYPERLINK l _Toc .重金属污染的主要缘由 PAGEREF _Toc h HYPERLINK l _Toc .确定污染源的位置 PAGEREF _Toc h HYPERLINK l _Toc 三、模型假设 PAGEREF _Toc h HYPERLINK l _Toc 四、符号阐明 PAGEREF _Toc h HYPERLINK l _Toc 五、模型的建立和求解 PAGEREF _Toc h HYPERLINK l _Toc .重金属元素的空间分布及污染程度模型 PAGEREF

5、 _Toc h HYPERLINK l _Toc .重金属元素的空间分布 PAGEREF _Toc h HYPERLINK l _Toc .重金属的污染程度分析 PAGEREF _Toc h HYPERLINK l _Toc .、 问题二的模型建立和求解 PAGEREF _Toc h HYPERLINK l _Toc .、重金属元素相关性分析 PAGEREF _Toc h HYPERLINK l _Toc .、土壤重金属主成分分析 PAGEREF _Toc h HYPERLINK l _Toc .重金属污染物传播特性及污染源位置确定模型 PAGEREF _Toc h HYPERLINK l _

6、Toc .污染物传播特性分析 PAGEREF _Toc h HYPERLINK l _Toc .模型建立与求解 PAGEREF _Toc h HYPERLINK l _Toc . 模型的叠加及求解 PAGEREF _Toc h HYPERLINK l _Toc .城市地质环境的演化模型的建立 PAGEREF _Toc h HYPERLINK l _Toc .问题三中模型的优点 PAGEREF _Toc h HYPERLINK l _Toc .问题三中模型的缺陷 PAGEREF _Toc h HYPERLINK l _Toc .城市地质环境的演化模型分析 PAGEREF _Toc h HYPER

7、LINK l _Toc .城市地质环境的演化方式 PAGEREF _Toc h HYPERLINK l _Toc 六、模型的检验 PAGEREF _Toc h HYPERLINK l _Toc 七、参考文献 PAGEREF _Toc h HYPERLINK l _Toc 八、附录 PAGEREF _Toc h 一、问题重述随着城市经济的快速开展和城市人口的不断添加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何运用查证获得的海量数据资料开展城市环境质量评价，研讨人类活动影响下城市地质环境的演化方式，日益成为人们关注的焦点。按照功能划分，城区普通可分为生活区、工业区

8、、山区、主干道路区及公园绿地域等，分别记为类区、类区、类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进展调查。为此，将所调查的城区划分为间距公里左右的网格子区域，按照每平方公里个采样点对表层土 厘米深度进展取样、编号，并用GPS记录采样点的位置。运用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件列出了采样点的位置、海拔高度及其所属功能区等信息，附件列出了种主要重金属元素在采样点处的浓度，附件列出了种主要重金属元素的背景值。现要求他们经过数学建模来完成

9、以下义务：() 给出种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。() 经过数据分析，阐明重金属污染的主要缘由。() 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。() 分析他所建立模型的优缺陷，为更好地研讨城市地质环境的演化方式，还应搜集什么信息？有了这些信息，如何建立模型处理问题？二、问题分析. 重金属元素的空间分布及不同区域重金属污染程度标题要求统计出种重金属元素在城区的空间分布，需求思索重金属元素在X，Y，Z三维空间的浓度分布特征。根据题中所给采样点的坐标和各重金属的浓度，用散乱数据插值法获得各重金属污染物浓度的空间分布，用matlab软件绘

10、制出重金属在该区的空间分布图及种重金属的浓度等高线在三维地形图曲面的投影图。同时可以经过单因子污染指数法和内梅罗综合指数法相互验证，得到城区内不同区域重金属的污染程度。.重金属污染的主要缘由根据附件给出的数据，他们分区域思索，确定每个区域里种重金属元素相互之间的相关性，相关性较显著的重金属元素可以以为于同一个污染源，由此可以大致确定每个区域污染源的种类。 由于每个区域中重金属的种类较多，相互之间又具有比较明显的相关性，所以可以在采用KMO法和Bartlett法对原始数据进展主成分分析适宜性检验的根底上，进一步思索运用主成分分析法，由此分析出每个区域影响比较显著的假设干种重金属，从而确定出重金属

11、污染的主要缘由。.确定污染源的位置由于重金属污染物主要是经过大气、土壤和水体进展传播，同时能够遭到酸雨淋融和冲刷三维影响，然后自然沉降吸附在土壤中。故遭到五方面要素影响，且影响要素相互关联，相互影响。经过建立对流方程建立模型，从中表达重金属传播的主要特征。在随时间的传播过程中，污染源所在地的浓度随时间的变化率较大。对重金属在土壤中的浓度进展时间和空间上的拟合，由此可以同对模型的反演确定污染源的大约位置。.重金属浓度和客观要素优化模型在研讨污染物传播的过程中，此题只给出了一个时间点上的数据，信息量明显缺乏，应搜集一段时间内的数据。由此建立模型可推出不同时间不同地点污染物的浓度。为确保模型的准确性

12、，还要思索到土壤浓度的PH值、地下水重金属的浓度、植物中重金属浓度等环境信息，以建立更真实的模型，对城市地质环境的演化做更准确的分析和预测。三、模型假设.固定时间内污染元素浓度坚持稳定。.取样点的浓度较好地反映了该地域污染物的浓度。.重金属在无穷空间分散，不计风力影响。.重金属的分散符合分散定律，即单位时间经过单法相面积的流量与它的浓度梯度成正比。.不思索不同区域内大气、土壤、水体等自然环境差别。.元素污染浓度在程度面与海拔高度的分布成独立相关性。.污染元素的分散只受程度面的分散影响。.污染元素在程度面得分散在各方向的分散速度一样。.污染源的位置大致固定，不随时间变化。四、符号阐明符号变量符号

13、阐明第个取样点污染指数第个取样点污染浓度丈量值第种重金属元素评价规范内梅罗污染综合指数空间任一点 在时辰的污染物浓度污染物分散系数悬浮泥沙对重金属的吸附速率重金属在悬浮泥沙中的解析速率重金属离子在土壤中固态和液态的分配系数污染物的总污染量分散系数重金属离子在土壤中迁移的阻滞因子混合沙的吸附量悬移质平均吸附量土壤介质的密度泥沙干容重断面平均水深五、模型的建立和求解.重金属元素的空间分布及污染程度模型.重金属元素的空间分布根据采样点各重金属元素的含量分布，运用matlab软件将该地域的三维地势及采样点在其上的综合空间分布图绘制如下 重金属As元素的空间分布图及海拔分布图 -由图-分析得，AS元素浓

14、度浓度较高部分在左图三维图中明显突出的三个柱体中，阐明这三个区域AS元素污染较为严重。同时思索浓度与海拔的关系，右图呈现各浓度主要分布在海拔高度为-m之间的范围内，且与海拔无特别关系。 重金属Cd元素的分布空间及海拔分布图-由图-分析得，Cd元素浓度较高部分在左图三维图中相对突出的几个柱体只可以，且该元素在城区内的污染较为严重。而由右图得在海拔高度为-m的范围内均有浓度分布，且浓度较高部分分布在海拔较低区域内。 重金属Cr元素的空间分布及海拔分布图-由图-分析得，Cr元素在城区总体分布较少，但有一个很明显的突出柱体区域。阐明污染主要集中在这一块。同时海拔图显示污染主要分布于-m海拔之间，且浓度

15、量大致在-g/g之间。 重金属元素Cu元素的空间分布及海拔分布图-由图-分析得，Cu元素在城区内的总体分布较少，高浓度主要集中于两个明显突出的柱体的区域，外表这些区域污染较为严重。右图得海拔图显示污染主要分布于-m海拔之间，且与海拔无特别关系。 重金属元素Hg元素的空间分布及海拔分布图-由图-分析得，Hg元素在城区内的总体分布较少，高浓度主要集中于三个明显突出的柱体的区域，阐明这些区域污染较严重。右图得海拔图显示污染主要分布于-m海拔之间，且与海拔无特别关系。 重金属元素Ni元素的空间分布及海拔分布图-由图-分析得，Ni元素在城区内的总体分布较少，高浓度主要集中于两个明显突出柱体的区域，外表这

16、些区域污染较为严重。右图得海拔图显示污染主要分布于-m海拔之间，且与海拔无特殊关系。 重金属Pb元素污染的空间分布及海拔分布图-由图-分析得，Pb元素浓度较高部分分布在左图中两个明显突出的柱体区域，阐明这两个区域Pb元素污染较严重。右图得在海拔高度为-m的范围内均有浓度分布，且浓度较高的部分大致分布在海拔较低区域内。 重金属Zn元素污染的空间分布及海拔分布图-由图-分析得，Zn元素浓度较高的部分在左图中几个突出柱体区域，阐明这几个区域Zn元素污染较为严重。由右图得在海拔高度为-m的范围内均有浓度分布，且与海拔无特殊关系。综上所述，由于各种重金属有不同的来源，比如交通污染，工业污染，生活渣滓污染

17、等，所以呈现出不同的区域分布。但从海拔方面来分析，其主要分布在-m的低海拔区域。.重金属的污染程度分析a.模型建立与求解他们运用地积累指数法建立模型。Muller曾提出地积累指数定量评价堆积物中的重金属污染程度，并规定了相应的污染程度级别划分规范表-，这种评价方法也可以用来评价土壤中的重金属的污染程度及其分级情况。其计算公式为： 式中Ci为重金属元素i在样本中的浓度，Bi为重金属元素i的背景值，K是为思索各地岩石差别能够会引起背景值的变动而取的系数普通取值K=.。表- Muller 地积累指数分级表根据附件表格给出的各重金属元素在多个采样点的浓度，同时他们也知道各个采样点对应的区域，所以他们可

18、以用Matlab计算出每种重金属在每个区域的平均浓度Ci。Bi取各元素背景值的平均值代入计算。那么可以算出各元素在各个区域的地积累指数Igeo。表-生活区重金属污染程度表-工业区重金属污染程度表-山区重金属污染程度表-主干道路区重金属污染程度表-公园绿地域重金属污染程度为了全面反映各污染物对土壤的不同作用，突出高浓度污染物对环境质量的影响，采用内梅罗综合污染指数法再次建立一个模型。内梅罗指数不仅思索到各种污染物的平均污染情况，并且强调最严重的污染物对污染指数的影响。因此，可采用单因子污染指数( PI) 和内梅罗综合污染指数( IPI) 对城区土壤重金属的污染程度进展综合评价。单因子污染指数法：

19、式中，PI 为重金属污染物i的污染指数;Ci为重金属污染物i的实测浓度，Si为重金属污染物i的国家规范值。内梅罗综合污染指数法:式中，Pmax为一个区域采样点中重金属单项污染指数最大值；Pave为一个区域一切重金属单项污染指数的平均值。IPI为一个区域土壤中重金属综合污染指数。评判规范如下：IPI，表示无污染； IPI，表示轻度污染； ，表示严重污染。由此可以进一步计算出内梅罗综合污染指数( IPI)。他们把城区分为个区域分别计算。计算过程中Ci取一个区域中重金属i的采样值的平均浓度，Si取为重金属i的国家规范值。经过Matlab编程计算出的结果如下表-所示。表-各区域综合污染程度.、 问题二

20、的模型建立和求解.、重金属元素相关性分析城市土壤重金属来源于成土母质和人类活动，同一来源的重金属之间存在着相关性，根据相关性可以判别土壤重金属污染来源能否一样。假设重金属之间存在显著的正相关，那么其来源能够一样，否那么来源能够不止一个。运用Matlab编程，可以计算出各个区域内种重金属相互之间的相关性，并对其相关性做出分析。表-生活区重金属元素相关性在生活区采集的土壤样品中，As和Ni，Cd和Pb呈现比较明显的相关性，相关系数分别是.和.，阐明As和Ni能够具有一样的污染源，Cd和Pb能够具有一样的污染源。但其他元素之间的相关性都不明显，所以生活区的重金属来源相对复杂，很难于某种污染源进展确认

21、。表-工业区重金属元素相关性在工业区采集的土壤样品中，Cd、Cr、Cu、Hg、Pb、Zn两两之间都有比较明显的相关性，阐明这六种重金属元素能够具有一样的污染源，As和Ni具有明显的相关性，能够具有一样的污染源。表-山区重金属元素相关性在山区采集的土壤样品中，Cd、Pb、Zn两两之间呈现明显的相关性，所以此三者能够有一样的污染源。As和Cu、Cu和Hg之间有明显的相关性Cr和Ni、Ni和Zn也有比较明显的相关性。表-交通区重金属元素相关性在主干道路区采集的土壤样品中，Cd、Cr、Cu、Ni、Pb、Zn两两之间具有明显的相关性，所以能够有一样的污染源。As、Hg与其本身外七种重金属的相关性都不明显

22、，两者与其他重金属的来源能够不同。表-公园区重金属元素相关性在公园绿地域采集的土壤样品中，As、Cr、Ni两两之间的相关性比较明显，能够有一样的污染源。Cu、Pb、Zn两两之间的相关性比较明显，能够有一样的污染源。.、土壤重金属主成分分析为了进一步明确该城区不同功能区土壤中重金属的污染来源, 利用Matlab软件对上述不同功能区土壤重金属的项目的进展主成分分析。在利用主成分分析进展源解析前,采用KMO法和Bartlett法对原始数据集进展主成分分析适宜性检验。经计算, 其KMO值为.大于. 的最小值要求，Bartlett检验统计的显著程度小于. ,因此可以为相关系数矩阵与单位阵有显著差别,以上

23、检验结果均阐明原始数据集适宜进展主成分分析。数据处置及分析结果见下表。表-生活区主成分分析 表-工业区主成分分析生活区的种重金属的全部信息可由个主成分表示。即对前个主成分进展分析曾经可以反映全部数据的大部分信息。结合表中数据可以看出:三个主成分中第一主成分的奉献率为.%，特点表现为因子变量在元素Pb、Cd上有较高的载荷。研讨结果阐明汽车尾气颗粒物中铅含量较高,在怠速形状下向大气的排放量分别为 mg/min。又知Pb和Cd在该区相关性明显，所以大致可以为交通污染是生活区的一个主要污染缘由。第二主成分的奉献率为.%，特点表现为因子变量在元素As、Zn上有较高的载荷，由于生活渣滓中会富含比较多的As

24、和Zn，所以生活区生活渣滓的堆积能够是重金属污染的主要缘由。第三主成分的奉献率为.%，特点表现为因子变量在元素Ni上有较高的载荷，普遍以为，Ni和Zn主要来源于轮胎的磨损，所以也阐明交通污染能够是生活区的污染来源之一。工业区的种重金属的全部信息可由个主成分表示。第一主成分的奉献率为.%，特点表现为因子变量在元素Cr、Cu、Hg、Pb、Zn上有较高的载荷,且这些重金属在该区相关性显著，可以以为来源一样，工业污染排放是工业区重金属污染的主要缘由。第二主成分的奉献率为.%，特点表现为因子变量在元素As、Hg上有较高的载荷, 这两种元素通常被视为燃煤和渣滓熄灭的标志元素，所以燃煤污染也是工业区重金属污

25、染的一个主要来源。表-山区主成分分析 表-主干道路区主成分分析山区的种重金属的全部信息可以由个主成分表示。第一主成分的奉献率为.%，特点表现为因子变量在元素Zn上有较高的载荷。第二主成分的奉献率为.%，特点表现为因子变量在元素Cd、Pb上有较高的载荷，第三主成分的奉献率为.%，特点表现为因子变量在元素Hg上有较高的载荷。综合思索，由于山区的污染程度较轻，所以重金属污染的主要缘由能够是山区富含较多的矿石。主干道路区的种重金属的全部信息可以由个主成分表示。第一主成分的奉献率为.%，特点表现为因子变量在元素Pb、Ni、Cu、Cd上有较高的载荷，所以交通污染中的汽车尾气排放和轮胎磨损能够是重金属污染的

26、主要缘由。第二主成分中Hg具有较高的载荷，所以也能够遭到工业污染如仪表、电气工业等的影响。第三主成分中As具有较高的载荷，所以能够遭到工业上排放As的工厂如化工、冶金业等的影响。表-公园绿地域主成分公园绿地域的种重金属的全部信息可以由个主成分表示。第一主成分的奉献率为.%，特点表现为因子变量在元素Cd、Cr、Pb、Zn上有较高的载荷，这些重金属元素主要来源可以以为是交通污染。 第三主成分中Hg具有较高的载荷，有研讨发现由于大量的荧光灯、汞灯等气体放电灯的广泛运用,使得绿地表层土壤中汞的含量远远超越背景值。所以荧光灯等的大量运用能够是Hg污染的主要缘由。.重金属污染物传播特性及污染源位置确定模型

27、.污染物传播特性分析重金属元素传播途径主要有水流传播、土壤浸透、大气分散等。在城区他们主要思索以大气分散方式传播。由此，他们经过分散方程建立模型，从中表达重金属传播的主要特性。.模型建立与求解A大气分散模型大气中重金属主要来源于工业消费、汽车尾气排放、汽车轮胎磨损产生的大量含有重金属的有害气体和粉尘。它们主要分布在工矿的周围及公路、铁路两侧。思索这些有害气体和粉尘在大气中分散的情况。将污染物开场释放的时辰记作，污染源选为坐标原点。时辰无穷空间中恣意一点的污染物浓度记为。根据假设，单位时间经过单位法向面积的流量 .是分散系数，表示梯度，负号表示由浓度高向浓度低的地方分散。调查空间域，的体积为，包

28、围的曲面为，的外法线向量为，那么在内经过的流量为 .而内污染物的增量为 .由质量守恒定律 .根据曲面积分的奥氏公式 .其中是散度记号。由.式再利用积分中值定理不难得到 .这是无界区域的抛物型偏微分方程。根据假设，初始条件为作用在坐标原点的点源函数，可记作 .表示污染物的总污染量，是单位强度的点源函数。方程.满足条件.的解为 .这个结果阐明，对于恣意时辰重金属浓度的等值面是球面，并且随着球面半径的添加的值是延续减少的；当或时B.土壤分散模型金属离子在多孔介质中的迁移取决于在水相中的分子分散、固相上的吸附、固相的紧缩比以及固相中的弯曲度，采用一维迁移方程计算金属元素在土壤中的迁移和浓度分布 ： .

29、其中为金属离子在溶液中的浓度，为分散系数，为水的通量密度，为吸附在土壤上金属离子的浓度，为时间，为被调查点到金属离子走漏源的间隔 ，为土壤介质的密度，为空隙率。金属离子在土壤介质中的吸附公式为 .其中是金属离子在土壤中固相和液相的分配系数。将式.代入式.： .其中是金属离子在土壤中迁移的阻滞因子： .金属离子在土壤中迁移的初始条件和边境条件为 .从式.可以得到以下分析解 ： .是误差补余函数， .在此不思索土壤密度受外部压力、地理位置的影响。假设污染源处重金属离子的相对浓度为，取样点重金属离子浓度均为与污染源重金属离子的相对浓度。在较小时，重金属离子在土壤中的相对浓度随间隔 的添加而迅速减小。

30、C.水体分散模型在此模型中，他们将问题理想化，以为水体传播仅仅是靠河流传播的。假设河流没有多大的改动。泥沙泥沙颗粒在水溶液中对重金属污染物的吸附遵照兰格茂尔吸附运动学方程 .在天然水体中，悬浮泥沙对重金属的吸附，和静态实验中的吸附本质并无不同，对部分河段可以为水环境化学条件根本一样。只是水流泥沙要素及水相浓度和泥沙吸附量与空间位置有关。对微小水体，可以为与静态实验一样。仍可用兰格茂尔吸附运动学方程，只需将对时间的微分改为偏微分，即对均匀沙有： .水流泥沙运动数学模拟研讨已有较长的历史，在水流泥沙运动确定后, 根据建立的重金属迁移转化方程, 吸附动力学方程及边境、初始条件就可确定河流重金属污染物

31、迁移转化过程, 构成河流一维重金属污染物迁移转化数学模型的根本方程如下：水流延续方程 .水流运动方程 .泥沙延续方程 .河床变形方程 .有，那么.假定悬浮泥沙的吸附量，对.式全断面积分得： .处理重金属迁移转换模型可将.至.式简化成： .根据投放泥沙静态吸附实验的结果，取b=.g/kg，k=. /(ppms)，k=. /s。重金属污染物的初始条件为： .；边境条件为： .在一样条件下，污染物的运动遵照以下方程： .用向后显式差分求解，结果如图-:图- 泥沙吸附对水相浓度的影响由上图可知：在不同时辰, 水相重金属浓度由于泥沙吸附, 峰值明显减少, 而且由于泥沙吸附需求时间, 因此河段上泥沙吸附量

32、的峰值总是滞后于水相浓度的峰值。污染物不受泥沙运动影响, 在输移过程中, 由于纵向综合紊动分散作用, 峰型变化一样条件下重金属污染物输移, 虽然图形类似, 但重金属污染物的峰值因泥沙吸附而下降很多, 特别是后半部减小更为明显, 且随时间添加明显加大。所以，重金属在水体中的传播特征是随着离污染源的间隔 添加而减少,速度也随之减小。D.自然沉降模型当固体处于流体中时，只需两者的密度有差别，那么在重力场中颗粒将在重力方向与流体做相对运动。直径为的球形颗粒遭到的重力为： .其中为颗粒密度。颗粒处于流体中，无论运动与否，都会遭到浮力。当流体处于重力场中，颗粒遭到的浮力等于： .分析颗粒沉降运动必需思索流

33、体对颗粒运动的阻力。两种阻力：包括表皮阻力和沉降阻力。当颗粒速度很小时，流体对球的运动阻力主要是粘性摩擦或表皮阻力。阻力大小的计算仿看管路阻力的计算，即以为阻力与相对运动速度的平方成正比，对于直径为球形颗粒： . 在重力场中，颗粒在流体中遭到重力、浮力和阻力。这些里会使颗粒产生一个加速度。根据牛顿第二定律：。当颗粒在流体中做匀速运动时， .现实上，颗粒从静止开场作沉降运动时，分为加速和匀速两个阶段。速度越大阻力越大，加速度越小；加速度为零时颗粒便做匀速运动，其速度称为沉降速度。普通而言，对小颗粒加速阶段时间很短，通常忽略，可以以为沉降过程是匀速的。令颗粒所受合力为零，便可解出沉降速度，解得：

34、.计算沉降速度首先要知道阻力系数，经过因次分析法可知它是颗粒与流体相对运动雷诺数的函数：而。大气层流区沉降所涉及的颗粒普通都很小，；颗粒在较安静的水体中沉降，即处于过渡区Allen区时，；颗粒处于湍流区牛顿区时，。E.酸雨淋溶对重金属元素的影响根据相关资料，在不同pH值的酸雨溶液中，酸性越强，重金属溶出率越高。元素Cu，Ni和Cr在较小粒径下小颗粒较容易溶出，这是由于颗粒物粒径较小，其比外表积越大，重金属可与酸液充分作用而较容易溶出；元素Zn在两粒径下溶出率相差较小。两粒径下的颗粒物中，元素Pb和Cr溶出率均随pH值降低而增高比较显著，重金属Cu和Zn在所测pH值范围内溶出率变化均较小。不同硫

35、氮比的酸雨溶液对重金属元素的溶出率也有一定的差别，整体上看，重金属元素的溶出会随着酸雨中含氮比例的添加而有所增大，阐明随着含氮比例的升高，酸雨的酸性会有所加强，对重金属的溶出才干也得以提高。由于并未给出该城区酸雨的相关信息，可以将酸雨淋溶作用作为重金属元素分散的参考要素，特别是淋溶作用明显的Cu，Ni，Cr等元素，很有能够在空间分布上遭到酸雨影响。. 模型的叠加及求解根据表层土壤重金属元素在采样点的浓度观测值，结合表层土壤重金属的传播特征，建立了关于偏微分方程的数学模型，求解重金属元素污染源的详细位置。记是时辰处的重金属污染物的浓度，建立一阶衰减的非稳定态的质量分散运输方程，即： .其中： 是

36、方向相应的分散系数；是衰减系数；为瞬时点源的质量；为污染源坐标。利用MATLAB软件编程，可知上述方程的解析解为： .针对每种重金属元素，参考元素浓度空间分布图，得到能够的污染源的个数，那么空间任一点在时辰的元素浓度： .根据上述重金属的传播特性模型，并分析附表 和附表 的数据值，结合第一问得等高线的分析方法，他们可以综合分析出污染源的位置。在处置数据时，明显看出有一些点的金属浓度与其他的浓度值不同，特别是 Hg 的最高浓度分布的那四个点，与其他位置的 Hg 浓度存在数量级的差别，运用概率的数理统计分析和整理方法，可知其为不正常点。剔除这几组不正常值，再用类似的方法分析其他的污染物，从而可以得

37、四处置了的数据，在上述传播特征模型的根底上，分析出污染源的位置如表-所示。表- 种重金属元素主要污染源位置重金属元素主要污染源位置As(,.),(.,.,.)Cd(,.)Cr(.,.,.)Cu(.,.,.)Hg(.,.,.),(,)Ni(.,.,.)Pb(.,.,.),(,)Zn(,.,.)表中数听阐明，本城区的污染主要是由于工业区、交通区和生活区的污染呵斥。污染源所在点即为该重金属浓度较高的点，普通为最高点和次高点为污染源。由于污染源处在不同的位置，而其受重力等传播特征的影响，会在传播中发生分散、降解等物理化学反响，其浓度在污染源的周围呈指数函数方式快速降低，构成较强的梯度差。传播的存在使城

38、区的重金属污染物在其他非污染源点也出现了较高的浓度，从而出现部分的轻度污染。.城市地质环境的演化模型的建立.问题三中模型的优点a. 模型分析了重金属污染在时间和空间上的传播特性，使其不局限于数据所给出的的单一时间点，更具推行性。.b. 重金属污染受多重要素影响，模型将大气分散模型、土壤分散模型、水体分散模型，以及自然沉降模型进展时间和空间上的拟合，结合酸雨淋溶对重金属元素的影响，使其更贴合实践，更具压服力。.问题三中模型的缺陷a. 仅分析一个时间点上的数据，且不思索该城区各取样点的土质，河流，降水量等环境要素的影响，模型简单化。b. 在模型中出现多种假设，将其理想化，能够与实践有差别。.城市地

39、质环境的演化模型分析在此题中，给出了采样点的位置、海拔高度及所属功能区、各个重金属元素在采样点处的浓度及相应的背景值。但分散传播模型所确定的偏微分方程不仅与空间位置相关，也与传播时间亲密相关，同时遭到外界环境的多重影响。因此为了更好地研讨地质环境的演化方式，他们应该搜集以下信息：a.土壤浓度的PH值土壤的PH值比土壤矿物学更具重要位置。由于土壤浓度的pH值不仅决议了各种土壤矿物的溶解度，而且影响着土壤溶液中各种离子在固相上的吸附程度。首先，随土壤溶液pH升高，各种重金属元素在土壤固相上的吸附量和吸附才干加强。其次，土壤溶液的pH影响土壤溶液中重金属元素离子活度。最后，土壤溶液pH对土壤中重金属

40、元素的生物有效性影响能够不是单一的递增关系。b.地下水重金属的浓度地下水重金属污染是污染物在地下水系统中进展各种物理、化学和生物过程的结果。在一定条件下,土壤中的重金属可随水的挪动而发生迁移,从而要挟人类的安康与其它动物的繁衍生息，因此积累在土壤中重金属的挪动性遭到人们广泛关注。普通以为,新土层对重金属有一定的吸附才干,重金属在土壤中的挪动性很有限,经过田间实地采集剖面土样的分析结果也阐明, 重金属主要积累在表层。土壤重金属的生物毒性不仅与其总量有关,更大程度上由其形状分布所决议,不同的形状产生不同的环境效应,直接影响到重金属的毒性、迁移及在自然界的循环。c.植物中重金属浓度植物对金属的吸收主要取决于植物从土壤中汲取金属以及向地上部运输金属的才干。植物对金属的吸收主要取决于自在态离子活度。所以重金属浓度与植物体内中重金属浓度有正相关性，经过植物中重金属浓度来评价污染程度。.城市地质环境的演化方式根据以上数据，先用指数法分别评价金属元素在土壤、PH、地下水、植物中的污染程度，然后采用层次分析法综合评价土壤污染程度，阐明地质环境的演化方式。再不同历