CN101852863B - Method for processing seismic data by using high-precision single-channel spectrum analysis technology - Google Patents
Method for processing seismic data by using high-precision single-channel spectrum analysis technology Download PDFInfo
- Publication number
- CN101852863B CN101852863B CN2009100814386A CN200910081438A CN101852863B CN 101852863 B CN101852863 B CN 101852863B CN 2009100814386 A CN2009100814386 A CN 2009100814386A CN 200910081438 A CN200910081438 A CN 200910081438A CN 101852863 B CN101852863 B CN 101852863B
- Authority
- CN
- China
- Prior art keywords
- data
- frequency
- processing
- frequency band
- seismic
- 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 discloses a method for processing seismic data by using high-precision single-channel spectrum analysis technology in petroleum geophysical exploration data processing. The method comprises the following steps of: labeling and selecting data, partitioning an analyzed time window segment into data segments, adding linear slopes at two ends of a time window and filling zero to the channel length of the primary data, performing Fourier transformation on the data of the disturbed m time window segments to obtain an amplitude spectrum, storing the amplitude spectrum into a seismic channel, determining an Nyquist frequency, arranging the calculated spectral values into a seismic data format to obtain a data frequency distribution rule, and selecting a processing method and a processing parameter according to the determined data frequency distribution rule to obtain a high-precision stratum section and determine the position of an oil-containing stratum. The method can quickly reflect the frequency characteristics of the seismic data at different processing stages, and can clearly display the problems so that the processing measure has pertinence.
Description
Technical field
The invention belongs to the geophysical prospecting for oil technology, relate to the method for in the seismic exploration data processing to high-density the data frequency spectrum being carried out express-analysis, be specifically related to utilize the method for single-channel spectrum analysis technical finesse mass seismic data.
Background technology
As everyone knows, in the geophysical prospecting for oil technology, the geological problem that improves seismic prospecting precision, solution continental sedimentary basin complexity is a systems engineering, relates to many aspects such as acquisition method, static correction technology, noise compact technique, meticulous imaging technique.The high density spatial Sampling techniques are the new developing direction of field of seismic exploration, for improving data signal to noise ratio (S/N ratio) and resolution, raising data quality, obtaining the effect that abundanter underground reflective information has uniqueness, therefore become one of focus of current seismic prospecting.The high density spatial Sampling techniques are meant implements some source excitation, some reception or small size combined reception, small distance between receivers or the observation of facet unit in the open air, seismic wave field is fully sampled, signal and noise are carried out " wide advance wide go out ", avoid coming to harm because of receiver pattern makes reflective information in the gatherer process, as far as possible the primitiveness of maintenance sampling wave field and rich.
But by the huge data volume that high density spatial sampling obtains, how adopting fast and effectively method that the data feature is analyzed is the problem that treatment people at first faces for the field; The problem that next faces be, in processing procedure, how treatment people is placed on the feature of correct understanding data, the problem that the objective analysis data exist with main energy, in this problem, how to take suitable technological means, select reasonable parameter that mass data is implemented to handle, the quality of mass data is improved constantly, make and handle the achievement needs of matter synthetic study contentedly; Moreover the problem that faces is how to grasp the problem of the quality of analysis means, Quality Control data or data map, and this to estimate the important symbol of a disposal system performance often, also is the embodiment of geological data processing enter complex art strength.Adopt multiple analytical technology to understand, grasp the feature of mass data, use quantitative analysis method that the processing quality of each link in the mass data processing process is carried out in time, monitored effectively, the details of observed data, can in time remove a hidden danger, increase work efficiency, be the problem of needing solution in the geophysical prospecting for oil technology badly always.
As everyone knows, frequency is one of key character of geological data.In order to solve top described technical matters, in disposal system commonly used at present, mutual single track or the multiple tracks spectral analysis method of many employings carries out frequency analysis to geological data, but in existing rarely seen mass data processing technology, exists various deficiencies or defective.Be that example is illustrated (referring to Fig. 1) with the ProMAX earthquake data processing system below.Interactive Spectral Analysis, promptly mutual spectral analysis method is composed estimation and demonstration to selected seismic trace; It adopts average power spectra, F-X spectrum and the average phase spectrum of Fourier transformation calculations seismic trace, and shows with alternant way.In the graphic structure shown in Fig. 1, the demonstration that mutual spectrum is estimated divides four zones: (a) input data (T-X) district; (b) corresponding road F-X spectrum; (c) multiple tracks power spectrum; (d) multiple tracks phase spectrum.The concrete grammar step is: (1) selects the input data; (2) delimit the time period to be analyzed; (3) import the parameter of Fourier conversion in the mode of windowing; (4) output is composed alternately.
For the mass seismic data of open-air high-density acquisition, conventional spectral analysis technology (as the ProMAX earthquake data processing system) mainly exists not enough in following several respects:
(1) data volume of handling when needs is big, and the data frequency feature is in the time, when direction in space there are differences, and the speed of conventional interactive mode spectrum output reduces, and is unfavorable for quick comparative analysis;
(2) F-X spectrum since the time window selectional restriction, precision is lower;
(3) fail to extract characteristic parameter on the frequency spectrum, can't realize quantizing comparative analysis.
In addition, because the singularity that geological data is handled, be that the geological data processing is that the specific aim of carrying out on the basis of overall understanding raw data holistic distribution characteristics is handled, the difficulty of data volume and processing accuracy requirement, DATA DISTRIBUTION feature, analysis is all different with other non-mass data processing methods with analyst coverage; Again because concerning seismic exploration to high-density, the data volume of collection increases several times again, more needs the data analysis technique that adapts with it.
Therefore, for the high density seismic exploration technique, be badly in need of at present a kind of method of suitable mass seismic data, with this improve the magnanimity earthquake processing accuracy with speed, reduce noise, take this to bring to obtain high-precision stratigraphic section and definite oil-bearing formation position by the identification active zone.
Summary of the invention
Deficiency at existing geological data system, the present invention is directed to the big characteristics of pre stack data amount in the high-density acquisition, the present invention proposes a kind of method of utilizing single-channel spectrum analysis technical finesse mass seismic data, and this method can be implemented geological data spectrum estimation fast and effectively and can finish the spectrum estimation of single track data as required.
May further comprise the steps according to the method for single-channel spectrum analysis technical finesse mass seismic data of utilizing of the present invention:
(1) survey line and emplacement are arranged in the preliminary in the open air area of determining that has oily to wish, use explosive source or ground vibroseis earthquake-wave-exciting, and with seismograph the propagation of seismic wave situation are noted;
(2) data that will handle are marked and choose: x (i) is pending seismic trace, and i is the sampling point sequence number, determines window Wa-Wb when to be analyzed according to whole district sign interval, objective interval and user's certain period of time;
(3) analyze the time window section on, choose 100ms with interior disturbance, and cut apart and form m data segment, m is 3～5 natural number;
(3) according to window length (number of samples of Wb-Wa+10 * n)/SI, the time window two ends add linear slope and zero padding is long to the road of raw data, wherein n is the random number between the 0-10, SI is a sampling interval;
The data of window section are done the Fourier conversion during (4) to the m after the disturbance, adopt formula (41) to calculate the spectral amplitude of single track data, and the data after the disturbance are asked on average, obtain the spectral amplitude A (f) of x (i);
In the formula: x
j(t
i) be the input seismic trace,
t
iBe the time,
J=1,2 ... m is the road sequence number after the random perturbation,
A (f) is the spectrum results estimated,
F is a frequency.
(5) the spectral amplitude A (f) with x (i) leaves in the seismic trace, determines the Nyquist frequency by sampling interval, by the end of the maximum frequency place;
(6) the spectrum value that will calculate is arranged as the geological data form, frequecy characteristic according to this geological data format checking raw data, can identify dissimilar wave frequency features rapidly and accurately with this, determine frequency values, obtain the data frequency regularity of distribution, and then improve the speed and the precision of frequency detecting;
(7) according to established data frequency distribution rule, select mass data processing method and definite rational processing parameter and flow process targetedly, improve the magnanimity earthquake processing accuracy with speed, reduce noise, take this to bring to obtain high-precision stratigraphic section and definite oil-bearing formation position by the identification active zone.
In the above-mentioned method of utilizing single-channel spectrum analysis technical finesse mass seismic data, the spectral amplitude A (f) of x (i) is left in the seismic trace, so just can handle, so that extract characteristic parameter data.
Data are carried out disposal route is: along the integral filtering or the medium filtering of frequency direction, avoid the influence of exceptional value to characteristic parameter.
A′
k(f)＝int?egralfilter[A
k(f)] (42)
A′
k(f)＝medianfilter[A
k(f)] (43)
In the formula: A '
k(f) be filtered data,
A
k(f) be the input data.
K is a sequence number
Being extracted as of characteristic parameter: the data after the Filtering Processing are carried out the extraction of characteristic parameter in different frequency band ranges; Frequency band range can be divided into the frequency band range, dominant frequency band scope, low-frequency range, high band of interference wave frequency band range, reflection wave etc., extracts the average energy (A in a certain frequency band as required
Eng) and crest frequency (f
Max), be equipped with key word and be stored in the trace header, obtain quantization characteristic curve or curved surface.
Average energy (the A of frequency band
Eng) curve can clearly reflect the energy horizontal change in the frequency band, analysing energy changes the corresponding relation that changes with near surface structure, helps treatment people to select testing site, confirmed test scheme.Crest frequency (f
Max) be the frequency of maximal value correspondence in the selected frequency band, the dominant frequency that this parameter makes treatment people can understand data apace distributes, use this parameter in the different links of data processing and can study data frequency Changing Pattern in the processing procedure, if three-dimensional data, the space distribution situation that can show express-analysis whole district data by the plane of this parameter, space-variant parameter in determining to handle and frequency variable element also can be monitored the variation of processing procedure medium frequency.
The present invention is the processing stage of high density spatial Sampling techniques indoor the time; utilize the numerical analysis method in the computing machine that noise is separated; make up neatly and superpose; thereby reach the purpose of suppressing noise, protection useful signal to greatest extent; for geologic interpretation provides high-resolution information, help meticulous reservoir description and reservoir prediction.
Brief Description Of Drawings
Fig. 1 is the signal graphic structure of ProMAX disposal system spectrum analysis;
Fig. 2 (a) is original single shot record local repressentation;
Fig. 2 (b) is the analysis of spectrum result of the inventive method;
Fig. 2 (c) is the analysis of spectrum result's of inventive method demonstration;
Fig. 3 (a) is raw data geophone offset collection demonstration altogether;
Fig. 3 (b) is that the common geophone offset collection after denoising and the amplitude processing shows;
Fig. 3 (c) is geophone offset collection demonstration altogether after static correction and the deconvolution;
Fig. 3 (d) is a stacked section;
Fig. 4 (a) is total to geophone offset collection analysis of spectrum for raw data;
Fig. 4 (b) is the analysis of spectrum of the common geophone offset collection after denoising and the amplitude processing;
Fig. 4 (c) is the analysis of spectrum that is total to the geophone offset collection after static correction and the deconvolution;
Fig. 4 (d) is the stacked section analysis of spectrum.
Embodiment
Below in conjunction with accompanying drawing, the method for utilizing single-channel spectrum analysis technical finesse mass seismic data that the present invention is proposed is described in detail.
At first, survey line and emplacement are arranged in the preliminary in the open air area of determining that has oily to wish, use explosive source or ground vibroseis earthquake-wave-exciting, and with seismograph the propagation of seismic wave situation are noted.Just, use the routine techniques means, gather the mass seismic data that acquisition will be handled in the open air.
The data that will handle are marked and choose: x (i) is pending seismic trace, and i is the sampling point sequence number, determines window Wa-Wb when to be analyzed according to whole district sign interval, objective interval and user's certain period of time.Wherein, the data that will handle also can be the geological data of acquired original, the geological data or the end result data of processing intermediate link.
Analyze the time window section on, get 100ms with interior disturbance, form m data segment, m is 3～5 natural number.((SI is a sampling interval to the number of samples of Wb-Wa+10 * n)/SI, and n is the random number between the 0-10 according to window length; Attention: n is the random number between the 0-10, is produced automatically by software, guarantees that disturbance is in 100ms; Wb-Wa is former definite (before the disturbance) window length, and middle "-" is minus sign, and the scope of Wa-Wb window when being meant, from Wa to Wb), the time window two ends add linear slope, and zero padding is long to the road of raw data; The data of window section are done the Fourier conversion during to the m after the disturbance, adopt formula (41) to calculate the spectral amplitude of single track data, and the data after the disturbance are asked on average, obtain the spectral amplitude A (f) of x (i);
In the formula: x
j(t
i) be the input seismic trace,
t
iBe the time,
J=1,2 ... m is the road sequence number after the random perturbation,
A (f) is the spectrum results estimated,
F is a frequency.
The result is left in the seismic trace, determine the Nyquist frequency, by the end of the maximum frequency place by sampling interval.
The spectrum value of calculating is arranged as the geological data form, checks the frequecy characteristic of raw data according to this, can identify dissimilar wave frequency features rapidly and accurately, determine frequency values.If the variation of research processing procedure medium frequency can adopt this method to compare, find that frequency band is unusual, improved the speed and the precision of frequency detecting.
The spectrum value of calculating is arranged as the geological data form, frequecy characteristic according to this geological data format checking raw data, as the dominant frequency of dominant frequency, frequency span, reflection wave frequency, ground roll, the dominant frequency of linear disturbance etc., compare with conventional method, can identify dissimilar wave frequency features rapidly and accurately on this method result of calculation, determine frequency values, obtain the data frequency regularity of distribution, and then improve the speed and the precision of frequency detecting.According to established data frequency distribution rule, select mass data processing method and definite rational processing parameter and flow process targetedly, improve the magnanimity earthquake processing accuracy with speed, reduce noise, take this to bring to obtain high-precision stratigraphic section and definite oil-bearing formation position by the identification active zone.
In the above-mentioned method of utilizing single-channel spectrum analysis technical finesse mass seismic data, result of calculation leaves in the seismic trace with the form of data, so just can handle data, so that extract characteristic parameter.
Disposal route:, avoid the influence of exceptional value to characteristic parameter along the integral filtering or the medium filtering of frequency direction.
A′
k(f)＝int?egratedfilter[A
k(f)] (42)
A′
k(f)＝medianfilter[A
k(f)] (43)
In the formula: A '
k(f) be filtered data,
A
k(f) be the input data.
K is a sequence number
The extraction of characteristic parameter: the data after the Filtering Processing are carried out the extraction of characteristic parameter in different frequency band ranges.Frequency band range can be divided into: the frequency band range of interference wave frequency band range, reflection wave, dominant frequency band scope, low-frequency range, high band etc., extract the average energy (A in a certain frequency band as required
Eng) and crest frequency (f
Max), be equipped with key word and be stored in the trace header, draw out quantization characteristic curve or curved surface.
Average energy (the A of certain frequency band
Eng) curve can clearly reflect the energy horizontal change in this frequency band, analysing energy changes the corresponding relation that changes with near surface structure, helps treatment people to select testing site, confirmed test scheme.Crest frequency (f
Max) be the frequency of maximal value correspondence in the selected frequency band, the dominant frequency that this parameter makes treatment people can understand data apace distributes, use this parameter in the different links of data processing and can study data frequency Changing Pattern in the processing procedure, if three-dimensional data, the space distribution situation that can show express-analysis whole district data by the plane of this parameter, space-variant parameter in determining to handle and frequency variable element also can be monitored the variation of processing procedure medium frequency.
The spectrum result that the geological data that uses the inventive method to obtain is handled has the following advantages:
(1) spectrum precision height is stored in the seismic trace in the mode of data and (notes the difference of sampling interval), is convenient to adopt the data comparative analysis in colored demonstration or visual way and original T-X territory;
(2) characteristic parameter of Ti Quing is sane, can represent the comprehensive characteristics of the data set of analyzing, and deposit in the trace header in the mode of numerical value, can curve, the mode of curved surface or color plane shows, so that, carry out synthetic study with other characteristic parameters (as face of land elevation curve, static correction value curve, energy distribution curve etc.) contrast.
(3) analysis of spectrum speed is fast, and after the analyst determined analyst coverage, parameter, the earth's surface-consistent spectrum was estimated and statistics is finished by Automatic Program, is applicable to the frequency analysis or the large-area three-dimensional data analysis of the mass data of high-density acquisition.
Below accompanying drawing is explained in detail.
In the graphic structure shown in Fig. 1, the demonstration that the mutual spectrum in the ProMax system is estimated divides four zones: (a) input data (T-X) district, and horizontal ordinate is CDP number (common midpoint number), ordinate is reflection interval (ms); (b) corresponding road F-X spectrum, horizontal ordinate are CDP number (common midpoint number), and ordinate is frequency (Hz); (c) multiple tracks power spectrum, horizontal ordinate are frequency (Hz), and ordinate is to divide the doubly power (dB) of expression; (d) multiple tracks phase spectrum, horizontal ordinate are frequency (Hz), and ordinate is phase place (angle °).The concrete grammar step is: (1) selects the input data; (2) delimit the time period to be analyzed; (3) import the parameter of Fourier conversion in the mode of windowing; (4) output is composed alternately.
Fig. 2 (a) is the demonstration of original single shot record at 1600-2200ms, horizontal ordinate FFID is a big gun number among the figure, CHAN is the record Taoist monastic name, ordinate is two-way time (Time (ms), Fig. 2 (b) is the result of this method after handling, and ordinate is frequency (Frequency (Hz)), as can be seen on little geophone offset road (Taoist monastic name scope 151-241), the ground roll energy is strong, and frequency distribution is at 5-18Hz.And on middle far offset trace, the energy of significant wave mainly is distributed in the 7-60HZ scope, and Fig. 2 (c) is that the colour of this method result shows, can more clearly understand the frequency distribution in the raw data.
Application example: the CZ area utilizes the frequency change in the spectral analysis method monitoring processing procedure
According to the requirement of this district's Data Processing task, seismic data will carry out meticulous pre-stack time migration to be handled, and requires section to have higher signal to noise ratio (S/N ratio) and resolution, can satisfy follow-up reservoir prediction and hydro carbons and detect the needs of handling.
Because the frequency fidelity that this district's geological tasks is handled data is had relatively high expectations, adopt the feature of spectral factorization method research reservoir in the seismic interpretation, so the analysis to frequecy characteristic is very important in the data processing.The conventional only variation of data quality as can be seen of analytical approach, but be difficult to the frequency distribution situation of analysis list track data, and analyst coverage is not comprehensive.After adopting this method that the data processing procedure has been carried out the all-the-way tracking monitoring, by diagnosis, analyzed the attenuation of interference wave (as ground roll, high frequency interference etc.), study the denoising means of being taked and whether useful signal has been had infringement, find the problem of denoising and adjusted processing parameter, guaranteed the precision that dominant frequency band is handled, for the good data basis has been established in the application of the Spectral Decomposition Technique in the geologic interpretation.
Fig. 3 (a) and (b), (c) are respectively the common geophone offset record (geophone offset is 2170m) after common geophone offset record, deconvolution and the static correction after original altogether geophone offset record, noise attentuation and the amplitude processing, and Fig. 3 (d) be the stacked section of correspondence.CDP is a common midpoint number among the figure, and ordinate is reflection interval (Time (ms)), and the display mode of this routine only can be from the variation of visually seeing the data frequency composition and the improve of signal to noise ratio (S/N ratio), but can not quantize to illustrate the change amount of each recording channel medium frequency.Fig. 4 (a) and (b), (c), (d) are the result with the corresponding this method of Fig. 3, horizontal ordinate is CDP number among the figure, ordinate is frequency (Frequency (Hz)), therefrom can be in original geophone offset record (Fig. 4 (a)) altogether, exist frequency to be higher than the high frequency interference of 80Hz on indivedual roads, the low-frequency disturbance that is lower than 10Hz is arranged on some road, and the signal of Mid Frequency 10-45Hz energy distribution in the horizontal is extremely inhomogeneous; Handle back (Fig. 4 (b)) through noise attentuation and earth's surface-consistent, the energy distribution in the data is more even, but dominant frequency band can not satisfy the explanation requirement of local area to thin reservoir at 15-50Hz, also needs to do the deconvolution processing and improves resolution; Fig. 4 (c) is the frequency spectrum of the common geophone offset record after static correction and the deconvolution, and frequency band has obtained widening, and the front end signal frequency can reach 80Hz, and energy distribution is even; Superpose on this basis, the frequency spectrum of stack result is shown in Fig. 4 (d), and stack has improved the signal to noise ratio (S/N ratio) of data, can find out with phasic property from the frequency spectrum obviously to strengthen, but the radio-frequency component of signal is had loss, and the above energy of 70Hz weakens to some extent.
Adopt the banded amount fractional analysis technology among the present invention, the frequecy characteristic that can reflect different disposal stage geological data fast, each link can clearly show the change of frequency content and the problem of existence, makes the treatment measures of taking have stronger specific aim.
As above-mentioned, the clear method that has described the present invention's proposition in detail.Although the present invention is described and explained to the preferred embodiments of the present invention in detail, but this area those skilled in the art be appreciated that, under the situation of the spirit and scope of the present invention that do not deviate from the claims definition, can in form and details, make multiple modification.
Claims (5)
1. method of utilizing single-channel spectrum analysis technical finesse geological data, it may further comprise the steps:
1) survey line and emplacement are arranged in the preliminary in the open air area of determining that has oily to wish, use explosive source or ground vibroseis earthquake-wave-exciting, and with seismograph the propagation of seismic wave situation are noted;
2) data that will handle are marked and choose: x (i) is pending seismic trace, and i is the sampling point sequence number, determines window scope Wa-Wb when to be analyzed according to whole district sign interval, objective interval and user's certain period of time; Window zero-time when Wa is, window termination time when Wb is, the ms of unit.
3) analyze the time window section on, choose 100ms with interior disturbance, and cut apart and form m data segment, m is 3～5 natural number;
3) according to window length (number of samples of Wb-Wa+10 * n)/SI, the time window two ends add linear slope and zero padding is long to the road of raw data, wherein n is the random number between the 0-10, SI is a sampling interval;
The data of window section are done the Fourier conversion during 4) to the m after the disturbance, adopt formula (41) to calculate the spectral amplitude of single track data, and the data after the disturbance are asked on average, obtain the spectral amplitude A (f) of x (i);
In the formula: x
j(t
i) be the input seismic trace,
t
iBe the time,
J=1,2 ... m is the road sequence number after the random perturbation,
A (f) is the spectrum results estimated,
F is a frequency.
5) the spectral amplitude A (f) with x (i) leaves in the seismic trace, determines the Nyquist frequency by sampling interval, by the end of the maximum frequency place;
6) the spectrum value that will calculate is arranged as the geological data form, according to the frequecy characteristic of this geological data format checking raw data, identify dissimilar wave frequency features with this, determine frequency values, obtain the data frequency regularity of distribution, and then improve the speed and the precision of frequency detecting;
7) according to established data frequency distribution rule, select data processing method and definite processing parameter and flow process, improve the earthquake processing accuracy with speed, reduce noise, take this to bring to obtain stratigraphic section and definite oil-bearing formation position by the identification active zone.
2. according to the process of claim 1 wherein that the spectral amplitude A (f) with x (i) leaves in the seismic trace, so just can handle, so that extract characteristic parameter to data.
3. according to the method for claim 2, wherein data are carried out disposal route and be:, avoid the influence of exceptional value characteristic parameter along the integral filtering or the medium filtering of frequency direction.
A′
k(f)＝int?egralfilter[A
k(f)]?(42)
A′
k(f)＝medianfilter[A
k(f)] (43)
In the formula: A '
k(f) be filtered data,
A
k(f) be the input data.
K is a sequence number
4. according to the method for claim 2, being extracted as of characteristic parameter wherein: the data after the Filtering Processing are carried out the extraction of characteristic parameter in different frequency band ranges; Frequency band range can be divided into frequency band range, dominant frequency band scope, low-frequency range, the high band of interference wave frequency band range, reflection wave, extracts the average energy A in a certain frequency band as required
EngWith crest frequency f
Max, be equipped with key word and be stored in the trace header, obtain quantization characteristic curve or curved surface.
5. according to the method for claim 4, the average energy A of its midband
EngCurve can clearly reflect the energy horizontal change in the frequency band, and analysing energy changes the corresponding relation that changes with near surface structure; Crest frequency f
MaxIt is the frequency of maximal value correspondence in the selected frequency band.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100814386A CN101852863B (en) | 2009-04-03 | 2009-04-03 | Method for processing seismic data by using high-precision single-channel spectrum analysis technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100814386A CN101852863B (en) | 2009-04-03 | 2009-04-03 | Method for processing seismic data by using high-precision single-channel spectrum analysis technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101852863A CN101852863A (en) | 2010-10-06 |
| CN101852863B true CN101852863B (en) | 2011-12-07 |
Family
ID=42804435
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100814386A Expired - Fee Related CN101852863B (en) | 2009-04-03 | 2009-04-03 | Method for processing seismic data by using high-precision single-channel spectrum analysis technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101852863B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103675915A (en) * | 2013-12-20 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for estimating stratum transverse relative quality factors based on seismic data |
| CN104422960A (en) * | 2013-09-02 | 2015-03-18 | 中国石油化工股份有限公司 | Seismic data fluid identification method based on self-adaption extraction of signal low-frequency intense anomaly |
| CN104536042A (en) * | 2014-12-24 | 2015-04-22 | 中国石油天然气股份有限公司 | Two-dimension superposed seismic data amplitude compensation method and device |
| CN105510963A (en) * | 2014-09-26 | 2016-04-20 | 中国石油化工股份有限公司 | Method for analyzing seismic exploration auxiliary data files |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323617B (en) * | 2011-06-13 | 2014-03-12 | 中国石油化工股份有限公司 | Merging processing method of 2D seismic data of complex surfaces |
| CN103217713B (en) * | 2012-01-19 | 2015-08-05 | 中国石油化工股份有限公司 | Oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND data optimization methods |
| CN103323876B (en) * | 2012-03-22 | 2015-08-19 | 中国石油天然气集团公司 | A kind of method determining the best low frequency sweep signal of vibroseis |
| GB2515218B (en) * | 2012-04-03 | 2017-08-16 | Total Sa | Computerized method of characterizing a geological region of interest and computer program implementing this method |
| CN103728661B (en) * | 2012-10-16 | 2016-08-03 | 中国石油化工股份有限公司 | A kind of high accuracy inverse Q filtering Processing Seismic Data |
| CN104820243B (en) * | 2015-04-20 | 2017-04-12 | 中国石油天然气股份有限公司 | Three-dimensional post-stack seismic data amplitude compensation method and device |
| CN107180391B (en) * | 2017-03-31 | 2023-03-24 | 中国电力科学研究院 | Wind power data span selection method and device |
| CN109031418A (en) * | 2018-07-17 | 2018-12-18 | 中国石油天然气股份有限公司 | A kind of earthquake collection quality analysis method and device |
| CN109212601A (en) * | 2018-08-31 | 2019-01-15 | 中国石油化工股份有限公司 | A kind of seismic data exception measuring point detection method |
| JP7236312B2 (en) * | 2019-04-04 | 2023-03-09 | 富士フイルムヘルスケア株式会社 | Ultrasound diagnostic device, signal processing device, and program |
| CN113126154A (en) * | 2019-12-30 | 2021-07-16 | 李智 | Surface wave spectrum analysis method for earthquake prediction |
| CN111765960B (en) * | 2020-07-23 | 2022-02-25 | 国网山西省电力公司太原供电公司 | Method for extracting seismic signals of OPGW (optical fiber composite overhead ground wire) optical cable based on distributed optical fiber sensing |
| CN112099086B (en) * | 2020-09-16 | 2022-03-29 | 中油奥博(成都)科技有限公司 | High-resolution optical fiber borehole seismic data deep frequency analysis method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998025161A2 (en) * | 1996-12-06 | 1998-06-11 | Amoco Corporation | Spectral decomposition for seismic interpretation |
| WO1998037437A1 (en) * | 1997-02-20 | 1998-08-27 | Schlumberger Canada Limited | A method of processing seismic data signals |
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| WO2004095073A2 (en) * | 2003-04-01 | 2004-11-04 | Exxonmobil Upstream Research Company | Shaped high frequency vibratory source |
| CN101221249A (en) * | 2007-01-10 | 2008-07-16 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic data quality analyzing technique |
-
2009
- 2009-04-03 CN CN2009100814386A patent/CN101852863B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998025161A2 (en) * | 1996-12-06 | 1998-06-11 | Amoco Corporation | Spectral decomposition for seismic interpretation |
| WO1998037437A1 (en) * | 1997-02-20 | 1998-08-27 | Schlumberger Canada Limited | A method of processing seismic data signals |
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| WO2004095073A2 (en) * | 2003-04-01 | 2004-11-04 | Exxonmobil Upstream Research Company | Shaped high frequency vibratory source |
| CN101221249A (en) * | 2007-01-10 | 2008-07-16 | 中国石油集团东方地球物理勘探有限责任公司 | Seismic data quality analyzing technique |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104422960A (en) * | 2013-09-02 | 2015-03-18 | 中国石油化工股份有限公司 | Seismic data fluid identification method based on self-adaption extraction of signal low-frequency intense anomaly |
| CN103675915A (en) * | 2013-12-20 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for estimating stratum transverse relative quality factors based on seismic data |
| CN103675915B (en) * | 2013-12-20 | 2016-07-06 | 中国石油天然气股份有限公司 | Based on the method and apparatus that seismic data estimates the laterally opposed quality factor in stratum |
| CN105510963A (en) * | 2014-09-26 | 2016-04-20 | 中国石油化工股份有限公司 | Method for analyzing seismic exploration auxiliary data files |
| CN105510963B (en) * | 2014-09-26 | 2018-03-06 | 中国石油化工股份有限公司 | A kind of method for analyzing seismic prospecting auxiliary data file |
| CN104536042A (en) * | 2014-12-24 | 2015-04-22 | 中国石油天然气股份有限公司 | Two-dimension superposed seismic data amplitude compensation method and device |
| CN104536042B (en) * | 2014-12-24 | 2017-01-25 | 中国石油天然气股份有限公司 | Two-dimension superposed seismic data amplitude compensation method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101852863A (en) | 2010-10-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101852863B (en) | Method for processing seismic data by using high-precision single-channel spectrum analysis technology | |
| Foti et al. | Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project | |
| CN100557464C (en) | Seismic prospecting horizon calibration method based on the prestack wave-field simulation | |
| Behm et al. | Retrieval of local surface wave velocities from traffic noise–an example from the La Barge basin (Wyoming) | |
| CN102305941B (en) | Method for determining stratum stack quality factor by direct scanning of prestack time migration | |
| CN100429527C (en) | Three-dimensional seismic data processing quality monitoring technology | |
| CN103605157B (en) | The method of decay near surface scattering wave | |
| CN103454685A (en) | Method and device for predicating sand body thicknesses through logging constraint wave impedance inversion | |
| CN101852864B (en) | Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology | |
| CN103728659A (en) | Method for improving underground karst detecting precision | |
| Gao et al. | Estimation of quality factor Q from the instantaneous frequency at the envelope peak of a seismic signal | |
| CN103257363A (en) | Method for detecting inclination angle of fissure in underground fissure type reservoir stratum | |
| CN103645499B (en) | Based on the earth surface consistency vibration amplitude compensation method of poststack reflected energy statistics | |
| Dal Moro et al. | Shear-wave velocity profiling according to three alternative approaches: a comparative case study | |
| CN103592682A (en) | Seismic wave field energy balancing method and processing device | |
| Song et al. | Utilization of multimode surface wave dispersion for characterizing roadbed structure | |
| CN102253414B (en) | Reservoir detecting method based on analysis of earthquake lines | |
| Boaga et al. | Shear wave structural models of Venice Plain, Italy, from time cross correlation of seismic noise | |
| CN102230973A (en) | Three-dimensional split-step Fourier viscous sound wave depth migration method | |
| CN104199087B (en) | Method and device for inverting sea water depth by use of data of underwater detector and land detector | |
| CN102385066A (en) | Pre-stack seismic quantitative imaging method | |
| CN104597485A (en) | Micro-fault detecting method and fault detecting device | |
| Kolkman-Quinn | Time-lapse VSP monitoring of CO2 sequestration at the CaMI Field Research Station | |
| Bouchaala et al. | Azimuthal Investigation of a Fractured Carbonate Reservoir | |
| CN104597497A (en) | Reservoir hydrocarbon prediction method based on prestack instantaneous frequency attribute analysis |
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: 20111207 |