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

Abstract

The application provides a three-dimensional nine-component wave field separation method, device, equipment and 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; determining a first angle for any wave source and any geophone, the first angle being the angle between the connecting line of the wave source and the geophone and the X axis in the X-Y coordinate system; acquiring relationship data based on the wave source type of the wave source, the relationship data being used to indicate the relationship between the first angle and the angles between the positive direction of the R component and the positive direction of the X component, the positive direction of the R component and the positive direction of the Y component, the positive direction of the T component and the positive direction of the X component, and the positive direction of the T component and the positive direction of the Y component of the seismic data acquired by the geophone of the wave source; and mapping the seismic data acquired by the geophone from the X-Y coordinate system to an R-T coordinate system based on the first angle and the relationship data, thereby realizing the wave field separation of the three-dimensional nine-component seismic data.

Description

Technical Field

[0001] This application relates to the field of seismic data processing technology, and in particular to a three-dimensional nine-component wavefield separation method, apparatus, equipment and medium. Background Technology

[0002] When acquiring multi-component seismic data, if the line connecting the shot point and the receiver point is oblique to the direction of the survey line, the SV wave and SH wave propagating to the receiver point will have projections on both the X and Y components. When preprocessing the seismic data, it is necessary to separate the wavefields of the SV wave and SH wave.

[0003] Currently, wavefield separation technology has been well developed for three-dimensional, three-component seismic data. However, with the successful development of S-wave mechanical seismic sources, a large amount of three-dimensional, nine-component seismic data has been acquired, but wavefield separation technology for this type of data has not yet been fully developed. Summary of the Invention

[0004] This application provides a method, apparatus, device, and medium for three-dimensional nine-component wavefield separation, achieving wavefield separation of three-dimensional nine-component seismic data. The technical solution is as follows:

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

[0006] Acquire three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X direction, and seismic data corresponding to SH-wave sources excited along the Y direction.

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

[0008] Establish an XY coordinate system, where 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.

[0009] For any receiver point of any wave source, based on the positions of the wave source and the receiver point in the XY coordinate system, a first angle is determined. The first angle is the angle between the line connecting the wave source and the receiver point in the XY coordinate system and the X-axis. Based on the wave source type, relational data is obtained. The relational data indicates the relationship between the first angle and the angles between the positive R-component and the positive X-component, the positive R-component and the positive Y-component, the positive T-component and the positive X-component, and the positive T-component and the positive Y-component of the seismic data acquired by the receiver point of the wave source. Based on the first angle and the relational data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

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

[0011] The first acquisition module is used to acquire three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X direction, and seismic data corresponding to SH-wave sources excited along the Y direction.

[0012] The determination module is used to determine the direction of the survey line based on the coordinates of any two receiver points on the unchanging observation line;

[0013] A module is established 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 the direction orthogonal to the direction of the survey line and satisfies the right-hand rule.

[0014] The separation module is used to determine a first angle for any receiver point of any wave source, based on the positions of the wave source and the receiver point in the XY coordinate system. The first angle is the angle between the line connecting the wave source and the receiver point in the XY coordinate system and the X-axis. Based on the wave source type, relational data is obtained. The relational data indicates the relationship between the first angle and the angles between the positive R-component and the positive X-component, the positive R-component and the positive Y-component, the positive T-component and the positive X-component, and the positive T-component and the positive Y-component of the seismic data acquired by the receiver point of the wave source. Based on the first angle and the relational data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0015] On the other hand, a computer device is provided, the computer device including a processor and a memory, the memory storing at least one piece of program code, the at least one piece of program code being loaded and executed by the processor to implement the three-dimensional nine-component wavefield 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 piece of program code is stored in the computer-readable storage medium, the at least one piece of program code being loaded and executed by a processor to implement the three-dimensional nine-component wavefield 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 piece of program code, the at least one piece of program code being loaded and executed by a processor to implement the three-dimensional nine-component wavefield separation method as described in any of the above implementations.

[0018] The beneficial effects of the technical solutions provided in this application include at least the following:

[0019] This application provides a three-dimensional nine-component wavefield separation method. By simply determining the angle between the line connecting the wave source and the receiver in the XY coordinate system and the X-axis, the angles between the positive directions of the R-component and X-component, the positive directions of the R-component and Y-component, the positive directions of the T-component and X-component, and the positive directions of the T-component and Y-component can be obtained from the relational data. Based on these angles, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system. This wavefield separation method is not only simple but also applicable to seismic data from various wave sources, achieving wavefield separation of three-dimensional nine-component seismic data. Attached Figure Description

[0020] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0021] Figure 1 This is a flowchart of a three-dimensional nine-component wavefield separation method provided in an embodiment of this application;

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

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

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

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

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

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

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

[0029] Figure 9 This is a schematic diagram of a survey line direction provided in an embodiment of this application;

[0030] Figure 10 This is a flowchart of a three-dimensional nine-component wavefield separation method provided in an embodiment of this application;

[0031] Figure 11 This is a schematic diagram of the positive R component and positive T component of seismic data from an SV wave source provided in an embodiment of this application;

[0032] Figure 12 This is a schematic diagram of the positive R component and positive T component of seismic data from an SH wave source provided in an embodiment of this application;

[0033] Figure 13 This is a schematic diagram of the X component of seismic data corresponding to a P-wave source provided in an embodiment of this application;

[0034] Figure 14 This is a schematic diagram of the Y component of seismic data corresponding to a P-wave source provided in an embodiment of this application;

[0035] Figure 15 This is a schematic diagram of the R component of seismic data corresponding to a P-wave source provided in an embodiment of this application;

[0036] Figure 16 This is a schematic diagram of the T component of seismic data corresponding to a P-wave source provided in an embodiment of this application;

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

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

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

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

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

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

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

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

[0045] Figure 25 This is a schematic diagram of the structure of a three-dimensional nine-component wave field separation device provided in an embodiment of this application;

[0046] Figure 26 This is a schematic diagram of the structure of a terminal provided in an embodiment of this application;

[0047] Figure 27 This is a schematic diagram of the structure of a server provided in an embodiment of this application. Detailed Implementation

[0048] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0049] The terms "first," "second," "third," and "fourth," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish different objects, not to describe a specific order. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or apparatuses.

[0050] Figure 1This is a flowchart illustrating a three-dimensional nine-component wavefield separation method provided in an embodiment of this application. This embodiment uses a computer device as the execution subject for illustrative purposes. See also... Figure 1 The method includes:

[0051] 101. The computer equipment acquires three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X-direction, and seismic data corresponding to SH-wave sources excited along the Y-direction.

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

[0053] The seismic data corresponding to a P-wave source is the seismic data acquired by the receiver after the P-wave is generated by the P-wave source. The seismic data corresponding to an SV-wave source generated along the X-direction is the seismic data acquired by the receiver after the SV-wave is generated by the SV-wave source along the X-direction. The seismic data corresponding to an SH-wave source generated along the Y-direction is the seismic data acquired by the receiver after the SH-wave is generated by the SH-wave source along the Y-direction.

[0054] The seismic data corresponding to the P-wave source includes data in the X, Y, and Z components; the seismic data corresponding to the SV-wave source excited along the X direction includes data in the X, Y, and Z components; and the seismic data corresponding to the SH-wave source excited along the Y direction includes data in the X, Y, and Z components. Therefore, the seismic data obtained in step 101 has a total of nine components, and thus 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 receiver points on the unchanging observation line.

[0056] In this context, a survey line is a straight line with a series of receiver points, and a constant observation line refers to a survey line where receiver points are evenly distributed. In some embodiments, receiver points are assigned their own numbers, and the survey line direction is the direction from smaller receiver points to larger receiver points on the survey line. Computer equipment can determine the survey line direction based on the coordinates of any two receiver points on the constant observation line.

[0057] For example, consider two geophone points A1 and A2 on an unchanging observation line, where the label of geophone point A1 is less than that of geophone point A2. The coordinates of geophone point A1 are (gx1, gy1), and the coordinates of geophone point A2 are (gx2, gy2). The difference between the x-coordinates of geophone points A1 and A2 is: dx = gx2 - gx1; the difference between the y-coordinates of geophone points A1 and A2 is: dy = gy2 - gy1; the distance between geophone points A1 and 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 inverse cosine function. Figure 2 In the coordinate system, the arrow indicates the direction of the survey line. When dx < 0 and dy > 0, as... Figure 3 As shown, the survey line direction angle = π - aa, where aa = acos(|dx / da|). When dx < 0 and dy < 0, as... Figure 4 As shown, the survey direction angle = π + aa, where aa = acos(|dx / da|). When dx > 0 and dy < 0, as... Figure 5 As shown, the survey direction angle = 2π - aa, where aa = acos(|dx / da|). When dx = 0 and dy > 0, as... Figure 6 As shown, the direction angle of the survey line is π / 2. When dx < 0 and dy = 0, as... Figure 7 As shown, the direction angle of the survey line is π. 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... Figure 9 As shown, the survey line direction angle is 0.

[0059] 103. The computer equipment establishes an XY coordinate system, where 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 this embodiment, the origin of the XY coordinate system is the location 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 location, 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 locations, an XY coordinate system needs to be established for the P-wave source based on its location, and the seismic data corresponding to the P-wave source needs to be processed based on this XY coordinate system; an XY coordinate system needs to be established for the SV-wave source excited along the X-direction based on its location, and the seismic data corresponding to the SV-wave source excited along the X-direction needs to be processed based on this XY coordinate system; and an XY coordinate system needs to be established for the SH-wave source excited along the Y-direction based on its location, and the seismic data corresponding to the SH-wave source excited along the Y-direction needs to be processed based on this XY coordinate system.

[0061] 104. For any receiver point of any wave source, the computer equipment determines a first angle based on the positions of the wave source and the receiver point in the XY coordinate system. This first angle is the angle between the line connecting the wave source and the receiver point in the XY coordinate system and the X-axis. Based on the wave source type, relational data is obtained. This relational data indicates the relationship between the first angle and the angles between the positive directions of the R and X components, the positive directions of the R and Y components, the positive directions of the T and X components, and the positive directions of the T and Y components of the seismic data acquired by the receiver point. Based on the first angle and the relational data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0062] In this embodiment, wavefield separation of seismic data can be achieved by mapping seismic data from the XY coordinate system to the RT coordinate system. When mapping seismic data from the XY coordinate system to the RT coordinate system, it is necessary to obtain the angles between the positive directions of the R component and the positive directions of the X component, the R component and the positive directions of the Y component, the T component and the positive directions of the X component, and the T component and the positive directions of the Y component in the seismic data acquired by the receiver in the XY coordinate system. In this embodiment, the relational data can represent the relationship between the angle between the line connecting the wave source and the receiver in the XY coordinate system and the X-axis, and the angles between the positive directions of the R component and the positive directions of the X component, the R component and the positive directions of the Y component, the T component and the positive directions of the X component, and the T component and the positive directions of the Y component in the seismic data acquired by the receiver in the XY coordinate system. Therefore, by determining the angle between the line connecting the wave source and the receiver in the XY coordinate system and the X-axis, we can obtain the angles between the positive directions of the R, Y, T, and X components of the seismic data acquired by the receiver from the wave source. After obtaining these angles, the seismic data can be mapped from the XY coordinate system to the RT coordinate system, achieving wavefield separation of the seismic data.

[0063] The three-dimensional nine-component wavefield separation method provided in this application only requires determining the angle between the line connecting the wave source and the receiver in the XY coordinate system and the X-axis. Then, through relational data, the angles between the positive directions of the R-component and X-component, the R-component and Y-component, the T-component and X-component, and the T-component and Y-component of the seismic data acquired by the receiver from the wave source can be obtained. Based on these angles, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system. This wavefield separation method is not only simple but also applicable to seismic data from various wave sources, achieving wavefield separation of three-dimensional nine-component seismic data.

[0064] In one possible implementation, relational data is obtained based on the wave source type, including:

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

[0066] The first relational data is represented as follows:

[0067]

[0068] Where 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, and θ is the first angle. In the first row and first column, θ represents the angle between the positive directions of the R component and the X component of the seismic data acquired from the receiver of the P-wave source. In the first row and second column, (90°-θ) represents the angle between the positive directions of the R component and the Y component of the seismic data acquired from the receiver of the P-wave source. In the second row and first column, (90°+θ) represents the angle between the positive directions of the T component and the X component of the seismic data acquired from the receiver of the P-wave source. In the second row and second column, θ represents the angle between the positive directions of the T component and the Y component of the seismic data acquired from the receiver of the P-wave source. cos represents the cosine function.

[0069] In one possible implementation, relational data is obtained based on the wave source type, including:

[0070] If the wave source type 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 and the quadrant to which the receiver point belongs in the XY coordinate system.

[0071] In one possible implementation, if the wave source type 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 and the quadrant to which the receiver point belongs in the XY coordinate system, relational data is obtained, including:

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

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

[0074] In one possible implementation, the second relational data is represented as:

[0075]

[0076] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the first row and second column, (180°-θ) represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction.

[0077] Third-party relation data is represented as follows:

[0078]

[0079] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive directions of the R and X components of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the x-direction. In the first row and second column, θ represents the angle between the positive directions of the R and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the x-direction. In the second row and first column, (180°-θ) represents the angle between the positive directions of the T and X components of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the x-direction. In the second row and second column, (90°-θ) represents the angle between the positive directions of the T and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the x-direction.

[0080] In one possible implementation, if the wave source type 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 and the quadrant to which the receiver point belongs in the XY coordinate system, relational data is obtained, including:

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

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

[0083] In one possible implementation, the fourth relation data is represented as:

[0084]

[0085] Where 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. In the first row and first column, (90°+θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the first row and second column, θ represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction.

[0086] The fifth relation data is represented as:

[0087]

[0088] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive directions of the R and X components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the first row and second column, θ represents the angle between the positive directions of the R and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the second row and first column, (180°-θ) represents the angle between the positive directions of the T and X components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the second row and second column, (90°-θ) represents the angle between the positive directions of the T and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction.

[0089] In one possible implementation, after mapping the seismic data acquired from the receiver points from the XY coordinate system to the RT coordinate system based on the first angle and relationship data to obtain the wavefield separation results, the method further includes:

[0090] The first angle is corrected according to the preset angle interval to obtain multiple reference angles;

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

[0092] Obtain the wavefield separation quality parameters of the wavefield separation result corresponding to the first angle and the wavefield separation quality parameters of the wavefield separation result corresponding to each reference angle;

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

[0094] All of the above-mentioned optional technical solutions can be combined in any way to form the optional embodiments of this application, and will not be described in detail here.

[0095] Figure 10 This is a flowchart illustrating a three-dimensional nine-component wavefield separation method provided in this application embodiment. This application embodiment uses a computer device as the execution subject for illustrative purposes. See [link to relevant documentation]. Figure 10 The method includes:

[0096] 1001. The computer equipment acquires three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X-direction, and seismic data corresponding to SH-wave sources excited along the Y-direction.

[0097] Step 1001 is the same as step 101 above, 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 receiver points on the unchanging observation line.

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

[0100] 1003. The computer equipment establishes an XY coordinate system, where 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.

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

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

[0103] In this embodiment, the angle between the line connecting the wave source and the receiver 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 receiver are (gx3, gy3). The difference in the abscissa between the wave source and the receiver is: dx = gx3; the difference in the ordinate between the wave source and the receiver is: dy = gy3; the distance between the wave source and the receiver is: The first angle is equal to acos(|dx / da|). Here, acos is the inverse cosine function.

[0104] 1005. If the wave source is a P-wave source, the computer device acquires the first relational data, and based on the first angle and the first relational data, maps the seismic data acquired by the receiver point from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0105] The first relational data is used to represent the relationship between the first angle and the angles between the positive directions of the R and X components, the positive directions of the R and Y components, the positive directions of the T and X components, and the positive directions of the T and Y components of the seismic data acquired from the receiver point of the wave source. In some embodiments, the first relational data is the 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 directions of the R component and the X component, (R, Y) represents the angle between the positive directions of the R component and the Y component, (T, X) represents the angle between the positive directions of the T component and the X component, and (T, Y) represents the angle between the positive directions of the T component and the Y component. Represents the X and Y components of seismic data. This represents the R and T components of the seismic data.

[0109] For a P-wave source, the X component of seismic data acquired from receivers in the first and fourth quadrants has positive polarity, while the Y component of seismic data acquired from receivers in the second and third quadrants has negative polarity; the Y component of seismic data acquired from receivers in the first and second quadrants has positive polarity, while the Y component of seismic data acquired from receivers in the third and fourth quadrants has negative polarity.

[0110] Based on the azimuth relationship between the P-wave source and the receiver and the vibration direction of the P-wave source, the positive directions of the R-component and Y-component of the seismic data from the P-wave source can be determined. The angles between the positive directions of the R-component and X-component, the R-component and Y-component, the positive directions of the T-component and X-component, and the positive directions of the T-component and Y-component can all be represented by a first angle. The above formula... It can be rewritten as:

[0111]

[0112] Where 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, and θ is the first angle. In the first row and first column, θ represents the angle between the positive directions of the R component and the X component of the seismic data acquired from the receiver of the P-wave source. In the first row and second column, (90°-θ) represents the angle between the positive directions of the R component and the Y component of the seismic data acquired from the receiver of the P-wave source. In the second row and first column, (90°+θ) represents the angle between the positive directions of the T component and the X component of the seismic data acquired from the receiver of the P-wave source. In the second row and second column, θ represents the angle between the positive directions of the T component and the Y component of the seismic data acquired from the receiver of the P-wave source. cos represents the cosine function.

[0113] In some embodiments, the first relational data is the formula described above.

[0114] 1006. If the wave source type is an SV wave source excited along the X direction, and the receiver point belongs to the first or third quadrant in the XY coordinate system, the computer device acquires the second relational data. Based on the first angle and the second relational data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0115] The positive directions of the R-component and T-component of seismic data are related to the positional relationship between the wave source and receiver, as well as 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 T-component of the seismic data acquired from receivers in different quadrants are as follows: Figure 11 As shown. Among them, Figure 11 In this context, S represents the wave source, and G represents the receiver. According to... Figure 11It can be seen that the positive directions of the R component and T component of seismic data acquired from receivers in the first and third quadrants are the same, and the positive directions of the R component and T component of seismic data acquired from receivers in the second and fourth quadrants are the same. If the wave source is an SV wave source excited along the X direction, and the receiver belongs to the first or third quadrant in the XY coordinate system, then the second relationship data is acquired. Based on the first angle and the second relationship data, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result. If the wave source is an SV wave source excited along the X direction, and the receiver belongs to the second or fourth quadrant in the XY coordinate system, then the third relationship data is acquired. Based on the first angle and the third relationship data, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

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

[0117]

[0118] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the first row and second column, (180°-θ) represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction.

[0119] 1007. If the wave source type is an SV wave source excited along the X direction, and the receiver point belongs to the second or fourth quadrant in the XY coordinate system, the computer device acquires the third relation data. Based on the first angle and the third relation data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

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

[0121]

[0122] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive R-component and the positive X-component of the seismic data acquired from the receiver point of the SV source excited along the X-direction in the second or fourth quadrant. In the first row and second column, θ represents the angle between the positive R-component and the positive Y-component of the seismic data acquired from the receiver point of the SV source excited along the X-direction in the second or fourth quadrant. In the second row and first column, (180°-θ) represents the angle between the positive T-component and the positive X-component of the seismic data acquired from the receiver point of the SV source excited along the X-direction in the second or fourth quadrant. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of the seismic data acquired from the receiver point of the SV source excited along the X-direction in the second or fourth quadrant.

[0123] 1008. If the wave source type is an SH wave source excited along the Y direction, and the receiver point belongs to the first or third quadrant of the XY coordinate system, the computer equipment acquires the fourth relation data. Based on the first angle and the fourth relation data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0124] The positive directions of the R-component and T-component of seismic data are related to the positional relationship between the wave source and receiver, as well as 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 acquired from receivers in different quadrants are as follows: Figure 12 As shown, according to Figure 12 It can be seen that the positive directions of the R component and T component of seismic data acquired from receivers in the first and third quadrants are the same, and the positive directions of the R component and T component of seismic data acquired from receivers in the second and fourth quadrants are the same. If the wave source is an SH wave source excited in the Y direction, and the receiver belongs to the first or third quadrant in the XY coordinate system, then the fourth relation data is acquired. Based on the first angle and the fourth relation data, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result. If the wave source is an SH wave source excited in the Y direction, and the receiver belongs to the second or fourth quadrant in the XY coordinate system, then the fifth relation data is acquired. Based on the first angle and the fifth relation data, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

[0125] according to Figure 12 It can be seen that the angles between the positive directions of the R and X components, the R and Y components, the T and X components, and the T and Y components of seismic data can all be represented by the first angle θ. For example, when the receiver point is located in the first or third quadrant, the angle between the positive directions of the R and X components is represented as (90°+θ), the angle between the R and Y components as θ, the angle between the T and X components as θ, and the angle between the T and Y components as (90°-θ). Therefore, the fourth relation data can be represented as:

[0126]

[0127] Where 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. In the first row and first column, (90°+θ) represents the angle between the positive R-component and the positive X-component of the seismic data acquired from the receiver point of the SH source excited along the Y-direction in the first or third quadrant. In the second row and second column, θ represents the angle between the positive R-component and the positive Y-component of the seismic data acquired from the receiver point of the SH source excited along the Y-direction in the first or third quadrant. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of the seismic data acquired from the receiver point of the SH source excited along the Y-direction in the first or third quadrant. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of the seismic data acquired from the receiver point of the SH source excited along the Y-direction in the first or third quadrant.

[0128] 1009. If the wave source type is an SH wave source excited along the Y direction, and the receiver point belongs to the second or fourth quadrant of the XY coordinate system, the computer equipment acquires the fifth relation data. Based on the first angle and the fifth relation data, the seismic data acquired by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

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

[0130]

[0131] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive directions of the R component and the X component of the seismic data acquired from the receiver point in the second or fourth quadrant of the SH source excited along the Y direction. In the second row and second column, θ represents the angle between the positive directions of the R component and the Y component of the seismic data acquired from the receiver point in the second or fourth quadrant of the SH source excited along the Y direction. In the second row and first column, (180°-θ) represents the angle between the positive directions of the T component and the X component of the seismic data acquired from the receiver point in the second or fourth quadrant of the SH source excited along the Y direction. In the second row and second column, (90°-θ) represents the angle between the positive directions of the T component and the Y component of the seismic data acquired from the receiver point in the second or fourth quadrant of the SH source excited along the Y direction.

[0132] As described in steps 1006 to 1009 above, if the wave source type is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then relational data is obtained based on the wave source type and the quadrant to which the receiver point belongs in the XY coordinate system. That is, the computer device obtains relational data based on the wave source type, including: if the wave source type is an SV wave source excited along the X direction or an SH wave source excited along the Y direction, then relational data is obtained based on the wave source type and the quadrant to which the receiver 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: Figure 13 As shown, the Y component of the seismic data corresponding to the P-wave source is as follows: Figure 14 As shown, after mapping the seismic data corresponding to the P-wave source 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: Figure 15 As shown, the T component of the seismic data corresponding to the P-wave source is as follows: Figure 16 As shown. According to Figures 13 to 16 It can be seen that after wavefield separation of the seismic data corresponding to the P-wave source, most of the energy is transferred to the R component, and a small amount of energy 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: Figure 17 As shown, the Y component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Figure 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: Figure 19 As shown, the T component of the seismic data corresponding to the SV wave source excited along the X direction is as follows: Figure 20 As shown. According to Figures 17 to 20As indicated by the box in the diagram, after wavefield 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, while a small amount of SH wave remains 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: Figure 21 As shown, the Y component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows: Figure 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: Figure 23 As shown, the T component of the seismic data corresponding to the SH wave source excited along the Y direction is as follows: Figure 24 As shown. According to Figures 21 to 24 As indicated by the box in the diagram, after wavefield 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 SH wave, while a small amount of SV wave remains on the T component.

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

[0137] When setting up an observation system for three-dimensional nine-component seismic data, human error can be introduced. For example, the X-component of the location of the receiver point may not strictly adhere to the direction of the survey line, thus resulting in an error in the first angle.

[0138] To achieve more accurate wavefield separation from seismic data, this application employs the maximum energy ratio angle scanning method to correct the wavefield separation results. Specifically, the maximum energy ratio angle scanning method involves correcting a first angle to obtain multiple reference angles, performing an angle scan on the first angle, and then determining the wavefield separation result with the highest 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 using the first angle as a reference and gradually adjusting the angle towards both ends. For example, if the first angle is 30° and the preset angle interval is 0.5°, using the first angle as a reference and gradually adjusting the angle towards both ends means adding 0.5° to the first angle 30° to obtain a reference angle of 30.5°, 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°, subtracting 0.5° from 29.5° to obtain a reference angle of 29°, and so on.

[0140] To obtain the wavefield separation result corresponding to each reference angle, replace the first angle with the reference angle and execute steps 1005 to 1009 again.

[0141] 1011. The computer device acquires the wave field separation quality parameters of the wave field separation result corresponding to the first angle and the wave field separation quality parameters of the wave field separation result corresponding to each reference angle. Based on the acquired wave field separation quality parameters, the wave field separation result with the highest wave field separation quality is obtained.

[0142] Wavefield separation quality parameters are parameters used to represent the quality of wavefield separation results. 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 seismic data is expressed as: Where top T represents the start time when the receiver receives the seismic wave, bottom T represents the end time when the receiver receives the seismic wave, i represents the i-th time, R(i) represents the R component of the seismic data at the i-th time, and T(i) represents the T component of the seismic data at the i-th time.

[0144] The three-dimensional nine-component wavefield separation method provided in this application only requires determining the angle between the line connecting the wave source and the receiver in the XY coordinate system and the X-axis. Then, through relational data, the angles between the positive directions of the R-component and X-component, R-component and Y-component, T-component and X-component, and T-component and Y-component of the seismic data acquired by the receiver from the wave source can be obtained. Based on these angles, the seismic data acquired by the receiver is mapped from the XY coordinate system to the RT coordinate system. This wavefield separation method is not only simple but also applicable to seismic data from various wave sources, achieving wavefield separation of three-dimensional nine-component seismic data.

[0145] In addition, the embodiments of this 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] Figure 25 This is a schematic diagram of the structure of a three-dimensional nine-component wavefield separation device provided in an embodiment of this application, as shown below. Figure 25 As shown, the device includes:

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

[0148] The determination module 2502 is used to determine the direction of the survey line based on the coordinates of any two receiver points on the unchanged observation line;

[0149] Module 2503 is established to establish an XY coordinate system, where 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.

[0150] The separation module 2504 is used to determine a first angle for any receiver point of any wave source, based on the positions of the wave source and receiver point in the XY coordinate system. The first angle is the angle between the line connecting the wave source and receiver point in the XY coordinate system and the X-axis. Based on the wave source type, relational data is obtained. The relational data is used to indicate the relationship between the first angle and the angles between the positive directions of the R component and the positive directions of the X component, the R component and the positive directions of the Y component, the T component and the positive directions of the X component, and the T component and the positive directions of the Y component of the seismic data obtained from the receiver point. Based on the first angle and the relational data, the seismic data obtained from the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

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

[0152] The first relational data is represented as follows:

[0153]

[0154] Where 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, and θ is the first angle. In the first row and first column, θ represents the angle between the positive directions of the R component and the X component of the seismic data acquired from the receiver of the P-wave source. In the first row and second column, (90°-θ) represents the angle between the positive directions of the R component and the Y component of the seismic data acquired from the receiver of the P-wave source. In the second row and first column, (90°+θ) represents the angle between the positive directions of the T component and the X component of the seismic data acquired from the receiver of the P-wave source. In the second row and second column, θ represents the angle between the positive directions of the T component and the Y component of the seismic data acquired from the receiver of the P-wave source. cos represents the cosine function.

[0155] In one possible implementation, the separation module 2504 is used to obtain relational data based on the source type of the wave source and the quadrant to which the receiver point belongs in the XY coordinate system if the wave source type 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 second relational data if the wave source type is an SV wave source excited along the X direction and the receiver point belongs to the first or third quadrant in the XY coordinate system; and to obtain third relational data if the wave source type is an SV wave source excited along the X direction and the receiver point belongs to the second or fourth quadrant in the XY coordinate system.

[0157] In one possible implementation, the second relational data is represented as:

[0158]

[0159] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the first row and second column, (180°-θ) represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SV source excited along the X-direction.

[0160] Third-party relation data is represented as follows:

[0161]

[0162] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the X-direction. In the first row and second column, θ represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the X-direction. In the second row and first column, (180°-θ) represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the X-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the second or fourth quadrant of an SV source excited along the X-direction.

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

[0164] In one possible implementation, the fourth relation data is represented as:

[0165]

[0166] Where 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. In the first row and first column, (90°+θ) represents the angle between the positive R-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the first row and second column, θ represents the angle between the positive R-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the second row and first column, θ represents the angle between the positive T-component and the positive X-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction. In the second row and second column, (90°-θ) represents the angle between the positive T-component and the positive Y-component of seismic data acquired from receivers in the first or third quadrant of an SH source excited along the Y-direction.

[0167] The fifth relation data is represented as:

[0168]

[0169] Where 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. In the first row and first column, (90°-θ) represents the angle between the positive directions of the R and X components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the first row and second column, θ represents the angle between the positive directions of the R and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the second row and first column, (180°-θ) represents the angle between the positive directions of the T and X components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction. In the second row and second column, (90°-θ) represents the angle between the positive directions of the T and Y components of seismic data acquired from receivers in the second or fourth quadrant of an SH source excited along the Y direction.

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

[0171] The correction module is used to correct the first angle according to a preset angle interval to obtain multiple reference angles;

[0172] The second acquisition module is used to acquire the wavefield separation results corresponding to each reference angle;

[0173] The third acquisition module is used to acquire the wave field separation quality parameters of the wave field separation result corresponding to the first angle and the wave field separation quality parameters of the wave field separation result corresponding to each reference angle.

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

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

[0176] In some embodiments, the computer device is provided as a terminal. Figure 26 This is a structural block diagram of a terminal provided in an embodiment of this application. The terminal 2600 includes a processor 2601 and a memory 2602.

[0177] Processor 2601 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. Processor 2601 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). Processor 2601 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, processor 2601 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the screen. In some embodiments, processor 2601 may also include an AI (Artificial Intelligence) processor, which is used to handle computational 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 high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory devices. In some embodiments, the non-transitory computer-readable storage media in the memory 2602 are used to store at least one program code, which is executed by the processor 2601 to implement the three-dimensional nine-component wavefield separation method provided in the method embodiments of this application.

[0179] In some embodiments, the terminal 2600 may optionally include a peripheral device interface 2603 and at least one peripheral device. The processor 2601, memory 2602, and peripheral device interface 2603 can be connected via a bus or signal line. Each peripheral device can be connected to the peripheral device interface 2603 via a bus, signal line, or circuit board. Specifically, the peripheral device includes at least one of the following: 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] Peripheral device interface 2603 can be used to connect at least one I / O (Input / Output) related peripheral device to processor 2601 and memory 2602. In some embodiments, processor 2601, memory 2602 and peripheral device interface 2603 are integrated on the same chip or circuit board; in some other embodiments, any one or two of processor 2601, memory 2602 and peripheral device interface 2603 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0181] Display screen 2605 is used to display a UI (User Interface). This UI may include graphics, text, icons, videos, and any combination thereof. When display screen 2605 is a touch display screen, it also has the ability to collect touch signals on or above its surface. These touch signals can be input as control signals to processor 2601 for processing. In this case, display screen 2605 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there may be one display screen 2605, which serves as the front panel of terminal 2600; in other embodiments, there may be at least two display screens 2605, respectively disposed on different surfaces of terminal 2600 or in a folded design; in still other embodiments, display screen 2605 may be a flexible display screen, disposed on a curved or folded surface of terminal 2600. Furthermore, display screen 2605 may also be configured as a non-rectangular, irregular shape, i.e., a non-rectangular screen. The display screen 2605 can be made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

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

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

[0184] In some embodiments, the computer device is provided as a server. Figure 27This is a schematic diagram of a server structure provided in an embodiment of this application. The server 2700 can vary considerably due to different configurations or performance. It may include one or more Central Processing Units (CPUs) 2701 and one or more memories 2702. The memory 2702 stores at least one line of program code, which is loaded and executed by the processor 2701 to implement the methods provided in the above-described method embodiments. Of course, the server may also have wired or wireless network interfaces, a keyboard, and input / output interfaces for input and output. The server may also include other components for implementing device functions, which will not be elaborated here.

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

[0186] This application also provides a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to implement the three-dimensional nine-component wavefield separation method as described in any of the above implementations.

[0187] This application also provides a computer program product, which includes at least one piece of program code that is loaded and executed by a processor to implement the three-dimensional nine-component wavefield separation method as described in any of the above implementations.

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

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

Claims

1. A three-dimensional nine-component wavefield separation method, characterized in that, The method includes: Acquire three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X direction, and seismic data corresponding to SH-wave sources excited along the Y direction. The direction of the survey line is determined based on the coordinates of any two detector points on the unchanging observation line. The unchanging observation line refers to a survey line with uniformly distributed detector points. Establish an XY coordinate system, where 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. For any detector point of any wave source, based on the positions of the wave source and the detector point in the XY coordinate system, a first angle is determined. The first angle is the angle between the line connecting the wave source and the detector point in the XY coordinate system and the X-axis. If the wave source is the P wave source, then the first relationship data is obtained, and the first relationship data is represented as follows: Where 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, and T represents the T component of the seismic data. From the first angle, The first line and These represent the angles between the positive directions of the R component and the positive directions of the X and Y components of the seismic data acquired by the receiver point of the P-wave source, respectively. The second row... and denoted by , respectively, represent the angles between the positive directions of the T component and the positive directions of the X and Y components of the seismic data acquired by the receiver point of the P-wave source, and cos represents the cosine function; If the wave source type is the SV wave source excited along the X direction, and the detector point belongs to the first or third quadrant in the XY coordinate system, obtain the second relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the first or third quadrant of the SV wave source excited along the X-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X-direction. If the wave source type is the SV wave source excited along the X direction, and the detector point belongs to the second or fourth quadrant in the XY coordinate system, obtain the third relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the second or fourth quadrant of the SV wave source excited along the X-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X-direction. If the wave source type is the SH wave source excited along the Y direction, and the detector point belongs to the first or third quadrant in the XY coordinate system, obtain the fourth relational data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the first or third quadrant of the SH wave source excited along the Y-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y-direction. If the wave source type is the SH wave source excited along the Y direction, and the detector point belongs to the second or fourth quadrant in the XY coordinate system, obtain the fifth relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component directions and the positive X-component directions, and the positive Y-component directions, respectively, of the seismic data acquired from receivers belonging to the second or fourth quadrant of the SH wave source excited along the Y-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y-direction. Based on the first angle and the corresponding relationship data, the seismic data obtained by the receiver point is mapped from the XY coordinate system to the RT coordinate system to obtain the wavefield separation result.

2. The method according to claim 1, characterized in that, The method further includes: The first angle is corrected according to a preset angle interval to obtain multiple reference angles; Obtain the wavefield separation results corresponding to each reference angle; Obtain the wavefield separation quality parameters of the wavefield separation result corresponding to the first angle and the wavefield separation quality parameters of the wavefield separation result corresponding to each reference angle; Based on the obtained wavefield separation quality parameters, the wavefield separation result with the highest wavefield separation quality is obtained.

3. A three-dimensional nine-component wavefield separation device, characterized in that, The device includes: The first acquisition module is used to acquire three-dimensional nine-component seismic data, which includes seismic data corresponding to P-wave sources, seismic data corresponding to SV-wave sources excited along the X direction, and seismic data corresponding to SH-wave sources excited along the Y direction. The determination module is used to determine the direction of the survey line based on the coordinates of any two detector points on the unchanging observation line, wherein the unchanging observation line refers to a survey line with uniformly distributed detector points. A module is established 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 the direction orthogonal to the direction of the survey line and satisfies the right-hand rule. The separation module is used to determine a first angle for any detector point of any wave source, based on the positions of the wave source and the detector point in the XY coordinate system. The first angle is the angle between the line connecting the wave source and the detector point in the XY coordinate system and the X-axis. The separation module is further configured to, if the wave source is the P-wave source, obtain first relational data, wherein the first relational data is represented as: Where 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, and T represents the T component of the seismic data. From the first angle, The first line and These represent the angles between the positive directions of the R component and the positive directions of the X and Y components of the seismic data acquired by the receiver point of the P-wave source, respectively. The second row... and denoted by , respectively, represent the angles between the positive directions of the T component and the positive directions of the X and Y components of the seismic data acquired by the receiver point of the P-wave source, and cos represents the cosine function; If the wave source type is the SV wave source excited along the X direction, and the detector point belongs to the first or third quadrant in the XY coordinate system, obtain the second relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the first or third quadrant of the SV wave source excited along the X-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the first quadrant or the third quadrant of the SV wave source excited along the X-direction. If the wave source type is the SV wave source excited along the X direction, and the detector point belongs to the second or fourth quadrant in the XY coordinate system, obtain the third relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the second or fourth quadrant of the SV wave source excited along the X-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the second quadrant or the fourth quadrant of the SV wave source excited along the X-direction. If the wave source type is the SH wave source excited along the Y direction, and the detector point belongs to the first or third quadrant in the XY coordinate system, obtain the fourth relational data, which is expressed as: in, The first line and These represent the angles between the positive R-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired from receivers belonging to the first or third quadrant of the SH wave source excited along the Y-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the first quadrant or the third quadrant of the SH wave source excited along the Y-direction. If the wave source type is the SH wave source excited along the Y direction, and the detector point belongs to the second or fourth quadrant in the XY coordinate system, obtain the fifth relationship data, which is expressed as: in, The first line and These represent the angles between the positive R-component directions and the positive X-component directions, and the positive Y-component directions, respectively, of the seismic data acquired from receivers belonging to the second or fourth quadrant of the SH wave source excited along the Y-direction. The second row... and These represent the angles between the positive T-component direction and the positive X-component direction, and the positive Y-component direction, respectively, of the seismic data acquired by the receiver points belonging to the second quadrant or the fourth quadrant of the SH wave source excited along the Y-direction. The separation module is further configured to map the seismic data acquired by the receiver point from the XY coordinate system to the RT coordinate system based on the first angle and the corresponding relationship data, thereby obtaining the wavefield separation result.

4. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one piece of program code, the at least one piece of program code being loaded and executed by the processor to implement the three-dimensional nine-component wavefield separation method as described in any one of claims 1 to 2.

5. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to implement the three-dimensional nine-component wavefield separation method as described in any one of claims 1 to 2.

Images