In the light of φ-mapping method and topological current theory, the effect of disclination lines on the free energy density of nematic liquid crystals is studied. It is pointed out that the total Frank free energy density can be divided into two parts. One is the distorted energy density of director field around the disclination lines. The other is the saddle-splay energy density, which is shown to be centralized at the disclination lines and to be topologically quantized in the unit of kπ /2 when the Jacobian determinant of the director field does not vanish at the singularities of the director field. The topological quantum numbers are determined by the Hopf indices and Brouwer degrees of the director field at the disclination lines, i.e., the disclination strengthes. When the Jacobian determinant vanishes, the generation, annihilation, intersection, splitting and merging processes of the saddle-splay energy density are detailed in the neighborhoods of the limit points and bifurcation points, respectively. It is shown that the disclination line with high topological quantum number is unstable and will evolve to the low topological quantum number states through the splitting process.