Three-dimensional nine-component wave field separation method, device, equipment and medium

Abstract

The invention provides a three-dimensional nine-component wave field separation method, device and equipment and a medium, and belongs to the technical field of seismic data processing. The method comprises the following steps: acquiring three-dimensional nine-component seismic data; establishing an X-Y coordinate system; for any detection point of any wave source, a first angle is determined, the first angle is the angle of an included angle between a connecting line of the wave source and the detection point in the X-Y coordinate system and the X axis, and relation data is acquired based on the wave source type of the wave source, the relation data is used for indicating the relation between the first angle and the included angle between the R component positive direction and the X component positive direction, the included angle between the R component positive direction and the Y component positive direction, the included angle between the T component positive direction and the X component positive direction and the included angle between the T component positive direction and the Y component positive direction of seismic data acquired by a detection point of the wave source; and based on the first angle and the relational data, mapping the seismic data acquired by the detection point from the X-Y coordinate system to the R-T coordinate system, thereby realizing wave field separation of the three-dimensional nine-component seismic data.

Description

Technical Field

[0001] The present application relates to the technical field of seismic data processing, and in particular to a three-dimensional nine-component wave field separation method, device, equipment and medium. Background Art

[0002] When collecting multi-component seismic data, when the line connecting the shot point and the detection point obliquely intersects the survey line direction, the SV wave and SH wave propagating to the detection point will be projected on both the X and Y components. When preprocessing the seismic data, it is necessary to separate the SV wave and SH wave wave fields.

[0003] At present, the wave field separation technology for three-dimensional three-component seismic data has been developed to maturity. However, with the successful development of S-wave mechanical seismic sources, a large amount of three-dimensional nine-component seismic data has been collected, but the wave field separation technology for three-dimensional nine-component seismic data has not yet been developed to maturity. Summary of the invention

[0004] The embodiment of the present application provides a three-dimensional nine-component wavefield separation method, device, equipment and medium, which realizes the wavefield separation of three-dimensional nine-component seismic data. The technical solution is as follows:

[0005] In one aspect, a three-dimensional nine-component wavefield separation method is provided, the method comprising:

[0006] Acquire three-dimensional nine-component seismic data, wherein the three-dimensional nine-component seismic data includes seismic data corresponding to a P wave source, seismic data corresponding to an SV wave source excited along an X direction, and seismic data corresponding to an SH wave source excited along a Y direction;

[0007] Determine the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line;

[0008] Establish an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is a direction orthogonal to the direction of the survey line and satisfies the right-hand rule;

[0009] For any detection point of any wave source, a first angle is determined based on the positions of the wave source and the detection point in the XY coordinate system, the first angle being the angle between the line connecting the wave source and the detection point and the X-axis in the XY coordinate system; based on the wave source type of the wave source, relationship data is obtained, the relationship data being used to indicate the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

[0010] On the other hand, a three-dimensional nine-component wave field separation device is provided, the device comprising:

[0011] A first acquisition module is used to acquire three-dimensional nine-component seismic data, wherein the three-dimensional nine-component seismic data includes seismic data corresponding to a P wave source, seismic data corresponding to an SV wave source excited along an X direction, and seismic data corresponding to an SH wave source excited along a Y direction;

[0012] A determination module, used to determine the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line;

[0013] An establishment module is used to establish an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is a direction orthogonal to the direction of the survey line and satisfies the right-hand rule;

[0014] A separation module is used to determine a first angle for any detection point of any wave source based on the positions of the wave source and the detection point in the XY coordinate system, wherein the first angle is the angle between the line connecting the wave source and the detection point and the X-axis in the XY coordinate system, and to obtain relationship data based on the wave source type of the wave source, wherein the relationship data is used to indicate the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

[0015] On the other hand, a computer device is provided, comprising a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the three-dimensional nine-component wave field separation method as described in any of the above implementations.

[0016] On the other hand, a computer-readable storage medium is provided, wherein at least one program code is stored in the computer-readable storage medium, and the at least one program code is loaded and executed by a processor to implement the three-dimensional nine-component wave field separation method as described in any of the above implementations.

[0017] On the other hand, a computer program product is provided, the computer program product comprising at least one program code, the at least one program code being loaded and executed by a processor to implement the three-dimensional nine-component wave field separation method as described in any of the above implementations.

[0018] The beneficial effects of the technical solution provided by the embodiments of the present application include at least:

[0019] The embodiment of the present application provides a three-dimensional nine-component wave field separation method. It only needs to determine the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis, and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source can be obtained through relational data. Based on these angles, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system. The wave field separation method is not only simple, but also applicable to seismic data of various wave sources, and realizes the wave field separation of three-dimensional nine-component seismic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 It is a flow chart of a three-dimensional nine-component wave field separation method provided in an embodiment of the present application;

[0022] Figure 2 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0023] Figure 3 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0024] Figure 4 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0025] Figure 5 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0026] Figure 6 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0027] Figure 7 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0028] Figure 8 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0029] Fig. 9 is a schematic diagram of a survey line direction provided in an embodiment of the present application;

[0030] Fig.10 It is a flow chart of a three-dimensional nine-component wave field separation method provided in an embodiment of the present application;

[0031] Fig.11 It is a schematic diagram of the positive direction of the R component and the positive direction of the T component of seismic data of an SV wave source provided in an embodiment of the present application;

[0032] Fig.12 It is a schematic diagram of the positive direction of the R component and the positive direction of the T component of seismic data of an SH wave source provided in an embodiment of the present application;

[0033] Fig.13 is a schematic diagram of an X component of seismic data corresponding to a P wave source provided in an embodiment of the present application;

[0034] Fig.14 is a schematic diagram of a Y component of seismic data corresponding to a P wave source provided in an embodiment of the present application;

[0035] Fig.15 is a schematic diagram of an R component of seismic data corresponding to a P wave source provided in an embodiment of the present application;

[0036] Fig.16 is a schematic diagram of a T component of seismic data corresponding to a P wave source provided in an embodiment of the present application;

[0037] Fig.17 is a schematic diagram of an X component of seismic data corresponding to an SV wave source provided in an embodiment of the present application;

[0038] Fig.18 is a schematic diagram of a Y component of seismic data corresponding to an SV wave source provided in an embodiment of the present application;

[0039] Fig.19 is a schematic diagram of an R component of seismic data corresponding to an SV wave source provided in an embodiment of the present application;

[0040] Fig. 20 is a schematic diagram of a T component of seismic data corresponding to an SV wave source provided in an embodiment of the present application;

[0041] Fig.21 is a schematic diagram of an X component of seismic data corresponding to an SH wave source provided in an embodiment of the present application;

[0042] Fig. 22 is a schematic diagram of a Y component of seismic data corresponding to an SH wave source provided in an embodiment of the present application;

[0043] Fig.23 is a schematic diagram of an R component of seismic data corresponding to an SH wave source provided in an embodiment of the present application;

[0044] Fig.24 is a schematic diagram of a T component of seismic data corresponding to an SH wave source provided in an embodiment of the present application;

[0045] Fig.25 It is a structural schematic diagram of a three-dimensional nine-component wave field separation device provided in an embodiment of the present application;

[0046] Fig.26 is a schematic diagram of the structure of a terminal provided in an embodiment of the present application;

[0047] Fig. 27 It is a structural diagram of a server provided in an embodiment of the present application. DETAILED DESCRIPTION

[0048] In order to make the objectives, technical solutions and advantages of the present application clearer, the implementation methods of the present application will be further described in detail below with reference to the accompanying drawings.

[0049] The terms "first", "second", "third" and "fourth" etc. in the specification and claims of the present application and the drawings are used to distinguish different objects, rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices.

[0050] Figure 1is a flow chart of a three-dimensional nine-component wave field separation method provided in an embodiment of the present application. The present application embodiment is illustrated by taking a computer device as an example. Figure 1 , the method comprising:

[0051] 101. A computer device acquires three-dimensional nine-component seismic data, wherein the three-dimensional nine-component seismic data includes seismic data corresponding to a P-wave source, seismic data corresponding to an SV-wave source excited along an X-direction, and seismic data corresponding to an SH-wave source excited along a Y-direction.

[0052] Among them, P waves are longitudinal waves, SH waves and SV waves are transverse waves. If the vibration of the transverse wave occurs in a vertical plane passing through the direction of wave propagation, it is called an SV wave, and if it occurs in a horizontal plane, it is called an SH wave.

[0053] The seismic data corresponding to the P wave source is the seismic data collected at the detection point after the P wave source excites the P wave. The seismic data corresponding to the SV wave source exciting along the X direction is the seismic data collected at the detection point after the SV wave source exciting along the X direction excites the SV wave. The seismic data corresponding to the SH wave source exciting along the Y direction is the seismic data collected at the detection point after the SH wave source exciting along the Y direction excites the SH wave.

[0054] Among them, the seismic data corresponding to the P wave source includes data on the X component, data on the Y component and data on the Z component; the seismic data corresponding to the SV wave source excited along the X direction includes data on the X component, data on the Y component and data on the Z component; the seismic data corresponding to the SH wave source excited along the Y direction includes data on the X component, data on the Y component and data on the Z component. Therefore, the seismic data acquired in the above step 101 has a total of nine components, and the seismic data is called three-dimensional nine-component seismic data.

[0055] 102. The computer equipment determines the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line.

[0056] The survey line is a straight line with a series of detection points, and the unchanged observation line means that the detection points set on the survey line are evenly distributed. In some embodiments, the detection points are set with their own labels, and the survey line direction is the direction from the small detection point on the survey line to the large detection point. The computer device can determine the survey line direction based on the coordinates of any two detection points on the unchanged observation line.

[0057] For example, take two detection points A1 and A2 on the unchanged observation line, where the number of detection point A1 is smaller than that of detection point A2, the coordinates of detection point A1 are (gx1, gy1), and the coordinates of detection point A2 are (gx2, gy2). The difference between the horizontal coordinates of detection point A1 and detection point A2 is: dx = gx2-gx1; the difference between the vertical coordinates of detection point A1 and detection point A2 is: dy = gy2-gy1; the distance between detection point A1 and detection point A2 is:

[0058] When dx>0 and dy>0, such as Figure 2 As shown, the survey line direction angle = acos(|dx / da|), where acos is the arc cosine function. Figure 2 The direction indicated by the arrow in the coordinate system is the direction of the survey line. When dx<0 and dy>0, Figure 3 As shown in the figure, the survey line direction angle = π-aa, where aa = acos(|dx / da|). When dx < 0 and dy < 0, as Figure 4 As shown in the figure, the survey line direction angle = π + aa, where aa = acos (|dx / da|). When dx>0 and dy<0, as Figure 5 As shown in the figure, the survey line direction angle = 2π-aa, where aa = acos(|dx / da|). When dx = 0 and dy> 0, as Figure 6 As shown, the survey line direction angle = π / 2. When dx < 0 and dy = 0, as Figure 7 As shown, the survey line direction angle = π. When dx = 0 and dy < 0, as Figure 8 As shown, the survey line direction angle = 1.5π. When dx>0 and dy=0, as Fig. 9 As shown, the survey line direction angle = 0.

[0059] 103. The computer device establishes an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is the direction orthogonal to the direction of the survey line and satisfies the right-hand rule.

[0060] In the embodiment of the present application, the origin of the XY coordinate system is the position of the wave source. When the P wave source, the SV wave source excited along the X direction, and the SH wave source excited along the Y direction are located at the same position, only one XY coordinate system needs to be established. When the P wave source, the SV wave source excited along the X direction, and the SH wave source excited along the Y direction are located at different positions, it is necessary to establish an XY coordinate system for the P wave source based on the position of the P wave source, and process the seismic data corresponding to the P wave source based on the XY coordinate system; based on the position of the SV wave source excited along the X direction, an XY coordinate system is established for the SV wave source excited along the X direction, and the seismic data corresponding to the SV wave source excited along the X direction is processed based on the XY coordinate system; based on the position of the SH wave source excited along the Y direction, an XY coordinate system is established for the SH wave source excited along the Y direction, and the seismic data corresponding to the SH wave source excited along the Y direction is processed based on the XY coordinate system.

[0061] 104. The computer device determines, for any detection point of any wave source, a first angle based on positions of the wave source and the detection point in an XY coordinate system, the first angle being the angle between a line connecting the wave source and the detection point and an X-axis in the XY coordinate system, and obtains relationship data based on a wave source type of the wave source, the relationship data being used to indicate a relationship between the first angle and an angle between a positive direction of an R component and a positive direction of an X component, an angle between a positive direction of an R component and a positive direction of a Y component, an angle between a positive direction of a T component and a positive direction of an X component, and an angle between a positive direction of a T component and a positive direction of a Y component of seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point are mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

[0062] In the embodiment of the present application, the wave field separation of the seismic data can be achieved by mapping the seismic data from the XY coordinate system to the RT coordinate system. When mapping the seismic data from the XY coordinate system to the RT coordinate system, it is necessary to obtain the angle between the positive direction of the R component of the seismic data and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component. In the embodiment of the present application, the relational data can represent the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis and the angle between the positive direction of the R component of the seismic data obtained by the detection point of the wave source and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component. Therefore, it is only necessary to determine the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis to obtain the angles between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source. After obtaining the angles between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source, the seismic data acquired by the detection point of the wave source can be mapped from the XY coordinate system to the RT coordinate system to achieve wave field separation of the seismic data.

[0063] The three-dimensional nine-component wave field separation method provided in the embodiment of the present application only needs to determine the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis, and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source can be obtained through relational data. Based on these angles, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system. The wave field separation method is not only simple, but also applicable to seismic data of various wave sources, and realizes the wave field separation of three-dimensional nine-component seismic data.

[0064] In a possible implementation, acquiring relational data based on the wave source type of the wave source includes:

[0065] If the wave source is a P wave source, obtaining first relationship data;

[0066] The first relational data is represented as:

[0067]

[0068] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, R represents the R component of the seismic data, T represents the T component of the seismic data, θ is the first angle, The θ in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the (90°-θ) in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, the (90°+θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the θ in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, and cos represents the cosine function.

[0069] In a possible implementation, acquiring relational data based on the wave source type of the wave source includes:

[0070] If the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, the relationship data is acquired based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system.

[0071] In a possible implementation, if the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system, the relationship data is obtained, including:

[0072] If the wave source type of the wave source is an SV wave source excited along the X direction, and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system, obtain the second relationship data;

[0073] If the wave source type of the wave source is an SV wave source excited along the X direction, and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system, the third relationship data is obtained.

[0074] In a possible implementation, the second relation data is represented as:

[0075]

[0076] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, (180°-θ) in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction;

[0077] The third relation data is represented as:

[0078]

[0079] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the x direction, the θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the x direction, the (180°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the x direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the x direction.

[0080] In a possible implementation, if the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system, the relationship data is obtained, including:

[0081] If the wave source type of the wave source is an SH wave source excited along the Y direction, and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system, fourth relationship data is obtained;

[0082] If the wave source type of the wave source is an SH wave source excited along the Y direction, and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system, the fifth relationship data is obtained.

[0083] In a possible implementation, the fourth relationship data is represented as:

[0084]

[0085] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. (90°+θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction;

[0086] The fifth relation data is represented as:

[0087]

[0088] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the (180°-θ) in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction.

[0089] In a possible implementation, based on the first angle and the relationship data, after mapping the seismic data acquired by the detection point from the XY coordinate system to the RT coordinate system and obtaining the wave field separation result, the method further includes:

[0090] Correcting the first angle according to a preset angle interval to obtain a plurality of reference angles;

[0091] Obtain the wavefield separation result corresponding to each reference angle;

[0092] Acquire a wavefield separation quality parameter of a wavefield separation result corresponding to the first angle and a wavefield separation quality parameter of a wavefield separation result corresponding to each reference angle;

[0093] Based on the obtained wavefield separation quality parameter, a wavefield separation result with the highest wavefield separation quality is obtained.

[0094] All the above optional technical solutions can be arbitrarily combined to form optional embodiments of the present application, which will not be described one by one here.

[0095] Fig.10 is a flowchart of a three-dimensional nine-component wave field separation method provided in an embodiment of the present application. The embodiment of the present application is illustrated by taking a computer device as an example. Fig.10 , the method comprising:

[0096] 1001. A computer device acquires three-dimensional nine-component seismic data, where the three-dimensional nine-component seismic data includes seismic data corresponding to a P-wave source, seismic data corresponding to an SV-wave source excited along an X-direction, and seismic data corresponding to an SH-wave source excited along a Y-direction.

[0097] The step 1001 is the same as the above step 101 and will not be described in detail here.

[0098] 1002. The computer equipment determines the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line.

[0099] The above step 1002 is the same as the above step 102 and will not be described in detail here.

[0100] 1003. The computer device establishes an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is a direction orthogonal to the direction of the survey line and satisfies the right-hand rule.

[0101] The above step 1003 is the same as the above step 1002, and will not be described in detail here.

[0102] 1004. The computer device determines a first angle for any detection point of any wave source based on positions of the wave source and the detection point in the XY coordinate system. The first angle is an angle between a line connecting the wave source and the detection point in the XY coordinate system and an X-axis.

[0103] In the embodiment of the present application, the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis is less than 90 degrees. For example, the coordinates of the wave source are (0, 0), and the coordinates of the detection point are (gx3, gy3). The difference between the horizontal coordinates of the wave source and the detection point is: dx = gx3; the difference between the vertical coordinates of the wave source and the detection point is: dy = gy3; the distance between the wave source and the detection point is: The first angle = acos(|dx / da|), where acos is the arc cosine function.

[0104] 1005. If the wave source is a P wave source, the computer device obtains first relationship data, and maps the seismic data obtained by the detection point from the XY coordinate system to the RT coordinate system based on the first angle and the first relationship data to obtain a wave field separation result.

[0105] The first relationship data is used to represent the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component of the seismic data acquired by the detection point of the wave source, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component. In some embodiments, the first relationship data is an affine transformation matrix equation of the P wave source.

[0106] In some embodiments, the affine transformation matrix equation of the P-wave source can be expressed as:

[0107]

[0108] Where cos is the cosine function, (R, X) represents the angle between the positive direction of the R component and the positive direction of the X component, (R, Y) represents the angle between the positive direction of the R component and the positive direction of the Y component, (T, X) represents the angle between the positive direction of the T component and the positive direction of the X component, and (T, Y) represents the angle between the positive direction of the T component and the positive direction of the Y component; represents the X and Y components of the seismic data, Represents the R component and T component of seismic data.

[0109] For the P-wave source, the X-component polarity of the seismic data obtained at the detection points in the first and fourth quadrants is positive, and the Y-component polarity of the seismic data obtained at the detection points in the second and third quadrants is negative; the Y-component polarity of the seismic data obtained at the detection points in the first and second quadrants is positive, and the Y-component polarity of the seismic data obtained at the detection points in the third and fourth quadrants is negative.

[0110] Based on the azimuth relationship between the P-wave source and the detection point and the vibration direction of the P-wave source, the positive direction of the R component and the positive direction of the Y component of the seismic data of the P-wave source can be determined, wherein the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component can all be expressed by the first angle. The above formula can be rewritten as:

[0111]

[0112] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, R represents the R component of the seismic data, T represents the T component of the seismic data, θ is the first angle, The θ in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the (90°-θ) in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, the (90°+θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the θ in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, and cos represents the cosine function.

[0113] In some embodiments, the first relationship data is the above formula

[0114] 1006. If the wave source type of the wave source is an SV wave source excited along the X direction, and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system, the computer device obtains the second relationship data, and maps the seismic data obtained by the detection point from the XY coordinate system to the RT coordinate system based on the first angle and the second relationship data to obtain a wave field separation result.

[0115] The positive directions of the R component and the T component of the seismic data are related to the positional relationship between the wave source and the detection point and the vibration direction of the wave source. When the wave source is an SV wave source excited along the X direction, the positive directions of the R component and the T component of the seismic data obtained by the detection points in different quadrants are as follows: Fig.11 As shown. Among them, Fig.11 The S in the equation represents the wave source, and the G represents the detection point. Fig.11It can be seen that the positive direction of the R component of the seismic data obtained by the detection points in the first and third quadrants is the same as the positive direction of the T component, and the positive direction of the R component of the seismic data obtained by the detection points in the second and fourth quadrants is the same as the positive direction of the T component. If the wave source is an SV wave source excited along the X direction, and the detection point belongs to the first or third quadrant in the XY coordinate system, the second relationship data is obtained, and the seismic data obtained by the detection point is mapped from the XY coordinate system to the RT coordinate system based on the first angle and the second relationship data to obtain the wave field separation result. If the wave source is an SV wave source excited along the X direction, and the detection point belongs to the second or fourth quadrant in the XY coordinate system, the third relationship data is obtained, and the seismic data obtained by the detection point is mapped from the XY coordinate system to the RT coordinate system based on the first angle and the third relationship data to obtain the wave field separation result.

[0116] according to Fig.11 It can be seen that the angle between the positive direction of the R component of the seismic data and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component can all be represented by the first angle θ. For example, when the detection point is located in the first quadrant or the third quadrant, the angle between the positive direction of the R component of the seismic data and the positive direction of the X component is expressed as (90°-θ), the angle between the positive direction of the R component and the positive direction of the Y component is expressed as (180°-θ), the angle between the positive direction of the T component and the positive direction of the X component is expressed as θ, and the angle between the positive direction of the T component and the positive direction of the Y component is expressed as (90°-θ). Therefore, the second relationship data can be expressed as:

[0117]

[0118] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained from the detection points in the first quadrant or the third quadrant of the SV wave source excited along the X direction, the (180°-θ) in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained from the detection points in the first quadrant or the third quadrant of the SV wave source excited along the X direction, the θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained from the detection points in the first quadrant or the third quadrant of the SV wave source excited along the X direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained from the detection points in the first quadrant or the third quadrant of the SV wave source excited along the X direction.

[0119] 1007. If the wave source type of the wave source is an SV wave source excited along the X direction, and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system, the computer device obtains the third relationship data, and maps the seismic data obtained by the detection point from the XY coordinate system to the RT coordinate system based on the first angle and the third relationship data to obtain a wave field separation result.

[0120] according to Fig.11 It can be seen that when the detection point is located in the second or fourth quadrant, the angle between the positive direction of the R component of the seismic data and the positive direction of the X component is expressed as (90°-θ), the angle between the positive direction of the R component and the positive direction of the Y component is expressed as θ, the angle between the positive direction of the T component and the positive direction of the X component is expressed as (180°-θ), and the angle between the positive direction of the T component and the positive direction of the Y component is expressed as (90°-θ). Therefore, the third relationship data can be expressed as:

[0121]

[0122] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the (180°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction.

[0123] 1008. If the wave source type of the wave source is an SH wave source excited along the Y direction, and the detection point belongs to the first quadrant or the third quadrant of the XY coordinate system, the computer device obtains the fourth relationship data, and based on the first angle and the fourth relationship data, maps the seismic data obtained by the detection point from the XY coordinate system to the RT coordinate system to obtain the wave field separation result.

[0124] The positive directions of the R component and T component of the seismic data are related to the positional relationship between the wave source and the detection point and the vibration direction of the wave source. When the wave source is an SH wave source excited along the Y direction, the positive directions of the R component and T component of the seismic data obtained by the detection points in different quadrants are as follows: Fig.12 As shown, according to Fig.12 It can be seen that the positive direction of the R component of the seismic data obtained by the detection points in the first and third quadrants is the same as the positive direction of the T component, and the positive direction of the R component of the seismic data obtained by the detection points in the second and fourth quadrants is the same as the positive direction of the T component. If the wave source is an SH wave source excited in the Y direction, and the detection point belongs to the first or third quadrant in the XY coordinate system, the fourth relationship data is obtained, and based on the first angle and the fourth relationship data, the seismic data obtained by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain the wave field separation result. If the wave source is an SH wave source excited in the Y direction, and the detection point belongs to the second or fourth quadrant in the XY coordinate system, the fifth relationship data is obtained, and based on the first angle and the fifth relationship data, the seismic data obtained by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain the wave field separation result.

[0125] according to Fig.12 It can be known that the angle between the positive direction of the R component of the seismic data and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component can all be represented by the first angle θ. For example, when the detection point is located in the first quadrant or the third quadrant, the angle between the positive direction of the R component of the seismic data and the positive direction of the X component is expressed as (90°+θ), the angle between the positive direction of the R component and the positive direction of the Y component is expressed as θ, the angle between the positive direction of the T component and the positive direction of the X component is expressed as θ, and the angle between the positive direction of the T component and the positive direction of the Y component is expressed as (90°-θ). Therefore, the fourth relationship data can be expressed as:

[0126]

[0127] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°+θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, the θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, the θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction.

[0128] 1009. If the wave source type of the wave source is an SH wave source excited along the Y direction, and the detection point belongs to the second quadrant or the fourth quadrant of the XY coordinate system, the computer device obtains the fifth relationship data, and based on the first angle and the fifth relationship data, maps the seismic data obtained by the detection point from the XY coordinate system to the RT coordinate system to obtain the wave field separation result.

[0129] according to Fig.12 It can be seen that when the detection point is located in the second quadrant or the fourth quadrant, the angle between the positive direction of the R component of the seismic data and the positive direction of the X component is expressed as (90°-θ), the angle between the positive direction of the R component and the positive direction of the Y component is expressed as θ, the angle between the positive direction of the T component and the positive direction of the X component is expressed as (180°-θ), and the angle between the positive direction of the T component and the positive direction of the Y component is expressed as (90°-θ). Therefore, the fifth relationship data can be expressed as:

[0130]

[0131] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the (180°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction.

[0132] According to the above steps 1006 to 1009, if the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then the relationship data is obtained based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system. That is, the computer device obtains the relationship data based on the wave source type of the wave source, including: if the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then the relationship data is obtained based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system.

[0133] In some embodiments, the X component of the seismic data corresponding to the P wave source is as follows: Fig.13 As shown, the Y component of the seismic data corresponding to the P wave source is as follows Fig.14 As shown, after the seismic data corresponding to the P wave source is mapped from the XY coordinate system to the RT coordinate system, the R component of the seismic data corresponding to the P wave source is as follows: Fig.15 As shown, the T component of the seismic data corresponding to the P wave source is as follows Fig.16 As shown. Figures 13 to 16 It can be seen that after the seismic data corresponding to the P-wave source are separated by wave field, most of the energy is transferred to the R component, and a small amount of residue remains on the T component.

[0134] In some embodiments, the X component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Fig.17 As shown in the figure, the Y component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Fig.18 As shown, after mapping the seismic data corresponding to the SV wave source excited along the X direction from the XY coordinate system to the RT coordinate system, the R component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Fig.19 As shown in the figure, the T component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Fig. 20 As shown. Figures 17 to 20As can be seen from the box mark, after the wave field separation of the seismic data corresponding to the SV wave source excited along the X direction, most of the energy is transferred to the R component, which is the SV wave, and there is a small amount of SH wave residue on the T component.

[0135] In some embodiments, the X component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows: Fig.21 As shown, the Y component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows Fig. 22 As shown, after mapping the seismic data corresponding to the SH wave source excited along the Y direction from the XY coordinate system to the RT coordinate system, the R component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows: Fig.23 As shown, the T component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows: Fig.24 As shown. Figure 21 to Figure 24 As can be seen from the box mark, after the seismic data corresponding to the SV wave source excited along the X direction are separated by wave field, most of the energy is transferred to the R component, which is the SH wave, and there is a small amount of SV wave residue on the T component.

[0136] 1010. The computer device corrects the first angle according to a preset angle interval to obtain multiple reference angles, and acquires a wave field separation result corresponding to each reference angle.

[0137] When arranging the observation system of three-dimensional nine-component seismic data, human errors will be introduced. For example, the X component of the placement of the detection point does not strictly follow the direction of the survey line, so there is an error in the first angle.

[0138] In order to more accurately perform wavefield separation on seismic data, the embodiment of the present application adopts the maximum energy ratio angle scanning method to correct the wavefield separation result. The maximum energy ratio angle scanning method is to correct the first angle to obtain multiple reference angles to perform angle scanning on the first angle, and then determine the wavefield separation result with the largest energy ratio among the wavefield separation results corresponding to the first angle and the multiple reference angles.

[0139] In some embodiments, correcting the first angle to obtain multiple reference angles means: taking the first angle as a reference, gradually adjusting the angle towards both ends. For example, the first angle is 30°, and the preset angle interval is 0.5°. Taking the first angle as a reference, gradually adjusting the angle towards both ends means adding 0.5° to the first angle 30° to obtain a reference angle of 30.5°, and subtracting 0.5° from the first angle 30° to obtain a reference angle of 29.5°; further adjustment means adding 0.5° to 30.5° to obtain a reference angle of 31°, and subtracting 0.5° from 29.5° to obtain a reference angle of 29°, and so on.

[0140] To obtain the wave field separation result corresponding to each reference angle, the first angle is replaced by the reference angle, and the above steps 1005 to 1009 are performed again.

[0141] 1011. The computer device obtains a wavefield separation quality parameter of a wavefield separation result corresponding to the first angle and a wavefield separation quality parameter of a wavefield separation result corresponding to each reference angle, and obtains a wavefield separation result with the highest wavefield separation quality based on the obtained wavefield separation quality parameters.

[0142] The wavefield separation quality parameter is a parameter used to indicate the quality of the wavefield separation result. In some embodiments, the wavefield separation quality parameter is the energy ratio of the R component to the T component of the seismic data, and the wavefield separation result with the largest energy ratio is determined as the final wavefield separation result.

[0143] In some embodiments, the energy ratio of the R component to the T component of the seismic data is expressed as: Among them, top T represents the start time when the detection point receives the seismic wave, bot T represents the end time when the detection point receives the seismic wave, i represents the i-th moment, R(i) represents the R component of the seismic data at the i-th moment, and T(i) represents the T component of the seismic data at the i-th moment.

[0144] The three-dimensional nine-component wave field separation method provided in the embodiment of the present application only needs to determine the angle between the line connecting the wave source and the detection point in the XY coordinate system and the X-axis, and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source can be obtained through relational data. Based on these angles, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system. The wave field separation method is not only simple, but also applicable to seismic data of various wave sources, and realizes the wave field separation of three-dimensional nine-component seismic data.

[0145] In addition, the embodiment of the present application can also correct the first angle according to a preset angle interval to determine the optimal wave field separation result, thereby improving the accuracy of the wave field separation result.

[0146] Fig.25 is a schematic diagram of the structure of a three-dimensional nine-component wave field separation device provided in an embodiment of the present application, such as Fig.25 As shown, the device comprises:

[0147] The first acquisition module 2501 is used to acquire three-dimensional nine-component seismic data, where the three-dimensional nine-component seismic data includes seismic data corresponding to a P wave source, seismic data corresponding to an SV wave source excited along an X direction, and seismic data corresponding to an SH wave source excited along a Y direction;

[0148] A determination module 2502, for determining the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line;

[0149] Establishing module 2503, used to establish an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the survey line direction, and the positive direction of the Y component is the direction orthogonal to the survey line direction and satisfies the right-hand rule;

[0150] The separation module 2504 is used to determine, for any detection point of any wave source, a first angle based on the positions of the wave source and the detection point in the XY coordinate system, where the first angle is the angle between the line connecting the wave source and the detection point and the X-axis in the XY coordinate system, and to obtain relationship data based on the wave source type of the wave source, where the relationship data is used to indicate the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

[0151] In a possible implementation, the separation module 2504 is configured to obtain first relationship data if the wave source is a P-wave source;

[0152] The first relational data is represented as:

[0153]

[0154] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, R represents the R component of the seismic data, T represents the T component of the seismic data, θ is the first angle, The θ in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the (90°-θ) in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, the (90°+θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained at the detection point of the P-wave source, the θ in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained at the detection point of the P-wave source, and cos represents the cosine function.

[0155] In a possible implementation, the separation module 2504 is used to obtain relationship data based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system if the wave source type of the wave source is an SV wave source excited along the X direction or an SH wave source excited along the Y direction.

[0156] In one possible implementation, the separation module 2504 is used to obtain the second relationship data if the wave source type of the wave source is an SV wave source excited along the X direction and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system; and to obtain the third relationship data if the wave source type of the wave source is an SV wave source excited along the X direction and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system.

[0157] In a possible implementation, the second relation data is represented as:

[0158]

[0159] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, (180°-θ) in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained from the detection points in the first or third quadrant of the SV wave source excited along the X direction;

[0160] The third relation data is represented as:

[0161]

[0162] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the (180°-θ) in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction.

[0163] In one possible implementation, the separation module 2504 is used to obtain the fourth relationship data if the wave source type of the wave source is an SH wave source excited along the Y direction and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system; and to obtain the fifth relationship data if the wave source type of the wave source is an SH wave source excited along the Y direction and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system.

[0164] In a possible implementation, the fourth relationship data is represented as:

[0165]

[0166] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. (90°+θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction;

[0167] The fifth relation data is represented as:

[0168]

[0169] Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, the (180°-θ) in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction.

[0170] In a possible implementation, the device further includes:

[0171] A correction module, used for correcting the first angle according to a preset angle interval to obtain a plurality of reference angles;

[0172] The second acquisition module is used to obtain the wave field separation result corresponding to each reference angle;

[0173] A third acquisition module is used to acquire a wavefield separation quality parameter of a wavefield separation result corresponding to the first angle and a wavefield separation quality parameter of a wavefield separation result corresponding to each reference angle;

[0174] The fourth acquisition module is used to acquire a wavefield separation result with the highest wavefield separation quality based on the acquired wavefield separation quality parameter.

[0175] It should be noted that: the three-dimensional nine-component wave field separation device provided in the above embodiment only uses the division of the above functional modules as an example when performing wave field separation. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the computer device is divided into different functional modules to complete all or part of the functions described above. In addition, the three-dimensional nine-component wave field separation device provided in the above embodiment and the three-dimensional nine-component wave field separation method embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

[0176] In some embodiments, the computer device is provided as a terminal. Fig.26 The terminal 2600 includes a processor 2601 and a memory 2602 .

[0177] The processor 2601 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 2601 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 2601 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 2601 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 2601 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0178] The memory 2602 may include one or more computer-readable storage media, which may be non-transitory. The memory 2602 may also include a high-speed random access memory, and a non-volatile memory, such as one or more disk storage devices, flash memory storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 2602 is used to store at least one program code, which is used to be executed by the processor 2601 to implement the three-dimensional nine-component wave field separation method provided in the method embodiment of the present application.

[0179] In some embodiments, the terminal 2600 may also optionally include: a peripheral device interface 2603 and at least one peripheral device. The processor 2601, the memory 2602 and the peripheral device interface 2603 may be connected via a bus or a signal line. Each peripheral device may be connected to the peripheral device interface 2603 via a bus, a signal line or a circuit board. Specifically, the peripheral device includes: at least one of a radio frequency circuit 2604, a display screen 2605, a camera 2606, an audio circuit 2607, a positioning component 2608 and a power supply 2609.

[0180] The peripheral device interface 2603 may be used to connect at least one peripheral device related to I / O (Input / Output) to the processor 2601 and the memory 2602. In some embodiments, the processor 2601, the memory 2602, and the peripheral device interface 2603 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 2601, the memory 2602, and the peripheral device interface 2603 may be implemented on a separate chip or circuit board, which is not limited in this embodiment.

[0181] The display screen 2605 is used to display a UI (User Interface). The UI may include graphics, text, icons, videos, and any combination thereof. When the display screen 2605 is a touch display screen, the display screen 2605 also has the ability to collect touch signals on the surface or above the surface of the display screen 2605. The touch signal can be input to the processor 2601 as a control signal for processing. At this time, the display screen 2605 can also be used to provide virtual buttons and / or virtual keyboards, also known as soft buttons and / or soft keyboards. In some embodiments, the display screen 2605 can be one, and the front panel of the terminal 2600 is set; in other embodiments, the display screen 2605 can be at least two, which are respectively set on different surfaces of the terminal 2600 or are folded; in some other embodiments, the display screen 2605 can be a flexible display screen, which is set on the curved surface or folded surface of the terminal 2600. Even, the display screen 2605 can also be set to a non-rectangular irregular figure, that is, a special-shaped screen. The display screen 2605 can be made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

[0182] The power supply 2609 is used to power various components in the terminal 2600. The power supply 2609 can be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 2609 includes a rechargeable battery, the rechargeable battery can support wired charging or wireless charging. The rechargeable battery can also be used to support fast charging technology.

[0183] Those skilled in the art will understand that Fig.26 The structure shown in the figure does not constitute a limitation on the terminal 2600, and the terminal 2600 may include more or fewer components than shown in the figure, or combine certain components, or adopt a different component arrangement.

[0184] In some embodiments, the computer device is provided as a server. Fig. 27It is a structural diagram of a server provided in an embodiment of the present application. The server 2700 may have relatively large differences due to different configurations or performances, and may include one or more processors (Central Processing Units, CPU) 2701 and one or more memories 2702, wherein the memory 2702 stores at least one program code, and the at least one program code is loaded and executed by the processor 2701 to implement the methods provided in the above-mentioned various method embodiments. Of course, the server may also have components such as a wired or wireless network interface, a keyboard, and an input and output interface for input and output, and the server may also include other components for implementing device functions, which will not be described in detail here.

[0185] The server 2700 is used to execute the steps executed by the server in the above method embodiment.

[0186] The embodiment of the present application further provides a computer-readable storage medium, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the three-dimensional nine-component wave field separation method as described in any of the above implementations.

[0187] The embodiment of the present application further provides a computer program product, which includes at least one program code, and the at least one program code is loaded and executed by a processor to implement the three-dimensional nine-component wave field separation method as described in any of the above implementations.

[0188] In some embodiments, the computer program involved in the embodiments of the present application may be deployed and executed on a computer device, or on multiple computer devices located at one location, or on multiple computer devices distributed at multiple locations and interconnected by a communication network. Multiple computer devices distributed at multiple locations and interconnected by a communication network may constitute a blockchain system.

[0189] The above are only optional embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.

Claims

1. A three-dimensional nine-component wave field separation method, It is characterized in that The method comprises: Acquire three-dimensional nine-component seismic data, wherein the three-dimensional nine-component seismic data includes seismic data corresponding to a P wave source, seismic data corresponding to an SV wave source excited along an X direction, and seismic data corresponding to an SH wave source excited along a Y direction; Determine the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line; Establish an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is a direction orthogonal to the direction of the survey line and satisfies the right-hand rule; For any detection point of any wave source, a first angle is determined based on the positions of the wave source and the detection point in the XY coordinate system, the first angle being the angle between the line connecting the wave source and the detection point and the X-axis in the XY coordinate system; based on the wave source type of the wave source, relationship data is obtained, the relationship data being used to indicate the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

2. The method according to claim 1, It is characterized in that The acquiring of relationship data based on the wave source type of the wave source comprises: If the wave source is the P wave source, obtaining first relationship data; The first relational data is expressed as: Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, R represents the R component of the seismic data, T represents the T component of the seismic data, θ is the first angle, The θ in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data acquired by the detection point of the P-wave source, the (90°-θ) in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data acquired by the detection point of the P-wave source, the (90°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data acquired by the detection point of the P-wave source, the θ in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the P-wave source, and cos represents the cosine function.

3. The method according to claim 1, It is characterized in that The acquiring of relationship data based on the wave source type of the wave source comprises: If the wave source type of the wave source is the SV wave source excited along the X direction or the SH wave source excited along the Y direction, relationship data is acquired based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system.

4. The method according to claim 3, It is characterized in that If the wave source type of the wave source is the SV wave source excited along the X direction or the SH wave source excited along the Y direction, acquiring relationship data based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system, includes: If the wave source type of the wave source is the SV wave source excited along the X direction, and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system, acquiring second relationship data; If the wave source type of the wave source is the SV wave source excited along the X direction, and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system, third relationship data is acquired.

5. The method according to claim 4, It is characterized in that The second relational data is expressed as: Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X direction, (180°-θ) in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X direction, θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X direction; The third relational data is represented as: Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, the (180°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X direction.

6. The method according to claim 3, It is characterized in that If the wave source type of the wave source is the SV wave source excited along the X direction or the SH wave source excited along the Y direction, acquiring relationship data based on the wave source type of the wave source and the quadrant to which the detection point belongs in the XY coordinate system, includes: If the wave source type of the wave source is the SH wave source excited along the Y direction, and the detection point belongs to the first quadrant or the third quadrant in the XY coordinate system, acquiring fourth relationship data; If the wave source type of the wave source is the SH wave source excited along the Y direction, and the detection point belongs to the second quadrant or the fourth quadrant in the XY coordinate system, fifth relationship data is acquired.

7. The method according to claim 6, It is characterized in that The fourth relational data is expressed as: Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, the θ in the second row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, the θ in the first row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction, and the (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y direction; The fifth relational data is expressed as: Wherein, X represents the X component of the seismic data, Y represents the Y component of the seismic data, θ is the first angle, cos represents the cosine function, R represents the R component of the seismic data, and T represents the T component of the seismic data. The (90°-θ) in the first row and the first column represents the angle between the positive direction of the R component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, θ in the first row and the second column represents the angle between the positive direction of the R component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, (180°-θ) in the second row and the first column represents the angle between the positive direction of the T component and the positive direction of the X component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction, and (90°-θ) in the second row and the second column represents the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data obtained by the detection points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y direction.

8. The method according to claim 1, It is characterized in that After mapping the seismic data acquired by the detection point from the XY coordinate system to the RT coordinate system based on the first angle and the relationship data to obtain a wave field separation result, the method further includes: Correcting the first angle according to a preset angle interval to obtain a plurality of reference angles; Obtain the wavefield separation result corresponding to each reference angle; Acquire a wavefield separation quality parameter of the wavefield separation result corresponding to the first angle and a wavefield separation quality parameter of the wavefield separation result corresponding to each reference angle; Based on the obtained wavefield separation quality parameter, a wavefield separation result with the highest wavefield separation quality is obtained.

9. A three-dimensional nine-component wave field separation device, It is characterized in that The device comprises: A first acquisition module is used to acquire three-dimensional nine-component seismic data, wherein the three-dimensional nine-component seismic data includes seismic data corresponding to a P wave source, seismic data corresponding to an SV wave source excited along an X direction, and seismic data corresponding to an SH wave source excited along a Y direction; A determination module, used to determine the direction of the survey line based on the coordinates of any two detection points on the unchanged observation line; An establishment module is used to establish an XY coordinate system, wherein the origin of the XY coordinate system is the position of the wave source, the positive direction of the X component is the direction of the survey line, and the positive direction of the Y component is a direction orthogonal to the direction of the survey line and satisfies the right-hand rule; A separation module is used to determine a first angle for any detection point of any wave source based on the positions of the wave source and the detection point in the XY coordinate system, wherein the first angle is the angle between the line connecting the wave source and the detection point and the X-axis in the XY coordinate system, and to obtain relationship data based on the wave source type of the wave source, wherein the relationship data is used to indicate the relationship between the first angle and the angle between the positive direction of the R component and the positive direction of the X component, the angle between the positive direction of the R component and the positive direction of the Y component, the angle between the positive direction of the T component and the positive direction of the X component, and the angle between the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the detection point of the wave source; based on the first angle and the relationship data, the seismic data acquired by the detection point is mapped from the XY coordinate system to the RT coordinate system to obtain a wave field separation result.

10. A computer device, It is characterized in that The computer device comprises a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the three-dimensional nine-component wave field separation method according to any one of claims 1 to 8.

11. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the three-dimensional nine-component wave field separation method according to any one of claims 1 to 8.

Images