随着国民经济的持续快速发展,国家“西部大开发”战略的逐步实施,对煤炭资源丰富的西部(新疆、宁夏、内蒙)地区进行开发利用,振兴西部经济十分有益.为了科学、更好地开发利用西部地区的煤炭资源,首先进行三维反射波地震勘探,以其查明煤层赋存范围、构造形态和储量(曹志勇等,2012;李维新等, 2009;孙文涛等,1997;赵殿栋,2009;张向林等,2007).
煤田地震勘探多年来取得很多进展,三维地震勘探及相应的3D资料处理得到飞速的发展 (崔若飞等,2008;马在田,1988;麻三怀等,2008;田春志等,2002),这是油气及矿产资源精细勘查的市场需求.虽然三维反射波地震勘探已成为煤炭开发前的常规性勘查技术工作(程建远等,2009;程建远等,2010;耿丽娟,2005.王伟等,2012),可是这些地区大多地形、地貌复杂多样,煤层构造较为复杂,干旱少雨,同时高山深谷相互切割,从而造成地震勘探激发与接收条件的差异较大,影响煤炭采区三维地震勘探方法的开展.为此,针对西部地区独特的地形、地貌与结构特点,对煤层反射波的激发、接收、数据采集条件和三维地震资料精细处理和解释进行探索性(李彦鹏等,2012;王建民等,2007;王赟等,2013;膝佃波等,2005),分析研究西部煤炭采区三维地震勘探方法技术及其效果是非常必要的.
本文针对新疆库车阿艾煤矿采区的地质地球物理概况以及该区三维地震勘探的难点,论述了高密度三维地震勘探实践中的地球物理勘探方法技术及勘探结果.
阿艾勘探区位于新疆维吾尔自治区的阿克苏地区,该区距库车县城约80余公里,行政区隶属阿克苏地区库车县阿艾乡管辖,217国道是贯穿天山南北的大通道,沿国道217 (独—库公路段)北行可抵达勘探区.全区二维、三维地震勘探面积为15.26 km2,其中三维地震勘探有效控制面积8.20 km2.
(1)查明测区内下1、下5、下10煤层的赋存形态(深度误差 不能大于1.0%) 及煤层厚度变化情况; 查明测区内主要煤 层露头位置,平面位置误差不大于30 m.
(2)查明测区内落差大于或等于5 m的断层,断层平面位置误 差不大于15 m,并解释3~5 m的断点;查明F2及F8 断层的赋存状态.
(3)查明测区内火烧区的各煤层下边界及火烧区内对各煤层的破坏程度.
(4) 利用瞬变电磁勘探手段探查火烧区、F2及F8断层的赋水性.
(1)地层:区内地层自下而上有:三叠系上统黄山街组(T3h)、侏罗系下统塔里奇克组(J1t)、阿合组(J1a),侏罗系中统克孜努尔组(J2k)和少量的第四系(Q4)
(2)构造:区内断裂不发育,小断层对煤层影响不大.测区南部有一条近东西走向的逆断层,倾角45°左右,倾向南,发育在测区西界外.
(3)煤层:全区可采煤层7层,煤层平均总厚29.38 m,所有煤层在地表均未见出露,煤层在地表均全部火烧.
其中全区主要可采煤层3层,均为厚煤层,埋藏深度为100~400 m.
测区位于塔里木盆地与天山山脉结合的过渡地带,属中低山区,表浅层地震地质条件的复杂多样,砂砾岩层覆盖全域;沿沟的断崖绝壁和基岩裸露的秃山石漠和地形切割细碎的地质条件,给地震资料采集带来很大困难.主要表现为:
(1)测区属低山区,沟壑纵横,道路稀少,设备搬运极为不便的特点,所以采用加拿大ARIES地震仪器,它具有灵活机动特点,适合野外适用.
(2)针对地表大部分基岩出露、冲沟发育,高差变化较大,炮点和检波器点的偏移无法避免的问题.采用山地钻成孔,保证炸药下到规定深度,并封孔激发,炮点移动要实测坐标;
(3)针对基岩裸露地段,检波器如何埋置问题.采用石膏固化法,把检波器粘在基岩上的方法解决,保证检波器埋置在设计位置;
(4)针对要查明落差5 m的断层及勘探区内火烧区范围,采集资料频宽要宽,高频成分要丰富的问题.采用60 Hz检波器接收及24次叠加次数的观测系统,保证纵横向叠加次数较为均匀,从而提高资料的分辨率和信噪比;
(5)针对煤层层数多,且埋深差异大问题.采用设计了不同的观测系统,保证了浅、中、深部目的层的叠加次数.另外,资料处理时制定了不同的处理流程,以保证各个目的层的信噪比和分辨率;
(6)针对区内火烧对煤层的破坏情况,利用高分辨三维地震、磁法相结合手段、多种参数,综合确定.
测区复杂地形条件下的地震资料采集,要有特定的施工 设备、采集参数、与之相适应的观测系统和行之有效的施工方法,才能确保较高质量的原始资料.为此在正式采用前我们进行了激发方式、采集方式的预研究.
对区内多处进行激发井深和药量点试验,井深分别为3 m、4 m、5 m、7 m、8 m、10 m、12 m共7种情况,药量分别为1.5 kg、2.0 kg、2.5 kg、3.0 kg共4种情况,考虑测区基本为 基岩孔及目的层埋深差异,试验结果井深为6 m,药量分别为1 kg(基岩裸露区),2 kg(第四系坡积物堆积区).其次进行了检波器类型及组合试验,试验结果为了在该区获得高质量原始记录,需要采用CDJ-60型高频检波器,和3串2并堆放.再次,对观测系统进行试验,采用采用加拿大ARIES地震仪器,束状8线4炮制中点激发,获得如下观测系统主要参数.
接收道数:8×72=576道(深部)、8×60=480道(中部)、8×36=288道(浅部);炮点网格:20 m×60 m(深部)、20 m×50 m(中部)、20 m×30 m(浅部);检波点网格:10 m×20 m;CDP网格:5 m×10 m;叠加次数: 6×4=24次.
上述观测系统可使记录获得较好的煤系地层反射波.
针对前述该区的资料采集条件分析,我们对资料处理流程中静校正、叠后去噪、偏移提频处理给予了重点关注(安西峰等,2010;董春晖等,2010;韩晓丽等,2008;兰海强等,2012;张正峰等,2013).
如图 1,所示,测区地形起伏较大,因此静校正工作是资料处理的重中之重,结合实际情况,采用绿山初至折射静校正的逐炮拾取阶段,保证所拾取的初至折射波来自于在全区较能连续追踪的同一层,建立精确的近地表模型.图 2 静校正前后的单炮记录,对比两个结果可以看出校正效果明显的.
为了得到准确可靠的偏移结果,提高叠加剖面的信噪比,去处随机噪音干扰是十分必要的,这样可以保证目地层在没有干扰下准确归位.对比一下可以看出目的层更加明显.如图3所示.
针对测区煤层多、层间距小的特点,为进一步提高分辨率,对剖面进行了提高频率处理,使有效波得到明显的提高,这样可以提高近距离煤层群分辨率间的距离(图4).
通过对新疆库车阿哎三维地震勘探资料的处理及对剖面的解释,对煤层、断层发育状况和主要煤层的分布形态及火烧区给予了解释,采用水平或平面图的方式显示具不连续性的断层、陷落柱、河道、地层边界等地质特征,来检查三维解释的正确性.如图5、图6、图7所示.
图5为下5煤层鸟瞰图.下5煤层较为稳定,能量强、连续性好,除部分地段可能受火烧影响外,基本全区赋存.
测区总体为一向南倾伏的弧状单斜构造,地层走向近东西为80°~280°,倾角为8°~20°,一般为10°左右.主采煤层下1、下5、下10赋存条件好,煤层顶底板岩性、岩相组合特征清楚,密度差异明显,煤层埋深较适中,能形成能量较强的反射波T1、T5、T10.另外,由于上部煤层的屏蔽作用,使得下部及中间薄煤层的反射波变弱,难于追踪识别.
测区北部由于煤层自燃后,烧变岩石较为破碎,裂隙发育,失去了形成层状地震反射波的条件. 煤层火烧区在时间剖面上的反映是反射波品质变差或突然消失,反射波变的杂乱无章. 如图6所示.
利用水平切片检查断层的纵、横向展布趋势,水平切片上同相轴的强度反映了反射波的强度,其宽度与地层的倾角大小有关.由于水平切片对断层具有较高的分辨能力,因此,充分利用它来识别断层,同时对断层组合的合理性及断层的延展情况进行检查.如图7所示.
在基岩裸露复杂地表条件区,利用“专用山地钻机成孔、高叠加次数、小炮检距、高频组合接收”采集技术,可获得高分辨率的原始资料,经过精细处理分析和综合对比解释研究,地质成果良好.实践证明,在新疆基岩裸露地区开展煤田三维地震勘探切实可行,我国西部煤炭勘探程度非常低,煤田地震勘探前景广阔.
[1] | An X F, Liu P, Zhang X. 2010.Surface consistent deconvolution effect and application in seismic exploration [J]. progress in geophysics.2010,25 (6): 2125-2129 |
[2] | Cheng J Y, Li X L, Zhang G Z, et al. 2009.Current status and outlook of seismic exploration applied underground in coal mine[J]. Progress in Exploration Geophysics (in Chinese), 32(2): 96-100. |
[3] | Cheng J Y, Zhao W, Cao D T, et al. 2010.Pondering on coalmine winning district 3d seismic correlation of prospecting and mining information effect analysis[J]. Coal Geology of China (in Chinese), 22(8): 67-72,82. |
[4] | Cao Z Y, Wang W W, Wang Y.2012.Numerial simulation and imaging of scattered wave of sunk pillar in coal seam.Chinese J.Geophys.(inChinese),(5):1749-1756,doi10.6038/j.issn.0001-5733. 05.032. |
[5] | Cui R F, Sun X K, Cui D W. 2008.Seismic inversion-new development in coal seismic survey[J]. Coal Geology of China (in Chinese), 20(6): 49-52. |
[6] | Dong chun-Hui, Zhang Jian-Feng.2010.Prestack time migration incorporated with redatuming.Chinese J. Geophys(in Chinese),53(10):2345-2441, doi:10.3969/j.issn.0001-5733.10.017. |
[7] | Geng Li juan.2005.3D seismic wave survey for coal layer in Hami desert area of Xinjiang[J]. Progress in Geophysics,20(2):393-398. |
[8] | Han Xiao-li,Yang Chang-chun,Ma San-huai,Qin Hong-guo .2008.Static of tomographic inversion by first breaks in complex areas [J]. Progress in Geophysics,23(2):475-483. |
[9] | Li W X,Wang H, Yao Zh,X, L Y K, Chang X. 2009.Constraint inversion of S-wave velocity and fluid replacement based on [J]. Progress in Geophysics. 52(3): 785-791. |
[10] | Li Yan-Peng,Ma Zai-Tian,Sun Peng-Yuan,Yang Hai-Shen.2012.Converted-wave static correction method for thick weathering area .ChineseJ.Geophys.(inChinese),55(2):614-621, doi10.3969/j.issn.0001-5733.02.024. |
[11] | Lan Hai-Qiang,Zhang Zhi,Xu Tao,Bai Zhi-Ming.2012.Effects due to the anisotropic stretching of the surface-fitting grid on the traveltime computation for irregular surface by the coordinate transforming method.ChineseJ.Geophys.(inChinese),(10):3355-3369,doi:10.6038/j.issn.0001-5733.2012.10.018. |
[12] | Ma Zai-Tian.1988.On the problems of 3-D seismic data processing .Chinese J. Geophys. (inChinese), 31(01): 99-107. |
[13] | Ma S H, Yang C C, Sun F L, Huang W C.2008. Kirchhoff prestack time migration technology applied in seismic data processing in complex structures [J]. progress in geophysics,23(3):754-760. |
[14] | Sun Wen-Tao ,Fang Zheng. 1997.Review and prospect of coal Geophysical technology in China.Chinese J.Geophys. (in Chinese), 40(S1):362-368. |
[15] | Tian C Z, Liu H. 2002.Study on seismic wave propagation for deep layer seismic exploration [J]. Progress in Geophysics, (in Chinese), 17(4):564-574. |
[16] | Wang Jian-Min,Chen Shu-Min,Su Mao-XinWang Yuan-Bo,Wang Li-Na,Guan Xin,Tian Gang,Shi Zhan-Jie.2007.A study of the near surface high-frequency compensation technology in 3-D seismic exploration .Chinese J. Geophys. (in Chinese), 50(6):1837-1843. |
[17] | Wang Yun,Zhang Yu-Gui,Xu Xiao-Kai.2013.Relationship between the maximum vitrinite reflectance and the elastic parameters of coal: A lab ultrasonic measurement of 6 metamorphic kinds of coals. Chinese Journal Geophysics, 56 (6): 2116-2122, doi :10. 6038/cjg20130631. |
[18] | Wang Wei,Gao Xing,Li Song-Ying,Yue Yong,Hu Guo-Ze,Li Yu-Ying.2012.Channelwavetomography method and itsapplicationincoalmineexploration :AnexamplefromHenanYimaMiningArea.Chinese J. Geophysics (in Chinese),55(3):1054-1062, doi:106038/j.issn.0001-5733.2012.03036. |
[19] | Xi T B, Wang P C, Wang Y, et al. 2005.Using prestack KirchhOff integral migration recognize small faults and low amplitude structure [J]. Progress in Exploration Geophysics,20(4):1035-1038. |
[20] | Zhao D D. 2009.Review and prospect on high-precision seismic exploration technique. Geophysical Prospecting for Petroleum (in Chinese), 48(5): 425-435. |
[21] | Zhang X L, Liu X R, Li J, et al.2007. Progress in oil productive monitoring technology in China[J]. Progress in Geophysics (in Chinese), 22(4): 1360-1363. |
[22] | Zhang Zheng-Feng,Guo Jian,Han Li-Guo,Xu Shan-Hui,Gong Xiang-Bo.2013.Transformation coordination converted wave static correction based on wave field continuation. Chinese Journal Geophysics, 56 (5): 1676-1685 ,doi :10. 6038/cjg20130525. |
[23] | 安西峰,刘平,张旭. 2010.地表一致性反褶积在地震勘探中的应用及效果[J]. 地球物理学进展.25(6):2125-2129. |
[24] | 程建远, 李锡龙, 张广忠,等. 2009.煤矿井下地震勘探技术应用现状与发展展望. 勘探地球物理进展, 32(2): 96-100. |
[25] | 程建远, 赵伟, 曹丁涛,等. 2010.煤矿采区三维地震探采对比效果的分析与思考[J]. 中国煤炭地质, 22(8): 67-72, 82. |
[26] | 曹志勇,王伟,王赟. 2012.煤层陷落柱散射波数值模拟与成像. 地球物理学报, 55(5):1749-1756,doi:10.6038/j.issn.0001-5733.2012.05.032 |
[27] | 崔若飞, 孙学凯, 崔大尉. 2008.地震反演-煤田地震勘探的新进展[J]. 中国煤田地质, 20(6): 49-52. |
[28] | 董春晖,张剑锋. 2010.结合基准面重建的叠前时间偏移方法[J]. 地球物理学报, 53(10):2435-2441,doi:10.3969/j.issn.0001-5733.2010.10.017. |
[29] | 耿丽娟. 2005.新疆哈密戈壁滩区煤层三维地震勘探实践 地球物理学进展,20(2):393-398. |
[30] | 韩晓丽,杨长春,麻三怀,等. 2008.复杂山区初至波层析反演静校正[J]. 地球物理学进展. 23(2):475-483 |
[31] | 李维新,王红,姚振兴,刘伊克,常旭. 2009.基于约束条件横波速度反演和流体替代[J]. 地球物理学进展, 52(3): 785-791. |
[32] | 李彦鹏,马在田,孙鹏远,杨海申.2012.厚风化层覆盖区转换波静校正方法. 地球物理学报,55(2):614-621,doi:10.3969/j.issn.0001-5733.2012.02.024. |
兰海强,张智,徐涛,白志明,2012.贴体网格各向异性对坐标变换法求解起伏地表下地震初至波走时的影响, 地球物理学报,(10):3355-3369,doi:10.6038/j.issn.0001-5733.2012.10.018. | |
[34] | 马在田. 1988.三维地震勘探数据处理的问题及其解决方法. 地球物理学报, 31(01): 99-107. |
[35] | 麻三怀,杨长春,孙福利,黄伟传. 2008.克希霍夫叠前时间偏移技术在复杂构造带地震资料处理中的应用[J].地球物理学进展. 23(3):754-760. |
[36] | 孙文涛,方正. 1997.我国煤田物探技术的回顾与展望. 地球物理学报, 40(S1):362-368. |
[37] | 田春志,刘洪.2002.深层地震勘探的地震波传播理论研究前景[J]. 地球物理学进展, 17(4):564-574. |
[38] | 王建民,陈树民,苏茂鑫,王元波,王丽娜,关昕,田钢,石战结. 2007.近地表高频补偿技术在三维地震勘探中的应用研究[J]. 地球物理学报, 50(6):1837-1843. |
[39] | 王赟, 张玉贵, 许小凯.2013.六种不同变质程度煤的最大镜质组反射率与弹性参数的关系, 地球物理学报, 56 (6): 2116-2122, doi :10.6038/cjg20130631. |
[40] | 王伟,高星,李松营,乐勇,胡国泽,历玉英.2012.槽波层析成像方法在煤田勘探中应用--以河南义马矿区为例. 地球物理学报,55(3):1054-1062,doi:10 6038/j.issn.0001-5733.2012.03036. |
[41] | 膝佃波,汪鹏程,王赟,等. 2005.利用叠前KirchhOff积分偏移识别小断裂与低幅度构造[J].地球物理学进展. 20(4):1035-1038. |
[42] | 赵殿栋.2009.高精度地震勘探技术发展回顾与展望. 石油物探, 48(5): 425-435. |
[43] | 张向林, 刘新茹, 李健,等. 2007.我国油气开发监测技术进展[J]. 地球物理学进展, 22(4): 1360-1363. |
[44] | 张正峰, 郭建, 韩立国, 徐善辉, 巩向博.2013.基于波场延拓的变换坐标转换波静校正. 地球物理学报,56 (5): 1676-1685 ,doi :10. 6038/cjg20130525. |