MULTI-MODE DISPERSION ENERGY IMAGING DEVICE AND METHOD FOR A FOUR-COMPONENT MARINE INTERFACE WAVE OF AN OCEAN BOTTOM SEISMOMETER

Abstract

The present invention provides a multi-mode dispersion energy imaging device and method for a four-component marine interface wave of an ocean bottom seismometer, belonging to the technical field of marine seismic exploration. The method includes the following steps: designing an marine interface wave artificial seismic observation system, designing a reasonable observation system according to the geological condition of the operation area to ensure the resolution of the imaging to perform the marine artificial source seismic operation; carrying out the data preprocessing of the seafloor surface wave, and then carrying out the three-component seismometer Scholte wave and the acoustic guided wave dispersion energy imaging, and the one-component hydrophone acoustic guided wave dispersion energy imaging; superposing and normalizing the three-component Scholte wave dispersion energy spectrum and the one-component acoustic guided wave dispersion energy spectrum. The device is implemented based on the method above.

Claims

1. A multi-mode dispersion energy imaging method for a four-component marine interface wave of an ocean bottom seismometer, comprising the following steps: step 101: designing an marine interface wave artificial seismic observation system, wherein the marine interface wave artificial seismic observation system comprises an ocean bottom seismometer deployment point and an artificial seismic source excitation point; carrying out marine seismic data acquisition according to the marine interface wave artificial seismic observation system, wherein the marine seismic data acquisition comprises ocean bottom seismometer deployment, air gun firing, ocean bottom seismometer retrieval, and ocean bottom seismometer original data retrieval; step 102: using excitation time and position navigation information of an air gun excitation seismic source to acquire a common receiving point gather excited by an air gun from an original continuous data sequence recorded by the ocean bottom seismometer while carrying out cross-section position correction and clock offset correction, performing trace equalization and band-pass filtering, exporting a result and saving the result into a SU or SEGY format; step 103: reading data of the common receiving point gather of the air gun excitation seismic source from a three-component seismometer, and performing attitude correction on the three-component seismometer to obtain three-component seismometer gather data in a shot line direction, a vertical shot line direction and a vertical direction; and performing dispersion energy imaging of a Scholte wave in the shot line direction and the vertical direction, respectively, using a phase shifting method to obtain a shot line dispersion energy spectrum of the Scholte wave and a vertical dispersion energy spectrum of the Scholte wave; step 104: reading the data of the common receiving point of the air gun excitation source from a hydrophone component, and using the phase shifting method to analyze dispersion energy of a marine acoustic guided wave of the hydrophone component to obtain a dispersion energy spectrum of the marine acoustic guided wave of the hydrophone component; step 105: superposing the shot line dispersion energy spectrum of the Scholte wave, the vertical dispersion energy spectrum of the Scholte wave and the dispersion energy spectrum of the marine acoustic guided wave to obtain a superposed marine interface wave dispersion energy spectrum, and carrying out normalization processing on the superposed marine interface wave dispersion energy spectrum to obtain a final marine interface wave dispersion energy map.

2. The multi-mode dispersion energy imaging method according to claim 1, wherein, in step 101, the ocean bottom seismometer for ocean bottom seismic observation and marine interface wave collection comprises: a four-component data acquisition system, wherein the four-component data acquisition system comprises the three-component seismometer and a one-component hydrophone; the three-component seismometer is a broadband seismometer, having a minimum frequency range of 1-250 Hz; the one-component hydrophone capable of recording low-frequency acoustic signals is an acoustic hydrophone, having a frequency range of 5-500 Hz; a data sampling rate is 500 SPS, the air gun excitation seismic source requires a total capacity of 3000 INCH.sup.3, and an excitation sub-wave band range width is greater than 60 Hz, wherein, a high frequency of the excitation sub-wave band range width is not higher than 120 Hz, and a low frequency of the excitation sub-wave band range width is not lower than 5 Hz.

3. The multi-mode dispersion energy imaging method according to claim 1, wherein, in step 101, a water depth of the marine interface wave artificial seismic observation system is no more than 100 m, and a maximum offset distance excited by the air gun excitation seismic source is no less than 2500 m; a submerged depth of the air gun is greater than 5 m, and an excitation point distance Δ.sub.x of the air gun is: ${\Delta }_x\le \frac{V_{ag}}{f_{\max }},$ wherein V.sub.ag is a minimum phase velocity of a marine interface wave of a detection area, and ƒ.sub.max is a maximum observable frequency of the marine interface wave, an operation process of the marine seismic data acquisition comprises: deploying the ocean bottom seismometer, after the ocean bottom seismometer sinks freely to an ocean bottom, performing an air gun seismic source excitation operation, and carrying out the ocean bottom seismometer retrieval.

4. The multi-mode dispersion energy imaging method according to claim 1, wherein, in step 103, the attitude correction of the three-component seismometer comprises: correcting a roll angle and a pitch angle; and correcting an azimuth angle; a correction method for the pitch angle and the roll angle is as follows:
V.sub.xyν=R.sub.θR.sub.ϕV.sub.xyz, wherein V.sub.xyν is a three-component seismic data column vector after the pitch angle and roll angle are corrected to a horizontal attitude, V.sub.xyz is an original three-component seismic data column vector, and R.sub.θ and R.sub.ϕ are rotation matrices of the pitch angle and the roll angle respectively, and expressions of the rotation matrices are: $R_{\theta }=\left[\begin{array}{ccc}\cos \theta & 0& -\sin \theta \\ 0& 1& 0\\ \sin \theta & 0& \cos \theta \end{array}\right],R_{\varphi }=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \varphi & -\sin \varphi \\ 0& \sin \varphi & \cos \varphi \end{array}\right],$ wherein θ is the pitch angle, and φ is the roll angle, a correction method for the azimuth angle is as follows:
V.sub.rtν=R.sub.γV.sub.xyν, wherein V.sub.rtν is a column vector after the azimuth angle is corrected, and R.sub.γ is an azimuth rotation matrix as follows: $R_{\gamma }=\left[\begin{array}{ccc}\cos \gamma & \sin \gamma & 0\\ -\sin \gamma & \cos \gamma & 0\\ 0& 0& 1\end{array}\right],$ wherein γ is the azimuth angle obtained by a direct wave polarization analysis.

5. The multi-mode dispersion energy imaging method according to claim 1, wherein, a method of calculating a dispersion energy spectrum of the marine interface wave by the phase shifting method comprises: performing Fourier transform on the common receiving point gather r.sub.i of a vertical component marine interface wave in a time-space domain:
R.sub.i(x.sub.k,ω)=FFT(r.sub.i), wherein i is a serial number of a trace, x.sub.k is an offset distance, and R.sub.i(x.sub.k, ω) is a frequency spectrum; performing normalization on an amplitude of the frequency spectrum to obtain a normalized frequency spectrum $T_i\left(x_k,\omega \right)=\frac{R_i\left(x_k,\omega \right)}{.Math.R_i\left(x_k,\omega \right).Math.},$ wherein |R.sub.i(x.sub.k, ω)| represents a modulo operation; when a phase velocity ν.sub.s and a frequency ω of the marine interface wave are given at a point, acquiring dispersion energy at the point, wherein the dispersion energy is: $D\left(v_s,\omega \right)=\frac{1}{N}\underset{1}{\overset{N}{.Math.}}{e}^{j\left(\frac{x_k\omega }{v_s}\right)}{T}_i\left(i=1,\mathrm{.Math.},N\right),$ wherein j is an imaginary unit, selecting a phase velocity range and a frequency range of the marine interface wave, and cyclically calculating the dispersion energy corresponding to all points in the phase velocity range and the frequency range at a specific phase velocity and a frequency interval to obtain an marine interface wave dispersion energy spectrum of the common receive point gather in the time-space domain, wherein for a four-component Scholte wave and acoustic guided wave dispersion energy spectrum calculation, the phase velocity range is 100 m/s-1800 m/s, a phase velocity interval is 1 m/s, and the frequency range is 0.8-150 Hz.

6. The multi-mode dispersion energy imaging method according to claim 1, wherein, superposition and normalization processing is performed on Scholte wave dispersion energy D.sub.r(ν.sub.s, ω) in the shot line direction, Scholte wave dispersion energy D.sub.z(ν.sub.s, ω) in the vertical direction, and acoustic guided wave dispersion energy D.sub.d(ν.sub.s, ω) of the hydrophone component by using the following equation:
D.sub.s(ν.sub.s,ω)=⅓{D.sub.r(ν.sub.s,ω)+D.sub.z(ν.sub.s,ω)+D.sub.d(ν.sub.s,ω)}, wherein D.sub.s(ν.sub.s, ω) is marine interface wave energy after being superposed and normalized.

7. A multi-mode dispersion energy imaging device for realizing the multi-mode dispersion energy imaging method according to claim 1, wherein, the multi-mode dispersion energy imaging device comprises: an ocean bottom seismometer raw data reading module configured for reading raw data in a binary format recorded by the ocean bottom seismometer; an marine interface wave data preprocessing module configured for converting the raw data format, extracting the common receiving point gather, and performing the attitude correction, the clock offset correction, the trace equalization and the band-pass filtering; an marine interface wave phase shifting method analysis module configured for calculating a dispersion energy spectrum of the common receiving point gather; an marine interface wave dispersion energy superposition module configured for superposing the shot line dispersion energy spectrum of the Scholte wave, the vertical dispersion energy spectrum of the Scholte wave and the dispersion energy spectrum of the acoustic guided wave, and performing normalization; an marine interface wave dispersion energy spectrum storage and exporting module configured for storing and exporting the final marine interface wave dispersion energy map.

8. The multi-mode dispersion energy imaging device according to claim 7, wherein, in step 101, the ocean bottom seismometer for ocean bottom seismic observation and marine interface wave collection comprises: a four-component data acquisition system, wherein the four-component data acquisition system comprises the three-component seismometer and a one-component hydrophone; the three-component seismometer is a broadband seismometer, having a minimum frequency range of 1-250 Hz; the one-component hydrophone capable of recording low-frequency acoustic signals is an acoustic hydrophone, having a frequency range of 5-500 Hz; a data sampling rate is 500 SPS, the air gun excitation seismic source requires a total capacity of 3000 INCH.sup.3, and an excitation sub-wave band range width is greater than 60 Hz, wherein, a high frequency of the excitation sub-wave band range width is not higher than 120 Hz, and a low frequency of the excitation sub-wave band range width is not lower than 5 Hz.

9. The multi-mode dispersion energy imaging device according to claim 7, wherein, in step 101, a water depth of the marine interface wave artificial seismic observation system is no more than 100 m, and a maximum offset distance excited by the air gun excitation seismic source is no less than 2500 m; a submerged depth of the air gun is greater than 5 m, and an excitation point distance Δ.sub.x of the air gun is: ${\Delta }_x\le \frac{V_{ag}}{f_{\max }},$ wherein V.sub.ag is a minimum phase velocity of a marine interface wave of a detection area, and ƒ.sub.max is a maximum observable frequency of the marine interface wave, an operation process of the marine seismic data acquisition comprises: deploying the ocean bottom seismometer, after the ocean bottom seismometer sinks freely to an ocean bottom, performing an air gun seismic source excitation operation, and carrying out the ocean bottom seismometer retrieval.

10. The multi-mode dispersion energy imaging device according to claim 7, wherein, in step 103, the attitude correction of the three-component seismometer comprises: correcting a roll angle and a pitch angle; and correcting an azimuth angle; a correction method for the pitch angle and the roll angle is as follows:
V.sub.xyν=R.sub.θR.sub.ϕV.sub.xyz, wherein V.sub.xyν is a three-component seismic data column vector after the pitch angle and roll angle are corrected to a horizontal attitude, V.sub.xyz is an original three-component seismic data column vector, and R.sub.θ and R.sub.ϕ are rotation matrices of the pitch angle and the roll angle respectively, and expressions of the rotation matrices are: $R_{\theta }=\left[\begin{array}{ccc}\cos \theta & 0& -\sin \theta \\ 0& 1& 0\\ \sin \theta & 0& \cos \theta \end{array}\right],$ $R_{\varphi }=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \varphi & -\sin \varphi \\ 0& \sin \varphi & \cos \varphi \end{array}\right],$ wherein θ is the pitch angle, and φ is the roll angle, a correction method for the azimuth angle is as follows:
V.sub.rtν=R.sub.γV.sub.xyν, wherein V.sub.rtν is a column vector after the azimuth angle is corrected, and R.sub.γ is an azimuth rotation matrix as follows: $R_{\gamma }=\left[\begin{array}{ccc}\cos \gamma & \sin \gamma & 0\\ -\sin \gamma & \cos \gamma & 0\\ 0& 0& 1\end{array}\right],$ wherein γ is the azimuth angle obtained by a direct wave polarization analysis.

11. The multi-mode dispersion energy imaging device according to claim 7, wherein, a method of calculating a dispersion energy spectrum of the marine interface wave by the phase shifting method comprises: performing Fourier transform on the common receiving point gather r.sub.i of a vertical component marine interface wave in a time-space domain:
R.sub.i(x.sub.k,ω)=FFT(r.sub.i), wherein i is a serial number of a trace, x.sub.k is an offset distance, and R.sub.i((x.sub.k,ω) is a frequency spectrum; performing normalization on an amplitude of the frequency spectrum to obtain a normalized frequency spectrum $T_i\left(x_k,\omega \right)=\frac{R_i\left(x_k,\omega \right)}{.Math.R_i\left(x_k,\omega \right).Math.},$ wherein |R.sub.i(x.sub.k,ω)| represents a modulo operation; when a phase velocity ν.sub.s and a frequency ω of the marine interface wave are given at a point, acquiring dispersion energy at the point, wherein the dispersion energy is: $D\left(v_s,\omega \right)=\frac{1}{N}\underset{1}{\overset{N}{.Math.}}{e}^{j\left(\frac{x_k\omega }{v_s}\right)}{T}_i\left(i=1,\mathrm{.Math.},N\right),$ wherein j is an imaginary unit, selecting a phase velocity range and a frequency range of the marine interface wave, and cyclically calculating the dispersion energy corresponding to all points in the phase velocity range and the frequency range at a specific phase velocity and a frequency interval to obtain an marine interface wave dispersion energy spectrum of the common receive point gather in the time-space domain, wherein for a four-component Scholte wave and acoustic guided wave dispersion energy spectrum calculation, the phase velocity range is 100 m/s-1800 m/s, a phase velocity interval is 1 m/s, and the frequency range is 0.8-150 Hz.

12. The multi-mode dispersion energy imaging device according to claim 7, wherein, superposition and normalization processing is performed on Scholte wave dispersion energy D.sub.r(ν.sub.s, ω) in the shot line direction, Scholte wave dispersion energy D.sub.z(ν.sub.s, ω) in the vertical direction, and acoustic guided wave dispersion energy D.sub.d(ν.sub.s, ω) of the hydrophone component by using the following equation:
D.sub.s(ν.sub.s,ω)=⅓{D.sub.r(ν.sub.s,ω)+D.sub.d(ν.sub.s,ω)+D.sub.d(ν.sub.s,ω)}, wherein D.sub.s(ν.sub.s, ω) is marine interface wave energy after being superposed and normalized.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] FIG. 1 is a flowchart of the method of the present invention;

[0021] FIG. 2 is a structural block diagram of the device of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0022] The specific embodiments of the present invention will be described below with reference to the drawings:

[0023] As shown in FIG. 1, the present invention provides a multi-mode dispersion energy imaging method for a four-component marine interface wave of an ocean bottom seismometer, including the following steps:

[0024] step 101: designing an marine interface wave artificial seismic observation system, which includes an ocean bottom seismometer deployment point and an artificial seismic source excitation point; carrying out marine seismic data acquisition according to the designed marine interface wave artificial seismic observation system, which includes ocean bottom seismometer deployment, air gun firing, ocean bottom seismometer retrieval, and ocean bottom seismometer original data retrieval;

[0025] step 102: using the excitation time and position navigation information of the artificial seismic source to acquire a common receiving point gather excited by an air gun from an original continuous data sequence recorded by an ocean bottom seismometer while carrying out cross-section position correction and clock offset correction, performing trace equalization and band-pass filtering, exporting a result and saving the result into SU or SEGY format;

[0026] step 103: reading the data of the common receiving point gather of the artificial seismic source from a three-component seismometer, and performing attitude correction on the three-component seismometer to obtain three-component seismometer gather data in a shot line direction, a vertical shot line direction and a vertical direction; and performing dispersion energy imaging of the Scholte wave in the shot line direction and the vertical direction respectively using a phase shifting method to obtain a shot line dispersion energy spectrum and a vertical dispersion energy spectrum, respectively, of the Scholte wave;

[0027] step 104: reading the data of the common receiving point of the artificial seismic source from a hydrophone component, and using the phase shifting method to analyze the dispersion energy of the marine acoustic guided wave of the hydrophone component to obtain the dispersion energy spectrum of the acoustic guided wave of the hydrophone component;

[0028] step 105: superposing the shot line dispersion energy spectrum of the Scholte wave, the vertical dispersion energy spectrum of the Scholte wave and the dispersion energy spectrum of the acoustic guided wave, and carrying out normalization processing on the superposed marine interface wave dispersion energy spectrum to obtain a final marine interface wave dispersion energy map.

[0029] In the above step 101, the ocean bottom seismometer for ocean bottom seismic observation and marine interface waves collection has the specification parameter requirements as follows: (1) a four-component data acquisition system, including a three-component seismometer and a one-component hydrophone; (2) the three-component seismometer is a broadband seismometer, having the minimum frequency band range of 1-250 Hz; (3) a frequency band range of an acoustic hydrophone capable of recording low-frequency acoustic signals is 5-500 Hz; (4) a data sampling rate meets 500 SPS, preferably 1000 SPS. As for the requirement for an air gun excitation seismic source, the total capacity of the air gun must meet the requirement of 3000 inch.sup.3, and an excitation sub-wave band range width is greater than 60 Hz, wherein, the high frequency is not higher than 120 Hz, and the low frequency is not lower than 5 Hz.

[0030] In the above step 101, the observation system is designed to have a water depth of no more than 100 m. The maximum offset distance excited by the air gun excitation seismic source is no less than 2500 m because the large offset distance increases the resolution of dispersion energy imaging. The submerged depth of the air gun is greater than 5 m, and it is recommended that the submerged depth of the air gun should be greater than 10 m in areas with greater water depth. An excitation point distance Δ.sub.x of the air gun needs to satisfy the following formula:

[00001]${\Delta }_x\le \frac{V_{\mathrm{ag}}}{f_{\max }},$

[0031] wherein V.sub.ag is the minimum phase velocity of the marine interface wave of the detection area, and f.sub.max is the maximum observable frequency of the marine interface wave. In the operation process, the ocean bottom seismometer is deployed first, and after the ocean bottom seismometer sinks freely to the ocean bottom, the air gun seismic source excitation operation is performed, and finally, the ocean bottom seismometer retrieval is carried out.

[0032] In the above step 103, the attitude correction of the three-component seismometer includes: correcting a roll angle and a pitch angle to ensure that the vertical component is perpendicular to the normal direction of the sea surface; and correcting an azimuth angle, so that the X component points to the positive direction of the shot line and the Y component is perpendicular to the positive direction of the shot line. The method for correcting the pitch angle and roll angle is as follows:

V.sub.xyν=R.sub.θR.sub.ϕV.sub.xyz,

[0033] wherein V.sub.xyν is a three-component seismic data column vector after the pitch angle and roll angle are corrected to a horizontal attitude, V.sub.xyz is an original three-component seismic data column vector, and R.sub.θ and R.sub.ϕ are the rotation matrices of the pitch angle and the roll angle respectively, and the expressions thereof are as follows:

[00002]$R_{\theta }=\left[\begin{array}{ccc}\cos \theta & 0& -\sin \theta \\ 0& 1& 0\\ \sin \theta & 0& \cos \theta \end{array}\right],R_{\varphi }=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \varphi & -\sin \varphi \\ 0& \sin \varphi & \cos \varphi \end{array}\right],$

[0034] wherein θ is the pitch angle, and φ is the roll angle.

[0035] A correction method for the azimuth angle is:

V.sub.rtν=R.sub.γV.sub.xyν,

[0036] wherein V.sub.rtν is the column vector after the azimuth angle is corrected, and R.sub.γ is the azimuth rotation matrix as follows

[00003]$R_{\gamma }=\left[\begin{array}{ccc}\cos \gamma & \sin \gamma & 0\\ -\sin \gamma & \cos \gamma & 0\\ 0& 0& 1\end{array}\right],$

[0037] wherein γ is the azimuth angle obtained by direct wave polarization analysis.

[0038] In the above steps 103 and 104, the method of calculating the dispersion energy spectrum of the marine interface wave by using the phase shifting method includes;

[0039] performing the Fourier transform on the common receiving point gather r.sub.i of the vertical component marine interface wave in the time-space domain:

R.sub.i(x.sub.k,ω)=FFT(r.sub.i),

[0040] wherein i is the serial number of the trace, x.sub.k is the offset distance, and R.sub.i(x.sub.k, ω) is the frequency spectrum;

[0041] then performing the normalization on the amplitude of the frequency spectrum to obtain the normalized frequency spectrum T.sub.i(x.sub.k, ω):

[00004]$T_i\left(x_k,\omega \right)=\frac{R_i\left(x_k,\omega \right)}{.Math.R_i\left(x_k,\omega \right).Math.},$

[0042] wherein |R.sub.i(x.sub.k, ω)| represents a modulo operation;

[0043] when the phase velocity ν.sub.s and frequency ω of the marine interface wave are given, acquiring the dispersion energy at the point as:

[00005]$D\left(v_s,\omega \right)=\frac{1}{N}\underset{1}{\overset{N}{.Math.}}{e}^{j\left(\frac{x_k\omega }{v_s}\right)}{T}_i\left(i=1,\mathrm{.Math.},N\right),$

[0044] wherein j is an imaginary unit,

[0045] selecting a certain phase velocity range and a frequency range of the marine interface wave, and cyclically calculating the dispersion energy corresponding to all points in the range at a specific phase velocity and a frequency interval to obtain the marine interface wave dispersion energy spectrum of the common receive point gather in the time-space domain.

[0046] For the four-component Scholte wave and acoustic guided wave dispersion energy spectrum calculation, the phase velocity range is preferably 100 m/s-1800 m/s, the phase velocity interval is preferably 1 m/s, and the frequency range is preferably 0.8-150 Hz.

[0047] In the above step 105, the superposition and normalization processing is performed on the Scholte wave dispersion energy D.sub.r(ν.sub.s, ω) in the shot line direction, the Scholte wave dispersion energy D.sub.z(ν.sub.s, ω) in the vertical direction, and the acoustic guided wave dispersion energy D.sub.d(ν.sub.s, ω) of the hydrophone by using the following equation:

D.sub.s(ν.sub.s,ω)=⅓{D.sub.r(ν.sub.s,ω)+D.sub.z(ν.sub.s,ω)+D.sub.d(ν.sub.s,ω)},

[0048] wherein D.sub.s(ν.sub.s, ω) is the marine interface wave energy after being superposed and normalized.

[0049] Based on the same inventive concept, the present invention also provides a multi-mode dispersion energy imaging device for a four-component marine interface wave of an ocean bottom seismometer.

[0050] FIG. 2 is a block diagram of a multi-mode dispersion energy imaging device for a four-component marine interface wave of an ocean bottom seismometer of the present invention, including the ocean bottom seismometer raw data reading module 201, the marine interface wave data preprocessing module 202, the marine interface wave phase shifting method analysis module 203, the marine interface wave dispersion energy superposition module 204, and the marine interface wave dispersion energy spectrum storage and exporting module 205. The structure of the device is illustrated as follows.

[0051] The ocean bottom seismometer raw data reading module 201 is configured for reading the raw artificial source seismic data recorded by the ocean bottom seismometer.

[0052] The marine interface wave data preprocessing module 202 is configured for converting the raw data format, extracting the common receiving point gather, and performing position correction, clock offset correction, trace equalization, band-pass filtering and attitude correction.

[0053] The marine interface wave phase shift method analysis module 203 is configured for calculating the dispersion energy spectrum of the common receiving point gather.

[0054] The marine interface wave dispersion energy superposition module 204 is configured for superposing the Scholte wave and the acoustic guided wave dispersion energy spectrum, and performing normalization.

[0055] The marine interface wave dispersion energy spectrum storage and exporting module 205 is configured for storing and exporting the normalized marine interface wave dispersion spectrum.

[0056] In the above module 202, for a attitude correction of the three-component seismometer, the correction includes: correcting the roll angle and the pitch angle are corrected to ensure that the vertical component is perpendicular to the normal direction of the sea surface; and correcting the azimuth angle so that the X component points to the positive direction of the shot line and the Y component is perpendicular to the positive direction of the shot line.

[0057] The method for correcting the pitch angle and the roll angle is as follows:

V.sub.xyν=R.sub.θR.sub.ϕV.sub.xyz,

[0058] wherein V.sub.xyν is the three-component seismic data column vector after the pitch angle and the roll angle are corrected to the horizontal attitude, V.sub.xyz is the original three-component seismic data column vector, and R.sub.θ and R.sub.ϕ are the rotation matrices of the pitch angle and the roll angle respectively, and the expressions thereof are as follows:

[00006]$R_{\theta }=\left[\begin{array}{ccc}\cos \theta & 0& -\sin \theta \\ 0& 1& 0\\ \sin \theta & 0& \cos \theta \end{array}\right],R_{\varphi }=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \varphi & -\sin \varphi \\ 0& \sin \varphi & \cos \varphi \end{array}\right],$

wherein θ is the pitch angle, and φ is the roll angle.
A correction method for the azimuth angle is:

V.sub.rtν=R.sub.γV.sub.xyν.

[0059] wherein V.sub.rtν is the column vector after the azimuth angle is corrected, and R.sub.γ is the azimuth rotation matrix as follows

[00007]$R_{\gamma }=\left[\begin{array}{ccc}\cos \gamma & \sin \gamma & 0\\ -\sin \gamma & \cos \gamma & 0\\ 0& 0& 1\end{array}\right],$

[0060] wherein γ is the azimuth angle obtained by direct wave polarization analysis.

[0061] In the above module 203, the method of calculating the dispersion energy spectrum of the marine interface wave by the phase shifting method includes:

[0062] performing Fourier transform on the common receiving point gather r.sub.i of the vertical component marine interface wave in the time-space domain:

R.sub.i(x.sub.k,ω)=FFT(r.sub.i),

[0063] wherein i is the serial number of the channel, x.sub.k is the offset distance, and R.sub.i((x.sub.k, ω) is the frequency spectrum;

[0064] then performing the normalization on the amplitude of the frequency spectrum to obtain the normalized frequency spectrum

[00008]$T_i\left(x_k,\omega \right)=\frac{R_i\left(x_k,\omega \right)}{.Math.R_i\left(x_k,\omega \right).Math.},$

[0065] wherein |R.sub.i(x.sub.k, ω)| represents a modulo operation;

[0066] when the phase velocity ν.sub.s and frequency ω of the marine interface wave are given, acquiring the dispersion energy at the point as:

[00009]$D\left(v_s,\omega \right)=\frac{1}{N}\underset{1}{\overset{N}{.Math.}}{e}^{j\left(\frac{x_k\omega }{v_s}\right)}{T}_i\left(i=1,\mathrm{.Math.},N\right),$

[0067] wherein j is an imaginary unit,

[0068] selecting a certain phase velocity range and a frequency range of the marine interface wave, and cyclically calculating the dispersion energy corresponding to all points in the range at a specific phase velocity and a frequency interval to obtain the marine interface wave dispersion energy spectrum of the common receive point gather in the time-space domain.

[0069] For the four-component Scholte wave and acoustic guided wave dispersion energy spectrum calculation, the phase velocity range is preferably 100 m/s-1800 m/s, the phase velocity interval is preferably 1 m/s, and the frequency range is preferably 0.8-150 Hz.

[0070] In the above module 205, superposition and normalization processing is performed on the Scholte wave dispersion energy D.sub.r(ν.sub.s, ω) in the shot line direction, the Scholte wave dispersion energy D.sub.z(ν.sub.s, ω) in the vertical direction, and the acoustic guided wave dispersion energy D.sub.d(ν.sub.s, ω) of the hydrophone by using the following equation:

D.sub.s(ν.sub.s,ω)=⅓{D.sub.r(ν.sub.s,ω)+D.sub.z(ν.sub.s,ω)+D.sub.d(ν.sub.s,ω)},

wherein D.sub.s(ν.sub.s, ω) is the marine interface wave energy after being superposed and normalized.

[0071] The present invention achieves the following technical effects: an ocean bottom seismometer four-component marine interface wave multi-mode dispersion energy imaging method and device realize data acquisition of the marine interface wave of the artificial source by using artificial seismic source excitation and an ocean bottom seismometer, and further realize multi-mode and multi-type marine interface wave dispersion energy spectrum comprehensive imaging by combining the dispersion energy spectrum imaging of the Scholte wave, the dispersion energy spectrum imaging of the acoustic guided wave, and the superposition of two types of the marine interface waves dispersion energy spectrum. This technical method is a basis for the extraction of high-mode dispersion curves in the marine interface wave analysis. The extreme value positions in the results map of the dispersion energy spectrum of this technology are utilized to extract the multi-mode marine interface wave dispersion curves, which provides more constraints for the inversion analysis on marine interface wave frequency dispersion curves and greatly promotes the development of marine interface wave detection technology.

[0072] The above are the preferred embodiments of the present invention. It should be pointed out that without departing from the principles of the present invention, those skilled in the art can make several improvements and modifications, and these improvements and modifications should also fall within the protection scope of the present invention.