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

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

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

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

含小断层煤层Rayleigh型槽波波场和频散分析

煤层隐伏小断层在煤田地质勘探期间能否查明,是影响煤矿安全生产的主要地质因素.煤层中传播的地震槽波,非常适用于探测煤田异常构造.本文利用谱元法模拟计算了含断层模型Rayleigh型地震槽波的产生、传播过程,对于直达槽波、反射槽波、透射槽波进行了波场特征分析,频散特征分析,频谱特征分析.根据不同小落差断层模型（垂直断距分别为1/4、1/2、3/4、1个煤厚;断层倾角分别为90°、60°、45°）,分析了Rayleigh型槽波传播过程中断层对反射槽波、透射槽波能量、频谱的影响效应.基于以上分析,对于小断层模型中的Rayleigh型地震槽波取得了一些规律性认识,并对Rayleigh型地震槽波的勘探应用做了一些探讨.     

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

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

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

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

瑞利面波地形效应的数值模拟研究

在地震工程中,经常观测到地震波受地表起伏的影响而产生地形效应.本文基于数值模拟,研究了地震远场区域起伏地形对地表瑞利面波传播影响和对应的地形效应特征.结果表明:起伏地形的转折点或面的散射或反射,产生反方向传播的反射RR波和向空间扩散传播的转换RS体波,导致通过地形的瑞利面波能量极大的减弱,地形总体表现出强的隔震效应.地形起伏与地震动的波长比越大,隔震效应越强.凹陷地形表现出远比隆起地形更大的隔震效应.起伏地形转折引起的地表张角变化,导致地震动能量分布密度改变,是地形出现局部放大或减弱特征的主要因素.地表张角小于180°,地震动放大,张角越小放大作用越明显.地表张角大于180°,地震动减弱,张角越大减弱作用越强.当地形对地震动产生放大作用时,地震动水平分量比垂直分量放大效果更明显,而当地形产生减弱作用时,相对于水平分量,垂直分量更容易被减弱.     

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

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

凹陷地形对Rayleigh面波传播影响的研究

利用局域离散波数法,研究了地形对近地表爆炸源产生的Rayleigh面波的影响,展示了Rayleigh面波穿越地形过程中波的传播进程,研究了穿越地形前后Rayleigh面波能量、频率的变化,并讨论了该变化与地形起伏程度的关系,从中可见,Rayleigh面波穿越地形时,激发出新的体波和面波,并且越陡的地形对面波传播的阻碍性越强,且面波高频部分能量损失的越多.     

朝鲜核爆的Rayleigh波震级测量

利用1995年至2009年中国东北及邻近地区11个宽频带台站记录到的77个地震事件、3个化学爆炸和2次朝鲜核爆的区域地震资料,标定该区域台网的Rayleigh波震级.通过对8~25 s 周期的垂直分量Rayleigh波形进行分析,获取基于最大振幅的面波震级.计算82个区域事件不同周期的台基响应,经过台基校正后取最大振幅的面波震级为事件震级.2006年和2009年两次朝鲜核爆的面波震级分别为2.93±0.19和3.62±0.21.将地震和核爆事件的面波震级M_s与体波震级m_b(Lg)进行比较,发现根据该区域台网的数据利用M_s-m_b识别方法无法鉴别朝鲜地区的核爆与地震.朝鲜核爆的面波震级相对较大,使M_s-m_b识别方法失效,其原因可能是源区介质的不均匀性、由核爆炸冲击引发的深部的拉伸破坏被抑制,或者是近爆源区存在张性的构造预应力.假定核爆可能的埋藏深度范围是0.01~1.0 km,用Rayleigh波震级估计朝鲜核爆的当量,对2006年和2009年核爆当量的估值范围分别为0.42~3.17 kt和2.06~15.53 kt.     

隧道反射地震超前预报波场分离方法的数值模拟研究

本文以有限元数值模拟方法为依托,依次从波场快照、地震记录、波场分离、工程实例等四个角度论证了地震波传播规律和波场分离特点.模拟了单界面介质条件下地震波在隧道围岩中传播的全波波场.研究了反射波时距曲线的线性相关系数.比较了τ-p变换和F-K变换的滤波特点.结合工程实例,讨论了这两类线性波场分离方法在提取隧道前方反射波中的可行性和准确性,研究结果表明:当隧道前方地层倾角由直立变为缓倾时,反射纵波和横波的时距关系都能很好的视为线性.地层倾角的改变主要影响反射波走时和视速度.因此,应用τ-p变换和F-K变换进行波场分离,提取隧道前方单独的反射纵波和横波是可行和有效的.但对于缓倾地层,线性变换波场分离的结果存在较大的走时和视速度误差.易误滤除反射信号,需进行偏移处理.     

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.     

0.9 m薄煤层SH型槽波频散特征及波形模式

在0.9 m薄煤层中使用放炮方法做微震震源的条件,通过对所采集到的薄煤层槽波信号进行分析,发现薄煤层槽波在频域中存在高频和低频两个独立并且不连续的波段.其中高频区大约以2000 Hz为中心,低频域的中心频率约为490 Hz,并且高频域和低频域的能量差异不显著.通过时频分析,可以清晰地看到高、低两个频域几乎在同一时刻触发,并且其小波相关系数在这两个域中的分布规律表现出一定的相似性.通过对0.9 m薄煤层槽波频散曲线的理论分析可知,现场观测到的Airy震相的频率及速度和其理论值较为接近.震源置于煤层中心,且炸药能量对顶底板的扰动,对第二阶对称波形模式下槽波高频部分的形成起着关键作用.在这一对称高阶波形模式下槽波的波速基本上和煤层顶底板中S波的波速一致.由于这两个触发的波形模式在时间域中具有相似的特征,建议在高频域和低频域同时发育较好的薄煤层槽波勘探中,可以采用在同一时间域中高低频相结合的方法提高利用槽波勘探分析的效果.     

A comparative analysis of seismic response of shallow buried underground structure under incident P,SV and Rayleigh waves

In this study, A time-domain seismic response analysis method and a calculation model of the underground structure that can realize the input of seismic P, SV and Rayleigh waves are established, based on the viscoelastic artificial boundary elements and the boundary substructure method for seismic wave input. After verifying the calculation accuracy, a comparative study on seismic response of a shallow-buried, double-deck, double-span subway station structure under incident P, SV and Rayleigh waves is conducted. The research results show that there are certain differences in the cross-sectional internal force distribution characteristics of underground structures under different types of seismic waves. The research results show that there are certain differences in the internal force distribution characteristics of underground structures under different types of seismic waves. At the bottom of the side wall, the top and bottom of the center pillar of the underground structure, the section bending moments of the underground structure under the incidences of SV wave and Rayleigh wave are relatively close, and are significantly larger than the calculation result under the incidence of P wave. At the center of the side wall and the top floor of the structure, the peak value of the cross-sectional internal force under the incident Rayleigh wave is larger than the calculation result under SV wave. In addition, the floor of the underground structure under Rayleigh waves vibrates in both the horizontal and vertical directions, and the magnification effect in the vertical direction is more significant. Considering that the current seismic research of underground structures mainly considers the effect of body waves such as the shear waves, sufficient attention should be paid to the incidence of Rayleigh waves in the future seismic design of shallow underground structures.     

槽波层析成像方法在煤田勘探中的应用——以河南义马矿区为例

在地下煤田的开发中,工作面内的小构造、异常体、煤层厚度变化等是需要解决的关键问题,而槽波探测则为这些问题的解决提供了重要的物探方法.本文对河南义马矿区11061工作面进行槽波透射法测量,在巷道显示的煤层厚度变化为1.5~8 m,从理论频散曲线分析速度与厚度关系,确定125 Hz频率槽波主要用于观测厚度约为2~5 m的煤层厚度变化;有效提取了684个频散曲线,并分别拾取了125 Hz时槽波群速度与走时,采用走时层析成像方法获得工作面内煤层速度、厚度以及高应力区分布特征,回采验证了结果的正确性.     

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.     

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

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

PARAMETER ESTIMATION AND FAULT DETECTION BY THREE-COMPONENT SEISMIC AND GEOELECTRICAL SURVEYS IN A COAL MINE*

Three-component seismic and geoelectrical in-mine surveys were carried out in Lyukobanya colliery near Miskolc, Hungary to determine the in situ petrophysical parameter distributions and to detect inhomogeneities in the coal seam. The seismic measurements comprise an underground vertical seismic profile, using body waves, and an in-seam seismic amplitude-depth distribution and transmission survey, using channel waves. The geoelectrical measurements are based on the drift- and seam-sounding method. Interval traveltime-, amplitude-, multiple-filter- and polarization analysis methods are applied to the seismic data. They lead to a five-layer model for the strata including the coal seam. The coal seam and two underlying beds act as a seismic waveguide. The layer sequence supports the propagation of both normal and leaky mode channel waves of the Love- and Rayleigh type. A calculation of the total reflected energy for each interface using Knott's energy coefficients shows that the velocity ranges of high reflection energy and of normal and leaky mode wavegroups coincide. The excitation of wavegroups strongly depends on the seismic source. A simultaneous inversion of a geoelectrical drift- and seam-sounding survey prevents misinterpretations of the seismic data by clearly identifying the low-velocity coal seam as a high-resistivity bed. Calculations of dispersion and sounding curves improve the resolution of the slowness and resistivity in each layer. Both diminished amplitudes and distortions in the polarization of transmission seismo-grams and decreasing resistivities in a geoelectrical pseudosection of the coal seam are related to an inhomogeneity. A calculation of synthetic seismograms for Love and Rayleigh channel waves with the finite-difference and the Alekseev-Mikhailenko method agrees well with the field data for the main features, i.e., particular arrivals in the wave train, wavegroups, velocities and symmetries or asymmetries. This in-mine experiment demonstrates that the simultaneous acquisition, processing and interpretation of seismic and geoelectrical data improve the lithological interpretation of petrophysical parameter distributions. Coal seam inhomogeneities can also be detected more reliably by the two independent surveys than by one alone.     

高频面波方法的若干新进展

面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中,由于该方法具有非侵入性、无损、高效及价格低的特点,越来越受到浅地表地球物理和地质工程学界的重视,视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质,尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量,从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少,这使得勒夫波的频散曲线比瑞雷波的简单.因此,勒夫波反演更稳定,非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰,并具有更高的分辨率,从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言,瑞雷波的探测深度是勒夫波的1.3~1.4倍;而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体,利用反卷积消除了地震子波对勒夫波波形的影响,通过更新每个块体的S波速度来拟合勒夫波波形,从而获得地下S波速度模型.该方法不基于水平层状模型假设,适用于任意二维介质模型.