CN110159263A - Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index - Google Patents
Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index Download PDFInfo
- Publication number
- CN110159263A CN110159263A CN201910460416.4A CN201910460416A CN110159263A CN 110159263 A CN110159263 A CN 110159263A CN 201910460416 A CN201910460416 A CN 201910460416A CN 110159263 A CN110159263 A CN 110159263A
- Authority
- CN
- China
- Prior art keywords
- gas
- well logging
- index
- logging
- surveys
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000012530 fluid Substances 0.000 title claims abstract description 19
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 17
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- NNPPMTNAJDCUHE-UHFFFAOYSA-N isobutane Chemical compound CC(C)C NNPPMTNAJDCUHE-UHFFFAOYSA-N 0.000 claims abstract description 13
- QWTDNUCVQCZILF-UHFFFAOYSA-N isopentane Chemical compound CCC(C)C QWTDNUCVQCZILF-UHFFFAOYSA-N 0.000 claims abstract description 11
- 239000001294 propane Substances 0.000 claims abstract description 9
- 239000001282 iso-butane Substances 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims description 13
- 241000209094 Oryza Species 0.000 claims description 12
- 235000007164 Oryza sativa Nutrition 0.000 claims description 12
- 235000009566 rice Nutrition 0.000 claims description 12
- 235000013847 iso-butane Nutrition 0.000 claims description 6
- AFABGHUZZDYHJO-UHFFFAOYSA-N dimethyl butane Natural products CCCC(C)C AFABGHUZZDYHJO-UHFFFAOYSA-N 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 4
- 239000013589 supplement Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 9
- 238000005259 measurement Methods 0.000 abstract description 4
- 150000001335 aliphatic alkanes Chemical class 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- OFBQJSOFQDEBGM-UHFFFAOYSA-N Pentane Chemical compound CCCCC OFBQJSOFQDEBGM-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012764 semi-quantitative analysis Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index, which comprises the following steps: a. respectively calculating ethane C collected by logging gas detector2Propane C3Isobutane iC4Isopentane iC5The gas power spectrum value of (1); b. superposing the gas measurement power spectrums obtained by calculation into a graph, and respectively calculating a superposed reverse sequence oil index and a superposed water index; c. and setting a discrimination mode to determine the oil content of the reservoir section and whether water is produced. The invention extends the traditional logging gas logging analysis method of a depth domain-time domain to a frequency domain, and establishes a reverse sequence superposition method based on the characteristic parameters of logging gas logging components according to the quantity relation among alkane components to judge the oil content and the water content.
Description
Technical field
The present invention relates to a kind of well logging gas to survey the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum, belongs to
The measurement of electro-mechanical arts industrial gas oil exploration and development scope log data, processing technology field.
Background technique
Currently, Oil Layer Reservoir flow net model is good at NORTH CHINA oil field and engineering services side, and well logging solution
It releases at present mainly based on scatterplot-crossplot method, also visual parameter curve method, is limited to above-mentioned analysis processing method
Characteristic, well logging Gas Logging Value is explained to determine to rely on and visually observes and manually sentence knowledge, subjective based on quantitative and semi-quantitative,
Accuracy of identification and accuracy rate be not high, and traditional logging explanation method is often " to have logging presentation just to analyze, have logging presentation just to solve
Release ", drop ply, misjudgement phenomenon are commonplace, need to research and develop it is a kind of can be important for log data is analyzed and is explained
Specific well section carries out the judgment method and technology of global analysis.
Retrieval discovery, the search result for having no " well logging Gas Logging Value " and " oil reservoir " while occurring utilize well logging Gas Logging Value
The patented technology for carrying out oil layer recognition is less, and " well logging Gas Logging Value " keyword is then visible " to be based on drilling fluid logging parameter Gas Logging Value
Logging explanation method (application number: 201810556205.6, the applying date: 2018-05-31, river celebrate the limited public affairs of drilling project
Department) ", " carry out the logging explanation method (application number: 201810553427.2 that reservoir sentences knowledge using triangle plate and probability graph road
The applying date: 2018-05-31, Chuan Qing drilling project Co., Ltd) " etc., above-mentioned means of interpretation be based on when m- Depth Domain benefit
Oil-gas Layer fluid identification is carried out with chart method, from frequency-domain analysis and does not utilize Gas Logging Value using spectrum analysis means;In addition, closing
Key word " well logging " and " Reservoir interpretation " then visible " gas detection logging Reservoir interpretation method (application based on oil reservoir quantitative judge plate
Number: 201510524973.X, the applying date: 2015-08-25, Tianjin Branch, China Sea Petroleum (China) have no right-refute
The state of returning) ", which states the explanation that reservoir fluid property identification is established according to the characteristic component ratio of well logging Gas Logging Value
Plate is the simple computation based on Depth Domain.Patent database search key " reservoir fluid ", entry is more, representative
Patent Case be " a kind of oil-gas-water layer well logging starry sky explanation evaluating method (and application number: 201710426025.1, the applying date: 2017-
06-08, China petroleum Group Pohai drilling Engineering Co., Ltd) ", which proposes a kind of utilization gas survey
Component " starry sky " means of interpretation carries out the scheme of oil-gas-water layer flow net model, tentatively solves point, line intersects, overlapping is blocked
The problems such as, integrated application well logging Gas Logging Value and geochemical logging TPI parameter etc., although having certain effect, but still in Depth Domain
It calculates and analyzes.To sum up, the oil reservoir well logging gas reservoir prediction method of mainstream at present, establish when m- Depth Domain on the basis of,
It usually only analyzes and high abnormal " logging presentation " layer of well logging gas survey occurs, there is dependence subjectivity to sentence knowledge, be easy to be limited to visually to see
Examine with subjective judgement and the risk of drop ply, erroneous judgement.
In addition, application number 201610962713.5, publication date 2017-02-15, disclose a kind of based on brill multi-parameter point
The oil gas water identification method and system of analysis, but be disadvantageous in that existing for it: (1), data source " color is described as according to its patent
Spectrum signal " and " chromatogram " do not limit mode and pattern (data format) that chromatographic signal obtains, patent claims
By " chromatographic signal acquisition chromatogram ", the processing of patent subsequent step and process object are chromatogram, the element of processing for description
Similarity of maps, figure correlation between chromatogram, data source are biased to instrument close beta and debugging.(2), tend to
Form comparison and descriptive analysis, main purpose are in Fu for carry out peak value or peak area to the chromatographic peak that chromatograph inside obtains
Leaf transformation calculates, and obtains frequency domain value, then calculates crosspower spectrum, emphasize chromatography to be calculated or chromatography to be explained for criterion numeral
According to correlation analysis, have a disadvantage in that the determination to normal data and specification do not elaborate in detail in the claims, give
The invention of people's one kind is the tentative imagination for being difficult to operate realization rather than practicability or the practical technique by experimental examination.
(3), its core concept is the similitude comparison of spectrogram, but setting for standard spectrogram and establishes clearly proposing,
The various for spectrogram situation and scatter plot of coefficient R are only provided, for the auxiliary judgment of resistivity curve in well logging etc.
There is no detailed standard and explanation, more generally.
Summary of the invention
It is an object of the invention to overcome the above problem of the existing technology, it is folded to provide a kind of well logging gas survey component frequency spectrum
Close the method that backward index determines properties of fluid in bearing stratum.Conventional depth domain-time-domain well logging gas is surveyed analysis method by the present invention
Frequency domain is expanded to, and establishes the backward overlapping side for surveying the characteristic parameter of component based on well logging gas according to quantitative relation between alkane component
Method differentiates oiliness and the property of water-bearing.
To achieve the above object, The technical solution adopted by the invention is as follows:
A kind of method that well logging gas surveys component frequency spectrum overlapping backward index judgement properties of fluid in bearing stratum, includes the following steps:
A, the collected ethane C of gas analyser for well logging is calculated separately2, propane C3, iso-butane iC4, isopentane iC5Gas measurement of power
Rate spectrum;
B, the gas power scale being calculated spectrum is overlapped into figure, calculates separately overlapping backward oil index and aqua index;
C, setting discrimination model determines Reservoir Section oiliness and whether produces water.
The step a is specifically included:
(1), it data prediction: adjusts well logging gas to be processed and surveys 2 required by parameter to spectrum analysisx, x ∈ Z+A data
Point, and trending processing and decentralization processing are carried out to data adjusted;
(2), spectra calculation: carrying out Fast Fourier Transform (FFT) for pretreated well logging gas measured data, obtains corresponding
The mould of sequence of complex numbers is calculated in sequence of complex numbers, draws power spectrum chart according to the mould of sequence of complex numbers.
Adjustment well logging gas to be processed survey parameter to set a data point include: set well logging gas survey parameter to be processed as
Ci, wherein i=1,2,3 ..., m, if m=2x+ y, and y is the minimum value enumerated in expression formula, as y=0, is directly surveyed to gas
Argument sequence is continued to employ, and as y < 0, then supplements 0 value of y item to argument sequence, as y > 0, then removes y data with existing point, will
Well logging gas to be processed is surveyed parameter and is adjusted to m=2xA data point.
With ethane C2For, it is described that trending is gone to handle are as follows: Ci=Ci-Ci-1;Decentralization processing are as follows: Ci=Ci-μc;Its
In, μcC is surveyed for well section gas to be processediThe average value of component, CiFor the Gas Logging Value ethane component of corresponding i-th of depth point.
The spectra calculation includes:
The pretreated well logging gas measured data that (1) is obtained carries out Fast Fourier Transform (FFT), obtains corresponding plural sequence
Column, if well logging gas measured data time series to be processed is x (t), wherein t=0,1,2,3 ..., n, n ∈ Z+If Fu to be asked
In leaf transformation complex result be X (k), formula is as follows:
The mould P (k) of calculated complex sequence, in addition, akFor the real part of P (k), bkFor the imaginary part of P (k), algorithm is as follows:
The step b is specifically included:
The gas surveys component C2、C3、iC4、iC5Corresponding power spectrum chart overlapped under same abscissa and present it is certain poor
Value, chooses C respectively3、C2Difference calculates overlapping backward oil index Δ Po, choose iC5、iC4Difference calculates overlapping backward aqua index Δ
Pw, wherein calculating formula is as follows:
Overlap backward oil index Δ PoCalculation is as follows:
Overlap backward aqua index Δ PwCalculation is as follows:
Wherein, j ∈ Z+, n is Reservoir Section thickness to be explained, max (f (C2)) it is Reservoir Section C to be explained2Power spectrum maximum number
Value, min (f (Cq)) it is Reservoir Section iC to be explained4Component power spectrum greatest measure, f (C3j) it is that jth rice corresponds to C3Power intensity values,
f(C2j) it is that jth rice corresponds to C2Power intensity values, f (Cpj) it is that jth rice corresponds to iC5Power intensity values, f (Cqj) it is that jth rice is corresponding
iC4Power intensity values.
The step c is specifically included:
Overlap backward oil index Δ PoWith overlapping backward aqua index Δ PwOil-gas possibility method of discrimination about reservoir to be explained
It is as follows:
Reservoir not oil-producing:
ΔPo=0
Oil reservoir (not aqueous):
ΔPo> 0
ΔPw< 0
Water layer:
ΔPo< 0
Or Δ Pw> 0.
Compared with prior art, using the present invention has the advantages that
(1) data acquisition source of the present invention is more: data source of the present invention is in including well logging chromatograph but being not limited to other
The collected gas measured data of analytical equipment, data format are that 1 depth point uniquely corresponds to the one-to-one " deep of 1 data
Degree evidence " is discrete data, and consistent with the chromatographic apparatus common data form of petroleum industry field application, has industry general
All over property and refering in particular to property, in addition, data source of the present invention can be not only restricted to the property of acquisition equipment, all gas that can collect surveys number
Inspection is surveyed according to and with well logging gas that Gas Logging Value measuring device that installation controlled condition described in invention or controlled factordisease influence obtains
In measured data is all contained in.
(2) trend data analysis of the present invention: the present invention is by carrying out the deletion of 2n data to discrete well logging Gas Logging Value data
Or cover, then carry out decentralization processing, then using Fast Fourier Transform (FFT) FFT by when m- Depth Domain on it is discrete
Gas measured data is transformed into frequency domain and obtains power spectrum chart;On this basis, calculating well logging Gas Logging Value is representative and boundary is anticipated
The propane C of justice3With ethane C2Frequency spectrum cumulative difference value (the more C of carbon number of component3The C less relative to carbon number2Increase
Amount), and the judgement section of oiliness and the property of water-bearing is judged in setting, the Land use systems of spectrogram are further specified, are given
The quantitative well logging Gas Logging Value of sxemiquantitative-judges the solution of reservoir fluid property, has operability.
(3) the technology of the present invention key point is that power spectrum chart establishes gentle survey characteristic value overlapping backward (oil, water) index calculating
Unlike the prior art with the differentiation of reservoir fluid property threshold, and there is continuity and specific aim.
(4) component characteristics parameter propane C is surveyed for gas for the first time3, ethane C2, iso-butane iC4And isopentane iC5Carry out power
Spectrum is calculated and is analyzed according to spectrogram folding method, and total hydrocarbon TG is divided into light component for boundary essentially according to propane by existing method
It is compared and analyzes with heavy constituent, m- Depth Domain when being confined to, analysis method relies on the naked eyes subjectivity such as plate and curve method and sentences
Knowledge method differentiates that difficulty is larger, is difficult to relationship between surveying component based on gas and carries out quantitative-semi-quantitative analysis;The present invention is directed to gas for the first time
Survey component characteristics parameter propane C3, ethane C2, iso-butane iC4And isopentane iC5It is calculated and is analyzed, added by spectrometer
The capacity volume variance of row power spectrum compares, and from frequency domain angle analysis well logging Gas Logging Value, and then quantitative-sxemiquantitative judges reservoir oiliness
And the property of water-bearing, and discrimination model is intuitive compared with conventional method, succinct.
Detailed description of the invention
Fig. 1 is that spectrogram backward of the present invention overlaps difference schematic diagram.
Specific embodiment
A kind of numerical value for well logging Gas Logging Value data of system of the present invention and descriptive means of interpretation.
Implementation step decomposes as follows:
A, the collected ethane C of gas analyser for well logging is calculated separately2, propane C3, iso-butane iC4, isopentane iC5Gas measurement of power
Rate spectrum;
B, the gas power scale being calculated spectrum is overlapped into figure, calculates separately overlapping backward oil index and aqua index;
C, setting discrimination model determines Reservoir Section oiliness and whether produces water.
The step a is specifically included:
(1), it data prediction: adjusts well logging gas to be processed and surveys 2 required by parameter to spectrum analysisx, x ∈ Z+A data
Point, and trending processing and decentralization processing are carried out to data adjusted;
(2), spectra calculation: carrying out Fast Fourier Transform (FFT) for pretreated well logging gas measured data, obtains corresponding
The mould of sequence of complex numbers is calculated in sequence of complex numbers, draws power spectrum chart according to the mould of sequence of complex numbers.
Adjustment well logging gas to be processed survey parameter to set a data point include: set well logging gas survey parameter to be processed as
Ci, wherein i=1,2,3 ..., m, if m=2x+ y, and y is the minimum value enumerated in expression formula, as y=0, is directly surveyed to gas
Argument sequence is continued to employ, and as y < 0, then supplements 0 value of y item to argument sequence, as y > 0, then removes y data with existing point, will
Well logging gas to be processed is surveyed parameter and is adjusted to m=2xA data point.
With ethane C2For, it is described that trending is gone to handle are as follows: Ci=Ci-Ci-1;Decentralization processing are as follows: Ci=Ci-μc;Its
In, μcC is surveyed for well section gas to be processediThe average value of component, CiFor the Gas Logging Value ethane component of corresponding i-th of depth point.
The spectra calculation includes:
Pretreated well logging gas measured data is subjected to Fast Fourier Transform (FFT), corresponding sequence of complex numbers is obtained, if wait locate
The well logging gas measured data time series of reason is x (t), wherein t=0,1,2,3 ..., n, n ∈ Z+If Fourier transformation to be asked is multiple
Number result is X (k), and formula is as follows:
The mould P (k) of calculated complex sequence, in addition, akFor the real part of P (k), bkFor the imaginary part of P (k), algorithm is as follows:
Ordinate is set as power P, and abscissa is set as frequency f, draws ethane C respectively2, propane C3, iso-butane iC4, it is different
Pentane iC5Gas power scale spectrogram.
2) spectrogram overlaps, and calculates overlapping backward index
The gas surveys component C2、C3、iC4、iC5Corresponding power spectrum chart overlapped under same abscissa and present it is certain poor
Value, chooses C respectively3、C2Difference calculates overlapping backward oil index Δ Po, choose iC5、iC4Difference calculates overlapping backward aqua index Δ
Pw, wherein calculating formula is as follows:
Overlap backward oil index Δ PoCalculation is as follows:
Overlap backward aqua index Δ PwCalculation is as follows:
Wherein, j ∈ Z+, n is Reservoir Section thickness to be explained, max (f (C2)) it is Reservoir Section C to be explained2Power spectrum maximum number
Value, min (f (Cq)) it is Reservoir Section iC to be explained4Component power spectrum greatest measure, f (C3j) it is that jth rice corresponds to C3Power intensity values,
f(C2j) it is that jth rice corresponds to C2Power intensity values, f (Cpj) it is that jth rice corresponds to iC5Power intensity values, f (Cqj) it is that jth rice is corresponding
iC4Power intensity values.
3) Reservoir Section oiliness is determined according to discrimination model and whether produces water
Overlap backward oil index Δ PoWith overlapping backward aqua index Δ PwOil-gas possibility method of discrimination about reservoir to be explained
It is as follows:
Reservoir not oil-producing:
ΔPo=0
Oil reservoir (not aqueous):
ΔPo> 0
ΔPw< 0
Water layer:
ΔPo< 0
Or Δ Pw> 0.
Claims (7)
1. a kind of well logging gas surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum, it is characterised in that: including
Following steps:
A, the collected ethane C of gas analyser for well logging is calculated separately2, propane C3, iso-butane iC4, isopentane iC5Gas power scale spectrum
Value;
B, the gas power scale being calculated spectrum is overlapped into figure, calculates separately overlapping backward oil index and aqua index;
C, setting discrimination model determines Reservoir Section oiliness and whether produces water.
2. well logging gas according to claim 1 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: the step a is specifically included:
(1), it data prediction: adjusts well logging gas to be processed and surveys 2 required by parameter to spectrum analysisx, x ∈ Z+A data point, and
Trending processing and decentralization processing are carried out to data adjusted;
(2), spectra calculation: pretreated well logging gas measured data is subjected to Fast Fourier Transform (FFT), obtains corresponding plural number
The mould of sequence of complex numbers is calculated in sequence, draws power spectrum chart according to the mould of sequence of complex numbers.
3. well logging gas according to claim 2 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: it includes: to set well logging gas to be processed to survey parameter that the adjustment well logging gas to be processed, which surveys parameter to a data point is set,
For Ci, wherein i=1,2,3 ..., m, if m=2x+ y, and y is the minimum value enumerated in expression formula, as y=0, directly to gas
It surveys argument sequence to continue to employ, as y < 0, then supplements 0 value of y item to argument sequence, as y > 0, then remove y data with existing point,
Well logging gas to be processed survey parameter is adjusted to m=2xA data point.
4. well logging gas according to claim 3 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: with ethane C2For, it is described that trending is gone to handle are as follows: Ci=Ci-Ci-1;Decentralization processing are as follows: Ci=Ci-μc;
Wherein, μcC is surveyed for well section gas to be processediThe average value of component, CiFor the Gas Logging Value ethane component of corresponding i-th of depth point.
5. well logging gas according to claim 4 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: the spectra calculation includes:
The pretreated well logging gas measured data that (1) is obtained carries out Fast Fourier Transform (FFT), obtains corresponding sequence of complex numbers, if
Well logging gas measured data time series to be processed is x (t), wherein t=0,1,2,3 ..., n, n ∈ Z+If Fourier to be asked becomes
Changing complex result is X (k), and formula is as follows:
The mould P (k) of calculated complex sequence, in addition, akFor the real part of P (k), bkFor the imaginary part of P (k), algorithm is as follows:
6. well logging gas according to claim 5 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: the step b is specifically included:
The gas surveys component C2、C3、iC4、iC5Corresponding power spectrum chart overlaps under same abscissa and presents certain difference, point
It Xuan Qu not C3、C2Difference calculates overlapping backward oil index Δ Po, choose iC5、iC4Difference calculates overlapping backward aqua index Δ Pw,
In, calculating formula is as follows:
Overlap backward oil index Δ PoCalculation is as follows:
Overlap backward aqua index Δ PwCalculation is as follows:
Wherein, j ∈ Z+, n is Reservoir Section thickness to be explained, max (f (C2)) it is Reservoir Section C to be explained2Power spectrum greatest measure,
min(f(Cq)) it is Reservoir Section iC to be explained4Component power spectrum greatest measure, f (C3j) it is that jth rice corresponds to C3Power intensity values, f
(C2j) it is that jth rice corresponds to C2Power intensity values, f (Cpj) it is that jth rice corresponds to iC5Power intensity values, f (Cqj) it is that jth rice corresponds to iC4
Power intensity values.
7. well logging gas according to claim 6 surveys the method that component frequency spectrum overlapping backward index determines properties of fluid in bearing stratum,
It is characterized by: the step c is specifically included:
Overlap backward oil index Δ PoWith overlapping backward aqua index Δ PwOil-gas possibility method of discrimination about reservoir to be explained is such as
Under:
Reservoir not oil-producing:
ΔPo=0
Oil reservoir (not aqueous):
ΔPo> 0
ΔPw< 0
Water layer:
ΔPo< 0
Or Δ Pw> 0.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910460416.4A CN110159263B (en) | 2019-05-30 | 2019-05-30 | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910460416.4A CN110159263B (en) | 2019-05-30 | 2019-05-30 | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110159263A true CN110159263A (en) | 2019-08-23 |
CN110159263B CN110159263B (en) | 2022-09-02 |
Family
ID=67630004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910460416.4A Active CN110159263B (en) | 2019-05-30 | 2019-05-30 | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110159263B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113159978A (en) * | 2021-05-07 | 2021-07-23 | 中国石油天然气集团有限公司 | Automatic analysis method for logging gas logging data |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040054478A1 (en) * | 2002-08-30 | 2004-03-18 | Brygynevych Valentyn Oleksandrovych | Method of wave diagnostics of the oil-and-gas-deposit |
US20060082367A1 (en) * | 2004-10-18 | 2006-04-20 | Edwards Carl M | Nuclear magnetic resonance method and apparatus for evaluating a characteristic of a region |
US20100027377A1 (en) * | 2006-02-24 | 2010-02-04 | Hannes Georges Zuercher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
CN102900433A (en) * | 2012-10-31 | 2013-01-30 | 中国石油集团川庆钻探工程有限公司 | Method for judging reservoir fluid by establishing gas logging chart through discriminant analysis |
CN104865614A (en) * | 2014-02-20 | 2015-08-26 | 中国石油化工股份有限公司 | Complicated reservoir fluid identification method based on variable skeleton parameter |
CN105804734A (en) * | 2016-04-12 | 2016-07-27 | 中国石油天然气集团公司 | Method for identifying thickened oil reservoir by utilizing nuclear magnetic resonance well logging |
CN106555587A (en) * | 2016-11-24 | 2017-04-05 | 中国石油集团渤海钻探工程有限公司 | A kind of gas surveys the difference profit recognition methods of hydrocarbon component longitudinal direction |
CN108798661A (en) * | 2018-06-11 | 2018-11-13 | 中国石油集团川庆钻探工程有限公司 | Method for identifying oil well reservoir and fluid-containing property by using logging gas logging component parameters |
CN109403958A (en) * | 2018-08-31 | 2019-03-01 | 中国石油集团川庆钻探工程有限公司 | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index |
CN109403960A (en) * | 2018-09-06 | 2019-03-01 | 中国石油集团川庆钻探工程有限公司 | Method for judging reservoir fluid properties by logging gas peak-logging state |
CN109798108A (en) * | 2019-03-06 | 2019-05-24 | 中国石油大学(北京) | The recognition methods of Deep tight reservoir-compact reservoir complex fluid property and device |
-
2019
- 2019-05-30 CN CN201910460416.4A patent/CN110159263B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040054478A1 (en) * | 2002-08-30 | 2004-03-18 | Brygynevych Valentyn Oleksandrovych | Method of wave diagnostics of the oil-and-gas-deposit |
US20060082367A1 (en) * | 2004-10-18 | 2006-04-20 | Edwards Carl M | Nuclear magnetic resonance method and apparatus for evaluating a characteristic of a region |
US20100027377A1 (en) * | 2006-02-24 | 2010-02-04 | Hannes Georges Zuercher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
CN102900433A (en) * | 2012-10-31 | 2013-01-30 | 中国石油集团川庆钻探工程有限公司 | Method for judging reservoir fluid by establishing gas logging chart through discriminant analysis |
CN104865614A (en) * | 2014-02-20 | 2015-08-26 | 中国石油化工股份有限公司 | Complicated reservoir fluid identification method based on variable skeleton parameter |
CN105804734A (en) * | 2016-04-12 | 2016-07-27 | 中国石油天然气集团公司 | Method for identifying thickened oil reservoir by utilizing nuclear magnetic resonance well logging |
CN106555587A (en) * | 2016-11-24 | 2017-04-05 | 中国石油集团渤海钻探工程有限公司 | A kind of gas surveys the difference profit recognition methods of hydrocarbon component longitudinal direction |
CN108798661A (en) * | 2018-06-11 | 2018-11-13 | 中国石油集团川庆钻探工程有限公司 | Method for identifying oil well reservoir and fluid-containing property by using logging gas logging component parameters |
CN109403958A (en) * | 2018-08-31 | 2019-03-01 | 中国石油集团川庆钻探工程有限公司 | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index |
CN109403960A (en) * | 2018-09-06 | 2019-03-01 | 中国石油集团川庆钻探工程有限公司 | Method for judging reservoir fluid properties by logging gas peak-logging state |
CN109798108A (en) * | 2019-03-06 | 2019-05-24 | 中国石油大学(北京) | The recognition methods of Deep tight reservoir-compact reservoir complex fluid property and device |
Non-Patent Citations (6)
Title |
---|
PING YIN等: "Application of element logging to lithologic identification of key horizons in Sichuan–Chongqing gas provinces", 《NATURAL GAS INDUSTRY》 * |
RUYUE WANG等: "Logging identification for the Lower Cambrian Niutitang shale reservoir in the Upper Yangtze region, China: A case study of the Cengong block, Guizhou Province", 《JOURNAL OF NATURAL GAS GEOSCIENCE》 * |
彭文春等: "气测数据校正在录井综合解释中的应用", 《录井工程》 * |
李战奎等: "气测录井储集层流体性质识别的新方法——正-异烷烃比值图板法", 《录井工程》 * |
赵卫平等: "最大熵谱分解技术在A油田薄砂储层厚度预测中的应用", 《工程地球物理学报》 * |
邵东波等: ""三低"高含水气藏录井解释评价方法", 《录井工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113159978A (en) * | 2021-05-07 | 2021-07-23 | 中国石油天然气集团有限公司 | Automatic analysis method for logging gas logging data |
CN113159978B (en) * | 2021-05-07 | 2024-06-07 | 中国石油天然气集团有限公司 | Automatic analysis method for logging gas logging data |
Also Published As
Publication number | Publication date |
---|---|
CN110159263B (en) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101936973B (en) | Method for rapidly classifying hydrocarbon oil with combined gas-phase chromatography-mass spectrometryer | |
CN110056348B (en) | Method and system for determining formation fluid composition and properties | |
CN106401579B (en) | Oil-gas-water identification method and system based on while-drilling multi-parameter analysis | |
CN107422381B (en) | A kind of earthquake low-frequency information fluid prediction method based on EEMD-ICA | |
NO20151108A1 (en) | Method to perform rapid formation fluid analysis | |
CN104142519A (en) | Mud rock crack oil deposit predicting method | |
CN108798661B (en) | Method for identifying oil well reservoir and fluid-containing property by using logging gas logging component parameters | |
CN105134190A (en) | Gas logging oil layer interpreting method based on oil layer quantitative recognition layout | |
CN107389648B (en) | A kind of identification of 3-D quantitative fluorescent characteristics peak and fluid type judgment method | |
CN109557578A (en) | A kind of reservoir gas-bearing property detection method and device | |
CN106447511A (en) | Heavy oil-based quantitative identification chart interpretation method with gas logging | |
Eltaleb et al. | Analysis of fracture injection tests using signal processing approach | |
CN111502650B (en) | Method for identifying condensate gas layer by using gas measurement derived parameters and application thereof | |
Yang et al. | Application of well logs integration and wavelet transform to improve fracture zones detection in metamorphic rocks | |
Masoudi et al. | Developing a method for identification of net zones using log data and diffusivity equation | |
CN104142517A (en) | Hydrocarbon detection method utilizing seismic data dynamic spectrum properties | |
CN109403958A (en) | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index | |
CN110159263A (en) | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index | |
CN110414169A (en) | A kind of fourier infrared gas detection logging method and device thereof | |
CN111257933B (en) | Novel method for predicting oil and gas reservoir based on low-frequency shadow phenomenon | |
CN102330553B (en) | MDT test dynamic spectrum fluid identification method | |
CN110017136B (en) | Water flooded layer identification and water production rate prediction method based on apparent water layer resistivity | |
Boulassel et al. | A new multifractal analysis-based for identifying the reservoir fluid nature | |
US8838393B2 (en) | Systems and methods for topographic analysis | |
CN113051518B (en) | Quantitative calculation method and device for gas-oil ratio, computer equipment and storage medium |
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 |