The following is the best probabilistic heuristic argument I know, which convinces me that either there is no blow up of the 3D Navier-Stokes equations for any initial data, or a blow up is very unlikely, or at least if there is a blow up it would be difficult to prove that it is a blow up:
See the paper "An Elementary Proof of the Existence and Uniqueness Theorem for the Navier-Stokes Equations" by Mattingly and Sinai. In that paper, they put the 3D equations in vorticity formulation and show what it looks like in Fourier space on page 12:
Let $\alpha=2$, $g_k=0$, $\mathcal Z = \mathbb Z^3$, and $k=(k_1,k_2,k_3) \in \mathbb Z^3$. Then they split the system into the equations for the "real and imaginary parts of $\{u_k\}_k$ and $\{\omega_k\}_k$. Letting $u_k(t)=u^{(1)}_k(t)+ i u^{(2)}_k(t)$, $\omega_k(t)=\omega_k^{(1)}(t)+i \omega_k^{(2)}(t)$,Letting and $g_k(t)=g^{(1)}_k(t)+ i g^{(2)}_k(t)$; we$$ \begin{split} u_k(t)&=u^{(1)}_k(t)+ iu^{(2)}_k(t),\\ \omega_k(t)&=\omega_k^{(1)}(t)+i \omega_k^{(2)}(t),\:\text{ and}\\ g_k(t)&=g^{(1)}_k(t)+ i g^{(2)}_k(t), \end{split} $$ we obtain \begin{align} \label{3DVorticityRealImag} \frac{d \omega_k^{(1)}(t)}{dt} = &2 \pi \sum_{\substack{l_1+l_2=k\\ l_1 , l_2 \in \mathcal{Z}}} \Bigl[ \bigl( u^{(1)}_{l_1}(t),k\bigr)\omega^{(2)}_{l_2}(t)+ \bigl(u^{(2)}_{l_1}(t),k\bigr)\omega^{(1)}_{l_2}(t) - \bigl( \omega^{(2)}_{l_1}(t),k\bigr)u^{(1)}_{l_2}(t) - \bigl( \omega^{(1)}_{l_1}(t),k\bigr)u^{(2)}_{l_2}(t) \Bigr] \notag \\ & \qquad -4 \pi^2 \nu |k|^\alpha \omega_k^{(1)}(t) + g_k^{(1)}(t) \\ \frac{d \omega_k^{(2)}(t)}{dt} = &-2 \pi \sum_{\substack{l_1+l_2=k\\ l_1 , l_2 \in \mathcal{Z}}} \Bigl[ \bigl( u^{(1)}_{l_1}(t),k\bigr)\omega^{(1)}_{l_2}(t)- \bigl(u^{(2)}_{l_1}(t),k\bigr)\omega^{(2)}_{l_2}(t) - \bigl( \omega^{(1)}_{l_1}(t),k\bigr)u^{(1)}_{l_2}(t) + \bigl( \omega^{(2)}_{l_1}(t),k\bigr)u^{(2)}_{l_2}(t) \Bigr] \notag \\ & \qquad + g_k^{(2)}(t) -4 \pi^2 \nu |k|^\alpha \omega_k^{(2)}(t) " \end{align}\begin{split} \label{3DVorticityRealImag} \frac{d \omega_k^{(1)}(t)}{dt} &= 2 \pi \sum_{\substack{l_1+l_2=k\\ l_1 , l_2 \in \mathcal{Z}}} \Bigl[ \bigl( u^{(1)}_{l_1}(t),k\bigr)\omega^{(2)}_{l_2}(t)+ \bigl(u^{(2)}_{l_1}(t),k\bigr)\omega^{(1)}_{l_2}(t) \\ &\qquad- \bigl(\omega^{(2)}_{l_1}(t),k\bigr)u^{(1)}_{l_2}(t) -\bigl( \omega^{(1)}_{l_1}(t),k\bigr)u^{(2)}_{l_2}(t) \Bigr]\\\\ &\qquad -4 \pi^2 \nu |k|^\alpha \omega_k^{(1)}(t) + g_k^{(1)}(t) \\\\ \frac{d \omega_k^{(2)}(t)}{dt} & = -2 \pi \sum_{\substack{l_1+l_2=k\\ l_1 , l_2 \in \mathcal{Z}}} \Bigl[ \bigl( u^{(1)}_{l_1}(t),k\bigr)\omega^{(1)}_{l_2}(t) - \bigl(u^{(2)}_{l_1}(t),k\bigr)\omega^{(2)}_{l_2}(t)\\ &\qquad - \bigl( \omega^{(1)}_{l_1}(t),k\bigr)u^{(1)}_{l_2}(t) + \bigl( \omega^{(2)}_{l_1}(t),k\bigr)u^{(2)}_{l_2}(t) \Bigr] \\\\ & \qquad + g_k^{(2)}(t) -4 \pi^2 \nu |k|^\alpha \omega_k^{(2)}(t) \end{split}
Notice that the convection part of the right hand side of the two real and imaginary part equations has an equal number of positive and negative terms. (Also note that each $u_k$ is bounded above by the square root of the total initial energy, which is a finite constant.) Because of this, it seems very difficult to find initial data which will prevent the terms of the convection part from approximately cancelling each other out at each time $t$. (The approximate cancellations would stop the enstrophy from going to infinity in finite time.) Also, the diffusion part of the right hand side of the equations can only help to prevent a blow up. Therefore, a blow up of the 3D Navier-Stokes equations seems unlikely, if not impossible.