Method for analyzing deformation of soil pore structure under uniaxial compaction

Abstract

The invention belongs to the technical field of soil science, and discloses a soil pore structure deformation analysis method under uniaxial compaction, which comprises the following steps: collecting a soil sample; performing CT scanning on the soil sample to obtain a first image slice; performing uniaxial compaction on the soil sample, and performing CT scanning on the compacted soil sample to obtain a second image slice; performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing the change of soil pore characteristic parameters before and after compaction; according to the first image slice and the second image slice, a pre-compaction binary image and a post-compaction binary image of the biological pore are obtained; and calculating the change of the characteristic parameters of the biological pores before and after compaction according to the binary image before compaction and the binary image after compaction, and evaluating the stability of the biological pores in the compaction process. And quantitative analysis can be carried out on the deformation of the soil pore structure under uniaxial compaction.

Description

Technical Field

[0001] The invention belongs to the technical field of soil science, and in particular relates to a method for analyzing soil pore structure deformation under uniaxial compaction. Background Art

[0002] The essence of soil compaction is the rearrangement of soil particles, which leads to an increase in soil bulk density and a decrease in porosity. The deformation of the pore structure caused by soil compaction further hinders water vapor transmission and root growth, and affects processes such as redox and greenhouse gas emissions in the soil. Soil compaction has become one of the main forms of soil degradation in the world. Soil compaction first affects the macropore structure of the soil, and the degree of pore deformation is related to factors such as applied stress, initial pore structure, and soil moisture. Among them, for soil pores with different initial structures, their deformation laws during the compaction process are different. For example, biological pores are relatively stable during the compaction process, and tubular biological pores can maintain their morphological stability under stress, maintain the ventilation and drainage functions of the soil, and have high application potential in alleviating compaction risks.

[0003] Traditional methods have difficulty in quantifying the deformation of soil pores during compaction and the stabilizing effect of biological pores. X-ray computed tomography has been widely used to quantify soil pore structure. However, the deformation of the microscopic pore structure of farmland soil under uniaxial compaction has not been systematically analyzed, especially the stability of tubular biological pores during compaction.

[0004] Therefore, it is not possible to quantitatively analyze the deformation of soil pore structure under uniaxial compaction. Summary of the invention

[0005] The purpose of the present invention is to provide a method for analyzing soil pore structure deformation under uniaxial compaction, which can quantitatively analyze the soil pore structure deformation under uniaxial compaction.

[0006] The first aspect of the present invention discloses a method for analyzing soil pore structure deformation under uniaxial compaction, comprising:

[0007] Collect soil samples;

[0008] Performing CT scanning on the soil sample to obtain a first image slice;

[0009] uniaxially compacting the soil sample, and performing CT scanning on the compacted soil sample to obtain a second image slice;

[0010] Performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing changes in soil pore characteristic parameters before and after compaction based on the binary image;

[0011] Acquire a pre-compacted binary image and a post-compacted binary image of biological pores according to the first image slice and the second image slice;

[0012] The changes of biological pore characteristic parameters before and after compaction are calculated based on the binary image before compaction and the binary image after compaction, and the stability of the biological pores during the compaction process is evaluated.

[0013] In some embodiments, before performing CT scanning on the soil sample, the method further comprises:

[0014] Wrap filter paper at the bottom of the soil sample and place it in a sand box for saturation. After the soil sample is saturated, adjust the suction value of the sand box to a preset value so that the soil sample is at the same suction.

[0015] In some embodiments, performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing changes in soil pore characteristic parameters before and after compaction based on the binary image, includes:

[0016] performing image preprocessing on the first image slice and the second image slice;

[0017] Selecting a central area in the preprocessed first image slice and the second image slice as a region of interest to obtain a target image;

[0018] Performing threshold segmentation on the target image using a single threshold segmentation method to obtain a binary image;

[0019] Based on the binary image, the changes of soil pore characteristic parameters before and after compaction are calculated, and the soil pore characteristic parameters include soil macroporosity, pore connectivity and soil pore size distribution.

[0020] In some embodiments, before selecting the central area in the preprocessed first image slice and the second image slice as the region of interest, the following preprocessing steps are also included:

[0021] The brightness and contrast of the first image slice and the second image slice are adjusted; and the first image slice and the second image slice are subjected to noise reduction by using median filtering.

[0022] In some embodiments, the calculation of changes in pore characteristic parameters before and after compaction based on the binary image includes:

[0023] Based on the binary image, the volume fraction module in Avizo software was used to calculate the change of soil macroporosity before and after compaction. The connectivity module in Avizo software was used to calculate the change of pore connectivity before and after compaction. The thickness module in Avizo software was used to calculate the change of soil pore size distribution before and after compaction.

[0024] In some embodiments, the soil sample is subjected to uniaxial compaction, and the compacted soil sample is subjected to CT scanning to obtain the second image slice, further comprising:

[0025] Drying and weighing the compacted soil sample to obtain the mass of the soil sample;

[0026] The changes in soil bulk density and soil porosity before and after compaction are calculated based on the sample mass and the deformation data of the soil sample before and after compaction.

[0027] In some embodiments, acquiring a pre-compacted binarized image and a post-compacted binarized image of biological pores according to the first image slice and the second image slice comprises:

[0028] Determining biological pores that meet requirements according to the first image slice and the second image slice;

[0029] The biological pores are separated from the first image slice and the second image slice by using the region growing method in the VG studio software, and the binary image before compaction and the binary image after compaction are obtained respectively.

[0030] In some embodiments, determining the biological pores that meet the requirements according to the first image slice and the second image slice includes:

[0031] In the first image slice, biological pores meeting the requirements are found by visual inspection according to the morphological characteristics of the biological pores, and corresponding biological pores are found in the second image slice according to the positions and trends of the biological pores in the first image slice.

[0032] In some embodiments, the biological pore characteristic parameters are volume, length and pore diameter of the biological pores, and the step of calculating the changes of the biological pore characteristic parameters before and after compaction based on the binary image before compaction and the binary image after compaction to evaluate the stability of the biological pores during the compaction process includes:

[0033] For the binarized image before compaction and the binarized image after compaction of the same biological pore, the three-dimensional structure of the biological pore before and after compaction is compared through the volume rendering module of the Avizo software, the change of the porosity of the biological pore before and after compaction is calculated through the volume fraction module of the Avizo software, the change of the pore size of the biological pore before and after compaction is calculated according to the three-dimensional structure through the thickness module of the Avizo software, the biological pore skeleton is generated through the skeleton module of the Avizo software, and the change of the biological pore length before and after compaction is calculated according to the biological pore skeleton.

[0034] In some embodiments, after generating the biological pore skeleton by the skeleton module of the Avizo software, the method further comprises:

[0035] The biological pores with an angle between the biological pore skeleton and the z-axis less than 30° are classified as longitudinal biological pores, and the biological pores with an angle between the biological pore skeleton and the z-axis greater than 60° are classified as transverse biological pores. The effect of the biological pore angle on pore deformation is analyzed.

[0036] The beneficial effect of the present invention is that, by performing CT scanning on soil samples before and after compaction, a three-dimensional image of the soil structure is obtained non-destructively, and the changes in soil macroporosity, pore connectivity and pore size distribution of the same soil sample before and after uniaxial compaction are calculated, and the influence of uniaxial compaction on the soil pore structure is quantified. By extracting the naturally formed biological pores in the soil, the changes in characteristic parameters related to the biological pores before and after compaction are calculated based on the binary images of the biological pores before and after compaction, the stability and deformation law of the biological pores in uniaxial compaction are quantitatively evaluated, and the quantitative analysis of the deformation law of soil pores under uniaxial compaction is comprehensively realized. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The accompanying drawings herein show specific examples of the technical solutions described in the present invention, and together with the specific implementation methods, constitute a part of the specification, and are used to explain the technical solutions, principles and effects of the present invention.

[0038] Unless otherwise specified or defined, the same reference numerals in different drawings represent the same or similar technical features, and the same or similar technical features may also be represented by different reference numerals.

[0039] Figure 1 It is a flow chart of a method for analyzing soil pore structure deformation under uniaxial compaction disclosed in an embodiment of the present invention;

[0040] Figure 2 is a comparison diagram of soil samples before and after compaction according to an embodiment of the present invention;

[0041] Figure 3 is a three-dimensional image of pore deformation under uniaxial compaction according to an embodiment of the present invention;

[0042] Figure 4 is a specific flow chart of step S400 in an embodiment of the present invention;

[0043] Figure 5 Schematic diagram of biological pores before and after compaction in an embodiment of the present invention.

[0044] Figure 6 It is a three-dimensional image of biological pores and their skeleton deformation under uniaxial compression according to an embodiment of the present invention. DETAILED DESCRIPTION

[0045] Unless otherwise specified or defined, all technical and scientific terms used herein have the same meaning as those generally understood by those skilled in the art. In the case of combining the technical solution of the present invention with realistic scenarios, all technical and scientific terms used herein may also have meanings corresponding to the purpose of implementing the technical solution of the present invention. The "first, second..." used herein is only used to distinguish the names and does not represent a specific quantity or order. The term "and / or" used herein includes any and all combinations of one or more related listed items.

[0046] It should be noted that when a component is considered to be "fixed to" another component, it can be directly fixed to the other component or there can be a central component; when an component is considered to be "connected to" another component, it can be directly connected to the other component or there can be a central component at the same time; when an component is considered to be "installed on" another component, it can be directly installed on the other component or there can be a central component at the same time. When an component is considered to be "set on" another component, it can be directly set on the other component or there can be a central component at the same time.

[0047] Unless otherwise specified or defined, the "said" and "the" used in this document refer to the technical features or technical contents mentioned or described before the corresponding position, and the technical features or technical contents may be the same as or similar to the technical features or technical contents mentioned therein. In addition, the terms "including" and "having" and any variations thereof used in this document are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units that are not listed, or may optionally include other steps or units inherent to these processes, methods, products or devices.

[0048] The present invention proposes a method for analyzing soil pore structure deformation under uniaxial compaction, which can calculate the changes in soil macroporosity, pore connectivity and soil pore size distribution before and after compaction, and quantitatively analyze the stability of biological pores during the compaction process, thereby realizing the quantitative analysis of the soil pore deformation law under uniaxial compaction.

[0049] The embodiment of the present invention takes farmland soil as an example, and uses non-destructive CT scanning technology to compare the pore deformation of the original soil of the same farmland before and after uniaxial compaction; and extracts the naturally formed biological pores in the field soil to achieve quantitative analysis of the stability of the biological pores during the compaction process. It should be emphasized that the analysis method of the present invention is not limited to farmland soil, and soil pore structure deformation analysis can also be performed on other types of soil.

[0050] In order to facilitate the understanding of the present invention, specific embodiments of the present invention will be described in more detail below with reference to the accompanying drawings.

[0051] like Figure 1 As shown, the method comprises the following steps:

[0052] Step S100: collecting soil samples;

[0053] The soil sample in this embodiment is an original soil sample of the tillage layer of the farmland. The collection process is as follows: a PC tube with a height of 5.5 cm and an outer diameter of 5 cm with one end ground into a blade is used to collect the original soil sample of the tillage layer of the field, and the original soil sample is wrapped with plastic wrap to facilitate CT scanning of the soil sample.

[0054] Step S200: performing CT scanning on the soil sample to obtain a first image slice;

[0055] In order to make the soil samples have the same suction force during uniaxial compaction, so as to compare the deformation laws of different soil samples, before CT scanning the soil samples, this embodiment wraps filter paper at the bottom of the soil samples and places them in a sandbox for saturation. After the soil samples are saturated, the suction force of the sandbox is adjusted to a preset value so that the soil samples have the same suction force value. Then, the soil samples are scanned and reconstructed using industrial CT, and grayscale image slices in TIFF format are obtained after reconstruction, i.e., the first image slices.

[0056] Step S300: performing uniaxial compaction on the soil sample, and performing CT scanning on the compacted soil sample to obtain a second image slice;

[0057] A uniaxial compaction test is performed on the soil sample using a universal testing machine. After the uniaxial compaction test, a CT scan is performed on the compacted soil sample. After reconstruction, a second image slice is obtained.

[0058] Specifically, a uniaxial compaction test is performed on the soil using a universal testing machine: a load greater than the soil pre-compression stress is set on the universal testing machine, for example, 200 kPa and 400 kPa loads are set to compact the soil sample, the load is applied for 4 hours, and finally the stress is adjusted back to 1 kPa for 0.5 hours to allow the soil to rebound naturally. The compacted soil sample is then CT scanned, and after reconstruction, a grayscale image slice of the compacted soil sample in TIFF format is obtained, i.e., the second image slice.

[0059] In addition, in this embodiment, after the soil sample is uniaxially compacted, Figure 2As shown, the deformation of the soil after stress removal is also measured to obtain the volume of the compacted soil, and the changes in soil bulk density and soil porosity before and after compaction are calculated. Specifically, the height of the compacted soil sample is first measured, and then the soil sample is dried and weighed, and then the changes in soil bulk density and soil porosity before and after compaction are calculated based on the sample mass and the measurement data.

[0060] The formula for calculating the bulk density of soil before compaction is: The calculation formula for the bulk density of soil after compaction is: The calculation formula for soil porosity before compaction is: The calculation formula for soil porosity after compaction is: Where m is the mass of the soil sample, in g; H is the height of the soil sample, in cm; d is the soil settlement after stress removal, in cm; r is the radius of the soil sample, in cm; ρ s is the particle density of the soil sample, that is, the ratio of the mass of the soil particles to their volume, in g / cm 3 .

[0061] Step S400: performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing changes in soil pore characteristic parameters before and after compaction based on the binary image;

[0062] The three-dimensional image of pore deformation under uniaxial compression is shown in Figure 3 As shown, for the first image slice and the second image slice, the soil pores are selected by Image software and Avizo software, the pore three-dimensional structure is rendered, and the changes in soil macroporosity and soil pore size distribution before and after compaction are calculated.

[0063] like Figure 4 As shown, the specific steps of this embodiment include:

[0064] Step S410: performing image preprocessing on the first image slice and the second image slice;

[0065] The first image slice and the second image slice are imported into ImageJ software, the brightness and contrast of the first image slice and the second image slice are adjusted by ImageJ software (Image-adjust-Brightness / Contrast), and the first image slice and the second image slice are subjected to noise reduction by using median filtering (Process-Filters-Median) to improve the quality of the target image, thereby improving the accuracy of the analysis;

[0066] Step S420: selecting a central area in the preprocessed first image slice and the second image slice as a region of interest to obtain a target image;

[0067] By selecting the central area as the region of interest (ROI) for further image processing, the influence of the PC tube wall during uniaxial compression can be reduced. Specifically, a certain range is selected at the center position of the first image slice and the center position of the second image slice to obtain the target image.

[0068] Step S430: using a single threshold segmentation method to perform threshold segmentation on the target image to obtain a binary image;

[0069] Specifically, the single threshold segmentation method in ImageJ software is used to perform threshold segmentation on the target image to obtain a binary image of the target image.

[0070] Step S440: Based on the binary image, calculate the changes in soil pore characteristic parameters before and after compaction, where the soil pore characteristic parameters include soil macroporosity, pore connectivity and soil pore size distribution.

[0071] Based on the binary image, the volume fraction module (Volume Fraction) in the Avizo software was used to calculate the change of soil macroporosity before and after compaction. The connectivity module (Axis Connecticity) in the Avizo software was used to calculate the change of pore connectivity before and after compaction. The thickness module (Thickness) in the Avizo software was used to calculate the change of soil pore size distribution before and after compaction.

[0072] Step S500: acquiring a binary image before compaction and a binary image after compaction of biological pores according to the first image slice and the second image slice;

[0073] Firstly, biological pores meeting the requirements are determined in the first image slice and the second image slice, for example, relatively complete, regularly shaped and other representative biological pores. Then, images corresponding to the biological pores are separated in the first image slice and the second image slice to obtain binary images of the biological pores before and after compaction.

[0074] In this embodiment, the specific process of determining the biological pores that meet the requirements in the first image slice and the second image slice is as follows: in the first image slice, the biological pores that meet the requirements are found by visual method based on the morphological characteristics of the biological pores (such as tubular shape and relatively continuous). For example, observe the soil sample and visually identify the relatively complete and regularly shaped biological pores in the grayscale image. According to the position and trend of the biological pores in the first image slice, the corresponding biological pores are found in the second image slice, that is, the same biological pores after compaction. For the biological pores before and after compaction, please refer to Figure 5 .

[0075] Then, based on the first image slice and the second image slice, the region growing method in VG studio is used, and the manual selection function is assisted to separate the representative soil biological pores from the first image slice and the second image slice, and obtain the binary image before compaction and the binary image after compaction, respectively. The region growing method of the VG studio tool can automatically identify the boundary between the target pore and the surrounding soil matrix, and can also be operated interactively by human and computer to avoid selecting the target biological pore and its adjacent non-biological pores at the same time, and finally obtain the binary image of the biological pore.

[0076] Step S600: Calculate the changes of biological pore characteristic parameters before and after compaction based on the binary image before compaction and the binary image after compaction, and evaluate the stability of biological pores during the compaction process.

[0077] The characteristic parameters of biological pores are usually the volume, length, pore diameter and angle of biological pores.

[0078] The three-dimensional image of biological pore deformation under uniaxial compression is shown in Figure 2. Figure 6 As shown in the figure, for the same biological pore before and after compaction, the changes in volume, length and pore diameter of the biological pore before and after compaction are calculated based on the binary images before and after compaction, respectively, to evaluate the stability of the original biological pores in the field.

[0079] Specifically, for the same biological pore in the binary image before and after compaction, the volume rendering module (Volume Rendering) of the Avizo software is used to generate and visualize the three-dimensional structure of the biological pore before and after compaction. The volume fraction module (Volume Fraction) of the Avizo software is used to calculate the change of the biological porosity before and after compaction according to the three-dimensional structure. The thickness module (Thickness) of the Avizo software is used to calculate the pore size distribution of the biological pore according to the three-dimensional structure, and the pore size weighted average is used as the pore size of the biological pore. The change of the pore size of the biological pore before and after compaction is calculated. The skeleton module (Skeleton) of the Avizo software is used to obtain the biological pore skeleton, and the change of the biological pore length before and after compaction is calculated based on the biological pore skeleton.

[0080] In this embodiment, after obtaining the biological pore skeleton through the skeleton module of the Avizo software, the biological pores are also classified to analyze the response law of the biological pore angle to uniaxial compaction. Specifically, considering that the pore stress deformation is related to the pore angle, the biological pores whose angle between the biological pore skeleton and the z-axis is less than 30° are classified as longitudinal biological pores, and the biological pores whose angle between the biological pore skeleton and the z-axis is greater than 60° are classified as transverse biological pores.

[0081] In some embodiments, the percentage of the biological pore characteristic parameter change before and after compaction is used instead of the absolute value change of the biological pore characteristic parameter to analyze the change pattern of the biological pore characteristics before and after compaction.

[0082] In summary, this embodiment compacts the soil sample under a specific matrix potential, and uses a universal testing machine to measure the settlement of the soil sample. The three-dimensional image of the soil structure is non-destructively obtained through CT scanning technology, and the soil particles and pores are separated through further image processing, and the changes in the soil pore structure parameters of the same soil sample before and after uniaxial compaction are calculated; biological pores at different angles are selected based on the visual method and regional growth method, and the stability and deformation law of biological pores in uniaxial compaction are quantitatively evaluated. Thus, through non-destructive CT scanning technology, the effect of uniaxial compaction on the porosity of the original soil in farmland is analyzed; the naturally formed biological pores in the farmland soil are extracted, and the quantitative analysis of the stability of the biological pores during the compaction process is realized, in order to provide a scientific basis for alleviating mechanical compaction in the field.

[0083] The purpose of the above embodiments is to exemplarily reproduce and deduce the technical solution of the present invention, and to fully describe the technical solution, purpose and effect of the present invention. Its purpose is to make the public understand the disclosed content of the present invention more thoroughly and comprehensively, and it does not limit the scope of protection of the present invention.

[0084] The above embodiments are not exhaustive enumerations of the present invention, and there may be multiple other implementations not listed. Any replacement and improvement made without violating the concept of the present invention shall fall within the protection scope of the present invention.

Claims

1. A method for analyzing soil pore structure deformation under uniaxial compaction, characterized in that: include: Collect soil samples; Performing CT scanning on the soil sample to obtain a first image slice; uniaxially compacting the soil sample, and performing CT scanning on the compacted soil sample to obtain a second image slice; Performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing changes in soil pore characteristic parameters before and after compaction based on the binary image; Acquire a pre-compacted binary image and a post-compacted binary image of biological pores according to the first image slice and the second image slice; The changes of biological pore characteristic parameters before and after compaction are calculated based on the binary image before compaction and the binary image after compaction, and the stability of the biological pores during the compaction process is evaluated.

2. The soil pore structure deformation analysis method under uniaxial compaction according to claim 1, characterized in that: Before performing CT scanning on the soil sample, the method further includes: Wrap filter paper at the bottom of the soil sample and place it in a sand box for saturation. After the soil sample is saturated, adjust the suction value of the sand box to a preset value so that the soil sample is at the same suction.

3. The soil pore structure deformation analysis method under uniaxial compaction according to claim 1, characterized in that: Performing image processing on the first image slice and the second image slice to obtain a binary image, and analyzing changes in soil pore characteristic parameters before and after compaction based on the binary image, including: performing image preprocessing on the first image slice and the second image slice; Selecting a central area in the preprocessed first image slice and the second image slice as a region of interest to obtain a target image; Performing threshold segmentation on the target image using a single threshold segmentation method to obtain a binary image; Based on the binary image, the changes of soil pore characteristic parameters before and after compaction are calculated, and the soil pore characteristic parameters include soil macroporosity, pore connectivity and soil pore size distribution.

4. The soil pore structure deformation analysis method under uniaxial compaction according to claim 3, characterized in that: Before selecting the central area as the region of interest in the preprocessed first image slice and the second image slice, the following preprocessing steps are also included: adjusting brightness and contrast of the first image slice and the second image slice; Median filtering is used to reduce noise on the first image slice and the second image slice.

5. The soil pore structure deformation analysis method under uniaxial compaction according to claim 3, characterized in that: The calculation of the change of soil pore characteristic parameters before and after compaction based on the binary image includes: Based on the binary image, the volume fraction module in Avizo software was used to calculate the change of soil macroporosity before and after compaction. The connectivity module in Avizo software was used to calculate the change of pore connectivity before and after compaction. The thickness module in Avizo software was used to calculate the change of soil pore size distribution before and after compaction.

6. The soil pore structure deformation analysis method under uniaxial compaction according to claim 1, characterized in that: The method further comprises: performing uniaxial compaction on the soil sample, performing CT scanning on the compacted soil sample, and obtaining a second image slice, and further comprising: drying and weighing the compacted soil sample to obtain the mass of the soil sample; The changes in soil bulk density and soil porosity before and after compaction are calculated based on the mass of the soil sample and the deformation data of the soil sample before and after compaction.

7. The soil pore structure deformation analysis method under uniaxial compaction according to claim 1, characterized in that: The step of obtaining a pre-compacted binary image and a post-compacted binary image of biological pores according to the first image slice and the second image slice comprises: Determining biological pores that meet requirements according to the first image slice and the second image slice; The biological pores are separated from the first image slice and the second image slice by using the region growing method in the VG studio software, and the binary image before compaction and the binary image after compaction are obtained respectively.

8. The soil pore structure deformation analysis method under uniaxial compaction according to claim 7, characterized in that: The step of determining the biological pores meeting the requirements according to the first image slice and the second image slice comprises: In the first image slice, biological pores meeting the requirements are found by visual inspection according to the morphological characteristics of the biological pores, and corresponding biological pores are found in the second image slice according to the positions and trends of the biological pores in the first image slice.

9. The method for analyzing soil pore structure deformation under uniaxial compaction according to any one of claims 1 to 8, characterized in that: The biological pore characteristic parameters are the volume, length and pore diameter of the biological pores, and the calculation of the changes of the biological pore characteristic parameters before and after compaction based on the binary image before compaction and the binary image after compaction to evaluate the stability of the biological pores during the compaction process includes: For the binarized image before compaction and the binarized image after compaction of the same biological pore, the three-dimensional structure of the biological pore before and after compaction is compared through the volume rendering module of the Avizo software, the change of the porosity of the biological pore before and after compaction is calculated through the volume fraction module of the Avizo software, the change of the pore size of the biological pore before and after compaction is calculated according to the three-dimensional structure through the thickness module of the Avizo software, the biological pore skeleton is generated through the skeleton module of the Avizo software, and the change of the biological pore length before and after compaction is calculated according to the biological pore skeleton.

10. The soil pore structure deformation analysis method under uniaxial compaction according to claim 9, characterized in that: After the biological pore skeleton is generated through the skeleton module of Avizo software, it also includes: The biological pores with an angle between the biological pore skeleton and the z-axis less than 30° are classified as longitudinal biological pores, and the biological pores with an angle between the biological pore skeleton and the z-axis greater than 60° are classified as transverse biological pores. The effect of the biological pore angle on pore deformation is analyzed.

Images