Angle gather extraction method and device based on construction direction constraint, electronic equipment and storage medium

Abstract

The invention discloses an angle gather extraction method and device based on construction direction constraint, electronic equipment and a storage medium. The method comprises the following steps: extracting an underground structure angle based on an earlier-stage migration imaging result or a speed model; calculating a shot point end Poynting vector based on an underground structure angle and a wave equation forward propagation wave field; calculating a normal vector of an underground reflection interface based on the underground structure angle; calculating an underground local reflection angle and a local azimuth angle based on the shot point end Poynting vector and the normal vector; an angle domain imaging condition is applied, corresponding reverse time migration is stored based on an underground local reflection angle and a local azimuth angle to obtain an imaging value, and a high-precision omni-directional angle domain common imaging point gather is formed. According to the method, the angle gather precision loss caused by inaccurate detection point wave field Poynting vector direction is avoided, and the precision of the omnibearing angle domain common imaging point gather is improved, so that the seismic imaging quality is remarkably improved.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to an angle gather extraction method, electronic device, storage medium and apparatus based on tectonic orientation constraints. Background Technology

[0002] Omnidirectional angular domain imaging is an effective method to eliminate imaging artifacts and improve image quality in complex media. Omnidirectional angular domain common imaging point gathers not only consider the influence of background velocity but also integrate important information provided by the azimuth angle reflection coefficient. It plays a crucial role in multiple fields, including migration velocity analysis, AVA (amplitude as a function of angle) / AVO (amplitude as a function of migration distance) analysis, anisotropy analysis, and high-precision imaging. High-quality imaging gathers can significantly reduce the risks of seismic exploration, thereby providing more reliable data support for geological exploration and resource development.

[0003] The key to extracting omnidirectional angle-domain common imaging point gathers lies in accurately calculating subsurface reflection angles. Currently, considering both computational efficiency and imaging quality, methods based on Poynting direction vectors for extracting local subsurface angles have been widely adopted, effectively enabling the calculation of omnidirectional angle gathers. Traditional Poynting direction vector extraction methods require calculating the Poynting vectors at the shot point and receiver point separately. Typically, the wavefield at the shot point is relatively simple, so the Poynting vector calculated at the shot point is more accurate; however, the wavefield at the receiver point is often exceptionally complex, so the calculation of the Poynting vector at the receiver point often suffers from instability, thus affecting the imaging quality of the entire angle gather.

[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 an angle gather extraction method, electronic device, storage medium, and apparatus based on structural orientation constraints. This invention aims to solve the problem of inaccurate angle calculation caused by the instability of the Poynting vector at the receiver point, especially in complex geological environments where the accuracy of omnidirectional angle gathers is severely affected, thus impacting the imaging quality of subsequent gathers. This invention avoids the calculation of the Poynting vector at the receiver point and improves the accuracy of omnidirectional angle domain common imaging point gathers.

[0006] To achieve the above objectives, this invention proposes an angle gather extraction method, electronic device, storage medium, and apparatus based on construction direction constraints.

[0007] According to a first aspect of the present invention, a method for extracting angle gathers based on construction direction constraints is proposed, comprising:

[0008] Subsurface structural angles are extracted based on previous migration imaging results or velocity models;

[0009] The Poynting vector at the shot point is calculated based on the underground structural angle and the forward propagation wave field of the wave equation.

[0010] Calculate the normal vector of the underground reflection interface based on the aforementioned underground structural angle;

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

[0012] By applying angular domain imaging conditions, imaging values ​​are obtained by storing the corresponding reverse time offsets based on the local underground reflection angles and local azimuth angles, forming a high-precision all-around angular domain common imaging point gather.

[0013] Optionally, the expression for the forward propagating wave field of the wave equation is:

[0014]

[0015] Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

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

[0017]

[0018] in, Represents the Poynting vector at the shot point. This represents the spatial derivative of the propagating wave field. This represents the time derivative of the propagating wave field.

[0019] Optionally, the expression for calculating the normal vector is:

[0020]

[0021] in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

[0022] Optionally, the expression for calculating the local reflection angle is:

[0023]

[0024] Where θ is the local reflection angle.

[0025] Optionally, the expression for calculating the local azimuth angle is:

[0026]

[0027] in, This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0028] Optionally, the omnidirectional angle domain common imaging point gather is:

[0029]

[0030] in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0031] According to a second aspect of the present invention, an angle gather extraction device based on construction direction constraints is proposed, comprising:

[0032] The extraction module is used to extract subsurface structural angles based on previous migration imaging results or velocity models;

[0033] The first calculation module is used to calculate the Poynting vector at the shot point based on the underground structural angle and the forward propagation wave field of the wave equation.

[0034] The second calculation module is used to calculate the normal vector of the underground reflection interface based on the underground structural angle;

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

[0036] The forming module is used to apply angular domain imaging conditions, and obtain imaging values ​​based on the corresponding reverse time offset stored in the underground local reflection angle and local azimuth angle, forming a high-precision all-round angular domain common imaging point gather.

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

[0038] At least one processor; and,

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

[0040] 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 angle gather extraction method based on construction direction constraints as described in any of the first aspects.

[0041] According to a fourth aspect of the invention, a non-transitory computer-readable storage medium is provided, which stores computer instructions for causing a computer to perform the angle gather extraction method based on construction direction constraints as described in any of the first aspects.

[0042] The beneficial effects of this invention are as follows: This invention calculates the Poynting vector at the shot point and the normal vector of the subsurface reflecting surface by extracting the directional information of the subsurface structure. Then, it uses these two vectors to calculate the local subsurface reflection angle and azimuth angle, achieving more accurate angle gather generation. This invention avoids the accuracy loss of angle gathers caused by inaccurate Poynting vector direction at the receiver wavefield, thus significantly improving the quality of seismic imaging. It provides strong support for the practical application of omnidirectional angle domain common imaging point gathers and promotes the development of seismic exploration technology. Through this invention, the efficiency and accuracy of geological structure analysis and resource exploration are expected to be significantly improved, thereby opening up broader application prospects in related fields. This invention can effectively improve gather accuracy and imaging quality without increasing additional computational load, thus significantly enhancing the practicality of omnidirectional angle gather technology and providing a fast and high-precision imaging tool for high-precision imaging of actual 3D seismic data.

[0043] 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

[0044] 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.

[0045] Figure 1 A flowchart illustrating the steps of an angle gather extraction method based on construction direction constraints according to the present invention is shown.

[0046] Figure 2A flowchart illustrating the steps of an angle gather extraction method based on construction direction constraints according to Embodiment 2 of the present invention is shown.

[0047] Figure 3 A schematic diagram of a velocity model according to Embodiment 2 of the present invention is shown.

[0048] Figure 4 A schematic diagram of the X-direction angle field according to Embodiment 2 of the present invention is shown.

[0049] Figure 5 A schematic diagram of the Y-direction angle field according to Embodiment 2 of the present invention is shown.

[0050] Figure 6 A partial schematic diagram of the angle gather results obtained by the angle gather extraction method based on construction direction constraints according to Embodiment 2 of the present invention is shown.

[0051] Figure 7 A partial schematic diagram of the angle gather result obtained by the conventional omnidirectional angle gather method based on the Poytning vector direction according to Embodiment 2 of the present invention is shown.

[0052] Figure 8 The second embodiment of the present invention is shown. Figure 5 A schematic diagram of the results of angle gather stacking imaging. Detailed Implementation

[0053] 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.

[0054] like Figure 1 As shown, an angle gather extraction method based on construction direction constraints according to the present invention includes:

[0055] Subsurface structural angles are extracted based on previous migration imaging results or velocity models;

[0056] Calculate the Poynting vector at the shot point based on the subsurface structural angle and the forward propagation wave field of the wave equation.

[0057] Calculate the normal vector of the underground reflection interface based on the underground structural angle;

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

[0059] By applying angular domain imaging conditions, imaging values ​​are obtained by storing the corresponding reverse time offsets based on the local underground reflection angles and local azimuth angles, forming a high-precision all-around angular domain common imaging point gather.

[0060] Specifically, omnidirectional angle gathering technology has promising applications in various aspects of seismic exploration. However, conventional omnidirectional angle gathering technology based on the Poynting vector direction cannot eliminate the problem of inaccurate calculation of the Poynting vector direction at the receiver point, thus making it difficult to obtain high-precision angle gathers. This invention extracts subsurface structural angles based on previous migration imaging results or velocity models, obtaining the X and Y dip angles of the structure. Then, the extracted X and Y dip angles of the structure, along with the velocity field and shot data, are input. During the forward propagation of the wave equation, the Poynting vector at the shot point is calculated. The normal vector of the subsurface reflecting interface is then calculated using the extracted X and Y dip angles of the structure. The local reflection angle and local azimuth angle are then calculated using the Poynting vector and the normal vector, respectively. Finally, the corresponding imaging values ​​are stored based on the local reflection angle and local azimuth angle, forming a high-precision omnidirectional angle domain common imaging point gather. This invention calculates the Poynting vector at the shot point and the normal vector of the subsurface reflecting surface by extracting directional information from subsurface structures. Then, it uses these two vectors to determine the local subsurface reflection angle and azimuth, achieving more accurate angle gather generation. This invention avoids the accuracy loss of angle gathers caused by inaccurate Poynting vector direction at the receiver wavefield, thus significantly improving the quality of seismic imaging. It provides strong support for the practical application of omnidirectional angle domain common imaging point gathers and promotes the development of seismic exploration technology. Through this invention, the efficiency and accuracy of geological structure analysis and resource exploration are expected to be significantly improved, opening up broader application prospects in related fields. This invention can effectively improve gather accuracy and imaging quality without increasing additional computation, thus significantly enhancing the practicality of omnidirectional angle gather technology and providing a fast and high-precision imaging tool for high-precision imaging of actual 3D seismic data.

[0061] In one example, the expression for the propagating wave field in the wave equation is:

[0062]

[0063] Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

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

[0065]

[0066] in, Represents the Poynting vector at the shot point. This represents the spatial derivative of the propagating wave field. This represents the time derivative of the propagating wave field.

[0067] In one example, the expression for calculating the normal vector is:

[0068]

[0069] in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

[0070] In one example, the expression for calculating the local reflection angle is:

[0071]

[0072] Where θ is the local reflection angle.

[0073] In one example, the expression for calculating the local azimuth angle is:

[0074]

[0075] in, This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0076] In one example, the omnidirectional angle domain common imaging point gather is:

[0077]

[0078] in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0079] 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.

[0080] Example 1

[0081] This embodiment provides a method for angle gather extraction based on construction direction constraints, including:

[0082] Subsurface structural angles are extracted based on previous migration imaging results or velocity models;

[0083] Calculate the Poynting vector at the shot point based on the subsurface structural angle and the forward propagation wave field of the wave equation.

[0084] Calculate the normal vector of the underground reflection interface based on the underground structural angle;

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

[0086] By applying angular domain imaging conditions, imaging values ​​are obtained by storing the corresponding reverse time offsets based on the local underground reflection angles and local azimuth angles, forming a high-precision all-around angular domain common imaging point gather.

[0087] The expression for the propagating wave field in the wave equation is:

[0088]

[0089] Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

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

[0091]

[0092] in, Represents the Poynting vector at the shot point. This represents the spatial derivative of the propagating wave field. This represents the time derivative of the propagating wave field.

[0093] The expression for calculating the normal vector is:

[0094]

[0095] in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

[0096] The expression for calculating the local reflection angle is:

[0097]

[0098] Where θ is the local reflection angle.

[0099] The expression for calculating the local azimuth angle is:

[0100]

[0101] in, This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0102] The omnidirectional angular domain common imaging point gather is:

[0103]

[0104] in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0105] Example 2

[0106] like Figure 2 As shown, this embodiment provides a method for angle gather extraction based on construction direction constraints, including:

[0107] (1) Read the speed model v0, such as Figure 3 As shown, the sfdip command in Madagascar is used to extract the angles of underground structures, obtaining the X-direction angle α and the Y-direction angle β, as shown below. Figure 4 and Figure 5 As shown;

[0108] (2) Input the X-direction angle α and Y-direction angle β obtained in step (1), as well as the offset velocity field v0, wavelet f(t) and shot data d. obs During the forward propagation of the wave equation, the Poynting vector at the shot point is calculated.

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

[0110]

[0111] Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

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

[0113]

[0114] in, Represents the Poynting vector at the shot point. This represents the spatial derivative of the propagating wave field. This represents the time derivative of the propagating wave field.

[0115] (3) Using the extracted structural direction X-direction angle α and Y-direction angle β, the normal vector of the underground reflecting interface is calculated.

[0116] Normal vector of the underground reflective interface:

[0117]

[0118] in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

[0119] (4) Using the two vectors obtained in steps (2) and (3), calculate the local underground reflection angle θ and the local azimuth angle.

[0120] The expression for calculating the local reflection angle is:

[0121]

[0122] Where θ is the local reflection angle.

[0123] In one example, the expression for calculating the local azimuth angle is:

[0124]

[0125] in, This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0126] (5) Applying the angle domain imaging conditions, the imaging value is obtained by storing the corresponding reverse time offset based on the angle obtained in step (4), and a high-precision all-around angle domain common imaging point gather is output, such as... Figure 6 As shown.

[0127] The omnidirectional angular domain common imaging point gather is:

[0128]

[0129] in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation.

[0130] Reverse propagation Tron:

[0131]

[0132] Where, d obs This is for artillery data.

[0133] Figure 6 This demonstrates the omnidirectional angle domain common imaging point gather obtained through the method of this embodiment. Figure 7 This result was obtained using a conventional omnidirectional angle gather technique based on the Poynting vector direction. Comparing these two images, we can clearly see a significant difference: Figure 6 The results show significantly less noise, clearer phase axes, and higher overall imaging accuracy. This indicates that the angle gather extraction method based on construction orientation constraints has a significant advantage in improving gather quality. Furthermore, Figure 8 Showing based on Figure 6 The results of the gather imaging in the image are shown. The imaging results are not only clear and discernible, but also further verify the accuracy and effectiveness of our method. These results demonstrate the application potential of this invention in omnidirectional angle gather technology, providing a more reliable imaging tool for seismic exploration.

[0134] Example 3

[0135] This embodiment proposes an angle gather extraction device based on construction direction constraints, including:

[0136] The extraction module is used to extract subsurface structural angles based on previous migration imaging results or velocity models;

[0137] The first calculation module is used to calculate the Poynting vector at the shot point based on the underground structural angle and the forward propagation wave field of the wave equation.

[0138] The second calculation module is used to calculate the normal vector of the underground reflection interface based on the underground structural angle;

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

[0140] The forming module is used to apply angular domain imaging conditions, store the corresponding reverse time offset based on the local underground reflection angle and local azimuth angle to obtain the imaging value, and form a high-precision all-round angular domain common imaging point gather.

[0141] The expression for the propagating wave field in the wave equation is:

[0142]

[0143] Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

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

[0145]

[0146] in, Represents the Poynting vector at the shot point. This represents the spatial derivative of the propagating wave field. This represents the time derivative of the propagating wave field.

[0147] The expression for calculating the normal vector is:

[0148]

[0149] in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

[0150] The expression for calculating the local reflection angle is:

[0151]

[0152] Where θ is the local reflection angle.

[0153] The expression for calculating the local azimuth angle is:

[0154]

[0155] in, This is a local azimuth angle. Represents the unit normal vector in the Z direction. Represents the unit normal vector in the X direction.

[0156] The omnidirectional angular domain common imaging point gather is:

[0157]

[0158] in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

[0159] Example 4

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

[0161] At least one processor; and,

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

[0163] The memory stores instructions that can be executed by the at least one processor, which enables the at least one processor to perform the angle gather extraction method based on construction direction constraints in Embodiment 1.

[0164] 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.

[0165] 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.

[0166] 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.

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

[0168] Example 5

[0169] This disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the angle gather extraction method based on construction direction constraints in Embodiment 1.

[0170] 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.

[0171] 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).

[0172] 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 method for extracting angle gathers based on construction direction constraints, characterized in that, include: Subsurface structural angles are extracted based on previous migration imaging results or velocity models; The Poynting vector at the shot point is calculated based on the underground structural angle and the forward propagation wave field of the wave equation. Calculate the normal vector of the underground reflection interface based on the aforementioned underground structural angle; Calculate the local underground reflection angle and local azimuth angle based on the Poynting vector and normal vector at the shot point; By applying angular domain imaging conditions, imaging values ​​are obtained by storing the corresponding reverse time offsets based on the local underground reflection angles and local azimuth angles, forming a high-precision all-around angular domain common imaging point gather.

2. The angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The expression for the propagating wave field of the wave equation is as follows: Among them, u f (x,t) represents the propagating wave field, v0(x) represents the velocity field, f(t) represents the wavelet, t represents time, and x represents three-dimensional space. For the Laplace operator.

3. The angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The expression for calculating the Poynting vector at the gun point is: in, Represents the Poynting vector at the shot point. The spatial derivative of the propagating wave field is represented by... This represents the time derivative of the propagating wave field.

4. The angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The expression for calculating the normal vector is: in, u represents the normal vector of the underground reflecting interface. nx u ny u nz Let α and β represent the three spatial components of the normal vector of the underground reflective interface, respectively, and let α and β be the X-direction angle and Y-direction angle of the underground structure, respectively.

5. The angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The expression for calculating the local reflection angle is: Where θ is the local reflection angle.

6. The angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The expression for calculating the local azimuth angle is: in, 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 angle gather extraction method based on construction direction constraints according to claim 1, characterized in that, The omnidirectional angle domain common imaging point gather is: in, Represents the common imaging point collection in the omnidirectional angular domain, u b (x,t) represents the inverted wave field of the wave equation. d obs This is for artillery data.

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 angle gather extraction method based on construction direction constraints 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 angle gather extraction method based on construction direction constraints as described in any one of claims 1-7.

10. An angle gather extraction device based on structural orientation constraints, characterized in that, include: The extraction module is used to extract subsurface structural angles based on previous migration imaging results or velocity models; The first calculation module is used to calculate the Poynting vector at the shot point based on the underground structural angle and the forward propagation wave field of the wave equation. The second calculation module is used to calculate the normal vector of the underground reflection interface based on the underground structural angle; The third calculation module is used to calculate the local underground reflection angle and local azimuth angle based on the Poynting vector and normal vector at the shot point. The forming module is used to apply angular domain imaging conditions, and obtain imaging values ​​based on the corresponding reverse time offset stored in the underground local reflection angle and local azimuth angle, forming a high-precision all-round angular domain common imaging point gather.

Images