地震波数值模拟方法、装置、介质及电子设备制造方法及图纸

本发明专利技术公开了地震波数值模拟方法、装置、介质及电子设备，该方法包括：生成全模型非结构网格；在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；删除希尔伯特空间填充曲线序列中的空序号；对计算范围内的网格单元构成的区域进行划分；基于划分后的区域进行地震波数值模拟计算。本发明专利技术利用基于空间填充曲线的区域划分，将网格划分为不同的计算区域，保持网格之间的相邻特征，保证各个子区域的网格单元数目相近，提高对非结构网格的划分效率，减少计算消耗，地震波数值模拟计算效率高。效率高。效率高。

【技术实现步骤摘要】
地震波数值模拟方法、装置、介质及电子设备


[0001]本专利技术石油勘探
，具体涉及一种基于空间填充曲线区域划分的地震波数值模拟方法、装置、介质及电子设备。

技术介绍

[0002]地震数值模拟是研究地震波传播特征，分析地质构造的可靠方法。地震数值模拟的方法主要有有限差分方法、有限元方法、伪谱法等。这些方法均依赖对模型的网格剖分。其中有限元方法具有精度较高，方法的适用性更强的特点。地震波数值模拟中的有限元网格通常是非结构性的四面体网格。
[0003]地震数值模拟所要面对的是非常巨大的数值模型。在地震勘探领域，整个模型的范围一般在几百立方千米，有限单元的网格数量在千万级到亿级的规模。因此模拟时通常使用并行计算的方法。目前最普遍的并行计算策略是将计算任务均分到不同计算单元，比如普通计算机的不同CPU核，或者超级计算的不同网络节点。每个计算单元会被分配到模型中一个固定区域的模拟任务。不同计算单元之间会在一定时间节点上进行计算结果的相互通信。上述并行任务的前提是对模型的区域划分。
[0004]对于模型的区域划分目前主要有两种实现方式本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于空间填充曲线区域划分的地震波数值模拟方法，其特征在于，包括：生成全模型非结构网格；在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；使用希尔伯特空间填充曲线填充所述计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；基于所述质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分；基于划分后的区域进行地震波数值模拟计算。2.根据权利要求1所述的基于空间填充曲线区域划分的地震波数值模拟方法，其特征在于，在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在所述全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单元。3.根据权利要求1所述的基于空间填充曲线区域划分的地震波数值模拟方法，其特征在于，所述建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系；其中，所述希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。4.根据权利要求3所述的基于空间填充曲线区域划分的地震波数值模拟方法，其特征在于，所述基于所述质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号包括：在希尔伯特空间填充曲线序列中，删除没有包含质心单元的希尔伯特单元对应的序号。5.根据权利要求1所述的基于空间填充曲线区域划分的地震...

【专利技术属性】
技术研发人员：韩德超，杨丽，佘德平，赵磊，司文朋，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：