CN105974092A - Method for full-dimension representation and analysis of dense reservoir pore throats - Google Patents

Method for full-dimension representation and analysis of dense reservoir pore throats Download PDF

Info

Publication number
CN105974092A
CN105974092A CN201610536741.0A CN201610536741A CN105974092A CN 105974092 A CN105974092 A CN 105974092A CN 201610536741 A CN201610536741 A CN 201610536741A CN 105974092 A CN105974092 A CN 105974092A
Authority
CN
China
Prior art keywords
throat
distribution frequency
rock sample
pore
data
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
CN201610536741.0A
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.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
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 Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201610536741.0A priority Critical patent/CN105974092A/en
Publication of CN105974092A publication Critical patent/CN105974092A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Solid-Sorbent Or Filter-Aiding Compositions (AREA)

Abstract

The invention discloses a method for full-dimension representation and analysis of dense reservoir pore throats. The method comprises the following steps of 1, data optimization: (a) preparing a dense oil rock sample; (b) using a high-pressure mercury injection apparatus Pore Master to perform mercury injection experiment on a rock sample; (c) using an automatic isothermal adsorption instrument to perform low-temperature nitrogen experiment at the temperature of 77.35K and relative pressure of 0.01-1MPa; (d) according to adsorption type pore throat distribution frequency, accumulating a distribution frequency chart, and performing interpolating on the accumulated distribution frequency chart; (e) connecting the mercury injection type pore throat distribution frequency and the converted adsorption type pore throat distribution frequency at 25nm; 2, data check. The representation method has the advantages that the steps are simple; the pore throat distribution rules of different reservoirs can be analyzed, and the more accurate and deeper evaluation is performed on different dimensions of space microstructures by crossing multiple orders of magnitude.

Description

A kind of compact reservoir pore throat full size characterizes and the method for analysis
Technical field
The invention belongs to fine and close oil exploration field, be specifically related to a kind of compact reservoir pore throat full size and characterize and analysis method.
Background technology
Affected with exploitation by North America densification oil successful exploration, the whole world fine and close oil fast development, China, Russia, etc. more than 10 states actively launch correlational study.The fine and close oil of China is widely distributed, and technically recoverable resources reaches 20 × 108~25×108T, its effective exploitation is most important to Energy Development in China.Grand celebration, long celebrating, Sichuan and huge port are the big fine and close oil producing regions of China the most typical four.One of marked feature of fine and close oil is the finest and close for reservoir, extensive development micron-nanometer level pore throat, and based on nanoscale pore throat.Therefore, fine sign compact reservoir micro/nano-scale pore throat system is the primary premise of this petroleum industry of effective exploitation " dark fund ".Owing to compact reservoir pore throat crosses over multiple orders of magnitude, radius as little as about 1nm, current research, in addition to by tradition pore throat measuring technology, also introduces such as high-pressure mercury, gas absorption, FIB-SEM, nanometer CT equally accurate and the higher technology of resolution.
High-pressure mercury technology belongs to the popularization of conventional pressure hydrargyrum technology, will enter mercury pressure and rise to 60000psi(413.8MPa on the basis of routine pressure hydrargyrum technology), pore throat radius test range expansion is to 1.8nm-500 μm.Tian Hua etc. (2012) utilize high-pressure mercury technology bound gas adsorption technology to attempt obtaining the full size distribution of fine and close rock sample aperture, but its full size collection of illustrative plates belongs to the simple concatenation of different experiments method, and the jumping characteristic of curve is big, the effectiveness after data fusion is also without strict checking.High-pressure mercury is combined by Huang Zhenkai etc. (2013) with gas absorption, fine and close blue sheep interior microscopic hole has been carried out qualitative and quantitative study, but the relative amount of different size hole is had in his research, do not provide the concrete distribution characteristics of single voids content.Although high-pressure mercury method can obtain the substantially distribution situation of pore communication and different scale hole, but it is high to enter mercury pressure, reaching more than 400MPa, relevant scholar thinks that it easily causes artificial crack, and the test for micro-pore still exists bigger error.Therefore, this technology the most also needs to further investigate further.
Gas absorption, common technology is low temperature nitrogen absorption, and test pore radius is between 1nm-50nm.Hole is divided into 3 classes by pore radius size by IUPAC (IUPAC): micropore (<1nm), mesoporous (1 ~ 25nm), grand hole (>25nm).Under low-temp low-pressure, the isothermal adsorption of N2 is possible not only to effectively reflect rock sample micropore and mesoporous pore-size distribution, and can calculate specific surface area by BET equation, can be distributed with calculated pore by BJH equation, by isothermal adsorption and the measurable pore shape of desorption curve.Low temperature nitrogen absorption method is used for compact reservoir pore configuration research by the most numerous scholars, but the pore-size distribution jumping characteristic that nitrogen adsorption obtains is big, real effectively in the range of mesoporous interval, need to be combined this research method with additive method could embody reservoir micropore structure feature comprehensively.
FIB-SEM and nanometer CT, in current research, are mainly used in two dimension fine description and three-dimensionalreconstruction digital cores, are different from the quantitative assessment of high-pressure mercury and gas-adsorption technique.The regularity of distribution for overall evaluation micro throat, Philip H. Nelson(2009) use optical microscope, scanning electron microscope, pressure hydrargyrum, low-angle neutron scattering and calculate the methods such as chemistry, have studied oil, water, alkane and clay mineral molecular size, compact reservoir Pore throat size, institute's Research scale crosses over 7 orders of magnitude, obtains famous " Nelson plate ".Thereafter, Zou Caineng etc. (2012) obtain (2013) such as the whole world typical case " the Nelson plate " of unconventionaloil pool reservoir nanoscale throat distribution, Yang Hua and obtain " the Nelson plate " of Ordos Basin prolongation group densification oil reservoir throat distribution.But this achievement in research can only provide certain oil district reservoir rock sample pore throat substantially distribution, it is impossible to quantifies Pore throat size and the regularity of distribution of content thereof.
For quantitatively characterizing compact reservoir pore throat overall distribution rule, full size distribution in aperture is defined by the present invention, and deeply contrast high-pressure mercury and low temperature nitrogen adsorption experiment data, syncaryon magnetic resonance and centrefuge experiment, preferably valid data, set up compact reservoir pore throat full size distribution characterizing method.
Summary of the invention
The present invention proposes a kind of compact reservoir pore throat full size and characterizes and analysis method, and the method can analyze different reservoir throat distribution rule, preferably instructs fine and close oil effective exploitation, reduces development cost, improve development efficiency.
The technical scheme is that and be achieved in that: a kind of compact reservoir pore throat full size characterizes and the method for analysis, comprise the steps: one, data preferred:
A () prepares fine and close oil rock sample, test rock sample quality by balance, by slide gauge test rock sample length and diameter, and calculates rock sample volume, calculates rock sample density by quality and volume, oozes instrument test rock sample porosity by hole;
B () uses high-pressure mercury instrument Pore Master to carry out rock sample pressing hydrargyrum experiment, and by the percentage composition of total void volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. mercury injection method pore throat distribution frequency;
C () uses automatic isothermal adsorption instrument, it is 77.35K in temperature, low temperature nitrogen adsorption experiment is carried out, by the percentage composition of total pore size volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. absorption method throat distribution frequency between relative pressure 0.01-1MPa;
D () obtains absorption method pore throat cumulative distribution frequency diagram according to absorption method throat distribution frequency, and at the enterprising row interpolation of cumulative distribution frequency diagram, interpolation method selects point-to-point transmission linear interpolation or polynomial interopolation, interpolation point selects the point of mercury injection method pore throat cumulative distribution below frequency diagram 100nm, cumulative distribution frequency after interpolation is transformed into throat distribution frequency, and making absorption method throat distribution frequency translation after converting is revised throat distribution frequency;
The revised throat distribution frequency data of (e) mesoporous interval employing, and between grand porose area, use mercury injection method pore throat distribution frequency data, two kinds of data carry out interpolation at the separation 25nm in mesoporous and grand hole and are connected, and are converted into cumulative distribution frequency diagram, i.e. comprehensive throat distribution frequency;
Two, data detection:
F () makes parallel rock sample, and to rock sample saturation water, carry out 0.14MPa, 0.28MPa, 1.38MPa and 2.76MPa centrifugal force the most successively, and front and back carry out NMR test the most centrifugal, by asking for the ratio of the semaphore after being centrifuged and the semaphore of saturated-water phase, the voids content of below below below below 1 μm, 0.5 μm, 0.1 μm, 0.05 μm can be obtained, such that it is able to draw NMR-to be centrifuged voids content cumulative distribution table;
G the comprehensive throat distribution frequency translation preferably gone out is comprehensive pore throat cumulative distribution frequency by (), and be centrifuged data with parallel rock sample NMR-and be plotted in same cumulative distribution table, sees whether that surmounting NMR-is centrifuged the point of cumulative distribution;Such as, Fig. 5 is centrifuged data with NMR-at 100nm point first and overlaps, illustrate that the data preferably gone out are higher in the cumulative amount of below 100nm, the most total cumulative amount is 102%, more than 100%, illustrate that below 100nm hole cumulative amount is inevitable higher, need to be modified.
H data are modified by (), correction step is: the following voids content of (1) pair radius 2nm retrains, owing to below 2nm hole great majority are invalid hole, the content of less hole can be arranged to 0 when retraining, cumulative distribution is controlled in the range of fair-sized;(2) if still deviation theory relation when following for radius 2nm content is tied to ultimate value 0, then check that abnormity point place is interval, row constraint of going forward side by side;As in Fig. 5,100nm point still surmounts NMR-centered plots, then need 50nm-100nm section is tested, find the abnormity point in this interval, suitably retrain;If deviation theory relation at 500nm, then check that 100nm-500nm is interval, by that analogy, checks one by one and revises.
(i) revised comprehensive pore throat cumulative distribution frequency is again converted to comprehensive throat distribution frequency and i.e. obtains pore throat full size distribution characteristics collection of illustrative plates;
(j) analysis to full size distribution characteristics collection of illustrative plates: the contribution plot of rock sample permeability is plotted in revised comprehensive throat distribution frequency, such as the increase with permeability, on the right side of different lithology densification rock sample pore-size distribution, crest is raised or is moved to right, illustrate that the contribution of permeability is affected substantially by macrovoid, as permeability contribution rate curve main peak is gradually deviated from throat distribution curve main peak, penetrating power is made to rely primarily on the contribution of a small amount of macrovoid, part is from nanoscale hole, this situation can make the reservoir exploitation difficulty will be very big, and substrate fluid is extremely difficult to be employed.
Beneficial effects of the present invention is as follows: this characterizing method step is simple, different reservoir throat distribution rule can be analyzed by the method, cross over multiple order of magnitude and the space microstructure of reservoir different scale is carried out the most deep evaluation, preferably instruct fine and close oil effective exploitation, reduce development cost, improve development efficiency.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, the accompanying drawing used required in embodiment or description of the prior art will be briefly described below, apparently, accompanying drawing in describing below is only some embodiments of the present invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 for being transformed into revised distribution frequency figure by absorption method cumulative frequency;
Fig. 2 is the coordinate diagram carrying out interpolation connection at 25nm;
Fig. 3 is centrifugal distribution and actual distribution theory relation schematic diagram;
Fig. 4 is the centrifugal accumulation of the NMR-before revising and comprehensive pore throat cumulative distribution frequency diagram;
Fig. 5 is the centrifugal accumulation of revised NMR-and comprehensive pore throat cumulative distribution frequency diagram;
Fig. 6 is the different permeability rank tight sand permeability contribution rate of long celebrating and pore throat full size distribution characteristics atlas analysis figure.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on the embodiment in the present invention, the every other embodiment that those of ordinary skill in the art are obtained under not making creative work premise, broadly fall into the scope of protection of the invention.
The technical scheme is that and be achieved in that:
A kind of compact reservoir pore throat full size characterizes and the method for analysis, comprise the steps: one, data preferred:
A () prepares fine and close oil rock sample, test rock sample quality by balance, by slide gauge test rock sample length and diameter, and calculates rock sample volume, calculates rock sample density by quality and volume, oozes instrument test rock sample porosity by hole;
B () uses high-pressure mercury instrument Pore Master to carry out rock sample pressing hydrargyrum experiment, and by the percentage composition of total void volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. mercury injection method pore throat distribution frequency;
C () uses automatic isothermal adsorption instrument, it is 77.35K in temperature, low temperature nitrogen adsorption experiment is carried out, by the percentage composition of total pore size volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. absorption method throat distribution frequency between relative pressure 0.01-1MPa;
Mercury injection method and absorption method evaluation are the interval volume sizes of different pore throat, can calculate the percentage composition (distribution frequency, shown in Fig. 1 left side) of total pore size volume shared by the volume of different pore throat intervals by parameters such as rock sample density, quality, porositys.The reason that absorption method data jumping characteristic is big be exactly its test point very little, cause single pore throat interval to widen, volume content change is acutely.And the measuring point of mercury injection method is intensive, one group of data is more than 80 measuring points, the most single pore throat interval is the narrowest, interval two ends Pore throat size difference is the least, according to differential thought, interval pore volume content can be equivalent to the pore volume content of a certain size pore throat, the pressure hydrargyrum throat distribution relative smooth that therefore Fig. 1 obtains, meet conventional understanding.
For the feature that low temperature nitrogen adsorpting data jumping characteristic is bigger, can be at its enterprising row interpolation of cumulative distribution frequency diagram, after the cumulative distribution frequency after interpolation is transformed into throat distribution frequency, jumping characteristic substantially weakens (as shown in Fig. 1 right side).Interpolation method is as shown in step d:
nullD () obtains absorption method pore throat cumulative distribution frequency diagram according to absorption method throat distribution frequency,And at the enterprising row interpolation of cumulative distribution frequency diagram,Interpolation method selects point-to-point transmission linear interpolation or polynomial interopolation,Because the 100nm data below point that absorption method selects,But the data point that absorption method obtains is fewer,Make data jumping characteristic big,To such an extent as to error is bigger,If error to be reduced,The data point that pressure hydrargyrum obtains can be inserted adsorption number strong point by interpolation algorithm,Absorption method data point is made to become close,Solve the error problem of absorption method,So interpolation point selects the pore throat point of mercury injection method pore throat cumulative distribution frequency diagram (being converted to by mercury injection method pore throat distribution frequency) below 100nm (because of the point contained in mercury injection method data,Illustrate to exist the pore throat of corresponding size inside actual rock sample,What simply its volume content and absorption test obtained has difference),Cumulative distribution frequency after interpolation is transformed into throat distribution frequency,Making absorption method throat distribution frequency translation after converting is revised throat distribution frequency;
Obtain due to absorption method and mercury injection method is all pore throat volume content, it is assumed that in test process, pore throat is not destroyed, then same the slightly larger person of volume content, it is believed that measuring accuracy is higher.Because the volume content that test obtains is often slightly smaller than true volume content, only constantly close to actual value, precision just improves constantly.By by revised throat distribution frequency and mercury injection method pore throat distribution frequency Data Comparison, equally find, revised throat distribution frequency distribution is apparently higher than mercury injection method pore throat distribution frequency data, but close to (hole is divided into 3 classes by pore radius size by IUPAC (IUPAC): micropore (<1nm), mesoporous (1 ~ 25nm), grand hole (>25nm) when of 25nm), data trend changes, revised throat distribution frequency drops to below mercury injection method pore throat distribution frequency data, so needing to use e to walk:
E () is at the revised throat distribution frequency data of mesoporous interval employing, and between grand porose area, use mercury injection method pore throat distribution frequency data, two kinds of data carry out interpolation at the separation 25nm in mesoporous and grand hole and are connected that (Fig. 2 just can illustrate the interpolation method at 25nm, the data between two peak values before and after sample point selection closes on 25nm during interpolation, then 25nm size pore throat is carried out linear or polynomial interopolation, after interpolation obtains the content of 25nm size pore throat, it is attached mesoporous with grand hole), i.e. can get comprehensive throat distribution frequency, as shown in Figure 2;
The data reasonability preferably gone out need to be checked further, main reason is that two aspects: (1) high-pressure mercury ruptures micro-pore and causes large hole gap volume content higher, shows as cumulative distribution often beyond 100%;(2) low temperature nitrogen absorption test is mainly for the above hole of radius 1nm, the hole that the adsorbance of below 1nm can be attached between 1-2nm when processing by instrumentation program, showing as mesoporous distribution jumping characteristic big, after being combined with mercury pressuring data, cumulative distribution is likely to occur the situation more than 100%.
Whether the focusing on of inspection examines single voids content and cumulative distribution beyond normal value, and is modified.Check that whether single hole is beyond normal value, it is impossible to realize;Only check whether cumulative distribution exceedes normal value, it is easy to accomplish, but problematic data point cannot be found.Set forth herein that employing nuclear magnetic resonance, NMR combines centrifugal test (i.e. NMR test) and tests the data preferably gone out.
The step of data detection is:
F () makes parallel rock sample, and rock sample saturation water (is made parallel rock sample and refers to that making crucial physical property (includes lithology, porosity, permeability, the degree of depth, oil district etc.) rock sample identical with step a), carry out 0.14MPa the most successively, 0.28MPa, 1.38MPa and 2.76MPa centrifugal force, and front and back carry out NMR test the most centrifugal, by asking for the ratio of the semaphore after being centrifuged and the semaphore of saturated-water phase, below 1 μm can be obtained, below 0.5 μm, below 0.1 μm, voids content below 0.05 μm, such that it is able to draw NMR-to be centrifuged voids content cumulative distribution table;Owing to the existence of irreducible water makes corresponding semaphore bigger than normal, cause the cumulative amount of same point than cumulative actual content higher (Fig. 3, cumulative actual content is theoretic content, it is an imaginary curve), but it is as the aperture big this difference of change to be gradually reduced, final two curves tend to overlapping at cumulative amount 100%, just because of NMR-is centrifuged hole content distribution and there is above-mentioned theory relation with actual distribution, just achieve and the data preferably gone out are carried out on cumulative distribution table segmentation inspection.
G the comprehensive throat distribution frequency translation preferably gone out is comprehensive pore throat cumulative distribution frequency by (), and be centrifuged data with parallel rock sample NMR-and be plotted in same cumulative distribution table, sees whether that surmounting NMR-is centrifuged the point of cumulative distribution;Such as, Fig. 4 is centrifuged data with NMR-at 100nm point first and overlaps, illustrate that the data preferably gone out are higher in the cumulative amount of below 100nm, the most total cumulative amount is 102%, more than 100%, illustrate that below 100nm hole cumulative amount is inevitable higher, need to be modified.
H data are modified by (), correction step is: the following voids content of (1) pair radius 2nm retrains, owing to below 2nm hole great majority are invalid hole, the content of less hole can be arranged to 0 when retraining, cumulative distribution is controlled in the range of fair-sized (i.e. the theoretical value 100% of cumulative distribution, if it exceeds partial invalidity data point will suitably be retrained by 100%);(2) if when following for radius 2nm content is tied to ultimate value 0 still deviation theory relation (theory relation refers to that 1. cumulative distribution frequency should be 100%, or infinite approach 100%;2. comprehensive pore throat cumulative distribution frequency should be in below NMR-centered plots, it is impossible to has data point to cross NMR-centered plots), then check that abnormity point place is interval, row constraint of going forward side by side;As in Fig. 5,100nm point still surmounts NMR-centered plots, then need 50nm-100nm section is tested, find the abnormity point in this interval, suitably retrain;If deviation theory relation at 500nm, then checking that 100nm-500nm is interval, by that analogy, checks one by one and revises, revised structure is as shown in Figure 5.
(i) revised comprehensive pore throat cumulative distribution frequency is again converted to comprehensive throat distribution frequency and i.e. obtains pore throat full size distribution characteristics collection of illustrative plates;
(j) analysis to full size distribution characteristics collection of illustrative plates: the contribution plot of rock sample permeability is plotted in revised comprehensive throat distribution frequency, such as the increase with permeability, on the right side of different lithology densification rock sample pore-size distribution, crest is raised or is moved to right, illustrate that the contribution of permeability is affected substantially by macrovoid, as permeability contribution rate curve main peak is gradually deviated from throat distribution curve main peak, penetrating power is made to rely primarily on the contribution of a small amount of macrovoid, part is from nanoscale hole, this situation can make the reservoir exploitation difficulty will be very big, and substrate fluid is extremely difficult to be employed.
Application example:
The different permeability rank tight sand single pore throat permeability contribution rate Variations of long celebrating shows (shown in Fig. 6), and permeability is contributed essentially from submicron order hole, on a small quantity from nanoscale.When permeability reduction is to about 0.03mD, permeability contribution rate curve main peak is gradually deviated from throat distribution curve main peak so that penetrating power relies primarily on the contribution of a small amount of macrovoid, and part is from nanoscale hole.Therefore, long celebrating below 0.03mD reservoir exploitation difficulty will be very big, and substrate fluid is extremely difficult to be employed.
By aperture full size scattergram modal data, can also be averaged the evaluation of the physical parameters such as radius, main flow venturi, sorting coefficient, affected owing to these parameters are existed by microcrack calculating when, low volume data statistical result often regularity is the strongest, jumping characteristic is big, but by the available regular understanding of the statistics of mass data, next step should further investigate the regularity of distribution of different oil districts this type of physical parameter of compact reservoir.Studied by this, can be that resource assessment preferably provides crucial physical parameter, lay the foundation for seepage flow mechanism research simultaneously.
The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.

Claims (1)

1. compact reservoir pore throat full size characterizes and the method for analysis, comprise the steps: one, data preferred:
A () prepares fine and close oil rock sample, test rock sample quality by balance, by slide gauge test rock sample length and diameter, and calculates rock sample volume, calculates rock sample density by quality and volume, oozes instrument test rock sample porosity by hole;
B () uses high-pressure mercury instrument Pore Master to carry out rock sample pressing hydrargyrum experiment, and by the percentage composition of total pore size volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. mercury injection method pore throat distribution frequency;
C () uses automatic isothermal adsorption instrument, it is 77.35K in temperature, low temperature nitrogen adsorption experiment is carried out, by the percentage composition of total pore size volume shared by rock sample density, quality, porosity calculation difference pore throat interval volume, i.e. absorption method throat distribution frequency between relative pressure 0.01-1MPa;
D () obtains absorption method pore throat cumulative distribution frequency diagram according to absorption method throat distribution frequency, and at the enterprising row interpolation of cumulative distribution frequency diagram, interpolation method selects point-to-point transmission linear interpolation or polynomial interopolation, interpolation point selects the point of mercury injection method pore throat cumulative distribution below frequency diagram 100nm, cumulative distribution frequency after interpolation is transformed into throat distribution frequency, and making absorption method throat distribution frequency translation after converting is revised throat distribution frequency;
The revised throat distribution frequency data of (e) mesoporous interval employing, and between grand porose area, use mercury injection method pore throat distribution frequency data, two kinds of data carry out interpolation at the separation 25nm in mesoporous and grand hole and are connected, and i.e. obtain comprehensive throat distribution frequency;
Two, data detection:
F () makes parallel rock sample, and to rock sample saturation water, carry out 0.14MPa, 0.28MPa, 1.38MPa and 2.76MPa centrifugal force the most successively, and front and back carry out NMR test the most centrifugal, by asking for the ratio of the semaphore after being centrifuged and the semaphore of saturated-water phase, the voids content of below below below below 1 μm, 0.5 μm, 0.1 μm, 0.05 μm can be obtained, such that it is able to draw NMR-to be centrifuged voids content cumulative distribution table;
G the comprehensive throat distribution frequency translation preferably gone out is comprehensive pore throat cumulative distribution frequency by (), and be centrifuged data with parallel rock sample NMR-and be plotted in same cumulative distribution table, sees whether that surmounting NMR-is centrifuged the point of cumulative distribution;Such as, Fig. 5 is centrifuged data with NMR-at 100nm point first and overlaps, illustrate that the data preferably gone out are higher in the cumulative amount of below 100nm, the most total cumulative amount is 102%, more than 100%, illustrate that below 100nm hole cumulative amount is inevitable higher, need to be modified;
H data are modified by (), correction step is: the following voids content of (1) pair radius 2nm retrains, owing to below 2nm hole great majority are invalid hole, the content of less hole can be arranged to 0 when retraining, cumulative distribution is controlled in the range of fair-sized;(2) if still deviation theory relation when following for radius 2nm content is tied to ultimate value 0, then check that abnormity point place is interval, row constraint of going forward side by side;As in Fig. 5,100nm point still surmounts NMR-centered plots, then need 50nm-100nm section is tested, find the abnormity point in this interval, suitably retrain;If deviation theory relation at 500nm, then check that 100nm-500nm is interval, by that analogy, checks one by one and revises;
(i) revised comprehensive pore throat cumulative distribution frequency is again converted to comprehensive throat distribution frequency and i.e. obtains pore throat full size distribution characteristics collection of illustrative plates;
(j) analysis to full size distribution characteristics collection of illustrative plates: the contribution plot of rock sample permeability is plotted in revised comprehensive throat distribution frequency, such as the increase with permeability, on the right side of different lithology densification rock sample pore-size distribution, crest is raised or is moved to right, illustrate that the contribution of permeability is affected substantially by macrovoid, as permeability contribution rate curve main peak is gradually deviated from throat distribution curve main peak, penetrating power is made to rely primarily on the contribution of a small amount of macrovoid, part is from nanoscale hole, this situation can make the reservoir exploitation difficulty will be very big, and substrate fluid is extremely difficult to be employed.
CN201610536741.0A 2016-07-08 2016-07-08 Method for full-dimension representation and analysis of dense reservoir pore throats Pending CN105974092A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610536741.0A CN105974092A (en) 2016-07-08 2016-07-08 Method for full-dimension representation and analysis of dense reservoir pore throats

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610536741.0A CN105974092A (en) 2016-07-08 2016-07-08 Method for full-dimension representation and analysis of dense reservoir pore throats

Publications (1)

Publication Number Publication Date
CN105974092A true CN105974092A (en) 2016-09-28

Family

ID=56951343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610536741.0A Pending CN105974092A (en) 2016-07-08 2016-07-08 Method for full-dimension representation and analysis of dense reservoir pore throats

Country Status (1)

Country Link
CN (1) CN105974092A (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106442268A (en) * 2016-10-31 2017-02-22 中国科学技术大学 Method for detecting pore size distribution of shale mesopores
CN106526079A (en) * 2016-10-27 2017-03-22 中国石油大学(北京) Method used for studying tight sandstone pore-throat structure dynamic change
CN106979917A (en) * 2017-05-17 2017-07-25 成都理工大学 A kind of characterizing method of shale gas RESERVOIR PORE STRUCTURE and the evaluation method of shale gas reservoir
CN108802073A (en) * 2017-05-05 2018-11-13 中国石油化工股份有限公司 A kind of rock electrical parameters acquisition methods and device based on digital cores
CN109060623A (en) * 2018-06-30 2018-12-21 西南石油大学 Tight sand pore size distribution combines characterizing method
CN109283114A (en) * 2018-11-07 2019-01-29 西安石油大学 A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat
CN109424363A (en) * 2017-08-30 2019-03-05 中国石油天然气股份有限公司 A kind of Fluid Identification Method based on pore throat character and resistivity
CN109839401A (en) * 2019-01-29 2019-06-04 太原理工大学 A kind of judgement and processing method of goaf seam area
CN110542752A (en) * 2019-10-09 2019-12-06 中国科学院地理科学与资源研究所 Slope runoff simulation device and using method thereof
CN110726778A (en) * 2018-07-17 2020-01-24 中国石油化工股份有限公司 Method and system for judging lower limit of flowing croup of shale oil
CN111189758A (en) * 2020-01-13 2020-05-22 闽江学院 Method for full-scale characterization of pore throats of compact sandstone or shale oil reservoir
CN111880234A (en) * 2020-06-24 2020-11-03 西南石油大学 Fine identification method for tight reservoir fractures based on conventional well logging
CN112461728A (en) * 2020-11-12 2021-03-09 成都理工大学 Method and system for determining whole core permeability and pore throat radius distribution
CN112577881A (en) * 2020-12-31 2021-03-30 煤炭科学技术研究院有限公司 Low-permeability coal seam multi-scale micro-pore structure fine characterization method
CN112964616A (en) * 2021-02-05 2021-06-15 福州大学 Method for determining connectivity of different-scale pores of cement-based material
CN113607621A (en) * 2021-07-30 2021-11-05 中国矿业大学 Characterization method for closed pores and connectivity of coal-series composite reservoir
CN113850689A (en) * 2021-09-24 2021-12-28 西南石油大学 Pore structure joint characterization method based on least square method and data processing method
CN115267141A (en) * 2022-09-26 2022-11-01 中国石油天然气股份有限公司 Reservoir comprehensive evaluation method and device
US11982609B1 (en) * 2022-11-04 2024-05-14 Southwest Petroleum University Method for determining lower radius limit of movable throat of shale

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030106995A1 (en) * 2001-12-12 2003-06-12 Exxonmobil Upstream Research Company Method for measuring absorbed and interstitial fluids
WO2012123882A3 (en) * 2011-03-15 2012-12-27 Schlumberger Technology B.V. Examination of porosity by nmr and intrusion porosimetry
CN105067795A (en) * 2015-07-29 2015-11-18 中国石油大学(北京) Method for evaluating microstructure of compact reservoir rock
CN105466830A (en) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 A pore-throat dimension distribution identification method for reservoir sandstone
WO2016065127A1 (en) * 2014-10-23 2016-04-28 Chevron U.S.A. Inc. A system and method of pore type classification for petrophysical rock typing
WO2016065251A1 (en) * 2014-10-23 2016-04-28 Aramco Services Company Measuring tensile strength of tight rock using electromagnetic heating
CN105547959A (en) * 2015-12-30 2016-05-04 中国石油天然气股份有限公司 Analysis method of natural gas storage and collection capability of matrix pores in dense bedrock

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030106995A1 (en) * 2001-12-12 2003-06-12 Exxonmobil Upstream Research Company Method for measuring absorbed and interstitial fluids
WO2012123882A3 (en) * 2011-03-15 2012-12-27 Schlumberger Technology B.V. Examination of porosity by nmr and intrusion porosimetry
WO2016065127A1 (en) * 2014-10-23 2016-04-28 Chevron U.S.A. Inc. A system and method of pore type classification for petrophysical rock typing
WO2016065251A1 (en) * 2014-10-23 2016-04-28 Aramco Services Company Measuring tensile strength of tight rock using electromagnetic heating
CN105067795A (en) * 2015-07-29 2015-11-18 中国石油大学(北京) Method for evaluating microstructure of compact reservoir rock
CN105466830A (en) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 A pore-throat dimension distribution identification method for reservoir sandstone
CN105547959A (en) * 2015-12-30 2016-05-04 中国石油天然气股份有限公司 Analysis method of natural gas storage and collection capability of matrix pores in dense bedrock

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖前华: "典型致密油区储层评价及渗流机理研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106526079A (en) * 2016-10-27 2017-03-22 中国石油大学(北京) Method used for studying tight sandstone pore-throat structure dynamic change
CN106442268A (en) * 2016-10-31 2017-02-22 中国科学技术大学 Method for detecting pore size distribution of shale mesopores
CN106442268B (en) * 2016-10-31 2018-11-27 中国科学技术大学 A kind of detection method of shale mesoporous pore size distribution
CN108802073A (en) * 2017-05-05 2018-11-13 中国石油化工股份有限公司 A kind of rock electrical parameters acquisition methods and device based on digital cores
CN108802073B (en) * 2017-05-05 2021-03-12 中国石油化工股份有限公司 Rock electrical parameter obtaining method and device based on digital rock core
CN106979917A (en) * 2017-05-17 2017-07-25 成都理工大学 A kind of characterizing method of shale gas RESERVOIR PORE STRUCTURE and the evaluation method of shale gas reservoir
CN106979917B (en) * 2017-05-17 2019-08-27 成都理工大学 A kind of characterizing method of shale gas RESERVOIR PORE STRUCTURE and the evaluation method of shale gas reservoir
CN109424363A (en) * 2017-08-30 2019-03-05 中国石油天然气股份有限公司 A kind of Fluid Identification Method based on pore throat character and resistivity
CN109424363B (en) * 2017-08-30 2021-11-02 中国石油天然气股份有限公司 Fluid identification method based on pore throat structure and resistivity
CN109060623A (en) * 2018-06-30 2018-12-21 西南石油大学 Tight sand pore size distribution combines characterizing method
CN110726778A (en) * 2018-07-17 2020-01-24 中国石油化工股份有限公司 Method and system for judging lower limit of flowing croup of shale oil
CN109283114A (en) * 2018-11-07 2019-01-29 西安石油大学 A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat
CN109839401A (en) * 2019-01-29 2019-06-04 太原理工大学 A kind of judgement and processing method of goaf seam area
CN109839401B (en) * 2019-01-29 2022-02-15 太原理工大学 Method for judging and processing goaf fracture development area
CN110542752A (en) * 2019-10-09 2019-12-06 中国科学院地理科学与资源研究所 Slope runoff simulation device and using method thereof
CN111189758A (en) * 2020-01-13 2020-05-22 闽江学院 Method for full-scale characterization of pore throats of compact sandstone or shale oil reservoir
CN111880234A (en) * 2020-06-24 2020-11-03 西南石油大学 Fine identification method for tight reservoir fractures based on conventional well logging
CN111880234B (en) * 2020-06-24 2022-09-30 西南石油大学 Fine identification method for tight reservoir fractures based on conventional well logging
CN112461728A (en) * 2020-11-12 2021-03-09 成都理工大学 Method and system for determining whole core permeability and pore throat radius distribution
CN112461728B (en) * 2020-11-12 2021-07-02 成都理工大学 Method and system for determining whole core permeability and pore throat radius distribution
CN112577881A (en) * 2020-12-31 2021-03-30 煤炭科学技术研究院有限公司 Low-permeability coal seam multi-scale micro-pore structure fine characterization method
CN112964616A (en) * 2021-02-05 2021-06-15 福州大学 Method for determining connectivity of different-scale pores of cement-based material
CN112964616B (en) * 2021-02-05 2022-05-10 福州大学 Method for determining connectivity of different-scale pores of cement-based material
CN113607621A (en) * 2021-07-30 2021-11-05 中国矿业大学 Characterization method for closed pores and connectivity of coal-series composite reservoir
WO2023005600A1 (en) * 2021-07-30 2023-02-02 中国矿业大学 Method for characterizing closed pores and connectivity of coal measure composite reservoir
US11781962B2 (en) 2021-07-30 2023-10-10 China University Of Mining And Technology Characterization method of closed pores and connectivity of coal measure composite reservoirs
CN113850689A (en) * 2021-09-24 2021-12-28 西南石油大学 Pore structure joint characterization method based on least square method and data processing method
US20230096906A1 (en) * 2021-09-24 2023-03-30 Southwest Petroleum University Joint characterization method of pore structures based on least square method and data processing method
US11668637B2 (en) * 2021-09-24 2023-06-06 Southwest Petroleum University Joint characterization method of pore structures based on least square method and data processing method
CN115267141A (en) * 2022-09-26 2022-11-01 中国石油天然气股份有限公司 Reservoir comprehensive evaluation method and device
US11982609B1 (en) * 2022-11-04 2024-05-14 Southwest Petroleum University Method for determining lower radius limit of movable throat of shale

Similar Documents

Publication Publication Date Title
CN105974092A (en) Method for full-dimension representation and analysis of dense reservoir pore throats
Liu et al. Experimental study on the pore structure fractals and seepage characteristics of a coal sample around a borehole in coal seam water infusion
CN108169099A (en) A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance
CN110927035A (en) Method for calculating saturation of irreducible tight sandstone bound water
Xiao et al. A full-scale characterization method and application for pore-throat radius distribution in tight oil reservoirs
CN109632594B (en) Method for characterizing pore throat characteristics of tight reservoir on basis of high-pressure mercury intrusion multi-scale
CN105574320A (en) Evaluation method for effective seepage capability of low-permeability sandstone reservoir
Meng et al. Rock fabric of tight sandstone and its influence on irreducible water saturation in Eastern Ordos Basin
CN108959767A (en) A kind of narrow river channel type gas reservoir difference well type condensate injury method for numerical simulation
CN114609010B (en) Method and device for measuring oil-water relative permeability of shale reservoir
CN107014728A (en) A kind of hole measuring method
Qian et al. Quantitative estimates of non-Darcy groundwater flow properties and normalized hydraulic aperture through discrete open rough-walled joints
Zhang et al. Experimental study of the pore structure characterization in shale with different particle size
CN110441209A (en) A method of rock permeability is calculated based on compact reservoir digital cores
Ng et al. Effects of soil suction on dilatancy of an unsaturated soil
Shen et al. The influence of pore structure on water flow in rocks from the Beibu Gulf oil field in China
Lei et al. A new method for physical simulation of flow characteristics of fluids in shale oil reservoirs
Zhao et al. Applicability comparison of nuclear magnetic resonance and mercury injection capillary pressure in characterisation pore structure of tight oil reservoirs
CN110984948A (en) Method for testing and evaluating injection capability of fracturing propping agent in fracture
CN111664825B (en) Method for improving accuracy of pore structure parameter measurement
Yao et al. Quantitative investigation on natural gas flooding characteristics in tight oil cores after fracturing based on nuclear magnetic resonance technique
Sun et al. Numerical model about the influence of cyclic load on the hydraulic conductivity of soft clay
CN114577677A (en) Reservoir heterogeneity-based compact reservoir nonlinear seepage characterization method
Wang et al. Micro-mechanical behaviour of artificially cemented sands under compression and shear
Jones et al. Correlation of radial flow-through and hollow cylinder dynamic pressurization test for measuring permeability

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20160928

RJ01 Rejection of invention patent application after publication