We investigate the interaction between $Omega$ baryons in the $^1S_0$ channel from 2+1 flavor lattice QCD simulations. On the basis of the HAL QCD method, the $OmegaOmega$ potential is extracted from the Nambu-Bethe-Salpeter wave function calculated on the lattice by using the PACS-CS gauge configurations with the lattice spacing $asimeq 0.09$ fm, the lattice volume $Lsimeq 2.9$ fm and the quark masses corresponding to $m_pi simeq 700$ MeV and $m_Omega simeq 1970$ MeV. The $OmegaOmega$ potential has a repulsive core at short distance and an attractive well at intermediate distance. Accordingly, the phase shift obtained from the potential shows moderate attraction at low energies. Our data indicate that the $OmegaOmega$ system with the present quark masses may appear close to the unitary limit where the scattering length diverges.