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

权利要求书

1.一种复杂空间曲面的拟合方法，其特征在于：包括如下步骤：

(1)获取工区范围及原始数据；

(2)读取曲面数据，划分大小网格，完成曲面原始数据投影；

(3)读取断层数据，直接绘制断层线或投影断层数据，并拟合断层面；

(4)采用空间网格化插值算法对所有的待插值点进行插值：

设层位数据投影的种子点记为B_i，所有种子点放入集合U_i＝{B₁，B₂，...，B_n}，1≤i≤n，n为原始的种子点个数；待插值点记为D_j，则所有待插值点的集合V＝{D₁，D₂，...，D_m}，1≤j≤m，m为带插值点的个数；记完成插值的点为D′_i，则所有已插值点集合U_k＝{D′₁，D′₂，...，D′_m}，0≤k≤m；

插值开始前，集合U_i中存放所有的原始种子点，集合U_k为空，k≥2，集合V中存放所有待插值点；

第一步，从集合U_i中取出一个种子点B_i，递归搜索临近点，判断种子点与临近点之间是否被断层分割：如果被分割，则该临近点不是待插值点；如果不被分割，则该临近点是待插值点，需要插值；

第二步，对需要插值的点，采用基于矩形网格多重剖分的快速搜索方法搜索与该点临近且不被断层分割的种子点，对该点进行插值；若找到的种子点数大于等于两个，用反比例加权法对D_j点的高程值进行估计；若找到的种子点数少于两个，直接将种子点的值赋于D_j；

第三步，验证插值是否正确，确保断层分割区域有多个不同高程值，且每两个值都是被断层分割的；

第四步，将插值后的点D_j标记为种子点，放入集合U_k，并且从集合V中删除D_j；

第五步，判断集合U_i是否为空，若集合U_i不为空，则跳至第一步；若集合U_i为空，再判断集合V是否为空，若集合V不为空，令U_i＝U_k，依次进行第一步到第四步并将插值后的点放入集合U_k+1；若集合V为空，则插值完毕；

(5)拟合空间曲面。

2.根据权利要求1所述的复杂空间曲面的拟合方法，其特征在于：所述种子点与临近点之间是否被断层分割的判断方法如下：

设B为种子点，D为待插值点，Z₁、Z₂分别为B、D上的层位高程值，F₁、F₂分别为B、D上的断层高程值；

①当B、D上都有断层点时，如果满足条件Z₁≥F₁且Z₂≤F₂，或者满足条件Z₁≤F₁且Z₂≥F₂，则B、D点被断层分割；

②当B上有断层点，D上没有断层点时，如果满足条件Z₁≥F≥Z₂，或者满足条件Z₁≤F≤Z₂，则B、D点被断层分割；

③当B上没有断层点，D上有断层点时，如果满足条件Z₁≥F≥Z₂，或者满足条件Z₁≤F≤Z₂，则B、D点被断层分割；

④其余情况下，B、D点不被断层分割。