We propose a method to use lattice QCD to compute the Borel transform of the vacuum polarization function appearing in the Shifman-Vainshtein-Zakharov (SVZ) QCD sum rule. We construct the spectral sum corresponding to the Borel transform from two-point functions computed on the Euclidean lattice. As a proof of principle, we compute the $s bar{s}$ correlators at three lattice spacings and take the continuum limit. We confirm that the method yields results that are consistent with the operator product expansion in the large Borel mass region. The method provides a ground on which the OPE analyses can be directly compared with non-perturbative lattice computations.