基于FPGA多芯片板卡全波形反演特定频率波场提取的方法技术

本发明专利技术公开了基于FPGA多芯片板卡全波形反演特定频率波场提取的方法，涉及到计算机技术领域，包括：S1、通过CPU加载速度场，观测系统数据，将观测系统数据传入FPGA计算板卡进行计算；S2、FPGA计算板卡中的主芯片按照确认的波场延拓过程，进行波场延拓计算；S3、主芯片上计算的波场值与CPU进行数据交互；S4、按照预先设置的频率因子将波场值分发给计算芯片，然后在若干计算芯片上同时对不同频率成分的波场信息进行提取，并将提取的数据信息传递给主芯片；S5、在完成S1至S4的过程的同时主芯片更新迭代时间，将时间向后推演，直到整个波场延拓完成。本发明专利技术将波场数据中特定成分提取步骤放在其它计算芯片上异步完成，提高了数据复用率，减少了数据的传输。减少了数据的传输。减少了数据的传输。

Method of extracting specific frequency wave field by full waveform inversion of multi chip board based on FPGA

【技术实现步骤摘要】
基于FPGA多芯片板卡全波形反演特定频率波场提取的方法


[0001]本专利技术涉及到计算机
，尤其涉及到基于FPGA多芯片板卡全波形反演特定频率波场提取的方法。

技术介绍

[0002]深水勘探速度建模和成像精度受到崎岖海底，海水速度，深层复杂构造，复杂绕射多次波等问题影响。全波形反演作为利用波形匹配计算模型扰动的反演方法，能够获取高分辨率的模型参数，提高偏移成像精度并降低勘探风险。提高反演算法的效率和精度，更加快速准确的获取地下介质特征，是油气勘探和复杂的工程地球物理探测永恒的目标。地震勘探是利用地震观测资料定量估计出地下介质的相关属性，为找寻油气资源提供依据，反演是油气勘探的主要手段之一。传统的反演方法，通常只使用地震响应中的部分信息，这就导致地震响应中的有些有用信息被忽略。随着油气勘探复杂程度的加深，叠前全波形反演(Full Waveform Inversion,FWI)方法作为一种能够客观反应地震波传播规律，适用于任意地质模型为研究地球深部构造、中深部找矿找油、浅地表地质环境调查以及工程建设等工作提供技术支持的方法越来越受到人们的重视。
[0003]全波形反演是基于逆时偏移技术发展起来的一种新的地震成像技术，不同的是逆势偏移技术反传的是地震记录，并在模型空间得到一次成像，模型并没用得到更新；而全波形反演技术反传的是地震记录的残差，模型不断的更新直至野外观测记录和模拟记录的匹配达到技术要求，最终的反演模型就是我们所能求取的最佳模型，但需要提供满足一定精度的初始模型。拥有相似的理论框架的逆时偏移技术的发展极大地促进了全波形反演技术的理论完善和实际应用。全波形反演包含两个过程即，正演过程与反演过程。正演是反演的重要组成部分，正演方法的选择决定着反演的精度和效率。因此，能否正确的选择正演方法决定了全波形反演的可行性，尤其是3D多炮点勘探情况下。三维情况下仅仅使用时间域算法进行全波形的正反演，其反演过程需要巨大的存储空间，耗费较长的CPU时间；直接求解法的正演过程对内存存储要求极高，炮集的LU分解处理需要的内存通常达到几百GB甚至TB数量级，并且难以实现并行，耗费时间长。为避免这一问题，在处理3D情况下的全波形反演将采取时间域正演方法。在正演得到时间域地震记录后灵活运用特定频率成分提取的方式，将时间域波场转换到频率域来实现频率域的全波形反演，这样我们不仅得到了频率域的波场值，还可以在时间域对数据进行所需要的预处理，并且一次正演可以得到多个所需频率的波场值。反演也可以在时间域或频率域实现。
[0004]而全波形反演对计算资源的要求是制约其发展的瓶颈之一，计算和存储量巨大这是目前需要解决的首要问题。地震勘探多年来一直是交叉领域，涉及到数学、物理学以及地球物理学等多学科。近年来随着计算机技术的发展，计算机的性能不断提高，资料处理能力也在逐年提高。3D全波形反演实际资料的应用开始陆续出现。目前，全波形反演的困难主要来源于巨大的计算量、庞大的存储空间需求以及问题高度不适定性所带来的计算不稳定性。

技术实现思路

[0005]本专利技术的目的在于提供一种基于FPGA多芯片板卡全波形反演特定频率波场提取的方法，用于解决上述技术问题。
[0006]本专利技术采用的技术方案如下：
[0007]基于FPGA多芯片板卡全波形反演特定频率波场提取的方法，包括：
[0008]S1、通过CPU加载速度场，观测系统数据，并将所述系统数据传入FPGA计算板卡进行计算；
[0009]S2、所述FPGA计算板卡中的主芯片按照确认的波场延拓过程，进行波场延拓计算；
[0010]S3、所述主芯片上计算的波场值与所述CPU进行数据交互；
[0011]S4、按照预先设置的频率因子将所述波场值分发给若干个计算芯片，然后在若干所述计算芯片上同时对不同频率成分的波场信息进行提取，并将提取的数据信息传递给所述主芯片；
[0012]S5、在完成S1至S4的过程的同时所述主芯片更新迭代时间，将时间向后推演，直到整个波场延拓完成。
[0013]作为优选，在S1中，所述系统数据包括不同频率成分的波场因子数据。
[0014]作为优选，在S2中，采用有限差分法进行波场延拓计算。
[0015]作为优选，在S3中，所述主芯片与所述CPU之间通过PCIE接口进行数据交互。
[0016]作为优选，在S3中，当进行数据交互后，需要在CPU平台上根据观测系统位置关系记录时刻的波场快照值，并将波场快照值写入磁盘。
[0017]作为优选，在S4中，当提取的数据信息传递给所述主芯片后，所述主芯片将提取的数据信息传递至所述CPU，并在CPU平台上将所述数据信息存盘。
[0018]上述技术方案具有如下优点或有益效果：
[0019](1)本专利技术中，通过利用FPGA的并行计算，在CPU侧预留时间域数据预处理模块，利用FPGA硬件逻辑电路实现算法的优势，采用流水式数据处理方式，将波场数据中特定成分提取步骤放在其他几颗计算芯片上异步完成，大大提高了数据的复用率，同时大大减少了数据的传输，通过这种架构大大提高了整个计算系统的并发能力，提高计算系统的吞吐能力，大大缩短了算法的计算时间；
[0020](2)本专利技术中的计算架构方式，使大尺度的三维全波形反演在频率反演的方法可以利用在时间域正演计算的波场值信息，大大利用FPGA针对地震正演封装的算子库，加大FPGA开发的bit的利用率；
[0021](3)本专利技术中，使得全波形反演所需的算力可以利用FPGA异构的方式满足，从而大大减少算力不足的问题；
[0022](4)本专利技术中，利用异构架构的方式可以灵活的在CPU平台上，在时间域对数据进行所需要的预处理，并且一次正演可以得到多个所需频率成分的波场值，对后续的反演在时间域或频率域实现提供了方便；
[0023](5)本专利技术中，将频率因子通过参数形式传入FPGA中，大大增强了整个系统的灵活性和可编程性；且通过测试，其加速性能也是普通加速方案的200倍以上，整个设计方案为地震勘探中的难题解决提供必要的算力支持，也为整个地震勘探速度建模提供了可行的算力支持。
附图说明
[0024]图1是本专利技术中的主芯片与计算芯片之间的架构图；
[0025]图2是本专利技术中的全波形反演特定频率成分波场信息提取模块流程图；
[0026]图3是本专利技术中基于FPGA多芯片板卡全波形反演特定频率波场提取的方法的流程图。
具体实施方式
[0027]下面将结合附图对本专利技术的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0028]在本专利技术的描述中，需要说明的是，如出现术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等，其所指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于FPGA多芯片板卡全波形反演特定频率波场提取的方法，其特征在于，包括：S1、通过CPU加载速度场，观测系统数据，并将所述系统数据传入FPGA计算板卡进行计算；S2、所述FPGA计算板卡中的主芯片按照确认的波场延拓过程，进行波场延拓计算；S3、所述主芯片上计算的波场值与所述CPU进行数据交互；S4、按照预先设置的频率因子将所述波场值分发给若干个计算芯片，然后在若干所述计算芯片上同时对不同频率成分的波场信息进行提取，并将提取的数据信息传递给所述主芯片；S5、在完成S1至S4的过程的同时所述主芯片更新迭代时间，将时间向后推演，直到整个波场延拓完成。2.如权利要求1所述的基于FPGA多芯片板卡全波形反演特定频率波场提取的方法，其特征在于，在S1中，所述系统数据包括不同频率成分的波场因子数据。3.如权利...

【专利技术属性】
技术研发人员：龙瞻，瞿雷，许康，曹腾杰，张子义，
申请(专利权)人：上海雪湖科技有限公司，
类型：发明
国别省市：