CN112305613B - Static correction method and device for converted transverse wave detector - Google Patents
Static correction method and device for converted transverse wave detector Download PDFInfo
- Publication number
- CN112305613B CN112305613B CN201910676852.5A CN201910676852A CN112305613B CN 112305613 B CN112305613 B CN 112305613B CN 201910676852 A CN201910676852 A CN 201910676852A CN 112305613 B CN112305613 B CN 112305613B
- Authority
- CN
- China
- Prior art keywords
- transverse wave
- superposition
- static correction
- wave
- processing
- 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.)
- Active
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 242
- 238000012937 correction Methods 0.000 title claims abstract description 233
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 160
- 238000006243 chemical reaction Methods 0.000 claims abstract description 57
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims description 30
- 238000004590 computer program Methods 0.000 claims description 10
- 239000010410 layer Substances 0.000 description 49
- 238000004364 calculation method Methods 0.000 description 14
- 238000003384 imaging method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 239000003550 marker Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011835 investigation Methods 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
- 238000010276 construction 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
- 239000000463 material 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/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 device for a converted transverse wave detector. The method comprises the following steps: determining the static correction value of the approximate conversion transverse wave detector according to the average time ratio and the static correction value of the longitudinal wave detector; performing gather extraction processing on the converted transverse wave data to obtain common asymptote gather data; carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data; carrying out the processing on the common asymptote gather data to obtain a first rotary transverse wave superposition profile; picking up a mark layer position on the longitudinal wave superposition section and the first rotary transverse wave superposition section, and stretching the longitudinal wave superposition section; taking the stretched longitudinal wave superposition profile as a first model channel to obtain a second rotary transverse wave superposition profile; repeating the processing to obtain a third round of converted transverse wave superposition profile and a residual static correction value; and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total.
Description
Technical Field
The invention relates to the technical field of converted transverse wave seismic data processing, in particular to a method and a device for static correction of a converted transverse wave detector.
Background
In the process of converting transverse wave seismic data, the most critical and fundamental processing links are the calculation of the static correction value of a converted transverse wave detector. According to the dynamic characteristics of the converted transverse wave, the converted transverse wave does not propagate in the fluid and only propagates in the rock framework, so that the propagation speed of the converted transverse wave is mainly determined by the rock framework. The longitudinal wave static correction value calculated by taking the shallow water surface as the top boundary of the high-speed layer is not suitable for converting transverse waves naturally. In addition, the transverse wave speed is smaller than the longitudinal wave speed, the underground travel time is different from the longitudinal wave, 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 amount is larger, and the calculation is more difficult.
The existing international mainstream technical method mainly avoids the method for calculating static correction by utilizing the first arrival pickup of cannons similar to longitudinal waves, and mainly comprises the following steps: step one: finishing longitudinal wave static correction value calculation; step two: the longitudinal wave data subjected to all static correction calculation are subjected to common-wave-detection-point superposition to form a longitudinal wave common-wave-detection-point superposition section; step three: on the basis of applying all longitudinal wave detector static correction values, performing superposition of the converted transverse wave common detector to form a superposition profile of the converted transverse wave common detector; step four: respectively picking up geological horizons corresponding to each other on the superposition profile of the longitudinal wave and the converted transverse wave common detection point; step five: and finally, carrying out mathematical calculation on the picked horizon to achieve the purpose of solving the static correction value of the converted wave detector.
However, the method has the application precondition that the geological structure in the work area is simple and gentle, the signal-to-noise ratio of the converted transverse wave data is high, and the same-phase axis of the mark layer can be continuously tracked in the common-detection-point superposition section. Otherwise, the horizon picking difficulty and the picking error easily occur in the human-computer interaction horizon picking process, so that larger errors occur in the process of calculating the static correction value. In addition, under the conditions of low signal-to-noise ratio and poor imaging quality of the converted transverse wave, the structure cannot be controlled, and then the situation that the structure does not conform to the actual geological condition and cannot meet the data processing requirement appears. Therefore, the application range of the prior art method is limited.
Disclosure of Invention
In order to solve the problems of limited application range, large error and the like in the current converted transverse wave seismic data processing technology, the embodiment of the invention provides a converted transverse wave detection point static correction method, which comprises the following steps:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel 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 profile;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first rotary transverse wave superposition section;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
taking the second rotary transverse wave superposition profile as a second 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 third round of converted transverse wave superposition profile and a third round of residual static correction value;
and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total amount.
Optionally, in an embodiment of the present invention, the method further includes: and determining the average time ratio of the corresponding layers of the longitudinal wave and 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 the static correction amount of the longitudinal wave detector according to the longitudinal wave data.
Optionally, in an embodiment of the present invention, before performing multi-round surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data to obtain a longitudinal wave superposition profile, the method includes: and carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
The embodiment of the invention also provides a device for correcting the static of the converted transverse wave detector, which comprises:
the approximate static correction module is used for determining the approximate conversion transverse wave detection point static correction according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the longitudinal wave detection point static correction;
the gather extraction processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave co-asymptote gather data;
the longitudinal wave superposition profile module is used for carrying out multi-wheel 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 profile;
the first wheel superposition profile module is used for carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first wheel rotation transverse wave superposition profile;
the profile stretching module is used for picking up a mark layer on the longitudinal wave superposition profile and the first rotary transverse wave superposition profile respectively, and stretching the longitudinal wave superposition profile by taking the mark layer as a control point;
the second wheel superposition profile module is used for taking the stretched longitudinal wave superposition profile as a first model channel, and carrying out 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 rotation transverse wave superposition profile;
the third wheel superposition profile module is used for taking the second wheel rotation transverse wave superposition profile as a second model track, and carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the conversion transverse wave common asymptote gather data to obtain a third wheel conversion transverse wave superposition profile and a third wheel residual static correction value;
and the static correction value determining module is used for superposing the static correction value of the approximate converted transverse wave detection point and the residual static correction value of the third wheel to obtain the final total static correction amount of the converted transverse wave detection point.
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 corresponding layers of the longitudinal wave and 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 module is used for determining the longitudinal wave detector static correction according to the longitudinal wave data.
Optionally, in an embodiment of the present invention, the apparatus further includes: and the gather sorting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting processing on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the following steps when executing the computer program:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel 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 profile;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first rotary transverse wave superposition section;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
taking the second rotary transverse wave superposition profile as a second 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 third round of converted transverse wave superposition profile and a third round of residual static correction value;
and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total amount.
The embodiment of the invention also provides a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel 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 profile;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first rotary transverse wave superposition section;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
taking the second rotary transverse wave superposition profile as a second 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 third round of converted transverse wave superposition profile and a third round of residual static correction value;
and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector 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 residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer dislocation, larger error and the like caused by the inaccuracy of the layer picked up in the low signal-to-noise ratio area. Therefore, the method has the advantages of achieving a good effect, better solving the problem of static correction of converted waves with an inclined structure and extremely poor signal-to-noise ratio, improving the imaging quality and precision of the converted waves and meeting the requirements of actual production.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions of the prior art, the drawings that are needed in the description of the embodiments will be briefly described below, it being obvious that the drawings in the description below are only some embodiments of the invention and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for static correction of a converted transverse wave detector in accordance with an embodiment of the present invention;
FIG. 2 is a flow chart of a method for converting transverse wave detector static correction in accordance with an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a static correction device for a converted transverse wave detector according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a method and a device for static correction of a converted transverse wave detector.
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Fig. 1 is a flowchart of a method for calibrating statics of a converted transverse wave detector according to an embodiment of the present invention, where the method includes:
step S1, determining a static correction value of an approximate conversion transverse wave detection point according to an average time ratio of a corresponding layer of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
step S2, conventional denoising, amplitude compensation, deconvolution and gather extraction processing are carried out on the converted transverse wave data to which the approximate converted transverse wave detection point static correction value is applied, so as to obtain converted transverse wave co-asymptote gather data;
s3, carrying out multi-wheel 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 profile;
s4, 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 first rotary transverse wave superposition section;
s5, respectively picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
s6, taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
s7, taking the second rotary transverse wave superposition profile as a second model channel, 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 third-wheel converted transverse wave superposition profile and a third-wheel residual static correction value;
and S8, superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total.
As an embodiment of the invention, the method further comprises: and determining the average time ratio of the corresponding layers of the longitudinal wave and the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed.
As an embodiment of the invention, the method further comprises: and determining the static correction amount of the longitudinal wave detector according to the longitudinal wave data.
In this embodiment, before performing multi-round surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data to obtain a longitudinal wave superposition profile, the method includes: and carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
In one embodiment of the present invention, the method of the present invention specifically comprises:
step 1: calculating a static correction value of the longitudinal wave detector;
step 2: counting the ratio gamma of the longitudinal wave speed and the transverse wave speed of the surface shallow layer of the work area;
step 3: calculating the average ratio alpha of the corresponding layers of the surface shallow longitudinal wave and the converted transverse wave in time according to the formula 1:
wherein T is ps To convert the shear wave horizon time, T p Is the longitudinal wave horizon time;
step 4: multiplying the static positive quantity of the longitudinal wave detector point by the average time ratio, and applying the static positive quantity to the converted transverse wave single-shot data;
step 5: performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave single shot data in the step 4 to obtain common asymptote gather (ACP) data;
step 6: performing conventional denoising, amplitude compensation and deconvolution and gather sorting on the longitudinal wave data to obtain common-center point gather (CMP) data;
step 7: performing multi-wheel surface consistency residual static correction processing on the longitudinal wave CMP data in the step 6, and performing superposition processing to obtain a longitudinal wave superposition profile;
step 8: carrying out multi-wheel residual static correction processing on the converted transverse wave ACP gather in the step 5, and carrying out superposition processing to obtain a first converted transverse wave superposition profile;
step 9: and respectively picking up the mark layers on the longitudinal wave superposition profile and the first rotary transverse wave superposition profile. Stretching the longitudinal wave superposition profile by using the picked mark layer as a control point to enable the longitudinal wave superposition profile to correspond to the converted wave superposition profile in time;
step 10: and (3) taking the stretched longitudinal wave superposition profile obtained in the step (9) as a model channel. Carrying out multi-wheel surface consistency residual static correction calculation on the converted transverse wave ACP gather generated in the step 4 again, and carrying out superposition processing to obtain a second wheel rotation transverse wave superposition profile;
step 11: and (3) taking the second-wheel converted transverse wave superposition profile obtained in the step (10) as a model channel, re-carrying out multi-wheel surface consistency residual static correction calculation on the converted transverse wave ACP channel set generated in the step (4), and carrying out superposition processing to obtain a third-wheel converted transverse wave superposition profile and a third-wheel residual static correction.
Step 12: and accumulating the static correction value in the step 4 and the third wheel residual static correction value in the step 11 as a converted transverse wave detector static correction value.
Fig. 2 is a flowchart of a method for correcting the statics of a converted transverse wave detector according to an embodiment of the present invention, and in particular, a flowchart of a method for correcting the statics of a converted transverse wave detector according to an embodiment of the present invention based on consistency of construction. In the general multi-component seismic data processing process, the first work done is the analysis of the data quality. By analyzing the data quality, whether the work area has poor earth surface receiving condition or has an inclined structure at the abdomen of the work area is judged, so that the converted transverse wave signal is weak, the signal to noise ratio is low, and the situation that the image cannot be accurately formed on the co-detection superposition section is caused. If not, adopting a traditional method; if so, the traditional conversion transverse wave static correction method is not suitable for the work area, and the method provided by the invention is needed to be adopted.
According to the method, first arrival picking of the longitudinal wave cannon is carried out on the acquired longitudinal wave single cannon data of the whole work area. And inputs the picked-up first arrival information into static correction software (e.g., geoEast, omega, geoVation, etc.). And calculating various static correction methods such as refraction static correction, chromatography static correction and the like by utilizing the advantages of respective software according to requirements. Applying different types of static correction values to longitudinal wave data, selecting static correction value with optimal application effect, and recording as ST p . At this time ST p Comprising two components, one being the detector component ST p R, one is shot component ST p S. For example, after the chromatographic static correction value is applied, compared with other static correction values, the longitudinal wave single shot reflection phase axis is smooth, the jitter phenomenon is weakened, and the superposition section phase axis is continuous. I.e. the amount of tomosynthesis correction is selected. And vice versa.
After the longitudinal wave static correction amount is determined, the ratio gamma of the longitudinal wave speed and the transverse wave speed of the surface shallow layer of the work area is counted according to geology, 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 there is some correlation between the two, γ is generally greater than 1, ranging from about 1 to 3.
In seismic data, the recording of seismic waves is on the ordinate of time, and thus the amount of static correction is also a time value. Which represents the amount of time that each seismic wave data trace has been shifted up or down. The difference in velocity necessarily varies in propagation time in the case where the propagation paths of the two waves are the same. Therefore, the average ratio alpha of the corresponding layers of the surface shallow longitudinal wave and the converted transverse wave in time is calculated by using the formula 1. For example, the ratio of the longitudinal and transverse wave velocities of the shallow layer of the work area is 2, and the time ratio of the longitudinal and transverse wave profile of the shallower reflecting layer, particularly the emitting layer near the ground surface, is 1.5. If a reflective layer X is about 1 second in longitudinal section, it is about 1.5 seconds in transverse section.
Wherein T is ps To convert the shear wave horizon time, T p Is the longitudinal wave horizon time;
after the time ratio is obtained, the selected longitudinal wave static correction amount ST is used p The detector component ST of (1) p Multiplying R by the time ratio α to obtain an approximate converted transverse wave detector static correction ST ps R, namely: ST (ST) ps _R＝α×ST p R; whereas static correction ST of converted transverse wave shot components ps S still adopts the static correction value ST of the longitudinal wave shot point component p S, i.e. ST ps _S＝ST p S. In the actual processing process, especially under the condition of inclined structure or low data signal-to-noise ratio, the approximate static correction amount is utilized to carry out superposition imaging processing, so that the phenomenon that the structure is inconsistent with the longitudinal wave, the horizon and the phase are dislocated easily occurs, and the actual geological condition is not met.
Thus, according to the flow, the optimal static correction amount ST is applied p Conventional noise attenuation of longitudinal wave dataSubtracting, earth surface consistency amplitude compensation, deconvolution, speed analysis and gather arrangement and selection to form a common center point (CMP) gather, and carrying out multi-round earth surface consistency residual static correction P_Res_st, post-stack denoising and the like to obtain the superposition profile STK1_P of the optimal longitudinal wave.
Also, the approximate static correction amount ST to be multiplied by the time ratio ps R is applied to the converted transverse wave data and is subjected to conventional noise attenuation, surface consistency amplitude compensation, deconvolution, velocity analysis and gather sorting to form a common Asymptote (ACP) gather. And then, based on the self data of the converted transverse wave, carrying out N rounds of residual static correction processing according to the convergence condition of the residual static correction iteration of the ground surface consistency to obtain the total residual static correction PS_Res_st1, completing post-stack denoising, and improving the signal to noise ratio so as to obtain the superimposed profile STK1_PS of the first rotary transverse wave.
Several marker horizons are picked up from shallow to deep on the optimal longitudinal wave stack profile and the first rotated transverse wave stack profile. For example 2 or 3, depending on the actual requirements of the material. From the above, it can be seen that the geological horizon time in the converted transverse wave is greater than the corresponding horizon time in the longitudinal wave. Stretching the longitudinal wave superposition profile by using the picked-up mark layer as a control point, and then carrying out interlayer detail matching on the basis of matching the mark layer, namely completing the time matching processing of the longitudinal wave and the transverse wave, so that the longitudinal wave and the transverse wave correspond to each other in time.
And then, taking the stretched longitudinal wave superposition profile as a model channel, and carrying out multi-round surface consistency residual static correction calculation on the converted transverse wave ACP channel set again to obtain the total residual static correction PS_Res_st2. And performing superposition and post-superposition denoising treatment to obtain a second rotary transverse wave superposition profile STK2_PS. The aim of this idea is to guarantee the requirements of consistency with the longitudinal wave structure.
And (3) taking the obtained second rotary transverse wave superposition profile STK2_PS as a model channel, and carrying out multi-round surface consistency residual static correction calculation again on the converted transverse wave ACP channel set to obtain the total residual static correction PS_Res_st3. And performing superposition processing to obtain a third-round converted transverse wave superposition profile STK3_PS.
When the above steps are completed, the approximate converted transverse wave is detectedWave point static correction value ST ps Sum of_R and the third residual static correction total PS_Res_st3 is used as the converted transverse wave detector static correction ST ps Final, as shown in equation 2. And the calculation of the static correction of the converted transverse wave detector is completed and output.
ST ps _Final＝ST ps R+ps_res_st3_2nd
In summary, 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 imaging effect of the data by carrying out multi-round surface consistency residual static correction iterative processing 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 residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer dislocation, larger error and the like caused by the inaccuracy of the layer picked up in the low signal-to-noise ratio area. Therefore, the method has the advantages of achieving a good effect, better solving the problem of static correction of converted waves with an inclined structure and extremely poor signal-to-noise ratio, improving the imaging quality and precision of the converted waves and meeting the requirements of actual production.
Fig. 3 is a schematic structural diagram of a static correction device for a converted transverse wave detector, according to an embodiment of the present invention, where the device includes:
the approximate static correction module 10 is configured to determine an approximate converted transverse wave detector static correction according to an average time ratio of a corresponding horizon of the longitudinal wave and the converted transverse wave and the longitudinal wave detector static correction;
the gather extraction processing module 20 is 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 detector static correction is applied, so as to obtain converted transverse wave co-asymptotic gather data;
the longitudinal wave superposition profile module 30 is used for performing multi-wheel 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 profile;
the first-wheel superposition profile module 40 is configured to perform multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, so as to obtain a first-wheel rotation transverse wave superposition profile;
a section stretching module 50, configured to pick up a marker layer on the longitudinal wave superposition section and the first rotated transverse wave superposition section, and stretch the longitudinal wave superposition section by using the marker layer as a control point;
the second wheel superposition profile module 60 is configured to use the stretched longitudinal wave superposition profile as a first model channel, perform multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and obtain a second wheel rotation transverse wave superposition profile;
a third-wheel superposition profile module 70, configured to use the second-wheel rotation transverse wave superposition profile as a second model track, and perform multi-wheel surface consistency residual static correction processing and superposition processing on the conversion transverse wave common asymptote gather data to obtain a third-wheel conversion transverse wave superposition profile and a third-wheel residual static correction value;
and the static correction value determining module 80 is configured to superimpose the approximate converted transverse wave detector static correction value and the third wheel residual static correction value, so as to obtain a final converted transverse wave detector 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 corresponding layers of the longitudinal wave and 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 module is used for determining the longitudinal wave detector static correction according to the longitudinal wave data.
In this embodiment, the apparatus further includes: and the gather sorting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting processing on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
Based on the same application conception as the method for correcting the converted transverse wave detector point, the invention also provides a device for correcting the converted transverse wave detector point. Because the principle of the converted transverse wave detector static correction device for solving the problem is similar to that of the converted transverse wave detector static correction method, the implementation of the converted transverse wave detector static correction device can refer to the implementation of the converted transverse wave detector 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 residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer dislocation, larger error and the like caused by the inaccuracy of the layer picked up in the low signal-to-noise ratio area. Therefore, the method has the advantages of achieving a good effect, better solving the problem of static correction of converted waves with an inclined structure and extremely poor signal-to-noise ratio, improving the imaging quality and precision of the converted waves and meeting the requirements of actual production.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the following steps when executing the computer program:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel 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 profile;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first rotary transverse wave superposition section;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
taking the second rotary transverse wave superposition profile as a second 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 third round of converted transverse wave superposition profile and a third round of residual static correction value;
and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total amount.
The embodiment of the invention also provides a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel 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 profile;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a first rotary transverse wave superposition section;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
taking the stretched longitudinal wave superposition profile as a first model channel, and performing multi-wheel surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data to obtain a second rotary transverse wave superposition profile;
taking the second rotary transverse wave superposition profile as a second 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 third round of converted transverse wave superposition profile and a third round of residual static correction value;
and superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total amount.
Based on the same application conception as the method for correcting the converted transverse wave detector point static, the invention also provides the computer equipment and the computer readable storage medium. Because the principle of the computer device and the computer readable storage medium for solving the problems is similar to that of a method for correcting the statics of the converted transverse wave detector, the implementation of the computer device and the computer readable storage medium can refer to the implementation of the method for correcting the statics of the converted transverse wave detector, and the repetition is omitted.
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 residual static correction calculation of the earth surface consistency, thereby avoiding the problems of layer dislocation, larger error and the like caused by the inaccuracy of the layer picked up in the low signal-to-noise ratio area. Therefore, the method has the advantages of achieving a good effect, better solving the problem of static correction of converted waves with an inclined structure and extremely poor signal-to-noise ratio, improving the imaging quality and precision of the converted waves and meeting the requirements of actual production.
Those of ordinary skill in the art will appreciate that all or a portion of the steps in implementing the methods of the above embodiments may be implemented by a program that instructs associated hardware, and the program may be stored on a computer readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (6)
1. A method for static correction of a converted transverse wave detector, the method comprising:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data, and carrying out noise attenuation on a superposition processing result to obtain a longitudinal wave superposition profile with high signal-to-noise ratio;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and carrying out noise attenuation on a superposition processing result to obtain a first rotary transverse wave superposition section with high signal-to-noise ratio;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
the stretched longitudinal wave superposition profile is used as a first model channel, the converted transverse wave common asymptote gather data is subjected to multi-round earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a second rotary transverse wave superposition profile with high signal to noise ratio is obtained;
the second wheel rotating transverse wave superposition section is used as a second model road, the conversion transverse wave common asymptote gather data is subjected to multi-wheel earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a third wheel of conversion transverse wave superposition section with high signal to noise ratio and a third wheel of residual static correction value are obtained;
superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total;
wherein the method further comprises: according to the ratio of the longitudinal wave speed to the transverse wave speed, the average time ratio of the corresponding layers of the longitudinal wave and the converted transverse wave is determined by using the following formula:
wherein T is ps To convert the shear wave horizon time, T p For longitudinal wave horizon time, γ is the ratio of longitudinal wave velocity to transverse wave velocity.
2. The method of claim 1, wherein the performing the multi-round surface consistency residual static correction and superposition on the longitudinal wave common center point gather data to obtain the longitudinal wave superposition profile comprises: and carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
3. A converted transverse wave detector static correction device, said device comprising:
the approximate static correction module is used for determining the approximate conversion transverse wave detection point static correction according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the longitudinal wave detection point static correction;
the gather extraction processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave co-asymptote gather data;
the longitudinal wave superposition profile module is used for carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data, and carrying out noise attenuation on the superposition processing result to obtain a longitudinal wave superposition profile with high signal-to-noise ratio;
the first wheel superposition profile module is used for carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and carrying out noise attenuation on a superposition processing result to obtain a first wheel rotation transverse wave superposition profile with high signal-to-noise ratio;
the profile stretching module is used for picking up a mark layer on the longitudinal wave superposition profile and the first rotary transverse wave superposition profile respectively, and stretching the longitudinal wave superposition profile by taking the mark layer as a control point;
the second wheel superposition profile module is used for taking the stretched longitudinal wave superposition profile as a first model channel, re-carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and carrying out noise attenuation on the superposition processing result to obtain a second wheel rotation transverse wave superposition profile with high signal to noise ratio;
the third wheel superposition profile module is used for taking the second wheel conversion transverse wave superposition profile as a second model track, re-carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the conversion transverse wave common asymptote gather data, and carrying out noise attenuation on a superposition processing result to obtain a third wheel conversion transverse wave superposition profile with high signal to noise ratio and a third wheel residual static correction value;
the static correction value determining module is used for superposing the static correction value of the approximate converted transverse wave detection point and the residual static correction value of the third wheel to obtain the final total static correction amount of the converted transverse wave detection point;
wherein the apparatus further comprises: the average time ratio module is used for determining the average time ratio of the corresponding layers of the longitudinal wave and the converted transverse wave according to the ratio of the longitudinal wave speed to the transverse wave speed by using the following formula:
wherein T is ps In order to switch the shear wave horizon time,T p for longitudinal wave horizon time, γ is the ratio of longitudinal wave velocity to transverse wave velocity.
4. A device according to claim 3, characterized in that the device further comprises: and the gather sorting processing module is used for carrying out conventional denoising, amplitude compensation, deconvolution and gather sorting processing on the longitudinal wave data to which the longitudinal wave detector static correction value is applied to obtain longitudinal wave common center point gather data.
5. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the following steps when executing the computer program:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data, and carrying out noise attenuation on a superposition processing result to obtain a longitudinal wave superposition profile with high signal-to-noise ratio;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and carrying out noise attenuation on a superposition processing result to obtain a first rotary transverse wave superposition section with high signal-to-noise ratio;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
the stretched longitudinal wave superposition profile is used as a first model channel, the converted transverse wave common asymptote gather data is subjected to multi-round earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a second rotary transverse wave superposition profile with high signal to noise ratio is obtained;
the second wheel rotating transverse wave superposition section is used as a second model road, the conversion transverse wave common asymptote gather data is subjected to multi-wheel earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a third wheel of conversion transverse wave superposition section with high signal to noise ratio and a third wheel of residual static correction value are obtained;
superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total;
wherein the steps further comprise: according to the ratio of the longitudinal wave speed to the transverse wave speed, the average time ratio of the corresponding layers of the longitudinal wave and the converted transverse wave is determined by using the following formula:
wherein T is ps To convert the shear wave horizon time, T p For longitudinal wave horizon time, γ is the ratio of longitudinal wave velocity to transverse wave velocity.
6. A computer readable storage medium having stored thereon a computer program, characterized in that the computer program when executed by a processor realizes the steps of:
determining a static correction value of the approximate conversion transverse wave detection point according to the average time ratio of the corresponding layers of the longitudinal wave and the conversion transverse wave and the static correction value of the longitudinal wave detection point;
performing conventional denoising, amplitude compensation, deconvolution and gather extraction processing on the converted transverse wave data to which the approximate converted transverse wave detector static correction value is applied to obtain converted transverse wave common asymptote gather data;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the longitudinal wave common center point gather data, and carrying out noise attenuation on a superposition processing result to obtain a longitudinal wave superposition profile with high signal-to-noise ratio;
carrying out multi-wheel earth surface consistency residual static correction processing and superposition processing on the converted transverse wave common asymptote gather data, and carrying out noise attenuation on a superposition processing result to obtain a first rotary transverse wave superposition section with high signal-to-noise ratio;
picking up a mark layer on the longitudinal wave superposition section and the first rotary transverse wave superposition section respectively, and stretching the longitudinal wave superposition section by using the mark layer as a control point;
the stretched longitudinal wave superposition profile is used as a first model channel, the converted transverse wave common asymptote gather data is subjected to multi-round earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a second rotary transverse wave superposition profile with high signal to noise ratio is obtained;
the second wheel rotating transverse wave superposition section is used as a second model road, the conversion transverse wave common asymptote gather data is subjected to multi-wheel earth surface consistency residual static correction processing and superposition processing again, noise attenuation is carried out on superposition processing results, and a third wheel of conversion transverse wave superposition section with high signal to noise ratio and a third wheel of residual static correction value are obtained;
superposing the approximate conversion transverse wave detector static correction value and the third wheel residual static correction value to obtain the final conversion transverse wave detector static correction total;
wherein the steps further comprise: according to the ratio of the longitudinal wave speed to the transverse wave speed, the average time ratio of the corresponding layers of the longitudinal wave and the converted transverse wave is determined by using the following formula:
wherein T is ps To convert the shear wave horizon time, T p For longitudinal wave horizon time, γ is the ratio of longitudinal wave velocity to transverse wave velocity.
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 CN112305613A (en) | 2021-02-02 |
| CN112305613B true 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 (7)
| 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 |
| 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 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9348049B2 (en) * | 2012-01-05 | 2016-05-24 | Cgg Services Sa | Simultaneous joint estimation of the P-P and P-S residual statics |
| US10386519B2 (en) * | 2015-12-18 | 2019-08-20 | 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 (7)
| 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 |
| 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 |
| 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.Exploration Geophysics.2014,全文. * |
| 杨海涛 ; 贺振华 ; 陈丹 ; 王栋 ; 王静 ; 王鸿燕 ; .转换波静校正技术在四川地区的应用研究.地球物理学进展.2017,(第04期),全文. * |
| 王淑荣等.转换波处理方法研究.石油地球物理勘探.2008,全文. * |
| 陈海峰等.联合纵波确定转换波接收点大的静校正量.地球物理学进展.2014,全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112305613A (en) | 2021-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111474584B (en) | Focusing superposition imaging method and system based on correlated seismic interference | |
| US10557956B2 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
| CN101614826B (en) | Method and device for realizing binning homogenization in three-dimensional seismic data processing | |
| Si et al. | Attenuation of random noise using denoising convolutional neural networks | |
| CN105182408A (en) | Manufacturing method and device for synthesizing earthquake record | |
| CN104090301B (en) | A kind of method for asking for three-D high frequency static correction value | |
| CN109917454A (en) | True earth's surface pre-stack depth migration imaging method and device based on double datum levels | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| CN107229071B (en) | A kind of subsurface structure inversion imaging method | |
| CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
| CN106353798A (en) | Multi-component joint Gaussian beam pre-stack reverse-time migration imaging method | |
| CN112305613B (en) | Static correction method and device for converted transverse wave detector | |
| CN106855639A (en) | Static correcting method and device | |
| CN104155688A (en) | High precision weighted stack method | |
| CN106338760B (en) | The relief surface offset method of error compensation | |
| CN115421195B (en) | Method, device, equipment and storage medium for generating velocity field in seismic exploration | |
| CN109581494B (en) | Pre-stack migration method and device | |
| CN116125530A (en) | Loess tableland seismic data processing method, loess tableland seismic data processing system and electronic equipment | |
| CN105445795A (en) | Method and system for analyzing velocity of tilted interface | |
| Li et al. | Progressive multitask learning for high-resolution prediction of reservoir elastic parameters | |
| CN113791447A (en) | Reflection wave chromatography inversion method guided by reflection structure | |
| CN107193047A (en) | Reflection static corr ections time-varying dynamic residual static corrections and device | |
| CN111596359B (en) | Inclined stratum wide line data correction imaging method and device | |
| CN111208568A (en) | Time domain multi-scale full waveform inversion method and system | |
| CN114415239B (en) | Earthquake reflection wave combined imaging method and system |
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 |