[发明专利]一种低渗气藏产能计算方法

说明书

技术领域

本发明涉及低渗气井的压裂测试技术，尤其涉及一种低渗气藏产能计算方法。

背景技术

目前的业务发展需要油藏能够准确的预测和评价压裂产能，正确的指导压裂工艺，减小无效压裂的次数，提高压裂服务的水平。

但目前还没有一种有效的方法和软件能够准确的描述裂缝形态，同时精确地预测压裂产能。主要的研究都是在常规的裂缝形态描述方法基础上进行改进和完善。近期大量的实践表明，这些现有技术存在网格极差变化大，网格模型在进行数值模拟计算时容易产生不收敛的问题。

发明内容

为了解决上述问题，本发明提出了一种低渗气藏产能计算方法，能够精细模拟压裂井的产能变化，计算出地层达到真实产能的时间和大小，节省了测试时间和费用。

为了达到上述目的，本发明提出了一种低渗气藏产能计算方法，该方法包括：

将油藏对象投影到二维平面上，获得油藏对象的二维油藏模型。

将二维油藏模型中的所有油藏对象表示成有序点集。

对有序点集进行德劳内Delaunay三角剖分，获得Delaunay三角形网。

基于Delaunay三角形网，生成Voronoi网格图。

根据Voronoi网格图计算低渗气藏产能。

优选地，油藏对象包括：油藏边界、层面、断层、井和水力压裂裂缝。

优选地，将油藏对象投影到二维平面上，获得油藏对象的二维油藏模型是指：

将油藏对象投影到二维平面上，边界转化为包络线，断层转化为多边形，井转化为点，斜度大于预定斜度阈值的井或者水平井转化为多边形，水力压裂裂缝转化为线段。

优选地，将二维油藏模型中的所有油藏对象表示成有序点集是指：将二维油藏模型中的包络线、多边形、线段以及点转化为二维平面上不同的有序点集。

优选地，该方法还包括：当油藏对象非均匀分布和/或油藏对象的有序点集中的散点个数超过预定的散点阈值时，增加二维油藏模型中的油藏对象的限定条件，将油藏对象的有序点集进行细分。

其中，增加二维油藏模型中的油藏对象的限定条件，将油藏对象的有序点集进行细分包括：增加每个和/或个别油藏对象的局部边界区域，在局部边界区域内部对该油藏对象的有序点集进行细分。

优选地，对有序点集进行德劳内Delaunay三角剖分，获得Delaunay三角形网的方法包括以下步骤：

61、遍历所述有序点集的所有散点，得出有序点集的包容盒，获得作为有序点集的凸壳的初始三角形并放入Delaunay三角形链表。

62、将有序点集中的一个散点插入Delaunay三角形链表，在Delaunay三角形链表中找出其外接圆包含插入的散点的三角形，删除影响三角形的公共边，将插入的散点与影响三角形的全部顶点连接起来，完成一个散点在Delaunay三角形链表中的插入。

63、根据Delaunay三角形优化准则对局部新形成的三角形进行优化，将优化后的三角形放入Delaunay三角形链表。

64、判断有序点集中是否还有未插入的散点，如果有则进行步骤62；如果没有则结束。

优选地，当引入地面的地性线和特征线时，对有序点集进行德劳内Delaunay三角剖分，获得Delaunay三角形网的方法包括以下步骤：

71、根据已有的地性线和特征线，形成控制边链表。

72、以控制边链表中的一个线段为基边，从有序点集中找出同基边的两端点的距离和最小的点，以该点为顶点，以基边为边，向外扩展一个三角形并放入三角形链表中。

73、判断控制边链表中是否还有未扩展的线段，如果有则进行步骤72；如果没有则进行步骤74。

74、依次将新形成的三角形的边作为基边，形成新的控制边链表，按照步骤72，对新的控制边链表中的所有的线段进行扩展。

75、判断新的控制边链表中是否还有未扩展的新形成的三角形，如果有则进行步骤74；如果没有则结束。

优选地，该方法还包括：当油藏对象非均匀分布和/或油藏对象的有序点集中的散点个数超过预定的散点阈值时，在对点集进行Delaunay三角剖分时，如果翻边后获得的三角形仍然不能满足Delaunay三角形的基本准则，则在该点集内部增加新的点作为辅助点，或增加新的点集作为辅助点集继续进行Delaunay三角剖分。

其中，在该点集内部增加新的点作为辅助点，或增加新的点集作为辅助点集采用随机增量构造法，随机增量构造法包括以下步骤：