基于MAPGIS的分形方法确定化探异常

摘要

地球化学元素的异常下限值确定是地球化学中重要的问题之一,目前还没有一个令人满意的具有科学依据的计算方法.传统的化探异常下限值计算是基于元素的地球化学分布呈正态分布或元素含量在空间上呈连续的变化这一假设为基础的,而事实上地球化学元素含量的空间分布是极其复杂的,研究表明,地球化学景观可能是一个具有低维吸引子的混沌系统,元素的地球化学背景值和异常具有各自独立的幂指数关系,由此导致了一种多重分形分布,因而可以利用元素的分形分布求出其异常下限.利用元素的分形分布求异常下限的几种常用方法均是以求取不同尺度r下,对应的N(r)数.如果用手工统计方法计算,计算过程简单但繁琐,地理信息系统(GIS)可以实现图形单元的动态查询和属性统计,这种特性很适合于统计不同尺度r下对应的N(r),因此,可以将分形方法与地理信息系统加以结合,在GIS平台上对地学图形和属性信息进行分析及统计处理,将原来比较繁杂的分形计算操作变得方便简捷,易于实现.本文以武汉中地信息公司开发的MAPGIS地理信息系统为例,对河北某地1∶5万Cu元素化探数据进行了处理,具体说明了以MAPGIS为工具,利用分形方法求取元素化探异常下限的过程,并将计算结果与传统异常下限计算方法所得的结果进行了对比,所得的异常下限值相近或一致,但由于传统的计算方法不仅要求对原始数据进行预处理,而且在下限值的确定上需要结合本地区的实际地质情况确定,因此人为干扰因素较大,而利用基于MAPGIS的分形方法计算元素异常下限值不仅操作简单,而且计算时由于不受元素特高值的影响,因而不用对原始数据进行预处理,人为干扰因素较少.通过对比得知该方法在实践工作中可行.