On the $\mathcal {R} $-boundedness of solution operator families for two-phase Stokes resolvent equations