时频图上仅反映出这个等效层的顶底反射界面。
V泥=3200m/s ρ泥=2323.09kg/m3 砂泥岩互层厚度=10米 V砂=3500m/s ρ砂=2375.72kg/m3
图3.41厚度为10m的砂泥岩等厚度互层模型及其反射界面的脉冲响应以及合成记录
图3.42 合成地震记录的时频分析结果(砂泥岩互层厚度都为10m)
图3.43为砂泥岩薄互层厚度为5m的六层砂泥岩等厚度互层模型及其反射界面的脉冲响应以及合成记录。图3.44对应合成地震记录的四种时频谱分析结果对比图。与图3.42的结果类似，双谱和三谱时频图上虽然反映出5m厚的砂体顶底反射界面，但对应的时间与实际砂体顶底反射界面的时间不同，说明，在砂体厚度比较薄时，上述三种时频分析方法也不能准确地分辨出砂体顶底反射界面。也进一步说明，一系列厚度为5m的砂泥岩薄互层等效于砂泥总厚度的单个层。合成记录以及时
频图上仅反映出这个等效层的顶底反射界面。
ρ泥=2323.09kg/m 3V泥=3200m/s 砂泥岩互层厚度=5米 V砂=3500m/s ρ砂=2375.72kg/m3
图3.43厚度为5m的砂泥岩等厚度互层模型及其反射界面的脉冲响应以及合成记录
图3.44 合成地震记录的时频分析结果(砂泥岩互层厚度都为5m)
3.4.6.3 楔形弹性介质模型
这里我们设计了一个楔形弹性介质模型，模型结构及弹性参数如图3.45(a)所示，对应的合成记录如图3.45(b)所示。计算时，雷克子波的主频为30Hz。合成记录上，34道以后的地震记录上，已经难以区分砂体的顶底反射界面。
V泥=3200m/s ρ泥=2323.1kg/m 3V砂=3500m/s ρ砂=2375.7kg/m 3砂泥岩厚度为从100米变化到0米，砂岩上覆泥岩厚度为100米，道间距为20米，无吸收衰减。 (a)
(b)
图3.45 楔状模型结构及对应的合成记录
图3.46为对应合成地震记录的频率分别为10Hz、20Hz、30Hz、40Hz、50Hz的双谱时频分析结果对比图，图中虚线为楔状砂体顶底界面反射时间。比较这些图可以看出，双谱时频分析结果在时间和频率上的分辨率能同时达到较高的分辨率。在频率为30Hz的双谱时频图上，CDP道号为50时，时频图仍然能较准确地分辨出砂体顶底反射界面，此时，砂体厚度大约1/8的地震波长，因此，双谱能够提高砂体厚度的分辨能力，这也可以从图3.47显示的各道时频图上看出来。
图3.46 对应合成地震记录的不同频率对应的双谱时频分析结果对比图
图3.47不同砂体厚度对应的合成地震道的时频图
3.5 实际地震记录的高阶谱计算
3.5.1 桩海地区某地震剖面
桩海地区位于胜利油区东北角，地质构造上属于埕岛—垦东断裂构造带的北段。
为了验证高阶谱对实际地震记录的含油气层位分辨能力，我们选取了该区一条过井地震剖面，如图3.48～3.49所示，产油层为Tg2层。我们采用40ms时窗进行双谱时频分析，图3.50显示出由双谱最大能量表示的地震记录剖面，图3.51显示出由频率为10Hz对应的双谱能量表示的地震记录剖面，图3.52显示出由频率为20Hz对应的双谱能量表示的地震记录剖面。图3.53～图3.56分别为不同地震道目的层段的双谱分析。图3.57～图3.59分别为20ms、40ms和60ms时窗的井旁地震道的双谱时频分布图。
埕北305 埕北302 埕北303 Tg1 Tg2
图3.48一条过井实际地震剖面
图3.49过井实际地震记录的彩色显示