获得地震资料的频率衰减梯度的方法和装置.pdf

上传人：奻奴

文档编号：6231337

上传时间：2019-05-23

格式：PDF

页数：9

大小：1.25MB

《获得地震资料的频率衰减梯度的方法和装置.pdf》由会员分享，可在线阅读，更多相关《获得地震资料的频率衰减梯度的方法和装置.pdf（9页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103558634 A (43)申请公布日 2014.02.05 CN 103558634 A (21)申请号 201310520856.7 (22)申请日 2013.10.29 G01V 1/28(2006.01) G01V 1/36(2006.01) (71)申请人中国石油集团川庆钻探工程有限公司地球物理勘探公司地址 610213 四川省成都市双流县华阳镇华阳大道一段一号 1-4 (72)发明人邹文张洞君陶正喜何光明 (74)专利代理机构北京铭硕知识产权代理有限公司 11286 代理人谭昌驰冯敏 (54) 发明名称获得地震资料的频率衰减梯。

2、度的方法和装置 (57) 摘要提供了一种获得地震资料的频率衰减梯度的方法和装置。所述获得地震资料的频率衰减梯度的方法包括：采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱；提取每个预定时间采样点的振幅谱的高频段数据；采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。 (51)Int.Cl. 权利要求书 2 页说明书 5 页附图 1 页 (19)中华人民共和国国家知识产。

3、权局 (12)发明专利申请权利要求书2页说明书5页附图1页 (10)申请公布号 CN 103558634 A CN 103558634 A 1/2 页 2 1. 一种获得地震资料的频率衰减梯度的方法，包括：采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱；提取每个预定时间采样点的振幅谱的高频段数据；采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。 2. 如权利要求 1 。

4、所述的获得地震资料的频率衰减梯度的方法，其中，三参数小波变换方法的三个参数包括小波的调制频率、能量延迟因子和能量衰减因子。 3. 如权利要求 1 所述的获得地震资料的频率衰减梯度的方法，其中，所述计算目标三维地震资料的时间频率域的振幅谱的步骤包括：计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间频率域的振幅谱。 4. 如权利要求 1 所述的获得地震资料的频率衰减梯度的方法，其中，高频段数据是最大振幅至奈奎斯特频率的数据。 5. 如权利要求 1 所述的获得地震资料的频率衰减梯度的方法，其中，指数衰减系数数组的长度与目标。

5、三维地震资料的预定时间采样点的长度相同。 6. 如权利要求 1 所述的获得地震资料的频率衰减梯度的方法，其中，所述去除指数衰减系数数组的背景衰减的步骤包括：采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减。 7. 如权利要求 6 所述的获得地震资料的频率衰减梯度的方法，所述采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减的步骤包括： (1) 设置循环计算的次数； (2) 采用最小二乘法对指数衰减系数数组进行四次多项式拟合以获得拟合的指数衰减系数数组； (3) 计算拟合的指数衰减系数数组与指数衰减系数数组的差值数组，寻找与差。

6、值数组的最大值对应的时间采样点，并采用与所述时间采样点对应的拟合的指数衰减系数数值来更新指数衰减系数数值以获得新的指数衰减数组； (4)针对新的指数衰减数组重复进行步骤(2)和(3)，直到重复次数达到设置的循环次数，以输出最后得到的指数衰减数组；以及 (5) 通过将指数衰减系数数组与最后得到的指数衰减数组进行差值计算来获得频率衰减梯度数组。 8. 一种获得地震资料的频率衰减梯度的装置，包括：振幅谱计算单元，采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱；高频段数据提取单元，提取每个预定时间采样点的振幅谱的高频段数据；曲线拟合计算单元。

7、，采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及权利要求书 CN 103558634 A 2 2/2 页 3 背景衰减去除单元，去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。权利要求书 CN 103558634 A 3 1/5 页 4 获得地震资料的频率衰减梯度的方法和装置技术领域 0001 本发明涉及勘探地球物理领域，更具体地，涉及一种获得地震资料的频率衰减梯度的方法和装置。背景技术 0002 地震波的衰减与砂岩含油气。

8、性存在一定的关系，当地震波穿过含油气砂岩时，能量会发生明显的高频衰减，即，振幅谱的高频段的能量发生明显地衰减。为了有效地利用高频衰减信息，现有技术中通过求取信号的高频段振幅谱的指数衰减系数来表征频率衰减的程度，即，频率衰减梯度分析技术，其中，高频段振幅谱的频率范围指最大振幅对应的频率至计算的最大频率，指数衰减系数越大，地震波的衰减程度越大。因此，常规的频率衰减梯度计算包含 2 个关键步骤：（1）计算信号的时间频率域的振幅谱；（2）计算高频段振幅谱的指数衰减系数。 0003 目前，频率衰减梯度的计算方法主要包括： Mitchell 等（1996。

9、）基于短时傅里叶变换方法和基于指数函数的线性拟合方法进行计算， EP-Tech 软件（2006）的 Image 模块采用三参数小波变换方法和两点线性拟合方法进行计算等等。 0004 目前虽然频率衰减梯度的计算方法有很多种，但是由于 EP-Tech 软件的 Image 模块的商业化程度高而被广泛应用于油气储层的预测研究。但是， EP-Tech 软件的 Image 模块由于采用三参数小波变换方法计算时间频率域的振幅谱，具有一定的缺陷：基于常规的连续小波变换得到的信号的时间频率域的振幅谱，采用了时间尺度信息转换，而在这种转换过程中每个尺度所对应的不是一个频率，而是一个。

10、频带，从而出现了能量交叠的情况，在高频时这种交叠更为严重，导致高频段振幅谱能量的衰减。EP-Tech 软件的 Image 模块虽然采用的是三参数小波，相对于常规的连续小波变换而言，其仅对小波函数进行了改进，而核心算法连续小波变换没有变，即，仍然存在高频段振幅谱的能量衰减现象，该衰减并不是油气储层引起的衰减，而是连续小波变换产生的衰减，从而会影响油气储层预测的可靠性。此外， EP-Tech 软件的 Image 模块采用两点线性拟合方法计算高频段振幅谱的衰减系数，其计算精度还有待提高。发明内容 0005 根据本发明的一方面，提供了一种获得地震资料的频率衰减梯度的。

11、方法，包括：采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱；提取每个预定时间采样点的振幅谱的高频段数据；采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。 0006 三参数小波变换方法的三个参数可包括小波的调制频率、能量延迟因子和能量衰减因子。说明书 CN 103558634 A 4 2/5 页 5 0007 所述计算目标三维地震资料的时间频率域的振幅谱的步骤可。

12、包括：计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间频率域的振幅谱。 0008 高频段数据可以是最大振幅至奈奎斯特频率的数据。 0009 指数衰减系数数组的长度可以与目标三维地震资料的预定时间采样点的长度相同。 0010 所述去除指数衰减系数数组的背景衰减的步骤可包括：采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减。 0011 所述采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减的步骤可包括： (1) 设置循环计算的次数； (2) 采用最小二乘法对指数衰减系数数组进行四次多。

13、项式拟合以获得拟合的指数衰减系数数组； (3) 计算拟合的指数衰减系数数组与指数衰减系数数组的差值数组，寻找与差值数组的最大值对应的时间采样点，并采用与所述时间采样点对应的拟合的指数衰减系数数值来更新指数衰减系数数值以获得新的指数衰减数组； (4) 针对新的指数衰减数组重复进行步骤 (2) 和 (3)，直到重复次数达到设置的循环次数，以输出最后得到的指数衰减数组；以及 (5) 通过将指数衰减系数数组与最后得到的指数衰减数组进行差值计算来获得频率衰减梯度数组。 0012 根据本发明的另一方面，提供了一种获得地震资料的频率衰减梯度的装置，包括：振幅谱计算单元，采。

14、用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱；高频段数据提取单元，提取每个预定时间采样点的振幅谱的高频段数据；曲线拟合计算单元，采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及背景衰减去除单元，去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。 0013 将在接下来的描述中部分阐述本发明另外的方面和 / 或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。附图说明 0014 通过下面结合附图进。

15、行的详细描述，本发明的上述和其它目的和特点将会变得更加清楚，其中： 0015 图 1 是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的方法的流程图； 0016 图 2 是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的装置的框图。具体实施方式 0017 现在，详细描述本发明的示例性实施例，其示例在附图中表示，其中，相同的标号始终表示相同的部件。 0018 图 1 是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的方法的流程图。说明书 CN 103558634 A 5 3/5 页 6 0019 参照图 1，在步骤 S110，设置计算。

16、的频率范围和三参数小波的三个参数。这里，更具体地，可设置计算的频率范围为开始计算的频率值 f1至结束计算的频率值 f2，三参数小波的三个参数可包括小波的调制频率、能量延迟因子和能量衰减因子。 0020 在步骤 S120，采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱，其中，所述振幅谱是二维数组。这里，更具体地，可计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间频率域的振幅谱。 0021 更具体地，仅作为示例而非限制，可通过以下等式来进行时频三参数小波变换方法的目标三维地震资料的时间频率域。

17、的振幅谱的计算： 0022 y(f,)=|Ws(,)|， =2f （1） 0023 在等式（1）中， y(f,) 表示时间频率域的振幅谱， |Ws(,)| 表示与 y(f,) 对应的频率振幅谱， f 表示频率，计算范围为（f1-f2），表示时间采样点，计算范围为（1-nt），表示角频率。 0024 0025 在等式（2）中， Ws(,)表示y(f,)的时频三参数小波变换， F(,)表示其小波变换，表示小波基的傅里叶变换，表示尺度， C表示小波函数容许函数，具体见以下等式（3）。 0026 0027 在等式（3）中， () 表示小波母函数的傅里叶变换。 0。

18、028 0029 在等式（4）中，表示小波母函数的小波变换公式，表示依赖于参数、的小波基函数，表示 () 的复共轭。 0030 0031 等式（5）表示连续小波变换定义，连续小波变换被定义为一系列小波与信号 f(t) 的褶积， f(t) 表示地震资料数据。说明书 CN 103558634 A 6 4/5 页 7 0032 0033 0034 0035 0036 等式（6）为三参数母小波表达式， e 为幂指数，等式（7）至等式（9）中的 k()、 p()、 q() 为计算三参数小波的中间系数，当等于 5.33，等于 0.0，等于 0.5 时，三参数小。

19、波就是 Morlet 小波。 0037 总体来说，等式（1）用于计算时频三参数小波变换的振幅谱，等式（2）至等式（4）用于计算时频三参数小波变换，等式（5）用于计算常规的三参数小波变换，等式（6）至等式（9）用于计算三参数小波。 0038 应该了解，除了上述列出的等式之外，本领域技术人员完全可采用其它的方法或算法来计算目标三维地震资料的时间频率域的振幅谱。 0039 在步骤 S130，针对目标三维地震资料的时间频率域的振幅谱，提取每个预定时间采样点的振幅谱（频率范围为 f1至 f2）的高频段数据（该数据的频率范围为预定时间采样点的振幅谱的最大。

20、值对应的频率至 f2），其中，所述预定时间采样点的振幅谱是一维数组，所述高频段数据是最大振幅至奈奎斯特频率的数据。 0040 在步骤 S140，采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组 q(nt)，其中，指数衰减系数数组的长度与目标三维地震资料的预定时间采样点的长度相同。 0041 更具体地，可通过以下等式来进行基于 Nelde-Mead 单纯形算法的曲线拟合计算： 0042 y=be-ax （10） 0043 0044 其中，等式（10）为指。

21、数拟合公式， a、 b 为拟合系数，等式（11）直接求解方程的最小值以获得与当前时间采样点对应的指数衰减系数，为拟合值， yk为实际值， E(a,b)为拟合值与实际值的误差。这里，上述 Nelde-Mead 单纯形算法为现有技术，在此不再赘述。说明书 CN 103558634 A 7 5/5 页 8 0045 在步骤 S150，去除指数衰减系数数组 q(nt) 的背景衰减，以获得频率衰减梯度数组 atn(nt)。这里，仅作为示例，可采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组 q(nt) 的背景衰减，更具体地，步骤 S150 可进一步包括。

22、： 0046 在步骤 S151，设置循环计算的次数 nc。 0047 在步骤 S152，采用最小二乘法对指数衰减系数数组 q(nt) 进行四次多项式拟合以获得拟合的指数衰减系数数组 qn(nt)，仅作为示例，拟合公式如等式 12 所示： 0048 f(x)=c1+c2x+c3x2+c4x3+c5x4 （12） 0049 其中， c1、 c2、 c3、 c4、 c5表示最小二乘法拟合出的系数， f(x) 表示拟合出的曲线，即，可采用等式 12 对指数衰减系数数组 q(nt) 进行四次多项式拟合，单本发明不限于此。 0050 在步骤 S153，计算 qn(nt) 与 q(nt)。

23、的差值数组 c(nt)，寻找与 c(nt) 的最大值对应的时间采样点 tmax，并采用与 tmax对应的拟合的指数衰减系数数值 qn(tmax) 来更新指数衰减系数数值 q(tmax) 以获得新的指数衰减数组 q (nt)。 0051 在步骤 S154，针对新的指数衰减数组 q (nt) 重复进行上述四次多项式拟合和指数衰减数值更新，直到重复次数达到设置的循环次数 nc，以输出最后得到的指数衰减数组 qn(nt) ； 0052 在步骤 S155，通过将指数衰减系数数组 q(nt) 与最后得到的指数衰减数组 qn(nt) 进行差值计算来获得频率衰减梯度数组 atn(nt)，。

24、即， atn(nt)=q(nt)-qn(nt)。 0053 图 2 是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的装置的框图。 0054 如图 2 所示，获得地震资料的频率衰减梯度的装置可包括振幅谱计算单元 210、高频段数据提取单元 220、曲线拟合计算单元 230、背景衰减去除单元 240。 0055 振幅谱计算单元 210 采用时频三参数小波变换方法来计算目标三维地震资料的时间频率域的振幅谱； 0056 高频段数据提取单元 220 提取每个预定时间采样点的振幅谱的高频段数据； 0057 曲线拟合计算单元 230 采用基于 Nelde-Mead 单纯形算法的非。

25、线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及 0058 背景衰减去除单元 240 去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。 0059 根据本发明的示例性实施例，通过采用时频三参数小波变换方法进行目标信号的时间频率域的振幅谱计算，不仅具有三参数小波的灵活性，而且可以获得真实的高频段的振幅谱，有助于提高频率衰减梯度计算的准确性；通过采用基于 Nelde-Mead 单纯形算法的非线性曲线拟合方法进行高频段振幅谱的指数曲线拟合计算，可以精确地拟合高频段振幅的衰减趋势，有助于提高频率衰减梯度计算的精度；通过采用基于多次曲线拟合迭代的背景衰减去除方法进行地震资料的背景衰减的去除，可以有效地提高背景衰减趋势的拟合精度，提高频率衰减梯度计算的可靠性。说明书 CN 103558634 A 8 1/1 页 9 图 1 图 2 说明书附图 CN 103558634 A 9 。

摘要
申请专利号：	CN201310520856.7	申请日：	2013.10.29
公开号：	CN103558634A	公开日：	2014.02.05
当前法律状态：	撤回	有效性：	无权
法律详情：	发明专利申请公布后的视为撤回IPC(主分类):G01V 1/28申请公布日:20140205\|\|\|实质审查的生效IPC(主分类):G01V 1/28申请日:20131029\|\|\|公开
IPC分类号：	G01V1/28; G01V1/36	主分类号：	G01V1/28
申请人：	中国石油集团川庆钻探工程有限公司地球物理勘探公司
发明人：	邹文; 张洞君; 陶正喜; 何光明
地址：	610213 四川省成都市双流县华阳镇华阳大道一段一号1-4
优先权：
专利代理机构：	北京铭硕知识产权代理有限公司 11286	代理人：	谭昌驰;冯敏
PDF完整版下载：	PDF下载

内容摘要

提供了一种获得地震资料的频率衰减梯度的方法和装置。所述获得地震资料的频率衰减梯度的方法包括：采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；提取每个预定时间采样点的振幅谱的高频段数据；采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。

权利要求书

权利要求书
1.  一种获得地震资料的频率衰减梯度的方法，包括：
采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；
提取每个预定时间采样点的振幅谱的高频段数据；
采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及
去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。

2.  如权利要求1所述的获得地震资料的频率衰减梯度的方法，其中，三参数小波变换方法的三个参数包括小波的调制频率、能量延迟因子和能量衰减因子。

3.  如权利要求1所述的获得地震资料的频率衰减梯度的方法，其中，所述计算目标三维地震资料的时间—频率域的振幅谱的步骤包括：
计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间—频率域的振幅谱。

4.  如权利要求1所述的获得地震资料的频率衰减梯度的方法，其中，高频段数据是最大振幅至奈奎斯特频率的数据。

5.  如权利要求1所述的获得地震资料的频率衰减梯度的方法，其中，指数衰减系数数组的长度与目标三维地震资料的预定时间采样点的长度相同。

6.  如权利要求1所述的获得地震资料的频率衰减梯度的方法，其中，所述去除指数衰减系数数组的背景衰减的步骤包括：
采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减。

7.  如权利要求6所述的获得地震资料的频率衰减梯度的方法，所述采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减的步骤包括：
(1)设置循环计算的次数；
(2)采用最小二乘法对指数衰减系数数组进行四次多项式拟合以获得拟合的指数衰减系数数组；
(3)计算拟合的指数衰减系数数组与指数衰减系数数组的差值数组，寻找与差值数组的最大值对应的时间采样点，并采用与所述时间采样点对应的拟合的指数衰减系数数值来更新指数衰减系数数值以获得新的指数衰减数组；
(4)针对新的指数衰减数组重复进行步骤(2)和(3)，直到重复次数达到设置的循环次数，以输出最后得到的指数衰减数组；以及
(5)通过将指数衰减系数数组与最后得到的指数衰减数组进行差值计算来获得频率衰减梯度数组。

8.  一种获得地震资料的频率衰减梯度的装置，包括：
振幅谱计算单元，采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；
高频段数据提取单元，提取每个预定时间采样点的振幅谱的高频段数据；
曲线拟合计算单元，采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及
背景衰减去除单元，去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。

说明书

说明书获得地震资料的频率衰减梯度的方法和装置
技术领域
本发明涉及勘探地球物理领域，更具体地，涉及一种获得地震资料的频率衰减梯度的方法和装置。
背景技术
地震波的衰减与砂岩含油气性存在一定的关系，当地震波穿过含油气砂岩时，能量会发生明显的高频衰减，即，振幅谱的高频段的能量发生明显地衰减。为了有效地利用高频衰减信息，现有技术中通过求取信号的高频段振幅谱的指数衰减系数来表征频率衰减的程度，即，频率衰减梯度分析技术，其中，高频段振幅谱的频率范围指最大振幅对应的频率至计算的最大频率，指数衰减系数越大，地震波的衰减程度越大。因此，常规的频率衰减梯度计算包含2个关键步骤：（1）计算信号的时间—频率域的振幅谱；（2）计算高频段振幅谱的指数衰减系数。
目前，频率衰减梯度的计算方法主要包括：Mitchell等（1996）基于短时傅里叶变换方法和基于指数函数的线性拟合方法进行计算，EP-Tech软件（2006）的Image模块采用三参数小波变换方法和两点线性拟合方法进行计算等等。
目前虽然频率衰减梯度的计算方法有很多种，但是由于EP-Tech软件的Image模块的商业化程度高而被广泛应用于油气储层的预测研究。但是，EP-Tech软件的Image模块由于采用三参数小波变换方法计算时间—频率域的振幅谱，具有一定的缺陷：基于常规的连续小波变换得到的信号的时间—频率域的振幅谱，采用了时间—尺度信息转换，而在这种转换过程中每个尺度所对应的不是一个频率，而是一个频带，从而出现了能量交叠的情况，在高频时这种交叠更为严重，导致高频段振幅谱能量的衰减。EP-Tech软件的Image模块虽然采用的是三参数小波，相对于常规的连续小波变换而言，其仅对小波函数进行了改进，而核心算法——连续小波变换没有变，即，仍然存在高频段振幅谱的能量衰减现象，该衰减并不是油气储层引起的衰减，而是连续小波变换产生的衰减，从而会影响油气储层预测的可靠性。此外，EP-Tech软件的Image模块采用两点线性拟合方法计算高频段振幅谱的衰减系数，其计算精度还有待提高。
发明内容
根据本发明的一方面，提供了一种获得地震资料的频率衰减梯度的方法，包括：采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；提取每个预定时间采样点的振幅谱的高频段数据；采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。
三参数小波变换方法的三个参数可包括小波的调制频率、能量延迟因子和能量衰减因子。
所述计算目标三维地震资料的时间—频率域的振幅谱的步骤可包括：计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间—频率域的振幅谱。
高频段数据可以是最大振幅至奈奎斯特频率的数据。
指数衰减系数数组的长度可以与目标三维地震资料的预定时间采样点的长度相同。
所述去除指数衰减系数数组的背景衰减的步骤可包括：采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减。
所述采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组的背景衰减的步骤可包括：(1)设置循环计算的次数；(2)采用最小二乘法对指数衰减系数数组进行四次多项式拟合以获得拟合的指数衰减系数数组；(3)计算拟合的指数衰减系数数组与指数衰减系数数组的差值数组，寻找与差值数组的最大值对应的时间采样点，并采用与所述时间采样点对应的拟合的指数衰减系数数值来更新指数衰减系数数值以获得新的指数衰减数组；(4)针对新的指数衰减数组重复进行步骤(2)和(3)，直到重复次数达到设置的循环次数，以输出最后得到的指数衰减数组；以及(5)通过将指数衰减系数数组与最后得到的指数衰减数组进行差值计算来获得频率衰减梯度数组。
根据本发明的另一方面，提供了一种获得地震资料的频率衰减梯度的装置，包括：振幅谱计算单元，采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；高频段数据提取单元，提取每个预定时间采样点的振幅谱的高频段数据；曲线拟合计算单元，采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及背景衰减去除单元，去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。
将在接下来的描述中部分阐述本发明另外的方面和/或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。
附图说明
通过下面结合附图进行的详细描述，本发明的上述和其它目的和特点将会变得更加清楚，其中：
图1是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的方法的流程图；
图2是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的装置的框图。
具体实施方式
现在，详细描述本发明的示例性实施例，其示例在附图中表示，其中，相同的标号始终表示相同的部件。
图1是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的方法的流程图。
参照图1，在步骤S110，设置计算的频率范围和三参数小波的三个参数。这里，更具体地，可设置计算的频率范围为开始计算的频率值f1至结束计算的频率值f2，三参数小波的三个参数可包括小波的调制频率ε、能量延迟因子β和能量衰减因子τ。
在步骤S120，采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱，其中，所述振幅谱是二维数组。这里，更具体地，可计算目标三维地震资料的每条测线的每个共深度点（CDP）或共中心点（CMP）道数据的时间—频率域的振幅谱。
更具体地，仅作为示例而非限制，可通过以下等式来进行时频三参数小波变换方法的目标三维地震资料的时间—频率域的振幅谱的计算：
y(f,τ)=|Ws(ω,τ)|，ω=2πf        （1）
在等式（1）中，y(f,τ)表示时间—频率域的振幅谱，|Ws(ω,τ)|表示与y(f,τ)对应的频率振幅谱，f表示频率，计算范围为（f1-f2），τ表示时间采样点，计算范围为（1-nt），ω表示角频率。
Ws(ω,τ)=1Cψ&Integral;-∞+∞Fω(σ,τ)ψ^(σω)e-iωtσdσσ5/2---(2)]]>
在等式（2）中，Ws(ω,τ)表示y(f,τ)的时频三参数小波变换，Fω(σ,τ)表示其小波变换，表示小波基的傅里叶变换，σ表示尺度，Cψ表示小波函数容许函数，具体见以下等式（3）。
Cψ=2π&Integral;-∞∞|ψ^(ω)|2ωdω<∞---(3)]]>
在等式（3）中，ψ(ω)表示小波母函数的傅里叶变换。
ψ^(σω)e-iωtσ=&Integral;-∞∞ψ(t-τσ)e-iωtdt---(4)]]>
在等式（4）中，表示小波母函数的小波变换公式，表示依赖于参数σ、τ的小波基函数，表示ψ()的复共轭。
Fω(σ,τ)=<f(t),ψσ,τ(t)>=&Integral;-∞∞f(t)1σψ^(t-τσ)dt---(5)]]>
等式（5）表示连续小波变换定义，连续小波变换被定义为一系列小波与信号f(t)的褶积，f(t)表示地震资料数据。
ψ(t)=e-τ(t-β)2{p(Λ)[cos(ϵt)-k(Λ)]+iq(Λ)sin(ϵt)}---(6)]]>
k(Λ)=e-ϵ2/4τ[cos(βϵ)+iq(Λ)p(Λ)sin(βϵ)]---(7)]]>
p(Λ)=[2τπ]1/4{4[e-ϵ22τ-e-3ϵ28τ]cos2(βϵ)+1-e-ϵ2/2τ}-1/2---(8)]]>
q(Λ)=[2τπ]1/4{4[e-ϵ22τ-e-3ϵ28τ]sin2(βϵ)+1-e-ϵ2/2τ}-1/2---(9)]]>
等式（6）为三参数母小波表达式，e为幂指数，等式（7）至等式（9）中的k(Λ)、p(Λ)、q(Λ)为计算三参数小波的中间系数，当ε等于5.33，β等于0.0，τ等于0.5时，三参数小波就是Morlet小波。
总体来说，等式（1）用于计算时频三参数小波变换的振幅谱，等式（2）至等式（4）用于计算时频三参数小波变换，等式（5）用于计算常规的三参数小波变换，等式（6）至等式（9）用于计算三参数小波。
应该了解，除了上述列出的等式之外，本领域技术人员完全可采用其它的方法或算法来计算目标三维地震资料的时间—频率域的振幅谱。
在步骤S130，针对目标三维地震资料的时间—频率域的振幅谱，提取每个预定时间采样点的振幅谱（频率范围为f1至f2）的高频段数据（该数据的频率范围为预定时间采样点的振幅谱的最大值对应的频率至f2），其中，所述预定时间采样点的振幅谱是一维数组，所述高频段数据是最大振幅至奈奎斯特频率的数据。
在步骤S140，采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组q(nt)，其中，指数衰减系数数组的长度与目标三维地震资料的预定时间采样点的长度相同。
更具体地，可通过以下等式来进行基于Nelde-Mead单纯形算法的曲线拟合计算：
y=be-ax         （10）
E(a,b)=Σk=1N(beaxk-yk)2---(11)]]>
其中，等式（10）为指数拟合公式，a、b为拟合系数，等式（11）直接求解方程的最小值以获得与当前时间采样点对应的指数衰减系数，为拟合值，yk为实际值，E(a,b)为拟合值与实际值的误差。这里，上述Nelde-Mead单纯形算法为现有技术，在此不再赘述。
在步骤S150，去除指数衰减系数数组q(nt)的背景衰减，以获得频率衰减梯度数组atn(nt)。这里，仅作为示例，可采用基于多次曲线拟合迭代的背景衰减去除方法去除指数衰减系数数组q(nt)的背景衰减，更具体地，步骤S150可进一步包括：
在步骤S151，设置循环计算的次数nc。
在步骤S152，采用最小二乘法对指数衰减系数数组q(nt)进行四次多项式拟合以获得拟合的指数衰减系数数组qn(nt)，仅作为示例，拟合公式如等式12所示：
f(x)=c1+c2x+c3x2+c4x3+c5x4          （12）
其中，c1、c2、c3、c4、c5表示最小二乘法拟合出的系数，f(x)表示拟合出的曲线，即，可采用等式12对指数衰减系数数组q(nt)进行四次多项式拟合，单本发明不限于此。
在步骤S153，计算qn(nt)与q(nt)的差值数组c(nt)，寻找与c(nt)的最大值对应的时间采样点tmax，并采用与tmax对应的拟合的指数衰减系数数值qn(tmax)来更新指数衰减系数数值q(tmax)以获得新的指数衰减数组q’(nt)。
在步骤S154，针对新的指数衰减数组q’(nt)重复进行上述四次多项式拟合和指数衰减数值更新，直到重复次数达到设置的循环次数nc，以输出最后得到的指数衰减数组qn(nt)；
在步骤S155，通过将指数衰减系数数组q(nt)与最后得到的指数衰减数组qn(nt)进行差值计算来获得频率衰减梯度数组atn(nt)，即，atn(nt)=q(nt)-qn(nt)。
图2是示出根据本发明示例性实施例的获得地震资料的频率衰减梯度的装置的框图。
如图2所示，获得地震资料的频率衰减梯度的装置可包括振幅谱计算单元210、高频段数据提取单元220、曲线拟合计算单元230、背景衰减去除单元240。
振幅谱计算单元210采用时频三参数小波变换方法来计算目标三维地震资料的时间—频率域的振幅谱；
高频段数据提取单元220提取每个预定时间采样点的振幅谱的高频段数据；
曲线拟合计算单元230采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法对高频段数据进行曲线拟合计算，以获得由与每个预定时间采样点对应的指数衰减系数构成的指数衰减系数数组；以及
背景衰减去除单元240去除指数衰减系数数组的背景衰减，以获得频率衰减梯度数组。
根据本发明的示例性实施例，通过采用时频三参数小波变换方法进行目标信号的时间—频率域的振幅谱计算，不仅具有三参数小波的灵活性，而且可以获得真实的高频段的振幅谱，有助于提高频率衰减梯度计算的准确性；通过采用基于Nelde-Mead单纯形算法的非线性曲线拟合方法进行高频段振幅谱的指数曲线拟合计算，可以精确地拟合高频段振幅的衰减趋势，有助于提高频率衰减梯度计算的精度；通过采用基于多次曲线拟合迭代的背景衰减去除方法进行地震资料的背景衰减的去除，可以有效地提高背景衰减趋势的拟合精度，提高频率衰减梯度计算的可靠性。