[发明专利]一种煤矿工作面透射地震体波与槽波分离方法

权利要求书

1.一种煤矿工作面透射地震体波与槽波分离方法，其特征在于，该方法包括如下步骤：

步骤1，读取煤矿工作面透射地震数据，按透射波观测系统设置道间距数组；

步骤2，将采样时间平方，对平方后的数组以最小间隔线性插值，并对地震数据线性插值，得到抛物线形式数据体；

步骤3，在频率域内对有效频段采用快速迭代收缩阈值法求解抛物线形式Radon变换；

步骤4，在时间-慢度域内对透射体波与槽波进行分离，对分离数据体应用Radon反变换并将采样时间开平方，即可得到分离后的体波和槽波。

2.根据权利要求1所述的一种煤矿工作面透射地震体波与槽波分离方法，其特征在于，所述步骤1中透射波观测系统道间距数组设置如下：

令道间距为△x，震源在另一巷道的投影点为零点，零点右边依次为为△x,2△x,3△x,...,零点左边依次为为-△x,-2△x,-3△x,...。

3.根据权利要求1所述的一种煤矿工作面透射地震体波与槽波分离方法，其特征在于，所述步骤3中具体求解过程如下：

步骤3.1，对插值后的地震数据体进行频谱分析，分析该数据体的有效频段；

步骤3.2，在频率域内，对有效频段的每一频率均采用快速迭代收缩阈值法求解稀疏反演问题：

式中，D为频率域内地震数据，M为频率域内Radon变换数据体，L为频率域Radon变换算子，λ是正则化因子；

令||M||₁＝φ₂，则

式中，k为迭代次数；

步骤3.3，判断有效频段是否处理完毕，如未完，则返回步骤3.2，如已计算完，则输出频率域内Radon变换数据体M，进行傅里叶反变换，得到时间-慢度域数据体。

4.根据权利要求1所述的一种煤矿工作面透射地震体波与槽波分离方法，其特征在于，所述步骤4中具体过程如下：

步骤4.1，由于透射体波传播速度大于透射槽波，在步骤3求得的透射波时间-慢度域分布区图中，分布在时间小，慢度小区域的为体波，时间大，慢度大区域的为槽波；

步骤4.2，对分离出的时间-慢度域数据进行傅里叶变换，然后利用频率域抛物线Radon反变换得到抛物线形式地震数据体，对数据体进行傅里叶反变换并对采样时间开平方，抽取与原数据体相对应时间的地震数据即得到分离后的透射体波和槽波数据体。