Abstract: The principle and method of solving three types of satellite gravity gradient boundary value problems by least-square are discussed in detail.And the kernel function expressions of least square solution of three geodetic boundary value problems with the observations are presented.From the results of recovery of gravity field using simulated gradient tensor data,a conclusion can be drawn that the satellite gravity gradient integral formula gotten from least square,used for recovering the gravity field,is valid and rigorous.