CN102176055B - Amplitude preserving treatment analysis and evaluation method - Google Patents
Amplitude preserving treatment analysis and evaluation method Download PDFInfo
- Publication number
- CN102176055B CN102176055B CN201110040379A CN201110040379A CN102176055B CN 102176055 B CN102176055 B CN 102176055B CN 201110040379 A CN201110040379 A CN 201110040379A CN 201110040379 A CN201110040379 A CN 201110040379A CN 102176055 B CN102176055 B CN 102176055B
- Authority
- CN
- China
- Prior art keywords
- segy
- geological data
- width
- cloth
- original
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the seismic data treatment field, in particular relating to a method for analyzing and evaluating the amplitude preservation performance of seismic treatment process. The technical scheme of the invention is as follows: inputting initial SEGY seismic data and SEGY seismic data subjected to amplitude preserving treatment; analyzing the seismic record difference; analyzing horizontal amplitude energy; analyzing horizontal frequency; analyzing the horizontal signal-to-noise ratio; and analyzing auto-correlation of seismic wavelets. The method provided by the invention has the following effect: amplitude and fidelity preservation is the key factor for the seismic data treatment quality, therefore, a manure method is needed to analyze and evaluate the amplitude preservation performance of the seismic data treatment process. In the method provided by the invention, the analysis of the amplitude preservation performance of the seismic records is mainly focused, and the comprehensive quality control can be carried out on the entire treatment process, especially on the intermediate process and the final results. The method has the advantages that the quality of the treatment process can be actually reflected through a quantized chart so as to obviate the direct loss to the effective component of the records, and processors can directly determine the amplitude preservation performance of the treatment process.
Description
Technical field
The invention belongs to the seismic data process field, is the method that a kind of assay earthquake treatment scheme is protected width of cloth property.
Prior art
Just at present, own to the treatment technology of structure imaging through comparative maturity.But under the overall background that the oil field prospecting degree is constantly deepened, seek the main flow that lithologic deposit and subtle reservoir become oil-gas exploration, the identification of lithologic trap is increasingly high to guarantor's width of cloth processing requirements of seismic data.For for the structure and the high-resolution seismic survey of lithologic deposit a little; How to monitor whole flow process and each parameter that seismic data is handled effectively, for the explanation personnel provide the seismic section of the underground form of true reflection, outstanding hydrocarbon-bearing pool characteristic most important.Monitoring means in the past seismic data work of treatment are less; It generally is the working experience that relies on treatment people; Rely on the correctness of the trend of seismic horizon in the visual determination section, be difficult to objective, actual treatment scheme and the processing parameter accurately of providing rationally, be difficult to judge guarantor's width of cloth property of each treatment step more; Directly influence the credibility of seismic section, caused explanation personnel's puzzlement.
Summary of the invention
The present invention protects unified evaluation criterion that the width of cloth handles and the disappearance problem of evaluation method proposes now in order to solve.Its purpose is in order to accomplish the accurate explanation of structure and lithologic body a little, to guarantee the levels of precision and the accuracy of imaging section.
The technical scheme of invention is:
1, the SEGY geological data after the original SEGY geological data of input and guarantor's width of cloth are handled.
2, seismologic record variance analysis
Handle back SEGY geological data with original SEGY geological data and guarantor's width of cloth and subtract each other, obtain the difference seismologic record, analyze the composition that whether comprises noise in effective signal content and the difference seismologic record in the difference seismologic record, guarantor's width of cloth property of judgment processing data.
3, along the analysis of layer amplitude energy
Select certain one deck position lineups; Analyze the amplitude energy curve along lineups; Analyze original SEGY geological data and handle the situation of change of the amplitude energy curve of back SEGY geological data, judge and protect the degree of agreement that the width of cloth is handled back SEGY geological data amplitude energy curve and original SEGY geological data with protecting the width of cloth.
4, along layer frequency analysis
Select certain one deck position lineups; Analyze the geological data frequecy characteristic along lineups; Analyze original SEGY geological data and the change of frequency situation of protecting width of cloth processing back SEGY geological data; Provide and protect the quantitative statistics information that the width of cloth is handled back SEGY geological data frequency and original SEGY geological data frequency, guarantor's width of cloth property of determination frequency disposal route.
5, along layer Analysis SNR
Select certain one deck position lineups; Analyze geological data signal to noise ratio (S/N ratio) characteristic along lineups; Analyze original SEGY geological data and the signal to noise ratio (S/N ratio) situation of change of protecting width of cloth processing back SEGY geological data; The quantitative statistics information that the width of cloth is handled back SEGY geological data and original SEGY geological data signal to noise ratio (S/N ratio) of protecting is provided, judges guarantor's width of cloth property of signal to noise ratio (S/N ratio) disposal route.
6, seismologic record Wavelet Auto-Correlated
Whether eliminated the influence of near surface in order to monitor processing links, the method intuitively of needing is analyzed the seismic event characteristic of the whole district, and statistics autocorrelation analysis technology is exactly a practicality and effective analytical approach.Because the wavelet of seismologic record is unknown; Therefore can utilize the auto-correlation of seismologic record to replace the auto-correlation (differing a scale-up factor between the two) of seismic wavelet, thereby indirect analysis different disposal method, different disposal parameter are to the influence of seismic wavelet characteristic.
The effect of invention: earthquake money is handled the quality of achievement, and protecting width of cloth fidelity is key issue.Therefore when EVALUATION DATA is handled quality, need the proven technique method guarantor's width of cloth property of processing procedure is carried out assay.The present invention lays particular emphasis on the analysis of protecting width of cloth property from seismologic record itself; This method can be carried out comprehensive quality control to the overall process of handling; Especially to middle process and net result; Come the true quality that reflects processing procedure through the demonstration that quantizes chart, avoid the effective constituent of record is caused direct loss, make treatment people guarantor's width of cloth property of determinating treatment course intuitively.
Description of drawings
Fig. 1 is a principle of work process flow diagram of the present invention.
Fig. 2 seismic section is protected width of cloth Treatment Analysis evaluation effect figure.
Fig. 3 is single big gun amplitude energy curve comparative analysis figure.
Fig. 4 signal-noise ratio of seismic sections two dimensional analysis figure.
Fig. 5 seismic section frequency plane analysis chart.
Embodiment
Protecting width of cloth Treatment Analysis evaluation method to the utilization of Shengli Oil Field typical earthquake data is:
1, the SEGY geological data after the original SEGY geological data of input and guarantor's width of cloth are handled.
2, seismologic record variance analysis
Handle back SEGY geological data with original SEGY geological data and guarantor's width of cloth and subtract each other, obtain the difference seismologic record, analyze the composition that whether comprises noise in effective signal content and the difference seismologic record in the difference seismologic record, guarantor's width of cloth property of judgment processing data.
3, along the analysis of layer amplitude energy
Select certain one deck position lineups; Analyze the amplitude energy curve along lineups; Analyze original SEGY geological data and handle the situation of change of the amplitude energy curve of back SEGY geological data, judge and protect the degree of agreement that the width of cloth is handled back SEGY geological data amplitude energy curve and original SEGY geological data with protecting the width of cloth.
4, along layer frequency analysis
Select certain one deck position lineups; Analyze the geological data frequecy characteristic along lineups; Analyze original SEGY geological data and the change of frequency situation of protecting width of cloth processing back SEGY geological data; Provide and protect the quantitative statistics information that the width of cloth is handled back SEGY geological data frequency and original SEGY geological data frequency, guarantor's width of cloth property of determination frequency disposal route.
5, along layer Analysis SNR
Select certain one deck position lineups; Analyze geological data signal to noise ratio (S/N ratio) characteristic along lineups; Analyze original SEGY geological data and the signal to noise ratio (S/N ratio) situation of change of protecting width of cloth processing back SEGY geological data; The quantitative statistics information that the width of cloth is handled back SEGY geological data and original SEGY geological data signal to noise ratio (S/N ratio) of protecting is provided, judges guarantor's width of cloth property of signal to noise ratio (S/N ratio) disposal route.
6, seismologic record Wavelet Auto-Correlated
Whether eliminated the influence of near surface in order to monitor processing links, the method intuitively of needing is analyzed the seismic event characteristic of the whole district, and statistics autocorrelation analysis technology is exactly a practicality and effective analytical approach.Because the wavelet of seismologic record is unknown; Therefore can utilize the auto-correlation of seismologic record to replace the auto-correlation (differing a scale-up factor between the two) of seismic wavelet, thereby indirect analysis different disposal method, different disposal parameter are to the influence of seismic wavelet characteristic.
According to the method described above; Finally obtain utilizing and protect width of cloth Treatment Analysis evaluation method to the seismic section (Fig. 1 Lower Half) before and after the treatment scheme; Carry out the detail analysis evaluation, obtained the curve map (Fig. 1 first half) of quantitative test, obtained good assay effect; Single big gun amplitude energy curve comparative analysis figure (Fig. 2), signal-noise ratio of seismic sections two dimensional analysis figure (Fig. 3) and seismic section frequency plane analysis chart a plurality of assay results such as (Fig. 4) have been obtained.
Claims (1)
1. protect width of cloth Treatment Analysis evaluation method for one kind, it is characterized in that:
1) original SEGY geological data of input and the SEGY geological data of protecting after the width of cloth is handled;
2) seismologic record variance analysis; Handling back SEGY geological data with original SEGY geological data and guarantor's width of cloth subtracts each other; Obtain the difference seismologic record, analyze the composition that whether comprises noise in effective signal content and the difference seismologic record in the difference seismologic record, guarantor's width of cloth property of judgment processing data;
3) along the analysis of layer amplitude energy; Select certain one deck position lineups; Analyze the amplitude energy curve along lineups; Analyze original SEGY geological data and handle the situation of change of the amplitude energy curve of back SEGY geological data, judge and protect the degree of agreement that the width of cloth is handled back SEGY geological data amplitude energy curve and original SEGY geological data with protecting the width of cloth;
4) along layer frequency analysis; Select certain one deck position lineups; Analyze the geological data frequecy characteristic along lineups; Analyze original SEGY geological data and handle the change of frequency situation of back SEGY geological data, provide and protect the quantitative statistics information that the width of cloth is handled back SEGY geological data frequency and original SEGY geological data frequency, guarantor's width of cloth property of determination frequency disposal route with protecting the width of cloth;
5) along layer Analysis SNR; Select certain one deck position lineups; Analyze geological data signal to noise ratio (S/N ratio) characteristic along lineups; Analyze original SEGY geological data and handle the signal to noise ratio (S/N ratio) situation of change of back SEGY geological data, the quantitative statistics information that the width of cloth is handled back SEGY geological data and original SEGY geological data signal to noise ratio (S/N ratio) of protecting is provided, judge guarantor's width of cloth property of signal to noise ratio (S/N ratio) disposal route with protecting the width of cloth;
6) the seismologic record Wavelet Auto-Correlated is analyzed, and utilizes the auto-correlation of seismologic record to replace the influence to the seismic wavelet characteristic of the auto-correlation indirect analysis different disposal method, different disposal parameter of seismic wavelet.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110040379A CN102176055B (en) | 2011-02-18 | 2011-02-18 | Amplitude preserving treatment analysis and evaluation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110040379A CN102176055B (en) | 2011-02-18 | 2011-02-18 | Amplitude preserving treatment analysis and evaluation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102176055A CN102176055A (en) | 2011-09-07 |
CN102176055B true CN102176055B (en) | 2012-10-17 |
Family
ID=44519257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110040379A Expired - Fee Related CN102176055B (en) | 2011-02-18 | 2011-02-18 | Amplitude preserving treatment analysis and evaluation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102176055B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467447A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Phased tendency energy matching seismic amplitude preservation evaluation method |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675915B (en) * | 2013-12-20 | 2016-07-06 | 中国石油天然气股份有限公司 | Method and device for estimating stratum transverse relative quality factor based on seismic data |
CN106547020B (en) * | 2015-09-17 | 2018-10-02 | 中国石油化工股份有限公司 | A kind of relative amplitude preserved processing method of seismic data |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634714A (en) * | 2009-08-26 | 2010-01-27 | 中国石油大学(华东) | Q structure establishment technology based on prestack seismic channel set earth layer |
CN101957455A (en) * | 2010-09-20 | 2011-01-26 | 中国海洋石油总公司 | Method of three-dimensional preserved-amplitude pre-stack time migration |
-
2011
- 2011-02-18 CN CN201110040379A patent/CN102176055B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101634714A (en) * | 2009-08-26 | 2010-01-27 | 中国石油大学(华东) | Q structure establishment technology based on prestack seismic channel set earth layer |
CN101957455A (en) * | 2010-09-20 | 2011-01-26 | 中国海洋石油总公司 | Method of three-dimensional preserved-amplitude pre-stack time migration |
Non-Patent Citations (1)
Title |
---|
郭树祥.地震资料保幅处理的讨论.《油气地球物理》.2009,第7卷(第1期),第1-3页. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467447A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Phased tendency energy matching seismic amplitude preservation evaluation method |
CN105467447B (en) * | 2014-09-04 | 2018-11-20 | 中国石油化工股份有限公司 | Width evaluation method is protected in the earthquake of phased trend energy match |
Also Published As
Publication number | Publication date |
---|---|
CN102176055A (en) | 2011-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kaklamanos et al. | Challenges in predicting seismic site response with 1D analyses: Conclusions from 114 KiK‐net vertical seismometer arrays | |
Ghofrani et al. | Implications of the 2011 M9. 0 Tohoku Japan earthquake for the treatment of site effects in large earthquakes | |
Liu et al. | A novel hydrocarbon detection approach via high-resolution frequency-dependent AVO inversion based on variational mode decomposition | |
Ge | Efficient mine microseismic monitoring | |
Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
Schorlemmer et al. | Probability of detecting an earthquake | |
CN101470211B (en) | Identification of q factor using s coda generated by microseism event | |
Régnier et al. | Influence of the VS profiles beyond 30 m depth on linear site effects: Assessment from the KiK‐net data | |
Lin et al. | Multi-station analysis of surface wave dispersion | |
US20100302905A1 (en) | Method for monitoring a subsoil zone, particularly during simulated fracturing operations | |
US8441891B2 (en) | Seismic data acquisition and processing quality control | |
Zhu et al. | Accurate identification of microseismic P-and S-phase arrivals using the multi-step AIC algorithm | |
Ji et al. | Effects of nonlinear soil behavior on Kappa (κ): Observations from the KiK‐Net database | |
Link et al. | An automatized XKS-splitting procedure for large data sets: Extension package for SplitRacer and application to the USArray | |
CN102176055B (en) | Amplitude preserving treatment analysis and evaluation method | |
Alan et al. | Precision of VS30 values derived from noninvasive surface wave methods at 31 sites in California | |
AU2013201862B2 (en) | Method for time-lapse wave separation | |
Taipodia et al. | Quantification of the resolution of dispersion image in active MASW survey and automated extraction of dispersion curve | |
Xue et al. | Detection of gas and water using HHT by analyzing P-and S-wave attenuation in tight sandstone gas reservoirs | |
Hallo et al. | Stochastic model to characterize high‐frequency ground motion at depth validated by KiK‐net vertical array data | |
Zheng et al. | Microseismic event denoising via adaptive directional vector median filters | |
Wang et al. | Falling weight deflectometer dispersion curve method for pavement modulus calculation | |
CN104570115B (en) | A kind of surface wave attenuation method and device | |
Cheng et al. | Q-estimation using seismic interferometry from vertical well data | |
Lv | Noise suppression of microseismic data based on a fast singular value decomposition algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121017 Termination date: 20200218 |