CN109403958B - Method for identifying formation fluid properties by using gas logging power spectrum equivalent index - Google Patents
Method for identifying formation fluid properties by using gas logging power spectrum equivalent index Download PDFInfo
- Publication number
- CN109403958B CN109403958B CN201811010807.8A CN201811010807A CN109403958B CN 109403958 B CN109403958 B CN 109403958B CN 201811010807 A CN201811010807 A CN 201811010807A CN 109403958 B CN109403958 B CN 109403958B
- Authority
- CN
- China
- Prior art keywords
- power
- peak
- gas
- logging
- power 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.)
- Active
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
Abstract
The invention discloses a method for identifying the property of formation fluid by using equivalent index of gas logging power spectrum, which comprises the following steps: a. calculating a power spectrum, and establishing the power spectrum according to the calculation result; b. partitioning the established power spectrogram; c. calculating equivalent identification indexes of the logging gas logging power spectrum according to the partitioned power spectrogram; d. and determining the reservoir fluid property according to the equivalent identification index according to a judgment rule. The method only carries out single-line operation on the aspects of gas measurement data time-depth-to-frequency domain conversion, power spectrogram partition, power spectrum equivalent identification index calculation and fluid judgment rules based on the equivalent identification index, and has the advantages of continuity and irreversibility in steps, and strong operability and reproducibility.
Description
Technical Field
The invention relates to a method for identifying the property of formation fluid by using equivalent index of gas logging power spectrum, belonging to the technical field of logging data measurement and processing in the exploration and development category of petroleum and natural gas industry in the electromechanical field.
Background
At present, for the logging gas measurement industry is generally limited in a time-depth domain, a visual observation method is adopted for overall trend observation, or data processing and analysis are carried out on a single point, and no effective and practical technical means is disclosed or patented in the field of overall analysis of a large well section. The more traditional processing ideas in the past are: "there is an abnormal response or the logging gas shows that it is processed locally, and it is not processed for the section without abnormality".
In 2015, a detailed method (version 1.0) for analyzing and processing logging gas measurement values by using a spectrum analysis technology and a meeting paper of an implementation process are published in the activity of gathering documents of the statement of the logging engineering in 2016 to the whole analysis technology of a large well section for logging gas logging started by the Chongqing drilling geological exploration and development institute, a time-series method-based gas logging interpretation method is published in 2016, an improved version (version 2.0) of spectrum analysis and fluid interpretation is formally disclosed in the report of the project conclusion of the Mizhongqing oil and gas field exploration institution in 2016 to describe the geological foundation and the implementation process of spectrum analysis implementation and application in more detail, but a fluid property identification part only provides a plurality of perceptual modes, so that the dependence on manual identification is high, and the method cannot be completely put into production and application.
In 3 months in 2017, a paper "new method for identifying gas reservoir gas-water layer by logging parameter in Longwanggio group gas reservoir gas-water layer in Anyue gas field Mixi zone" published by Chongwanggio institute for Chongqing drilling geological exploration and development in Chuan details application and practice conditions of processing logging gas logging data in a high-wear structure in Chuan Yu gas zone by a spectrum analysis method, but contents such as fluid property identification are not published. It is noted that shanghai Shenkai corporation focusing on the development and production of chromatographs in 2017 published a method and system for identifying oil, gas and water based on multi-parameter analysis while drilling, patent application No. 201610962713.5, published time: 2017-02-15, the invention patent proposes that data are subjected to spectrum analysis in a chromatograph, a standard spectrogram needs to be established for comparison and analysis, and the core is correlation coefficient calculation and correlation degree judgment. In practice, it is inconvenient and difficult to implement for the production application unit of logging gas logging, so it is considered to use other simpler and more general methods to quickly process and judge the logging gas measurement value, and then identify the properties of the drilling reservoir fluid.
Disclosure of Invention
The invention aims to overcome the problems in the prior art and provides a method for identifying the property of formation fluid by using an equivalent index of a gas logging power spectrum. The method only carries out single-line operation on the aspects of gas measurement data time-depth-to-frequency domain conversion, power spectrogram partition, power spectrum equivalent identification index calculation and fluid judgment rules based on the equivalent identification index, and has the advantages of continuity and irreversibility in steps, and strong operability and reproducibility.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a method for identifying the properties of formation fluid by using equivalent indexes of gas logging power spectrums is characterized by comprising the following steps:
a. calculating a power spectrum, and establishing the power spectrum according to the calculation result;
b. partitioning the established power spectrogram;
c. calculating equivalent identification indexes of the logging gas logging power spectrum according to the partitioned power spectrogram;
d. and determining the reservoir fluid property according to the equivalent identification index according to a judgment rule.
The step a specifically comprises the following steps:
a1, preparing data;
a2, fast Fourier transform and power spectrum calculation: logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence;
a3, establishing a power spectrogram: calculating a2 to obtain a plurality of series of modules, setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrum.
The a1 data preparation:
setting the gas logging parameters to be processed to form a time sequence XtWherein t is 1,2,3, …, m is 2x+ y and y is the minimum value listed in the above expression, leaving it directly available for the above gasometric parameter sequence when y is 0, when y is<When 0, then supplement y 0 values to the parameter sequence, when y>When 0, the y existing data points are removed, and in sum, the logging gas logging parameters to be processed are adjusted to m 2xA data point;
the data is then detrended, i.e. the
X′t＝Xt-Xt-1
Followed by a de-centering treatment which is carried out,
wherein, muxIs X 'corresponding to the selected well section'tAverage value of (a).
The a2 fast Fourier transform and power spectrum calculation:
logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence:
through fast fourier transform, the complex number sequence is obtained as:
y＝aj+ibj。
establishing the a3 power spectrum:
calculating a2 to obtain a modulus of the complex sequence, the calculation method is: a + B ═ a ═ B' ═ B-
Setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrogram according to the following formula;
in said step b, a drill rod length range (T) is seta,Tb) As the characteristic period of the single aftereffect peak, the corresponding frequency interval (f) is inversely calculateda,fb) The power spectrogram is divided into two intervals of a display peak and a non-display peak by taking the frequency fa as the starting frequency and the demarcation frequency.
The frequency interval of the display peak is (0, f)a) The frequency interval of the non-display peak is (f)aInfinity); on the basis, the non-display peak is further divided into three intervals of a single peak, a hysteresis peak and other peaks, and the corresponding frequency intervals are respectively (f)a,fb)、(fbN), (n, ∞), where n is 3.5min/ROP (average).
In the step c, calculating an equivalent identification index figure (f) of the logging gas logging power spectrum, hereinafter referred to as f (f), in the following way:
wherein P isx′、fx' the power value between the display peak regions and the corresponding frequency value thereof are respectively; f. ofd' and fd' are the maximum power value and the corresponding frequency of a single peak interval respectively; sigma Px′·fx' is the product sum of all power values and corresponding frequency values in the display peak interval; sigma Pd′·fd' is the product sum of all power values and corresponding frequency values in a single peak interval; max (P)x) To display the maximum value of the power between the peak regions, fxFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksc) Is the maximum value of power between hysteresis peak regions, fcFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksd) Is the maximum value of the power between the single peak regions, fdThe frequency value corresponding to the maximum power value of the peak interval is displayed.
In the step d, reservoir fluid properties are determined by using a descriptive morphological analysis method, which specifically comprises the following steps:
describing morphological analysis, observing and judging the power value, the maximum power value and the peak area of a displayed peak and a non-displayed peak including a single peak, a delayed peak and other peaks by a power spectrogram after the partition, wherein the judgment rule is as follows:
a. gas layer: the maximum value of the power between the peak areas is over 100, and the part of the power reaches over 200; the power and the advantages between the display peak areas are obvious, and the power sum between the non-display peak areas is far lower than the display peak areas;
b. gas-water homostorey: the maximum value of the display peak power is much higher or closer to the non-display peak, but the power value is lower, usually less than 100.
c. Water layer: the power value of the display peak is lower than or not higher than 10, the total power value is smaller than that of the non-display peak, the maximum powers of the display peak and the non-display peak are relatively close, and no order difference exists;
d. dry layer: the peak power maximum is shown to be lower than the non-shown peak, or the total power value is lower than the non-shown peak, the power intensity is also lower, typically lower than 1.
In the step d, the reservoir fluid properties are determined by adopting a numerical calculation method, which specifically comprises the following steps:
numerical zoning to determine fluid properties
Establishing equivalent identification index F based on logging gas logging power spectrum(f)Discrimination pattern of numerical intervals:
gas layer: equivalent index F(f)∈[1000，∞)；
Gas-water homostorey: equivalent index F(f)∈[100，1000)；
Water layer: equivalent index F(f)∈[10，100)；
Dry layer: equivalent index F(f)∈[-∞，10)。
The invention has the advantages that:
first, principle and mode of partitioning logging gas logging data spectrogram
The prior art does not find the condition of carrying out power spectrum analysis on logging gas logging data and carrying out partition recognition on a power spectrogram according to the periodic occurrence rule of drilling hardness factors. The invention completely shows the principle, scheme and implementation steps of partitioning the logging gas power spectrogram, provides a partitioning effect graph, and has a certain effect on distinguishing and recognizing the influence factors and useful geological information according to the test application.
Second, logging gas logging power spectrum equivalent identification index calculation method
The invention provides a logging gas logging power spectrum equivalent identification index calculation method based on observation and analysis of a power spectrum, which is initiated domestically and internationally, analyzes and calculates a power spectrum partition diagram comprehensively, and represents and depicts different modes of power spectrum response.
Third, judging the implementation scheme of the gas layer, the water layer, the gas-water layer and the dry layer
The invention provides two discrimination modes of morphological analysis and numerical value partition on the basis of calculating the equivalent recognition index of the logging gas logging power spectrum, recognizes and judges the logging gas logging power spectrum from two aspects of a visual observation method and a data partition, has the capability of recognizing a gas layer, a gas-water layer, a water layer and a dry layer, utilizes and mines the gas logging data of a non-display well section which is ignored in the past, and mixes the gas logging data with the display data for comprehensive analysis, thereby improving the utilization rate of logging gas logging data. No relevant processing and analysis of logging gas measurements is known in the prior art.
Comparing with application number 201610962713.5 (hereinafter referred to as comparison document):
1. comparing the calculation methods:
(1) data source and data type are different
The data source of the comparison file is described as a ' chromatographic signal ' and a ' chromatogram ' according to the patent thereof, the mode and the style (data format) of obtaining the chromatographic signal are not limited, the patent claims describe obtaining the chromatogram ' through the chromatographic signal, the processing and processing objects of the subsequent steps of the patent are the chromatograms, the processing elements are graphic similarity and graphic correlation among the chromatograms, and the data source is biased to be tested and debugged inside the instrument.
The data of the invention is derived from gas logging data which comprises a logging chromatograph but is not limited to gas logging data acquired by other analysis equipment, the data format of the invention is 'depth data' with 1 depth point uniquely corresponding to 1 data, the data is discrete data, the data format is consistent with the common data format of a chromatographic instrument applied in the field of the petroleum industry, and the invention has industrial universality and specificity.
(2) The core algorithm is different in concept and implementation steps
The comparison document mainly aims to perform Fourier transform calculation of peak values or peak areas on chromatographic peaks acquired inside a chromatograph to obtain frequency domain values, then calculate cross-power spectrums and emphasize correlation analysis of chromatograms to be calculated or interpreted on standard data.
According to the method, 2n data deletion or bit supplement is carried out on discrete logging gas measurement data, then decentralization processing is carried out, then fast Fourier transform FFT (MICROSOFT EXCEL or SAS software is adopted in a test) is used for converting the discrete gas measurement data on a time-depth domain into a frequency domain to obtain a power spectrogram, on the basis, according to the actual situation that a drill stem is used while drilling, a drill stem length range (Ta, Tb) is set as a characteristic period of a single post-effect peak, a corresponding frequency interval (fa, fb) is inversely calculated, and the power spectrogram is divided into a display peak interval and a non-display peak interval by taking frequency fa as an initial frequency and a boundary frequency. Wherein, the frequency interval of the display peak is (0, fa), and the frequency interval of the non-display peak is (fa, infinity); on the basis, the non-display peak is further divided into three sections of a single peak, a hysteresis peak and other peaks, corresponding frequency sections are (fa, fb), (fb, n) and (n, infinity), wherein n is 3.5min/ROP (average), and the main reason of 3.5min is that 3.5min is the normal pipeline delay time, and can also be used as the boundary point of obtaining the hysteresis gas measurement value and managing the internal delay and high-frequency disturbance of a gas measuring instrument.
2. The two methods are different:
the core idea of the comparison file is the similarity comparison of spectrograms, but the setting and establishment of a standard spectrogram are not explicitly proposed in the claims, various spectrogram conditions and scatter diagrams of a correlation coefficient R are only given in a patent attached atlas, and detailed standards and descriptions are not provided for auxiliary judgment of resistivity curves and the like in logging in the claims, so that the comparison file is more general.
The invention only relates to the field of logging engineering, and the three steps of establishing and partitioning a technical key point power spectrogram, calculating an equivalent index and judging the property of a fluid are different from the three steps. The method only carries out single-line operation on the aspects of gas measurement data time-depth-to-frequency domain conversion, power spectrogram partition, power spectrum equivalent identification index calculation and a fluid judgment rule based on the equivalent identification index, the steps have coherence and irreversibility, the power spectrum partition mainly depends on a time sequence periodicity analysis idea on the basis of geological experience analysis of the gas measurement data periodicity fluctuation reason in the drilling process, the power spectrum equivalent identification method mainly carries out calculation and analysis based on total power comparison and comparison ideas of all partitions, and finally the fluid property judgment rule is set.
3. Comparison of difficulty or Effect
The comparison file has the characteristics of a certain range and mixed description words, and comprises the engineering fields of logging, well logging and the like, and meanwhile, the selection, preparation and manufacturing elements of the standard characteristic spectrogram for judging the correlation of the core technology are not mentioned, and certain difficulty is brought to the implementation and reproduction.
The method does not depend on the existing sample analysis, carries out definite numerical value interval judgment according to a set judgment rule, carries out period-frequency spectrum analysis according to the gas measurement periodic fluctuation principle in the drilling process, is convenient to understand and apply, and has the characteristics of reproducibility and popularization and application.
Drawings
FIG. 1 is a schematic view of a power spectrum partition according to the present invention;
FIG. 2 is a schematic diagram of calculation of equivalent identification index of logging gas power spectrum of the present invention.
Detailed Description
Example 1
A method for identifying the properties of formation fluid by using equivalent indexes of gas logging power spectrum comprises the following steps:
a. calculating a power spectrum, and establishing the power spectrum according to the calculation result;
b. partitioning the established power spectrogram;
c. calculating equivalent identification indexes of the logging gas logging power spectrum according to the partitioned power spectrogram;
d. and determining the reservoir fluid property according to the equivalent identification index according to a judgment rule.
The step a specifically comprises the following steps:
a1, preparing data;
a2, fast Fourier transform and power spectrum calculation: logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence;
a3, establishing a power spectrogram: calculating a2 to obtain a plurality of series of modules, setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrum.
The a1 data preparation:
setting the gas logging parameters to be processed to form a time sequence XtWherein t is 1,2,3, …, m is 2x+ y and y is the minimum value listed in the above expression, leaving it directly available for the above gasometric parameter sequence when y is 0, when y is<When 0, then supplement y 0 values to the parameter sequence, when y>When 0, the y existing data points are removed, and in sum, the logging gas logging parameters to be processed are adjusted to m 2xA data point;
the data is then detrended, i.e. the
Xt′＝Xt-Xt-1
Followed by a de-centering treatment which is carried out,
wherein, muxCorresponding X for the selected well sectiont' average value.
The a2 fast Fourier transform and power spectrum calculation:
logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence:
through fast fourier transform, the complex number sequence is obtained as:
y＝aj+ibj。
establishing the a3 power spectrum:
calculating a2 to obtain a modulus of the complex sequence, the calculation method is: a + B ═ a ═ B' ═ B-
Setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrogram according to the following formula;
in said step b, a drill rod length range (T) is seta,Tb) As the characteristic period of the single aftereffect peak, the corresponding frequency interval (f) is inversely calculateda,fb) The power spectrogram is divided into two intervals of a display peak and a non-display peak by taking the frequency fa as the starting frequency and the demarcation frequency.
The frequency interval of the display peak is (0, f)a) The frequency interval of the non-display peak is (f)aInfinity); on the basis, the non-display peak is further divided into three intervals of a single peak, a hysteresis peak and other peaks, and the corresponding frequency intervals are respectively (f)a,fb)、(fbN), (n, ∞), where n is 3.5min/ROP (average).
In the step c, calculating an equivalent identification index figure (f) of the logging gas logging power spectrum, hereinafter referred to as f (f), in the following way:
wherein P isx′、fx' the power value between the display peak regions and the corresponding frequency value thereof are respectively; f. ofd' and fd' are the maximum power value and the corresponding frequency of a single peak interval respectively; sigma Px′·fx' is the product sum of all power values and corresponding frequency values in the display peak interval; sigma Pd′·fd' is the product sum of all power values and corresponding frequency values in a single peak interval; max (P)x) To display the maximum value of the power between the peak regions, fxFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksc) Is the maximum value of power between hysteresis peak regions, fcFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksd) Is the maximum value of the power between the single peak regions, fdThe frequency value corresponding to the maximum power value of the peak interval is displayed.
In the step d, reservoir fluid properties are determined by using a descriptive morphological analysis method, which specifically comprises the following steps:
describing morphological analysis, observing and judging the power value, the maximum power value and the peak area of a displayed peak and a non-displayed peak including a single peak, a delayed peak and other peaks by a power spectrogram after the partition, wherein the judgment rule is as follows:
a. gas layer: the maximum value of the power between the peak areas is over 100, and the part of the power reaches over 200; the power and the advantages between the display peak areas are obvious, and the power sum between the non-display peak areas is far lower than the display peak areas;
b. gas-water homostorey: the maximum value of the display peak power is much higher or closer to the non-display peak, but the power value is lower, usually less than 100.
c. Water layer: the power value of the display peak is lower than or not higher than 10, the total power value is smaller than that of the non-display peak, the maximum powers of the display peak and the non-display peak are relatively close, and no order difference exists;
d. dry layer: the peak power maximum is shown to be lower than the non-shown peak, or the total power value is lower than the non-shown peak, the power intensity is also lower, typically lower than 1.
In the step d, the reservoir fluid properties are determined by adopting a numerical calculation method, which specifically comprises the following steps:
numerical zoning to determine fluid properties
Establishing equivalent identification index F based on logging gas logging power spectrum(f)Discrimination pattern of numerical intervals:
gas layer: equivalent index F(f)∈[1000，∞)；
Gas-water homostorey: equivalent index F(f)∈[100，1000)；
Water layer: equivalent index F(f)∈[10，100)；
Dry layer: equivalent index F(f)∈[-∞，10)。
Example 2
The implementation steps of the invention can be decomposed into the following four items: 1) calculating a power spectrum and establishing a power spectrum; 2) partitioning a power spectrogram; 3) calculating equivalent identification indexes of the logging gas logging power spectrum; 4) and determining the reservoir fluid property according to the equivalent identification index according to a judgment rule.
1) Power spectrum calculation and power spectrum establishment
(1) Data preparation
Setting the gas logging parameters to be processed to form a time sequence XtWherein t is 1,2,3, …, m is 2x+ y and y is the minimum value listed in the above expression, leaving it directly available for the above gasometric parameter sequence when y is 0, when y is<When 0, then supplement y 0 values to the parameter sequence, when y>When 0, the y existing data points are removed, and in sum, the logging gas logging parameters to be processed are adjusted to m 2xA data point.
The data is then detrended, i.e. the
Xt′＝Xt-Xt-1
Followed by a de-centering treatment which is carried out,
wherein, muxCorresponding X for the selected well sectiont' average value.
(2) Fast Fourier transform and power spectrum calculation
The logging gas logging data array prepared by the data obtained in the step (1)
And carrying out fast Fourier transform to obtain a corresponding complex sequence.
Through fast fourier transform, the complex number sequence is obtained as:
y＝aj+ibj
(3) power spectrogram establishment
Calculating the modulus of the complex sequence obtained in the step 2, wherein the calculation method comprises the following steps: a + B ═ a ═ B' ═ B-
The power spectrum can be plotted with the ordinate set to power P and the abscissa set to frequency f.
2) Power spectrogram partitioning
Setting a drill rod length range (T)a,Tb) As the characteristic period of the single aftereffect peak, the corresponding frequency interval (f) is inversely calculateda,fb) The power spectrogram is divided into two intervals of a display peak and a non-display peak by taking the frequency fa as the starting frequency and the demarcation frequency. Wherein the frequency interval showing the peak is (0, f)a) The frequency interval of the non-display peak is (f)aInfinity); on the basis, the non-display peak is further divided into three intervals of a single peak, a hysteresis peak and other peaks, and the corresponding frequency intervals are respectively (f)a,fb)、(fbN) and (n, ∞), wherein n is 3.5min/ROP (average), and the main reason for 3.5min is that 3.5min is a typical line delay time, and may beSo as to be used as a demarcation point for obtaining a delayed gas measurement value and managing internal delay and high-frequency disturbance by a gas measuring instrument.
As shown in fig. 1, the power spectrum is thus divided into displayed peak intervals, non-displayed peak intervals (single peak interval, delayed peak interval and other peak intervals).
3) Calculating equivalent identification indexes of the logging gas logging power spectrum:
as shown in fig. 2, a core parameter of the equivalent identification method, i.e., the equivalent identification index fig. of the logging gas power spectrum, (f), is calculated as follows:
wherein P isx′、fx' the power value between the display peak regions and the corresponding frequency value thereof are respectively; f. ofd' and fd' are the maximum power value and the corresponding frequency of a single peak interval respectively; sigma Px′·fx' is the product sum of all power values and corresponding frequency values in the display peak interval; sigma Pd′·fd' is the product sum of all power values and corresponding frequency values in a single peak interval; max (P)x) To display the maximum value of the power between the peak regions, fxFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksc) Is the maximum value of power between hysteresis peak regions, fcFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksd) Is the maximum value of the power between the single peak regions, fdThe frequency value corresponding to the maximum power value of the peak interval is displayed.
4) Determining reservoir fluid properties according to decision rules:
the invention judges the fluid property according to the equivalent identification index of the logging gas power spectrum calculated in the step 3, and the judgment scheme is divided into two types: one is a descriptive morphological analysis method, and the other is a numerical calculation method.
First, describe the analysis of morphology
Describing morphological analysis, and observing and judging the power value, the maximum power value and the peak area of a display peak and a non-display peak (comprising a single peak, a hysteresis peak and other peaks) by a power spectrogram after the partition is finished. The decision rule is as follows:
a. gas layer: the maximum value of the power between the peak areas is over 100, and the part of the power reaches over 200; the power and the advantages between the display peak areas are obvious, and the power sum between the non-display peak areas is far lower than the display peak areas;
b. gas-water homostorey: the maximum value of the display peak power is much higher or closer to the non-display peak, but the power value is lower, usually less than 100.
c. Water layer: the power value of the display peak is lower than or not higher than 10, the total power value is smaller than that of the non-display peak, the maximum powers of the display peak and the non-display peak are relatively close, and no order difference exists;
d. dry layer: the peak power maximum is shown to be lower than the non-shown peak, or the total power value is lower than the non-shown peak, the power intensity is also lower, typically lower than 1.
Two, numerical value partition determination of fluid properties
Establishing a logging-based gas logging power spectrum equivalent identification index F according to the response characteristics of partial wells to the power spectrogram(f)Discrimination pattern of numerical intervals:
gas layer: equivalent index F(f)∈[1000，∞)；
Gas-water homostorey: equivalent index F(f)∈[100，1000)；
Water layer: equivalent index F(f)∈[10，100)；
Dry layer: equivalent index F(f)∈[-∞，10)。
Claims (7)
1. A method for identifying the properties of formation fluid by using equivalent indexes of gas logging power spectrums is characterized by comprising the following steps:
a. calculating a power spectrum, and establishing the power spectrum according to the calculation result;
b. partitioning the established power spectrogram;
c. calculating equivalent identification indexes of the logging gas logging power spectrum according to the partitioned power spectrogram;
d. determining the property of the reservoir fluid according to the equivalent identification index according to a judgment rule;
in the step c, calculating an equivalent identification index figure (f) of the logging gas logging power spectrum, hereinafter referred to as f (f), in the following way:
wherein P isx′、fx' the power value between the display peak regions and the corresponding frequency value thereof are respectively; pd' and fd' power values between single peak regions and corresponding frequencies thereof, respectively; sigma Px′·fx' is the product sum of all power values and corresponding frequency values in the display peak interval; sigma Pd′·fd' is the product sum of all power values and corresponding frequency values in a single peak interval; max (P)x) To display the maximum value of the power between the peak regions, fxFor indicating the frequency value, Max (P), corresponding to the maximum value of the power between the peaksc) Is the maximum value of power between hysteresis peak regions, fcFor the frequency corresponding to the maximum value of the power in the lag peak interval, Max (P)d) Is the maximum value of the power between the single peak regions, fdThe frequency corresponds to the maximum value of the power among the single peak areas;
in the step d, reservoir fluid properties are determined by using a descriptive morphological analysis method, which specifically comprises the following steps:
describing morphological analysis, observing and judging the power value, the maximum power value and the peak area of a displayed peak and a non-displayed peak including a single peak, a delayed peak and other peaks by a power spectrogram after the partition, wherein the judgment rule is as follows:
(1) and a gas layer: the maximum value of the power between the peak areas is over 100, and the part of the power reaches over 200; the power and the advantages between the display peak areas are obvious, and the power sum between the non-display peak areas is far lower than the display peak areas;
(2) gas-water layer: the maximum value of the display peak power is far higher or more close to the non-display peak, but the power value is lower, usually less than 100;
(3) and a water layer: the power value of the display peak is lower than or not higher than 10, the total power value is smaller than that of the non-display peak, the maximum powers of the display peak and the non-display peak are relatively close, and no order difference exists;
(4) and a dry layer: the maximum value of the display peak power is lower than that of the non-display peak, or the total power value is lower than that of the non-display peak, and the power intensity is also lower and is usually lower than 1;
in the step d, the reservoir fluid properties are determined by adopting a numerical calculation method, which specifically comprises the following steps:
numerical zoning to determine fluid properties
Establishing equivalent identification index F based on logging gas logging power spectrum(f)Discrimination pattern of numerical intervals:
gas layer: equivalent index F(f)∈[1000，∞)；
Gas-water homostorey: equivalent index F(f)∈[100，1000)；
Water layer: equivalent index F(f)∈[10，100)；
Dry layer: equivalent index F(f)∈[-∞，10)。
2. The method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 1, wherein: the step a specifically comprises the following steps:
a1, preparing data;
a2, fast Fourier transform and power spectrum calculation: logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence;
a3, establishing a power spectrogram: calculating a2 to obtain a plurality of series of modules, setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrum.
3. The method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 2, wherein: the a1 data preparation:
setting the gas logging parameters to be processed to form a time sequence XtWherein t is 1,2,3, …, m is 2x+ y, where y is the minimum value listed in the above expression, when y is 0, the gas logging parameter sequence is directly reserved, when y is less than 0, y 0 values are supplemented to the parameter sequence, when y is greater than 0, y existing data points are removed, and in sum, the logging gas logging parameter to be processed is adjusted to m is 2xA data point;
the data is then detrended, i.e. the
Xt′＝Xt-Xt-1
Followed by a de-centering treatment which is carried out,
X″t＝Xt′-μx
wherein, muxFor selecting corresponding X of well sectiont' average value.
4. The method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 3, wherein: the a2 fast Fourier transform and power spectrum calculation:
logging gas logging data array prepared by using data obtained by a1
Performing fast Fourier transform to obtain a corresponding complex sequence:
through fast fourier transform, the complex number sequence is obtained as:
y＝aj+ibj。
5. the method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 4, wherein: establishing the a3 power spectrum:
calculating a2 to obtain a modulus of the complex sequence, the calculation method is:
setting the ordinate as power P and the abscissa as frequency f, and drawing a power spectrogram according to the following formula;
6. the method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 5, wherein: in said step b, a drill rod length range (T) is seta，Tb) As the characteristic period of the single aftereffect peak, the corresponding frequency interval (f) is inversely calculateda，fb) The power spectrogram is divided into two intervals of a display peak and a non-display peak by taking the frequency fa as the starting frequency and the demarcation frequency.
7. The method of identifying formation fluid properties using gas log power spectrum equivalence index of claim 6, wherein: the frequency interval of the display peak is (0, f)a) The frequency interval of the non-display peak is (f)aInfinity); on the basis, the non-display peak is further divided into three intervals of a single peak, a hysteresis peak and other peaks, and the corresponding frequency intervals are respectively (f)a，fb)、(fbN), (n, ∞), where n is 3.5min/ROP, ROP being the rate of penetration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811010807.8A CN109403958B (en) | 2018-08-31 | 2018-08-31 | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811010807.8A CN109403958B (en) | 2018-08-31 | 2018-08-31 | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109403958A CN109403958A (en) | 2019-03-01 |
| CN109403958B true CN109403958B (en) | 2021-10-22 |
Family
ID=65463689
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811010807.8A Active CN109403958B (en) | 2018-08-31 | 2018-08-31 | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109403958B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110159263B (en) * | 2019-05-30 | 2022-09-02 | 中国石油集团川庆钻探工程有限公司 | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index |
| CN110725679B (en) * | 2019-09-12 | 2022-10-11 | 中国石油集团川庆钻探工程有限公司 | Method for judging fracturing modification potential of unconventional oil and gas reservoir by utilizing rock stratum fracture index |
| CN113159978A (en) * | 2021-05-07 | 2021-07-23 | 中国石油天然气集团有限公司 | Automatic analysis method for logging gas logging data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102169188A (en) * | 2010-12-15 | 2011-08-31 | 中国海洋石油总公司 | Method for surveying oil and gas based on Morlet spectrum |
| CN103510946A (en) * | 2012-06-19 | 2014-01-15 | 中国石油集团长城钻探工程有限公司 | Method for evaluating reservoir fluid property through gas logging data |
| CN104295291A (en) * | 2014-07-31 | 2015-01-21 | 中国石油集团长城钻探工程有限公司 | Method for evaluating reservoir fluid property through gasometry component ratio |
| CN105041305A (en) * | 2015-06-18 | 2015-11-11 | 中国海洋石油总公司 | Method for oil-gas-water identification based on physical property and gas logging parameter intersection |
| CN106447511A (en) * | 2016-09-09 | 2017-02-22 | 中国海洋石油总公司 | Heavy oil-based quantitative identification chart interpretation method with gas logging |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110231099A1 (en) * | 2006-12-05 | 2011-09-22 | Charles Elkins | Portable gas monitor |
-
2018
- 2018-08-31 CN CN201811010807.8A patent/CN109403958B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102169188A (en) * | 2010-12-15 | 2011-08-31 | 中国海洋石油总公司 | Method for surveying oil and gas based on Morlet spectrum |
| CN103510946A (en) * | 2012-06-19 | 2014-01-15 | 中国石油集团长城钻探工程有限公司 | Method for evaluating reservoir fluid property through gas logging data |
| CN104295291A (en) * | 2014-07-31 | 2015-01-21 | 中国石油集团长城钻探工程有限公司 | Method for evaluating reservoir fluid property through gasometry component ratio |
| CN105041305A (en) * | 2015-06-18 | 2015-11-11 | 中国海洋石油总公司 | Method for oil-gas-water identification based on physical property and gas logging parameter intersection |
| CN106447511A (en) * | 2016-09-09 | 2017-02-22 | 中国海洋石油总公司 | Heavy oil-based quantitative identification chart interpretation method with gas logging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109403958A (en) | 2019-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109403958B (en) | Method for identifying formation fluid properties by using gas logging power spectrum equivalent index | |
| US10995606B2 (en) | Well integrity analysis using sonic measurements over depth interval | |
| AU2011380514B2 (en) | Statistical reservoir model based on detected flow events | |
| US9519865B2 (en) | Apparatus and methods of analysis of pipe and annulus in a wellbore | |
| Aliev et al. | Robust technology and system for management of sucker rod pumping units in oil wells | |
| US20080162098A1 (en) | Method and apparatus for multi-dimensional data analysis to identify rock heterogeneity | |
| WO2014165487A2 (en) | Cement evaluation | |
| Yang et al. | Application of well logs integration and wavelet transform to improve fracture zones detection in metamorphic rocks | |
| Cipolla et al. | Evaluating Stimulation Effectiveness with Permanent Optical Fiber and Sealed Wellbore Pressure Monitoring: Bakken Shale Case Study | |
| CN109738955B (en) | Metamorphic rock lithology comprehensive judgment method based on component-structure classification | |
| CN110017136B (en) | Water flooded layer identification and water production rate prediction method based on apparent water layer resistivity | |
| CN109339760B (en) | Horizontal well section multi-cluster fracturing fracture number diagnosis method | |
| CN113988142B (en) | Tunnel lining cavity acoustic identification method based on convolutional neural network | |
| Mohaghegh et al. | Predicting well-stimulation results in a gas-storage field in the absence of reservoir data with neural networks | |
| CN110159263B (en) | Method for judging reservoir fluid properties by logging gas logging component spectrum superposition reverse order index | |
| CN111582528A (en) | Inter-well connectivity discrimination method based on fracture prediction and dynamic response | |
| CN111488666A (en) | Gas reservoir horizontal well productivity prediction model establishing and predicting method and device | |
| CN113775326B (en) | Method and device for evaluating movable water saturation, electronic equipment and medium | |
| CN1538166A (en) | Comprehensive detecting method for quality of well cementation | |
| CN112943224B (en) | Method for calculating dynamic liquid level of heavy oil well | |
| Alwarda et al. | Automated procedure for quantifying ISIP and friction losses from stage by stage hydraulic fracture treatment falloff data | |
| CN110761778B (en) | Fracturing fracture monitoring method and system based on sound waves | |
| CN113720952A (en) | Method, device, equipment and medium for generating chart for reservoir interpretation and evaluation | |
| Shaposhnikov et al. | Advanced techniques in integrated cement evaluation | |
| CN114818516B (en) | Intelligent prediction method for corrosion form profile of shaft |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhao Lei Inventor after: Li Xianghua Inventor after: Xu Wen Inventor before: Zhao Lei |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |