城区表层土壤重金属污染分析模型南开大学数学建模.docx

《城区表层土壤重金属污染分析模型南开大学数学建模.docx》由会员分享，可在线阅读，更多相关《城区表层土壤重金属污染分析模型南开大学数学建模.docx（30页珍藏版）》请在三一办公上搜索。

1、城区表层土壤重金属污染分析模型摘要 本文针对城市表层土壤重金属污染的分析与评价，在综合考虑了各种不定因素的影响条件下，综合运用模糊数学对各种重金属的污染程度给出判定。利用基础几何与回归分析结合的方法建立了确定污染源的数学模型。并用数学软件对数据进行很好的处理，得到较为准确的数字信息，解决了所需解决的问题。模型一中，对于重金属元素的空间分布及该城区内不同区重金属的污染程度的问题，文中根据单个元素在不同地域的污染程度，准确画出各种重金属元素的空间分布图。依据土壤重金属污染程度模糊性，采用模糊隶属函数模型推导得出数据，辅以加权平均模型，得出各个区的污染级别。生活区、工业区、山区、主干道路区及公园绿地

2、区的污染级别分别时三级、三级、二级、三级、三级，从而断言其污染程度。并通过相关性系数分析法来判断元素之间的相关关系，并结合重金属元素的空间分布图说明重金属污染的主要原因。模型二中，对于重金属元素的传播特征的理解是它的分布特点和动态传播规律。元素分布的特点是通过模型一中的重金属元素的空间分布结合它们之间的相关性来描述的。元素的动态变化趋势函数是通过基于基本几何和区域散点线性关系得出的，即：通过取As重金属污染浓度空间分布图中的一块图形进行定性的分析得出潜意识矢量方程。通过计算得到污染源的浓度为。再通过将计算得出的浓度带入动态变化趋势函数得到坐标。这个点的坐标就是污染源的坐标。模型三，模型三基于物

3、力、物理化学、流体力学等基础学科的知识建立了非保守性污染物质在土壤中迁移转化一维控制方程。通过公式的推倒，得到方程： 得到的结论是:污染物浓度在时空上变化规律的近似解析解在形式上与正态概率分布函数相同的。 关键字 模糊数学 相关性系数分析法 隶属函数 一维弥漫方程 渗透方程1问题重述11 问题的提出城市规划原理中，对于城区功能区的划分一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（010 厘米深度

4、）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。通过数据分析，说明重金属污染的主要原因。(2) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置，并分析所建立模型的优缺点。(3)为了更好地研究城市地质环境的演变模式，还应收集什么信息？以及有了这些信息，如何建立模型解决问题？12 问题的分析针对问题一 不同的功能区对土壤的影响大小是

5、不同的。要评定8种重金属在该城区的空间分布和对该城区不同区域的污染程度，仅用数据描述过于抽象。因此采用matlab 成图技术直接观察这8种主要重金属元素在该城区的空间分布，这样比仅用数据描述更为形象。在评定重金属元素污染程度的时候，由于考虑到元素毒性问题，不能仅凭元素所占含量来划分它们的污染程度，因此必须对元素进行权重的衡量。要说明重金属元素污染的主要原因，需要结合它们的空间分布图来说明。不同元素的污染源不同，但是它们之间存在着一定的相关关系，从这种关系里找出它们污染的主要原因。 针对问题二 要得到传播方程，必须从污染物浓度随传播路径的变化规律中找到一般的方程。它们的关系应类似于的关系。得到传

6、播方程后，对给定的区域取一些离散的点，并做出回归方程。由方程进一步推倒出污染源的位置。针对问题三 由于问题二中所建立的方程过于简陋，在现实生活中运用不大，只能作为一个理想的模型。问题三所要解决的就是，实现方程与实际生活关系。能运用方程解决现与实相关的问题。2模型的假设与符号的约定21 模型的假设1 每个监测点一种重金属污染浓度平均值就是该区污染的评价值2 背景值作为该区的污染评价标准3 元素与元素之间不发生化学反应，它们之间的性质在短时期内不会改变22 符号的约定 第种污染物的检测浓度 第种污染物评价标准浓度 该重金属含量的实测值 第个因子的权重 某种重金属元素对一级重金属环境质量的隶属度函数

7、 某种重金属元素对二级重金属环境质量的隶属度函数 某种重金属元素对三级重金属环境质量的隶属度函数 该重金属对应于一级土壤重金属环境质量状况的标准值 该重金属对应于二级土壤重金属环境质量状况的标准值 该重金属对应于三级土壤重金属环境质量状况的标准值U 污染物评价因素集合V 评价等级集合组成一个模糊矩阵,即关系模糊矩阵该指标对应的各土壤重金属环境质量级别的标准值3 模型的建立与求解31 建模过程模型1.1对于重金属元素在该城区的空间分布，本文采用matlab 绘制出实验的319个样点的GPS坐标，生成相应的散点图，录入各点重金属浓度，并以此浓度作图生成等值线图，最终得到重金属空间分布格局图。具体见

8、附录。 8种主要重金属元素在该城区的空间分布图图一 土壤重金属As含量（g/g）空间分布图图二 土壤重金属Cd含量（ng/g）空间分布图图三 土壤重金属Cr含量（g/g）空间分布图图四 土壤重金属Cu含量（g/g）空间分布图图五 土壤重金属Hg含量（ng/g）空间分布图图六 土壤重金属Ni含量（g/g）空间分布图图七 土壤重金属Pb含量（g/g）空间分布图图八 土壤重金属Zn含量（g/g）空间分布图注：点（.）、 小圆（o） 、叉子符（x）、 加号（+）、星号（*），分别为生活区、工业区、山区、交通区、公园绿地区区的标志。模型1-2 对题目所给的附录2表格中，筛选出各金属元素的区域分布情况，并

9、把单个区域所存在的元素浓度取平均值得到各检测点浓度表。同时，把附录3中的背景值浓度做5个划分，分别为无污染、尚污染、轻污染、中污染和重污染。表一 城区各功能区监测点评价浓度表As (g/g)Cd (ug/g)Cr (g/g)Cu (g/g)Hg (ug/g)Ni (g/g)Pb (g/g)Zn (g/g)一区6.27050.2899669.01849.4030.09304118.34269.106237.01二区7.25140.3931153.409127.540.6423619.81293.041277.93三区4.04410.1523238.9617.3170.04095615.45436

10、.55673.294四区5.7080.3600158.05462.2150.4468217.61763.534242.85五区6.26370.2805443.63630.1920.1149915.2960.709154.24评价等级一级（背景值）二级三级污染程度无污染尚污染轻污染中污染重污染As (g/g)1.82.73.64.55.4Cd (ug/g)00701013016019Cr (g/g)1322314049Cu (g/g)69.613.216.820.4Hg (ug/g)00190027003500430051Ni (g/g)4.78.512.316.119.9Pb (g/g)192

11、5313743Zn (g/g)4155698397表二 该地区土地污染等级评价标准由于考虑到元素的毒性，必须对各个元素进行权重的确定，确定权重的公式这里采用反映土壤各种重金属元素相对含量大小的加权法 该计算权重的方法在一定程度上可以反映污染超标的轻重对因子权重的影响 该方法计算权重的一般公式为： 并用matlab编写权重矩阵。表三: 由matlab计算得到元素毒性权重表As (g/g)Cd (ug/g)Cr (g/g)Cu (g/g)Hg (ug/g)Ni (g/g)Pb (g/g)Zn (g/g)0.088170.112910.11270.189450.134560.075480.11284

12、0.173880.046390.069640.039680.222550.422720.03710.069120.092770.117850.122920.131850.137630.122770.131810.123720.111440.051630.090180.060980.153480.415730.046640.066740.114610.106640.132270.086270.140190.201370.076190.120030.13701为了进行模糊运算需要确定隶属度函数并以隶属度来描述土壤污染状况的模糊界线使用表2的值作为各个隶属度函数的拐点,用分段的直线函数模拟,可以得到

13、不同PH 值下各种重金属污染状况的隶属度函数. 一般来说 ,各种土壤重金属环境质量状况的隶属度函数可以用下面三个分段函数表示某种重金属元素对一级重金属环境质量的隶属度函数：某种重金属元素对二级重金属环境质量的隶属度函数：某种重金属元素对三级重金属环境质量的隶属度函数：加权平均模型式中为最终评价结果对应于第j个等级的隶属度，为对应的权重，为模糊关系矩阵中的对应元素，为参评因子个数，为所划分的等级数，类似的计算得到评价向量由于该模型计箅结果已经自动归一化，所以以该集合中最大值所对应的级别作为最终评价结果3。用计算模糊矩阵：对于矩阵，有5个矩阵，分别有，模糊矩阵复合运算模型的评价：评价向量的计算：计

14、算各采样点的关系模糊矩阵和权重模糊矩阵，分别用单因素决定模型和加权平均模型进行评价。 对评价结果进行归一化处理，得到各采样点的评价向量评价分值的计算：根据土壤重金属环境质量各级别的控制意义，取土壤重金属环境质量一级水平的分值为100,二级水平的分值为80三级水平的分值为60,按照公式(4) ,取k 值为2,计算在两种评价模型下各采样点的评价分值。从评价分值可以直观的看出各采样点土壤重金属环境质量的优劣，可比性大大增强。表四 各区污染程度分析表评价模型评价结果一区二区三区四区五区加权平均模型一级隶属00000二级隶属0.01540.000430360.606230.0140110.071929三

15、级隶属0.98450.999570.386880.985990.92807评价等级三级三级二级三级三级污染程度严重污染严重污染中度污染重污染重污染评价分值60.3160.0171.7160.2861.44321模型解释依据表四的结果，按照最大隶属的原则，可以评价出各个区的土壤重金属环境质量级别。但需要注意的是，当某个区具有两个最大隶属度时，需要考虑次大值更靠近哪个最大值，以该最大值对应的级别为最终评价结果。模型优缺点：加权评价模型主要适用于各个参评因素超标情况接近，即不存在单因素否决，且评价出发点为希望评价体现不同参评因子对土壤重金属环境质量的综合影响的情况。而单因素决定模型适用于个别参评重金

16、属元素超标过大，严重影响土壤环境质量，评价出发点为希望评价体现单因素否决的情况。根据实际情况，采用加权评价模型能够更准确的评价环境污染程度，但是对于山区的评价有待商榷。因为环境没有过多的人类活动的时刻影响，所以外来的重金属沉降积累很容易占较大的权重。由此，评价结果缺乏可信度。33重金属污染的主要原因的分析使用中的相关性系数分析法，对8种重金属元素，按功能区分类，进行两两之间的相关性指数分析2。 表五 土壤中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 10.255*0.189*0.160*0.0640.317*0.290*0.247*Cd 10.352*0.397*0

17、.265*0.329*0.660*0.431*Cr 10.532*0.1030.716*0.383*0.424*Cu 10.417*0.495*0.520*0.387*Hg 10.1030.298*0.196*Ni 10.307*0.436*Pb 10.494*Zn 1*. 在 .01 水平（双侧）上显著相关。观察发现，7种元素两两之间的相关性十分明显。除了之间的相关性为显著相关外，其余的重金属之间都达到了极显著相关关系。将重金属元素两两为一组，按照相关系数R值的由高到低来排列，得出前7位的顺序：Ni-CrPb-CdCu-CrPb-CuZn-Pb，它们之间的相关系数分别为：0.716，0.66

18、0，0.520，0.494，0.436。其余重金属间的相关系数都小于04。由此表明：Ni-Cr Pb-Cd Cu-Cr Pb-Cu Zn-Pb。这7组重金属之间的相关性实验中，对土壤重金属分布起到主导作用。相关性分析可以帮助揭示重金属元素的污染来源，按土壤所在的不同功能区分类，将该城区中的319个样点的重金属数据进行两两之间的相关性指数分析。表六 生活区中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 10.381*0.2380.531*0.2930.605*0.450*-0.017Cd 10.349*0.499*0.397*0.2830.802*0.346*Cr 10

19、.376*0.1500.527*0.416*0.412*Cu 10.1980.434*0.502*0.238Hg 10.2110.340*0.242Ni 10.300*0.334*Pb 10.328*Zn 1*. 在 0.05 水平（双侧）上显著相关。*. 在 .01 水平（双侧）上显著相关。观测发现，生活区土壤金属元素中与As含量与Ni、Pb含量，Cd含量与Cu、Pb含量，Cr含量与Ni、Pb、Zn含量，Cu含量与Ni、Pb含量, 达到了O01水平的极显著正相关。相关系数分别是0.605，0.531,0.450，0.499，0.802，0.527，0.416，0.412，0.434，0.50

20、2。As含量与Cd含量，Cd含量与Cr、Zn含量，Cr含量与含量含量，Hg含量与Pb含量, Ni含量与Pb、Zn含量, Pb含量与Zn含量, 达到了O.05水平的显著正相关。由于Cd和Pb主要来自交通污染源，而Cr、Zn来自自非交通污染源。这说明了在生活区，以工业污染源为主的重金属元素之间有着明显的交互作用，交通污染源为主的重金属元素之间也有着明显的交互作用，但这两种污染源之间的交互作用并不明显。表七 工业区中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 10.3290.380*0.1530.1810.690*0.395*0.518*Cd 10.541*0.566*0

21、.533*0.489*0.829*0.754*Cr 10.920*0.902*0.698*0.675*0.695*Cu 1.983*0.503*0.670*0.622*Hg 10.479*0.612*0.590*Ni 10.578*0.634*Pb 10.739*Zn 1*. 在 0.05 水平（双侧）上显著相关。*. 在 .01 水平（双侧）上显著相关。观测发现，工业区土壤金属元素中，As含量与与Ni含量，Cd含量与Pb、Zn含量，Cr含量与Cu、Hg、Ni、Pb、Zn含量，Cu含量与Hg、Pb、Zn含量，Hg含量与Pb含量，Ni含量与Zn含量，Pb含量与Zn含量，达到了O01水平的极显著正

22、相关。相关系数分别是：0.690，0.829，0.754，0.920，0.902，0.698，0.675，0.695，0.983，0.670，0.622，0.612，0.634，0.739。As含量与Cr、Pb含量，达到0.05水平的显著正相关。相关系数是：0.380，0.395。这说明了，在工业区，无论是交通污染源内部和工业污染源内部，还是交通污染源与工业污染源之间，都存在着明显的交互作用；体现出在这一区域范围内，2种污染源对土壤重金属的共同作用。表八 山区中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 1-0.291*0.1130.527*0.0750.078-0

23、.205-0.176Cd 10.0660.0900.246*0.0490.766*0.606*Cr 10.364*0-.0060.945*0.1070.627*Cu 10.505*0.358*0.1220.252*Hg 1-.0450.2260.170Ni 10.0280.629*Pb 10.590*Zn 1*. 在 0.05 水平（双侧）上显著相关。*. 在 .01 水平（双侧）上显著相关。观测发现，山区土壤重金属元素中，As含量与Cu含量，Cd含量与Pb、Zn含量，Cr含量与Ni、Zn含量，Cu含量与Hg含量，Ni含量与Zn含量，Pb含量与Zn含量，达到了O01水平的极显著正相关。相关系数

24、是：0.527，0.766，0.606，0.945，0.627，0.505，0.629，0.590。Cd元素与Hg元素，Cu元素与Zn元素，达到了0.05水平的显著正相关。相关系数是：0.246，0.252。由于Cr与Ni、Zn主要来自于汽车运输，而其他元素没有极显著相关或着存在负相关关系。且有六种元素的实际值在背景值的区间里。故而可显示出，山区土壤污染源主要是交通运输，而受工业污染源的影响较小。表九 交通区中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 10.1210.1390.092-0.0040.228*0.0600.188*Cd 10.373*0.424*0.

25、211*0.351*0.615*0.294*Cr 10.894*0.0120.869*0.428*0.395*Cu 1.0320.886*0.506*0.432*Hg 10.0400.266*0.118Ni 10.396*0.503*Pb 10.482*Zn 1*. 在 .01 水平（双侧）上显著相关。*. 在 0.05 水平（双侧）上显著相关。观测发现，交通区土壤重金属元素中，Cd含量与Cu、Pb含量，Cr含量与Cu、Ni、Pb含量，Cu含量与Ni、Pb、Zn含量，Ni含量与Zn含量，Pb含量与Zn含量，达到了O01水平的极显著正相关。相关系数分别是 ：0.424，0.615，0.894，0

26、.869，0.428，0.886，0.506，0.432，0.503，0.482。As含量与Zn含量，Cd含量与Hg含量，达到了0.05水平的显著正相关。相关系数是 ：0.188，0.211。由于Cd、Pb、Cr、Ni、Zn元素主要来自于交通污染源，而且Cd与Pb显著相关较大，据此推测可能来自于工业污染源。Hg与其他的元素出现负相关关系，表明Hg存在着不同的污染源。表十 公园绿化区中8种元素的相关性 As Cd Cr Cu Hg Ni Pb Zn As 10.358*0.689*0.1070.1760.691*0.2650.285Cd 10.564*0.500*0.0540.433*0.598

27、*0.712*Cr 10.357*0.0230.739*0.397*0.509*Cu 10.1360.2670.756*0.521*Hg 1-0.0480.389*0.063Ni 10.1680.298Pb 10.748*Zn 1*. 在 0.05 水平（双侧）上显著相关。*. 在 .01 水平（双侧）上显著相关。观测发现，公园绿化区土壤重金属元素中，As含量与Cr、Ni含量，Cd含量与Zn含量，Cr含量与Ni含量，Cu含量与Pb含量，Pb含量与Zn含量，达到了O01水平的极显著正相关。相关系数分别是：0.689，0.691，0.712，0.739，0.756，0.748。As含量与Cd含量，

28、Cr含量与Cu、Pb含量，Hg含量与Pb含量，达到了O05水平的显著正相关。相关系数分别是：0.358，0.357，0.397，0.389。由于Cr、Ni、Cd、Zn、Pb来自于交通污染，且有以交通污染源为主的Cu、Pb、Zn之间存在着明显的交互作用，显示出居民区的土壤重金属受到交通污染源的控制，受工业污染源的影响较小。Ni与其他元素出现负相关关系，说明Hg还存在着其他污染源。Hg与其他元素没有相关性，其含量与公园绿化区的背景值最为接近，说明Hg它很有可能来源于自然环境中的岩石风化。331 土壤重金属污染评价综合以上信息分析得出：这八种土壤重金属含量分布比较均匀。将全部重金属元素与背景值相比较

29、知道均比背景值大出许多，此地存在中污染区和重污染区。在生活区，重金属元素比重差异不是太大(由权重表可知)，但含量却是不小，均大于背景值，但Zn比重最大，占了17.388 。在工业区，重金属元素比重差异很大，尤是Hg元素的比重最大，占了42.272 。通过资料查询知晓，Hg主要来自于工业废水，故而可知工业区的废水排放很不符合环境标准。观察浓度表，所有元素含量都大大超过背景值。在山区，各元素权重最为均匀，可是依然处于污染状态。在交通区，Hg元素的权重亦是最大，是41.573，来自于交通污染源的各元素交互作用明显，含量也是不可小觑，其值都超过了背景值。在公园绿地区来自于交通污染源的重金属元素和来自于

30、非交通污染源的重金属元素，或者各自内部之间都存在明显的交互作用，体现出不同污染源对土壤重会属分布的共同作用。由此对比该地区土地污染等级评价标准，可得生活区、工业区、交通区、公园绿地区都达到了重污染水平，而山区也是中污染水平。 综合八种重金属元素含量的大小，不难得出五个区的污染程度排序：工业区 生活区 交通区 公园绿地区 山区 。反映了重金属污染从地区内部到远郊递减的规律。34 模型二341 传播特征的分析本文对题中传播特征的解释是在这片地域的重金属污染物的分布特点和动态变化趋势。1.分布特点：对于这种元素，它的分布主要集中于工厂较多的工业园区，和道路交通较为发达的地区。对于这种元素，它的分布主

31、要集中于交通要道地区附近，并且有扩散的潜能，如在公园河山区土壤中也出现了这种元素，在从元素的空间分布图来看，工业区主要建在山区附件。因此，山区也有也存在了这种元素。类似，公园绿地和居民区都存在较多的重金属元素Cd。对于Cr这种元素，它的分布主要集中于交通地区，和工业园区。Cr主要用于合金(如供铬钢用的铬铁合金或如镍铬电阻合金)及电镀(如汽车减震器和内部装璜以及切削工具)，电镀和制造特种钢。而Cr的来源主要以铬铁矿或铬铅矿形式存在。因此，工业地区可能较多。对于Cu这种元素，它的分布主要集中于工业区，由于Cu的用途之广，所以该城区的污染浓度非常大对于Hg这种元素，它的分布主要集中于工业地区，以及对

32、于汞矿的开采与汞的冶炼。对于Ni这种元素，它的分布主要集中在交通，工业园区，分布面积广泛。对于Pb这种元素，它的分布同Cd，是重工业地区污染排放物之一对于Zn这种元素，它的分多存在于交通地区。2一般动态变化趋势 同一污染源，不同元素的迁移变化趋势是不同的。它是一个类似于的函数，即三维、浓度、距离关于平面弥漫速度的关系。在不考虑特殊的情况下，重金属浓度越高，它的扩散速度越大，扩散的距离就越远。同时考虑到地形影响的情况，在海拔高的地方，污染物迁移速度变慢，污染较轻（此处无污染源）。海拔低的地方，污染物迁移速度快。342 模型假设1城区的重金属污染是一种连续的稳定场2污染源以点污染源为主 3未知的污

33、染源的扩散方式是同心圆模式4污染物的浓度变化方式不受高度影响，并且以相同的梯度减小5 对于单个污染物浓度来说，它的最大点就是该地区的污染源6 不考虑高度对污染物扩散方式的影响343符号的约定 横坐标 纵坐标 各重金属的浓度 能量极大点 行走极限点 行走极限点到能量极大点的距离 浓度变化 颇高344 模型建立对第一问中，取As重金属污染浓度空间分布图中的一块图形进行定性和定量的分析。对第一问中，取As重金属污染浓度空间分布图中的一块图形进行定性和定量的分析。 由于污染源的不确定性，假设该地区的污染源为。在这一扩散过程中，选取浓度为背景值的点,再选取周围大于背景值的点,其中,为我们选取点中浓度比较

34、大的一个或某几个点。分别作出它们的线型关系。由于污染源是以等浓度梯度扩散的一种源，且不考虑高度的影响，因此当距离为时，这点的浓度就是。如下示意图所示：上图为扩散的理想模型，背景值点是某一元素达到一级标准的点，我们把这个点称为行走极限点，记为。当这个元素到达这个点的时候，它就不在向外扩散了。而位置浓度扩散点就是我们要找的污染源 现假设污染源的扩散能力极强，足以穿透的地面，我们把这些点称为能量极大点，记为。随着污染物的迁移，逐渐的增大，浓度逐渐变小，这是一个及其漫长的过程，我们把这种过程成为弥漫时间。 现在，由已知的一系列点对其做一维线性回归方程,称这个方程为的潜意识一维矢量方程，的过错称为记忆过

35、程。我们就知道了这个方程的任意的坐标。： 在问题一中的元素空间分布图中取任意元素的局部等浓梯度变化图，找出该元素的背景值浓度，以该点为点向比它浓度更大的地方做的潜意识一维矢量方程，其中,是有界域的值。任意带入该界域的值，称这对点为的原始最终点。当该点以最接近点的方式进行移动时。污染物浓度扩散的分布是均等分布的，因此它的浓度和上图所切割的理想面积成正比，点形成的浓度梯度角为（L以米计算）。当取时，,面积；当取时，,面积 当取时， ，面积面积和浓度的关系为点的浓度那么把分为个等分的段，那么每段的浓度为,某一点的浓度为到此我们确定了简单的一维弥漫方程：做法如下图：当确定了的潜意识一维矢量方程时，在界

36、域内给定x点，就可由简单的一维弥漫方程得到它的浓度，和相比：如果时,我们则认为过程不够彻底，并继续朝矢量方程的方向前进。直到时，我们则认为的记忆过程基本完成。若并确认这点为所要找的污染源，这种通过随机取点来找污染源的方式所找的污染源是不准确的点污染源，而已有一定半径的面污染源。模型运用现以As的部分图来说明：作出上图的一维矢量方程 任意给出多个点的坐标 ,带入一维弥漫方程 As的背景值点为，由浓度对比图看出该的横坐标为5200，并从5200开始，以直线向左前进，通过ab到达c得到，通过到达得到,最后计算出， ，,所以最后得到的浓度为。本文认为这点就为基本的污染源点。最后通过这点反求出的值为，带

37、入上面的矢量方程得到，所以污染源的坐标是。考虑到污染源的模糊性，以为半径左的模糊污染源。并在这个理想的模糊污源里用以上的方法继续寻找最终的理想点污染源。35 模型三城市地质环境演变模式：城市地质环境同自然地质环境的相同点在于它们拥有相同的水圈、大气圈和岩层圈，由于重金属元素对大气圈和岩层圈的影响几乎可以忽略，我们就认为城市地质环境的演变中，重金属元素只影响水圈。模型二中没有考虑到重金属元素的对于大地的渗透作用，实际上，重金属元素的渗透也是不可忽视的。如果他们过多被排放，就会透过渗透作用到地下水层，严重的影响了地下水层的水质。因此我们要控制排放量，通过减少排放量从而减少它对地下水层的污染。模型三

38、基于物力、物理化学、流体力学等基础学科的知识建立了非保守性污染物质在土壤中迁移转化一维控制方程。还需要收集的数据还要：Ss为比贮水量 Ci 污染组分i在地下水中的浓度 介质体积含水率，对于饱和土壤，则是介质有效孔隙度Dij 弥散系数张量vi 方向的实际渗流速度或平均孔隙流速，m/s； n个不同反应中第个反应的溶解污染组分i的产率xi 笛卡儿坐标 模型建立 在稳定均匀流及x轴与流速向量方向一致的情况下，如果不考虑给排水和生物降解作用，并忽略弥散系数与空间和时间的关系，可由 得到非保守性污染物质（与土壤介质发生吸附、沉淀或生化反应等作用）在饱和土壤中迁移转化的一维控制方程： 对于有机物污染问题，吸

39、附作用是影响污染物迁移转化的主要因素。如果把吸附和解吸完全看作一个可逆过程，则： 如线性平衡吸附方程能够很好地描述介质的吸附特性，那么上式中的可表示为： 代入计算整理后可得： 令，其中为阻滞系数，表示多孔介质对于非保守性物质的阻滞能力，则上式可简化为： 在一维、稳定流场和连续注入定污染物浓度（C0）时，方程的定解条件如下： 0 0 0借助于Laplace变换，并代入上述边界条件，方程的解析解为： 在10的情况下，式（2-19）中右端第二项比第一项小的多，可以忽略不计，这时污染物浓度在时空上的变化规律可以用以下近似解表示： 经过适当变量替换后，式（2-20）可以表示为： 式中和分别为余误差函数和

40、正态分布函数，即： 由上述求解可以看出，污染物浓度在时空上变化规律的近似解析解在形式上与正态概率分布函数相同，所以目前分析室内一维、稳定流场和连续注入定示踪剂浓度的实验数据，以确定多孔介质弥散系数的方法主要有正态分布函数法和反函数法。但是正态分布函数法存在明显缺点，而反函数法在考虑示踪剂和多孔介质之间的吸附与解吸作用时不能够应用。为此，本研究借用反函数变换思路，在利用保守性物质确定土壤介质的有效孔隙度基础上，通过相关计算确定非保守性物质BTEX在介质中的弥散系数D和阻滞系数Rd6。模型的评价与改进模型一的优点：在元素浓度分布图中，能直观的反映单个元素在该区的污染程度；在重金属污染的评价中，考虑到元素毒性的特点，建立判定元素毒性的权重矩阵，清晰的判断了该区的污染程度；在重金属污染程度的影响中，利用了Person相关性系数分析法，合理的说明了元素之间的相关性