CN117784223A - Ray tracing method, device, equipment and computer readable storage medium - Google Patents
Ray tracing method, device, equipment and computer readable storage medium Download PDFInfo
- Publication number
- CN117784223A CN117784223A CN202211151458.8A CN202211151458A CN117784223A CN 117784223 A CN117784223 A CN 117784223A CN 202211151458 A CN202211151458 A CN 202211151458A CN 117784223 A CN117784223 A CN 117784223A
- Authority
- CN
- China
- Prior art keywords
- speed model
- gridding
- initial
- ray
- ray path
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000004590 computer program Methods 0.000 claims description 17
- 238000010586 diagram Methods 0.000 description 22
- 238000012545 processing Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Abstract
The present disclosure relates to a ray tracing method, apparatus, device and computer readable storage medium, the method comprising: determining an initial ray in an initial meshing speed model; updating the initial gridding speed model according to a preset scaling ratio, and determining the path of the initial ray in the updated initial gridding speed model. According to the method and the device, the initial grid speed model is scaled in multiple scales, so that the grid structure far away from the initial ray path is changed into the grid structure adjacent to the initial ray path, ray tracing in a remote geological structure is realized, a tracing blind area in a conventional piecewise iteration method is eliminated, a ray tracing result is closer to a real condition, and the accuracy of ray tracing is improved.
Description
Technical Field
The present disclosure relates to the field of data processing technologies, and in particular, to a ray tracing method, apparatus, device, and computer readable storage medium.
Background
Ray tracing is a method for calculating the propagation path of an artificially simulated seismic wave by utilizing the difference of elasticity and density of underground media, and plays a key role in the processes of seismic exploration observation system design, seismic exploration offset imaging, microseism positioning, tomography and inversion.
Currently, the piecewise iterative method is one of the most dominant methods of ray tracing. The piecewise iterative method firstly needs to determine the positions of a vibration source point and a detection point in a preset area, takes a connecting line between the two points as an initial path, and carries out iterative updating on the initial path point by point according to a refraction Law (Snell's Law). However, in the region where the medium changes severely, the method cannot consider the medium far from the initial path, so that the ray tracing in the complex medium has a dead zone, the accuracy of the ray tracing is lower, and the obtained result is greatly different from the real propagation path.
Disclosure of Invention
In order to solve the above technical problems, the present disclosure provides a ray tracing method, apparatus, device and computer readable storage medium, so as to improve the accuracy of ray tracing.
In a first aspect, an embodiment of the present disclosure provides a ray tracing method, including:
determining an initial ray in an initial meshing speed model;
updating the initial gridding speed model according to a preset scaling ratio, and determining the path of the initial ray in the updated initial gridding speed model.
In some embodiments, updating the initial gridding speed model according to a preset scaling, and determining a ray path of the initial ray in the updated initial gridding speed model includes:
performing thinning treatment on the initial gridding speed model to obtain a first gridding speed model, wherein the grid line spacing in the first gridding speed model is larger than that in the initial gridding speed model;
updating the position of each intersection point of the initial ray path and the grid lines in the first gridding speed model to obtain a plurality of first ray path points, and determining a first ray path;
encrypting the first gridding speed model to obtain a second gridding speed model, wherein the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the first gridding speed model;
and updating the position of each intersection point of the grid lines in the first ray path and the second meshing speed model to obtain a plurality of second ray path points, and determining a second ray path.
In some embodiments, performing thinning processing on the initial meshing speed model to obtain a first meshing speed model, including:
determining the number of intersection points of the initial ray paths in a preset area and grid lines in the initial gridding speed model;
and if the number of the intersection points is larger than a preset threshold value, performing thinning treatment on the initial meshing speed model to obtain a first meshing speed model.
In some embodiments, after determining the second ray path, the method further comprises:
rotating the second ray path by a preset angle towards a preset direction to obtain a rotated second ray path;
covering the rotated second ray path by a third gridding speed model, wherein the grid line spacing in the third gridding speed model is the same as the grid line spacing in the initial gridding speed model;
updating the position of each intersection point of the rotated second ray path and the grid lines in the third meshing speed model to obtain a plurality of third ray path points, and determining a third ray path;
and rotating the third ray path by a preset angle towards the reverse direction of the preset direction to obtain a fourth ray path.
In some embodiments, after deriving the fourth ray path, the method further comprises:
if the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the initial gridding speed model, encrypting the second gridding speed model;
and updating the rotated fourth ray path according to the encrypted second meshing speed model and the third meshing speed model.
In a second aspect, an embodiment of the present disclosure provides a ray tracing apparatus, including:
a first determining module for determining an initial ray in an initial gridding speed model;
and the second determining module is used for updating the initial gridding speed model according to a preset scaling ratio and determining the path of the initial ray in the updated initial gridding speed model.
In a third aspect, an embodiment of the present disclosure provides an electronic device, including:
a memory;
a processor; and
a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement the method according to the first aspect.
In a fourth aspect, embodiments of the present disclosure provide a computer-readable storage medium having stored thereon a computer program for execution by a processor to implement the method of the first aspect.
In a fifth aspect, the disclosed embodiments also provide a computer program product comprising a computer program or instructions which, when executed by a processor, implement a ray tracing method as described above.
According to the ray tracing method, the device, the equipment and the computer readable storage medium, the grid structure far away from the initial ray path is changed into the grid structure adjacent to the initial ray path by performing multi-scale scaling on the initial grid speed model, so that ray tracing in a remote geological structure is realized, a tracing blind area in a conventional piecewise iteration method is eliminated, a ray tracing result is closer to a real condition, and the ray tracing accuracy is improved.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description, serve to explain the principles of the disclosure.
In order to more clearly illustrate the embodiments of the present disclosure or the solutions in the prior art, the drawings that are required for the description of the embodiments or the prior art will be briefly described below, and it will be obvious to those skilled in the art that other drawings can be obtained from these drawings without inventive effort.
FIG. 1 is a flow chart of a ray tracing method provided in an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of an initial meshing speed model according to an embodiment of the present disclosure;
FIG. 3 is a flowchart of a ray tracing method according to another embodiment of the present disclosure;
FIG. 4 is a schematic diagram of a first meshing speed model provided by an embodiment of the present disclosure;
FIG. 5 is a first ray path schematic provided by an embodiment of the present disclosure;
FIG. 6 is a schematic diagram of a second meshing speed model provided by an embodiment of the present disclosure;
FIG. 7 is a second ray path schematic provided by an embodiment of the present disclosure;
FIG. 8 is a schematic diagram of a first forward model provided by an embodiment of the present disclosure;
FIG. 9 is a schematic diagram of a first forward model ray tracing result provided by an embodiment of the present disclosure;
FIG. 10 is a flowchart of a ray tracing method according to another embodiment of the present disclosure;
FIG. 11 is a schematic diagram of a rotated second ray path provided by an embodiment of the present disclosure;
FIG. 12 is a schematic view of a second ray path covered by a third meshing speed model provided by an embodiment of the present disclosure;
FIG. 13 is a third ray path schematic provided by an embodiment of the present disclosure;
FIG. 14 is a fourth ray path schematic provided by an embodiment of the present disclosure;
FIG. 15 is a schematic diagram of a second forward model provided by an embodiment of the present disclosure;
FIG. 16 is a diagram illustrating a second forward model ray tracing result provided by an embodiment of the present disclosure;
FIG. 17 is a schematic diagram of a ray tracing apparatus according to an embodiment of the disclosure;
fig. 18 is a schematic structural diagram of an electronic device according to an embodiment of the disclosure.
Detailed Description
In order that the above objects, features and advantages of the present disclosure may be more clearly understood, a further description of aspects of the present disclosure will be provided below. It should be noted that, without conflict, the embodiments of the present disclosure and features in the embodiments may be combined with each other.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present disclosure, but the present disclosure may be practiced otherwise than as described herein; it will be apparent that the embodiments in the specification are only some, but not all, embodiments of the disclosure.
Embodiments of the present disclosure provide a ray tracing method, which is described below in connection with specific embodiments.
Fig. 1 is a flowchart of a ray tracing method according to an embodiment of the disclosure. The method can be applied to the initial meshing speed model shown in fig. 2. As shown in fig. 2, an initial meshing speed model is established in the electronic device with the data processing function, and the initial meshing speed model is a rectangular grid with equal size and can be regarded as meshing the earth medium space. The propagation velocity and other attributes of the seismic waves are defined in the grid or on the nodes, and the propagation of the seismic waves in the initial grid velocity model can be simulated through a piecewise iteration method, so that ray tracing is completed. The electronic device with the data processing function can be a smart phone, a palm computer, a tablet computer, a wearable device with a display screen, a desktop computer, a notebook computer, an integrated machine and the like. It can be appreciated that the ray tracing method provided by the embodiments of the present disclosure may also be applied in other scenarios.
The ray tracing method shown in fig. 1 is described below with reference to the initial gridding velocity model shown in fig. 2, and includes the following specific steps:
s101, determining an initial ray in an initial meshing speed model.
Seismic surveying refers to a geophysical prospecting method for deducing the nature and morphology of subsurface rock formations by observing and analyzing the propagation law of seismic waves generated by artificial earthquakes in subsurface media based on the differences in elasticity and density of the subsurface media. Seismic waves refer to vibration propagated around by a seismic source, and refer to elastic waves radiated around from the source, and ray tracing is a method for calculating the propagation path of the seismic waves.
The gridding velocity model is a parameterized model that provides a model basis for a piecewise iterative ray tracing algorithm by defining respective seismic wave propagation velocities within or on a plurality of equally sized rectangular grids. In the initial gridding velocity model shown in FIG. 2, given the propagation velocity of the seismic wave at each point within the grid, P1 is determined to be the source point, P16 is the detector point, and the point P1 and point P16 are connected to obtain an initial ray path.
S102, updating the initial meshing speed model according to a preset scaling ratio, and determining a path of the initial ray in the updated initial meshing speed model.
After the initial rays in the initial grid speed model are determined, the positions of all path points on the initial ray path can be corrected through a piecewise iterative method, so that ray tracing is completed. For example, the intersections of the initial ray paths between the source point P1 and the detector point P16 and the grid lines in the initial meshing speed model are the path points P2 to P15 of the initial rays.
Because the piecewise iteration method is to carry out iterative update on ray paths between every two path points based on the seismic wave propagation speeds of media corresponding to adjacent grids, the influence of the media corresponding to grids at a position far from the initial ray paths on ray tracking cannot be considered in an initial gridding speed model with denser grids. Thus, the initial meshing speed model is updated according to the preset scaling. The method comprises the following steps: the initial gridding speed model is amplified according to a preset proportion, namely the grid density in the initial gridding speed model is reduced, so that a medium corresponding to a grid far away from an initial ray path in the initial gridding speed model becomes a medium corresponding to a grid where the initial ray path is located or a grid adjacent to the grid where the initial ray path is located in the initial gridding speed model after the grid density is reduced, and the initial ray is tracked in the initial gridding speed model after the grid density is reduced, so that a new ray path is obtained. And then reducing the initial gridding speed model with reduced grid density according to a preset proportion, namely restoring the grid density to be the same as that of the initial gridding speed model, carrying out ray tracing on new rays in the restored initial gridding speed model again, and finally determining the path of the initial rays in the updated initial gridding speed model. Optionally, the density of the grids in the initial gridding speed model after the grid density is reduced can be gradually increased, and the ray tracing of the iteration segment by segment is performed after the grid density is increased each time until the grid density in the initial gridding speed model is restored to be the same as the initial gridding speed model, and the path of the initial ray in the updated initial gridding speed model is determined.
Embodiments of the present disclosure initialize a velocity model by determining an initial ray in the initial meshing velocity model; updating the initial gridding speed model according to a preset scaling ratio, determining the path of the initial ray in the updated initial gridding speed model, and scaling the initial gridding speed model in a multi-scale mode to change a grid structure far away from the initial ray path into a grid structure adjacent to the initial ray path, so that ray tracing in a remote geological structure is realized, a tracing blind area in a conventional piecewise iteration method is eliminated, a ray tracing result is more approximate to a real condition, and the accuracy of ray tracing is improved.
Fig. 3 is a flowchart of a ray tracing method according to another embodiment of the disclosure, as shown in fig. 3, the method includes the following steps:
s301, determining an initial ray in an initial meshing speed model.
S302, performing thinning treatment on the initial gridding speed model to obtain a first gridding speed model, wherein the grid line spacing in the first gridding speed model is larger than that in the initial gridding speed model.
Fig. 4 is a schematic diagram of a first meshing speed model according to an embodiment of the disclosure. As shown in fig. 4, the initial gridding velocity model shown in fig. 2 is subjected to thinning processing, that is, the initial gridding velocity model is scaled up, the initial gridding velocity model is changed from the initial 9*9 initial gridding velocity model into a 3*3 first gridding velocity model, wherein a point A1 and a point A4 respectively correspond to a seismic source point P1 and a detection point P16 in the initial gridding velocity model, and a connection line between the points A1 and A4 is an initial ray path. The intersection points A2 and A3 of the connecting line of the points A1 and A4 and the grid lines in the first grid speed model are the path points of the initial ray path in the first grid speed model, and the point A2 and the point A3 respectively correspond to the path points P8 and P9 in the initial grid speed model.
S303, updating the position of each intersection point of the initial ray path and the grid lines in the first meshing speed model to obtain a plurality of first ray path points, and determining a first ray path.
Fig. 5 is a schematic diagram of a first ray path provided by an embodiment of the present disclosure. And updating the positions of A2 and A3 according to the seismic wave propagation speeds of the initial ray path on two sides of the grid line where the path points in the first grid speed model are positioned by using a piecewise iteration method to obtain first ray path points B2 and B3 shown in FIG. 5, and further determining the first ray path shown in FIG. 5. Point B1 and point B4 correspond to source point P1 and detector point P16, respectively, in the initial gridding velocity model.
S304, encrypting the first gridding speed model to obtain a second gridding speed model, wherein the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the first gridding speed model.
Fig. 6 is a schematic diagram of a second meshing speed model according to an embodiment of the disclosure. As shown in fig. 6, the first gridding speed model shown in fig. 4 is encrypted, that is, the scale of the first gridding speed model is reduced, the first gridding speed model is changed from 3*3 to 5*5 to a second gridding speed model, wherein the point C1 and the point C8 correspond to the source point P1 and the detector point P16 in the initial gridding speed model respectively, and the connection line between the point C1 and the point C8 is the initial ray path. And the intersection points C2-C7 of the grid lines in the first ray path and the second grid speed model between the point C1 and the point C8 are the path points of the first ray path in the second grid speed model.
S305, updating the position of each intersection point of the grid lines in the first ray path and the second meshing speed model to obtain a plurality of second ray path points, and determining a second ray path.
Fig. 7 is a schematic diagram of a second ray path provided by an embodiment of the present disclosure. And updating the positions of the points C2-C7 according to the seismic wave propagation speeds of the first ray path on two sides of the grid line where the path points in the second grid speed model are located by using a piecewise iteration method to obtain second ray path points D2-D7 shown in FIG. 7, and further determining a second ray path shown in FIG. 7. Point D1 and point D8 correspond to source point P1 and detector point P16, respectively, in the initial gridding velocity model.
And repeatedly encrypting the gridding speed model according to the steps, and iteratively updating the ray paths in the gridding speed model until the grid line spacing of the encrypted gridding speed model is the same as that of the initial gridding speed model, and determining the paths of the initial rays in the updated initial gridding speed model.
Fig. 8 is a schematic diagram of a first forward model according to an embodiment of the disclosure. In geophysical exploration, theoretical values (mathematical simulation) of the geological body are obtained through calculation by constructing a mathematical model according to the shape, the occurrence and the physical property data of the geological body, or values (physical simulation) of geophysical effects generated by observing the model by constructing a physical model, the model required by the forward modeling is the forward model, as shown in fig. 8, the first forward model provided by the embodiment of the disclosure is a two-dimensional model with depth of 5000 meters and length of 8000 meters, the surface of a subsurface medium is at the depth of 2000 meters, the propagation speed V1 of a shallow layer is 2000m/S at the depth of more than 2000 meters, the propagation speed V2 of a deep layer is 5000m/S, a seismic source point S1 (2000,1500) and a detection point R1 (6000,1500) are set, and ray tracing is performed in the first forward model according to the method described in the embodiment of the disclosure.
According to the embodiment of the disclosure, the initial grid speed model is subjected to multi-scale thinning and encryption, so that the grid structure far away from the initial ray path is changed into the grid structure adjacent to the initial ray path, ray tracing in a remote geological structure is realized, a tracing blind area in a conventional piecewise iteration method is eliminated, a ray tracing result is closer to a real condition, and the accuracy of ray tracing is improved. In addition, as the ray path points in the diluted gridding speed model correspond to the path points of the initial ray path in the embodiment of the disclosure, the new path points obtained by using difference value calculation can still keep the attribute values in the original model when the gridding speed model is encrypted later, the accuracy of ray tracing is improved, the calculated amount in the data processing process is reduced, and the efficiency of the ray tracing method is improved.
FIG. 10 is a flowchart of a ray tracing method according to another embodiment of the present disclosure, as shown in FIG. 10, the method includes the following steps
S1001, determining an initial ray in an initial meshing speed model.
S1002, determining the number of intersection points of the initial ray paths in a preset area and grid lines in the initial grid speed model.
And S1003, if the number of the intersection points is larger than a preset threshold value, performing thinning processing on the initial meshing speed model to obtain a first meshing speed model.
The predetermined region may be a region of uniform medium in the earth medium space in which the seismic wave propagation velocity is uniform. And the gridding speed model is thinned until only the ray path points smaller than or equal to a preset threshold value exist in a medium uniform area, so that the accuracy of the path point positions after the iterative calculation section by section can be ensured to the greatest extent.
And S1004, updating the position of each intersection point of the initial ray path and the grid lines in the first meshing speed model to obtain a plurality of first ray path points, and determining a first ray path.
S1005, carrying out encryption processing on the first meshing speed model to obtain a second meshing speed model.
S1006, updating the position of each intersection point of the grid lines in the first ray path and the second gridding speed model to obtain a plurality of second ray path points, and determining a second ray path.
Specifically, the implementation principles and the specific implementation processes of S1004 to 1006 are consistent, and are not described herein.
S1007, rotating the second ray path by a preset angle towards a preset direction to obtain a rotated second ray path.
Fig. 11 is a schematic diagram of a rotated second ray path provided by an embodiment of the present disclosure. After the second ray paths in the second meshing speed model are obtained, the second ray paths are rotated by a preset angle, and the rotated second ray paths and the corresponding second meshing speed model are obtained as shown in fig. 11. During the real seismic wave propagation process, a reverse wave may occur, and the corresponding ray path is turned back inside a certain grid, that is, after entering the grid from one grid line, the ray exits the grid from the same grid line again to enter the next adjacent grid. The piecewise iteration method is based on the refraction law, and cannot solve the problem that waves are emitted from the same interface after entering from the interface, so that the piecewise iteration method cannot track the return waves.
As shown in fig. 11, the second ray path is rotated 45 degrees in the clockwise direction, and the rotated second ray path and the corresponding second meshing speed model are obtained. It will be appreciated that the rotation angle and direction provided in the embodiments of the present disclosure are merely examples, and are not limiting of the present disclosure, and may be adjusted accordingly in actual situations according to the size, parameters, complexity, and the like of the meshing speed model.
And S1008, covering the rotated second ray path through a third gridding speed model, wherein the grid line spacing in the third gridding speed model is the same as the grid line spacing in the initial gridding speed model.
Fig. 12 is a schematic diagram of a second ray path covered by a third meshing speed model according to an embodiment of the disclosure. As shown in fig. 12, a third gridding speed model is used to cover the second ray paths and the corresponding second gridding speed model, wherein the directions of the gridlines in the third gridding speed model are horizontal or vertical directions, and the grid line intervals are the same as those in the initial gridding speed model. The point E1 and the point E17 correspond to the source point P1 and the detector point P16 in the initial gridding velocity model, respectively, and the intersection points E2 to E16 of the rotated second ray path and the gridlines in the third gridding velocity model are the path points of the second ray path in the third gridding velocity model.
S1009, updating the position of each intersection point of the rotated second ray path and the grid line in the third gridding speed model to obtain a plurality of third ray path points, and determining a third ray path.
Fig. 13 is a schematic diagram of a third ray path provided by an embodiment of the present disclosure. And updating the positions of E2 to E16 according to the seismic wave propagation speeds of the second ray path on two sides of the grid line where the path points in the third gridding speed model are positioned by using a piecewise iteration method to obtain F2 to F16 shown in FIG. 13, wherein the point F1 and the point F17 correspond to a seismic source point P1 and a detection point P16 in the initial gridding speed model respectively, and further determining the third ray path shown in FIG. 13.
S1010, rotating the third ray path by a preset angle towards the opposite direction of the preset direction to obtain a fourth ray path.
Fig. 14 is a fourth ray path schematic provided by an embodiment of the present disclosure. And rotating the third ray path by 45 degrees in the anticlockwise direction to obtain a fourth ray path shown in 14 and a second gridding speed model corresponding to the fourth ray path. The point G1 and the point G8 correspond to the source point P1 and the detector point P16 in the initial gridding velocity model, respectively, and the points G2 to G7 are intersections of the fourth ray path and the second gridding velocity model, that is, the path points of the fourth ray path in the second gridding velocity model.
And repeatedly encrypting the gridding speed model according to the steps, and iteratively updating the ray paths in the gridding speed model until the grid line spacing of the encrypted gridding speed model is the same as that of the initial gridding speed model. Specifically, each pair of gridding speed models is subjected to encryption operation once, the ray paths in the gridding speed models are rotated after iteration section by section, the ray paths are rotated to the original angle in the opposite direction after iteration section by section again, and the next encryption operation is continued.
Fig. 15 is a schematic diagram of a second forward model according to an embodiment of the disclosure. The second forward model provided in the embodiment of the present disclosure is a two-dimensional model with a depth of 5000 meters and a length of 8000 meters, the propagation speed of the seismic wave increases linearly along the arrow direction, the speed value V1 at (0, 0) is 2000m/S, the speed value V2 at (5000 ) is 5000m/S, the source point S2 (2000,500) and the detector point R2 (6000,500) are set, and ray tracing is performed in the second forward model according to the method described in the above embodiment. Fig. 16 is a schematic diagram of a second forward model ray tracing result according to an embodiment of the disclosure. As shown in fig. 16, the solid curve represents the real ray path, the dash-dot line represents the ray tracing result obtained by the conventional piecewise iteration method, the dotted curve represents the ray tracing result obtained by the ray tracing method provided by the embodiment of the present disclosure, and it can be seen that the ray tracing result obtained by the ray tracing method provided by the embodiment of the present disclosure is closer to the real ray path.
According to the embodiment of the disclosure, on the basis of carrying out multi-scale thinning and encryption on the initial grid speed model, ray tracking in a remote geological structure is realized, the direction of grid lines is changed by obtaining the multi-directional grid speed model, rays which do not meet the refraction law under the condition that the refraction wave and the like are in the original grid direction are adjusted to a state capable of being iteratively updated section by section, so that ray tracking on the refraction wave and the like is realized, dead zones in the ray tracking process are further reduced, the tracking effect of the ray tracking in the complex geological structure is improved, and the data processing effect of the seismic exploration and other works of the complex area is obviously improved.
Fig. 17 is a schematic structural diagram of a ray tracing apparatus according to an embodiment of the disclosure. The ray tracing apparatus may be an electronic device having a data processing function as described in the above embodiment, or the ray tracing apparatus may be a part or component in the electronic device having a data processing function. The ray tracing apparatus provided in the embodiment of the present disclosure may execute the processing flow provided in the embodiment of the ray tracing method, as shown in fig. 17, the ray tracing apparatus 170 includes: a first determination module 171, a second determination module 172; wherein the first determining module 171 is configured to determine an initial ray in the initial meshing speed model; the second determining module 172 is configured to update the initial meshing speed model according to a preset scaling, and determine a path of the initial ray in the updated initial meshing speed model.
Optionally, the second determining module 172 includes a thinning unit 1721, a path determining unit 1722, and an encrypting unit 1723; the thinning unit 1721 is configured to perform thinning processing on the initial gridding speed model to obtain a first gridding speed model, where a grid line pitch in the first gridding speed model is greater than a grid line pitch in the initial gridding speed model; the path determining unit 1722 is configured to update a position of each intersection point of the initial ray path and the grid line in the first meshing speed model to obtain a plurality of first ray path points, and determine a first ray path; the encryption unit 1723 is configured to encrypt the first meshing speed model to obtain a second meshing speed model, where the grid line pitch in the second meshing speed model is smaller than the grid line pitch in the first meshing speed model; the path determining unit 1722 is further configured to update a position of each intersection point of the first ray path and the grid line in the second meshing speed model, obtain a plurality of second ray path points, and determine a second ray path.
Optionally, the thinning unit 1721 is further configured to determine the number of intersections between the initial ray paths and grid lines in the initial meshing speed model in a preset area; and if the number of the intersection points is larger than a preset threshold value, performing thinning treatment on the initial meshing speed model to obtain a first meshing speed model.
Optionally, the second determining module 172 further includes a rotating unit 1724, configured to rotate the second ray path by a preset angle towards a preset direction, to obtain a rotated second ray path; covering the rotated second ray path by a third gridding speed model, wherein the grid line spacing in the third gridding speed model is the same as the grid line spacing in the initial gridding speed model; the path determining unit 1722 is further configured to update a position of each intersection point of the rotated second ray path and the grid line in the third meshing speed model to obtain a plurality of third ray path points, and determine a third ray path; the rotating unit 1724 is further configured to rotate the third ray path by a preset angle in a direction opposite to the preset direction, so as to obtain a fourth ray path.
Optionally, the encrypting unit 1723 is further configured to encrypt the second meshing speed model if the grid line pitch in the second meshing speed model is smaller than the grid line pitch in the initial meshing speed model; the path determination unit 1722 is further configured to update the rotated fourth ray path according to the encrypted second meshing speed model and the third meshing speed model.
The ray tracing apparatus of the embodiment shown in fig. 17 may be used to implement the technical solution of the embodiment of the ray tracing method described above, and its implementation principle and technical effects are similar, and will not be described herein again.
Fig. 18 is a schematic structural diagram of an electronic device according to an embodiment of the disclosure. The electronic device may be an electronic device having a data processing function as described in the above embodiment. The electronic device provided in the embodiment of the present disclosure may execute the processing flow provided in the embodiment of the ray tracing method, as shown in fig. 18, the electronic device 180 includes: memory 181, processor 182, computer programs, and communication interface 183; wherein the computer program is stored in the memory 181 and configured to perform the ray tracing method as described above by the processor 182.
In addition, the embodiment of the present disclosure also provides a computer-readable storage medium having stored thereon a computer program that is executed by a processor to implement the ray tracing method described in the above embodiment.
Furthermore, embodiments of the present disclosure provide a computer program product comprising a computer program or instructions which, when executed by a processor, implement a ray tracing method as described above.
Furthermore, the disclosed embodiments also provide a computer program product comprising a computer program or instructions which, when executed by a processor, implements a vehicle voice control method as described above.
Computer program code for carrying out operations of the present disclosure may be written in one or more programming languages, including, but not limited to, an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).
It should be noted that in this document, relational terms such as "first" and "second" and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The foregoing is merely a specific embodiment of the disclosure to enable one skilled in the art to understand or practice the disclosure. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the disclosure. Thus, the present disclosure is not intended to be limited to the embodiments shown and described herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (10)
1. A method of ray tracing, the method comprising:
determining an initial ray in an initial meshing speed model;
updating the initial gridding speed model according to a preset scaling ratio, and determining the path of the initial ray in the updated initial gridding speed model.
2. The method of claim 1, wherein updating the initial meshing speed model at a preset scale, determining a ray path of the initial ray in the updated initial meshing speed model, comprises:
performing thinning treatment on the initial gridding speed model to obtain a first gridding speed model, wherein the grid line spacing in the first gridding speed model is larger than that in the initial gridding speed model;
updating the position of each intersection point of the initial ray path and the grid lines in the first gridding speed model to obtain a plurality of first ray path points, and determining a first ray path;
encrypting the first gridding speed model to obtain a second gridding speed model, wherein the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the first gridding speed model;
and updating the position of each intersection point of the grid lines in the first ray path and the second meshing speed model to obtain a plurality of second ray path points, and determining a second ray path.
3. The method of claim 2, wherein the thinning the initial meshing speed model to obtain a first meshing speed model includes:
determining the number of intersection points of the initial ray paths in a preset area and grid lines in the initial gridding speed model;
and if the number of the intersection points is larger than a preset threshold value, performing thinning treatment on the initial meshing speed model to obtain a first meshing speed model.
4. The method of claim 2, wherein after determining the second ray path, the method further comprises:
rotating the second ray path by a preset angle towards a preset direction to obtain a rotated second ray path;
covering the rotated second ray path by a third gridding speed model, wherein the grid line spacing in the third gridding speed model is the same as the grid line spacing in the initial gridding speed model;
updating the position of each intersection point of the rotated second ray path and the grid lines in the third meshing speed model to obtain a plurality of third ray path points, and determining a third ray path;
and rotating the third ray path by a preset angle towards the reverse direction of the preset direction to obtain a fourth ray path.
5. The method of claim 4, wherein after the fourth ray path is obtained, the method further comprises:
if the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the initial gridding speed model, encrypting the second gridding speed model;
and updating the rotated fourth ray path according to the encrypted second meshing speed model and the third meshing speed model.
6. A ray tracing apparatus, comprising:
a first determining module for determining an initial ray in an initial gridding speed model;
and the second determining module is used for updating the initial gridding speed model according to a preset scaling ratio and determining the path of the initial ray in the updated initial gridding speed model.
7. The apparatus of claim 6, wherein the second determination module is further configured to:
performing thinning treatment on the initial gridding speed model to obtain a first gridding speed model, wherein the grid line spacing in the first gridding speed model is larger than that in the initial gridding speed model;
updating the position of each intersection point of the initial ray path and the grid lines in the first gridding speed model to obtain a plurality of first ray path points, and determining a first ray path;
encrypting the first gridding speed model to obtain a second gridding speed model, wherein the grid line spacing in the second gridding speed model is smaller than the grid line spacing in the first gridding speed model;
and updating the position of each intersection point of the grid lines in the first ray path and the second meshing speed model to obtain a plurality of second ray path points, and determining a second ray path.
8. The apparatus of claim 7, wherein the second determination module is further configured to:
rotating the second ray path by a preset angle towards a preset direction to obtain a rotated second ray path;
covering the rotated second ray path by a third gridding speed model, wherein the grid line spacing in the third gridding speed model is the same as the grid line spacing in the initial gridding speed model;
updating the position of each intersection point of the rotated second ray path and the grid lines in the third meshing speed model to obtain a plurality of third ray path points, and determining a third ray path;
and rotating the third ray path by a preset angle towards the reverse direction of the preset direction to obtain a fourth ray path.
9. An electronic device, comprising:
a memory;
a processor; and
a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement the method of any one of claims 1-5.
10. A computer readable storage medium, on which a computer program is stored, which computer program, when being executed by a processor, implements the method according to any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211151458.8A CN117784223A (en) | 2022-09-21 | 2022-09-21 | Ray tracing method, device, equipment and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211151458.8A CN117784223A (en) | 2022-09-21 | 2022-09-21 | Ray tracing method, device, equipment and computer readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117784223A true CN117784223A (en) | 2024-03-29 |
Family
ID=90387692
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211151458.8A Pending CN117784223A (en) | 2022-09-21 | 2022-09-21 | Ray tracing method, device, equipment and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117784223A (en) |
-
2022
- 2022-09-21 CN CN202211151458.8A patent/CN117784223A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rasht‐Behesht et al. | Physics‐informed neural networks (PINNs) for wave propagation and full waveform inversions | |
| Givoli | Non-reflecting boundary conditions | |
| CN110945386B (en) | Angle domain co-imaging gather based on reverse time migration modeling | |
| He et al. | Reparameterized full-waveform inversion using deep neural networks | |
| EP2966602A1 (en) | Regional stress inversion using frictional faults | |
| EP3039458B1 (en) | Unified application framework for finite-difference modeling | |
| CN108983285B (en) | moment tensor-based multi-seismic source wave field simulation method and device | |
| JP7142968B2 (en) | FULL WAVEFORM INVERSION METHOD, APPARATUS AND ELECTRONICS | |
| US11169287B2 (en) | Method and system for automated velocity model updating using machine learning | |
| Zheng et al. | A method for identifying three-dimensional rock blocks formed by curved fractures | |
| CN116774292B (en) | Seismic wave travel time determining method, system, electronic equipment and storage medium | |
| WO2016001697A1 (en) | Systems and methods for geologic surface reconstruction using implicit functions | |
| Gadylshin et al. | Numerical dispersion mitigation neural network for seismic modeling | |
| CN109709602B (en) | Remote detection sound wave offset imaging method, device and system | |
| Haghighat et al. | En-DeepONet: An enrichment approach for enhancing the expressivity of neural operators with applications to seismology | |
| CN117784223A (en) | Ray tracing method, device, equipment and computer readable storage medium | |
| Potter et al. | Ordered line integral methods for solving the eikonal equation | |
| CN113221409B (en) | Two-dimensional numerical simulation method and device for acoustic waves with coupled finite elements and boundary elements | |
| GB2468023A (en) | Seismic attributes for structural analysis | |
| Zhang et al. | Modeling and visualization of tsunamis | |
| Gillberg | Fast and accurate front propagation for simulation of geological folds | |
| Gandham | High performance high-order numerical methods: applications in ocean modeling | |
| Pletzer et al. | Mimetic interpolation of vector fields on Arakawa C/D grids | |
| US8494778B2 (en) | Variable grid for finite difference computation | |
| Betteridge | Efficient elliptic solvers for higher order DG discretisations on modern architectures and applications in atmospheric modelling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |