The fuzzy logic method has been used widely in civil and earthquake engineering, but there is no comprehensive point of view for utilizing fuzzy approach in order to obtain ground motion prediction equations (GMPEs) for probabilistic seismic hazard analysis (PSHA). Hence, fuzzy magnitude–distance method as a new approach for choosing GMPEs in the process of PSHA, is developed in this research through the selection of the ruling peak ground acceleration (PGA) of each common cell (the combined cell of earthquake intensity and site to source distance). The presented method reduces the need for engineering judgments in seismic analysis based on a newly developed benchmark. It enables designers to not only determine the range of acceptable fuzzy results but also introduces a concept which ensures the selection of initial well-suited GMPEs for the analysis.