The feasibility of performing hydraulic fracturing can be assessed by analyzing the effect of skin on the overall deliverability of a reservoir. Aspects of the disclosure provide a method to determine equivalent skin in hydraulic fractured system through a semi analytical correlation obtained by numerical modeling. Two 3D models are considered: a base case without fracture and a complex case with fractures. The base case is simulated by varying matrix permeabilities while the fractured model is run with different sets of varying fracture's lengths, widths, and permeability as well as matrix permeabilities. The two cases are comparatively analyzed and skin is determined. A generalized correlation for determination of skin is derived through application of multivariate regression technique on obtained skins and their corresponding variables.