The multipole expansion is a good approximation to allow us to deal with complicated shapes, by adding terms and “complicacies” step by step: point charge + sphere + ellipsoid + ….
Making use of a perturbation theory without multipole expansion-like terms would make the description of the system (mathematically speaking) challenging. Having to deal with all the distance vectors in a complicated shape system is nearly impossible.