一维合成地震记录

合集下载

人工合成地震记录作业

人工合成地震记录作业

人工合成地震记录程序设计(一)、人工合成地震记录原理:地震记录上看到的反射波波形是地震子波在地下各反射界面上发生反射时形成的。

反射波的振幅有大有小(决定于界面反射系数的绝对值)、极性有正有负(取决于反射系数的正负)、到达时间有先有后(取决于反射界面的深度)的地震反射子波叠加的结果。

如果地震子波的波形用S (t )表示,地震剖面的反射系数为双程垂直反射时间t 的函数,用R (t )表示,那么反射波地震记录形成的物理过程在数学上就可以用S (t )的R (t )的褶积表示,即某一时刻的反射波地震记录f (t )是:)()()(t R t S t f *=其离散形式为:))(()()(1t m n R t m S t n f M m ∆-⋅∆=∆∑=如果大地为多层介质,在地面记录长度内可接收的反射波地震记录为:))(()()(11t m n R t m S t n f Mm N n ∆-⋅∆=∆∑∑== 式中,n 为合成地震记录的采样序号,n =1,2,3...N ;N 为合成一道地震记录的采样点数;m =1,2,3...M ,为离散子波的采样点数;△t 为采样间隔。

这种褶积模型将地震波的实际传播过程进行了简化:1、在合成地震记录的过程中没有考虑大地的吸收作用,所有薄层的反射波都与地震子波的形式相同,只是振幅和符号不同。

2、假设地震波垂直入射到界面上,并原路径返回。

3、假设地层横向是均匀的,在深度(纵向)方向上假设密度为常数,只是速度发生变化。

4、不考虑地震波在传播过程中的透射损失。

(二)、人工合成地震记录的方法1、 反射系数序列在有速度测井资料的情况下,可以用速度曲线代替波阻抗曲线,计算反射系数序列。

在没有速度资料的情况下,可根据干扰波调查剖面分析的结果设计地质模型。

如设计的地质模型如图a 所示,图中H 为层厚度,V 为层速度,根据下式计算反射系数: 11)(--+-=N N N N N V V V V H R 式中H 为反射界面的深度,N 为反射层序号,随深度变化的反射系数序列如图b 所示。

地震资料解释之合成记录制作

地震资料解释之合成记录制作

How To Do Troubleshooting?
Master the basic knowledge & skills ! Think about your trouble ! Always use your head to think and always use
your hands to do it !
地震记录 合成记录
声波曲线
最终得到一个较好的合成地震记录
五、怎样判别合成记录的好坏
相关系数:通过求取合成记录与地震的相关系数来进行定量的 判断。一般在目的层段内的时窗范围内,相关系数都应在0.6之 上,井过断层时除外。
相关系数0.8
与工区内平均速度对比 对做完合成地震记录的井按一定间隔提取其时间-深度对, 与工区的平均速度曲线相对比,应比较一致。
某工区内井的时深关系与平均速度对比图 (红的为平均速度)
剩余记录
地震剖面
合成记录
剩余记录
剩余记录 越少,合 成记录越 精确。
测井曲线
多元标定技术
地震记录
合成记录
测井曲线
综合利用电性曲线、岩性剖面及钻井地质分层的多种信息,精确标定出 地下某地质体的顶底界面,分析储层的横向变化,提高了标定的准确和 形象性,将其与对该区的整体构造与地质认识相对照,分析其是否符合 该区地质沉积规律,也可对合成地震记录起到检查的作用。
自然电位
多元标定技术在合成记录中的应用
横向对比法
通过任意连井线,将测井曲线(波阻抗、自然电位等)投在地震剖 面上,根据地震标准反射层与测井响应一致原则,使井与井之间的 地震地质标定的关系保持一致,从而检验空间地震地质标定的合理 性和一致性。
内容提要
13 为什么要制作合成地震记录 23 合成地震记录原理 3 制作合成地震记录的关键技术 43 体会与认识

Geoframe地震合成记录制作(绝对原版)

Geoframe地震合成记录制作(绝对原版)

Stretch—sqeeze:局部拉伸收缩。微调。 中建选择拉伸点,左键实现拉伸 Bulkshift:对曲线进行整体漂移。粗略对应 地震层位。左键选好始终线实现漂移 时深曲线校正。前提是:有 checkshot。没有的话,用声波曲 线形成后保存checkshot
时深关系曲线编辑(velocity survey):
4. 子波的选取或提取。
选取的基本思路是:先用默认的雷克子波进行观察,看效果。如果不行,自己提取子 波。参数有频率、相位、波长、极性等,提取方法有自相关(默认)和维纳—莱文森 混相子波提取方法。
常用的极性判别方法:
1、单轨、双轨剖面判断法:在正极性剖面上,正反
射系数界面,如基岩顶面、海底、火成岩顶面等,表现为单轨 强峰;而负反射系数界面,如大套油页岩、煤层顶面,表现为 双轨强峰。在负极性剖面上,特征相反。
3dv—define—borehole appearance—synthetic
点亮display synthetics。 进行参数设置。
注意:点亮了沿井显示,则不
能在3dv—tools—synthetics上 移动(见菜单讲解大全)
Update –close。必要时可在marker 、curve将层位测井曲线投上
波长 极性 相位
5.对synthetic进行调整,使其与实际的地震剖面进行匹配。
点右键,选borehole投井,对井右键选appearance投synthetic和曲线 点右键,选content对子波类型选择
Post— correlation 投相关系 数。右键 Maximun correlation Values查看
2.合成记录(synthetic)制作界面的进入
IEXS---application---synthetics

姚姚-地震数值模拟

姚姚-地震数值模拟

(四)地震数值模拟的发展
最简单的地震数值模拟就是一维合成地震 记录形成。它是利用声测井资料得到反射系数 序列,然后与不同的地震子波进行褶积计算得 到各种一维合成地震记录。其目的主要是要进 行波阻抗反演。
后来发展了射线追踪方法,然后又发展波 动方程数值模拟方法。
目前,全三维弹性波波动方程数值模拟已 经步入了实用阶段,考虑更为复杂介质情况的 全三维波动方程数值模拟也有可能了。
逆断层射线追踪数值模拟一炮的射线路径图 正断层射线追踪数值模拟一炮的射线路径图
(一)射线追踪数值模拟的发展与应用
传统射线追踪数值模拟主要是使用两点射线追踪 的试射法和弯曲法。
目前,射线追踪数值模拟在两个方面得到发展。 一是在计算方法上,发展了适应性更强的程函方程法、 波前重建法、最短路径法等;二是在克服只能得到运 动学特征的问题上,发展了近轴射线追踪、动力学射 线追踪、高斯束射线追踪等方法。
(二)地震数值模拟的应用
地震数值模拟在地震勘探、开发中的应用 范围非常广泛,归纳起来最主要有如下三方面 的应用:
1、在地震波场和地震勘探方法的理论研究 中有重要的作用,利用它可以研究各种复杂地 质条件下的地震波场响应特征和验证新方法的 可靠性。
2、在地震资料采集、处理和解释中均起着 重要的辅助作用(如观测系统的设计、解释结 果的验证、处理方法的使用等)。现在我们面 临着越来越复杂的地质对象,可以说,离开它 地震方法寸步难行。
垂直分量(左)
水平分量(右)
双层介质波动方程数值模拟一炮地面记录
射线追踪数值模拟只能模拟地震场的运动 学特征,但是运算速度快,而且提供的射线图 十分直观,在生产实际中还是很受青睐的,其 方法也比较成熟。但是它存在盲区、焦点等难 以解决的问题,特别是缺少地震波的动力学特 征,需要改进。

地震子波 数字信号实验报告

地震子波 数字信号实验报告
实验分析
根据所学知识对实验结果进行分析;
地震子波由震源激发,在地层中传播,因为在沉积地层中,每层介质的物理性质不相同,从而使得地震波的传播速度也不相同。当地震波传播到两层介质的分界面时,会发生反射,由于每层介质的反射系数不同,所以反射波的能量也不相同,检波器接收到不同时刻的、不同能量的反射波,形成一个地震记录。
subplot(2,2,4),plot(f2,'k'),title('最小相位子波地震记录')
f1(i)=0;
f2(i)=0;
end
fori=1:n+m-1%地震记录长度
forj=1:m%地层深度
ifi-j>0&i-j<=200%满足雷克子波的时间序列长度
a1(j)=r(j)*w1(i-j);
a2(j)=r(j)*w2(i-j);
f1(i)=f1(i)+a1(j);
f2(i)=f2(i)+a2(j);
地震子波波形显示及一维地震合成记录
姓名:杨肖迪学号:050422009040专业:地球信息科学与技术2009级
实验目的
1.认识地震子波(以雷克子波为例),对子波的波形有直观的认识。
2.利用褶积公式合成一维地震记录。
实验步骤
1.雷克子波
(零相位子波)
(最小相位子波)
其中 代表子波的中心频率, 代表子波宽度,随着 的增大,子波能量后移,当 =7时,最小相位子波可视为混合相位子波,这里 = 25 Hz, =4;
end
end
end
subplot(2,2,1),plot(w1,'k'),axis([0,200,-1,1]),title('零相位子波')

VSP测井资料在合成地震记录制作中的应用研究

VSP测井资料在合成地震记录制作中的应用研究

[收稿日期]2019-03-15[基金项目]国家科技重大专项(2016ZX05002-002)[作者简介]张建华(1986 ),男,山西绛县人,中国石化胜利油田分公司勘探开发研究院工程师,硕士,主要从事地球物理方法研究㊂doi :10.3969/j.issn.1673⁃5935.2019.03.001V SP 测井资料在合成地震记录制作中的应用研究张建华(中国石化胜利油田分公司勘探开发研究院,山东东营257015)[摘 要] 合成地震记录的制作分为计算波阻抗曲线㊁对波阻抗曲线进行时深转换㊁计算反射系数曲线㊁将反射系数序列与子波进行褶积4个步骤㊂对波阻抗曲线进行时深转换一般采用方法是对声波测井曲线进行简单的积分,但存在两个问题:①在井的最浅部很难进行声波测井,因为缺少浅部的声波测井数据,导致求得的时深关系存在绝对误差;②声波测井数据存在误差,当对测井数据沿着几千米的深度积分时,误差也在累积㊂通常情况下制作的合成地震记录与实际井旁地震道进行对比时需要进行整体漂移及局部拉伸压缩㊂而VSP 测井资料可以提供精确的时深关系,利用VSP 测井资料对波阻抗曲线进行时深转换,制作的合成地震记录不需要整体漂移,也不需要进行局部拉伸压缩就可以与井旁地震道很好吻合㊂[关键词] 合成地震记录;VSP 测井;声波测井;时深转换[中图分类号]P631 [文献标识码]A [文章编号]1673⁃5935(2019)03⁃0001⁃05 地震解释首先要进行层位标定,确定地震反射层位和地质层位之间的关系㊂目前最常用的层位标定方法是合成地震记录标定㊂合成地震记录制作过程中需要对波阻抗曲线进行时深转换,时深转换一般采用的方法是对声波测井曲线进行简单的积分[1]㊂这样制作的合成地震记录在与井旁地震道进行对比时,需要依据标志层的对比分析进行整体的漂移,但在标志层对准的情况下,经常会存在其他的地震反射同相轴无法对准问题㊂笔者对利用声波测井资料计算时深关系的原理进行深入分析,阐明产生这种现象的原因,并研究利用VSP 测井资料的时深关系对波阻抗曲线进行时深转换,进而制作合成地震记录的方法,并对其效果进行对比分析㊂1 合成地震记录基本原理合成地震记录可以看作是一种一维的地震模型正演,即认为地层界面是水平的,射线在界面上是正入射的,射线路径是垂直的㊂这种近似尽管忽略了地震响应随入射角的变化效应,但在通常情况下是合理的㊂当入射角确实对地震响应的影响较大时,必须使用近炮检距部分地震道集叠加而不是全部叠加,来与合成地震记录进行对比㊂合成地震记录的制作通常可以分计算波阻抗曲线㊁对波阻抗曲线进行时深转换㊁计算反射系数曲线㊁将反射系数序列与子波褶积生成合成地震记录4个步骤:(1)计算波阻抗曲线㊂波阻抗是地层密度和纵波速度的乘积,分别可以由密度测井资料和声波测井资料得到㊂如果缺少密度测井资料,考虑到密度的变化远远小于纵波速度的变化,因此可以近似地假定密度不变,即以速度曲线代替波阻抗曲线来计算反射系数,所产生的误差一般情况下可以忽略㊂也可以利用加德纳(Gardner)公式由声波测井资料近似计算得到密度㊂ρ=0.31v 0.25.(1)式中,ρ为密度,g /cm 3;v 为纵波速度,m /s㊂如果缺少声波测井资料,可以利用法斯特(Faust)公式由电阻率测井资料近似计算得到纵波速度,因为岩层速度和岩层电阻率都随岩石孔隙度增加而变小㊂这个公式的适用范围是深度大于12019年9月中国石油大学胜利学院学报Sep.2019第33卷 第3期Journal of Shengli College China University of Petroleum Vol.33 No.3200~300m,自然电位曲线上没有特殊的峰值,而且地层水的矿化度变化很小㊂v=KH1/6R1/6.(2)式中,H为深度,m;R为电阻率,Ω㊃m;K为一个与岩石性质有关的常数㊂(2)对波阻抗曲线进行时深转换㊂将以深度为函数的波阻抗曲线转化为以双程旅行时为函数的波阻抗曲线㊂这就需要知道井的时深关系㊂通常情况下,通过对声波测井曲线进行简单的积分来得到时深关系[2]㊂t H=∫H0τ(h)d h.(3)式中,t H为深度H对应的旅行时,μs;τ(h)为声波时差曲线,μs/m㊂(3)计算反射系数曲线㊂反射系数是地震波垂直入射时,反射振幅与入射振幅的比值㊂地下第n个地层界面处纵波的零炮检距反射系数由如下公式算出:R n=(A n+1-A n)/(A n+1+A n).(4)式中,R n为地下第n个地层界面的反射系数,无量纲;A n和A n+1分别是第n个界面上下两个地层的波阻抗,kg/(m2㊃s)㊂(4)将反射系数序列与子波进行褶积生成合成地震记录㊂地震子波的选取通常有两种方法,一种方法是使用理论子波,如雷克子波,另一种方法是从地震数据中提取子波㊂f(t)=S(t)*R(t).(5)式中,f(t)为合成地震记录,S(t)为地震子波,R(t)为反射系数㊂2 利用声波测井资料进行时深转换存在问题 在上述合成地震记录制作的第2个步骤中,通常通过对声波测井曲线进行简单积分得到时深关系,进而对波阻抗曲线进行时深转换㊂实际计算时,可以采取一种比较简单的方法实现(3)式的计算:假设某口井的声波测井曲线起始深度为d1,在深度上有d1㊁d2㊁d3㊁ ㊁d n-1㊁d n共计n个采样点,可以认为它们将地层分成了c1㊁c2㊁ ㊁c n-1共计n-1层(图1),把每一层当成是均质的,且认为它们的速度分别为1/τ2㊁1/τ3㊁ ㊁1/τn,则每一层的双程旅行时分别为2(d2-d1)τ2㊁2(d3-d2)τ3㊁ ㊁2(d n-d n-1)τn,这样某一深度点所对应的双程旅行时即为其上各层双程旅行时的总和,即t n=∑n k=22(d k-d k-1)τk.(6)式中,t n为深度d n对应的双程旅行时,μs;d k㊁d k-1为声波测井第k个㊁第k-1个采样点对应的深度,m;τk为深度d k对应的声波时差,μs/m㊂图1 利用声波时差计算时深关系的方法由式(6)即可较为方便得到井点的时深关系,但是它存在以下两方面的问题:(1)式(6)计算的双程旅行时是地下某一深度到测井起始深度d1的双程旅行时,而地震数据的双程旅行时是地下某一深度到地震基准面的双程旅行时㊂如图2所示,要得到地下某一深度到地震基准面的双图2 T0求取示意图程旅行时,需要将地下某一深度到测井起始深度的双程旅行时加上测井起始深度至地震基准面之间的双程旅行时T0㊂T0由两部分组成,一部分是地震基准面至地表的双程旅行时T01,可以由地震基准面与地表的高程差和替换速度求得;另一部分是地表至测井起始深度的双程旅行时T02,由于缺少这一深度段的声波时差数据,T02无法直接求取㊂这样,由于缺少井2第33卷 中国石油大学胜利学院学报 2019年 第3期的最浅部的声波时差测井数据而使得T0难以直接的准确求取,从而无法直接得到基于地震基准面的时深关系[3]㊂所以一般情况下制作的合成地震记录无法直接与井旁地震道进行对比分析㊂对于T0难以直接准确求取这个问题,可以采用标志层对比的方法来解决㊂通过分析井点的速度特征,选取一个明显的速度界面作为标志层,在地震剖面上找到该标志层对应的反射同相轴,地震剖面上该标志层对应的时间值与声波时差计算出的该标志层对应的时间值之差即为T0[4⁃5]㊂这样,根据标志层的对比分析,对合成地震记录进行整体的漂移,就可以得到井点基于地震基准面的时深关系,以上就是利用合成地震记录标定得到基于地震基准面的时深关系及平均速度的原理㊂(2)声波时差曲线本身的误差会使计算的时深关系不准确㊂声波时差曲线受井径㊁地层厚度㊁ 周波跳跃”现象的影响,可能无法准确反映地层的时差值[6],进而由声波时差计算得到的时深关系的精度也受到影响㊂这样,在进行合成地震记录对比时,虽然标志层通过整体漂移实现了准确的对应(图3中反图3摇利用声波时差数据进行时深转换的合成地震记录标定3张建华:VSP测井资料在合成地震记录制作中的应用研究1250~1325ms处的3个强轴),但其他的某些射同相轴还是会对不准(图3中725和875ms处的强轴)㊂这时需要对合成地震记录进行局部的拉伸或压缩,使其他的反射同相轴也能准确对应㊂但这种局部拉伸压缩的方法不能滥用,不能通过这种局部拉伸压缩的方法将合成地震记录上所有的同相轴都与井旁地震道的同相轴一一对准㊂因为合成地震记录和实际地震道中都有可能存在不可靠的同相轴㊂声波时差的误差不只影响时深关系,还影响波阻抗㊁反射系数的计算,从而影响合成地震记录中某些同相轴的可靠性㊂井旁地震反射中也会出现多次波等不能真实反映地层界面的同相轴㊂除此之外,还有一些其他的原因,如地震采集处理过程中的某些因素等[7],使得合成地震记录中的同相轴无法与实际地震道中的同相轴完全一一对应㊂3 利用VSP时深关系进行时深转换如果在合成地震记录制作时,该口井有VSP测井资料,利用VSP测井提供的时深关系代替由声波时差计算的时深关系对波阻抗曲线进行时深转换会取得更好的对应效果㊂VSP测井是将测井检波器放入井中,记录下从井口至井下各深度处直达波的传播时间,通过井口地表海拔与地震基准面的高程差和替换速度即可换算出地震基准面至井下各深度处的双程旅行时,不存在像声波测井那样浅部无法获取数据的问题,又由于时间值是直接读得的,因此时深关系误差较小㊂因此,VSP测井是地震与测井信息之间建立联系的最佳方法㊂图4是利用VSP 时深关系进行图4 利用VSP时深关系进行时深转换的合成地震记录4第33卷 中国石油大学胜利学院学报 2019年 第3期时深转换得到的合成地震记录,与图3相比,首先不需要进行整体的漂移,1250~1325ms 处的标志层(3个强轴)对应的比较准确;其次,不需要进行局部的拉伸压缩,其他的反射同相轴(725ms 和875ms处的强轴)就有较好的对应关系㊂对比利用声波时差计算的时深关系和VSP 时深关系也可以看出,在1250~1325ms 标志层处两者较为吻合,越向浅部两者的差别越大㊂图5 声波时差计算的时深关系与VSP 时深关系对比4 结 论在合成地震记录的制作过程中需要将波阻抗曲线由深度域转换到时间域㊂通常利用声波时差测井数据计算时深关系来进行这种时深转换㊂但由于井的浅部缺少声波时差数据,无法直接得到基于地震基准面的时深关系,同时声波时差曲线本身存在误差,会使计算的时深关系不准确㊂而由VSP 测井资料得到的时深关系不存在像声波测井那样浅部无法获取数据的问题,又由于VSP 测井测量的井口至井下各深度处直达波的传播时间值是直接读得的,时深关系误差较小㊂因此,利用VSP 时深关系对波阻抗曲线进行时深转换制作的合成地震记录不需要进行整体的漂移,也不需要进行局部的拉伸和压缩,便可以实现与井旁地震道的标志层和大部分反射同相轴较好的对应㊂[参考文献][1] 程远锋,王延光.三维地震解释[M].北京:石油工业出版社,2013:46⁃51.[2] 陆基孟.地震勘探原理(上)[M].东营:石油大学出版社,1993:207⁃209.[3] 张建华,周进,周川,等.利用声波测井资料求取平均速度方法研究[J].低渗透油气田,2015(1):74⁃77.[4] 浦世照,李世宏,牟中海.地震地质综合层位标定方法及其应用[J].石油地质与工程,2007,21(4):27⁃29.[5] 张永华,陈萍,赵雨晴,等.基于合成记录的综合层位标定技术[J].石油地球物理勘探,2004,38(1):92⁃96.[6] 丁次乾.矿场地球物理[M].东营:石油大学出版社,2002:90⁃93.[7] 郝涛,朱伟.准噶尔盆地腹部高精度合成地震记录制作技术研究分析[J].天然气勘探与开发,2005,28(4):16⁃19.[责任编辑] 董大伟5张建华:VSP 测井资料在合成地震记录制作中的应用研究。

地球物理正反演理论(正演部分)

地球物理正反演理论(正演部分)

正演理论方法
正演理论方法
•地震波场正演数值计算
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
正演理论方法
•算例演示’s Principle
Diffractions
正演理论方法
Huygen’s Principle
Diffractions
正演理论方法
•地震波与岩石物性
•岩石性质与地震波速度 •与岩石性质关系
正演理论方法
•地震波与岩石物性
•岩石性质与地震波速度 •与孔隙度关系
正演理论方法
Diffractions
正演理论方法
Huygen’s Principle
Diffractions
正演理论方法
Huygen’s Principle
Diffractions
正演理论方法
Huygen’s Principle
Diffractions
正演理论方法
Huygen’s Principle
Diffractions
•速度和密度资料的获取 •地震子波的选取
正演理论方法
正演理论方法
•一维模型计算(人工合成地震记录)
•计算合成地震记录
正演理论方法
正演理论方法
正演理论方法
•一维模型计算(人工合成地震记录)
•考虑透过系数合成地震记录(自学) •考虑多次波和透过系数合成地震记录(自学)
正演理论方法
•波动方程克希霍夫积分解

一维合成地震记录

一维合成地震记录

2700
平 均 2600 速 度 2500

) 2400 2300 2500
3000
方案(a)的 方案(b)的
3500 深度(m)
3平均速度分布图
2正确解释方案
4.Q补偿剖面上的合成 地震记录解释 27
解决特殊地质问题
28
S1井气砂层改为含水砂岩和泥岩时 地震响应及与地震剖面的比较
29
VSP资料应用
7
时深曲线
地质 分层
声波曲线
反射系数 地震子波 合成记录 地震子波 合成记录 地震剖面
正极性
负极性
自然电位
地质 分层
8
单井标定—确定井中测量的地层与地震反射的对应关系
官11
Es1
T2
T3
T4
T6 T6’
9
具体制作和使用中需要考虑的因素
(1)制作前
目的,合成地震记录要作什么用,目的层范围
有所估计,有所预计
2、合成地震记录
1
合成地震记录史话:
始于褶积模型的概念,Norman Ricker 1940年, 1953年的经典著作中雏形。褶积模型和CVL的出 现,1950年Peterson用声波测井曲线和一套光电模 型装置首次实现了合成地震记录,示波仪显示, 令人鼓舞。其后,随着数字革命的开始,进行数 字运算。
了解地震剖面的类型,2D,3D,叠加或偏移,
波场类型,频率,井位,反射面貌,断层,要有
良好的反射,剖面质量较好,可能的子波估计,
排列长度,处理流程
测井曲线,AC,SP,井径,电阻率,感应等,
仪器类型,所关心的层位有没有界面,测井曲线
的编辑、取值、异常剔除、异常标记
处理模块
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

车407VSPLOG与地震剖面对比图
微电极
0
声波时差
Μs/M 200
微电极
0 2 4 ΩM 400
声波时差
Μs/M 300 0
微电极
2 4 ΩM
声波时差
400 Μs/M 300
2ΩM 300
拐点
(a)
(b)
图5-13 不同厚度地层的声波速度测井曲线
(c)
Y120井漂移曲线
声波测井曲线
漂移前后合成记录与地震剖面的对比
·声波测井曲线的环境校正 ----泥 浆 浸 泡 对 泥 岩 测 井 速 度 的 影响
3个关键参数: 声波 密度 子波
声波测井曲线的漂移校正 在 输 入 计 算 机 之 前 ,首 先 参 考 其 他 测 井曲线对声波曲线进行编辑解释,以消 除明显的畸变影响或对可能的畸变范围 有所预见。如果有检查炮资料,用检查 炮方法来调整声波测井曲线是一种好方 法。
改进的措施: (1)测井曲线的调整 (2)采样 (3)褶积模型
S1井气砂层改为含水砂岩和泥岩时 地震响应及与地震剖面的比较
VSP资料应用
垂直地震剖面法基本概念
垂直地震剖面是相对于地面地震剖面而言的,它 是在井中观测地震波场,将井下检波仪置于深井中不 同深度记录地面震源所产生的地震信号。 常规地震资料的纵向刻度是以时间(秒或毫秒) 为单位的,即说是时间域的资料,要用于指导勘探 开发得转换到空间域,即需要进行时间-深度转换, 这也常常是解释中的基本问题和难题,而由于VSP 方法是在井下测量,其观测资料中同时具有时间和 深度信息,因此在地震解释中VSP资料有独特的作 用,可用VSP平均速度和VSOLOG进行地质层位的 识别和标定。另外利用VSP对井孔附近的地层和构 造的细节变化进行研究。
RT
AC
?
尖灭砂体在地震记录和合成记录上的响应示意图
~ ~
~ ~
~ ~
地震记录
合成地震记录
改进合成地震纪录的质量 从以上层位匹配流程中可以看出, 声波 合成记录与地震剖面来源不同, 在实际应用 中往往会出现两者难以较好的匹配的现象, 这就需要对其中之一反复修改。 在地震资料 出站后,仅仅为了标定一个层位,不可能再 进行资料的再处理,要使两者较好的匹配, 只有从合成记录的制作过程出发, 对流程中 每一个需要注意的焦点环节问题反复考虑, 合理修改。
初步解释结果(陷阱)
正确的解释方案
Y50-01井合成地震记录 解释分析
2700
平 2600 均 速 度 2500

方案(a)的 方案(b)的
m/s
) 2400 2300 2500
3000
3500
深度(m)
1初步解释方案(陷阱)
3平均速度分布图
2正确解释方案
4.Q补偿剖面上的合成 地震记录解释
解决特殊地质问题
实际 上 合 成 记 录 制作 流 程 是 一 个 相 当 复杂,每个参数的选择和提取都需要仔细斟 酌,制作过程中既要考虑到测井资料的真实 性,同时也要分析地震资料在处理过程中所 用的处理方法,只有这样,两者才能更好地 匹配。任何一个参数的失真,都会严重影响 合成记录制作的质量。
精细合成地震记录:参考其他测井信息对声波测井曲 线进行层状模型解释;用检查炮数据调整声波测井曲线; 选择与地震剖面相适应的褶积模型(包括反射系数和地震 子波两个方面)。 合成地震记录的检验分析:用平均速度分析方法和层 序分析方法来检验合成地震记录标定的正确与否。平均 速度方法是先作出不同解释方案的平均速度分布图,然 后与已知的平均速度进行比较识别陷阱。层序波形分析 法是先确定测区内有代表性、地质含义清楚的反射界面 的波形特征,在解释合成地震记录时以该界面的波形为 对比依据。
雷克脉冲
e
t
2 2
sin( 2 f p t

2
)
fb=1.12fp
N=4
带宽约1.5OCT
辛克子波
2 f 4 sin( 2 f 4 t ) 2 f 4 t

2 f 1 sin( 2 f 1 t ) 2 f 1t
fb=0.7f4
带宽任调
CGG
2 ( 6 t f p 0 .5 ) e
子波的选取
理论子波
子 波
时变子波
统计子波
井旁地震道提取子波
原始子波
常见的几种零相位子波
名称 经典雷克子波 解析式 波形 频谱形状 与量频关系 备注 Ricker(1953) 带宽约1.5OCT
1 2 (
fpt )
2
e
f pt )
2
fp
fb=1.28fp
fb=1.20fp
N=3
N 2 f p /
如果考虑密度变化,其反射系数 k = (4198× 2.21- 4297× 2.53) / (4198× 2.21+4297× 2.53) =- 0.08 密度对反射系数影响率: ( -0.08- -0.012)/-0.08=85%
由此可见,密度对合成记录的影 响程度非常显著。 据此,我们可以在灰岩和砂岩速 度接近时,利用二者的密度差异进行 层位的区分与标定。
18
下 泥 岩 层 段 : 泥 岩 厚 度 93m(2847 - 2958m), 平 均 测 井 速 度 3000m/s, 岩 芯 测 试 平 均 速 度 4050m/s, △ T= △ H /△ V = 93/3000-93/4050=0.008s =8ms 上 泥 岩 层 段 : 泥 岩 厚 度 132m(2670 - 2820m), 平 均 测 井 速 度 2960m/s, 岩 芯 测 试 平 均 速 度 4050m/s, △ T= △ H /△ V =132/2960-132/4050-2960 =0.012s =12ms
地震剖面给出的是时间 !
时变子波 合成 记录
地震剖面
GM
SP
RT
AC
测井资料给出了地层埋深和层速度!
合成地震记录标定
合成地震记录制造流程
大古67井深度域合成地震记录
首先对时 间域的子 波进行时 深转换, 以得到深 度域的子 波,然后 利用与深 度域的反 射系数褶 积,求得 合成记录。
Ed底 Es3底 Ek底
ห้องสมุดไป่ตู้
时深曲线 地质 分层
声波曲线
反射系数
地震子波 合成记录 地震子波 合成记录 地震剖面 正极性 负极性
自然电位 地质 分层
单井标定—确定井中测量的地层与地震反射的对应关系
官11
Es1
T2 T3
T4
T6 T6’
具体制作和使用中需要考虑的因素
(1)制作前 目的,合成地震记录要作什么用,目的层范围 有所估计,有所预计 了解地震剖面的类型,2D,3D,叠加或偏移, 波场类型,频率,井位,反射面貌,断层,要有 良好的反射,剖面质量较好,可能的子波估计, 排列长度,处理流程 测井曲线,AC,SP,井径,电阻率,感应等, 仪器类型,所关心的层位有没有界面,测井曲线 的编辑、取值、异常剔除、异常标记 处理模块 前人的经验
如果考虑密度变化,其反射系数 k = (5159× 2.61- 3747× 2.30) / (5159× 2.61+3747× 2.30) = 0.22 密度对反射系数影响率: (0.22- 0.16)/0.22=27%
史 111 井 在 3467m 存 在 一 个 波 阻 抗 界 面,上为泥灰岩,下为细砂岩,岩芯实验 测 定 物 理 参 数 为 : 泥 灰 岩 速 度 4297m/s , 密 度 2.53; 细 砂 岩 速 度 4198m/s , 密 度 2.21。 如果不考虑密度变化,其反射系数 k = (4198ρ - 4297ρ ) / (4198ρ +4297ρ ) = - 0.012
上为 泥 岩 , 下为 砂 岩 , 岩 芯 实 验 测 定 物 理 参 数 为 : 泥 岩 速 度 3747m/s , 密 度 2.30; 砂 岩 速 度 5159m/s , 密 度 2.61。 如 果 不 考 虑 密 度 变 化 , ρ 为 定 数 , 其 反射系数 k = (5159ρ - 3747ρ ) / (5159ρ +3747 ρ ) = 0.16
校正后
密度对总反射系数的影响 在合成地震纪录的制作过程中往往不太 注意密度对反射系数的影响,只是由声波转 换 而 成 , 实 际 上 反 射 系 数 由 两 者 共 同 决 定 的。在某些情况下,密度的影响是显著的。
史南地区史 育层段,于
105 井 在 沙 三 中 浊 积 岩 发
3269m 存 在 一 个 波 阻 抗 界 面 ,
2
2
6 t
2
f
2 p
fb=fp
与经典Ricker极为相 似,只是显频等于主 频。
统计法提取子波
极性
合成地震记录的应用陷阱分析 合成地震记录一般都不是从地表作起 的 , 其 本 身 也 就 很 难 与 地 震 剖 面 完 全 吻 合,且在使用时在大致的层段用相位对比 的方法,选择较为吻合的方案,这样就可 能使解释落入陷阱。可用平均分析方法和 层序或突变性界面的波形对比方法识别陷 阱。
合成地震记录与地震剖面相关性因素分析
• 与地震剖面有关的: • 野外方法 • 处理过程,叠加、偏 移、滤波、反褶积等 • 与合成记录有关的: • 井的,泥浆,井斜 • 测井的 • 褶积模型,反射系数、 地震子波 • 两者对比: • 频率一速度弥散效应 • 传播路径一各向异性 效应 • 分辨率CDP叠加与 AVO 效应
(2)在制作过程中 采样率 初始 T0 V0 B 褶积模型 算子选择,算子长度,频率,类型,相位 放大显示 波形分析 极性
(3)制作后--怎样用好合成地震记录
用途 评价 陷阱分析 波形分析
相关文档
最新文档