Multi-node parallel output angle gather extraction method and device based on adaptive decomposition, electronic equipment and storage medium

Abstract

The invention discloses a multi-node parallel output angle gather extraction method and device based on adaptive decomposition, electronic equipment and a storage medium. The method comprises the following steps: calculating the number of nodes for angle gather output; decomposing along the Y direction based on the node number, and calculating the output range of each node; calculating a shot point end Poynting vector and a detection point end Poynting vector in the forward transmission process and the backward transmission process of the wave equation by adopting a finite difference algorithm, and further calculating an underground local reflection angle and a local azimuth angle; and respectively applying an angle domain imaging condition based on the output range of each node, storing corresponding reverse time migration according to the local reflection angle and the local azimuth angle to obtain an imaging value, obtaining a local angle domain common imaging point gather of each node, and reducing the local angle domain common imaging point gather as an overall angle domain common imaging point gather. According to the method, the computer memory requirement required by the angle gather can be effectively reduced, the calculation efficiency is improved, and the calculation cost is reduced.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to a method for extracting multi-node parallel output angle gathers based on adaptive decomposition, an electronic device, a storage medium, and an apparatus. Background Technology

[0002] Omnidirectional angular domain imaging can effectively eliminate artifacts produced by imaging in complex media and significantly improve image quality. In addition to subsurface imaging information, omnidirectional angular gathers also contain important parameters such as azimuth, dip, and subduction angle, providing fundamental data support for subsequent gather processing and applications. Angle gathers are typically considered as at least a five-dimensional data volume, meaning they require a very large amount of memory and place extremely high demands on computer memory. While angle gathers contain a wealth of seismic information, the high memory requirements for extracting the five-dimensional data volume significantly limit the application of this technology in production. Especially when processing actual three-dimensional seismic data, the computational demands often exceed the capacity of conventional computing resources. This not only increases the computational cost of the technology but also further restricts its promotion and application in the industry.

[0003] Therefore, current wave equation-based angle gather extraction techniques are usually limited to processing two-dimensional seismic data or small three-dimensional seismic data volumes, which greatly limits the widespread application of this technique in actual production.

[0004] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0005] The purpose of this invention is to propose a multi-node parallel output angle gather extraction method, electronic device, storage medium and apparatus based on adaptive decomposition, which solves the problems of large computer memory consumption and high computational cost in omnidirectional angle gather calculation, effectively reduces the computer memory requirement for angle gathers, improves computational efficiency and reduces computational cost.

[0006] To achieve the above objectives, this invention proposes a multi-node parallel output angle gather extraction method, electronic device, storage medium, and apparatus based on adaptive decomposition.

[0007] According to a first aspect of the present invention, a multi-node parallel output angle gather extraction method based on adaptive decomposition is proposed, comprising:

[0008] The number of nodes used for the output of the angle gather is calculated based on the angle gather dimension.

[0009] Based on the number of nodes, decompose along the Y direction and calculate the output range of each node;

[0010] The Poynting vector at the shot point during forward propagation of the wave equation is calculated using the finite difference algorithm based on the migration velocity field and seismic wavelet.

[0011] Based on the offset velocity field and shot data, the finite difference algorithm is used to calculate the Poynting vector at the receiver end during the back propagation of the wave equation.

[0012] Calculate the local underground reflection angle and local azimuth angle based on the Poynting vector at the shot point and the Poynting vector at the receiver point.

[0013] Based on the output range of each node, angle domain imaging conditions are applied respectively. The imaging value is obtained by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle. The local angle domain common imaging point gather corresponding to each node is obtained.

[0014] All the local angle domain common imaging point gathers are normalized to the global angle domain common imaging point gathers.

[0015] Optionally, the expression for calculating the number of nodes is:

[0016] N node =(nx*ny*nz*na*nzai*4) / M max ;

[0017] Where nzai is the azimuth angle, na is the subduction angle or inclination angle, and M max The maximum usable memory for a node is denoted by nx, ny, and nz, which represent the number of points in the three-dimensional space along the X, Y, and Z directions, respectively.

[0018] Optionally, the expression for calculating the output range of each node is:

[0019]

[0020] y i =(ny / N) node )*i,i<ny-1

[0021] y ny-1 =ny-(ny / N) node )*i;

[0022] y ny =ny-1

[0023] in, y represents the output range of the i-th node. i y i+1Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

[0024] Optionally, the expression for calculating the Poynting vector at the shot point is:

[0025]

[0026] in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x,t) represents the forward propagating wave field of the wave equation. t represents time, x represents three-dimensional space, and f(t) represents the seismic wavelet. This represents the Laplace operator.

[0027] Optionally, the calculation expression for the Poynting vector at the detector point is:

[0028]

[0029] in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0030] Optionally, the calculation expressions for the local reflection angle and the local azimuth angle are as follows:

[0031] Local reflection angle:

[0032] Local azimuth angle:

[0033] Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0034] Optionally, the local angle domain common imaging point gather is:

[0035]

[0036] in, Let be the common imaging point gather of the angular domain of the i-th node.

[0037] According to a second aspect of the present invention, a multi-node parallel output angle gather extraction device based on adaptive decomposition is proposed, comprising:

[0038] The first calculation module is used to calculate the number of nodes for the angle gather output based on the angle gather dimension;

[0039] The second calculation module is used to decompose the number of nodes along the Y direction and calculate the output range of each node.

[0040] The third calculation module is used to calculate the Poynting vector at the shot point during the forward propagation of the wave equation based on the migration velocity field and the seismic wavelet using the finite difference algorithm.

[0041] The fourth calculation module is used to calculate the Poynting vector at the receiver end during the back propagation of the wave equation using the finite difference algorithm based on the offset velocity field and the shot data.

[0042] The fifth calculation module is used to calculate the local underground reflection angle and local azimuth angle based on the Poynting vector at the shot point and the Poynting vector at the receiver point.

[0043] The forming module is used to apply angle domain imaging conditions based on the output range of each node, and obtain imaging values ​​by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle, thereby forming a local angle domain common imaging point gather for each node.

[0044] The reduction module is used to reduce all the local angle domain common imaging point gathers to the global angle domain common imaging point gathers.

[0045] According to a third aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0046] At least one processor; and,

[0047] A memory communicatively connected to the at least one processor; wherein,

[0048] The memory stores instructions executable by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the multi-node parallel output angle gather extraction method based on adaptive decomposition as described in any of the first aspects.

[0049] According to a fourth aspect of the present invention, a non-transitory computer-readable storage medium is provided, the non-transitory computer-readable storage medium storing computer instructions for causing a computer to perform the multi-node parallel output angle gather extraction method based on adaptive decomposition as described in any of the first aspects.

[0050] The beneficial effects of this invention are as follows: This invention calculates the number of nodes required for angle gather output by using the angle gather size and the available memory size of the input nodes. Then, based on the number of output nodes, it calculates the range of the angle gather output by each node. Next, it uses a finite difference algorithm to calculate the Poynting vectors of the forward and reverse propagating wavefields of the wavefield equations, thereby calculating the local subsurface reflection angle and azimuth. The corresponding imaging values ​​are stored according to the angles, forming a high-precision local angle gather, achieving adaptive region decomposition of the angle gather and parallel processing of different nodes. The decomposed angle gathers are output in parallel at different nodes. After all data calculations are completed, the angle gathers output in parallel from different nodes are normalized to the overall angle gather, achieving overall extraction of the angle gather. This invention effectively reduces the computer memory requirements for angle gathers, improves computational efficiency, and reduces computational costs.

[0051] The system of the present invention has other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0052] The above and other objects, features and advantages of the present invention will become more apparent from the accompanying drawings, in which like reference numerals generally denote like parts.

[0053] Figure 1 A flowchart illustrating the steps of the multi-node parallel output angle gather extraction method based on adaptive decomposition according to the present invention is shown.

[0054] Figure 2 A flowchart illustrating the steps of a multi-node parallel output angle gather extraction method based on adaptive decomposition according to Embodiment 2 of the present invention is shown.

[0055] Figure 3 A schematic diagram of the reverse-time offset model of the slit hole model according to Embodiment 2 of the present invention is shown.

[0056] Figure 4 A schematic diagram showing the calculation range of each node output of the slit hole model according to Embodiment 2 of the present invention is shown.

[0057] Figure 5 A schematic diagram of the angle gather obtained by the method of this embodiment according to Embodiment 2 of the present invention is shown.

[0058] Figure 6 A schematic diagram of an angle gather obtained without using the method of this embodiment is shown in Embodiment 2 of the present invention.

[0059] Figure 7 The following is an embodiment of the present invention. Figure 5 A schematic diagram of the offset result obtained by image stacking of angle gathers.

[0060] Figure 8 The following is an embodiment of the present invention. Figure 6 A schematic diagram of the offset result obtained by image stacking of angle gathers.

[0061] Figure 9 The following is an embodiment of the present invention. Figure 6 and Figure 7 A schematic diagram of the difference between the offset results. Detailed Implementation

[0062] The invention will now be described in more detail with reference to the accompanying drawings. While preferred embodiments of the invention are shown in the drawings, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

[0063] like Figure 1 As shown, a multi-node parallel output angle gather extraction method based on adaptive decomposition according to the present invention includes:

[0064] The number of nodes used for the output of the angle gather is calculated based on the angle gather dimension.

[0065] Decompose the nodes along the Y direction based on the number of nodes and calculate the output range of each node;

[0066] The Poynting vector at the shot point during forward propagation of the wave equation is calculated using the finite difference algorithm based on the migration velocity field and seismic wavelet.

[0067] The Poynting vector at the receiver end during the back propagation of the wave equation is calculated using the finite difference algorithm based on the offset velocity field and shot data.

[0068] Calculate the local underground reflection angle and local azimuth based on the Poynting vector at the shot point and the Poynting vector at the receiver point;

[0069] Based on the output range of each node, the angle domain imaging conditions are applied respectively. The imaging value is obtained by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle. The local angle domain common imaging point gather corresponding to each node is obtained.

[0070] All local angle domain common imaging point gathers are normalized to global angle domain common imaging point gathers.

[0071] Specifically, this invention calculates the number of nodes required for angle gather output by using the size of the angle gather and the available memory size of the input nodes. Then, based on the number of output nodes, it calculates the range of the angle gather output by each node. Next, it uses a finite difference algorithm to calculate the Poynting vectors of the forward and reverse propagating wavefields of the wavefield equation, thereby calculating the local subsurface reflection angle and azimuth. The corresponding imaging values ​​are stored according to the angles, forming a high-precision local angle gather, achieving adaptive region decomposition and parallel processing of different nodes. The decomposed angle gathers are output in parallel at different nodes. After all data calculations are completed, the angle gathers output in parallel from different nodes are normalized to the overall angle gather, achieving overall angle gather extraction. This invention effectively reduces the computer memory requirements for angle gathers, improves computational efficiency, reduces computational costs, and effectively avoids the bottleneck problem caused by excessive memory requirements. It significantly enhances the potential for the application of wave equation-based angle gather extraction technology in 3D actual seismic survey areas, providing more stable and reliable technical support for subsequent angle gather applications.

[0072] In one example, the expression for calculating the number of nodes is:

[0073] N node =(nx*ny*nz*na*nzai*4) / M max ;

[0074] Where nzai is the azimuth angle, na is the subduction angle or inclination angle, and M max The maximum usable memory for a node is denoted by nx, ny, and nz, which represent the number of points in the three-dimensional space along the X, Y, and Z directions, respectively.

[0075] In one example, the expression for calculating the output range of each node is:

[0076]

[0077] y i =(ny / N) node )*i,i<ny-1

[0078] y ny-1 =ny-(ny / N) node )*i;

[0079] y ny =ny-1

[0080] in, y represents the output range of the i-th node. i y i+1 Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

[0081] In one example, the expression for calculating the Poynting vector at the shot point is:

[0082]

[0083] in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x,t) represents the forward propagating wave field of the wave equation. t represents time, x represents three-dimensional space, and f(t) represents the seismic wavelet. This represents the Laplace operator.

[0084] In one example, the expression for calculating the Poynting vector at the receiver point is:

[0085]

[0086] in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0087] In one example, the expressions for calculating the local reflection angle and the local azimuth angle are as follows:

[0088] Local reflection angle:

[0089] Local azimuth angle:

[0090] Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0091] In one example, the local angle domain common imaging point gather is:

[0092]

[0093] in, Let be the common imaging point gather of the angular domain of the i-th node.

[0094] The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but this is not intended to limit the invention. It should be noted that, unless otherwise specified, the embodiments and features described in the embodiments of the present invention can be combined with each other.

[0095] Example 1

[0096] This embodiment provides a multi-node parallel output angle gather extraction method based on adaptive decomposition, including:

[0097] The number of nodes used for the output of the angle gather is calculated based on the angle gather dimension.

[0098] Decompose the nodes along the Y direction based on the number of nodes and calculate the output range of each node;

[0099] The Poynting vector at the shot point during forward propagation of the wave equation is calculated using the finite difference algorithm based on the migration velocity field and seismic wavelet.

[0100] The Poynting vector at the receiver end during the back propagation of the wave equation is calculated using the finite difference algorithm based on the offset velocity field and shot data.

[0101] Calculate the local underground reflection angle and local azimuth based on the Poynting vector at the shot point and the Poynting vector at the receiver point;

[0102] Based on the output range of each node, the angle domain imaging conditions are applied respectively. The imaging value is obtained by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle. The local angle domain common imaging point gather corresponding to each node is obtained.

[0103] All local angle domain common imaging point gathers are normalized to global angle domain common imaging point gathers.

[0104] The expression for calculating the number of nodes is:

[0105] N node =(nx*ny*nz*na / nzai*4) / M max ;

[0106] Where nzai is the azimuth angle, na is the subduction angle or inclination angle, and M maxThe maximum usable memory for a node is denoted by nx, ny, and nz, which represent the number of points in the three-dimensional space along the X, Y, and Z directions, respectively.

[0107] The expression for calculating the output range of each node is:

[0108]

[0109] y i =(ny / N) node )*i,i<ny-1

[0110] y ny-1 =ny-(ny / N) node )*i;

[0111] y ny =ny-1

[0112] in, y represents the output range of the i-th node. i y i+1 Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

[0113] The expression for calculating the Poynting vector at the shot point is:

[0114]

[0115] in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x,t) represents the forward propagating wave field of the wave equation. t represents time, x represents three-dimensional space, and f(t) represents the seismic wavelet. This represents the Laplace operator.

[0116] The expression for calculating the Poynting vector at the receiver point is:

[0117]

[0118] in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x,t) represents the inverted wave field of the wave equation. dobs For artillery data.

[0119] The expressions for calculating the local reflection angle and the local azimuth angle are as follows:

[0120] Local reflection angle:

[0121] Local azimuth angle:

[0122] Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0123] The local angle domain common imaging point gather is:

[0124]

[0125] in, Let be the common imaging point gather of the angular domain of the i-th node.

[0126] Example 2

[0127] like Figure 2 As shown, this embodiment provides a multi-node parallel output angle gather extraction method based on adaptive decomposition, including:

[0128] First, based on the input angle gather size and the available node memory size, the number of nodes required for angle gather output is calculated. Then, based on the number of output nodes, the range of the angle gather output by each node is calculated. Next, the Poynting vectors of the forward and reverse propagating wavefields of the wavefield equations are calculated separately using the finite difference algorithm, and the local subsurface reflection angle and azimuth angle are calculated. Then, the corresponding imaging values ​​are stored according to the angles to form high-precision angle gathers. Finally, these angle gathers are output in parallel at different nodes according to the previous adaptive decomposition scheme. After all data calculations are completed, the angle gathers output in parallel from different nodes are reduced to the overall angle gather. The specific implementation process includes the following steps:

[0129] (1) Input the offset velocity field v0, such as Figure 3 As shown, and the gun data d obs The seismic wavelet f(t) and the angle gather dimensions, including the number of points in three-dimensional space nx, ny, nz, azimuth nzai, subtended angle or dip angle na, and the maximum available memory M for each node. max Calculate the number of nodes N used for angle gather output. node ;

[0130] Number of nodes required: Nnode =(nx*ny*nz*na*nzai*4) / M max .

[0131] (2) Based on the number of nodes used for angle gather output calculated in step (1), decompose along the Y direction and calculate the calculation range of each node output, such as... Figure 4 As shown;

[0132]

[0133] y i =(ny / N) node )*i,i<ny-1

[0134] y ny-1 =ny-(ny / N) node )*i;

[0135] y ny =ny-1

[0136] in, y represents the output range of the i-th node. i y i+1 Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

[0137] (3) Read the input offset velocity field v0 and the seismic wavelet f(t), and calculate the Poynting vector at the shot point during the forward propagation of the wave equation using the finite difference algorithm;

[0138] Forward propagating wave field according to the wave equation:

[0139] Poynting vector at the gun point:

[0140] in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x,t) represents the forward propagation wave field of the wave equation, t represents time, x represents three-dimensional space, and f(t) represents the seismic wavelet. This represents the Laplace operator.

[0141] (4) Read the input offset velocity field v0 and shot data d obsThe Poynting vector at the receiver point is calculated using the finite difference algorithm during the back propagation of the wave equation.

[0142] Wave field propagated back through the wave equation:

[0143] Poynting vector at the detector point:

[0144]

[0145] in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x,t) represents the inverted wave field of the wave equation, d obs For artillery data.

[0146] (5) Calculate the local subsurface reflection angle θ and azimuth using the Poynting vector at the shot point and the Poynting vector at the receiver point.

[0147] Local reflection angle:

[0148] Local azimuth angle:

[0149] Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0150] (6) Based on the output range of each node obtained in step (2), apply the angle domain imaging conditions respectively, and store the corresponding reverse time offset according to the angle obtained in step (5) to obtain the imaging value. Each node obtains a local angle domain common imaging point gather.

[0151] (7) Reduce the local corner gathers obtained from each output node to form the final global corner gather, such as... Figure 5 As shown.

[0152] Figure 5 The angle gathers obtained by the method of this embodiment are shown, and Figure 6 The angle gather is produced directly without regional decomposition, and the two are consistent. Figure 7 and Figure 8 They are respectively Figure 5 and Figure 6 An offset profile is produced by overlaying angle gathers with imaging. Figure 9The difference between the two offset profiles is 0, which clearly shows that the angle gather obtained by the method of this embodiment is completely consistent with the angle gather calculated normally, thus demonstrating the correctness of the method of this embodiment. Figure 4 As can be seen from the angle gather calculation table, this invention can transform a single-node angle gather with a large memory footprint into multiple nodes, each with a smaller memory footprint. This demonstrates that this invention can effectively reduce the memory usage of a single node, increase the application potential of angle gathers in actual production, and provide a more reliable imaging tool for seismic exploration.

[0153] Example 3

[0154] This embodiment provides a multi-node parallel output angle gather extraction device based on adaptive decomposition, including:

[0155] The first calculation module is used to calculate the number of nodes for the angle gather output based on the angle gather dimension;

[0156] The second calculation module is used to decompose along the Y direction based on the number of nodes and calculate the output range of each node.

[0157] The third calculation module is used to calculate the Poynting vector at the shot point during the forward propagation of the wave equation based on the migration velocity field and the seismic wavelet using the finite difference algorithm.

[0158] The fourth calculation module is used to calculate the Poynting vector at the receiver end during the back propagation of the wave equation using the finite difference algorithm based on the offset velocity field and shot data.

[0159] The fifth calculation module is used to calculate the local subsurface reflection angle and local azimuth angle based on the Poynting vector at the shot point and the Poynting vector at the receiver point.

[0160] The forming module is used to apply angle domain imaging conditions based on the output range of each node, and obtain imaging values ​​by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle, thus forming the local angle domain common imaging point gather for each node.

[0161] The reduction module is used to reduce all local angle domain common imaging point gathers to global angle domain common imaging point gathers.

[0162] The expression for calculating the number of nodes is:

[0163] N node =(nx*ny*nz*na*nzai*4) / M max ;

[0164] Where nzai is the azimuth angle, na is the subduction angle or inclination angle, and M maxThe maximum usable memory for a node is denoted by nx, ny, and nz, which represent the number of points in the three-dimensional space along the X, Y, and Z directions, respectively.

[0165] The expression for calculating the output range of each node is:

[0166]

[0167] y i =(ny / N) node )*i,i<ny-1

[0168] y ny-1 =ny-(ny / N) node )*i;

[0169] y ny =ny-1

[0170] in, y represents the output range of the i-th node. i y i+1 Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

[0171] The expression for calculating the Poynting vector at the shot point is:

[0172]

[0173] in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x,t) represents the forward propagating wave field of the wave equation. t represents time, x represents three-dimensional space, and f(t) represents the seismic wavelet. This represents the Laplace operator.

[0174] The expression for calculating the Poynting vector at the receiver point is:

[0175]

[0176] in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x,t) represents the inverted wave field of the wave equation. dobs For artillery data.

[0177] The expressions for calculating the local reflection angle and the local azimuth angle are as follows:

[0178] Local reflection angle:

[0179] Local azimuth angle:

[0180] Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0181] The local angle domain common imaging point gather is:

[0182]

[0183] in, Let be the common imaging point gather of the angular domain of the i-th node.

[0184] Example 4

[0185] This disclosure also provides an electronic device, which includes:

[0186] At least one processor; and,

[0187] A memory communicatively connected to the at least one processor; wherein,

[0188] The memory stores instructions that can be executed by the at least one processor, which enables the at least one processor to perform the multi-node parallel output angle gather extraction method based on adaptive decomposition in Embodiment 1.

[0189] An electronic device according to embodiments of the present disclosure includes a memory and a processor. The memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0190] The processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0191] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0192] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0193] Example 5

[0194] This disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the multi-node parallel output angle gather extraction method based on adaptive decomposition in Embodiment 1.

[0195] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0196] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0197] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A multi-node parallel output angle gather extraction method based on adaptive decomposition, characterized in that, include: The number of nodes used for the output of the angle gather is calculated based on the angle gather dimension. Based on the number of nodes, decompose along the Y direction and calculate the output range of each node; The POynting vector at the shot point during forward propagation of the wave equation is calculated using the finite difference algorithm based on the migration velocity field and seismic wavelet. Based on the offset velocity field and shot data, the finite difference algorithm is used to calculate the POynting vector at the receiver end during the back propagation of the wave equation. Calculate the local underground reflection angle and local azimuth angle based on the POynting vector at the shot point and the POynting vector at the receiver point. Based on the output range of each node, angle domain imaging conditions are applied respectively. The imaging value is obtained by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle. The local angle domain common imaging point gather corresponding to each node is obtained. All the local angle domain common imaging point gathers are normalized to the global angle domain common imaging point gathers.

2. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The expression for calculating the number of nodes is: N node =(nx*ny*nz*na*nzai*4) / M max ; Where nzai is the azimuth angle, na is the subduction angle or inclination angle, and M max The maximum usable memory for a node is denoted by nx, ny, and nz, which represent the number of points in the three-dimensional space along the X, Y, and Z directions, respectively.

3. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The expression for calculating the output range of each node is: in, Let yi and yi' represent the output range of the i-th node. +1 Let y represent the start and end points of the i-th node along the Y-axis. ny-1 y ny These represent the start and end points of the last node along the Y-axis, respectively.

4. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The expression for calculating the POynting vector at the shot point is: in, Let v0(x) represent the Poynting vector at the gun point, and v0(x) represent the offset velocity field. The spatial derivative of the propagating wave field is represented by... u represents the time derivative of the propagating wave field. f (x, t) represents the forward propagating wave field of the wave equation. t represents time, x represents three-dimensional space, and T(t) represents the seismic wavelet. This represents the Laplace operator.

5. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The expression for calculating the Poynting vector at the detector point is: in, This represents the Poynting vector at the receiver point. The spatial derivative of the backpropagating wave field is represented. The time derivative of the backpropagating wavefield, u b (x, t) represents the inverted wave field of the wave equation. d obs For artillery data.

6. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The calculation expressions for the local reflection angle and the local azimuth angle are as follows: Local reflection angle: Local azimuth angle: Where θ is the local reflection angle. This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

7. The multi-node parallel output angle gather extraction method based on adaptive decomposition according to claim 1, characterized in that, The local angle domain common imaging point gather is: in, Let be the common imaging point gather of the angular domain of the i-th node.

8. An electronic device, characterized in that, The electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the multi-node parallel output angle gather extraction method based on adaptive decomposition as described in any one of claims 1-7.

9. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions for causing a computer to perform the multi-node parallel output angle gather extraction method based on adaptive decomposition as described in any one of claims 1-7.

10. A multi-node parallel output angle gather extraction device based on adaptive decomposition, characterized in that, include: The first calculation module is used to calculate the number of nodes for the angle gather output based on the angle gather dimension; The second calculation module is used to decompose the number of nodes along the Y direction and calculate the output range of each node. The third calculation module is used to calculate the Poynting vector at the shot point during the forward propagation of the wave equation based on the migration velocity field and the seismic wavelet using the finite difference algorithm. The fourth calculation module is used to calculate the Poynting vector at the receiver end during the back propagation of the wave equation using the finite difference algorithm based on the offset velocity field and the shot data. The fifth calculation module is used to calculate the local underground reflection angle and local azimuth angle based on the Poynting vector at the shot point and the Poynting vector at the receiver point. The forming module is used to apply angle domain imaging conditions based on the output range of each node, and obtain imaging values ​​by storing the corresponding reverse time offset according to the local reflection angle and local azimuth angle, thereby forming a local angle domain common imaging point gather for each node. The reduction module is used to reduce all the local angle domain common imaging point gathers to the global angle domain common imaging point gathers.

Images