一种基于压裂模拟结果的多段压裂水平井产能模拟方法与流程

1.本发明属于石油开采领域，更具体地，涉及一种基于压裂模拟结果的多段压裂水平井产能模拟方法、电子设备及介质。


背景技术：

2.页岩气是一种致密低渗的非常规气藏，其开发必须依赖多段压裂水平井才能实现商业化开发。为了降低开发风险，必须对多段压裂水平井的产能进行准确预测。数值模拟技术是预测多段压裂水平井产能的重要技术手段，在页岩气开发中广泛应用。准确预测多段压裂水平井产能的基础是准确刻画压后裂缝网络的分布。压裂模拟是表征压后改造裂缝分布的重要手段。
3.目前，常用的压裂模拟软件mangrove、fracman均可以模拟压裂后复杂的裂缝网格。在完成压裂模拟之后，常规的数值模拟软件会将压裂模拟结果进行粗化，然后采用双孔双渗或者双孔单渗模型进行模拟。粗化之后不能保留压裂模拟获得的裂缝真实形态，也不能考虑支撑剂分布造成裂缝内部开度和渗透率的不同，因此产能模拟的精度有所损失。嵌入式离散裂缝是近十年研发出来的新型裂缝模拟方法，与双孔双渗或者双孔单渗相比，嵌入式离散裂缝模型可以准确地表征不同走向的裂缝，模拟裂缝内的流动，模拟精度较双孔双渗或者双孔单渗有明显优势，但目前并没有嵌入式离散裂缝方法与压裂模拟软件模拟相结合的技术。
4.因此，期待发明一种多段压裂水平井产能模拟方法，能够有效解决现有技术中在数值模拟过程中无法完整保留裂缝的空间几何形态的问题。


技术实现要素：

5.本发明的目的是提出一种多段压裂水平井产能模拟方法，以解决现有技术中在数值模拟过程中无法完整保留裂缝的空间几何形态的问题。
6.为了实现上述目的，本发明提供一种基于压裂模拟结果的多段压裂水平井产能模拟方法，包括：
7.步骤1：通过压裂模拟获得裂缝片，通过所述裂缝片的子裂缝片确定所述裂缝片的顶点，并利用基质网格得到所述裂缝片的新子裂缝片；
8.步骤2：基于所述新子裂缝片，对所述裂缝片的物性进行重构，得到所述新子裂缝片的物性；
9.步骤3：基于所述新子裂缝片的物性，对所述多段压裂水平井进行产能模拟，得到产能模拟结果。
10.可选地，所述步骤1包括：
11.步骤11：通过所述压裂模拟获得裂缝片f，其中，所述裂缝片f为矩形，包括多个三角形的子裂缝片；
12.步骤12：对每个所述子裂缝片的三个顶点坐标进行编码；
13.步骤13：分别对所述步骤12中所有子裂缝片的顶点的x坐标、y坐标和z坐标进行排序，并将所述x坐标排序结果中的最大值和最小值记为x
max
和x
min
，将所述y坐标排序结果中的最大值和最小值记为y
max
和y
min
，将所述z坐标排序结果中的最大值和最小值记为z
max
和z
min
；
14.步骤14：基于所述x
max
、x
min
、y
max
、y
min
、z
max
和z
min
，对所述步骤12中所有子裂缝片的顶点的坐标进行检索，确定所述裂缝片f的顶点fp1、fp2、fp3和fp4；
15.步骤15：将所述裂缝片f的所述顶点fp1、fp2、fp3和fp4嵌入至所述基质网格中并与所述基质网格进行拓扑计算，得到多个新子裂缝片。
16.可选地，所述步骤14包括：
17.如果x
max
＝x
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4；
18.如果y
max
＝y
min
，将检索出包含坐标x
max
、x
min
、z
max
和z
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4；
19.如果x
max
≠x
min
且y
max
≠y
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4。
20.可选地，所述步骤15包括：
21.将所述裂缝片f的所述顶点fp1、fp2、fp3和fp4嵌入到所述基质网格中并与所述基质网格进行拓扑计算；
22.基于拓扑计算结果，利用所述基质网格将所述裂缝片f切割为多个所述新子裂缝片。
23.可选地，所述步骤2包括：
24.获取每个所述子裂缝片fi的三角网格物性ai；其中，i＝1，2，
…
，n且n为所述裂缝片f包括所述子裂缝片的个数；
25.计算每个所述新子裂缝片f
′j的中心到所述子裂缝片f1、f2、
…
、fn的中心的距离d
j1
、d
j2
、
…
、d
jn
，其中，j＝1，2，
…
，m且m为所述新子裂缝片的个数；
26.基于所述三角网格物性a1、a2、
…
、an以及每组所述距离d
j1
、d
j2
、
…
、d
jn
，通过公式计算每个所述新子裂缝片f
′j的物性a
′j。
27.可选地，所述公式为
[0028][0029]
可选地，所述步骤3包括：
[0030]
计算每个所述新子裂缝片的面积；
[0031]
基于每个所述新子裂缝片的面积和对应的所述新子裂缝片的物性，计算所述基质网格与每个所述新子裂缝片之间的传导率；
[0032]
基于所述基质网格与每个所述新子裂缝片之间的传导率，分别计算基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导并构建jacobian矩阵；
[0033]
采用高斯消去法对所述jacobian矩阵进行求解，得到所述产能模拟结果。
[0034]
一种电子设备，所述电子设备包括：
[0035]
存储器，存储有可执行指令；
[0036]
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0037]
一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0038]
本发明的有益效果在于：
[0039]
本发明首先通过压裂模拟获得裂缝片，通过裂缝片的子裂缝片确定裂缝片的顶点，并利用基质网格得到裂缝片的新子裂缝片，其次，基于新子裂缝片，对裂缝片的物性进行重构，得到新子裂缝片的物性，最后基于新子裂缝片的物性，对多段压裂水平井进行产能模拟，得到产能模拟结果；本发明应用压裂模拟的结果来构建嵌入式离散裂缝，通过这种方法完整保留了裂缝的空间几何形态，然后进行多段压裂水平井的产能模拟，从而更加准确的模拟页岩气多段压裂水平井的产能。
[0040]
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
[0041]
通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。
[0042]
图1示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的流程图。
[0043]
图2示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的裂缝片的结构示意图。
[0044]
图3示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的裂缝片嵌入至基质网格的结构示意图。
[0045]
图4示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的基质-裂缝jacobian矩阵示意图。
[0046]
图5示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的压裂模拟获得的压裂缝网分布示意图。
[0047]
图6示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的嵌入式离散裂缝的模型示意图。
[0048]
图7示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的嵌入式离散裂缝的多段压裂水平井的预测产能与实测产能的对比图。
具体实施方式
[0049]
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
[0050]
根据本发明的一种基于压裂模拟结果的多段压裂水平井产能模拟方法，包括：
[0051]
步骤1：通过压裂模拟获得裂缝片，通过裂缝片的子裂缝片确定裂缝片的顶点，并
利用基质网格得到裂缝片的新子裂缝片；
[0052]
步骤2：基于新子裂缝片，对裂缝片的物性进行重构，得到新子裂缝片的物性；
[0053]
步骤3：基于新子裂缝片的物性，对多段压裂水平井进行产能模拟，得到产能模拟结果。
[0054]
具体地，本发明首先通过压裂模拟获得裂缝片，通过裂缝片的子裂缝片确定裂缝片的顶点，并利用基质网格得到裂缝片的新子裂缝片，其次，基于新子裂缝片，对裂缝片的物性进行重构，得到新子裂缝片的物性，最后基于新子裂缝片的物性，对多段压裂水平井进行产能模拟，得到产能模拟结果；本发明应用压裂模拟的结果来构建嵌入式离散裂缝，通过这种方法完整保留了裂缝的空间几何形态，然后进行多段压裂水平井的产能模拟，从而更加准确的模拟页岩气多段压裂水平井的产能。
[0055]
进一步地，基质网格为现有技术，本领域技术人员很容易获得，计算基质网格12条棱与裂缝片f的交点，将所有交点连接成封闭的多边形即为该基质网格切割裂缝片获得的子裂缝片，采用此方法依次计算获得所有网格与裂缝片f得到的新子裂缝片。
[0056]
在一个示例中，步骤1包括：
[0057]
步骤11：通过压裂模拟获得裂缝片f，其中，裂缝片f为矩形，包括多个三角形的子裂缝片；
[0058]
步骤12：对每个子裂缝片的三个顶点坐标进行编码；
[0059]
步骤13：分别对步骤12中所有子裂缝片的顶点的x坐标、y坐标和z坐标进行排序，并将x坐标排序结果中的最大值和最小值记为x
max
和x
min
，将y坐标排序结果中的最大值和最小值记为y
max
和y
min
，将z坐标排序结果中的最大值和最小值记为z
max
和z
min
；
[0060]
步骤14：基于x
max
、x
min
、y
max
、y
min
、z
max
和z
min
，对步骤12中所有子裂缝片的顶点的坐标进行检索，确定裂缝片f的顶点fp1、fp2、fp3和fp4；
[0061]
步骤15：将裂缝片f的顶点fp1、fp2、fp3和fp4嵌入至基质网格中并与基质网格进行拓扑计算，得到多个新子裂缝片。
[0062]
在一个示例中，步骤14包括：
[0063]
如果x
max
＝x
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4；
[0064]
如果y
max
＝y
min
，将检索出包含坐标x
max
、x
min
、z
max
和z
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4；
[0065]
如果x
max
≠x
min
且y
max
≠y
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4。
[0066]
在一个示例中，步骤15包括：
[0067]
将裂缝片f的顶点fp1、fp2、fp3和fp4嵌入到基质网格中并与基质网格进行拓扑计算；
[0068]
基于拓扑计算结果，利用基质网格将裂缝片f切割为多个新子裂缝片。
[0069]
在一个示例中，步骤2包括：
[0070]
获取每个子裂缝片fi的三角网格物性ai；其中，i＝1，2，
…
，n且n为裂缝片f包括子裂缝片的个数；
[0071]
计算每个新子裂缝片f
′j的中心到子裂缝片f1、f2、
…
、fn的中心的距离d
j1
、d
j2
、
…
、djn
，其中，j＝1，2，
…
，m且m为新子裂缝片的个数；
[0072]
基于三角网格物性a1、a2、
…
、an以及每组距离d
j1
、d
j2
、
…
、d
jn
，通过公式计算每个新子裂缝片f
′j的物性a
′j。
[0073]
在一个示例中，公式为
[0074][0075]
在一个示例中，步骤3包括：
[0076]
计算每个新子裂缝片的面积；
[0077]
基于每个新子裂缝片的面积和对应的新子裂缝片的物性，计算基质网格与每个新子裂缝片之间的传导率；
[0078]
基于基质网格与每个新子裂缝片之间的传导率，分别计算基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导并构建jacobian矩阵；
[0079]
采用高斯消去法对jacobian矩阵进行求解，得到产能模拟结果。
[0080]
具体地，基质网格与新子裂缝片的传导率t的公式为
[0081][0082]
其中，k
nnc
为基质渗透率和新子裂缝片渗透率的调和平均数，a
nnc
为基质渗透率和新子裂缝片的交接面面积，d
nnc
为基质网格中心和新子裂缝片形心间距离。
[0083]
进一步地，根据渗流物理中两点间流动的达西公式计算基质间流动项、基质-裂缝间流动项和裂缝-裂缝间流动项，并将上述流动项分别对压力和含气饱和度求偏导获得基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导，将上述基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导代入质量守恒方程构建jacobian矩阵。
[0084]
一种电子设备，电子设备包括：
[0085]
存储器，存储有可执行指令；
[0086]
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0087]
一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0088]
实施例1
[0089]
如图1所示，一种基于压裂模拟结果的多段压裂水平井产能模拟方法，包括：
[0090]
步骤1：通过压裂模拟获得裂缝片，通过裂缝片的子裂缝片确定裂缝片的顶点，并利用基质网格得到裂缝片的新子裂缝片；
[0091]
步骤2：基于新子裂缝片，对裂缝片的物性进行重构，得到新子裂缝片的物性；
[0092]
步骤3：基于新子裂缝片的物性，对多段压裂水平井进行产能模拟，得到产能模拟结果。
[0093]
具体实施方式为
[0094]
步骤1包括：
[0095]
步骤11：如图1所示，通过压裂模拟获得裂缝片f，其中，裂缝片f为矩形，包括四个三角形的子裂缝片f1、f2、f3和f4；
[0096]
步骤12：对每个子裂缝片的三个顶点坐标进行编码，其中，所述子裂缝片f1的三个顶点坐标按逆时针方向编为f1p1(x
f1p1
,y
f1p1
,z
f1p1
)、f1p2(x
f1p2
,y
f1p2
,z
f1p2
)和f1p3(x
f1p3
,y
f1p3
,z
f1p3
)，所述子裂缝片f2的三个顶点坐标按逆时针方向编为f2p1(x
f2p1
,y
f2p1
,z
f2p1
)、f2p2(x
f2p2
,y
f2p2
,z
f2p2
)和f2p3(x
f2p3
,y
f2p3
,z
f2p3
)，所述子裂缝片f3的三个顶点坐标按逆时针方向编为f3p1(x
f3p1
,y
f3p1
,z
f3p1
)、f3p2(x
f3p2
,y
f3p2
,z
f3p2
)和f1p3(x
f3p3
,y
f3p3
,z
f3p3
)，所述子裂缝片f4的三个顶点坐标按逆时针方向编为f4p1(x
f4p1
,y
f4p1
,z
f4p1
)、f4p2(x
f4p2
,y
f4p2
,z
f4p2
)和f4p3(x
f4p3
,y
f4p3
,z
f4p3
)；
[0097]
步骤13：分别对步骤12中所有子裂缝片的顶点的x坐标、y坐标和z坐标进行排序，并将x坐标排序结果中的最大值和最小值记为x
max
和x
min
，将y坐标排序结果中的最大值和最小值记为y
max
和y
min
，将z坐标排序结果中的最大值和最小值记为z
max
和z
min
；
[0098]
步骤14：基于x
max
、x
min
、y
max
、y
min
、z
max
和z
min
，对步骤12中所有子裂缝片的顶点的坐标进行检索，确定裂缝片f的顶点fp1、fp2、fp3和fp4；
[0099]
步骤15：如图3所示，将裂缝片f的顶点fp1、fp2、fp3和fp4嵌入至基质网格中并与基质网格进行拓扑计算，得到新子裂缝片f
′1和f
′2。
[0100]
其中，步骤14包括：
[0101]
如果x
max
＝x
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4；
[0102]
如果y
max
＝y
min
，将检索出包含坐标x
max
、x
min
、z
max
和z
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4；
[0103]
如果x
max
≠x
min
且y
max
≠y
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为裂缝片f的顶点fp1、fp2、fp3和fp4。
[0104]
其中，步骤15包括：
[0105]
将裂缝片f的顶点fp1、fp2、fp3和fp4嵌入到基质网格中并与基质网格进行拓扑计算；
[0106]
基于拓扑计算结果，利用基质网格将裂缝片f切割为新子裂缝片f
′1和f
′2。
[0107]
步骤2包括：
[0108]
分别获取子裂缝片f1的三角网格物性a1、子裂缝片f2的三角网格物性a2、子裂缝片f3的三角网格物性a3、子裂缝片f4的三角网格物性a4；
[0109]
分别计算新子裂缝片f
′1的中心到子裂缝片f1、f2、f3和f4的中心的距离d
11
、d
12
、d
13
和d
14
，并分别计算新子裂缝片f
′2的中心到子裂缝片f1、f2、f3和f4的中心的距离d
21
、d
22
、d
23
和d
24
；
[0110]
基于三角网格物性a1、a2、a3和a4以及距离d
11
、d
12
、d
13
和d
14
，通过第一公式计算新子裂缝片f
′1的物性a
′1，并基于三角网格物性a1、a2、a3和a4以及距离d
21
、d
22
、d
23
和d
24
，通过第二公式计算新子裂缝片f
′2的物性a
′2。
[0111]
其中，第一公式为
[0112][0113]
第二公式为
[0114][0115]
步骤3包括：
[0116]
计算每个新子裂缝片的面积；
[0117]
基于每个新子裂缝片的面积和对应的新子裂缝片的物性，计算基质网格与每个新子裂缝片之间的传导率；
[0118]
基于基质网格与每个新子裂缝片之间的传导率，分别计算基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导并构建jacobian矩阵，如图4所示；
[0119]
采用高斯消去法对jacobian矩阵进行求解，得到产能模拟结果。
[0120]
具体地，图5示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的压裂模拟获得的压裂缝网分布示意图，由图5可知，压裂缝网分布非常复杂，图6示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的嵌入式离散裂缝的模型示意图，通过图5和图6对比可知，本发明可以准确地反映出裂缝的分布及裂缝的属性，从而保证了数值模拟的精度；图7示出了根据本发明的一个实施例的一种基于压裂模拟结果的多段压裂水平井产能模拟方法的嵌入式离散裂缝的多段压裂水平井的预测产能与实测产能的对比图，由图7可知，本发明可以准确地预测页岩气多段压裂水平井的产能。
[0121]
实施例2
[0122]
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0123]
根据本公开实施例的电子设备包括存储器和处理器。
[0124]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0125]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0126]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0127]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0128]
实施例3
[0129]
本公开提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述基于压裂模拟结果的多段压裂水平井产能模拟方法。
[0130]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指
令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0131]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0132]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

技术特征：
1.一种基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，包括：步骤1：通过压裂模拟获得裂缝片，通过所述裂缝片的子裂缝片确定所述裂缝片的顶点，并利用基质网格得到所述裂缝片的新子裂缝片；步骤2：基于所述新子裂缝片，对所述裂缝片的物性进行重构，得到所述新子裂缝片的物性；步骤3：基于所述新子裂缝片的物性，对所述多段压裂水平井进行产能模拟，得到产能模拟结果。2.根据权利要求1所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述步骤1包括：步骤11：通过所述压裂模拟获得裂缝片f，其中，所述裂缝片f为矩形，包括多个三角形的子裂缝片；步骤12：对每个所述子裂缝片的三个顶点坐标进行编码；步骤13：分别对所述步骤12中所有子裂缝片的顶点的x坐标、y坐标和z坐标进行排序，并将所述x坐标排序结果中的最大值和最小值记为x
max
和x
min
，将所述y坐标排序结果中的最大值和最小值记为y
max
和y
min
，将所述z坐标排序结果中的最大值和最小值记为z
max
和z
min
；步骤14：基于所述x
max
、x
min
、y
max
、y
min
、z
max
和z
min
，对所述步骤12中所有子裂缝片的顶点的坐标进行检索，确定所述裂缝片f的顶点fp1、fp2、fp3和fp4；步骤15：将所述裂缝片f的所述顶点fp1、fp2、fp3和fp4嵌入至所述基质网格中并与所述基质网格进行拓扑计算，得到多个新子裂缝片。3.根据权利要求2所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述步骤14包括：如果x
max
＝x
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4；如果y
max
＝y
min
，将检索出包含坐标x
max
、x
min
、z
max
和z
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4；如果x
max
≠x
min
且y
max
≠y
min
，将检索出包含坐标x
max
、x
min
、y
max
和y
min
的四个顶点作为所述裂缝片f的顶点fp1、fp2、fp3和fp4。4.根据权利要求2所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述步骤15包括：将所述裂缝片f的所述顶点fp1、fp2、fp3和fp4嵌入到所述基质网格中并与所述基质网格进行拓扑计算；基于拓扑计算结果，利用所述基质网格将所述裂缝片f切割为多个所述新子裂缝片。5.根据权利要求2所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述步骤2包括：获取每个所述子裂缝片f
i
的三角网格物性a
i
；其中，i＝1，2，
…
，n且n为所述裂缝片f包
括所述子裂缝片的个数；计算每个所述新子裂缝片f
′
j
的中心到所述子裂缝片f1、f2、
…
、fn的中心的距离d
j1
、d
j2
、
…
、d
jn
，其中，j＝1，2，
…
，m且m为所述新子裂缝片的个数；基于所述三角网格物性a1、a2、
…
、a
n
以及每组所述距离d
j1
、d
j2
、
…
、d
jn
，通过公式计算每个所述新子裂缝片f
′
j
的物性a
′
j
。6.根据权利要求5所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述公式为7.根据权利要求2所述的基于压裂模拟结果的多段压裂水平井产能模拟方法，其特征在于，所述步骤3包括：计算每个所述新子裂缝片的面积；基于每个所述新子裂缝片的面积和对应的所述新子裂缝片的物性，计算所述基质网格与每个所述新子裂缝片之间的传导率；基于所述基质网格与每个所述新子裂缝片之间的传导率，分别计算基质间流动偏导、基质-裂缝间流动偏导和裂缝-裂缝间流动偏导并构建jacobian矩阵；采用高斯消去法对所述jacobian矩阵进行求解，得到所述产能模拟结果。8.一种电子设备，其特征在于，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现根据权利要求1-7中任一项所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。9.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的基于压裂模拟结果的多段压裂水平井产能模拟方法。

技术总结
本发明公开了一种基于压裂模拟结果的多段压裂水平井产能模拟方法，包括：步骤1：通过压裂模拟获得裂缝片，通过裂缝片的子裂缝片确定裂缝片的顶点，并利用基质网格得到裂缝片的新子裂缝片；步骤2：基于新子裂缝片，对裂缝片的物性进行重构，得到新子裂缝片的物性；步骤3：基于新子裂缝片的物性，对多段压裂水平井进行产能模拟，得到产能模拟结果。本发明应用压裂模拟的结果来构建嵌入式离散裂缝，通过这种方法完整保留了裂缝的空间几何形态，然后进行多段压裂水平井的产能模拟，从而更加准确的模拟页岩气多段压裂水平井的产能。拟页岩气多段压裂水平井的产能。拟页岩气多段压裂水平井的产能。


技术研发人员：戴城 方思冬 刘华 王卫红 胡小虎 王妍妍
受保护的技术使用者：中国石油化工股份有限公司石油勘探开发研究院
技术研发日：2022.03.25
技术公布日：2023/10/8