System.ArgumentOutOfRangeException: 索引和长度必须引用该字符串内的位置。 参数名: length 在 System.String.Substring(Int32 startIndex, Int32 length) 在 zhuanliShow.Bind() 一种黏声介质地震成像方法、装置、电子设备及介质制造方法及图纸_技高网

一种黏声介质地震成像方法、装置、电子设备及介质制造方法及图纸

技术编号:41749496 阅读:5 留言:0更新日期:2024-06-21 21:34
本申请公开了一种黏声介质地震成像方法、装置、电子设备及介质。成像方法包括:获得当前炮数据不同时刻的正传波场,所述正传波场不含拟微分项;基于所述正传波场,依据波动方程,求解当前炮数据在不同时刻的反传波场;对于同一时刻,对所述正传波场与所述反传波场进行相关成像;对于传播中的每一个时刻都进行相关成像求和,输出当前炮集的成像结果,对所有炮集循环上述步骤得到最终成像结果。

【技术实现步骤摘要】

本专利技术涉及地震波场偏移成像领域,更具体地,涉及一种黏声介质地震成像方法、装置、电子设备及介质


技术介绍

1、在黏声介质覆盖探区如海上气云或烟囱,存在由衰减引起的振幅能量衰减、频散相位畸变、走时误差以及高频成分信号的损失,严重影响着常规逆时偏移的成像精度,特别地会影响衰减区下覆地层的构造成像信息,这些地层的吸收衰减效应也会影响着后续的油藏预测与定量解释的正确性,而这种效应需要在岩石特性表征和avo分析中加以考虑。因此,对偏移成像进行黏声介质的衰减补偿具有重要的意义。反q滤波技术在实际产生中得到广泛应用,其算法本质上是对叠后数据进行空变或时变的滤波处理,由于数据的干涉作用该方法不适用于远偏移距且对于强衰减区域也不适用。相对于信号滤波方法,基于偏移成像的算法在理论上要先进很多,其中基于射线的qpsdm偏移方法通过沿射线路径实现吸收衰减的补偿,从而改善了地震资料的avo属性,然而由于基于高频近似,较难处理多波至问题,在复杂介质模型中的应用受到限制。基于波动理论的单程波偏移方法以波动方法为理论基础,可以实现波场更准确的传播,其算法是在频域中进行实现,而反映地层吸收衰减效应的q值与频率密切相关,因此,其实现较为容易高效,但受理论的限制单程波算法对于高陡构造不能准确成像。对于复杂地质构造,逆时偏移算法被认为是目前成像精度最高的成像方法,相较于克希霍夫和单程波具有显著的成像优势,现阶段是深度域成像的常规技术,但是对于基于黏声介质的逆时偏移成像技术的研究并不多。目前基于逆时偏移的q补偿算法还是基于伪谱法进行实现,其采用的是全局传播算子,其特点是求解容易,精度较高,但高维傅里叶正反变换使得该方法计算效率非常低,需要付出巨大的计算量和存储,并且不合适大规模的并行化处理,所以此类方法在实际生产应用中并没有得到广泛的应用与推广。

2、因此,有必要开发一种黏声介质地震成像方法、装置、电子设备及介质。

3、公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现思路

1、本专利技术提出了一种黏声介质地震成像方法、装置、电子设备及介质,其在时间空间域通过有限差分的方式高效地实现了黏声介质方程的数值求解。

2、第一方面,本公开实施例提供了一种黏声介质地震成像方法,包括:

3、获得当前炮数据不同时刻的正传波场,所述正传波场不含拟微分项;

4、基于所述正传波场,依据波动方程,求解当前炮数据在不同时刻的反传波场;

5、对于同一时刻,对所述正传波场与所述反传波场进行相关成像;

6、对于传播中的每一个时刻都进行相关成像求和,输出当前炮集的成像结果,对所有炮集循环上述步骤得到最终成像结果。

7、作为本公开实施例的一种具体实现方式,对当前炮数据利用下式的时空域黏声介质波动方程,求解所述正传波场,

8、

9、

10、其中σ为应力场,v为介质速度,γ为辅助变量,τσ,τξ为应变参数,k为衰减模型参数阶数,为x,z方向的偏导数。

11、作为本公开实施例的一种具体实现方式,所述时空域黏声介质波动方程通过下式的频散关系推导而出,

12、

13、其中mc,mr为复模量,w为频率,τσ,τξ为应变参数,k为衰减模型参数阶数。

14、作为本公开实施例的一种具体实现方式,所述正传波场与所述反传波场利用下式进行相关成像:

15、

16、其中sc(t,x),rc(t,x)表示t时刻空间位置x=(x,z)处考虑了地层吸收补偿的炮点正传波场与检波点反传波场,t为波场延拓最大时间。

17、第二方面,本公开实施例还提供了一种黏声介质地震成像装置,包括:正传波场模块,用于获得当前炮数据不同时刻的正传波场,所述正传波场不含拟微分项;

18、反传波场模块,用于基于所述正传波场,依据波动方程,求解当前炮数据在不同时刻的反传波场;

19、成像模块,用于对于同一时刻,对所述正传波场与所述反传波场进行相关成像;对于传播中的每一个时刻都进行相关成像求和,输出当前炮集的成像结果,进而得到所有炮集的最终成像结果。

20、作为本公开实施例的一种具体实现方式,所述正传波场模块对当前炮数据利用下式的时空域黏声介质波动方程,求解所述正传波场,

21、

22、

23、其中σ为应力场,v为介质速度,γ为辅助变量,τσ,τξ为应变参数,k为衰减模型参数阶数,为x,z方向的偏导数。

24、作为本公开实施例的一种具体实现方式,,所述时空域黏声介质波动方程通过下式的频散关系推导而出,

25、

26、其中mc,mr为复模量,w为频率,τσ,τξ为应变参数,k为衰减模型参数阶数。

27、作为本公开实施例的一种具体实现方式,所述正传波场与所述反传波场利用下式进行相关成像:

28、

29、其中sc(t,x),rc(t,x)表示t时刻空间位置x=(x,z)处考虑了地层吸收补偿的炮点正传波场与检波点反传波场,t为波场延拓最大时间。

30、第三方面,本公开实施例还提供了一种电子设备,该电子设备包括:

31、存储器,存储有可执行指令;

32、处理器,所述处理器运行所述存储器中的所述可执行指令,以实现所述的黏声介质地震成像方法。

33、第四方面,本公开实施例还提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现所述的黏声介质地震成像方法。

34、其有益效果在于:

35、本专利技术基于时空域的局部微分算法进行实现,克服了巨大的计算量和海量存储问题,解决了拟微分算子的求解难题,通过常规的高阶有限差分的数值求解算法进行实现,具有高效的计算能力,为该算法的实际工业化应用提供了技术支撑。

36、本专利技术的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本专利技术的特定原理。

本文档来自技高网...

【技术保护点】

1.一种黏声介质地震成像方法,其特征在于,包括:

2.根据权利要求1所述的黏声介质地震成像方法,其中,对当前炮数据利用下式的时空域黏声介质波动方程,求解所述正传波场,

3.根据权利要求2所述的黏声介质地震成像方法,其中,所述时空域黏声介质波动方程通过下式的频散关系推导而出,

4.根据权利要求1所述的黏声介质地震成像方法,其中,所述正传波场与所述反传波场利用下式进行相关成像:

5.一种黏声介质地震成像装置,其特征在于,包括:

6.根据权利要求5所述的黏声介质地震成像装置,其中,所述正传波场模块对当前炮数据利用下式的时空域黏声介质波动方程,求解所述正传波场,

7.根据权利要求6所述的黏声介质地震成像装置,其中,所述时空域黏声介质波动方程通过下式的频散关系推导而出,

8.根据权利要求5所述的黏声介质地震成像装置,其中,所述正传波场与所述反传波场利用下式进行相关成像:

9.一种电子设备,其特征在于,所述电子设备包括:

10.一种计算机可读存储介质,其特征在于,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现权利要求1-4中任一项所述的黏声介质地震成像方法。

...

【技术特征摘要】

1.一种黏声介质地震成像方法,其特征在于,包括:

2.根据权利要求1所述的黏声介质地震成像方法,其中,对当前炮数据利用下式的时空域黏声介质波动方程,求解所述正传波场,

3.根据权利要求2所述的黏声介质地震成像方法,其中,所述时空域黏声介质波动方程通过下式的频散关系推导而出,

4.根据权利要求1所述的黏声介质地震成像方法,其中,所述正传波场与所述反传波场利用下式进行相关成像:

5.一种黏声介质地震成像装置,其特征在于,包括:

6.根据权利要求5所述的黏声介质地震成像装置,其中,所述正...

【专利技术属性】
技术研发人员:许璐李博肖建恩
申请(专利权)人:中国石油化工股份有限公司
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1