We study thermodynamic properties of 2 + 1 flavor QCD applying the Small Flow-time eXpansion (SFtX) method based on the gradient flow. The method provides us with a general way to compute correctly renormalized observables irrespective of explicit violation of symmetries due to the regularization, such as the Poincaré and chiral symmetries on the lattice. We report on the status of our on-going project to compute the energy-momentum tensor and the chiral condensate at the physical point with improved Wilson quarks, extending our previous study with slightly heavy u and d quarks. We also report on our test of two-loop matching coefficients recently calculated by Harlander et al., revisiting the case of QCD with slightly heavy u and d quarks. Our results suggest that the SFtX method is powerful in extracting physical observables on the lattice.