体曲率获取方法和设备的制作方法.docx

体曲率获取方法和设备的制作方法
体曲率获取方法和设备的制作方法
本发明提供了用于获取体曲率的方法和设备，其中，所述方法包括：获取三维地震数据的每个地震数据采样点的线方向、道方向和时间/深度方向的梯度矢量；基于获取的梯度矢量建立地震数据采预测方法则是一种既直接又现实的裂缝预测方法，并且在国内外部分油气田也有成功应用的先例。综合比较而言，基于裂缝纵波地震响应的预测方法更为实际可行。
[0008]在基于裂缝纵波地震响应的裂缝预测中，包括基于叠前地震资料的裂缝预测和基于叠后地震资料的裂缝预测方法两大类，其中后者又主要包括相干、曲率等方法。
[0009]目前，曲率计算都是以倾角和方位角作为输入，且主要采用复数道分析算法或离散倾角扫描算法。复数道分析算法由于只是用单点的倾角信息，没有平滑掉背景噪声，因此效果较差；而离散倾角扫描算法的效果虽好，但是计算速度太慢。
【发明内容】
[0010]本发明目的在于提供体曲率获取方法和设备。
[0011]为实现上述目的，根据本发明的一方面，提供了一种获取体曲率的方法，包括:A)获取三维地震数据的每个地震数据采样点的线方向、道方向和时间/深度方向的梯度矢量；B)基于获取的梯度矢量建立地震数据采样点的梯度结构张量；C)对梯度结构张量的元素进行平滑；D)基于平滑后的梯度结构张量获取P波属性和Q波属性，其中，P波属性是一种反应地下岩层沿线方向的变化率的属性，Q波属性是一种反应地下岩层沿道方向的变化率的属性；E)基于P波属性和Q波属性获取体曲率。
[0012]可选地，步骤C)中可通过下面的高斯滤波器对梯度结构张量的元素进行平滑:
【权利要求】
，包括: A)获取三维地震数据的每个地震数据采样点的线方向、道方向和时间/深度方向的梯度矢量； B)基于获取的梯度矢量建立地震数据采样点的梯度结构张量； C)对梯度结构张量的元素进行平滑； D)基于平滑后的梯度结构张量获取P波属性和Q波属性，其中，P波属性是一种反应地下岩层沿线方向的变化率的属性，Q波属性是一种反应地下岩层沿道方向的变化率的属性； E)基于P波属性和Q波属性获取体曲率。
，其特征在于，步骤C)中通过下面的高斯滤波器对梯度结构张量的元素进行平滑:
G(Xi’ Vj,z;;σ) = cxp[(.v； + \? +?)/(2σ~)] 其中，Xp Yj和Zj表示沿着X，y和ζ轴的待平滑的梯度结构张量的元素为中心
的第j道，σ是平滑时窗。
，其特征在于，基于平滑后的梯度结构张量获取P波属性和Q波属性是对平滑后的梯度结构张量进行特征值分解得到P波属性和Q波属性， 其中，P波属性是进行特征值分解后的最大特征值所对应的特征向量的第一个元素和第二个元素、第三个元素的差的比值；Q波属性是进行特征值分解后的最大特征值所对应的特征向量的第二个元素和第二个元素、第三个元素的差的比值。
，其特征在于，步骤E)包括: 基于P波属性和Q波属性获取下面的反射面的二次曲面表达式的系数: ζ (X，y) =ax2+by2+cxy+dx+ey+f 其中，a的值是P波属性沿X方向的偏导数值的一半，b的值是Q波属性沿y方向的偏导数值的一半，c是P波属性沿X方向的偏导数