基于MATLAB实现地震速度切片成图.pdfVIP

2008年 第 3 基于MATL~J3实现地震速度切片成图 周艳萍 刘鸿福 (太原理工大学 030024) 【摘要】野外采集的天然地震数据一般都庞大而复杂，我们无法从这些地震数据中来判断出隐 含的地质现象和规律，所以有必要将这些经过预处理的天然地震数据运用科学可视化技术直观 地显示出来。本文针对天然地震数据的特点，运用MATLAB强大数学计算和图像可视化功能绘制 出新疆地区地下 300km内八层的水平速度切片图，从而可以直观地了解地下各个区域的速度变 化情况，为地质研究人员研究新疆地质特征提供参考。 【关键词】科学计算可视化；MATLAB；地震数据：水平速度切片图 1引言 科学计算可视化n¨(VisualizationinScientificComputing)是指运用计算机图形 学和图像处理技术，将科学计算过程中的计算结果数据转换为图形及图像在屏幕上显示出来 并进行交互处理的理论、方法和技术。它的运用领域也不断的扩大，目前主要于医学图像、 地质学、气象学、空间探索、建筑设计等流域 。 本文是利用MATLAB强大的数据处理和图像显示功能，通过整理和分析新疆富蕴一库尔勒的 天然地震数据，并根据这些地震数据来绘制地下各层水平速度切片，依此对了解新疆地下结 构以及天山山脉的成因、板块运动、构造单元的演化、断裂特征、地震的成因及分布规律等 地质特征具有重要的意义，为地质解释人员提供参考。 2gtATLAB平台的介绍与资料的收集 ~,TLAB是MathWorks公司于 1984推 出的一款数学软件，是用于工程科学计算的高效率 的高级语言。经过 MathWorks公司的不断完善，时至今 日，MATLAB不仅在数值科学计算上 独 占鳌头，而且图形可视化和图像处理上也毫不逊色。MATLAB可以通过简单的函数调用便 可显示二维、三维乃至四维的图形，而且通过对图形的线型、立面、色彩、渲染、光线、视 角等的控制，可把数据表现的淋漓尽致 。本文正是利用它强大数据处理和可视化功能优点， 绘制了新疆地区的水平速度切片和空间三维速度切片。 本文所用的天然地震数据是摆放在横跨天山南北，基本上沿富蕴一库尔勒一线摆开，中 部地段分布较密的48个台站接收的P波到时数据。台站最北位于富蕴附近(北纬 47。8．28)， 最南位于库尔勒 以南 (北纬 4O。57．44)，南北范围长度接近 700公里，东西宽约 400公 里 。并通过对大量的地震事件 的筛选，挑取震源 4．5级以上并且信号 良好的地震 P波到时 数据，并利用 ACH算法所获得不同深度的八层地震数据。其数据最终形式如表 1，数据格 式是．dat，有三列变量组成：第一列是经度坐标，第二列是纬度坐标，第三列是速度特性值。 21 工程她质计算机应用 2008年 第3期 总51期 表 I 新疆天然地震数据 (部分) 经度坐标 维度坐标 速度值 87．756783 44．505299． 5．899659 87．815041 44．489922 5．914009 87．873291 44．474541 5．927305 87．931541 44．459164 5．938489 87．989792 44．443783 5．947208 88．048042 44．428402 5．951290 88．106293 44．413025