基于LMI的模糊增益调参控制器的设计.docVIP

黄金玉，等：地下水调查点特殊分布的评价方法解析—以长江三角洲地区为例 黄金玉，等：地下水调查点特殊分布的评价方法解析—以长江三角洲地区为例 ——————————————— 收稿日期：2016-4-26；修回日期： 基金项目：中国地质调查局项目“长江三角洲地区地下水污染调查评价”（1212010634400）。 作者简介: 黄金玉(1982-),女,本科,高级工程师，主要从事GIS应用与开发、水文地质环境方面的研究。E-mail：windtones@163.com 地下水调查点特殊分布的评价方法解析 ——以长江三角洲地区为例 黄金玉，姜月华，苏晶文，张泰丽 (南京地质调查中心 江苏 南京 210016) 摘 要：分析现使用GIS软件网格等值线成图特点，根据调查评价点分布不均匀的特征，结合Kriging网格等值线分析方法，对Kriging等值线绘制方法和运算函数进行解析，总结出更精确更适用的分析方法，对调查点具有不均匀特征的地区进行地下水污染程度及变化趋势的综合评价。依据分析结果开发GPSE-Data系统，使分析步骤更简便快捷运用最终结论和方法，以长江三角洲地区为例，进行地下水污染评价。 关键词: 地下水;分析评价; Kriging;长江三角洲;GPSE-Data系统 中图分类号:X143 　　　 文献标志码: A 分析方法现状 使用GIS软件对调查区进行评价，现有的等值线分析方法主要针对平均分布（如1km*1km）且分析范围为规则图形，但在做地下水污染调查时，无法做到均匀采集调查点数据，而调查区域范围也需要根据项目具体要求和结合其他调查和工作手段。因此，如果在使用不均匀分布调查点和不规则调查区域范围内仍然使用原始等值线分析方法，就会使分析结果与实际情况有较大偏差。所以，要根据此种情况对分析公式和计算参数进行分析和修正，使结果更符合实际结果。 评价方法分析 现有评价方法分析 等值线分析方法是Kriging分析方法，参与计算的参数包括横坐标数组、纵坐标数组以及分析指标值数组， 进行单指标质量评价，计算函数是： _tnKrigingGridingToBuf(double * tmpXPtr , double * tmpYPtr , double * tmpZPtr , long allDatNum , DemInfoStrcT * mHdInf , SearchDatStrcT * PtSdat , KrigingParmStrcT * PtKrgParm , double * tZptr)。其中，函数中参数 double *tmpXPtr,*tmpYPtr,*tmpZPtr及long allDatNum所代表的意义分别是对应离散数据的横坐标X值、纵坐标Y值、需要进行质量评价的单指标Z值以及总点数。 均匀分布调查点分析时，点数只作为总数量进入计算，但是针对不均匀点，点数就要用于计算，而调查点数量是变化的，无法提前预判，也不可能固定为某值。 另外函数中DemInfoStrcT *mHdInf 是网格化参数，DemInfoStrcT是网格高程数据信息结构.，这个数据结构是 typedef struct { long nx,ny; double xmin,xmax; double ymin,ymax; double zmin,zmax; } DemInfoStrcT (其中：zmin代表X向步长,zmax代表Y向步长)，如针对平均分布调查点，调查点间隔是等值，则X和Y的两方向步长就是调查时密度。但是，如调查点数据不是均匀分布的，各点之间间距是不均匀的，有的间距超过了均值，有的间距很小，如采用默认步长，则对间距较大的数据间数据就可能无法参与分析，结果会与实际情况有很大出入。 新研发的系统采用方法 由于分析调查点总数变化且需要保存参与计算，因此，在新开发的GPSE-Data系统中针对不同数量值的变量上选择动态数组方式，根据分析的数据量定义存储变量的大小，这种处理方法就合理的解决了两个问题。 计算关键参数X和Y向步长，通过分析后从向步长值出发，采用的是用网格线数来推算向步长的方式。由向步长（即网格间距）推算出网格线数值，通常来说网格线越密集越接近调查点区域范围，所以如果通过循环获取调查点ID值以及坐标值，计算出各个相邻点间距值，最终计算出所有间距的平均值设定为向步长的默认值，而网格线数值就分别通过公式计算后设置，算法是： mHdInf.nx=(mHdInf.xmax-mHdInf.xmin)/mHdInf.zmin和mHdInf.ny =(mHdInf.ymax-mHdInf.ymin