## 4 The Influence of Different Factors on the Development of Temperature Field Patterns  
A quantitative analysis was conducted along path DGJ-3 perpendicular to the direction of the groundwater flow and path DGJ-4 along the direction of the groundwater flow. Figure 6 shows a schematic diagram of the path at X = 10 m. \begin{table}
\begin{tabular}{c c c c c c c} \hline \hline \multirow{2}{*}{**Soil Layer**} & \multirow{2}{*}{\begin{tabular}{c} **Density** \\ **Kg/m3** \\ \end{tabular} } & \multirow{2}{*}{\begin{tabular}{c} **Moisture** \\ **Content** \\ **\%** \\ \end{tabular} } & \multirow{2}{*}{\begin{tabular}{c} **Thermal** \\ **Conductivity** \\ **\%** \\ \end{tabular} } & \multirow{2}{*}{\begin{tabular}{c} **Specific Heat** \\ **J/(kg\cdot K)** \\ \end{tabular} } & \multirow{2}{*}{\begin{tabular}{c} **Latent Heat of** \\ **Phase Change** \\ **KJ/m3** \\ \end{tabular} } & \multirow{2}{*}{\begin{tabular}{c} **Permeability** \\ **Coefficientcm** \\ **M/d** \\ \end{tabular} } \\ \hline \multirow{3}{*}{**Artificial fill soil**} & \multirow{3}{*}{\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(1600\)} & \multirow{3}{*}{\(16.05\)} & \multirow{3}{*}{\(1.398\)} & \multirow{3}{*}{\(1420\)} & \multirow{3}{*}{\(36.32\)} & \multirow{3}{*}{\(2.97\times 10^{-2}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \multirow{3}{*}{**Medium gravel sand**} & \multirow{3}{*}{\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(2000\)} & \multirow{3}{*}{\(13.54\)} & \multirow{3}{*}{\(1.217\)} & \multirow{3}{*}{\(1410\)} & \multirow{3}{*}{\(31.48\)} & \multirow{3}{*}{\(2.86\times 10^{1}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \multirow{3}{*}{**Silly clay with silt**} & \multirow{3}{*}{\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(1820\)} & \multirow{3}{*}{\(47.6\)} & \multirow{3}{*}{\(1.485\)} & \multirow{3}{*}{\(1510\)} & \multirow{3}{*}{\(41.26\)} & \multirow{3}{*}{\(1.98\times 10^{-4}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \multirow{3}{*}{**Fine clay**} & \multirow{3}{*}{\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(2010\)} & \multirow{3}{*}{\(26.37\)} & \multirow{3}{*}{\(1.537\)} & \multirow{3}{*}{\(1560\)} & \multirow{3}{*}{\(56.24\)} & \multirow{3}{*}{\(1.50\times 10^{-4}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \multirow{3}{*}{**Medium sand**} & \multirow{3}{*}{\
\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(2020\)} & \multirow{3}{*}{\(17.92\)} & \multirow{3}{*}{\(1.756\)} & \multirow{3}{*}{\(1360\)} & \multirow{3}{*}{\(32.87\)} & \multirow{3}{*}{\(3.98\times 10^{0}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \multirow{3}{*}{**Gravelly clay**} & \multirow{3}{*}{
\begin{tabular}{c} **Unfrozen soil** \\ **Frozen soil** \\ \end{tabular} } & \multirow{3}{*}{\(1890\)} & \multirow{3}{*}{\(31.98\)} & \multirow{3}{*}{\(1.442\)} & \multirow{3}{*}{\(1580\)} & \multirow{3}{*}{\(48.25\)} & \multirow{3}{*}{\(5.78\times 10^{-4}\)} \\ \cline{3-3} \cline{5-6}  & & & & & & \\ \hline \hline \end{tabular}
\end{table} 
Table 1: Thermal physical parameters of each soil layer. Figure 4: Grid Partition after Increasing the Smallest Cell Size. Figure 5: Boundary conditions of new pipe curtain freezing method. In order to study the influence of different seepage velocities on the temperature field of the new pipe curtain freezing method under seepage conditions, different seepage velocities were set according to the step length for multiple numerical simulations. The initial seepage velocity was 0.27 m/d, increasing by 1 m/d each time.