2. 河南义马煤业集团股份有限公司,河南义马 472300;
3. 吉林大学, 长春 130029;
4. 东北石油大学地球科学学院,辽宁大庆 163318
2. Henan Yima Coal Group Co.,LTD, Yima Henan 472300,China;
3. Jilin University, Changchun 130029,China;
4. Northeast Petroleum University, Daqing Liaoning 163318, China
在煤矿开采过程中常因遇到未探明的异常地质体使开采技术复杂化、开采成本增加、可开采煤炭储量降低,甚至导致灾难性事故频发,并使井下开采人员的生命安全受到威胁,造成巨大经济损失.因此在开采之前,探明更小的地质构造,对煤矿的开采计划、储量估计、危害评价、以及矿水管理等至关重要.而地面地震勘探具有不可克服的缺点,如高频成分随传播距离的增加而很快衰减、激发频率低、层间多次波干扰严重等,致使其分辨率精度不高,为解决这些问题,在实际工作中,普遍采用槽波探测技术[1-4].
槽波是在煤层中传播的一种导波(围限波或者管波),其传播距离远、能量强、波形特征易于识别,具有明显的频散特征.槽波地震勘探就是利用此特性来探测煤层内异常体[5-6].事实证明,这种方法不仅能够有效地确定异常体位置,而且探测范围也是其他物探方法无法比拟的.一般而言,反射方法探测的距离一般为煤层厚度的120 倍,而透射法探测距离能达到煤层厚度的300倍[7].
槽波探测技术首先应用于断层的识别,它可以提供断层的走向、类型,延伸程度与方向[8].由于槽波探测技术提供非直接信息,因此易于引起误差,误差大小与煤层的槽导性以及探测距离有直接关系.其探测结果的准确性、可靠程度随探测距离增加而减小.另外地质结构的复杂性也降低了对数据处理结果的解释精度.为了增加解释的可靠性,应做进一步的详细探测,如辅以其他物探方法,或增加探测精度[9-12].尽管如此,槽波对断层的识别率可达90%.
层析成像一般用易于识别的直达P 波或者折射波到时对地下结构进行速度结构成像[13-20].然而在槽波探测中,直达P 波反映的是围岩速度信息,但后续的槽波具有频散特性,与煤层的结构-煤厚变化、夹矸层、密度、速度、异常地质体等有直接关系.一般而言,在煤层内传播的槽波速度较低,煤层与顶底板围岩之间波阻抗差异大,在此条件下形成的槽波频散特征明显[21],长波长、低频成分波传播速度快,先于高频成分波到达,易于识别.
槽波波列上有一个特殊震相---埃里震相,其频率高、能量强,表现在群速度曲线上的极值点位置.槽波波列包含多个具有不同相速度的相位,拾取单一相速度到时比较困难,因此槽波走时成像,一般按一组波列包络极值点的走时进行反演计算[22-23].其中埃里震相就是这些波组包络中速度最低、能量最强、衰减较慢的一组波列.故埃里震相信息特别适用于走时反演,以了解煤层内的异常地质体信息.
对Love波和Rayleigh 波进行频散分析,可通过多次滤波技术[24-25]、时窗移动法[26]、Gabor矩阵分析[27]、小波变换[28]、以及基于多道处理技术的F-K 域变换[29]和拉东变换[30]等方法计算得到的.在频散分析曲线上能够容易地识别出埃里震相.
对每道数据进行滤波以增强埃里震相的能量.提取基阶模式曲线埃里震相到时进行反演计算,类似于P波到时层析成像反演.分别提取Love波和Rayleigh波埃里震相走时进行反演,这两类波反演得到煤层内速度扰动信息,结合围岩信息与显示的地质结构进行综合解释[31].
煤层厚度变化对Love 波频散曲线影响明显[2, 32-35],随着煤厚增大,Love波所含主要频段向低频移动,因此主频段槽波携带煤厚信息,应用层析成像方法获得速度分布与煤层厚度分布.
运用层析成像进行反演计算的算法有:反投影法(BG),求解逆矩阵的共轭梯度法(CG)与最小平方法(LSQR),代数迭代法(ART),瞬时迭代法(SIRT)等[36-41].我们选用瞬时迭代算法,主要是因为:(1)SIRT 适用于稀疏、不规则、信噪比低的实测数据[42];(2)易于加入先验信息,适用于弯曲射线;(3)方法稳健,但是该方法计算量大,耗时长[24].
本文通过对共检波点道集数据分析,得到围岩纵波波速、频率变化信息、巷道反射槽波以及具有完整槽波的路径信息;采用移动时窗法计算频散曲线并提取125Hz震相走时,应用于11061工作面透射层析成像,圈定了高速异常区;通过速度与厚度拟合曲线,获得工作面内煤层厚度分布图,为井下煤矿安全、高效开采设计提供了更详细、合理的物理参数.
2 层析成像方法原理射线是几何光学中导出的概念,它是用射线长度的仿射参数来描述的.在层析成像中要使用射线理论,首先需要对地震波的传播做几何光学近似,近似的条件为波速的空间变化在一个波长范围内很小,这实际上就是取波动方程的高频近似解.当用射线来描述地震波的传播时,由震源i到接收点j的理论走时Tij可用线积分表示为
(1) |
式中u(x,y,z)表示波的慢度(速度的倒数),V(x,y,z)表示波速,receive表示接收点位置,source表示震源位置,dl表示沿射线路径的线元.根据费马原理和斯涅尔定律,在各向同性、完全弹性的分层介质模型中,可以追踪出地震波射线路径.
在层析成像中,对介质进行网格离散化(图 1),积分公式转变为离散化形式
(2) |
其中tkcalc、dkj、sj、N分别表示第k条射线的走时,第k条射线经过第j面元的长度,第j面元的慢度,射线经过的网格剖分点数.
所有射线路径走时方程可以写成矩阵形式:
(3) |
其中T为M维时间向量,S为N维慢度向量,D为(M×N)矩阵算子,利用瞬时迭代重构法求解,得到速度分布图像.
瞬时迭代重构法求解过程:
(1) 用方程(2)对给定的初始模型计算走时;
(2) 计算射线经过每一个面元的慢度校正量;
(3) 将上述结果加上初始慢度产生新的慢度;
(4) 重新计算走时,并与实际走时结果进行比较,直到满足最小误差条件.
在槽波层析成像中,从频散曲线中拾取埃里相走时(速度)或某一频率震相走时(速度)信息,用于反演计算.
3 实例分析 3.1 施工设计与数据采集为验证上述理论,我们将槽波透射法应用到河南义马矿区11061工作面的实际探测中,以预测煤层厚度变化以及压力异常区(图 2).工作面宽度135m, 勘探面长度420 m, 煤层厚度变化1.5~8 m, 显示的煤层中没有夹矸层,揭露的断层区域用红色椭圆标注,排瓦斯钻孔位于黄色椭圆范围内.需要探测煤层厚度分布情况.观测系统布置如下:炮线安置于轨道巷(红色星号位置),检波器安置于皮带巷(绿色箭头指示位置).图中上下部分为皮带巷与轨道巷揭露的煤层厚度变化分布.
激发震源为炸药类型,药量200g, 炮距10 m, 共36炮,采用2分量检波器接收,检波点距为20m, 共22个检波点.
3.2 数据分析原始数据分析:将共炮集数据抽成共接收点集数据,并经过滤波处理,只有编号06(见图 3)、07、08、09检波器可以采集到完整的煤层波信息,由此可以推测此处的地质异常体比较小.从编号10 检波器Y分量数据分析得到P波速度,从X分量观测到槽波频率变化情况(见图 4).
频散曲线分析:根据模型分析对比[28, 43-45],可得到该工作面的一些参数,工作面巷道揭露的煤层厚度变化剧烈,厚度在1.5~8m 之间;1000m/s的速度值区域表示该区域煤层厚度大于或等于5 m 左右;速度大于2200m/s的区域表示该区域煤层厚度小于或等于2.0 m.该速度值与围岩和薄煤层的S波的速度相近(图 5).根据对理论频散曲线的分析,确定以125Hz槽波的理论到时进行层析成像反演计算.
图 6为从实际数据提取的频散曲线,这些数据来自不同的炮点检波点(S20-G05,S29-G06,S26-G14).从图中可确定埃里相速度及125 Hz槽波速度与走时.
作为一个实例,我们选择125Hz频率用于观测厚度约为2~5m 的煤层的厚度变化(图 7).从图 7所示的速度分布分析,红色部分(高速度区)表示此处有一个与工作面掘进方向近乎垂直、近乎西北至东南走向、宽约为150 m 的较大地质异常体.分析表明这些断层的延伸范围不会太大或者工作面内的断层落差与煤层厚度比小于0.5.
煤层厚度的变化对槽波频散特征影响十分明显,随着厚度的增加,勒夫波主频段降低,同时对应于同一频率波速降低,即同一频率,波速随厚度增加而减小.按照巷道内显示的煤层厚度与层析成像结果,采用4阶多项式拟合速度-厚度关系(图 8),并转换成煤层厚度分布图(图 9).从图 8 可看到从巷道内显示的煤层厚度与速度分布对应关系,速度范围800~2400m/s, 厚度变化范围1.5~8m, 变化趋势为随速度增加煤层厚度减薄.从图 9 可看到煤层厚度分布情况为中间薄、两端厚,从右至左的变化是,在右侧煤层厚度达6 m, 到工区中部,煤层变薄,在左侧煤层厚度又增加到5m 左右.
由于沿着巷道煤厚变化剧烈,进行精确的层析成像十分困难,成像使用的125Hz频率是一个折衷的办法,它不适用于煤厚较大的区域(煤厚大于5m),对煤厚小于2m 的区域也不完全适用,仅对2~5m之间变化的区域比较可信.
图 7中是以速度等值线分布的形式显示的层析成像结果,1000m/s的速度值指的是煤厚大于等于5m 的区域,5~8 m 或8 m 以上的区域则无法区分,只有选择小于100 Hz的频率才可能区分.煤厚小于等于2.5m 区域波速值高于2000m/s, 这个速度接近横波在围岩中的传播速度,在煤层较薄的区域波速不会明显地改变,针对较薄的煤层需要选择较高的频率才能获得较高的分辨率.成像结果里显示了一个宽达150 m 走向为NW-SE 的、与工作面走向垂直的高速异常区,这个异常区可能与工作面东边界的断层有关.由较高的速度变化可以看出东部边界非常明显,高速区域常常显示的是煤层结构变化较大的地方或反应该区域具有较高应力分布,这个断层的规模不是很大,在工作面内部落差只有煤厚的1/2左右;另外煤层顶板或底板的岩性可能会有变化,从勘探区域的末端煤厚又开始变厚.在实际回采中也验证了本研究成果的准确性.
当然,影响成像质量的因素很多,如煤厚的变化使槽波在其中的传播像经过透镜体一样,其传播路线发生弯曲,造成以直射线层析成像结果存在误差;另外,煤厚变化较小的区域也无法识别,这不仅与所选频率成分有关,亦与勘探选择的精度,比如炮检距及其坐标有关,因此,需要根据研究问题的特殊性选择施工方案和参数,以期达到更好的效果.
[1] | Evison F F. A coal seam as a guide for seismic energy. Nature , 1955, 176(4495): 1224-1225. DOI:10.1038/1761224a0 |
[2] | Krey T C. Channel waves as a tool of applied geophysics in coal mining. Geophysics , 1963, 28(5): 701-714. DOI:10.1190/1.1439258 |
[3] | 孙瑞霞, 张雯霁, 黄晓玲. 煤田地震勘探的新方法——槽坡地震勘探. 河北建筑科技学院学报 (自然科学版) , 1998, 15(3): 59–72. Sun R X, Zhang W J, Huang X L. A new method of seismic prospecting in coalfield. Journal of Hebei Institute of Architectural Science and Technology (in Chinese) , 1998, 15(3): 59-72. |
[4] | 刘天放, 潘冬明, 李德春, 等. 槽波地震勘探. 北京: 中国矿业大学出版社, 1994 . Liu T F, Pan D M, Li D C, et al. In-Seam Seismic Exploration (in Chinese). Beijing: China Mining University Press, 1994 . |
[5] | Buchanan D J. The propagation of attenuated SH channel waves. Geophys. Prosp. , 1978, 26(1): 16-8l. DOI:10.1111/gpr.1978.26.issue-1 |
[6] | Krey T C, Arnetzl H H, Knecht M. Theoretical and practical aspects of absorption in the application of in-seam seismic coal exploration. Geophysics , 1982, 47(12): 1645-1656. DOI:10.1190/1.1441314 |
[7] | Lagsse P E, Mason I M. Guided modes in coal seams and their application to underground seismic surveying.// Proc. IEEE Ultrasonic Symposium. Los Angeles, 1975: 64-67. |
[8] | Kerner C, Dressen L. The influence of dirt bands and faults on the propagation of Love seam waves. J. Geophys. , 1985, 57: 77-89. |
[9] | 卫金善, 张晋武. 综合勘探方法在成庄矿井地质构造探测中的应用. 中国煤田地质 , 2002, 14(4): 19–21. Wei J S, Zhang J W. Application of comprehensive prospecting method on search of geological structure in Chengzhuang mine. Coal Geology of China (in Chinese) , 2002, 14(4): 19-21. |
[10] | 程建远, 李淅龙, 张广忠, 等. 煤矿井下地震勘探技术应用现状与发展展望. 勘探地球物理进展 , 2009, 32(2): 76–82. Cheng J Y, Li X L, Zhang G Z, et al. Current status and outlook of seismic exploration applied underground in coal mine. Progress in Exploration Geophysics (in Chinese) , 2009, 32(2): 76-82. |
[11] | Luo X, Hatherly P, Duncan F M, et al. Microseismic monitoring of roof reaction to high wall mining //Aziz N, Indraratna B. Proceedings International Conference on Geo mechanics / Ground Control in Mining and Underground Construction. Australia: CSRIO, 1998: 1059-1065. |
[12] | Phllips W S, Rutledge J T, House L S, et al. Induced microearthquake patterns in hydrocarbon and geothermal reservoirs: six case studies. Pure and Applied Geophysics , 2002, 159(1-3): 345-369. |
[13] | Anderson D L. Global mapping of the uppermantle by surface wave tomography.// Midterm Report of the International Lithosphere Program in the Volume, The Composition, Structure and Dynamics of the Lithosphere-Asthenosphere System, 1987: 89-97. |
[14] | Chapman C H, Pratt R G. Traveltime tomography in anisotropic media – I. theory. Geophys. J. Int. , 1992, 109(1): 1-19. DOI:10.1111/gji.1992.109.issue-1 |
[15] | Lo T W, Inderwiesen P L. Fundamentals of seismic tomography. Society of Exploration Geophysicists , 1994: 1-178. |
[16] | Zhang J, Toks?z M N. Nonlinear refraction traveltime tomography. Geophysics , 1998, 63(5): 1726-1737. DOI:10.1190/1.1444468 |
[17] | Zhao D P, Hasegawa A. P wave tomographic imaging of the crust and upper mantle beneath the Japan islands. J. Geophys. Res. , 1993, 98(B3): 4333-4353. DOI:10.1029/92JB02295 |
[18] | 房立华, 吴建平, 吕作勇. 华北地区基于噪声的瑞利面波群速度层析成像. 地球物理学报 , 2009, 52(3): 663–671. Fang L H, Wu J P, Lü Z Y. Rayleigh wave group velocity tomography from ambient seismic noise in North China. Chinese J. Geophys. (in Chinese) , 2009, 52(3): 663-671. DOI:10.1002/cjg2.v52.3 |
[19] | 何正勤, 叶太兰, 丁志峰. 华北东北部的面波相速度层析成像研究. 地球物理学报 , 2009, 52(5): 1233–1242. He Z Q, Ye D L, Ding Z F. Surface wave tomography for the phase velocity in the northeastern part of North China. Chinese J. Geophys. (in Chinese) , 2009, 52(5): 1233-1242. |
[20] | 黄国娇, 白超英. 二维复杂层状介质中地震多波走时联合反演成像. 地球物理学报 , 2010, 53(12): 2972–2981. Huang G J, Bai C Y. Simultaneous inversion with multiple traveltimes within 2-D complex layered media. Chinese J. Geophys. (in Chinese) , 2010, 53(12): 2972-2981. |
[21] | 刘玉忠. 采区槽波探测的适应性分类. 煤田地质与勘探 , 1997, 25(4): 54–57. Liu Y Z. The adaptability classification of in seam wave surveying in mining districts. Coal Geology & Exploration (in Chinese) , 1997, 25(4): 54-57. |
[22] | 杨元海. 槽波速度谱相对透射系数图法. 煤田地质与勘探 , 1993, 21(4): 53–56. Yang Y H. The method for relative transmittance map of seam-wave velocity spectrum. Coal Geology & Exploration (in Chinese) , 1993, 21(4): 53-56. |
[23] | Wilson R G. In-seam seismic detection of geological structures. Geoexploration , 1987, 24(4-5): 295-300. DOI:10.1016/0016-7142(87)90003-2 |
[24] | Krajewski C, Dresen L, Gelbke C, et al. Iterative tomographic methods to locate seismic low-velocity anomalies: a model study. Geophysical Prospecting , 1989, 37(7): 717-751. DOI:10.1111/gpr.1989.37.issue-7 |
[25] | Buchanan D J, Jackson P J. Dispersion relation extraction by multi-trace analysis. Bull. Seis. Soc. Am. , 1983, 73(2): 391-404. |
[26] | Cox K B, Mason I M. Velocity analysis of the SH-channel wave in the Schwalbach seam at Ensdorf Colliery. Geophysical Prospecting , 1988, 36(3): 298-317. DOI:10.1111/gpr.1988.36.issue-3 |
[27] | R?der D, Schott W, Dresen L, et al. Calculation of dispersion curves and amplitude-depth distributions of love channel waves in horizontally-layered media. Geophysical Prospecting , 1985, 33(6): 800-816. DOI:10.1111/gpr.1985.33.issue-6 |
[28] | 杨真, 冯涛, WangS G. 0.9 m薄煤层SH型槽波频散特征及波形模式. 地球物理学报 , 2010, 53(2): 442–449. Yang Z, Feng T, Wang S G. Dispersion characteristics and wave shape mode of SH channel wave in a 0. Chinese J. Geophys. (in Chinese) , 2010, 53(2): 442-449. |
[29] | 张碧星, 鲁来玉. 用频率-波数法分析瑞利波频散曲线. 工程地球物理学报 , 2005, 2(4): 245–255. Zhang B X, Lu L Y. Investigation on the dispersion curves of rayleigh wave by Frequency-Wavenumber analysis method. Chinese Journal of Engineering Geophysics (in Chinese) , 2005, 2(4): 245-255. |
[30] | 潘冬明, 胡明顺, 崔若飞, 等. 基于拉东变换的瑞雷面波频散分析与应用. 地球物理学报 , 2010, 53(11): 2760–2766. Pan D M, Hu M S, Cui R F, et al. Dispersion analysis of Rayleigh surface waves and application based on Radon transform. Chinese J. Geophys. (in Chinese) , 2010, 53(11): 2760-2766. |
[31] | Bodoky T, Hermann L, Dianiska L. Processing of the in-seam seismic transmission measurements // 47th Annual EAEG Meeting in Budapest. Hungary, 1985. |
[32] | 程久龙. Love型槽波理论地震图的计算. 山东矿业学院学报 , 1994, 13(4): 349–353. Cheng J L. Calculation of theoretical seismograms of Love channel waves. Journal of Shandong Mining Institute (in Chinese) , 1994, 13(4): 349-353. |
[33] | 潘冬明, 刘天放. 拉夫型槽波的有限差分合成. 中国煤田地质 , 1990, 2(2): 54–59. Pan D M, Liu T F. Love channel wave finite difference synthesis. Coal Geology of China (in Chinese) , 1990, 2(2): 54-59. |
[34] | 杨文强. 槽波地震勘探的数学模型研究. 地质与勘探 , 2001, 37(3): 58–60. Yang W Q. Modeling research of channel wave seismic exploration. Geology and Prospecting (in Chinese) , 2001, 37(3): 58-60. |
[35] | 刘天放, 程久龙, 潘冬明, 等. 槽波的吸收衰减. 煤炭学报 , 1993, 18(5): 83–86. Liu T F, Cheng J L, Pan D M, et al. Channel wave attenuation by medium absorption. Journal of China Coal Society (in Chinese) , 1993, 18(5): 83-86. |
[36] | Herman G T. Image Reconstruction from Projections: The Fundamentals of Computerized Tomography. New York: Academic Press, 1980 . |
[37] | Hestenes M R, Stiefel E. Methods of conjugate gradients for solving linear systems. Nat. Bur. Standards J. Res. , 1952, 49(6): 409-436. DOI:10.6028/jres.049.044 |
[38] | 崔岩, 王彦飞, 杨长春. 带先验知识的波阻抗反演正则化方法研究. 地球物理学报 , 2009, 52(8): 2135–214l. Cui Y, Wang Y F, Yang C C. Regularizing method with a priori knowledge for seismic impedance inversion. Chinese J. Geophys. (in Chinese) , 2009, 52(8): 2135-214l. |
[39] | Paige C C, Saunders M A. LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Trans. Math. Softw. , 1982. |
[40] | Gordon R, Bender R, Herman G T. Algebraic Reconstruction Techniques (ART) for three-dimensional electron microscopy and x-ray photography. J. Theor. Biol. , 1970, 29(3): 471-481. DOI:10.1016/0022-5193(70)90109-8 |
[41] | Gilbert P E C. Iterative methods for the three-dimensional reconstruction of an object from projections. J. Theor. Biol. , 1972, 36(1): 105-117. DOI:10.1016/0022-5193(72)90180-4 |
[42] | Miranda F. Geophysical tomographic reconstruction of travel-time and amplitude anomalies. University College of Swansea, Department of Earth Sciences , 1989. |
[43] | Hu Y R, McMechan G A. Imaging mining hazards within coalbeds using prestack wave equation migration of in-seam seismic survey data: A feasibility study with synthetic data. Journal of Applied Geophysics , 2007, 63(1): 24-34. DOI:10.1016/j.jappgeo.2007.03.002 |
[44] | Xu Y Z, Cui R F, Huang W C, et al. Reflectivity forward modeling and a CSSI method seismic inversion study of igneous intrusive area, coked area, and gas-enriched area located within a coal seam. Mining Science and Technology , 2009, 19(4): 457-462. |
[45] | Ge M, Wang H, Hardy H R Jr, et al. Void detection at an anthracite mine using an in-seam seismic method. International Journal of Coal Geology , 2008, 73(3-4): 201-212. DOI:10.1016/j.coal.2007.05.004 |