Abstract: Exploring the characterization of singular black hole spacetimes, we study the relation between energy density, curvature invariants, and geodesic completeness using a quadratic f(R) gravity theory coupled to an anisotropic fluid. Working in a metric-affine approach, our models and solutions represent minimal extensions of general relativity (GR) in the sense that they rapidly recover the usual Reissner-Nordstrm solution from near the inner horizon outwards. The anisotropic fluid helps modify only the innermost geometry. Depending on the values and signs of two parameters on the gravitational and matter sectors, a breakdown of the correlations between the finiteness/ divergence of the energy density, the behavior of curvature invariants, and the (in) completeness of geodesics is obtained. We find a variety of configurations with and without wormholes, a case with a de Sitter interior, solutions that mimic nonlinear models of electrodynamics coupled to GR, and configurations with up to four horizons. Our results raise questions regarding what infinities, if any, a quantum version of these theories should regularize.