基于地震分级敏感属性融合的深层薄互储层定量表征方法
【技术领域】
[0001] 本发明涉及地球物理勘探及综合研究领域，尤其涉及到一种基于地震分级敏感属 性融合方法预测深层薄互储层。
【背景技术】
[0002] 地震属性是表征地震波几何形态、运动学特征、动力学特征和统计特征的物理量， 在当前的石油勘探综合研究中是进行储层预测的一个重要方面。由于局限于使用单一属性 预测储层，预测结果往往多解性很强，这种多解性问题可以通过多属性融合方法加以改 善。在测井资料的基础上，研究地震属性特征，分析各属性对储层的敏感性，应用多属性 融合方法将几种属性融合于一体，再利用井位计算出各种地震属性的融合比重，预测有利 储层的发育位置。该思路有机地结合了各属性的优点，能有效降低储层预测结果的多解 性。但随着地层深度的加深，产状的变化以及储层厚度的变薄，多属性融合的预测精度会 大大降低。
【发明内容】
[0003] 本发明的目的是针对现行技术对深层薄互层储层预测精度较低的问题，提供一种 利用地震分级（区带分级和时窗分级）的基于地震分级敏感属性融合的深层薄互储层定量 表征方法，有效提高了针对深层薄互储层地震属性分析的效率，大大提高了预测精度。特别 适用于储层埋藏深、断块复杂、储层厚度薄、边界难以确定的地区。
[0004] 本发明的目的可通过如下技术措施来实现：
[0005] 基于地震分级敏感属性融合的深层薄互储层定量表征方法：大时窗段的地质背 景下，针对目的层段的地质目标分区域提取的合理小时窗段属性，对有利储层进行精细雕 亥IJ，最终与大背景融合，得到针对地质目标的属性预测图。该方法包括两个关键步骤：一是 大地质背景下，利用常规地震多属性相关降维、针对目的层段优选目标敏感属性。二是基于 地震分级对属性采取先分级：确定不同背景下目标区域合理的小时窗段；再融合：分序级 小时窗参数下敏感属性在融合定量预测研究区内储层范围。其中所述的分级和融合采用实 钻资料约束下的支持向量机（SVM)算法实现。
[0006] 发明的具体技术方案：
[0007] 1、结合井震高频层序划分的沉积背景，对常规地震多属性相关降维、对目的层段 进行属性优选。完成常规多属性融合图，预测深层薄互储层。
[0008] 2.在常规多属性预测储层的基础上，对全区储层进行分级处理。大范围大范围的 区域，储层的产状变化较大，大时窗的预测不能满足小范围砂体的精细雕刻，因此对全区进 行分级处理，包括：区带分级和时窗分级。区带分级主要是利用测井、钻井资料，结合沉积特 征，在已钻井的基础上，将需要详细刻画的储层确定出来，在平面上划分出区域；时窗分级 的原则是在砂体分布比较集中的地方可以用同一个小时窗，在不集中的地方可以进一步的 划分目标区域再确定时窗，直到能把砂体精细刻画。但小时窗的范围必须在大时窗的范围 内，不能超出，主要目的是为了更细致的描述砂体的边界。
[0009] 3、在确定分级的基础上，对各个分级区带和分级时窗内的储层进行分级预测。其 算法包含两个阶段：训练阶段和预测识别阶段。训练阶段主要是为了建立基于支持向量机 (SVM)的多属性情况识别模型，主要步骤为：对地震属性训练、集中属性信息进行归一化： 提取特征向量。建立地震属性特征向量训练集：由此为基础建立支持向量机（SVM)分类模 型；识别阶段主要步骤为：对要分级的地震属性信息进行归一化；提取其特征向量。其中分 级建立支持向量机（SVM)识别模型过程中还加入了实钻井的约束：提取属性特征向量时加 入钻井储层参数的约束，是预测结果与实钻井吻合度更高。
[0010] 4、对分区带分时窗的属性进行第二次融合，得到全区分级属性融合属性。运算中 对优选属性随机取若干采样点，进行分类函数f(x)的运算，计算分类阀值，按照最优化原 贝1J，即该阀值能将属性中的不同分类正确地分开，并且折中考虑使样本错分率最小和分类 间隔最大。不断修改分类阀值，直到不同分类中的样本点误差都收敛至可接受的范围内。这 样得到的融合结果避免了多属性融合的多解性，使预测结果更加精确。
[0011] 5.完成的基于地震分级敏感算法大背景下的属性分类结果：针对分类结果进一 步细分目标进行小时窗的属性预测，通过属性优化的手段与反应背景的原始属性信息再融 合，得到地震分级属性融合预测图。
[0012] 6.利用地震分级属性融合储层预测图件结合实钻资料的储层参数，半定量化的预 测出深层薄储层的展布范围，使预测结果与实钻井吻合度更高，更精确，结果可推广应用。
[0013] 上述方案步骤（3)识别阶段中对要分级的地震属性信息进行归一化，提取其特征 向量；将该向量作为支持向量机（SVM)输入进行分类函数的运算、每个支持向量对应一个 分类阀值b，对分类阀值b低于特定数值的的归为一类，依此确定分类界限；提取属性特征 向量时加入钻井储层参数的约束，通过井点附近地震属性信息与大背景下地震属性比对 分析，不断计算收敛，优化识别建模过程中特征向量的提取，使模型更接近于实际地质情 况，降低多解性。
[0014] 上述方案步骤（4)对分区带分时窗的属性进行第二次融合，得到全区分级属性融 合属性，对优选属性随机取若干采样点，进行分类函数f (X)的运算，计算分类阀值，按照最 优化原则，即该阀值能将属性中的不同分类正确地分开，并且折中考虑使样本错分率最小 和分类间隔最大，不断修改分类阀值，直到不同分类中的样本点误差都收敛至可接受的范 围内；在该过程中，以分级模型为督导进行训练，得到若干类基于分类函数的特征向量，每 类特征向量在地震属性上都对应一个分类集合，至此分类过程结束；整个过程以支持向量 机（SVM)为核心算法，算法的实现方式如下：
[0015] Η为分类线，满足线性方程：w*x+b = 0, HI、H2分别是过各类中离分类线最近的 样本且平行于分类线Η的直线，它们之间的距离称做分类间隔；所谓的最优分类线就是要 求分类线不但能将两类正确分开，即训练错误率为〇,而且使分类间隔2/ || w ||最大；针对 线性可分集合（Xl，yi)对分类线方程进行标准化处理后满足yi(w*xi+b)-l = 1，一η ye {-1，1}，使分类间隔2/11 W II最大等价与使II W II2最小，满足前式且使II W II2/2最小 的分类面称为最优分类面，H1、H2上的样本点称为支持向量，求取最优分类面也就是求取线 性方程中最佳（w, b)。即在yi (w*xi+b)_l彡0i = 1，…n y e {-1，1}的约束下求|| w || 2/2的最小值；定义如下的Lagrange函数：
[0016]
[0017] 其中ai彡0为拉格朗日系数；求II w II 2/2的极小值就是对w, b求拉氏函数的极小 值，对L函数做微分并令其等于0,拉氏函数及转换为：在Σ yA = 0, a;彡0, i = 1, ···]!约 束下求拉格朗日系数ai对函数：
[0018]
[0019] 其极大值最优解满足aJyJw^Xi+tO-l) = 0, i = 1,…η,解中a;不为零的样本 即支持向量，也就是H1、H2上的样本点，因此W可简化为：
[0020]
[0021] 求解后得到的最优化分类函数为：
[0022]
[0023] 其中sgn为分类函数、b为分类阀值，X为给定样本、xi, yi为线性可分样本集；对 于未知样