井间地震弹性波传播规律研究.docxVIP

井间地震弹性波传播规律研究 与地面地震相比，井间地震勘探具有以下特点：区域内不受低速带的影响，地震波传播距离短，源激发频率高（通常100.100m），接收点和源点可以接近研究目标，接收到的地震记录频率高（地震波的10倍或10倍），能量强，信噪比高。波场的类型很多，复杂的波散射波能量很强。可以高区分波、射波、散射波、导波、破断波、转换波和反射波等波型。因此, 井间地震是一种有效的地球物理技术, 它不但在地震学研究, 煤、金属与非金属等固体矿产资源, 工程地质, 环境地质和地质灾害勘探中得到卓有成效的应用, 而且在石油天然气勘探与开发中获得成功应用 由于井间地震波场类型多, 而且各种类型波之间互相叠加或干涉, 形成了非常复杂的弹性波场, 这对于井间地震波场识别与分离以及成像带来了很大困难。研究井间地震弹性波场特征和传播规律将具有非常重要的理论和实际意义。 井间地震与地面地震的激发、接收方式不同, 因此, 在地下介质中它们的波场特征和波场成分均存在差异。张碧星等 井间地震波场正演模拟是研究复杂地层中井间地震波传播特征的有效手段。目前, 地震波场数值模拟方法可以归纳为3大类, 即地震波动方程数值解法, 积分方程法和射线追踪法。对于一些复杂的本构方程, 由于积分方程法和射线追踪法不满足假设条件, 因而限制了其应用。地震波动方程数值解法自20世纪60年代以来得到了飞速发展, 已形成了目前具有有限差分法、有限元法、虚谱法等各种数值模拟方法的现代地震波动方程数值模拟技术。这些解法各有优缺点。对于复杂构造、复杂地质体和复杂岩性地震波场数值模拟而言, 交错网格高阶有限差分法的综合性能 (占内存大小、模拟精度、计算效率和并行算法实现) 最好, 并且最具实用性 笔者将采用弹性波方程交错网格高阶有限差分数值模拟方法, 重点研究基础性模型的井间地震弹性波场特征及传播规律, 为井间地震资料处理、成像及解释提供理论依据。 1 反射p波入射角的差异 为了分析固-固界面上的地震波场, 设计如图1所示的2层水平层状模型。模型大小为100 m×200 m。网格大小为1 m。数值模拟中, 模型四周为 完全匹配层 (PML) 吸收边界 图2所示为2层水平层状模型井间地震不同炮点弹性波场记录。从图中可以看出:①当震源为P波源时, 在固-固界面上的井间地震弹性波场主要包括直达P波, 反射P波, 反射PS波 (转换波) , 透射P波, 透射PS波 (转换波) , 首波 (P波和S波) 和面波等;②随着炮号的增大, 即P波入射角的增大 (其变化范围为十几度至90°) , 井间地震弹性波场发生了明显变化, 其中, 反射P波的振幅变化最大。反射P波、反射PS波、透射P波、透射PS波的振幅变化与根据Zoeppritz方程所计算的在弹性界面 (速度界面为3 800 m/5 600 m) 上, 当P波入射时, P波、PSV波的反射和透射系数相吻合 (对比图3) 。 为了便于进一步分析P波入射角小于临界角和大于临界角时弹性界面上反射波场特征, 从图2的12炮记录中抽出了反射点在界面中心处的共反射点记录 (图4) , 共12道。记录的入射角见表1。 从图4可以清楚地看出:①当入射角增大但小于临界角 (42.7°) 时 (见图4中1～4道) , 反射P波振幅增幅不大, 而反射PS波 (转换波) 增幅比较大;②当入射角大于临界角后, 反射P波增幅快速增大接近直达波振幅 (见图4中6～10道) , 而反射PS波振幅在第6道达到最大, 之后其振幅逐渐减小 (图4中7～10道) ;③对比1～4道和5～10可以看出, 当小于临界角和大于临界角的反射P波相位发生了明显变化 (直达波为最小相位) , 即从最小相位到0相位再到最大相位。图4中的第4道反射P波的相位接近0相位, 而反射PS波 (转换波) 无此现象 (由于上层纵波速度大于下层横波速度而不会出现临界反射之故) 。 图4b的第9～11道上的低频波成分中, 第10道上面波最强, 离开界面能量迅速减小。该面波是当震源位于界面上或附近时, 沿界面传播的瑞利面波传播到接收器一侧后, 在界面处重辐射 (二次源) 所形成的一类面波 (或管波) 。 综上可得:①井间地震入射角变化范围广 (二十几度至90°) , 反射波多为大角反射;②当P波入射角大于临界角后, 反射P波的相位发生了明显变化 (与入射波相比, 或与小于临界角反射波相比) 。这正好解释了复数反射系数的物理含义。 2 高速夹层水平夹层模型 导波在井间地震勘探中占有重要地位。下面通过数值模型模拟结果来分析高速层和低速层中导波的不同特性。 低速夹层水平层状模型取泊松比为0.25, 模型大小为200 m×200 m, 夹层厚度为30 m, 网格大小为1 m (图5) 。数值模拟中, 模型四周为PML吸收边界