CN108802822B - 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 - Google Patents
方位各向异性介质中的保幅直接叠前时间偏移方法及装置 Download PDFInfo
- Publication number
- CN108802822B CN108802822B CN201810606448.6A CN201810606448A CN108802822B CN 108802822 B CN108802822 B CN 108802822B CN 201810606448 A CN201810606448 A CN 201810606448A CN 108802822 B CN108802822 B CN 108802822B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- azimuth angle
- seismic data
- migration
- wide
- 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.)
- Expired - Fee Related
Links
- 238000013508 migration Methods 0.000 title claims abstract description 149
- 230000005012 migration Effects 0.000 title claims abstract description 149
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000003384 imaging method Methods 0.000 claims abstract description 109
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000000523 sample Substances 0.000 claims description 7
- 230000015654 memory Effects 0.000 description 11
- 239000000243 solution Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 239000010410 layer Substances 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005526 G1 to G0 transition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000007983 Tris buffer Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- 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
- G01V2210/512—Pre-stack
-
- 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/58—Media-related
- G01V2210/586—Anisotropic media
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
技术领域
本发明涉及油气资源勘探中反射地震资料处理与成像技术领域，尤其涉及一种方位各向异性介质中的保幅直接叠前时间偏移方法及装置。
背景技术
地震勘探技术是地球物理探测方法中的重要手段。地震勘探是指人工激发所引起的弹性波利用地下介质弹性和密度的差异，通过观测和分析人工地震产生的地震波在地下的传播规律，推断地下岩层的性质和形态的地球物理勘探方法。目前偏移处理是地震资料处理中最重要的环节，其中叠前时间偏移是复杂构造成像最有效的方法之一，能适应纵横向速度变化较大的情况，适用于大倾角的偏移成像，在建立偏移速度模型、计算效率、灵活性等方面具有很多明显优势。
影响叠前时间偏移成像效果的因素包括：偏移速度、地震走时计算方法、偏移计算的偏移孔径、计算偏移幅值时的权系数、偏移算法实现流程。其中，走时计算与偏移速度共同决定了反射波能够正确归位，偏移孔径及其应用方式决定了偏移和偏移算法的计算量，权系数决定了成像幅值能否正确反映地下实际界面的物性参数变化，偏移算法实现流程对偏移的计算效率和存储需求有重要影响。随着对勘探精度要求的提高以及勘探深度的加深，考虑地层各向异性的研究成为叠前时间偏移成像的发展方向。HTI(HorizontalTransverse Isotropy)介质是描述各向同性介质中分布着一组平行的定向排列的垂直裂隙所构成的各向异性介质模型，它属于方位各向异性，地震波在该类介质中传播，速度随方向变化的特性不仅表现在随着相位角的变化而变化，而且随着观测方位的变化而变化。HTI介质代表了方位各向异性的实际地层(存在断裂、裂痕的地层)，在HTI介质中会出现横波分裂现象，当横波在HTI介质中按照一定方向传播时，将会分裂为一快一慢的两个横波，可见地震波在该类介质中的传播情况较为复杂。
随着地震采集技术的提高，宽方位角地震数据得到了普遍应用，因此有必要发展相应的方位各向异性三维叠前时间偏移技术。现有的叠前时间偏移方法在计算地震波走时、幅值并未考虑方位各向异性的影响，因此对宽方位角地震数据进行处理时出现成像问题。
发明内容
本发明提供一种方位各向异性介质中的保幅直接叠前时间偏移方法及装置，实现基于方位各向异性的情况下宽方位角地震数据的保幅成像。
本发明第一方面提供一种方位各向异性介质中的保幅直接叠前时间偏移方法，包括：
获取宽方位角地震数据，所述宽方位角地震数据指示方位各向异性特征；
根据预设方位角度将所述宽方位角地震数据进行分组，得到多个方位角区间对应的宽方位角地震数据；
通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
根据每个方位角区间对应的所述第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，所述偏移参数包括走时和振幅值；
根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
将所有方位角区间对应的第二偏移成像结果进行叠加，得到所述宽方位角地震数据的保幅偏移成像结果。
可选的，所述根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，包括：
根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
可选的，所述获取宽方位角地震数据，包括：
对地震数据进行方位角统计，绘制所述地震数据的方位角信息图；所述地震数据包括所有地震道的炮点和检波点的坐标；
根据所述方位角信息图，确定所述地震数据为宽方位角地震数据。
可选的，所述根据所述方位角信息图，确定所述地震数据为宽方位角地震数据，包括：
根据所述方位角信息图确定所述地震数据的方位角的纵横比；
判断所述地震数据的方位角的纵横比是否满足预设纵横比，若是，则确定所述地震数据为宽方位角地震数据。
可选的，所述根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果，包括：
根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。
本发明第二方面提供一种方位各向异性介质中的保幅直接叠前时间偏移装置，包括：
获取模块，用于获取宽方位角地震数据，所述宽方位角地震数据指示方位各向异性特征；
分组模块，用于根据预设方位角度将所述宽方位角地震数据进行分组；得到多个方位角区间对应的宽方位角地震数据；
计算成像模块，用于通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
确定模块，用于根据每个方位角区间对应的所述第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
计算模块，还用于根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，所述偏移参数包括走时和振幅值；
所述计算成像模块，用于根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
所述计算成像模块，还用于将所有方位角区间对应的第二偏移成像结果进行叠加，得到所述宽方位角地震数据的保幅偏移成像结果。
可选的，所述计算模块，具体用于根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
可选的，所述保幅直接叠前时间偏移装置还包括：
图像生成模块，用于对地震数据进行方位角统计，绘制所述地震数据的方位角信息图；所述地震数据包括所有地震道的炮点和检波点的坐标；
所述确定模块，还用于根据所述方位角信息图，确定所述地震数据为宽方位角地震数据。
可选的，所述确定模块，用于根据所述方位角信息图确定所述地震数据的方位角的纵横比；
判断所述地震数据的方位角的纵横比是否满足预设纵横比，若是，则确定所述地震数据为宽方位角地震数据。
可选的，所述计算成像模块，具体用于根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。
本发明第三方面提供一种方位各向异性介质中的保幅直接叠前时间偏移装置，包括：
存储器；
处理器；以及
计算机程序；
其中，所述计算机程序存储在所述存储器中，并被配置为由所述处理器执行以实现如本发明第一方面任一项所述的方法。
本发明第四方面提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行以实现如本发明第一方面任一项所述的方法。
本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移方法及装置，将获取到的宽方位角地震数据按照预设方位角度进行分组，通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的第一偏移成像结果，并确定每个方位角区间对应的偏移速度场；根据每个偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数；根据每个方位角区间对应的偏移参数得到得到每个方位角区间对应的第二偏移成像结果；将所有方位角区间对应的第二偏移成像结果进行叠加，最终得到宽方位角地震数据的保幅偏移成像结果。由于每个方位角区间对应的偏移成像结果考虑了方位各向异性对地震波走时和幅值的影响，实现各向异性介质中的保幅偏移成像。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移方法的流程示意图；
图2为本发明实施例提供的方位角信息图；
图3为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置的结构示意图一；
图4为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置的结构示意图二；
图5为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置的硬件结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
本发明实施例提供一种方向各向异性介质中的保幅直接叠前时间偏移方法，用于解决宽方位角地震数据带来的方位各向异性问题，能够有效消除同相轴层间抖动、不聚焦等问题，识别裂缝方位，为油气开发提供地球物理依据。
下面结合具体实施例对本发明提供的保幅直接叠前时间偏移方法做详细说明。
图1为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移方法的流程示意图，如图1所示，本实施例的保幅直接叠前时间偏移方法包括以下步骤：
S101、获取宽方位角地震数据，宽方位角地震数据指示方位各向异性特征；
本实施例中，根据接收到多个地震道的地震数据首先进行方位角统计，在确定地震数据为宽方位角地震数据时，进行方位各向异性偏移成像。也就是说，在确定地震数据具有各向异性特征的情况下才进行后续方位各向异性偏移成像。
针对某一地震道的地震数据，根据该地震道的炮点坐标和检波点坐标，确定该地震道的方位角。示例性的，
某地震道的炮点坐标为(xs,ys)，检波点坐标为(xr,yr)，则根据公式(1)计算地震道的方位角φ：
通过计算各地震道的方位角，绘制地震数据的方位角信息图；根据方位角信息图，确定地震数据是否为宽方位角地震数据。具体的，
根据方位角信息图确定地震数据的方位角的纵横比，若地震数据的方位角的纵横比满足预设纵横比，则确定地震数据为宽方位角地震数据。例如，预设纵横比为1:1，不同应用场景的预设纵横比不同，对此本实施例不作具体限定，满足行业规定值即可。图2为本发明实施例提供的方位角信息图，由图2可知，根据上述方法统计的地震数据的方位角的纵横比基本满足预设纵横比1:1，可以确定该地震数据为宽方位角地震数据。
S102、根据预设方位角度将宽方位角地震数据进行分组，得到多个方位角区间对应的宽方位角地震数据；
在获取宽方位角地震数据之后，根据预设方位角度将宽方位角地震数据进行分组，得到多个方位角区间对应的宽方位角地震数据，示例性的，预设方位角度设置为15°，每15°为一个方位角区间，方位角区间包括：0°-15°，15°-30°，30°-45°等。可以理解，在每个方位角区间内包括多个地震道的地震数据。
S103、通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
具体来说，根据S102的分组，将每个方位角区间对应的宽方位角地震数据进行各向同性叠前时间偏移，得到每个方位角区间对应的走时和振幅值第一偏移成像结果，该第一偏移成像结果并未考虑方向异性对地震波幅值的影响。
S104、根据每个方位角区间对应的第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
本实施例基于不同方位角区间，以各向同性叠前时间偏移处理得到的第一偏移成像结果为基础，进行速度分析与速度场重建，也就是说，以各向同性叠前时间偏移处理所用速度场为背景场进行速度分析与速度场重建，确定每个方位角区间对应的偏移速度场，可以理解，每个方位角区间的偏移速度场不同，通过对每个方位角区间对应的第一偏移成像结果进行速度扫描，可以得到新的各向异性的偏移速度场。
S105、根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，偏移参数包括走时和振幅值；
本实施例的方位各向异性公式包括方位方位各向异性走时公式和方位各向异性幅值公式，具体的，
根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
需要指出的是，偏移参数中的走时和振幅值是在新的各向异性偏移速度场下计算得到的。上述计算过程，以单程波理论和稳相点原理为基础，通过提出新的各向异性等效参数得到方位各向异性介质中的地震波走时和幅值。
下面对方位各向异性介质中地震波走时、幅值的计算原理进行详细说明。
各向异性介质中，P波传播的相速度公式可表示为：
其中，V(θ)表示相速度，θ表示相角度，VP0表示P波垂向速度，VS0表示S波垂向速度，ε和δ表示各项异性参数，由于VS0对相速度的影响几乎可以忽略不计，令VS0为0，则(2)式写成：
相速度与波数-频率的关系可表示为：
cosθ＝V(θ)Ky/ω (5)
将(4)(5)(6)代入(3)，整理得到：
又，
其中，Kx、Ky、Kz为波数K在直角坐标系x、y、z方向上的分量，Kr为波数K在球坐标系上的分量。
将(8)(9)代入(7)，替换公式(7)中的波数K，并开根号，从而得到：
另(10)＝(11)，从而对应得到新的各向异性等效参数：
其中，式(10)至(13)中的φ为方位角。
由深度偏移的相移法可得，在波数－频率域，单个炮点或检波点的地震波场的走时(时间深度)表示为：
其中，ΔTi是各层介质用单程旅行时表达的厚度，n是某一深度包含的层数，是该实际深度的时间深度(单程旅行时)，vnmo(i)是各层介质的各向异性动校正速度，ηi是各层介质的各向异性参数，px是射线参数，j是单位虚数，ω是频率，f(ω)是炮点或检波点的时域信号的傅立叶变换。
将新的各向异性等效参数(12)(13)代入(14)并作空间付氏反变换，得空间－频率域波场，表示为：
式(15)可用稳相点原理求得渐进解为
式(16)中，定义
其中，是的零点，求得代入(18)、(19)得地震波的走时和振幅值。
通过上述计算过程，可以得到每个地震道的走时和幅值。
本实施例中，方位各向异性走时公式具体由式(18)表示，方位各向异性幅值公式具体由式(19)表示。
S106、根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
本实施例中，根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。具体的，
利用波动方程叠前深度偏移的反褶积成像条件保幅成像。设震源是一时间脉冲，则对单个地震道，有成像结果：
其中，F′(t)是f(ω)对应的时域函数的一阶导数，As和ts是炮点至成像点的幅值和走时，Ar和tr是检波点至成像点的幅值和走时。
将S105的每个方位角区间对应的偏移参数分别代入式(20)得到每个方位角区间对应的第二偏移成像结果。
S107、将所有方位角区间对应的第二偏移成像结果进行叠加，得到宽方位角地震数据的保幅偏移成像结果。
基于S106得到的每个方位角区间对应的第二偏移成像结果，将每个方位角区间的第二偏移成像结果在空间上进行叠加，得到宽方位角地震数据的总的保幅偏移成像结果。
现有的真振幅偏移技术只是在地下介质较为简单的情况下(各向同性介质中)可以得到较好的应用，当地质情况较为复杂时，将无法得到较好的保幅成像结果。
本实施例根据成像点处的偏移速度和各向异性等效参数，计算总走时tr+ts和权系数Ar/As，在地震信号的一阶导数上拾取ts+tr时刻的值并乘上权系数Ar/As，得到该地震道在该成像点的偏移幅值，具有较好的保幅效果权，其中权系数Ar/As保证了偏移成像中正确补偿了地震波在非均匀各向异性介质中传播的几何扩散效应。
本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移方法，将获取到的宽方位角地震数据按照预设方位角度进行分组，通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的第一偏移成像结果，并确定每个方位角区间对应的偏移速度场；根据每个偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，根据每个方位角区间对应的偏移参数得到每个方位角区间对应的第二偏移成像结果；将所有方位角区间对应的第二偏移成像结果进行叠加，最终得到宽方位角地震数据的保幅偏移成像结果。由于每个方位角区间对应的偏移成像结果考虑了方位各向异性对地震波走时和幅值的影响，实现各向异性介质中的保幅偏移成像。
图3为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置的结构示意图一，如图3所述，本实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置30，包括：
获取模块31，用于获取宽方位角地震数据，所述宽方位角地震数据指示方位各向异性特征；
分组模块32，用于根据预设方位角度将所述宽方位角地震数据进行分组；得到多个方位角区间对应的宽方位角地震数据；
计算成像模块33，用于通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
确定模块34，用于根据每个方位角区间对应的所述第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
计算模块35，还用于根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，所述偏移参数包括走时和振幅值；
所述计算成像模块33，用于根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
所述计算成像模块33，还用于将所有方位角区间对应的第二偏移成像结果进行叠加，得到所述宽方位角地震数据的保幅偏移成像结果。
可选的，所述计算模块35，具体用于根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
图4为本发明提供的方位各向异性介质中的保幅直接叠前时间偏移装置的结构示意图二，在图3所示装置的基础上，如图4所示，所述保幅直接叠前时间偏移装置30还包括：
图像生成模块36，用于对地震数据进行方位角统计，绘制所述地震数据的方位角信息图；所述地震数据包括所有地震道的炮点和检波点的坐标；
所述确定模块34，还用于根据所述方位角信息图，确定所述地震数据为宽方位角地震数据。
可选的，所述确定模块34，用于根据所述方位角信息图确定所述地震数据的方位角的纵横比；
判断所述地震数据的方位角的纵横比是否满足预设纵横比，若是，则确定所述地震数据为宽方位角地震数据。
可选的，所述计算成像模块33，具体用于根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。
本实施例提供的保幅直接叠前时间偏移装置，可以执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
图5为本发明实施例提供的方位各向异性介质中的保幅直接叠前时间偏移装置的硬件结构图，如图5所示，本实施例提供的保幅直接叠前时间偏移装置50，包括：
存储器51；
处理器52；以及
计算机程序；
其中，计算机程序存储在存储器51中，并被配置为由处理器52执行以实现如前述任一项方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
可选的，存储器51既可以是独立的，也可以跟处理器52集成在一起。
当存储器51是独立于处理器52之外的器件时，连通井组渗流通道的监测装置50还包括：
总线53，用于连接存储器51和处理器52。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器52执行以实现如上方法实施例中保幅直接叠前时间偏移装置所执行的各个步骤。
其中，存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)、随机存储记忆体(Random Access Memory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard DiskDrive，缩写：HDD)或固态硬盘(Solid-State Drive，SSD)等；存储介质还可以包括上述种类的存储器的组合。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种方位各向异性介质中的保幅直接叠前时间偏移方法，其特征在于，包括：
获取宽方位角地震数据，所述宽方位角地震数据指示方位各向异性特征；
根据预设方位角度将所述宽方位角地震数据进行分组，得到多个方位角区间对应的宽方位角地震数据；
通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
根据每个方位角区间对应的所述第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，所述方位各向异性公式包括方位各向异性走时公式和方位各向异性幅值公式，所述偏移参数包括走时和振幅值；
根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
将所有方位角区间对应的第二偏移成像结果进行叠加，得到所述宽方位角地震数据的保幅偏移成像结果。
2.根据权利要求1所述的方法，其特征在于，所述根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，包括：
根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
3.根据权利要求1所述的方法，其特征在于，所述获取宽方位角地震数据，包括：
对地震数据进行方位角统计，绘制所述地震数据的方位角信息图；所述地震数据包括所有地震道的炮点和检波点的坐标；
根据所述方位角信息图，确定所述地震数据为宽方位角地震数据。
4.根据权利要求3所述的方法，其特征在于，所述根据所述方位角信息图，确定所述地震数据为宽方位角地震数据，包括：
根据所述方位角信息图确定所述地震数据的方位角的纵横比；
判断所述地震数据的方位角的纵横比是否满足预设纵横比，若是，则确定所述地震数据为宽方位角地震数据。
5.根据权利要求1所述的方法，其特征在于，所述根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果，包括：
根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。
6.一种方位各向异性介质中的保幅直接叠前时间偏移装置，其特征在于，包括：
获取模块，用于获取宽方位角地震数据，所述宽方位角地震数据指示方位各向异性特征；
分组模块，用于根据预设方位角度将所述宽方位角地震数据进行分组；得到多个方位角区间对应的宽方位角地震数据；
计算成像模块，用于通过各向同性叠前时间偏移方法进行偏移成像，得到多个方位角区间对应的宽方位角地震数据的第一偏移成像结果；
确定模块，用于根据每个方位角区间对应的所述第一偏移成像结果，确定每个方位角区间对应的偏移速度场；
计算模块，还用于根据每个方位角区间对应的偏移速度场和方位各向异性公式，计算每个方位角区间对应的偏移参数，所述方位各向异性公式包括方位各向异性走时公式和方位各向异性幅值公式，所述偏移参数包括走时和振幅值；
所述计算成像模块，用于根据每个方位角区间对应的偏移参数，得到每个方位角区间对应的第二偏移成像结果；
所述计算成像模块，还用于将所有方位角区间对应的第二偏移成像结果进行叠加，得到所述宽方位角地震数据的保幅偏移成像结果。
7.根据权利要求6所述的装置，其特征在于，
所述计算模块，具体用于根据每个方位角区间对应的偏移速度场和方位各向异性走时公式，计算每个方位角区间对应的走时；以及
根据每个方位角区间对应的偏移速度场和方位各向异性幅值公式，计算每个方位角区间对应的振幅值。
8.根据权利要求6所述的装置，其特征在于，所述装置还包括：
图像生成模块，用于对地震数据进行方位角统计，绘制所述地震数据的方位角信息图；所述地震数据包括所有地震道的炮点和检波点的坐标；
所述确定模块，还用于根据所述方位角信息图，确定所述地震数据为宽方位角地震数据。
9.根据权利要求8所述的装置，其特征在于，
所述确定模块，用于根据所述方位角信息图确定所述地震数据的方位角的纵横比；
判断所述地震数据的方位角的纵横比是否满足预设纵横比，若是，则确定所述地震数据为宽方位角地震数据。
10.根据权利要求6所述的装置，其特征在于，
所述计算成像模块，具体用于根据每个方位角区间对应的偏移参数和波动方程叠前深度偏移的反褶积成像条件，得到每个方位角区间对应的第二偏移成像结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810606448.6A CN108802822B (zh) | 2018-06-13 | 2018-06-13 | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810606448.6A CN108802822B (zh) | 2018-06-13 | 2018-06-13 | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108802822A CN108802822A (zh) | 2018-11-13 |
| CN108802822B true CN108802822B (zh) | 2019-05-24 |
Family
ID=64085816
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810606448.6A Expired - Fee Related CN108802822B (zh) | 2018-06-13 | 2018-06-13 | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108802822B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109725354B (zh) * | 2018-11-20 | 2020-07-10 | 中国石油天然气集团有限公司 | 各向异性速度建模方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957455A (zh) * | 2010-09-20 | 2011-01-26 | 中国海洋石油总公司 | 三维保幅叠前时间偏移方法 |
| CN102141633A (zh) * | 2010-12-10 | 2011-08-03 | 中国科学院地质与地球物理研究所 | 各向异性三维叠前时间偏移方法 |
| CN102540250A (zh) * | 2010-12-08 | 2012-07-04 | 同济大学 | 基于方位保真角度域成像的裂缝型油气储层地震探测方法 |
| CN102455436B (zh) * | 2010-11-02 | 2013-10-30 | 中国石油大学(北京) | 有限方位角条件下压噪叠前纵波各向异性裂缝检测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7830747B2 (en) * | 2006-11-01 | 2010-11-09 | Pgs Geophysical As | Method for multi-azimuth prestack time migration for general heterogeneous, anisotropic media |
-
2018
- 2018-06-13 CN CN201810606448.6A patent/CN108802822B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957455A (zh) * | 2010-09-20 | 2011-01-26 | 中国海洋石油总公司 | 三维保幅叠前时间偏移方法 |
| CN102455436B (zh) * | 2010-11-02 | 2013-10-30 | 中国石油大学(北京) | 有限方位角条件下压噪叠前纵波各向异性裂缝检测方法 |
| CN102540250A (zh) * | 2010-12-08 | 2012-07-04 | 同济大学 | 基于方位保真角度域成像的裂缝型油气储层地震探测方法 |
| CN102141633A (zh) * | 2010-12-10 | 2011-08-03 | 中国科学院地质与地球物理研究所 | 各向异性三维叠前时间偏移方法 |
Non-Patent Citations (1)
| Title |
|---|
| 海上多方位地震资料处理关键技术;陈礼;《石油天然气学报（江汉石油学院学报）》;20130731;第35卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108802822A (zh) | 2018-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| Verdon et al. | Microseismic monitoring using a fiber-optic distributed acoustic sensor array | |
| CA2964893C (en) | Structure tensor constrained tomographic velocity analysis | |
| US20120275267A1 (en) | Seismic Data Processing | |
| US9541659B2 (en) | Noise removal from 3D seismic representation | |
| US9733371B2 (en) | Creating seismic images using expanded image gathers | |
| CN106556861B (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
| US10215869B2 (en) | System and method of estimating anisotropy properties of geological formations using a self-adjoint pseudoacoustic wave propagator | |
| US11733417B2 (en) | Quality control and preconditioning of seismic data | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| AU2017366758A1 (en) | Diving wave illumination using migration gathers | |
| CN103852785B (zh) | 地层各向异性的评价方法 | |
| US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| Zhang et al. | Automated microseismic event location by amplitude stacking and semblance | |
| CN108802822B (zh) | 方位各向异性介质中的保幅直接叠前时间偏移方法及装置 | |
| CN109143352A (zh) | 一种各向异性介质地震反射特征方程建立方法 | |
| CN111474580A (zh) | 一种基于炮检距矢量片的方位角道集提取方法和系统 | |
| CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
| CN109425892B (zh) | 地震子波的估计方法及系统 | |
| CN110579798B (zh) | 一种等反射角度间隔的地震采集观测方法及系统 | |
| US11255993B2 (en) | Variable aperture estimation using bottom-up ray tracing | |
| CN110824555B (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
| US20160238727A1 (en) | Removal of an Estimated Acquisition Effect from a Marine Survey Measurement | |
| CN111323819B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190524 Termination date: 20210613 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |