煤巷三维槽波全波场数值模拟与波场特征分析（英文）

当前煤巷地震槽波超前探测数值模拟主要以二维波场研究为主,难以精确模拟三维全波场特征和全空间观测方式接收的地震记录。本文应用三维弹性波动方程的交错网格有限差分算法模拟了在巷道迎头前方煤层中以纵波震源激发的三维全波场。采用取平均值的办法实现了交错网格地震波场模拟的振动速度三分量同点记录。查明了三维顶底板对称煤巷模型和非对称模型中传播的波型及其传播特征。研究发现巷道迎头前方煤层内Rayleigh型槽波能量较强,Love型槽波能量较弱,煤层内巷道迎头反射Rayleigh型槽波在煤层顶、底板中泄露能量较小,在煤层中传播距离较远;在巷道周围存在巷道面波和折射头波;在煤巷两侧煤层内接收的地震记录中Rayleigh型槽波能量较强,受头波和面波的干扰较弱,而Love型槽波能量较弱,受到相对强能量头波和面波干扰,难以识别;增加接收点在煤层中的深度可以有效降低面波对Rayleigh型槽波的干扰,但不能有效减弱面波对Love型槽波的干扰。基于以上波场认识,发现传统的采煤工作面Love型槽波构造探测方法难以适用于掘进煤巷构造超前探测,Rayleigh型槽波可以用于掘进煤巷构造超前探测。     

煤矿井下槽波三维数值模拟及频散分析

采用交错网格高阶有限差分法编制了地震波场三维正演模拟软件,设计了基于镜像法原理处理煤矿井下近水平和起伏巷道特殊空间的算法;模拟了煤矿井下含巷道和不含巷道情况下煤层中传播的地震波场,并分析其频散特征.结果发现:由于巷道的影响,巷道壁上产生很强的巷道振型槽波,煤层中则出现了以Love型为主的槽波,据此分析了实际槽波记录的形成机理,研究结果对今后煤矿井下巷道地震超前探测和工作面弹性波透视等具有重要的理论意义和实际价值.     

平行巷道条带异常体在透射槽波层析成像中造成的假象

近年来透射法槽波地震勘探技术在煤矿勘探领域取得广泛应用,该技术依据槽波频散特征,对特定频率下槽波走时进行层析速度反演成像,在煤层厚度探测方面效果明显.由于煤矿巷道施工的特殊性,震源和检波器布置在巷帮煤层中,观测系统多采用两边式或三边式,导致槽波传播角度有限,容易对层析成像造成误差.在层析正演中采用最短路径与射线弯曲法联合,兼顾走时计算的精度和效率,在层析反演中进行正则化约束,利用平滑和阻尼因子提高算法的精度.通过对大量模型进行正反演发现,当煤层中存在垂直巷道异常体或局部异常体时,透射槽波层析反演不存在假象;当存在平行巷道的条带异常体时,层析结果出现交叉状速度异常假象.这是由于透射槽波采集得到的走时不受异常体位置影响,高值区在平面上表现为交叉型,导致层析反演出现假象.实际勘探施工条件允许情况下,可在煤层工作面四周布置震源和检波器消除这类假象.     

井下槽波地震勘探——预防煤矿灾害的一种地球物理方法

当在煤层中激发地震波时,由于顶、底板围岩速度明显高于煤层中地震波的传播速度,当地震波入射角大于临界角时会发生全反射,经过多次全反射混合叠加,在煤层中形成槽波.槽波是一种围限波,它最大的特点是频散现象,槽波在煤层中传播时还具有速度低、衰减弱等特点,因此可以用来探测煤层中的采空区、断层等地质构造.本文初步阐述了槽波地震勘探技术的研究现状,槽波的形成,槽波的特点、槽波地震勘探的工作方法以及应用实例等内容,并基于其特点对槽波地震勘探技术作出了展望.     

透射法典型槽波数据波场分析

煤矿井下巷道煤层中采集的槽波地震数据中包含多种类型地震波,对这些波组认识存在争议.通过正演模拟获得的槽波波场过于理想化,与现场采集槽波数据存在一定偏差.针对煤矿实际采集的典型槽波数据,采用时频域极化滤波方法,利用槽波数据中的纵波在波传播方向能量最强、垂直方向能量最弱的特点,将目前国内煤矿常用的德国SUMMIT槽波地震仪水平双分量检波器接收到的信号,校正到平行波传播方向和垂直波传播方向,有利于判定地震波极化运动特征.经过槽波波场分离,得到三个特征明显的波组,按照时间到达先后顺序,认为分别是折射纵波、瑞利型槽波和勒夫型槽波.其中折射纵波传播速度最快,在共炮点道集上表现为双曲线特征;瑞利型槽波主频最低,水平平行分量和垂直分量都能接收到;勒夫型槽波振动方向垂直波传播方向,只在检波器水平垂直分量有显示.     

地下全空间瞬变电磁技术在煤矿巷道掘进头的连续跟踪超前探测

文中介绍了利用地下全空间瞬变电磁(UWTEM)技术在某煤矿-250排水巷道掘进过程中进行的连续跟踪超前探测试验。根据初探结果,采用不共面同轴偶极观测方式时,UWTEM在巷道掘进头前方的有效探测范围约为70m,巷道每掘进50m进行一次跟踪探测。试验表明,利用该技术可以有效地超前预测煤矿巷道掘进头前方30~50m范围的富水异常,为及时采取水害治理措施提供依据     

反射波超前成像预报系统及其应用

矿山及煤矿井巷工程建设需要超前探测技术的有力支持,目前使用的超前反射波探测技术需进一步发展.本文提出了巷道掘进条件下有限空间三维反射波数据采集系统及其现场工作技术方法,形成了一套具有自主知识产权适用于矿山巷道掘进地质条件的反射波超前成像预报技术系统,包括现场硬件系统和处理软件.并应用实际生产中,为矿山、隧道、地铁等地下工程提供了有力的技术保障.     

矿井瞬变电磁法超前探测合成孔径成像研究

在矿井瞬变电磁超前探测数据处理和解释中,为了更好地识别巷道前方含水体电性分界面,利用扩散场与波动场间的函数关系,将矿井瞬变电磁法超前探测中瞬变电磁波场转换为拟地震波场,实现波场转换.同时,为了解决波场转换后信号幅值偏弱,以增加信噪比利于解释,利用雷达的合成孔径方法对波场转换后的拟地震信号进行相关叠加处理.将矿井瞬变电磁法巷道超前探测的同一测点不同探测方向的数据作为孔径数据进行合成孔径成像,以突出弱异常而提高信噪比,结合应用实例验证了方法的效果.研究表明,矿井瞬变电磁法合成孔径成像方法能够提高电性界面信息的信噪比和界面几何分辨率,对矿井瞬变电磁法精细超前探测含水构造具有重要的理论意义和应用价值.     

透射与反射槽波联合探测小构造应用研究——以山西龙泉矿区为例

矿井勘探方法中,槽波地震是探测小构造最重要的一种方法,小构造的定量探测对于工作面的设计和开采具有重要的意义,通过透反射联合勘探并以揭露信息进行约束能够取得准确的结果.对山西龙泉矿4203工作面进行槽波透反射勘探,巷道的煤厚为6.8 m,巷道揭露的断层f29落差为7 m,f30落差为11.8 m,利用椭圆切线法和透射层析成像信息得出槽波的解释成果,从构造发育规律和三维地震解释的断层的数量、延伸长度、走向进行统计分析,合理推断槽波探测断层在工作面内外分布及延伸长度,并通过钻探得到了验证,得出在有效利用先验地质资料条件下矿井槽波地震透反射法具有广泛的应用前景.     

掘进巷道直流电法四极装置超前探测的模拟及应用研究

为研究矿井直流电法四极装置在掘进巷道超前探测中的效果,本文以矿井直流电法超前探测原理为基础,通过分析COMSOL Multiphysics点电源正演模型模拟值与理论值的误差,验证了COMSOL Multiphysics研究直流电法超前探测的可靠性与可行性;建立巷道超前探测三维地电模型,通过正演模拟,分析了巷道空腔对数据的影响,对比了二极、三极装置与四极装置的探测效果及对不同距离低阻体的响应特征;基于相似性原理构建物理模拟实验,绘制了超前探测视电阻率等值线图,研究了四极装置对迎头前方不同形状、距离异常体的实际探测效果.研究结果表明：直流电法四极装置能够对迎头前方低阻体产生异常响应,探测数据剔除巷道空腔影响后,能够精准识别出异常体的横向位置;并且相较于板状异常体,四极装置对球状异常体的超前探测效果更好．通过将直流电法四极装置超前探测应用于工程实例中,经钻探验证与探测效果一致,说明该技术具有很强的实用性和推广性．     

STUDIES OF ROADWAY MODES IN A COAL SEAM BY DISPERSION AND POLARIZATION ANALYSIS: A CASE HISTORY*

Field data from two-component in-seam seismic measurements are used to study roadway modes and their interaction with reflected seam waves. Using the multiple-filter technique to investigate the dispersion behaviour of the different waves, it can be shown that the roadway modes disperse very similarly to the related transmitted seam waves. However, because of the free surface of the coal face, the dispersion curves of the roadway modes show a velocity reduction and a slight shift to lower frequencies compared to those of the related transmitted seam waves. Polarization analysis using hodograms, rectilinearity and polarization angle confirms these results. The parameters found by polarization analysis can be used to design polarization filters which help to separate roadway modes and reflected events in the traveltime range of superposition in the presented field case.     

煤巷小构造Rayleigh型槽波超前探测数值模拟

对煤巷小构造地震波场进行了数值模拟研究,分析了层状煤层中地震波的传播特征.研究表明:(1)在煤巷迎头前方煤层内以纵波震源激发的Rayleigh型槽波相对于体波能量较强,波列较长,波速较低.(2)沿煤层传播的Rayleigh型槽波在小构造面上产生Rayleigh型槽波反射波,反射Rayleigh型槽波垂直分量相对于水平分量能量较强.沿煤层反向传播的反射Rayleigh型槽波在煤巷迎头面上转换为沿煤巷底板传播的Rayleigh面波.沿煤巷底板可以接收到能量较强的反射Rayleigh型槽波产生的Rayleigh面波,其可以作为超前探测小构造面的特征波.在地震记录上反射Rayleigh型槽波产生的Rayleigh面波波至最迟,在时间域与其他波列时间间隔较大,其垂直分量能量相对于水平分量较强,在地震记录上容易识别.(3)在相同的地质条件下应用反射地震超前探测方法,标志煤巷迎头前方存在小构造面的反射地震波能量较弱,受煤巷顶、底板界面和采煤迎头面的强反射波干扰,在地震记录中难以识别.     

矿井全空间孔中瞬变电磁响应积分方程法数值模拟研究

基于地面-钻孔瞬变电磁法,将发射源布设于煤矿巷道中,而将接收探头放置在钻孔内,从上至下依次观测,形成矿井巷道-钻孔瞬变电磁法.为了研究巷道全空间条件下巷道-钻孔瞬变电磁场的响应特征,选取煤层底板受水害威胁的代表性地层为研究对象,建立煤层底板存在含水低阻地质异常体的三维数学模型,采用积分方程法进行数值模拟,结果表明：相对高阻的煤层侧帮对孔内垂向感应电动势的影响较小,而对孔内水平感应电动势除浅部约30 m内影响相对较大外,对深部影响较小.当煤层底板含水平低阻板状异常体时,在对异常体的纵向分辨率和异常响应的延续时间方面,孔内感应电动势的水平分量相对优于垂向分量,但垂向感应电动势的幅值强于水平感应电动势.因此,实际观测时,不仅要观测孔内感应电动势的垂向分量,也要观测水平分量.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

煤层陷落柱散射波数值模拟与成像

煤层陷落柱是煤田勘探开发中常见的一种典型的非均匀地质体.由于来自陷落柱的反射信号少、反射能量弱,使得基于反射波原理的常规地震成像方法难以有效识别陷落柱.本文以散射波理论为基础,采用数值模拟方法,研究了陷落柱的散射波场特征,研究表明地面接收的波场中含有来自陷落柱陡倾角界面的散射波场.通过共散射点道集波场的模拟,可以清晰地识别散射波,获得地下散射点和非均匀地质体的信息,判断散射点的位置,从而勾画出不均匀地质体的形态.采用等效偏移距假设抽取共散射点道集,在此基础上进行叠前偏移,对陷落柱成像;模拟与实际数据成像结果对比表明此方法能够合理地提取散射点的散射波场信息,对陷落柱形态及内部结构准确成像,是一种有效的煤田陷落柱成像方法.     

IN-SEAM SEISMIC LOVE WAVE SCATTERING MODELED BY THE FINITE ELEMENT METHOD*

The dynamic finite-element method allows frequency-dependent reflection and transmission coefficients to be computed for Love waves scattered by two-dimensional inhomogeneities in coal seams. Clean faults of zero hade angle show a throw-dependent cross-over frequency in reflection spectra, and throw-dependent conversion from fundamental to first higher mode energy in transmission spectra. Oblique faults show mode conversion in reflection spectra which is relatable to the fault hade angle by simple Huygens-theory models. Thin fracture zones or dykes normal to the seam show a reflection maximum when thickness of the zone is of order one quarter of the seam wave wavelength. Published field data from two known faults and a dyke are compared with the modeling results and support the belief that broad-band seam-wave data are capable of characterizing a seam discontinuity (throw, hade angle, dyke thickness) as well as locating it. Development of such procedures will require extensions to existing field practice and processing. This research was funded in part by the National Energy Research Development and Demonstration Council of Australia. The authors thank The Broken Hill Proprietary Company Limited for permission to quote from company case histories.     

Numerical simulations of full-wave fields and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling twodimensional wave fields and therefore cannot accurately simulate three-dimensional (3-D) full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity–stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity V_x, V_y, and V_z for the same node in 3-D staggered-grid finite difference models by calculating the average value of V_y, and V_z of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.     

单裂缝岩石标本的弹性波频响特征初步分析

通过分析单裂缝岩石标本中弹性波到时分布、波震相转换和波振幅衰减等规律，进一步用波场及频谱量刻划介质横向不均匀性。Ｐ波震相遇裂缝端部后出现“之”字形波峰，分叉成折射Ｐ波、绕射Ｐ波，它的特征也依赖于发射源和接收点相对于裂缝的位置，当波以大角度斜跨裂缝和发射点在约裂缝半长度范围时波到时和振幅差别更大。对其研究有助于认识弹性波场在含断层结构中的复杂图象特征、在野外勘探中提高地震资料解释、断裂构造识别能力和提高弹性波在非完整岩石实验中的声发射定位等能力     

Numerical study of seismic scattering and waveguide excitation in faulted coal seams

Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam.