strata操作
- 格式：doc
- 大小：841.00 KB
- 文档页数：23
一、 反演简介
地震道的褶积模型如下：
地震道 = 子波 * 反射系数 + 噪声 上式中 * 代表褶积 上式成立的假设条件是：
⏹ 地震道是零偏移距的叠后数据 ⏹ 没有多次波 ⏹ 没有A VO 效应
⏹ 噪音是随机的白噪音, 与地震道不相关, 即没有相关噪音 ⏹ 子波是常量, 没有时变性
⏹ 地震数据已经做过偏移, 每个地震道只与其下方的反射系数序列有
关
第i 个界面的反射系数定义如下：
二、 软件操作
1、 建数据库（database ）
正演建模
输入:
处理:
输出:
反演建模
1111-i i i i i i i
i i V V r V V ρρρρ++++=+density of layer i velocity of layer i
i i V ρ==
设定目录（给的数据库名字不要给后缀）
点击OK，完成。
2、加载地震数据（只有3D数据）
点击右下角
新建地震工区，若工区已建好的则选第二或第三项）
给地震工区名（给地震工区名时，也不要加后缀）
单击OK，出现以下窗口：
点Open seismic from file找到地震数据存放文件：
NEXT》
现有数据起始位一般是INLINE：9 XLINE：21 X：73 Y：77
若是不清楚的，可单击上图detail specification…按扭查看：
此处要注意后面三项的设置，一般分别为：stacked data(4)\meters\根据实际情况。
点NEXT》继续查看，查看完毕进行参数设置。
设置完毕后，点击NEXT》，将会自动扫描出以下结果：检查是否正确，确定后单击OK。
地震数据加载完毕。
3、井数据加载
文件格式栏，现在有的数据选有底色项。
OK。
最后一项填如深度值所在列数。
OK。
将井名填好NEXT》
填写井坐标NEXT》
点击display file…弹出右边信息，找到需要加载的数据（如DEN 等），DEN横向顺序为2，在左边窗口中找到对应行第二行，选择对应曲线类型名称density,并在第二格将名称改为习惯简写DEN，并将usage栏打勾。
以此类推，将所有需要的曲线设置完毕，点OK，井曲线加载完毕。
4、井分层加载
可通过文件加载(tops from file…)和手动加载（tops from table…）
OK。
单击OK。
完成加载。
下面是手工加载：
选中要加载的井，手工敲入有底色表格中的内容。
文件加载有适用于批量加井，手工主要用于个别井添加情况。
4、加载地震解释层位及层位拾取加载选择inport horizon…
拾取选择Pick horizon…
若层位文件如上排列，则对应以下参数设置：inline 和xline分别存放在1和2列，从第四列开始为有效数据，地层的顶底数据分别放在3、4列。
设置好参数，OK。
地震解释层位加载完毕。
若是没有解释好的地震分层，则必须自己拾取，选择Pick horizon…，鼠标左键拾取。
5、合成记录，井震对比
单击edit,弹出以下窗口
单击弹出下面对话框：选择地震道，一般选默认，井旁第一道。
OK，
首先由地震统计得一子波：
设置选用地震数据范围（区域和时窗）：
NEXT》
OK。
得到如下子波我wave1：
将子波设置为wave1，然后用鼠标对应拾取合成记录（蓝）和井旁地震道（红）相似特征明显点，点stretch按扭拉伸调整，使之尽量相关，测井地质分层（tops）要和地震解释层位（horizons）对应。
点stretch按扭时会弹出如下对话框：
注意，黑色曲线为原始声波时差，红色为拉伸压缩后的曲线，两者不能偏离太远，若太远则改变了原始的时深关系太多，产生人为的错误改变。
从以下窗口可以设置计算时窗，进一步提高相关系数。
然后再根据井资料提取子波。
选井：
设时窗：
用默认参数，生成子波wave2：
再将子波设成wave2进行拉伸调整。
重复以上过程（提取子波，拉伸调整），直到相关系数较高，合成道和地震道很好拟合。
6、模型建立
给模型名字：
选择参与建模的井：
生成模型：
7、反演
反演模块分四个：
从上到下，依次是：基于模型反演、有限带宽反演、稀疏脉冲反演、神经网络反演。
(1)、基于模型反演
设置反演数据体的范围：
设置参数：
重要的参数是：反演选项、最大波阻抗变化、平均块的尺寸、迭代的次数。
该参数控制着如何约束反演过程：
基于模型反演使以下的目标函数达到最小化：
J = weight1 x (T - W*r) + weight2 x (M - H*r)
其中：
T = 地震道
W = 子波
r = 最终的反射系数
M = 初始估算模型的波阻抗值
H = 积分因子, 它与最终的反射系数褶积产生最终的波阻抗
* = 褶积
目标函数由两部分组成：
使第一部分(T - W*r)最小化，使得解与真实的地震道相拟合。
使第二部分(M - H*r)最小化,使用指定的块大小使得解与初始估算波阻抗相拟合。
这两部分通常是相互矛盾的。
这些权weight1 和weight2决定着两部分所占的比例。
随机反演中目标函数如上所示. 权由以下参数决定：
模型约束值即为目标函数中weight2 的值。
令这个值为0 则地震道占主导地位。
令这个值为1 则初始估算模型占主导地位。
因为最终的模型极有可能与初始估算模型偏离较远, 所以这叫做“软约束”，但是随着偏离值变大，也会加大对目标函数的惩罚值。
在约束反演中，第二项完全去掉了。
但是，该算法还是能够在特定的范围内通过最大波阻抗变化值来约束最终的波阻抗值。
因为波阻抗值不能超过固定的边界, 所以被称为“硬约束”。
实际工作中, 基于最大波阻抗变化值比模型约束更有意义, 约束反演优于随机反演。
平均块大小：
这个参数控制着最终结果的分辨率。
按该尺寸将初始估算模型划分为一系列均匀的块体
最终的反演结果可能会改变块体的大小，但块体的个数是不变的。
这意味着一些块体变大, 另一些块体变小, 但是总的平均值保持不变。
采用小的块体(例如2ms)会提高分辨率, 但增加的细节可能来自于初始估算模型。
采用小的块体总能提高输入地震道和最终合成记录的拟合程度。
迭代次数：
因为STRATA通过一系列的迭代达到收敛，那迭代次数就决定着收敛的程度。
实际中，一般迭代3次以后都可以达到收敛。
当然了, 多迭代几次除了影响运行时间外不会有什么不好的影响。
达到收敛所需的迭代次数取决于反演中所用的块体的大小。
细的分块需要更多次数的迭代。
（2）、稀疏脉冲反演
前面部分参数和基于模型反演相同，以下参数决定了每个地震道上所需的脉冲数：
最大脉冲个数：
该参数设置每道所允许的最大脉冲个数。
系统缺省为与窗口中的总样点数是相同的, 也就是说正常情况下, 该参数不让改变。
脉冲检测门槛：
每个脉冲加上去时, 它的振幅要与当前所检测到的脉冲的平均振幅相比较，当该振幅小于平均振幅的特定百分比时, 程序则停止添加该脉冲。
（3）带限反演
带限反演唯一的一个参数是：
约束的高切频率：这个参数是应用到初始模型上的滤波器的高切频率，最终反演结果的低频分量来自于初始模型。
初始模型中的高于这个值以上的所有频率成分都滤掉。
递归反演道中所有低于这个值的频率成分也滤掉。
这两个结果然后加在一起。
8、色标调整
反演完后，可更改色标，使显示效果达到最佳，目标层更明显：
点击左上角眼睛图标，弹出右边窗口，可设置不同不同波阻抗值段的显示颜色。
9、生成切片图件
水平切片图：。
