Heteroclinic intersections between Invariant Circles of   Volume-Preserving Maps
  We develop a Melnikov method for volume-preserving maps with codimension one invariant manifolds. The Melnikov function is shown to be related to the flux of the perturbation through the unperturbed invariant surface. As an example, we compute the Melnikov function for a perturbation of a three-dimensional map that has a heteroclinic connection between a pair of invariant circles. The intersection curves of the manifolds are shown to undergo bifurcations in homology 
H. E. Lomeli and J. D. Meiss