Earthquake‐related hydrochemical changes in thermal springs have been widely observed; however, quantitative modeling of the reactive transport process is absent. In the present study, we apply reactive transport simulation to capture the hydrochemical responses in a thermal spring following the Wenchuan Ms 8.0 and Lushan Ms 7.0 earthquakes. We first constrain deep reservoir geothermal fluid compositions and temperature by multicomponent geothermometry, and then a reactive geochemical transport model is constructed to reproduce the hydrochemical evolution process. The results show that the recharge from the shallow aquifer increases gradually until it reaches a peak because of the permeability enhancement caused by the Lushan earthquake, which may be the mechanism to explain the earthquake‐related hydrochemical responses. In contrast to the postseismic effect of the Wenchuan earthquake, the chemical evolution can be considered as hydrochemical anomalies related to the Lushan earthquake. This study proves that the efficient simulation of reactive transport processes is useful for investigating earthquake‐related signals in hydrochemical time series.