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

1.本发明涉及核磁共振测井数据处理技术领域，特别是涉及一种岩石物理参数的确定方法、装置及介质。


背景技术：

2.核磁共振测井是评价地下油气资源的一种重要地球物理测井方法，根据核磁共振测井数据反演得到的核磁共振谱可以进一步得到储层的孔隙度、饱和度和渗透率等岩石物理参数，因此，反演核磁共振谱的精度直接影响岩石物理参数的准确性。
3.核磁共振数据反演是一个病态问题，测量数据微小的噪声波动就会影响反演核磁共振谱的结果。目前反演核磁共振谱的方法采用正则化方法，在构建目标函数时采用残差l2范数，导致在对非常规储层低信噪比的核磁共振测井数据反演时其解的稀疏性较差，从而导致核磁共振谱的分辨率较差，影响后续核磁共振测井数据在非常规储层勘探中的应用效果。
4.因此，如何提高非常规储层的核磁共振测井数据反演精度是本领域技术人员亟需要解决的。


技术实现要素：

5.本发明的目的是提供一种岩石物理参数的确定方法、装置及介质，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
6.为解决上述技术问题，本发明提供一种岩石物理参数的确定方法，应用于油气勘探的核磁共振测井数据，包括：
7.获取多维核磁共振测井数据；
8.根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；
9.对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
10.在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
11.将最优反演解逆变换得到多维核磁共振谱；
12.根据多维核磁共振谱确定岩石物理参数。
13.优选地，根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，包括：
14.根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据；
15.将多维数据进行降维处理得到降维数据，其中降维数据包括降维后的多维核磁共振测井数据、待反演解和多维核磁共振测井数据的核矩阵；
16.将降维后的多维核磁共振测井数据、多维核磁共振测井数据的核矩阵和待反演解
进行范数处理得到残差l1范数和正则化l1范数；
17.将残差l1范数、正则化l1范数和正则化参数组合得到反演目标函数。
18.优选地，对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解，包括：
19.获取当前正则化参数、当前拉格朗日乘子、当前迭代次数、近端参数和预设迭代次数；
20.根据当前正则化参数得到当前惩罚参数；
21.根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数，其中拉格朗日乘子首次处理为0；
22.根据反演解函数确定当前反演解，其中当当前反演解存在负项时，将负项对应的数据置0得到新的当前反演解，当前反演解首次处理时为0；
23.判断当前迭代次数和当前反演解是否满足预设条件；
24.若是，则输出当前反演解；
25.若否，则返回至根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤。
26.优选地，预设条件为当前迭代次数达到预设迭代次数或当前反演解与上一个当前反演解之间的l1范数解小于阈值。
27.优选地，在多个正则化参数中选取最优正则化参数，包括：
28.设置正则化参数的参数范围；
29.根据对数均匀分布规则在参数范围内中选取各正则化参数；
30.将各正则化参数输入至反演目标函数得到对应的反演解；
31.根据各反演解计算各正则化参数对应的目标函数残差l1范数；
32.在各目标函数残差l1范数中选取最小l1范数并根据最小l1范数对应的正则化参数作为最优正则化参数。
33.优选地，根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据，包括：
34.根据第一类fredholm积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据。
35.优选地，反演解大于或等于0且正则化参数大于0。
36.为解决上述技术问题，本发明还提供一种岩石物理参数的确定装置，包括：
37.第一获取模块，用于获取多维核磁共振测井数据；
38.建立模块，用于根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；
39.迭代模块，用于对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
40.第二获取模块，用于在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
41.变换模块，用于将最优反演解逆变换得到多维核磁共振谱；
42.确定模块，用于根据多维核磁共振谱确定岩石物理参数。
43.为解决上述技术问题，本发明还提供一种岩石物理参数的确定装置，包括：
44.存储器，用于存储计算机程序；
45.处理器，用于执行计算机程序时实现如上述岩石物理参数的确定方法的步骤。
46.为解决上述技术问题，本发明还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如上述岩石物理参数的确定方法的步骤。
47.本发明提供的一种岩石物理参数的确定方法，应用于油气勘探的核磁共振测井数据，包括获取多维核磁共振测井数据；根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；将最优反演解逆变换得到多维核磁共振谱；根据多维核磁共振谱确定岩石物理参数。该方法通过反演目标函数包括残差l1范数和正则化l1范数以实现基于双稀疏约束，同时在多个正则化参数中选取最优正则化参数以根据建立的反演目标函数确定最优反演解，从而根据最优反演解逆变换得到多维核磁共振谱，进而根据多维核磁共振谱确定岩石物理参数。基于双稀疏约束的反演目标函数相较于现有的反演目标函数包括残差l2范数在对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其得到的核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
48.另外，本发明还提供了一种岩石物理参数的确定装置及介质，具有如上述岩石物理参数的确定方法相同的有益效果。
附图说明
49.为了更清楚地说明本发明实施例，下面将对实施例中所需要使用的附图做简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
50.图1为本发明实施例提供的一种岩石物理参数的确定方法的流程图；
51.图2为本发明实施例提供的一种岩石物理参数的确定装置的结构图；
52.图3为本发明实施例提供的另一种岩石物理参数的确定装置的结构图；
53.图4为数值模拟实验构造的含有三峰的t
1-t2谱的示意图；
54.图5为基于正演的数据信噪比40的t
1-t2数据的示意图；
55.图6为基于正演的数据信噪比10的t
1-t2数据的示意图；
56.图7为基于brd方法反演信噪比为40的数据得到的核磁共振t
1-t2谱的示意图；
57.图8为本发明提供的方法反演信噪比为40的数据得到的核磁共振t
1-t2谱的示意图；
58.图9为基于brd方法反演信噪比为10的数据得到的核磁共振t
1-t2谱的示意图；
59.图10为本发明提供的方法反演信噪比为10的数据得到的核磁共振t
1-t2谱的示意图。
具体实施方式
60.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下，所获得的所有其他实施例，都属于本发明保护范围。
61.本发明的核心是提供一种岩石物理参数的确定方法、装置及介质，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
62.为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。
63.需要说明的是，本发明提供的岩石物理参数的确定方法，适用于一维核磁共振测井数据和多维核磁共振测井数据，主要确定反演的核磁共振谱以确定岩石物理参数，反演的核磁共振谱的精度直接影响岩石物理参数的准确性。
64.图1为本发明实施例提供的一种岩石物理参数的确定方法的流程图，应用于油气勘探的核磁共振测井数据，如图1所示，该方法包括：
65.s11：获取多维核磁共振测井数据；
66.s12：根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；
67.s13：对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
68.s14：在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
69.s15：将最优反演解逆变换得到多维核磁共振谱；
70.s16：根据多维核磁共振谱确定岩石物理参数。
71.可以理解的是，以二维核磁共振测井数据为例，采集核磁共振测井数据，在采集的过程中数据为t
1-t2测量数据，弛豫过程是在射频脉冲信号解除后被激发的氢原子核逐步释放能量，其相位和能级都恢复到初始状态，磁化矢量回到平衡状态的过程。弛豫包括横向弛豫和纵向弛豫，为两个独立的过程。纵向弛豫是高能的质子释放能量回到低能级的过程，一般用t1值来表述纵向弛豫的快慢。横向弛豫能量过程是质子内部之间发生，为自旋弛豫，一般用t2值表述横向弛豫快慢。
72.以二维核磁共振t
1-t2谱的反演为例，通常是将二维数据转换至一维数据，再进行反演，将反演后的一维数据解逆变换到二维空间以实现反演。对于多维核磁共振的反演过程是将多维数据降维处理得到一维数据，进而反演，将反演后的一维数据解逆变换到多维空间以实现反演。
73.获取多维核磁共振测井数据，根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，可以理解的是，反演目标函数中的正则化参数是一个变量。现有的反演目标函数中的残差项为l2范数，其本发明实施例中反演目标函数的残差项为l1范数，正则化项也为li范数，两个l1范数以实现双稀疏约束，l1范数的模型假设参数是服从拉普拉斯分布，l2范数的模型假设参数是服从高斯分布。l1范数可以得到稀疏解，l1范数比l2范数能够产生更加稀疏的谱。
74.从而对不同正则化参数时反演目标函数进行迭代求解，以得到反演解。在反演目标函数中由多维核磁共振测井数据、采集参数和正则化参数得到，则后续根据选定的正则化参数输入至反演目标函数中，进行迭代求解。需要说明的是，本实施例中的迭代求解对于迭代方法不做具体限定，只要能得到最终的迭代解即可。
75.可以理解的是，正则化参数存在多个且不相同，根据步骤s13中对不同正则化参数时反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解，进而在不同正则化参数中选取最优正则化参数以确定对应的最优反演解。对于最优正则化参数的选取规则不做具体限定，可以基于l1范数选取，也可以基于l2范数选取，由于l1范数比l2范数其解的稀疏性较好，作为优选地实施例，根据l1范数选取最优正则化参数。
76.将确定的最优反演解作逆变换得到多维核磁共振谱，进而根据多维核磁共振谱确定岩石物理参数。其通过多维核磁共振谱确定岩石物理参数采用何种方法不做具体限定，本发明主要基于确定多维核磁共振谱，使其反演的多维核磁共振谱的精度得到提高才可确定岩石物理参数更为准确。
77.本发明实施例提供的一种岩石物理参数的确定方法，应用于油气勘探的核磁共振测井数据，包括获取多维核磁共振测井数据；根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；将最优反演解逆变换得到多维核磁共振谱；根据多维核磁共振谱确定岩石物理参数。该方法通过反演目标函数包括残差l1范数和正则化l1范数以实现基于双稀疏约束，同时在多个正则化参数中选取最优正则化参数以根据建立的反演目标函数确定最优反演解，从而根据最优反演解逆变换得到多维核磁共振谱，进而根据多维核磁共振谱确定岩石物理参数。基于双稀疏约束的反演目标函数相较于现有的反演目标函数包括残差l2范数在对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其得到的核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
78.在上述实施例的基础上，步骤s12中的根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，包括：
79.根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据；
80.将多维数据进行降维处理得到降维数据，其中降维数据包括降维后的多维核磁共振测井数据、待反演解和多维核磁共振测井数据的核矩阵；
81.将降维后的多维核磁共振测井数据、多维核磁共振测井数据的核矩阵和待反演解进行范数处理得到残差l1范数和正则化l1范数；
82.将残差l1范数、正则化l1范数和正则化参数组合得到反演目标函数。
83.具体地，以二维数据为例，将多维核磁共振测井数据转换为积分方程形式，其公式如下：
[0084][0085]
其中，b(tw,t)为核磁共振t
1-t2测量数据，tw为等待时间，f(t1,t2)为待求的t
1-t2谱，也就是待求的二维核磁共振谱，e(tw,t)为噪声数据。
[0086]
进一步将积分方程的形式转换为矩阵形式，通过以下公式表示：
[0087][0088]
其中，b对应上述核磁共振t
1-t2测量数据，f对应待求的t
1-t2谱，代表t
1-t2两个矩阵，e代表噪声数据。
[0089]
进一步，将二维矩阵形式的多维数据进行降维处理得到降维数据，其公式如下：
[0090]
b＝kf e
[0091]
其中，b＝vect(b)为核磁共振测井算子，f＝vect(f)为待反演解，e＝vect(e)为噪声算子，为两个矩阵的kronecker积，也就是核磁共振测井数据的核矩阵，vect()是将矩阵所有列按顺序叠加成一个向量的算子。因此，二维核磁共振数据反演可以看成是已知k和b，求解f的过程。
[0092]
将以上的降维数据，其中降维数据包括降维后的多维核磁共振测井数据、待反演解和多维核磁共振测井数据的核矩阵进行范数处理得到残差l1范数||kf-b||1和正则化l1范数α||f||1。将残差l1范数和正则化l1范数合并得到反演目标函数，其公式如下：
[0093][0094]
其中，α为正则化参数。
[0095]
本实施例提供的根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，通过反演目标函数包括残差l1范数和正则化l1范数以实现基于双稀疏约束。基于双稀疏约束的反演目标函数相较于现有的反演目标函数中的残差项为l2范数在对非常规储层低信噪比的核磁共振测井数据反演时得到的解的稀疏性较好，其得到的核磁共振谱的分辨率较高，提高非常规储层的核磁共振测井数据的反演精度以确定岩石物理参数。
[0096]
在上述实施例的基础上，步骤s13中的对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解，包括：
[0097]
获取当前正则化参数、当前拉格朗日乘子、当前迭代次数、近端参数和预设迭代次数；
[0098]
根据当前正则化参数得到当前惩罚参数；
[0099]
根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数，其中拉格朗日乘子首次处理为0；
[0100]
根据反演解函数确定当前反演解，其中当当前反演解存在负项时，将负项对应的数据置0得到新的当前反演解，当前反演解首次处理时为0；
[0101]
判断当前迭代次数和当前反演解是否满足预设条件；
[0102]
若是，则输出当前反演解；
[0103]
若否，则返回至根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤。
[0104]
获取当前正则化参数α，当前拉格朗日乘子ω，当前迭代次数iter，近端参数τ＝0.01，预设迭代次数iter
max
。对于预设迭代次数为经验值，可根据实际情况设定。
[0105]
根据当前正则化参数得到当前惩罚参数ρ，其公式如下：
[0106]
ρ＝2500/α
[0107]
根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格
朗日函数变形。具体地，将反演目标函数改写为以下公式：
[0108][0109]
根据公式构造增广拉格朗日函数，其公式如下：
[0110][0111]
进而得到反演解函数公式如下：
[0112]fiter
＝sign(z
iter
)
·
max(|z
iter
|-τ/ρ,0)
[0113]
其中，z
iter
＝f
iter-1-τ(k
t
(kf
iter-1-e
iter-b-ω
iter-1
/ρ))；
[0114]eiter
＝sign(kf
iter-1-b-ω
iter-1
/ρ)
·
max(|kf
iter-1-b-ω
iter-1
/ρ|-1/αρ,0)。
[0115]
拉格朗日乘子函数公式如下：
[0116]
ω
iter
＝ω
iter-1-ρ(kf
iter
e
iter-b)
[0117]
其中，拉格朗日乘子首次处理为0。
[0118]
根据上述的反演解函数确定当前反演解，若当前反演解存在负项时，则将负项对应的数据置0，得到新的当前反演解，对应首次处理时为0。
[0119]
判断当前迭代次数和当前反演解是否满足预设条件，若满足，则退出循环，输出当前反演解。若不满足，继续返回至根据反演目标函数、当前拉格朗日乘子、当前惩罚参数和近端参数进行增广拉格朗日函数变形得到反演解函数和拉格朗日乘子函数的步骤直到当前迭代次数和当前反演解满足预设条件输出当前反演解为止。
[0120]
本发明实施例提供的对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解以构造增广拉格朗日函数的方法进行迭代求解，实现数值优化。
[0121]
在上述实施例中，预设条件为当前迭代次数达到预设迭代次数或当前反演解与上一个当前反演解之间的l1范数解小于阈值。
[0122]
当前迭代次数达到预设迭代次数时iter》iter
max
，或者当前反演解与上一个当前反演解之间的关系确定的范数解小于阈值时，也就是相邻两次迭代结果的差值小于阈值，说明比较靠近，其结果较为收敛稳定，可以退出循环。
[0123]
相邻两次迭代结果的差值小于阈值通过如下公式：
[0124][0125]
需要说明的是，本发明实施例的预设条件也可以设定当前迭代次数达到预设迭代次数或当前反演解与上一个当前反演解之间的l2范数解小于阈值，不做具体限定，根据实际情况设定。
[0126]
本实施例提供的预设条件为当前迭代次数达到预设迭代次数或当前反演解与上一个当前反演解之间的l1范数解小于阈值，实现迭代求解。
[0127]
在上述实施例的基础上，步骤s14中的在多个正则化参数中选取最优正则化参数，包括：
[0128]
设置正则化参数的参数范围；
[0129]
根据对数均匀分布规则在参数范围内中选取各正则化参数；
[0130]
将各正则化参数输入至反演目标函数得到对应的反演解；
[0131]
根据各反演解计算各正则化参数对应的目标函数残差l1范数；
[0132]
在各目标函数残差l1范数中选取最小l1范数并根据最小l1范数对应的正则化参数作为最优正则化参数。
[0133]
具体地，正则化参数的参数范围设定为0.02-200之间，进而在参数范围内选取各正则化参数，由于参数范围较大，可以根据对数均匀分布规则选取，还可以根据参数的精度进而设置间隔等，不做具体限定。
[0134]
将各正则化参数输入至反演目标函数得到对应的反演解，根据各反演解计算各正则化参数对应的目标函数残差l1范数；在各目标函数残差l1范数中选取最小l1范数并根据最小l1范数对应的正则化参数作为最优正则化参数。
[0135]
需要说明的是，选取最小l1范数解对应的正则化参数，除了根据目标函数残差li范数求得，还可以根据其他方式选取，不做具体限定。
[0136]
本实施例提供的根据各反演解计算各正则化参数对应的目标函数残差l1范数；在各目标函数残差l1范数中选取最小l1范数并根据最小l1范数对应的正则化参数作为最优正则化参数。节省数据处理流程，简化选取过程，提高数据处理效率。
[0137]
在上述实施例的基础上，根据积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据，包括：
[0138]
根据第一类fredholm积分方程将多维核磁共振测井数据转换为矩阵形式的多维数据。
[0139]
具体地，将多维核磁共振测井数据转换为积分方程形式，对于何种积分方程形式不做具体限定，为了得到稳定的数值解，其通常采用正则化方法求解，也就是后续的输出处理过程中采取的正则化迭代求解。
[0140]
本实施例提供的根据第一类fredholm积分方程将核磁共振测井数据转换为矩阵形式的多维数据，使其后续求得的反演解更加稳定。
[0141]
在上述实施例的基础上，反演解大于或等于0，且正则化参数大于0，作为限定的条件，使其求得的反演解不存在负项，趋于稳定。
[0142]
上述详细描述了岩石物理参数的确定方法对应的各个实施例，在此基础上，本发明还公开与上述方法对应的岩石物理参数的确定装置，图2为本发明实施例提供的一种岩石物理参数的确定装置的结构图。如图2所示，岩石物理参数的确定装置包括：
[0143]
第一获取模块11，用于获取多维核磁共振测井数据；
[0144]
建立模块12，用于根据多维核磁共振测井数据、采集参数和正则化参数建立反演目标函数，其中反演目标函数包括残差l1范数和正则化l1范数；
[0145]
迭代模块13，用于对不同正则化参数的反演目标函数进行迭代求解确定不同正则化参数对应反演目标函数的反演解；
[0146]
第二获取模块14，用于在多个正则化参数中选取最优正则化参数并根据最优正则化参数获取对应的最优反演解；
[0147]
变换模块15，用于将最优反演解逆变换得到多维核磁共振谱；
[0148]
确定模块16，用于根据多维核磁共振谱确定岩石物理参数。
[0149]
由于装置部分的实施例与上述的实施例相互对应，因此装置部分的实施例请参照
上述方法部分的实施例描述，在此不再赘述。
[0150]
对于本发明提供的一种岩石物理参数的确定装置的结构图的介绍请参照上述方法实施例，本发明在此不再赘述，其具有上述岩石物理参数的确定方法相同的有益效果。
[0151]
图3为本发明实施例提供的另一种岩石物理参数的确定装置的结构图，如图3所示，该装置包括：
[0152]
存储器21，用于存储计算机程序；
[0153]
处理器22，用于执行计算机程序时实现岩石物理参数的确定方法的步骤。
[0154]
本实施例提供的岩石物理参数的确定装置可以包括但不限于智能手机、平板电脑、笔记本电脑或者台式电脑等。
[0155]
其中，处理器22可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器22可以采用数字信号处理器(digital signal processor，dsp)、现场可编程门阵列(field－programmable gate array，fpga)、可编程逻辑阵列(programmable logic array，pla)中的至少一种硬件形式来实现。处理器22也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称中央处理器(central processing unit，cpu)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器22可以集成有图像处理器(graphics processing unit，gpu)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器22还可以包括人工智能(artificial intelligence，ai)处理器，该ai处理器用于处理有关机器学习的计算操作。
[0156]
存储器21可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器21还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。本实施例中，存储器21至少用于存储以下计算机程序211，其中，该计算机程序被处理器22加载并执行之后，能够实现前述任一实施例公开的岩石物理参数的确定方法的相关步骤。另外，存储器21所存储的资源还可以包括操作系统212和数据213等，存储方式可以是短暂存储或者永久存储。其中，操作系统212可以包括windows、unix、linux等。数据213可以包括但不限于岩石物理参数的确定方法所涉及到的数据等等。
[0157]
在一些实施例中，岩石物理参数的确定装置还可包括有显示屏23、输入输出接口24、通信接口25、电源26以及通信总线27。
[0158]
领域技术人员可以理解，图3中示出的结构并不构成对岩石物理参数的确定装置的限定，可以包括比图示更多或更少的组件。
[0159]
处理器22通过调用存储于存储器21中的指令以实现上述任一实施例所提供的岩石物理参数的确定方法。
[0160]
对于本发明提供的一种岩石物理参数的确定装置的结构图的介绍请参照上述方法实施例，本发明在此不再赘述，其具有上述岩石物理参数的确定方法相同的有益效果。
[0161]
进一步的，本发明还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器22执行时实现如上述岩石物理参数的确定方法的步骤。
[0162]
可以理解的是，如果上述实施例中的方法以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明
的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可以存储程序代码的介质。
[0163]
对于本发明提供的一种计算机可读存储介质的介绍请参照上述方法实施例，本发明在此不再赘述，其具有上述岩石物理参数的确定方法相同的有益效果。
[0164]
作为一种优选地实施例，以二维核磁共振t
1-t2谱反演为例，进行数值模拟实验验证本发明提供的岩石物理参数的确定方法。
[0165]
数值模拟实验首先构造三峰的t
1-t2谱，三个峰对应的(t1，t2)值分别为(8，5)ms，(50，30)ms和(160，150)ms。设置40组等待时间[0.2,0.5,0.8,1.0,1.5,2.0,3.0,4.0,5.0,8.0,10,12,15,18,20,25,28,30,40,50,60,70,80,90,100,120,150,180,200,220,250,280,300,400,500,600,700,800,900,1000]ms。回波间隔为0.2ms，回波个数为4096。根据设置的采集参数，正演计算得到t
1-t2数据，并加入不同信噪比(40和10)的噪声，得到对应的含噪声的t
1-t2数据。再分别利用核磁共振数据反演方法(butler-reeds-dawson，brd)和提出的方法反演含噪声的t
1-t2数据，得到对应的t
1-t2谱。
[0166]
图4为数值模拟实验构造的含有三峰的t
1-t2谱的示意图。图5为基于正演的数据信噪比40的t
1-t2数据的示意图，图6为基于正演的数据信噪比10的t
1-t2数据的示意图。图5、图6为基于正演的含噪声的t
1-t2数据。将图5、图6中的回波数据个数压缩到400，再分别利用brd方法和提出的方法反演压缩后的回波数据，得到不用信噪比条件下核磁共振t
1-t2谱，图7为基于brd方法反演信噪比为40的数据得到的核磁共振t
1-t2谱的示意图；图8为本发明提供的方法反演信噪比为40的数据得到的核磁共振t
1-t2谱的示意图；图9为基于brd方法反演信噪比为10的数据得到的核磁共振t
1-t2谱的示意图；图10为本发明提供的方法反演信噪比为10的数据得到的核磁共振t
1-t2谱的示意图。
[0167]
当信噪比为40时，brd方法和提出的方法反演的t
1-t2谱与模型均方根误差rmse分别为3.38
×
10-5
和3.23
×
10-5
，反演t
1-t2谱与模型一致性较好。当信噪比为10时，brd方法反演的t
1-t2谱中不同流体信号峰发生了重叠，与模型均方根误差rmse为5.92
×
10-5
，而本发明提供的方法反演的t
1-t2谱与模型均方根误差rmse为4.05
×
10-5
，提出的方法反演的t
1-t2谱精度高于brd方法反演的t
1-t2谱。本发明提供的方法反演的t
1-t2谱受数据噪声影响较小。
[0168]
以上对本发明所提供的一种岩石物理参数的确定方法、装置及介质进行了详细介绍。说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
[0169]
还需要说明的是，在本说明书中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”“包含”或者其任何其他变体意
在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。