一种地震相干层析成像方法技术

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

An seismic coherence tomography method

The invention discloses a method of seismic coherence tomography, which belongs to the technical field of seismic exploration, the method includes using dynamic time warping method based on time window envelope coherence, acquiring seismic gather events coherent transit time; according to the acquired events coherent time difference, combined with the iterative method of nonlinear tomography ray path difference. Get the degree of access to the slow field; slowness field reciprocal, generating velocity field. The invention utilizes the travel time information of the reflected wave to estimate the depth of the velocity field, and obtains a high precision velocity field by combining the high-precision coherent time difference acquisition method.

【技术实现步骤摘要】
一种地震相干层析成像方法
本专利技术涉及地震勘探
，尤其涉及一种地震相干层析成像方法。
技术介绍
准确获取地下介质的速度参数是地震勘探的核心问题之一，为勘探开发深层超深层油气，必须获得深层地质介质的精细速度场。现有使用天然地震初至走时层析成像反演获得的速度场精细度不够，精度不高；使用人工反射地震走时数据测算地下介质速度，需要有测井或VSP资料约束标定，无法准确获得井未钻达层位或无井地区的介质速度；使用人工反射地震的层析成像方法，因巨大的计算需求和高度的非线性，无法适用于深层条件。因此，以上获取深层地质介质的精细速度场的方式均无法满足深层超深层油气勘探需求。
技术实现思路
为了解决现有技术的问题，本专利技术实施例提供了一种地震相干层析成像方法。所述技术方案如下：一种地震相干层析成像方法，包括：利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差；根据获取到的地震炮道集中的同相轴相干时差，结合射线路径差的非线性层析迭代方法，获取慢度场；对获取到的慢度场求倒数，生成速度场。进一步的，利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差，包括：在本文档来自技高网...

【技术保护点】
一种地震相干层析成像方法，其特征在于，包括：利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差；根据获取到的地震炮道集中的同相轴相干时差，结合射线路径差的非线性层析迭代方法，获取慢度场；对获取到的慢度场求倒数，生成速度场。

【技术特征摘要】
1.一种地震相干层析成像方法，其特征在于，包括：利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差；根据获取到的地震炮道集中的同相轴相干时差，结合射线路径差的非线性层析迭代方法，获取慢度场；对获取到的慢度场求倒数，生成速度场。2.根据权利要求1所述的地震相干层析成像方法，其特征在于，利用基于时窗包络相干性的动态时间规整方法，获取地震炮道集中的同相轴相干时差，包括：在第is个地震炮集中，第ig个检波器和其相邻的检波器的地震记录分别为dis,ig(it)和dis,ig+1(it)；其中，it为时间采样点号，it取整数，it范围为1≤it≤lt，lt为地震道总的时间采样点；之后，进行Hilbert变换，得到，其中，dh1和dh2分别为dis,ig和dis,ig+1的Hilbert变换结果，iτ为计算过程中的临时变量，Lh为Hilbert算子长度，h为Hilbert算子，h为：包络信号e1(it)和e2(it)分别为：为了弱化振幅因素对相干时差提取的影响，对包络信号e1(it)和e2(it)进行如下处理，e1'(it)＝e1(it)/ae2'(it)＝e2(it)/b其中，e1'(it)和e2'(it)为e1(it)和e2(it)处理后的结果，a和b分别为e1(it)和e2(it)时间序列上的振幅绝对值的最大值；计算e1'(it)和e2'(it)在局部时间窗范围内的对齐误差D(it1,it2)，其中，it1和it2分别为e1'(it)和e2'(it)在计算过程中的时间采样点号，it1和it2为整数，取值范围为1≤it1≤lt，1≤it2≤lt；Lw为计算对齐误差时的半时窗长度；利用动态时间规整方法中的二阶对称公式，计算递归累积误差和A(it1,it2)，A(it1,it2)＝D(it1,it2)+g其中，第一个地震道的it1时刻与第二个地震道的it2时刻的相干时差为Δt(it1)，Δtis,ig(it1)＝(it2-it1)*dt其中，dt为地震道采样时间间隔；在每一行上搜索相干性最高的元素时需要遵循以下条件：①当某行上存在多个最小值时，取最靠近对角线的元素；②相邻两个时间点的相干时差需满足：it1*dt+Δtis,ig(it1)≤(it1+1)*dt+Δtis,ig(it1+1)；在同一炮集中，将上述方法应用于每一组相邻地震道，即可获得每一个地震道的相干时差；另外，在后续的层析成像步骤中，应用的相干时差通常是某一地震同相轴的相干时差；因此，需要沿同相轴提取相干时差；在第is个炮集中的第ig个地震道上给定一个同相轴的种子点it3，其对应的相干时差为Δtis,ig(it3)，那么其临近的下一地震道ig+1的同相轴点it4为，it4＝it3+Δtis,ig(it3)/dt这一地震道的同相轴点处的相干时差为Δtis,ig+1(it4)；同样地，可以用相同的方式追踪下一道的同相轴点；在该炮集中的所有地震道上进行追踪，即可获得该同相轴对应的所有相干时差。3.据权利要求2...

【专利技术属性】
技术研发人员：曹俊兴，胡江涛，何晓燕，蒋旭东，廖小芳，吕雪松，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川,51