A method for determining the unfrozen water content in the three-dimensional pores of frozen soil under freeze-thaw cycles

Abstract

The present application is suitable for frozen soil engineering technical field, and provides a method for determining the unfrozen water content in three-dimensional pore of frozen soil under freeze-thaw cycle. The present application focuses on the thawing process of frozen soil under freeze-thaw cycle, combines the adsorption and capillary action of water, establishes the quantitative relationship between the unfrozen water content and temperature, medium size and pore size based on three-dimensional pore structure, divides the change of unfrozen water content into three states, and deduces the expression formula for each state according to the micro-geometric relationship; when the values of temperature, medium size and pore size and other variables are known, the unfrozen water content can be quickly and accurately determined by judging the state of unfrozen water and substituting into the corresponding formula, which effectively improves the calculation efficiency and accuracy. The present application not only provides a new idea for the characterization of unfrozen water content in frozen soil during freeze-thaw cycle, but also assists the groundwater geophysical work in frozen soil area, especially has important application value for detecting the groundwater content in frozen soil area by using magnetic resonance technology.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of frozen soil engineering, and particularly relates to a method for determining the unfrozen water content in three-dimensional pores of frozen soil under freeze-thaw cycles. BACKGROUND

[0002] As a temperature-sensitive geologic body, the physical properties of the active layer frozen soil on the plateau change with temperature. The freeze-thaw process of the frozen soil causes the conversion of water-ice phase, which causes the interaction among the pores, ice and water in the frozen soil medium, accompanied by a series of freeze-thaw cycle phenomena such as pore frost heaving, ice freezing and unfrozen water migration, resulting in the extremely unstable pore structure, unfrozen water content and ice occurrence state of the frozen soil. In this case, it is of great significance to accurately obtain the relationship between the unfrozen water content and the pore distribution and temperature change.

[0003] A method for obtaining the unfrozen water content in rocks under freeze-thaw cycles is disclosed in Chinese patent CN110806422A. The method prepares a rock column sample, saturates the rock column sample with formation water, obtains a T2 spectrum curve through nuclear magnetic resonance testing, converts the T2 spectrum curve into a nuclear magnetic pore throat distribution curve, and then substitutes it into an integral formula to obtain the relationship curves between the unfrozen water content and temperature at the time of melting and freezing, respectively. This method has the advantages of high precision, high reliability and easy access to the unfrozen water content. A method for calculating the unfrozen water content of frozen soil considering the initial water content is disclosed in Chinese patent CN115047032A. The method obtains a theoretical formula for the unfrozen water content in frozen soil considering the initial water content by the steps of soil sampling, parameter measurement and substitution of the calculation coefficient into the function, and uses the obtained formula to calculate the unfrozen water content in the soil freeze-thaw process. This method has the advantages of simple expression, convenient use and fast calculation of the unfrozen water content in the soil freezing and thawing process.

[0004] However, both of the above methods have certain limitations in actual application. For the method disclosed in patent CN110806422A, it needs to prepare a rock sample first, and then perform a series of conversion calculations after nuclear magnetic resonance testing. The overall process is relatively cumbersome, and the calculation process does not consider the geometric relationship between the pore and the unfrozen water system in detail. The method disclosed in patent CN115047032A focuses on the resistivity model, measures the change of the soil resistivity in the freeze-thaw process, and analyzes the relationship between the change of the resistivity and the change of the unfrozen water content and temperature. However, resistivity measurement is an indirect measurement mode, and the measurement result is easily disturbed by various external factors. At the same time, this method does not consider the changes in the form of unfrozen water caused by adsorption and capillary action, and the changes in the form of unfrozen water will directly cause the deviation of the calculation result of the unfrozen water content and affect the accuracy of the calculation.

[0005] In view of the above deficiencies of the prior art, the present application provides a method for determining the unfrozen water content in the three-dimensional pores of frozen soil under freeze-thaw cycles. SUMMARY

[0006] The purpose of the present application is to provide a method for determining the unfrozen water content in the three-dimensional pores of frozen soil under freeze-thaw cycles, aiming to solve the problems raised in the background art.

[0007] The purpose of the present application is achieved by the following technical solutions:

[0008] A method for determining the unfrozen water content in the three-dimensional pores of frozen soil under freeze-thaw cycles, comprising the following steps:

[0009] Sample acquisition and parameter determination: acquire a frozen soil sample and determine its basic physical parameters;

[0010] Temperature measurement: perform a freeze-thaw cycle process on the frozen soil sample and monitor its temperature change in real time;

[0011] State determination: calculate the numerical value of the capillary water meniscus curvature radius r and determine the unfrozen water state of the frozen soil sample at the current temperature; the states include:

[0012] State 1: when r=0, only adsorbed water exists in the system;

[0013] State 2: when 0

[0014] State 3: when r> the inscribed sphere radius of the pore model, the capillary water is connected;

[0015] Unfrozen water content calculation: according to the determined state, use the corresponding unfrozen water content calculation formula to calculate the unfrozen water content under the corresponding state.

[0016] Further, in the state determination step, the capillary water meniscus curvature radius r is calculated by the following formula:

[0017] ;

[0018] Wherein: is the capillary water contact angle; is the ice-water interface free energy coefficient; is the capillary water pressure.

[0019] Further, the capillary water pressure is calculated by the total pressure of water and the adsorbed water pressure , and the relationship is as follows:

[0020] ;

[0021] Wherein: adsorbed water pressure Molecular force Electrostatic force And structural force Together constitute the relationship as follows:

[0022] .

[0023] Further, in the state determination step, the inscribed sphere radius of the pore model is solved according to the geometric relationship as (R= (3V) / (4S) ) (1) ) , wherein The medium average radius of the sample.

[0024] Further, in the state 1, the unfrozen water content is equal to the adsorbed water content, wherein the adsorbed water volume is equal to the product of its surface area and adsorbed water thickness d The unfrozen water content is solved according to the geometric relationship as (Wu= (S*H) / (3V) ) (2) :

[0025] ;

[0026] Wherein, The density of water.

[0027] Further, in the state 2, the unfrozen water content is equal to the sum of the adsorbed water content and the capillary water content, and the unfrozen water content is solved according to the geometric relationship as (Wu= (S*H) / (3V) + (S*H) / (2R) ) (3) :

[0028] ;

[0029] Wherein, the upper limit of the integral is solved according to the geometric relationship as (H= (2R) / (2R-1) ) (4) , The distance from the vertex of the capillary water meniscus to the bottom edge, The cross-sectional area of the capillary water, The angle set in the geometric relationship.

[0030] Further, in the state 3, the unfrozen water volume is equal to the difference between the total pore volume and the ice volume, and the unfrozen water content is solved according to the geometric relationship as (Wu= (V- (4 / 3) *Pi*R^3) / (4 / 3) *Pi*R^3) (5) :

[0031] .

[0032] Compared with the prior art, the beneficial effects of the present application are:

[0033] The application focuses on the thawing process of frozen soil under freeze-thaw cycle conditions, fully considers the influence of water adsorption and capillary action on the form of unfrozen water, and based on the dense distribution structure of the medium, a technical solution is constructed based on three-dimensional pore structure: on the one hand, the quantitative relationship between the content of unfrozen water and temperature, medium size, pore size is established, the change of the content of unfrozen water in the thawing process is divided into three states, and the expression formula is derived according to the micro-geometric relationship of each state; on the other hand, when the temperature, medium size and pore size and other variable values are known, the state of unfrozen water can be determined and substituted into the corresponding formula to quickly and accurately determine the content of unfrozen water, which effectively improves the calculation efficiency and accuracy. The application not only provides a new idea for the characterization of the content of unfrozen water in frozen soil in the freeze-thaw cycle process, but also assists the groundwater geophysical work in the frozen soil area, especially has important application value for detecting the content of groundwater in the frozen soil area by using magnetic resonance technology. BRIEF DESCRIPTION OF DRAWINGS

[0034] Figure 1 The method flowchart of the application.

[0035] Figure 2 The three-dimensional medium tight distribution schematic diagram; wherein (a) is the left and right two-equal-angle axial side view; (b) is the front view.

[0036] Figure 3 The unfrozen water state transformation schematic diagram in the thawing process; wherein (a) is the initial thawing (no capillary water); (b) is the formation of capillary water, before intersection; (c) is the intersection of capillary water critical; (d) is after the intersection of capillary water; (e) is the end of thawing.

[0037] Figure 4 The unfrozen water state determination flowchart.

[0038] Figure 5 The solution diagram of the meniscus curvature radius r of capillary water. DETAILED DESCRIPTION

[0039] In order to have a clearer understanding of the technical features, objectives and beneficial effects of the application, the technical solutions of the application will be described in detail as follows, but it cannot be understood as limiting the scope of the application.

[0040] The specific implementation of the application will be described in detail in combination with specific embodiments.

[0041] This invention provides a method for determining the unfrozen water content in the three-dimensional pores of frozen soil under freeze-thaw cycles. In freeze-thaw cycle studies, the core variables of frozen soil samples include pore distribution, temperature, and unfrozen water content. Temperature is a variable that is easily measured directly, while pore distribution is relatively stable and unrelated to temperature changes. Unfrozen water content, however, is influenced by both pore distribution and temperature, and these two factors jointly determine the dynamic changes in unfrozen water content. Therefore, the core objective of this invention is to establish a quantitative relationship among pore distribution, temperature, and unfrozen water content. By using directly measurable temperature combined with known pore distribution parameters, the unfrozen water content can be accurately calculated. The specific method is as follows (flowchart shown). Figure 1 (as shown)

[0042] I. Sample Acquisition and Parameter Measurement;

[0043] Obtain frozen soil samples and determine their basic physical parameters (mean radius of the medium of the sample). Average pore radius );

[0044] II. Temperature Measurement;

[0045] The frozen soil sample was subjected to a sufficient freeze-thaw cycle (placed for a set duration in a temperature-controlled instrument). During this process, the internal medium distribution of the frozen soil tended to be uniform and densely distributed. Figure 2 As shown in (a) and (b), the temperature changes of the frozen soil samples were monitored in real time during the freeze-thaw cycle and subsequent thawing analysis, and the current temperature was recorded. data.

[0046] III. Determination of State and Calculation of Unfrozen Water Content;

[0047] Taking the melting process as an example, initially only adsorbed water exists in the system, and the unfrozen water content equals the adsorbed water content, which is defined as state 1. As the temperature rises, melting proceeds, the thickness of the adsorbed water changes, and capillary water appears. At this point, the unfrozen water content equals the sum of the adsorbed water content and the capillary water content, which is defined as state 2. When the temperature continues to rise, melting continues until the capillary water reaches contact. At this point, the unfrozen water volume equals the difference between the total pore volume and the ice volume, which is defined as state 3. The above state changes are as follows: Figure 3 As shown in (a)-(e), the methods for calculating the unfrozen water content differ for the three states.

[0048] Regarding the determination of the three states during melting mentioned above, this invention classifies the states based on the radius of curvature r of the capillary manifold surface: when r = 0, the system is in state 1; when 0 < r < the radius of the inscribed sphere of the pore model, the system is in state 2; when r > the radius of the inscribed sphere of the pore model, the system is in state 3. The above-mentioned unfrozen water state determination process is as follows:Figure 4 The inner sphere radius of the pore model can be solved according to the geometric relationship as follows: ) wherein is the medium average radius of the sample. The derivation of r is as follows:

[0049] Step 1: Derivation of the curvature radius r of the capillary water meniscus (including the calculation of secondary correlation parameters);

[0050] (1) The curvature radius r of the capillary water meniscus can be solved by the following formula:

[0051] Formula 1: ;

[0052] wherein is the capillary water contact angle, which approaches 0°; is the ice-water interface free energy coefficient; is the capillary water pressure.

[0053] (2) The capillary water pressure in Formula 1 can be solved by the following formula:

[0054] Formula 2: ;

[0055] Formula 3: ;

[0056] wherein is the total pressure of water; represents the adsorbed water pressure, which is composed of molecular force , electrostatic force and structural force .

[0057] (3) The molecular force , electrostatic force and structural force in Formula 3 can be solved by the following formula:

[0058] Formula 4: ;

[0059] Formula 5: ;

[0060] Formula 6: ;

[0061] wherein is the Hamaker constant (J); is the adsorbed water thickness; is the relative dielectric constant of liquid water; is the vacuum dielectric constant; is the potential difference of the interaction surface; For surface correlation coefficient; For characteristic thickness of hydration layer.

[0062] (4) The adsorbed water thickness described by formula 4, formula 5, formula 6 , can be solved by the following formula:

[0063] Formula 7: ;

[0064] Wherein: And The parameter related to the interfacial interaction potential energy; , The melting point temperature, The current temperature; The density of water; The latent heat of phase change; The medium average radius of the sample.

[0065] (5) The melting point temperature described by formula 7 , can be solved by the following formula:

[0066] Formula 8: ;

[0067] Wherein: 273.15K; The density of ice; The average pore radius.

[0068] (6) The total pressure of water described by formula 2 , can be solved by the following formula:

[0069] Formula 9: ;

[0070] Wherein: And The total pressure of ice and the density of ice respectively.

[0071] Step 2: The independent variable of the capillary water meniscus curvature radius r is summarized;

[0072] Through the hierarchical derivation of step 1, it is finally determined that the independent variable of r is only related to the current temperature , the average pore radius And the medium average radius of the sample R (see Figure 5 The solving logic is shown).

[0073] Step 3: Calculation of unfrozen water content;

[0074] (1) Calculate the unfrozen water content of state 1 When r = 0, the system is in state 1, and the unfrozen water content is equal to the adsorbed water content, wherein the adsorbed water volume is equal to the product of the surface area of the adsorbed water in contact with the medium and the adsorbed water thickness The unfrozen water content in state 1 is obtained according to a geometric relationship as follows:

[0075] Equation 10:

[0076] (2) Calculation of the unfrozen water content in state 2 When 0 < r < (r0 ) , the system is in state 2, and the unfrozen water content is equal to the sum of the adsorbed water content and the capillary water content, and the unfrozen water content is obtained according to a geometric relationship as follows:

[0077] Equation 11:

[0078] wherein the is obtained by a geometric relationship is the distance from the vertex of the capillary water meniscus to the bottom edge; is the capillary water cross-sectional area; is an angle set in the geometric relationship.

[0079] The distance from the vertex of the capillary water meniscus to the bottom edge and the capillary water cross-sectional area , according to a geometric relationship, can be obtained by the following equation:

[0080] Equation 12:

[0081] Equation 13:

[0082] wherein and (mutual, same as , a set parameter in the geometric relationship solving process, and does not have physical meaning) are reciprocal, and according to a geometric relationship, the following can be obtained:

[0083] Equation 14:

[0084] Equation 15:

[0085] (3) Calculation of the unfrozen water content in state 3 When r > (r0 ) , the system is in state 3, at this time, the ice is completely wrapped by water, and the unfrozen water volume is equal to the difference between the total pore volume and the ice volume, and the unfrozen water content is obtained according to a geometric relationship as follows: ​​​​​​​​​​

[0086] Formula 16: .

[0087] The above are only preferred embodiments of the present application, it should be noted that for those skilled in the art, without departing from the concept of the present application, can also make several variations and improvements, these should also be considered as the protection scope of the present application, these will not affect the effect and practicality of the present application.

Claims

1. A method for determining the unfrozen water content in the three-dimensional pore of frozen soil under freeze-thaw cycles, characterized by, The method comprises the following steps: Sample acquisition and parameter determination: acquire frozen soil samples and determine their basic physical parameters; Temperature measurement: perform freeze-thaw cycle process on the frozen soil samples and monitor their temperature changes in real time; State determination: calculate the numerical value of the capillary water meniscus curvature radius r and determine the unfrozen water state of the frozen soil samples at the current temperature; the states include: State 1: when r = 0, only adsorbed water exists in the system; State 2: when 0 < r < the inscribed sphere radius of the pore model, adsorbed water and capillary water coexist in the system; State 3: when r > the inscribed sphere radius of the pore model, capillary water is in contact; Unfrozen water content calculation: according to the determined state, the corresponding unfrozen water content calculation formula is used to calculate the unfrozen water content under the corresponding state; In the state determination step, the capillary water meniscus curvature radius r is obtained by calculation through the following formula: ； wherein: is the capillary water contact angle; is the ice water interfacial free energy coefficient; is the capillary water pressure; In said state 1, the unfrozen water content is equal to the adsorbed water content, where the adsorbed water volume is equal to the product of its surface area and the adsorbed water thickness d, which is solved from the geometric relationship is: ； wherein, is the density of water, is the medium average radius of the sample; In the state 2, the unfrozen water content is equal to the sum of the adsorbed water content and the capillary water content, and the unfrozen water content is solved according to the geometric relationship is: ； Wherein, according to the geometric relationship, the integral upper limit is solved , The distance from the vertex of the capillary water bend to the bottom edge, The capillary water cross-sectional area, The angle set in the solution of the geometric relationship; In the state 3, the unfrozen water volume is equal to the difference between the total pore volume and the ice volume, and the unfrozen water content is solved according to the geometric relationship is: 。 2. The method of determining unfrozen water content in frozen soil three-dimensional pores under freeze-thaw cycles according to claim 1, characterized in that, The capillary water pressure By the total pressure of water With the adsorbed water pressure The relationship is calculated as follows: ； Where: adsorbed water pressure Molecular forces Electrostatic forces And structural forces Together comprise the following relationship: 。 3. The method of determining unfrozen water content in frozen soil three-dimensional pores under freeze-thaw cycles according to claim 1, characterized in that, In the state determination step, the inscribed sphere radius of the pore model is solved according to a geometric relationship as ) wherein is the medium average radius of the sample.

Images