You haven't specified what sorts of explanations you would find "nice," so the following collection of thoughts might not be helpful, but here goes (or it might be helpful, but contain only things you already know). [By the way, this answer only attempts to "explain" why $V_K(i) = \Delta(-1)$ -- it does not address the additional question of why they both yield the knot determinant. That's a great question; I guess one "explanation" is simply that the invariants are determined by their skein relation.]

Firstly, a reference I like, in addition to the one you mention by Sartori, is this paper by Kauffman and Saleur. In that paper, many (though not all) of the theorems and computations are for the general Lie superalgebra $\mathfrak{gl}(m|n)$, where $m,n \in \mathbb{Z} >= 0$ are not necessarily equal. The Jones and Alexander cases are each specializations of this general case, with $(m,n) = (2,0)$ or $(1,1)$, respectively, so therefore a general strategy for a certain kind of "explanation" would be: observe that for certain specializations of the quantum parameter $q$, the generally-defined $\mathfrak{gl}(m|n)$ quantum link invariants, and even much of the structure of the quantum group $U_q(\mathfrak{gl}(m|n))$, depends only on $m+n$.

Here are two specific instantiations of the strategy. For the first, see page 299 (= page 7) of the paper I linked to. There the authors point out that for $q = i$, the quantum $R$ matrices for $U_q(\mathfrak{gl}(m|n))$ and $U_q(\mathfrak{gl}(m+n))$ agree. This is the essential ingredient for showing the relevant polynomial invariants agree at $q=i$, but to be complete one would want to go through the details of how those invariants are defined, in terms of the $R$ matrix.

The second is along those lines: the first link invariant defined in the Kauffman and Saleur paper has this skein relation (see p. 302):

$$ q^{n-m}P^+ - q^{-(n-m)}P^- = (q-q^{-1})P_{||} $$

(hopefully my notation is clear enough). However, one can see, by passing from a given knot to a split (1,1) tangle, say, that this invariant $P$ is forced to be zero for all knots if the quantum dimension $D_q$ of $U_q(\mathfrak{gl}(m|n))$ is zero. However, one can define another invariant $P':= P/D_q,$ which will be non-zero. The quantum dimension is $$ D_q = \frac{q^{n-m} - q^{-(n-m)}}{q-q^{-1}},$$ which is zero for all $q$ if $n=m$. Anyway, I haven't done it carefully myself, but I think it should be fairly straightforward to go from the above skein relation to the skein relation satisfied by $P'$ at $q = i$, and see that this only depends on $m+n$. (A final step would be to relate these normalizations of the invariants to those in the literature).

Let me know if I've said anything confusing or wrong, and if this helps in any way; or if you were imagining something different in terms of an explanation.