CN1325937C - Static correcting method during converted wave delaying - Google Patents
Static correcting method during converted wave delaying Download PDFInfo
- Publication number
- CN1325937C CN1325937C CNB2004100885637A CN200410088563A CN1325937C CN 1325937 C CN1325937 C CN 1325937C CN B2004100885637 A CNB2004100885637 A CN B2004100885637A CN 200410088563 A CN200410088563 A CN 200410088563A CN 1325937 C CN1325937 C CN 1325937C
- Authority
- CN
- China
- Prior art keywords
- wave
- shear
- converted
- compressional
- velocity
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000012937 correction Methods 0.000 claims abstract description 27
- 230000005540 biological transmission Effects 0.000 claims description 5
- 230000000644 propagated effect Effects 0.000 claims description 4
- 238000011835 investigation Methods 0.000 claims description 3
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 230000005570 vertical transmission Effects 0.000 claims description 2
- 239000000523 sample Substances 0.000 abstract 3
- 230000005284 excitation Effects 0.000 abstract 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Abstract
The present invention relates to a converted wave static correction method which is a set of complete converted wave static correction method. The present invention receives by adopting longitudinal wave excitation three-components, demodulator probes corresponding to a Z component and an X component are mutually interrelated, mutual correlation results of the demodulator probes are superimposed to obtain a mutual correlation stacked section, the first arrival difference of the two components are determined according to the peak position of the section, and the converted wave delay time is calculated by the first arrival difference of the two components; then, a transverse wave mould close to the ground is established according to the converted wave delay time, the accurate transverse wave datum plane static correction amount of the demodulator probes is resolved finally, and the obtained static correction amount is a result for correcting a datum plane.
Description
Technical field
The present invention relates to the geophysical prospecting for oil data processing technique, is a kind of when utilizing up transmission transformed wave and compressional wave to ask for converted wave delaying in the first arrival time difference that waste mantle causes, sets up the shear wave near-surface model, calculates the method for shear wave statics.
Background technology
Transformed wave seismic prospecting is a kind of special method of exploration.It excites compressional wave by blast or ground vibroseis in the well, and compressional wave is propagated downwards in the stratum, and behind the arrival interface, part energy is converted to shear wave (SV ripple and SH ripple).Compressional wave, SV ripple and three components of SH ripple are upwards propagated respectively and are arrived the face of land.The face of land receives with three-component seismometer, to obtain complete seismic wave field information.
The influence different with compressional wave (P ripple), that shear wave (S ripple) is not risen and fallen by the near surface water table.This mainly is because the sediment of saturated with fluid has bigger breathing power than dried deposit, modulus of shearing be subjected to water saturated influence very little.After up S ripple entered water table aquifer, velocity variations was very little, continued to keep low speed; And the low velocity layer (LVL) of up P ripple ends at water table.Therefore, the near-surface low-velocity layer variable thickness of P ripple and S ripple causes, and makes that almost it doesn't matter between the static correction value of compressional wave and shear wave.Normally incorrect with compressional wave static correction value linear-apporximation shear wave statics.Shear wave statics on the same position usually than the big 2-10 of compressional wave doubly.
Existing " utilizing up transmission transformed wave to calculate the method for shear wave static correction " patent of invention supposes that up compressional wave and shear wave are nearly vertical-path propagation in low velocity layer (LVL), therefore utilize the time difference of transformed wave first arrival and compressional wave first arrival, add the static correction value of compressional wave, as the static correction value of transformed wave, obtained effect preferably.But in fact, there are certain angle in upward traveling wave and vertical direction, and the time difference is different during with converted wave delaying, directly utilizes relative time error as relative static correction value, still has certain static correction value problem.
In addition, many static correction operators are found the solution shot point and receiver static correction amount simultaneously.But for the P-SV ripple, a demand is separated the static correction value of geophone station.This just requires to suppress shot statics and finds the solution.Therefore, must adopt static correcting method targetedly.
Summary of the invention
The present invention utilizes the first arrival time difference of transformed wave and compressional wave, the notion when proposing converted wave delaying.Set up the shear wave surface-level model during then according to converted wave delaying, thereby calculate the datum static correction amount of shear wave.
At the problem that exists in the background technology, the present invention is according to the relation of compressional wave and transformed wave raypath, the relation when deriving converted wave delaying, when compressional wave and the transformed wave first arrival time difference and compressional wave postpone.According to this relation, when at first asking for converted wave delaying, set up the shear wave near-surface model again, finally ask for shear wave datum static correction amount.Concrete steps are as follows:
When 1, calculating converted wave delaying (Fig. 1)
For double-layer structure, velocity of longitudinal wave is respectively VP1 and VP2, and shear wave velocity is V
S1, A and A ' they are shot point, D is an acceptance point.Excited by the A point, the D point receives, and compressional wave refraction wave path is ABCD, and α is critical angle, and the velocity of propagation of AB section and CD section is VP1, and BC section velocity of propagation is VP2; Transformed wave slides with VP2 speed along the interface, arrives the E point and produces transformed wave, arrives the D point along ray ED with speed VS1.
Theoretical according to refraction:
Because VP1＞VS1, so α＞β, the more close D point of transfer point E C point is as Fig. 1.
From refraction path as can be seen, for any one shot point A or A ', compressional wave is identical to the mistiming Δ t of acceptance point D with transformed wave.Can obtain according to raypath:
Δt＝[T
AB+T
BE+T
ED]-[T
AB+T
BC+T
CD]
＝[T
A′B′+T
B′E+T
ED]-[T
A′B′+T
B′C+T
CD] (1)
＝T
CE+T
ED-T
CD
Transformed wave and compressional wave time difference Δ t can be converted into following form:
Δt＝T
CE+T
ED-T
CD＝T
ED-T
EF+T
EF+T
CE-T
CD
(2)
＝(T
ED-T
EF)-(T
CD-T
CF)
Order: d
P=T
CD-T
CF, d
S=T
ED-T
EF
When obtaining converted wave delaying thus two parts are formed: when the first arrival time difference of transformed wave and compressional wave and compressional wave postpone, both when being converted wave delaying:
d
S＝Δt+d
P (3)
Wherein, d
SWhen being defined as converted wave delaying, d
pWhen postponing for compressional wave, Δ t is the first arrival time difference of transformed wave and compressional wave.
2, set up the shear wave near-surface model, determine thickness h or calculation of near surface shear velocity VS1
Can obtain calculation of near surface shear velocity according to near-surface investigation, obtain the calculation of near surface shear velocity of the whole district after the interpolation; Perhaps, obtain calculation of near surface shear velocity indirectly according to velocity of longitudinal wave aspect ratio parameter.During according to converted wave delaying, calculation of near surface shear velocity and compressional wave refraction wave speed, can calculate the thickness on top layer, formula is as follows:
Otherwise, because compressional wave and transformed wave occur in same interface, therefore can be from the thickness of compressional wave data computation shear wave near-surface model, the calculation of near surface shear velocity of calculation of parameter near-surface model such as during then according to converted wave delaying, formula is as follows:
3, calculate shear wave statics
According to near-surface model, calculate shear wave datum static correction amount.Computing formula is as follows:
Wherein, h
iBe the thickness of each layer of near-surface model, V
iBe the shear wave velocity of each layer, H
dBe datum elevation, H
gBe bottommost layer interface elevation, V
sBe the shear wave replacement velocity, τ is the vertical transmission time of shot point seismic event from the shaft bottom to ground.
The invention has the advantages that:
The present invention is a kind of complete converted wave statics method of transformed wave seismic prospecting.When it postpones according to compressional wave and first arrival time difference when calculating converted wave delaying, set up the shear wave near-surface model during according to converted wave delaying, finally calculate shear wave statics.The static correction value that this method provides is the precise results that is corrected to reference field.
Fig. 2 is the stacked section contrast before and after the static correction of transformed wave geological data.The shear wave static correction problem is better solved as can be seen, and the lineups continuity obviously improves, and it is smooth and continuous that the lineups of distortion become, and especially for the relatively more outstanding left side of static correction problem, effect is remarkable.
Description of drawings
Fig. 1 is compressional wave refraction wave and transformed wave refraction wave raypath synoptic diagram
Fig. 2 is before the static correction and static correction post-detection point stacked section comparison diagram
Fig. 3 is Z component and X component simple crosscorrelation stacked section
Fig. 4 is the shear wave near-surface model
Embodiment
Embodiment 1
Because the transformed wave and the compressional wave first arrival time difference are constant on the common detector gather, (first arrival of X component is the energy of up transmission converted shear wave by geophone station X component and Z component, the Z component first arrival is the energy of up compressional wave refraction) in giving the timing window scope, carry out simple crosscorrelation, and geophone station cross correlation results stack obtained the simple crosscorrelation stacked section, determine the first arrival time difference Δ t (Fig. 3) of transformed wave and compressional wave according to energy peak.D when postponing with corresponding geophone station compressional wave by Δ t
pD when sum is calculated the converted wave delaying of this geophone station
s
During according to the known calculation of near surface shear velocity of shear wave near-surface investigation, calculate weathering thickness, set up shear wave near-surface model (Fig. 4) according to (4) formula.Calculate the benchmark of shear wave and static correction value according to (6) formula then.
Embodiment 2
Because the transformed wave and the compressional wave first arrival time difference are constant on the common detector gather, in giving the timing window scope, carry out simple crosscorrelation by geophone station X component and Z component, and geophone station cross correlation results stack obtained the simple crosscorrelation stacked section, determine the first arrival time difference Δ t of transformed wave and compressional wave according to energy peak.D when postponing with corresponding geophone station compressional wave by Δ t
pD when sum is calculated the converted wave delaying of this geophone station
s
By compressional wave data computation near-surface model thickness, calculate calculation of near surface shear velocity according to (5) formula, set up the shear wave near-surface model.Calculate the datum static correction amount of shear wave then according to (6) formula.
Claims (1)
1, a kind of static correcting method during converted wave delaying, it is to utilize up transmission transformed wave to calculate the method for shear wave statics in the seismic prospecting, excite compressional wave by blast or ground vibroseis in the well, compressional wave is propagated downwards in the stratum, after arriving the interface, part energy is converted to shear wave (SV ripple and SH ripple), compressional wave, SV ripple and three components of SH ripple are upwards propagated respectively and are arrived the face of land, the face of land receives with three-component seismometer, obtain complete seismic wave field information, it is characterized in that: (d when postponing according to first arrival time difference (Δ t) of transformed wave and compressional wave and compressional wave
P) (d when calculating converted wave delaying
s), d
S=Δ t+d
P, set up the shear wave surface-level model during then according to converted wave delaying, calculate shear wave datum static correction amount at last,
Specific practice is as follows:
At first, utilize the first arrival (X component) of the up transmission converted shear wave of common detector gather and the first arrival (Z component) of up compressional wave refraction, in giving the timing window scope, carry out simple crosscorrelation, and geophone station cross correlation results stack obtained the simple crosscorrelation stacked section, determine the first arrival time difference Δ t of transformed wave and compressional wave according to energy peak;
The second, when calculating converted wave delaying
For double-layer structure, velocity of longitudinal wave is respectively V
P1And V
P2Shear wave velocity is V
S1, A and A ' they are shot point, and D is an acceptance point, is excited by the A point, and the D point receives, and compressional wave refraction wave path is ABCD, and α is critical angle, and the velocity of propagation of AB section and CD section is V
P1, BC section velocity of propagation is V
P2Transformed wave along the interface with V
P2Speed slides, and arrives the E point and produces transformed wave, arrives the D point along ray ED with speed VS1,
Theoretical according to refraction:
Because V
P1＞V
S1So, α＞β, the more close D point of transfer point EC point,
For any one shot point A or A ', compressional wave is identical to the mistiming Δ t of acceptance point D with transformed wave, can obtain according to raypath:
Δt＝[T
AB+T
BE+T
ED]-[T
AB+T
BC+T
CD]
＝[T
A′B′+T
B′E+T
ED]-[T
A′B′+T
B′C+T
CD]
＝T
CE+T
ED-T
CD
Transformed wave and compressional wave time difference Δ t can be converted into following form:
Δt＝T
CE+T
ED-T
CD＝T
ED-T
EF+T
EF+T
CE-T
CD
＝(T
ED-T
EF)-(T
CD-T
CF)
Order: d
P=T
CD-T
CF, d
S=T
ED-T
EF
When obtaining converted wave delaying thus two parts are formed: when the first arrival time difference of transformed wave and compressional wave and compressional wave postpone, both when being converted wave delaying:
d
S＝Δt+d
P
Wherein, d
SWhen being defined as converted wave delaying, d
PWhen postponing for compressional wave, Δ t is the first arrival time difference of transformed wave and compressional wave;
The 3rd, set up the shear wave near-surface model, determine thickness h or calculation of near surface shear velocity V
S1
Can obtain calculation of near surface shear velocity V according to near-surface investigation
S1, obtain the calculation of near surface shear velocity of the whole district after the interpolation; Perhaps, obtain calculation of near surface shear velocity indirectly, d during according to converted wave delaying according to velocity of longitudinal wave aspect ratio parameter
S, calculation of near surface shear velocity V
S1With compressional wave refraction wave speed V
S2, calculating the thickness h on top layer, formula is as follows:
According to the thickness h of compressional wave data computation shear wave near-surface model, d during then according to converted wave delaying
SWith compressional wave refraction wave speed V
S2, the calculation of near surface shear velocity V of calculating near-surface model
S1, computing formula is as follows:
At last, calculate shear wave statics
According to near-surface model, calculate shear wave datum static correction amount, computing formula is as follows:
Wherein, h
iBe the thickness of each layer of near-surface model, V
iBe the shear wave velocity of each layer, H
dBe datum elevation, H
gBe bottommost layer interface elevation, V
sBe the shear wave replacement velocity, τ is the vertical transmission time of shot point seismic event from the shaft bottom to ground.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100885637A CN1325937C (en) | 2004-11-09 | 2004-11-09 | Static correcting method during converted wave delaying |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100885637A CN1325937C (en) | 2004-11-09 | 2004-11-09 | Static correcting method during converted wave delaying |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1773310A CN1773310A (en) | 2006-05-17 |
| CN1325937C true CN1325937C (en) | 2007-07-11 |
Family
ID=36760370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004100885637A Expired - Fee Related CN1325937C (en) | 2004-11-09 | 2004-11-09 | Static correcting method during converted wave delaying |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1325937C (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101363916B (en) * | 2007-08-10 | 2012-05-30 | 中国石油集团东方地球物理勘探有限责任公司 | High precision refraction static correction data inversion method |
| CN101393270B (en) * | 2007-09-21 | 2011-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | Continuum area surface layer high precision static correction method |
| CN101598810B (en) * | 2008-06-04 | 2011-06-22 | 中国石油天然气集团公司 | Method for processing refraction static correction border |
| CN101598811B (en) * | 2008-06-04 | 2011-04-20 | 中国石油天然气集团公司 | Method for calculating shot-static correction for two-dimensional vertical seismic section data |
| CN101609165B (en) * | 2008-06-19 | 2011-05-25 | 中国石油天然气集团公司 | Refractive wave static correction method based on relief surface |
| CN103630934B (en) * | 2012-08-23 | 2016-12-21 | 中国石油天然气集团公司 | A kind of method determining shear wave statics that converted wave geophone station is big |
| CN104133244B (en) * | 2013-05-02 | 2017-03-08 | 中国石油化工股份有限公司 | A kind of near surface static corrections processing method of focus shear wave exploration |
| CN103293554B (en) * | 2013-05-07 | 2014-12-24 | 东北石油大学 | Statics correction method based on seismic data |
| CN103412333B (en) * | 2013-08-26 | 2015-11-18 | 郑鸿明 | A kind of static correction base level defining method |
| CN103487834B (en) * | 2013-09-09 | 2016-07-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Converted wave common geophone stack static correcting method |
| CN104656139B (en) * | 2013-11-21 | 2018-04-06 | 中国石油天然气集团公司 | The static correcting method and equipment of a kind of converted wave |
| CN105093319B (en) * | 2014-05-14 | 2017-11-28 | 中国石油化工股份有限公司 | Ground micro-seismic static correcting method based on 3D seismic data |
| CN105093316B (en) * | 2014-05-14 | 2017-10-13 | 中国石油化工股份有限公司 | Declined geological data static correcting method along the surface array of well track search |
| CN104570111B (en) * | 2015-01-21 | 2016-03-02 | 中国矿业大学(北京) | Altogether attitude Dao Ji position angle is analyzed and bearing calibration and device |
| CN111965701B (en) * | 2019-05-20 | 2023-09-26 | 中国石油天然气集团有限公司 | Near-surface structure inversion method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4405999A (en) * | 1980-12-31 | 1983-09-20 | Mobil Oil Corporation | Method for collecting and generating composite trace signals with improved signal to noise ratios |
| US6442490B1 (en) * | 2000-06-30 | 2002-08-27 | Pgs Americas, Inc. | Vertical cable time processing |
| GB2387227A (en) * | 2002-04-01 | 2003-10-08 | Pgs Americas Inc | Determining statics by using pilot trace defined by CDP, azimuth and offset information |
-
2004
- 2004-11-09 CN CNB2004100885637A patent/CN1325937C/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4405999A (en) * | 1980-12-31 | 1983-09-20 | Mobil Oil Corporation | Method for collecting and generating composite trace signals with improved signal to noise ratios |
| US6442490B1 (en) * | 2000-06-30 | 2002-08-27 | Pgs Americas, Inc. | Vertical cable time processing |
| GB2387227A (en) * | 2002-04-01 | 2003-10-08 | Pgs Americas Inc | Determining statics by using pilot trace defined by CDP, azimuth and offset information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1773310A (en) | 2006-05-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1325937C (en) | Static correcting method during converted wave delaying | |
| CN102053261B (en) | Method for processing seismic data | |
| CN102967882B (en) | The modeling method of the interval velocity model on stratum | |
| CN101630016B (en) | Method for improving imaging quality of vertical seismic profile | |
| US20020075759A1 (en) | Static correction method for exploration seismic data using first arrivals of seismic waves | |
| US6188964B1 (en) | Method for using global optimization to the estimation of surface-consistent residual statics | |
| CN106226818A (en) | Seismic data processing technique and device | |
| CN108196305B (en) | Mountain land static correction method | |
| CN103630934B (en) | A kind of method determining shear wave statics that converted wave geophone station is big | |
| CN101363916A (en) | High precision refraction static correction data inversion method | |
| CN103116184A (en) | Method for calculating transverse wave statics correction value in surface survey | |
| CN102073067A (en) | Method for increasing automatic pickup efficiency of geological data first-motion waves | |
| CN102692648A (en) | Refraction wave residual static correction method based on cannon first arrival | |
| US8139441B2 (en) | Method for bispectral picking of anelliptical NMO correction parameters | |
| US20060062083A1 (en) | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration | |
| AU2004220765A1 (en) | Method for stable estimation of anisotropic parameters for P-wave prestack imaging | |
| CN109444959B (en) | Full-frequency high-precision layer velocity field establishment method | |
| CN101598808A (en) | A kind of method that improves image quality of seismic data | |
| CN104122582B (en) | The method that accurately seismic velocity is asked for using stack velocity | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| CN107817516A (en) | Near surface modeling method and system based on preliminary wave information | |
| CN102043167B (en) | Curved earth surface seismic data velocity analysis method | |
| CN104977615A (en) | Model-statistics-pickup-based multiple suppression method of deep-sea OBC data | |
| CN104914469A (en) | Static correction method and device for converted S-wave | |
| CN101598810B (en) | Method for processing refraction static correction border |
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 |
Granted publication date: 20070711 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |