第二章分辨率与信噪比的基本概念Word下载.docx

1、t （ 1）2.3f式中，f*是地震子波的视频率（或主频）。从这个公式出发，只要剖面上主频偏高，就认为 其时间分辨率高。如果把层速度V乘上tR，就得到厚度分辨率公式.*A V tR 丸Zr 二-Z - （ 2）2 4.6入为视波长，z-成为可分辨厚度。按这个公式，如果视频率为 30Hz,速度为3000m/ms,可分辨厚度为 22m这两个公式假定了地震子波是理想的雷克子波。 用楔形模型可以证明： 在上述分辨厚度时，相邻两个雷克子波的过零点互相重合，此时叠加的合成波形在两个波峰之间出现一个波谷，波 谷的谷底振幅为零，两个波峰刚刚分开，如图 11 （b）。一般来说，我们可以接受这两个公式，但请注意图

2、11 （ b）右方波形的振幅已经不真实，由 1.000下降到0.595。此外，如果实际子波不是零相位的，并且反褶积并没有把它压缩成一个正峰，则可分辨厚度比此公式肯定要差。例 如图11 （d）及（e）分别为最小相位子波及混合相位子波的例子，当子波错开半个视周期时， 其合成波形看不到相应的两个波峰，很难分辨谁是谁。所以，瑞利（Rayleigh ）准则指出：“一个反射波的分辨率的极限是 1/4波长。”注意|和曲乍量卞Hl应的带通f-诡/ K|UT和这里说的极限两字。公式（2）的Zr实际上也接近于四分之一视波长， 因而Knapp把这四分之一视波长定义称为“薄层分辨率”。由于实际地震子波并不是零相位的雷

3、克子波，公式（ 1 ）及（2）并不严格。所以我们应该进一步讨论一下对不同的子波具有不同频宽时，它们的特点是什么？三、三种相位特性带通子波的分辨率我用零相位、最小相位及混合相位三类带通子波，绘了一张不同子波的波形图，见图 12。从此图来看，第一，子波的分辨率首先与频带宽度有关。例如左边两个子波，由于频带很窄， 小于一个倍频程，则产生连续振动。此时，虽然视波长不是很大，但连续相位互相干涉，其分 辨率肯定不能再用公式（1）来表达。第二，不论什么样的子波，随着频带的加宽，子波的连续 振动相位数逐渐减少，并且主峰慢慢变窄、变尖锐，也就是分辨率逐渐变高，混合相位的子波 也是如此。第三，在频带窄的时候，最小

4、相位及混合相位子波与零相位的波形差别不大。但是 当频宽在一个倍频程以上的时候差别就大了。 例如零相位的10 80及10 160Hz的子波分辨率相当好，一个正峰十分突出，而最小相位子波变成正负相间的一个强头部，后面跟着一连串的 振荡波形。而混合相位子波，则表现为一个高频的头部和后面拖着一个低频的长长的尾巴！分辨率就降低了。所以说相同频带范围的子波中，以零相位子波的分辨率为最高。以后我还要提 到，实际子波都是混合相位的。四、频带上下限与分辨率的关系现在介绍一下频带上下限与分辨率的关系。1.绝对频宽如果我们将带通子波通频带的下限称为 fl,上限为f2,则fl-f2称为绝对频宽B，即绝对频宽B = f

5、 2 一 f 1 （ 3）对于零相位子波，绝对频宽决定了包络的形态。 Knapp在他的文章中对此作出了数学证明，他称之为“频移过程中子波包络的不变性”。例如1040Hz的子波包络与30 60Hz的包络完全 一样，因此他认为二者的分辨率也一样。见图 13上方，左右两个子波的 B都等于30,它们的包络就一样，30 60Hz的子波振动相位数却增加了两倍！2.相对频宽f2与f相除称为相对频宽 R即r f2R （ 4）f 1相对频宽还有一种表达方法，即采用倍频程数f 2OCT -3.32log 10（）一个倍频程就是音乐中的一个“高八度” ，也就是频率翻一翻。对于两个零相位带通子波而言，只要以倍频程表示

6、的相对频宽一样，那么，它们的波形就 是一样的，只是波形的胖瘦不一样，如图 13,左边上下两个子波，它们的相对频宽 R都等于4，于是它们的波形不变，只是包络缩放了一倍。通带为 20 80Hz的子波波形与1040的子波，只是在时间轴上压缩了一倍。 Knapp指出倍频程一样，波形一样时，还是瘦的子波分辨率最高。所以他认为分辨率不能用倍频程来衡量，只能用绝对频宽来衡量。相对频宽决定了子波的振动相位数。图 14是我统计的图12中零相位子波的连续相位数及主峰视周期宽度的示意图。 从图中可以看到零相位子波窄于一个倍频程时， 连续相位迅速增多，一个倍频程时相位数约等于 1.5个（此时用负极性显示，可见两个波峰

7、） 。所以，如果我们在地震剖面上看到大多数反射波有两到三个连续相位时，就可以说：它的分辨率不高，频宽不够一 个倍频程（或者是相位谱不好）。图丄4母相牠了-波的畅宽与HI位數蚕卞同關2：乾墓Number tjf phase- t即 &pTcmimnt pericd of a sere- phase wnvetetrefers tD different hfind - withh3.视频率通频带的中心频率fc决定了视频率fp （或称主频），即1/ 丄、fp 二 fc （ f 2 f 1） （6）2主频就是我们眼睛看到的剖面上同相轴的胖瘦程度。以上便是对带通子波的一般分析结果。五、关于分辨率的争议关

8、于分辨率的讨论其实还存在着不少争议。例如， Widess 认为在没有噪声的情况下，反射 波的可分辨厚度应该定为 1/8 视波长。又如 Knapp 指出，如果研究孤立的薄地层，其实不一定 需要有足够的带宽，例如大学物理教科书中就有使用单色光和相干原理测量肥皂泡薄膜厚度的 实验。Knapp 除了把四分之一视波长称为“薄层分辨率”以外，又划分了另外一种分辨率 韵律层的分辨率 ，即对薄层状的周期性重复的沉积韵律层来说， 可以用高频来检测薄互层的厚度， 这种韵律层产生一种对频率高次谐波的调谐作用。 我记得过去 50 年代老的地震教科书上也有一 段关于测定薄互层厚度的反射调谐振幅的文章。但是，以上测肥皂泡

9、厚度的例子必须是一层单 独的薄膜，而韵律层则必须是许多互层周期地作简单重复，它们都是在特殊情况下的特例。在 这些特例中可能不需要频带很宽，也不需要低频成分。 但是对我们地下的实际沉积剖面来说， 我认为：无疑地我们需要一个较宽的频带，才能解决分辨率问题。因为地的砂岩、泥岩、灰岩 的沉积往往不是一个“孤立”的层，也不是简单的周期重复，而是千变万化的。没有足够的频 宽是不能准确描述它的。六、不同频宽子波对不同厚度砂岩层的反映我认为好的分辨率不能光看绝对频宽是多少， 而还要看所占的频段在哪里 。例如 40 160Hz 的宽频（ B=120 ）不一定比 1080Hz（B=70 ）的好。因为前者缺乏 10

10、40Hz 低频段，会使厚 砂层在波阻抗反演后形成砂层中央的阻抗值下陷。为了说明不同频宽对描述砂岩分布的能力的 情况，我作了如下的试验。图 15 是不同频宽的零相位带通子波（即 Sinc 子波）反映砂层的反 射波形图。左方设计了从厚度为 3m到60m的7个砂层模型的反射系数，假定采样率为 2ms,层速度为3000m/s。右方各道是不同频带的子波与反射系数的褶积结果。这些带通子波的形态 可参看图20（ a）。从此图来看：第条曲线 40 50Hz的子波频带太窄，连续振动多，肯定不好；曲线及频带最宽， 为5 135Hz及10 160Hz,它们的反射波形都接近为反射系数曲线， 曲线及为10-23HZ及1

11、0 33Hz，频带虽然不太宽，但定性地尚可以反映出存在着 7个砂层。现在，我们把这个图的各个道积分一下，变成了相对波阻抗的形式，见图 16。左边第一道就是 7个砂层的波阻抗模型。现在我们来看右边哪一条曲线描述这个模型最好。显然, 55Hz及 10160Hz 两条曲线最好。右边曲线 1080Hz 也不错,只要砂层都是孤立的,互相不挨,就都能分辨，但它对6m到3m的砂层描绘得过分胖而没有棱角了。 再右边10 40Hz及10 23Hz曲线进一步圆滑,没有棱角,正峰宽度几乎相等,厚度也就成一致的了,同时负峰明显增强,10 23Hz的曲线连厚砂层也似乎变薄了。曲线 40 50Hz更不行，子波连续振动 3

12、.5个相位，分辨率最差, 看不清与砂层的关系。 而曲线及对厚砂层也反映很差。 40160Hz（B=120Hz）按分辨率定义来说分辨率很高,但由于缺乏低频,反映 10m 以上砂层的能力是很差的,此曲线对21m至60m砂层都出现中央波阻抗值下陷。 所以，绝对频宽很宽的子波如果缺少了低频， 也不是好事！七、不同的频率有不同的用处我认为不同的频率有不同的作用,其实道理很简单,在时间域,当子波的主瓣宽度（半周 期）和砂层的时间厚度相一致时,褶积后,输出振幅达到最强,否则振幅要变弱。这种效应相 当于一种对不同砂层厚度的“滤波器” 。L510123 sis SIS SI 口 31 K-r _!:-= - -

13、 WMLU一 H - - b -陌胡 为亘體狀雪苓盡直击邕T.慕屛尋琲霄冋弟器匕 産卿淙请虽satrid 烹 Ene -HycTed tlBdf- E. hxwurnhcMC:嗚cimi 诗一Ac 昨工田陰 嶋耳靈浅3翻盒耳eEsi辭隈弼筋曲辛虞即丞孟谯再囲注两彖 sejnfkronfryrred modelsverzl.on i !?!? nf 三 PA/R bnnd WHVeT茁显然，被增强的砂层厚度大致为 1/4视波长（即时间厚度为半周期），所以有4f表2增强砂层的厚度与加强频率的关系增强砂层的厚度（m）7537.52518.812.59.47.56.24.7加强频率102030405060708090例如，你希望搞清楚厚度为 60m的砂层，按表1, fi必须达到10Hz,否则如图16曲线，20 66Hz的曲线就不能反映出 60m的砂层，但是它还可以勉强反映 30m的砂层。又例如要使6m的砂层正确地反映出它的宽度和波阻抗数值，按附表 2则需要高频f2达120Hz，