地震波数值模拟方法、装置、介质及电子设备与流程

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


背景技术：

2.地震数值模拟是研究地震波传播特征，分析地质构造的可靠方法。地震数值模拟的方法主要有有限差分方法、有限元方法、伪谱法等。这些方法均依赖对模型的网格剖分。其中有限元方法具有精度较高，方法的适用性更强的特点。地震波数值模拟中的有限元网格通常是非结构性的四面体网格。
3.地震数值模拟所要面对的是非常巨大的数值模型。在地震勘探领域，整个模型的范围一般在几百立方千米，有限单元的网格数量在千万级到亿级的规模。因此模拟时通常使用并行计算的方法。目前最普遍的并行计算策略是将计算任务均分到不同计算单元，比如普通计算机的不同cpu核，或者超级计算的不同网络节点。每个计算单元会被分配到模型中一个固定区域的模拟任务。不同计算单元之间会在一定时间节点上进行计算结果的相互通信。上述并行任务的前提是对模型的区域划分。
4.对于模型的区域划分目前主要有两种实现方式：一是现将模型进行人为划分区域，之后对各个子区域生成网格；二是先将全模型进行网格剖分，之后对网格进行区域划分。其中方法二通过算法控制区域划分的子区域大小，可以使每个子区域的网格数量相当，有利于并行计算。
5.对于有限元方法，因为它采用非结构网格，网格的排序一般也没有规律性，模型内部网格单元数目的分布与模型的速度大小分布及几何形态直接相关。目前已有一些开源程序包可以实现非结构网格的区域划分功能。如metis，使用的是多层次递归二分切分法、多层次k路切分法以及多约束划分机制。但在面对地震模型的模拟任务时，该软件需要消耗巨大的计算资源，区域划分的效率较低。
6.在地震勘探领域，当地面观测系统是随炮滚动的方式时，每次模拟需要从模型中切出一部分进行。因此，每次模拟算法运行前都需要对当前切出的网格进行区域划分，在三维多炮计算时，因此，这个过程将更加耗时和耗费资源，计算效率较低。
7.因此需要研发一种节省计算消耗，计算效率高的地震波数值模拟方法。


技术实现要素：

8.本发明的目的是提出一种节省计算消耗，计算效率高的地震波数值模拟方法。
9.第一方面，本发明提供一种基于空间填充曲线区域划分的地震波数值模拟方法，包括：生成全模型非结构网格；在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；使用希尔伯特空间填充曲线填充所述计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；基于所述质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的
空序号；基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分；基于划分后的区域进行地震波数值模拟计算。
10.可选的，在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在所述全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单元。
11.可选的，所述建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系；其中，所述希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
12.可选的，所述基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
13.可选的，所述全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
14.第二方面，本发明还提供一种电子设备，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现上述基于空间填充曲线区域划分的地震波数值模拟方法。
15.第三方面，本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述基于空间填充曲线区域划分的地震波数值模拟方法。
16.第四方面，本发明还提供一种基于空间填充曲线区域划分的地震波数值模拟装置，包括：全模型非结构网格生成模块，生成全模型非结构网格；质心单元建立模块，在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；对应关系建立模块，使用希尔伯特空间填充曲线填充所述计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；删除空序号模块，基于所述质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；划分模块，基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分；计算模块，基于划分后的区域进行地震波数值模拟计算。
17.可选的，在所述全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在所述全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的并进行编号，建立质心单元。
18.可选的，所述建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系；其中，所述希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
19.可选的，所述基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希
尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
20.可选的，所述全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
21.本发明的有益效果在于：本发明的基于空间填充曲线区域划分的地震波数值模拟方法利用基于空间填充曲线的区域划分，将网格划分为不同的计算区域，保持网格之间的相邻特征，保证各个子区域的网格单元数目相近，提高对非结构网格的划分效率，减少计算消耗，基于划分后的区域进行数值模拟计算，地震波数值模拟计算效率高。
22.本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
23.通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
24.图1示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的流程图。
25.图2示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的全模型非结构网格示意图。
26.图3示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的当前炮在全模型非结构网格中获取的对应计算范围示意图。
27.图4示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的不同阶的hilbert空间填充曲线的生成和曲线的结果示意图。
28.图5示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的质心单元与希尔伯特空间填充曲线序列的对应关系图。
29.图6示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的区域划分后的二维和三维示意图。
30.图7示出了根据本发明的一个实施例的基于空间填充曲线区域划分的地震波数值模拟装置的框图。
31.附图标记说明
32.102、全模型非结构网格生成模块；104、质心单元建立模块；106、对应关系建立模块；108、删除空序号模块；110、划分模块；112、计算模块。
具体实施方式
33.下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
34.本发明提供一种基于空间填充曲线区域划分的地震波数值模拟方法，包括：生成全模型非结构网格；在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的
质心，建立质心单元；使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；基于质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分；基于划分后的区域进行地震波数值模拟计算。
35.具体的，生成全模型非结构网格，该网格在二维是三角形单元，三维是四面体单元，这一步不特定网格生成方法，只需生成二维三角形或三维四面体网格即可。在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元和质心坐标，获得质心单元，使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，获取质心单元与希尔伯特空间填充曲线序列的对应关系；因为非结构性单元在空间上没有明显的规律性，从而会造成在某些hilbert空间填曲线的单元中不包含质心坐标，因此在归拢相邻质心时，将这些空单元剔除。对删除空序号后的希尔伯特空间填充曲线序列进行切割，从而对与其对应的计算范围内的网格单元构成的区域进行划分，基于划分后的区域进行地震波数值模拟计算。
36.在地震勘探的随炮滚动观测方式中，每一炮的计算范围不同，需要重新按上述方法做区域划分，基于划分后的区域进行地震波数值模拟计算。
37.根据示例性的实施方式，基于空间填充曲线区域划分的地震波数值模拟方法利用基于空间填充曲线的区域划分，将网格划分为不同的计算区域，保持网格之间的相邻特征，保证各个子区域的网格单元数目相近，提高对非结构网格的划分效率，减少计算消耗，基于划分后的区域进行数值模拟计算，地震波数值模拟计算效率高。
38.作为可选方案，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单元。
39.具体的，根据观测系统模拟需求，检索出当前模拟炮的计算范围，以及该计算范围在全模型非结构网格中内的网格单元，对当前计算范围内的网格单元，用网格单元的质心坐标代替网格单元，获得质心单元。这一步的目的是减少地震模拟的计算消耗，因为将全模型的网格进行模拟会有很大一部分的区域并不需要模拟，造成资源的浪费。基于观测系统将主要区域检索出来，减少计算量，同时也不会影响模拟结果。
40.作为可选方案，建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系。
41.具体的，获取希尔伯特空间填充曲线中每个希尔伯特单元包含的质心单元的编号，将希尔伯特单元与包含的质心单元编号对应，获得质心单元与希尔伯特空间填充曲线序列的对应关系。
42.作为可选方案，希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
43.具体的，确定hilbert(希尔伯特)空间填充曲线的最小阶数，使hilbert空间填充曲线的每个最小的单元最多只有一个质心坐标。
44.作为可选方案，基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内
的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
45.具体的，按照设定的区域划分数目，切割一维hilbert空间填充曲线，即切割删除空序号后的希尔伯特空间填充曲线序列，同时实现对与希尔伯特空间填充曲线序列对应的质心单元代表的网格单元构成的区域进行划分。
46.作为可选方案，全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
47.第二方面，本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述基于空间填充曲线区域划分的地震波数值模拟方法。
48.第三方面，本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述基于空间填充曲线区域划分的地震波数值模拟方法。
49.第四方面，本发明还提供一种基于空间填充曲线区域划分的地震波数值模拟装置，包括：全模型非结构网格生成模块，生成全模型非结构网格；质心单元建立模块，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；对应关系建立模块，使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；删除空序号模块，基于质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；划分模块，基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分；计算模块，基于划分后的区域进行地震波数值模拟计算。
50.具体的，生成全模型非结构网格，该网格在二维是三角形单元，三维是四面体单元，这一步不特定网格生成方法，只需生成二维三角形或三维四面体网格即可。在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元和质心坐标，获得质心单元，使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，获取质心单元与希尔伯特空间填充曲线序列的对应关系；因为非结构性单元在空间上没有明显的规律性，从而会造成在某些hilbert空间填曲线的单元中不包含质心坐标，因此在归拢相邻质心时，将这些空单元剔除。对删除空序号后的希尔伯特空间填充曲线序列进行切割，从而对与其对应的计算范围内的网格单元构成的区域进行划分，基于划分后的区域进行地震波数值模拟计算。
51.在地震勘探的随炮滚动观测方式中，每一炮的计算范围不同，需要重新按上述方法做区域划分，基于划分后的区域进行地震波数值模拟计算。
52.根据示例性的实施方式，基于空间填充曲线区域划分的地震波数值模拟方法利用基于空间填充曲线的区域划分，将网格划分为不同的计算区域，保持网格之间的相邻特征，保证各个子区域的网格单元数目相近，提高对非结构网格的划分效率，减少计算消耗，基于划分后的区域进行数值模拟计算，地震波数值模拟计算效率高。
53.作为可选方案，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单
元。
54.具体的，根据观测系统模拟需求，检索出当前模拟炮的计算范围，以及该计算范围在全模型非结构网格中内的网格单元，对当前计算范围内的网格单元，用网格单元的质心坐标代替网格单元，获得质心单元。这一步的目的是减少地震模拟的计算消耗，因为将全模型的网格进行模拟会有很大一部分的区域并不需要模拟，造成资源的浪费。基于观测系统将主要区域检索出来，减少计算量，同时也不会影响模拟结果。
55.作为可选方案，建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系。
56.具体的，获取希尔伯特空间填充曲线中每个希尔伯特单元包含的质心单元的编号，将希尔伯特单元与包含的质心单元编号对应，获得质心单元与希尔伯特空间填充曲线序列的对应关系。
57.作为可选方案，希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
58.具体的，确定hilbert空间填充曲线的最小阶数，使hilbert曲线的每个最小的单元最多只有一个质心坐标。
59.作为可选方案，基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
60.具体的，按照设定的区域划分数目，切割一维hilbert空间填充曲线，即切割删除空序号后的希尔伯特空间填充曲线序列，同时实现对与希尔伯特空间填充曲线序列对应的质心单元代表的网格单元构成的区域进行划分。
61.作为可选方案，全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
62.实施例一
63.图1示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的流程图。图2示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的全模型非结构网格示意图。图3示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的当前炮在全模型非结构网格中获取的对应计算范围示意图。图4示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的不同阶的hilbert空间填充曲线的生成和曲线的结果示意图。图5示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的质心单元与希尔伯特空间填充曲线序列的对应关系图。图6示出了根据本发明的一个实施例的一种基于空间填充曲线区域划分的地震波数值模拟方法的区域划分后的二维和三维示意图。
64.结合图1、图2、图3、图4、图5和图6所示，该基于空间填充曲线区域划分的地震波数值模拟方法，包括：
65.步骤1：生成全模型非结构网格；
66.其中，全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
67.步骤2：在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；
68.其中，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单元。
69.步骤3：使用希尔伯特空间填充曲线填充计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；
70.其中，建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系。
71.其中，希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
72.步骤4：基于质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；
73.步骤5：基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分；
74.其中，基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
75.步骤6：基于划分后的区域进行地震波数值模拟计算。
76.以一个二维复杂地质模型的全模型非结构网格为例，总网格数为871921个单元。模型大小为x＝42.62km，z＝10km。所有网格单元均为三角形单元，单元排列没有特定方向，且单元大小也没有一定规则。根据观测系统模拟需求，检索出当前模拟炮的计算范围，以及该计算范围内对应的全模型非结构中的网格单元。如图3展示了观测范围为17km的观测系统对应的第m炮和第k炮的模型切取范围，这是一个随炮滚动的过程，每一炮的模拟都需要重新进行，在图2和图3中横坐标position表示点，depth表示位置。其中k炮的模拟范围切出的单元数为354924，仅为全模型的40％左右。在计算范围内，用单元的质心坐标代替网格单元。图4展示了不同阶数的hilbert曲线，可以看到每个阶数的曲线其实都对应了一定排列的方形小单元，这里称为hilbert曲线的单元。确定hilbert空间填充曲线的最小阶数，使hilbert曲线的每个最小的单元最多只有一个质心坐标。这一步中以三角形单元的质心为切分标准，不断提高hilbert曲线的阶数，使得每个hilbert曲线的单元中最多包含一个质心。使用hilbert空间填充曲线填充待划分区域，建立质心单元与hilbert曲线序列的对应关系。如图5展示了用hilbert曲线填充待划分区域后，每个hilbert单元内包含的三角形质心的单元编号，由此便建立了质心单元与希尔伯特空间填充曲线序列对应关系。剔除hilbert空间填充曲线上的空序号，剔除掉图5中不包含质心单元的hilbert单元后，就将二维的单元转换成了具有空间相邻特征的一维序列。按照设定的区域划分数目，切割一维hilbert空间填充曲线，从而对与其对应的计算范围内的质心单元构成的区域进行划分。图6展示了区域划分的二维和三维示意图。按照划分后的区域进行地震波数值模拟计算，在图
6中横坐标position表示点，depth表示位置，x-distance表示x方向的距离，y-distance表示y方向的距离。
77.实施例二
78.图7示出了根据本发明的一个实施例的基于空间填充曲线区域划分的地震波数值模拟装置的框图。
79.如图7所示，该基于空间填充曲线区域划分的地震波数值模拟装置，包括：
80.全模型非结构网格生成模块102，生成全模型非结构网格；
81.质心单元建立模块104，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元；
82.对应关系建立模块106，使用希尔伯特空间填充曲线填充所述计算范围内的网格单元构成的空间，建立质心单元与希尔伯特空间填充曲线序列的对应关系；
83.删除空序号模块108，基于所述质心单元与希尔伯特空间填充曲线序列的对应关系，删除希尔伯特空间填充曲线序列中的空序号；
84.划分模块110，基于删除空序号后的希尔伯特空间填充曲线序列，对所述计算范围内的网格单元构成的区域进行划分；
85.计算模块112，基于划分后的区域进行地震波数值模拟计算。
86.其中，在全模型非结构网格中获取当前模拟炮对应计算范围内的网格单元的质心，建立质心单元包括：根据观测系统的模拟需求，在全模型非结构网格中获取当前模拟炮对应的计算范围内的网格单元；确定每个网格单元的质心并进行编号，建立质心单元。
87.其中，建立质心单元与希尔伯特空间填充曲线序列的对应关系包括：基于每个希尔伯特单元包含的质心单元的编号，建立质心单元与希尔伯特空间填充曲线序列的对应关系。
88.其中，希尔伯特空间填充曲线的每个希尔伯特单元中至多包含一个质心单元。
89.其中，基于删除空序号后的希尔伯特空间填充曲线序列，对计算范围内的网格单元构成的区域进行划分包括：根据预设的区域划分数量，对删除空序号后的希尔伯特空间填充曲线序列进行切分，从而对与其对应的计算范围内的网格单元构成的区域进行划分。
90.其中，全模型非结构网格在二维空间中是三角形单元，在三维空间中是四面体单元。
91.实施例三
92.本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述基于空间填充曲线区域划分的地震波数值模拟。
93.根据本公开实施例的电子设备包括存储器和处理器。
94.该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
95.该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力
的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
96.本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
97.有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
98.实施例四
99.本公开提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述基于空间填充曲线区域划分的地震波数值模拟。
100.根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
101.上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
102.以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。