In this answer I assume you really mean "set" rather than "multiset" (as we might see more typically in statistics).
One measure of set-similarity is Jaccard similarity $J(A,B) = {{|A \cap B|}\over{|A \cup B|}}$ (also called the Jaccard index). That is, the number in both sets divided by the number in either set.
Correspondingly, the Jaccard dissimilarity between two sets is $1-J(A,B)$.
We could generalize the Jaccard similarity to more than two sets readily enough.
If $A$ is a collection of sets $A_1, A_2,...,A_n$, then $J(A) = {{|\bigcap_i A_i|}\over{|\bigcup_i A_i|}}$ and then perhaps define a measure of dissimilarity (taking on some sense of "variability") as its complement, $1-J(A)$.
(However numerous other similarity measures exist, as whuber points out; it depends what you want to measure)
You mention using entropy (by which I assume you mean something like cross-entropy).
To work with cross entropy you'd need to assign some sort of probabilities to the elements.
If the sets were finite, and one were to define the probabilities uniformly that might work, but
cross entropy is also not symmetric; you'd presumably want a symmetric measure (you could perhaps add the two cross entropies $d(p,q)=H(p,q)+H(q,p)$).
then you'd need to generalize to more than two sets, possibly by summing all the pairwise $d$s. However
I don't think this would be especially satisfactory as it stands since it's not 0 when the sets are the same.
Related to it but better still might be the symmetrized Kullback-Liebler divergence. Again you would need to generalize to multiple sets.
Hopefully these give you some ideas. You should probably look around some of the other similarity and dissimilarity indices that already exist.