[发明专利]一种适用于煤层底板复杂陷落柱的全波形反演方法

说明书

本发明公开了一种适用于煤层底板复杂陷落柱的全波形反演方法，先获取震源的观测数据，然后经过计算获得震源的正传波场和伴随源的反传波场，然后确定速度场模型的梯度表达式，接着根据正传波场和反传波场及梯度计算表达式，计算得到速度场模型参数的梯度，有利于获得准确的速度场模型参数的更新方向，最终通过对速度场模型参数进行迭代更新，直至数据残差达到设定阈值或迭代更新次数达到预设值，停止迭代更新确定速度场模型参数，此时采用该参数的速度场模型进行全波形反演。因此本发明能对煤层底板下方条件复杂的陷落柱进行精确定位及形态探测，从而为后续开采提供数据支撑。

技术领域

本发明涉及一种适用于煤层底板复杂陷落柱的全波形反演方法，属于地震勘探技术领域。

背景技术

陷落柱是煤系下伏灰岩中溶洞顶部覆岩发生塌陷所形成的陷落体,是我国华北型煤田中一种典型的隐伏垂向构造。因具有随机性及不确定性，其诱发的突水隐蔽性强、速度快、水量大、破坏性强，对煤矿安全生产的危害极大。地震三维勘探在绿色煤炭资源丰富的中西部地区受地表地形以及浅部上组煤的开采和多层采空区及地表沉陷的影响，陷落柱的探测精度受到了很大限制，矿井地震勘探具有距离目标体近、地震波能量和高频成分衰减少、地震波传播路径干扰少、分辨率高等特点。因此，开展井下针对陷落柱探测具有重要意义。

目前工作面回采过程中针对煤层陷落柱的探测以槽波最为典型。当陷落柱发育至煤层时，会破坏煤层的连续性，利用槽波的地震波属性进行探测，典型的探测方式是利用工作面其中一条巷道布置检波器，另外一条巷道进行震源激发进行CT探测。自Evison在新西兰首次于煤层中激发并且记录到了槽波，随后众多学者开展了大量研究，目前，槽波勘探已成为煤层工作面回采面内构造探测的常规工作。

在槽波发育的地区，槽波勘探作为主要的煤层勘探方法，但是，若煤层相对较薄或者煤层槽波不发育时，槽波则不能作为探测手段，此时发现沿着顶底板传播的折射波具有陷落柱探测可能性。因此，当单一类型的地震波成像精度受限，开展全波场信息的综合使用可以提升成像分辨率。为了获得煤田高精度的速度模型，全波形反演方法在煤田中的应用已经有了一些初探。有学者将全波形反演方法应用到煤田多参数反演中，并且全波形反演结果与测井曲线的吻合度比较高。还有一些学者在煤矿井下针对煤层内地质构造高精度成像问题，发展了槽波全波形反演方法，相比于常规的槽波成像方法，利用全波形反演可以获得更高精度的反演结果。但是目前研究仅限于合成数据的参数反演，同时现场陷落柱条件复杂，目前的全波形反演方法无法对条件复杂多个陷落柱的精确反演，因此如何能对条件复杂的陷落柱进行精确的全波形反演是亟需解决的问题。

发明内容

针对上述现有技术存在的问题，本发明提供一种适用于煤层底板复杂陷落柱的全波形反演方法，能对煤层底板下方条件复杂的陷落柱进行精确定位及形态探测，从而为后续开采提供数据支撑。

为了实现上述目的，本发明采用的技术方案是：一种适用于煤层底板复杂陷落柱的全波形反演方法，该方法的具体步骤为：

A、在煤层工作面一侧的巷道壁相互等间距布设多个震源；在煤层工作面另一侧的巷道壁相互等间距布设多个检波器，多个检波器均通过连接主线与采集基站及采集主机连接，使检波器、采集基站及采集主机形成地震观测系统；激发多个震源，通过地震观测系统接收观测数据并获得各个地震子波；

B、先构建基于卷积型的目标函数，具体定义如下：

式中，d是观测数据，u是合成数据，x_r是检波器的位置参数，*是时间卷积算子，x_ref表示提取参考道的位置参数；表示欧式范数，m＝(v_p，v_s)为速度场模型参数，密度为常数；

其中合成数据中的正传波场通过弹性波动方程合成，具体公式如下：

L[m(x)]u(x,t；m)＝s(x_s,t) (2)