Seismic surface wave tomography methods, devices, media, and program products

Abstract

The application discloses a kind of seismic surface wave tomography method, equipment, medium and program product, method includes: first, the seismic surface wave tomography study area is divided into uniform spherical equilateral triangle grid;For each first spherical equilateral triangle grid to be divided, the corresponding seismic ray density of first spherical equilateral triangle grid is calculated, the first spherical equilateral triangle grid with greater density is divided, and the second spherical equilateral triangle grid of completed division is obtained;Seismic ray tracing and inversion are carried out based on the second spherical equilateral triangle grid, and the imaging result of the seismic surface wave tomography study area is formed.The embodiment of the application realizes the matching of seismic ray distribution density and tomography pixel by the seismic surface wave tomography of variable size grid based on seismic ray distribution driving, and effectively improves the actual resolution capability of seismic surface wave tomography.

Description

Technical Field

[0001] The present invention relates to the technical field of seismic tomography, and in particular to a seismic surface wave tomography method, equipment, medium and program product. Background Art

[0002] Seismic surface wave tomography (SWT) uses the principles of medical CT to image the Earth's internal structure. The difference lies in the fact that the "X-rays" used are seismic surface waves traveling between the earthquake source and the observation station (hereinafter referred to as seismic rays for ease of understanding). Similarities exist in that both methods divide the scanned area into many small cells, determining the properties of each cell through inversion calculations to achieve imaging of the scanned area. Therefore, to perform SWT in a given area, sufficient earthquakes, observation stations, and corresponding seismic rays must be present.

[0003] In medical CT, X-rays are actively generated by humans, ensuring uniform coverage of the scanned lesion, so the size of the divided units can also be uniform. However, for seismic surface wave tomography, the location of earthquakes is uncontrollable. While the location of station installation is generally controllable, it is occasionally restricted by the geographical environment and cannot be uniformly deployed, making it difficult for seismic rays to pass evenly through the study area. However, most current seismic tomography methods typically use a uniform square grid with a fixed spacing that matches the average data coverage level of the study area, ignoring the uneven distribution of data. The advantage of this grid definition method is that it is simple and easy to perform routine processing such as calculating surface wave travel times, adding regularization constraints, and testing the detection plate resolution. The disadvantage is that it weakens the lateral resolution ability in areas with advantageous data distribution or strengthens the lateral resolution ability in areas with disadvantageous data distribution.

[0004] The existing technology usually divides the study area into Figure 2 The equidistant square grid shown in Figure 1 is used to trace each ray according to the grid definition and establish the propagation equation for each ray, for example: Figure 2 The propagation time equation of the ray is , the propagation time equation of ray B is , each ray can establish the above linear equation. Assuming there are m rays, we can obtain a system of equations consisting of m equations and n unknowns (n=30): , where the propagation time (T) is known and the propagation path ( ) we can see that the grid speed ( ) to be determined; finally, by inverting the above linear equations, the velocity of each grid is obtained ( ).

[0005] From the above existing technical solutions, it can be seen that the definition of the equally spaced square grid is independent of the ray distribution, which means that the grid velocity obtained by inversion is uniformly distributed, which can be understood as the uniform distribution of the imaged pixels. Figure 3 It can be seen that in actual situations, the distribution of seismic rays is very uneven, and the resolution of tomography is closely related to the distribution of seismic rays. The greater the density of ray distribution, the more effective linear equations can be established, the more reliable the grid speed of the inversion calculation, and the higher the resolution. The grid size should be reduced to increase the pixels; on the contrary, the smaller the density of seismic ray distribution, the less reliable the inversion grid speed, and the lower the resolution. The grid size should be increased to reduce the pixels.

[0006] Therefore, existing technologies ignore the correlation between seismic ray distribution and imaging resolution. Imaging with a uniformly sized grid (i.e., fixed pixels) will underestimate the actual resolution capability in areas with good ray coverage and overestimate the actual resolution capability in areas with poor ray coverage, resulting in imaging distortion in local areas.

[0007] Therefore, it is urgent to invent a variable-size grid seismic surface wave tomography method to solve the problem that the existing seismic surface wave tomography method ignores the correlation between seismic ray distribution and imaging resolution, which leads to underestimation of the actual resolution capability in areas with good ray coverage and overestimation of the actual resolution capability in areas with poor ray coverage, resulting in imaging distortion in local areas. Summary of the Invention

[0008] In view of this, embodiments of the present invention provide a seismic surface wave tomography method, device, medium, and program product, which at least partially solve the problems existing in the prior art.

[0009] Other features and advantages of the present invention will become apparent from the following detailed description, or may be learned in part by practice of the present invention.

[0010] In order to achieve the above objectives, the embodiments of the present invention provide the following technical solutions:

[0011] According to a first aspect of an embodiment of the present invention, a seismic surface wave tomography method is provided, the method comprising:

[0012] The seismic surface wave tomography study area is divided into uniform spherical equilateral triangle grids;

[0013] For each first spherical equilateral triangle grid to be segmented, calculating the seismic ray density corresponding to the first spherical equilateral triangle grid;

[0014] dividing the first spherical equilateral triangular grid according to the seismic ray density corresponding to the first spherical equilateral triangular grid to obtain a divided second spherical equilateral triangular grid;

[0015] Seismic ray tracing and inversion are performed based on the equilateral triangle grid of the second sphere after segmentation to form the imaging results of the seismic surface wave tomography study area.

[0016] Furthermore, calculating the seismic ray density corresponding to the first spherical equilateral triangle grid includes:

[0017] Based on a uniform spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each spherical equilateral triangle grid;

[0018] Obtaining, based on the propagation path lengths of the seismic rays in the equilateral triangular grids on the spherical surface, a statistical count of the number of times the first equilateral triangular grid is traversed by different seismic rays;

[0019] The seismic ray density corresponding to the first spherical equilateral triangle grid is obtained according to the number of times the first spherical equilateral triangle grid is passed through by different seismic rays.

[0020] Furthermore, the first spherical equilateral triangular mesh is segmented according to the seismic ray density corresponding to the first spherical equilateral triangular mesh to obtain a segmented second spherical equilateral triangular mesh, including:

[0021] For each first spherical equilateral triangle mesh to be segmented, determining whether a ray density corresponding to the first spherical equilateral triangle mesh is greater than a preset density threshold;

[0022] If the ray density corresponding to the first spherical equilateral triangle mesh is less than or equal to a preset density threshold, directly using the first spherical equilateral triangle mesh as the second spherical equilateral triangle mesh completed by segmentation;

[0023] If the ray density corresponding to the first spherical equilateral triangle mesh is greater than a preset density threshold, dividing the first spherical equilateral triangle mesh into four first spherical equilateral triangle meshes of equal area;

[0024] For the first spherical equilateral triangle grid generated by equal division, a loop is performed to calculate the seismic ray density corresponding to the first spherical equilateral triangle grid, and it is determined whether the ray density corresponding to the first spherical equilateral triangle grid is greater than a preset density threshold.

[0025] Furthermore, seismic ray tracing and inversion are performed based on the segmented second spherical equilateral triangle mesh to form the imaging results of the seismic surface wave tomography study area, including:

[0026] Perform seismic ray tracing based on each second spherical equilateral triangle grid to obtain ray tracing results;

[0027] Constructing a ray propagation linear equation system based on the ray tracing results and performing inversion processing to obtain grid velocities corresponding to each second spherical equilateral triangle grid;

[0028] The imaging results of the seismic surface wave tomography study area are formed according to the grid velocities corresponding to all the equilateral triangle grids of the second sphere.

[0029] Furthermore, seismic ray tracing is performed based on each second spherical equilateral triangle grid to obtain ray tracing results, including:

[0030] Based on the second spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each second spherical equilateral triangle grid.

[0031] Furthermore, a ray propagation linear equation group is constructed based on the ray tracing results and inversion processing is performed to obtain the mesh velocities corresponding to each second spherical equilateral triangle mesh, including:

[0032] According to the propagation path length of each seismic ray in each equilateral triangle grid of the second sphere and the corresponding propagation time, a linear equation group of seismic ray propagation velocity is constructed;

[0033] The linear equations for the seismic ray propagation velocity are inverted and calculated to obtain the grid velocities corresponding to the equilateral triangle grids on the second spherical surface.

[0034] According to a second aspect of an embodiment of the present invention, there is provided a device, comprising: a processor and a memory;

[0035] The memory is used to store one or more program instructions;

[0036] The processor is used to run one or more program instructions to execute the steps of a seismic surface wave tomography method as described in any one of the above items.

[0037] According to a third aspect of an embodiment of the present invention, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps of a seismic surface wave tomography method as described in any one of the above items are implemented.

[0038] According to a fourth aspect of an embodiment of the present invention, a computer program product is provided, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer implements the steps of a seismic surface wave tomography method as described in any one of the above items.

[0039] The seismic surface wave tomography method, device, medium, and program product provided by the embodiments of the present invention have the following advantages:

[0040] The existing technology based on uniform square grid ray tracing is usually simplified to straight line propagation on a plane, involving only plane geometry theory. The embodiment of the present invention adopts a ray tracing method of seismic rays passing through a spherical grid, involving complex spherical geometry theory, which is more in line with the actual situation of seismic surface waves propagating along the earth's surface than the existing technology.

[0041] Based on the characteristic that a spherical equilateral triangular grid can be easily divided infinitely into four equilateral triangular grids of equal area, the embodiment of the present invention determines the number of grid divisions by the grid ray density. A grid with a low ray density is no longer divided, while a grid with a high ray density can be divided once or multiple times, thereby achieving the purpose of matching the ray density with the grid size, and ultimately achieving the technical effect of matching the seismic ray distribution density with the tomographic imaging pixels. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] To more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for the embodiments or the description of the prior art. Obviously, the drawings described below are merely exemplary, and those skilled in the art can derive other implementation drawings based on the provided drawings without inventive effort.

[0043] Figure 1 A schematic flow chart of a seismic surface wave tomography method provided by an embodiment of the present invention;

[0044] Figure 2 This is a schematic diagram of the effect of grid division of the study area in the prior art of the technical field to which the embodiment of the present invention belongs;

[0045] Figure 3 Schematic diagrams comparing seismic ray distribution in the technical field to which the embodiments of the present invention pertain, wherein the left diagram shows an ideal uniform ray distribution, and the right diagram shows an actual non-uniform ray distribution;

[0046] Figure 4 A schematic diagram of the effect of a uniform spherical equilateral triangle grid on a study area of ​​a seismic surface wave tomography method provided by an embodiment of the present invention;

[0047] Figure 5 A schematic diagram illustrating the effect of seismic ray density on a spherical equilateral triangle grid of uniform size in a seismic surface wave tomography method provided by an embodiment of the present invention;

[0048] Figure 6 A schematic diagram of the effect of a non-uniform spherical equilateral triangle mesh after being segmented once in a seismic surface wave tomography method provided by an embodiment of the present invention;

[0049] Figure 7A schematic diagram comparing tomographic imaging effects of a seismic surface wave tomography method provided by an embodiment of the present invention;

[0050] Figure 8 A schematic diagram of segmenting a spherical equilateral triangle grid for a seismic surface wave tomography method provided by an embodiment of the present invention. DETAILED DESCRIPTION

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

[0052] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that the numbers used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0053] Figure 1 The figure shows a flow chart of a seismic surface wave tomography method according to an embodiment of the present invention.

[0054] like Figure 1 As shown, the seismic surface wave tomography method according to an embodiment of the present invention may include step S100, step S200, step S300 and step S400.

[0055] In step S100 , the seismic surface wave tomography study area is divided into uniform spherical equilateral triangle grids.

[0056] Figure 4 A schematic diagram showing the effect of a uniform spherical equilateral triangle grid on a study area of ​​a seismic surface wave tomography method provided by an embodiment of the present invention is shown.

[0057] Next, in step S200 , for each first spherical equilateral triangle mesh to be segmented, the seismic ray density corresponding to the first spherical equilateral triangle mesh is calculated.

[0058] Specifically, the above steps include:

[0059] Based on a uniform spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each spherical equilateral triangle grid.

[0060] According to the propagation path lengths in each spherical equilateral triangular grid, the number of times the first spherical equilateral triangular grid is penetrated by different seismic rays is obtained by counting. According to the number of times the first spherical equilateral triangular grid is penetrated by different seismic rays, the seismic ray density corresponding to the first spherical equilateral triangular grid is converted.

[0061] Figure 5 A schematic diagram showing the seismic ray density effect of a first spherical equilateral triangle grid of uniform size in a seismic surface wave tomography method provided by an embodiment of the present invention is shown. Figure 5 The ray density on the right side of the study area is greater than that on the left side.

[0062] Next, in step S300, the first spherical equilateral triangular mesh is segmented according to the seismic ray density corresponding to the first spherical equilateral triangular mesh to obtain a segmented second spherical equilateral triangular mesh.

[0063] Specifically, the above steps include:

[0064] For each first spherical equilateral triangle mesh to be segmented, it is determined whether a ray density corresponding to the first spherical equilateral triangle mesh is greater than a preset density threshold.

[0065] If the ray density corresponding to the first spherical equilateral triangle mesh is less than or equal to the preset density threshold, the first spherical equilateral triangle mesh is directly used as the second spherical equilateral triangle mesh completed by segmentation.

[0066] If the ray density corresponding to the first spherical equilateral triangle mesh is greater than the preset density threshold, the first spherical equilateral triangle mesh is divided into four first spherical equilateral triangle meshes of the same area. Figure 8 A schematic diagram of segmenting a spherical equilateral triangle grid in a seismic surface wave tomography method provided by an embodiment of the present invention is shown.

[0067] For the first spherical equilateral triangle grid generated by the equal division, a loop is performed to calculate the seismic ray density corresponding to the first spherical equilateral triangle grid, and it is determined whether the ray density corresponding to the first spherical equilateral triangle grid is greater than a preset density threshold.

[0068] Figure 6 A schematic diagram showing the effect of a non-uniform spherical equilateral triangle grid after being segmented once in a seismic surface wave tomography method provided by an embodiment of the present invention is shown.

[0069] pass Figure 6 It can be seen that after performing non-uniform grid segmentation using a seismic surface wave tomography method provided by an embodiment of the present invention, the ray density of different regions can be well reflected.

[0070] Next, in step S400, seismic ray tracing and inversion are performed based on the segmented second spherical equilateral triangle mesh to form an imaging result of the seismic surface wave tomography study area.

[0071] Specifically, the above steps include:

[0072] Based on the final segmented second spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each second spherical equilateral triangle grid.

[0073] According to the propagation path length of each seismic ray in each second spherical equilateral triangle grid and the corresponding propagation time, a linear equation group of seismic ray propagation velocity is constructed for each ray and each grid attribute.

[0074] The linear equations of seismic ray propagation velocity are inverted to obtain the grid velocities corresponding to each equilateral triangle grid on the second sphere.

[0075] The imaging results of the seismic surface wave tomography study area are formed according to the grid velocities corresponding to all the equilateral triangle grids of the second sphere.

[0076] Figure 7 A schematic diagram showing a comparison of tomographic imaging effects of a seismic surface wave tomography method provided by an embodiment of the present invention is shown.

[0077] Figure 7 (a) Tomography results based on uniform-sized grids.

[0078] Figure 7 (b) is the tomographic imaging result after a non-uniform grid segmentation using an embodiment of the present invention. Figure 7 (b) contains two different grid sizes (pixels), Figure 7 (b) The lower left area has a large grid size (low pixels) Figure 7 (a) The ray sparse area, Figure 7 The upper right area of ​​(b) has a small grid size (low pixels) corresponding to Figure 7 (a) Ray-dense area.

[0079] Figure 7 (c) To utilize the embodiment of the present invention in Figure 7 (b) The tomographic imaging results after non-uniform grid segmentation are performed again. Figure 7(c) Includes three different grid sizes (pixels). It can be seen that the large, medium and small grids (pixels) are distributed in the lower left, middle and upper right areas respectively, and the image resolution is greatly improved.

[0080] In addition, an embodiment of the present invention also provides a device, which includes: a processor and a memory; the memory is used to store one or more program instructions; the processor is used to run one or more program instructions to execute the steps of a seismic surface wave tomography method as described above.

[0081] In addition, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the seismic surface wave tomography method described above are implemented.

[0082] In addition, an embodiment of the present invention further provides a computer program product, which includes computer program instructions. When the computer program instructions are executed by a processor, the steps of the seismic surface wave tomography method described above are implemented.

[0083] The seismic surface wave tomography method, device, medium, and program product provided by the embodiments of the present invention have the following advantages:

[0084] Seismic ray tracing based on a spherical equilateral triangular grid is realized. The existing technology of ray tracing based on a uniform square grid is usually simplified to straight-line propagation on a plane, involving only plane geometry theory. The embodiment of the present invention adopts a ray tracing method of seismic rays passing through a spherical grid, involving complex spherical geometry theory, which is more in line with the actual situation of seismic surface waves propagating along the earth's surface than the existing technology.

[0085] Based on the characteristic that a spherical equilateral triangular grid can be easily divided infinitely into four equilateral triangular grids of equal area, the embodiment of the present invention determines the number of grid divisions by the grid ray density. A grid with a low ray density is no longer divided, while a grid with a high ray density can be divided once or multiple times. This effectively optimizes the actual lateral resolution capability, thereby achieving the purpose of matching the ray density distribution density with the grid size, and ultimately achieving the technical effect of matching the seismic ray distribution density with the tomographic imaging pixels.

[0086] In the embodiments of the present invention, the processor can be an integrated circuit chip with signal processing capabilities. The processor can be a general-purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA), or other programmable logic device, discrete gate or transistor logic device, or discrete hardware component. It can implement or execute the various methods, steps, and logic block diagrams disclosed in the embodiments of the present invention. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the methods disclosed in the embodiments of the present invention can be directly implemented and executed by a hardware decoding processor, or by a combination of hardware and software modules within the decoding processor. The software module can be located in a storage medium well-known in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, or registers. The processor reads information from the storage medium and, in conjunction with its hardware, completes the steps of the aforementioned method. The storage medium can be a memory, for example, volatile memory or non-volatile memory, or can include both volatile and non-volatile memory. Non-volatile memory can be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), or flash memory. Volatile memory can be random access memory (RAM), which is used as an external cache. By way of example and not limitation, many forms of RAM are available, such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous dynamic random access memory (SDRAM), double data rate synchronous dynamic random access memory (DDRSDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link dynamic random access memory (SLDRAM), and direct rambus random access memory (DRRAM).The storage media described in the embodiments of the present invention are intended to include, but are not limited to, these and any other suitable types of memory. Those skilled in the art will appreciate that in one or more of the above examples, the functions described in the present invention can be implemented using a combination of hardware and software. When software is used, the corresponding functions can be stored in a computer-readable medium or transmitted as one or more instructions or codes on a computer-readable medium. Computer-readable media include computer storage media and communication media, wherein communication media includes any medium that facilitates the transmission of computer programs from one place to another. The storage medium can be any available medium that can be accessed by a general or special-purpose computer. Although the present invention has been described in detail above using general instructions and specific embodiments, it will be apparent to those skilled in the art that modifications or improvements can be made based on the present invention. Therefore, these modifications or improvements made without departing from the spirit of the present invention are within the scope of protection claimed in the present invention.

[0087] The above description is only a preferred embodiment of the present invention and does not limit the present invention in any form. Those skilled in the art can make some simple modifications, equivalent changes or modifications based on the technical content disclosed above, which all fall within the scope of protection of the present invention.

Claims

1. A seismic surface wave tomography method, characterized in that: The method comprises: The seismic surface wave tomography study area is divided into uniform spherical equilateral triangle grids; For each first spherical equilateral triangle grid to be segmented, calculating the seismic ray density corresponding to the first spherical equilateral triangle grid; The first spherical equilateral triangular mesh is segmented according to the seismic ray density corresponding to the first spherical equilateral triangular mesh to obtain a segmented second spherical equilateral triangular mesh, including: For each first spherical equilateral triangle mesh to be segmented, determining whether a ray density corresponding to the first spherical equilateral triangle mesh is greater than a preset density threshold; If the ray density corresponding to the first spherical equilateral triangle mesh is less than or equal to a preset density threshold, directly using the first spherical equilateral triangle mesh as the second spherical equilateral triangle mesh completed by segmentation; If the ray density corresponding to the first spherical equilateral triangle mesh is greater than a preset density threshold, dividing the first spherical equilateral triangle mesh into four first spherical equilateral triangle meshes of equal area; For the first spherical equilateral triangle grid generated by the equal division, looping to calculate the seismic ray density corresponding to the first spherical equilateral triangle grid, and determining whether the ray density corresponding to the first spherical equilateral triangle grid is greater than a preset density threshold; Seismic ray tracing and inversion are performed based on the equilateral triangle grid of the second sphere after segmentation to form the imaging results of the seismic surface wave tomography study area.

2. A seismic surface wave tomography method according to claim 1, characterized in that: Calculating the seismic ray density corresponding to the first spherical equilateral triangle grid includes: Based on a uniform spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each spherical equilateral triangle grid; Obtaining, based on the propagation path lengths of the seismic rays in the equilateral triangular grids on the spherical surface, a statistical count of the number of times the first equilateral triangular grid is traversed by different seismic rays; The seismic ray density corresponding to the first spherical equilateral triangle grid is obtained according to the number of times the first spherical equilateral triangle grid is passed through by different seismic rays.

3. The seismic surface wave tomography method according to claim 1, characterized in that: Seismic ray tracing and inversion are performed based on the segmented second sphere equilateral triangle grid to generate the imaging results of the seismic surface wave tomography study area, including: Perform seismic ray tracing based on each second spherical equilateral triangle grid to obtain ray tracing results; Constructing a ray propagation linear equation group based on the ray tracing results and performing inversion processing to obtain grid velocities corresponding to each second spherical equilateral triangle grid; The imaging results of the seismic surface wave tomography study area are formed according to the grid velocities corresponding to all the equilateral triangle grids of the second sphere.

4. A seismic surface wave tomography method according to claim 3, characterized in that: Seismic ray tracing is performed based on each second spherical equilateral triangle grid to obtain ray tracing results, including: Based on the second spherical equilateral triangle grid, ray tracing is performed on each seismic ray in the study area to obtain the propagation path length of each seismic ray in each second spherical equilateral triangle grid.

5. The seismic surface wave tomography method according to claim 4, characterized in that: According to the ray tracing results, a ray propagation linear equation group is constructed and inversion processing is performed to obtain the grid velocity corresponding to each second spherical equilateral triangle grid, including: According to the propagation path length of each seismic ray in each equilateral triangle grid of the second sphere and the corresponding propagation time, a linear equation group of seismic ray propagation velocity is constructed; The linear equations for the seismic ray propagation velocity are inverted and calculated to obtain the grid velocities corresponding to the equilateral triangle grids on the second spherical surface.

6. A seismic surface wave tomography device, characterized in that: The device includes: a processor and a memory; The memory is used to store one or more program instructions; The processor is used to run one or more program instructions to execute the steps of the seismic surface wave tomography method according to any one of claims 1 to 5.

7. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the seismic surface wave tomography method according to any one of claims 1 to 5 are implemented.

8. A computer program product, characterized in that The computer program product comprises computer program instructions, which, when executed by a processor, implement the steps of a seismic surface wave tomography method according to any one of claims 1 to 5.

Images