A three-dimensional fluid model of a double-driver negative hydrogen ion source for CFETR neutral beam injection is developed. In this model, the magnetic filter field is generated by 16 permanent magnets, which are surrounded by a soft iron. In order to accurately describe the transportation of charged species in the presence of strong magnetic field, both the electron magnetization and ion magnetization are taken into account, and the accuracy of the model has been proved by comparing with experimental data. By employing this model, the spatial distributions of the plasma parameters have been investigated, and three methods are proposed to optimize the symmetry at the bottom of the expansion region of a double-driver source. The results indicate that by adjusting the power of Driver I while keeping the power of Driver II, the symmetry of the electron density and negative hydrogen ion density could be improved. Besides, the addition of partition causes the symmetry of the electron density and electron temperature becomes better, but the influence on the negative hydrogen ion density is limited. Finally, the application of magnetic shield can not only improve the symmetry of the electron density and negative hydrogen ion density, but also increase their densities at the bottom of the expansion region.