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

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

HTI煤层介质槽波波场与频散特征初步研究

煤层内裂隙较为发育,具有明显的各向异性.目前槽波理论研究以各向同性介质为主,对HTI介质中槽波及其频散性质研究很少.本文以弱各向异性、含垂直裂隙HTI煤层介质为研究对象,研究了HTI煤层介质中的三维槽波波场,采用交错网格高阶有限差分法模拟槽波,推导了三层水平层状HTI煤层介质的Love型槽波理论频散公式和振幅深度分布,分析了HTI各弹性参数对频散曲线的影响.HTI介质和各向同性介质基阶Love槽波频散曲线差异较小,高阶较大;煤厚主要影响Airy相频率,而Airy相速度不变;煤层vs对Airy相速度影响很大;煤层γ对基阶Love槽波影响很小,高阶稍大.各波偏振方向不再与波的传播方向平行或垂直,而是呈一定夹角.利用基阶Love槽波频散曲线推测裂隙发育较为困难,可利用高阶频散曲线.     

Analysis of quality factors for Rayleigh channel waves

To facilitate investigation of the effect of imperfect elastic dissipation on thepropagation of Rayleigh-type channel waves and use of their quality factors in investigationsof the properties of coal seams, a simple method for calculating the quality factor QR isproposed in this paper. Introduction of complex velocities into the dispersion function allowscalculation of the dispersion function of Rayleigh-type channel waves in coal seams. By thecontrol variable method, we analyzed changes in QR with changes in coal seam thickness andP- and S-wave Q-factors within the coal seam and adjacent rock layers. The numerical resultsshow that the trend of the QR curve is consistent with the group velocity curve. The minimumQR value occurs at the Airy phase frequency; the Airy phase frequency decreases as coal seamthickness increases. The value of QR increases with increasing Qs2 （quality factor for S wavein coal seam）. We can compensate for the absorption of Rayleigh-type channel waves usingthe computed QR curve. Inversion of the QR curve can also be used to predict the thicknessesand litholoeies of coal seams.     

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

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

CALCULATION OF DISPERSION CURVES AND AMPLITUDE-DEPTH DISTRIBUTIONS OF LOVE CHANNEL WAVES IN HORIZONTALLY-LAYERED MEDIA*

We present dispersion curves, and amplitude-depth distributions of the fundamental and first higher mode of Love seam waves for two characteristic seam models. The first model consists of four layers, representing a coal seam underlain by a root clay of variable thickness. The second model consists of five layers, representing coal seams containing a dirt band with variable position and thickness. The simple three-layer model is used for reference. It is shown that at higher frequencies, depending on the thickness of the root clay and the dirt band, the coal layers alone act as a wave guide, whereas at low frequencies all layers act together as a channel. Depending on the thickness, and position of the dirt band and the root clay, in the dispersion curves of the group velocity, secondary minima grow in addition to the absolute minima. Furthermore, the dispersion curves of the group velocity of the two modes can overlap. In all these cases, wave groups in addition to the Airy phase of the fundamental mode (propagating with minimum group velocity) occur on the seismograms recorded in in-seam seismic surveys, thus impeding their interpretation. Hence, we suggest the estimation of the dispersion characteristics of Love seam waves in coal seams under investigation preceding actual field surveys. All numerical calculations were performed using a fast and stable phase recursion algorithm.     

基于背景噪声初步研究河北及邻区的剪切波速度结构

本文根据2010年1～12月河北及邻区的83个宽频地震仪12个月连续噪声记录, 分析了河北及邻区瑞利面波的群速度频散曲线并反演了主要分区内的典型路径剪切波速度结构。 首先采用多重滤波方法提取了台站对5～50 s的面波群速度频散曲线, 然后用Herrmann线性反演方法反演了剪切波速度结构。 结果表明, 群速度频散曲线及剪切波速度分布特征与地表地质和构造特征表现出较好的相关性, 清晰地揭示了地壳内部的横向速度变化。 在短周期(8～20 s), 拥有较厚的沉积层的平原地区表现为明显的低速特征, 而隆起地区则表现为较高的群速度分布特征: 随着周期的增加(>20 s)速度的特征有所改变, 30 s之后由于受地壳厚度的影响, 两者的速度差异逐渐缩小, 在中下地壳波速度随深度增加而增大。 同时, 研究发现在研究区域内西南方向的噪声源占主导作用。     

THE INFLUENCE OF AN ASYMMETRY IN THE SEQUENCE ROCK/COAL/ROCK ON THE PROPAGATION OF RAYLEIGH SEAM WAVES*

Small offsets in hard coal seams can be detected with the aid of seam (channel) waves. Transmission and reflection of seam waves depend, among other parameters, upon the symmetry properties of the sequence rock/coal/rock. Two typical unsymmetrical sequences are found in European coal deposits: (i) coal seams with roof and floor of differing acoustic impedance and (ii) coal seams interlayered with rock and soil. Two-dimensional analog models with appropriate impedance contrasts are used to study the effect of the unsymmetrical layers upon the propagation of Rayleigh seam waves. Data analysis is based upon amplitude measurements both parallel and perpendicular to the layers and dispersion curves. The effect of unsymmetrical roof (rock 1) and floor (rock 2) was studied with models containing homogeneous coal seams. Leaky mode wave groups with phase velocities (c_R) in the range between the SV-wave velocities (β_{r1, 2}) of the two rock materials, i.e. β_r1≥c_R > β_r2, form a characteristic part of the Rayleigh seam wave signal. Using Knott's energy coefficient calculations it is shown that in that range energy leakage into the surrounding rock by refracted SV-waves is restricted to only one of the two interfaces, namely coal/rock 2. At the other interface, coal/rock 1, all waves are totally reflected. Thus, the high amplitudes of these leaky mode wave groups are explained by “quasi-normal mode” features. The influence of a dirt bed on wave propagation was studied in models where the roof and the floor have the same elastic properties. The maximum thickness of the dirt bed did not exceed 20% of the total seam thickness. The effect of the bed's location within the seam was also investigated. For all recorded normal-mode wave groups either the total seam or the coal layers could be regarded as wave guides. This was shown by the fact that the phases could be associated with the phase velocity dispersion curves calculated for the symmetrical sequence rock/coal/rock. These curves are relevant under the condition that the thickness of the coal layer assumed under the calculation coincides with the thickness of the effective wave guide of the respective wave groups. Wave groups guided in the total seam are not influenced by either the thickness or the position of the dirt bed. On the other hand, for wave groups guided in the coal layers, the quotient of signal amplitudes in the coal layers is influenced by the position of the dirt bed.     

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

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

Automatic pickup of arrival time of channel wave based on multi-channel constraints

Accurately detecting the arrival time of a channel wave in a coal seam is very important for in-seam seismic data processing. The arrival time greatly affects the accuracy of the channel wave inversion and the computed tomography (CT) result. However, because the signal-to-noise ratio of in-seam seismic data is reduced by the long wavelength and strong frequency dispersion, accurately timing the arrival of channel waves is extremely difficult. For this purpose, we propose a method that automatically picks up the arrival time of channel waves based on multi-channel constraints. We first estimate the Jaccard similarity coefficient of two ray paths, then apply it as a weight coefficient for stacking the multichannel dispersion spectra. The reasonableness and effectiveness of the proposed method is verified in an actual data application. Most importantly, the method increases the degree of automation and the pickup precision of the channel-wave arrival time.     

套管井偶极弯曲波频散向高频偏移的特性

对套管井偶极弯曲模式波的频散特性进行了系统的数值考察、实例对比和分析.发现套管井弯曲波频散曲线随地层横波速度的降低,特别是地层横波速度小于2000 m·s^-1以下,会迅速移向高频区,偶极弯曲波基础模式主频散区（或截止频率）可出现在13 kHz以上,以致超出了现行低频偶极子声波测井仪的激发与接收频带,这是一过去没有被研究者注意到的现象,并进一步被现场实例所证实.研究表明控制套管井弯曲波频散曲线主频散区位置的主要是钢套管的厚度和地层横波速度.对地层横波速度大于井孔流体声速的快速地层,在钢套管壁厚一定（8 mm）的情况下,频散曲线主频散区可移至11 kHz以上,可能出现的最大可能频域位置是同一井孔内径,井外全钢时的频散曲线上等于、小于地层横波速度那一段,这对各种地层和套管参数都是适用的.对地层横波速度小于等于井孔流体声速（1500 m·s^-1）的慢速地层,弯曲波频散曲线随地层横波速度的降低移向高频区的特点更为明显,可能移至16 kHz以上;而套管厚度的影响,也比快速地层大的多,对地层横波速度小于1380 m·s^-1的慢速地层,无论用多高的频率激发,都不能在现行使用的各类套管井（壁厚6~12 mm）中用偶极声波测井仪测到弯曲模式波.     

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

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

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.     

Shear Wave Velocity Structure of the Sinai Peninsula from Rayleigh Wave Analysis

The lithospheric structure of the Sinai Peninsula is shown by means of nine shear velocity profiles for depths ranging from zero to 50 km, determined from the Rayleigh wave analysis. The traces of 30 earthquakes, which occurred from 1992 to 1999 in and around the study area, have been used to obtain Rayleigh wave dispersion. These earthquakes were registered by a broadband station located in Egypt (KEG station). The dispersion curves were obtained for periods between 3 and 40 s, by digital filtering with a combination of MFT and TVF filtering techniques. After that, all seismic events were grouped in source zones to obtain a dispersion curve for each source-station path. These dispersion curves were inverted according to generalized inversion theory, to obtain shear wave velocity models for each source-station path, which is the main goal of this study. The shear velocity structure obtained for the Sinai Peninsula is shown through the shear velocity distributions with depth. These results agree well with the geology and other geophysical results, previously obtained from seismic and gravity data. The obtained velocity models suggest the existence of lateral and vertical heterogeneity. The shear velocity increases generally with depth for all paths analyzed in the study area. Nevertheless, in some paths a small low velocity channel in the upper or lower crust occurs. Along these profiles, it is found that the crustal structure of the Sinai Peninsula consists of three principal layers: upper crust with a sedimentary layer and lower crust. The upper crust has a sedimentary cover of 2 km thick with an average S-velocity of 2.53 km/s. This upper crust has a variable thickness ranging from 12 to 18 km, with S-wave velocity ranging from 3.24 to 3.69 km/s. The Moho discontinuity is located at a depth of 30 km, which is reflected by a sharp increase in the S-velocity values that jump from 3.70–4.12 to 4.33–4.61 km/s.     

Shear wave crustal velocity model of the Western Bohemian Massif from Love wave phase velocity dispersion

We propose a new quantitative determination of shear wave velocities for distinct geological units in the Bohemian Massif, Czech Republic (Central Europe). The phase velocities of fundamental Love wave modes are measured along two long profiles (~200 km) crossing three major geological units and one rift-like structure of the studied region. We have developed a modified version of the classical multiple filtering technique for the frequency-time analysis and we apply it to two-station phase velocity estimation. Tests of both the analysis and inversion are provided. Seismograms of three Aegean Sea earthquakes are analyzed. One of the two profiles is further divided into four shorter sub-profiles. The long profiles yield smooth dispersion curves; while the curves of the sub-profiles have complicated shapes. Dispersion curve undulations are interpreted as period-dependent apparent velocity anomalies caused both by different backazimuths of surface wave propagation and by surface wave mode coupling. An appropriate backazimuth of propagation is found for each period, and the dispersion curves are corrected for this true propagation direction. Both the curves for the long and short profiles are inverted for a 1D shear wave velocity model of the crust. Subsurface shear wave velocities are found to be around 2.9 km/s for all four studied sub-profiles. Two of the profiles crossing the older Moldanubian and Teplá-Barrandian units are characterized by higher velocities of 3.8 km/s in the upper crust while for the Saxothuringian unit we find the velocity slightly lower, around 3.6 km/s at the same depths. We obtain an indication of a shear wave low velocity zone above Moho in the Moldanubian and Teplá-Barrandian units. The area of the Eger Rift (Teplá-Barrandian–Saxothuringian unit contact) is significantly different from all other three units. Low upper crust velocities suggest sedimentary and volcanic filling of the rift as well as fluid activity causing the earthquake swarms. Higher velocities in the lower crust together with weak or even missing Moho implies the upper mantle updoming.     

Calculation of Rayleigh-wave phase velocities due to models with a high-velocity surface layer

Rayleigh-wave phase velocities have been utilized to determine shear (S)-wave velocities in near-surface geophysics since early 1980s. One of the key steps is to calculate theoretical dispersion curves of an earth model. When the S-wave velocity of the surface layer is higher than some of the layers below, however, the Rayleigh-wave phase velocity in a high-frequency range calculated by existing algorithms approaches the lowest S-wave velocity among the layers above the half-space, rather than a value related to the S-wave velocity of the surface layer. According to our numerical modeling results based on wave equation, trends of the Rayleigh-wave dispersive energy approach about a 91% of the S-wave velocity of the surface layer at a high-frequency range when its wavelength is much shorter than the thickness of the surface layer, which cannot be fitted by a dispersion curve calculated by existing algorithms. We propose a method to calculate Rayleigh-wave phase velocities of models with a high-velocity surface layer by considering its penetration depth. We build a substituted model that only contains the layer with the lowest S-wave velocity among the layers above the half-space and the layers above it. We use the substituted model to replace the original model to calculate phase velocities when the Rayleigh-wave wavelength is not long enough to penetrate the lowest S-wave velocity layer. Several synthetic models are used to verify fitness between the dispersion curve calculated by our proposed method and the trend of the highest dispersive energy. Examples of inversion also demonstrate high accuracy of using our method as the forward calculation method during the inversions.     

常温压条件下六种变质程度煤的超声弹性特征

煤弹性是反映煤的物质组分和结构的重要力学特征之一.在钻孔与测井的约束下,运用地震勘探获取煤层的弹性特征以反映其物性等,对于煤炭井工开采和煤层气储层评价及开发具有重要的工程意义,而其中,煤样的超声实验是实现地震反演煤层物性的基础.鉴于此,本文针对中国义马、阜康、淮南、平顶山、鹤壁和焦作6个矿区采集的6种不同变质程度的原煤样30块,在常温常压条件下分别进行了沿煤层走向、倾向及垂直层理三个方向煤样的实验室超声波测量.测试结果显示:煤样纵横波速度在走向、倾向、垂向三个方向上依次减小,存在各向异性,且P波速度的平均各向异性强于S波;品质因子与弹性模量在三个方向上也存在较大差异,且S波的品质因子大于P波的品质因子;弹性模量除泊松比外,均小于一般的沉积岩.通过本实验与分析进一步证明了:Gardener与Castagna公式不适用于中国煤田的煤岩弹性描述,并给出了精度更高的经验公式.     

含流体砂岩地震波频散实验研究

为了研究孔隙流体对不同渗透率岩石地震波速度的影响,在实验室利用跨频带岩石弹性参数测试系统得到了应变幅值10-6的2~2000Hz频段下的地震波速度和1 MHz频率下的超声波速度,利用差分共振声谱法得到了频率600Hz岩石干燥和完全饱水情况下岩石声学参数.实验表明,在低饱和度下,致密砂岩在地震和超声频段下没有明显的频散;在高饱和度下纵波速度的频散变得明显.从干燥到完全水饱和条件,不同频率测量的致密砂岩的体积模量随岩石孔隙度增高而降低,且体积模量的变化量受岩石微观孔隙结构的影响较大.高孔、高渗砂岩无论在低含水度下还是在高含水饱和度下频散微弱,并且在地震频段下围压对于岩石纵横波速度的影响要大于频率的影响.高孔、高渗砂岩和致密砂岩不同含水饱和度下的频散差异可应用于储层预测,油气检测等方面,同时该研究可以更好地帮助理解岩石的黏弹性行为,促进岩石物理频散理论的发展,提高地震解释的精度.     

SEISMIC MODELLING WITH CHANNEL WAVES IN SEAM STRUCTURES INFLUENCED BY MYLONITE ZONES1

Channel waves generated in coal-seams and their reflections from discontinuities are widely used to indicate the tectonic and stratigraphic features of coal deposits, resulting in greater efficiency and safety in coal-mining. In the mining area of Ibbenbüren (F.R.G.) seam structures sometimes contain so-called mylonite zones, which are crushed coal deposits capable of binding gas. If mining hits a mylonite zone this would probably not only reduce output of the mine, but could even cause gas explosions. To investigate the influence of a mylonite zone on the propagation of channel waves, Rayleigh channel wave measurements for 2D analogue models were performed and synthetic seismograms of Love channel waves were calculated. Analogue modelling of the mylonite zone using Rayleigh seam waves within the seam was carried out using a perforation technique. Calculations were made to obtain an estimate of velocity reduction due to perforation. The results agree well with velocity values measured up to a perforation of 25% in a 2D epoxy resin model. Reflected channel wave energy was found by applying dispersion analysis in the case where the impedance reduction between the mylonite seam structure and the undisturbed seam was approximately 5%. The horizontal width of the mylonite structure was detectable from the time lag between reflected channel wave signals from both in-seam borders of the mylonite zone. Resolution of two discrete borders was possible for a width of 1.5 λ's. The influence of a vertical fault, positioned within the mylonite zone, could only poorly be resolved. Numerical model investigations of Love seam waves were concerned mainly with the variation of the horizontal width of the mylonite zone. Mylonite zones with dimensions of the order of 0.4 λ's allow definite statements about their widths from dispersion and spectral analyses. For zones with smaller widths down to 0.2 λ's, it was found that reflectivity and transmissivity analyses give a qualitative criterion for distinguishing a mylonite structure surrounding a fault from a pure fault.     

基于背景噪声研究大别-苏鲁及其邻区的瑞雷波群速度结构

大别-苏鲁造山带是中国大陆东部地区最重要的构造带之一. 为了研究该地区上中下地壳的速度结构, 选用国家数字地震台网和中国区域地震数据台网5省(山东、 安徽、 江苏、 河南和湖北)连续两年(2009年5月—2011年5月)的垂直向地震记录, 进行背景噪声互相关处理, 叠加得到了台站对间的面波经验格林函数. 采用多重滤波法提取了近4000条频散曲线, 并反演得到了研究区10—25 s的瑞雷波群速度分布结果. 通过分析大别-苏鲁及其邻区的瑞雷波群速度结构图像, 发现不同构造块体具有不同的瑞雷波群速度结构: ① 研究区内的郯庐断裂带及其周边地区包括鲁西地块和胶北地块上中下地壳均表现出明显的高速异常, 可能是在拉张环境下岩浆岩上涌, 导致高速、 高密度的变质岩在地壳富集而形成; ② 苏鲁高压变质带的瑞雷波速度在10—25 s周期内明显高于其它地区, 其上中下地壳均表现出较高的群速度结构特征, 认为苏鲁高压变质带至少延伸到下地壳, 而大别造山带在10 s时表现出高速特征, 但在15—25 s没有明显的高速特征, 故无法从其结果中判断大别高压变质带的垂向延伸范围; ③ 华北板块上中下地壳均表现为低速特征, 体现了研究区内华北板块的大陆地壳减薄特征.      

ON THE ABSORPTION-DISPERSION CHARACTERISTICS OF CHANNEL WAVES PROPAGATING IN COAL SEAMS OF VARYING THICKNESS1

The WKB-method is used for the derivation of both the complex dispersion relation and displacement functions for Love channel-waves that propagate in a coal seam of varying thickness. The constant Q-model is used to describe the anelastic friction. With numerical solutions of the absorption-dispersion relation, the influence of thickness changes on the phase velocity and absorption coefficient of Love seam-waves is analysed at various frequencies. It is shown that the changes in the seam thickness can be optimally detected around the average Airy-phase frequency. An equivalence is pointed out between the wave guide structures: homogeneous with varying seam thickness and horizontally inhomogeneous with constant seam thickness.