In　this　paper,　the　problem　of　the　global　L-2　stability　for　large　solutions　to　the　nonhomogeneous　incompressible　Navier-Stokes　equations　in　3D　bounded　or　unbounded　domains　is　studied.　By　delicate　energy　estimates　and　under　the　suitable　condition　of　the　large　solutions,　it　shows　that　if　the　initial　data　are　small　perturbation　on　those　of　the　known　strong　solutions,　the　large　solutions　are　stable.