CN117706618A

CN117706618A - R-T component calculation method and device for OBS seismic data

Info

Publication number: CN117706618A
Application number: CN202211115712.9A
Authority: CN
Inventors: 乌云巴特尔; 杨丽莹; 张济东
Original assignee: Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd; China Petroleum and Chemical Corp
Current assignee: Sinopec Petroleum Geophysical Exploration Technology Research Institute Co ltd; China Petroleum and Chemical Corp
Priority date: 2022-09-14
Filing date: 2022-09-14
Publication date: 2024-03-15

Abstract

The disclosure relates to the technical field of R-T component calculation of OBS seismic data, and provides an R-T component calculation method and device of OBS seismic data. The method comprises the following steps: acquiring a plurality of seismic trace data, each seismic trace data comprising a first arrival signal; projecting the amplitude polarity of the first arrival signal to a target coordinate system, and calculating the horizontal component azimuth angle of the target OBS; calculating a polarization angle based on the projected coordinates of the first arrival signal and the horizontal component azimuth of the OBS; and rotating an X-Y component in the horizontal component of the seismic channel data of the corresponding seismic channels to an R-T direction based on the polarization angle to obtain R-T components of the seismic channel data of the plurality of seismic channels. The embodiments of the present disclosure may accurately obtain a polarization angle and accurately rotate data in a seismic trace from a horizontal component to an R-T component based on the polarization angle.

Description

R-T component calculation method and device for OBS seismic data

Technical Field

The disclosure relates to the technical field of marine seismic exploration, in particular to an R-T component calculation method and device of OBS seismic data.

Background

Multi-component exploration technology begins in the 70 s of the 20 th century in oil and gas exploration, and exploration modes include transverse wave exploration, terrestrial three-dimensional three-component exploration, marine four-component exploration and the like. At present, offshore multicomponent exploration mainly adopts OBS seismic exploration, which is longitudinal wave excitation and simultaneously receives longitudinal wave and converted wave information. In the current data processing links of offshore OBS multi-component exploration, horizontal component rotation is an indispensable link.

In OBS (Ocean bottom seismometer, ocean bottom seismograph) seismic exploration, seismic wave energy excited by a shot on the sea surface mainly propagates on a vertical plane of shot and wave detection point connecting lines, the connecting line direction of the shot and the wave detection points is called Radial (R), the direction perpendicular to the Radial direction is called tangential (T), and converted waves are mainly recorded on the Radial direction (R) and the tangential direction (T). In OBS seismic exploration, a detector consists of an X-component receiver that receives primarily shear wave information, a Y-component receiver that receives primarily compressional wave information, a vertical Z-component receiver that receives primarily compressional wave information, and a P-component that receives primarily compressional wave components propagating in the sea. The vibrations received by the component of the detector X, Y are primarily a function of the converted transverse wave energy propagating radially and tangentially. Therefore, the premise of the research of transverse waves is that the data in the X, Y direction needs to be rotated to the radial direction and the tangential direction.

The effect of the R-T rotation depends on the accuracy of the polarization angle. In OBS exploration, OBS is freely put into the seabed after sea GPS positioning, is influenced by sea waves and the like, and the X, Y component azimuth of the seabed of the OBS is inconsistent with the azimuth when being put, even if a part of compass with the azimuth recorded during OBS design is influenced by a marine magnetic field, the azimuth recorded by the compass is inaccurate, so that a polarization angle cannot be accurately obtained, and data in X, Y directions cannot be accurately rotated to radial and tangential (R-T) directions.

Disclosure of Invention

In view of this, the embodiments of the present disclosure provide a method and a device for calculating R-T components of OBS seismic data, so as to solve the problem in the prior art that the polarization angle cannot be accurately obtained, and thus the data in the X, Y direction cannot be accurately rotated to the radial and tangential (R-T) directions.

In a first aspect of an embodiment of the present disclosure, a method for calculating an R-T component of OBS seismic data is provided, including:

acquiring seismic channel data of a plurality of seismic channels corresponding to a plurality of shot points in a preset radius range by taking a point of a target OBS vertically mapped to a sea level as a center, wherein the seismic channel data of each seismic channel comprises a first arrival signal;

projecting the amplitude polarity of the first arrival signal of the seismic channel data to a target coordinate system, and calculating the horizontal component azimuth angle of the target OBS, wherein the target coordinate system takes the OBS as an origin coordinate, the positive east direction of the origin coordinate is an X-axis positive direction, and the positive north direction of the origin coordinate is a Y-axis positive direction;

calculating a polarization angle of a seismic trace corresponding to each first-arrival signal based on the projection coordinates of each first-arrival signal and the horizontal component azimuth of the OBS;

and rotating the X-Y component in the horizontal component of the seismic channel data of the corresponding seismic channel based on each polarization angle to obtain the R-T component of the seismic channel data of the plurality of seismic channels.

In some embodiments, before projecting the amplitude polarity of the first arrival signal of the seismic trace data to the target coordinate system, the method further comprises:

and picking up the first arrival time of the first arrival signal of each seismic channel data, and acquiring the amplitude polarity of the first arrival signal at the first arrival time.

In some embodiments, the amplitude polarity of each of the first arrival signals includes an X-component amplitude polarity and a Y-component amplitude polarity;

the projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system, calculating the horizontal component azimuth of the target OBS, including:

calculating a first included angle between a projection coordinate of a first arrival signal with positive Y component amplitude polarity and an origin of the target coordinate system and an X axis of the target coordinate system to obtain a plurality of first included angles;

calculating a first included angle maximum value and a first included angle minimum value in the plurality of first included angles;

calculating a second included angle between a projection coordinate of a first arrival signal with positive X component amplitude polarity and an origin of the target coordinate system and a Y axis of the target coordinate system to obtain a plurality of second included angles;

calculating a second included angle maximum value and a second included angle minimum value in the plurality of second included angles;

And calculating the angle difference between the maximum value of the first included angle and the minimum value of the second included angle, or calculating the angle difference between the minimum value of the first included angle and the maximum value of the second included angle, so as to obtain the horizontal component azimuth angle of the target OBS.

In some embodiments, the calculating the polarization angle of the seismic trace corresponding to each first arrival signal based on the projected coordinates of each first arrival signal and the horizontal component azimuth of the OBS comprises:

calculating an offset azimuth of the seismic channel corresponding to each first arrival signal based on the projection coordinates of each first arrival signal;

the polarization angle of each seismic trace is calculated based on the horizontal component azimuth of the OBS and the offset azimuth of each seismic trace.

In some embodiments, the calculating the offset azimuth of the seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal includes:

substituting the projection coordinates of each first arrival signal and the original point coordinates of the target coordinate system into a first calculation formula, and calculating the offset azimuth of the seismic channel corresponding to each first arrival signal, wherein the first calculation formula is as follows:

wherein alpha represents the offset azimuth, x ₁ Representing the X-coordinate, y-coordinate representing the origin coordinate of the target coordinate system ₁ Representing the coordinate system representing the targetY coordinate, x of origin coordinate ₂ X-coordinate, y-coordinate representing projection coordinates of first arrival signal representing shot point ₂ And a Y coordinate representing the projection coordinates of the first arrival signal representing the shot.

In some embodiments, the calculating the polarization angle for each seismic trace based on the OBS horizontal component azimuth and the offset azimuth for each seismic trace comprises:

substituting the OBS horizontal component azimuth and the offset azimuth of each seismic trace into a second calculation formula, and calculating the polarization angle of each seismic trace, wherein the second calculation formula is as follows:

θ＝α-β-π/2

where θ represents the polarization angle of the seismic trace, α represents the offset azimuth, and β represents the azimuth of the horizontal component of the OBS.

In some embodiments, the rotating the X-Y component of the horizontal component of the trace data of the corresponding trace based on each polarization angle to obtain the R-T component of the trace data of the plurality of traces includes:

the polarization angle of each seismic channel and X and Y components in horizontal components of seismic channel data are led into a preset third calculation formula, the horizontal components are rotated into R-T components, and the third calculation formula is as follows:

wherein U is _X Representing the X component, U, of the horizontal components _Y Representing the Y component, U, of the horizontal components _R Representing the R component, U, of the R-T components _T Representing the T component of the R-T components, and θ representing the polarization angle of each seismic trace.

In a second aspect of embodiments of the present disclosure, there is provided an R-T component calculation apparatus for OBS seismic data, comprising:

the system comprises an acquisition module, a prediction module and a prediction module, wherein the acquisition module is used for acquiring seismic channel data of a plurality of seismic channels corresponding to a plurality of shot points in a preset radius range by taking a point of a target OBS (on-board) vertically mapped to a sea surface as a center, wherein the seismic channel data of each seismic channel comprises a first arrival signal;

the first calculation module is used for projecting the amplitude polarity of the first arrival signal of the seismic channel data to a target coordinate system and calculating the horizontal component azimuth angle of the target OBS, wherein the target coordinate system takes the OBS as an origin coordinate, the positive east direction of the origin coordinate is an X-axis positive direction, and the positive north direction of the origin coordinate is a Y-axis positive direction;

the second calculation module is used for calculating the polarization angle of the seismic channel corresponding to each first arrival signal based on the projection coordinates of each first arrival signal and the horizontal component azimuth of the OBS;

and the rotation calculation module is used for rotating X-Y components in horizontal components of the seismic channel data of the corresponding seismic channels based on each polarization angle to obtain R-T components of the seismic channel data of the plurality of seismic channels.

wherein alpha represents the offset azimuth, x ₁ Representing the X-coordinate, y-coordinate representing the origin coordinate of the target coordinate system ₁ Representing Y-coordinate, x representing the origin coordinate of the target coordinate system ₂ X-coordinate, y-coordinate representing projection coordinates of first arrival signal representing shot point ₂ And a Y coordinate representing the projection coordinates of the first arrival signal representing the shot.

θ＝α-β-π/2

In a third aspect of the disclosed embodiments, there is provided an electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the steps implemented by the processor when executing the computer program comprising:

θ＝α-β-π/2

In a fourth aspect of the disclosed embodiments, there is provided a computer-readable storage medium storing a computer program which, when executed by a processor, performs steps comprising:

θ＝α-β-π/2

Advantageous effects

Compared with the prior art, the beneficial effects of the embodiment of the disclosure at least comprise: calculating a horizontal component azimuth of the target OBS by projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system, and calculating a polarization angle of the seismic trace corresponding to each first arrival signal based on the projection coordinate of each first arrival signal; the polarization angle can be accurately obtained, and the R-T component is accurately obtained by rotating the data in the seismic trace from the horizontal component based on the polarization angle.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the drawings that are required for the embodiments or the description of the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are only the embodiments, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

FIG. 1 is a schematic illustration of one scenario of a method for R-T component computation of OBS seismic data provided in accordance with an embodiment of the present disclosure;

FIG. 2 is a flow chart of a second embodiment of a method for computing R-T components of OBS seismic data provided in accordance with an embodiment of the present disclosure;

FIG. 3a is a simplified schematic diagram of a shot, OBS and associated coordinate system of an R-T component calculation method for OBS seismic data provided in accordance with an embodiment of the present disclosure;

FIG. 3b is a schematic diagram of a distribution of first arrival signal polarity with respect to the X component (X axis) for an R-T component calculation method for OBS seismic data provided in accordance with an embodiment of the disclosure;

FIG. 3c is a schematic diagram of a distribution of first arrival signal polarity with respect to the Y component (Y axis) for an R-T component calculation method for OBS seismic data provided in accordance with an embodiment of the disclosure;

FIG. 3d is a schematic waveform diagram of the X, Y, R and T components of the first arrival signal polarity of an R-T component calculation method for OBS seismic data provided in accordance with an embodiment of the disclosure;

FIG. 4 is a flow chart of an embodiment III of another method for R-T component calculation of OBS seismic data provided in accordance with an embodiment of the present disclosure;

FIG. 5 is a simplified block diagram of an R-T component calculation device for OBS seismic data provided in accordance with an embodiment of the present disclosure;

fig. 6 is a schematic diagram of an electronic device provided according to an embodiment of the present disclosure.

Detailed Description

Embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While certain embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete. It should be understood that the drawings and embodiments of the present disclosure are for illustration purposes only and are not intended to limit the scope of the present disclosure.

It should be further noted that, for convenience of description, only a portion relevant to the present disclosure is shown in the drawings. Embodiments of the present disclosure and features of embodiments may be combined with each other without conflict.

It should be noted that the terms "first," "second," and the like in this disclosure are merely used to distinguish between different systems, devices, modules, or units and are not intended to limit the order or interdependence of functions performed by such systems, devices, modules, or units.

It should be noted that references to "one", "a plurality" and "a plurality" in this disclosure are intended to be illustrative rather than limiting, and those of ordinary skill in the art will appreciate that "one or more" is intended to be understood as "one or more" unless the context clearly indicates otherwise.

The names of messages or information interacted between the various devices in the embodiments of the present disclosure are for illustrative purposes only and are not intended to limit the scope of such messages or information.

The present disclosure will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.

Embodiment one:

fig. 1 is a schematic diagram of an application scenario of an R-T component calculation method of OBS seismic data according to a first embodiment of the present disclosure.

In the application scenario of fig. 1, first, a computing device 101 may acquire seismic trace data 102 of a plurality of seismic traces corresponding to a plurality of shots within a preset radius range with a point where a target OBS is vertically mapped to a sea level as a center, where the seismic trace data of each of the seismic traces includes a first arrival signal 103;

secondly, projecting the amplitude polarity of a first arrival signal 103 of the seismic trace data to a target coordinate system, and calculating a horizontal component azimuth 104 of the target OBS, wherein the target coordinate system takes the OBS as an origin coordinate, the positive east direction of the origin coordinate is an X-axis positive direction, and the positive north direction of the origin coordinate is a Y-axis positive direction;

again, the computing device 101 may calculate a polarization angle 106 of the seismic trace corresponding to each first arrival signal based on the projected coordinates 105 of each first arrival signal and the horizontal component azimuth 104 of the OBS;

Finally, the computing device 101 may rotate the X-Y component of the horizontal component of the trace data for the corresponding trace to the R-T direction based on each polarization angle 106, resulting in the R-T component 107 of the trace data for the plurality of traces.

The computing device 101 may be hardware or software. When the computing device is hardware, the computing device may be implemented as a distributed cluster formed by a plurality of servers or terminal devices, or may be implemented as a single server or a single terminal device. When the computing device is embodied as software, it may be installed in the hardware devices listed above. It may be implemented as a plurality of software or software modules, for example, for providing distributed services, or as a single software or software module. The present invention is not particularly limited herein.

It should be understood that the number of computing devices in fig. 1 is merely illustrative. There may be any number of computing devices, as desired for an implementation.

Embodiment two:

with continued reference to FIG. 2, a flow 200 of a second embodiment of a method of R-T component calculation of OBS seismic data in accordance with the present disclosure is shown. The method may be performed by the computing device 101 in fig. 1. The method for calculating the R-T component of the OBS seismic data comprises the following steps:

Step 201, obtaining seismic channel data of a plurality of seismic channels corresponding to a plurality of shots within a preset radius range by taking a point of a target OBS vertically mapped to a sea level as a center, wherein the seismic channel data of each seismic channel comprises a first arrival signal.

In some alternative implementations, the execution body of the R-T component calculation method of the OBS seismic data (such as the calculation device 101 shown in fig. 1) may connect to the target device through a wired connection or a wireless connection, and then acquire seismic trace data of a plurality of seismic traces corresponding to a plurality of shots within a preset radius range with a point where the target OBS is vertically mapped to the sea level as a center, where the seismic trace data of each of the seismic traces includes a first arrival signal.

The target OBS may refer to an OBS for the current calculation. In calculating relevant data for an OBS, first, it is necessary to determine seismic trace data for a plurality of seismic traces corresponding to a plurality of shots used for calculation. Each seismic trace data includes a first arrival signal. The first arrival signal may refer to a first segment of seismic trace data (waveform data) corresponding to each shot detected by the OBS. The preset radius range may be set as desired, for example, 10 meters, 20.3 meters, etc., and is not limited herein. The more the number of shots, the higher the accuracy of the subsequent calculation, and the number of the maps is set according to the accuracy requirement, which is not limited herein. However, when setting shots, the distribution of the shots of the seismic channels within the preset radius range is required to be as uniform as possible, and the distances of each shot within the range are as consistent as possible; secondly, the shot points in the radius range can be covered as comprehensively as possible. The shots are distributed in all directions, and the smaller the shot distance is, the higher the accuracy of horizontal component azimuth is.

It should be noted that the wireless connection may include, but is not limited to, 3G/4G/5G connection, wiFi connection, bluetooth connection, wiMAX connection, zigbee connection, UWB (ultra wideband) connection, and other now known or later developed wireless connection.

And 202, projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system, and calculating the horizontal component azimuth angle of the target OBS, wherein the target coordinate system takes the OBS as an origin coordinate, the positive east direction of the origin coordinate is an X-axis positive direction, and the positive north direction of the origin coordinate is a Y-axis positive direction.

In some embodiments as shown in fig. 3a, S1 and S2 may represent 2 shots, the upper sea level coordinate system, i.e. the point on the sea surface where the OBS is mapped perpendicularly to the sea level, is the origin coordinate, the forward direction of the origin coordinate is the E-axis, the north direction of the origin coordinate is the N-axis forward coordinate system, and the lower coordinate system may refer to the target coordinate system, i.e. the coordinate system below the sea level where the OBS is the origin coordinate, the forward direction of the origin coordinate is the X-axis, and the north direction of the origin coordinate is the Y-axis forward coordinate system. Obviously, the N axis corresponds to the Y axis, and the directions of the N axis and the Y axis are the same; the E axis corresponds to the X axis, and the directions of the E axis and the X axis are the same. The Z axis is a direction perpendicular to the sea level and is used to represent each direction, and there are no other uses in this disclosure and will not be described in any great detail. The amplitude polarity of the first-arrival signal may refer to the signal amplitude of the first-arrival signal as it arrives at the OBS, as determined based on the X-axis and Y-axis data. The amplitude polarity of the first arrival signal may be obtained by a detector in the OBS.

As shown in fig. 3b and 3c, projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system may refer to displaying the acquired amplitude polarity data of the first arrival signal of each seismic trace to the target coordinate system, so that the execution subject may calculate the horizontal component azimuth of the target OBS based on the following steps:

in the first step, the executing body may calculate a first included angle between a first line (i.e. a line between a coordinate point corresponding to any plus sign in fig. 3b and the origin in the graph) between the projection coordinates of the first arrival signal with positive Y component amplitude polarity and the origin of the target coordinate system, and the first axis of the target coordinate system, to obtain a plurality of first included angles. The projection coordinates of the first arrival signal may refer to coordinate points represented by either + or-numbers in fig. 3b or 3 c. The first included angle may refer to an angle between the first connecting line and the X-axis, and the magnitude of the first included angle may be calculated based on a clockwise (or counterclockwise) direction, from 0 degrees to 360 degrees.

And a second step, the executing body may calculate a first included angle maximum value and a first included angle minimum value in the plurality of first included angles.

And thirdly, the execution main body can calculate a second included angle between the projection coordinate of the first arrival signal with positive amplitude polarity of each X component and the origin of the target coordinate system and the Y axis of the target coordinate system to obtain a plurality of second included angles. The second included angle may refer to an angle of the first connection line with the Y-axis, which may be calculated based on a clockwise (or counter-clockwise) direction, from 0 degrees to 360 degrees. It should be noted that the first angle and the second angle may be calculated based on the clockwise (or counterclockwise) direction, but the first angle and the second angle are calculated in the same manner, i.e. calculated in the clockwise direction or calculated in the counterclockwise direction.

And step four, the executing body can calculate a second included angle maximum value and a second included angle minimum value in the plurality of second included angles.

And fifthly, the executing body may calculate an angle difference between the maximum value of the first included angle and the minimum value of the second included angle, or calculate an angle difference between the minimum value of the first included angle and the maximum value of the second included angle, so as to obtain a horizontal component azimuth angle of the target OBS.

Note that, the first arrival signal having positive amplitude polarity may be selected for calculation, or the projection coordinates of the first arrival signal having negative amplitude polarity may be selected for calculation, and the calculation results may be the same, and the first arrival signal and the projection coordinates may be selected as needed, and the calculation is not particularly limited.

Step 203, calculating the polarization angle of the seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal and the horizontal component azimuth of the OBS.

In some embodiments, the executing entity may calculate the polarization angle of the seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal and the horizontal component azimuth of the OBS. The polarization angle may be calculated based on the prior art, or may be calculated based on the following calculation method:

The first step, the executing body may substitute the projection coordinates of each first arrival signal and the origin coordinates of the target coordinate system into a first calculation formula, and calculate the offset azimuth of the seismic trace corresponding to each first arrival signal, where the first calculation formula is:

The second step, the executing body may substitute the azimuth of the horizontal component of the OBS and the azimuth of the offset of each seismic trace into a second calculation formula, and calculate the polarization angle of each seismic trace, where the second calculation formula is:

θ＝α-β-π/2

And 204, rotating X-Y components in horizontal components of the seismic channel data of the corresponding seismic channels based on each polarization angle to obtain R-T components of the seismic channel data of the plurality of seismic channels.

In some embodiments, the executing entity may rotate the X-Y component of the horizontal component of the trace data of the corresponding trace to the R-T direction by each polarization angle based on a manner in the prior art, to obtain the R-T components of the trace data of the plurality of traces.

In some optional implementations of some embodiments, the executing body may import the polarization angle of each seismic trace and the X component and the Y component in the horizontal component of the seismic trace data into a preset third calculation formula, and rotate the horizontal component into an R-T component, where the third calculation formula is:

The beneficial effects of one of the implementation manners of the above embodiments of the disclosure include at least: calculating a horizontal component azimuth of the target OBS by projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system, and calculating a polarization angle of the seismic trace corresponding to each first arrival signal based on the projection coordinate of each first arrival signal; the polarization angle based on which the data in the seismic trace is rotated from the horizontal component to the R-T component can be obtained accurately.

The correctness and effectiveness of the method is demonstrated below on the basis of a specific embodiment:

the following table is a built numerical model, a 2.5D HTI (Horizontal Transverse Isotropy, transverse isotropy) three-layer horizontal lamellar medium model is built, wave detection points are arranged on the seabed, namely within the range of 1000X1000m at the bottom of a first layer, the wave distance is 25m, the line distance is 25m, the surface layer of the model is blasted, the second layer is an HTI layer, multi-component data is received, and the X component azimuth of the wave detection point is in the direction of the right east.

Where ε may indicate the intensity of the longitudinal wave anisotropy, γ may indicate the intensity of the transverse wave anisotropy, δ may indicate the transition parameters of the longitudinal wave propagation along the symmetry axis and the vertical symmetry axis, α may indicate the HTI medium symmetry axis orientation, and Φ may indicate the tilt angle of the symmetry axis.

Referring to fig. 3d, fig. 3d includes a (1) th portion representing an X component in shot set data obtained by forward modeling of the three-layer horizontal layered medium model, a (2) th portion representing a Y component in shot set data obtained by forward modeling of the three-layer horizontal layered medium model, a (3) th portion representing a rotated R component, and a (4) th portion representing a rotated T component, wherein an abscissa of each portion is a seismic trace number, an ordinate is time, and waveforms from top to bottom in each portion are respectively corresponding to a direct wave, a reflected longitudinal wave, and a converted transverse wave in the drawing.

As can be seen from fig. 3d, the energy of the direct wave corresponding to the X, Y component is relatively strong, so that the display is relatively clear. While the R component corresponds to a stronger direct wave energy and the T component corresponds to a weaker direct wave energy (shown unclear).

The X, Y component corresponds to a stronger energy of the reflected longitudinal wave, so that the display is clearer. While the R component corresponds to a stronger reflected longitudinal wave energy and the T component corresponds to a weaker reflected longitudinal wave energy (not clearly shown).

The energy of the converted transverse wave corresponding to the X, Y component and the R, T component is relatively strong, so that the display is relatively clear.

The distribution accords with the characteristics of the HTI medium, so that the conversion effect is good.

Embodiment III:

with continued reference to FIG. 4, a flow 400 of a third embodiment of a method of R-T component calculation of OBS seismic data in accordance with the present disclosure is shown, which may be performed by the computing device 101 of FIG. 1. The R-T component calculation method of the OBS seismic data comprises the following steps:

step 401, obtaining seismic channel data of a plurality of seismic channels corresponding to a plurality of shots within a preset radius range by taking a point of a target OBS vertically mapped to a sea level as a center, wherein the seismic channel data of each seismic channel comprises a first arrival signal.

Step 402, picking up first arrival time of first arrival signals of each seismic channel data, and obtaining amplitude polarity of the first arrival signals at the first arrival time.

In some embodiments, the executing body may pick up a first arrival time of a first arrival signal of each of the seismic trace data, and obtain an amplitude polarity of the first arrival signal at the first arrival time. Because of the large signal quantity, it is necessary to first determine the first arrival time of the first arrival signal of each seismic trace data, and then determine the amplitude polarity of the first arrival signal corresponding to the first arrival time from the acquired signals based on the first arrival time.

And step 403, projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system, and calculating the horizontal component azimuth angle of the target OBS, wherein the target coordinate system takes the OBS as an origin coordinate, the positive east direction of the origin coordinate is an X-axis positive direction, and the positive north direction of the origin coordinate is a Y-axis positive direction.

Step 404, calculating an offset azimuth of the seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal.

Step 405, calculating a polarization angle of each seismic trace based on the horizontal component azimuth of the OBS and the offset azimuth of each seismic trace.

And step 406, rotating the X-Y component in the horizontal component of the seismic channel data of the corresponding seismic channel based on each polarization angle to obtain the R-T components of the seismic channel data of the plurality of seismic channels.

In some alternative implementations, the specific implementation of steps 401 and 403-406 and the technical effects thereof may refer to the steps in the second embodiment corresponding to fig. 2, which are not described herein.

All the above optional solutions may be combined arbitrarily to form an optional embodiment of the present application, which is not described here in detail.

Embodiment four:

The following are device embodiments of the present disclosure that may be used to perform method embodiments of the present disclosure. For details not disclosed in the embodiments of the apparatus of the present disclosure, please refer to the embodiments of the method of the present disclosure.

With further reference to FIG. 5, as an implementation of the method described above for each of the above figures, the present disclosure provides an embodiment of an R-T component calculation device for OBS seismic data, which corresponds to the embodiment described above with respect to FIG. 2.

As shown in fig. 5, the R-T component calculation apparatus 500 of the OBS seismic data of the present embodiment includes:

the acquiring module 501 is configured to acquire seismic trace data of a plurality of seismic traces corresponding to a plurality of shots within a preset radius range with a point of vertical mapping of a target OBS to a sea level as a center, where the seismic trace data of each seismic trace includes a first arrival signal;

the first calculating module 502 is configured to project an amplitude polarity of a first arrival signal of the seismic trace data to a target coordinate system, and calculate a horizontal component azimuth of the target OBS, where the target coordinate system uses the OBS as an origin coordinate, a forward direction of the origin coordinate is an X-axis forward direction, and a forward north direction of the origin coordinate is a Y-axis forward direction;

a second calculating module 503, configured to calculate a polarization angle of a seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal and the horizontal component azimuth of the OBS;

The rotation calculation module 504 is configured to rotate an X-Y component in the horizontal component of the trace data of the corresponding trace based on each polarization angle, so as to obtain R-T components of the trace data of the plurality of traces.

In some optional implementations of some embodiments, before projecting the amplitude polarity of the first arrival signal of the seismic trace data to the target coordinate system, the method further includes:

In some alternative implementations of some embodiments, the amplitude polarity of each of the first arrival signals includes an X-component amplitude polarity and a Y-component amplitude polarity;

In some optional implementations of some embodiments, the calculating a polarization angle of the seismic trace corresponding to each first arrival signal based on the projected coordinates of each first arrival signal and the horizontal component azimuth of the OBS includes:

In some optional implementations of some embodiments, the calculating an offset azimuth of the seismic trace corresponding to each first arrival signal based on the projection coordinates of each first arrival signal includes:

In some alternative implementations of some embodiments, the calculating the polarization angle for each seismic trace based on the OBS horizontal component azimuth and the offset azimuth for each seismic trace includes:

θ＝α-β-π/2

In some optional implementations of some embodiments, the rotating the X-Y component of the horizontal component of the trace data of the corresponding trace based on each polarization angle, to obtain the R-T component of the trace data of the plurality of traces, includes:

wherein U is _X Representing the X component, U, of the horizontal components _Y Representing the Y component, U, of the horizontal components _R Representing the R component, U, of the R-T components _T Representing the T component of the R-T components, and θ representing the polarization angle of each seismic trace. It will be appreciated that the modules described in the apparatus 500 correspond to the various steps in the method described with reference to fig. 2. Thus, the operations, features and resulting benefits described above with respect to the method are equally applicable to the apparatus 500 and the modules contained therein, and are not described in detail herein.

Fifth embodiment:

as shown in fig. 6, the electronic device 600 may include a processing means (e.g., a central processing unit, a graphics processor, etc.) 601, which may perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 602 or a program loaded from a storage means 608 into a Random Access Memory (RAM) 603. In the RAM 603, various programs and data required for the operation of the electronic apparatus 600 are also stored. The processing device 601, the ROM 602, and the RAM 603 are connected to each other through a bus 604. An input/output (I/O) interface 605 is also connected to bus 604.

In general, the following devices may be connected to the I/O interface 605: input devices 606 including, for example, a touch screen, touchpad, keyboard, mouse, camera, microphone, accelerometer, gyroscope, and the like; an output device 607 including, for example, a Liquid Crystal Display (LCD), a speaker, a vibrator, and the like; storage 608 including, for example, magnetic tape, hard disk, etc.; and a communication device 609. The communication means 609 may allow the electronic device 600 to communicate with other devices wirelessly or by wire to exchange data. While fig. 6 shows an electronic device 600 having various means, it is to be understood that not all of the illustrated means are required to be implemented or provided. More or fewer devices may be implemented or provided instead. Each block shown in fig. 6 may represent one device or a plurality of devices as needed.

In particular, according to the present embodiment, the process described above with reference to the flowcharts may be implemented as a computer software program. For example, the present embodiment includes a computer program product comprising a computer program embodied on a computer readable medium, the computer program containing program code for performing the method shown in the flowchart. In the present embodiment, the computer program can be downloaded and installed from a network through the communication means 609, or installed from the storage means 608, or installed from the ROM 602. When the computer program is executed by the processing means 601, the following steps may be performed:

θ＝α-β-π/2

It should be noted that, in this embodiment, the computer readable medium may be a computer readable signal medium or a computer readable storage medium, or any combination of the two. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples of the computer-readable storage medium may include, but are not limited to: an electrical connection having 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. In this embodiment, a computer-readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present embodiment, however, the computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, with the computer-readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: electrical wires, fiber optic cables, RF (radio frequency), and the like, or any suitable combination of the foregoing.

In some implementations, the clients, servers may communicate using any currently known or future developed network protocol, such as HTTP (HyperText Transfer Protocol ), and may be interconnected with any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network ("LAN"), a wide area network ("WAN"), the internet (e.g., the internet), and peer-to-peer networks (e.g., ad hoc peer-to-peer networks), as well as any currently known or future developed networks.

The computer readable medium may be embodied in the apparatus; or may exist alone without being incorporated into the electronic device. The computer readable medium carries one or more programs which, when executed by the electronic device, cause the electronic device to perform the steps of:

The computer program code for carrying out operations of the present embodiments may be written in one or more programming languages, including an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).

The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

The modules described in the present embodiment may be implemented by software or hardware. The described modules may also be provided in a processor, for example, as:

The device comprises an acquisition module, a first calculation module, a second calculation module and a conversion calculation module. For example, the acquisition module may also be described as "a module that acquires seismic trace data of a plurality of seismic traces corresponding to a plurality of shots".

The functions described above herein may be performed, at least in part, by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: a Field Programmable Gate Array (FPGA), an Application Specific Integrated Circuit (ASIC), an Application Specific Standard Product (ASSP), a system on a chip (SOC), a Complex Programmable Logic Device (CPLD), and the like.

The foregoing description is only of the preferred embodiments of the present disclosure and description of the principles of the technology being employed. It will be appreciated by those skilled in the art that the scope of the invention in the embodiments of the present disclosure is not limited to the specific combination of the above technical features, but encompasses other technical features formed by any combination of the above technical features or their equivalents without departing from the spirit of the invention. Such as the above-described features, are mutually substituted with (but not limited to) the features having similar functions disclosed in the embodiments of the present disclosure.

Claims

1. A method for computing an R-T component of OBS seismic data, comprising:

2. The method of claim 1, wherein projecting the amplitude polarity of the first arrival signal of the seismic trace data to a target coordinate system is preceded by:

3. The method of claim 1, wherein the amplitude polarity of each of the first arrival signals comprises an X-component amplitude polarity and a Y-component amplitude polarity;

4. The method of claim 1, wherein the calculating a polarization angle of the seismic trace corresponding to each first arrival signal based on the projected coordinates of each first arrival signal and the horizontal component azimuth of the OBS comprises:

5. The method of claim 4, wherein calculating the offset azimuth of the seismic trace corresponding to each first arrival signal based on the projected coordinates of each first arrival signal comprises:

6. The method of claim 4, wherein said calculating a polarization angle for each seismic trace based on the OBS horizontal component azimuth and an offset azimuth for each seismic trace comprises:

θ＝α-β-π/2

7. The method of claim 4, wherein rotating the X-Y component of the horizontal component of the trace data of the corresponding trace based on each polarization angle results in the R-T component of the trace data of the plurality of traces, comprising:

8. An R-T component calculation apparatus for OBS seismic data for offshore exploration, comprising:

9. An electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any of claims 1 to 7 when the computer program is executed.

10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method according to any one of claims 1 to 7.