CN117688743A - Porous medium generation method - Google Patents
Porous medium generation method Download PDFInfo
- Publication number
- CN117688743A CN117688743A CN202311682469.3A CN202311682469A CN117688743A CN 117688743 A CN117688743 A CN 117688743A CN 202311682469 A CN202311682469 A CN 202311682469A CN 117688743 A CN117688743 A CN 117688743A
- Authority
- CN
- China
- Prior art keywords
- growth
- cluster
- probability
- core
- initial
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000012010 growth Effects 0.000 claims abstract description 82
- 239000011435 rock Substances 0.000 claims abstract description 26
- 238000009826 distribution Methods 0.000 claims abstract description 19
- 230000003698 anagen phase Effects 0.000 claims abstract description 6
- 239000011148 porous material Substances 0.000 abstract description 18
- 238000002910 structure generation Methods 0.000 abstract description 8
- 238000012800 visualization Methods 0.000 description 4
- 239000002245 particle Substances 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a porous medium generation method, which comprises the following steps: (1) Determining initial growth core distribution probability C d Porosity of rock V r And a rock gauge number kn, determining a rock initial core number kk= (kn) C in two dimensions d Initial growth probability P for each growth core with eight different growth directions i Each of the kk initial growth cores is a cluster, each cluster having a different code and a different probability of additional growth P c The total growth probability per cluster was p=p i *P c The method comprises the steps of carrying out a first treatment on the surface of the (2) The eight directions of the initial growth core are given a random probability number D which is (0, 1) and obeys normal distribution when the initial growth core starts growing i (3) repeating step (2) until the volume fraction of the growth phase reaches a given value, i.e., a given volume fraction V r The method comprises the steps of carrying out a first treatment on the surface of the (4) And finally generating a set of CSGS for the cluster structure when the growth stops. The invention randomly grows in four parametersThe improved cluster structure generation set based on the algorithm more truly simulates the pore structure of the rock.
Description
Technical Field
The invention provides a porous medium generation method, and belongs to the technical field of porous media.
Background
Because of the complexity and randomness of porous media, it is very difficult, if not impossible, to fully describe their microstructure. We can only obtain some statistical-based average information such as average porosity or pore size distribution. In fact, a more accurate prediction of porous media transport characteristics requires a more detailed description of the morphology of the porous media as a whole, including geometric characteristics (such as particle or pore shape), volume, and topological properties (such as pore deformation and connectivity). Random growth of the media is the simplest method of constructing an artificial porous media.
The prior art is the four parameter random growth algorithm (QSGS) (Wang, m., et al, mesoscopic predictions of the effective thermal conductivity for microscale random porous media. Phys Rev E Stat Nonlin Soft Matter Phys,2007.75 (3 pt 2): p.036702), which comprises the following steps:
(1) Setting a growth core distribution probability C d The first growth cores are randomly defined in the growth grid with values not greater than the volume fraction of the growth phase. Each growth point in the grid will be assigned a random number by a uniform distribution function within (0, 1), each random number not greater than C d Will be selected as a growth core.
(2) Each growth core has eight different growth directions on the two-dimensional plane, according to the growth probability D of each given growth direction i (i represents the direction), each growth direction of each growth nucleus of the growth stage is given a random probability number of (0, 1) if the random probability number of the direction is larger than the growth probability D thereof i The point where the direction connects with the growth nucleus becomes a new growth nucleus.
(3) Repeating the process (2) until the volume fraction of the growth phase reaches a given value, i.e. a given volume fraction P 2 (i.e. its porosity epsilon in the rock).
(4) Finally, when growth stops, the set of quad structure generation sets (QSTS) is obtained.
The rock pores in the nature are disordered and often have no regularity, but the rock pores simulated by the prior art scheme can find the regularity, and the growth probability D in each direction i Under the same condition, the rock pores grown by the algorithm are more regular and do not accord with the condition of real rock pores. Probability of growth D in a certain direction i When the probability of growth is much larger than that of other, the rock pores grown by the algorithm are extremely high, more rock grows in the direction with high probability, and the porosity is smaller in the direction. The other direction is the opposite. In summary, the method of this technique does not simulate the pore structure of real rock well.
Disclosure of Invention
The invention aims to improve the technology, adds a concept of growing clusters into a four-parameter random growth algorithm (QSPS), and names a Cluster Structure Generation Set (CSGS), so that the improved algorithm can simulate the relatively real rock pore structure.
The specific technical scheme is as follows:
a porous medium generation method comprising the steps of:
(1) Determining initial growth core distribution probability C d Porosity of rock V r And a rock gauge number kn, determining a rock initial core number kk= (kn) C in two dimensions d Each growth core has eight different growth directions, the initial growth probability P of these eight directions i (i represents the direction) is the same, each growth core in the kk initial growth cores is a cluster, each cluster has different codes and different additional growth probabilities P c I.e. the total growth probability per cluster is p=p i *P c ；
(2) The eight directions of the initial growth core are given a random probability number D which is (0, 1) and obeys normal distribution when the initial growth core starts growing i If the total cluster growth probability P corresponding to the growth nucleus is larger than the random probability D i The direction and the growthThe point of nuclear connection becomes a new growth nucleus, and the growth nucleus still belongs to the growth cluster, and the cluster is endowed with independent codes;
(3) Repeating the step (2) until the volume fraction of the growth phase reaches a given value, i.e. a given volume fraction V r (i.e., its porosity epsilon in the rock);
(4) And finally generating a set of CSGS for the cluster structure when the growth stops.
The method can effectively and truly simulate the pore structure of the rock, and can obviously show that the pore structure simulated by the visualization of the previous quadruple structure generation set (QGS) is more realistic in the visualization of the Cluster Structure Generation Set (CSGS), and can find that the pore structure generated by logarithmic distribution is more realistic than that generated by normal distribution in the visualization of the Cluster Structure Generation Set (CSGS), so that the improved Cluster Structure Generation Set (CSGS) based on a four-parameter random growth algorithm (QGS) more truly simulates the pore structure of the rock.
Drawings
FIG. 1 is a CSGS flow chart of the invention;
FIG. 2 is a set of CSGS with an embodiment having kn of 250 and a normal distribution of cluster probabilities;
FIG. 3 is a set of CSGS with an embodiment having a kn of 250 and a cluster probability of lognormal distribution;
FIG. 4 is a set of CSGS with an embodiment having a kn of 350 and a normal distribution of cluster probabilities;
FIG. 5 is a CSGS set with an embodiment having a kn of 350 and a cluster probability of lognormal distribution;
FIG. 6 is a set of CSGS with an embodiment having a kn of 350 and a normal distribution of cluster probabilities;
FIG. 7 is a CSGS set with an embodiment having a kn of 350 and a cluster probability of lognormal distribution;
fig. 8 is a three-dimensional CSGS set with kn of 200 and a normal distribution of cluster probabilities for an example.
Detailed Description
A porous medium generation method comprising the steps of:
(1) Determining initial growth core distribution probability C d Porosity of rock V r And a rock gauge number kn, determining a rock initial core number kk= (kn) C in two dimensions d Each growth core has eight different growth directions, the initial growth probability P of these eight directions i (i represents the direction) is the same, each growth core in the kk initial growth cores is a cluster, each cluster has different codes and different additional growth probabilities P c I.e. the total growth probability per cluster is p=p i *P c ；
(2) The eight directions of the initial growth core are given a random probability number D which is (0, 1) and obeys normal distribution when the initial growth core starts growing i If the total cluster growth probability P corresponding to the growth nucleus is larger than the random probability D i The point where the direction connects with the growth nuclei becomes a new growth nuclei and the growth nuclei remain in the growth cluster giving the cluster an independent code;
(3) Repeating the step (2) until the volume fraction of the growth phase reaches a given value, i.e. a given volume fraction V r (i.e., its porosity epsilon in the rock);
(4) And finally generating a set of CSGS for the cluster structure when the growth stops.
For the above algorithm of Cluster Structure Generation Set (CSGS), the present invention establishes an algorithm flow chart, as shown in fig. 1, and performs visualization of simulated pore structure by matlab procedure, and in addition, the present embodiment also performs logarithmic random cluster probability, as shown in fig. 1 to 7, and performs upgrading of three-dimensional pore structure based on two-dimensional plane pore structure, as shown in fig. 8.
Claims (2)
1. A method of generating a porous medium, comprising the steps of:
(1) Determining initial growth core distribution probability C d Porosity of rock V r And a rock gauge number kn, determining a rock initial core number kk= (kn) C in two dimensions d Each of the raw materialsThe long cores all have eight different growth directions, the initial growth probability P of these eight directions i All the same, i represents the direction, each growth core in the kk initial growth cores is a cluster, each cluster has different codes and different additional growth probabilities P c I.e. the total growth probability per cluster is p=p i *P c ；
(2) The eight directions of the initial growth core are given a random probability number D which is (0, 1) and obeys normal distribution when the initial growth core starts growing i ；
(3) Repeating the step (2) until the volume fraction of the growth phase reaches a given value, i.e. a given volume fraction V r ；
(4) And finally generating a set of CSGS for the cluster structure when the growth stops.
2. The method according to claim 1, wherein in the step (2), if the total cluster growth probability P corresponding to the growth nuclei is larger than the random probability D thereof i The point where the direction connects with the growth nuclei becomes a new growth nuclei and the growth nuclei remain in the growth cluster giving the cluster an independent code.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311682469.3A CN117688743A (en) | 2023-12-08 | 2023-12-08 | Porous medium generation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311682469.3A CN117688743A (en) | 2023-12-08 | 2023-12-08 | Porous medium generation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117688743A true CN117688743A (en) | 2024-03-12 |
Family
ID=90133020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311682469.3A Pending CN117688743A (en) | 2023-12-08 | 2023-12-08 | Porous medium generation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117688743A (en) |
-
2023
- 2023-12-08 CN CN202311682469.3A patent/CN117688743A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Solé et al. | Are rainforests self-organized in a critical state? | |
| KR20120095869A (en) | Method and apparatus for encoding a mesh model, encoded mesh model, and method and apparatus for decoding a mesh model | |
| CN105499575A (en) | Design and manufacturing method of porous grid structure material | |
| CN107680131B (en) | Method for rapidly determining volume size of porous medium characterization unit | |
| Bang Pham et al. | Rapid freeze prototyping technique in bio‐plotters for tissue scaffold fabrication | |
| CN117688743A (en) | Porous medium generation method | |
| CN113191044B (en) | Topological optimization design method of single-material porous structure | |
| Popovych et al. | Loss of coherence in a system of globally coupled maps | |
| Machlus et al. | Correlated rigidity percolation in fractal lattices | |
| Zhao et al. | Percolation laws of a fractal fracture-pore double medium | |
| Steiner et al. | A cellular model for the evolutionary development of lightweight material with an inner structure | |
| CN115859402A (en) | Plant fiber mulching film structure simulation method and system based on improved four-parameter random growth method | |
| CN104485970B (en) | The building method of solid size rate, the pattern matrix of multi code Rate of Chinese character QC LDPC codes | |
| Masood et al. | Development of scaffold building units and assembly for tissue engineering using fused deposition modelling | |
| Li et al. | Coupling control of pore size and spatial distribution in bone scaffolds based on a random strategy for additive manufacturing | |
| CN111122270A (en) | Silicate high-porosity artificial core and preparation method thereof | |
| CN115630461A (en) | Randomized broken stone particle generation method considering three-dimensional broken stone particle form boundary | |
| Van Segbroeck et al. | Stochastic simulation of the chemoton | |
| CN113539374A (en) | Method, device, medium and apparatus for generating protein sequence of high-thermal-stability enzyme | |
| CN115114787B (en) | Complex organization mode characterization method for reservoir fractal pore structure | |
| CN113177334A (en) | Solid finite element modeling method and system for multi-cell TPMS porous structure | |
| CN115544836B (en) | Optimization method for evolution of fluid-solid common regulation structure | |
| CN111310777A (en) | Method and system for acquiring target category number in K-means algorithm | |
| Chandler et al. | Realistic computer models of aerogels | |
| Tunbridge et al. | An L‐systems approach to the modelling of fungal growth |
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 |