[发明专利]三维地震数据动态拉伸时差校正方法及系统

说明书

本发明公开了一种三维地震数据动态拉伸时差校正方法及系统，包括：1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；2)对三维时移量数据进行平滑处理；3)利用经平滑处理的时移量数据对于第二组三维地震数据逐线进行时差校正。应用本发明提出的三维地震数据动态拉伸时差校正方法，通过能量一致性处理，提高了二维时移量时差校正的准确度，通过逐线处理，将二维时移量时差校正方法扩展到三维，使之更好的服务于实际应用之中。

技术领域

本发明属于地震数据处理领域，更具体地，涉及一种三维地震数据动态拉伸时差校正方法及系统。

背景技术

地震数据处理中有时需要对不同数据进行时差(深度差做深时转换后可转换为时差)校正，比如模拟数据与观测数据的匹配、不同偏移距(或传播角度)成像剖面的对齐、P波成像剖面与S波成像剖面的标定和剩余时差校正等。当数据间的时移误差比较大，或时移误差随着空间或时间变化时，时差校正是比较困难的。美国科罗拉多矿院Dave Hale提出了一种动态时间拉伸(dynamic time warping，DTW)方法，可以有效对一维、二维数据进行动态拉伸，在PP波剖面与PS波剖面标定对齐方面得到较好的应用。

发明人发现，如果两个数据能量关系存在差异，直接利用上面提到的二维方法无法得到好的时移校正量，影响时差校正的效果，另外在实际应用中，我们往往需要对三维数据进行时差校正，因此有必要提出一种适用于三维地震数据动态拉伸时差校正方法。

发明内容

本发明面向三维地震数据时差校正需求，针对二维动态时间拉伸存在的问题，提出一种三维地震数据动态拉伸时差校正方法，目的是形成一种处理效果更好的时差校正方法，用于成像时差校正和剖面对齐。

根据本发明的一方面，提出了一种三维地震数据动态拉伸时差校正方法，该方法包括：

1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；

2)对所述三维时移量数据进行平滑处理；

3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，

其中，步骤3)包括：

3.1)读取第二组三维地震数据的其中一条线数据D₂(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据T_s(nx,nz)；

3.2)对于所述第二组三维地震数据的线数据D₂(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量T_s(nx,nz)进行时差校正：

3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。

优选地，所述第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个所述维度的样点数分别表示为ny、nx和nz。

优选地，步骤1)包括：

1.1)读取第二组三维地震数据其中一条线数据D₂(nx,nz)和相应的第一组三维地震数据的线数据D₁(nx,nz)；

1.2)计算所述第二组三维地震数据的线数据D₂(nx,nz)相对于所述第一组三维地震数据的线数据D₁(nx,nz)的每个点的时移量，将所述时移量按线号位置存放在三维数组中；

1.3)读取第二组三维地震数据和第一组三维地震数据的下一条线数据，继续进行时移量计算，直到第二组三维地震数据的所有线数据完成时移量计算，获得所述三维时移量数据T(ny,nx,nz)。