Layered spacing evaluation method and device, electronic equipment and storage medium

Abstract

One or more embodiments of the invention provide a stratigraphic spacing evaluation method and device, electronic equipment and a storage medium. The method comprises the steps that a stratum triangulation network model of a target area is acquired, a reference stratum and a target stratum in the triangulation network model are determined, and the target stratum is represented by a plurality of triangular faces; constructing an octree representation of the triangular surface, wherein the octree comprises a plurality of leaf nodes; constructing a ray in a preset direction by taking a discrete point in the reference datum stratum as a starting point; leaf nodes intersecting with the rays in the octree are determined; performing intersection test on the triangular patches in the leaf nodes, and determining target triangular patches intersected with the rays in the leaf nodes; in response to determining that the intersection point of the ray and the target triangular patch is located in the target triangular patch, obtaining the stratum spacing between the reference stratum and the target stratum according to the coordinates of the starting point and the intersection point; and generating a stratigraphic spacing visualization model according to the stratigraphic spacing, and evaluating the stratigraphic spacing according to the stratigraphic spacing visualization model.

Description

TECHNICAL FIELD

[0001] One or more embodiments of the present disclosure relate to the technical field of geotechnical engineering, and in particular to a stratum interval evaluation method and device, electronic equipment and storage medium. BACKGROUND

[0002] It should be noted that the above introduction to the technical background is only for the convenience of clearly and completely describing the technical solutions of the present disclosure, and for the convenience of understanding by those skilled in the art. The above technical solutions cannot be considered as known to those skilled in the art merely because they are described in the background section of the present disclosure.

[0003] Coal resources occupy a dominant position in China's energy structure and are one of the most widely used and important energies in the world, providing an important driving force for global economic growth for a long time. Many safety problems and hidden dangers are involved in the process of coal mining, and the safe mining of coal is still an important foundation for energy security and economic stability.

[0004] In view of the research on coal mine safety accidents, China's traditional coal mine safety management is changing from "event center theory" to "risk center theory", from "feedback response-based" to "feedforward intervention-based", from "experience-driven" to "data-driven" coal mine safety risk management mode.

[0005] Among them, the accurate calculation of stratum interval is crucial for accurately evaluating roof stability, predicting water hazard risk, and designing a reasonable mine exploitation scheme. The traditional stratum interval evaluation method mainly relies on the measurement of two-dimensional cross-section, which cannot fully reflect the real distribution of strata in three-dimensional space. With the development of three-dimensional geological modeling technology, the stratum interval calculation method based on three-dimensional geometric model has gradually become a research hotspot.

[0006] However, the three-dimensional stratum interval evaluation method in the related art has the problems of limited calculation accuracy and low efficiency, and limited accuracy. It is of great significance to develop a high-precision and high-efficiency stratum interval evaluation method to realize more accurate address evaluation for coal mine safety production and rational development of resources. SUMMARY

[0007] Therefore, the purpose of one or more embodiments of the present disclosure is to propose a stratum interval evaluation method, device, electronic equipment and storage medium to solve the problems in the background art.

[0008] In order to achieve the above purpose, one or more embodiments of the present disclosure provide a stratum interval evaluation method, which comprises:

[0009] obtaining a stratigraphic triangulation model of a target area, and determining a reference base stratum and a target stratum in the triangulation model, the target stratum being represented by a plurality of triangular faces;

[0010] constructing an octree representation of the triangular faces, the octree comprising a plurality of leaf nodes;

[0011] constructing a ray along a preset direction starting from a discrete point in the reference base stratum;

[0012] determining a leaf node in the octree intersecting the ray;

[0013] performing an intersection test on the triangular patches within the leaf node to determine a target triangular patch within the leaf node intersecting the ray;

[0014] in response to determining that an intersection point of the ray and the target triangular patch is located within the target triangular patch, obtaining a stratigraphic interval between the reference base stratum and the target stratum according to coordinates of the starting point and the intersection point;

[0015] generating a stratigraphic interval visualization model according to the stratigraphic interval, and performing stratigraphic interval evaluation according to the stratigraphic interval visualization model.

[0016] Optionally, the constructing the octree representation of the triangular faces comprises:

[0017] constructing a bounding box according to spatial positions of all the triangular faces of the target stratum, and taking the bounding box as a root node;

[0018] dividing the root node until a number of triangular patches within a leaf node obtained by the division is less than or equal to a preset threshold;

[0019] dividing the triangular faces into at least one leaf node according to spatial positions of the triangular faces, each of the triangular faces belonging to one or more leaf nodes;

[0020] wherein the dividing the root node is iteratively performed until a number of triangular patches within a leaf node obtained by the division is less than or equal to a preset threshold:

[0021] determining a center point coordinate of each intermediate node in a current round;

[0022] dividing the intermediate node according to the center point coordinate to obtain child nodes of the intermediate node, each of the intermediate nodes being divided into eight child nodes;

[0023] taking the child nodes of the intermediate node as intermediate nodes in a next round.

[0024] Optionally, the determining the leaf node intersected by the ray in the octree comprises:

[0025] determining whether the ray intersects the root node according to the start point, the direction of the ray, and the spatial position of the root node;

[0026] in response to determining that the ray intersects the root node, recursively detecting the leaf node intersected by the ray;

[0027] wherein the following steps are iteratively performed to recursively detect the leaf node intersected by the ray:

[0028] determining the child node of the intermediate node intersected by the ray according to the start point, the direction of the ray, and the position coordinates of the intermediate node in the current round;

[0029] taking the child node as the intermediate node in the next round.

[0030] Optionally, the following steps are performed to determine the intersection point of the ray and the triangular patch:

[0031] constructing a plane equation of the triangular patch and a ray equation of the ray, the plane equation being constructed according to any vertex of the triangular patch;

[0032] obtaining an intersection point parameter according to the plane equation and the ray equation;

[0033] judging whether the intersection point is valid according to the intersection point parameter;

[0034] in response to determining that the intersection point is valid, calculating the intersection point of the ray and the triangular patch through intersection point coordinates.

[0035] Optionally, the following steps are performed to determine whether the intersection point is located within the triangular patch:

[0036] obtaining a first basis vector, a second basis vector, and a third basis vector according to vertex coordinates of the triangular patch and intersection point coordinates of the intersection point, the first basis vector being obtained according to vertex coordinates of a first vertex and a second vertex of the triangular patch, the second basis vector being obtained according to vertex coordinates of the first vertex and a third vertex of the triangular patch, and the third basis vector being obtained according to vertex coordinates of the first vertex of the triangular patch and the intersection point coordinates;

[0037] obtaining a parameter of a barycentric coordinate according to the first basis vector, the second basis vector, and the third basis vector;

[0038] in response to the barycentric coordinate satisfying validity detection, determining that the intersection point is located within the triangular patch.

[0039] Optionally, the stratum triangulation model is obtained by performing the following steps:

[0040] Obtaining a stratum point cloud data set of the target area;

[0041] Constructing a Delaunay triangulation grid according to the stratum point cloud data set;

[0042] Processing the Delaunay triangulation grid to obtain an optimized triangulation grid model;

[0043] Calculating a normal vector of each triangular face in the optimized triangulation grid model, and constructing the stratum triangulation model.

[0044] Optionally, the method further comprises:

[0045] Determining a stratum interval interval to which the stratum interval belongs and a corresponding mapping color according to the stratum interval and a stratum interval interval mapping table;

[0046] Rendering the stratum triangulation model according to the mapping color corresponding to the stratum interval and the start point coordinates of the start point, to obtain a stratum interval visualization model.

[0047] Based on the same inventive concept, one or more embodiments of the present disclosure also provide a stratum interval evaluation device, comprising:

[0048] An obtaining module configured to obtain a stratum triangulation model of a target area, and determine a reference datum stratum and a target stratum in the triangulation model, the target stratum being represented by a plurality of triangular faces;

[0049] A first constructing module configured to construct an octree representation of the triangular faces, the octree comprising a plurality of leaf nodes;

[0050] A second constructing module configured to construct a ray along a preset direction with a discrete point in the reference datum stratum as a start point;

[0051] A first positioning module configured to determine a leaf node in the octree that intersects with the ray;

[0052] A second positioning module configured to perform an intersection test on the triangular faces within the leaf node, and determine a target triangular face within the leaf node that intersects with the ray;

[0053] A calculating module configured to, in response to determining that an intersection point of the ray and the target triangular face is located within the target triangular face, obtain a stratum interval between the reference datum stratum and the target stratum according to coordinates of the start point and the intersection point.

[0054] An evaluation module is configured to generate a stratigraphic interval visualization model according to the stratigraphic intervals, and to perform stratigraphic interval evaluation according to the stratigraphic interval visualization model.

[0055] Based on the same inventive concept, one or more embodiments of the present disclosure further provide an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the stratigraphic interval evaluation method according to any one of the above when executing the program.

[0056] Based on the same inventive concept, one or more embodiments of the present disclosure further provide a non-transitory computer-readable storage medium, which stores computer instructions for causing the computer to execute the stratigraphic interval evaluation method according to any one of the above.

[0057] As can be seen from the above, the stratigraphic interval evaluation method provided by one or more embodiments of the present disclosure, after obtaining the stratigraphic triangulation network model of the target area, determining the reference datum stratum and the target stratum, constructing the octree representation of the triangular surface, emitting a preset direction ray from the reference datum stratum discrete point as the starting point, using the octree to quickly locate the leaf node intersected by the ray, performing intersection test on the triangular patches in the leaf node, when the intersection point of the ray and the target triangular patch is located in the patch, calculating the stratigraphic interval according to the coordinates of the starting point and the intersection point, and then generating a visualization model to assist stratigraphic interval evaluation. Its advantage lies in that the ray tracing complexity is optimized from O(n) to O(log n+k) (n represents the total number of triangular patches) with the help of octree spatial index, which greatly improves the stratigraphic interval calculation efficiency of large-scale complex geological model, combines the barycentric coordinate method to ensure the intersection point detection accuracy, and three-dimensionally visualizes the obtained stratigraphic interval, providing an integrated solution from data modeling to decision support in the fields of petroleum exploration and engineering geology.

[0058] The stratigraphic interval evaluation device, the electronic device, and the computer-readable storage medium provided by the present disclosure can all implement the steps of the stratigraphic interval evaluation method described above, and therefore also have the beneficial effects of the geological evaluation method described above. BRIEF DESCRIPTION OF DRAWINGS

[0059] In order to more clearly illustrate the technical solutions in the one or more embodiments of the present disclosure or the prior art, brief introductions will be given below to the drawings needed to be used in the embodiments or prior art descriptions. Obviously, the drawings in the following description are only one or more embodiments of the present disclosure, and other drawings can also be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0060] Figure 1 FIG. 1 is a flowchart of the stratigraphic interval evaluation method according to one or more embodiments of the present disclosure.

[0061] Figure 2 A structural schematic diagram of a stratum interlayer spacing evaluation device for one or more embodiments of the present disclosure;

[0062] Figure 3 An electronic device hardware structural schematic diagram for one or more embodiments of the present disclosure. DETAILED DESCRIPTION

[0063] In order to make the objectives, technical solutions and advantages of the present disclosure clearer, the present disclosure is further described in detail below with reference to specific embodiments and drawings.

[0064] It should be noted that, unless otherwise defined, technical terms or scientific terms used in one or more embodiments of the present disclosure should be understood as the general meaning understood by those skilled in the art to which the present disclosure belongs. The terms "first", "second", and similar terms used in one or more embodiments of the present disclosure do not represent any order, quantity, or importance, but are only used to distinguish different components. The terms "include", "contain", and similar terms mean that the elements or objects before the terms cover the elements or objects listed after the terms and their equivalents, and do not exclude other elements or objects. The terms "connect" or "connected" and similar terms are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. The terms "up", "down", "left", "right", and the like are only used to represent relative positional relationships, and when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0065] For the convenience of understanding, the technical terms such as reference datum stratum involved in the embodiments of the present application are briefly described below.

[0066] The reference datum stratum is a specific stratum unit selected as a spatial position reference standard in geological evaluation, which is represented as a discrete point set or a grid surface with continuous spatial form in a three-dimensional stratum triangulation model, and is often used as the starting point of ray emission to calculate the spatial distance from the target stratum. It is often selected as a landmark stratum with stable lithology, wide distribution, or clear age (such as a regional sedimentary interface), and the spatial relationship with the target stratum is quantified to provide a reference for geological structure analysis, resource evaluation, etc. For example, in oil exploration, a regional stable distribution of mudstone interlayer is often used as a reference datum stratum to calculate the thickness distribution of the overlying oil layer.

[0067] The stratigraphic triangulation model is a three-dimensional geometric model that discretely expresses the spatial form of the target area stratum through a series of connected triangular facets. The model can accurately depict the complex forms such as the relief and folds of the stratum by collecting the discrete point coordinates of the stratum interface and constructing the spatial topological relationship of the stratum surface in a triangular topological structure. The model not only retains the geometric characteristics (such as the strike and dip angle) of the stratum, but also reflects the spatial distribution law of the stratum through the spatial distribution of the triangular facets. It is a basic data structure for interlayer spacing calculation and structural analysis in geological assessment, and is commonly used in the digital expression and analysis of underground stratum form in the fields of three-dimensional geological modeling and resource exploration.

[0068] The triangular facet is the basic geometric unit of the stratigraphic triangulation model, which is formed by connecting three ordered vertices in three-dimensional space. Each vertex contains accurate coordinate information (x, y, z), and the topological relationship of the three edges forms a planar polygon with a normal vector (which can be calculated by vertex cross multiplication). In geological modeling, the triangular facet is used to discretely express the local form of the stratum interface. Multiple triangular facets are spliced together to form a continuous stratum surface, which can depict the relief, fractures and other detailed features of the stratum, and provide a basic geometric carrier for ray tracing, interlayer spacing calculation and other processes. It is a core element for realizing the digitization of stratum spatial form and geological assessment analysis.

[0069] As described in the background section, the three-dimensional stratigraphic interlayer spacing evaluation method in the related art has the problems of limited accuracy and low efficiency.

[0070] The three-dimensional stratigraphic interlayer spacing evaluation method in the related art is mostly based on grid division or interpolation estimation, which has limited calculation accuracy and low efficiency.

[0071] Especially in complex geological structure areas, due to the large changes in stratum strike and dip angle, the traditional method is difficult to accurately express the spatial relationship between the strata.

[0072] Therefore, the present disclosure proposes a stratigraphic interlayer spacing evaluation method, which constructs an octree index of the stratigraphic triangulation model, emits a ray from a reference base stratum discrete point and performs intersection detection with a target stratum triangular facet, calculates the interlayer spacing and generates a visualization model, and realizes efficient stratigraphic interlayer spacing evaluation.

[0073] Reference Figure 1 The stratigraphic interlayer spacing evaluation method of one or more embodiments of the present disclosure includes the following steps:

[0074] Step S101: Obtain the stratigraphic triangulation model of the target area, and determine the reference base stratum and the target stratum in the triangulation model. The target stratum is represented by a plurality of triangular facets.

[0075] Step S102: constructing an octree representation of the triangular facets, the octree including a plurality of leaf nodes;

[0076] Step S103: constructing a ray along a preset direction from a discrete point in the reference datum formation;

[0077] Step S104: determining a leaf node in the octree that intersects the ray;

[0078] Step S105: performing an intersection test on the triangular facets within the leaf node to determine a target triangular facet within the leaf node that intersects the ray;

[0079] Step S106: in response to determining that an intersection point of the ray and the target triangular facet is located within the target triangular facet, obtaining a formation interval between the reference datum formation and the target formation according to coordinates of the starting point and the intersection point;

[0080] Step S107: generating a formation interval visualization model according to the formation interval, and performing formation interval evaluation according to the formation interval visualization model.

[0081] In some implementations of the present disclosure, the formation triangular mesh model can be constructed based on formation point cloud data to obtain a Delaunay triangular mesh.

[0082] Specifically, in some implementations of the present disclosure, the step of constructing the formation triangular mesh model can be specifically:

[0083] First step, obtaining a formation point cloud data set, the formation including at least a reference datum formation and a target formation.

[0084] Second step, pre-processing the formation point cloud data set to obtain an optimized point cloud data set. In some embodiments of the present disclosure, the pre-processing can include denoising processing and down-sampling processing.

[0085] Third step, constructing a Delaunay triangular mesh according to the optimized point cloud data set.

[0086] The Delaunay triangular mesh is a special kind of triangular subdivision. In three-dimensional geological modeling, it constructs a triangular mesh based on a discrete point set to ensure that no other discrete points are contained in the circumcircle of any triangle (empty circle property), thereby forming a high-quality mesh that maximizes the minimum angle and avoids the appearance of narrow triangles. This mesh not only accurately fits the spatial form of the formation interface, but also has good topological structure and geometric properties, and can be effectively applied to formation surface reconstruction, spatial interpolation, and layer interval calculation scenarios, providing a regularized and high-quality basic mesh framework for three-dimensional geological model construction and analysis, and is an important technical means for realizing digital expression of formations.

[0087] The fourth step is to optimize the Delaunay triangular mesh to obtain an optimized triangular mesh model. In some embodiments of this disclosure, the optimization may include removing non-popular edges and smoothing processes.

[0088] The fifth step is to calculate the normal vector of each optimized triangular mesh model and construct a complete stratigraphic triangular mesh model.

[0089] In developing this disclosure, the applicant discovered that stratigraphic triangulation models typically contain a large number of triangular facets. Even if someone skilled in the art were to conceive of calculating stratigraphic spacing based on a stratigraphic triangulation model, for large geological models, the efficiency of a brute-force search method that traverses all triangular facets of the target stratum to determine the target triangular facet is extremely low.

[0090] Therefore, the technical solution disclosed herein introduces octree space partitioning technology to improve computational efficiency.

[0091] An octree is a hierarchical data structure in three-dimensional space. Its core principle is to recursively decompose the three-dimensional space into eight equal-sized cubic subspaces, forming a tree-like topology. Each node represents a spatial region, with the root node covering the entire target space and leaf nodes storing specific geometric objects (such as triangular faces). The advantage of this structure is that it can organize complex three-dimensional data hierarchically according to spatial location, significantly reducing the scope of data retrieval and computation by quickly locating the target region.

[0092] In this disclosed technical solution, the bounding box of the target stratigraphic triangular facet is used as the root node. The space is recursively divided according to a preset threshold, and the triangular facets are assigned to all intersecting child nodes, achieving spatialized data storage. During ray tracing, the intersection of the ray with the octree node is determined first. Precise intersection tests are performed only on the triangular facets within the intersecting nodes, avoiding a global traversal of all facets. This approach optimizes the interlayer spacing calculation efficiency from traditional linear time complexity to logarithmic complexity, making it particularly suitable for the rapid analysis of large-scale stratigraphic models and providing an efficient data processing foundation for geological assessment.

[0093] In this implementation, the leaf nodes are used to store triangular faces. The number of triangular faces stored in each leaf node does not exceed a preset threshold, and a triangular face may belong to one or more leaf nodes.

[0094] In some embodiments of this disclosure, the method for constructing the octree representation of the triangular face described above may include:

[0095] Step S201: Construct a bounding box based on the spatial positions of all the triangular faces of the target stratum, and use the bounding box as the root node;

[0096] Step S202: Divide the root node until the number of triangular faces in the resulting leaf nodes is less than or equal to a preset threshold.

[0097] Step S203: Based on the spatial position of the triangular facets, divide the triangular facets into at least one of the leaf nodes, with each triangular facet belonging to one or more of the leaf nodes;

[0098] The following steps are performed iteratively to divide the root node until the number of triangular faces in the resulting leaf nodes is less than or equal to a preset threshold:

[0099] Step S2021: Determine the coordinates of the center point of each intermediate node in the current round;

[0100] Step S2022: Divide the intermediate nodes based on the coordinates of the center point to obtain the child nodes of the intermediate nodes. Each intermediate node is divided into eight child nodes.

[0101] Step S2023: Use the child nodes of the above intermediate nodes as the intermediate nodes of the next round.

[0102] In other words, in each round, all existing nodes are checked to determine the number of triangular faces in each node. When the number of triangular faces in any node exceeds a preset value, the space of all nodes in the current round is partitioned, with each node divided into eight parts.

[0103] For example, in one embodiment of this disclosure, the root node is first spatially divided into 8 parts to obtain 8 intermediate nodes. It is determined that at least one intermediate node has a greater than a preset value of triangular facets, and the division continues, with each intermediate node further divided into 8 parts, resulting in 64 intermediate nodes. It is then determined that the number of triangular facets stored in these 64 intermediate nodes is less than the preset value, and these 64 intermediate nodes are designated as leaf nodes.

[0104] Thus, in the octree representation of this embodiment, there is one root node, eight intermediate nodes below the root node, and 64 leaf nodes below the level of the intermediate nodes.

[0105] In another embodiment of this disclosure, the root node is first spatially divided into 8 parts, resulting in 8 intermediate nodes. If at least one intermediate node has a greater than a preset value of triangular facets, the division continues, with each intermediate node further divided into 8 parts, resulting in 64 intermediate nodes. If at least one of these 64 intermediate nodes still has a greater than the preset value of triangular facets, the division continues, with each intermediate node further divided into 8 parts, resulting in 4096 intermediate nodes. If it is determined that the number of triangular facets stored in all 4096 intermediate nodes is less than the preset value, these 4096 intermediate nodes are designated as leaf nodes.

[0106] Thus, in the octree representation of this embodiment, there is one root node, eight intermediate nodes below the root node, 64 intermediate nodes below the level of these eight intermediate nodes, and 4096 leaf nodes below the level of these 64 intermediate nodes.

[0107] By employing octree spatial partitioning, the computational complexity of intersection detection can be reduced from O(n) to an average of O(log n). In actual testing, for a stratigraphic model containing 100,000 triangular patches, the computation time was reduced by approximately 95%, from tens of seconds in the traditional method to less than one second, significantly improving the efficiency of calculating stratigraphic interlayer spacing distribution.

[0108] The rays constructed in step S103 are used to perform intersection tests with the triangular facets of the target stratum's octree index. By determining the coordinates of the effective intersection points, the interlayer spacing between the reference stratum and the target stratum is calculated, providing key data for generating a visualization model.

[0109] In some implementations of this disclosure, the starting point of the ray (a discrete point in the reference stratum) can be represented as P. s (x s ,y s ,z s ).

[0110] In some implementations of this disclosure, to obtain the definition formula for the ray direction, we can first define the ray as pointing to a simulated point P. t (x s ,y s ,z s +Δz), where Δz represents the height increment between the simulated point and the discrete point in the reference stratum.

[0111] Thus, the direction of the ray can be defined as:

[0112] D = P t -P s =(x s ,y s ,z s +Δz)-(x s ,y s ,z s )=(0,0,Δz);

[0113] Among them, P s (x s ,y s ,z s () indicates the starting point.

[0114] Furthermore, the ray emitted from a discrete point (starting point) of the reference stratum can be defined as:

[0115] R(t)=P s +t·D,t≥0;

[0116] In the technical solution disclosed herein, the triangular facets that intersect with the aforementioned rays in the target stratum are defined as target triangular facets, and the stratum spacing is determined based on the distance between the starting point and the target triangular facets.

[0117] To achieve rapid positioning of target triangular facets, the technical solution disclosed herein introduces octree space partitioning technology.

[0118] In some implementations of this disclosure, after introducing octree space partitioning technology, the process of determining the target triangular facet is the above steps S104 to S105.

[0119] Regarding step S104, determining the leaf nodes in the octree that intersect with the ray, some embodiments of this disclosure can be refined as follows:

[0120] Step S401: Based on the starting point and direction of the ray and the spatial position of the root node, determine whether the ray intersects with the root node.

[0121] Step S402: In response to determining that the above ray intersects with the root node, recursively detect the leaf nodes that intersect with the above ray;

[0122] The process iteratively executes the following steps to recursively detect the leaf nodes that intersect with the aforementioned rays:

[0123] Step S4021: Based on the starting point and direction of the ray and the position coordinates of the intermediate node in the current round, determine the child node of the intermediate node that intersects with the ray.

[0124] Step S4022: Use the above child nodes as the intermediate nodes for the next round.

[0125] The above steps can be understood as follows: when determining the leaf nodes that intersect with the ray, it is first necessary to determine whether the ray intersects with the root node. If the ray does not intersect with the root node, it means that the ray does not intersect with any leaf nodes within the root node, and no further calculation is needed. Conversely, if the ray intersects with the root node, it means that the ray must intersect with at least one leaf node within the root node, and the calculation continues.

[0126] In some implementations of this disclosure, leaf nodes intersecting with rays are determined using recursive algorithms to improve computational efficiency.

[0127] That is, after determining that the ray intersects with the root node, the intermediate nodes that intersect with the ray in the next level are determined iteratively until the leaf nodes that intersect with the ray are determined in the last round.

[0128] In some implementations of this disclosure, the intermediate nodes (or leaf nodes) of each level are sorted according to their spatial positions to obtain a sequence of intermediate nodes (or leaf nodes).

[0129] In one embodiment of this disclosure, in the first iteration, it is calculated whether the ray intersects with the 8 intermediate nodes, and it is determined that the ray intersects with the 5th intermediate node of this level. In the second iteration, it is calculated whether the ray intersects with the 8 leaf nodes under the aforementioned 5th intermediate node, and it is determined that the ray intersects with the 7th and 8th leaf nodes among the 8 leaf nodes.

[0130] In another embodiment of this disclosure, in the first iteration of the first round, it is calculated whether the ray intersects with the 8 intermediate nodes, and it is determined that the ray intersects with the 3rd intermediate node of this level. In the second iteration, it is calculated whether the ray intersects with the 8 intermediate nodes under the aforementioned 3rd intermediate node, and it is determined that the ray intersects with the 4th and 5th intermediate nodes of this level. In the third iteration, it is calculated whether the ray intersects with the 16 leaf nodes under the aforementioned 4th and 5th intermediate nodes, and it is determined that the ray intersects with the 8th leaf node under the 4th intermediate node and the 1st leaf node under the 5th intermediate node.

[0131] The intersection detection in step S105 is to determine whether the intersection point between the ray and the triangular facet in the leaf node is valid, including whether the intersection point is inside the triangular facet and not on the edge of the triangular facet.

[0132] In an implementation of this disclosure, the above steps may include:

[0133] Based on the normal vectors of the triangular faces within the intersecting leaf nodes, the coordinate representation of the intersection point is obtained;

[0134] Determine whether the above intersection points are within the triangular facet using the centroid coordinate method.

[0135] Specifically, in the implementation of this disclosure, the equation of the plane containing each triangular facet is defined as: a(x-x1)+b(y-y1)+c(z-z1)=0, using the normal vector N=(a,b,c) of each triangular facet.

[0136] The above plane equation can be expanded as: ax + by + cz + d = 0, where d = -(ax1 + by1 + cz1).

[0137] Substitute the ray equation from the above into the plane equation: The intersection point can be obtained: P = (x s +t·d x ,y s +t·d y ,z s +t·d z ), where dx d y and d z It is a component of the ray direction vector, d z =Δz.

[0138] According to the centroid coordinate method, the vertices of the triangle face can be defined as V1, V2 and V3, and the intersection point can be expressed as: P = V1 + u(V2-V1) + v(V3-V1).

[0139] Thus, u and v can be represented as:

[0140]

[0141] If the following conditions are met: 0≤u≤1, 0≤v≤1, and u+v≤1, then the intersection point is inside the triangular facet; otherwise, the intersection point is invalid.

[0142] Based on the starting point coordinates and the distance between the effective intersection points, the stratigraphic distance between the reference stratum and the target stratum can be determined. The formula for calculating this distance is as follows:

[0143] Some implementations of this disclosure, including visualization of stratigraphic interlayer data, may include:

[0144] Step S701: Based on the above stratigraphic spacing and stratigraphic spacing interval mapping table, determine the stratigraphic spacing interval to which the above stratigraphic spacing belongs and the corresponding mapping color;

[0145] Step S702: Render the above-mentioned stratigraphic triangulation model according to the mapping color corresponding to the above-mentioned stratigraphic spacing and the starting coordinate of the starting point, so as to obtain the above-mentioned stratigraphic spacing visualization model.

[0146] In some embodiments of this disclosure, the above-mentioned stratigraphic spacing range can be: 0-5 meters, 5-15 meters, 15-25 meters, and greater than 25 meters.

[0147] Each zone is configured with a different display method, differing in color or style. Users can intuitively understand the stratigraphic spacing in different areas through the display of different stratigraphic spacings.

[0148] Some implementations disclosed herein, after determining the stratigraphic spacing interval to which the stratigraphic spacing belongs, can also determine the corresponding stratigraphic spacing assessment report, such as the hidden key disaster-causing factors, making it convenient for users to view the corresponding stratigraphic spacing assessment results.

[0149] Some implementations of this disclosure can also perform the following analyses and evaluations based on the aforementioned stratigraphic spacing data. Some implementations of this disclosure can also simultaneously display the analysis and evaluation results in a stratigraphic spacing visualization model.

[0150] Specifically, some embodiments of this disclosure may include the analysis and evaluation of stratigraphic interlayer data as follows:

[0151] Based on stratigraphic interlayer data, statistical characteristic values ​​are calculated, which may include one or more of the following: minimum value, maximum value, average value, and standard deviation.

[0152] Based on the above statistical characteristic values, abnormal regions in the target area are identified. These abnormal regions may include abnormally thin regions and abnormally thick regions.

[0153] Based on the acquired geological structural data, the aforementioned anomalous areas can be analyzed to assess potential geological risks.

[0154] Based on the aforementioned geological risk data, a stratigraphic spacing analysis report is generated to provide a basis for subsequent stratigraphic spacing assessment and geological assessment.

[0155] It is understandable that this method can be executed by any device, equipment, platform, or cluster of devices with computing and processing capabilities.

[0156] It should be noted that the methods of one or more embodiments of this disclosure can be executed by a single device, such as a computer or server. The methods of this embodiment can also be applied in a distributed scenario, where multiple devices cooperate to complete the process. In such a distributed scenario, one of these devices may execute only one or more steps of the methods of one or more embodiments of this disclosure, and the multiple devices will interact with each other to complete the method described.

[0157] It should be noted that the above description pertains to specific embodiments of this disclosure. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps described in the claims may be performed in a different order than those shown in the embodiments and may still achieve the desired results. Furthermore, the processes depicted in the drawings do not necessarily require a specific or sequential order to achieve the desired results. In some embodiments, multitasking and parallel processing are possible or may be advantageous.

[0158] Based on the same inventive concept, corresponding to any of the methods in the above embodiments, this disclosure also provides a geological assessment device. For example... Figure 2 As shown, the stratigraphic interlayer assessment device includes:

[0159] The acquisition module 11 is configured to acquire a stratigraphic triangulation model of the target area and determine the reference stratigraphic layer and the target stratigraphic layer in the triangulation model, wherein the target stratigraphic layer is represented by multiple triangular faces;

[0160] The first construction module 12 is configured to construct an octree representation of the triangular face, the octree including multiple leaf nodes;

[0161] The second construction module 13 is configured to construct a ray along a preset direction, starting from a discrete point in the reference stratum;

[0162] The first positioning module 14 is configured to determine the leaf nodes in the octree that intersect with the ray;

[0163] The second positioning module 15 is configured to perform an intersection test on the triangular facets within the leaf node to determine the target triangular facets within the leaf node that intersect with the ray.

[0164] The calculation module 16 is configured to, in response to determining that the intersection of the ray and the target triangular facet is located within the target triangular facet, obtain the stratigraphic distance between the reference stratum and the target stratum based on the coordinates of the starting point and the intersection point;

[0165] The evaluation module 17 is configured to generate a formation spacing visualization model based on the formation spacing and to evaluate the formation spacing based on the formation spacing visualization model.

[0166] For ease of description, the above apparatus is described in terms of its functions, divided into various modules. Of course, when implementing one or more embodiments of this disclosure, the functions of each module can be implemented in one or more software and / or hardware.

[0167] The apparatus described above is used to implement the corresponding methods in the foregoing embodiments and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0168] Figure 3 This embodiment illustrates a more specific hardware structure of an electronic device, which may include a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, memory 1020, input / output interface 1030, and communication interface 1040 are interconnected internally via the bus 1050.

[0169] The processor 1010 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this disclosure.

[0170] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 1020 can store the operating system and other applications. When the technical solutions provided in the embodiments of this disclosure are implemented by software or firmware, the relevant program code is stored in the memory 1020 and is called and executed by the processor 1010.

[0171] The input / output interface 1030 is used to connect input / output modules to realize information input and output. Input / output modules can be configured as components within the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Input devices may include keyboards, mice, touchscreens, microphones, various sensors, etc., while output devices may include displays, speakers, vibrators, indicator lights, etc.

[0172] The communication interface 1040 is used to connect a communication module (not shown in the figure) to enable communication between this device and other devices. The communication module can communicate via wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0173] Bus 1050 includes a pathway for transmitting information between various components of the device, such as processor 1010, memory 1020, input / output interface 1030, and communication interface 1040.

[0174] It should be noted that although the above-described device only shows the processor 1010, memory 1020, input / output interface 1030, communication interface 1040, and bus 1050, in specific implementations, the device may also include other components necessary for normal operation. Furthermore, those skilled in the art will understand that the above-described device may only include the components necessary for implementing the embodiments of this disclosure, and not necessarily all the components shown in the figures.

[0175] The electronic devices described above are used to implement the corresponding methods in the foregoing embodiments and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0176] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. Information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by a computing device.

[0177] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of this disclosure (including the claims) is limited to these examples; within the framework of this disclosure, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of one or more embodiments of this disclosure as described above, which are not provided in detail for the sake of brevity.

[0178] Additionally, to simplify the description and discussion, and to avoid obscuring one or more embodiments of this disclosure, the provided drawings may or may not show well-known power / ground connections to integrated circuit (IC) chips and other components. Furthermore, the apparatus may be shown in block diagram form to avoid obscuring one or more embodiments of this disclosure, and this also takes into account the fact that the details of implementation of these block diagram apparatuses are highly dependent on the platform on which one or more embodiments of this disclosure will be implemented (i.e., these details should be fully understood by those skilled in the art). While specific details (e.g., circuitry) have been set forth to describe exemplary embodiments of this disclosure, it will be apparent to those skilled in the art that one or more embodiments of this disclosure may be implemented without these specific details or with variations thereof. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0179] Although this disclosure has been described in conjunction with specific embodiments thereof, many substitutions, modifications, and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0180] This disclosure includes one or more embodiments intended to cover all such substitutions, modifications, and variations falling within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of one or more embodiments of this disclosure should be included within the scope of protection of this disclosure.

Claims

1. A geological assessment method, characterized in that, include: Obtain a stratigraphic triangulation model of the target area, and determine the reference stratigraphic layer and the target stratigraphic layer in the triangulation model. The target stratigraphic layer is represented by multiple triangular faces. Construct an octree representation of the triangular face, wherein the octree includes multiple leaf nodes; Construct a ray along a preset direction, starting from a discrete point in the reference stratum; Determine the leaf nodes in the octree that intersect with the ray; Intersection tests are performed on the triangular facets within the leaf node to determine the target triangular facets within the leaf node that intersect with the ray; In response to determining that the intersection of the ray and the target triangular facet is located within the target triangular facet, the stratigraphic distance between the reference stratum and the target stratum is obtained based on the coordinates of the starting point and the intersection point; A formation spacing visualization model is generated based on the formation spacing, and the formation spacing is evaluated based on the formation spacing visualization model.

2. The method according to claim 1, characterized in that, The construction of the octree representation of the triangular face includes: Construct a bounding box based on the spatial positions of all triangular faces of the target stratum, and use the bounding box as the root node; The root node is divided until the number of triangular facets in the resulting leaf nodes is less than or equal to a preset threshold. Based on the spatial position of the triangular facet, the triangular facet is divided into at least one leaf node, and each triangular facet belongs to one or more leaf nodes; The root node is divided iteratively until the number of triangular faces within the resulting leaf nodes is less than or equal to a preset threshold: Determine the coordinates of the center point of each intermediate node in the current round; The intermediate nodes are divided based on the coordinates of the center point to obtain the child nodes of the intermediate nodes, and each intermediate node is divided into eight child nodes; The child nodes of the intermediate node are used as the intermediate nodes in the next round.

3. The method according to claim 2, wherein determining the leaf nodes in the octree that intersect with the ray comprises: Based on the starting point and direction of the ray, and the spatial position of the root node, determine whether the ray intersects with the root node; In response to determining that the ray intersects the root node, the leaf nodes that intersect the ray are recursively detected; The process involves iteratively performing the following steps to recursively detect the leaf nodes that intersect with the ray: Based on the starting point and direction of the ray, and the position coordinates of the intermediate node in the current round, determine the child nodes of the intermediate node that intersect with the ray; The child node will be used as the intermediate node in the next round.

4. The method according to claim 1, characterized in that, The following steps are performed to determine the intersection point of the ray and the triangular facet: Construct the plane equation of the triangular facet and the ray equation of the ray, wherein the plane equation is constructed based on any vertex of the triangular facet; Based on the plane equation and the ray equation, the intersection parameters are obtained; Based on the intersection parameters, determine whether the intersection is valid; In response to determining that the intersection point is valid, the intersection point between the ray and the triangular facet is calculated using the intersection point coordinates.

5. The method according to claim 4, characterized in that, The following steps are performed to determine whether the intersection point lies within the triangular facet: Based on the vertex coordinates of the triangular facet and the intersection coordinates of the intersection point, a first basis vector, a second basis vector, and a third basis vector are obtained. The first basis vector is obtained based on the vertex coordinates of the first and second vertices of the triangular facet. The second basis vector is obtained based on the vertex coordinates of the first and third vertices of the triangular facet. The third basis vector is obtained based on the vertex coordinates of the first vertex of the triangular facet and the intersection coordinates. Based on the first basis vector, the second basis vector, and the third basis vector, the parameters of the centroid coordinates are obtained; In response to the validity check being satisfied by the centroid coordinates, it is determined that the intersection point is located within the triangular facet.

6. The method according to claim 1, characterized in that, The stratigraphic triangulation model is obtained by performing the following steps: Obtain the stratigraphic point cloud dataset of the target area; Based on the aforementioned stratigraphic point cloud dataset, a Delaunay triangular mesh is constructed; The Delaunay triangular mesh is processed to obtain an optimized triangular mesh model; Calculate the normal vector of each triangle face in the optimized triangular mesh model, and construct the stratigraphic triangular mesh model.

7. The method according to claim 1, characterized in that, Also includes: Based on the stratigraphic spacing and stratigraphic spacing interval mapping table, determine the stratigraphic spacing interval to which the stratigraphic spacing belongs and the corresponding mapping color; The stratigraphic triangulation model is rendered based on the mapping color corresponding to the stratigraphic spacing and the starting coordinates of the starting point to obtain the stratigraphic spacing visualization model.

8. A geological assessment device, characterized in that, include: The acquisition module is configured to acquire a stratigraphic triangulation model of the target area and determine the reference stratigraphic layer and the target stratigraphic layer in the triangulation model, wherein the target stratigraphic layer is represented by multiple triangular faces; The first building module is configured to build an octree representation of the triangular face, the octree including multiple leaf nodes; The second construction module is configured to construct a ray along a preset direction, starting from a discrete point in the reference stratum; The first positioning module is configured to determine the leaf nodes in the octree that intersect with the ray; The second positioning module is configured to perform an intersection test on the triangular facets within the leaf node to determine the target triangular facets within the leaf node that intersect with the ray. The calculation module is configured to, in response to determining that the intersection of the ray and the target triangular facet is located within the target triangular facet, obtain the stratigraphic distance between the reference stratum and the target stratum based on the coordinates of the starting point and the intersection point; The evaluation module is configured to generate a formation spacing visualization model based on the formation spacing and to evaluate the formation spacing based on the formation spacing visualization model.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executed by the processor, characterized in that, When the processor executes the computer program, it implements the method as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions for causing the computer to perform the method of any one of claims 1 to 7.

Images