A kind of method of quantitative assessment reservoir diagenetic evolutionary process porosity Spatio-temporal Evolution
Technical field
The present invention relates to petroleum exploration field, and specifically, the present invention relates to a kind of quantitative assessment reservoir diagenetic to develop
The method of journey porosity Spatio-temporal Evolution.
Background technology
Petroleum resources has become the strategic resource closely bound up with national economy, social development and national security.Nothing
By the conventional gas and oil for being traditional, technology relative maturity, or the emerging, shale gas that has a high potential, coal bed gas etc. are non-
Conventional resource, equal preservation is in reservoir rock hole or crack.To effectively be exploited to it, it is necessary to clear and definite oil gas into
Rule is hidden, the physical property such as porosity and permeability of reservoir condition is correctly evaluated.
The distribution of petroleum resources and physical property are mainly by deposition and diagenetic joint effect and control:1. deposition is made
After suitable place is reached with the material for referring to be carried by moving medium, precipitation, the process of accumulation occurs since condition changes
Effect, deposition be influence reservoir quality " congenital " factor, decide the spatial distribution of reservoir and original physical property;2. into
Rock effect refers to after deposits to the whole process of rock bonding, continuous due to overlying sediments thing after deposits
Accumulation, deposit experienced different type and the diagenesis of varying strength during burying, formation, guarantor to reservoir pore space
Deposit, destroy and the physical property of reservoir plays particularly important effect.Therefore, diagenesis be influence reservoir quality " day after tomorrow " because
Element, decides the final physical property of reservoir, and the property such as the porosity of reservoir, permeability is the knot that diagenesis transforms it now
Fruit.
Diagenetic process and its complexity, according to typical diagenesis event, can be divided into multiple continuous periods, each period is equal
The conditions such as corresponding different buried depth, temperature, pressure, at the same also tectonic movement, external fluid intrusion, primary rock producing hydrocarbon into etc. it is a variety of
The influence of factor, spatially also shows as heterogeneous characteristic.The diagenesis mainly occurred in reservoir has compaction, glue
Knot effect, metasomasis, crystallization, eluviation, aquation and biological chemistry action etc..Whole diagenetic process is substantially all
Carried out in the environment for have water, in very long diagenetic process, complicated underground aqueous solution and country rock mineral Long Term Contact are inevitable
A series of physics, chemical reaction can occur, cause mineral that dissolving or precipitation occurs, trigger reservoir porosity to change.Therefore, urgently
Under the conditions for the treatment of diagenesis different in clear and definite diagenetic process, influence of the interior physics, chemical action occurred of reservoir to porosity,
Exactly recover the temporal-spatial erosion process of porosity during diagenesis, be oil gas so that the distribution of Favorable Reservoir can accurately be evaluated
Direction is found out in exploration, improves production efficiency.
On the research of diagenetic process porosity, on the one hand current technology scheme is to be based on the geologic informations such as well point,
By means such as theory analysis, petrology test and physical analogys, simulation sediment filling history, buried history, heat source Evolutionary History, recover
Calcite cementation history, silicious gangue history, the cementing history of clay mineral, spectra dissolving history etc., determine the Typical Diagenetic of different times
The change of effect, intensity and reservoir porosity, the methods of test based on petrology, determines that the original hole of reservoir retains and secondary hole
The situation of gap increase and decrease, belongs to qualitative or sxemiquantitative analysis, what is obtained is the integrated curved of porosity, is not had with time-evolution
Quantitative analysis and evaluation, fail accurately to provide the porosity concrete numerical value in some time or some region in diagenetic process.
On the other hand, there is a small amount of document to be based on geochemistry theoretical research diagenesis at present, with calculation of thermodynamics or
Numerical simulation means, evaluate transformation of a certain diagenesis to porosity, but these researchs are all to stick to single typical case
Diagenesis, such as dolomitization, illitization effect, belongs to mechanism study, can not predict actual diagenetic process mesoporous
The change of porosity.However, actual diagenesis includes multiple diagenesis, process and its complexity.At present, there has been no a kind of technology
Scheme can be based on the comprehensive knowledges such as geology, geochemistry and computer, quantitatively evaluate and predict reservoir in diagenetic process
Porosity with space-time differentiation.
The content of the invention
It is an object of the invention to provide a kind of side of quantitative assessment reservoir diagenetic evolutionary process porosity Spatio-temporal Evolution
Method;This method can utilize the means that experiment and numerical simulation are combined, and replicate reservoir diagenetic evolutionary process exactly, and quantitative
Recover differentiation curve of the porosity with the time.
The present invention is combined with numerical simulation technology by systematicness experiment, establish diagenesis evolutionary process porosity when
Empty evolvement simulation system, thus during the very long complicated diagenesis evolution of clear and definite reservoir reservoir porosity development law, it is quantitative
Evolving trend and distribution characteristics spatially of the reservoir porosity with the time are evaluated, discloses the Forming Mechanism and sky of Favorable Reservoir
Between Distribution Characteristics, which will point the direction for the exploration of petroleum resources, improve the production efficiency of oil gas.
Data of the technical program based on petrology test analysis etc. qualitatively, fully use geochemistry numerical simulation skill
Art, realizes the quantitative evaluation that porosity is developed with the time and spatially.
In order to achieve the above object, on the one hand, the present invention provides a kind of quantitative assessment reservoir diagenetic evolutionary process porosity
The method of Spatio-temporal Evolution, described method includes following steps:
(1) Diagenetic Sequence and porosity variation tendency of reservoir are determined, the diagenesis event of identification reservoir experience, chooses device to hole
The influential diagenesis event of porosity (such as organic acid filling, carbon dioxide intrusion, abnormal high temperature high pressure), delimit multiple typical cases
, continuous diagenesis period (each diagenesis period there are a variety of possible diagenesis conditions, the speed filled such as organic acid, filling
Amount etc.), establish one or more possible diagenesis evolution processes;
(2) the diagenesis evolution process speculated respectively in simulation steps (1), by the analog result of each diagenesis evolution process with
The Diagenetic Sequence and porosity variation tendency of reservoir reality are contrasted, and select analog result therein consistent with real data
Diagenesis drilled process, it is determined that the diagenesis condition in diagenesis evolution process, diagenesis period and each diagenesis period that reservoir is undergone;
(3) according to the diagenesis bar of identified diagenesis evolution process, diagenesis period and each diagenesis period in step (2)
Part, carries out the high temperature and pressure diagenesis laboratory experiment in different diagenesis periods, test the front and rear fluid of each diagenesis period experiment into
Point, rock sample component and pore structure;
(4) it is completely the same with the experiment condition in step (3), the diagenesis numerical simulation in different diagenesis periods is carried out,
The analog result in each diagenesis period is contrasted with experimental result, correction simulation Minerals dynamics calculation parameter and Bi Biao
Area;
(5) simulated time (5-8 days) in each diagenesis period in step (3), (4) is extended to the actual diagenesis of reservoir
Time scale (hundreds thousand of years to millions of years), based on corrected analog parameter in step (4), carries out actual diagenesis evolution mistake
The numerical simulation (millions of years) of journey, by the reservoir diagenetic sequence in analog result and step (1) and porosity variation tendency into
Row contrast, calibration model parameter, establishes quantitative simulation system of the porosity with time-evolution again;
(6) based on the geology in regional extent and diagenesis data, the quantitative mould with established porosity with time-evolution
Plan system, builds threedimensional model, and the evolutionary process of reservoir in simulated domain scope, contrasts with real data, finally establish hole
Degree realizes porosity spatially with the quantitative evaluation and prediction of time with the quantitative simulation system of Spatio-temporal Evolution.
Some specific embodiments according to the present invention, wherein, step (1) is to be based on petrology and geochemistry, with reference to rock
Sample thin section identification, rock sample X ray diffracting data, fluid inclusion and maturity analysis test, determine reservoir Diagenetic Sequence and
Porosity variation tendency, the authigenetic clay rim of identification reservoir experience, divides typical, continuous diagenesis period, establishes diagenesis
Evolutionary process.
Some specific embodiments according to the present invention, wherein, step (1) is by indoor petrology test analysis, fluid
Inclusion enclave test, isotope test, specify the material resource condition and thermal evolution history of reservoir, based on Mineral Transformation relation and different diagenesis
The test data of timing node, determines the Diagenetic Sequence and porosity variation tendency of reservoir, thus it is speculated that the diagenesis evolution process of reservoir,
And according to the influence degree to porosity, every kind of possible diagenesis evolution process is generalized as continuous several diagenesis periods.
Some specific embodiments according to the present invention, wherein, step (1) is influence degree of the basis to porosity, will be every
The possible diagenesis evolution process of kind is generalized as continuous 4-6 period.
Some specific embodiments according to the present invention, wherein, step (1) is when delimiting 4-6 typical, continuous diagenesis
Phase.
Some specific embodiments according to the present invention, wherein, step (1) is to establish the possible diagenesis evolution process of 3-5 kinds.
Some specific embodiments according to the present invention, wherein, step (2) is to use simulation softward, simulation by periods step (1)
Middle speculated diagenesis evolution process.
Some specific embodiments according to the present invention, wherein, the simulation softward of step (2) is the migration of multiphase flow reactive solute
Simulation softward.
The multiphase flow reactive solute transported simulation software can be presently commercially available conventional corresponding software, for example
TOUGHREACT, TOUCHSTONE etc..
Some specific embodiments according to the present invention, wherein, the analog result of step (2) includes rock mineral composition, ore deposit
Thing precipitates/dissolving trend and porosity in one or more.
Some specific embodiments according to the present invention, wherein, step (3) is filled with indoor diagenesis physical simulation experiment
Put, carry out the high-temperature high-voltage reaction experiment in different diagenesis periods, test the front and rear fluid composition of each period experiment, rock sample
Component and pore structure.
Some specific embodiments according to the present invention, wherein, each period experimental period of step (3) is 5-8 days, common 20-
48 days.
Some specific embodiments according to the present invention, wherein, step (3) is exemplary ion in measurement fluid pH value and fluid
Concentration.
Some specific embodiments according to the present invention, wherein, step (3) described exemplary ion includes K+、Ca2+、Na+、Mg2 +、Al3+, Fe (containing all valence states), Cl-, Si (full Si measure) and SO4 2-In one or more combinations.
Wherein it is understood that in the art, the concentration mensuration of Fe refers to the iron ion total amount of each valence state;Si's is dense
Degree measure refers to the measure of full Si in aqueous solution.
Some specific embodiments according to the present invention, wherein, step (3) test rock sample component includes carrying out rock sample
Following test:1. X-ray diffraction and spectrum analysis;2. scanning electron microscope and thin section identification;3. casting body flake;4. X-ray micron
CT。
Some specific embodiments according to the present invention, wherein, step (4) is to be based on the condition of laboratory experiment in step (3),
The numerical simulation in different diagenesis periods is carried out, analog result and experimental result are contrasted, corrects analog parameter.
Some specific embodiments according to the present invention, wherein, step (4) carries out the diagenesis numerical value in different diagenesis periods
Simulation, the temperature in each diagenesis period, pressure, diagenesis period corresponding with step (3) time it is consistent, will be each
The analog result in diagenesis period is contrasted with experimental result, correction simulation Minerals dynamics calculation parameter and specific surface area.
Some specific embodiments according to the present invention, wherein, the analog result of step (4) includes the opposite of each mineral and contains
Amount, the concentration of each ion and porosity.
In conclusion the present invention provides a kind of side of quantitative assessment reservoir diagenetic evolutionary process porosity Spatio-temporal Evolution
Method.The method of the present invention has the following advantages that:
(1) on diagenetic research it was in the past the petrology test analysis means that are based on mostly, it is qualitative or semi-quantitatively
Study Diagenetic Sequence and porosity develops.However, geochemistry process plays huge effect in diagenetic process, especially
There is direct influence on reservoir porosity, it might even be possible to say the evolution process for determining porosity, thus be accordingly used in the calculating earth
The thermodynamics and kinetics theory of chemical action is particularly important.Geochemistry numerical simulation technology has been applied to diagenesis by the present invention
In the research of process, with reference to means such as calculation of thermodynamics, laboratory experiment, petrology test analysis and natural analogies, it is mutually authenticated,
Differentiation curve of the porosity with the time, but also the geology diagenesis feature in the range of calmodulin binding domain CaM, evaluation and prediction are not only recovered
The Spatio-Temporal Change Characteristics of regional extent internal porosity, study diagenetic process porosity in terms of time-domain and spatial domain two
Characteristics of Evolution, thoroughly realizes the quantification of the evaluation of diagenetic process porosity and prediction, at home and abroad there has been no this research.
(2) diagenetic process and its complexity, are by multiple stages, multiple subprocess collective effects, influence each other between them,
Close ties, the research on diagenesis geochemistry numerical simulation in the past, is often confined to single diagenesis, such as white clouds
Petrifaction, illitization effect etc..The porosity that the technical program is considered under all diagenesis joint effects develops
Feature, it is consistent with the diagenetic process of actual reservoir, more direct, more valuable information can be provided for oil-gas exploration.
(3) the technical scheme is that making first with the diagenesis of comprehensive geochemistry synthesis in theory study of many complicated
With consistent with actual diagenesis.
(4) technical solution of the present invention, after preliminary generalization (first step), adds numerical simulation work, that is,
Indoors test before corrected repeatedly by numerical simulation means, determine diagenetic process, diagenesis period and diagenesis condition, so again into
Row laboratory test, success rate are high.This step is very crucial, considerably increases the success rate of the technical solution.
Brief description of the drawings
Fig. 1 is 1 process flow diagram flow chart of embodiment;
Fig. 2 is the schematic diagram of water-rock interaction model during 1 diagenesis evolution of embodiment;
Fig. 3 contrasts for 1 analog result of embodiment and experimental data;
Fig. 4 is differentiation schematic diagram of 1 porosity of embodiment with Geologic Time;
Fig. 5 is the space inhomogeneous distribution schematic diagram of 1 porosity of embodiment and typical mineral.
Embodiment
Below by way of the specific embodiment implementation process that the present invention will be described in detail and the beneficial effect produced, it is intended to which help is read
Reader more fully understand the present invention essence and feature, not as to this case can practical range restriction.
Embodiment 1
Research background
Northern Ordos gas field possesses huge petroleum resources amount, and Shihezi Formation is main force's gas pay.By very long multiple
Miscellaneous Geological Evolution, Shihezi Formation reservoir are subjected to a series of strong diagenesis, and average pore now is only 0.06,
Mean permeability is only 1mD, belongs to the tight sandstone reservoir of low porosity and low permeability, it is stronger non-that diagenetic transformation has reservoir
Homogenieity, causes the exploitation of compact sandstone gas to encounter the problem of " hypotonic, low pressure, low yield, low abundance ".At present, Shihezi Formation
There is dispute in the Forming Mechanism of tight sandstone reservoir, reservoir occurs fine and close time and process and do not know still, that is, porosity
Evolution process it is unclear, cause the production efficiency of compact sandstone gas low., can quantitative assessment diagenesis evolution by the technology of the present invention
Process porosity specifies the mechanism of reservoir densification, so as to be carried for the exploitation of compact sandstone gas with the differentiation in time and space
For instructing.Its detailed process is as follows:
1) tested by a series of indoor petrology test analysis, fluid inclusion test, isotope etc., it specify that stone
The Diagenetic Sequence and porosity variation tendency of box group reservoir, the test number based on Mineral Transformation relation and different diagenesis nodes
According to, thus it is speculated that the diagenesis evolution process of Shihezi Formation reservoir has 4 kinds of possibility, according to influence journey of the geochemistry process to porosity
Degree, 6 continuous typical periods (as shown in Figure 1) are generalized as by every kind of possible diagenesis evolution process:
1. compacting → early stage shallow burial diagenesis → abnormal high temperature high pressure → CO2The buried Tibetan of filling → organic acid filling → late period
Diagenesis;
2. compacting → early stage shallow burial diagenesis → CO2The buried Tibetan of filling → abnormal high temperature high pressure → organic acid filling → late period
Diagenesis;
3. compacting → early stage shallow burial diagenesis → organic acid filling → abnormal high temperature high pressure → CO2Filling → late period buried Tibetan
Diagenesis;
4. compacting → early stage shallow burial diagenesis → abnormal high temperature high pressure → organic acid filling → CO2Filling → late period buried Tibetan
Diagenesis;
Meanwhile the diagenesis condition in each diagenesis period also has a variety of possibility, such as organic acid, CO2The filling speed of filling, fill
Fluence and filling duration etc..
2) multiphase flow reactive solute transported simulation software TOUGHREACT, the 4 kinds of diagenesis evolutions more than simulation speculated are used
Process, every kind of diagenesis evolution process are divided into 6 period simulations, and each period corresponds to a variety of diagenesis condition (organic acid, CO2Fill
Filling speed, charging amount and the filling duration of note).By various analog results and the Diagenetic Sequence and porosity variation tendency of reservoir
Contrasted, determine the main geochemical reaction that Shihezi Formation reservoir is undergone, Mineral Transformation relation such as formula (1)-(4) institute
Show.After Shihezi Formation reservoir experienced initial compacting and early diagenesis buries, with hydrocarbon charge, a large amount of CO are experienced2Enter
Invade, the change to porosity is larger, subsequent stratum there occurs high-temperature and high-pressure conditions mutation, then with the second of oil gas
Filling, organic acid injection, finally experienced late period buried Tibetan diagenesis evolution process.In short, worked by the numerical simulation of this step,
Learn diagenesis evolution process that Shihezi Formation reservoir undergone with speculate in step 1) the 2. kind may be consistent, and other 3 kinds
The trend that possible diagenesis evolution process simulation results are showed with measured data in step 1) is inconsistent.Meanwhile pass through number
Value simulation work, it is determined that each diagenesis condition such as the temperature of diagenetic stage, pressure.
CaCO3(calcite)+H+→Ca2++HCO3 -
Feldspathoid+H2O+H+→(K+,Na+,Ca2+)+SiO2+ kaolinite
Kaolinite+K+→ illite+H2O+H+
Smectite+K++Al3+→ illite+Na++Ca2++Fe3++Mg2++Si4+
3) it is real in laboratory, the high-temperature high-voltage reaction that this 6 diagenesis periods have been carried out with Petrogenetic Simulation system experimental device
Test, first 8 days period, equal 6 days of other periods.Wherein, initial rock sample component and fluid concentrations are according to the research in step 1
Area's petrology test result is counter to be pushed away, and Tables 1 and 2 is that the fluid of initial rock sample component and diagenesis evolution process is dense respectively
Degree.The major parameter such as the temperature in each stage, pressure is taken from step 2 in experiment, as shown in table 3.
The relative amount of each mineral during 1 initial rock sample of table is prepared
The concentration (mg/L) of each fluid during table 2 is tested
The experiment condition in each stage during 3 diagenesis evolution of table
The rock sample and water sample after each diagenesis period reaction are collected, water sample is cooled to room temperature, filters, pH value determination and solution
Middle exemplary ion (K+、Ca2+、Na+、Mg2+、Al3+、Fe、Cl-, Si and SO4 2-) concentration;Following test has been carried out to rock sample:①X
X ray diffraction (XRD) and spectrum analysis (XRF), test rock sample Minerals and constituent content.2. scanning electron microscope (SEM) and thin slice mirror
It is fixed, observe corrosion or the Precipitation of rock surface mineral.3. casting body flake, measures and calculates the Areal porosity of rock sample.4. X-ray
Micron CT, scans the pore structure and face characteristic of rock sample.
4) condition based on laboratory experiment, carries out the numerical simulation of diagenesis evolution process (including different diagenesis periods), mould
Type schematic diagram is as shown in Fig. 2, analog result and experimental result are contrasted, as shown in figure 3, correction simulation mineral dynamics meter
Calculate parameter and specific surface area.
5) extend the simulated time in each diagenesis period, carry out the diagenesis evolution process mould of the actual rock-forming time scale of reservoir
Intend, totally 280 million years.Analog result and measured data are contrasted, calibration model parameter, establishes porosity and drilled with the time
The quantitative simulation system of change.The differentiation curve of porosity is as shown in figure 4, porosity passes through compaction, by 30% drop at initial stage
To 15%, early diagenesis period porosity is held essentially constant, with CO2The filling of gas, mineral dissolution, porosity increase, but
After stopping injection, and temperature and pressure increase, cause a large amount of carbonate rocks to precipitate, porosity is drastically down to less than 10%, storage
Layer occurs fine and close, that is to say, that is filled with the first time of oil gas, CO2Densification of the intrusion of gas to reservoir plays pass
Key acts on, and the organic acid filling in later stage and the influence into Tibetan to porosity are little.
6) based on the material resource and diagenesis data in the range of research region, threedimensional model is built, is stored up in simulated domain scope
The diagenesis evolution process of layer, analog result are as shown in Figure 5.The result shows that analog result is consistent with region measured data, show
It is entirety based on the dissolving of feldspathoid, is converted into illite, and generates ankerite and dawsonite precipitation.But
Starting mineral and porosity, permeability it is heterogeneous to cause diagenesis to be shown as on region heterogeneous, final porosity is in
Existing heterogeneous state, there occurs different degrees of densification for reservoir.According to the distribution of the porosity of simulation spatially, Neng Gouping
Valency Favorable Reservoir.
Technical solution of the present invention is combined by serial experiment with simulation, it is proposed that complete set, science, system
The method for evaluating porosity, can establish quantitative simulation system of the porosity with Spatio-temporal Evolution, not only recover very long diagenesis evolution
The differentiation curve of process porosity, and can realize the quantitative evaluation and prediction of porosity spatially, both found out storage
The densification mechanism of layer, and can point the direction for the exploration of the resources such as oil gas, improve production efficiency.