The carbonate reservoir of Ordovician in the sixth—seventh blocks of the Tahe Oilfield is a kind of special reservoir, which is affected by the karstification, and having strong heterogeneity. Consequently, it is important to depict the shape and 3D distribution of this reservoir for the oil and gas field exploration and development. The method of “multipoint geostatistics” was first proposed for sixth—seventh blocks in the Tahe Oilfield from the point of vug reservoir development rule and mode in the vertical, and on the basis of the special corresponding relationship between vug reservoir and seismic attributes. According to the corresponding relationship between vug reservoir and seismic attributes, the training image was built base on the special attribute. Taking drilling, seismic and well logging interpretation as hard date, and using the cave development probability cube from seismic wave impedance as software date to restrain the interwell. The result shows that the method not only can display the shape of vug reservoir, but also reproduce the complex structure and spatial distribution of the reservoir, and the simulation results has consistence with the reservoir development rule in the vertical.