CN110244369A - Reservoir constraint and movable fluid distribution determination method, apparatus and system - Google Patents
Reservoir constraint and movable fluid distribution determination method, apparatus and system Download PDFInfo
- Publication number
- CN110244369A CN110244369A CN201910572371.XA CN201910572371A CN110244369A CN 110244369 A CN110244369 A CN 110244369A CN 201910572371 A CN201910572371 A CN 201910572371A CN 110244369 A CN110244369 A CN 110244369A
- Authority
- CN
- China
- Prior art keywords
- fluid
- relaxation time
- constraint
- distributed data
- lateral relaxation
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Abstract
This specification embodiment discloses a kind of constraint of reservoir and movable fluid distribution determination method, apparatus and system, the method includes obtaining the nuclear magnetic resonance echo data in target work area, the original lateral relaxation time distributed data that inverting obtains the target work area is carried out to the echo data；It whether include that short relaxation peak determines type belonging to the original lateral relaxation time distributed data according to the original lateral relaxation time distributed data；The original lateral relaxation time distributed data is handled according to the corresponding proportionality coefficient computation model of the type, obtains the proportionality coefficient of constraint fluid and movable fluid；The proportionality coefficient for being utilized respectively the constraint fluid and movable fluid handles the original lateral relaxation time distributed data, obtains the constraint fluid distrbution data and movable fluid distributed data in the target work area.Using each embodiment of this specification, the continuous and quantitative characterization that fluid and movable fluid distribution are fettered to reservoir can be accurately realized.
Description
Technical field
The present invention relates to log data processing technical fields in petroleum exploration and development, particularly, are related to a kind of reservoir constraint
With movable fluid distribution determination method, apparatus and system.
Background technique
Constraint fluid and movable fluid distribution can directly reflect micropore structure feature and percolation ability of rock etc.
Petrophysical property.It is effective to fetter fluid and movable fluid distribution characterization in Oil And Gas Exploration And Development, especially movably
Fluid distrbution characterization is to carry out the important evidence of evaluating reservoir, capability forecasting and reservoir Efficient Development.Therefore, Study In Reservoir beam
It ties up fluid and movable fluid distribution characterizing method is of great significance.
Low-field nuclear magnetic resonance (NMR) measures obtained lateral relaxation time (T2) it is distributed the stream that can be characterized in blowhole
Body distribution.Laboratory respectively to fluid saturated rocks and centrifugation after rock carry out NMR experiment measurement be NMR technology characterization constraint and
A kind of effective means of movable fluid distribution.NMR measurement is carried out to rock core first in the case where being saturated fluid state, obtains saturation fluid
T2Distribution.Then the movable fluid being centrifuged out under the conditions of reservoir pressure in rock core carries out the rock core under constraint fluid state
NMR measurement obtains constraint fluid T2Distribution.It is saturated fluid T2Distribution subtracts constraint fluid T2Distribution, can be obtained movable fluid T2
Distribution.But NMR centrefuge experiment can only measure limited core sample, cannot achieve to down-hole formation constraint and can
The continuous characterization of dynamic fluid distrbution.Although NMR well logging can provide down-hole formation continuous T2Distribution, but due to fluid in stratum
The complexity of distribution, there are no a kind of effective methods so far can be based on the T that NMR logs well2Distribution is accurate to realize underground
The continuous and quantitative characterization of stratum constraint and movable fluid distribution.
Summary of the invention
This specification embodiment be designed to provide a kind of constraint of reservoir and movable fluid distribution determination method, device and
System can accurately realize the continuous characterization that fluid and movable fluid distribution are fettered to reservoir.
This specification a kind of constraint of reservoir is provided and movable fluid distribution determination method, apparatus and system be include such as lower section
What formula was realized:
A kind of constraint of reservoir and movable fluid distribution determination method, comprising:
The nuclear magnetic resonance echo data for obtaining target work area carries out inverting to the echo data and obtains the target work area
Original lateral relaxation time distributed data；
It whether include that short relaxation peak determines the original transverse relaxation according to the original lateral relaxation time distributed data
Type belonging to Annual distribution data；
The original lateral relaxation time distributed data is carried out according to the corresponding proportionality coefficient computation model of the type
Processing obtains the proportionality coefficient of constraint fluid and movable fluid；
The proportionality coefficient of the constraint fluid and movable fluid is utilized respectively to the original lateral relaxation time distribution number
According to being handled, the constraint fluid distrbution data and movable fluid distributed data in the target work area are obtained.
In another embodiment of the method that this specification provides, the type includes the first kind, the second class, accordingly
, type belonging to the determination original lateral relaxation time distributed data, comprising:
If the original lateral relaxation time distributed data includes short relaxation peak, it is determined that the original lateral relaxation time
Type belonging to distributed data is the first kind；
If the original lateral relaxation time distributed data does not include short relaxation peak, it is determined that when the original transverse relaxation
Between type belonging to distributed data be the second class.
It is described according to the corresponding proportionality coefficient of the type in another embodiment of the method that this specification provides
Computation model handles the original lateral relaxation time distributed data, comprising:
If type belonging to original lateral relaxation time distributed data is the first kind, following first proportionality coefficient meters are utilized
Model is calculated to handle the original lateral relaxation time distributed data:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤ε1
Wherein, αi=i Δ α, Δ α > 0, i=1,2 ..., N, αiIndicate i-th of α value, N indicates the maximum value of α
Number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid proportionality coefficient
Function shape, pksIndicate that short relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data, ε1It is
A kind of lateral relaxation time is distributed threshold value, Pb(T2,αi) it is αiCorresponding constraint fluid proportional coefficient；
Determine the corresponding α of minimum i for meeting above-mentioned first proportionality coefficient computation modeli, it is denoted as αopt1；
By αopt1Corresponding Pb(T2,αopt1) and Pm(T2,αopt1) it is identified as constraint fluid proportional coefficient and movable stream
Body proportionality coefficient, wherein Pm(T2,αopt1)=1-Pb(T2,αopt1)。
It is described according to the corresponding proportionality coefficient of the type in another embodiment of the method that this specification provides
Computation model handles the original lateral relaxation time distributed data, comprising:
If type belonging to original lateral relaxation time distributed data is the second class, following second proportionality coefficient meters are utilized
Model is calculated to handle the original lateral relaxation time distributed data:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤ε2
Wherein, αj=(N+1-j) Δ α, Δ α > 0, j=1,2 ..., N, αjIndicate j-th of α value, N indicates the maximum of α
Value number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid ratio
The function shape of example coefficient, pklIndicate that long relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data,
ε2Threshold value, P are distributed for the second class lateral relaxation timem(T2,αi) it is αiCorresponding movable fluid proportionality coefficient；
Determine the corresponding α of minimum j for meeting above-mentioned second proportionality coefficient computation modelj, it is denoted as αopt2；
By αopt2Corresponding Pb(T2,αopt2) and Pm(T2,αopt2) it is identified as constraint fluid proportional coefficient and movable stream
Body proportionality coefficient, wherein Pb(T2,αopt2)=1-Pm(T2,αopt2)。
In another embodiment of the method that this specification provides, the ratio system of the constraint fluid and movable fluid
Number includes the exponential function using natural constant the bottom of as.
In another embodiment of the method that this specification provides, the ratio system of the constraint fluid and movable fluid
Number includes:
Wherein, Pb(T2, α) and it is constraint fluid proportional coefficient distribution function, Pm(T2, α) and it is that movable fluid proportionality coefficient is distributed
Function, α is form factor, to control Pb(T2, α) and Pm(T2, α) function shape.
This specification provide the method another embodiment in, it is described be utilized respectively the constraint fluid and movably
The proportionality coefficient function of fluid handles corresponding sorted original lateral relaxation time distributed data, comprising:
fb(T2)=Pb(T2,αopt)·f(T2)
fm(T2)=Pm(T2,αopt)·f(T2)
Wherein, fb(T2) it is constraint fluid lateral relaxation time distributed data, Pb(T2,αopt) it is constraint fluid proportional system
Number, fm(T2) it is movable fluid lateral relaxation time distributed data, Pm(T2,αopt) it is movable fluid proportionality coefficient, f (T2) it is original
Beginning lateral relaxation time distributed data.
On the other hand, this specification embodiment also provides a kind of constraint of reservoir and movable fluid is distributed determining device, described
Device includes:
Data acquisition module carries out the echo data anti-for obtaining the nuclear magnetic resonance echo data in target work area
Drill the original lateral relaxation time distributed data for obtaining the target work area；
Determination type module, for whether including that short relaxation peak determines according to the original lateral relaxation time distributed data
Type belonging to the original lateral relaxation time distributed data；
Proportionality coefficient determining module is used for according to the corresponding proportionality coefficient computation model of the type to the original transverse direction
Relaxation time distributed data is handled, and the proportionality coefficient of constraint fluid and movable fluid is obtained；
Fluid distrbution determining module, for being utilized respectively the proportionality coefficient of the constraint fluid and movable fluid to the original
Beginning lateral relaxation time distributed data is handled, and the constraint fluid distrbution data and movable fluid point in the target work area are obtained
Cloth data.
On the other hand, this specification embodiment also provides a kind of constraint of reservoir and movable fluid is distributed and determines equipment, including
Processor and memory for storage processor executable instruction, realized when described instruction is executed by the processor include with
Lower step:
The nuclear magnetic resonance echo data for obtaining target work area carries out inverting to the echo data and obtains the target work area
Original lateral relaxation time distributed data；
It whether include that short relaxation peak determines the original transverse relaxation according to the original lateral relaxation time distributed data
Type belonging to Annual distribution data；
The original lateral relaxation time distributed data is carried out according to the corresponding proportionality coefficient computation model of the type
Processing obtains the proportionality coefficient of constraint fluid and movable fluid；
The proportionality coefficient of the constraint fluid and movable fluid is utilized respectively to the original lateral relaxation time distribution number
According to being handled, the constraint fluid distrbution data and movable fluid distributed data in the target work area are obtained.
On the other hand, this specification embodiment also provides a kind of constraint of reservoir and movable fluid is distributed the system of determination, described
System includes at least one processor and the memory for storing computer executable instructions, and the processor executes described instruction
The step of Shi Shixian above-mentioned any embodiment the method.
The reservoir constraint and movable fluid distribution determination method, device and be that this specification one or more embodiment provides
System, can be by previously according to lateral relaxation time T2Distribution characteristics it is different by T2Classify, and building is different respectively
The T of type2The corresponding data processing model of distribution.In actual use, original T is first determined2Classification belonging to distributed data, so
Afterwards, it is handled using the corresponding data processing model of the type, to determine the proportionality coefficient shape of constraint fluid and movable fluid
Formula.Recycle proportionality coefficient to original T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.It utilizes
Each embodiment of this specification can accurately realize the continuous and quantitative table that fluid and movable fluid distribution are fettered to reservoir
Sign.
Detailed description of the invention
In order to illustrate more clearly of this specification embodiment or technical solution in the prior art, below will to embodiment or
Attached drawing needed to be used in the description of the prior art is briefly described, it should be apparent that, the accompanying drawings in the following description is only
The some embodiments recorded in this specification, for those of ordinary skill in the art, in not making the creative labor property
Under the premise of, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is that the process of a kind of constraint of reservoir and movable fluid distribution determination method embodiment that this specification provides is illustrated
Figure；
Fig. 2 is the constraint fluid predicted conventional sandstone sample in one embodiment that this specification provides
Proportionality coefficient and movable fluid proportionality coefficient distribution schematic diagram；
Fig. 3 is the constraint stream predicted tight sand sample in another embodiment that this specification provides
Body proportionality coefficient and movable fluid proportionality coefficient distribution schematic diagram；
Fig. 4 is the constraint fluid T of the conventional sandstone sample in another embodiment that this specification provides2Distribution signal
Figure；
Fig. 5 is the constraint fluid T of the tight sand sample in another embodiment that this specification provides2Distribution signal
Figure；
Fig. 6 is the movable fluid T of the conventional sandstone sample in another embodiment that this specification provides2Distribution signal
Figure；
Fig. 7 is the movable fluid T of the tight sand sample in another embodiment that this specification provides2Distribution signal
Figure；
Fig. 8 is the modular structure that a kind of reservoir constraint that this specification provides and movable fluid are distributed determining device embodiment
Schematic diagram.
Specific embodiment
In order to make those skilled in the art more fully understand the technical solution in this specification, below in conjunction with this explanation
Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete
Site preparation description, it is clear that described embodiment is only specification a part of the embodiment, instead of all the embodiments.Based on saying
Bright book one or more embodiment, it is obtained by those of ordinary skill in the art without making creative efforts all
The range of this specification example scheme protection all should belong in other embodiments.
NMR centrefuge experiment can only measure limited core sample, cannot achieve to down-hole formation constraint and movable
The continuous characterization of fluid distrbution.Although NMR well logging can provide down-hole formation continuous T2Distribution, but due to fluid point in stratum
The complexity of cloth, there are no a kind of effective methods so far can be based on the T that NMR logs well2Distribution is accurate with realizing underground
The continuous and quantitative characterization of layer constraint and movable fluid distribution.
Correspondingly, this specification embodiment provides a kind of constraint of reservoir and movable fluid distribution determination method, Ke Yitong
It crosses previously according to lateral relaxation time T2Distribution characteristics it is different by T2Classify, and constructs different types of T respectively2
The corresponding data processing model of distribution.In actual use, original T is first determined2Then classification belonging to distributed data utilizes
The corresponding data processing model of the type is handled, and determines the proportionality coefficient form of constraint fluid and movable fluid, is recycled
Proportionality coefficient is to original T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.Utilize this specification
Each embodiment can accurately realize the continuous and quantitative characterization that fluid and movable fluid distribution are fettered to reservoir.
Fig. 1 is that a kind of reservoir constraint that this specification provides and movable fluid distribution determination method embodiment process are shown
It is intended to.Although being based on present description provides as the following examples or method operating procedure shown in the drawings or apparatus structure
Less behaviour after may include routinely or without creative labor more in the method or device or part merging
Make step or modular unit.In the step of there is no necessary causalities in logicality or structure, the execution of these steps is suitable
Sequence or the modular structure of device are not limited to this specification embodiment or execution shown in the drawings sequence or modular structure.The side
Device, server or the end product of method or modular structure in practice are in application, can be according to shown in embodiment or attached drawing
Method or modular structure carry out sequence execution or it is parallel execute (such as parallel processor or multiple threads environment,
It even include the implementation environment of distributed treatment, server cluster).
Specific one embodiment is as shown in Figure 1, the reservoir constraint and movable fluid that this specification provides are distributed determination side
In one embodiment of method, the method may include:
S102: obtaining the nuclear magnetic resonance echo data in target work area, carries out inverting to the echo data and obtains the mesh
Mark the original lateral relaxation time distributed data in work area.
The nuclear magnetic resonance log echo data in target work area can be acquired, then, the echo data is carried out at inverting
Reason obtains the original lateral relaxation time T in target work area2Distributed data.It, such as can be using singular value point in some embodiments
(Singular Value Decomposition, SVD), BRD algorithm (Butler-Reeds-Dawson, BRD) etc. are solved to described
Echo data carries out inverting, obtains original lateral relaxation time T2Distributed data.
Whether S104: including that short relaxation peak determines the original transverse direction according to the original lateral relaxation time distributed data
Type belonging to relaxation time distributed data.
It can be according to whether determining type belonging to the original lateral relaxation time distributed data comprising short relaxation peak.It is short
The division limits of relaxation can sets itself according to actual needs.It, can be with T in some embodiments of this specification2Value is
30ms is divided, and is short relaxation data less than or equal to 30ms, is long relaxation data greater than 30ms.
It, can be by the original lateral relaxation time T in some embodiments2Distributed data is based on whether include short relaxation peak
It is divided into two classes, is respectively set as the first kind and the second class, wherein the first kind is the T comprising short relaxation peak2Distributed data, second
Class is the T not comprising short relaxation peak2Distributed data.
Correspondingly, if the original lateral relaxation time distributed data includes short relaxation peak, it is determined that the original transverse direction
Type belonging to relaxation time distributed data is the first kind；If the original lateral relaxation time distributed data does not include short relaxation
Peak, it is determined that type belonging to the original lateral relaxation time distributed data is the second class.
S106: according to the corresponding proportionality coefficient computation model of the type to the original lateral relaxation time distributed data
It is handled, obtains the proportionality coefficient of constraint fluid and movable fluid.
It can be previously according to different types of T2The data characteristics of distributed data constructs corresponding proportionality coefficient computation model,
It, can the corresponding proportionality coefficient computation model pair of the type according to belonging to original lateral relaxation time distributed data when practical application
The original lateral relaxation time distributed data is handled, and the proportionality coefficient of constraint fluid and movable fluid is obtained.According to not
The T of same type2The data characteristics of distributed data constructs computation model respectively, and the ratio of constraint fluid and movable fluid can be improved
The accuracy that example coefficient determines, and then improve T2It is distributed the accuracy of characterization result.
In one embodiment of this specification, it is described according to the corresponding proportionality coefficient computation model of the type to the original
Beginning lateral relaxation time distributed data is handled, and may include:
If type belonging to original lateral relaxation time distributed data is the first kind, following first ratios system can use
Number computation model handles the original lateral relaxation time distributed data:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤ε1 (1)
Wherein, αi=i Δ α, Δ α > 0, i=1,2 ..., N, αiIndicate i-th of α value, N indicates the maximum value of α
Number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid proportionality coefficient
Function shape, pksIndicate that short relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data, ε1It is
A kind of lateral relaxation time is distributed threshold value, Pb(T2,αi) it is αiCorresponding constraint fluid proportional coefficient；
Determine the corresponding α of minimum i for meeting above-mentioned first proportionality coefficient computation modeli, it is denoted as α opt1；
By α opt1Corresponding Pb(T2,αopt1) and Pm(T2,αopt1) it is identified as constraint fluid proportional coefficient and movable stream
Body proportionality coefficient, wherein Pm(T2,αopt1)=1-Pb(T2,αopt1)。
In some embodiments, first kind lateral relaxation time is distributed threshold epsilon1It is demarcated using core experiment.Institute
It states short relaxation peak-to-peak value and seeks function pksFunction can be sought for the peak value provided in Matlab software, it is of course also possible to use
Other kinds of peak value seeks function, here without limitation.
In another embodiment of this specification, it is described according to the corresponding proportionality coefficient computation model of the type to described
Original lateral relaxation time distributed data is handled, and may include:
If type belonging to original lateral relaxation time distributed data is the second class, following second ratios system can use
Number computation model handles the original lateral relaxation time distributed data:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤ε2 (2)
Wherein, αj=(N+1-j) Δ α, Δ α > 0, j=1,2 ..., N, αjIndicate j-th of α value, N indicates the maximum of α
Value number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid ratio
The function shape of example coefficient, pklIndicate that long relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data,
ε2Threshold value, P are distributed for the second class lateral relaxation timem(T2,αi) it is αiCorresponding movable fluid proportionality coefficient；
Determine the corresponding α of minimum j for meeting above-mentioned second proportionality coefficient computation modelj, it is denoted as αopt2；
By αopt2Corresponding Pb(T2,αopt2) and Pm(T2,αopt2) it is identified as constraint fluid proportional coefficient and movable stream
Body proportionality coefficient, wherein Pb(T2,αopt2)=1-Pm(T2,αopt2)。
In some embodiments, the second class lateral relaxation time is distributed threshold epsilon2Core experiment can also be used to be marked
It is fixed.The long relaxation peak-to-peak value seeks function pklFunction can be sought for the peak value provided in Matlab software, can also used
Other kinds of peak value seeks function, here without limitation.
Original T is determined using model provided by the above embodiment2The constraint fluid and movable fluid proportionality coefficient of distribution, can
To further increase constraint fluid and movable fluid T2It is distributed the accuracy of continuous and quantitative characterization.
In one or more embodiment of this specification, the proportionality coefficient of the constraint fluid and movable fluid be can wrap
Include the exponential function using natural constant the bottom of as.Preferably, the proportionality coefficient function of the constraint fluid and movable fluid can wrap
It includes:
Wherein, Pb(T2, α) and it is constraint fluid proportional coefficient distribution function, Pm(T2, α) and it is that movable fluid proportionality coefficient is distributed
Function, α is form factor, to control Pb(T2, α) and Pm(T2, α) function shape.The optimal value of α can be according to above-mentioned meter
Model (1) or (2) is calculated to determine, thus the more accurate coefficient function of certainty ratio really Pb(T2, α) and Pm(T2, α) value.
S108: the proportionality coefficient of the constraint fluid and movable fluid is utilized respectively to the original lateral relaxation time point
Cloth data are handled, and the constraint fluid distrbution data and movable fluid distributed data in the target work area are obtained.
The proportionality coefficient of constraint fluid and movable fluid that above-mentioned steps determine be can use to the original transverse relaxation
Annual distribution data are handled, and the constraint fluid distrbution data and movable fluid distributed data are obtained.Some embodiments
In, it can be in the following manner to original lateral relaxation time T2Distributed data is handled:
fb(T2)=Pb(T2,αopt)·f(T2) (5)
fm(T2)=Pm(T2,αopt)·f(T2) (6)
Wherein, fb(T2) it is constraint fluid T2Distributed data, fm(T2) it is movable fluid T2Distributed data, Pb(T2,αopt) be
Fetter fluid proportional coefficient, Pm(T2,αopt) it is movable fluid proportionality coefficient.Belonging to original lateral relaxation time distributed data
When type is the first kind, αoptValue is αopt1, when the type belonging to the original lateral relaxation time distributed data is the second class,
αoptValue is αopt2。
The scheme provided based on the above embodiment, this specification provide a kind of scheme using above-described embodiment also to characterize
Reservoir fetters the specific example of fluid and movable fluid distribution situation, as follows:
Certain research 19 pieces of area conventional sandstone sample and certain research 19 pieces of area tight sand sample are collected, NMR experiment is carried out and surveys
Amount obtains the T under rock sample original (saturation fluid state), constraint fluid state and movable fluid state2Distribution.It needs to illustrate
It is the original T of conventional sandstone sample used in this example2It is distributed while including I class and II class T2Distribution, tight sand sample used are former
Beginning T2Distribution only includes I class T2Distribution.
Then, conventional sandstone sample and tight sand sample are carried out using this specification scheme provided by the above embodiment
Fluid and movable fluid forecast of distribution are fettered, as follows:
Step 1: acquisition conventional sandstone sample and the corresponding nuclear magnetic resonance echo data of tight sand sample, and to number of echoes
Original T is obtained according to inverting is carried out2Distribution.
Step 2: by T2Distribution is divided into two classes according to whether comprising short relaxation peak: if original T2Distribution includes short relaxation peak, then
For I class T2Distribution；If original T2Distribution does not include short relaxation peak, then is II class T2Distribution.
Step 3: construction is using natural constant as the exponential function at bottom, to characterize constraint and movable fluid proportionality coefficient point
Cloth, construction constraint and movable fluid proportionality coefficient distribution function can be with are as follows:
Wherein, Pb(T2, α) and it is constraint fluid proportional coefficient distribution function, Pm(T2, α) and it is that movable fluid proportionality coefficient is distributed
Function, α is form factor, to control Pb(T2, α) and Pm(T2, α) function shape.
Step 4: according to the exponential function of construction and original T2Type belonging to distribution, respectively to constraint and movable fluid T2
Distribution is characterized.
S401: for I class T2Distribution seeks constraint and movable fluid proportionality coefficient distribution function according to following formula:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤ε1
Wherein, αi=i Δ α, Δ α > 0, i=1,2 ..., N.pksIndicate that short relaxation peak-to-peak value seeks function, f (T2) be
Original T2Distribution, ε1For I class T2It is distributed threshold value, can be demarcated by core experiment.Meet the corresponding α of minimum i of the formulaiFor
Optimal α, is denoted as αopt, corresponding Pb(T2,αopt) and Pm(T2,αopt) it is the original T2The constraint of distribution and movable fluid ratio
Example coefficient function.
S402: for II class T2Distribution seeks constraint and movable fluid proportionality coefficient distribution function according to following formula:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤ε2
Wherein, αj=(N+1-j) Δ α, Δ α > 0, j=1,2 ..., N.pklIndicate that long relaxation peak-to-peak value seeks function, ε2
For II class T2It is distributed threshold value, can be demarcated by core experiment.Meet the corresponding α of minimum j of the formulajFor optimal α, it is denoted as
αopt, corresponding Pb(T2,αopt) and Pm(T2,αopt) it is the original T2The constraint of distribution and movable fluid proportionality coefficient function.
Step 5: according to the constraint and movable fluid proportionality coefficient function, seeking constraint and movable stream according to following formula
Body T2Distribution:
fb(T2)=Pb(T2,αopt)·f(T2)
fm(T2)=Pm(T2,αopt)·f(T2)
Wherein, fb(T2) it is constraint fluid T2Distribution, fm(T2) it is movable fluid T2Distribution.
Fig. 2 indicates the constraint stream obtained using the scheme of above-described embodiment to conventional sandstone sample prediction (predicted)
Body proportionality coefficient is distributed (solid line) and movable fluid proportionality coefficient distribution (dotted line).Fig. 3 indicates to utilize the scheme of above-described embodiment
The constraint fluid proportional coefficient distribution (solid line) and the distribution of movable fluid proportionality coefficient predict tight sand sample are (empty
Line).Wherein, the abscissa in Fig. 2 and Fig. 3 is T2Value, ordinate are sample number into spectrum, and " C " indicates that conventional sandstone, " T " indicate to cause
Close sandstone.For conventional sandstone sample, ε used in this example1=0.050, ε2=0.035；For tight sand sample, this example
ε used1=0.001.
Fig. 4 indicates that the experiment measurement (exp) of conventional sandstone sample fetters fluid T2It is distributed (solid line) and according to above-mentioned implementation
The constraint fluid T of the method prediction of example2It is distributed (dotted line).Fig. 5 indicates the experiment measurement constraint fluid T of tight sand sample2Distribution
(solid line) and the constraint fluid T predicted according to the method for above-described embodiment2It is distributed (dotted line).
Fig. 6 indicates that the experiment of conventional sandstone sample measures movable fluid T2It is distributed (solid line) and the side according to above-described embodiment
The movable fluid T of method prediction2It is distributed (dotted line).Fig. 7 indicates that the experiment of tight sand sample measures movable fluid T2It is distributed (solid line)
With the movable fluid T predicted according to the method for above-described embodiment2It is distributed (dotted line).
Analyses and comparison Fig. 4 to Fig. 7 is it is found that the reservoir constraint fluid and movable fluid point provided according to this specification embodiment
The constraint fluid and movable fluid T that cloth characterizing method is predicted2Distribution is almost the same with experimental measurements, shows to utilize this
The scheme that specification embodiment provides, which can be realized accurately and effectively, determines reservoir constraint and the continuous of movable fluid distribution situation
Scale sign.
It should be noted that examples detailed above is only schematically illustrated by taking conventional sandstone and tight sand as an example, this explanation
The scheme of book embodiment is still applicable in other types reservoir, is not construed as limiting to this.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Specifically it is referred to
The description of aforementioned relevant treatment related embodiment, does not do repeat one by one herein.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
The reservoir constraint and movable fluid distribution determination method that this specification one or more embodiment provides, can pass through
Previously according to lateral relaxation time T2Distribution characteristics it is different by T2Classify, and constructs different types of T respectively2Point
Data processing model corresponding to cloth.In actual use, original T is first determined2Then classification belonging to distributed data utilizes this
The corresponding data processing model of type is handled, in the form of determining the proportionality coefficient of constraint fluid and movable fluid.It recycles
Proportionality coefficient is to original T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.Utilize this specification
Each embodiment can accurately realize the continuous and quantitative characterization that fluid and movable fluid distribution are fettered to reservoir.
Based on reservoir constraint described above and movable fluid distribution determination method, this specification one or more embodiment
A kind of constraint of reservoir is also provided and movable fluid is distributed determining device.The device may include having used this specification implementation
System, software (application), module, component, server of example the method etc. simultaneously combine the necessary device for implementing hardware.It is based on
Same innovation thinking, the device in one or more embodiments that this specification embodiment provides is as described in the following examples.
Since the implementation that device solves the problems, such as is similar to method, the implementation of the specific device of this specification embodiment can be joined
See the implementation of preceding method, overlaps will not be repeated.Used below, term " unit " or " module " may be implemented pre-
Determine the combination of the software and/or hardware of function.Although device described in following embodiment is preferably realized with software,
The realization of the combination of hardware or software and hardware is also that may and be contemplated.Specifically, provided Fig. 8 shows specification
A kind of modular structure schematic diagram of reservoir constraint and movable fluid distribution determining device embodiment, as shown in figure 8, described device can
To include:
Data acquisition module 202 can be used for obtaining the nuclear magnetic resonance echo data in target work area, to the echo data
Carry out the original lateral relaxation time distributed data that inverting obtains the target work area；
Whether determination type module 204 can be used for according to the original lateral relaxation time distributed data including short relaxation
Henan peak determines type belonging to the original lateral relaxation time distributed data；
Proportionality coefficient determining module 206 can be used for the corresponding proportionality coefficient computation model of the type to described original
Lateral relaxation time distributed data is handled, and the proportionality coefficient of constraint fluid and movable fluid is obtained；
Fluid distrbution determining module 208 can be used for being utilized respectively the proportionality coefficient of the constraint fluid and movable fluid
The original lateral relaxation time distributed data is handled, obtain the target work area constraint fluid distrbution data and can
Dynamic fluid distrbution data.
In another embodiment of this specification, the determination type module 204 may include:
First kind determination unit, if can be used for the original lateral relaxation time distributed data includes short relaxation peak,
Then determine that type belonging to the original lateral relaxation time distributed data is the first kind；
Second Type determination unit, if can be used for the original lateral relaxation time distributed data not comprising short relaxation
Peak, it is determined that type belonging to the original lateral relaxation time distributed data is the second class.
In another embodiment of this specification, the proportionality coefficient determining module 206 may include:
First processing units, if can be used for type belonging to original lateral relaxation time distributed data is the first kind,
The original lateral relaxation time distributed data is handled using following first proportionality coefficient computation models:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤ε1
Wherein, αi=i Δ α, Δ α > 0, i=1,2 ..., N, αiIndicate i-th of α value, N indicates the maximum value of α
Number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid proportionality coefficient
Function shape, pksIndicate that short relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data, ε1It is
A kind of lateral relaxation time is distributed threshold value, Pb(T2,αi) it is αiCorresponding constraint fluid proportional coefficient；
First best factors determination unit is determined for meeting the minimum i of above-mentioned first proportionality coefficient computation model
Corresponding αi, it is denoted as αopt1；
First proportionality coefficient determination unit, can be used for αopt1Corresponding Pb(T2,αopt1) and Pm(T2,αopt1) true respectively
It is set to constraint fluid proportional coefficient and movable fluid proportionality coefficient, wherein Pm(T2,αopt1)=1-Pb(T2,αopt1)。
In another embodiment of this specification, the proportionality coefficient determining module 206 may include:
The second processing unit, if can be used for type belonging to original lateral relaxation time distributed data is the second class,
The original lateral relaxation time distributed data is handled using following second proportionality coefficient computation models:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤ε2
Wherein, αj=(N+1-j) Δ α, Δ α > 0, j=1,2 ..., N, αjIndicate j-th of α value, N indicates the maximum of α
Value number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid ratio
The function shape of example coefficient, pklIndicate that long relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data,
ε2Threshold value, P are distributed for the second class lateral relaxation timem(T2,αi) it is αiCorresponding movable fluid proportionality coefficient；
Second best factors determination unit is determined for meeting the minimum j of above-mentioned second proportionality coefficient computation model
Corresponding αj, it is denoted as αopt2；
Second proportionality coefficient determination unit, can be used for αopt2Corresponding Pb(T2,αopt2) and Pm(T2,αopt2) true respectively
It is set to constraint fluid proportional coefficient and movable fluid proportionality coefficient, wherein Pb(T2,αopt2)=1-Pm(T2,αopt2)。
In another embodiment of this specification, the fluid distrbution determining module 208 be can be also used for according to following public affairs
The determination of formula progress distributed data:
fb(T2)=Pb(T2,αopt)·f(T2)
fm(T2)=Pm(T2,αopt)·f(T2)
Wherein, fb(T2) it is constraint fluid lateral relaxation time distributed data, Pb(T2,αopt) it is constraint fluid proportional system
Number, fm(T2) it is movable fluid lateral relaxation time distributed data, Pm(T2,αopt) it is movable fluid proportionality coefficient, f (T2) it is original
Beginning lateral relaxation time distributed data.
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
The reservoir constraint and movable fluid that this specification one or more embodiment provides are distributed determining device, can pass through
Previously according to lateral relaxation time T2Distribution characteristics it is different by T2Classify, and constructs different types of T respectively2Point
Data processing model corresponding to cloth.In actual use, original T is first determined2Then classification belonging to distributed data utilizes this
The corresponding data processing model of type is handled, in the form of determining the proportionality coefficient of constraint fluid and movable fluid.It recycles
Proportionality coefficient is to original T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.Utilize this specification
Each embodiment can accurately realize the continuous and quantitative characterization that fluid and movable fluid distribution are fettered to reservoir.
Method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program
It collects and records on a storage medium, the storage medium can be read and be executed with computer, realize this specification embodiment institute
The effect of description scheme.Therefore, this specification also provides a kind of constraint of reservoir and movable fluid is distributed and determines equipment, including processing
The memory of device and storage processor executable instruction, when described instruction is executed by the processor realize the following steps are included:
The nuclear magnetic resonance echo data for obtaining target work area carries out inverting to the echo data and obtains the target work area
Original lateral relaxation time distributed data；
It whether include that short relaxation peak determines the original transverse relaxation according to the original lateral relaxation time distributed data
Type belonging to Annual distribution data；
The original lateral relaxation time distributed data is carried out according to the corresponding proportionality coefficient computation model of the type
Processing obtains the proportionality coefficient of constraint fluid and movable fluid；
The proportionality coefficient of the constraint fluid and movable fluid is utilized respectively to the original lateral relaxation time distribution number
According to being handled, the constraint fluid distrbution data and movable fluid distributed data in the target work area are obtained.
It should be noted that equipment described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
The storage medium may include the physical unit for storing information, usually by after information digitalization again with benefit
The media of the modes such as electricity consumption, magnetic or optics are stored.It may include: that letter is stored in the way of electric energy that the storage medium, which has,
The device of breath such as, various memory, such as RAM, ROM；The device of information is stored in the way of magnetic energy such as, hard disk, floppy disk, magnetic
Band, core memory, magnetic bubble memory, USB flash disk；Using optical mode storage information device such as, CD or DVD.Certainly, there are also it
Readable storage medium storing program for executing of his mode, such as quantum memory, graphene memory etc..
The constraint of reservoir described in above-described embodiment and movable fluid, which are distributed, determines equipment, can be by previously according to lateral relaxation
Henan time T2Distribution characteristics it is different by T2Classify, and constructs different types of T respectively2The corresponding data of distribution
Handle model.In actual use, original T is first determined2Then classification belonging to distributed data utilizes the corresponding data of the type
Processing model is handled, in the form of determining the proportionality coefficient of constraint fluid and movable fluid.Recycle proportionality coefficient to original
T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.It, can using each embodiment of this specification
Accurately to realize the continuous and quantitative characterization for fettering fluid and movable fluid distribution to reservoir.
This specification also provides a kind of constraint of reservoir and movable fluid is distributed the system of determination, and the system can be individual
Fluid distrbution feature determines system, can also apply in a variety of Reservoir Analysis systems.The system can be individually clothes
Business device also may include the service for having used one or more the methods or one or more embodiment devices of this specification
Device cluster, system (including distributed system), software (application), practical operation device, logic gates device, quantum computer
Deng and combine it is necessary implement hardware terminal installation.Reservoir constraint and movable fluid be distributed the system of determination may include to
A few processor and the memory for storing computer executable instructions, the processor is realized above-mentioned when executing described instruction
The step of method described in any one or multiple embodiments.
It should be noted that system described above can also include others according to the description of method or Installation practice
Embodiment, concrete implementation mode are referred to the description of related method embodiment, do not repeat one by one herein.
The constraint of reservoir described in above-described embodiment and movable fluid are distributed the system of determination, can be by previously according to lateral relaxation
Henan time T2Distribution characteristics it is different by T2Classify, and constructs different types of T respectively2The corresponding data of distribution
Handle model.In actual use, original T is first determined2Then classification belonging to distributed data utilizes the corresponding data of the type
Processing model is handled, in the form of determining the proportionality coefficient of constraint fluid and movable fluid.Recycle proportionality coefficient to original
T2Distributed data is handled, and the T of constraint fluid and movable fluid is obtained2Distribution.It, can using each embodiment of this specification
Accurately to realize the continuous and quantitative characterization for fettering fluid and movable fluid distribution to reservoir.
It should be noted that this specification device or system described above according to the description of related method embodiment also
It may include other embodiments, concrete implementation mode is referred to the description of embodiment of the method, does not go to live in the household of one's in-laws on getting married one by one herein
It states.All the embodiments in this specification are described in a progressive manner, and same and similar part is mutual between each embodiment
Mutually referring to each embodiment focuses on the differences from other embodiments.Especially for hardware+program
For class, storage medium+program embodiment, since it is substantially similar to the method embodiment, so be described relatively simple, it is related
Place illustrates referring to the part of embodiment of the method.
Although mentioning SVD inverting, BRD inverting etc. in this specification embodiment content to obtain, definition, interaction, calculate, judgement
Deng operation and data description, still, this specification embodiment is not limited to comply with standard data model/template or sheet
Situation described in specification embodiment.Certain professional standards or the practice processes described using customized mode or embodiment
On embodiment modified slightly also may be implemented above-described embodiment it is identical, it is equivalent or it is close or deformation after it is anticipated that reality
Apply effect.Using the embodiment of the acquisitions such as these modifications or deformed data acquisition, storage, judgement, processing mode, still may be used
To belong within the scope of the optional embodiment of this specification.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used
Think personal computer, laptop computer, vehicle-mounted human-computer interaction device, cellular phone, camera phone, smart phone, individual
Digital assistants, media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or
The combination of any equipment in these equipment of person.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when specification one or more, it can also be with
The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Installation practice described above is only
It is only illustrative, for example, in addition the division of the unit, only a kind of logical function partition can have in actual implementation
Division mode, such as multiple units or components can be combined or can be integrated into another system or some features can be with
Ignore, or does not execute.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical
Some interfaces are crossed, the indirect coupling or communication connection of device or unit can be electrical property, mechanical or other forms.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method or equipment of element.
It will be understood by those skilled in the art that this specification one or more embodiment can provide as method, system or calculating
Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or
The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or
It is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to magnetic disk storage,
CD-ROM, optical memory etc.) on the form of computer program product implemented.
This specification one or more embodiment can computer executable instructions it is general on
It hereinafter describes, such as program module.Generally, program module includes executing particular task or realization particular abstract data type
Routine, programs, objects, component, data structure etc..This this specification one can also be practiced in a distributed computing environment
Or multiple embodiments, in these distributed computing environments, by being held by the connected remote processing devices of communication network
Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage equipment
In storage medium.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ",
The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material
Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term
Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description
Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other,
Those skilled in the art can be by different embodiments or examples described in this specification and different embodiments or examples
Feature is combined.
The foregoing is merely the embodiments of this specification, are not limited to this specification.For art technology
For personnel, this specification can have various modifications and variations.It is all made any within the spirit and principle of this specification
Modification, equivalent replacement, improvement etc., should be included within the scope of the claims of this specification.
Claims (10)
1. a kind of reservoir constraint and movable fluid distribution determination method characterized by comprising
The nuclear magnetic resonance echo data for obtaining target work area carries out the original that inverting obtains the target work area to the echo data
Beginning lateral relaxation time distributed data；
It whether include that short relaxation peak determines the original lateral relaxation time according to the original lateral relaxation time distributed data
Type belonging to distributed data；
The original lateral relaxation time distributed data is handled according to the corresponding proportionality coefficient computation model of the type,
Obtain the proportionality coefficient of constraint fluid and movable fluid；
Be utilized respectively it is described constraint fluid and movable fluid proportionality coefficient to the original lateral relaxation time distributed data into
Row processing, obtains the constraint fluid distrbution data and movable fluid distributed data in the target work area.
2. the method according to claim 1, wherein the type includes the first kind, the second class, correspondingly, institute
It states and determines type belonging to the original lateral relaxation time distributed data, comprising:
If the original lateral relaxation time distributed data includes short relaxation peak, it is determined that the original lateral relaxation time distribution
Type belonging to data is the first kind；
If the original lateral relaxation time distributed data does not include short relaxation peak, it is determined that the original lateral relaxation time point
Type belonging to cloth data is the second class.
3. according to the method described in claim 2, it is characterized in that, described calculate mould according to the corresponding proportionality coefficient of the type
Type handles the original lateral relaxation time distributed data, comprising:
If type belonging to original lateral relaxation time distributed data is the first kind, mould is calculated using following first proportionality coefficients
Type handles the original lateral relaxation time distributed data:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤ε1
Wherein, αi=i Δ α, Δ α > 0, i=1,2 ..., N, αiIndicate i-th of α value, N indicates the maximum value number of α, Δ
α indicates the value interval of α, and α is form factor, to control the function of constraint fluid proportional coefficient and movable fluid proportionality coefficient
Shape, pksIndicate that short relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data, ε1For first kind cross
Threshold value, P are distributed to the relaxation timeb(T2,αi) it is αiCorresponding constraint fluid proportional coefficient；
Determine the corresponding α of minimum i for meeting above-mentioned first proportionality coefficient computation modeli, it is denoted as αopt1；
By αopt1Corresponding Pb(T2,αopt1) and Pm(T2,αopt1) it is identified as constraint fluid proportional coefficient and movable fluid ratio
Coefficient, wherein Pm(T2,αopt1)=1-Pb(T2,αopt1)。
4. according to the method described in claim 2, it is characterized in that, described calculate mould according to the corresponding proportionality coefficient of the type
Type handles the original lateral relaxation time distributed data, comprising:
If type belonging to original lateral relaxation time distributed data is the second class, mould is calculated using following second proportionality coefficients
Type handles the original lateral relaxation time distributed data:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤ε2
Wherein, αj=(N+1-j) Δ α, Δ α > 0, j=1,2 ..., N, αjIndicate j-th of α value, N indicates the maximum value of α
Number, Δ α indicate the value interval of α, and α is form factor, to control constraint fluid proportional coefficient and movable fluid ratio system
Several function shapes, pklIndicate that long relaxation peak-to-peak value seeks function, f (T2) it is original lateral relaxation time distributed data, ε2For
Second class lateral relaxation time is distributed threshold value, Pm(T2,αi) it is αiCorresponding movable fluid proportionality coefficient；
Determine the corresponding α of minimum j for meeting above-mentioned second proportionality coefficient computation modelj, it is denoted as αopt2；
By αopt2Corresponding Pb(T2,αopt2) and Pm(T2,αopt2) it is identified as constraint fluid proportional coefficient and movable fluid ratio
Coefficient, wherein Pb(T2,αopt2)=1-Pm(T2,αopt2)。
5. the method according to claim 3 or 4, which is characterized in that the proportionality coefficient of the constraint fluid and movable fluid
Including using natural constant as the exponential function at bottom.
6. according to the method described in claim 5, it is characterized in that, the proportionality coefficient packet of the constraint fluid and movable fluid
It includes:
Wherein, Pb(T2, α) and it is constraint fluid proportional coefficient distribution function, Pm(T2, α) and it is movable fluid proportionality coefficient distribution function,
α is form factor, to control Pb(T2, α) and Pm(T2, α) function shape.
7. the method according to claim 1, wherein described be utilized respectively the constraint fluid and movable fluid
Proportionality coefficient function handles corresponding sorted original lateral relaxation time distributed data, comprising:
fb(T2)=Pb(T2,αopt)·f(T2)
fm(T2)=Pm(T2,αopt)·f(T2)
Wherein, fb(T2) it is constraint fluid lateral relaxation time distributed data, Pb(T2,αopt) it is constraint fluid proportional coefficient, fm
(T2) it is movable fluid lateral relaxation time distributed data, Pm(T2,αopt) it is movable fluid proportionality coefficient, f (T2) it is original cross
To relaxation time distributed data.
8. a kind of constraint of reservoir and movable fluid are distributed determining device, which is characterized in that described device includes:
Data acquisition module carries out inverting to the echo data and obtains for obtaining the nuclear magnetic resonance echo data in target work area
Obtain the original lateral relaxation time distributed data in the target work area；
Determination type module, for whether including described in short relaxation peak determines according to the original lateral relaxation time distributed data
Type belonging to original lateral relaxation time distributed data；
Proportionality coefficient determining module is used for according to the corresponding proportionality coefficient computation model of the type to the original transverse relaxation
Annual distribution data are handled, and the proportionality coefficient of constraint fluid and movable fluid is obtained；
Fluid distrbution determining module, for being utilized respectively the proportionality coefficient of the constraint fluid and movable fluid to the original cross
It is handled to relaxation time distributed data, obtains the constraint fluid distrbution data and movable fluid distribution number in the target work area
According to.
9. a kind of reservoir constraint and movable fluid are distributed and determine equipment, which is characterized in that handled including processor and for storage
The memory of device executable instruction, when described instruction is executed by the processor realize the following steps are included:
The nuclear magnetic resonance echo data for obtaining target work area carries out the original that inverting obtains the target work area to the echo data
Beginning lateral relaxation time distributed data；
It whether include that short relaxation peak determines the original lateral relaxation time according to the original lateral relaxation time distributed data
Type belonging to distributed data；
The original lateral relaxation time distributed data is handled according to the corresponding proportionality coefficient computation model of the type,
Obtain the proportionality coefficient of constraint fluid and movable fluid；
Be utilized respectively it is described constraint fluid and movable fluid proportionality coefficient to the original lateral relaxation time distributed data into
Row processing, obtains the constraint fluid distrbution data and movable fluid distributed data in the target work area.
10. a kind of reservoir constraint and movable fluid are distributed the system of determination, which is characterized in that the system comprises at least one processing
Device and the memory for storing computer executable instructions, the processor are realized in claim 1-7 when executing described instruction
The step of any one the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572371.XA CN110244369B (en) | 2019-06-28 | 2019-06-28 | Reservoir constraint and movable fluid distribution determination method, device and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572371.XA CN110244369B (en) | 2019-06-28 | 2019-06-28 | Reservoir constraint and movable fluid distribution determination method, device and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110244369A true CN110244369A (en) | 2019-09-17 |
| CN110244369B CN110244369B (en) | 2020-10-13 |
Family
ID=67890077
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910572371.XA Active CN110244369B (en) | 2019-06-28 | 2019-06-28 | Reservoir constraint and movable fluid distribution determination method, device and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110244369B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114046145A (en) * | 2021-11-26 | 2022-02-15 | 中国石油大学(北京) | Reservoir fluid identification and saturation determination method and device |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1869733A (en) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | Method for determining nuclear magnetic resonance logging T2 spectral T2 end value |
| CN102608664A (en) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | Method and device for obtaining transverse relaxation time spectrum by depth-dimension nuclear magnetic resonance inversion |
| CN104330433A (en) * | 2014-10-28 | 2015-02-04 | 中国石油天然气股份有限公司 | Method and device for obtaining transverse relaxation time distribution of target storage layer |
| CN104932027A (en) * | 2015-05-06 | 2015-09-23 | 中国石油大学(北京) | Reservoir classification method based on nuclear magnetic resonance logging |
| CN104991280A (en) * | 2015-07-02 | 2015-10-21 | 中国石油天然气股份有限公司 | Perforated stratum selection method and apparatus based on microscopic pore throat structure index |
| CN106249306A (en) * | 2016-10-12 | 2016-12-21 | 贵州大学 | Shale pore structure detection method based on nuclear magnetic resonance, NMR |
| CN106324688A (en) * | 2016-09-23 | 2017-01-11 | 中国石油大学(北京) | Reservoir irreducible water saturation determining method and device |
| CN106772645A (en) * | 2016-12-15 | 2017-05-31 | 中国石油大学(北京) | Nuclear magnetic resonance data inversion method and device based on the constraint of general prior information |
| CN107843611A (en) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | Low permeability sandstone reservoir moveable gel nuclear magnetic resonance parameter characterizes new method |
| CN108956678A (en) * | 2018-06-11 | 2018-12-07 | 西南石油大学 | A kind of T based on nuclear magnetic resonance log2Compose sensitive parameter extracting method |
| CN109030311A (en) * | 2018-07-16 | 2018-12-18 | 西南石油大学 | Based on nuclear magnetic resonance T2Compose the pore structure classification and recognition methods of sensitive parameter |
| CN109594971A (en) * | 2018-12-21 | 2019-04-09 | 中国石油天然气集团有限公司 | Flow net model method based on nuclear magnetic resonance log enhanced diffustion gas-bearing formation recognition factor |
-
2019
- 2019-06-28 CN CN201910572371.XA patent/CN110244369B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1869733A (en) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | Method for determining nuclear magnetic resonance logging T2 spectral T2 end value |
| CN102608664A (en) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | Method and device for obtaining transverse relaxation time spectrum by depth-dimension nuclear magnetic resonance inversion |
| CN104330433A (en) * | 2014-10-28 | 2015-02-04 | 中国石油天然气股份有限公司 | Method and device for obtaining transverse relaxation time distribution of target storage layer |
| CN104932027A (en) * | 2015-05-06 | 2015-09-23 | 中国石油大学(北京) | Reservoir classification method based on nuclear magnetic resonance logging |
| CN104991280A (en) * | 2015-07-02 | 2015-10-21 | 中国石油天然气股份有限公司 | Perforated stratum selection method and apparatus based on microscopic pore throat structure index |
| CN107843611A (en) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | Low permeability sandstone reservoir moveable gel nuclear magnetic resonance parameter characterizes new method |
| CN106324688A (en) * | 2016-09-23 | 2017-01-11 | 中国石油大学(北京) | Reservoir irreducible water saturation determining method and device |
| CN106249306A (en) * | 2016-10-12 | 2016-12-21 | 贵州大学 | Shale pore structure detection method based on nuclear magnetic resonance, NMR |
| CN106772645A (en) * | 2016-12-15 | 2017-05-31 | 中国石油大学(北京) | Nuclear magnetic resonance data inversion method and device based on the constraint of general prior information |
| CN108956678A (en) * | 2018-06-11 | 2018-12-07 | 西南石油大学 | A kind of T based on nuclear magnetic resonance log2Compose sensitive parameter extracting method |
| CN109030311A (en) * | 2018-07-16 | 2018-12-18 | 西南石油大学 | Based on nuclear magnetic resonance T2Compose the pore structure classification and recognition methods of sensitive parameter |
| CN109594971A (en) * | 2018-12-21 | 2019-04-09 | 中国石油天然气集团有限公司 | Flow net model method based on nuclear magnetic resonance log enhanced diffustion gas-bearing formation recognition factor |
Non-Patent Citations (3)
| Title |
|---|
| 丁业娇 等: ""利用核磁共振原始回波数据确定可动水饱和度的方法"", 《2015中国地球科学联合学术年会论文集》 * |
| 代齐全 等: ""基于核磁共振新参数的致密砂岩储层孔隙结构特征"", 《石油学报》 * |
| 谢然红 等: ""Advanced fluid-typing methods for NMR logging"", 《PET. SCI》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114046145A (en) * | 2021-11-26 | 2022-02-15 | 中国石油大学(北京) | Reservoir fluid identification and saturation determination method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110244369B (en) | 2020-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106372402B (en) | The parallel method of fuzzy region convolutional neural networks under a kind of big data environment | |
| Wang et al. | Tectonic–geochemical exploration modeling for characterizing geo-anomalies in southeastern Yunnan district, China | |
| CN103590827B (en) | Based on the compact clastic rock natural gas well PRODUCTION FORECASTING METHODS of Reservoir Classification | |
| CN104380260B (en) | Utilize the reservoir simulation of scalable grid computing | |
| CN107229768B (en) | Slopereliability parameter acquiring method and device based on fuzzy classification technology | |
| CN105891089B (en) | Reservoir permeability determines method and device | |
| CN109736794A (en) | Fractured pore type tight sandstone reservoir feature determines method, apparatus and system | |
| CN104101902B (en) | Seismic properties clustering method and device | |
| CN107797154B (en) | A kind of electric imaging logging image scale method and device | |
| CN109636094A (en) | Reservoir in low-permeability oilfields kind identification method and device | |
| CN105205239A (en) | Method and device for modeling reservoir physical property parameter | |
| CN106355571A (en) | Determination method and device for dolostone reservoir quality | |
| CN109298464A (en) | Tight sandstone reservoir Diagenetic Facies Logging Identification Method and device | |
| CN110244369A (en) | Reservoir constraint and movable fluid distribution determination method, apparatus and system | |
| Bhattacharyya et al. | A novel methodology for fast reservoir simulation of single-phase gas reservoirs using machine learning | |
| Pan et al. | Study on reservoir characteristics and evaluation methods of altered igneous reservoirs in Songliao Basin, China | |
| CN107991710B (en) | A kind of reservoir pores distribution acquiring method and device | |
| CN106682347B (en) | Discontinuous Deformation Analysis Method based on rock mass discontinuity random strength | |
| Altomare et al. | Using clouds for smart city applications | |
| CN105224777B (en) | A kind of oil gas unit investigation and prospecting degree recognition methods | |
| Hox et al. | Multigroup and multilevel approaches to measurement equivalence | |
| CN109308381A (en) | A kind of method, apparatus and system of determining oil volume factor | |
| CN108761533B (en) | A kind of method, apparatus and system of determining P-S wave velocity ratio | |
| CN108897061A (en) | A kind of method, apparatus and system of determining sandstone reservoir ratio | |
| CN109100802A (en) | Reservoir determines method and apparatus |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |