A method for analyzing a rock sample includes segmenting a digital image volume corresponding to the rock sample, to associate voxels in the digital image volume with pore space or solid material. A distance transform is applied to each pore space voxel. The distance transform assigns a distance value to the pore space voxel specifying distance from the pore space voxel to a solid material voxel. Drainage is numerically simulated by, for a pore space, selecting each distance value assigned to a pore space voxel that is greater than a predetermined threshold value to represent a radius of a sphere of a non-wetting fluid introduced into the pore space. The sphere is centered at the pore space voxel corresponding to the distance value. The digital image volume is numerically analyzed to characterize a material property of the rock sample at a non-wetting fluid saturation produced by the drainage.