[发明专利]一种地震相干层析成像方法

说明书

本发明公开了一种地震相干层析成像方法，属于地震勘探技术领域，所述方法包括利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差；根据获取到的同相轴相干时差，结合射线路径差的非线性层析迭代方法，获取慢度场；对获取到的慢度场求倒数，生成速度场。本发明利用反射波的走时信息实现深层速度场的估计，结合高精度的相干时差获取方法，获得高精度的速度场。

技术领域

本发明涉及地震勘探技术领域，尤其涉及一种地震相干层析成像方法。

背景技术

准确获取地下介质的速度参数是地震勘探的核心问题之一，为勘探开发深层超深层油气，必须获得深层地质介质的精细速度场。

现有使用天然地震初至走时层析成像反演获得的速度场精细度不够，精度不高；使用人工反射地震走时数据测算地下介质速度，需要有测井或VSP资料约束标定，无法准确获得井未钻达层位或无井地区的介质速度；使用人工反射地震的层析成像方法，因巨大的计算需求和高度的非线性，无法适用于深层条件。

因此，以上获取深层地质介质的精细速度场的方式均无法满足深层超深层油气勘探需求。

发明内容

为了解决现有技术的问题，本发明实施例提供了一种地震相干层析成像方法。所述技术方案如下：

一种地震相干层析成像方法，包括：

利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差；

根据获取到的地震炮道集中的同相轴相干时差，结合射线路径差的非线性层析迭代方法，获取慢度场；

对获取到的慢度场求倒数，生成速度场。

进一步的，利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差，包括：

在第is个地震炮集中，第ig个检波器和其相邻的检波器的地震记录分别为d_is,ig(it)和d_is,ig+1(it)；

其中，it为时间采样点号，it取整数，it范围为1≤it≤lt，lt为地震道总的时间采样点；

之后，进行Hilbert变换，得到，

其中，d_h1和d_h2分别为d_is,ig和d_is,ig+1的Hilbert变换结果，iτ为计算过程中的临时变量，Lh为Hilbert算子长度，h为Hilbert算子，h为：

包络信号e₁(it)和e₂(it)分别为：

进一步的，利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差，还包括：

为了弱化振幅因素对相干时差提取的影响，对包络信号e₁(it)和e₂(it)进行如下处理，

e₁'(it)＝e₁(it)/a

e₂'(it)＝e₂(it)/b

其中，e₁'(it)和e₂'(it)为e₁(it)和e₂(it)处理后的结果，a和b分别为e₁(it)和e₂(it)时间序列上的振幅绝对值的最大值。