CN102540252B - High-precision median stacking method on basis of cross-correlation - Google Patents

High-precision median stacking method on basis of cross-correlation Download PDF

Info

Publication number
CN102540252B
CN102540252B CN201110422335.9A CN201110422335A CN102540252B CN 102540252 B CN102540252 B CN 102540252B CN 201110422335 A CN201110422335 A CN 201110422335A CN 102540252 B CN102540252 B CN 102540252B
Authority
CN
China
Prior art keywords
sampling point
stacking
collection
value
trace
Prior art date
Application number
CN201110422335.9A
Other languages
Chinese (zh)
Other versions
CN102540252A (en
Inventor
罗红明
何光明
陈爱萍
曹中林
张华�
刘鸿
吕文彪
白静
Original Assignee
中国石油集团川庆钻探工程有限公司地球物理勘探公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 中国石油集团川庆钻探工程有限公司地球物理勘探公司 filed Critical 中国石油集团川庆钻探工程有限公司地球物理勘探公司
Priority to CN201110422335.9A priority Critical patent/CN102540252B/en
Publication of CN102540252A publication Critical patent/CN102540252A/en
Application granted granted Critical
Publication of CN102540252B publication Critical patent/CN102540252B/en

Links

Abstract

The invention provides a high-precision median stacking method on the basis of cross-correlation. The method comprises the following steps: (1) stacking and filtering a trace gather subjected to dynamic correction treatment to obtain a model trace; (2) calculating a cross-correlation time difference of the model trace and one seismic trace of the current trace gather and correcting the time difference; (3) calculating a cross-correlation time difference of the model trace and a next seismic trace of the current gather and correcting the time difference until the time difference correction of all the seismic traces of the current gather is completed; (4) removing abnormal values in sampling points which take part in stacking, calculating an arithmetic mean value on the residual sampling points and removing the sampling point values of which symbols are opposite to those of the arithmetic mean value; (5) generating a weight of each sampling point which takes part in stacking in the residual sampling points; (6) unitizing the weight of each sampling point of the current gather, which takes part in stacking; (7) according to the unitized weights, completing the weighted stacking of each sampling point; and (8) repeatedly executing the steps (1) to (7) and carrying out weighted stacking of each sampling point of a next gather until the weighted stacking of each sampling point of all gathers is completed.

Description

High precision intermediate value stacking method based on simple crosscorrelation

Technical field

The present invention relates to petroleum gas acquiring seismic exploration data field, specifically, relate to a kind of high precision intermediate value stacking method based on simple crosscorrelation, the field seismic data that is mainly used in petroleum gas seismic prospecting is processed explanation.

Background technology

Complex area stacking image method research is the medium-term and long-term existence of oil earth rock gas physical prospecting in the frontier nature research topic solving not yet completely at present, and research has at present obtained a lot of theories and practical application achievement.Current research concentrates on following two aspects mostly:

1, how research obtains the more preferably method of normal moveout correction road collection

In seismic data processing practically, the lineups of the same reflection point of normal moveout correction road collection generally cannot strictly be smoothed, often by considering that coefficient of anisotropy or high-order term realize the road collection of offset distance far away and even up, but in seismic data processing practically, because the factor on earth's surface cannot be eliminated completely, and the horizontal change of shallow-layer speed makes nearly offset distance or middle offset distance still can have the phenomenon that cannot smooth conventionally, these means of shaking impenetrably propagation realize smoothing, so by tectonic model road, the correlation theory of utilizing mathematics time realize the processing of residual static correction in window, obtain more desirable road collection.

Above-mentioned research mainly comprises two parts content: the one, improve the computational accuracy of stack velocity, main method is to utilize anisotropy velocity analysis to obtain stack velocity and anisotropic parameter etc., thereby realize normal moveout correction road collection in, offset distance far away correctly playbacks; The 2nd, improve the precision that normal moveout correction is calculated, as utilize high-order normal moveout correction method etc. realize moving correcting track collection in, offset distance far away correctly playbacks.Because the influence factor of complex structural area seismic velocity is a lot, although said method has improved the precision of velocity analysis to a certain extent, can not ensure to obtain accurate stacking velocity field.

2, utilize stacking method more flexibly to realize the high precision imaging to road collection

These class methods are main study hotspots at present, have obtained some achievements in research, have obtained application comparatively widely.For example, but the research of this respect mainly still respectively lays particular emphasis on some aspects, and, Coherent addition is just considered the correlation computations of collection, the rejecting to exceptional value is just considered in smart stacking.These methods can only be exported the data after stack, cannot take into account in many aspects simultaneously, therefore, conventionally have to accept or reject to some extent, will inevitably affect like this precision of imaging.Even if said method adopts weighting processing, its weights calculate also very complicated, apply inaccurately, and do not consider that some seismic trace that participates in stack should reject, and can more be conducive to the stacking image of section after rejecting.

Although anisotropy normal moveout correction and high-order normal moveout correction method have advantages of that precision is high, they are not suitable for complex structure imaging, and existing stacking method function singleness, has little significance for complex area practical operation.

Summary of the invention

The problem existing for prior art, the invention provides a kind of high precision intermediate value stacking method based on simple crosscorrelation, described method comprises: (1) superposes to normal moveout correction road after treatment collection, and the data to stack are carried out filtering, the model trace using filtered data as the seismic trace correction calculation for current road collection; (2) when given in window, the simple crosscorrelation time difference between the seismic trace that the participation of computation model road and current road collection is calculated, and the seismic trace that described participation is calculated carries out TEC time error correction; (3) the simple crosscorrelation time difference between next seismic trace that the participation of computation model road and current road collection is calculated, and next seismic trace of described participation calculating is carried out to TEC time error correction, until complete the TEC time error correction of all seismic traces of the participation calculating of current road collection; (4) reject the exceptional value in the sampling point that participates in stack, and remaining sampling point is asked to arithmetic mean, reject the sample value contrary with arithmetic mean symbol; (5) in remaining sampling point, produce the weights that participate in each sampling point superposeing; (6) weights of each sampling point of the participation stack to current road collection are normalized; (7), according to the weights after the normalization of calculating, complete the weighted stacking of each sampling point; (8) repeated execution of steps (1), to step (7), is descended the weighted stacking of each sampling point of collection together, until complete the weighted stacking of each sampling point of all roads collection, to realize the accurately image of underground complex structure.

In step (1), normal moveout correction road after treatment collection can be common midpoint gather.

In step (5), can adopt the anti-distance weighted mode of sampling point amplitude and average, in remaining sampling point, produce the weights of each sampling point that participates in stack.

Can calculate according to formula below the weights of each sampling point:

W i = 1 ( x i - A ) n

Wherein, the arithmetic mean that A is sampling point, x ifor current sample value, n is index, the sequence number that i is each sampling point.

In step (4), the exceptional value in sampling point can be the extreme value in sampling point.

High precision intermediate value stacking method based on simple crosscorrelation according to the present invention is processed and imaging applicable to the structure of Complex Mountain, there is the features such as easy, counting yield is high and imaging effect is good of calculating, in Complex Mountain seismic data is processed, have broad application prospects.

Brief description of the drawings

In conjunction with the drawings, from the description of the following examples, the present invention these and/or other side and advantage will become clear, and are easier to understand, wherein:

Fig. 1 is according to the process flow diagram of the high precision intermediate value stacking method based on simple crosscorrelation of the present invention.

Embodiment

In the present invention, normal moveout correction common midpoint gather after treatment is superposeed and obtains collection data, then these data are carried out to filtering processing, the model trace using the data that obtain as the seismic trace correction calculation for current road collection.By calculating the simple crosscorrelation time difference between the seismic trace (can be described as and calculate road) that the participation of this model trace and current road collection calculates in window when given, the seismic trace that participates in calculating is carried out to TEC time error correction.Equally, other seismic trace that the participation of current road collection is calculated also adopts this given time window, and the calculating simple crosscorrelation time difference is also proofreaied and correct the simple crosscorrelation time difference, thereby completes the TEC time error correction of all seismic traces that the participation of current road collection is calculated.

Then, adopt the mode of smart stacking, reject the exceptional value in the sampling point that participates in stack, and remaining sampling point is asked to arithmetic mean, reject the sample value contrary with mean value symbol.In remaining sampling point, (employing sampling point amplitude and average) anti-distance weighted mode, produces the weights that participate in each sampling point superposeing, and finally realizes weighted stacking.

Can descend the TEC time error correction of collection together by similar mode, until complete the TEC time error correction of all roads collection, thus realize the more accurate imaging of underground complex structure.

Below, describe the high precision intermediate value stacking method based on simple crosscorrelation in detail with reference to Fig. 1.Fig. 1 is according to the process flow diagram of the high precision intermediate value stacking method based on simple crosscorrelation of the present invention.

With reference to Fig. 1, in step 101, normal moveout correction road after treatment collection is superposeed, and the data of stack are carried out to filtering, the model trace using filtered data as the seismic trace correction calculation for current road collection.Preferably, normal moveout correction road after treatment collection is common midpoint gather.

In step 102, when given in window, the simple crosscorrelation time difference between the seismic trace that the participation of computation model road and current road collection is calculated, and the seismic trace that described participation is calculated carries out TEC time error correction.

In step 103, the simple crosscorrelation time difference between next seismic trace that the participation of computation model road and current road collection is calculated, and next seismic trace that described participation is calculated carries out TEC time error correction.

In step 104, determine whether the TEC time error correction of all seismic traces of the participation calculating of current road collection.

If do not complete the TEC time error correction of all seismic traces of the participation calculating of current road collection, turn back to step 103, proceed the TEC time error correction of other seismic trace of the participation calculating of current road collection.

The TEC time error correction of all seismic traces that the participation of current road collection is calculated if be over, in step 105, rejects the exceptional value in the sampling point that participates in stack, and remaining sampling point is asked to arithmetic mean, rejects the sample value contrary with arithmetic mean symbol.Preferably, the exceptional value in sampling point can be the extreme value in sampling point.

In step 106, in remaining sampling point, produce the weights of each sampling point that participates in stack.Preferably, can adopt the anti-distance weighted mode of sampling point amplitude and average, in remaining sampling point, produce the weights of each sampling point that participates in stack.

The weights computing formula of each sampling point is as follows:

W i = 1 ( x i - A ) n

Wherein, the arithmetic mean that A is sampling point, x ifor current sample value, n is index, the sequence number that i is each sampling point.

In step 107, the weights of each sampling point to the participation stack with collection (current road collection) are normalized.

In step 108, according to the weights after the normalization of calculating, complete the weighted stacking of each sampling point.

In step 109, repeated execution of steps 101-108, descends the weighted stacking of each sampling point of collection together.

In step 110, determine whether the weighted stacking of each sampling point of all roads collection.

If do not complete the weighted stacking of each sampling point of all roads collection, turn back to step 101, proceed the weighted stacking of each sampling point of other road collection.

If completed the weighted stacking of each sampling point of all roads collection, the method stops, and finally realizes the weighted stacking of each sampling point of all roads collection, to realize the accurately image of underground complex structure.

For example, for a series of sample value, if each sample value equals the mean value of these a series of sample value, the superposition value of these sample value just equals its sample value.For example, if there are four sample value: 5,5,5,5, its superposition value M=(5+5+5+5)/4=5; For example, if wherein certain sample value is exceptional value (four sample value is respectively 5,5,9,5), its superposition value M '=(5+5+9+5)/4=6, therefore can not reflect original actual value M=5, now calculate the weights (wherein n=2) of each sampling point according to formula described above:

W i = 1 ( x i - A ) n

The weights of each sampling point calculating are respectively 1,1,1/9,1, its normalized value is respectively 9/28,9/28,1/28,9/28, therefore the weighted stacking value of these sampling points is: 5 × 9/28+5 × 9/28+9 × 1/28+5 × 9/28=5.142, more approaches actual value 5 than M '=6 like this.Therefore,, even also there are indivedual exceptional values after by foregoing step 105 rejecting abnormalities value, also can, by the interference of the weighted stacking rejecting abnormalities value of each sampling point, be convenient to realize the more accurate imaging of underground complex structure.

High precision intermediate value stacking method based on simple crosscorrelation according to the present invention can be had and had the following advantages:

1, adopt directly stack production model road, result of calculation is more accurate, and counting yield also improves greatly;

2, utilize model trace and calculate the calculating of the simple crosscorrelation time difference and the mode of proofreading and correct between road, computational accuracy is also largely increased.

High precision intermediate value stacking method based on simple crosscorrelation according to the present invention is specially adapted to the structure of Complex Mountain and processes and imaging, there is the features such as easy, counting yield is high and imaging effect is good of calculating, in Complex Mountain seismic data is processed, have broad application prospects.

Although the present invention is described in detail and shows with reference to its exemplary embodiment, but will be understood by those skilled in the art that, in the case of not departing from the spirit and scope of the present invention that are defined by the claims, can carry out to it various changes of form and details.

Claims (3)

1. the high precision intermediate value stacking method based on simple crosscorrelation, comprising:
(1) normal moveout correction road after treatment collection is superposeed, and the data of stack are carried out to filtering, the model trace using filtered data as the seismic trace correction calculation for current road collection;
(2) when given in window, the simple crosscorrelation time difference between the seismic trace that the participation of computation model road and current road collection is calculated, and the seismic trace that described participation is calculated carries out TEC time error correction;
(3) the simple crosscorrelation time difference between next seismic trace that the participation of computation model road and current road collection is calculated, and next seismic trace of described participation calculating is carried out to TEC time error correction, until complete the TEC time error correction of all seismic traces of the participation calculating of current road collection;
(4) reject the exceptional value in the sampling point that participates in stack, and remaining sampling point is asked to arithmetic mean, reject the sample value contrary with arithmetic mean symbol;
(5) in remaining sampling point, produce the weights that participate in each sampling point superposeing;
(6) weights of each sampling point of the participation stack to current road collection are normalized;
(7), according to the weights after the normalization of calculating, complete the weighted stacking of each sampling point;
(8) repeated execution of steps (1), to step (7), is descended the weighted stacking of each sampling point of collection together, until complete the weighted stacking of each sampling point of all roads collection, to realize the accurately image of underground complex structure,
Wherein, in step (5), adopt the anti-distance weighted mode of sampling point amplitude and average, in remaining sampling point, produce the weights of each sampling point that participates in stack,
Wherein, calculate the weights of each sampling point according to formula below:
W i = 1 ( x i - A ) n
Wherein, the arithmetic mean that A is sampling point, x ifor current sample value, n is index, the sequence number that i is each sampling point.
2. high precision intermediate value stacking method according to claim 1, wherein, in step (1), normal moveout correction road after treatment collection is common midpoint gather.
3. high precision intermediate value stacking method according to claim 2, wherein, in step (4), the exceptional value in sampling point is the extreme value in sampling point.
CN201110422335.9A 2011-12-15 2011-12-15 High-precision median stacking method on basis of cross-correlation CN102540252B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110422335.9A CN102540252B (en) 2011-12-15 2011-12-15 High-precision median stacking method on basis of cross-correlation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110422335.9A CN102540252B (en) 2011-12-15 2011-12-15 High-precision median stacking method on basis of cross-correlation

Publications (2)

Publication Number Publication Date
CN102540252A CN102540252A (en) 2012-07-04
CN102540252B true CN102540252B (en) 2014-08-06

Family

ID=46347562

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110422335.9A CN102540252B (en) 2011-12-15 2011-12-15 High-precision median stacking method on basis of cross-correlation

Country Status (1)

Country Link
CN (1) CN102540252B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3324216A4 (en) * 2016-04-01 2018-10-17 Petrochina Company Limited Method and device for processing seismic data

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103728665B (en) * 2012-10-12 2016-06-08 中国石油化工股份有限公司 A kind of fine residual static correction method
CN103399346B (en) * 2013-08-15 2016-04-06 电子科技大学 A kind of well shake associating impedance initial value modeling method
CN104155688A (en) * 2014-08-13 2014-11-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 High precision weighted stack method
CN104181588B (en) * 2014-08-15 2016-08-31 中国石油集团川庆钻探工程有限公司地球物理勘探公司 A kind of method building high-quality model trace
CN105425295B (en) * 2015-12-18 2017-11-14 中国石油集团东方地球物理勘探有限责任公司 Geological data is superimposed weights acquisition methods and stage by stage stacking method stage by stage
CN109254322B (en) * 2017-07-14 2020-08-11 中国石油天然气股份有限公司 Method for establishing residual static correction model channel based on heredity

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5995907A (en) * 1998-02-05 1999-11-30 Geoquest Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5995907A (en) * 1998-02-05 1999-11-30 Geoquest Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
互相关时差分析技术及其应用-以胜利油田三维地震资料连片处理为例;李继光等;《石油物探》;20100131;第49卷(第1期);全文 *
李继光等.互相关时差分析技术及其应用-以胜利油田三维地震资料连片处理为例.《石油物探》.2010,第49卷(第1期),23-29.
横向寻优叠加提高地震资料成像精度;葛川庆等;《勘探地球物理进展》;20050430;第28卷(第2期);全文 *
葛川庆等.横向寻优叠加提高地震资料成像精度.《勘探地球物理进展》.2005,第28卷(第2期),122-124.

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3324216A4 (en) * 2016-04-01 2018-10-17 Petrochina Company Limited Method and device for processing seismic data

Also Published As

Publication number Publication date
CN102540252A (en) 2012-07-04

Similar Documents

Publication Publication Date Title
CN102508293B (en) Pre-stack inversion thin layer oil/gas-bearing possibility identifying method
CN103713315B (en) A kind of seismic anisotropy parameter full waveform inversion method and device
CN103376464B (en) A kind of inversion method for stratigraphic quality factor
CN103034123B (en) Based on the parallel robot control method of kinetic parameters identification
CN102890290B (en) Prestack depth migration method under a kind of relief surface condition
CN103954295B (en) A kind of step-recording method based on acceleration transducer
CN102998704B (en) Geophysical exploration seismic data processing method
CN105589100B (en) A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method
CN102692645B (en) Method for performing joint inversion on P-wave and S-wave velocity ratio of reservoir by utilizing P-wave and converted wave data
CN100545679C (en) The earthquake poststack forward modeling method of relief surface
CN101251604B (en) Method for analyzing and NMO correcting two parameters transformation wave speed
CN102053261B (en) Method for processing seismic data
CN104599292B (en) A kind of anti-noise moving object detection algorithm decomposed based on low-rank matrix
Lebedev et al. Mapping the Moho with seismic surface waves: a review, resolution analysis, and recommended inversion strategies
CN104155691B (en) Converted wave anisotropic velocity analysis method and device
CN102914791B (en) Kirchhoff prestack time migration method for processing seismic data of undulating surface
CN104391323B (en) A kind of method utilizing lower wave number composition in reflected wave information inversion speed field
CN104237940A (en) Diffracted wave imaging method and diffracted wave imaging device based on dynamical features
CN104123463A (en) Time domain identification method of random dynamic loads
CN102841375A (en) Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition
CN103064115B (en) Method for matching longitudinal wave and converted wave of radial parameter region
CN103399300B (en) Wave packet superposition microseism ground location method
CN103323615A (en) Mobile terminal and method for calculating walking speed through acceleration sensor
CN101776768B (en) Anisotropy speed analysis and dynamic correction method
CN102937721B (en) Limited frequency tomography method for utilizing preliminary wave travel time

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180129

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

Effective date of registration: 20180129

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

TR01 Transfer of patent right

Effective date of registration: 20200918

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.