大时窗搜索在初至波收获中的应用.docxVIP

大时窗搜索在初至波收获中的应用 0 确定初至波一般准则 目前，地震初至波的自动采集方法主要分为以下几类：一是基于地震记录的多通道相关属性的相关方法；第二种是基于地震记录的即时特性的能量比法和改进算法。第三种是基于地震初至波信息的人工智能方法。第四，参考数据处理方法，将地震记录分类到边缘检测方法，并使用初始值。 在综合考虑以上各类初至波拾取方法优缺点的基础上, 文中选择在天然地震相识别中应用广泛的赤迟信息准则 (AIC) 方法进行地震初至波拾取。考虑到在使用AIC法拾取地震初至波过程中首先需要确定在给定的计算时窗内存在地震初至波, 以往研究者通过能量比法得到该时窗非常费时, 尤其是当地震数据体比较大时, 逐道搜索将会导致自动拾取速度很慢, 故结合炮检波点分布关系及初至波走时关系, 在约束炮初至波精确拾取的基础上, 估算预测出其他炮的初至波所在时窗, 在该时窗内使用AIC方法精确拾取地震初至波。提出基于空变时窗约束的地震初至波自动拾取方法 1 时窗约束方法 地震初至波自动拾取思路为:在约束炮初至波精确拾取的基础上, 使用空变时窗约束方法预估其他炮地震初至波所在时窗, 然后在该时窗内使用赤迟信息准则法确定初至时间, 再根据初至时间与偏移距分布关系剔除拾取误差较大的初至时间, 最后对每个初至时间进行微调得到初至波波峰对应的时刻。 1.1 初至波预测时间t 根据野外采集工作的认识, 地理位置相近放炮, 接收的单炮地震资料上的初至波形态基本一致。故分析观测系统, 选定可以对其他炮都有约束效果的炮作为约束炮, 根据观测系统信息, 对所有炮按照距离约束炮空间位置由近及远的顺序排序, 在相邻的两炮上相同检波点上, 上一炮的初至波形态与下一炮的初至波形态基本一致 考虑到当野外放炮顺序与炮点空间位置排序顺序差异较大时, 地震记录与炮点空间位置排序结果顺序不一致。若按以上方法进行计算时需要对地震记录进行排序, 当地震记录数据量大时, 使用以上约束方法反而容易出错, 故课题中引入表示地震初至波趋势的初至波预测时间。 初至波中虽然包含直达波、折射波、反射波、转换波等, 但目前地震勘探中初至波主要以折射波和直达波为主 式中:t 所以, 地震初至波的预测时间应该包括作为约束炮上相应道约束道的初至时间和被约束道与约束道的时差, 约束道初至时间已知, 故预测时间计算即可转化为相邻道直达波或者折射波时差计算。起伏地表层状介质折射初至射线路径如图1所示。 如图1所示, 对于分别在相邻的2个炮点S 折射波时差可表示为: 根据三角关系, 式 (3) 也可表示为: 当地表近水平, 反射层为水平层状介质时, O 而当地表及地下构造复杂时, 炮点空间距离与O 1.2 ar-aic取值法 在地震记录中, 由于噪声信号和初至波信号具有不同的统计特性, 可分别看作是一个平稳过程, 在二者分界处, 噪声信号和初至波信号统计特性差别较大, 在最小平方意义下这2种信号的拟合度最差, 因此该分界位置可解释为相位起跳位置即为地震初至波, 这个计算过程被定义为AR-AIC拾取法 式中, var表示数据序列的方差;N表示每道采样点数;k表示逐个取计算时窗内的样点。 在1.1小节中得到的预测时窗内依次从第一个点开始, 利用公式 (5) 可直接计算出每个采样点的AIC值。在初至预测时窗内AIC值最小的点对应的时刻即为精确拾取的地震初至波。 1.3 时间-偏移距随机分布特征 地震初至波按时间-偏移距排列时呈现线性分布, 拾取正确的初至波时间-偏移距会分布在一条或者多条线性束上, 而当初至波拾取错误时, 该初至波在时间-偏移距分布图上会呈现随机分布的特征。故可根据该特征将偏离线性束较远的点直接去掉, 不参与后续的静校正计算。 1.4 采样点的确定 以上自动拾取到的地震初至波是起跳位置, 而在资料处理过程中多拾取的是初至波的波峰, 可使用希尔伯特变换, 根据模型约束后的初至时间计算地震初至波的波峰位置 对于第i道地震数据, 根据前一步拾取的初至时间, 找到初至波前后的M个采样点, 其中 式中N表示单道地震数据的采样点数。 根据以下公式求这些采样点的希尔伯特变换: 式中:h (t) 为x (t) 对应的希尔伯特变换;x (t) 是t时刻地震道的振幅值。 然后寻找满足公式 (8) 的点, 即为波峰N 式中:cosθ (t) 为瞬时相位函数;a (t) 为瞬时包络。 根据以上算法, 计算每道上地震初至波的波峰。 2 空变时窗约束正确处理初至时间 采用某地二维活断层地震勘探数据对文中所提方法的有效性进行测试。该二维活断层资料共有134炮, 最小文件号为63, 测线地表高差起伏较大, 高差达450 m。 首先选定文件号63的单炮作为约束炮, 对SPS文件及约束炮信息进行综合分析, 最终按照空间距离信息对所有炮点进