复杂断块油藏产量劈分方法与流程

本发明涉及油藏开放，尤其涉及一种复杂断块油藏产量劈分方法。

背景技术：

1、复杂断块的油藏一般表现为油层多、井段长、层间差异大、油水关系复杂、油砂体面积小、物性较差、平面非均质性强等地质特征；地质条件、措施部署、工程工艺等因素也会影响油藏层间、平面动用情况，使得差异更加明显。而对于油藏开发评价、小层动用状况分析、合理开发调整政策制定等都需要合理、准确的小层油藏产量劈分的数据支撑。

2、产量劈分一直都是油田开发中非常困扰技术人员的重要问题。目前，一般采用的是数值模拟方法和静态产量劈分方法。对于数值模拟方法而言，虽然其模拟结果的可信度较高，但是，该方法需要前期的地质建模做支撑，而且该方法涉及的数据量和计算量都较为巨大、拟合周期长、操作难度高，此外，该方法对软硬件的要求也很高，需要强大的人力、物力为基础，资源投入很大。对于静态产量劈分方法而言，其一般是运用油水井的静态地质参数来对动态数据进行劈分，其考虑的影响因素较为单一，准确度不高，并且，使用静态数据来劈分动态数据本身就极为不合理，特别是在含水率的劈分问题上。

技术实现思路

1、本发明实施例提供一种复杂断块油藏产量劈分方法，旨在解决或者部分解决如何合理、准确地计算复杂断块油藏的产量劈分的问题。

2、为了解决上述技术问题，本发明是这样实现的：

3、第一方面，本发明提供了一种复杂断块油藏产量劈分方法，所述方法包括：

4、注采井网格划分：对含油面积区域的注采关系和网格进行划分，并获取网格对应的属性；所述属性包括孔隙度、渗透率、孔隙体积、储量和含油饱和度；

5、液量权重系数计算：计算对应井段的纵向权重系数；

6、劈分液量、注水量计算各方向上的注入倍数：基于液量权重系数计算步骤所得的纵向权重系数，确定每个注水井在不同小层、不同方向上的劈分量和注入倍数；

7、计算各方向上含水率：查询注入倍数与含水率的关系表，获取含水率和对应的平均含水饱和度，并引入拟合系数拟合含水率；

8、统计单井产水量：计算单个采油井的总产水量，并与实际产水量进行对比，对拟合系数进行修正，并基于修正后的拟合系数，在注入倍数与含水率的关系表中重新获取注入倍数，再次计算总产水量，直至总产水量与实际产水量的误差值达到第一阈值；

9、计算含水饱和度：基于统计单井产水量步骤所得的拟合系数，重新计算平均含水饱和度，并对应更新纵向权重系数；

10、输出各方向一定时间内产油量、产水量和注水量：根据计算含水饱和度步骤所得的平均含水饱和度，以及含水率和注入倍数，确定注水井和每个采油井组成的油水井在注采方向上所对应的产注劈分量，并计算在注采方向上的累积产油量和累积注水量。

11、可选地，在注采井网格划分步骤中，还包括将含油面积区域按照逐层且间隔一定时间的顺序对注采关系和网格进行划分。

12、可选地，对注采关系和网格进行划分的步骤具体为：

13、在规则注采井网区域中，以注水井为中心，并以每相邻两口一线采油井的连线中点与该注水井的连线作为划分线，进行注采关系和网格的划分；

14、在不规则注采井网区域中，注水井向外扩展三分之一井距，采油井向外扩展二分之一到三分之一井距，并以每相邻两口一线采油井的连线中点与对应的注水井的连线作为划分线，进行注采关系和网格的划分；或者，以砂体尖灭边界、有效储层边界、封闭断层边界或区域边界中的一者或多者作为划分边界，进行注采关系和网格的划分。

15、可选地，当含油面积区域内存在补孔、封层和开关井措施中的任一者时，对应按照一定的时间间隔进行数据更新；

16、其中，压裂后的渗透率按照10倍原渗透率的数值进行计算，酸化后的渗透率按照3倍原渗透率的数值进行计算。

17、可选地，在液量权重系数计算步骤中，

18、当产油井和注水井的全井段均有产吸剖面的测试数据时，将测试剖面的权重系数定义为对应井段的纵向权重系数，其中，产吸剖面为产液剖面和吸水剖面；

19、当无产吸剖面的测试数据时，对应井段的纵向权重系数按照下述公式进行计算：

20、

21、式中，wijk为注水井i与采油井j之间的连通层k层的权重系数；

22、k为渗透率；

23、h为层厚。

24、可选地，在液量权重系数计算步骤中，

25、当产油井和注水井的部分井段有产吸剖面的测试数据时，先按照无产吸剖面测试数据时的纵向权重系数计算方法计算有测试数据的产吸剖面对应的层位的纵向权重系数和无测试数据的产吸剖面对应的层位的纵向权重系数，并劈分有测试数据的产吸剖面和无测试数据的产吸剖面对应层段的液量；

26、然后将有测试数据的产吸剖面对应的井段，以及无测试数据的产吸剖面对应的井段按照测试数据的权重分别将对应液量劈分到相应的小层上。

27、可选地，劈分液量、注水量计算各方向上的注入倍数步骤具体为：

28、按照权重系数，按照下列计算式来确定每个注水井在各方向上的劈分量和注入倍数：

29、

30、

31、

32、

33、式中，为生产井i在采油井j的第k连通方向上的采液量的劈分量；

34、为采油井j的采液量；

35、wijk为注水井i与采油井j之间的连通层k层的权重系数；

36、为注入井i在生产井j的第k连通方向上的注水量的劈分量；

37、为注入井i的注水量；

38、为t时刻注入井i在生产井j的第k连通方向上的累积注水量；为单方向的累积注水量；

39、pv为注入井i在生产井j的第k连通方向上的空隙体积；

40、pvi为注入井i在生产井j第k连通方向上的注入倍数。

41、可选地，计算各方向上含水率步骤具体为：

42、根据劈分液量、注水量计算各方向上的注入倍数步骤计算得到的注入倍数，查询注入倍数与含水率的关系表，获取含水率和对应的平均含水饱和度，引入拟合系数∝，对注入倍数进行调整，拟合含水率采用下述公式进行计算：

43、fw＝f1(pvi,∝)；

44、

45、式中，fw为含水率；

46、为平均含水饱和度；

47、pvi为注入倍数；

48、f1和f2为调整过程中第一次和第二次查表获取的含水率；

49、∝为拟合系数，初值取1。

50、可选地，所述计算含水饱和度步骤具体为：

51、基于所得拟合系数，重新计算所述平均含水饱和度，并根据下式更新所述权重系数：

52、

53、式中，wijk为注入井i在生产井j的第k连通方向上的液量权重系数；

54、k为井间的平均渗透率，单位为md；

55、h为井间的平均厚度，单位为m；

56、pi为注水井的井底流压，单位为mpa；

57、pj为生产井的井底流压，单位为mpa；

58、l为注采井间距离，单位为m；

59、为平均含水饱和度；

60、为油相相渗，为的函数；

61、为水相相渗，为的函数。

62、可选地，所述方法在计算含水饱和度步骤和输出各方向一定时间内产油量、产水量和注水量步骤之间还包括迭代步骤：

63、重复液量权重系数计算步骤至计算含水饱和度步骤，进行数值迭代，直到纵向权重系数的变化值达到预期值或达到最大迭代步数；

64、其中，预期值和最大迭代步数均为预设值。

65、第二方面，本发明实施例提供了一种电子设备，电子设备包括：

66、至少一个处理器；以及，

67、与至少一个处理器通信连接的存储器；其中，

68、存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器能够执行本发明实施例的第一方面提供的方法步骤。

69、第三方面，本发明实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时可以实现本发明实施例第一方面提供的方法步骤。

70、本发明实施例包括以下优点：

71、通过上述技术方案，本发明的方法结合不同时刻下的不同注采对应关系可以对单井单层单方向上的纵向权重系数进行合理且准确的计算，能够将一定时间内的产注量劈分到各个方向上，最后累加得到的各个方向上的累积产油量和累积产水量能够更加准确，并能够契合于实际情况。

72、在本公开的方法中，将数值模拟图版、相渗曲线与含水饱和度挂钩，力求将注采方向上的液量、含水率进行合理劈分。本发明将小层含水饱和度作为一项重要影响因素，通过相渗曲线将油相相对渗透率以及水相相对渗透率与小层含水饱和度相结合，将油水相对渗透率看作为小层含水饱和度的函数，通过数值模拟图版将其建立函数关系。综合考虑小层渗透率、有效厚度、小层含水饱和度、生产压差、油水井距、地层对油水井的贡献、油水粘度等参数对产量劈分的影响以及相关增产资料，并结合剖面测试资料对劈分结果进行修正，应用本发明的劈分方法对油藏产量劈分，可以实现劈分精确度高、适用性强，可满足实际应用需要。