基于三维射线追踪的采空区三维正演数值模拟分析.docxVIP

基于三维射线追踪的采空区三维正演数值模拟分析 在采矿和开采过程中，地下煤矿和采石场的开采留下了许多空洞和空洞，形成了一定范围内的采空区，破坏了现有的平衡地下系统。岩层在覆岩的应力作用下发生变形、断裂、位移和冒落,并在这一过程中,重新达到应力平衡,在采空区上方形成“三带”:冒落带、裂隙带及弯曲带。未采空煤层与围岩的速度和密度差异较大,是一个明显的波阻抗界面,能够形成较强的反射波。但在其被采空后,本层、相邻区和上部区域的波阻抗发育复杂变化,在地震时间剖面上表现为反射波通向周的弯曲、减弱、直到消失以及低频响应增强,可以通过地震勘探的手段加以识别和圈定。 目前煤田三维地震勘探中,受采集成本所限和强调施工效率,通常采用常规束状正交型观测系统,其评价依据侧重于面元大小、覆盖次数、最大偏移距、偏移孔径、常规炮检距及方位角分布等参数与属性的分析,在地层倾角较小、构造发育简单、规则采空区分布等地震地质条件简单地区,这种观测系统较为适用。 实际生产中,由于小煤窑私挖乱采、越界开采等因素造成了许多历史遗留问题,导致地下采空区分布毫无规律,冒落、位移情况多变。在这种复杂情况下,使用常规束状正交型观测系统已难以满足业主的地质任务和要求,其适用性存在疑问。在当前基于目标探测的高精度构造勘探、岩性解释研究的煤田三维勘探技术发展趋势下,对观测系统设计也提出了更高的要求,需要从采集脚印压制、面元属性、照明属性以及成像角度作更全面细致的分析和评价。笔者采用三维地震数据模拟的方法,基于多种属性分析,讨论了不同观测系统类型及参数在采空区探测上的适用性。 1 实际地质资料及速度模型 三维地质模型主要依据山西省沁水盆地煤田中南部典型煤矿实际地质资料建立,便于研究,将煤层简化为单一层,采空区为范围为350 m×1 200 m的不规则体,模型参数见表1,速度模型见图1、图2。 2 比较观测系统 2.1 线5参数对比 为了科学对比不同观测系统类型及参数对采空区探测效果的影响,在覆盖次数(20次)、面元尺寸(5×10 m)、接收道数(640道)、道距(10 m)、炮点距(20 m)、炮线距(80 m)等参数相同,最大偏移距变化不大(452～576 m)的基础上,分别建立了代表常规束状观测系统的8线8炮制与宽方位直线型代表性的8线5炮制观测系统,对比其不同激发方式情况下(正交式、砖墙式、45°斜交式)各类面元属性及基于目标地质体正演照明效果(见图3、表2)。 2.2 目标层深度、偏移距与煤质地层深度的关系 面元属性包括覆盖次数、方位角和偏移距等信息,由于目的煤层埋深在410～510m,因此在限制偏移距大小约等于目的层深度的基础上,在三维观测系统上圈定了一个子集,在该范围内对以上属性加以对比分析。 2.2.1 观测系统方位角随偏移距的变化 常规束状观测系统属于窄方位观测系统,直线型观测系统属于宽方位观测系统。后者与前者相比,横向(crossline)滚动距离变小,减少了采集脚印现象,压制了横向散射噪声,可明显改善静校正效果,增加了子集数据记录道的实际有效覆盖次数。 如图4、5所示,使用直线型观测系统方位角随偏移距分布均匀,在限制偏移距的情况下,仍旧可以保证在各个角度都有较多的检波点接收信息,对于复杂采空区等不规则目标地质体来说,成像效果更好。另外,较宽的方位分布还有利于垂直测线方位信息的获得,在地下采空区延伸方向未知的情况下,可保留更完整的异常信息。 2.2.2 面元的覆盖次数 限制偏移距覆盖次数代表了面元上炮检距的分布情况,如果面元上缺失一些炮检距值,则其限制偏移距覆盖次数就会下降。在一般情况下大部分面元的限制偏移距覆盖次数较常规覆盖次数低,反映了面元上真正对叠加起作用的覆盖次数,代表了面元的属性。显然,限制偏移距覆盖次数越高,对叠加真正的贡献越大。 图6为六种不同观测系统的限制偏移距覆盖次数对比。纵向对比显示,8线5炮宽方位直线型普遍比8线8炮窄方位束状型覆盖次数高;横向对比显示,45°斜交式观测系统覆盖次数分布最均匀,但高覆盖次数面元较少,正交式与砖墙式覆盖次数分布差异不大,但高覆盖次数面元较多。因此,对于信噪比高、高差变化不大、地表一致性条件较好的复杂采空区,选用8线5炮45°斜交式观测系统效果最佳。 2.2.3 偏距对比分析 偏移距分布变化主要反映子集内的偏移距分布的均匀性。面元间偏移距分布变化差异越小的观测系统越符合水平叠加理论的假设条件,越有利于叠加速度求取和随机噪声压制。 观察图7中六种不同观测系统的偏移距分布变化,纵向对比显示,整体而言,8线5炮宽方位直线型比8线8炮窄方位束状型的偏移距分布变化小,但是45°斜交式不如8线8炮型偏移距分布变化均匀;横线对比显示,45°斜交式效果最好,砖墙式其次,正交式最差。因此,从资料处理角度考虑,选用8线8炮型45°斜交式观测