地球物理数据网格化方法的选取 - 图文-.docVIP

第34卷第1期物探与化探Vo.l34,No.1 2010年2月GEOPHYSI CALGEOCHE M ICAL EXPLORAT I ON F eb.,2010地球物理数据网格化方法的选取 刘兆平1,2,杨进1,2,武炜3 (1.中国地质大学地球物理与信息技术学院,北京100083;2.教育部地下信息探测技术与仪器 重点实验室,北京100083;3.中国科学院地质与地球物理研究所,北京100029 摘要:绘制地球物理数据等值线图时,需要根据客观环境特征和数据本身的特点,选择合适的网格化方法。通过实例,分别介绍了加权反距离插值法、克里格法、最小曲率法、最近邻点法、多项式回归法、径向基函数法、带线性插值的三角剖分法等常用网格化方法的选取方法、适用范围及参数设置等使用技巧。 关键词:地球物理数据;等值线图;网格化方法;克里格法 中图分类号:P631文献标识码:A文章编号:1000-8918(201001-0093-05 在地球物理数据处理中,等值线的绘制是重要的组成部分。通过各种勘探手段获取的数据要客观地描述或表征研究对象,须进行数据网格化。而不同的网格化方法有不同的效果,因此,结合实际资料特征,分析和研究不同的网格化方法、适应范围、及参数设置,对于客观、正确地描述研究对象有重要的意义。 在众多的商业化绘图软件中,美国GOLDE N软件公司的Surf er软件,以其方便、直观、快捷、安装简单、对系统要求低等优点得到广大地球物理工作者的青睐。利用Surf er强大的插值功能,根据不同的数据特征,科学地选择插值方法,灵活地进行参数设置,方便快捷地绘制各种常用的等值线图件,不仅可以提高工作效率,避免人为误差,改善技术成果的质量,还可被地球物理工作人员用于对特定研究对象的多方位分析,准确合理掌握研究对象的数据特征及其所反映的地质地球物理概况。 1网格化方法的特征及应用条件 绘制地球物理数据等值线图时,首先需要对离散数据进行网格化处理。所谓网格化是指通过一定的插值方法,将稀疏的、不规则分布的数据插值加密为规则分布的数据,以适合绘图的需要。Surf er绘图软件中提供了多种网格化方法,包括:加权反距离法(i n verse distance to a po wer、克里格法(K ri g ing、最小曲率法(m ini m um c urvature、改进谢别德法(modified Shepards method、自然邻点法(natural neighbor、最近邻点法(nearest ne i g hbor、多项式回归法(polyno m ial regression、径向基函数法(rad i a l basis f uncti o n、带线性插值的三角剖分法(triangu la2 tion/li n er i n terpolation、移动平均法(moving average 、数据度量法(data metrics和局部多项式方法(l o2 cal pol y no m ia l[1-2]。同一组数据采用不同的网格化方法,产生的绘图效果不同。因而根据研究内容的环境特征和离散数据原有的特点,选择合适的网格化方法,是正确研究、分析目标对象的基本保证。 克里格法最初是由南非金矿地质学家克里格根据南非金矿的具体情况提出的计算矿产储量的方法:按照样品与待估块段的相对空间位置和相关程度来计算块段品位及储量,并使估计误差为最小。后来,法国学者马特隆对克里格法进行了详细的研究,使之公式化和合理化。克里格方法的基本原理是根据相邻变量的值(如若干样品元素含量值,利用变差函数所揭示的区域化变量的内在联系来估计空间变量数值。该方法总是尽可能地去描述原数据所隐含的趋势特征,以区域化变量理论为基础,以变差函数为主要工具,在保证研究对象的估计值满足无偏性条件和最小方差条件的前提下求得估计值。例如,对于高值数据点会使之沿某一/脊0分布,而不围绕该点孤立插值,不形成/公牛眼0等值线。克里格法极为灵活,广泛地应用于各个科学领域,适于各种类型的离散数据,网格化精度高,是极佳的网格化方法。但该方法随着数据量的增大,计算速度较慢,因此最适合数量小于250个点数据的网格化,对 收稿日期:2009-05-13 基金项目:国家自然科学基金项目物探与化探34卷 于250~1000个数据点,效果也不错。换言之,克里格插值法是一种特定的滑动加权平均法。 最小曲率法是目前地学界流行的一种网格化方法。它采用迭代的方法逐次求取网格节点数据,其插值面类似一个薄的、线性弹性形变板,该/板0经过所有的数据点,且每个数据点具有最小曲率。由于最小曲率法采用全区的数据进行网格化,因而比较适合于数据分布不均匀的情况。在尽可能体现原数据的同时,最小曲率法产生尽可能的光滑曲面,绘制的图件比较