1.一种基于最小曲率法插值的含断层的等值线图自动生成方法，其特征在于，包括：步骤1，网格化，确定研究区域，并在所述研究区域内建立网格；

步骤2，三角化，将每个所述网格以其中一条对角线划分为两个三角形，并将整个所述研究区域建立一个互不交叉的三角形列表；

步骤3，数据的预处理，建立点列表、边列表和三角形列表，同时记录边类型；

步骤4，网格点的插值，采用最小曲率法插值方法对每个所述网格点插值；

步骤5，等值线追踪，采用虚拟点技术进行所述等值线追踪；

步骤6，填充区域的搜索，建立基于边框搜索线和断层搜索线形成的环结构，按照逆时针方向遍历搜索线寻找所述填充区域；

步骤7，等值线数据的格式化，所述等值线数据生成后，将所述等值线数据和所述填充区域的数据格式化后保存；

步骤3中所述的边类型包括边框边、内部边、断层边和虚拟边，所述边框边是指边位于边框上；所述断层边是指边有断层通过；所述虚拟边是指边本身不存在，随着等值线的追踪动态创建；所述内部边是指除所述边框边、所述断层边和所述虚拟边的边；

步骤4中所述网格点插值主要包括：

(1)所述网格点的迭代方程的建立；

(2)所述网格点的初始值设置；

(3)所述网格点的迭代方程的迭代求解；

不同所述网格点的迭代方程包括：分别是远离边框的网格点迭代方程，如公式a；边框外的虚拟点的迭代方程；靠近断层的虚拟点的迭代方程；受原始数据点约束的网格点的迭代方程，如公式d；和断层端点附近的网格点的迭代方程；

所述公式a，如下：

公式(a)

其中，α为网格的横纵比 Δx为网格矩形的长，Δy为网格矩形的宽，TI为内部张力，默认为0；

所述边框外的虚拟点的迭代方程，包括：框外的第一层虚拟点的迭代方程，如公式b-1；边框外的第二层虚拟点的迭代方程，如公式b-2；边框外的位于角落的虚拟点的迭代方程，如公式b-3： 公式(b-1)

2 2

z-20＝z20+α(z11+z1-1-z-11-z-1-1)-2(1+α)(z10-z-10) 公式(b-2)z-1-1＝z1-1+z-11-z11 公式(b-3)所述靠近断层的虚拟点的迭代方程:所述虚拟点所属的每个网格点都有对应的迭代方程加权平均，得到所述虚拟点的迭代方程；

公式(d)

所述断层端点附近的网格点的迭代方程：虚拟点A的迭代方程由所述虚拟点A所属的网格点的迭代方程和与所述虚拟点A有连通关系的网格点的迭代方程加权平均得到；

其中，z01，z0-1，z10，z-10表示目标网格点的四个内环点的高程值，z11，z-11，z1-1，z-1-1表示目标网格点的四个中环点的高程值，z02，z0-2，z20，z-20表示目标网格点的四个外环点的高程值；

TB表示边框张力，取值范围是0到1；

z1、z2、z3和z4表示相对于原始数据点和目标网格点的临近网格点的高程值，z5表示原始数据点的高程值；

b1、b2、b3、b4表示临近网格点相对于目标网格点的膨胀因子，b5表示原始数据点相对于目标网格点的膨胀因子；

在建立所述断层端点附近的网格点的迭代方程时，通过周围所述网格点加权平均来设置；

建立边框附近的所述网格点的迭代方程和断层附近的所述网格点迭代方程时，当出现邻近网格点不足的情况时，附加两排虚拟网格点，通过所述虚拟网格点所属的所述网格点的迭 代加权平均来建立所述虚拟网格点的迭代方程；

所述网格点的初始值的设置方法如下：以受所述原始数据点约束的所述网格点为种子，向上、下、左、右以及左上、右上、左下、右下八个方向扩散，在扩散过程中查看对应的边是否为断层边，如果为断层边，则说明被断层阻断，朝此方向的扩散结束；循环每个种子点的扩散，每个所述网格点都会得到一个关于所述种子点和扩散次数的列表，根据所述列表，通过反距离加权法计算所述网格点的初始值；

所述等值线追踪采用基于三角形的基本原理，同时采用所述虚拟点技术，拓展断层单侧数据的连续性，在所述等值线追踪过程中，时刻检查所述等值线是否和断层相交，并选择所述等值线和断层相交的点做为停靠点。

2.根据权利要求1所述的基于最小曲率法插值的含断层的等值线图自动生成方法，其特征在于，所述研究区域为包含原始数据点和断层线的最大矩形，所述原始数据点包含的属性有横坐标x、纵坐标y、高程值z。

3.根据权利要求1所述的基于最小曲率法插值的含断层的等值线图自动生成方法，其特征在于，建立所述网格点的迭代方程时，首先建立以所述网格点为中心的中心区域，在所述中心区域内，通过将多个所述原始数据点加权平均，使得每个所述网格点的矩形区域至多有一个原始数据点与之对应。