CN109212600A - Seismic data evaluation method and system based on resolution capability - Google Patents
Seismic data evaluation method and system based on resolution capability Download PDFInfo
- Publication number
- CN109212600A CN109212600A CN201710548357.7A CN201710548357A CN109212600A CN 109212600 A CN109212600 A CN 109212600A CN 201710548357 A CN201710548357 A CN 201710548357A CN 109212600 A CN109212600 A CN 109212600A
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- lithology
- resolution capability
- thin interbed
- 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.)
- Pending
Links
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/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Abstract
The present invention provides a kind of seismic data evaluation method and system based on resolution capability, this method comprises: collecting lithologic log data and corresponding speed data, and establishes lithology thin interbed model according to lithologic log data and speed data；It is reflection coefficient sequence by lithology thin interbed model conversation；When purpose layer region is opened, window analyzes frequency and extracts seismic wavelet；The seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic reflection record of the seismic wavelet；Obtained seismic reflection is recorded and is compared with lithology thin interbed model；By obtaining the evaluation result of the seismic data resolution capability to the two comparative analysis.Seismic data evaluation method of the invention is with strong points, and this method is the supplement of common seismic data evaluation, more efficient can intuitively evaluate seismic data.
Description
Technical field
The present invention relates to oil gas field of geophysical exploration, comment more particularly, to the seismic data based on resolution capability
Valence method and system.
Background technique
The oil exploration in east old area at present develops, large-scale self-contained oil-gas reservoir towards increasingly Subtle reservoir exploration
All advantageously, it has been found that and control the small layer of sand of remaining oil distribution, small fault block, due to being limited by seismic data resolution, bore
Well success rate is not high always, referred to as restricts the important obstruction of remaining exploration activity in next step.
It is the key points and difficulties studied at present for further increasing craven fault, the imaging capability of thin sand body and resolution ratio.
At present for the resolution ratio and acquisition quality of raising earthquake-capturing data, generally requires and carry out applying for stringent high quality in field
Work selects the wave detector that can excite the excitation condition of wideband and selection that can receive wideband, but to earthquake after having acquired
The evaluation of data is very crucial factor, because the various measures taken in collection process are provided to improve earthquake money
The frequency bandwidth and signal-to-noise ratio of material, but how effect needs a reasonable method and is judged, not so can not know
Whether the technical measures that field is taken are effective.
It is all the frequency bandwidth and dominant frequency of analysis of data in the past in evaluation, but these are the concept of a frequency,
It is difficult to improve how many on earth to one resolution capability of people, how thin of reservoir can be differentiated on earth, therefore previous method is to people
Feeling it is not intuitive.
Summary of the invention
Present invention aim to address cannot intuitive assay seismic data point in current seismic data assessment technique
The problem of distinguishing ability improves the reasonability and science of Mountain Seismic data evaluation, a correct conclusion is provided, so as to guide life
It produces.
The seismic data evaluation method based on resolution capability that the present invention provides a kind of, solves the above technical problem, leads to
It crosses and establishes thin interbed model and carry out convolution with the seismic wavelet of extraction, more can intuitively analyze the evaluation energy of seismic data
Power.
According to an aspect of the present invention, a kind of seismic data evaluation method based on resolution capability, this method packet are provided
It includes:
Lithologic log data and corresponding speed data are collected, and establishes lithology according to lithologic log data and speed data
Thin interbed model；
It is reflection coefficient sequence by lithology thin interbed model conversation；
When purpose layer region is opened, window analyzes frequency and extracts seismic wavelet；
The seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic reflection of the seismic wavelet
Record；
Obtained seismic reflection is recorded and is compared with lithology thin interbed model；
By obtaining the evaluation result of the seismic data resolution capability to the two comparative analysis.
It further, is reflection coefficient sequence by lithology thin interbed model conversation using formula I=(V1-V2)/(V1+V2)
Column,
In formula, I is reflection coefficient, and V1 is upper interval velocity, and V2 is lower interval velocity.
Further, lithology thin interbed model is established according to lithologic log data and speed data.
Further, obtained seismic reflection is recorded and is compared with lithology thin interbed model.
Further, by determining that the seismic wavelet of certain frequency range can differentiate the thin interbed sand of much thickness
Body obtains the evaluation result of the seismic data resolution capability.For example, being obtained in earthquake record according to simulation, between two wave crests
Apart from size, be determined to differentiate the protective layer sand body of much thickness.
According to another aspect of the present invention, a kind of seismic data evaluation system based on resolution capability, the system packet are provided
It includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Lithologic log data and corresponding speed data are collected, and establishes lithology according to lithologic log data and speed data
Thin interbed model；
It is reflection coefficient sequence by lithology thin interbed model conversation；
When purpose layer region is opened, window analyzes frequency and extracts seismic wavelet；
The seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic reflection of the seismic wavelet
Record；
Obtained seismic reflection is recorded and is compared with lithology thin interbed model；
By obtaining the evaluation result of the seismic data resolution capability to the two comparative analysis.
The present invention is to overcome previous data evaluation method not to an important supplement of previous earthquake data evaluation technology
It can be to the defect of the concept of one resolution capability of people.By establishing thin interbed model and carrying out convolution with the seismic wavelet of extraction,
The evaluating ability of seismic data more can be intuitively analyzed,
Detailed description of the invention
Disclosure illustrative embodiments are described in more detail in conjunction with the accompanying drawings, the disclosure above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in disclosure illustrative embodiments, identical reference label
Typically represent same parts.
Fig. 1 shows the techniqueflow chart of the method for the present invention.
Fig. 2 shows lithology model and different frequency convolution results, and wherein 30Hz is the convolution result of real data wavelet.
Fig. 3 is the earthquake record and lithology thin interbed that the seismic wavelet (30Hz in Fig. 2) for the real data extracted obtains
The comparing result of model
Specific embodiment
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that this disclosure will be more thorough and complete, and can be by the disclosure
Range is completely communicated to those skilled in the art.
The present invention provides a kind of seismic data evaluation method based on resolution capability.This method is for current seismic data
Evaluation is without intuitive effective aiming at the problem that resolution capability is evaluated, the new evaluation method of one kind of proposition.It is directed at present
The evaluation of seismic data is mainly spectrum analysis, it is therefore an objective to and dominant frequency, which to be evaluated, than old data improves how many, frequency band has widened how many,
Only give everybody the concept of a frequency size on frequency band, but these data can not intuitively illustrate point of the seismic data
It distinguishes that ability, this method are the supplements of common seismic data evaluation, more efficient intuitively seismic data can be evaluated.
Method of the invention can be achieved by the following technical measures:
Step 1, to the seismic data region of needs assessment, first collection lithologic log data and corresponding speed number
According to, and lithology thin interbed model is established according to lithologic log data and speed data；
It step 2, is reflection coefficient sequence by lithology thin interbed model conversation；
Step 3, the seismic data of needs assessment is subjected to spectrum analysis, window is analyzed frequency and mentioned when purpose layer region is opened
Take seismic wavelet；
Step 4, the seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the earthquake of the seismic wavelet
Wave reflection record；
Step 5, obtained seismic reflection is recorded and is compared with lithology thin interbed model；
Step 6, by determining that the seismic wavelet of the frequency range can differentiate much thickness to the two comparative analysis
Thin interbed sand body obtains the evaluation result of the seismic data resolution capability.
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
To enable objects, features and advantages of the present invention to be clearer and more comprehensible, hereafter cooperate institute's accompanying drawings, elaborate as
Under:
Fig. 1 is a kind of flow chart of seismic data evaluation method based on resolution capability of the invention.
It is the lithology thin interbed model established by lithologic log data shown in Fig. 2, and extracts real data wavelet
(30Hz) carries out convolution result and compares with realistic model, analyzes the resolution capability of real data, the convolution of 15Hz and 40Hz
As a result for being compared with 30Hz convolution result.
Fig. 3 is the earthquake record and lithology thin interbed that the seismic wavelet (30Hz in Fig. 2) for the real data extracted obtains
The comparing result of model, arrow both is designated as corresponding preferable position in figure, it can be seen that corresponding wave crest in waveform is thin
Certain big jacket layer positions in alternating layers, and relatively thin layer position is not shown in earthquake record, from two neighboring wave crest most
Small distance analysis, the seismic wavelet minimum can differentiate the layer of sand of 12m thickness.
As shown in Figure 1, in step 101, to the seismic data region of needs assessment, first collection lithologic log data
And corresponding speed data, and lithology thin interbed model is established according to lithologic log data and speed data.Such as shown in Fig. 2,
Middle Fig. 2A is the lithologic log data for certain well being collected into, and includes that detailed depth information is corresponding in the lithologic log data
Lithology and lithology thickness.Geological model is established in Tesseral software according to the lithologic log data, and defeated in the model
Enter the corresponding velocity information in each rock stratum, as shown in Figure 2 B.Process enters step 102.
It is reflection coefficient using formula I=(V1-V2)/(V1+V2), I, V1 is upper interval velocity, under V2 is in step 102
Lithology thin interbed model conversation is reflection coefficient sequence by interval velocity.For example, after having established model, the sgy number of output model
According to, and be entered into MATLAB, reflection series model is converted to using above-mentioned formula.Enter step 103.
In step 103, the seismic data of needs assessment is subjected to spectrum analysis, window analyzes frequency when purpose layer region is opened
And extract seismic wavelet.Process enters step 104.
In step 104, the seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic wavelet
Seismic reflection record.In the processing system, such as PROMAX, extract the wavelet of single-shot interval of interest, and the wavelet extracted with
The thin interbed reflectivity model established carries out convolution, completes convolution operation in the processing system.Will reflection series model with
The seismic wavelet of extraction carries out convolution and obtains seismic reflection record, as shown in Figure 2 C.Process enters step 105.
In step 105, obtained seismic reflection is recorded and is compared with lithology thin interbed model.Such as Fig. 2 B and figure
2C, Fig. 2 B are lithology thin interbed model, and Fig. 2 C is the seismic reflection record that different frequency obtains, with the earthquake note of 40Hz simulation
For record, the wave crest of seismic wave and the layer position of lithology thin interbed model are compared, if it is possible to can be right with regard to explanation on
It is differentiated this layer of position.From the point of view of the faults comparison of 15Hz, 30Hz, 40H earthquake record simulated and lithology thin interbed model, frequency
Rate is higher, corresponding better, higher to the layer bit resolution of lithology thin interbed model.Process enters step 106.
In step 106, by determining that the seismic wavelet of the frequency range can differentiate much thickness to the two comparative analysis
The thin interbed sand body of degree, obtains the evaluation result of the seismic data resolution capability.For example, obtained in earthquake record according to simulation,
The distance between two wave crests size is determined to differentiate the protective layer sand body of much thickness.Process terminates.
The presently disclosed embodiments is described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or technological improvement to the technology in market for best explaining each embodiment, or lead this technology
Other those of ordinary skill in domain can understand each embodiment disclosed herein.
Claims (10)
1. a kind of seismic data evaluation method based on resolution capability, which is characterized in that this method comprises:
Lithologic log data and corresponding speed data are collected, and it is thin mutually according to lithologic log data and speed data to establish lithology
Layer model；
It is reflection coefficient sequence by lithology thin interbed model conversation；
When purpose layer region is opened, window analyzes frequency and extracts seismic wavelet；
The seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic reflection note of the seismic wavelet
Record；
Obtained seismic reflection is recorded and is compared with lithology thin interbed model；
By obtaining the evaluation result of the seismic data resolution capability to the two comparative analysis.
2. the seismic data evaluation method according to claim 1 based on resolution capability, which is characterized in that utilize formula I
Lithology thin interbed model conversation is reflection coefficient sequence by=(V1-V2)/(V1+V2),
In formula, I is reflection coefficient, and V1 is upper interval velocity, and V2 is lower interval velocity.
3. the seismic data evaluation method according to claim 1 based on resolution capability, which is characterized in that recorded according to lithology
Well data and speed data establish lithology thin interbed model.
4. the seismic data evaluation method according to claim 1 based on resolution capability, which is characterized in that the ground that will be obtained
Seismic wave reflection record is compared with lithology thin interbed model.
5. the seismic data evaluation method according to claim 1 based on resolution capability, which is characterized in that by determining
The seismic wavelet of certain frequency range can differentiate the thin interbed sand body of much thickness, obtain commenting for the seismic data resolution capability
Valence result.
6. a kind of seismic data evaluation system based on resolution capability, which is characterized in that the system includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Lithologic log data and corresponding speed data are collected, and it is thin mutually according to lithologic log data and speed data to establish lithology
Layer model；
It is reflection coefficient sequence by lithology thin interbed model conversation；
When purpose layer region is opened, window analyzes frequency and extracts seismic wavelet；
The seismic wavelet of extraction and reflection coefficient sequence are subjected to convolution operation, obtain the seismic reflection note of the seismic wavelet
Record；
Obtained seismic reflection is recorded and is compared with lithology thin interbed model；
By obtaining the evaluation result of the seismic data resolution capability to the two comparative analysis.
7. the seismic data evaluation method according to claim 6 based on resolution capability, which is characterized in that utilize formula I
Lithology thin interbed model conversation is reflection coefficient sequence by=(V1-V2)/(V1+V2),
In formula, I is reflection coefficient, and V1 is upper interval velocity, and V2 is lower interval velocity.
8. the seismic data evaluation method according to claim 6 based on resolution capability, which is characterized in that recorded according to lithology
Well data and speed data establish lithology thin interbed model.
9. the seismic data evaluation method according to claim 6 based on resolution capability, which is characterized in that the ground that will be obtained
Seismic wave reflection record is compared with lithology thin interbed model.
10. the seismic data evaluation method according to claim 6 based on resolution capability, which is characterized in that pass through determination
The seismic wavelet of certain frequency range can differentiate the thin interbed sand body of much thickness out, obtain the seismic data resolution capability
Evaluation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710548357.7A CN109212600A (en) | 2017-07-06 | 2017-07-06 | Seismic data evaluation method and system based on resolution capability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710548357.7A CN109212600A (en) | 2017-07-06 | 2017-07-06 | Seismic data evaluation method and system based on resolution capability |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109212600A true CN109212600A (en) | 2019-01-15 |
Family
ID=64992399
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710548357.7A Pending CN109212600A (en) | 2017-07-06 | 2017-07-06 | Seismic data evaluation method and system based on resolution capability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109212600A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859160A (en) * | 2019-11-12 | 2021-05-28 | 中国石油天然气股份有限公司 | Method and device for predicting thickness of thin interbed sand body |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7525873B1 (en) * | 2004-07-23 | 2009-04-28 | Bf Corporation North America Inc. | Seismic inversion of conditioned amplitude spectra |
| US20130144534A1 (en) * | 2011-12-03 | 2013-06-06 | Geokinetics Acquisition Company | Multi-Component Spectral Warping Analysis for Thin Bed Resolution |
| CN105182408A (en) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | Manufacturing method and device for synthesizing earthquake record |
-
2017
- 2017-07-06 CN CN201710548357.7A patent/CN109212600A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7525873B1 (en) * | 2004-07-23 | 2009-04-28 | Bf Corporation North America Inc. | Seismic inversion of conditioned amplitude spectra |
| US20130144534A1 (en) * | 2011-12-03 | 2013-06-06 | Geokinetics Acquisition Company | Multi-Component Spectral Warping Analysis for Thin Bed Resolution |
| CN105182408A (en) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | Manufacturing method and device for synthesizing earthquake record |
Non-Patent Citations (4)
| Title |
|---|
| HONGLIU ZENG ET AL.: "How thin is a thin bed? An alternative perspective", 《THE LEADING EDGE》 * |
| 蔡希玲等: "砂泥岩薄互层分辨率的理论分析", 《石油物探》 * |
| 邓小力等: "薄层速度分析", 《石油地球物理勘探》 * |
| 邝少荣等: "美法地震资料处理及解释概况", 《石油地球物理勘探》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859160A (en) * | 2019-11-12 | 2021-05-28 | 中国石油天然气股份有限公司 | Method and device for predicting thickness of thin interbed sand body |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104216008B (en) | Downhole fracturing microseismic event identification method | |
| CN104142519B (en) | Mud rock crack oil deposit predicting method | |
| US6952649B2 (en) | Petroleum exploration and prediction apparatus and method | |
| CN105974407B (en) | A kind of Ground Penetrating Radar subterranean layer position detecting method | |
| CN107678064B (en) | Real-time extraction method for sound wave time difference | |
| CN106154324A (en) | Down-hole micro-seismic event automatic identifying method based on multiple tracks scanning superposition | |
| CN105005080A (en) | Method for identifying stratigraphic trap pinch-out line by using amplitude ratio attribute | |
| CN104360382A (en) | Method for detecting oil and gas by aid of stacked seismic data | |
| CN106401579A (en) | While-drilling multi-parameter analysis-based oil-gas-water identification method and system | |
| CN104265277B (en) | Method for extracting stratum sound velocity by using tube wave and stratum sound wave interference principle | |
| CN105353411A (en) | Description method of reservoir sand bodies on ultra-deep coal layer | |
| CN106526678A (en) | Reflection acoustic logging wave field separation method and device | |
| Imrie | The application of pattern recognition and machine learning to determine cement channeling & bond quality from azimuthal cement bond logs | |
| CN109212600A (en) | Seismic data evaluation method and system based on resolution capability | |
| CN107870359B (en) | Micro-seismic event recognition methods and device | |
| CN105064987A (en) | Interpretation and evaluation method for water layer recognition by logging while drilling Q parameter | |
| CN104977602A (en) | Control method and apparatus for earthquake data acquisition construction | |
| US20040114464A1 (en) | Petroleum exploration and prediction apparatus and method | |
| RU2233001C2 (en) | Method for improvement of seismic signal produced by drill bit with the use of drilling parameters | |
| CN114167515B (en) | Lithologic trap effectiveness identification method | |
| CN112282739B (en) | Method for identifying rock debris scatterer in measurement while drilling diameter | |
| Oktariena et al. | Gabor Deconvolution as Preliminary Method to Reduce Pitfall in Deeper Target Seismic Data | |
| CN112180465B (en) | Well-seismic combined stratum spread determination method | |
| KR101415309B1 (en) | Determination method and apparatus of the optimum distance of upward continuation for stripping | |
| CN113835123B (en) | Seismic acquisition parameter analysis method based on geological target pre-stack migration imaging |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190115 |
|
| RJ01 | Rejection of invention patent application after publication |