ee[v1_?NumericQ, v2_?NumericQ, v3_?NumericQ] :=
CountRoots[\[Sigma]q33[\[Omega], v1, v2, v3], {\[Omega], 0,
Infinity}];
Show[Table[RegionPlot3D[
ee[v1, v2, v3] == i, {v1, 0, 50}, {v2, 0, 50}, {v3, 0, 100},
Mesh -> None, LabelStyle -> (FontSize -> 14),
PlotStyle -> Directive[GrayLevel[ i/10 + 1/10, 3/10]],
BoundaryStyle -> None], {i, 0, 4}]]