CN108169099A - A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance - Google Patents

A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance Download PDF

Info

Publication number
CN108169099A
CN108169099A CN201810042413.4A CN201810042413A CN108169099A CN 108169099 A CN108169099 A CN 108169099A CN 201810042413 A CN201810042413 A CN 201810042413A CN 108169099 A CN108169099 A CN 108169099A
Authority
CN
China
Prior art keywords
aperture
magnetic resonance
nuclear magnetic
gas reservoir
shale
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
Application number
CN201810042413.4A
Other languages
Chinese (zh)
Inventor
王亮
贾玉泽
杨琴琴
黄林林
谢韦峰
朱明月
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201810042413.4A priority Critical patent/CN108169099A/en
Publication of CN108169099A publication Critical patent/CN108169099A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Dispersion Chemistry (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation methods based on nuclear magnetic resonance, include the following steps:Rock core is collected;Duplicate Samples are drilled through, are carried out from oil suction, the measurement of water nuclear magnetic resonance experiment;Comparative analysis Duplicate Samples oil, water nuclear magnetic resonance T 2 spectrum difference, determine distribution of the different wetting porosity type on nuclear magnetic resonance T 2 spectrum;Shale gas reservoir full aperture distribution curve is obtained according to high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption;Further, aperture and corresponding T2 times crossplot are obtained;According to different aperture type aperture and the crossplot of corresponding T2 values, point porosity type establishes the quantitative calculation in aperture.The advantage of the invention is that:The technology can quantify and calculate shale gas reservoir pore space full aperture distribution curve;Meanwhile nuclear-magnetism measurement is quick, simple, lossless, it is stronger compared to high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption practicability;Compared to conventional method, result of calculation is more accurate.

Description

A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance
Technical field
The present invention relates to shale gas RESERVOIR PORE STRUCTURE assessment technique field, more particularly to a kind of page based on nuclear magnetic resonance Rock gas reservoir pore structure quantitative calculation method.
Background technology
Shale gas refers mainly to the natural gas in shale matrix pores or crack with ADSORPTION STATE or free state preservation.Shale gas It is the unconventional petroleum resources that another is important after compact sandstone gas and coal bed gas.At present, shale gas reservoir pore space is evaluated The experimental method of structure mainly includes:Scanning electron microscope, high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption, nuclear magnetic resonance etc. are real It tests;
(1) scanning electron microscope, pressure mercury, nitrogen adsorption, carbon dioxide adsorption
Scanning electron microscope combination argon ion polishing technology can carry out qualitative observation, but be difficult to quantitatively characterizing shale to pore structure The pore structure characteristic of gas reservoir.High-pressure mercury, nitrogen adsorption, carbon dioxide adsorption experiment can obtain quantitative response reservoir pore space The pore size distribution curve of structure.However, because of the difference of its experimental principle, experiment condition, mercury injection method mainly characterizes>The hole of 50nm Diameter, nitrogen adsorption methods mainly characterize the aperture of 2~50nm, and carbon dioxide adsorption method mainly characterizes<The aperture of 2nm;Therefore, not It can obtain the pore size distribution curve of objective characterisation reservoir micropore structure.
(2) nuclear magnetic resonance evaluation pore structure
Nuclear magnetic resonance logging data evaluation RESERVOIR PORE STRUCTURE has its unique advantage;Nuclear magnetic resonance evaluates conventional sandstone or carbonic acid The Method And Principle of rock salt RESERVOIR PORE STRUCTURE is as follows:
By NMR relaxation mechanism it is found that the lateral relaxation time T2 observed can be expressed as
In formula:Volume (freedom) relaxation times of the T2B for fluid, ms;D is diffusion coefficient, μm2/ms;G is magnetic field gradient, gauss/cm;TE is echo sounding, ms;R is pore radius, um;ρ is the lateral surfaces relaxation strength of rock, μm/s, by hole The influence of surface wettability, i.e., it is related with porosity type.Fs is the pore shape factor, for spherical hole, Fs=3;And For column pipeline, Fs=2.
For diffusion relaxation item (Section 3 on the right of equation), when using smaller echo sounding (TE=0.06ms) and uniformly During magnetic field observation, then either saturated water or saturated oils condition, diffusion relaxation item can be ignored.Above formula becomes:
Normal conditions, fluid volume relaxation T2B are more much bigger (T2B than the T2 relaxation times>>T2).Therefore, equal sign in formula 1 The right first item Volume Relaxation part can ignore to get to:
R=ρ Fs×T2 (4)
R=C × T2 (5)
In formula:C=ρ Fs.The key of nuclear-magnetism evaluation reservoir pores r is the accurately relationship of determining C or determining formula (5) Formula;
At present, mainly include three classes using the method for nuclear magnetic resonance characterization reservoir micropore structure:Linear function scale Method, J functions and SDR models couplings method, classification subsection power function method.He Yudan et al. is it was verified that linear function scale method structure The capillary pressure curve built and experiment capillary pressure curve are larger in aperture larynx fractional error, and effect is undesirable;Classification segmentation power Function method overcomes linear function scale method aperture throat point effect is undesirable the shortcomings that;But this method needs known store up The porosity of layer, permeability could build capillary pressure curve after classifying to reservoir;Since shale gas reservoir properties are poor, The dyscalculia of porosity and permeability, similar to classification subsection power function method, J functions and SDR models coupling methods fail to solve The determining problem of J function category standards so that the utilization of the method is also restrained.However the above method is for sandstone or carbon The conventional reservoirs such as hydrochlorate, it is less to the research of shale gas reservoir micropore structure nuclear-magnetism characterization.
For the evaluation of shale gas RESERVOIR PORE STRUCTURE, Li Jun et al. in periodical《Logging technique》On delivered an article 《Shale gas rock core nuclear magnetic resonance T2 and aperture size quantitative relationship》, it is believed that C values are 47 μm/s.But for shale gas reservoir Complex Rock component, porosity type, the value not should be definite value, i.e., do not consider that surface is relaxed in shale different wetting hole component type The difference of Henan rate.
Invention content
The present invention in view of the drawbacks of the prior art, provides a kind of shale gas RESERVOIR PORE STRUCTURE based on nuclear magnetic resonance and determines Computational methods are measured, can effectively solve the problem that the above-mentioned problems of the prior art.
In order to realize more than goal of the invention, the technical solution that the present invention takes is as follows:
A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance, includes the following steps:
Step 1, rock core is collected;
Step 2, Duplicate Samples are drilled through, are carried out from oil suction, the measurement of water nuclear magnetic resonance experiment;
Drill through cylindric shale oil and gas reservoir rock core;The rock core is cut from centre, it is smooth to be prepared as end face, similar length One group of Duplicate Samples, two pieces of core numbers are respectively a, b;Rock core a self-primings brine, another piece of rock core b are selected from oil suction;According to rock Sample nuclear magnetic resonance parameter tests specifications of surveys, carries out Nuclear Magnetic Resonance Measurement to two blocks of rock cores, obtains the nuclear magnetic resonance of two blocks of rock cores T2 is composed;
Step 3, comparative analysis Duplicate Samples oil, water nuclear-magnetism T2 resonance spectrum differences, determine different wetting porosity type in core Distribution in magnetic resonance T2 spectrums;
Step 4, shale gas reservoir full aperture distribution curve is obtained according to high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption;
Step 5, aperture and corresponding nuclear magnetic resonance T2 time crossplots are obtained;
It is descending by aperture by shale full aperture distribution curve, the corresponding saturation degree in each aperture is added up, is obtained To accumulative saturation degree and the change curve in aperture;It is descending by T2 values to full water nuclear magnetic resonance T 2 spectrum, every T2 values are corresponded to Saturation degree added up, obtain accumulation saturation degree and T2 change curve.
It is bent in the change curve in accumulative saturation degree and aperture and the variation of accumulation saturation degree and T2 using equisaturation principle The T2 values corresponding to the aperture under same accumulation saturation degree are found in line;According to the corresponding T2 values in aperture, cross plot is done;Base In step 3, the porosity type carried out to different T2 values divides, further by porosity type difference, by the friendship in cross plot It can put and be divided into three parts;Different aperture type (organic hole, inorganic hole, microcrack) is presented different between T2 times and aperture Variation tendency, the difference of variation tendency may be due to the shape of hole corresponding to different aperture type and pore surface relaxation Caused by the difference of Henan rate.
Step 6, according in step 5, the crossplot of different aperture type aperture and corresponding T2 times, a point porosity type are built The quantitative calculation in vertical aperture;
To different aperture type aperture and T2 values, using exponential function form, aperture and the relational expression of T2 are fitted, has been shown in Table 1。
Table 1
Based on nuclear magnetic resonance T 2 spectrum, according to the formula in table 1, by porosity type difference, different T2 values ranges is used Different calculation formula, is calculated aperture in table 1.
Preferably, cylindric shale oil and gas reservoir core diameter is 2.54cm or 3.81cm in step 2, length is more than 6cm。
Preferably, rock core a self-primings brine is the Nacl solution of 40000ppm in step 2.
Preferably, in step 2 rock core b from oil suction be dodecane.
Preferably, rock core self-priming time is 48 hours in step 2.
Further, step 4 is as follows:1) shale is analyzed using high-pressure mercury method, with reference to Washburn equations obtain the pore volume in each aperture in the first pore diameter range;2) using nitrogen adsorption methods to the first powdery shale It is analyzed, the pore volume in each aperture in the second pore diameter range is obtained with reference to BJH models;3) using carbon dioxide adsorption method pair Second powdery shale is analyzed, and the pore volume in each aperture in third pore diameter range is obtained with reference to DFT models;4) such as the first hole Model is overlapped between diameter range, the second pore diameter range and third pore diameter range, each hole in overlapping range is obtained using weighted mean method The pore volume of diameter.Obtain shale full aperture distribution curve.
Compared with prior art the advantage of the invention is that:
(1) conventional means of the characterization such as high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption shale gas pore-size distribution and side There is limitation in method, i.e., mercury injection method mainly characterizes in the range of aperture characterization>The aperture of 50nm, the main table of nitrogen adsorption methods The aperture of 2~50nm is levied, carbon dioxide adsorption method mainly characterizes<The aperture of 2nm;The technology can quantify and calculate shale gas reservoir Hole full aperture distribution curve;Meanwhile nuclear-magnetism measurement is quick, simple, lossless, compared to high-pressure mercury, nitrogen adsorption, titanium dioxide Carbon adsorption practicability is stronger.
(2) technology considers the shadow that different aperture type and its wetting sex differernce in shale characterize nuclear-magnetism pore-size distribution It rings, compared to conventional method, result of calculation is more accurate.
Description of the drawings
Fig. 1 is the nuclear magnetic resonance T 2 spectrum figure of two blocks of rock cores of the embodiment of the present invention;
Fig. 2 is shale full aperture distribution map of the embodiment of the present invention;
Fig. 3 is the reversed accumulation curve in aperture of shale samples of embodiment of the present invention H10;
Fig. 4 is the T2 distribution reversal accumulation curves of shale samples of embodiment of the present invention H10;
Fig. 5 is shale samples of embodiment of the present invention different pore size and corresponding T2 times cross plot;
Fig. 6 is predicted for LM1 of embodiment of the present invention nuclear magnetic resonance with surveying pore-size distribution comparison diagram;
Fig. 7 is predicted for LM2 of embodiment of the present invention nuclear magnetic resonance with surveying pore-size distribution comparison diagram;
Fig. 8 is predicted for LM3 of embodiment of the present invention nuclear magnetic resonance with surveying pore-size distribution comparison diagram;
Fig. 9 is predicted for LM4 of embodiment of the present invention nuclear magnetic resonance with surveying pore-size distribution comparison diagram;
Figure 10 is predicted for H10 of embodiment of the present invention nuclear magnetic resonance with surveying pore-size distribution comparison diagram.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention more comprehensible, develop simultaneously embodiment referring to the drawings, right The present invention is described in further details.
A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance, includes the following steps:
Step 1, rock core is collected;
Step 2, Duplicate Samples are drilled through, are carried out from oil suction, the measurement of water nuclear magnetic resonance experiment;
Drill through a diameter of 2.54cm (or 3.81cm), length is more than the cylindric shale oil and gas reservoir rock cores of 6cm;By the rock core It is cut from centre, is prepared as that end face is smooth, one group of Duplicate Samples of similar length (two pieces of core numbers are respectively a, b);Select rock Heart a self-primings brine (the Nacl solution of 40000ppm), another piece of rock core b are from oil suction (dodecane);Rock core self-priming time is small for 48 When.According to rock sample nuclear magnetic resonance parameter experiment specifications of surveys (SY/T 6490-2007), nuclear magnetic resonance survey is carried out to two blocks of rock cores Amount, obtains the nuclear magnetic resonance T 2 spectrum of two blocks of rock cores, as shown in Figure 1.
Step 3, comparative analysis Duplicate Samples oil, water nuclear magnetic resonance T 2 spectrum difference, determine different wetting porosity type in core Distribution in magnetic resonance T2 spectrums;
Duplicate Samples rock core is compared from oil suction, the difference of Absorb Water nuclear magnetic resonance T 2 spectrum, thinks that rock core nuclear-magnetism T2 spectrums exist accordingly Spectral peak (P1) between 0.01~0.3ms is mainly related with the response of proton in organic matter, that is, corresponds to organic hole (glossy wet); Spectral peak (P2) between 0.3~10ms is mainly by (clay mineral intergranular pore and intracrystalline pore, pyrite intracrystalline pore) hydrogen in inorganic hole The response of core causes, that is, corresponds to inorganic hole (water wetting);>The spectral peak (P3) of 10ms corresponds to microcrack (mixing wetting).
Step 4, shale gas reservoir full aperture distribution curve is obtained according to high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption;
It is as follows:1) shale is analyzed using high-pressure mercury method, first is obtained with reference to washburn equations The pore volume in each aperture in pore diameter range;2) the first powdery shale is analyzed using nitrogen adsorption methods, with reference to BJH models Obtain the pore volume in each aperture in the second pore diameter range;3) the second powdery shale is divided using carbon dioxide adsorption method Analysis, the pore volume in each aperture in third pore diameter range is obtained with reference to DFT models;4) such as the first pore diameter range, the second aperture model It encloses and is overlapped model between third pore diameter range, the pore volume in each aperture in overlapping range is obtained using weighted mean method.Foundation Step and method shown in the patent obtain shale full aperture distribution curve, as shown in Figure 2.
Step 5, aperture and corresponding nuclear magnetic resonance T2 time crossplots are obtained
For petrophysics angle, what Fig. 2 pore size distribution curves reflected is the hole under a certain pore size control The distribution of volume.What nuclear magnetic resonance T2 distributions equally reflected is the porosity component under a certain aperture.Pore size distribution curve and T2 The meaning that distribution curve is characterized is the same, but since the cloth points of pore size distribution curve and T2 distribution curves have differences, because This is difficult to directly compare by form, and corresponding T2 values are obtained from pore size distribution curve.
For this purpose, by shale full aperture distribution curve in Fig. 2, to each aperture corresponding saturation degree descending by aperture Added up, obtain the change curve in accumulative saturation degree and aperture, such as Fig. 3;To in Fig. 1 full water nuclear-magnetism T2 compose, by T2 values by It arrives greatly small, the corresponding saturation degree of every T2 values is added up, obtain accumulation saturation degree and the change curve of T2, such as Fig. 4.
Using " equisaturation " principle, the T2 corresponding to the aperture under same accumulation saturation degree is found in Fig. 3 and Fig. 4 Value.According to the corresponding T2 values in aperture, do cross plot and obtain Fig. 5 based on the hole in step 3, carried out to different T2 values Plotted point in cross plot 5 further by porosity type difference, is divided into three parts by Type division.Different aperture class in Fig. 5 Different variation tendencies are presented in type (organic hole, inorganic hole, microcrack) between T2 times and aperture, the difference of variation tendency may It is caused by the shape of hole corresponding to different aperture type and the difference of pore surface relaxation rate.
Step 6, according in step 5, the crossplot of different aperture type aperture and corresponding T2 values, a point porosity type is established The quantitative calculation in aperture
To different aperture type aperture in Fig. 5 and T2 values, using exponential function form, it is fitted aperture and the relationship of T2 Formula is shown in Table 1.
1 shale gas reservoir pores transformation model of table
Based on nuclear magnetic resonance T 2 spectrum, according to the formula in table 1, i.e., by porosity type difference, different T2 value ranges are adopted With different calculation formula, aperture can be calculated.Fig. 6~Figure 10 shows that the pore-size distribution of nuclear-magnetism prediction and experiment measure aperture The identical property of distribution is preferable;The shale gas RESERVOIR PORE STRUCTURE quantitative calculation method that this patent proposes is reliable.
Those of ordinary skill in the art will understand that the embodiments described herein, which is to help reader, understands this hair Bright implementation, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.Ability The those of ordinary skill in domain can be made according to these technical inspirations disclosed by the invention it is various do not depart from essence of the invention its Its various specific deformation and combination, these deformations and combination are still within the scope of the present invention.

Claims (6)

1. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance, it is characterised in that including following step Suddenly:
Step 1, rock core is collected;
Step 2, Duplicate Samples are drilled through, are carried out from oil suction, the measurement of water nuclear magnetic resonance experiment;
Drill through cylindric shale oil and gas reservoir rock core;The rock core from centre is cut, is prepared as that end face is smooth, the one of similar length Group Duplicate Samples, two pieces of core numbers are respectively a, b;Rock core a self-primings brine, another piece of rock core b are selected from oil suction;According to rock sample core Magnetic resonance parameters test specifications of surveys, carry out Nuclear Magnetic Resonance Measurement to two blocks of rock cores, obtain the nuclear magnetic resonance T 2 spectrum of two blocks of rock cores;
Step 3, comparative analysis Duplicate Samples oil, water nuclear magnetic resonance T 2 spectrum difference, determine that different wetting porosity type is total in nuclear-magnetism The distribution shaken in T2 spectrums;
Step 4, shale gas reservoir full aperture distribution curve is obtained according to high-pressure mercury, nitrogen adsorption, carbon dioxide adsorption;
Step 5, aperture and corresponding nuclear magnetic resonance T2 time crossplots are obtained;
It is descending by aperture by shale full aperture distribution curve, the corresponding saturation degree in each aperture is added up, is tired out Count the change curve in saturation degree and aperture;It is descending by T2 values to full water nuclear magnetic resonance T 2 spectrum, it is corresponding to every T2 values full Added up with degree, obtain accumulation saturation degree and the change curve of T2;
Using equisaturation principle, in accumulative saturation degree and the change curve in aperture and the change curve of accumulation saturation degree and T2 Find the T2 values corresponding to the aperture under same accumulation saturation degree;According to the corresponding T2 values in aperture, cross plot is done;Based on step In rapid 3, the porosity type carried out to different T2 values divides, further by porosity type difference, by the plotted point in cross plot It is divided into three parts;Organic hole, inorganic hole and microcrack are included for porosity type, presented between T2 times and aperture different Variation tendency, the difference of variation tendency may be due to the shape of hole corresponding to different aperture type and pore surface relaxation Caused by the difference of rate;
Step 6, according in step 5, the crossplot of different aperture type aperture and corresponding T2 times, a point porosity type establish hole The quantitative calculation of diameter;
To different aperture type aperture and T2 values, using exponential function form, aperture and the relational expression of T2 are fitted, has been shown in Table 1;
Table 1
Based on nuclear magnetic resonance T 2 spectrum, according to the formula in table 1, by porosity type difference, to different T2 values ranges using in table 1 Different calculation formula, is calculated aperture.
2. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance according to claim 1, It is characterized in that:Cylindric shale oil and gas reservoir core diameter is 2.54cm or 3.81cm in step 2, length is more than 6cm.
3. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance according to claim 1, It is characterized in that:Rock core a self-primings brine is the Nacl solution of 40000ppm in step 2.
4. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance according to claim 1, It is characterized in that:In step 2 rock core b from oil suction be dodecane.
5. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance according to claim 1, It is characterized in that:Rock core self-priming time is 48 hours in step 2.
6. a kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance according to claim 1, It is characterized in that:Step 4 is as follows:1) shale is analyzed using high-pressure mercury method, with reference to washburn equations Obtain the pore volume in each aperture in the first pore diameter range;2) the first powdery shale is analyzed using nitrogen adsorption methods, tied It closes BJH models and obtains the pore volume in each aperture in the second pore diameter range;3) using carbon dioxide adsorption method to the second powdery page Rock is analyzed, and the pore volume in each aperture in third pore diameter range is obtained with reference to DFT models;4) such as the first pore diameter range, the Model is overlapped between two pore diameter ranges and third pore diameter range, the pore-body in each aperture in overlapping range is obtained using weighted mean method Product.Obtain shale full aperture distribution curve.
CN201810042413.4A 2018-01-17 2018-01-17 A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance Pending CN108169099A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810042413.4A CN108169099A (en) 2018-01-17 2018-01-17 A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810042413.4A CN108169099A (en) 2018-01-17 2018-01-17 A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance

Publications (1)

Publication Number Publication Date
CN108169099A true CN108169099A (en) 2018-06-15

Family

ID=62514990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810042413.4A Pending CN108169099A (en) 2018-01-17 2018-01-17 A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance

Country Status (1)

Country Link
CN (1) CN108169099A (en)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109085108A (en) * 2018-09-14 2018-12-25 重庆科技学院 A method of evaluation slippery water intrusion volume is to shale permeability impact effect
CN109283114A (en) * 2018-11-07 2019-01-29 西安石油大学 A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat
CN109520894A (en) * 2018-12-29 2019-03-26 西南石油大学 A kind of shale organic matter is forced from pipette prediction technique
CN109856176A (en) * 2019-03-27 2019-06-07 中国石油大学(华东) The full-scale hole quantitatively characterizing method of compact reservoir based on NMR and LTNA
CN109887614A (en) * 2019-01-22 2019-06-14 中国石油天然气股份有限公司 Hydraulic fracture analysis method and device
CN109884109A (en) * 2019-03-22 2019-06-14 中国石油大学(华东) A kind of shale hole water adsorption amount, free amount and its evaluation method of distribution
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
CN110410058A (en) * 2019-06-20 2019-11-05 中石化石油工程技术服务有限公司 A method of correction core experiment result scale two dimension nuclear magnetic resonance log
CN110687612A (en) * 2019-09-17 2020-01-14 中国石油天然气股份有限公司 Shale oil analysis method and device for continuously representing content of adsorbed oil and free oil
CN111398122A (en) * 2020-04-03 2020-07-10 中国矿业大学 Comprehensive characterization method for heterogeneity characteristics of full-scale pore structure of shale
CN111537544A (en) * 2020-06-08 2020-08-14 长安大学 Improve nuclear magnetic resonance T2Conversion method for spectral characterization of dense reservoir pore size distribution precision
CN111624146A (en) * 2019-02-28 2020-09-04 中国石油化工股份有限公司 Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir
CN112198185A (en) * 2019-07-08 2021-01-08 中国石油化工股份有限公司 Method for quantitatively characterizing shale oil occurrence characteristics
CN112284999A (en) * 2020-10-26 2021-01-29 中国石油大学(华东) Sandstone pore size distribution determination method and application thereof
CN112304999A (en) * 2020-09-16 2021-02-02 宜宾学院 Quantitative statistics method for pore characteristics of scanning electron microscope of shale micro-nano pores
CN112577881A (en) * 2020-12-31 2021-03-30 煤炭科学技术研究院有限公司 Low-permeability coal seam multi-scale micro-pore structure fine characterization method
CN112782477A (en) * 2019-11-11 2021-05-11 中国石油化工股份有限公司 Method and system for measuring electrical response characteristics of rock core in different wetting states
CN112946005A (en) * 2021-02-02 2021-06-11 中国石油大学(华东) Shale microcrack evaluation 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
CN113075102A (en) * 2021-03-11 2021-07-06 中国地质大学(北京) Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time
CN113405962A (en) * 2020-03-16 2021-09-17 中国石油天然气股份有限公司 Storage space evaluation method and device for carbonate reservoir
CN113533156A (en) * 2021-06-30 2021-10-22 西安石油大学 Identification method for microscopic pore structure characteristics and multi-type pore fluid of shale oil reservoir
CN113850689A (en) * 2021-09-24 2021-12-28 西南石油大学 Pore structure joint characterization method based on least square method and data processing method
CN115078210A (en) * 2022-05-05 2022-09-20 中国科学院地质与地球物理研究所 Shale pore structure testing method
CN115406921A (en) * 2022-08-27 2022-11-29 西安石油大学 Method for determining optimal degassing temperature of low-maturity hydrocarbon source rock
CN115753866A (en) * 2022-11-17 2023-03-07 重庆科技学院 Method for quantitatively characterizing shale pore structure
CN117110349A (en) * 2023-06-28 2023-11-24 东北石油大学 Quantitative evaluation method for Dan Duo-scale reservoir space wettability of shale oil reservoir rock

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141637A (en) * 2010-01-28 2011-08-03 中国石油天然气股份有限公司 Method for continuously quantitative evaluation of pore structures of reservoir strata by utilizing nuclear magnetic resonance well logging data
CN104237284A (en) * 2014-10-17 2014-12-24 西安石油大学 Nuclear magnetic resonance detection method for brittle and hard shale microcrack damage
CN104634718A (en) * 2015-03-05 2015-05-20 中国石油大学(华东) Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance
CN104697915A (en) * 2015-03-20 2015-06-10 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Shale micropore size and fluid distribution analysis method
CN104819923A (en) * 2015-05-17 2015-08-05 西南石油大学 Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance
CN104990851A (en) * 2015-06-23 2015-10-21 西南石油大学 Novel experimental research method for shale sensitivity
CN105241913A (en) * 2015-10-10 2016-01-13 西安石油大学 Nuclear magnetic resonance quantitative analysis method for rock micro-crack damage variable
CN105445161A (en) * 2015-11-16 2016-03-30 中国石油大学(北京) Characterizing method for pore volumes of full apertures of shale
CN105466830A (en) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 A pore-throat dimension distribution identification method for reservoir sandstone
CN106053315A (en) * 2016-06-06 2016-10-26 中国石油大学(华东) Method for classifying bioclastic limestone reservoir pore structures

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141637A (en) * 2010-01-28 2011-08-03 中国石油天然气股份有限公司 Method for continuously quantitative evaluation of pore structures of reservoir strata by utilizing nuclear magnetic resonance well logging data
CN104237284A (en) * 2014-10-17 2014-12-24 西安石油大学 Nuclear magnetic resonance detection method for brittle and hard shale microcrack damage
CN104634718A (en) * 2015-03-05 2015-05-20 中国石油大学(华东) Calibration method for representing dense sandstone pore size distribution by adopting nuclear magnetic resonance
CN104697915A (en) * 2015-03-20 2015-06-10 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Shale micropore size and fluid distribution analysis method
CN104819923A (en) * 2015-05-17 2015-08-05 西南石油大学 Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance
CN104990851A (en) * 2015-06-23 2015-10-21 西南石油大学 Novel experimental research method for shale sensitivity
CN105241913A (en) * 2015-10-10 2016-01-13 西安石油大学 Nuclear magnetic resonance quantitative analysis method for rock micro-crack damage variable
CN105445161A (en) * 2015-11-16 2016-03-30 中国石油大学(北京) Characterizing method for pore volumes of full apertures of shale
CN105466830A (en) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 A pore-throat dimension distribution identification method for reservoir sandstone
CN106053315A (en) * 2016-06-06 2016-10-26 中国石油大学(华东) Method for classifying bioclastic limestone reservoir pore structures

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
李军等: "《利用核磁共振技术确定有机孔与无机孔孔径分布—以四川盆地涪陵地区志留系龙马溪组页岩气储层为例》", 《石油与天然气地质》 *
李军等: "《页岩气岩心核磁共振T2与孔径尺寸定量关系》", 《测井技术》 *
李海波等: "《核磁共振T2谱换算孔隙半径分布方法研究》", 《波谱学杂志》 *
王凯等: "《基于二氧化碳和液氮吸附、高压压汞和低场核磁共振的煤岩多尺度孔径表征》", 《中国矿业》 *
肖佃师等: "《联合核磁共振和恒速压汞方法测定致密砂岩孔喉结构》", 《石油勘探与开发》 *

Cited By (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109085108A (en) * 2018-09-14 2018-12-25 重庆科技学院 A method of evaluation slippery water intrusion volume is to shale permeability impact effect
CN109283114A (en) * 2018-11-07 2019-01-29 西安石油大学 A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat
CN109520894A (en) * 2018-12-29 2019-03-26 西南石油大学 A kind of shale organic matter is forced from pipette prediction technique
CN109887614A (en) * 2019-01-22 2019-06-14 中国石油天然气股份有限公司 Hydraulic fracture analysis method and device
CN111624146A (en) * 2019-02-28 2020-09-04 中国石油化工股份有限公司 Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir
CN109884109B (en) * 2019-03-22 2022-08-02 中国石油大学(华东) Evaluation method for shale pore water adsorption amount, free amount and distribution of shale pore water adsorption amount and free amount
CN109884109A (en) * 2019-03-22 2019-06-14 中国石油大学(华东) A kind of shale hole water adsorption amount, free amount and its evaluation method of distribution
CN109856176B (en) * 2019-03-27 2022-03-08 中国石油大学(华东) NMR and LTNA-based compact reservoir full-size pore quantitative characterization method
CN109856176A (en) * 2019-03-27 2019-06-07 中国石油大学(华东) The full-scale hole quantitatively characterizing method of compact reservoir based on NMR and LTNA
CN110410058A (en) * 2019-06-20 2019-11-05 中石化石油工程技术服务有限公司 A method of correction core experiment result scale two dimension nuclear magnetic resonance log
CN112198185A (en) * 2019-07-08 2021-01-08 中国石油化工股份有限公司 Method for quantitatively characterizing shale oil occurrence characteristics
CN112198185B (en) * 2019-07-08 2023-11-28 中国石油化工股份有限公司 Method for quantitatively representing occurrence characteristics of shale oil
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
US11086044B2 (en) 2019-09-17 2021-08-10 Petrochina Company Limited Shale oil analysis method and apparatus for continuously characterizing saturation of adsorbed oil and free oil
CN110687612B (en) * 2019-09-17 2020-09-08 中国石油天然气股份有限公司 Shale oil analysis method and device for continuously representing content of adsorbed oil and free oil
CN110687612A (en) * 2019-09-17 2020-01-14 中国石油天然气股份有限公司 Shale oil analysis method and device for continuously representing content of adsorbed oil and free oil
CN112782477A (en) * 2019-11-11 2021-05-11 中国石油化工股份有限公司 Method and system for measuring electrical response characteristics of rock core in different wetting states
CN112782477B (en) * 2019-11-11 2024-05-14 中国石油化工股份有限公司 Method and system for measuring electrical response characteristics of rock core in different wetting states
CN113405962B (en) * 2020-03-16 2023-02-10 中国石油天然气股份有限公司 Storage space evaluation method and device for carbonate reservoir
CN113405962A (en) * 2020-03-16 2021-09-17 中国石油天然气股份有限公司 Storage space evaluation method and device for carbonate reservoir
CN111398122A (en) * 2020-04-03 2020-07-10 中国矿业大学 Comprehensive characterization method for heterogeneity characteristics of full-scale pore structure of shale
CN111537544A (en) * 2020-06-08 2020-08-14 长安大学 Improve nuclear magnetic resonance T2Conversion method for spectral characterization of dense reservoir pore size distribution precision
CN111537544B (en) * 2020-06-08 2023-04-07 长安大学 Improve nuclear magnetic resonance T 2 Conversion method for spectral characterization of pore size distribution precision of tight reservoir
CN112304999B (en) * 2020-09-16 2023-09-15 宜宾学院 Quantitative statistical method for shale micro-nano pore characteristics of scanning electron microscope
CN112304999A (en) * 2020-09-16 2021-02-02 宜宾学院 Quantitative statistics method for pore characteristics of scanning electron microscope of shale micro-nano pores
CN112284999A (en) * 2020-10-26 2021-01-29 中国石油大学(华东) Sandstone pore size distribution determination method and application thereof
CN112577881A (en) * 2020-12-31 2021-03-30 煤炭科学技术研究院有限公司 Low-permeability coal seam multi-scale micro-pore structure fine characterization method
CN112946005A (en) * 2021-02-02 2021-06-11 中国石油大学(华东) Shale microcrack evaluation method and application thereof
CN112946005B (en) * 2021-02-02 2022-06-07 中国石油大学(华东) Shale microcrack evaluation method and application thereof
CN113075102A (en) * 2021-03-11 2021-07-06 中国地质大学(北京) Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time
CN113075102B (en) * 2021-03-11 2022-02-01 中国地质大学(北京) Method for establishing mathematical model of relation between spontaneous imbibition amount of porous medium and time
CN113008927A (en) * 2021-04-01 2021-06-22 中海石油(中国)有限公司 Conversion method of nuclear magnetic resonance T2 spectrum and pore distribution of coal seam
CN113533156A (en) * 2021-06-30 2021-10-22 西安石油大学 Identification method for microscopic pore structure characteristics and multi-type pore fluid of shale oil reservoir
CN113533156B (en) * 2021-06-30 2024-02-02 西安石油大学 Shale oil reservoir microscopic pore structure characteristic and multi-class pore fluid identification method
CN113850689A (en) * 2021-09-24 2021-12-28 西南石油大学 Pore structure joint characterization method based on least square method and data processing method
CN115078210A (en) * 2022-05-05 2022-09-20 中国科学院地质与地球物理研究所 Shale pore structure testing method
CN115406921A (en) * 2022-08-27 2022-11-29 西安石油大学 Method for determining optimal degassing temperature of low-maturity hydrocarbon source rock
CN115753866A (en) * 2022-11-17 2023-03-07 重庆科技学院 Method for quantitatively characterizing shale pore structure
CN115753866B (en) * 2022-11-17 2024-04-09 重庆科技学院 Method for quantitatively representing shale pore structure
CN117110349A (en) * 2023-06-28 2023-11-24 东北石油大学 Quantitative evaluation method for Dan Duo-scale reservoir space wettability of shale oil reservoir rock
CN117110349B (en) * 2023-06-28 2024-05-28 东北石油大学 Quantitative evaluation method for wettability of Dan Duo-scale reservoir space of shale oil reservoir rock

Similar Documents

Publication Publication Date Title
CN108169099A (en) A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance
CN107727679B (en) A kind of characterization Deep Carbonate Rocks petrophysics characterization method
Wenlian et al. Evaluation of the stress sensitivity in tight reservoirs
US10996155B2 (en) Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability
CN103926186B (en) A kind of method that quantitative evaluation water filling affects throat distribution
Nabawy et al. Pore-throat characterization in highly porous and permeable sandstones
CN106153662A (en) The measuring method of rock core stress sensitivity
CN109030311A (en) Based on nuclear magnetic resonance T2Compose the pore structure classification and recognition methods of sensitive parameter
Cheng et al. Nuclear magnetic resonance features of low-permeability reservoirs with complex wettability
CN103926267A (en) Method for quantitatively evaluating pore throat change degree in stress sensitive process
CN102253069A (en) Method and device for determining permeability according to nuclear magnetic resonance T2 spectrum
CN104819923A (en) Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance
CN111007230B (en) Method for quantitatively evaluating oil content of low-porosity compact oil reservoir of continental-phase lake basin
US11035769B2 (en) Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability
Wang et al. Fractal analysis of tight sandstone petrophysical properties in unconventional oil reservoirs with NMR and rate-controlled porosimetry
US20240027379A1 (en) Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids
Wang et al. A novel model of predicting Archie’s cementation factor from nuclear magnetic resonance (NMR) logs in low permeability reservoirs
WO2020097037A1 (en) Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability
MAO et al. An improved method of using NMR T2 distribution to evaluate pore size distribution
Liu et al. A new method for determining tight sandstone permeability based on the characteristic parameters of the NMR T 2 distribution
Xu et al. Reducing residual oil saturation: underlying mechanism of imbibition in oil recovery enhancement of tight reservoir
CN109594971A (en) Flow net model method based on nuclear magnetic resonance log enhanced diffustion gas-bearing formation recognition factor
RU2542998C1 (en) Gas permeability anisotropy laboratory measurement on whole cores
Yanjie et al. Lower limit of tight oil flowing porosity: Application of high-pressure mercury intrusion in the fourth Member of Cretaceous Quantou Formation in southern Songliao Basin, NE China
CN110031496A (en) A method of evaluation compact reservoir movable fluid distribution characteristics

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180615