The prime objective of the article is to present the spatial characteristics of stress distribution around a shallowly buried spherical Karst cave containing fill materials in limestone strata. Firstly, considering the external load effects, stress field in the Earth’s crust, internal filling, and the Karst landform characteristics in China, a spatial axial-symmetrical model was established. Concurrently, combining available work and the concept of elasticity, the boundary conditions were determined. Subsequently, Love displacement method was introduced, in addition to the expressions of stress components were gained. The diagram characteristics of each stress component were summarized, which are affected by various influencing factors. Finally, in order to prove the rationality of the general solution, a comparison between numerical simulation result and theoretical calculation result was carried out, as well as the maximum error between theoretical calculation value and numerical simulation value is less than 5%. Thus, the analytical solution could represent the spatial characteristics of stress distribution around a shallowly buried spherical Karst cave containing fill materials in limestone strata.