En este trabajo, investigamos el comportamiento de las trayectorias de baja energía en el marco dinámico del problema restringido de 4 cuerpos elíptico espacial, desarrollado utilizando el formalismo hamiltoniano. Introduciendo transformaciones canónicas, la función hamiltoniana en la vecindad del punto de libración colineal (o ), puede expresarse como una suma de tres integrales locales de movimiento de segundo orden, que proporcionan una descripción topológica compacta de los tránsitos de baja energía, capturas y órbitas de puntos de libración cuasiperiódicas, además de términos de orden superior que representan perturbaciones. El problema de los denominadores pequeños se aplica entonces al orden tres de la función hamiltoniana transformada, para identificar los efectos de la resonancia orbital de los primarios sobre las órbitas cuasiperiódicas. Se determinan soluciones estacionarias para estos términos resonantes, correspondientes a órbitas cuasiperiódicas que existen en presencia de resonancia orbital. El modelo propuesto se aplica al sistema Júpiter-Europa-Io, determinando órbitas cuasiperiódicas en los alrededores de Júpiter-Europa considerando la resonancia orbital 2:1 entre Europa e Io.