CN109655394A - A kind of nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters - Google Patents
A kind of nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters Download PDFInfo
- Publication number
- CN109655394A CN109655394A CN201811574964.1A CN201811574964A CN109655394A CN 109655394 A CN109655394 A CN 109655394A CN 201811574964 A CN201811574964 A CN 201811574964A CN 109655394 A CN109655394 A CN 109655394A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- permeability
- pore throat
- spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 120
- 239000011148 porous material Substances 0.000 title claims abstract description 103
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 96
- 238000001228 spectrum Methods 0.000 title claims abstract description 75
- 238000004364 calculation method Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000002790 cross-validation Methods 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000013507 mapping Methods 0.000 claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 17
- 238000013528 artificial neural network Methods 0.000 claims description 10
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 10
- 229910052753 mercury Inorganic materials 0.000 claims description 10
- 239000011435 rock Substances 0.000 claims description 9
- 210000000867 larynx Anatomy 0.000 claims description 8
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 210000002569 neuron Anatomy 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 239000004576 sand Substances 0.000 abstract description 2
- 230000005311 nuclear magnetism Effects 0.000 description 11
- 230000035515 penetration Effects 0.000 description 5
- 238000009738 saturating Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 235000016838 Pomo dAdamo Nutrition 0.000 description 2
- 244000003138 Pomo dAdamo Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Dispersion Chemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The present invention discloses the nuclear magnetic resonance T 2 spectrum calculation of permeability under a kind of pore throat character restriction on the parameters, comprising: S1: converts pseudo- capillary pressure curve for nuclear magnetic resonance T 2 spectrum curve；S2: pore throat character parameter is calculated according to the pseudo- capillary pressure curve；S3: pore throat character parameter, the NMR logging data being calculated after pseudo- capillary pressure curve are converted using the core analysis permeability data and nuclear magnetic resonance T 2 spectrum of acquisition, establishes pore throat character parameter, nuclear magnetic resonance T 2 spectrum-permeability corresponding relation data body；S4: selected non-linear map establishes Permeability Prediction model；S5: cross validation method is used, debugging optimization is carried out to the Permeability Prediction model, establishes Nonlinear Mapping relational model, and form final network model；S6: node-by-node algorithm is carried out according to the final network model, obtains the nuclear magnetic resonance T 2 spectrum permeability under pore throat character restriction on the parameters；At least one order of magnitude of the computational accuracy of tight sand permeability can be improved in the present invention.
Description
Technical field
The present invention relates to petroleum geology exploration technology and logging evaluation technical fields, and in particular to a kind of pore throat character parameter
Nuclear magnetic resonance T 2 spectrum calculation of permeability under constraint.
Background technique
During oil-gas exploration and development, usually acquisition absolute permeability (abbreviation permeability) method have core test and
The methods of calculated using well-log information, wherein the permeability of core test is the most accurate, is a kind of direct permeability acquisition side
Method is commonly used to the permeability of scale well-log information calculating.But the permeability of core test is limited by sampled point, acquisition
Permeability value is discontinuous.Therefore, well-log information permeameter is established by the relationship between research permeability and log response parameter
Model is calculated, permeability is calculated using well-log information and just seems very necessary.
There are two kinds using the thinking that Conventional Logs obtain permeability, one is establish permeability and various reservoirs ginseng
Number, the statistical model between log response parameter, for example, permeability and porosity, natural gamma relative value statistical model, point
Penetration rate model, the model of neural computing permeability etc. that flow unit is established；Another method is to pass through rock physics
Relationship between model foundation permeability and reservoir characteristic parameter, such as Wyllie-Rose permeability log interpretation model,
Timur permeability log interpretation model etc..Two kinds of methods complement one another, but Wyllie-Rose equation, Timur equation are hypotonic
In saturating reservoir permeability logging evaluation, the permeability ratio of permeability and the core test analysis of calculating is compared with error is larger.With equivalent
Based on rock constituents are theoretical, similarity principle is migrated according to charge migration and fluid molecule, on effective conductive porosity basis
On, the permeability log interpretation model between effective flowing porosity and permeability is established, model accuracy is better than Wyllie-Rose
Equation and Timur equation.
With going deep into for domestic oil-gas exploration and development, Low-porosity And Low Permeability Reservoir receives more concerns.Low porosity and low permeability storage
The pore throat character for collecting layer is complicated, is affected to permeability, and the above method for calculating permeability is not able to satisfy production needs.
Nuclear magnetic resonance log can measure a variety of formation informations such as porosity, pore throat character, be increasingly being applied to low porosity and low permeability
In the evaluation of permeability of reservoir.Classical nuclear magnetic resonance log penetration rate model rely primarily on nuclear magnetic resonance with relaxation behavior and
Diffusion property, it is established that statistical relational expression.By analyzing the correlation of its parameter with permeability, final permeability is obtained
Model simultaneously calculates continuous permeability curve with it.There is Coates model with the classical model that nuclear magnetic resonance log calculates permeability
With SDR model.Coates model establishes model using movable fluid, constraint fluid, nuclear-magnetism porosity；SDR model utilizes nuclear-magnetism
Porosity, T2 geometric mean establish model.
But application effect of the classical NMR Permeability Models in Low-porosity And Low Permeability Reservoir computing permeability is not very managed
Think, main cause is that low porosity and low permeability reservoir heterogeneity is serious, and relevant parameter cannot sufficiently reflect irreducible water saturation, pore throat
The features such as structure；Pore throat character parameter pair is weakened again using the set a song to music calculation of permeability of line morphology of nuclear magnetic resonance log T2
The important function of computing permeability.
Summary of the invention
The technical problem to be solved in the present invention is that the nuclear magnetic resonance T 2 spectrum under providing a kind of pore throat character restriction on the parameters seeps
Saturating rate calculation method.
The technical solution adopted by the present invention to solve the technical problems is: the nuclear-magnetism under a kind of pore throat character restriction on the parameters is total
The T2 that shakes composes calculation of permeability, comprising the following steps:
S1: pseudo- capillary pressure curve is converted by nuclear magnetic resonance T 2 spectrum curve；
S2: pore throat character parameter, including replacement pressure, maximum pore throat half are calculated according to the pseudo- capillary pressure curve
Diameter, hole throat median radius, hole throat mean value, primary flow bore larynx radius average value and sorting coefficient；
S3: it is counted after converting pseudo- capillary pressure curve using the core analysis permeability data and nuclear magnetic resonance T 2 spectrum of acquisition
Obtained pore throat character parameter, NMR logging data establishes pore throat character parameter, nuclear magnetic resonance T 2 spectrum-permeability pair
Answer relation data body (Him, T2ij, Pij, ki) (m=1 ..., l；I=1 ..., n), wherein H refers to pore throat character parameter value,
T2 refers to nuclear magnetic resonance T 2 spectrum time component Discrete value, and P refers to the discrete value of T2 spectral amplitude component, and k refers to that permeability value, i refer to acquisition
Core analysis permeability sample spot serial number, m refer to pore throat character parameter serial number, and j refers to component serial number, and l is each depth point hole Adam's apple
The number of structure parameter, n are the core test permeability sample spot number of acquisition；
S4: selected non-linear map establishes the nuclear magnetic resonance T 2 spectrum Permeability Prediction mould under pore throat character restriction on the parameters
Type；
S5: cross validation method is used, to the nuclear magnetic resonance T 2 spectrum Permeability Prediction under the pore throat character restriction on the parameters
Model carries out debugging optimization, establishes Nonlinear Mapping relational model, and forms final network model；
S6: node-by-node algorithm is carried out according to the final network model, obtains the nuclear magnetic resonance under pore throat character restriction on the parameters
T2 composes permeability.
Preferably, in step sl, pseudo- capillary pressure curve is converted by nuclear magnetic resonance T 2 spectrum curve, including walked as follows
It is rapid:
Step S11: calculating scale conversion between nuclear magnetic resonance T 2 spectrum and capillary pressure curve according to similarity principle is
Number C is converted by formula PC=C/T2 linear relationship and is obtained pseudo- capillary pressure curve；Or
Step S12: it is calculated between nuclear magnetic resonance T 2 spectrum and capillary pressure curve laterally according to two dimension segmentation equal-area method
Calibration factor and longitudinal calibration factor in large and small aperture, are calculated pseudo- capillary pressure curve.
Preferably, in step s 2,
The replacement pressure refers to that nonwetting phase initially enters the pressure of rock sample maximum venturi, and the calculation method used is from most
Low pressure point starts, more adjacent pressure spot into mercury saturation degree, starting point of the two into mercury saturation degree difference greater than 1% is row
Drive pressure；
The maximum pore throat radius refers to the corresponding pore throat radius of the replacement pressure；
The hole throat median radius refers to the corresponding pore throat radius value of 50% mercury saturation degree；
The pore constriction mean value is to describe the mean place of experimental data value, that is, indicates the average bit of full throat distribution
It sets；
The primary flow bore larynx radius average value refer to permeability contribution margin it is accumulative up to 95% when pore throat radius it is average
Value；
The sorting coefficient is the parameter for reflecting Pore throat size uniform level.
Preferably, in step s3, the core test permeability data and pore throat character parameter, nuclear magnetic resonance log money
Material is after carrying out core Location by permeability measured by core experiment room, with nuclear magnetic resonance T 2 spectrum and nuclear-magnetism T2 spectrum conversion puppet
The corresponding relationship for the pore throat character parameter being calculated after capillary pressure curve.
Preferably, in step s 4, the non-linear map uses BP neural network method.
Preferably, for the BP neural network method using single hidden layer network, input layer unit number is the nuclear magnetic resonance T2
Spectral component discrete value n, output layer unit number are 1, and One hidden layer neuron number is calculated with following formula and obtained: n1=(n+1)1/2+a,
Wherein, a is the constant between 1~10.
Preferably, in step s 5, the cross validation method is in given modeling sample, and most of sample is used for
Nonlinear Mapping relational model is established, stays fraction sample for model testing, and seek the prediction error of the fraction sample
Quadratic sum.
Preferably, the sample that most of sample is 70%, the sample that the fraction sample is 30%.
Preferably, during establishing the Nonlinear Mapping relational model, using iterative algorithm constantly to described 70%
Sample calculated；After each iteration is completed, the data of 15% sample cross verifying are assessed；When described 70%
Sample computational accuracy when being difficult to improve, error minimum is calculated using the data that wherein 15% sample carries out cross validation
When corresponding network model be final network model, and using the data of last 15% test to the final network model progress
Test.
Preferably, in step s 6, node-by-node algorithm is carried out to full well section using the final network model, obtains hole Adam's apple
Nuclear magnetic resonance T 2 spectrum permeability under structure restriction on the parameters.
Implement the invention has the following advantages: the nuclear-magnetism under a kind of pore throat character restriction on the parameters proposed by the present invention is total
The T2 that shakes composes calculation of permeability, solves and weakens pore throat using the calculation of permeability of nuclear magnetic resonance T 2 spectrum tracing pattern
The problem that structural parameters cause the permeability precision for thinking poorly of hole low permeability reservoir relatively low, it had both considered nuclear magnetic resonance T 2 spectrum and includes
All pore structural informations, also strengthen effect of the pore throat character parameter in computing permeability.And it does not need to carry out a large amount of
Rock core nuclear magnetic resonance experiment, and have higher precision, and be easily achieved, can satisfy production requirement.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples, in attached drawing:
Fig. 1 is the flow chart of the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters of the present invention；
Fig. 2 is that the present invention utilizes nuclear magnetic resonance T 2 spectrum, the nuclear magnetic resonance T 2 spectrum of 30 sampled points to be converted into capillary pressure
The penetration rate model precision schematic diagram of median radius and core analysis the permeability building calculated after curve；
Fig. 3 is the nuclear magnetic resonance T 2 spectrum calculation of permeability well logging money under pore throat character restriction on the parameters proposed by the present invention
Material processing achievement schematic diagram.
Specific embodiment
As shown in Figure 1, being the nuclear magnetic resonance T 2 spectrum computing permeability side under a kind of pore throat character restriction on the parameters of the invention
Method causes to evaluate for solving to weaken pore throat character parameter using the calculation of permeability of nuclear magnetic resonance T 2 spectrum tracing pattern
The relatively low problem of the permeability precision of low porosity and low permeability reservoir, comprising the following steps:
S1: pseudo- capillary pressure curve is converted by nuclear magnetic resonance T 2 spectrum curve；In this step, by nuclear magnetic resonance T 2 spectrum
Curve is converted into pseudo- capillary pressure curve, including following method:
Step S11: calculating scale conversion between nuclear magnetic resonance T 2 spectrum and capillary pressure curve according to similarity principle is
Number C is converted by formula PC=C/T2 linear relationship and is obtained pseudo- capillary pressure curve；Or
Step S12: it is calculated between nuclear magnetic resonance T 2 spectrum and capillary pressure curve laterally according to two dimension segmentation equal-area method
Calibration factor and longitudinal calibration factor in large and small aperture, are calculated pseudo- capillary pressure curve.
S2: pore throat character parameter, including replacement pressure, maximum pore throat half are calculated according to the pseudo- capillary pressure curve
Diameter, hole throat median radius, hole throat mean value, primary flow bore larynx radius average value and sorting coefficient；In step s 2,
The replacement pressure refers to that nonwetting phase initially enters the pressure of rock sample maximum venturi, and the calculation method used is from most
Low pressure point starts, more adjacent pressure spot into mercury saturation degree, starting point of the two into mercury saturation degree difference greater than 1% is row
Drive pressure；
The maximum pore throat radius refers to the corresponding pore throat radius of the replacement pressure；
The hole throat median radius refers to the corresponding pore throat radius value of 50% mercury saturation degree；
The pore constriction mean value is to describe the mean place of experimental data value, that is, indicates the average bit of full throat distribution
It sets；
The primary flow bore larynx radius average value refer to permeability contribution margin it is accumulative up to 95% when pore throat radius it is average
Value；
The sorting coefficient is the parameter for reflecting Pore throat size uniform level.
In the present embodiment, pore constriction mean value is calculated using following formula (1), (2), (3), primary flow bore larynx radius is averaged
Value and sorting coefficient
R in formulaM, Rz, Sp be respectively pore constriction mean value, primary flow bore larynx radius average value and sorting coefficient；ΔS,i
For the mercury saturation degree of the i-th interval measure；riFor throat radius corresponding with Δ Si.
S3: it is calculated after converting pseudo- capillary pressure curve using core analysis permeability data and nuclear-magnetism the T2 spectrum of acquisition
Pore throat character parameter, the NMR logging data arrived establishes the corresponding pass of pore throat character parameter, nuclear magnetic resonance T 2 spectrum-permeability
It is data volume (Him, T2ij, Pij, ki) (m=1 ..., l；I=1 ..., n), this data volume are that known core analysis is seeped
The nuclear magnetic resonance T 2 spectrum data of same depth point and turned by the depth point nuclear magnetic resonance T 2 spectrum after saturating rate data and core Location
Change the pore throat character parameter composition calculated after capillary pressure curve into.Wherein, H refers to that pore throat character parameter value, T2 refer to that nuclear-magnetism is total
Component Discrete value between T2 time spectrum of shaking, P refer to the discrete value of T2 spectral amplitude component, and k refers to that permeability value, i refer to that the core test of acquisition seeps
Saturating rate sample spot serial number, m refer to pore throat character parameter serial number, and j refers to component serial number, and l is of each depth point pore throat character parameter
Number, n are the core test permeability sample spot number of acquisition.
It is to be appreciated that in step s3, the core analysis permeability data and pore throat character parameter, nuclear magnetic resonance are surveyed
Well data is to turn after carrying out core Location by permeability measured by core experiment room with nuclear magnetic resonance T 2 spectrum and nuclear-magnetism T2 spectrum
Change the corresponding relationship for the pore throat character parameter being calculated after pseudo- capillary pressure curve.
S4: selected non-linear map establishes the nuclear magnetic resonance T 2 spectrum Permeability Prediction mould under pore throat character restriction on the parameters
Type；In the present embodiment, in step s 4, the non-linear map uses BP neural network method, further, the BP
For neural network method using single hidden layer network, input layer unit number is the nuclear magnetic resonance T 2 spectrum component Discrete value n, output layer list
First number is 1, and One hidden layer neuron number is calculated with following formula and obtained: n1=(n+1)1/2+ a, wherein a is normal between 1~10
Number.
It is to be appreciated that establishing Permeability Prediction model based on BP neural network non-linear map.BP neural network
Method is a kind of Multi-layered Feedforward Networks by the training of error backpropagation algorithm, using steepest descent method, by backpropagation come
The constantly weight and threshold value of adjustment network keeps the error sum of squares of network minimum.Since neural network is too sensitive, it is necessary to big
Sample is modeled the model that can just obtain with stronger generalization ability, because the method requires the quantity of sample more as far as possible, i.e. n value
It is big as far as possible, since BP neural network carries out global optimizing using gradient descent algorithm, it is contemplated that gradient descent method is easily trapped into
The problem of local optimum, (seeps input (pore throat character parameter and sampled point nuclear magnetic resonance T 2 spectrum) and output when establishing model
Saturating rate) it is normalized, to improve the generalization ability of model.Such as 3 pore throat character parameters, 30 sampled point nuclear magnetic resonance T2
33 × 10 × 1 network structure can be used to be modeled for spectrum.
S5: cross validation method is used, to the nuclear magnetic resonance T 2 spectrum Permeability Prediction under the pore throat character restriction on the parameters
Model carries out debugging optimization, establishes Nonlinear Mapping relational model, and forms final network model；It is to be appreciated that using handing over
Verification method is pitched, model debugging is carried out, to obtain best result.
In step s 5, the cross validation method is in given modeling sample, and most of sample is non-for establishing
Linear mapping relation model stays fraction sample for model testing, and seeks the Prediction sum squares of the fraction sample.
To a preferred embodiment, the sample that the major part sample is 70%, the sample that the fraction sample is 30%.
During establishing the Nonlinear Mapping relational model, using iterative algorithm constantly to described 70% sample into
Row calculates；After each iteration is completed, the data of 15% sample cross verifying are assessed；When described 70% sample meter
Precision is calculated when being difficult to improve, it is corresponding when being calculated error minimum using the data that wherein 15% sample carries out cross validation
Network model is final network model, and is tested using the data of last 15% test the final network model.
S6: node-by-node algorithm is carried out according to the final network model, obtains the nuclear magnetic resonance under pore throat character restriction on the parameters
T2 composes permeability.
Specifically, in step s 6, full well section is counted point by point using the final network model obtained in step s 5
It calculates, obtains the nuclear magnetic resonance T 2 spectrum permeability under pore throat character restriction on the parameters.
Further, model involved in step S5 and step S6 of the present invention, is all with the Permeability Prediction in step S4
Based on model, debugging optimization is carried out by cross validation method etc., forms final network model, utilizes the final network model
Node-by-node algorithm is carried out to full well section, the nuclear magnetic resonance T 2 spectrum permeability under pore throat character restriction on the parameters can be obtained.
The present invention passes through the nuclear-magnetism using same depth point after the core analysis permeability data and core Location acquired
The T2 that resonates composes data, and is converted by the depth point nuclear magnetic resonance T 2 spectrum pore throat character number calculated after capillary pressure curve
According to establishing pore throat character parameter, nuclear magnetic resonance T 2 spectrum and permeability corresponding relation data body.Using the BP nerve of None-linear approximation
Network algorithm, and optimal pore throat character, T2 are established by the optimization of the normalized of parameter, cross validation and algorithm parameter
The Nonlinear Mapping of fractions distribution and permeability obtains the nuclear magnetic resonance T 2 spectrum tracing pattern in the case where pore throat character constrains and calculates infiltration
The model of rate realizes reservoir permeability Continuous plus.The present invention compared with traditional nuclear magnetic resonance logging data calculates penetration rate model,
At least one order of magnitude of the computational accuracy of tight sand permeability can be improved.
Fig. 2 is that the present invention utilizes nuclear magnetic resonance T 2 spectrum, the nuclear magnetic resonance T 2 spectrum of 30 sampled points to be converted into capillary pressure
The penetration rate model precision schematic diagram of median radius and core analysis the permeability building calculated after curve.Wherein, data volume
151 data are shared, for the picture left above to return the precision sentenced after 70% sample, top right plot is the number of 15% cross validation in figure
According to lower-left figure is the data of 15% test, and bottom-right graph is the prediction result of last all sample points.From the results, it was seen that building
Mode coherence coefficients have reached R=0.95, and cross validation precision has reached R=0.69, and measuring accuracy has reached R=0.78.It is found that
The model accuracy is higher, and the permeability value of calculating is more reliable.
As shown in figure 3, being the nuclear magnetic resonance T 2 spectrum computing permeability side under pore throat character restriction on the parameters proposed by the present invention
Method Well Data Processing achievement, wherein first is gamma ray curve in figure, and second is depth, and third road is total for nuclear-magnetism
The T2 that shakes composes form, and the 4th scatterplot is core analysis permeability, and the 4th curve is the permeability value that the present invention calculates.It can be bright
Aobvious to find out, calculating permeability using the nuclear magnetic resonance T 2 spectrum of pore throat character restriction on the parameters is that effectively, it does not both need to carry out big
The rock core nuclear magnetic resonance experiment of amount, and have higher precision, and be easily achieved, it can satisfy production requirement.
Implement the invention has the following advantages: the nuclear-magnetism under a kind of pore throat character restriction on the parameters proposed by the present invention is total
The T2 that shakes composes calculation of permeability, solves and is weakened using the set a song to music calculation of permeability of line morphology of nuclear magnetic resonance log T2
The problem that pore throat character parameter causes the permeability precision for thinking poorly of hole low permeability reservoir relatively low, it both considers nuclear magnetic resonance T 2 spectrum
All pore structural informations for including also strengthen effect of the pore throat character parameter in computing permeability.And it does not need to carry out
A large amount of rock core nuclear magnetic resonance experiment, and have higher precision, and be easily achieved, it can satisfy production requirement.
It should be understood that above embodiments only express the preferred embodiment of the present invention, description is more specific and detailed
Carefully, but it cannot be understood as limitations on the scope of the patent of the present invention；It should be pointed out that for the common skill of this field
For art personnel, without departing from the inventive concept of the premise, above-mentioned technical characterstic can be freely combined, can also be done
Several modifications and improvements out, these are all within the scope of protection of the present invention；Therefore, all to be done with scope of the invention as claimed
Equivalents and modification, should belong to the covering scope of the claims in the present invention.
Claims (10)
1. the nuclear magnetic resonance T 2 spectrum calculation of permeability under a kind of pore throat character restriction on the parameters, which is characterized in that including following
Step:
S1: pseudo- capillary pressure curve is converted by nuclear magnetic resonance T 2 spectrum curve；
S2: pore throat character parameter, including replacement pressure, maximum pore throat radius, hole are calculated according to the pseudo- capillary pressure curve
Gap throat median radius, hole throat mean value, primary flow bore larynx radius average value and sorting coefficient；
S3: it is calculated after converting pseudo- capillary pressure curve using the core analysis permeability data and nuclear magnetic resonance T 2 spectrum of acquisition
Pore throat character parameter, the NMR logging data arrived establishes the corresponding pass of pore throat character parameter, nuclear magnetic resonance T 2 spectrum-permeability
It is data volume (Him, T2ij, Pij, ki) (m=1 ..., l；I=1 ..., n), wherein H refers to that pore throat character parameter value, T2 refer to
Nuclear magnetic resonance T 2 spectrum time component Discrete value, P refer to the discrete value of T2 spectral amplitude component, and k refers to that permeability value, i refer to the rock core of acquisition
Permeability sample spot serial number is analyzed, m refers to pore throat character parameter serial number, and j refers to component serial number, and l is each depth point pore throat character ginseng
Several numbers, n are the core test permeability sample spot number of acquisition；
S4: selected non-linear map establishes the nuclear magnetic resonance T 2 spectrum Permeability Prediction model under pore throat character restriction on the parameters；
S5: cross validation method is used, to the nuclear magnetic resonance T 2 spectrum Permeability Prediction model under the pore throat character restriction on the parameters
Debugging optimization is carried out, Nonlinear Mapping relational model is established, and forms final network model；
S6: node-by-node algorithm is carried out according to the final network model, obtains the nuclear magnetic resonance T 2 spectrum under pore throat character restriction on the parameters
Permeability.
2. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 1, special
Sign is, in step sl, converts pseudo- capillary pressure curve for nuclear magnetic resonance T 2 spectrum curve, includes the following steps:
Step S11: calculating scale conversion coefficient C between nuclear magnetic resonance T 2 spectrum and capillary pressure curve according to similarity principle,
It is converted by formula PC=C/T2 linear relationship and obtains pseudo- capillary pressure curve；Or
Step S12: lateral scale between nuclear magnetic resonance T 2 spectrum and capillary pressure curve is calculated according to two dimension segmentation equal-area method
Coefficient and longitudinal calibration factor in large and small aperture, are calculated pseudo- capillary pressure curve.
3. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 1, special
Sign is, in step s 2,
The replacement pressure refers to that nonwetting phase initially enters the pressure of rock sample maximum venturi, and the calculation method used is from minimal pressure
Force starts, more adjacent pressure spot into mercury saturation degree, starting point of the two into mercury saturation degree difference greater than 1% is that row drives pressure
Power；
The maximum pore throat radius refers to the corresponding pore throat radius of the replacement pressure；
The hole throat median radius refers to the corresponding pore throat radius value of 50% mercury saturation degree；
The pore constriction mean value is to describe the mean place of experimental data value, that is, indicates the mean place of full throat distribution；
The primary flow bore larynx radius average value refer to permeability contribution margin it is accumulative up to 95% when pore throat radius average value；
The sorting coefficient is the parameter for reflecting Pore throat size uniform level.
4. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 1, special
Sign is that in step s3, the core analysis permeability data and pore throat character parameter, NMR logging data are to pass through
After permeability measured by core experiment room carries out core Location, pseudo- capillary is converted with nuclear magnetic resonance T 2 spectrum and nuclear magnetic resonance T 2 spectrum
The corresponding relationship for the pore throat character parameter being calculated after pipe pressure curve.
5. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 1, special
Sign is that in step s 4, the non-linear map uses BP neural network method.
6. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 5, special
Sign is, the BP neural network method using single hidden layer network, input layer unit number be the nuclear magnetic resonance T 2 spectrum component from
Value n is dissipated, output layer unit number is 1, and One hidden layer neuron number is calculated with following formula and obtained: n1=(n+1)1/2+ a, wherein a is
Constant between 1~10.
7. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 1, special
Sign is that in step s 5, the cross validation method is in given modeling sample, and most of sample is non-thread for establishing
Property mapping relations model, stay fraction sample for model testing, and seek the Prediction sum squares of the fraction sample.
8. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 7, special
Sign is, the sample that the major part sample is 70%, the sample that the fraction sample is 30%.
9. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 8, special
Sign is, during establishing the Nonlinear Mapping relational model, using iterative algorithm constantly to described 70% sample into
Row calculates；After each iteration is completed, the data of 15% sample cross verifying are assessed；When described 70% sample meter
Precision is calculated when being difficult to improve, it is corresponding when being calculated error minimum using the data that wherein 15% sample carries out cross validation
Network model is final network model, and is tested using the data of last 15% test the final network model.
10. the nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters according to claim 9, special
Sign is, in step s 6, carries out node-by-node algorithm to full well section using the final network model, obtains pore throat character parameter about
Nuclear magnetic resonance T 2 spectrum permeability under beam.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811574964.1A CN109655394B (en) | 2018-12-21 | 2018-12-21 | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811574964.1A CN109655394B (en) | 2018-12-21 | 2018-12-21 | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109655394A true CN109655394A (en) | 2019-04-19 |
| CN109655394B CN109655394B (en) | 2022-05-27 |
Family
ID=66115655
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811574964.1A Active CN109655394B (en) | 2018-12-21 | 2018-12-21 | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109655394B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111950193A (en) * | 2020-07-15 | 2020-11-17 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| CN112966365A (en) * | 2021-02-04 | 2021-06-15 | 中海石油(中国)有限公司 | Ultra-low permeability condensate gas reservoir retrograde condensation injury evaluation method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2272997A1 (en) * | 1997-03-11 | 1998-09-17 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived with nmr |
| CN102253069A (en) * | 2011-04-20 | 2011-11-23 | 中国石油天然气股份有限公司 | Method and device for determining permeability according to nuclear magnetic resonance T2 spectrum |
| CN103279647A (en) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | Method for calculating tight sandstone reservoir permeability based on hole feature parameters |
| CN103884633A (en) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | Method and device for confirming rock permeability |
| CN108181219A (en) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | A kind of method and its device based on carbonate porosity structure prediction permeability |
| CN108288092A (en) * | 2018-01-09 | 2018-07-17 | 中国海洋石油集团有限公司 | A method of obtaining tight sand permeability using nuclear magnetic resonance T 2 spectrum form |
-
2018
- 2018-12-21 CN CN201811574964.1A patent/CN109655394B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2272997A1 (en) * | 1997-03-11 | 1998-09-17 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived with nmr |
| CN102253069A (en) * | 2011-04-20 | 2011-11-23 | 中国石油天然气股份有限公司 | Method and device for determining permeability according to nuclear magnetic resonance T2 spectrum |
| CN103279647A (en) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | Method for calculating tight sandstone reservoir permeability based on hole feature parameters |
| CN103884633A (en) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | Method and device for confirming rock permeability |
| CN108181219A (en) * | 2017-11-15 | 2018-06-19 | 中国石油天然气股份有限公司 | A kind of method and its device based on carbonate porosity structure prediction permeability |
| CN108288092A (en) * | 2018-01-09 | 2018-07-17 | 中国海洋石油集团有限公司 | A method of obtaining tight sand permeability using nuclear magnetic resonance T 2 spectrum form |
Non-Patent Citations (2)
| Title |
|---|
| 中国石油天然气总公司劳资局: "《油层物性实验工》", 30 November 1997, 石油工业出版社 * |
| 成志刚等: "基于储层孔喉特征参数计算致密砂岩渗透率的新方法", 《测井技术》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111950193A (en) * | 2020-07-15 | 2020-11-17 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| WO2022011893A1 (en) * | 2020-07-15 | 2022-01-20 | 中海油田服务股份有限公司 | Reservoir-based modeling method and device for pore network model |
| CN111950193B (en) * | 2020-07-15 | 2023-03-14 | 中海油田服务股份有限公司 | Modeling method and device of pore network model based on reservoir |
| CN112966365A (en) * | 2021-02-04 | 2021-06-15 | 中海石油(中国)有限公司 | Ultra-low permeability condensate gas reservoir retrograde condensation injury evaluation method |
| CN112966365B (en) * | 2021-02-04 | 2023-09-12 | 中海石油(中国)有限公司 | Method for evaluating reverse condensation injury of ultralow condensation gas reservoir |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109655394B (en) | 2022-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bauer et al. | Sources of uncertainty in shear stress and roughness length estimates derived from velocity profiles | |
| Halford et al. | User guide for the drawdown-limited, multi-node well (MNW) package for the US Geological Survey's modular three-dimensional finite-difference ground-water flow model, versions MODFLOW-96 and MODFLOW-2000 | |
| Li et al. | Three‐dimensional geostatistical inversion of flowmeter and pumping test data | |
| CN105182424B (en) | A kind of method and apparatus based on patchy saturation quantitative forecast reservoir porosity | |
| CN108303752A (en) | Glutenite effective reservoir conventional logging quantitative identification method | |
| CN105426612A (en) | Stratum component optimization determining method and device | |
| CN106154351A (en) | A kind of evaluation method of low porosity permeability reservoir permeability | |
| CN108596780B (en) | A kind of Reservoir type division methods and system based on multiple information | |
| CN104698492B (en) | A kind of method for calculating abnormal formation pressure | |
| CN106223942A (en) | A kind of Conglomerate Reservoir shale content computational methods based on Well logging curve reconstruction | |
| CN109509111A (en) | The prediction technique and system of prospect pit strata pressure | |
| CN104727813A (en) | Porosity measuring method for resource storage commensalism type stratums | |
| CN104932027A (en) | Reservoir classification method based on nuclear magnetic resonance logging | |
| CN107808032A (en) | Shale gas reservoir water saturation computational methods and system | |
| CN109655394A (en) | A kind of nuclear magnetic resonance T 2 spectrum calculation of permeability under pore throat character restriction on the parameters | |
| CN111812736A (en) | Method for evaluating gas content of compact sandstone anhydrous gas reservoir | |
| CN105844011B (en) | A kind of calculation of permeability based on capillary model | |
| CN107762483A (en) | A kind of Fluid Identification Method of coefficient correlation and envelope size based on log | |
| Zhou et al. | Determination of pore size distribution in tight gas sandstones based on Bayesian regularization neural network with MICP, NMR and petrophysical logs | |
| CN116859478B (en) | Groundwater simulation method and system based on transient electromagnetic imaging | |
| Sureshjani et al. | Estimating reservoir permeability distribution from analysis of pressure/rate transient data: A regional approach | |
| CN103698494B (en) | Method and device for determining saturation degree of hydrocarbon in lithologic trap | |
| CN112069444B (en) | Method and computer for calculating reservoir well testing permeability by using well logging data | |
| CN111965720B (en) | Method for acquiring hydraulic conductivity based on ground-well combination | |
| Hadibeik et al. | Petrophysical properties of unconventional low-mobility reservoirs (shale gas and heavy oil) by using newly developed adaptive testing approach |
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 |