48 results about "Submarine earthquake" patented technology

The invention relates to the technical field of earthquake prediction, in particular to a submarine seismic oscillation prediction method and device based on earthquake exploration, and the method comprises the following steps: selecting a monitoring region, collecting the submarine topographic data of the monitoring region, and drawing a topographic map; selecting a plurality of monitoring points according to the topographic map; geomagnetic field information, electric wave information, micro-seismic activity information and seabed deformation information of each monitoring point are collected and recorded; converting the collected information into acoustic signals and transmitting the acoustic signals to a sea surface buoy, and converting the acoustic signals into electric signals by the sea surface buoy and transmitting the electric signals to a ground data center; by optimizing site selection of monitoring points, integrating geomagnetic-electric wave-microseismic-deformation multi-source data, realizing high-efficiency data transmission by using an acoustic-electric staged communication chain, and finally realizing early recognition of earthquake precursor in combination with threshold analysis, the method can significantly improve the coverage range, data precision and early warning timeliness of seafloor earthquake monitoring, and has good application prospects. And a technical support is provided for marine geological disaster prevention and control.

One or more embodiments of the invention provide a Scholte wave frequency dispersion curve adaptive inversion method and a related device, and the method comprises the steps: obtaining and processing the original seismic data of a target seabed, and obtaining a Scholte wave frequency dispersion curve; constructing an initial shear wave velocity model through linear inversion based on a Schlte wave frequency dispersion curve; and performing fitting optimization on the Schlte wave frequency dispersion curve based on the initial shear wave velocity model to obtain a two-dimensional shear wave velocity profile. According to the technical scheme, the problem that a traditional linear inversion method depends on an initial model is solved, the calculation efficiency is remarkably improved, the problem of local extremum is avoided, and an automatic high-precision inversion solution suitable for batch data processing is provided for submarine seismic exploration.

The invention provides a submarine earthquake optical fiber monitoring data processing method and device, equipment and a medium, which can be applied to the technical field of communication. The method comprises the following steps: acquiring original optical signal data generated by a submarine seismic signal; preprocessing the original optical signal data to obtain initial monitoring data; performing seismic event identification on the initial monitoring data to obtain seismic event data; and compressing the seismic event data to obtain compressed monitoring data. According to the method provided by the invention, a series of data processing such as preprocessing, filtering and denoising, seismic event identification and compression processing are carried out on the optical fiber monitoring data in sequence, so that the efficiency and precision of accessing the optical fiber monitoring data to the seismic monitoring platform are improved.

The invention provides a seabed free site accurate calibration method combining HVSR and 2D-GIT. According to the method, on the basis of submarine seismological observation data, HVSR and 2D-GIT are adopted to obtain submarine site transfer functions respectively, and a submarine free site is accurately calibrated through difference analysis of the site transfer functions; comprising the steps of seabed seismic data acquisition and preprocessing, HVSR spectrum ratio extraction, two-dimensional earth crust attenuation model modeling, 2D-GIT inversion and data fusion, conjoint analysis and environmental parameter correction, site type evaluation and dynamic parameter calibration. The characteristic that the HVSR technology is obviously influenced by local landforms when the horizontal transfer function of the seabed site is evaluated and the characteristic that the 2D-GIT technology can accurately evaluate the horizontal transfer function of each type of seabed site are jointly utilized, and analysis results of the two methods are compared. The method achieves the evaluation of the type of the seabed site and the more accurate and reliable positioning and characteristic evaluation of the free seabed site, effectively makes up for the limitation of a single technology, and provides a comprehensive and accurate type of the seabed site.

The invention provides a seismic wave field analytical simulation method, which comprises the following steps of: constructing a seabed interface analytical equation according to a displacement potential function and a boundary condition of reflected / transmitted waves; obtaining a first function relationship; obtaining a second function relationship; converting the seabed interface analytical equation by using the first function relationship and the second function relationship to obtain a function relationship between the reflection / transmission coefficient and the influence factor, and determining the function relationship as a third function relationship; a plurality of influence factors are obtained, and reflection / transmission coefficients corresponding to the influence factors are obtained based on the third function relation; and obtaining seabed seismic wave field data through forward modeling based on the influence factors and the reflection / transmission coefficients which are in one-to-one correspondence. According to the seismic wave field analysis simulation method, the obtained seabed seismic wave field is highly similar to real data, and a theoretical model can be provided for prediction and attenuation of seabed multiples and ghost waves.

The invention relates to the technical field of underwater monitoring, and provides a seabed earthquake observation base station which comprises a seabed base, a seismograph, a floating ball, a locking piece and a protective cylinder, the protective cylinder is fixedly installed on the seabed base, the floating ball is connected with the seismograph through a tying rope, the seismograph is installed in the protective cylinder in a pluggable mode, and the protective cylinder is provided with a lock hole; the seismograph has a first state and a second state in the laying process; in the first state, the floating ball is fixed to the seismograph, the locking piece penetrates through the lock hole, and the seismograph is locked in the protective cylinder. And in the second state, the locking piece retreats from the protective cylinder, the seismograph is pulled out of the protective cylinder under the action of external force, and the cable is released so that the floating ball pulls the seismograph to sink into the landfill pit of the seabed. The seismograph can stably sink into the landfill pit of the seabed, the inclination angle of the seismograph is small, the laying difficulty of the whole seabed seismological observation base station is small, and the monitoring effectiveness of the seismograph is improved.

The invention discloses a deep residual network-based submarine seismic signal denoising enhancement method. The method comprises the following steps of decomposing an original signal into a plurality of sub-band signals through a wavelet packet; performing feature extraction on each sub-band signal by using a deep residual convolutional network in combination with a frequency domain attention mechanism; network hyper-parameters are optimized by using an improved marine predator algorithm, and each sub-band feature map is fused; enhancing the signal through short-time Fourier transform and a dynamic gating threshold; periodically, adaptively and finely adjusting network parameters; and adjusting the network according to the signal-to-noise ratio gain, the weak seismic event detection success rate and the positioning error index feedback. According to the invention, the signal-to-noise ratio of the seismic signal and the event detection precision in the complex marine environment are obviously improved.

The application discloses a method and device for extracting Scholte wave attenuation coefficients, a system and a storage medium, and comprises the following steps: collecting seabed seismic data of an observation area; obtaining a pretreated common-receiver gather based on the seabed seismic data of the observation area; based on the CRG gather data in the SU format pretreated by the seabed seismograph, using high-resolution forward and inverse transformation, separating and extracting the base-order CRG gather with energy focusing and single order from the CRG gather coupled and stacked by multiple-order Scholte waves; based on the base-order CRG seismic gather, calculating the attenuation coefficients of any two gathers, and after the cyclic calculation of all gather pairs is completed, the average attenuation coefficient model corresponding to the gather pair is obtained. By using the technical scheme of the application, the Scholte wave dispersion order can be accurately separated, and the base-order attenuation coefficient can be specially extracted.

The application relates to the technical field of underwater detection, in particular to an underwater vehicle autonomous detection system and method based on seabed seismic wave signals. The system comprises a three-axis seismic wave signal acquisition module, a signal superposition amplification module, a signal filtering module, a signal envelope extraction module, an environmental noise integration module and a threshold comparison determination module. The acquisition module is orthogonally arranged to receive three-axis seismic wave signals and convert the signals into charge signals, the superposition amplification module completes signal conversion and superposition, the filtering module suppresses noise, the envelope extraction module extracts a detection value, the noise integration module obtains a noise integral value, the comparison module compares the two values, and when the detection value is greater than the integral value, an underwater vehicle existing signal is output. The application uses seabed seismic waves as a detection source, has strong anti-interference performance, works in a passive mode, is concealed and has low power consumption, can accurately capture weak seismic wave signals of the underwater vehicle, realizes autonomous detection of the underwater vehicle, and is suitable for long-term unattended deployment on the seabed.

The application belongs to the technical field of seismic data acquisition, and specifically discloses a kind of underwater seismic data continuous acquisition method, system, equipment and medium, the system includes the main control module built-in watchdog mechanism, sensor, acquisition auxiliary module, data storage module and power supply module;The method comprises: the main control module starts data acquisition task, sensor collects seabed seismic data, and is stored to data storage module;During data acquisition, the main control module sends the watchdog to send dog feeding signal regularly, if not on time to send dog feeding signal, then the watchdog forces to restart the main control module, and the main control module restarts data acquisition task through acquisition auxiliary module.The application solves the problem of data interruption and time loss caused by program exception of seabed node, ensures the continuity and integrity of marine seismic data acquisition, and improves the reliability and economic benefit of marine oil and gas exploration operation.The application is suitable for continuous acquisition of underwater seismic data.

The present application provides a method, device, equipment and medium for determining the optimal earth radial parameters. Among them, the original observation data of a number of seafloor stations in the epicentral area of submarine earthquakes are obtained; the seafloor station coordinate position solution data and the observed values of the displacement time series are determined based on the original observation data of each seafloor station; a Burgers viscoelastic earth model is constructed, and the reset dislocation Love number and the complex Green's function are determined based on the Burgers viscoelastic earth model; based on the reset dislocation Love number and the complex Green's function, the calculated values of the displacement time series are determined according to the observed values of the displacement time series; the optimal transient viscosity parameter and the optimal steady-state viscosity parameter are determined based on the observed values and the calculated values of the displacement time series, and the transient viscosity parameter and the optimal steady-state viscosity parameter are jointly composed of the optimal earth radial parameters. By adopting the above scheme, the earth radial parameters can be accurately determined to ensure the correct understanding and accurate prediction of the earthquake process.

The present application provides an ocean seismic node acquisition device, belonging to the technical field of marine geophysical exploration. Inside the pressure-resistant inner housing of the device, there are a three-component geophone, a power supply and a circuit unit. Among them, the circuit unit is used to control the piezoelectric geophone to detect the seismic longitudinal wave pressure in seawater, control the three-component geophone to pick up the submarine seismic signals, and control the acoustic transducer for submarine positioning. The power supply is used to provide electrical energy for other components in the device. After the device is placed on the seabed, the device can independently perform detection work as a single node. After a certain period of time, the data stored in the acquisition circuit in the device can be exported through an external socket. Through this device, seismic data corresponding to 4 acquisition channels measured by the piezoelectric geophone and the three-component geophone can be obtained. Since there is no need to use a transmission cable to transmit data to the host system in real time, the cost is relatively low.

A system and method for deploying and retrieving a plurality of ocean bottom seismic nodes to and from the seabed by using a floating platform. The system comprises an autonomous underwater vehicle (AUV) (101) coupled to a node skid (131) that is configured to handle a plurality of ocean bottom seismic nodes. The AUV and coupled skid may be lowered to and raised from the seabed and a platform in a garage. The AUV may comprise one or more fuel cells in addition to or as a replacement of traditional rechargeable batteries. The system may utilize an unmanned underwater vehicle (UUV) or unmanned surface vehicle (USV) to monitor and track AUV deployment and to provide acoustic communications between the surface platform and the AUV.

The application provides a deep-sea T-phase calculation method and system based on a seafloor seismic source mechanism, and the method comprises the following steps: determining calculation physical quantities in an elastic multilayer ocean crust medium and a fluid deep-sea sound channel for a deep-sea seafloor seismic T-phase propagation process, and establishing an acoustic wave-seismic wave combined elliptical differential equation in a two-dimensional cylindrical coordinate system; the elliptical differential equation is converted into an acoustic wave-seismic wave combined parabolic differential equation, a horizontal recursive equation in a root form of the acoustic wave-seismic wave combined parabolic differential equation is given, an initial field horizontal recursive equation in the root form suitable for a seafloor seismic event is established, the horizontal recursive equation in the root form is converted into an N p step loop recursive equation, numerical discrete processing is performed on operator action in the loop recursive equation, a representation form of calculation quantities of virtual points in a discrete process is given, and a seismic wave field in an ocean crust medium and an acoustic wave field in a deep-sea sound channel are obtained according to the loop recursive equation.

The embodiment of the invention provides a multi-point submarine seismic oscillation generation method and related equipment. The time domain non-stationarity of seismic oscillation and the spatial variability of complex terrain can be considered when multi-point seismic oscillation is synthesized. The method comprises the following steps: generating one-dimensional seismic waveform data according to three-direction seismic waveform data, random noise and a target single-point random seismic oscillation model; splicing the one-dimensional seismic waveform data with random data satisfying Dirichlet distribution to obtain spliced multipoint seismic oscillation data; according to the size of the semi-space uniform isotropic geologic body, the position of the fault plane and the size of the fault plane, establishing a three-dimensional site model containing the fault; according to the fault plane parameters corresponding to different moment magnitudes and the three-dimensional site model, performing simulation analysis of seismic oscillation through a spectral element method to obtain multi-point seismic oscillation data of each measuring point in the measuring point set; and inputting the spliced multi-point seismic oscillation data and the multi-point seismic oscillation data into a multi-point seismic oscillation data generation model to obtain multi-point seismic oscillation.

Embodiments of the present application provide a kind of multi-point seabed seismic generation method and related equipment, can be synthesized multi-point seismic time-domain non-stationary and the spatial variability of complex topography when considering seismic.The method comprises: according to three direction seismic waveform data, random noise and target single-point random seismic model generation one-dimensional seismic waveform data;One-dimensional seismic waveform data is spliced with random data satisfying Dirichlet distribution, to obtain spliced multi-point seismic data;According to the size of semi-space homogeneous isotropic geologic body, the position of fault plane and the size of the fault plane, establish three-dimensional site model containing fault;According to the fault plane parameters corresponding to different moment magnitude and three-dimensional site model, the simulation analysis of seismic is carried out by spectral element method, to obtain the multi-point seismic data of each measuring point in measuring point set;Spliced multi-point seismic data and multi-point seismic data are input into multi-point seismic data generation model, to obtain multi-point seismic.

The invention relates to the technical field of underwater monitoring, and provides a submarine seismograph base station and system, the submarine seismograph base station comprises a seabed base, a seismograph, a floating ball, a locking piece and a protective cylinder, the protective cylinder is fixedly installed on the seabed base, the floating ball is connected with the seismograph through a tether, the seismograph is installed in the protective cylinder in a pluggable mode, and the protective cylinder is provided with a lock hole; the seismograph has a first state and a second state in the laying process; in the first state, the floating ball is fixed to the seismograph, the locking piece penetrates through the lock hole, and the seismograph is locked in the protective cylinder. And in the second state, the locking piece retreats from the protective cylinder, the seismograph is pulled out of the protective cylinder under the action of external force, and the cable is released so that the floating ball pulls the seismograph to sink into the landfill pit of the seabed. The seismograph can stably sink into the landfill pit of the seabed, the inclination angle of the seismograph is small, the laying difficulty of the whole seabed seismological observation base station is small, and the monitoring effectiveness of the seismograph is improved.

The invention discloses an efficient fluctuating seabed seismic wave numerical simulation method and a related device, and belongs to the technical field of seismology. The method comprises the following steps: constructing a push-up interface label-vector coupling wave equation; constructing a seismic wave operator recursion process of the push-up fluctuation interface based on the push-up interface mark-vector coupling wave equation; performing model reconstruction and interpolation through the original elastic parameter model to obtain a Kaiser window sinchinger function filtering speed model; and substituting the Kaiser window sinchinger function filtering speed model into the standard-vector coupling wave equation, and carrying out numerical simulation and outputting a seismic wave field in combination with a seismic wave operator recursion process of an up-pushing fluctuating interface. The method successfully solves the problems of stepped scattering and inaccurate travel time in seismic wave numerical simulation under the fluctuating seabed condition, does not need to greatly increase extra calculation amount, and can accurately simulate the propagation process of seismic waves in the marine environment.

The invention relates to the field of seismic data analysis, in particular to a seismic data noise suppression method and system based on deep learning. The method comprises the following steps: acquiring original seafloor seismic data of a target area and marine environment parameters which are in space-time correlation with the original seafloor seismic data; the marine environment parameters are mapped to generate a space-time continuous physical guidance noise field, the generation process of the physical guidance noise field is constrained by physical consistency, and the physical guidance noise field has low-frequency, coherent and smooth physical attributes of flow-induced noise; and taking the physical guidance noise field as physical guidance information, adaptively separating an effective seismic signal from the original seabed seismic data, and reconstructing and outputting a seabed seismic data set after noise suppression. Through combination of physical reasoning and deep learning, the model can distinguish signals like noise and noise of the image signals, and the high fidelity and generalization ability of seismic data noise reduction are improved.

The embodiment of the invention provides a construction system and a construction method for a submarine earthquake monitoring system. The system comprises a wind power platform, a workboat, an underwater robot and a seabed monitoring base station, wherein the wind power platform is used for arranging a winch, a lifting belt, a beam embracing clamp, a chain block and a pulley; the workboat is used for loading an underwater robot, a water entry bridge, a protection pipe, a seabed monitoring base station, a submarine cable, a winch, a sling, a beam embracing clamp, a chain block, a pulley, a detacher, an acoustic releaser and an operator. The underwater robot is used for monitoring drawing of the protection pipe; the seabed monitoring base station is used for monitoring seabed seismic data. By utilizing the existing structure and facilities of the existing wind power platform, a special monitoring platform or foundation does not need to be additionally built, the cost of foundation construction materials, manpower and time is effectively reduced, and a large amount of repeated cost input is avoided; and the monitoring of the pulling of the protection pipe can be completed by releasing the underwater robot, so that the labor cost is reduced, and the safety is improved.

It is proposed a seismic data acquisition device (400) intended to be placed on an ocean bottom floor, comprising a polymeric casing (412) defining a chamber that houses at least art of a data acquisition system (440, 444, 445); and a metallic device (414) in which the polymeric casing (412) is trapped, the metallic device (414) comprising two metallic beams (4141, 4142) that extend on opposite sides of the polymeric casing (412). It is also proposed a method for assembling such a device and a corresponding method for seabed seismic data acquisition.

The application relates to the technical field of seismic exploration, and provides a seabed seismic acquisition system, which comprises an acquisition station and an acquisition device. The acquisition device comprises a pressure cabin, a seismic detector arranged on the pressure cabin, the seismic detector being used for acquiring seismic reflection wave data, the pressure cabin having a first end and a second end opposite to each other, a first cone body, one end of the first cone body being connected with the first end of the pressure cabin, the other end of the first cone body extending away from the pressure cabin and being used for being inserted into the seabed, and a buffer reset assembly, two opposite ends of the buffer reset assembly being connected with the acquisition station and the pressure cabin respectively, the buffer reset assembly being used for making the pressure cabin and the acquisition station first approach each other and then move away from each other. The seabed seismic acquisition system can improve the coupling effect between the seabed seismic acquisition system and the seabed.

The present invention relates to a method for rapidly calculating the submarine seismic wave field induced by an axisymmetric volume sound source in water, comprising the steps of: constructing an equivalent sound source intensity - sound pressure relationship equation; constructing a series form of the Green's function; calculating the reflection coefficient under an elastic seabed; constructing a vibration velocity - equivalent sound source relationship equation; calculating the equivalent sound source intensity; constructing an axisymmetric finite element coupling calculation model for the equivalent source intensity; constructing an axisymmetric finite element coordinate system in cylindrical coordinates; establishing an axisymmetric finite element corresponding to the structure; establishing an axisymmetric finite element corresponding to the seawater; establishing an axisymmetric finite element corresponding to the acoustic - solid coupling boundary surface; establishing a matching infinite element for truncating the infinite outer domain; establishing an acoustic - vibration coupling dynamic equation under axisymmetric conditions; constructing an axisymmetric finite element coupling calculation model for the equivalent source intensity; and calculating the submarine seismic wave induced by an axisymmetric body target in water. The present invention realizes rapid and high - precision calculation of the submarine seismic wave field induced by an axisymmetric volume sound source in water.

A method of modeling a subsurface area of a seabed using fewer seismic sensors than a seismic sensors that generate sufficient data to permit a desired model of the subsurface area. Seismic datasets are received from the fewer seismic sensors. The fewer seismic sensors are located at points defined in relation to the seabed. Boundary conditions of a seismic wavefield at a fluid-solid interface between the seabed and water above the seabed are identified. An estimation, from the seismic datasets and the boundary conditions, is made of one or more virtual datasets. Each of the virtual datasets includes corresponding estimated data for additional points defined in relation to the seabed. The points and the additional points, combined, are sufficient to permit the desired model of the subsurface area. A subsurface model of the subsurface area is generated by modeling the subsurface area using the seismic datasets and the virtual datasets.

An embodiment of the present invention provides a bottom-seated submarine observation station for marine natural disaster research, which relates to the field of underwater operations. It includes: a protective cabin, a bracket, a mounting frame, a protective tube, a seabed seismic detector, and a bottom-seating mechanism; the bracket is mounted on the lower surface of the protective cabin; the mounting frame is mounted on the inner side of the bracket; the protective tube is mounted at the center of the mounting frame; the outer wall of the protective tube is provided with a plurality of through grooves equidistantly along the circumference; the seabed seismic detector is mounted in the protective tube, and the seabed seismic detector is connected to the lower surface of the protective cabin via a traction rope; two bottom-seating mechanisms are symmetrically mounted on the bracket; the bottom-seating mechanism can automatically unfold when it contacts the seabed, increasing the contact area between the device and the seabed. The device can naturally unfold the landing plate after landing, increasing the contact area with the seabed, and improve the positioning capability of the device by interlocking rivets with the seabed.

The application discloses a seabed slope model box for simulating seabed moving water pressure waves, which comprises a seabed slope model chamber, a moving water pressure wave device, a power equipment chamber and a weight lowering device. The weight lowering device can lift the moving water pressure wave device assembly to hit the water body at the bottom of the model box, so as to stimulate the moving water pressure wave generated by the fault dislocation caused by the seabed earthquake, and then the real marine water environment simulation around the seabed slope in the earthquake process is realized. The application has the advantages of simple structure, simple driving force provided by the weight stack, manual remote control of the moving water pressure wave stimulation time, and indoor seabed slope test conditions which are more in line with the actual working conditions and easy to observe the results.

The invention belongs to the technical field of seismic data acquisition, and particularly discloses an underwater seismic data continuous acquisition method, system, device and medium, the system comprises a main control module with a built-in watchdog mechanism, a sensor, an acquisition auxiliary module, a data storage module and a power supply module; the method comprises the steps that a main control module starts a data acquisition task, and a sensor acquires submarine seismic data and stores the submarine seismic data to a data storage module; in the data acquisition process, the main control module sends a dog feeding signal to the watchdog regularly, if the dog feeding signal is not sent on time, the watchdog forcibly restarts the main control module, and the main control module restarts the data acquisition task through the acquisition auxiliary module. According to the method, the problems of data acquisition interruption and time loss caused by abnormal programs of seabed nodes are solved, the continuity and integrity of marine seismic data acquisition are ensured, and the reliability and economic benefits of marine oil and gas exploration operation are improved. The method is suitable for continuous acquisition of underwater seismic data.

The invention provides an optimal earth radial parameter determination method and device, equipment and a medium, and the method comprises the steps: obtaining the original observation data of a plurality of seabed stations in a seabed earthquake epicentral region; based on the original observation data of each seabed station, determining the coordinate position resolving data of the seabed station and the observation value of the displacement time sequence; constructing a Burgers viscoelastic body earth model, and determining a reset error Love number and a complex green function based on the Burgers viscoelastic body earth model; determining a calculated value of the displacement time sequence according to the observed value of the displacement time sequence based on the wrong reset Love number and a complex green function; and determining an optimal transient viscosity parameter and an optimal steady-state viscosity parameter based on the observed value and the calculated value of the displacement time sequence, and jointly forming the optimal earth radial parameter by the transient viscosity parameter and the optimal steady-state viscosity parameter. By adopting the scheme, the radial parameters of the earth can be accurately determined, and correct understanding and accurate prediction of the earthquake process are ensured.

The application provides a seismic wave field analytical simulation method, comprising the following steps: constructing a seabed interface analytical equation according to a displacement potential function of a reflection / transmission wave and a boundary condition; obtaining a first function relationship; obtaining a second function relationship; transforming the seabed interface analytical equation by using the first function relationship and the second function relationship to obtain a function relationship formula of a reflection / transmission coefficient and an influence factor, and determining the function relationship formula as a third function relationship; obtaining a plurality of influence factors, and obtaining the reflection / transmission coefficients corresponding to the influence factors based on the third function relationship; and obtaining seabed seismic wave field data through forward modeling based on the one-to-one corresponding influence factors and reflection / transmission coefficients. The seabed seismic wave field obtained by the above seismic wave field analytical simulation method is highly similar to the real data, and can provide a theoretical model for the prediction and attenuation of seabed multiple waves and ghost waves.

The present application relates to the technical field of geophysical exploration, and particularly relates to a method and system for monitoring submarine earthquakes. Two interference optical fiber geological sensing rings with different physical parameters are deployed on a submarine to be monitored, and a phase demodulation interferometry method is used to obtain an observation signal vector; the observation signal vector is decomposed into a common seismic signal component and independent sparse interference components of each channel; an observation model is constructed based on a sparse dictionary transformation matrix, a wavelet domain sparse coefficient vector of each component, an observation signal vector, a sensing matrix and measurement noise; then, a joint optimization model is constructed according to the observation model, combined with a data error constraint term, a structured sparse constraint term and a physical prior constraint term, and solved to obtain a target sparse coefficient vector of each component; finally, a submarine earthquake monitoring result is obtained according to the target sparse coefficient vector of the seismic signal component, and precise and efficient submarine earthquake monitoring is achieved.