The　objective　of　this　paper　is　to　study　the　well-posedness　of　solutions　for　the　three　dimensional　planetary　geostrophic　equations　of　large-scale　ocean　circulation　with　additive　noise.　Since　strong　coupling　terms　and　the　noise　term　create　some　difficulties　in　the　process　of　showing　the　existence　of　weak　solutions,　we　will　first　show　the　existence　of　weak　solutions　by　the　monotonicity　methods　when　the　initial　data　satisfies　some　＂regular＂　condition.　For　the　case　of　general　initial　data,　we　will　establish　the　existence　of　weak　solutions　by　taking　a　sequence　of　＂regular＂　initial　data　and　proving　the　convergence　in　probability　as　well　as　some　weak　convergence　of　the　corresponding　solution　sequences.　Finally,　we　establish　the　existence　of　weak　D-pullback　mean　random　attractors　in　the　framework　developed　in　[11,　25].