A soil pore mass transfer simulation method, device, equipment and storage medium

Abstract

Embodiments of the present application provide a soil pore mass transfer simulation method, device, equipment and storage medium. The method comprises acquiring a soil sample, CT scanning the soil sample, and acquiring continuous two-dimensional slices of the soil sample. The continuous two-dimensional slices are subjected to pore extraction and boundary extraction to obtain pore extraction data and boundary extraction data, and a three-dimensional model of the soil sample is constructed according to the pore extraction data and the boundary extraction data. Finite element analysis is performed on the three-dimensional model to obtain a discretized finite element grid, and the properties of the finite element grid include internal pressure field, flow field and mass transfer field. A coupling calculation formula is established according to the properties of the finite element grid, and the mass transfer of the next time of the finite element grid is predicted according to the mass transfer model, the environmental data and the coupling calculation formula. In this way, effective soil pore mass transfer analysis and prediction can be realized, and a basis is provided for subsequent soil management and environmental protection.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to the field of soil pore mass transfer simulation, and particularly relate to a soil pore mass transfer simulation method, device, equipment and storage medium. BACKGROUND

[0002] The pore structure in soil is one of the most important physical properties of soil, which not only determines the basic mechanical characteristics of soil, but also directly affects the migration process of water, air, various solutes and particles in soil. Therefore, the research on soil pore structure has important theoretical significance and practical application value. In recent years, the research on soil pore structure at home and abroad mainly focuses on the basic level of quantitative description and geometric characterization of pore characteristics, such as the determination of parameters such as porosity, pore size distribution and connectivity. However, there is still a lack of systematic research on how to establish an effective model based on these pore characteristics to further analyze the influence of pore structure on mass transfer process (such as water migration, solute diffusion) and the migration behavior of heavy metal and other pollutants. Especially in the prediction of heavy metal migration and distribution, the existing research is still limited to qualitative analysis, lacks accurate model support, and cannot fully reveal the internal relationship between pore structure and pollutant migration.

[0003] Therefore, how to carry out mass transfer simulation for soil pore structure and carry out effective mass transfer analysis and prediction is a problem to be solved at present. SUMMARY

[0004] According to the embodiments of the present application, a soil pore mass transfer simulation method, device, equipment and storage medium are provided, which can realize effective soil pore mass transfer analysis and prediction, and greatly improve the efficiency of soil environment governance.

[0005] In a first aspect of the present application, a soil pore mass transfer simulation method is provided. The method comprises:

[0006] obtaining a soil sample, performing CT scanning on the soil sample, and obtaining continuous two-dimensional slices of the soil sample;

[0007] performing pore extraction and boundary extraction on the continuous two-dimensional slices to obtain pore extraction data and boundary extraction data, and constructing a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data;

[0008] performing finite element analysis on the three-dimensional model to obtain a discrete finite element grid, and the properties of the finite element grid include internal pressure field, flow field and mass transfer field;

[0009] establishing a coupling calculation formula according to the properties of the finite element grid, and predicting the mass transfer of the next time of the finite element grid according to the mass transfer model, environmental data and the coupling calculation formula.

[0010] In a possible implementation, the CT scanning adopts a micro-CT scanning method.

[0011] In a possible implementation, the pore extraction adopts Avizo software to perform pore extraction, to obtain pore extraction data.

[0012] In a possible implementation, the boundary extraction includes:

[0013] The continuous two-dimensional slices are subjected to vectorization processing by using an R2V vector diagram processing tool, to obtain vectorization data of the continuous two-dimensional slices.

[0014] The boundary extraction data are obtained by performing boundary extraction on the vectorization data by using CAD.

[0015] In a possible implementation, the coupling calculation formula includes a first coupling formula, a second coupling formula, and a third coupling formula.

[0016] The first coupling formula is a coupling formula between an internal pressure field and a flow field established according to Darcy's law.

[0017] The second coupling formula is a coupling formula between a flow field and a mass transfer field established according to a diffusion equation.

[0018] The third coupling formula is a coupling formula between heat conduction and a flow field established according to an energy balance equation.

[0019] In a possible implementation, the first coupling formula is:

[0020]

[0021] wherein v is a flow velocity vector of the flow field (m / s), μ is a dynamic viscosity of the fluid in the flow field (Pa·s), K is a permeability of the soil, is a pressure (Pa), ρ is a fluid density (kg / m 3 ), and g is a gravity acceleration vector.

[0022] The second coupling formula is:

[0023]

[0024] wherein C is a heavy metal concentration, and D is a diffusion coefficient.

[0025] The third coupling formula is:

[0026]

[0027] wherein T is a temperature, c p is a specific heat capacity under constant pressure, k T is a thermal conductivity of the soil.

[0028] In a possible implementation, the mass transfer model is pre-trained by historical environmental data and historical mass transfer data of the finite element grid;

[0029] The mass transfer model is constrained by the coupling calculation formula during the training process.

[0030] In a second aspect of the present application, a soil pore mass transfer simulation device is provided. The device comprises:

[0031] The scanning module is configured to obtain a soil sample, perform CT scanning on the soil sample, and obtain continuous two-dimensional slices of the soil sample.

[0032] The construction module is configured to perform pore extraction and boundary extraction on the continuous two-dimensional slices, obtain pore extraction data and boundary extraction data, and construct a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data.

[0033] The finite element module is configured to perform finite element analysis on the three-dimensional model, and obtain a discrete finite element grid. The properties of the finite element grid include an internal pressure field, a flow field, and a mass transfer field.

[0034] The prediction module is configured to establish a coupling calculation formula according to the properties of the finite element grid, and predict the mass transfer of the next time of the finite element grid according to the mass transfer model, the environmental data, and the coupling calculation formula.

[0035] In a third aspect of the present application, an electronic device is provided. The electronic device comprises a memory and a processor. The memory stores a computer program. When the processor executes the program, the method described above is implemented.

[0036] In a fourth aspect of the present application, a computer readable storage medium is provided. The computer readable storage medium stores a computer program. When the program is executed by a processor, the method according to the first aspect of the present application is implemented.

[0037] The soil pore mass transfer simulation method provided by the embodiments of the present application comprises the following steps: obtaining a soil sample, performing CT scanning on the soil sample to obtain continuous two-dimensional slices of the soil sample, performing pore extraction and boundary extraction on the continuous two-dimensional slices to obtain pore extraction data and boundary extraction data, constructing a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data, performing finite element analysis on the three-dimensional model to obtain a discretized finite element grid, the attributes of the finite element grid including an internal pressure field, a flow field and a mass transfer field, establishing a coupling calculation formula according to the attributes of the finite element grid, and predicting the mass transfer of the next moment of the finite element grid according to the mass transfer model, environmental data and the coupling calculation formula. That is, the continuous two-dimensional slices of the soil sample are obtained by using the CT scanning technology, the complex structure inside the soil sample can be accurately captured, and the actual situation of the soil is more truly reflected. Moreover, the three-dimensional model is constructed through pore extraction and boundary extraction, which provides a basis for subsequent finite element analysis of the soil sample. Then, the mass transfer of the next moment of the finite element grid is predicted by using the mass transfer model, environmental data and the coupling calculation formula, the physical conditions in the actual environment are coupled and reflected through the coupling calculation of multiple physical fields, and then the mass transfer model is constrained to make a more accurate description, thereby providing effective decision support for soil environment monitoring and pollution prevention.

[0038] It should be understood that the content described in the summary section is not intended to limit the key or important features of the embodiments of the present application, nor to limit the scope of the present application. Other features of the present application will become apparent through the following description. BRIEF DESCRIPTION OF DRAWINGS

[0039] The above and other features, advantages and aspects of the embodiments of the present application will become more apparent by describing in detail the following embodiments with reference to the attached drawings in which:

[0040] Figure 1 A flow chart of the soil pore mass transfer simulation method according to the embodiments of the present application;

[0041] Figure 2 An example diagram of pore extraction and boundary extraction according to the embodiments of the present application;

[0042] Figure 3 A structural schematic diagram of a mass transfer model according to the embodiments of the present application;

[0043] Figure 4 A block diagram of a soil pore mass transfer simulation device according to the embodiments of the present application;

[0044] Figure 5 A structural schematic diagram of a terminal device or a server suitable for implementing the embodiments of the present application. DETAILED DESCRIPTION

[0045] In order to make the objects, technical solutions and advantages of the embodiments of the present disclosure clearer, the following will be combined with the drawings in the embodiments of the present disclosure to make a clear and complete description of the technical solutions in the embodiments of the present disclosure. Obviously, the described embodiments are part of the embodiments of the present disclosure, rather than all the embodiments. Based on the embodiments in the present disclosure, all other embodiments obtained by a person of ordinary skill in the art without creative labor fall within the protection scope of the present disclosure.

[0046] In addition, the term "and / or" in this paper is only to describe the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B, which can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in this paper generally represents that the front and rear associated objects are in an "or" relationship.

[0047] Figure 1 A flowchart of a method for simulating soil pore mass transfer according to an embodiment of the present disclosure is shown. Referring to Figure 1 , the method comprises:

[0048] S101, acquiring a soil sample, CT scanning the soil sample, and acquiring continuous two-dimensional slices of the soil sample.

[0049] In this embodiment, the continuous two-dimensional slices of the soil sample are obtained by CT scanning, which provides high-resolution continuous two-dimensional slice images, which helps to capture the details of the internal microstructure of the soil.

[0050] Optionally, the CT scanning adopts a micro-CT scanning method.

[0051] In this embodiment, Micro-CT (Micro-Computed Tomography) is a kind of microcomputer tomography technology, which uses X-ray to non-destructively image the object, has high resolution, and can obtain detailed feature information of the soil sample. The quality and clarity of the scanning can be adaptively selected as needed in the Micro-CT scanning.

[0052] S102, pore extraction and boundary extraction are performed on the continuous two-dimensional slices, pore extraction data and boundary extraction data are obtained, and a three-dimensional model of the soil sample is constructed according to the pore extraction data and the boundary extraction data.

[0053] In this embodiment, the specific pore extraction data and boundary extraction data are extracted from the continuous two-dimensional slices, which provides an accurate data basis for the three-dimensional modeling of the soil sample, and the soil structure is more intuitively understood by using three-dimensional modeling, which is beneficial to subsequent mass transfer simulation and analysis and prediction.

[0054] Optionally, the pore extraction is performed by using Avizo software to obtain the pore extraction data.

[0055] As shown in Figure 2 Fig. 1 is an example diagram of pore extraction and boundary extraction according to an embodiment of the present application, wherein the left (1) is an example diagram of pore extraction data.

[0056] In this embodiment, Avizo is a three-dimensional visualization and analysis software. For the continuous two-dimensional slice data obtained by micro-CT scanning, Avizo can perform denoising, contrast enhancement, threshold segmentation and other processing operations, and then extract the soil pore region to generate pore extraction data.

[0057] Optionally, the boundary extraction includes:

[0058] The continuous two-dimensional slice is vectorized by using an R2V vector diagram processing tool to obtain vectorized data of the continuous two-dimensional slice.

[0059] The boundary extraction data is obtained by performing boundary extraction on the vectorized data by using CAD.

[0060] As shown in Figure 2 Fig. 2 is an example diagram of pore extraction and boundary extraction according to an embodiment of the present application, wherein the right (2) is an example diagram of boundary extraction data.

[0061] In this embodiment, R2V (Raster to Vector) is a software tool that can convert raster images into vector graphics. In the soil pore mass transfer simulation in the present application, R2V can convert the continuous two-dimensional slice obtained by micro-CT scanning into vectorized data. CAD (Computer-Aided Design) is a computer-aided design software, which is widely used in engineering design, drawing and other fields. In the soil pore mass transfer simulation in the present application, the CAD software can obtain the boundary data of the continuous two-dimensional slice from the vectorized data processed by R2V, realizing the effective combination of R2V software and CAD software.

[0062] S103, performing finite element analysis on the three-dimensional model to obtain a discretized finite element mesh, and the properties of the finite element mesh include an internal pressure field, a flow field and a mass transfer field.

[0063] The finite element analysis can be performed by using ABAQUS, ANSYS and MSC, which are widely used finite element analysis (FEA, Finite Element Analysis) software. The internal pressure field is the pressure distribution of the fluid inside the finite element mesh of the soil sample, the flow field is the flow condition of the fluid in the finite element mesh of the soil sample, and the mass transfer field is the transport process of the heavy metal and other substances in the finite element mesh of the soil sample.

[0064] In this embodiment, the complex soil geometry is finely discretized and simulated by finite element analysis, and more accurate physical field distribution is provided for the finite element grid of the soil sample.

[0065] In S104, a coupling calculation formula is established according to the properties of the finite element grid, and the mass transfer of the next time of the finite element grid is predicted according to the mass transfer model, the environmental data and the coupling calculation formula.

[0066] In this embodiment, the mass transfer and distribution of the substance in the soil can be more comprehensively predicted by combining the mass transfer model, the environmental data and the coupling calculation formula.

[0067] Optionally, the coupling calculation formula includes a first coupling formula, a second coupling formula and a third coupling formula.

[0068] The first coupling formula is a coupling formula between the internal pressure field and the flow field established according to Darcy's law.

[0069] The second coupling formula is a coupling formula between the flow field and the mass transfer field established according to the diffusion equation.

[0070] The third coupling formula is a coupling formula between the heat conduction and the flow field established according to the energy balance equation.

[0071] The first coupling formula is a coupling formula between the internal pressure field and the flow field established according to Darcy's law, which aims to describe the basic relationship between the movement of the fluid in the soil pore and the pressure distribution. In the soil pore, the movement of the fluid is affected by the pore structure and the permeability, and Darcy's law is a basic and practical empirical law in fluid mechanics for describing the flow behavior of the fluid in the porous medium, which relates the flow rate of the fluid, the pressure gradient, the viscosity of the fluid and the permeability of the porous medium. It has great significance for establishing the coupling relationship between the internal pressure field and the flow field in the soil pore. The second coupling formula is a coupling formula between the flow field and the mass transfer field established according to the diffusion equation, which aims to describe the relationship between the fluid flow in the soil pore and the transport of the heavy metal or other substance in the medium. The diffusion equation, also known as the convection-diffusion equation, is commonly used to describe the natural diffusion process of the substance (such as solute, heat, etc.) caused by the concentration gradient, as well as the convection transport caused by the fluid movement. It has important application value in the field of soil pore mass transfer simulation. In addition, since the fluid flow in the soil pore affects the heat distribution, and the change of the heat distribution also affects the properties of the fluid, and further affects the flow of the fluid, the third coupling formula between the heat conduction and the flow field established according to the energy balance equation is also added in the present application, which further considers the influence of multiple factors on the soil pore mass transfer when coupling the movement of the fluid in the soil pore.

[0072] In this embodiment, according to Darcy's law, diffusion equation and energy balance equation, the coupling formulas of internal pressure field, flow field and mass transfer field are respectively established, and a complete simulation framework of multi-physical field coupling is constructed.

[0073] Optionally, the first coupling formula is:

[0074]

[0075] Wherein, v is the flow velocity vector of the flow field (m / s), μ is the dynamic viscosity of the fluid in the flow field (Pa·s), K is the permeability of the soil, is the pressure (Pa), ρ is the fluid density (kg / m 3 ), g is the gravity acceleration vector;

[0076] The second coupling formula is:

[0077]

[0078] Wherein, C is the heavy metal concentration, D is the diffusion coefficient;

[0079] The third coupling formula is:

[0080]

[0081] Wherein, T is the temperature, c p is the specific heat capacity under constant pressure, k T is the thermal conductivity of the soil.

[0082] Wherein, some parameters in the coupling calculation formula can be set as follows:

[0083] Property Value Unit Property Group Density 2200 Kg / m 3 ]] Basic Kinematic Viscosity 1e-3 Pa·s Basic Porosity 0.68 1 Basic Permeability 1e-5 m 2 ]]> Basic Heat Capacity at Constant Pressure 650 J / (kg·K) Basic Thermal Conductivity 0.42 W / (m·K) Basic

[0084] In addition, the coupling calculation formula can also include an effective permeability coupling formula, which calculates the effective permeability k eff of the finite element grid through porosity and fluid residual saturation, and the calculation formula is as follows:

[0085] k eff = μ (1-S r ) m ,

[0086] Wherein, μ is the porosity, S r is the fluid residual saturation, and m is an empirical index that can be adjusted according to actual needs. Through the effective porosity coupling formula, more physical field scenarios can be further fitted, so that the final mass transfer prediction is more accurate.

[0087] In this embodiment, the coupling calculation formula is not limited to the internal pressure field, flow field and mass transfer field, and can be additionally and individually set according to the actual physical scene.

[0088] Optionally, the mass transfer model is obtained by pre-training based on historical environmental data and historical mass transfer data of the finite element grid.

[0089] The mass transfer model is constrained by the coupling calculation formula during the training process.

[0090] The historical environmental data includes but is not limited to temperature, humidity and wind speed, and the historical mass transfer data includes but is not limited to the historical average heavy metal concentration in the finite element grid. Examples of the historical environmental data and the historical mass transfer data of the finite element grid are shown as follows:

[0091]

[0092] Figure 3 The structure diagram of the mass transfer model according to the embodiment of the present application is shown as follows: Figure 3

[0093] The mass transfer model includes an input layer, a PINN layer, a fully connected layer and an output layer. The input layer is used for standardizing and preprocessing the historical environmental data and the historical mass transfer data of the finite element grid to be input, so as to conform to the input format of the mass transfer model. The PINN layer (Physics-Informed Neural Networks) is used for feature extraction of the historical environmental data and the historical mass transfer data of the finite element grid after standardization and preprocessing, to obtain characteristic information. The Physics-Informed Neural Networks is a machine learning model combining data-driven and physical laws, which trains the neural network by adding the constraint of the physical scene in the loss function, so that the model not only fits the observed data, but also follows the known physical laws. In the present application, the coupling calculation formula is added during the training process to constrain the model, so as to fit multiple different physical field environments. The fully connected layer is used for mapping the characteristic information of multiple finite element grids output by the PINN layer to a lower dimensional space and performing nonlinear transformation, to prepare for the final prediction result. The output layer is responsible for integrating and predicting the mapping data output by the fully connected layer to obtain the final soil pore mass transfer prediction information. In addition, the loss function used for pre-training of the mass transfer model is the adaptive sum of the loss between each physical field value obtained by the coupling calculation formula and the label value, and the adaptive parameter can be set as needed or learned during the model training process. The cross-entropy loss is a basic loss function for solving the mass transfer model prediction according to the probability distribution.

[0094] ​In a possible implementation, there is a soil sample, after continuous two-dimensional slices are obtained through CT scanning, the continuous two-dimensional slices are subjected to pore extraction and boundary extraction, and a three-dimensional model is constructed according to the pore extraction data and the boundary extraction data. Then, finite element analysis is performed on the three-dimensional model, for example, the three-dimensional model is divided into five finite element grids, and then the current environmental information and the average heavy metal concentration of the five finite element grids are input into the mass transfer model, and the overall heavy metal concentration prediction value of the soil sample corresponding to the five finite element grids at the next moment is obtained through the mass transfer model. Then, according to the overall heavy metal concentration prediction value, it can be further judged whether the heavy metal concentration of the soil sample at the next moment exceeds the standard, and preventive measures can be taken in advance. Further, there is a heavy metal threshold h, when the overall heavy metal concentration prediction value is greater than the heavy metal threshold h, the current soil heavy metal concentration is high, and measures such as applying soil conditioner, adjusting irrigation plan or source control need to be taken to control the soil environment.

[0095] In this embodiment, the PINN physical information neural network combining data driving and physical laws is selected, and the coupling calculation constraint of multiple physical fields is added, so that the model can better fit the actual environment, and the accuracy of the soil pore mass transfer simulation is improved.

[0096] According to the embodiments of the present disclosure, the following technical effects are achieved:

[0097] 1. By combining CT scanning technology, pore extraction and boundary extraction, three-dimensional reconstruction of the soil sample is realized, providing a basis for subsequent finite element analysis and soil pore mass transfer simulation prediction.

[0098] 2. The coupling calculation formula comprehensively considers the pressure field, flow field and mass transfer field in the soil sample, optimizes the accuracy of the mass transfer model by considering the interaction between multiple physical fields, and establishes a calculation model that can accurately describe the soil pore mass transfer process.

[0099] 3. The coupling calculation formula constrains the mass transfer model, making the simulation result closer to the actual situation, effectively utilizing the interaction and influence between different physical processes in the soil, and realizing accurate prediction of the migration and analysis of heavy metals in the soil.

[0100] It should be noted that, for the foregoing method embodiments, in order to simply describe, they are all described as a series of action combinations, but those skilled in the art should know that the present application is not limited to the action order described, because according to the present application, certain steps can be performed in other order or simultaneously. Secondly, those skilled in the art should know that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily required by the present application.

[0101] The above is the introduction of the method embodiment, and the following further describes the scheme of the present application through the device embodiment.

[0102] Figure 4 A block diagram of a soil pore mass transfer simulation device according to an embodiment of the present application is shown, as Figure 4 shown, includes:

[0103] The scanning module 401 is configured to acquire a soil sample, perform CT scanning on the soil sample, and acquire continuous two-dimensional slices of the soil sample.

[0104] The construction module 402 is configured to perform pore extraction and boundary extraction on the continuous two-dimensional slices, acquire pore extraction data and boundary extraction data, and construct a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data.

[0105] The finite element module 403 is configured to perform finite element analysis on the three-dimensional model, and acquire a discretized finite element grid, the attributes of the finite element grid including an internal pressure field, a flow field, and a mass transfer field.

[0106] The prediction module 404 is configured to establish a coupling calculation formula according to the attributes of the finite element grid, and predict the mass transfer of the next time of the finite element grid according to the mass transfer model, the environmental data, and the coupling calculation formula.

[0107] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the described modules can refer to the corresponding process in the foregoing method embodiments, which will not be described herein.

[0108] Figure 5 A structural schematic diagram of a terminal device or a server suitable for implementing an embodiment of the present application is shown.

[0109] As Figure 5 shown, the terminal device or the server includes a central processing unit (CPU) 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage portion 508 to a random access memory (RAM) 503. In the RAM 503, various programs and data required for the operation of the terminal device or the server are also stored. The CPU 501, the ROM 502, and the RAM 503 are connected to each other through a bus 504. An input / output (I / O) interface 505 is also connected to the bus 504.

[0110] The following components are connected to the I / O interface 505: an input part 506 including a keyboard, a mouse, etc.; an output part 507 including a display such as a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage part 508 including a hard disk, etc.; and a communication part 509 including a network interface card such as a LAN card, a modem, etc. The communication part 509 performs communication processing via a network such as the Internet. A drive 510 is also connected to the I / O interface 505 as necessary. A removable media 511 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc. is attached to the drive 510 as necessary, so that a computer program read out therefrom is installed in the storage part 508 as necessary.

[0111] In particular, the above method flow steps can be implemented as a computer software program in accordance with embodiments of the present application. For example, embodiments of the present application include a computer program product comprising a computer program carried on a machine-readable medium, the computer program containing program code for executing the methods illustrated in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via the communication part 509, and / or installed from the removable media 511. When the computer program is executed by the central processing unit (CPU) 501, the above-described functions defined in the system of the present application are executed.

[0112] It should be noted that the computer-readable medium shown in 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,—but not limited to—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 or flash memory), optical fiber, portable compact disk 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 computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. Computer-readable signal media can also be any computer-readable medium other than computer-readable storage media, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wireless, wire, optical fiber, RF, etc., or any suitable combination thereof.

[0113] 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. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing 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 the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, 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.

[0114] The units or modules described in the embodiments of the present application can be implemented in the form of software or in the form of hardware. The units or modules described can also be arranged in a processor. In some cases, the names of the units or modules do not constitute a limitation on the units or modules themselves.

[0115] As another aspect, the present application also provides a computer readable storage medium, which can be included in the electronic device described in the above embodiments, or can exist separately without being assembled into the electronic device. The computer readable storage medium stores one or more programs, and the programs are used by one or more processors to execute the methods described in the present application.

[0116] The above description is merely preferred embodiments of the present application and a description of the principles of the technology used. Those skilled in the art should understand that the scope of the application described in the present application is not limited to the technical solutions formed by the specific combinations of the above technical features, and should also cover other technical solutions formed by the combinations of the above technical features or their equivalent features without departing from the above application concept. For example, the technical solutions formed by the mutual replacement of the above features and the technical features applied in the present application (but not limited to) having similar functions.

Claims

1. A method of simulating mass transfer in soil pores, characterized by, The method comprises the following steps: obtaining a soil sample, performing CT scanning on the soil sample, and obtaining continuous two-dimensional slices of the soil sample; performing pore extraction and boundary extraction on the continuous two-dimensional slices to obtain pore extraction data and boundary extraction data, and constructing a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data; performing finite element analysis on the three-dimensional model to obtain a discretized finite element grid, and the properties of the finite element grid include an internal pressure field, a flow field, and a mass transfer field; establishing a coupling calculation formula according to the properties of the finite element grid, and predicting the mass transfer of the next time of the finite element grid according to a mass transfer model, environmental data, and the coupling calculation formula; the coupling calculation formula comprises a first coupling formula, a second coupling formula, and a third coupling formula; the first coupling formula is a coupling formula between the internal pressure field and the flow field established according to Darcy's law; the second coupling formula is a coupling formula between the flow field and the mass transfer field established according to a diffusion equation; the third coupling formula is a coupling formula between heat conduction and the flow field established according to an energy balance equation; the first coupling formula is: ， wherein is the flow velocity vector (m / s) of the flow field, is the dynamic viscosity of the fluid in the flow field (m ) K is the permeability of the soil, is the pressure (Pa), is the fluid density (kg / m 3 ), is the gravitational acceleration vector; the second coupling formula is: ， wherein C is the concentration of heavy metals, and D is the diffusion coefficient; the third coupling formula is: ， where T is the temperature, Cp is the specific heat capacity at constant pressure, k is the thermal conductivity of the soil.

2. The method of soil pore mass transfer simulation of claim 1, wherein, the CT scanning adopts a micro-CT scanning method.

3. The method of simulating mass transfer in soil pores according to claim 1, wherein, The pore extraction adopts Avizo software to perform pore extraction and obtain the pore extraction data.

4. The method of soil pore mass transfer simulation of claim 1, wherein, The boundary extraction comprises the following steps: performing vectorization processing on the continuous two-dimensional slices by using an R2V vector diagram processing tool to obtain vectorization data of the continuous two-dimensional slices; performing boundary extraction on the vectorization data by using CAD to obtain the boundary extraction data.

5. The method of soil pore mass transfer simulation of claim 1, wherein, The mass transfer model is pre-trained by historical environmental data and historical mass transfer data of the finite element grid; the mass transfer model is constrained by the coupling calculation formula during the training process.

6. A soil pore mass transfer simulation apparatus, characterized by, The method comprises the following steps: a scanning module for obtaining a soil sample, performing CT scanning on the soil sample, and obtaining continuous two-dimensional slices of the soil sample; a construction module for performing pore extraction and boundary extraction on the continuous two-dimensional slices to obtain pore extraction data and boundary extraction data, and constructing a three-dimensional model of the soil sample according to the pore extraction data and the boundary extraction data; a finite element module for performing finite element analysis on the three-dimensional model to obtain a discretized finite element grid, and the properties of the finite element grid include an internal pressure field, a flow field, and a mass transfer field; a prediction module for establishing a coupling calculation formula according to the properties of the finite element grid, and predicting the mass transfer of the next time of the finite element grid according to a mass transfer model, environmental data, and the coupling calculation formula; the coupling calculation formula comprises a first coupling formula, a second coupling formula, and a third coupling formula; the first coupling formula is a coupling formula between the internal pressure field and the flow field established according to Darcy's law; the second coupling formula is a coupling formula between the flow field and the mass transfer field established according to a diffusion equation; the third coupling formula is a coupling formula between heat conduction and the flow field established according to an energy balance equation; The third coupling formula is a coupling formula between heat conduction and the flow field established according to an energy balance equation; The first coupling formula is: ， wherein is the flow velocity vector (m / s) of the flow field, is the dynamic viscosity of the fluid in the flow field (m ), K is the permeability of the soil, is the pressure (Pa), is the fluid density (kg / m 3 ), is the gravitational acceleration vector; The second coupling formula is: ， Wherein, C is the heavy metal concentration, D is the diffusion coefficient; The third coupling formula is: ， where T is temperature, Cp is the specific heat capacity at constant pressure, k is the thermal conductivity of the soil.

7. An electronic device comprising a memory and a processor, said memory having stored thereon a computer program, characterized in that, The processor implements the method of any one of claims 1-5 when executing the computer program.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the method of any one of claims 1-5.

Images