Carbonate stimulation design constitutes many uncertainties and challenges. Reservoir heterogeneity, in addition to the formation damage, complicates the attempt to understand matrix acidizing stimulation design. Without proper reservoir characterization and modeling, stimulation treatment uncertainty cannot be predicted nor evaluated, which poses difficulties for treatment evaluation and improvement processes. In practice, complete and representative data necessary for treatment design may not available.In a mature carbonate reservoir, the degree of uncertainties is increasing significantly. The additional uncertainties come from the different hydrocarbon recovery methods across the formation which leads to different resistance to flow for the stimulation fluid that sometimes is not taken into account in the stimulation design.The increasing number of mature carbonates reservoirs in the world calls for a methodology to evaluate the uncertainties in mature carbonate formations that consider recovery uncertainties, pressure uncertainties, and formation damage uncertainties. As operators are developing hydrocarbon potential in more complex carbonate formations, conventional matrix acidizing may not be the best treatment to achieve good production results. This paper will present the workflow that could optimize carbonate stimulation design, taking into account the uncertainties in mature fields. Stimulation strategy should not only be based on the knowledge of petrophysical data, but also on other key parameters that influence design. Several scenarios for those uncertainties will be simulated. The outcomes of this simulation scenario will be proposed for stimulation fluid selection and optimization.