第十章高分辨率地震勘探剖面的解释

1、第十章高分辨率地震勘探剖面的解释高分辨率的剖面应该如何来作解释？这个问题过去没有人认真考虑过，不少研究单位直接用高分辨率的剖面去作砂层追踪，或者使用其振幅参数，用调谐曲线作厚度计算。这是有问题的。我作了不少理论记录，发现反射记录上的同相轴实质上是追踪着反射系数而不是追踪着砂岩。在分辨率不高的常规剖面上，这个矛盾暴露得不明显。 当有效频带展宽到 80Hz以上，主频达50Hz以上时，反射剖面愈来愈不像砂层剖面。而当主频达1OOHZA上时，反射剖面实际上是接近为反射系数剖面，人们很难看懂它。只有把这种剖面转换成波阻抗剖面，或者简单地从上到下连续积分累加，才能真实地反映砂层的变化情况。第一节一个完整砂

2、岩体的模型试看我作的理论地震记录：图 83 (a)是设计的一个地下砂层分布的剖面。此剖面是声 速剖面。图中黑色的是砂层，每个道的基线就是泥岩基线，它相当于速度为3000m/s。两个道之间的距离相当于速度增量为100Om/s。所以砂岩的速度在 3400- 390Om/s之间变化。采样率为1ms,设地层平均速度为 300Om/s,则每个采样点相当于 1.5m厚度。此剖面的右方砂层较多、较厚，相当于一个砂岩体的根部。而向左砂层逐渐减薄，层数增加，相当于砂岩体的前缘，(三角洲相)。再向左边进入浅湖相，砂岩薄而且层数稀少。因 而此模型可以表达一个完整的砂岩体的情况。整个剖面的总厚度是 200ms,相当于

3、地下300nl在右方单层砂岩最厚的达 33m,到左方最薄的可以为 1.5m (一个采样点)。我们从Gardner经验公式出发即：总体密度P = 0.31v°.25( 69)v为纵波层速度，就可以推算出波阻抗p. v。从而得到反射系数曲线如图83 (b)。图83 (c)是同一个模型用另外一种显示方式，内容与(4)相同。它是以曲线的形式表达的声速模型。在图 83 (d)中列出了我们作试验中所采用的不同的子波形态(见注释)。让我们先来看看零相位带通子波的理论模型。将子波褶积反射系数模型图83 (b)就可得到如下理论的合成地震剖面。图84 (a)是10 20HZH个倍频程的带通子波(即 Si

4、nC子波)的例子。这个频档大致 相当于我们新生界盆地中常规剖面3s左右深层的反射频带。图的右方就画出了这个子波的形态。与图83 (a)答案相比，这个例子的分辨率太差了。基本上看不清地下砂层的分布情 况。将答案透明薄膜图 83 (a)上移20msM,可以看到左方 30m 33m及27m三个厚砂层在 图84 (a)上有三个波峰就正好吻合。(上移的原因是由于反射波从砂层顶开始起跳)。但其它的砂层都对不上。图 84 (a)中凡是有强振幅处，一般地反映它附近有地下较厚的砂层。 而薄砂层全部没有反映。图84 (b)是零相位带通子波 10-40Hz的两个倍频程的情况，子波绘于右方。这是我 们常规地震剖面1s

5、到3s之间典型的通频带。注释：(本书将图83 (a)及(b)用透明薄膜印了一份红色的波形图，是为了让读者可以剪 下来蒙在其它几幅图上自己去作对比。)i, tiffirifril4 ILL 尸 p J多大手 Lrrf;,二卜 AL L， 二；Q I v -， -1a b . 二匚匚匚匚 WTT-4弓 &LLL1JL. JWZJM. thrtar d，TW-I一 ) 二二二, 1 二一 rfT rctjrr. 4 my r 4H弓一二J一 二. .二 一 J"=rrb :上 J UTf>一riku尸. - -7 二二三：二：'1r二,根据力速模型推亮.假定密度与速度

6、之间有G.Mnc,公式.P=0.3|P0图83(a) Theoretical sand body- Sonic logs model with sample rate 1 ms.Thickness of sandstone varies from I. 5m. to 33 meters.(b) Reflection coefficient model derived from sonic logs usingGardner'、empirical formula.m m3m 324733£? L/u Jr户，LEL A =丁Larm;丁LQrmxr3 ;l JaL n La N

7、 ,二L刍彳丫L , L-cr-;.;,二 '.-K l"l :h/L,二，i<rns.L ;多.；w;yl. I.d产 LZ、，、L二，；士,工j>广 LJ.L .JC: ；.>JLr LVJ w 匚厂 nLaMs(EJl 2二Gr'LJ之二7 卜.5.hkrAL L X ,>k: %Lrv%，LA 之1工r工1at CQrA之三.h白丁f7:止；eL*AYF广言户广广r . A虻、wrf /之"rkz叼勿。7, 二.一号吟ju,uNr:七互-土"匕，尸厂” h/I L.W.r'J frr :；rrr三it-LJe

8、r-A4£.； :匚>丁 HJ'rkccJ/.J, r ,>?r:qeLA产二 £匚 J幺 I NG为ZJJJCUUU 二L j 二、>>二匚尸二sLUUVLrG7mt, JLQ£lr.4/匚I- dqNkL,c2二teJ; 艺"亍户乏力七支二2歹Kkr ,云二->尸一七七.； 之二1' z ，GZ 上，：/,r ( i-lc>h><rcr 1子波经常通滤波后 的形态B328 35l70l90HZ 带通滤波»i克劳德F波%相位带通子波般小相位的帚通f ffi雷克于波算通频带恒为12

9、5必L 一 一”，一n力相位校正及 极性次正后的子波Irmo.采样率全为ImOI-6OHZo180 Hl0140 千来 fhj.kJ(d)试验中使用的各神子波图83地下个完整砂岩体的声波速度模型(0)Same as Fig. 83 (a) different display mode (wiggled trace).(d) The wavelets used in this test. Ricker's. Klauder*s & Min. phase,-零相位带通f波/ -10-20+的理论合成她震道别由 这肿只有一个倍赖程的低糖记求是很睢反映心必的变化的I。.15卜，零相位带

10、通f波/ =1。-小用,的理论介成她震迤副血这是我们通常所树到的臃带为10 402的反射记录.它不能细致地反映薄砂片伸.它大致说明“砂奴介何处(a) Syutluuic mmsouc refl. srcfion ： liand width 10 20Hz. (h) liand width 卜)10Hz.十二-J-q程二至空孑全海吟笆秋残常爱足*烷律行而料会盅工装后4金皿豪域尼把笈裆零型liiZH图三居.唐娈.+f 5'二.一一. 4 _ . _ / a *_* _l£l丁 一 r C r ft r u T ，, 下 T , _ 5?f*MYYVviAY7724W* ,.丁 1

11、*寸布*廿isrMrVf忙汽筑、建:她山婢簿/这一 十 JxaXiAaZaAJXlV将此图84 (b)与答案83 (a)核对时，应该将 83 (a)上移约18mso右方33m及27m 砂层有较好的反映。其它厚砂层也都有些反映，但反映不真实。例如我们打X号之处就是地 下有10来米砂层而反射波形反而变为波谷或者压在过零点上。图的最左边1 8道在浅湖相沉积中的6m砂层却有着良好的反映，而且振幅相当的强。这便是为什么通常在大套泥岩的 湖相沉积之中追踪薄砂层常常能够取得成功的道理。从第8道至第20道之间有三个倾斜的同相轴是复合叠加所造成的“ 视同相轴”。在我们模型中其实并没有倾斜的地层，这种现象也值得我

12、们在解释工作中加以注意！我在陆相沉积地震地层学若干问题一文中曾经指出这种视同相轴有时会造成假的叠瓦、上超或下超现象。图84 (c)是带通子波108OHz三个倍频程的情况。这相当于目前石油勘探(中层反 射)的高分辨率的地震剖面。此剖面的主频约65Hz左右。这是通过努力将子波化为零相位之后才有的效果。右边是它的子波形态。看到这条剖面使人十分纳闷，除了左边6m砂层有着良好的反映之外，其它的地方是很难理解的！再看分辨率更高的例子。图84 (d)是带通子波1016OHz四个倍频程的、零相位的、理想的高分辨率剖面。我作完这个理论记录后自己都感到吃惊！除了左边浅湖相的6m砂层跃然纸上以外，其它的砂层连影子也

13、不见了。第二节追踪砂层必须用波阻抗剖面或积分地震道看了半天，才明白这张图其实并不是反映着砂层，而是反映着反射系数序列。我将图83 (b)蒙在图84 (d)上，一点不错！(这次你不需要移动，只要对准图框)。几乎每一根反射系数的棍棍都有相应的反映。所以我得到如下的结论：(1)我们反射剖面的波形并不是追踪着砂层，而是追踪着反射系数。(2)反射剖面上的波形胖瘦程度与子波主瓣的胖瘦程度是一致的。子波主瓣瘦的就反映 不了厚砂层；子波胖的反映不了薄砂层。当砂层厚度与子波主瓣一样胖瘦时，振幅就猛的加强，这便是振幅的“调谐”作用。(3)当砂层厚度大于主瓣宽度的时候，砂层的顶、底界分成为两个反射：一个反射为正，另

14、一个反射为负。但由于解释人员不能通过肉眼直观地识别它，所以看不懂它与砂层的联系。不少人拿着108OHz的高分辨率剖面去追踪砂层，以为同相轴中断处便是砂层尖灭，这样的解释是十分危险的。例如图84 (c)中右边，我打X的地方的波峰同相轴就根本不是砂层的尖灭。那么，怎样才能看懂地下的砂层分布情况呢？答案是：你必须把反射波形转化为波阻抗剖面才能理解它。最简单的办法是把反射波形从上到下来一个积分，并去除其直流成分。便称为积分地震道。大家知道反射系数R=f匕匕(70)：2： 2，：T 1当波阻抗反差不大时，P2y2 一 P1V1 = AP",且设 0为P2V2与P1V1的平均值。则有因此对反射系

15、数取积分便近似有1 :可Rdt2"= -ln P、2(71)所以反射系数的积分便正比于波阻抗Pv的自然对数，这是一种最简单的相对波阻抗概O现在我们把图84 (a) (b) (c) (d)都来一次积分，就得到图85 (a) (b*) (c) (d*)四张。先看图85 (d*)带通子波10 160Hz经过积分以后，将它与图 83 (a)的答案相比，简直太好了：每个厚度为 4.5m以上的砂层都反映良好。 3m的砂层也有较弱的反映。（如右边 第43道上）。甚至在左边湖盆中有一层 1.5m的砂层也能追踪（第173道）。如果我们今后五 年左右能够争取频宽达到 10160Hz,这样的剖面会给我们带

16、来极大的用处。箔底君假L事王朝凌忸Egg等奔詹等ir.y三案七减£追初.华阑*£m跄专喜图硅闫会h ,萼魅全壁埠科 -运岳馆底埠石定黑赛理中晶第S南金8t、H替：三H j搭卜但隼我等忡；二时阳为的百性旨酒苜Eg意於案就申F国南澳出H心隆一二一郭运理技，"屋蟹年£日51上比 显露朝逮盥威军一归建倒中露定后屋相去=Mr挡善装冬舜卧T二LTktjL 11 上 "上1-1 * h*_*4_L :，f* *。_认/尸上上上:；1 /*BT - BBH"1TAw<taiK_*AaAa>>aAa )/1Am I*. .JIk-vA

17、yAi .d。. * a k- 2-人. ，' .* - 一号联三sr留把 一.NTdE 一 n>lH-nx(3工)上crg乩LrJC口<LrEHLU.sJ5!管 MT1£rPJ 二上:厂-3y京知谜割金濠总堂上”取决回图85 (c)是10-8OHz的积分地震道。它的分辨率差一些，但是和答案83 (a)相比，凡是6m以上的砂层都有了正确的反映，这是当前我们的技术所允许的高分辨率成果。. * . . - .一 一 、.、 、 、 . . 一 、一图85 (b )是10 4OHz的积分地震道，它相当于目前常规剖面转为波阻抗之后的情况。厚于12m的砂层才能得到正确的反映

18、，只有左边6m的砂层非常突出，可以分辨。在第 15 20道附近，积分道剖面上也还有倾斜的轴，这也是分辨率不够高所引起的“视同相轴倾斜 现象”。图85 (b*)的另外一个现象便是你与答案83 (a)对比时不需要再向上移动图框了，因. 、 一 . . . . . * 一 .、 一 .为经过积分后，砂层的位置都已经还原到它应该在的位置了。图 85 (a)的情况也是如此。这也是积分地震道的好处。 . . * . . . . . . . . . . .从图85( a *)来看，分辨率太低了是没有办法的。即使积分后，也不能很好地分辨出砂层来，但它比不积分的还是稍好些，至少厚砂层的位置对齐了。所以我打趣的向

19、人说： “高分辨率的地震反射叠偏剖面似乎在讲我们听不怪的外国话。 需要一个翻译，这个翻译便是波阻抗或积分地震道”。因而我认为高分辨率勘探的最终成果，交到解释人员手里的剖面不应该是一般的水平叠加或者叠偏剖面，而应该是波阻抗剖面或者是积分地震道剖面。否则解释人员很可能作出错误的解释。“采集处理丁高分辨率地震资料后，如果不作波阻抗或积分道剖面，就好比农夫辛辛苦苦种了庄稼而不去收获！ ”。 我这句话并不太夸大。五年前我访问了美国某家大公司，该公司在墨西哥湾海上寻找非背斜油藏，一直取得成功。其成功决窍有三：亮点、积分地震道和AVO术。他们公司多年来始终规定：凡是送到解释人员手中的剖面必须是积分地震道剖面

20、。我国东部地区早已进入非背斜隐蔽油藏的搜索阶段，所以我积极宣传对砂岩储集层搞积分地震道。当然， 有条件能够作绝对波阻抗当然更好，但是它要花较多的时间和精力。与绝对波阻抗相比，积分道的优点是：( l )递推时累计误差反而小；( 2)计算简单，不需要反射系数的标定； ( 3)没有井的控制也能作。缺点是它不知道波阻抗的绝对值。最近， 胜利油田等单位试作了积分地震道，取得了很好的成效。如果在大庆油田作，效果会更好！第三节 不能丢失低频请再看两个雷克子波的例子。图86 (a)二是主频为40Hz的雷克子波，图86 ( b)是主频为80Hz的雷克子波。右方都画出了其子波波形。从这两张图来看，显然可以看出还是

21、子 波的胖瘦决定了反射剖面的胖瘦。请注意：图 86( a) 子波的胖瘦程度和图84( c) 1O 8OHz子波的主瓣宽度几乎是接近的。图86 (b)的也和图84 (d)的接近。但雷克子波的频宽仅有1.25个二倍频程。(-3dB点的f2/fi只有2.36倍)。所以主频40Hz的图86 (a)的分辨率 要比带通10 80Hz差。(抱歉我把雷克子波的理论记录都画成波形加变面积，而前面图84(a) (b) (c) (d)都画成波谷也涂黑的“密点显示” 了。可能会造成视觉上的差别，不过波形是正确的)。图87 (a*)及图87 (b*)是它们相应的波阻抗(积分地震道)。将它们与图85 (c*)及图85 (

22、d*)比较就显得逊色多了。图 87 (a*)中打“ ？”的地方是与答案图 83 (a)不符之 处。图87 (b*)中打“ ？”的地方是厚砂层变成中央下陷，甚至似乎分为两个砂层的情况。 . . . ._ , . * . _ 这是因为80Hz主频的雷克子波缺乏低频所造成。但图 87 (b)对反映薄砂层还是不错的。再作一个试验，将一条剖面的低频成分滤去，看其结果。我们将图84 (b)的剖面用一个带通滤波器 BP (283570 9OHz对它作一次滤波，于是得到图 88 (a)。滤波后，视 周期的确小了，主频似乎偏高了。但实际上频带更窄，分辨率更差。表现为有连续振动的多相值、其滤波后的实际子波形态画在

23、最右边。图 88 (b)是它的积分地震道。拿它去和图85(b*)相比，效果就差多了：左方6m砂层由一层变为三层；第 12道至第18道似乎地下存在着不少厚砂层，其实不然。右方也出现了许多与答案83 (a)相矛盾之处，我把这种滤去了有效的低频成分而表面上看来视周期变小的现象叫作“假分辨率WLX- V 工-二工三.f 卜三，二 M 专一方 u:=m H=rJxA士石上 tu 曲(C痛忘留之会武理史cL乏箕更叼官,雄4三豪学三笑品区工笑工工雷拿端工二或己至'swe'MT.IEM签$孤我三-嬴三_三号3餐未建京落fi二乾带金祟分.（广）物为40H,甯克子波反射记录的积分地应道 反演近假波

24、阳抗，.反映虱6m的阶层是不 储的！薄眇品没有反映！打-? ”的四处是行砂层分义的而反演曲我没有反映出家匕糠妥三、a克f波, AA、,、A A A-AJk A Akji , Ag4LAAMA卜1 AU A 2LrA卜 R A EPIbEm JUL A A A> Em£.A A A A =4 JhAAAFAAAAA A AAr rA£ArA,A A.AAAr( ArLKAAlAALa Akr一 AlbEA A尸L匚4 r aia 卜 aial a. A.AkAk aa! rAAAA aal 广laaaeg rA a a A-rrrrA a a ALAr £ A a a a rrA > E r Ar, LE-:. trk A>rA、rr A上较为HOH,而克夕波反射氾录的根分施震道（反演近似波阻抗）反映海砂必较好：.对M砂公会造成判断错 «由，甯克f波