This work describes simulation of well productivity reduction due to near-well flow restrictions. The flow restrictions considered include: mechanical skin, partial completion, altered formation due to residual completion fluid or liquid dropout, and non-Darcy flow effects. It addresses two specific issues. The first is how to predict changes in well performance induced by changes in near-well parameters, as needed in history matching and well test analysis. A subset of this issue is how to estimate the impact of using well parameters to simplify simulation models such as partial completion skin and non-Darcy flow skin, instead of using more layers and solving the non-Darcy flow equation. The second issue is how to accurately preserve well productivity in a large grid block, as in the case of well productivity reduced by liquid dropout in a gas condensate reservoir model.To resolve these issues, this work uses direct simulation and heuristic total skin equations developed for one-dimensional and two-dimensional multilayered cases with crossflow. The equation correctly represents the effects of nonlinear interactions among near-well flow restrictions on well productivity, resulting in an effective total skin that is much larger than the conventional linear sum of individual skin components. Numerical examples show the total skin equation can be used to predict well performance and to accurately preserve well productivity in a field-scale gas condensate reservoir model. The latter use of the equation shows a feedback mechanism, wherein the total skin is used to adjust the coarse-grid well productivity.