摘要:高精度的合成地震记录是井震联合统层的桥梁。从合成记录根本原理反射角度出发,提出了基于地震波理论,采用非褶积模型的万能高精度合成地震记录制作新方法,即子波反射叠加合成地震记录制作法。这种方法可以制作任意多、任意厚度的合成记录。这种任意厚度就是实际地下的任意地质界面。
加入收藏
传统人工合成地震记录制作是基于褶积理论模型的[1]。设地震子波是s(t),地震记录f(t)是地震子波s(t)与反射系数r(t)褶积的结果。
一说到褶积,就意味着反射系数必须是等间隔的,而等间隔的反射系数与实际地质界面并不是一一对应的,这样制作的合成地震记录的精度就会受影响。
在此基础上,提出了基于地震波理论,采用非褶积模型的万能高精度合成地震记录制作新方法——子波反射叠加合成地震记录制作法。这种方法可以制作任意多、任意厚度的合成记录。这种任意厚度就是实际地下的任意地质界面。
1、理论基础
地震记录上看到的反射波波形是地震子波在地下各反射界面上发生反射时形成的。反射波的振幅有大有小(决定于界面反射系数的绝对值)[2]、极性有正有负(取决于反射系数的正负)、到达时间有先有后(取决于反射界面的深度)的地震反射子波叠加的结果。
应用反射叠加这一最基础的原理,按一定的采样率制作合成地震记录。本文模型的采样间隔均为nms。
假定子波为雷克子波,不考虑波损与透射的影响,第K层的厚度为HK,速度为VK,那么子波到达第K层的时间为:
子波在T时刻发生了反射且在T时刻刚好到达地表h=0处,则第K个反射界面T时刻在地表H=0处的振幅值有贡献,其贡献值取决于T与第K个界面双程旅行时之差,应用这个差值确定在子波中的位置,利用子波的解析式求出子波在该差值时间(任意的、连续的)上的位移再乘上界面的反射系数即为第K个界面对合成记录T处的贡献值,将所有界面的贡献值累计求和即为合成记录总的位移值。那么第K层的反射在T时刻能够到达地表并对地表的振幅值有贡献其条件为:
式中:WL——子波的波长。
我们已经知道在T=2TK,则T时刻将要接收到第K层的反射;如果在T时刻2TK≤T≤2TK+WL,那么第K层对地表反射一直有贡献(如图1所示)。
图1反射法合成记录制作原理图
2、模型建立
2.1任意厚度模型
首先建立一套多层随机厚度模型[3],包括受反射波叠加影响而不好分辨的砂泥薄互层,也包括能很好分辨的砂泥厚层,通过反射法利用40Hz雷克子波对该模型做正演,如图2所示,反射法正演结果与该模型匹配较好。
2.2楔状薄互层模型
图2反射法制作合成记录与反射系数对应关系
首先建立一个在泥岩中嵌入单层砂岩薄层楔状模型,两层泥岩间的砂岩厚度由零逐渐增大到波长的四分之一时间厚度处,即砂岩的厚度在最厚处(H)等于雷克子波的1/4波长(调谐厚度),其余厚度均小于1/4波长。我们在给定相同子波条件下,应用褶积与反射法两种方法对模型进行正演。
利用褶积对模型做的正演,如图4a所示,受调谐厚度作用的影响,模型地震响应特征在砂体最厚处振幅幅度达到最大,随后呈现幅度逐渐减薄的状态,并在模型1/4处砂体分辨达到极限;利用反射法对模型做正演,如图4b所示,其地震振幅响应也总体呈现减弱的状态,但变化比较平稳,在对于砂岩分辨能力,在模型1/2处分辨达到极限;通过对比分析发现,对于薄层砂岩分辨,反射法的正演分辨能力高于褶积正演方法的分辩能力。
图3砂岩楔状模型
图4单层砂岩薄层楔状模型的地震响应
3、应用实例
在某区块,利用反射法制作合成记录,与井间地震道进行对比。本次也是利用声波和密度曲线制作合成记录,分以下几部分完成:
(1)对声波和密度曲线进行多次的圆滑处理;
(2)对声波和密度进行中值滤波,去除曲线上的奇异点;
(3)对声波和密度按一定深度间隔重新采样,采样间隔一般为实际工区薄层的分辨率,这样做是因为一般测井曲线的采样间隔没有做到地质与地震界面的结合;
(4)将采样深度域数据转换成时间域数据;
(5)时间域数据按制作合成地震记录的时间间隔进行重采样,采样间隔一般与地震道数据保持一致;
(6)计算波阻抗,按时间间隔计算反射系数;
(7)与给定子波进行反射叠加制作合成地震记录。
某井目的层在深层,根据其分辨率深度域重采样的间隔给定9m。合成记录制作中所使用的子波为零相位雷克子波,从井旁地震道提取子波,确定其主频为25Hz。如图5所示,该井利用反射法与褶积法做合成记录,其中曲线I为褶积法所做合成记录,曲线Ⅱ为反射法所做合成记录,经对比分析,认为反射法合成记录与反射系数对应关系较好。为其它区块利用该法开展工作提供了理论依据。
图5工区实际地震合成记录对比
4、结论
(1)反射法可制作任意厚度、任意多层的合成记录,其地震响应与模型匹配较好;
(2)反射法与褶积法对比,其对薄层砂岩有更好的分辨能力;
(3)在实际工作应用中,该方法应用效果较好,为其它区块井震联合统层提供了很好的参考依据。
参考文献:
[1]陆基孟.地震勘探原理[M].东营:石油大学出版社,1993.
[2]靳玲,苏桂芝,刘桂兰,朱丕月,李磊,胡艳萍.合成地震记录制作的影响因素及对策[J].石油物探,2004,43(3):267-271.
[3]徐伟慕,郭平,胡天跃.薄互层调谐与分辨率分析[J].石油地球物理勘探,2013,48(5):750-757.
孙丽梅.高精度反射法制作合成记录及其应用[J].西部探矿工程,2020,32(05):61-63.
基金:国家科技重大专项(2011ZX05023-006-002).
分享:
近年来,各种分类或聚类算法应用于地球物理资料处理和反演,主要包括模糊C均值(FCM)聚类[1]、支持向量机(SVM)[2-3]和神经网络(NN)[4-5]等。其中,SVM和NN算法均属于监督学习分类方法,需要足够数量的样本维持结果的稳定性;而FCM聚类算法是一种无监督学习方法,具有样本量需求少的优势[6-7],更适用于地球物理方法联合反演。
2024-08-05咸潮是发生在滨海河口地区的一种水文现象,它是由于太阳和月球对地表海水的引力作用,海水倒灌进入内陆河道,导致水体变咸。咸潮危害巨大,例如影响水厂供水安全及部分水生生物的正常生存。因此,加强对咸潮的监测,并据此及时制定应对措施,对当地防灾减灾和城市供水保障等具有重要意义。
2024-01-02在新时期背景下,对国土空间规划功能进行准确定位,全面落实各项规划方案,不仅可以保护周围生态环境,而且能够实现统筹治理目标。结合现阶段我国国土空间规划功能定位现状可知,环境修复难度不断增加,各项管控措施落实效果较差,严重影响生态环境质量。鉴于此,本文重点分析新形势下国土空间规划功能定位和实施途径。
2021-06-01日益严重的环境污染、自然资源无节制开发,不仅开始制约人类社会的发展,同时威胁人类的生存和繁衍,如何改变这种情况,成为当今人类亟待解决的问题。基于卫星遥感的国土空间监测的发展与应用可以帮助人类全面调查全球自然资源和生态环境,确保人类更合理的开发、分配、利用自然资源,使实现可持续发展成为可能。
2021-05-04文章围绕着国土空间规划中对地理信息大数据的有效应用展开充分讨论,首先简要介绍了地理信息大数据的概述,然后表明了国土空间规划的特点,最后着重探讨了地理信息大数据在国土空间规划中的具体应用,包括提供合理数据信息、优化国土空间规划、强化信息的数据源、明确空间具体位置,以此完善地理信息的内容,并形成重要的数据,促进国土空间规划的良好发展。
2021-04-27在当前这一社会与经济快速发展背景下,地震预警已经成为了人民群众关心的问题,因为地震为社会与经济带来的威胁远超于其他地质灾害,再加上地震的前期预兆不明显,所以就需要利用合理的方式进行监测预警,进而有效提高防控效率,降低灾害为人民群众生命财产带来的威胁。那么在当前这一社会背景下,信息技术成为了各行业广泛应用的技术之一。
2021-02-19受测量条件限制,激光扫描技术往往需要架设多站点进行扫描,并对获取的点云数据进行配准以消除不同站点的相对位置变化。目前点云配准算法主要分粗配准和精配准2个阶段。迭代最近点算法(iterativeclosestpoint,ICP)[1]为目前应用最广泛且十分有效的一种精配准方法,但计算成本大、易陷入局部最小值。
2020-12-05高铁列车可当做一种全新的震源类型,其确定的长度和荷载,固定的行驶路线以及近乎匀速运动的特性,使之具有可重复震源的特征,另外,地震仪记录的高铁震源记录,具有宽频带分立谱的特征,为高铁高架桥附近地表结构和物性高精度探测提供了条件.因此,近年来,以高铁列车震源为基础的高铁地震学已成为了一个新的研究热点.
2020-12-04第三极不同地区(如季风区和西风影响区)和不同冰川表面(如老雪区、新雪区、裸冰区)雪冰中高分辨率的碳质组分的观测与分析势在必行,特别是消融期开始后碳质组分在短时间尺度内的浓度、黑碳粒径以及有机碳组成变化对于模型的改进至关重要,可为我们全面认识碳质组分对青藏高原不同地区冰川消融影响的机制提供科学依据。
2020-10-24喜马拉雅地区的冰川是南亚地区许多大江大河的源头,如印度河、恒河和布拉马普特河(上游中国境内河段称雅鲁藏布江)。喜马拉雅高大的山体阻挡来自印度洋和阿拉伯海的湿润水汽,南坡温暖湿润,北坡寒冷干旱,剧烈的海拔梯度差异使得在短距离空间内水热急速变化,为各种动植物生存繁衍提供了场所,成为全球生物多样性最丰富的地区之一。
2020-10-24我要评论
期刊名称:地震地质
期刊人气:1475
主管单位:中国地震局
主办单位:中国地震局地质研究所,中国地震局地质研究所
出版地方:北京
专业分类:地质
国际刊号:0253-4967
国内刊号:11-2192/P
邮发代号:82-809
创刊时间:1979年
发行周期:双月刊
期刊开本:16开
见刊时间:一年半以上
影响因子:0.435
影响因子:1.417
影响因子:0.700
影响因子:0.925
影响因子:1.760
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!