含采空区层状岩体微震信号传播规律试验研究

微震信号在层状岩体中传播规律的研究对于实现震源准确定位具有重要的现实意义。在天然岩体中,微震信号在岩层中传播的速度,不仅受其内部因素即岩石自身物理性质的影响,还受外部因素如地质结构面、断层及采空区等条件的影响。通过室内试验的方法对微震信号在不同岩层及采空区处的传播规律进行研究,结果表明：波速随着传播距离的增加而逐渐衰减,且岩层密度越大衰减越慢;穿过的结构面数量越多,微震信号传播速度衰减比例越大;能量的衰减与波速的衰减一致。采空区断面横截面积越大,微震信号通过的用时越长;采空区周围岩石密度、弹性模量等物理性质的数值越大,微震信号通过时衰减越少;边际谱分析结果表明：断层对频率高的微震信号阻碍作用更强。     

含采空区层状岩体微震信号传播规律试验研究

微震信号在层状岩体中传播规律的研究对于实现震源准确定位具有重要的现实意义。在天然岩体中,微震信号在岩层中传播的速度,不仅受其内部因素即岩石自身物理性质的影响,还受外部因素如地质结构面、断层及采空区等条件的影响。通过室内试验的方法对微震信号在不同岩层及采空区处的传播规律进行研究。结果表明：波速随着传播距离的增加而逐渐衰减,且岩层密度越大衰减越慢;穿过的结构面数量越多,微震信号传播速度衰减比例越大;能量的衰减与波速的衰减一致。采空区断面横截面积越大,微震信号通过的用时越长;采空区周围岩石密度、弹性模量等物理性质的数值越大,微震信号通过时衰减越少。边际谱分析结果表明：断层对频率高的微震信号阻碍作用更强。     

基于变分模态分解的采空区“三带”微震信号能量衰减规律

为探明微震信号能量在采空区“三带”结构中的衰减规律,拟开展采空区覆岩相似模型试验,采集人工激发微震波经由采空区结构传播的微震信号,通过变分模态分解（variationalmodedecomposition,VMD）处理微震信号,获取各频率下模态分量。针对采空区微震信号在VMD下各模态分量中心频率与能量之间的关系展开分析。根据中心频率法确定微震信号最佳模态数量,并计算微震信号欠分解状态、最佳分解状态、过分解状态下各分量能量;对各震源下信号最佳分解状态时各模态分量能量与中心频率分布关系进行拟合,分析在“三带”结构中,微震信号不同传播状态下各结构层对信号能量影响作用。研究结果表明：（1）在VMD过程中,人工激发震动信号有效模态数量在6～11范围内,微震信号能量随模态数量变化明显。（2）采用幂函数可实现对微震信号模态能量与频率关系的拟合,且拟合状态良好（决定系数大于0.9）,其中低频模态分量包含能量占信号总能量近50%;采用高斯函数可以拟合震源各分量能量在频域上的分布表现,拟合状态较好,且表现出高斯单峰特征。（3）微震信号穿越采空区“三带”结构,微震信号能量随震源位置与传感器距离增加而减小,同...     

粘弹性介质地震波传播特征及反射特征研究

理想弹性介质并不能解决许多复杂的实际问题。在实际介质中,地震波的传播过程存在吸收作用,使地震波能量发生了损耗、频带变窄、相位延迟、频率降低,尤其是高频部份。这种特性可以用传播矢量和衰减矢量来表示,衰减矢量越大,介质的吸收越强。这里通过引入广义平面波、纵横波品质因子以及衰减角,利用应力应变关系和连续性条件,推导了粘弹性介质精确的Zoeppritz方程,并通过求解精确的Zoeppritz方程,分析了粘弹性介质中传播矢量、衰减矢量、纵波反射系数、转换波反射系数与频率、衰减角、入射角、品质因子、速度等物理量之间的关系,详细地了解了粘弹性介质中地震波传播特征及反射特征,为粘弹性介质中储层预测问题的研究提供了理论基础。     

黏弹介质槽波超前探测数值模拟与波场分析

目前常用槽波超前探测掘进工作面前方隐伏断层,而实际煤层具有黏弹性,对槽波有吸收衰减作用。为研究黏弹煤层介质中槽波超前探测的波场特征和传播规律,基于Kelvin-Voigt一阶速度-应力方程,建立了三维含断层煤系地质模型,进行了三维有限差分数值模拟。结果表明：黏弹介质中反射槽波能量衰减大,其传播与衰减特征比完全弹性介质更符合实际煤层反射槽波传播特征。煤层Q值越大槽波超前探测效果越好,煤层Q值过小时,槽波超前探测效果差,可利用x分量的反射横波与z分量的P-S波超前探测。断层落差小于煤厚时,y和z分量的槽波超前探测效果都较好;落差大于煤厚时,可利用y分量进行槽波超前探测。断层夹角大于60°时,可利用y和z分量进行槽波超前探测;夹角小于60°时,可利用x和y分量进行槽波超前探测。研究结果可为煤矿槽波超前探测提供理论支持。     

节理面透射模型及其隔振性能研究

振动波入射到节理界面上,与软弱断层一样会产生波场的分解,能够起到隔振作用。根据线性变形节理的透反射规律研究发现,隔振后各点的波场可由透射P波与透射S波叠加而成,其能量分配服从Zoeppritz方程,不连续位移造成了透射波的衰减;法向刚度增加透射纵波系数会增加,而切向刚度与透射纵波无关,法向与切向刚度对透射S波均有影响,刚度增加透射S波振幅逐渐降低,逐步逼近于0;由于透射后波场相互干涉,不同波场延时不同,透射合成波与入射波相比,振幅降低,持续时间增加,波形尾段有畸变;总能量系数与节理面的刚度无关,而随入射角增加不断降低。     

基于粘滞性声波方程的吸收补偿方法

通常地层介质都是粘滞性的,由于地层的吸收作用,地震波在向下传播过程中能量被衰减,特别是高频能量的衰减,使反射信号频带宽度变窄。为了补偿被吸收的能量,笔者采用粘滞声波方程波场延拓方法进行能量吸收和衰减补偿。首先,从粘滞性波动方程波场延拓函数理论出发,推导了波场延拓函数稳定性条件;其次,对其有效性进行了测试和分析;最后,通过设计数值模型进行正演模拟对比试验,证明了基于粘滞性声波方程的吸收补偿方法能够很好地对粘性介质吸收衰减进行能量补偿,改善处理效果。     

陡倾结构金属矿山采空区围岩破坏机制研究

在地下金属矿山开采中会造成采空区围岩变形破坏,影响地下采矿安全。采用现场调查、位移监测、微震监测,并结合理论分析等方法对金山店东区地下采空区的围岩破坏机制进行研究。考虑水平构造应力和崩落岩体的影响,基于极限平衡理论和能量法分别得到上下盘采空区围岩破坏机制,即倾倒滑移破坏和溃曲滑移破坏。上盘和下盘围岩在陡倾结构面的切割作用下,分别形成反倾结构和顺倾结构。在水平构造应力作用下,上盘陡倾结构岩体向采空区侧发生倾倒破坏,岩体变形破坏进而导致F4断层活化。随着采矿的深入,破坏的岩体数量和范围不断增加,形成一个深部滑动面,并穿过F4断层。下盘采空区围岩在水平构造应力作用下诱发了矿体边界处F1断层滑移,同时引起顺倾结构岩柱发生溃曲-滑移破坏,沿着节理面形成滑动面。     

不同饱和度土层分界面上剪切波的反射与透射

基于多相孔隙介质弹性理论,给出了非饱和土中不同弹性波的传播方程。根据分界面上的边界条件,建立了各势函数波幅值之间的关系式,讨论了入射剪切波在不同饱和度土层分界面上的反射与透射问题。在无限空间非饱和土体中存在3种压缩波和1种剪切波,因此,当剪切波传播到不同饱和度的非饱和土层分界面上将分别在上、下土层激发产生4种反射波和4种透射波。推导出不同反射波和透射波的振幅比例系数和能量比例系数的理论表达式,并且在此基础上进行数值分析。在数值算例中分别研究了各反射波与透射波的能量比例系数（即能量反射率和能量透射率）受入射频率、入射角度以及上、下土层土体饱和度变化的影响情况。计算结果表明：各能量反射率和能量透射率不仅与入射角和入射频率有关,而且其受上、下土层饱和度变化的影响也同样不能忽视。     

应力波斜入射对岩体节理端部的影响

岩体节理在受到应力波扰动时端部受力会发生不同程度的连续性变化,为具体分析这种动态变化与应力波入射角度之间的关系,引用岩石非线性法向本构关系与线性切向本构关系的组合模型以及相应的在P波斜入射节理的应力波传播方程,结合岩石断裂力学中Ⅰ、Ⅱ型裂纹尖端应力和位移场的相关理论,得出组合形式下的节理端部应力场和位移场随质点振速变化的计算公式。通过不同入射角的模拟脉冲信号作用对端部应力位移的变化分析：应力波斜入射节理导致节理端部上、下两侧应力与位移场非对称分布,随着入射波质点振速的增减变化,应力集中位置会出现变化;计算并整理模型中节理端部上、下两侧0.005 m位置的数据,节理法向刚度由入射波质点振速带来的变化直接影响到应力波的透射与反射,进而会导致端部的应力与位移出现滞后效应;节理端部横向位移值与入射角度并非单调变化,而竖向位移会随着入射角的增大呈现下降趋势。     

Analysis of Stochastic Seismic Wave Interaction with a Slippery Rock Fault

Stochastic seismic wave interaction with a slippery rock fault is studied, based on the principle of conservation of momentum at the wave fronts along the fault. By using the displacement discontinuity method, the wave propagation equations are derived for incident longitudinal-(P-) and shear-(S-) waves, respectively. This is an extension of the study by Li and Ma (2010) for blast-induced wave propagation across a linear rock joint. Stochastic seismic waves are generated from a frequency spectrum and used to analyze the seismic wave interaction with a rock fault having a Coulomb-slip behavior. Parametric studies are carried out to investigate the effect of the intensity and impinging angle of the incident seismic waves on wave propagation across a slippery rock fault. Results show that the transmission of the incident P-wave is almost not affected by the fault, on the contrary, this is not the case for an incident S-wave, due to the occurrence of a relative slip which is related to the impinging angle of the incident S-wave. A quantitative study is presented which is of help in understanding the propagation and attenuation laws of seismic waves in discontinuous rock masses.     

采动微地震波传播与衰减特性研究

为分析采动地震波煤岩响应及其波场特性,构建顶板岩层和煤层2种传播介质接收采动地震波的多通道微震监测系统,研究2种介质采动地震波的传播和衰减特征。研究结果表明：顶板岩层接收到采动地震波的平均波速、加速度均高于煤层,2种介质中的地震波加速度均与爆破装药量成正相关,与传播距离、介质密度负相关。采动微地震波能量、信号时长随传播距离分别呈指数、线性衰减,近爆心处煤层传播微地震波能量、信号时长均高于顶板,超过一定距离后顶板高于煤层;采动微地震波传播距离有限,顶板的有效传播距离高于煤层,地震波有效传播距离与爆破装药量线性正相关。采动地震波低频部分的品质因子高于高频部分,顶板岩层的品质因子高于煤层。研究成果有助于地震波传播介质选取、检波器布设优化。     

南川地区溶洞及采空区地震资料针对性处理方法

为解决地震勘探中地下溶洞及采空区位置地震资料成像效果差的难题,利用正演模拟结合实际资料分析,总结该类区域地震资料的频率、能量和信噪比特征,在此基础上,提出一套针对性的处理方法,包括微测井约束层析静校正、弱信号提取与补偿、叠前五维数据规则化和基于构造约束的网格层析速度建模。通过在南川地区的应用,地震资料信噪比和波组连续性得到明显改善。证明了此方法可以有效提高地震勘探中溶洞及采空区位置地震资料品质。      

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Regional hazard prediction of rock bursts using microseismic energy attenuation tomography in deep mining

Rock burst prediction is a worldwide challenge that we have long tried to overcome. This study tentatively proposed a method to regionally predict rock burst hazards using microseismic energy attenuation. To verify the feasibility of the proposal, first, the mechanism of microseismic energy propagation and attenuation in rock medium was explored, and dominant attenuation characteristics of microseismic waves were analyzed. Second, a spatial attenuation model of microseismic energy was established, and the average energy attenuation coefficient for each wave path was defined. A 3D seismic energy attenuation inversion algorithm was put forward, and the corresponding computation matrix was developed. Third, a continuous microseismic field investigation was carried out in a deep coal mine. Seismic energy attenuation coefficient was confirmed using the calibrated focus position and energy determination. Based on data discretization processing, energy attenuation inversion and tomography, potential rock burst hazard regions were strictly zoned in mining areas. Finally, regional prediction results obtained from the microseismic energy attenuation were compared with the direct measurement results obtained from the classical drilling dust method to verify the reliability of proposed approach. It turns out that rock burst hazard regions predicted by the microseismic energy attenuation agreed well with the objective hazardous situations. Seismic energy attenuation coefficient is a significant evaluation factor that directly mirrors the inelastic performance of rock medium. Energy attenuation coefficient threshold used for determining the rock burst hazard regions was 3.0 km^?1. Reliability of the seismic energy attenuation inversion and tomography was closely related to the spatial distribution of microseisms in a localized region. The optimum spatial density of microseisms was 0.2 m^?3. Regional rock burst prediction using microseismic energy attenuation is an effective approach for revealing potential hazardous regions in deep mining conditions. This approach improves the pertinence of geological hazard prevention and provides a beforehand reference for targeted hazard management.     

复杂断层构造电磁波衰减系数层析成像模拟研究

传统矿井电磁波透视层析成像技术只考虑电磁波的传播距离而忽视传播的方向性,造成探测复杂构造时的层析成像分辨率较低。依据电磁波在煤岩介质中的传播规律推导出了与入射角有关的新振幅衰减常数。引用新振幅衰减常数并基于代数重建算法建立了层析成像数学模型,编制了程序软件。利用新振幅衰减常数层析成像算法分别对平行和交叉断层进行数值模拟。通过与传统衰减常数层析成像的效果进行对比,结果表明:改进后的新振幅衰减常数层析成像效果更好,对断层的位置、走向、延展长度等情况反映清晰准确;在介质背景衰减场不易确定的情况下,绝对衰减层析成像方式优于相对衰减层析成像方式;平行关系的断层构造分辨率要高于交叉关系的断层构造。通过在某采煤工作面的电磁波层析成像中成功圈定了一组隐伏的平行断层构造实例,说明采用新衰减常数的层析成像技术效果提高明显。     

岩体夹层应力波能量演化及应力响应特征分析

夹层是常见的地质结构,在地震或爆破荷载作用下,夹层对应力波的传播以及对岩体的响应具有重要的影响。以往对应力波在夹层介质中传播的研究多集中于夹层对应力波的隔振或透射性能,而对应力波在夹层中的多次折、反射过程中能量的演化规律缺乏讨论,对夹层介质的应力响应与破坏没有开展较好地分析。通过理论方法对应力波在夹层内部传播过程中的能量系数的变化规律进行了研究,分析了岩体介质波阻抗和应力波入射角对夹层内外介质中累积波动能量系数的影响规律,以及平面型边坡中软弱夹层的应力响应特征和动态稳定性。研究发现,应力波在夹层内部往复反射过程中,夹层内剩余应力波能量随折、反射发生次数呈指数曲线下降,第4次折、反射后产生的应力波能量可以忽略;夹层内外介质中应力波的累积能量系数的差异随着夹层内外介质波阻抗的相对差异的增大而增大。在平面谐波入射下,边坡内部的夹层中的剪应力和抗剪强度呈波动变化;相对P波,SV波入射会产生较高水平的剪应力,对边坡稳定性影响最大。且SV波入射时,边坡的安全系数对夹层的倾角变化更为敏感,随倾角增大而迅速降低。     

Numerical study on S-wave transmission across a rough,filled discontinuity

This paper presents a numerical simulation of S-wave propagation across a rough, filled discontinuity using the universal distinct element code (UDEC). The ability of UDEC to simulate a stress wave across a smooth and planar discontinuity filled with an elastic material is validated through comparisons with analytical solutions. Next, the effect of the plastic deformation of the fill on the wave propagation is investigated. The model is extended to further study S-wave propagation across a filled discontinuity with rough interfaces, which is described using a sawtooth. The transmission coefficient defined by the energy is used to measure the wave attenuation. Finally, a parametric study is conducted to investigate the influences of the filled thickness, asperity angle, and incident amplitude on the transmission waves and transmission coefficients. The asperity angle and filled thickness together determine the transmitted waveform and transmission coefficient. The transmitted wave may be cut off when the incident wave amplitude exceeds a threshold value. The transmission coefficient decreases with a different trend with the incident wave amplitude increasing when the asperity angle varies. Compared with planar discontinuity, a filled discontinuity with rough interfaces is more sensitive to the amplitude of the incident wave. The causes of these phenomena are analyzed in detail. In addition, the deformation of the fill material is strongly related to the wave attenuation.