买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：西南石油大学

摘要：本发明公开了一种有效拆解核磁共振弛豫时间谱的方法，步骤S1、根据岩心核磁共振T2谱中孔隙和溶洞的分布特征构造孔隙和溶洞两部分的概率密度函数；S2、令核磁谱面积和为1，在孔洞分布的区间内随机生成n个数据点，并通过筛选峰值将所有数据点划分为多个数据集；S3、计算似然函数；S4、计算每个数据点属于孔隙或溶洞的先验概率；S5、利用先验概率更新迭代法中的各参数；S6、基于更新后的参数，计算新的对数似然性，并更新迭代法中的期望Ep和Ev，期望Ep和Ev；重复步骤S2至步骤S5，直至满足收敛条件，迭代终止。本发明将核磁共振测试与迭代法相结合拆解核磁共振谱，提升了谱形分布的准确度。

主权项：1.一种有效拆解核磁共振弛豫时间谱的方法，其特征在于，包括以下步骤：S1、根据岩心核磁共振T2谱中孔隙和溶洞的分布特征构造孔隙和溶洞两部分的概率密度函数；S2、数据初始化；首先计算初始对数似然选择初始估计，并对每组峰幅度归一化处理，令核磁谱面积和为1，在孔洞分布的区间内随机生成n个数据点，并通过筛选峰值将所有数据点划分为多个数据集；数据集的个数等于核磁共振谱峰值的数量，每个数据集的分布范围用孔隙度分量出现的频次来表征；S3、计算似然函数；将核磁共振弛豫时间看作有限个独立的数据集组成，引入孔隙部分的期望和方差，溶洞部分的期望和方差，求解分布函数的最大似然估计；使用迭代法时，假设数据点都是独立的，似然函数由概率密度函数给出，似然函数公式如下； 式中：n为数据点的数量；wp为数据点属于孔隙部分的概率；wv为数据点属于溶洞部分的概率；t2i为第i个数据点，i＝1，2…n；为孔隙部分的密度函数，为溶洞部分的密度函数；0代表迭代次数；S4、确定期望Ep和Ev，期望Ep和Ev分别是根据当前参数计算的每个数据点来自孔隙或溶洞的概率，每个数据点属于孔隙或溶洞的先验概率计算公式如下： 式中：Tip为第i个数据点属于孔隙的先验概率；Tiv为第i个数据点属于溶洞的先验概率；np和nv分别为孔隙和溶洞的期望矩阵；m为迭代次数；l为数据点序号；K表示高斯模型的个数；为第l个数据点的权重矩阵；S5、利用先验概率更新迭代法中的各参数；S6、基于更新后的参数，依据下述公式计算新的对数似然性，并更新迭代法中孔隙和溶洞的概率wp和wv，期望Ep和Ev，方差∑p和∑v；重复步骤S2至步骤S5，直至满足收敛条件迭代终止；最终得到迭代后的孔隙和溶洞的期望Ep和Ev、方差∑p和∑v； 式中，m和m+1均为迭代次数。

全文数据：

权利要求：

百度查询： 西南石油大学 一种有效拆解核磁共振弛豫时间谱的方法