Optimized storage method, system and equipment for Gaussian beam travel time table and medium

Abstract

The invention relates to the technical field of storage, and discloses a Gaussian beam travel time table optimization storage method, system and device and a medium, and the method comprises the steps: obtaining a starting shot point in a Gaussian beam travel time table and the initial direction of a Gaussian beam; performing ray tracing based on the initial shot point and the initial direction of the Gaussian beam to obtain a propagation path of a central ray of the Gaussian beam; determining the effective width of the Gaussian beam based on the propagation path of the central ray of the Gaussian beam and the starting shot point; based on the effective width of the Gaussian beam, determining the minimum coverage rectangular grid range of the Gaussian beam; performing local extraction on the Gaussian beam travel time table based on the minimum coverage rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table; and storing the optimized Gaussian beam travel time table. According to the invention, the storage cost of the Gaussian beam travel time table can be greatly reduced.

Description

Technical Field

[0001] The present invention relates to the field of storage technology, and in particular to a Gaussian beam travel time table optimization storage method, system, equipment and medium. Background Art

[0002] Gaussian beam migration is a widely used migration imaging method in the industry. During its computational implementation, it is necessary to calculate and store Gaussian beam travel time tables for different propagation directions underground. The traditional Gaussian beam travel time table storage method stores the entire underground imaging grid space, which consumes a large amount of data storage space, seriously affecting the computational efficiency and scalability of Gaussian beam migration imaging technology.

[0003] The traditional Gaussian beam travel time table is defined based on the overall underground imaging grid space. The technical implementation process is as follows: first, the central path of the Gaussian beam is calculated by ray tracing; second, the effective width of the Gaussian beam is calculated based on the central ray, and the distribution range of the Gaussian beam travel time table is determined on this basis; next, the overall imaging grid containing the Gaussian beam is extracted and stored.

[0004] Since the effective range of the Gaussian beam is significantly smaller than the overall imaging grid, the traditional method will waste a lot of storage space, seriously affecting the computational efficiency and scalability of Gaussian beam migration imaging technology. Summary of the invention

[0005] The purpose of the present invention is to provide a Gaussian beam travel time table optimization storage method, system, device and medium to solve the problem that traditional methods waste a lot of storage space and seriously affect the calculation efficiency and scalability of Gaussian beam migration imaging technology.

[0006] In order to achieve the above-mentioned invention object, the technical solution adopted by the present invention is: In a first aspect, the present invention provides a method for optimizing storage of a Gaussian beam travel time table, the method comprising: Obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; Ray tracing is performed based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam; Based on the propagation path and the starting shot point of the central ray of the Gaussian beam, the effective width of the Gaussian beam is determined; Based on the effective width of the Gaussian beam, determine the minimum coverage rectangular grid range of the Gaussian beam; Based on the minimum coverage rectangular grid range of the Gaussian beam, the Gaussian beam travel time table is locally extracted to obtain the optimized Gaussian beam travel time table; The optimized Gaussian beam travel time table is stored.

[0007] Preferably, ray tracing is performed based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam, including: At the starting shot point, ray tracing is performed according to the initial direction of the Gaussian beam, and the kinematic ray tracing partial differential equation is obtained; The kinematic ray tracing partial differential equation is numerically solved based on the Runge-Kutta algorithm to obtain the propagation path of the central ray of the Gaussian beam.

[0008] Preferably, determining the effective width of the Gaussian beam based on the propagation path and the starting shot point of the central ray of the Gaussian beam comprises: Extract any point on the propagation path of the central ray of the Gaussian beam as an unknown point; Obtain the velocity of the unknown point, and determine the kinematic ray parameters of the unknown point based on the velocity of the unknown point; Obtain the seismic wave propagation travel time from the Gaussian beam travel time table; Determine the effective width of the unknown point based on the kinematic ray parameters of the unknown point, the travel time of seismic wave propagation, and the starting shot point; The effective width of the Gaussian beam is obtained by combining the effective widths corresponding to all points on the propagation path of the central ray of the Gaussian beam.

[0009] Preferably, the calculation expression of the kinematic ray parameters of the unknown point is: ; In the formula, represents the kinematic ray parameters of the unknown point P, represents the velocity of the unknown point P, and dt represents the time sampling interval of ray tracing.

[0010] Preferably, the calculation expression of the effective width of the unknown point is: ; In the formula, represents the effective width of the unknown point P, B represents the starting shot point, T represents the propagation time of the seismic wave, represents the kinematic ray parameters of the unknown point P, Represents the kinematic ray parameters of the starting shot point.

[0011] Preferably, based on the effective width of the Gaussian beam, determining the minimum coverage rectangular grid range of the Gaussian beam includes: Traverse all points on the propagation path of the central ray of the Gaussian beam, and extract the horizontal range and vertical range corresponding to each point according to the effective width corresponding to the point; Based on the horizontal range and vertical range corresponding to all points, the minimum coverage rectangular grid range of the Gaussian beam is determined.

[0012] Preferably, storing the optimized Gaussian beam travel time table includes: In the optimized Gaussian beam travel time table, a two-dimensional array is constructed based on the horizontal range and vertical range corresponding to each point; The optimized Gaussian beam travel time table is stored based on a two-dimensional array.

[0013] In a second aspect, the present invention provides a Gaussian beam travel time table optimization storage system, which is used to implement the above-mentioned Gaussian beam travel time table optimization storage method, and the system includes: A data acquisition module, used to obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; A path tracing module is used to perform ray tracing based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam; A width calculation module, used to determine the effective width of the Gaussian beam based on the propagation path and the starting shot point of the central ray of the Gaussian beam; A range calculation module, used for determining the minimum coverage rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam; A range extraction module is used to perform local extraction of the Gaussian beam travel time table based on the minimum coverage rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table; The travel time table storage module is used to store the optimized Gaussian beam travel time table.

[0014] In a third aspect, the present invention provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned Gaussian beam travel time table optimization storage method when executing the computer program.

[0015] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned Gaussian beam travel time table optimization storage method.

[0016] The beneficial effects of the present invention are mainly embodied in: The present invention can quickly determine the effective width of the Gaussian beam according to the propagation path of the central ray of the Gaussian beam and the starting shot point. Based on the effective width of the Gaussian beam, the minimum coverage rectangular grid range of the Gaussian beam can be dynamically determined. The Gaussian beam travel time table is extracted and stored based on the minimum coverage rectangular grid range of the Gaussian beam, which can greatly reduce the storage cost of the Gaussian beam travel time table. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following specific implementations, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings: Figure 1 It is a flow chart of a Gaussian beam travel time table optimization storage method provided by one embodiment of the present invention; Figure 2 is a schematic diagram of the travel time of a Gaussian beam in a Gaussian beam travel time table provided in one embodiment of the present invention; Figure 3 is a schematic diagram of the effective width of a Gaussian beam provided by one embodiment of the present invention; Figure 4 is a schematic diagram of the minimum coverage rectangular grid range of a Gaussian beam provided by an embodiment of the present invention; Figure 5 It is a comparative schematic diagram of the storage space of the optimized travel time table provided by one embodiment of the present invention; Figure 6 It is a schematic diagram of an optimized extraction result of a Gaussian beam travel time table with different initial directions provided by an embodiment of the present invention; Figure 7 is a schematic diagram of another optimized extraction result of a Gaussian beam travel time table with a different initial direction provided by an embodiment of the present invention; Figure 8 It is a block diagram of a Gaussian beam travel time table optimization storage system provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0018] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the present invention will be briefly introduced below in combination with the drawings and the description of the embodiments or the prior art. Obviously, the following description of the structure of the drawings is only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. It should be noted that the description of these embodiments is used to help understand the present invention, but does not constitute a limitation of the present invention.

[0019] Embodiment 1 Figure 1 FIG. 1 is a flow chart of a method for optimizing storage of Gaussian beam travel time tables provided in one embodiment of the present invention. Figure 1 As shown, this embodiment provides a method for optimizing storage of a Gaussian beam travel time table, the method comprising: Step S10: Obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; in this embodiment, Figure 2 As shown, Figure 2is the travel time of a Gaussian beam in the Gaussian beam travel time table. The starting shot point and the initial direction of the Gaussian beam can be extracted from the travel time.

[0020] Step S20: performing ray tracing based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam.

[0021] As a further optimization of this embodiment, ray tracing is performed based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam, including: Step S201: At the starting shot point, ray tracing is performed according to the initial direction of the Gaussian beam to obtain a kinematic ray tracing partial differential equation.

[0022] Step S202: numerically solve the kinematic ray tracing partial differential equation based on the Runge-Kutta algorithm to obtain the propagation path of the central ray of the Gaussian beam.

[0023] In this embodiment, the Runge-Kutta algorithm is a classic numerical method for solving initial value problems of ordinary differential equations (ODEs), and its core idea is to improve the calculation accuracy by weighted averaging of multi-stage slopes.

[0024] In this embodiment, the kinematic ray tracing partial differential equation is a conventional equation in the art, and its specific function expressions are not described one by one in this embodiment.

[0025] After solving step S201 to step S202, we can get Figure 2 The blue curve part of the propagation path.

[0026] Step S30: Determine the effective width of the Gaussian beam based on the propagation path of the central ray of the Gaussian beam and the starting shot point.

[0027] In this embodiment, the traditional Gaussian beam width calculation method requires extremely time-consuming dynamic ray tracing. To ensure calculation efficiency, the effective width of the Gaussian beam is calculated as follows: Step S301: extract any point on the propagation path of the central ray of the Gaussian beam as an unknown point.

[0028] Step S302: Obtain the velocity of the unknown point, and determine the kinematic ray parameters of the unknown point based on the velocity of the unknown point; wherein the calculation expression of the kinematic ray parameters of the unknown point is: ; In the formula, represents the kinematic ray parameters of the unknown point P, represents the velocity of the unknown point P, and dt represents the time sampling interval of ray tracing.

[0029] Step S303: Obtain the seismic wave propagation travel time of the Gaussian beam travel time table.

[0030] Step S304: Determine the effective width of the unknown point based on the kinematic ray parameters of the unknown point, the travel time of seismic wave propagation, and the starting shot point; wherein the calculation expression of the effective width of the unknown point is: ; In the formula, represents the effective width of the unknown point P, B represents the starting shot point, T represents the propagation time of the seismic wave, represents the kinematic ray parameters of the unknown point P, Represents the kinematic ray parameters of the starting shot point.

[0031] Step S305: combining the effective widths corresponding to all points on the propagation path of the central ray of the Gaussian beam to obtain the effective width of the Gaussian beam.

[0032] In this embodiment, the velocity of the unknown point P, the travel time of seismic wave propagation and other parameters can be obtained by kinematic ray tracing partial differential equations, without the need for time-consuming dynamic ray tracing, which can effectively improve the calculation efficiency of the Gaussian beam travel time table.

[0033] After the processing of steps S301 to S305, the effective width of the Gaussian beam is obtained as follows: Figure 3 As shown, Figure 3 The three red curves in the figure are used to represent the range of the effective width of the three points, and the arc shape is used to represent that the amplitude of the Gaussian beam gradually decays along the direction of the arc, showing a trend of being high in the middle and decaying on both sides.

[0034] Step S40: Determine the minimum coverage rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam.

[0035] As a further optimization of this embodiment, based on the effective width of the Gaussian beam, the minimum coverage rectangular grid range of the Gaussian beam is determined, including: Step S401: traverse all points on the propagation path of the central ray of the Gaussian beam, and extract the horizontal range and vertical range corresponding to each point according to the effective width corresponding to the point. The horizontal range and vertical range are as follows: Figure 4 As shown, in Figure 4 In the example, the horizontal range is (x_min, x_max) and the vertical range is (z_min, z_max).

[0036] Step S402: Based on the horizontal range and vertical range corresponding to all points, determine the minimum coverage rectangular grid range of the Gaussian beam; Figure 4 As shown, the minimum coverage rectangular grid range of the Gaussian beam is Figure 4The part within the dotted line.

[0037] Step S50: performing local extraction on the Gaussian beam travel time table based on the minimum coverage rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table.

[0038] In the example of this book, according to the minimum coverage rectangular grid range of the Gaussian beam, the number of horizontal grid points is x_max-x_min+1 and the number of vertical grid points is z_max-z_min+1. Then, the Gaussian beam travel time table is locally extracted to obtain the following: Figure 5 The orange area shown is the storage space of the optimized travel time table ( Figure 5 The orange area in the middle is significantly smaller than the overall underground grid space ( Figure 5 medium grey area).

[0039] Step S60: storing the optimized Gaussian beam travel time table.

[0040] As a further optimization of this embodiment, the optimized Gaussian beam travel time table is stored, including: Step S601: construct a two-dimensional array based on the horizontal range and vertical range corresponding to each point in the optimized Gaussian beam travel time table; wherein the two-dimensional array is a surface consisting of m points in the horizontal direction and n points in the vertical direction in the Gaussian beam travel time table.

[0041] Step S602: storing the optimized Gaussian beam travel time table based on a two-dimensional array.

[0042] In this embodiment, by applying the above steps S10 to S60 to Gaussian beams of different initial directions, the optimized storage of the Gaussian beam travel time table can be achieved, effectively reducing the storage cost of the Gaussian beam travel time table; Figure 6 and Figure 7 The optimized extracted results of Gaussian beam travel time tables for two other different initial directions show that they are significantly smaller than the overall imaging grid in the underground.

[0043] Therefore, the present invention can quickly determine the effective width of the Gaussian beam according to the propagation path of the central ray of the Gaussian beam and the starting shot point. Based on the effective width of the Gaussian beam, the minimum coverage rectangular grid range of the Gaussian beam can be dynamically determined. The Gaussian beam travel time table is extracted and stored based on the minimum coverage rectangular grid range of the Gaussian beam, which can greatly reduce the storage cost of the Gaussian beam travel time table.

[0044] Embodiment 2 Figure 8 1 is a block diagram of a Gaussian beam travel time table optimization storage system provided by an embodiment of the present invention. Figure 8As shown, this embodiment provides a Gaussian beam travel time table optimization storage system, the system is used to implement the Gaussian beam travel time table optimization storage method in embodiment 1, the system includes: A data acquisition module, used to obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; A path tracing module is used to perform ray tracing based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam; A width calculation module, used to determine the effective width of the Gaussian beam based on the propagation path and the starting shot point of the central ray of the Gaussian beam; A range calculation module, used for determining the minimum coverage rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam; A range extraction module is used to perform local extraction of the Gaussian beam travel time table based on the minimum coverage rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table; The travel time table storage module is used to store the optimized Gaussian beam travel time table.

[0045] This embodiment further provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned Gaussian beam travel time table optimization storage method when executing the computer program.

[0046] This embodiment further provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned Gaussian beam travel time table optimization storage method is implemented.

[0047] The present invention can quickly determine the effective width of the Gaussian beam according to the propagation path of the central ray of the Gaussian beam and the starting shot point. Based on the effective width of the Gaussian beam, the minimum coverage rectangular grid range of the Gaussian beam can be dynamically determined. The Gaussian beam travel time table is extracted and stored based on the minimum coverage rectangular grid range of the Gaussian beam, which can greatly reduce the storage cost of the Gaussian beam travel time table.

[0048] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0049] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A system that specifies the functions of a box or boxes.

[0050] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A Gaussian beam travel time table optimization storage method, characterized in that: The method comprises: Obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; Ray tracing is performed based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam; Based on the propagation path and the starting shot point of the central ray of the Gaussian beam, the effective width of the Gaussian beam is determined; Based on the effective width of the Gaussian beam, determine the minimum coverage rectangular grid range of the Gaussian beam; Based on the minimum coverage rectangular grid range of the Gaussian beam, the Gaussian beam travel time table is locally extracted to obtain the optimized Gaussian beam travel time table; The optimized Gaussian beam travel time table is stored.

2. The Gaussian beam travel time table optimization storage method according to claim 1, characterized in that: Ray tracing is performed based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam, including: At the starting shot point, ray tracing is performed according to the initial direction of the Gaussian beam, and the kinematic ray tracing partial differential equation is obtained; The kinematic ray tracing partial differential equation is numerically solved based on the Runge-Kutta algorithm to obtain the propagation path of the central ray of the Gaussian beam.

3. The Gaussian beam travel time table optimization storage method according to claim 1, characterized in that: Based on the propagation path and starting shot point of the central ray of the Gaussian beam, the effective width of the Gaussian beam is determined, including: Extract any point on the propagation path of the central ray of the Gaussian beam as an unknown point; Obtain the velocity of the unknown point, and determine the kinematic ray parameters of the unknown point based on the velocity of the unknown point; Obtain the seismic wave propagation travel time from the Gaussian beam travel time table; Determine the effective width of the unknown point based on the kinematic ray parameters of the unknown point, the travel time of seismic wave propagation, and the starting shot point; The effective width of the Gaussian beam is obtained by combining the effective widths corresponding to all points on the propagation path of the central ray of the Gaussian beam.

4. The Gaussian beam travel time table optimization storage method according to claim 3, characterized in that: The calculation expression of the kinematic ray parameters of the unknown point is: ； In the formula, represents the kinematic ray parameters of the unknown point P, represents the velocity of the unknown point P, and dt represents the time sampling interval of ray tracing.

5. The Gaussian beam travel time table optimization storage method according to claim 4, characterized in that: The calculation expression of the effective width of the unknown point is: ； In the formula, represents the effective width of the unknown point P, B represents the starting shot point, T represents the propagation time of the seismic wave, represents the kinematic ray parameters of the unknown point P, Represents the kinematic ray parameters of the starting shot point.

6. The Gaussian beam travel time table optimization storage method according to claim 3, characterized in that: Based on the effective width of the Gaussian beam, determine the minimum coverage rectangular grid range of the Gaussian beam, including: Traverse all points on the propagation path of the central ray of the Gaussian beam, and extract the horizontal range and vertical range corresponding to each point according to the effective width corresponding to the point; Based on the horizontal range and vertical range corresponding to all points, the minimum coverage rectangular grid range of the Gaussian beam is determined.

7. The Gaussian beam travel time table optimization storage method according to claim 6, characterized in that: The optimized Gaussian beam travel time table is stored, including: In the optimized Gaussian beam travel time table, a two-dimensional array is constructed based on the horizontal range and vertical range corresponding to each point; The optimized Gaussian beam travel time table is stored based on a two-dimensional array.

8. A Gaussian beam travel time table optimization storage system, used to implement the Gaussian beam travel time table optimization storage method according to any one of claims 1 to 7, characterized in that: The system comprises: A data acquisition module, used to obtain the starting shot point and the initial direction of the Gaussian beam in the Gaussian beam travel time table; A path tracing module is used to perform ray tracing based on the starting shot point and the initial direction of the Gaussian beam to obtain the propagation path of the central ray of the Gaussian beam; A width calculation module, used to determine the effective width of the Gaussian beam based on the propagation path and the starting shot point of the central ray of the Gaussian beam; A range calculation module, used for determining the minimum coverage rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam; A range extraction module is used to perform local extraction of the Gaussian beam travel time table based on the minimum coverage rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table; The travel time table storage module is used to store the optimized Gaussian beam travel time table.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the Gaussian beam travel time table optimization storage method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the Gaussian beam travel time table optimization storage method described in any one of claims 1 to 7 is implemented.

Images