MATLAB线性拟合测定古田地区震动持续时间震级.docVIP

MATLAB线性拟合测定古田地区 震动持续时间震级 段刚 （福建省地震局监测中心） 摘要：本文介绍了振动持续时间震级，并利用MATLAB线性拟合来确定公式中的未知参数，讨论振动持续时间震级的优缺点及适用性。 主题词： 振动持续时间 震级 MATLAB 线性拟合 引言 震级是表征地震强弱的测度，是地震的基本参数之一，常用的震级有近震震级，面波震级，体波震级和振动持续时间震级，另外，对于巨大地震描述还有矩震级和谱震级M T 。、、均属于里氏震级系统，它是建立在单一频率波振幅定震级的基础上，振幅大小表现了震源体内地震所释放的能量大小。、及M T 则是利用其他因素来测定震级。 1. 振动持续时间震级 在实际观测中发现这样一个重要现象：同一个地震，当各台比较稳定时，各台记录的振动持续时间也比较稳定。一般说来，地震越大，振动的持续时间越长；地震越小，振动的持续时间越短。振动持续时间在一定距离内几乎和震中距无关。既然振动的持续时间与地震的强弱有如此的依赖关系，就可以用振动的持续时间作为一个新的测度去测定地震震级，建立如下关系式： （1） 式中： F-P 是表示从P波起始，到振动衰减至相当于干扰背景水平的2倍时为止（表示为F）振动所经历的时间，以秒为单位，，是待定系数，很小，在震中距小于600km时可以忽略不计，此时（1）式可写为 （2） 式中表明与 F-P 为线性函数关系，表明F-P 越大，就越大，这个公式基本反映了上述振动持续时间和地震强弱的依赖关系。 2. 待定系数确定 2.1. 原理 在现实中经常遇到这样的问题，一个函数并不是以某个数学表达式的形式给出，而是以一些自变量与因变量的对应表给出，在研究这一系列成对的数据（x,y、x,y、…、x.y）之间的相互关系时，通常把因变量y与自变量x表达成线性回归关系式，称为经验表达式，之所以叫经验是因为我们无法知道两者的因果数学表达式，回归关系式存在一定偏差，只是偏差很小罢了。为了减小误差的影响，往往要增加测量次数，出现了多次测量，方程的个数多于未知量个数，这时就不能由方程直接算出唯一解，要让经验方程最好的描述实际关系，就要使得观测点和估计点的距离的平方和达到最小（之所以要平方是考虑到误差可正可负直接相加可能相互抵消），这就是最小二乘法的思想。利用此方法求解线性回归关系式 2 中的待定系数与。 2.2 资料选取 由于与地震仪器的放大倍数或灵敏度有关，与仪器的频带宽度有关，所以在选取资料时应选取性能稳定期间的地震，P是震动开始时间，可以是，也可以是；F是从P波起始，到振动衰减至相当于干扰背景水平的2倍时为止，且以后10s内不再出现2倍于背景振幅的尾波时的时间， F-P 取各台平均值；为对比方便，以震级为基准，即取 。 基于以上要求，选取古田水口 GTSK 、闽清雄江 MQXJ 、南平樟湖 NPZH 、古田玉田 GTYT 四个台资料（这四个台在定位古田地区地震时满足精度要求，且记录良好），将四个台振动持续时间的平均值作为一次地震的平均振动持续时间 D ,四个台震级的平均值作为此次地震的平均震级 M ，共选取10次地震，如表一： 表一：古田地区平均震级与平均振动持续时间 序号 1 2 3 4 5 6 7 8 9 10 D s 27.49 31.71 34.76 44.31 47.66 56.21 57.59 58.68 62.16 88.53 M 0.3 0.4 0.4 0.6 0.6 0.8 1.0 1.1 1.2 1.6 2.3 MATLAB数据处理 MATLAB拥有强大的数据处理能力，在多学科领域有广泛的应用，利用其数据拟合能力可以计算出参数的最佳估计值，提高了效率和精度。根据表一的数据可先绘制D-M关系图，MATLAB程序如下： D [27.49,31.71,34.76,44.31,47.66,56.21,57.59,58.68,62.16,88.53]; M [0.3,0.4,0.4,0.6,0.6,0.8,1.0,1.1,1.2,1.6]; plot D,M,-* 得出的图形如图一所示。 D-M关系大致为一直线，故可以利用线性拟合进行参数确定。 图一：D-M关系图 由（1.2）式可知，需要先将表一中的D数据进行取对数运算，为进一步处理做准备。MATLAB程序如下： D [27.49,31.71,34.76,44.31,47.66,56.21,57.59,58.68,62.16,88.53]; i [1:10]; p log10 D i 得出结果如表二： 表二：古田地区平均震级与平均振动持续时间对数值 序号 1 2 3 4 5 6 7 8 9 10 lgD 1.439 1.501 1.541 1.647 1.678 1.749 1.760 1.769 1.794