A method for optimizing fracture treatments includes interpreting a nominal pump schedule corresponding to a nominal value for each fracture control parameter. The method further includes interpreting environmental variables, and interpreting probability distributions for each of the environmental variables that is uncertain. The method further includes defining an objective function such as a net present value of each fracture treatment over a 365 day period following the fracture treatment. The method includes determining an optimal value for each fracture control parameter according to the objective function by determining the fracture control parameter values that yield the best mean net present value given the variability in the environmental variables as described by their probability distributions.