[发明专利]一种基于侵蚀-膨胀算法的数字岩心两相流模拟结果定量表征方法及其应用

摘要

本发明涉及一种基于侵蚀‑膨胀算法的数字岩心两相流模拟结果定量表征方法及其应用，本发明利用侵蚀‑膨胀算法提取数字岩心的孔隙网络模型，实现数字岩心中两相流动直接模拟结果的定量化表征。本发明准确描述了每一个孔隙和喉道的两相流体的饱和度，描述了不同直径的孔隙或喉道中两相流体的饱和度，避免了孔隙网络模型流动模拟误差大的缺点，同时又借鉴了孔隙网络模型的孔隙表征方法，将直接流动模拟的结果进行定量化描述；不同直径孔隙和喉道的饱和度分布对认识数字岩心中的两相流动具有重要意义，根据结果能够表征流体在数字岩心内部的流动能力，得到限制流体流动的孔隙和喉道的直径等。

主权项

1.一种基于侵蚀‑膨胀算法的数字岩心两相流模拟结果定量表征方法，其特征在于，具体步骤包括：(1)获取数字岩心的几何结构数据文件及两相流动模拟结果文件；所述几何结构数据文件包括每个像素点及该像素点对应的几何结构，几何结构包括孔隙和固体；所述两相流动模拟结果文件是指数字岩心中每个像素点位置的两相流体的饱和度；所述两相包括水相、油相；(2)利用侵蚀‑膨胀算法提取步骤(1)所述数字岩心的孔隙网络模型，所述孔隙网络模型包括多个孔隙及多个喉道，每个孔隙包括多个孔隙像素点，每个喉道包括多个喉道像素点，获取每个孔隙的半径及体积，并获取每个喉道的半径及体积；(3)求取每个孔隙的两相流体的饱和度：每个孔隙像素点的两相流体的饱和度，并求和取平均值，该平均值即为该孔隙的两相流体的饱和度；求取每个喉道的两相流体的饱和度：求取每个喉道包括的每个喉道像素点的两相流体的饱和度，并求和取平均值，该平均值即为该喉道的两相流体的饱和度；(4)利用体积平均方法，获取每个直径相同的孔隙的两相流体的平均饱和度并获取每个直径相同的喉道的两相流体的平均饱和度即数字岩心两相流模拟结果定量表征结果；具体是指：式(I)、(II)中，i表示具有相同直径的孔隙的个数，j表示具有相同直径的喉道的个数，VPi表示相同直径的孔隙中第i个孔隙的体积，S1i表示第i个孔隙中第一相流体的饱和度，S2i表示第i个孔隙中第二相流体的饱和度；VTj表示相同直径的喉道中第j个喉道的体积，S1j表示第j个喉道中第一相流体的饱和度，S2j表示第j个喉道中第二相流体的饱和度。