本申请提供了一种全波形反演方法、装置、电子设备及存储介质；所述方法包括：设定震源，并构建所述震源对应的目标函数；基于所述震源对应的模拟数据，确定所述震源对应的震源波场；基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场；基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度；根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值。如此，能够智能地进行全波形反演，提高了全波形反演的精度和效率。提高了全波形反演的精度和效率。提高了全波形反演的精度和效率。
【技术实现步骤摘要】
一种全波形反演方法、装置、电子设备及存储介质
[0001]本申请涉及全波形反演技术，尤其涉及一种全波形反演方法、装置、电子设备及存储介质。
技术介绍
[0002]现有的大多数全波形反演是使用流变模型作为控制方程，由于流变模型表征的模型参数相对较多，且振幅衰减和相位频散作用相互耦合，现有的全波形反演的精度较低，且全波形反演的计算过程的十分耗时。人们更希望减少全波形反演的计算过程的时间并提高全波形反演的精度。
[0003]因此，如何智能地进行全波形反演，以提高全波形反演的精度和效率是一直追求的目标。
技术实现思路
[0004]本申请实施例提供了一种全波形反演方法、装置、电子设备及存储介质。
[0005]根据本申请的第一方面，提供了一种全波形反演方法，该方法包括：设定震源，并构建所述震源对应的目标函数；基于所述震源对应的模拟数据，确定所述震源对应的震源波场；基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场；基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度；根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值。
[0006]根据本申请一实施方式，所述基于所述震源对应的模拟数据，确定所述震源对应的震源波场，包括：基于所述模拟数据，确定对应的第一稳定补偿算子；所述第一稳定补偿算子包括振幅补偿算子、高频稳定算子、相位频散算子、相位校正算子中的至少之一；基于所述第一稳定补偿算子，确定所述模拟数据对应的波场延拓算子；基于所述波场延拓算子和所述震源，以时间正向的方式求解补偿方程，得的所述震源对应的正传波场；对所述正传波场以时间反向的方式进行重构，得到震源波场。
[0007]根据本申请一实施方式，所述基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场，包括：基于所述模拟数据，确定对应的第二稳定补偿算子；所述第二稳定补偿算子包括振幅补偿算子、相位频散算子中的至少之一；基于所述第二稳定补偿算子，求解衰减方程，得到所述震源对应的粘声波场；基于所述粘声波场和所述观测数据，确定伴随震源；基于所述伴随震源，确定所述逆时伴随波场。
[0008]根据本申请一实施方式，所述基于所述伴随震源，确定所述逆时伴随波场，包括：基于所述伴随震源，确定对应的第三稳定补偿算子；所述第三稳定补偿算子包括振幅补偿算子、高频稳定算子、相位频散算子、相位校正算子中的至少之一；基于所述第三稳定补偿算子，确定所述伴随震源对应的伴随算子；基于所述伴随算子和所述伴随震源，以时间反向的方式求解伴随方程，得到所述震源对应的逆时伴随波场。
[0009]根据本申请一实施方式，所述基于所述震源波场和所述逆时伴随波场，确定所述
目标函数的补偿梯度，包括：将所述震源波场和所述逆时伴随波场互相关，得到所述目标函数的补偿梯度。
[0010]根据本申请一实施方式，所述根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值之后，所述方法还包括：将最后一次迭代更新后的参数作为所述目标函数对应的全波形反演结果。
[0011]根据本申请一实施方式，所述振幅补偿算子、高频稳定算子、相位频散算子、相位校正算子均包括固定阶数的拉普拉斯算子。
[0012]根据本申请的第二方面，提供了一种全波形反演装置，所述全波形反演装置包括：设定模块，用于设定震源，并构建所述震源对应的目标函数；第一确定模块，用于基于所述震源对应的模拟数据，确定所述震源对应的震源波场；第二确定模块，用于基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场；梯度确定模块，用于基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度；迭代模块，用于根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值。
[0013]根据本申请的第三方面，提供了一种电子设备，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行本申请所述的方法。
[0014]根据本申请的第四方面，提供了一种存储有计算机指令的非瞬时计算机可读存储介质，所述计算机指令用于使所述计算机执行本申请所述的方法。
[0015]本申请实施例的方法，设定震源，并构建所述震源对应的目标函数；基于所述震源对应的模拟数据，确定所述震源对应的震源波场；基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场；基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度；根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值。如此，能够智能地进行全波形反演，提高了全波形反演的精度和效率。
[0016]需要理解的是，本申请的教导并不需要实现上面所述的全部有益效果，而是特定的技术方案可以实现特定的技术效果，并且本申请的其他实施方式还能够实现上面未提到的有益效果。
附图说明
[0017]通过参考附图阅读下文的详细描述，本申请示例性实施方式的上述以及其他目的、特征和优点将变得易于理解。在附图中，以示例性而非限制性的方式示出了本申请的若干实施方式，其中：在附图中，相同或对应的标号表示相同或对应的部分。
[0018]图1示出了本申请实施例提供的全波形反演方法的处理流程示意图一；图2示出了本申请实施例提供的全波形反演方法的处理流程示意图二；图3示出了本申请实施例提供的全波形反演方法的处理流程示意图三；
图4示出了本申请实施例提供的全波形反演方法的处理流程示意图四；图5示出了本申请实施例提供的全波形反演方法的处理流程示意图五；图6示出了本申请实施例提供的全波形反演方法的一种应用场景图；图7示出了本申请实施例提供的全波形反演方法的另一种应用场景图；图8示出了本申请实施例提供的全波形反演装置的一种可选示意图；图9示出了本申请实施例提供的电子设备的组成结构示意图。
具体实施方式
[0019]为使本申请的目的、特征、优点能够更加的明显和易懂，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而非全部实施例。基于本申请中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
[0020]在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解，“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。
[0021]对本申请实施例提供的全波形反演方法中的处理流程进行说明。参见图1，图1是本申请实施例提供的全波形反演方法的处理流程示意图一，将结合图1示出的步骤S101
‑
S105进行说明。
[0022]步骤S101，设定震源，并构建震源对应的目标函数本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种全波形反演方法，其特征在于，所述方法包括：设定震源，并构建所述震源对应的目标函数；基于所述震源对应的模拟数据，确定所述模拟数据对应的第一稳定补偿算子；所述第一稳定补偿算子包括振幅补偿算子、高频稳定算子、相位频散算子、相位校正算子中的至少之一；基于所述第一稳定补偿算子，确定所述模拟数据对应的波场延拓算子；基于所述波场延拓算子和所述震源，以时间正向的方式求解补偿方程，得到所述震源对应的正传波场；对所述正传波场以时间反向的方式进行重构，得到所述震源对应的震源波场；基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场；基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度；根据所述补偿梯度对所述目标函数的初始参数进行迭代更新，直至所述迭代更新的次数达到预设值。2.根据权利要求1所述的方法，其特征在于，所述基于所述震源对应的模拟数据和观测数据，确定所述震源对应的逆时伴随波场，包括：基于所述模拟数据，确定对应的第二稳定补偿算子；所述第二稳定补偿算子包括振幅补偿算子、相位频散算子中的至少之一；基于所述第二稳定补偿算子，求解衰减方程，得到所述震源对应的粘声波场；基于所述粘声波场和所述观测数据，确定伴随震源；基于所述伴随震源，确定所述逆时伴随波场。3.根据权利要求2所述的方法，其特征在于，所述基于所述伴随震源，确定所述逆时伴随波场，包括：基于所述伴随震源，确定对应的第三稳定补偿算子；所述第三稳定补偿算子包括振幅补偿算子、高频稳定算子、相位频散算子、相位校正算子中的至少之一；基于所述第三稳定补偿算子，确定所述伴随震源对应的伴随算子；基于所述伴随算子和所述伴随震源，以时间反向的方式求解伴随方程，得到所述震源对应的逆时伴随波场。4.根据权利要求1所述的方法，其特征在于，所述基于所述震源波场和所述逆时伴随波场，确定所述目标函数的补偿梯度，包括：将所述震源波场和所述逆时伴随波场互相关，得到所...
【专利技术属性】
技术研发人员：王宁，石颖，王维红，李松龄，关珊，袁红旗，
申请(专利权)人：东北石油大学三亚海洋油气研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。