CN109655912A - 三维地震数据动态拉伸时差校正方法及系统 - Google Patents

三维地震数据动态拉伸时差校正方法及系统 Download PDF

Info

Publication number
CN109655912A
CN109655912A CN201710935821.8A CN201710935821A CN109655912A CN 109655912 A CN109655912 A CN 109655912A CN 201710935821 A CN201710935821 A CN 201710935821A CN 109655912 A CN109655912 A CN 109655912A
Authority
CN
China
Prior art keywords
seismic data
group
data
shift amount
time shift
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710935821.8A
Other languages
English (en)
Other versions
CN109655912B (zh
Inventor
段心标
刘定进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710935821.8A priority Critical patent/CN109655912B/zh
Publication of CN109655912A publication Critical patent/CN109655912A/zh
Application granted granted Critical
Publication of CN109655912B publication Critical patent/CN109655912B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

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

Description

三维地震数据动态拉伸时差校正方法及系统
技术领域
本发明属于地震数据处理领域，更具体地，涉及一种三维地震数据动态拉伸时差校正方法及系统。
背景技术
地震数据处理中有时需要对不同数据进行时差(深度差做深时转换后可转换为时差)校正，比如模拟数据与观测数据的匹配、不同偏移距(或传播角度)成像剖面的对齐、P波成像剖面与S波成像剖面的标定和剩余时差校正等。当数据间的时移误差比较大，或时移误差随着空间或时间变化时，时差校正是比较困难的。美国科罗拉多矿院Dave Hale提出了一种动态时间拉伸(dynamic time warping，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经过时差校正后的剖面。
图6示出了根据本发明的一个实施例的三维时移量数据的一个时间切片。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例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详细说明根据本发明的三维地震数据动态拉伸时差校正方法的具体步骤。
步骤1，分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据。
在一个示例中，第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个维度的样点数分别表示为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)所有时窗的时移量计算结束后，求取不同时窗计算的时移量的平均值，作为时窗重叠部分的时移量。
具体地，两个二维数据能量关系不一致时，需要进行分时窗能量一致性处理，然后再计算时移量。
首先，分时窗设置两个参数，一是窗口长度nw，二是滑动步长nl，且窗口长度减滑动步长为重叠部分长度。
其次，对于一个时窗内的两个地震数据W1(nx,nw)和W2(nx,nw)而言，分别统计所有点的平方和，分别记为E1和E2，然后对第二个地震数据进行能量处理。
在一个示例中，对时窗内的第二组三维地震数据的线数据进行能量一致性处理的公式为：
其中，j,k为二维数据中一点的序号，W2(j,k)为在时窗内第二组三维地震数据的线数据，E1为第一组三维地震数据的时窗内线数据的能量，E2为第二组三维地震数据的时窗内线数据的能量。
然后，能量一致性处理后，按照背景技术中提到的动态时间拉伸(dynamic timewarping，DTW)方法计算两个二维数据的时移量。
最后，所有时窗的时移量计算结束后，对时窗重叠部分，求取不同时窗计算的时移量的平均值，作为重叠部分的时移量。这样即可得到两个二维数据剖面间的时移量。
步骤2，对三维时移量数据进行平滑处理。
在一个示例中，对三维时移量数据的每一点而言，计算周围一定范围内所有点时移量的平均值，作为该点的时移量，平滑处理后的三维时移量数据表示为：
其中，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，利用经平滑处理的时移量数据对于第二组三维地震数据进行时差校正，其中，步骤3)包括：
3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；
3.2)对于第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正；
3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。
在一个示例中，步骤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。
本实施例面向三维地震数据时差校正需求，针对两个数据能量关系不一致时已有二维动态拉伸方法存在的问题，提供了一种更适用于实际数据处理的三维地震数据动态拉伸时差校正方法，对于能量关系差异较大的两个三维地震数据先进行能量一致性处理，提高了二维时移量时差校正的准确度，从而获得较高精度的动态时移量，进而将两个数据校正为一致，通过逐线处理，将二维时移量时差校正方法扩展到三维，使之更好的服务于实际应用之中。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
某区分别用两种不同偏移方法得到两个偏移成像结果，两者存在细微的时差。为了便于对比分析，对两个三维数据进行时差校正。
图2示出了在上述应用示例中的偏移方法1的成像剖面。图3示出了在上述应用示例中的偏移方法2的成像剖面。图4示出了在上述应用示例中的时移量剖面。图5示出了在上述应用示例中的偏移方法2经过时差校正后的剖面。
给出两组三维地震数据，图2-图3是某一条线的两个不同的成像剖面，图2是第一组三维地震数据的其中一条线数据，图3是与之相应的第二组三维地震数据的线数据，对比可以看出它们之间的时差。
首先计算图3所示的线数据相对于图2线数据的每个点的时移量，计算出的时移量如图4所示，逐线处理获得第二组三维时移量数据；然后，第二组三维时移量数据进行平滑处理；最后，利用经平滑处理的时移量数据对于第二组三维地震数据进行时差校正。图5是图3所示的线数据经过动态拉伸时差校正后的剖面，可以看出它与图2所示的线基本不再存在时差问题。图6给出的是三维时移量数据的一个时间切片，其中x方向样点数为701，y方向样点数分别为221。
本应用示例面向三维地震数据时差校正需求，针对两个数据能量关系不一致时已有二维动态拉伸方法存在的问题，提供了一种更适用于实际数据处理的三维地震数据动态拉伸时差校正方法，对于能量关系差异较大的两个三维地震数据先进行能量一致性处理，提高了二维时移量时差校正的准确度，从而获得较高精度的动态时移量，进而将两个数据校正为一致，通过逐线处理，将二维时移量时差校正方法扩展到三维，使之更好的服务于实际应用之中。
实施例2
根据本发明的实施例，提供了一种三维地震数据动态拉伸时差校正系统，系统包括：存储器，存储有计算机可执行指令；处理器，运行存储器上的计算可执行指令时，处理器实现以下步骤：
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)。
本实施例面向三维地震数据时差校正需求，针对两个数据能量关系不一致时已有二维动态拉伸方法存在的问题，提供了一种更适用于实际数据处理的三维地震数据动态拉伸时差校正方法，对于能量关系差异较大的两个三维地震数据先进行能量一致性处理，提高了二维时移量时差校正的准确度，从而获得较高精度的动态时移量，进而将两个数据校正为一致，通过逐线处理，将二维时移量时差校正方法扩展到三维，使之更好的服务于实际应用之中。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

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所述的三维地震数据动态拉伸时差校正方法，其中，平滑处理后的三维时移量数据表示为：
其中，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。
7.根据权利要求1所述的三维地震数据动态拉伸时差校正方法，其中，步骤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。
8.一种三维地震数据动态拉伸时差校正系统，其特征在于，所述系统包括：
存储器，存储有计算机可执行指令；
处理器，运行所述存储器上的计算可执行指令时，所述处理器实现以下步骤：
1)分别逐线读取第一组三维地震数据和第二组三维地震数据，计算两组三维数据的每个点的时移量，获得三维时移量数据；
2)对所述三维时移量数据进行平滑处理；
3)利用经平滑处理的时移量数据对于所述第二组三维地震数据进行时差校正，
其中，步骤3)包括：
3.1)读取第二组三维地震数据的其中一条线数据D2(nx,nz)和相应的经平滑处理后的三维时移量数据的线数据Ts(nx,nz)；
3.2)对于所述第二组三维地震数据的线数据D2(nx,nz)的每一个点，利用相应的经平滑处理后的三维时移数据量Ts(nx,nz)进行时差校正；
3.3)读取第二组三维地震数据和经平滑处理后的三维时移量数据的下一条线数据进行时差校正，直到第二组三维地震数据的所有线数据完成校正。
9.根据权利要求8所述的三维地震数据动态拉伸时差校正系统，其中，所述第一组三维地震数据和第二组三维地震数据的维度包括：线、道和深度点，每个所述维度的样点数分别表示为ny、nx和nz。
10.根据权利要求8所述的三维地震数据动态拉伸时差校正系统，其中，步骤1)包括：
1.1)读取第二组三维地震数据其中一条线数据D2(nx,nz)和相应的第一组三维地震数据的线数据D1(nx,nz)；
1.2)计算所述第二组三维地震数据的线数据D2(nx,nz)相对于所述第一组三维地震数据的线数据D1(nx,nz)的每个点的时移量，将所述时移量按线号位置存放在三维数组中；
1.3)读取第二组三维地震数据和第一组三维地震数据的下一条线数据，继续进行时移量计算，直到第二组三维地震数据的所有线数据完成时移量计算，获得所述三维时移量数据T(ny,nx,nz)。
CN201710935821.8A 2017-10-10 2017-10-10 三维地震数据动态拉伸时差校正方法及系统 Active CN109655912B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710935821.8A CN109655912B (zh) 2017-10-10 2017-10-10 三维地震数据动态拉伸时差校正方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710935821.8A CN109655912B (zh) 2017-10-10 2017-10-10 三维地震数据动态拉伸时差校正方法及系统

Publications (2)

Publication Number Publication Date
CN109655912A true CN109655912A (zh) 2019-04-19
CN109655912B CN109655912B (zh) 2021-05-25

Family

ID=66109433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710935821.8A Active CN109655912B (zh) 2017-10-10 2017-10-10 三维地震数据动态拉伸时差校正方法及系统

Country Status (1)

Country Link
CN (1) CN109655912B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112649873A (zh) * 2019-10-10 2021-04-13 中国石油化工股份有限公司 地震数据分频动校正处理方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014172565A1 (en) * 2013-04-17 2014-10-23 Drilling Info, Inc. System and method for automatically correlating geologic tops
CN105182420A (zh) * 2015-10-13 2015-12-23 中国石油天然气集团公司 一种动态匹配动校正方法
CN106291707A (zh) * 2015-06-01 2017-01-04 中国石油化工股份有限公司 三维地震资料多块拼接处理时差校正方法和装置
CN106353818A (zh) * 2016-08-12 2017-01-25 中国石油天然气集团公司 一种地震数据时差校正的方法及系统
CN107179552A (zh) * 2016-03-11 2017-09-19 中国石油化工股份有限公司 一种基于波形动态匹配的子波拉伸校正处理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014172565A1 (en) * 2013-04-17 2014-10-23 Drilling Info, Inc. System and method for automatically correlating geologic tops
CN106291707A (zh) * 2015-06-01 2017-01-04 中国石油化工股份有限公司 三维地震资料多块拼接处理时差校正方法和装置
CN105182420A (zh) * 2015-10-13 2015-12-23 中国石油天然气集团公司 一种动态匹配动校正方法
CN107179552A (zh) * 2016-03-11 2017-09-19 中国石油化工股份有限公司 一种基于波形动态匹配的子波拉伸校正处理方法
CN106353818A (zh) * 2016-08-12 2017-01-25 中国石油天然气集团公司 一种地震数据时差校正的方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAVE HALE: "Dynamic warping of seismic images", 《GEOPHYSICS》 *
蒋雪珍等: "基于动态图像变形的PP与PS波层位直接匹配", 《地球物理学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112649873A (zh) * 2019-10-10 2021-04-13 中国石油化工股份有限公司 地震数据分频动校正处理方法及系统
CN112649873B (zh) * 2019-10-10 2024-04-09 中国石油化工股份有限公司 地震数据分频动校正处理方法及系统

Also Published As

Publication number Publication date
CN109655912B (zh) 2021-05-25

Similar Documents

Publication Publication Date Title
CN106504276B (zh) 非局部立体匹配方法
US8125493B2 (en) Generating a texture from multiple images
US8331686B2 (en) Method and apparatus for extracting line segments from image using incremental pixel extension
CN101655982B (zh) 基于改进Harris角点的图像配准方法
CN104463894A (zh) 一种多视角三维激光点云全局优化整体配准方法
CN102721376A (zh) 一种大视场三维视觉传感器的标定方法
EP2277001B1 (en) Method and system for positioning an object on an optical sensor system for acquiring a topography thereof
CN106780610B (zh) 一种位置定位方法及装置
CN103080973A (zh) 立体图像处理装置及立体图像处理方法
Winkler et al. Actin filament tracking in electron tomograms of negatively stained lamellipodia using the localized radon transform
CN109655912A (zh) 三维地震数据动态拉伸时差校正方法及系统
CN115272080A (zh) 基于图像拼接的全局变形测量方法和系统
Qin et al. A curvature constraint Exemplar-based image inpainting
Yin et al. Initializing and accelerating Stereo-DIC computation using semi-global matching with geometric constraints
Ma et al. Variable subset DIC algorithm for measuring discontinuous displacement based on pixel-level ZNCC value distribution map
JP2007322404A (ja) 画像処理装置及びその処理方法
US6961459B2 (en) Three-dimensional data processing device, three-dimensional data processing method, and program providing medium
CN106228593B (zh) 一种图像密集匹配方法
CN110376644A (zh) 叠前地震数据拼接方法及装置
CN115727836A (zh) 一种基于激光反光板和里程计的融合定位方法及系统
AU2021201735B2 (en) Representing an interior of a volume
CN112070810B (zh) 定位方法、可移动设备及计算机可读存储介质
CN109782216A (zh) 一种二维空间谱峰值搜索的简易方法
Wöhler et al. Monocular 3D scene reconstruction at absolute scale
CN116934956A (zh) 三维模型的构建方法、装置、设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant