补充2第12章全国地下水污染调查评价与综合研究20141104p73周金龙.pptVIP

将调查指标分类分级； 建立指标等级体系； 明确层级评价的主要内容； 梯级结果集成实现综合表达。 强调几点 调查：要控制主要污染源和污染类型。 评价：要深入解读污染发生发展过程， 深入研究背景，突破简单统计。 综合研究：以专题研究入手，典型深化； 通过系统编图，解读主成果； 集成数据库和各区成果。 I 西、北江三角洲冲积平原 孔隙水环境单元 II 东江三角洲冲积平原孔隙 水环境单元 III 广花平原区孔隙水环境 单元 IV 高要冲积平原孔隙水环 境单元 V 新兴-鹤山台地丘陵区裂 隙水环境单元 VI 潭江河谷平原、台地、 丘陵区裂隙孔隙水-裂隙 水环境单元 VII 四会低山丘陵基岩裂隙 水环境单元 VIII 从化-增城低山丘陵区 块状基岩裂隙水环境 单元 IX 东江谷地平原孔隙水环 境单元 X 东部山地丘陵区基岩 裂隙水环境单元 XI 海积平原环境单元 地下水环境单元划分 根据地下水的形成条件（如地形地貌、地质构造特征、含水层岩性等），全面考虑地下水赋存环境的径流条件、氧化还原条件、酸碱度等因素，对研究区进行地下水环境单元划分，将水文地球化学作用具有共同性的总体进行归类。 剔除异常值 异常值检验 方法 异常值 浓度 分位数 检验分布类型 分布类型 公式 正态分布 对数正态分布 偏态分布 以中位数为集中值，以百分位数为背景区间 背景值计算 横坐标为 对数值 背景值计算实例 地下水环境单元 检出率（%） 分布类型 背景集中特征值（mg/L） 背景值范围 （mg/L） Ⅰ 100 偏态 38.96 13.90-72.01 Ⅱ 100 偏态 29.24 8.69-56.65 Ⅲ 100 偏态 27.15 4.92-53.00 Ⅳ 100 偏态 26.32 5.26-38.61 Ⅴ 100 正态 27.73 11.63-43.83 Ⅵ 100 偏态 16.84 4.70-57.97 Ⅶ 100 正态 14.51 7.47-21.55 Ⅷ 100 对数正态 14.89 6.93-32.01 Ⅸ 100 偏态 20.18 3.51-43.87 Ⅹ 100 偏态 14.04 3.51-65.39 Ⅺ 100 偏态 28.32 6.22-68.31 以氯离子为例 (mg/L) 计算公式： K阶趋势面方程： 为常数，i=1/2(k+1)(k+2)-1 二、趋势面分析法 趋势面分析：利用数学模型来模拟（或拟合）地理 数据的空间分布及区域性变化趋势。 多项式函数polynomial.m function yhat=polynomial(bb,x) x1=x(:,1); x2=x(:,2); a0=bb(1); a1=bb(2); a2=bb(3); … yhat=a0+a1*x1+a2*x2+a3*x1.*x2+a4*x1.^2+a5*x2.^2; %二次多项式%yhat=a0+a1*x1+a2*x2+a3*x1.*x2+a4*x1.^2+a5*x2.^2+a6*(x1.^2).*x2+a7*(x2.^2).*x1+a8*x1.^3+ a9*x2.^3 %三次多项式 … End 拟合函数ninlinfit.m(为MATLAB自带函数，代码不再列出) 主函数surface fitting.m x=[输入]; y=[输入]; beta0=[1,1,1,1,1,1,1,1,1,1]; %给参数赋起始值，进行迭代计算 [beta,r,J]=nlinfit(x,y,@polynomial,beta0) %调用非线性最小二乘法拟合函数 yhat = polynomial(beta,x) %返回拟合值 r1=r; Sd=r1*r y1=yhat-mean(y); Sr=y1*y1 NC=Sr/(Sr+Sd) n=(输入); % 输入采样点数 K=（输入） % 输入除去常数项的多项式项数 F=(Sr/K)/(Sd/(n-K-1)) % 计算F统计量 运用国际通用计算软件Matlab编制趋势面分析计算程序。 二、趋势面分析法 三、预筛选法（Pre-selection method） 监测数据 初选原则： 样品离子平衡相对 误差不得超过10%； 取样深度明确； 含水层岩性明确； 不得为地热含水层、 高盐量含水层； 利用中值平均法消