该研究是作为中、低放废物浅地层处置安全评价方法研究的一部分,而进行的非饱和带水气二相流方程及两相运移方程数值计算方法方面的研究.对于水气二相流方程,采用Garlerkin有限元方法进行数值模拟计算,对于其中的时间导数采用向后差分格式.并在已有的关于二相流数值计算方法的基础上,引入了离散Newton-Raphson迭代方法,替代Binning研究中的Picard迭代方法,用于线性化经过空间和时间离散后得到的非线性有限元方程组,从而将迭代方法的收敛阶由原来的线性(一阶)收敛提高到平方(二阶)收敛.但是在计算时发现,如果每次迭代的初值与实际值不接近的话,离散Newton方法易于出现数值解的发展现象,而改进的Picard方法则不会出现这种情况,所以迭代的初值的选择是用好离散Newton方法的关键,该文中是将前一个时刻的预测值作为后一时刻迭代的初值来进行计算的.同时在求解两相的有限元方程组时,采用两相方程并行迭代的方法,这与将两相方程联立来求解相比,大大地节省了计算时内存的占用量.