CN103217713A - Method for optimizing seismic velocity analysis data in oil and gas exploration - Google Patents
Method for optimizing seismic velocity analysis data in oil and gas exploration Download PDFInfo
- Publication number
- CN103217713A CN103217713A CN2012100171167A CN201210017116A CN103217713A CN 103217713 A CN103217713 A CN 103217713A CN 2012100171167 A CN2012100171167 A CN 2012100171167A CN 201210017116 A CN201210017116 A CN 201210017116A CN 103217713 A CN103217713 A CN 103217713A
- Authority
- CN
- China
- Prior art keywords
- data
- frequency
- oil
- gas exploration
- analysis
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for optimizing seismic velocity analysis data in oil and gas exploration. The method for optimizing the seismic velocity analysis data in the oil and gas exploration includes the following steps: velocity trace gather sorting is carried out on the seismic data, and an ultra trace gather used for velocity analysis is sorted out; sub-band abnormal amplitude attenuation is carried out; F-X-Y domain three-dimensional prestack random noise attenuation is carried out; frequency filtering is carried out, so that energy of a superiority frequency band of reflection waves can be enhanced; and wavelet shaping deconvolution processing is carried out to improve resolution. The method for optimizing the seismic velocity analysis data in the oil and gas exploration solves the problems that in the prior art, focusing of an energy group in velocity analysis is poor and precision is low, and has the advantages that effective reflection energy can be enhanced, the quality of the velocity analysis data is improved, the signal-to-noise ratio and the resolution of the velocity analysis trace gather are improved, and the focusing property of a velocity spectrum can be improved.
Description
Technical field
The present invention relates to the seismic data process field, particularly relate to the data-optimized method of a kind of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND.
Background technology
Seismic imaging requires more and more higher to the analysis precision of speed, present conventional speeds analytical approach, the main very big diagnostic method of poststack energy that adopts, the velocity analysis process is not carried out the specific aim optimization process to participating in velocity analysis ground geological data, so that make speed energy out-focus, the velocity accuracy of mid-deep strata is not enough, has influenced ground quality imaging precision.Though developed the velocity analysis method that four item velocity analysiss are smoothed collection as far as possible early stage, can't solve mid-deep strata speed energy focusing problem.At mid-deep strata velocity spectrum energy group out-focus, problem that the velocity analysis precision is low, do not see corresponding velocity analysis method.And, when producing velocity spectrum usually, all be the result who directly utilizes overall data process.When overall data process, consider factors such as amplitude, frequency in length and breadth to fidelity, when using the noise compact technique, denoising effect usually can be weaker, and is bigger to the velocity spectrum influence.We have invented a kind of new data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND for this reason, at the velocity analysis data have that signal to noise ratio (S/N ratio) is low, a little less than the usable reflection, cause the relatively poor problem of velocity spectrum energy group focusing, to participating in the data of velocity analysis, to improve road collection signal to noise ratio (S/N ratio) and velocity spectrum quality is main target, adopt disposal route and the parameter different, can access better velocity spectrum with overall data.
Summary of the invention
The purpose of this invention is to provide one and can plant the signal to noise ratio (S/N ratio) of raising velocity analysis road collection and the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND of resolution.
Purpose of the present invention can realize by following technical measures: the data-optimized method of this oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND comprises: step 1, geological data is carried out the collection sorting of speed road, and sub-elect the channel set that is used for velocity analysis; Step 2, the geological data that step 1 is generated carries out the decay of branch frequency band anomalous amplitude; Step 3, the geological data that step 2 is generated carries out the three-dimensional prestack random noise decay in F-X-Y territory; Step 4, the geological data that step 3 is generated carries out frequency filtering, and the dominant frequency band energy of reflection wave is strengthened; And step 5, the geological data that step 4 is generated carries out wavelet shaping deconvolution processing, improves resolution.
Purpose of the present invention also can realize by following technical measures:
The data-optimized method of this oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND also is included in before the step 1, and input is through processed conventionally this geological data.
In step 2, the geological data that step 1 is generated carries out frequency division band scanning analysis earlier, determine the noise band scope after, carry out the decay of branch frequency band anomalous amplitude again.
In step 2, when carrying out the decay of branch frequency band anomalous amplitude, window when setting up a plurality of hyperbolic curve according to stack velocity, window and usable reflection axle basically identical when making this a plurality of hyperbolic curve, these a plurality of hyperbolic curve window lengths are the 200-500 millisecond, and are overlapped between the window in the time of up and down, overlapping scope 30%-50%, the geological data that step 1 is generated is divided into a plurality of minutes frequency band data, then the data timesharing window that contains noise band is carried out the amplitude energy statistics by following formula:
In the formula
EBe energy,
ABe amplitude,
iBe the sample period,
kBe Taoist monastic name,
fBe frequency band,
tFor the time window,
NBe total sample, medium filtering road number is set
nWith threshold value, according to the time in the window statistical value with
nThe average of road statistics compares, and just uses medium filtering compacting anomalous amplitude noise when the ratio of statistical value and average surpasses threshold value.
In step 3, when the geological data that step 2 is generated carries out the three-dimensional prestack random noise decay of F-X-Y territory, with the CMP axle of every line as the longitudinal axis, the single track of concentrating with the CMP road of the big minispread of offset distance is as transverse axis, form the plane of both direction and constitute 3-D data volume with the time, this 3-D data volume is done the Fourier transform of time orientation, obtain the data in F-X-Y territory, for each frequency content, utilize rectangle predictive operator of plural least square principle design, make corresponding to the 3D seismic data of this frequency content after this rectangle predictive operator prediction the prediction residual energy of its output
Q(
f) minimum, be about to
Q(
f) this rectangle predictive operator is asked local derviation and made it is zero, the matrix equation that can to obtain a coefficient be the Hermite matrix
R'·P=R
Wherein,
R'Be the autocorrelative Hermite matrix of F-X-Y territory geological data multiple tracks,
RBe F-X-Y territory geological data multiple tracks auto-correlation array,
PIt is the array of each component of this rectangle predictive operator, solve each component of this rectangle predictive operator, obtain this rectangle predictive operator again corresponding to this frequency content, with this rectangle predictive operator the data in the F-X-Y territory of this frequency content are carried out plural two-dimentional convolution, just obtained the result in F-X-Y territory after the three-dimensional denoising, the anti-Fourier transform in working frequency territory obtains three-dimensional random sound attenuation result.
In step 4, the geological data that step 3 is generated carries out the advantage spectrum analysis earlier, determine the dominant frequency band of reflection wave after, carry out frequency filtering again, the dominant frequency band energy of reflection wave is strengthened.
In step 4, when carrying out frequency filtering, the low frequency value of the corresponding dominant frequency band of low cut-off frequency of bandpass filter, the low frequency zone of transition adopts the 9dB raising frequency, the high frequency value of the corresponding dominant frequency band of high cut-off frequency of bandpass filter, high-frequency transition region adopts the 36dB frequency reducing.
The data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND among the present invention, at first from processed conventionally geological data, sub-elect the geological data that participates in velocity analysis, by optimizing sound attenuation technology parameter, improve the signal to noise ratio (S/N ratio) of velocity analysis data, adopt the wavelet shaping deconvolution to improve the resolution of velocity analysis road collection, utilize processing means such as frequency constraint, dominant frequency band reinforcement, strengthen the energy of usable reflection, improve the quality of velocity analysis data.The invention has the beneficial effects as follows, improve the signal to noise ratio (S/N ratio) and the resolution of velocity analysis road collection, the focusing of velocity spectrum is improved.
Description of drawings
Fig. 1 is the process flow diagram of the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND of the present invention;
Fig. 2 is the synoptic diagram of the road collection before and after the compound multiple domain sound attenuation in one embodiment of the invention;
Fig. 3 is the synoptic diagram of the road collection before and after the optimization process in one embodiment of the invention;
Fig. 4 is the synoptic diagram of the velocity spectrum before and after the optimization process in one embodiment of the invention.
Embodiment
For above-mentioned and other purposes, feature and advantage of the present invention can be become apparent, cited below particularlyly go out preferred embodiment, and cooperate appended graphicly, be described in detail below.
As shown in Figure 1, Fig. 1 is the process flow diagram of the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND of the present invention.In step 101, input is through processed conventionally geological data.Flow process enters into step 102.
In step 102, the geological data that step 101 is imported carries out the collection sorting of speed road, sub-elects the channel set that is used for velocity analysis.Flow process enters into step 103.
In step 103, the geological data that step 102 is generated carries out branch frequency band scanning analysis, determines the noise band scope, carries out the decay of branch frequency band anomalous amplitude again.Divide frequency band anomalous amplitude damped system as follows:
Window when setting up a plurality of hyperbolic curve according to stack velocity, window and usable reflection axle basically identical when making, window length 200-500 millisecond, overlapped between the window in the time of up and down, overlapping scope 30%-50%.Geological data is divided into a plurality of minutes frequency band data, then the data timesharing window that contains noise band is carried out the amplitude energy statistics by following formula:
,
In the formula
EBe energy,
ABe amplitude,
iBe the sample period,
kBe Taoist monastic name,
fBe frequency band,
tFor the time window,
NBe total sample.Medium filtering road number is set
nWith threshold value, according to the time in the window statistical value with
nThe average of road statistics compares, and just uses medium filtering compacting anomalous amplitude noise when the ratio of statistical value and average surpasses threshold value.Flow process enters into step 104.
In step 104, the geological data that step 103 is generated carries out the three-dimensional prestack random noise decay in F-X-Y territory.As the longitudinal axis, and the single track of concentrating with the CMP road of the big minispread of offset distance forms the plane of both direction and constitutes 3-D data volume with the time as transverse axis with the CMP axle of every line.Data volume is done the Fourier transform of time orientation hereto, obtains the data in F-X-Y territory.For each frequency content, can utilize rectangle predictive operator of plural least square principle design, make 3D seismic data corresponding to this frequency content after this rectangle predictive operator prediction, the prediction residual energy of its output
Q(
f) minimum, be about to
Q(
f) the rectangle predictive operator is asked local derviation and made it is zero, the matrix equation that can to obtain a coefficient be the Hermite matrix
R'·P=R
Wherein:
R'It is the autocorrelative Hermite matrix of F-X-Y territory geological data multiple tracks;
RIt is F-X-Y territory geological data multiple tracks auto-correlation array;
PIt is the array of each component of rectangle predictive operator.
According to the characteristics of Hermite matrix, can solve each component of this rectangle predictive operator apace, obtain rectangle predictive operator again corresponding to this frequency content.With this operator the data in the F-X-Y territory of this frequency content are carried out plural two-dimentional convolution, just obtained the result in F-X-Y territory after the three-dimensional denoising.Then, the anti-Fourier transform in working frequency territory obtains final three-dimensional random sound attenuation result.Flow process enters into step 105.
In step 105, the geological data that step 104 is generated carries out the advantage spectrum analysis, and the dominant frequency band scope of determining reflection wave is that amplitude is-frequency band range of 18dB correspondence, carries out bandpass filtering again.When carrying out bandpass filtering, the low frequency value of the corresponding dominant frequency band of low cut-off frequency of bandpass filter, the low frequency zone of transition adopts the 9dB raising frequency, the high frequency value of the corresponding dominant frequency band of high cut-off frequency of bandpass filter, high-frequency transition region adopts the 36dB frequency reducing.Flow process enters into step 106.
In step 106, the geological data that step 105 is generated carries out wavelet shaping deconvolution processing.The desired output of deconvolution is a zero-phase wavelet, selects time delay by experiment, adopts the high conduct output result of resolution.Flow process finishes.
Arrive shown in Figure 4 as Fig. 2, Fig. 2 is the synoptic diagram of the road collection before and after the compound multiple domain sound attenuation in one embodiment of the invention, Fig. 3 is the synoptic diagram of the road collection before and after the optimization process in one embodiment of the invention, and Fig. 4 is the synoptic diagram of the velocity spectrum before and after the optimization process in one embodiment of the invention.From Fig. 2 to Fig. 4 as can be seen, by analyzing the seismic data characteristic analysis, after technological synthesiss such as the compound multiple domain sound attenuation of The data (as shown in Figure 2), the dominant frequency band energy to the participation velocity analysis strengthened, raising resolution are handled, the quality of road collection is obviously improved (as shown in Figure 3), the velocity spectrum quality that forms obviously improves, the energy group of velocity spectrum focuses on, the focal region downward continuation more than 1 second, improved the precision (as shown in Figure 4) of SEISMIC VELOCTTY ANALYSIS AND.
Claims (7)
1. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND is characterized in that, the data-optimized method of this oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND comprises:
Step 1 is carried out the collection sorting of speed road to geological data, sub-elects the channel set that is used for velocity analysis;
Step 2, the geological data that step 1 is generated carries out the decay of branch frequency band anomalous amplitude;
Step 3, the geological data that step 2 is generated carries out the three-dimensional prestack random noise decay in F-X-Y territory;
Step 4, the geological data that step 3 is generated carries out frequency filtering, and the dominant frequency band energy of reflection wave is strengthened; And
Step 5, the geological data that step 4 is generated carries out wavelet shaping deconvolution processing, improves resolution.
2. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 1 is characterized in that, the data-optimized method of this oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND also is included in before the step 1, and input is through processed conventionally this geological data.
3. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 1, it is characterized in that in step 2, the geological data that step 1 is generated carries out frequency division band scanning analysis earlier, after determining the noise band scope, carry out the decay of branch frequency band anomalous amplitude again.
4. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 3, it is characterized in that, in step 2, when carrying out the decay of branch frequency band anomalous amplitude, window when setting up a plurality of hyperbolic curve according to stack velocity, window and usable reflection axle basically identical when making this a plurality of hyperbolic curve, these a plurality of hyperbolic curve window lengths are the 200-500 millisecond, overlapped between the window in the time of up and down, overlapping scope 30%-50%, the geological data that step 1 is generated is divided into a plurality of minutes frequency band data, then the data timesharing window that contains noise band is carried out the amplitude energy statistics by following formula:
,
In the formula
EBe energy,
ABe amplitude,
iBe the sample period,
kBe Taoist monastic name,
fBe frequency band,
tFor the time window,
NBe total sample, medium filtering road number is set
nWith threshold value, according to the time in the window statistical value with
nThe average of road statistics compares, and just uses medium filtering compacting anomalous amplitude noise when the ratio of statistical value and average surpasses threshold value.
5. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 1, it is characterized in that, in step 3, when the geological data that step 2 is generated carries out the three-dimensional prestack random noise decay of F-X-Y territory, with the CMP axle of every line as the longitudinal axis, the single track of concentrating with the CMP road of the big minispread of offset distance forms the plane of both direction and constitutes 3-D data volume with the time, to this three-dimensional data as transverse axis
Body is done the Fourier transform of time orientation, obtain the data in F-X-Y territory,, utilize rectangle predictive operator of plural least square principle design for each frequency content, make corresponding to the 3D seismic data of this frequency content after this rectangle predictive operator prediction the prediction residual energy of its output
Q(
f) minimum, be about to
Q(
f) this rectangle predictive operator is asked local derviation and made it is zero, the matrix equation that can to obtain a coefficient be the Hermite matrix
R'·P=R
Wherein,
R'Be the autocorrelative Hermite matrix of F-X-Y territory geological data multiple tracks,
RBe F-X-Y territory geological data multiple tracks auto-correlation array,
PIt is the array of each component of this rectangle predictive operator, solve each component of this rectangle predictive operator, obtain this rectangle predictive operator again corresponding to this frequency content, with this rectangle predictive operator the data in the F-X-Y territory of this frequency content are carried out plural two-dimentional convolution, just obtained the result in F-X-Y territory after the three-dimensional denoising, the anti-Fourier transform in working frequency territory obtains three-dimensional random sound attenuation result.
6. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 1, it is characterized in that, in step 4, the geological data that step 3 is generated carries out the advantage spectrum analysis earlier, after determining the dominant frequency band of reflection wave, carry out frequency filtering again, the dominant frequency band energy of reflection wave is strengthened.
7. the data-optimized method of oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND according to claim 6, it is characterized in that, in step 4, when carrying out frequency filtering, the low frequency value of the corresponding dominant frequency band of low cut-off frequency of bandpass filter, the low frequency zone of transition adopts the 9dB raising frequency, the high frequency value of the corresponding dominant frequency band of high cut-off frequency of bandpass filter, and high-frequency transition region adopts the 36dB frequency reducing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210017116.7A CN103217713B (en) | 2012-01-19 | 2012-01-19 | Oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND data optimization methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210017116.7A CN103217713B (en) | 2012-01-19 | 2012-01-19 | Oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND data optimization methods |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103217713A true CN103217713A (en) | 2013-07-24 |
CN103217713B CN103217713B (en) | 2015-08-05 |
Family
ID=48815665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210017116.7A Active CN103217713B (en) | 2012-01-19 | 2012-01-19 | Oil-gas exploration SEISMIC VELOCTTY ANALYSIS AND data optimization methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103217713B (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954993A (en) * | 2014-04-17 | 2014-07-30 | 中国石油天然气集团公司 | Scale domain multichannel filtering method and device based on seismic signal continuous wavelet transformation |
CN104133239A (en) * | 2014-07-29 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for defining earthquake single-shot record target stratum time window |
CN105319586A (en) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | Velocity spectrum three-dimensional filtering optimization velocity analytical method |
CN106125134A (en) * | 2016-07-04 | 2016-11-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Based on the geological data signal-noise ratio computation method of window during hyperbolic |
CN106338769A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Seismic data denoising method and system |
CN107807390A (en) * | 2016-09-09 | 2018-03-16 | 中国石油化工股份有限公司 | The processing method and system of geological data |
CN109116425A (en) * | 2018-10-31 | 2019-01-01 | 中国石油化工股份有限公司 | Utilize the method for the frequency spectrum design filter removal noise of back wave |
CN109212603A (en) * | 2018-11-08 | 2019-01-15 | 北京邮电大学 | A kind of Earthquake signal detection algorithm based on GS transformation filtering and EMD denoising |
CN109343119A (en) * | 2018-11-10 | 2019-02-15 | 北京科胜伟达石油科技股份有限公司 | The main energy relative amplitude preserved processing method of low signal-to-noise ratio seismic data useful signal |
CN110873899A (en) * | 2018-08-30 | 2020-03-10 | 中国石油化工股份有限公司 | Seismic data super gather generation method based on offset distance grouping |
CN110895351A (en) * | 2018-09-12 | 2020-03-20 | 中国石油化工股份有限公司 | Method for suppressing secondary interference |
CN114002736A (en) * | 2021-09-07 | 2022-02-01 | 中国矿业大学 | Seismic exploration multi-frequency data fusion method based on weight deconvolution |
CN114089420A (en) * | 2020-08-24 | 2022-02-25 | 中国石油化工集团有限公司 | Method for quantitatively detecting crosstalk interference |
CN116388696A (en) * | 2023-06-06 | 2023-07-04 | 苏州华兴源创科技股份有限公司 | Noise signal generation method and device |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2230861A (en) * | 1989-03-06 | 1990-10-31 | Amoco Corp | Velocity analysis in seismic data processing |
CN101419293A (en) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | Method for enhancing earthquake data signal-to-noise ration |
CN101852863A (en) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | Method for processing seismic data by using high-precision single-channel spectrum analysis technology |
CN102096103A (en) * | 2010-12-03 | 2011-06-15 | 中国石油天然气集团公司 | Velocity analysis method for low signal-to-noise ratio seismic data |
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
-
2012
- 2012-01-19 CN CN201210017116.7A patent/CN103217713B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2230861A (en) * | 1989-03-06 | 1990-10-31 | Amoco Corp | Velocity analysis in seismic data processing |
US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
CN101419293A (en) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | Method for enhancing earthquake data signal-to-noise ration |
CN101852863A (en) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | Method for processing seismic data by using high-precision single-channel spectrum analysis technology |
CN102096103A (en) * | 2010-12-03 | 2011-06-15 | 中国石油天然气集团公司 | Velocity analysis method for low signal-to-noise ratio seismic data |
Non-Patent Citations (2)
Title |
---|
刘晓峰等: "三维AVO处理方法及应用", 《油气地质与采收率》, vol. 10, no. 03, 25 June 2003 (2003-06-25), pages 36 - 38 * |
车廷信等: "地震数据处理中的高频速度分析研究", 《世界地质》, vol. 22, no. 02, 30 June 2003 (2003-06-30), pages 191 - 195 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954993A (en) * | 2014-04-17 | 2014-07-30 | 中国石油天然气集团公司 | Scale domain multichannel filtering method and device based on seismic signal continuous wavelet transformation |
CN105319586A (en) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | Velocity spectrum three-dimensional filtering optimization velocity analytical method |
CN105319586B (en) * | 2014-07-04 | 2018-02-02 | 中国石油化工股份有限公司 | A kind of normal-moveout spectrum three-dimensional filtering optimal velocity analysis |
CN104133239A (en) * | 2014-07-29 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for defining earthquake single-shot record target stratum time window |
CN106338769A (en) * | 2015-07-07 | 2017-01-18 | 中国石油化工股份有限公司 | Seismic data denoising method and system |
CN106125134A (en) * | 2016-07-04 | 2016-11-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Based on the geological data signal-noise ratio computation method of window during hyperbolic |
CN107807390B (en) * | 2016-09-09 | 2019-08-23 | 中国石油化工股份有限公司 | The processing method and system of seismic data |
CN107807390A (en) * | 2016-09-09 | 2018-03-16 | 中国石油化工股份有限公司 | The processing method and system of geological data |
CN110873899A (en) * | 2018-08-30 | 2020-03-10 | 中国石油化工股份有限公司 | Seismic data super gather generation method based on offset distance grouping |
CN110873899B (en) * | 2018-08-30 | 2021-12-28 | 中国石油化工股份有限公司 | Seismic data super gather generation method based on offset distance grouping |
CN110895351A (en) * | 2018-09-12 | 2020-03-20 | 中国石油化工股份有限公司 | Method for suppressing secondary interference |
CN109116425B (en) * | 2018-10-31 | 2021-02-26 | 中国石油化工股份有限公司 | Method for removing noise by using frequency spectrum design filter of reflected wave |
CN109116425A (en) * | 2018-10-31 | 2019-01-01 | 中国石油化工股份有限公司 | Utilize the method for the frequency spectrum design filter removal noise of back wave |
CN109212603A (en) * | 2018-11-08 | 2019-01-15 | 北京邮电大学 | A kind of Earthquake signal detection algorithm based on GS transformation filtering and EMD denoising |
CN109343119A (en) * | 2018-11-10 | 2019-02-15 | 北京科胜伟达石油科技股份有限公司 | The main energy relative amplitude preserved processing method of low signal-to-noise ratio seismic data useful signal |
CN114089420A (en) * | 2020-08-24 | 2022-02-25 | 中国石油化工集团有限公司 | Method for quantitatively detecting crosstalk interference |
CN114089420B (en) * | 2020-08-24 | 2024-05-24 | 中国石油化工集团有限公司 | Method for quantitatively detecting crosstalk interference |
CN114002736A (en) * | 2021-09-07 | 2022-02-01 | 中国矿业大学 | Seismic exploration multi-frequency data fusion method based on weight deconvolution |
CN116388696A (en) * | 2023-06-06 | 2023-07-04 | 苏州华兴源创科技股份有限公司 | Noise signal generation method and device |
CN116388696B (en) * | 2023-06-06 | 2023-08-15 | 苏州华兴源创科技股份有限公司 | Noise signal generation method and device |
Also Published As
Publication number | Publication date |
---|---|
CN103217713B (en) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103217713A (en) | Method for optimizing seismic velocity analysis data in oil and gas exploration | |
CN100349011C (en) | Method for suppressing background noise not related to excitation source in seismic data enhancement | |
CN101598803B (en) | Method for directly obtaining stacked section of converted wave | |
CN108037526A (en) | Reverse-time migration method based on all-wave wave field VSP/RVSP seismic datas | |
CN101382598B (en) | Pressing method for true 3-d seismics data linear noise | |
CN102109612A (en) | Seismic wave absorption and attenuation compensation method | |
CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
CN109633752B (en) | Offshore towing cable data self-adaptive ghost wave compression method based on three-dimensional fast Radon transformation | |
CN102692647B (en) | Stratum oil-gas possibility prediction method with high time resolution | |
CN105445801B (en) | A kind of processing method for eliminating 2-d seismic data random noise | |
CN103605157B (en) | Method for attenuating near-surface scattered waves | |
CN104808245A (en) | Gather optimized processing method and device thereof | |
CN103728660A (en) | Multi-channel matching tracking method based on seismic data | |
CN104345341A (en) | Region constraint-based frequency band division energy seismic surface wave processing method | |
CN104142518A (en) | Method for analyzing pre-stack time migration response of seismological observation system | |
CN103645497A (en) | Emergence angle based multi-component wave field separation method | |
CN101852864B (en) | Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology | |
CN103576197A (en) | Method for extracting converted wave angle channel set | |
CN103076626B (en) | A kind of wave field purification treatment method | |
CN102103215B (en) | Method for suppressing surface waves of three-dimensional high-density seismic prospecting records before stack | |
CN104635264B (en) | Pre-stack seismic data processing method and device | |
CN106873031B (en) | A kind of 3 D seismic observation system vertical resolution quantitative analysis evaluation method | |
CN103592684B (en) | A kind of massive seismic data compression method keeping space attribute information and device | |
CN105093282A (en) | Energy replacement surface wave suppressing method based on frequency constraint | |
CN104133244B (en) | A kind of near surface static corrections processing method of focus shear wave exploration |
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 |