专利名称：基于地震子波衰减谱的油气检测方法及装置的制作方法
技术领域：
本发明涉及石油地球物理勘探中的油气检测领域，特别是涉及一种基于地震子波衰减谱的油气检测方法及装置。
背景技术：
目前的研究和实际应用表明，当地质体中的孔隙充填有油气时，地震反射吸收加大，地震高频能量吸收衰减加剧，含油气地层吸收系数可比相同岩性不含油气地层的吸收系数高出一个甚至几个数量级。基于上述原理，Mitchell等人(Mitchell J Τ, Derzhi N， Lichma Ε. Energy absorption analysis :A case study. Expanded Abstracts of 66th SEG Mtg, 1996)提出一种计算地震信号能量衰减的分析方法(即EAA技术)，来进行地质体的油气检测。该方法的核心是求解信号谱的高频指数衰减系数，具体计算时以一系列小时窗对地震道连续作谱分析，进而计算衰减系数。发明人在实现本发明的过程中发现，上述方法存在以下三个缺点一、地震信号是由子波与反射系数褶积，再加上噪音组成的，直接对地震道进行谱分析，会受到地震噪音和反射系数的严重影响，导致地震信号谱估算不准，误差较大；二、采用傅里叶变换法计算信号谱，要求较大的时窗，但是针对较小的地质体，过大的时窗会导致计算结果产生较大的误差；三、采用一系列的单时窗对地震信号进行高频衰减分析，导致含油气识别的对比性不强，不够直观。

发明内容
本发明的目的是为了克服上述背景技术的不足，提供一种基于地震子波衰减谱的油气检测方法及装置，能够避免受到地震噪音和反射系数的影响，实现准确、敏感地识别油气，使含油气识别具有较强的对比性，更加直观。一方面，为达上述目的，本发明提供的基于地震子波衰减谱的油气检测方法，包括以下步骤采用独立成分分析方法对地震数据进行去噪处理；对去噪处理后的地震数据进行目的层顶、底地震层位追踪；沿目的层顶、底开上、下两个时窗；对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；根据预设的设计低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；对上、下两个时窗内的地震子波谱进行差值计算，得到上、下两个时窗内的地震子
5波衰减谱；通过分析地震子波衰减谱来检测地层的含油气性。另一方面，为达上述目的，本发明实施例还提供了一种基于地震子波衰减谱的油气检测装置，所述装置包括去噪处理单元，用于采用独立成分分析方法对地震数据进行去噪处理；地震层位追踪单元，用于对去噪处理后的地震数据进行目的层顶、底地震层位追踪；时窗选取单元，用于沿目的层顶、底开上、下两个时窗；频谱分析单元，用于对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；复倒谱分析单元，用于对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；复倒谱分离单元，用于根据预设的低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；复倒谱逆变换单元，用于对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；地震子波衰减谱获取单元，用于对上、下两个时窗内的地震子波谱进行差值计算， 得到上、下两个时窗内的地震子波衰减谱；地层含油气性检测单元，用于通过分析地震子波衰减谱来检测地层的含油气性。本发明实施例的有益效果在于本发明实施例的方法和装置，能够克服地震噪音和反射系数的影响，并且通过上下时窗地震子波谱差值处理，使得含油气识别更加直观，具有较强对比性，能够实现准确、敏感地识别油气，为地质目标的勘探开发提供支持。


为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本发明实施例的油气检测方法的整体流程图2为本发明实施例的油气检测方法的具体流程图3为本发明实施例的沿目的层顶、底开上、下两个时窗的示意图4为本发明实施例的上、下时窗地震子波衰减谱；
图5为本发明实施例的地震频谱能量图6为本发明实施例的4 口单井井点处地震道上下时窗地震子波频谱图
图7本发明实施例地的震子波衰减谱切片分析示意图8本发明实施例地的震子波衰减谱切片平面图9为本发明实施例的油气检测装置的整体功能框图10为本发明实施例的油气检测装置的具体功能框图。
具体实施例方式为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明实施例中，首先对地震数据进行去噪处理，尽可能地消除地震噪音影响，然后求取沿目的层的上、下时窗的地震子波衰减谱，再通过地震子波衰减谱来检测地层的含油气性。图1为本发明实施例的油气检测方法的整体流程图。如图1所示，该方法包括如下步骤101、采用独立成分分析方法对地震数据进行去噪处理；102、对去噪处理后的地震数据进行目的层顶、底地震层位追踪；103、沿目的层顶、底开上、下两个时窗；104、对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；105、对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；106、根据预设的低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；107、对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；108、对上、下两个时窗内的地震子波谱进行差值计算，得到上、下两个时窗内的地
震子波衰减谱；109、通过分析所述地震子波衰减谱来检测地层的含油气性。本发明实施例图1所示的方法，能够克服地震噪音和反射系数的影响，并且通过上下时窗地震子波谱差值处理，使得含油气识别更加直观，具有较强对比性，能够实现准确、敏感地识别油气，为地质目标的勘探开发提供支持。图2为本发明实施例的油气检测方法的具体流程图，如图2所示，本发明实施例提供的基于地震子波衰减谱的油气检测方法，具体包括以下步骤(1)地震数据准备。(2)去噪处理采用独立成分分析方法对地震数据进行去噪处理，以提高地震数据的信噪比，尽可能地消除地震噪音影响。(3)测井数据准备。(4)合成记录测井数据与地震数据结合开展合成记录工作，把深度域的地质层位标定到时间域的地震剖面上。(5)目的层顶、底层位追踪在合成记录制作的基础上，开展准确的目的层顶、底地震层位追踪。(6)选择已知井点的联井剖面在三维地震数据体上提取已知井的联井数据剖面。三位地震数据体是在工作区进行野外数据采集，对采集进行数据处理而得到的。(7)沿目的层顶、底开上、下时窗在已知井的联井数据剖面上，沿目的层的顶、底开上、下两个时窗，上、下时窗大小都为t0，参见图3所示，其中图3为本发明实施例沿目的层顶、底开上、下两个时窗的示意图。(8)对上、下时窗地震数据进行小波变换分频应用小波变换分频技术，对上、下时窗内的地震数据进行频谱分析，得到上、下时窗内地震数据的振幅谱分布。(9)进行复倒谱转换对上、下时窗内的地震数据振幅谱进行复倒谱转换，得到上、下时窗地震数据振幅谱的复倒谱。(10)设计低通时窗地震子波的振幅谱是光滑的，其复倒谱集中分布在原点
附近(分布范围Fl F2)，而反射系数的振幅谱是震荡的，其复倒谱分布在远离原点
处。设计一个能够分离地震子波振幅谱的低通时窗。低通时窗具体设计做法建立函数 A,F\<f<F2;
F( f) = \将F(f)与复倒谱相乘。复倒谱分布范围是Fl F2。F1，F2指
地震子波振幅谱的复倒谱的分布范围。(11)分离出上、下时窗地震子波复倒谱在对上、下时窗地震数据振幅谱进行复倒谱变换的基础上，结合该低通时窗，分离出上、下时窗地震子波振幅谱的复倒谱。(12)获得上、下时窗地震子波谱对上、下时窗地震子波振幅谱的复倒谱进行逆变换，得到上、下时窗地震子波谱。(13)获得上、下时窗地震子波衰减谱对上、下两个时窗内地震子波谱进行差值计算，得到上、下两个时窗内的地震子波衰减谱。(14)对比已知钻井与地震子波衰减谱地震子波衰减谱是用变密度显示的，如图 4所示，图4为上、下时窗地震子波衰减谱，图4左上图是地震剖面，横坐标代表cdp (common depth point，共深度点)位置(每个cdp对应一个大地坐标位置)，纵坐标代表时间；图4 左下是地震子波衰减谱，横坐标代表cdp位置，纵坐标代表频率(0-100赫兹)。图4右是地震子波衰减谱的色标，不同灰度代表不同程度的衰减，深色代表衰减最大。在地震子波衰减谱图中(图4左下)，可以看到在cdp79附近(yl2井点处，黑线所示)在48赫兹左右存在一个明显的衰减(深色团块处)。从图5示出的地震频谱能量图中可以看出地震资料主频 25赫兹左右，40-60赫兹属于地震频谱的高频段。也就是说在cdp79位置附近，地震子波通过目的层段，在高频段存在一个明显的衰减(必须是高频段存在衰减)，说明该处含油气性的可能性很高。在一实施例中是通过各口单井的统计，来确定衰减谱大于预设的衰减能量值的区域认为是可能的含油气区。具体地，对井点处上、下时窗地震子波谱进行统计分析， 如图6所示，横坐标代表频率，纵坐标代表归一化能量，在水井和干井(y203、y3井)处，在高频段(40-60hz)上、下时窗子波谱在高频段微弱衰减，在油气井处(y2、y201井)，上、下时窗子波谱在高频段呈现大于0. 36db (归一化后)的衰减。即衰减谱大于0. 36db处为可能的含油气区。其中，db为衰减能量单位。这个量化的值是经过已经钻井统计得出的，不同区域值不同，视实际统计结果而定。(15)选择最优时窗对已知钻井和地震子波衰减谱进行对比分析，如果在油气井处，地震子波衰减谱在高频段出现严重衰减，即衰减程度大于预设的衰减能量值时，说明时窗大小to选择合适。如果已知钻井情况与地震子波衰减谱不匹配，则重新调整时窗大小 to to 士 η* At (η = 1,2,3-, Δ t 是时窗步长)，按照上述(7) (8) (9) (10) (11) (12) (13) (14)步骤重新进行地震子波衰减谱的计算，直到所有已知钻井与地震子波衰减谱匹配，最终确定最合适的时窗大小。(16)全区地震子波衰减谱计算在选择最优时窗基础上，运用(7)⑶(9) (10) (11) (12) (13) (14)步骤，对全区的三维地震数据体进行地震子波衰减谱计算。(17)地震子波衰减谱切片分析地震子波衰减谱是一个三维数据体，图6为地震子波衰减谱切片分析示意图，如图7所示，X，y代表大地坐标，ζ代表的是频率。三维数据体数据可以表示成(x，y，频率，子波衰减谱值)，切片就是固定频率值，提取(x，y，子波衰减谱值)。在三维数据体上，按垂直于频率轴(ζ轴)进行切片，得到不同频率的子波衰减平面图。(18)有利含油气区预测选择高频的地震子波衰减谱切片平面图来进行含油气性的平面预测，衰减严重的区域就是可能的含油气区域。以塔里木盆地塔北地区作为实例说明，图8为地震子波衰减谱切片平面图，如图8所示，y202、y201、y2_4、y204、y2、y206、 y3、y4、y2_hl、yl2井点处地震子波衰减能量大于0. 36db，说明y202、y201、y2_4、y204、 12、y206、y3、y4、y2_hl、yl2井点处是有利的产油气区，实际钻探上述10 口井均产工业油气流，预测结果与实际钻探的结果相符。y203、yl两口井井点处的地震子波衰减能量小于 0. 36db，预测两口井点处不产油气流，实际钻探显示y203井产水，yl井产油，只有yl井的预测结果与实际钻探情况不符。并且可以看到存在西南方向的有利区带(在黑色虚线范围内)。y206字样符号代表钻井井位位置。本发明实施例能够克服地震噪音和反射系数的影响，并且通过上下时窗地震子波谱差值处理，使得含油气识别更加直观，具有较强对比性，能够开展更加准确、敏感的油气检测，特别是对碳酸盐岩裂缝型储层。在塔里木盆地的塔中、塔北地区，碳酸盐岩裂缝型储层的油气检测中，预测结果与实钻结果吻合率达到91. %。本发明实施例还提供了一种基于地震子波衰减谱的油气检测装置。图9为本发明实施例的油气检测装置的整体功能框图，如图9所示，该装置包括去噪处理单元302，用于采用独立成分分析方法对地震数据进行去噪处理；地震层位追踪单元304，用于对去噪处理后的地震数据进行目的层顶、底地震层位追踪；时窗选取单元306，用于沿目的层顶、底开上、下两个时窗；频谱分析单元308，用于对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；复倒谱分析单元310，用于对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；复倒谱分离单元312，用于根据预设的低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；复倒谱逆变换单元314，用于对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；地震子波衰减谱获取单元316，用于对上、下两个时窗内的地震子波谱进行差值计算，得到上、下两个时窗内的地震子波衰减谱；地层含油气性检测单元318，用于通过分析地震子波衰减谱来检测地层的含油气性。
图10为本发明实施例的油气检测装置的具体功能框图，如图10所示，可选地，所述装置还可以包括合成记录制作单元320，用于开展合成记录制作，将地震层位转换到时间域，并标定到地震剖面上。可选地，所述装置还可以包括联井剖面选择单元322，用于选择一条联串已知井的联井剖面。可选地，所述频谱分析单元308，具体可以是采用小波变换分频技术对所述上、下两个时窗内的地震数据进行频谱分析。可选地，所述装置还可以包括时窗大小确定单元324，用于确定最合适的上、下时窗大小。可选地，所述时窗大小确定单元324，具体可以用于将已知钻井情况和地震子波衰减谱进行对比分析，若油气井处的地震子波衰减谱在高频段出现大于预设的衰减能量值的衰减时，则表明上、下时窗大小选择合适。例如，该预设的衰减能量值可以为0. 36db，但本发明实施例不以此为限。可选地，所述时窗大小确定单元324，还可以用于若已知钻井情况与地震子波衰减谱不匹配，则重新调整上、下时窗大小，直到所有已知钻井与地震子波衰减谱匹配，最终确定合适的上、下时窗大小。具体地，所述地层含油气性检测单元318包括频带分析单元326，用于对上、下两个时窗内的地震子波衰减谱进行频带分析，确定所述地震子波衰减谱中的高频段；切片分析单元328，用于对所述地震子波衰减谱的高频段执行切片分析，根据切片分析预测可能含油气区域。可选地，所述切片分析单元3 具体包括切片模块330，用于在地震子波衰减谱对应的三维数据体上，按垂直于频率轴进行切片，得到不同频率的地震子波衰减谱切片平面图；预测模块332，用于选择高频的地震子波衰减谱切片平面图来进行含油气性的平面预测，确定衰减大于预设的衰减能量值的区域为可能的含油气区域。本发明实施例的油气检测装置的工作过程已在前面的方法实施例中详述，故不再赘述。本发明实施例的油气检测装置的优点在于，能够克服地震噪音和反射系数的影响，并且通过上下时窗地震子波谱差值处理，使得含油气识别更加直观，具有较强对比性， 能够开展更加准确、敏感的油气检测，特别是对碳酸盐岩裂缝型储层。在塔里木盆地的塔中、塔北地区，碳酸盐岩裂缝型储层的油气检测中，预测结果与实际钻井结果吻合率达到 91. 7%。本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory, ROM)或随机存储记忆体(Random Access Memory, RAM)等。以上实施例仅用以说明本发明实施例的技术方案，而非对其限制；尽管参照前述
10实施例对本发明实施例进行了详细的说明，本领域的普通技术人员应当理解其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例各实施例技术方案的精神和范围。
权利要求
1.一种基于地震子波衰减谱的油气检测方法，其特征在于，所述方法包括以下步骤采用独立成分分析方法对地震数据进行去噪处理；对去噪处理后的地震数据进行目的层顶、底地震层位追踪；沿目的层顶、底开上、下两个时窗；对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；根据预设的低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；对上、下两个时窗内的地震子波谱进行差值计算，得到上、下两个时窗内的地震子波衰减谱；通过分析所述地震子波衰减谱来检测地层的含油气性。
2.根据权利要求1所述的基于地震子波衰减谱的油气检测方法，其特征在于所述进行目的层顶、底地震层位追踪之前还包括以下步骤开展合成记录制作，将地震层位转换到时间域，并标定到地震剖面上。
3.根据权利要求1所述的基于地震子波衰减谱的油气检测方法，其特征在于在所述沿目的层顶、底开上、下两个时窗之前还包括步骤选择一条联串已知井的联井地震数据剖
4.根据权利要求1所述的基于地震子波衰减谱的油气检测方法，其特征在于所述对上、下两个时窗内的地震数据进行频谱分析，是采用小波变换分频技术。
5.根据权利要求1所述的基于地震子波衰减谱的油气检测方法，其特征在于所述通过分析地震子波衰减谱来检测地层的含油气性包括步骤对上、下两个时窗内的地震子波衰减谱进行频带分析，确定所述地震子波衰减谱中的高频段；对所述地震子波衰减谱的高频段执行切片分析，根据切片分析预测可能含油气区域。
6.根据权利要求5所述的基于地震子波衰减谱的油气检测方法，其特征在于在对所述地震子波衰减谱的高频段执行切片分析之前，还包括步骤确定合适的上、下时窗大小。
7.根据权利要求6所述的基于地震子波衰减谱的油气检测方法，其特征在于所述确定合适的上、下时窗大小包括将已知钻井情况和地震子波衰减谱进行对比分析，若油气井处的地震子波衰减谱在高频段出现大于预设的衰减能量值的衰减时，则表明上、下时窗大小选择合适。
8.根据权利要求7所述的基于地震子波衰减谱的油气检测方法，其特征在于所述确定合适的上、下时窗大小还包括若已知钻井情况与地震子波衰减谱不匹配，则重新调整上、下时窗大小，直到所有已知钻井与地震子波衰减谱匹配，最终确定合适的上、下时窗大
9.根据权利要求5所述的基于地震子波衰减谱的油气检测方法，其特征在于所述对所述地震子波衰减谱的高频段执行切片分析，根据切片分析预测可能含油气区域包括步骤在地震子波衰减谱对应的三维数据体上，按垂直于频率轴进行切片，得到不同频率的地震子波衰减谱切片平面图；选择高频的地震子波衰减谱切片平面图来进行含油气性的平面预测，确定衰减大于预设的衰减能量值的区域为可能的含油气区域。
10.一种基于地震子波衰减谱的油气检测装置，其特征在于，所述装置包括 去噪处理单元，用于采用独立成分分析方法对地震数据进行去噪处理；地震层位追踪单元，用于对去噪处理后的地震数据进行目的层顶、底地震层位追踪； 时窗选取单元，用于沿目的层顶、底开上、下两个时窗；频谱分析单元，用于对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布； 复倒谱分析单元，用于对所述频谱分布进行复倒谱分析，得到上、下时窗地震数据振幅谱的复倒谱；复倒谱分离单元，用于根据预设的低通时窗从所述地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；复倒谱逆变换单元，用于对所述地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；地震子波衰减谱获取单元，用于对上、下两个时窗内的地震子波谱进行差值计算，得到上、下两个时窗内的地震子波衰减谱；地层含油气性检测单元，用于通过分析地震子波衰减谱来检测地层的含油气性。
11.根据权利要求10所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述装置还包括合成记录制作单元，用于开展合成记录制作，将地震层位转换到时间域，并标定到地震剖面上。
12.根据权利要求10所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述装置还包括联井剖面选择单元，用于选择一条联串已知井的联井剖面。
13.根据权利要求10所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述频谱分析单元，具体是采用小波变换分频技术对所述上、下两个时窗内的地震数据进行频谱分析。
14.根据权利要求10所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述地层含油气性检测单元包括频带分析单元，用于对上、下两个时窗内的地震子波衰减谱进行频带分析，确定所述地震子波衰减谱中的高频段；切片分析单元，用于对所述地震子波衰减谱的高频段执行切片分析，根据切片分析预测可能含油气区域。
15.根据权利要求14所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述装置还包括时窗大小确定单元，用于确定合适的上、下时窗大小。
16.根据权利要求15所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述时窗大小确定单元，具体用于将已知钻井情况和地震子波衰减谱进行对比分析，若油气井处的地震子波衰减谱在高频段出现大于预设的衰减能量值的衰减时，则表明上、下时窗大小选择合适。
17.根据权利要求16所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述时窗大小确定单元，还用于若已知钻井情况与地震子波衰减谱不匹配，则重新调整上、下时窗大小，直到所有已知钻井与地震子波衰减谱匹配，最终确定合适的上、下时窗大小。
18.根据权利要求14所述的基于地震子波衰减谱的油气检测装置，其特征在于，所述切片分析单元包括切片模块，用于在地震子波衰减谱对应的三维数据体上，按垂直于频率轴进行切片，得到不同频率的地震子波衰减谱切片平面图；预测模块，用于选择高频的地震子波衰减谱切片平面图来进行含油气性的平面预测， 确定衰减大于预设的衰减能量值的区域为可能的含油气区域。
全文摘要
本发明实施例提供了一种基于地震子波衰减谱的油气检测方法及装置，该方法包括采用独立成分分析方法对地震数据进行去噪处理；对去噪处理后的地震数据进行目的层顶、底地震层位追踪；沿目的层顶、底开上、下两个时窗；对上、下两个时窗内的地震数据进行频谱分析，得到其频谱分布；对频谱分布进行复倒谱分析，得到地震数据振幅谱的复倒谱；根据预设的低通时窗从地震数据振幅谱的复倒谱中分离出地震子波振幅谱的复倒谱；对地震子波振幅谱的复倒谱进行逆变换，得到上、下两个时窗内的地震子波谱；对地震子波谱进行差值计算，得到地震子波衰减谱；通过分析地震子波衰减谱来检测地层的含油气性。本发明能实现准确、敏感地识别油气。
文档编号G01V1/30GK102305943SQ20111020684
公开日2012年1月4日 申请日期2011年7月22日 优先权日2011年7月22日
发明者刘建新, 刘晓梅, 姚清洲, 孙东, 孙勤华, 张虎权, 潘建国, 王宏斌, 王锦喜 申请人:中国石油天然气股份有限公司