The present paper addresses the existence of J2 invariant relative orbits with arbitrary relative magnitude over the infinite time using the Routh reduction and Poincare techniques in the J2 Hamiltonian problem. The current research also proposes a novel numerical searching approach for J2 invariant relative orbits from the dynamical system point of view. A new type of Poincare mapping is defined from different central manifolds of the pseudo-circular orbits (parameterized by the Jacobi energy E, the polar component of momentum Hz and the measure of distance Δr between the fixed point and its central manifolds) to the nodal periods Td and the drifts of longitude of the ascending node during one period (ΔΩ), which differs from Koon et al.’s (AIAA 2001) definition on central manifolds parameterized by the same fixed point. The Poincare mapping is surjective because it compresses the three-dimensional variables into two-dimensional images, and the mapping degenerates into a bijective mapping in consideration of the fixed points. An iteration algorithm to the degenerated bijective mapping is proposed from the continuation procedure to perform the ergodic representation of E- and Hz-contour maps on the space of Td–ΔΩ. For the surjective mapping with Δr ≠ 0, different pseudo-circular or elliptical orbits may share the same images. Hence, the inverse surjective mapping may achieve non-unique variables from a single image, which makes the generation of J2 invariant relative orbits possible. The pseudo-circular or elliptical orbits generated from the surjective mapping will be defined in different meridian planes. Hence, the critical contribution of the present paper is the assignment of J2 invariant relative orbits to different invariant parameters E and Hz depending on the E- and Hz-contour map, which will hold J2 invariant relative orbits for extended durations. To investigate the high-order nonlinearity neglected by previous studies, a formation configuration with a large magnitude of 500 km is successfully generated from the theory developed in the present work, which is beyond the scope of the linear conditions of J2 invariant relative orbits. Therefore, the existence of J2 invariant relative orbit with an arbitrary relative magnitude over the infinite time is achieved from the dynamical system point of view.