当前位置：文档之家 › 合成地震记录

合成地震记录

  • 格式：docx
  • 大小：14.20 KB
  • 文档页数：5

下载文档原格式

  / 5

合集下载

合成地震记录制作

合成地震记录制作

我们知道计算合成地震记录的基本原理是,合成地震记录=子波与反射系数的褶积所以需要子波和反射系数.但是用于计算的数据一般是深度域的,要转换到时间域来必须有时深关系.所以.需要的数据:时间/深度关系数据:checkshot或者DT,用于计算反射系数的数据,一般是DT和密度(RHO B).基本步骤:1, 加载数据:如果是斜井的话,加载井斜,计算出SSTVD,设置成Prefered DS(deviation survey);如果有来自VSP或者其他可信渠道的时深关系的话加载进来,叫checkshot,就是时间,深度关系对,用于提供时深关系;加载DT,RHOB曲线;2,数据质量检查:查看checkshot数据覆盖范围,和品质;查看DT,RHOB曲线的品质,如果不好需要用well-edit或者synthetics里带的一些功能进行编辑.DT,RHOB曲线应该是做过Depth match,需要拼接的话是splice好的.3,制作合成地震记录:点击Post,依次选择时深关系,声波曲线,密度曲线(如果没有密度曲线或者品质不好也可以使用经验公式来代替),声波阻抗,反射系数,子波,合成地震记录,地震数据.软件完全是根据原理走的,如果时深关系没选,后续工作无法开展,如果没有DT,密度,就无法生成声波阻抗和反射系数...软件自带有Ricker30经验子波.如果效果不好可以自己提取子波,也可以使用时变子波.4,对比合成地震记录和井旁道实际地震记录,通过bulkshift或者拉伸压缩来调整时深关系.有时需要用c heckshot来校正DT.一般可能先使用Ricker30子波试一下,看看大致情况,如果效果不好,再尝试提取子波.这是一个反复实验的过程.合成地震记录的品质和制作的数据来源的品质有关,对比的好坏和实际地震数据的品质也有关系.总是实际情况总是复杂的.。

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

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

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_合成反射地震记录 [兼容模式]

实验1_合成反射地震记录 [兼容模式]

地震反射
地震反射的振幅谱特征
复杂地质模型的反射地震记地震记录，N=256
振幅谱
相位谱
雷克子波（Ricker Wavelet）
w t 1 2 f t e
2 2 2 m
2 2 2 f m t
其中：fm表示雷克子波的主频。
Ricker子波，fm=60Hz，N=31， ts=1ms
简单地质模型的反射地震记 录合成及谱分析：
含双反射界面的地质模型 Time Thickness=18ms
合成反射地震记录及谱分析基础实验
在地震勘探中，接收到的地震信号可看成是发射的震 源子波受大地滤波后得到的结果，因此，反射地震信号 （记录）可表示为地震反射系数与震源子波的卷积，即
s t r t w t
其中：s(t)为反射地震信号或地震记录，r(t)表示地震 反射系数，w(t)为震源子波，t表示时间，根据傅里叶变换 的卷积定理，在频率域，地震信号的频谱可表示为地震反 射系数与震源子波两者频谱的乘积，即
S R W
其中：S(ω) 、R(ω) 和W(ω)分别表示为地震记录、地震反射系数和震源子 波的频谱。
对于地震波垂直入射的情况，不产生转换波，地震纵 波的反射系数可简化表示为
r
V p 2 2 V p 1 1 V p 2 2 V p1 1
其中，V,ρ分别表示地层的纵波速度和密度，下标1，2分 别表示相邻的上一地层和下一地层。可简化表示为 在人工合成地震记录时，常利 用雷克子波作为震源子波，如下 式

褶积模型合成地震记录matlab

褶积模型合成地震记录matlab

褶积模型合成地震记录matlab
在地震学中，褶积模型（Convolution Model）是一种常见的地
震记录合成方法，它可以用于模拟地震波在地下介质中的传播和记
录地震仪上的地震信号。

使用Matlab进行褶积模型的合成可以通过
以下步骤实现：
1. 定义地震波模型，首先，你需要定义地震波的源项，包括地
震波的类型（P波或S波）、震源机制、震源位置等。

这些参数将
用于生成地震波的初始模型。

2. 定义介质模型，接下来，你需要定义地下介质的模型，包括
介质的速度结构、密度结构等。

这些参数将影响地震波在介质中的
传播。

3. 生成地震记录，使用定义的地震波模型和介质模型，利用褶
积运算来合成地震记录。

在Matlab中，你可以使用卷积函数（conv）来进行褶积运算，将地震波模型与介质响应进行卷积，得到地震记录。

4. 可视化和分析，最后，你可以对合成的地震记录进行可视化
和分析，以便进一步研究地震波在地下介质中的传播特性。

需要注意的是，在实际操作中，还需要考虑到地震波的衰减效应、边界条件等因素，以及对地震记录进行滤波、去噪等后处理步骤。

同时，合成地震记录的精确度也取决于所用的地震波模型和介
质模型的准确性。

总之，在Matlab中进行褶积模型的合成需要对地震波和地下介
质有深入的理解，同时结合Matlab强大的数学计算和可视化功能，
可以实现高效、准确的地震记录合成。

希望这些信息能够帮助到你。

提高合成地震记录精度方法探究

提高合成地震记录精度方法探究

提高合成地震记录精度方法探究合成地震记录(简称合成记录)是连接地震资料与测井资料的纽带,是联合高分辨率的测井信息与区域性的地震信息的桥梁,其精度直接影响到地质层位的准确标定。

目前油气勘探工作越来越向隐蔽性油气藏发展,目标尺度越来越小,对合成地震记录提出了更高的要求[1]。

由于合成记录的制作过程中存在诸多的制约因素,合成地震记录与实际地震剖面往往并不能完全一致。

本文将着重从原理上分析这些制约因素,并总结相应的技术对策。

1、合成地震记录的制作原理合成记录的制作是一个简化的一维正演的过程,合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果。

F(t) =S(t)×R(t)合成地震记录制作的一般流程是:由速度和密度测井曲线计算得到反射系数,将反射系数与提取的地震子波进行褶积得到初始合成地震记录。

根据较精确的速度场对初始合成地震记录进行校正,再与井旁地震道匹配调整,得到最终合成地震记录。

(图1)图1 东部油田某井合成地震记录2、合成记录的制作过程中的质量控制合成地震记录制作主要包括计算反射系数、提取地震子波和匹配调整时深关系三个环节,如何对这三个环节进行有效的质量控制成为决定最终精度的关键,在实际制作过程中依照下列方法进行相应的质量控制[2]。

2.1反射系数由于测井曲线在测量过程中受所使用仪器、施工单位以及环境的影响,所以在制作合成地震记录之前应该对密度和声波曲线进行标准化处理。

测井曲线的标准化处理是一个非常重要的工作,精确的测井曲线可以提高反射系数的准确性和可信度。

整理测井曲线可以从以下几个方面进行: 对测井曲线进行标准化处理是为了消除曲线中的一些压制正常的波峰、波谷的奇异值。

首先计算标准井标准层段的测井曲线平均值,然后计算每口井目的层段测井曲线平均值,求出校正差,再从测井曲线中减去校正量,即完成标准化校正(图2)。

图2 标准化曲线流程图反射系数由速度和密度测井数据计算得到,反映上、下两层介质的波阻抗差。

合成地震记录

应用合成地震记录来标定地震层位是地震资料解释中非常重要的手段，也是将地震资料与测井资料相结合的一条纽带。

它最终使抽象的地震数据与实际的地质模型连接起来，为地震资料解释的可靠性提供了依据。

合成记录的精度将直接影响到地震地质层位标定的准确性，因此，提高合成记录的精度就成了地震层位标定的首要问题。

１合成记录的方法原理１．１合成地震记录制作的一般方法一般而言，人工合成地震记录，是利用声波和密度测井资料求取一反射系数序列，再将这一反射系数序列与某一子波反褶积得到结果。

Ｓ（ｔ）　＝　Ｒ（ｔ）　＊　Ｗ（ｔ） （１）式中 Ｓ（ｔ）　——　合成地震记录； Ｒ（ｔ）　——　反射系数序列； Ｗ（ｔ）　——　地震子波。

上式表明，合成记录的好坏与反射系数序列的求取和子波的选择有着密切的关系。

反射系数序列的准确性和精确程度又与测井资料（声波、密度）的采集、处理等过程密切相关；子波的选择，则要考虑子波的长度、相位、频率等诸多因素。

在实际工作中，所得到的结果往往不尽人意［１］，主要表现在：（１）　合成地震记录与井旁地震道附近的地震剖面层位不吻合现象较多，或者说同相轴吻合的时窗长度有限；（２）　合成地震记录与井旁地震道附近的地震剖面能量不吻合现象较多，或者说同相轴“胖瘦”程度吻合有限；（３）　合成地震记录与井旁地震道附近的地震剖面存在一定的时移。

其原因主要在于：①子波受地质条件变化的影响，难以给得恰到好处；②深—时转换存在误差；③褶积模型并不能完全准确地反应地震记录；④实际地震记录存在噪声。

１．２实用优化方法１．２．１校正测井数据首先对测井数据进行校正，对反射系数序列进行非均匀采样［２，３］。

１．２．２选择合适的子波（１）子波的类型。

常用的子波有两类，一是典型子波，如Ｒｉｃｈｅｒ、Ｔｒａｐｅｒｉｏｄ子波等；二是提取子波，提取子波一般有维纳—莱文森混相位子波提取法和自相关子波提取法两种［４，５］。

从剖面提取的实际子波制作的合成记录，虽然其合成地震记录层位精细标定应用研究＊洪余刚 陈景山 代宗仰 李凌峰（西南石油学院资源与环境学院，四川省成都市６１０５００）摘 要：通过对合成记录制作的一般方法进行分析，结合研究区实际地质、地震资料，提出合成记录的制作在层位标定中的实用优化方法，强调了子波的提取方法和子波相位引起的偏差。

GeoFrame合成记录

GeoFrame 合成地震记录培训和训练手册 版本 4.2斯伦贝谢科技 斯伦贝谢科技（ 科技（中国） 中国）公司 SIS 部门 2010 年 3 月BreakthroughTeam PerformanceSchlumberger版权声明© 2008 Schlumberger。

保留所有权利。

未经 Schlumberger 公司 SIS 部门（5599San Felipe, Suite 100, Houston, TX 77056-2722）的 事先书面许可，禁止任何人对本手册的任何部分进行复制或保存，如保存于可重新获 取的系统内、转化为其它存储格式、通过电子或机械手段（如复印和录音）。

免责声明本产品的使用通过许可证来控制。

Schlumberger 公司对于本产品，不作任何保证、声 明、暗示或规定，并声明不限于任何商业或特殊用途。

Schlumberger 保留在不另行通 知的情况下，修订本手册中信息的权利。

商标信息除非另有说明，本手册中使用的软件应用标志是斯伦贝谢公司的注册商标。

某些其它 产品和产品名称为其各自公司或组织的商标或注册商标。

目 录第一章 合成地震记录的基本概念 .......................................... .......................................... 1一、介绍 ............................................................................................................................................................ 1 二、为何与如何使用地震合成记录？ ............................................................................................................ 2 三、总结 ............................................................................................................................................................ 2第二章 数据要求、 数据要求、典型工作流程 .......................................... .......................................... 3一、数据要求 .................................................................................................................................................... 3 二、典型工作流程 ............................................................................................................................................ 3 三、练习 1......................................................................................................................................................... 5 四、总结 ............................................................................................................................................................ 6第三章 理解 GeoFrame 中的基准面和高程 中的基准面和高程 ................................... ................................... 7一、概念介绍 .................................................................................................................................................... 7 二、练习 2....................................................................................................................................................... 10 三、总结 .......................................................................................................................................................... 11第四章 准备开始 ....................................................... ....................................................... 12一、介绍 .......................................................................................................................................................... 12 二、主窗口 – 常规特征 .................................................................................................................................. 12 三、Layout（布局） ....................................................................................................................................... 13 四、Panels （图版） ...................................................................................................................................... 14 五、MB1 （鼠标左键）用法 ......................................................................................................................... 14 六、MB3 （鼠标右键）用法 ......................................................................................................................... 15 七、数据依存关系 .......................................................................................................................................... 16 八、练习 3....................................................................................................................................................... 17 九、总结 .......................................................................................................................................................... 17第五章 数据质控和编辑 ................................................. ................................................. 18一、介绍 .......................................................................................................................................................... 18 二、质量控制（Quality Control） ................................................................................................................. 18 三、曲线数据 .................................................................................................................................................. 18 四、声波 .......................................................................................................................................................... 18 五、密度 .......................................................................................................................................................... 19 六、曲线编辑 .................................................................................................................................................. 19 七、深度校正 .................................................................................................................................................. 20 八、拼接 .......................................................................................................................................................... 20 九、平抑尖峰 .................................................................................................................................................. 20 十、曲线方波化 .............................................................................................................................................. 20 十一、Checkshot ............................................................................................................................................. 21十二、井斜（Deviation Surveys） ................................................................................................................. 21 十三、地震数据 .............................................................................................................................................. 21 十四、练习 4................................................................................................................................................... 22 十五、总结 ...................................................................................................................................................... 23第六章 声波阻抗和反射系数序列 声波阻抗和反射系数序列 ......................................... ......................................... 24一、介绍 .......................................................................................................................................................... 24 二、声波阻抗（AI）....................................................................................................................................... 24 三、反射系数（RC）序列 .............................................................................................................................. 24 四、RC Parameters（反射系数参数） .......................................................................................................... 26 五、RC 类型 .................................................................................................................................................... 27 六、RC 采样..................................................................................................................................................... 27 七、Offset RC （偏移距 RC）......................................................................................................................... 27 八、常见问题 .................................................................................................................................................. 28 九、RC Modeling （反射系数模型）............................................................................................................. 28 十、练习 5....................................................................................................................................................... 29 十一、总结 ...................................................................................................................................................... 29第七章时间时间-深度关系 ................................................. ................................................. 30一、介绍 .......................................................................................................................................................... 30 二、Checkshot vs 声波曲线 .......................................................................................................................... 30 三、VSP（垂直地震剖面）相关的一个注释 ................................................................................................ 31 四、编辑时间-深度关系 ................................................................................................................................. 31 五、整体偏移（BulkShift）............................................................................................................................ 31 六、拉伸和压缩（Stretch and Squeeze） ..................................................................................................... 31 七、Checkshot 编辑器 .................................................................................................................................... 32 八、调整校正速度 .......................................................................................................................................... 32 九、保存时间-深度曲线 ................................................................................................................................. 33 十、练习 6....................................................................................................................................................... 33 十一、总结 ...................................................................................................................................................... 33第八章声波曲线校正 .................................................. .................................................. 34一、介绍 .......................................................................................................................................................... 34 二、工作流程 .................................................................................................................................................. 34 三、Velocity Survey （速度校正）窗口 ........................................................................................................ 34 四、偏移曲线 .................................................................................................................................................. 35 五、校正声波曲线 .......................................................................................................................................... 36 六、练习 7....................................................................................................................................................... 37 七、总结 .......................................................................................................................................................... 37第九章 子波提取 ....................................................... ....................................................... 38一、介绍 .......................................................................................................................................................... 38二、使用一个存在的子波 .............................................................................................................................. 38 三、从其它工区输入一个子波 ...................................................................................................................... 39 四、理论子波 .................................................................................................................................................. 39 五、子波提取 .................................................................................................................................................. 40 六、统计提取算法 .......................................................................................................................................... 41 七、确定提取算法 .......................................................................................................................................... 42 八、统计参数 .................................................................................................................................................. 46 九、时变子波 .................................................................................................................................................. 47 十、编辑子波 .................................................................................................................................................. 47 十一、结论 ...................................................................................................................................................... 48 十二、练习 8................................................................................................................................................... 48 十三、练习 9 （IESX 用户参看） ................................................................................................................. 49第十章 相关工具 ....................................................... ....................................................... 51一、介绍 .......................................................................................................................................................... 51 二、准备开始 .................................................................................................................................................. 51 三、相关道显示 .............................................................................................................................................. 51 四、相关偏差 .................................................................................................................................................. 53 五、编辑时间-深度关系 ................................................................................................................................. 53 六、练习 10..................................................................................................................................................... 54 七、总结 .......................................................................................................................................................... 54第十一章 保存和使用合成记录数据 ....................................... ....................................... 55一、介绍 .......................................................................................................................................................... 55 二、结果分析 .................................................................................................................................................. 55 三、保存合成记录 .......................................................................................................................................... 56 四、在 Charisma 中显示 .................................................................................................................................. 57 五、在 IESX 中显示 ......................................................................................................................................... 57 六、练习 11..................................................................................................................................................... 58 七、总结 .......................................................................................................................................................... 58Schlumberger合成地震记录的基本概念第一章一、介绍合成地震记录的基本概念什么是合成地震记录？合成地震记录是一个在井眼位置的地震道模型。

合成地震记录业务流程

合成地震记录业务流程一、准备工作。

咱得先把相关的数据都找齐喽。

比如说，测井数据那是相当重要的。

就像我们找宝藏得有个地图一样，测井数据就是我们合成地震记录的地图。

这里面包括声波测井曲线、密度测井曲线等。

这些数据就像是一个个小零件，缺了哪个都不行。

而且呀，我们还得确保这些数据的准确性，如果数据错了，那就好比做菜的时候盐当成了糖，做出来的东西肯定不对味。

另外呢，我们还需要有一些地质分层信息，这个就像是房子的框架结构，能让我们清楚地知道不同地层的情况，知道在哪个地层该怎么操作。

二、选择合适的子波。

子波就像是合成地震记录的画笔。

有好多不同类型的子波可以选呢。

我们得根据实际的地质情况和研究目的来挑。

如果是比较简单的地层结构，可能选个简单点的子波就够用啦。

但要是地层情况很复杂，就像一个超级复杂的迷宫一样，那我们就得找个功能强大、能适应复杂情况的子波。

这时候就得花点心思去对比不同的子波，看看哪个画出来的“画”（也就是合成的地震记录）最符合我们对这个地下情况的预期。

三、计算反射系数。

这一步就像是在做数学题，不过是很有趣的那种。

我们要根据前面准备好的测井数据，像声波和密度这些，来计算反射系数。

反射系数就像是镜子的反射率一样，它能告诉我们地震波在不同地层界面上反射的情况。

这个计算可不能马虎，要是算错了，那合成出来的地震记录就会像一个歪歪扭扭的积木塔，一点都不稳定也不准确。

我们得仔仔细细地按照公式来算，就像小心翼翼地搭积木一样，一块都不能搭错。

四、合成地震记录。

好啦，前面的工作都做好了，就到了最激动人心的合成地震记录这一步啦。

我们把选好的子波和计算好的反射系数放在一起，就像把颜料和画笔放在一起准备画画一样。

然后通过一些算法，让它们相互作用，就像魔法一样，一个地震记录就慢慢合成出来了。

这时候我们就像一个小魔法师，看着自己的作品一点点呈现出来。

不过呢，这时候还不能掉以轻心，我们还得检查这个合成出来的地震记录是不是合理。

五、验证与调整。

2合成地震记录


如果考虑密度变化，其反射系数 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
2
6 t f p
2
2
fb=fp
与经典Ricker极为相 似，只是显频等于主 频。
统计法提取子波
极性
合成地震记录的应用陷阱分析 合成地震记录一般都不是从地表作起 的，其本身也就很难与地震剖面完全吻 合，且在使用时在大致的层段用相位对比 的方法，选择较为吻合的方案，这样就可 能使解释落入陷阱。可用平均分析方法和 层序或突变性界面的波形对比方法识别陷 阱。
微电极
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井漂移曲线
声波测井曲线
漂移前后合成记录与地震剖面的对比
·声波测井曲线的环境校正 ----泥浆浸泡对泥岩测井速度的 影响
如果考虑密度变化，其反射系数 k ＝(4198×2.21－4297×2.53) / (4198×2.21+4297×2.53) =－0.08 密度对反射系数影响率： ( -0.08－-0.012)/-0.08=85％
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
R(ilayer) = (z2-z1) / (z2+z1);
end
%
%计算各反射界面所对应的时间tlength
%
tlength(1) = 2*h(1)/V(1);
for ilayer = 2:n-1
tlength(ilayer) = tlength(ilayer-1) + 2*h(ilayer)/V(ilayer);
end
reply = input('请输入Ricker子波的频率f和采样间隔dt(Defalt=40 0.004):','s');
if isempty(reply)
f = 40;%子波频率f和采样间隔dt
dt = 0.004;
else
clear a;
a = sscanf(reply,'%f',[2 1]);
deltT(i) = t(i) - t(1);
delay(i) = floor(deltT(i)/dt);
end
%
%形成地震记录SS不明白！！！！
%
for i = 1:ntrace
SS(1:2*nsample,i) = i*1;
for j = delay(i)+1:2*nsample
SS(j,i) = S(j-delay(i))+1*i;
plot(SS,Tss,'b');
hold on;
plot(Q,Qt,'r-');
set(gca,'XAxisLocation','top');
set(gca,'YDir','reverse');
title('地震合成记录');
for i = 1:length(wavelet);
Tw(i) = i*dt;
end
%subplot(2,2,2);
plot(Tw,wavelet);%应该是能画出雷克子波
title('Ricker子波');
%
%褶积后得到单道记录S
%
S = conv(RR,wavelet);
for i = 1:length(S)
Ts(i) = i*dt;
end
%subplot(2,2,3);
plot(S,Ts);
set(gca,'XAxisLocation','top');
set(gca,'YDir','reverse');
title('单道记录');
%
%计算多层反射界面是的平均速度Vav,等效深度h0以及获得偏移距deltX
f = a(1);
dt = a(2);
end
%
%计算各反射界面所对应的采样点数nR
%
nsample = floor(tlength(n-1)/dt);
for ilayer = 1:n-1
nR(ilayer) = floor(tlength(ilayer)/dt);
end
%
%形成反射系数序列RR
%
RR(1:2*nsample) = 0;%？这个地方反射系数的长度应该是nsample/2
a = sscanf(reply,'%f',[1 1]);
deltX = a;
end
%
%计算延迟时间delay完全不知道算的啥意思！！！！！@
%
ntrace = 10;
t(1) = sqrt(deltX^2 + 4*h0^2) / Vav;
deltT(1) = 0;
for i = 2:ntrace
t(i) = sqrt((deltX*i)^2 + 4*h0^2) / Vav;
for ilayer = 1:n-1
RR(nR(ilayer)) = R(ilayer); %只有在有界面的地方反射系数才有值
end
%subplot(2,2,1);
stem(RR);
title('反射系数序列');
%
%形成一个Ricker子波wavelet
%
wavelet = ricker(f,dt);
%
for i = 1:n-1
v(i) = V(i);
end
Vav = sum(h) / sum(h/v);
h0 = sum(h);
reply = input('请输入偏移距deltX(Defalt = 400):','s');
if isempty(reply)
deltX = 400;
else
clear a;
V = a(1,:);
dens = a(2,:);
h = a(3,:);
end
%
%计算反射系数R
%
for ilayer = 1:n-1
z1(ilayer) = V(ilayer) * dens(ilayer);
z2(ilayer) = V(ilayer+1) * dens(ilayer+1);%各层反射系数R
%地震合成记录
%日期：07.07.19
%
clc
clear
reply = input('请输入层数n(Default=5):','s');%层数为n
if isempty(reply)
n = 5;
else
n = sscanf(reply,'%f',[1 1]);
end
reply = input...
('请输入各层速度、密度及层厚(Defaul=[600 1000 1500 2000 2500;1500 1800 2000 2500 3000;500 700 400 300]):','s');
end
end
% for i = 1:ntrace
% S(1:nsample,i) = i*1;
% delay = 10*(i-1);
% for j = delay+1:nsample
% S(j,i) = SS(j-delay)+1*i;
% end
% end
for i = 1:length(SS)
Tss(i) = dt*i;
if isemp00];
dens = [1500 1800 2000 2500 3000];%速度和密度v和den
h = [500 700 400 300];
else
clear a;
a = sscanf(reply,'%f',[3 n]);
end
for i = 1:nsample
if abs(S(i))>0
bg=i;
break;
end
end
for i = 1:ntrace
Qt(i) = dt*(delay(i)+bg);
end
for i = 1:ntrace
Q(i) = SS(delay(i)+bg,i);
end
%subplot(2,2,4);