A ray tracing area analysis method and apparatus
By constructing a uniformly symmetric box in the 3D geological model and determining the ray intersection points, the problem of low ray tracing calculation efficiency is solved, realizing efficient and accurate ray tracing analysis and meeting the real-time analysis needs of the geological model.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA NAT PETROLEUM CORP
- Filing Date
- 2023-11-21
- Publication Date
- 2026-06-26
Smart Images
Figure CN120028831B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of geophysical exploration, and in particular to a method and apparatus for ray tracing regional analysis. Background Technology
[0002] Ray tracing is currently the most widely used forward modeling technique for seismic waves in practical applications. This technique is based on fundamental kinematic theories of seismic wave propagation and computer technology. In real-world applications, ray tracing forward modeling of 3D geological models can provide strong technical support for optimizing seismic observation schemes in complex areas. Ray tracing can obtain not only the travel time information of seismic waves but also the ray trajectory, which has significant theoretical and practical value for research in the field of inversion tomography.
[0003] The ray tracing analysis method for geological models is based on the technique of finding the intersection points of rays with the three-dimensional geological model. By obtaining the intersection points of rays in various strata of the complex three-dimensional geological model and the endpoints of ray segments on the surface of the geological model, the path of rays emitted from different points in the strata and the imaging range on the surface can be determined, providing technical support for the exploration of complex three-dimensional geological models. Summary of the Invention
[0004] This invention reveals that when using conventional techniques for finding intersections between rays and 3D geological models for seismic ray tracing, the ray tracing process is time-consuming and computationally inefficient when the number of triangles or rays in the 3D geological model is large. This makes it difficult to meet the needs of real-time ray tracing analysis based on geological models. Therefore, how to quickly and accurately find the intersections of rays in the 3D geological model to improve project construction efficiency has become an urgent problem to be solved.
[0005] In view of the above problems, the present invention is proposed to provide a ray tracing region analysis method and apparatus that overcomes or at least partially solves the above problems.
[0006] In a first aspect, embodiments of the present invention provide a ray tracing region analysis method, comprising:
[0007] Obtain a three-dimensional geological model of the area to be explored; the three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple homosymmetric boxes are pre-constructed in each geological block. Each homosymmetric box includes triangles in the triangular mesh that correspond to the homosymmetric box.
[0008] The initial ray is obtained based on the selected ray origin and the preset deflection angle of the formation normal at the relative ray origin;
[0009] Determine the symmetrical box within the geological block where the ray originates and intersects with the initial ray;
[0010] Based on the symmetrical boxes that intersect with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block, and obtain the set of intersection points.
[0011] Select intersection points that meet the preset conditions from the set of intersection points, and use them as the path endpoint of the initial ray in the geological block where the ray originates, as well as the ray origin of the adjacent geological blocks of the geological block where the ray originates;
[0012] If the endpoint of the path is not located on the surface of the 3D geological model, the next ray is determined based on the stratum velocity of the two adjacent geological blocks and the ray origin of the adjacent geological blocks. The process continues to determine the symmetrical box in the current geological block where the ray origin is located that intersects with the initial ray until the intersection point with the geological blocks in the 3D geological model, including the surface, is obtained, thus obtaining the endpoint of the ray.
[0013] In some alternative embodiments, the process of constructing multiple homosymmetric boxes within a geological block includes:
[0014] Based on the shape parameters of the triangles included in the triangular mesh of the geological block, determine the bounding box shape parameters of the geological block;
[0015] Based on the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or based on the shape parameters of the bounding box and the preset division ratio, the bounding box is divided into multiple initial symmetrical boxes.
[0016] Determine whether the number of triangles included in the initial symmetrical box meets the preset number requirement; if not, adjust the shape parameters of the initial symmetrical box or adjust the division ratio, and then execute the step of dividing the bounding box; if yes, obtain multiple symmetrical boxes included in the geological block.
[0017] In some alternative embodiments, the triangle shape parameters include vertices; the bounding box shape parameters include the length, width, and height of the bounding box; and the initial symmetric box shape parameters include the length, width, and height of the initial symmetric box.
[0018] In some alternative embodiments, determining the number of triangles included in the initial homosymmetric box includes:
[0019] Determine the total number of initial homosymmetric boxes for the geological blocks;
[0020] Determine the total number of triangles in the triangulation of the geological block;
[0021] The number of triangles in the triangulation of the geological blocks included in each initial homosymmetric box can be determined using the following formula:
[0022] triPerCell = totalTri / cells, where triPerCell represents the number of triangles in the triangulation of the geological block included in each initial homosymmetric box, totalTri represents the total number of triangles in the triangulation of the geological block, and cells represents the total number of initial homosymmetric boxes.
[0023] Adjust the shape parameters of the initial symmetrical box, including:
[0024] Adjust the length, width, and height of the initial symmetrical box using the following formulas: L1 = L / n, W1 = W / n, H1 = H / n, where L1, W1, and H1 are the adjusted length, width, and height of the initial symmetrical box, respectively; L, W, and H are the original length, width, and height of the initial symmetrical box; and n is the adjustment ratio.
[0025] In some optional embodiments, determining the total number of initial homosymmetric boxes for the geological blocks includes:
[0026] The number of initial symmetrical boxes in the X direction of the geological block is determined based on the length of the bounding box and the length of each initial symmetrical box.
[0027] The number of initial symmetrical boxes in the Y direction for a geological block is determined based on the width of the bounding box and the width of each initial symmetrical box.
[0028] The number of initial symmetrical boxes in the Z direction of the geological block is determined based on the height of the enclosing box and the height of each initial symmetrical box.
[0029] Multiply the number of initial homosymmetric boxes in the X, Y, and Z directions for each geological block, and obtain the total number of initial homosymmetric boxes constructed from the geological block based on the result of the multiplication.
[0030] In some alternative embodiments, determining the perfectly symmetric box within the geological block where the ray originates and intersects with the initial ray includes:
[0031] Determine the geological block where the ray originates;
[0032] Determine all homosymmetric boxes contained in the geological block where the ray originates, and establish a set of homosymmetric boxes;
[0033] Determine the homosymmetric boxes that intersect the initial ray with the set of homosymmetric boxes, and obtain the homosymmetric boxes that intersect with the initial ray.
[0034] In some alternative embodiments, determining the geological block where the ray originates includes:
[0035] Based on the triangulation where the ray originates, obtain geological blocks B1 and B2 that are adjacent to the triangulation.
[0036] The position coordinates of the ray starting point are finely adjusted according to the preset adjustment rules;
[0037] Using the finely adjusted starting point of the ray as the endpoint, construct a vertical segment b1 that extends beyond geological block B1 and a vertical segment b2 that extends beyond geological block B2.
[0038] Determine whether the number of intersections between perpendicular segment b1 and geological block B1 and the number of intersections between perpendicular segment b2 and geological block B2 are odd. The geological block with an odd number of intersections is the geological block where the ray originates.
[0039] In some optional embodiments, based on the symmetrical boxes intersecting the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, resulting in a set of intersection points, including:
[0040] Traverse the homosymmetric boxes that intersect with the initial ray;
[0041] Determine whether the initial ray intersects with a triangle in the triangular mesh contained in an intersecting symmetrical box;
[0042] If so, determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block, and add the determined intersection points to the intersection point set;
[0043] After traversing all the homogeneous symmetric boxes that intersect with the initial ray, we obtain the set of intersection points.
[0044] In some optional embodiments, determining the intersection of the initial ray with the triangles in the triangulation of the current geological block includes:
[0045] Given the three vertices p1, p2, and p3 of the triangle, calculate the side vectors v1 and v2 of the triangle.
[0046] Determine the cross product xv of the edge vectors based on the edge vectors;
[0047] The ray segment vector vr is determined based on the starting point r1 and the ending point r2 of the initial ray;
[0048] Determine vector v3 based on the vertex p1 of the triangle and the starting point r1 of the initial ray;
[0049] The dot product d1 of the two vectors is determined by the cross product xv and the vector vr.
[0050] The dot product d2 of vectors v3 and xv is determined based on the vectors v3 and xv.
[0051] The point in the triangle is determined by the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle;
[0052] Determine the area A of the triangle formed by vertices p1, p2, and p3, the area a1 of the triangle formed by vertices p1, p2, and zp, the area a2 of the triangle formed by vertices p2, p3, and zp, and the area a3 of the triangle formed by p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then point zp in the triangle is the intersection of the initial ray and the triangle in the triangular mesh of the current geological block.
[0053] In some optional embodiments, the above method further includes:
[0054] Based on the endpoint of the rays determined by multiple initial rays passing through the ray origin with a preset deflection angle, the imaging range of the ground surface included in the three-dimensional geological model is obtained by emitting rays with the ray origin as the starting point and the preset deflection angle as the direction.
[0055] In some optional embodiments, multiple initial rays passing through a preset deflection angle at the ray origin include:
[0056] Based on the selected ray starting point and the preset deflection angle of the stratum normal at the relative ray starting point, multiple initial rays are obtained by using the preset deflection angle of the stratum normal as the ray direction and distributing them evenly around the circumference at preset intervals.
[0057] Secondly, embodiments of the present invention provide a ray tracing region analysis device, comprising:
[0058] The data acquisition module is used to acquire a three-dimensional geological model of the area to be explored. The three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple homosymmetric boxes are pre-constructed in each geological block. Each homosymmetric box includes triangles in the triangular mesh that correspond to the homosymmetric box.
[0059] The ray determination module is used to obtain an initial ray based on the selected ray origin and the preset deflection angle of the formation normal at the relative ray origin; and to determine the next ray based on the formation velocity of two adjacent geological blocks and the ray origin of the adjacent geological blocks.
[0060] The intersection point determination module is used to determine the symmetrical boxes in the geological block where the ray originates and intersects with the initial ray; based on the symmetrical boxes intersecting with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, resulting in an intersection point set; from the intersection point set, the intersection points that meet preset conditions are selected as the path endpoint of the initial ray in the geological block where the ray originates, and the ray origin of the adjacent geological block of the geological block where the ray originates; if the path endpoint is not located on the surface of the 3D geological model, the ray determination module is notified and, based on the next ray determined by the ray determination module, the intersection point determination module is notified to continue the step of determining the symmetrical boxes in the current geological block where the ray originates and intersects with the initial ray, until the intersection points with the geological blocks including the surface in the 3D geological model are obtained, thus obtaining the ray endpoint.
[0061] In some optional embodiments, the above-described apparatus further includes:
[0062] The imaging range determination module is used to obtain the imaging range of the ground surface included in the three-dimensional geological model based on the ray endpoint determined by multiple initial rays passing through the ray origin.
[0063] This invention also provides a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement a ray tracing region analysis method.
[0064] This invention also provides a terminal device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement a ray tracing region analysis method.
[0065] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:
[0066] The method provided in this invention divides the three-dimensional geological model of the area to be explored into triangular facets, and constructs a uniformly symmetric box based on the facet division. An initial ray is constructed based on the selected ray origin and a preset deflection angle relative to the stratigraphic normal at the ray origin. Intersections are performed using the uniformly symmetric boxes to determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block, resulting in a set of intersection points. This eliminates the need for intersection calculations on the entire geological block, significantly reducing the computational load and improving processing speed and efficiency. Intersection points meeting preset conditions are selected from the set of intersection points as the path endpoint of the initial ray in the geological block where the ray origin is located. The path endpoint of the current geological block is used as the ray origin for the next adjacent geological block, and the next ray is constructed. Intersection processing continues for the next adjacent geological block based on the uniformly symmetric boxes included in the adjacent geological blocks, until an intersection point is obtained with a geological block in the three-dimensional geological model including the surface, thus obtaining the ray endpoint. This method, when determining the intersection points of rays with various strata and geological blocks, first identifies the homosymmetric boxes through which the ray passes and then calculates the intersection points. For homosymmetric boxes that the ray does not pass through, there is no need to calculate the intersection points, which greatly reduces the amount of data processing required for intersection calculations and improves the processing speed and efficiency. It can solve the problems of high time consumption, poor computational efficiency, and difficulty in meeting the needs of real-time ray tracing analysis based on geological models in existing technologies, improves the accuracy and efficiency of ray tracing, and provides strong technical support for seismic ray tracing analysis.
[0067] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.
[0068] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description
[0069] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:
[0070] Figure 1 This is a flowchart of the ray tracing region analysis method in Embodiment 1 of the present invention.
[0071] Figure 2 This is a flowchart of the ray tracing region analysis method in Embodiment 2 of the present invention;
[0072] Figure 3 This is a schematic diagram of a three-dimensional geological model in Embodiment 2 of the present invention;
[0073] Figure 4 This is a schematic diagram of a geological block in the three-dimensional geological model of Embodiment 2 of the present invention;
[0074] Figure 5 This is a schematic diagram of a perfectly symmetrical box representing a geological block in Embodiment 2 of the present invention;
[0075] Figure 6 This is a schematic diagram illustrating the determination of the intersection point between the initial ray and the triangle in Embodiment 2 of the present invention;
[0076] Figure 7 This is a schematic diagram of ray tracing in Embodiment 2 of the present invention;
[0077] Figure 8 This is a schematic diagram of the imaging range of the initial ray from a ray starting point on the ground surface in Embodiment 2 of the present invention;
[0078] Figure 9 This is a schematic diagram of the imaging range of the initial rays from multiple ray initiation points on the earth's surface in Embodiment 2 of the present invention;
[0079] Figure 10 This is a structural diagram of the ray tracing region analysis device in Embodiment 2 of the present invention. Detailed Implementation
[0080] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
[0081] Ray tracing technology plays a crucial role in seismic exploration, with wide applications in seismic location, geophysical tomography, and other seismic digital processing and inversion. The speed and accuracy of ray tracing are key to solving practical production problems. However, existing ray tracing techniques suffer from high tracing time and poor computational efficiency when dealing with a large number of rays, making it difficult to meet the needs of real-time ray tracing analysis based on geological models. To address these issues, this invention provides a ray tracing regional analysis method that overcomes the shortcomings of existing technologies and provides strong technical support for exploration in complex geological areas.
[0082] Example 1
[0083] Embodiment 1 of the present invention provides a method for ray tracing region analysis, the process of which is as follows: Figure 1 As shown, it includes the following steps:
[0084] S101: Obtain a three-dimensional geological model of the area to be explored; the three-dimensional geological model includes multiple layers and multiple geological blocks, each layer includes multiple triangular meshes, each geological block is a closed body enclosed by triangular meshes from multiple layers, each triangular mesh includes multiple triangles, and multiple symmetrical boxes are pre-constructed in each geological block, each symmetrical box includes the triangle in the triangular mesh corresponding to the symmetrical box;
[0085] S102: Based on the selected ray starting point and the preset deflection angle of the formation normal at the relative ray starting point, the initial ray is obtained;
[0086] S103: Determine the symmetrical box in the geological block where the ray originates and intersects with the initial ray;
[0087] S104: Based on the symmetrical boxes that intersect with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, determine the intersection points of the initial ray and the triangles in the triangulation of the current geological block, and obtain the set of intersection points.
[0088] S105: Select intersection points that meet the preset conditions from the set of intersection points, and use them as the path endpoint of the initial ray in the geological block where the ray originates, as well as the ray origin of the adjacent geological blocks of the geological block where the ray originates.
[0089] S106: If the endpoint of the path is not located on the surface of the 3D geological model, determine the next ray based on the stratum velocity of the two adjacent geological blocks and the ray starting point of the adjacent geological blocks. Continue to execute the step of determining the symmetrical box in the current geological block where the ray starting point is located that intersects with the initial ray until the intersection point with the geological blocks in the 3D geological model including the surface is obtained, and the endpoint of the ray is obtained.
[0090] The method provided in this invention divides the three-dimensional geological model of the area to be explored into triangular facets, and constructs a uniformly symmetric box based on the facet division. An initial ray is constructed based on the selected ray origin and a preset deflection angle relative to the stratigraphic normal at the ray origin. Intersections are performed using the uniformly symmetric boxes to determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block, resulting in a set of intersection points. This eliminates the need for intersection calculations on the entire geological block, significantly reducing the computational load and improving processing speed and efficiency. Intersection points meeting preset conditions are selected from the set of intersection points as the path endpoint of the initial ray in the geological block where the ray origin is located. The path endpoint of the current geological block is used as the ray origin for the next adjacent geological block, and the next ray is constructed. Intersection processing continues for the next adjacent geological block based on the uniformly symmetric boxes included in the adjacent geological blocks, until an intersection point is obtained with a geological block in the three-dimensional geological model including the surface, thus obtaining the ray endpoint. This method, when determining the intersection points of rays with various strata and geological blocks, first identifies the homosymmetric boxes through which the ray passes and then calculates the intersection points. For homosymmetric boxes that the ray does not pass through, there is no need to calculate the intersection points again, which greatly reduces the amount of data processing required for intersection calculations and improves the processing speed and efficiency. It can solve the problems of high time consumption, poor computational efficiency, and difficulty in meeting the needs of real-time ray tracing analysis based on geological models in existing technologies, improves the accuracy and efficiency of ray tracing, and provides strong technical support for seismic ray tracing analysis.
[0091] Example 2
[0092] Embodiment 2 of the present invention provides a ray tracing region analysis method. After determining the ray endpoint based on the method provided in Embodiment 1, the ray starting point is further determined as the starting point, and the imaging range of the ray emitted at a preset deflection angle is determined. The process is as follows: Figure 2 As shown, it includes the following steps:
[0093] Step S201: Obtain a three-dimensional geological model of the area to be explored.
[0094] Step S202: Based on the selected ray starting point and the preset deflection angle of the stratum normal at the relative ray starting point, obtain the initial ray.
[0095] Step S203: Determine the symmetrical box in the geological block where the ray originates and intersects with the initial ray.
[0096] Step S204: Based on the symmetrical boxes that intersect with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, determine the intersection points of the initial ray and the triangles in the triangulation of the current geological block, and obtain the set of intersection points.
[0097] Step S205: Select intersection points that meet the preset conditions from the set of intersection points, as the path endpoint of the initial ray in the geological block where the ray originates, and the ray origin of the adjacent geological block of the geological block where the ray originates.
[0098] Step S206: Determine whether the endpoint of the path is located on the surface of the three-dimensional geological model. If not, proceed to step S207; if so, proceed to step S208.
[0099] Step S207: Based on the formation velocity of two adjacent geological blocks and the ray starting point of the adjacent geological blocks, determine the next ray and continue to execute step S203.
[0100] Step S208: Take the endpoint of the path located on the surface of the three-dimensional geological model as the endpoint of the ray.
[0101] Step S209: Based on the endpoint of the rays determined by multiple initial rays passing through the ray origin at a preset deflection angle, obtain the imaging range of the ground surface included in the three-dimensional geological model by rays emitted from the ray origin at the preset deflection angle.
[0102] The specific implementation process of each step in Embodiment 1 and Embodiment 2 above is described in detail below.
[0103] Preferably, in steps S101 and S201 above, the three-dimensional geological model of the area to be explored is described by a structure of "volume → block → surface → triangular mesh → triangle". The physical property interfaces between different geological blocks are composed of triangular meshes. The three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple symmetrical boxes are pre-constructed in each geological block. Each symmetrical box includes a triangle in the triangular mesh that corresponds to the symmetrical box.
[0104] Figure 3 This is a schematic diagram of a three-dimensional geological model. Figure 3 The three-dimensional geological model contains six strata, and the closed area enclosed by the strata is called a geological block. Figure 4 This is a schematic diagram of a geological block in a three-dimensional geological model. Figure 4 The top face of the geological block is composed of 7 triangular meshes, and the sides and front are each composed of one triangular mesh. Correspondingly, the other 3 faces of the geological block are also composed of triangular meshes. Each triangular mesh contains several triangular facets. The regular cube outside the geological block is the bounding box of the geological block.
[0105] Preferably, multiple homosymmetric boxes are constructed within the geological block, including:
[0106] Based on the shape parameters of the triangles included in the triangular mesh of the geological block, determine the bounding box shape parameters of the geological block;
[0107] Based on the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or based on the shape parameters of the bounding box and the preset division ratio, the bounding box is divided into multiple initial symmetrical boxes.
[0108] Determine whether the number of triangles included in the initial symmetrical box meets the preset number requirement; if not, adjust the shape parameters of the initial symmetrical box or adjust the division ratio, and then execute the step of dividing the bounding box; if yes, obtain multiple symmetrical boxes included in the geological block.
[0109] Constructing multiple uniformly symmetric boxes for a geological block can greatly reduce the amount of computation. Instead of performing intersection calculations for triangles in the entire geological block, the intersection calculations are performed on the triangles contained in the uniformly symmetric boxes that intersect with the initial ray. This can improve computational efficiency and accuracy when encountering geological blocks with a large number of triangles.
[0110] Optionally, the shape parameters of the triangles and the bounding boxes can be selected as needed. The triangle shape parameters may include vertices; the bounding box shape parameters may include the length, width, and height of the bounding box; and the initial symmetrical box shape parameters may include the length, width, and height of the initial symmetrical box. For example, in this embodiment, if the number of triangles included in each initial symmetrical box is less than 10, the length, width, and height of the initial symmetrical box are adjusted. After adjustment, the number of triangles included in the adjusted symmetrical box is determined again until the number of triangles in the initial symmetrical box is greater than 10. Then, the adjustment of the shape parameters of the initial symmetrical box is stopped, and the symmetrical box divided at this time is taken as the final symmetrical box.
[0111] Figure 5 This is a schematic diagram of a perfectly symmetrical box representing a geological block. Figure 5 The bounding box of the geological block was divided into several initial homosymmetric boxes according to the partitioning rules. Each homosymmetric box contains several triangles. Some homosymmetric boxes may contain triangles from multiple triangulations, for example... Figure 5 The first column on the left contains four initial symmetrical boxes, each containing triangles from multiple triangulations. The rays in the diagram pass through some of these initial symmetrical boxes.
[0112] Preferably, determining the number of triangles included in the initial homosymmetric box includes:
[0113] Determine the total number of initial homosymmetric boxes for the geological blocks;
[0114] Determine the total number of triangles in the triangulation of the geological block;
[0115] The number of triangles in the triangulation of the geological blocks included in each initial homosymmetric box can be determined using the following formula:
[0116] triPerCell = totalTri / cells, where triPerCell represents the number of triangles in the triangulation of the geological block included in each initial homosymmetric box, totalTri represents the total number of triangles in the triangulation of the geological block, and cells represents the total number of initial homosymmetric boxes.
[0117] Preferably, adjusting the shape parameters of the initial symmetrical box includes:
[0118] In this embodiment, the shape parameters of the initial symmetrical box are adjusted using the following formulas: L1 = L / n, W1 = W / n, H1 = H / n, where L1, W1, and H1 are the adjusted length, width, and height of the initial symmetrical box, respectively; L, W, and H are the original length, width, and height of the initial symmetrical box; and n is the adjustment ratio. For example, in this embodiment, the initial symmetrical box is adjusted by a ratio of 2. The ratio parameter can be selected according to the required level of detail. The shape parameters of the initial symmetrical box can be selected as needed and adjusted according to preset adjustment rules.
[0119] Preferably, determining the total number of initial homosymmetric boxes for the geological blocks includes:
[0120] The number of initial symmetrical boxes in the X direction of the geological block is determined based on the length of the bounding box and the length of each initial symmetrical box.
[0121] The number of initial symmetrical boxes in the Y direction for a geological block is determined based on the width of the bounding box and the width of each initial symmetrical box.
[0122] The number of initial symmetrical boxes in the Z direction of the geological block is determined based on the height of the enclosing box and the height of each initial symmetrical box.
[0123] Multiply the number of initial homosymmetric boxes for each geological block in the X, Y, and Z directions. Based on the result of this multiplication, obtain the total number of initial homosymmetric boxes constructed from the geological block. See also Figure 5 As shown, the initial homosymmetric box of this geological block is 7*4*1, and a total of 28 initial homosymmetric boxes are divided. Accordingly, the total number of homosymmetric boxes can be calculated according to the above method.
[0124] Preferably, in steps S102 and S202 above, the initial ray is oriented at a preset deflection angle, starting from the selected ray origin, and attempts to emit the initial ray towards the surface of the strata in the three-dimensional geological model. It should be noted that the ray in this invention represents the propagation path of seismic waves moving through media with different velocities, and this propagation path is represented by a virtual ray.
[0125] Preferably, steps S103 and S203 include:
[0126] Determine the geological block where the ray originates;
[0127] Determine all homosymmetric boxes contained in the geological block where the ray originates, and establish a set of homosymmetric boxes;
[0128] Determine the homosymmetric boxes that intersect the initial ray with the set of homosymmetric boxes, and obtain the homosymmetric boxes that intersect with the initial ray.
[0129] See Figure 5 As shown, with Figure 5 The ray in the image is the initial ray. Based on the set of homosymmetric boxes in the geological block, determine the four homosymmetric boxes in the set of homosymmetric boxes that intersect with the initial ray.
[0130] Preferably, the geological block where the ray originates is located includes:
[0131] Based on the triangulation where the ray originates, obtain geological blocks B1 and B2 that are adjacent to the triangulation.
[0132] The position coordinates of the ray starting point are finely adjusted according to the preset adjustment rules;
[0133] Using the finely adjusted starting point of the ray as the endpoint, construct a vertical segment b1 that extends beyond geological block B1 and a vertical segment b2 that extends beyond geological block B2.
[0134] Determine whether the number of intersections between perpendicular line segment b1 and geological block B1, and the number of intersections between perpendicular line segment b2 and geological block B2 are odd. The geological block with an odd number of intersections is the geological block where the ray originates.
[0135] In this embodiment, the adjustment rule can be to adjust the coordinates of the ray starting point in the Z direction, for example, by adjusting the coordinates of the ray starting point in the Z direction upwards by 0.25. Other adjustment rules can be designed according to actual exploration needs.
[0136] Preferably, steps S104 and S204 include:
[0137] Traverse the homosymmetric boxes that intersect with the initial ray;
[0138] Determine whether the initial ray intersects with a triangle in the triangular mesh contained in an intersecting symmetrical box;
[0139] If so, determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block, and add the determined intersection points to the intersection point set;
[0140] After traversing all the homogeneous symmetric boxes that intersect with the initial ray, we obtain the set of intersection points.
[0141] In flat strata, the initial ray intersects with a triangle in the triangulation of a perfectly symmetric box. Where strata are concave or convex, folded regions form between the strata. See [reference needed]. Figure 4 As shown, Figure 4 The top layer is not a smooth layer; there are bulges and depressions. At this time, wrinkled areas will be formed between the layers. The symmetrical box containing the wrinkled areas will contain triangles of multiple triangular meshes. Therefore, there will be multiple intersection points between the initial ray and the triangles. Add the intersection points to the set to form the intersection point set.
[0142] Preferably, determining the intersection points of the initial ray and the triangles in the triangulation of the current geological block includes:
[0143] Given the three vertices p1, p2, and p3 of the triangle, calculate the side vectors v1 and v2 of the triangle.
[0144] Determine the cross product xv of the edge vectors based on the edge vectors;
[0145] The ray segment vector vr is determined based on the starting point r1 and the ending point r2 of the initial ray;
[0146] Determine vector v3 based on the vertex p1 of the triangle and the starting point r1 of the initial ray;
[0147] The dot product d1 of the two vectors is determined by the cross product xv and the vector vr.
[0148] The dot product d2 of vectors v3 and xv is determined based on the vectors v3 and xv.
[0149] The point in the triangle is determined by the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle;
[0150] Determine the area A of the triangle formed by vertices p1, p2, and p3, the area a1 of the triangle formed by vertices p1, p2, and zp, the area a2 of the triangle formed by vertices p2, p3, and zp, and the area a3 of the triangle formed by p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then point zp in the triangle is the intersection of the initial ray and the triangle in the triangular mesh of the current geological block.
[0151] For example, in this embodiment, if the value of (a1+a2+a3) / A is not greater than 1.000001, zp can be determined as the intersection point of the initial ray and the triangle in the triangulation of the current geological block. The preset area threshold can then be adjusted to meet actual exploration needs. See the schematic diagram for determining the intersection point of the initial ray and the triangle. Figure 6 As shown in the figure, the required parameters are marked.
[0152] Preferably, in steps S105 and S205 above, it is necessary to select intersection points that meet preset conditions from the intersection point set. For example, in this embodiment, the point closest to the ray starting point in the intersection point set is selected as the path endpoint of the initial ray in the geological block where the ray starting point is located, as well as the ray starting point of the adjacent geological block of the geological block where the ray starting point is located. See Figure 4 As shown, for folded stratigraphic regions, there are multiple intersections between rays and triangulation networks in geological blocks. It is necessary to select intersections that meet the requirements. In this case, selection rules can be set. For example, the intersection closest to the ray origin can be selected. Of course, other rules can also be used to select intersections that meet other conditions.
[0153] The following is a method for determining the point closest to the origin of a ray in the set of intersection points:
[0154] Calculate the distances between the starting point of the ray and the intersection points in the set of intersection points in turn, and find the smallest distance among them. The intersection point corresponding to the smallest distance is the intersection point that meets the conditions.
[0155] Optionally, in steps S106 and S206 above, after determining the path endpoint, the determined path endpoint is judged to determine whether it is located on the surface of the three-dimensional geological model. If it is not located there, the next ray is determined based on the stratum velocity of the two adjacent geological blocks and the ray starting point of the adjacent geological blocks. Step S106 continues to execute step S103, and step S206 continues to execute S203 until the determined path endpoint is located on the surface of the three-dimensional geological model. The path endpoint is then used as the ray endpoint, and the process of finding the intersection point ends.
[0156] Because the medium in each stratum is non-uniform and the strata velocities are different, according to the principle of wave transmission, the initial ray will be transmitted when it passes through each stratum. Therefore, the direction of the initial ray propagating to the next stratum needs to be corrected based on the strata velocities of the two adjacent geological blocks. The intersection point of the ray with the next stratum is traced again. The steps of determining the symmetrical box that intersects with the initial ray in the geological block where the ray originates are repeated until the intersection point is traced close enough to or reaches the surface of the geological model. By fitting the intersection point in each stratum, a complete transmitted ray can be obtained, thus tracing the ray path of the initial ray in the strata.
[0157] Figure 7This is a schematic diagram of ray tracing. Figure 7 Each point represents the intersection of the initial ray in the strata. In the three-dimensional geological model of this embodiment, there are 6 strata, so the process of finding the intersection points needs to be performed 5 times. Each pair of intersection points forms a ray segment. Since the medium of the strata is not uniform and the strata velocity is different, the direction of the initial ray needs to be continuously corrected. Finally, all intersection points and ray segments are integrated to obtain a transmission ray.
[0158] Preferably, in step S209 above, the multiple initial rays passing through the preset deflection angle of the ray origin include:
[0159] Based on the selected ray starting point and the preset deflection angle of the stratum normal at the relative ray starting point, multiple initial rays are obtained by using the preset deflection angle of the stratum normal as the ray direction and distributing them evenly around the circumference at preset intervals.
[0160] Figure 8 A schematic diagram of the imaging range of an initial ray originating from a given ray on the Earth's surface. Figure 8 As can be seen, multiple initial rays are emitted from the selected ray starting point at a preset deflection angle. These initial rays are evenly distributed circumferentially with the preset deflection angle of the stratum normal as the ray direction and at preset interval angles. Through the ray area tracking analysis method of the present invention, the intersection point of each initial ray in each stratum can be determined; the transmitted ray is determined based on the intersection point; and the imaging range of the selected ray starting point on the ground surface included in the three-dimensional geological model is obtained based on the ray endpoint determined by multiple initial rays, thereby providing efficient and accurate technical support for seismic exploration in complex areas.
[0161] Figure 9 This is a schematic diagram showing the imaging range of the initial rays from multiple selected ray initiation points on the Earth's surface. Figure 9 As can be seen, not all initial rays will intersect with the ground surface. At the boundary of a geological block, there may be some initial rays that cannot reach the ground surface. In actual exploration, it is only necessary to obtain the endpoints of those rays that can reach the ground surface, thereby determining the imaging range of multiple initial rays emitted from the ray origin on the ground surface.
[0162] Based on the same inventive concept, embodiments of the present invention also provide a ray tracing region analysis device, which can be installed in a device having computer instruction processing capabilities, and the structure of the device is as follows. Figure 10 As shown, it includes:
[0163] The data acquisition module 11 is used to acquire a three-dimensional geological model of the area to be explored. The three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple symmetrical boxes are pre-constructed in each geological block. Each symmetrical box includes triangles in the triangular mesh corresponding to the symmetrical box.
[0164] The ray determination module 12 is used to obtain an initial ray based on the selected ray origin and the preset deflection angle of the formation normal at the relative ray origin; and to determine the next ray based on the formation velocity of two adjacent geological blocks and the ray origin of the adjacent geological blocks.
[0165] The intersection point determination module 13 is used to determine the symmetrical boxes in the geological block where the ray originates and intersects with the initial ray; based on the symmetrical boxes intersecting with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, resulting in an intersection point set; from the intersection point set, the intersection points that meet the preset conditions are selected as the path endpoint of the initial ray in the geological block where the ray originates, and the ray origin of the adjacent geological block where the ray originates; if the path endpoint is not located on the surface of the three-dimensional geological model, the ray determination module is notified and, based on the next ray determined by the ray determination module, the intersection point determination module is notified to continue the step of determining the symmetrical boxes in the current geological block where the ray originates and intersects with the initial ray, until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained, thus obtaining the ray endpoint.
[0166] Preferably, the above-mentioned device further includes:
[0167] The imaging range determination module 14 is used to obtain the imaging range of the ground surface included in the three-dimensional geological model based on the ray endpoint determined by multiple initial rays passing through the ray origin.
[0168] Preferably, the ray determination module 12 is also used to obtain multiple initial rays based on the selected ray starting point and the preset deflection angle of the stratum normal at the relative ray starting point, and the multiple initial rays are evenly distributed in a circle at preset interval angles; based on the stratum velocity of two adjacent geological blocks and the ray starting point of the adjacent geological blocks, the next ray that propagates between the multiple initial rays in different strata is determined.
[0169] Regarding the ray tracing region analysis device in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated here.
[0170] The method and apparatus described in this invention acquire a three-dimensional geological model of the area to be explored and establish a descriptive structure to construct a uniformly symmetric box of the geological blocks in the geological model; obtain an initial ray based on the ray origin selected in the strata of the geological model and the preset deflection angle of the stratum normal at the relative ray origin; determine the uniformly symmetric boxes in the geological block where the ray origin is located that intersect with the initial ray; determine the intersection points of the initial ray with the triangles in the triangulation of the current geological block to obtain a set of intersection points; select intersection points that meet preset conditions from the set of intersection points as the path endpoint of the initial ray in the geological block where the ray origin is located, and the ray origin of the next adjacent geological block; if the path endpoint is not located on the surface of the three-dimensional geological model, determine the next ray based on the stratum velocity of the two adjacent geological blocks and the ray origin of the adjacent geological blocks, and continue to execute the step of obtaining stratum intersection points until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained, thus obtaining the ray endpoint; determine the imaging range of the selected ray origin on the surface according to the ray endpoint, and determine the ray path of the initial ray in the three-dimensional geological model according to the intersection points of the initial ray with each stratum. The above-described method and apparatus overcome the problem of high time consumption in existing technologies for ray tracing, effectively improving the efficiency and accuracy of ray tracing analysis based on three-dimensional geological models, and providing technical support for exploration in complex areas.
[0171] Unless otherwise specifically stated, terms such as processing, calculation, operation, determination, display, etc., may refer to the actions and / or processes of one or more processing or computing systems or similar devices that represent the manipulation and conversion of data representing physical (e.g., electronic) quantities within the registers or memory of the processing system into other data similarly representing physical quantities within the memory, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals can be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description can be represented by voltage, current, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.
[0172] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process may be rearranged without departing from the scope of this disclosure. The appended method claims provide elements of various steps in an exemplary order and are not intended to limit the scope to the specific order or hierarchy described.
[0173] In the detailed description above, various features are combined together in a single embodiment to simplify this disclosure. This approach to disclosure should not be construed as reflecting an intention that embodiments of the claimed subject matter require more features than are explicitly stated in each claim. Rather, as reflected in the appended claims, the invention is presented with fewer features than all of the features in a single disclosed embodiment. Therefore, the appended claims are hereby explicitly incorporated into the detailed description, with each claim representing a separate preferred embodiment of the invention.
[0174] Those skilled in the art will also understand that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments herein can be implemented as electronic hardware, computer software, or a combination thereof. To clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described in terms of their functionality. Whether such functionality is implemented as hardware or software depends on the specific application and the design constraints imposed on the overall system. Those skilled in the art can implement the described functionality in alternative ways for each specific application; however, such implementation decisions should not be construed as departing from the scope of this disclosure.
[0175] The steps of the methods or algorithms described in conjunction with the embodiments herein can be directly embodied in hardware, software modules executed by a processor, or a combination thereof. The software modules can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disks, removable disks, CD-ROMs, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor, enabling the processor to read information from and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and storage medium can reside in an ASIC. The ASIC can reside in a user terminal. Alternatively, the processor and storage medium can exist as discrete components in the user terminal.
[0176] For software implementation, the techniques described in this application can be implemented using modules (e.g., procedures, functions, etc.) that perform the functions described in this application. This software code can be stored in memory units and executed by a processor. The memory units can be implemented within the processor or outside the processor; in the latter case, they are communicatively coupled to the processor via various means, as is well known in the art.
[0177] The foregoing description includes examples of one or more embodiments. It is certainly impossible to describe all possible combinations of components or methods in order to describe the above embodiments, but those skilled in the art will recognize that further combinations and arrangements of the various embodiments are possible. Therefore, the embodiments described herein are intended to cover all such changes, modifications, and variations that fall within the scope of the appended claims. Furthermore, the term "comprising" as used in the specification or claims is interpreted in a manner similar to the term "including," as interpreted when used as a conjunction in the claims. Additionally, the use of any term "or" in the specification of the claims is intended to mean "non-exclusive or."
Claims
1. A ray tracing region analysis method, characterized in that, include: Obtain a three-dimensional geological model of the area to be explored; The three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple symmetrical boxes are pre-constructed in each geological block. Each symmetrical box includes triangles in the triangular mesh corresponding to the symmetrical box. An initial ray is obtained based on the selected ray origin and the preset deflection angle relative to the formation normal at the ray origin; Determine the symmetrical box within the geological block where the ray originates and intersects with the initial ray; Based on the symmetrical boxes that intersect with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, resulting in a set of intersection points. If the initial ray intersects with a triangle in the triangulation of the intersecting symmetrical boxes, the intersection point zp of the initial ray and the triangle in the triangulation of the current geological block is determined. The intersection point zp satisfies the following condition: the value of (a1+a2+a3) / A is not greater than a preset area threshold, where A is the area of the triangle enclosed by vertices p1, p2, and p3, a1 is the area of the triangle enclosed by vertices p1, p2, and zp, a2 is the area of the triangle enclosed by vertices p2, p3, and zp, and a3 is the area of the triangle enclosed by vertices p1, p3, and zp. Select intersection points that meet the preset conditions from the set of intersection points, and use them as the path endpoint of the initial ray in the geological block where the ray originates, as well as the ray origin of the adjacent geological block to the geological block where the ray originates. If the endpoint of the path is not located on the surface of the three-dimensional geological model, the next ray is determined based on the stratum velocity of the two adjacent geological blocks and the ray origin of the adjacent geological blocks. The step of determining the symmetrical box in the current geological block where the ray origin is located that intersects with the initial ray is continued until the intersection point with the geological blocks including the surface in the three-dimensional geological model is obtained, and the endpoint of the ray is obtained.
2. The method as described in claim 1, characterized in that, The process of constructing multiple homosymmetric boxes within a geological block includes: Based on the shape parameters of the triangles included in the triangular mesh of the geological block, determine the bounding box shape parameters of the geological block; Based on the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or based on the shape parameters of the bounding box and the preset division ratio, the bounding box is divided into multiple initial symmetrical boxes. Determine whether the number of triangles included in the initial symmetrical box meets the preset number requirement; if not, adjust the shape parameters of the initial symmetrical box or adjust the division ratio, and then execute the step of dividing the bounding box; if yes, obtain multiple symmetrical boxes included in the geological block.
3. The method as described in claim 2, characterized in that, The triangle shape parameters include vertices; the bounding box shape parameters include the length, width, and height of the bounding box; the initial symmetrical box shape parameters include the length, width, and height of the initial symmetrical box.
4. The method as described in claim 2, characterized in that, Determine the number of triangles included in the initial homosymmetric box, including: Determine the total number of initial homosymmetric boxes for the geological blocks; Determine the total number of triangles in the triangulation of the geological block; The number of triangles in the triangulation of the geological blocks included in each initial homosymmetric box can be determined using the following formula: triPerCell = totalTri / cells, where triPerCell represents the number of triangles in the triangulation of the geological block included in each initial homosymmetric box, totalTri represents the total number of triangles in the triangulation of the geological block, and cells represents the total number of initial homosymmetric boxes. The adjustment of the shape parameters of the initial symmetrical box includes: Adjust the length, width, and height of the initial symmetrical box using the following formulas: L1=L / n, W1=W / n, H1=H / n, where L1, W1, and H1 are the adjusted length, width, and height of the initial symmetrical box, respectively; L, W, and H are the original length, width, and height of the initial symmetrical box; and n is the adjustment ratio.
5. The method according to claim 4, characterized in that, The total number of initial homosymmetric boxes for determining the geological blocks includes: The number of initial symmetrical boxes in the X direction of the geological block is determined based on the length of the bounding box and the length of each initial symmetrical box. The number of initial symmetrical boxes in the Y direction for a geological block is determined based on the width of the bounding box and the width of each initial symmetrical box. The number of initial symmetrical boxes in the Z direction for a geological block is determined based on the height of the bounding box and the height of each initial symmetrical box. Multiply the number of initial homosymmetric boxes in the X, Y, and Z directions for each geological block, and obtain the total number of initial homosymmetric boxes constructed from the geological block based on the result of the multiplication.
6. The method according to claim 1, characterized in that, Determining the symmetrical box within the geological block where the ray originates and intersects with the initial ray includes: Determine the geological block where the ray originates; Determine all homosymmetric boxes contained in the geological block where the ray originates, and establish a set of homosymmetric boxes; Determine the homosymmetric boxes that intersect the initial ray with the set of homosymmetric boxes, thus obtaining the homosymmetric boxes that intersect with the initial ray.
7. The method according to claim 6, characterized in that, The geological block where the ray origin is located includes: Based on the triangulation where the ray originates, obtain geological blocks B1 and B2 that are adjacent to the triangulation. The position coordinates of the ray starting point are finely adjusted according to the preset adjustment rules; Using the finely adjusted starting point of the ray as the endpoint, construct a vertical segment b1 that extends beyond geological block B1 and a vertical segment b2 that extends beyond geological block B2. Determine whether the number of intersections between perpendicular segment b1 and geological block B1 and the number of intersections between perpendicular segment b2 and geological block B2 are odd. The geological block with an odd number of intersections is the geological block where the ray originates.
8. The method as described in claim 1, characterized in that, Based on the symmetrical boxes that intersect with the initial ray, and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, resulting in a set of intersection points, including: Traverse the homosymmetric boxes that intersect with the initial ray; Determine whether the initial ray intersects with a triangle in the triangular mesh contained in an intersecting symmetrical box; If so, determine the intersection point of the initial ray with the triangle in the triangulation of the current geological block, and add the determined intersection point to the intersection point set; After traversing all the homogeneous symmetric boxes that intersect with the initial ray, we obtain the set of intersection points.
9. The method as described in claim 8, characterized in that, Determining the intersection points of the initial ray and the triangles in the triangulation of the current geological block includes: Calculate the side vectors v1 and v2 of the triangle based on its three vertices p1, p2, and p3. Determine the cross product xv of the edge vectors based on the edge vectors; The ray segment vector vr is determined based on the initial ray's starting point r1 and ending point r2. Determine vector v3 based on the vertex p1 of the triangle and the starting point r1 of the initial ray; The dot product d1 of the two is determined based on the cross product xv and the vector vr; The dot product d2 of vector v3 and the cross product xv is determined based on the vector v3 and the cross product xv. The point in the triangle is determined by the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle; Determine the area A of the triangle formed by vertices p1, p2, and p3, the area a1 of the triangle formed by vertices p1, p2, and zp, the area a2 of the triangle formed by vertices p2, p3, and zp, and the area a3 of the triangle formed by p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then point zp in the triangle is the intersection of the initial ray and the triangle in the triangular mesh of the current geological block.
10. The method according to claim 1, characterized in that, Also includes: Based on the endpoints of multiple initial rays that pass through the preset deflection angle of the ray origin, the imaging range of the ground surface included in the three-dimensional geological model is obtained by emitting rays with the ray origin as the starting point and the preset deflection angle as the direction.
11. The method according to claim 10, characterized in that, Multiple initial rays passing through the preset deflection angle of the ray origin include: Based on the selected ray starting point and the preset deflection angle of the stratum normal relative to the ray starting point, multiple initial rays are obtained by using the preset deflection angle of the stratum normal as the ray direction and distributing them evenly in a circle at preset intervals.
12. A ray tracing region analysis device, characterized in that, include: The data acquisition module is used to acquire a three-dimensional geological model of the area to be explored. The three-dimensional geological model includes multiple layers and multiple geological blocks. Each layer includes multiple triangular meshes. Each geological block is enclosed by triangular meshes from multiple layers. Each triangular mesh includes multiple triangles. Multiple symmetrical boxes are pre-constructed in each geological block. Each symmetrical box includes triangles in the triangular mesh that correspond to the symmetrical box. A ray determination module is used to obtain an initial ray based on a selected ray origin and a preset deflection angle relative to the formation normal at the ray origin. The next ray is determined based on the stratigraphic velocity of two adjacent geological blocks and the ray origin of the adjacent geological blocks; The intersection point determination module is used to determine the symmetrical boxes in the geological block where the ray origin is located that intersect with the initial ray; based on the symmetrical boxes that intersect with the initial ray and the corresponding triangles in the triangulation of the intersecting symmetrical boxes, the intersection points of the initial ray and the triangles in the triangulation of the current geological block are determined, and the intersection point set is obtained. Select intersection points that meet preset conditions from the set of intersection points as the path endpoint of the initial ray in the geological block where the ray originates, and as the ray origin of the adjacent geological block of the geological block where the ray originates; if the path endpoint is not located on the surface of the three-dimensional geological model, notify the ray determination module and, based on the next ray determined by the ray determination module, notify the intersection point determination module to continue executing the step of determining the symmetrical box that intersects with the initial ray in the current geological block where the ray originates, until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained, and the ray endpoint is obtained; If the initial ray intersects with a triangle in the triangular mesh contained in the intersecting symmetrical box, then the intersection point zp of the initial ray and the triangle in the triangular mesh of the current geological block is determined. The intersection point zp satisfies the following condition: the value of (a1+a2+a3) / A is not greater than the preset area threshold, where A is the area of the triangle enclosed by vertices p1, p2, and p3, a1 is the area of the triangle enclosed by vertices p1, p2, and zp, a2 is the area of the triangle enclosed by vertices p2, p3, and zp, and a3 is the area of the triangle enclosed by triangles p1, p3, and zp.
13. The apparatus according to claim 12, characterized in that, Also includes: The imaging range determination module is used to obtain the imaging range of the ground surface included in the three-dimensional geological model based on the ray endpoint determined by multiple initial rays passing through the ray origin.
14. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the ray tracing region analysis method according to any one of claims 1-11.
15. A terminal device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements a ray tracing region analysis method according to any one of claims 1-11.