CN111060970B - Reservoir constrained deconvolution method and application - Google Patents
Reservoir constrained deconvolution method and application Download PDFInfo
- Publication number
- CN111060970B CN111060970B CN201911397428.3A CN201911397428A CN111060970B CN 111060970 B CN111060970 B CN 111060970B CN 201911397428 A CN201911397428 A CN 201911397428A CN 111060970 B CN111060970 B CN 111060970B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- wavelet
- information
- data
- processing
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Abstract
The invention provides a reservoir constrained deconvolution method, comprising the steps of: processing seismic data to obtain migration pure wave data; step two, wavelet processing is carried out on the offset pure wave data obtained in the step one, and a plurality of wavelet results are obtained; processing the reservoir information of the well logging to obtain reservoir information with the most complete layer position for serving as a basis of the wavelet obtained in the processing step two; and step four, processing the wavelet result obtained in the step two and the reservoir data obtained in the step three, and determining the wavelet reflecting the reservoir characteristics to obtain the standard wavelet. The invention has the beneficial effects that: a technical scheme based on a well drilling reservoir constraint deconvolution method is established, and the method mainly comprises the steps of obtaining reasonable wavelets by applying well logging information, enabling the wavelets after deconvolution to be consistent with wavelets of known well reservoir information, establishing deconvolution parameters by taking the wavelets as a standard, completing reservoir constraint deconvolution, and reducing errors of obtained data results and geological conditions.
Description
Technical Field
The invention relates to the technical field of reservoirs, in particular to a reservoir constrained deconvolution method and application.
Background
In the seismic exploration in the petroleum industry, the seismic data processing is to convert the seismic data collected in the field into result data capable of reflecting the information of underground strata, wherein the relevant application of the deconvolution technical method has the following characteristics: 1) the deconvolution technique is one of the key technical methods, mainly solving the complex wave problem existing in the data and improving the authenticity of the data; 2) the application process of the deconvolution technology comprises the steps of designing parameters of a plurality of correlation factor lengths, calculation step lengths and time window lengths, completing a plurality of corresponding results by a computer, selecting a 'better' parameter from the results, and applying the parameter to production to complete the deconvolution processing work.
The well control deconvolution technology is a relatively advanced technology in theory, and is mainly used for horizon calibration because more logging information cannot be utilized, reservoir information in well data is not involved, and the actual production application effect is not ideal.
At present, a commonly used deconvolution method has uncertainty in parameter selection, in actual deconvolution processing, reference standards are few, personal experience is mostly used as reference, the selected parameters have no certainty, and reservoir information in well logging data is not used, and the reservoir information is generally referred to as non-reservoir constrained deconvolution.
The deficiencies in the prior art are summarized as follows:
1. the choice of deconvolution parameters varies from person to person and lacks standardization and certainty;
2. the selection of the deconvolution parameters has no standard, only one group of parameters is selected from test results, and the subjective factors are more;
3. lack of methods for parameter validation;
4. determining parameters without a fixed flow, and having multi-solution;
5. the end result is often in error from the geological situation.
Disclosure of Invention
The invention overcomes the defects in the prior art and provides a reservoir constrained deconvolution method.
The purpose of the invention is realized by the following technical scheme.
The reservoir constrained deconvolution method and application thereof comprise the following steps:
processing seismic data to obtain migration pure wave data;
step two, wavelet processing is carried out on the offset pure wave data obtained in the step one, and a plurality of wavelet results are obtained;
processing the reservoir information of the well logging to obtain reservoir information with the most complete layer position for serving as a basis for processing the wavelets obtained in the second step;
and step four, processing the wavelet result obtained in the step two and the reservoir information obtained in the step three, and determining the wavelet reflecting reservoir characteristics to obtain the standard wavelet.
Further, the offset pure wave data of the first step is data without filtering, channel equalization and gain.
Furthermore, the wavelet processing in the second step is to perform resolution improvement, phase adjustment and denoising processing on the offset pure wave data.
Further, the method for processing the offset pure wave data in the second step specifically comprises:
step one, respectively finishing the processing of a plurality of different increments for increasing the frequency of 10-50Hz, and outputting a plurality of sets of results;
secondly, denoising the structure in the first step;
and thirdly, performing phase conversion processing on the denoising result in the second step.
Further, the method for obtaining the most complete reservoir data in the third step comprises the following steps:
firstly, obtaining logging information;
and secondly, obtaining the most complete reservoir data through the comparison analysis of the horizon by using the information of the first step.
Further, the method for obtaining the logging information in the first step comprises the following steps:
1) extracting frequency information for the interval for the reservoir section;
2) and extracting the horizon of the well data, extracting the thickness and lithology characteristics of the reservoir stratum and solving the root mean square speed of the stratum to obtain the information of the logging data.
Further, the specific method of the horizon control analysis in the second step is as follows:
1) performing seismic horizon picking and description of reservoir boundary and thickness on the result seismic data which is not subjected to wavelet processing in the step one, and calibrating stratums which are not easy to track or lack;
2) projecting the horizon and reservoir information obtained in the step 1) to data processed by different wavelets, and obtaining the most complete data of the horizon through comparison and analysis of the horizon.
Further, the step four-way passes through the data of the bandwidth, the dominant frequency, the energy and the phase to determine the standard wavelet.
Further, the method for determining the standard wavelet in the fourth step includes the following steps:
1) calibrating according to the information containing the stratum characteristics in the logging information and the result data in the second step, and determining wavelets matched with the logging information;
2) selecting wavelets meeting the stability of both continuity and traceability in a reservoir in the step 1);
3) and finally determining the wavelet which can reflect the reservoir characteristics most in the step 2) by combining the reservoir characteristics of the frequency-enhanced data according to the logging information of the reservoir.
And applying the standard wavelets obtained by the reservoir constraint deconvolution method to reservoir constraint deconvolution.
The invention has the beneficial effects that:
a technical scheme based on a well drilling reservoir constraint deconvolution method is established, and the method mainly comprises the steps of obtaining reasonable wavelets by applying well logging information, enabling the wavelets after deconvolution to be consistent with wavelets of known well reservoir information, establishing deconvolution parameters by taking the wavelets as a standard, completing reservoir constraint deconvolution, and reducing errors of obtained data results and geological conditions. The method has the advantages of reducing the influence of human factors, having good production and application effects, and solving the problems of different choice of deconvolution parameters, lack of standard property and certainty.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a schematic diagram of a process for well location calibration;
FIG. 3 is actual drilling information for an Rx well;
Detailed Description
The technical solution of the present invention is further illustrated by the following specific examples.
A method of reservoir constrained deconvolution as shown in fig. 1, comprising the steps of:
processing seismic data to obtain migration pure wave data, wherein the processing mainly completes coordinate loading of field data, removal of various noises, energy compensation, deconvolution, speed analysis, residual static correction and migration imaging processing to obtain preliminary data; the offset pure wave data which is directly provided can also be used, and the data is not added with filtering, channel equalization and gain;
step two, wavelet processing is carried out on the offset pure wave data obtained in the step one, and a plurality of wavelet results are obtained; the step is mainly used for improving resolution, phase adjustment and denoising wavelet processing, when the resolution is improved, processing of improving a plurality of different increments of 10-50Hz frequency is respectively completed, a plurality of sets of results are output, appropriate denoising is carried out, meanwhile, processing such as phase conversion can be carried out, and the strength characteristics of seismic data and well data are more consistent;
processing the reservoir information of the well logging to obtain reservoir information with the most complete layer position for serving as a basis for processing the wavelets obtained in the second step;
mainly aiming at a reservoir zone, extracting frequency information and layer position information of the zone, then completing layer position extraction of well information, extraction of reservoir thickness and lithology characteristics and solving of root mean square velocity of stratum to obtain information of standard well-passing information; selecting the seismic data which are not subjected to wavelet processing by utilizing the information, carrying out seismic horizon pickup and description of reservoir boundary and thickness, calibrating stratums which are difficult to track or lose, projecting the horizon and reservoir information onto the data which are subjected to different wavelet processing, and selecting the data with the most complete horizon as the wavelet basis of data processing through comparison and analysis of the horizon;
and step four, processing the wavelet result obtained in the step two and the reservoir data obtained in the step three, finally establishing data of frequency bandwidth, dominant frequency, energy and phase of data processing, and determining the wavelet reflecting reservoir characteristics to obtain the standard wavelet. The method comprises the following steps:
1) well position calibration: calibrating according to the information containing the stratum characteristics in the logging information and the result data in the second step, and determining wavelets matched with the logging information; the calibration process is shown in FIG. 2, where well location calibration parameters 5 are a better set of results;
2) and (3) horizon tracking: selecting wavelets meeting the stability of both continuity and traceability in a reservoir in the step 1);
3) extracting reservoir wavelet: and finally determining the wavelet which can reflect the reservoir characteristics most in the step 2) by combining the reservoir characteristics of the frequency-enhanced data according to the logging information of the reservoir.
When 2) the horizon tracking does not screen out the conforming sub-waves, the reason is that the problem of continuity is solved by adopting a phase adjustment mode in the first step and the second step of the method aiming at the lack of the horizon information, the existing horizon information is searched on another well aiming at the lack of the horizon information, then the frequency of the horizon information is filtered to obtain the required horizon information, and the discontinuous information and missing layer information adjusting methods related to the adjustment of the 2) horizon tracking are the prior art, and the invention does not technically improve the horizon tracking.
Application of standard wavelets: in actual production, standard wavelets can be directly input in a wavelet deconvolution program to finish production; and other deconvolution parameters can be used for adjustment, so that the output wavelet is close to the standard wavelet of the reservoir, and the deconvolution application of reservoir constraint is completed.
Fig. 3 is the actual drilling information for the Rx well, and fig. 3 shows that there are not so many strong wave impedance interfaces, and the results achieved by reservoir constrained deconvolution fit well with the well information.
While one embodiment of the present invention has been described in detail, the description is only a preferred embodiment of the present invention and should not be taken as limiting the scope of the invention. All equivalent changes and modifications made within the scope of the present invention shall fall within the scope of the present invention.
Claims (9)
1. A method of reservoir constrained deconvolution, characterized by: the method comprises the following steps:
processing seismic data to obtain migration pure wave data;
step two, wavelet processing is carried out on the offset pure wave data obtained in the step one, and a plurality of wavelet results are obtained;
processing the reservoir information of the well logging to obtain reservoir information with the most complete layer position for serving as a basis for processing the wavelets obtained in the second step;
step four, processing the wavelet result obtained in the step two and the reservoir data obtained in the step three, and determining a wavelet reflecting reservoir characteristics to obtain a standard wavelet;
the method for processing the offset pure wave data in the second step specifically comprises the following steps:
firstly, respectively finishing the processing of a plurality of different increments of 10-50Hz frequency and outputting a plurality of sets of results;
secondly, denoising the structure in the first step;
and thirdly, performing phase conversion processing on the denoising result obtained in the second step.
2. A method of reservoir constrained deconvolution according to claim 1, characterized by: and the offset pure wave data in the first step is data without filtering, channel equalization and gain.
3. A method of reservoir constrained deconvolution according to claim 1, characterized by: and in the second step, the wavelet processing is to perform phase adjustment and denoising processing on the offset pure wave data, so that the resolution is improved.
4. A method of reservoir constrained deconvolution according to claim 1, characterized by: the method for obtaining the most complete reservoir data in the third step comprises the following steps:
firstly, obtaining logging information;
and secondly, obtaining the most complete reservoir data through the comparison analysis of the horizon by using the logging data information of the first step.
5. A method of reservoir constrained deconvolution according to claim 4, characterized by: the method for acquiring the logging information in the first step comprises the following steps:
1) extracting frequency information of each section aiming at the reservoir section;
2) the information of the logging information is finally obtained by extracting the horizon, the reservoir thickness and the lithology characteristics of the well information and solving the root mean square speed of the stratum.
6. A method of reservoir constrained deconvolution according to claim 4, characterized by: the specific method for the control analysis of the horizon in the second step comprises the following steps:
1) seismic horizon picking and description of reservoir boundary and thickness are carried out on the seismic data which are not subjected to wavelet processing in the step one, and strata which are not easy to track or lack are calibrated;
2) projecting the horizon and reservoir information obtained in the step 1) to data processed by different wavelets, and obtaining the most complete reservoir data through comparison and analysis of the horizon.
7. A method of reservoir constrained deconvolution according to claim 1, characterized by: and determining the standard wavelet through the frequency bandwidth, the main frequency, the energy and the phase data in the fourth step.
8. A method of reservoir constrained deconvolution according to claim 7, characterized by: the method for determining the standard wavelet in the fourth step comprises the following steps:
1) calibrating according to the information of the stratum characteristics contained in the logging information and the result of the wavelet in the second step, and determining the wavelet matched with the logging information;
2) selecting wavelets meeting the stability of both continuity and traceability in a reservoir in the step 1);
3) and finally determining the wavelet which can reflect the reservoir characteristics most in the step 2) by combining the reservoir characteristics of the frequency-enhanced data according to the logging information of the reservoir.
9. Use of standard wavelets obtained using the method of reservoir constrained deconvolution of any of claims 1-8 to achieve reservoir constrained deconvolution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911397428.3A CN111060970B (en) | 2019-12-30 | 2019-12-30 | Reservoir constrained deconvolution method and application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911397428.3A CN111060970B (en) | 2019-12-30 | 2019-12-30 | Reservoir constrained deconvolution method and application |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111060970A CN111060970A (en) | 2020-04-24 |
| CN111060970B true CN111060970B (en) | 2022-06-03 |
Family
ID=70305169
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911397428.3A Active CN111060970B (en) | 2019-12-30 | 2019-12-30 | Reservoir constrained deconvolution method and application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111060970B (en) |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2933499B1 (en) * | 2008-07-03 | 2010-08-20 | Inst Francais Du Petrole | METHOD OF JOINT INVERSION OF SEISMIC DATA REPRESENTED ON DIFFERENT TIME SCALES |
| FR2957682B1 (en) * | 2010-03-16 | 2012-03-23 | Inst Francais Du Petrole | METHOD FOR MONITORING A GEOLOGICAL STORAGE SITE OF GAS BY STRATIGRAPHIC REVERSAL OF SEISMIC DATA |
| CN102998703B (en) * | 2012-11-23 | 2015-07-08 | 中国石油大学(北京) | Method and device for conducting reservoir prediction and based on earth surface consistency deconvolution |
| CN104656142B (en) * | 2013-11-19 | 2017-06-06 | 中国石油天然气集团公司 | One kind is using vertical seismic profiling (VSP) and the united seismic layer labeling method of well logging |
| CN103954992B (en) * | 2014-04-03 | 2017-04-05 | 中国石油天然气股份有限公司 | A kind of the Method of Deconvolution and device |
| CN104181589A (en) * | 2014-08-20 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Nonlinear deconvolution method |
| BR112017004124A2 (en) * | 2014-09-10 | 2017-12-05 | Halliburton Energy Services Inc | integrity evaluation method of joints adhering wellbore casing to a wellbore, and wellbore system. |
| CN105044772B (en) * | 2015-08-06 | 2017-09-26 | 中国石油天然气集团公司 | A kind of method and apparatus for realizing time-varying spectrum analog deconvolution |
| CN105182408B (en) * | 2015-08-28 | 2017-09-01 | 中国石油天然气集团公司 | The preparation method and device of a kind of synthetic seismogram |
| CN106405654A (en) * | 2016-10-26 | 2017-02-15 | 成都理工大学 | Seismic spectrum imaging method based on deconvolution generalized S transform |
| US11360233B2 (en) * | 2017-09-12 | 2022-06-14 | Schlumberger Technology Corporation | Seismic image data interpretation system |
| CN109100784B (en) * | 2018-06-08 | 2020-04-28 | 恒泰艾普（北京）能源科技研究院有限公司 | Three-dimensional VSP source detection interchange full-wave-field imaging method |
| CN109283581B (en) * | 2018-11-20 | 2020-08-07 | 成都理工大学 | Reservoir gas content evaluation method based on depth domain seismic wave frequency dispersion analysis |
| CN109736791A (en) * | 2019-01-02 | 2019-05-10 | 中国石油大港油田勘探开发研究院 | Second interpretation technical method towards reservoir in a kind of oil and gas development |
-
2019
- 2019-12-30 CN CN201911397428.3A patent/CN111060970B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111060970A (en) | 2020-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106597532B (en) | Pre-stack seismic data frequency band expanding method combining well data and horizon data | |
| CN109425896B (en) | Dolomite oil and gas reservoir distribution prediction method and device | |
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| CN103837893B (en) | Acoustic logging bearing calibration based on VSP data | |
| WO2018107904A1 (en) | Method for precisely inverting young's modulus and poisson's ratio | |
| WO2017024702A1 (en) | Inversion system for ray elastic parameter | |
| US11230922B2 (en) | Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations | |
| CN107784159B (en) | Method for determining anisotropic coefficient of reservoir resistivity | |
| CN105700013B (en) | The method for building shale gas dessert productivity model | |
| CN106443774A (en) | Method for improving pre-stack depth migration imaging precision of irregular earth surface | |
| CN108427140A (en) | A method of being used for small scale fracture and cave reservoir seismic recognition | |
| CN111060970B (en) | Reservoir constrained deconvolution method and application | |
| CN108254785B (en) | Reservoir determination method and device | |
| CN109657258B (en) | Method and device for processing mechanical characteristic profile of stratum rock | |
| CN105277981A (en) | Non-consistent time-lapse seismic bin matching method based on wave field continuation compensation | |
| CN109884701B (en) | Geologic body scattering angle guiding depth imaging method | |
| CN114966828A (en) | Seismic reflection characteristic enhancing method for fractured solution reservoir | |
| CN112147700A (en) | Low-frequency model construction method and system for speed abnormal area | |
| CN113960658B (en) | Logging constraint speed modeling method and device based on geological seismic model | |
| CN107561580B (en) | The method for establishing initial geological model based on the more wells of multilayer | |
| CN112433247B (en) | While-drilling adjusting method and device for position of stratum to be drilled | |
| CN113945978B (en) | Depth domain seismic data optimization method and device | |
| CN110824565B (en) | Construction method and system of space-variant wavelet body | |
| CN111827966B (en) | Multi-well acoustic logging curve consistency processing method and device and storage medium | |
| CN111060966B (en) | First arrival picking method for static correction technology |
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 |