三维地震数据动态拉伸时差校正方法及系统技术方案

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

Dynamic Tensile Time Difference Correction Method and System for 3-D Seismic Data

The invention discloses a method and system for dynamic stretching time difference correction of three-dimensional seismic data, which includes: 1) reading the first group of three-dimensional seismic data and the second group of three-dimensional seismic data line by line, calculating the time shift of each point of the two groups of three-dimensional data, obtaining three-dimensional time shift data; 2) smoothing the three-dimensional time shift data; 3) using the smoothed time shift data for processing. The second group of three-dimensional seismic data is corrected by time difference line by line. By applying the dynamic stretching time difference correction method of the three-dimensional seismic data proposed by the invention, the accuracy of the two-dimensional time difference correction is improved through energy consistency processing. By line-by-line processing, the two-dimensional time difference correction method is extended to three-dimensional, so as to better serve the practical application.

【技术实现步骤摘要】
三维地震数据动态拉伸时差校正方法及系统
本专利技术属于地震数据处理领域，更具体地，涉及一种三维地震数据动态拉伸时差校正方法及系统。
技术介绍
地震数据处理中有时需要对不同数据进行时差(深度差做深时转换后可转换为时差)校正，比如模拟数据与观测数据的匹配、不同偏移距(或传播角度)成像剖面的对齐、P波成像剖面与S波成像剖面的标定和剩余时差校正等。当数据间的时移误差比较大，或时移误差随着空间或时间变化时，时差校正是比较困难的。美国科罗拉多矿院DaveHale提出了一种动态时间拉伸(dynamictimewarping，DTW)方法，可以有效对一维、二维数据进行动态拉伸，在PP波剖面与PS波剖面标定对齐方面得到较好的应用。专利技术人发现，如果两个数据能量关系存在差异，直接利用上面提到的二维方法无法得到好的时移校正量，影响时差校正的效果，另外在实际应用中，我们往往需要对三维数据进行时差校正，因此有必要提出一种适用于三维地震数据动态拉伸时差校正方法。
技术实现思路
本专利技术面向三维地震数据时差校正需求，针对二维动态时间拉伸存在的问题，提出一种三维地震数据动态拉伸时差校正方法，目的是形成一种处理效果更好的时差校正方法，用于成像时差校正和剖面对齐。根据本专利技术的一方面，提出了一种三维地震数据动态拉伸时差校正方法，该方法包括：1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；2)对所述三维时移量数据进行平滑处理；3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，其中，步骤3)包括：3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；3.2)对于所述第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正：3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。优选地，所述第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个所述维度的样点数分别表示为ny、nx和nz。优选地，步骤1)包括：1.1)读取第二组三维地震数据其中一条线数据D2(nx,nz)和相应的第一组三维地震数据的线数据D1(nx,nz)；1.2)计算所述第二组三维地震数据的线数据D2(nx,nz)相对于所述第一组三维地震数据的线数据D1(nx,nz)的每个点的时移量，将所述时移量按线号位置存放在三维数组中；1.3)读取第二组三维地震数据和第一组三维地震数据的下一条线数据，继续进行时移量计算，直到第二组三维地震数据的所有线数据完成时移量计算，获得所述三维时移量数据T(ny,nx,nz)。优选地，在步骤1)中，通过以下步骤计算时移量：a)设置窗口长度为nw，滑动步长为nl的时窗；b)在所述设置时窗内，利用第一组三维地震数据的线数据，对第二组三维地震数据的线数据进行能量一致性处理；c)计算同一时窗内经能量一致性处理后的第二组三维地震数据的线数据相对于所述第一组三维地震数据的线数据的时移量；d)所有时窗的时移量计算结束后，求取不同时窗计算的时移量的平均值，作为时窗重叠部分的时移量。优选地，对时窗内的第二组三维地震数据的线数据进行能量一致性处理的公式为：其中，j,k为二维数据中一点的序号，W2(j,k)为在时窗内第二组三维地震数据的线数据，E1为第一组三维地震数据的时窗内线数据的能量，E2为第二组三维地震数据的时窗内线数据的能量。优选地，平滑处理后的三维时移量数据表示为：其中，i,j,k为三维数据中一点的序号，M是每个方向的平滑范围的半径，若i+n＜1，则令(i+n)的值为1，若i+n＞ny，则令(i+n)的值为ny，若j+m＜1，则令(j+m)的值为1，若j+m＞nx，则令(j+m)的值为ny，若k+p＜1，则令(k+p)的值为1，若k+p＞nz，则令(k+p)的值为ny。优选地，步骤3.2)中的时差校正后的第二组三维地震数据的线数据表示为：D'2(j,k)＝D2(j,k1)(3)k1＝k+Ts(j,k)(4)其中，j,k为二维数据中一点的序号，Ts(j,k)为经平滑处理后的三维时移量数据，若k1＜1，则令k1的值为1，若k1＞nz，则令k1的值为nz。根据本专利技术的另一方面，提出了一种三维地震数据动态拉伸时差校正系统，所述系统包括：存储器，存储有计算机可执行指令；处理器，运行所述存储器上的计算可执行指令时，所述处理器实现以下步骤：1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；2)对所述三维时移量数据进行平滑处理；3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，其中，步骤3)包括：3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；3.2)对于所述第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正：3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。优选地，所述第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个所述维度的样点数分别表示为ny、nx和nz。优选地，步骤1)包括：1.1)读取第二组三维地震数据其中一条线数据D2(nx,nz)和相应的第一组三维地震数据的线数据D1(nx,nz)；1.2)计算所述第二组三维地震数据的线数据D2(nx,nz)相对于所述第一组三维地震数据的线数据D1(nx,nz)的每个点的时移量，将所述时移量按线号位置存放在三维数组中；1.3)读取第二组三维地震数据和第一组三维地震数据的下一条线数据，继续进行时移量计算，直到第二组三维地震数据的所有线数据完成时移量计算，获得所述三维时移量数据T(ny,nx,nz)。本专利技术的有益效果在于：面向三维地震数据时差校正需求，针对两个数据能量关系不一致时已有二维动态拉伸方法存在的问题，提供了一种更适用于实际数据处理的三维地震数据动态拉伸时差校正方法，对于能量关系差异较大的两个三维地震数据先进行能量一致性处理，提高了二维时移量时差校正的准确度，从而获得较高精度的动态时移量，进而将两个数据校正为一致，通过逐线处理，将二维时移量时差校正方法扩展到三维，使之更好的服务于实际应用之中。本专利技术的其它特征和优点将在随后具体实施方式部分予以详细说明。附图说明通过结合附图对本专利技术示例性实施方式进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的三维地震数据动态拉伸时差校正方法的流程图。图2示出了根据本专利技术的一个实施例的偏移方法1的成像剖面。图3示出了根据本专利技术的一个实施例的偏移方法2的成像剖面。图4示出了根据本专利技术的一个实施例的时移量剖面。图5示出了根据本专利技术的一个实施例的偏移方法2本文档来自技高网...

【技术保护点】
1.一种三维地震数据动态拉伸时差校正方法，其特征在于，该方法包括：1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；2)对所述三维时移量数据进行平滑处理；3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，其中，步骤3)包括：3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；3.2)对于所述第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正；3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。

【技术特征摘要】
1.一种三维地震数据动态拉伸时差校正方法，其特征在于，该方法包括：1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；2)对所述三维时移量数据进行平滑处理；3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，其中，步骤3)包括：3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；3.2)对于所述第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正；3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。2.根据权利要求1所述的三维地震数据动态拉伸时差校正方法，其中，所述第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个所述维度的样点数分别表示为ny、nx和nz。3.根据权利要求1所述的三维地震数据动态拉伸时差校正方法，其中，步骤1)包括：1.1)读取第二组三维地震数据其中一条线数据D2(nx,nz)和相应的第一组三维地震数据的线数据D1(nx,nz)；1.2)计算所述第二组三维地震数据的线数据D2(nx,nz)相对于所述第一组三维地震数据的线数据D1(nx,nz)的每个点的时移量，将所述时移量按线号位置存放在三维数组中；1.3)读取第二组三维地震数据和第一组三维地震数据的下一条线数据，继续进行时移量计算，直到第二组三维地震数据的所有线数据完成时移量计算，获得所述三维时移量数据T(ny,nx,nz)。4.根据权利要求3所述的三维地震数据动态拉伸时差校正方法，其中，在步骤1)中，通过以下步骤计算时移量：a)设置窗口长度为nw，滑动步长为nl的时窗；b)在所述设置时窗内，利用第一组三维地震数据的线数据，对第二组三维地震数据的线数据进行能量一致性处理；c)计算同一时窗内经能量一致性处理后的第二组三维地震数据的线数据相对于所述第一组三维地震数据的线数据的时移量；d)所有时窗的时移量计算结束后，求取不同时窗计算的时移量的平均值，作为时窗重叠部分的时移量。5.根据权利要求4所述的三维地震数据动态拉伸时差校正方法，其中，对时窗内的第二组三维地震数据的线数据进行能量一致性处理的公式为：其中，j,k为二维数据中一点的序号，W2(j,k)为在时窗内第二组三维地震数据的线数据，E1为第一组三维地震数据的时窗内线数据的能量，E2为第二组三维地震数据的时窗内线数据的能量。6.根据权利要求1所述的三维地震数据动态拉伸时...

【专利技术属性】
技术研发人员：段心标，刘定进，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11