The invention relates to a numerical
simulation method for multiphase fluid stratum seepage, which belongs to the technical field of numerical
simulation and comprises the following steps: S1, generating an unstructured grid by adopting a mode of combining a
leading edge propulsion method and
Delaunay triangulation; s2, performing dynamic adjustment through an adaptive grid
encryption technology according to the change of the seepage
physical quantity; s3, establishing a dynamic
phase change model, and obtaining a function of representing the
phase change rate as temperature, pressure and concentration of each phase; s4, the
mass and energy change of each phase in the
phase change process is calculated in real time through a kinetic equation; s5,
processing dynamic evolution of a phase interface by adopting a phase field method, and describing distribution conditions of different phases through phase field variables; and S6,
coupling the
grid system obtained by dividing the multi-scale unstructured grid with the dynamic phase
change model through a numerical solution
algorithm. According to the method, the complex geometrical shape of the stratum can be better fitted through multi-scale unstructured grid division, and then the change of the phase change process along with time and space can be tracked in real time through the dynamic phase
change model.