重非水相（DNAPLs）是地下水常见的有机污染物，理解其运移机制对于污染物修复具有重要意义。为解释实验中多孔介质润湿性对DNAPL的流动速度及残余饱和度的影响，使用多相流相场法，在构造的二维孔隙中模拟DNAPL液滴的下落过程。结果表明：模拟能准确刻画多相流界面的不稳定性，重现与实验室尺度相似的现象。文章提出的双界面模型相较于以往的微观模型能够更好的解释介质润湿性对DNAPL运移的影响。DNAPL的下落速度受前后界面共同影响，而饱和度主要由后界面决定。同时，润湿性对运移的影响大小受制于多孔介质的非均质性。

Dense Non-aqueous Phase Liquids (DNAPLs) are common pollutants in groundwater. Understanding the migration mechanisms is of great significance. To explain the effect of wettability of porous media on the flow rate and residual saturation of DNAPL in the experiment, we simulated the DNAPL migration via solving the basic fluid mechanics equations at pore scale and reproduced the
behaviors similar to the experimental results. The results show that the multi-interface model can better explain the migration mechanism and reproduce a phenomenon similar to the laboratory scale. The falling speed is affected by both advancing interface and rear interface while the residual saturation is mainly determined by the rear interface. In addition, we demonstrate that the effect of wettability on motion works through the heterogeneity of porous media.