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

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

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

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

In-mine geoelectric investigations for detecting tectonic disturbances in coal seam structures

The methods of in-mine seam-sounding and transillumination (geoelectric tomography) for the detection of tectonic disturbances of coal seams were developed by the Department of Geophysics of the University of Miskolc in the 1970–80’s with the effective support of the former “Borsod” Coal Mines Ltd. The paper gives an overview about the theory of seam-sounding and a special geoelectric tomographic inversion, and introduces the in-mine geoelectric seam-sounding and transillumination measurement systems using vertical electrode dipoles. In the second part the paper, the results of an in-mine geoelectric measurement are presented, which was carried out in order to detect tectonic disturbances of the Miocene aged coal seams situated in Slovakia. As results of the geophysical investigation, the authors forecasted the tectonic features in the coal seam. The company confirmed the results by independent information about seam disturbances and tectonic features arising from the excavation of the investigated area.     

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波的波速一致.由于这两个触发的波形模式在时间域中具有相似的特征,建议在高频域和低频域同时发育较好的薄煤层槽波勘探中,可以采用在同一时间域中高低频相结合的方法提高利用槽波勘探分析的效果.     

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.     

Fault proximity surveys in coal using guided seismic waves—Numerical modelling

Guided seismic waves are generated at lateral discontinuities in coal seams as an integral part of the mechanics of scattering downgoing body waves from a surface energy source.Numerical and laboratory model studies undertaken for both thin and thick coal seam waveguides show that small faults generate recognizable diffraction patterns by convertingP andS waves into channel wave arrivals at in-seam geophones on a walkawayVSP profile.These results suggest that in-seam geophones can be used during subsequent surface reflection surveys so as to improve mine profitability by optimising longwall layout in faulted ground.     

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.     

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

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

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.     

MODELLING CHANNEL WAVES WITH SYNTHETIC SEISMOGRAMS IN AN ANISOTROPIC IN-SEAM SEISMIC SURVEY1

In-seam seismic surveys with channel waves have been widely used in the United Kingdom and elsewhere to map coal-seams and to detect anomalous features such as dirt bands, seam thinning and thickening, and particularly in-seam faulting. Although the presence of cleat-induced anisotropy has been recognized in the past, almost all previous analyses have assumed homogeneous isotropic or transversely isotropic coal-seams. Channel waves, however, exhibit properties which cannot be fully explained without introducing anisotropy into the coal-seam. In particular, Love-type channel waves are observed for recording geometries where, in a homogeneous isotropic or transversely isotropic structure, the source would not be expected to excite transverse motion. Similarly, modes of channel-wave propagation display the coupled three-component motion of generalized modes in anisotropic substrates, which would not be expected for Rayleigh and Love wave motion in isotropy or in transversely isotropic media with azimuthal isotropy. We model the observed in-seam seismic channel waves with synthetic seismograms to gain an understanding of the effects of cleat-induced anisotropy on the behaviour of channel waves. The results show a reasonable good match with the observations in traveltime, relative amplitudes, dispersion characteristics and particle motions. We demonstrate that anisotropy in the surrounding country rocks contributes significantly to the coupling of channel wave particle motion, although its effect is not as strong as the anisotropy in the coal-seam. We conclude that the effects of cleat- and stress-induced anisotropy are observed and can be modelled with synthetic seismograms, and that anisotropy must be taken into account for the detailed interpretation of channel waves.     

A CROSS-HOLE AND FACE-TO-BOREHOLE IN-SEAM SEISMIC EXPERIMENT AT INVINCIBLE COLLIERY,AUSTRALIA*

An experimental cross-hole and face-to-borehole in-seam seismic survey was carried out at Invincible Colliery in the Western Coalfield of New South Wales. Objectives were to determine propagation characteristics of the Lithgow seam and to establish the infrastructure for seismic mapping hole-to-hole in Australia. The seam supports leaky P-, S- and P-SV-modes. These modes propagate with group velocities (at 60 Hz) of 3.1, 1.5 and 1.2 km/s respectively. Particle motion polarization is well developed, as is dispersion of the SH-mode. Attenuation rates are high. The seam is lossy (Q of approximately 20). Two prominent structures were mapped by mode conversion. One is believed to be a fracture zone, the other a zone of intense roof thrusting. The old workings and a minor strike-slip fault, which intersected raypaths, were found to be relatively transparent to P- and S-waves at 60 Hz. Telemetry delay and shot-break timing errors of the exploder box are significant. The resulting traveltime scatter is reduced by means of a least-squares “statics” procedure. The group velocity estimation algorithm (based on Fourier transform) yields dispersion characteristics which can be matched with theoretical results for a simple model of a coal seam waveguide. The experiment demonstrates the capability to retrieve in-seam seismic data of diagnostic quality over an appreciable distance (2 km). The experience gained in both survey layout and data processing will be beneficial to future seam wave surveying of Australian coal mines.     

Instantaneous polarization filtering focused on suppression of surface waves

The different characteristics of polarization of body and Rayleigh waves make it possible to separate these two types of waves by their characteristics and suppress the latter. The moving time-window analysis often is used in polarization filtering but it is difficult to determine a suitable time-window length, resulting in some problems, such as complex eigenvalues and non-convergence. For overcoming these disadvantages, in this paper, we introduce the concept of complex-trace analysis and conduct de-noise processing to suppress undesirable surface waves by instantaneous polarization analysis in the case of horizontal and vertical component seismic recordings from the Hauinan coal mine. The performance of the method is illustrated by examples with synthetic and field data and its effectiveness to remove surface waves from multi-component seismic data is demonstrated.     

瞬时偏振滤波技术压制X分量中的面波(英文)

瑞雷面波与体波在偏振上的差异使得利用偏振特征分离这两种波场成为可能。为克服滑动时窗法在应用中存在的一些问题,本文将工程地震和天然地震信号分析中的复数道技术与瞬时偏振分析技术引进到能源领域多分量地震勘探转换波记录中的面波压制处理。理论模拟与野外实测数据的处理试验显示了较好的压制面波、保留有效信号的去噪效果,证明了该技术值得在多分量数据处理中推广应用。     

CROSS-BOREHOLE OBSERVATION OF MODE CONVERSION FROM BOREHOLE STONELEY WAVES TO CHANNEL WAVES AT A COAL LAYER1

Conversion of borehole Stoneley waves to channel waves was observed in data from a seismic cross-borehole experiment conducted between wellbores penetrating a thin coal layer at 2022 m depth, near Rifle, Colorado. Traveltime moveout observations show that borehole Stoneley waves underwent partial conversion to channel waves at the coal layer. The channel waves were detected directly in an adjacent borehole 35 m away at receiver positions within the coal. Stoneley waves, subsequently produced by partial conversion of channel waves, were also detected at receiver positions located up to 50 m above and below the coal layer in the adjacent borehole. We infer the channel wave to be the first-higher Rayleigh mode by comparing the observed group velocity with theoretically derived dispersion curves. Identifying the conversion between borehole and stratigraphically guided waves is significant because coal penetrated by multiple wells may be detected without placing a transmitter or receiver in the coal itself.     

SEISMIC MODELLING OF SEAM WAVES EXCITED BY ENERGY TRANSMISSION INTO A SEAM1

The traditional method of exciting channel waves in coal deposits underground consists of firing explosive sources in a mid-seam position generating seam waves of the Rayleigh and Love type. We investigate various source positions and excitation mechanisms within the bedrock structure surrounding the seam to evaluate the effects of source positions adjacent to the seam. The investigation is based on analogue and numerical modelling of half- and full-space cases, for which the excitation and the nature of Rayleigh channel waves are examined. In the analogue modelling, sources, located from mid-seam out into the bedrock, along the edge of a 2D plate model, excited channel waves through a conversion of the free surface Rayleigh wave at the edge of the plate. The excited channel wave belongs to the normal mode range. Frequency-wavenumber analysis shows that the symmetric 2nd mode of the channel wave is excited with frequencies comparable to the forcing frequency of the source signal. The polarization changes from retrograde to prograde, as the wave develops from the front to the rear of the seam, respectively. The amplitude-depth distribution resembles that of an ordinarily excited seam wave, for the symmetric component. However, the antisymmetric component does not show the characteristic change of sign in amplitudes across the mid-seam axis. Numerical modelling with sources located in the bedrock (full-space case) shows that relocating the source away from the seam lowers the frequency content of the excited channel wave. Based on these investigations, the influence of a lower-frequency source signal on the excitation of the channel wave is examined in an analogue experiment. Sources are sited in the bedrock adjacent to the seam at three locations. A lower-frequency wavelet is calculated for each source location from the results obtained in the numerical analysis. For comparison, a higher-frequency wavelet is also used which is known to be optimal for this model geometry when excited by a mid-seam source location. It is found that in two cases the use of the lower-frequency wavelet improves the channel wave excitation, while no amplification is achieved in one case.     

Dispersion features of transmitted channel waves and inversion of coal seam thickness

In-seam seismic survey currently is a hot geophysical exploration technology used for the prediction of coal seam thickness in China. Many studies have investigated the relationship between the group velocity of channel wave at certain frequency and the actual thickness of exposed coal beds. But these results are based on statistics and not universally applicable to predict the thickness of coal seams. In this study, we first theoretically analyzed the relationship between the depth and energy distribution of multi-order Love-type channel waves and found that when the channel wave wavelength is smaller than the thickness of the coal seam, the energy is more concentrated, while when the wavelength is greater than the thickness, the energy reduces linearly. We then utilized the numerical simulation technology to obtain the signal of the simulated Love-type channel wave, analyzed its frequency dispersion, and calculated the theoretical dispersion curves. The results showed that the dispersion characteristics of the channel wave are closely related to the thickness of coal seam, and the shear wave velocity of the coal seam and its surrounding rocks. In addition, we for the first time realized the joint inversion of multi-order Love-type channel waves based on the genetic algorithm and inversely calculated the velocities of shear wave in both coal seam and its surrounding rocks and the thickness of the coal seam. In addition, we found the group velocity dispersion curve of the single-channel transmitted channel wave using the time–frequency analysis and obtained the phase velocity dispersion curve based on the mathematical relationship between the group and phase velocities. Moreover, we employed the phase velocity dispersion curve to complete the inversion of the above method and obtain the predicted coal seam thickness. By comparing the geological sketch of the coal mining face, we found that the predicted coal seam thickness is in good agreement with the actual thickness. Overall, adopting the channel wave inversion method that creatively uses the complete dispersion curve can obtain the shear wave velocities of the coal and its surrounding rocks, and analyzing the depth of the abruptly changed shear wave velocity can accurately obtain the thickness of the coal seam. Therefore, our study proved that this inversion method is feasible to be used in both simulation experiments and actual detection.     

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.     

DISPERSION CALCULATIONS FOR SH AND P-SV WAVES IN MULTILAYERED COAL SEAMS*

Recent work by Räder, Schott, Dresen and Rüter has provided a means of calculating dispersion and amplitude distribution curves for SH channel waves in multilayered media. An alternative calculation procedure is presented and the P-SV wave type is also included. Numerical results are given for a coal seam containing a band of dirt of either higher or lower seismic velocities than the coal itself.     

THE INFLUENCE OF OBLIQUELY DIPPING DISCONTINUITIES ON THE USE OF RAYLEIGH CHANNEL WAVES FOR THE IN-SEAM SEISMIC REFLECTION METHOD *

The geological sequence rock-coal-rock represents a seismic low-velocity channel. Channel waves generated in a coal seam and their reflections from discontinuities can be used for proving the minability of the seam. To investigate the process of reflection, two-dimensional models of the sequence rock-coal-rock have been investigated by means of the ultrasonic transducer technique. Two-dimensional models have the advantage that the wave field can be observed at an arbitrarily chosen point of the model plate. Thus, by means of these models the direct and the reflected Rayleigh channel wave can be observed along their path of propagation. From the geophysical point of view the various types of discontinuities of a coal seam can be divided into two basic types: one is restricted to the seam, the other includes in addition a fractured zone in the adjoining rock. The investigation of the symmetrical Rayleigh channel wave reflected by a discontinuity yields the following results: For dip angles γ between approximately 90° and approximately 60° the reflectivity is virtually independent of the type of discontinuity. This does not hold for the limiting case of γ= 90° (vertical dip) for which the reflectivity increases with increasing influence of the fractured zone. For dip angles γ between approximately 60° and approximately 40° the reflectivity is still independent of the fractured zone but the shape of the reflected wave deteriorates with increasing influence of the fractured zone. For dip angles γ below approximately 40° the reflected wave deteriorates such that the application of the in-seam seismic reflection method will be difficult or even impossible. The conversion of the direct wave of the symmetrical fundamental mode into a reflected wave of the antisymmetrical fundamental mode has been observed.