Nanoqueous phase liquid(NAPL) simulator is a powerful and popular mathematical model for modeling the flow and transport of non-aqueous phase liquids in subsurface,but the testing of its feasibility under water table fluctuation has received insufficient attention.The feature in a column test was tested through two cycles of water table fluctuation.The sandy medium in the column was initially saturated,and each cycle of water table fluctuation consisted of one water table falling and one rising,resulting in a drainage and an imbibition of the medium,respectively.It was found that the difference between the simulated and measured results in the first drainage of the column test was minor.However,with the propagation of the water table fluctuations,the simulation errors increased,and the simulation accuracy was not acceptable except for the first drainage in the two fluctuation cycles.The main reason was proved to be the estimation method of residual saturation used in this simulator.Also,based on the column tests,it was assumed that the resulting residual saturation from an incomplete imbibition process was a constant,with a value equal to that of the residual value resulting from the main imbibition process.The results obtained after modifying NAPL simulator with this assumption were found to be more accurate in the first cycle of water table fluctuation,but this accuracy decreased rapidly in the second one.It is concluded that NAPL simulator is not adequate in the case of LNAPL migration under water table fluctuation in sandy medium,unless a feasible assumption to estimate residual saturation is put forward.
The hysteresis of saturation-capillary pressure （S-p） relations was investigated in a fine sandy medium under consecutive drainage-imbibition cycles, which resulted from scheduled water level fluctuations. A drainage-imbibition cycle starts with a drainage process and ends with an imbibition process in sequence. The saturation and capillary pressure were measured online with time domain reflectometry （TDR） probes and T5 tensiometers, respectively. Results show that the relation between the degree of hysteresis and the number of the drainage-imbibition cycles is not obvious. However, the degree decreases with the increase of the initial water saturation of the imbibition processes in these drainage-imbibition cycles. The air-entry pressure of a sandy medium is also found to be constant, which is independent of the drainage-imbibition cycles and the initial water saturation of the drainage process. In all the imbibition processes, parameter a of the van Genuchten （VG） model decreases with the increase of the initial water saturation, which corresponds positively to the magnitude of the hysteresis.
