Accurate quantification of propeller excitation is essential for effective vibro-acoustic analysis and mitigation in marine vessels. This study introduces a Bayesian statistical framework that indirectly reconstruct uncertain forces and system parameters by fusing available prior information, forward modeling and measured vibration responses through Bayesian theorem. The Markov Chain Monte Carlo (MCMC) algorithm, enhanced with Polynomial Chaos Kriging (PC-Kriging) surrogates, is employed to facilitate efficient estimation of likelihood functions and robust uncertainty propagation, providing posterior probabilistic estimates and credible intervals for the identified forces. Additionally, the framework incorporates model class selection based on model evidence to determine the most plausible empirical spectrum model for the propeller broadband thrust spectrum. Specifically, this study proposes and evaluates two candidate empirical spectrum models, the Ornstein-Uhlenbeck (OU)-Gaussian and OU-Weibull models, which effectively capture the frequency- amplitude characteristics of unsteady thrust. Numerical results demonstrate substantial parameter uncertainty reduction, accurate thrust spectrum reconstruction, and effective quantification of statistical properties. Moreover, the OU-Gaussian model is revealed to be more plausible for characterizing the propeller thrust spectrum and response prediction. This framework may offer a practical solution for propeller excitation identification, enabling rapid vibro-acoustic performance assessments of marine vessels.