CN106446396B - For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability - Google Patents
For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability Download PDFInfo
- Publication number
- CN106446396B CN106446396B CN201610833467.3A CN201610833467A CN106446396B CN 106446396 B CN106446396 B CN 106446396B CN 201610833467 A CN201610833467 A CN 201610833467A CN 106446396 B CN106446396 B CN 106446396B
- Authority
- CN
- China
- Prior art keywords
- carbonate
- solution cavity
- strata
- crack
- 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.)
- Expired - Fee Related
Links
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
Abstract
Provided by the present invention for determining solution cavity distribution in crack to the method for Carbonate Reservoir Permeability, it is related to oilfield prospecting developing technology, this method carries out flow simulating according to the crack on Carbonate Reservoir stratum and the true fractue spacing parameter and solution cavity distribution parameter of solution cavity, it does not need fracture and solution cavity carries out any simplified processing, and it is not necessarily to consider the border issue between crack solution cavity and matrix using Lattice Boltzmann Method unified model, all grids are uniformly processed, so that can obtain accurate result while calculation amount greatly simplifies again.Simultaneously, it is distributed according to different cracks and solution cavity, by simulate obtain considering crack and solution cavity after entire carbonate strata improve after matrix permeability, the quality for characterizing the permeability improvement of different cracks and solution cavity distribution to the carbonate strata from the ratio of intrinsic permeability with the matrix permeability after improvement realizes the distribution of crack solution cavity to the quantum chemical method of Carbonate Reservoir Permeability.
Description
Technical field
The present invention relates to oilfield prospecting developing technologies, more particularly to one kind is for determining solution cavity distribution in crack to carbonate
The method and apparatus that reservoir permeability influences.
Background technique
The percolation law of Carbonate Reservoir is influenced by its complex geological condition, reservoir include a large amount of crack and
Solution cavity.In order to study the flowing law under Carbonate Reservoir complex geological condition, need to study different cracks and solution cavity development
Influence to reservoir permeability, convenient for the research of the percolation law of development Carbonate Reservoir.
At present in the numerical simulation of Carbonate Reservoir, usually crack and solution cavity are regarded as a kind of continuous medium, use
Double medium model or multi-dielectric model are simulated, and the truth in this and crack and solution cavity noncontinuity is not met；
Another part scholar carries out numerical simulation using discrete fractures model, only accounts for the influence in crack in this simulation, and
Fracture has carried out dimensionality reduction and has simplified processing, and cannot consider the influence of solution cavity；In short, the Numerical-Mode of current Carbonate Reservoir
The true distribution situation of crack and solution cavity cannot be considered in quasi- to carry out analog study.
Summary of the invention
The present invention provides a kind of method and apparatus being distributed for determining crack solution cavity to Carbonate Reservoir Permeability,
It is intended to reduce the distribution of simulation fracture solution cavity to the difficulty of Carbonate Reservoir Permeability, simplifies the distribution of crack solution cavity to carbonate
The numerical simulation calculation amount that reservoir permeability influences improves the distribution of crack solution cavity to the Numerical-Mode of Carbonate Reservoir Permeability
Quasi- precision.
On the one hand, the present invention provides a kind of for determining solution cavity distribution in crack to the side of Carbonate Reservoir Permeability
Method, which comprises
According to seismic data and log data, the formation parameter of carbonate strata and the crack of the carbonate strata are obtained
Distribution parameter and solution cavity distribution parameter, wherein the formation parameter of the carbonate strata includes the stratum of the carbonate strata
Area, intrinsic permeability, the fractue spacing parameter of the carbonate strata include fracture dip, fracture half-length and fracture aperture,
The solution cavity distribution parameter includes solution cavity center position coordinates and solution cavity area；
According to the fluid sample and log data of the carbonate strata, obtain the carbonate strata fluid density and
Temperature；
According to equilibrium state distribution function, the first fluid speed of the carbonate strata is determined, wherein the equilibrium-like
State distribution function isWherein R is gas constant, u, ρ and T difference
For first fluid speed, density and the temperature of the carbonate strata, eiFor grid speed, ωiFor weight coefficient；
According to formulaDetermine the second fluid speed of the carbonate strata, wherein R is gas
Constant, τ are that there are relationships between slack time, with viscosity: v=(τ -0.5) RT；
According to the first fluid speed and the second fluid speed, determine in the carbonate strata horizontal direction
Mean fluid velocity；
Using Lattice Boltzmann unified model according to the first fluid speed, the mean fluid velocity, the carbon
Pressure difference between the entrance and exit on hydrochlorate stratum determines the matrix permeability of the carbonate strata；
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to store up carbonate
The infiltrative influence value of layer.
Optionally, described according to the first fluid speed and the second fluid speed, determine the carbonate strata
Mean fluid velocity in horizontal direction, comprising:
According to formulaDetermine the mean fluid velocity in the carbonate strata horizontal directionWherein, u is the first fluid speed of the carbonate strata, and u' is the second fluid speed of the carbonate strata, and τ is
Slack time.
Optionally, described to utilize Lattice Boltzmann unified model according to the first fluid speed, the mean flow
Speed, the carbonate strata entrance and exit between pressure difference, determine the matrix permeability of the carbonate strata, comprising:
Grid dividing is carried out to the carbonate strata using orthogonal grid；
For the carbonate strata after grid division, according to the fractue spacing parameter and solution cavity of the carbonate strata
Distribution parameter determines that the type of each grid, the type of the grid include matrix, crack and solution cavity；
For Lattice Boltzmann unified model, the filtrational resistance in crack and the filtrational resistance of solution cavity are set as unlimited
It is small, the filtrational resistance R of matrix1According to formula R1=vk-1It determines；
Using Lattice Boltzmann unified model, according to formulaDetermine the matrix permeability of the carbonate strata
K, whereinFor the mean fluid velocity, u is the first fluid speed, and Δ p is the entrance of the carbonate strata and goes out
Pressure difference between mouthful, l are the zone length of the carbonate strata.
Optionally, the ratio by between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity
To the infiltrative influence value of Carbonate Reservoir, comprising:
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to store up carbonate
The infiltrative influence value of layer；
Judge that crack solution cavity is distributed to the infiltrative improvement of Carbonate Reservoir according to the size of the influence value,
In, the ratio is bigger, and solution cavity distribution in crack is better to the infiltrative improvement of Carbonate Reservoir.
On the other hand, the embodiment of the present invention provides a kind of for determining solution cavity distribution in crack to Carbonate Reservoir permeability shadow
Loud device, described device include:
First obtains module, for obtaining formation parameter and the institute of carbonate strata according to seismic data and log data
State the fractue spacing parameter and solution cavity distribution parameter of carbonate strata, wherein the formation parameter of the carbonate strata includes institute
Layer spreads, the intrinsic permeability of carbonate strata are stated, the fractue spacing parameter of the carbonate strata includes fracture dip, splits
Seam half is long and fracture aperture, the solution cavity distribution parameter include solution cavity center position coordinates and solution cavity area；
Second acquisition module obtains the carbonic acid for the fluid sample and log data according to the carbonate strata
The density and temperature of salt formation fluid；
First determining module, for determining the first fluid speed of the carbonate strata according to equilibrium state distribution function
Degree, wherein the equilibrium state distribution function isWherein R is gas
Body constant, u, ρ and T are respectively first fluid speed, density and the temperature of the carbonate strata, eiFor grid speed, ωiFor
Weight coefficient；
Second determining module, for according to formulaDetermine the second fluid of the carbonate strata
Speed, wherein R is gas constant, and τ is that there are relationships between slack time, with viscosity: v=(τ -0.5) RT；
Third determining module, for determining the carbonic acid according to the first fluid speed and the second fluid speed
Mean fluid velocity in the horizontal direction of salt stratum；
4th determining module, for using Lattice Boltzmann unified model according to the first fluid speed, described flat
Pressure difference between the entrance and exit of fluid velocity, the carbonate strata, determines the matrix permeability of the carbonate strata；
5th determining module, it is molten for the ratio between the matrix permeability and the intrinsic permeability to be determined as crack
Hole is distributed to the infiltrative influence value of Carbonate Reservoir.
Optionally, the third determining module is specifically used for:
According to formulaDetermine the mean fluid velocity in the carbonate strata horizontal directionWherein, u is the first fluid speed of the carbonate strata, and u' is the second fluid speed of the carbonate strata, and τ is
Slack time.
Optionally, the 4th determining module is specifically used for:
Grid dividing is carried out to the carbonate strata using orthogonal grid；
For the carbonate strata after grid division, according to the fractue spacing parameter and solution cavity of the carbonate strata
Distribution parameter determines that the type of each grid, the type of the grid include matrix, crack and solution cavity；
For Lattice Boltzmann unified model, the filtrational resistance in crack and the filtrational resistance of solution cavity are set as unlimited
It is small, the filtrational resistance R of matrix1According to formula R1=vk-1It determines；
Using Lattice Boltzmann unified model, according to formulaDetermine the matrix infiltration of the carbonate strata
Rate k, whereinFor the mean fluid velocity, u is the first fluid speed, Δ p be the carbonate strata entrance and
Pressure difference between outlet, l are the zone length of the carbonate strata.
Optionally, the 5th determining module is specifically used for:
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to store up carbonate
The infiltrative influence value of layer；
Judge that crack solution cavity is distributed to the infiltrative improvement of Carbonate Reservoir according to the size of the influence value,
In, the ratio is bigger, and solution cavity distribution in crack is better to the infiltrative improvement of Carbonate Reservoir.
The method provided in an embodiment of the present invention for being used to determine the distribution of crack solution cavity to Carbonate Reservoir Permeability, root
Flow simulating is carried out according to the crack on Carbonate Reservoir stratum and the true fractue spacing parameter and solution cavity distribution parameter of solution cavity, no
Fracture and solution cavity is needed to carry out any simplified processing, and using Lattice Boltzmann Method unified model without considering crack
Border issue between solution cavity and matrix, all grids are uniformly processed, so that can obtain standard while calculation amount greatly simplifies again
True result.Meanwhile an intrinsic permeability is given for the carbonate strata, then according to different cracks and solution cavity point
Cloth, by simulate obtain considering crack and solution cavity after entire carbonate strata improve after matrix permeability, with the base after improvement
Matter permeability, which characterizes different cracks and solution cavity distribution from the ratio of intrinsic permeability, improves the permeability of the carbonate strata
Quality, realize the distribution of crack solution cavity to the quantum chemical method of Carbonate Reservoir Permeability.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention without any creative labor, may be used also for those of ordinary skill in the art
To obtain other drawings based on these drawings.
Fig. 1 is provided in an embodiment of the present invention for determining solution cavity distribution in crack to the side of Carbonate Reservoir Permeability
The flow diagram of method；
Fig. 2 is provided in an embodiment of the present invention for determining solution cavity distribution in crack to the dress of Carbonate Reservoir Permeability
The structural schematic diagram set.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Term " includes " in description and claims of this specification and above-mentioned attached drawing and " having " and they
Any deformation, it is intended that cover it is non-exclusive include, for example, containing the process, method of a series of steps or units, being
System, product or equipment those of are not necessarily limited to be clearly listed step or unit, but may include be not clearly listed or
For the intrinsic other step or units of these process, methods, product or equipment.
It is provided in an embodiment of the present invention to be used to determine that the distribution of crack solution cavity to the method for Carbonate Reservoir Permeability
To be executed by terminal device.Terminal device in the embodiment of the present invention can be handheld device including display screen, mobile unit,
Wearable device calculates equipment and various forms of user equipmenies (User Equipment；Referred to as: UE), mobile station
(Mobile Station；Referred to as: MS) and terminal (terminal) etc..Exemplary, the terminal device of the embodiment of the present invention can be with
It is desktop computer, server etc..
It is provided in an embodiment of the present invention to be used to determine that solution cavity distribution in crack to the method for Carbonate Reservoir Permeability, to be used
During oil-gas exploration, evaluation crack solution cavity distribution is on the infiltrative influence of Carbonate Reservoir, according to Carbonate Reservoir
The crack of layer and the true fractue spacing parameter and solution cavity distribution parameter of solution cavity carry out flow simulating, do not need fracture and solution cavity
Any simplified processing is carried out, and using Lattice Boltzmann Method unified model without considering between crack solution cavity and matrix
Border issue, all grids are uniformly processed, so that can obtain accurate result while calculation amount greatly simplifies again.Meanwhile it is right
An intrinsic permeability is given in the carbonate strata, is then distributed according to different cracks and solution cavity, is examined by simulation
Matrix permeability after entire carbonate strata improves behind worry crack and solution cavity is permeated with the matrix permeability after improvement with intrinsic
The ratio of rate is distributed the quality improved to the permeability of the carbonate strata to characterize different cracks and solution cavity, and it is molten to realize crack
Hole is distributed the quantum chemical method to Carbonate Reservoir Permeability.
Technical solution of the present invention is described in detail with specific embodiment below.These specific implementations below
Example can be combined with each other, and the same or similar concept or process may be repeated no more in some embodiments.
Fig. 1 is provided in an embodiment of the present invention for determining solution cavity distribution in crack to the side of Carbonate Reservoir Permeability
The flow diagram of method.As shown in Figure 1, the offer of the embodiment of the present invention one is used to determine that crack solution cavity to be distributed to Carbonate Reservoir
The method of Permeability the following steps are included:
Step S101: according to seismic data and log data, the formation parameter and the carbonate of carbonate strata are obtained
The fractue spacing parameter and solution cavity distribution parameter on stratum.
During the realization of step S101, the integrally texture in the oil field is usually obtained using method of seismic exploration first
Make, specifically, obtain the geological structure data in the oil field by method of seismic exploration, so by computer carry out data processing and
Draw the geologic structure diagram in the oil field.It should be noted that the specific implementation process of method of seismic exploration, those skilled in the art can
With reference to the prior art.
Certainly, during the realization of step S101, the stratum of carbonate strata can also be obtained by way of well logging
The fractue spacing parameter and solution cavity distribution parameter of parameter and the carbonate strata.But also can simultaneously using method of seismic exploration and
Well logging method obtains the formation parameter of the carbonate strata and the fractue spacing parameter of the carbonate strata and solution cavity distribution parameter, this
Inventive embodiments are not specifically limited in this embodiment.
It should be noted that the formation parameter of the carbonate strata includes the layer spreads of carbonate strata, intrinsic infiltration
Rate, the fractue spacing parameter of the carbonate strata include fracture dip, fracture half-length and fracture aperture, the carbonate strata it is molten
Hole distribution parameter includes solution cavity center position coordinates and solution cavity area.
Step S102: according to the fluid sample and log data of the carbonate strata, the carbonate strata stream is obtained
The density and temperature of body.
Specifically, can obtain the fluid sample of the carbonate strata by way of drilling extracting core, in turn, use is close
Degree meter measures the density of the carbonate strata；In addition in drilling process, according to fluid in the available stratum of temperature sensor
Temperature.
Step S103: according to equilibrium state distribution function, the first fluid speed of the carbonate strata is determined.
Specifically, according to the first fluid speed, density of carbonate strata in equilibrium state distribution function and temperature three
Between relationship, determine the first fluid speed of the carbonate strata.
Further, which isIts
Middle R is gas constant, and u, ρ and T are respectively first fluid speed, density and the temperature of the carbonate strata, eiFor grid speed,
ωiFor weight coefficient.
It should be noted that for weight coefficient ωi, the embodiment of the present invention is not specifically limited, those skilled in the art
It can refer to prior art setting.It is exemplary, for the Lattice Boltzmann unified model of two-dimensional D2Q9, weight coefficient ωi
For ω0=4/9, work as i=1, when 2,3,4, ωi=1/9, work as i=5, when 6,7,8, ωi=1/36.
It should be noted that for grid speed eiSpecific value, the embodiment of the present invention is not specifically limited, this field
Technical staff can refer to the prior art.It is exemplary, for the Lattice Boltzmann unified model of two-dimensional D2Q9, grid speed
eiFor e0=(0,0), e1,3=(± 1,0), e2,4=(0, ± 1), e5,7=(± 1, ± 1) and
Step S104: according to formulaDetermine the second fluid speed of the carbonate strata.
Specifically, after obtaining the first fluid speed of the carbonate strata, according to first fluid speed and second
Relationship between body speed calculates the second fluid speed of the carbonate strata.Further, according to formulaDetermine the second fluid speed of the carbonate strata；Wherein, R is gas constant, is one in states of matter
The physical constant of each thermodynamic function is linked in equation, in general, R=8.3138462；τ is slack time, for refer to object by
Power deformation, the time needed for material restores normal condition after external force releases, there are relationships between slack time τ and viscosity: v=
(τ-0.5)RT。
Step S105: according to the first fluid speed and the second fluid speed, the carbonate strata water is determined
Square upward mean fluid velocity.
Specifically, according to first fluid speed and second fluid speed, the carbonate strata is calculated in the horizontal direction
Mean fluid velocity.Further, according to formulaDetermine the carbonate strata in the horizontal direction
Mean fluid velocityWherein, u is the first fluid speed of the carbonate strata, and u' is the second of the carbonate strata
Body speed, τ are slack time.
Step S106: using Lattice Boltzmann unified model according to the first fluid speed, mean flow speed
It spends, the pressure difference between the entrance and exit of the carbonate strata, determines the matrix permeability of the carbonate strata.
Specifically, grid dividing is carried out using orthogonal grid first against the carbonate strata, after grid division
The carbonate strata determines the kind of each grid according to the fractue spacing parameter of the carbonate strata and solution cavity distribution parameter
Class, the type of the grid include matrix, crack and solution cavity, exemplary, according to the fractue spacing parameter of the carbonate strata and molten
Hole distribution parameter determines that each grid is matrix or crack either solution cavity；After the type that grid has been determined, for lattice
Sub- Boltzmann unified model sets infinitely small for the filtrational resistance of the filtrational resistance in its crack and solution cavity, the seepage flow of matrix
Resistance R1According to formula R1=vk-1It determines；Grid hereby graceful unified model is finally utilized, according to formulaDetermine the carbonic acid
The matrix permeability k on salt stratum, whereinFor mean fluid velocity, u is first fluid speed, and Δ p is the carbonate strata
Pressure difference between entrance and exit, l are the zone length of the carbonate strata.
Step S107: the ratio between the matrix permeability and the intrinsic permeability is determined as the distribution pair of crack solution cavity
The infiltrative influence value of Carbonate Reservoir.
Specifically, the ratio between the matrix permeability and intrinsic permeability of determining above-mentioned acquisition, and then the ratio is determined
It is distributed for crack solution cavity to the infiltrative influence value of Carbonate Reservoir, crack solution cavity point is then judged according to the size of the influence value
Cloth is to the infiltrative improvement of Carbonate Reservoir, wherein the ratio is bigger, and crack solution cavity is distributed to Carbonate Reservoir permeability
Improvement it is better.
The method provided in an embodiment of the present invention for being used to determine the distribution of crack solution cavity to Carbonate Reservoir Permeability, root
Flow simulating is carried out according to the crack on Carbonate Reservoir stratum and the true fractue spacing parameter and solution cavity distribution parameter of solution cavity, no
Fracture and solution cavity is needed to carry out any simplified processing, and using Lattice Boltzmann Method unified model without considering crack
Border issue between solution cavity and matrix, all grids are uniformly processed, so that can obtain standard while calculation amount greatly simplifies again
True result.Meanwhile an intrinsic permeability is given for the carbonate strata, then according to different cracks and solution cavity point
Cloth, by simulate obtain considering crack and solution cavity after entire carbonate strata improve after matrix permeability, with the base after improvement
Matter permeability, which characterizes different cracks and solution cavity distribution from the ratio of intrinsic permeability, improves the permeability of the carbonate strata
Quality, realize the distribution of crack solution cavity to the quantum chemical method of Carbonate Reservoir Permeability.
Fig. 2 is provided in an embodiment of the present invention a kind of for determining solution cavity distribution in crack to Carbonate Reservoir Permeability
Device structural schematic diagram.As shown in Fig. 2, provided in an embodiment of the present invention for determining that carbonate is stored up in solution cavity distribution in crack
Layer Permeability device include:
First obtains module 210, for according to seismic data and log data, obtain carbonate strata formation parameter and
The fractue spacing parameter and solution cavity distribution parameter of the carbonate strata, wherein the formation parameter of the carbonate strata includes
Layer spreads, the intrinsic permeability of the carbonate strata, the fractue spacing parameter of the carbonate strata include fracture dip,
Fracture half-length and fracture aperture, the solution cavity distribution parameter include solution cavity center position coordinates and solution cavity area；
Second acquisition module 220 obtains the carbon for the fluid sample and log data according to the carbonate strata
The density and temperature of hydrochlorate formation fluid；
First determining module 230, for determining the first fluid of the carbonate strata according to equilibrium state distribution function
Speed, wherein the equilibrium state distribution function isWherein R is
Gas constant, u, ρ and T are respectively first fluid speed, density and the temperature of the carbonate strata, eiFor grid speed, ωi
For weight coefficient；
Second determining module 240, for according to formulaDetermine the second of the carbonate strata
Body speed, wherein R is gas constant, and τ is that there are relationships between slack time, with viscosity: v=(τ -0.5) RT；
Third determining module 250, for determining the carbon according to the first fluid speed and the second fluid speed
The upward mean fluid velocity of hydrochlorate water flooding square；
4th determining module 260, for using Lattice Boltzmann unified model according to the first fluid speed, described
Mean fluid velocity, the carbonate strata entrance and exit between pressure difference, determine the carbonate strata matrix infiltration
Rate；
5th determining module 270, for the ratio between the matrix permeability and the intrinsic permeability to be determined as splitting
Solution cavity distribution is stitched to the infiltrative influence value of Carbonate Reservoir.
Optionally, third determining module 250 is specifically used for:
According to formulaDetermine the mean fluid velocity in the carbonate strata horizontal directionWherein, u is the first fluid speed of the carbonate strata, and u' is the second fluid speed of the carbonate strata, and τ is pine
It relaxes the time.
Optionally, the 4th determining module 260 is specifically used for:
Grid dividing is carried out to the carbonate strata using orthogonal grid；
For the carbonate strata after grid division, according to the fractue spacing parameter and solution cavity of the carbonate strata
Distribution parameter determines that the type of each grid, the type of the grid include matrix, crack and solution cavity；
For Lattice Boltzmann unified model, the filtrational resistance in crack and the filtrational resistance of solution cavity are set as unlimited
It is small, the filtrational resistance R of matrix1According to formula R1=vk-1It determines；
Using Lattice Boltzmann unified model, according to formulaDetermine the matrix permeability of the carbonate strata
K, whereinFor the mean fluid velocity, u is the first fluid speed, and Δ p is the entrance of the carbonate strata and goes out
Pressure difference between mouthful, l are the zone length of the carbonate strata.
Optionally, the 5th determining module 270 is specifically used for:
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to store up carbonate
The infiltrative influence value of layer；
Judge that crack solution cavity is distributed to the infiltrative improvement of Carbonate Reservoir according to the size of the influence value,
In, the ratio is bigger, and solution cavity distribution in crack is better to the infiltrative improvement of Carbonate Reservoir.
It should be understood that provided by the above embodiment for determining solution cavity distribution in crack to Carbonate Reservoir permeability shadow
Loud device, only the example of the division of the above functional modules, in practical application, can according to need and will be above-mentioned
Function distribution is completed by different functional modules, i.e., the internal structure of equipment is divided into different functional modules, with complete with
The all or part of function of upper description.In addition, provided by the above embodiment for determining that carbonate is stored up in solution cavity distribution in crack
The device of layer Permeability with for determining solution cavity distribution in crack to the embodiment of the method category of Carbonate Reservoir Permeability
In same design, specific implementation process is detailed in embodiment of the method, and which is not described herein again.
The device provided in an embodiment of the present invention for being used to determine the distribution of crack solution cavity to Carbonate Reservoir Permeability, root
Flow simulating is carried out according to the crack on Carbonate Reservoir stratum and the true fractue spacing parameter and solution cavity distribution parameter of solution cavity, no
Fracture and solution cavity is needed to carry out any simplified processing, and using Lattice Boltzmann Method unified model without considering crack
Border issue between solution cavity and matrix, all grids are uniformly processed, so that can obtain standard while calculation amount greatly simplifies again
True result.Meanwhile an intrinsic permeability is given for the carbonate strata, then according to different cracks and solution cavity point
Cloth, by simulate obtain considering crack and solution cavity after entire carbonate strata improve after matrix permeability, with the base after improvement
Matter permeability, which characterizes different cracks and solution cavity distribution from the ratio of intrinsic permeability, improves the permeability of the carbonate strata
Quality, realize the distribution of crack solution cavity to the quantum chemical method of Carbonate Reservoir Permeability.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above-mentioned each method embodiment can lead to
The relevant hardware of program instruction is crossed to complete.Program above-mentioned can be stored in a computer readable storage medium.The journey
When being executed, execution includes the steps that above-mentioned each method embodiment to sequence；And storage medium above-mentioned include: ROM, RAM, magnetic disk or
The various media that can store program code such as person's CD.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement；And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (8)
1. a kind of for determining solution cavity distribution in crack to the method for Carbonate Reservoir Permeability, which is characterized in that the side
Method includes:
According to seismic data and log data, the formation parameter of carbonate strata and the fractue spacing of the carbonate strata are obtained
Parameter and solution cavity distribution parameter, wherein the formation parameter of the carbonate strata include the carbonate strata layer spreads,
Intrinsic permeability, the fractue spacing parameter of the carbonate strata include fracture dip, fracture half-length and fracture aperture, described molten
Hole distribution parameter includes solution cavity center position coordinates and solution cavity area；
According to the fluid sample and log data of the carbonate strata, the density and temperature of the carbonate strata fluid are obtained
Degree；
According to equilibrium state distribution function, the first fluid speed of the carbonate strata is determined, wherein the equilibrium state point
Cloth function isWherein R is gas constant, and u, ρ and T are respectively institute
State first fluid speed, density and the temperature of carbonate strata, eiFor grid speed, ωiFor weight coefficient；
According to formulaDetermining the second fluid speed of the carbonate strata, wherein R is gas constant,
τ is that there are relationships between slack time, with viscosity: v=(τ -0.5) RT, v is viscosity coefficient；
According to the first fluid speed and the second fluid speed, being averaged in the carbonate strata horizontal direction is determined
Fluid velocity；
Using Lattice Boltzmann unified model according to the first fluid speed, the mean fluid velocity, the carbonate
Pressure difference between the entrance and exit on stratum determines the matrix permeability of the carbonate strata；
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to seep Carbonate Reservoir
The influence value of permeability.
2. the method according to claim 1, wherein described according to the first fluid speed and the second
Body speed determines the mean fluid velocity in the carbonate strata horizontal direction, comprising:
According to formulaDetermine the mean fluid velocity in the carbonate strata horizontal directionIts
In, u is the first fluid speed of the carbonate strata, and u' is the second fluid speed of the carbonate strata, when τ is relaxation
Between.
3. the method according to claim 1, wherein described utilize Lattice Boltzmann unified model according to
First fluid speed, the mean fluid velocity, the carbonate strata entrance and exit between pressure difference, determine the carbonic acid
The matrix permeability on salt stratum, comprising:
Grid dividing is carried out to the carbonate strata using orthogonal grid；
For the carbonate strata after grid division, it is distributed according to the fractue spacing parameter and solution cavity of the carbonate strata
Parameter determines that the type of each grid, the type of the grid include matrix, crack and solution cavity；
For Lattice Boltzmann unified model, the filtrational resistance in crack and the filtrational resistance of solution cavity are set as infinitely small, base
The filtrational resistance R of matter1According to formula R1=vk-1It determines；
Using Lattice Boltzmann unified model, according to formulaDetermine the matrix permeability k of the carbonate strata,
In,For the mean fluid velocity, u is the first fluid speed, and △ p is between the entrance and exit of the carbonate strata
Pressure difference, l be the carbonate strata zone length.
4. the method according to claim 1, wherein described by the matrix permeability and the intrinsic permeability
Between ratio be determined as crack solution cavity be distributed to the infiltrative influence value of Carbonate Reservoir, comprising:
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to seep Carbonate Reservoir
The influence value of permeability；
Judge solution cavity distribution in crack to the infiltrative improvement of Carbonate Reservoir according to the size of the influence value, wherein institute
It is bigger to state ratio, solution cavity distribution in crack is better to the infiltrative improvement of Carbonate Reservoir.
5. a kind of for determining solution cavity distribution in crack to the device of Carbonate Reservoir Permeability, which is characterized in that the dress
It sets and includes:
First obtains module, for obtaining the formation parameter and the carbon of carbonate strata according to seismic data and log data
The fractue spacing parameter and solution cavity distribution parameter on hydrochlorate stratum, wherein the formation parameter of the carbonate strata includes the carbon
Layer spreads, the intrinsic permeability on hydrochlorate stratum, the fractue spacing parameter of the carbonate strata include fracture dip, crack half
Long and fracture aperture, the solution cavity distribution parameter include solution cavity center position coordinates and solution cavity area；
Second acquisition module obtains the carbonic acid saline land for the fluid sample and log data according to the carbonate strata
The density and temperature of layer fluid；
First determining module, for determining the first fluid speed of the carbonate strata according to equilibrium state distribution function,
In, the equilibrium state distribution function isWherein R is that gas is normal
Number, u, ρ and T are respectively first fluid speed, density and the temperature of the carbonate strata, eiFor grid speed, ωiFor weight
Coefficient；
Second determining module, for according to formulaDetermine the second fluid speed of the carbonate strata,
Wherein, R is gas constant, and τ is that there are relationships between slack time, with viscosity: v=(τ -0.5) RT, v is viscosity coefficient；
Third determining module, for determining the carbonic acid saline land according to the first fluid speed and the second fluid speed
Mean fluid velocity in layer horizontal direction；
4th determining module, for utilizing Lattice Boltzmann unified model according to the first fluid speed, the mean flow
Body speed, the carbonate strata entrance and exit between pressure difference, determine the matrix permeability of the carbonate strata；
5th determining module, for the ratio between the matrix permeability and the intrinsic permeability to be determined as crack solution cavity point
Cloth is to the infiltrative influence value of Carbonate Reservoir.
6. device according to claim 5, which is characterized in that the third determining module is specifically used for:
According to formulaDetermine the mean fluid velocity in the carbonate strata horizontal directionIts
In, u is the first fluid speed of the carbonate strata, and u' is the second fluid speed of the carbonate strata, when τ is relaxation
Between.
7. device according to claim 5, which is characterized in that the 4th determining module is specifically used for:
Grid dividing is carried out to the carbonate strata using orthogonal grid；
For the carbonate strata after grid division, it is distributed according to the fractue spacing parameter and solution cavity of the carbonate strata
Parameter determines that the type of each grid, the type of the grid include matrix, crack and solution cavity；
For grid, hereby graceful unified model, the filtrational resistance in crack and the filtrational resistance of solution cavity are set as infinitely small, matrix
Filtrational resistance R1According to formula R1=vk-1It determines；
Using Lattice Boltzmann unified model, according to formulaDetermine the matrix permeability k of the carbonate strata,
In,For the mean fluid velocity, u is the first fluid speed, and △ p is between the entrance and exit of the carbonate strata
Pressure difference, l be the carbonate strata zone length.
8. device according to claim 5, which is characterized in that the 5th determining module is specifically used for:
Ratio between the matrix permeability and the intrinsic permeability is determined as the distribution of crack solution cavity to seep Carbonate Reservoir
The influence value of permeability；
Judge solution cavity distribution in crack to the infiltrative improvement of Carbonate Reservoir according to the size of the influence value, wherein institute
It is bigger to state ratio, solution cavity distribution in crack is better to the infiltrative improvement of Carbonate Reservoir.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610833467.3A CN106446396B (en) | 2016-09-20 | 2016-09-20 | For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610833467.3A CN106446396B (en) | 2016-09-20 | 2016-09-20 | For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106446396A CN106446396A (en) | 2017-02-22 |
| CN106446396B true CN106446396B (en) | 2019-04-16 |
Family
ID=58166334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610833467.3A Expired - Fee Related CN106446396B (en) | 2016-09-20 | 2016-09-20 | For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106446396B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107169192B (en) * | 2017-05-10 | 2020-06-09 | 中国石油天然气股份有限公司 | Method and device for identifying size of karst cave while drilling |
| CN110529105B (en) * | 2018-05-23 | 2022-06-03 | 中国石油天然气股份有限公司 | Design method and design device for development of multi-medium water-reservoir horizontal well |
| CN112329353B (en) * | 2020-09-25 | 2022-11-01 | 中国石油天然气股份有限公司 | Method and device for determining channeling flow, computer equipment and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009068550A1 (en) * | 2007-11-27 | 2009-06-04 | Basf Se | Self-incompatible solvent |
| CN104948176A (en) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | Method for identifying carbonate reservoir fractures based on permeability increasing rate |
| CN105350960A (en) * | 2015-12-07 | 2016-02-24 | 西南石油大学 | Method of determining fractured horizontal well crack parameters of low-permeability anisotropic gas reservoir |
| CN105793520A (en) * | 2013-10-17 | 2016-07-20 | 西格纳化学有限责任公司 | Alkali metal silicides for high temperature wettability improvement in carbonate reservoirs |
| CN105891083A (en) * | 2016-04-05 | 2016-08-24 | 中国石油大学(华东) | Digital core two-phase flowing simulation result quantitative representation method based on erosion-expansion algorithm and application of method |
-
2016
- 2016-09-20 CN CN201610833467.3A patent/CN106446396B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009068550A1 (en) * | 2007-11-27 | 2009-06-04 | Basf Se | Self-incompatible solvent |
| CN105793520A (en) * | 2013-10-17 | 2016-07-20 | 西格纳化学有限责任公司 | Alkali metal silicides for high temperature wettability improvement in carbonate reservoirs |
| CN104948176A (en) * | 2015-05-08 | 2015-09-30 | 西南石油大学 | Method for identifying carbonate reservoir fractures based on permeability increasing rate |
| CN105350960A (en) * | 2015-12-07 | 2016-02-24 | 西南石油大学 | Method of determining fractured horizontal well crack parameters of low-permeability anisotropic gas reservoir |
| CN105891083A (en) * | 2016-04-05 | 2016-08-24 | 中国石油大学(华东) | Digital core two-phase flowing simulation result quantitative representation method based on erosion-expansion algorithm and application of method |
Non-Patent Citations (2)
| Title |
|---|
| 利用格子Boltzmann方法计算页岩渗透率;张磊等;《中国石油大学学报（自然科学版）》;20140220;第38卷(第1期);全文 |
| 基于格子玻尔兹曼方法的碳酸盐岩数字岩心渗流特征分析;王晨晨等;《中国石油大学学报（自然科学版）》;20121220;第36卷(第6期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106446396A (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wen et al. | Use of border regions for improved permeability upscaling | |
| Kaufmann et al. | Karst aquifer evolution in fractured rocks | |
| AU2018201250A1 (en) | Determining geomechanics completion quality | |
| CN105510993B (en) | The buried extrusion pressing type complexity gypseous-salt rock identification of foreland basin and distribution forecasting method | |
| CN106951660A (en) | A kind of marine clastics horizontal well reservoir log interpretation method and device | |
| CN108603402A (en) | The variation of capillary pressure and relative permeability in porous media caused by mineral deposits and dissolving is modeled and predicted | |
| Guihéneuf et al. | Insights about transport mechanisms and fracture flow channeling from multi-scale observations of tracer dispersion in shallow fractured crystalline rock | |
| CN106446396B (en) | For determining that crack solution cavity is distributed the method and apparatus influenced on reservoir permeability | |
| NO339162B1 (en) | Fracture network characterization method | |
| CN104914470B (en) | Carbonate rock fractured cave reservoir reserves bearing calibration | |
| Griggs et al. | Ground‐water flow dynamics and development strategies at the atoll scale | |
| Santos et al. | Characterization of natural fracture systems: Analysis of uncertainty effects in linear scanline results | |
| CN105019892B (en) | A method of simulation fracture hole type Reservoir Body electric logging response | |
| CN107366534A (en) | It is roughened the determination method and apparatus of permeability | |
| CN104834007B (en) | The method that carbonate rock fractured cave type reservoir filling operation is calculated during seismic inversion | |
| CN110259439A (en) | For evaluating the method and device of carbonate rock fault sealing property Evolutionary History | |
| CN110094196A (en) | A kind of carbonate rock open-hole horizontal well segmentation acid fracturing effect evaluation method | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| Knochenmus et al. | Descriptions of anisotropy and heterogeneity and their effect on ground-water flow and areas of contribution to public supply wells in a karst carbonate aquifer system | |
| Fabbri et al. | Estimation of hydraulic conductivity using the slug test method in a shallow aquifer in the Venetian Plain (NE, Italy) | |
| CN107609284A (en) | A kind of method of quantitative assessment Landslide hazards | |
| Bierkens et al. | Upscaling hydraulic conductivity: theory and examples from geohydrological studies | |
| Hambrey et al. | Structural evolution during cyclic glacier surges: 1. Structural glaciology of Trapridge Glacier, Yukon, Canada | |
| CN103675945B (en) | A kind of method and apparatus of the saturation degree for determining hole type reservoir | |
| CN106644839B (en) | Rock mass transporting capability determination method for parameter and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190416 Termination date: 20190920 |