A method, system, medium and device for constructing a three-dimensional network fracture model

Abstract

The application provides a three-dimensional network fracture model construction method, system, medium and equipment, relates to the three-dimensional network fracture simulation technical field, determines basic parameters according to a network fracture geometric model, and establishes a three-dimensional network fracture model, carries out inert solute transport simulation and reactive solute transport simulation on the basis of the groundwater flow field of the three-dimensional network fracture model, to obtain reactive solute transport parameters, and calculates the spatial distribution of the reactive solute in the three-dimensional network fracture model according to the reactive solute transport parameters, so that the transmission state of the reactive solute in the fracture medium is accurately obtained and analyzed, and scientific and reasonable data support is provided for the control of the pollutant transport and reaction process.

Description

Technical Field

[0001] This application relates to the field of three-dimensional network fracture simulation technology, specifically to a method, system, medium, and device for constructing a three-dimensional network fracture model. Background Technology

[0002] Fractured media are key natural channels for groundwater storage and migration, and their fluid flow is primarily controlled by the characteristics of the fracture structure. However, in practical engineering applications, cost constraints and the complexity of fractured media structures make it difficult to characterize them precisely. The equivalent continuous medium method is often used to generalize bedrock as an equivalent porous medium. However, this equivalent method has significant limitations and cannot meet the needs of in-depth research into the seepage laws of fractured media.

[0003] Unlike inert solute transport mechanisms, reactive solute transport is a complex process integrating water flow, solute transport, and chemical reactions, and it is particularly prominent in fractured media with significant heterogeneity. Reactive solute transport in fractured media exhibits anomalous characteristics, such as premature arrival of products, concentration tailing, and average concentrations lower than expected. These phenomena lead to practical engineering problems such as prolonged remediation cycles and lower-than-expected remediation efficiency. Currently, the understanding of the processes and laws governing reactive solute migration is incomplete, especially in structurally complex network fractured systems, where existing pollutant migration prediction models generally have limitations. Therefore, a research scheme for reactive solute transport mechanisms in network fractures is urgently needed. Summary of the Invention

[0004] To address the aforementioned technical problems, this application is proposed. Embodiments of this application provide a method, system, medium, and device for constructing a three-dimensional network fracture model.

[0005] According to one aspect of this application, a method for constructing a three-dimensional network fracture model is provided, comprising: determining basic parameters of the three-dimensional network fracture model based on a network fracture geometric model; wherein the basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures; establishing the three-dimensional network fracture model based on the basic parameters; wherein the groundwater parameters of the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability; solving the groundwater flow field of the three-dimensional network fracture model based on the three-dimensional network fracture model; and conducting inertial analysis on the groundwater flow field. Inert solute transport simulations are performed to determine the inert solute transport parameters of the three-dimensional network fracture model. These inert solute transport parameters include initial concentration, infiltration concentration of the inert solute, and flow rate. Based on these inert solute transport parameters, reactive solute transport simulations are performed on the three-dimensional network fracture model to determine the reactive solute transport parameters of the model. These reactive solute transport parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. Based on these reactive solute transport parameters, the spatial distribution of the reactive solute in the three-dimensional network fracture model is calculated.

[0006] In one embodiment, determining the basic parameters of the three-dimensional network fracture model based on the network fracture geometry model includes: determining the probability density distribution model of the basic parameters based on the network fracture geometry model; wherein, the fracture location satisfies a normal distribution, the fracture trace length satisfies a uniform distribution, and the fracture orientation and the fracture dip angle satisfy a Fisher distribution.

[0007] In one embodiment, establishing the three-dimensional network fracture model based on the basic parameters includes: setting the initial temperature condition of the three-dimensional network fracture model to a constant temperature value, setting the model inlet boundary to a constant pressure boundary, setting the outlet boundary to a pressureless boundary, and setting the remaining fracture walls to fluxless boundaries.

[0008] In one embodiment, establishing the three-dimensional network fracture model based on the basic parameters includes: setting the fracture medium as a porous medium with a porosity of 1, and using the Brinkman equation to describe the fluid flow characteristics in the fracture, thereby obtaining the three-dimensional network fracture model.

[0009] In one embodiment, the step of conducting an inert solute transport simulation on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model includes: conducting an inert solute transport simulation on the groundwater flow field, and using the convection-dispersion equation to describe the solute transport phenomenon to obtain the inert solute transport parameters.

[0010] In one embodiment, the step of performing reactive solute transport simulation on the three-dimensional network fracture model based on the inert solute transport parameters, and determining the reactive solute transport parameters of the three-dimensional network fracture model, includes: the reactive solute transport characteristic equation of the three-dimensional network fracture model is: ; in, u For fluid velocity, C This refers to the solute concentration. D The molecular diffusion coefficient is... R It is a chemical reaction constant. For the Laplace operator.

[0011] In one embodiment, calculating the spatial distribution of the reactive solute in the three-dimensional network fracture model based on the reactive solute transport parameters includes: calculating the reaction rate of the reactive solute in the three-dimensional network fracture model based on the concentration value of the reactive solute, so as to obtain the spatial distribution of the reactive solute in the three-dimensional network fracture model.

[0012] According to another aspect of this application, a system for constructing a three-dimensional network fracture model is provided, comprising: a basic parameter determination module, used to determine the basic parameters of the three-dimensional network fracture model based on a network fracture geometric model; wherein the basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures; a fracture model establishment module, used to establish the three-dimensional network fracture model based on the basic parameters; wherein the groundwater parameters of the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability; a flow field solution module, used to solve the groundwater flow field of the three-dimensional network fracture model based on the three-dimensional network fracture model; and an inertia parameter determination module, used to determine the groundwater flow field of the three-dimensional network fracture model. Inert solute transport simulations are conducted on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model. These parameters include initial concentration, infiltration concentration of the inert solute, and flow velocity. A reaction parameter determination module is used to perform reactive solute transport simulations on the three-dimensional network fracture model based on the inert solute transport parameters to determine the reactive solute transport parameters of the three-dimensional network fracture model. These parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. A spatial distribution calculation module is used to calculate the spatial distribution of reactive solutes in the three-dimensional network fracture model based on the reactive solute transport parameters.

[0013] According to another aspect of this application, a computer-readable storage medium is provided, the storage medium storing a computer program for performing any of the methods described above.

[0014] According to another aspect of this application, an electronic device is provided, comprising: a processor; a memory for storing processor-executable instructions; the processor being configured to perform any of the methods described above.

[0015] This application provides a method, system, medium, and device for constructing a three-dimensional network fracture model. The method involves determining the basic parameters of the three-dimensional network fracture model based on a geometric model of the network fractures. These basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures. Based on these basic parameters, a three-dimensional network fracture model is established. The groundwater parameters of the three-dimensional network fracture model include initial conditions, boundary conditions, and network fracture permeability. Based on the three-dimensional network fracture model, the groundwater flow field of the three-dimensional network fracture model is solved. Inert solute transport simulation is conducted on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model. These inert solute transport parameters include initial concentration, infiltration concentration of inert solute, and flow velocity. Based on the inert solute transport parameters, reactive... Solute transport simulation was used to determine the reactive solute transport parameters in a three-dimensional network fracture model. These parameters included initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. Based on these parameters, the spatial distribution of reactive solutes within the three-dimensional network fracture model was calculated. Basic parameters were determined according to the network fracture geometry model, and a three-dimensional network fracture model was established. Inert solute transport simulation and reactive solute transport simulation were conducted based on the groundwater flow field within the three-dimensional network fracture model to obtain the reactive solute transport parameters. The spatial distribution of reactive solutes within the three-dimensional network fracture model was then calculated based on these parameters, thus accurately identifying and analyzing the transport state of reactive solutes in the fractured medium. This provides scientific and reasonable data support for the control of pollutant transport and reaction processes. Attached Figure Description

[0016] The above and other objects, features, and advantages of this application will become more apparent from the more detailed description of the embodiments of this application in conjunction with the accompanying drawings. The drawings are provided to further illustrate the embodiments of this application and form part of the specification. They are used together with the embodiments of this application to explain the application and do not constitute a limitation thereof. In the drawings, the same reference numerals generally represent the same components or steps.

[0017] Figure 1 This is a flowchart illustrating a method for constructing a three-dimensional network fracture model provided in an exemplary embodiment of this application.

[0018] Figure 2 This is a schematic diagram of the structure of a three-dimensional network fracture model provided in an exemplary embodiment of this application.

[0019] Figure 3This is a schematic diagram of the mesh partitioning structure of a three-dimensional network fracture model provided in an exemplary embodiment of this application.

[0020] Figure 4 This is a schematic diagram of the reactive solute penetration curve of a three-dimensional network fracture model provided in an exemplary embodiment of this application.

[0021] Figure 5 This is a schematic diagram of the structure of a three-dimensional network fracture model construction system provided in an exemplary embodiment of this application.

[0022] Figure 6 This is a structural diagram of an electronic device provided in an exemplary embodiment of this application. Detailed Implementation

[0023] Hereinafter, exemplary embodiments according to this application will be described in detail with reference to the accompanying drawings. Obviously, the described embodiments are merely some embodiments of this application, and not all embodiments of this application. It should be understood that this application is not limited to the exemplary embodiments described herein.

[0024] Figure 1 This is a flowchart illustrating a method for constructing a three-dimensional network fracture model provided in an exemplary embodiment of this application. Figure 1 As shown, the method for constructing this three-dimensional network fracture model includes the following steps: Step 110: Determine the basic parameters of the three-dimensional network fracture model based on the network fracture geometry model.

[0025] The basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures. The network fracture geometry model is a mathematical model used to quantitatively describe and characterize the spatial distribution, morphology, and connectivity of fracture (crack) systems in rock masses, concrete, or other porous media. It mainly includes the following parameters: (1) Individual fracture scale parameters: Occurrence: dip and dip angle (spatial orientation).

[0026] Size: Usually expressed as trace length (length seen from the outcrop), radius, or diameter.

[0027] Opening width: The width of the fissure, which directly affects the fluid flow capacity.

[0028] Shape: Often simplified to a circular, oval, or polygonal disc.

[0029] (2) Network system parameters: Linear density: The number of cracks passing through a unit length of measuring line.

[0030] Areal density: The total length of the crack traces per unit area.

[0031] Volume density: The number or total area of ​​cracks per unit volume.

[0032] Spacing: The vertical distance between adjacent cracks in the same group of cracks.

[0033] Number of groups: The number of dominant fracture groups obtained from attitude clustering analysis.

[0034] Connectivity: The degree to which fissures intersect and connect to form a through path is the key factor determining permeability.

[0035] Step 120: Based on the basic parameters, establish a three-dimensional network fracture model.

[0036] The groundwater parameters for the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability. This application imports the basic parameters of the three-dimensional network fracture model into simulation software for three-dimensional network fracture modeling and determines the groundwater parameters required for the model.

[0037] Step 130: Solve the groundwater flow field of the three-dimensional network fracture model.

[0038] The groundwater flow field refers to the spatial distribution of water head (or pressure), velocity, and direction as groundwater moves through the pores and fractures of soil and rock masses. This application solves for the groundwater flow field of the three-dimensional network fracture model based on the three-dimensional network fracture model and the corresponding groundwater parameters.

[0039] Step 140: Conduct inert solute transport simulation on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model.

[0040] Inert solutes (or tracers) do not chemically react with the solid matrix during transport (e.g., adsorption, desorption, degradation), but are only controlled by physical processes. By conducting inert solute transport simulations, we focus on understanding the complex physical transport mechanisms in fractured media. The inert solute transport parameters include initial concentration, infiltration concentration of the inert solute, and flow velocity. This application, based on the groundwater flow field, adds a rare matter transport interface to the simulation software to conduct inert solute transport simulations and determines the initial concentration, infiltration concentration of the inert solute, flow velocity, and other inert solute transport parameters in the model.

[0041] Step 150: Based on the inert solute transport parameters, perform reactive solute transport simulation on the three-dimensional network fracture model to determine the reactive solute transport parameters of the three-dimensional network fracture model.

[0042] The reactive solute transport parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. Based on the simulation of inert solute transport, a chemical reaction interface is added to the simulation software to simulate reactive solute transport, determining the initial concentration, chemical reaction frequency factor, activation energy, chemical reaction rate, and other reactive solute transport parameters.

[0043] Step 160: Calculate the spatial distribution of reactive solutes in the three-dimensional network fracture model based on the reactive solute transport parameters.

[0044] The core of reactive transport simulation lies in coupling chemical reaction equations with convection-dispersion equations. These reactions alter solute concentrations, which in turn affect their transport behavior. Therefore, after determining the reactive solute transport parameters in the three-dimensional network fracture model, this application further calculates the spatial distribution of reactive solutes within the model, thereby reducing the research difficulty of reactive solute transport and determining how pollutants (such as heavy metals, nuclides, and organic solvents) migrate, transform, and are ultimately retained or degraded in fractured rock masses.

[0045] This application provides a method for constructing a three-dimensional network fracture model. The method involves determining the basic parameters of the three-dimensional network fracture model based on a geometric model of the network fractures. These basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures. Based on these basic parameters, a three-dimensional network fracture model is established. The groundwater parameters of the three-dimensional network fracture model include initial conditions, boundary conditions, and network fracture permeability. Based on the three-dimensional network fracture model, the groundwater flow field of the three-dimensional network fracture model is solved. An inert solute transport simulation is conducted on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model. These inert solute transport parameters include initial concentration, infiltration concentration of inert solute, and flow velocity. Based on the inert solute transport parameters, a reactive solute transport simulation is performed on the three-dimensional network fracture model. This study aims to determine the reactive solute transport parameters in a three-dimensional network fracture model. These parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. Based on these parameters, the spatial distribution of reactive solutes within the three-dimensional network fracture model is calculated. Basic parameters are determined according to the network fracture geometry model, and a three-dimensional network fracture model is established. Inert solute transport simulations and reactive solute transport simulations are conducted based on the groundwater flow field within the three-dimensional network fracture model to obtain the reactive solute transport parameters. Furthermore, the spatial distribution of reactive solutes within the three-dimensional network fracture model is calculated based on these parameters, thereby accurately understanding and analyzing the transport state of reactive solutes in the fractured medium. This provides scientific and reasonable data support for the control of pollutant transport and reaction processes.

[0046] In one embodiment, the specific implementation of step 110 above may be: based on the network fracture geometry model, determine the probability density distribution model of the basic parameters; wherein, the fracture location satisfies a normal distribution, the fracture trace length satisfies a uniform distribution, and the fracture orientation and fracture dip angle satisfy a Fisher distribution.

[0047] Specifically, the location of the crack follows a normal distribution, and the corresponding probability density function is expressed as: ; in, This is the lower bound of the sample values. For the mathematical expectation value, The standard deviation is given if and only if When the time is right, it follows a standard normal distribution.

[0048] The length of the crack trace follows a uniform distribution, and the corresponding probability density function is expressed as: ; in, The maximum value in the range of possible values. It represents the minimum value within the range of possible values.

[0049] The fracture orientation and dip angle follow a Fisher distribution, and the corresponding probability density function is expressed as: ; in, and These represent the degrees of freedom for the fracture orientation and fracture dip angle, respectively. When both values ​​approach infinity, it is considered a normal distribution.

[0050] In one embodiment, the specific implementation of step 120 above may be as follows: the initial temperature condition of the three-dimensional network fracture model is set to a constant temperature value, the model inlet boundary is set to a constant pressure boundary, the outlet boundary is set to a pressureless boundary, and the remaining fracture walls are set to fluxless boundaries.

[0051] In the three-dimensional network fracture model, it is assumed that the fluid and solute only flow and migrate within the fracture medium, without considering the permeation or diffusion with the matrix. The initial temperature condition of the entire three-dimensional network fracture model is set to a constant 293.15K. The heat transfer effect during solute migration is ignored. The model inlet boundary is set as a constant pressure boundary. The seepage state within the model is changed by changing the pressure value at the inlet. The outlet boundary is set as a pressureless boundary, and the remaining fracture walls are set as fluxless boundaries. The influence of the fracture matrix is ​​ignored, and only the seepage and mass transfer process within the fracture is considered.

[0052] In one embodiment, step 120 can be implemented by setting the fracture medium as a porous medium with a porosity of 1, and using the Brinkman equation to describe the fluid flow characteristics in the fracture, thereby obtaining a three-dimensional network fracture model.

[0053] By oversimplifying the fractured medium as a porous medium with a porosity close to 1, the fluid flow characteristics in the fracture are described using the Brinkman equation, whose governing equation is as follows: ; in, Permeability coefficient, The fluid dynamic viscosity coefficient, The viscosity coefficient is the internal viscosity of the porous medium. Indicates fluid velocity. For pressure, Indicates the effect of gravity. The Laplace operator is specifically mathematically represented as follows: ; ; in, x, y, z These are the components in the three directions of the vertical coordinate system.

[0054] In one embodiment, step 140 can be implemented by: conducting an inert solute transport simulation on a groundwater flow field, and using a convection-dispersion equation to describe the solute transport phenomenon to obtain inert solute transport parameters.

[0055] By simulating the transport of inert solutes in a groundwater flow field and using the convection-dispersion equation to describe the solute transport phenomenon, the following equation was obtained: ; in, The fluid velocity; C This refers to the solute concentration. D The molecular diffusion coefficient; For the Laplace operator.

[0056] In one embodiment, step 150 can be specifically implemented as follows: the reactive solute transport characteristic equation of the three-dimensional network fracture model is: ; in, For fluid velocity, C This refers to the solute concentration. D The molecular diffusion coefficient is... R It is a chemical reaction constant. For the Laplace operator.

[0057] In one embodiment, step 160 can be implemented by: calculating the reaction rate of the reactive solute in the three-dimensional network fracture model based on the concentration value of the reactive solute, so as to obtain the spatial distribution of the reactive solute in the three-dimensional network fracture model.

[0058] A bimolecular reaction is the most basic form of chemical reaction, with the reaction form A + B = C. It is a first-order irreversible chemical reaction. The reaction rates of reactants A and B and product C are set as follows: ; ; ; in, , , These are the concentration values ​​of reactants A and B, and product C, respectively. The reaction rate constant is temperature-dependent and is calculated using the following formula: ; in, For the frequency factor (e.g., set to 1.106 m) 3 / (mol·s)), The activation energy (e.g., set to 30.103 J / mol). This is the universal gas constant (e.g., set to 8.314 J / (mol·K)). This refers to the local temperature.

[0059] The following simulation experiment uses the three-dimensional network fracture model constructed above to illustrate this. Figure 2 As shown, five sets of network fracture geometric models with different geometric structures were set up, denoted as PL1, PL2, PL3, PL4, and PL5, respectively. PL1 has a total of 20 fractures, with two sets of fractures having an average orientation of 45° and 135°, an average fracture length of 10 cm, and an average dip angle of 45°. Control groups (PL2, PL3, PL4, and PL5) with different numbers, lengths, orientations, and dip angles were set up based on PL1. Specific geometric parameter settings are shown in Table 1, and the schematic diagram of the geometric model mesh is shown below. Figure 3 As shown (using PL3 as an example).

[0060] Table 1. Geometric parameter settings for network fractures The numerical simulation of inert solute transport is divided into two stages. In the first stage, it is assumed that the material filling the network fractures is water, with an initial concentration of 0, and the inflow solute concentration at the inlet is set to 1 mol / m. 3The solute was allowed to completely fill the fissures in the first stage; in the second stage, the initial concentration of the solute in the fissures was changed to 1 mol / m. 3 The concentration of solute flowing into the inlet is set to 0 to simulate the solute displacement process.

[0061] The numerical simulation of reactive solute transport assumes that the network gaps are filled with substance A, and sets... c A The initial concentration was 1 mol / m 3 Substance B is continuously introduced at the inlet at a concentration of 1 mol / m³. 3 The reaction and migration between solutes are simulated using the chemical reaction form A+B=C. The specific values ​​of other parameters in the numerical model are shown in Table 2.

[0062] Table 2 Numerical Model Parameter Settings like Figure 4 As shown in the solute concentration penetration curve in the network fracture model, with the injection of solute B, it undergoes a bimolecular reaction with the existing solute A in the model to generate a new substance C. Then, as solute A is gradually displaced by solute B, the concentration of product C decreases, eventually being completely displaced. The entire process is accompanied by a gradual decrease in the concentration of solute A and a continuous increase in the concentration of solute B, while the concentration of reactant C does not reach the ideal concentration for complete reaction, with a peak concentration of only 0.2-0.4 mol / m³. 3 Furthermore, the breakthrough curve of the product concentration C exhibits a clear non-normal distribution, and the curve shows a clear non-Fick transport phenomenon.

[0063] Figure 5 This is a schematic diagram of the structure of a three-dimensional network fracture model construction system provided in an exemplary embodiment of this application. Figure 5As shown, the three-dimensional network fracture model construction system 50 includes: a basic parameter determination module 51, used to determine the basic parameters of the three-dimensional network fracture model based on the network fracture geometric model; wherein the basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures; a fracture model establishment module 52, used to establish the three-dimensional network fracture model based on the basic parameters; wherein the groundwater parameters of the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability; a water flow field solution module 53, used to solve the groundwater flow field of the three-dimensional network fracture model based on the three-dimensional network fracture model; and an inertia parameter determination module 54, used to determine the inertia parameter in the groundwater flow field. Inert solute transport simulations are conducted on the drainage flow field to determine the inert solute transport parameters of the three-dimensional network fracture model. These parameters include initial concentration, infiltration concentration of the inert solute, and flow velocity. A reaction parameter determination module 55 is used to simulate reactive solute transport on the three-dimensional network fracture model based on the inert solute transport parameters, determining the reactive solute transport parameters of the three-dimensional network fracture model. These parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. A spatial distribution calculation module 56 is used to calculate the spatial distribution of reactive solutes in the three-dimensional network fracture model based on the reactive solute transport parameters.

[0064] This application provides a system for constructing a three-dimensional network fracture model. A basic parameter determination module 51 determines the basic parameters of the three-dimensional network fracture model based on the network fracture geometric model. These basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures. A fracture model establishment module 52 establishes the three-dimensional network fracture model based on these basic parameters. The groundwater parameters of the three-dimensional network fracture model include initial conditions, boundary conditions, and network fracture permeability. A flow field solution module 53 solves for the groundwater flow field of the three-dimensional network fracture model. An inert parameter determination module 54 simulates inert solute transport in the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model. These inert solute transport parameters include initial concentration, inert solute infiltration concentration, and flow velocity. A reaction parameter determination module 55 determines the reaction parameters based on the inert solute transport parameters. Reactive solute transport was simulated on a three-dimensional network fracture model to determine the reactive solute transport parameters of the three-dimensional network fracture model. These parameters included initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate. A spatial distribution calculation module 56 calculated the spatial distribution of reactive solutes in the three-dimensional network fracture model based on these parameters. Basic parameters were determined according to the network fracture geometry model, and a three-dimensional network fracture model was established. Based on the groundwater flow field of the three-dimensional network fracture model, inert solute transport simulations and reactive solute transport simulations were conducted to obtain the reactive solute transport parameters. The spatial distribution of reactive solutes in the three-dimensional network fracture model was calculated based on these parameters, thereby accurately understanding and analyzing the transport state of reactive solutes in the fractured medium, providing scientific and reasonable data support for the control of pollutant transport and reaction processes.

[0065] In one embodiment, the basic parameter determination module 51 can be further configured to: determine the probability density distribution model of the basic parameters based on the network crack geometry model; wherein the crack location satisfies a normal distribution, the crack trace length satisfies a uniform distribution, and the crack orientation and crack dip angle satisfy a Fisher distribution.

[0066] In one embodiment, the fracture model establishment module 52 can be further configured to: set the initial temperature condition of the three-dimensional network fracture model to a constant temperature value, set the model inlet boundary to a constant pressure boundary, set the outlet boundary to a pressureless boundary, and set the remaining fracture walls to fluxless boundaries.

[0067] In one embodiment, the fracture model building module 52 can be further configured to: set the fracture medium as a porous medium with a porosity of 1, and use the Brinkman equation to describe the fluid flow characteristics in the fracture, thereby obtaining a three-dimensional network fracture model.

[0068] In one embodiment, the inert parameter determination module 54 can be further configured to: conduct inert solute transport simulation on the groundwater flow field, and use the convection-dispersion equation to describe the solute transport phenomenon to obtain inert solute transport parameters.

[0069] In one embodiment, the reaction parameter determination module 55 can be further configured such that the reactive solute transport characteristic equation of the three-dimensional network fracture model is: ; in, u For fluid velocity, C This refers to the solute concentration. D The molecular diffusion coefficient is... R It is a chemical reaction constant. For the Laplace operator.

[0070] In one embodiment, the spatial distribution calculation module 56 can be further configured to: calculate the reaction rate of the reactive solute in the three-dimensional network fracture model based on the concentration value of the reactive solute, so as to obtain the spatial distribution of the reactive solute in the three-dimensional network fracture model.

[0071] Below, for reference Figure 6 This application describes an electronic device according to embodiments thereof. The electronic device may be either or both of a first device and a second device, or a standalone device independent of them, which may communicate with the first device and the second device to receive acquired input signals from them.

[0072] Figure 6 A block diagram of an electronic device according to an embodiment of this application is illustrated.

[0073] like Figure 6 As shown, the electronic device 10 includes one or more processors 11 and memory 12.

[0074] The processor 11 may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device 10 to perform desired functions.

[0075] The memory 12 may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, and the processor 11 may execute the program instructions to implement the methods of the various embodiments of this application described above and / or other desired functions. Various contents such as input signals, signal components, and noise components may also be stored in the computer-readable storage medium.

[0076] In one example, the electronic device 10 may also include an input device 13 and an output device 14, which are interconnected via a bus system and / or other forms of connection mechanism (not shown).

[0077] When the electronic device is a standalone device, the input device 13 can be a communication network connector for receiving the collected input signals from the first device and the second device.

[0078] In addition, the input device 13 may also include, for example, a keyboard, a mouse, etc.

[0079] The output device 14 can output various information to the outside, including determined distance information, direction information, etc. The output device 14 may include, for example, a display, a speaker, a printer, and a communication network and its connected remote output devices, etc.

[0080] Of course, for the sake of simplicity, Figure 6 Only some of the components of the electronic device 10 relevant to this application are shown in this illustration; components such as buses, input / output interfaces, etc., are omitted. In addition, the electronic device 10 may include any other suitable components depending on the specific application.

[0081] In addition to the methods and apparatus described above, embodiments of this application may also be computer program products, which include computer program instructions that, when executed by a processor, cause the processor to perform the steps in the methods according to various embodiments of this application described in the "Exemplary Methods" section above.

[0082] The computer program product can be written in any combination of one or more programming languages ​​to perform the operations of the embodiments of this application. The programming languages ​​include object-oriented programming languages ​​such as Java and C++, as well as conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on the user's computing device, partially on the user's computing device, as a standalone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0083] Furthermore, embodiments of this application may also be computer-readable storage media storing computer program instructions thereon, which, when executed by a processor, cause the processor to perform the steps in the methods according to various embodiments of this application described in the "Exemplary Methods" section above.

[0084] The computer-readable storage medium may be any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable 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.

[0085] The basic principles of this application have been described above with reference to specific embodiments. However, it should be noted that the advantages, benefits, and effects mentioned in this application are merely examples and not limitations, and should not be considered as essential features of each embodiment of this application. Furthermore, the specific details disclosed above are for illustrative and facilitative purposes only, and are not limitations. These details do not limit the application to the necessity of employing the aforementioned specific details for implementation.

[0086] The block diagrams of devices, apparatuses, devices, and systems involved in this application are merely illustrative examples and are not intended to require or imply that they must be connected, arranged, or configured in the manner shown in the block diagrams. As those skilled in the art will recognize, these devices, apparatuses, devices, and systems can be connected, arranged, and configured in any manner. Words such as “comprising,” “including,” “having,” etc., are open-ended terms meaning “including but not limited to,” and are used interchangeably with them. The terms “or” and “and” as used herein refer to the terms “and / or,” and are used interchangeably with them unless the context clearly indicates otherwise. The term “such as” as used herein refers to the phrase “such as but not limited to,” and is used interchangeably with it.

[0087] It should also be noted that in the apparatus, equipment, and methods of this application, the components or steps can be disassembled and / or recombined. These disassemblies and / or recombinations should be considered as equivalent solutions of this application.

[0088] The above description of the disclosed aspects is provided to enable any person skilled in the art to make or use this application. Various modifications to these aspects will be readily apparent to those skilled in the art, and the general principles defined herein can be applied to other aspects without departing from the scope of this application. Therefore, this application is not intended to be limited to the aspects shown herein, but rather to be accorded the widest scope consistent with the principles and novel features disclosed herein.

[0089] The above description has been given for purposes of illustration and description. Furthermore, this description is not intended to limit the embodiments of this application to the forms disclosed herein. Although numerous exemplary aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations thereof.

Claims

1. A method for constructing a three-dimensional network fracture model, characterized in that, include: Based on the network fracture geometry model, the basic parameters of the three-dimensional network fracture model are determined; wherein, the basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures; Based on the aforementioned basic parameters, a three-dimensional network fracture model is established; wherein, the groundwater parameters of the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability; Based on the aforementioned three-dimensional network fracture model, the groundwater flow field of the three-dimensional network fracture model is solved; Inert solute transport simulations were conducted on the groundwater flow field to determine the inert solute transport parameters of the three-dimensional network fracture model; wherein, the inert solute transport parameters include the initial concentration, the infiltration concentration of the inert solute, and the flow velocity; Based on the inert solute transport parameters, reactive solute transport simulations are performed on the three-dimensional network fracture model to determine the reactive solute transport parameters of the three-dimensional network fracture model; wherein, the reactive solute transport parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate; Based on the aforementioned reactive solute transport parameters, the spatial distribution of reactive solutes in the three-dimensional network fracture model is calculated.

2. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The basic parameters for determining the three-dimensional network fracture model based on the network fracture geometry model include: Based on the network fracture geometry model, the probability density distribution model of the basic parameters is determined; wherein, the fracture location satisfies a normal distribution, the fracture trace length satisfies a uniform distribution, and the fracture orientation and the fracture dip angle satisfy a Fisher distribution.

3. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The establishment of the three-dimensional network fracture model based on the basic parameters includes: The initial temperature conditions of the three-dimensional network fracture model are set to a constant temperature value, the model inlet boundary is set to a constant pressure boundary, the outlet boundary is set to a pressureless boundary, and the remaining fracture walls are set to fluxless boundaries.

4. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The establishment of the three-dimensional network fracture model based on the basic parameters includes: The fracture medium is set as a porous medium with a porosity of 1, and the fluid flow characteristics in the fracture are described using the Brinkman equation to obtain the three-dimensional network fracture model.

5. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The process of conducting inert solute transport simulations on the groundwater flow field and determining the inert solute transport parameters of the three-dimensional network fracture model includes: Inert solute transport was simulated in the groundwater flow field, and the convection-dispersion equation was used to describe the solute transport phenomenon to obtain the inert solute transport parameters.

6. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The step of performing reactive solute transport simulations on the three-dimensional network fracture model based on the inert solute transport parameters, and determining the reactive solute transport parameters of the three-dimensional network fracture model, includes: The reactive solute transport characteristic equation of the three-dimensional network fracture model is as follows: ； in, u For fluid velocity, C This refers to the solute concentration. D The molecular diffusion coefficient is... R It is a chemical reaction constant. For the Laplace operator.

7. The method for constructing a three-dimensional network fracture model according to claim 1, characterized in that, The calculation of the spatial distribution of reactive solutes in the three-dimensional network fracture model based on the reactive solute transport parameters includes: Based on the concentration value of the reactive solute, the reaction rate of the reactive solute in the three-dimensional network fracture model is calculated to obtain the spatial distribution of the reactive solute in the three-dimensional network fracture model.

8. A system for constructing a three-dimensional network fracture model, characterized in that, include: The basic parameter determination module is used to determine the basic parameters of the three-dimensional network fracture model based on the network fracture geometry model; wherein, the basic parameters include fracture location, fracture trace length, fracture orientation, fracture dip angle, and number of fractures; The fracture model establishment module is used to establish the three-dimensional network fracture model based on the basic parameters; wherein, the groundwater parameters of the three-dimensional network fracture model include: initial conditions, boundary conditions, and network fracture permeability; The water flow field solution module is used to solve the groundwater flow field of the three-dimensional network fracture model based on the three-dimensional network fracture model. An inert parameter determination module is used to conduct inert solute transport simulation on the groundwater flow field and determine the inert solute transport parameters of the three-dimensional network fracture model; wherein, the inert solute transport parameters include initial concentration, infiltration concentration of inert solute, and flow velocity; The reaction parameter determination module is used to perform reactive solute transport simulation on the three-dimensional network fracture model based on the inert solute transport parameters, and to determine the reactive solute transport parameters of the three-dimensional network fracture model; wherein, the reactive solute transport parameters include initial concentration, chemical reaction frequency factor, activation energy, and chemical reaction rate; The spatial distribution calculation module is used to calculate the spatial distribution of reactive solutes in the three-dimensional network fracture model based on the reactive solute transport parameters.

9. A computer-readable storage medium, characterized in that, The storage medium stores a computer program for performing the method described in any one of claims 1-7.

10. An electronic device, characterized in that, include: processor; Memory used to store the processor's executable instructions; The processor is used to execute the method described in any one of claims 1-7.

Images