高密度地震的资料品质综合分析技术的研究.docVIP

高密度地震的资料品质综合分析技术的研究 　　 【摘 要】在地震勘探中，对采集数据的定量分析是评价采集质量和制定处理流程、参数选择的依据，也是分析资料能否满足勘探需要的重要手段。与原始地震资料品质有关的属性包括有效信号能量、频率、噪音能量、频率、资料信噪比等方面。常规采集中，质量监控的手段主要有对记录的分析和现场处理剖面监控。当面对超万道的高密度采集资料时，常规的监控分析技术不论在技术方面还是时效方面均不能适应。在这种情况下，展开高密度采集质量监控分析技术的研究就很有必要。 　　【关键词】能量;频率;信噪比;定量分析 　　0 引言 　　随着地震勘探技术的不断深入，野外采集地震道数越来越多，特别是近几年来万道地震仪的广泛应用，使地震数据量成倍增长，地震处理（数据输入、定义观测系统、解编记录、滤波等）工作量越来越大。按照常规的处理手段不论在技术方面还是时效方面均已不能满足处理要求，不能及时有效的进行反馈，严重阻碍野外生产质量、生产进度，面对高密度地震数据在正式确定处理流程之前需要定性、定量地对地震资料进行评价，对数据中的激发能量、激发子波、高频干扰和信噪比等参数进行定量、系统、直观地分析为科学地制定后续处理流程以及快速、有效地实现对处理数据的全面质量控制提供依据。 　　1 综合能量分析及应用 　　实际地震勘探中激发能量受各种因素的影响而发生变化，传播路径、储层物性、地表岩性、激发井深和爆炸药量等因素均会对激发能量造成影响。激发能量的变化又将直接影响激发频率和激发子波的变化。因此在地震资料处理前有效地分析各炮的激发能量变化，依据近地表资料和静校正分析结果为后续地震资料处理提供炮集质量的变化和近地表变化的信息。该技术为快速分析各炮激发能量提供了手段。 　　1.1 方法原理 　　利用希尔伯特变换能有效地从地震信号中提取复杂信号的瞬时参数一一瞬时振幅、瞬时频率、瞬时相位。设一个解析信号可表示为依赖于时间变化的复变量，即 　　Ut=Xt+jYt（1） 　　式中：X（t）为信号本身;Y（t）为它的正交。正交是记录信号的900相移。对X（t）进行希尔伯特反变换就可以得到 　　Yt=■Xt（2） 　　代入式（1）得 　　Ut=Xt+j■Xt=δt+j■Xt（3） 　　U（t）的指数表达式为 　　Ut=Rte■（4） 　　式中：R（t）为瞬时振幅;？?Jt为瞬时相位。 　　希尔伯特瞬时振幅变换能使地震剖面上能量包络更加清晰，该技术利用希尔伯特变换得到的振幅包络对属性参数进行统计分析的功能。在统计激发能量时，将分析时窗分成7个子时窗分别统计，最后求取平均值。 　　1.2 能量分析技术应用 　　1.2.1 野外激发能量监控 　　从罗家全区能量统计平面图可以清楚地看出全工区各炮反射能量的变化。全区中部和南部初至能量强，即激发能量强，通过叠加后能量分布比较均衡。根据分析结果，处理人员不仅可以全面了解原始资料的品质，而且可以通过平面监测与监控炮相结合的方式全程监控整个处理过程。 　　1.2.2 干扰能量分析 　　地震采集中经常遇到的干扰波主要分为规则干扰和环境噪音两类。噪音是影响资料质量的主要因素，直接影响了剖面的成像效果。根据干扰能量分析图，可以快捷有效的分析全区或整束线的噪音分布情况，查找干扰源，排除干扰，提高资料质量。同时，还可以根据噪音的能量分布特征来压制噪音，提高信噪比。 　　1.2.3 振幅处理效果监控 　　在资料处理中，振幅一致性处理是一个重要环节，振幅处理效果如何，如果逐炮分析，既费时，又不直观。借助综合能量分析的方法，可以方便直观的进行处理效果监控。 　　2 信噪比综合定量分析 　　剖面法信噪比估算方法基于地震数据存在相干性的原理。它以小时窗作为研究对象，并滑动外推进行整条测线的信噪比计算，抛开常规信噪比估算用一个时窗内信号和噪音总能量的计算方法，因此对全测线资料信噪比的描述更加直观。它的估算前提有两个，一个是信号的相关性，二个是噪音的随机分布。剖面法信噪比估算理论上基本可以分为三个过程： 　　2.1 地震信号的相关性判断及能量估算 　　假设相邻两道中地震波分别用f1（t）和f2（t）表示，计算两道的互相关就可表为： 　　Φτ=■f■tf■tdt 　　Φ=■（5） 　　两道相关性决定于相关系数： 　　根据？?J值可以进行信号或噪音的判断。给定一个门槛值，假设？?J的值大于这个门槛值，则认为地震道相干，作为有效信号统计;假设？?J的值小于这个门槛值，则认为地震道不相干，作为噪音统计。 　　在信噪比分析时窗内，将数据体记为： 　　D=d■M×N（6） 　　式中，M为时窗内采样点数，N为时窗内地震道数，i=1，2，3，…，M，j=1，2，3，…，N。 　　假设在整个时窗内地震