一种天然源煤槽地震波勘探方法

测量装置的制造及其应用技术1.本发明涉及一种煤层地震波勘探方法，具体是一种天然源煤槽地震波勘探方法。背景技术：2.在现有的矿井地球物理勘探方法中，地震槽波勘探是具有发展潜力和应用前景的物探手段方法之一。它具有探测距离远、分辨精度高、抗电干扰能力强、波形特征较易识别以及最终成果直观的优点。其震源是关键部分，直接影响地震信号质量。槽波地震勘探的常用震源有炸药震源、锤击震源、雷管震源、夯击震源、皮带机震源、掘进机震源、采煤机震源等。目前，矿井震源以炸药震源为主。3.由于煤炭的开采过程中会伴随着煤与瓦斯突出、瓦斯爆炸、透水、顶板垮落等多种灾害事故。其中煤与瓦斯突出灾害一直是世界矿井领域的巨大难题之一，也是煤矿生产安全的最大威胁之一。随着采掘深度的不断增加，煤岩应力条件更加复杂，高瓦斯突出问题会愈加凸显出来。这些直接导致以炸药震源为主的槽波勘探在高瓦斯突出矿井中并不适用。锤击、夯击等震源相对操作简单，但其能量弱，可探测距离较短；掘进机、采煤机震源逐渐成为研究热点，但尚处于试验研究阶段，不适用于槽波勘探。因此如何找到一种新的非炸药震源槽波勘探方法，使其能在保证对煤层进行精准探测的前提下，具有观测系统布设方便使其每次探测效率高，且探测范围较广，是本行业的研究方向之一。技术实现要素：4.针对上述现有技术存在的问题，本发明提供一种天然源煤槽地震波勘探方法，利用天然源槽波进行地质探测，使其能在保证对煤层进行精准探测的前提下，具有观测系统布设方便使其每次探测效率高，且探测范围较广。5.为了实现上述目的，本发明采用的技术方案是：一种天然源煤槽地震波勘探方法，具体步骤为：6.a、观测系统布设及数据采集：在巷道靠近工作面的侧壁中部沿巷道走向呈一排布设多个检波器，且各个检波器之间的间距相等，将各个检波器均与地震监测仪连接组成观测系统；观测系统的各个检波器在一段时间内持续进行天然源槽波数据的采集，并将采集的数据反馈给地震监测仪，完成一次数据采集过程；7.b、数据处理：采用扩展空间自相关法对步骤a采集的天然源槽波数据进行处理，获得天然源槽波数据的频散谱，沿着频散谱的能量脊能提取频散点，最终根据提取各个频散点得到天然源槽波数据的频散曲线；8.c、频散曲线的解释：采用波长深度转换系数法对步骤b获得的频散曲线进行解释，从而能得出巷道周围煤层及岩层的纵横波速度剖面图，最后根据煤层及岩层的纵横波速度剖面图确定周围煤层及岩层的地质信息。9.本发明利用天然源槽波探查巷道周围煤层及围岩的地质信息的原理是：目前行业公知，煤层中存在着频率低、发生时间长的天然震源，无论是地表各类人工或机械活动，还是下部工作面中的采煤工作区等机械运作、人工活动、岩层断裂都会因振动产生地震波。本技术的发明人研究发现，当这些地震波在地层传播过程中若以一定角度由巷道顶底板的岩土体入射到煤层内时，如图1所示，多种地震波会在煤层内发生相互干涉并形成天然源槽波，天然源槽波在巷道周围煤层及围岩内传播后进入巷道内，并且发明人进一步研究发现，利用天然源槽波进行地质勘探相较于人工锤击、夯击等震源，其具有探查范围大、能量较强、可探测距离远、操作简单、不影响生产和可持续监测的显著优势；并且由于煤巷地质条件沿巷道走向变化不大的情况下，天然源槽波的频散特征与煤层厚度、煤层纵横波速度以及围岩纵横波速度相关。因此如能提取出天然源槽波的频散特征曲线，并确定巷道周围煤层及围岩距离巷道的不同深度位置中不同波长的煤层纵横波速度以及围岩纵横波速度情况，则能得出巷道周围煤层及围岩不同深度各自的纵横波速度剖面图，最后根据煤层及围岩各自的纵横波速度剖面图确定周围煤层及围岩的地质信息。10.进一步，所述步骤b具体为：11.扩展空间自相关法是保持频率不变，改变圆周半径，计算自相关系数并与贝塞尔函数拟合，求得自相关系数随距离的变化关系.不同半径下频率为f的空间自相关系数s定义为：[0012][0013][0014]上式中，r0n表示所改变的半径，γ0为初始半径，对应第一类零阶贝塞尔函数j0，n代表所改变的半径数，e表示拟合误差.通过拟合计算的空间自相关系数s0n(f,r0n)能导出槽波相速度；[0015]①选择一个提取频率，以线性排列的检波器中点为圆心点o，各个检波器与圆心点o之间的距离为半径，从而对不同半径的各个检波器进行数据提取；[0016]②对提取到的数据进行滤波处理；[0017]③分别计算不同半径圆心点o与圆周上第i个检波器采集微动信号的自功率谱s(r,ω)、s(r,ω,θi)和互功率谱soi(r,ω,θi)，进而得到该频率下的spac系数ρ(r,ω)；[0018]④根据上述的公式将spac系数与贝塞尔函数进行拟合，求得spac系数随着距离的变化关系；再根据最小误差原则将其与第一类零阶贝塞尔函数j0进行拟合，求得该频率下的瑞雷波的相速度c(ω0)；再选择一个提取频率并重复步骤①至④，从而能依次求出每一频率下瑞雷波的相速度，进而能计算获取频散谱，沿着频散谱的能量脊能提取频散点，最终根据提取各个频散点得到天然源槽波数据的频散曲线。[0019]进一步，所述步骤a中检波器的数量为12个，各个检波器之间的间距为10m，每次数据采集的时长为120s，数据的采样间隔为1ms。[0020]进一步，所述步骤c具体为：[0021]波长深度转换系数法中基于经典的槽波理论分析，计算出槽波波长与其勘探深度的转换关系为h＝k×λ；其中k为波长深度转换系数，h为勘探深度，λ为槽波波长；由于天然源槽波的频散曲线特征与煤层厚度、煤层纵横波速度以及围岩纵横波速度相关，为了获取不同地层深度的k值，在该巷道进行地层探查之前，先在该矿区采用多种勘探方式进行地质探测获取该矿区巷道周围煤层及围岩距离巷道的不同深度位置中不同波长的煤层纵横波速度以及围岩纵横波速度情况，从而确定巷道周围煤层及围岩不同深度位置的深度转化系数k；最后采用波长深度转换系数法对天然源槽波的频散曲线进行解释，从而得出巷道周围煤层及围岩不同深度各自的纵横波速度剖面图，最后煤层及围岩各自的纵横波速度剖面图确定周围煤层及围岩的地质信息。[0022]与现有技术相比，本发明先在巷道内布设观测系统，然后利用观测系统进行天然源槽波数据的采集，由于发明人研究发现天然源槽波的频散特征与煤层厚度、煤层纵横波速度以及围岩纵横波速度相关。因此如能提取出天然源槽波的频散特征曲线，并确定巷道周围煤层及围岩距离巷道的不同深度位置中不同波长的煤层纵横波速度以及围岩纵横波速度情况，则能得出巷道周围煤层及围岩不同深度各自的纵横波速度剖面图，最后根据煤层及围岩各自的纵横波速度剖面图确定周围煤层及围岩的地质信息。接着本发明采用扩展空间自相关法提取出天然源槽波的频散特征曲线；并采用波长深度转换系数法对获得的频散曲线进行解释，从而能得出巷道周围煤层及围岩不同深度各自的纵横波速度剖面图，最后根据煤层及围岩不同深度各自的纵横波速度剖面图确定周围煤层及围岩的地质信息。因此本发明利用天然源槽波进行地质探测，使其能在保证对煤层进行精准探测的前提下，具有观测系统布设方便使其每次探测效率高，且探测范围较广。附图说明[0023]图1是本发明中天然源槽波形成的示意图；[0024]图2是本发明观测系统的布置平面图；[0025]图3是本发明中观测系统接收不同方向天然源槽波的示意图；[0026]图4是本发明采集的天然源槽波的数据图[0027]图5是本发明对图4采用扩展空间自相关法获得的频散谱结果图。具体实施方式[0028]下面将对本发明作进一步说明。[0029]如图2所示，本发明的具体步骤为：[0030]a、观测系统布设及数据采集：在巷道靠近工作面的侧壁中部沿巷道走向呈一排布设12个检波器，且各个检波器之间的间距均为10m，将各个检波器均与地震监测仪连接组成观测系统；如图3所示，观测系统的各个检波器在120s时间内持续进行天然源槽波数据的采集，数据的采样间隔为1ms，并将采集的数据反馈给地震监测仪，完成一次数据采集过程；[0031]b、数据处理：如图4所示，采用扩展空间自相关法对步骤a采集的天然源槽波数据进行处理，获得天然源槽波数据的频散谱，如图5所示，沿着频散谱的能量脊能提取频散点，最终根据提取各个频散点得到天然源槽波数据的频散曲线，具体为：[0032]扩展空间自相关法是保持频率不变，改变圆周半径，计算自相关系数并与贝塞尔函数拟合，求得自相关系数随距离的变化关系.不同半径下频率为f的空间自相关系数s定义为：[0033][0034][0035]上式中，r0n表示所改变的半径，γ0为初始半径，对应第一类零阶贝塞尔函数j0，n代表所改变的半径数，e表示拟合误差.通过拟合计算的空间自相关系数s0n(f,r0n)能导出槽波相速度；[0036]①选择一个提取频率，以线性排列的检波器中点为圆心点o，各个检波器与圆心点o之间的距离为半径，从而对不同半径的各个检波器进行数据提取；[0037]②对提取到的数据进行滤波处理；[0038]③分别计算不同半径圆心点o与圆周上第i个检波器采集微动信号的自功率谱s(r,ω)、s(r,ω,θi)和互功率谱soi(r,ω,θi)，进而得到该频率下的spac系数ρ(r,ω)；[0039]④根据上述的公式将spac系数与贝塞尔函数进行拟合，求得spac系数随着距离的变化关系；再根据最小误差原则将其与第一类零阶贝塞尔函数j0进行拟合，求得该频率下的瑞雷波的相速度c(ω0)；再选择一个提取频率并重复步骤①至④，从而能依次求出每一频率下瑞雷波的相速度，进而能计算获取频散谱，沿着频散谱的能量脊能提取频散点，最终根据提取各个频散点得到天然源槽波数据的频散曲线。[0040]c、频散曲线的解释：采用波长深度转换系数法对步骤b获得的频散曲线进行解释，从而能得出巷道周围煤层及岩层的纵横波速度剖面图，最后根据煤层及岩层的纵横波速度剖面图确定周围煤层及岩层的地质信息，具体为：[0041]波长深度转换系数法中基于经典的槽波理论分析，计算出槽波波长与其勘探深度的转换关系为h＝k×λ；其中k为波长深度转换系数，h为勘探深度，λ为槽波波长；由于天然源槽波的频散曲线特征与煤层厚度、煤层纵横波速度以及围岩纵横波速度相关，为了获取不同地层深度的k值，在该巷道进行地层探查之前，先在该矿区采用多种勘探方式进行地质探测获取该矿区巷道周围煤层及围岩距离巷道的不同深度位置中不同波长的煤层纵横波速度以及围岩纵横波速度情况，从而确定巷道周围煤层及围岩不同深度位置的深度转化系数k；最后采用波长深度转换系数法对天然源槽波的频散曲线进行解释，从而得出巷道周围煤层及围岩不同深度各自的纵横波速度剖面图，最后煤层及围岩各自的纵横波速度剖面图确定周围煤层及围岩的地质信息。[0042]以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。









图片声明：本站部分配图来自人工智能系统AI生成,觅知网授权图片,PxHere摄影无版权图库。本站只作为美观性配图使用,无任何非法侵犯第三方意图,一切解释权归图片著作权方,本站不承担任何责任。如有恶意碰瓷者,必当奉陪到底严惩不贷!




内容声明：本文中引用的各种信息及资料（包括但不限于文字、数据、图表及超链接等）均来源于该信息及资料的相关主体（包括但不限于公司、媒体、协会等机构）的官方网站或公开发表的信息。部分内容参考包括:(百度百科,百度知道,头条百科,中国民法典,刑法,牛津词典,新华词典,汉语词典,国家院校,科普平台)等数据,内容仅供参考使用,不准确地方联系删除处理！本站为非盈利性质站点,发布内容不收取任何费用也不接任何广告!




免责声明：我们致力于保护作者版权，注重分享，被刊用文章因无法核实真实出处，未能及时与作者取得联系，或有版权异议的，请联系管理员，我们会立即处理，本文部分文字与图片资源来自于网络，部分文章是来自自研大数据AI进行生成,内容摘自(百度百科,百度知道,头条百科,中国民法典,刑法,牛津词典,新华词典,汉语词典,国家院校,科普平台)等数据,内容仅供学习参考,不准确地方联系删除处理!的,若有来源标注错误或侵犯了您的合法权益，请立即通知我们，情况属实，我们会第一时间予以删除，并同时向您表示歉意,谢谢!