We apply the truncated Wigner function method to the model of degenerated optical parametric oscillators (DOPOs) with mutual optical couplings. First we found the origin of the quantum entanglement in the two DOPOs with mutual optical coupling. With injecting the squeezed vacuum mixed into the open port of the optical path, as a result, we found that the two DOPOs are more entangled.