Commonly, capillary pressure - saturation - relative permeability (P-c - S - K-r) relationships are obtained by means of laboratory experiments carried out on soil samples that are up to 10 - 12 cm long. In obtaining these relationships, it is implicitly assumed that the soil sample is homogeneous. However, it is well known that even at such scales, some micro-heterogeneities may exist. These heterogeneous regions will have distinct multiphase flow properties and will affect saturation and distribution of wetting and non-wetting phases within the soil sample. This, in turn, may affect the measured two-phase flow relationships. In the present work, numerical simulations have been carried out to investigate how the variations in nature, amount, and distribution of sub-sample scale heterogeneities affect P-c - S - K-r relationships for dense non-aqueous phase liquid (DNAPL) and water flow. Fourteen combinations of sand types and heterogeneous patterns have been defined. These include binary combinations of coarse sand imbedded in fine sand and vice versa. The domains size is chosen so that it represents typical laboratory samples used in the measurements of P-c - S - K-r curves. Upscaled drainage and imbibition P-c - S - K-r relationships for various heterogeneity patterns have been obtained and compared in order to determine the relative significance of the heterogeneity patterns. Our results show that for micro-heterogeneities of the type shown here, the upscaled P-c - S curve mainly follows the corresponding curve for the background sand. Only irreducible water saturation ( in drainage) and residual DNAPL saturation ( in imbibition) are affected by the presence and intensity of heterogeneities.