一种岩石物理参数的确定方法、装置及介质.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN114994776A(43)申请公布日2022.09.02(21)申请号202210581419.5(22)申请日2022.05.26(71)申请人中国石油大学（北京）地址102299北京市昌平区府学路18号(72)发明人郭江峰谢然红(74)专利代理机构北京集佳知识产权代理有限公司11227专利代理师刘珂(51)Int.Cl.G01V3/32(2006.01)G01V3/38(2006.01)权利要求书2页说明书11页附图5页(54)发明名称一种岩石物理参数的确定方法、装置及介质(57)摘要本发明公开了一种岩石物理参数的确定方法、装置及介质，适用于核磁共振测井数据处理技术领域。根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其反演目标函数包括残差L1范数和正则化L1范数以实现基于双稀疏约束，同时在多个正则化参数中选取最优正则化参数以确定最优反演解，根据最优反演解逆变换得到多维核磁共振谱的反演，根据反演结果确定岩石物理参数。基于双稀疏约束的反演目标函数相较于现有的反演目标函数的残差项L2范数对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。CN114994776ACN114994776A权利要求书1/2页1.一种岩石物理参数的确定方法，其特征在于，应用于油气勘探的核磁共振测井数据，包括：获取多维核磁共振测井数据；根据所述多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中所述反演目标函数包括残差L1范数和正则化L1范数；对不同所述正则化参数的反演目标函数进行迭代求解确定不同所述正则化参数对应反演目标函数的反演解；在多个所述正则化参数中选取最优正则化参数并根据所述最优正则化参数获取对应的最优反演解；将所述最优反演解逆变换得到多维核磁共振谱；根据所述多维核磁共振谱确定所述岩石物理参数。2.根据权利要求1所述的岩石物理参数的确定方法，其特征在于，所述根据所述多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，包括：根据积分方程将所述多维核磁共振测井数据转换为矩阵形式的多维数据；将所述多维数据进行降维处理得到降维数据，其中所述降维数据包括降维后的多维核磁共振测井数据、待反演解和所述多维核磁共振测井数据的核矩阵；将所述降维后的多维核磁共振测井数据、所述多维核磁共振测井数据的核矩阵和所述待反演解进行范数处理得到所述残差L1范数和所述正则化L1范数；将所述残差L1范数、所述正则化L1范数和所述正则化参数组合得到所述反演目标函数。3.根据权利要求2所述的岩石物理参数的确定方法，其特征在于，所述对不同所述正则化参数的反演目标函数进行迭代求解确定不同所述正则化参数对应反演目标函数的反演解，包括：获取当前正则化参数、当前拉格朗日乘子、当前迭代次数、近端参数和预设迭代次数；根据所述当前正则化参数得到当前惩罚参数；根据所述反演目标函数、所述当前拉格朗日乘子、所述当前惩罚参数和所述近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数，其中所述拉格朗日乘子首次处理为0；根据所述反演解函数确定当前反演解，其中当所述当前反演解存在负项时，将所述负项对应的数据置0得到新的所述当前反演解，所述当前反演解首次处理时为0；判断所述当前迭代次数和所述当前反演解是否满足预设条件；若是，则输出所述当前反演解；若否，则返回至所述根据所述反演目标函数、所述当前拉格朗日乘子、所述当前惩罚参数和所述近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤。4.根据权利要求3所述的岩石物理参数的确定方法，其特征在于，所述预设条件为所述当前迭代次数达到所述预设迭代次数或所述当前反演解与上一个所述当前反演解之间的L1范数解小于阈值。5.根据权利要求1所述的岩石物理参数的确定方法，其特征在于，在多个所述正则化参2CN114994776A权利要求书2/2页数中选取所述最优正则化参数，包括：设置所述正则化参数的参数范围；根据对数均匀分布规则在所述参数范围内中选取各正则化参数；将各所述正则化参数输入至所述反演目标函数得到对应的所述反演解；根据各所述反演解计算各所述正则化参数对应的目标函数残差L1范数；在各所述目标函数残差L1范数中选取最小L1范数并根据所述最小L1范数对应的所述正则化参数作为所述最优正则化参数。6.根据权利要求2所述的岩石物理参数的确定方法，其特征在于，所述根据积分方程将所述多维核磁共振测井数据转换为矩阵形式的多维数据，包括：根据第一类Fredholm积分方程将所述多维核磁共振测井数据转换为矩阵形式的所述多维数据。7.根据权利要求1至6任意一项所述的岩石物