In fractured-vuggy reservoirs， triple media consisting of fractures， vugs， or pores significantly affects oil and gas flow.However， traditional well testing interpretation methods have been proven to be ineffective in accurately describing the flow characteristics of the fluid. This paper introduced a numerical well testing method based on the fast marching method. This approach employed an embedded discrete fracture model to describe the triple media of fractured-vuggy reservoirs. The Eikonal equation was established to describe the fluid transmission process between different media and solved with the fast marching method. Then， the control equation under a one-dimensional diffusive time of flight coordinate system could be established. On this basis， two reservoir models were established according to the geological characteristics of different combinations of pores， fractures， and vugs： the vuggy reservoir model and the fractured-vuggy reservoir model. The results indicate that the size of vugs and their relative position to the wellbore lead to differences in the characteristics of the typical curve in vuggy reservoirs. When the vugs are connected to the wellbore， the pressure dynamic is influenced by different vug sizes， and the typical curve exhibits four distinct stages： wellbore storage， radial flow， late-time decline， and boundary flow. In contrast， when the vugs are not connected to the wellbore， the relative position of the vugs to the wellbore affects the pressure dynamic， and the typical curve experiences five stages： wellbore storage，radial flow， boundary flow， radial flow， and boundary flow. In the fractured-vuggy reservoirs， the typical curve undergoes four stages： wellbore storage， early-time decline， late-time decline， and boundary flow. The effective communication between vugs and natural fractures is a key factor that affects the productivity of oil reservoirs. Moreover， vugs serve as an advantageous space for fluid storage and have a greater impact on pressure dynamic changes during the whole production cycles of the reservoirs.Furthermore， the proposed numerical well testing method has been successfully applied to the dynamic analysis of actual reservoirs，demonstrating its applicability and effectiveness. Compared with traditional methods， this approach can more accurately characterize the fluid flow behavior in complex fractured-vuggy reservoirs， providing valuable insights for dynamic evaluation and development of similar reservoirs.