[发明专利]基于测井曲线的岩性识别方法及装置

权利要求书

1.一种基于测井曲线的岩性识别方法，其特征在于，包括：

获取目标井的多条测井曲线；

标准化所述测井曲线；

对标准化后的测井曲线进行主成分分析，以获取主成分测井数据，所述对标准化后的测井曲线进行主成分分析包括：获取初始测井数据的协方差矩阵，所述初始测井数据为标准化后的测井曲线的测井数据；确定所述协方差矩阵的特征值及特征向量；对所述特征值进行从大到小排序，并根据公式确定大于预设阈值时的主成分测井曲线个数；其中，λ_i为第i个特征值，M为测井曲线个数，L为大于预设阈值时的主成分测井曲线个数；根据公式X_NL＝Y_NM*U_ML确定主成分测井数据；其中，X_NL为主成分测井数据，Y_NM为初始测井数据，U_ML为对初始测井数据的协方差矩阵求取的特征向量；

基于预设的聚类算法对所述主成分测井数据进行聚类，获得所述目标井的岩性识别结果。

2.如权利要求1所述的基于测井曲线的岩性识别方法，其特征在于，所述标准化所述测井曲线，包括：

对所述多条测井曲线进行标准差标准化处理。

3.如权利要求1所述的基于测井曲线的岩性识别方法，其特征在于，所述预设的聚类算法包括核K-均值聚类算法。

4.如权利要求3所述的基于测井曲线的岩性识别方法，其特征在于，对于所述目标井中的每口井，所述基于预设的聚类算法对所述主成分测井数据进行聚类，包括：

a、按照N/K的间隔取K个数据点作为初始聚类中心点，并根据所述初始聚类中心点对所述主成分测井数据进行初始分类；其中，N为主成分测井数据个数，K为初始分类数；

b、根据公式将主成分测井数据从输入空间映射到高维特征空间，求取核函数矩阵；其中，KM_i,j为核函数矩阵第i行第j列元素，x_i为第i个主成分测井数据，x_j为第j个主成分测井数据，κ(x_i,x_j)为第i个主成分测井数据与第j个主成分测井数据的核函数，φ(x_i)为第i个主成分测井数据到高维特征空间的映射，φ(x_j)为第j个主成分测井数据在高维特征空间的映射，＜φ(x_i),φ(x_j)＞为第i个主成分测井数据到高维空间的映射与第j个主成分测井数据到高维空间的映射的内积，σ为核参数，L为主成分测井曲线个数，x_il为第i个主成分测井数据的第l个分量，x_jl为第j个主成分测井数据的第l个分量；

c、统计每类的数据样本个数count_i；

d、根据公式p_i＝count_i/N确定每类所占的概率；其中p_i为第i类在N中所占的比例；

e、根据公式确定高维特征空间中每个主成分测井数据x_n与各类c_i的频率差异测度；其中，Dist_i为高维特征空间中每个主成分测井数据φ(x_n)与φ(c_i)类中心的频率差异测度，φ(x_n)为第n个主成分测井数据到高维特征空间的映射，c_i为第i类，φ(c_i)为c_i类中心到高维空间的映射，λ为惩罚因子；

f、确定每个主成分测井数据的最小Dist_i，并将每个主成分测井数据归类至与其最小Dist_i所对应的类别中；

g、更新聚类中心点，并重复c～f，直至满足预设的迭代终止条件；

h、输出每个主成分测井数据所属类别，并根据测井曲线特征，给各类别赋予岩性。