CN110389378A - Seismic data zero phase bearing calibration - Google Patents
Seismic data zero phase bearing calibration Download PDFInfo
- Publication number
- CN110389378A CN110389378A CN201910343233.4A CN201910343233A CN110389378A CN 110389378 A CN110389378 A CN 110389378A CN 201910343233 A CN201910343233 A CN 201910343233A CN 110389378 A CN110389378 A CN 110389378A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- seismic
- zero phase
- lineups
- bearing calibration
- 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
- 238000003908 quality control method Methods 0.000 claims abstract description 24
- 238000012937 correction Methods 0.000 claims abstract description 18
- 238000011160 research Methods 0.000 claims abstract description 4
- 238000003325 tomography Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 239000004058 oil shale Substances 0.000 description 1
- 230000005502 phase rule Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 239000009671 shengli Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of seismic data zero phase bearing calibration, which includes: step 1: selecting suitable seismic data lineups as phasing with reference to lineups on research work area seismic profile;Step 2: selecting suitable Quality Control seismic profile or seismic data cube;Step 3: the corresponding lineups equalized amplitude energy of the different correction angles of calculating;Step 4: for the equalized amplitude energy recorded in step 3, selecting the corresponding angle of maximum equalized amplitude energy as the zero phase of seismic data cube and correct angle;Step 5: the zero phase calculated in step 4 correction angle being applied to seismic data cube, to carry out the zero phase correction of seismic data.The seismic data zero phase bearing calibration clear principle, it is easy to operate, it is particularly suitable for seismic geology interpretation technical staff.
Description
Technical field
The present invention relates to Seismic Data Processing Technique fields, especially relate to a kind of seismic data zero phase correction side
Method.
Background technique
Not only resolution ratio with higher and the seismic data calibration and Seismic Sedimentology explanation of zero phase seismic data
Basis.Currently used seismic data dephasing processing technique and means, generally all require in the seism processing stage
Wavelet extraction and deconvolution processing are carried out, processing method and process are complicated, and quality control difficulty is big, and uncertain factor is more, ground
Shake data interpretation personnel are difficult to grasp.Thus we have invented a kind of new seismic data zero phase bearing calibration, solve with
Upper technical problem.
Summary of the invention
Most collect the object of the present invention is to provide energy when one kind is based on seismic data zero phase and neutralizes the highest rule of resolution ratio
Rule, may be implemented the seismic data zero phase bearing calibration of the dephasing of seismic data.
The purpose of the present invention can be achieved by the following technical measures: seismic data zero phase bearing calibration, earthquake money
Material zero phase bearing calibration includes: step 1: selecting suitable seismic data lineups as phase on research work area seismic profile
Bit correction refers to lineups;Step 2: selecting suitable Quality Control seismic profile or seismic data cube;Step 3: calculating different corrections
The corresponding lineups equalized amplitude energy of angle;Step 4: for the equalized amplitude energy recorded in step 3, selecting maximum balanced
The corresponding angle of amplitude energy corrects angle as the zero phase of seismic data cube;Step 5: the zero phase that will be calculated in step 4
It corrects angle and is applied to seismic data cube, to carry out the zero phase correction of seismic data.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, seismic data lineups selection standard is to select contact relation simple, continuous-stable in transverse direction,
Change apparent lineups with upper and lower formation lithology or impedance characteristic.
In step 2, selected seismic profile or seismic data cube are to be mutated these far from tomography, stratal pinch-outs, lithology
Influence the special geobody of seismic data lineups stability.
In step 3, in step 2 select Quality Control seismic profile or seismic data cube within the scope of certain angle into
Row phasing, and record the Quality Control seismic profile of corresponding correction angle or Quality Control seismic data cube correspond to step 1 select it is same
To the equalized amplitude energy of axis.
In step 3, selected angle correct range is -90 degree to+90 degree.
Step 3 includes:
Step a: record Quality Control seismic profile or Quality Control seismic data cube each seismic channel correspond to step 1 select it is in the same direction
The peak swing energy of axis;
Step b: the peak swing of all seismic channels can be measured into arithmetic mean or take root mean square as the equal of the lineups
Weigh amplitude energy.
Seismic data zero phase bearing calibration in the present invention, according to " "ockham's razor" principle ", which is known as " such as nothing
Necessity does not increase entity ", i.e., " simple and effective principle ", when seismic data is zero phase, seismic reflection lineups energy most collects
In, Energy maximum value position and reflecting interface position consistency;The offset of phase will lead to the dispersion of reflected energy, lineups energy
It is subjected to displacement at amount maximum value with geology reflecting interface.It is highest most to collect neutralization resolution ratio for energy when based on seismic data zero phase
Rule, if the method that can be judged at the maximum position of lineups energy that can be referred to can be found in seismic data,
It can realize the dephasing of seismic data.With by extracting wavelet, seismic data dephasing is realized using convolution operator
It is compared etc. fine wavelet shaping method, the seismic data zero phase bearing calibration clear principle, easy to operate, data is accurate, easily
In realization, it is particularly suitable for seismic geology interpretation technical staff.
Detailed description of the invention
Fig. 1 is the flow chart of a specific embodiment of seismic data zero phase bearing calibration of the invention;
Fig. 2 is different angle Ricker wavelet in a specific embodiment of the invention to different type reservoir forward result figure;
Fig. 3 is the t6 lineups selected in a specific embodiment of the invention and the schematic diagram of phasing Quality Control section;
Fig. 4 is phase calibration and equalized amplitude energy value relational graph in a specific embodiment of the invention;
Fig. 5 is dephasing front and back contrast of seismic profile figure in a specific embodiment of the invention;
Fig. 6 is dephasing front and back strata slicing comparison diagram in a specific embodiment of the invention.
Specific embodiment
In order to which the technical problems, technical solutions and beneficial effects solved by the present invention is more clearly understood, below in conjunction with
Accompanying drawings and embodiments, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein is only used
To explain the present invention, it is not intended to limit the present invention.
As shown in FIG. 1, FIG. 1 is the flow charts of seismic data zero phase bearing calibration of the invention.
Step 101: selecting suitable seismic data lineups to refer to as phasing on research work area seismic profile
Lineups.Seismic data lineups selection standard is that selection contact relation is simple, continuous-stable and upper and lower formation rock in transverse direction
Property or impedance characteristic change apparent lineups.
Step 102: selecting suitable Quality Control seismic profile or seismic data cube;Selected seismic profile or seismic data cube
Need the special geobody for influencing seismic data lineups stability far from tomography, stratal pinch-outs, lithology mutation etc. as far as possible.
Step 103: in step 102 select Quality Control seismic profile or seismic data cube within the scope of certain angle into
Row phasing, and record the Quality Control seismic profile of corresponding correction angle or Quality Control seismic data cube correspond to step 1 select it is same
To the equalized amplitude energy of axis.Selected angle correct range is typically selected to -90 degree to+90 degree.
The step includes:
Step a: record Quality Control seismic profile or Quality Control seismic data cube each seismic channel correspond to step 101 select it is same
To the peak swing energy of axis;
Step b: the peak swing of all seismic channels can be measured into arithmetic mean or take root mean square as the equal of the lineups
Weigh amplitude energy.
Step 104: for the equalized amplitude energy recorded in step 103, selecting the corresponding angle of maximum equalized amplitude energy
It spends and corrects angle as the zero phase of the seismic data cube.
Step 105: the zero phase calculated in step 104 correction angle being applied to the seismic data cube, so that realizing should
The zero phase of seismic data corrects.
The following are apply a specific embodiment of the invention:
Embodiment 1
Forward modeling is carried out to the combination of different type reservoir with Ricker wavelet, from attached drawing 2 it can be seen that when seismic data is
When zero phase, seismic reflection lineups energy is most concentrated, and Energy maximum value is in reflecting interface position consistency.This is ceiling capacity
Theoretical basis has been established in the correction of method zero phase.
Processing step (as shown in Fig. 1) carries out zero phase to the practical work area seismic data of Shengli Oil Field according to the present invention
Correction.
Step 1: selecting practical work area seismic data t6 lineups to refer to as phasing based on lineups Criterion of Selecting
Lineups, the lineups have apparent Seismogeological Significance (as shown in Figure 3) as the lower oil shale top interface of sand three;
Step 2: the Quality Control seismic profile (as shown in Figure 3) for selecting the lineups far from tomography stable;
Step 3: carrying out phase school within the scope of -90 degree to an angle of 90 degrees degree for the Quality Control seismic profile selected in step 2
Just, and the equalized amplitude energy for corresponding to the Quality Control seismic profile t6 axis in the same direction of correction angle is recorded;
Step 4: for the equalized amplitude energy recorded in step 3, selected maximum equalized amplitude energy pair (as shown in Figure 4)
15 degree of the angle the answered zero phases as the seismic data cube correct angle;
Step 5: the zero phase calculated in step 4 correction angle being applied to the seismic data cube, to realize the earthquake
The zero phase of data corrects, and seismic data resolution is significantly improved (as shown in Figure 5), while zero phase after zero phase correction
Change front and back peak swing energy strata slicing comparison, amplitude attribute is remarkably reinforced, and transitional region is more clear, convenient for explaining (such as
Shown in Fig. 6).
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.
Claims (6)
1. seismic data zero phase bearing calibration, which is characterized in that the seismic data zero phase bearing calibration includes:
Step 1: selecting suitable seismic data lineups as phasing with reference to lineups on research work area seismic profile;
Step 2: selecting suitable Quality Control seismic profile or seismic data cube;
Step 3: the corresponding lineups equalized amplitude energy of the different correction angles of calculating;
Step 4: for the equalized amplitude energy recorded in step 3, selecting the corresponding angle of maximum equalized amplitude energy as ground
The zero phase for shaking data volume corrects angle;
Step 5: the zero phase calculated in step 4 correction angle being applied to seismic data cube, to carry out the zero of seismic data
Phasing.
2. seismic data zero phase bearing calibration according to claim 1, which is characterized in that in step 1, seismic data
Lineups selection standard is that selection contact relation is simple, continuous-stable in transverse direction, and formation lithology or impedance characteristic become up and down
Change apparent lineups.
3. seismic data zero phase bearing calibration according to claim 1, which is characterized in that in step 2, selected ground
Shake section or seismic data cube are to be mutated these far from tomography, stratal pinch-outs, lithology to influence seismic data lineups stability
Special geobody.
4. seismic data zero phase bearing calibration according to claim 1, which is characterized in that in step 3, for step
The Quality Control seismic profile or seismic data cube selected in 2 carry out phasing within the scope of certain angle, and record corresponding correction
The Quality Control seismic profile or Quality Control seismic data cube of angle correspond to the equalized amplitude energy of the axis in the same direction of step 1 selection.
5. seismic data zero phase bearing calibration according to claim 4, which is characterized in that in step 3, selected angle
Spending correcting range is -90 degree to+90 degree.
6. seismic data zero phase bearing calibration according to claim 4, which is characterized in that step 3 includes:
Step a: record Quality Control seismic profile or Quality Control seismic data cube each seismic channel correspond to the axis in the same direction that step 1 selects
Peak swing energy;
Step b: can measure arithmetic mean for the peak swing of all seismic channels or take equilibrium vibration of the root mean square as the lineups
Width energy.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910343233.4A CN110389378B (en) | 2019-04-26 | 2019-04-26 | Seismic data zero phase correction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910343233.4A CN110389378B (en) | 2019-04-26 | 2019-04-26 | Seismic data zero phase correction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110389378A true CN110389378A (en) | 2019-10-29 |
CN110389378B CN110389378B (en) | 2021-03-26 |
Family
ID=68285269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910343233.4A Active CN110389378B (en) | 2019-04-26 | 2019-04-26 | Seismic data zero phase correction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110389378B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112925020A (en) * | 2021-01-21 | 2021-06-08 | 滨州学院 | Parametric array type shallow profile intelligent energy compensation method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110069581A1 (en) * | 2008-08-11 | 2011-03-24 | Christine E Krohn | Removal of Surface-Wave Noise In Seismic Data |
CN103217716A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Method of fast correction processing of phase of seismic data |
US20160202376A1 (en) * | 2015-01-13 | 2016-07-14 | Bp Corporation North America Inc. | Stack Ghost Suppression |
CN106842297A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Borehole restraint unstable state method for correcting phase |
-
2019
- 2019-04-26 CN CN201910343233.4A patent/CN110389378B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110069581A1 (en) * | 2008-08-11 | 2011-03-24 | Christine E Krohn | Removal of Surface-Wave Noise In Seismic Data |
CN103217716A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Method of fast correction processing of phase of seismic data |
US20160202376A1 (en) * | 2015-01-13 | 2016-07-14 | Bp Corporation North America Inc. | Stack Ghost Suppression |
CN106842297A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Borehole restraint unstable state method for correcting phase |
Non-Patent Citations (4)
Title |
---|
S.LEVY,等: "AUTOMATIC PHASE CORRECTION OF COMMON-MIDPOINT STACKED DATA", 《GEOPHYSICS》 * |
方中于: "海上油气田三维地震的整体综合研究", 《中国博士学位论文全文数据库》 * |
曹孟起,等: "统计法同态反褶积", 《石油地球物理勘探》 * |
李庆忠: "《走向精确勘探的道路》", 31 December 1993 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112925020A (en) * | 2021-01-21 | 2021-06-08 | 滨州学院 | Parametric array type shallow profile intelligent energy compensation method |
CN112925020B (en) * | 2021-01-21 | 2022-08-05 | 滨州学院 | Parametric array type shallow profile intelligent energy compensation method |
Also Published As
Publication number | Publication date |
---|---|
CN110389378B (en) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570125B (en) | A kind of method utilizing well data to improve image taking speed model accuracy | |
US6388947B1 (en) | Multi-crosswell profile 3D imaging and method | |
US6826484B2 (en) | 3D prestack time migration method | |
CA2558000C (en) | Method, media, and signals for processing seismic data to obtain a velocity field | |
US8082107B2 (en) | Methods and computer-readable medium to implement computing the propagation velocity of seismic waves | |
CN102879799B (en) | Multi-direction seismic energy gradient difference carbonate karst cave type reservoir identification method | |
CN108333628A (en) | Elastic wave least square reverse-time migration method based on regularization constraint | |
CN102944896B (en) | The modelling static correcting method of surface survey data | |
CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
CN105956238A (en) | Palaeogeomorphic restoration method and apparatus | |
DE68909461T2 (en) | Methods for seismic data processing. | |
CN105974479A (en) | Chromatographic 2D/3D anisotropic depth domain velocity modeling method for GPU (Graphics Processing Unit) spatial grids | |
CN106547020B (en) | A kind of relative amplitude preserved processing method of seismic data | |
US20040148148A1 (en) | Well planning using seismic coherence | |
NO318304B1 (en) | Procedure for three-dimensional inversion of seismic data | |
CN110389378A (en) | Seismic data zero phase bearing calibration | |
Wang et al. | Least-squares reverse time migration with inaccurate velocity model | |
Perrone et al. | Linearized wave-equation migration velocity analysis by image warping | |
CN104597496B (en) | A kind of three dimensions method for homing of 2-d seismic data medium velocity | |
Verwest | Seismic migration in elliptically anisotropic media | |
CN108957554A (en) | Seismic inversion method in a kind of geophysical exploration | |
CA2815211A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
CN102466818A (en) | Method for imaging anisotropic medium by utilizing cross-well seismic data | |
WO2015199757A1 (en) | System and method for seismic imaging of a complex subsurface | |
CN111538080B (en) | Method of seismic imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |