Radial basis functions quasi-interpolation is very useful tool for the numerical solution of differential equations, since it possesses shape-preserving and high-order approximation properties. Based on multiquadric quasi-interpolations, this study suggests a meshless symplectic procedure for KdV equation. The method has a number of advantages over existing approaches including no need to solve a resultant full matrix, accuracy and ease of implementation. We also present a theoretical framework to show the conservativeness and convergence of the proposed method. As the numerical experiments show, it not only offers a high order accuracy but also has a good property of long-time tracking capability.