1.一种致密储层三维空间孔喉定量分级评价方法,其特征在于,所述方法包括:
1)利用avizo软件构建数字岩心模型并提取孔隙数据和喉道数据;
2)根据配位数和所述孔隙数据将孔隙划分为有效孔隙和无效孔隙;
3)将有效孔隙、无效孔隙和喉道数据分别进行BOX-COX的变换;
4)分别进行无量纲化处理;
5)处理后的有效孔隙、无效孔隙和喉道数据进行K-均值聚类分析,并将分类结果分别存储到数据库1、数据库2和数据库3中;
6)将数据库1、数据库2和数据库3的数据库整合分析,并相结合对致密储层内部孔喉进行细分。
2.根据权利要求1所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,所述步骤3)中BOX-COX变换具体包括:分别对数据进行BOX-COX的变换,BOX-COX公式为:其中,y为反应变量,c为了使y+c>0的常数,λ为变换参数。
3.根据权利要求2所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,利用极大似然估计法估算λ,极大然估计法的公式为:(λ)
式中,y 代表BOX-COX的公式, 代表y1,y2...yn的几何平均值,n为观测值,I为系数因子,X为特征向量。
4.根据权利要求1所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,步骤4)中无量纲化处理的公式为:式中,xi为准备处理的数据; 为对应数据列的均;σ为标准差;yi为处理后的数据。
5.根据权利要求1所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,在步骤4)之后步骤5)之前还包括对处理后的喉道数据进行灰关联分析,具体包括:分别选出每一个孔隙连通多个喉道的喉道数据,将孔隙数据作为母序列用X0(k)表示,将该孔隙对应的多个喉道数据作为比较序列用Xi(k)表示,进行灰关联分析,关联度值最大的作为表征该孔隙的喉道数据,灰关联分析公式为:式中,Δi(k)=|X0(k)-Xi(k)|代表比较序列与参考序列对应元素的绝对差值,y为分辨系数。
6.根据权利要求5所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,所述步骤5)中K-均值聚类分析具体包括以下步骤:
5.1根据储层特点选择K个初始点,将每类样品集合记为C(l);
5.2计算出某一个孔隙数据或者喉道数据到K类重心的欧式距离;
5.3将样品移动到与K类中心距离最近的C(l);
5.4若此数据类别发生改变,则需要重新计算各类重心;
5.5对选取孔隙或者喉道的数据继续按照这样的方式迭代下去直到达到选择的迭代次数或者迭代稳定为止。
7.根据权利要求6所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,所述步骤5.2中K类中心的距离公式为:式中,d(xi,xj)为欧氏距离;xi为第i个样品;xir为第i个样品的第r个特征参数。
8.根据权利要求6所述一种致密储层三维空间孔喉定量分级评价方法,其特征在于,所述步骤5.3中重心计算公式为:式中,l=1,2,...,K。