CN110658558A - 吸收衰减介质叠前深度逆时偏移成像方法及系统 - Google Patents
吸收衰减介质叠前深度逆时偏移成像方法及系统 Download PDFInfo
- Publication number
- CN110658558A CN110658558A CN201910912183.7A CN201910912183A CN110658558A CN 110658558 A CN110658558 A CN 110658558A CN 201910912183 A CN201910912183 A CN 201910912183A CN 110658558 A CN110658558 A CN 110658558A
- Authority
- CN
- China
- Prior art keywords
- wave field
- medium
- time
- imaging
- value
- 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.)
- Pending
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/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/51—Migration
- G01V2210/512—Pre-stack
Abstract
一种吸收衰减介质叠前深度逆时偏移成像方法及系统,方法包括:步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值;步骤2:基于L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;步骤3:根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;步骤5:获得黏声介质波场正传算子和黏声介质波场反传算子,根据正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。本发明在保证成像精度条件下大幅提高叠前深度逆时偏移成像的效率。
Description
技术领域
本发明涉及地震波场偏移成像技术领域,具体涉及一种吸收衰减介质叠前深度逆时偏移成像方法及系统。
背景技术
地震波在地下介质中传播时,介质的非完全弹性会对地震波场造成吸收衰减作用。这种介质本身固有的吸收特性通常用品质因子来描述,它与介质内部的结构特征以及饱和度、孔隙度、渗透率等因素密切相关。利用地层吸收性质可以用来预测岩性、砂泥岩分布,在有利条件下还可以直接用来预测石油和天然气的存在。因此,地震波的吸收特性分析研究,在油气、水资源勘探等工程应用领域,具有十分重要的意义。
另外,随着油气勘探的深入及勘探目标复杂程度的增加,构造成像已不能满足勘探的需求,对岩性成像的需求越来越迫切。在保证成像位置准确的情况下,为了获得高精度保幅成像,需要校正黏滞性对成像的影响,实现真振幅成像。大部分真振幅成像算法是建立在完全弹性介质条件的假设下,然而实际地下介质,特别是在近地表及油气储层,多为一种黏弹性介质。地震波在黏弹性介质中的传播主要表现为速度频散与振幅衰减。不考虑黏滞性的叠前成像算法不仅会使成像位置发生偏离,而且还会引起成像振幅的欠估计,严重影响甚至误导随后的地震数据处理、解释等工作。
对由黏滞性引起的地震波吸收衰减的补偿,常规的方法是进行反Q滤波以增强地震波特别是高频成分的能量,达到提高分辨率的效果。但常规的反Q滤波方法大部分基于层状介质假设,虽能够对能量进行一定补偿,但在复杂介质条件下,该方法并不符合地震波的传播规律,不是一种精确的补偿方法;另一种方法是在成像的过程中对地震波进行能量补偿及相位校正,即反Q偏移方法。目前已有的反Q偏移方法主要有三类,即基于射线理论、基于单程波方程与基于双程波方程的反Q偏移方法。基于射线理论的成像方法由于基于高频近似,较难处理多波至问题,在复杂介质模型中的应用受到限制;基于波动理论的单程波反Q偏移方法对陡倾角构造成像有限制,不能适应复杂地质构造。近年来许多学者将逆时偏移扩展到黏声介质中,通过层析反演技术建立精确的Q模型,并通过Q-RTM(逆时偏移成像)技术实现地层的非弹性吸收衰减补偿,这也是近年来高分辨率地震勘探的重要尝试。
由于Q值是关于频率的函数,传统Q-RTM的实现一般在频率域进行,其特点是求解容易,精度较高,然而高维傅里叶正反变换使得该方法计算效率非常低,不能进行大规模应用。因此,期待开发一种能够保证成像精度且效率更高的吸收衰减介质逆时偏移成像方法,以使该算法能适应于大规模数据应用。
发明内容
本发明一方面提出了一种吸收衰减介质叠前深度逆时偏移成像方法,包括:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据所述拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据所述正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
优选地,根据以下公式(1)计算所述积分函数值:
其中,F表示积分函数,w为当前计算频率。
优选地,所述拟合参数表达式为:
其中,wa、wb为有效频带范围的边界值,Q0表示当前Q值。
优选地,所述时间空间域波场正传控制方程和时间空间域波场反传控制方程分别如公式(8)和(9)所示:
其中,P为波场,v为介质速度,t为旅行时,Ψ为辅助变量,s为震源。
优选地,利用高阶有限差分算法在时间空间域分别求解所述时间空间域正传控制方程和时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子。
优选地,所述成像条件如公式(10)所示:
其中,Ic(x)表示吸收衰减介质中归一化互相关成像结果,Sc(t,x),Rc(t,x)分别表示t时刻空间位置x处的炮点正传波场与检波点反传波场,T为波场延拓最大时间。
本发明另一方面提供一种吸收衰减介质叠前深度逆时偏移成像方法系统,所述系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据所述拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据所述正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
优选地,根据以下公式(1)计算所述积分函数值:
其中,F表示积分函数,w为当前计算频率;
所述拟合参数表达式为:
其中,wa、wb为有效频带范围的边界值,Q0表示当前Q值。
优选地,所述时间空间域波场正传控制方程和时间空间域波场反传控制方程分别如公式(8)和(9)所示:
其中,P为波场,v为介质速度,t为旅行时,Ψ为辅助变量,s为震源。
优选地,所述成像条件如公式(10)所示:
其中,Ic(x)表示吸收衰减介质中归一化互相关成像结果,Sc(t,x),Rc(t,x)分别表示t时刻空间位置x处的炮点正传波场与检波点反传波场,T为波场延拓最大时间。
本发明的有益效果在于,基于衰减介质标准线性体模型,实现了对Q值的多参数高精度拟合,基于拟合参数从频散关系出发在时间空间域重新推导吸收衰减介质条件下波场正传控制方程与带吸收衰减补偿的波场反传控制方程,新的控制方程不再含有拟微分算子项,因此可以利用时间空间域高阶有限差分算法进行数值求解,能够在保证成像精度条件下大幅提高吸收衰减介质Q-RTM成像算法的计算效率,使得该算法能适应大规模实际数据应用。本方法还可利用CPU/GPU协同计算进一步提高计算效率。
本发明具有其它特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的附图标记通常代表相同部件。
图1显示根据本发明示例性实施例的吸收衰减介质叠前深度逆时偏移成像方法的流程图;
图2a和图2b分别显示本发明示例性实施例中的异常体二维速度模型和理论Q值模型;
图3a和图3b分别显示本发明示例性实施例中的某时刻二维声波介质波场快照与黏声介质波场快照;
图4显示图3a中的第700道声波波形和图3b中的黏声介质波形对比;
图5a和图5b分别显示对图3a的二维声波介质进行正演的单炮和对图3b的黏声介质进行正演的单炮;
图6和图7分别显示吸收衰减介质常规RTM结果和Q-RTM偏移结果;
图8a和图8b分别显示本发明示例性实施例中的三维速度模型和理论Q值模型;
图9a和图9b分别显示本发明示例性实施例中的三维声波介质波场快照和三维黏声介质波场快照;
图10a和图10b分别是图9a和图9b的内部剖视图;
图11a和图11b分别是图9a和图9b的内部切片图;
图12a和图12b分别显示对图9a和图9b所示波场进行正演的单炮;
图13显示常规RTM的成像结果;
图14显示根据本发明的吸收衰减介质叠前深度逆时偏移成像方法的成像结果;
图15a和图15b分别显示常规RTM成像对应的频谱和根据本发明的吸收衰减介质叠前深度逆时偏移成像方法对应的频谱。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
本发明一方面提供一种吸收衰减介质叠前深度逆时偏移成像方法,如图1所示,该方法包括:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述l个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
实际地下介质为非完全弹性介质,即为黏弹性介质。地震波在黏弹性介质中传播存在的吸收衰减作用主要表现为速度频散与振幅衰减,这种介质本身固有的吸收特性通常用品质因子Q来描述。为了更真实地描述波场的传播规律,进行地震波的正演与成像的时候,需要考虑地下黏性介质对波场的影响。传统吸收衰减介质Q-RTM方法由于其波场传播方程中含有拟微分算子,不能进行时间空间域的有限差分高效求解,因此在频率域实现,然而其高维傅里叶正反变换使得该计算方法效率非常低,不适用于三维大规模计算,在实际数据应用中存在瓶颈。
本发明实施例基于衰减介质标准线性体模型,实现了对Q值的多参数高精度拟合,基于拟合参数从频散关系出发在时间空间域重新推导吸收衰减介质条件下波场正传控制方程与带吸收衰减补偿的波场反传控制方程,新的控制方程不再含有拟微分算子项,因此可以利用时间空间域高阶有限差分算法进行数值求解,能够在保证成像精度条件下大幅提高吸收衰减介质Q-RTM成像算法的计算效率,使得该算法能适应大规模实际数据应用。。
在常规地震数据的频带范围内,Q值近似不随频率变化,基于广义标准线性体模型(SLS),通过引入辅助变量对常Q(不随频率变化)模型进行拟合,推导能利用有限差分法进行高效求解的非拟微分波场传播控制方程是提高黏声介质正演模拟效率的有效途径。
以下具体描述根据本发明实施例的基于黏声方程的正演模拟方法的步骤。
在步骤1中,在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L。在实际应用中,可根据经验值给定L个应变拟合参数τσl的初始值。
在步骤2中,在有效频带范围内,基于L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值。
其中,积分函数F如公式(1)所示:
其中,F表示积分函数,w为当前计算频率。
在步骤3中,基于积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ。
具体地,在地震数据的有效频带范围内,依据给定的Q模型构建拟合目标函数,通过使拟合目标函数的值最小,可得到Q值的L个广义标准线性体拟合后的优化参数,其中拟合目标函数为:
其中,J为拟合目标函数,wa、wb为有效频带范围的边界值,w为当前计算频率,Q0和Q分别为当前Q值以及实际拟合得到的Q值。
使拟合目标函数的值最小,即令拟合目标函数的导数为零,如以下公式(3)所示:
已知实际拟合得到的Q值Q与应变拟合参数τσl、应力拟合参数τ存在如下关系:
将公式(1)和公式(4)代入公式(3)中,即有:
综合上述,得到应力拟合参数τ的表达式为:
在地震数据有效频带范围内,给定L个应变拟合参数τσl的初始值,通过(6)式实现了将Q值用多参数的Taup值进行表征,便于从频率波数域频散关系推导出新的波场控制方程。
在步骤4中,针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体,拟合参数数据体包括对应于每个Q值的应变拟合参数τσl(l=1…L)和应力拟合参数τ。
在步骤5中,根据拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场。
黏声介质频散关系描述了波的传播规律,其表达式如公式(7)所示:
其中,w为频率,V为相速度,c0,w0为给定的初始速度与频率,cp为实际介质速度,γ为与衰减因子相关的参数,其表达式为γ=arctan(1/Q)/π。
基于拟合参数数据体,推导公式(7)对应的时间空间域波场正传控制方程,如公式(8)所示:
其中,p为波场,v为介质速度,t为旅行时,L为拟合标准线性体的个数,ψ为辅助变量,s为震源。
由吸收衰减介质补偿原理可知,Q-RTM对应的波场反传算子为带吸收衰减补偿的稳定逆时传播算子,即保证振幅的补偿同时保持相位频散关系不变,由频散关系推导其反传控制方程为:
利用拟合后的多参数表征Q值,得到的时间空间域波场正传控制方程(8)和反传控制方程(9)与常规声波波动方程逆时偏移控制方程类似,不含有拟微分项,可以利用高阶有限差分算法在时空域进行高效求解计算。利用高阶有限差分算法在时空域求解时间空间域正传控制方程和时间空间域反传控制方程,即获得黏声介质波场正传算子和黏声介质波场反传算子,根据正传算子和反传算子进行波场的正演模拟,即获得炮点正传波场与检波点反传波场。
在步骤6中,利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
具体地,在常密度吸收衰减介质中,逆时偏移通过求解波动方程实现波场外推,其炮点正传波场与检波点带吸收衰减补偿的反传波场进行零延迟互相关成像,如公式(10)所示:
其中,Ic(x)表示吸收衰减介质中归一化互相关成像结果,Sc(t,x),Rc(t,x)分别表示t时刻空间位置x=(x,y,z)处带补偿的炮点正传波场与检波点反传波场,T为波场延拓最大时间。
成像条件公式(10)实现了成像结果振幅的补偿与相位频散校正,消除了地震波吸收衰减对成像结果的影响,改善成像分辨率,提高整体成像质量。
本发明实施例还提供一种吸收衰减介质叠前深度逆时偏移成像系统,所述系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据所述拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据所述正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
优选地,根据以下公式(1)计算所述积分函数值:
其中,F表示积分函数,w为当前计算频率;
所述拟合参数表达式为:
其中,wa、wb为有效频带范围的边界值,Q0表示当前Q值。
优选地,所述时间空间域波场正传控制方程和时间空间域波场反传控制方程分别如公式(8)和(9)所示:
其中,P为波场,v为介质速度,t为旅行时,Ψ为辅助变量,s为震源。
优选地,所述成像条件如公式(10)所示:
其中,Ic(x)表示吸收衰减介质中归一化互相关成像结果,Sc(t,x),Rc(t,x)分别表示t时刻空间位置x处的炮点正传波场与检波点反传波场,T为波场延拓最大时间。
实施例
以下参考附图描述根据本发明实施例的吸收衰减介质叠前深度逆时偏移成像方法的效果。
为了验证根据本发明实施例的基于黏声方程的正演模拟方法的正演效果,以异常体二维速度模型和复杂三维速度模型进行测试。图2a和图2b分别显示本发明示例性实施例中的异常体二维速度模型和理论Q值模型。图3a和图3b分别显示本发明示例性实施例中某时刻的二维声波介质波场快照与黏声介质波场快照,从图3b中可见黏声介质波场的振幅能量衰减,相位发生频散。图4显示图3a中的第700道声波波形和图3b中的黏声介质波形,其中带有圆点的实线表示声波波形,光滑实线表示黏声介质波形。图5a和图5b分别显示对图3a的二维声波介质进行正演的单炮和对图3b的黏声介质进行正演的单炮。图6和图7分别显示吸收衰减介质常规RTM结果和Q-RTM偏移结果,在图7中,振幅异常被消除,畸变相位得到校正。
为了验证本发明对复杂模型的正演效果,图8a和图8b分别显示本发明示例性实施例中的三维速度模型和理论Q值模型;图9a和图9b分别显示本发明示例性实施例中的三维声波介质波场快照和三维黏声介质波场快照;图10a和图10b分别是图9a和图9b的内部剖视图;图11a和图11b分别是图9a和图9b的内部切片图,波场快照直观地反映了地震波在声波介质与黏性介质中的不同传播规律与特征。图12a和图12b分别显示对图9a和图9b所示波场进行正演的单炮,由于衰减因子Q的存在导致了黏性介质中正演数据的主频向低频端移动,频带变窄,其会影响偏移成像特别是中深层成像效果。
图13显示常规RTM的成像结果;图14显示根据本发明的吸收衰减介质叠前深度逆时偏移成像方法的成像结果;图15a和图15b分别显示常规RTM成像对应的频谱和根据本发明的吸收衰减介质叠前深度逆时偏移成像方法对应的频谱。
不同模型测试结果对比本发明与常规声波介质的成像结果,验证了本发明的正确性与有效性。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种吸收衰减介质叠前深度逆时偏移成像方法,其特征在于,包括:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据所述拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据所述正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
5.根据权利要求4所述的吸收衰减介质叠前深度逆时偏移成像方法,其特征在于,利用高阶有限差分算法在时间空间域分别求解所述时间空间域正传控制方程和时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子。
7.一种吸收衰减介质叠前深度逆时偏移成像系统,其特征在于,所述系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值,其中L表示拟合标准线性体的个数,l=1…L;
步骤2:在所述有效频带范围内,基于所述L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;
步骤3:基于所述积分函数值,根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;
步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;
步骤5:根据所述拟合参数数据体、时间空间域正传控制方程、时间空间域反传控制方程,获得黏声介质波场正传算子和黏声介质波场反传算子,根据所述正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;
步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910912183.7A CN110658558A (zh) | 2019-09-25 | 2019-09-25 | 吸收衰减介质叠前深度逆时偏移成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910912183.7A CN110658558A (zh) | 2019-09-25 | 2019-09-25 | 吸收衰减介质叠前深度逆时偏移成像方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110658558A true CN110658558A (zh) | 2020-01-07 |
Family
ID=69039091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910912183.7A Pending CN110658558A (zh) | 2019-09-25 | 2019-09-25 | 吸收衰减介质叠前深度逆时偏移成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110658558A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115201896A (zh) * | 2022-02-17 | 2022-10-18 | 成都理工大学 | 吸收衰减介质逆时偏移方法、装置、成像方法、介质 |
CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105652321A (zh) * | 2015-12-30 | 2016-06-08 | 中国石油大学(华东) | 一种粘声各向异性最小二乘逆时偏移成像方法 |
US20160170059A1 (en) * | 2014-12-16 | 2016-06-16 | Pgs Geophysical As | Visco-acoustic reverse-time migration using pseudo-analytical method |
US20190033479A1 (en) * | 2017-07-25 | 2019-01-31 | Advanced Geophysical Technology Inc. | Memory efficient q-rtm computer method and apparatus for imaging seismic data |
-
2019
- 2019-09-25 CN CN201910912183.7A patent/CN110658558A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160170059A1 (en) * | 2014-12-16 | 2016-06-16 | Pgs Geophysical As | Visco-acoustic reverse-time migration using pseudo-analytical method |
CN105652321A (zh) * | 2015-12-30 | 2016-06-08 | 中国石油大学(华东) | 一种粘声各向异性最小二乘逆时偏移成像方法 |
US20190033479A1 (en) * | 2017-07-25 | 2019-01-31 | Advanced Geophysical Technology Inc. | Memory efficient q-rtm computer method and apparatus for imaging seismic data |
Non-Patent Citations (7)
Title |
---|
JOAKIM O. BLANCH ET AL.: "Modeling of a constant Q: Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique", 《GEOPHYSICS》 * |
LIU HONG-WEI ET AL.: "THE ALGORITHM OF HIGH ORDER FINITE DIFFERENCE PRE-STACK REVERSE TIME MIGRATION AND GPU IMPLEMENTATION", 《CHINESE JOURNAL OF GEOPHYSICS》 * |
TONG ZHOU ET AL.: "An efficient local operator-based Q-compensated reverse time migration algorithm with multistage optimization", 《GEOPHYSICS》 * |
WENYI HU ET AL.: "An efficient Q-RTM algorithm based on local differentiation operators", 《2016 SEG INTERNATIONAL EXPOSITION AND 86TH ANNUAL MEETING》 * |
周彤等: "一种黏声波方程逆时偏移成像中的衰减补偿方法", 《地球物理学报》 * |
李振春等: "《地震叠前成像理论与方法》", 30 September 2011, 中国石油大学出版社 * |
许璐等: "粘声介质逆时偏移成像方法研究", 《SPG/SEG 南京 2020 年国际地球物理会议》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115201896A (zh) * | 2022-02-17 | 2022-10-18 | 成都理工大学 | 吸收衰减介质逆时偏移方法、装置、成像方法、介质 |
CN115793058A (zh) * | 2023-02-08 | 2023-03-14 | 成都理工大学 | 局部路径频变复旅行时的计算方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2009282330B2 (en) | Estimation of soil properties using waveforms of seismic surface waves | |
CN105467444B (zh) | 一种弹性波全波形反演方法及装置 | |
Borisov et al. | Application of 2D full-waveform inversion on exploration land data | |
CN106054244B (zh) | 截断时窗的低通滤波多尺度全波形反演方法 | |
US10345466B2 (en) | Memory efficient Q-RTM computer method and apparatus for imaging seismic data | |
CN112327358B (zh) | 一种粘滞性介质中声波地震数据正演模拟方法 | |
CN110658558A (zh) | 吸收衰减介质叠前深度逆时偏移成像方法及系统 | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
Fu et al. | Multiscale phase inversion for 3D ocean‐bottom cable data | |
Mu et al. | Modeling of pure visco-qP-wave propagation in attenuating tilted transversely isotropic media based on decoupled fractional Laplacians | |
Liu et al. | An analytic signal-based accurate time-domain viscoacoustic wave equation from the constant-Q theory | |
Zhao et al. | Viscoacoustic prestack reverse time migration based on the optimal time-space domain high-order finite-difference method | |
CN111257930B (zh) | 一种黏弹各向异性双相介质区域变网格求解算子 | |
Wang et al. | Q full-waveform inversion based on the viscoacoustic equation | |
CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
Zhang et al. | 2D and 3D Q-compensated image-domain least-squares reverse time migration through the hybrid point spread functions and the hybrid deblurring filter | |
Sun et al. | The stability problem of reverse time migration for viscoacoustic VTI media | |
CN113866823A (zh) | 一种粘声各向异性介质中的正演成像方法 | |
You et al. | Q-compensated wavefield depth extrapolation-based migration using a viscoacoustic wave equation | |
Zhou et al. | High angle prestack depth migration with absorption compensation | |
CN112558146A (zh) | 基于黏声方程的正演模拟方法及系统 | |
CN108680957B (zh) | 基于加权的局部互相关时频域相位反演方法 | |
Guo et al. | Dynamic convolution-based misfit function for time domain full waveform inversion | |
Mu et al. | A simple and high-efficiency viscoacoustic reverse time migration calculated by finite difference | |
Bai et al. | Numerical simulation of seismic wavefields in TTI media using the rotated staggered-grid compact finite-difference scheme |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200107 |
|
RJ01 | Rejection of invention patent application after publication |