Let us leave scalar field aside for now and discuss just the geometry, which is the same as that of an extreme Reissner–Nordström (ERN) solution.
The ERN metric is type D of Petrov classification, so vanishing of Weyl invariant implies that the Weyl tensor is zero at the horizon (and not just null). Under suitable choice of NP tetrad there is only one Weyl scalar of NP formalism:
$$
\Psi_2 =\frac{m(m-r)}{r^4},
$$
which describes “Coulomb-like” gravitational field of an isolated charge. To clarify the physical meaning of this structure of Weyl tensor one could consider the stresses induced within some extended body (“gravitational compass”) by tidal accelerations of its different parts as was discussed by Szekeres here, see case (c) in this image:
So the vanishing of Weyl tensor means that the stresses from tidal accelerations on such a device would vanish when it approaches the black hole horizon.
Note, that for a non-extreme RN solutions there is still a hypersurface where Weyl tensor vanishes located at the minimum of $g_{00}$ metric component, $r=Q^2/M$ in Schwarzschild coordinates. For non-extreme RN black hole this hypersurface is spacelike and located between inner and outer horizons, while for over-extreme RN solutions it is timelike and serves as a boundary between regions where gravitational attraction toward the central singularity at larger $r$ changes to gravitational repulsion at smaller $r$.
Vanishing of Weyl tensor on the horizon also means that the near-horizon geometry is conformally flat. This near horizon limit of ERN geometry is Bertotti–Robinson spacetime $AdS_2 \times S^2$.
Weyl tensor also vanishes for large distances $r\to \infty$ and there is in fact a discrete conformal isometry that exchanges extremal black hole horizons (future and past, $\mathcal{H}^+$ and $\mathcal{H}^-$) and future and past null infinities, $\mathcal{I}^+$ and $\mathcal{I}^-$ discovered by Couch and Torrence in 1984. With ERN metric written in “isotropic” coordinates:
$$ds^2=−\left(1 +\frac{Q}{r}\right)^{−2}dt^2+\left(1 +\frac{Q}{r}\right)^2(dr^2+r^2dΩ_2^2),$$ where horizon is located at $r=0$, this conformal isometry is just an inversion:
$$
\hat{r}=\frac{Q^2}{r},\qquad d\hat{s}^2=\frac{Q^2}{r^2}ds^2,
$$
The new metric is the same as ERN but with a new radial coordinate $\hat{r}$. In recent works such as this this symmetry has been used to relate conserved quantities at horizon and at null infinity.
Now, let us discuss the scalar field of the BBMB solution. To me the intuitive reason why scalar field must diverge at the horizon is simple: the spatial distance to the horizon as measured by an integral $\int \sqrt{|g_{rr}|} dr$ along $t=\mathrm{const}$ hypersurface is infinite. And since $\phi$ is dependent only on $r$ and its SE tensor remains bounded on the horizon the field variation must also be infinite as it reaches horizon. Note, that scalar field singularity at the horizon is rather mild from the point of view of both neutral particles and particles carrying scalar charge as has been demonstrated by Bekenstein. Also, the existence of enhanced $AdS_2$ near-horizon symmetry, which maps time independent scalar field onto time dependent but still having a static SE tensor, suggests to me that there must also be a dynamical version of BBMB black hole with time dependent scalar field that sources precisely the same geometry and has a nicer behavior across horizon. There is probably a connection with “Aretakis charges” describing polynomial growth of scalar field perturbations near ERN horizon. While Aretakis' original work concerns test fields, there must also be a nonlinear version of the results.
One could also look for a version of Couch–Torrence inversion symmetry relating both metric (such as asymptotic values of NP scalars) and scalar field quantities at horizon and at spatial infinity.
All this suggests that there must be some interesting physics happening with conformal scalar field near the horizon of BBMB black hole.