该文对复杂地下洞室结构的三维渗流场进行了数值分析.首先,针对传统变网格法的缺陷,用节点虚流量法求解渗流自由面.在对有自由面及渗出面渗流问题的有限元法求解时,只采用一可包括结构全区域的固定单元网格;在迭代计算的过程中,网格的范围及剖分形式始终不变,各单元的渗透矩阵只需建立一次,总体渗透矩阵也只需形成及分解一次,从而克服了传统变网格法的缺陷,也克服了其它固定网格算法在求解复杂渗流问题时有时不易收敛的现象.其次,该文借助于排水子结构技术严密地模拟了排水孔的实际渗流行为,并联合并用节点虚流量法及排水子结构技术,一方面把排水子结构纳入节点虚流量法中,使后者能适用于求解含有众多排水孔时的复杂渗流场;另一方面把节点虚流量法中的渗流虚域与实域及虚单元与实单元的概念引入排水子结构,按节点虚流量法原理严密地计算排水子结构穿过或不穿过渗流自由面时排水子结构的出口传导矩阵和右端项,从而在理论上和算法上都严密地解决了排水孔穿过或不穿过渗流自由面时渗流场的求解问题.再次,该文将上述算法运用到对地下洞室群区渗流场行为的模拟中,成功地求解了糯扎渡水电站地下厂房洞室群围岩中的复杂渗流场.