The current methods to optimize the contact structure among main transformer stations in distribution network based on power supply capability, only considers the interconnections of contact corridors among main transformers, but ignores the interconnections among feeders. In view of the shortages of existing methods, a method based on the Pareto optimality of the main transformers contact structure optimization is proposed. First, in order to improve total supply capacity of the distribution area and reduce construction cost which takes geographic factors into account, a multi-objective optimization model based on the feeder interconnection matrix for the contact structure of main transformers between substations is built. Furthermore, the NSGA-Ⅱ algorithm is used to solve the model and the obtained individuals of Pareto frontier are taken as alternative options. To sort the alternative options, the variation coefficient method and the TOPSIS method are used and the optimal contact structure planning is screened out. The resulting optimization not only contains the contact corridors between the main transformers, but also shows the composition of the number and location of contact branches clearly. The validity and effectiveness of the proposed method is verified by the analysis of a practical example. This work is supported by National Natural Science Foundation of China (No. 51177009).