（4）渗透矩阵调整法^[5-7]

渗流场有饱和区域和非饱和区域之分，饱和区中空隙介质完全被水充满，水的流动满足达西定律，非饱和区的水流动是复杂的两相问题。

本文借助Bathe(1979年)提出的单元渗透矩阵调整法，把整个渗流域分为2个子域R₁和R₂。在浸润线以下的子域R₁中，各点有大于零的流速，在浸润线以上的子域R2中，各点流速为零，其渗透系数降低至很小的值，取为饱和区渗透系数的1/1000^[8]。

（5）非稳定渗流程序结构^[9-14]

海堤渗流为非稳定渗流，本文海堤渗流分析的初始温度场选择如下：由t=0时稳态渗流计算的渗流场水头值作为下一潮位的瞬态计算初始条件。以后，每个潮位瞬态渗流计算的初值是取上一潮位的渗流场水头分布。

由于选取了整点时刻的一系列潮位高程作为输入，每个潮位的瞬态计算时，还需改变上游坝坡的水头边界条件，该值取该时刻潮位高程。每一时段计算中，采取迭代方式反复调整饱和非饱和区及逸出点位置。