A Method, System, Device and Medium for Optimizing the Storage of Gaussian Beam Travel-Time Tables

Abstract

The present invention relates to the field of storage technology, and discloses an optimized storage method, system, device and medium for Gaussian beam travel time tables. The method includes: obtaining the starting shot point in the Gaussian beam travel time table and the initial direction of the Gaussian beam; 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; 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; determining the minimum covering rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam; locally extracting the Gaussian beam travel time table based on the minimum covering 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. The present invention can significantly reduce the storage cost of the Gaussian beam travel time table.

Description

Technical Field

[0001] The present invention relates to the field of storage technology, and particularly to an optimized storage method, system, device and medium for Gaussian beam travel time tables. Background Art

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

[0003] The traditional Gaussian beam travel time table is defined based on the entire underground imaging grid space. The technical implementation process is as follows: First, calculate the central path of the Gaussian beam through ray tracing; Second, calculate the effective width of the Gaussian beam based on the central ray, and on this basis, determine the spread range of the Gaussian beam travel time table; Next, extract and store the entire imaging grid containing the Gaussian beam.

[0004] Since the effective range of the Gaussian beam is significantly smaller than the entire imaging grid, the traditional method will waste a large amount of storage space and seriously affect the calculation efficiency and scalability of the Gaussian beam migration imaging technology. Summary of the Invention

[0005] The purpose of the present invention is to provide an optimized storage method, system, device and medium for Gaussian beam travel time tables, so as to solve the problem that the traditional method will waste a large amount of storage space and seriously affect the calculation efficiency and scalability of the Gaussian beam migration imaging technology.

[0006] To achieve the above-mentioned invention purpose, the technical solution adopted by the present invention is:

[0007] In the first aspect, the present invention provides an optimized storage method for Gaussian beam travel time tables, and the method includes:

[0008] Obtain the starting shot point in the Gaussian beam travel time table and the initial direction of the Gaussian beam;

[0009] 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;

[0010] 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;

[0011] Determine the minimum covering rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam;

[0012] Perform local extraction on the Gaussian beam travel time table based on the minimum covering rectangular grid range of the Gaussian beam to obtain the optimized Gaussian beam travel time table;

[0013] Store the optimized Gaussian beam travel time table.

[0014] 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:

[0015] At the starting shot point, ray tracing is performed according to the initial direction of the Gaussian beam to obtain the kinematic ray tracing partial differential equation;

[0016] Based on the Runge-Kutta algorithm, numerical solution is performed on the kinematic ray tracing partial differential equation to obtain the propagation path of the central ray of the Gaussian beam.

[0017] Preferably, based on the propagation path of the central ray of the Gaussian beam and the starting shot point, the effective width of the Gaussian beam is determined, including:

[0018] Extract any point on the propagation path of the central ray of the Gaussian beam as the unknown point;

[0019] Obtain the velocity of the unknown point, and based on the velocity of the unknown point, determine the kinematic ray parameters of the unknown point;

[0020] Obtain the travel time of seismic wave propagation in the Gaussian beam travel time table;

[0021] Based on the kinematic ray parameters of the unknown point, the travel time of seismic wave propagation, and the starting shot point, determine the effective width of the unknown point;

[0022] Combine 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.

[0023] Preferably, the calculation expression for the kinematic ray parameters of the unknown point is:

[0024] ;

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

[0026] Preferably, the calculation expression for the effective width of the unknown point is:

[0027] ;

[0028] In the formula, represents the effective width of the unknown point P, B represents the starting shot point, T represents the travel time of seismic wave propagation, represents the kinematic ray parameters of the unknown point P, Represents the kinematic ray parameters of the starting shot point.

[0029] Preferably, based on the effective width of the Gaussian beam, determine the minimum covering rectangular grid range of the Gaussian beam, including:

[0030] Traverse all points on the propagation path of the central ray of the Gaussian beam, and extract the corresponding horizontal range and vertical range of each point according to the effective width corresponding to each point;

[0031] Based on the horizontal range and vertical range corresponding to all points, determine the minimum covering rectangular grid range of the Gaussian beam.

[0032] Preferably, store the optimized Gaussian beam travel time table, including:

[0033] In the optimized Gaussian beam travel time table, construct a two-dimensional array based on the horizontal range and vertical range corresponding to each point;

[0034] Store the optimized Gaussian beam travel time table based on the two-dimensional array.

[0035] In a second aspect, the present invention provides a Gaussian beam travel time table optimization storage system for implementing the above-mentioned Gaussian beam travel time table optimization storage method. The system includes:

[0036] A data acquisition module for acquiring the starting shot point in the Gaussian beam travel time table and the initial direction of the Gaussian beam;

[0037] A path tracking module for 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;

[0038] A width calculation module for 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;

[0039] A range calculation module for determining the minimum covering rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam;

[0040] A range extraction module for locally extracting the Gaussian beam travel time table based on the minimum covering rectangular grid range of the Gaussian beam to obtain the optimized Gaussian beam travel time table;

[0041] A travel time table storage module for storing the optimized Gaussian beam travel time table.

[0042] In a third aspect, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-mentioned Gaussian beam travel time table optimization storage method.

[0043] Fourthly, the present invention provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the above-mentioned optimized storage method for Gaussian beam travel time table is implemented.

[0044] The beneficial effects of the present invention are mainly reflected in:

[0045] Based on the propagation path of the central ray of the Gaussian beam and the starting shot point, the present invention can quickly determine the effective width of the Gaussian beam. Within the effective width of the Gaussian beam, the minimum covering rectangular grid range of the Gaussian beam can be dynamically determined. Based on the minimum covering rectangular grid range of the Gaussian beam, the Gaussian beam travel time table is extracted and stored, which can greatly reduce the storage cost of the Gaussian beam travel time table. Description of the Drawings

[0046] The 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 embodiments, they are used to explain the embodiments of the present invention, but do not constitute a limitation to the embodiments of the present invention. In the drawings:

[0047] Figure 1 is a flowchart of the optimized storage method for Gaussian beam travel time table provided by an embodiment of the present invention;

[0048] Figure 2 is a travel time schematic diagram of a Gaussian beam in the Gaussian beam travel time table provided by an embodiment of the present invention;

[0049] Figure 3 is a schematic diagram of the effective width of the Gaussian beam provided by an embodiment of the present invention;

[0050] Figure 4 is a schematic diagram of the minimum covering rectangular grid range of the Gaussian beam provided by an embodiment of the present invention;

[0051] Figure 5 is a comparative schematic diagram of the storage space of the optimized travel time table provided by an embodiment of the present invention;

[0052] Figure 6 is a schematic diagram of the optimized extraction result of the Gaussian beam travel time table with different initial directions provided by an embodiment of the present invention;

[0053] Figure 7 is a schematic diagram of the optimized extraction result of the Gaussian beam travel time table with another different initial direction provided by an embodiment of the present invention;

[0054] Figure 8 is a block diagram of the Gaussian beam travel time table optimized storage system provided by an embodiment of the present invention. Detailed Embodiments

[0055] To more clearly illustrate the technical solutions in the embodiments of the present invention or in the prior art, the present invention will be briefly introduced below in combination with the accompanying drawings and the descriptions of the embodiments or the prior art. Obviously, the following descriptions of the structures of the accompanying drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings. It should be noted here that the descriptions of these embodiments are used to help understand the present invention, but do not constitute a limitation to the present invention.

[0056] Embodiment 1

[0057] Figure 1 is a flowchart of an optimized storage method for Gaussian beam travel time tables provided by an embodiment of the present invention. As Figure 1 shown, this embodiment provides an optimized storage method for Gaussian beam travel time tables, and the method includes:

[0058] Step S10: Obtain the starting shot point in the Gaussian beam travel time table and the initial direction of the Gaussian beam; in this embodiment, as Figure 2 shown, Figure 2 is the travel time of a Gaussian beam in the Gaussian beam travel time table, and the starting shot point and the initial direction of the Gaussian beam can be extracted from the travel time.

[0059] Step S20: 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.

[0060] As a further optimization of this embodiment, 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 includes:

[0061] Step S201: At the starting shot point, perform ray tracing according to the initial direction of the Gaussian beam to obtain the kinematic ray tracing partial differential equation.

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

[0063] In this embodiment, the Runge-Kutta algorithm is a class of classical numerical methods for solving the initial value problems of ordinary differential equations (ODEs). Its core idea is to improve the calculation accuracy by means of weighted averaging of multi-stage slopes.

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

[0065] After the solution of steps S201 to S202, the following can be obtained Figure 2The propagation path of the blue curve part in

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

[0067] In this embodiment, the traditional method for calculating the Gaussian beam width requires extremely time-consuming dynamic ray tracing. To ensure the calculation efficiency, the calculation method for the effective width of the Gaussian beam is as follows:

[0068] Step S301: Extract any point on the propagation path of the central ray of the Gaussian beam as the unknown point.

[0069] 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; where the calculation expression for the kinematic ray parameters of the unknown point is:

[0070] ;

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

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

[0073] Step S304: Determine the effective width of the unknown point based on the kinematic ray parameters of the unknown point, the travel time of the seismic wave, and the starting shot point; where the calculation expression for the effective width of the unknown point is:

[0074] ;

[0075] In the formula, represents the effective width of the unknown point P, B represents the starting shot point, T represents the travel 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.

[0076] Step S305: Combine 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.

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

[0078] After the processing of steps S301 to S305, the obtained effective width of the Gaussian beam is asFigure 3 As shown Figure 3 The three red curves in are used to characterize the range of the effective width of three points, and the arc shape is used to characterize 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.

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

[0080] As a further optimization of this embodiment, determining the minimum covering rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam includes:

[0081] Step S401: Traverse all points on the propagation path of the central ray of the Gaussian beam, and extract the corresponding horizontal range and vertical range according to the effective width corresponding to each point. The horizontal range and vertical range are as Figure 4 shown. In Figure 4 , the horizontal range is (x_min, x_max), and the vertical range is (z_min, z_max).

[0082] Step S402: Determine the minimum covering rectangular grid range of the Gaussian beam based on the horizontal range and vertical range corresponding to all points; as Figure 4 shown, the minimum covering rectangular grid range of the Gaussian beam is Figure 4 the part within the dashed line in.

[0083] Step S50: Locally extract the travel-time table of the Gaussian beam based on the minimum covering rectangular grid range of the Gaussian beam to obtain an optimized travel-time table of the Gaussian beam.

[0084] In this example, according to the minimum covering rectangular grid range of the Gaussian beam, a horizontal grid point number of x_max - x_min + 1 and a vertical grid point number of z_max - z_min + 1 are established, and then the travel-time table of the Gaussian beam is locally extracted to obtain the orange area as shown in Figure 5 . The storage space of the optimized travel-time table ( Figure 5 the orange area in) is significantly smaller than the overall underground grid space ( Figure 5 the gray area in).

[0085] Step S60: Store the optimized travel-time table of the Gaussian beam.

[0086] As a further optimization of this embodiment, storing the optimized travel-time table of the Gaussian beam includes:

[0087] Step S601: In the optimized travel-time table of the Gaussian beam, construct a two-dimensional array based on the horizontal range and vertical range corresponding to each point; among them, the two-dimensional array is a plane, which is composed of m points in the horizontal direction and n points in the vertical direction in the travel-time table of the Gaussian beam.

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

[0089] In this embodiment, by applying the above steps S10 to S60 to Gaussian beams with 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 are the optimized extraction results of the Gaussian beam travel time tables for two other Gaussian beams with different initial directions, and it can be seen that they are significantly smaller than the overall imaging grid underground.

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

[0091] Embodiment 2

[0092] Figure 8 is a block diagram of an optimized storage system for a Gaussian beam travel time table provided by an embodiment of the present invention. As Figure 8 shown, this embodiment provides an optimized storage system for a Gaussian beam travel time table. The system is used to implement the optimized storage method of the Gaussian beam travel time table in Embodiment 1. The system includes:

[0093] A data acquisition module, configured to acquire the starting shot point in the Gaussian beam travel time table and the initial direction of the Gaussian beam.

[0094] A path tracking module, configured 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.

[0095] A width calculation module, configured to 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.

[0096] A range calculation module, configured to determine the minimum covering rectangular grid range of the Gaussian beam based on the effective width of the Gaussian beam.

[0097] A range extraction module, configured to locally extract the Gaussian beam travel time table based on the minimum covering rectangular grid range of the Gaussian beam to obtain an optimized Gaussian beam travel time table.

[0098] A travel time table storage module, configured to store the optimized Gaussian beam travel time table.

[0099] This embodiment also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above-mentioned Gaussian beam travel time table optimization storage method is implemented.

[0100] This embodiment also 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.

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

[0102] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can 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 code.

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

[0104] The above are only the embodiments of the present application and are not used to limit the present application. For those skilled in the art, various changes and modifications can be made to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall 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