**4. Application of the Calculation Method**

In this section, the calculation method is applied to investigate the influence of overconsolidation effect on thermal pore water pressure using MATLAB software. The consolidation pressure is fixed at 300 kPa, and thermal loadings are taken as 35 ◦C, 50 ◦C, 65 ◦C, and 80 ◦C, respectively. Furthermore, the required parameters for numerical calculation are the same as those in Table 1. The evolution law of thermal pore water pressure with *OCR* under different thermal loadings is shown in Figure 3.

**Figure 3.** Evolution of the thermal pore water pressure with *OCR* under different thermal loadings.

From Figure 3, it can be observed that the value of thermal pore water pressure drops sharply with increasing *OCR*. Meanwhile, the difference of thermal pore water pressure under different thermal loadings also decreases gradually. This is mainly because the thermal pore water pressure is generated by the difference in thermal expansion characteristics of soil particles and pore water, and the degree of overconsolidation of soil can reduce this difference of thermal expansion. In addition, it can be found from Figure 3 that when thermal loading is large and *OCR* is small, the value of thermal pore water pressure can exceed the effective stress of the soil. In other words, the soil may suffer thermal damage. In this case, if *OCR* is properly increased, the thermal pore water pressure will be significantly reduced, thus avoiding the occurrence of thermal damage.
