One way to see this is to consider the fact that the vector space of rank (3,3) completely antisymmetric tensors ($ \Lambda_3^3(R^3) $) has dimension one (it's just a linear algebra exercise). Then define the tensor:
$$
M_{ijk}^{lmn} = \delta_i^{[l} \, \delta_j^m \, \delta_k^{n]} = \frac{1}{3!} \sum_{\sigma \in S_3} sgn(\sigma) \, \delta_i^{\sigma(l)} \, \delta_j^{\sigma(m)} \, \delta_k^{\sigma(n)}
$$
where we are summing over all the permutations of three numbers $\sigma$, and $sgn(\sigma)$ denotes the sign of the permutation. It's worthy to note that
$$
M_{ijk}^{lmn}=\frac{1}{3!}\begin{vmatrix} \delta_i^l & \delta_i^m & \delta_i^n\\
\delta_j^l & \delta_j^m & \delta_j^n \\ \delta_k^l & \delta_k^m & \delta_k^n \end{vmatrix}
$$
by the Leibniz formula of the determinant (http://en.wikipedia.org/wiki/Leibniz_formula_for_determinants). So we have $M \in \Lambda_3^3(R^3) $. Since $M \neq 0$, $M$ is a basis for the space $ \Lambda_3^3(R^3) $. Now consider the tensor
$$
\epsilon_{ijk} \, \epsilon^{lmn} = B_{ijk}^{lmn}
$$
Since $B \in \Lambda_3^3(R^3) $ and $M$ is a basis, there exists a constant $k$ such that
$$
B_{ijk}^{lmn} = k \, M_{ijk}^{lmn} \implies \epsilon_{ijk} \, \epsilon^{lmn} = k \, \delta_i^{[l} \, \delta_j^m \, \delta_k^{n]}
$$
Now to determine $k$ contract $\epsilon_{lmn}$ on both sides and use the fact that $\epsilon_{lmn} \, \epsilon^{lmn} =3!$ (since you sum $3!$ terms equal to one)
$$
\epsilon_{ijk} \, \epsilon^{lmn} \, \epsilon_{lmn} = 3! \, \epsilon_{ijk} = k \, \delta_i^{[l} \, \delta_j^m \, \delta_k^{n]} \epsilon_{lmn} = k \, \delta_i^{l} \, \delta_j^m \, \delta_k^{n} \, \epsilon_{[lmn]} = k \, \delta_i^{l} \, \delta_j^m \, \delta_k^{n} \, \epsilon_{lmn} = k \, \epsilon_{ijk}
$$
So we finally get $k=3!$ and
$$
\epsilon_{ijk} \, \epsilon^{lmn}=\begin{vmatrix} \delta_i^l & \delta_i^m & \delta_i^n\\
\delta_j^l & \delta_j^m & \delta_j^n \\ \delta_k^l & \delta_k^m & \delta_k^n \end{vmatrix}
$$
And you can get the identity you want by contracting. The same argument could be used in any dimension to show that
$$
\epsilon_{i_1,\dots,i_n} \, \epsilon^{j_1,\dots,j_n} = \begin{vmatrix} \delta_{i_1}^{j_1} & \dots & \delta_{i_1}^{j_n}\\
\vdots & & \vdots \\ \delta_{i_n}^{j_1} & \dots & \delta_{i_n}^{j_n} \end{vmatrix}
$$