CN110244369B - Reservoir constraint and movable fluid distribution determination method, device and system - Google Patents
Reservoir constraint and movable fluid distribution determination method, device and system Download PDFInfo
- Publication number
- CN110244369B CN110244369B CN201910572371.XA CN201910572371A CN110244369B CN 110244369 B CN110244369 B CN 110244369B CN 201910572371 A CN201910572371 A CN 201910572371A CN 110244369 B CN110244369 B CN 110244369B
- Authority
- CN
- China
- Prior art keywords
- fluid
- relaxation time
- distribution data
- time distribution
- transverse 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.)
- Active
Links
Images
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the specification discloses a reservoir constraint and movable fluid distribution determination method, a device and a system, wherein the method comprises the steps of obtaining nuclear magnetic resonance echo data of a target work area, and carrying out inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area; determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not; processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid; and processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
Description
Technical Field
The invention relates to the technical field of logging data processing in oil exploration and development, in particular to a method, a device and a system for determining reservoir constraint and movable fluid distribution.
Background
The distribution of the bound fluid and the mobile fluid can directly reflect the rock physical properties such as the micro-pore structure characteristics and the seepage capability of the rock. In the exploration and development of oil and gas resources, the distribution characterization of effective bound fluids and movable fluids, particularly the distribution characterization of the movable fluids, is an important basis for reservoir evaluation, capacity prediction and efficient reservoir development. Therefore, the method for researching reservoir bound fluid and mobile fluid distribution characterization is of great significance.
Transverse relaxation by low field Nuclear Magnetic Resonance (NMR) measurementsTime (T)2) The distribution may characterize the fluid distribution in the rock pores. Laboratory NMR experimental measurements on saturated fluid rocks and centrifuged rocks, respectively, are an effective means for NMR techniques to characterize the distribution of bound and mobile fluids. Firstly, carrying out NMR measurement on a rock core in a saturated fluid state to obtain a saturated fluid T2And (4) distribution. Centrifuging out movable fluid in the core under the condition of reservoir pressure, and performing NMR measurement on the core in the bound fluid state to obtain bound fluid T2And (4) distribution. Saturated fluid T2Distribution minus bound fluid T2Distributing to obtain movable fluid T2And (4) distribution. However, NMR centrifugation experiments can only measure limited core samples and cannot achieve continuous characterization of downhole formation restriction and mobile fluid distribution. Although NMR logging may provide a continuous T of the formation downhole2Distribution, but due to the complexity of fluid distribution in the formation, there has not been an effective method to date based on T obtained from NMR logging2The distribution accurately realizes the continuous quantitative characterization of the underground formation constraint and the movable fluid distribution.
Disclosure of Invention
The embodiments of the present specification aim to provide a reservoir bound and mobile fluid distribution determination method, device and system, which can accurately implement continuous characterization of reservoir bound fluid and mobile fluid distribution.
The specification provides a reservoir constraint and movable fluid distribution determination method, a device and a system, which are realized by the following modes:
a reservoir tiedown and mobile fluid distribution determination method, comprising:
acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not;
processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid;
and processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
In another embodiment of the method provided in this specification, the types include a first type and a second type, and accordingly, the determining the type to which the original transverse relaxation time distribution data belongs includes:
if the original transverse relaxation time distribution data contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a first type;
and if the original transverse relaxation time distribution data do not contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a second type.
In another embodiment of the method provided in this specification, the processing the original transverse relaxation time distribution data according to the scale factor calculation model corresponding to the type includes:
if the type of the original transverse relaxation time distribution data is a first type, processing the original transverse relaxation time distribution data by using a first proportional coefficient calculation model as follows:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiA corresponding bound fluid proportionality coefficient;
determining α corresponding to the minimum i that satisfies the first scale factor calculation modeliIs recorded as αopt1；
α will be mixedopt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1)。
In another embodiment of the method provided in this specification, the processing the original transverse relaxation time distribution data according to the scale factor calculation model corresponding to the type includes:
if the type of the original transverse relaxation time distribution data is a second type, processing the original transverse relaxation time distribution data by using a second proportionality coefficient calculation model as follows:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklRepresents the long relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiA corresponding mobile fluid proportionality coefficient;
determining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2；
α will be mixedopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)。
In another embodiment of the method provided herein, the proportionality coefficients for the bound and mobile fluids comprise an exponential function based on a natural constant.
In another embodiment of the method provided herein, the proportionality coefficients of the confining fluid and the mobile fluid comprise:
wherein, Pb(T2α) is a bound fluid proportionality coefficient distribution function, Pm(T2α) is a movable fluid proportionality coefficient distribution function and α is a form factor for controlling Pb(T2α) and Pm(T2α).
In another embodiment of the method provided in this specification, the processing the corresponding classified raw transverse relaxation time distribution data by using the proportionality coefficient functions of the bound fluid and the mobile fluid respectively includes:
fb(T2)＝Pb(T2,αopt)·f(T2)
fm(T2)＝Pm(T2,αopt)·f(T2)
wherein f isb(T2) For confining the transverse relaxation time distribution data of the fluid, Pb(T2,αopt) To bound the fluid proportionality coefficient, fm(T2) For transverse relaxation time distribution data of the mobile fluid, Pm(T2,αopt) Is a mobile fluid proportionality coefficient, f (T)2) Raw transverse relaxation time distribution data.
In another aspect, embodiments of the present description also provide a reservoir tiedown and mobile fluid distribution determination apparatus, the apparatus comprising:
the data acquisition module is used for acquiring nuclear magnetic resonance echo data of a target work area and inverting the echo data to obtain original transverse relaxation time distribution data of the target work area;
the type determining module is used for determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contains a short relaxation peak or not;
a proportionality coefficient determining module, configured to process the original transverse relaxation time distribution data according to the proportionality coefficient calculation model corresponding to the type, so as to obtain proportionality coefficients of a bound fluid and a movable fluid;
and the fluid distribution determining module is used for processing the original transverse relaxation time distribution data by utilizing the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain the bound fluid distribution data and the movable fluid distribution data of the target work area.
In another aspect, embodiments of the present specification further provide a reservoir tiedown and mobile fluid distribution determination apparatus, including a processor and a memory for storing processor-executable instructions that when executed by the processor implement steps including:
acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not;
processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid;
and processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
In another aspect, embodiments of the present description also provide a reservoir tiedown and mobile fluid distribution determination system, the system comprising at least one processor and a memory storing computer-executable instructions that, when executed by the processor, implement the steps of the method of any of the embodiments described above.
The method, the device and the system for determining reservoir constraint and movable fluid distribution provided by one or more embodiments of the specification can be used for determining the reservoir constraint and movable fluid distribution according to the transverse relaxation time T in advance2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2The classification to which the distribution data belongs is then processed using a data processing model corresponding to that type to determine the form of the scaling coefficients for the bound and mobile fluids. Reuse of scaling factor vs. original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a schematic flow chart diagram of an embodiment of a reservoir tiedown and mobile fluid distribution determination method provided herein;
figure 2 is a schematic illustration of a binding fluid proportionality coefficient and mobile fluid proportionality coefficient distribution predicted for a conventional sandstone sample in one embodiment provided herein;
fig. 3 is a schematic diagram of the distribution of the scale factor of the confining fluid and the scale factor of the mobile fluid predicted from the tight sandstone sample in another embodiment provided in the present specification;
figure 4 is a tethering fluid T of a conventional sandstone sample in another embodiment provided herein2A schematic distribution diagram;
figure 5 is a confining fluid T of a tight sandstone sample in another embodiment provided in the present specification2A schematic distribution diagram;
figure 6 is a mobile fluid T of a conventional sandstone sample in another embodiment provided herein2A schematic distribution diagram;
figure 7 is a mobile fluid T of a tight sandstone sample in another embodiment provided herein2A schematic distribution diagram;
fig. 8 is a schematic block diagram of an embodiment of a reservoir tiedown and mobile fluid distribution determination apparatus provided herein.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in one or more embodiments of the present specification will be clearly and completely described below with reference to the drawings in one or more embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments of the specification, and not all embodiments. All other embodiments obtained by a person skilled in the art based on one or more embodiments of the present specification without making any creative effort shall fall within the protection scope of the embodiments of the present specification.
NMR centrifugal experiments can only measure limited core samples, and continuous characterization of underground formation constraint and mobile fluid distribution cannot be realized. Although NMR logging may provide a continuous T of the formation downhole2Distribution, but due to the complexity of fluid distribution in the formation, there has not been an effective method to date based on T obtained from NMR logging2The distribution accurately realizes the continuous quantitative characterization of the underground formation constraint and the movable fluid distribution.
Accordingly, the embodiments of the present specification provide a reservoir restriction and mobile flowThe volume distribution determination method can be carried out in advance according to the transverse relaxation time T2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2Classifying the distributed data, processing by the data processing model corresponding to the type, determining the proportional coefficient form of the bound fluid and the movable fluid, and using the proportional coefficient to process the original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
Fig. 1 is a schematic flow chart of an embodiment of the reservoir tiedown and mobile fluid distribution determination method provided in the present specification. Although the present specification provides the method steps or apparatus structures as shown in the following examples or figures, more or less steps or modules may be included in the method or apparatus structures based on conventional or non-inventive efforts. In the case of steps or structures which do not logically have the necessary cause and effect relationship, the execution order of the steps or the block structure of the apparatus is not limited to the execution order or the block structure shown in the embodiments or the drawings of the present specification. The described method or module structure can be executed in sequence or in parallel according to the embodiments or the method or module structure shown in the drawings (for example, in the environment of parallel processors or multi-thread processing, or even in the environment of distributed processing and server cluster) when the method or module structure is applied to a device, a server or an end product in practice.
In one embodiment of the reservoir tiedown and mobile fluid distribution determination method provided herein, as shown in fig. 1, the method may comprise:
s102: acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area.
Nuclear magnetic resonance logging echo data of target work area can be collectedThen, the echo data is inverted to obtain the original transverse relaxation time T of the target work area2And distributing the data. In some embodiments, for example, Singular Value Decomposition (SVD), BRD (Butler-streams-Dawson, BRD) and the like may be used to perform inversion on the echo data to obtain the original transverse relaxation time T2And distributing the data.
S104: and determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contains short relaxation peaks.
The type to which the raw transverse relaxation time distribution data belongs can be determined from whether short relaxation peaks are included. The division limit of the short relaxation can be set according to the actual requirement. In some embodiments of the present description, T may be2The values are divided into 30ms, short relaxation data if the values are less than or equal to 30ms, and long relaxation data if the values are greater than 30 ms.
In some embodiments, the original transverse relaxation time T may be determined2The distribution data is divided into a first class and a second class based on whether the distribution data contains a short relaxation peak, wherein the first class is T containing a short relaxation peak2Distribution data, second type T without short relaxation peaks2And distributing the data.
Correspondingly, if the original transverse relaxation time distribution data contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a first type; and if the original transverse relaxation time distribution data do not contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a second type.
S106: and processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid.
Can be preset according to different types of T2The data characteristics of the distribution data construct a corresponding proportional coefficient calculation model, and when the proportional coefficient calculation model is actually applied, the proportional coefficient calculation model corresponding to the type of the original transverse relaxation time distribution data can be usedAnd processing the original transverse relaxation time distribution data to obtain the proportionality coefficients of the bound fluid and the movable fluid. According to different types of T2The data characteristics of the distributed data are used for respectively constructing calculation models, so that the accuracy of determining the proportionality coefficients of the bound fluid and the movable fluid can be improved, and the T is further improved2The distribution characterizes the accuracy of the results.
In an embodiment of the present specification, the processing the original transverse relaxation time distribution data according to the scale factor calculation model corresponding to the type may include:
if the type of the original transverse relaxation time distribution data is a first type, the original transverse relaxation time distribution data can be processed by using the following first scale factor calculation model:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1(1)
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiA corresponding bound fluid proportionality coefficient;
determining α corresponding to the minimum i that satisfies the first scale factor calculation modeliRecorded as α opt1；
α opt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1)。
In some embodiments, the first type of lateral relaxationThreshold of distribution of the time of relaxation1The calibration can be performed by using a core experiment. The short relaxation peak-to-peak value is used for solving a function pksThe peak extracting function provided in Matlab software may be used, and of course, other types of peak extracting functions may be used, which are not limited herein.
In another embodiment of the present specification, the processing the original transverse relaxation time distribution data according to the scale factor calculation model corresponding to the type may include:
if the type of the original transverse relaxation time distribution data is the second type, the original transverse relaxation time distribution data can be processed by using the following second proportionality coefficient calculation model:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2(2)
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklRepresents the long relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiA corresponding mobile fluid proportionality coefficient;
determining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2；
α will be mixedopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)。
In some embodiments, the second type of transverse relaxation time distribution threshold2Can also be marked by using a rock core experimentAnd (4) determining. The long relaxation peak-to-peak value is used for solving a function pklThe peak extracting function provided in Matlab software may be used, and other types of peak extracting functions may also be used, which are not limited herein.
Determining the original T Using the model provided in the above embodiments2The proportion coefficient of the distributed bound fluid and the movable fluid can further improve the bound fluid and the movable fluid T2Accuracy of continuous quantitative characterization of distribution.
In one or more embodiments of the present disclosure, the proportionality coefficients for the bound and mobile fluids may include an exponential function based on a natural constant. Preferably, the scaling factor function for the bound fluid and the mobile fluid may include:
wherein, Pb(T2α) is a bound fluid proportionality coefficient distribution function, Pm(T2α) is a movable fluid proportionality coefficient distribution function and α is a form factor for controlling Pb(T2α) and Pm(T2α) the optimal value of α can be determined from the above-described calculation model (1) or (2) to more accurately determine the scale factor function Pb(T2α) and Pm(T2α).
S108: and processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
The raw transverse relaxation time distribution data may be processed using the proportionality coefficients of the bound fluid and the mobile fluid determined in the above steps to obtain the bound fluid distribution data and the mobile fluid distribution data. In some embodiments, the following can be said to applyOriginal transverse relaxation time T2And (3) processing the distribution data:
fb(T2)＝Pb(T2,αopt)·f(T2) (5)
fm(T2)＝Pm(T2,αopt)·f(T2) (6)
wherein f isb(T2) To confine fluid T2Distribution data, fm(T2) Is a mobile fluid T2Distribution data, Pb(T2,αopt) To constrain the fluid proportionality coefficient, Pm(T2,αopt) When the type of the original transverse relaxation time distribution data is of the first type, αoptValue of αopt1α when the type to which the original transverse relaxation time distribution data belongs is of the second classoptValue of αopt2。
Based on the solutions provided by the above embodiments, the present specification further provides a specific example of applying the solutions of the above embodiments to characterize reservoir-bound fluid and mobile fluid distribution, as follows:
collecting 19 conventional sandstone samples in a certain research area and 19 compact sandstone samples in the certain research area, and carrying out NMR experimental measurement to obtain original rock samples (in a saturated fluid state), a bound fluid state and T in a movable fluid state2And (4) distribution. Note that the conventional sandstone sample used in this example had the original T2The distribution includes both class I and class II T2Distribution, the original T of the tight sandstone sample used2The distribution containing only class I T2And (4) distribution.
Then, the scheme provided by the above embodiment of the present specification is used to predict the distribution of the confining fluid and mobile fluid of the conventional sandstone sample and the tight sandstone sample, as follows:
step 1: collecting nuclear magnetic resonance echo data corresponding to a conventional sandstone sample and a compact sandstone sample, and performing inversion on the echo data to obtain an original T2And (4) distribution.
Step 2: will T2According to whether the package is distributedThe short relaxation-containing peaks are divided into two categories: if original T2The distribution contains a short relaxation peak, then is class I T2Distributing; if original T2The distribution does not contain a short relaxation peak, and is a type II T2And (4) distribution.
And step 3: constructing an exponential function with a natural constant as a base to characterize the distribution of the confinement and mobile fluid proportionality coefficients, wherein the confinement and mobile fluid proportionality coefficient distribution function can be constructed by:
wherein, Pb(T2α) is a bound fluid proportionality coefficient distribution function, Pm(T2α) is a movable fluid proportionality coefficient distribution function and α is a form factor for controlling Pb(T2α) and Pm(T2α).
And 4, step 4: according to the constructed exponential function and the original T2The type of distribution being, respectively, bound and mobile fluid T2And (5) characterizing the distribution.
S401: for class I T2And (3) calculating a distribution function of the proportion coefficients of the bound fluid and the movable fluid according to the following formula:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N。pksRepresents the short relaxation peak-to-peak solving function, f (T)2) Is original T2The distribution of the water content is carried out,1is a class I T2Distribution threshold, which can be calibrated by core experiments α corresponding to the minimum i satisfying this formulaiIs optimum α and is recorded as αoptCorresponding to Pb(T2,αopt) And Pm(T2,αopt) Is the original T2Distributed bundleBound and mobile fluid scaling factor functions.
S402: for class II T2And (3) calculating a distribution function of the proportion coefficients of the bound fluid and the movable fluid according to the following formula:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N。pklRepresenting a long relaxation peak-to-peak solving function,2is class II T2Distribution threshold, which can be calibrated by core experiments α corresponding to the minimum j satisfying this formulajIs optimum α and is recorded as αoptCorresponding to Pb(T2,αopt) And Pm(T2,αopt) Is the original T2Distribution constraints and mobile fluid proportionality coefficient functions.
And 5: calculating the constraint and mobile fluid T according to the following formula according to the constraint and mobile fluid proportionality coefficient function2Distribution:
fb(T2)＝Pb(T2,αopt)·f(T2)
fm(T2)＝Pm(T2,αopt)·f(T2)
wherein f isb(T2) To confine fluid T2Distribution, fm(T2) Is a mobile fluid T2And (4) distribution.
Figure 2 shows the bounding fluid proportionality coefficient distribution (solid line) and the mobile fluid proportionality coefficient distribution (dashed line) predicted (predicted) for a conventional sandstone sample using the scheme of the above example. Figure 3 shows the binding fluid proportionality coefficient distribution (solid line) and mobile fluid proportionality coefficient distribution (dashed line) predicted for a tight sand sample using the protocol of the example described above. Wherein the abscissa in FIGS. 2 and 3 is T2The values, ordinate, are sample numbers, "C" for regular sandstone and "T" for tight sandstone. For conventional sandstone samples, the example used1＝0.050，2＝0.035; for tight sandstone samples, the example uses1＝0.001。
Figure 4 shows the experimental measurement (exp) of a conventional sandstone sample bound fluid T2Distribution (solid line) and bound fluid T predicted according to the method of the above embodiment2Distribution (dashed line). Figure 5 shows an experimental measurement of a tight sandstone sample tie-up fluid T2Distribution (solid line) and bound fluid T predicted according to the method of the above embodiment2Distribution (dashed line).
Figure 6 shows the experimentally measured mobile fluid T of a conventional sandstone sample2Distribution (solid line) and mobile fluid T predicted according to the method of the above embodiment2Distribution (dashed line). Figure 7 shows experimentally measured mobile fluid T of tight sandstone samples2Distribution (solid line) and mobile fluid T predicted according to the method of the above embodiment2Distribution (dashed line).
As can be seen from the analysis and comparison of fig. 4 to fig. 7, the reservoir bound fluid and mobile fluid distribution characterization method provided in the embodiments of the present specification predicts the bound fluid and mobile fluid T2The distribution is basically consistent with the experimental measurement result, which shows that the scheme provided by the embodiment of the specification can accurately and effectively realize continuous quantitative characterization on the reservoir constraint and movable fluid distribution conditions.
It should be noted that the above examples are only schematically illustrated by taking conventional sandstone and tight sandstone as examples, and the solution of the embodiment of the present disclosure is still applicable to other types of reservoirs, and is not limited thereto.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. For details, reference may be made to the description of the related embodiments of the related processing, and details are not repeated herein.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
The reservoir constraint and mobile fluid distribution determination method provided by one or more embodiments of the specification can be determined according to the transverse relaxation time T in advance2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2The classification to which the distribution data belongs is then processed using a data processing model corresponding to that type to determine the form of the scaling coefficients for the bound and mobile fluids. Reuse of scaling factor vs. original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
Based on the reservoir constraint and mobile fluid distribution determination method, one or more embodiments of the present specification further provide a reservoir constraint and mobile fluid distribution determination apparatus. The apparatus may include systems, software (applications), modules, components, servers, etc. that utilize the methods described in the embodiments of the present specification in conjunction with hardware implementations as necessary. Based on the same innovative conception, embodiments of the present specification provide an apparatus as described in the following embodiments. Since the implementation scheme of the apparatus for solving the problem is similar to that of the method, the specific implementation of the apparatus in the embodiment of the present specification may refer to the implementation of the foregoing method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated. Specifically, fig. 8 is a schematic block diagram illustrating an embodiment of a reservoir-bound and mobile-fluid-distribution determining apparatus provided in the specification, where, as shown in fig. 8, the apparatus may include:
the data acquisition module 202 may be configured to acquire nuclear magnetic resonance echo data of a target work area, and perform inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
a type determining module 204, configured to determine a type to which the original transverse relaxation time distribution data belongs according to whether the original transverse relaxation time distribution data contains a short relaxation peak;
a proportionality coefficient determining module 206, configured to process the original transverse relaxation time distribution data by using the proportionality coefficient calculation model corresponding to the type to obtain proportionality coefficients of a bound fluid and a movable fluid;
the fluid distribution determining module 208 may be configured to process the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the mobile fluid, respectively, to obtain bound fluid distribution data and mobile fluid distribution data of the target work area.
In another embodiment of the present specification, the type determining module 204 may include:
a first type determination unit, configured to determine that the type of the original transverse relaxation time distribution data is a first type if the original transverse relaxation time distribution data includes a short relaxation peak;
the second type determining unit may be configured to determine that the type to which the original transverse relaxation time distribution data belongs is of a second type if the original transverse relaxation time distribution data does not include a short relaxation peak.
In another embodiment of the present disclosure, the scaling factor determining module 206 may include:
the first processing unit may be configured to, if the type to which the original transverse relaxation time distribution data belongs is a first type, process the original transverse relaxation time distribution data by using the following first scale factor calculation model:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiA corresponding bound fluid proportionality coefficient;
a first optimal factor determination unit for determining α corresponding to the minimum i satisfying the first scale factor calculation modeliIs recorded as αopt1；
A first scaling factor determination unit, which may be used to determine αopt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1)。
In another embodiment of the present disclosure, the scaling factor determining module 206 may include:
the second processing unit may be configured to, if the type of the original transverse relaxation time distribution data is a second type, process the original transverse relaxation time distribution data by using the following second scale factor calculation model:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklTo representLong relaxation peak-to-peak function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiA corresponding mobile fluid proportionality coefficient;
a second optimization factor determination unit for determining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2；
A second scaling factor determination unit, which may be used to determine αopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)。
In another embodiment of the present description, the fluid distribution determination module 208 may be further configured to determine the distribution data according to the following formula:
fb(T2)＝Pb(T2,αopt)·f(T2)
fm(T2)＝Pm(T2,αopt)·f(T2)
wherein f isb(T2) For confining the transverse relaxation time distribution data of the fluid, Pb(T2,αopt) To bound the fluid proportionality coefficient, fm(T2) For transverse relaxation time distribution data of the mobile fluid, Pm(T2,αopt) Is a mobile fluid proportionality coefficient, f (T)2) Raw transverse relaxation time distribution data.
It should be noted that the above-described apparatus may also include other embodiments according to the description of the method embodiment. The specific implementation manner may refer to the description of the related method embodiment, and is not described in detail herein.
The reservoir confinement and mobile fluid distribution determination device provided in one or more embodiments of the present specification can be determined by previously determining the reservoir confinement and mobile fluid distribution according to the lateral directionRelaxation time T2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2The classification to which the distribution data belongs is then processed using a data processing model corresponding to that type to determine the form of the scaling coefficients for the bound and mobile fluids. Reuse of scaling factor vs. original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
The method or apparatus provided by the present specification and described in the foregoing embodiments may implement service logic through a computer program and record the service logic on a storage medium, where the storage medium may be read and executed by a computer, so as to implement the effect of the solution described in the embodiments of the present specification. Accordingly, the present specification also provides a reservoir tiedown and mobile fluid distribution determination apparatus comprising a processor and a memory storing processor-executable instructions which when executed by the processor implement steps comprising:
acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not;
processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid;
and processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
It should be noted that the above description of the apparatus according to the method embodiment may also include other embodiments. The specific implementation manner may refer to the description of the related method embodiment, and is not described in detail herein.
The storage medium may include a physical device for storing information, and typically, the information is digitized and then stored using an electrical, magnetic, or optical media. The storage medium may include: devices that store information using electrical energy, such as various types of memory, e.g., RAM, ROM, etc.; devices that store information using magnetic energy, such as hard disks, floppy disks, tapes, core memories, bubble memories, and usb disks; devices that store information optically, such as CDs or DVDs. Of course, there are other ways of storing media that can be read, such as quantum memory, graphene memory, and so forth.
The reservoir confinement and mobile fluid distribution determination apparatus according to the above embodiment may be determined by pre-determining the transverse relaxation time T2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2The classification to which the distribution data belongs is then processed using a data processing model corresponding to that type to determine the form of the scaling coefficients for the bound and mobile fluids. Reuse of scaling factor vs. original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
The present specification also provides a reservoir tiedown and mobile fluid distribution determination system that may be a single fluid distribution characterization system or may be applied in a variety of reservoir analysis systems. The system may be a single server, or may include a server cluster, a system (including a distributed system), software (applications), an actual operating device, a logic gate device, a quantum computer, etc. using one or more of the methods or one or more of the example devices of the present specification, in combination with a terminal device implementing hardware as necessary. The reservoir tiedown and mobile fluid distribution determination system may comprise at least one processor and memory storing computer executable instructions which when executed by the processor implement the steps of the method described in any one or more of the embodiments above.
It should be noted that the above-mentioned system may also include other implementation manners according to the description of the method or apparatus embodiment, and specific implementation manners may refer to the description of the related method embodiment, which is not described in detail herein.
The reservoir confinement and mobile fluid distribution determination system of the above embodiment may be determined by pre-determining the transverse relaxation time T2Difference in distribution characteristics will T2The distributions are classified and different types of T are constructed respectively2The corresponding data processing model is distributed. When in actual use, the original T is determined2The classification to which the distribution data belongs is then processed using a data processing model corresponding to that type to determine the form of the scaling coefficients for the bound and mobile fluids. Reuse of scaling factor vs. original T2Processing the distribution data to obtain T of the bound fluid and the movable fluid2And (4) distribution. By utilizing the embodiments of the specification, the continuous quantitative characterization of the distribution of the reservoir bound fluid and the mobile fluid can be accurately realized.
It should be noted that, the above-mentioned apparatus or system in this specification may also include other implementation manners according to the description of the related method embodiment, and a specific implementation manner may refer to the description of the method embodiment, which is not described herein in detail. The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class, storage medium + program embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and for the relevant points, refer to the partial description of the method embodiment.
Although the description of operations and data such as SVD inversion, BRD inversion, etc. acquisition, definition, interaction, computation, judgment, etc. are referred to in the context of embodiments of the present specification, embodiments of the present specification are not limited to what must be consistent with a standard data model/template or described in embodiments of the present specification. Certain industry standards, or implementations modified slightly from those described using custom modes or examples, may also achieve the same, equivalent, or similar, or other, contemplated implementations of the above-described examples. The embodiments using these modified or transformed data acquisition, storage, judgment, processing, etc. may still fall within the scope of the alternative embodiments of the present description.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. One typical implementation device is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a vehicle-mounted human-computer interaction device, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, when implementing one or more of the present description, the functions of each module may be implemented in one or more software and/or hardware, or a module implementing the same function may be implemented by a combination of multiple sub-modules or sub-units, etc. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may therefore be considered as a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method or apparatus that comprises the element.
As will be appreciated by one skilled in the art, one or more embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, one or more embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, one or more embodiments of the present description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
One or more embodiments of the present description may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. One or more embodiments of the present specification can also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description of the specification, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
The above description is only an example of the present specification, and is not intended to limit the present specification. Various modifications and alterations to this description will become apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present specification should be included in the scope of the claims of the present specification.
Claims (7)
1. A reservoir tiedown and mobile fluid distribution determination method, comprising:
acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not; wherein the types comprise a first type and a second type; correspondingly, the determining the type of the original transverse relaxation time distribution data includes: if the original transverse relaxation time distribution data contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a first type; if the original transverse relaxation time distribution data do not contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a second type;
processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid, wherein the proportional coefficients comprise: if the type of the original transverse relaxation time distribution data is a first type, processing the original transverse relaxation time distribution data by using a first proportional coefficient calculation model as follows:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiDetermining α corresponding to the minimum i satisfying the first scale factor calculation modeliIs recorded as αopt1A cover αopt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1) (ii) a If the type of the original transverse relaxation time distribution data is a second type, processing the original transverse relaxation time distribution data by using a second proportionality coefficient calculation model as follows:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklRepresents the long relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiDetermining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2A cover αopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)；
And processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
2. The method of claim 1, wherein the proportionality coefficients for the bound and mobile fluids comprise an exponential function based on a natural constant.
3. The method of claim 2, wherein the proportionality coefficients for the bound and mobile fluids comprise:
wherein, Pb(T2α) is a bound fluid proportionality coefficient distribution function, Pm(T2α) is a movable fluid proportionality coefficient distribution function and α is a form factor for controlling Pb(T2α) and Pm(T2α).
4. The method of claim 1, wherein the processing of the raw transverse relaxation time distribution data using the proportionality coefficients of the bound and mobile fluids, respectively, comprises:
fb(T2)＝Pb(T2,αopt)·f(T2)
fm(T2)＝Pm(T2,αopt)·f(T2)
wherein f isb(T2) For confining the transverse relaxation time distribution data of the fluid, Pb(T2,αopt) To bound the fluid proportionality coefficient, fm(T2) For transverse relaxation time distribution data of the mobile fluid, Pm(T2,αopt) Is a mobile fluid proportionality coefficient, f (T)2) Raw transverse relaxation time distribution data.
5. A reservoir tiedown and mobile fluid distribution determination apparatus, the apparatus comprising:
the data acquisition module is used for acquiring nuclear magnetic resonance echo data of a target work area and inverting the echo data to obtain original transverse relaxation time distribution data of the target work area;
the type determining module is used for determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contains a short relaxation peak or not; wherein the types comprise a first type and a second type; correspondingly, the determining the type of the original transverse relaxation time distribution data includes: if the original transverse relaxation time distribution data contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a first type; if the original transverse relaxation time distribution data do not contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a second type;
a proportionality coefficient determining module, configured to process the original transverse relaxation time distribution data according to the proportionality coefficient calculation model corresponding to the type, so as to obtain proportionality coefficients of a bound fluid and a movable fluid, where the proportionality coefficient determining module is configured to: if the type of the original transverse relaxation time distribution data is a first type, processing the original transverse relaxation time distribution data by using a first proportional coefficient calculation model as follows:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiDetermining α corresponding to the minimum i satisfying the first scale factor calculation modeliIs recorded as αopt1A cover αopt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1) (ii) a If the type of the original transverse relaxation time distribution data is a second type, processing the original transverse relaxation time distribution data by using a second proportionality coefficient calculation model as follows:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklRepresents the long relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiDetermining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2A cover αopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)；
And the fluid distribution determining module is used for processing the original transverse relaxation time distribution data by utilizing the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain the bound fluid distribution data and the movable fluid distribution data of the target work area.
6. A reservoir tiedown and mobilizable fluid distribution determination device comprising a processor and a memory for storing processor-executable instructions that when executed by the processor implement steps comprising:
acquiring nuclear magnetic resonance echo data of a target work area, and performing inversion on the echo data to obtain original transverse relaxation time distribution data of the target work area;
determining the type of the original transverse relaxation time distribution data according to whether the original transverse relaxation time distribution data contain short relaxation peaks or not; wherein the types comprise a first type and a second type; correspondingly, the determining the type of the original transverse relaxation time distribution data includes: if the original transverse relaxation time distribution data contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a first type; if the original transverse relaxation time distribution data do not contain short relaxation peaks, determining that the type of the original transverse relaxation time distribution data belongs to a second type;
processing the original transverse relaxation time distribution data according to the proportional coefficient calculation model corresponding to the type to obtain the proportional coefficients of the bound fluid and the movable fluid, wherein the proportional coefficients comprise: if the type of the original transverse relaxation time distribution data is a first type, processing the original transverse relaxation time distribution data by using a first proportional coefficient calculation model as follows:
|pks[Pb(T2,αi)·f(T2)]-pks[f(T2)]|≤1
wherein, αi＝i·Δα，Δα>0，i＝1,2,…,N，αiRepresenting the ith α value, N the maximum number of α values, Δ α the value interval of α, α the shape factor, for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pksRepresents the short relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,1is a transverse relaxation time distribution threshold of the first type, Pb(T2,αi) Is αiA corresponding bound fluid proportionality coefficient; meter for determining first scale factor satisfying the aboveα for minimum i of computational modeliIs recorded as αopt1A cover αopt1Corresponding Pb(T2,αopt1) And Pm(T2,αopt1) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pm(T2,αopt1)＝1-Pb(T2,αopt1) (ii) a If the type of the original transverse relaxation time distribution data is a second type, processing the original transverse relaxation time distribution data by using a second proportionality coefficient calculation model as follows:
|pkl[Pm(T2,αj)·f(T2)]-pkl[f(T2)]|≤2
wherein, αj＝(N+1-j)·Δα，Δα>0，j＝1,2,…,N，αjRepresenting the jth α value, N representing the maximum number of values α, Δ α representing the value interval α, α being a shape factor for controlling the function shape of the bound fluid proportionality coefficient and the movable fluid proportionality coefficient, pklRepresents the long relaxation peak-to-peak solving function, f (T)2) For the original transverse relaxation time distribution data,2is a transverse relaxation time distribution threshold of the second type, Pm(T2,αi) Is αiDetermining α corresponding to the minimum j satisfying the second scale factor calculation modeljIs recorded as αopt2A cover αopt2Corresponding Pb(T2,αopt2) And Pm(T2,αopt2) Respectively determined as a bound fluid proportionality coefficient and a mobile fluid proportionality coefficient, wherein Pb(T2,αopt2)＝1-Pm(T2,αopt2)；
And processing the original transverse relaxation time distribution data by using the proportionality coefficients of the bound fluid and the movable fluid respectively to obtain bound fluid distribution data and movable fluid distribution data of the target work area.
7. A reservoir tiedown and mobile fluid distribution determination system, the system comprising at least one processor and memory storing computer executable instructions which when executed by the processor implement the steps of the method of any one of claims 1 to 4.
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 CN110244369A (en) | 2019-09-17 |
| CN110244369B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114046145B (en) * | 2021-11-26 | 2023-06-16 | 中国石油大学(北京) | Reservoir fluid identification and saturation determination method and device |
Citations (9)
| 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 |
| 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 |
| 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104330433B (en) * | 2014-10-28 | 2016-08-03 | 中国石油天然气股份有限公司 | A kind of method and device obtaining the distribution of purpose reservoir T2 |
| CN106249306A (en) * | 2016-10-12 | 2016-12-21 | 贵州大学 | Shale pore structure detection method based on nuclear magnetic resonance, NMR |
| CN109594971B (en) * | 2018-12-21 | 2022-07-05 | 中国石油天然气集团有限公司 | Fluid property identification method based on nuclear magnetic resonance logging enhanced diffusion gas layer identification factor |
-
2019
- 2019-06-28 CN CN201910572371.XA patent/CN110244369B/en active Active
Patent Citations (9)
| 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 |
| 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 |
| 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 |
Non-Patent Citations (3)
| Title |
|---|
| "Advanced fluid-typing methods for NMR logging";谢然红 等;《Pet. Sci》;20111231(第2（2011）期);第163-169页 * |
| "利用核磁共振原始回波数据确定可动水饱和度的方法";丁业娇 等;《2015中国地球科学联合学术年会论文集》;20151231;第1984第1段-1987页第1段 * |
| "基于核磁共振新参数的致密砂岩储层孔隙结构特征";代齐全 等;《石油学报》;20160731;第37卷(第7期);第887-897页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110244369A (en) | 2019-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Slope reliability analysis by updated support vector machine and Monte Carlo simulation | |
| Fayazi et al. | Estimation of pseudo relative permeability curves for a heterogeneous reservoir with a new automatic history matching algorithm | |
| Kim et al. | Dependence of pore-to-core up-scaled reaction rate on flow rate in porous media | |
| Jerauld et al. | Validation of a workflow for digitally measuring relative permeability | |
| CN110244369B (en) | Reservoir constraint and movable fluid distribution determination method, device and system | |
| CN108227033B (en) | Method and device for predicting porosity of sandstone reservoir and computer storage medium | |
| Yeh et al. | Inspection of two sophisticated models for sand based on generalized plasticity: Monotonic loading and Monte Carlo analysis | |
| loannidis et al. | Computer generation and application of 3-D model porous media: from pore-level geostatistics to the estimation of formation factor | |
| Chen et al. | Numerical simulation of unstable suction transients in unsaturated soils: the role of wetting collapse | |
| Sanford | Uncertainties in sediment erodibility estimates due to a lack of standards for experimental protocols and data interpretation | |
| Suzuki et al. | Characterization of 3D printed fracture networks | |
| CN106980009A (en) | Shale oil reservoir evaluation methods and device | |
| Klemin et al. | Digital rock technology for quantitative prediction of acid stimulation efficiency in carbonates | |
| Manasipov et al. | Automated interpretation tool for synchronous history matching of multiple SCAL experiments with advance Nurbs representations of relevant functions | |
| Boek | Pore scale simulation of flow in porous media using lattice-Boltzmann computer simulations | |
| CN109800508A (en) | The calculation method and terminal device of the empty top plate thickness at rock-socketed piles end | |
| Hamid et al. | Enhanced assessment of stress dependent permeability of sandstone rock | |
| CN114486978A (en) | Method and system for quantitatively analyzing residual oil of rock core | |
| Satti et al. | From computerized tomography to CFD: production flow prediction in realistic perforation tunnels | |
| CN108874735B (en) | Method and device for determining pressure of ancient fluid in sedimentary basin | |
| Ferreira et al. | Use of emulator methodology for uncertainty reduction quantification | |
| Boek et al. | Multiphase flow and reactive transport at the pore scale using lattice-Boltzmann computer simulations | |
| Meshgi et al. | Analysing tension infiltrometer data from sloped surface using two‐dimensional approximation | |
| Li et al. | Conditional stochastic moment equations for uncertainty analysis of flow in heterogeneous reservoirs | |
| Badarch et al. | A free-surface immersed-boundary lattice Boltzmann method for flows in porous media |
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 |