针对裂隙岩体渗流非连续、非均质和各向异性的特点，提出裂隙岩体非稳定渗流的复合单元算法，给出详细的计算格式和数值处理的关键技术，包括复合单元的前处理、求解自由面边界流量项的积分方法以及溢出边界的判断准则。该算法的特点是：(1)可准确地离散岩体内任意分布的裂隙；(2)可考虑岩块的渗透性及其与裂隙间的水量交换；(3)可与传统有限单元计算程序无缝整合；(4)复合单元的拓扑信息由裂隙面和常规有限单元的交切计算生成，由于事先未考虑裂隙的存在，极大地减轻了人工剖分网格的工作量；(5)对裂隙岩体非稳定渗流过程中裂隙与岩块渗流不同步的现象具有较强的描述能力，并能提供近似于有限单元法的计算精度。通过典型砂槽试验验证了该算法的正确性和有效性。最后，对裂隙岩体进行非稳定渗流分析，计算结果很好地反映了岩体内水压力的传递过程，实现了对裂隙与岩块因水力特性不同而导致的不同于整体流向的局部流动的模拟。