Horizontal component azimuth correction method, device and equipment of ocean bottom seismograph and medium

Abstract

The invention discloses a horizontal component azimuth correction method, device and equipment of a sea bottom seismograph and a medium. The method comprises the following steps: acquiring target direct water wave signals corresponding to a plurality of target shot point positions on a target survey line; for each target shot point position, target direct water wave total energy corresponding to the target direct water wave signal is determined; determining a plurality of candidate horizontal component orientations in a preset orientation range, and determining a target horizontal component orientation corresponding to the total energy of each target direct water wave from the plurality of candidate horizontal component orientations; and determining a target rotation angle according to a target orientation in the target horizontal component orientations corresponding to the plurality of target shot point positions, and correcting the initial horizontal component orientation according to the target rotation angle to obtain a reference horizontal component orientation. According to the scheme, the direct water wave maximum energy orientation is determined by using the direct water wave signal in the OBS horizontal component orientation so as to realize efficient and accurate correction of the horizontal component orientation, and the stability and practicability of horizontal component orientation correction are enhanced.

Description

Technical Field

[0001] The present invention relates to the technical field of marine seismic exploration, and particularly to a method, device, equipment and medium for horizontal component azimuth correction of an ocean bottom seismometer. Background Art

[0002] With the increasing attention of global scientists to the research on the deep structure of the ocean, the ocean bottom seismometer (OBS) has become a key tool for revealing marine geological tectonic activities. The OBS is a multi-component seismic observation system that deploys geophones on the seabed and realizes automatic recovery through acoustic signals, and has broad application prospects in marine geophysical exploration and research.

[0003] However, due to the inevitable rotation during the deployment process and the influence of the water flow in the sea, the directionality of the geophones in the OBS will inevitably change. The horizontal component direction of the OBS is of great significance for research such as natural earthquakes; in particular, for seismic receiver functions and converted shear waves, horizontal component rotation is fundamental. Therefore, in order to improve the basic instrument parameters of the OBS and subsequent seismological research, it is necessary to correct the horizontal component direction of the OBS.

[0004] In the related art, the maximum eigenvector value of the converted shear wave is solved by using the cross-energy matrix, so as to determine the polarization direction for rotation. Due to the large uncertainty in the detection of converted shear waves, such as the inability to detect converted shear waves or the poor quality of the detected converted shear waves, it is difficult to ensure the accuracy and efficiency of the horizontal component direction correction of the OBS, and the stability and practicability are limited. Summary of the Invention

[0005] The present invention provides a method, device, equipment and medium for horizontal component azimuth correction of an ocean bottom seismometer, which uses the direct water wave signals recorded in the horizontal component azimuth of the OBS to determine the maximum energy azimuth of the direct water wave, so as to achieve efficient and accurate correction of the horizontal component azimuth, and at the same time enhance the stability and practicability of the horizontal component azimuth correction.

[0006] According to one aspect of the present invention, there is provided a method for horizontal component azimuth correction of an ocean bottom seismometer, the method comprising:

[0007] Obtaining target direct water wave signals corresponding to a plurality of target shot point positions on a target survey line; wherein, the target direct water wave signal refers to the direct water wave phase received by the ocean bottom seismometer in the initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line;

[0008] For each of the target shot point positions, determine the total energy of the target direct water wave signal corresponding thereto; wherein, the total energy of the target direct water wave is the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth.

[0009] Determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine, from the plurality of candidate horizontal component azimuths, the target horizontal component azimuth corresponding to each of the total energies of the target direct water waves; wherein, the target horizontal component azimuth is the horizontal component azimuth with the maximum direct water wave energy in the target azimuth, and the target azimuth is the first azimuth or the second azimuth.

[0010] Determine a rotation angle as the target rotation angle according to the target azimuth among the target horizontal component azimuths corresponding to the plurality of target shot point positions, and correct the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth.

[0011] According to another aspect of the present invention, there is provided a horizontal component azimuth correction device for a marine seismograph, the device comprising:

[0012] A target direct water wave signal acquisition module, configured to acquire target direct water wave signals corresponding to a plurality of target shot point positions on a target survey line; wherein, the target direct water wave signal is a direct water wave phase received by the marine seismograph in the initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line.

[0013] A target direct water wave total energy determination module, configured to, for each of the target shot point positions, determine the total energy of the target direct water wave signal corresponding thereto; wherein, the total energy of the target direct water wave is the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth.

[0014] A target horizontal component azimuth determination module, configured to determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine, from the plurality of candidate horizontal component azimuths, the target horizontal component azimuth corresponding to each of the total energies of the target direct water waves; wherein, the target horizontal component azimuth is the horizontal component azimuth with the maximum direct water wave energy in the target azimuth, and the target azimuth is the first azimuth or the second azimuth.

[0015] A horizontal component azimuth correction module, configured to determine a rotation angle as the target rotation angle according to the target azimuth among the target horizontal component azimuths corresponding to the plurality of target shot point positions, and correct the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth.

[0016] According to another aspect of the present invention, there is provided an electronic device, which includes:

[0017] at least one processor; and,

[0018] a memory communicatively connected to the at least one processor; wherein,

[0019] the memory stores a computer program executable by the at least one processor, and when the computer program is executed by the at least one processor, the at least one processor is enabled to execute the horizontal component azimuth correction method of the seafloor seismograph according to any embodiment of the present invention.

[0020] According to another aspect of the present invention, there is provided a computer-readable storage medium storing computer instructions for causing a processor to implement the horizontal component azimuth correction method of the seafloor seismograph according to any embodiment of the present invention when executed.

[0021] The technical solution of the embodiment of the present invention first obtains target direct water wave signals corresponding to multiple target shot point positions on a target survey line; wherein, the target direct water wave signal refers to the direct water wave phase received by the seafloor seismograph at an initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line; then for each target shot point position, the total energy of the target direct water wave corresponding to the target direct water wave signal is determined respectively; wherein, the total energy of the target direct water wave refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth; furthermore, multiple candidate horizontal component azimuths are determined within a preset azimuth range, and the target horizontal component azimuth corresponding to each total energy of the target direct water wave is determined from the multiple candidate horizontal component azimuths; wherein, the target horizontal component azimuth refers to the horizontal component azimuth with the maximum direct water wave energy at the target azimuth, and the target azimuth is the first azimuth or the second azimuth; then a rotation angle is determined as the target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to the multiple target shot point positions, and the initial horizontal component azimuth is corrected according to the target rotation angle to obtain a reference horizontal component azimuth. This technical solution uses the direct water wave signals recorded on the OBS horizontal component azimuth to determine the azimuth with the maximum direct water wave energy, thereby achieving efficient and accurate correction of the horizontal component azimuth, and at the same time enhancing the stability and practicability of the horizontal component azimuth correction.

[0022] It should be understood that the content described in this part is not intended to identify the key or important features of the embodiments of the present invention, nor is it used to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] To more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the accompanying drawings required for the description of the embodiments. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings.

[0024] Figure 1 It is a flowchart of a method for horizontal component azimuth correction of a seafloor seismograph according to Embodiment 1 of the present invention;

[0025] Figure 2A It is a top view schematic diagram of the position of an OBS and the position of a seismic source according to Embodiment 1 of the present invention;

[0026] Figure 2B It is a schematic diagram of a target direct water wave signal according to Embodiment 1 of the present invention;

[0027] Figure 3A It is a single-shot energy search result graph according to Embodiment 1 of the present invention;

[0028] Figure 3B It is a schematic diagram of a target fitting curve according to Embodiment 1 of the present invention;

[0029] Figure 4 It is a flowchart of a method for horizontal component azimuth correction of a seafloor seismograph according to Embodiment 2 of the present invention;

[0030] Figure 5 It is a horizontal component profile before and after correction according to Embodiment 2 of the present invention;

[0031] Figure 6 It is a schematic structural diagram of a horizontal component azimuth correction device of a seafloor seismograph according to Embodiment 3 of the present invention;

[0032] Figure 7 It is a schematic structural diagram of an electronic device for implementing the horizontal component azimuth correction method of a seafloor seismograph in an embodiment of the present invention. Detailed implementation manners

[0033] In order to enable those skilled in the art to better understand the solutions of the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0034] It should be noted that the terms "first", "second", "target", etc. in the description, claims and above-mentioned drawings of the present invention are used to distinguish similar objects, and do not necessarily have to be used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present invention described here can be implemented in an order other than those illustrated or described here. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device that includes a series of steps or units does not have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products or devices.

[0035] Embodiment 1

[0036] Figure 1 FIG. is a flowchart of a method for horizontal component azimuth correction of a seafloor seismograph provided in Embodiment 1 of the present invention. This embodiment is applicable to the situation of efficiently and accurately correcting the horizontal component azimuth of an OBS. This method can be executed by a horizontal component azimuth correction device of a seafloor seismograph. The horizontal component azimuth correction device of the seafloor seismograph can be implemented in the form of hardware and / or software, and the horizontal component azimuth correction device of the seafloor seismograph can be configured in an electronic device with data processing capabilities. As Figure 1 shown, the method includes:

[0037] S110, obtaining target direct water wave signals corresponding to multiple target shot point positions on a target survey line.

[0038] Among them, the target survey line can refer to a set of a series of shot points set in a specific direction in the exploration area for collecting seismic wave data, that is, multiple shot points are set on the target survey line, and the target survey line is a straight line. The target shot point position can refer to the position where seismic waves are artificially excited and can be used to represent the source position.

[0039] Among them, the target direct water wave signal refers to the direct water wave phase received by the ocean bottom seismograph (OBS) in the initial horizontal component azimuth, and the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other. It should be noted that the OBS can be used to record multi-component seismic wave signals, where the multi-components include a vertical component and a horizontal component. Specifically, the vertical component refers to the component perpendicular to the seabed surface and can be used to record the longitudinal wave field; the horizontal component refers to the component parallel to the seabed surface and includes two mutually perpendicular components, which can be used to record the propagation characteristics of shear waves and longitudinal waves in the horizontal direction. The horizontal component azimuth can be used to characterize the azimuth where the horizontal component is located and includes a first azimuth and a second azimuth that are orthogonal to each other. Exemplarily, the first azimuth and the second azimuth can be the east direction and the north direction in the geographical azimuth respectively. The initial horizontal component azimuth can refer to the OBS horizontal component azimuth before calibration. The direct water wave phase can refer to the attenuated seismic wave signal of the original seismic wave signal emitted from the target shot point position along the shortest propagation path reaching the OBS detector. The target direct water wave signal is the direct water wave phase received by the OBS detector in the initial horizontal component azimuth, specifically including the direct water wave phase in the first azimuth and the direct water wave phase in the second azimuth.

[0040] Figure 2A This is a top view schematic diagram of the OBS position and the source position provided by Embodiment 1 of the present invention. As Figure 2A shown, the black dashed line represents the target survey line, the red mark represents the OBS position, the yellow mark represents the target shot point position (i.e., the source position), and the blue line represents the shortest path of signal transmission. Specifically, it can be set that the exploration ship sails along the target survey line and fires shots with an air gun at multiple positions in an equal-time or equal-spacing manner during the sailing process, thereby emitting corresponding seismic wave signals. When the air gun fires a shot, the on-board GPS will record the accurate position information of the air gun, and the air gun firing position is the corresponding shot point position, whereby multiple target shot point positions on the target survey line can be determined. Figure 2B This is a schematic diagram of a target direct water wave signal provided by Embodiment 1 of the present invention. As Figure 2B shown, SHE and SHN respectively represent the first azimuth and the second azimuth in the initial horizontal component azimuth.

[0041] In this embodiment, multiple shot points within the preset range of the OBS can be selected as target shot points, or multiple shot points with a preset number near the OBS can be selected as target shot points. Exemplarily, the preset range can be a circular area within 5 km from the OBS. Among them, the positions and firing times of the target shot points can be accurately obtained. The OBS can calculate the signal reception time when the detector receives the direct water wave phase according to the firing time of the target shot point and the seawater depth, so as to locate which target shot point position the received signal comes from.

[0042] S120. For each target shot point position, respectively determine the total energy of the target direct water wave signal corresponding thereto.

[0043] Among them, the total energy of the target direct water wave refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth. It can be understood that since the horizontal component azimuth includes the first azimuth and the second azimuth that are orthogonal to each other, correspondingly, the target direct water wave signal includes the direct water wave phase in the first azimuth and the direct water wave phase in the second azimuth. Therefore, the total energy of the target direct water wave signal in the entire horizontal direction can be determined based on the target direct water wave signals corresponding to the first azimuth and the second azimuth as the total energy of the target direct water wave.

[0044] In this embodiment, optionally, for each target shot point position, respectively determining the total energy of the target direct water wave signal corresponding thereto includes: for each target shot point position, respectively performing time integration on the amplitude of the target direct water wave signal in the first azimuth of the initial horizontal component azimuth to obtain the first direct water wave energy; for each target shot point position, respectively performing time integration on the amplitude of the target direct water wave signal in the second azimuth of the initial horizontal component azimuth to obtain the second direct water wave energy; determining the total energy of the target direct water wave signal corresponding thereto according to the sum of the squares of the first direct water wave energy and the second direct water wave energy.

[0045] Specifically, based on the relationship between the signal amplitude and the signal energy, the corresponding signal energy can be obtained by performing time integration on the signal amplitude. In this embodiment, for each target shot point position, respectively performing time integration on the amplitudes of the target direct water wave signals in the first azimuth and the second azimuth of the initial horizontal component azimuth to obtain the first direct water wave energy and the second direct water wave energy, and then calculating the sum of the squares of the first direct water wave energy and the second direct water wave energy as the total energy of the target direct water wave signal corresponding thereto.

[0046] Through such a setting, this solution can quickly and accurately determine the total energy of the target direct water wave signal corresponding to each target shot point position by using the relationship between the signal amplitude and the signal energy.

[0047] S130. Determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine the target horizontal component azimuth corresponding to each total energy of the target direct water wave from the plurality of candidate horizontal component azimuths.

[0048] Among them, the preset azimuth range can refer to the azimuth search range preset according to actual needs, and can be used to represent the angular deviation range between the candidate horizontal component azimuth and the initial horizontal component azimuth. Exemplarily, the preset azimuth range can be set to 0-180 degrees or 0-360 degrees. Among them, the candidate horizontal component azimuth can refer to the horizontal component azimuth obtained by rotating a certain angle with the initial horizontal component azimuth as the reference. Among them, the target horizontal component azimuth refers to the horizontal component azimuth with the maximum direct water wave energy at the target azimuth, and the target azimuth is the first azimuth or the second azimuth. It can be understood that the target horizontal component azimuth can be the horizontal component azimuth with the maximum direct water wave energy at the first azimuth, or the horizontal component azimuth with the maximum direct water wave energy at the second azimuth. It should be noted that the sum of the direct water wave energies at the first azimuth and the second azimuth among different candidate horizontal component azimuths is constant, that is, the total target direct water wave energy. Therefore, if the total target direct water wave energy is decomposed into different candidate horizontal component azimuths, according to the Pythagorean theorem of a triangle, when the direct water wave energy of one azimuth of the candidate horizontal component azimuth increases, the direct water wave energy of the other orthogonal azimuth will inevitably decrease.

[0049] In this embodiment, optionally, determining a plurality of candidate horizontal component azimuths within the preset azimuth range includes: determining a plurality of candidate rotation angles within the preset azimuth range based on a preset azimuth interval; and rotating the initial horizontal component azimuth according to each candidate rotation angle to obtain a plurality of candidate horizontal component azimuths.

[0050] Among them, the preset azimuth interval can refer to the azimuth search interval preset according to actual needs. Exemplarily, the preset azimuth interval can be set to 0.5 degrees, 1 degree or 2 degrees, specifically depending on the accuracy required for azimuth search. The candidate rotation angle can be used to represent the rotation angle of the initial horizontal component azimuth determined based on the preset azimuth interval within the preset azimuth range. The candidate horizontal component azimuth can refer to the horizontal component azimuth obtained by rotating the initial horizontal component azimuth by the candidate rotation angle.

[0051] Exemplarily, assuming that the preset azimuth range is 0-360 degrees and the preset azimuth interval is 1 degree, then the candidate rotation angles can be determined to be 1-360 degrees, and the candidate horizontal component azimuths are the horizontal component azimuths obtained by rotating the initial horizontal component azimuth by 1-360 degrees respectively.

[0052] In this embodiment, after determining a plurality of candidate horizontal component orientations, the target horizontal component orientation corresponding to the total energy of each target direct water wave can be determined from the plurality of candidate horizontal component orientations. Optionally, determining the target horizontal component orientation corresponding to the total energy of each target direct water wave from the plurality of candidate horizontal component orientations includes: for each target shot point position, decomposing the total energy of the target direct water wave onto the first orientation and the second orientation of each candidate horizontal component orientation; and taking the candidate horizontal component orientation with the maximum direct water wave energy in the target orientation as the target horizontal component orientation corresponding to the total energy of the target direct water wave.

[0053] Specifically, after determining a plurality of candidate horizontal component orientations, the total energy of the target direct water wave corresponding to each target shot point position can be respectively decomposed onto the first orientation and the second orientation of each candidate horizontal component orientation, and the candidate horizontal component orientation with the maximum direct water wave energy in the target orientation (i.e., the first orientation or the second orientation) can be selected as the target horizontal component orientation corresponding to the total energy of the target direct water wave. It can be understood that for the plurality of candidate horizontal component orientations, the direct water wave energy in the first orientation or the second orientation of the target horizontal component orientation is the maximum.

[0054] Figure 3A FIG. 1 is a single-shot energy search result diagram provided in Embodiment 1 of the present invention. As Figure 3A shown, SHR and SHT respectively represent the first orientation and the second orientation in the candidate horizontal component orientations, that is, SHR and SHT respectively correspond to SHE and SHN, the abscissa and the ordinate respectively represent the candidate rotation angle and the direct water wave energy, and the preset orientation range is 0-360 degrees. Figure 3A FIG. 2 shows the energy search result corresponding to the target shot point position with shot number 2420. Among them, the shot number can be used to uniquely identify the shot point position. According to Figure 3A it can be known that taking SHR as an example, the direct water wave energy corresponding to 20 degrees is the maximum. Therefore, the candidate horizontal component orientation corresponding to 20 degrees can be taken as the target horizontal component orientation of the target shot point position 2420.

[0055] Through such a setting in this solution, by decomposing the total energy of the target direct water wave corresponding to each target shot point position onto each candidate horizontal component orientation, the target horizontal component orientation corresponding to the total energy of the target direct water wave under each target shot point position can be quickly and accurately located according to the decomposition result.

[0056] S140. Determine a rotation angle as the target rotation angle according to the target orientation in the target horizontal component orientations corresponding to the plurality of target shot point positions, and correct the initial horizontal component orientation according to the target rotation angle to obtain the reference horizontal component orientation.

[0057] Among them, the target rotation angle may refer to the angle that the initial horizontal component azimuth actually needs to be rotated. In this embodiment, after determining multiple target horizontal component azimuths, a rotation angle can be determined as the target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to multiple target shot positions. Optionally, determining a rotation angle as the target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to multiple target shot positions includes: determining a target fitting curve according to the target azimuth in the target horizontal component azimuth corresponding to each target shot position; determining a target curve boundary value according to the change trend of the target fitting curve, and determining the target curve boundary value as the target rotation angle.

[0058] Specifically, first, a smooth curve is used to connect the target azimuths in the target horizontal component azimuths corresponding to each target shot position to obtain a target fitting curve, and then the boundary value of the target fitting curve is determined according to the change trend of the target fitting curve as the target curve boundary value, and the target curve boundary value is determined as the target rotation angle. Figure 3B It is a schematic diagram of a target fitting curve provided by Embodiment 1 of the present invention. As Figure 3B shown, the abscissa and ordinate respectively represent the shot number and the rotation angle corresponding to the target azimuth in the target horizontal component azimuth, and the blue curve represents the target fitting curve. From Figure 3B it can be seen that the target curve boundary value is 75 degrees marked by the red line, and at this time, the target rotation angle can be determined as 75 degrees.

[0059] Correspondingly, correcting the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth includes: rotating the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth, so that the target azimuth corresponding to the reference horizontal component azimuth is parallel to the direction where the target survey line is located.

[0060] In this embodiment, after determining the target curve boundary value as the target rotation angle, the initial horizontal component azimuth can be rotated by the target rotation angle to obtain a reference horizontal component azimuth, so that the target azimuth corresponding to the reference horizontal component azimuth is parallel to the direction where the target survey line is located. At this time, the correction task of the OBS horizontal component azimuth can be completed.

[0061] It should be noted that when the seismic wave signal emitted from the target shot position travels along the shortest path (see Figure 2AWhen the blue line in [[]] is transmitted, the direct water wave energy at the target azimuth in the corresponding horizontal component azimuth is the largest. At this time, the horizontal component azimuth is the target horizontal component azimuth corresponding to the target shot point position. For the shot points on the target survey line that are infinitely far away from the OBS (i.e., the two ends of the target survey line are infinitely far away), the angle between the target survey line and the shortest signal transmission path is almost 0 degrees. At this time, it can be considered that the direction of the target survey line is parallel to the direction of the shortest signal transmission path. Based on the change trend of the target fitting curve, it can be known that the target curve boundary value (i.e., the target rotation angle) is the angle that the target horizontal component azimuth corresponding to the shot point positions at both ends of the target survey line that are infinitely far away needs to rotate. That is to say, for the shot point positions at both ends of the target survey line that are infinitely far away, after rotating the initial horizontal component azimuth by the target rotation angle, the reference horizontal component azimuth is obtained, and the target azimuth corresponding to the reference horizontal component azimuth is the direction of the shortest signal transmission path of the signal of this shot point. Since the direction of the shortest signal transmission path of the signal of this shot point is parallel to the direction of the target survey line, at this time, the target azimuth corresponding to the reference horizontal component azimuth is parallel to the direction of the target survey line.

[0062] Through such a setting in this solution, the target rotation angle corresponding to the shot points at both ends of the target survey line that are infinitely far away can be determined by using the change trend of the target fitting curve. Rotating the initial horizontal component azimuth according to the target rotation angle makes the target azimuth corresponding to the rotated reference horizontal component azimuth parallel to the direction of the target survey line, so as to achieve the rapid and accurate correction of the horizontal component azimuth of the OBS.

[0063] The technical solution of the embodiment of the present invention is as follows: First, obtain the target direct water wave signals corresponding to multiple target shot point positions on the target survey line; wherein, the target direct water wave signal refers to the direct water wave phase received by the ocean bottom seismograph in the initial horizontal component azimuth, and the horizontal component azimuth includes the first azimuth and the second azimuth that are orthogonal to each other, and the target survey line is a straight line; then, for each target shot point position, determine the total energy of the target direct water wave corresponding to the target direct water wave signal respectively; wherein, the total energy of the target direct water wave refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth; furthermore, determine multiple candidate horizontal component azimuths within the preset azimuth range, and determine the target horizontal component azimuth corresponding to each total energy of the target direct water wave from the multiple candidate horizontal component azimuths; wherein, the target horizontal component azimuth refers to the horizontal component azimuth with the maximum direct water wave energy in the target azimuth, and the target azimuth is the first azimuth or the second azimuth; then, determine a rotation angle as the target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to the multiple target shot point positions, and correct the initial horizontal component azimuth according to the target rotation angle to obtain the reference horizontal component azimuth. This technical solution uses the direct water wave signals recorded in the horizontal component azimuth of the OBS to determine the azimuth with the maximum direct water wave energy, thereby realizing the efficient and accurate correction of the horizontal component azimuth, and at the same time enhancing the stability and practicality of the horizontal component azimuth correction.

[0064] Embodiment 2

[0065] Figure 4 The flowchart of a method for correcting the horizontal component azimuth of an ocean bottom seismograph provided by Embodiment 2 of the present invention is shown. This embodiment is optimized based on the above embodiment. The specific optimization is as follows: After obtaining the target direct water wave signals corresponding to multiple target shot point positions on the target survey line, the method further includes: determining the signal-to-noise ratio corresponding to the target direct water wave signal; if the signal-to-noise ratio is lower than the preset threshold, the corresponding target direct water wave signal is removed.

[0066] As Figure 4 shown, the method of this embodiment specifically includes the following steps:

[0067] S210, obtain the target direct water wave signals corresponding to multiple target shot point positions on the target survey line, and determine the signal-to-noise ratio corresponding to the target direct water wave signal.

[0068] Wherein, the target direct water wave signal refers to the direct water wave phase received by the ocean bottom seismograph in the initial horizontal component azimuth, and the horizontal component azimuth includes the first azimuth and the second azimuth that are orthogonal to each other, and the target survey line is a straight line.

[0069] It should be noted that if the seismic wave signal emitted by the air gun is interfered by other signals (such as noise) during the transmission process, it may cause the target direct water wave signal received by the OBS detector to be chaotic, thus affecting the accuracy of the azimuth correction of the OBS horizontal component. Therefore, after obtaining the target direct water wave signals corresponding to multiple target gunpoint positions on the target survey line, the signal-to-noise ratio corresponding to the target direct water wave signal can also be calculated, and the degree of interference of the target direct water wave signal can be reflected through the signal-to-noise ratio. It can be understood that the lower the signal-to-noise ratio, the more severely the signal is interfered.

[0070] S220, if the signal-to-noise ratio is lower than the preset threshold, then the corresponding target direct water wave signal is excluded.

[0071] Among them, the preset threshold can refer to the signal-to-noise ratio reference value set according to actual requirements. For example, the preset threshold can be set to 10. In this embodiment, if it is determined that the signal-to-noise ratio corresponding to the target direct water wave signal is lower than the preset threshold, it indicates that the target direct water wave signal is severely interfered. At this time, in order to ensure the accuracy and availability of the target direct water wave signal, the target direct water wave signal with a signal-to-noise ratio lower than the preset threshold can be excluded.

[0072] S230, for each target gunpoint position, respectively determine the total energy of the target direct water wave corresponding to the target direct water wave signal.

[0073] Among them, the total energy of the target direct water wave refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth.

[0074] S240, determine multiple candidate horizontal component azimuths within the preset azimuth range, and determine the target horizontal component azimuth corresponding to each total energy of the target direct water wave from the multiple candidate horizontal component azimuths.

[0075] Among them, the target horizontal component azimuth refers to the horizontal component azimuth with the maximum direct water wave energy in the target azimuth, and the target azimuth is the first azimuth or the second azimuth.

[0076] S250, determine a rotation angle as the target rotation angle according to the target azimuth among the target horizontal component azimuths corresponding to multiple target gunpoint positions, and correct the initial horizontal component azimuth according to the target rotation angle to obtain the reference horizontal component azimuth.

[0077] Among them, the specific implementation manners of S230 - S250 can refer to the relevant descriptions in Embodiment 1, and will not be elaborated here.

[0078] In the technical solution of the embodiment of the present invention, after obtaining the target direct water wave signals corresponding to multiple target shot point positions on the target survey line, the signal-to-noise ratio corresponding to the target direct water wave signals is determined; if the signal-to-noise ratio is lower than a preset threshold, the corresponding target direct water wave signals are excluded, thereby ensuring the accuracy and availability of the target direct water wave signals, which helps to improve the accuracy of horizontal component azimuth correction.

[0079] In this embodiment, optionally, after correcting the initial horizontal component azimuth according to the target rotation angle to obtain the reference horizontal component azimuth, the method further includes: determining the target characteristic parameters of the initial horizontal component azimuth and the reference horizontal component azimuth; performing quality control verification on the reference horizontal component azimuth according to the change of the target characteristic parameters.

[0080] Among them, the target characteristic parameters can be used as an important basis for quality control verification of the horizontal component azimuth correction result. Exemplarily, the target characteristic parameters can be direct water wave energy or converted shear wave intensity, etc.

[0081] In this embodiment, in order to verify the effectiveness of the horizontal component azimuth correction, a quality control verification link is specially set up, and the effectiveness of the horizontal component azimuth correction result is tested by using the target characteristic parameters of the horizontal component azimuth before and after correction. Specifically, after correcting the initial horizontal component azimuth according to the target rotation angle to obtain the reference horizontal component azimuth, the target characteristic parameters of the initial horizontal component azimuth and the reference horizontal component azimuth are respectively determined, and then the quality control verification is performed on the reference horizontal component azimuth according to the change of the target characteristic parameters, so as to determine whether the horizontal component azimuth correction is effective. For example, other shot points on the target survey line except the target shot point can be used to perform quality control verification on the reference horizontal component azimuth by comparing the change of the direct water wave energy on the horizontal component azimuth before and after correction. If the total energy of the target direct water wave corresponding to this shot point is concentrated on a certain azimuth in the reference horizontal component azimuth compared with the initial horizontal component azimuth, it indicates that the horizontal component azimuth correction effect is good. In addition, any shot point on the target survey line can be used to perform quality control verification on the reference horizontal component azimuth by comparing the change of the converted shear wave intensity corresponding to the horizontal component azimuth before and after correction. If the converted shear wave intensity corresponding to the reference horizontal component azimuth increases compared with the initial horizontal component azimuth, it indicates that the horizontal component azimuth correction is effective.

[0082] Figure 5 This is a horizontal component profile before and after correction provided for the second embodiment of the present invention. Among them, the left side is the horizontal component profile before correction, and the right side is the horizontal component profile after correction. The abscissa and ordinate respectively represent the offset and time. As Figure 5 shown, the converted shear wave phase can be clearly seen from the horizontal component (OBS-SHR) profile after correction, which can prove the effectiveness of the horizontal component azimuth correction.

[0083] With such a setting, this solution can utilize the changes in the target feature parameters before and after the horizontal component azimuth correction to perform quality control verification on the reference horizontal component azimuth, so as to verify the effectiveness of the horizontal component azimuth correction.

[0084] Embodiment III

[0085] Figure 6 FIG. is a schematic structural diagram of a horizontal component azimuth correction device for a marine seismograph provided in Embodiment III of the present invention. This device can execute the horizontal component azimuth correction method for a marine seismograph provided in any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method. As Figure 6 shown, the device includes:

[0086] A target direct water wave signal acquisition module 310, configured to acquire target direct water wave signals corresponding to a plurality of target shot point positions on a target survey line; wherein, the target direct water wave signal refers to the direct water wave phase received by the marine seismograph at the initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line;

[0087] A target direct water wave total energy determination module 320, configured to respectively determine the target direct water wave total energy corresponding to the target direct water wave signal for each of the target shot point positions; wherein, the target direct water wave total energy refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth;

[0088] A target horizontal component azimuth determination module 330, configured to determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine the target horizontal component azimuth corresponding to each of the target direct water wave total energies from the plurality of candidate horizontal component azimuths; wherein, the target horizontal component azimuth refers to the horizontal component azimuth with the maximum direct water wave energy at the target azimuth, and the target azimuth is the first azimuth or the second azimuth;

[0089] A horizontal component azimuth correction module 340, configured to determine a rotation angle as the target rotation angle according to the target azimuth among the target horizontal component azimuths corresponding to the plurality of target shot point positions, and correct the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth.

[0090] Optionally, the target direct water wave total energy determination module 320 is specifically configured to:

[0091] For each of the target shot point positions, respectively perform time integration on the amplitudes of the target direct water wave signals on the first azimuth in the initial horizontal component azimuth to obtain the first direct water wave energy;

[0092] For each of the target shot point positions, respectively, the amplitude of the target direct water wave signal at the second azimuth in the initial horizontal component azimuth is time-integrated to obtain the second direct water wave energy;

[0093] The target direct water wave total energy corresponding to the target direct water wave signal is determined according to the square sum of the first direct water wave energy and the second direct water wave energy.

[0094] Optionally, the target horizontal component azimuth determination module 330 is configured to:

[0095] Determining a plurality of candidate rotation angles within the preset orientation range based on a preset orientation interval;

[0096] The initial horizontal component orientation is rotated according to each candidate rotation angle to obtain a plurality of candidate horizontal component orientations.

[0097] Optionally, the target horizontal component azimuth determination module 330 is further configured to:

[0098] For each of the target shot point positions, decomposing the target direct water wave total energy into the first and second azimuths of each of the candidate horizontal component azimuths;

[0099] The candidate horizontal component azimuth with the maximum direct water wave energy on the target azimuth is used as the target horizontal component azimuth corresponding to the total energy of the direct water wave of the target.

[0100] Optionally, the horizontal component azimuth correction module 340 is used to:

[0101] Determine a target fitting curve according to the target azimuth in the target horizontal component azimuth corresponding to each of the target shot point positions;

[0102] Determine a target curve boundary value according to a change trend of the target fitting curve, and determine the target curve boundary value as a target rotation angle;

[0103] The initial horizontal component azimuth is rotated according to the target rotation angle to obtain a reference horizontal component azimuth, so that the target azimuth corresponding to the reference horizontal component azimuth is parallel to the direction of the target survey line.

[0104] Optionally, the device further includes: a data processing module, configured to:

[0105] After obtaining target direct water wave signals corresponding to a plurality of target shot point positions on the target survey line, determining a signal-to-noise ratio corresponding to the target direct water wave signals;

[0106] If the signal-to-noise ratio is lower than a preset threshold, the corresponding target direct water wave signal is eliminated.

[0107] Optionally, the device further includes: a quality control verification module, configured to:

[0108] After correcting the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth, determine target characteristic parameters of the initial horizontal component azimuth and the reference horizontal component azimuth;

[0109] Perform quality control verification on the reference horizontal component azimuth according to the change of the target characteristic parameters.

[0110] The horizontal component azimuth correction device of a marine seismograph provided by an embodiment of the present invention can execute a horizontal component azimuth correction method of a marine seismograph provided by any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method.

[0111] Embodiment 4

[0112] Figure 7 FIG. shows a schematic structural diagram of an electronic device 10 that can be used to implement an embodiment of the present invention. The electronic device is intended to represent various forms of digital computers, such as, a laptop computer, a desktop computer, a workbench, a personal digital assistant, a server, a blade server, a mainframe computer, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as, a personal digital processor, a cellular phone, a smart phone, a wearable device (such as a helmet, glasses, a watch, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present invention described and / or claimed herein.

[0113] As Figure 7 shown, the electronic device 10 includes at least one processor 11, and a memory communicatively connected to at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc. Among them, the memory stores a computer program executable by at least one processor. The processor 11 can execute various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. The input / output (I / O) interface 15 is also connected to the bus 14.

[0114] Multiple components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disc, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0115] The processor 11 can be various general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the horizontal component azimuth correction method of the subsea seismograph.

[0116] In some embodiments, the horizontal component azimuth correction method of the subsea seismograph can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as the storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed onto the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the horizontal component azimuth correction method of the subsea seismograph described above can be executed. Alternatively, in other embodiments, the processor 11 can be configured to execute the horizontal component azimuth correction method of the subsea seismograph in any other suitable way (e.g., by means of firmware).

[0117] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems on a chip (SOCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include: implemented in one or more computer programs, the one or more computer programs can be executed and / or interpreted on a programmable system including at least one programmable processor, the programmable processor can be a special or general-purpose programmable processor, can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit the data and instructions to the storage system, the at least one input device, and the at least one output device.

[0118] A computer program for implementing the method of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general purpose computer, a special purpose computer, or other programmable data processing apparatus, such that the computer programs, when executed by the processor, cause the functions / operations specified in the flowchart and / or block diagram to be implemented. The computer program may be executed entirely on the machine, partly on the machine, as a stand-alone software package partly on the machine and partly on a remote machine, or entirely on the remote machine or server.

[0119] In the context of the present invention, a computer-readable storage medium may be a tangible medium that can contain, or store a computer program for use in or in connection with an instruction execution system, apparatus, or device. The computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. Alternatively, the computer-readable storage medium may be a machine-readable signal medium. More specific examples of the machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0120] In order to provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the electronic device. Other kinds of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0121] The systems and techniques described herein can be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system can be interconnected with each other by digital data communication in any form or medium (e.g., a communication network). Examples of communication networks include: local area network (LAN), wide area network (WAN), blockchain network, and the Internet.

[0122] A computing system can include a client and a server. The client and the server are generally remote from each other and typically interact through a communication network. The client-server relationship is created by computer programs that run on the respective computers and have a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or a cloud host, which is a host product in the cloud computing service system, and solves the defects of difficult management and weak business scalability existing in traditional physical hosts and VPS services.

[0123] It should be understood that various forms of the processes shown above can be used, with steps reordered, added, or deleted. For example, the steps recited in the present invention can be executed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of the present invention can be achieved, and this is not limited herein.

[0124] The above specific embodiments do not constitute a limitation on the protection scope of the present invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims

1. A method for correcting the horizontal component azimuth of a seafloor seismograph, characterized in that: The method comprises: Obtaining target direct water wave signals corresponding to multiple target shot point positions on the target survey line; wherein the target direct water wave signal refers to the direct water wave phase received by the seafloor seismograph in the initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line; For each of the target shot point positions, the target direct water wave total energy corresponding to the target direct water wave signal is determined respectively; wherein the target direct water wave total energy refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth; Determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine a target horizontal component azimuth corresponding to each target direct water wave total energy from the plurality of candidate horizontal component azimuths; wherein the target horizontal component azimuth refers to the horizontal component azimuth with the largest direct water wave energy on the target azimuth, and the target azimuth is the first azimuth or the second azimuth; A rotation angle is determined as a target rotation angle according to a target azimuth in the target horizontal component azimuths corresponding to the multiple target shot point positions, and the initial horizontal component azimuth is corrected according to the target rotation angle to obtain a reference horizontal component azimuth.

2. The method according to claim 1, characterized in that For each of the target shot point positions, the target direct water wave total energy corresponding to the target direct water wave signal is determined respectively, including: For each of the target shot point positions, respectively, the amplitude of the target direct water wave signal at the first azimuth in the initial horizontal component azimuth is time-integrated to obtain the first direct water wave energy; For each of the target shot point positions, respectively, the amplitude of the target direct water wave signal at the second azimuth in the initial horizontal component azimuth is time-integrated to obtain the second direct water wave energy; The target direct water wave total energy corresponding to the target direct water wave signal is determined according to the square sum of the first direct water wave energy and the second direct water wave energy.

3. The method according to claim 1, characterized in that A plurality of candidate horizontal component orientations are determined within a preset orientation range, including: Determining a plurality of candidate rotation angles within the preset orientation range based on a preset orientation interval; The initial horizontal component orientation is rotated according to each candidate rotation angle to obtain a plurality of candidate horizontal component orientations.

4. The method according to claim 3, characterized in that Determining the target horizontal component azimuth corresponding to the total energy of the direct water wave of each target from the multiple candidate horizontal component azimuths includes: For each of the target shot point positions, decomposing the target direct water wave total energy into the first and second azimuths of each of the candidate horizontal component azimuths; The candidate horizontal component azimuth with the maximum direct water wave energy on the target azimuth is used as the target horizontal component azimuth corresponding to the total energy of the direct water wave of the target.

5. The method according to claim 1, characterized in that Determining a rotation angle as a target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to the multiple target shot point positions includes: Determine a target fitting curve according to the target azimuth in the target horizontal component azimuth corresponding to each of the target shot point positions; Determine a target curve boundary value according to a change trend of the target fitting curve, and determine the target curve boundary value as a target rotation angle; Accordingly, the initial horizontal component azimuth is corrected according to the target rotation angle to obtain a reference horizontal component azimuth, including: The initial horizontal component azimuth is rotated according to the target rotation angle to obtain a reference horizontal component azimuth, so that the target azimuth corresponding to the reference horizontal component azimuth is parallel to the direction of the target survey line.

6. The method according to any one of claims 1 to 5, characterized in that: After obtaining target direct water wave signals corresponding to multiple target shot point positions on the target survey line, the method further includes: Determine the signal-to-noise ratio corresponding to the target direct water wave signal; If the signal-to-noise ratio is lower than a preset threshold, the corresponding target direct water wave signal is eliminated.

7. The method according to claim 6, characterized in that After correcting the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth, the method further includes: Determining target characteristic parameters of the initial horizontal component orientation and the reference horizontal component orientation; The reference horizontal component orientation is quality-controlled and verified according to the change of the target characteristic parameter.

8. A horizontal component azimuth correction device for a seafloor seismograph, characterized in that: The device comprises: A target direct water wave signal acquisition module is used to acquire target direct water wave signals corresponding to multiple target shot point positions on a target survey line; wherein the target direct water wave signal refers to the direct water wave phase received by the seafloor seismograph in the initial horizontal component azimuth, the horizontal component azimuth includes a first azimuth and a second azimuth that are orthogonal to each other, and the target survey line is a straight line; A target direct water wave total energy determination module is used to determine the target direct water wave total energy corresponding to the target direct water wave signal for each target shot point position; wherein the target direct water wave total energy refers to the total energy of the direct water wave phase corresponding to the initial horizontal component azimuth; A target horizontal component azimuth determination module is used to determine a plurality of candidate horizontal component azimuths within a preset azimuth range, and determine the target horizontal component azimuth corresponding to each target direct water wave total energy from the plurality of candidate horizontal component azimuths; wherein the target horizontal component azimuth refers to the horizontal component azimuth with the largest direct water wave energy on the target azimuth, and the target azimuth is the first azimuth or the second azimuth; The horizontal component azimuth correction module is used to determine a rotation angle as a target rotation angle according to the target azimuth in the target horizontal component azimuths corresponding to the multiple target shot point positions, and to correct the initial horizontal component azimuth according to the target rotation angle to obtain a reference horizontal component azimuth.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and, a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can perform the horizontal component azimuth correction method for the seafloor seismograph according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the horizontal component azimuth correction method of an ocean bottom seismograph according to any one of claims 1 to 7 when executed.

Images