Petrel中体属性及层属性提取
第一步： 体属性提取
1. 2. 3. 4.
在home?window? 3D window。
input?seismic?地震数据体上点右键菜单下选volume attributes。 弹出volume attributes窗口。
在volume attributes窗口下可以一次提一个属性，也可以一次提多个属性。
⑴点击add a new attribute to chain图标添加一行，可以点多次添加多行。 ⑵从category种类列可以选择所提属性属于哪一类，这样可以有效过滤属性；在
attributes列选择所提的属性，realize列选择是否产生实体，一般先不选realize，产生一个虚体出来，看效果好的情况下，再realize成实体属性数据。
⑶Input右边初始数据体已经输入进来，output右边空格处可以不输入名字，计算过程中会自动命名，命名规则是：原始地震数据体名［属性名］，提取属性位于原始地震数据体下面，如下：
1
⑷在方差体上点右键菜单选择insert time slice intersection创建时间切片，勾选方差体下的inline，xline，Z显示到3D window中，如下：
⑸Home?inspector，在窗口中出现inspector工具，鼠标切换成select模式，在任何一个地震剖面上点一下，激活inspector窗口信息，到inspector?attribute parameters下，展开看到属性参数的调整界面，可以实时调整参数，看窗口中属性的变化，直到确认一个好的属性显示效果。
⑹设置地震属性存储的统一路径：参照下面Petrel工区规范设置： SeisAttribute存放说明
2
存放个人crop的实体地震体，提取的实体地震属性体，目录按\\survey\\个人存放 可以在Petrel的File-> system -> system setting 下的seismic setting设置路径
只需要选到<项目组>_petrel\\SeisAttribute，后面跟上%s\\%u，每次做realize的时候，如果“survey \\个人”目录不存在，会自动建立目录。
对于这些数据体，项目结束后，要及时删掉数据体，需要到目录下对数据体进行删除。
⑺生成实体属性数据：在方差属性体上点右键选realize，到弹出界面下点击realize即可。
第二步：层属性提取
5. Petrel地震层面属性简介
地震层面属性是储层描述过程中反应储层岩性物性特征的重要参考，计算结果是地
质层段对应的地震时窗内的属性平均值；注意当时窗厚度为零，属性值则为瞬时值。时窗的定义是准确计算的根本。由于振幅类属性比较好理解，本文着重对时窗的定义和波形类的几个属性通过实例进行介绍，帮助读者加深认识。
首先看看属性分类，目前Petrel版本可计算多达50多个层面属性，以下仅列出地震解释人员常用的若干振幅属性和常被忽视的几个波形属性： ⑴振幅类
Average energy – 平均能量（功率） Average magnitude – 平均绝对值振幅 Average negative amplitude – 平均负振幅
3
Average peak value – 平均波峰值（包括局部峰值） Average trough value – 平均波谷值（包括局部谷值） Extract value – 提取沿层瞬时值 Maximum amplitude – 最大振幅
Maximum magnitude – 最大绝对值振幅 Mean amplitude – 算术平均振幅 Minimum amplitude – 最小振幅 RMS amplitude – 均方根振幅 Sum of Amplitude – 振幅值累加 ⑵波形类
首先了解地震波形的几个术语，波形从波谷过渡到波峰的转折点叫S型零值点，从波峰到波谷的转折点叫Z型零值点，两个零值点之间零线和波形线包含的区域叫Loop，把一个Loop等面积划分，上半部分叫Upper Loop， 下半部分叫Lower Loop
下面要经常提到这些术语：
Loop Kurtosis – Kurtosis在统计学中可对数据分布的直方图包络线进行描述，如数据分布是松散（包络线胖）的其K值就小，如果是紧凑的（包络线瘦）其K值就大。换句话说这个参数反映直方图曲线分布的尖锐度。同理，用于地震波形上它可描述波形的胖瘦情况。 Loop Asymmetry – 波形非对称性描述，把LOOP面积等分上下两部分，值大则说明上半部分面积大，头重脚轻；值小则头轻脚重。
Loop duration – 两个零值线之间的距离，波形胖则值大。 Upper loop area – 上半部分的面积。 Lower loop area – 下半部分的面积。 ⑶其它类
Arc length – 波形弧长，也可以间接看作是波形属性，值大有可能是振幅大，也可能是复波的表现。
Number of zero crossing – 时窗内有几个零值点，间接说明频率。
Standard deviation of amplitude – 振幅值的标准方差，时窗内各样点振幅值分布差异。 Positive to negative ratio – 正负振幅比率，主瓣旁瓣的大小比例。
4
⑷非属性类
Make offset horizons – 按照时窗的定义参数，生成时窗顶底层面，可用于后续显示参考也可以用于定义新的时窗。
Isochron thickness – 两层之间的等时厚度，计算等时间图的好地方（Petrel中也可以用surface calculator 计算等厚或等时间图）
Cipher attributes – 提取层面上下相邻同相轴的波峰和波谷的值，以及对应的时间值，或者上下相邻峰谷与参考层面之间的等时间厚度值 6. 层面属性计算界面及应用实例
⑴时窗上下限的指定
大多数碎硝岩储层单元相对较薄，从四分之波长到一个波长不等，考虑到计算的稳定性一般需要一定数量的样点，所以计算属性所需要的时窗厚度一般至少为半个波长到一个波长，亦即一个波峰或波谷，或二者的组合。如果目的层段的顶解释为波峰，按照具体情况时窗就要向下开半个周期至下面的波谷或紧邻的下一个波峰；如果目的层的顶追踪的是波峰上面的零值线，那么时窗就要往下开到波峰下面的零值线或者再下面一个零值线；见下图. 总而言之，时窗最好有一个完整的波峰或波谷，或者两个组合。
下面这个例子提取均方根（RMS）振幅，沿单层Houston上下各开30毫秒开时窗，左面为纯平行时窗没落实到波谷上，右面在上下30毫秒范围内先寻找最大波谷位置再计算RMS振幅。其结果有一定的差异，见下图红框内。右面的RMS高值连片性明显高于左面的平行时窗结果。
5
Petrel中体属性与层属性提取操作流程 - 图文
