[发明专利]使用各向异性参数来生成时间偏移图像道集的方法和系统

权利要求书

1.一种用于储层表征和解释的方法，所述方法包括：

在勘探区域上设置包括多条炮点线和多条接收器线的网格；

沿着所述多条炮点线布置多个炮点，其中每条炮点线上均布置有多个炮点；

获取由所述多个炮点生成的叠前地震数据，其中所述叠前地震数据包括具有水平偏移轴和垂直时间轴的共炮点道集；

使用基于长偏移距P波数据的常规非椭圆率各向异性参数分析来预处理初始非椭圆率各向异性参数；

使用基于近偏移距P波数据的常规双曲线时差分析或者基于长偏移距P波数据的预偏移非双曲线时差分析来对基线正常时差速度进行预处理；

将所获取的叠前地震数据、所获取的初始非椭圆率各向异性参数和所获取的基线正常时差速度存储到存储资源中；

使用所存储的叠前地震数据、所存储的初始非椭圆率各向异性参数和所存储的基线正常时差来计算各向异性基尔霍夫时间偏移局部图像；

从所计算出的各向异性基尔霍夫时间偏移局部图像中生成具有水平偏移轴和垂直时间轴的远偏移距图像道集；

将所生成的具有水平偏移轴和垂直时间轴的远偏移距图像道集存储到存储资源中；

使用所存储的远偏移距图像道集来计算多个相似度参数；

从所计算出的多个相似度参数中生成最终相似度体积；

将最终相似度体积存储到存储资源中；

从所存储的最终相似度体积的最大值中选择更新的非椭圆率各向异性参数；

将更新的非椭圆率各向异性参数存储到存储资源中；

计算更新的非椭圆率各向异性参数和初始非椭圆率各向异性参数之间的扰动变化值；

输入用户定义的公差值；

将所存储的更新的非椭圆率各向异性参数设置为等于初始非椭圆率各向异性参数；

重复如下步骤：计算各向异性基尔霍夫时间偏移局部图像；生成具有水平偏移轴和垂直时间轴的远偏移距图像道集；将所生成的具有水平偏移轴和垂直时间轴的远偏移距图像道集存储到存储资源中；计算多个相似度参数；生成最终相似度体积；存储最终相似度体积；选择更新的非椭圆率各向异性参数；存储更新的非椭圆率各向异性参数；计算扰动变化值，直到所计算出的扰动变化值小于或等于所输入的用户定义的公差值，或者直到所生成的远偏移距图像道集平坦；以及

生成最终非椭圆率各向异性参数。

2.根据权利要求1所述的方法，其特征在于，使用所生成的远偏移距图像道集来计算多个相似度参数还包括步骤：

从所生成的远偏移距图像道集中选择对应于所获取的初始非椭圆率各向异性参数和所获取的基线正常时差速度的远偏移距图像道集；

从所有生成的远偏移距图像道集中根据下式来计算所选择的远偏移距图像道集的偏移时间的算法，作为所获取的初始非椭圆率各向异性参数的函数：

其中，t(x)是偏移时间，t₀是零偏移处图像点的偏移时间，x是震源和接收器之间的偏移，V_nmo是有效正常时差速度，η₀是初始的有效非椭圆率参数，η是更新的有效非椭圆率参数；

使用所计算出的远偏移距图像道集的偏移时间根据下式来计算相似度模型：

其中，q是偏移图像道集中的道振幅，l是用于平滑相似度谱估计的窗口长度；和

重复选择远偏移距图像道集的步骤，计算所选择的远偏移距图像道集的偏移时间的算法的步骤，以及计算相似度模型的步骤，直到从所存储的远偏移距图像道集中的所有远偏移距图像道集都已经被选择了为止。

3.根据权利要求2所述的方法，其特征在于，所述相似度S通过扫描从-0.2到0.2的η值来生成。

4.一种用于储层表征和解释的方法的计算系统，所述计算系统包括：

设置在勘探区域上的包括多条炮点线和多条接收器线的网格，其中，每条炮点线上均布置有多个炮点，每条接收器线上均布置有多个接收器；

用于发送和接收来自勘探区域的叠前地震数据的遥测系统；

用于存储数据的存储资源；

计算机系统输出设备；和

与遥测系统、存储资源和计算机系统输出设备耦合的非暂时性计算机可读存储设备，其被编程为执行以下操作：

获取叠前地震数据，其中所述叠前地震数据包括具有水平偏移轴和垂直时间轴的共炮点道集；

使用基于长偏移距P波数据的常规非椭圆率各向异性参数分析来预处理初始非椭圆率各向异性参数；

使用基于近偏移距P波数据的常规双曲线时差分析或者基于长偏移距P波数据的预偏移非双曲线时差分析来对基线正常时差速度进行预处理；

将所获取的叠前地震数据、所获取的初始非椭圆率各向异性参数和所获取的基线正常时差速度存储到存储资源中；

使用所存储的叠前地震数据、所存储的初始非椭圆率各向异性参数和所存储的基线正常时差来计算各向异性基尔霍夫时间偏移局部图像；

从所计算出的各向异性基尔霍夫时间偏移局部图像中生成具有水平偏移轴和垂直时间轴的远偏移距图像道集；

将所生成的具有水平偏移轴和垂直时间轴的远偏移距图像道集存储到存储资源中；

使用所存储的远偏移距图像道集来计算多个相似度参数；

从所计算出的多个相似度参数中生成最终相似度体积；

将最终相似度体积存储到存储资源中；

从所存储的最终相似度体积的最大值中选择更新的非椭圆率各向异性参数；

将更新的非椭圆率各向异性参数存储到存储资源中；

计算更新的非椭圆率各向异性参数和初始非椭圆率各向异性参数之间的扰动变化值；

输入用户定义的公差值；

将所存储的更新的非椭圆率各向异性参数设置为等于初始非椭圆率各向异性参数；

重复如下步骤：计算各向异性基尔霍夫时间偏移局部图像；生成具有水平偏移轴和垂直时间轴的远偏移距图像道集；将所生成的具有水平偏移轴和垂直时间轴的远偏移距图像道集存储到存储资源中；计算多个相似度参数；生成最终相似度体积；存储最终相似度体积；选择更新的非椭圆率各向异性参数；存储更新的非椭圆率各向异性参数；计算扰动变化值，直到所计算出的扰动变化值小于或等于所输入的用户定义的公差值，或者直到所生成的远偏移距图像道集平坦；以及

生成最终非椭圆率各向异性参数。