地震属性计算公式.docx

1．层序类层属性大于门槛值的百分比(Percent Greater than Threshold)该类属性用于分析储层内的同相轴，如由很高的值集中于数据引起的振幅异常。计算公式为：该属性主要用于分析地层的延伸，海进和海退垂直序列层序会在高振幅砂岩面和低振幅页岩面之间产生。通过该属性，可以确定这些垂直变化和绘出横向变化的范围图。同样，它可以帮助区分出整合基底(高振幅)、丘状起伏基底(较低振幅)和杂乱反射基底(低振幅)之间的不同。小于门槛值的百分比(Percent Less than Threshold)该类属性用于分析储层内的同相轴，如由很低的值集中于数据引起的振幅异常。计算公式为：该属性主要用于地层走向方面的研究。在特定的第三纪盆地内，三角洲层序是从富含砂，高均方根振幅，到富含页岩前三角洲或深海平原里面的低振幅来划分的。这些油页岩比率的变化通过看图中的小于门槛值的百分比就可以很容易确定。同样，它可以帮助区分出整合基底(高振幅)、丘状起伏基底(较低振幅)和杂乱反射基底(低振幅)之间的不同。吸收系数(Absorption Coefficient)吸收系数是用来表示地震波振幅A沿传播距离的衰减系数，即：其中，为波的传播距离，为起始振幅。吸收系数与地震波速度之间存在明显的对比关系，高速的岩石，吸收系数低；低速的岩石，吸收系数高。吸收系数如同速度一样，频数异常现象较弱。2.地震纹理属性(1)地震数据灰度化设三维地震数据中的一个地震体素点为X(x,y,z,a)，其中x，y，z代表了线号、道号和时间；a代表振幅值。设定某个灰度阶数G，用将地震数据a转化为灰度数据g，式中Amin和Amax所有地震数据中的最小值和最大值；灰度阶数G决定了地震数据的粗化程度。(2)灰度共生矩阵灰度共生矩阵定义为三维地震数据中的范围内，两个距离为、方向为的地震体素在地震数据中出现的概率P，记为由于矩阵只能表示一个面，所以每次计算分别取x，y，z中的两个方向，地震体素简记为X(x,y,g)，i和j代表某位置的灰度值。N表示的范围内地震体素总数。应取、、和四个方向，分别计算出纹理属性值后，再求四个值的平均作为输出的纹理属性值。考虑到地震数据中相邻道的连续性问题，取。归一化后阶的灰度共生矩阵M有如下形式其中：当，时，；当，时，。考虑到地震数据中不同地质目标的不同尺度问题，数据分析窗对地震纹理的影响也很大，也是一个值得讨论的参数。图2演示了从图1左上角取出大小的4阶灰度数据后计算灰度共生矩阵的结果，可以发现矩阵M中的值在对角线集中的现象。图1地震道振幅数据灰度化示意图图2不同方向的灰度共生矩阵（未归一化）示意图(a)；(b)；(c)；(d)(3)4种灰度共生矩阵的统计量特征纹理熵：纹理均质性：纹理对比度：纹理能量：3.曲率属性计算曲率离不开求导，因此求导的方法决定了曲率的计算方法，这里主要介绍差分法计算的倾角(Dip Angle)和最大正曲率(Most Positive Curvature)。为计算如式(1)所表达的曲面上某一点的曲率，用其周围网格点的值对局部二次曲面进行最小二乘法拟合。若采用 3×3 网格单元作逼近（图 3），由一阶和二阶导数的定义，得式中各系数（式（2）~式（6））。图3用于差分计算的网格单元 (1) (2) (3) (4) (5) (6)式中：是图 5 中所示层面网格结点 1至9 的值；、是网格结点间的距离。在此基础上，可计算各种曲率属性。(1)倾角(2)最大正曲率4.频谱类层属性有效带宽(Effective Bandwidth)有效带宽是自相关函数的极值与双边自相关函数值的总和之比。计算公式为：式中，是时窗内样点数据的零延迟自相关，代表道两旁所有自相关的总和，是采样间隔。该属性反映地层的横向相似性，宽的带宽反映非均匀，复杂地层；窄的带宽反映简单或圆滑的反射特征，可认为是均质的地层。带宽能帮助我们在数据体中识别噪声区，有噪声的数据体比没有噪声的数据体有很明显宽的带宽。同时，它可用于指示沉积环境的变化。提取时窗应大于2倍主波长，一般为50~200ms.波形长度(Arc Length)波形长度是指时窗内地震波振动曲线的长度。计算公式为：式中，代表时窗内样点个数，代表采样间隔，代表时窗内第i个样点的振幅值。该属性反映高振幅高频率和高振幅低频率之间与低振幅低频率和低振幅高频率之间的差别。比如页岩和砂岩的界面，一般有一些突变和高阻抗的反差，波形长度就用于页岩层序和含砂量较高的层序之间的识别，带宽越小，它就越接近总绝对值振幅。它类似于反映反射层的均质性。提取时应选取较大的时窗，一般为50~100ms.平均零交叉点数(Average Zero Crossing)平均零交叉点数也称为平均零交叉点频率，计算公式为：式中，代表时窗内过零点的样点个数，、代表