To study the grey photovoltaic soliton families in two-photon photovoltaic photorefractive materials under open-circuit conditions, the dynamical evolution equation and the numerical solution of the grey photovoltaic soliton families were established. The existence curves of the grey photovoltaic soliton families were simulated by the numerical method. The results show that the soliton families are derived under the assumption that the multiple carrier beams share the same polarization, wavelength, and are mutually incoherent. Such grey photovoltaic soliton families reduce to grey-grey photovoltaic soliton pairs when they contain only two components. Furthermore, such grey soliton families will reduce to dark soliton families or soliton pairs when the greyness of soltion takes zero. Relevant examples were presented where the photovoltaic photorefractive material was LiNbO3.