Abstarct: Vuggy porosity, which is a porosity type in carbonate rocks, significantly effects on permeability, pressure drop and recovery factor, then; its identification and modeling is important for reservoir characterization and history matching. Conventionally, Image logs, RCAL and SCAL tests are main instruments to determine vugs, which all are costy and usually unavailable. Therefore developing a generalized approach that distinguishes Vugs will be useful. In this paper, a combined Parzen-wavelet-baxsed algorithm is developed to identify Vuggy zones in wavelet coefficient domain, using gamma ray (GR), neutron porosity (NPHI), bulk density (RHOB) and sonic (DT) logs. Compatibility between results of core tests and this newly developed method, confirms applicability of the proposed method.