2. 中国石油天然气股份有限公司华北油田分公司, 河北任丘062552
2. PetroChina Huabei Oilfield Company, Renqiu 062552, China
二连盆地朝克乌拉凹陷位于苏尼特隆起东部, 是白垩系沉积盆地, 油气资源丰富。本区火成岩多期次喷发, 新生界的巨厚玄武岩广泛分布于地表; 受新生代的地壳变动和风化剥蚀作用, 地形起伏剧烈, 地震地质条件十分复杂, 激发接收条件极差, 地震资料能量弱、信噪比低、成像困难, 曾一度被称为地震“勘探禁区”。华北油田1984—2005年曾先后采用井炮、可控震源开展过多次采集技术攻关, 但始终没有取得突破性进展。2015年, 华北油田在火成岩覆盖区边缘钻遇优质烃源岩, 证明本区具有巨大油气勘探潜力, 2016年再次对本区开展了地震勘探技术攻关。
表层火成岩覆盖区地震勘探属世界级技术难题, 已经引起了勘探学家的普遍关注和高度重视。安学勇等[1]对朝克乌拉凹陷地震采集方法进行了研究总结, 认为火成岩覆盖区“地震波散射严重和难以产生反射波这两种因素无法克服, 而地震波屏蔽严重和激发接收困难这两种因素, 可通过采取适当的措施, 降低其对地震资料的影响”。张爱印[2]对宣化—下花园煤田火成岩覆盖区地震勘探技术进行了研究, 认为采用大药量激发、高覆盖接收、精细静校正和精细去噪方法, 能够取得较好的效果。裴正林等[3]采用P波广角反射消除火成岩高速层对地震波的屏蔽作用。吴希光等[4-6]对复杂地区地震资料低信噪比的原因及对策进行了研究总结, 认为在复杂地表区, 强散射干扰是导致记录信噪比低的根本原因。李灿苹[7]对散射波特征与非均匀地质体对应关系进行了详细分析研究, 认为地震波的散射程度与非均匀地质体的尺度和波长的比值有关。本文通过对二连盆地朝克乌拉凹陷火成岩覆盖区地震地质特点的分析, 在总结以往相关研究成果的基础上, 阐明了造成该区地震原始资料品质差的主要原因和地震勘探的技术难点, 提出了一套以高精度可控震源低频激发、宽线、长排列、高覆盖接收为主的地震资料采集技术和以综合静校正、多域去噪和高精度速度分析为主的处理方法, 展示了该方法技术的应用效果, 希望对类似地区的地震勘探工作具有借鉴意义。
1 难点分析 1.1 地震地质条件研究区为朝克乌拉凹陷火成岩覆盖区, 位于该凹陷朝三次洼中部以西, 区内发育有下白垩系巴彦花群的赛汉组、腾二段、腾一段、阿尔善组地层, 各层底界反射能量中等, 连续性较差。其中腾二段、腾一段和阿尔善组地层为该区良好的储集层, 也是该区勘探的主要目的层, 勘探深度为800~3 000 m。区内火成岩覆盖于近地表, 有的出露地表, 有的隐伏于近地表之下, 表面粗糙、成不规则块状, 内部多孔, 分布零散、破碎, 其间有第四系黄土不同程度填充(图 1)。根据以往地质露头、电法等资料推测凹陷内表层火成岩夹沉积岩总厚度东部一般在50 m左右; 中部厚度变化较大, 在50~200 m; 西部火成岩普遍较厚, 大部分地区在200 m以上(图 2), 分布极不均匀。根据表层调查结果, 本区近地表低速层速度一般为650~1 100 m/s, 高速层速度一般为2 100~3 000 m/s, 最高速度可达4 500 m/s(图 3)。
2016年以前, 在研究区共开展过4次地震勘探。图 4展示了1998年和2005年相邻两条火成岩区二维地震剖面。1998年利用AMG P23可控震源施工, 主要采集参数:振动台次为4台×8次、扫描长度为10 s、扫描频率为13~64 Hz、道距为50 m、覆盖次数为30次。在火成岩覆盖区地震记录信噪比低, 有效信息微弱, 剖面上(图 4a)难以见到有效目的层反射, 缺少低频信息。2005年利用井炮施工, 采用25 m固定井深、24 kg药量、25 m道距、180次覆盖, 地震剖面(图 4b)效果无明显改善。
根据本区地震地质条件及以往勘探经验, 认为本区地震勘探主要有以下难点。
难点一:散射波干扰严重, 提高地震资料信噪比困难。从研究区典型的一条原始地震记录(图 5)上可以看出, 散射干扰遍布整个地震记录, 尤其是近道散射能量强, 远道能量衰减快, 面波、有效波都发生了强烈的散射, 很难见到有效反射波同相轴。杂乱分布的不规则块状火成岩与地表(图 1)以及埋藏于近地表的火成岩与第四系填充的黄土之间都存在不规则波阻抗界面, 地震波在传播过程中遇到不规则波阻抗界面发生散射, 主要表现在以下3个方面。
1) 由炮点下行的地震波在通过近地表巨厚火成岩覆盖区时, 遇到不规则波阻抗界面发生散射, 使下行波到达目的层反射界面的能量降低, 激发效果变差。
2) 由反射界面上行的反射波通过近地表巨厚火成岩覆盖区时, 遇到不规则波阻抗界面发生散射, 使上行反射波到达接收点的能量降低, 接收效果变差, 上行波的散射产生的上行波被接收后降低了资料信噪比。
3) 面波在近地表传播时, 遇到地表不规则波阻抗界面产生散射。由于面波能量一般较强, 强面波散射严重影响地震记录信噪比, 如图 5所示的近道强干扰即是面波散射的结果。
难点二:巨厚火成岩屏蔽作用严重。用水钻对近地表火成岩覆盖区钻孔发现, 地表以下火成岩呈块状杂乱堆积, 孔隙发育, 漏水严重。由于表层大量孔隙的存在, 使之与下伏地层形成很强的波阻抗界面, 激发产生的地震入射波到达这一界面时, 产生的反射波能量强, 透射波能量小, 造成目的层反射能量弱, 巨厚火成岩严重的屏蔽作用, 影响目的层成像效果。
难点三:静校正问题突出。火成岩覆盖区地表起伏剧烈(图 6), 从CK16-11线近地表结构调查模型(图 7)可以看出, 表层速度、厚度横向变化大, 静校正问题严重, 影响成像效果。
从上述难点中可知, 强散射波干扰是造成研究区地震资料信噪比低的最主要原因; 巨厚火成岩散射和屏蔽作用是造成目的层反射能量弱的主要原因; 静校正问题突出进一步影响了研究区目的层的成像效果。
2 技术对策针对上述难点, 在地震资料采集环节采用可控震源低频激发、宽线、长排列、高覆盖接收提高原始资料品质; 在地震资料处理环节, 采用多域去噪、综合静校正和高精度速度分析技术提高剖面成像效果。
2.1 可控震源低频激发技术在地震勘探中, 介质的非均匀性是针对地震波波长而言的[7-8], 一定尺度的非均质体对短波长的地震波而言是非均匀介质, 而对长波长而言可等效为均匀介质, 波长越长, 被等效为均匀介质的地质体越大。根据这一理论, 低频长波长信号穿透非均匀介质的能力强, 散射作用弱; 相反, 高频信号穿透能力弱, 散射作用强。图 8对比了二维不规则界面不同主频正演模拟地震记录的结果。由图 8b可见, 不规则界面L1的散射波很弱, 几乎看不到连续的散射波; 由图 8c和图 8d可以看出, 大量连续的散射波形成背景干扰, 同时使L2界面反射产生扭曲, 这表明低频信号较高频信号抗散射能力强。另外, 低频信号下传的能量较高频信号强, 有利于克服巨厚火成岩屏蔽作用, 提高深层信号能量。因此, 用低频激发技术降低表层巨厚火成岩散射和屏蔽作用的影响, 提高原始资料品质是可行的。
由于工区表层火成岩与黄土混杂、疏松, 井炮激发钻井困难, 成功率低, 激发效果差, 以往井炮采集结果也证实了这一点。EV56低频可控震源具有吨位大、起始扫描频率低、扫描频带宽的特点, 并且其振板设计特殊, 与地面耦合效果和弹性波转换效能较常规震源有明显提高, 因此我们施工时采用EV56低频可控震源, 扫描频率为1.5~64.0 Hz。图 9为可控震源不同扫描频率的叠加剖面。从图 9a可以看出, 扫描频率为12~64 Hz所得的剖面从浅到深难以见到有效地层反射信息, 从图 9b中可以看出, 扫描频率为1.5~64.0 Hz所得的剖面尽管分辨率较低, 但可以清晰见到主要目的层的反射信息。
根据地震波传播原理, 当入射角大于临界角时发生广角反射[9-11], 广角反射的能量一般大于常规反射的能量(图 10)。
采用广角反射增大炮检距的观测方式, 即可克服高速层屏蔽作用, 同时由于广角反射信号能量强, 有利于低信噪比地区资料的获得。广角反射临界角计算公式为:
$ \theta = \sin{ ^{ - 1}}({v_1}/{v_2}) $ | (1) |
式中:θ为临界角; v1为第1层层速度; v2为第2层层速度。
临界排列长度计算公式为:
$ {X_{ \max }} \ge 2h \times \tan \theta $ | (2) |
式中:Xmax为炮检距; h为第1层厚度。
图 11是根据朝克乌拉凹陷的表层和深层资料建立的正演模拟速度模型, 根据模型计算主要目的层腾二段底(T6)发生广角反射的炮检距为1 600 m, 阿尔善组顶(T8)发生广角反射的炮检距为3 500 m, 阿尔善组底(T11)发生广角反射的炮检距为4 000 m。
由火成岩覆盖区(左支)4 000 m炮检距的正演模拟地震记录(图 12)可见, 能产生阿尔善组底以上地层的广角反射, 地震记录正演结果与理论计算结果一致。
图 13是20~40 Hz分频地震记录, 可以清楚地看到广角反射现象, 且广角反射较远炮检距反射能量更强。处理地震资料时可采用少切除或不切除和高阶动效正处理技术, 以更好地利用广角反射信息。
宽线观测具有组合效应, 可以有效压制侧面干扰。高覆盖技术是地震勘探提高信噪比和目的层能量的主要技术。火成岩覆盖区受表层火成岩影响, 资料信噪比极低(如图 4), 宽线、高覆盖结合是有效解决目的层能量和信噪比的有效技术。
图 14显示了火成岩覆盖区采用不同覆盖次数的照明效果。从图 14中可以看出, 随着覆盖次数的增加, 目的层的照明效果越来越好。
定量分析不同覆盖次数叠加剖面, 可以得到火成岩覆盖区下伏目的层信噪比随覆盖次数的变化曲线, 结果如图 15所示。从图 15可以看出, 随着覆盖次数的增加, 目的层的信噪比得到提高, 当覆盖次数大于1 000次时, 信噪比增大幅度不明显。
本次朝克乌拉火成岩覆盖区地震采集技术攻关采用了4线4炮120道960次覆盖的观测系统。对比不同覆盖次数叠加剖面(图 16)可以看出, 随着覆盖次数增加资料信噪比得到明显提升。960次覆盖剖面构造特征清楚, 能够满足火成岩区勘探需求。
复杂地表区静校正是地震资料处理的关键技术, 目前单一的静校正方法很难或不可能解决全区静校正问题, 需多种方法联合使用。朝克乌拉火成岩覆盖地区表层有巨厚火成岩, 难以用小折射、微测井等常规近地表调查方法获得精确的近地表结构, 模型法静校正技术不适宜该区。另外, 该区表层结构复杂, 地表起伏剧烈, 表层速度、岩性横向变化快, 火成岩无稳定折射层, 难以划分出统一折射层面, 折射静校正效果不好。因此, 在该区采用高程静校正、层析静校正和综合静校正处理技术。
对比图 17a和图 17b可以看出, 除局部区域高程静校正叠加剖面好于层析静校正叠加剖面外, 整体上层析静校正叠加剖面好于高程静校正叠加剖面。而综合了高程静校正与层析静校正获得的剖面(图 17c)兼有了两者的优势, 整体效果较好。
火成岩覆盖区干扰波类型较多, 除常见的面波、线性干扰、异常振幅和随机噪声干扰外, 在近偏移距区还存在强能量散射干扰(图 18)。
这些干扰波的存在使反褶积过程的自相关函数变差, 影响反褶积的效果。地震资料处理时依照“先低频噪声后高频噪声、先线性噪声后随机噪声、先强能量噪声后弱能量噪声”的去除原则, 根据各种干扰波的自身特点, 分别在炮域、共检波点域、CMP域或OVT(共偏移距矢量片)域等进行分域去噪, 逐步去除各种噪声, 提高资料信噪比。多域去噪的技术流程如图 19所示。
火成岩覆盖区的强能量散射干扰是本次去噪的重点。在原始地震记录上, 近偏移距区散射干扰的能量明显大于远偏移距区(图 20), 特别是在0~1 250 m偏移距内散射波能量约为无散射干扰区道能量的70倍, 因此, 必须对散射波能量进行衰减。常规的地表一致性异常振幅去除技术, 主要是基于多道统计单道压制的原理, 其应用前提是时窗中要衰减的噪声的振幅大于有效信号的振幅。很显然, 火成岩覆盖区的炮集记录中近偏移距区不能满足该应用条件。为此, 采用OVT域处理技术, 先将炮集记录按OVT矢量片的方式进行划分, 并对OVT域道集进行随机重排, 改变强能量干扰道的分布规律, 然后进行地表一致性异常振幅处理, 从而实现对近偏移距区强能量散射干扰的有效压制(图 21)。
通过保幅、保真的渐进式去噪, 各种干扰得到很好的压制, 有效波成分得到突出, 成像剖面整体信噪比得到明显改善(图 22)。
在低信噪比地区, 经常会由于速度分析精度低造成剖面成像质量变差, 甚至会出现假构造现象。BANCROFT等[12]提出了基于共散射点道集等效偏移距的叠前时间偏移速度分析方法, 即依据地震波旅行时的双平方根方程, 采用叠前克希霍夫积分偏移原理, 将地震道按产生的散射点, 在给定的偏移距范围内抽道集, 称为共散射点道集(CSP道集)。CSP道集与CMP道集相比具有较高的信噪比, 利用CSP道集进行速度分析, 具有提高覆盖次数的优势, 能获取高精度的均方根速度场, 从而提高地震剖面质量[13-18]。
图 23是相同点CMP道集和CSP道集速度谱对比结果, 可以看出后者比前者信噪比高, 能量团更集中, 大道集中浅层能看到明显的双曲线同相轴, 速度解释可靠性更高。
图 24是采用某二维测线的CMP道集和CSP道集进行速度分析建立的速度模型。从图 24a可以看出, 浅层速度横向变化平缓, 而深层速度横向变化相对剧烈, 这种现象与本区地层发育和速度变化规律都不符, 而从图 24b可以看出, 由浅到深速度的变化更符合本区地层发育和速度变化规律。图 25是与图 24对应的叠加剖面, 可以看出, 图 25b所示剖面波组特征清晰, 更符合本区地层发育规律。
本次攻关野外采集采用EV56低频高精度可控震源, 振动台次为2台1次, 扫描长度为14 s, 扫描频率为1.5~64.0 Hz, 驱动幅度为65%, 采用4线4炮120道的观测系统, 道距为50 m, 炮点距为50 m, 覆盖次数为960次。地震资料处理过程中在完成前文介绍的综合静校正、叠前系列去噪的基础上, 开展基于散射波的速度分析, 获得了高质量的速度场。以此开展真地表叠前时间偏移成像攻关, 取得了较为理想的效果。图 26是朝克乌拉火成岩覆盖区相距2 km的老新成果剖面对比结果。从图 26可以看出, 新的叠前时间偏移剖面信噪比明显提高, 主要目的层反射清楚, 地质现象清晰, 剖面质量较老剖面有很大提高。
朝克乌拉凹陷火成岩覆盖区地震勘探技术的再次攻关研究结果表明:
1) 近地表散射是造成本区地震资料反射能量弱、信噪比低的最主要原因, 严重的静校正问题进一步影响了研究区目的层的成像效果。可控震源低频激发是取得攻关突破的关键技术, 宽线、长排列、高覆盖接收是提高原始资料信噪比的有效技术手段;
2) 多域去噪处理技术是提高地震成像剖面信噪比的关键, 基于散射波理论的高精度速度分析技术是提高低信噪比资料速度分析精度的有效措施。
本次地震勘探攻关虽然取得了突破性进展, 所采用的配套技术对类似地区具有重要借鉴意义, 但还有一些技术难题有待进一步攻关研究:
1) 成像剖面上主要目的层发育及展布清楚, 能够满足构造解释需要, 但火成岩散射严重, 对地震资料的分辨率和信噪比的影响还没有完全消除, 需进一步深入研究;
2) 最新地质研究成果表明, 研究区不仅近地表被火成岩覆盖, 深层火成岩也十分发育, 提高成果剖面的品质还需进一步研究深层火成岩成像效果影响因素及先进技术, 比如基于散射波的成像处理技术等。
致谢: 本文在撰写过程中得到了范国增老专家的指导与帮助, 在此表示衷心感谢![1] |
安学勇, 李六五, 于培峰. 二连盆地火成岩地区地震勘探采集方法研究与应用[J]. 石油物探, 2004, 43(2): 171-175. AN X Y, LI L W, YU P F. Research and application of seismic exploration and acquisition methods in igneous rock area of Erlian Basin[J]. Geophysical Prospecting for Petroleum, 2004, 43(2): 171-175. DOI:10.3969/j.issn.1000-1441.2004.02.017 |
[2] |
张爱印. 火山岩覆盖区三维地震勘探方法探究[J]. 中国煤炭地质, 2010, 22(6): 55-58. ZHANG A Y. Probe into 3D-seismic prospecting in igneous rock covered regions[J]. Coal Geology of China, 2010, 22(6): 55-58. DOI:10.3969/j.issn.1674-1803.2010.06.14 |
[3] |
裴正林, 牟永光. 火成岩区地震波传播规律研究[J]. 石油物探, 2004, 43(5): 433-437. PEI Z L, MOU Y G. Study of seismic wave propagation regularity in igneous rock area[J]. Geophysical Prospecting for Petroleum, 2004, 43(5): 433-437. DOI:10.3969/j.issn.1000-1441.2004.05.004 |
[4] |
吴希光, 李亚林, 张孟, 等. 复杂地区地震资料低信噪比的原因及对策(一)[J]. 天然气工业, 2012, 32(1): 27-32. WU X G, LI Y L, ZHANG M, et al. Causes of and solutions to a low signal-to noise ratio of seismic data incomplex areas, Part Ⅰ[J]. Natural Gas Industry, 2012, 32(1): 27-32. DOI:10.3787/j.issn.1000-0976.2012.01.005 |
[5] |
吴希光, 李志荣, 李正佳, 等. 复杂地区地震资料低信噪比的原因及对策(二)[J]. 天然气工业, 2012, 32(2): 38-42. WU X G, LI Z R, LI Z J, et al. Causes of and solutions to a low signal-to noise ratio of seismic data incomplex areas Part Ⅱ.[J]. Natural Gas Industry, 2012, 32(2): 38-42. DOI:10.3787/j.issn.1000-0976.2012.02.008 |
[6] |
吴希光, 李志荣, 张晓斌, 等. 复杂地区地震资料低信噪比的原因及对策(三)[J]. 天然气工业, 2012, 32(3): 57-61. WU X G, LI Z R, ZHANG X B, et al. Causes of and solutions to a low signal-to noise ratio of seismic data incomplex areas, Part Ⅲ[J]. Natural Gas Industry, 2012, 32(3): 57-61. DOI:10.3787/j.issn.1000-0976.2012.03.013 |
[7] |
李灿苹.散射波特征与非均匀地质体对应关系研究[D].北京: 中国地质大学, 2006 LI C P.A study on corresponding relationship between scattering wave characteristics and heterogeneous geologic bodies[J].Beijing: China University of Geosciences, 2006 |
[8] |
范小平, 杨从杰, 李清河. 介质非均匀性参数对散射波包络的影响[J]. 地球物理学进展, 2013, 28(2): 687-694. FAN X P, YANG C J, LI Q H. The effect of inhomogeneiety parameters on scattering wave fnvelopes[J]. Progress in Geophysics, 2013, 28(2): 687-694. |
[9] |
周芝旭, 董幼瑞. 火成岩区广角-转换波地震勘探方法研究[M]. 北京: 石油工业出版社, 2005: 1-103. ZHOU Z X, DONG Y R. Wide-angle-converted wave seismic exploration method in igneous rock area[M]. Beijing: Petroleum Industry Press, 2005: 1-103. |
[10] |
徐文君, 於文辉, 胡中平. 广角反射波的特征及正演模拟[J]. 石油地球物理勘探, 2006, 41(4): 390-395. XU W J, YU W H, HU Z P. Feature and forward simulation of wide-angle reflection[J]. Oil Geophysical Prospecting, 2006, 41(4): 390-395. DOI:10.3321/j.issn:1000-7210.2006.04.007 |
[11] |
王熙明, 王志, 贺振华, 等. 广角反射方式下观测系统设计与优化[J]. 物探化探计算技术, 2002, 24(4): 298-303. WANG X M, WANG Z, HE Z H, et al. The observation system design and optimization of the wide angle reflection[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2002, 24(4): 298-303. DOI:10.3969/j.issn.1001-1749.2002.04.003 |
[12] |
BANCROFT J C, GEIGER H D. Equivalent offsets and CRP gathers for prestack migration[J]. Expanded Abstracts of 64th Annual Internat SEG Mtg, 1994, 672-675. |
[13] |
刘开隆.地震散射波EOM成像技术应用研究[D].西安: 长安大学, 2015 LIU K L.An application research of seismic scattering wave EOM imaging method[D].Xi'an: Chang'an University, 2015 |
[14] |
王伟, 尹军杰, 刘学伟, 等. 等效偏移距方法及应用[J]. 地球物理学报, 2007, 50(6): 1823-1830. WANG W, YIN J J, LIU X W, et al. Equivalent offset method and its application[J]. Chinese Journal of Geophysics, 2007, 50(6): 1823-1830. DOI:10.3321/j.issn:0001-5733.2007.06.024 |
[15] |
徐基祥.山地地震勘探中近地表散射波压制技术[R].SPG/SEG北京2016国际地球物理会议, 北京, 2016 XU J X.The suppression technique of near-surface scattered wave in mountain seismic exploration[R].SPG/SEG Beijing 2016 International Geophysical Conference, Beijing, 2016 |
[16] |
王要森.地震散射波正演模拟及去噪方法研究[D].青岛: 中国石油大学(华东), 2012 WANG Y S.Research of scattered wave forward simulation and denoising technology[D].Qingdao: China University of Petroleum(East China), 2012 http://kns.cnki.net/KCMS/detail/Error.aspx |
[17] |
宋翔宇, 李振春, 周卿, 等. 共散射点道集映射噪声压制方法及其应用[J]. 石油物探, 2013, 52(5): 524-529. SONG X Y, LI Z C, ZHOU Q, et al. Mapping noise suppression in common scattering point gather and its application[J]. Geophysical Prospecting for Petroleum, 2013, 52(5): 524-529. DOI:10.3969/j.issn.1000-1441.2013.05.012 |
[18] |
王勇, 朱亚平, 杨慧珠. 共散射点成像及其在低信噪比三维地震数据处理中的应用[J]. 石油地球物理勘探, 2000, 35(1): 20-26. WANG Y, ZHU Y P, YANG H Z. Common scattering point imaging and its application to low S/N 3-D seismic data processing[J]. Oil Geophysical Prospecting, 2000, 35(1): 20-26. DOI:10.3321/j.issn:1000-7210.2000.01.003 |