Analytical Solution of Steady-State Temperature Field of Single Freezing Pipe under Action of Seepage Field

To accurately describe the distribution law of the temperature field formed by a single freezing pipe under the action of a seepage field, the shape of the freezing front was simplified using a segmentation-equivalent method.)e analytical solution of the steadystate temperature field was derived, and the accuracy was verified using a physical model test. Combined with the results of the model test and the calculation results of the analytical solution, the distribution law of the freezing temperature field formed by a single pipe under different seepage velocities was analyzed. It was found that compared with the no flow rate, when the seepage velocity was 3, 6, and 9m/day, the frozen area was reduced from 17.97×10mm to 15.77×10, 3.84×10, and 3.05×10mm, respectively. )e proportion of frozen area below −5°C increased from 39.43% to 40.19%, 49.84%, and 51.52%, respectively. )e average freezing temperature field reduced from −5.78 to −5.86, −7.31, and −7.50°C, respectively. As the seepage velocity increased, the frozen area formed by a single pipe decreased while the proportion of the low-temperature zone increased and the average temperature of the temperature field decreased.


Introduction
e artificial formation freezing method was first proposed by Poetsch in the 19th century and was originally used in the shaft construction process in soft and water-rich soil layers. During the implementation of this method, an external refrigeration unit first reduces the low-temperature refrigerant to a negative temperature (the temperature of the brine or alcohol is generally reduced to −40°C, and the temperature of the liquid nitrogen can be reduced to −190°C), and a low-temperature refrigerant is subsequently pumped through each freezing pipe. e low-temperature refrigerant in the freezing pipe exchanges heat with the surrounding soil layer through the pipe wall, thereby freezing the soil layer [1][2][3][4][5]. is method has a little influence on the surrounding environment, and the formed frozen wall exhibits a good sealing performance and high strength as a temporary supporting structure. erefore, this method has gradually developed into one of the main methods used in underground construction in complex strata. In recent years, in the construction of subway tunnels in coastal cities, such as Tianjin, Guangzhou, Shenzhen, and Ningbo, the freezing method has been widely used in the soil reinforcement of subway communication channels and the soil reinforcement of the import and export of shield tunnels [6][7][8][9]. Due to the influence of marine sedimentary environments, seawater impregnation, and tidal systems, the groundwater in coastal urban strata has high salinity and flow rate. When the artificial ground freezing method is adopted in these types of strata, a large difference exists in the freezing temperature field distribution law in the nonseepage flow field.
Aiming to understand the influence of the seepage field on the artificial freezing temperature field, scholars worldwide have conducted considerable related research. Many experimental studies have been conducted. Wang et al. used a model test method to experimentally study the effect of the groundwater flow rate on the freezing temperature field under freezing conditions created using liquid nitrogen (temperatures −57 and −80°C, respectively); the results showed that the groundwater with a flow rate of 10 m/d or more had a significant effect on the freezing by use of liquid nitrogen [10]. Zhou et al. carried out an orthogonal model test of double-pipe freezing and studied the effects of the seepage velocity and pipe spacing on the closure time of the frozen wall in saturated sand layers. Furthermore, they studied the development behavior of the temperature fields upstream and downstream [11]. Pimentel et al. summarized the results of the formation behavior of the freezing temperature field under the action of groundwater. e problem of the loss of freezing capacity was fully considered, and a new test device was designed based on the predecessor devices. Large-scale model tests of artificial ground freezing under seepage conditions with flow rates of 0, 1, 1.5, 2.0, and 2.1 m/d were carried out. Based on the test results, several analytical solutions of the closure time of the frozen wall were discussed [12]. Huang et al. studied the formation behavior of the single-pipe freezing temperature field under different flow rates (the maximum flow rate was 2 m/d) [13]. Based on the similarity theory, Li et al. established an interlaced double-row pipe freezing model test system under the action of seepage and conducted orthogonal tests on the main influencing factors for the formation of the frozen wall in the seepage stratum [14].
Several researchers have studied this problem mathematically or numerically. Victor derived a mathematical expression of the freezing temperature field of a single freezing pipe with seepage, which is the only analytical expression that accounts for the influence of the seepage velocity on the temperature field [15]. Harlan first proposed a coupled model of hydrothermal migration, and many scholars have used this result as the basis of further studies [16]. Trupak and Bakholdin proposed an analytical solution for the frozen temperature field formed by single pipe, single-row, and double-row freezing pipes [17,18]. Sanger and Sayles studied the distribution formula of freezing temperature field formed by a single-row pipe [19]. Tobe and Akimoto derived an analytical solution to a multipiped frozen temperature field [20]. Lai et al. obtained the governing differential equations for the problem with coupled temperature and seepage fields [21]. Yang and Pi established a mathematical model for the development of the freezing front of a single freezing pipe using the theory of heat transport in porous media and Darcy's law and analyzed the variation of the temperature field and groundwater flow field during the freezing process [22]. Xu obtained a coupled mathematical model of the temperature-seepage field including phase changes of the lowtemperature rock mass and applied the model to the analysis of engineering problems [23]. Hu et al. established a set of methods for solving the steady-state temperature field of artificial ground freezing based on the principle of superposition of potential functions. Combining this principle with mathematical methods, a series of analytical results for the steady-state temperature field were obtained [24][25][26][27][28]. In order to benchmark flow and energy transport models that include pore water phase change, Kurylyk et al.
presented a detailed derivation of the Lunardini solution and found that the accuracy of Lunardini solution was directly proportional to the Stefan number [29]. Vitel et al.  constructed a hydrothermal numerical model consistent  with thermodynamics to simulate the artificial ground  freezing of saturated nondeformable porous media under  seepage conditions. is numerical model has been well verified by the results of three-dimensional ground freezing experiments under high seepage velocity conditions [30][31][32]. Marwan et al. used the "ant colony algorithm" to optimize the layout of freezing pipes under the action of low-flow-rate groundwater, thereby shortening the closure time of the frozen wall and making the strength of the frozen wall more uniform. is method led to new ideas for the optimal design of the layout of freezing pipes [33]. Scheidegger et al. proposed a combination of one-dimensional heat conduction model and two-dimensional thermal convection model including freeze-thaw action to simulate permafrost development in complex geological environments and made it possible to evaluate the influence of convective heat flow and load on the state of glaciers [34]. Based on the heat transfer and seepage theory of the finite element method in porous media, Hu et al. built a fully coupled numerical model to simulate changes in temperature and groundwater flow fields; the ability of the model to solve the phase transition problem was verified by the classical analysis solution of freezing temperature field, and the factors that influence the frozen wall crossing time are also analyzed [35]. Grenier et al. verified a series of the numerical codes being developed for coupled nonlinear freeze-thaw systems, and this verification effort provided a basis for addressing more general challenges related to numerical simulation of cold-region hydrogeological processes [36]. Huang et al. developed a hydrothermal coupling model to simulate the effect of water flow on the freezing process by considering the water/ice phase transition and combined the model with the Nelder-Mead simplex method based on the COMSOL Multiphysics platform to optimize the positions of the freezing pipes around a circular tunnel [13,37,38]. Alzoubi et al. constructed a numerical model of the hydrothermal coupling using the enthalpy-porosity method and verified the rationality of the model through experiments. Based on the model, the factors affecting the development of the freezing temperature field were studied, and the influence of the seepage field on the freezing temperature field was analyzed in depth by introducing the concept of "Heatlines" [39,40].
After the artificial formation freezing enters the stable freezing stage, the temperature drop rate of the freezing temperature field becomes very slow; therefore, the transient temperature field at this time can be approximated by the steady-state temperature field [24][25][26][27][28]. When a seepage field exists in the stratum, after entering the stable freezing phase, the convective heat transfer between the water and the frozen soil on the freezing front and the heat transfer of the frozen pipe cancel each other out under the condition that the seepage velocity and the water temperature are constant. e temperature field at this time is close to the steady-state temperature field; therefore, it is also possible to use a steady-state temperature field instead of a transient temperature field for analysis.
Considerable progress has occurred in the analytical calculation of the steady-state temperature field of artificial freezing in nonseepage strata [24][25][26][27][28], but an analytical solution of the steady-state freezing temperature field under the action of a seepage field is lacking. e study of the freezing temperature field formed by single pipe is the basis of artificial frozen wall research and has important research value. In this paper, the analytical expression of the steadystate temperature field of single freezing pipe under the action of seepage field was derived, and the accuracy of the analytical solution was verified by physical model test. Based on the formulas provided in this paper, the temperature field formed by single freezing pipe under different seepage velocity was calculated.

Basic Assumption.
Under the action of the seepage field, the water flow changes the heat distribution of the entire frozen area through the convective heat transfer, and the freezing range of the downstream area is obviously larger than the upstream area; therefore, the shape of the freezing front is no longer a circle, but a shape close to the "cardioid" [12,15].
To facilitate the description of the distribution law of the freezing temperature field, a rectangular coordinate system was constructed as shown in Figure 1. e center of the freezing pipe was taken as the coordinate origin; the axis passing through the origin and parallel to the water flow direction was taken as the x-axis, and the axis passing through the origin and perpendicular to the water flow direction was taken as the y-axis. e intersections I(ξ V 1 , 0) and III(ξ V 3 , 0) of the freezing front and the x-axis represent the extended range of the upstream and downstream of the freezing front under the action of the seepage field, respectively; intersection point II(ξ V 2 , 0) represents the extended range of the side of the freezing front. Point M(x 0 , y 0 ) is any point within the freezing front of the right side of the x-axis (x > 0). M′(x 1 , y 1 ) is the intersection of the freezing front and straight line passing through the origin (0, 0) and the point M(x 0 , y 0 ). N(x 2 , y 2 ) is any point within the freezing front of the left side of the x-axis (x < 0) and N′(x 3 , y 3 ) is the intersection of the freezing front and straight line passing through the origin (0, 0) and the point N(x 2 , y 2 ).
According to the characteristics of the steady-state freezing temperature field of a single pipe under the seepage field, before performing analytical solution derivation, the following hypotheses are proposed: (1) e flow rate, direction, and temperature of the seepage field remain unchanged. (2) e porous medium is uniform (particle size, porosity, and thermal conductivity are equal at all locations). (3) e temperature of the freezing pipe wall is constant.
(4) After entering the stable freezing phase, the freezing temperature field becomes a steady-state temperature field. (5) As shown in Figure 1, taking the y-axis as the dividing line, the shape of the freezing fronts on both sides of the dividing line is equivalent to half an ellipse. And the shape of the freezing front can be expressed by the following equation:

Derivation of Analytical Solution.
According to the first law of thermodynamics, the equation for two-dimensional heat conduction is expressed as [41] zq where q x and q y are the heat flux densities along the x and y axes, respectively. According to Fourier's law, q x and q y can be expressed as [41] where k and T are the thermal conductivity and temperature of the soil, respectively. It is assumed that the soil is isotropic in thermophysics, so k was the same in all directions. And the following equation can be obtained [41]: e polar form of equation (5) is [41] In the steady-state freezing temperature field formed by a single freezing pipe, the temperature of the freezing pipe and the freezing front are stabilized at T f and T 0 , respectively. erefore, the boundary conditions are as follows: where T 0 is the freezing temperature of the soil, T f is the temperature of the freezing pipe, r 0 is the radius of the freezing pipe, and ξ is the distance from the freezing front to the center of the freezing pipe. Under the action of flowing water, the shape of the freezing front is no longer a circle, and ξ changes with position (x, y) and seepage velocity v. e freezing front is considered to consist of countless points. For any point where the distance from the freezing front to the freezing pipe is ξ, the point can also be regarded as a point on the circle with the freezing pipe as the center and the radius of ξ, as shown in Figure 1. According to the calculation principle of steady-state freezing temperature field of a single pipe [17,26,28] and combining equations (6) and (7), the following equation is obtained: When x ≥ 0, the straight line passing through the point M(x 0 , y 0 ) and the origin (0, 0) can be expressed as e expressions of the intersection M′(x 1 , y 1 ) of the straight line and the freezing front can be obtained by combining equations (2) and (9): Substitution of equation (11) into equation (8) yields e temperature expression at any point in the frozen region (x ≥ 0) is obtained by simplifying the following equation: When x < 0, the straight line passing through the point N(x 2 , y 2 ) and the origin (0, 0) can be expressed as e expressions of the intersection N′(x 3 , y 3 ) of the straight line and the freezing front can be obtained by combining equations (2) and (14): e distance from N′(x 3 , y 3 ) to the coordinate origin (0, 0) is Substituting equation (16) into equation (8) and considering the arbitrariness of N(x 2 , y 2 ), the temperature expression at any point in the frozen area (x < 0) is
Based on equations (13) and (17), the formula for the isotherm of any temperature T can be derived as At the wall of the freezing pipe, T � T f , equation (20) is simplified as 1 At the freezing front T � T 0 , equation (20) is simplified as 1 e expression for the extension radius ξ i of the isotherm of any temperature (T) in equation (20) is e frozen area surrounded by any isotherm can be expressed as e frozen area in any temperature range (T, T − Δ) is

Model Test Design.
To verify the rationality of the analytical solution, a large-scale physical model test system was designed. e whole test system consists of a porous medium test zone, a freezing system, a seepage field simulation system, and a data acquisition system, as shown in Figure 2. e test plane was set at 500 mm intervals in the box, the test plane consisted of 7 axes labeled A-G, and 13 measurement points were arranged on each axis. e spacing between the measurement points was 50 mm and that between adjacent axes was 200 mm. e arrangement of measurement points is shown in Figure 3. In the test, only the freezing pipe in the middle position is in the working state.
In this study, a total of 4 sets of tests were conducted, the control parameters of the seepage field in the test are shown in Table 1.
e wall temperature of the inlet pipe was taken as the main reference, and its temperature was controlled at −30°C. To reach this temperature and account for the cooling loss in the pipeline, the alcohol temperature and flow rate of the refrigeration unit were set to −32°C and 2.5 m 3 /h, respectively.

Analysis of Steady-State Temperature Field of Model Test.
e spatial distribution laws of the temperature of each measuring point (D1−13) on the D axis were plotted as graphs, as shown in Figure 4. Figure 4 shows that when no seepage occurred in the saturated porous medium, the temperatures at the measurement points located at the symmetric positions on the D axis were basically the same at the same time, indicating that the positions of the measurement points were relatively accurate and the sand layer was relatively uniform. e spacing between adjacent two curves in the late stage of freezing gradually decreased, indicating that the temperature drop of the freezing temperature field became slower after entering the relatively stable stage in the late stage of freezing. e curve with the freezing time of 24 h and the curve of 28 h basically coincided, so it can be considered that the freezing temperature field entered the stable freezing stage when the freezing time was 28 h.
When v � 3 m/day, the temperatures at the measurement points located in the downstream area were lower than those at the measurement points in the upstream area, indicating that the convective heat transfer of the water flow changed Advances in Civil Engineering the spatial distribution of the heat of the frozen area. e water was cooled as it flowed through the upstream region, the temperature difference between the water flow and the porous medium decreased, and the convective heat transfer weakened when flowing in the downstream region. us, the temperature drop rate at the measurement point in the downstream region was higher in the early stage of freezing.
As the freezing time increased, the downstream area gradually entered the stable freezing stage and the frozen range no longer changed, whereas the temperatures at the measurement points in the upstream area continued to decrease. erefore, the temperature difference between the upstream and downstream regions decreased during the late stage of freezing. e parts of the curves below 0°C with freezing times of 36 h and 40 h basically coincided, so it can be considered that the freezing temperature field entered the stable freezing stage when the freezing time was 40 h with the seepage velocity being 3 m/d.
When v � 6 m/day, the temperature of D1-D4 measurement points in the upstream area was maintained at about 15°C, indicating that under this flow rate condition, the freezing pipe could not affect the area beyond 150 mm upstream. e temperatures at the measurement points in the downstream region gradually decreased as the distance from the freezing pipe decreased; the shorter the distance from the freezing pipe, the larger the temperature gradient of the adjacent measurement points. e upstream part of the temperature distribution curves at different times was almost coincident, whereas the curves of the downstream part gradually decreased. e two curves with freezing times of 20 h and 24 h were coincident, indicating that the freezing temperature field entered a stable freezing stage when the freezing time was 24 h with a seepage velocity of 6 m/day.
When v � 9 m/day, the temperature at D1-D5 measurement points in the upstream area was about the same during the initial stage of freezing, and only the temperature at D5 slightly decreased during the late freezing stage. is showed that under this flow rate condition, the influence of the freezing pipe on the area outside the 100 mm range was weak. e temperatures at the measurement points in the downstream area gradually decreased as the distance from the freezing pipe decreased; the shorter the distance from the freezing pipe, the larger the temperature gradient of the adjacent measurement points. e three temperature curves for 16, 20, and 24 h coincided, indicating that the freezing temperature field entered a stable freezing stage when the freezing time was 20 h with a seepage velocity of 9 m/day. By comparing the temperature curves under several flow conditions, it can be found that in the initial freezing period, the larger the seepage velocity, the lower the temperature in the downstream region, and in the late freezing period, the  Advances in Civil Engineering smaller the seepage velocity, the lower the temperature in the downstream region. e reasons for this phenomenon are as follows: during the construction of the artificial ground freezing method, when there is a seepage field in the local layer, due to the convective heat transfer effect, the water flow will release heat when it passes through the freezing front, and the temperature of the water flow will decrease; when low-temperature water flows to downstream area where the temperature is higher, it absorbs heat and the temperature of the downstream formation decreases. In the initial stage of freezing, the difference in freezing area is small. e greater the seepage velocity, the more intense the convective heat transfer effect and the lower the temperature in the downstream area. As the freezing time increases, the difference between the frozen area at high seepage velocity and the frozen area at low seepage velocity increases. When the seepage velocity in the ground layer is small, the contact surface between the water flow and the freezing front is larger, which results in lower temperature of the water flowing to the downstream area than that with high seepage velocity. erefore, the lower the seepage velocity in the formation, the lower the temperature in the downstream area in the late stage of freezing.
According to the temperature test results, the extended radius ξ i of the freezing front of the steady-state temperature field in the three directions is shown in  law of the extended radius ξ i with the flow rate is shown in Figure 5, where t is the time required for the freezing temperature field to enter the stable freezing section. ξ 1 , ξ 2 , and ξ 3 are the extended radii of the freezing front facing the upstream, both sides, and the downstream direction of the freezing pipe, respectively. Table 2 and Figure 5 show that the extended radius ξ i of the freezing front of the steady-state freezing temperature field had three significant phases with the change in the seepage velocity v. When the seepage velocity v was lower than 3 m/day, ξ i changed little with the seepage velocity, and the differences between the extended radii ξ 1 , ξ 2 , and ξ 3 in each direction were small. When v increased from 3 to 6 m/ day, ξ 1 , ξ 2 , and ξ 3 sharply reduced from 220, 220, and 240 mm to 65, 108, and 169 mm, and the reduction ratio was 70.45%, 50.90%, and 29.58%, respectively. e gap between ξ 1 , ξ 2 , and ξ 3 at this stage gradually increased as the flow rate increased. When the seepage velocity v was greater than 6 m/ day, the expansion radius ξ i of the frozen area in three directions was further reduced, but the overall variation was small. e main reasons for the above phenomena are as follows: the development law of the freezing temperature field under the action of a seepage field was the result of the interaction of heat conduction and convection heat transfer, whereas the steady-state freezing temperature field was a state in which the two effects reached equilibrium. When the seepage velocity was low, the convective heat transfer effect was weak; therefore, when v increased from 0 to 3 m/day, the change in the extended radius ξ i of the freezing front of the steady-state freezing temperature field was small. e convective heat transfer was enhanced with the increase in the seepage velocity, and the heat transfer between the soil layer and the freezing pipe far from the freezing pipe was offset by convective heat transfer. erefore, when v increased from 3 to 6 m/day, the frozen area drastically reduced. When the frozen area reduced to a position close to the freezing pipe, the heat conduction in the frozen area became very strong. At this time, the increase in the seepage velocity had less influence on the frozen area. erefore, when the seepage velocity was larger than 6 m/day, the change in ξ i with seepage velocity was small.

Rationality Verification of Analytical Solution.
e calculation results and the test results of the temperature on the axis parallel to the direction of water flow and perpendicular to the direction of water flow under different seepage speed conditions were compared, as shown in Figures 6(a) and 6(b), respectively.
By comparison, we found that the formula calculation results had a high degree of agreement with the test results. erefore, a mathematical description of the steady-state temperature field under the action of seepage field can be realized based on this formula.

Calculation and Analysis of Steady-State Temperature Field of Single Freezing Pipe under Action of Seepage Field
e test results only partially reflect the distribution law of the freezing temperature field due to the limited number of measuring points in the test. However, based on the analytical solution of the steady-state freezing temperature field under the seepage field that was derived from this study, the distribution law of the freezing temperature field was comprehensively described. Based on the formula derived from this study, the overall distribution law of the freezing temperature field under different seepage velocities was obtained, as shown in Figure 7. Figure 7 shows that the expansion of the radius of the frozen area in all directions was equal without a seepage field, and the temperature of the frozen area dropped sharply as the distance from the freezing pipe decreased. When the seepage velocity was 3 m/day, after entering the stable freezing phase, the frozen zone upstream of the freezing pipe and on both sides of the freezing pipe was slightly reduced compared to that without seepage. It showed that the seepage field with velocity of 3 m/d would reduce the temperature drop rate of the freezing temperature field formed by single pipe, but it did not have a significant effect on the temperature distribution of the steady-state temperature field compared to that without seepage. When the seepage velocity was 6 m/day, the convective heat transfer between the water flow and the frozen zone intensified, so the frozen zone area drastically reduced. Since the upstream    Advances in Civil Engineering 9 position was affected by the frontal impact of the water flow, the frozen area of upstream position was significantly smaller than that of the frozen areas downstream and on both sides of freezing pipe. When the seepage velocity was 9 m/day, the range of the frozen area further decreased, but the overall distribution law of the temperature field was close to that with a seepage velocity of 6 m/day. Using equations (23) and (24), the expansion of any one isotherm in the frozen temperature field and the area enclosed by any isotherm can be calculated, respectively. e calculation results are shown in Figure 8. Figure 8(a) shows that without a seepage field, the area surrounded by the 0°C isotherm was 17.97 × 10 4 mm 2 . As the temperature decreased, the area enclosed by the isotherm decreased sharply, and the area enclosed by the −5°C isotherm was 7.09 × 10 4 mm 2 , accounting for 39.45% of the total frozen area, whereas the area enclosed by the −25°C isotherm had an area of 0.06 × 10 4 mm 2 , which was only 0.31% of the total frozen area. When a seepage field existed in the sand layer, the change law of the area enclosed by the isotherm at different temperatures was the same as that with no seepage field; the area enclosed by the isotherm decreased sharply as the temperature decreased. However, due to the convective heat transfer between the seepage field and the temperature field, the frozen area formed by a single freezing pipe was smaller than that without a flow rate. When v was 3 m/day, the area enclosed by the 0°C isotherm was 15.77 × 10 4 mm 2 , which was 12.24% lower than that without a flow rate. However, when v was 6 and 9 m/day, the area enclosed by the 0°C isotherm was only 3.84 × 10 4 and 3.05 × 10 4 mm 2 , respectively, which are 78.63% and 83.03% lower than that without flow rate, respectively.
is indicated that the seepage with a velocity greater than 6 m/day had a significant effect on the extended range of the steady-state temperature field of the single freezing pipe. e frozen area between any temperature range was calculated using equation (23). e proportional relationship between the frozen area S (T,T−Δ) of the temperature  range (T, T − Δ) and the whole frozen area S T�0 under different flow rate conditions was calculated, as shown in Figure 9. e average temperature of the entire frozen wall was calculated using the weighted average algorithm as shown in Table 3. Figure 9 and Table 3 show that compared with the no flow rate, when the seepage velocities were 3, 6, and 9 m/day, the proportion of frozen zone in the temperature range of 0 to −4°C was reduced from 60.57% to 59.81%, 50.16%, and 48.48%, respectively; the proportion of frozen areas in the temperature range of −5 to −9, −10 to −14, −15 to −19, and −20 to −24°C increased from 24.14%, 9.62%, 3.83%, and 1.53% to 24.31%, 9.89%, 4.02%, and 1.64%; 25.51%, 13.25%, 6.81%, and 3.55%; and 25.60%, 13.75%, 7.39%, and 3.98%, respectively. e average temperature of the freezing temperature field was reduced from −5.78°C to −5.86, −7.31, and −7.50°C, respectively. is showed that as the seepage velocity increased, the frozen area formed by single pipe decreased; the proportion of the low-temperature zone increased; and the average temperature of the temperature field decreased.

Conclusions
In this paper, the shape of freezing front was simplified using a segmentation-equivalent method. e analytical expression of the steady-state temperature field of single freezing pipe under the action of seepage field was derived, and the accuracy of the analytical solution was verified using a physical model test. e verification results showed that the calculated results of the analytical solution were highly consistent with the experimental results, so the steady-state temperature field under the seepage field could be mathematically described based on the formulas derived in this paper. e results of the model test and calculation can be summarized as follows:   (1) e model test results showed that when the seepage velocity was 0, 3, 6, and 9 m/day, the time required for the temperature field formed by a single freezing pipe to enter the stable freezing phase was 28, 40, 24, and 20 h, respectively, and the frozen area of the steady-state freezing temperature field was 17.97 × 10 4 , 15.77 × 10 4 , 3.84 × 10 4 , and 3.05 × 10 4 mm 2 , respectively.
is indicated that seepage with a velocity greater than 6 m/day had a significant effect on the extended range of the steadystate temperature field of a single freezing pipe.
(2) e calculation results of the formulas showed that compared with no flow rate, when the seepage velocity was 3, 6, and 9 m/day, the proportion of frozen area below −5°C increased from 39.43% to 40.19%, 49.84%, and 51.52%, respectively, and the average temperature of the freezing temperature field was reduced from −5.78 to −5.86, −7.31, and −7.50°C, respectively. is showed that as the seepage velocity increased, the area of the single-pipe freezing zone decreased, the proportion of the low-temperature zone increased, and the average temperature of the temperature field decreased.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.