发明内容
本发明的示例性实施例的目的在于克服现有技术中的上述的和/或其他的问题。因此,本发明的示例性实施例提供了一种基于D-S证据理论的地震属性融合方法,通过将多种地震属性融合,以达到在快速实现属性融合的同时,有效提高地质岩性信息的准确性、进而有效地提高对储层空间展布形态、分布范围和储层岩石物理特征参数的描述与预测。
根据本发明的示例性实施例,提供了一种基于D-S证据理论的地震属性融合方法,所述方法包括下述步骤:步骤A、导入待融合的地震属性;步骤B、对导入的地震属性的各幅值分别进行归一化处理;步骤C、确定用于识别数据的识别框架;步骤D、利用确定的识别框架来描述每个属性数据;步骤E、将每个属性数据进行分块,以对每个属性数据的各个块构造证据体;步骤F、分别计算单证据体作用下识别框架中各命题的信度区间;步骤G、利用D-S合成规则计算所有证据体联合作用下的基本概率赋值和信度区间;步骤H、将D-S融合后的概率值作为加权系数,与属性相结合,从而实现地震属性的融合。
在步骤B,可以通过下式进行归一化处理:
在式(1)中,Xi为第i种属性的幅值,为第i种属性归一化处理后的幅值。
待融合的地震属性可以为均方根振幅地震属性和相干体切片地震属性。在步骤C,可以确定识别框架,以利用熵、平均梯度、点锐度来识别属性数据。
在步骤D,可以通过下述算法来描述每个属性数据:
基于熵值的评价算法
H=-∑pi lnpi
(2)
在式(2)中,pi为灰度值等于i的像素数与图像总像素数之比,ln为自然对数。
基于平均梯度的评价算法:
在式(3)中,M和N为图像的长度和宽度,g(i,j)为图像在(i,j)点的像素值。
基于点锐度的评价算法:
在式(4)中,m、n为图像的长和宽,df为灰度变化幅值,dx为像元间的距离增量。
在步骤E,可以将每个属性数据等分成n×m块,可以将对每个属性数据的各个块用熵、平均梯度、点锐度来评价所得到三个评价因子作为一个证据体。
在步骤F,可以将每个属性相同位置的块的证据体的值分别排序,可以以最小的值作为标准,根据证据体评价因子的值分别赋予概率值。
在步骤F,可以将属性数据评判为好,较好,一般这三档,并可以对应地赋予概率值0.6、0.4、0.2。
在步骤F,可以对每种评价算法所得的基本概率值进行归一化处理:
在式(5)中,i表示第i种属性,n指示属性数据分块所得的第n块,pni为第n块的第i个属性所得概率值,l为总的属性个数。
在步骤G,可以通过三个评价算法联合作用,将每个算法评价的概率融合成一个概率:
其中,
在式(6)中,K表示了证据冲突的额度,m1、m2、m3为三种评价算法得到的基本概率值,X,Y,W为三种评价算法,Z为属性数据相同的部分。
在步骤H,可以通过下式得到融合后的系数:
在式(7)中,Z
n为第n块融合后的系数,P
ni为融合后的系数,
为第i个属性的第n块数据。
本发明的示例性实施例还提供了一种基于D-S证据理论的地震多属性融合方法,通过对源属性数据进行特征提取,构造证据以获得基本概率值,将基本概率值在D-S证据理论指导下进行融合得到融合后的概率值并将其作为加权系数,进而将属性数据各部分小块与加权数相结合从而实现多种地震属性的融合。因此,可以解决利用单一地震属性信息解释复杂地质问题时单一地震属性信息所带来的多解性问题。
本发明的示例性实施例利用D-S证据理论将多种地震属性数据进行合成处理,具有简捷、快速可靠,效率准确度高,处理能力强等特点,对数据量较大的地震属性的融合效果尤为明显,有效提高地质岩性信息的准确性、进而有效地提高对储层空间展布形态、分布范围和储层岩石物理特征参数的描述与预测。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法与直接加权融合有同样的执行效率。与传统的直接加权融合方法相比,根据本发明的示例性实施例的算法本身就可以确定加权因子,不需要先验信息和人为干预,即使属性差异很大,融合的效果也比单一的加权融合好,能够把地质岩性信息更准确地表现出来,进而有效提高对储层空间展布形态、分布范围和储层岩石物理特征参数的描述与预测。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法效率准确度高,处理能力强,对数据量较大的地震属性的融合效果尤为明显,即使属性差异很大,融合的效果也比单一的加权融合好,能综合多种地震属性信息,将有效的地质信息、储层信息及含流体特征信息更准确地表现出来。根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法适用范围广,可以通过调整评价算法来适应不同的情况,适用于各种地震属性数据的融合。由于该方法不需要先验信息和人为干预,使用起来方便快捷,且综合了更多有效的地震属性有效信息,因此将在地质岩性分析、储层空间的描述与预测、含油气检测等领域有良好的应用前景。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法拓展了常规的属性数据融合方法,使其在地震属性分析,地震储层预测、岩性与流体识别取得更一步的良好效果。采用根据本发明的示例性实施例的方法可以融合多种地震属性,进行地震属性分析,地震储层预测、岩性与流体识别等研究,弥补使用单一地震属性信息参数解决地质现象存在的多解性缺陷。
具体实施方式
下文中,将参照附图来详细描述本发明的示例性实施例。然而,本发明的示例性实施例可以以许多不同的形式来实施,且不应该被局限于这里阐述的示例。相反,提供这些示例使得本公开将是彻底并完整的,并将向本领域技术人员充分地传达本发明的范围。
图1为根据本发明示例性实施例的基于D-S证据理论的地震属性融合方法的应用原理示意图,图2为根据本发明示例性实施例的基于D-S证据理论的地震属性融合方法的流程示意图。
参照图1和图2,根据本发明的示例性实施例,首先,可以导入待融合的地震属性,并可以在对导入的待融合的地震属性的幅值分别规范到[0,1]区间(归一化处理)后,将其作为常量输入。接下来,可以对输入的属性数据进行特征提取,构造证据,以获得基本概率值。在一个本发明的示例性实施例中,可以将输入的一个属性数据等分成n×m块小数据块,可以对每个小数据块进行基于熵、梯度和平均梯度的评价,并可以将得到的评价因子与其他属性对应的评价因子作比较。然后,可以以比较后的结果作为依据来求取基本概率值。可以将基本概率值通过D-S证据理论来融合得到融合后的概率值,可以将得到的融合后的概率值作为加权系数。最后,可以将属性数据的各个分得的小数据块块与加权数相结合,从而实现多种地震属性的融合。
具体地讲,根据本发明的一个实施例的基于D-S证据理论的地震属性融合方法可以包括下述步骤。
步骤A:导入待融合的地震属性;
步骤B:属性幅值的归一化处理,即,对导入的地震属性的各幅值分别进行归一化处理;
步骤C:在深入分析决策问题的基础上,构造系统识别框架;
步骤D:针对属性数据,构建基于识别框架的证据体(例如,熵值、点锐度、平均梯度值);
步骤E:将每个属性数据进行分块,搜集每个属性中各个块的证据体资料;
步骤F:综合每个属性的对应块的证据体,分别计算单证据体作用下识别框架中各命题的信度区间;
步骤G:利用D-S合成规则计算所有证据体联合作用下的基本概率赋值和信度区间;
步骤H:将D-S融合后的概率值作为加权系数,与属性相结合实现多地震属性融合。
下面,将参照图3a-图3c,以均方根振幅地震属性和相干体切片地震属性的融合为例,具体地说明根据本发明的一些示例性实施例的基于D-S证据理论的地震属性融合方法。图3为根据本发明示例性实施例的融合属性,其中,图3a为均方根振幅地震属性,图3b为相干体切片地震属性,图3c为融合的属性。
步骤A
在步骤A,可以导入均方根振幅属性和相干体切片属性。
步骤B
在步骤B,可以进行属性幅值的归一化处理。例如,可以对在步骤A导入的地震属性的幅值分别进行归一化处理。可以通过下式进行归一化处理:
在式(1)中,X
i为第i种属性的幅值,
为第i种属性归一化处理后的幅值。
步骤C
在步骤C,可以确定识别框架,利用熵、平均梯度、点锐度来识别属性数据。
步骤D
在步骤D,可以构造证据体,通过熵、平均梯度、点锐度来描述每个属性数据。示例性的算法如下:
基于熵值的评价算法:
H=-∑pi lnpi
(2)
在式(2)中,pi为灰度值等于i的像素数与图像总像素数之比。
基于平均梯度的评价算法:
在式(3)中,M和N为图像的长度和宽度,g(i,j)为图像在(i,j)点的像素值。
基于点锐度的评价算法:
在式(4)中,m、n为图像的长和宽,df为灰度变化幅值,dx为像元间的距离增量。式(4)可描述为:逐个对图像中的每点取8邻域点与之相减。先求8个差值的加权和加权的大小取决于距离,距离远,则权小,如45°和135°方向的差值需乘以
,再将所有点所得值相加除以像素总个数。其中df/dx较边缘锐度算法之所以由平方变为绝对值的原因主要是,整幅图像的8邻域灰度差值和已足够反映图像的灰度分布情况。若加上平方的放大作用,反而会因噪声等众多因素而造成该参数的波动。式(4)可以理解为对图像每点周围的灰度扩散程度的统计,即扩散程度越剧烈,其值越大,图像越清晰。从另一个角度看,该算法可以近似地等同于对图像点扩展函数能量分布情况的统计。这里,属性数据与图像数据一样,是一个二维数据,将属性数据归一化后,可以将属性数据看做图像数据。
步骤E
在步骤E,可以取得证据体。例如,可以将每个属性数据等分成n×m块,可以对每个属性数据的各个块构造证据体,即,可以将对每个块用熵、平均梯度、点锐度来评价所得到三个评价因子作为一个证据体;
步骤F
在步骤F,可以进行基本概率值的赋值。例如,可以将每个属性相同位置的块的证据体的值分别排序,以最小的值作为标准,根据证据体评价因子的值赋予三种概率值,将属性数据评判为好、较好、一般这三档,对应的概率值为0.6,0.4,0.2。为了便于融合,最后对每种评价算法所得的基本概率值进行归一化处理。公式如下:
在式(5)中,i指示第i种属性,n指示属性数据分块所得的第n块,pni为第n块的第i个属性所得概率值,l为总的属性个数。
步骤G
在步骤G,可以实现属性融合。例如,可以通过三个评价算法联合作用,将每个算法评价的概率最终融合成一个概率,如下式所示:
其中,
在式(6)中,K表示了证据冲突(Evidential Conflict)的额度,m1、m2、m3为三种评价算法得到的基本概率值,X,Y,W为三种评价算法,Z为属性数据相同的部分。
步骤H
在步骤H,可以实现属性对步骤G所得每个属性的各个块的融合系数,则属性融合的对应幅值输出为各属性的各个幅值与其对应的D-S证据理论的融合系数的乘积之和。例如,假设Zn为融合后的属性数据,则有:
在式(7)中,Zn为第n块融合后的属性数据,Pni为融合后的系数,为第i个属性的第n块数据。
图3a为均方根属性数据的显示图,图3b为相干体属性数据的显示图。可以看到图3a和图3b中由于属性单一,存在很多冗余信息。然而,如通过上述方法融合得到的图3c所示,通过重叠部分可以看出,大大降低了多解性,突出了感兴趣区域。
本发明的示例性实施例利用D-S证据理论将多种地震属性数据进行合成处理,具有简捷、快速可靠,效率准确度高,处理能力强等特点,对数据量较大的地震属性的融合效果尤为明显,有效提高地质岩性信息的准确性、进而有效地提高对储层空间展布形态、分布范围和储层岩石物理特征参数的描述与预测。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法与直接加权融合有同样的执行效率。与传统的直接加权融合方法相比,根据本发明的示例性实施例的算法本身就可以确定加权因子,不需要先验信息和人为干预,即使属性差异很大,融合的效果也比单一的加权融合好,能够把地质岩性信息更准确地表现出来,进而有效提高对储层空间展布形态、分布范围和储层岩石物理特征参数的描述与预测。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法效率准确度高,处理能力强,对数据量较大的地震属性的融合效果尤为明显,即使属性差异很大,融合的效果也比单一的加权融合好,能综合多种地震属性信息,将有效的地质信息、储层信息及含流体特征信息更准确地表现出来。根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法适用范围广,可以通过调整评价算法来适应不同的情况,适用于各种地震属性数据的融合。由于该方法不需要先验信息和人为干预,使用起来方便快捷,且综合了更多有效的地震属性有效信息,因此将在地质岩性分析、储层空间的描述与预测、含油气检测等领域有良好的应用前景。
根据本发明的示例性实施例的基于D-S证据理论的地震属性融合方法拓展了常规的属性数据融合方法,使其在地震属性分析,地震储层预测、岩性与流体识别取得更一步的良好效果。采用根据本发明的示例性实施例的方法可以融合多种地震属性,进行地震属性分析,地震储层预测、岩性与流体识别等研究,弥补使用单一地震属性信息参数解决地质现象存在的多解性缺陷。
虽然已经示出并描述了本发明的示例性实施例的示例,但是本领域技术人员应该理解的是,本发明的示例性实施例不限于此,在不脱离根据权利要求所限定的本发明的精神和范围的情况下,可以对这些示例性实施例进行各种修改。