Fractured rock mass seepage process simulation method and device, computer equipment and medium
Patent Information
- Application Number
- CN202310487973.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-04-24
- Publication Date
- 2026-09-22
- Estimated Expiration
- 2043-04-24
AI Technical Summary
[0004]本申请实施例的目的在于提出一种裂隙岩体渗流过程模拟方法，以解决相关技术中对岩体渗流模拟的精准性低的问题
[0036]通过获取待测裂隙岩体试样的待测图像，对待测图像进行归一化处理，得到灰度图像，其中，灰度图像包括各个归一化后的像素点的强度值，根据预设的数值计算模型，构造灰度图像的变化函数，其中，每个变化函数描述每个归一化后的像素点的强度值，强度值与测裂隙岩体试样的渗流特性参数具有对应关系，根据变化函数设置待测图像中的每个像素点的渗流特性参数，基于预设的计算参数，对待测图像中的每个像素点的渗流特性参数进行求解，得到待测裂隙岩体试样的裂隙岩体渗流过程模拟结果，由于灰度图像上的每个像素点有不同的强度值也就是灰度值，即每个像素点可以代表岩体的一处，这样就可以通过变化函数来描述待测图像上每个像素点的渗流特性参数，也就是获取裂隙岩体每一处的渗流特性参数，进而基于预设的计算参数和每个像素点的渗流特性参数来实现裂隙岩体渗流过程的模拟，即使得裂隙岩体的每一处都可以被精确模拟，提高了模拟的精度。
Smart Images
Figure CN116630243B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of fractured rock mass seepage technology, and in particular to a method, apparatus, computer equipment and storage medium for simulating the seepage process in fractured rock mass. Background Technology
[0002] Seepage in fractured rock masses is a common engineering geological phenomenon, posing significant risks to engineering construction due to the resulting safety hazards. 35%-40% of hydropower dam instability, over 90% of mountain slope failures, and 60%-65% of coal mine tunnel water inrush accidents are related to seepage in fractured rock masses. Oil and gas resource extraction also involves seepage in fractured rock masses. Although numerical methods in related technologies can simulate the entire process of seepage in fractured rock masses under various seepage conditions—that is, the parameters set in these technologies assume the rock mass to be homogeneous, meaning that porosity, permeability, density, and compressibility are set completely uniformly at all points in the rock mass—the inherent complexity of natural rock masses, with numerous pores and fractures, limits the accuracy of these technologies in describing and simulating the seepage process in fractured rock masses.
[0003] Since describing and simulating the flow distribution of fluids in fractures and assessing the seepage characteristics of fractured rock masses are of great significance for maintaining the safety and stability of rock mass engineering and the effective development and utilization of resources, improving the accuracy of the simulation of seepage processes in fractured rock masses is an urgent technical problem to be solved. Summary of the Invention
[0004] The purpose of this application is to propose a method for simulating the seepage process in fractured rock mass, so as to solve the problem of low accuracy in the simulation of rock mass seepage in related technologies.
[0005] To address the aforementioned technical problems, this application provides a method for simulating seepage processes in fractured rock masses, comprising the following steps:
[0006] Obtain the image of the fractured rock mass sample to be tested;
[0007] The image to be tested is normalized to obtain a grayscale image, which includes the intensity value of each normalized pixel.
[0008] Based on the preset numerical calculation model, a variation function of the grayscale image is constructed, where each variation function describes the intensity value of each normalized pixel; the intensity value corresponds to the seepage characteristic parameters of the fractured rock mass sample.
[0009] The percolation characteristic parameters of each pixel in the image to be tested are set according to the change function;
[0010] Based on preset calculation parameters, the seepage characteristic parameters of each pixel in the image to be tested are solved to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested.
[0011] In some implementations, the image to be tested is a nuclear magnetic resonance image. Acquiring the image of the fractured rock mass sample to be tested includes:
[0012] The fractured rock mass sample to be tested was saturated.
[0013] Obtain nuclear magnetic resonance images of the fractured rock mass sample after saturation treatment.
[0014] In some embodiments, the above method further includes:
[0015] Based on the preset intensity value range and the intensity value of each pixel in the image to be tested, the region type to which the fractured rock mass sample to be tested belongs is determined.
[0016] In some implementations, the intensity value range includes a first range, a second range, and a third range, where the first range is (0.5-1], the second range is (0.25-0.5], and the third range is (0-0.25];
[0017] Based on the preset intensity value range and the intensity value of each pixel in the image to be tested, the region type to which the fractured rock mass sample to be tested belongs is determined, including:
[0018] When the strength value is within the first interval, the region type to which the fractured rock mass sample belongs is determined to be the main fracture region.
[0019] When the strength value is within the second range, the region type of the fractured rock mass sample to be tested is determined to be a secondary fracture region.
[0020] When the strength value is within the third interval, the region type of the fractured rock mass sample to be tested is determined to be the rock matrix region.
[0021] In some implementations, the image to be tested is normalized to obtain a grayscale image, including:
[0022] Obtain the intensity values of the three channels corresponding to each pixel in the image to be tested;
[0023] The intensity values of each pixel are averaged across the three channels to obtain the intensity value of each pixel in the grayscale image.
[0024] In some implementations, the preset calculation parameters include permeability, porosity, rock density, solid compressibility, fluid density, fluid compressibility, and fluid dynamic viscosity coefficient.
[0025] To address the aforementioned technical problems, this application provides a device for simulating the seepage process in fractured rock mass. The device includes:
[0026] The acquisition module is used to acquire the image of the fractured rock mass sample to be tested.
[0027] The normalization module is used to normalize the image to be tested to obtain a grayscale image, wherein the grayscale image includes the intensity value of each normalized pixel.
[0028] The construction module is used to construct the variation function of the grayscale image according to the preset numerical calculation model. Each variation function describes the intensity value of each normalized pixel. The intensity value has a corresponding relationship with the seepage characteristic parameters of the fractured rock mass sample.
[0029] The setting module is used to set the percolation characteristic parameters of each pixel in the image under test according to the change function;
[0030] The solution module is used to solve the seepage characteristic parameters of each pixel in the image to be tested based on preset calculation parameters, so as to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested.
[0031] In some embodiments, the apparatus for simulating seepage processes in fractured rock masses includes:
[0032] The region determination module is used to determine the region type of the fractured rock mass sample to be tested based on the preset intensity value range and the intensity value of each pixel in the image to be tested.
[0033] To address the aforementioned technical problems, this application also provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described method for simulating the seepage process in fractured rock mass.
[0034] To address the aforementioned technical problems, this application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described method for simulating the seepage process in fractured rock mass.
[0035] Compared with related technologies, the embodiments of this application have the following main advantages:
[0036] By acquiring the image of the fractured rock mass sample to be tested, the image is normalized to obtain a grayscale image. The grayscale image includes the intensity values of each normalized pixel. Based on a preset numerical calculation model, a variation function for the grayscale image is constructed. Each variation function describes the intensity value of each normalized pixel, and the intensity values correspond to the seepage characteristic parameters of the fractured rock mass sample. The seepage characteristic parameters of each pixel in the image are set according to the variation function. Based on preset calculation parameters, the seepage characteristic parameters of each pixel in the image are... The calculation is performed to obtain the simulation results of the seepage process of the fractured rock mass sample. Since each pixel on the grayscale image has a different intensity value, that is, a grayscale value, each pixel can represent a part of the rock mass. Thus, the seepage characteristic parameters of each pixel on the image can be described by the variation function, that is, the seepage characteristic parameters of each part of the fractured rock mass can be obtained. Then, based on the preset calculation parameters and the seepage characteristic parameters of each pixel, the simulation of the seepage process of the fractured rock mass can be realized, so that each part of the fractured rock mass can be accurately simulated, thus improving the simulation accuracy. Attached Figure Description
[0037] To more clearly illustrate the solutions in this application, the accompanying drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0038] Figure 1 This is an exemplary system architecture diagram to which this application can be applied;
[0039] Figure 2 This is a schematic flowchart of the method for simulating the seepage process in fractured rock mass provided in the embodiments of this application;
[0040] Figure 3 This is a schematic diagram of the fractured rock mass sample to be tested after saturation treatment, provided in the embodiments of this application;
[0041] Figure 4-a , Figure 4-b and Figure 4-c These are the strength value distribution diagrams of the fractured rock mass samples under test after saturation treatment in the first interval, the second interval, and the third interval, respectively, in the embodiments of this application.
[0042] Figure 5 This refers to the cracked specimen mesh model involved in the embodiments of this application;
[0043] Figure 6 This is a simulation result diagram of seepage in fractured rock mass involved in the embodiments of this application;
[0044] Figure 7 This is a diagram of experimental results used to verify simulation results in the embodiments of this application;
[0045] Figure 8 This is a comparison chart of the error between the experimental results and numerical results involved in the embodiments of this application;
[0046] Figure 9 This is a schematic diagram of the structure of an embodiment of the seepage process simulation device for fractured rock mass provided in this application;
[0047] Figure 10 This is a schematic diagram of the structure of one embodiment of the computer device provided in this application. Detailed Implementation
[0048] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; the terminology used herein in the specification of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application; the terms "comprising" and "having," and any variations thereof, in the specification, claims, and foregoing drawings of this application, are intended to cover non-exclusive inclusion. The terms "first," "second," etc., in the specification, claims, or foregoing drawings of this application are used to distinguish different objects, not to describe a particular order.
[0049] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.
[0050] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings.
[0051] like Figure 1 As shown, system architecture 100 may include terminal devices 101, 102, and 103, a network 104, and a server 105. Network 104 serves as the medium for providing communication links between terminal devices 101, 102, and 103 and server 105. Network 104 may include various connection types, such as wired or wireless communication links, or fiber optic cables, etc.
[0052] Users can use terminal devices 101, 102, and 103 to interact with server 105 via network 104 to receive or send messages, etc.
[0053] Terminal devices 101, 102, and 103 can be various electronic devices with displays and web browsing capabilities, including but not limited to smartphones, tablets, laptops, and desktop computers.
[0054] Server 105 can be a server that provides various services, such as a backend server that supports the pages displayed on terminal devices 101, 102, and 103.
[0055] It should be noted that the simulation method for seepage process in fractured rock mass provided in the embodiments of this application is generally executed by a server / terminal device, and correspondingly, the simulation device for seepage process in fractured rock mass is generally set in the server / terminal device.
[0056] It should be understood that Figure 1 The number of terminal devices, networks, and servers shown is merely illustrative. Depending on implementation needs, any number of terminal devices, networks, and servers can be included.
[0057] In the embodiments of this application, such as Figure 2 As shown, Figure 2 This is a flowchart illustrating the method for simulating seepage processes in fractured rock masses provided in this application embodiment. The specific implementation of the method for simulating seepage processes in fractured rock masses includes:
[0058] S201: Obtain the image of the fractured rock mass sample to be tested.
[0059] The image to be tested can represent an image of a saturated fractured rock mass sample obtained after saturation processing. This image can be a nuclear magnetic resonance image.
[0060] In some implementations, acquiring the image of the fractured rock mass sample to be tested includes:
[0061] The fractured rock mass sample to be tested was saturated.
[0062] Obtain nuclear magnetic resonance images of the fractured rock mass sample after saturation treatment.
[0063] Specifically, the saturation treatment process includes: first, drying the fractured rock mass sample to be tested at a constant temperature of 105℃ for 24 hours; then, placing the dried sample in a YZK-2 type rock vacuum saturation apparatus and evacuating it under vacuum at -0.1 MPa for 10 hours; finally, saturating it with water at a pressure of 10 MPa for 48 hours. Nuclear magnetic resonance (NMR) images of the fractured rock mass sample are then acquired using a nuclear magnetic resonance imaging (NMR) device.
[0064] like Figure 3 As shown, Figure 3 This is a schematic diagram of the saturated treated fractured rock mass sample provided in the embodiments of this application. Figure 3It can be seen that the image of the saturated fractured rock mass sample to be tested obtained by the nuclear magnetic resonance imaging device shows that there is no seepage inside the saturated fractured rock mass sample to be tested.
[0065] S202: Normalize the image to be tested to obtain a grayscale image, wherein the grayscale image includes the intensity value of each normalized pixel.
[0066] The normalization process involves normalizing the intensity values of the three basic colors in the image of the fractured rock mass sample to be tested, resulting in a grayscale image. The intensity value of each normalized pixel in the grayscale image can be defined as the normalized feature value, i.e., the I value.
[0067] In some implementations, the image to be tested is normalized to obtain a grayscale image, including:
[0068] Obtain the intensity values of the three channels corresponding to each pixel in the image to be tested;
[0069] The intensity values of each pixel are averaged across the three channels to obtain the intensity value of each pixel in the grayscale image.
[0070] Specifically, the intensity values R, G, and B of the three channels (red, green, and blue) for each pixel can be represented by the range 0 to 255, respectively, to obtain the intensity value I of each pixel in the grayscale image: The obtained intensity value I is then processed to [0~1], which means that the intensity value range after normalization is changed from [0~255] to [0~1].
[0071] The intensity value corresponds to the seepage characteristic parameters of the fractured rock mass sample being tested. Here, the fractured rock mass sample being tested refers to a sample after saturation treatment. Seepage characteristic parameters can include permeability and porosity (or porosity). A correspondence is established between the I value obtained after normalization and the permeability and porosity of the saturated fractured rock mass sample. For example, a grayscale image of the saturated fractured rock mass sample is obtained using low-field nuclear magnetic resonance imaging (NMR). Different fracture regions are delineated based on different intensity values in the grayscale image, and different fracture regions have different porosities and permeabilities.
[0072] In some embodiments, the above method further includes:
[0073] Based on the preset intensity value range and the intensity value of each pixel in the image to be tested, the region type to which the fractured rock mass sample to be tested belongs is determined.
[0074] Among them, the region type refers to the different crack regions divided according to different intensity values.
[0075] In some implementations, the intensity value range includes a first range, a second range, and a third range, where the first range is (0.5-1), the second range is (0.25-0.5), and the third range is (0-0.25), as shown below. Figure 4-a , Figure 4-b and Figure 4-c The figure shown is a distribution diagram of the strength values of the fractured rock mass sample after saturation treatment in the first interval, the second interval, and the third interval, respectively, in the embodiments of this application.
[0076] Based on the preset intensity value range and the intensity value of each pixel in the image to be tested, the region type to which the fractured rock mass sample to be tested belongs is determined, including:
[0077] When the strength value is within the first interval, the region type to which the fractured rock mass sample belongs is determined to be the main fracture region.
[0078] When the strength value is within the second range, the region type of the fractured rock mass sample to be tested is determined to be a secondary fracture region.
[0079] When the strength value is within the third interval, the region type of the fractured rock mass sample to be tested is determined to be the rock matrix region.
[0080] The region with an intensity value (I value) between 0.5 and 1 is designated as the primary fracture region, the region between 0.25 and 0.5 as the secondary fracture region, and the region between 0 and 0.25 as the relatively intact rock matrix region. For example, the permeability of the primary fracture region is 1 × 10⁻⁶. -11 -1×10 -10 The porosity is 0.7-0.9; the permeability of the main fracture region is 1×10⁻⁶. -12 -1×10 -11 The porosity is 0.4-0.7; the permeability of the main fracture region is 1×10⁻⁶. -13 -1×10 -10 The porosity is 0.7-0.9.
[0081] The correspondence between the I value obtained after normalization and the permeability and porosity of the fractured rock mass sample after saturation treatment is shown in Table 1 below.
[0082] Table 1. Permeability and porosity of fractured samples
[0083]
[0084] S203: Construct a grayscale image variation function based on a preset numerical calculation model, wherein each variation function describes the intensity value of each normalized pixel.
[0085] Specifically, the grayscale image is imported into COMSOL Multiphysics, a multiphysics simulation software. This software, based on the finite element method, simulates the seepage process of fluid (water) in rock mass by solving partial differential equations. A variation function AF is constructed regarding the pixel position (x, y), which characterizes the normalized intensity value of the nuclear magnetic resonance image of the fractured rock mass sample. During the modeling process of the fractured rock mass sample, its geometric model needs to be meshed using a numerical calculation model, such as... Figure 5 As shown, the numerical calculation model can be the fractured sample grid model involved in the embodiments of this application, and then physical parameters are set, including the porosity, permeability, density, compressibility, etc. of the rock mass. Due to the heterogeneity of the rock mass, the porosity, permeability, density, and compressibility at each point of the rock mass may be different. Therefore, ideally, each grid after division should have its own physical parameters. In the past, when simulating the seepage process of fractured rock mass, the physical parameters were set as if the rock mass were homogeneous, that is, the porosity, permeability, density, and compressibility at each point of the rock mass were set to be completely consistent. In comparison, the fractured rock mass seepage process simulated in the embodiments of this application is more accurate.
[0086] S204: Set the percolation characteristic parameters for each pixel in the image to be tested according to the change function.
[0087] The porosity and permeability of each pixel in the image of the fractured rock mass sample after saturation processing are set using a variation function AF. For example, porosity and permeability can be expressed as ep(im(x,y)) and ka(im(x,y)), respectively, where x and y represent the pixel positions. Since each pixel can represent a location in the fractured rock mass, the porosity and permeability at each location in the fractured rock mass can be determined. When setting physical parameters after meshing the simulation of the seepage process in the fractured rock mass using multiphysics simulation software, a variation function with respect to the pixel position (x,y) is constructed. This variation function allows setting the porosity and permeability represented by each pixel in the image, thereby improving the simulation accuracy.
[0088] S205: Based on preset calculation parameters, the seepage characteristic parameters of each pixel in the image to be tested are solved to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested.
[0089] The preset calculation parameters include permeability, porosity, rock density, solid compressibility, fluid density, fluid compressibility, and fluid dynamic viscosity coefficient.
[0090] Specifically, in unsaturated seepage, groundwater level and saturation are not constants but vary over time. To solve the Richards equations, an unsaturated seepage interface is needed. Here, the Richards equations refer to the parabolic differential equations that describe the unsteady flow in unsaturated porous media such as soil, derived from the generalized Darcy's law and the continuity equation. Solving the Richards equations using an unsaturated seepage interface allows for the determination of seepage characteristic parameters. In the unsaturated seepage interface, groundwater level and saturation are treated as two variables, not two parameters. Numerical methods are typically used to simulate the Richards equations in the unsaturated seepage interface. Commonly used numerical methods include the finite difference method and the finite element method. In the finite difference method, the Richards equations can be discretized into a system of difference equations, which can then be solved numerically. In the finite element method, the Richards equations can be expressed as spatial integrals, which can then be solved numerically.
[0091] For example, the preset calculation parameters are shown in Table 2 below. The seepage characteristic parameters of each pixel in the image to be tested are calculated to obtain the flow velocity distribution of the fractured rock mass sample at different time steps, as shown below. Figure 6 As shown, Figure 6 This is a simulation result diagram of seepage in fractured rock mass involved in the embodiments of this application. Figure 6 The simulation results show the effect of different seepage velocities of the fractured rock mass sample under test at 20s, 40s, 60s, 80s and 100s.
[0092] Table 2 Calculation parameters for seepage model of fractured specimens
[0093]
[0094] In this embodiment, the fluid type is assumed to be water, a Newtonian fluid, flowing at an isothermal and isodense density. Sandstone is selected as the sample for the fractured rock mass. Due to the small model size, deformation of the fractured rock mass caused by fluid flow is not considered. The upper and lower boundaries of the model are impermeable boundaries with no flow. The seepage inlet is at the left end, and the outlet is at the right end. Under a constant pressure of 0.5 MPa, the fluid flows into the computational model from the inlet, and the outlet seepage pressure is set to 0. A transient solver is selected for calculation. The transient solver typically uses numerical methods to discretize time and uses a system of discrete-time differential equations to solve the dynamic process. Commonly used transient solvers include the Euler method, the fourth-order Runge-Kutta method, the low-order Runge-Kutta method, and the high-order Runge-Kutta method.
[0095] To verify the reliability of the simulation results, a seepage test was conducted on the fractured rock samples under the same seepage conditions using a low-field nuclear magnetic resonance (NMR) spectrometer and a seepage displacement imaging device. NMR imaging was performed every 20 seconds, and the results are as follows: Figure 7 As shown, Figure 7 This is a diagram of experimental results used to verify simulation results in the embodiments of this application. Figure 8 This is a comparison chart of the errors between the experimental results and numerical results involved in the embodiments of this application. From the comparison of experimental values and simulated values, the overall simulation results better simulate the entire seepage process of real fractured sandstone. The fluid mainly flows through the main fractures with a relatively high flow velocity, while the flow velocity is lower in the relatively intact rock matrix area. The simulation results reproduce the entire process of the fluid flowing freely in the main fractures and gradually filling the micropores of the rock mass.
[0096] By acquiring the image of the fractured rock mass sample to be tested, the image is normalized to obtain a grayscale image. The grayscale image includes the intensity values of each normalized pixel. Based on a preset numerical calculation model, a variation function for the grayscale image is constructed. Each variation function describes the intensity value of each normalized pixel, and the intensity values correspond to the seepage characteristic parameters of the fractured rock mass sample. The seepage characteristic parameters of each pixel in the image are set according to the variation function. Based on preset calculation parameters, the seepage characteristic parameters of each pixel in the image are... The calculation is performed to obtain the simulation results of the seepage process of the fractured rock mass sample. Since each pixel on the grayscale image has a different intensity value, that is, a grayscale value, each pixel can represent a part of the rock mass. Thus, the seepage characteristic parameters of each pixel on the image can be described by the variation function, that is, the seepage characteristic parameters of each part of the fractured rock mass can be obtained. Then, based on the preset calculation parameters and the seepage characteristic parameters of each pixel, the simulation of the seepage process of the fractured rock mass can be realized, so that each part of the fractured rock mass can be accurately simulated, thus improving the simulation accuracy.
[0097] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. This computer program can be stored in a computer-readable storage medium, and when executed, it can include the processes of the embodiments of the methods described above. The aforementioned storage medium can be a non-volatile storage medium such as a magnetic disk, optical disk, or read-only memory (ROM), or random access memory (RAM).
[0098] It should be understood that although the steps in the flowcharts of the accompanying figures are shown sequentially as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the accompanying figures may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the sub-steps or stages of other steps.
[0099] Further reference Figure 9 As a response to the above Figure 2 The implementation of the method shown in this application provides an embodiment of a device for simulating the seepage process in fractured rock mass. This device embodiment is similar to... Figure 2 Corresponding to the method embodiments shown, this device can be specifically applied to various electronic devices.
[0100] like Figure 9 The diagram shown is a structural schematic of an embodiment of the fractured rock mass seepage process simulation device provided in this application. The device further includes: an acquisition module 91, a normalization module 92, a construction module 93, a setting module 94, and a solution module 95.
[0101] The acquisition module 91 is used to acquire the image of the fractured rock mass sample to be tested.
[0102] The normalization module 92 is used to normalize the image to be tested to obtain a grayscale image, wherein the grayscale image includes the intensity value of each normalized pixel.
[0103] The construction module 93 is used to construct the variation function of the grayscale image according to the preset numerical calculation model, wherein each variation function describes the intensity value of each normalized pixel; the intensity value has a corresponding relationship with the seepage characteristic parameters of the fractured rock mass sample.
[0104] Setting module 94 is used to set the seepage characteristic parameters of each pixel in the image to be tested according to the change function;
[0105] The solver module 95 is used to solve the seepage characteristic parameters of each pixel in the image to be tested based on preset calculation parameters, so as to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested.
[0106] In some implementations, the preset calculation parameters in the solver module 95 include permeability, porosity, rock density, solid compressibility, fluid density, fluid compressibility, and fluid dynamic viscosity coefficient.
[0107] In some embodiments, the apparatus for simulating seepage processes in fractured rock masses includes:
[0108] The region determination module is used to determine the region type of the fractured rock mass sample to be tested based on the preset intensity value range and the intensity value of each pixel in the image to be tested.
[0109] In some implementations, the image to be tested is an MRI image, and the acquisition module 91 is further used for:
[0110] The fractured rock mass sample to be tested was saturated.
[0111] Obtain nuclear magnetic resonance images of the fractured rock mass sample after saturation treatment.
[0112] The intensity value range includes a first range, a second range, and a third range. The first range is (0.5-1], the second range is (0.25-0.5], and the third range is (0-0.25].
[0113] In some implementations, the region determination module is also used for:
[0114] When the strength value is within the first interval, the region type to which the fractured rock mass sample belongs is determined to be the main fracture region.
[0115] When the strength value is within the second range, the region type of the fractured rock mass sample to be tested is determined to be a secondary fracture region.
[0116] When the strength value is within the third interval, the region type of the fractured rock mass sample to be tested is determined to be the rock matrix region.
[0117] In some implementations, the normalization module 92 is also used for:
[0118] Obtain the intensity values of the three channels corresponding to each pixel in the image to be tested;
[0119] The intensity values of each pixel are averaged across the three channels to obtain the intensity value of each pixel in the grayscale image.
[0120] Regarding the simulation device for seepage process in fractured rock mass in the above embodiments, the specific methods of operation of each module have been described in detail in the embodiments of the relevant method, and will not be elaborated here.
[0121] To address the aforementioned technical problems, embodiments of this application also provide a computer device. Please refer to [link / reference needed]. Figure 10 , Figure 10 This is a basic structural block diagram of the computer device in this embodiment.
[0122] The computer device 10 includes a memory 101, a processor 102, and a network interface 103 that are interconnected via a system bus. It should be noted that only the computer device 10 with components 101-103 is shown in the figure; however, it should be understood that it is not required to implement all the shown components, and more or fewer components can be implemented alternatively. Those skilled in the art will understand that the computer device described here is a device capable of automatically performing numerical calculations and / or information processing according to pre-set or stored instructions, and its hardware includes, but is not limited to, microprocessors, application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), digital signal processors (DSPs), embedded devices, etc.
[0123] The computer device can be a desktop computer, laptop, handheld computer, or cloud server, etc. The computer device can interact with the user via a keyboard, mouse, remote control, touchpad, or voice control.
[0124] The memory 101 includes at least one type of readable storage medium, including flash memory, hard disk, multimedia card, card-type memory (e.g., SD or D-type fractured rock mass seepage process simulation memory), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, optical disk, etc. In some embodiments, the memory 101 may be an internal storage unit of the computer device 10, such as the hard disk or memory of the computer device 10. In other embodiments, the memory 101 may also be an external storage device of the computer device 10, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., equipped on the computer device 10. Of course, the memory 101 may also include both the internal storage unit and the external storage device of the computer device 10. In this embodiment, the memory 101 is typically used to store the operating system and various application software installed on the computer device 10, such as program code for a method to simulate seepage processes in fractured rock masses. Furthermore, the memory 101 can also be used to temporarily store various types of data that have been output or will be output.
[0125] In some embodiments, the processor 102 may be a central processing unit (CPU), controller, microcontroller, microprocessor, or other data processing chip. The processor 102 is typically used to control the overall operation of the computer device 10. In this embodiment, the processor 102 is used to run program code stored in the memory 101 or process data, for example, to run program code for the simulation method of seepage process in fractured rock mass.
[0126] The network interface 103 may include a wireless network interface or a wired network interface, which is typically used to establish communication connections between the computer device 10 and other electronic devices.
[0127] This application also provides another embodiment, namely, a computer-readable storage medium storing a simulation program for the seepage process of fractured rock mass, which can be executed by at least one processor to perform the steps of the fractured rock mass seepage process simulation method as described above.
[0128] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes several instructions to cause a terminal device (which may be a mobile phone, computer, server, air conditioner, or network device, etc.) to execute the methods described in the various embodiments of this application.
[0129] Obviously, the embodiments described above are only some embodiments of this application, not all embodiments. The accompanying drawings show preferred embodiments of this application, but do not limit the patent scope of this application. This application can be implemented in many different forms; rather, the purpose of providing these embodiments is to provide a more thorough and comprehensive understanding of the disclosure of this application. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing specific embodiments, or make equivalent substitutions for some of the technical features. Any equivalent structures made using the content of this application's specification and drawings, directly or indirectly applied to other related technical fields, are similarly within the scope of patent protection of this application.
Claims
1. A method for simulating seepage processes in fractured rock masses, characterized in that, The method includes: The fractured rock mass sample to be tested was saturated. The nuclear magnetic resonance image of the fractured rock mass sample to be tested was obtained after saturation treatment; The nuclear magnetic resonance image is normalized to obtain a grayscale image, wherein the grayscale image includes the intensity value of each normalized pixel. Based on a preset intensity value range and the intensity value of each pixel in the grayscale image, the region type to which the fractured rock mass sample belongs is determined. The intensity value range includes a first range (0.5-1), a second range (0.25-0.5), and a third range (0-0.25). When the intensity value is within the first range, the region type to which the fractured rock mass sample belongs is determined to be a primary fracture region; when the intensity value is within the second range, the region type to which the fractured rock mass sample belongs is determined to be a secondary fracture region; and when the intensity value is within the third range, the region type to which the fractured rock mass sample belongs is determined to be a rock matrix region. Based on a preset numerical calculation model, a variation function for the pixel position is constructed, wherein each variation function describes the intensity value of each normalized pixel; the intensity value corresponds to the seepage characteristic parameters of the fractured rock mass sample, wherein the seepage characteristic parameters include permeability and porosity; The percolation characteristic parameters of each pixel in the nuclear magnetic resonance image are set according to the change function; Based on preset calculation parameters and using an unsaturated seepage interface to solve the Richards equation, the seepage characteristic parameters of each pixel in the nuclear magnetic resonance image are solved to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested. The preset calculation parameters include permeability, porosity, rock mass density, solid compressibility, fluid density, fluid compressibility, and fluid dynamic viscosity coefficient.
2. The method for simulating seepage processes in fractured rock masses according to claim 1, characterized in that, The normalization process of the nuclear magnetic resonance image to obtain a grayscale image includes: Obtain the intensity values corresponding to the three channels of each pixel in the nuclear magnetic resonance image; The intensity values of the three channels corresponding to each pixel are averaged to obtain the intensity value of each pixel in the grayscale image.
3. A device for simulating seepage processes in fractured rock masses, characterized in that, The device for simulating the seepage process in fractured rock mass includes: The acquisition module is used to saturate the fractured rock mass sample to be tested; and to acquire the nuclear magnetic resonance image of the fractured rock mass sample after saturation. The normalization module is used to normalize the nuclear magnetic resonance image to obtain a grayscale image, wherein the grayscale image includes the intensity value of each normalized pixel. The region determination module is used to determine the region type of the fractured rock mass sample to be tested based on a preset intensity value range and the intensity value of each pixel in the grayscale image. The intensity value range includes a first range (0.5-1), a second range (0.25-0.5), and a third range (0-0.25). When the intensity value is within the first range, the region type of the fractured rock mass sample to be tested is determined to be a primary fracture region; when the intensity value is within the second range, the region type is determined to be a secondary fracture region; and when the intensity value is within the third range, the region type is determined to be a rock matrix region. The construction module is used to construct a variation function of pixel position according to a preset numerical calculation model, wherein each variation function describes the intensity value of each normalized pixel; the intensity value has a corresponding relationship with the seepage characteristic parameters of the fractured rock mass sample, wherein the seepage characteristic parameters include permeability and porosity; The setting module is used to set the percolation characteristic parameters of each pixel in the nuclear magnetic resonance image according to the change function; The solution module is used to solve for the seepage characteristic parameters of each pixel in the nuclear magnetic resonance image based on preset calculation parameters and by solving the Richards equation using an unsaturated seepage interface, so as to obtain the simulation results of the seepage process of the fractured rock mass sample to be tested. The preset calculation parameters include permeability, porosity, rock mass density, solid compressibility, fluid density, fluid compressibility, and fluid dynamic viscosity coefficient.
4. A computer device comprising a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the method for simulating the seepage process in fractured rock mass as described in any one of claims 1 to 3.
5. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the steps of the method for simulating the seepage process in fractured rock mass as described in any one of claims 1 to 3.