一种密度曲线校正方法及系统的制作方法

一种密度曲线校正方法及系统的制作方法
【专利摘要】本申请提供一种密度曲线校正方法及系统。所述方法包括：采集原始测井数据；根据所述原始测井数据，获取测井曲线，所述测井曲线包括伽玛曲线、声波曲线、电阻率曲线以及密度曲线；根据所述测井曲线，并结合实际钻井岩性信息确定密度曲线异常值段；利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到拟合表达式；根据所述拟合表达式，对异常值段的密度曲线进行校正。本申请实施例提供的一种密度曲线校正方法及系统，利用多曲线拟合的方法，可以确定其他测井曲线和密度曲线之间的关系，并以数学形式表示出来。那么对于存在异常值的密度曲线段，就可以利用其他测井曲线与密度曲线之间的关系，快速准确地对异常值进行校正。
【专利说明】-种密度曲线校正方法及系统

【技术领域】
[0001] 本申请设及地球物理勘探【技术领域】，特别设及地球物理勘探中的测井技术，具体 讲的是一种密度曲线校正方法及系统。

【背景技术】
[0002] 测井技术是地球物理勘探技术的一个分支，可W应用物理学的原理解决地质问 题，并已在石油、天然气、金属矿、煤田、工程机水文地质等许多方面得到应用。测井技术可 W在已经钻探好的井眼中，将各种激发设备和接收设备放置在一根钢管内。其中，激发设备 和接收设备可W保持几米的距离，并通过电缆连接到地面上。在地面上，通过匀速下降电缆 和匀速提升电缆，可W灵活地变动激发设备和接收设备在井中的位置。在井中，激发设备可 W进行信息激发，该些激发的信息穿过井壁地层，由接收设备接收并通过电缆传输到地面 存储设备中。
[000引在测井技术中，密度测井可W进行井周围地层密度的测量。密度测井可W在已经 钻探好的井中布置测井设备。测井设备包括位于屏蔽体内的放射性源和放射性探测器。由 放射性源发出的伽玛射线射向井壁地层，该些伽玛射线可W看作为高速粒子，在地层中与 电子碰撞。每次碰撞时伽玛射线传递能量给电子而失去一部分能量，能量减少后的伽玛射 线继续前进并可W由放射性探测器接收。伽玛射线能量的变化可W作为地层密度的指示被 记录下来，该样就形成了密度曲线。密度曲线是制作地震合成记录必不可少的曲线，是地球 物理地震地质研究中很重要的参数。但是密度测井对井眼质量要求高，对于扩径、不平整井 壁均需要进行校正，否则其值往往不能精确的反映地层的密度特征，难W实现对储层进行 精确预测。测井密度曲线的校正就是利用各种办法对存在异常值的密度曲线进行校正，达 到曲线能够精确的反映井眼周围地层的地球物理特征的目的。
[0004] 目前对密度曲线的校正，通常包括W下几种方法；一种是Garner方程校正方法。 该方法利用了的测井数据进行统计分析，得到纵波测井曲线与密度测井曲线之间的对应关 系。然而该种方法不具备普遍性，基于大量墨西哥湾钻井数据得到的该种对应关系并不一 定适合其他地区的地质情况，因此难W实现对测井密度曲线的精确校正。还有一种是岩石 物理模拟正演方法。该种方法需要把测井地层评价处理得到的各矿物含量、孔隙度、饱和度 等参数收集全，然后选用岩石物理模型进行纵、横波速度的建模，再通过微调骨架点参数， 使模型数据和实测数据达到很高的相关性。该种方法处理过程中需要反复试验，才能达到 最佳效果，需要大量的人力和设备资源。并且该种方法需要大量的岩屯、参数，对岩屯、参数不 全的井，密度曲线校正仍然不能进行。


【发明内容】

[0005] 为解决上述问题，本申请实施例提供一种密度曲线校正方法及系统，W达到能够 准确进行密度曲线校正的目的。
[0006] 本申请实施例提供的一种密度曲线校正方法及系统是该样实现的：
[0007] 一种密度曲线校正方法，包括：
[000引采集原始测井数据；
[0009] 根据所述原始测井数据，获取测井曲线，所述测井曲线包括伽玛曲线、声波曲线、 电阻率曲线W及密度曲线；
[0010] 根据所述测井曲线，并结合实际钻井岩性信息确定密度曲线异常值段；
[0011] 利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到拟合表达式；
[0012] 根据所述拟合表达式，对异常值段的密度曲线进行校正。
[0013] 一种密度曲线校正系统，包括：原始测井数据采集单元，测井曲线获取单元，异常 值段确定单元，拟合单元，校正单元，其中：
[0014] 所述原始测井数据采集单元，用来采集原始测井数据；
[0015] 所述测井曲线获取单元，用来根据所述原始测井数据，获取测井曲线，所述测井曲 线包括伽玛曲线、声波曲线、电阻率曲线W及密度曲线；
[0016] 所述异常值段确定单元，用来根据所述测井曲线，并结合实际钻井岩性信息确定 密度曲线异常值段；
[0017] 所述拟合单元，用来利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合， 得到拟合表达式；
[0018] 所述校正单元，用来根据所述拟合表达式，对异常值段的密度曲线进行校正。
[0019] 本申请实施例提供的一种密度曲线校正方法及系统，利用多曲线拟合的方法，在 井眼环境好的井段，可W对其他测井曲线进行估计和预测，确定该些测井曲线和密度曲线 之间的关系，并W数学形式表示出来。那么对于存在异常值的密度曲线段，就可W利用其他 测井曲线与密度曲线之间的关系，快速准确的对异常值进行校正。

【专利附图】

【附图说明】
[0020] 图1为本申请实施例提供的一种密度曲线校正方法的流程图；
[0021] 图2为本申请实施例中伽玛-密度的交汇示意图；
[0022] 图3为本申请实施例中经过校正后的伽玛-密度交汇示意图；
[0023] 图4为本申请另一实施例提供的一种密度曲线校正方法的流程图；
[0024] 图5为本申请实施例提供的一种密度曲线校正系统的功能模块图；
[0025] 图6为本申请另一实施例提供的一种密度曲线校正系统的功能模块图；
[0026] 图7为本申请实施例提供的一种密度曲线校正系统中测井曲线获取单元的功能 模块图。

【具体实施方式】
[0027] 本申请实施例提供一种密度曲线校正方法及系统。
[002引为了使本【技术领域】的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都应当属于本申请保护 的范围。
[0029] 图1为本申请实施例提供的一种密度曲线校正方法的流程图。如图1所示，所述 方法包括W下步骤：
[0030] S100 ;采集原始测井数据。
[0031] 原始测井数据可W用过野外测井施工来进行采集。在已经钻探好的井眼中，可W 将各种激发设备和接收设备放置在一根钢管内。其中，激发设备和接收设备可W保持几米 的距离，并通过电缆连接到地面上。在地面上，通过匀速下降电缆和匀速提升电缆，可W灵 活地变动激发设备和接收设备在井中的位置。在井中，激发设备可W进行信息激发，该些激 发的信息穿过井壁地层，由接收设备接收并通过电缆传输到地面存储设备中。接收设备接 收到的数据就可W作为原始测井数据被记录在数字磁带上。
[0032] S200;根据所述原始测井数据，获取测井曲线，所述测井曲线包括伽玛曲线、声波 曲线、电阻率曲线W及密度曲线。
[0033] 原始测井数据记录在数字磁带上后，可W进一步地对记录在磁带上的原始测井数 据进行处理。该些原始测井数据包括了井壁周围地层的电阻率、声波值、伽玛值W及密度等 物理信息。测井数据的处理可W将该些记录在磁带上的物理信息转换为不同的测井曲线， 具体是通过W下两个步骤来实现的：
[0034] S201 ;对测井数据进行检查与预处理。
[0035] 首先可W将记录在数字磁带上的测井数据打印出来，从而进一步地检查测井数据 文件的鉴别号、深度、采样间距、采样数据是否合理准确，并改变记录内容的记录格式，使之 成为便于计算机使用的文件格式。接着可W对数据进行深度校正、平滑滤波处理、环境校正 等预处理，从而得到采样间距一致、深度对齐、数据正确的测井数据。
[0036] S202 ;根据处理后的测井数据，获取测井曲线。
[0037] 针对步骤S201预处理后的测井资料，可W运用各种测井分析软件对数据进行自 动处理解释，从而获得钻井中地层的各种物理信息，并最终W连续曲线图的方式显示出来。 例如，可W利用Geo化ame或者express等测井软件，对测井数据进行处理解释，从而可W将 钻井中地层的各种物理信息W连续曲线图的方式显示出来。本申请实施例通过对测井数据 的处理，可W得到钻井中地层的伽玛曲线、声波曲线、电阻率曲线W及密度曲线。
[003引 S300;根据所述测井曲线，并结合实际钻井岩性信息确定密度曲线异常值段。
[0039] 在本申请实施例中通过对测井数据的处理，可W得到钻井中地层的伽玛曲线、声 波曲线、电阻率曲线W及密度曲线。地层密度的大小与地层的岩性、电性、孔隙、流体等因素 均有着一定的关系。其中，地层的岩性特征可W由伽玛曲线中的伽玛值进行表征；地层的电 性特征则可W很好地通过电阻率曲线进行解释；地层的孔隙和流体特征可W通过地层的速 度来进行分辨，而地层的速度可W直接由地层的声波曲线进行求解。因此地层的密度曲线 与该地层的伽玛曲线、声波曲线W及电阻率曲线也存在着某种确定性或者统计性的关系。 另外，伽玛曲线、声波曲线W及电阻率曲线由于其测试的深度较广，从而不易失真。那么便 可W通过横向比较伽玛曲线、声波曲线、电阻率曲线与密度曲线的对应关系，从而对密度曲 线进行分析，确定出密度曲线的异常值段。
[0040] 现在W伽玛曲线为例来解释如何通过横向对比伽玛曲线和密度曲线，来确定出密 度曲线的异常值段。通常而言，泥岩的密度要大于砂岩的密度，泥岩的放射性也要大于砂 岩的放射性。泥岩和砂岩的密度特征可W通过两者的密度曲线进行表征，而泥岩和砂岩的 放射性特征也可w通过伽玛曲线进行表征。那么泥岩的测井曲线特征就应该为高密度高伽 玛，砂岩的测井曲线特征就应该为低密度低伽玛。通过对两者的伽玛曲线和密度曲线进行 横向分析，如果出现高伽玛低密度或低伽玛高密度的区域，则可W视为异常值段。该异常值 段的密度曲线就需要进行进一步的校正。图2为本申请实施例中伽玛-密度的交汇示意图。 如图2所示，横坐标代表地层的伽玛值，纵坐标代表该地层的密度，伽玛-密度交汇示意图 中的点就代表地层的采样点对应的伽玛值和密度。从图2中可W看出，绝大多数地层采样 点对应的伽玛值和密度均符合高密度高伽玛W及低密度低伽玛的特性。但是在图2中还有 部分离散的地层采样点具有高伽玛低密度或低伽玛高密度的特性。例如位于黑色曲线框内 的地层采样点则明显具有高密度低伽玛的特性。该些具有异常伽玛-密度特性的地层采样 点就可W被认为是异常值段。该异常值段的密度曲线就需要进行进一步的校正。
[0041] S400;利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到拟合表达 式。
[0042] 在确定出密度曲线的异常值段后，可W对正常值的数据进行处理。在正常值的数 据中，可W利用地层的伽玛曲线、声波曲线W及电阻率曲线对密度曲线进行拟合，得到拟合 表达式。具体拟合的过程如下所述：
[0043] 由于地层的密度曲线与该地层的伽玛曲线、声波曲线W及电阻率曲线存在着某种 确定性或者统计性的关系，那么就可W利用多曲线拟合，通过伽玛曲线、声波曲线W及电阻 率曲线去估算密度曲线。首先可W对正常的伽玛曲线、声波曲线W及电阻率曲线进行计算， 得到伽玛曲线、声波曲线W及电阻率曲线和密度测井曲线值之间的关系，并把该种关系用 一个数学公式来表示。该数学关系式可W表示为：
[0044] den = aX+bY+cZ+d
[0045] 其中，den代表地层的密度值，a, b，c，d为待求取的系数，X代表地层的伽玛值，Y 代表地层的声波值，Z代表地层的电阻率值。
[0046] 从上述公式可W看出，地层的密度值可W通过地层的伽玛值、声波值W及电阻率 值来表示。现在可W在地层压力、沉积环境W及岩性特征相同的区域，选取4组正常的地层 密度值、伽玛值、声波值W及电阻率值带入上述公式，得到四元一次方程组。
[0047] 通过对该四元一次方程组的求解，便可W得到a, b，C，d该四个待求取的系数。将 该四个已求取的系数带入上述公式后，便可W得到利用伽玛曲线、声波曲线W及电阻率曲 线针对密度曲线的拟合表达式。
[0048] S500 ;根据所述拟合表达式，对异常值段的密度曲线进行校正。
[0049] 步骤S400中通过正常值的数据拟合出的密度曲线拟合表达式能够反映出密度曲 线与伽玛曲线、声波曲线W及电阻率曲线之间的联系。那么对于密度曲线的异常值段，便可 W利用伽玛值、声波值W及电阻率值通过该拟合表达式进行计算，得到理论的密度值，并用 该理论的密度值替换异常值段的密度值，达到对异常值段的密度值进行校正的目的。图3 为本申请实施例中经过校正后的伽玛-密度交汇示意图。通过对比图2和图3可W看出， 原本位于异常值段的地层采样点通过拟合表达式的计算，已经全部得到校正。从图3中可 W看出，校正后的伽玛-密度交汇示意图中的地层采样点均符合高密度高伽玛W及低密度 低伽玛的特性。
[0化0] 在本申请的另一实施例中，在对异常值段的密度曲线进行校正后，还对该校正的 准确性进行了确定，并对校正的结果进行进一步的修正。图4为本申请另一实施例提供的 一种密度曲线校正方法的流程图。如图4所示，所述方法还包括：
[0化1] S600 ;对所述校正后的密度曲线进行井震标定。
[0052] 在对密度曲线进行校正后，可W利用校正后的密度曲线，在井震标定软件里面，首 先把声波曲线和校正后的密度曲线进行乘积运算，得到波阻抗曲线。接着可W由波阻抗曲 线直接计算得到反射系数序列，计算的公式如下所示：
[005引 Ri 二 7山I ~ 7' ' i = l，2,3wN 之i+l十
[0化4] 其中，而表示反射系数序列，Z i为第i层的波阻抗，N为地下反射界面的数目。
[0055] 接着可W对测井数据通过反權积，求取子波。利用子波与上述的反射系数序列进 行權积得到人工合成地震记录。由于对密度曲线进行校正后，测井数据也将得到校正，该样 最终得到的人工合成地震记录将会受到密度曲线校正的影响。将该人工合成地震记录与实 际地震记录进行对比，可W根据对比结果再次校正密度曲线，使得人工合成地震记录与实 际地震记录达到最佳匹配，从而能够对校正的结果进行进一步修正，完成井震标定。
[0056] 本申请实施例还提供一种密度曲线校正系统。图5为本申请实施例提供的一种密 度曲线校正系统的功能模块图。如图5所示，所述系统包括：
[0化7] 原始测井数据采集单元1，用来采集原始测井数据；
[005引测井曲线获取单元2,用来根据所述原始测井数据，获取测井曲线，所述测井曲线 包括伽玛曲线、声波曲线、电阻率曲线W及密度曲线；
[0059] 异常值段确定单元3,根据所述测井曲线，确定密度曲线异常值段；
[0060] 拟合单元4,利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到拟合 表达式；
[0061] 校正单元5,根据所述拟合表达式，对异常值段的密度曲线进行校正。
[0062] 此外，本申请一优选实施例还可W对图5所示的系统做进一步的补充。图6为本 申请另一实施例提供的一种密度曲线校正系统的功能模块图。如图6所示，所述系统除了 包含图5所示系统的功能模块外，还包括：
[0063] 井震标定单元6,用来对所述校正后的密度曲线进行井震标定。
[0064] 图7为本申请实施例提供的一种密度曲线校正系统中测井曲线获取单元的功能 模块图。如图7所示，测井曲线获取单元2具体包括：
[00化]预处理模块201，用来对测井数据进行检查与预处理；
[0066] 测井曲线获取模块202,用来根据处理后的测井数据，获取测井曲线。
[0067] 通过上述本申请提供的实施例可W看出，本申请利用多曲线拟合的方法，在井眼 环境好的井段，可W对其他测井曲线进行估计和预测，确定该些测井曲线和密度曲线之间 的关系，并W数学形式表示出来。那么对于存在异常值的密度曲线段，就可W利用其他测井 曲线与密度曲线之间的关系，快速准确的对异常值进行校正。
[0068] 在20世纪90年代，对于一个技术的改进可W很明显地区分是硬件上的改进（例 如，对二极管、晶体管、开关等电路结构的改进）还是软件上的改进（对于方法流程的改 进）。然而，随着技术的发展，当今的很多方法流程的改进已经可W视为硬件电路结构的直 接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路 结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑 器件（Programm油le Logic Device, PLD)(例如现场可编程口阵列（Field Programm油le Gate Array, FPGA))就是该样一种集成电路，其逻辑功能由用户对器件编程来确定。由设 计人员自行编程来把一个数字系统"集成"在一片PLD上，而不需要请巧片制造厂商来设 计和制作专用的集成电路巧片2。而且，如今，取代手工地制作集成电路巧片，该种编程也 多半改用"逻辑编译器（logic compiler)"软件来实现，它与程序开发撰写时所用的软件 编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描 述语言（Hardware Description Language,皿L)，而皿L也并非仅有一种，而是有许多种， 女日 ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、 JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是 VHDL(Ve;ry-Hi曲-Speed Integrated Circuit Hardware Description Language)与 Verilog2。本领域技术人员 也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路 中，就可W很容易得到实现该逻辑方法流程的硬件电路。
[0069] 控制器可W按任何适当的方式实现，例如，控制器可W采取例如微处理器或处理 器W及存储可由该（微）处理器执行的计算机可读程序代码（例如软件或固件）的计算 机可读介质、逻辑口、开关、专用集成电路（Application Specific Integrated Cir州it, ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于W下微控制 器；ARC 625D、Atmel AT91SAM、Microchip PIC18巧服20 W及 Silicone L油S C8051F320， 存储器控制器还可W被实现为存储器的控制逻辑的一部分。
[0070] 本领域技术人员也知道，除了 W纯计算机可读程序代码方式实现控制器W外，完 全可W通过将方法步骤进行逻辑编程来使得控制器W逻辑口、开关、专用集成电路、可编程 逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此该种控制器可W被认为是一种 硬件部件，而对其内包括的用于实现各种功能的装置也可W视为硬件部件内的结构。或者 甚至，可W将用于实现各种功能的装置视为既可W是实现方法的软件模块又可W是硬件部 件内的结构。
[0071] 上述实施例阐明的系统、装置、模块或单元，具体可W由计算机巧片或实体实现， 或者由具有某种功能的产品来实现。
[0072] 为了描述的方便，描述W上装置时W功能分为各种单元分别描述。当然，在实施本 申请时可W把各单元的功能在同一个或多个软件和/或硬件中实现。
[0073] 通过W上的实施方式的描述可知，本领域的技术人员可W清楚地了解到本申请可 借助软件加必需的通用硬件平台的方式来实现。基于该样的理解，本申请的技术方案本质 上或者说对现有技术做出贡献的部分可软件产品的形式体现出来，该计算机软件产品 可W存储在存储介质中，如R0M/RAM、磁碟、光盘等，包括若干指令用W使得一台计算机设备 (可W是个人计算机，服务器，或者网络设备等）执行本申请各个实施例或者实施例的某些 部分所述的方法。
[0074] 本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部 分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实 施例而言，由于其基本相似于方法实施例，所w描述的比较简单，相关之处参见方法实施例 的部分说明即可。
[0075] 本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、月良 务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置 顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括W上任何系统或设备 的分布式计算环境等等。
[0076] 本申请可W在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序 模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组 件、数据结构等等。也可W在分布式计算环境中实践本申请，在该些分布式计算环境中，由 通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可W 位于包括存储设备在内的本地和远程计算机存储介质中。
[0077] 虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和 变化而不脱离本申请的精神，希望所附的权利要求包括该些变形和变化而不脱离本申请的 精神。
【权利要求】
1. 一种密度曲线校正方法，其特征在于，包括： 采集原始测井数据； 根据所述原始测井数据，获取测井曲线，所述测井曲线包括伽玛曲线、声波曲线、电阻 率曲线以及密度曲线； 根据所述测井曲线，并结合实际钻井岩性信息确定密度曲线异常值段； 利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到拟合表达式； 根据所述拟合表达式，对异常值段的密度曲线进行校正。
2. 如权利要求1所述的一种密度曲线校正方法，其特征在于，所述方法还包括： 对所述校正后的密度曲线进行井震标定。
3. 如权利要求1或2所述的一种密度曲线校正方法，其特征在于，所述根据所述原始测 井数据，获取测井曲线具体包括： 对测井数据进行检查与预处理； 根据处理后的测井数据，获取测井曲线。
4. 一种密度曲线校正系统，其特征在于，包括：原始测井数据采集单元，测井曲线获取 单元，异常值段确定单元，拟合单元，校正单元，其中： 所述原始测井数据采集单元，用来采集原始测井数据； 所述测井曲线获取单元，用来根据所述原始测井数据，获取测井曲线，所述测井曲线包 括伽玛曲线、声波曲线、电阻率曲线以及密度曲线； 所述异常值段确定单元，用来根据所述测井曲线，并结合实际钻井岩性信息确定密度 曲线异常值段； 所述拟合单元，用来利用伽玛曲线、声波曲线、电阻率曲线对密度曲线进行拟合，得到 拟合表达式； 所述校正单元，用来根据所述拟合表达式，对异常值段的密度曲线进行校正。
5. 如权利要求4所述的一种密度曲线校正系统，其特征在于，所述系统还包括： 井震标定单元，用来对所述校正后的密度曲线进行井震标定。
6. 如权利要求4或5所述的一种密度曲线校正系统，其特征在于，所述测井曲线获取单 元具体包括： 预处理模块，用来对测井数据进行检查与预处理； 测井曲线获取模块，用来根据处理后的测井数据，获取测井曲线。
【文档编号】G01V11/00GK104502996SQ201410814789
【公开日】2015年4月8日 申请日期:2014年12月23日 优先权日:2014年12月23日
【发明者】邹义, 杨洋, 张泉, 徐博, 马培领, 冯许魁, 刘永雷 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司