CN109443867B - The method that the physical parameter of a kind of pair of tight rock is continuously detected - Google Patents
The method that the physical parameter of a kind of pair of tight rock is continuously detected Download PDFInfo
- Publication number
- CN109443867B CN109443867B CN201811258760.7A CN201811258760A CN109443867B CN 109443867 B CN109443867 B CN 109443867B CN 201811258760 A CN201811258760 A CN 201811258760A CN 109443867 B CN109443867 B CN 109443867B
- Authority
- CN
- China
- Prior art keywords
- core sample
- water
- rock
- curve
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 43
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 83
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 19
- 239000000344 soap Substances 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 13
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 12
- 239000000835 fiber Substances 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 230000035699 permeability Effects 0.000 claims description 50
- 238000005259 measurement Methods 0.000 claims description 19
- 230000001186 cumulative effect Effects 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 16
- 238000004519 manufacturing process Methods 0.000 claims description 15
- 239000011148 porous material Substances 0.000 claims description 12
- 238000002474 experimental method Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 7
- 238000005325 percolation Methods 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 4
- 238000001035 drying Methods 0.000 claims description 3
- 229920001971 elastomer Polymers 0.000 claims description 3
- 239000011521 glass Substances 0.000 claims description 3
- 239000004816 latex Substances 0.000 claims description 3
- 229920000126 latex Polymers 0.000 claims description 3
- 238000002386 leaching Methods 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000006837 decompression Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 239000008398 formation water Substances 0.000 claims description 2
- 238000005498 polishing Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000005514 two-phase flow Effects 0.000 claims description 2
- 239000011800 void material Substances 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 43
- 239000003350 kerosene Substances 0.000 description 10
- 238000012986 modification Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002427 irreversible effect Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
- G01N2001/2873—Cutting or cleaving
Landscapes
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Sampling And Sample Adjustment (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses the methods that the physical parameter to tight rock is continuously detected, it include: selection rock sample, rock sample is dried, then its size and weight is measured, take water flooding, and infiltrate the fiber containing hole, then rock core is rolled repeatedly on fiber, soak outer surface of the rock core in addition to both ends of the surface uniformly, then the good Nuclear Magnetic Resonance for putting it into debugging obtains the T2 spectral curve of sample, it further takes out core sample and dries, saturated process is being done to rock core using water flooding, and measure T2 spectral curve again with Nuclear Magnetic Resonance；Core sample is removed and placed in core holding unit again, apply displacement straightening and reaches suitable displacement pressure difference to it, after the system stabilizes, parameter needed for accurate recording, then core sample is continued to measure using soap film flowmeter, its T2 spectral curve is finally measured again, and collected data are subjected to analytical calculation, obtains required many kinds of parameters.The present invention has accurate reliable, the time saving convenient and fast beneficial effect of data.
Description
Technical field
The present invention relates to oil-gas field development technical field, the physical parameter of specifically a kind of pair of tight rock is continuously examined
The method of survey.
Background technique
As the exploitation of conventional gas and oil resource is close to saturation, the technical research of untraditional reservoir is increasingly deepened,
On the basis of this, fine and close oil gas has become the hot spot of petroleum exploration and development in recent years.It is the porosity of reservoir rock, permeability, opposite
The physical parameters such as permeability and saturation degree are the basic indexs for describing fine and close oil-gas reservoir reservoir characteristic.Rock basic physical properties parameter
It is accurate to determine, it is of great significance to oil-gas reservoir economy, Efficient Development.
The measurement of rock sample porosity needs to obtain total volume, the skeleton volume of rock sample and the pore volume of rock sample of rock sample.
Currently, pressing its measurement method, it is broadly divided into following a few classes:
(1) gas-liquids method: geometric measurement method is applied in total volume measurement.Regular rock sample very simple for shape, such as
Cylindrical body, and there is no notch or groove, the diameter and length of rock core directly can be directly measured with vernier caliper, then rock core
Volume can be calculated by cylindrical body.Skeleton volume is then measured with Boyle law.The total volume of rock sample subtracts skeleton volume, is exactly
The pore volume of rock sample, finally, obtaining the porosity of rock sample with the ratio calculation of pore volume and total volume.
(2) it is saturated kerosene method: weighing the aerial weight of rock sample with electronic balance.Saturation kerosene is vacuumized, so
Weight and kerosene specific gravity of the rock sample in kerosene, air after weighing saturation kerosene afterwards.To be saturated the rock sample after kerosene in air
Difference with the weight in kerosene calculates the total volume of rock sample divided by the specific gravity of kerosene；To be saturated kerosene and rock sample in air
Weight difference divided by kerosene specific gravity, calculate pore volume；The hole of rock sample is finally obtained divided by the total volume of rock sample with pore volume
Porosity.
(3) using Boyle law as foundation, the relationship of pressure and volume, the totality of measure and calculation rock sample gas method: are utilized
Two in product, skeleton volume and pore volume, the porosity of rock sample is calculated with porosity defined formula.
When measuring the permeability of sample, steady state method is generally used.And in steady state method measurement rock sample permeability process, usually
Be according to contrived experiment with Darcy's law, wherein rate of discharge is the important parameter for obtaining rock sample permeability, at present frequently with
Following several method:
(1) flowmeter method: one group of range in parallel electronic flowmeter not of uniform size in outlet end is flowed through by measuring medium
The size of sample flow selects suitable flowmeter.
(2) soap-film discharge method: the quantity tube of a certain specification is connected in outlet end, in quantity tube one end soap-film, passes through acquisition
The time that soap-film flows through certain scale calculates flow.
When measuring sample saturation degree, centrifugal process or displacement process are mostly used:
(1) centrifugal process: the rock sample of abundant saturation water flooding is put into a centrifuge, the weight difference of centrifugation front and back rock sample is passed through
Rock sample saturation degree is acquired in conjunction with weight method.
(2) displacement process: the rock sample of abundant saturation water flooding is put into driving device, rock sample is driven by gas drive water
It replaces, rock sample is taken out when the weight of displacement to rock sample is no longer changed, and acquire rock sample saturation degree in conjunction with weight method.
But detection method on the market still has following deficiency at present: (1) porosity comes according to Boyle law mostly
Measurement, what this method obtained is absolute porosity, and in production more concerned with effecive porosity, gas method is unable to measure rock
Effecive porosity, at present frequently with method be saturation weight method, but operation link is more, many times weigh rock core weight, by
Human factor is larger.(2) permeability mostly uses gas to survey and surveys two classes with liquid, is easy to produce kirschner effect in gas survey, causes to measure
As a result bigger than normal.When liquid surveys method measurement permeability, especially often it is difficult to inject liquid into rock core when measuring tight rock, usually
Higher production pressure is needed, this requires equipment voltage endurance capabilities with higher, and compact rock core is extremely difficult to surely
Determine state.(3) it then mostly uses centrifugal process or displacement process to obtain for saturation degree, but since tight rock interstitial space is small, passes
The metering device of system is difficult accurately to measure the variation of fluid miniscule volume in rock core.(4) these physical parameters were measured in the past all
Different experiments method is needed, a parameter is measured and just needs to apply rock core certain confining pressure, and reverse cyclic loadings confining pressure can make rock
Stone deforms, and it is irreversible that these, which deform some and are reversible some, if caused after measuring a parameter to rock core
Irreversible deformation, the parameter of subsequent measurement just do not meet actual conditions.
Summary of the invention
In view of the above-mentioned problems, the invention discloses the method that the physical parameter of a kind of pair of tight rock is continuously detected,
Having for rock is accurately calculated first with the variation characteristic and the variation of nuclear magnetic resonance area under the curve of rock sample weight before and after self-priming
Porosity is imitated, while remaining water saturation is accurately obtained in conjunction with nuclear magnetic resonance using cold store enclosure measurement Gas And Water Relative Permeability
Degree and corresponding gas phase relative permeability, finally by the ratio of the gas phase permeability and gas phase relative permeability that are measured in experiment
The absolute permeability of rock can be obtained in value.The disadvantage that the present invention overcomes previous experimental methods is single, metering is inaccurate etc., can be effective
Measure the basic physical properties parameter of tight rock.The present invention has accurate reliable, the time saving convenient and fast beneficial effect of data.
Technical scheme is as follows:
The method that the physical parameter of a kind of pair of tight rock is continuously detected, includes the following steps:
S1, selection need to carry out the rock sample of detection zone, carry out cutting polishing to it, then carry out washing oil and the desalinization of soil by flooding or leaching；
S2, rock sample is dried, then measures its size, and record its length L, diameter D, quality m；
S3, configuration sample enough water floodings, it is made to meet the requirement of the underground water density p of institute's detection zone；
S4, the fiber containing hole, and the stratum water-soaked fiber obtained with step S3 are taken, then by rock core in fiber
On repeatedly roll repeatedly, soak outer surface of the rock core in addition to both ends of the surface uniformly, as 1# core sample, record weight at this time
And it is recorded as m1；
S5, Nuclear Magnetic Resonance has been debugged, has set parameter, be then put into 1# core sample, and obtained the T2 spectrum of sample
Curve；
S6,1# core sample is taken out, drying measures its weight and is recorded as m3；
S7, saturated process is done to rock core using water flooding, as 2# core sample, 2# core sample is taken out after being saturated,
It weighs its weight and is recorded as m4；
S8, it repeats step S5 and is put into 2# core sample, measure T2 spectral curve；
S9, core sample 2 is taken out, and puts it into core holding unit, and applied displacement straightening and reach suitable to it
Displacement pressure difference, and timing, after the system stabilizes, when accurate recording water breakthrough time, cumulative water production, cumulative gas, displacement
Between and displacement pressure；
S10, when being driven to the weight that balance measures and being no longer changed, test the gas phase permeability under remaining water state,
2# core sample is continued to measure using soap film flowmeter, obtains the sampling flow q as unit of mL/ming；
After the completion of S11, test, 2# rock core sample is taken out, is weighed and is recorded as m5, then repeatedly step S5, measurement
Its T2 spectral curve；
S12, collected data are subjected to analytical calculation, obtain porosity, saturation degree, participates in saturation degree, is remaining watery
Gas phase permeability, rock absolute permeability under state, and draw phase percolation curve.
Further, in the step S4, step S4 is repeated, separated in time repeatedly weighs rock core, directly
Until the mass conservation of rock sample, the weight recorded at this time is m2, then carry out step S5.
It further, can scrambled record gas production, water yield at experiment initial stage in step S9；Gradually with water yield
It reduces, the record time can be extended.
Further, the particular content of step S10 is: when being driven to the weight that balance measures and being no longer changed (continuous 4
Hour is constant), the gas phase permeability under remaining water state is tested, the pipe for measuring bottle entrance is unloaded, is connect with hose
The air inlet of upper soap-film glass tube injects appropriate soap lye in the latex water dropper of soap film flowmeter, opens sampler, adjusts sampling
The flow of device gently pinches rubber head after waiting samplers to stablize sampling one minute, increase soap lye and generates a soap-film, writes down soap with stopwatch
Film is by certain volume required time, the sampling flow q being converted into as unit of mL/ming, repeatedly, it is averaged.
Compared with prior art, the invention has the beneficial effects that:
1, the present invention sufficiently combine each physical parameter definition, Porosimeter combine at last rock core weight change feature and
The area of T2 spectral curve is calculated using the definition of porosity, and the prior art is typically all to be counted just for T2 spectral curve
It calculates, calculated result accuracy is poor；
2, what is studied at present is tight rock, and tight rock pore radius fine pore space is much smaller than conventional rock, tradition
Metering device be unable to measure the variation of subtle Fluid Volume in tight rock, this uses nuclear magnetic resonance T 2 spectrum curve detection rock
The variation of middle fluid volume greatly improves the precision of experimental data；
3, method energy one-shot measurement multi-group data of the invention can greatly save detection time, improve detection efficiency,
And without repeatedly pressurizeing to rock core, be able to maintain the unification of rock core shape, it is ensured that experimental data it is accurate reliable.
Detailed description of the invention
Fig. 1 is the flow chart that the present invention tests.
Fig. 2 water saturation schematic diagram；
Fig. 3 is air water phase percolation curve；
Fig. 4 is irreducible water saturation schematic diagram；.
In figure:
1 valve a, 2 valve b, 3 valve c, 4 high-pressure nitrogen bottles, 5 first pressure controllers, 6 core holding units, 7 measurement bottles,
8 electronic balances, 9 water volume measurement pipes, 10 control valves, 11 thermometers, 12 data acquisition control platforms, 13 confining pressures pump, 14 second pressures
Force controller, 15 wet flow indicators, 16 valve d, 17 second pressure gauges.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The invention discloses a kind of methods of physical parameter for measuring tight rock, and this method is by traditional measurement and nuclear-magnetism phase
In conjunction with accurately calculating rock first with variation characteristic and nuclear magnetic resonance the area under the curve variation of rock sample weight before and after self-priming
Effecive porosity, while residual water in conjunction with nuclear magnetic resonance is accurately obtained using cold store enclosure measurement Gas And Water Relative Permeability
Saturation degree and corresponding gas phase relative permeability, finally by the gas phase permeability and gas phase relative permeability measured in experiment
Ratio the absolute permeability of rock can be obtained.The disadvantage that the present invention overcomes previous experimental methods is single, metering is inaccurate etc., energy
The effectively basic physical properties parameter of measurement tight rock.The present invention has accurate reliable, the time saving convenient and fast beneficial effect of data.
In petroleum works field, nuclear magnetic resonance is widely used in core analysis in laboratory with quick, efficient, lossless advantage
With live well logging study, cardinal principle is: when applying the radiofrequency signal of certain frequency to the fluid in external magnetic field,
The phenomenon that proton in fluid absorbs the energy of radio frequency, generates resonance, the nuclear proton of hydrogen atom is by non-equilibrium state to equilibrium-like
The relaxation process of state is transverse relaxation and the combination of longitudinal relaxation, and corresponding is lateral relaxation time T2 and longitudinal relaxation time
T1.But when measuring petrophysical parameter with Nuclear Magnetic Resonance, operating procedure is more complex, and professional requirement is high, especially
When measuring the physical parameter of tight rock, it is difficult directly to measure exact value.
The method that the physical parameter of a kind of pair of tight rock is continuously detected, includes the following steps:
1. experimental procedure
(1) washing oil and the desalinization of soil by flooding or leaching.
(2) length, diameter and the dry weight of measurement drying rock sample, is denoted as L, D and m respectively.
(3) water flooding, precise measurement water flooding density p are filteredw。
(4) one fiber containing hole is soaked with water flooding, rock core is rolled on fiber, guarantee that rock core removes both ends of the surface
Except outer surface uniformly soak.
(5) weight of record at this time is m1, Nuclear Magnetic Resonance is opened, according to setting constants of nuclear magnetic resonance is illustrated, sample is put into
In Nuclear Magnetic Resonance, the T2 spectral curve of sample is obtained.
(6) step (4) being repeated, separated in time repeatedly weighs rock core, until the mass conservation of rock sample,
The weight of record at this time is m2。
(7) rock core is put into Nuclear Magnetic Resonance, after constants of nuclear magnetic resonance is set, obtains the T2 spectral curve of rock.
(8) it is dried 48 hours after taking out rock sample, re-measures its weight and be denoted as m3。
(9) sample is done into saturated process (vacuumizing saturation simulation water flooding, saturation time is 48h or more), after being saturated
Sample is taken out, example weight is weighed and is denoted as m4.Nuclear Magnetic Resonance is opened, according to setting constants of nuclear magnetic resonance is illustrated, by the sample of saturation
Product are put into Nuclear Magnetic Resonance, obtain the T2 spectral curve of sample.
(10) driving device is connected first, after whether check device there are the failures such as leak, gas leakage, is opened as shown in Figure 1
Valve a1, valve b2, valve c3 and valve d16, empty entire experimental provision, be then shut off valve a1, valve b2, valve c3 and
The rock core of saturation is fitted into core holding unit 6 by valve d16, opens valve b2, is pumped 13 pairs of rock samples by confining pressure and is applied centainly
Confining pressure and internal pressure.
(11) valve a1, valve c3 are opened, system pressure is applied to reasonable displacement pressure difference Pm, use time controller
11 setting some times were recorded, after the system stabilizes, accurate recording water breakthrough time, cumulative water production, cumulative gas production, drive
For time and displacement pressure, the initial stage of experiment can scrambled record gas production, water yield；With gradually decreasing for water yield, can prolong
The long record time.
(12) when being driven to the weight that balance measures and being no longer changed (continuous 4 hours constant) remaining water state, is tested
Under gas phase permeability: by measure bottle entrance pipe unload, the air inlet of soap-film glass tube is connected with hose, in soap
The appropriate soap lye of injection in the latex water dropper of film flowmeter, opens sampler, adjusts the flow of sampler, and samplers is waited to stablize sampling
Rubber head is gently pinched after one minute, increase soap lye and generates a soap-film, writes down soap-film by changing the time required to certain volume with stopwatch
The sampling flow q being counted as unit of mL/ming, repeatedly, it is averaged.
(13) after the completion of entire measurement process, internal pressure and confining pressure is unloaded, rock sample is weighed, weight at this time is recorded
For m5, sample is put into Nuclear Magnetic Resonance, the T2 spectral curve of sample is obtained.
2. data processing
2.1 porosity calculation
As shown in Fig. 2, the area that curve S1 is surrounded is A1, T2 spectral curve when rock sample sucks a part of water, curve are indicated
The area that S2 is surrounded is A2, and the T2 spectral curve of the rock sample after indicating abundant Absorb Water, the water saturation of rock sample is denoted as Δ at this time
Sw1。
The area that the curve of front and back two in T2 spectrum obtained by Nuclear Magnetic Resonance surrounds, can obtain saturation degree Δ Sw1Such as
Shown in formula (1):
It can get the active porosity volume of rock sample by the ratio of rock sample weight difference and stratum water density before and after self-priming:
Total pore volume of rock sample at this time are as follows:
According to the definition of porosity, the porosity of rock sample is total pore size volume VpWith the ratio of apparent volume V:
Wherein,
2.2 saturation computation
2.2.1 degree of saturation is verified
The active porosity volume of rock sample can be obtained by the ratio between the weight change of rock sample before and after sucking water and stratum water density
It is as follows:
Then the effecive porosity of rock sample is that the ratio between active porosity volume and apparent volume are shown below:
And if only ifWhen, rock sample saturation sufficiently, otherwise will repeat step (9), so that rock sample is sufficiently full
With until meeting the condition.
2.2.2 irreducible water saturation calculates
As shown in figure 4, the area that curve S4 is surrounded is A4, and the T2 spectral curve of the rock sample after indicating saturation water flooding, curve
The area that S3 is surrounded is A3, indicates the T2 spectral curve of the rock sample after gas drive water is stablized, the water saturation of rock sample as can at this time
Dynamic fluid saturation.
The area that the curve of front and back two in T2 spectrum obtained by Nuclear Magnetic Resonance surrounds can show that movable fluid is saturated
Degree is as shown in formula (6):
And because
In formula, Vq=VpΔSw2, as can the fluid that goes out of displacement, remaining water saturation can be obtained by being substituted into dissolve in formula (7)
Degree are as follows:
ΔSow=1- Δ Sw2 (8)
2.3 computing permeability
2.3.1 phase percolation curve is drawn
(1) value being adapted to the cumulative fluid yield values measured under rock sample ambient outlet pressures power under rock sample average pressure:
Vi=Δ Vwi+Vi-1+[2Pa/(Δp+2Pa)]ΔVgi (9)
In formula:
Vi--- cumulative fluid yield, mL；
Vi-1--- the above cumulative fluid yield of a bit, mL；
ΔVwi--- the water increment of some time interval, mL；
ΔVgi--- the gas increment of some time interval measured under ambient outlet pressures power, mL；
Δ P --- displacement pressure difference, MPa；
Pa--- atmospheric pressure when measurement, MPa.
(2) cumulative gas production ∑ V is drawng, cumulative water production ∑ VwWith the relation curve of accumulation injection length ∑ t.
(3) it is uniformly taken a little on curve, obtains corresponding gas production △ V in △ t at certain time intervalsgiWith △ Vwi, press
Following formula is calculated.
In formula:
Sg,av--- average gas saturation, %；
Vw--- accumulation outlet water, mL；
qgi--- gas flow when two-phase flow, mL/s；
qg--- gas flow when single-phase flow, mL/s；
fg--- void fraction (decimal)；
fw--- moisture content (decimal)；
μg--- injection gas viscosity, mPa s；
μw--- the viscosity of the simulated formation water of saturation rock sample, mPa s；
C --- decompression volume factor (decimal)；
P3--- rock sample inlet pressure (absolute), MPa；
P4--- rock sample outlet pressure (absolute), MPa.
(4) relation curve of Gas And Water Relative Permeability and water saturation is drawn, as shown in Figure 3.
2.3.2 the gas phase permeability under remaining water state calculates
It is as follows that gas phase permeability can be measured by soap film flowmeter:
In formula:
Kg--- gas permeability, μm2；
qg--- water flow, mL/s；
μg--- the viscosity of gas under measuring temperature, mPa s；
L --- rock sample length, cm；
A --- rock sample sectional area, cm2；
P1--- rock sample inlet pressure (gauge pressure), MPa；
P2--- rock sample outlet pressure (gauge pressure), MPa；
2.3.3 rock absolute permeability calculates
By the way that corresponding gas phase permeability value K under irreducible water saturation can be obtained on phase percolation curverg, and because residual water is full
It is equal to the ratio of gas phase permeability and rock absolute permeability with the gas phase permeability under degree, therefore the absolute of rock can be calculated
Shown in permeability such as formula (12):
The above described is only a preferred embodiment of the present invention, be not intended to limit the present invention in any form, though
So the present invention has been disclosed as a preferred embodiment, and however, it is not intended to limit the invention, any technology people for being familiar with this profession
Member, without departing from the scope of the present invention, when the technology contents using the disclosure above make a little change or modification
For the equivalent embodiment of equivalent variations, but anything that does not depart from the technical scheme of the invention content, according to the technical essence of the invention
Any simple modification, equivalent change and modification to the above embodiments, all of which are still within the scope of the technical scheme of the invention.
The above described is only a preferred embodiment of the present invention, be not intended to limit the present invention in any form, though
So the present invention has been disclosed as a preferred embodiment, and however, it is not intended to limit the invention, any technology people for being familiar with this profession
Member, without departing from the scope of the present invention, when the technology contents using the disclosure above make a little change or modification
For the equivalent embodiment of equivalent variations, but anything that does not depart from the technical scheme of the invention content, according to the technical essence of the invention
Any simple modification, equivalent change and modification to the above embodiments, all of which are still within the scope of the technical scheme of the invention.
Claims (3)
1. the method that the physical parameter of a kind of pair of tight rock is continuously detected, which comprises the steps of:
The core sample in the region that S1, selection are detected, carries out cutting polishing to it, then carries out washing oil and the desalinization of soil by flooding or leaching；
S2, core sample is dried, then measures its size, and record its length L, diameter D, quality m；
S3, configuration sample enough water floodings, it is made to meet the water flooding density p of institute's detection zonewIt is required that；
S4, the fiber containing hole, and the stratum water-soaked fiber obtained with step S3 are taken, then by core sample in fiber
On repeatedly roll repeatedly, soak outer surface of the core sample in addition to both ends of the surface uniformly, as 1# core sample, record at this time
Weight is simultaneously recorded as m1；
S5, Nuclear Magnetic Resonance has been debugged, has set parameter, be then put into 1# core sample, and obtained the T2 spectral curve of sample,
It is named as curve S1；Step S4 is repeated, separated in time repeatedly weighs core sample, until the quality of core sample
Until constant, the weight recorded at this time is m2；Core sample is put into Nuclear Magnetic Resonance, after constants of nuclear magnetic resonance is set, obtains rock
The T2 spectral curve of heart sample, is named as curve S2；
S6,1# core sample is taken out, drying measures its weight and is recorded as m3；
S7, saturated process is done to core sample using water flooding, as 2# core sample, 2# core sample is taken out after being saturated,
It weighs its weight and is recorded as m4；
S8, Nuclear Magnetic Resonance has been debugged, has set parameter, 2# core sample has been put into, measured T2 spectral curve, be named as curve
S3；
S9,2# core sample is taken out, and puts it into core holding unit, and applied displacement straightening and reach suitable displacement to it
Pressure difference, and timing, after the system stabilizes, the accurate recording water breakthrough time, cumulative water production, cumulative gas, the displacement time and
Displacement pressure；
S10, when being driven to the weight that balance measures and being no longer changed, test the gas phase permeability under remaining water state, use
Soap film flowmeter continues to measure to 2# core sample, obtains the gas sampling flow q as unit of mL/sg；
After the completion of S11, test, 2# core sample is taken out, is weighed and is recorded as m5, Nuclear Magnetic Resonance has then been debugged, has been set
Good parameter, 2# core sample is put into, and is measured its T2 spectral curve, is named as curve S4；
S12, collected data are subjected to analytical calculation, obtain the gas phase under porosity, irreducible water saturation, remaining water state
Permeability, rock absolute permeability, and draw phase percolation curve；Specifically include following data processing step:
S121, porosity calculation
The area that curve S1 is surrounded is A1, and the area that curve S2 is surrounded is A2, and the water saturation of core sample is denoted as Δ at this time
Sw1, Δ Sw1The formula as shown in formula (1) calculates:
It can get the active porosity body of core sample by the ratio of core sample weight difference and stratum water density before and after self-priming
Product:
In formula, ρwRepresent stratum water density；
Total pore volume of core sample at this time are as follows:
According to the definition of porosity, the porosity of core sample is total pore size volume VpWith the ratio of apparent volume V:
Wherein,
S122, saturation computation
(1) degree of saturation is verified
The active porosity of core sample can be obtained by the ratio between the weight change of core sample before and after sucking water and stratum water density
Volume is as follows:
Then the effecive porosity of core sample is that the ratio between active porosity volume and apparent volume are shown below:
And if only ifWhen, core sample saturation sufficiently, otherwise will repeat step S7, so that core sample is abundant
Saturation is until meeting the condition；
(2) irreducible water saturation calculates
The area that curve S4 is surrounded is A4, and the area that curve S3 is surrounded is A3, and the water saturation of core sample as can at this time
Dynamic fluid saturation, calculates shown in moveable gel such as formula (6):
And because
In formula, Vq=VpΔSw2, as can displacement go out fluid, irreducible water saturation can be obtained by being substituted into formula (7) are as follows:
ΔSow=1- Δ Sw2 (8)
S123, computing permeability
(1) phase percolation curve is drawn
The cumulative fluid yield values measured under core sample ambient outlet pressures power are adapted to core sample average pressure by the first step
Under value:
Vi=Δ Vwi+Vi-1+[2Pa/(ΔP+2Pa)]ΔVgi (9)
In formula:
Pa --- atmospheric pressure when measurement, MPa；
Vi--- cumulative fluid yield, mL；
Vi-1--- the above cumulative fluid yield of a bit, mL；
ΔVwi--- the water increment of some time interval, mL；
ΔVgi--- the gas increment of some time interval measured under ambient outlet pressures power, mL；
Δ P --- displacement pressure difference, MPa；
Second step draws cumulative gas production ∑ Vg, cumulative water production ∑ VwWith the relation curve of accumulation injection length ∑ t；
Third step uniformly takes a little on curve, obtains corresponding gas production △ V in △ t at certain time intervalsgiWith water yield △
Vwi, it is calculated as the following formula:
In formula: Krg--- gas phase relative permeability；
Krw--- water phase relative permeability；
K --- absolute permeability value；
A --- core sample sectional area, cm2；
△Vi--- some time interval adds up fluid-withdrawal rate, mL；
P3--- core sample inlet pressure, absolute pressure, MPa；
P4--- core sample outlet pressure, absolute pressure, MPa；
C --- decompression volume factor, decimal；
Sg,av--- average gas saturation, %；
Vw--- accumulation outlet water, mL；
qgi--- gas flow when two-phase flow, mL/s；
qga————Gas flow when single-phase flow, mL/s；
fg--- void fraction, decimal；
fw--- moisture content, decimal；
μg--- injection gas viscosity, mPas；
μw--- the viscosity of the simulated formation water of saturated core sample, mPas；
4th step, the relation curve for drawing Gas And Water Relative Permeability and water saturation；
(2) gas phase permeability under remaining water state calculates
It is as follows that gas phase permeability can be measured by soap film flowmeter:
In formula: Kg--- gas phase permeability, μm2；
qg--- gas sampling flow, mL/s；
μg--- the viscosity of gas, mPas under measuring temperature；
L --- core sample length, cm；
A --- core sample sectional area, cm2；
P1--- core sample inlet pressure, gauge pressure, MPa；
P2--- core sample outlet pressure, gauge pressure, MPa；
(3) rock absolute permeability calculates
By the way that corresponding gas phase relative permeability value K under irreducible water saturation can be obtained on phase percolation curverg, and because residual water is full
It is equal to the ratio of gas phase permeability and rock absolute permeability with the gas phase permeability under degree, therefore calculates the Absolute permeation of rock
Shown in rate K such as formula (12):
2. the method that the physical parameter of a kind of pair of tight rock according to claim 1 is continuously detected, feature exist
In in step S9, in experiment initial stage scrambled record gas production, water yield；With gradually decreasing for water yield, when extending record
Between.
3. the method that the physical parameter of a kind of pair of tight rock according to claim 1 is continuously detected, feature exist
In the particular content of step S10 is: when being driven to the weight that balance measures and being no longer changed, testing under remaining water state
The pipe for measuring bottle entrance is unloaded, the air inlet of soap-film glass tube is connected with hose, in soap-film stream by gas phase permeability
The appropriate soap lye of injection in the latex water dropper of meter, opens sampler, adjusts the flow of sampler, and samplers is waited to stablize one point of sampling
Rubber head is gently pinched after clock, increase soap lye and generates a soap-film, writes down soap-film by being converted into the time required to certain volume with stopwatch
Sampling flow q as unit of mL/sg, repeatedly, it is averaged.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811258760.7A CN109443867B (en) | 2018-10-26 | 2018-10-26 | The method that the physical parameter of a kind of pair of tight rock is continuously detected |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811258760.7A CN109443867B (en) | 2018-10-26 | 2018-10-26 | The method that the physical parameter of a kind of pair of tight rock is continuously detected |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109443867A CN109443867A (en) | 2019-03-08 |
| CN109443867B true CN109443867B (en) | 2019-08-09 |
Family
ID=65547681
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811258760.7A Active CN109443867B (en) | 2018-10-26 | 2018-10-26 | The method that the physical parameter of a kind of pair of tight rock is continuously detected |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109443867B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20230133700A1 (en) * | 2021-10-29 | 2023-05-04 | Halliburton Energy Services, Inc. | Nuclear magnetic resonance (nmr) fluid substitution using machine learning |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110715888B (en) * | 2019-05-30 | 2022-04-22 | 西南石油大学 | Method for measuring size distribution change of sulfur deposit pores of high-sulfur-content gas reservoir core |
| CN110296931B (en) * | 2019-08-13 | 2022-06-21 | 中国石油大学(华东) | Characterization method and system for oil-water relative permeability information of tight sandstone |
| CN112782477B (en) * | 2019-11-11 | 2024-05-14 | 中国石油化工股份有限公司 | Method and system for measuring electrical response characteristics of rock core in different wetting states |
| CN111337408B (en) * | 2020-03-27 | 2022-02-08 | 西南石油大学 | Method for testing rock crack porosity by using low-field nuclear magnetic resonance equipment |
| CN111638158A (en) * | 2020-07-17 | 2020-09-08 | 西南石油大学 | Compact sandstone gas-water phase permeability testing device and method based on capacitance method |
| CN111879678B (en) * | 2020-07-31 | 2022-06-07 | 西南石油大学 | Self-priming method-based experimental method for gas-water relative permeability of tight sandstone |
| CN112505084B (en) * | 2020-11-30 | 2022-10-25 | 中国石油大学(华东) | Evaluation model, evaluation method and application for improving shale oil mobility through gas injection |
| CN112505085B (en) * | 2021-02-05 | 2021-04-09 | 西南石油大学 | Method for measuring porosity effective stress coefficient based on nuclear magnetic resonance |
| CN113834840B (en) * | 2021-09-24 | 2024-05-14 | 西安工程大学 | Method for testing core imbibition efficiency |
| CN114858683B (en) * | 2022-05-20 | 2023-03-31 | 西南石油大学 | Method for evaluating pore structure change of gas storage core after salt precipitation |
| CN114894832A (en) * | 2022-05-27 | 2022-08-12 | 中国石油大学(华东) | Shale oil CO injection 2 Experimental device and method for evaluating electromagnetic heating auxiliary energy increasing and yield increasing effects |
| CN115078438A (en) * | 2022-06-19 | 2022-09-20 | 西南石油大学 | Method for establishing pore network model based on nuclear magnetic resonance test digital core |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105866009A (en) * | 2016-05-30 | 2016-08-17 | 中国石油大学(北京) | Method and device for calculating effective porosity of compact oil storage layer |
| CN106483057A (en) * | 2016-09-30 | 2017-03-08 | 西安石油大学 | A kind of method of quantitative assessment ultra-deep reservoir movable fluid and its application |
-
2018
- 2018-10-26 CN CN201811258760.7A patent/CN109443867B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105866009A (en) * | 2016-05-30 | 2016-08-17 | 中国石油大学(北京) | Method and device for calculating effective porosity of compact oil storage layer |
| CN106483057A (en) * | 2016-09-30 | 2017-03-08 | 西安石油大学 | A kind of method of quantitative assessment ultra-deep reservoir movable fluid and its application |
Non-Patent Citations (3)
| Title |
|---|
| 基于核磁共振T2图谱测量岩心裂缝孔隙度;宋南希 等;《勘探开发》;20170228(第2期);第127页 * |
| 基于核磁共振技术研究页岩自发渗吸过程;蒙冕模 等;《特种油气藏》;20151031;第22卷(第5期);第137-140页 * |
| 致密气藏气水相对渗透率理论及实验分析;莫邵元 等;《天然气地球科学》;20151130;第26卷(第11期);第2149-2153页 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20230133700A1 (en) * | 2021-10-29 | 2023-05-04 | Halliburton Energy Services, Inc. | Nuclear magnetic resonance (nmr) fluid substitution using machine learning |
| US11828901B2 (en) * | 2021-10-29 | 2023-11-28 | Halliburton Energy Services, Inc. | Nuclear magnetic resonance (NMR) fluid substitution using machine learning |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109443867A (en) | 2019-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109443867B (en) | The method that the physical parameter of a kind of pair of tight rock is continuously detected | |
| US10837278B2 (en) | Device and method for evaluating gas-water two-phase flow fracture conductivity of fractured shale | |
| CN110346258B (en) | Method for measuring relative permeability of oil phase of compact rock | |
| CN105910971B (en) | The simultaneous measuring method of rich organic matter compact rock core gas permeability and diffusion coefficient | |
| CN208171813U (en) | A kind of multi-functional permeability test device | |
| CN103954544B (en) | A kind of polymer control water increases experimental provision and the experimental technique of gas effect assessment | |
| CN109001243B (en) | Method and device for evaluating dynamic water lock effect of coal by adopting low-field nuclear magnetic resonance | |
| CN109030292B (en) | novel method for determining wettability of compact rock | |
| CN110793901B (en) | High-temperature high-pressure gas reservoir permeability flow rate sensitivity test method considering bound water | |
| CN206410978U (en) | A kind of tight rock gas phase relative permeability measurement apparatus | |
| CN210264648U (en) | Multifunctional rock core displacement device | |
| CN109283029B (en) | A kind of method, apparatus and clay preparing instrument measuring clay bound water and mechanics parameter | |
| CN104297126B (en) | low permeability reservoir gas seepage starting pressure gradient measuring device and measuring method | |
| CN112986097B (en) | Experimental measurement method for determining relative permeability curve of tight reservoir steady state method | |
| CN105910941B (en) | The test method of content of unfrozen water in frozen earth based on pressure plate apparatus | |
| CN105928832B (en) | The assay method of capillary viscosity analyzer and fluid viscosity | |
| CN107271347A (en) | The method that lignite effecive porosity is measured based on low-field nuclear magnetic resonance technology | |
| CN110296931A (en) | A kind of characterizing method and system of tight sand oil-water relative permeability information | |
| CN208155804U (en) | Compact reservoir rock electrical measurement based on nuclear magnetic resonance | |
| CN106501151A (en) | A kind of shale aperture measurement device and method based on imbibition and ion diffusion property | |
| CN108444890B (en) | Unsteady state titration device and method for testing medium and high permeability core liquid permeability | |
| CN113866069A (en) | Shale core permeability experiment device and method | |
| CN109799177A (en) | A kind of device and method multiple groups rock sample Non-Darcy Flow in Low Permeability Reservoir test while measured | |
| CN109580454A (en) | A method of compact reservoir Fluid Sensitivity is tested with pressure oscillation method | |
| CN209821099U (en) | Multifunctional compact gas reservoir dynamic parameter joint measurement device based on nuclear magnetic resonance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |