CN106546525A - The method and apparatus for setting up three-dimensional penetration rate model - Google Patents
The method and apparatus for setting up three-dimensional penetration rate model Download PDFInfo
- Publication number
- CN106546525A CN106546525A CN201610944944.3A CN201610944944A CN106546525A CN 106546525 A CN106546525 A CN 106546525A CN 201610944944 A CN201610944944 A CN 201610944944A CN 106546525 A CN106546525 A CN 106546525A
- Authority
- CN
- China
- Prior art keywords
- target reservoir
- fracture
- density region
- well
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
Abstract
The embodiment of the present application provides a kind of method and apparatus for setting up three-dimensional penetration rate model, wherein, the method is comprised the following steps：Obtain the result data of target reservoir；According to the result data of target reservoir, the stress field simulation geological model of target reservoir is set up；According to the result data and target reservoir Image Logging Data of target reservoir, by stress field simulation geological model, target reservoir Stress Field Distribution data are calculated；According to target reservoir Stress Field Distribution data, target reservoir Image Logging Data and target reservoir well test data, solution obtains correcting permeability；Three-dimensional penetration rate model is set up according to correction permeability.As the program considers the impact of microcrack, and three-dimensional penetration rate model is set up according to correction permeability, therefore solve that three-dimensional penetration rate model precision present in existing modeling method is low, the technical problem that error is big, realize the technique effect of the three-dimensional penetration rate model precision for improving set up.
Description
Technical field
The application is related to oil reservoir development technical field, more particularly to a kind of method and dress for setting up three-dimensional penetration rate model
Put.
Background technology
In oil reservoir development field, when researching and analysing to target reservoir, often well logging first according to target reservoir
Data and experimental result data, set up the three-dimensional penetration rate model with regard to target reservoir, then by the three-dimensional penetration rate model pair
Target reservoir launches specific reservoir engineering research and reservoir numerical simulation research.
At present, the main modeling method for adopting is generally that, according to well logging permeability, foundation is oozed with regard to the three-dimensional of target reservoir
Saturating rate model.However, in the specific implementation, in the geologic reservoir of target reservoir, often development has many microcracks, and this is slightly
Crack can affect the migration of the fluid in porous media, and then cause in actual development, the true permeability meeting of geologic reservoir
Than logging well, permeability is much higher.Therefore, cause by existing modeling method, the three-dimensional directly set up according to well logging permeability is oozed
Thoroughly rate model often has that precision is low, the technical problem that error is big.When studying to ultra-permeable reservior, existing modeling side
Above-mentioned technical problem produced by method is particularly evident.
For the problems referred to above, effective solution is not yet proposed at present.
The content of the invention
The embodiment of the present application provides a kind of method and apparatus for setting up three-dimensional penetration rate model, to solve existing modeling side
The technical problem that model accuracy is low, error is big that method is present.
The embodiment of the present application provides a kind of method for setting up three-dimensional penetration rate model, including：
Obtain the result data of target reservoir；
According to the result data of the target reservoir, the stress field simulation geological model of target reservoir is set up；
According to the well-log information of the result data and target reservoir of the target reservoir, by the stress field simulation geology
Model, is calculated target reservoir Stress Field Distribution data；
According to the well testing of the target reservoir Stress Field Distribution data, the well-log information of the target reservoir and target reservoir
Data, solution obtain correcting permeability；
According to the correction permeability, the three-dimensional penetration rate model of the target reservoir is set up.
In a preferred embodiment, the result data includes：Rock mechanics relevant parameter, stress field size and should
The horizontal biggest principal stress direction in the field of force.
In a preferred embodiment, the result data for obtaining target reservoir, including：
Obtain the drill cores sample of target reservoir；
Rock triaxial stress test is carried out to the drill cores sample of the target reservoir, the rock mechanics is obtained related
Parameter；
Acoustic Emission of Rock experiment is carried out to the drill cores sample of the target reservoir, the stress field size is obtained；
According to the well-log information of the target reservoir, the level of stress field biggest principal stress direction is obtained.
In a preferred embodiment, the result data according to the target reservoir, sets up the stress field of target reservoir
Simulation geological model, including：
According to the result data of the target reservoir, divided by map grids, set up the stress field of the target reservoir
Simulation geological model.
In a preferred embodiment, according to the target reservoir Stress Field Distribution data, the survey of the target reservoir
The well test data of well data and target reservoir, solution obtain correcting permeability, including：
According to the target reservoir Stress Field Distribution data, the gash fracture rate of each measuring point in target reservoir is calculated
Distributed constant value；
According to the gash fracture rate of each measuring point in the target reservoir Stress Field Distribution data and the target reservoir point
Cloth parameter value, is calculated the fracture spacing predictive value of each measuring point in target reservoir；
According to the fracture spacing predictive value of each measuring point in the target reservoir, the target reservoir is divided into into multiple splitting
Seam development density region；
According to the well-log information and well test data of the target reservoir, solution obtains each measuring point in the target reservoir
Correction permeability.
In a preferred embodiment, it is described according to the target reservoir Stress Field Distribution data, it is calculated target
The gash fracture rate distributed constant value of each measuring point in oil reservoir, including：
According to the target reservoir Stress Field Distribution data, according to below equation, it is calculated in the target reservoir each
The gash fracture rate distributed constant value of individual measuring point：
In(j)=τn(j)/[τn](j)
It(j)=σt(j)/[σt](j)
I(j)=α(j)It(j)+β(j)In(j)
Wherein, [τn](j)To number the shearing strength of the rock of the measuring point for being j, C(j)To number the direct stress of the measuring point for being j
The shearing strength of rock, σ when being zeron(j)To number the plane of fracture shear stress of the measuring point for being j,It is j's for numbering
The coefficient of internal friction of measuring point, In(j)To number the crack shear fracture rate of the measuring point for being j, τn(j)To number the shear stress of the measuring point for being j
Value, It(j)To number the crack tension fracture rate of the measuring point for being j, σt(j)To number the plane of fracture tensile stress size of the measuring point for being j,
[σt](j)To number the tensile strength of rock of the measuring point for being j, α(j)To number the tension fissure ratio of the measuring point for being j, β(j)For numbering
For the shear crack ratio of the measuring point of j, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j.
In a preferred embodiment, it is described according to the target reservoir Stress Field Distribution data and the target reservoir
In each measuring point gash fracture rate distributed constant value, be calculated the fracture spacing predictive value of each measuring point in target reservoir,
Including：
According to the gash fracture rate of each measuring point in the target reservoir Stress Field Distribution data and the target reservoir point
Cloth parameter value, according to below equation, is calculated the fracture spacing predictive value of each measuring point in the target reservoir：
I(j)＜ Ic,ξ(j)=a1(j)I(j) 2+a2(j)
I(j)≥Ic,ξ(j)=a1(j)I(j) 2+a2(j)+a3(j)(A(j)+a4(j))
Wherein, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j, IcFor I(j)Marginal value, ξ(j)To compile
Number for j measuring point fracture spacing predictive value, a1(j)、a2(j)、a3(j)、a4(j)To number the undetermined constant of the measuring point for being j, A(j)For
Numbering is the energy value of the measuring point of j.
In a preferred embodiment, the well-log information and well test data according to the target reservoir, solves
The correction permeability of each measuring point in the target reservoir, including：
According to the well-log information of the target reservoir, during solution obtains the plurality of fracture development density region, each crack is sent out
Educate the meansigma methodss of the well logging permeability in density region；
According to the target reservoir well test data, solution obtains each fracture development in the plurality of fracture development density region
Well testing permeability in density region；
According to the meansigma methodss of the well logging permeability in each fracture development density region in the plurality of fracture development density region
With well testing permeability value, solution obtains in the plurality of fracture development density region correction coefficient in each fracture development density region；
According to each fracture development density in the well-log information and the plurality of fracture development density region of the target reservoir
Correction coefficient in area, solution obtain the correction permeability of each measuring point in the target reservoir.
In a preferred embodiment, the well-log information according to the target reservoir, solution obtain the plurality of crack
The meansigma methodss of the well logging permeability in development density region in each fracture development density region, including：
According to the well-log information of the target reservoir, according to below equation, solution obtains the plurality of fracture development density
The meansigma methodss of the well logging permeability in area in each fracture development density region：
Wherein,For the well logging permeability meansigma methodss in the F fracture spacing development density region in target reservoir, K
(φ(j),MZ(j),Vsh(j))(j)Well logging for the measuring point that the numbering in the i-th substratum in the F fracture spacing development density region is j is oozed
Saturating rate, φ(j)、MZ(j)、Vsh(j)、VshNumbering in respectively the F fracture spacing development density region in the i-th substratum is the measuring point of j
Porosity, median grain diameter, shale content, hi(F)For the i-th small thickness in the F fracture spacing development density region, h(F)For
Oil reservoir actual production interval thickness in F fracture spacing development density region, niFor in the F fracture spacing development density region i-th
The measuring point number included in substratum, N(F)For the little number of plies included in the F fracture spacing development density region.
In a preferred embodiment, the well test data according to the target reservoir, solution obtain the plurality of crack
Well testing permeability in development density region in each fracture development density region, including：
According to the well test data of the target reservoir, according to below equation, solution obtains the plurality of fracture development density
Well testing permeability in area in each fracture development density region：
Wherein, KT(F)For the well testing permeability in the F fracture spacing development density region in target reservoir, Q(F)For target
Well test rate in oil reservoir in the F fracture spacing development density region, μ(F)For the ground in the F fracture spacing development density region
Layer viscosity of crude, re(F)For the drainage radius in the F fracture spacing development density region, rw(F)It is close for the F fracture spacing development
Oil well radius in degree area, Δ P(F)For the producing pressure differential in the F fracture spacing development density region, h(F)It is close for the F crack
Well testing reservoir thickness in degree development density region.
In a preferred embodiment, according in each fracture development density region in the plurality of fracture development density region
Well logging permeability meansigma methodss and well testing permeability value, solution obtains each crack in the plurality of fracture development density region and sends out
Correction coefficient in density region is educated, including：
According to the meansigma methodss of the well logging permeability in each fracture development density region in the plurality of fracture development density region
With well testing permeability value, according to below equation, solution obtains each fracture development density in the plurality of fracture development density region
Correction coefficient in area：
Wherein, the λ(F)For correction coefficient in the F fracture spacing development density region, KT(F)Send out for the F fracture spacing
Well testing permeability in density region is educated,For permeability meansigma methodss of logging well in the F fracture spacing development density region.
In a preferred embodiment, the well-log information and the plurality of fracture development density according to the target reservoir
Correction coefficient in each fracture development density region in area, solution obtain the correction permeability of each measuring point in the target reservoir,
Including：
According to each fracture development density in the well-log information and the plurality of fracture development density region of the target reservoir
Correction coefficient in area, according to below equation, solution obtains the correction permeability of each measuring point in the target reservoir：
Wherein, the PERM(j)In target reservoir numbering be j measuring points correction permeability, λ(F)It is that j measuring points are located for numbering
The F fracture spacing development density region in correction coefficient, K (φ(j),MZ(j),Vsh(j)) it is that the well logging for numbering the measuring point for being j is oozed
Saturating rate.
Based on identical inventive concept, the embodiment of the present application additionally provides a kind of device for setting up three-dimensional penetration rate model,
Including：
Acquisition module, for obtaining the result data of target reservoir；
First sets up module, for the result data according to the target reservoir, sets up the stress field of the target reservoir
Simulation geological model；
Computing module, for the well-log information of result data and target reservoir according to the target reservoir, by described
Stress field simulation geological model, is calculated target reservoir Stress Field Distribution data；
Correction module, for according to the target reservoir Stress Field Distribution data, the target reservoir well-log information and mesh
Mark oil reservoir well test data, solution obtain the correction permeability of each measuring point in target reservoir；
Second sets up module, for the correction permeability according to each measuring point in the target reservoir, sets up the target
The three-dimensional penetration rate model of oil reservoir.
In the embodiment of the present application, shadow of the microcrack to permeability is considered during three-dimensional penetration rate model is set up
Ring, first permeability of logging well is corrected, the three-dimensional penetration rate model of target reservoir is set up further according to correction permeability.Because examining
Impact of the microcrack to permeability has been considered, so as to the technology that precision is low, error is big is solved present in existing modeling method
Problem, has reached the purpose for improving three-dimensional penetration rate model accuracy.
Description of the drawings
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
Accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise for not paying creative labor
Under, can be with according to these other accompanying drawings of accompanying drawings acquisition.
Fig. 1 is the process chart of the method for the three-dimensional penetration rate model of foundation according to the embodiment of the present application；
Fig. 2 is the composition structure chart of the device of the three-dimensional penetration rate model of foundation according to the embodiment of the present application；
Fig. 3 is that the crack that the method, device of the three-dimensional penetration rate model of foundation provided using the embodiment of the present application is obtained is broken
Split rate distributed constant Distribution value prognostic chart；
Fig. 4 is that the crack that the method, device of the three-dimensional penetration rate model of foundation provided using the embodiment of the present application is obtained is close
Degree forecast of distribution figure.
Specific embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement
Example is only some embodiments of the present application, rather than the embodiment of whole.Based on the embodiment in the application, this area is common
The every other embodiment obtained under the premise of creative work is not made by technical staff, should all belong to the application protection
Scope.
In view of existing modeling method, the geologic reservoir of target reservoir is set up according to well logging permeability directly typically
Three-dimensional penetration rate model.And directly using the permeability of well logging, due to not accounting for the microcrack developed in reservoir to reality
The impact of permeability, causes set up three-dimensional penetration rate model inaccurate, often has that precision is low, the technology that error is big is asked
Topic.For produce above-mentioned technical problem basic reason, the application take into account target reservoir geologic reservoir in microcrack to oozing
The impact of rate thoroughly, and by being corrected permeability of logging well according to well test data, obtain more accurately correcting permeability, then
According to correction permeability, three-dimensional penetration rate model is set up.Such that it is able to solve modeling accuracy present in existing modeling method it is poor,
The big technical problem of error, realizes the technique effect of the accuracy of the three-dimensional penetration rate model for improving set up.
Based on above-mentioned thinking thinking, this application provides a kind of method for setting up three-dimensional penetration rate model.Refer to Fig. 1.
The method of the three-dimensional penetration rate model of foundation that the application is provided, may comprise steps of.
Step 11：Obtain the result data of target reservoir.
In the present embodiment, according to practical study needs, the above results data can include：The related ginseng of rock mechanics
Number, stress field size and level of stress field biggest principal stress direction.Certainly, except above-mentioned cited rock mechanics relevant parameter,
Outside stress field size and level of stress field biggest principal stress direction, according to the situation that is embodied as and actual demand, the above results number
According to other corresponding geologic datas can also be included.
In one preferred embodiment, in order to obtain result data, first target reservoir can be sampled, then leads to
Cross carries out the corresponding result data of test acquisition to sample.Specifically can implement according to following steps：
S1：Obtain the drill cores sample of target reservoir.
S2：Rock triaxial stress test is carried out to the drill cores sample of the target reservoir, the rock mechanics is obtained
Relevant parameter.
S3：Acoustic Emission of Rock experiment is carried out to the drill cores sample of the target reservoir, the stress field size is obtained.
S4：According to the well-log information of the target reservoir, the level of stress field biggest principal stress direction is obtained.
Step 12：According to the result data of the target reservoir, the stress field simulation Geological Model of the target reservoir is set up
Type.
In one preferred embodiment, in order to make full use of result data, stress field simulation geological model is set up, can
With the result data according to the target reservoir, divided by map grids, set up the stress field simulation ground of the target reservoir
Matter model.Wherein, above-mentioned map grids are divided, and can be specifically to carry out stress and strain model according to mesh generation principle, set up three-dimensional
Regular grid.
Step 13：According to the result data and target reservoir well-log information of the target reservoir, by the stress field mould
Intend geological model, be calculated target reservoir Stress Field Distribution data.
In the present embodiment, according to practical study demand and being embodied as situation, target reservoir stress described here
Field distribution data can specifically include data below：The shearing strength of the rock of each measuring point, the direct stress of each measuring point is
The shearing strength of rock when zero, the plane of fracture shear stress of each measuring point, the coefficient of internal friction of each measuring point, each measuring point
Crack shear fracture rate, the value of shearing of each measuring point, the crack tension fracture rate of each measuring point, the plane of fracture tensile stress of each measuring point
Size, the tensile strength of rock of each measuring point, the tension fissure ratio of each measuring point, the shear crack ratio of each measuring point, each survey
Gash fracture rate distributed constant value of point etc..
In one preferred embodiment, in order to make full use of the well logging of the result data of target reservoir, target reservoir
Data and the stress field simulation geological model having built up, obtain target reservoir Stress Field Distribution data.When being embodied as, can
With：According to the geological model set up, using rock mechanics parameters (Young's moduluss, Poisson's ratio in the result data of above-mentioned acquisition
And the density of rock etc.) target reservoir numerical simulation for stress field research is carried out using Finite Element Method, with SuperSAP
(Super structural analysis＆design program, structural analyses design the abbreviation of program) relative program, choosing
Take above-mentioned parameter value and be adjusted calculating, obtain target reservoir Stress field features figure.Certainly, what deserves to be explained is obtaining here
Target reservoir stress characteristics figure is intended merely to more intuitively react target reservoir stress situation, is easy to this area research worker to read
The related stress data of analysis.The application is not restricted for whether drawing target reservoir Stress field features figure.During construction, can be with
As the case may be, real needs, flexibly choose whether to draw target reservoir Stress field features figure.
In the above-described embodiment, the SuperSAP for being used is a kind of according to structure finite element establishment, for big
Volume or large area object carry out a kind of program software of stress field or Temperature calculating.Certainly, it should be noted that the application
When being embodied as, it is possible to use SuperSAP assists calculating to try to achieve target reservoir Stress Field Distribution data, it is also possible to use other
The similar software assistance obtains required target reservoir Stress Field Distribution data.In this regard, the application is not construed as limiting.
Step 14：According to the target reservoir Stress Field Distribution data, the target reservoir well-log information and target reservoir
Well test data, solution obtain correcting permeability.
In one preferred embodiment, provide to make full use of the result data of target reservoir, target reservoir to log well
The target reservoir Stress Field Distribution data of material sum, obtain in target reservoir that each fracture development is close in multiple fracture development density regions
The correction permeability of each measuring point in degree area.Specifically can perform according to the following steps：
S1：According to the target reservoir Stress Field Distribution data, gash fracture rate distributed constant value is calculated.
S2：According to the target reservoir Stress Field Distribution data and the gash fracture rate distributed constant value, it is calculated
Fracture spacing predictive value.
S3：According to the fracture spacing predictive value, the target reservoir is divided into into multiple fracture development density regions.
S4：According to the target reservoir well-log information and well test data, solution obtains the plurality of fracture development density region
In in each fracture development density region each measuring point correction permeability.
In one preferred embodiment, in order to be calculated rupture rate according to target reservoir Stress Field Distribution data,
Specifically can solve according to below equation, obtain the distributed constant value of the rupture rate of each measuring point：
In(j)=τn(j)/[τn](j)
It(j)=σt(j)/[σt](j)
I(j)=α(j)It(j)+β(j)In(j)
Wherein, [τn](j)To number the shearing strength of the rock of the measuring point for being j, C(j)To number the direct stress of the measuring point for being j
The shearing strength of rock, σ when being zeron(j)To number the plane of fracture shear stress of the measuring point for being j,It is j's for numbering
The coefficient of internal friction of measuring point, In(j)To number the crack shear fracture rate of the measuring point for being j, τn(j)To number the shear stress of the measuring point for being j
Value, It(j)To number the crack tension fracture rate of the measuring point for being j, σt(j)To number the plane of fracture tensile stress size of the measuring point for being j,
[σt](j)To number the tensile strength of rock of the measuring point for being j, α(j)To number the tension fissure ratio of the measuring point for being j, β(j)For numbering
For the shear crack ratio of the measuring point of j, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j.
In one preferred embodiment, in order to according to target reservoir Stress Field Distribution data and the distribution of gash fracture rate
Parameter value, is calculated fracture spacing predictive value, specifically can be solved according to equation below, obtains each measuring point fracture spacing
Predictive value：
I(j)＜ Ic,ξ(j)=a1(j)I(j) 2+a2(j)
I(j)≥Ic,ξ(j)=a1(j)I(j) 2+a2(j)+a3(j)(A(j)+a4(j))
Wherein, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j, IcFor I(j)Marginal value, ξ(j)To compile
Number for j measuring point fracture spacing predictive value, a1(j)、a2(j)、a3(j)、a4(j)To number the undetermined constant of the measuring point for being j, A(j)For
Numbering is the energy value of the measuring point of j.
In one preferred embodiment, in order to according to target reservoir well-log information and well test data, solution is obtained
State the correction permeability of each measuring point in each fracture development density region in multiple fracture development density regions, specifically can according to
Lower step is performed：
S4_1：According to the target reservoir well-log information, during solution obtains the plurality of fracture development density region, each splits
The meansigma methodss of the well logging permeability in seam development density region.
S4_2：According to the target reservoir well test data, during solution obtains the plurality of fracture development density region, each splits
Well testing permeability in seam development density region.
S4_3：According to the flat of the well logging permeability in each fracture development density region in the plurality of fracture development density region
Average and well testing permeability value, solution obtain the correction in the plurality of fracture development density region in each fracture development density region
Coefficient.
S4_4：According to each fracture development in the well-log information and the plurality of fracture development density region of the target reservoir
Correction coefficient in density region, solution obtain each survey in each fracture development density region in the plurality of fracture development density region
The correction permeability of point.
In one preferred embodiment, in order to according to target reservoir well-log information, solution obtains above-mentioned multiple cracks
The meansigma methodss of the well logging permeability in development density region in each fracture development density region, specifically can be according to the target reservoir
Well-log information, according to below equation, solution is obtained in the plurality of fracture development density region in each fracture development density region
The meansigma methodss of well logging permeability：
Wherein,For the well logging permeability meansigma methodss in the F fracture spacing development density region in target reservoir, K
(φ(j),MZ(j),Vsh(j))(j)Well logging for the measuring point that the numbering in the i-th substratum in the F fracture spacing development density region is j is oozed
Saturating rate, φ(j)、MZ(j)、Vsh(j)、VshNumbering in respectively the F fracture spacing development density region in the i-th substratum is the measuring point of j
Porosity, median grain diameter, shale content, hi(F)For the i-th small thickness in the F fracture spacing development density region, h(F)For
Oil reservoir actual production interval thickness in F fracture spacing development density region, niFor in the F fracture spacing development density region i-th
The measuring point number included in substratum, N(F)For the little number of plies included in the F fracture spacing development density region.
In one preferred embodiment, in order to make full use of target reservoir well test data, above-mentioned multiple cracks are obtained
Well testing permeability in development density region in each fracture development density region, specifically can provide according to the target reservoir well testing
Material, according to below equation, solution obtains the well testing in the plurality of fracture development density region in each fracture development density region and oozes
Saturating rate：
Wherein, KT(F)For the well testing permeability in the F fracture spacing development density region in target reservoir, Q(F)For target
Well test rate in oil reservoir in the F fracture spacing development density region, μ(F)For the ground in the F fracture spacing development density region
Layer viscosity of crude, re(F)For the drainage radius in the F fracture spacing development density region, rw(F)It is close for the F fracture spacing development
Oil well radius in degree area, Δ P(F)For the producing pressure differential in the F fracture spacing development density region, h(F)It is close for the F crack
Well testing reservoir thickness in degree development density region.
In one preferred embodiment, in order to close according to each fracture development in above-mentioned multiple fracture development density regions
The meansigma methodss and well testing permeability value of the well logging permeability in degree area, obtain each crack in above-mentioned multiple fracture development density regions
Correction coefficient in development density region, specifically can be according to the survey in each fracture development density region in multiple fracture development density regions
The meansigma methodss and well testing permeability value of well permeability, according to below equation, solution is obtained in the plurality of fracture development density region
Correction coefficient in each fracture development density region：
Wherein, the λ(F)For correction coefficient in the F fracture spacing development density region, KT(F)Send out for the F fracture spacing
Well testing permeability in density region is educated,For permeability meansigma methodss of logging well in the F fracture spacing development density region.
In one preferred embodiment, in order to the well-log information according to target reservoir and the plurality of fracture development are close
Correction coefficient in each fracture development density region in degree area, in obtaining the plurality of fracture development density region, each fracture development is close
The correction permeability of each measuring point in degree area.Specifically can be close according to the well-log information of target reservoir and the plurality of fracture development
Correction coefficient in each fracture development density region in degree area, according to below equation, solution obtains the plurality of fracture development density
In area in each fracture development density region each measuring point correction permeability：
Wherein, the PERM(j)In target reservoir numbering be j measuring points correction permeability, λ(F)It is that j measuring points are located for numbering
The F fracture spacing development density region in correction coefficient, K (φ(j),MZ(j),Vsh(j)) it is that the well logging for numbering the measuring point for being j is oozed
Saturating rate.
Step 15：According to the correction permeability, the three-dimensional penetration rate model of the target reservoir is set up.
In one preferred embodiment, in order to set up the three-dimensional penetration rate model of target reservoir, specifically can basis
Permeability after correction, using Geologic modeling software, is analyzed by variogram, sets up target reservoir three-dimensional penetration rate model.
It should be noted that the height of permeability can be seen substantially with crack by setting up in the three-dimensional penetration rate model for obtaining
Density size and show distributional difference spatially.
In the embodiment of the present application, compared to existing modeling method, the application is by logging well according to well test data correction
Permeability, obtains correcting permeability, sets up three-dimensional penetration rate model further according to correction permeability, so as to solve existing modeling side
The three-dimensional technical problem that penetration rate model precision is low, error is big set up present in method, realizes raising and sets up three-dimensional and ooze
The technique effect of saturating rate model accuracy.
Based on same inventive concept, a kind of device for setting up three-dimensional penetration rate model in the embodiment of the present invention, is additionally provided,
As described in the following examples.Due to the principle of the device solve problem it is similar to the method for setting up three-dimensional penetration rate model, because
The enforcement of this device may refer to the enforcement of the method for setting up three-dimensional penetration rate model, repeats part and repeats no more.It is following to be made
, term " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although following examples
Described device is preferably with software realizing, but hardware, or the realization of the combination of software and hardware be also may be simultaneously
It is contemplated.Fig. 2 is referred to, is a kind of composition structure chart of the device of the three-dimensional penetration rate model of foundation of the embodiment of the present application,
The device can include：Acquisition module 201, first are set up module 202, computing module 203, correction module 204 and second and set up
Module 205, is specifically described to the structure below.
Acquisition module 201, for obtaining the result data of target reservoir, wherein, the result data includes：Rock mechanics
Relevant parameter, stress field size and level of stress field biggest principal stress direction；
First sets up module 202, for the result data according to the target reservoir, sets up the stress of the target reservoir
Field stimulation geological model；
Computing module 203, for the result data according to the target reservoir and target reservoir well-log information, by described
Stress field simulation geological model, is calculated target reservoir Stress Field Distribution data；
Correction module 204, for according to the target reservoir Stress Field Distribution data, the target reservoir well-log information and
Target reservoir well test data, solution obtain correcting permeability；
Second sets up module 205, for according to the correction permeability, setting up the three-dimensional permeability mould of the target reservoir
Type.
In a specific application scenarios, the method for the three-dimensional penetration rate model of foundation for for example proposing using the application/
Device carries out the analysis and research of target reservoir to somewhere, specifically can perform according to the following steps.
S1：Carry out stress field simulation, prediction target reservoir fracture development and the regularity of distribution
S1-1, admission target reservoir actual well drilled rock core, carry out rock triaxial stress test, obtain crucial rock mechanics
Relevant parameter is (rock porosity, rock density, tensile strength, Young's moduluss, Poisson's ratio, divergence cone angle, internal friction angle, cohesion
Power)；Carry out Acoustic Emission of Rock experiment, obtain stress field stress intensity.
S1-2, analysis target reservoir well-log information, obtain level of stress field biggest principal stress direction.
S1-3, integration objective reservoir geology achievement in research and experimental result, are divided using map grids, set up three-dimensional regular
Grid, so as to set up the geological model of target reservoir stress field simulation.
S1-4, according to set up geological model, using above-mentioned acquisition rock mechanics parameters (Young's moduluss, Poisson's ratio with
And the density of rock etc.) target reservoir numerical simulation for stress field research is carried out using Finite Element Method, it is related with SuperSAP
Program, chooses above-mentioned parameter value and is adjusted calculating, obtain target reservoir Stress field features figure.
S1-5, calculating gash fracture rate distributed constant value
The gash fracture rate distributed constant value of different drill cores is calculated one by one using formula 1, and using the geology set up
Model, by the method for the stochastic modeling of Geologic modeling software, works out gash fracture rate distributed constant Distribution value prognostic chart.Specifically
Fig. 3 can be referred to, stain well location numbering different with digital expression in figure.
Total gash fracture rate distributed constant value：IAlways=α It+βIn
In formula：A, β are the tension fissure and shear crack ratio of core observation statistics.
Crack shear fracture rate：In=τn/[τn]
In formula：τnFor value of shearing, unit：MPa；[τn] for the shearing strength of rock, unit：MPa.
In formula：C is cohesive strength, is the shearing strength of direct stress rock when being zero, unit：MPa；For internal friction angle, list
Position：°；For coefficient of internal friction.
Crack tension fracture rate：It=σt/[σt]
In formula：σtFor plane of fracture tensile stress size；[σt] for tensile strength of rock.
S1-6, calculating prediction fracture spacing
According to binary method prediction fracture spacing method (formula 2) prediction target area fracture spacing, and made according to fracture spacing
Fracture spacing forecast of distribution figure, can specifically refer to Fig. 4.
I ＜ Ic, β=a1I2+a2
I ＞ Ic, β=a1I2+a2+a3(A+a4)
In formula：β is fracture spacing predictive value；I is rupture value；a1, a2, a3, a4For four undetermined constants；IcFor rupture value I
Marginal value；A-energy value.
S2：Carry out reservoir performance analysis, obtain target reservoir Dynamic Permeability
Analysis target reservoir dynamic monitoring information, especially well test data, using percolation hydraulic theory, try to every mouth well
Well data carries out well test analysis (formula 3), objective of interpretation oil reservoir every mouth well performance permeability.
In formula：KTFor Reservoir behavior permeability, unit：μm2；Q is oil well output, unit：cm3/ s (underground value)；μ is ground
Layer viscosity of crude, unit：mPa·s；reFor drainage radius, unit：cm；rwFor oil well radius, unit：cm；Δ P is production pressure
Difference, unit：MPa, h are reservoir thickness, unit：cm.
S3：Well testing permeability is set up with well logging permeability corresponding relation
It is according to the different fracture development density regions set up and well testing well location corresponding relation, close according to different fracture developments
Degree area, sets up different well testing permeabilities and well logging permeability corresponding relation, and then obtains well testing permeability calibration model.Specifically
Process is as follows：
S-1, the permeability that oil reservoir actual production interval is calculated using the mean value method of well logging permeability
In formula：For oil reservoir actual production interval well testing permeability meansigma methodss, unit：μm2；K(φ,MZ,Vsh) for i-th
The well logging permeability of substratum jth sampled point, unit：μm2；φ、Mz、VshThe respectively porosity of the i-th substratum jth sampled point, grain
Degree intermediate value, shale content；hiFor the i-th small thickness, unit：m；H be oil reservoir actual production interval thickness, unit：m；n1For i-th
The well logging permeability points included in substratum；n2For the little number of plies is included in pressure measurement interval.
S3-2, the calibration model in order to set up well testing permeability and well logging permeability, introduce correction coefficient λ, set up different
Well testing permeability with well logging permeability corresponding relation.
S3-3, using correction coefficient λ, fracture density development not same district, using different correction coefficients, is set up well logging and is oozed
The calibration model of saturating rate.
Table 1 can be specifically referred to, table 1 is the school between the Dynamic Permeability and well logging permeability obtained using said method
Positive coefficient.It is very different with well logging permeability by Dynamic Permeability is can be seen that in table.Wherein, both difference maximums is reachable
47 times, minimum also has 6 times, typically all changes between 10～25 times.
1 well log interpretation permeability dynamic calibration coefficient of table
Note：Oil volume factor is calculated in table and takes 1.112, viscosity of crude takes 6.7mPas.
Reservoir behavior test permeability can reflect the seepage flow characteristics of oil reservoir, closer to the true permeability of oil reservoir.Using life
Dynamic Permeability is produced to permeability timing of logging well, well logging permeability (PERM) dynamic calibration model can be set up, i.e.,：
PERM=λ × K (φ, Mz,Vsh)
S4：Using the well logging permeability of correction, permeability threedimensional model is set up.
According to the permeability after correction, using Geologic modeling software, analyzed by variogram, set up target reservoir three-dimensional
Penetration rate model, can see that in three-dimensional penetration rate model the height of permeability is substantially showed with the density size in crack
Go out distributional difference spatially.By being embodied as, can mainly show as at following 2 points：
1st, in the area that local cracks density is larger, permeability is also higher, and this is mainly crack and permeability is generated very
Big contribution；
2nd, the less regional reservoir permeability of fracture spacing is also less, is mainly affected by matrix permeability, less crack
Also less is contributed to local permeability.
The impact to permeating that fracture spacing is produced is coincide with local area exploitation reality very much.In the big block of fracture spacing,
Initial stage of development Yield compari@is good, and development late stage water logging is extremely serious, and this is that due to higher in these block permeabilities, fluid is along splitting
Seam can smoothly flow into pit shaft, but due to a large amount of presence in crack, injection water is easy to along crack tonguing, so as to result in oil well
Water breakthrough is very fast, the phenomenon that later stage water logging is serious.And in the little relatively low-permeability block of fracture development density, although oilwell produced fluid amount
Relatively low, but rate of water cut increase is also relatively slow, injection water phenomenon of advancing by leaps and bounds is not obvious.This not only absolutely proves built fractured model
Correctness.Also having embodied impact of the crack to ultra-permeable reservior local penetration rate can be fairly obvious.Therefore, it can be seen that sharp
The method or apparatus of the three-dimensional penetration rate model of foundation provided with the application can more accurately set up reservoir permeability model,
Reliable basis are provided for oil reservoir development.
As can be seen from the above description, the method for the three-dimensional penetration rate model of foundation that the embodiment of the present application is provided and
Device first obtains the result data of target reservoir compared to existing modeling method, and setting up target reservoir according to result data should
The geological model of field of force simulation, obtains supplemental characteristic by the geological model of target reservoir stress field simulation；In view of in reservoir
Impact of the microcrack to true permeability, is corrected to permeability of logging well according to supplemental characteristic and with reference to well test data, obtains
Correction permeability；Three-dimensional penetration rate model is set up according to correction permeability.The application is by setting up three-dimensional according to correction permeability
Penetration rate model, existing sets up according to well logging permeability that the precision that three-dimensional penetration rate model is present is low, error is big so as to solve
Technical problem, realize the technique effect for improving the accuracy for setting up three-dimensional penetration rate model；Additionally, the application is also by dividing
Correction coefficient that Ji Suan be in different fracture spacing development areas, it is close to difference according to the correction coefficient of different fracture spacing development areas
Well logging permeability in degree development area is targetedly corrected, and is further increased and is set up the accurate of three-dimensional penetration rate model
Degree.
Although this application provides the method operating procedure as described in embodiment or flow chart, is based on routine or noinvasive
The means of the property made can include more or less operating procedures.The step of enumerating in embodiment order is only numerous steps
A kind of mode in execution sequence, does not represent unique execution sequence.When device or client production in practice is performed, can
To perform or executed in parallel (such as at parallel processor or multithreading according to embodiment or method shown in the drawings order
The environment of reason, even distributed data processing environment).Term " including ", "comprising" or its any other variant are intended to contain
Lid nonexcludability is included, so that a series of process, method, product or equipment including key elements not only will including those
Element, but also including other key elements being not expressly set out, or also include for this process, method, product or equipment
Intrinsic key element.In the absence of more restrictions, it is not excluded that including the process of the key element, method, product or
Also there are other identical or equivalent elements in person's equipment.
Device that above-described embodiment is illustrated or module etc., specifically can be realized by computer chip or entity, or by having
There is the product of certain function to realize.For convenience of description, it is divided into various modules with function when describing apparatus above to retouch respectively
State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when the application is implemented,
The module for realizing same function can be realized by the combination of multiple submodule etc..Device embodiment described above is only
Schematically, for example, the division of the module, only a kind of division of logic function can have other drawing when actually realizing
Point mode, such as multiple module or components can with reference to or be desirably integrated into another system, or some features can be ignored,
Or do not perform.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, it is complete
Entirely can by by method and step carry out programming in logic cause controller with gate, switch, special IC, may be programmed
The form of logic controller and embedded microcontroller etc. is realizing identical function.Therefore this controller is considered one kind
Hardware component, and the device for realizing various functions included to its inside can also be considered as the structure in hardware component.Or
Person even, can be used for realizing that the device of various functions be considered as not only being the software module of implementation method but also being hardware
Structure in part.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes execution particular task or realizes the routine of particular abstract data type, program, object, group
Part, data structure, class etc..The application is put into practice in a distributed computing environment can also, in these distributed computing environment,
Task is performed by the remote processing devices connected by communication network.In a distributed computing environment, program module can
With positioned at including in the local and remote computer-readable storage medium including storage device.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.The application can be used for crowd
In more general or special purpose computing system environments or configuration.For example：Personal computer, server computer, handheld device or
Portable set, laptop device, multicomputer system, set based on the system of microprocessor, set top box, programmable electronics
Standby, network PC, minicomputer, mainframe computer, including the distributed computing environment etc. of any of the above system or equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application.
Claims (13)
1. a kind of method for setting up three-dimensional penetration rate model, it is characterised in that include：
Obtain the result data of target reservoir；
According to the result data of the target reservoir, the stress field simulation geological model of target reservoir is set up；
According to the well-log information of the result data and target reservoir of the target reservoir, by the stress field simulation Geological Model
Type, is calculated target reservoir Stress Field Distribution data；
Provided according to the well testing of the target reservoir Stress Field Distribution data, the well-log information of the target reservoir and target reservoir
Material, solution obtain correcting permeability；
According to the correction permeability, the three-dimensional penetration rate model of target reservoir is set up.
2. method according to claim 1, it is characterised in that the result data includes：Rock mechanics relevant parameter, should
Field of force size and level of stress field biggest principal stress direction.
3. method according to claim 2, it is characterised in that the result data of the acquisition target reservoir, including：
Obtain the drill cores sample of target reservoir；
Rock triaxial stress test is carried out to the drill cores sample of the target reservoir, the related ginseng of the rock mechanics is obtained
Number；
Acoustic Emission of Rock experiment is carried out to the drill cores sample of the target reservoir, the stress field size is obtained；
According to the well-log information of the target reservoir, the level of stress field biggest principal stress direction is obtained.
4. method according to claim 2, it is characterised in that according to the result data of the target reservoir, set up target
The stress field simulation geological model of oil reservoir, including：
According to the result data of the target reservoir, divided by map grids, set up the stress field simulation of the target reservoir
Geological model.
5. method according to claim 2, it is characterised in that according to target reservoir Stress Field Distribution data, described
The well test data of the well-log information and target reservoir of target reservoir, solution obtain correcting permeability, including：
According to the target reservoir Stress Field Distribution data, the gash fracture rate distribution of each measuring point in target reservoir is calculated
Parameter value；
According to the gash fracture rate distribution ginseng of each measuring point in the target reservoir Stress Field Distribution data and the target reservoir
Numerical value, is calculated the fracture spacing predictive value of each measuring point in target reservoir；
According to the fracture spacing predictive value of each measuring point in the target reservoir, the target reservoir is divided into into multiple cracks and is sent out
Educate density region；
According to the well test data of the well-log information and the target reservoir of the target reservoir, solution is obtained in the target reservoir
The correction permeability of each measuring point.
6. method according to claim 5, it is characterised in that described according to the target reservoir Stress Field Distribution data,
The gash fracture rate distributed constant value of each measuring point in target reservoir is calculated, including：
According to the target reservoir Stress Field Distribution data, according to below equation, each survey in the target reservoir is calculated
The gash fracture rate distributed constant value of point：
In(j)=τn(j)/[τn](j)
It(j)=σt(j)/[σt](j)
I(j)=α(j)It(j)+β(j)In(j)
Wherein, [τn](j)To number the shearing strength of the rock of the measuring point for being j, C(j)Direct stress for the measuring point that numbering is j is zero
The shearing strength of Shi Yanshi, σn(j)To number the plane of fracture shear stress of the measuring point for being j,For the measuring point that numbering is j
Coefficient of internal friction, In(j)To number the crack shear fracture rate of the measuring point for being j, τn(j)To number the value of shearing of the measuring point for being j,
It(j)To number the crack tension fracture rate of the measuring point for being j, σt(j)To number the plane of fracture tensile stress size of the measuring point for being j, [σt](j)
To number the tensile strength of rock of the measuring point for being j, α(j)To number the tension fissure ratio of the measuring point for being j, β(j)For the survey that numbering is j
The shear crack ratio of point, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j.
7. method according to claim 5, it is characterised in that it is described according to the target reservoir Stress Field Distribution data and
In the target reservoir, the gash fracture rate distributed constant value of each measuring point, is calculated the crack of each measuring point in target reservoir
Density prediction value, including：
According to the gash fracture rate distribution ginseng of each measuring point in the target reservoir Stress Field Distribution data and the target reservoir
Numerical value, according to below equation, is calculated the fracture spacing predictive value of each measuring point in the target reservoir：
I(j)＜ Ic,ξ(j)=a1(j)I(j) 2+a2(j)
I(j)≥Ic,ξ(j)=a1(j)I(j) 2+a2(j)+a3(j)(A(j)+a4(j))
Wherein, I(j)To number the gash fracture rate distributed constant value of the measuring point for being j, IcFor I(j)Marginal value, ξ(j)It is j for numbering
Measuring point fracture spacing predictive value, a1(j)、a2(j)、a3(j)、a4(j)To number the undetermined constant of the measuring point for being j, A(j)For numbering
For the energy value of the measuring point of j.
8. method according to claim 5, it is characterised in that the well-log information according to the target reservoir and described
The well test data of target reservoir, solution obtain the correction permeability of each measuring point in the target reservoir, including：
According to the well-log information of the target reservoir, during solution obtains the plurality of fracture development density region, each fracture development is close
The meansigma methodss of the well logging permeability in degree area；
According to the target reservoir well test data, solution obtains each fracture development density in the plurality of fracture development density region
Well testing permeability in area；
According to meansigma methodss and the examination of the well logging permeability in each fracture development density region in the plurality of fracture development density region
Well permeability value, solution obtain in the plurality of fracture development density region correction coefficient in each fracture development density region；
According in each fracture development density region in the well-log information and the plurality of fracture development density region of the target reservoir
Correction coefficient, solution obtain the correction permeability of each measuring point in the target reservoir.
9. method according to claim 8, it is characterised in that according to the well-log information of the target reservoir, solution is obtained
The meansigma methodss of the well logging permeability in the plurality of fracture development density region in each fracture development density region, including：
According to the well-log information of the target reservoir, according to below equation, solution is obtained in the plurality of fracture development density region
The meansigma methodss of the well logging permeability in each fracture development density region：
Wherein,For the well logging permeability meansigma methodss in the F fracture spacing development density region in target reservoir, K (φ(j),
MZ(j),Vsh(j))(j)The numbering developed in density region in the i-th substratum for the F fracture spacing is the well logging permeability of the measuring point of j,
φ(j)、MZ(j)、Vsh(j)、VshNumbering in respectively the F fracture spacing development density region in the i-th substratum is the hole of the measuring point of j
Porosity, median grain diameter, shale content, hi(F)For the i-th small thickness in the F fracture spacing development density region, h(F)For F
Oil reservoir actual production interval thickness in fracture spacing development density region, niFor the i-th substratum in the F fracture spacing development density region
Interior included measuring point number, N(F)For the little number of plies included in the F fracture spacing development density region.
10. method according to claim 8, it is characterised in that according to the well test data of the target reservoir, solution is obtained
Well testing permeability in the plurality of fracture development density region in each fracture development density region, including：
According to the well test data of the target reservoir, according to below equation, solution is obtained in the plurality of fracture development density region
Well testing permeability in each fracture development density region：
Wherein, KT(F)For the well testing permeability in the F fracture spacing development density region in target reservoir, Q(F)For in target reservoir
Well test rate in the F fracture spacing development density region, μ(F)For the in-place oil in the F fracture spacing development density region
Viscosity, re(F)For the drainage radius in the F fracture spacing development density region, rw(F)For in the F fracture spacing development density region
Oil well radius, Δ P(F)For the producing pressure differential in the F fracture spacing development density region, h(F)For the F fracture spacing development
Well testing reservoir thickness in density region.
11. methods according to claim 8, it is characterised in that according in the plurality of fracture development density region each split
The meansigma methodss and well testing permeability value of the well logging permeability in seam development density region, solution obtain the plurality of fracture development density
Correction coefficient in each fracture development density region in area, including：
According to meansigma methodss and the examination of the well logging permeability in each fracture development density region in the plurality of fracture development density region
Well permeability value, according to below equation, solution is obtained in the plurality of fracture development density region in each fracture development density region
Correction coefficient：
Wherein, the λ(F)For correction coefficient in the F fracture spacing development density region, KT(F)It is close for the F fracture spacing development
Well testing permeability in degree area,For permeability meansigma methodss of logging well in the F fracture spacing development density region.
12. methods according to claim 8, it is characterised in that according to the well-log information of the target reservoir and described many
Correction coefficient in each fracture development density region in individual fracture development density region, solution obtain each measuring point in the target reservoir
Correction permeability, including：
According in each fracture development density region in the well-log information and the plurality of fracture development density region of the target reservoir
Correction coefficient, according to below equation, solution obtains the correction permeability of each measuring point in the target reservoir：
Wherein, the PERM(j)In target reservoir numbering be j measuring points correction permeability, λ(F)For numbering be j measuring points be located the
Correction coefficient in F fracture spacing development density region, K (φ(j),MZ(j),Vsh(j)) it is the well logging permeability for numbering the measuring point for being j.
13. a kind of devices for setting up three-dimensional penetration rate model, it is characterised in that include：
Acquisition module, for obtaining the result data of target reservoir；
First sets up module, for the result data according to the target reservoir, sets up the stress field simulation of the target reservoir
Geological model；
Computing module, for the well-log information of result data and target reservoir according to the target reservoir, by the stress
Field stimulation geological model, is calculated target reservoir Stress Field Distribution data；
Correction module, for according to the target reservoir Stress Field Distribution data, the target reservoir well-log information and target oil
Well test data is hidden, solution obtains the correction permeability of each measuring point in target reservoir；
Second sets up module, for the correction permeability according to each measuring point in the target reservoir, sets up the target reservoir
Three-dimensional penetration rate model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610944944.3A CN106546525A (en) | 2016-11-02 | 2016-11-02 | The method and apparatus for setting up three-dimensional penetration rate model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610944944.3A CN106546525A (en) | 2016-11-02 | 2016-11-02 | The method and apparatus for setting up three-dimensional penetration rate model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106546525A true CN106546525A (en) | 2017-03-29 |
Family
ID=58393659
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610944944.3A Pending CN106546525A (en) | 2016-11-02 | 2016-11-02 | The method and apparatus for setting up three-dimensional penetration rate model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106546525A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107132170A (en) * | 2017-04-10 | 2017-09-05 | 中国石油天然气股份有限公司 | A kind of determination method and apparatus of reservoir stress sensitivity |
| CN107153038A (en) * | 2017-06-08 | 2017-09-12 | 铁道第三勘察设计院集团有限公司 | Stratum osmotic coefficient quickly determines probe and its application method |
| CN107366534A (en) * | 2017-08-10 | 2017-11-21 | 中国石油天然气股份有限公司 | It is roughened the determination method and apparatus of permeability |
| CN108843300A (en) * | 2018-06-25 | 2018-11-20 | 中国石油天然气股份有限公司 | It is a kind of complexity porous media in determine mainstream channel type method and device |
| CN109145375A (en) * | 2018-07-18 | 2019-01-04 | 中国石油天然气股份有限公司 | The bearing calibration of fracture permeabgility model and device |
| CN112818501A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for correcting carbonate reservoir static permeability based on dynamic monitoring data |
| CN113586042A (en) * | 2021-07-09 | 2021-11-02 | 中国石油天然气股份有限公司 | Method and device for determining dynamic permeability of porous carbonate rock |
-
2016
- 2016-11-02 CN CN201610944944.3A patent/CN106546525A/en active Pending
Non-Patent Citations (3)
| Title |
|---|
| YONG-CHAO XUE ET AL.: "A new method of building permeability model in low-permeability reservoir numerical simulation", 《NATURAL SCIENCE》 * |
| 薛永超: "肯基亚克盐下油田应力场模拟及裂缝预测", 《复杂油气藏》 * |
| 陈子光: "《岩石力学性质与构造应力场》", 30 June 1986 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107132170B (en) * | 2017-04-10 | 2019-09-06 | 中国石油天然气股份有限公司 | A kind of determination method and apparatus of reservoir stress sensitivity |
| CN107132170A (en) * | 2017-04-10 | 2017-09-05 | 中国石油天然气股份有限公司 | A kind of determination method and apparatus of reservoir stress sensitivity |
| CN107153038A (en) * | 2017-06-08 | 2017-09-12 | 铁道第三勘察设计院集团有限公司 | Stratum osmotic coefficient quickly determines probe and its application method |
| CN107366534A (en) * | 2017-08-10 | 2017-11-21 | 中国石油天然气股份有限公司 | It is roughened the determination method and apparatus of permeability |
| CN107366534B (en) * | 2017-08-10 | 2020-08-11 | 中国石油天然气股份有限公司 | Method and device for determining coarsening permeability |
| CN108843300A (en) * | 2018-06-25 | 2018-11-20 | 中国石油天然气股份有限公司 | It is a kind of complexity porous media in determine mainstream channel type method and device |
| CN108843300B (en) * | 2018-06-25 | 2022-03-01 | 中国石油天然气股份有限公司 | Method and device for determining type of main flow channel in complex porous medium |
| CN109145375A (en) * | 2018-07-18 | 2019-01-04 | 中国石油天然气股份有限公司 | The bearing calibration of fracture permeabgility model and device |
| CN109145375B (en) * | 2018-07-18 | 2022-02-01 | 中国石油天然气股份有限公司 | Crack permeability model correction method and device |
| CN112818501A (en) * | 2019-11-15 | 2021-05-18 | 中国石油天然气股份有限公司 | Method for correcting carbonate reservoir static permeability based on dynamic monitoring data |
| CN112818501B (en) * | 2019-11-15 | 2022-11-04 | 中国石油天然气股份有限公司 | Method for correcting carbonate reservoir static permeability based on dynamic monitoring data |
| CN113586042A (en) * | 2021-07-09 | 2021-11-02 | 中国石油天然气股份有限公司 | Method and device for determining dynamic permeability of porous carbonate rock |
| CN113586042B (en) * | 2021-07-09 | 2023-09-26 | 中国石油天然气股份有限公司 | Method and device for determining dynamic permeability of pore type carbonate rock |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106546525A (en) | The method and apparatus for setting up three-dimensional penetration rate model | |
| Zhang et al. | Investigating hydraulic fracturing complexity in naturally fractured rock masses using fully coupled multiscale numerical modeling | |
| Lei et al. | A new approach to upscaling fracture network models while preserving geostatistical and geomechanical characteristics | |
| Shahbazi et al. | A review of existing methods used to evaluate the hydraulic conductivity of a fractured rock mass | |
| US10961826B2 (en) | Subsurface modeler workflow and tool | |
| CN104677734A (en) | Multi-scale information integration tight formation brittleness index measuring method | |
| CN103913772B (en) | Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter | |
| Fu et al. | Predicting hydraulic fracture trajectory under the influence of a mine drift in EGS Collab Experiment I | |
| CN104865610B (en) | A kind of numerical representation method method of deep reservoir rock core | |
| Zhou et al. | Numerical modeling and investigation of fracture propagation with arbitrary orientation through fluid injection in tight gas reservoirs with combined XFEM and FVM | |
| CN109655903A (en) | Rammell S-Wave Velocity Predicted Method and system | |
| Suo et al. | Experimental study on mechanical and anisotropic properties of shale and estimation of uniaxial compressive strength | |
| White et al. | Code comparison study fosters confidence in the numerical simulation of enhanced geothermal systems | |
| Earnest et al. | Investigating the role of hydromechanical coupling on flow and transport in shallow fractured-rock aquifers | |
| Burnell et al. | Geothermal supermodels: the next generation of integrated geophysical, chemical and flow simulation modelling tools | |
| Wang et al. | Gas shale hydraulic fracturing: a numerical investigation of the fracturing network evolution in the Silurian Longmaxi formation in the southeast of Sichuan Basin, China, using a coupled FSD approach | |
| Cornet et al. | Nonlinear viscoelastic closure of salt cavities | |
| Zhou et al. | Optimizing the construction parameters of modified zipper fracs in multiple horizontal wells | |
| CN105929461B (en) | A kind of sound state rock mechanics parameters correction system | |
| Fernandes et al. | Analytical Model to Effective Permeability Loss Monitoring and Boundary Effects Identification in Oil Wells with Finite Hydraulic Fractures | |
| CN109211745A (en) | A kind of restoration methods of the evolutionary process rich in organic matter mud shale porosity | |
| CN102436014B (en) | Method for evaluating multi-parameter three-dimensional geologic structure of underground water sealed cave depot | |
| CN116738794A (en) | Rock physical numerical simulation method, device, equipment and medium for pore fracture medium | |
| CN110863814A (en) | Method and device for determining single-well subsection specific fluid production index of ultra-thick type cuttings limestone oil reservoir | |
| Andersson et al. | Site investigations: Strategy for rock mechanics site descriptive model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170329 |