1.一种剩余油分散程度表征方法,其特征在于,包括:步骤S1:建立待测区域的油藏数值模拟模型,之后进行生产历史拟合;步骤S2:根据所述步骤S1所获取的参数,利用剩余油丰度公式计算不同时间步的剩余油丰度;步骤S3:利用邻域边界追踪算法识别所述步骤S2中高丰度的连续油相,高丰度的连续油相即剩余油斑块,并获取各所述剩余油斑块的面积、周长及所述剩余油斑块的总块数,并由此计算平均斑块面积MPS、斑块密度PD以及平均形状指数MSI;步骤S4:将所述步骤S3中得到的不同时间步的所述平均斑块面积MPS、所述斑块密度PD和所述平均形状指数MSI组成指标矩阵,然后通过熵权法计算得到不同时间步上的剩余油分散度;其中,在所述步骤S3中,所述平均形状指数MSI的计算公式为:其中,E为剩余油斑块总周长,单位为:m;A为剩余油总面积,单位为:m 2 ;在所述熵权法计算的过程中,需要对正向指标所述斑块密度PD和所述平均形状指数MSI归一化处理以及对负向指标所述平均斑块面积MPS归一化处理。
2.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S1中,通过获取待测区域的地质属性参数、流体分布参数、注采井生产参数和注采控制参数,以建立所述油藏数值模拟模型。
3.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S1中,所述生产历史拟合包括区块历史拟合和单井历史拟合。
4.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S2中,所述剩余油丰度公式为:剩余油储量丰度公式: 或剩余油可采储量丰度公式: 其中,J o1 为剩余油储量丰度,单位为:10 4 t/km 2 ;J o2 为剩余油可采储量丰度,单位为:10 4 t/km 2 ;h为储层厚度,单位为:m;φ为孔隙度;S o 为含油饱和度;S or 为残余油饱和度;ρ o 为原油密度,单位为:g/cm 3 ;B o 为原油体积系数。
5.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S3中,所述平均斑块面积MPS的计算公式为:其中,A为剩余油总面积,单位为:m 2 ;N为剩余油斑块块数。
6.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S3中,所述斑块密度PD的计算公式为:其中,A为剩余油总面积,单位为:m 2 ;N为剩余油斑块块数。
7.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S4中,所述熵权法的赋权的计算过程包括以下步骤:假设有m个时间点,n个评价指标,x ij 是第i个时间点第j个指标的值(i=1,…,m;j=1,2,…,n);步骤S41:指标归一化处理,其中,正向指标所述斑块密度PD和所述平均形状指数MSI的归一化公式为 负向指标所述平均斑块面积MPS的归一化公式为 步骤S42:计算第j项指标下第i个样本值占该指标的比重,公式为 步骤S43:计算第j项指标的熵值,公式为 其中,k=1/lnm,当p ij =0时, 步骤S44:计算信息熵冗余度,公式为d j =1-e j ;步骤S45:计算各指标的权重,公式为 步骤S46:计算剩余油分散度,公式为 。
8.根据权利要求1所述的剩余油分散程度表征方法,其特征在于,在所述步骤S1中,所述油藏数值模拟模型包括网格属性、油藏属性和流体属性。
9.根据权利要求8所述的剩余油分散程度表征方法,其特征在于,所述网格属性包括网格数目和网格大小;所述油藏属性包括油藏深度、孔隙度和原始地层压力;所述流体属性包括原油密度、原油粘度和体积系数。