CN114114417B - 地震成像方法及装置 - Google Patents

地震成像方法及装置 Download PDF

Info

Publication number
CN114114417B
CN114114417B CN202010861535.3A CN202010861535A CN114114417B CN 114114417 B CN114114417 B CN 114114417B CN 202010861535 A CN202010861535 A CN 202010861535A CN 114114417 B CN114114417 B CN 114114417B
Authority
CN
China
Prior art keywords
offset
point
imaging
target
sub
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.)
Active
Application number
CN202010861535.3A
Other languages
English (en)
Other versions
CN114114417A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010861535.3A priority Critical patent/CN114114417B/zh
Publication of CN114114417A publication Critical patent/CN114114417A/zh
Application granted granted Critical
Publication of CN114114417B publication Critical patent/CN114114417B/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/34—Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
    • 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
    • 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
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • 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/51—Migration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

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

本申请提供了一种地震成像方法及装置，该方法包括：获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。本申请能够提高复杂高陡构造地层的地震成像的准确性，进而能够提高油气藏勘探的精度。

Description

地震成像方法及装置
技术领域
本申请涉及地震数据处理技术领域，尤其涉及一种地震成像方法及装置。
背景技术
在油气藏勘探领域，高陡复杂构造及其油气藏十分重要，高陡复杂构造的准确成像是高陡构造及其油气藏地震勘探面临的关键问题之一；随着地震勘探向地质条件复杂地区的发展，复杂构造的准确成像成为亟需解决的问题。
在野外，地震数据大都是在不同偏移距的接收点上采集的。由于地下构造的复杂性，地震波在传播过程中，随着偏移距的不同，传播路径及其经过的地层在横向和纵向上均存在着差异，有时该差异还很大，这必然会导致波传播速度的差异。常规叠加中所采用的道集是由不同偏移距组合而成的，叠加时大都是对整个共中心点道集采用相同的速度曲线，忽略了近道、远道和非双曲线等不同因素而造成的动态时差，以及不同偏移距道对动校正响应的差异，使校正后的道集内仍存在着剩余时差，从而影响共深度点的同相叠加。同时，由于地质构造的复杂性，基于水平层状介质的共中心点叠加技术存在先天缺陷，主要是共中心点道集不再完全是共中心点下方的反射，而是来自不同构造部位的反射。对于高陡复杂构造的成像，共中心点道集不再完全是共中心点下方的反射，是无法准确求解高陡构造速度和实现同向叠加。因此，对于高陡复杂构造，共中心点道集不再完全是共中心点下方的反射，求解高陡构造速度和叠加成像的准确性较低。
发明内容
针对现有技术中的问题，本申请提出了一种地震成像方法及装置，能够提高复杂高陡构造地层的地震成像的准确性，进而能够提高油气藏勘探的精度。
为了解决上述技术问题，本申请提供以下技术方案：
第一方面，本申请提供一种地震成像方法，包括：
获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；
从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；
从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
进一步地，所述获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱，包括：获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据；根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
进一步地，所述分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果，包括：生成各个所述分偏移距组段各自对应的分偏移距速度谱；根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图；基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
进一步地，所述从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，包括：从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
进一步地，所述基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图，包括：根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数；应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数；根据所述目标逐点函数，生成所述目标二维地震测线对应的目标速度场；基于所述目标速度场，得到所述目标叠加成像剖面图。
进一步地，在所述从所述共中心点道集地震数据中提取多个分偏移距组段之后，还包括：分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果；从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
第二方面，本申请提供一种地震成像装置，包括：
获取模块，用于获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；
第一叠加成像模块，用于从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；
目标叠加成像模块，用于从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
进一步地，所述获取模块，包括：获取单元，用于获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据；生成单元，用于根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
进一步地，所述第一叠加成像模块，包括：生成速度谱单元，用于生成各个所述分偏移距组段各自对应的分偏移距速度谱；生成剖面图单元，用于根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图；第一叠加成像单元，用于基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
进一步地，所述目标叠加成像模块包括：获取剖面单元，用于从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
进一步地，所述目标叠加成像模块包括：生成逐点速度单元，用于根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数；生成目标逐点速度单元，用于应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数；生成目标速度场单元，用于根据所述目标逐点函数，生成所述目标二维地震测线对应的目标速度场；目标叠加成像单元，用于基于所述目标速度场，得到所述目标叠加成像剖面图。
进一步地，所述的地震成像装置，还包括：第二叠加成像模块，用于分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果；地震成像模块，用于从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
第三方面，本申请提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的地震成像方法。
第四方面，本申请提供一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现所述的地震成像方法。
由上述技术方案可知，本申请提供一种地震成像方法及装置。其中，该方法包括：获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图，能够提高复杂高陡构造地层的地震成像的准确性，进而能够提高油气藏勘探的精度；具体地，能够避免高陡构造地层叠加速度和中心点叠加多解的问题，可以很好的获得准确的高陡复杂构造成像，能够提高油气藏勘探的精度和地震剖面信噪比。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是现有技术中平缓地层的叠加速度分析的原理示意图；
图2是现有技术中高陡构造地层的叠加速度分析的原理示意图；
图3是本申请实施例中地震成像方法的流程示意图；
图4是本申请实施例中地震成像方法中步骤101和步骤102的流程示意图；
图5是本申请实施例中地震成像方法中步骤201至步骤203的流程示意图；
图6是本申请另一实施例中地震成像方法的流程示意图；
图7是本申请又一实施例中地震成像方法的流程示意图；
图8是本申请实施例中地震成像方法中步骤210和步骤220的流程示意图；
图9是本申请具体应用实例中的偏移距范围0-10000m的全偏移距数据速度谱和MVFS速度扫描叠加段之间的比较示意图；
图10示出了本申请具体应用实例中的全偏移距数据叠加剖面；
图11示出了本申请具体应用实例中的0-4000m偏移距数据叠加剖面；
图12示出了本申请具体应用实例中的4000-7000m偏移距数据叠加剖面；
图13示出了本申请具体应用实例中的7000-10000m偏移距数据叠加剖面；
图14是本申请具体应用实例中的偏移距范围4000-10000m的远偏移距数据速度谱和MVFS速度扫描叠加段之间的比较示意图；
图15示出了现有技术中的叠加成像剖面；
图16示出了本申请具体应用实例中的叠加成像剖面；
图17是本申请实施例中地震成像装置的结构示意图；
图18是本申请实施例的电子设备的系统构成示意框图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
目前，在共中心点的道集上，以动校正能否将反射时距曲线校平和叠加能量是否增强为判别标准，是所有地震数据处理中叠加成像普遍采用的技术，但是，对于高陡复杂构造的成像，共中心点道集不再完全是共中心点下方的反射，无法准确求解高陡构造速度和实现同向叠加成像。
现有技术中叠加速度分析的原理如图1所示，r1和r2表示不同的反射波，当提前假设条件就是地层为水平层状介质时，其地下反射点与中心点位置在同一个位置，其获得的反射波同向轴符合双曲线的特征，这样获得的速度和叠加是一一对应的，这样就可以得到准确叠加成像。但是当地层存在倾角(即有高陡构造存在)的情况下，如图2所示，r1和r2表示不同的反射波，反射点位置与中心点位置将不在是同一个位置点，而且从图2可以看出，偏移距包括：检波点1相对于炮点1的偏移距，中心点为中心点1；检波点2相对于炮点2的偏移距，中心点为中心点2；偏移距越大，中心点与地下反射点位置差异就越大，这样就不能简单的利用常规的经典理论用于高陡构造的速度分析及叠加成像。如图2所示，复杂构造越陡，采集到地层的地震信息就更多的集中于中远偏移距上，因此，为了提高求取复杂高陡构造地层的速度与成像精度，本申请考虑从改变现有的地震成像方式入手，在地震数据处理阶段，利用不同偏移距包含不同地层的地震波优势信息的特征，建立对应的速度谱，可以避免速度多解性的干扰，准确求取高陡构造地层和整个剖面的速度，这样就可以使高陡构造准确成像并提高地震剖面信噪比。本申请通过确定高陡构造的地震信息主要集中的优势偏移距组段，利用优势偏移距组段进行速度分析，能够准确得到高陡构造地层的速度，进而应用此速度再进行高陡复杂构造区的叠加成像，能够提高高陡构造地层的成像精度。
具体通过下述各个实施例进行说明。
为了提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度，本实施例提供一种执行主体是地震成像装置的地震成像方法，该地震成像装置包括但不限于服务器，如图3所示，该方法具体包含有如下内容：
步骤100：获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱。
具体地，可以通过观测系统定义、层析静校正、异常振幅衰减、地表一致性反褶积和剩余静校正等处理方式，得到目标二维地震测线的叠加剖面和全偏移距的共中心点道集地震数据；可以根据所述共中心点道集地震数据生成该目标二维地震测线的全偏移距速度谱。
其中，二维地震勘探需要在地面上布置一条条的测线，沿各条测线进行地震勘探施工，该测线即为二维地震测线；所述全偏移距为所述目标二维地震测线对应的全部偏移距。
步骤200：从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果。
具体地，所述分偏移距组段可以是近中偏移距组段或偏移距组段等，可以根据预先设定的偏移距组段划分标准提取各个所述分偏移距组段，所述偏移距组段划分标准可以根据实际需要设置，也可以根据以下步骤生成所述偏移距组段划分标准：获得动校正后的共中心点道集地震数据和全偏移距速度谱，依据动校正后的共中心点道集地震数据的拉平程度，以及不同偏移距数据包含不同深度地层信息的优势程度，还有高陡构造成像的位置及上覆地层倾角大小等因素，确定相应的偏移距组段划分标准。
步骤300：从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
具体地，可以根据预设的成像效果评价规则，确定最佳的第一叠加成像结果，该成像效果评价规则可根据实际需要进行设置。
其中，步骤300可以分为：步骤310：从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段；步骤320：基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
为了进一步提高获取全偏移距速度谱的准确性，进而提高接下来应用该全偏移距速度谱和分偏移距速度谱实现地震成像的可靠性，在本申请一个实施例中，参见图4，步骤100包含有：
步骤101：获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据。
步骤102：根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
为了准确获取避免高陡构造地层对应的多个叠加成像结果，进而提高接下来应用所述叠加成像结果确定高陡构造地层的成像效果最佳的叠加成像剖面的准确性，在本申请一个实施例中，参见图5，步骤200所述的分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果，包含有：
步骤201：生成各个所述分偏移距组段各自对应的分偏移距速度谱。
步骤202：根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图。
可以理解的是，所述分偏移距速度谱和叠加成像剖面图为一一对应关系。
步骤203：基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
为了避免高陡构造地层叠加速度和中心点叠加多解的问题，获取高陡构造地层的成像效果最佳的叠加成像剖面，进而提高地震成像的准确性，在本申请一个实施例中，参见图6，步骤310包含有：
步骤311：从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
为了进一步提高生成目标叠加成像剖面图的准确性，在本申请一个实施例中，参见图7，步骤320包含有：
步骤321：根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数。
步骤322：应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数。
具体地，将所述第一逐点速度函数和第二逐点速度函数进行合并，合并过程中，以第一逐点速度函数为母本，把分偏移距组段对应的叠加成像剖面图中高陡构造地层的第一逐点速度函数替换为第二逐点速度函数，形成目标逐点速度函数。
步骤323：根据所述目标逐点函数，生成所述目标二维地震测线对应的目标速度场。
步骤324：基于所述目标速度场，得到所述目标叠加成像剖面图。
具体地，可以根据所述目标速度场的速度，对全偏移距进行叠加成像，得到所述目标叠加成像剖面图。
为了进一步提高地震成像的准确性，同时考虑高陡构造地层和平缓地层的速度对地震成像的影响，本申请一个实施例中，参见图8，在步骤200所述的从所述共中心点道集地震数据中提取多个分偏移距组段之后，还包括：
步骤210：分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果。
具体地，所述分偏移距组段和第二叠加成像结果是一一对应关系。
步骤220：从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
为了进一步说明本方案，本申请提供一种图像成像方法的具体应用实例，在地震数据采集时，分布于不同偏移距段的检波器会接收到来自不同深度位置的地震波信息，尤其在存在高陡复杂构造的地区，其下伏构造的地震波信号更多的会被分布于中、远偏移距的检波器接收到。基于此，利用这一特征来提高高陡构造的速度分析准确率和叠加成像精度，具体包括以下步骤：
S11：获取常规的全偏移距CMP共中心点道集地震数据。
S12：将步骤S11的全偏移距CMP共中心点道集地震数据生成全偏移距速度谱并进行速度v1拾取，获得应用v1动校正后的全偏移距CMP超道集数据和对应的速度谱，依据动校正后的CMP超道集的拉平程度，以及高陡构造地层成像的位置及上覆地层倾角大小等，确定相应的偏移距组段划分标准。
具体地，将步骤S11的全偏移距CMP共中心点道集地震数据生成全偏移距速度谱并进行速度拾取，获得动校正后的全偏移距CMP超道集和速度谱，依据动校正后的CMP超道集的拉平程度，以及不同偏移距数据包含不同深度地层信息的优势程度，还有高陡构造成像的位置及上覆地层倾角大小等，确定相应的偏移距组段划分标准。
S13：依据步骤S12的偏移距组划分标准，将常规的全偏移距CMP共中心点道集地震数据按偏移距大小分选为近中及中远偏移距组段，得到不同偏移距组段的地震数据，即分偏移距CMP道集数据，同时生成相应的速度谱，在相应的速度谱上拾取速度，得到相应的分偏移距组段的叠加成像速度。
具体地，依据步骤S12的偏移距组划分标准，将常规的全偏移距CMP共中心点道集地震数据分选为近中及中远的偏移距组段，得到不同偏移距组段的地震数据，并生成相应的速度谱，并在相应的速度谱上拾取得到相应的分偏移距组段的叠加成像速度。
S14：由步骤S13得到的分偏移距CMP道集数据和相应的叠加速度，通过叠加方式获得相应的近中及中远偏移距的叠加成像剖面。然后通过对比分析，确定高陡构造地层的叠加成像是否为最佳，如果为最佳叠加成像，则确定偏移距分组范围及叠加成像速度，否则返回至步骤S12和步骤S13，直至确定最优的偏移距分组范围和获得准确叠加速度为止。
S15：对利用由步骤S14得到的最优近中偏移距叠加与中远偏移距叠加进行对比分析，近中偏移距的叠加成像在平缓地层成像效果较好，而高陡构造地层域的中远偏移距叠加成像则更好。依据这一认识，在叠加剖面上确定平缓地层和高陡构造地层成像优势范围。
S16：利用由步骤S15确定的平缓地层与高陡构造地层的成像优势范围，相应的选取各自的成像速度，就是把促使高陡构造地层和平缓地层精确的成像优势范围的对应速度选取出来，而后合成插值为一个全地震剖面成像的速度场。
S17：在获得最优叠加成像速度场的基础上，在步骤S11获得的全偏移距CMP道集上，进行高精度叠加，这样可以避免由于地层复杂引起的速度多解性，同时又可获得全剖面的准确地震成像。
具体地，以不同偏移距组段包含不同深度地层的地震信号优势程度为依据，首先对全偏移距的地震道数据进行速度分析，依据动校正后的CMP超道集的拉平程度，以及高陡构造成像的位置及上覆地层倾角大小等因素，确定相应的偏移距组段初步划分标准。依据初步划分标准，将常规的CMP道集地震数据分选为近中及中远的偏移距组段，并在相应的速度谱上拾取得到相应的分偏移距组段的叠加速度的对应叠加成像，然后通过对比分析，确定高陡构造的成像是否为最佳，如果为最佳叠加成像，确定偏移距分组及叠加成像速度，否则重新分组，直至确定最优的偏移距分组和准确叠加速度为止。基于此，通过反复迭代，获得最终全剖面叠加速度和最终的准确成像。能够避免高陡构造地层叠加速度和中心点叠加多解的问题，可以获得准确的高陡复杂构造成像结果，进而能够提高油气藏勘探的精度。
从上述描述可知，本应用实例提供的地震成像方法，利用互相关方法及全偏移距CMP道集数据建立速度谱，并分析速度谱获得初始叠加速度；由于地震叠加成像基于水平层状和双曲线的假设，于是基于叠加能量相干最大的原则，利用初始速度及不同偏移距段的数据进行分偏移距叠加，进而确定包含高陡构造地层成像信息的有效偏移距组段；在确定高陡构造叠加成像的优势偏移距组段基础上，利用包含高陡构造地震反射信息的偏移距组段建立速度谱，拾取获得高陡构造下伏地层的准确速度；把平缓地层的叠加速度与准确的高陡叠加成像速度进行构造约束下的速度融合，进而形成整个工区的叠加速度场；在准确的叠加速度场基础上，利用全偏移距数据进行叠加处理并获得最终成像数据。本应用实例通过不同偏移距信息包含不同的地球物理信息为基础，实现了高陡构造下伏地层的速度准确求解和准确成像，从而提高高陡构造的成像质量。在不同地震地质信息基础上，利用不同偏移距包含不同优势的地球物理信息，尽可能的减少速度多解性的问题，以便求解出准确的速度从而刻画高陡构造的成像。
为了进一步说明本方案地震成像过程的实用性和可行性，本申请提供一种地震成像方法的具体应用实例，在塔里木盆地的库车地区，已发现多个大型油气田，是塔里木盆地探明油气储量最多的地区之一。本区因地表复杂(如山体陡峭或沟壑纵横等)和地下构造复杂(如逆冲推覆或断块破碎等)，采集的地震数据信噪比极低。且由于高陡构造、复杂断裂带与刺穿盐体交织发育，形成异常复杂的地下地质结构，呈现很强的速度横向变化和陡倾角地层分布，地震数据成像异常困难。鉴于库车坳陷地质目标储层埋藏较深，并与刺穿盐体交织发育，叠前地震数据理的重点是远偏移距上的广角反射信号的恢复和成像速度的准确求取。这样才能获得高陡构造下的较好成像。下面以库车东秋地区的一条二维测线为例。
S21：二维地震测线进行常规处理，通过观测系统定义、层析静校正、异常振幅衰减、地表一致性反褶积和剩余静校正等处理方式中的至少一种，得到该测线的叠加剖面和全偏移距CMP道集数据。
S22：利用步骤S21得到的全偏移距CMP道集数据生成该二维测线的速度谱，参见图9，可以通过偏移距范围0-10000m的全偏移距数据速度谱和MVFS速度扫描叠加段之间的比较和分析，拾取得到一个叠加的初始速度，获得对应得逐点速度TV函数V1，然后通过插值方法得到一个全速度场，在图9和图14中，数据速度谱即Semblance速度谱的横坐标为速度，单位：米/秒，纵坐标为时间，单位：秒；MVFS速度扫描叠加段对应的横坐标为道，纵坐标为时间，单位：秒。
S23：将步骤S21的得到的全偏移距CMP道集数据进行偏移距分组，分成近中远的偏移距组段，而后利用步骤S22得到的全偏移距初始速度场进行近中远偏移距组段的叠加，参见图10，得到全偏移距(0-10000m)数据叠加剖面，参见图11，得到0-4000m偏移距数据叠加剖面，参见图12，得到4000-7000m偏移距数据叠加剖面，参见图13，得到7000-10000m偏移距数据叠加剖面，从而分析不同偏移距组段对高陡构造主体部位的成像贡献率，其中，图10至图13，以及图15和图16的横坐标为CMP号，纵坐标为时间，单位：秒。
S24：由步骤S23确定的最优复杂构造主体部位叠加成像的偏移距组合，在最优偏移距组合基础上，利用自相关的求取速度谱的方法再生成一个分偏移距的速度谱，如图14所示，可以通过偏移距范围4000-10000m的远偏移距数据速度谱即Semblance速度谱和MVFS速度扫描叠加段之间的比较和分析，拾取速度。
S25：在步骤S24获得的速度谱基础上，拾取速度，获得对应的逐点速度TV函数V2，同时与步骤S22由全偏移距数据获得的逐点速度TV函数V1进行合并，合并过程中，以V1速度函数为母本，把剖面中高陡构造地层的V1速度函数替换为V2，形成最终逐点速度TV函数V3，逐点速度TV函数V1至逐点速度TV函数V3可以为逐点速度TV函数对V1至逐点速度TV函数对V3。
S26：利用由步骤S25得到逐点速度TV函数V3，生成整个剖面的最终速度场，而后利用最新的速度对全偏移距进行叠加成像，如图16所示，可以看到本具体应用实例的叠加成像剖面相较于图15中的常规方法叠加成像剖面，在主体构造部分信噪比明显的提高。
从软件层面来说，为了提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度，本申请提供一种用于实现所述地震成像方法中全部或部分内容的地震成像装置的实施例，参见图17，所述地震成像装置具体包含有如下内容：
获取模块10，用于获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱。
第一叠加成像模块20，用于从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果。
目标叠加成像模块30，用于从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
在本申请一个实施例中，所述获取模块包括：
获取单元，用于获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据。
生成单元，用于根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
在本申请一个实施例中，所述第一叠加成像模块，包括：
生成速度谱单元，用于生成各个所述分偏移距组段各自对应的分偏移距速度谱；
生成剖面图单元，用于根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图。
第一叠加成像单元，用于基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
在本申请一个实施例中，所述目标叠加成像模块包括：
获取剖面单元，用于从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
在本申请一个实施例中，所述目标叠加成像模块包括：
生成逐点速度单元，用于根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数。
生成目标逐点速度单元，用于应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数。
生成目标速度场单元，用于根据所述目标逐点函数，生成所述目标二维地震测线对应的目标速度场。
目标叠加成像单元，用于基于所述目标速度场，得到所述目标叠加成像剖面图。
在本申请一个实施例中，所述的地震成像装置还包括：
第二叠加成像模块，用于分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果。
地震成像模块，用于从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
本说明书提供的地震成像装置的实施例具体可以用于执行上述地震成像方法的实施例的处理流程，其功能在此不再赘述，可以参照上述地震成像方法实施例的详细描述。
由上述描述可知，本申请提供的地震成像方法及装置，能够提高复杂高陡构造地层的地震成像的准确性，进而能够提高油气藏勘探的精度；具体地，能够避免高陡构造地层叠加速度和中心点叠加多解的问题，可以很好的获得准确的高陡复杂构造成像，能够提高油气藏勘探的精度和地震剖面信噪比。
从硬件层面来说，为了提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度，本申请提供一种用于实现所述地震成像方法中的全部或部分内容的电子设备的实施例所述电子设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现所述地震成像装置以及用户终端等相关设备之间的信息传输；该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例用于实现所述地震成像方法的实施例及用于实现所述地震成像装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
图18为本申请实施例的电子设备9600的系统构成的示意框图。如图18所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图18是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
在本申请一个或多个实施例中，地震成像功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：
步骤100：获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱。
步骤200：从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果。
步骤300：从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
从上述描述可知，本申请的实施例提供的电子设备，能够提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度。
在另一个实施方式中，地震成像装置可以与中央处理器9100分开配置，例如可以将地震成像装置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现地震成像功能。
如图18所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图18中所示的所有部件；此外，电子设备9600还可以包括图18中没有示出的部件，可以参考现有技术。
如图18所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
上述描述可知，本申请的实施例提供的电子设备，能够提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度。
本申请的实施例还提供能够实现上述实施例中的地震成像方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的地震成像方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤100：获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱。
步骤200：从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果。
步骤300：从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
从上述描述可知，本申请实施例提供的计算机可读存储介质，能够提高复杂高陡构造地层的地震成像的准确性，进而提高油气藏勘探的精度。
本申请中上述方法的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。相关之处参见方法实施例的部分说明即可。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本申请中应用了具体实施例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。

Claims (12)

1.一种地震成像方法，其特征在于，包括：
获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；
从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；
从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图；
所述基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图，包括：
根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数；
应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数；
根据所述目标逐点速度函数，生成所述目标二维地震测线对应的目标速度场；
基于所述目标速度场，得到所述目标叠加成像剖面图；
所述应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数，包括：将所述第一逐点速度函数和第二逐点速度函数进行合并，合并过程中，以第一逐点速度函数为母本，把分偏移距组段对应的叠加成像剖面图中高陡构造地层的第一逐点速度函数替换为第二逐点速度函数，形成目标逐点速度函数。
2.根据权利要求1所述的地震成像方法，其特征在于，所述获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱，包括：
获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据；
根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
3.根据权利要求1所述的地震成像方法，其特征在于，所述分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果，包括：
生成各个所述分偏移距组段各自对应的分偏移距速度谱；
根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图；
基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
4.根据权利要求3所述的地震成像方法，其特征在于，所述从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，包括：
从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
5.根据权利要求1所述的地震成像方法，其特征在于，在所述从所述共中心点道集地震数据中提取多个分偏移距组段之后，还包括：
分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果；
从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
6.一种地震成像装置，其特征在于，包括：
获取模块，用于获取目标二维地震测线针对于全偏移距的共中心点道集地震数据和全偏移距速度谱；
第一叠加成像模块，用于从所述共中心点道集地震数据中提取多个分偏移距组段，分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的高陡构造地层的第一叠加成像结果；
目标叠加成像模块，用于从各个所述分偏移距组段中，获取最佳的第一叠加成像结果对应的第一分偏移距组段，基于该第一分偏移距组段的分偏移距速度谱和所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图；
所述目标叠加成像模块包括：
生成逐点速度单元，用于根据所述全偏移距速度谱得到第一逐点速度函数，以及根据所述第一分偏移距组段的分偏移距速度谱，得到所述高陡构造地层对应的第二逐点速度函数；
生成目标逐点速度单元，用于应用所述第一逐点速度函数和第二逐点速度函数，生成目标逐点速度函数；
生成目标速度场单元，用于根据所述目标逐点速度函数，生成所述目标二维地震测线对应的目标速度场；
目标叠加成像单元，用于基于所述目标速度场，得到所述目标叠加成像剖面图；
所述生成目标逐点速度单元，具体用于将所述第一逐点速度函数和第二逐点速度函数进行合并，合并过程中，以第一逐点速度函数为母本，把分偏移距组段对应的叠加成像剖面图中高陡构造地层的第一逐点速度函数替换为第二逐点速度函数，形成目标逐点速度函数。
7.根据权利要求6所述的地震成像装置，其特征在于，所述获取模块，包括：
获取单元，用于获取目标二维地震测线针对于所述全偏移距的共中心点道集地震数据；
生成单元，用于根据所述共中心点道集地震数据，生成所述全偏移距速度谱。
8.根据权利要求6所述的地震成像装置，其特征在于，所述第一叠加成像模块，包括：
生成速度谱单元，用于生成各个所述分偏移距组段各自对应的分偏移距速度谱；
生成剖面图单元，用于根据所述分偏移距速度谱，生成所述第一分偏移距组段的叠加成像剖面图；
第一叠加成像单元，用于基于所述叠加成像剖面图，得到所述高陡构造地层的第一叠加成像结果，所述叠加成像剖面与第一叠加成像结果为一一对应关系。
9.根据权利要求8所述的地震成像装置，其特征在于，所述目标叠加成像模块包括：
获取剖面单元，用于从各个所述叠加成像剖面中，获取所述高陡构造地层的成像效果最佳的叠加成像剖面，将该叠加成像剖面对应的分偏移距组段作为所述第一分偏移距组段。
10.根据权利要求6所述的地震成像装置，其特征在于，还包括：
第二叠加成像模块，用于分别应用各个所述分偏移距组段，得到所述目标二维地震测线对应的平缓地层的第二叠加成像结果；
地震成像模块，用于从各个所述分偏移距组段中，确定最佳的第二叠加成像结果对应的第二分偏移距组段，基于所述第一分偏移距组段和第二分偏移距组段各自的分偏移距速度谱，以及所述全偏移距速度谱，生成所述全偏移距对应的目标叠加成像剖面图。
11.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至5任一项所述的地震成像方法。
12.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述指令被执行时实现权利要求1至5任一项所述的地震成像方法。
CN202010861535.3A 2020-08-25 2020-08-25 地震成像方法及装置 Active CN114114417B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010861535.3A CN114114417B (zh) 2020-08-25 2020-08-25 地震成像方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010861535.3A CN114114417B (zh) 2020-08-25 2020-08-25 地震成像方法及装置

Publications (2)

Publication Number Publication Date
CN114114417A CN114114417A (zh) 2022-03-01
CN114114417B true CN114114417B (zh) 2024-03-01

Family

ID=80373789

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010861535.3A Active CN114114417B (zh) 2020-08-25 2020-08-25 地震成像方法及装置

Country Status (1)

Country Link
CN (1) CN114114417B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076623A (zh) * 2011-10-25 2013-05-01 中国石油化工股份有限公司 一种基于叠前相干的裂缝检测方法
CN104614762A (zh) * 2014-12-17 2015-05-13 中国石油天然气股份有限公司 疏松砂岩气藏边界确定方法及装置
CN107179546A (zh) * 2017-05-11 2017-09-19 中国石油天然气股份有限公司 一种确定速度场的方法及装置
CN109917461A (zh) * 2019-03-04 2019-06-21 中国石油天然气集团有限公司 一种提高侏罗系地层及基岩成像质量的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076623A (zh) * 2011-10-25 2013-05-01 中国石油化工股份有限公司 一种基于叠前相干的裂缝检测方法
CN104614762A (zh) * 2014-12-17 2015-05-13 中国石油天然气股份有限公司 疏松砂岩气藏边界确定方法及装置
CN107179546A (zh) * 2017-05-11 2017-09-19 中国石油天然气股份有限公司 一种确定速度场的方法及装置
CN109917461A (zh) * 2019-03-04 2019-06-21 中国石油天然气集团有限公司 一种提高侏罗系地层及基岩成像质量的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
分偏移距分频地震属性裂缝预测技术在页岩气裂缝预测中的应用;吕财;《江汉石油职工大学学报》;第31卷(第3期);第1-4页 *
牟书令 等.《中国海油气勘探理论技术与实践》.北京地质出版社,2009,第425-434页. *

Also Published As

Publication number Publication date
CN114114417A (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN105182408A (zh) 一种合成地震记录的制作方法和装置
CN112415591A (zh) 绕射波的成像方法、装置、电子设备及存储介质
CN103454681B (zh) 评价三维地震观测系统成像效果的方法和设备
CN109765615A (zh) 一种地层品质因子反演方法及装置
CN103018771A (zh) 确定拖缆形状的装置和方法
CN106125139B (zh) 一种三维地震数据处理方法及系统
Brodic et al. Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging
CN104199103A (zh) 转换波静校正方法及装置
CN111551989B (zh) 一种透射槽波成像方法、设备及计算机可读存储介质
Hung et al. Effects of failure of the ocean‐bottom seismograph leveling system on receiver function analysis
CN105093296A (zh) 一种优化观测系统的方法及装置
CN114114417B (zh) 地震成像方法及装置
WO2024119967A1 (zh) 一种多激发点瞬态面波勘探方法、设备及存储介质
CN111158050B (zh) 数据采集系统、方法及隧道地震波超前预报方法
CN111077577B (zh) 井地联合储层描述方法及装置
CN112255679B (zh) 地震资料绕射深度偏移处理方法及装置
CN112817042B (zh) 沉积储层中油气识别方法及装置
CN105445795A (zh) 倾斜界面的速度分析方法及其系统
WO2020001353A1 (zh) 消除井旁界面方位不确定性的反射横波测井系统及方法
CN111691876A (zh) 一种利用声波测井对邻井成像的方法、装置及存储介质
CN113325476B (zh) 地震资料静校正方法及装置
CN112505748B (zh) 炮域反射波拾取方法及装置
CN112888970A (zh) 确定声波慢度的方法和装置
CN112764107B (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