Motivated by the possible large annihilation contributions implied by recent CDF and LHCb measurements on nonleptonic annihilation B-meson decays, and the refined experimental measurements on hadronic B-meson decays, we study the strength of annihilation contributions within QCD factorization (QCDF) in this paper. With the available measurements of two-body B u,d,s → ππ, π K , KK decays, a comprehensive fit on the phenomenological parameters Xwhich are used to parameterize the endpoint singularity in annihilation amplitudes is performed with the statistical χ 2 approach. It is found that (1) flavor symmetry breaking effects are hardly to be distinguished between X i A,s and X i A,d due to the large experimental errors and theoretical uncertainties, where X i A,s and X i A,d are related to the nonfactorization annihilation contributions in B s and B u,d decay, respectively. So X i A,s X i A,d is a good approximation by now. (2) In principle, parameter X f A which is related to the factorization annihilation contributions and independent of the initial state can be regarded as the same variable for B u,d,s decays. (3) Numerically, two solutions are found, one is (ρ i A , φ i A [ • ]) = (2.98 +1.12 −0.86 , −105 +34 −24 ) and (ρ f A , φ f A [ • ]) = (1.18 +0.20 −0.23 , −40 +11 −8 ), the other is (ρ i A , φ i A [ • ]) = (2.97 +1.19 −0.90 , −105 +32 −24 ) and (ρ f A , φ f A [ • ]) = (2.80 +0.25 −0.21 , 165 +4 −3 ). Obviously, nonfactorization annihilation parameter X i A is generally unequal to factorization annihilation parameter X f A , which differs from the traditional treatment. With the fitted parameters, all results for observables of B u,d,s → ππ, π K , KK decays are in good agreement with experimental data.