基于matlab的地震活动性分析.doc

Matlab在地震活动性图像分析中的应用 1），李红光2） 1）河北省地震局 2）中国地震应急搜救中心 摘要：地震活动性分析是地震预测、地震工程的一个重要依据，地震活动性分析又多是通过图像来表现。Matlab是一种简单易学、强大的计算功能和编程可视化的计算机语言。本文用Matlab语言编程，实现了地震统计区内地震的快速选取，并根据这些选中的地震进行地震活动性分析。 关键词：Matlab语言；地震活动性 引言 地震活动性研究就是通过分析一定震级区间内的地震时间、空间的分布特征，探讨其物理含义，进而对地震发生的规律进行科学总结。通过地震活动性研究，可对地壳介质非均匀性和运动形态有宏观的了解和总体把握，因此可服务于地震预测和地球动力过程等研究[1]。在地震安全性评价中，通过地震活动性分析，为工程场地一定时间内的地震活动性趋势和地震环境做出评价，为合理划分潜源区和确定其地震活动性参数提供依据[2]。地震活动图像的分析方法很多，有简单的图像描述法，如地震震中分布、蠕变曲线、M-T图等；也有采用统计参数表征地震活动时空图像特征的方法，如b值、地震活动度S、地震能流密度、地震强度因子MF分布等。Matlab具有强大的计算能力、计算结果可视化和编程效率高的优势，它是地震活动性分析的一个有力工具。 Matlab是1984年由美国MathWorks公司推出的荣誉产品。早在20世纪80年代中期，Matlab就在我国出现，大规模流行时再90年代中期以后。现在Matlab已被广泛应用在科学研究、工程计算等方面。Matlab采用全新的数据类型和面向对象编程技术，采用了新控制流和函数结构，特别是包含很多常用的子函数，非计算机专业人员非常容易用Matlab来实现很复杂的计算程序。并且Maltab提供了图像处理功能，可以很方便的生成图形。在地震活动性数字图像分析中，用Matlab可以很简单、方便的实现研究人员的思想。 Matlab在地震预测、地震工程、地震波记录处理和地震定位等方面得到了应用。王凤基于Matlab的BP预测模型进行地震前兆预测的研究[3]；焦旭霞利用Matlab分析地震日常记录[4]；李晋平利用Matlab实现了地震数据的小波变换[5]；李敬利用Malab中的数字信号处理工具箱进行数字地震记录中干扰波的去除[6]。Matlab以它的简单易学、强大的工具箱和结果可视化等优势已经在地震行业被广泛应用。 本文通过编写Matlab程序，快速画出任意地震统计区M-T图和应变释放曲线图来说明Matlab在地震活动性分析中的优势。 一 计算地震与地震统计区关系 画M-T图和应变释放曲线图首先要选择地震统计区内的所有符合要求的地震目录。地震统计区是CPSHA方法中特有的概念，区内大小地震的发生具有统计上的协调性，满足G-R关系。而且地震统计区的边界是决定地震样本的不同归属，具有很重要的意义[7]。所以在地震安全性评价中为了准确的统计处地震统计区地震活动性参数，我们必须根据地震统计区的边界精确的判定出地震的归属。 我们可以把地震统计区看作一个任意的多边形，地震震中作为一个点，判断地震是否发生在地震统计区内，可以转换为判断一个点是否在多边形内的问题。用计算机实现这个问题，不仅大大的加快了工作效率，还提高了结果的精度。 二 判断点与多边形的关系理论 点与多边形位置关系的判断问题是计算机图形学中的基本问题，也是多年来国内外的诸多学者比较关注的问题，他们针对点在多边形内的检测问题作了广泛的研究，提出了许多算法以提高计算效率[8]。归纳起来常用的方法有3种：叉积判断法、夹角之和检验法与射线法．其中，叉积判断法要同时对全部边进行叉积运算，角度法则要使用复杂的三角运算，因此二者的计算量都非常大．目前，工程上应用最多的还是射线法，这种方法简单可靠。本文采用射线法来解决这个问题，下面简单的介绍射线法的原理。 地震统计区是一组首尾相连的线段的有限集合，线段的端点就是地震统计区的顶点，所以用有序的顶点序列P１，P２，? Pｎ，Ｐｎ＋１表示地震统计区，其中Ｐｎ＋１＝Ｐ１.每两个相邻的点连接起来就是地震统计区的一个边。 射线法判断点与多边形关系的关键是通过判断点引一条射线，求射线与多边形边的交点个数，如果交点个数为偶数则点在多边形外，如果为奇数则点在多边形内．如图1所示，有5个判断点分别对应5条水平向右的射线，它们与多边形相交的位置各不相同：如 = 1 \* GB3 ①-⑥ 只与多边形的边相交； = 2 \* GB3 ②-⑦ 和 = 3 \* GB3 ③-⑧ 与多边形的顶点相交；射线 = 4 \* GB3 ④-⑨ 和 = 5 \* GB3 ⑤-⑩ 与多边形的某些边重合．不同情况下交点个数的求解方法不同，下面分析不同情况下交点个数的求解过程． １.一般情况 如图1中的 = 1 \* GB3