ContourPlot3D[(x^2)+(y^2)+(z^2)==4,{x,-3,3},{y,-3,3},{z,-3,3},Mesh->False,ContourStyle->{LightBlue,Opacity[0.4]}]

cil=ContourPlot3D[(x^2)+(y^2)==1,{x,-3,3},{y,-2,2},{z,-2,2},Mesh->False,ContourStyle->{LightBlue,Opacity[0.4]}]

elipse=ContourPlot3D[(x^2)+(y^2)==1,{x,-3,3},{y,-2,2},{z,-2,2},ContourStyle->{LightBlue,Opacity[0.4]},MeshFunctions->(#3-#1&),Mesh->{{0}},MeshStyle->{Blue,Thick}]

eixoz=ParametricPlot3D[{0,0,z},{z,-2,2},PlotStyle->{Blue,Thick}];

Show[elipse,eixoz]

circ1=ParametricPlot[{Cos[t],Sin[t]},{t,0,2π},PlotStyle->{Blue,Thick},GridLines->Automatic]

circ2=ParametricPlot[{2Cos[t],2Sin[t]},{t,0,2π},PlotStyle->{Blue,Thick},GridLines->Automatic];

Show[circ2,circ1]

anel=ParametricPlot[{r Cos[t],r Sin[t]},{t,0,2π},{r,1,2},Mesh->True]

Show[circ2,circ1,anel]

RegionPlot[1≤(x^2)+(y^2)≤4,{x,-2,2},{y,-2,2},GridLines->Automatic,BoundaryStyle->{Thick,Blue}]

circ1esp=ParametricPlot3D[{Cos[t],Sin[t],0},{t,0,2π},PlotStyle->{Blue,Thick}]

Show[cil,circ1esp,eixoz]

paracirc=ParametricPlot3D[{Cos[t],Sin[t],z},{t,0,2π},{z,-2,2},PlotStyle->{LightGray,Opacity[0.4]}]

Show[paracirc,circ1esp,eixoz]

pisa=ParametricPlot3D[{Cos[t],Sin[t]+0.3z,z},{t,0,2π},{z,-2,2},PlotStyle->{LightGray,Opacity[0.4]}]

Show[pisa,circ1esp,eixoz]

pisaeq=ContourPlot3D[(x^2)+(((y-0.3 z))^2)==1,{x,-3,3},{y,-2,2},{z,-2,2},Mesh->False,ContourStyle->{LightBlue,Opacity[0.4]}]

Show[pisaeq,circ1esp,eixoz]