Extracting flexibility of heterogeneous deferrable loads via polytopic projection approximation