Step 1: Use the divergence formula.
For
\[
\vec{F}=P\hat{i}+Q\hat{j}+R\hat{k},
\]
\[
\nabla\cdot\vec{F}
=
\frac{\partial P}{\partial x}
+
\frac{\partial Q}{\partial y}
+
\frac{\partial R}{\partial z}.
\]
Here,
\[
P=2x^{2}z,\qquad
Q=xy^{2}z,\qquad
R=3yz^{2}.
\]
Step 2: Differentiate and substitute the point.
\[
\frac{\partial P}{\partial x}=4xz,
\]
\[
\frac{\partial Q}{\partial y}=2xyz,
\]
\[
\frac{\partial R}{\partial z}=6yz.
\]
Thus,
\[
\nabla\cdot\vec{F}
=
4xz+2xyz+6yz.
\]
At
\[
(x,y,z)=(1,1,1),
\]
\[
\nabla\cdot\vec{F}
=
4+2+6
=
12.
\]
Hence,
\[
\boxed{12}
\]
Therefore,
\[
\boxed{(D)}
\]
is the correct answer.