Stable and efficient seismic wave propagation vector determination method and device

Abstract

The invention discloses a stable and efficient seismic wave propagation vector determination method and device. The method comprises the following steps: solving a sound wave fluctuation equation according to a finite difference method to obtain seismic wave field data; calculating a spatial domain gradient of a seismic wave field at a certain moment, a time domain integral of the seismic wave field at a certain moment and a spatial domain gradient of the time domain integral of the seismic wave field at a certain moment based on the seismic wave field data; and finally determining a seismic wave propagation vector according to the seismic wave field data, the spatial domain gradient of the seismic wave field at a certain moment, the time domain integral of the seismic wave field at a certain moment and the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment. Compared with the traditional Poynting vector, the method has the advantages that the oscillation phenomenon of the traditional Poynting vector is eliminated, and the calculation complexity is relatively low.

Description

Technical Field

[0001] The present invention belongs to the technical field of seismic exploration, and in particular relates to a stable and efficient method and device for determining a seismic wave propagation vector. Background Art

[0002] The seismic wave propagation vector, also known as the Poynting vector, is an energy flux density vector. In the prior art, the Poynting vector represents the propagation direction of the energy flux density in the entire seismic wave field. The propagation direction obtained is also the propagation direction of the entire wave field, which can be used to calculate the earthquake propagation direction. Because its calculation efficiency is much higher than other methods, it is widely used in the field of seismic exploration.

[0003] The existing method for determining the Poynting vector contains a simple harmonic function term, which leads to oscillation in the determination of the Poynting vector. The determined wave field propagation direction is not accurate enough. In addition, the Poynting vector needs to be preprocessed before the next calculation, which consumes additional computing resources.

[0004] Therefore, how to reduce the computational complexity of the stable Poynting vector and avoid the oscillation phenomenon that exists when the traditional Poynting vector is determined is a technical problem to be solved by those skilled in the art. Summary of the Invention

[0005] The purpose of the present invention is to solve the technical problems in the prior art that the calculation of a stable Poynting vector is relatively complicated, or that an oscillation phenomenon occurs when the Poynting vector is determined.

[0006] To achieve the above technical objectives, the present invention provides, on the one hand, a stable and efficient method for determining a seismic wave propagation vector, the method comprising: The acoustic wave equation is solved using the finite difference method to obtain seismic wave field data; Calculating, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The seismic wave propagation vector is determined according to the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

[0007] Furthermore, the seismic wavefield data is specifically determined by the following formula: ; Where, is the seismic wave field data, is the plane wave amplitude, is the angular frequency term, is the wave propagation direction vector, is the propagation speed of the wave in the medium.

[0008] Furthermore, the spatial domain gradient of the seismic wave field at a certain moment is specifically determined by the following formula: ; Where, is the spatial gradient of the seismic wave field at a certain moment, is the partial differential symbol, is the spatial position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector.

[0009] Furthermore, the time domain integral of the seismic wave field at a certain moment is specifically determined by the following formula: ; Where, is the time domain integral of the seismic wave field at a certain moment, is the wave field data, is the integral term, is the time direction position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector, is the spatial position, is the constant term obtained from the integral calculation.

[0010] Furthermore, the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment is specifically determined by the following formula: ; Where, is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, is the time domain integral of the seismic wave field at a certain moment, is the spatial position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector.

[0011] Furthermore, the seismic wave propagation vector is determined based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment, specifically by the following formula: ; Where, is the seismic wave propagation vector, is the time domain integral of the seismic wave field at a certain moment, is the spatial gradient of the seismic wave field at a certain moment, is the seismic wave field data, is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, is the time direction position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector, For spatial location.

[0012] On the other hand, the present invention also provides a stable and efficient device for determining a seismic wave propagation vector, the device comprising: A solution module, used to solve the acoustic wave equation according to the finite difference method to obtain seismic wave field data; A first determination module is configured to calculate, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The second determination module is used to determine the seismic wave propagation vector based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

[0013] The present invention provides a stable and efficient method and device for determining seismic wave propagation vectors. Compared with existing technologies, this method first solves the acoustic wave equation using the finite difference method to obtain seismic wavefield data; then, based on the seismic wavefield data, calculates the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; finally, determines the seismic wave propagation vector based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment. Compared with the traditional Poynting vector, this method eliminates the oscillation phenomenon existing in the traditional Poynting vector and has lower computational complexity. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] In order to more clearly illustrate the embodiments of this specification or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0015] Figure 1 FIG2 is a flow chart of a stable and efficient method for determining a seismic wave propagation vector provided by an embodiment of this specification; Figure 2 FIG2 is a schematic diagram of the structure of a stable and efficient device for determining a seismic wave propagation vector provided in an embodiment of this specification; Figure 3 Shown is a schematic diagram of the seismic wave field; Figure 4 Shown is a schematic diagram of the horizontal component of the seismic wave propagation vector determined by this scheme; Figure 5 Shown is a schematic diagram of the vertical component of the seismic wave propagation vector determined by this scheme; Figure 6 Shown is a schematic diagram of the horizontal component of the Poynting vector determined by the conventional scheme; Figure 7 Shown is a schematic diagram of the vertical component of the Poynting vector determined by the conventional scheme; Figure 8 The figure shows the direction of the Poynting vector determined by the conventional scheme; Figure 9 Shown is a schematic diagram of the direction of the seismic wave propagation vector determined by the solution of this application. DETAILED DESCRIPTION

[0016] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without making creative efforts shall fall within the scope of protection of this application.

[0017] like Figure 1Shown is a flow chart of a stable and efficient method for determining the seismic wave propagation vector provided by an embodiment of this specification. Although this specification provides the method operation steps or device structures shown in the following embodiments or drawings, the method or device may include more or fewer operation steps or module units after partial merger based on routine or no creative labor. In the steps or structures where there is no necessary causal relationship logically, the execution order of these steps or the module structure of the device is not limited to the execution order or module structure shown in the embodiments or drawings of this specification. When the method or module structure is applied in actual devices, servers or terminal products, it can be executed sequentially or in parallel according to the method or module structure shown in the embodiments or drawings (for example, in a parallel processor or multi-threaded processing environment, or even in a distributed processing or server cluster implementation environment).

[0018] The stable and efficient method for determining the seismic wave propagation vector provided in the embodiments of this specification can be applied to various terminal devices, such as Figure 1 As shown, the method specifically includes the following steps: Step S101: Solve the acoustic wave equation according to the finite difference method to obtain seismic wave field data.

[0019] Specifically, seismic data is obtained, and then the acoustic wave equation is determined through forward simulation. The acoustic wave equation is as follows: ; Where, is the velocity field, is the seismic wave field, is the source function, , , are spatial horizontal, spatial vertical and temporal positions respectively, is the spatial direction vector, is the spatial coordinate of the earthquake source. The wave field data of the seismic wave field is obtained by solving the acoustic wave equation.

[0020] In the embodiment of the present application, the above-mentioned acoustic wave equation is solved by the finite difference method, and the seismic wave field data is specifically determined by the following formula: ; Where, is the seismic wave field data, is the plane wave amplitude, is the angular frequency term, is the wave propagation direction vector, is the propagation speed of the wave in the medium.

[0021] Step S102: Calculate the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment based on the seismic wavefield data.

[0022] Specifically, by solving the finite difference method, we can obtain the high-precision spatial gradient of the wave field, that is, the spatial gradient of the seismic wave field, where of The component calculation formula is: , of The component calculation formula is: , where is the difference order, is the finite difference coefficient. To balance the calculation accuracy and efficiency, this method uses the 10th order finite difference accuracy in the test. 、 The calculation data of the two items have been obtained in the forward simulation process, and reusing these two data can save the amount of calculation.

[0023] In the embodiment of the present application, the spatial domain gradient of the seismic wave field at a certain moment is specifically determined by the following formula: ; Where, is the spatial gradient of the seismic wave field, is the partial differential symbol.

[0024] Specifically, the time domain integral of the wave field at a certain moment It can be calculated from the wave field before that moment, The calculation formula is: Then the finite difference method is used to calculate The spatial gradient ,in of The component calculation formula is: , of The component calculation formula is: , the test also uses 10th order precision finite difference.

[0025] In the embodiment of the present application, the time domain integral of the seismic wave field at a certain moment is specifically determined by the following formula: ; Where, is the time domain integral of the seismic wave field at a certain moment, is the integral term, is the constant term obtained from the integral calculation.

[0026] The spatial domain gradient of the time domain integral of the seismic wave field at a certain moment is determined by the following formula: ; Where, is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, It is the time domain integral of the seismic wave field at a certain moment.

[0027] Step S103: determining a seismic wave propagation vector according to the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

[0028] In the embodiment of the present application, the seismic wave propagation vector is determined based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment, specifically by the following formula: ; Where, is the seismic wave propagation vector, is the time domain integral of the seismic wave field at a certain moment, is the spatial gradient of the seismic wave field, is the seismic wave field data, It is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment.

[0029] Specifically, the process of determining the seismic wave propagation vector in this application only needs to include a constant term at the end and the unit direction vector , does not involve trigonometric functions. Therefore, this method has two advantages in calculating the wave field propagation direction: first, since there are no trigonometric functions, this method does not have the oscillation phenomenon that occurs in the Poynting vector determined by conventional methods, which significantly improves the calculation stability of the seismic wave propagation vector; second, the method of this application only uses the unit direction vector Determine the size ratio between the vectors, and the coefficients of the calculated seismic wave propagation vectors are , the proportion of each component of the vector remains unchanged, which facilitates the direct calculation of the wave field propagation direction and improves the calculation accuracy of the seismic wave propagation vector. Figure 3 The diagram below shows the seismic wave field. Figure 4 The diagram shows the horizontal component of the seismic wave propagation vector determined by this scheme, as shown in Figure 5 The diagram shows the vertical component of the seismic wave propagation vector determined by this scheme, as shown in Figure 6 The diagram shows the horizontal component of the Poynting vector determined by the conventional method, as shown in Figure 7 The figure shows the schematic diagram of the vertical component of the Poynting vector determined by the conventional scheme. In the figure, the horizontal axis Distance / km represents the distance in kilometers, and the vertical axis Depth / km represents the depth in kilometers. Figures 4 to 7 It can be seen that after the conventional scheme determines the Poynting vector, there is a wavefront oscillation phenomenon in both the horizontal and vertical components. However, the seismic wave propagation vector determined by the present application scheme has no oscillation phenomenon in both the horizontal and vertical components. Figure 8 The following is a schematic diagram of the direction of the Poynting vector determined by the conventional scheme, as shown in Figure 9 The figure shows the direction of the seismic wave propagation vector determined by the present application scheme. Figures 3 to 9 In the image, the crosshair in the center represents the earthquake source, which is represented by Figure 8 and Figure 9 It can be seen that the direction of the seismic wave propagation vector determined by the present application solution is more accurate in the details than that determined by the conventional solution.

[0030] Based on the above-mentioned stable and efficient method for determining the propagation vector of seismic waves, one or more embodiments of this specification also provide a platform and terminal for determining the propagation vector of seismic waves in a stable and efficient manner. The platform or terminal may include devices, software, modules, plug-ins, servers, clients, etc. that use the methods described in the embodiments of this specification and are combined with necessary hardware implementation devices. Based on the same innovative concept, the systems in one or more embodiments provided in the embodiments of this specification are as described in the following embodiments. Since the implementation scheme and method for solving the problem of the system are similar, the implementation of the specific system in the embodiments of this specification can refer to the implementation of the aforementioned method, and the repetitions will not be repeated. The terms "unit" or "module" used below can be a combination of software and / or hardware that implements predetermined functions. Although the system described in the following embodiments is preferably implemented in software, hardware and a combination of software and hardware are also possible and conceived.

[0031] Specifically, Figure 2 This is a schematic diagram of the module structure of an embodiment of a stable and efficient seismic wave propagation vector determination device provided in this specification, such as Figure 2 As shown, the stable and efficient seismic wave propagation vector determination device provided in this specification includes: A solution module 201 is used to solve the acoustic wave equation according to the finite difference method to obtain seismic wave field data; A first determining module 202 is configured to calculate, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The second determination module 203 is used to determine the seismic wave propagation vector based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

[0032] It should be noted that the above-mentioned system may also include other implementation methods according to the description of the corresponding method embodiment. The specific implementation methods can refer to the description of the above-mentioned corresponding method embodiment, and will not be described one by one here.

[0033] An embodiment of the present application further provides an electronic device, including: processor; a memory for storing instructions executable by the processor; The processor is configured to execute the method provided in the above embodiment.

[0034] The electronic device provided in the embodiment of the present application stores executable instructions of the processor in a memory. When the processor executes the executable instructions, it can first solve the acoustic wave equation according to the finite difference method to obtain seismic wave field data; then calculate the spatial domain gradient of the seismic wave field at a certain moment, the time domain integral of the seismic wave field at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment based on the seismic wave field data; finally, determine the seismic wave propagation vector based on the seismic wave field data, the spatial domain gradient of the seismic wave field at a certain moment, the time domain integral of the seismic wave field at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment. Compared with the traditional Poynting vector, the oscillation phenomenon existing in the traditional Poynting vector is eliminated, and the computational complexity is lower.

[0035] The foregoing description of this specification describes specific embodiments. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in an order different from that described in the embodiments and still achieve the desired results. Furthermore, the processes depicted in the accompanying drawings do not necessarily require the specific order shown or the sequential order to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0036] The methods or devices described in the above embodiments of this specification can implement business logic through computer programs and record them on storage media. The storage media can be read and executed by a computer to achieve the effects of the solutions described in the embodiments of this specification, such as: The acoustic wave equation is solved using the finite difference method to obtain seismic wave field data; Calculating, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The seismic wave propagation vector is determined according to the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

[0037] The storage medium may include a physical device for storing information, typically digitizing the information and then storing it in a medium utilizing electrical, magnetic, or optical means. Examples of such storage media include: devices that store information electrically, such as various types of memory devices like RAM and ROM; devices that store information magnetically, such as hard disks, floppy disks, magnetic tapes, magnetic core memories, bubble memories, and USB flash drives; and devices that store information optically, such as CDs and DVDs. Of course, other types of readable storage media exist, such as quantum memories and graphene memories.

[0038] The embodiments of this specification are not limited to those that must comply with industry communication standards, standard computer resource data update and data storage rules, or the situations described in one or more embodiments of this specification. Certain industry standards or slightly modified implementation plans based on the implementation described in the embodiments using custom methods or embodiments can also achieve the same, equivalent, or similar implementation effects as the above embodiments, or the expected implementation effects after deformation. The embodiments obtained by applying these modified or deformed data acquisition, storage, judgment, processing methods, etc. can still fall within the scope of the optional implementation plans of the embodiments of this specification.

[0039] The controller can be implemented in any suitable manner. For example, the controller can take the form of a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro)processor, logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers. Examples of controllers include, but are not limited to, the following microcontrollers: ARC 625D, Atmel ATMEL AT91SAM, Microchip PIC18F26K20, and Silicon Labs C8051F320. The memory controller can also be implemented as part of the control logic of the memory. Those skilled in the art will also appreciate that, in addition to implementing the controller in pure computer-readable program code, the controller can also be implemented in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded microcontrollers by logically programming the method steps. Therefore, such a controller can be considered a hardware component, and the means for implementing various functions included therein can also be considered as structures within the hardware component. Alternatively, the means for implementing various functions can be considered as both software modules implementing the method and structures within the hardware component.

[0040] The device embodiments described above are merely illustrative. For example, the division of units described is merely a logical functional division. Actual implementations may employ alternative divisions, such as combining or integrating multiple units or plug-ins into another system, or omitting or disabling certain features. Furthermore, the coupling or direct coupling or communication connection shown or discussed between devices or units may be through interfaces, or indirect coupling or communication connection between devices or units may be electrical, mechanical, or otherwise.

[0041] These computer program instructions can also be loaded onto a computer or other programmable resource data updating device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for executing on the computer or other programmable device to implement the process. Figure 1 a process or multiple processes and / or boxes Figure 1 The steps for the function specified in one or more boxes.

[0042] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from the other embodiments. In particular, since the system embodiments are generally similar to the method embodiments, their description is relatively simple, and relevant parts can be referenced to the partial description of the method embodiments. Throughout this specification, reference to the terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of this specification. In this specification, the schematic representations of these terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and integrate the different embodiments or examples, and features of different embodiments or examples, described in this specification, without conflict.

[0043] Those skilled in the art will appreciate that the embodiments described herein are intended to help readers understand the principles of the present invention, and it should be understood that the scope of protection of the present invention is not limited to such specific descriptions and embodiments. Those skilled in the art can make various other specific variations and combinations based on the technical teachings disclosed in the present invention without departing from the essence of the present invention, and such variations and combinations are still within the scope of protection of the present invention.

Claims

1. A stable and efficient method for determining seismic wave propagation vector, characterized in that: The method comprises: The acoustic wave equation is solved using the finite difference method to obtain seismic wave field data; Calculating, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The seismic wave propagation vector is determined according to the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

2. The stable and efficient method for determining the seismic wave propagation vector according to claim 1, wherein: The seismic wave field data is specifically determined by the following formula: ； Where, is the seismic wave field data, is the plane wave amplitude, is the angular frequency term, is the wave propagation direction vector, is the propagation speed of the wave in the medium.

3. The stable and efficient method for determining the seismic wave propagation vector according to claim 1, wherein: The spatial domain gradient of the seismic wave field at a certain moment is specifically determined by the following formula: ； Where, is the spatial gradient of the seismic wave field at a certain moment, is the partial differential symbol, is the spatial position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector.

4. The stable and efficient method for determining the seismic wave propagation vector according to claim 1, wherein: The time domain integral of the seismic wave field at a certain moment is specifically determined by the following formula: ； Where, is the time domain integral of the seismic wave field at a certain moment, is the wave field data, is the integral term, is the time direction position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector, is the spatial position, is the constant term obtained from the integral calculation.

5. The stable and efficient method for determining seismic wave propagation vector according to claim 1, wherein: The spatial domain gradient of the time domain integral of the seismic wave field at a certain moment is specifically determined by the following formula: ； Where, is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, is the time domain integral of the seismic wave field at a certain moment, is the spatial position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector.

6. The stable and efficient method for determining seismic wave propagation vector according to claim 1, wherein: The seismic wave propagation vector is determined based on the seismic wave field data, the spatial domain gradient of the seismic wave field at a certain moment, the time domain integral of the seismic wave field at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, specifically by the following formula: ； Where, is the seismic wave propagation vector, is the time domain integral of the seismic wave field at a certain moment, is the spatial gradient of the seismic wave field at a certain moment, is the seismic wave field data, is the spatial domain gradient of the time domain integral of the seismic wave field at a certain moment, is the time direction position, is the plane wave amplitude, is the angular frequency term, is the propagation speed of the wave in the medium, is the time direction position, is the wave propagation direction vector, For spatial location.

7. A stable and efficient device for determining seismic wave propagation vector, characterized in that: The device comprises: A solution module, used to solve the acoustic wave equation according to the finite difference method to obtain seismic wave field data; A first determination module is configured to calculate, based on the seismic wavefield data, a spatial domain gradient of the seismic wavefield at a certain moment, a time domain integral of the seismic wavefield at a certain moment, and a spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment; The second determination module is used to determine the seismic wave propagation vector based on the seismic wavefield data, the spatial domain gradient of the seismic wavefield at a certain moment, the time domain integral of the seismic wavefield at a certain moment, and the spatial domain gradient of the time domain integral of the seismic wavefield at a certain moment.

Images