《地震属性计算公式》由会员分享，可在线阅读，更多相关《地震属性计算公式（12页珍藏版）》请在金锄头文库上搜索。
1、1层序类层属性1）大于门槛值的百分比（Percent Greater than Threshold）该类属性用于分析储层内的同相轴，如由很高的值集中于数据引起的振幅异常。计算公式为：PGT =大于门槛值的样点数总的样点数x 100%该属性主要用于分析地层的延伸，海进和海退垂直序列层序会在高振幅砂岩面和低振幅 页岩面之间产生。通过该属性，可以确定这些垂直变化和绘出横向变化的范围图。同样，它 可以帮助区分出整合基底（高振幅）、丘状起伏基底（较低振幅）和杂乱反射基底（低振幅）之间的 不同。2）小于门槛值的百分比（Percent Less than Threshold）该类属性用于分析储层内的同相轴，
2、如由很低的值集中于数据引起的振幅异常。计算公式为：PLT =小于门槛值的样点数总的样点数x100%该属性主要用于地层走向方面的研究。在特定的第三纪盆地内，三角洲层序是从富含砂， 高均方根振幅，到富含页岩前三角洲或深海平原里面的低振幅来划分的。这些油页岩比率的 变化通过看图中的小于门槛值的百分比就可以很容易确定。同样，它可以帮助区分出整合基 底（高振幅）、丘状起伏基底（较低振幅）和杂乱反射基底（低振幅）之间的不同。3）吸收系数（Absorption Coefficient）吸收系数是用来表示地震波振幅A沿传播距离的衰减系数，即：A = A0 e-g其中，x为波的传播距离，A0为起始振幅。吸收系数
3、与地震波速度之间存在明显的对 比关系，高速的岩石，吸收系数低；低速的岩石，吸收系数高。吸收系数如同速度一样，频 数异常现象较弱。2. 地震纹理属性（1）地震数据灰度化设三维地震数据中的一个地震体素点为X（x,y,z,a），其中x，y，z代表了线号、道号和时间；a代表振幅值。设定某个灰度阶数G,用a-Ag = minx GA - A maxmin将地震数据a转化为灰度数据g,式中A .和A 所有地震数据中的最小值和最大值；灰minmax度阶数G决定了地震数据的粗化程度。(2) 灰度共生矩阵灰度共生矩阵定义为三维地震数据中W xW的范围内，两个距离为5二：Ax2 +Ay2、xy方向为0的地震体素在
4、地震数据中出现的概率p,记为P(i, j,5,0) = (x, y),(x + Ax, y + Ay川f (x, y) = i, f (x + Ax, y + Ay) = jx, y = 0,1,., N -1由于矩阵只能表示一个面，所以每次计算分别取X，y，z中的两个方向，地震体素简记 为X(x,y,g),i和j代表某位置的灰度值。N表示W xW的范围内地震体素总数。0应取0、xy45、90和135。四个方向，分别计算出纹理属性值后，再求四个值的平均作为输出的纹理属性值。考虑到地震数据中相邻道的连续性问题，取5 = 1。归一化后Gx G阶的灰度共生矩阵M有如下形式PPP0,00,10,G-1
5、PPP1,01,11, G -11RPPPi,0i,1i,G-1PPPG-1,0G-1,1G-1,G-1其中：当 0= 0，90 时，R = 2 x W x (W -1)；当 0 = 45，135 时，xyR = 2x(W -1)x(W -1)。考虑到地震数据中不同地质目标的不同尺度问题，数据分析窗xyW xW 对地震纹理的影响也很大，也是一个值得讨论的参数。 xy图2演示了从图1左上角取出3x3大小的4阶灰度数据后计算灰度共生矩阵的结果， 可以发现矩阵 M 中的值在对角线集中的现象。图1地震道振幅数据灰度化示意图jjDIVSD 12n图2不同方向的灰度共牛矩阵(未归一化)示意图(a) 。；
6、(b) 90。； (c) 45。；初135。(3) 4 种灰度共牛矩阵的统计量特征纹理熵：Entropy三审 P log Pij ij11 + (i - j )2Pijij纹理均质性：Homogeneity三Hi纹理对比度：Contrast三审(i - j)2 Pij ij纹理能量:Energy三审P 2ijij3. 曲率属性计算曲率离不开求导，因此求导的方法决定了曲率的计算方法，这里主要介绍差分法计 算的倾角(Dip Angle)和最大正曲率(Most Positive Curvature)。为计算如式(1)所表达的曲面上某一点的曲率，用其周围网格点的值对局部二次曲面进行最小二乘法拟合。若采
7、用3X3网格单元作逼近(图3),由一阶和二阶导数的定义，得 式中各系数(式(2)式(6)。图3用于差分计算的网格单元z (x, y) = ax 2 + by 2 + cxy + dx + ey + f(1)1 d 2 z z + z + z + z + z + z z + z + za = 134679 2582 dx 212Ax 26Ax 242)71 d 2 z z + z + z + z + z + z z + z + zb = = 1 2 3 7 8 9 4 5 62 dy 212 Ay 26Ay 2(3)d2 zz + z z zc = 3719dxdy4AxAy了 dz z + z
8、 + z z z zd = 369147dx6Ax 2dz z + z + z z z ze =23789dy6Ay 2(5)式中：幷 z9是图5中所示层面网格结点1至9的值；Ax、Ay是网格结点间的距离。在此基础上，可计算各种曲率属性。(1)倾角dip = tan-i( d2 + e2)(2)最大正曲率Kpos =(a + b)+Ka b)2 + c 24. 频谱类层属性1) 有效带宽(Efective Bandwidth) 有效带宽是自相关函数的极值与双边自相关函数值的总和之比。计算公式为EB = r (0)/T 工 M(n)式中，r(o)是时窗内样点数据的零延迟自相关，工r(n)代表道两
9、旁所有自相n= M关的总和， T 是采样间隔。该属性反映地层的横向相似性，宽的带宽反映非均匀，复杂地层； 窄的带宽反映简单或圆滑的反射特征，可认为是均质的地层。带宽能帮助我们在数据体中识 别噪声区，有噪声的数据体比没有噪声的数据体有很明显宽的带宽。同时，它可用于指示沉 积环境的变化。提取时窗应大于2倍主波长，一般为50200ms.2) 波形长度(Arc Length)波形长度是指时窗内地震波振动曲线的长度。计算公式为：AL = - 迓、，(x (i +1) x (i)2 +T 2NT Vi=1式中，N代表时窗内样点个数，T代表采样间隔，x(i)代表时窗内第i个样点的振幅 值。该属性反映高振幅高
10、频率和高振幅低频率之间与低振幅低频率和低振幅高频率之间的差 别。比如页岩和砂岩的界面，一般有一些突变和高阻抗的反差，波形长度就用于页岩层序和 含砂量较高的层序之间的识别，带宽越小，它就越接近总绝对值振幅。它类似于反映反射层 的均质性。提取时应选取较大的时窗，一般为50100ms.3) 平均零交叉点数(Average Zero Crossing)平均零交叉点数也称为平均零交叉点频率，计算公式为N 1AZC =2 (f t )21式中， N 代表时窗内过零点的样点个数， t 、t 代表时窗内样点第一次和最后一次过 u12零点的反射时间。该属性相对于平均瞬时频率更稳定、更灵敏，是测量频率成分的另一种
11、方 法。时窗应大于2倍主波长，一般为50200ms.5. 振幅类层属性提取1) 均方根振幅(RMS Amplitude)：ARMS式中，N为时窗内样点个数，X：为时窗内第i个样点的振幅值。该属性可以指示孤立的 或极值振幅异常，追踪三角洲河道和含气砂岩等。一般情况时窗选取为50100ms.2) 总能量(Total Energy):A =丄工C x)Var N i式中，xi代表时窗内第i个样点的振幅值，N为时窗内样点个数。该属性适于描述层 序地层内岩性变化和含油气砂岩的变化。一般情况下，时窗选取为50100ms.3) 振幅的方差(Variance in Amplitude)：A 二丄 Y (x -
12、 X)Var Ni式中，N为时窗内样点个数，xi为时窗内第i个样点的振幅值，x是时窗内样点振幅值 的平均值。该属性是振幅相对摆动量的量度。提取时，时窗一般选取为50100ms.4) 平均能量(Average Energy)：E =丄2Avg N ii=1式中，N为时窗内样点个数，无为时窗内第i个样点的振幅值。该属性可用来分析有意 义的区段/层位的振幅异常，是一个检测亮点或暗点的关键属性。提取时，时窗一般选取为 50100ms.5) 最大绝对值振幅(Max Absolute Amplitude)：A = max (x.|)0 i N 1式中，X.为分析时窗内第i个样点的振幅值，N为分析时窗内样点
13、个数。该属性用于识 别岩性变化和含气砂岩，检测具体反射层内的振幅异常。提取时，一般选取较小的时窗： 20 100ms.6) 平均绝对值振幅(Average Absolute Amplitude)：A =迓 lx IAvg N ii=1式中，N为时窗内样点个数，x.为时窗内第i个样点的振幅值。该属性能识别层序地层 内整合的、丘状的或杂乱堆积的地层间差异。提取时，时窗一般选取为50100ms.7) 总振幅(Total Amplitude)：ATt =YI x.li=1式中，xi为时窗内第.个样点的振幅值，N为时窗内样点个数。该属性适于描述层序地 层内岩性变化和含油气砂岩的变化。提取时，一般选取较小
14、的时窗： 20100ms.8) 振幅的斜度(Skew in Amplitude):A二丄工(x - X)Skew N i式中，N为时窗内样点个数，xi为时窗内第i个样点的振幅值，x是时窗内样点振幅值的平均值。它是振幅相对摆动量的变化率。一般选取较小的时窗：20100ms.9) 振幅的峰态(Kurtosis in Amplitude)： A二丄工(x - X)Kurt Ni式中， N 为时窗内样点个数， xi 为时窗内第 i 个样点的振幅值， x 是时窗内样点振幅值的平均值。该属性是振幅相对摆动量的量度。用来确定记录相对平稳(波的到达相当少)的数据部分。提取时，一般选取较小的时窗： 20100ms.6. 地震纹理分析1) 能量(Energy)Energy = W p (z, j )2( 0-7)ij能量是图像灰度分布均匀程度和纹理粗细的一个度量。当图像较细致、均匀 时，能量值较大，最大时为1，表明区域内图像灰度分布完全均匀；反之，当图 像灰度分布很不均匀、表面呈现出粗糙特性时，此时能量值较小。2) 对比度(Contrast)f=W (i - j)2 p (i, j)
