Ray tracing method, device, apparatus and computer readable storage medium
By scaling the initial meshed velocity model at multiple scales, the problem of blind spots in ray tracing in complex media is solved by the segmented iterative method, resulting in more accurate ray tracing results and more efficient data processing.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- ENN (TIANJIN) ENERGY TECH CO LTD
- Filing Date
- 2022-09-21
- Publication Date
- 2026-06-05
AI Technical Summary
In regions with drastic changes in the medium, existing technologies using iterative ray tracing methods have blind spots, resulting in low ray tracing accuracy and significant discrepancies between the results and the actual propagation path.
By scaling the initial meshed velocity model at multiple scales, including thinning and densification, and adjusting the mesh structure, the ray paths in distant geological structures can be accurately tracked, eliminating the tracking blind spots in the conventional piecewise iterative method.
It improves the accuracy of ray tracing, making the results closer to reality, reduces the amount of computation in the data processing process, and improves the efficiency of the method.
Smart Images
Figure CN117784223B_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of data processing technology, and in particular to a ray tracing method, apparatus, device, and computer-readable storage medium. Background Technology
[0002] Ray tracing is a method that uses the differences in elasticity and density of underground media to calculate the propagation path of artificially simulated seismic waves. It plays a key role in the design of seismic exploration observation systems, seismic exploration migration imaging, microseismic location, tomography, and other processing and inversion processes.
[0003] Currently, the step-by-step iterative method is one of the most important methods for ray tracing. This method first requires determining the positions of the source and receiver points within a predetermined region, using the line connecting these two points as the initial path, and then iteratively updating the initial path point by point according to Snell's Law. However, in regions with drastic changes in the medium, this method cannot account for media far from the initial path, resulting in blind spots in ray tracing within complex media. This leads to lower accuracy in ray tracing, and the obtained results differ significantly from the actual propagation path. Summary of the Invention
[0004] To address the aforementioned technical problems, this disclosure provides a ray tracing method, apparatus, device, and computer-readable storage medium to improve the accuracy of ray tracing.
[0005] In a first aspect, embodiments of this disclosure provide a ray tracing method, including:
[0006] Determine the initial ray in the initial meshed velocity model;
[0007] The initial meshed velocity model is updated according to a preset scaling ratio to determine the path of the initial ray in the updated initial meshed velocity model.
[0008] In some embodiments, updating the initial meshed velocity model according to a preset scaling ratio and determining the ray path of the initial ray in the updated initial meshed velocity model includes:
[0009] The initial meshed velocity model is thinned to obtain a first meshed velocity model, wherein the grid line spacing in the first meshed velocity model is greater than the grid line spacing in the initial meshed velocity model.
[0010] The position of each intersection point between the initial ray path and the grid line in the first meshed velocity model is updated to obtain multiple first ray path points, and the first ray path is determined.
[0011] The first meshed velocity model is densified to obtain a second meshed velocity model, wherein the grid line spacing in the second meshed velocity model is smaller than that in the first meshed velocity model.
[0012] The position of each intersection point between the first ray path and the grid line in the second meshed velocity model is updated to obtain multiple second ray path points, and the second ray path is determined.
[0013] In some embodiments, the initial meshed velocity model is thinned to obtain a first meshed velocity model, including:
[0014] Determine the number of intersection points between the initial ray path and the grid lines in the initial meshed velocity model within the preset area;
[0015] If the number of intersection points is greater than a preset threshold, the initial meshed velocity model is thinned to obtain a first meshed velocity model.
[0016] In some embodiments, after determining the second ray path, the method further includes:
[0017] The second ray path is rotated in a preset direction by a preset angle to obtain the rotated second ray path;
[0018] The rotated second ray path is covered by a third meshed velocity model, wherein the mesh line spacing in the third meshed velocity model is the same as the mesh line spacing in the initial meshed velocity model;
[0019] The position of each intersection point between the rotated second ray path and the grid line in the third meshed velocity model is updated to obtain multiple third ray path points, and the third ray path is determined.
[0020] The third ray path is rotated by a preset angle in the opposite direction to the preset direction to obtain the fourth ray path.
[0021] In some embodiments, after obtaining the fourth ray path, the method further includes:
[0022] If the grid line spacing in the second meshed velocity model is smaller than the grid line spacing in the initial meshed velocity model, then the second meshed velocity model is densified.
[0023] The rotated fourth ray path is updated based on the encrypted second mesh velocity model and the third mesh velocity model.
[0024] Secondly, embodiments of this disclosure provide a ray tracing device, comprising:
[0025] The first determining module is used to determine the initial rays in the initial meshing velocity model;
[0026] The second determining module is used to update the initial meshed velocity model according to a preset scaling ratio and determine the path of the initial ray in the updated initial meshed velocity model.
[0027] Thirdly, embodiments of this disclosure provide an electronic device, including:
[0028] Memory;
[0029] Processor; and
[0030] Computer programs;
[0031] The computer program is stored in the memory and configured to be executed by the processor to implement the method as described in the first aspect.
[0032] Fourthly, embodiments of this disclosure provide a computer-readable storage medium having a computer program stored thereon, the computer program being executed by a processor to implement the method described in the first aspect.
[0033] Fifthly, embodiments of this disclosure also provide a computer program product, which includes a computer program or instructions that, when executed by a processor, implement the ray tracing method as described above.
[0034] The ray tracing method, apparatus, device, and computer-readable storage medium provided in this disclosure, by scaling the initial meshed velocity model at multiple scales, transforms mesh structures that are far from the initial ray path into mesh structures that are adjacent to the initial ray path, thereby enabling ray tracing in distant geological structures. This eliminates the tracing blind spots in conventional step-by-step iterative methods, making the ray tracing results closer to reality and improving the accuracy of ray tracing. Attached Figure Description
[0035] The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure.
[0036] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the accompanying drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0037] Figure 1 A flowchart of the ray tracing method provided in this embodiment of the disclosure;
[0038] Figure 2 A schematic diagram of an initial meshing velocity model provided in an embodiment of this disclosure;
[0039] Figure 3 A flowchart of a ray tracing method provided in another embodiment of this disclosure;
[0040] Figure 4 This is a schematic diagram of a first meshed velocity model provided in an embodiment of the present disclosure;
[0041] Figure 5 This is a schematic diagram of the first ray path provided in an embodiment of the present disclosure;
[0042] Figure 6 This is a schematic diagram of the second meshed velocity model provided in an embodiment of this disclosure;
[0043] Figure 7 This is a schematic diagram of the second ray path provided in an embodiment of the present disclosure;
[0044] Figure 8 A schematic diagram of the first forward model provided in this embodiment of the disclosure;
[0045] Figure 9 A schematic diagram of ray tracing results for the first forward model provided in this embodiment of the disclosure;
[0046] Figure 10 A flowchart of a ray tracing method provided in another embodiment of this disclosure;
[0047] Figure 11 This is a schematic diagram of the rotated second ray path provided in an embodiment of the present disclosure;
[0048] Figure 12 A schematic diagram of the second ray path after the third meshed velocity model is covered, as provided in an embodiment of this disclosure;
[0049] Figure 13 A schematic diagram of the third ray path provided in the embodiments of this disclosure;
[0050] Figure 14 A schematic diagram of the fourth ray path provided in the embodiments of this disclosure;
[0051] Figure 15 This is a schematic diagram of the second forward model provided in an embodiment of the present disclosure;
[0052] Figure 16 A schematic diagram of the ray tracing results of the second forward model provided in this embodiment of the disclosure;
[0053] Figure 17 This is a schematic diagram of the structure of the ray tracing device provided in the embodiments of this disclosure;
[0054] Figure 18 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this disclosure. Detailed Implementation
[0055] To better understand the above-mentioned objectives, features, and advantages of this disclosure, the solutions disclosed herein will be further described below. It should be noted that, unless otherwise specified, the embodiments and features described herein can be combined with each other.
[0056] Numerous specific details are set forth in the following description in order to provide a full understanding of this disclosure, but this disclosure may also be implemented in other ways different from those described herein; obviously, the embodiments in the specification are only some, and not all, of the embodiments of this disclosure.
[0057] This disclosure provides a ray tracing method, which will be described below with reference to specific embodiments.
[0058] Figure 1 A flowchart illustrating a ray tracing method provided in an embodiment of this disclosure. This method can be applied to... Figure 2 In the initial meshed velocity model shown. For example... Figure 2 As shown, an initial meshed velocity model is established in an electronic device with data processing capabilities. This initial meshed velocity model consists of rectangular grids of equal size, which can be viewed as a grid division of the Earth's medium space. Seismic wave propagation velocities and other attributes are defined within the grid or at the nodes. The propagation of seismic waves within the initial meshed velocity model can be simulated using a piecewise iterative method, thereby completing ray tracing. The electronic device with data processing capabilities can be a smartphone, PDA, tablet, wearable device with a display screen, desktop computer, laptop computer, all-in-one computer, etc. It is understood that the ray tracing method provided in this embodiment can also be applied to other scenarios.
[0059] The following is combined Figure 2 The initial meshed velocity model shown is for Figure 1 The ray tracing method shown is described below, and the specific steps of this method are as follows:
[0060] S101. Determine the initial ray in the initial meshed velocity model.
[0061] Seismic surveying is a geophysical exploration method that infers the properties and morphology of underground rock strata by observing and analyzing the propagation patterns of seismic waves generated by artificial earthquakes in the underground medium based on the differences in elasticity and density of the underground medium. Seismic waves are vibrations that propagate from the earthquake source in all directions; they are elastic waves that radiate outward from the source. Ray tracing is a method for calculating the propagation path of seismic waves.
[0062] The gridded velocity model is a parametric model that defines the seismic wave propagation velocity within or at nodes of multiple equally sized rectangular grids, providing a model basis for the piecewise iterative ray tracing algorithm. In, for example... Figure 2 In the initial gridded velocity model shown, given the propagation velocity of the seismic wave at each point within the grid, P1 is determined as the source point and P16 as the receiver point. Connecting point P1 and point P16 yields the initial ray path.
[0063] S102. Update the initial meshed velocity model according to a preset scaling ratio, and determine the path of the initial ray in the updated initial meshed velocity model.
[0064] After determining the initial ray in the initial meshed velocity model, the positions of each path point on the initial ray path can be corrected by a piecewise iterative method, thereby completing ray tracing. For example, the intersection of the initial ray path between the source point P1 and the receiver point P16 and the grid lines in the initial meshed velocity model is the path points P2 to P15 of the initial ray.
[0065] Since the segmented iterative method updates the ray path between every two path points based on the seismic wave propagation velocity of the medium corresponding to adjacent grids, it cannot account for the influence of the medium corresponding to grids far from the initial ray path on ray tracing in the initially gridded velocity model with a relatively dense grid. Therefore, the initial gridded velocity model is updated according to a preset scaling ratio. Specifically, the initial gridded velocity model is first enlarged according to a preset ratio, that is, the grid density in the initial gridded velocity model is reduced, so that the medium corresponding to the grids far from the initial ray path in the initial gridded velocity model becomes the medium corresponding to the grid where the initial ray path is located or the grid adjacent to the grid where the initial ray path is located in the initial gridded velocity model with reduced grid density. Ray tracing of the initial ray is then performed in the initial gridded velocity model with reduced grid density to obtain a new ray path. Then, the initial gridded velocity model with reduced grid density is scaled down according to a preset ratio, that is, its grid density is restored to the same as the initial gridded velocity model. Ray tracing of the new ray is performed again in the restored initial gridded velocity model, and finally the path of the initial ray in the updated initial gridded velocity model is determined. Optionally, the mesh density in the initial meshed velocity model after the mesh density reduction can be increased successively. After each increase in mesh density, ray tracing is performed iteratively until the mesh density is restored to the same as the initial meshed velocity model, and the path of the initial ray in the updated initial meshed velocity model is determined.
[0066] This embodiment of the disclosure determines an initial ray in an initial meshed velocity model; updates the initial meshed velocity model according to a preset scaling ratio, determines the path of the initial ray in the updated initial meshed velocity model, and by scaling the initial meshed velocity model at multiple scales, mesh structures far from the initial ray path become mesh structures adjacent to the initial ray path, thereby realizing ray tracing in distant geological structures, eliminating the tracking blind spots in conventional piecewise iterative methods, making the ray tracing results closer to the real situation, and improving the accuracy of ray tracing.
[0067] Figure 3 A flowchart of a ray tracing method provided in another embodiment of this disclosure is shown below. Figure 3 As shown, the method includes the following steps:
[0068] S301. Determine the initial ray in the initial meshed velocity model.
[0069] S302. The initial meshed velocity model is thinned to obtain a first meshed velocity model, wherein the grid line spacing in the first meshed velocity model is greater than the grid line spacing in the initial meshed velocity model.
[0070] Figure 4 This is a schematic diagram of a first meshed velocity model provided in an embodiment of this disclosure. Figure 4 As shown, it will be as follows Figure 2 The initial meshed velocity model shown is thinned out, i.e., its scale is enlarged from the original 9*9 initial meshed velocity model to a 3*3 first meshed velocity model. Points A1 and A4 correspond to the source point P1 and receiver point P16 in the initial meshed velocity model, respectively. The line connecting points A1 and A4 is the initial ray path. The intersections A2 and A3 of the line connecting points A1 and A4 with the grid lines in the first meshed velocity model are the path points of the initial ray path in the first meshed velocity model. Points A2 and A3 correspond to path points P8 and P9 in the initial meshed velocity model, respectively.
[0071] S303. Update the position of each intersection point between the initial ray path and the grid line in the first meshed velocity model to obtain multiple first ray path points and determine the first ray path.
[0072] Figure 5 This is a schematic diagram of the first ray path provided in an embodiment of this disclosure. Using a segment-by-segment iterative method, the positions of A2 and A3 are updated based on the seismic wave propagation velocities on both sides of the grid line where the path point of the initial ray path is located in the first gridded velocity model, resulting in the following... Figure 5 The first ray path points B2 and B3 shown are further determined as follows: Figure 5The first ray path is shown. Points B1 and B4 correspond to the source point P1 and receiver point P16 in the initial meshed velocity model, respectively.
[0073] S304. The first meshed velocity model is densified to obtain a second meshed velocity model, wherein the grid line spacing in the second meshed velocity model is smaller than the grid line spacing in the first meshed velocity model.
[0074] Figure 6 This is a schematic diagram of the second meshed velocity model provided in an embodiment of this disclosure. Figure 6 As shown, it will be as follows Figure 4 The first gridded velocity model shown is refined by reducing its scale from a 3x3 model to a 5x5 model. Points C1 and C8 correspond to the source point P1 and receiver point P16 in the initial gridded velocity model, respectively. The line connecting points C1 and C8 forms the initial ray path. The intersections C2 to C7 of the first ray path between points C1 and C8 with the grid lines in the second gridded velocity model are the path points of the first ray path in the second gridded velocity model.
[0075] S305. Update the position of each intersection point of the first ray path and the grid line in the second meshed velocity model to obtain multiple second ray path points and determine the second ray path.
[0076] Figure 7 This is a schematic diagram of the second ray path provided in an embodiment of this disclosure. Using a segment-by-segment iterative method, the positions of points C2 to C7 are updated based on the seismic wave propagation velocities on both sides of the grid lines where the path points of the first ray path are located in the second gridded velocity model, resulting in the following... Figure 7 The second ray path points D2 to D7 shown are further determined as follows: Figure 7 The second ray path is shown. Points D1 and D8 correspond to the source point P1 and receiver point P16 in the initial meshed velocity model, respectively.
[0077] Repeat the above steps to refine the meshed velocity model and iteratively update the ray paths within it until the mesh line spacing of the refined meshed velocity model is the same as that of the initial meshed velocity model. Then determine the path of the initial ray in the updated initial meshed velocity model.
[0078] Figure 8This is a schematic diagram of a first forward model provided in an embodiment of this disclosure. Forward modeling is (in geophysical exploration research, based on the shape, occurrence, and physical properties of a geological body, calculating its theoretical values by constructing a mathematical model (mathematical simulation), or observing the numerical values of the geophysical effects produced by the model by constructing a physical model (physical simulation). The model required for forward modeling is called a forward model. For example... Figure 8 As shown, the first forward model provided in this embodiment is a two-dimensional model with a depth of 5000 meters and a length of 8000 meters. The subsurface medium interface is located at a depth of 2000 meters. The propagation velocity of seismic waves above a depth of 2000 meters (shallow layer) is 2000 m / s, and the propagation velocity of seismic waves below a depth of 2000 meters (deep layer) is 5000 m / s. The source point S1 (2000, 1500) and the receiver point R1 (6000, 1500) are set, and ray tracing is performed in the first forward model according to the method described in the above embodiment. Figure 9 A schematic diagram of ray tracing results for a first forward model provided in an embodiment of this disclosure. (See diagram below.) Figure 9 As shown, the dashed line represents the ray tracing result obtained by the conventional step-by-step iterative method, and the solid line represents the ray tracing result obtained by the ray tracing method provided in this embodiment. It can be seen that the ray tracing result obtained by the ray tracing method provided in this embodiment avoids the perception blind zone generated near the interface of the underground medium, and the ray tracing result obtained is closer to the path in the real situation than the ray tracing result obtained by the conventional step-by-step iterative method.
[0079] This embodiment of the disclosure performs multi-scale thinning and densification of the initial meshed velocity model, transforming mesh structures far from the initial ray path into those adjacent to it. This enables ray tracing in distant geological structures, eliminating the tracing blind spots inherent in conventional piecewise iterative methods. The ray tracing results are thus closer to reality, improving accuracy. Furthermore, since the ray path points in the thinned meshed velocity model correspond to the path points of the initial ray path, the new path points obtained through interpolation calculations during subsequent densification of the meshed velocity model retain the attribute values from the original model. This improves accuracy while reducing computational load during data processing, thereby increasing efficiency.
[0080] Figure 10 A flowchart of a ray tracing method provided in another embodiment of this disclosure is shown below. Figure 10 As shown, the method includes the following steps.
[0081] S1001. Determine the initial ray in the initial meshed velocity model.
[0082] S1002. Determine the number of intersection points between the initial ray path and the grid lines in the initial meshed velocity model within the preset area.
[0083] S1003. If the number of intersection points is greater than a preset threshold, the initial meshed velocity model is thinned to obtain a first meshed velocity model.
[0084] The preset region can be a homogeneous area in the Earth's medium where seismic wave propagation speed is uniform. Thinning the meshed velocity model to a homogeneous region with only ray path points less than or equal to a preset threshold ensures the accuracy of path point locations after segment-by-segment iterative calculations to the greatest extent possible.
[0085] S1004. Update the position of each intersection point between the initial ray path and the grid line in the first meshed velocity model to obtain multiple first ray path points and determine the first ray path.
[0086] S1005. The first gridded velocity model is encrypted to obtain the second gridded velocity model.
[0087] S1006. Update the position of each intersection point of the first ray path and the grid line in the second meshed velocity model to obtain multiple second ray path points and determine the second ray path.
[0088] Specifically, the implementation principles and specific implementation processes of S1004 to S1006 are the same, and will not be repeated here.
[0089] S1007. Rotate the second ray path in a preset direction by a preset angle to obtain the rotated second ray path.
[0090] Figure 11 This is a schematic diagram of the rotated second ray path provided in an embodiment of this disclosure. After obtaining the second ray path in the second meshed velocity model, the second ray path is rotated by a preset angle to obtain, as shown below. Figure 11 The diagram shows the rotated second ray path and the corresponding second gridded velocity model. In actual seismic wave propagation, folded waves may occur, where the corresponding ray path "folds back" within a certain grid. That is, after entering the grid from one grid line, it leaves the grid again from the same grid line and enters the next adjacent grid. The piecewise iterative method, based on the law of refraction, cannot solve the problem of waves entering from an interface and then exiting from the same interface. Therefore, the piecewise iterative method cannot track folded waves.
[0091] like Figure 11As shown, the second ray path is rotated 45 degrees clockwise to obtain the rotated second ray path and the corresponding second meshed velocity model. It is understood that the rotation angle and direction provided in this embodiment are merely examples and not limitations of this disclosure. In practice, adjustments can be made based on the size, parameters, and complexity of the meshed velocity model.
[0092] S1008. The rotated second ray path is covered by a third meshed velocity model, wherein the mesh line spacing in the third meshed velocity model is the same as the mesh line spacing in the initial meshed velocity model.
[0093] Figure 12 This is a schematic diagram of the second ray path after being covered by the third meshed velocity model, as provided in an embodiment of this disclosure. Figure 12 As shown, a third gridded velocity model is used to cover the second ray path and its corresponding second gridded velocity model. The grid lines in the third gridded velocity model should be horizontal or vertical, and the grid line spacing should be the same as that in the initial gridded velocity model. Points E1 and E17 correspond to the source point P1 and receiver point P16 in the initial gridded velocity model, respectively. The intersection points E2 to E16 of the rotated second ray path and the grid lines in the third gridded velocity model are the path points of the second ray path in the third gridded velocity model.
[0094] S1009. Update the position of each intersection point between the rotated second ray path and the grid line in the third meshed velocity model to obtain multiple third ray path points and determine the third ray path.
[0095] Figure 13 This is a schematic diagram of the third ray path provided in an embodiment of this disclosure. Using a segment-by-segment iterative method, the positions of E2 to E16 are updated based on the seismic wave propagation velocities on both sides of the grid lines where the path points of the second ray path are located in the third gridded velocity model, resulting in the following... Figure 13 As shown in F2~F16, points F1 and F17 correspond to the source point P1 and receiver point P16 in the initial meshed velocity model, respectively. Further determination is needed as follows... Figure 13 The third ray path is shown.
[0096] S1010. Rotate the third ray path in the opposite direction of the preset direction by a preset angle to obtain the fourth ray path.
[0097] Figure 14This is a schematic diagram of the fourth ray path provided in an embodiment of this disclosure. The third ray path is rotated 45 degrees counterclockwise to obtain the fourth ray path and its corresponding second gridded velocity model, as shown in Figure 14. Points G1 and G8 correspond to the source point P1 and receiver point P16 in the initial gridded velocity model, respectively. G2 to G7 are the intersection points of the fourth ray path and the second gridded velocity model, i.e., the path points of the fourth ray path in the second gridded velocity model.
[0098] Repeat the above steps to refine the meshed velocity model, iteratively updating the ray paths until the mesh line spacing of the refined meshed velocity model is the same as that of the initial meshed velocity model. Specifically, after each refinement operation on the meshed velocity model, the ray paths are iterated segment by segment and then rotated. After another segment-by-segment iteration, they are rotated in the opposite direction to the original angle, and the next refinement operation is performed.
[0099] Figure 15 This is a schematic diagram of the second forward model provided in this embodiment. The second forward model provided in this embodiment is a two-dimensional model with a depth of 5000 meters and a length of 8000 meters. The seismic wave propagation velocity increases linearly along the direction of the arrow. The velocity value V1 at (0,0) is 2000 m / s, and the velocity value V2 at (5000,5000) is 5000 m / s. The source point S2 (2000,500) and the receiver point R2 (6000,500) are set. Ray tracing is performed in the second forward model according to the method described in the above embodiment. Figure 16 This is a schematic diagram of the ray tracing results of the second forward model provided in an embodiment of this disclosure. Figure 16 As shown, the solid curve represents the actual ray path, the dotted line represents the ray tracing result obtained by the conventional step-by-step iterative method, and the dashed curve represents the ray tracing result obtained by the ray tracing method provided in this embodiment. It can be seen that the ray tracing result obtained by the ray tracing method provided in this embodiment is closer to the actual ray path.
[0100] This embodiment of the present disclosure, by performing multi-scale thinning and densification of the initial gridded velocity model to achieve ray tracing in distant geological structures, obtains a multi-directional gridded velocity model, changes the direction of the grid lines, and adjusts rays such as reflection waves that do not satisfy the law of refraction under the original grid direction to a state that can be updated segment by segment iteratively. This enables ray tracing of reflection waves, etc., further reduces blind spots in the ray tracing process, improves the ray tracing effect in complex geological structures, and significantly improves the data processing effect of seismic exploration and other work in complex areas.
[0101] Figure 17This is a schematic diagram of the structure of a ray tracing device provided in an embodiment of this disclosure. The ray tracing device can be an electronic device with data processing capabilities as described in the above embodiments, or it can be a component or assembly within that electronic device. The ray tracing device provided in this disclosure can execute the processing flow provided in the ray tracing method embodiments, such as... Figure 17 As shown, the ray tracing device 170 includes: a first determining module 171 and a second determining module 172; wherein, the first determining module 171 is used to determine the initial ray in the initial meshed velocity model; the second determining module 172 is used to update the initial meshed velocity model according to a preset scaling ratio, and determine the path of the initial ray in the updated initial meshed velocity model.
[0102] Optionally, the second determining module 172 includes a thinning unit 1721, a path determining unit 1722, and an encryption unit 1723; wherein, the thinning unit 1721 is used to thin the initial meshed velocity model to obtain a first meshed velocity model, wherein the grid line spacing in the first meshed velocity model is greater than the grid line spacing in the initial meshed velocity model; the path determining unit 1722 is used to update the position of each intersection point of the initial ray path and the grid lines in the first meshed velocity model to obtain multiple first ray path points and determine the first ray path; the encryption unit 1723 is used to encrypt the first meshed velocity model to obtain a second meshed velocity model, wherein the grid line spacing in the second meshed velocity model is less than the grid line spacing in the first meshed velocity model; the path determining unit 1722 is also used to update the position of each intersection point of the first ray path and the grid lines in the second meshed velocity model to obtain multiple second ray path points and determine the second ray path.
[0103] Optionally, the thinning unit 1721 is further configured to determine the number of intersections between the initial ray path and the grid lines in the initial meshed velocity model within a preset area; if the number of intersections is greater than a preset threshold, the initial meshed velocity model is thinned to obtain a first meshed velocity model.
[0104] Optionally, the second determining module 172 further includes a rotation unit 1724, used to rotate the second ray path in a preset direction by a preset angle to obtain a rotated second ray path; cover the rotated second ray path with a third meshed velocity model, wherein the grid line spacing in the third meshed velocity model is the same as the grid line spacing in the initial meshed velocity model; the path determining unit 1722 is also used to update the position of each intersection point of the rotated second ray path and the grid lines in the third meshed velocity model to obtain multiple third ray path points and determine the third ray path; the rotation unit 1724 is also used to rotate the third ray path in the opposite direction of the preset direction by a preset angle to obtain a fourth ray path.
[0105] Optionally, the encryption unit 1723 is further configured to encrypt the second meshed velocity model when the mesh line spacing in the second meshed velocity model is less than the mesh line spacing in the initial meshed velocity model; the path determination unit 1722 is further configured to update the rotated fourth ray path according to the encrypted second meshed velocity model and the third meshed velocity model.
[0106] Figure 17 The ray tracing device shown in the embodiment can be used to execute the technical solution of the above-described ray tracing method embodiment. Its implementation principle and technical effect are similar, and will not be described again here.
[0107] Figure 18 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this disclosure. The electronic device can be an electronic device with data processing capabilities as described in the above embodiments. The electronic device provided in this disclosure can execute the processing flow provided in the ray tracing method embodiments, such as… Figure 18 As shown, the electronic device 180 includes: a memory 181, a processor 182, a computer program, and a communication interface 183; wherein the computer program is stored in the memory 181 and is configured to be executed by the processor 182 using the ray tracing method described above.
[0108] In addition, this disclosure also provides a computer-readable storage medium having a computer program stored thereon, the computer program being executed by a processor to implement the ray tracing method described in the above embodiments.
[0109] Furthermore, this disclosure also provides a computer program product, which includes a computer program or instructions that, when executed by a processor, implement the ray tracing method described above.
[0110] Furthermore, this disclosure also provides a computer program product, which includes a computer program or instructions that, when executed by a processor, implement the vehicle voice control method described above.
[0111] Computer program code for performing the operations of this disclosure can be written in one or more programming languages or a combination thereof, including but not limited to object-oriented programming languages such as Java, Smalltalk, and C++, as well as conventional procedural programming languages such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).
[0112] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.
[0113] The above description is merely a specific embodiment of this disclosure, enabling those skilled in the art to understand or implement it. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of this disclosure. Therefore, this disclosure is not to be limited to the embodiments described herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims
1. A ray tracing method, characterized in that, The method includes: Determine the initial ray in the initial meshed velocity model; The initial meshed velocity model is thinned to obtain a first meshed velocity model, wherein the grid line spacing in the first meshed velocity model is greater than the grid line spacing in the initial meshed velocity model. The position of each intersection point between the initial ray path and the grid line in the first meshed velocity model is updated to obtain multiple first ray path points, and the first ray path is determined. The first meshed velocity model is densified to obtain a second meshed velocity model, wherein the grid line spacing in the second meshed velocity model is smaller than that in the first meshed velocity model. The position of each intersection point between the first ray path and the grid line in the second meshed velocity model is updated to obtain multiple second ray path points, and the second ray path is determined. The second ray path is rotated in a preset direction by a preset angle to obtain the rotated second ray path; The rotated second ray path is covered by a third meshed velocity model, wherein the mesh line spacing in the third meshed velocity model is the same as the mesh line spacing in the initial meshed velocity model; The position of each intersection point between the rotated second ray path and the grid line in the third meshed velocity model is updated to obtain multiple third ray path points, and the third ray path is determined. The third ray path is rotated by a preset angle in the opposite direction to the preset direction to obtain the fourth ray path.
2. The method according to claim 1, characterized in that, The initial meshed velocity model is thinned to obtain a first meshed velocity model, including: Determine the number of intersection points between the initial ray path and the grid lines in the initial meshed velocity model within the preset area; If the number of intersection points is greater than a preset threshold, the initial meshed velocity model is thinned to obtain a first meshed velocity model.
3. The method according to claim 1, characterized in that, After obtaining the fourth ray path, the method further includes: If the grid line spacing in the second meshed velocity model is smaller than the grid line spacing in the initial meshed velocity model, then the second meshed velocity model is densified. The rotated fourth ray path is updated based on the encrypted second mesh velocity model and the third mesh velocity model.
4. A ray tracing device, characterized in that, include: The first determining module is used to determine the initial rays in the initial meshing velocity model; The second determining module is used to thin out the initial meshed velocity model to obtain a first meshed velocity model, wherein the grid line spacing in the first meshed velocity model is greater than the grid line spacing in the initial meshed velocity model. The position of each intersection point between the initial ray path and the grid lines in the first meshed velocity model is updated to obtain multiple first ray path points, and the first ray path is determined; the first meshed velocity model is densified to obtain a second meshed velocity model, wherein the grid line spacing in the second meshed velocity model is smaller than that in the first meshed velocity model; The position of each intersection point between the first ray path and the grid line in the second meshed velocity model is updated to obtain multiple second ray path points, and the second ray path is determined. The second determining module is also used for: The second ray path is rotated in a preset direction by a preset angle to obtain the rotated second ray path; The rotated second ray path is covered by a third meshed velocity model, wherein the mesh line spacing in the third meshed velocity model is the same as the mesh line spacing in the initial meshed velocity model; The position of each intersection point between the rotated second ray path and the grid line in the third meshed velocity model is updated to obtain multiple third ray path points, and the third ray path is determined. The third ray path is rotated by a preset angle in the opposite direction to the preset direction to obtain the fourth ray path.
5. An electronic device, characterized in that, include: Memory; processor; as well as Computer programs; The computer program is stored in the memory and configured to be executed by the processor to implement the method as described in any one of claims 1-3.
6. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the method as described in any one of claims 1-3.