一种储层岩石渗透率预测方法、装置、设备及存储介质

本发明涉及地震岩石物理领域，特别涉及一种储层岩石渗透率预测方法、装置、设备及存储介质。

背景技术：

1、随着油气勘探开发程度的不断深入，面向的对象也越来越复杂。对储层参数的预测的难度越来越高。孔隙度和饱和度影响岩石的弹性性质，如岩石的模量、速度及密度。但是渗透率对这些参数并不敏感，现有技术只是认为渗透率会影响地震波在岩石中传播的速度频散及衰减，但是也缺乏定量的表征。因此渗透率的预测，一直是研究的重点与难点。

2、当前获得渗透率的办法主要有岩芯测量、测井和地震预测等。岩芯测量得到的渗透率精确度高，但是从地下钻取岩芯成本高且范围小；测井只能反映井点的储层渗透性，具有局限性；地震预测具有地震数据量丰富且成本较低较易得到的优势，因此研究较多，得到储层渗透性与地震频散关系很密切，但是还只是定性还没有做到定量预测。

3、因此，如何在低成本条件下，定量预测岩石渗透率是当前亟需解决的问题。

技术实现思路

1、有鉴于此，本发明的目的在于提供一种储层岩石渗透率预测方法、装置、设备及存储介质，解决了现有技术中无法实现低成本、高效率的定量预测岩石渗透率的问题。

2、为解决上述技术问题，本发明提供了一种储层岩石渗透率预测方法，包括：

3、利用岩石声电渗综合全自动测量系统得到岩芯在不同压力下的渗透率；

4、利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下频变体积模量的第一拟合函数，所述频变体积模量为不同频率的对数对应的体积模量；

5、根据所述不同压力下的渗透率和所述第一拟合函数确定渗透率与频变体积模量对数频率导数的关系函数；所述频变体积模量对数频率导数为所述频变体积模量对频率的对数求导得到的导数；

6、从地震数据中反演获取目标存储岩石的体积模量，并基于所述关系函数对所述目标储层岩石进行渗透率预测。

7、可选的，所述利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下频变体积模量的第一拟合函数，包括：

8、利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下的频变体积模量散点数据；

9、利用最小二乘法对所述不同压力下的频变体积模量散点数据进行拟合，得到不同压力下对应的所述第一拟合函数。

10、可选的，所述利用最小二乘法对所述不同压力下的频变体积模量散点数据进行拟合，得到不同压力下对应的所述第一拟合函数，包括：

11、利用三次函数拟合所述不同压力下的频变体积模量散点数据，得到不同压力下对应的所述第一拟合函数，公式为：

12、fj(logf)＝aj(logf)3+bj((logf))2+cj(logf)+dj

13、kji＝fj(logfi)

14、logfi表示第i个频率fi的对数；fj(logf)表示第j个压力pj对应的第一拟合函数；κji表示压力为pj时，第i个频率fi的体积模量；aj、bj、cj、dj表示fj(logf)的系数，当公式中lj最小时对应的aj、bj、cj、dj数值作为所述fj(logf)的系数；lj表示在第j个压力下，每个频率下所述第一拟合函数得到的值与所述频变体积模量散点数据的误差的平方和；x表示地震波的个数。

15、可选的，所述根据所述不同压力下的渗透率和所述第一拟合函数确定渗透率与频变体积模量对数频率导数的关系函数，包括：

16、将所述第一拟合函数对所述频率的对数求导数，得到不同压力下频变体积模量导数的第二拟合函数fj'，公式为

17、利用不同压力下的所述第二拟合函数fj'和所述不同压力下的渗透率κj确定所述渗透率与频变体积模量导数的关系函数。

18、可选的，所述利用所述第二拟合函数和所述不同压力下的渗透率确定所述渗透率与频变体积模量导数的关系函数，包括：

19、确定所述第二拟合函数与所述不同压力下的渗透率在各个压力、各个频率的各个交点；

20、利用二次函数分别拟合同频率下的各个交点，得到各频率对应的第三拟合函数，所述第三拟合函数为所述关系函数；所述关系函数以所述渗透率为横坐标，所述频变体积模量对数频率导数为纵坐标。

21、可选的，所述利用二次函数分别拟合同频率下的各个交点，得到各频率对应的第三拟合函数，包括：

22、利用二次函数拟合同压力同频率下的各个交点，公式为：

23、mi(κ)＝diκ2+eiκ+gi

24、mi(κ)表示第i个频率对应的第三拟合函数；dj、ej、gj为mi(κ)的系数，当公式中的zi最小时对应的dj、ej、gj数值作为mi(κ)的系数；zi表示在第i个频率时，每个渗透率下所述第三拟合函数得到的值与所述频变体积模量对数频率导数的误差的平方和；κfldji表示在第i个频率下，第j个频变体积模量对数频率导数；mi(κj)表示第i个频率下，第j个第三拟合函数值；n表示压力的个数。

25、本发明还提供了一种储层岩石渗透率预测装置，包括：

26、渗透率实测模块，用于利用岩石声电渗综合全自动测量系统得到岩芯在不同压力下的渗透率；

27、频变体积模量拟合函数获取模块，用于利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下频变体积模量的第一拟合函数，所述频变体积模量为不同频率的对数对应的体积模量；

28、关系函数确定模块，用于根据所述不同压力下的渗透率和所述第一拟合函数确定渗透率与频变体积模量对数频率导数的关系函数；所述频变体积模量对数频率导数为所述频变体积模量对频率的对数求导得到的导数；

29、渗透率预测模块，用于从地震数据中反演获取目标存储岩石的体积模量，并基于所述关系函数对所述目标储层岩石进行渗透率预测。

30、可选的，所述频变体积模量拟合函数获取模块，包括：

31、散点数据获取单元，用于利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下的频变体积模量散点数据；

32、第一拟合单元，用于利用最小二乘法对所述不同压力下的频变体积模量散点数据进行拟合，得到不同压力下对应的所述第一拟合函数。

33、本发明还提供了一种储层岩石渗透率预测设备，包括：

34、存储器，用于存储计算机程序；

35、处理器，用于执行所述计算机程序时实现上述的储层岩石渗透率预测方法。

36、本发明还提供了一种存储介质，所述存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器加载并执行时，实现上述的储层岩石渗透率预测方法。

37、可见，本发明通过利用岩芯声电渗综合全自动测量系统得到岩芯在不同压力下的渗透率；利用储层岩石跨频段地震弹性参数测试系统和应力应变法得到岩芯在不同压力下频变体积模量的第一拟合函数，频变体积模量为不同频率的对数对应的体积模量；根据不同压力下的渗透率和第一拟合函数确定渗透率与频变体积模量对数频率导数的关系函数；频变体积模量对数频率导数为频变体积模量对频率的对数求导得到的导数；从地震数据中反演获取目标存储岩石的体积模量，并基于关系函数对目标储层岩石进行渗透率预测。本发明从岩石物理实验出发，寻找地震波频散规律，得到频变体积模量对数频率导数与渗透率变化之间的定量关系，解决了岩芯测量成本高和效率低的问题，以及测井局限性问题。同时还考虑了不同压力和频率对渗透率的影响，为后续利用定量关系函数进行储层参数预测新理论推导提供了基础。

38、此外，本发明还提供了一种储层岩石渗透率预测装置、设备及存储介质，同样具有上述有益效果。