全国大学生数学建模竞赛A题论文——城市表层土壤重金属污染分析

1、2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置

2、报名号的话）： 所属学校（请填写完整的全名）： 中国人民大学 参赛队员 (打印并签名) ：1. 于焕杰 2. 陈姝含 3. 毕雅枫 指导教师或指导教师组负责人 (打印并签名)： 日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响

3、也日益突出。另一方面，重金属污染是破坏土壤环境的重要因素，直接或间接危害到人体健康。本文以给定的城市为例，对土壤中的Hg、Pb、Cd、As、Zn、Cr、Ni和Cu这8种重金属进行分析，运用了统计分析方法、Muller地积累指数法、因子分析法、数据模拟等方法，研究了重金属在该城区的空间分布特征，进行相应的污染程度评价，在此基础上进行污染因素的定量分析，分析污染源并确定其位置，主要内容和结论如下：（1）通过相关数据处理软件对数据进行标准化、插值和拟合，以直观的分布图配合地积累指数评价得出8种重金属元素在该城区的空间分布，并对其污染程度做出评价。结果表明：大多数重金属污染主要集中在生活区，工业区和交

4、通区的交叉地段。Hg和Cu在整个城区的污染最为严重，Zn、Pb和Cd也为主要的污染成分。工业区的各项污染指标均高于其他区域。（2）因子分析法有效地将来自相似污染源的重金属元素归类为四个主因子，再对照重金属元素来源有效地分析重金属污染的原因。过程中，对因子分析法的可行性，因子的成份和污染因素都进行了有效、细致地分析。最后得出主因子1（Cr和Ni）和主因子2（Cd、Pb和Zn）主要来源于电镀厂、电子厂、塑料厂等工厂排放的废水和汽车尾气，同时因子2沿交通带分布；主因子3中Hg的广泛来源和易迁移性时其在整个区域广泛分布。（3）运用传播特征分析和因子分析的结果，估测出若干个污染源位置所在的范围区域，利用

5、SPSS软件拟合出适配的对数形式的复合函数，确定污染源的具体坐标。本文拟出了四个主因子的九个污染源的具体坐标。（4）通过对所建立的各模型优劣性的分析和整个土壤重金属污染分析过程的思考，引入污染物进入研究区的速率，建立预测研究区未来一定时期内各时刻土壤中指定污染物的含量的土壤环境质量的预测模型。【关键词】 重金属污染 Muller地积累指数法 因子分析法 数据拟合 MATLAB编程 目录摘要（1）关键词（1）1. 问题重述及分析 1.1 问题重述（3） 1.2 问题分析（3）2. 模型假设及符号说明 2.1 模型假设（4） 2.2 符号说明（4）3. 表层土壤重金属污染分析 3.1 任务一 3.

6、1.1 利用MATLAB绘图、插值和拟合（5） 3.1.2 利用Muller地积累指数评价土壤重金属污染 3.1.2.1 模型一的建立（） 3.1.2.2 模型一的求解分析（） 3.1.3 结论（10） 3.2 任务二 3.2.1 模型二：用因子分析法进行数据处理（10） 3.2.2 重金属污染的原因（12） 3.3 任务三 3.3.1 重金属污染物传播特征分析（） 3.3.2 模型三的建立和求解分析（） 3.3.3 模型三的误差分析（）4. 任务四：模型评价和拓展 4.1 模型评价 （13） 4.2 模型拓展 （13）5. 参考文献（14）6附件 6.1 附录一（14） 6.2 附录二（14

7、）6.3 附录三（） §1问题重述及分析1.1 问题重述随着城市的发展，人们对重金属资源的需求越来越大，而生产、加工的过程中产生的重金属废物也越来越多，渐渐导致土壤中重金属含量超标，这势必会对生态环境造成一定的影响和破坏。按照不同的土地的功能可以将城区分为五大功能区：生活区、工业区、山区、主干道及公园绿地区，不同的区域环境受人类活动影响的程度不同。现要求对某城市城区土壤地址环境进行调查。根据已给出的该城区表层土壤取样点的位置、海拔高度、所属功能区、重金属元素浓度及重金属元素背景值等数据，请通过数学建模的方法研究解决土壤重金属污染分析问题，需要完成一下任务：给出8种主要重金属元素在该城

8、区的空间分布，并分析该城区内不同区域重金属的污染程度；通过数据分析，说明重金属污染的主要原因；分析重金属污染物的传播特征，建立模型，确定污染源的位置；分析所建立模型的优缺点。为了更好地研究城市地质环境的演变模式，讨论我们还应该收集的信息点，以及利用信息如何建立模型解决问题。1.1 问题分析（1）为了讨论该城区内不同区域重金属的污染程度，我们首先需要对附件中所给的数据按功能区进行筛选和分类处理。然后根据数据，利用MATLAB软件绘制出取样点点位坐标和插值后整个城区各个重金属元素的污染评价，可以粗略地得出8中重金属元素在该城区的空间分布。若要进行较为准确的定量分析，还需要采用适合的指数评价方法，本

9、文中采用的是Muller地积累指数评价方法。（2）为了分析重金属污染的原因，即确定重金属污染物的来源是什么，需要充分利用采样点的各项数据指标，但是如果分别分析8中重金属污染物在五大功能区的污染源，不仅工作量大，造成污染来源的重复确定，而且也不符合该城区各个功能区相互贯穿的具体实际情况。所以，我们需要寻求某种可行的方法将关联性大的重金属元素归类处理。（3）通常情况下，污染物的传播都存在一定的规律，我们需要通过查找文献和观察分析采样点的数据找出重金属污染物传播特征，假设模型，然后通过对数据的分类，利用MATLAB软件拟合数据，得到能反映重金属污染物传播规律的函数，才能有效地确定污染源位置。（4）研

10、究城市地质环境的演变模式，我们初步设想为研究污染物含量随时间推移的变化。对于应该收集的信息，我们可能要在研究出以上几个求污染程度、污染源位置等问题以后才能具体地设计及建立有效模型。§2模型的假设与符号说明2.1模型的假设（1）本文中土壤重金属污染分析的研究，假设是在忽略地区的特殊性下进行，即忽略各个地区的差异，采用统一的评价标准。（2）假设调查的污染物浓度中只含有少数的高峰值，使数据的分布呈现高度的偏倚性。忽略重金属污染物空间分布的极大变异性和偏倚性。（3）2.2符号说明Igeo：Muller地积累指数；Cn：重金属元素n在土壤中的含量；Bn：重金属元素n的背景值；Ct：土壤中重金属

11、元素在t单位时间后的含量；C0：重金属的起始含量；V：每个单位时间内进入土壤的重金属含量；P：特定重金属在研究区土壤中的单位时间残留率；t：及为需要预测的时间变量；p、m：任务三中模型进行数据拟合时需要的参数；x、y：表示采样点的横、纵坐标；a、b：任务三中模型代表污染源横、纵坐标的参数；F：因子分析中指定因子的因子变量得分。§3表层土壤重金属污染分析3.1 任务一3.1.1利用MATLAB绘图、插值和拟合针对问题中附件1的相关数据，我们利用MATLAB软件绘制出该城区的三维立体地形图（如图3-1），让我们较为直观的了解该城区的地形分布情况，以便接下来各任务的分析。图3-1 该城区的

12、四维立体地形和8种重金属元素综合浓度指标图 由于在调查的过程中，是将所考察的城区划分为间距1公里左右的网格子区域，再按照每平方公里1个采样点对表层土壤（010厘米深度）进行取样、编号，所以根据问题中附件1所给的采样点坐标数据和分区情况，利用MATLAB制图，得到该城区五大功能区的模拟分布情况（如图3-2）：图3-2 该城区五大功能区的模拟分布情况接下来，为了粗略地由图示分析出8种重金属元素在该城区的空间分布和不同区域的污染程度，我们利用问题中附件2中的8种主要重金属元素的浓度，利用MATLAB进行插值和拟合后，得到各个重金属元素在整个城区区域范围内的污染评价（如图3-3）：AsCrCdCuHg

13、NiPbZn 图3-3 各个重金属元素在整个区域内的污染评价对照图3-2和图3-3，可以看到各个重金属元素在该城区各区域的表层土壤的富集情况，我们可以得到如下分析：在山区和公园绿地区，8种重金属元素污染程度都较轻；工业区的As、Cd、Hg和Pb含量均较其他区域要高；交通区的As、Cr、Cd、Cu、Hg和Zn污染评价较高；生活区的Cr、Cu、Hg和Zn也呈现很高的态势；Ni含量在整个区域造成的污染水平均较低。由此我们也可以验证知工业“三废”的排放和人们生活中产生的废弃物是造成土壤重金属污染的原因。3.1.2 利用Muller地积累指数评价土壤重金属污染3.1.2.1 模型一的建立以上由MATLA

14、B插值拟合的图让我们直观的了解8种重金属元素的污染情况，但是只能粗略地确定污染区域，并不能有效地分析不同功能区的污染程度。接下来，我们将建立模型进行较为准确的定量分析和测试。地累积指数法地累积指数法是德国海德堡大学沉积物研究所的科学家Muller于1969年提出的，在欧洲被广泛采用。其计算式如下： （I）试中：Cn是元素n在土壤中的含量；Bn是普通页岩中该元素的地球化学背景值（本文中用问题附件中已给出的背景值）；K为考虑各地岩石差异可能会引起背景值的变动而取的系数（一般取值为1.5，本文即取定为1.5）。Igeo值为Muller地积累指数。分级标准根据Igeo值将污染等级分为6级，对应污染程度

15、为无污染至极强污染（如表一）：Muller地积累指数分级标准地积累指数Igeo5<Igeo<104<Igeo<53<Igeo<42<Igeo<31<Igeo<20<Igeo<1Igeo<=0分级6543210污染程度极严重污染强-极严重污染强污染中等-强污染中等污染轻度-中等污染无污染 表一 Muller地积累指数分级标准该城区不同功能区表层土壤重金属平均浓度根据对该城区不同功能区表层土壤采样中重金属浓度的分析测试， 通过计算分析，得到各个功能区土壤重金属的平均浓度（如表二）：各区域重金属平均浓度 (g/g)AsCd

16、 CrCuHg NiPbZn生活区6.270.2899669.0249.40.0930418.3469.11237.01工业区7.250.3931153.41127.540.6423619.8193.04277.93山区4.040.1523238.9617.320.0409615.4536.5673.29交通区5.710.3600158.0562.210.4468217.6263.53242.85公园绿地区6.260.2805443.6430.190.1149915.2960.71154.24 表二 五大功能区重金属元素平均浓度土壤重金属元素背景值问题附件中已给出不同土壤重金属元素的背景值（如

17、表三）：8种重金属元素的背景值 (g/g)AsCdCrCuHgNiPbZn3.60.133113.20.03512.33169 表三 8种重金属元素的背景值3.1.2.2 模型一的求解分析利用表二的各个功能区的重金属平均浓度和计算式（I），可以得到各个功能区中8种重金属浓度的地积累指数（见表四）：不同功能区土壤重金属浓度Muller地积累指数As Cd Cr Cu Hg Ni Pb Zn 生活区0.21550.57240.56981.3190.8255-0.00860.57171.1953工业区0.4251.01150.19992.68743.6130.10261.00061.4251山区-0

18、.4186-0.3564-0.2552-0.1931-0.3581-0.256-0.347-0.4979交通区0.08050.88460.32011.65163.0893-0.06640.45021.2304公园绿地区0.21320.5247-0.09160.60861.1311-0.2710.38470.5755表四 五大功能区土壤重金属浓度的Muller地积累指数将表四与表一的分级标准对照，我们即可得出不同功能区土壤重金属浓度的Muller地积累指数评价（如表五和图3-），进而分析重金属元素在该城区的空间分布和各个城区的污染程度。不同功能区土壤重金属浓度Muller地积累指数评价AsCdC

19、rCuHgNiPbZn分级分级分级分级分级分级分级分级生活区11121012工业区12134122山区00000000交通区11124012公园绿地区11012011表五 五大功能区土壤重金属浓度的Muller地积累指数评价图3-4 五大功能区土壤重金属浓度的Muller地积累指数评价柱状图分析从以上的表格和柱状图3-4中可以看出：在山区中，8种重金属元素含量均较低，无污染情况；在其余的四个功能区中，Hg污染非常严重，在工业区和交通区中均已经达到了强污染；Ni在各个功能区中含量很低，只有工业区中的有轻度-中污染；As、Cd、Cr、Pb、Zn在四个功能区中均有轻度-中或者中污染；而Cu污染也较严

20、重，且在各个功能区中差异较大，Cu浓度大小为：工业区（127.54g/g）>交通区（62.21g/g）>生活区（49.4g/g）>公园绿地区（30.19g/g）>山区（17.32）。3.1.3 结论按照Muller地积累指数评价法，8种重金属元素在不同的生态功能区中的污染水平表现出了较大的差异。如在工业区中，8种主要的重金属元素均存在不同程度的污染，而山区则暂时没有受到重金属的污染；另一方面，不同的生态功能区，8种重金属元素也呈现出较大的差异。比如，各个功能区中Ni污染几乎不存在，而Hg、Cu污染严重，有的已达到强污染水平，污染问题严重，不可忽视。同时，我们看到的，工业

21、区中Hg为强污染，Cu达到中等-强污染水平，Cd、Pb、Zn达到中等污染水平；交通区中Hg、Cu、Zn污染严重，As、Cd、Cr、也有一定污染；生活区中Cu和Zn污染严重；山区土壤环境良好；该城区表层土壤中Ni含量很少这与前面由MATLAB软件插值得到的污染评价是相符的，很大程度上验证了我们模型的适用性和准确性。3.2 任务二 模型二：用因子分析法进行数据处理因子分析因子分析是由Charles Spearman在1904年首次提出，其主要目的是将具有相近的各个变量置于一个公因子之下，即将相关比较密切的几个变量归在同一类中，每一类变量就成为一个因子，以较少的几个因子反映原资料的大部分信息。因子分

22、析过程该城区表层土壤取样点重金属含量的数据特征完全符合因子分析的要求，本文中，我们以Hg、Pb、Cd、As、Zn、Cr、Ni和Cu这8种重金属元素指标做因子分析。这样，我们在解释该城区重金属污染的原因时，就可以通过讨论综合因子来解释各个指标变化异常时的情况。我们首先对该城区表层土壤取样点重金属元素含量的数据进行标准化处理，标准化的公式为Xij =（XijXj）/j，其中Xij为第i个样本的第j个指标值，而Xj和j分别为j指标的均值（本文中为所给的背景值平均值）和标准差。标准化后的数据见本文附录一。接下来，我们运用SPSS11.0软件进行因子分析，分析结果如下（如表六、表七）（详细的因子分析过程

23、见附录二）：成份初始特征值旋转前旋转后合计方差的 %累积 %合计方差的 %合计方差的 %累积 %13.5644.544.53.5644.52.09726.21226.21221.1514.37758.8771.1514.3772.06725.83552.04730.96512.06370.9410.96512.0631.25815.72167.76840.7689.59680.5370.7689.5961.02112.76980.53750.5787.2287.75660.4325.39993.15670.3013.76996.92480.2463.076100表六 特征值和累计贡献率旋转前旋

24、转后项目F1F2F3F4F1F2F3F4As_0.426-0.20.6810.5510.1340.1730.0260.969Cd_0.7110.2810.282-0.3220.1310.8510.1310.096Cr_0.735-0.444-0.303-0.0460.8780.2420.0340.009Cu_0.7560.125-0.3650.1370.5770.330.546-0.025Hg_0.4080.673-0.2970.4490.0030.1460.9420.035Ni_0.723-0.515-0.190.1370.8770.1540.0360.218Pb_0.7640.3140.2

25、37-0.2480.170.840.2280.117Zn_0.699-0.0370.123-0.2410.3990.6280.020.096表七 因子载荷矩阵利用因子分析所得到的4个因子可作出各个因子在空间分布的等直线图，能更直观地说明各个因子在空间平面上的分布特征（如图3-5）：图3-5 主因子在平面空间的等值线分布图 根据因子载荷矩阵，主因子1的成分为Cr、Ni的组合，这说明这两种这两种重金属污染物可能是来自同一来源或者相似的来源；主因子2是Cd、Pb和Zn，表明这几者来自相似来源；主因子3是Cu和Hg；主因子4是As。主因子2和3的方差累计贡献率达到41.556%，这说明Hg、Cu、Zn

26、、Cd和Pb五种重金属的方差累计贡献率要超过第一主成分Cr和Ni。这可以解释Hg、Cu、Zn、Cd和Pb为该城区土壤主要的重金属污染元素，这与前文验证的结果是完全一致的。3.2.2 重金属污染的原因土壤中重金属的来源是多途径的，总体上分为自然来源和人为因素。自然因素主要是成土母质本身含有重金属，不同的成土母质和形成过程所形成的土壤中重金属污染含量差异较大。前人的文献中总结了我国土壤重金属污染的主要来源（如表八）：（加文献：陈怀满）中国土壤重金属的主要来源来源重金属矿产开采、冶炼、加工排放的废气、废水和废渣Cr, Hg, As, Pb, Ni, Mo煤和石油燃烧过程中排放的飘尘Cr, Hg, A

27、s, Pb电镀工业废水Cr, Cd, Ni, Pb, Cu, Zn塑料、电池、电子工业排放的废水Hg, Cd, Pb, Ni, Zn汞工业排放的废水Hg染料、化工制革工业排放的废水Cr, Cd汽车尾气Pb农药、化肥As, Cu, Cd表八 中国土壤重金属的主要来源通过因子分析结果参照五大功能区的模拟分布情况及土壤重金属的主要来源，认为该城区土壤中重金属污染物可能来源于铜汞等重金属冶炼和加工排放的废气、废水和废渣，农业活动过程中施用的农药、化肥，以及汽车尾气。从前面的主因子的平面空间分布图可以看出，主因子1（Cr和Ni）主要分布于生活区、交通区和工业区交叉地段，可能来源于钢铁厂，电镀厂的工业废水

28、和居民使用的废旧电池；主因子2（Cd、Pb和Zn）主要分三分布在交通带和工业区上，说明Cd、Pb和Zn来源于相似的污染源，主要来源于电镀厂、电子厂、塑料厂等工厂排放的废水。而Pb的主要来源是汽车尾气，该因子沿交通带分布正说明了这一点；主因子4（Hg）在生活区、交通区、工业区和公园绿地区中都有较大范围的分布，前文的污染程度分析中也有说明这一点，因为天然本地情况下，汞在大气、土壤和水体中均有分布，汞的迁移转化也较易在水、陆、空之间发生。它主要来自氯碱、塑料、电池、电子等工业排放的废水。3.3 任务三§4 任务四：模型评价和拓展4.1 模型评价41.1优点:（1）. 本文中各模型实用性好，

29、理论性强，不仅采用直观的立体、平面图形，还有较为准确的定量分析的表格、模型。且各模型间能够相互辅助论证和有效验证。（2）. 在任务一中，我们采用Muller地积累指数法进行污染评价，地积累指数很好地对该城区土壤重金属污染程度进行较为全面的评价，且充分注意了人为活动对环境的影响，非常有助于我们对该区重金属污染分布、程度和成因的研究。（3）. 在任务二中，用因子分析法建立的模型有效地归类出相关性较大的重金属成分，避免了分别讨论各个重金属元素所需要的繁杂过程，也更符合重金属污染的实际情况。（4）. 缺点：（1）.理论上，地积累指数计算式中的K值，是为了考虑地岩石差异可能引起背景值的变动而取的系数。但

30、是我们的模型中用已给定的背景值，可能造成背景值与系数不是最匹配的，影响该城区重金属污染整体分级程度的判断。（2）. 对于该城区内各种污染物的来源，还需进一步补充该城区的现场实际调查和监测分析材料，才能进一步验证上述推断的合理性。（3）.在污染源的位置研究上，本文的研究大多采用统计分析方法结合在一些被怀疑为污染源的附近采点研究，没有非常有效地考虑空间位置，不能扩展为大区域范围内的处理。模型的假设不够充分，对模型的检验也不够详细。4.2模型拓展对城市土壤调查采样可以帮助我们了解所研究城市的地质环境现状，但如果是进行高频率重复采样和分析，不仅费时费力，从经济上考虑也是不现实的。所以，为了更好的研究城

31、市土壤重金属污染水平，研究地质环境演变模式。在分析出城市如今的重金属污染程度和污染源状况后，我们还需要了解重金属污染物的传播、迁移规律，以及各种污染物在未来一段时间轴上的变化趋势，以便我们及时采取防治措施。假设我们想要知道某个污染场址未来一定时间内各时刻土壤中指定污染物的含量，或者污染源附近某个还未受污染的场址多长时间内会受影响，在已知重金属起始含量和距离污染源距离等情况下，必须还要有该场址重金属污染物含量时间变化的量，也就是污染物进入该场址的速率（V）和重金属元素在土壤中降解和残留的残留率（P）。然后利用预测模型推导（计算式如下）： （III）式中：Ct是土壤中某重金属元素在t单位时间后的含

32、量；C0是该重金属的起始含量；V为每个单位时间进入土壤的量；P为该重金属在该地土壤中的单位时间残留率；t及为需要预测的时间变量。、 而经过查阅资料，我们知道：重金属污染物进入土壤后，通过土壤对污染物质的物理吸附、过滤阻留、胶体的物理化学吸附、化学沉淀、生物吸收等过程，使其不断在土壤中积累。此外，土壤重金属污染物质可通过稀释、扩散、分解以及植物吸收等作用而净化。总之，在土壤中重金属污染物的累积和净化是同时进行的，两者处于一定的相对动态平衡状态。所以，在获取V和P的值是，还需要考虑这些信息，仍然需要一个复杂的定量分析的过程。值得注意的是，该模型要求试验的环境无突变情景，且不考虑重金属元素随作物被带

33、出该场址土壤、重金属向土壤下层迁移等情况，只用P计算重金属在土壤中的降解。§5参考文献1韩中庚.数学建模竞赛-获奖论文精选与点评M.北京：科学出社,2007.2王文波.数学建模及其基础知识详解M.武汉：武汉大学出版社，2006.3王庚，王敏生.现代数学建模方法M.北京:科学出版社,2006.4赵静,但琦.数学建模与数学实验M.北京:高等教育出版社,2007.5姜启源,谢金星,叶俊.数学模型M.北京:高等教育出版社,2003.6冯杰,黄力伟,王勤,尹成义.数学建模原理与案例M.北京:科学出版社,2007.§6 附录6.1附录标准化数据编号AsCdCrCuHgNiPbZn14.

34、710.791.482.0428.881.550.730.2422.592.381.583.25280.610.812.431.4416.854.1659.617.621.8910.8933.5743.292.410.256.8163.04-0.2518.272.9953.0613.173.1528.9895.632.0823.1646.93611.6432.14.1182.06125.634.1867.364.1275.934.667.28.7810.751.455.326.9886.6931.2128.29698.691683.137.7458.4496.3594.2333.136.353

35、8.461995.633.5523.5661.27105.6929.835.1816.1211.71.8632.7394.31112.5914.64.929.1612.53.0313.4236126.195.231.449.0516.631.955.2911.01132.3662.826.6412.91.094.776.11140.9920.837.6523.53263.3211.6318.17152.1221.624.8522.0236.881.9513.220.19168.7230.487.0934.3820.882.6373.5838.07172.8215.117.1461.25152.

36、3211.0817.29184.717.262.9623.3222.251.926.6417.02194.233.080.017.0313.44-0.9616.274.64205.4416.1379.2732.615.135.3932.9467.46212.1218.596.0616.5720.751.7612.0323.88226.9336.1298.87375.461034.2625.08124.96236.436.5315.725.475.885.059.9825.7424-0.332.742.137.0613-0.4568.37250.54-0.01-0.46-0.860.71-1.1

37、9-0.36-0.06262.8210.515.7442.8313.621.6113.3722.91271.97.881.2614.2919.131.346.028.12280.16.21.65.792.50.085.319.12920.39.824.7112.92185.634.088.787.623016.4216.677.328.2276.255.9216.56111.22317.716.843.6724.4919.384.2121.9439.0132-0.118.8811.9312.7116.883.1310.13201.75333.06-0.62-0.2412.7721.48-0.1

38、6-0.5-0.53342.1212.886.7721.0719.251.9715.1230.41350.997.6525.2430.575.250.1612.5856.0736-0.1111.119.0216.5910.370.474.938.88372.124.242.634.398.620.249.4511.02380.773.43-0.27.45209.91-0.661.271.89390.7710.644.0210.68350.8912.6616.12403.293.131.013.33114.370.820.213.454114.4210.977.099.44233.1326.72

39、10.44424.236.915.1465.466.872.0511.4310.07432.590.541.562.596.371.290.861.83441.2114.6230.5615.89108.131.188.119.12451.4410.31-0.2114.279.25-0.457.2210.67461.90.731.83.1213.131.321.081.68470.7713.22.718.957.25211.6121.49482.123.282.034.93.372.792.492.36491.2115.634.2330.2710.751.9711.9525.575042.822

40、.157.4915.132.242.7923.93511.915.691.528.8620.845.818.38521.447.111.2212.675.630.398.0114.77531.445.792.7814.633.537.39.95263.7540.327.140.489.638.130.7417.125.14550.16.460.465.961.880.425.7513.1556-0.549.542.999.5716.5-0.1114.212.36570.326.10.3173.512.50.456.2116.19580.32-0.112.11-0.132.321.68-0.27

41、0.0659-0.773.44-0.682.346.34-1.042.372.2260-0.330.050.55-0.531.11-0.33-0.16-0.5611.445.463.2510.512.250.087.4914.96620.54-0.1-0.383-0.63-0.14.411.1963-0.774.11.6616.21-0.360.5525.5264-0.984.961.385.393.66-0.633.734.7765-0.5420.712.4611.28-0.89-0.484.886.1466-0.3311.932.245.930.31-0.393.944.0767-0.33

42、12.450.864.790.17-1.343.792.42682.820.451.192.66.740.051.19-0.24690.10.693.181.61-0.093.341.241.43700.990.270.396.211.920.08-0.93-1.1871-0.54-2.02-0.9-1.57-1.03-1.22-1.49-1.14724.962.131.25.0622.280.342.424.17730.17.250.751.88-1.03-0.63.910.87740.14.680.2412.4-1.160.325.173.5175-0.110.860.26-0.310.9

43、9-0.25-0.49-1.2576-0.98-0.27-1.22-1.52-1.3-1.11-0.62-0.79775.444.381.663.232.051.52.071.1178-2.032.5215.824.060.8516.244.1211.4979-1.1911.290.675.231.25-0.045.717.86800.11.010.813.58-1.160.521.520.99812.823.231.2414.961.791.060.551828.21-0.672.7915.522.992.211.31-0.41833.0613.42.9519.619.51.897.0816

44、.38429.4820.4623.61-0.891.416.968.07850.328.131.5710.2700.714.526.3186-0.772.56-1.5-0.33-0.49-1.793.80.7387-1.822.13-1.440.83-1.03-1.41.551.83883.762.621.114.264.060.7610.471.1189-0.774.50.524.5-1.3-0.814.923.85902.5925.221.34.141.112.3310.616.73912.592.430.551.02-2.50.582.471.46924.237.1914.2716.06

45、6.251.539.7222.28930.77-1.15-0.951.16.34-0.57-0.51-0.13941.211.260.241.63.120.970.41-0.95952.361.162.054.1111.430.341.110.51963.5212.090.766.2313.440.2913.853.57970.991.690.691.53-0.89-0.10.91-0.3698-0.546.57-0.562.93-1.16-0.665.474.9399-0.777.171.0538.61-1.430.854.517.271000.7714.593.259.851.3715.469.931014.48-1.431.650.45-0.630.61-0.29-1.441022.59-0.041.840.86-2.761.410.421.62103-1.45.471.3-0.265.540.858.742.18104-0.777.990.5910.91-1.56-1.131.841.081052.362.125.483.