Method, system, device and storage medium for obtaining soil saturated hydraulic conductivity

Abstract

The application discloses a soil saturated permeability coefficient obtaining method, system, device and storage medium, belongs to the field of geotechnical engineering, hydrology and environmental geology, and the method comprises the steps of: calculating the effective particle size of soil based on particle size distribution data; obtaining soil porosity; and obtaining the saturated permeability coefficient based on solid-liquid two-phase fluid mechanics equations. By incorporating porosity, the effective particle size of particle distribution and water flow density and other parameters are embodied, the synergistic effect of multiple parameters is quantified and comprehensive influence is achieved, the single parameter dependency of traditional empirical formula is broken through; meanwhile, the formula is derived based on fluid mechanics theory, the simplification defect of the theoretical model is broken through, the calculation efficiency is significantly improved on the basis of ensuring the theoretical rigor. By incorporating porosity and effective particle size based on particle size distribution, and the square form of the two, the difference in the structure of non-homogeneous soil is amplified, the difference in the water flow path is reflected in combination with the fluid mechanics theory, so that the difference in the seepage characteristics is accurately quantified.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of geotechnical engineering, hydrogeology and environmental geology, and particularly relates to a soil saturated permeability coefficient obtaining method, system, device and storage medium. BACKGROUND

[0002] The soil saturated permeability coefficient is a key parameter for characterizing the seepage capacity of porous media, and has important engineering significance in the fields of slope stability evaluation, groundwater pollution prevention and control, and underground engineering seepage prevention design. However, the saturated permeability coefficient determination technology generally faces systematic defects such as excessive simplification of theoretical models, long test cycle, and single parameter characterization, which seriously restricts the fine design and risk control ability of seepage related engineering. The existing theoretical prediction models are mostly based on empirical formulas (such as Hazen formula, Kozeny-Carman equation), and the core defect of the traditional empirical formula is that it only takes the effective particle size as a single variable. In addition, the traditional empirical formula derivation does not strictly follow the fluid mechanics theory, and it is difficult to characterize the nonlinear seepage characteristics under high hydraulic gradient. SUMMARY

[0003] The purpose of the present application is to provide a soil saturated permeability coefficient obtaining method, system, device and storage medium to solve the problem of single parameter characterization of traditional empirical formula, non-strict adherence to fluid mechanics theory in traditional empirical formula derivation, and difficulty in characterizing nonlinear seepage characteristics under high hydraulic gradient.

[0004] The present application is implemented in the following way: a soil body saturated permeability coefficient obtaining method comprises:

[0005] calculating the effective particle size of the soil based on the particle size distribution data;

[0006] obtaining the soil porosity;

[0007] obtaining the saturated permeability coefficient based on the solid-liquid two-phase fluid mechanics equation.

[0008] Optionally, in some embodiments of the present application, the particle size distribution data includes the particle size value in the soil particle size distribution and the mass fraction corresponding to the particle size in the soil particle size distribution .

[0009] Optionally, in some embodiments of the present application, based on the particle size distribution data, the effective particle size is calculated according to the specific surface area of the soil particles, and the reciprocal of the effective particle size is equal to the sum of the specific surface area of the soil particles, that is, the calculation formula of the effective particle size is as follows: ;

[0010] In the formula: is the effective particle size of the soil, unit: ; The particle size value in soil particle size distribution, unit: ; This represents the mass fraction corresponding to the particle size in the soil particle size distribution, in dimensionless form. The group number refers to the group 1, group 2, ..., group n, which are different particle sizes and their corresponding mass fractions. The unit is dimensionless. n is the total number of soil particle size distribution groups, which represents the number of groups obtained after dividing the soil according to the particle size range. The unit is dimensionless.

[0011] Optionally, in some embodiments of this application, the formula for calculating soil porosity is as follows:

[0012] ;

[0013] In the formula: Soil dry density, unit: ; Soil porosity, unit: dimensionless; The intrinsic density of soil particles, in units of: Its value is 2650 ; The drying quality of soil collected in the field, unit: ; Total volume of soil collected in the field, unit: .

[0014] Optionally, in some embodiments of this application, the method for obtaining soil porosity is as follows:

[0015] Provide soil dry density; the formula for calculating soil dry density is as follows:

[0016] ;

[0017] In the formula:

[0018] Soil dry density, unit: ;

[0019] The drying quality of soil collected in the field, unit: ;

[0020] Total volume of soil collected in the field, unit: ;

[0021] Provide soil drying density formula ;

[0022] Will Substitute In this context, the formula for calculating soil porosity is obtained.

[0023] Optionally, in some embodiments of the application, the formula for calculating the saturated permeability coefficient is as follows:

[0024] ;

[0025] In the formula:

[0026] ρw is the density of pure water, with a value of 1000, unit: ;

[0027] g is the acceleration of gravity, with a value of 9.81, unit: ;

[0028] θ is the soil porosity, unit: dimensionless;

[0029] D is the effective particle size of the soil, unit: ;

[0030] ηw is the dynamic viscosity coefficient of pure water, with a value of 0.001, unit: .

[0031] Optionally, in some embodiments of the application, the formula for calculating the saturated permeability coefficient is obtained as follows:

[0032] The continuity equation of the solid phase is provided:

[0033] ;

[0034] The energy conservation equation of the solid phase is:

[0035] ;

[0036] The continuity equation of the liquid phase is:

[0037] ;

[0038] The energy conservation equation of the liquid phase is:

[0039] ;

[0040] In the formula:

[0041] ρs is the inherent density of soil particles, with a value of 2650, unit: ;

[0042] ρw is the density of pure water, with a value of 1000, unit: ;

[0043] porosity of soil, dimensionless;

[0044] velocity field of soil particle movement, unit: ;

[0045] velocity field of water flow movement, unit: ;

[0046] macroscopic stress tensor of soil, unit: ;

[0047] stress tensor of water flow, unit: ;

[0048] time, unit: ;

[0049] divergence operator;

[0050] drag coefficient, unit: , and the calculation formula is as follows:

[0051] ;

[0052] In the formula:

[0053] dynamic viscosity coefficient of clear water, taking 0.001, unit: ;

[0054] effective particle size of soil, unit: ;

[0055] porosity of soil, dimensionless;

[0056] stress tensor of clear water flow is:

[0057] ;

[0058] In the formula:

[0059] velocity field of water flow movement, unit: ;

[0060] pressure field of water flow, unit: ;

[0061] is the kinematic viscosity of water, value is 0.001, unit: ;

[0062] is the unit tensor, unit: dimensionless;

[0063] is the gradient operator;

[0064] is the transpose operator;

[0065] Assume that the soil remains in a static state when the water flow moves in the soil, at this time the flow velocity field of the soil particle movement ; the soil porosity is assumed to remain unchanged, the formula and the formula are naturally satisfied; again assume that the water flow in the soil is a constant uniform flow state, the term and the term in the formula are both zero, at this time the left side of the equation is equal to zero, ignoring the dynamic viscosity of the water flow, i.e. the term of is zero; finally, bring and into to obtain:

[0066] ,

[0067] Simplify to get: ;

[0068] In the formula:

[0069] is the flow velocity field of the water flow movement, unit: ;

[0070] is the pressure field of the water flow, unit: ;

[0071] is the density of water, value is 1000, unit: ;

[0072] is the acceleration of gravity, value is 9.81, unit: ;

[0073] is the soil porosity, unit: dimensionless;

[0074] is the dynamic viscosity coefficient of clear water, taking a value of 0.001, and the unit is: ;

[0075] is the position head of water flow in soil, and the unit is: ;

[0076] is the saturated permeability coefficient of soil;

[0077] is the gradient operator;

[0078] so as to obtain the saturated permeability coefficient is

[0079] ;

[0080] In the formula:

[0081] is the density of clear water, taking a value of 1000, and the unit is: ;

[0082] is the acceleration of gravity, taking a value of 9.81, and the unit is: ;

[0083] is the soil porosity, and the unit is dimensionless;

[0084] is the effective particle size of soil, and the unit is: ;

[0085] is the dynamic viscosity coefficient of clear water, taking a value of 0.001, and the unit is: .

[0086] Correspondingly, the embodiment of the application also provides a soil saturated permeability coefficient obtaining system, comprising:

[0087] an effective particle size of soil obtaining module, which calculates the effective particle size of soil based on particle size distribution data;

[0088] a soil porosity obtaining module, which obtains the soil porosity;

[0089] a saturated permeability coefficient obtaining module, which obtains the saturated permeability coefficient based on solid-liquid two-phase fluid mechanics equations.

[0090] Correspondingly, the embodiment of the present application also provides a computer device, comprising a storage and a processor, wherein the storage stores a computer program, and the computer program is executed by the processor to make the processor execute the steps of the above method.

[0091] Correspondingly, the embodiment of the present application also provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to make the processor execute the steps of the above method.

[0092] In summary, due to the adoption of the above technical solutions, the beneficial effects of the present application are:

[0093] Based on the theory of fluid mechanics, the present application obtains parameters required for numerical calculation, including soil porosity, particle size distribution characteristics, water flow density, water flow dynamic viscosity, gravity acceleration, etc., and deduces the saturated permeability coefficient by comprehensively considering the soil porosity, particle size distribution characteristics and other parameters, so as to solve the problems of single parameter representation of traditional empirical formula, non-strict compliance with the theory of fluid mechanics and difficulty in representing nonlinear seepage characteristics under high hydraulic gradient. The method for obtaining the saturated permeability coefficient of the present application quantitatively considers the comprehensive influence of multiple parameters such as porosity, effective particle size reflecting particle distribution and water flow density, thereby breaking through the single parameter dependency of traditional empirical formula; meanwhile, the formula is derived based on the theory of fluid mechanics, thereby breaking through the simplification defect of theoretical model, significantly improving the calculation efficiency on the basis of ensuring the theoretical rigor. The saturated permeability coefficient calculation formula of the present application considers the porosity and the effective particle size based on the particle size distribution, and the two parameters exist in the form of square, which has an amplification effect on the structural differences of heterogeneous soil (such as clay soil and organic soil), and reflects the water flow path difference based on the theory of fluid mechanics, thereby accurately quantifying the seepage characteristic differences, and providing high-precision parameter support for slope anti-seepage design optimization, groundwater pollution prevention and control and geological disaster risk assessment.

[0094] Compared with the limitation of traditional formula relying on single particle size parameter, the present application selects the porosity and the effective particle size based on the particle size distribution to represent the pore structure and particle distribution respectively, and quantifies the synergistic effect of the two parameters in the form of square and combination of numerator and denominator in the saturated permeability coefficient formula. The saturated permeability coefficient formula can capture the key differences of complex soil, is sensitive to parameter differences, and has universal theory, thereby being applicable to various types of complex soil, significantly improving the efficiency and reliability of engineering anti-seepage design, groundwater simulation and geological disaster prevention, and having wide engineering application value. BRIEF DESCRIPTION OF DRAWINGS

[0095] Figure 1 The flow chart of the method for obtaining the saturated permeability coefficient of soil body of the present application is shown in the figure;

[0096] Figure 2 The particle size distribution curve of a certain soil sample of the present application is shown in the figure. DETAILED DESCRIPTION

[0097] In order to make the objectives, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and should not be used to limit the present application.

[0098] The technical solutions of the present application are as follows:

[0099] Referring to Figure 1 The embodiments of the present application provide a method for obtaining a saturated permeability coefficient of soil, comprising:

[0100] S01, calculating an effective particle size of the soil based on particle size distribution data;

[0101] S02, obtaining a soil porosity;

[0102] S03, obtaining the saturated permeability coefficient based on solid-liquid two-phase fluid mechanics equations.

[0103] In the S01,

[0104] In some embodiments, the particle size distribution data includes particle size values in the soil particle size distribution and mass fractions corresponding to the particle sizes in the soil particle size distribution .

[0105] It can be understood that the inherent densities of all soil particles are consistent.

[0106] Further, based on the particle size distribution data, the effective particle size is calculated according to the specific surface area of the soil particles, and the reciprocal of the effective particle size is equal to the sum of the specific surface areas of the soil particles, that is, the calculation formula of the effective particle size is as follows: ;

[0107] In the formula, is the effective particle size of the soil, the unit is ; is the particle size value in the soil particle size distribution, the unit is ; is the mass fraction corresponding to the particle size in the soil particle size distribution, the unit is dimensionless; is a group number, which sequentially represents the grouping of different particle sizes and corresponding mass fractions of the first group, the second group, …, and the nth group, the unit is dimensionless; n is the total number of groupings of the soil particle size distribution, representing the number of groupings obtained after the soil is divided according to the particle size range, the unit is dimensionless.

[0108] It can be understood that The group number is used to identify the soil particle size distribution group, sequentially referring to the group with different particle sizes and corresponding mass fractions of group 1, group 2... group n, enabling the calculation of each particle size distribution group. The unit is dimensionless. n is the total number of soil particle size distribution groups, representing the number of groups obtained after dividing the soil according to the particle size range, clarifying the total number of particle size distribution groups participating in the effective particle size calculation. The unit is dimensionless.

[0109] In S02:

[0110] In some embodiments, the formula for calculating soil porosity is as follows:

[0111] ;

[0112] In the formula: Soil dry density, unit: ; Soil porosity, unit: dimensionless; The intrinsic density of soil particles, in units of: Its value is 2650 ; The drying quality of soil collected in the field, unit: ; Total volume of soil collected in the field, unit: .

[0113] Furthermore, the method for obtaining soil porosity is as follows:

[0114] Provide soil dry density; the formula for calculating soil dry density is as follows:

[0115] ;

[0116] In the formula:

[0117] Soil dry density, unit: ;

[0118] The drying quality of soil collected in the field, unit: ;

[0119] Total volume of soil collected in the field, unit: ;

[0120] Provide soil drying density formula ;

[0121] Will Substitute In this context, the formula for calculating soil porosity is obtained.

[0122] Exemplarily, the soil dry density is calculated by weighing the soil original sample of a fixed volume after drying and collecting in the field.

[0123] In the S03, the formula for calculating the saturated permeability coefficient is as follows:

[0124] In some embodiments, the formula for calculating the saturated permeability coefficient is as follows:

[0125] ;

[0126] In the formula, the formula for calculating the saturated permeability coefficient is as follows:

[0127] is the density of pure water, taking the value of 1000, unit: ;

[0128] is the acceleration of gravity, taking the value of 9.81, unit: ;

[0129] is the soil porosity, unit: dimensionless;

[0130] is the effective particle size of the soil, unit: ;

[0131] is the dynamic viscosity coefficient of pure water, taking the value of 0.001, unit: .

[0132] Further, the formula for calculating the saturated permeability coefficient is obtained as follows:

[0133] The continuity equation of the solid phase is provided:

[0134] ;

[0135] The energy conservation equation of the solid phase is:

[0136] ;

[0137] The continuity equation of the liquid phase is:

[0138] ;

[0139] The energy conservation equation of the liquid phase is:

[0140] ;

[0141] In the formula, the formula for calculating the saturated permeability coefficient is as follows:

[0142] is the inherent density of soil particles, taking the value of 2650, unit: ;

[0143] For clear water density, take the value of 1000, unit: ;

[0144] For soil porosity, unit: dimensionless;

[0145] For the flow velocity field of soil particle movement, unit: ;

[0146] For the flow velocity field of water flow movement, unit: ;

[0147] For the soil macroscopic stress tensor, unit: ;

[0148] For the stress tensor of water flow, unit: ;

[0149] For time, unit: ;

[0150] For the divergence operator;

[0151] For the drag coefficient, unit: , the calculation formula is as follows:

[0152] ;

[0153] In the formula:

[0154] For the dynamic viscosity coefficient of clear water, take the value of 0.001, unit: ;

[0155] For the effective particle size of soil, unit: ;

[0156] For the soil porosity, unit: dimensionless;

[0157] The stress tensor of clear water flow is:

[0158] ;

[0159] In the formula:

[0160] For the flow velocity field of water flow movement, unit: ;

[0161] is the pressure field of water flow, unit: ;

[0162] is the dynamic viscosity coefficient of water, taking the value of 0.001, unit: ;

[0163] is the unit tensor, unit: dimensionless;

[0164] is the gradient operator;

[0165] is the transpose operator;

[0166] Assuming that the soil remains static when the water flow moves in the soil, the flow velocity field of soil particles at this time ; assuming that the soil porosity remains unchanged, the formula and the formula are naturally satisfied; assuming that the water flow in the soil is a constant and uniform flow state, the term and the term in the formula are both zero at this time, the left side of the equation is equal to zero, ignoring the dynamic viscosity of the water flow, i.e. the term of is zero; finally, bring and into to obtain:

[0167] ,

[0168] Simplify to get: ;

[0169] In the formula:

[0170] is the flow velocity field of water flow, unit: ;

[0171] is the pressure field of water flow, unit: ;

[0172] is the density of water, taking the value of 1000, unit: ;

[0173] is the acceleration of gravity, taking the value of 9.81, unit: ;

[0174] is the soil porosity, unit: dimensionless;

[0175] is the effective particle size of the soil, unit: ;

[0176] is the clear water dynamic viscosity coefficient, taking a value of 0.001, unit: ;

[0177] is the position water head of the water flow in the soil, unit: ;

[0178] is the soil saturated permeability coefficient;

[0179] is the gradient operator;

[0180] so as to obtain the saturated permeability coefficient is

[0181] ;

[0182] In the formula:

[0183] is the clear water density, taking a value of 1000, unit: ;

[0184] is the gravity acceleration, taking a value of 9.81, unit: ;

[0185] is the soil porosity, unit: dimensionless;

[0186] is the effective particle size of the soil, unit: ;

[0187] is the clear water dynamic viscosity coefficient, taking a value of 0.001, unit: .

[0188] In a second aspect, the embodiments of the present application provide a soil saturated permeability coefficient obtaining system, comprising:

[0189] an effective particle size of soil obtaining module, which calculates the effective particle size of the soil based on the particle size distribution data;

[0190] a soil porosity obtaining module, which obtains the soil porosity;

[0191] The saturated permeability coefficient obtaining module obtains the saturated permeability coefficient based on solid-liquid two-phase fluid mechanics equations.

[0192] In the effective particle size obtaining module of the soil:

[0193] In some embodiments, the particle size data in the particle size distribution of the soil includes particle size values in the particle size distribution of the soil and mass fractions corresponding to the particle sizes in the particle size distribution of the soil .

[0194] It can be understood that the inherent densities of all soil particles are consistent.

[0195] Further, based on the particle size data, the effective particle size is calculated according to the specific surface area of the soil particles, and the reciprocal of the effective particle size is equal to the sum of the specific surface areas of the soil particles, that is, the calculation formula of the effective particle size is as follows: ;

[0196] In the formula: is the effective particle size of the soil, the unit is: ; is the particle size value in the particle size distribution of the soil, the unit is: ; is the mass fraction corresponding to the particle size in the particle size distribution of the soil, the unit is dimensionless; is a group number, which sequentially refers to the grouping of different particle sizes and corresponding mass fractions of the first group, the second group, …, and the nth group, the unit is dimensionless; n is the total number of groupings of the particle size distribution of the soil, which represents the number of groupings obtained after the soil is divided according to the particle size range, and the total amount of the particle size distribution groupings participating in the calculation of the effective particle size, the unit is dimensionless.

[0197] It can be understood that, is a group number, which is used to identify the particle size distribution groupings of the soil, sequentially refers to the grouping of different particle sizes and corresponding mass fractions of the first group, the second group, …, and the nth group, realizes the one-by-one calculation of each particle size distribution grouping, the unit is dimensionless; n is the total number of groupings of the particle size distribution of the soil, which represents the number of groupings obtained after the soil is divided according to the particle size range, and the total amount of the particle size distribution groupings participating in the calculation of the effective particle size, the unit is dimensionless.

[0198] In the soil porosity obtaining module:

[0199] In some embodiments, the calculation formula of the soil porosity is as follows:

[0200] ;

[0201] In the formula: is the dry density of the soil, the unit is: ; is the porosity of the soil, the unit is dimensionless; is the inherent density of the soil particles, the unit is The value is 2650 ; The oven-dried mass of the soil collected in the field, unit: ; The total volume of the soil collected in the field, unit: .

[0202] Further, the method for obtaining the soil porosity is as follows:

[0203] The soil dry density is provided, and the calculation formula of the soil dry density is as follows:

[0204] ;

[0205] In the formula:

[0206] The soil dry density, unit: ;

[0207] The oven-dried mass of the soil collected in the field, unit: ;

[0208] The total volume of the soil collected in the field, unit: ;

[0209] The soil oven-dried density formula is provided ;

[0210] Bring Into , the calculation formula of the soil porosity is obtained.

[0211] Exemplarily, the soil dry density is calculated by weighing the oven-dried original soil sample of a fixed volume collected in the field.

[0212] In the saturated permeability coefficient obtaining module:

[0213] In some embodiments, the calculation formula of the saturated permeability coefficient is as follows:

[0214] ;

[0215] In the formula:

[0216] The fresh water density, the value is 1000, unit: ;

[0217] The gravitational acceleration, the value is 9.81, unit: ;

[0218] Porosity of soil, unit: dimensionless;

[0219] Effective particle size of soil, unit: ;

[0220] Dynamic viscosity coefficient of clear water, value is 0.001, unit: .

[0221] Further, the calculation formula of the saturated permeability coefficient is obtained as follows:

[0222] The continuity equation of solid phase is provided:

[0223] ;

[0224] The energy conservation equation of solid phase is:

[0225] ;

[0226] The continuity equation of liquid phase is:

[0227] ;

[0228] The energy conservation equation of liquid phase is:

[0229] ;

[0230] In the formula:

[0231] Intrinsic density of soil particles, value is 2650, unit: ;

[0232] Density of clear water, value is 1000, unit: ;

[0233] Porosity of soil, unit: dimensionless;

[0234] Flow velocity field of soil particle movement, unit: ;

[0235] Flow velocity field of water flow movement, unit: ;

[0236] Macroscopic stress tensor of soil, unit: ;

[0237] Stress tensor of water flow, unit: ;

[0238] is the time, unit: ;

[0239] is the divergence operator;

[0240] is the drag coefficient, unit: , the calculation formula is as follows:

[0241] ;

[0242] In the formula:

[0243] is the dynamic viscosity coefficient of clear water, the value is 0.001, unit: ;

[0244] is the effective particle size of soil, which is calculated by formula (1), unit: ;

[0245] is the soil porosity, unit: dimensionless;

[0246] Stress tensor of clear water flow is:

[0247] ;

[0248] In the formula:

[0249] is the flow velocity field of water flow, unit: ;

[0250] is the pressure field of water flow, unit: ;

[0251] is the dynamic viscosity coefficient of clear water, the value is 0.001, unit: ;

[0252] is the unit tensor, unit: dimensionless;

[0253] is the gradient operator;

[0254] is the transpose operator;

[0255] Assuming that the water flow moves in the soil, the soil remains in a static state, at this time the flow velocity field of soil particles ; the latter assuming soil porosity ; the former keeping unchanged ; the latter assuming soil porosity ; the former keeping unchanged ; the latter assuming soil porosity ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil

[0256] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil

[0257] ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0258] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0259] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0260] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0261] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0262] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0263] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil

[0264] ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0265] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0266] ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil

[0267] ; the former keeping unchanged ; the latter assuming constant and uniform flow in the soil ; the former keeping unchanged

[0268] Thus, the saturated permeability coefficient is obtained For

[0269] ;

[0270] In the formula:

[0271] is the density of clean water, and the value is 1000, the unit is: ;

[0272] is the acceleration of gravity, and the value is 9.81, the unit is: ;

[0273] is the soil porosity, the unit is dimensionless;

[0274] is the effective particle size of the soil, the unit is: ;

[0275] is the dynamic viscosity coefficient of clean water, and the value is 0.001, the unit is: .

[0276] In a third aspect, the present application provides a computer device, comprising a storage and a processor, the storage stores a computer program, and the computer program is executed by the processor to make the processor execute the steps of the soil saturated permeability coefficient obtaining method as described above.

[0277] The computer device can be a desktop computer, a notebook computer, a palm computer, a cloud server, or the like. The computer device can interact with a user through a keyboard, a mouse, a remote controller, a touchpad, a voice control device, or the like.

[0278] The memory includes at least one type of readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card-type memory (e.g., an SD or D interface display memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, etc. In some embodiments, the memory can be an internal storage unit of the computer device, such as a hard disk or a memory of the computer device. In other embodiments, the memory can also be an external storage device of the computer device, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the computer device. Of course, the memory can also include both the internal storage unit and the external storage device of the computer device. In the present embodiment, the memory is commonly used to store an operating system and various application software installed on the computer device, such as a program code of the soil saturated permeability coefficient obtaining method, etc. In addition, the memory can also be used to temporarily store various data that have been output or will be output.

[0279] The processor in some embodiments can be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chip. The processor is commonly used to control the overall operation of the computer device. In the present embodiment, the processor is used to run the program code or process data stored in the memory, such as running the program code of the soil saturated permeability coefficient obtaining method.

[0280] In a fourth aspect, the present application provides a computer readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the steps of the soil saturated permeability coefficient obtaining method as described above.

[0281] In the fourth aspect, the computer readable storage medium stores an interface display program, which can be executed by at least one processor to cause the at least one processor to perform the steps of the soil saturated permeability coefficient obtaining method as described above.

[0282] Through the description of the above embodiments, those skilled in the art can clearly understand that the above-mentioned example method can be realized by means of software and the necessary general hardware platform, of course, it can also be realized by hardware, but in many cases, the former is a better embodiment. Based on such understanding, the technical solutions of the present application can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes a plurality of instructions for causing a terminal device (which can be a mobile phone, computer, server or network device, etc.) to execute the soil saturated permeability coefficient obtaining method described in the embodiments of the present application.

[0283] The application will be further described below in combination with application examples.

[0284] Application Example

[0285] The volume of the sampler is 0.001 V , the original soil sample is collected from the slope of a certain debris flow ditch in China, and the volume of the sample is 1.70 kg, and the dried mass of the soil sample is 1.65 kg after drying. M The density of the water flow is 1000 , , the inherent density of the soil particles is 2650 , ; The dynamic viscosity of the water flow is ; The gravitational acceleration (unit: ) is 9.81 ; Figure 2 The particle size distribution curve of the soil is used to obtain the effective particle size required for calculating the saturated permeability coefficient, and is obtained by sieving the sampled soil. The above data is substituted into the formula as follows:

[0286] 1. According to the particle size distribution curve of the original soil sample collected from the slope of a certain debris flow ditch in China (as shown in Figure 2 ), the effective particle size is calculated:

[0287] ;

[0288] 2. The soil porosity is calculated according to the volume of the sampled soil V and the dried mass of the soil M :

[0289] ;

[0290] 3. The effective particle size and the soil porosity calculated are brought into the calculation of the saturated permeability coefficient :​

[0291] .

[0292] The above descriptions are only the preferred embodiments of the present application, not intended to limit the present application, and any modification, equivalent replacement and improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. Method for obtaining the saturated permeability coefficient of a soil, characterized in that, The method comprises: calculating the effective grain diameter of the soil based on the grain size distribution data; obtaining the soil porosity; obtaining the saturated permeability coefficient based on the solid-liquid two-phase fluid mechanics equation; The calculation formula of the saturated permeability coefficient is as follows: ； In the formula: For clear water density, take the value of 1000, unit: ; g is the gravitational acceleration, having a value of 9.81, units: ; porosity of the soil, unit: dimensionless; effective diameter of soil, unit: ; Cf is the clear water dynamic viscosity coefficient, with a value of 0.001, unit: .

2. The method for obtaining the soil saturated permeability coefficient according to claim 1, characterized in that, The particle size data of the particle size distribution of the soil includes the particle size value and the mass fraction corresponding to the particle size in the particle size distribution of the soil .

3. The method for obtaining the soil saturated permeability coefficient according to claim 2, characterized in that, Based on the grain size distribution data, the effective grain diameter is calculated according to the specific surface area of soil particles, and the reciprocal of the effective grain diameter is equal to the sum of the specific surface area of soil particles, that is, the calculation formula of the effective grain diameter is as follows: ； In the formula: is the effective particle size of the soil, unit: ; is the particle size value in the soil particle size distribution, unit: ; is the mass fraction corresponding to the particle size in the soil particle size distribution, unit: dimensionless; is the grouping number, which sequentially represents the grouping of the first group, the second group, …, the nth group of different particle sizes and corresponding mass fractions, unit: dimensionless; n is the total number of groupings of the soil particle size distribution, representing the number of groupings obtained after the soil is divided according to the particle size range, unit: dimensionless.

4. The method for obtaining the soil saturated permeability coefficient according to claim 1, characterized in that, The calculation formula of the soil porosity is as follows: ； wherein: is the dry density of the soil, unit: ; is the porosity of the soil, unit: dimensionless; is the inherent density of the soil particles, unit: with a value of 2650 ; is the oven-dried mass of the field-collected soil, unit: ; is the total volume of the field-collected soil, unit: .

5. The method for obtaining the soil saturated permeability coefficient according to claim 4, characterized in that, The method for obtaining the soil porosity is as follows: The soil dry density is provided, and the calculation formula of the soil dry density is as follows: ； In the formula: soil dry density, unit: ; To dry the quality of the soil collected in the field, unit: ; Total volume of soil collected in the field, unit: ; Providing soil dry density formula ; Bringing into , the formula for calculating soil porosity is obtained.

6. The method of obtaining the soil saturated hydraulic conductivity according to claim 1, characterized in that, The calculation formula of the saturated permeability coefficient is obtained as follows: The continuity equation of the solid phase is provided: ； The energy conservation equation of the solid phase is provided: ； The continuity equation of the liquid phase is provided: ； The energy conservation equation of the liquid phase is provided: ； In the formula: Soil particle inherent density, valued at 2650, units: ; For clear water density, take the value of 1000, unit: ; for soil porosity, unit: dimensionless; is the flow velocity field of the soil particle motion, units: ; is the velocity field of the water flow movement, with units of: ; is the soil macroscopic stress tensor, units: ; is the stress tensor of the water flow, units: ; Time, units: ; is the divergence operator; For the drag coefficient, units: The formula is as follows: ； In the formula: Cw is the clear water dynamic viscosity coefficient, with a value of 0.001, unit: ; effective diameter of soil, unit: ; for soil porosity, unit: dimensionless; Stress tensor of a clear water flow is: ； In the formula: is the velocity field of the water flow movement, with units of: ; is the pressure field of the water flow, in units of: ; Cw is the clean water dynamic viscosity coefficient, with a value of 0.001, unit: ; unit tensor, unit: dimensionless; is the gradient operator; is the transpose operator; Assume that the water flow moves in the soil while the soil keeps static, the flow velocity field of soil particles at this time ; the latter assumes that the soil porosity remains unchanged, the formula and the formula naturally satisfy; again assume that the water flow in the soil is a constant uniform flow state, the term and the term in the formula are zero, at this time the left side of the equation is equal to zero, ignore the water flow dynamic viscosity, that is, the term of the formula is zero; finally, bring and into to get: ， Simplifying gives: ; In the formula: is the velocity field of the water flow movement, with units of: ; is the pressure field of the water flow, in units of: ; For clear water density, take the value of 1000, unit: ; g is the gravitational acceleration, having a value of 9.81, units: ; for soil porosity, unit: dimensionless; effective diameter of soil, unit: ; Cw is the clear water dynamic viscosity coefficient, with a value of 0.001, unit: ; Water head for the position of the water flow in the soil, unit: ; Ks is the saturated soil permeability coefficient; is the gradient operator; Thus, the saturated permeability coefficient is ； In the formula: For clear water density, take the value of 1000, unit: ; g is the gravitational acceleration, having a value of 9.81, units: ; Porosity of the soil, unit: dimensionless; effective diameter of soil, unit: ; Cf is the clear water dynamic viscosity coefficient, with a value of 0.001, unit: .

7. A system for obtaining the saturated hydraulic conductivity of a soil, characterized in that, The method comprises: The effective grain diameter obtaining module calculates the effective grain diameter of the soil based on the grain size distribution data; The soil porosity obtaining module obtains the soil porosity; The saturated permeability coefficient obtaining module obtains the saturated permeability coefficient based on the solid-liquid two-phase fluid mechanics equation; The calculation formula of the saturated permeability coefficient is as follows: ； In the formula: For clear water density, take the value of 1000, unit: ; g is the gravitational acceleration, having a value of 9.81, units: ; for soil porosity, unit: dimensionless; effective diameter of soil, unit: ; Cf is the clear water dynamic viscosity coefficient, with a value of 0.001, unit: .

8. Computer device, characterized in that The computer program is stored in the storage and executed by the processor, so that the processor executes the steps of the method according to any one of claims 1-6.

9. A computer readable storage medium, characterized in that, The computer program is stored in the storage and executed by the processor, so that the processor executes the steps of the method according to any one of claims 1-6.

Images