We present an investigation of disordered photonic crystals (PhCs) based on the combination of photonic Wannier functions with the concept of the coherent potential approximation (CPA). In particular, we provide the theoretical foundation of a real-space cluster CPA that is causal, enforces the proper symmetries of the effective medium, and includes effects of multiple scattering of the same and nearby defects, which is essential for strong defects. Based on this, we present results for the density of states of disordered PhCs for different types of disorder. Our results are thus relevant to such diverse areas as random lasing and the analysis of fabricational imperfections in PhCs.