We introduce ensembles of repelling charged particles restricted to a ball in a non-archimedean field (such as the $p$-adic numbers) with interaction energy between pairs of particles proportional to the logarithm of the ($p$-adic) distance between them. In the {em canonical ensemble}, a system of $N$ particles is put in contact with a heat bath at fixed inverse temperature $beta$ and energy is allowed to flow between the system and the heat bath. Using standard axioms of statistical physics, the relative density of states is given by the $beta$ power of the ($p$-adic) absolute value of the Vandermonde determinant in the locations of the particles. The partition function is the normalizing constant (as a function of $beta$) of this ensemble, and we identify a recursion that allows this to be computed explicitly in finite time. Probabilities of interest, including the probabilities that specified subsets will have a prescribed occupation number of particles, and the conditional distribution of particles within a subset given a prescribed occupation number, are given explicitly in terms of the partition function. We then turn to the {em grand canonical ensemble} where both the energy and number of particles are variable. We compute similar probabilities to those in the canonical ensemble and show how these probabilities can be given in terms the canonical and grand canonical partition functions. Finally, we briefly consider the multi-component ensemble where particles are allowed to take different integer charges, and we connect basic properties of this ensemble to the canonical and grand canonical ensembles.