CN115963565B - Seabed multiple attenuation method and device based on hyperbolic vector median filter - Google Patents
Seabed multiple attenuation method and device based on hyperbolic vector median filter Download PDFInfo
- Publication number
- CN115963565B CN115963565B CN202310098387.8A CN202310098387A CN115963565B CN 115963565 B CN115963565 B CN 115963565B CN 202310098387 A CN202310098387 A CN 202310098387A CN 115963565 B CN115963565 B CN 115963565B
- Authority
- CN
- China
- Prior art keywords
- wave
- hyperbolic
- point
- submarine
- correction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 74
- 239000013598 vector Substances 0.000 title claims abstract description 70
- 238000012937 correction Methods 0.000 claims abstract description 100
- 238000003384 imaging method Methods 0.000 claims abstract description 65
- 238000001228 spectrum Methods 0.000 claims abstract description 44
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 25
- 230000005012 migration Effects 0.000 claims abstract description 15
- 238000013508 migration Methods 0.000 claims abstract description 15
- 230000002238 attenuated effect Effects 0.000 claims abstract description 11
- 238000012876 topography Methods 0.000 claims abstract description 9
- 239000013535 sea water Substances 0.000 claims description 22
- 238000001514 detection method Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 4
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000011410 subtraction method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- 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/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a seabed multiple attenuation method and device based on a hyperbolic vector median filter, and belongs to the technical field of seismic waves. The method comprises the following steps: performing water depth correction on the seismic channel according to P, S wave offset speed and submarine topography data; pre-stack time migration is carried out on the amplitudes of all the sample points in the seismic channel after the water depth correction, and non-hyperbolic correction is carried out on the PS wave travel time of all the sample points according to the submarine incident point positions, so that a PS wave non-hyperbolic correction co-imaging gather is obtained; determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition of focusing velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points; and carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window. By the method, the first-order and Gao Jiehai bottom multiples on the common imaging gather can be effectively attenuated, and obvious efficiency advantages are achieved.
Description
Technical Field
The embodiment of the application relates to the technical field of seismic waves, in particular to a seabed multiple attenuation method and device based on a hyperbolic vector median filter.
Background
Attenuation of the ocean bottom multiples is an important link in seismic data processing, and particularly for ocean bottom seismic data, free surface multiples develop severely, which can cause imaging artifacts of the seismic data if not attenuated. As shown in FIG. 1, sea surface seismic waves originate from shot point O at an acoustic velocity v W Incident on the sea floor at a W point and a P wave velocity v P Reflected wave mode conversion propagating to interface C point at S wave velocity v S The wave field to be damped is first-order multiple and higher-order multiple indicated by other paths.
The existing attenuation method of the submarine multiples comprises the following steps: the method comprises a prediction deconvolution method, a method for suppressing multiple based on the time difference characteristic of primary reflected wave and multiple, a prediction subtraction method based on wave equation and a median filtering subtraction method, wherein the median filtering method has wide application in free surface multiple attenuation and can be implemented on shot gathers, common detection point gathers or imaging gathers.
However, since the wave pattern of the submarine multiples can be converted for multiple times in the propagation process, and the difficulty in identifying the S-type multiples is high in shot gather records and common-detector gather records, especially the high-order S-type multiples, the conventional median filtering method is difficult to directly use for attenuating the submarine multiples, and the attenuation effect is not ideal.
Disclosure of Invention
The embodiment of the application provides a seabed multiple attenuation method and device based on a hyperbolic vector median filter, which are used for attenuating seabed multiple on a co-imaging gather.
In order to solve the technical problems, the application is realized as follows:
in a first aspect, an embodiment of the present application provides a method for attenuating a plurality of seafloor waves based on a hyperbolic vector median filter, including:
carrying out water depth correction on the seismic channel according to P, S wave migration velocity and submarine topography data to obtain the seismic channel after water depth correction and submarine incident point positions corresponding to various points of the seismic channel;
performing prestack time migration on the amplitudes of all the sample points in the seismic channel after the water depth correction, and performing non-hyperbolic correction on the PS wave travel time of all the sample points according to the submarine incident point positions to obtain a PS wave non-hyperbolic correction co-imaging gather;
determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition of focusing velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points;
and carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, wherein after vector median filtering, the submarine multiples developed on the PS wave non-hyperbolic correction co-imaging gather are attenuated.
In a second aspect, an embodiment of the present application further provides a device for attenuating a submarine multiple based on a hyperbolic vector median filter, including:
the water depth correction module is used for carrying out water depth correction on the seismic channel according to the P, S wave offset speed and the submarine topography data to obtain the seismic channel after the water depth correction and the submarine incident point positions corresponding to various points of the seismic channel;
the non-hyperbolic correction module is used for carrying out pre-stack time migration on the amplitudes of all the points in the seismic channel after the water depth correction, and carrying out non-hyperbolic correction on the PS wave travel time of all the points according to the submarine incident point position to obtain a PS wave non-hyperbolic correction co-imaging gather;
the hyperbolic time window construction module is used for determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition that velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points are focused;
and the multiple attenuation module is used for carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, and attenuating the submarine multiple developed on the PS wave non-hyperbolic correction co-imaging gather after the vector median filtering.
In a third aspect, embodiments of the present application further provide a computer device, where the computer device includes a processor and a memory, where the memory stores at least one computer program, and the at least one computer program is loaded and executed by the processor to implement the above-mentioned method for attenuating a subsea multiple based on a hyperbolic vector median filter.
In a fourth aspect, embodiments of the present application further provide a computer readable storage medium, where at least one computer program is stored, where the computer program is loaded and executed by a processor to implement the above-mentioned method for ocean bottom multiple attenuation based on a hyperbolic vector median filter.
The technical scheme that this application provided can include following beneficial effect:
according to the embodiment of the application, the seismic channel is subjected to water depth correction according to P, S wave offset speed and submarine topography data; pre-stack time migration is carried out on the amplitudes of all the sample points in the seismic channel after the water depth correction, non-hyperbolic correction is carried out on the PS wave travel time of all the sample points according to the submarine incident point positions, a PS wave non-hyperbolic correction common imaging gather is obtained, and the primary reflection PS wave and submarine multiple can be effectively distinguished.
Further, determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition of focusing velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points; according to the hyperbolic time window, vector median filtering is carried out on the PS wave non-hyperbolic correction co-imaging gather, so that first-order and Gao Jiehai bottom multiples on the co-imaging gather can be attenuated at one time, and the method has obvious efficiency advantages.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the application and together with the description, serve to explain the principles of the application.
FIG. 1 shows a schematic diagram of primary reflected waves and first order multiple modes received subsea;
FIG. 2 shows a schematic flow chart of a method for attenuating a plurality of ocean bottom waves based on a hyperbolic vector median filter according to an embodiment of the present application;
FIG. 3 is a schematic flow chart of another method for attenuating a plurality of ocean bottom waves based on a hyperbolic vector median filter according to an embodiment of the present application;
FIG. 4 shows a schematic diagram of a non-hyperbolic corrected PS wave CIG and its velocity spectrum provided by an embodiment of the present application;
FIG. 5 shows an effect plot of non-hyperbolic corrected seafloor PS-wave co-imaging gathers versus velocity spectrum provided by embodiments of the present application;
FIG. 6 shows an effect diagram of a submarine PS-wave co-imaging gather and velocity spectrum suppressing S-multiples provided by an embodiment of the present application;
fig. 7 shows a schematic structural diagram of a submarine multiple attenuation device based on a hyperbolic vector median filter according to an embodiment of the present application;
FIG. 8 shows a schematic structural diagram of another submarine multiple attenuation device based on a hyperbolic vector median filter according to an embodiment of the present application;
fig. 9 shows a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples are not representative of all implementations consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with some aspects of the present application as detailed in the accompanying claims.
FIG. 1 shows a schematic diagram of primary reflected wave and first order multiples received at the seafloor, as shown in the figure, with paths O W C D being the imaged target wavefield and the wavefield requiring attenuation being the first order multiples and higher order multiples indicated by the other paths. These multiples are free surface multiples that propagate in the sea in a P-wave mode, and the propagation between the reflective interface and the sea floor can be either a P-wave mode or an S-wave mode; but eventually returns from the reflective interface to the hydrophone in S-wave mode and may be referred to as S-wave multiples.
Assume that the horizontal distances between OW, WC, CD are x, respectively 1 、x 2 、x 3 The method comprises the steps of carrying out a first treatment on the surface of the The vertical distance z between the W point and the sea surface W The vertical distance from the C point to the seabed is z C The method comprises the steps of carrying out a first treatment on the surface of the PS-wave travel time t of OWCD segment ray path od The method comprises the following steps:
wherein, the earthquake wave travel time t in the sea water w The method comprises the following steps:
here, v P For P-wave velocity, v S Is the S-wave velocity.
Typically, P-wave velocity and sea water depth are known information before PS-wave imaging, and S-wave velocity needs to be iteratively updated in the co-imaging gather (Common Image Gather, CIG).
The embodiment of the application provides a seabed multiple attenuation method based on a hyperbolic vector median filter, which realizes attenuation of S-shaped multiple while updating S-wave speed, as shown in fig. 2, and the method 200 specifically comprises the following steps:
s201: carrying out water depth correction on the seismic channel according to P, S wave migration velocity and submarine topography data to obtain the seismic channel after water depth correction and submarine incident point positions corresponding to various points of the seismic channel;
s202: performing prestack time migration on the amplitudes of all the sample points in the seismic channel after the water depth correction, and performing non-hyperbolic correction on the PS wave travel time of all the sample points according to the submarine incident point positions to obtain a PS wave non-hyperbolic correction co-imaging gather;
s203: determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition of focusing velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points;
s204: and carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, wherein after vector median filtering, the submarine multiples developed on the PS wave non-hyperbolic correction co-imaging gather are attenuated.
The inventors found that the time difference of the primary reflected PS wave has a hyperbolic characteristic capable of highlighting the distinction of the primary reflected PS wave from the S-shaped ocean bottom multiples. Based on this, the embodiment of the application proposes to perform the depth correction on the seismic trace, perform the pre-stack time migration on the amplitudes of all the sample points in the seismic trace after the depth correction, and perform the non-hyperbolic correction on the PS wave travel time of each sample point according to the submarine incident point position, so as to obtain the PS wave non-hyperbolic correction co-imaging gather, so that the primary reflection PS wave can be effectively distinguished from the submarine multiple.
Further, according to the velocity spectrum corresponding to the PS wave co-imaging gather after non-hyperbolic correction (namely, the PS wave non-hyperbolic correction co-imaging gather), a hyperbolic time window for vector median filtering is constructed by sample points, and after the hyperbolic vector median filtering is carried out on the PS wave co-imaging gather, first-order and high-order S-type multiples can be attenuated at one time.
In one possible embodiment, depth correction of seismic traces from seismic wave migration velocity and seafloor terrain data comprises: acquiring the position coordinates, offset distance, sea water depth at the shot point end, P wave offset speed and S wave offset speed of the shot point and the wave detector of the seismic channel; determining the position coordinates of the submarine incident point according to the offset distance, the sea water depth of the shot point end, the P wave offset speed and the S wave offset speed; determining a second horizontal distance between the submarine incident point and the PS wave reflection point and a third horizontal distance between the PS wave reflection point and the detection point according to the first horizontal distance between the gun point and the submarine incident point; calculating a sum of the first horizontal distance, the second horizontal distance and the third horizontal distance; resetting the position coordinates of the submarine incident point by adopting a dichotomy under the condition that the difference value between the sum value and the offset exceeds a preset threshold value until the calculated difference value between the sum value and the offset is smaller than the preset threshold value; according to the reset position coordinates of the submarine incident points, the seismic wave travel time in the sea water is solved, and according to the seismic wave travel time in the sea water, the time of each sample point of the seismic channel is corrected.
In a specific implementation, as shown in fig. 3, first a grid of offset imaging is set, based on which PS wave co-imaging gathers will be generated by pre-stack time offset. Inputting an earthquake channel, and acquiring relevant parameters of the earthquake channel, wherein the method specifically comprises the following steps of: position coordinates of shot point and wave detector, offset distance l and sea water depth z of shot point end w P-wave offset velocity v P S-wave offset velocity v S 。
Here, the S-wave shift velocity v S Is the initial S-wave offset speed, and is updated in the CIG in a subsequent iteration.
Correcting the influence of sea depth on seismic wave travel time by sample points through the related parameters of the seismic channels, namely moving shot point O to seabed W point, and respectively determining x 1 、x 2 、x 3 Is of a size of (a) and (b).
According to the offset distance, the sea water depth of the shot point end, the P wave offset speed and the S wave offset speed, the position coordinates of the submarine incident point are determined, and the specific formula is as follows:
wherein,,
here, T 0 Is a double journey travel time when the seismic wave offset is 0.
Further, a second horizontal distance x between the incident point of the sea bottom and the reflection point of the PS wave can be calculated 2 And a third horizontal distance x between the PS wave reflection point and the detection point 3 The specific calculation formula is as follows:
wherein,,
then there are:
l'＝x 1 +x 2 +x 3 ，
comparing the sum value l' with the offset distance l, and resetting the position coordinates of the submarine incident point by adopting a dichotomy under the condition that the difference value between the sum value and the offset distance exceeds a preset threshold value until the calculated difference value between the sum value and the offset distance is smaller than the preset threshold value, wherein the difference value is understood to be close enough to the l.
According to twoX determined by division 1 Solving the seismic wave travel time t in the sea water for the input seismic channel sample by sample W Correcting the time of each sample point of the seismic channel according to the seismic wave travel time in the sea water so as to eliminate the influence of the water depth on the time of each sample point of the input seismic channel, so that the PS wave travel time t WD The method meets the following conditions:
the water depth correction of the seismic channel is completed through the steps, and x corresponding to each sample point of the seismic channel is output 1 The method comprises the steps of carrying out a first treatment on the surface of the If the sea bottom has the problem of uneven fluctuation, the shot point and the wave detection point can be corrected to the sea bottom datum plane by adopting a static correction method.
Further, according to the submarine incident point position, performing non-hyperbolic correction on PS wave travel time of each sample point, including: determining the horizontal distance between the submarine incident point and the detection point according to the submarine incident point position; and correcting the PS wave travel time of each sample point to a hyperbolic time difference position according to the horizontal distance, the vertical distance of the submarine and PS wave reflection interface and the reflected wave imaging speed.
In specific implementation, for the seismic trace after the water depth correction, moving the amplitude into the CIG by the prestack time offset sample by sample, and correcting the PS wave travel time to the hyperbolic time difference position by the following formula, namely finishing the non-hyperbolic correction of the PS wave CIG:
wherein t is WD For PS wave travel time, Δt WD Is a non-hyperbolic moveout correction amount.
As shown in fig. 4, the velocity spectrum energy bin of the non-hyperbolically corrected primary reflected PS wave is focused, and the multiple S waves form a lower velocity energy bin, which would occur in a lower velocity interval if there were higher order multiple S waves. On the velocity spectrum of fig. 4, primary reflected wave imaging velocity points are picked up(T 0 ,v c ) Here, v c For a double travel time T when the seismic wave offset is 0 in the velocity spectrum 0 Corresponding speeds. In one possible implementation, in the case that the velocity spectrum energy cluster corresponding to the primary reflection PS wave velocity point is not focused, the S wave offset velocity is updated according to the P wave offset velocity and the velocity of the picked up primary reflection PS wave, and the calculation formula is as follows:
a new PS wave non-hyperbolic corrected CIG is then generated whose focusing of velocity spectrum energy clusters of the primary reflected PS wave will be improved.
Here, whether the velocity spectrum energy bolus is focused may be determined empirically, or the focusing property of the velocity spectrum energy bolus may be determined by calculating specific values such as average amplitude energy and cross-correlation coefficient, which is not particularly limited.
Further, the velocity spectrum of the PS wave CIG regenerated after updating the S wave velocity is picked up (T 0 ,v C ) As shown by the grey dotted line in fig. 4, a continuous hyperbolic time window t is generated over the CIG:
where y is the offset of the CIG.
And carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, wherein the vector median filtering comprises the following steps: acquiring a target seismic channel, taking the target seismic channel as a center, and selecting a preset number of seismic channels from the PS wave non-hyperbolic correction co-imaging channel set to form a vector set; calculating the sum of the distances between the amplitude vector of the target sample point and the amplitude vectors of other sample points except the target sample point in the hyperbolic time window, and outputting the amplitude vector with the minimum sum of the distances between the amplitude vector and the amplitude vector of all the sample points.
In an embodiment, the seismic traces (offset y) of the input CIG are aligned with each other toTaking the same CIG as the center, selecting 1+2n seismic channel composition vector sets (n is a set numerical value), taking the channel spacing as delta y, and calculating the amplitude vector of each sample point in a hyperbolic time window by sample point
Amplitude vector +.>
Wherein the sum of the distances from all sample amplitude vectors and the smallest amplitude vector is taken as the output of the hyperbolic vector median filter:
wherein,,
y 1 ,y 2 ∈{y-nΔy,...,y-Δy,y,y+Δy,...,y+nΔy}.
after hyperbolic vector median filtering, the S-shaped multiple developed on the PS-wave co-imaging gather can be attenuated, wherein the multiple comprises a first-order multiple and a high-order multiple.
As shown in fig. 5, after the co-imaging gather of the submarine PS wave is subjected to non-hyperbolic correction, the energy of the PS wave reflected once on the velocity spectrum is well focused, such as a white dotted line in fig. 5, which is a picked-up (T0, vC) point; before multiple S-wave compression, the first order (at the solid white box) on the PS-wave co-imaging gather is developed in comparison with the high order (at the dashed white box) multiple S-waves; there are distinct multiples of the S-wave wavefield on the co-imaging gather, as indicated by the black arrow. After the method provided by the embodiment of the application is adopted to attenuate multiple S waves, as shown in fig. 6, energy groups on a velocity spectrum are mainly focused at the primary reflection S waves, the multiple S waves of each order are effectively attenuated, and the multiple S waves are suppressed on a common imaging gather.
Fig. 7 shows a schematic structural diagram of a device for attenuating a submarine multiple based on a hyperbolic vector median filter according to an embodiment of the present application, where the device 700 includes:
the water depth correction module 710 is configured to perform water depth correction on the seismic trace according to the P, S wave offset speed and the submarine topography data, so as to obtain a seismic trace after the water depth correction and a submarine incident point position corresponding to each sample point of the seismic trace;
the non-hyperbolic correction module 720 is configured to perform pre-stack time migration on amplitudes of various points in the seismic trace after the depth correction, and perform non-hyperbolic correction on PS wave travel times of the various points according to the submarine incident point positions, so as to obtain a PS wave non-hyperbolic correction co-imaging gather;
the hyperbolic time window construction module 730 is configured to determine a velocity spectrum corresponding to the PS-wave non-hyperbolic correction co-imaging gather, pick up a primary reflection PS-wave velocity point from the velocity spectrum, and construct a hyperbolic time window under the condition that a velocity spectrum energy cluster corresponding to the primary reflection PS-wave velocity point is focused;
and the multiple attenuation module 740 is configured to perform vector median filtering on the PS-wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, and after the PS-wave non-hyperbolic correction co-imaging gather is subjected to vector median filtering, the developed seafloor multiple is attenuated.
The apparatus 700 provided in the embodiment of the present application may perform the method described in fig. 2 and implement the functions of the embodiment shown in fig. 2, which are not described herein.
Fig. 8 shows a schematic structural diagram of another submarine multiple attenuation device based on a hyperbolic vector median filter according to an embodiment of the present application, where the device 800 includes:
the water depth correction module 810 is configured to perform water depth correction on the seismic trace according to the P, S wave offset speed and the submarine topography data, so as to obtain a seismic trace after the water depth correction and a submarine incident point position corresponding to each sample point of the seismic trace;
the non-hyperbolic correction module 820 is used for performing pre-stack time migration on the amplitudes of various points in the seismic trace after the water depth correction, and performing non-hyperbolic correction on the PS wave travel time of various points according to the submarine incident point position to obtain a PS wave non-hyperbolic correction co-imaging trace set;
the hyperbolic time window construction module 830 is configured to determine a velocity spectrum corresponding to the PS-wave non-hyperbolic correction co-imaging gather, pick up a primary reflection PS-wave velocity point from the velocity spectrum, and construct a hyperbolic time window under the condition that a velocity spectrum energy cluster corresponding to the primary reflection PS-wave velocity point is focused;
the multiple attenuation module 840 is configured to perform vector median filtering on the PS-wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, and after the PS-wave non-hyperbolic correction co-imaging gather is subjected to vector median filtering, the developed seafloor multiple is attenuated.
And the speed updating module 850 is configured to update the S-wave offset speed according to the P-wave offset speed and the speed of the picked up primary reflected PS-wave when the velocity spectrum energy cluster corresponding to the primary reflected PS-wave velocity point is not focused.
Wherein, the water depth correction module 810 is specifically configured to:
acquiring the position coordinates, offset distance, sea water depth at the shot point end, P wave offset speed and S wave offset speed of the shot point and the wave detector of the seismic channel;
determining the position coordinates of the submarine incident point according to the offset distance, the sea water depth of the shot point end, the P wave offset speed and the S wave offset speed;
determining a second horizontal distance between the submarine incident point and the PS wave reflection point and a third horizontal distance between the PS wave reflection point and the detection point according to the first horizontal distance between the gun point and the submarine incident point;
calculating a sum of the first horizontal distance, the second horizontal distance and the third horizontal distance;
resetting the position coordinates of the submarine incident point by adopting a dichotomy under the condition that the difference value between the sum value and the offset exceeds a preset threshold value until the calculated difference value between the sum value and the offset is smaller than the preset threshold value;
according to the reset position coordinates of the submarine incident points, the seismic wave travel time in the sea water is solved, and according to the seismic wave travel time in the sea water, the time of each sample point of the seismic channel is corrected.
The non-hyperbolic correction module 820 is specifically configured to:
determining the horizontal distance between the submarine incident point and the detection point according to the submarine incident point position;
and correcting the PS wave travel time of each sample point to a hyperbolic time difference position according to the horizontal distance, the vertical distance of the submarine and PS wave reflection interface and the reflected wave imaging speed.
The multiple attenuation module 840 is specifically configured to:
acquiring a target seismic channel, taking the target seismic channel as a center, and selecting a preset number of seismic channels from the PS wave non-hyperbolic correction co-imaging channel set to form a vector set;
calculating the sum of the distances between the amplitude vector of the target sample point and the amplitude vectors of other sample points except the target sample point in the hyperbolic time window, and outputting the amplitude vector with the minimum sum of the distances between the amplitude vector and the amplitude vector of all the sample points.
The apparatus 800 provided in this embodiment of the present application may perform the methods described in the foregoing method embodiments, and implement the functions and beneficial effects of the methods described in the foregoing method embodiments, which are not described herein again.
Fig. 9 shows a schematic diagram of a hardware structure of a computer device for executing the embodiments of the present application, and referring to the figure, at a hardware level, the computer device includes a processor, and optionally includes an internal bus, a network interface, and a memory. The Memory may include a Memory, such as a Random-Access Memory (RAM), and may further include a non-volatile Memory (non-volatile Memory), such as at least 1 disk Memory. Of course, the computer device may also include hardware required for other services.
The processor, network interface, and memory may be interconnected by an internal bus, which may be an industry standard architecture (Industry Standard Architecture, ISA) bus, a peripheral component interconnect standard (Peripheral Component Interconnect, PCI) bus, or an extended industry standard architecture (Extended Industry Standard Architecture, EISA) bus, among others. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, only one bi-directional arrow is shown in the figure, but not only one bus or one type of bus.
And a memory for storing the program. In particular, the program may include program code including computer-operating instructions. The memory may include memory and non-volatile storage and provide instructions and data to the processor.
The processor reads the corresponding computer program from the nonvolatile memory into the memory and then runs to form a device for locating the target user on a logic level. A processor executing the program stored in the memory, and specifically executing: the embodiments shown in fig. 2-3 disclose the method and implement the functions and advantages of the methods described in the foregoing method embodiments, which are not described in detail herein.
The methods disclosed above in the embodiments of fig. 2-3 of the present application may be implemented in or by a processor. The processor may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software. The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; but also digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present application may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the embodiments of the present application may be embodied directly in hardware, in a decoded processor, or in a combination of hardware and software modules in a decoded processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in a memory, and the processor reads the information in the memory and, in combination with its hardware, performs the steps of the above method.
The computer device may also execute the methods described in the foregoing method embodiments, and implement the functions and beneficial effects of the methods described in the foregoing method embodiments, which are not described herein.
Of course, in addition to software implementation, the computer device of the present application does not exclude other implementation, such as a logic device or a combination of software and hardware, that is, the execution subject of the following process is not limited to each logic unit, but may also be hardware or a logic device.
The embodiments of the present application further provide a computer readable storage medium storing one or more programs, which when executed by an electronic device including a plurality of application programs, cause the electronic device to execute the method disclosed in the embodiments shown in fig. 1-2 and implement the functions and benefits of the methods described in the foregoing method embodiments, which are not described herein.
The computer readable storage medium includes Read-Only Memory (ROM), random access Memory (Random Access Memory RAM), magnetic disk or optical disk, etc.
Further, embodiments of the present application also provide a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, implement the following flow: the embodiments shown in fig. 1-2 disclose the method and implement the functions and advantages of the methods described in the foregoing method embodiments, which are not described herein.
In summary, the foregoing description is only a preferred embodiment of the present application, and is not intended to limit the scope of the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
The system, apparatus, module or unit set forth in the above embodiments may be implemented in particular by a computer chip or entity, or by a product having a certain function. One typical implementation is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises the element.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments.
Claims (10)
1. A method for attenuating a plurality of seafloor waves based on a hyperbolic vector median filter, comprising:
carrying out water depth correction on the seismic channel according to P, S wave migration velocity and submarine topography data to obtain the seismic channel after water depth correction and submarine incident point positions corresponding to various points of the seismic channel;
performing prestack time migration on the amplitudes of all the sample points in the seismic channel after the water depth correction, and performing non-hyperbolic correction on the PS wave travel time of all the sample points according to the submarine incident point positions to obtain a PS wave non-hyperbolic correction co-imaging gather;
determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition of focusing velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points;
and carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, wherein after vector median filtering, the submarine multiples developed on the PS wave non-hyperbolic correction co-imaging gather are attenuated.
2. The method of claim 1, wherein depth correcting the seismic traces from the seismic migration velocity and the seafloor terrain data comprises:
acquiring the position coordinates, offset distance, sea water depth at the shot point end, P wave offset speed and S wave offset speed of the shot point and the wave detector of the seismic channel;
determining the position coordinates of the submarine incident point according to the offset distance, the sea water depth of the shot point end, the P wave offset speed and the S wave offset speed;
determining a second horizontal distance between the submarine incident point and the PS wave reflection point and a third horizontal distance between the PS wave reflection point and the detection point according to the first horizontal distance between the gun point and the submarine incident point;
calculating a sum of the first horizontal distance, the second horizontal distance and the third horizontal distance;
resetting the position coordinates of the submarine incident point by adopting a dichotomy under the condition that the difference value between the sum value and the offset exceeds a preset threshold value until the calculated difference value between the sum value and the offset is smaller than the preset threshold value;
according to the reset position coordinates of the submarine incident points, the seismic wave travel time in the sea water is solved, and according to the seismic wave travel time in the sea water, the time of each sample point of the seismic channel is corrected.
3. The method of claim 1, wherein performing non-hyperbolic correction of PS-wave travel times for each sample based on the seafloor entry point location comprises:
determining the horizontal distance between the submarine incident point and the detection point according to the submarine incident point position;
and correcting the PS wave travel time of each sample point to a hyperbolic time difference position according to the horizontal distance, the vertical distance of the submarine and PS wave reflection interface and the reflected wave imaging speed.
4. The method of claim 1, wherein vector median filtering the PS wave non-hyperbolic corrected co-imaging gather according to the hyperbolic time window comprises:
acquiring a target seismic channel, taking the target seismic channel as a center, and selecting a preset number of seismic channels from the PS wave non-hyperbolic correction co-imaging channel set to form a vector set;
calculating the sum of the distances between the amplitude vector of the target sample point and the amplitude vectors of other sample points except the target sample point in the hyperbolic time window, and outputting the amplitude vector with the minimum sum of the distances between the amplitude vector and the amplitude vector of all the sample points.
5. The method as recited in claim 1, further comprising:
and under the condition that the velocity spectrum energy group corresponding to the primary reflection PS wave velocity point is not focused, updating the S wave offset velocity according to the P wave offset velocity and the velocity of the picked primary reflection PS wave.
6. A hyperbolic vector median filter-based ocean bottom multiple attenuation device, comprising:
the water depth correction module is used for carrying out water depth correction on the seismic channel according to the P, S wave offset speed and the submarine topography data to obtain the seismic channel after the water depth correction and the submarine incident point positions corresponding to various points of the seismic channel;
the non-hyperbolic correction module is used for carrying out pre-stack time migration on the amplitudes of all the points in the seismic channel after the water depth correction, and carrying out non-hyperbolic correction on the PS wave travel time of all the points according to the submarine incident point position to obtain a PS wave non-hyperbolic correction co-imaging gather;
the hyperbolic time window construction module is used for determining a velocity spectrum corresponding to the PS wave non-hyperbolic correction co-imaging gather, picking up primary reflection PS wave velocity points from the velocity spectrum, and constructing a hyperbolic time window under the condition that velocity spectrum energy groups corresponding to the primary reflection PS wave velocity points are focused;
and the multiple attenuation module is used for carrying out vector median filtering on the PS wave non-hyperbolic correction co-imaging gather according to the hyperbolic time window, and attenuating the submarine multiple developed on the PS wave non-hyperbolic correction co-imaging gather after the vector median filtering.
7. The method of claim 6, wherein the water depth correction module is specifically configured to:
acquiring the position coordinates, offset distance, sea water depth at the shot point end, P wave offset speed and S wave offset speed of the shot point and the wave detector of the seismic channel;
determining the position coordinates of the submarine incident point according to the offset distance, the sea water depth of the shot point end, the P wave offset speed and the S wave offset speed;
determining a second horizontal distance between the submarine incident point and the PS wave reflection point and a third horizontal distance between the PS wave reflection point and the detection point according to the first horizontal distance between the gun point and the submarine incident point;
calculating a sum of the first horizontal distance, the second horizontal distance and the third horizontal distance;
resetting the position coordinates of the submarine incident point by adopting a dichotomy under the condition that the difference value between the sum value and the offset exceeds a preset threshold value until the calculated difference value between the sum value and the offset is smaller than the preset threshold value;
according to the reset position coordinates of the submarine incident points, the seismic wave travel time in the sea water is solved, and according to the seismic wave travel time in the sea water, the time of each sample point of the seismic channel is corrected.
8. The method of claim 6, wherein the non-hyperbolic correction module is specifically configured to:
determining the horizontal distance between the submarine incident point and the detection point according to the submarine incident point position;
and correcting the PS wave travel time of each sample point to a hyperbolic time difference position according to the horizontal distance, the vertical distance of the submarine and PS wave reflection interface and the reflected wave imaging speed.
9. The method of claim 6, wherein the multiple attenuation module is specifically configured to:
acquiring a target seismic channel, taking the target seismic channel as a center, and selecting a preset number of seismic channels from the PS wave non-hyperbolic correction co-imaging channel set to form a vector set;
calculating the sum of the distances between the amplitude vector of the target sample point and the amplitude vectors of other sample points except the target sample point in the hyperbolic time window, and outputting the amplitude vector with the minimum sum of the distances between the amplitude vector and the amplitude vector of all the sample points.
10. The method as recited in claim 6, further comprising:
and the speed updating module is used for updating the S-wave offset speed according to the P-wave offset speed and the speed of the picked primary reflection PS wave under the condition that the speed spectrum energy group corresponding to the primary reflection PS wave speed point is not focused.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310098387.8A CN115963565B (en) | 2023-01-20 | 2023-01-20 | Seabed multiple attenuation method and device based on hyperbolic vector median filter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310098387.8A CN115963565B (en) | 2023-01-20 | 2023-01-20 | Seabed multiple attenuation method and device based on hyperbolic vector median filter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115963565A CN115963565A (en) | 2023-04-14 |
| CN115963565B true CN115963565B (en) | 2023-06-16 |
Family
ID=85905038
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310098387.8A Active CN115963565B (en) | 2023-01-20 | 2023-01-20 | Seabed multiple attenuation method and device based on hyperbolic vector median filter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115963565B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005026776A1 (en) * | 2003-09-16 | 2005-03-24 | Geosystem S.R.L. | Wide-offset-range pre-stack depth migration method for seismic exploration |
| CN105974405A (en) * | 2016-05-04 | 2016-09-28 | 哈尔滨工业大学 | Amplitude weighting-based ground penetrating radar back projection imaging method |
| WO2018071628A1 (en) * | 2016-10-13 | 2018-04-19 | Downunder Geosolutions (America) Llc | Method for the attenuation of multiple refelections in shallow water settings |
| CN108828668A (en) * | 2018-03-27 | 2018-11-16 | 中国石油天然气集团有限公司 | A kind of pre-stack time migration data processing method and device |
| WO2020168406A1 (en) * | 2019-02-20 | 2020-08-27 | Eneva S.A. | Method for improving 2d seismic acquisition |
| CN112327362A (en) * | 2020-10-30 | 2021-02-05 | 中国海洋大学 | Submarine multiple prediction and tracking attenuation method in velocity domain |
| CN113093280A (en) * | 2021-04-07 | 2021-07-09 | 青岛海洋地质研究所 | Equal-floating correction method for virtual reflection travel-time cable based on coherent function control |
| CN115144899A (en) * | 2022-06-24 | 2022-10-04 | 中国地质大学（北京） | Rugged seabed OBN elastic wave combined deflection imaging method and device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9477000B2 (en) * | 2013-01-11 | 2016-10-25 | Cgg Services Sa | System and method for the removal of shallow water multiples using a hybrid multi-channel prediction method |
| CN111474584B (en) * | 2020-05-29 | 2023-05-05 | 核工业北京地质研究院 | Focusing superposition imaging method and system based on correlated seismic interference |
-
2023
- 2023-01-20 CN CN202310098387.8A patent/CN115963565B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005026776A1 (en) * | 2003-09-16 | 2005-03-24 | Geosystem S.R.L. | Wide-offset-range pre-stack depth migration method for seismic exploration |
| CN105974405A (en) * | 2016-05-04 | 2016-09-28 | 哈尔滨工业大学 | Amplitude weighting-based ground penetrating radar back projection imaging method |
| WO2018071628A1 (en) * | 2016-10-13 | 2018-04-19 | Downunder Geosolutions (America) Llc | Method for the attenuation of multiple refelections in shallow water settings |
| CN108828668A (en) * | 2018-03-27 | 2018-11-16 | 中国石油天然气集团有限公司 | A kind of pre-stack time migration data processing method and device |
| WO2020168406A1 (en) * | 2019-02-20 | 2020-08-27 | Eneva S.A. | Method for improving 2d seismic acquisition |
| CN112327362A (en) * | 2020-10-30 | 2021-02-05 | 中国海洋大学 | Submarine multiple prediction and tracking attenuation method in velocity domain |
| CN113093280A (en) * | 2021-04-07 | 2021-07-09 | 青岛海洋地质研究所 | Equal-floating correction method for virtual reflection travel-time cable based on coherent function control |
| CN115144899A (en) * | 2022-06-24 | 2022-10-04 | 中国地质大学（北京） | Rugged seabed OBN elastic wave combined deflection imaging method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115963565A (en) | 2023-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3017322B1 (en) | Predicting multiples in survey data | |
| CN107193045B (en) | A kind of seismic data processing technique and device | |
| US20090245022A1 (en) | Dual-wavefield multiple attenuation | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| WO1997011390A2 (en) | Multiple suppression in geophysical data | |
| CN108445532B (en) | A kind of Depth Domain inverse migration method and device | |
| CN115963565B (en) | Seabed multiple attenuation method and device based on hyperbolic vector median filter | |
| MXPA06011579A (en) | Generalized 3d surface multiple prediction. | |
| CA3117329A1 (en) | Seismic random noise attenuation | |
| CN111239832B (en) | 3D3C-VSP imaging processing method, device and equipment | |
| CN117452478A (en) | Reconstruction method and device for seismic compressed sampling data | |
| CN111337973B (en) | Seismic data reconstruction method and system | |
| CN109471173B (en) | Residual static correction method and device | |
| CN111781644B (en) | Linear interference attenuation method and device for shallow and medium stratum seismic data | |
| CN111948714B (en) | Noise reduction method and noise reduction device for uplink wave in submarine node seismic data | |
| CN106597534A (en) | Method and device for separating related multiples of free surfaces of different orders | |
| WO2002065159A1 (en) | Method of seismic imaging using direct travel time computing | |
| CN113031071B (en) | Static correction method and device for long wavelength of earthquake wave | |
| CN113640882B (en) | Acquisition footprint noise removal method, electronic device and computer readable storage medium | |
| CN112698394B (en) | Method and device for determining size of design surface element of two-dimensional observation system | |
| CN112394410B (en) | Ultra-shallow imaging information processing method and device | |
| CN115963543A (en) | Seismic data processing method and device, electronic equipment and storage medium | |
| Yao et al. | Weak Signals Recovery of Downhole DAS with Scale-Weighted Non-Local Selective Attention | |
| CN116027413A (en) | Background velocity field reconstruction method based on reflected wave waveform inversion gradient preprocessing | |
| CN112180448A (en) | Residual static correction method and device and electronic equipment |
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 |