勘探重力实验matlab,MATLAB在重力图制作中的应用

1引言地球重力场是反映地球物质分布特征的物理场,制约地球及其空间任何物体的运动,与空间技术发展密切相关,是建设数字地球或数字中国的基础物理场信息。关于重力场的表现和描述主要有地球重力场模型、大地水准面、重力图(重力异常图)等,重力图是地球物理勘查工作的一项基础性成果,可供地质、石油、煤炭等部门和测绘系统使用。MATLAB主要由MATLAB主程序、Simulink动态系统仿真和MATLAB工具箱(Toolbox)三大部分组成[1]。MATLAB的一个主要优点之一就是可以使计算结果可视化,而重力图是重力场变化的主要表现之一,是重力场异常的可视化成果,所以应用MATLAB的这个特点可以使计算得到的重力异常值可视化得到重力图。本文在讨论了重力图的一些知识和MATLAB的特点后将MATLAB和重力图的制作结合起来。2重力图的相关知识重力图即重力异常图,是重力异常等值线图。如果在测区内进行了若干点的重力测量,则可以按式:(g0-r0)空=g+1g-0(1)式中(g0-r0)空表示空间重力异常,g为地面点的重力值,1g为空间改正,0为正常重力值。计算各点的空间重力异常,然后将重力点按其经、纬度展绘到图纸上,将图上空间异常值相等的各点联成曲线,称为等值线,这样的图称为空间异常图[2]。如图1所示。重力异常是由于地球质量分布不均匀造成的重力场中各点的重力矢量g和正常重力矢量0的数量之差。即:g=g-0(2)重力异常反映了地壳内部物质的分布情况,可供间接了解地壳的内部构造。反映区域地质构造的异常称为区域重力异常,反映局部地质体的异常称为局部重力异常。因此重力图是进行重力勘探的主要工具之一。区域重力调查的数据、成果及其编制的重力图件,是图1某区域的空间重力异常图行地球物理研究的基础资料,又有长期利用的价值。这区域重力资料可以用于研究沉积层到上地幔的结构和构特征,划分大地构造单元,研究区域重力场与矿产分布关系,划出成矿远景区,为探讨成矿规律,进行矿产预提供地球物理资料;区域重力资料,也可为研究地球形、发展空间技术、天然地震的预报等提供重要的重力场息[3]。因此,编制重力图尤其是大比例尺的重力图是一非常有意义的工作。目前主要的重力改正有空间改正、布格改正、地形改正、壳均衡改正等,同时进行了这些改正得到对应的重力异常分为空间重力异常、布格重力异常、地改后的重力异常和均衡力异常,同时有这些异常得到了相应的空间重力异常图、布重力异常图、地改后的重力异常图和均衡重力异常图。用MATLAB制作重力图在一些试验、工程测量及科学计算中,经常会得到大的离散数据点,称作采样点。在对这些数据利用之前,对其进行分析和处理,因为对于大量的原始数据,很难接从中找出内在的规律,而把这些数据用各种形式的图表示出来,即把数据所反映的很多内在规律直观地展示人们面前。因此,数据可视化是人们研究科学、认识世不可缺少的重要手段。正是基于上面所述的需要,MATLAB在数据可视化方提供了强大的功能,它可以把数据用二维、三维乃至四图形表现出来。对于定义在x-y-z坐标系上的四维可视,MATLAB的色彩控制和动画等命令可为四维甚至更高的表现提供现实手段,即主要用色彩、切片图和等位线进行四维表现。在MATLAB中提供了函数slice和函数ontourslice来绘制三维物体切片图,用于三元函数R=(x,y,z)的四维可视化。slice主要是通过三维坐标点上的色变化把图形的表现能力扩展到四维。MATLAB通过对图形的线型、立面、色彩、渲染、光以及视角等属性的处理,即将重力测量的数据经计算得的重力异常数据的特性表现出来,主要