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 simulation system for porous permeable rocks.
Background
The geothermal resource is a clean renewable energy source with large reserve and good stability, and has great significance for energy conservation, emission reduction, global warming response and haze treatment. However, due to the uncertainty of the distribution of geothermal resources and the migration of underground water, before the development and utilization of geothermal resources, the local underground heat flow and water flow conditions need to be effectively analyzed and predicted so as to improve the survey accuracy and reduce the development risk, thereby reducing the time investment and the capital investment of the whole process of the development and utilization of geothermal resources.
The numerical simulation is one of important means for predicting underground heat flow and water flow, can be used for preliminarily testing the local geological condition before large-scale exploration, and the simulation degree of the constructed numerical model is one of key factors for determining the effectiveness of the testing condition. Two properties exist for underground rock: one is permeability, i.e. the fluid has some mobility in the rock; the other is porosity, i.e. solid particles of a certain size can pass through the rock. However, at present, there are two main ways of dealing with underground rock in modeling: one is to reduce the rock model to a model with certain porosity and permeability, but only through fluid and not through solid particles; the other is to simplify the rock model into a solid compact model with a series of pores, and the solid compact model is a solid body without porosity and permeability except the pores, and any fluid and solid particles cannot pass through the solid body. It can be seen that the above two models can only characterize a single aspect of the underground liquid-solid physical action, and still have certain defects in the simulation degree.
Disclosure of Invention
In order to reflect the migration characteristics of underground fluid more truly and fully represent the physical properties of the underground rock, namely porosity (certain solid particles can pass through) and permeability (fluid can pass through at a certain speed), a certain degree of fusion and improvement on the existing model are necessary.
In order to achieve the purpose, the invention adopts the technical scheme that:
the invention provides a simulation method of porous permeable rock, which comprises the following steps:
step 1, collecting rock samples and fluid samples of rock stratums corresponding to a research area;
step 2, obtaining the porosity, permeability and 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 according to the rock model and the fluid model to obtain a porous permeable rock model.
Preferably, in step 2, the porosity, permeability and average pore diameter of the rock sample obtained in step 1 are obtained by the following specific method:
carrying out pore permeability test analysis on the rock sample obtained in the step 1 to obtain the porosity and permeability of the rock sample;
and (3) acquiring a pore size distribution characteristic curve of the rock sample by using a gas adsorption analysis method, and calculating the average pore size of the rock sample according to the pore size distribution characteristic curve.
Preferably, in step 2, the solid content ratio and the average particle size 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 fraction distribution rule of the solid particles in the fluid sample by using a solid particle size fraction analysis method, and calculating the average particle size of the solid particles in the fluid sample according to the size fraction 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 specific method is as follows:
the pore type rock model comprises a matrix, wherein a plurality of pore channels arranged in a matrix form are formed in the matrix, and the pore channels are communicated along any direction of the matrix;
setting the porosity of the matrix to be 0 according to the permeability and the porosity characteristics of the rock sample; the diameter of the pore canal is the average pore diameter of the rock sample;
calculating the pore channel spacing of the pore channel according to the porosity and the average pore diameter of the rock sample;
and calculating the permeability of the matrix according to the permeability of the rock sample.
Preferably, the porosity type rock model is a cubic structure of L × L × L, wherein L > 10d, and d is the average pore diameter of the rock sample.
Preferably, in step 3, the fluid model is constructed according to the 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 round ball models, the particle size distribution rule of each round ball 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 round ball models is consistent with the average particle size of the solid particles in the fluid sample;
the solid phase flow has the same proportion as the proportion of the solid components in the fluid sample;
the parameters of the liquid phase stream 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:
and introducing the rock model and the fluid model into fluid calculation software, and establishing a geometric model, wherein the geometric model is a porous permeable rock model.
A simulation system of porous rock is capable of operating the simulation method of porous rock, and comprises 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 rock strata corresponding to the research area;
the data processing unit is used for acquiring the porosity, the permeability and the average pore diameter of the rock sample; obtaining the solid component ratio and the average particle size of solid particles in the fluid sample;
the modeling unit is used for constructing a pore type rock model according to the 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:
the invention provides a simulation method of a porous permeable rock, which is characterized in that on the basis of determining a research object, a rock sample and a fluid sample of a rock stratum corresponding to a research area are collected, relevant parameters of the rock sample and the fluid sample are respectively obtained, model parameters such as pore canal spacing, matrix permeability and the like are determined according to the measured rock parameter and fluid parameter, a rock model and a fluid model are respectively constructed, and then a porous permeable rock model is obtained according to the rock model and the fluid model; the rock model related by the invention is a honeycomb rock model considering the matrix seepage characteristic and the pore channel blocking condition at the same time; the related fluid model comprises two different phase states of a liquid phase and a solid phase; meanwhile, the rock model and the fluid model are unified, so that the real migration state of underground water in the mining and irrigation process can be simulated in a high-simulation manner, data support and technical reserve are provided for further underground fluid research, and the method has better popularization significance in the field of exploration and development of the geothermal mining and irrigation analysis stage.
Drawings
FIG. 1 is a schematic structural view of the present invention;
FIG. 2 is a schematic flow diagram of the present invention.
Detailed Description
In order to make up for the lack of effective modeling methods for porous permeable rocks, the invention aims to provide a simulation method for porous permeable rocks. The method is mainly used for collecting rock samples and fluid samples of corresponding regions and rock stratums on the basis of determining a research object, wherein the rock samples and the fluid samples are required to meet relevant requirements. Secondly, carrying out pore permeability 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 using a gas adsorption analysis means; and (3) analyzing the components of the collected fluid sample, determining the solid component content a%, and determining the size fraction distribution rule of the solid particles by using a solid particle size fraction analyzing means. Thirdly, determining the pore canal spacing l and the matrix permeability k according to the measured rock parameter and the fluid parametermAnd (5) waiting for model parameters, and constructing a rock model and a fluid model. And finally, introducing the rock model and the fluid model into fluid calculation software for subsequent operation, and providing data support and technical reserve for the migration state analysis of the underground fluid under the actual geological condition.
Referring to fig. 1 and 2, embodiments of the present invention will be described in detail below with reference to the accompanying drawings and examples.
The invention discloses a simulation method of porous permeable rock, which comprises the following operation steps:
step 1, determining a research object, and collecting a rock sample and a fluid sample.
After a certain research area and rock stratum are determined as a research object, rock samples of the corresponding area and rock stratum are collected. When the sample is collected, the rock sample with the regional representation is selected as much as possible, and the overall levels of lithology, microstructure, physical parameters and the like in a certain region can be reflected. The rock sample should be at least 5cm x 10cm in size to allow for subsequent correlation test analysis. At the same time, a fluid sample should also be collected near the rock, in an amount of at least 100ml, and to ensure that no significant compositional contamination occurs during transport.
And 2, acquiring basic parameters of the rock and the surrounding fluid.
Performing a hole in a collected rock samplePerforming seepage test analysis to obtain the porosity phi and the permeability k of the rock sample; and use of N2/CO2And 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, analyzing the components of the collected fluid sample to determine the solid content in the fluid sample to be a%; and determining the size fraction distribution rule of the solid particles in the fluid sample by using a solid particle size fraction analysis means, and calculating the average particle size D of the solid particles in the fluid sample according to the rule. And the parameters are recorded in sequence, and a data base is laid for building a rock model and a fluid model.
And 3, establishing a rock model according to the rock parameters.
In order to better trace both the permeability and porosity of the rock, a "honeycomb" model can be created in which the cells 2 allow the free passage of fluid but to a certain extent limit the passage of solid particles in the fluid, while the matrix 1 has a certain percolation capacity allowing the passage of fluid according to darcy's law but preventing the passage of solid particles.
In the "honeycomb" model, the porosity is determined by the rock sample
And average pore diameter d, calculating the pore space l in the model:
In addition, in the "honeycomb" model, since two parts of the matrix 1 and the cell channels 2 are involved at the same time, the overall permeability k of the actual rock is converted into the permeability k of the matrixmThe 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, and k isfreeIs the free fluid permeability. From this, the permeability k of the matrix in the model can be derivedm。
And constructing a pore type rock model by using the obtained parameters. The rock model is an L multiplied by L cube (L is a certain length which is more than 10 d), and through holes are arranged in the cube along a certain direction, and the size of the holes is
The channel spacing is l. Porosity of matrix model
Is 0 and has a permeability of k
m(ii) a The pore canal is a blank model.
And 4, establishing a fluid model according to the fluid parameters.
Considering that two phases, namely a liquid phase and a solid phase, are involved in the fluid, both parts should be considered in the fluid model. The solid phase flow 4 is simplified into a series of spherical ball models, the particle size distribution rule of the solid phase flow is consistent with the analysis result of the solid particle size fraction, and the average particle size is D; the liquid phase stream 3 should then be in accordance with the relevant parameters of the measured fluid.
And 5, introducing the rock model and the fluid model into fluid calculation software for subsequent operation.
And (4) establishing a geometric model in software such as ANSYS based on the solid model of the porous permeable rock and the fluid model of the liquid-solid two-phase flow constructed in the steps 3 and 4, and calculating the motion state of the fluid in the porous permeable rock according to actual boundary conditions and action conditions 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 effective modeling method of the porous permeable rock, can provide a 'honeycomb' rock model which considers the matrix seepage characteristic and the pore blocking condition at the same time, can convert the actual rock parameters into 'honeycomb' rock model parameters, can also provide a fluid model which considers two different phase states of liquid phase, solid phase and the like, unifies the rock model and the fluid model, simulates the real migration state of the underground water in the process of mining and irrigating with high simulation, provides data support and technical reserve for further underground fluid research, and has better popularization significance in the field of exploration and development in the stage of geothermal mining and irrigating analysis.
The above description is only exemplary of the present invention and should not be construed as limiting the scope of the invention, so that the substitution of equivalent elements or the equivalent changes and modifications made in accordance with the scope of the present invention should be covered thereby.