Sub-azimuth anisotropic first arrival ray tracing method, device, equipment and medium

Abstract

The invention provides a sub-azimuth anisotropic first-arrival ray tracing method, device and equipment and a medium, and belongs to the field of first-arrival ray tracing. The method comprises the steps that grid discrete processing is carried out on an HTI medium, an HTI rectangular area is obtained, the HTI rectangular area comprises a plurality of rectangular units, and each rectangular unit comprises a plurality of rays; determining a group angle and a length of each ray in the rectangular unit of the receiving point; inputting the preset work area data and the group angle of each ray corresponding to the receiving point into a pre-constructed sub-azimuth anisotropic velocity initial model to obtain a ray group velocity of a rectangular unit of the receiving point; determining travel time of the receiving point rectangular unit according to the ray group velocity of the receiving point rectangular unit and the length of each ray; and based on a wavefront scanning search algorithm, tracking a first arrival ray path according to travel time of a receiving point rectangular unit.

Description

Technical Field

[0001] The present invention relates to the field of first-arrival ray tracing, and in particular to an azimuthal anisotropic first-arrival ray tracing method, an azimuthal anisotropic first-arrival ray tracing device, an electronic device and a computer-readable storage medium. Background Art

[0002] With the continuous increase in exploration efforts, the exploration targets are getting deeper and deeper, the height difference of the work area is getting larger and larger, and the lateral changes in shallow velocity are becoming more dramatic; with the continuous development of new acquisition technologies, the proportion of three-dimensional seismic data acquired with large offsets and wide azimuths is increasing.

[0003] The anisotropic characteristics are reflected in these data. In many work areas, it has been found that there are obvious differences in the shallow velocity models inverted from first arrivals at different azimuths. The differences in the low-velocity part are particularly obvious and directly affect the imaging effect.

[0004] However, in the prior art, the ray path obtained only by isotropic first-arrival ray tracing cannot accurately characterize the characteristics of the azimuthally anisotropic velocity field, resulting in large differences in the results of tomographic inversion in different azimuths and an inability to better describe the near-surface velocity distribution.

[0005] How to better utilize large-offset, wide-azimuth seismic data to complete the anisotropy investigation and modeling of refraction waves, eliminate the influence of surface anisotropy, and establish a more accurate near-surface velocity model is an issue that needs to be addressed urgently. Summary of the invention

[0006] The purpose of the present invention is to provide a method, device, equipment and medium for tracing azimuthally anisotropic first-arrival rays, so as to at least solve the problem that the prior art cannot accurately characterize the characteristics of azimuthally anisotropic velocity fields.

[0007] In order to achieve the above object, on the one hand, the present invention provides an azimuthal anisotropic first-arrival ray tracing method, the method comprising: Performing grid discretization processing on the HTI medium to obtain an HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays; Determine the group angle and length of each ray within the rectangular unit of the receiving point; The preset work area data and the group angle of each ray corresponding to the receiving point are input into the pre-constructed azimuthal anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameter and vertical velocity; Determine the travel time of the rectangular unit at the receiving point according to the ray group velocity of the rectangular unit at the receiving point and the length of each ray; Based on the wavefront scanning search algorithm, the first arrival ray path is traced according to the travel time of the rectangular unit of the receiving point.

[0008] Preferably, based on a wavefront scanning search algorithm, tracing the first arrival ray path according to the travel time of the rectangular unit of the receiving point includes: Perform wavefront scanning on the HTI rectangular area and determine the minimum travel time of all rectangular units based on the travel time of the rectangular units at the receiving point; Search for the first arrival ray path based on the minimum travel time of all rectangular cells.

[0009] Preferably, the method further comprises: constructing an initial azimuthal anisotropic velocity model, comprising: Determine the anisotropic Thomsen parameters of VTI media; Performing coordinate axis conversion on the HTI medium and the VTI medium to obtain conversion results, wherein the conversion results include: anisotropic Thomsen parameters of the HTI medium and a conversion relationship between the HTI medium and the VTI medium; The initial model of azimuthal anisotropic velocity is constructed based on the conversion results.

[0010] Preferably, the azimuthal anisotropic velocity initial model for:

[0011] ;

[0012] In the formula, is the vertical speed, is the conversion relationship between HTI medium and VTI medium, is the first anisotropy parameter of the VTI medium, which is used to characterize the longitudinal anisotropy strength. is the second anisotropy parameter of VTI medium.

[0013] Preferably, the conversion relationship between the HTI medium and the VTI medium is: ;

[0014] In the formula, The axis of symmetry is axial, is the azimuth of the HTI medium, is the angle between the incident ray and the vertical.

[0015] Preferably, the anisotropic Thomsen parameters of the HTI medium include: a first anisotropic parameter of the HTI medium, and the calculation expression of the first anisotropic parameter of the HTI medium is: ;

[0016] In the formula, is the first anisotropic parameter of HTI medium, is the first anisotropy parameter of the VTI medium.

[0017] Preferably, the anisotropic Thomsen parameter of the HTI medium further includes: a second anisotropic parameter of the HTI medium, and the calculation expression of the second anisotropic parameter of the HTI medium is: ;

[0018] In the formula, is the second anisotropy parameter of HTI medium, is the second anisotropy parameter of the VTI medium, f is the frequency;

[0019] in, , where is the shear wave velocity.

[0020] On the other hand, the present invention further provides an azimuthally anisotropic first-arrival ray tracing device for implementing the above-mentioned azimuthally anisotropic first-arrival ray tracing method, the device comprising: A discrete processing module, used for performing grid discrete processing on the HTI medium to obtain an HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays; A ray determination module, used to determine the group angle and length of each ray within the rectangular unit of the receiving point; The velocity determination module is used to input the preset work area data and the group angle of each ray corresponding to the receiving point into the pre-built azimuth anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameter and vertical velocity; A travel time determination module, used to determine the travel time of the receiving point rectangular unit according to the ray group velocity of the receiving point rectangular unit and the length of each ray; The ray tracing module is used to trace the first arrival ray path according to the travel time of the rectangular unit of the receiving point based on the wavefront scanning search algorithm.

[0021] On the other hand, the present invention also 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 azimuthal anisotropic first-arrival ray tracing method when executing the computer program.

[0022] On the other hand, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned azimuthal anisotropic first-arrival ray tracing method.

[0023] Through the above technical solution, the present invention has at least the following technical effects: The present invention can improve the adaptability of the first-arrival wave ray tracing technology. In the collection area of ​​large offset and wide azimuth data, it has the characteristics of azimuthal anisotropy, can more accurately describe the ray path of the first-arrival wave, and plays a key role in near-surface anisotropy modeling.

[0024] Other features and advantages of the present invention will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] 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 an azimuthal anisotropic first-arrival ray tracing method provided by one embodiment of the present invention; Figure 2 is a schematic diagram of a HTI rectangular area provided by one embodiment of the present invention; Figure 3 It is a schematic diagram of moving from one node to another node under a uniform medium in a HTI rectangular area provided by an embodiment of the present invention; Figure 4 It is a block diagram of an azimuthal anisotropic first-arrival ray tracing device provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0026] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here is only used to illustrate and explain the embodiment of the present invention, and is not used to limit the embodiment of the present invention.

[0027] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by technicians in the technical field to which this application belongs; the terms used herein are only for the purpose of describing specific embodiments and are not intended to limit this application; the terms "including" and "having" in the specification and claims of this application and the above-mentioned figure descriptions and any variations thereof are intended to cover non-exclusive inclusions.

[0028] In the description of the embodiments of the present application, the technical terms "first", "second", etc. are only used to distinguish different objects, and cannot be understood as indicating or implying relative importance or implicitly indicating the number, specific order or primary and secondary relationship of the indicated technical features. In the description of the embodiments of the present application, the meaning of "multiple" is more than two, unless otherwise clearly and specifically defined.

[0029] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0030] Embodiment 1 Figure 1 is a flow chart of an azimuthal anisotropic first-arrival ray tracing method provided by one embodiment of the present invention, such as Figure 1 As shown, this embodiment provides an azimuthal anisotropic first-arrival ray tracing method, the method comprising: Step S101: performing grid discretization processing on the HTI medium to obtain a HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays.

[0031] In this embodiment, the HTI medium is discretized into (m*n*l) adjacent rectangular units of the same size through a grid, and the velocity is uniform in each unit, such as Figure 2 and Figure 3 shown.

[0032] In the subsequent azimuthal anisotropic ray tracing process, starting from a known travel time point (such as the source point), the minimum travel time and ray direction information are gradually calculated according to the Fermat principle; let Q be the set of known travel time points q, p be the unknown travel time point adjacent to it, t q ,t p are the minimum travel times at points q and p, respectively, and t pq is the minimum travel time from q to p (local travel time), r is the position of the secondary source of p (the intersection of the ray passing through point p and the wavefront), at this time: ;

[0033] According to Huygens' principle, q only needs to traverse the boundary of Q (i.e., the wavefront points). When the minimum travel time of all neighboring wavefront points is found, these points become new wavefront points.

[0034] Step S102: Determine the group angle and length of each ray in the rectangular unit of the receiving point.

[0035] exist Figure 3 In the geometric relationship of the ray path at the receiving point, the group angle corresponding to the ray path can be calculated ; The length of the ray is the length of the ray experienced in the current grid (the rectangular unit corresponding to the receiving point) .

[0036] Step S103: Input the preset work area data and the group angle of each ray corresponding to the receiving point into the pre-constructed azimuthal anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameters and vertical velocity.

[0037] In this embodiment, the method further includes: constructing an initial azimuthal anisotropic velocity model, including: Step a01: Determine the anisotropic Thomsen parameters of the VTI medium; Step a02: performing coordinate axis conversion on the HTI medium and the VTI medium to obtain a conversion result, wherein the conversion result includes: anisotropic Thomsen parameters of the HTI medium and a conversion relationship between the HTI medium and the VTI medium; Step a03: Construct an initial azimuthal anisotropic velocity model based on the conversion results.

[0038] Among them, the initial model of azimuthal anisotropic velocity for:

[0039] (1);

[0040] In the formula, is the vertical speed, is the conversion relationship between HTI medium and VTI medium, is the first anisotropy parameter of the VTI medium, which is used to characterize the longitudinal anisotropy strength. is the second anisotropy parameter of VTI medium.

[0041] The conversion relationship between the HTI medium and the VTI medium is: (2);

[0042] In the formula, The axis of symmetry is axial, is the azimuth of the HTI medium, is the angle between the incident ray and the vertical.

[0043] The anisotropic Thomsen parameters of the HTI medium include: a first anisotropic parameter of the HTI medium, and the calculation expression of the first anisotropic parameter of the HTI medium is: (3);

[0044] In the formula, is the first anisotropic parameter of HTI medium, is the first anisotropy parameter of the VTI medium.

[0045] The anisotropic Thomsen parameter of the HTI medium further includes: a second anisotropic parameter of the HTI medium, and the calculation expression of the second anisotropic parameter of the HTI medium is: (4);

[0046] In the formula, is the second anisotropy parameter of HTI medium, is the second anisotropy parameter of the VTI medium, f is the frequency;

[0047] in, , where is the shear wave velocity.

[0048] In this embodiment, the anisotropic parameters of the HTI medium include: the second anisotropic parameter of the HTI medium and the first anisotropic parameter of the HTI medium. According to the actual working condition data of the receiving point, the preset working area data can be obtained, that is, the azimuth angle of the HTI medium, the axial direction of the symmetry axis, the anisotropic parameter of the HTI medium and the vertical velocity can be obtained; When the anisotropy parameters of the HTI medium, the azimuth angle of the HTI medium, the axial direction of the symmetry axis, the anisotropy parameters of the HTI medium and the vertical velocity are known, the group angle corresponding to the ray path will be calculated. As According to formula (1), the group velocity of the current grid can be calculated. The group velocity of each ray in the same grid unit is the same.

[0049] Step S104: Determine the travel time of the receiving point rectangular unit according to the ray group velocity of the receiving point rectangular unit and the length of each ray.

[0050] Therefore, according to the ray length experienced in the current grid (the rectangular unit corresponding to the receiving point) Divided by the group velocity of the current grid, which can be calculated according to formula (1), the travel time of all wavefront points on the boundary of the current grid can be obtained. Since the current grid takes into account the rectangular cells adjacent to the grid, the travel time of the wavefront neighbor corresponding to each wavefront point is the minimum travel time among the travel times of all wavefront points in the grid cells.

[0051] Step S105: Based on the wavefront scanning search algorithm, the path of the first arrival ray is tracked according to the travel time of the rectangular unit of the receiving point.

[0052] In this embodiment, based on the wavefront scanning search algorithm, the first arrival ray path is tracked according to the travel time of the rectangular unit of the receiving point, including: Step b01: perform wavefront scanning on the HTI rectangular area, and determine the minimum travel time of all rectangular units according to the travel time of the rectangular units at the receiving point; Step b02: Search for the first arrival ray path based on the minimum travel time of all rectangular cells.

[0053] Among them, the wavefront scanning of the HTI rectangular area is: Let the source point (i.e. Figure 2 The travel time of the excitation position in the source is zero, and the travel time of other nodes except the source point is infinite. First, select a scanning center O, and take O as the center, and the square with r as half the length as the scanning square. Now let the r variable increase by 1, gradually expand, and scan the entire HTI rectangular area. The nodes on the four boundaries of each scanning square are wavefront points, and then the minimum travel time of the wavefront neighboring points is calculated for each wavefront point, and then shrink from the boundary of the HTI rectangular area to the scanning center, and scan from the outside to the inside to calculate the travel time of the first arrival wave and the ray path on all nodes.

[0054] When searching for the first arrival ray path, since each rectangular unit has a certain wavefront point travel time and position on its edge, the ray path from any point to the source point can be found. Figure 2 ), then take the receiving point as the starting point and find the corresponding secondary source coordinates in sequence until the source point ends. Since each secondary source coordinate records the minimum travel time, the searched path is the minimum travel time path, that is, the first arrival ray path.

[0055] In this embodiment, according to the group velocity at any position and in combination with the first arrival ray path, the minimum travel time for any shot inspection pair to pass through the anisotropic velocity field and anisotropic field path can be calculated, so as to prepare data for anisotropic multi-parameter inversion.

[0056] This embodiment can improve the adaptability of the first-arrival wave ray tracing technology. For the near-surface with azimuthal anisotropy in the collection area of ​​large offset and wide azimuth data, it can more accurately describe the ray path of the first-arrival wave, and plays a key role in near-surface anisotropy modeling.

[0057] Embodiment 2 Figure 4 is a block diagram of an azimuthal anisotropic first-arrival ray tracing device provided by one embodiment of the present invention, such as Figure 4 As shown, based on the same inventive concept as in the first embodiment, this embodiment further provides an azimuthally anisotropic first-arrival ray tracing device, which is used to implement the azimuthally anisotropic first-arrival ray tracing method in the first embodiment, and the device includes: A discrete processing module, used for performing grid discrete processing on the HTI medium to obtain an HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays; A ray determination module, used to determine the group angle and length of each ray within the rectangular unit of the receiving point; The velocity determination module is used to input the preset work area data and the group angle of each ray corresponding to the receiving point into the pre-built azimuth anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameter and vertical velocity; A travel time determination module, used to determine the travel time of the receiving point rectangular unit according to the ray group velocity of the receiving point rectangular unit and the length of each ray; The ray tracing module is used to trace the first arrival ray path according to the travel time of the rectangular unit of the receiving point based on the wavefront scanning search algorithm.

[0058] Based on the same inventive concept as in the first embodiment, this embodiment further 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 azimuthal anisotropic first-arrival ray tracing method when executing the computer program.

[0059] Based on the same inventive concept as in the first embodiment, 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 azimuthal anisotropic first-arrival ray tracing method is implemented.

[0060] This embodiment can improve the adaptability of the first-arrival wave ray tracing technology. For the near-surface with azimuthal anisotropy in the collection area of ​​large offset and wide azimuth data, it can more accurately describe the ray path of the first-arrival wave, and plays a key role in near-surface anisotropy modeling.

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

[0062] 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 device that provides the functions specified in a block or multiple blocks.

[0063] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

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

[0065] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0066] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0067] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0068] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0069] 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 method for tracing azimuthally anisotropic first-arrival rays, characterized in that: The method comprises: Performing grid discretization processing on the HTI medium to obtain an HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays; Determine the group angle and length of each ray within the rectangular unit of the receiving point; The preset work area data and the group angle of each ray corresponding to the receiving point are input into the pre-constructed azimuthal anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameter and vertical velocity; Determine the travel time of the rectangular unit at the receiving point according to the ray group velocity of the rectangular unit at the receiving point and the length of each ray; Based on the wavefront scanning search algorithm, the first arrival ray path is traced according to the travel time of the rectangular unit of the receiving point.

2. The method according to claim 1, characterized in that Based on the wavefront scanning search algorithm, the first arrival ray path is tracked according to the travel time of the rectangular unit of the receiving point, including: Perform wavefront scanning on the HTI rectangular area and determine the minimum travel time of all rectangular units based on the travel time of the rectangular units at the receiving point; Search for the first arrival ray path based on the minimum travel time of all rectangular cells.

3. The method according to claim 1, characterized in that The method further comprises: constructing an initial azimuthal anisotropic velocity model, comprising: Determine the anisotropic Thomsen parameters of VTI media; Performing coordinate axis conversion on the HTI medium and the VTI medium to obtain conversion results, wherein the conversion results include: anisotropic Thomsen parameters of the HTI medium and a conversion relationship between the HTI medium and the VTI medium; The initial model of azimuthal anisotropic velocity is constructed based on the conversion results.

4. The method according to claim 3, characterized in that The initial model of azimuthal anisotropic velocity for: ； In the formula, is the vertical speed, is the conversion relationship between HTI medium and VTI medium, is the first anisotropy parameter of the VTI medium, which is used to characterize the longitudinal anisotropy strength. is the second anisotropy parameter of VTI medium.

5. The method according to claim 4, characterized in that The conversion relationship between the HTI medium and the VTI medium is: ； In the formula, The axis of symmetry is axial, is the azimuth of the HTI medium, is the angle between the incident ray and the vertical.

6. The method according to claim 4, characterized in that The anisotropic Thomsen parameters of the HTI medium include: a first anisotropic parameter of the HTI medium, and the calculation expression of the first anisotropic parameter of the HTI medium is: ； In the formula, is the first anisotropic parameter of HTI medium, is the first anisotropy parameter of the VTI medium.

7. The method according to claim 6, characterized in that The anisotropic Thomsen parameter of the HTI medium also includes: a second anisotropic parameter of the HTI medium, and the calculation expression of the second anisotropic parameter of the HTI medium is: ； In the formula, is the second anisotropy parameter of HTI medium, is the second anisotropy parameter of the VTI medium, f is the frequency; in, , where is the shear wave velocity.

8. An azimuthally anisotropic first-arrival ray tracing device, used to implement the azimuthally anisotropic first-arrival ray tracing method according to any one of claims 1 to 7, characterized in that: The device comprises: A discrete processing module, used for performing grid discrete processing on the HTI medium to obtain an HTI rectangular area, wherein the HTI rectangular area includes a plurality of rectangular units, and each rectangular unit includes a plurality of rays; A ray determination module, used to determine the group angle and length of each ray within the rectangular unit of the receiving point; The velocity determination module is used to input the preset work area data and the group angle of each ray corresponding to the receiving point into the pre-built azimuth anisotropic velocity initial model to obtain the ray group velocity of the rectangular unit of the receiving point; wherein the preset work area data includes: HTI medium azimuth angle, symmetry axis axial direction, HTI medium anisotropy parameter and vertical velocity; A travel time determination module, used to determine the travel time of the receiving point rectangular unit according to the ray group velocity of the receiving point rectangular unit and the length of each ray; The ray tracing module is used to trace the first arrival ray path according to the travel time of the rectangular unit of the receiving point based on the wavefront scanning search algorithm.

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 azimuthal anisotropic first-arrival ray tracing method described in 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 azimuthal anisotropic first-arrival ray tracing method described in any one of claims 1 to 7 is implemented.

Images