CN105259573B - 一种地震沉积体信息的获取方法及装置 - Google Patents
一种地震沉积体信息的获取方法及装置 Download PDFInfo
- Publication number
- CN105259573B CN105259573B CN201510599636.7A CN201510599636A CN105259573B CN 105259573 B CN105259573 B CN 105259573B CN 201510599636 A CN201510599636 A CN 201510599636A CN 105259573 B CN105259573 B CN 105259573B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- frequency
- information
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000013049 sediment Substances 0.000 title abstract 4
- 238000001228 spectrum Methods 0.000 claims abstract description 26
- 230000008021 deposition Effects 0.000 claims abstract description 14
- 238000000605 extraction Methods 0.000 claims description 26
- 230000037237 body shape Effects 0.000 claims description 22
- 239000000284 extract Substances 0.000 claims 2
- 238000005516 engineering process Methods 0.000 abstract description 17
- 208000035126 Facies Diseases 0.000 abstract description 13
- 238000004458 analytical method Methods 0.000 abstract description 12
- 230000009466 transformation Effects 0.000 abstract 1
- 239000004576 sand Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 241001074085 Scophthalmus aquosus Species 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000007491 morphometric analysis Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地震沉积体信息的获取方法及装置，涉及石油勘探技术领域，方法包括：获取地震道数据，并将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据；从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；根据峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示峰值频率数据体的频率信息；根据峰值振幅数据体提取沉积体形态信息，并在二维色标上显示峰值振幅数据体的振幅信息。本发明能够解决现有技术中的利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
Description
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种地震沉积体信息的获取方法及装置。
背景技术
目前，随着地震技术的发展，人们开始利用地震资料来研究地层和沉积相，特别是三维地震资料的出现，为人们利用地震波形和属性研究沉积提供了极大的便利，并相继发展了地震相分析技术和地震沉积学分析技术。
地震相分析技术是通过对研究区目的层的地震反射外形、地震反射的内部结构、顶底接触关系、振幅、视周期、连续性、平滑性以及反射特征的横向变化等特征参数进行分析，得到一张地震相分类图，然后结合研究区的区域沉积特征判断其沉积类型。然而地震相分析技术只能分析百米厚的地层，而且无法区分不同岩性及其厚度。地震沉积学分析技术是通过地震切片上振幅的强弱变化和形态特征分析沉积分布。然而，影响地震波振幅与能量的因素有许多，在表层条件、激发—接收条件不稳定时，如果资料处理中保真度差,时间剖面上反射波振幅或能量的变化不能准确地反映薄层厚度变化。并且这种线性关系只有在一定的厚度、反射系数及波长条件下才能成立。因此，传统的利用地震振幅切片研究沉积相形态和厚度的地震沉积学分析技术往往在应用上会遇到存在多解性的问题。
当前，利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度。
发明内容
本发明实施例提供一种地震沉积体信息的获取方法及装置，以解决现有技术中的利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
为了达到上述目的，本发明采用如下技术方案：
一种地震沉积体信息的获取方法，包括：
获取地震道数据，并将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据；
从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；
根据所述峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示所述峰值频率数据体的频率信息；所述二维色标包括一X轴，所述X轴用于表示所述频率信息；
根据所述峰值振幅数据体提取沉积体形态信息；并在所述二维色标上显示所述峰值振幅数据体的振幅信息；所述二维色标还包括一Y轴，所述Y轴用于表示所述振幅信息。
具体的，所述时窗信息为一窗函数：
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
具体的，所述将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据，包括：
根据公式：
将所述地震道数据进行行变速广义S变换；
其中，h(t)为所述地震道数据；S(τ,f)为所述h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
具体的，根据所述峰值振幅数据体提取沉积体形态信息，包括：
根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于反映所述沉积体形态信息。
具体的，根据所述峰值频率数据体提取沉积体厚度信息，包括：
根据公式：
确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t'；
其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值；fo为雷克子波的峰值频率。
具体的，在一二维色标上显示所述峰值频率数据体的频率信息，包括：
根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息。
具体的，在所述二维色标上显示所述峰值振幅数据体的振幅信息，包括：
根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li；
其中：Li为在Y轴第i个亮度值；NL为颜色个数；所述Li的取值范围为0至1。
一种地震沉积体信息的获取装置，包括：
时频谱数据生成单元，用于获取地震道数据，并将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据；
读取单元，用于从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；
沉积体厚度信息提取单元，用于根据所述峰值频率数据体提取沉积体厚度信息；
频率信息显示单元，用于在一二维色标上显示所述峰值频率数据体的频率信息；所述二维色标包括一X轴，所述X轴用于表示所述频率信息；
沉积体形态信息提取单元，用于根据所述峰值振幅数据体提取沉积体形态信息；
振幅信息显示单元，用于在所述二维色标上显示所述峰值振幅数据体的振幅信息；所述二维色标还包括一Y轴，所述Y轴用于表示所述振幅信息。
具体的，所述时频谱数据生成单元中的时窗信息为一窗函数：
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
另外，所述时频谱数据生成单元，具体用于：
根据公式：
将所述地震道数据进行行变速广义S变换；
其中，h(t)为所述地震道数据；S(τ,f)为所述h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
此外，所述沉积体形态信息提取单元，具体用于：
根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于反映所述沉积体形态信息。
另外，所述沉积体厚度信息提取单元，具体用于：
根据公式：
确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t'；
其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值；fo为雷克子波的峰值频率。
另外，所述频率信息显示单元，具体用于：
根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息。
此外，所述振幅信息显示单元，具体用于：
根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li；
其中：Li为在Y轴第i个亮度值；NL为颜色个数；所述Li的取值范围为0至1。
本发明实施例提供的一种地震沉积体信息的获取方法及装置，能够获取地震道数据，从而能够将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据；之后，从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；根据峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示峰值频率数据体的频率信息；根据峰值振幅数据体提取沉积体形态信息，并在二维色标上显示峰值振幅数据体的振幅信息；可见通过地震道数据的变速广义S变换，并进行沉积体信息提取，进而在二维色标上显示频率信息和振幅信息，可以实现同时准确确定沉积形态和沉积厚度信息，且同时显示频率信息和振幅信息的效果，这样在同一切片上既能反映沉积形态，又能说明沉积体厚度。本发明能够现有技术中的利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种地震沉积体信息的获取方法的流程图；
图2为本发明实施例中的一个J1井柱状示意图；
图3为本发明实施例中的J1井对应井旁地震道的变速广义S变换的时频谱的示意图；
图4为本发明实施例中的过J1井的地震剖面的示意图；
图5为本发明实施例中的对应S砂岩的峰值振幅切片的示意图；
图6为本发明实施例中的对应S砂岩的峰值频率切片的示意图；
图7为本发明实施例中的采用二维色标显示的组合图像；
图8为本发明实施例提供的一种地震沉积体信息的获取装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供的一种地震沉积体信息的获取方法，包括：
步骤101、获取地震道数据，并将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据。
步骤102、从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体。
步骤103、根据峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示峰值频率数据体的频率信息。
其中，二维色标包括一X轴，X轴用于表示频率信息。
步骤104、根据峰值振幅数据体提取沉积体形态信息，并在二维色标上显示峰值振幅数据体的振幅信息。
其中，二维色标还包括一Y轴，Y轴用于表示振幅信息。
本发明实施例提供的一种地震沉积体信息的获取方法，能够获取地震道数据，从而能够将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据；之后，从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；根据峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示峰值频率数据体的频率信息；根据峰值振幅数据体提取沉积体形态信息，并在二维色标上显示峰值振幅数据体的振幅信息；可见通过地震道数据的变速广义S变换，并进行沉积体信息提取，进而在二维色标上显示频率信息和振幅信息，可以实现同时准确确定沉积形态和沉积厚度信息，且同时显示频率信息和振幅信息的效果，这样在同一切片上既能反映沉积形态，又能说明沉积体厚度。本发明能够现有技术中的利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
在一实施例中，上述步骤101中的将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据，可以是：
设函数h(t)∈L2(R)，其中h(t)为地震道数据；L2(R)表示能量有限函数空间。这样，h(t)的S变换可以定义为：
在该S变换中，窗函数为：
该S变换的窗函数满足条件窗函数的宽度与频率成反比，低频端的时窗较宽，频率分辨率较高但时间分辨率很低；高频端时窗较窄，时间分辨率较高但频率分辨率很低；仅具有固定的时频分辨率。
在本实施例中，对S变换的窗函数进行改造，引入λ和p两个调节参数，则步骤101中的时窗信息为一窗函数：
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
则该改造后的窗函数对应的变速广义S变换为：
其中，h(t)为地震道数据；S(τ,f)为h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
则上述步骤101可以根据该公式(4)将地震道数据进行行变速广义S变换。
在实际应用中，可通过调节λ和p为两个调节参数，使在低频端保持频率分辨率的同时，相对得提高时间分辨率；高频端保持时间分辨率的同时，相对提高频率分辨率；解决了普通的S变换在地震道数据处理中的应用限制。通过对地震道数据作变速广义S变换，可以得到地震道数据对应的时频谱数据，为下一步提取沉积体信息提供了数据基础。
另外，上述步骤104中的根据峰值振幅数据体提取沉积体形态信息，可以包括如下过程：
根据峰值振幅数据体，在待测的目的层段提取等时切片。该等时切片用于反映沉积体形态信息。
由于变速广义S变换具有高质量的时频分布和较强的薄层探测能力，因此，和基于常规90°相位数据体切片相比，基于峰值振幅数据体提取的等时切片能反映更丰富的沉积体形态信息。
在一实施例中，峰值频率对反射系数大小和极性变化不敏感，只随着薄层厚度的增加而减小，当薄层顶底反射系数比趋近-1时，峰值频率与薄层厚度呈现良好的线性关系。对于输入子波为零相位雷克子波，反射复合波峰值频率与薄层时间厚度存在如下关系：
其中，薄层顶部反射系数为r1、薄层底部反射系数为r2；fo为雷克子波的峰值频率；fp为峰值频率；t'为沉积体薄层的双程旅行时间厚度。
通过公式(5)，可以公式两边均除以r2 2得到：
其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值。
由此可见，地震记录的峰值频率fp只与薄层的双程旅行时间t'和薄层顶底反射系数比R有关，而与顶底反射系数绝对值大小无关。
在实际应用中，对应具体数据，fo为确定值。且对应薄层沉积，其顶底反射系数比R趋近‐1，薄层时间厚度t'与峰值频率fp增量为良好的线性关系，即地震峰值频率随薄层厚度的减小而增加。因此，与峰值振幅数据体中提取的等时切片相对应，在峰值频率数据体中提取对应层位的等时切片，通过等时切片中峰值频率值的大小可以确定对应沉积体的相对厚度大小。
另外，上述步骤103中的在一二维色标上显示峰值频率数据体的频率信息，可以包括：
根据二维色标的X轴的色相H、饱和度S和亮度L来表示频率信息。此处的二维色标的X轴的色相H、饱和度S和亮度L来表示频率信息与一维色标(X轴)的表示方式相同。
另外，上述步骤104中的在二维色标上显示峰值振幅数据体的振幅信息，可以包括：
根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li；
其中：Li为在Y轴第i个亮度值；NL为颜色个数；Li的取值范围为0至1。
从而，在二维色标的Y轴上，通过色相H、饱和度S和亮度Li可以表示振幅信息。
为了使本领域的技术人员更好的了解本发明实施例，下面列举一个具体实例来对本发明实施例进行说明。
如图2所示的J1井柱状图，以青山口组二三段底部S砂岩为例说明。S砂岩在研究区内发育的时间深度范围在1.4-1.8s，故在此时窗内做变速广义S变换。如图3所示是J1井对应井旁地震道的变速广义S变换的时频谱，时频分辨率基本达到平衡，有效反映原始信息。
为使提取的沉积体信息准确反映目标砂岩S，首先，应在S砂岩上下选择稳定的地震地质等时同相轴，研究区S砂岩所在的青山口组有4个比较连续的同相轴，其中，T2、TX对应最大洪泛面，发育稳定的泥岩段(如图4所示)，是理想的地震地质等时同相轴。以T2、TX为顶底约束，制作对应S砂岩的地层切片，并提取对应的峰值频率和峰值振幅。图5为对应S砂岩的峰值振幅切片，从峰值振幅切片上可以看到，研究区发育三角洲前缘分流河道沉积，河道形态清晰。J1、J2、J3井对应的峰值振幅值相近，分别为2483、2950、2390，井中实际砂岩厚度分别为8.3m、6.1m、6.3m，振幅和砂岩厚度对应关系差。图6为对应S砂岩的峰值频率切片，其河道形态不如峰值振幅切片清晰，J1、J2、J3井对应的峰值频率分别为45hz、60hz、56hz，低频对应厚砂体、高频对应薄砂体，两者有很好的对应关系。
之后，可以采用二维色标将峰值振幅和峰值频率信息显示在同一张地震切片上，不仅可以使我们对研究区三角洲沉积的形态有清晰的认识，也可以对三角洲砂体的厚度有定量的认识，从而帮助我们确定物源方向和进行井位部署。图7是采用二维色标显示的组合图像，色标的X轴对应峰值频率，Y轴对应峰值振幅。在组合图像中，分流河道的形态由峰值振幅控制，砂体厚度由峰值频率表示，其中，J1井对应色标为暗灰色(J1标记所在位置)，说明峰值频率值为45hz，J2井和J3井对应色标为灰白色，说明峰值频率值分别为60hz、56hz。
可见，本发明实施例可以实现沉积体形态和沉积体厚度的显示。解决了利用地震相分析技术和地震沉积学分析技术很难同时确定沉积形态和沉积厚度的难点。在本发明指导下生成的切片(如图7所示)有效识别了沉积形态，准确刻画了河道的平面展布，指出了分流河道砂体发育较厚的区域，预测结果与已有钻井结果较为吻合，如J1井钻遇砂体较厚，实际钻遇砂体8.3m，J2、J3井钻遇砂体稍薄，实际钻遇砂体分别为6.1m、6.3m。本发明的应用结果为以后的钻井部署提供了可靠依据。
对应于上述图1对应的方法实施例，本发明实施例还提供一种地震沉积体信息的获取装置，如图8所示，包括：
时频谱数据生成单元81，可以获取地震道数据，并将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据。
读取单元82，可以从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体。
沉积体厚度信息提取单元83，可以根据峰值频率数据体提取沉积体厚度信息。
频率信息显示单元84，可以在一二维色标上显示峰值频率数据体的频率信息；二维色标包括一X轴，X轴用于表示频率信息。
沉积体形态信息提取单元85，可以根据峰值振幅数据体提取沉积体形态信息。
振幅信息显示单元86，可以在二维色标上显示峰值振幅数据体的振幅信息；二维色标还包括一Y轴，Y轴用于表示振幅信息。
具体的，时频谱数据生成单元81中的时窗信息为一窗函数：
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
另外，时频谱数据生成单元81，具体可以根据公式：
将地震道数据进行行变速广义S变换。
其中，h(t)为地震道数据；S(τ,f)为h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
此外，沉积体形态信息提取单元85，具体可以根据峰值振幅数据体，在待测的目的层段提取等时切片；等时切片用于反映沉积体形态信息。
另外，沉积体厚度信息提取单元83，具体可以根据公式：
确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t'。其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值；fo为雷克子波的峰值频率。
另外，该频率信息显示单元84，具体可以根据二维色标的X轴的色相H、饱和度S和亮度L来表示频率信息。
此外，振幅信息显示单元86，具体可以根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li。
其中：Li为在Y轴第i个亮度值；NL为颜色个数；Li的取值范围为0至1。
值得说明的是，本发明实施例提供的地震沉积体信息的获取装置的具体实现方式可以参见上述图1对应的方法实施例，此处不再赘述。
本发明实施例提供的一种地震沉积体信息的获取装置，能够获取地震道数据，从而能够将地震道数据根据一时窗信息进行变速广义S变换，生成地震道数据对应的时频谱数据；之后，从时频谱数据中读取峰值频率和峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；根据峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示峰值频率数据体的频率信息；根据峰值振幅数据体提取沉积体形态信息，并在二维色标上显示峰值振幅数据体的振幅信息；可见通过地震道数据的变速广义S变换，并进行沉积体信息提取，进而在二维色标上显示频率信息和振幅信息，可以实现同时准确确定沉积形态和沉积厚度信息，且同时显示频率信息和振幅信息的效果，这样在同一切片上既能反映沉积形态，又能说明沉积体厚度。本发明能够现有技术中的利用地震相分析技术和地震沉积学分析技术难以同时准确确定沉积形态和沉积厚度的问题。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种地震沉积体信息的获取方法，其特征在于，包括：
获取地震道数据，并将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据；
从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；
根据所述峰值频率数据体提取沉积体厚度信息，并在一二维色标上显示所述峰值频率数据体的频率信息；所述二维色标包括一X轴，所述X轴用于表示所述频率信息；其中，根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息；
根据所述峰值振幅数据体提取沉积体形态信息；并在所述二维色标上显示所述峰值振幅数据体的振幅信息；所述二维色标还包括一Y轴，所述Y轴用于表示所述振幅信息；具体地，根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li；其中：Li为在Y轴第i个亮度值；NL为颜色个数；所述Li的取值范围为0至1。
2.根据权利要求1所述的地震沉积体信息的获取方法，其特征在于，所述时窗信息为一窗函数：
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mo>|</mo>
<mi>f</mi>
<msup>
<mo>|</mo>
<mi>p</mi>
</msup>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<msup>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>p</mi>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
3.根据权利要求2所述的地震沉积体信息的获取方法，其特征在于，所述将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据，包括：
根据公式：
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mo>|</mo>
<mi>f</mi>
<msup>
<mo>|</mo>
<mi>p</mi>
</msup>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>p</mi>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
将所述地震道数据进行行变速广义S变换；
其中，h(t)为所述地震道数据；S(τ,f)为所述h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
4.根据权利要求3所述的地震沉积体信息的获取方法，其特征在于，根据所述峰值振幅数据体提取沉积体形态信息，包括：
根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于反映所述沉积体形态信息。
5.根据权利要求4所述的地震沉积体信息的获取方法，其特征在于，根据所述峰值频率数据体提取沉积体厚度信息，包括：
根据公式：
<mrow>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>&lsqb;</mo>
<msup>
<mi>R&pi;t</mi>
<mo>&prime;</mo>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>p</mi>
</msub>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mo>&lsqb;</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>p</mi>
</msub>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t'；
其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值；fo为雷克子波的峰值频率。
6.一种地震沉积体信息的获取装置，其特征在于，包括：
时频谱数据生成单元，用于获取地震道数据，并将所述地震道数据根据一时窗信息进行变速广义S变换，生成所述地震道数据对应的时频谱数据；
读取单元，用于从所述时频谱数据中读取峰值频率和所述峰值频率对应的峰值振幅，分别形成峰值频率数据体和峰值振幅数据体；
沉积体厚度信息提取单元，用于根据所述峰值频率数据体提取沉积体厚度信息；
频率信息显示单元，用于在一二维色标上显示所述峰值频率数据体的频率信息；所述二维色标包括一X轴，所述X轴用于表示所述频率信息；其中，根据二维色标的X轴的色相H、饱和度S和亮度L来表示所述频率信息；
沉积体形态信息提取单元，用于根据所述峰值振幅数据体提取沉积体形态信息；
振幅信息显示单元，用于在所述二维色标上显示所述峰值振幅数据体的振幅信息；所述二维色标还包括一Y轴，所述Y轴用于表示所述振幅信息；具体地，根据X轴的亮度值L，确定二维色标的Y轴的亮度值Li；其中：Li为在Y轴第i个亮度值；NL为颜色个数；所述Li的取值范围为0至1。
7.根据权利要求6所述的地震沉积体信息的获取装置，其特征在于，所述时频谱数据生成单元中的时窗信息为一窗函数：
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mo>|</mo>
<mi>f</mi>
<msup>
<mo>|</mo>
<mi>p</mi>
</msup>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<msup>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>p</mi>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，t为窗函数表示的时间；f为频率；λ和p为两个调节参数。
8.根据权利要求7所述的地震沉积体信息的获取装置，其特征在于，所述时频谱数据生成单元，具体用于：
根据公式：
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mo>|</mo>
<mi>f</mi>
<msup>
<mo>|</mo>
<mi>p</mi>
</msup>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>p</mi>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
将所述地震道数据进行行变速广义S变换；
其中，h(t)为所述地震道数据；S(τ,f)为所述h(t)的时频谱；τ和f分别表示时频谱中的时间和频率。
9.根据权利要求8所述的地震沉积体信息的获取装置，其特征在于，所述沉积体形态信息提取单元，具体用于：
根据所述峰值振幅数据体，在待测的目的层段提取等时切片；所述等时切片用于反映所述沉积体形态信息。
10.根据权利要求9所述的地震沉积体信息的获取装置，其特征在于，所述沉积体厚度信息提取单元，具体用于：
根据公式：
<mrow>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>&lsqb;</mo>
<msup>
<mi>R&pi;t</mi>
<mo>&prime;</mo>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>p</mi>
</msub>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mo>&lsqb;</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>p</mi>
</msub>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
确定峰值频率fp对应的沉积体薄层的双程旅行时间厚度t'；
其中，R为薄层顶部反射系数r1与薄层底部反射系数r2的比值；fo为雷克子波的峰值频率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510599636.7A CN105259573B (zh) | 2015-09-18 | 2015-09-18 | 一种地震沉积体信息的获取方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510599636.7A CN105259573B (zh) | 2015-09-18 | 2015-09-18 | 一种地震沉积体信息的获取方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105259573A CN105259573A (zh) | 2016-01-20 |
| CN105259573B true CN105259573B (zh) | 2017-12-19 |
Family
ID=55099334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510599636.7A Active CN105259573B (zh) | 2015-09-18 | 2015-09-18 | 一种地震沉积体信息的获取方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105259573B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106054247B (zh) * | 2016-05-25 | 2020-09-29 | 中国石油集团东方地球物理勘探有限责任公司 | 基于转换波地震数据求取高精度反射系数的方法 |
| CN109709607B (zh) * | 2018-12-07 | 2020-07-10 | 中国石油天然气股份有限公司 | 一种预测薄层砂岩储层厚度方法及装置 |
| CN110456415A (zh) * | 2019-07-17 | 2019-11-15 | 中国石油大港油田勘探开发研究院 | 一种基于局部峰值频率的陆相薄储层解释方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100165791A1 (en) * | 2008-09-05 | 2010-07-01 | Statoilhydro Asa | Method for quantitatively making a thickness estimate of thin geological layers based on seismic reflection signals in the frequency domain |
| US20130144534A1 (en) * | 2011-12-03 | 2013-06-06 | Geokinetics Acquisition Company | Multi-Component Spectral Warping Analysis for Thin Bed Resolution |
| CN103454676B (zh) * | 2013-08-09 | 2016-12-28 | 中国石油天然气股份有限公司 | 一种非线性等时地层切片的确定方法及系统 |
| CN103454678B (zh) * | 2013-08-12 | 2016-02-10 | 中国石油天然气股份有限公司 | 一种地震切片等时性的确定方法及系统 |
-
2015
- 2015-09-18 CN CN201510599636.7A patent/CN105259573B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105259573A (zh) | 2016-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101551466B (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
| CN104570067B (zh) | 一种地球物理勘探中相控地震反演方法 | |
| Tsuji et al. | VP/VS ratio and shear-wave splitting in the Nankai Trough seismogenic zone: Insights into effective stress, pore pressure, and sediment consolidation | |
| CN106054248A (zh) | 一种基于大面积致密储层地震岩石物理反演方法 | |
| Pugin et al. | Multicomponent high-resolution seismic reflection profiling | |
| CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
| CN104316965B (zh) | 一种裂缝方位和强度的预测方法及系统 | |
| CN107728204A (zh) | 基于叠前纵波各向异性的裂缝预测方法及系统 | |
| CN101276001B (zh) | 地下非均匀介质地震探测复杂性定量评估方法 | |
| CN102866426B (zh) | 一种利用avo大角度道集分析岩体油气信息的方法 | |
| CN105259573B (zh) | 一种地震沉积体信息的获取方法及装置 | |
| Odum et al. | VS 30 and Spectral Response from Collocated Shallow, Active‐, and Passive‐Source VS Data at 27 Sites in Puerto Rico | |
| CN107728205B (zh) | 一种地层压力预测方法 | |
| Yuan et al. | 6D phase-difference attributes for wide-azimuth seismic data interpretation | |
| CN104216013B (zh) | 基于宽方位角资料的c3相干体的方法 | |
| CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
| CN104570086A (zh) | 一种在共偏移距共方位角域进行叠前裂缝预测的方法 | |
| Iqbal et al. | Near surface velocity estimation using GPR data: Investigations by numerical simulation, and experimental approach with AVO response | |
| Liu et al. | Combined petrophysics and 3D seismic attributes to predict shale reservoirs favourable areas | |
| Bruno | Seismic Exploration Methods for Structural Studies and for Active Fault Characterization: A Review | |
| Li et al. | Seismic wave field anomaly identification of ultra-deep heterogeneous fractured-vuggy reservoirs: A case study in Tarim Basin, China | |
| CN105093300B (zh) | 一种地质体边界识别方法及装置 | |
| CN102914797A (zh) | 一种获得地层各向异性系数的方法及装置 | |
| Nardone et al. | Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations | |
| Li et al. | SV-P extraction and imaging for far-offset vertical seismic profile data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |