Oxygen vacancy formation and migration energies in the pyrochlores Gd2(ZrxTi1-x)2O7 were computed using atomistic energy minimization methods, with particular attention to the role of cation antisite disorder. The interatomic potentials recognize oxygens in distinct crystallographic sites (polarizations) as distinct species, and permit assessment of the directional anisotropy in vacancy migration energies. Oxygen vacancy migration is determined by O48f - O48f transitions along the shortest edges of the TiO6 octahedra, in agreement with previous results. Oxygen vacancies migrate towards disorder to reduce local stresses. In contrast to previous work, the VO8a and VO8a + IO8b Frenkel defects are favored over the VO8f and VO48f + IO8b Frenkel. The transition VO8a -> VO48f is energetically favored for the ordered structure, and activates the three dimensional O48f - O48f vacancy migration network by filling empty 8a sites with 48f oxygens. Local regions of order and moderate disorder must both be present for optimum oxygen conductivity. High disorder retards vacancy migration.