[发明专利]基于多GPU的数字岩心相渗曲线的计算方法、装置及计算机设备

权利要求书

1.一种基于多GPU的数字岩心相渗曲线的计算方法，其特征在于，包括：

对在网格空间中根据微岩心所构建的数字岩心，获取并存储各孔隙所在孔隙点的网格数据；

对所述数字岩心进行油驱替饱和水模拟，在所述油驱替饱和水模拟中根据所述网格数据进行数值计算，获得各个时刻的油相参数和水相参数，直至所述油驱替饱和水模拟达到稳定状态；

以所述稳定状态时的油相参数和水相参数作为初始参数，对所述数字岩心进行水驱替油模拟，在所述模拟中，根据所述初始参数以及所述网格数据计算得到为所述微岩心所预测的相渗曲线，所述相渗曲线是根据在所述水驱替油模拟过程中各个时刻的含水饱和度、油相渗透率和水相渗透率所获得的。

2.根据权利要求1所述的方法，其特征在于，所述对在网格空间中根据微岩心所构建的数字岩心，获取并存储各孔隙所在孔隙点的网格数据之前，还包括：

获取对所述微岩心进行扫描所得到的岩心三维灰度图像；

对所述岩心三维灰度图像进行图像识别，获得所述微岩心的孔隙骨架三维数据；

在网格空间根据所述孔隙骨架三维数据进行所述微岩心的三维构建，得到所述数字岩心。

3.根据权利要求1所述的方法，其特征在于，所述对在网格空间中根据微岩心所构建的数字岩心，获取并存储各孔隙所在孔隙点的网格数据，包括：

根据所述数字岩心中各模型单元所对应的结构类型标示值，获得孔隙所在的模型单元；

根据所述孔隙所在的模型单元，获得孔隙所在的孔隙点，所述孔隙点是所述孔隙所在的模型单元在所述网格空间中所占据的网格点；

在为所述孔隙点所分配的孔隙数据组中写入网格数据，所述网格数据包括所述孔隙点的坐标、初始化配置的分布函数和密度函数；

将所述孔隙数组存储于所对应分配的GPU中，以在所述油驱替饱和水模拟和所述水驱替油模拟过程中通过调用所述GPU中的孔隙数组进行对应参数的计算。

4.根据权利要求3所述的方法，其特征在于，所述孔隙数据中还包括直接邻居孔隙索引信息，所述将所述孔隙数组存储于所对应分配的GPU中之前，还包括：

根据在所述孔隙所在的模型单元中与所述孔隙点直接相邻的若干直接邻居孔隙点，构建所述孔隙点所对应孔隙数组与所述直接邻居孔隙点所对应孔隙数据的索引；

根据所构建的索引生成直接邻居孔隙索引信息，并将所生成的直接邻居孔隙索引信息写入所述孔隙所对应孔隙数组中。

5.根据权利要求1所述的方法，其特征在于，所述对所述数字岩心进行油驱替饱和水模拟，在所述油驱替饱和水模拟中根据所述网格数据进行数值计算得到各个时刻的油相参数和水相参数，直至所述油驱替饱和水模拟达到稳定状态之前，还包括：

对所述数字岩心进行压汞模拟，得到孔喉半径所占据的网格数量；

判断所述所占据的网格数量是否低于设定的阈值；

若为是，则对所述网格空间中的网格进行加密；若为否，则执行对所述数字岩心进行油驱替饱和水模拟的步骤。

6.根据权利要求1所述的方法，其特征在于，所述数字岩心包括入口和出口，所述对所述数字岩心进行油驱替饱和水模拟，在所述油驱替饱和水模拟中根据所述网格数据进行数值计算得到各个时刻的油相参数和水相参数，直至所述油驱替饱和水模拟达到稳定状态，包括：

为所述数字岩心配置所述入口的油相压力P1和所述出口的水相压力P2，其中P1＞P2；

在所述入口处的油相压力P1与所述出口处的水相低压P2所形成压力差的驱动下，进行所述油驱替饱和水模拟；

在所述油驱替饱和水模拟过程中，根据所述网格数据进行数值计算，得到各个时刻各孔隙点的油相参数和水相参数；

若计算所得到的油相参数和水相参数在设定时间内保持不变，则确定所述油驱替饱和水模拟达到稳定状态。