CN112305613A - Static correction method and device for converted transverse wave detection point - Google Patents
Static correction method and device for converted transverse wave detection point Download PDFInfo
- Publication number
- CN112305613A CN112305613A CN201910676852.5A CN201910676852A CN112305613A CN 112305613 A CN112305613 A CN 112305613A CN 201910676852 A CN201910676852 A CN 201910676852A CN 112305613 A CN112305613 A CN 112305613A
- Authority
- CN
- China
- Prior art keywords
- wave
- static correction
- converted
- round
- superposition
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 249
- 238000012937 correction Methods 0.000 title claims abstract description 244
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000001514 detection method Methods 0.000 title abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 146
- 239000000523 sample Substances 0.000 claims abstract description 76
- 238000006243 chemical reaction Methods 0.000 claims abstract description 26
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 238000004590 computer program Methods 0.000 claims description 11
- 239000010410 layer Substances 0.000 description 45
- 238000004364 calculation method Methods 0.000 description 14
- 238000003384 imaging method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a static correction method and a static correction device for a converted transverse wave demodulator probe. The method comprises the following steps: determining an approximate converted shear wave detection point static correction value according to the average time ratio and the longitudinal wave detection point static correction value; carrying out gather extraction processing on the converted transverse wave data to obtain common asymptote gather data; performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common-center gather data; processing the common asymptote gather data to obtain a first round-robin transverse wave superposition profile; picking up a mark layer position on the longitudinal wave superposition section and the first cycle-converted transverse wave superposition section, and stretching the longitudinal wave superposition section; taking the stretched longitudinal wave superposition section as a first model channel to obtain a second-wheel converted transverse wave superposition section; repeating the processing to obtain a third converted shear wave superposition section and a residual static correction value; and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
Description
Technical Field
The invention relates to the technical field of converted transverse wave seismic data processing, in particular to a converted transverse wave demodulator probe static correction method and device.
Background
In the processing of converted transverse wave seismic data, the most key and basic processing link is the calculation of the static correction value of a converted transverse wave demodulator probe. According to the dynamic characteristics of the converted transverse wave, the converted transverse wave does not propagate in the fluid but only propagates in the rock skeleton, so that the propagation speed of the converted transverse wave is mainly determined by the rock skeleton. The longitudinal wave static correction value calculated by taking the shallow water surface as the top boundary of the high-speed layer is naturally not suitable for converting the transverse wave. In addition, the transverse wave velocity is smaller than the longitudinal wave velocity, the underground travel time is different from the longitudinal wave velocity, and the absorption attenuation influence of the earth surface on the transverse wave is larger than that of the longitudinal wave, so that the converted wave static correction value is larger, and the calculation is more difficult.
The existing international mainstream technical method mainly avoids a method for calculating static correction by utilizing cannon first arrival pickup similar to longitudinal waves, and mainly comprises the following steps: the method comprises the following steps: completing the calculation of the longitudinal wave static correction value; step two: carrying out common-detection-wave point superposition on the longitudinal wave data subjected to static correction calculation to form a longitudinal wave common-detection-wave point superposition section; step three: on the basis of applying all longitudinal wave demodulator probe static correction values, carrying out converted transverse wave common detector point superposition to form a converted transverse wave common detector point superposition section; step four: picking up geological horizons corresponding to each other on the superposed profiles of the common wave detection points of the longitudinal waves and the converted transverse waves respectively; step five: and finally, performing mathematical calculation on the picked horizon to achieve the purpose of solving the static correction value of the converted wave detection point.
However, the method is applicable to the precondition that the geological structure in the work area is simple and smooth, the signal-to-noise ratio of the converted transverse wave data is high, and the in-phase axis of the marker horizon is continuously traceable in the common detection point superposition section. Otherwise, the horizon picking difficulty and picking error are easy to occur in the human-computer interaction horizon picking process, so that a large error occurs in the process of calculating the static correction value. In addition, under the conditions of low signal-to-noise ratio of converted transverse waves and poor imaging quality, the structure cannot be controlled, and then the situation that the converted transverse waves are not accordant with the actual geological situation and cannot meet the data processing requirement occurs. Thus, the application range of the prior art methods is limited.
Disclosure of Invention
In order to solve the problems of limited application range, large error and the like in the existing converted transverse wave seismic data processing technology, the embodiment of the invention provides a converted transverse wave demodulator probe static correction method, which comprises the following steps:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
Optionally, in an embodiment of the present invention, the method further includes: and determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
Optionally, in an embodiment of the present invention, the method further includes: and determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
Optionally, in an embodiment of the present invention, before the performing multiple rounds of surface-consistent residual static correction processing and superposition processing on the data of the longitudinal wave common midpoint gather to obtain a longitudinal wave superposition profile, the method includes: and carrying out conventional denoising, amplitude compensation, deconvolution and trace set arrangement and selection on the longitudinal wave data to which the longitudinal wave demodulator point static correction value is applied to obtain longitudinal wave common center point trace set data.
The embodiment of the invention also provides a static correction device for the converted shear wave detection point, which comprises:
the approximate static correction value module is used for determining an approximate converted transverse wave demodulator probe static correction value according to the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave and the longitudinal wave demodulator probe static correction value;
a gather extraction processing module, configured to perform conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator node static correction value is applied, so as to obtain converted transverse wave common asymptote gather data;
the longitudinal wave superposition section module is used for carrying out multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common-center-point gather data to obtain a longitudinal wave superposition section;
the first round of superimposed section module is used for carrying out multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote gather data to obtain a first round of converted shear wave superimposed section;
the profile stretching module is used for picking up a mark layer on the longitudinal wave superposition profile and the first wheel-converted transverse wave superposition profile respectively and stretching the longitudinal wave superposition profile by using the mark layer as a control point;
the second-round superposed section module is used for taking the stretched longitudinal wave superposed section as a first model channel, and carrying out multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second-round converted transverse wave superposed section;
the third round of superimposed section module is used for taking the second round of converted shear wave superimposed section as a second model channel, and performing multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote channel set data to obtain a third round of converted shear wave superimposed section and a third round of residual static correction value;
and the static correction value determining module is used for superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
Optionally, in an embodiment of the present invention, the apparatus further includes: and the average time ratio module is used for determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
Optionally, in an embodiment of the present invention, the apparatus further includes: and the longitudinal wave static correction value module is used for determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
Optionally, in an embodiment of the present invention, the apparatus further includes: and the channel set arranging and selecting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and channel set arranging and selecting processing on the longitudinal wave data to which the longitudinal wave demodulator node static correction value is applied to obtain longitudinal wave common center point channel set data.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the following steps are implemented:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
The invention effectively utilizes the longitudinal wave data to control the structural characteristics of the converted transverse wave, and utilizes the self data of the converted transverse wave to carry out the residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer position error, larger error and the like caused by the inaccuracy of picking up the layer position in the low signal-to-noise ratio area. Therefore, the method has the advantages of obtaining a better effect, better solving the problem of converted wave static correction of an inclined structure and extremely poor signal-to-noise ratio, improving the converted wave imaging quality and precision and meeting the requirements of actual production.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without inventive labor.
FIG. 1 is a flow chart of a converted shear wave demodulator probe statics method according to an embodiment of the present invention;
FIG. 2 is a flow chart of a method for converted shear wave demodulator probe statics in accordance with an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a converted shear wave demodulator probe static correction device according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a static correction method and device for a converted transverse wave detection point.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of a converted shear wave demodulator probe statics correction method according to an embodiment of the present invention, where the method includes:
step S1, determining approximate converted transverse wave demodulator probe static correction value according to the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave and the longitudinal wave demodulator probe static correction value;
step S2, carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied, so as to obtain converted transverse wave common asymptote gather data;
step S3, multi-round earth surface consistency residual static correction processing and superposition processing are carried out on the longitudinal wave common center gather data to obtain a longitudinal wave superposition section;
step S4, multi-round earth surface consistency residual static correction processing and superposition processing are carried out on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition section;
step S5, picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
step S6, taking the stretched longitudinal wave superposition section as a first model road, and performing multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second wheel converted transverse wave superposition section;
step S7, taking the second-wheel converted shear wave superposition section as a second model channel, and performing multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote channel set data to obtain a third-wheel converted shear wave superposition section and a third-wheel residual static correction value;
and step S8, superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
As an embodiment of the present invention, the method further comprises: and determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
As an embodiment of the present invention, the method further comprises: and determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
In this embodiment, before performing multiple rounds of surface-consistent residual static correction processing and superposition processing on the longitudinal-wave common-center gather data to obtain a longitudinal-wave superposition profile, the method includes: and carrying out conventional denoising, amplitude compensation, deconvolution and trace set arrangement and selection on the longitudinal wave data to which the longitudinal wave demodulator point static correction value is applied to obtain longitudinal wave common center point trace set data.
In a specific embodiment of the present invention, the method specifically includes:
step 1: calculating a longitudinal wave demodulator probe static correction value;
step 2: counting the ratio gamma of the longitudinal wave velocity and the transverse wave velocity of the earth surface superficial layer of the work area;
and step 3: calculating the average ratio alpha of the corresponding positions of the longitudinal waves and the converted transverse waves of the superficial layer of the earth surface in time according to the formula 1:
wherein, TpsFor converting transverse wave horizon time, TpIs longitudinal wave horizon time;
and 4, step 4: multiplying the static positive quantity of the longitudinal wave demodulator probe by the average time ratio and applying the product to the converted shear wave single shot data;
and 5: carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted shear wave single shot data in the step 4 to obtain common asymptote gather (ACP) data;
step 6: carrying out conventional denoising, amplitude compensation, deconvolution and trace gather arrangement and selection processing on longitudinal wave data to obtain common center point trace gather (CMP) data;
and 7: performing multi-round earth surface consistency residual static correction processing on the longitudinal wave CMP data obtained in the step 6, and performing superposition processing to obtain a longitudinal wave superposition section;
and 8: performing multi-round residual static correction processing on the converted transverse wave ACP channel set in the step 5, and performing superposition processing to obtain a first round converted transverse wave superposition section;
and step 9: and picking up a mark layer on the longitudinal wave superposition section and the first wheel converted transverse wave superposition section respectively. Using the picked mark layer as a control point to stretch the longitudinal wave superposition section so as to make the longitudinal wave superposition section correspond to the converted wave superposition section in time;
step 10: and (4) taking the stretched longitudinal wave superposition section obtained in the step (9) as a model channel. Performing multi-round residual static correction calculation on the converted transverse wave ACP channel set generated in the step (4) again, and performing superposition processing to obtain a second round converted transverse wave superposition section;
step 11: and (4) taking the second-round converted shear wave superposed section obtained in the step (10) as a model channel, performing multi-round earth surface consistency residual static correction calculation again on the converted shear wave ACP channel set generated in the step (4), and performing superposition processing to obtain a third-round converted shear wave superposed section and a third-round residual static correction value.
Step 12: and accumulating the static correction amount in the step 4 and the third round of residual static correction amount in the step 11 to be used as a converted shear wave demodulator probe static correction amount.
Fig. 2 is a flowchart of a converted shear wave demodulator probe statics correction method according to an embodiment of the present invention, and particularly shows a flowchart of a low snr converted wave demodulator probe statics correction method based on configuration consistency according to an exemplary embodiment of the present invention. In the general multi-component seismic data processing process, the first completed work is the analysis work of the data quality. Whether the work area has poor earth surface receiving conditions or whether the ground abdomen of the work area has an inclined structure or not is judged through analyzing the data quality, so that the conditions that the converted transverse wave signal is weak, the signal to noise ratio is low, and accurate imaging cannot be performed on a common detection wave superposition section are caused. If not, adopting a traditional method; if so, the traditional converted shear wave static correction method is not suitable for the work area, and the method provided by the invention needs to be adopted.
According to the method, firstly, the collected longitudinal wave single-shot data of the whole work area is subjected to longitudinal wave cannon first arrival pickup. And inputting the picked first arrival information into static correction software (such as GeoEast, Omega, GeoVation, etc.). And calculating a plurality of static correction methods such as refraction static correction and chromatography static correction by using the advantages of respective software according to needs. Different types of static correction values are applied to the longitudinal wave data, and the static correction value with the best application effect is selected from the different types of static correction values and is recorded as STp. At this time, STpComprising two components, one being a detection point component STpR, one is the shot component STpAnd (4) S. For example, after applying tomographic statics correction, the longitudinal wave single shot reflection in-phase axis is smoothed over other statics correctionThe jitter phenomenon is weakened, and the superposed section has continuous in-phase axes. I.e. selecting the amount of the chromatographic static correction. And vice versa.
After the compressional wave static correction value is determined, the ratio gamma of the compressional wave speed and the shear wave speed of the superficial layer of the earth surface of the work area is counted by working according to geology and earth surface investigation, such as old data analysis, field surface investigation and the like. Since the longitudinal wave velocity is greater than the transverse wave velocity and the two have a certain correlation, γ is generally greater than 1 and ranges from about 1 to 3.
In seismic data, the record of seismic waves is plotted against time, and thus the static correction is also a time value. Which represents the amount of time each trace of seismic data moves up or down. In the case where the propagation paths of the two waves are the same, the difference in velocity inevitably differs in propagation time. Therefore, the average ratio α of the corresponding horizon of the longitudinal waves and the converted transverse waves of the superficial earth surface over time is calculated by using the formula 1. For example, if the shallow compressional-shear velocity ratio of the work area is 2, the time ratio of the longitudinal and shear wave cross sections of the shallow reflecting layer, especially the emitting layer near the earth surface, is 1.5. If a reflecting layer X is about 1 second on the longitudinal wave section, it is about 1.5 seconds on the converted shear wave section.
Wherein, TpsFor converting transverse wave horizon time, TpIs longitudinal wave horizon time;
after obtaining the time ratio, using the selected compressional static correction STpThe detected wave point component ST inpMultiplying R by the time ratio alpha to obtain an approximate converted shear wave demodulator probe static correction value STpsR, namely: ST (ST)ps_R＝α×STpR; while static correction ST of transverse wave shot components is convertedpsS still adopts the static correction value ST of the longitudinal wave shot point componentpS, i.e. STps_S＝STpAnd (4) S. During actual processing, especially in the case of inclined structure or low data signal-to-noise ratio, superposition imaging is carried out by using the approximate static correction valueAnd the phenomena of inconsistent structure and longitudinal waves and dislocation of the horizon and the phase are easy to occur, and the actual geological condition is not met.
Thus, according to the flow, the optimum static correction amount ST is appliedpThe longitudinal wave data is subjected to conventional noise attenuation, earth surface consistency amplitude compensation, deconvolution, velocity analysis and gather sorting to form a common center point (CMP) gather, and multiple rounds of earth surface consistency residual static correction P _ Res _ st, post-stack denoising and other processing are carried out to obtain an optimal longitudinal wave superposition section STK1_ P.
Similarly, the approximate static correction amount ST multiplied by the time ratiopsR is applied to converted shear wave data, and conventional noise attenuation, earth surface consistency amplitude compensation, deconvolution, speed analysis and trace gather sorting are carried out to form a common Asymptote (ACP) trace gather. And then, on the basis of the data of the converted transverse wave, carrying out N-round residual static correction processing according to the convergence condition of the earth surface consistency residual static correction iteration to obtain a residual static correction total amount PS _ Res _ st1, completing post-stack denoising, and improving the signal-to-noise ratio so as to obtain a superposition section STK1_ PS of the first-round converted transverse wave.
Picking up several marker horizons from shallow to deep on the optimal compressional and first-round converted shear stacking profiles. For example 2 or 3, according to the actual requirement of the data. From the above, the geologic horizon time in the converted shear waves is greater than the corresponding horizon time in the longitudinal waves. And stretching the longitudinal wave superposition section by using the picked mark layer as a control point, and then performing interlayer detail matching on the basis of matching the mark layer, namely completing the longitudinal and transverse wave time matching processing to ensure that the two correspond in time.
And then, taking the stretched longitudinal wave superposed section as a model channel, and performing multi-round earth surface consistency residual static correction calculation on the converted transverse wave ACP channel set again to obtain a residual static correction total amount PS _ Res _ st 2. And performing superposition and post-superposition denoising to obtain a second round-conversion transverse wave superposition section STK2_ PS. The purpose of this idea is to guarantee the requirement of consistency with the longitudinal wave structure.
And taking the obtained second-round converted shear wave superposition section STK2_ PS as a model channel, and carrying out multi-round earth surface consistency residual static correction calculation again on the converted shear wave ACP channel set to obtain a residual static correction total amount PS _ Res _ st 3. And performing superposition processing to obtain a third converted transverse wave superposition section STK3_ PS.
When the above steps are completed, the approximate converted shear wave detection point statics correction amount ST is adjustedpsR and the third remaining static correction total amount PS _ Res _ ST3 are added as a converted shear wave detection point static correction amount STpsFinal, as shown in formula 2. The calculation of the static correction of the converted shear wave demodulator probe is completed and output.
STps_Final＝STpsR + PS _ Res _ st3 equation 2
In conclusion, the invention effectively utilizes the longitudinal wave data to control the structural characteristics of the converted transverse wave, and improves the signal-to-noise ratio and the imaging effect of the data by multiple rounds of residual static correction iterative processing of the earth surface consistency on the basis of the self data of the converted transverse wave.
The invention effectively utilizes the longitudinal wave data to control the structural characteristics of the converted transverse wave, and utilizes the self data of the converted transverse wave to carry out the residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer position error, larger error and the like caused by the inaccuracy of picking up the layer position in the low signal-to-noise ratio area. Therefore, the method has the advantages of obtaining a better effect, better solving the problem of converted wave static correction of an inclined structure and extremely poor signal-to-noise ratio, improving the converted wave imaging quality and precision and meeting the requirements of actual production.
Fig. 3 is a schematic structural diagram of a converted shear wave demodulator probe statics correction apparatus according to an embodiment of the present invention, where the apparatus includes:
an approximate static correction value module 10, configured to determine an approximate converted shear wave demodulator probe static correction value according to an average time ratio of the corresponding layers of the longitudinal waves and the converted shear waves and a longitudinal wave demodulator probe static correction value;
a gather extraction processing module 20, configured to perform conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator node static correction value is applied, so as to obtain converted transverse wave common asymptote gather data;
the longitudinal wave superposition section module 30 is used for performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common-center-point gather data to obtain a longitudinal wave superposition section;
a first round of superimposed section module 40, configured to perform multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote gather data to obtain a first round of converted shear wave superimposed section;
a profile stretching module 50, configured to pick up a mark level on the longitudinal wave superposition profile and the first round-conversion transverse wave superposition profile, respectively, and stretch the longitudinal wave superposition profile by using the mark level as a control point;
a second-round superimposed section module 60, configured to use the stretched longitudinal wave superimposed section as a first model road, and perform multi-round earth-surface consistency residual static correction processing and superimposing processing on the converted shear wave common asymptote gather data to obtain a second-round converted shear wave superimposed section;
a third round of superimposed section module 70, configured to use the second round of converted shear wave superimposed section as a second model road, and perform multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote gather data to obtain a third round of converted shear wave superimposed section and a third round of residual static correction value;
and a static correction value determining module 80, configured to superimpose the approximate converted transverse wave detection point static correction value and the third remaining static correction value to obtain a final converted transverse wave detection point static correction total amount.
As an embodiment of the present invention, the apparatus further comprises: and the average time ratio module is used for determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
As an embodiment of the present invention, the apparatus further comprises: and the longitudinal wave static correction value module is used for determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
In this embodiment, the apparatus further includes: and the channel set arranging and selecting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and channel set arranging and selecting processing on the longitudinal wave data to which the longitudinal wave demodulator node static correction value is applied to obtain longitudinal wave common center point channel set data.
Based on the same application concept as the converted transverse wave demodulator probe static correction method, the invention also provides the converted transverse wave demodulator probe static correction device. Because the principle of solving the problems of the converted transverse wave detection point static correction device is similar to that of the converted transverse wave detection point static correction method, the implementation of the converted transverse wave detection point static correction device can refer to the implementation of the converted transverse wave detection point static correction method, and repeated parts are not repeated.
The invention effectively utilizes the longitudinal wave data to control the structural characteristics of the converted transverse wave, and utilizes the self data of the converted transverse wave to carry out the residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer position error, larger error and the like caused by the inaccuracy of picking up the layer position in the low signal-to-noise ratio area. Therefore, the method has the advantages of obtaining a better effect, better solving the problem of converted wave static correction of an inclined structure and extremely poor signal-to-noise ratio, improving the converted wave imaging quality and precision and meeting the requirements of actual production.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the following steps are implemented:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
The invention also provides the computer equipment and a computer readable storage medium based on the same application concept as the converted shear wave demodulator probe static correction method. Because the principle of solving the problems of the computer equipment and the computer-readable storage medium is similar to that of a converted transverse wave detection point static correction method, the implementation of the computer equipment and the computer-readable storage medium can refer to the implementation of the converted transverse wave detection point static correction method, and repeated parts are not repeated.
The invention effectively utilizes the longitudinal wave data to control the structural characteristics of the converted transverse wave, and utilizes the self data of the converted transverse wave to carry out the residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer position error, larger error and the like caused by the inaccuracy of picking up the layer position in the low signal-to-noise ratio area. Therefore, the method has the advantages of obtaining a better effect, better solving the problem of converted wave static correction of an inclined structure and extremely poor signal-to-noise ratio, improving the converted wave imaging quality and precision and meeting the requirements of actual production.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by relevant hardware instructed by a program, and the program may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method of converted shear wave demodulator probe statics correction, the method comprising:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
2. The method of claim 1, further comprising: and determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
3. The method of claim 1, further comprising: and determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
4. The method of claim 3, wherein before performing multiple rounds of surface-consistent residual static correction and superposition on the compressional common-midpoint gather data to obtain compressional superposition profiles, the method comprises: and carrying out conventional denoising, amplitude compensation, deconvolution and trace set arrangement and selection on the longitudinal wave data to which the longitudinal wave demodulator point static correction value is applied to obtain longitudinal wave common center point trace set data.
5. A converted shear wave demodulator probe statics correction apparatus, comprising:
the approximate static correction value module is used for determining an approximate converted transverse wave demodulator probe static correction value according to the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave and the longitudinal wave demodulator probe static correction value;
a gather extraction processing module, configured to perform conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator node static correction value is applied, so as to obtain converted transverse wave common asymptote gather data;
the longitudinal wave superposition section module is used for carrying out multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common-center-point gather data to obtain a longitudinal wave superposition section;
the first round of superimposed section module is used for carrying out multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote gather data to obtain a first round of converted shear wave superimposed section;
the profile stretching module is used for picking up a mark layer on the longitudinal wave superposition profile and the first wheel-converted transverse wave superposition profile respectively and stretching the longitudinal wave superposition profile by using the mark layer as a control point;
the second-round superposed section module is used for taking the stretched longitudinal wave superposed section as a first model channel, and carrying out multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second-round converted transverse wave superposed section;
the third round of superimposed section module is used for taking the second round of converted shear wave superimposed section as a second model channel, and performing multi-round earth surface consistency residual static correction processing and superimposed processing on the converted shear wave common asymptote channel set data to obtain a third round of converted shear wave superimposed section and a third round of residual static correction value;
and the static correction value determining module is used for superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
6. The apparatus of claim 5, further comprising: and the average time ratio module is used for determining the average time ratio of the longitudinal wave to the corresponding layer of the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
7. The apparatus of claim 5, further comprising: and the longitudinal wave static correction value module is used for determining a longitudinal wave demodulator probe static correction value according to the longitudinal wave data.
8. The apparatus of claim 7, further comprising: and the channel set arranging and selecting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and channel set arranging and selecting processing on the longitudinal wave data to which the longitudinal wave demodulator node static correction value is applied to obtain longitudinal wave common center point channel set data.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the computer program implements the steps of:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of:
determining approximate converted transverse wave demodulator probe static correction values according to the average time ratio of the longitudinal waves to the corresponding layer positions of the converted transverse waves and the longitudinal wave demodulator probe static correction values;
carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave demodulator point static correction value is applied to obtain converted transverse wave common asymptote gather data;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common midpoint gather data to obtain a longitudinal wave superposition section;
performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a first round converted shear wave superposition profile;
respectively picking up a mark layer on the longitudinal wave superposition section and the first wheel-converted transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposed section as a first model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second round converted transverse wave superposed section;
taking the second round converted shear wave superposition section as a second model road, and performing multi-round earth surface consistency residual static correction processing and superposition processing on the converted shear wave common asymptote gather data to obtain a third round converted shear wave superposition section and a third round residual static correction value;
and superposing the approximate conversion transverse wave demodulator probe static correction value and the third round of residual static correction value to obtain the final conversion transverse wave demodulator probe static correction total amount.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910676852.5A CN112305613B (en) | 2019-07-25 | 2019-07-25 | Static correction method and device for converted transverse wave detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910676852.5A CN112305613B (en) | 2019-07-25 | 2019-07-25 | Static correction method and device for converted transverse wave detector |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112305613A true CN112305613A (en) | 2021-02-02 |
| CN112305613B CN112305613B (en) | 2024-03-01 |
Family
ID=74330104
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910676852.5A Active CN112305613B (en) | 2019-07-25 | 2019-07-25 | Static correction method and device for converted transverse wave detector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112305613B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3539982A (en) * | 1969-03-03 | 1970-11-10 | Texas Instruments Inc | Method and apparatus for consistent static correction of common depth point seismic signals |
| US4881209A (en) * | 1988-10-12 | 1989-11-14 | Mobil Oil Corporation | Method of processing converted mode seismic data |
| US20130176823A1 (en) * | 2012-01-05 | 2013-07-11 | Cggveritas Services Sa | Simultaneous joint estimation of the p-p and p-s residual statics |
| CN103630934A (en) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | Method for determining large transverse wave static correction value of converted wave demodulator probe |
| CN104155694A (en) * | 2014-08-28 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Residual static correction method of reflection converted shear wave common geophone stack section |
| CN104570110A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Multi-component data joint speed analysis method based on longitudinal and horizontal wave matching |
| CN104914469A (en) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | Static correction method and device for converted S-wave |
| CN106443791A (en) * | 2016-11-11 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for calculating transverse wave residual static correction quantity of tilted stratum or anisotropic formation stratum |
| US20170176617A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
-
2019
- 2019-07-25 CN CN201910676852.5A patent/CN112305613B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3539982A (en) * | 1969-03-03 | 1970-11-10 | Texas Instruments Inc | Method and apparatus for consistent static correction of common depth point seismic signals |
| US4881209A (en) * | 1988-10-12 | 1989-11-14 | Mobil Oil Corporation | Method of processing converted mode seismic data |
| US20130176823A1 (en) * | 2012-01-05 | 2013-07-11 | Cggveritas Services Sa | Simultaneous joint estimation of the p-p and p-s residual statics |
| CN103630934A (en) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | Method for determining large transverse wave static correction value of converted wave demodulator probe |
| CN104570110A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Multi-component data joint speed analysis method based on longitudinal and horizontal wave matching |
| CN104914469A (en) * | 2014-03-13 | 2015-09-16 | 中国石油化工股份有限公司 | Static correction method and device for converted S-wave |
| CN104155694A (en) * | 2014-08-28 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Residual static correction method of reflection converted shear wave common geophone stack section |
| US20170176617A1 (en) * | 2015-12-18 | 2017-06-22 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
| CN106443791A (en) * | 2016-11-11 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for calculating transverse wave residual static correction quantity of tilted stratum or anisotropic formation stratum |
Non-Patent Citations (4)
| Title |
|---|
| PU WANG等: "Converted-wave imaging technology and application for complex structures" * |
| 杨海涛;贺振华;陈丹;王栋;王静;王鸿燕;: "转换波静校正技术在四川地区的应用研究", no. 04 * |
| 王淑荣等: "转换波处理方法研究" * |
| 陈海峰等: "联合纵波确定转换波接收点大的静校正量" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112305613B (en) | 2024-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107505651B (en) | Seismic first break and back wave combine slope chromatography imaging method | |
| US9013956B2 (en) | Method and system for seismic imaging and earth modeling using beam tomography | |
| US10557956B2 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
| CN104237940B (en) | A kind of diffraction wave imaging method based on dynamic characteristic and device | |
| CN109917454B (en) | True earth surface prestack depth migration imaging method and device based on double reference surfaces | |
| CN102590862B (en) | Prestack time migration method for compensating absorptive attenuation | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| CN107219554B (en) | The automatic obtaining method of the Value of residual static correction of land seismic data | |
| CN105182408A (en) | Manufacturing method and device for synthesizing earthquake record | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| CN104199103B (en) | Static correction method and static correction device for converted wave | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| US20160161619A1 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| CN110907995B (en) | Reverse time migration method and device for VSP seismic data in well | |
| CN113740901A (en) | Land seismic data full-waveform inversion method and apparatus based on complex undulating surface | |
| CN104977615B (en) | A kind of multiple ripple drawing method of deep water OBC data based on modeling statistics pickup | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN109655890B (en) | Depth domain shallow-medium-deep layer combined chromatography inversion speed modeling method and system | |
| Duret et al. | Near-surface velocity modeling using a combined inversion of surface and refracted P-waves | |
| CN103645499A (en) | Earth surface consistency amplitude compensation method based on overlaid reflection wave energy statistics | |
| CN104597485B (en) | A kind of miniature fault detection method and fault detection device | |
| CN102798888B (en) | Method for calculating velocity ratio of longitudinal wave to transverse wave by using non-zero wellhead distance data | |
| CN104570090B (en) | The extraction of full waveform inversion noise filter operator and the method filtered using its noise | |
| CN104155688A (en) | High precision weighted stack method | |
| CN105445795B (en) | The velocity analysis method and its system of tilted interface |
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 |