CN113296153B - 轴对称介质各向异性参数确定方法及装置 - Google Patents
轴对称介质各向异性参数确定方法及装置 Download PDFInfo
- Publication number
- CN113296153B CN113296153B CN202010111495.0A CN202010111495A CN113296153B CN 113296153 B CN113296153 B CN 113296153B CN 202010111495 A CN202010111495 A CN 202010111495A CN 113296153 B CN113296153 B CN 113296153B
- Authority
- CN
- China
- Prior art keywords
- longitudinal wave
- wave speed
- parameter
- initial
- speed adjusting
- 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
Links
Images
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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/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/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
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
技术领域
本发明涉及地震资料处理技术领域,尤其涉及一种轴对称介质各向异性参数确定方法及装置。
背景技术
传播介质的弹性特征随方向变化而变化,致使地震波传播的速度值与传播方向有关,称这种物理现象为“各向异性”。随着新的地震资料采集技术和处理算法的进步,地震各向异性对地震数据的定量处理和解释产生巨大的影响。在许多地震数据处理实例中,不考虑各向异性的影响,很可能会导致错误的处理和解释。因此,对于存在各向异性的复杂地质区域,根据地震波数据反演地下介质的速度和各向异性参数,从而建立地下介质的速度和各向异性模型成为当前油气勘探与开发的一项重要的任务。
具有垂向对称轴的横向各向同性(VTI)介质是一种最简单、最普遍的各向异性介质。在均匀的轴对称介质中,地震波沿水平方向传播时体现了横向各向同性,沿其它方向传播时体现各向异性,即速度与方向有关。并且随着各向异性程度的增加,地震波的振幅和走时与各向同性介质中的相比有更大的差别。通过引入Thomsen各向异性参数,能够更简单直观地描述轴对称介质中速度和各向异性参数估算的问题,不仅扩展了传统的勘探地震数据处理技术,也是油气勘探的需要。
根据广义胡克定律,从弹性模量矩阵来说,某一方向上各向同性介质共有5个独立常数,即c11,c13,c33,c44和c66。轴对称介质弹性模量矩阵如下:
其中,vp表示纵波垂直方向传播速度;ε表示纵波各向异性参数,描述了P波水平速度和垂向速度的微小差别;δ表示纵波速度调节参数,影响垂直介质对称轴方向的纵波速度大小;ρ表示地下介质密度。
因此,只要得到vp、δ和ε这3个参数,就可以采用各向异性介质的弯曲射线追踪方法计算轴对称介质叠前深度偏移的旅行时,从而进行轴对称介质叠前深度偏移。
目前,对轴对称介质的叠前深度偏移,已经总结出高效准确的各向异性射线追踪方法,但对于各向异性参数的准确估计,尚无比较成熟的技术措施。对于轴对称介质来说,vp、δ和ε这3个参数互相耦合,单纯利用地震资料不能同时反演得到,且存在着非唯一性。目前工业界常用的做法是:通过测井或VSP约束层析反演等方法迭代得到轴对称介质的各项同性速度vp0,使得叠前深度偏移成像效果最佳,通过公式(5)井震地层厚度对比确定δ,利用公式(6)确定各向异性速度vp,在已知vp和δ的情况下,由地震层析反演确定ε。
其中,δ表示纵波速度调节参数,h表示地震地层厚度,h0表示实际地层厚度,vp0表示各向同性速度,vp表示各向异性速度。
但实际应用中,经常遇到工区内井数量比较少、地震资料的信噪较低等不利因素,导致各向异性参数估计精度不高,主要表现在以下两个方面:
各向同性速度求取精度不够,各向同性偏移剖面成像深度偏差较大,这种误差在利用公式(5)计算后将传递给δ,导致δ值异常,最终得到的δ值不符合实际岩石物理规律;
利用测井或VSP速度强约束各向同性速度场vp0,得到各向异性速度场vp,如果不更新各向异性参数场δ进行叠前深度偏移,将出现共成像点道集中同相轴不平、成像质量下降的情况;其中,共成像点道集中同相轴不平反映了轴对称介质各向异性参数不准确;如果更新各向异性参数场δi,那么新生成的δi与原有的δ0如何合成应用,目前尚无有效技术支撑。
因此,现有技术中求取的轴对称介质各向异性参数准确度不高,不符合实际岩石物理规律,且以此轴对称介质各向异性参数进行的叠前深度偏移,会出现共成像点道集中同相轴不平、成像质量下降的情况。
发明内容
本发明实施例提供一种轴对称介质各向异性参数确定方法,用以提高求取的轴对称介质各向异性参数的准确度,进而保证叠前深度偏移的共成像点道集中同相轴拉平,提高叠前深度偏移的成像质量,该方法包括:
获取轴对称介质的地震资料,根据所述地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
根据所述初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
根据所述纵波速度调节参数和所述各向同性纵波速度,得到初始各向异性纵波速度;
循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据所述初始各向异性纵波速度和所述纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数;
将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;
根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数,包括:
根据所述新的地震地层厚度和所述实际地层厚度,得到井震误差;
根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值;
根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数;
按照如下公式,根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差;h0表示实际地层厚度;
按照如下公式,根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数;δ0表示初始纵波速度调节参数。
本发明实施例还提供一种轴对称介质各向异性参数确定装置,用以提高求取的轴对称介质各向异性参数的准确度,进而保证叠前深度偏移的共成像点道集中同相轴拉平,提高叠前深度偏移的成像质量,该装置包括:
初始数据获取模块,用于获取轴对称介质的地震资料,根据所述地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
纵波速度调节参数计算模块,用于根据所述初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
初始各向异性纵波速度计算模块,用于根据所述纵波速度调节参数和所述各向同性纵波速度,得到初始各向异性纵波速度;
轴对称介质各向异性参数确定模块,用于循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据所述初始各向异性纵波速度和所述纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数;
将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;
所述轴对称介质各向异性参数确定模块具体用于:
根据所述新的地震地层厚度和所述实际地层厚度,得到井震误差;
根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值;
根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数;
按照如下公式,根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差;
按照如下公式,根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述轴对称介质各向异性参数确定方法。
本发明实施例也提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述轴对称介质各向异性参数确定方法的计算机程序。
本发明实施例中,通过获取轴对称介质的地震资料,根据地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;根据初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度;循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:根据初始各向异性纵波速度和纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据新的地震地层厚度和实际地层厚度得到新的纵波速度调节参数;将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;根据岩石物理允许范围,得到纵波速度调节参数,使轴对称介质各向异性参数符合实际岩石物理规律,提高了求取的轴对称介质各向异性参数的准确度;通过多次重复循环过程,使得最终的各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,从而提高叠前深度偏移的成像质量。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中轴对称介质各向异性参数确定方法的示意图。
图2-图7为本发明一具体应用实施中不同各向异性参数下的速度场及偏移剖面示意图。
图8为本发明实施例中轴对称介质各向异性参数确定装置的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了解决轴对称介质各向异性参数准确度不高的问题,本发明实施例提供了一种轴对称介质各向异性参数确定方法,用以提高求取的轴对称介质各向异性参数的准确度,进而保证叠前深度偏移的共成像点道集中同相轴拉平,提高叠前深度偏移的成像质量,如图1所示,该方法包括:
步骤101:获取轴对称介质的地震资料,根据地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
步骤102:根据初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
步骤103:根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度;
步骤104:循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
步骤105:根据初始各向异性纵波速度和纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据新的地震地层厚度和实际地层厚度得到新的纵波速度调节参数;
步骤106:将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数。
由图1所示流程可以得知,本发明实施例中,通过根据岩石物理允许范围,得到纵波速度调节参数,使轴对称介质各向异性参数符合实际岩石物理规律,提高了求取的轴对称介质各向异性参数的准确度;通过多次重复循环过程,使得最终的各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,从而提高叠前深度偏移的成像质量。
具体实施时,首先获取轴对称介质的地震资料,根据地震资料,反演迭代求取轴对称介质的各向同性纵波速度,并根据得到的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度。
其次,根据初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数。具体过程,包括:
根据初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数;
根据岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数。
其中,根据初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数,具体包括:
获取测井资料,根据获取的测井资料,计算得到实际地层厚度;
根据实际地层厚度和初始地震地层厚度,求取初始纵波速度调节参数。
具体实施时,按照如下公式,根据实际地层厚度和初始地震地层厚度,求取初始纵波速度调节参数:
其中,δ0表示初始纵波速度调节参数;h表示初始地震地层厚度;h0表示实际地层厚度。
上述岩石物理允许范围是指实际岩石物理规律下所允许的纵波速度调节参数范围,可根据实际各向异性参数的统计结果总结得到,例如,Thomsen(1986)曾经总结了多种常见的沉积岩的各向异性参数统计结果,归纳了截至当时所有公开测得的沉积岩各向异性参数,其研究结果表明大多数岩石属于小到中等程度各向异性,其|δ|<0.2。表1给出了一组VTI类型的岩石对应的Thomsen参数。
表1 VTI类型岩石的Thomsen参数
Thomsen参数 | ε | δ | γ |
砂岩 | 0.110 | -0.035 | 0.255 |
泥岩 | 0.034 | 0.211 | 0.046 |
疏松砂岩 | 0.097 | 0.091 | 0.051 |
泥质页岩 | 0.189 | 0.204 | 0.175 |
页岩 | 0.255 | -0.05 | 0.480 |
具体实施例中,根据岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数,包括:
判断初始纵波速度调节参数是否超出岩石物理允许范围,若超出岩石物理允许范围,调整初始纵波速度调节参数到岩石物理允许范围内,得到纵波速度调节参数;
若未超出岩石物理允许范围,将初始纵波速度调节参数作为纵波速度调节参数。
得到纵波速度调节参数后,根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度。具体实施时,按照如下公式进行计算:
其中,vp表示初始各向异性纵波速度;vp0表示各向同性纵波速度;δ表示纵波速度调节参数。
得到初始各向异性纵波速度后,循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平且井震误差符合预设标准,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据初始各向异性纵波速度和纵波速度调节参数,反演迭代得到新的各向异性纵波速度,根据纵波速度调节参数和反演得到的新的各向异性纵波速度进行各向异性叠前深度偏移得到新的地震地层厚度,根据新的地震地层厚度和实际地层厚度得到纵波速度调节参数的校正值;
根据纵波速度调节参数的校正值得到纵波速度调节参数,将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数。
具体实施时,根据新的地震地层厚度和实际地层厚度得到新的纵波速度调节参数包括:
根据新的地震地层厚度和实际地层厚度,得到井震误差;
根据井震误差和实际地层厚度,计算得到纵波速度调节参数的校正值;
根据纵波速度调节参数和纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数。
其中,按照如下公式,根据井震误差和实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差;Δh=h′-h0,h′表示新的地震地层厚度;i表示反演迭代的重复次数。
推导过程如下:
令井震误差(通常为正)Δh带入式(7)有:
按照如下公式,根据纵波速度调节参数和纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数。
可以理解的是,上述计算公式仅为举例,实施时可以对该公式进行变形,或采用其它公式或方法,本领域技术人员可以理解,变形后的这些公式或方法均落入本发明的保护范围,实施例中不再赘述。
下面给出一具体实例说明本发明实施例如何确定轴对称介质各向异性参数。
如图2所示,为通过层析反演迭代求取最优成像的各向同性速度场及各向同性偏移剖面,可以看到共成像点偏移道集基本拉平,成像效果较为理想,但由于资料信噪比相对较低,这种道集拉平程度的判定存在人为干扰,所求取得各项同性速度场必然存在一定误差,而这部分误差将会累计到井震深度误差中。
图3为根据测井与偏移剖面深度差通过式(5)、(6)求取的各向异性速度场Vp及各向异性参数场δ0,可以看到,虽然偏移成像质量未发生改变,井震深度误差减小,但各向异性速度场与测井速度差异较大,同时δ场存在明显异常,超出正常的岩石物理允许范围,说明这种井震误差并不是完全由各向异性参数δ0引起,一定程度上还累加了各向同性速度不准导致的误差。
图4为利用测井速度对各向同性速度场进行强约束得到的各向异性速度场和偏移剖面,此时各向异性参数场δ人为给定区域经验常数(δ=0.05),可以看到,虽然各向异性速度与测井速度吻合度较高,井震深度误差减小,但是出现共成像点偏移道集不平,偏移剖面成像质量明显下降的情况,说明这种井震速度误差也不完全是由各向异性速度的误差造成的,一定程度上还受各向异性参数估算精度不足的影响。
图5为通过各向异性速度场Vp和初始人为设置为合理的各向异性参数场δ(|δ|<0.2)进行偏移的结果,可以看到井震深度误差减小,偏移成像质量基本保持略有下降。
图6为通过本发明实施例中提供的轴对称介质各向异性参数确定方法,更新迭代图5中的各向异性速度场Vp和各向异性场δ后的偏移结果,可以看到利用本发明实施例中的更新迭代求取方法,在保证道集拉平的同时,井震深度误差减小,与钻井曲线吻合度较高,各向异性速度场Vp合理,并且各向异性参数场δ估值在岩石物理允许范围内。
四川盆地南缘外高山区-云贵高原YS工区,地震资料信噪比低,页岩储层裂缝发育,各向异性复杂。通过本发明实施例中的轴对称介质各向异性参数确定方法的应用,最终各向异性速度场Vp合理,各向异性参数场δ估值在岩石物理允许范围内,井震误差较小,较各向同性偏移剖面成像质量有明显提高,效果对比图如图7所示。
上述具体应用的实施仅为举例,其余实施方式不再一一赘述。
基于同一发明构思,本发明实施例还提供一种轴对称介质各向异性参数确定装置,由于轴对称介质各向异性参数确定装置所解决问题的原理与轴对称介质各向异性参数确定方法相似,因此轴对称介质各向异性参数确定装置的实施可以参见轴对称介质各向异性参数确定方法的实施,重复之处不再赘述,具体结构如图8所示:
初始数据获取模块801,用于获取轴对称介质的地震资料,根据地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
纵波速度调节参数计算模块802,用于根据初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
初始各向异性纵波速度计算模块803,用于根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度;
轴对称介质各向异性参数确定模块804,用于循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据初始各向异性纵波速度和纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据新的地震地层厚度和实际地层厚度得到新的纵波速度调节参数;
将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数。
具体实施例中,纵波速度调节参数计算模块802包括:
初始参数求取单元,用于根据初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数;
参数调整单元,用于根据岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数。
其中,初始参数求取单元具体用于:
根据获取的测井资料,计算得到实际地层厚度;
根据实际地层厚度和初始地震地层厚度,求取初始纵波速度调节参数。
具体地,按照如下公式,根据实际地层厚度和初始地震地层厚度,求取初始纵波速度调节参数:
其中,δ0表示初始纵波速度调节参数;h表示初始地震地层厚度;h0表示实际地层厚度。
参数调整单元具体用于:
判断初始纵波速度调节参数是否超出岩石物理允许范围,若超出岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数;
若未超出岩石物理允许范围,将初始纵波速度调节参数作为纵波速度调节参数。
具体实施时,初始各向异性纵波速度计算模块803具体用于:按照如下公式,根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度:
其中,vp表示初始各向异性纵波速度;vp0表示各向同性纵波速度;δ表示纵波速度调节参数。
具体实施例中,轴对称介质各向异性参数确定模块804具体用于:
根据新的地震地层厚度和实际地层厚度,得到井震误差;
根据井震误差和实际地层厚度,计算得到纵波速度调节参数的校正值;
根据纵波速度调节参数和纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数。
具体地,按照如下公式,根据井震误差和实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差。
具体地,根据纵波速度调节参数和纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述轴对称介质各向异性参数确定方法。
本发明实施例还提供一种计算机可读存储介质,存储有执行上述轴对称介质各向异性参数确定方法的计算机程序。
综上所述,本发明实施例提供的轴对称介质各向异性参数确定方法及装置具有如下优点:
通过获取轴对称介质的地震资料,根据地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;根据初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;根据纵波速度调节参数和各向同性纵波速度,得到初始各向异性纵波速度;循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:根据初始各向异性纵波速度和纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据新的地震地层厚度和实际地层厚度得到新的纵波速度调节参数;将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;根据岩石物理允许范围,得到纵波速度调节参数,使轴对称介质各向异性参数符合实际岩石物理规律,提高了求取的轴对称介质各向异性参数的准确度;通过多次重复循环过程,使得最终的各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平且井震误差符合预设标准,从而提高叠前深度偏移的成像质量。
本领域内的技术人员应明白,本发明的实施例可提供为方法、装置、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种轴对称介质各向异性参数确定方法,其特征在于,包括:
获取轴对称介质的地震资料,根据所述地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
根据所述初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
根据所述纵波速度调节参数和所述各向同性纵波速度,得到初始各向异性纵波速度;
循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据所述初始各向异性纵波速度和所述纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数;
将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;
根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数,包括:
根据所述新的地震地层厚度和所述实际地层厚度,得到井震误差;
根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值;
根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数;
按照如下公式,根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差;h0表示实际地层厚度;
按照如下公式,根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数;δ0表示初始纵波速度调节参数。
2.如权利要求1所述的方法,其特征在于,根据所述初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数,包括:
根据所述初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数;
根据所述岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数。
3.如权利要求2所述的方法,其特征在于,根据所述初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数,包括:
根据获取的测井资料,计算得到实际地层厚度;
根据所述实际地层厚度和所述初始地震地层厚度,求取初始纵波速度调节参数。
5.如权利要求2所述的方法,其特征在于,根据岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数,包括:
判断所述初始纵波速度调节参数是否超出岩石物理允许范围,若超出岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数;
若未超出岩石物理允许范围,将所述初始纵波速度调节参数作为纵波速度调节参数。
7.一种轴对称介质各向异性参数确定装置,其特征在于,包括:
初始数据获取模块,用于获取轴对称介质的地震资料,根据所述地震资料,反演迭代求取轴对称介质的各向同性纵波速度,进行叠前深度偏移得到初始地震地层厚度;
纵波速度调节参数计算模块,用于根据所述初始地震地层厚度、实际地层厚度和岩石物理允许范围,得到纵波速度调节参数;
初始各向异性纵波速度计算模块,用于根据所述纵波速度调节参数和所述各向同性纵波速度,得到初始各向异性纵波速度;
轴对称介质各向异性参数确定模块,用于循环执行以下步骤,直至各向异性叠前深度偏移的共成像点偏移道集中同相轴拉平,确定新的纵波速度调节参数为轴对称介质各向异性参数:
根据所述初始各向异性纵波速度和所述纵波速度调节参数,反演迭代得到新的各向异性纵波速度,进行各向异性叠前深度偏移得到新的地震地层厚度,根据所述新的地震地层厚度和所述实际地层厚度得到新的纵波速度调节参数;
将新的纵波速度调节参数和新的各向异性纵波速度作为下一次循环的初始各向异性纵波速度和纵波速度调节参数;
所述轴对称介质各向异性参数确定模块具体用于:
根据所述新的地震地层厚度和所述实际地层厚度,得到井震误差;
根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值;
根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数;
按照如下公式,根据所述井震误差和所述实际地层厚度,计算得到纵波速度调节参数的校正值:
其中,δi表示纵波速度调节参数的校正值;i表示反演迭代的重复次数;Δh表示井震误差;h0表示实际地层厚度;
按照如下公式,根据所述纵波速度调节参数和所述纵波速度调节参数的校正值,累加求和得到新的纵波速度调节参数:
其中,δ表示新的纵波速度调节参数;n表示反演迭代的重复总次数;δ0表示初始纵波速度调节参数。
8.如权利要求7所述的装置,其特征在于,所述纵波速度调节参数计算模块包括:
初始参数求取单元,用于根据所述初始地震地层厚度和实际地层厚度,求取初始纵波速度调节参数;
参数调整单元,用于根据实际岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数。
9.如权利要求8所述的装置,其特征在于,所述初始参数求取单元具体用于:
根据获取的测井资料,计算得到实际地层厚度;
根据所述实际地层厚度和所述初始地震地层厚度,求取初始纵波速度调节参数。
11.如权利要求8所述的装置,其特征在于,所述参数调整单元具体用于:
判断所述初始纵波速度调节参数是否超出岩石物理允许范围,若超出岩石物理允许范围,调整初始纵波速度调节参数,得到纵波速度调节参数;
若未超出岩石物理允许范围,将所述初始纵波速度调节参数作为纵波速度调节参数。
13.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010111495.0A CN113296153B (zh) | 2020-02-24 | 2020-02-24 | 轴对称介质各向异性参数确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010111495.0A CN113296153B (zh) | 2020-02-24 | 2020-02-24 | 轴对称介质各向异性参数确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113296153A CN113296153A (zh) | 2021-08-24 |
CN113296153B true CN113296153B (zh) | 2023-04-25 |
Family
ID=77317812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010111495.0A Active CN113296153B (zh) | 2020-02-24 | 2020-02-24 | 轴对称介质各向异性参数确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113296153B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100135115A1 (en) * | 2008-12-03 | 2010-06-03 | Chevron U.S.A. Inc. | Multiple anisotropic parameter inversion for a tti earth model |
MX2018002023A (es) * | 2015-09-18 | 2018-04-13 | Halliburton Energy Services Inc | Estimacion basada en inversion global de parametros de anisotropia para medios ortorrombicos. |
CA3008206C (en) * | 2016-01-15 | 2022-08-02 | Landmark Graphics Corporation | Semblance-based anisotropy parameter estimation using isotropic depth-migrated common image gathers |
CN109725354B (zh) * | 2018-11-20 | 2020-07-10 | 中国石油天然气集团有限公司 | 各向异性速度建模方法及系统 |
CN109696704B (zh) * | 2019-01-28 | 2020-05-01 | 中国海洋石油集团有限公司 | 一种基于纵波阻抗约束的地震各向异性δ建模方法 |
-
2020
- 2020-02-24 CN CN202010111495.0A patent/CN113296153B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113296153A (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN106443774A (zh) | 一种提高起伏地表叠前深度偏移成像精度的方法 | |
CN104459782A (zh) | 使用薄层的沿层层析反演地震速度建模方法以及建模单元 | |
CN103616723A (zh) | 基于avo特征的crp道集真振幅恢复方法 | |
CN109946742A (zh) | 一种TTI介质中纯qP波地震数据模拟方法 | |
CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
CN105626056B (zh) | 一种确定目标区三维地层孔隙压力的方法和装置 | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
CN113296153B (zh) | 轴对称介质各向异性参数确定方法及装置 | |
CN106932823A (zh) | 一种光滑基准面选择技术 | |
CN111435174B (zh) | 强反射地区地震资料振幅补偿方法及装置 | |
CN111399037B (zh) | 高速顶界面提取的方法和装置 | |
CN111381280A (zh) | 预测储层含烃饱和度的方法及装置 | |
CN115586573A (zh) | 一种致密砂岩储层的动态约束物性参数地震反演方法 | |
CN113138412A (zh) | 深层页岩气孔隙度地震预测方法及装置 | |
CN111624662B (zh) | 速度场校正方法、装置和设备 | |
CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
CN112444848B (zh) | 一种全波形反演的方法及系统 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
NO20190217A1 (en) | Correcting biases in microseismic-event data | |
CN111198401B (zh) | 基于多地球物理参数相关性的测井曲线校正方法与装置 | |
CN113156494B (zh) | 面元等效覆盖次数的计算方法及装置 | |
CN116931076A (zh) | 一种全节点地震资料多约束深度域速度建模方法 | |
CN113267808B (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 |