The thermodynamics of the vortex lattice of high-temperature superconductors has been studied by solving the generalized Ginzburg-Landau equations derived microscopically. Our numerical simulation indicates that the structure of the vortex lattice is oblique at the temperature far away from the transition temperature Tc, where the mixed s-dx2-ya state is expected to have the lowest energy. Whereas, very close to Tc, the dx2-ya wave is slightly lower energetically, and a triangular vortex lattice recovers. The coexistence and the coupling between the s and d waves would account for the unusual dynamic behaviours such as the upward curvature of the upper critical field curve Hc2(T), as observed in dc magnetization measurements on single-crystal YBa2Cu307 samples.