本发明专利技术公开了一种岩石物理参数的确定方法、装置及介质，适用于核磁共振测井数据处理技术领域。根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其反演目标函数包括残差L1范数和正则化L1范数以实现基于双稀疏约束，同时在多个正则化参数中选取最优正则化参数以确定最优反演解，根据最优反演解逆变换得到多维核磁共振谱的反演，根据反演结果确定岩石物理参数。基于双稀疏约束的反演目标函数相较于现有的反演目标函数的残差项L2范数对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。井数据的反演精度以确定岩石物理参数。井数据的反演精度以确定岩石物理参数。
【技术实现步骤摘要】
一种岩石物理参数的确定方法、装置及介质
[0001]本专利技术涉及核磁共振测井数据处理
，特别是涉及一种岩石物理参数的确定方法、装置及介质。
技术介绍
[0002]核磁共振测井是评价地下油气资源的一种重要地球物理测井方法，根据核磁共振测井数据反演得到的核磁共振谱可以进一步得到储层的孔隙度、饱和度和渗透率等岩石物理参数，因此，反演核磁共振谱的精度直接影响岩石物理参数的准确性。
[0003]核磁共振数据反演是一个病态问题，测量数据微小的噪声波动就会影响反演核磁共振谱的结果。目前反演核磁共振谱的方法采用正则化方法，在构建目标函数时采用残差L2范数，导致在对非常规储层低信噪比的核磁共振测井数据反演时其解的稀疏性较差，从而导致核磁共振谱的分辨率较差，影响后续核磁共振测井数据在非常规储层勘探中的应用效果。
[0004]因此，如何提高非常规储层的核磁共振测井数据反演精度是本领域技术人员亟需要解决的。
技术实现思路
[0005]本专利技术的目的是提供一种岩石物理参数的确定方法、装置及介质，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
[0006]为解决上述技术问题，本专利技术提供一种岩石物理参数的确定方法，应用于油气勘探的核磁共振测井数据，包括：
[0007]获取多维核磁共振测井数据；
[0008]根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差L1范数和正则化L1范数；
[0009]对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
[0010]在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
[0011]将最优反演解逆变换得到多维核磁共振谱；
[0012]根据多维核磁共振谱确定岩石物理参数。
[0013]优选地，根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，包括：
[0014]根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据；
[0015]将多维数据进行降维处理得到降维数据，其中降维数据包括降维后的多维核磁共振测井数据、待反演解和多维核磁共振测井数据的核矩阵；
[0016]将降维后的多维核磁共振测井数据、多维核磁共振测井数据的核矩阵和待反演解
进行范数处理得到残差L1范数和正则化L1范数；
[0017]将残差L1范数、正则化L1范数和正则化参数组合得到反演目标函数。
[0018]优选地，对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解，包括：
[0019]获取当前正则化参数、当前拉格朗日乘子、当前迭代次数、近端参数和预设迭代次数；
[0020]根据当前正则化参数得到当前惩罚参数；
[0021]根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数，其中拉格朗日乘子首次处理为0；
[0022]根据反演解函数确定当前反演解，其中当当前反演解存在负项时，将负项对应的数据置0得到新的当前反演解，当前反演解首次处理时为0；
[0023]判断当前迭代次数和当前反演解是否满足预设条件；
[0024]若是，则输出当前反演解；
[0025]若否，则返回至根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤。
[0026]优选地，预设条件为当前迭代次数达到预设迭代次数或当前反演解与上一个当前反演解之间的L1范数解小于阈值。
[0027]优选地，在多个正则化参数中选取最优正则化参数，包括：
[0028]设置正则化参数的参数范围；
[0029]根据对数均匀分布规则在参数范围内中选取各正则化参数；
[0030]将各正则化参数输入至反演目标函数得到对应的反演解；
[0031]根据各反演解计算各正则化参数对应的目标函数残差L1范数；
[0032]在各目标函数残差L1范数中选取最小L1范数并根据最小L1范数对应的正则化参数作为最优正则化参数。
[0033]优选地，根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据，包括：
[0034]根据第一类Fredholm积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据。
[0035]优选地，反演解大于或等于0且正则化参数大于0。
[0036]为解决上述技术问题，本专利技术还提供一种岩石物理参数的确定装置，包括：
[0037]第一获取模块，用于获取多维核磁共振测井数据；
[0038]建立模块，用于根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差L1范数和正则化L1范数；
[0039]迭代模块，用于对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
[0040]第二获取模块，用于在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
[0041]变换模块，用于将最优反演解逆变换得到多维核磁共振谱；
[0042]确定模块，用于根据多维核磁共振谱确定岩石物理参数。
[0043]为解决上述技术问题，本专利技术还提供一种岩石物理参数的确定装置，包括：
[0044]存储器，用于存储计算机程序；
[0045]处理器，用于执行计算机程序时实现如上述岩石物理参数的确定方法的步骤。
[0046]为解决上述技术问题，本专利技术还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如上述岩石物理参数的确定方法的步骤。
[0047]本专利技术提供的一种岩石物理参数的确定方法，应用于油气勘探的核磁共振测井数据，包括获取多维核磁共振测井数据；根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差L1范数和正则化L1范数；对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；将最优反演解逆变换得到多维核磁共振谱；根据多维核磁共振谱确定岩石物理参数。该方法通过反演目标函数包括残差L1范数和正则化L1范数以实现基于双稀疏约束，同时在多个正则化参数中选取最优正则化参数以根据建立的反演目标函数确定最优反演解，从而根据最优反演解逆变换得到多维核磁共振谱，进而根据多维核磁共振谱确定岩石物理参数。基于双稀疏约束的反演目标函数相较于现有的反演目标函数包括残差L2范数在对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其得到的核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
[0048]另外本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种岩石物理参数的确定方法，其特征在于，应用于油气勘探的核磁共振测井数据，包括：获取多维核磁共振测井数据；根据所述多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中所述反演目标函数包括残差L1范数和正则化L1范数；对不同所述正则化参数的反演目标函数进行迭代求解确定不同所述正则化参数对应反演目标函数的反演解；在多个所述正则化参数中选取最优正则化参数并根据所述最优正则化参数获取对应的最优反演解；将所述最优反演解逆变换得到多维核磁共振谱；根据所述多维核磁共振谱确定所述岩石物理参数。2.根据权利要求1所述的岩石物理参数的确定方法，其特征在于，所述根据所述多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，包括：根据积分方程将所述多维核磁共振测井数据转换为矩阵形式的多维数据；将所述多维数据进行降维处理得到降维数据，其中所述降维数据包括降维后的多维核磁共振测井数据、待反演解和所述多维核磁共振测井数据的核矩阵；将所述降维后的多维核磁共振测井数据、所述多维核磁共振测井数据的核矩阵和所述待反演解进行范数处理得到所述残差L1范数和所述正则化L1范数；将所述残差L1范数、所述正则化L1范数和所述正则化参数组合得到所述反演目标函数。3.根据权利要求2所述的岩石物理参数的确定方法，其特征在于，所述对不同所述正则化参数的反演目标函数进行迭代求解确定不同所述正则化参数对应反演目标函数的反演解，包括：获取当前正则化参数、当前拉格朗日乘子、当前迭代次数、近端参数和预设迭代次数；根据所述当前正则化参数得到当前惩罚参数；根据所述反演目标函数、所述当前拉格朗日乘子、所述当前惩罚参数和所述近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数，其中所述拉格朗日乘子首次处理为0；根据所述反演解函数确定当前反演解，其中当所述当前反演解存在负项时，将所述负项对应的数据置0得到新的所述当前反演解，所述当前反演解首次处理时为0；判断所述当前迭代次数和所述当前反演解是否满足预设条件；若是，则输出所述当前反演解；若否，则返回至所述根据所述反演目标函数、所述当前拉格朗日乘子、所述当前惩罚参数和所述近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤。4...
【专利技术属性】
技术研发人员：郭江峰，谢然红，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。