Abstract: A mathematic model for gas flow through naturally fractured?vuggy carbonate reservoirs with complex fracture?vug connection is developed, by which the patterns of the flow rate changes in each fracture and the pressure changes in each vug during pressure decline and buildup stages are gained in this paper. Taking the reservoir with six macroscopic vugs and five fractures as an example, the computation indicates that during pressure decline stage, the pressure decline in vugs near the production well is much faster than that in vugs far away from the production well, and the vugs near the well contribute most gas flow in the initial pressure decline stage and the vugs far away from the well join in gradually later; during pressure buildup stage, the pressure in vugs and production in fractures vary widely because of the complexity of fracture?vug connection. In some of fractures backflow may appear in late stage. Also, this paper presents the influences of vug size, fracture length and fracture width on the changes of pressure in vugs and gas production in fractures