【技术实现步骤摘要】
一种核磁共振数据反演方法、装置、存储介质及设备
本申请涉及油气勘探中的核磁共振测井数据处理
，特别涉及一种核磁共振数据反演方法、装置、存储介质及设备。
技术介绍
核磁共振测井利用氢原子核与磁场相互作用发生的共振现象来实现井下油气储层信息的观测，是一种有效的探测岩石物理性质的方法。核磁共振测井测量地层的原始数据是成千上万个回波组成的回波串，通常将测量得到的回波串反演成核磁共振谱用于储层岩石物理参数的计算，从而反演核磁共振谱的精度直接关系着计算的岩石物理参数的准确性。所以，研究高精度的核磁共振数据反演方法非常重要。目前，已有的核磁共振数据反演方法包括正则化方法、迭代方法以及智能优化方法。这些方法在反演核磁共振数据时，目标函数都包括残差L2范数这一项，导致反演得到的解的稀疏性差，从而得到的核磁共振反演谱的分辨率较低。因此，亟需一种核磁共振数据反演方法，以提高核磁共振反演谱的分辨率。
技术实现思路
本申请实施例的目的是提供一种核磁共振数据反演方法、装置、存储介质及设备，以提高核磁共振反演谱的分辨率。...
【技术保护点】
1.一种核磁共振数据反演方法，其特征在于，包括：/n获取核磁共振回波数据；/n基于向量的Lp范数，建立目标函数；其中，所述向量满足预设的约束方程，所述约束方程通过所述核磁共振回波数据确定，所述Lp范数中的参数p的取值范围为大于等于0且小于等于1；/n根据所述目标函数的解，生成所述核磁共振回波数据的反演谱。/n
【技术特征摘要】
1.一种核磁共振数据反演方法，其特征在于，包括：
获取核磁共振回波数据；
基于向量的Lp范数，建立目标函数；其中，所述向量满足预设的约束方程，所述约束方程通过所述核磁共振回波数据确定，所述Lp范数中的参数p的取值范围为大于等于0且小于等于1；
根据所述目标函数的解，生成所述核磁共振回波数据的反演谱。
2.根据权利要求1所述的方法，其特征在于，所述目标函数为：
其中，为向量的Lp范数；向量s为所述目标函数的解，向量s中的所有元素都大于等于0；p为Lp范数的参数；
所述约束方程为：
Ks＝y
其中，向量y为不同时刻的回波幅度组成的向量，所述不同时刻的回波幅度通过所述核磁共振回波数据得到；K为核矩阵，所述核矩阵K中的元素为：
其中，t为回波幅度衰减时间，T2为横向弛豫时间。
3.根据权利要求2所述的方法，其特征在于，利用以下步骤得到所述目标函数的解：
设置参数ε、阈值tol、Lp范数中的参数p、最大迭代次数itmax，并将迭代次数it的初始值设置为1；所述参数ε和所述阈值tol的取值范围为大于0；
计算向量s(0)＝K\y，并确定向量s(0)中的元素个数n；
计算初始残差
进行第it次迭代计算，所述第it次迭代计算包括：计算第it次迭代后的解s(it)，以及第it次迭代后的残差r(it)；所述第it次迭代后的解s(it)，以及第it次迭代后的残差r(it)通过如下方法得到：
s(it)＝WKT/(KWKT)y
其中，对角矩阵W＝diag(1./w)，向量w＝[w1,w2,…wi,…wn]，在第1次迭代时，为所述向量s(0)中的第i个元素；在迭代次数it的值大于等于2时，为前一次迭代后的解s(it-1)中的第i个元素；
在|r(it)-r|＜tol；或，|r(it)-r|≥tol且it＝itmax的情况下，令s＝s(it)，得到所述向量s。
4.根据权利要求3所述的方法，其特征在于，还包括：
在|r(it)-r|≥tol且it＜itmax的情况下，令s＝s(it)、r＝r(it)；
将迭代次数it增加1后，进行第it次迭代计算。
5.根据权利要求3所述的方法，其特征在于，所述令s＝s(it)，得到向量s，包括：
令s...
【专利技术属性】
技术研发人员：郭江峰，谢然红，肖立志，徐陈昱，谷明宣，吴勃翰，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。