Simulation method and system for porous permeable rock
Technical Field
The invention relates to the technical field of geothermal exploration and development, in particular to a simulation method and a system for porous permeable rock.
Background
The geothermal resource is a clean renewable energy source with large reserves and good stability, and has great significance for saving energy, reducing emission, coping with global warming and treating haze. However, due to uncertainty of geothermal resource distribution and groundwater migration conditions, prior to geothermal resource development and utilization, effective analysis and prediction of local groundwater heat flow and water flow conditions are required to improve survey accuracy and reduce development risk, thereby reducing time investment and capital investment for the whole process of geothermal development and utilization.
Numerical simulation is one of the important means for predicting subsurface heat flow and water flow, and can be used for preliminary bottoming of the local geological condition before large-scale exploration is not performed, and the simulation degree of the constructed numerical model is one of the key factors for determining whether bottoming condition is effective or not. Underground rock has two properties at the same time: one is permeability, i.e. the fluid has a certain mobility in the rock; the other is porosity, i.e. solid particles of a certain size can pass through the rock. However, at present, in modeling, there are two main ways of treating the subsurface rock: one is to reduce the rock model to a model with a certain porosity and permeability, but only through fluid, not through solid particles; the other is to simplify the rock model into a solid compact model with a series of pore channels, and the solid compact model is a solid body without porosity and permeability except the pore channels, and the solid body cannot pass through any fluid and solid particles. It can be seen that the above two models can only characterize a single aspect of the subsurface liquid-solid physical effect, and still have certain drawbacks in terms of simulation.
Disclosure of Invention
In order to more truly reflect the migration characteristics of underground fluid, fully characterize the physical properties of the underground rock, such as porosity (certain solid particles can pass through) and permeability (certain rate of fluid can pass through), and therefore the existing model needs to be fused and improved to a certain extent, the invention aims to provide a simulation method and system of the porous permeable rock, and the defect of low simulation precision of the action simulation of the geothermal fluid and the rock at present is overcome.
In order to achieve the above purpose, the invention adopts the following technical scheme:
the invention provides a simulation method of porous permeable rock, which comprises the following steps:
step 1, collecting rock samples and fluid samples of corresponding rock strata in a research area;
step 2, acquiring the porosity, the permeability and the average pore diameter of the rock sample obtained in the step 1; acquiring the solid component ratio and the average particle size of solid particles in the fluid sample obtained in the step 1;
step 3, constructing a pore type rock model according to the rock parameters obtained in the step 2; constructing a fluid model according to the fluid parameters;
and 4, constructing a porous permeable rock model according to the rock model and the fluid model.
Preferably, in step 2, the porosity, permeability and average pore size of the rock sample obtained in step 1 are obtained by the following specific methods:
performing pore-penetration test analysis on the rock sample obtained in the step 1 to obtain the porosity and the permeability of the rock sample;
and obtaining a pore size distribution characteristic curve of the rock sample by using a gas adsorption analysis method, and calculating according to the pore size distribution characteristic curve to obtain the average pore size of the rock sample.
Preferably, in step 2, the solid component ratio and the average particle diameter of the solid particles in the fluid sample obtained in step 1 are obtained by:
analyzing the components of the collected fluid sample to obtain the solid component ratio in the fluid sample;
and determining the size distribution rule of the solid particles in the fluid sample by using a solid particle size analysis method, and calculating the average particle size of the solid particles in the fluid sample according to the size distribution rule of the solid particles.
Preferably, in step 3, a pore type rock model is constructed according to the rock parameters obtained in step 2, and the concrete method is as follows:
the pore type rock model comprises a matrix, wherein a plurality of pore channels which are arranged in a matrix form are formed in the matrix, and the pore channels penetrate through the matrix along any direction;
setting the porosity of the matrix to 0 according to the permeability and the pore characteristics of the rock sample; the pore diameter is the average pore diameter of the rock sample;
calculating to obtain pore space of pore according to porosity and average pore diameter of the rock sample;
and calculating the permeability of the matrix according to the permeability of the rock sample.
Preferably, the pore rock model is a cubic structure of lxl×l, where L > 10d, d is the average pore size of the rock sample.
Preferably, in step 3, a fluid model is constructed according to fluid parameters, and the specific method is as follows:
the fluid model comprises a solid phase flow and a liquid phase flow, wherein the solid phase flow comprises a plurality of sphere models, the particle size distribution rule of each sphere model is consistent with the particle size distribution rule of solid particles in the fluid sample obtained in the step 2, and the average particle size of the plurality of sphere models is consistent with the average particle size of the solid particles in the fluid sample;
the solid phase flow has the same ratio as the solid component in the fluid sample;
the parameters of the liquid phase flow are consistent with the parameters of the fluid sample.
Preferably, in step 4, a porous permeable rock model is constructed according to the rock model and the fluid model, and the specific method is as follows:
the rock model and the fluid model are imported into fluid computing software to build a geometric model, which is a porous permeable rock model.
A simulation system of porous permeable rock, which can run the simulation method of porous permeable rock, comprising a data acquisition unit, a data processing unit and a modeling unit, wherein:
the data acquisition unit is used for acquiring rock samples and fluid samples of corresponding rock strata in the research area;
the data processing unit is used for acquiring the porosity, the permeability and the average pore diameter of the rock sample; acquiring the solid component ratio and the average particle size of solid particles in a fluid sample;
the modeling unit is used for constructing a pore type rock model according to rock parameters; constructing a fluid model according to the fluid parameters; and constructing a porous permeable rock model according to the rock model and the fluid model.
Compared with the prior art, the invention has the beneficial effects that:
according to the simulation method of the porous permeable rock, on the basis of determining a research object, rock samples and fluid samples corresponding to rock layers in a research area are collected, relevant parameters of the rock samples and the fluid samples are obtained respectively, model parameters such as pore space, matrix permeability and the like are determined according to the measured rock parameters and the measured fluid parameters, a rock model and a fluid model are built respectively, and then the porous permeable rock model is built according to the rock model and the fluid model; the rock model related by the invention is a honeycomb rock model which simultaneously considers the seepage characteristic of a matrix and the blocking condition of a pore canal; the fluid model comprises two different phases of liquid phase and solid phase; meanwhile, the rock model and the fluid model are unified, the real migration state of underground water in the process of mining and irrigation can be simulated in a high-simulation mode, data support and technical reserve are provided for further underground fluid research, and the method has good popularization significance in the field of exploration and development in the geothermal mining and irrigation analysis stage.
Drawings
FIG. 1 is a schematic diagram of the structure of the present invention;
fig. 2 is a schematic flow chart of the present invention.
Detailed Description
In order to compensate for the lack of an effective modeling method of the porous permeable rock, the invention aims to provide a simulation method of the porous permeable rock. The invention mainly collects rock samples and fluid samples of corresponding areas and rock strata on the basis of determining research objects, wherein the rock samples and the fluid samples meet related requirements. Secondly, carrying out pore permeation test analysis on the collected rock sample to obtain the porosity phi and the permeability k of the rock sample, and calculating to obtain the average pore diameter d of the rock by utilizing a gas adsorption analysis means; and (3) carrying out component analysis on the collected fluid sample, determining the solid component proportion a%, and determining the particle size distribution rule of the solid particles by using a solid particle size analysis means. Again, from the measured rock parameters and fluid parameters, the cell spacing l, matrix permeability k are determined m And model parameters, constructing a rock model and a fluid model. Finally, the rock model and the fluid model are led into fluid computing software to carry out subsequent operation, so as to analyze the migration state of the underground fluid under the actual geological conditionProviding data support and technical reserves.
Referring to fig. 1 and 2, embodiments of the present invention will be described in detail with reference to the drawings and examples.
The invention relates to a simulation method of porous permeable rock, which comprises the following operation steps:
and step 1, determining a research object, and collecting a rock sample and a fluid sample.
After a region of interest and a formation are determined to be the subject, rock samples are collected for the corresponding region and formation. When the sample is collected, the rock sample with the area representation is selected as much as possible, and the integral level of lithology, microstructure, physical parameters and the like in a certain area can be reflected. The rock sample should be at least 5cm by 10cm in size to allow for subsequent correlation test analysis. At the same time, a fluid sample should be collected near the rock, the fluid sample collection should be at least 100ml, and ensure that no significant contamination of components occurs during transport.
And 2, acquiring basic parameters of the rock and surrounding fluid.
Performing pore-penetration test analysis on the collected rock sample to obtain the porosity phi and the permeability k of the rock sample; and utilize N 2 /CO 2 And analyzing the pore size distribution characteristic curve of the rock by using an equal gas adsorption analysis means, and calculating the average pore size d of the rock according to the pore size distribution characteristic curve. Meanwhile, carrying out component analysis on the collected fluid sample to determine the solid component in the fluid sample to account for a percent of a%; and determining the particle size distribution rule of the solid particles in the fluid sample by using a solid particle size analysis means, and calculating the average particle size D of the solid particles in the fluid sample according to the particle size distribution rule. And sequentially recording the parameters, and laying a data foundation for building a rock model and a fluid model.
And 3, building a rock model according to the rock parameters.
To better describe both rock permeability and porosity, a "honeycomb" model may be created in which the channels 2 allow free passage of fluid but restrict to some extent the passage of solid particles in the fluid, while the matrix 1 has a certain percolation capacity allowing passage of fluid according to darcy's law but preventing passage of solid particles.
In the "honeycomb" model, according to the rock sample porosity
And (3) calculating the pore space/in the model by averaging the pore diameter d:
In addition, in the "honeycomb" model, since two parts of the matrix 1 and the channels 2 are involved at the same time, the overall permeability k of the actual rock needs to be converted into the permeability k of the matrix m The calculation formula is as follows:
wherein k is the permeability of the rock sample, l is the pore space, d is the average pore diameter of the rock sample, k free Is free fluid permeability. From this, the matrix permeability k in the model can be derived m 。
And constructing a pore type rock model by using the obtained parameters. The rock model is a cube with L multiplied by L (L is a certain length of more than 10 d), and through holes are arranged in the cube along a certain direction, and the size of the holes is
The spacing of the pore channels is l. Matrix model porosity->
0, permeability k
m The method comprises the steps of carrying out a first treatment on the surface of the The pore canal is a blank model.
And 4, establishing a fluid model according to the fluid parameters.
Since two phases, liquid and solid, are involved in the fluid, these two parts should be considered simultaneously in the fluid model. The solid phase flow 4 is simplified into a series of sphere models, the particle size distribution rule of the solid phase flow is consistent with the analysis result of the particle size of the solid particles, and the average particle size is D; the liquid phase stream 3 should then be consistent with the parameters associated with the measured fluid.
And 5, introducing the rock model and the fluid model into fluid computing software for subsequent operation.
Based on the solid model of the porous permeable rock and the fluid model of the liquid-solid two-phase flow constructed in the step 3 and the step 4, a geometric model is built in ANSYS and other software, and the motion state of the fluid in the porous permeable rock is calculated according to actual boundary conditions and action conditions so as to provide data support and technical reserve for the migration state of the underground fluid under actual geological conditions.
The invention makes up the lack of an effective modeling method of porous permeable rock, can provide a honeycomb rock model which simultaneously considers the matrix seepage characteristic and the pore canal blocking condition, can convert the actual rock parameters into the honeycomb rock model parameters, can also provide a fluid model which considers two different phases such as liquid phase and solid phase, unifies the rock model and the fluid model, simulates the real migration state of groundwater in the mining and filling process in a high-simulation manner, provides data support and technical reserve for further underground fluid research, and has better popularization significance in the exploration and development field of the geothermal mining and filling analysis stage.
The foregoing description of the embodiments of the invention should not be taken as limiting the scope of the invention, and therefore, all changes and modifications that come within the meaning and range of equivalency of the claims are intended to be embraced therein.