Irreversible electroporation (IRE), a relatively new energy-based tumor ablation technology, has shown itself in the last decade to be able to safely ablate tumors with favorable clinical outcomes, yet little work has been done on optimizing the IRE protocol to variously sized tumors. Incomplete tumor ablation has been shown to be the main reason leading to the local recurrence and thus treatment failure. The goal of the present study was to develop a general optimization approach to optimize the IRE protocol for cervical tumors in different sizes, while minimizing the damage to normal tissues. This kind of approach can lay a foundation for future personalized treatment of IRE. First, a statistical IRE cervical tumor death model was built using previous data in our group. Then, a multi-objective optimization problem model was built, in which the decision variables are five IRE-setting parameters, namely the pulse strength (U), the length of active tip (H), the number of pulses delivered in one round between a pair of electrodes (A), the distance between electrodes (D), and the number of electrodes (N). The domains of the decision variables were determined based on the clinical experience. Finally, the problem model was solved by using NSGA-II algorithm to give respective optimal protocol for three sizes of cervical tumors. The results of the study demonstrate this approach can theoretically provide the optimal IRE protocol for different sizes of tumors and may be generalizable to other types, sizes, and locations of tumors.