We report on the calculation of the two-loop penguin amplitudes in non-leptonic B decays in the framework of QCD factorization. We discuss the computation of this genuine two-loop, two-scale problem and provide details on the matching from QCD onto SCET, the evaluation of the master integrals, and the convolution of the hard scattering kernels with the light-cone distribution amplitude of the light meson. Preliminary results on the size of the two-loop correction are given.