气体水合物反射层的速度结构研究.pdfVIP

维普资讯 第 1z卷 第6期 天 然 气 地 球 科 学 V01．12 N0．6 2001j 12月 NATURAL GAS GEOSCIENcE Dec． 2001 气体水合物反射层的速度结构研究 摘 要 地震反射删面的成像特征反映了海底模拟反射层(BSRs)与海底底面平 行特征，人们通常称谣反射屡就是甲烷水舍物稳定带的基底 本文根据地震反射数 据的波形反演技术，分析了用于估算BSR层的反射速度姑构特点和水合物层厚 度。通过大洋钻探计划(ODP)验证结果表 明，BSR层下部的沉积物孔隙层 中舍有 大约 30m厚的游离甲烷气 水舍物层和其下部的气体代表着 巨大的甲烷气藏，它 特有着重要的经济意义和影响全球气候变化 关键词 气体水舍物 反射层 速度姑构 BSR 波形反演技术 海底模拟反射层(BSRs)通常出现在大洋底都和相邻大陆边缘沉积物层的上部几百米 深处。最大可能出现在汇聚型大陆边缘加积堆积的棱柱体沉积物 中。BSR被广泛解释为 甲 烷水合物稳定带存在的基底 。在这样的环境下，甲烷水合物存在的稳定性主要受温度控制， 倒如 ，BSR常被用来估计地热梯度和热流产生的原 因。然而 ，BSRs最重要的意义是 由于在 水合物的内部和下部沉积层 中圈闭大量的甲烷气体 ，这具有重要的经济价值，同时水合物 的 形成和分解对全球气候变化也起着重要的作用 。 现在很多研究者用多道地震反射数据对 BSR上部 的水合物饱和层与其下部气体饱和 层的物理特性进行对 比，以建立BSR模型 但是，近来一些研究者认为，天然气的存在与否 并不需要波阻抗对 比，因为波阻抗对 比的产生是由BSR上部水合物沉积层和其下部的水饱 和带相互作用的结果，这个观点和建议是建立在收集了加拿大温哥华岛Cascadia沉降带的 1987年高质量地震反射数据并进行了详细研究的基础上(图1)。研究区的BSR反射层延伸 几万米。然而，常规合成地震剖面不能区分这两种模式，而通过钻井能达到区分的目的 ，如 1992年 1O月大洋钻探计划(ODP)164测线的钻探结果就证实了这点。为了识别BSR反射 层进行全球性钻撵是不允许的，也是不现实的。因此，用地震反射数据的波形反演能够区分 这两种模式并可提供 BSR下部游离气带的厚度 研究所使用的数据(图中圆点所示)来 自测线 89一O8，图中同时也标出了钻井位置 。 波形反演是最小化样品问的现场地震数据和台成地震数据的差值。由于波形差值函数 是非线性的，只有使用蒙特卡洛(随机)方法才可保证达到整体最小值(寻求一个单一的模 型)。如果未知函数的数据不是太大，这一方法将是可行的。同通常的地震反演情况一样，使 用蒙特卡洛方法在整个模型空间对单一函数进行最优化在经济上是相 当昂贵 的。然而，这个 问题可以用少参数函数 的蒙特卡洛技术得到解决 。这一技术使地震反射数据转换亏损得到 了补偿。我们首次使用随机检索方法估算速度变化的长波波长，然后用线性局部方法估算密 度和速度变化的短波波长。 我们将该反演技术应用于 Hyndman和 Spence地 区VI一5井的部分数据，该处海底和 BSR平行(图2)。反演分为三步： 第一步 ．我们通过在 120ms反射层的最大能量来估算三个突出的反射层 的速度平方 维普资讯 No．6 气体水舍物反射层的速度姑构研究 39 根。已确定的反射层的双程反射时间 为 2．23ms(图3)，其反射界面为BSR 层底界。我们计算出速度平方根的上 限和下限，用能量 曲率和最少能量来 求出单层层速度 。 从 图4F获得 的最佳反演模型， 时问轴上的箭头指示的第一步和第三 步反演的反射层。 第二步，同时对三个速度 的层速 度范围问的模型空间进行研究。开始 选用 5O个随机模型，而且每次迭代中 有～参数是随机变化的。如果能量随 着随机变化而增加，则该迭代模型即 为～个新模型。经 100次迭代后，该模 田1 加拿大盈哥华岛jI上Vl