Abstract: The MQ-B(multiquadric-Biharmonic method) proposed by Hardy as a geometry fitting method was widely used in geodesy,geophysics,etc.In this paper,we firstly expanded a MQ-B method that suitable for large scale velocity field establishment on the Reference ellipsoid.The method we expanded use the geodesicto establish the kernel function,the way we establish theKernel Functionis more reasonable than the traditional B and L differential kernel function as to Hardy's theory.Secondly,we find that this method can be extended in vector fieldas if we can calculate the projection matrix,andthe form of projection matrix is shown in this paper.Finally,this new geodesic Kernel Function MQ-B method was tested through an experiment.It shows that this new method has a same precisionasthe traditional B and L differentialkernel function method in the region of china.However,this new method has a more rational theoretical accountthanthe traditional B and L differentialkernel function method,It can be used to establish the large scale velocity field.