This paper proposes a novel meshless numerical simulation method for fractured-vuggy reservoirs. Initially, model discretization involves extracting fracture and vug feature nodes along the contours of the reservoir bodies and in alignment with fracture orientations, tailored to the characteristics of the fractured-vuggy oil deposit. Subsequently, a connection element model is established within the influence domain. Based on the seepage equations of the connection system, a computational method for the parameters of the connection elements (connection transmissibility and connection volume) is defined. In addition, a sequential method is employed to solve for the pressure and saturation, achieving rapid prediction of the reservoir's production dynamics. On this basis, an automatic history matching algorithm is utilized for real-time fitting of oil–water dynamic indicators, inverting the characteristic parameters of the connection element model, and quantitatively characterizing the injection–production connection elements. The research findings indicate that the method is capable of rapidly fitting and predicting the dynamics of oil reservoir production. Furthermore, conceptual model examples have substantiated that it achieves comparable computational efficiency and superior computational accuracy when compared to the traditional finite difference (volume) method under the same nodal conditions. Additionally, the node arrangement for fractures and vugs is comparatively flexible, and it can ensure the integrity of the flow paths with a limited number of nodes to enhance predictive accuracy. Therefore, this method can effectively meet the rapid prediction requirements for fractured-vuggy reservoirs and also provides a novel meshless computational approach for the numerical simulation of such reservoirs.