A reconstruction algorithm is proposed, based on multi-dictionary learning (MDL), to improve the reconstruction quality of acoustic tomography for complex temperature fields. Its aim is to improve the under-determination of the inverse problem by the sparse representation of the sound slowness signal (i.e., reciprocal of sound velocity). In the MDL algorithm, the K-SVD dictionary learning algorithm is used to construct corresponding sparse dictionaries for sound slowness signals of different types of temperature fields; the KNN peak-type classifier is employed for the joint use of multiple dictionaries; the orthogonal matching pursuit (OMP) algorithm is used to obtain the sparse representation of sound slowness signal in the sparse domain; then, the temperature distribution is obtained by using the relationship between sound slowness and temperature. Simulation and actual temperature distribution reconstruction experiments show that the MDL algorithm has smaller reconstruction errors and provides more accurate information about the temperature field, compared with the compressed sensing and improved orthogonal matching pursuit (CS-IMOMP) algorithm, which is an algorithm based on compressed sensing and improved orthogonal matching pursuit (in the CS-IMOMP, DFT dictionary is used), the least square algorithm (LSA) and the simultaneous iterative reconstruction technique (SIRT).