第７卷第２期 ２０１０年４月 工程 球物理季赧 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＥＮＧＩＮＥＥＲＩＮＧ ＧＥＯＰＨＹＳＩＣＳ Ｖｏ１．７，Ｎｏ．２ Ａｐｒ．，２０１０ 文章编号：１６７２—７９４Ｏ（２Ｏ１Ｏ）Ｏ２一Ｏ１９０—０７ ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７２—７９４０．２０１０．０２．０１２ 基于小波域从叠前地震数据提取 地层品质因子Ｑ 袁恩辉，顾汉明 （中国地质大学地球物理与空问信息学院，武汉４３００７４） 摘 要：结合小波理论，利用地震波在黏弹性介质中的吸收衰减传播规律，在小波域以修正的复Ｍｏｒｌｅｔ小波 为基本小波，实现了从叠前炮集和共中心点道集地震数据，提取等效地层品质因子，然后通过剥层法估算层品 质因子。文中采用地面地震正演记录和实际海上地震资料试算，验证了本文提出的方法的正确性．结果表 明，方法正确可行，可靠有效，最后提出改善效果的措施和实际应用中应注意的问题。 关键词：品质因子；共中心点道集；小波变换；均方根品质因子；层品质因子 中图分类号：Ｐ６３１ 文献标识码：Ａ 收稿日期：２ｏｌｏ一０４一ｏ６ Ｅｓｔｉｍａｔｉｏｎ ｏｆ Ｑｕａｌｉｔｙ Ｆａｃｔｏｒ Ｑ ｆｒｏｍ Ｐｒｅ－－ｓｔａｃｋ Ｓｅｉｓｍｉｃ Ｄａｔａ ｉｎ Ｗａｖｅｌｅｔ Ｄｏｍａｉｎ Ｙｕａｎ Ｅｎｈｕｉ， Ｇｕ Ｈａｎｍｉｎｇ （Ｔｈｅ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ＆Ｇｅｏｍａｔｉｃｓ．Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｇｅｏｓｃｉｅｎｃｅｓ，Ｗｕｈａｎ ４３００７４，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｉｓ ｐａｐｅｒ ｕｓｅｓ ａｎ ｅｎｅｒｇｙ ａｔｔｅｎｕａｔｉｏｎ ｆｏｒｍｕｌａ ｏｆ ｓｅｉｓｍｉｃ ｗａｖｅｓ ｉｎ ｔｈｅ ｗａｖｅｌｅｔ ｓｃａｌｅ ｄｏｍａｉｎ ｆｒｏｍ ｔｈｅ ｗａｖｅｌｅｔ ｔｈｅｏｒｙ ａｎｄ ｓｅｉｓｍｉｃ ｗａｖｅ ｐｒｏｐａｇａｔｉｏｎ ｅｑｕａｔｉｏｎ ｉｎ ａｎ ｅｌａｓｔｉｃ ｍｅｄｉｕｍ．Ａｎｄ ｓｕｂｓｅｑｕｅｎｔｌｙ ａ ｆｏｒｍｕｌａ ｉｓ ｏｂｔａｉｎｅｄ ｔｈｒｏｕｇｈ ｗｈｉｃｈ ｃａｎ ｃａｌｃｕｌａｔｅ ｔｈｅ ｆａｃｔｏｒ ｉｎ ｔｈｅ ｗａｖｅｌｅｔ ｄｏｍａｉｎ．Ｕｓｉｎｇ ｔｈｅ ａｍｅｎｄｅｄ ｃｏｍｐｌｅｘ Ｍｏｒｌｅｔ ｗａｖｅｌｅｔ ａｓ ａｎａｌｙｚｉｎｇ ｗａｖｅｌｅｔｓ，ａ ｎｅｗ ｍｅｔｈｏｄ ｆｏｒ ｅｓｔｉｍａｔｉｎｇ Ｑ—ｆａｃｔｏｒｓ ｆｒｏｍ ＣＭＰ ｇａｔｈｅｒｓ ａｎｄ ｓｈｏｔ ｇａｔｈｅｒｓ ｉｎ ｔｈｅ ｗａｖｅｌｅｔ ｄｏｍａｉｎ ｉｓ ａｄｏｐｔｅｄ ｔｏ ｅｓｔｉｍａｔｅｓ ｔｈｅ Ｑ—ｆａｃｔｏｒｓ ｗｉｔｈ Ｉａｙｅｒ—ｓｔｒｉｐｐｉｎｇ ａｐｐｒｏａｃｈ．Ｂｙ ｓｙｎｔｈｅｔｉｃ ｓｅｉｓｍｉｃ ｒｅｃｏｒｄ ａｎｄ ｒｅａｌ ｓｅａ ｓｅｉｓｍｉｃ ｄａｔａ，‘ｗｅ ｅｓｔｉｍａｔｅ ｔｈｅ Ｑ ｗｈｉｃｈ ａｒｅ ｎｅａｒｌｙ ｔｈｅ ｓａｍｅ ａｓ ｍｏｄｅｌ ｐａｒａｍｅｔｅｒ ａｎｄ ｐｒｏｖｅｄ ｔｈｅ ｍｅｔｈｏｄ ｗｅ ａｄｏｐｔｅｄ ｉｓ ｃｏｒｒｅｃｔ．Ｄｅｓｉｒｅｄ ｒｅｓｕｌｔｓ ｈａｖｅ ｂｅｅｎ ａ— ｃｈｉｅｖｅｄ．Ｔｈｅ ｆｉｎａｌ ｒｅｃｏｍｍｅｎｄａｔｉｏｎｓ ｏｎ ｉｍｐｒｏｖｅｍｅｎｔ ｍｅａｓｕｒｅｓ ａｎｄ ｔｈｅ ｐｒａｃｔｉｃａｌ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｒｅｓｕｌｔｓ ｓｈｏＵｌｄ ｂｅ ａｔｔｅｎｔｉｏｎ． Ｋｅｙ ｗｏｒｄｓ：ｑｕａｌｉｔｙ ｆａｃｔｏｒｓ；ＣＭＰ ｇａｔｈｅｒｓ； ｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍａｔｉｏｎ；ｍｅａｎ。——ｓｑｕａｒｅ ——ｒｏｏｔ ｑｕａｌｉｔｙ ｆａｃｔｏｒｓ；ｉｎｔｅｒｖａｌ ｑｕａｌｉｔｙ ｆａｃｔｏｒｓ 波在其中传播时，介质中质点振动的能量因质点 １ 引 言 间摩擦等作用，有一部分转化为热能。波前扩散、 地层吸收和散射使得地震波能量随着传播距离的 实际地层介质是粘弹性非均质性介质，地震 增加而迅速衰减，特别是地层对地震波的高频成 基金项目：国家重点基础研究发展计划（９７３计划）（编号：２００９ＣＢ２１９４０３）资助。 作者简介：袁恩辉（１９７５一），男，工程师。硕士研究生，主要从事油气地球物理方法研究和油气储层预测。Ｅ－－ｍａｉｌ：ｙｕａｎｅｈ＠１６３．ｃｏｍ 顾汉明（１９６３一）・男，教授，博士生导师，主要从事地球物理勘探及资料处理解释的教学与研究。 第２期 袁恩辉等：基于小波域从叠前地震数据提取地层品质因子Ｑ １９１ 分的吸收，越到深处，对高频成分的吸收越厉害， 振幅按指数规律衰减，使地震波的深部分辨率降 ２．２小波尺度域能量衰减方程 对粘弹性介质中地震波传播方程进行连续小 波变换，小波采用（９）式的修正Ｍｏｒｌｅｔ小波，就可 以得到小波尺度域中的能量衰减密度方程，有关 低，地震子波的振幅谱和相位谱产生畸变。研究、 估算和补偿地下介质对地震波的吸收是叠前地震 资料处理的关键，它对提高地震成像精度、合理地 解释ＡｖＯ效应和正确地反演介质物理属性等都 是极其重要的。它能明显地提高地震资料的分辨 推导过程可以参阅文献１－６３，这里直接给出公式： Ｅ（ｃｃ，，ｚ）：２ａ√＿ｃ一 Ｕ（ ，０） Ｐ ， －－（ｏａ—ｍ） ／ｒ （５） 率，特别是中深部资料的分辨率，而要做好补偿， 式中 （２＞０）为地震波传播距离，【，（ ，Ｏ）为 一０ 关键在于准确地估算地层品质因子。 关于估算品质因子（Ｑ值）的方法，以前的学 者都做了大量的研究，主要有频谱比法，质心频移 法等，但大都借助于傅里叶变换，具有傅里叶变换 固有的缺陷限制。由于小波变换具有良好的局部 化性质，具有傅里叶变换所不具有的优势，很快 就被引入地震资料的处理中，在前人研究的基础 上，本文提出在小波域采用类似“频谱比”的方法， 由叠前炮集和ＣＭＰ道集估算Ｑ值。需要注意的 是，这里的“频率”与尺度有关，不是直接傅里叶变 换中的频率， ２理论基础 ２．１连续小波变换的定义 连续小波变换实信号，（ｆ）对于解析小波（￡） 的连续小波变换为 （硼 ，）（ ）：Ｊ口 ｆ ，㈤ （ ） （１） 式中￡为时间，窗函数 （ ）称为母小波，其上 面的短横“一”表示取复数共轭（下同），参数口，ｂ 分别为尺度因子和平移因子，（叫 ，）（口，６）称为 小波变换系数．母小波 （ ）的傅氏变换为 （０）， 母小波应当是绝对可积及平方可积，它满足允许 条件 Ｊｒ ＊ Ｉ∞ｌ ｄ∞＜＋∞ （２） 本文选择复数型小波的Ｍｏｒｌｅｔ小波： （￡）一Ｔ：－ ／ Ｐ拥 Ｐ一 。 （３） 式中７ｎ为调制频率， ≥５，为了使小波与地震子 波的波形特征相接近，可对Ｍｏｒｌｅｔ小波进行修 正，得到修正后的Ｍｏｒｌｅｔ小波 （ｆ）＝丌＿ ／ ‘８一‘ ‘／。 （４） 式中Ｃ为调幅因子，它控制了小波函数的长度，参 数的具体选择见文献Ｅｓｌ。 时的震源波场，ｃ（ｏＪ）为相速度。 对角频率，假设震源子波为理想的脉冲源，即 己，（ｃｃ，，０）。＝１，且忽略由频散效应引起的衰减，可 以推导出小波尺度域能量衰减公式 Ｅｄ：ｃ Ｐ （６） 式中Ｅａ为尺度能量．由于品质因子Ｑ》１，选取 合适的Ｃ，使得Ｃ２／Ｑ ａ 足够小，忽略上式中的平 方项，则式（６）简化成 Ｅ＝ｆ一 （７） 如果在ｔ 时刻的尺度能量谱为Ｅ。 ，在ｔ 时 刻的尺度能量谱为Ｅ ，则可以得到 Ｅｌ一 ：一 Ｃ （８） 对式（８）取对数，然后可用一条直线拟合， 其斜率是有效品质因子的函数。利用该式估计的 Ｑ值是地层的等效品质因子Ｑ值，要求取单套地 层的Ｑ值，有 —Ｑ 一 Ｑ口， 。＋ Ｑ （９）… 式中，ＱＤ，，，Ｑ 分别是反射时间ｔ。，ｔ。时刻的等 效品质因子，Ｑ。为ｔ。与ｔ 之间的品质因子． 接着可以采取剥层法，再利用张长军提出的 方法，采用类似于层速度分析的方法，得出计算层 Ｑ值的公式如下 （１０） 式中，Ｑ 是层ｉ在不同偏移距计算的Ｑ值的平 均，ｔ。． 是反射界面ｉ的零偏移距旅行时。 ３正演模型计算 ３．１水平地质模型的炮集模拟记录计算 图１是设计的吸收衰减地质模型，模拟海上 地震勘探，模型参数如表１所示。图２是利用粘 弹性波动方程模拟产生的单炮记录，采用和野外 一致的观测系统参数进行模拟，雷克子波主频 １９２ 工程地球物理学报（Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｎｇｉｎｅｅｒｉｎｇ Ｇｅｏｐｈｙｓｉｃｓ） ∞ Ｏ Ｏ Ｏ Ｏ 第７卷 ｌ ｌ １ ｌ ０ ４５Ｈｚ，检波器间距为３０ｍ，炮间距１２０ｍ，最小偏 移距６０ｍ，排列长度３６００ｍ，在海面以下１０ｍ处 ０．２ 激发，网格１０ｍ×１０ｍ。 ０．４ 表１模型参数 ０．６ Ｔａｂｌｅ １ Ｍｏｄｅｌ ｐａｒａｍｅｔｅｒｓ ０．８ 介质 ｍ 品质因子Ｑ 舶 厦 １ ‘ 留 Ｉ．２ Ｉ．４ Ｉ．∞ ６ ∞：。∞ 如 ｘ／ｍ １．８ ２ ５Ｏ ｌＯＯ １５Ｏ 频率ｆ／Ｈｚ 图４炮集记录第６３道及其小波时频谱 图１地质模型 Ｆｉｇ．４ Ａ ｔｒａｃｅ ６３ ｏｆ ｓｈｏｔ ｇａｔｈｅｒｓ ａｎｄ ｔｈｅ Ｆｉｇ．１ Ｅｘａｍｐｌｅ ｏｆ ｆｏｕｒ－－ｌａｙｅｒ ｍｏｄｅｌ ｔｉｍｅ－ｆｒｅｑｕｅｎｃｙ ｄｉｓｔｒｉｂｕｔｉｏｎ 图２是正演模拟产生的地震记录，图３是经 过球面扩散补偿，动校正后的炮集记录，图４是 经过动校正后的炮集记录，图３中第５３道数据的 时域波形，采用前面提到的修正的Ｍｏｒｌｅｔ小波作 为母小波，对该道信号进行小波变换，在小波域 计算得到的尺度振幅包络剖面（图４右图）。其中 调制频率优一１００，ｃ一３０．（下面的数值试验中， 点号 都采用同样的小波参数），尺度和频率有一定对应 关系，本文已将尺度换算成了对应的频率。 图２单炮地震记录 Ｆｉｇ．２ Ｓｈｏｔ ｇａｔｈｅｒｓ 图５是进行线性回归估计算第１层Ｑ值时， 拟合选取的“伪尺度”（时间／尺度）范围为４～６ （图５（ａ）），对应的频率为６０￣９０Ｈｚ（图５（ｂ）），利 用公式（８）估计Ｑ值，得到反射时间ｔ 一Ｏｍｓ到ｔ。 一７８４ｍｓ处，即地震波在海水中传播吸收衰减的 品质因子Ｑ＝５０９．４５，与实际海水很少的吸收衰 减情况比较吻合，但达不到理论给定的２００００。 图６是进行线性回归估计算第１层Ｑ值时， 拟合选取的“伪尺度”（时间／尺度）范围（图６（ａ）） 在整个频带范围（图６（ｂ））内的分布，从中可以看 点号 出尺度振幅包络对数曲线直线反映了模拟的地震 图３经过动校正后的炮集记录 波在海水还是有吸收衰减的，而且符合高频衰减 Ｆｉｇ．３ Ｓｈｏｔ ｇａｔｈｅｒｓ ａｆｔｅｒ ＮＭＯ 的快，低频衰减的慢的吸收衰减规律。 ： 第２期 瓣靛苍删 一 一 一 一 一籁靛 栏 籁茛盛删箍 袁恩辉等：基于小波域从叠前地震数据提取地层品质因子Ｑ 一一 一 一一 一瓣窝 哪媪 １９３ ４ ５ ６ ７ ８ ＂ ９６ 、 ～， １ 、、！ ＼● 、 — 、 ＼ —、 （ａ） 计算结果Ｑ＝５Ｏ９．４５４１ （ｂ） 图５线性回归 Ｆｉｇ．５ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉｏｎ ｏｆ Ｑ ｆａｃｔｏｒ ＿ 、 。．． 』 ．． － …一’翟ｒ … （ａ） （ｂ） 图６线性回归 Ｆｉｇ．６ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉｏｎ ｏｆ Ｑ ｆａｃｔｏｒ 下面是对该ＣＭＰ道集进行计算，提取的品 质因子（表２），从中可以看出，估算的结果与实际 模型给定的理论值吻合得比较好。 计算结果如表２。 表２计算得到品质因子与地质模型参数对比 Ｔａｂｌｅ ２ Ｑｕａｌｉｔｙ ｆａｃｔｏｒ ｃａｌｃｕｌａｔｅｄ ｂｙ ｃｏｍｐａｒｉｓｏｎ ｗｉｔｈ ｔｈｅ ｇｅｏｌｏｇｉｃａｌ ｍｏｄｅｌ ｐａｒａｍｅｔｅｒｓ ３．２起伏地质模型的ＣＭＰ模拟记录计算 图７是设计的带吸收衰减起伏地质模型，模 拟海上地震勘探，模型参数、地层Ｑ值、地层密 度、各层速度如表１所示。图８是利用粘弹性波 动方程模拟产生的单炮记录，采用和野外一致的 观测系统参数进行模拟，其中激发震源采用雷克 子波主频５ＯＨｚ，检波器间距为３０ｍ，炮间距 １２０ｍ，最小偏移距６０ｍ，排列长度３６００ｒｎ，在海面 以下５ｍ处激发，网格ｌＯｒｅ×ｌＯｍ。 表３模型参数 Ｔａｂｌｅ ３ Ｍｏｄｅｌ ｐａｒａｍｅｔｅｒｓ ｘ／ｍ Ｏ 量Ｉ ２ 图７地质模型 Ｆｉｇ．７ Ｅｘａｍｐｌｅ ｏｆ ｆｏｕｒ—ｌａｙｅｒ ｍｏｄｅｌ ０．Ｏ０ ０．Ｏ０ ０．２５ Ｏ．２５ ０．５０ ０．５Ｏ ∞ ０．７５ Ｏ．７５ １．Ｏ０ １．Ｏ０ １．２５ １．２５ １．５０ １．５０ １．７５ １．７５ 点号 图８单炮地震记录 Ｆｉｇ．８ Ｓｈｏｔ ｇａｔｈｅｒｓ 图８是正演模拟产生的地震记录，图９（ａ）是 ｅｄｐ号为５１０￣５１４的共中心点道集速度谱，图９ （ｂ）是共中心点５１０道集，图９（ｃ）是ＣＭＰ５１０中 第八道数据的时域波形，采用前面提到的修正的 Ｍｏｒｌｅｔ小波作为母小波，对其进行小波变换，在 小波域计算的尺度振幅包络剖面。 工程地球物理学报（ｃｈｉ？？竺 ！：：！！ ！ 竺！ ：！：竺！！竺！！竺竺 ———— １９４ Ｃ ＣＬ＾ｕ 正 集 墓 畦 地 Ｆ 速 度 第卟 八 数 、Ｊ ｂ 口 ｑ Ｍ 陂 蝶划 顺 螽４３・８ 表４计算得到品质因子与地质模型参数 Ｔａｂｌｅ ４ Ｑｕａｌｉｔｙ￣ａｃｔｏｒ ｃａｌｃｕｌａｌｅｄ ｐ ｍ 。 ｗｉｔｈ ｔｈｅ ｇｅｏｌｏｇｉｃａｌ 捉４３．６ 警 嘲４３．４ ４３． ４３ 度Ⅱ ０度 小 振 。晏 波 幅 时 Ａ 燕 窝 甾 嘲 越 魁 艚 刹 面 ４海上实际资料计算 图１０线性回归 Ｆｉｇ．１０ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉ。ｎ。ｆ Ｑ ｆａｃｔｏｒ 振 图１１海上单炮记录 Ｆｉｇ．１１ Ｓｈｏｔ ｇａｔｈｅｒｓ 第２期 袁恩辉等：基于小波域从叠前地震数据提取地层品质因子Ｑ １９５ 图１２海上ＣＭＰ道集 Ｆｉｇ．１２ ＣＭＰ ｇａｔｈｅｒｓ 海洋实际信号 小波时频图 望 厘 蓄 １ Ｆ １ ＿ ｒ Ｉ ｊ ｜ ＇｜｜ ｒ －葺ｉ ｝ ３ ｊｉ ｒ ５Ｏ ｌＯＯ ｌ５Ｏ 频率ｆｌＨｚ 图１３海上ＣＭＰ道集第６３道及其小波时频谱 Ｆｉｇ．１３ ＣＭＰ ｇａｔｈｅｒｓ ａｎｄ ｔｈｅ ｔｉｍｅ－－ｆｒｅｑｕｅｎｃｙ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｔｒａｃｅ ６３ 图１４是进行线性回归估计算第１层Ｑ值 时，拟合选取的“伪尺度”（时间／尺度）范围为２０ ～２５，对应的频率为１４５～１８５Ｈｚ，（海上地震采 集加高切为２０６Ｈｚ）进行线性拟合Ｑ值的结果。 利用公式（８）估计Ｑ值，得到反射时间ｔ ＝０ｍｓ 到ｔ：一２２１９ｍｓ处，即地震波在海水中传播吸收衰 减的品质因子Ｑ＝１２９７．６２，与实际海水吸收衰 籁 韶 捌 疆 ｔｉｍｅ／ａ （ａ） 图１４ 线性回归 Ｆｉｇ．１４ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉｏｎ ｏｆ Ｑ ｆａｃｔｏｒ 、＼ｋ 、 －： ｛ 、 ｑ ‘、～ 卜 ｔｉｍｅ／ａ （ａ） 图１５ 线性回归 Ｆｉｇ．１５ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉｏｎ ｏｆ Ｑ ｆａｃｔｏｒ 减很少的情况比较吻合。 图１５是进行线性回归估计算第１层Ｑ值 时，拟合选取的“伪尺度”（时间／尺度）范围（图１５ （ａ））在整个频带范围（图１５（ｂ））内的分布，从中 可以看出尺度振幅包络对数曲线直线反映了地震 波在海水还是有吸收衰减的，而且符合高频衰减 的快，低频衰减的慢的规律，只不过非常的少。 图１６是进行线性回归估计算第７层Ｑ值 时，拟合选取的“伪尺度”（时间／尺度）范围为２ｌ ～３（图１ ６（ａ）），对应的频率为７ ０～１ ０ ０ Ｈｚ １９６ 工程地球物理学报（Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｎｇｉｎｅｅｒｉｎｇ Ｇｅｏｐｈｙｓｉｃｓ） 第７卷 介质中的吸收衰减传播规律，在小波域以修正的 籁 靛 盛 一 翟 、 、 复Ｍｏｒｌｅｔ小波为基本小波，实现了从叠前炮集和 —冉 共中心点道集地震数据，提取等效地层品质因子， 卜、 — 一 然后通过剥层法估算层品质因子。文中采用地面 地震正演记录和实际海上地震资料试算，验证了 本文提出的方法的正确性．结果表明，方法正确 可行，可靠有效，最后提出改善效果的措施和实际 应用中应注意的问题。 ２）计算层Ｑ值，沿目的反射层进行，注意选 择资料信噪较高地震资料，尽量避免多次波，严重 绕射波的影响，研究中发现这些干扰波对品质因 子的提取影响特别大，尤其是针对起伏较大的地 图１６ 线性回９ｊ 层，更是如此。 Ｆｉｇ．１ ６ Ｌｉｎｅａｒ ｒｅｇｒｅｓｓｉｏｎ ｏｆ Ｑ ｆａｃｔｏｒ ３）在实际计算层Ｑ值时，要注意做好静校 正、扩散补偿，速度分析等预处理工作，在此基础 （图１６（ｂ）），进行线性拟合Ｑ值的结果，计算得到 上进行估算Ｑ值，效果更佳，建议在剩余静校正 第７层岩层的品质因子Ｑ一１３７．２。 后进行。 对其中第６６道做变换，进行了品质因子提 ４）本文针对的都是厚层的地质模型，而且起 取，估算的Ｑ值曲线如图１７所示，由于该区岩性 伏曲率不太大的情况，至于薄层或者起伏很大的 主要为沉积岩类，故可以根据品质因子做大体推 地层的情况还有待于进一步的研究。 断，结果如下： 参考文献： Ｏ～２２１９ｍｓ，Ｑｌ＝５００以上，为海水；２２１９～ ２２７５ｍｓ，０．２＝２Ｏ．３，为第四系海底松散物，可能以 ［１］Ｓｈｅｒｉｆｆ Ｒ Ｅ，Ｇｅｌｄａｒｔ Ｌ Ｐ．Ｅｘｐｌｏｒａｔｉｏｎ Ｓｅｉｓｍｏｌｏｇｙ 黏土之类为主的物质；２２７５～２５４０ｍｓ，Ｑ一１２８．１， （Ｓｅｃｏｎｄ Ｅｄｉｔｉｏｎ）Ｉ－Ｍ３．Ｎｅｗ Ｙｏｒｋ：Ｃａｍｂｒｉｄｇｅ Ｕｎｉ— 该层较厚，可能为泥岩粉砂岩互层；２５４０￣２６４８ｍｓ， ｖｅｒｓｉｔｙ Ｐｒｅｓｓ，１９９５． ＱＩ＝５０．６，以粉砂泥岩为主；２６４８￣２７８２ｍｓ， ： ［２］Ｈａｕｇｅ Ｐ Ｓ．Ｍｅａｓｕｒｅｍｅｎｔｓ ｏｆ ａｔｔｅｎｕａｔｉｏｎ ｆｒｏｍ ｖｅｒｔｉ— ２７．２，粉砂质泥岩、泥岩互层为主；２７８２￣２８６０ｍｓ， ｃａｌ ｓｅｉｓｍｉｃ ｐｒｏｆｉｌｅｓｌ，Ｊ－］．Ｇｅｏｐｈｙｓｉｃｓ，１９８１，４６（１１）： １５４８～１５５８． Ｑ６—４６．７，可能是中砂、粉砂岩互层；２８６０～ ［３］Ｓｔａｉｎｓ ｂｙ Ｓ Ｄ，Ｗｏｒｔｈｉｎｇｔｏｎ ＭＨ．Ｑ ｅｓｔｉｍａｔｉｏｎ ｆｒｏｍ ３３２９ｍｓ， 一１３７．２，砾状砂岩中粗砂岩等。 ｖｅｒｔｉｃａｌ ｓｅｉｓｍｉｃ ｐｒｏｆｉｌｅ ｄａｔａ ａｎｄ ａｎｏｍａｌｏｕｓ ｖａｒｉａｔｉｏｎｓ ｉｎ ｔｈｅ Ｎｏｒｔｈ Ｓｅａ口］．Ｇｅｏｐｈｙｓｉｃｓ，１９８５，５０（４）：６１５￣６２６． ［４３ Ｒｉｃｋｅｒ Ｎ．Ｔｈｅ ｆｏｒｍ ａｎｄ ｎａｔｕｒｅ ｏｆ ｓｅｉｓｍｉｃ ｗａｖｅｓ ａｎｄ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ｏｆ ｓｅｉｓｍｏｇｒａｍｌ，Ｊ３．Ｇｅｏｐｈｙｓｉｃｓ，１９４０，５ （４）：３４８～３６６． ［５］高静怀，汪文秉，朱光明，等．地震资料处理中小波函数 的选取研究［Ｊ］．地球物理学报，１９９６，３９（３）：３８９￣３９７． ［－６３李宏兵，赵文智，曹宏，等．小波尺度域含气储层地震波 衰减特征［Ｊ］．地球物理学报，２００４，４７（５）：８９２￣８９８． ［７３高静怀，杨森林．利用零偏移ＶＳＰ资料估计介质品 质因子方法研究Ｉ－Ｊ－Ｉ．地球物理学报，２００７，５０（４）： 图１７ 估算的品质因子 １１９８～１２０９． Ｆｉｇ．１ ７ Ｅｓｔｉｍａｔｅｄ ｖａｌｕｅｓ ｏｆ Ｑ ｆａｃｔｏｒ ，１８３陈爱萍，邹文，刘天佑．共中心点道集中品质因子的 估计及应用Ｉ－Ｊ３．勘探地球物理进展，２００４，２７（１）： ５ 结 论 ３２￣４４． ［９］赵伟，葛艳．利用零偏移距ＶＳＰ资料在小波域计算介质 １）本文结合小波理论，利用地震波在粘弹性 Ｑ值Ｄ］．地球物理学报，２００８，５１（４）：１２０２￣１２０８．