[发明专利]基于线形自适应灰阶化理论的地震属性突变界线求取方法

权利要求书

1.基于线形自适应灰阶化理论的地震属性突变界线求取方法，其方法步骤为：

步骤一：用地震数据计算得到一个二维层属性数据；

步骤二：将层属性上所有点按值从大到小排列到一个表中；

步骤三：设置一个空的二维数组，按照表中点原来所在的二维数组的坐标，对这个新的二维数组对应位置的点周围的点进行检查；

步骤四：检查过程中，判断该点周围情况：如果周围没有被检查过的点，则赋予该点一个新值k，可以视作一个区域编号，同时k＝k+1,作为下一个新区域的编号；如果周围有一个检查过的点，或者有多个但所有临近点区域编号都相同，那么该点赋予邻近点相同的区域编号；如果周围有多个检查过的点但区域编号不同，那么该点赋予临近点中最小的那个区域编号，同时，在整个二维数组中，搜索与其它区域编号的相同的点，均赋予此次的最小区域编号，并且记录每个区域编号被更改过的点的数目；

步骤五：如果单个区域编号更改过的点数目大于一个数n，那么储存该点对应的属性值，作为一个属性值域的分割点；

步骤六：回到步骤三，直到所有点都被检查过，最后会得到一个属性值分割点数组，用该数组对地震属性值域进行分割，并将分割好的值区间从小到大赋予一个整数编号，然后将层属性每个点的值转化为其在值域中对应值区间的编号，地震二维层属性数据即转化为了灰阶数据；

步骤七：以N_x×N_y为窗口扫描前一步得到的灰阶数据，计算灰度共生矩阵，从位置(x,y)灰度为i的点出发，统计与其距离为δ＝(Δx²+Δy²)^1/2、灰度为j的点同时出现的概率P

对于地震属性，考虑各向异性，θ取值为0°、45°、90°和135°，距离δ取1；对于0°和90°灰度共生矩阵，∑_i∑_jP(i，j，θ)＝2×N_x×(N_y-1)，对于45°和135°灰度共生矩阵，∑_i∑_jP(i，j，θ)＝2×(N_x-1)×(N_y-1)；

再对四个矩阵内元素分别进行归一化，使每个矩阵内元素的和为1

其中∑_i∑_jP(i，j，θ)可以由在灰阶数据上计算灰度共生矩阵的窗口大小N_x、N_y确定；

步骤八：在步骤七得到的灰度共生矩阵基础上，计算该点变化值，计算公式为

C(x，y)即为点(x,y)变化值结果；

步骤九：当灰阶数据所有点都完成步骤七、八的扫描计算后，得到了变化值组成的二维数据，其极大值点所连成的线即为该层属性中的突变界线。

2.根据权利要求1所述的基于线形自适应灰阶化理论的地震属性突变界线求取方法，其特征在于：步骤四中记录的每个区域编号被更改过的点的数目，在步骤五首先通过该数目是否超过n，判断在检测点附近存在的点是否包含于一个有效区域中，即该区域点数是否达到最小区域限制条件，进而确定该点是否为分割了不同区域的突变界线中的点，这样可以有效避开随机噪音形成的伪突变界线。

3.根据权利要求1所述的基于线形自适应灰阶化理论的地震属性突变界线求取方法，其特征在于：步骤五中的n为最小区域所具有的点数，随计算过程动态变化；如果检查过程中，间隔很短连续触发步骤四，则n增加；一旦储存值的条件触发，n急速减少，但不小于一个最小值。

4.根据权利要求1所述的基于线形自适应灰阶化理论的地震属性突变界线求取方法，其特征在于：其中在步骤七得到的灰度共生矩阵基础上，计算该点变化值，还可以采用如下公式：

5.根据权利要求4所述的基于线形自适应灰阶化理论的地震属性突变界线求取方法，其特征在于：其中灰度共生矩阵在得到该点四个不同θ的灰度共生矩阵的变化值之后，取四个中最大的一个作为该点最终的变化值结果：

C(x，y)＝MAX{C(x，y，θ)}。