[发明专利]一种分层级虚网格获取方法、装置及系统

说明书

本申请实施方式公开了一种分层级虚网格获取方法、装置及系统，分层级虚网格获取方法包括：创建基层矩形网格；根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。

技术领域

本申请涉及虚网格技术领域，特别涉及一种分层级虚网格获取方法、装置及系统。

背景技术

等值线图是地质研究工作中的一种重要的表示方法，它通过插值一些离散的点来生成平滑的地质结构图。等值线图的插值方法有很多，比较经典的方法有反距离加权、样条插值法、克里金方法、径向基函数法、最小曲率法等等。当插值约束只有离散点的情况下，这些经典方法比较简单实用，但是当需要添加如断层、断距、相控、趋势面等复杂约束的时候，这些方法实施起来就会比较困难，而由J.L.Mallet教授提出的离散光滑插值(DSI)方法却能够克服这一困难，可以非常方便的在插值过程中添加各种复杂约束。

DSI插值方法的数学模型是一个最优化模型，一般采用迭代方法对它进行求解。当采用迭代方法来求解的时候，如果其所依赖的网格模型和约束的数量级较大，那么其运算速度会比较慢，而且插值效果也会较差。

DSI插值方法所基于的网格主要包括两种：矩形网格和三角网格。对于矩形网格，因为其规则性，很容易实施分层级的划分，然而对复杂约束，例如断层线的符合程度却比较差。对于三角网格，虽然对复杂约束的符合程度非常好，但是因为它的任意性和不规则性，对于分层级的划分却非常困难。

发明内容

本申请实施方式的目的是提供一种分层级虚网格获取方法、装置及系统，针对网格数量级较大且具有复杂几何结构的模型，建立一种分层级虚网格获取方案，该方案获得的分层级虚网格能够加速DSI方法迭代求解的速度，改善插值效果。

为实现上述目的，本申请实施方式提供一种分层级虚网格获取方法，包括：

创建基层矩形网格；

根据分层级次数在所述基层矩形网格上创建每一层级的矩形网格；其中，每一层级的矩形网格是上一层级每个矩形网格十字劈分后的结果；

在第n层级的矩形网格上创建混合网格；其中，n为分层级次数的最大值；

遍历每一层级的每个矩形网格，获得当前层级的矩形网格包含的混合网格，对所述混合网格进行分类，根据分类结果创建每一层级的虚网格。

优选地，还包括：

遍历每一层级的虚网格，对每一层的虚网格进行迭代求解。

优选地，在第n层级的矩形网格上创建混合网格的步骤包括：

利用约束线对所述第n层级的矩形网格进行劈分，使得在无所述约束线的地方矩形网格保持不变；在有所述约束线的地方，每个矩形网格被劈分成三角网格，并分裂所述约束线两侧的三角网格，获得混合网格。

优选地，根据分类结果创建每一层级的虚网格的步骤包括：

遍历每一层级的矩形网格，确定当前层级的矩形网格中的混合网格内是否包含三角网格；

如果当前层级的矩形网格中的混合网格内均为矩形网格，则创建一个虚网格，所述虚网格的四个顶点是当前层级的矩形网格的四个顶点；

如果当前层级的矩形网格中的混合网格内包含三角网格，根据混合网格内三角网格和矩形网格的连通性来创建虚网格。

优选地，根据混合网格内三角网格和矩形网格之间的连通性创建虚网格的步骤包括：