My course notes give the Taylor series up to degree 2 as $$f(x+h,y+k)\approx f(x,y)+\frac{\partial f}{\partial x}h+\frac{\partial f}{\partial y}k+\frac{1}{2}\frac{\partial^2 f}{\partial x^2}h^2+\frac{\partial^2 f}{\partial x \partial y}hk+\frac{1}{2}\frac{\partial^2 f}{\partial y^2}k^2.$$
It looks like the formula could be something like $$f(x+h,y+k)=\sum_{i=0}^\infty \left(\sum_{j=0}^i \left(\frac{\partial f^i}{\partial (x^{i-j}y^j)}\cdot\frac{h^{i-j}k^j}{i!} \right)\right),$$
except, in degree 2, the factorial doesn't seem to apply equally to the mixed term, and I can't see what the pattern would be without more terms.
It is indeed applied to the miwed term, but to see it plainly, you have to use a symbolic binomial expansion: $$f(x+h,y+k)=\sum_{n=0}^\infty \frac 1{n!}\left(\frac{\partial f} {\partial x}\,h+\frac{\partial f}{\partial y}\,k\right)^{\mkern -5mu(n)},$$ where we set $$\left(\frac{\partial f} {\partial x}\,h+\frac{\partial f}{\partial y}\,k\right)^{\mkern -5mu(n)}=\sum_{i=0}^n\binom ni\frac{\partial^n\mkern-2mu f}{\partial x^i\,\partial y^{n-i}}\,h^i k^{n-i}.$$
For instance, the terms of degree $2$ are
$$\frac1{2!}\biggl(\frac{\partial f}{\partial x}\,h+\frac{\partial f}{\partial y}\,k\biggr)^{\mkern -5mu(2)}= \frac12\biggl(\frac{\partial^2\mkern-2mu f}{\partial x^i}\,h^2 + 2\frac{\partial^2\mkern-2mu f}{\partial x\,\partial y } \,h \, k + \frac{\partial^2\mkern-2mu f}{\partial y^2}\,k^{2}\biggr),$$ which are exactly what you have, after a simplification.