Soil pore water quantification method, device, electronic equipment and storage medium

Abstract

The application relates to the technical field of soil structure detection, in particular to a soil pore water quantification method and device, an electronic device and a storage medium. The method comprises the following steps: filling a preset gas into undisturbed soil to obtain undisturbed soil in a gas-liquid equilibrium state; determining a three-dimensional image of a gas phase, a liquid phase and a soil matrix phase in the undisturbed soil; and determining water distribution and water content of the undisturbed soil pores according to the three-dimensional image. According to the technical scheme, the water distribution recognition stability and accuracy of the soil pores can be effectively ensured without changing the water structure and natural state of the soil pores.

Description

Technical Field

[0001] This application relates to the field of soil structure detection technology, specifically to a method, apparatus, electronic device, and storage medium for quantifying soil pore moisture. Background Technology

[0002] Computed tomography (CT) combined with 3D visualization software can be used for non-destructive observation of soil pore structure. However, because the X-ray attenuation coefficients of water and soil pores are similar and their grayscale ranges overlap significantly, the accuracy of pore identification in water-bearing soils is low.

[0003] Currently, liquid tracers such as iodides are used to detect soil structure, but this liquid tracer method changes the soil's moisture distribution and surface tension, disrupts its natural state, and cannot be repeated. Summary of the Invention

[0004] In view of the above problems, the embodiments of this application provide a method, apparatus, electronic device and storage medium for quantifying soil pore moisture, which can effectively ensure the stability and accuracy of soil pore moisture distribution identification without changing the moisture structure and natural state within the soil pores.

[0005] According to one aspect of the embodiments of this application, a method for quantifying soil pore moisture is provided, the method comprising: The soil to be tested is obtained by filling the undisturbed soil with a preset gas to obtain the soil under gas-liquid equilibrium state; Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

[0006] In some embodiments of this application, the step of filling the undisturbed soil with a preset gas to obtain the soil to be tested in a gas-liquid equilibrium state includes: The undisturbed soil is placed in a sealed cavity with an air inlet valve; wherein the sealed cavity is a PVC pipe or an acrylic pipe. The preset gas is controlled to be introduced into the sealed cavity through the air inlet valve of the sealed cavity to obtain the soil to be tested in a gas-liquid equilibrium state.

[0007] In some embodiments of this application, the step of filling the undisturbed soil with a preset gas to obtain the soil to be tested in a gas-liquid equilibrium state further includes: According to preset filling conditions, the original soil is filled with preset gas; wherein, the preset filling conditions include maintaining preset pressure, preset filling rate, preset temperature and preset humidity, and continuously filling the preset gas concentration for a first preset duration; The original soil filled with the preset gas is left to stand for a second preset time to obtain the soil to be tested in a gas-liquid equilibrium state.

[0008] In some embodiments of this application, the preset pressure range is 0.05-0.20 MPa, the preset filling rate range is 30-50 mL / min, the preset temperature range is 15-30 °C, the preset humidity range is 0-10%, and the preset concentration range is 99%-100%.

[0009] In some embodiments of this application, the volume fraction of carbon dioxide in the preset gas is greater than or equal to 99%.

[0010] In some embodiments of this application, the method further includes: The soil samples were obtained under gas-liquid equilibrium by introducing a preset gas into the undisturbed soil at different times. Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

[0011] In some embodiments of this application, determining the three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested includes: The three-dimensional voxel data is obtained by scanning the soil to be tested using a scanning device, and the three-dimensional voxel data is converted into grayscale histogram data. The grayscale histogram data is divided into three phases to determine the gas phase data, liquid phase data, and soil matrix phase data in the grayscale histogram data; wherein, the grayscale range of the gas phase data is greater than 0 and less than or equal to 3500, the grayscale range of the liquid phase data is greater than 3500 and less than or equal to 8000, and the grayscale range of the soil matrix phase data is greater than 8000. The grayscale histogram data after the three phases are divided is processed by grayscale image processing and then reconstructed to obtain three-dimensional images of the gas phase, liquid phase, and soil matrix phase; wherein, the grayscale image processing includes at least one of boundary refinement processing, connectivity analysis processing, artifact region clipping processing, Gaussian smoothing processing, and voxel erosion processing.

[0012] According to another aspect of the embodiments of this application, a soil pore moisture quantification device is provided, the device comprising: The gas filling module is used to fill the undisturbed soil with a preset gas to obtain the soil to be tested under gas-liquid equilibrium. The image determination module is used to determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The result output module is used to determine the moisture distribution and water content of the undisturbed soil pores based on the three-dimensional image.

[0013] According to another aspect of the embodiments of this application, an electronic device is provided, comprising: Controller; The memory is used to store one or more programs that, when executed by the controller, enable the controller to implement the soil pore moisture quantification method described above.

[0014] According to another aspect of the embodiments of this application, a computer-readable storage medium is provided, wherein a computer program is stored in the computer program, which, when run on a soil pore moisture quantification device / electronic device, causes the soil pore moisture quantification device / electronic device to perform the steps of the soil pore moisture quantification method as described above.

[0015] In the soil pore moisture quantification method provided in this application embodiment, a preset gas is introduced into the undisturbed soil to obtain the test soil in a gas-liquid equilibrium state. This allows the introduced preset gas to fully dissolve in the undisturbed soil, thereby significantly altering the X-ray attenuation characteristics of the undisturbed soil containing water voxels. The grayscale of the soil pores containing water voxels shows a measurable increase, while the grayscale of the soil pores without water remains basically unchanged, thus presenting high contrast. Based on this, a three-dimensional image of the gas phase, liquid phase, and soil matrix phase in the test soil can be determined using a scanning device. Finally, the moisture distribution and water content of the undisturbed soil pores are accurately calculated based on the determined three-dimensional image. Through the technical solution of this application, the clear separation of the liquid and gas phases is achieved through the enhanced mechanism of gas combining with pore water in the soil. It also ensures the reversibility and repeatability of the measurement process, effectively avoiding the chemical and physical disturbances that affect the original moisture state caused by using liquid tracers to detect soil structure. It can effectively ensure the stability and accuracy of soil pore moisture distribution identification without changing the moisture structure and natural state within the soil pores.

[0016] The above description is merely an overview of the technical solutions of the embodiments of this application. In order to better understand the technical means of the embodiments of this application and to implement them in accordance with the contents of the specification, and to make the above and other objects, features and advantages of the embodiments of this application more obvious and understandable, specific implementation methods of this application are described below. Attached Figure Description

[0017] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application. It is obvious that the drawings described below are merely some embodiments of this application, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort.

[0018] Figure 1 A schematic flowchart of one embodiment of the soil pore moisture quantification method of this application is shown.

[0019] Figure 2 This application shows Figure 1 A flowchart of an embodiment of S100.

[0020] Figure 3 A schematic diagram of the gas filling device in the soil pore moisture quantification method of this application is shown.

[0021] Figure 4 This application shows Figure 1 A flowchart of another embodiment of S100.

[0022] Figure 5 The diagram shows CT grayscale images before and after carbon dioxide injection in the soil pore moisture quantification method of this application.

[0023] Figure 6 The diagram shows a three-dimensional reconstruction of the pore size and moisture distribution enhanced by potassium iodide-barium chloride in the relevant technology.

[0024] Figure 7 A three-dimensional reconstruction of water distribution enhanced by potassium iodide-barium chloride is shown in the related technology.

[0025] Figure 8 The paper presents a three-dimensional reconstruction of the soil pore moisture distribution enhanced by carbon dioxide gas filling in the soil moisture quantification method of this application.

[0026] Figure 9 The paper presents a three-dimensional reconstruction of water distribution enhanced by introducing carbon dioxide gas in the soil pore water quantification method of this application.

[0027] Figure 10 A schematic diagram of an embodiment of the soil pore moisture quantification device of this application is shown.

[0028] Figure 11 A schematic diagram of an embodiment of the electronic device provided in this application is shown.

[0029] Explanation of reference numerals in the attached figures: 1-Pressure reducing valve; 2-Flow meter; 3-Gas source cylinder; 4-Inlet valve; 5-Sealed cavity. Detailed Implementation

[0030] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0031] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0032] The flowcharts shown in the accompanying drawings are merely illustrative and do not necessarily include all content and operations / steps, nor do they necessarily have to be performed in the described order. For example, some operations / steps can be broken down, while others can be combined or partially combined; therefore, the actual execution order may change depending on the specific circumstances.

[0033] In this application, "multiple" refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A alone, A and B simultaneously, or B alone. The character " / " generally indicates that the preceding and following related objects have an "or" relationship.

[0034] First, a brief overview of the background of this application is provided below: Soil pore structure and its moisture distribution characteristics are key parameters for studying soil physical properties, water transport patterns, and the plant root growth environment. Currently, X-ray CT combined with 3D visualization software is widely used for the observation and analysis of soil microstructure. In some technical solutions, undisturbed soil samples are directly CT scanned to obtain 3D grayscale voxel data of the soil matrix, pore water, and pore gas. To distinguish between the liquid and gas phases, some related technical solutions employ the method of injecting liquid tracers such as potassium iodide into the soil. The high X-ray attenuation characteristics of heavy metal elements are used to enhance the grayscale value of water-bearing areas, and then image segmentation algorithms are used to identify moisture distribution and calculate water content.

[0035] However, in these related schemes, because the hydrogen and oxygen elements in natural soil pore water have X-ray attenuation coefficients that are quite similar to those of the soil matrix, there is a significant overlap in the grayscale ranges between water-containing and water-free voxels, making it difficult to achieve accurate phase separation. Furthermore, introducing liquid tracers alters the chemical composition and surface tension of soil pore water, disturbing the original water distribution. This process is irreversible, making it impossible to repeatedly measure the same sample to monitor the dynamic changes in water content.

[0036] In view of this, to address the aforementioned problems, this application proposes a method for quantifying soil pore moisture, which effectively ensures the stability and accuracy of soil pore moisture distribution identification without altering the moisture structure and natural state within the soil pores. The subject executing this method can be a terminal device, server, or other processing device with computing capabilities. The terminal device can be a user equipment (UE), computer, mobile device, user terminal, terminal, personal digital assistant (PDA), computing device, wearable device, etc. In some possible implementations, this method can be implemented by a processor calling computer-readable instructions stored in memory.

[0037] Specifically, please refer to, for example Figure 1 As shown, the soil pore water quantification method of this embodiment can be applied to farmland irrigation regulation, quantitative characterization of pore water, and dynamic monitoring of seepage and evaporation processes; the method includes the following steps: S100. The undisturbed soil is filled with a preset gas to obtain the soil to be tested under gas-liquid equilibrium. The term "original soil" refers to natural soil, which may be, but is not limited to, typical red soil, lateritic red soil, kaolin, etc. No specific limitation is made here. The choice is based on the actual application scenario. For example, in this embodiment, a typical red soil column with a length of 0-20cm can be selected.

[0038] The preset gas can refer to a gaseous contrast agent that is injected into the undisturbed soil and dissolved in the soil pore water. The core component of this preset gas is carbon dioxide (CO2), with a volume fraction greater than or equal to 99%, meaning that high-purity carbon dioxide is used as the main gas component. Other gas components less than 1% can include nitrogen, oxygen, etc. In this embodiment, if the preset gas contains a large amount of low-solubility impurities such as nitrogen and oxygen, the change in soil pore water density will be insignificant, thus failing to produce a distinguishable grayscale difference in CT images. Therefore, this embodiment uses high-purity carbon dioxide to ensure that the preset gas has sufficient dissolving capacity, improving the grayscale contrast between the moisture-free soil pores and the liquid phase. That is, CO2 dissolves into the pore water until the partial pressure of the gas phase and the concentration of the dissolved liquid phase reach an approximate equilibrium state described by Henry's Law. At this point, a stable gas-liquid equilibrium state is formed inside the soil, resulting in a significant increase in the X-ray attenuation coefficient in the soil pore water, while the grayscale of the moisture-free soil pores remains basically unchanged. This allows for high contrast, thereby achieving non-invasive, high-precision moisture quantification without disturbing the soil's moisture content.

[0039] For example, in practical operation, industrial-grade high-purity carbon dioxide with a purity of 99.90% can be selected and injected at a rate of 30-50 mL / min using a flow meter to ensure that almost all the gas entering the soil pores is carbon dioxide. This high-concentration gas injection fully utilizes the high solubility of carbon dioxide in water, causing a significant increase in the density of pore water after absorbing carbon dioxide. This, in turn, enables a clear distinction of grayscale values ​​between water-containing and water-free pore voxels, achieving non-invasive, high-precision moisture quantification.

[0040] In some embodiments of this application, please refer to, for example... Figure 2 and Figure 3 As shown, in order to illustrate the specific implementation of filling the undisturbed soil with a preset gas to obtain the soil under gas-liquid equilibrium, this application also includes at least S111-S112, the steps of which are as follows: S111. Place the undisturbed soil into a sealed cavity with an air inlet valve; S112. Control the preset gas to be filled into the sealed cavity through the air inlet valve of the sealed cavity to obtain the soil to be tested in a gas-liquid equilibrium state.

[0041] The sealed cavity can refer to a container used to hold the original soil and isolate it from the external environment. The sealed cavity can be, but is not limited to, PVC (polyvinyl chloride) pipe or plexiglass pipe. The choice is made according to the actual application scenario, and no specific limitation is made here.

[0042] The intake valve is a one-way or two-way control interface set on the sealed cavity. It is used to connect the external air source with the internal space of the cavity and is the intake channel for controlling the injection of carbon dioxide gas.

[0043] For example, a soil column with an inner diameter of 2 cm and a height of 3 cm is selected. A PVC pipe with a wall thickness of 1 mm is chosen as the sealing cavity. The collected undisturbed red soil column is pushed into the pipe, and ventilated end caps with annular sealing rings are installed at both ends of the sealing cavity. One end cap integrates a miniature needle valve as an air inlet valve. Please refer to [example]. Figure 3 As shown in the figure, this diagram illustrates the structure of the gas injection device in the soil pore moisture quantification method. The diagram shows the connection between the sealed cavity 5, the gas source cylinder 3, the pressure reducing valve 1, and the flow meter 2 in the gas injection device. The sealed cavity 5 can be connected to the inlet valve 4 via a flexible conduit, forming a closed gas circulation path. Controlled injection refers to adjusting the pressure and flow rate of the gas source cylinder to allow the preset gas to flow through the inlet valve at a stable rate into the sealed cavity. During this process, the preset gas first replaces the original air in the macropores of the soil, and then diffuses into the micropores under the drive of the partial pressure difference and dissolves in the pore water.

[0044] In this embodiment, the gas source cylinder 3 is opened, and the output pressure is maintained at a set value through the pressure reducing valve 1. The gas flow rate is monitored by the flow meter 2, allowing high-purity preset gas to be continuously injected into the sealed cavity 5 through the inlet valve 4. For example, carbon dioxide gas is controlled to be injected into a PVC sealed cavity containing saturated water-containing soil at a rate of 40 mL / min, and the inflation is continued for 60 seconds to ensure that the pore air is completely replaced. Then, the inlet valve is closed, and the cavity is left to stand for 24 hours. Through this embodiment, the undisturbed soil is placed in a sealed cavity of a specific material, and combined with the precise control of the inlet valve, the efficient injection and uniform distribution of preset gas are achieved. On this basis, the controlled inflation process and the post-inflation settling equilibrium mechanism work together to promote the preset gas to reach a stable dissolution equilibrium in the pore water, thereby significantly amplifying the grayscale contrast of the liquid and gas phases in CT images without changing the natural water content of the soil.

[0045] In some embodiments of this application, please refer to, for example... Figure 4 As shown, in order to illustrate the process of filling undisturbed soil with a preset gas and allowing it to stand to obtain the soil under gas-liquid equilibrium, this application also includes at least S121-S122, the steps of which are as follows: S121. Inflate the undisturbed soil with a preset gas according to preset inflating conditions; wherein, the preset inflating conditions include maintaining a preset pressure, preset inflating rate, preset temperature and preset humidity, and continuously inflating a preset gas concentration for a first preset duration. S122. The original soil filled with the preset gas is left to stand for a second preset time to obtain the soil to be tested in a gas-liquid equilibrium state.

[0046] Among them, the preset filling conditions can refer to a series of environmental and operational parameters set in order to cause the preset gas to reach a dissolution equilibrium approximately described by Henry's Law in the soil pore water. The preset filling conditions include maintaining preset pressure, preset filling rate, preset temperature and preset humidity for a first preset duration of continuous filling of preset gas at a preset concentration.

[0047] In some embodiments of this application, to illustrate the specific parameter limitations for obtaining the test soil in a gas-liquid equilibrium state by inflating undisturbed soil with a preset gas, this embodiment includes at least: inflating the undisturbed soil with a preset gas according to preset inflating conditions, the preset inflating conditions including maintaining a preset pressure, a preset inflating rate, a preset temperature, and a preset humidity for a first preset duration of inflating a preset gas concentration; the preset pressure range in the preset inflating conditions is 0.05-0.20 MPa, the preset inflating rate range is 30-50 mL / min, the preset temperature range is 15-30℃, the preset humidity range is 0-10%, and the preset concentration range is 99%-100%; and allowing the undisturbed soil inflated with the preset gas to stand for a second preset duration to obtain the test soil in a gas-liquid equilibrium state.

[0048] The first preset duration can refer to the duration of time during which preset gas is filled into the sealed cavity, so that enough carbon dioxide is filled into the sealed cavity. The first preset duration can be, but is not limited to, 1 minute, 1.5 minutes, 2 minutes, etc., depending on the actual application scenario. This is only an example and is not specifically limited.

[0049] The second preset duration refers to the duration for which the sealed cavity remains closed and the soil sample is kept in a static environment after inflation stops. This second preset duration can be, but is not limited to, 20 hours, 24 hours, 30 hours, etc., depending on the specific application scenario. This is merely an example and not a specific limitation. By setting the second preset duration, the principle of molecular diffusion is utilized to allow sufficient time for the high-concentration gas already injected into the pores to dissolve in the pore water until the chemical potential between the gas and liquid phases reaches dynamic equilibrium, achieving a gas-liquid equilibrium state. During the settling process, the gas concentration dissolved in the pore water gradually stabilizes, causing a measurable change in the X-ray attenuation coefficient of the water-containing voxel, while the anhydrous pore gas phase maintains a lower attenuation coefficient, thus forming a clear grayscale difference in subsequent CT scans.

[0050] For example, this application can control the preset filling rate to 40 mL / min, the preset pressure to 0.10 MPa, the preset temperature to 25°C, the preset humidity to 7%, and the first preset duration to 60 seconds through a flow meter; through this multi-parameter coordinated control of the air filling method, high-concentration carbon dioxide gas can be uniformly filled into the soil pores without damaging the original soil structure.

[0051] Through the embodiments of this application, by continuously filling a preset gas for a first preset time according to preset filling conditions, and then allowing it to stand for a second preset time after filling, a synergistic effect of uniform gas distribution and deep dissolution in soil pores is achieved. Specifically, the continuous filling process, maintaining preset pressure, preset filling rate, preset temperature, and preset humidity, not only rapidly replaces the original air in the pores but also establishes a high partial pressure environment and high concentration gradient conducive to gas dissolution. On this basis, the second preset time of standing allows the gas to fully dissolve in the pore water until a gas-liquid equilibrium is reached, significantly improving the accuracy of subsequent three-phase differentiation based on CT scanning and the reliability of moisture quantification results.

[0052] S200. Determine the three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested.

[0053] After obtaining the soil sample in a gas-liquid equilibrium state, it is necessary to determine the three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil sample. This can be achieved by scanning the soil sample using scanning equipment such as CT scanners, thereby obtaining the three-dimensional voxel data of the soil sample. This data is then used for subsequent operations such as conversion, three-phase segmentation, grayscale processing, and three-dimensional image reconstruction based on the three-dimensional voxel data of the soil sample.

[0054] In some embodiments of this application, in order to illustrate the scanning process of the soil to be tested by the scanning device, this application further includes at least the following steps: obtaining three-dimensional voxel data by scanning the soil to be tested with the scanning device, and converting the three-dimensional voxel data into grayscale histogram data; The grayscale histogram data is divided into three phases to determine the gas phase data, liquid phase data, and soil matrix phase data in the grayscale histogram data; After performing grayscale image processing on the grayscale histogram data after the three phases were divided, and then reconstructing, three-dimensional images of the gas phase, liquid phase, and soil matrix phase were obtained.

[0055] The scanning equipment can refer to an X-ray CT scanner, which penetrates the soil sample to be tested and collects projection data from different angles. After processing by a reconstruction algorithm, it generates three-dimensional voxel data containing spatial location and density information. Converting the scanned three-dimensional voxel data into grayscale histogram data can involve statistically analyzing the frequency of each grayscale value in the entire three-dimensional voxel dataset to form a distribution map with grayscale value on the x-axis and frequency on the y-axis. This map visually displays the grayscale aggregation characteristics within the sample. For the grayscale distribution maps of CO2-filled and non-CO2-filled samples in this application, please refer to [reference needed]. Figure 5 As shown.

[0056] For example, the soil sample filled with CO2 in this application was scanned using an X-ray CT device. Before scanning, the CT device underwent flat-field and dark-field correction, and then the parameters of the CT scanning device were set. The parameters of the CT scanning device can be set as follows: Tube voltage: 100kV; tube current: 150μA; exposure time: 400ms; number of projections: 1800; voxel resolution: 0.1mm; reconstruction algorithm: Feldkamp back projection (FBP).

[0057] After the CT scan is completed, three-dimensional voxel data (.tiff sequence) is obtained. Then, the CT voxel data is imported into three-dimensional image visualization and analysis software and converted into grayscale histogram data. The grayscale histogram data is divided into three phases to determine the gas phase data, liquid phase data and soil matrix phase data in the grayscale histogram data. Finally, the gas phase data, liquid phase data and soil matrix phase data after the three-phase division are processed by grayscale image and reconstructed to obtain three-dimensional images of the gas phase, liquid phase and soil matrix phase.

[0058] In this embodiment, grayscale image processing refers to a series of image processing steps performed after the three-phase division to further optimize image quality, eliminate noise interference, and correct geometric errors. Grayscale image processing may include at least one of boundary refinement, connectivity analysis, artifact region cropping, Gaussian smoothing, and voxel erosion. Boundary refinement is used to eliminate jagged artifacts at the segmentation boundary, making the phase interface smoother and more continuous; connectivity analysis is used to determine the adjacency relationship between voxels in three-dimensional space, identify and mark interconnected pore networks or isolated water droplets, thereby distinguishing between real structures and random noise; artifact region cropping is specifically designed to identify and remove non-real signals such as annular bright or dark bands generated near the sample container wall, avoiding their interference with volume calculation; Gaussian smoothing reduces high-frequency noise in the image through convolution operations, improving the signal-to-noise ratio; and voxel erosion is used to shrink object boundaries and remove false connections or excessively large voxel blocks caused by partial volume effects. For example, when processing red soil sample data, a bright ring with abnormal grayscale was found near the tube wall. This ring was removed by artifact region cropping to prevent it from being misidentified as soil matrix. Simultaneously, for tiny, isolated water films, connectivity analysis was used to confirm whether they belonged to effective pore water; if they were merely single-voxel noise, they were discarded. Through this series of collaborative processing steps, the final generated three-dimensional images of the gas phase, liquid phase, and soil matrix phases preserved the true microstructural features, significantly improving the clarity of the phase interfaces and the accuracy of the structure. This provided a high-quality model foundation for subsequent accurate calculations of the volume fraction of each phase and water distribution parameters.

[0059] In the aforementioned scheme, the three-phase division can refer to the process of discretizing continuous data in a grayscale histogram into subsets representing three different physical components based on a preset grayscale threshold range. Gas phase data can refer to a set of voxels with grayscale values ​​greater than 0 and less than or equal to 3500, primarily representing air pores without water filling. Since the X-ray attenuation coefficient of carbon dioxide gas is similar to that of air and much lower than that of water and solids, its grayscale value remains in a low range. Liquid phase data can refer to a set of voxels with grayscale values ​​greater than 3500 and less than or equal to 8000, representing pore water dissolved with carbon dioxide. The dissolution of carbon dioxide significantly increases the electron density and X-ray attenuation capability of water, causing a measurable increase in its grayscale value compared to pure water or an unaerated state, thus forming a clear boundary with the gas phase data. Soil matrix phase data can refer to a set of voxels with grayscale values ​​greater than 8000, representing the mineral particles and organic matter solid skeleton in the soil. Due to its highest density and strongest X-ray attenuation, it occupies the highest grayscale range.

[0060] To distinguish between the gas phase, liquid phase, and soil matrix phase in soil, the aforementioned threshold ranges are determined based on a calibration model established using standard reference samples. Specifically, under the same scanning conditions, the lower limit of the matrix phase is determined by scanning a dry soil sample with 0% moisture content. Since only the soil matrix and pores exist at this point, the gray values ​​greater than 8000 in the scanned image are identified as the soil matrix, while those less than 8000 represent pores. In water-containing soil, since the image gray values ​​less than 8000 represent gas pores and water within the pores, they cannot be distinguished. Therefore, the gas-liquid boundary is determined by comparing saturated soil samples before and after aeration, thereby determining the gray value ranges of the gas phase, liquid phase, and soil matrix phase.

[0061] Based on the standard reference sample, a calibration model was established to determine the standard values, which, after calibration, yielded the first equation for the linear gray-level-attenuation coefficient relationship: ; in, The linear attenuation coefficient is expressed in cm. -1 ; Grayscale value; The scaling factor can represent the attenuation factor of the CT system. Convert to grayscale The scaling ratio; The constant offset can represent the setting when the attenuation coefficient is... The baseline grayscale value when it is 0; and All calibrations are based on the CT system and the scanned soil. For example, in the embodiments of this application, the calibration is based on the CT system scanning of the aqueous phase before and after CO2 dissolution. It can be 3, calibrated. It can be 4.

[0062] Based on the aforementioned linear gray-level-attenuation coefficient relationship, the second relationship for the gray-level difference component is obtained as follows: ; The effective attenuation coefficient of the aqueous phase before CO2 dissolution is: The effective attenuation coefficient of the aqueous phase after dissolving CO2 The grayscale difference component is calculated as follows: ; Due to the measurement Then the corresponding It has approximately 1500–2000 grayscale units, thus achieving a stable distinction between the aqueous phase and the gas phase.

[0063] Through the embodiments of this application, the grayscale increment brought about by the enhanced carbon dioxide dissolution, combined with the calibrated three-phase grayscale range, solves the problem of indistinguishability caused by the overlap of liquid and gas phase grayscale in traditional methods, ensuring the accuracy of the initial classification. Furthermore, by introducing a combination of various image processing techniques such as boundary refinement, connectivity analysis, artifact cropping, Gaussian smoothing, and voxel erosion, the system noise and container edge artifacts generated during the scanning process are effectively removed, the geometric continuity and topological connectivity of the phase interface are optimized, and the classification error caused by partial volume effects is eliminated.

[0064] S300. Determine the moisture distribution and water content of the undisturbed soil pores based on the three-dimensional image.

[0065] In conjunction with the aforementioned embodiments, after obtaining three-dimensional images of the gas phase, liquid phase, and soil matrix phase using a CT system, the voxel number of each phase in the undisturbed soil can be determined, thereby determining the water distribution and moisture content of the undisturbed soil pores.

[0066] In this embodiment, moisture distribution refers to the spatial location, connectivity, and morphological characteristics of water in the soil pore network, while water content refers to the volume ratio of water per unit volume or unit mass of soil. Moisture distribution and water content are calculated through voxel-level statistical analysis of the generated three-dimensional images. Specifically, the number of voxels belonging to the liquid phase data in the three-dimensional image is counted, and this number is compared with the total number of voxels or the number of voxels in the soil matrix phase to obtain the volumetric water content. Simultaneously, three-dimensional visualization rendering technology is used to visually display the continuous path of liquid phase voxels in the pore space and their contact interfaces with other phases. For example, based on the reconstructed three-dimensional model, the pore water volume fraction of a certain red soil sample can be calculated to be 23.4%, and it can be identified that water is mainly distributed in large pores with a diameter greater than 0.5 mm, while exhibiting a discontinuous distribution in micropores.

[0067] In this embodiment, the third relationship for the volume fraction of each phase is calculated using a voxel-level statistical method: ; in, Let i be the number of phase i (gas, water, solid). It is the total prime number.

[0068] Through the embodiments of this application, high-resolution three-dimensional voxel data are obtained by CT scanning, and three-dimensional images of the gas phase, liquid phase and soil matrix phase are accurately constructed by combining grayscale histogram analysis and image segmentation technology; then, based on the three-dimensional images, voxel statistics and spatial analysis are performed to achieve non-destructive and high-precision quantification of the distribution morphology and water content of soil pore water.

[0069] In some embodiments of this application, the method for quantifying soil pore moisture includes the following steps: The soil samples were obtained under gas-liquid equilibrium by introducing a preset gas into the undisturbed soil at different times. Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

[0070] In this embodiment, multiple discrete or continuous time points are selected on the time axis for the same undisturbed soil sample, and the aeration and equilibration operations are repeatedly performed. At each selected time, the undisturbed soil is placed in a sealed cavity with an air inlet valve, and a preset gas is controlled to be introduced into the sealed cavity through the air inlet valve to achieve a gas-liquid equilibrium state. After gas-liquid equilibrium is achieved at each time, the soil under test is scanned using a scanning device to obtain the unique three-dimensional voxel data at that time, which is then converted into grayscale histogram data. The data is then divided into three phases according to a pre-calibrated grayscale threshold range, and grayscale image processing is performed on the divided data. Since the testing standards are consistent at each time, the generated series of three-dimensional images are highly comparable and can intuitively show the subtle changes in pore structure over time.

[0071] In this embodiment, based on the reconstructed 3D images at each time point, the spatial distribution characteristics of soil moisture and quantitative water content indices at that time point are calculated and output. Since the aeration, imaging, and calculation processes are performed separately at different times, a set of water distribution data and water content values ​​arranged in a time series can be output. This reflects the static water content state at a single moment and reveals the dynamic laws governing the migration, transformation, and flux changes of water in soil pores. For example, by comparing the 3D images at time t1 and time t2, the evaporation rate or infiltration depth of water within the time interval can be quantified, providing precise dynamic data support for the study of soil hydraulic properties.

[0072] Through the embodiments of this application, the process of air balancing, three-dimensional imaging and quantitative calculation is repeatedly executed at different times, realizing long-term, continuous and non-destructive moisture monitoring on the same undisturbed soil sample.

[0073] Based on the above embodiments, combined with, for example Figures 5-9 Comparing CT images of soil without CO2 injection and soil with CO2 injection reveals that the grayscale overlap between the gas and aqueous phases is approximately 40% in the un-CO2-injected soil, decreasing to 15% in the CO2-injected soil, resulting in an approximately 35% improvement in voxel recognition accuracy. Grayscale differentiation is improved from... By definition, the method of infusing CO2 into soil improves grayscale discrimination by approximately 60%. This result indicates that CO2, as a gaseous contrast agent, can effectively amplify the grayscale difference between the water and gas phases, ensuring the stability and accuracy of moisture distribution identification.

[0074] Based on this, the preset pressure of CO2 filling was changed to 0.05 MPa, 0.10 MPa, 0.15 MPa and 0.20 MPa. The results showed that when the preset pressure was below 0.05 MPa, CO2 penetration was incomplete and the pore contrast was insufficient; when the preset pressure exceeded 0.20 MPa, micro-disturbances may occur. Therefore, the optimal gas pressure range was 0.10–0.15 MPa, at which the grayscale separation was the greatest and the sample structure was stable.

[0075] To verify the non-invasive advantage of this application, soil infused with CO2 was compared with that infused with a potassium iodide-barium chloride solution (10%+10%) tracer method. The results showed that the liquid tracer produced a significant difference in grayscale after scanning, thus quantifying soil pore moisture. However, due to the need to add the liquid tracer, the calculated pore water content increased by approximately 8% compared to the initial value, indicating a significant disturbance effect. In contrast, the soil moisture content before and after CO2 infusion showed no significant difference. Therefore, this application can maintain the natural moisture state, avoiding chemical or capillary changes, and is suitable for long-term monitoring.

[0076] The following describes an embodiment of the soil pore moisture quantification device of this application, which can be used to execute the soil pore moisture quantification method in the above embodiments of this application. For details not disclosed in the device embodiments of this application, please refer to the embodiments of the soil pore moisture quantification method described above in this application.

[0077] Please see Figure 10 As shown, Figure 10 A schematic diagram of an embodiment of the soil pore moisture quantification device applied to an electronic device according to this application is shown. The soil pore moisture quantification device 400 includes a gas injection module 410, an image determination module 420, and a result output module 430; The gas filling module 410 is used to fill the undisturbed soil with a preset gas to obtain the soil to be tested under gas-liquid equilibrium. Image determination module 420 is used to determine three-dimensional images of the gas phase, liquid phase and soil matrix phase in the soil to be tested; The result output module 430 is used to determine the moisture distribution and water content of the undisturbed soil pores based on the three-dimensional image.

[0078] In some embodiments of this application, based on the foregoing scheme, the gas filling module 410 places the undisturbed soil into a sealed cavity with an air inlet valve; wherein, the sealed cavity is a PVC pipe or an plexiglass pipe; The preset gas is controlled to be introduced into the sealed cavity through the air inlet valve of the sealed cavity to obtain the soil to be tested in a gas-liquid equilibrium state.

[0079] In some embodiments of this application, based on the foregoing scheme, the gas filling module 410 is further configured to fill the undisturbed soil with a preset gas according to preset filling conditions; wherein, the preset filling conditions include a first preset duration for continuously filling a preset gas concentration at a preset pressure, preset filling rate, preset temperature and preset humidity. The original soil filled with the preset gas is left to stand for a second preset time to obtain the soil to be tested in a gas-liquid equilibrium state.

[0080] The preset pressure range is 0.05-0.20 MPa, the preset filling rate range is 30-50 mL / min, the preset temperature range is 15-30℃, the preset humidity range is 0-10%, and the preset concentration range is 99%-100%.

[0081] In some embodiments of this application, based on the foregoing scheme, the volume fraction of carbon dioxide in the preset gas is greater than or equal to 99%.

[0082] In some embodiments of this application, based on the aforementioned scheme, the gas filling module 410 is also used to fill the undisturbed soil with a preset gas at different times to obtain the soil to be tested in a gas-liquid equilibrium state. Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

[0083] In some embodiments of this application, based on the aforementioned scheme, the image determination module 420 is further configured to scan the soil to be tested using a scanning device to obtain three-dimensional voxel data, and convert the three-dimensional voxel data into grayscale histogram data. The grayscale histogram data is divided into three phases to determine the gas phase data, liquid phase data, and soil matrix phase data in the grayscale histogram data; wherein, the grayscale range of the gas phase data is greater than 0 and less than or equal to 3500, the grayscale range of the liquid phase data is greater than 3500 and less than or equal to 8000, and the grayscale range of the soil matrix phase data is greater than 8000. The grayscale histogram data after the three phases are divided is processed by grayscale image processing and then reconstructed to obtain three-dimensional images of the gas phase, liquid phase, and soil matrix phase; wherein, the grayscale image processing includes at least one of boundary refinement processing, connectivity analysis processing, artifact region clipping processing, Gaussian smoothing processing, and voxel erosion processing.

[0084] In the soil pore moisture quantification device provided in this application embodiment, the undisturbed soil is filled with a preset gas through the gas filling module 410 to obtain the test soil in a gas-liquid equilibrium state. This allows the filled preset gas to fully dissolve in the undisturbed soil, thereby significantly changing the X-ray attenuation characteristics in the undisturbed soil containing water voxels. The gray value of the soil pores containing water voxels shows a measurable increase, while the gray value of the soil pores without water remains basically unchanged, thus presenting a high contrast. Based on this, the three-dimensional image of the gas phase, liquid phase, and soil matrix phase in the test soil can be determined by the image determination module 420. Finally, the moisture distribution and water content of the undisturbed soil pores are accurately calculated based on the three-dimensional image determined by the result output module 430. The technical solution of this application achieves clear separation of the liquid and gas phases through the enhanced mechanism of gas binding to pore water in soil. It also ensures the reversibility and repeatability of the measurement process, effectively avoiding the chemical and physical disturbances that affect the original moisture state caused by using liquid tracers to detect soil structure. It can effectively ensure the stability and accuracy of soil pore moisture distribution identification without changing the moisture structure and natural state in the soil pores.

[0085] It should be noted that the soil pore moisture quantification device 400 provided in the above embodiments and the soil pore moisture quantification method provided in the aforementioned embodiments belong to the same concept. The specific way in which each module and unit performs its operation has been described in detail in the method embodiments, and will not be repeated here.

[0086] Figure 11 The diagram illustrates the structure of an embodiment of the electronic device of this application, and also shows the structure of a computer system suitable for implementing the electronic device of this application. The specific embodiments of this application do not limit the specific implementation of the electronic device.

[0087] Please see Figure 11As shown, the electronic device includes: a controller; and a memory for storing one or more programs, which, when executed by the controller, perform the aforementioned method for quantifying soil pore moisture.

[0088] Please continue reading. Figure 11 As shown, the computer system 500 of this electronic device includes a Central Processing Unit (CPU) 501, which can perform various appropriate actions and processes, such as executing the methods described in the above embodiments, based on programs stored in Read-Only Memory (ROM) 502 or programs loaded from storage portion 508 into Random Access Memory (RAM) 503. The RAM 503 also stores various programs and data required for system operation. The CPU 501, ROM 502, and RAM 503 are interconnected via a bus 504. An input / output (I / O) interface 505 is also connected to the bus 504.

[0089] The following components are connected to I / O interface 505: an input section 506 including a keyboard, mouse, etc.; an output section 507 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and speakers, etc.; a storage section 508 including a hard disk, etc.; and a communication section 509 including a network interface card such as a LAN (Local Area Network) card, modem, etc. The communication section 509 performs communication processing via a network such as the Internet. A drive 510 is also connected to I / O interface 505 as needed. A removable medium 511, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., is installed on drive 510 as needed so that computer programs read from it can be installed into storage section 508 as needed.

[0090] Specifically, according to embodiments of this application, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of this application include a computer program product comprising a computer program carried on a computer-readable medium, the computer program including a computer program for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication section 509, and / or installed from removable medium 511. When the computer program is executed by central processing unit (CPU) 501, it performs various functions defined in the system of this application.

[0091] Another aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for quantifying soil pore moisture as described above. This computer-readable storage medium may be included in the electronic device described in the above embodiments, or it may exist independently and not incorporated into the electronic device.

[0092] Another aspect of this application provides a computer program product or computer program comprising at least one executable instruction that, when executed on a soil pore moisture quantification device / electronic device, causes the soil pore moisture quantification device / electronic device to perform the soil pore moisture quantification method as described below: The soil to be tested is obtained by filling the undisturbed soil with a preset gas to obtain the soil under gas-liquid equilibrium state; Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

[0093] In the soil pore moisture quantification method provided in this application embodiment, a preset gas is introduced into the undisturbed soil to obtain the test soil in a gas-liquid equilibrium state. This allows the introduced preset gas to fully dissolve in the undisturbed soil, thereby significantly altering the X-ray attenuation characteristics of the undisturbed soil containing water voxels. The grayscale of the soil pores containing water voxels shows a measurable increase, while the grayscale of the soil pores without water remains basically unchanged, thus presenting high contrast. Based on this, a three-dimensional image of the gas phase, liquid phase, and soil matrix phase in the test soil can be determined using a scanning device. Finally, the moisture distribution and water content of the undisturbed soil pores are accurately calculated based on the determined three-dimensional image. Through the technical solution of this application, the clear separation of the liquid and gas phases is achieved through the enhanced mechanism of gas combining with pore water in the soil. It also ensures the reversibility and repeatability of the measurement process, effectively avoiding the chemical and physical disturbances that affect the original moisture state caused by using liquid tracers to detect soil structure. It can effectively ensure the stability and accuracy of soil pore moisture distribution identification without changing the moisture structure and natural state within the soil pores.

[0094] It should be noted that the computer-readable medium shown in the embodiments of this application can be a computer-readable signal medium or a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, optical fiber, portable compact disc read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this application, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In this application, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying a computer-readable computer program. The transmitted data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. The computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The computer program contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to wireless, wired, etc., or any suitable combination thereof.

[0095] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. Each block in a flowchart or block diagram may represent a module, segment, or portion of code, which contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0096] The units described in the embodiments of this application can be implemented in software or hardware, and the described units can also be located in a processor. The names of these units do not necessarily limit the specific unit itself.

[0097] According to one aspect of the embodiments of this application, a computer system is also provided, including a Central Processing Unit (CPU), which can perform various appropriate actions and processes based on a program stored in read-only memory (ROM) or a program loaded from storage into random access memory (RAM), such as performing the methods described above. Various programs and data required for system operation are also stored in the RAM. The CPU, ROM, and RAM are interconnected via a bus. Input / output (I / O) interfaces are also connected to the bus.

[0098] The following components are connected to the I / O interface: input components including keyboards, mice, etc.; output components including cathode ray tubes (CRTs), liquid crystal displays (LCDs), and speakers; storage components including hard drives; and communication components including network interface cards such as LAN (Local Area Network) cards and modems. The communication components perform communication processing via networks such as the Internet. Drives are also connected to the I / O interface as needed. Removable media, such as disks, optical discs, magneto-optical discs, semiconductor memories, etc., are installed on the drive as needed so that computer programs read from them can be installed into the storage components as required.

[0099] The above content is merely a preferred exemplary embodiment of this application and is not intended to limit the implementation of this application. Those skilled in the art can easily make corresponding modifications or alterations based on the main concept and spirit of this application. Therefore, the scope of protection of this application should be determined by the scope of protection claimed in the claims.

[0100] In practice, the collection and processing of data in this application should strictly comply with the requirements of relevant national laws and regulations, obtain the informed consent or separate consent of the data subject, and carry out subsequent data use and processing within the scope of laws and regulations and the authorization of the data subject.

Claims

1. A method for quantifying soil pore moisture, characterized in that, The method includes: The soil to be tested is obtained by filling the undisturbed soil with a preset gas to obtain the soil under gas-liquid equilibrium state; Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

2. The method for quantifying soil pore moisture according to claim 1, characterized in that, The process of filling undisturbed soil with a preset gas to obtain soil under gas-liquid equilibrium includes: The undisturbed soil is placed in a sealed cavity with an air inlet valve; wherein the sealed cavity is a PVC pipe or an acrylic pipe. The preset gas is controlled to be introduced into the sealed cavity through the air inlet valve of the sealed cavity to obtain the soil to be tested in a gas-liquid equilibrium state.

3. The method for quantifying soil pore moisture according to claim 1, characterized in that, The process of filling the undisturbed soil with a preset gas to obtain the soil to be tested under gas-liquid equilibrium also includes: According to preset filling conditions, the original soil is filled with preset gas; wherein, the preset filling conditions include maintaining preset pressure, preset filling rate, preset temperature and preset humidity, and continuously filling the preset gas concentration for a first preset duration; The original soil filled with the preset gas is left to stand for a second preset time to obtain the soil to be tested in a gas-liquid equilibrium state.

4. The method for quantifying soil pore moisture according to claim 3, characterized in that, The preset pressure range is 0.05-0.20 MPa, the preset filling rate range is 30-50 mL / min, the preset temperature range is 15-30℃, the preset humidity range is 0-10%, and the preset concentration range is 99%-100%.

5. The method for quantifying soil pore moisture according to claim 1, characterized in that, The volume fraction of carbon dioxide in the preset gas is greater than or equal to 99%.

6. The method for quantifying soil pore moisture according to claim 1, characterized in that, The method further includes: The soil samples were obtained under gas-liquid equilibrium by introducing a preset gas into the undisturbed soil at different times. Determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The moisture distribution and water content of the undisturbed soil pores are determined based on the three-dimensional image.

7. The method for quantifying soil pore moisture according to claim 1, characterized in that, The process of determining the three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested includes: The three-dimensional voxel data is obtained by scanning the soil to be tested using a scanning device, and the three-dimensional voxel data is converted into grayscale histogram data. The grayscale histogram data is divided into three phases to determine the gas phase data, liquid phase data, and soil matrix phase data in the grayscale histogram data; wherein, the grayscale range of the gas phase data is greater than 0 and less than or equal to 3500, the grayscale range of the liquid phase data is greater than 3500 and less than or equal to 8000, and the grayscale range of the soil matrix phase data is greater than 8000. The grayscale histogram data after the three phases are divided is processed by grayscale image processing and then reconstructed to obtain three-dimensional images of the gas phase, liquid phase, and soil matrix phase; wherein, the grayscale image processing includes at least one of boundary refinement processing, connectivity analysis processing, artifact region clipping processing, Gaussian smoothing processing, and voxel erosion processing.

8. A device for quantifying soil pore moisture, characterized in that, The device includes: The gas filling module is used to fill the undisturbed soil with a preset gas to obtain the soil to be tested under gas-liquid equilibrium. The image determination module is used to determine three-dimensional images of the gas phase, liquid phase, and soil matrix phase in the soil to be tested; The result output module is used to determine the moisture distribution and water content of the undisturbed soil pores based on the three-dimensional image.

9. An electronic device, characterized in that, include: Controller; A memory for storing one or more programs that, when executed by a controller, cause the controller to implement the soil pore moisture quantification method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when run on a soil pore moisture quantification device / electronic device, causes the soil pore moisture quantification device / electronic device to perform the steps of the soil pore moisture quantification method as described in any one of claims 1 to 7.

Images