In this paper, we perform the first application of the hybrid method (exact low modes plus stochastically estimated high modes) for all-to-all propagators to the HAL QCD method. We calculate the HAL QCD potentials in the I = 2 ππ scattering in order to see how statistical fluctuations of the potential behave under the hybrid method. All of the calculations are performed with the 2+1 flavor gauge configurations on a 16 3 × 32 lattice at the lattice spacing a ≈ 0.12 fm and m π ≈ 870 MeV. It is revealed that statistical errors for the potential are enhanced by stochastic noises introduced by the hybrid method, which, however, are shown to be reduced by increasing the level of dilutions, in particular, that of space dilutions. From systematic studies, we obtain a guiding principle for a choice of dilution types/levels and a number of eigenvectors to reduce noise contamination to the potential while keeping numerical costs reasonable. We also confirm that we can obtain the scattering phase shifts for the I = 2 ππ system by the hybrid method within a reasonable numerical cost; these phase shifts are consistent with the result obtained with the conventional method. The knowledge that we obtain in this study will become useful for the investgation of hadron resonances that require quark annihilation diagrams such as the ρ meson by the HAL QCD potential with the hybrid method. arXiv:1904.09549v2 [hep-lat]