[发明专利]复杂空间曲面的拟合方法

摘要

本发明公开了一种复杂空间曲面的拟合方法，通过将平面的断层投影拟合成空间中的断层面，利用断层信息作为约束条件对网格点进行插值，以消除断层对网格点插值的影响，解决了含有任意复杂断层(包括正断层和逆断层)的非连续曲面网格化问题，并通过矩形网格多重剖分的快速搜索算法实现快速网格化，可实现含有任意断层的复杂空间曲面拟合，该方法具有高效性。

主权项

一种复杂空间曲面的拟合方法，其特征在于：包括如下步骤：(1)获取工区范围及原始数据；(2)读取曲面数据，划分大小网格，完成曲面原始数据投影；(3)读取断层数据，直接绘制断层线或投影断层数据，并拟合断层面；(4)采用空间网格化插值算法对所有的待插值点进行插值：设层位数据投影的种子点记为Bi，所有种子点放入集合Ui＝{B1，B2，...，Bn}，1≤i≤n，n为原始的种子点个数；待插值点记为Dj，则所有待插值点的集合V＝{D1，D2，...，Dm}，1≤j≤m，m为带插值点的个数；记完成插值的点为D′i，则所有已插值点集合Uk＝{D′1，D′2，...，D′m}，0≤k≤m；插值开始前，集合Ui中存放所有的原始种子点，集合Uk为空，k≥2，集合V中存放所有待插值点；第一步，从集合Ui中取出一个种子点Bi，递归搜索临近点，判断种子点与临近点之间是否被断层分割：如果被分割，则该临近点不是待插值点；如果不被分割，则该临近点是待插值点，需要插值；第二步，对需要插值的点，采用基于矩形网格多重剖分的快速搜索方法搜索与该点临近且不被断层分割的种子点，对该点进行插值；若找到的种子点数大于等于两个，用反比例加权法对Dj点的高程值进行估计；若找到的种子点数少于两个，直接将种子点的值赋于Dj；第三步，验证插值是否正确，确保断层分割区域有多个不同高程值，且每两个值都是被断层分割的；第四步，将插值后的点Dj标记为种子点，放入集合Uk，并且从集合V中删除Dj；第五步，判断集合Ui是否为空，若集合Ui不为空，则跳至第一步；若集合Ui为空，再判断集合V是否为空，若集合V不为空，令Ui＝Uk，依次进行第一步到第四步并将插值后的点放入集合Uk+1；若集合V为空，则插值完毕；(5)拟合空间曲面。