Determination of the relationship among capillary pressure, saturation and interfacial area: a pore unit assembly approach