This thesis investigates ultrasonic wave propagation in unconsolidated sands in the presence of different pore fluids. Laboratory experiments have been conducted in the sub-MHz range using quartz sand fully saturated with one or two liquids. Elastic wave propagation in unconsolidated granular material is computed with different numerical models: in one-dimension a scattering model based on an analytical propagator solution, in two dimensions a numerical approach using the boundary integral equation method, in three dimensions the local flow model (LFM), the combined Biot and squirt flow theory (BISQ) and the dynamic composite elastic medium theory (DYCEM). The combination of theoretical and experimental analysis yields a better understanding of how wave propagation in unconsolidated sand is affected by (a) homogeneous phase distribution; (b) inhomogeneous phase distribution, (fingering, gas inclusions); (c) pore fluids of different viscosity; (d) wettabilities of a porous medium. The first study reveals that the main ultrasonic P-wave signatures, as a function of the fraction on nonaqueous-phase liquids in initially water-saturated sand samples, can be explained by a 1-D scattering model. The next study investigates effects of pore fluid viscosity on elastic wave propagation, in laboratory experiments conducted with sand samples saturated with fluids of different viscosities. The last study concentrates on the wettability of the grains and its effect on elastic wave propagation and electrical resistivity.