基于连续小波变换目标处理技术在储层预测中的应用-尹继尧.pdf

《基于连续小波变换目标处理技术在储层预测中的应用-尹继尧.pdf》由会员分享，可在线阅读，更多相关《基于连续小波变换目标处理技术在储层预测中的应用-尹继尧.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第55卷第3期2016年5月石油物探GEOPHYSICAL PROSPE(?nNG FOR PETROLEUMVoL 55，No3May，2016尹继尧，钟磊，张吉辉，等基于连续小波变换目标处理技术在储层预测中的应用EJ石油物探，2016，55(3)：433440YIN Jiyao，ZHONG Lei，ZHANG Jihui，et a1Target processing by continuous wavelet transform coefficients applied to reservoir predictionJGeophysical Prospecting for Petroleu

2、m，2016，55(3)：433440基于连续小波变换目标处理技术在储层预测中的应用尹继尧，钟磊，张吉辉，赵俊，张田(中国石油天然气股份有限公司新疆油田分公司勘探开发研究院，新疆克拉玛依834000)摘要：随着油气藏勘探开发的不断深入，预测隐蔽性复杂岩性油气藏的空间展布规律变得越来越重要。采用基于傅里叶算法的连续小波变换信号分解和重构技术，对目的层段进行地震资料解释性目标处理，以提高储层预测精度。薄层模型正演模拟和两个实例应用结果表明，该方法在保持地震资料信噪比的基础上能有效提高地震资料中具有地质意义的弱信号，突出薄互层内部细节，其结果有利于精细刻画地层岩性油气藏的空间展布，提高薄储层预测精度

3、。关键词：连续小波变换；弱信号；目标处理；储层预测中图分类号：P631 文献标识码：A文章编号：10001441(2016)03-043308 DOI：103969jissn10001441201603014Target processing by continuous wavelet transform coefficientsapplied to reservoir predictionYIN Jiyao，ZHONG Lei，ZHANG Jihui，ZHAO Jun，ZHANG Tian(Petroleum Exploration and Development Research Inst

4、itute，PetroChina Xingjiang Oilfield Company，Karamay 834000，China)Abstract：As the deepening of hydrocarbon exploration and development，the prediction of the spatial distribution of subtlecomplex lithologic reservoir becomes more and more importantTherefore，w。e adopted the signal decomposition and rec

5、onstruction technique by continuous wavelet transform based on Fourier algorithm to carry out interpretative target processingfor seismic data and tO improve the precision of reservoir predictionThe forward modeling and two application cases indicatethat the method can effectively strengthen the wea

6、k signals with geological meaning while preserve the SN of seismic data，which highlights the inner details of thin interbedded layersMoreover，the results are beneficial for the fine characterizationon the spatial distribution of stratigraphic and lithologic reservoirs and improve the precision of re

7、servoir predictionKeywords：continuous wavelet transform(CWT)，weak signal，target processing，reservoir prediction随着油气藏勘探开发的深人，物探技术正经历着由构造油气藏勘探到地层、岩性油气藏勘探和油藏地球物理描述的发展方向转变，研究任务也由简单的构造成像转为隐蔽性复杂岩性预测，甚至储层孔渗特征描述、油藏静态描述和动态监测等。相应的解释需求发展方向也由叠后走向叠前，从单一走向综合。勘探和评价开发的精细化对储层预测的精度要求越来越高，更需要精细物探技术的支持和处理解释的一体化，对地震基础资料

8、精细程度的要求也在不断提高。但是由于受检波器、信号处理、野外采集环境等因素影响，地震资料主频通常偏低，一些地层的厚度远低于当前地震勘探的分辨范围。因此，如何在现有地震资料基础上针对目的层进一步提高地震资料品质，对于更加精确预测储层收稿日期：20150702；改回日期：20151030。作者简介：尹继尧(1981一)，男，博士，高级工程师，主要从事地震地质综合研究。基金项目：中国博士后科学基金项目(2013M542414)资助。This research is financially supported by China Postdoctoral Science Foundation(Grant

9、 No2013M542414)万方数据434 石油物探 第55卷的空间展布规律是一项具有重要意义的工作。目前，常规目标处理通常需要兼顾浅、中、深层地震剖面的整体效果，往往无法有效突出由于薄互层等影响造成的有用地质弱信号，这对后续地震解释中的反演和属性提取等储层预测工作带来非常大的影响，因此，必须高精度处理解释一体化，井震结合，切实提高目的层的有用地质弱信号，提高储层预测精度。近年来迅速发展的小波分析方法是一种具有多分辨功能的时频分析方法，被誉为数学显微镜，在石油勘探、图像处理和语音分析等众多领域都有广泛应用。路鹏飞等1利用能量准则和相似系数将小波变换与振幅谱补偿相结合，在保持较高信噪比前提下，

10、给出了薄互层高分辨率分析方法。高静怀等2采用定性分析与Monte Carlo仿真相结合的方法，研究了在物理小波为基本小波的变换域随机噪声局部小波功率谱的统计特性，利用假设检验给出了在给定置信度下确定有效信号在小波变换域内能量分布空间的方法，通过在有效信号能量分布空间重建信号来压制噪声。王真理3在研究薄互层的储层预测问题时提出特征分辨率和储层特征信号的概念，同时利用二进正交小波变换，结合相似系数法提高薄互层分辨率。李春峰等4根据小波变换特征提出奇异性指数，能够准确反映数据中奇点位置和奇异性强度，可以作为一种精确刻画地层边界的地震属性。熊晓军等5利用广义S变换通过提取单频剖面或单频曲线方法，研究如

11、何提高地震波主频。岳文正等6根据小波变换特点提出一种识别储层流体新方法。罗永健等7基于小波变换多分辨率特性和调焦作用，提出利用小波在薄层地震记录上检测结构特性和估算薄层厚度新方法。高静怀等8提出三参数小波，通过选择这3个参数，可使其最适合于所给定问题，并通过对三参数小波的参数加以约束，获得新近似解析小波，论证了通过恰当选择参数，三参数小波不仅适合于分析包含慢变频率和振幅分量信号，而且也适合于包含快变分量信号；并以三参数小波为分析小波，提出了一种用于薄互层的地震资料分析方法。此外，在地震资料去除面波干扰、混合相位子波模拟、噪声去除及储层预测方面小波变换也得到了广泛的应用睁13。上述给出的小波分析

12、方法，大部分都是在离散小波变换基础上进行的信号分析和重构，无法针对所需要的细节信号进行任意尺度和任意空间上的操作和分析，也不能精确针对特定频率段的子波特征增强弱信号。由于调谐频率等存在，有用地质弱信号往往存在于特定频率段中，必须进行连续小波变换的重构，而由于连续小波为冗余变换，存在着逆变换的定义非唯一性，不能直接进行逆变换。基于以上原因，本文采用基于傅里叶算法的连续小波变换进行信号的分解和重构，达到连续小波变换逆变换的目的，同时有效结合目标层段的地震反射特征和地质特征，并且井震结合，有针对性地进行解释性目标处理，从而提高地震资料的品质，提高储层预测精度。1 方法和原理连续小波变换(CWT)表达

13、式145为：cb，删一D去驴(学)出(1)式中：a为缩放因子(对应于频率信息)；6为平移因子(对应于时空信息)；曲为小波函数(又叫基本小波或母小波)；“*”表示复共轭。连续小波变换是一种类似于信号基于频率的滤波处理，将公式(1)改写成逆傅里叶变换，得到如下公式16：C(a，6；厂()，妒()一麦I fGo)口(妒(黝)。ej曲do) (2)式中：函()和于(叫)分别为子波和信号的傅里叶变换。逆连续小波变换通常为双积分形式(doubleintegral form)，为使之能进行唯一逆变换，需将其变为单积分形式(singleintegral form)，逆CWT单积分公式为：f(t)一2Re士序M

14、舭)da(3)通过(3)式我们可以依靠增强CWT中包含有用地质弱信号的特定尺度信号到原始信号中，突出目的层的有效信息，从而达到提高目的层储层预测精度的目的。目前有许多小波函数可供选择，这些函数具有不同的时频特征，根据研究目的选择合适的分析小波函数就显得非常重要。Morlet小波具有较好频率分析能力，很多学者在使用小波变换进行相关地震属性分析时，大都采用Morlet小波函数，因此本文选取该小波函数进行研究。基于连续小波变换解释性目标处理流程如图1所示，具体过程和步骤如下：1)对目的层段进行地质分析，确定目的层段万方数据第3期 尹继尧等基于连续小波变换目标处理技术在储层预测中的应用435岩电特征，

15、井震结合精细标定，特别对含油层段反复确认标定，综合确定目的层段储层地震地质特征；2)对原始地震资料进行全频段子波分解和重构，首先对有效频带范围内的子波以窗长5 Hz，步长10 Hz滑动扫描方案进行分解和重构，把重构后的子波段叠加到原始资料中得到多个新数据体(按频宽575 Hz计算，将会得到14个新数据体)，对该批新数据体再进行新一轮标定、属性分析、反演等验证分析，选取标定效果最好、反演和属性分析结果最符合区域基本地质认识和储层特征描述的数据体为最终数据体，最后再对该数据体所对应的子波频率段以窗长1 Hz，步长5 Hz进行精细扫描，选取最优方案。图1 基于连续小波变换解释性目标处理流程经过上述步

16、骤处理后的结果既保留了原始地震资料所有信息，同时又有效突出目的层段储层地震响应特征，处理结果为目的层段的地震解释和储层预测提供了更好的基础地震资料。2正演分析为检验基于傅里叶算法连续变换解释性目标处理技术对提高储层预测精度的有效性，设计模型如图2所示，图中黄色代表砂体，上部为一套厚度10-20 m的薄砂层，速度2 800 ms，薄砂层上部速度为2 400 ms，下部速度为3 200 ms，底下为两套呈叠置关系的高速透镜砂体，厚度为010 m，速度为3 600 ms，两者之间距离最小3 m，最大6 m。利用50 Hz Ricker子波合成地震记录，所得结果如图2中的波形所示。从图2可以看到，由于

17、分辨率较低，使得薄砂层中部出现复波，底下两套透镜砂体由于间隔小于A4，形成砂泥岩薄互层，使得底下一套透镜砂体中间部分的顶界面由于受到上层透镜砂体的强反射及折射影响，反射信息较弱，在原始合成地震记录上无法单独对这两套透镜体进行精确识别和描述，而薄层单砂体边界的精细刻画对于优化油田开发方案和井位部署具有非常重要的作用和意义。利用基于傅里叶算法连续小波变换对信号进行子波特征分解和重构，为比较本文所使用的连续小波变换结果和常规离散小波变换结果的差异，同时对原始合成地震记录进行离散小波变换信号重构，所得结果如图3所示。其中，图3a为原始合成地震记录，图3b为基于傅里叶算法连续小波重构后的合成地震信号，图

18、3c为离散小波变换重构后的合成地震信号。对比图3a，图3b和图3c可知，基于傅里叶算法连续小波变换信号重构能够很好地分解地质模型中薄砂层的复波，同时突出了底部薄互层的内部细节，底下透镜砂体顶界面弱反射信息得到有效加强，为砂体展布刻画和储层预测提供了更好品质的地震资料，且保持了较好的地震信噪比和保幅性；而离散小波变换重构后得到的合成地震信号，虽然同样也能区分复波和突出相应薄互层内部细节，但不能有效控制特定频率段子波重构。因此，在信号重构的过程中，如要获得和连续小波变换相同分辨率，不可避免地会带来更多高频噪声万方数据436 石油物探 第55卷2040g6080lOO图2 薄层地质模型和合成地震记录

19、叠合结果道号0 20 30 40 50、一、-：；：道号0 20 30 40 50“、h、huu、“t、：。f、：fjji：(LLLLkhn“I“(二(【：o，M仍留；刀，：!40、6080100道号O 20 30a b c图3不同类型小波变换解释性目标处理效果对比a原始合成地震记录b基于傅里叶变换算法连续小波变换地震记录重构效果；c基于离散小波变换地震记录重构效果干扰，使得资料的信噪比降低，且保幅性能也相对较差。分析上述正演结果可知，本文所使用的基于傅里叶算法连续小波变换地震信号重构能提高目标层段储层预测精度。3实际资料分析31金龙井区实例准噶尔盆地西北缘金龙井区二叠系上乌尔禾组与佳木河组为

20、区域性地层不整合接触，上乌尔禾组地层自东南向西北逐层超覆沉积于二叠系佳木河组和石炭系地层之上，根据岩电特征自下而上可分为乌一段、乌二段和乌三段，其中乌一段、乌二段为主要的目的层段。乌一段和乌二段属于砂岩与泥质粉砂岩互层的河口砂坝，储层相对较薄，为岩性构造油藏，油水关系较为复杂。该工区实际三维叠后地震资料目的层为二叠系上乌尔禾组，主频仅为25 Hz，采样间隔2 ms，通过井震结合，确定了需要分辨的目标层段储层厚度。根据前述原理和方法，针对目的层段进行基于连续小波变换子波特征弱信号增强处理，得到重新处理后的三维数据体。图4为连续小波变换信号重构前、后连井地震剖面。从图4可以看出：解释性目标处理后，

21、地震剖面同相轴连续性增强，如图中黑色虚线所示(黑色虚线为百口泉组顶界面)，处理后使得该界面横向更易解释；从图4虚线椭圆中可以看到，很多复波得到分解，薄互层内部细节更加突出，断点更加清晰；此外，从图4b频谱图上可以看出，基于连续小波变换解释性目标处理后，在不改变原始频率有效频宽和低频情况下，提高了资料的主频。为进一步评价解释性目标处理对提高储层预测精度效果，本文对处理前、后数据体在相同井控和反演控制参数条件下进行稀疏脉冲反演，结果如图5所示。图5中白色、蓝色偏泥岩，红、黄、绿色偏砂岩。图5剖面中曲线为自然电阻率测井曲线，对比处理前、后反演结果，可以清楚看到，解释性目加g越卷加量、万方数据第3期

22、尹继尧等基于连续小波变换目标处理技术在储层预测中的应用 4372道号ab图4金龙井区某连续小波变换信号重构前(a)、后(b)的连井地震剖面(左下图为剖面对应的频谱分析结果)2 5002 600l 2 700：2 8002 9002 5002 6002 700：2 8002 900a道号b图5在相同井控和反演控制参数条件下数据体处理前(a)、后(b)的稀疏脉冲反演结果万方数据438 石油物探 第55卷标处理后分辨率得到了明显的提高，砂体响应关系更好，并且砂体展布更符合地质实际。另外，试油结果表明，位于高部位的X24井为含油水层，而位于低部位的X25井为油水同层，从目标处理后的结果明显可以看出，X

23、24井含油水层段属于单独一套透镜砂体(图中红色箭头所示)。可见，解释目标处理突出了有效地质弱信号，提高了储层预测精度，对该区域的油水关系复杂问题进行了有效解释。本文处理后的数据以及相应反演结果已作为基础资料用于该研究区域内油藏目标评价工作，为后续多口评价井的布设和开发部署方案优化提供了重要参考依据。32玛北斜坡区实例玛北斜坡区主要目的层段为三叠系百口泉组，位于扇三角洲前缘相带，砂体纵、横向变化快，叠置特征突出，有利相带边界和相对优质储层的精细刻画难度大。三维叠后地震资料采样间隔2 ms，百lZl泉组主频仅为30 Hz，对该数据进行解释性目标处理，得到连续小波变换重新处理后的三维数据体。图6为连

24、续小波变换信号重构前、后连井地震剖面。从图6可以看到：明显增强了有用弱信号，目的层百口泉组内幕地震反射信息更加丰富，地层接触关系更加清晰，更易于刻画砂体叠置特征；很多复波被分解，薄互层内部的细节更突出；处理后的地震剖面同相轴连续性增强。相对于原始叠后资料，处理后的地震数据在保持原有波组关系的基础上获得了更高的分辨率，复波被分解，地震资料的信噪比也有所保持，处理后的地震数据基本上具有振幅相对保持的特点，有利于分辨和解释薄互层。 哆一”_-_一 M 卅， 一M二p 。- 一 l_r j Pr一。 一、，：=；p_r，一一 ，一-：-!?皇一d+，-!-= P一r一 。-，r一蚺”一飞一一茎至豸。

25、-二：一忖：二一 。：蕊 o，r-_一，一，一=”， ， r哩y一一：J 一一哆 巳!v；”1多 r一；乡1 t口f 一 _-p一Z 一_r。、，：y” 誊一一|-i；”噬=_， pV_r夕 -一，-，一 ， t，_一，一：， 乞b图6玛北斜坡区某连续小波变换信号重构前(a)、后(b)连井地震剖面为进一步评价处理后地震资料的效果和保真性，对处理前、后两套数据体分别进行单井标定和平面地震相分析比较。图7为M3井地震资料目标处理前、后的标定结果。图7中红色框内为目的层段，可以明显看到，目标处理后井旁地震道和合成记录波组关系对应更好，突出了薄互储层的反射弱信号，表明采用本文方法处理后，在保持信噪比前

26、提下，提高了有用地质信号，为精细储层研究奠定了更高品质的资料基础。此外，通过综合测井、地质分析，该区域三叠系百口泉组物源主要来自东北方向的夏子街物源(图8a)。对比图8a，图8b和图8c可知，处理后平均波峰振幅属性所揭示的地震相和沉积相更加吻合，表明处理后突出的弱信号为具有地质意义的有效信号，在此基础上再进行叠前反演等储层预测研究，将会更好地提高该区域的储层预测精度。咖啪枷姗2222咖m瑚姗2222万方数据第3期 尹继尧等基于连续小波变换目标处理技术在储层预测中的应用 439一 _处理前道号416 1 418 1 420 l 422 1 416J P I lII、兰心。姜二三声广一，百、广、梦

27、，_)j，) )l、广、r、ll，-，夕I l I l I l lIl谚谚沙萝声潍一，)，J，)、I，l J I l I II lI I Jl l抄，【I，ll，J，J l，lk、k、I l IJ lJ萼豳睁缈声)F爹ji净P峥薄萼芗W淞，声，JJl，IJ J J Ik电阻率p(nm)T，kT1b3T，也T。b隧隘迸L图7 M3井地震资料目标处理前、后标定结果_平原辫状河道 I平原砂质碎屑流-前缘砂质碎屑流_I：I：口前缘近岸水下河道口前缘远岸水下河道 口滩坝河口坝口平原河道间 口前缘水下分流问湾口前缘席状砂 圄前扇三角洲、滨浅湖泥a C图8地震资料目标处理前、后地震属性对比a玛北斜坡区百口泉

28、组二段沉积体系分布b百1：3泉组二段平均波峰振幅(处理前)；c百13泉组二段平均波峰振幅(处理后)4 结束语本文采用基于傅里叶算法连续小波变换进行信号分解和重构，达到连续小波变换逆变换的目的。根据目标层段地震反射特征和地质特征，井震结合，针对目的层段进行解释性目标处理，提高了储层预测精度。分析了正演和离散小波变换结果，表明本文方法能很好地分解地质模型中薄层复波，突出砂泥薄互层内部细节，同时能够较好保持地震信噪比和保幅性能；而离散小波变换处理结果降低了资料的信噪比。实例分析结果表明，处理后的数吕、j姗姗伽姗一：瑚咖2222222咖咖啪姗22222万方数据440 石油物探 第55卷据突出了有效地质

29、弱信号，为更好地刻画砂体边界和提高储层预测精度打下了坚实的资料基础。参考文献1路鹏飞，郭爱华，朱德献，等薄互层高分辨率分析方法研究J石油天然气学报，2010，32(3)：6972LU P F，GU0 A H，ZHU D X，et al，The Techniqueof high resolution analysis in thin inter-beddingsJJournal of Oil and Gas Technology，2010，32(3)：69722高静怀，毛剑，满蔚仕，等叠前地震资料噪声衰减的小波域方法研究J地球物理学报，2006，49(4)：11551163GA0 J H，J，M

30、AN W S，et a1On the denoising method of prestack seismic data in wavelet domainl-JChinese Journal of Geophysics，2006，49(4)：115511633王真理地震波场数值模拟的高效算法研究(博士后出站报告)R北京：中国科学院地球物理研究所，1997：4546WANG Z LThe research on efficient algorithms ofthe numerical simulation of seismic wave field(Postdoctoral report)RB

31、eijing：Geophysical ResearchInstitute of China，1997：45464李春峰，LINER C基于小波多尺度分析的奇性指数：一种新地震属性J地球物理学报，2005，48(4)：882-888LI C F，LINER CSingularity exponent from wavele-t based muhiscale analysis：a new seismic attributeI-JChinese Journal of Geophysics，2005，48(4)：8828885熊晓军，贺振华，黄德济，等广义S变换在提高地震高分辨率处理中的应用JI勘探

32、地球物理进展，2006，29(6)：415-418X10NG X J，HE Z H，HUANG D J，et al。Application of generalized S transforlTl in seismic high resolution processing I-JProgess in Exploration Geophysics，2006，29(6)：415-4186岳文正，陶果小波变换在识别储层流体性质中的应用J地球物理学报，2003，46(6)：863869YUE W Z，TAO GA method for recognition of fluid property in

33、reservoirs using；vavelet transformationrJ-Chinese Journal of Geophysics，2003，46(6)：8638697罗永健，俞根苗，张守宏小波变换在薄层地震信号分析中的应用J西安交通大学学报，2001，35(8)：848-850LUO Y J，YU G M，ZHANG S HWavelet transform analysis of thin-bed seismic signalJJournalof Xian Jiaotong University，2001，35(8)：848-8508高静怀，万涛，陈文超，等三参数小波及其在地震

34、资料分析中的应用J地球物理学报，2006，49(6)：1802一18lZGAO J H，WAN T，CHEN W C，et a1Three parameter wavelet and its applications to seismic data processingl,JChinese Journal of Geophysics，2006，49(6)：180218129张华，潘冬明，张兴岩二维小波变换在去除面波干扰中的应用J石油物探，2007，46(2)：147150ZHANG H，PAN D M，ZHANG X YApplication of2：D wavelet transfclrma

35、tion in elimination surfacewave interferenceJGeophysical Prospecting forPetroleum，2007，46(2)：1471。50101王西文，高静怀，李幼铭。高分辨地震资料处理中的导数小波函数的构造J石油物探，2000，39(2)：6471WANGX W，GAO J H，LI Y M Construction ofderivative wavelet function in high resolution seismicdata processingl,JGeophysical Prospecting for Petrol

36、eum，2000，39(2)：64711-11彭才，朱仕军，孙建库，等小波变换域KL变换及其去噪效果分析J石油物探，2007，46(2)：147150：112114PENG C，ZHU S J，SUN J K，et a1K-L transfornlation in wavelet conbersion domain and the analysis ofde-noise effectJ1Geophysical Prospecting for Petroleum，2007，46(2)：112114f12李曙光，徐天吉，甘其刚，等频率域小波变换分频处理在川西地震勘探中的应用J石油物探，2010，4

37、9(5)：500503LI S G，XU T J，GAN Q G，et a1Application of frequency division processing on frequency-domainwavelet transformaion in deep Chuanxi DepressionjGeophysical Prospecting for Petroleum，2010，49(5)：50050313陈广军，张善文，李建明，等小波分析技术在薄砂岩储集层描述中的应用以埕岛地区馆上段为例J石油物探，2002，41(1)：9599CHEN G J，ZHANG S W，LI J M，et

38、a1Applicationof wavelet analysis to the characterization of thinsandstone reservoirs：an example from the uppermember of Guantao Formation in Chengdao arealJGeophysical Prospecting for Petroleum，2002，41(1)：959914DAUBECHIES L Ten lectures on waveletsMPhiladelphia，PA：Society for Industrial and AppliedMathematics，1992：1-377E15FARCE M Wavelet transforms and their applicationto turbulenceJAnnual Review of Fluid Mechanics，1992，24：39545716TORRENCE C，COMPO G PA practical guide towavelet analysisJBulletin of the American Meteorological Society，1998，79：6178(编辑：顾石庆)万方数据