Operator product expansion at finite temperature

We extend an earlier, configuration space method to find the Wilson coefficients of operators appearing in the short distance expansion of thermal correlation functions of different quark bilinears. Considering all the different correlation functions, there arise, up to dimension four, two new operators, in addition to the two appearing already in the vacuum correlation functions. They would contribute substantially to the QCD sum rules, when the temperature is not too low.