叠前逆时偏移成像的方法及装置制造方法及图纸

本发明专利技术提供了一种叠前逆时偏移成像的方法及装置，其中方法包括：在孔径范围内增加随机速度边界，生成随机边界速度模型；采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并且从最大接收时间反向传播地表接收到的单炮数据，将单炮数据反向传播回地下；将所述正演波场和单炮数据进行相关，形成图像。本发明专利技术提供的方法，能够提高叠前逆时偏移成像的速度，减小对于存储空间的需求。

【技术实现步骤摘要】

本专利技术实施例涉及成像
，尤其涉及一种叠前逆时偏移成像的方法及装置。
技术介绍
随着油气田勘探程度的不断深化，油气勘探日益受到关注。采用传统的偏移方法 无法圈定盐丘边界对盐下构造进行成像，而逆时偏移(Reverse Time Migration，简称RTM) 成像方法采用双波动方程，可以精确地描述波的传播过程，利用多次波、回转波等对陡倾角 以及反转构造进行准确成像，应用有限的偏移方法进行复杂地质构造成像。现有的叠前逆时偏移成像方法包括如下的步骤读取孔径内的速度模型，生成震源子波，利用有限差分方法正演震源波场到最大 接收时间，记录不同时刻的正演波场；读取相同时刻的震源波场，并同时将单炮数据规则化，从最大接收时间反向传播 地表接收的单炮数据，反向传播回地下；应用相关的成像技术，将震源波场和单炮数据进行相关，形成图像。现有技术的这种成像方法中，正演震源波场是正向传播的，而单炮数据是反向传 输的，所以，无法直接将获取到的震源波场和单炮数据直接进行相关来形成图像，而必须将 震源波场和单炮数据中的一个先进行存储，然后再与另一个通过相关形成图像，这样，就需 要巨大的存储空间。并且，由于需要先存储震源波场或者单炮数据，然后再将震源波场和单 炮数据进行相关形成图像，所以，叠前逆时偏移成像的速度慢。
技术实现思路
本专利技术提供一种叠前逆时偏移成像的方法及装置，用以解决现有技术中需要的数 据存储空间大、成像速度慢的问题，减小逆时偏移成像所需的存储空间的需求，提高逆时偏 移成像的速度。本专利技术提供了一种叠前逆时偏移成像的方法，包括读取预设参数，并读取孔径内的速度模型；在孔径范围内增加随机速度边界，生成随机边界速度模型；采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接收时间反 向传播最后两个正演波场，并且从最大接收时间反向传播地表接收到的单炮数据，将所述 单炮数据反向传播回地下；将所述正演波场和单炮数据进行相关，形成图像。本专利技术还提供了一种叠前逆时偏移成像的装置，包括预处理模块，用于读取预设参数，并读取孔径内的速度模型；生成模块，用于在孔径范围内增加随机边界速度，生成随机边界速度模型；处理模块，与所述生成模块连接，用于采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并且从最大接收时间反向 传播地表接收到的单炮数据，反向传播回地下；相关模块，与所述处理模块连接，用于将所述正演波场和单炮数据进行相关，形成 图像。本专利技术提供的叠前逆时偏移成像的方法及装置，采用随机边界速度模型正演震源 波场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并从最大接收时间反 向传播地表接收的单炮数据，将单炮数据反向传播回地下；正演波场和从地表接收的单炮 数据均为反向传播，这样，在一个时刻可以同时获得正演波场和单炮数据的波场分布，可以 将各个时刻获得的正演波场和单炮数据相关后获得图像。而无需先存储正演波场或者单炮 数据的波场分布数据，这样就能够减小对于数据存储空间的需求，并且逆时偏移成像速度 快。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以 根据这些附图获得其他的附图。图1所示为本专利技术叠前逆时偏移成像的方法实施例一的流程图；图2所示为本专利技术叠前逆时偏移成像的方法实施例二的流程图；图3所示为本专利技术叠前逆时偏移成像的装置的结构示意图；图4所示为本专利技术实施例中涉及到的观测系统示意图。具体实施例方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例 中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员 在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。如图1所示为本专利技术叠前逆时偏移成像的方法实施例一的流程图，包括步骤101、读取预设参数，并读取孔径内的速度模型。预设参数可以包括偏移参数 和单炮数据等。步骤102、在孔径范围内增加随机速度边界，生成随机边界速度模型。步骤103、采用随机边界速度模型，正演震源波场到最大接收时间，从最大接收时 间反向传播最后两个正演波场，并且从最大接收时间反向传播地表接收的单炮数据，反向 传播回地下。现有技术中采用的速度模型具有吸收边界，采用这种速度模型正演波场时，正演 波场不能反向传播，如果反向传播会导致波场失真。在本专利技术的步骤103中，采用的是步骤102中基于孔径范围内增加的随机速度边 界生成的随机边界模型，正演波场可以反向传播，而不会导致正演波场失真。步骤104、将正演波场和单炮数据进行相关，形成图像。正演波场和单炮数据是随着时间变化的，每一个时刻都有相应的波场分布。在步骤104中，具体可以是将每一时刻的 正演波场和单炮数据进行相关。在步骤104中，由于步骤103中的正演波场和单炮数据都是反向传播的，这样可以 在每一个时刻同时得到正演波场和单炮数据，无需将正演波场或单炮数据每一个时刻的波 场分布先存储后再进行相关，从而减小了数据存储空间的需求。本专利技术实施例一提供的方法，采用随机边界速度模型正演震源波场到最大接收时 间，从最大接收时间反向传播最后两个正演波场，并从最大接收时间反向传播地表接收的 单炮数据，反向传播回地下；正演波场和从地表接收的单炮数据均为反向传播，这样，在一 个时刻可以同时获得正演波场和单炮数据的波场分布，可以将各个时刻获得的正演波场和 单炮数据相关后获得图像。而无需先存储正演波场或者单炮数据的波场分布数据，这样就 能够减小对于数据存储空间的需求，并且逆时偏移成像速度快。如图1所示的实施例中，步骤103可以由图形处理器(Graphic Processing Unit, 简称GPU)执行。GPU是从中央处理器(Central Processing Unit，简称CPU)衍生出来的一个概 念，它是显卡的心脏，和显卡板载内存成为一个子系统，共同决定计算机系统的图形处理性 能。GPU可以支持多线程处理。在本专利技术实施例一中，是通过将正演波场和单炮数据相关来形成图像的。相关条 件如公式(1)所示image(x,z)= YjS(XiZJ)RiX^J) ( J )time公式(1) φ, image (χ, ζ)是成像结果，S (x，z，t)是正演波场，R(x，z，t)是单炮数 据反向传播的波场。公式(1)所示的相关条件是一种互相关成像条件，该相关条件容易实现，便于并 行，不存在稳定性问题，并且可以处理多波至问题，不会丢失波场信息。对于一个观测系统而言，一个点可以用(X，ζ)表示，波场信息的坐标定位可以通 过网格点定义。波场传播时以及应用公式(1)所示的条件成像时，空间中的每个网格点都 是解耦的，即彼此之间相互独立，所有的网格点可以并行计算，并行粒度很小，这样每个网 格点的计算就可以通过GPU来进行。GPU与CPU相比，在处理这种小粒度并行的计算时，更 具有优势。在步骤103中，可以由GPU采用随机边界速度模本文档来自技高网...

【技术保护点】
一种叠前逆时偏移成像的方法，其特征在于，包括：　　读取预设参数，并读取孔径内的速度模型；　　在孔径范围内增加随机速度边界，生成随机边界速度模型；　　采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并且从最大接收时间反向传播地表接收到的单炮数据，将所述单炮数据反向传播回地下；　　将所述正演波场和单炮数据进行相关，形成图像。

【技术特征摘要】
一种叠前逆时偏移成像的方法，其特征在于，包括读取预设参数，并读取孔径内的速度模型；在孔径范围内增加随机速度边界，生成随机边界速度模型；采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并且从最大接收时间反向传播地表接收到的单炮数据，将所述单炮数据反向传播回地下；将所述正演波场和单炮数据进行相关，形成图像。2.根据权利要求1所述的方法，其特征在于，采用所述随机边界速度模型，正演震源波 场到最大接收时间，从最大接收时间反向传播最后两个正演波场，并且从最大接收时间反 向传播地表接收的单炮数据，反向传播回地下，包括由图形处理器采用所述随机边界速度模型，正演震源波场到最大接收时间，从最大接 收时间反向传播最后两个正演波场，并且由图形处理器从最大接收时间反向传播地表接收 的单炮数据，将所述单炮数据反向传播回地下。3.根据权利要求2所述的方法，其特征在于，由图形处理器采用所述随机边界速度模 型，正演震源波场到最大接收时间，并且由图形处理器从最大接收时间反向传播地表接收 的单炮数据，将所述单炮数据反向传播回地下，包括由图形处理器采用所述随机边界速度模型，将波场值分别映射到图形处理器的各个线 程上，正演震源波场到最大接收时间，并且由图形处理器将单炮数据分别映射到图形处理 器的各个线程上，分别在各个线程上将从最大接收时间反向传播地表接收的单炮数据，将 所述单炮数据反向传播回地下。4.根据权利要求1-3中任一权利要求所述的方法，...

【专利技术属性】
技术研发人员：李博，
申请(专利权)人：中国科学院地质与地球物理研究所，北京吉星吉达科技有限公司，
类型：发明
国别省市：11[]