In the case of Mersenne primes $p=2^n-1$ it is known (by elementary factorisation) that the exponent $n$ is necessarily a prime. Therefore, in this special case $\lim_{k \rightarrow} \frac{|S_k|}{k}$. To the best of my knowledge the general case is not known. The closest results in the literature that I am aware of are to be found in Hooley's book (Applications of sieve methods to the theory of numbers. Cambridge Tracts in Mathematics, No. 70. 1976) [http://www.ams.org/mathscinet-getitem?mr=404173] As far as I remember, (without the book at hand), Hooley needs quite strong assumptions (one of these GRH) to say more. However, he studies the related cases of Cullen primes, $p=n 2^n+1$, and here he shows, unconditionally, that the corresponding limit is $0$. (This case is easier as one has more control over the distribution in residue classes). There is some further work by Rieger (Rieger, G. J. Über Primzahlen und dünne Folgen. Arch. Math. (Basel) 28 (1977), no. 6, 600–602. http://www.ams.org/mathscinet-getitem?mr=447169 ). The fact that there are very few results on this topic reflects the fact that such thin sequences lead to problems where current methods do not quite work.