Heat kernel asymptotics for Laplace type operators and matrix KdV hierarchy