一种基于直方图匹配算法的属性连片处理方法
【技术领域】
[0001] 本发明属于油气勘探开发中地震属性处理领域，具体涉及一种基于直方图匹配算 法的属性连片处理方法。
【背景技术】
[0002] 在地震资料采集过程中，由于受地质勘探目标、经济能力、勘探技术、勘探周期等 因素的影响与制约，相邻区块间获得的地震资料往往分属于不同时间采集、处理的成果，后 期为了全面掌握区域构造、岩性分布特征，需要对不同区块的地震资料进行连片解释，但由 于不同时期采集、处理资料的差异性，在连片解释前需要进行不同区块地震资料的叠前连 片处理工作，以解决地震资料品质不一致等问题。
[0003] 针对不同区块地震资料的连片处理，往往需要处理人员在工区网格重新划分的基 础上，开展叠前地震资料的地表一致性振幅补偿、预测反褶积、全区统一速度分析、地表一 致性剩余静校正、面元均化处理等。当地震资料差异较大时，该处理流程工作量极大，虽然 能够得到较好的地震资料品质用于后续的属性提取及解释，但不能很好满足油田开发所要 求的时效性。
【发明内容】
[0004] 本发明的目的在于解决上述现有技术中存在的难题，提供一种基于直方图匹配算 法的属性连片处理方法，解决不同区块由于不同时期采集、处理地震资料不同所造成计算 的属性数值差异问题。通过不同区块属性数值的去差异化处理，实现不同区块间地震属性 数据的较好一致性，最终完成不同区块属性的连片拼接，为地震资料的连片解释提供地质 依据，满足地震资料连片解释快速处理的要求。
[0005] 本发明是通过以下技术方案实现的：
[0006] -种基于直方图匹配算法的属性连片处理方法，包括：
[0007] (1)参考属性数据的均衡化处理：针对参考区域范围内的原始地震资料的属性数 据yy分布特征，即目标直方图，由此计算变换函数R(yu);
[0008] (2)原始属性数据均衡化处理：针对需要进行均衡化变换区域范围内的原始地震 资料的属性数据Xu分布特征，即原始直方图，由此计算变换函数T( Xl> ]);
[0009] (3)对参考属性数据的变换函数R(yi,,)求逆，得到R 1G1,,);
[0010] (4)原始属性数据直方图匹配处理：利用参考变换函数的逆函数R1G1,,)进行原 始属性数据均衡化处理；
[0011] (5)基于直方图匹配处理后的属性数据进行空间直接拼接，得到连片属性数据结 果。
[0012] 所述步骤（1)是这样实现的：
[0013] 利用公式
1计算变换函数，其中，i、j分别为该属性的行号和 列号，k为该数据yl,,的最大值，i是数据yl,,分布的直方图函数，该函数通过统计数据 yl, 值的分布区间得到。
[0014] 所述步骤（2)是这样实现的：
[0015] 利用公式
计算变换函数，其中，i、j分别为该属性的行号和 列号，k为该数据X1,,的最大值，ha是数据X1,,分布的直方图函数，该函数通过统计数据X 1, 值的分布区间得到。
[0016] 所述步骤（4)是这样实现的：
[0017] 将原始属性数据变换函数T (X1,,)代入到R 1G1, ,)中，实现原始属性数据均衡化处 理。
[0018] 所述步骤（5)是这样实现的：
[0019] 根据直方图匹配处理后的不同的区块属性数据，按坐标直接把属性数据显示在同 一张图像上。
[0020] 与现有技术相比，本发明的有益效果是：塔河油田不同区块地震资料采集、处理结 果的不一致性，造成相同地震属性值域分布的差异，原始的不同区块地震属性集中显示的 平面结果，其在属性能量上及地质现象刻画能力上存在较大的差别，基于直方图匹配处理 后的五连片地震属性，其多区块属性的一致性得到进一步的提升，在拼接位置无明显属性 差异，属性过渡自然平滑，譬如一间房组顶面（T74)相干属性比原始分区块相干属性更加 能够反映区域上不同尺度断裂分布特征。证明该方法对属性一致性处理有效，在塔河油田 连片性地震属性处理中取得了较明显的效果。
【附图说明】
[0021] 图1本发明方法的步骤框图
[0022] 图2-1目标直方图
[0023] 图2-2原始直方图
[0024] 图2-3基于图2-1的目标数据直方图对图2-2的原始数据直方图进行直方图匹配 后的直方图
[0025] 图3-1强振幅聚类属性连片处理前的图片
[0026] 图3-2强振幅聚类属性连片处理后的图片
[0027] 图4-1相干断裂检测属性连片（局部）处理前的图片
[0028] 图4-2相干断裂检测属性连片（局部）处理后的图片。
【具体实施方式】
[0029] 下面结合附图对本发明作进一步详细描述：
[0030] 本发明针对塔河油田主体区不同区块间不同时期采集及所处理地震资料的不一 致性问题，造成分区块提取的地震平面属性难以开展连片解释。通过开发基于直方图匹配 算法对不同区块地震资料平面属性的再处理，解决由于属性能量分布差异所造成的平面属 性不一致问题，获得的连片属性具有较好的一致性，提高基于地震属性连片识别断裂及古 水系分布特征的能力。
[0031] 直方图匹配算法来自于图像处理领域，通过将图像直方图以标准图像的直方图为 标准作变换，使两图像的直方图相同或近似，从而使两幅图像具有类似的色调和反差，达到 原始图像所要求的效果。
[0032] (1)技术原理
[0033] 本发明针对塔河油田主体区五块高精度三维地震资料的连片解释问题，借助于图 像处理技术中的直方图匹配算法，通过将原始数据分布的直方图以目标数据分布的直方图 为标准作变换，最终使得两数据分布的直方图相同或近似，从而使两类数据具有类似的显 示形式，实现五连片属性匀化处理。本发明方法可以用于任意数量的连片处理解释，只要利 用一块三维资料作为标准，其它区块资料都按直方图匹配算法向标准资料靠拢即可实现匀 化处理。
[0034] 实际应用中以五块三维数据中的一块地震属性数据为标准，建立属性分布标准直 方图，通过数据分布特征映射其余区块属性数据分布特征，使得其余属性数据的分布特征 与标准直方图分布特征相一致，从而实现连片地震属性数据图形显示及地质现象识别能力 上的一致性。
[0035] 其基本原理为：
[0036] 其中yi> ]表示参考区域范围内地震属性数据，其属性值的直方图均衡化变换函数 表示为
得到的结果为参考属性值变换后属性均质Zb。
[0037] 相应地，针对需要进行均衡化变换区域范围内的地震属性数据X1,,，其属性值的直 方图均衡化变换可表示为
[0038] 若要使属性数据X1,,直方图均衡化变换后的属性均质等于Zb，则利用参考属性数 据直方图均衡化变换函数R(Yu)的逆函数来实现，即利用RlTU^)]。这样，由上式所变 换后的属性数据 Xlu和参考属性数据yu将具有相似的数据