CN115310026B - Bentonite expansive force prediction method and system considering ionic hydration energy - Google Patents

Bentonite expansive force prediction method and system considering ionic hydration energy Download PDF

Info

Publication number
CN115310026B
CN115310026B CN202211243933.4A CN202211243933A CN115310026B CN 115310026 B CN115310026 B CN 115310026B CN 202211243933 A CN202211243933 A CN 202211243933A CN 115310026 B CN115310026 B CN 115310026B
Authority
CN
China
Prior art keywords
bentonite
potential
dimensionless
ratio
expansive force
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211243933.4A
Other languages
Chinese (zh)
Other versions
CN115310026A (en
Inventor
孙海泉
王立忠
洪义
国振
李玲玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hainan Institute of Zhejiang University
Original Assignee
Hainan Institute of Zhejiang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hainan Institute of Zhejiang University filed Critical Hainan Institute of Zhejiang University
Priority to CN202211243933.4A priority Critical patent/CN115310026B/en
Publication of CN115310026A publication Critical patent/CN115310026A/en
Application granted granted Critical
Publication of CN115310026B publication Critical patent/CN115310026B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Processing Of Solid Wastes (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention relates to the fields of geotechnical engineering and soil mechanics, in particular to a bentonite expansive force prediction method and a system considering ionic hydration energy, wherein the bentonite expansive force prediction method considering ionic hydration energy comprises the following steps: acquiring the porosity ratio of the bentonite to be detected; obtaining the compensation dimensionless potential of the dimensionless potential and the ion hydration energy of the bentonite to be detected by utilizing the pore ratio; constructing a bentonite expansive force prediction model according to the dimensionless potential and the compensation dimensionless potential; and predicting the expansive force of the bentonite to be detected through the bentonite expansive force prediction model. According to the invention, the bentonite expansive force can be accurately predicted by utilizing the porosity ratio of the bentonite through combining the physical and chemical theory of the soil and the contribution of charge potential between bentonite layers to the expansive force of the bentonite. The invention provides a material selection method with high reliability and high precision for selecting the barrier buffer material for nuclide migration in a high-radioactive nuclear waste disposal library.

Description

Bentonite expansive force prediction method and system considering ionic hydration energy
Technical Field
The invention relates to the fields of geotechnical engineering and soil mechanics, in particular to a bentonite expansive force prediction method and system considering ionic hydration energy.
Background
The bentonite has the characteristics of high expansion force, low permeability, high adsorbability and the like, so the bentonite is preferentially applied to a high-radioactive nuclear waste disposal warehouse as a barrier buffer material for nuclide migration in various countries in the world. The swelling property of bentonite is due to its special soil structure. In order to maximize the swelling and sealing properties of bentonite, bentonite powder having a certain water content is previously compacted into a bentonite compact having a certain initial dry density. These compacted bentonite samples can be placed between the nuclear waste disposal tank and the underground surrounding rock as a barrier material, and when groundwater is immersed in the bentonite, the bentonite swells in the presence of water, generating expansive force, preventing the immersion of groundwater, and if the nuclear waste disposal tank leaks, the bentonite also prevents nuclides from migrating to the outside. The swelling force of bentonite is related to various factors such as the compaction degree (namely dry density), the water content, the mineral components and the like of the bentonite. If the expansive force of the selected compacted bentonite sample is too large, damage can be caused to surrounding rocks, and further infiltration of underground water is induced; if the expansion force of the sample is too small, the sample cannot play the role of a barrier buffer nuclear waste disposal tank; therefore, bentonite with expansive force within a certain range needs to be selected as a buffer material, but the prediction of the expansive force of the bentonite is based on an empirical method, so that the applicability is low, the prediction result is inaccurate and is not uniform; meanwhile, the traditional prediction of the swelling force of the bentonite based on the theory of physical and chemical aspects of soil only considers the contribution of charge potential among clay sheets to the swelling force of the bentonite, does not consider the contribution of hydration energy of ions in the bentonite after hydration to the swelling force of the bentonite, and has great deviation between the prediction result and the actual result. At present, a bentonite expansive force prediction method based on a pore ratio, which considers the soil body structure characteristics of bentonite and fully considers the influence of energy after ion hydration on the expansive force of the bentonite, has no solution.
Disclosure of Invention
In view of the defects of the prior art, in a first aspect, the invention provides a bentonite expansive force prediction method considering ionic hydration energy, which comprises the following steps: acquiring the porosity ratio of the bentonite to be detected; obtaining the compensation dimensionless potential of the dimensionless potential and the ion hydration energy of the bentonite to be detected by utilizing the pore ratio; constructing a bentonite expansive force prediction model according to the dimensionless potential and the compensation dimensionless potential; and predicting the expansive force of the bentonite to be detected through the bentonite expansive force prediction model. According to the invention, the bentonite expansive force can be accurately predicted by utilizing the pore ratio of the bentonite through the physical and chemical theory of the soil and the contribution of charge potential between bentonite layers to the expansive force of the bentonite. The invention provides a material selection method with high reliability and high precision for selecting the barrier buffer material for nuclide migration in a high-radioactive nuclear waste disposal library.
Optionally, the obtaining of the porosity ratio of the bentonite to be detected includes the following steps: obtaining the dry density of the bentonite to be detected; and obtaining the porosity ratio of the bentonite to be detected by utilizing the dry density, wherein the porosity ratio meets the following formula:
Figure DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE002
the porosity ratio of the bentonite to be measured is shown,
Figure DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure DEST_PATH_IMAGE004
the weight of the bentonite to be tested is shown,
Figure DEST_PATH_IMAGE005
the dry density of the bentonite to be tested is shown.
Optionally, the obtaining of the compensated dimensionless potential of the bentonite to be tested and the dimensionless potential of the ion hydration energy by using the porosity ratio includes the following steps: constructing an expansive force double electric layer model according to the physical characteristics of the bentonite to be detected; acquiring a function relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using an expansive force double electric layer model; obtaining the dimensionless potential of the bentonite to be detected according to the functional relation between the dimensionless potential and the pore ratio; acquiring a functional relation between the compensation dimensionless potential of the ion hydration energy and the pore ratio according to the microstructure change of the bentonite to be detected after the ion hydration; and acquiring the compensated non-dimensional potential of the ion hydration energy according to the functional relation between the compensated non-dimensional potential and the pore ratio.
Optionally, the expansion force electric double layer model satisfies the following formula:
Figure DEST_PATH_IMAGE007
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE008
the force of the expansion is indicated by,
Figure DEST_PATH_IMAGE009
the ion concentration of the bentonite pore water is shown,
Figure DEST_PATH_IMAGE010
the values of the boltzmann constants are expressed,
Figure DEST_PATH_IMAGE011
and u represents a dimensionless potential.
Optionally, the dimensionless potential as a function of void ratio satisfies the following equation:
Figure 100002_DEST_PATH_IMAGE012
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE013
and
Figure DEST_PATH_IMAGE014
the coefficient of fit is represented by the value of,
Figure 523372DEST_PATH_IMAGE002
it is indicated that the ratio of the pores,
Figure DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure 962443DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure 393031DEST_PATH_IMAGE004
the gravity of the bentonite to be tested is shown, and S is the specific surface area of the bentonite particles to be tested.
Optionally, the functional relationship of the compensated dimensionless potential to void ratio satisfies the following equation:
Figure DEST_PATH_IMAGE016
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE017
representing a compensated non-dimensional potential that is,
Figure 298670DEST_PATH_IMAGE002
the void ratio is indicated.
Optionally, the obtaining a functional relationship between the dimensionless potential and the pore ratio of the bentonite to be tested by using the expansive force electric double layer model includes the following steps: setting a plurality of expansion force values; obtaining dimensionless potential and bentonite layer spacing corresponding to the expansion force value through the expansion force double electric layer model; acquiring a functional relation between the space of the bentonite layer and the pore ratio; acquiring a pore ratio corresponding to the expansion force value by utilizing the functional relation between the space of the bentonite layer and the pore ratio; fitting a functional relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using the dimensionless potential and the pore ratio corresponding to the multiple groups of expansion force values.
Optionally, the functional relationship of the dimensionless potential to the porosity ratio comprises: the dimensionless potential of the german Bavaria bentonite as a function of the void ratio satisfies the following formula:
Figure DEST_PATH_IMAGE018
wherein, the first and the second end of the pipe are connected with each other,
Figure 100002_DEST_PATH_IMAGE019
representing the dimensionless potential of Bavaria bentonite in germany,
Figure DEST_PATH_IMAGE020
represents the porosity ratio of the German Bavaria bentonite,
Figure 316174DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure DEST_PATH_IMAGE021
represents the specific gravity of German Bavaria bentonite,
Figure DEST_PATH_IMAGE022
representing the severity of Bavaria bentonite in germany,
Figure DEST_PATH_IMAGE023
represents the specific surface area of German Bavaria bentonite particles; the functional relationship of the dimensionless potential and the porosity ratio of Spanish S-2 bentonite satisfies the following formula:
Figure DEST_PATH_IMAGE024
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE025
represents the dimensionless potential of spanish S-2 bentonite,
Figure DEST_PATH_IMAGE026
represents the porosity ratio of Spain S-2 bentonite,
Figure 932838DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure DEST_PATH_IMAGE027
represents the specific gravity of Spain S-2 bentonite,
Figure DEST_PATH_IMAGE028
indicating the severity of spanish S-2 bentonite,
Figure DEST_PATH_IMAGE029
representing the specific surface area of the Spanish S-2 bentonite particles.
Optionally, the bentonite expansive force prediction model satisfies the following formula:
Figure DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 329184DEST_PATH_IMAGE008
the force of the expansion is indicated by,
Figure 264779DEST_PATH_IMAGE009
the ion concentration of the bentonite pore water is shown,
Figure 582628DEST_PATH_IMAGE010
the values of the boltzmann constants are expressed,
Figure 989601DEST_PATH_IMAGE011
Figure DEST_PATH_IMAGE032
denotes a natural constant, u denotes a dimensionless potential,
Figure 506033DEST_PATH_IMAGE017
representing the compensated dimensionless potential.
In a second aspect, the present invention further provides a bentonite expansive force prediction system considering ionic hydration energy, which includes an input device, a processor, a memory and an output device, where the input device, the processor, the memory and the output device are connected with each other, where the memory is used for storing a computer program, the computer program includes program instructions, and the processor is configured to call the program instructions to execute the bentonite expansive force prediction method considering ionic hydration energy according to the first aspect of the present invention. The system is compact in structure and high in applicability, operation efficiency is greatly improved, and prediction of simulated bentonite expansive force is quickly and accurately realized by considering the bentonite expansive force prediction method of the ionic hydration energy.
Drawings
FIG. 1 is a flow chart of a method for predicting swelling force of bentonite according to the invention, taking ionic hydration energy into consideration;
FIG. 2 is a comparison graph of the experimental test value, the predicted value without considering the hydration energy and the predicted value considering the hydration energy of the German Bavaria bentonite obtained by the invention;
FIG. 3 is a comparison graph of experimental test values, predicted values without and with regard to hydration energy of Spain S-2 bentonite obtained in the present invention;
fig. 4 is a structural diagram of a bentonite swelling force prediction system in consideration of ionic hydration energy according to the invention.
Detailed Description
Specific embodiments of the present invention will be described in detail below, and it should be noted that the embodiments described herein are merely illustrative and are not intended to limit the present invention. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be apparent to those of ordinary skill in the art that: it is not necessary to employ these specific details to practice the present invention. In other instances, well-known circuits, software, or methods have not been described in detail so as not to obscure the present invention.
Throughout the specification, reference to "one embodiment," "an embodiment," "one example," or "an example" means: the particular features, structures, or characteristics described in connection with the embodiment or example are included in at least one embodiment of the invention. Thus, the appearances of the phrases "in one embodiment," "in an embodiment," "one example" or "an example" in various places throughout this specification are not necessarily all referring to the same embodiment or example. Furthermore, the particular features, structures, or characteristics may be combined in any suitable combination and/or sub-combination in one or more embodiments or examples. Further, those of ordinary skill in the art will appreciate that the illustrations provided herein are for illustrative purposes and are not necessarily drawn to scale.
Referring to fig. 1, in one embodiment, the present invention provides a method for predicting swelling force of bentonite considering ion hydration energy, including the following steps:
s1, obtaining the porosity ratio of the bentonite to be detected.
The pore ratio is the ratio of the pore volume of the bentonite soil body to the volume of the soil solid, and the larger the pore ratio is, the more pores are, the less dense the soil is. In practical engineering, the relation between the dry density (or porosity ratio) and the expansive force of bentonite is generally used to characterize the expansive force trend of bentonite, and the dry density and the porosity ratio of the bentonite can be converted into each other through the specific gravity of the bentonite.
In an optional embodiment, the obtaining of the porosity ratio of the bentonite to be measured in step S1 includes the following steps: obtaining the dry density of the bentonite to be detected; and obtaining the porosity ratio of the bentonite to be detected by utilizing the dry density, wherein the porosity ratio meets the following formula:
Figure 222316DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 89778DEST_PATH_IMAGE002
the porosity ratio of the bentonite to be measured is shown,
Figure 798977DEST_PATH_IMAGE005
the dry density of the bentonite to be tested is shown,
Figure 373177DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure 650575DEST_PATH_IMAGE004
and the gravity of the bentonite to be detected is represented, and in detail, the gravity of the bentonite to be detected represents the product of the density of the bentonite to be detected and the gravity acceleration.
S2, obtaining the compensation dimensionless potential of the dimensionless potential and the ion hydration energy of the bentonite to be detected by utilizing the pore ratio.
In an optional embodiment, the step S2 of obtaining the compensated non-dimensional potential of the non-dimensional potential and the ion hydration energy of the bentonite to be detected by using the porosity ratio includes the following steps: constructing an expansive force double electric layer model according to the physical characteristics of the bentonite to be detected; acquiring a function relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using an expansive force double electric layer model; obtaining the dimensionless potential of the bentonite to be detected according to the functional relation between the dimensionless potential and the pore ratio; acquiring a functional relation between the compensation dimensionless potential of the ion hydration energy and the pore ratio according to the microstructure change of the bentonite to be detected after the ion hydration; and acquiring the compensated non-dimensional potential of the ion hydration energy according to the functional relation between the compensated non-dimensional potential and the pore ratio.
In an alternative embodiment, the expansion force electric double layer model satisfies the following formula:
Figure DEST_PATH_IMAGE034
wherein the content of the first and second substances,
Figure 83962DEST_PATH_IMAGE008
the force of the expansion is indicated by,
Figure 737797DEST_PATH_IMAGE009
the ion concentration of the bentonite pore water is shown,
Figure 166504DEST_PATH_IMAGE010
denotes the Boltzmann constant, T denotes the temperature in Kelvin, room temperature 25 deg.C
Figure 972393DEST_PATH_IMAGE011
Wherein
Figure DEST_PATH_IMAGE035
The unit of (d) is kelvin, and u represents a dimensionless potential. The dimensionless potential varies in magnitude for different types of bentonite. The expansion force double electric layer model does not consider the contribution of hydration energy of ions in the bentonite after hydration to the expansion force of the bentonite, and meanwhile, for different types of bentonite, the expansion force double electric layer model is used and substituted into corresponding parameters for operation, so that the corresponding function relation between the dimensionless potential and the pore ratio is obtained.
In an optional embodiment, the obtaining a functional relationship between the dimensionless potential and the void ratio of the bentonite to be tested by using the expansive force double electric layer model includes the following steps: setting a plurality of expansion force values; obtaining dimensionless electric potential and bentonite layer distance corresponding to the expansion force value through the expansion force double electric layer model; acquiring a functional relation between the space of the bentonite layer and the pore ratio; acquiring a pore ratio corresponding to the expansion force value by utilizing the functional relation between the space of the bentonite layer and the pore ratio; fitting a functional relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using the dimensionless potential and the pore ratio corresponding to the multiple groups of expansion force values.
In this embodiment, through the expansion force double electric layer model, the dimensionless electric potential and the bentonite layer interval corresponding to the expansion force value can be obtained through what is proposed by Olpehen in "thermomonynamics of interlayer absorption of water in clays
Figure DEST_PATH_IMAGE037
Figure DEST_PATH_IMAGE039
The following equation is obtained in combination with the bentonite properties:
Figure DEST_PATH_IMAGE041
Figure 407922DEST_PATH_IMAGE039
using this formula, the value of y is found, where y = z, where y represents the dimensionless potential at a location x from the bentonite layer, where x =0, and z represents the value of y where x =0, the value of z is obtained,
Figure DEST_PATH_IMAGE042
is a distance function between bentonite layers, K represents a double electric layer coefficient, B represents ion exchange energy, S represents the specific surface area of the bentonite particles to be measured,
Figure DEST_PATH_IMAGE043
which represents the dielectric constant of a vacuum,
Figure DEST_PATH_IMAGE044
which represents the relative dielectric constant of the material,
Figure 475235DEST_PATH_IMAGE009
the ion concentration of the bentonite pore water is shown,
Figure 820766DEST_PATH_IMAGE010
the values of the boltzmann constants are expressed,
Figure 66065DEST_PATH_IMAGE011
Figure DEST_PATH_IMAGE045
represents the dielectric constant of the fluid; and substituting the z value into the u value obtained by substituting the expansion force double-layer model into the set expansion force value according to the z value in the Olpehen article
Figure DEST_PATH_IMAGE047
Double-layer dielectric function equation of one-piece combined bentonite
Figure DEST_PATH_IMAGE048
Obtaining the bentonite layer interval corresponding to the set expansion force value, wherein the double-layer dielectric function equation
Figure DEST_PATH_IMAGE049
Which represents the amount of the cell charge,
Figure DEST_PATH_IMAGE050
represents the average number of ionic valences in the bentonite. In this example, the bentonite layer spacing as a function of void ratio satisfies the following equation:
Figure DEST_PATH_IMAGE051
wherein, in the process,
Figure 192153DEST_PATH_IMAGE002
it is expressed as a ratio of the pores,
Figure 187790DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure 997615DEST_PATH_IMAGE004
the gravity of the bentonite to be measured is shown, and S represents the specific surface area of the bentonite particles to be measured. In actual engineering, basic parameters such as specific gravity, particle specific surface area and the like of the bentonite to be tested can be obtained by testing the basic physical properties of the bentonite.
In this embodiment, the dimensionless potential as a function of void ratio satisfies the following equation:
Figure 990978DEST_PATH_IMAGE012
wherein the content of the first and second substances,
Figure 542045DEST_PATH_IMAGE013
and
Figure 433385DEST_PATH_IMAGE014
the coefficient of fit is represented by the value of,
Figure 691191DEST_PATH_IMAGE002
it is expressed as a ratio of the pores,
Figure 917773DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure DEST_PATH_IMAGE052
the specific gravity of the bentonite to be measured is shown,
Figure 34764DEST_PATH_IMAGE004
the gravity of the bentonite to be measured is shown, and S represents the specific surface area of the bentonite particles to be measured. Fitting coefficient of different types of bentonite
Figure 309888DEST_PATH_IMAGE013
And
Figure 687780DEST_PATH_IMAGE014
different. The function relation of the dimensionless potential and the pore ratio is combined with the function relation of the bentonite layer space and the pore ratio, so that an empirical expression of the dimensionless potential and the bentonite layer space can be obtained, and the empirical expression meets the following formula:
Figure DEST_PATH_IMAGE053
wherein, in the step (A),
Figure DEST_PATH_IMAGE054
the bentonite layer spacing is indicated. The method comprises the steps of setting a plurality of expansion force values, obtaining corresponding dimensional potential and bentonite layer distance by utilizing an expansion force double-electric-layer model, selecting a fitting type by utilizing curve fitting software, and obtaining a corresponding fitting curve and fitting coefficients by inputting specific values of the dimensionless potential and the bentonite layer distance. The predicted expansion force on the fitted curve and the actually measured expansion force are now at low dry density (dry density less than 1.5 g/cm) 3 ) The error is small but at high dry density (greater than 1.5 g/cm) 3 ) The error is large.
In an optional embodiment, the step of obtaining a functional relationship between the dimensionless potential and the pore ratio of the bentonite to be tested by using the expansive force electric double layer model in the previous embodiment obtains the following functional relationship between the dimensionless potential and the pore ratio, including: the dimensionless potential of the german Bavaria bentonite as a function of the void ratio satisfies the following formula:
Figure 960629DEST_PATH_IMAGE018
wherein the content of the first and second substances,
Figure 679816DEST_PATH_IMAGE019
representing the dimensionless potential of Bavaria bentonite in germany,
Figure 555368DEST_PATH_IMAGE020
represents the porosity ratio of the German Bavaria bentonite,
Figure 459870DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure 700358DEST_PATH_IMAGE021
represents the specific gravity of the German Bavaria bentonite,
Figure 713314DEST_PATH_IMAGE022
representing the severity of Bavaria bentonite in germany,
Figure 251611DEST_PATH_IMAGE023
the specific surface area of German Bavaria bentonite particles is expressed; the functional relationship of the dimensionless potential and the porosity ratio of Spanish S-2 bentonite satisfies the following formula:
Figure 338516DEST_PATH_IMAGE024
wherein, the first and the second end of the pipe are connected with each other,
Figure 77802DEST_PATH_IMAGE025
represents the dimensionless potential of spanish S-2 bentonite,
Figure 922261DEST_PATH_IMAGE026
represents the porosity ratio of Spain S-2 bentonite,
Figure 874037DEST_PATH_IMAGE015
the electric double layer coefficient is expressed,
Figure 81027DEST_PATH_IMAGE027
represents the specific gravity of Spanish S-2 bentonite,
Figure 351734DEST_PATH_IMAGE028
indicating the severity of spanish S-2 bentonite,
Figure 808123DEST_PATH_IMAGE029
representing the specific surface area of the Spanish S-2 bentonite particles. In the embodiment, dimensionless potentials corresponding to German Bavaria bentonite and Spain S-2 bentonite are obtained, and a foundation is laid for predicting the expansive force of German Bavaria bentonite and Spain S-2 bentonite by a subsequent bentonite expansive force prediction model.
In an optional embodiment, when ions in the bentonite are hydrated in water, a water film is formed around the ions, so that the space between bentonite layers is increased, and energy is released after the ions are hydrated, so that the microstructure of the bentonite is influenced, and when the swelling force of the bentonite considers the hydration energy of the ions, the functional relation between the compensation dimensionless potential and the pore ratio meets the following formula:
Figure 766851DEST_PATH_IMAGE016
wherein the content of the first and second substances,
Figure 766031DEST_PATH_IMAGE017
a compensated non-dimensional potential is represented,
Figure 581541DEST_PATH_IMAGE002
indicating the void ratio.
And S3, constructing a bentonite expansive force prediction model according to the dimensionless potential and the compensation dimensionless potential.
In an alternative embodiment, the bentonite expansive force prediction model satisfies the following formula:
Figure DEST_PATH_IMAGE056
wherein, the first and the second end of the pipe are connected with each other,
Figure 384281DEST_PATH_IMAGE008
the force of the expansion is indicated by,
Figure 943438DEST_PATH_IMAGE009
the ion concentration of the bentonite pore water is shown,
Figure 859441DEST_PATH_IMAGE010
the values of the boltzmann constants are expressed,
Figure 721218DEST_PATH_IMAGE011
Figure 417779DEST_PATH_IMAGE032
denotes a natural constant, in this example
Figure DEST_PATH_IMAGE057
Value taking
Figure DEST_PATH_IMAGE058
And u represents a non-dimensional potential,
Figure 669375DEST_PATH_IMAGE017
representing the compensated dimensionless potential. The bentonite expansive force prediction model considers the contribution of a hydraulic coupling field of hydration energy of ions in bentonite after hydration to the bentonite expansive force, and then compensates the dimensionless potential by utilizing the contribution, so that the prediction error of the traditional bentonite expansive force is reduced, and the predicted value and the measured value have better consistency.
And S4, predicting the expansive force of the bentonite to be detected through the bentonite expansive force prediction model.
In an alternative embodiment, referring to fig. 2 and 3, predicted values of swelling force of Bavaria bentonite in germany and predicted values of swelling force of Bavaria bentonite in germany are obtained through steps S1 to S4, respectivelySpain S-2 Bentonite expansive force prediction value, wherein, the dimensionless potential of Germany Bavaria Bentonite is
Figure 111989DEST_PATH_IMAGE018
The dimensionless potential of Spanish S-2 bentonite is
Figure 472563DEST_PATH_IMAGE024
Respectively, as a dimensionless potential of Bavaria bentonite in Germany of
Figure DEST_PATH_IMAGE060
The dimensionless potential of Spanish S-2 bentonite is
Figure DEST_PATH_IMAGE062
Wherein, in the process,
Figure DEST_PATH_IMAGE063
representing the space between germany Bavaria bentonite layers,
Figure DEST_PATH_IMAGE064
representing spanish S-2 bentonite layer spacing. Fig. 2 is an experimental test value, a predicted value without considering the water energy and a predicted value considering the water energy of Bavaria bentonite in germany, wherein the predicted value considering the water energy is obtained by the method of the present invention. Fig. 3 is an experimental test value, a predicted value without considering water energy, and a predicted value considering water energy of spanish S-2 bentonite, wherein the predicted value considering water energy is obtained by the method of the present invention. As shown in fig. 2 and 3, the ion hydration energy is considered when predicting the bentonite in the invention, so that the consistency degree of the predicted value and the experimental test value is higher than the consistency degree of the predicted value without considering the ion hydration energy, that is, the bentonite expansive force prediction model of the invention has good prediction performance on the expansive force of the bentonite with low dry density and high dry density.
According to the invention, the bentonite expansive force can be accurately predicted by utilizing the porosity ratio of the bentonite through combining the physical and chemical theory of the soil and the contribution of charge potential between bentonite layers to the expansive force of the bentonite. Meanwhile, the invention provides a material selection method with high reliability and high precision for selecting the barrier buffer material for nuclide migration in the high radioactive nuclear waste disposal library.
Referring to fig. 4, the present invention further provides a bentonite expansive force prediction system considering ion hydration energy, which includes an input device, a processor, a memory and an output device, where the input device, the processor, the memory and the output device are connected to each other, where the memory is used for storing a computer program, the computer program includes program instructions, and the processor is configured to call the program instructions to execute the bentonite expansive force prediction method considering ion hydration energy according to the first aspect of the present invention. The system is compact in structure and high in applicability, the operation efficiency is greatly improved, and prediction of simulated bentonite expansion force is quickly and accurately realized by considering the bentonite expansion force prediction method of the ionic hydration energy.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the present invention, and they should be construed as being included in the following claims and description.

Claims (6)

1. The method for predicting the swelling force of the bentonite by considering the ionic hydration energy is characterized by comprising the following steps of:
the method for obtaining the porosity ratio of the bentonite to be detected comprises the following steps:
obtaining the dry density of the bentonite to be detected;
and obtaining the porosity ratio of the bentonite to be detected by utilizing the dry density, wherein the porosity ratio meets the following formula:
Figure 904381DEST_PATH_IMAGE001
wherein, the first and the second end of the pipe are connected with each other,
Figure 861973DEST_PATH_IMAGE002
the porosity ratio of the bentonite to be tested is shown,
Figure 333405DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure 387949DEST_PATH_IMAGE004
the weight of the bentonite to be tested is shown,
Figure 237700DEST_PATH_IMAGE005
expressing the dry density of the bentonite to be detected;
and acquiring the compensation dimensionless potential and the ion hydration energy of the bentonite to be detected by utilizing the pore ratio, wherein the functional relation between the dimensionless potential and the compensation dimensionless potential and the pore ratio respectively meets the following formula:
Figure 10484DEST_PATH_IMAGE006
Figure 20028DEST_PATH_IMAGE007
wherein u represents a dimensionless potential,
Figure 335603DEST_PATH_IMAGE008
representing the compensated non-dimensional electrical potential,
Figure 201928DEST_PATH_IMAGE009
and
Figure 665270DEST_PATH_IMAGE010
the fitting coefficient is expressed as a ratio of,
Figure 665456DEST_PATH_IMAGE002
it is expressed as a ratio of the pores,
Figure 694592DEST_PATH_IMAGE011
the electric double layer coefficient is expressed,
Figure 200660DEST_PATH_IMAGE003
the specific gravity of the bentonite to be measured is shown,
Figure 557823DEST_PATH_IMAGE004
representing the gravity of the bentonite to be detected, and S represents the specific surface area of the bentonite particles to be detected;
according to the dimensionless potential and the compensation dimensionless potential, a bentonite expansive force prediction model is constructed, the bentonite expansive force prediction model has good prediction performance on the expansive force of the bentonite with low dry density and high dry density, and the bentonite expansive force prediction model meets the following formula:
Figure DEST_PATH_IMAGE012
wherein the content of the first and second substances,
Figure 971487DEST_PATH_IMAGE013
the force of the expansion is indicated by,
Figure 746807DEST_PATH_IMAGE015
the ion concentration of the bentonite pore water is shown,
Figure 220513DEST_PATH_IMAGE016
which represents the boltzmann constant, represents,
Figure 392869DEST_PATH_IMAGE017
Figure 220011DEST_PATH_IMAGE018
denotes a natural constant, u denotes a dimensionless potential,
Figure DEST_PATH_IMAGE019
representing the compensated dimensionless potential;
and predicting the expansive force of the bentonite to be tested through the bentonite expansive force prediction model, so that the consistency degree of the predicted value and the experimental test value is higher than that of the predicted value without considering the ion hydration energy.
2. The method for predicting the expansive force of bentonite according to claim 1, wherein the step of obtaining the dimensionless potential and the compensated dimensionless potential of the ionic hydration energy of the bentonite to be tested by using the porosity ratio comprises the following steps:
constructing an expansive force double electric layer model according to the physical characteristics of the bentonite to be detected;
acquiring a function relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using an expansive force double electric layer model;
obtaining the dimensionless potential of the bentonite to be detected according to the functional relation between the dimensionless potential and the pore ratio;
acquiring a functional relation between the compensation dimensionless potential of the ion hydration energy and the pore ratio according to the microstructure change of the bentonite to be detected after the ion hydration;
and acquiring the compensated non-dimensional potential of the ion hydration energy according to the functional relation between the compensated non-dimensional potential and the pore ratio.
3. The method for predicting swelling force of bentonite according to claim 2, wherein the swelling force electric double layer model satisfies the following formula:
Figure 145110DEST_PATH_IMAGE020
wherein the content of the first and second substances,
Figure 992980DEST_PATH_IMAGE013
the force of the expansion is indicated by,
Figure 714949DEST_PATH_IMAGE014
the ion concentration of the bentonite pore water is shown,
Figure 673677DEST_PATH_IMAGE016
representing boltzmann's constant, T =298.15, u represents a dimensionless potential.
4. The method for predicting the swelling force of the bentonite according to claim 2, wherein the step of obtaining the functional relation between the dimensionless potential and the pore ratio of the bentonite to be tested by using a swelling force double electric layer model comprises the following steps:
setting a plurality of expansion force values;
obtaining dimensionless electric potential and bentonite layer distance corresponding to the expansion force value through the expansion force double electric layer model;
acquiring a functional relation between the space of the bentonite layer and the pore ratio;
acquiring a pore ratio corresponding to the expansion force value by utilizing the functional relation between the space of the bentonite layer and the pore ratio;
fitting a functional relation between the dimensionless potential and the pore ratio of the bentonite to be detected by using the dimensionless potential and the pore ratio corresponding to the multiple groups of expansion force values.
5. The method of predicting bentonite expansive force according to claim 4, wherein the function of the dimensionless potential and the porosity ratio comprises:
the dimensionless potential of the German Bavaria bentonite is in function of the porosity ratio and satisfies the following formula:
Figure 141699DEST_PATH_IMAGE021
wherein the content of the first and second substances,
Figure 222787DEST_PATH_IMAGE022
representing the dimensionless potential of Bavaria bentonite in germany,
Figure 789642DEST_PATH_IMAGE023
represents the porosity ratio of the German Bavaria bentonite,
Figure 348799DEST_PATH_IMAGE011
the electric double layer coefficient is expressed,
Figure 264802DEST_PATH_IMAGE024
represents the specific gravity of German Bavaria bentonite,
Figure 126579DEST_PATH_IMAGE025
representing the severity of Bavaria bentonite in germany,
Figure 557561DEST_PATH_IMAGE026
represents the specific surface area of German Bavaria bentonite particles;
the function relation of the dimensionless potential and the porosity ratio of Spanish S-2 bentonite satisfies the following formula:
Figure 858092DEST_PATH_IMAGE027
wherein the content of the first and second substances,
Figure 81132DEST_PATH_IMAGE028
represents the dimensionless potential of spanish S-2 bentonite,
Figure 972864DEST_PATH_IMAGE029
represents the porosity ratio of Spain S-2 bentonite,
Figure 94404DEST_PATH_IMAGE011
expressing electric double layer coefficient,
Figure 136309DEST_PATH_IMAGE030
Represents the specific gravity of Spanish S-2 bentonite,
Figure 558063DEST_PATH_IMAGE031
indicating the severity of spanish S-2 bentonite,
Figure 89539DEST_PATH_IMAGE032
representing the specific surface area of the Spain S-2 bentonite particles.
6. A bentonite expansive force prediction system considering ionic hydration energy, which is characterized by comprising an input device, a processor, a memory and an output device, wherein the input device, the processor, the memory and the output device are connected with each other, the memory is used for storing a computer program, the computer program comprises program instructions, and the processor is configured to call the program instructions to execute the bentonite expansive force prediction method considering ionic hydration energy according to any one of claims 1 to 5.
CN202211243933.4A 2022-10-12 2022-10-12 Bentonite expansive force prediction method and system considering ionic hydration energy Active CN115310026B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211243933.4A CN115310026B (en) 2022-10-12 2022-10-12 Bentonite expansive force prediction method and system considering ionic hydration energy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211243933.4A CN115310026B (en) 2022-10-12 2022-10-12 Bentonite expansive force prediction method and system considering ionic hydration energy

Publications (2)

Publication Number Publication Date
CN115310026A CN115310026A (en) 2022-11-08
CN115310026B true CN115310026B (en) 2022-12-30

Family

ID=83867610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211243933.4A Active CN115310026B (en) 2022-10-12 2022-10-12 Bentonite expansive force prediction method and system considering ionic hydration energy

Country Status (1)

Country Link
CN (1) CN115310026B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115656481B (en) * 2022-12-14 2023-03-28 海南浙江大学研究院 Method and system for representing water holding capacity of soil microstructure

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007209903A (en) * 2006-02-09 2007-08-23 Shimizu Corp Plug structure in facility for burying and disposing waste
JP2014048099A (en) * 2012-08-30 2014-03-17 Shimizu Corp Behavior prediction method of water absorption expansive clay
CN106680330A (en) * 2017-01-22 2017-05-17 东南大学 Method for carrying out on-site evaluation on expansion property of expansive soil by utilizing electrical resistivity
CN107014730A (en) * 2017-03-15 2017-08-04 中国矿业大学 A kind of method for simulating true groundwater erosion Nuclear Waste Repository padded coaming
CN109187213A (en) * 2018-09-21 2019-01-11 同济大学 Concrete-high-pressure solid bentonite combined system condition simulation experiment method and device
CN110849924A (en) * 2019-12-02 2020-02-28 武汉科技大学 Bentonite quality evaluation method
CN112903969A (en) * 2021-01-26 2021-06-04 同济大学 Expansive soil strength evaluation model considering physical and chemical effects
CN114169154A (en) * 2021-11-25 2022-03-11 欣旺达电动汽车电池有限公司 Expansive force prediction method, expansive force prediction system, electronic device, and storage medium

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7022398B2 (en) * 2003-09-18 2006-04-04 Amcol International Corporation Moisture-impervious water-swellable clay-containing “water-stop” composition containing a water-penetrable coating
CN106815391A (en) * 2016-08-29 2017-06-09 中国辐射防护研究院 A kind of bentonite Radionuclide Migration appraisal procedure for being based on two Scale Models
CN112033891B (en) * 2020-08-03 2023-09-12 南京交通职业技术学院 Novel expansive force measuring method for expansive soil
CN112146983B (en) * 2020-09-24 2023-11-03 大连民族大学 Dimensionless soil body compression coefficient representation method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007209903A (en) * 2006-02-09 2007-08-23 Shimizu Corp Plug structure in facility for burying and disposing waste
JP2014048099A (en) * 2012-08-30 2014-03-17 Shimizu Corp Behavior prediction method of water absorption expansive clay
CN106680330A (en) * 2017-01-22 2017-05-17 东南大学 Method for carrying out on-site evaluation on expansion property of expansive soil by utilizing electrical resistivity
CN107014730A (en) * 2017-03-15 2017-08-04 中国矿业大学 A kind of method for simulating true groundwater erosion Nuclear Waste Repository padded coaming
CN109187213A (en) * 2018-09-21 2019-01-11 同济大学 Concrete-high-pressure solid bentonite combined system condition simulation experiment method and device
CN110849924A (en) * 2019-12-02 2020-02-28 武汉科技大学 Bentonite quality evaluation method
CN112903969A (en) * 2021-01-26 2021-06-04 同济大学 Expansive soil strength evaluation model considering physical and chemical effects
CN114169154A (en) * 2021-11-25 2022-03-11 欣旺达电动汽车电池有限公司 Expansive force prediction method, expansive force prediction system, electronic device, and storage medium

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
双电层理论与高庙子膨润土的体变特征;叶为民等;《岩土力学》;20090731;第30卷(第7期);第1899-1903页 *
高压实膨润土膨胀力预测研究;陈永贵 等;《同济大学学报(自然科学版)》;20181231;第46卷(第12期);第1628-1636页 *

Also Published As

Publication number Publication date
CN115310026A (en) 2022-11-08

Similar Documents

Publication Publication Date Title
Zhang et al. Deformation response of an existing tunnel to upper excavation of foundation pit and associated dewatering
Zhang et al. A new soil-water characteristic curve model for unsaturated loess based on wetting-induced pore deformation
Villar et al. Influence of dry density and water content on the swelling of a compacted bentonite
Tarantino et al. Making unsaturated soil mechanics accessible for engineers: Preliminary hydraulic–mechanical characterisation & stability assessment
Chijimatsu et al. Experiment and validation of numerical simulation of coupled thermal, hydraulic and mechanical behaviour in the engineered buffer materials
Mendes et al. Influence of initial water content on the mechanical behavior of unsaturated sandy clay soil
Gerard et al. Study of the soil–atmosphere moisture exchanges through convective drying tests in non‐isothermal conditions
Fattah et al. Influence of soil suction on swelling pressure of bentonite-sand mixtures
Xiong et al. Soil–water–air fully coupling finite element analysis of slope failure in unsaturated ground
CN115310026B (en) Bentonite expansive force prediction method and system considering ionic hydration energy
Li et al. Volume change behavior of two compacted clayey soils under hydraulic and mechanical loadings
Thomas et al. Water infiltration into a large-scale in-situ experiment in an underground research laboratory
Pedarla et al. Evaluation of swell behavior of expansive clays from internal specific surface and pore size distribution
Al-Dakheeli et al. Interrelationship between elastic deformation and soil-water characteristic curve of expansive soils
Bruch Laboratory study of evaporative fluxes in homogeneous and layered soils
Ma et al. Constitutive model of unsaturated soils considering the effect of intergranular physicochemical forces
Gallé et al. Gas permeability of unsaturated cement-based materials: application of a multi-scale network model
Jia et al. Modified Duncan-Chang constitutive model for modeling supported excavations in granular soils
Pham et al. Effect of temperature on the soil–water retention characteristics in unsaturated soils: Analytical and experimental approaches
Lee et al. Suction and water uptake in unsaturated compacted bentonite
Chen et al. Coupled thermohydromechanical modeling of the full-scale in situ test “Prototype repository”
Molina-Gomez et al. Internal erosion of a 0–5 mm crushed sand in a rigid wall-permeameter: Experimental methods and results
Zhang et al. Investigation of suction effects due to stress release with compacted MX80 bentonite
Jian-Zhong et al. One-dimensional viscous behavior of clay and its constitutive modeling
Bharat et al. Prediction of compressibility behavior of clayey soils of different plasticity for containment applications at large consolidation pressures

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant