Ray tracing area analysis method and device

Abstract

The invention discloses a ray tracing area analysis method and device. The method comprises the following steps: acquiring a three-dimensional geologic model of a to-be-explored area, establishing a description structure of the model, and constructing a symmetrical box; obtaining an initial ray based on the ray starting point and a preset deflection angle of the stratum normal at the relative ray starting point; determining an aligned symmetric box intersected with the initial ray; determining intersection points of the initial ray and triangles in the triangulation network in the intersected symmetrical box to obtain an intersection point set, selecting intersection points meeting preset conditions, continuing to execute the process of determining the intersection points of the ray and triangles in the triangulation network of the geological block until the intersection points of the ray and the geological block including the earth surface in the three-dimensional geological model are obtained, and obtaining a ray end point; and determining the imaging range of the ray starting point on the ground surface according to the ray end points. The calculation amount of solving the intersection point between the ray and the three-dimensional geologic model can be reduced, the processing time is shortened, the accuracy of solving the intersection point is effectively improved, and the efficiency of the actual exploration process is improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration, and particularly to a ray tracing area analysis method and device. Background Art

[0002] Seismic wave ray tracing is a seismic wave forward modeling technology that is currently most widely used in actual production. This technology is based on some kinematic basic theories of seismic wave propagation and computer technology. In actual application scenarios, the ray tracing forward modeling technology of a three-dimensional geological model can provide strong technical support for the optimal design of seismic observation schemes in complex areas. Ray tracing can not only obtain the travel time information of seismic waves, but also obtain ray trajectories, which has important theoretical significance and practical value for the research in the field of inverse tomography.

[0003] The ray tracing analysis method of a geological model is based on the technology of finding the intersection points of rays and a three-dimensional geological model. By obtaining the intersection points of rays in each formation of the three-dimensional complex geological model and the ray endpoints on the surface of the geological model of the ray segments, the paths of rays emitted from different points in the formation and the imaging range on the surface are determined, providing technical support for the exploration of the three-dimensional complex geological model. Summary of the Invention

[0004] The present invention discovers that when using the conventional technology of finding the intersection points of rays and a three-dimensional geological model for seismic wave ray tracing, when the number of triangles in the three-dimensional geological model layer is large or the number of rays is huge, the time consumption of ray tracing is high, the calculation efficiency is poor, and it is difficult to meet the requirements of real-time analysis of ray tracing based on the geological model. How to quickly and accurately obtain the intersection points of rays in the three-dimensional geological model and improve the 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 area analysis method and device that overcome the above problems or at least partially solve the above problems.

[0006] In a first aspect, an embodiment of the present invention provides a ray tracing area analysis method, including:

[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 of multiple layers to form a closed body, each triangular mesh includes multiple triangles, and multiple quasi-symmetric boxes are pre-constructed in each geological block, and each quasi-symmetric box includes the triangles in the triangular mesh corresponding to the quasi-symmetric box;

[0008] Based on the selected ray starting point and a preset deviation angle relative to the formation normal line at the ray starting point, obtain an initial ray;

[0009] Determine the symmetrical box in the geological block where the ray starts that intersects the initial ray;

[0010] According to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulated network included in the intersecting symmetrical box, the intersection point of the initial ray and the triangle in the triangulated network of the current geological block is determined to obtain an intersection point set;

[0011] Selecting an intersection point that meets a preset condition from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located;

[0012] If the end point of the path is not located on the surface of the three-dimensional geological model, the next ray is determined based on the formation velocity of the two adjacent geological blocks and the starting point of the ray in the adjacent geological block, and the step of determining the symmetrical box intersecting with the initial ray in the current geological block where the ray starting point is located is continued until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the end point of the ray is obtained.

[0013] In some optional embodiments, the process of constructing a plurality of symmetrical boxes in a geological block includes:

[0014] Determine the bounding box shape parameters of the geological block according to the triangle shape parameters included in the triangulated network in the geological block;

[0015] According to the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or according to the shape parameters of the bounding box and the preset division ratio, the bounding box is divided into a plurality of 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 perform the step of dividing the bounding box; if so, obtain multiple symmetrical boxes included in the geological block.

[0017] In some optional 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 aligned symmetric box shape parameters include the length, width and height of the bounding box.

[0018] In some optional embodiments, determining the number of triangles included in the initial symmetrical box includes:

[0019] Determine the total number of initial symmetric boxes of the geological block;

[0020] Determine the total number of triangles in the triangulation network of the geological block;

[0021] The number of triangles in the triangulation of the geological block included in each initial symmetric box is determined using the following formula:

[0022] triPerCell=totalTri / cells, where triPerCell represents the number of triangles in the triangulation network of the geological block included in each initial symmetrical box, totalTri represents the total number of triangles in the triangulation network of the geological block, and cells represents the total number of initial symmetrical boxes;

[0023] Adjust the shape parameters of the initial symmetrical box, including:

[0024] Adjust the length, width, and height of the initial symmetrical box by the following formula: L 1 =L / n,W 1 =W / n,H 1 =H / n, where L 1 , W 1 , H 1 are the length, width and height of the initial symmetrical box after adjustment, L, W and H are the length, width and height of the initial symmetrical box before adjustment, and n is the adjustment ratio value.

[0025] In some optional embodiments, determining the total number of initial symmetrical boxes of the geological block includes:

[0026] Determine the number of initial symmetrical boxes of the geological block in the X direction according to the length of the bounding box and the length of each initial symmetrical box;

[0027] Determine the number of initial symmetrical boxes of the geological block in the Y direction according to the width of the bounding box and the width of each initial symmetrical box;

[0028] Determine the number of initial symmetrical boxes of the geological block in the Z direction according to the height of the enclosing sum and the height of each initial symmetrical box;

[0029] The number of initial aligned symmetric boxes of each geological block in the X, Y, and Z directions is multiplied, and the total number of multiple initial aligned symmetric boxes constructed by the geological block is obtained according to the multiplication result.

[0030] In some optional embodiments, determining an symmetrical box in a geological block where a ray starting point is located and intersecting with the initial ray comprises:

[0031] Determine the geological block where the ray starts;

[0032] Determine all the symmetrical boxes contained in the geological block where the ray starting point is located, and establish a set of symmetrical boxes;

[0033] Determine the symmetric boxes in the set of symmetric boxes that intersect with the initial ray to obtain the symmetric boxes that intersect with the initial ray.

[0034] In some alternative embodiments, determining the geological block where the ray starting point is located includes:

[0035] According to the triangulation where the ray starting point is located, obtain geological blocks B1 and B2 adjacent to the triangulation;

[0036] Fine-tune the position coordinates of the ray starting point according to a preset adjustment rule;

[0037] Taking the position of the ray starting point after fine-tuning as an endpoint, construct a perpendicular segment b1 that extends beyond geological block B1 and a perpendicular segment b2 that extends beyond geological block B2;

[0038] Judge whether the number of intersection points between the perpendicular segment b1 and geological block B1 and the number of intersection points between the perpendicular segment b2 and geological block B2 are odd, and determine that the geological block with an odd number of intersection points is the geological block where the ray starting point is located.

[0039] In some alternative embodiments, according to the symmetric boxes that intersect with the initial ray and the corresponding triangles in the triangulation included in the intersecting symmetric boxes, determine the intersection points of the initial ray and the triangles in the triangulation of the current geological block to obtain an intersection point set, including:

[0040] Traverse the symmetric boxes that intersect with the initial ray;

[0041] Judge whether the initial ray intersects with the triangles in the triangulation included in the intersecting symmetric boxes;

[0042] If so, determine the intersection points of the initial ray and the triangles in the triangulation of the current geological block, and add the determined intersection points to the intersection point set;

[0043] After traversing the symmetric boxes that intersect with the initial ray, obtain the intersection point set.

[0044] In some alternative embodiments, determining the intersection points of the initial ray and the triangles in the triangulation of the current geological block includes:

[0045] According to the three vertices p1, p2, p3 of the triangle, calculate the side vectors v1 and v2 of the triangle;

[0046] Determine the cross product xv of the side vectors according to the side vectors;

[0047] Determine the ray segment vector vr according to the starting point r1 and the ending point r2 of the initial ray;

[0048] Determine the vector v3 according to the vertex p1 of the triangle and the starting point r1 of the initial ray;

[0049] Determine the dot product d1 of the two based on the cross product xv and the vector vr;

[0050] Determine the dot product d2 of the vector v3 and the cross product xv;

[0051] Determine the point in the triangle according to the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle;

[0052] Determine the area A of the triangle enclosed by vertices p1, p2, and p3, the area a1 of the triangle enclosed by vertices p1, p2, and zp, the area a2 of the triangle enclosed by vertices p2, p3, and zp, and the area a3 of the triangle p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then the point zp in the triangle is the intersection of the initial ray and the triangle in the triangulation of the current geological block.

[0053] In some optional embodiments, the above method further includes:

[0054] Based on the ray endpoints determined by multiple initial rays passing through the ray starting points with preset deflection angles, an imaging range of the ground surface included in the three-dimensional geological model of the rays emitted with the ray starting points as the starting points and with the preset deflection angles as the directions is obtained.

[0055] In some optional embodiments, the multiple initial rays passing through the ray starting point with a preset deflection angle include:

[0056] Based on the selected ray starting point and the preset deflection angle of the formation normal relative to the ray starting point, a plurality of initial rays are obtained by taking the preset deflection angle of the formation normal as the ray direction and uniformly distributing the rays around the circumference at preset interval angles.

[0057] In a second aspect, an embodiment of the present invention provides a ray tracing region analysis device, comprising:

[0058] A 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 levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and each geological block has multiple symmetrical boxes pre-constructed, each symmetrical box includes triangles in the triangulated network corresponding to the symmetrical box;

[0059] A ray determination module is used to obtain an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point; and to determine a next ray based on stratum velocities of two adjacent geological blocks and ray starting points of adjacent geological blocks;

[0060] The intersection determination module is used to determine the symmetrical box intersecting with the initial ray in the geological block where the ray starting point is located; determine the intersection of the initial ray and the triangle in the triangulation network of the current geological block according to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulation network included in the intersecting symmetrical box, and obtain an intersection set; select the intersection that meets the preset conditions from the intersection set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located; if the path end point 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 determination module, and continue to execute the step of determining the symmetrical box intersecting with the initial ray in the current geological block where the ray starting point is located, until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the ray end point is obtained.

[0061] In some optional embodiments, the above device 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 of the ray emitted from the ray starting point with a preset deflection angle as the starting point based on the ray end point determined by multiple initial rays passing through the ray starting point.

[0063] An embodiment of the present invention further provides a computer storage medium, in which computer executable instructions are stored. When the computer executable instructions are executed by a processor, a ray tracing area analysis method is implemented.

[0064] An embodiment of the present invention further provides a terminal device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements a ray tracing area analysis method when executing the program.

[0065] The beneficial effects of the above technical solution provided by the embodiment of the present invention include at least:

[0066] The method provided by the embodiment of the present invention divides the three-dimensional geological model of the exploration area into triangular face elements, constructs an symmetrical box based on the face element division, constructs an initial ray based on the selected ray starting point and the preset deflection angle relative to the stratum normal at the ray starting point, performs intersection based on the symmetrical box, determines the intersection point of the initial ray with the triangle in the triangulation network of the current geological block, and obtains an intersection point set without performing intersection calculation for the entire geological block, thereby greatly reducing the amount of calculation in the intersection process and improving the processing speed and efficiency; selects an intersection point that meets the preset conditions from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, constructs the next ray with the path end point of the current geological block as the ray starting point of the next adjacent geological block, and continues to perform the intersection processing of the next adjacent geological block based on the symmetrical box included in the adjacent geological block until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained, thereby obtaining the ray end point. When obtaining the intersection points of rays with various strata and geological blocks, the method first determines the symmetrical boxes that the rays pass through and then calculates the intersection points. There is no need to perform intersection calculations for the symmetrical boxes that the rays do not pass through, which greatly reduces the amount of data processing 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 analysis of ray tracing based on geological models in the prior art, improves the accuracy and analysis efficiency of ray tracing, and provides strong technical support for seismic wave ray tracing analysis.

[0067] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings.

[0068] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0070] Figure 1 Flow chart of the ray tracing regional analysis method in Embodiment 1 of the present invention

[0071] Figure 2 is a flow chart of the ray tracing regional analysis method in Embodiment 2 of the present invention;

[0072] Figure 3 This is a schematic diagram of a three-dimensional geological model in the second embodiment of the present invention;

[0073] Figure 4 A schematic diagram of geological blocks in a three-dimensional geological model in the second embodiment of the present invention;

[0074] Figure 5 It is a schematic diagram of a symmetrical box of a geological block in the second embodiment of the present invention;

[0075] Figure 6 A schematic diagram of determining the intersection point of an initial ray and a triangle in Embodiment 2 of the present invention;

[0076] Figure 7 This is a schematic diagram of ray tracing in the second embodiment of the present invention;

[0077] Figure 8 This is a schematic diagram of the imaging range of an initial ray from a ray starting point on the ground surface in the second embodiment of the present invention;

[0078] Fig. 9 Schematic diagram of the imaging range of the initial rays of multiple ray starting points on the ground surface in the second embodiment of the present invention;

[0079] Fig.10 It is a structural diagram of the ray tracing regional analysis device in the second embodiment of the present invention. DETAILED DESCRIPTION

[0080] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0081] Ray tracing technology occupies a very important position in seismic exploration. It is widely used in other seismic digital processing and inversion such as seismic positioning and geophysical tomography. The speed and accuracy of ray tracing are the key to solving practical production problems. However, the ray tracing technology used in the prior art has the problems of high tracing time and poor calculation efficiency when there are a large number of rays. It is difficult to meet the needs of real-time analysis of ray tracing based on geological models. In order to solve the problems existing in the prior art, an embodiment of the present invention provides a ray tracing regional analysis method, which can solve the defects in the prior art and provide strong technical support for the exploration of complex geological areas.

[0082] Embodiment 1

[0083] Embodiment 1 of the present invention provides a ray tracing regional analysis method, the process of which is as follows: Figure 1 As shown, the following steps are included:

[0084] S101: Acquire a three-dimensional geological model of the area to be explored; the three-dimensional geological model includes multiple levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and multiple symmetrical boxes are pre-constructed in each geological block, each symmetrical box includes triangles in the triangulated network corresponding to the symmetrical box;

[0085] S102: obtaining an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point;

[0086] S103: determining an symmetrical box in the geological block where the starting point of the ray is located and intersecting with the initial ray;

[0087] S104: determining the intersection points of the initial ray and the triangles in the triangulation network of the current geological block according to the symmetrical boxes intersecting with the initial ray and the corresponding triangles in the triangulation network included in the symmetrical boxes intersecting with the initial ray, and obtaining an intersection point set;

[0088] S105: Selecting an intersection point that meets a preset condition from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located;

[0089] S106: If the end point of the path is not located on the surface of the three-dimensional geological model, the next ray is determined based on the formation velocity of the two adjacent geological blocks and the starting point of the ray in the adjacent geological block, and the step of determining the symmetrical box that intersects with the initial ray in the current geological block where the ray starting point is located is continued until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the end point of the ray is obtained.

[0090] The method provided by the embodiment of the present invention divides the three-dimensional geological model of the exploration area into triangular face elements, constructs an symmetrical box based on the face element division, constructs an initial ray based on the selected ray starting point and the preset deflection angle relative to the stratum normal at the ray starting point, performs intersection based on the symmetrical box, determines the intersection point of the initial ray with the triangle in the triangulation network of the current geological block, and obtains an intersection point set without performing intersection calculation for the entire geological block, thereby greatly reducing the amount of calculation in the intersection process and improving the processing speed and efficiency; selects an intersection point that meets the preset conditions from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, constructs the next ray with the path end point of the current geological block as the ray starting point of the next adjacent geological block, and continues to perform the intersection processing of the next adjacent geological block based on the symmetrical box included in the adjacent geological block until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained, thereby obtaining the ray end point. When obtaining the intersection points of rays with various strata and geological blocks, the method first determines the symmetrical boxes that the rays pass through and then calculates the intersection points. There is no need to perform intersection calculations for the symmetrical boxes that the rays do not pass through, which greatly reduces the amount of data processing 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 analysis of ray tracing based on geological models in the prior art, improves the accuracy and analysis efficiency of ray tracing, and provides strong technical support for seismic wave ray tracing analysis.

[0091] Embodiment 2

[0092] Embodiment 2 of the present invention provides a ray tracing area analysis method. After determining the ray end point 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 as follows: Figure 2 As shown, the following steps are included:

[0093] Step S201: Acquire a three-dimensional geological model of the area to be explored.

[0094] Step S202: obtaining an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point.

[0095] Step S203: Determine an symmetrical box in the geological block where the starting point of the ray is located that intersects with the initial ray.

[0096] Step S204: According to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulation network included in the symmetrical box intersecting with the initial ray, the intersection points of the initial ray and the triangles in the triangulation network of the current geological block are determined to obtain a set of intersection points.

[0097] Step S205: Selecting intersection points that meet preset conditions from the intersection point set as the path end points of the initial ray in the geological block where the ray starting point is located, and the ray starting points of the geological blocks adjacent to the geological block where the ray starting point is located.

[0098] Step S206: Determine whether the end point of the path is located on the surface of the three-dimensional geological model. If not, execute step S207; if so, execute step S208.

[0099] Step S207: Based on the formation velocities of two adjacent geological blocks and the ray starting points of the adjacent geological blocks, determine the next ray and continue to execute step S203.

[0100] Step S208: taking the end point of the path located on the surface of the three-dimensional geological model as the end point of the ray.

[0101] Step S209: Based on the ray endpoints determined by multiple initial rays passing through the ray starting point with a preset deflection angle, an imaging range of the ground surface included in the three-dimensional geological model of the rays emitted from the ray starting point with the preset deflection angle as the direction is obtained.

[0102] The specific implementation process of each step in the above-mentioned embodiment 1 and embodiment 2 is described in detail below.

[0103] Preferably, in the above steps S101 and S201, the three-dimensional geological model of the area to be explored is described using the structure of "volume -> block -> surface -> triangulated network -> triangle", and the physical interfaces between different geological blocks are composed of triangulated networks, wherein the three-dimensional geological model includes multiple levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and multiple symmetrical boxes are pre-constructed in each geological block, each symmetrical box includes a triangle in the triangulated network corresponding to the symmetrical box.

[0104] Figure 3 This is a schematic diagram of the three-dimensional geological model. Figure 3 The three-dimensional geological model in contains six strata, and the closed area surrounded by the strata is called a geological block. Figure 4 This is a schematic diagram of geological blocks in the three-dimensional geological model. Figure 4 The top surface of the geological block in the figure is composed of 7 triangulated meshes, and the side and front surfaces are composed of one triangulated mesh respectively. Correspondingly, the surfaces in the other three directions of the geological block are also composed of triangulated meshes, each of which contains several triangular facets. The regular cube outside the geological block is the bounding box of the geological block.

[0105] Preferably, a plurality of symmetrical boxes are constructed in the geological block, including:

[0106] Determine the bounding box shape parameters of the geological block according to the triangle shape parameters included in the triangulated network in the geological block;

[0107] According to the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or according to the shape parameters of the bounding box and the preset division ratio, the bounding box is divided into a plurality of 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 perform the step of dividing the bounding box; if so, obtain multiple symmetrical boxes included in the geological block.

[0109] Constructing multiple symmetrical boxes for the geological block can greatly reduce the amount of calculation. There is no need to perform intersection calculations on the triangles in the entire geological block. Instead, the intersection calculations are performed on the triangles contained in the symmetrical boxes that intersect with the initial rays. This can improve the calculation efficiency and accuracy when encountering intersection points with geological blocks with a large number of triangles.

[0110] Optionally, the shape parameters of the triangle and the shape parameters of the bounding box can be selected as needed, and the above-mentioned triangle shape parameters may include vertices; the above-mentioned bounding box shape parameters may include the length, width and height of the bounding box; and the above-mentioned 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, and after the 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 used as the final symmetrical box.

[0111] Figure 5 This is a schematic diagram of a symmetrical box of a geological block. Figure 5 The bounding box of the geological block is divided into several initial symmetrical boxes according to the division rules. Each symmetrical box includes several triangles. Some symmetrical boxes may contain triangles from multiple triangulated networks, such as Figure 5 There are four initially aligned symmetrical boxes in the first column on the left. Each box includes multiple triangles in the triangulated network. The rays in the figure pass through some of the initially aligned symmetrical boxes.

[0112] Preferably, determining the number of triangles included in the initial symmetrical box comprises:

[0113] Determine the total number of initial symmetric boxes of the geological block;

[0114] Determine the total number of triangles in the triangulation network of the geological block;

[0115] The number of triangles in the triangulation of the geological block included in each initial symmetric box is determined using the following formula:

[0116] triPerCell=totalTri / cells, where triPerCell represents the number of triangles in the triangulation network of the geological block included in each initial symmetrical box, totalTri represents the total number of triangles in the triangulation network of the geological block, and cells represents the total number of initial symmetrical 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 by the following formula: 1 =L / n,W 1 =W / n,H 1 =H / n, where L 1 , W 1 , H 1 are the length, width and height of the adjusted initial symmetrical box, L, W and H are the length, width and height of the initial symmetrical box before adjustment, and n is the adjustment ratio value. For example, in this embodiment, the initial symmetrical box is adjusted with an adjustment ratio of 2, and the ratio parameter can be selected according to the fine division requirements. The shape parameters of the initial symmetrical box can be selected as needed and adjusted according to the preset adjustment rules.

[0119] Preferably, determining the total number of initial symmetrical boxes of the geological block comprises:

[0120] Determine the number of initial symmetrical boxes of the geological block in the X direction according to the length of the bounding box and the length of each initial symmetrical box;

[0121] Determine the number of initial symmetrical boxes of the geological block in the Y direction according to the width of the bounding box and the width of each initial symmetrical box;

[0122] Determine the number of initial symmetrical boxes of the geological block in the Z direction according to the height of the enclosing sum and the height of each initial symmetrical box;

[0123] Multiply the number of initial aligned symmetric boxes of each geological block in the X, Y, and Z directions, and obtain the total number of multiple initial aligned symmetric boxes constructed by the geological block based on the multiplication result. Figure 5 As shown, the initial symmetrical box of the geological block is 7*4*1, and a total of 28 initial symmetrical boxes are divided. Accordingly, the total number of symmetrical boxes can be calculated according to the above method.

[0124] Preferably, in the above steps S102 and S202, the initial ray is directed at a preset deflection angle, starting from the selected ray starting point, and attempts to emit the initial ray to the stratum surface of the three-dimensional geological model. It should be noted that the ray in the present invention is the propagation path of the seismic wave moving in a medium with different velocities, and the propagation path is represented by a virtual ray.

[0125] Preferably, the above step S103 and step S203 include:

[0126] Determine the geological block where the ray starts;

[0127] Determine all the symmetrical boxes contained in the geological block where the ray starting point is located, and establish a set of symmetrical boxes;

[0128] Determine the symmetrical box that intersects the initial ray with the symmetrical box set, and obtain the symmetrical box that intersects the initial ray.

[0129] See also Figure 5 As shown, Figure 5 The ray in is the initial ray. According to the set of symmetrical boxes in the geological block, four symmetrical boxes in the set of symmetrical boxes that intersect with the initial ray are determined.

[0130] Preferably, determining the geological block where the ray starting point is located includes:

[0131] According to the triangulated network where the starting point of the ray is located, the geological blocks B1 and B2 adjacent to the triangulated network are obtained;

[0132] Fine-tune the position coordinates of the ray starting point according to the preset adjustment rules;

[0133] Taking the finely adjusted ray starting point position as the endpoint, construct the vertical segment b1 beyond the geological block B1 and the vertical segment b2 beyond the geological block B2;

[0134] It is determined whether the number of intersections between the vertical line segment b1 and the geological block B1 and the number of intersections between the vertical line segment b2 and the geological block B2 are odd numbers, and the geological block with an odd number of intersections is determined to be the geological block where the ray starting point is located.

[0135] In this embodiment, the adjustment rule may be to adjust the coordinate of the ray starting point in the Z direction, for example, to adjust the coordinate of the ray starting point in the Z direction upward by 0.25. Other adjustment rules may be designed according to actual exploration requirements.

[0136] Preferably, the above step S104 and step S204 include:

[0137] Traverse the flush symmetric box that intersects the initial ray;

[0138] Determine whether the initial ray intersects with the triangle in the triangulated network contained in the intersecting symmetrical box;

[0139] If so, determine the intersection of the initial ray and the triangle in the triangulation network of the current geological block, and add the determined intersection to the intersection set;

[0140] After traversing all the symmetrical boxes that intersect with the initial ray, we get the set of intersection points.

[0141] In flat areas, the initial ray and the triangles in the triangulated network of the symmetrical box have an intersection point. When encountering concave or uplifted strata, wrinkled areas will form between the layers of the strata, see Figure 4 As shown, Figure 4 The top layer is not a smooth layer, but has ridges and depressions. At this time, wrinkle areas will be formed between the layers. The symmetrical box containing the wrinkle areas will contain multiple triangles of the triangulated network. Therefore, there will be multiple intersections between the initial rays and the triangles. The intersections are added to the set to form an intersection set.

[0142] Preferably, determining the intersection of the initial ray with a triangle in the triangulation network of the current geological block comprises:

[0143] According to the three vertices p1, p2, and p3 of the triangle, calculate the edge vectors v1 and v2 of the triangle;

[0144] Determine the cross product xv of the edge vector based on the edge vector;

[0145] Determine the ray segment vector vr 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] Determine the dot product d1 of the two based on the cross product xv and the vector vr;

[0148] Determine the dot product d2 of the vector v3 and the cross product xv;

[0149] Determine the point in the triangle according to the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle;

[0150] Determine the area A of the triangle enclosed by vertices p1, p2, and p3, the area a1 of the triangle enclosed by vertices p1, p2, and zp, the area a2 of the triangle enclosed by vertices p2, p3, and zp, and the area a3 of the triangle p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then the point zp in the triangle is the intersection of the initial ray and the triangle in the triangulation 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 of the initial ray and the triangle in the triangulation of the current geological block, and the preset area threshold can be adjusted to meet the actual exploration needs. Figure 6 As shown, the various parameters that need to be used have been marked in the figure.

[0152] Preferably, in the above steps S105 and S205, it is necessary to select an intersection point that meets the preset conditions from the intersection point set. For example, in this embodiment, the point closest to the ray starting point and the intersection point set is selected as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located. Figure 4 As shown, for the stratigraphic fold area, there are multiple intersections between the ray and the triangulated network in the geological block. It is necessary to select the intersection that meets the requirements. At this time, you can set the selection rules, such as selecting the intersection closest to the starting point of the ray. Of course, you can also select intersections that meet other conditions according to other rules.

[0153] Here's how to determine the point in the set of intersection points that is closest to the starting point of the ray:

[0154] Calculate the distances between the starting point of the ray and the intersection points in the intersection point set in sequence, and find the minimum distance among them. The intersection point corresponding to the minimum distance is the intersection point that meets the conditions.

[0155] Optionally, in the above steps S106 and S206, after the end point of the path is determined, the determined end point of the path is judged to determine whether it is located on the surface of the three-dimensional geological model. If not, the next ray is determined based on the formation velocity of the two adjacent geological blocks and the ray starting point of the adjacent geological block. Step S106 continues to execute step S103, and step S206 continues to execute S203 until the determined end point of the path is located on the surface of the three-dimensional geological model. The path end point is used as the ray end point to end the process of obtaining the intersection point.

[0156] Since the medium in each stratum is non-uniform and the stratum velocities are different, according to the principle of wave transmission, the initial ray will be transmitted when passing through each stratum. Therefore, the direction of the initial ray propagation to the next stratum should be corrected according to the stratum velocities of the two adjacent geological blocks, and the intersection of the ray with the next stratum should be traced again. The steps of determining the symmetrical box that intersects with the initial ray in the geological block where the ray starting point is located are repeated until the intersection that is close enough to or reaches the ground surface of the geological model is traced. By fitting the intersection points in each stratum, a complete transmission ray can be obtained, thereby tracing the ray path of the initial ray in the stratum.

[0157] Figure 7It is a schematic diagram of ray tracing. Figure 7 Each point in it represents the intersection point of the initial ray in the formation. In the three-dimensional geological model of this embodiment, there are 6 formations, so 5 processes of finding intersection points are required. A ray segment is formed between every two intersection points. Since the medium of the formation is non-uniform and the formation velocities are different, the direction of the initial ray needs to be continuously corrected. Finally, all intersection points and ray segments are integrated to obtain a transmitted ray.

[0158] Preferably, in the above step S209, multiple initial rays with a preset deflection angle at the ray starting point include:

[0159] Based on the selected ray starting point and the preset deflection angle relative to the formation normal at the ray starting point, multiple initial rays are obtained by circularly distributing them at equal intervals with the preset deflection angle of the formation normal as the ray direction.

[0160] Figure 8 It is a schematic diagram of the imaging range of the initial rays of a ray starting point on the ground surface. From Figure 8 it can be seen that with the selected ray starting point as the starting point of multiple initial rays, multiple are emitted according to the preset deflection angle. These initial rays are circularly distributed at equal intervals with the preset deflection angle of the formation normal as the ray direction. Through the ray region tracking and analysis method of the present invention, the intersection point of each initial ray in each formation can be determined; the transmitted ray can be determined according to the intersection points; according to the ray end points determined by multiple initial rays, the imaging range of the selected ray starting point on the ground surface included in the three-dimensional geological model can be obtained, thereby providing efficient and accurate technical support for seismic exploration in complex regions.

[0161] Fig. 9 It is a schematic diagram of the imaging range of the initial rays of multiple selected ray starting points on the ground surface. In Fig. 9 it can be seen that not necessarily all initial rays will have intersection points with the ground surface. At the boundary position of the geological block, there will be some situations where the initial rays cannot reach the ground surface. In actual exploration, only the ray end points that can reach the ground surface need to be obtained, so as to determine the imaging range of multiple initial rays emitted from the ray starting point on the ground surface.

[0162] Based on the same inventive concept, the embodiment of the present invention also provides a ray tracing region analysis device. This device can be set in a device with the ability to process computer instructions. The structure of this device is as Fig.10 shown and 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 levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and each geological block has multiple symmetrical boxes pre-constructed, each symmetrical box includes triangles in the triangulated network corresponding to the symmetrical box;

[0164] The ray determination module 12 is used to obtain an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point; determine a next ray based on stratum velocities of two adjacent geological blocks and ray starting points of adjacent geological blocks;

[0165] The intersection determination module 13 is used to determine the symmetrical box intersecting with the initial ray in the geological block where the ray starting point is located; determine the intersection of the initial ray and the triangle in the triangulation network of the current geological block according to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulation network included in the intersecting symmetrical box, and obtain an intersection set; select the intersection that meets the preset conditions from the intersection set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located; if the path end point 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 determination module, and continue to execute the step of determining the symmetrical box intersecting with the initial ray in the current geological block where the ray starting point is located, until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the ray end point is obtained.

[0166] Preferably, the above device further comprises:

[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 of the ray emitted from the ray starting point with a preset deflection angle based on the ray end point determined by multiple initial rays passing through the ray starting point.

[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 deviation angle of the formation normal relative to the ray starting point, and the multiple initial rays are evenly distributed around the circle at a preset interval angle; based on the formation velocity of two adjacent geological blocks and the ray starting point of the adjacent geological blocks, determine the next ray of the multiple initial rays propagating between different formations.

[0169] Regarding the ray tracing area analysis device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0170] The above method and device of the embodiment of the present invention obtains a three-dimensional geological model of the area to be explored and establishes a description structure to construct an symmetrical box of the geological block of the geological model; obtains an initial ray based on a ray starting point selected in the stratum of the geological model and a preset deflection angle relative to the stratum normal at the ray starting point; determines an symmetrical box in the geological block where the ray starting point is located that intersects with the initial ray; determines the intersection point between the initial ray and a triangle in the triangulation network of the current geological block to obtain an intersection point set; selects an intersection point that meets preset conditions from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located and the ray starting point of the next adjacent geological block; if the path end point is not located on the surface of the three-dimensional geological model, determines the next ray based on the stratum velocity of the two adjacent geological blocks and the ray starting point of the adjacent geological block, and continues to perform the step of obtaining the stratum intersection point until the intersection point with the geological block including the surface in the three-dimensional geological model is obtained to obtain the ray end point; determines the imaging range of the selected ray starting point on the ground surface according to the ray end point, and determines the ray path of the initial ray in the three-dimensional geological model according to the intersection point of the initial ray with each stratum. The above method and device overcome the problem of high time consumption of ray tracing in the prior art, can effectively improve the efficiency of ray tracing analysis based on three-dimensional geological models, improve the accuracy of ray tracing, and provide technical support for the exploration of complex areas.

[0171] Unless otherwise specifically stated, terms such as processing, computing, calculating, determining, displaying, etc. may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, which operate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may 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 may be represented by voltages, currents, 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 can be rearranged without departing from the scope of protection of the present disclosure. The attached method claims present the elements of the various steps in an exemplary order and are not intended to be limited to the specific order or hierarchy described.

[0173] In the above detailed description, various features are grouped together in a single embodiment to simplify the disclosure. This method of disclosure should not be interpreted as reflecting an intention that the embodiments of the claimed subject matter require more features than are clearly stated in each claim. On the contrary, as reflected in the appended claims, the invention is in a state of having less than all the features of the disclosed individual embodiments. Therefore, the appended claims are hereby expressly incorporated into the detailed description, with each claim standing on its own as a separate preferred embodiment of the invention.

[0174] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithmic steps described in conjunction with the embodiments herein can all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, various illustrative components, blocks, modules, circuits, and steps are generally described above around their functions. Whether such functions are implemented as hardware or software depends on specific applications and the design constraints imposed on the entire system. A skilled person can implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present disclosure.

[0175] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0176] For software implementation, the techniques described in this application can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in this application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0177] The above description includes examples of one or more embodiments. Of course, it is impossible to describe all possible combinations of components or methods for the purpose of describing the above embodiments, but it should be recognized by those skilled in the art that the various embodiments may be further combined and arranged. Therefore, the embodiments described herein are intended to cover all such changes, modifications and variations that fall within the scope of protection of the appended claims. In addition, with respect to the term "comprising" used in the specification or claims, the word is covered in a manner similar to the term "including", just as "including," is explained as a transitional word in the claims. In addition, any term "or" used in the specification of the claims is intended to mean "non-exclusive or".

Claims

1. A ray tracing regional analysis method, It is characterized in that include: Obtain a three-dimensional geological model of the area to be explored; The three-dimensional geological model includes multiple levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by the triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and multiple symmetrical boxes are pre-constructed in each geological block, each symmetrical box includes triangles in the triangulated network corresponding to the symmetrical box; Obtaining an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point; Determine an symmetrical box in the geological block where the starting point of the ray is located and intersects with the initial ray; According to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulated network included in the symmetrical box intersecting with the initial ray, the intersection point of the initial ray and the triangle in the triangulated network of the current geological block is determined to obtain an intersection point set; Selecting an intersection point that meets a preset condition from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located; If the end point of the path is not located on the surface of the three-dimensional geological model, the next ray is determined based on the formation velocity of the two adjacent geological blocks and the starting point of the ray in the adjacent geological block, and the step of determining the symmetrical box intersecting with the initial ray in the current geological block where the starting point of the ray is located is continued until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the end point of the ray is obtained.

2. The method according to claim 1, It is characterized in that The process of constructing multiple symmetrical boxes in a geological block includes: Determine the bounding box shape parameters of the geological block according to the triangle shape parameters included in the triangulated network in the geological block; According to the shape parameters of the bounding box and the predetermined shape parameters of the initial symmetrical box, or according to the shape parameters of the bounding box and a preset division ratio, the bounding box is divided into a plurality of 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 perform the step of dividing the bounding box; if so, obtain multiple symmetrical boxes included in the geological block.

3. The method according to claim 2, It is characterized in that The triangle shape parameters include vertices; the bounding box shape parameters include the length, width and height of the bounding box; and the initial symmetrical box shape parameters include the length, width and height of the bounding box.

4. The method according to claim 2, It is characterized in that Determine the number of triangles included in the initial symmetrical box, including: Determine the total number of initial symmetric boxes of the geological block; Determine the total number of triangles in the triangulation network of the geological block; The number of triangles in the triangulation of the geological block included in each initial symmetric box is determined using the following formula: triPerCell=totalTri / cells, where triPerCell represents the number of triangles in the triangulation network of the geological block included in each initial symmetrical box, totalTri represents the total number of triangles in the triangulation network of the geological block, and cells represents the total number of initial symmetrical boxes; The adjusting of the shape parameters of the initial symmetrical box includes: Adjust the length, width, and height of the initial symmetrical box by the following formula: L 1 =L / n,W 1 =W / n,H 1 =H / n, where L 1 , W 1 , H 1 are the length, width and height of the initial symmetrical box after adjustment, L, W and H are the length, width and height of the initial symmetrical box before adjustment, and n is the adjustment ratio value.

5. The method according to claim 4, It is characterized in that The total number of initial symmetrical boxes of the geological block is determined, including: Determine the number of initial symmetrical boxes of the geological block in the X direction according to the length of the bounding box and the length of each initial symmetrical box; Determine the number of initial symmetrical boxes of the geological block in the Y direction according to the width of the bounding box and the width of each initial symmetrical box; Determine the number of initial symmetrical boxes of the geological block in the Z direction according to the height of the bounding box and the height of each initial symmetrical box; The number of initial aligned symmetric boxes of each geological block in the X, Y, and Z directions is multiplied, and the total number of multiple initial aligned symmetric boxes constructed by the geological block is obtained according to the multiplication result.

6. The method according to claim 1, It is characterized in that Determining a symmetrical box in the geological block where the starting point of the ray is located and intersecting with the initial ray comprises: Determine the geological block where the ray starts; Determine all the symmetrical boxes contained in the geological block where the ray starting point is located, and establish a set of symmetrical boxes; An symmetrical box intersecting the initial ray with the symmetrical box set is determined to obtain an symmetrical box intersecting the initial ray.

7. The method according to claim 6, It is characterized in that The step of determining the geological block where the ray starting point is located comprises: According to the triangulated network where the starting point of the ray is located, the geological blocks B1 and B2 adjacent to the triangulated network are obtained; Fine-tune the position coordinates of the ray starting point according to the preset adjustment rules; Taking the finely adjusted ray starting point position as the endpoint, construct the vertical segment b1 beyond the geological block B1 and the vertical segment b2 beyond the geological block B2; It is determined whether the number of intersections between the vertical line segment b1 and the geological block B1 and the number of intersections between the vertical line segment b2 and the geological block B2 are odd numbers, and the geological block with an odd number of intersections is determined to be the geological block where the ray starting point is located.

8. The method according to claim 1, It is characterized in that According to the symmetrical box intersecting with the initial ray and the corresponding triangle in the triangulated network included in the intersecting symmetrical box, the intersection point of the initial ray and the triangle in the triangulated network of the current geological block is determined to obtain an intersection point set, including: Traverse the flush symmetric box that intersects the initial ray; Determine whether the initial ray intersects with the triangle in the triangulated network contained in the intersecting symmetrical box; If yes, determine the intersection of the initial ray and the triangle in the triangulation network of the current geological block, and add the determined intersection to the intersection set; After traversing all the symmetrical boxes that intersect with the initial ray, we get the set of intersection points.

9. The method according to claim 8, It is characterized in that The step of determining the intersection of the initial ray and the triangle in the triangulation network of the current geological block comprises: Calculate the edge vectors v1 and v2 of the triangle according to the three vertices p1, p2, and p3 of the triangle; Determine the cross product xv of the edge vector based on the edge vector; Determine the ray segment vector vr based on the starting point r1 and the ending point r2 of the initial ray; Determine vector v3 based on the vertex p1 of the triangle and the starting point r1 of the initial ray; Determine the dot product d1 of the two according to the cross product xv and the vector vr; Determine the dot product d2 of the vector v3 and the cross product xv; Determine the point in the triangle according to the following formula: zp = r1 + (-d2 / d1) * vr, where zp represents a point in the triangle; Determine the area A of the triangle enclosed by vertices p1, p2, and p3, the area a1 of the triangle enclosed by vertices p1, p2, and zp, the area a2 of the triangle enclosed by vertices p2, p3, and zp, and the area a3 of the triangle p1, p3, and zp. If the value of (a1+a2+a3) / A is not greater than the preset area threshold, then the point zp in the triangle is the intersection of the initial ray and the triangle in the triangulation of the current geological block.

10. The method according to claim 1, It is characterized in that Also includes: Based on the ray end points determined by multiple initial rays passing through the ray starting point with a preset deflection angle, an imaging range of the ground surface included in the three-dimensional geological model of the rays emitted with the ray starting point as the starting point and with the preset deflection angle as the direction is obtained.

11. The method according to claim 10, It is characterized in that A plurality of initial rays passing through the ray starting point with a preset deflection angle include: Based on the selected ray starting point and the preset deflection angle relative to the formation normal at the ray starting point, a plurality of initial rays are obtained by taking the preset deflection angle of the formation normal as the ray direction and uniformly distributing the rays around the circumference at preset interval angles.

12. A ray tracing area analysis device, It is characterized in that include: A 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 levels and multiple geological blocks, each level includes multiple triangulated networks, each geological block is surrounded by triangulated networks of multiple levels to form a closed body, each triangulated network includes multiple triangles, and each geological block has multiple symmetrical boxes pre-constructed, each symmetrical box includes triangles in the triangulated network corresponding to the symmetrical box; A ray determination module, configured to obtain an initial ray based on a selected ray starting point and a preset deflection angle relative to a stratum normal at the ray starting point; Determining a next ray based on the formation velocities of two adjacent geological blocks and the ray starting points of the adjacent geological blocks; An intersection determination module is used to determine an symmetrical box in the geological block where the starting point of the ray is located that intersects with the initial ray; according to the symmetrical box that intersects with the initial ray and the corresponding triangle in the triangulation network included in the intersecting symmetrical box, determine the intersection of the initial ray and the triangle in the triangulation network of the current geological block to obtain an intersection set; Selecting an intersection point that meets a preset condition from the intersection point set as the path end point of the initial ray in the geological block where the ray starting point is located, and the ray starting point of the geological block adjacent to the geological block where the ray starting point is located; If the end point of the path 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 determination module is notified to continue the step of determining the symmetrical box that intersects with the initial ray in the current geological block where the starting point of the ray is located until the intersection with the geological block including the surface in the three-dimensional geological model is obtained, and the end point of the ray is obtained.

13. The device according to claim 12, It is 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 of the ray emitted from the ray starting point with a preset deflection angle as the direction based on the ray end point determined by multiple initial rays passing through the ray starting point.

14. A computer storage medium, It is characterized in that The computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, a ray tracing area analysis method as described in any one of claims 1-11 is implemented.

15. A terminal device, It is characterized in that include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, a ray tracing area analysis method as described in any one of claims 1 to 11 is implemented.

Images