[发明专利]基于地质构造散点数据自动求取含烃面积的方法

权利要求书

1.一种基于地质构造散点数据自动求取含烃面积的方法，其特征在于，包括如下步骤：

第一步：利用道格拉斯-普克抽稀算法对原始构造散点进行抽稀，去掉冗余的点；

第二步：依据第一步抽稀后的散点，求取该散点的外包矩形，作为构造等值线网格的边界；

第三步：求取所有散点间距集合，从集合中获得最小间距，取整后作为网格的纵横间距(或者说网格步长)；

第四步：基于扫描线算法原理，求抽稀后的等值线与各条网格线的交点集合；

第五步：基于所有交点集合，依次计算交点所在单元格的顶点值；先利用加权平均法计算单元格左上顶点，再依据左上顶点所在四个单元格的边，进一步利用线性关系计算各边的另一个顶点，如此递归，直到完成所有交点所在单元格的顶点值；

第六步：将第五步获得的所有单元格顶点作为样本点，利用反距离插值算法对未赋值的单元格顶点赋值，得到构造等值线网格；

第七步：根据流体界面深度，在第六步中得到的构造网格中追踪该深度的平面位置，进而得到含烃面积。

2.根据权利要求1所述的基于地质构造散点数据自动求取含烃面积的方法，其特征在于，

所述第一步是利用现有的道格拉斯-普克抽稀算法(Douglas–Peucker algorithm)对原始构造散点进行抽稀，去掉冗余的点；

所述第二步是依据抽稀后的散点，求取其外包矩形，作为构造等值线网格的边界是依次取出抽稀后的点，比较其X，Y坐标，将获得的最小X、最大X、最小Y、最大Y坐标作为外包矩形四条边的坐标值；

所述第四步的基于扫描线算法原理，求抽稀后的等值线与各条网格线的交点集合是取出每条等值线的线段，求取线段两个端点内的横、纵网格线，利用相似三角形快速计算线段与横、纵坐标的交点，即为网格的交点；

第六步：把所有获得的单元格顶点作为样本点，利用反距离插值算法对未赋值的单元格顶点赋值，得到构造等值线网格，对任意一个单元格顶点A，以三个样本点为例，分别为样本点B、样本点C和样本点D，获得顶点A赋值的反距离插值算法为：

设：B1＝1/(A到B的距离)，C1＝1/(A到C的距离)，D1＝1/(A到D的距离)，单元格A顶点的值Z＝B的Z值*B1/(B1+C1+D1)+C的Z值*C1/(B1+C1+D1)+D的Z值*D1/(B1+C1+D1)；

第七步：根据流体界面深度，在得到的构造网格中追踪该深度的平面位置，进而得到含烃面积是：在网格中提取所有深度大于流体界面深度的网格顶点及构造等值线，根据等值线两侧构造点的深度，将一侧为大于界面深度、一侧为小于界面深度的线提取出来，对其做闭合处理，得到含烃面积。