(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202110914265.2 (22)申请日 2021.08.10 (71)申请人 林磊 地址 241000 安徽省芜湖市无为县赫店镇 留桥行政村 鲁林自然村24 号 申请人 蒋恕　钟志 (72)发明人 林磊　蒋恕　钟志　 (74)专利代理 机构 北京神州信德知识产权代理 事务所(普通 合伙) 11814 代理人 刘真 (51)Int.Cl. E21B 47/00(2012.01) G06F 30/27(2020.01) G06Q 50/02(2012.01) (54)发明名称 一种基于二维核磁共振测井数据识别页岩 油储层流体成分及相对体积的方法 (57)摘要 本发明公开了一种基于二维核磁共振测井 数据识别页岩油储层流体成分及相对体积的方 法， 包括以下步骤： S1.构建T1 ‑T2二维核磁学习 样本库； S2.将学习样本库中的每一个T1 ‑T2二维 数据拼接组成一个大型二维矩阵； S3.初始化源 数， 得到初始化源矩阵和初始化系数矩阵； S4.对 初始化的系数矩阵进行主成分分析， 确定最佳源 数； S5.将非负矩阵分解的源数设置为最佳源数， 得到源矩阵和系数矩阵； S6.确定每类流体质心 对应的T1和T2坐标； S7.对数据库中每一个T1 ‑T2 图数据中的像素点按照距离以每类流体质心欧 式距离的远近对流体类型进行划分； S8.根据流 体类型划分结果确定每类流体的总信号强度， 根 据信号强度比值确定 每类流体的相对体积。 权利要求书2页 说明书7页 附图4页 CN 114000864 A 2022.02.01 CN 114000864 A 1.一种基于二维核磁共振测井数据识别页岩油储层流体成分及相对体积的方法， 其特 征在于该 方法包括以下步骤： S1.对T1‑T2二维核磁共振测井数据按采样点提取数据， 构建T1 ‑T2二维核磁学习样本 库； S2.将学习样本库中的每一个T1 ‑T2二维数据向量化为一列向量， 再拼接组成一个大型 二维矩阵； S3.初始化源数， 对由S1步骤中形成的样本库中的数据所建立的大型二维矩阵进行非 负矩阵分解， 得到初始 化源矩阵和初始 化系数矩阵； 其中， 源 数为页岩储层中流体类型的个 数； S4.对初始化的系数矩阵进行主成分 分析， 确定最佳源数； S5.将非负矩阵分解的源数设置为最佳源数， 再次对样本库建立的大型矩阵进行非负 矩阵分解， 得到源矩阵和系数矩阵； S6.源矩阵中每一列代表一类流体类型对应的T1 ‑T2图， 从而确定每类流体质心对应的 T1和T2坐标； S7.对数据库中每一个T1 ‑T2图数据中的像素点按照距离以每类流体质心欧式距离的 远近对流体 类型进行划分； S8.根据流体类型划分结果确定每类流体的总信号强度， 根据信号强度比值确定每类 流体的相对体积。 2.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S1中， 对应一个的二维矩阵， 按采样点从上至下逐个提取T1 ‑T2二维 核磁共振测井数据中的T1 ‑T2图像， 构成学习样本库， 以提高学习样本库的构建速度和效 率。 3.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S2中， 每一幅T1 ‑T2图都是一个64 ×64的二维矩阵， 将每一幅T1 ‑T2图 的二维矩阵都向量化为一个4096维的列向量， 再将所有的列向量在列的维度上拼接为大型 的二维矩阵。 4.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于在S3中， S2 中构建的大型矩阵V， 可以近似为源矩阵W和系数矩阵H进 行矩阵乘法得到的结果， 如公式 V≈WH所示， 并利用非负矩阵分解的无监督学习算法求解源矩阵W和系数矩阵H， 通过公 式 C＝‖V‑WH‖2定义二范数损失函数， C为损失值， 在W＞0,H＞0的约束下通过乘法更新准则 的梯度下降原理求解损失函数的最小值； 其中源矩阵W的每一列是一种流体类型的T1 ‑T2图 的列向量形式， 系数矩阵H的每一列是样本库中合成T1 ‑T2图的的一组系数， 给定不同的分 解源数r既流体 类别数， 会得到不同的分解结果； 在S3中， 初始化给定的源数r＝8。 5.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S4中， 对S 3中获得的初始 化系数矩阵进 行主成分分析， 得到每个主成 分所携带 的解释方差的估计， 给定一个累计解释方差阈值， 超过该累计方差阈值的最小主 成分数就是最佳源数。权　利　要　求　书 1/2 页 2 CN 114000864 A 26.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S 5中， 将待学习的大型二 维矩阵V的分解源数设置为S4得到的最佳分 解源数4， 对矩阵V进 行二次非负矩阵分解， 得到r＝4对应的源矩阵和系数矩阵， 将源矩阵每 一列提取 出来恢复成64 ×64的二维矩阵。 7.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S6中， 计算S 5得到的四种源流体类型T1 ‑T2图中信号 强度大于给定信 号强度下限的像素点所对应的T1和T2值， 计算每种源流体类型的质心， 质心的计算公式如 下： 式中， T1C为某一类流体的质心T1轴的坐标； T2C为某一类流体 的质心T2轴的坐标； T1i， T2i分别为T1 ‑T2图中信号强度大于给定信号强度下限的第i个 像素点所对应的T1和T2值。 8.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S7中， 计算数据库每一幅T1 ‑T2图中的信号强度大于信号强度下限 0.05的像素点距离S 6中得到的四种流体类型的质心的欧式距离， 将其分给距离质心欧式距 离最小的这 一类流体， 完成流体 类型的划分。 9.如权利要求1所述的基于二维核磁共振测井数据识别页岩油储层流体成分及相对体 积的方法， 其特征在于S8中， 根据S7中流体类型划分 的结果， 计算数据库中每 幅T1‑T2图中 各类流体的信号强度总和， 然后计算不同类型流体的总信号强度比值， 该比值就是每类流 体的相对体积， 完成不同流体 类型的相对体积计算。权　利　要　求　书 2/2 页 3 CN 114000864 A 3