Numerical simulation method to determine unsaturated hydraulic parameters of fracture

According to the fractal dimension, the disturbed distance between upper and lower fracture surfaces and the ratio of area contact (the proportion of contact area to whole area of fracture) of fracture, the successive random addition method is used to generated a rough-walled fracture with position-dependent aperture, Based on the parallel-plate approximation for suitably small cell in the fracture plane, the fracture is divided into many equally spaced cells and numerical simulations are done for the process of unsaturated seepage flow during drainage and imbibition. In the light of the simulated results of water saturation and unsaturated permeability under different capillary pressure, drainage and imbibition curves arc obtained and relationships between unsaturated permeability and capillary pressure during drainage and imbibition are also presented. The qualitative comparison of the simulated results and experimental results shows that the numerical simulation method in this paper is feasible and reliable to determine the unsaturated hydraulic parameters of rock fracture.