二维地震过障碍观测系统模式及其参数设计.pdf

3 四川省双流市中和镇石油物探研究中心 , 610212 本文于 2005 年 3 月 25 日收到 , 修改稿于同年 9 月 9 日收到。 采集技术 二维地震过障碍观测系统模式及其参数设计 梁顺军 3 ①② 肖 敏 ① 唐 怡 ① 蔡友宏 ① 张延充 ① 刘定锦 ① ①四川石油地调处物探研究中心 ; ②新疆塔里木油田勘探开发研究院 摘 要 梁顺军 , 肖敏 , 唐怡 , 蔡友宏 , 张延充 , 刘定锦 . 二维地震过障碍观测系统模式及其参数设计 . 石油地球物理勘探 , 2006, 41 1 1～ 7 在二维地震勘探中 , 为了避免因地表障碍物使地震反射剖面出现间断现象 , 需要改变观测系统设计 , 跨越 障碍物 江河、城镇等 , 以保证反射同相轴能连续追踪对比。在以往地震勘探中 , 有时因过障碍观测系统设计不 合理 , 不但增加勘探费用 , 而且影响采集质量。因此只有科学合理地设计过障碍观测系统模式和参数 , 才能保证 地震测线顺利通过障碍物 , 得到较好的障碍物下方的地震资料并降低勘探成本。为此本文在调查和研究大量野 外实际的过障碍观测系统的基础上 , 总结归纳成三大类过障碍观测系统模式 , 并分析其特点和应用范围。采用 图解法 , 推导出主要模式的跨越宽度、最小炮检距及最大炮检距等参数公式。采用这套过障碍观测模式不仅可 以得到和障碍区两侧等质量的地震资料 , 而且可以有效地降低勘探成本。 关键词 山地 二维地震 测线 障碍物 观测系统 模式 参数 1 问题的提出 地震资料采集中的地表障碍物泛指无法按常规 方法布置炮点或无法埋置检波器的地段 , 如江河湖 泊、厂矿城镇等。在这种情况下需要采用过障碍的观 测系统。以往地震勘探野外过障碍观测系统的形式 可概括为“炮点偏离”、“折线避开”、“加密井位”、“加 大最小炮检距”及“排列移动”等。这些形式虽然在生 产中发挥了积极作用 \[ 1, 2 \], 但也存在不足 其一是为 了保证障碍段的 CM P 点达到满覆盖次数 , 障碍段两 侧加炮太多 , 使部分 CM P 点叠加次数是设计的几 倍 , 造成浪费 ; 其二是过障碍最小炮检距太大 , 在地 震剖面上形成很深的“ W ”形缺口 , 不仅使障碍段浅 层资料缺失 , 而且造成过障碍段的中深层地震反射 特征和正常段存在明显差异 , 影响构造解释和储层 图 1 过障碍时采用不合理观测系统得到的地震剖面 2006 年 2 月 石油地球物理勘探 第 41 卷 第 1 期 研究 图 1。为此 , 笔者认为应当改进过障碍观测系 统的设计 , 使其更趋合理和实用 , 因此有必要从理论 上对过障碍观测系统进行分析和归类 , 以便建立一 种基本模式 , 方便指导野外生产。 本文以采用 120 道中间放炮的 6 次覆盖观测系 统为例建立模式 , 但也适用于其他仪器道数和发炮 方式。采用图解法进行数学公式推导 , 确定过障碍观 测系统参数。通常正常中间 或双边 放炮二维观测 系统有两种表示方法 图 2 , 但意义是一样的 其一 是分离法 , 左、右支排列分别表示在列线图下、上区 域 , 视图直观 , 但图幅面积大 , 野外施工携带不方便 ; 其二是合并法 , 左、右支排列均表示在列线图上区 域 , 图幅面积小 , 携带方便。为了行文和公式推导方 便 , 本文采用分离法。 图 2 二维观测系统中间和双边放炮两种图示法 a分离列线图 ; b合并列线图 2 二维地震过障碍观测系统图示 为了便于理解过障碍观测系统 , 采用图示法来 表述其特点、术语和定义 图 3。文中的变量均以道 数表示。 1 正常排列 L 及过障碍排列 L Z 正常排列 L 是指在无障碍情况下按正常观测系统设计炮点与接 收点的对应关系。过障碍排列 L Z 是指跨越障碍的排 列 , 其炮点和接收点位于障碍区两侧 , 相对于正常观 测系统而言 , 炮点与接收点的对应关系已发生变化。 其实该观测系统相当于用延长时距曲线 加大炮检 距 方法来获得障碍物下面的深层反射资料。显然 , 因受初至波和浅层折射波的严重干扰 , 此法对于浅 层反射资料叠加成像困难 , 所以在资料处理中应对 浅层资料进行切除 图 1。 2障碍宽度 Z 及过障碍观测系统跨越宽度 Z k 障碍宽度 Z 是指野外实测的障碍物宽度。 Z k 是指过 障碍观测系统所能跨越障碍物的最大宽度 , 简称跨越 宽度 , 它反映了过障碍观测系统跨越障碍的能力。只 有当 Z≤ Z k 时 , 观测系统才能顺利地通过障碍物。 3正常最小炮检距 L及过障碍最小炮检距 Lz 在正常排列情况下 , L为最小炮检距 , 又称偏移距 ; Lz 为在一个过障碍排列上炮点与排列第一道之间 距离。通常 , L≤ Lz。 4正常最大炮检距 X 及过障碍最大炮检距 X z X 为在正常排列情况下 , 炮点到最远一道的距 离。在过障碍排列情况下 , X z 为炮点到最远一道的 距离。显然 , X ≤ X z。 5 正常炮点距 D 及加密炮点距 d D 为按正 常排列设计的炮点距。 d 为按过障碍排列设计的炮 点距 , 一般障碍两侧加密井位 炮点 , 可以等间隔或 不等间隔 , 通常 d 小于正常炮点距的整数倍 , 其最 图 3 二维观测系统过障碍模式图示法 2 石 油 地 球 物 理 勘 探 2006 年 小值为一个道间隔。一般 , d≤ D。 6不满排列三角区 在邻近障碍的排列区段 , 每一炮接收道数小于正常设计的接收道数 , 在观测 系统图上表现形式如三角形。其宽度一般为 [ N 1 - 1 L]x 1 7 其他参数 x 为道间距 ; N 、 n 分别为仪器 道数和覆盖次数 ; N 1、 N 2 分别为右支和左支排列的 接收道数 , 对于中间发炮情况有 N 1 N 2 N 2。在 野外生产中炮检距大小常用道数和距离来表示 , 其 关系为距离 \[道数 - 1 L]道间距。 3 过障碍观测模式特点 依据障碍物能否放炮、能否埋置检波器及障碍物 宽度与排列的关系 , 过障碍观测系统可分为三大类。 3. 1 不能放炮、能埋检波器 3. 1. 1 中间放炮对称排列观测模式 图 4 为中间放炮对称排列模式。其特点是炮检 距不变 , 在障碍两侧加密的井位数要满足障碍段的 覆盖次数要求 , 此种观测系统模式适用于跨越较小 的障碍物和高分辨率地震勘探。将该模式的左支或 右支排列去掉后 , 就变成单边放炮过障碍模式。随着 地震仪器道数的成倍增加 , 单边和双边放炮观测系 统的炮检距变得太大 , 因此现代地震资料采集中一 般不采用这种方式。 图 4 中间放炮对称排列观测模式示意图 3. 1. 2 中间放炮不对称排列观测模式 图 5 为中间放炮不对称排列模式。当障碍宽度 较大时 , 其丢井个数 W n 满足下式 W n ≥ n - 1 L D 2 若沿用图 4 中的观测模式 , 地震剖面就要出现间断 或降低覆盖次数 , 所以必须采用图 5 的观测模式 , 即 以加大右支或左支排列接收道数的办法来满足障碍 段的覆盖次数。该模式与图 4 的放炮方式相同 , 不同 的是右支和左支共炮点线的接收道数不等 , 即所谓 排列转换 , 例如 120 道对称排列中间放炮 , 可转换成 50～ 70 道、 40～ 80 道、 30～ 90 道等不对称排列形 式。通过排列转换 , 加大了炮检距 , 使得跨越障碍能 力增强。 图 5 中间放炮不对称排列观测模式示意图 3. 2 不能放炮、不能埋检波器 其特点是在障碍区内不能埋检波器也不能打井 放炮 , 如江河湖泊和厂矿城镇 , 即只能丢炮和丢道 , 为此可考虑采用如下三种办法。 3. 2. 1 单边移动排列观测模式 图 6 为单边移动排列模式。单边移动是指障碍 一侧每放一炮后炮点向前或向后移动若干个加密炮 点距 d , 而障碍的另一侧排列不动 , 因此每炮过障碍 排列的最小和最大炮检距是不等的。 图 6 单边移动排列观测模式示意图 3. 2. 2 双边移动排列观测模式 图 7 为双边移动排列观测模式。双边移动是指 障碍一侧每放一炮后 , 炮点和障碍另一侧排列向同 一方向移动若干个相同加密炮点距 d , 故每炮过障 碍排列的最小及最大炮检距是相等的。 图 7 双边移动排列观测模式示意图 3. 2. 3 “ X”形排列观测模式 图 8 为“ X”形排列观测模式 , 障碍段内丢道和丢 炮 , 过障碍排列的接收道数少于正常排列接收道数 , 但过障碍最小和最大炮检距分别等于正常最小和最 3 第 41 卷 第 1 期 梁顺军等 二维地震过障碍观测系统模式及其参数设计 大炮检距。该模式适用于高分辨率地震勘探和跨越 小型障碍物。 图 8 “ X”形排列观测模式示意图 需要指出的是 , 把上述某一种基本模式重复使 用或将几种模式组合应用 , 以获得连续多个障碍物 地下的反射资料 , 故将这类观测模式称为复合观测模 式。限于篇幅 , 不再赘述。 3. 3 二维、三维组合观测模式 当障碍宽度大于正常最大炮检距或一个排列长 度时 , 称为大型障碍 , 这时地震剖面将出现间断 ; 或 仅使用二维观测系统难以得到障碍物下面的地震资 料时 , 可采用二维、三维观测系统组合模式来解决。 如图 9 所示 , 其中三维观测系统采用多线条带法 , 由 相隔一定距离的多组接收条带与多组炮线正交而构 成。炮点和接收线对称分布于“剖面输出线” 沿障 碍物“长轴”方向 两侧。障碍物两端仍按正常二维 观测系统施工。 图 9 二维、三维组合观测模式示意图 4 过障碍观测模式的参数设计 过障碍观测模式的主要参数有跨越宽度、最小 炮检距和最大炮检距。在下列公式推导及图示中涉 及 B 和 h 表示的有关变量均是中间变量 , 不是观测 系统参数。 4. 1 不能放炮、能埋检波器 4. 1. 1 中间放炮对称排列观测模式 1跨越宽度 Z k, 如图 10 所示 , 其表达式为 B \[2 N 1 - 1 L ]x 2 \[N - 2 2L]x 2 3 又有 B Z k n - 1 d x 4 化简整理式 3和式 4得跨越宽度 Z k \[N - 2 2L - 2 n - 1 d ]x 2 5 由式 5 可知 , 加密炮点距 d 越小 , 跨越宽度 Z k 越 大 , 反之亦然。通常 d 为 x 的整数倍 , 小于正常炮点 距 D。当障碍宽度 Z 较小时 , d 可取大些 , 使炮点均匀 分布。另外 , 覆盖次数 n 较小时 , d 也可取大些 , 反之 亦然。 2 过障碍最小炮检距 , 其表达式为 Lz L 6 3过障碍最大炮检距 , 其表达式为 X z X 7 4 石 油 地 球 物 理 勘 探 2006 年 图 10 中间放炮对称排列观测模式参数计算示意图 4. 1. 2 单边放炮观测模式 图 11 为单边放炮观测模式参数计算图示 , 其表 达式为 B \[N - 1 L]x 2 8 又有 B Z k n - 1 d x 9 化简整理式 8和式 9 得跨越宽度 Z k \[N - 1 L - 2 n - 1 d ]x 2 10 其余参数的确定与对称排列中间放炮模式相同。 图 11 单边放炮观测模式参数计算示意图 4. 1. 3 正常炮点距两端单边放炮观测模式 一般单边放炮观测系统是在测线的一端 小 或 大 桩号放炮 , 大 或小 桩号接收 向一个方向放炮 , 故可称为一端单边放炮观测系统 \[ 2 \]; 而两端单边放 炮观测是在测线的两端 , 一端小桩号放炮大桩号接 收 , 而另一端大桩号放炮小桩号接收 , 向相对方向 放炮。 对于跨越打井困难段 , 一般过障碍观测系统采 取的方法是加密井位或加大最小炮检距 , 这样不利 于复杂构造和高分辨率地震勘探 , 影响浅层资料的 效果。而采用正常炮点距两端单边发炮观测模式可 弥补上述过障碍模式的缺陷。以跨越高陡构造顶部 灰岩出露区为例 图 12 , 在野外施工中 , 排列布设在 反射界面上倾方向 , 下倾方向激发 , 上倾方向接收。 由两翼的单边发炮观测系统向构造顶部逐步会合 , 共炮点线依次交会“搭篷”。在不改变正常观测系统 参数的条件下 , 构造顶部灰岩出露段可不打井 , 也可 获得满覆盖次数。其跨越宽度 Z k 表达式为 Z k \[N - 1 L - n - 1D ]x 11 把式 11的计算结果列成表 1, 与常规单边放炮 观测系统比较 , 可以看出该过障碍观测系统有以下 特点。 1降低勘探成本 由表 1 可知 x 越大 , 跨越 宽度 Z k 越大 , 即跨越障碍能力越强 ; 覆盖次数越高 , 少打炮井数越多 , 例如 N 120, L 4. 5, n 20, x 40m , 跨越障碍宽度为 2660m , 障碍范围内少打炮 井 22 口。 图 12 正常炮点距两端单边放炮模式参数计算示意图 5 第 41 卷 第 1 期 梁顺军等 二维地震过障碍观测系统模式及其参数设计 表 1 正常炮点距两端单边放炮观测模式 在障碍范围内少打炮井数一览表 覆盖次数 道间距 m 跨越宽度道数 距离 m 障碍范围内少打炮井数 施工参数 10 10 20 30 40 69. 5 695 1390 2085 2780 12 单边发炮 N 120 L 4. 5 D 6 20 10 20 30 40 66. 5 665 1330 1995 2660 22 单边发炮 N l20 L 4. 5 D 3 2按正常观测系统设计施工 没有加密井位 , 障碍段内仍获得满覆盖次数 , 有利于得到较好的构 造顶部浅层资料。 3 有利于选择激发条件 在陡构造两翼有利 于选择砂泥岩条件进行激发 , 易得到较好的单炮 记录。 4. 2 不能放炮、不能埋检波器 4. 2. 1 单边移动观测模式 图 13 为单边移动观测模式参数计算图示 , 其主 要参数表达式如下。 图 13 单边移动观测模式参数计算示意图 1跨越宽度 Z k \[N - 1 - 2 n - 1 d - L]x 12 2过障碍最小炮检距 Lz Z 13 3过障碍最大炮检距 X z \[Z N - 1 n - 1 d ]x 14 4. 2. 2 双边移动观测模式 图 14 为双边移动观测模式参数计算图示 , 其主 要参数表达式如下 1跨越宽度 Z k \[N - 1 - n - 1 d - L]x 15 2过障碍最小炮检距 Lz Z n - 1 d x 16 图 14 双边移动模式参数计算示意图 3 过障碍最大炮检距 X z Z \[N - 1 n - 1 d ]x 17 比较式 12 和式 15 可知 , 双边移动跨越宽度大于 单边移动的跨越宽度 , 其差为 n - 1 d x 也就是说 , 双边移动观测模式跨越障碍能力比单边 移动观测模式要强一些。 4. 2. 3 “ X”形观测模式 图 15 为“ X”形模式参数计算图示 , 其中 B 1 \[N 1 - 1 L]x 2 18 B 2 \[N 1 - 1 L n - 1 d ]x 19 又有 B 3 B 4 \[N 1 - 1 L]x 2 n - 1 d x 20 B B 1 B 2 \[3 N 1- 1 L ]x 2 n- 1 d x 21 B Z k 2B 3 Z k \[N 1 - 1 L 2 n - 1 d ]x 22 化简整理式 21和式 22得跨越宽度 Z k \[ N 1- 1 L 2 - n - 1 d ]x 23 其余参数确定与对称排列中间放炮模式相同。 在不能放炮、不能埋检波器过障碍观测系统的条件 下 , 比较式 12、式 15 和式 23 可知 ,“ X”形模式 的跨越宽度最小。 图 15 “ X”形模式参数计算示意图 6 石 油 地 球 物 理 勘 探 2006 年 5 过障碍观测系统设计的基本要求 1尽量减小过障碍最小炮检距 , 使其降低到与 正常观测系统相同的水平。这是因为在观测系统设 计中 , 正常最小炮检距是按勘探目的层设计的 , 当测 线要跨越障碍物时 , 一般都要加大最小炮检距 , 若大 于目的层埋藏深度 , 就要以“牺牲”浅层资料质量为 代价来获得深层资料。 2选择良好的激发和接收条件。因为加大最小 炮检距后 , 增加了反射波射线路径长度 , 降低了反射 波能量和信噪比 , 故应尽量选择良好的激发和接收 条件 , 保证过障碍排列记录有足够的反射能量。 3因地制宜 , 灵活多变。在山地条件下 , 因地表 地质条件复杂多变 , 每个障碍物又有各自的特点 , 因 此设计过障碍观测系统应从实际出发。图 16 为通过 多个障碍物 果园、稻田、茶园及河流沙坝 时的观测 系统设计方案 , 其中过障碍观测系统连续采用中间放 炮对称排列观测模式 , 进行扩展和优化组合 , 穿越多 个无法放炮障碍区 , 笔者认为此设计方案是合理的。 图 16 川东云安场构造 A 测线观测系统 6 结束语 综上所述 , 二维地震过障碍观测系统设计不能 停留在定性阶段 , 有可能通过建立规范的观测模式 及科学的参数计算和实例分析 , 为设计野外过障碍 观测系统提供新的思路和方法。从本文提供的图示 和分析来看 , 只要针对具体的障碍分布和现有设备 情况 , 总可以科学合理地设计出与障碍区段相匹配 的过障碍观测系统 , 从而有可能用较低的费用获得 较好的勘探效果。 参 考 文 献 \[1 \] 邓志文 . 高陡逆掩推覆构造地区地震观测系统研究 . 石 油物探 , 2000, 41 2 127～ 131 \[2 \] 周绪文 . 反射波地震勘探方法 . 北京 石油工业出版社 , 1989, 102～ 126 本文编辑 刘英 7 第 41 卷 第 1 期 梁顺军等 二维地震过障碍观测系统模式及其参数设计 作 者 介 绍 梁 顺军 1957 年生 ; 1985 年和 1997 年先后毕业于西南石油 学院石油地质专业和物探专业 函授 。现在四川石油地球 物理勘探公司研究中心工作 , 并借聘到塔里木油田勘探开 发研究院 , 从事地震资料综合解释和研究工作 , 已发表和 宣读论文 10 余篇。 吴 艳辉 1981 年生 ; 2003 年毕业于内蒙古师范大学数学系 , 2003 年至今就读于大连理工大学数学系。目前在东方地 球物理公司物探技术研究中心从事技术协作。 杨 海申 工程师 , 1965 年出生 ; 1983 年毕业于石油物探学 校。长期从事于野外采集及静校正 , 现在东方公司物探技 术研究中心从事静校正方法研究和采集软件开发工作。 裴 正林 副研究员 , 1962 年生 ; 2000 年于中国地质大学 北 京 获博士学位 ; 2003 年于石油大学 北京 博士后出站。 主要从事井间地震层析成像、地震波传播理论及数值模拟 方法、小波变换应用研究。 马 德堂 讲师 , 1965 年生 ; 1991 年毕业于陕西师范大学数学 专业 , 获硕士学位。毕业后一直在西安地质学院从事数学 教学与研究工作 , 现为长安大学应用地球物理所工程地质 专业在职博士生 , 研究方向为波动方程正反演。 段 云卿 副教授 , 1958 年生 ; 1982 年毕业于华东石油学院物 探专业。发表论文多篇 , 现在中国地质大学 北京 从事教 学和地震资料处理及解释方法研究工作。 冉 建斌 高级工程师 , 1966 年生 ; 1988 年毕业于石油大学 东营 , 1997 年获石油大学 北京 煤田、油气地质与勘探 专业硕士学位 , 现在中国地质大学 北京 攻读博士学位。 先后主持和参与 30 多项生产和科研课题 , 获部级和局级 科技进步奖 10 项 , 发表论文 10 余篇。 贾 晓峰 博士后 , 1979 年生 ; 2000 年毕业于北京大学地球物 理学专业 , 获学士学位 ; 2005 年获北京大学固体地球物理 学博士学位。发表论文 5 篇 , 现在美国加州大学做博士后 研究 , 研究领域为勘探地球物理学。 闫 玉魁 高级工程师 , 1963 年生 ; 1984 年毕业于西南石油学 院地球物理勘探专业 , 获学士学位。长期从事油田勘探管 理工作 , 发表论文多篇 , 现在中油股份吐哈油田分公司勘 探事业部工作。 李 昌彪 博士生 ; 研究方向为计算智能与数据挖掘。目前在 西安交通大学电子信息工程学院攻读博士学位。 许 丽 高级工程师 , 1964 年生 ; 1986 年毕业于华东石油学 院勘探系。现在胜利油田地质科学研究院从事地震地质综 合研究 , 同时为中国石油大学 东营 在读博士生。 郭 平 高级工程师 , 1967 年生 ; 1989 年毕业于中国地质大 学物探专业。现在中油辽河油田勘探开发研究院计算所从 事地震资料处理及方法研究工作 , 同时在南京大学攻读硕 士学位。 王 西文 教授级高级工程师 , 1956 年生 ; 1982 年毕业于西安 地质学院 , 1987 年在该校硕士研究生毕业后留校任教 , 1989 年调入西北石油地质研究所 , 1997 年进入中国科学 院攻读博士学位。现在中石油勘探开发研究院西北分院地 球物理研究所工作。 金 龙 1978 年生 ; 2001 年毕业于大庆石油学院 , 现为中国 石油大学 北京 资源与信息学院博士研究生。研究兴趣包 括时移地震资料处理与分析、信息新技术在地震勘探开发 中的应用等。 吴 磊 1974 年生 ; 1997 年毕业于西南石油学院 , 现为中国 石油大学 北京 在读博士生 , 从事油气地质综合研究。 武 丽 工程师 , 1974 年生 ; 2000 年毕业于中国地质大学 北京 油气田开发专业 , 获硕士学位。现在中国石化石油 勘探开发研究院处理与解释中心工作 , 主要从事地震资料 解释与储层预测研究 , 同时在中国地质大学 北京 能源学 院矿产普查与勘探专业攻读博士学位。 王 在民 高级工程师 , 1943 年生 ; 1967 年毕业于南京大学数 学系数学专业。 1978 年北京大学数学系地震数据处理进 修班进修 , 现在新疆油田公司石油勘探开发研究院地球物 理研究所处理中心 , 从事方法研究和程序设计工作。 杨 辉 高级工程师 , 1961 年生 ; 2000 年毕业于同济大学 , 获博士学位 ; 2003 年从石油大学 北京 博士后流动站出 站。曾主持和参与多项省 部 级科研项目 , 获部级科技进 步奖两次 , 发表论文 20 余篇。现在中国石油勘探开发研究 院从事地质、物探综合研究工作。 底 青云 研究员 , 1964 年生 ; 1990 年毕业于长春地质学院地 球物理系 , 获硕士学位 ; 1990～ 1999 年在中国科学院地球 物理研究所工作 , 并于 1998 年获博士学位。目前主要从事 电磁波正、反演研究 , 合著专著 1 本 , 发表论文近 40 篇。 丁 建荣 高级工程师 , 1970 年生 ; 1993 年毕业于成都地质学 院应用地球物理系 , 获学士学位 ; 2004 年同济大学海洋与 地球科学学院毕业 , 获硕士学位。现在中石化江苏油田勘 探处从事勘探工作 , 并在中科院地质与地球物理所攻读博 士学位。 沈 建国 教授 , 1963 年生 ; 1983 年华东石油学院测井专业毕 业 , 2000 年获中国科学院声学研究所博士学位 , 2002 年从 清华大学物理学博士后流动站出站。现为天津大学教授 , 专门从事声波测井方法原理、资料处理和技术研究。 李 明杰 高级工程师 , 1964 年生 ; 1990 年获西安地质学院构 造地质学硕士学位。长期从事石油地质综合研究工作 , 发 表论文 10 多篇 , 现任东方地球物理公司研究院地研中心 总地质师 , 同时在中国地质大学 北京 攻读博士学位。 2006 年 2 月 石油地球物理勘探 第 41 卷 第 1 期