This question sounds like an exercise: Split the sequence into I sequences of iid random variables. Apply the classical SLLN to each sequence. Recombine.
Tom: Of course it is true with exponential decay of the corellation function, but it is not easy.
The essential difficulty is that one wants to reduce the SLLN to an exponentially growing subsequence of N's. In the classical case, this is done by a martingale inequality. (Prob of a supremum of a martingale is dominated by the probability at end of the martingale.)
Once one moves away from an implicit martingale structure, then tricks have to be employed---of which the most obvious is that if the sequence of random variables is bounded, then obviously you can reduce to an exponentially growing subsequence. This point is much of the content of the paper of Lyons cited already.
Not sure that this would appear in a text book however. My sense is that these considerations are well-known.