[发明专利]一种基于地震波的地质预报信息处理方法及其系统

权利要求书

1.一种基于地震波的地质预报信息处理方法，其特征在于，包括以下步骤：

S1、通过地震数据采集仪器获取M个含有N个分量数据的地震数据，每个分量数据含有L个采样点数；

S2、依次合并M个地震数据中相同位置的分量数据，得到N个数据文件，每个数据文件包含M个分量数据；

S3、分别将N个数据文件中的M个分量数据依次进行零漂校正、低通滤波和非对称放大，得到N个非对称放大后的地震数据文件；

S4、对N个非对称放大后的地震数据依次进行道间均衡，得到N个道间均衡后的地震数据文件；

S5、分别对N个道间均衡后的地震数据文件中的M个分量数据进行初至拾取，得到N个初至拾取后的地震数据文件；

S6、分别对N个初至拾取后的地震数据文件中的M个分量数据进行分量融合，得到N个1分量数据的地震数据文件；

S7、根据N个1分量数据的地震数据文件，通过采样点数规则得出顶部满足条件点的集合TS和底部满足条件点的集合BS；

根据每个采样点数在所述顶部满足条件点的集合TS和底部满足条件点的集合BS出现的次数，再对所述L个采样点数进行归类，通过判断规则对归类后的L个采样点数进行判断，得出地质情况的地质预报信息。

2.根据权利要求1所述的基于地震波的地质预报信息处理方法，其特征在于，通过所述地质情况的预报信息，通过地质情况的距离公式计算L个采样点数到掌子面的距离d，得出所述地质情况的预报信息距离掌子面的地质预报信息；

其中，距离公式为d＝2×p×V/H；

p为L个采样点数中的一个采样点数，H为地震数据的采样频率、V为地震数据的预设速度。

3.根据权利要求1所述的基于地震波的地质预报信息处理方法，其特征在于，所述步骤S1中，地震数据包括波形数据、采样频率H、数据类型和预设速度V，且M个地震数据中的N个分量数据、采样频率H、采样点数和数据类型相同。

4.根据权利要求4所述的基于地震波的地质预报信息处理方法，其特征在于，所述步骤S3中，所述步骤S3中，所述零漂校正的具体步骤为：

在N个数据文件中的M个分量数据中，设第i个文件数据文件中的第j个分量波形的数据为S_i，j＝Array[i，j]；

S_i，j包含L个采样点数{S_i，j[1]，S_i，j[2]，...，S_i，j[L-1]，S_i，j[L]}；

按照S_i，j第int(L*0.2)个点至第L个采样点数的平均值进行矫正，第j个分量数据中的L个采样点数均减去即得出零漂校正后的第i个数据文件中的第j个分量数据；

所述第i个数据文件中的第j个分量数据S_i，j的低通滤波采用巴特沃斯滤波器。

5.根据权利要求5所述的基于地震波的地质预报信息处理方法，其特征在于，所述步骤S3中，所述第i个数据文件中的第j个分量数据S_i，j的非对称放大的具体步骤为：

首先对第i个数据文件中的第j个分量数据S_i，j的L个采样点数取绝对值得到S′_i，j，设S′_i，j从第1点至第int(L/5)采样点数的平均值与第int(L*0.8)点至L点平均值的比值为R，在区间[1，R]取等间隔的L个数R_i，j，将S_i，j与R_i，j对应位置的数字相乘，得到非对称放大后的分量数据，即S_i，j[p]＝S_i，j[p]*R_i，j[p]。