本发明提供了一种模拟石油在地下流动状态的方法，属于石油勘探技术领域，包括：在利用颜色梯度格子Boltzmann多相流模型模拟地下石油的流动状态时：利用中心差分法，求解颜色梯度格子Boltzmann多相流模型进出口边界格点f(i,j)的虚拟临近格点的网格变量；利用f(i,j)的虚拟临近格点的网格变量，求出颜色梯度格子Boltzmann多相流模型的连续表面力模型中的流体相场梯度▽ρN和界面曲率κ；其中，流体相场梯度▽ρN用于确定流体颜色梯度相场，界面曲率κ用于确定两相界面张力。该方法解决了现有技术中，利用颜色梯度格子Boltzmann多相流模型模拟地下石油的流动状态时，进出口边界处流体格点处物理量存在在某些计算方向上缺乏对应的信息，使得边界格点计算时失真的问题。