基于matlab的地震活动性分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Matlab在地震活动性图像分析中的应用
1),李红光2)
1)河北省地震局
2)中国地震应急搜救中心
摘要:地震活动性分析是地震预测、地震工程的一个重要依据,地震活动性分析又多是通过图像来表现。Matlab是一种简单易学、强大的计算功能和编程可视化的计算机语言。本文用Matlab语言编程,实现了地震统计区内地震的快速选取,并根据这些选中的地震进行地震活动性分析。
关键词:Matlab语言;地震活动性
引言
地震活动性研究就是通过分析一定震级区间内的地震时间、空间的分布特征,探讨其物理含义,进而对地震发生的规律进行科学总结。通过地震活动性研究,可对地壳介质非均匀性和运动形态有宏观的了解和总体把握,因此可服务于地震预测和地球动力过程等研究[1]。在地震安全性评价中,通过地震活动性分析,为工程场地一定时间内的地震活动性趋势和地震环境做出评价,为合理划分潜源区和确定其地震活动性参数提供依据[2]。地震活动图像的分析方法很多,有简单的图像描述法,如地震震中分布、蠕变曲线、M-T图等;也有采用统计参数表征地震活动时空图像特征的方法,如b值、地震活动度S、地震能流密度、地震强度因子MF分布等。Matlab具有
强大的计算能力、计算结果可视化和编程效率高的优势,它是地震活动性分析的一个有力工具。
Matlab是1984年由美国MathWorks公司推出的荣誉产品。早在20世纪80年代中期,Matlab就在我国出现,大规模流行时再90年代中期以后。现在Matlab已被广泛应用在科学研究、工程计算等方面。M atlab采用全新的数据类型和面向对象编程技术,采用了新控制流和函数结构,特别是包含很多常用的子函数,非计算机专业人员非常容易用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种:叉积判断法、夹角之和检验法与射线法.其中,叉积判断法要同时对全部边进行叉积运算,角度法则要使用复杂的三角运算,因此二者的计算量都非常大.目前,工程上应用最多的还是射线法,这种方法简单可靠。本文采用射线法来解决这个问题,下面简单的介绍射线法的原理。
地震统计区是一组首尾相连的线段的有限集合,线段的端点就是地震统计区的顶点,所以用有序的顶点序列P1,P2,⋯ Pn,Pn+1表示地震统计区,其中Pn+1=P1.每两个相邻的点连接起来就是地震统计区的一个边。
射线法判断点与多边形关系的关键是通过判断点引一条射线,
求射线与多边形边的交点个数,如果交点个数为偶数则点在多边形外,如果为奇数则点在多边形内.如图1所示,有5个判断点分别对应5条水平向右的射线,它们与多边形相交的位置各不相同:如①-⑥只与多边形的边相交;②-⑦和③-⑧与多边形的顶点相交;射线④-⑨和⑤-⑩与多边形的某些边重合.不同情况下交点个数的求解方法不同,下面分析不同情况下交点个数的求解过程.1.一般情况
如图1中的①-⑥只与多边形的边相交,不通过多边形的顶点,属于射线与边相交的一般情况,此时可以很容易得到交点个数。
2.射线通过多边形顶点的情况。射线通过多边形顶点的情况有2种:I、通过单个或多个互不相邻的顶点;II、通过2个或多个相邻顶点,即射线与多边形的边重合。
如图1中所示的射线②-⑦和③-⑧属于第1种情况,即射线②-⑦通过单个顶点,射线③-⑧通过2个不相邻顶点。射线②-⑦通过多边形的顶点P2,该顶点对应的2条边都位于顶点的上侧;射线③-⑧通过多边形的顶点P11和P4其中P11对应的2条边都位于顶点的下侧,P4对应的2条边则分别位于顶点的两侧.由图1可以看出:射线②-⑦通过顶点P2并与边P1P12和P3P4各有一个交点,如果把P2算作一个交点,则射线与多边形的交点个数等于3,根据射线法的奇偶判断原则,点②位于多边形内。但从图1中可以看出点②位于多边形外,因此判断有误;射线③-⑧通过P4和P11顶点,点③应该在多边形外,而从图1可看出,点③在多边形内,由此看出,当射线通过多
边形的顶点时,交点个数的计算与射线通过的顶点所对应边的情况有关。当射线通过多边形的顶点时,如果该顶点的前后2个相邻顶点位于射线的同侧,则认为在相交的顶点处有2个交点或者没有交点,因为断定点是否在多边形内的主要依据是射线与多边形边交点个数的奇偶性,所以认为这种顶点不作为交点,如顶点P2。如果该顶点的前后2个相邻顶点位于射线的异侧,则认为在相交处的顶点就是一个顶点,如顶点P4。
射线通过相邻顶点的情况如图1所示的射线④-⑨和⑤-⑩,这2条射线分别通过多边形边P5P6和P8P9,此时由于线段与射线重合,所以可以分别把P5、P6和P8、P9各作为一个交点,此后的交点个数的判断与射线通过单个顶点时的判断方法相同。
3.判断点在多边形的边上的情况。如果判断点与边的下侧端点连线的斜率等于该边斜率,则判断点在多边形的边上。此外,如果判断点与多边形的某个顶点重合,判断点在多边形的边上。