We firstly present two popular space discretization models of the Nonlinear Schrödinger Equation in dark solitons motion : the Direct-Discrete model and the Ablowitz-Ladik model. Applying the midpoint scheme to the space discretization models, we get two time-space discretization models: the Crank-Nicolson method and the New-Difference method. Secondly, we demonstrate that the solutions of the two space discretization models converge to the solution of the Nonlinear Schrödinger Equation. Also, we prove that the convergence order of the two time-space discretization models are O(h2+τ2) in discrete L2-norm error estimates. Finally, the numerical experiments agree well with the proven theoretical results.