I'll just add a few more references to what Ben McKay mentioned.
Phillip Griffiths and I wrote a paper on invariants of parabolic equations in 1 space variable (Characteristic cohomology of differential systems, II: Conservation laws for a class of parabolic equations, Duke Mathematical Journal 78 (1995), 531–676. MR1334205) that explains how the method of equivalence can be used to determine invariants of parabolic second order PDE. We use it to, for example, classify up to equivalence the parabolic equations that admit non-trivial conservation laws.
As Ben McKay mentioned, Jeanne Clelland's thesis solves the equivalence problem for parabolic equations in two space variables and uses it to classify the parabolic equations that admit infinitely many independent conservation laws. (Jeanne was a student of mine.)
More recently, another student of mine, Ben McMillan, worked on the equivalence problem for parabolic equations in an arbitrary number of space variables, with a view towards classifying the ones that admit a `large' number of conservation laws. He wrote two papers that are on the arXiv, arXiv:1810.00458 and arXiv:1810.02346, that you may find interesting.