The accurate identification of lithofacies is indispensable for reservoir parameter prediction. In recent years, the application of multivariate statistical methods has gained more and more attention in petroleum geology. In terms of the identification for lithofacies, the commonly used multivariate statistical methods include discriminant analysis and cluster analysis. Fisher and Bayesian discriminant analyses are two different discriminant analysis methods, which include intrinsic advantages and disadvantages. Given the discriminant efficiency of different methods, calculation cost, difficulty in the degree of determining the parameters, and the ability to analyze statistical characteristics of data, we put forward a new method combined with seismic information to classify reservoir lithologies and pore fluids. This method integrates the advantages of Fisher discrimination, the kernel function, and Bayesian discrimination. First, we analyze training data and search a projection direction. Then, data are transformed through Fisher transformation according to this direction and different kinds of facies can be distinguished more efficiently by exploiting transformed data than by using primitive data. Subsequently, using the kernel function estimates the conditional probability density function of the transformed variable. A classifier is constructed based on Bayesian theory. Then, the pending data are input to the classifier and the solution whose posteriori probability reaches the maximum is extracted as the predicted result at each grid node. An a posteriori probability distribution of predicted lithofacies can be acquired as well, from which interpreters can evaluate the uncertainty of the results. The ultimate goal of this study is to provide a novel and efficient lithofacies discriminant method. Tests on model and field data indicate that our method can obtain more accurate identification results with less uncertainty compared with conventional Fisher approaches and Bayesian methods.