Relevance. In the process of formation, carbonate reef massifs undergo a large number of secondary changes: leaching, fracturing, recrystallization, dolomitization, which has a significant impact on the distribution of filtration-capacity properties of reservoirs. High intensity of diagenetic changes in carbonates can form zones with anomalously high filtration-capacity properties – supercollectors. The complication of the deposit by the presence of supercollectors is a significant factor of uncertainty, which largely determines the development of the object, since it is the main zone of fluid filtration. This carries, on the one hand, high risks of premature watering of producing wells and formation of zones of residual unrecoverable reserves, on the other hand, the presence of a supercollector significantly increases the productivity of wells, which allows keeping high production levels, so the problem of identifying and predicting supercollectors is relevant. Aim. Development of an approach for supercollector identification and modeling on the example of the N field. Methods. Statistical methods, geological and hydrodynamic modeling. Results and conclusions. The authors have analyzed core data for the presence and distribution of secondary transformations at the depth of wells, resulting in a three-dimensional cube of the intensity of catagenetic changes. Further, they identified the intervals with anomalously high filtration-capacitance properties using a set of different-scale studies (core, geophysical well surveys, field-geophysical survey) and constructed the curves of probability of supercollector presence in the well section. Using different boundary values of supercollector allocation on the probability curves (minimum – P90, maximum – P10 and optimum – P50), the authors obtained three variants of supercollector realization in the reservoir volume. At distribution of supercollectors in the interwell space the cube of intensity of secondary transformations was used as a trend. Based on these realizations, three permeability distribution arrays were calculated, which were further used in the hydrodynamic model. On the basis of reproduction of the field development history the authors compared the obtained cubes. The best convergence with actual production was obtained for realization P50. Deviations on accumulated liquid and oil production amounted to –9,6 and –7,7%, respectively.