Grand canonical Monte Carlo (GCMC) simulations of binary Lennard-Jones mixtures in the zeolite silicalite have been used to predict the adsorption of CH4 and CF4 mixtures as a function of gas phase composition, total pressure, and temperature. For single components and mixtures, predictions of adsorption isotherms and isosteric heats are in good agreement with experiment at room temperature. Within the experimental pressure range of 0 to 17 bar, the mixtures are well described by the ideal adsorbed solution (IAS) theory. For very high loading, deviations from IAS theory appear. The configurations generated in the simulation were used to calculate sorbate-zeolite interaction energy distributions for different types of siting locations within the zeolite pores. These distributions display a pore shape related energetic heterogeneity in different regions of silicalite. Near saturation at a total loading of 12 molecules per unit cell, the shape of the observed energy distribution is relatively independent of the composition in the pore. Nevertheless, the energetic heterogeneity is responsible for a mild segregation in the adsorbed mixtures, with methane adsorbed preferentially in the silicalite zigzag channels and CF4 preferentially in the straight channels.