THE SCATTERING OF SH-CHANNEL WAVES BY A FAULT IN A COAL SEAM*

The imaging of faults in coal seams by the in-seam seismic method has now become standard practice. In the UK over 300 surveys have been undertaken and the technique is now part of the exploration arsenal of colliery planners. From these users comes the pressure for two major improvements, namely an increase in range and target identification. This paper is directed towards the latter problem. It has long been recognized that the reflected channel waves must contain information on the fault structure that caused the reflection, and model experiments have been undertaken to investigate the reflection process. Only recently, however, have attempts been made to quantify the reflection process. Calculations using both the finite-difference and finite-element techniques have been carried out, and estimates of the reflection coefficient as a function of frequency have been obtained. The object of this paper is to extend these considerations by calculating analytically the scattering matrix of an SH-channel wave after interaction at a fault plane. The scattering matrix is calculated as a function of frequency, hade angle, and fault throw. The method employed is based on the decomposition of the incident SH-channel wave into Fourier components, the calculation of plane wave reflection and transmission coefficients within the constraints of geometrical acoustics, and finally the synthesis of the scattering matrix by application of the Helmholtz-Kirchhoff integral. The calculation throughout is restricted to normal modes.     

LOVE-TYPE SEAM-WAVES IN WASHOUT MODELS OF COAL SEAMS*

The propagation of Love seam-waves across washouts of coal seams was studied by calculating synthetic seismograms with a finite-difference method. Seam interruption, seam end and seam thinning models were investigated. The horizontal offset, the dip of the discontinuities and the degree of erosion served as variable parameters. Maximum displacement amplitudes, relative spectral amplitudes and phase and group slowness curves were extracted from the synthetic seismograms. Both seam interruption and seam thinning reduce the maximum displacement amplitudes of the transmitted Love seam-waves. The degree of amplitude reduction depends on the horizontal offset and the degree of erosion. It is four times greater for a total seam interruption than for an equivalent seam thinning with a horizontal offset of four times the seam thickness. In a seam cut vertically, the impedance contrast between the coal and the washout filling determines the maximum displacement amplitudes of the reflected Love seam-waves. They diminish by a maximum factor of four in oblique interruption zone discontinuities with a dip of maximum 27°, and by a maximum factor of ten in a seam thinning with a degree of erosion of at least 22%. The analysis of the relative spectral amplitudes indicates a preferential transmission of those phases with frequencies below, and a preferential reflection of those phases with frequencies above the first mode Airy-phase. The relative spectral amplitudes of the reflected Love seam-waves show a distinct interference pattern of the waves reflected at both interruption zone discontinuities. The dispersion analysis reveals a flattening of the phase and group slowness curves with increasing frequencies, horizontal offset and degrees of erosion. These results imply that a detection of washouts in-mine will be possible in a frequency range including at least the first mode Airy-phase. An interference pattern and a flattening of the dispersion curve indicate a washout rather than other seam obstructions and leads to an estimate of the washout dimension.     

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

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

Applying the tilt‐depth and contact‐depth methods to the magnetic anomalies of thin dykes

A new method for the calculation of the depth, location, and dip of thin dykes from pole‐reduced magnetic data is introduced. The depth can be obtained by measuring the distance between chosen values of a tilt angle that is based upon the ratio of the magnetic field and its Hilbert transform over the dyke. Alternatively, it can be obtained from the horizontal derivative of the ratio of the Hilbert transform of the field to the field itself, over the dyke. The latter method also allows the dip of the dyke to be estimated from the gradient of the depth estimates.     

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.     

地震反射法中薄煤层分辨能力的研究

本文以非线性地震层状构造理论为基础,对薄层分辨率进行了研究,研究中采用理论计算和合成地震记录,并用超声地震模型试验加以验证。结果表明,薄层可检测性准则,应是它的复合反射波的振幅不小于围岩中较强反射波的振幅。由Farr提出的可检测薄层的极限厚度为λ/12（λ为薄层中地震波波长）,是在薄层顶、底界面反射系数很小时的近似值。指出,当薄层顶、底界面反射系数较大时,如煤层,其可检测性高于Farr定义,且随界面反射系数的增大或频率的提高而提高;适当地提高地震反射波的频率,就能得到厚2-3m薄煤层的大振幅、高信噪比的反射波。上述结论均用煤田地震勘探实例证实。     

Finding the strike direction of fractures using GPR

GPR reflection energy varies with antenna orientation relative to the strike and dip of the reflector. This directional dependence of GPR responses was investigated through numerical experiments and was used to estimate the azimuth of fractures and joints. Three antenna configurations were considered in this study: perpendicular-broadside (YY mode), parallel-broadside (XX mode) and cross-polarization (YX mode). The reflection energy in the cross-polarization mode shows a shape characteristic similar to the strike, regardless of the dip angle. Those in the other two modes show quite different amplitudes from the strike, depending on the dip angle.  We have developed a strike-direction-finding scheme using data obtained from the three different modes for the same survey line. The azimuth angle of each reflector was displayed in colour on the GPR profile. This scheme was applied to a field survey at a granite quarry in southern Korea. The GPR profiles showed different images of the reflectors depending on the antenna configuration. The estimated azimuths of reflectors obtained using our scheme matched fairly well with those of known fractures and joints.     

三河-平谷8.0级地震区地壳结构和活动断裂研究——利用单次覆盖深反射和浅层地震剖面

跨1679年三河-平谷8.0级地震区完成的单次覆盖深地震反射剖面和浅层地震反射剖面,揭示了三河-平谷地震区的地壳结构和断裂的深、浅构造特征.结果表明,该区地壳以TWT6～7 s左右的强反射带为界分为上地壳和下地壳,上地壳厚约18～21 km,下地壳厚约13～15 km.剖面揭示的地壳深断裂和浅部活动断裂具有上下一致的对...     

Surface stresses associated with arrested dykes in rift zones

Many theoretical models predict that arrested dykes may generate major grabens at rift-zone surfaces. Arrested dyke tips in eroded rift zones, however, are normally not associated with major grabens or normal faults that could be generated by dyke-induced stresses ahead of the tips, and normal faults and grabens tend to be less common in those parts of eroded rift zones where dykes are comparatively abundant. Similarly, there are feeder dykes, as well as dykes arrested a few metres below the surface, that do not generate faults or grabens at the surface. Here I propose that this discrepancy between theoretical models and field observations may be explained by the mechanical layering of the crust. Numerical models presented here show that abrupt changes in Young's moduli, layers with high dyke-normal compressive stresses (stress barriers), and weak, horizontal contacts have large effects on the dyke-induced stress fields. For the models considered, the surface tensile stresses induced by arrested dykes are normally too small to lead to significant fault or graben formation at the rift-zone surface. The only significant dyke-induced surface tensile stresses (2 MPa) in these models are for a dyke tip arrested at 1 km depth below the surface of a rift zone with a weak contact at 400 m depth and subject to extension. That tensile stress, however, peaks above the ends of the weak horizontal contact, which, in the model considered, occur at distances of 4 km to either side of the dyke, and shows no simple relation to the depth to the dyke tip. Thus, for a layered crust with weak contacts, straightforward inversion of surface geodetic data to infer dyke geometries may result in unreliable results.Editorial responsibility: A. Woods     

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

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

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.     

低纬哨声通过低电离层传播的透射特征——全波解

本文利用全波解方法研究了低纬哨声通过低电离层传播的过程。文中主要计算了哨声波的透射损失和反射系数随投射角与方位角的变化,并结合哨声多站同时观测结果进行了简要的讨论。全波解数值计算结果表明,透射损失随投射角的变化具有南北不对称性,当波矢沿磁力线方向时透射损失最小,且这种不对称性随磁倾角的减小而明显加剧。计算结果还表明,下行哨声的内反射系数基本上不随投射角而变化;但随波频率的升高而减小,并伴有数个相间的极大和极小值,这可能导致哨声回波频谱畸变。     

Analytical study of the reflection and transmission coefficient of the submarine interface

The analytical study of the reflection and transmission coefficient of the seafloor interface is essential for the characterization of the ocean bottom in marine seismic exploration. Based on the boundary conditions of the seafloor interface, the analytical expression of the reflection and transmission coefficient at the submarine interface is derived in this study by using the steady-state wave solution of the elastic wave in a homogeneous, isotropic medium. With this analytical expression, the characteristics of the reflection and transmission coefficient at the submarine interface are analysed and discussed using critical angles. The results show that the change in the reflection and transmission coefficient with the incidence angle presents a “segmented” characteristic, in which the critical angle is the dividing point. The amplitude value and phase angle of the coefficient at the submarine interface change dramatically at the critical angle, which is related to the P- and S-wave velocities in the seabed layer. Compared with the stiff seabed, the soft seabed has a larger P-wave critical angle and an absence of the converted S-wave critical angle, owing to the low P- and S-wave velocities in the solid seabed layer. By analysing and discussing the special changes that occur in the coefficient values at the critical angle, the reflection and transmission characteristics of the different incident angles are obtained. Synthetic models of both stiff and soft seafloors are provided in this study to verify the analytical results. Finally, we compared our synthetic results with real data from the Gulf of Mexico, which enabled the validation of our conclusions.     

Seismic Rheological Model and Reflection Coefficients of the Brittle–Ductile Transition

It is well established that the upper—cooler—part of the crust is brittle, while deeper zones present ductile behaviour. In some cases, this brittle–ductile transition is a single seismic reflector with an associated reflection coefficient. We first develop a stress–strain relation including the effects of crust anisotropy, seismic attenuation and ductility in which deformation takes place by shear plastic flow. Viscoelastic anisotropy is based on the eigenstrain model and the Zener and Burgers mechanical models are used to model the effects of seismic attenuation, velocity dispersion, and steady-state creep flow, respectively. The stiffness components of the brittle and ductile media depend on stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the reflection and transmission coefficients of a single brittle–ductile interface and of a ductile thin layer. The PP scattering coefficient has a Brewster angle (a sign change) in both cases, and there is substantial PS conversion at intermediate angles. The PP coefficient is sensitive to the layer thickness, unlike the SS coefficient. Thick layers have a well-defined Brewster angle and show higher reflection amplitudes. Finally, we compute synthetic seismograms in a homogeneous medium as a function of temperature.     

Reflection and transmission of plane SH-waves in an initially stressed inhomogeneous anisotropic magnetoelastic medium

This study investigates the reflection and transmission of plane SH-waves in two semi-infinite anisotropic magnetoelastic media. The lower half-space is considered as initially stressed and inhomogeneous. The density of lower half-space is taken exponentially varying with depth. The solutions for half-spaces are obtained analytically. The expressions for reflection and transmission coefficient are obtained in the closed form subject to continuity conditions at the interfaces of anisotropic magnetoelastic half-spaces and the Snell’s law. It is found that these coefficients depend on the initial stress, inhomogeneity parameter, the magnetoelastic coupling parameter, and the angle at which wave crosses the magnetic field of the half-spaces. Numerical computations are performed for these coefficients for a specific model of two different anisotropic magnetoelastic half-spaces. The numerical results are illustrated by the graph of reflection and transmission coefficient versus the angle of incidence. In general, as the initial stress increases the reflection and transmission coefficient increases, the affect is more prominent for more than 10 GPa. Inhomogeneity in the density of the material also increases the reflection and transmission coefficient. The anisotropic magnetoelastic parameter and the angle at which the wave crosses the magnetic field for both the half-spaces have a quite significant effect on the reflection and transmission coefficient.     

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

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

Detection of Magnetically Susceptible Dyke Swarms in a Fresh Coastal Aquifer

Groundwater constitutes the main source of freshwater in Shalatein, on the western coast of the Red Sea, in Egypt. The fresh aquifer of Shalatein is intensively dissected by shallow and deep faults associated with the occurrence of dykes and/or dyke swarms. In this context, synthesis of electrical resistivity, ground magnetics, and borehole data was implemented to investigate the freshwater aquifer condition, locate the intrusive dykes and/or dyke swarms, and demarcate the potential freshwater zones. Nine Schlumberger VES’s with maximum current electrode half-spacing (AB/2) of 682 m were conducted. The subsurface was successfully delineated by general four layers. The fresh aquifer of the Quaternary and Pre-Quaternary alluvium sediments was effectively demarcated with true resistivities ranged from 30 to 105 Ωm and thickness ranged between 20 and 60 m. A ground magnetic survey comprised 35 magnetic profiles, each 7 km in length. Magnetic data interpretation of the vertical derivatives (first and second order), downward continuation (100 m), apparent susceptibility (depth of 100 m), and wavelength filters (Butterworth high-pass of wavelengths <100 m and Band-Pass of wavelengths 30–100 m) successfully distinguished the near surface structure with five major clusters of dyke swarms, whereas filters of the upward continuation (300 m) and Butterworth low-pass (wavelengths >300 m) clearly reflected the deep-seated structure. The computed depth by the 3D Euler deconvolution for geological contacts and faults (SI = 0) ranged from 14 to 545 m, whereas for dyke and sill (SI = 1), it ranged from 10 to 1,095 m. The western part of the study area is recommended as a potential freshwater zone as it is characterized by depths >100 m to the top of the dykes, higher thickness of the fresh aquifer (45–60 m), depths to the top of the fresh aquifer ranging from 25 to 40 m, and higher resistivities reflecting better freshwater quality (70–105 Ωm).     

Exploration and Study of Deep Crustal Structure in the Quanzhou Basin and Its Adjacent Area

INTRODUCTIONQuanzhou Citylocates on the southeast coast of Fujian Province .Tectonically,it locates in thevolcanic fault depression zone in East Fujian betweenthe Wuyi-Daiyunfolded doming-upregion of theSouth China block and the depression zone of the Taiwan Straits . Being on the north segment of thesoutheast seismic belt of China , many destructive earthquakes inthe history affected the area and theregional seismicityis closelyrelated withfaulting.Inthe Quanzhou basin,large-scalelo…     

武清-北京-赤城二维地壳结构和构造

对武清—北京—赤城剖面段的人工地震测深资料,利用常规地震测深方法进行处理与计算,目的是查明该地区二维地壳速度结构和构造的基本特征、断层的空间展布特征和深浅部地质构造环境等.结果表明,不同区域的地壳速度结构在纵向和横向上均具有明显的非均匀性.地壳呈层状结构,地壳厚度由武清的31.0 km向北西至延庆逐渐加厚为40.0 km.在北京下方壳内界面起伏较大,莫霍面在北京的北西方向有一急剧加深,在约20 km范围内M面起伏变化达5.0 km. 根据地震波动力学与运动学特征,以及二维地壳速度等值线和地震界面起伏变化等特征,沿剖面推测了5条地壳浅部断裂.     

基于地震活动性资料估计海原断裂倾角

本文以海原断裂带区域活动构造为基础,将海原断裂划分为西、中、东3段.基于1999年12月26日至2010年7月26日间的精定位小震目录,估计了海原断裂带各段的倾角.考虑海原断裂沿走向可能存在南倾与北倾两种情况,将倾角的范围设置为.首先运用网格搜索法确定了平面断层模型,其次以特征深度节点为基础数据运用多项式构建了曲面断层模型.结果表明：当进行平面拟合时,海原断裂西段与中段、东段的倾向不同,西段为南倾,其倾角值为71°,而中段、东段为北倾,其倾角值分别为72°、65°,各段的倾角值均由地表以下8 km地震资料确定.当进行曲面拟合时,在8 km深度以内海原断裂西段、中段、东段的倾角均处于80°左右,即接近陡立.西段的倾角在深度为9 km处出现转换,之后倾角接近陡立;中段的倾角在深度为16 km处出现转换,之后倾角逐渐减小,当深度为18 km时倾角为30°;东段的倾角在深度为11 km处出现转换,倾角为42°,在深度为16 km处出现第二次转换,倾角为55°,之后倾角逐渐减小.结合震源机制解和大地测量观测资料反演拟合的合理性,验证了本文所估计倾角的可靠性.