A new Bayesian approach for multiple satellite faults detection and exclusion is proposed by introducing a classification variable to each satellite observation. If we treat this classification variable as random and assume a prior distribution for it, then a rule for satellite fault detection and exclusion based on the posterior probabilities of the classification variables is constructed under the framework of Bayesian hypothesis testing. Secondly, the Gibbs sampler is introduced to compute the posterior probabilities of the classification variables. Then the implementation for a Bayesian Receiver Autonomous Integrity Monitoring (RAIM) algorithm is designed with the Gibbs sampler. Finally, different schemes are designed to evaluate the performance of the new Bayesian RAIM algorithm in the case of multiple faults. We compare the method in this paper with the Range Consensus (RANCO) method. Experiments illustrate that the proposed algorithm in this paper is capable of detecting and eliminating multiple satellite faults, and the probability of correctly detecting faults is high.