By employing the perturbative QCD (pQCD) factorization approach, we calculate the full leading and the partial next-to-leading order (NLO) contributions to the seven B →πη(') and η(')η(') decays. For B+ →π+η(') decays, the pQCD predictions for their decay rates agree very well with the data after the inclusion of the small NLO contributions. For neutral decays, the pQCD predictions are also consistent with the experimental upper limits and can be tested by the LHC experiments. The measured value of ACdiPr(π±η)=-19±7% can also be accommodated by the pQCD approach.