1.一种多尺度堵剂封堵后油井产能的预测方法,其特征在于,所述预测方法包括:将封堵后的高含水压裂直井划分为多个渗流区域,所述多个渗流区域包括第一渗流区域和第二渗流区域,所述第一渗流区域包括区域1、区域3和区域4,所述区域1为主裂缝区域,所述区域3为所述主裂缝区域的最远端至半井距构成的区域,所述区域4为所述第二渗流区域至排距构成的区域;确定第二渗流区域封堵后在当前时间步的渗透率 和孔隙度 ;根据边界条件确定每个渗流区域的基本渗流关系式,以确定在当前时间步下每个渗流区域在拉普拉斯空间中的油相拟压力值 和产量;确定下一个时间步的每个渗流区域在封堵后的初始平均压力值 ;根据下一个时间步的初始平均压力值 确定出下一个时间步的封堵后的高含水压裂直井的含油饱和度 ,并确定出每个第一渗流区域在下一个时间步的渗透率 ;根据下一个时间步的含油饱和度 和所述第一渗流区域中的主裂缝区域在下一个时间步的油相拟压力值 预测封堵后的油井在下一个时间步的油井产量;其中,所述根据边界条件确定每个渗流区域的基本渗流关系式,以确定在当前时间步下每个渗流区域在拉普拉斯空间中的油相拟压力值 和产量包括:根据边界条件确定每个渗流区域的基本渗流关系式,将第一关系式和第二关系式代入每个基本渗流关系式进行化简,以确定在当前时间步下每个渗流区域在拉普拉斯空间中的油相拟压力值 和产量;所述第一关系式为公式(3): (3) ;其中, 为渗流区域n封堵后在当前时间步的孔隙度, 为与每个渗流区域对应的封堵前在当前时间步的孔隙度, 为岩石的压缩系数, 为渗流区域n封堵后在当前时间步的压力值, 为油藏在当前时间步的初始压力值;所述第二关系式为公式(4): (4);其中, , 为渗流区域n在当前时间步的油相的相对渗透率,n=1,2,3,4, 为渗流区域n封堵后在当前时间步的渗透率, 为渗流区域n的油相拟压力, 为原油的粘度, 为原油的体积系数;其中,所述根据边界条件确定每个渗流区域的基本渗流关系式包括:所述区域4的基本渗流方程为公式(5): (5);所述区域3的基本渗流方程为公式(6): (6);所述第二渗流区域的基本渗流方程为公式(7): (7);所述区域1的基本渗流方程为公式(8): (8);其中, 为所述区域2的基质中的压力值, 为所述区域2的裂缝中的压力值, 为所述区域3压力值,D为分形维数, 为连通指数, 为垂直于主缝方向堵剂封堵长度, 为距水力裂缝距离的变量, 为原油的密度, 为所述主裂缝区域的宽, 为所述主裂缝区域的半长。
2.根据权利要求1所述的多尺度堵剂封堵后油井产能的预测方法,其特征在于,所述确定第二渗流区域封堵后在当前时间步的渗透率 和孔隙度 包括,根据公式(1)确定所述第二渗流区域封堵后在当前时间步的渗透率 : (1); 为所述第二渗流区域封堵前在当前时间步的渗透率, 为所述第二渗流区域封堵后在当前时间步的渗透率, 为垂直于主缝方向堵剂封堵长度,即堵剂运移距离, 为距水力裂缝距离的变量,D为分形维数, 为连通指数。
3.根据权利要求1所述的多尺度堵剂封堵后油井产能的预测方法,其特征在于,所述确定第二渗流区域封堵后在当前时间步的渗透率 和孔隙度 包括,根据公式(2)确定所述第二渗流区域封堵后在当前时间步的孔隙度 : (2);其中, 为垂直于主缝方向堵剂封堵长度,即堵剂运移距离, 为距水力裂缝距离的变量,D为分形维数, 为与所述第二渗流区域封堵前在当前时间步的孔隙度, 为所述第二渗流区域封堵后在当前时间步的孔隙度。
4.根据权利要求1所述的多尺度堵剂封堵后油井产能的预测方法,其特征在于,所述确定下一个时间步的封堵后的高含水压裂直井的含油饱和度 包括,根据公式(9)计算出下一个时间步的封堵后的高含水压裂直井的含油饱和度 : (9)其中, ,N为原始地层储量, 为采出油量, 为半井距, 为束缚水饱和度, 为所述主裂缝区域的宽, 为所述主裂缝区域的半长, 为排距。
5.根据权利要求1所述的多尺度堵剂封堵后油井产能的预测方法,其特征在于,所述根据下一个时间步的含油饱和度 和所述第一渗流区域中的主裂缝区域在下一个时间步的油相拟压力值 预测封堵后的油井在下一个时间步的油井产量包括:根据公式(10)计算出所述封堵后的油井在下一个时间步在拉普拉斯空间内的油井产量 :对所述油井产量 进行stehfest数值反演,以得到所述封堵后的油井在下一个时间步的油井产量q; (10);其中, 为在拉普拉斯空间内所述第一渗流区域中的主裂缝区域在下一个时间步的油相拟压力值,s为拉普拉斯变换中的变量。
6. 根据权利要求5所述的多尺度堵剂封堵后油井产能的预测方法,其特征在于,根据公式(11)计算出 : (11)其中, , 为等式代号, 为井筒半径, 为系数,利用边界条件求出, 为第一类虚宗量Bessel函数,可直接调用, 为系数,利用边界条件求出, 为第二类虚宗量Bessel函数,可直接调用, 为下一个时间步的封堵后的高含水压裂直井的含油饱和度, 为所述区域1在下一个时间步的渗透率, 为所述第二渗流区域在当前时间步的油相的相对渗透率, 为所述主裂缝区域的宽; 的表达式为以下公式(12)所示: (12)其中, , ,D为分形维数, 为连通指数, 和 是直接通过边界条件化简得到的参数, 、 以及 分别为第一类、第二类虚宗量Bessel函数,n、 和 均为等式代号; 、 、n的表达式分别如以下公式(13)、公式(14)和公式(15)所示: (13); (14); (15);其中, , 为变换化简后的中间量, 为原油的粘度, 为垂直于主缝方向堵剂封堵长度,即堵剂运移距离, 为所述主裂缝区域的半长, 为半井距, , 为岩石的压缩系数, 为渗流区域n封堵后在当前时间步的孔隙度。
7.一种处理器,其特征在于,被配置成执行根据权利要求1至6中任意一项所述的多尺度堵剂封堵后油井产能的预测方法。
8.一种多尺度堵剂封堵后油井产能的预测装置,其特征在于,包括根据权利要求7所述的处理器。