>[0136] 标准化处理单元1103,用于对待预测区域的叠前地震数据体进行标准化处理；
[0137] 动力学参数确定单元1104,用于标准化后的叠前地震数据体作为初始数据，利用 所述海色矩阵进行反演迭代，获得待预测区域的动力学参数；其中，所述动力学参数包括： 纵波速度Vp、横波速度Vs和密度P ;
[0138] 转换单元1105,用于对所述动力学参数进行转换，得到待预测区域的弹性参数；
[0139] 流体敏感因子确定单元1106,用于利用所述待预测区域的弹性参数确定待预测区 域的流体敏感因子的取值情况；
[0140] 预测单元1107,用于根据所述待预测区域的流体敏感因子的取值情况预测含油气 储层分布。
[0141] 实施例：
[0142] 区域叠前三维地震处理面积312km2。目前区内完钻14 口井，在平湖组和宝石组钻 遇油气层，并获得工业油流，证明研究区具有良好的勘探前景。叠前研究区（孔雀亭井区） 位于平北地区的北部，区内完钻7 口井，其中有4 口井具有横波资料。针对西湖凹陷平湖组 目前勘探现状，本案例主要从叠前地震资料入手，将丰富的叠前速度密度信息和弹性参数 与先进的流体识别技术结合起来，识别含油储层和含气储层。
[0143] 由于研究区平湖组埋深较深，为了研究目的层的物性情况，我们选取KQT-2井平 湖组进行统计分析，如图2所示，为本实施例的KQT-2井平湖组测井曲线统计分析图。从 KQT-2井可以看出在平湖组4114-4128m、4150-4169m处的气层对应泊松比为0. 05-0. 23,储 层孔隙度为11-15 %，渗透率为3-8md，含气饱和度为30 % -70 %。可知，在研究区的深层目 的层处，若泊松比、纵横波速度比、拉梅系数乘密等参数能较好的指示储层和流体，则在相 应深度处的储层物性条件也较好。
[0144] 如图3所示，为本实施例的叠前地震资料分析图之一；如图4所示，为本实施例的 叠前地震资料分析图之二；如图5所示，为本实施例的叠前地震资料分析图之三；如图6所 示，为本实施例的叠前地震资料分析图之四。受覆盖次数影响，浅层地震资料信噪比较低， 同相轴连续性变差。同一条剖面振幅能量和频率均有差异，AVA特征明显。对过KQT-I井 的角道集剖面分析发现：同一条地震测线不同角度剖面的振幅、频率均有差异，随角度增大 振幅能量减弱，AVA变化明显反映了岩性和流体的差异；有利于叠前多参数反演和流体的 检测。
[0145] 如图7所示，为本实施例的不同角道集放大分析图。经临界角分析可知利用近似 的公式只有在角度小于30°的情况下，反演结果与实际结果才会吻合得较好，因此此次反 演选择了前三个角道集作为反演的输入。
[0146] 如图8所示，为本实施例的实际地震资料的含油性特征图，如图9所示，为本实施 例的实际地震资料的AVO特征图。图8与图9相比较可知，随角度增加地震波振幅、频率逐 渐衰减，尤其在15°到20°衰减最明显。同一剖面、相同时窗，振幅能量和频率均有差异， 具有明显的AVA特征。
[0147] 如图10所示，为本实施例的过KQT-2井的M值剖面图。通过叠前近似公式反演得 到的动力学参数数据体，将动力学参数数据体转换为弹性参数数据体，根据计算得到的弹 性参数数据体，排除异常值后，计算M值数据体，分析M值范围，根据测井解释成果，分析不 同类型的储层M值范围，再根据这个阈值，分析和预测本区储层的含油气性。结合图中过 KQT-2井的M值剖面和测井解释成果可以看出，水层的M值要比气层的M值偏高。水层的M 值在0. 23-0. 30之间而气层的M值在0. 08-0. 22之间，根据这一阈值，可以判断全区储层的 含油气性。
[0148] 下表2为叠前储层预测吻合率数据，通过井吻合率对比发现，通过计算M值得到的 剖面，与井的吻合率较高，说明该方法真实可信。
[0151] 以上所述的【具体实施方式】，对本发明的目的、技术方案和有益效果进行了进一步 详细说明，所应理解的是，以上所述仅为本发明的【具体实施方式】而已，并不用于限定本发明 的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含 在本发明的保护范围之内。
【主权项】
1. 一种基于多参数梯度向量和海色矩阵的油气检测方法，其特征在于，包括： 利用Aki-Richards近似公式建立叠前反演目标函数； 对所述叠前反演目标函数进行泰勒展开，得到海色矩阵；其中，所述海色矩阵中的元素 值通过对所述Aki-Richards近似公式求导得到； 对待预测区域的叠前地震数据体进行标准化处理； 标准化后的叠前地震数据体作为初始数据，利用所述海色矩阵进行反演迭代，获得待 预测区域的动力学参数；其中，所述动力学参数包括：纵波速度Vp、横波速度Vs和密度P; 利用动力学参数与弹性学参数之间的关系式对所述动力学参数进行转换，得到待预测 区域的弹性参数； 利用所述待预测区域的弹性参数确定待预测区域的流体敏感因子的取值情况； 根据所述待预测区域的流体敏感因子的取值情况预测含油气储层分布。2. 如权利要求1所述的方法，其特征在于，所述流体敏感因子的表达式为：其中，Μ表示流体敏感因子；λ、μ均表示弹性参数。3.如权利要求1或2所述的方法，其特征在于，所述弹性参数的表达式为： λ=p(Vp2-2Vs2) μ= 0· 5*(PVp2_λ) =pVs2 其中，λ、μ均表示弹性参数，Vp表示纵波速度，￥；3表示横波速度，p表示密度。4.一种基于多参数梯度向量和海色矩阵的油气检测装置，其特征在于，包括： 叠前反演目标函数建立单元，用于利用Aki-Richards近似公式建立叠前反演目标函 数； 海色矩阵确定单元，用于对所述叠前反演目标函数进行泰勒展开，得到海色矩阵；其 中，所述海色矩阵中的元素值通过对所述Aki-Richards近似公式求导得到； 标准化处理单元，用于对待预测区域的叠前地震数据体进行标准化处理； 动力学参数确定单元，用于标准化后的叠前地震数据体作为初始数据，利用所述海色 矩阵进行反演迭代，获得待预测区域的动力学参数；其中，所述动力学参数包括：纵波速度 Vp、横波速度Vs和密度P; 转换单元，用于对所述动力学参数进行转换，得到待预测区域的弹性参数； 流体敏感因子确定单元，用于利用所述待预测区域的弹性参数确定待预测区域的流体 敏感因子的取值情况； 预测单元，用于根据所述待预测区域的流体敏感因子的取值情况预测含油气储层分 布。5.如权利要求4所述的装置，其特征在于，所述流体敏感因子确定单元确定的流体敏 感因子的表达式为：其中，Μ表示流体敏感因子；λ、μ均表示弹性参数。6.如权利要求4或5所述的装置，其特征在于，所述转换单元得到的弹性参数的表达式 为：其中，λ、μ均表示弹性参数，Vp表示纵波速度，￥；3表示横波速度，p表示密度。
【专利摘要】本发明涉及一种基于多参数梯度向量和海色矩阵的油气检测方法及装置，其中，方法包括：利用Aki-Richards近似公式建立叠前反演目标函数；对所述叠前反演目标函数进行泰勒展开，得到海色矩阵；其中，所述海色矩阵中的元素值通过对所述Aki-Richards近似公式求导得到；对待预测区域的叠前地震数据体进行标准化处理；标准化后的叠前地震数据体作为初始数据，利用所述海色矩阵进行反演迭代，获得待预测区域的动力学参数；对所述动力学参数进行转换，得到待预测区域的弹性参数；利用所述待预测区域的弹性参数确定待预测区域的流体敏感因子的取值情况；根据所述待预测区域的流体敏感因子的取值情况预测含油气储层分布。
【IPC分类】G01V1/30, G01V1/28
【公开号】CN105301644
【申请号】CN201510794474
【发明人】黄捍东, 董月霞, 刘洪昌, 向坤
【申请人】中国石油大学(北京), 中国石油天然气股份有限公司冀东油田分公司
【公开日】2016年2月3日
【申请日】2015年11月18日