CN111537544A - Improve nuclear magnetic resonance T2Conversion method for spectral characterization of dense reservoir pore size distribution precision - Google Patents
Improve nuclear magnetic resonance T2Conversion method for spectral characterization of dense reservoir pore size distribution precision Download PDFInfo
- Publication number
- CN111537544A CN111537544A CN202010532575.3A CN202010532575A CN111537544A CN 111537544 A CN111537544 A CN 111537544A CN 202010532575 A CN202010532575 A CN 202010532575A CN 111537544 A CN111537544 A CN 111537544A
- Authority
- CN
- China
- Prior art keywords
- pore
- nuclear magnetic
- magnetic resonance
- size distribution
- pore size
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 33
- 238000012512 characterization method Methods 0.000 title claims abstract description 18
- 230000003595 spectral effect Effects 0.000 title claims abstract description 13
- 238000006243 chemical reaction Methods 0.000 claims abstract description 44
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 33
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 33
- 238000002474 experimental method Methods 0.000 claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 239000008398 formation water Substances 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 4
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 4
- 239000012530 fluid Substances 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 238000003825 pressing Methods 0.000 abstract description 12
- 238000007796 conventional method Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 8
- 238000004836 empirical method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000002791 soaking Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 150000002696 manganese Chemical class 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Geophysics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Remote Sensing (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geochemistry & Mineralogy (AREA)
- Dispersion Chemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for improving nuclear magnetic resonance T2The conversion method for the spectral characterization of the precision of the pore size distribution of the tight reservoir comprises the following steps: the method comprises the following steps: performing nuclear magnetic resonance test on a compact reservoir sample of saturated formation water to obtain nuclear magnetic resonance T2A spectrum; step two: performing a constant-speed mercury-pressing experiment on the dry and compact reservoir sample to obtain mercury-pressing pore size distribution data; step three: carrying out sectional processing on the mercury intrusion pore size distribution data; step four: obtaining T using a cylindrical pore model2The conversion relation with the aperture size; step five: respectively using the conversion relation in the step four to obtain the T of each section2Conversion method from aperture size to obtain complete T2A conversion method for spectral characterization of tight reservoir pore size distribution. The method improves the accuracy of pore size distribution obtained by nuclear magnetic resonance experiments.
Description
Technical Field
The invention relates to a reservoir prediction technology, in particular to a nuclear magnetic resonance characterization method of a compact reservoir pore structure, and belongs to the field of petroleum exploration and development.
Background
One of the main features of tight reservoirs compared to conventional reservoirs is the small pore throat space, which makes migration of fluids therein difficult. Therefore, accurate evaluation of the pore structure of the compact reservoir is one of the preconditions for evaluating the potential, the exploitability and the productivity of compact oil and gas resources, nuclear magnetic resonance is a more common test means, and the method needs to test T in a saturated water state2Converted to a pore radius. The advantages are that (1) the pore size distribution of the fluid can be definitely bound by combining a centrifugal or drying means; (2) by combining the relaxation spectrum of the saturated manganese sample, the pore size distribution occupied by the oil phase can be determined; (3) the control of the throat on the bound fluid or residual oil can be determined by combining a constant-speed mercury pressing experiment.
Usually will T2Conversion to pore size distribution is an experimental approach incorporating mercury injection. The conventional conversion method is derived by the formula: t is2=CrtWhere C is the conversion coefficient, rtPore size obtained for mercury injection. The above methods are described in the prior art of CN104634718A, CN106249306A and the like. However, this method has a problem that T is2And rtThe linear relationship of (a) and (b) does not fit well to both. This is more pronounced in tight reservoirs, probably due to the complex pore throat structure and smaller pore size distribution. For this reason, some researchers have obtained empirical formulas based on conventional methods:wherein n is a constant. But the formula is notCan be obtained through mathematical derivation, and simultaneously lacks of physical significance. Moreover, a common conversion method is based on actually measured mercury intrusion pore size distribution data, if the actually measured data are densely distributed in a certain mercury intrusion saturation section, the conversion formula can take more consideration of the weight of the actually measured data, so that the fitting degree of the actually measured data sparse section is poor, and the utilization of T is caused2The conversion to pore size distribution is less effective.
Thus, the present invention proposes to utilize T2The new method for converting the pore size distribution also fully considers the cylindrical pore model on the basis of common derivation. On the one hand, the transverse relaxation time T is not well influenced by the linear relation2And the pore size distribution conversion is carried out, and on the other hand, the proposed model is more in line with the physical significance. In addition, the processing of the mercury intrusion pore size distribution data by the method can be helpful for improving the accuracy of pore size distribution obtained by nuclear magnetic resonance experiments and evaluating the potential of compact oil and gas resources relatively accurately.
Disclosure of Invention
The invention aims to provide a new material which has certain physical significance and is beneficial to improving nuclear magnetic resonance T2The conversion method for the spectrum characterization of the pore size distribution precision of the compact reservoir relatively accurately evaluates the potential of compact oil and gas resources.
In order to achieve the above purpose, the embodiments of the present invention provide the following technical solutions: improve nuclear magnetic resonance T2The conversion method for the spectral characterization of the precision of the pore size distribution of the tight reservoir comprises the following steps:
the method comprises the following steps: performing nuclear magnetic resonance test on a compact reservoir sample of saturated formation water to obtain nuclear magnetic resonance T2A spectrum;
step two: performing a constant-speed mercury-pressing experiment on the dry and compact reservoir sample to obtain mercury-pressing pore size distribution data;
step three: carrying out sectional processing on the mercury intrusion pore size distribution data;
step four: obtaining T using a cylindrical pore model2The conversion relation with the aperture size;
step five: using the data of step three and step four respectivelyConverting the relation to obtain T of each segment2Conversion method from aperture size to obtain complete T2A conversion method for spectral characterization of tight reservoir pore size distribution.
Further, the segmentation processing method in the third step is as follows: and according to the sequence from the small aperture to the large aperture, performing difference calculation on the aperture size data of different test points, and selecting the aperture size data with the aperture difference within an order of magnitude range as a data segment.
Further, the conversion relationship of the step four is as follows:
t expressed using KST model2Relationship between pore volume and pore fluid volume:
the volume of the lamellar fluid in the cylindrical pore model is:
the surface area of the cylindrical pores is:
St=2πrt·l (3)
the formula (2) and the formula (3) are brought into the formula (1) to obtain
Order to
Wherein the content of the first and second substances,T2representing transverse relaxation time, T2sDenotes the relaxation, V, of the particle surfacesDenotes the pore fluid volume, V denotes the pore volume, Vs tThe volume of the thin layer fluid in the cylindrical pore model is shown, h represents the thickness of the thin layer fluid, rtDenotes the pore radius,/. denotes the cylinder height, StExpressed as the surface area of the cylindrical pores, p2Is the relaxation rate; S/V is the pore specific surface; fsIs a form factor, dimensionless, to spherical pores, Fs3; for cylindrical pipes, Fs2. Relaxation rate ρ2Pore shape factor FsCan be approximately regarded as a constant, so C is a constant value;
by bringing (5), (6) and (7) into (4), T can be obtained2Conversion relationship with aperture size:
further, the C and h values of each segment can be obtained by the least square principle of each segment of data, so that the T of each segment is obtained2A method of conversion from pore size;
further, complete T is obtained2A conversion method for spectral characterization of tight reservoir pore size distribution.
Compared with the prior art, the invention has the following characteristics and advantages:
1. the invention fully considers the pore structure model of the rock, changes the common plane physical model for nuclear magnetic resonance into the cylindrical pore model consistent with mercury intrusion experiment, and improves the nuclear magnetic resonance data T2And accuracy of the correspondence between pore radii obtained from mercury intrusion data.
2. The method fully considers the characteristics of mercury intrusion data distribution, and reduces fitting errors caused by the data by using piecewise fitting.
Drawings
FIG. 1 is a sectional view of mercury intrusion pore size distribution data
FIG. 2 is a conversion process from a planar physical model to a cylindrical pore model
FIG. 3 shows T obtained by the conversion method of the present invention2Fitting relation with pore radius
FIG. 4 is a comparison result of the pore size distribution and mercury intrusion pore size distribution calculated by the conversion method of the present invention
FIG. 5 is a graph of T obtained empirically2Fitting relation with pore radius
FIG. 6 is a comparison of pore size distribution and mercury intrusion pore size distribution calculated by empirical method
FIG. 7 shows T obtained by a conventional method2Fitting relation with pore radius
FIG. 8 is a comparison of pore size distribution and mercury intrusion pore size distribution calculated by a conventional method
Detailed Description
The details of the present invention can be more clearly understood in conjunction with the accompanying drawings and the description of the embodiments of the present invention. However, the specific embodiments of the present invention described herein are for the purpose of illustration only and are not to be construed as limiting the invention in any way. Any possible variations based on the present invention may be conceived by the skilled person in the light of the teachings of the present invention, and these should be considered to fall within the scope of the present invention.
Example 1
Take four tight sandstone samples from the west of the slope of the Isanza, Ordos basin, China as an example; nuclear magnetic resonance and constant-rate mercury intrusion experiments were performed at the institute of hydrodynamics, the chinese academy of sciences.
The nuclear magnetic resonance experiment adopts a RecCore2500 type low-field nuclear magnetic resonance instrument, the resonance frequency is 2.38MHz, the number of echoes is 2048, the scanning frequency is 128, the waiting time is 5000ms, the echo interval is 0.6ms, and the gain is 50. Vacuumizing the rock sample at 120 ℃ for at least 24 hours, weighing, soaking in simulated formation water (total mineralization is 25000mg/L) for at least 24 hours, removing surface water of the rock sample by using micro-wet filter paper, weighing the rock sample and carrying out first nuclear magnetic resonance measurement to obtain T2Spectra.
The constant-speed mercury-pressing experiment adopts an ASPE-730 type constant-speed mercury-pressing experiment device produced by Coretest company, the contact angle is 140 degrees, and the surface tension is 485 dyne/cm. Preparing a rock sample into a cylindrical rock core with the diameter of 1cm and the length of 1cm, vacuumizing at 120 ℃ for at least 24 hours, soaking in a mercury solution, feeding mercury into the rock core at a constant speed of 0.00005ml/min, and finishing the test when the pressure reaches 6.2055 MPa.
Further, the fitting accuracy is affected by considering the distribution of the data points, because the more concentrated the data is, the better the fitting effect is, and the less distributed the data points is, the less dispersed the fitting effect is, the less desirable.
For compact sandstone, the large aperture is relatively small, and the small aperture is relatively small, and meanwhile, the pressure points of mercury intrusion test are intensively distributed in the high-pressure area corresponding to the small aperture, so that the data are reflected on mercury intrusion data, and the data points of the small aperture are more and more intensively distributed, while the data points of the large aperture are less and more dispersedly distributed.
According to the invention, the difference value calculation is carried out on the pore sizes of the mercury intrusion data of different test points according to the sequence from the small pore size to the large pore size, and the result shows that: the pore size differences increase in turn, indicating that the larger the pore size, the fewer data points per pore size range (as shown in FIG. 1, sample J-1 for example). Therefore, starting with the minimum aperture difference, the corresponding aperture range with an aperture difference within an order of magnitude range is selected as a data segment. The problem that sparse large aperture data are rarely considered in the fitting process of the data can be avoided to a certain extent by carrying out segmentation transformation on the aperture.
Next, T was obtained using a cylindrical pore model2A translation with aperture size. The specific process is as follows:
t expressed using KST model2Relationship between pore volume and pore fluid volume:
a planar physical model (volume of thin layer fluid is V) in a conventional method or an empirical methods) Replacement is by a cylindrical pore model (as shown in FIG. 2), the circleThe volume of the thin layer fluid in the cylindrical pores is:
the surface area of the cylindrical pores is:
St=2πrt·l (3)
the formula (2) and the formula (3) are brought into the formula (1) to obtain
Order to
By bringing (5), (6) and (7) into (4), T can be obtained2A conversion relationship with pore radius;
next, using the linear least squares principle, the C and h values are solved. The process of finding C and h values is based on the previous study by the inventors of the core T disclosed in the Effects of hole-through structure on gas permeability in the light and store resistance properties of the Upper Triassic Yanchang formation in the WesternOrdos Basin, China2A method for converting the spectrum transverse relaxation time into pore radius distribution.
Finally, the C and h values of different sections are obtained by adopting the same method for other sections of data, so that the complete rock core T is obtained2A method for converting spectrogram characterization pore size distribution. The calculation results are shown in Table 1The following steps:
t obtained by the above method2The fit to the pore radius is shown in FIG. 3 (sample T-6 as an example):
the comparison of the pore radius distribution of mercury intrusion and the pore radius distribution calculated by the above method is shown in FIG. 4 (sample T-6 is taken as an example)
Comparative example 1
Taking an empirical method as an example, namely:logarithms are taken at two sides, and the values of C and n are obtained by adopting the least square principle.
The calculation results are shown in table 2:
empirically obtained T2The fit to the pore radius is shown in FIG. 5 (sample T-6 as an example):
the comparison of the empirically calculated pore radius distribution with the mercury intrusion pore radius distribution is shown in FIG. 6 (sample T-6 is taken as an example):
comparative example 2
By way of example in a conventional manner, i.e. n is 1, T2=CrtAnd obtaining the C value by adopting a least square principle.
The calculation results are shown in table 3:
t obtained by conventional methods2The fit to the pore radius is shown in FIG. 7 (sample T-6 as an example):
the comparison of the pore radius distribution calculated by the conventional method with the mercury intrusion pore radius distribution is shown in FIG. 8 (taking sample T-6 as an example):
comparing tables 1, 2 and 3, and, at the same time, comparing FIGS. 3, 5 and 7, the results show that T is the conversion method obtained by the present invention2The value and the aperture fitting coefficient are both above 0.99, the fitting effect is better in intuition and is obviously better than that of a conventional method and an empirical method, and the method for predicting the aperture radius of the nuclear magnetic signal by utilizing the conversion relation obtained by the method is shown to be capable of better predicting the aperture distribution of the nuclear magnetic signal to a certain extent.
Further, comparing FIGS. 4, 6 and 8, the results show that the NMR T obtained by the conversion method of the present invention2Compared with the radius of the constant-speed mercury pressing pore throat, the pore diameter converted by the spectrum has better corresponding relation between the pore diameter and the constant-speed mercury pressing pore throat. Nuclear magnetic resonance T obtained by conventional method2Compared with the radius of the pore throat of the constant-speed mercury intrusion pore, the integral difference of the prediction result of the conventional formula method is larger. And the nuclear magnetic resonance T obtained by an empirical method2Compared with the radius of a constant-speed mercury-pressing pore throat, the pore diameter of spectrum conversion is compared with that of a constant-speed mercury-pressing pore throat, although the corresponding relation of a small pore diameter part and a large pore diameter part is good in a prediction result of an empirical formula method, the corresponding relation of the large pore diameter part is poor, the accuracy of the empirical formula method for large pore diameter prediction is low, and the side surface illustrates that the mercury-pressing data are processed in a segmented mode, so that the prediction accuracy is improved, and particularly the prediction for a large pore diameter section.
Therefore, the invention will T2The spectrum and aperture conversion relation uses a cylindrical model, and the mercury intrusion data is processed in a segmented manner, so that the nuclear magnetic resonance T with higher precision is obtained2A conversion method for spectral characterization of tight reservoir pore size distribution.
The present invention is not limited to the above embodiments, and in particular, various features described in different embodiments can be arbitrarily combined with each other to form other embodiments, and the features are understood to be applicable to any embodiment except the explicitly opposite descriptions, and are not limited to the described embodiments.
Claims (6)
1. Improve nuclear magnetic resonance T2The conversion method for the spectral characterization of the precision of the pore size distribution of the tight reservoir is characterized by comprising the following steps of:
the method comprises the following steps: performing nuclear magnetic resonance test on a compact reservoir sample of saturated formation water to obtain nuclear magnetic resonance T2A spectrum;
step two: carrying out mercury injection experiment on the dry compact reservoir sample to obtain mercury injection pore size distribution data;
step three: carrying out sectional processing on the mercury intrusion pore size distribution data;
step four: obtaining T using a cylindrical pore model2The conversion relation with the aperture size;
step five: respectively using the conversion relation in the step four to obtain the T of each section2Conversion method from aperture size to obtain complete T2A conversion method for spectral characterization of tight reservoir pore size distribution.
2. The method of claim 1 for increasing nuclear magnetic resonance T2The conversion method for the spectrum characterization of the pore size distribution precision of the tight reservoir is characterized in that the segmented processing method in the step three comprises the following steps: and according to the sequence from the small aperture to the large aperture, performing difference calculation on the aperture size data of different test points, and selecting the aperture size data with the aperture difference within an order of magnitude range as a data segment.
3. The method of claim 1 for increasing nuclear magnetic resonance T2The conversion method for the spectral characterization of the pore size distribution precision of the tight reservoir is characterized in that the conversion relation obtaining method in the fourth step is as follows:
t expressed using KST model2Relationship between pore volume and pore fluid volume:
the volume of the lamellar fluid in the cylindrical pore model is:
the surface area of the cylindrical pores is:
Si=2πrt·l (3)
the formula (2) and the formula (3) are brought into the formula (1) to obtain
Order to
Wherein, T2Representing transverse relaxation time, T2sDenotes the relaxation, V, of the particle surfacesDenotes the pore fluid volume, V denotes the pore volume, Vs iThe volume of the thin layer fluid in the cylindrical pore model is shown, h represents the thickness of the thin layer fluid, rtDenotes the pore radius,/. denotes the cylinder height, SiExpressed as the surface area of the cylindrical pores, p2Is the relaxation rate; S/V is the pore specific surface; fsIs a form factor, dimensionless, to spherical pores, Fs3; for cylindrical pipes, Fs2. Relaxation rate ρ2Pore shape factor FsCan be approximately regarded as a constant, so C is a constant value;
by bringing (5), (6) and (7) into (4), T can be obtained2Conversion relationship with pore radius:
4. the method of claim 3 for increasing nuclear magnetic resonance T2The conversion method for the spectral characterization of the precision of the pore size distribution of the tight reservoir is characterized in that the C and h values of each section are obtained by the data of each section through the least square principle, so that the T of each section is obtained2And the conversion method between the pore size and the pore size.
5. The method of claim 1 for increasing nuclear magnetic resonance T2The conversion method for the spectrum characterization of the pore size distribution precision of the tight reservoir is characterized in that the tight reservoir is one of tight sandstone, shale or mudstone.
6. The method of claim 1 for increasing nuclear magnetic resonance T2The conversion method for the spectral characterization of the pore size distribution precision of the compact reservoir is characterized in that the nuclear magnetic resonance in the step one is low-field nuclear magnetic resonance, and the mercury intrusion experiment in the step two is constant-speed mercury intrusion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010532575.3A CN111537544B (en) | 2020-06-08 | 2020-06-08 | Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010532575.3A CN111537544B (en) | 2020-06-08 | 2020-06-08 | Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111537544A true CN111537544A (en) | 2020-08-14 |
CN111537544B CN111537544B (en) | 2023-04-07 |
Family
ID=71978299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010532575.3A Active CN111537544B (en) | 2020-06-08 | 2020-06-08 | Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111537544B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112284999A (en) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | Sandstone pore size distribution determination method and application thereof |
CN113008927A (en) * | 2021-04-01 | 2021-06-22 | 中海石油(中国)有限公司 | Conversion method of nuclear magnetic resonance T2 spectrum and pore distribution of coal seam |
CN114544466A (en) * | 2022-04-27 | 2022-05-27 | 中国林业科学研究院木材工业研究所 | Method for measuring transverse relaxation rate of pore surface of wood cell cavity |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104634718A (en) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance |
CN106249306A (en) * | 2016-10-12 | 2016-12-21 | 贵州大学 | Shale pore structure detection method based on nuclear magnetic resonance, NMR |
CN108169099A (en) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance |
CN109916943A (en) * | 2019-04-25 | 2019-06-21 | 西安石油大学 | Rock core nuclear magnetic resonance T2Spectrogram relaxation time and pressure mercury pore throat radius conversion method |
CN110231272A (en) * | 2019-07-09 | 2019-09-13 | 中国地质大学(北京) | Tight sand aperture and nuclear magnetic resonance T2It is worth the determination method and system of transformational relation |
CN110618158A (en) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | Method for constructing capillary pressure curve of rock core by utilizing nuclear magnetic resonance information |
-
2020
- 2020-06-08 CN CN202010532575.3A patent/CN111537544B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104634718A (en) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance |
CN106249306A (en) * | 2016-10-12 | 2016-12-21 | 贵州大学 | Shale pore structure detection method based on nuclear magnetic resonance, NMR |
CN108169099A (en) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance |
CN109916943A (en) * | 2019-04-25 | 2019-06-21 | 西安石油大学 | Rock core nuclear magnetic resonance T2Spectrogram relaxation time and pressure mercury pore throat radius conversion method |
CN110231272A (en) * | 2019-07-09 | 2019-09-13 | 中国地质大学(北京) | Tight sand aperture and nuclear magnetic resonance T2It is worth the determination method and system of transformational relation |
CN110618158A (en) * | 2019-10-28 | 2019-12-27 | 中石化石油工程技术服务有限公司 | Method for constructing capillary pressure curve of rock core by utilizing nuclear magnetic resonance information |
Non-Patent Citations (1)
Title |
---|
HEXIN HUANG 等: "Effects of pore-throat structure on gas permeability in the tight sandstone reservoirs of the Upper Triassic Yanchang formation in the Western Ordos Basin, China", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112284999A (en) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | Sandstone pore size distribution determination method and application thereof |
CN113008927A (en) * | 2021-04-01 | 2021-06-22 | 中海石油(中国)有限公司 | Conversion method of nuclear magnetic resonance T2 spectrum and pore distribution of coal seam |
CN114544466A (en) * | 2022-04-27 | 2022-05-27 | 中国林业科学研究院木材工业研究所 | Method for measuring transverse relaxation rate of pore surface of wood cell cavity |
CN114544466B (en) * | 2022-04-27 | 2022-07-05 | 中国林业科学研究院木材工业研究所 | Method for measuring transverse relaxation rate of pore surface of wood cell cavity |
Also Published As
Publication number | Publication date |
---|---|
CN111537544B (en) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111537544B (en) | Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir | |
Liu et al. | Study on the seepage characteristics of coal based on the Kozeny-Carman equation and nuclear magnetic resonance experiment | |
Jian-Ping et al. | A QUANTITATIVE EVALUATION METHOD OF LOW PERMEABLE SANDSTONE PORE STRUCTURE BASED ON NUCLEAR MAGNETIC RESONANCE (NMR) LOGGING: A CASE STUDY OF Es4 FORMATION IN THE SOUTH SLOPE OF DONGYING SAG | |
CN110133035B (en) | Method for quantitatively representing pore throat structure of sandstone reservoir by coupling nuclear magnetic resonance with constant-speed mercury intrusion | |
CN105866002B (en) | A kind of accurate dice NMR porosity test method | |
CN108694264B (en) | Method for determining permeability of shale gas reservoir | |
CN110595953B (en) | Experimental test device and method for shale mixing wettability | |
CN109856176B (en) | NMR and LTNA-based compact reservoir full-size pore quantitative characterization method | |
CN113075102B (en) | Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time | |
CN105910941B (en) | The test method of content of unfrozen water in frozen earth based on pressure plate apparatus | |
CN110108616A (en) | Coal sample centrifugal process loss of signal compensation method in a kind of pair of nuclear magnetic resonance experiment | |
CN115078210B (en) | Shale pore structure testing method | |
Wang et al. | Multifractal characteristics of shale and tight sandstone pore structures with nitrogen adsorption and nuclear magnetic resonance | |
CN104142290A (en) | Rock sample bulk volume and variable density determination apparatus and method thereof | |
CN106897531B (en) | Quantitative evaluation method for permeability of low-permeability limestone reservoir | |
Li et al. | The description of shale reservoir pore structure based on method of moments estimation | |
CN110309611B (en) | Gas-water two-phase seepage rule prediction method and system based on gas-water thickness distribution | |
Li et al. | Controlling effects of pore‐throat structure and fractal characteristics on the physical properties of ultra‐low permeability sandstone reservoirs: A case study of the sixth member of the Yanchang Formation in the Xiaojiahe area, Ordos Basin | |
Teng et al. | Pore structure characterization based on NMR experiment: A case from the Shanxi Formation tight sandstones in the Daning-Jixian area, eastern Ordos Basin | |
CN111579584A (en) | Method for identifying gas enrichment area in coal based on infrared thermal imaging technology | |
CN113984590B (en) | Method for calculating space tortuosity and gas diffusion coefficient distribution of heterogeneous rock | |
CN108846218A (en) | Predict the method and system of the inner ratio surface area of shale reservoir different type hole | |
CN113984589B (en) | Method for calculating rock tortuosity and gas diffusion coefficient | |
Song et al. | Pore fractal characteristics of lignite at different temperatures based on mercury intrusion test | |
Song et al. | Estimation of oil saturation via pseudo capillary pressure curve from nuclear magnetic resonance log data in tight conglomerate reservoirs |
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 |