A BVP approach to the reduction of spaceborne gradiometry: theory and simulations