Tahe Oilfield is a marine carbonate fractured and vuggy reservoir with a very special type of reservoir. The fractured and vuggy reservoirs have diverse types of reservoirs, complex spatial structures, and fluid flow types including not only seepage, but also fracture flow and cave flow. To solve the problem of fluid flow in this complex medium, this paper establishes a mathematical model of dual porosity and dual permeability black oil, and presents a joint simultaneous solution method of fracture and matrix for the first time for the numerical simulation of three-dimensional two-phase micro compressible fluid in a fractured vuggy reservoir. The finite volume dispersion is applied, the control volume is a polyhedron, and the entire grid is an arbitrary unstructured grid. This method is used to carry out numerical simulation research on a unit that has been proved to be connected, preliminarily complete the establishment of the reservoir numerical model, and the matching of the production performance history of the unit and single well and the remaining oil distribution law match the actual well. It will lay a foundation for the numerical simulation of fracture cavity units of different reservoir models in the next stage.