CN117745979B - Three-dimensional fracture-pore coupling network simulation generation method and system - Google Patents
Three-dimensional fracture-pore coupling network simulation generation method and system Download PDFInfo
- Publication number
- CN117745979B CN117745979B CN202410189990.1A CN202410189990A CN117745979B CN 117745979 B CN117745979 B CN 117745979B CN 202410189990 A CN202410189990 A CN 202410189990A CN 117745979 B CN117745979 B CN 117745979B
- Authority
- CN
- China
- Prior art keywords
- pore
- fracture
- dimensional
- domain
- random
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 200
- 230000008878 coupling Effects 0.000 title claims abstract description 61
- 238000010168 coupling process Methods 0.000 title claims abstract description 61
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000004088 simulation Methods 0.000 title claims abstract description 45
- 239000011435 rock Substances 0.000 claims abstract description 47
- 239000007787 solid Substances 0.000 claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 24
- 238000010276 construction Methods 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000011800 void material Substances 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 8
- 230000006399 behavior Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 239000003344 environmental pollutant Substances 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 231100000719 pollutant Toxicity 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the technical field of fracture-pore coupling, and particularly discloses a three-dimensional fracture-pore coupling network simulation generation method and system, wherein the method comprises the following steps: constructing a three-dimensional entity model; acquiring fracture characteristics in actual geological conditions, and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; acquiring pore characteristics in actual geological conditions, and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; in the three-dimensional solid model, eliminating cracks at the overlapping part of the cracks and the pores, and constructing a random three-dimensional crack-pore coupling network; and performing gridding treatment on the random three-dimensional fracture-pore coupling network, and establishing a hybrid physical numerical model to realize rock mechanics analysis and simulation. The method can quickly, accurately and efficiently construct the multi-scale random three-dimensional fracture-pore coupling network, and accurately reflect the real distribution condition inside the rock.
Description
Technical Field
The invention relates to the technical field of fracture-pore coupling, in particular to a three-dimensional fracture-pore coupling network simulation generation method and system.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
The fissures and voids often coexist in the rock and soil mass, and their distribution has a critical impact on subsurface fluid migration, geologic hazard assessment, water resource management, and the like. However, the conventional method is mostly limited to independent analysis of the cracks and the pores, and interaction between the cracks and the pores cannot be considered, so that in an actual system, fluid exchange and pressure transfer equivalent between the cracks and the pores have important influence on the overall seepage behavior. Ignoring these effects may lead to inaccurate predictions of the percolation behavior, thereby reducing the accuracy and reliability of the model. In practical applications, this may lead to errors in the prediction of seepage behavior, affecting the reliability of engineering design and decisions. In addition, the fissures and pores often constitute complex porous media systems, and independent analysis does not provide a comprehensive understanding of the overall system behavior, and such limitations may lead to misunderstanding of fluid transport, contamination diffusion, or rock mechanics behavior.
In the prior art, the research on fracture-pore coupling is mainly focused on the construction of a seepage coupling equation of a fictitious two-dimensional model, and the research on the construction of a fracture-pore coupling numerical simulation model is lacked; in addition, in nature, the cracks and the pores are three-dimensional, and the two-dimensional crack-pore coupling cannot reflect the real distribution condition inside the rock.
Although the prior art discloses that the CT scanning and nuclear magnetic imaging technology can be utilized to obtain the scanning images of the cracks and the pores through scanning the real sample for numerical simulation analysis, the three-dimensional crack-pore coupling model has high complexity, and meanwhile, the interaction between the cracks and the pores needs to be considered, so that the numerical simulation analysis method has higher complexity of numerical calculation and higher calculation cost for the three-dimensional crack-pore coupling model.
Disclosure of Invention
In order to solve the problems, the invention provides a three-dimensional fracture-pore coupling network simulation generation method and system, wherein a random three-dimensional fracture network and a random three-dimensional pore network are coupled and numerically processed by a numerical simulation method, so that a multi-scale random three-dimensional fracture-pore coupling network can be quickly, accurately and efficiently constructed, and further, a data body is numerically processed, so that a real three-dimensional fracture-pore structure in a rock can be accurately generated.
In some embodiments, the following technical scheme is adopted:
a three-dimensional fracture-pore coupling network simulation generation method comprises the following steps:
Constructing a three-dimensional entity model based on actual geological conditions;
Acquiring fracture characteristics in actual geological conditions, and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; wherein the fracture characteristics include at least fracture occurrence, fracture length, and fracture number;
acquiring pore characteristics in actual geological conditions, and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; the pore characteristics include at least pore shape, porosity, and pore size;
In the three-dimensional solid model, eliminating cracks at the overlapping part of the cracks and the pores, and constructing a random three-dimensional crack-pore coupling network;
and performing gridding treatment on the random three-dimensional fracture-pore coupling network, and establishing a hybrid physical numerical model to realize rock mechanics analysis and simulation.
The method comprises the following steps of generating a random three-dimensional fractured rock mass model in the three-dimensional solid model through a discrete fracture network, wherein the specific process is as follows:
randomly generating a space coordinate point in the three-dimensional solid model, and generating a crack based on the acquired crack occurrence and crack length characteristics by taking the space coordinate point as a center; repeating the fracture generation process, thereby generating a fracture network; the fracture occurrence comprises average trend and average inclination angle of the fracture, and the fracture length is characterized by minimum trace length, maximum trace length or trace length average value of the fracture.
Generating a random three-dimensional pore rock mass model in the three-dimensional entity model through a discrete pore network, wherein the specific process comprises the following steps of:
Randomly generating a space coordinate point in the three-dimensional solid model, and generating a pore based on the acquired pore shape and pore size characteristics by taking the space coordinate point as a center; and determining the number of pores through the porosity characteristics, and repeating the pore generation process to generate a pore network.
In other embodiments, the following technical solutions are adopted:
a three-dimensional fracture-pore coupled network simulation generation system, comprising:
the three-dimensional entity model construction module is used for constructing a three-dimensional entity model based on actual geological conditions;
The random three-dimensional fracture rock mass model construction module is used for acquiring fracture characteristics in actual geological conditions and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; wherein the fracture characteristics include at least fracture occurrence, fracture length, and fracture number;
The random three-dimensional pore rock mass model construction module is used for acquiring pore characteristics in actual geological conditions and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; the pore characteristics include at least pore shape, porosity, and pore size;
the random three-dimensional crack-pore coupling network module is used for removing cracks at the overlapping part of the cracks and the pores in the three-dimensional solid model to construct a random three-dimensional crack-pore coupling network;
And the model processing module is used for carrying out gridding processing on the random three-dimensional fracture-pore coupling network and establishing a hybrid physical numerical model so as to realize rock mechanics analysis and simulation.
In other embodiments, the following technical solutions are adopted:
a terminal device comprising a processor and a memory, the processor for implementing instructions; the memory is used for storing a plurality of instructions adapted to be loaded by the processor and to perform the three-dimensional fracture-void coupled network simulation generation method described above.
In other embodiments, the following technical solutions are adopted:
A computer readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the three-dimensional fracture-void coupled network simulation generation method described above.
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the invention, the random three-dimensional fracture network and the random three-dimensional pore network are coupled and numerically processed, so that the multi-scale random three-dimensional fracture-pore coupling network can be quickly, accurately and efficiently constructed, the real distribution condition in the rock can be accurately reflected, the numerically processed data body is further numerically processed, the mixed physical model is constructed, and the data body is cut at any angle and grid point information is stored.
(2) According to the seepage physical process of the fluid in the cracks, the pores and the coupling network, the corresponding flow equation is established, model parameters can be set and verified through the existing seepage simulation software, simulation calculation of the hybrid physical model is carried out, and analysis and simulation of rock mechanics can be realized.
Additional features and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional fracture-pore coupled network simulation generation method in an embodiment of the invention;
FIG. 2 is a schematic representation of a random three-dimensional pore morphology in an embodiment of the present invention;
FIG. 3 is a schematic diagram of a random three-dimensional fracture network model in an embodiment of the invention;
FIG. 4 is a schematic diagram of a random three-dimensional pore network model in accordance with an embodiment of the present invention;
FIG. 5 is a schematic diagram of the crack and pore coupling in an embodiment of the present invention;
FIG. 6 is a schematic diagram of a random three-dimensional fracture-pore coupled network model in accordance with an embodiment of the present invention;
FIG. 7 is a numerical gridding schematic diagram of a random three-dimensional fracture-pore coupling model in an embodiment of the present invention;
FIG. 8 is a projection of a random three-dimensional slot-aperture coupling model onto a plane at a height cross-section in an embodiment of the invention.
Detailed Description
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the application. Unless defined otherwise, 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 belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
Example 1
In one or more embodiments, a three-dimensional fracture-pore coupled network simulation generation method is disclosed, and in combination with fig. 1, the method specifically includes the following steps:
(1) Constructing a three-dimensional entity model based on actual geological conditions;
in the embodiment, the shape and size data of the solid model are determined based on the actual geological space where the cracks and the pores are located; the shape of the solid model can be a cube or a cylinder; wherein the dimension data of the cube comprises a length Lc, a width Wc and a height Hc; the dimensional data of the cylinder includes a radius Ro, a height Ho. The specific size data is determined according to actual needs. In addition, more complex model shapes can be built according to actual geological conditions.
(2) Acquiring fracture characteristics in actual geological conditions, and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network;
In this embodiment, the fracture characteristics in the actual geological conditions may be obtained by means of CT scanning, etc.; the fracture characteristics include at least fracture morphology, fracture length, and fracture number; the fracture morphology includes an average strike T a and an average dip D a; the fracture length is characterized by a fracture minimum trace length I min, a maximum trace length I max, or an average of fracture trace lengths I avg; the number of cracks N c is set by grouping, and the crack distribution form may be a random uniform distribution or a normal distribution.
And, depending on the discrete fracture network, the fracture shape generated may be linear, disc-shaped, or octahedral.
The generation process of the random three-dimensional fractured rock mass model comprises the following steps: randomly generating a space coordinate point in the three-dimensional solid model, and generating a crack based on the acquired crack occurrence and crack length characteristics by taking the space coordinate point as a central point and a starting point; the fracture generation process is repeated to generate a fracture network, as shown in fig. 3.
(3) Acquiring pore characteristics in actual geological conditions, and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network;
In this embodiment, the type of multi-scale pores may be microscopic fissures, macroscopic fissures, or dissolved pores, and the pore characteristics include at least pore shape, porosity, and pore size; with reference to fig. 2, the shape of the pores is mainly small sphere, and the change of the pore shape is realized by controlling whether the pores intersect or not; the number of voids in the molded body can be set by the porosity Po; the pore size is controlled by the pore radius R p.
In this embodiment, the process of generating a random three-dimensional pore rock mass model in a three-dimensional solid model through a discrete pore network is: randomly generating a space coordinate point in the three-dimensional solid model, and generating a pore based on the acquired pore shape and pore size characteristics by taking the space coordinate point as a center; the pore generation process is repeated to generate a network of pores, as shown in fig. 4.
(4) In the three-dimensional solid model, eliminating cracks at the overlapping part of the cracks and the pores, and constructing a random three-dimensional crack-pore coupling network;
specifically, in connection with fig. 5, the process of removing the fissures in the portion overlapping the aperture is as follows:
Random pores and pores are generated in a solid model with the same size, and the starting point coordinates are (0, 0);
defining a three-dimensional entity model as a domain 1, defining a random fracture network model as a domain 2, and defining a random pore network model as a domain 3;
Respectively storing the entity of the domain 1, the crack in the domain 2 and the pore in the domain 3 as three sets of space coordinate data;
determining intersection boundary surfaces of the domains 2 and 3 through space coordinate data overlap detection, and defining the intersection boundary surfaces as a domain 4;
Recording a dataset (x i,yi,zi) in the fracture network beyond the boundary surface of the intersection, defined as domain 5;
The domain 5 portion (x m-i,ym-i,zm-i) of the domain 2 dataset is eliminated, i.e., the fracture network after the intra-pore fracture is eliminated, defined as domain 6.
The process of constructing the random three-dimensional fracture-pore coupling network is as follows;
the elimination of the domain 4 portion (x n-i,yn-i,zn-i) of the domain 3 dataset, i.e., the elimination of the pore network at the intersection boundary, is defined as domain 7.
And carrying out contour merging on the positions of the intersecting boundary surfaces of the domain 6 and the domain 7 to obtain a coupled random three-dimensional fracture-pore coupling network.
(5) And performing gridding treatment on the random three-dimensional fracture-pore coupling network, and establishing a hybrid physical numerical model to realize rock mechanics analysis and simulation.
In this embodiment, the numerical gridding method of the random data body may be a three-sided gridding dividing method or a four-sided gridding dividing method, and the randomly generated fissures and pores are subjected to gridding encryption processing; the fissures and the pores generally have different dimensions and geometric characteristics, and the detail characteristics of the fissures and the pores can be better captured by respectively gridding, so that the simulation precision is improved. The grid encryption with different degrees is carried out on the areas with different scales and characteristics, so that the calculation amount can be reduced while the precision is kept, and the simulation efficiency is improved. Cracks and pores may have a local influence on the seepage or mechanical properties. The separate encryption grids may better capture local effects, ensuring a more accurate simulation of the changes and effects of these regions, without ignoring or fading these critical effects in the overall simulation. The separate meshing encryption of the fissures and the pores is also beneficial to subsequent processing and analysis. The results at different lattice densities can provide more information of different dimensions, helping to understand the complexity of the porous media in depth.
In addition, the data body can be cut at any angle through the cutting subdivision of the data grid body, and the cutting surface is mapped to a two-dimensional plane for numerical simulation of a two-dimensional fracture-pore network; according to the subsequent analysis requirement, the endpoint coordinate data of the grid points, the coordinates of the three-dimensional cracks and the central points of the pores in the grid and the permeability of the three-dimensional cracks and the central points of the pores can be derived in batches and stored.
Wherein the mesh encryption processing of the randomly generated cracks or pores comprises:
Performing gridding treatment on the coupling network by adopting a three-side grid or four-side grid method; random cracks or pores in the field 6 or the field 7 are selected respectively, and the maximum value and the minimum value of the unit grid quantity range are changed respectively, so that the grid encryption of the single cracks or the single pores can be realized.
In this embodiment, different flow equations are defined for the fields 6 and 7 according to different physical processes of the fissures and the pores; by setting model parameters and verifying, simulation calculation of the hybrid physical model can be performed.
Specifically, assuming that the underground aquifer, particularly the karst, has fissures as the main water guide channels, and the pores have water storage as the main;
according to the physical process of seepage of the fluid in the fracture, the flow equation of the definition domain 6 is as follows:
;
according to the physical process of seepage of the fluid in the pores, the flow equation of the domain 7 is defined as follows:
;
according to the physical process of seepage of fluid in the fracture and pore coupling network, the flow equation of the definition domain 4 is as follows:
;
Wherein, in the flow equation, K is the permeability coefficient; h is pore water head; h f is a slit head; c is a proportionality constant; s s f is the crack water storage rate; s s is the pore water storage rate.
As a specific embodiment, model shape and size data are defined, the model shape is set to be a rectangular parallelepiped, and the model size data are defined as follows: the length L c of the model was 10mm, the width W c was 5mm, and the height H c was 5mm. Meanwhile, the lower left corner of the model is set as the origin coordinate, i.e., (0, 0).
Setting fracture characteristics, namely setting 5 groups of fractures in total, wherein the average trend T a of each group of fractures is respectively set to be 90 degrees, 0 degrees, 30 degrees, 45 degrees, 60 degrees, 15 degrees and 75 degrees, the average inclination angle D a is respectively set to be 0 degrees, 90 degrees, 60 degrees, 10 degrees, 50 degrees, 30 degrees and 72 degrees, and the standard deviation of trend and trend is 0; the crack distribution is randomly generated through a normal distribution rule, and the crack length is characterized through an average trace length I avg and is respectively set to be 6.05mm, 5.35mm, 8.72mm, 10.86mm, 8.04mm, 12.60mm and 9.43mm; the number of cracks N c was set to 30 in total. FIG. 3 is a random three-dimensional fracture network generated from the above data.
Setting pore characteristics, and selecting small spherical pores; for clarity of the drawing, the number of pores was set to 300; the internal radius R p of the aperture was set to 0.05mm and the radius standard deviation was set to 0.05. Fig. 4 is a multi-scale random three-dimensional pore network generated from the above data.
And converting the generated random fracture network and the random pore network into vector data, taking a union set by using the lower left corner coordinates (0, 0) for coupling, and removing the fractures in the pores to realize the construction of the coupling network of the random three-dimensional fracture and the pores. FIG. 6 is a graph of a resulting multi-scale random three-dimensional fracture-pore coupled network.
The generated multiscale random three-dimensional fracture-pore coupling network is subjected to numerical grid division through a free tetrahedral grid, as shown in fig. 7, and the randomly generated fracture and pore are subjected to grid encryption treatment; in addition, horizontal cutting dissection was performed on the position of the data grid body h=3.6 cm, and the cut surface was mapped to a two-dimensional plane. Fig. 8 is a mapping of a cut surface on a two-dimensional plane.
The embodiment couples and digitizes the random three-dimensional fracture network and the random three-dimensional fracture network, can quickly, accurately and efficiently construct a random three-dimensional fracture-pore coupling network, further digitizes a data body, performs arbitrary angle cutting and lattice point information storage on the data body, and provides a model foundation for aspects such as rock-soil body seepage research, mechanics simulation, pollution transmission prediction, resource exploration and development, decision support and the like, for example:
① Seepage simulation: the model can more accurately simulate seepage behavior in the porous medium, and takes effects of mutual seepage, pressure transmission and the like between cracks and pores into consideration. This has important applications in groundwater flow, fluid migration in hydrocarbon reservoirs, and the like.
② Rock mechanics simulation: the three-dimensional fracture-pore coupling network model can simulate deformation and stress transmission in the rock and predict the influence of the fracture on the mechanical properties of the rock. The method has important significance for geotechnical engineering in engineering, underground storage design and the like.
③ Pollution transmission prediction: for environmental engineering and groundwater resource management, the model can help predict the transmission and diffusion of pollutants in an underground medium, and the influence of cracks and pores is considered, so that the influence range and the propagation path of the pollutants can be accurately evaluated.
④ Resource exploration and development: in the field of oil, gas and other resource exploration, the model is helpful for better understanding the characteristics of underground reservoirs and fluid migration rules and guiding the development and utilization of resources.
⑤ Optimizing design and decision support: by simulating the coupling network of the cracks and the pores, more accurate data can be provided, and support is provided for decisions such as engineering design, environment management, risk assessment and the like.
Example two
In one or more embodiments, a three-dimensional fracture-pore coupled network simulation generation system is disclosed, comprising:
the three-dimensional entity model construction module is used for constructing a three-dimensional entity model based on actual geological conditions;
The random three-dimensional fracture rock mass model construction module is used for acquiring fracture characteristics in actual geological conditions and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; wherein the fracture characteristics include at least fracture occurrence, fracture length, and fracture number;
The random three-dimensional pore rock mass model construction module is used for acquiring pore characteristics in actual geological conditions and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; the pore characteristics include at least pore shape, porosity, and pore size;
the random three-dimensional crack-pore coupling network module is used for removing cracks at the overlapping part of the cracks and the pores in the three-dimensional solid model to construct a random three-dimensional crack-pore coupling network;
And the model processing module is used for carrying out gridding processing on the random three-dimensional fracture-pore coupling network and establishing a hybrid physical numerical model so as to realize rock mechanics analysis and simulation.
The specific implementation manner of each module is the same as that in the first embodiment, and will not be described in detail.
Example III
In one or more embodiments, a terminal device is disclosed that includes a server including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the three-dimensional fracture-void coupled network simulation generation method of embodiment one when executing the program. For brevity, the description is omitted here.
It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general purpose processors, digital signal processors DSP, application specific integrated circuits ASIC, off-the-shelf programmable gate array FPGA or other programmable logic device, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include read only memory and random access memory and provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory. For example, the memory may also store information of the device type.
In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software.
Example IV
In one or more embodiments, a computer-readable storage medium is disclosed, in which are stored a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the three-dimensional fracture-void coupled network simulation generation method of embodiment one.
While the foregoing description of the embodiments of the present invention has been presented in conjunction with the drawings, it should be understood that it is not intended to limit the scope of the invention, but rather, it is intended to cover all modifications or variations within the scope of the invention as defined by the claims of the present invention.
Claims (8)
1. The three-dimensional fracture-pore coupling network simulation generation method is characterized by comprising the following steps of:
Constructing a three-dimensional entity model based on actual geological conditions;
Acquiring fracture characteristics in actual geological conditions, and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; wherein the fracture characteristics include at least fracture occurrence, fracture length, and fracture number;
acquiring pore characteristics in actual geological conditions, and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; the pore characteristics include at least pore shape, porosity, and pore size;
In the three-dimensional solid model, eliminating cracks at the overlapping part of the cracks and the pores, and constructing a random three-dimensional crack-pore coupling network;
The specific process for removing the cracks at the overlapping part of the holes is as follows: defining a three-dimensional entity model as a domain 1, defining a random three-dimensional fracture rock mass model as a domain 2, and defining a random three-dimensional pore rock mass model as a domain 3; storing the fissures in the domain 2 and the pores in the domain 3 as a space coordinate data set respectively; determining overlapping coordinates of the domain 2 and the domain 3 through overlapping detection of space coordinate data, determining an intersection boundary surface of the overlapping coordinates, and defining the intersection boundary surface as a domain 4; recording a data set exceeding the boundary surface of the intersecting line in the fracture network, and defining the data set as a domain 5; eliminating a domain 5 part in the domain 2 data set to obtain a crack of which the part overlapped with the pore is removed, and defining the crack as a domain 6;
The construction of the random three-dimensional fracture-pore coupling network comprises the following specific processes: eliminating a domain 4 part in the domain 3 data set to obtain a pore network with intersection boundary surfaces removed, and defining the pore network as a domain 7; carrying out contour merging on the positions of the intersecting boundary surfaces of the domain 6 and the domain 7 to obtain a random three-dimensional fracture-pore coupling network;
and performing gridding treatment on the random three-dimensional fracture-pore coupling network, and establishing a hybrid physical numerical model to realize rock mechanics analysis and simulation.
2. The method for generating the three-dimensional fracture-pore coupling network simulation according to claim 1, wherein a random three-dimensional fracture rock mass model is generated in the three-dimensional solid model through a discrete fracture network, and the method comprises the following specific steps of:
randomly generating a space coordinate point in the three-dimensional solid model, and generating a crack based on the acquired crack occurrence and crack length characteristics by taking the space coordinate point as a center; repeating the fracture generation process, thereby generating a fracture network; the fracture occurrence comprises average trend and average inclination angle of the fracture, and the fracture length is characterized by minimum trace length, maximum trace length or trace length average value of the fracture.
3. The method for generating the three-dimensional fracture-pore coupling network simulation according to claim 1, wherein the random three-dimensional pore rock mass model is generated in the three-dimensional solid model through a discrete pore network, and the method comprises the following specific steps of:
Randomly generating a space coordinate point in the three-dimensional solid model, and generating a pore based on the acquired pore shape and pore size characteristics by taking the space coordinate point as a center; and determining the number of pores through the porosity characteristics, and repeating the pore generation process to generate a pore network.
4. The method for generating the three-dimensional fracture-pore coupling network simulation according to claim 1, wherein the method for generating the three-dimensional fracture-pore coupling network simulation is characterized by establishing a hybrid physical numerical model, and specifically comprises the following steps:
according to the physical process of seepage of the fluid in the fracture, the flow equation of the definition domain 6 is as follows:
;
according to the physical process of seepage of the fluid in the pores, the flow equation of the domain 7 is defined as follows:
;
according to the physical process of seepage of fluid in the fracture and pore coupling network, the flow equation of the definition domain 4 is as follows:
;
Wherein, 、/>、/>The permeability coefficients in the x, y and z directions are respectively; h is pore water head; h f is a slit head; c is a proportionality constant; s s f is the crack water storage rate; s s is pore water storage rate; t is time.
5. The method for generating the three-dimensional fracture-pore coupling network simulation according to claim 1, wherein the random three-dimensional fracture-pore coupling network is subjected to gridding treatment, and the method is specifically as follows:
And performing gridding treatment on the coupling network by adopting a three-sided grid or four-sided grid method, and performing grid encryption treatment on randomly generated cracks or pores by adjusting the number range of all cracks or pore unit grids.
6. A three-dimensional fracture-void coupled network simulation generation system, comprising:
the three-dimensional entity model construction module is used for constructing a three-dimensional entity model based on actual geological conditions;
The random three-dimensional fracture rock mass model construction module is used for acquiring fracture characteristics in actual geological conditions and generating a random three-dimensional fracture rock mass model in the three-dimensional solid model through a discrete fracture network; wherein the fracture characteristics include at least fracture occurrence, fracture length, and fracture number;
The random three-dimensional pore rock mass model construction module is used for acquiring pore characteristics in actual geological conditions and generating a random three-dimensional pore rock mass model in the three-dimensional solid model through a discrete pore network; the pore characteristics include at least pore shape, porosity, and pore size;
The random three-dimensional crack-pore coupling network module is used for removing cracks at the overlapping part of the cracks and the pores in the three-dimensional solid model to construct a random three-dimensional crack-pore coupling network; the specific process for removing the cracks at the overlapping part of the holes is as follows: defining a three-dimensional entity model as a domain 1, defining a random three-dimensional fracture rock mass model as a domain 2, and defining a random three-dimensional pore rock mass model as a domain 3; storing the fissures in the domain 2 and the pores in the domain 3 as a space coordinate data set respectively; determining overlapping coordinates of the domain 2 and the domain 3 through overlapping detection of space coordinate data, determining an intersection boundary surface of the overlapping coordinates, and defining the intersection boundary surface as a domain 4; recording a data set exceeding the boundary surface of the intersecting line in the fracture network, and defining the data set as a domain 5; eliminating a domain 5 part in the domain 2 data set to obtain a crack of which the part overlapped with the pore is removed, and defining the crack as a domain 6; the construction of the random three-dimensional fracture-pore coupling network comprises the following specific processes: eliminating a domain 4 part in the domain 3 data set to obtain a pore network with intersection boundary surfaces removed, and defining the pore network as a domain 7; carrying out contour merging on the positions of the intersecting boundary surfaces of the domain 6 and the domain 7 to obtain a random three-dimensional fracture-pore coupling network;
And the model processing module is used for carrying out gridding processing on the random three-dimensional fracture-pore coupling network and establishing a hybrid physical numerical model so as to realize rock mechanics analysis and simulation.
7. A terminal device comprising a processor and a memory, the processor for implementing instructions; a memory for storing a plurality of instructions adapted to be loaded by a processor and to perform the three-dimensional fracture-void coupled network simulation generation method of any of claims 1-5.
8. A computer readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the three-dimensional fracture-void coupled network simulation generation method of any of claims 1-5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410189990.1A CN117745979B (en) | 2024-02-21 | 2024-02-21 | Three-dimensional fracture-pore coupling network simulation generation method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410189990.1A CN117745979B (en) | 2024-02-21 | 2024-02-21 | Three-dimensional fracture-pore coupling network simulation generation method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117745979A CN117745979A (en) | 2024-03-22 |
CN117745979B true CN117745979B (en) | 2024-05-03 |
Family
ID=90251213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410189990.1A Active CN117745979B (en) | 2024-02-21 | 2024-02-21 | Three-dimensional fracture-pore coupling network simulation generation method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117745979B (en) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787220A (en) * | 2016-04-22 | 2016-07-20 | 山东科技大学 | Coal bed high-pressure water injection fracturing-flow seeping value simulation method |
CN106127777A (en) * | 2016-06-27 | 2016-11-16 | 中山大学 | A kind of three dimensions crack separation identification and characterizing method |
CN109887083A (en) * | 2019-01-29 | 2019-06-14 | 中国石油集团测井有限公司西南分公司 | A kind of method for building up of Fracture-Pore dual media coupling penetration rate model |
CN110853138A (en) * | 2019-11-21 | 2020-02-28 | 科吉思石油技术咨询(北京)有限公司 | Dual-medium carbonate rock pore-fracture dual-network model construction method |
WO2020215524A1 (en) * | 2019-04-24 | 2020-10-29 | 山东科技大学 | Method for predicting coal porosity and permeability parameters based on fractal theory and ct scanning |
CN112132966A (en) * | 2020-09-29 | 2020-12-25 | 成都理工大学 | Shale fracture network connectivity characterization method based on topological structure |
CN113075728A (en) * | 2021-02-26 | 2021-07-06 | 河海大学 | Method for establishing compact sandstone multi-scale three-dimensional rock physical drawing board |
WO2022011894A1 (en) * | 2020-07-15 | 2022-01-20 | 中海油田服务股份有限公司 | Convolutional neural network-based modeling method and device for pore network model |
CN115524275A (en) * | 2022-10-18 | 2022-12-27 | 中国水利水电科学研究院 | Fractured rock mass permeability tensor determination method considering rock mass seepage and fracture seepage |
CN115937467A (en) * | 2023-03-13 | 2023-04-07 | 山东科技大学 | Method and system for dividing random three-dimensional fracture in upscale model grid |
CN116882221A (en) * | 2023-09-08 | 2023-10-13 | 山东科技大学 | Geothermal exploitation numerical simulation method and system based on three-dimensional fracture type thermal storage model |
JP7418767B1 (en) * | 2022-11-02 | 2024-01-22 | 浙大城市学院 | Three-dimensional fracture network rock model characterization method with multiscale heterogeneity |
CN117557874A (en) * | 2023-10-16 | 2024-02-13 | 重庆大学 | Intelligent prediction method for coal rock fatigue deformation fracture expansion |
-
2024
- 2024-02-21 CN CN202410189990.1A patent/CN117745979B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787220A (en) * | 2016-04-22 | 2016-07-20 | 山东科技大学 | Coal bed high-pressure water injection fracturing-flow seeping value simulation method |
CN106127777A (en) * | 2016-06-27 | 2016-11-16 | 中山大学 | A kind of three dimensions crack separation identification and characterizing method |
CN109887083A (en) * | 2019-01-29 | 2019-06-14 | 中国石油集团测井有限公司西南分公司 | A kind of method for building up of Fracture-Pore dual media coupling penetration rate model |
WO2020215524A1 (en) * | 2019-04-24 | 2020-10-29 | 山东科技大学 | Method for predicting coal porosity and permeability parameters based on fractal theory and ct scanning |
CN110853138A (en) * | 2019-11-21 | 2020-02-28 | 科吉思石油技术咨询(北京)有限公司 | Dual-medium carbonate rock pore-fracture dual-network model construction method |
WO2022011894A1 (en) * | 2020-07-15 | 2022-01-20 | 中海油田服务股份有限公司 | Convolutional neural network-based modeling method and device for pore network model |
CN112132966A (en) * | 2020-09-29 | 2020-12-25 | 成都理工大学 | Shale fracture network connectivity characterization method based on topological structure |
CN113075728A (en) * | 2021-02-26 | 2021-07-06 | 河海大学 | Method for establishing compact sandstone multi-scale three-dimensional rock physical drawing board |
CN115524275A (en) * | 2022-10-18 | 2022-12-27 | 中国水利水电科学研究院 | Fractured rock mass permeability tensor determination method considering rock mass seepage and fracture seepage |
JP7418767B1 (en) * | 2022-11-02 | 2024-01-22 | 浙大城市学院 | Three-dimensional fracture network rock model characterization method with multiscale heterogeneity |
CN115937467A (en) * | 2023-03-13 | 2023-04-07 | 山东科技大学 | Method and system for dividing random three-dimensional fracture in upscale model grid |
CN116882221A (en) * | 2023-09-08 | 2023-10-13 | 山东科技大学 | Geothermal exploitation numerical simulation method and system based on three-dimensional fracture type thermal storage model |
CN117557874A (en) * | 2023-10-16 | 2024-02-13 | 重庆大学 | Intelligent prediction method for coal rock fatigue deformation fracture expansion |
Non-Patent Citations (5)
Title |
---|
一种新的基于拓扑结构特征的微裂隙-孔隙空间描述方法;王鑫;蒋泽云;姚军;张琦;张国浩;张臻;段利亚;;科学通报;20180228(第Z1期);全文 * |
工程科技Ⅰ辑.2023,全文. * |
工程科技Ⅱ辑.2023,全文. * |
张元.随机溶孔-裂隙网络生成及其渗透特性研究.基础科学 * |
郭家繁.孔隙—裂隙双重介质渗流特性与模型研究.基础科学 * |
Also Published As
Publication number | Publication date |
---|---|
CN117745979A (en) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Makedonska et al. | Particle tracking approach for transport in three-dimensional discrete fracture networks: Particle tracking in 3-D DFNs | |
Ledoux | On the validation of solids represented with the international standards for geographic information | |
Xu et al. | A connectivity index for discrete fracture networks | |
CN106886682B (en) | Random walking particle tracking method for solute transport numerical simulation in single fracture | |
CN107655908A (en) | Method and device for constructing digital core | |
CN108921945B (en) | Pore network model construction method combining centering axis and solid model | |
CN115937467B (en) | Method and system for dividing random three-dimensional fracture in upscale model grid | |
CN116030218A (en) | Tetrahedral mesh dividing method, device, system and storage medium | |
Davy et al. | DFN, why, how and what for, concepts, theories and issues | |
Yao et al. | A Novel Numerical Model for Fluid Flow in 3D Fractured Porous Media Based on an Equivalent Matrix‐Fracture Network | |
Krotz et al. | Variable resolution Poisson-disk sampling for meshing discrete fracture networks | |
Hashemi Beni et al. | Toward 3D spatial dynamic field simulation within GIS using kinetic Voronoi diagram and Delaunay tetrahedralization | |
CN114065593A (en) | Method and device for underground water coupling simulation in multi-scale fracture matrix system | |
CN117745979B (en) | Three-dimensional fracture-pore coupling network simulation generation method and system | |
CN107886573B (en) | Slope three-dimensional finite element grid generation method under complex geological conditions | |
Fu et al. | Spatial topology identification of three-dimensional complex block system of rock masses | |
Ren et al. | Methods for the modelling of concrete mesostructures: a critical review | |
CN116740285A (en) | Gridding method after image reconstruction | |
Tu et al. | Comparative investigation of parallel spatial interpolation algorithms for building large-scale digital elevation models | |
CN115908733A (en) | Real-time segmentation and three-dimensional visualization method for angular point grid geological model | |
Lv et al. | The application of a complex composite fractal interpolation algorithm in the seabed terrain simulation | |
Jia et al. | Modeling of complex geological body and computation of geomagnetic anomaly | |
CN114187423A (en) | Surrounding rock fracture reconstruction method, electronic equipment and storage medium for three-dimensional simulation test | |
Kang et al. | A precise modeling method of three‐dimensional discrete fracture network based on rectangular joint model | |
CN116882255B (en) | Method and system for randomly generating porous medium model based on Fourier series |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |