CN113567485A - Method and device for determining surface relaxation rate of compact sandstone - Google Patents
Method and device for determining surface relaxation rate of compact sandstone Download PDFInfo
- Publication number
- CN113567485A CN113567485A CN202010353466.5A CN202010353466A CN113567485A CN 113567485 A CN113567485 A CN 113567485A CN 202010353466 A CN202010353466 A CN 202010353466A CN 113567485 A CN113567485 A CN 113567485A
- Authority
- CN
- China
- Prior art keywords
- rotating speeds
- different rotating
- saturated oil
- oil core
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
The invention provides a method and a device for determining the surface relaxation rate of compact sandstone, wherein the method comprises the following steps: obtaining a first nuclear magnetic T2 spectrum cumulative integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotating speeds; a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds; determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds; determining the pore radius of the saturated oil core sample corresponding to different rotating speeds; and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds. The method can determine the surface relaxation rate of the compact sandstone, and has strong operability and high accuracy.
Description
Technical Field
The invention relates to the technical field of increasing the recovery ratio of a tight sandstone oil reservoir, in particular to a method and a device for determining the surface relaxation rate of tight sandstone.
Background
The pore type of the compact sandstone reservoir is mainly micro-nano grade, and the fluid seepage characteristic is complex and difficult to quantitatively describe. The low-field nuclear magnetic T2 spectrum is an effective means for determining the macroscopic statistical law of the microfluid in the pores, and the relaxation time T2 is closely related to the pore distribution and can reflect the pore distribution characteristics to a certain extent.
In a uniformly distributed magnetic field, the low field NMR relaxation time T is not taken into account for the influence of diffusion relaxation and free relaxation (negligible compared to the influence of surface relaxation)2The following relationship can be established with the pore radius:
in the formula, T2Is the relaxation time, ms; rho is the surface relaxation rate, mu m/s; s is the surface area of the core in cm2(ii) a V is pore volume, cm3(ii) a R is the pore radius, cm; c is a constant, C ═ 1, 2, 3 are used for the flat plate model, the capillary bundle model and the spherical model, respectively, and for the tight sandstone, the capillary bundle model is generally selected, i.e., C ═ 2.
Therefore, the relaxation time T2 can be converted into a pore radius by calculating the surface relaxation rate, and the pore distribution law can be characterized by using a T2 spectrum. For the same core sample, the surface relaxation rate calculation methods can be roughly divided into three categories: (1) direct calculation using low-field nuclear magnetic test results, i.e. using the mode of CPMG or IR pulse sequence and limited diffusion pulse sequence test (such as pulse field gradient and pulse field gradient excitation echo); (2) combining the measurement results of the pore radius (including mercury porosimetry and nitrogen adsorption) and the relaxation time (low-field nuclear magnetic method) to calculate the surface relaxation rate; (3) and (3) combining the test results of specific surface area (including nitrogen adsorption, cation exchange capacity and imaging analysis) and relaxation time (low-field nuclear magnetic method) to calculate the surface relaxation rate.
The first method is a direct measurement method, and is not suitable for unconventional tight sandstone, because in the rock core with a large number of micro-nano pores, the diffusion relaxation phenomenon has obvious influence, the rapid diffusion relaxation of hydrogen protons is more complex, and the influence on the test result is larger. The second and third methods are indirect measurement methods, and the measurement results depend on which measurement means is used to evaluate the porosity, surface area, pore distribution, and the like of the tight sandstone. The method suitable for evaluating the pore distribution of the tight sandstone is a high-pressure mercury intrusion measurement technology, and the surface relaxation rate can be solved by combining the pore radius measured by the high-pressure mercury intrusion measurement and the relaxation time measured by low-field nuclear magnetism. However, the high-pressure mercury intrusion test process is time-consuming and labor-consuming, the operation flow is complex, and the core is permanently damaged.
Disclosure of Invention
The embodiment of the invention provides a method for determining the surface relaxation rate of compact sandstone, which is used for determining the surface relaxation rate of compact sandstone and has strong operability and high accuracy, and the method comprises the following steps:
obtaining a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds, wherein the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds.
The embodiment of the invention provides a device for determining the surface relaxation rate of compact sandstone, which is used for determining the surface relaxation rate of the compact sandstone and has strong operability and high accuracy, and the device comprises:
the device comprises a T2 spectrum obtaining module, a sampling module and a data processing module, wherein the T2 spectrum obtaining module is used for obtaining a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds, the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
the quasi-T2 cut-off value obtaining module is used for determining the quasi-T2 cut-off values of the saturated oil core sample corresponding to different rotating speeds according to a first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds;
the pore radius determining module is used for determining the pore radius of the saturated oil core sample corresponding to different rotating speeds;
and the compact sandstone surface relaxation rate determining module is used for determining the compact sandstone surface relaxation rate according to the simulated T2 cut-off value and the pore radius of the saturated oil core sample corresponding to different rotating speeds.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the method for determining the surface relaxation rate of the tight sandstone when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the method for determining the surface relaxation rate of tight sandstone.
In the embodiment of the invention, a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained, wherein the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds; determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds; determining the pore radius of the saturated oil core sample corresponding to different rotating speeds; and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds. In the process, a first nuclear magnetic T2 spectrum accumulated integral curve and a second nuclear magnetic T2 spectrum accumulated integral curve are obtained by directly adopting a low-field nuclear magnetic measurement and high-speed centrifugation process, so that a simulated T2 cut-off value of a saturated oil core sample is determined, and after the pore radii of the saturated oil core sample corresponding to different rotating speeds are determined, the surface relaxation rate of the compact sandstone is finally obtained.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
fig. 1 is a flow chart of a method of determining the surface relaxation rate of tight sandstone in an embodiment of the invention;
fig. 2 is a detailed flowchart of a method for determining the surface relaxation rate of tight sandstone according to an embodiment of the present invention;
3-6 are schematic diagrams of the accumulated signal amplitude of the second nuclear magnetic T2 spectrum corresponding to the compact sandstone plunger samples A11-A14 with different rotating speeds in the embodiment of the invention;
fig. 7-10 are schematic diagrams of fitting functions corresponding to tight sandstone plunger samples a11-a14, respectively, in an example of the invention;
11-14 are comparative illustrations of pore diameters obtained with different methods for tight sandstone plug samples A11-A14, respectively, in accordance with an embodiment of the present invention;
figure 15 is a schematic diagram of an apparatus for determining the surface relaxation rate of tight sands in an embodiment of the present invention;
FIG. 16 is a diagram of a computer device in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
In the description of the present specification, the terms "comprising," "including," "having," "containing," and the like are used in an open-ended fashion, i.e., to mean including, but not limited to. Reference to the description of the terms "one embodiment," "a particular embodiment," "some embodiments," "for example," etc., means that a particular feature, structure, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. The sequence of steps involved in the embodiments is for illustrative purposes to illustrate the implementation of the present application, and the sequence of steps is not limited and can be adjusted as needed.
Fig. 1 is a flowchart of a method for determining a surface relaxation rate of tight sandstone according to an embodiment of the present invention, and as shown in fig. 1, the method includes:
101, obtaining a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds, wherein the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
102, determining simulated T2 cut-off values of the saturated oil core sample corresponding to different rotating speeds according to a first nuclear magnetic T2 spectrum accumulated integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulated integral curves corresponding to different rotating speeds;
103, determining the pore radii of the saturated oil core samples corresponding to different rotating speeds;
and 104, determining the surface relaxation rate of the compact sandstone according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds.
According to the method provided by the embodiment of the invention, a first nuclear magnetic T2 spectrum accumulated integral curve and a second nuclear magnetic T2 spectrum accumulated integral curve are obtained by directly adopting a low-field nuclear magnetic measurement and high-speed centrifugation process, so that the quasi T2 cutoff value of a saturated oil core sample is determined, and after the pore radii of the saturated oil core sample corresponding to different rotating speeds are determined, the surface relaxation rate of the tight sandstone is finally obtained.
In specific implementation, in step 101, when performing low-field nuclear magnetic measurement, a low-field nuclear magnetic device may be selected, and when performing high-speed centrifugation, a high-speed refrigerated centrifuge may be selected, and the second nuclear magnetic T2 spectrum cumulative integral curve corresponding to the set rotation speed is obtained by centrifuging at the set rotation speed for 2 hours. And then, increasing the rotation speed of the high-speed refrigerated centrifuge, and repeating the steps until the measured second nuclear magnetic T2 spectrum cumulative integral curve hardly changes (for example, the change of the T2 spectrum cumulative integral area is less than 3%), so that a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotation speeds can be obtained.
In one embodiment, determining quasi-T2 cut-off values of saturated oil core samples corresponding to different rotating speeds according to a first nuclear magnetic T2 spectrum cumulative integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotating speeds includes:
and determining a T2 value corresponding to the intersection point of the first nuclear magnetic T2 spectrum cumulative integral curve and the second nuclear magnetic T2 spectrum cumulative integral curve corresponding to different rotating speeds as the quasi-T2 cut-off value of the saturated oil core sample corresponding to different rotating speeds.
In the above embodiment, for each second nuclear magnetic T2 spectrum cumulative integral curve, the straight line segment of the second nuclear magnetic T2 spectrum cumulative integral curve is extended reversely to intersect with the first nuclear magnetic T2 spectrum cumulative integral curve, and the T2 value at the intersection point is the pseudo T2 cutoff value. The method for determining the quasi-T2 cut-off value is simple and convenient and has strong operability.
In one embodiment, determining pore radii of saturated oil core samples corresponding to different rotation speeds comprises:
obtaining centrifugal forces corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment;
and determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds.
In the embodiment, according to the principle of determining the pore radius, during high-speed centrifugation, the centrifugal force is equal to the capillary force of the saturated oil core sample, that is, the centrifugal force corresponding to different rotating speeds is the capillary force corresponding to different rotating speeds, so that the pore radius of the saturated oil core sample corresponding to different rotating speeds is determined according to the capillary force corresponding to different rotating speeds.
In one embodiment, the centrifugal forces corresponding to different rotation speeds when the saturated oil core sample is subjected to high-speed centrifugation are obtained by the following formula:
wherein, PcAs a centrifugal forceIn MPa; l is the length of the saturated oil core sample, and the unit is cm; reThe rotation radius of a saturated oil core sample is in cm; delta rho is the density difference of oil and gas two phases, and the unit is g/cm 3; n is the rotating speed of the centrifuge, and the unit is r/min;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds by adopting the following formula:
wherein, PciCapillary force in MPa; sigma is the surface tension of oil gas, and the unit is mN/m; θ is the wetting angle in °; r is the pore radius in cm.
In one embodiment, determining the surface relaxation rate of the tight sandstone according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds comprises the following steps:
determining surface relaxation rates corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to the different rotating speeds;
determining a fitting function of the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample according to the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the simulated T2 cut-off value as the compact sandstone surface relaxation rate of the saturated oil core sample when the surface relaxation rate is zero based on the fitting function.
In the above embodiment, a plurality of quasi-T2 cut-off values and a plurality of surface relaxivity values may be determined by means of linear fitting, to form a fitting function, and the surface relaxivity of the fitting function when the quasi-T2 cut-off value (i.e., the coordinate axis of the T2 value) is zero is the tight sandstone surface relaxivity of the saturated oil core sample, i.e., the final tight sandstone surface relaxivity.
In one embodiment, the surface relaxation rates corresponding to different rotating speeds are determined according to the quasi-T2 cutoff values and the pore radii of the saturated oil core samples corresponding to the different rotating speeds by adopting the following formula:
wherein, T2,pcutoffIs a quasi-T2 cutoff value in ms; rho is the surface relaxation rate and has the unit of mu m/s; r is the pore radius in cm; and C is 2.
In addition, a T2 spectrum can be converted into pore radius distribution according to the surface relaxation rate of the compact sandstone, the pore radius distribution is compared with pore distribution results determined by an average value method and a high-pressure mercury pressing method, when the pore radius distribution results are compared, large samples belonging to the same compact sandstone plunger can be sliced, the first part adopts the average value method to obtain the pore radius distribution, the second part adopts the high-pressure mercury pressing method to obtain the pore radius distribution, and the third part adopts the method provided by the invention to obtain the pore radius distribution. Through comparison, the pore radius distribution obtained by the embodiment of the invention is determined to be highly consistent with the pore radius distribution obtained by an average value method and basically consistent with the general trend of a high-pressure mercury porosimetry method, so that the finally determined surface relaxation rate calculation result is real and reliable.
After the pore radius distribution is obtained by adopting an average value method, the formula for calculating the surface relaxation rate by adopting the average value method is as follows:
wherein, T2LMThe mean value of the logarithm of the relaxation time is ms; rpIs the average pore radius in μm; r isiThe radius of the pore at the ith point is the unit of mu m; siIs the i-th point mercury saturation in%.
Based on the above embodiment, the present invention provides the following embodiment to illustrate a detailed flow of the method for determining the surface relaxation rate of tight sandstone, and fig. 2 is a detailed flow chart of the method for determining the surface relaxation rate of tight sandstone according to the embodiment of the present invention, which includes:
and step 207, determining the surface relaxation rate when the simulated T2 cut-off value is zero as the compact sandstone surface relaxation rate of the saturated oil core sample based on the fitting function.
Of course, it is understood that other variations of the above detailed flow can be made, and all such variations are intended to fall within the scope of the present invention.
A specific example is given below to illustrate the specific application of the method proposed by the present invention.
Selecting four compact sandstone plunger samples with the length of 632 of the Ordos basin extension group main force development layer, wherein the coring depth is 2075-2225m, the lithology of the samples is mainly quartz feldspar sandstone, and the clay mineral composition comprises illite (45.8%), chlorite (55.4%) and an illite/montmorillonite mixed layer (8.1%);
after oil washing and drying treatment, cutting a section (with the length of 1.8-2.0cm) at the end face of the compact sandstone plunger sample for high-pressure mercury intrusion test, taking another section (with the length of 1-1.2cm) for contact angle test (contact angle method), and performing vacuum saturated oil treatment (vacuumizing for 48 hours, the confining pressure of 30MPa and the kerosene saturation for 5 days) on the rest (including 4 parts of A11-A14 and the length of 3.6-3.8cm) by using a vacuum pressurizing saturation device;
measuring a compact sandstone plunger sample A11-A14 after the compact sandstone saturated oil treatment by using a MicroMR12-025V small-size nuclear magnetic resonance analyzer to obtain a first nuclear magnetic T2 spectrum accumulation integral curve corresponding to each part of the compact sandstone plunger sample, then placing each part of the compact sandstone plunger sample in a CSC-12(S) super core high-speed refrigerated centrifuge, centrifuging for 2 hours at the rotating speed of 3000rpm/min, taking out, and measuring a second nuclear magnetic T2 spectrum accumulation integral curve corresponding to the rotating speed of 3000 rpm/min;
and increasing the rotating speed to 4000rpm, and repeating the steps until the measured morphology of the second nuclear magnetic T2 spectrum cumulative integration curve is hardly changed (the change of the T2 spectrum cumulative integration area is less than 3%), wherein in the figures 3-6, the graphs are respectively a schematic diagram of the cumulative signal amplitude of the second nuclear magnetic T2 spectrum corresponding to different rotating speeds of the compact sandstone plunger samples A11-A14 in the embodiment of the invention, and the cumulative signal amplitude is the cumulative integration curve, and in the figures 3-6, the rotating speeds corresponding to the cumulative integration curves from top to bottom are respectively 0rpm, 1000rpm, 3000rpm, 5000rpm, 6000rpm, 7000rpm, 8000rpm and 9000 rpm.
For each part of the compact sandstone plunger sample, for each second nuclear magnetic T2 spectrum accumulation integral curve, reversely extending a straight line segment of the second nuclear magnetic T2 spectrum accumulation integral curve, intersecting with the first nuclear magnetic T2 spectrum accumulation integral curve, wherein a T2 value at an intersection point is a quasi-T2 cut-off value;
and (4) calculating the pore radius of each part of the compact sandstone plunger sample corresponding to different rotating speeds according to the formula (2) and the formula (3).
Then calculating the surface relaxation rate of each part of the tight sandstone plunger sample corresponding to different rotating speeds by adopting a formula (4), and determining a fitting function of the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample according to the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample corresponding to different rotating speeds, wherein FIGS. 7-10 are schematic diagrams of the fitting functions corresponding to the tight sandstone plunger samples A11-A14 in the embodiment of the invention respectively; and determining the surface relaxation rate of the simulated T2 cut-off value as the tight sandstone surface relaxation rate of the saturated oil core sample based on the fitting function, so as to determine the final tight sandstone surface relaxation rate of the tight sandstone plunger sample A11-A14.
Fig. 11-14 are schematic diagrams respectively comparing pore diameters of the tight sandstone plunger samples a11-a14 obtained by different methods in the examples of the invention, wherein the method of the invention is a T2-like cut-off value method in fig. 11-14, and table 1 shows the results of surface relaxivity of the tight sandstone plunger samples a11-a14 obtained by different methods, and the error between the two results is within an acceptable range.
TABLE 1 results of the mean value method and the surface relaxation rate of the method of the invention
In summary, in the method provided by the embodiment of the present invention, a first nuclear magnetic T2 spectrum cumulative integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotation speeds are obtained, where the saturated oil core sample is a tight sandstone plunger sample after vacuum pumping saturated oil treatment, and the first nuclear magnetic T2 spectrum cumulative integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds; determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds; determining the pore radius of the saturated oil core sample corresponding to different rotating speeds; and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds. In the process, a first nuclear magnetic T2 spectrum accumulated integral curve and a second nuclear magnetic T2 spectrum accumulated integral curve are obtained by directly adopting a low-field nuclear magnetic measurement and high-speed centrifugation process, so that a simulated T2 cut-off value of a saturated oil core sample is determined, and after the pore radii of the saturated oil core sample corresponding to different rotating speeds are determined, the surface relaxation rate of the compact sandstone is finally obtained.
The embodiment of the invention also provides a device for determining the surface relaxation rate of the compact sandstone, and the principle of the device is similar to that of the method for determining the surface relaxation rate of the compact sandstone, and the device is not repeated here.
Fig. 15 is a schematic diagram of an apparatus for determining the surface relaxation rate of tight sandstone according to an embodiment of the present invention, as shown in fig. 15, the apparatus includes:
the T2 spectrum obtaining module 1501 is configured to obtain a first nuclear magnetic T2 spectrum cumulative integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotation speeds, where the saturated oil core sample is a tight sandstone plunger sample subjected to vacuumized saturated oil treatment, and the first nuclear magnetic T2 spectrum cumulative integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
a quasi-T2 cutoff value obtaining module 1502 for determining quasi-T2 cutoff values of the saturated oil core sample corresponding to different rotation speeds according to a first nuclear magnetic T2 spectrum cumulative integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotation speeds;
the pore radius determining module 1503 is used for determining pore radii of the saturated oil core samples corresponding to different rotating speeds;
and the compact sandstone surface relaxation rate determining module 1504 is used for determining the compact sandstone surface relaxation rate according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds.
In an embodiment, the pseudo T2 cutoff obtaining module 1502 is specifically configured to:
and determining a T2 value corresponding to the intersection point of the first nuclear magnetic T2 spectrum cumulative integral curve and the second nuclear magnetic T2 spectrum cumulative integral curve corresponding to different rotating speeds as the quasi-T2 cut-off value of the saturated oil core sample corresponding to different rotating speeds.
In one embodiment, the aperture radius determination module 1503 is specifically configured to:
obtaining centrifugal forces corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment;
and determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds.
In one embodiment, the aperture radius determination module 1503 is specifically configured to:
the centrifugal force corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment is obtained by adopting the following formula:
wherein, PcIs the centrifugal force, in MPa; l is the length of the saturated oil core sample, and the unit is cm; reThe rotation radius of a saturated oil core sample is in cm; delta rho is the density difference of oil and gas two phases, and the unit is g/cm 3; n is the rotating speed of the centrifuge, and the unit is r/min;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds by adopting the following formula:
wherein, PciCapillary force in MPa; sigma is the surface tension of oil gas, and the unit is mN/m; θ is the wetting angle in °; r is the pore radius in cm.
In an embodiment, the tight sandstone surface relaxation rate determination module 1504 is specifically configured to:
determining surface relaxation rates corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to the different rotating speeds;
determining a fitting function of the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample according to the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the simulated T2 cut-off value as the compact sandstone surface relaxation rate of the saturated oil core sample when the surface relaxation rate is zero based on the fitting function.
In an embodiment, the tight sandstone surface relaxation rate determination module 1504 is specifically configured to:
determining the surface relaxation rates of the compact sandstone corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds by adopting the following formula:
wherein, T2,pcutoffIs a quasi-T2 cutoff value in ms; rho is the surface relaxation rate and has the unit of mu m/s; r is the pore radius in cm; and C is 2.
In summary, in the apparatus provided in the embodiment of the present invention, a first nuclear magnetic T2 spectrum cumulative integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotation speeds are obtained, where the saturated oil core sample is a tight sandstone plunger sample after vacuum pumping saturated oil treatment, and the first nuclear magnetic T2 spectrum cumulative integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds; determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds; determining the pore radius of the saturated oil core sample corresponding to different rotating speeds; and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds. In the process, a first nuclear magnetic T2 spectrum accumulated integral curve and a second nuclear magnetic T2 spectrum accumulated integral curve are obtained by directly adopting a low-field nuclear magnetic measurement and high-speed centrifugation process, so that a simulated T2 cut-off value of a saturated oil core sample is determined, and after the pore radii of the saturated oil core sample corresponding to different rotating speeds are determined, the surface relaxation rate of the compact sandstone is finally obtained.
An embodiment of the present application further provides a computer device, fig. 16 is a schematic diagram of the computer device in the embodiment of the present invention, where the computer device is capable of implementing all steps in the method for determining a surface relaxation rate of tight sandstone, and the electronic device specifically includes the following contents:
a processor (processor)1601, a memory (memory)1602, a communication Interface (Communications Interface)1603, and a bus 1604;
the processor 1601, the memory 1602 and the communication interface 1603 are configured to communicate with each other via the bus 1604; the communication interface 1603 is used for realizing information transmission among relevant devices such as server-side devices, detection devices and user-side devices;
the processor 1601 is configured to call a computer program in the memory 1602, and the processor executes the computer program to implement all the steps of the method for determining the tight sandstone surface relaxation rate in the above embodiment.
Embodiments of the present application also provide a computer readable storage medium, which can implement all the steps of the method for determining tight sandstone surface relaxation rate in the above embodiments, and the computer readable storage medium stores a computer program, which when executed by a processor implements all the steps of the method for determining tight sandstone surface relaxation rate in the above embodiments.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention 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.
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.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (14)
1. A method of determining the surface relaxation rate of tight sandstone, comprising:
obtaining a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds, wherein the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
determining the quasi-T2 cutoff values of the saturated oil core samples corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the compact sandstone according to the simulated T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds.
2. The method for determining the surface relaxation rate of tight sandstone according to claim 1, wherein the determining the pseudo-T2 cutoff values of the saturated oil core sample corresponding to different rotating speeds according to the first nuclear magnetic T2 spectrum cumulative integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum cumulative integral curves corresponding to different rotating speeds comprises:
and determining a T2 value corresponding to the intersection point of the first nuclear magnetic T2 spectrum cumulative integral curve and the second nuclear magnetic T2 spectrum cumulative integral curve corresponding to different rotating speeds as the quasi-T2 cut-off value of the saturated oil core sample corresponding to different rotating speeds.
3. The method for determining the surface relaxation rate of tight sandstone according to claim 2, wherein determining the pore radius of the saturated oil core sample corresponding to different rotating speeds comprises:
obtaining centrifugal forces corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment;
and determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds.
4. The method for determining the surface relaxation rate of tight sandstone according to claim 3, wherein the centrifugal force corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugation is obtained by adopting the following formula:
wherein, PcIs the centrifugal force, in MPa; l is the length of the saturated oil core sample, and the unit is cm; reThe rotation radius of a saturated oil core sample is in cm; delta rho is the density difference of oil and gas two phases, and the unit is g/cm 3; n is the rotating speed of the centrifuge, and the unit is r/min;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds by adopting the following formula:
wherein, PciCapillary force in MPa; sigma is the surface tension of oil gas, and the unit is mN/m; θ is the wetting angle in °; r is the pore radius in cm.
5. The method for determining the surface relaxation rate of tight sandstone according to claim 3, wherein the determining the surface relaxation rate of tight sandstone according to the pseudo-T2 cut-off value and the pore radius of the saturated oil core sample corresponding to different rotating speeds comprises the following steps:
determining surface relaxation rates corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to the different rotating speeds;
determining a fitting function of the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample according to the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the simulated T2 cut-off value as the compact sandstone surface relaxation rate of the saturated oil core sample when the surface relaxation rate is zero based on the fitting function.
6. The method for determining the surface relaxation rate of tight sandstone according to claim 5, wherein the surface relaxation rate of tight sandstone corresponding to different rotating speeds is determined according to the quasi-T2 cutoff value and the pore radius of the saturated oil core sample corresponding to different rotating speeds by adopting the following formula:
wherein, T2,pcutoffIs a quasi-T2 cutoff value in ms; rho is the surface relaxation rate and has the unit of mu m/s; r is the pore radius in cm; and C is 2.
7. An apparatus for determining the surface relaxation rate of tight sandstone, comprising:
the device comprises a T2 spectrum obtaining module, a sampling module and a data processing module, wherein the T2 spectrum obtaining module is used for obtaining a first nuclear magnetic T2 spectrum accumulation integral curve of a saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds, the saturated oil core sample is a compact sandstone plunger sample subjected to vacuumizing saturated oil treatment, and the first nuclear magnetic T2 spectrum accumulation integral curve is obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample; the second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds are obtained by performing low-field nuclear magnetic measurement on the saturated oil core sample subjected to high-speed centrifugal treatment at different rotating speeds;
the quasi-T2 cut-off value obtaining module is used for determining the quasi-T2 cut-off values of the saturated oil core sample corresponding to different rotating speeds according to a first nuclear magnetic T2 spectrum accumulation integral curve of the saturated oil core sample and a plurality of second nuclear magnetic T2 spectrum accumulation integral curves corresponding to different rotating speeds;
the pore radius determining module is used for determining the pore radius of the saturated oil core sample corresponding to different rotating speeds;
and the compact sandstone surface relaxation rate determining module is used for determining the compact sandstone surface relaxation rate according to the simulated T2 cut-off value and the pore radius of the saturated oil core sample corresponding to different rotating speeds.
8. The apparatus for determining densified sandstone surface relaxation rate of claim 7, wherein the pseudo-T2 cutoff value obtaining module is specifically configured to:
and determining a T2 value corresponding to the intersection point of the first nuclear magnetic T2 spectrum cumulative integral curve and the second nuclear magnetic T2 spectrum cumulative integral curve corresponding to different rotating speeds as the quasi-T2 cut-off value of the saturated oil core sample corresponding to different rotating speeds.
9. The apparatus for determining densified sandstone surface relaxation rate of claim 8, wherein the pore radius determination module is specifically configured to:
obtaining centrifugal forces corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment;
and determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds.
10. The apparatus for determining densified sandstone surface relaxation rate of claim 9, wherein the pore radius determination module is specifically configured to:
the centrifugal force corresponding to different rotating speeds when the saturated oil core sample is subjected to high-speed centrifugal treatment is obtained by adopting the following formula:
wherein, PcIs the centrifugal force, in MPa; l is the length of the saturated oil core sample, and the unit is cm; reThe rotation radius of a saturated oil core sample is in cm; delta rho is the density difference of oil and gas two phases, and the unit is g/cm 3; n is the rotating speed of the centrifuge, and the unit is r/min;
determining the pore radius of the saturated oil core sample corresponding to different rotating speeds according to the centrifugal force corresponding to different rotating speeds by adopting the following formula:
wherein, PciCapillary force in MPa; sigma is the surface tension of oil gas, and the unit is mN/m; θ is the wetting angle in °; r is the pore radius in cm.
11. The apparatus for determining compacted sandstone surface relaxation rate of claim 9, wherein the compacted sandstone surface relaxation rate determination module is specifically configured to:
determining surface relaxation rates corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to the different rotating speeds;
determining a fitting function of the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample according to the quasi-T2 cutoff value and the surface relaxation rate of the saturated oil core sample corresponding to different rotating speeds;
and determining the surface relaxation rate of the simulated T2 cut-off value as the compact sandstone surface relaxation rate of the saturated oil core sample when the surface relaxation rate is zero based on the fitting function.
12. The apparatus for determining compacted sandstone surface relaxation rate of claim 11, wherein the compacted sandstone surface relaxation rate determination module is specifically configured to:
determining the surface relaxation rates of the compact sandstone corresponding to different rotating speeds according to the quasi-T2 cut-off values and the pore radii of the saturated oil core samples corresponding to different rotating speeds by adopting the following formula:
wherein, T2,pcutoffIs a quasi-T2 cutoff value in ms; rho is the surface relaxation rate and has the unit of mu m/s; r is the pore radius in cm; and C is 2.
13. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 6 when executing the computer program.
14. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010353466.5A CN113567485A (en) | 2020-04-29 | 2020-04-29 | Method and device for determining surface relaxation rate of compact sandstone |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010353466.5A CN113567485A (en) | 2020-04-29 | 2020-04-29 | Method and device for determining surface relaxation rate of compact sandstone |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113567485A true CN113567485A (en) | 2021-10-29 |
Family
ID=78158305
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010353466.5A Pending CN113567485A (en) | 2020-04-29 | 2020-04-29 | Method and device for determining surface relaxation rate of compact sandstone |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113567485A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160139291A1 (en) * | 2014-11-13 | 2016-05-19 | Colorado School Of Mines | Surface Relaxivity Calculation using Nuclear Magnetic Resonance (NMR) Measurement, Three Dimensional (3D) Rock Model and NMR Response Simulation |
| CN106990131A (en) * | 2017-02-21 | 2017-07-28 | 上海大学 | A kind of nano-fluid changes the fast appraisement method of core wettability |
| CN110108616A (en) * | 2019-04-24 | 2019-08-09 | 中国地质大学（北京） | Coal sample centrifugal process loss of signal compensation method in a kind of pair of nuclear magnetic resonance experiment |
-
2020
- 2020-04-29 CN CN202010353466.5A patent/CN113567485A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160139291A1 (en) * | 2014-11-13 | 2016-05-19 | Colorado School Of Mines | Surface Relaxivity Calculation using Nuclear Magnetic Resonance (NMR) Measurement, Three Dimensional (3D) Rock Model and NMR Response Simulation |
| CN106990131A (en) * | 2017-02-21 | 2017-07-28 | 上海大学 | A kind of nano-fluid changes the fast appraisement method of core wettability |
| CN110108616A (en) * | 2019-04-24 | 2019-08-09 | 中国地质大学（北京） | Coal sample centrifugal process loss of signal compensation method in a kind of pair of nuclear magnetic resonance experiment |
Non-Patent Citations (3)
| Title |
|---|
| YUN JIANG ET AL.: "Experimental Study on Spontaneous Imbibition under Confining Pressure in Tight Sandstone Cores Based on Low-Field Nuclear Magnetic Resonance Measurements", 《ENERGY & FUELS》, vol. 32, no. 3, pages 3152 - 3162 * |
| 任会康;王安民;李昌峰;曹代勇;魏迎春;: "基于核磁共振技术的低阶煤储层孔隙特征研究", 煤炭科学技术, no. 04, pages 143 - 148 * |
| 贾辰: "基于核磁共振测量岩芯孔隙的研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》, no. 2, pages 005 - 291 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103884633B (en) | A kind of method of definite rock permeability and device | |
| EP3593159B1 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr | |
| US10156531B2 (en) | Measuring connectivity between different pore types in porous media | |
| CN109030292B (en) | novel method for determining wettability of compact rock | |
| US20060132131A1 (en) | Method of measuring rock wettability by means of nuclear magnetic resonance | |
| CN105891089B (en) | Reservoir permeability determines method and device | |
| US20240027379A1 (en) | Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids | |
| US10545105B2 (en) | Estimating formation properties using saturation profiles | |
| CN110927035A (en) | Method for calculating saturation of irreducible tight sandstone bound water | |
| US20050216223A1 (en) | Method of evaluating the capillary pressure curve of an underground deposit rocks based on rock cuttings measurements | |
| CN107991335A (en) | Tight sand Water-blocking damage evaluation test method | |
| CN103257151A (en) | Method for quantitative evaluation on pore throat usage law in oil and gas secondary migration process | |
| CN106644879B (en) | A kind of method and device of determining rock core different aperture component permeability contribution margin | |
| CN104990854A (en) | Method and device for determining irreducible water saturation | |
| CN113075102A (en) | Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time | |
| CN113567485A (en) | Method and device for determining surface relaxation rate of compact sandstone | |
| CN110031496A (en) | A method of evaluation compact reservoir movable fluid distribution characteristics | |
| CN113218834A (en) | Experimental device and method for quantitatively describing seepage damage of dense gas fracturing fluid and reservoir | |
| CN115420766B (en) | Based on T 2 Method for measuring distribution of residual oil in sandstone water injection oil extraction by spectrum analysis | |
| CN116183458A (en) | Shale oil effective porosity determination method | |
| CN112147050A (en) | Method and device for measuring permeability of oil and gas reservoir rock dynamic liquid | |
| CN113376069A (en) | Method and device for determining physical property lower limit of tight reservoir | |
| US20190033238A1 (en) | Estimating formation properties using saturation profiles | |
| CN111693427B (en) | Analysis method for mobility of oil and gas reservoir fluid | |
| CN109298008B (en) | Method, device and system for analyzing oil-water content in porous medium displacement process |
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 |