专利名称：一种检测地震信号的控制方法及装置的制作方法
技术领域：
本发明涉及地震信号采集领域，尤其涉及一种检测地震信号的控制方法及装置。
背景技术：
地震勘探主要分采集、处理和解释三大部分，地震数据采集时需要拾取地震信号，地震波初至的拾取广泛应用在表层调查方法和地质构造的探测中，初至的拾取对矿物资源和石油资源的开发，都起着重要的作用。在地震观测中，地震事件的检测识别及初至到时拾取是一项重要而又繁琐的工作，当前出现了一种基于地震波形峭度和偏斜度的检测方法。对于一个有限长度的离散实数序列{xj ,其偏斜度S (skewness)和峭度K(kurtosis,又译峰度)定义为
权利要求
1.一种检测地震信号的控制方法，其特征在于，包括 获取地震波的N个样点的离散序列IxJ，当前时窗长度为M ； 计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和；将时窗向前移动一个样点，计算新时窗内样点的各值，是在上一个时窗的值基础上，减去上一个时窗第I个样点的值，并加上当前时窗最后一个样点的值，依次类推，得到所有时窗N个样点的峭度和偏斜度的数值；并通过得到的所有时窗N个样点的峭度和偏斜度的数值得到峭度曲线和偏斜度曲线； 通过得到的地震波的峭度曲线和偏斜度曲线，判断地震波到时的位置。
2.如权利要求I所述的方法，其特征在于， 所述得到所有时窗N个样点的峭度K的数值的步骤，包括 计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和，其中，将当前时窗长度为M的时窗内离散序列的算术和、平方和、立方和、4次方和分别记为Smj, j=l, 2，3，4,公式如下所示
3.如权利要求I所述的方法，其特征在于， 所述得到所有时窗N个样点的偏斜度S的数值的步骤，包括 计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和，其中，将当前时窗长度为M的时窗内离散序列的算术和、平方和、立方和分别记为smy j=l, 2，3，公式如下所示
4.如权利要求I所述的方法，其特征在于， 所述通过得到的地震波的峭度曲线和偏斜度曲线，判断地震波到时的位置的步骤，包括 通过地震波形峭度和偏斜度检测地震波到时来判断峭度曲线和偏斜度曲线的极大值，根据极大值点前曲线斜率最大的位置判断地震波到时的位置。
5.一种检测地震信号的控制装置，其特征在于，包括 获取地震波样点模块，用于获取地震波的N个样点的离散序列IxJ，当前时窗长度为M ； 计算模块，用于计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和；将时窗向前移动一个样点，计算新时窗内样点的各值，是在上一个时窗的值基础上，减去上一个时窗第I个样点的值，并加上当前时窗最后一个样点的值，依次类推，得到所有时窗N个样点的峭度和偏斜度的数值； 存储模块，用于存储数据 '及判断模块，用于通过得到的所有时窗N个样点的峭度和偏斜度的数值得到峭度曲线和偏斜度曲线，并通过得到的地震波的峭度曲线和偏斜度曲线，判断地震波到时的位置。
6.如权利要求5所述的装置，其特征在于， 所述计算模块得到所有时窗N个样点的峭度K的数值，是指 所述计算模块计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和、4次方和，其中，将当前时窗长度为M的时窗内离散序列的算术和、平方和、立方和、4次方和分别记为sm,-，j=l, 2，3，4,公式如下所示
7.如权利要求5所述的装置，其特征在于， 所述计算模块得到所有时窗N个样点的偏斜度S的数值，是指 所述计算模块计算当前时窗长度为M的时窗内样点的算术和、平方和、立方和，其中，将当前时窗长度为M的时窗内离散序列的算术和、平方和、立方和分别记为Sn^，j=l, 2，3，公式如下所示
8.如权利要求5所述的装置，其特征在于， 所述判断模块通过得到的地震波的峭度曲线和偏斜度曲线，判断地震波到时的位置，是指所述判断模块通过地震波形峭度和偏斜度检测地震波到时来判断峭度曲线和偏斜度曲线的极大值，根据极大值点前曲线斜率最大的位置判断地震波到时的位置。
全文摘要
一种检测地震信号的控制方法及系统，包括获取地震波的N个样点的离散序列；计算当前时窗内样点的算术和、平方和、立方和、4次方和；将时窗向前移动一个样点，计算新时窗内样点的各值，是在上一个时窗的值基础上，减去上一个时窗第1个样点的值，并加上当前时窗最后一个样点的值，依次类推，得到所有时窗N个样点的峭度和偏斜度的数值得到峭度曲线和偏斜度曲线，判断地震波到时的位置。本发明解决了当前传统的基于地震波形峭度和偏斜度的检测方法存在运算量过大导致效率低的问题。
文档编号G01V1/22GK102707312SQ20121020475
公开日2012年10月3日 申请日期2012年6月15日 优先权日2012年6月15日
发明者刘劲松 申请人:中国科学院地质与地球物理研究所