[发明专利]复杂空间的曲面重构方法

权利要求书

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

第一步、原始断层数据预处理：包括断层插值、拟合断层面和根据层位原始三维数据生成关联层位的断层多边形；其中断层插值和拟合断层面是先确定一个断层包络范围，利用克里金算法直接插值，然后根据插值数据拟合成断面；所述断层多边形是断层与层位在空间上相交而成的一个多边形结构，其边界由断层线组成，断层线分为上盘线和下盘线，上下盘线组成一个三维空间中的闭合多边形；

第二步、原始层位数据预处理：去除违法的点数据；

第三步、层位插值：

(1)初始化关联的断层多边形：获取与将要插值的层位关联的所有断层多边形，初始化层位插值的约束断层多边形边界集合；

(2)初始化网格属性：确定网格参数，将关联断层多边形边界按照网格离散为三维控制点，并将关联断层多边形投影到网格，确定各个网格的属性；

(3)初始化种子点数据：获取将要插值的层位原始数据和关联断层多边形边界的三维控制点，转换为离散的三维种子点数据结构，获取每个种子点的断层关联属性，然后将种子点按照(x，y)坐标投影到各个网格中，定义种子点结构，获得初始化种子点的属性；

(4)初始化插值点数据：在插值之前对网格的每一道进行插值点初始化，确定网格的每一道需要插值的层位点数目，以及每个插值点的断层关联属性，定义插值点结构，获得初始化插值点的属性；

(5)层位插值：根据第(4)步得到的插值点数据，通过搜索插值点周围的合法种子点来对该插值点进行插值计算；在种子点稀疏的层位中，保留无法计算Z值的插值点，在所有能依靠种子点进行插值计算的插值点处理完毕后，将插值完成的点当作新的种子点来进行层位补偿插值，直到所有插值点均处理完毕；

第四步、层位拟合成面：根据事先确定的边界约束条件将层位与断层相邻部分的边界插值点与断层多边形的上下盘线段连接起来，得到一个层面和断层完全无缝连接的面结构。

2.根据权利要求1所述的复杂空间的曲面重构方法，其特征在于：所述合法种子点的判断规则为：设待插值点为I_i，关联的断层集合uF_i和lF_i；搜索到的种子点为S_j，关联的断层集合uF_i和lF_j；定义与插值点和种子点的水平连线I_iS_j相交的断层多边形边界集合Boundary，其中uP为断层上盘边界集合，lP为断层下盘边界集合，up_k为断层F_k关联插值层位的多边形上盘边界，lp_k为断层F_k关联插值层位的多边形下盘边界，m为插值层位关联的断层多边形总数：

Boundary＝{uP，lP}

uP={upi1,upi2,L,upiM},0≤iM≤m]]>

lP={lpk1,lpk2,L,lpkN},0≤kN≤m]]>

A.若目则种子点S_j合法；

B.若存在断层F_k∈uF_i且F_k∈lF_j，或者F_k∈lF_i且F_k∈uF_j，则种子点S_j不合法；

C.若Boundary中存在垂直断层边界，则种子点S_j不合法；

D.若Boundary中同时存在某个断层的上盘边界和下盘边界，则种子点S_j不合法；

E.若uF_j中最后一个断层或者lF_j中第一个断层为正断层，且Boundary中包含该正断层的边界，则种子点S_j不合法；

F.若Boundary集合不满足D和E中的条件，同时有uF_i＝uF_j且lF_i＝lF_j，则种子点S_j合法；

G.若存在断层F_k∈uF_i且lp_k∈F_k且lp_k∈lP，或者F_k∈lF_i且up_k∈F_k且up_k∈uP，则种子点S_j不合法；

H.若存在断层F_k∈uF_j且lp_k∈F_k且lp_k∈lP，或者F_k∈lF_j且up_k∈F_k且up_k∈uP，则种子点S_j不合法；

I.其它情况，则判定种子点S_j合法。