In the framework of the k_T-factorization QCD approach we consider the production of b quark pairs in pbar p collisions at the Fermilab Tevatron. We investigate the dependence of the b quark, B meson and decay muon differential cross sections on the different forms of unintegrated gluon distributions. The analysis also covers the azimuthal correlations between the b and bar b quarks and their decay muons. Our theoretical results agree well with recent data taken by the D0 and CDF collaborations at Tevatron. Finally, we present our predictions for muon-muon and muon-jet cross sections at the Tevatron and CERN LHC conditions.