各向异性渗流渗流的数值计算时，当主渗透方向与坐标轴交叉时，系数矩阵的对角不占优往往使计算难以进行。本文从各向异性渗流控制方程本身出发，将由于主渗透方向与坐标轴不一致而产生的、易导致矩阵对角不占优的交叉项移至源汇项中，采用延迟迭代方法逐步逼近原方程的解，增强计算的稳定性，同时提高了求解结果的精度，较好的解决了带有交叉项的各向异性渗流数值求解问题。

In seepage flow simulation, the directions of principal permeation are not always coincided with the axises of global computational coordinates because of complex geologic structure of the natural porous medium. When the traditional method is applied, it can result in an illness discreted coefficient matrix and incorrect solution as a result of the fact that it is not a diagonal-dominance matrix. On the basis of anisotropic governing equation, a new method to divide the permeability tensor into two parts, diagonal elements and non-diagonal ones, then a differed correction method was employed to approach anisotropic model. The cross-partial derivative terms introduced by the anisotropic orientations were taken as source/sinks terms which can be easily handled by differed correction method. This method improves the diagonal dominance condition and raises the numerical accuracy and numerical stability.