Semiempirical Correlation betweenP-WaveVelocity andThermal Conductivity of Frozen Silty Clay Soil

Xi’an Research Institute Co. Ltd., China Coal Technology and Engineering Group Corp., Xi’an, 710077, Shaanxi, China Shaanxi Key Laboratory of Coal Mine Water Hazard Prevention and Control Technology, Xi’an, 710077, Shaanxi, China School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, China School of Geology and Environment, Xi’an University of Science and Technology, Xi’an 710054, China


Introduction
Seasonal frozen regions are widely distributed in the civilized world. e structure, physical and mechanical properties, and thermal physical properties of the soil undergo substantial changes due to the freeze-thaw cycle caused by temperature changes throughout the year [1]. In addition, the thawing temperature decreases as the higher pressure is applied on the ice or soil. e thawing of the seasonal frozen soil in the cold regions occurs because of the substructures and superstructures, which are contacting the frozen soil. Subsequently, the microfissures can be generated due to the volume shrinkage of frozen soil. e bearing capacity of foundation decreases, and the consolidation settlement of the superstructures constructed in the cold regions occurs accordingly. When the engineering structure is built on the seasonally frozen soil, the interaction between engineering structure and frozen soil is therefore studied deeply. Especially, the temperature field, frost heave, and consolidation and settlement of frozen soil need to be studied in depth [2,3]. Soil thermal conductivity is an important heat transfer property which shows the ability of soil to conduct heat [4,5].
ere are many factors that affect the thermal conductivity such as mineral composition, moisture content, porosity, and temperature. In addition, thermal conductivity is an important parameter in the calculation of the temperature field [6,7]. e analysis and evaluation of soil thermal conductivity during the freeze-thaw cycle have become a focus in engineering and theoretical research. e estimated temperature distribution and freezing degree are essential to the structural health of buildings [8].
In the last few decades, some experimental and theoretical approaches have been proposed in the literature to determine the soil thermal conductivity [9][10][11][12][13]. For example, the weighted average thermal conductivity of each component in the soil matrix was proposed in the De Vries model [14] and a correlation between thermal conductivity and degree of saturation which incorporates a variable kappa that accounts for the soil type effect was deduced by Cote and Konrad [15]. Lu et al. [16] carried out a series of tests with a thermo-time domain reflectometry probe on twelve types of natural soils which ranged from sand and silt to clay loam and established a thermal conductivity model by examining a new relationship between thermal conductivity and saturation. Ultrasonic technique has been used to characterize the constituents and structure of soils [17]. For instance, researchers used the ultrasonic technique to determine physical and mechanical properties of frozen soils [18]. e main factors affecting soil P-wave velocity are soil mineralogy, saturation, porosity, cementation, and temperature. e fact that both thermal conductivity and P-wave velocity are affected by the same parameters may help to determine their correlation [13,19]. In addition, the elastic waves of soil can be determined easily and quickly, and therefore determining the soil thermal conductivity from soil elastic wave is possible and timely. However, very few studies have attempted to generalize the correlation between thermal conductivity and elastic wave [20,21].
In this study, to thoroughly examine the variations in thermal conductivity and P-wave velocity of silty clay soil, experiments on soil samples with different water contents at different temperatures are conducted and an association model for these two physical parameters is established. e applicability of the proposed theoretical equations for thermal conductivity and P-wave velocity of silty clay soil is also validated.

A Theoretical Model for P-Wave Velocity and Thermal Conductivity of Frozen Soil
Few research studies have attempted to develop theoretical models to determine the thermal conductivity and P-wave velocity of frozen soils. Moreover, generalizations of the relationships between thermal conductivity and P-wave velocity of frozen soils have yet to be derived. e thermal conductivity of frozen soil depends on the thermal conductivity of the single component (solid particles, water, ice, and air) and upon their volume fractions. According to the contributions of four-phase substances, the series model simplifies the structure of soils [22]. If the direction of the heat flow is parallel to the direction of each component (Figure 1), the thermal conductivity of frozen soils can be determined by the weighted arithmetic mean method as follows: where λ and S are the thermal conductivity and cross section area of frozen soil, respectively, and the subscripts s, w, i, and a represent solid particles, pore water, ice, and air, respectively (S s + S w+ S i+ S a � 1). It is suggested that the change in volume of water when it is frozen can be neglected. e unfrozen water content of the frozen soil at a certain temperature is w u . en, where n, S r , and w are the porosity, saturation, and gravimetric water content of the frozen soil. Equation (1) can be described as Since n(1− S r )<1 and λ a ≪ λ s , λ a n(1 − S r ) can be neglected. Equation (3) can be written as e porosity and saturation are constants for a soil with a given water content, and λ s , λ w , and λ i are constants. us, equation (4) can be described as where a � λ s (1 − n) + nS r λ i and b � nS r (λ w − λ i ) are constants. Previous studies have demonstrated that the wave velocity can be used to determine the unfrozen water content [23,24]. Timur [25] proposed a model to estimate the wave velocity in a frozen soil, which is written as where v, v s , v w , and v i are the effective, solid, water, and ice wave velocities and ϕ w , ϕ i , and ϕ s are the proportion of the phases (ϕ w + ϕ i + ϕ s � 1). en,

Heat flux
Solid particles Water Ice Gas Equation (6) can be described as Since v s , v w , and v i are constants, equation (7) can be described as constants for a soil with a given water content. erefore, the unfrozen water content can be described as Incorporating equation (10) into equation (5), the correlation between thermal conductivity (λ) and wave velocity (v) is expressed as where to the soil porosity and soil saturation and are constants for a soil with a given water content. It can be therefore observed that the thermal conductivity of frozen soil is inversely proportional to wave velocity.

Preparation of Soil
Samples. e silty clay soil used for the determination of thermal conductivity and P-wave velocity was collected from typical eluvial deposits in farmland in Xuzhou city, Jiangsu Province, China. e basic physical properties and the grain-size distribution of the experimental silty clay soil are shown in Table 1 and Figure 2. e main composition is quartz and kaolinite. First, the soil samples with different gravimetric water contents of 10%, 12%, 15%, 18%, and 20% were put in a stainless steel cylindrical container (Φ5 × 10 cm). e detailed production process of soil samples is given by Geng et al. [26]. After compaction, the samples are cured naturally for 24 h and then demoulded. en, the soil samples were subjected to eleven target temperatures of 10°C, 5°C, 2°C, 0°C, − 2°C, − 4°C, − 6°C, − 8°C, − 10°C, − 15°C, and − 20°C in a BPH-B incubator (Shanghai Yiheng Scientific Instruments Co., Ltd, China) for 24 h to ensure uniform cooling after reaching the target temperatures. e thermal conductivity and P-wave velocity were then measured.

Test Instruments and Method.
e thermal conductivity of soil samples was measured through the test hole on the incubator using a DRE-2C thermal conductivity tester which uses the transient plane source method. e details on the DRE-2C thermal conductivity tester are given in the study by Sun and Lü [13]. en, the P-wave velocity was measured using a RSM-SY5 acoustic wave detector (see Figure 3). e sampling interval was 0.5 μs. Vaseline petroleum jelly was used as a medium to maximize acoustic coupling and smeared onto the transducer and the surface of the sample.

Test Results and Analysis
4.1. Analysis of P-Wave Velocity. When a soil sample is subjected to a low temperature, the phase transition of water and the development of cracks will lead to changes in the number of pores and pore structure, which will change their P-wave velocities. e freezing temperature of soils is affected by the soil texture, structure, moisture, and salinity. Ultrasonic p-wave velocity has been used in previous studies to determine the unfrozen water content and freezing temperature of frozen soil [23,24]. e temperature at which the P-wave velocity rapidly increases can be considered as the freezing temperature. e variations in the P-wave velocity of the samples with different water contents are shown in Figure 4. It can be observed that the freezing temperature of the silty clay samples is about − 2°C. At a temperature range of 10°C to − 2°C or Stage I, the P-wave velocity is almost constant. e P-wave velocity of the samples with a water content of 10% is about 750 m/s while the P-wave velocity of those with a water content of 15% and 18% is close to 1000 m/s. When the temperature is low and ranges from − 2°C to − 20°C, the P-wave velocity varies and can be described as a rapid increase followed by a slower increase and then gradually becomes constant. e specific rates of increase are shown in Table 2. At a temperature range of − 2°C to − 6°C or Stage II, it can be seen that samples with a higher water content show a greater increase in their P-wave velocity. However, samples with a low water content show a greater increase in their P-wave velocity from − 6°C to − 20°C or at Stages III and IV. Figure 5 shows the relationship between thermal conductivity and temperature of the silty clay samples with different water contents. It can be seen that the variations in thermal conductivity are divided into four stages. e temperature ranges from 10°C to − 2°C in Stage I, in which the thermal conductivity of samples with a water content of 10%, 12%, 15%, 18%, and 20% is basically constant at 1.37, 1.45, 1.69, 1.68, and 1.66 W·m − 1 K − 1 , respectively. At a temperature range of − 2°C to − 6°C or at Stage II, the thermal conductivity substantially increases with reduced temperature. However, the thermal conductivity quickly decreases as the temperature declines for the soil samples with a high water content (15∼20%) from − 4°C to − 6°C. At a temperature range of − 6°C to − 10°C or Stage III, the thermal conductivity of all of the samples shows a slow trend of decline with decreased temperature.

Shock and Vibration
When the temperature is lower than − 10°C or at Stage IV, the thermal conductivity remains constant on the whole at 1.84, 1.95, 2.1, 2.15, and 2.25 W·m − 1 K − 1 with a water content of 10%, 12%, 15%, 18%, and 20%, respectively. It is obvious that the thermal conductivity at − 20°C is about 0.5 W·m − 1 K − 1 higher than that at 10°C.

Discussion
e unfrozen water content in frozen soil has a great influence on the physical and mechanical properties of soil. Previous research studies have shown that the unfrozen water content is rapidly reduced from − 2°C to − 5°C and then is reduced slowly [27][28][29]. When the temperature is below − 20°C, the unfrozen water content is approximately the same even though the initial water content is different. In our experiment, the unfrozen water content greatly varies between − 2°C and − 6°C. e variations in the thermal conductivity and P-wave velocity can be examined as four phases: pretransformation, intense transformation, transition, and permanently frozen phases.
(1) Pretransformation (temperature range of 10°C to − 2°C): e degree that the samples are frozen is very low, and there is no substantial change in their internal structure. e thermal conductivity decreases slowly with decrease in temperature due to the deceleration of molecular heat movement. e P-wave velocity is basically constant.     clay therefore changes from a particle assembly of individual grains to a solid composite that consists of particles and solid water [30]. Meanwhile, since the thermal conductivity and P-wave velocity of water is lower than ice, the thermal conductivity and P-wave velocity of the samples increase substantially. In addition, the volume expansion of water due to the phase transition separates the solid particles [24]. Cracks therefore propagate on the surface of the soil samples due to frost heave. At a temperature of − 6°C, frost heave cracks are greatly generated for the samples with a higher water content as shown in Figure 6(a). Consequently, the thermal conductivity slightly reduces with decreased temperature. e thermal conductivity for the samples with a water content of 10% and 12% increases as there are few cracks. Soil samples with high water content have more frost heave cracks, so the volume fraction of air in samples is higher. erefore, the P-wave velocity for the samples with a water content of 18% and 20% is lower than that for the samples with a water content of 15%.
(3) Transition (temperature range of − 6°C to − 10°C): A small amount of water continues to freeze in soil samples. erefore, the thermal conductivity is slowly decreased with the development and extension of frost heave cracks as shown in Figure 6(b). e P-wave velocity is gradually increased, and the rate of increase in the P-wave velocity is less than that in the intense transformation phase. (4) Permanently frozen (temperature range of − 10°C to − 20°C): In general, soil samples have a low unfrozen water content, and thus the frost heave cracks are rarely generated. e thermal conductivity remains constant, and the P-wave velocity is slightly increased. Figure 7 shows the variations in the generalized correlation between thermal conductivity and P-wave  Shock and Vibration velocity of frozen silty clay soil from − 4°C to − 20°C. e thermal conductivity has an inversely proportional relationship with P-wave velocity. e statistical parameters A and B for the various fits are listed in Table 3. e correlation coefficients of the equations fitted by the experimental data are high. e results of the experiments are consistent with the inversely proportional function of the theoretical correlation for the thermal conductivity and the P-wave velocity. e parameters A and B are mainly related to the water content of samples. It can be observed from Table 3 that the parameter A is first decreased with water content in general and then is increased when the water content exceeds 15%. e fitting parameter B is first increased with water content and then is decreased when the water content exceeds 15%. Determining soil thermal conductivity through correlation to its P-wave velocity seems to be a simple and effective solution.
e generalized correlation between thermal conductivity and P-wave velocity for the frozen silty clay soil is studied in the paper. Future research needs to be conducted on the different types of soil, and more testing data need to be analyzed.

Conclusions
In summary, this paper has examined the variations in thermal conductivity and P-wave velocity for samples with different water contents under low temperatures based on the formation of frozen soil and analyzed the correlation between the various parameters.
e main outcomes are summarized as follows: (1) e variations in thermal conductivity and P-wave velocity with temperature can be summarized as four stages: pretransformation phase, intense transformation phase, transition phase, and permanently frozen phase. e thermal conductivity and P-wave velocity all rapidly increase in intense transformation phase due to the formation of ice which bonds the solid particles. e thermal conductivity of soil samples with a water content of 15%, 18%, and 20% is slightly reduced at − 6°C due to the rapid growth of the frost heave cracks.
(2) e correlation between thermal conductivity and P-wave velocity of frozen silty clay soil is well described by the equation λ � A + B/v. e experimental results verify the validity of the theoretical equation for the thermal conductivity and the P-wave velocity. e findings provide the possibility for estimating the thermal conductivity of frozen soil from the P-wave velocity in practice.

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