CN111753414B - Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law - Google Patents
Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law Download PDFInfo
- Publication number
- CN111753414B CN111753414B CN202010546917.7A CN202010546917A CN111753414B CN 111753414 B CN111753414 B CN 111753414B CN 202010546917 A CN202010546917 A CN 202010546917A CN 111753414 B CN111753414 B CN 111753414B
- Authority
- CN
- China
- Prior art keywords
- model
- seepage
- continuous medium
- grid
- permeability
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses a seepage characteristic simulation method for a rock progressive destruction process based on the law of Darcy, which comprises the following steps: establishing a particle model and a continuous medium model which correspond to each other; carrying out secondary development on a seepage module in the continuous medium model to enable the seepage module to reflect the non-Darcy seepage characteristics; calculating the non-Darcy seepage characteristics of the rock at different damage stages in the continuous medium model and introducing the pore water pressure of each grid into the particle model; calculating the mesoscopic failure condition of the rock in the particle model to update the permeability in the continuous medium model; and repeating the processes of seepage characteristic calculation, pore water pressure introduction, particle model deformation calculation and permeability updating until the particle model is destroyed. Through the macro-microscopic bidirectional coupling alternate calculation, the non-Darcy seepage characteristics in the rock progressive destruction process are finally obtained. The method considers the local permeability change caused by the crack development in the rock progressive destruction process, and can accurately reflect the seepage characteristics in the rock progressive destruction process.
Description
Technical Field
The invention relates to a seepage characteristic simulation method, in particular to a seepage characteristic simulation method for a rock progressive failure process based on the law of Darcy, and belongs to the field of rock mechanics.
Background
The continuous accumulation of local microscopic damage often leads to the instability of the overall macroscopic deformation of the engineering rock mass. Under the condition of seepage, the local damage of the engineering rock mass can cause the local permeability characteristic to change obviously, the change can greatly promote the evolution process of the local damage, and the overall stability of the engineering rock mass is seriously threatened. Therefore, if the influence of local damage is neglected in the design of the rock mass engineering and only the macroscopic mechanical response is considered, the reliability of the design is greatly reduced.
In the progressive destruction process of the rock, the local permeability is obviously increased due to the development and the expansion of the fracture, so that the seepage characteristics of the rock are obviously heterogeneous and anisotropic, a more obvious non-Darcy flow effect is caused, and the mechanical properties of the rock are obviously influenced. Therefore, the seepage mechanical characteristics in the rock progressive failure process need to be researched, the seepage mechanical responses at different failure stages are revealed, and support is provided for the design of rock engineering.
At present, the influence of the development of the microscopic fracture on the seepage characteristics is less considered in the numerical research on the rock seepage mechanics, the research on the seepage is mainly based on Darcy's law rather than the more accurate non-Darcy's law, and the research result is difficult to accurately express the seepage mechanics characteristics in the rock progressive destruction process.
Disclosure of Invention
The invention aims to: aiming at the problems that the influence of the development of the microscopic fractures is less considered in the rock seepage characteristic simulation at present and the simulation method is mainly based on Darcy's law, the invention provides a seepage characteristic simulation method for the rock progressive destruction process based on the non-Darcy's law.
The technical scheme is as follows: the invention relates to a seepage characteristic simulation method of a rock progressive destruction process based on the non-Darcy law, which comprises the following steps:
(1) Establishing a rock particle model and a continuous medium model with the same size, extracting grids of the continuous medium model, introducing the grids into the particle model, and grouping particles of the continuous medium model;
(2) Introducing equivalent permeability to carry out secondary development on a seepage module of the continuous medium model, so that the seepage module can reflect the non-Darcy seepage characteristics;
(3) Calculating the non-Darcy seepage in the continuous medium model until the seepage is balanced, and recording the non-Darcy seepage characteristics of the continuous medium model at different rock damage stages; seepage characteristics comprise pore water pressure, seepage speed and the like;
(4) Extracting the pore water pressure of each grid in the continuous medium model, introducing the pore water pressure into the particle model, establishing a seepage topological structure of the particle model, and applying the pore water pressure with corresponding size to the particles in each grid;
(5) Carrying out deformation calculation under corresponding pore water pressure by using a particle model, determining the development and expansion conditions of the pore fracture at the stage, and then calculating and deriving the permeability of each deformed grid;
(6) Calling a continuous medium model in the particle model, and updating the permeability of the corresponding grid in the continuous medium model by using the permeability calculated by the particle model;
(7) Repeating the steps (3) to (6) until the particle model is destroyed;
(8) After the particle model is damaged, the step (3) is executed again, and seepage characteristics of the model during damage are recorded;
(9) And (4) sorting the seepage characteristics of each deformation stage of the particle model, and completing the seepage characteristic simulation of the rock progressive destruction process based on the law of Darcy.
Specifically, in the step (1), the method of extracting the mesh of the continuous medium model, introducing the mesh into the particle model, and grouping the particles comprises the following steps:
the continuous medium model is subjected to mesh division, vertex coordinates of each mesh are extracted and are derived in a Geometry format, each Geometry occupies a Geometry set, the name of the set is the same as the serial number of the meshes in the continuous medium model, then the Geometry is led into the particle model, particles in the particle model are grouped, and the grouping name is the same as the name of the Geometry set.
In the step (2), the method for performing secondary development on the seepage module of the continuous medium model by introducing the equivalent permeability may be:
converting the formula Forchheimer to the form of Darcy formula:
wherein k is e Equivalent permeability;
and (4) carrying out secondary development on the seepage module by using the above formula, and updating the permeability of each unit by using the continuous medium model in real time.
Preferably, the step (4) specifically includes the following steps (4.1) to (4.3):
(4.1) sequentially traversing the grids in the continuous medium model, and deriving the pore water pressure and the grid serial number of each grid into a pore water pressure file in the traversing process;
(4.2) reading the information in the pore water pressure file in the particle model into an array;
and (4.3) establishing a seepage topological structure of the particle model by using a 'tube domain method', determining the pore water pressure applied to the particles according to the grids where the particles are located, wherein the direction of the pore water pressure is the connecting line direction of the domain and the center of the particles.
The step (5) may specifically include the following steps (5.1) to (5.5):
(5.1) carrying out deformation calculation on the particle model;
(5.2) establishing inscribed circles of each grid as measurement circles of the porosity by taking the center of each grid as the center of a circle, and directly extracting the porosity measured by each measurement circle after deformation calculation is finished;
(5.3) grouping the cracks according to the positions of the cracks relative to the grids, then carrying out statistics on the cracks in each grid, and determining the crack development condition in each grid;
(5.4) calculating the permeability of each grid according to the porosity and fracture development conditions in the grids;
(5.5) exporting the group name and permeability of each mesh to a permeability file.
In the step (6), the method for updating the permeability of the corresponding grid in the continuous medium model comprises the following steps: and reading the permeability of each grid derived by the particle model in the continuous medium model, and updating the permeability of each grid of the continuous medium model.
In the step (9), the process of sorting the seepage characteristics of the particle model in each deformation stage is as follows: and (3) corresponding the strain of each deformation stage of the particle model to the seepage characteristic of the continuous medium model when the seepage is balanced under the strain one by one.
Has the advantages that: compared with the prior art, the invention has the advantages that: the simulation method of the invention can express the non-Darcy effect of seepage by carrying out secondary development on the seepage module in the continuous medium model, then the non-Darcy seepage characteristics of the rock at different damage stages are calculated in the continuous medium model and led into the particle model, the mesoscopic damage condition of the rock is obtained in the particle model, the permeability of each grid in the continuous medium model is updated, and the non-Darcy seepage characteristics in the progressive damage process of the rock are finally obtained by the macroscopic and mesoscopic bidirectional coupling alternate calculation of the continuous medium model and the particle model. The method considers the local permeability change caused by the crack development in the rock progressive destruction process, so that the non-Darcy seepage characteristics of the rock can accurately reflect the influence of the rock progressive destruction.
Drawings
FIG. 1 is a flow chart of a method for simulating seepage characteristics of a rock progressive failure process based on the law of Darcy of the invention;
FIG. 2 is a model of a particle and a model of a continuous medium used in the examples;
FIG. 3 is a diagram illustrating a case where particle models are grouped using a continuous medium model mesh in the embodiment;
FIG. 4 is a percolation topology in an embodiment;
FIG. 5 shows the arrangement of measurement circles in the example;
FIG. 6 is an initial water pressure distribution of the continuous medium model in the example;
FIG. 7 is the fracture development of the particle model at a deformation stage in the example;
FIG. 8 is a permeability distribution of a continuous medium model at a deformation stage in the embodiment;
FIG. 9 is the pore water pressure distribution of the continuous medium model in a certain deformation stage in the example.
Detailed Description
The technical solution of the present invention is further explained with reference to the drawings and the embodiments.
Referring to fig. 1, the method for simulating the non-darcy seepage characteristics based on the rock progressive failure process of the invention comprises the following steps:
(1) Establishing a rock particle model and a continuous medium model with the same size, extracting grids of the continuous medium model, introducing the grids into the particle model, and grouping particles of the grids;
specifically, in the step (1), the method of extracting the mesh of the continuous medium model, introducing the mesh into the particle model, and grouping the particles comprises the following steps:
the continuous medium model is subjected to mesh division, vertex coordinates of each mesh are extracted and are derived in a Geometry format, each Geometry occupies a Geometry set, the name of the set is the same as the serial number of the meshes in the continuous medium model, then the Geometry is led into the particle model, particles in the particle model are grouped, and the grouping name is the same as the name of the Geometry set.
(2) Introducing equivalent permeability to carry out secondary development on a seepage module of the continuous medium model, so that the seepage module can reflect the non-Darcy seepage characteristics;
specifically, the method for carrying out secondary development on the seepage module of the continuous medium model by introducing the equivalent permeability in the step (2) comprises the following steps:
converting the formula Forchheimer to the form of Darcy formula:
wherein k is e Equivalent permeability;
the seepage module is developed for the second time by using the above formula, so that the continuous medium model can update the permeability of each unit in real time.
(3) And performing non-Darcy seepage calculation in the continuous medium model until the seepage reaches balance, and recording the seepage characteristics of the continuous medium model at the stage, wherein the seepage characteristics comprise pore water pressure, seepage speed and the like.
(4) Extracting the pore water pressure of each grid in the continuous medium model, introducing the pore water pressure into the particle model, establishing the seepage topological structure of the particle model, and applying the pore water pressure with corresponding size to the particles in each grid.
The method specifically comprises the following steps (4.1) - (4.3):
(4.1) sequentially traversing the grids in the continuous medium model, and deriving the pore water pressure and the grid serial number of each grid into a pore water pressure file in the traversing process;
(4.2) reading the information in the pore water pressure file in the particle model into an array;
and (4.3) establishing a seepage topological structure of the particle model by using a 'tube domain method', determining the pore water pressure applied to the particles according to the grids where the particles are located, wherein the direction of the pore water pressure is the connecting line direction of the domain and the center of the particles.
(5) And (3) carrying out deformation calculation under corresponding pore water pressure by using the particle model, determining the development and expansion conditions of the pore fracture at the stage, and then calculating and deriving the permeability of each grid after deformation.
The method comprises the following steps (5.1) - (5.5):
(5.1) carrying out deformation calculation on the particle model;
(5.2) establishing inscribed circles of each grid as measurement circles of the porosity by taking the center of each grid as the center of a circle, and directly extracting the porosity measured by each measurement circle after deformation calculation is finished;
(5.3) grouping the cracks according to the positions of the cracks relative to the grids, then carrying out statistics on the cracks in each grid, and determining the crack development condition in each grid;
(5.4) calculating the permeability of each grid according to the porosity and fracture development conditions in the grids;
(5.5) exporting the group name and permeability of each mesh to a permeability file.
(6) And calling the continuous medium model in the particle model, and updating the permeability of the corresponding grid in the continuous medium model by using the permeability calculated by the particle model.
Specifically, the continuous medium model is called in the particle model, and then the permeability file derived from the particle model is read in the continuous medium model, so that the permeability of each grid of the continuous medium model is updated.
(7) Repeating the steps (3) to (6) until the particle model is destroyed;
(8) After the particle model is damaged, the step (3) is executed again, and seepage characteristics of the model during damage are recorded;
(9) And (4) sorting the seepage characteristics of each deformation stage of the particle model, and completing the seepage characteristic simulation of the rock progressive destruction process based on the law of Darcy.
The process of sorting out seepage characteristics of each deformation stage of the particle model comprises the following steps: and (3) corresponding the strain of each deformation stage of the particle model to the seepage characteristic of the continuous medium model when the seepage is balanced under the strain one by one.
Examples
The seepage characteristic simulation method of the rock progressive failure process based on the non-Darcy law is explained by taking a simulated seepage mechanical test as an example.
1. Respectively establishing a particle model and a continuous medium model, wherein the models are standard rock samples with the same size, and the established particle model and the established continuous medium model are shown in figure 2.
2. Extracting grids in the continuous medium model into the particle model, and grouping the particles in the particle model, wherein the grouping name is the same as the grid serial number, and the grouping condition is shown in fig. 3.
3. And (3) establishing a particle model seepage topological structure, wherein the seepage topological structure is shown in a figure 4.
4. And establishing an inscribed circle as a measurement circle of the porosity by taking the center of each grid as the center of the circle, wherein the established measurement circle is as shown in figure 5.
5. And carrying out secondary development on the seepage module of the continuous medium model by using the equivalent permeability, wherein the equation used in the secondary development is as follows:
wherein k is e Is the equivalent permeability.
6. And (3) carrying out seepage calculation by using the continuous medium model to obtain seepage characteristics in the model, wherein the seepage characteristics comprise pore water pressure distribution, seepage speed and the like, and the pore water pressure distribution in the initial stage is shown in figure 6.
7. And introducing the pore water pressure of each grid in the continuous medium model into the particle model, applying the pore water pressure to each grouped particle according to the specified size, and then performing deformation calculation on the particle model at a certain time step to obtain the development and expansion conditions of the pore fracture in the model.
8. And calculating the permeability of each grid of the particle model after the deformation calculation is finished according to the development and expansion conditions of the pore fracture in the particle model, and then exporting the permeability and the grid sequence number of each grid to a file.
9. And calling the continuous medium model in the particle model, and updating the permeability of the corresponding grid in the continuous medium model by using the permeability calculated by the particle model.
10. And (5) circularly executing 6-9 until the calculation of the particle model reaches a termination condition, performing seepage calculation in the continuous medium model during the damage again, and finishing the calculation result to complete the simulation of the non-Darcy seepage characteristics based on the rock progressive damage process. The fracture development at a certain stage is shown in figure 7, the permeability distribution at the stage is shown in figure 8, and the pore water pressure distribution is shown in figure 9.
Claims (1)
1. A seepage characteristic simulation method for a rock progressive failure process based on the non-Darcy law is characterized by comprising the following steps:
(1) Establishing a rock particle model and a continuous medium model with the same size, extracting grids of the continuous medium model, introducing the grids into the particle model, and grouping particles of the continuous medium model; the process is as follows: dividing meshes of the continuous medium model, extracting vertex coordinates of each mesh and deriving the coordinates in a Geometry format, wherein each Geometry occupies a Geometry set, the name of the set is the same as the serial number of the meshes in the continuous medium model, introducing the Geometry into the particle model, and grouping particles in the particle model, wherein the grouping name is the same as the name of the Geometry set;
(2) Introducing equivalent permeability to carry out secondary development on a seepage module of the continuous medium model, so that the seepage module can reflect the non-Darcy seepage characteristics; the process is as follows: converting the formula Forchheimer to the form of Darcy formula:
wherein k is e Equivalent permeability;
carrying out secondary development on the seepage module by using the above formula, and updating the permeability of each unit by using the continuous medium model in real time;
(3) Calculating the non-Darcy seepage in the continuous medium model until the seepage reaches balance, and recording the non-Darcy seepage characteristics of the continuous medium model at different damage stages of the rock;
(4) Extracting the pore water pressure of each grid in the continuous medium model, introducing the pore water pressure into the particle model, establishing a seepage topological structure of the particle model, and applying the pore water pressure with corresponding size to the particles in each grid; the method comprises the following steps (4.1) - (4.3):
(4.1) sequentially traversing the grids in the continuous medium model, and deriving the pore water pressure and the grid serial number of each grid into a pore water pressure file in the traversing process;
(4.2) reading the information in the pore water pressure file in the particle model into an array;
(4.3) establishing a seepage topological structure of the particle model by using a 'tube domain method', determining the pore water pressure applied to the particles according to the grids where the particles are located, wherein the direction of the pore water pressure is the connecting line direction of the domain and the center of the particles;
(5) Carrying out deformation calculation under corresponding pore water pressure by using a particle model, determining the development and expansion conditions of the pore fracture at the stage, and then calculating and deriving the permeability of each deformed grid; comprises the following steps (5.1) to (5.5):
(5.1) carrying out deformation calculation on the particle model;
(5.2) establishing inscribed circles of each grid as measurement circles of the porosity by taking the center of each grid as the center of a circle, and directly extracting the porosity measured by each measurement circle after deformation calculation is finished;
(5.3) grouping the fractures according to the positions of the fractures relative to the grids, then counting the fractures in each grid, and determining the fracture development condition in each grid;
(5.4) calculating the permeability of each grid according to the porosity and fracture development conditions in the grids;
(5.5) exporting the group name and permeability of each mesh into a permeability file;
(6) Calling a continuous medium model in the particle model, and updating the permeability of the corresponding grid in the continuous medium model by using the permeability calculated by the particle model, wherein the process is as follows: reading the permeability of each grid derived from the particle model in the continuous medium model, and updating the permeability of each grid of the continuous medium model;
(7) Repeating the steps (3) to (6) until the particle model is destroyed;
(8) After the particle model is damaged, the step (3) is executed again, and seepage characteristics of the model during damage are recorded;
(9) And (3) sorting seepage characteristics of the particle model in each deformation stage, and enabling the strain magnitude of each deformation stage of the particle model to correspond to the seepage characteristics of the continuous medium model in seepage balance under the strain one by one, so as to complete seepage characteristic simulation of the rock progressive destruction process based on the non-Darcy law.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010546917.7A CN111753414B (en) | 2020-06-16 | 2020-06-16 | Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010546917.7A CN111753414B (en) | 2020-06-16 | 2020-06-16 | Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111753414A CN111753414A (en) | 2020-10-09 |
| CN111753414B true CN111753414B (en) | 2023-03-24 |
Family
ID=72675362
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010546917.7A Active CN111753414B (en) | 2020-06-16 | 2020-06-16 | Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111753414B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112362556B (en) * | 2020-11-13 | 2024-03-29 | 重庆大学 | Method for obtaining continuous function of permeability coefficient of coal mine mining stable region |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104021277A (en) * | 2014-05-14 | 2014-09-03 | 河海大学 | Numerical analysis method for piping phenomenon |
| CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
| CN110361312A (en) * | 2019-07-05 | 2019-10-22 | 河海大学 | The determination method of permeability and porosity relationship during rock seepage liquefaction |
-
2020
- 2020-06-16 CN CN202010546917.7A patent/CN111753414B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104021277A (en) * | 2014-05-14 | 2014-09-03 | 河海大学 | Numerical analysis method for piping phenomenon |
| CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
| CN110361312A (en) * | 2019-07-05 | 2019-10-22 | 河海大学 | The determination method of permeability and porosity relationship during rock seepage liquefaction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111753414A (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111338300B (en) | Physical simulation method and system of production line based on digital twins | |
| CN110704888B (en) | Unconventional oil and gas reservoir encrypted well volume fracturing construction parameter optimization design method | |
| CN106988739B (en) | Shale reservoir fracturing fracture identification and interpretation evaluation method | |
| CN102081030B (en) | Geotechnical model test system based on macro-micro mechanics as well as elaboration test method | |
| US11067970B1 (en) | Method for designing production line based on digital twin | |
| RU2010122059A (en) | SIMULATION OF THE BREAK OF THE COLLECTOR | |
| CN111753414B (en) | Seepage characteristic simulation method for rock progressive destruction process based on non-Darcy law | |
| CN115983062B (en) | High arch dam seismic damage assessment method and system based on finite element model correction | |
| CN112326524A (en) | Rock pore permeability measurement method based on CT scanning image | |
| CN111123374A (en) | Ground penetrating radar full-waveform inversion method based on matched filtering | |
| CN104134013A (en) | Wind turbine blade modal analysis method | |
| CN101701935B (en) | Refractory brick inside quality online detection device | |
| CN114324580A (en) | Intelligent knocking detection method and system for structural defects | |
| CN103574291A (en) | Tube bursting positioning method based on artificial immunity system | |
| CN106934729B (en) | Building detection and identification method and device | |
| CN115470675A (en) | Prediction method and system for crack propagation path of heavy-load pressure vessel | |
| CN110441204B (en) | Digital core simulation-based compact reservoir fracturing fluid damage digital evaluation method | |
| CN112001055B (en) | Low-amplitude thin oil reservoir water content prediction method based on microstructure | |
| Sladojevic et al. | Investigation of the influence of aerodynamic coupling on response levels of mistuned bladed discs with weak structural coupling | |
| CN109992841B (en) | Blasting block space omnibearing partition coupling numerical simulation method | |
| CN112558179A (en) | Intelligent inversion method of far field stress state | |
| CN111474905A (en) | Parameter drift fault diagnosis method in manufacturing process of electromechanical product | |
| CN113886992B (en) | Digital twin modeling method based on multi-fidelity data | |
| CN117421939B (en) | Shale oil fracture system simulation agent method based on track piecewise linearization | |
| CN109472083B (en) | Civil engineering structure damage detection method based on statistical process control |
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 |