A three-dimensional sectioning method for generalized triangular prisms with vertical side edges

Abstract

The present invention discloses a three-dimensional sectioning method of a generalized triangular prism with vertical side edges, which relates to the technical field of three-dimensional geological information models. The method includes obtaining a VGTP unit that needs to be sectioned, determining an effective sectioning plane, and sectioning and segmenting the unit. Specifically, the VGTP unit that needs to be sectioned is determined by the intersection relationship between all sectioning planes and the VGTP unit, and then the effective sectioning plane is determined by the inclusion relationship of the sectioned body. Finally, a series of vertical sections are used to segment the sectioned VGTP unit to form a new VGTP unit. The sectioning of the entire three-dimensional geological information model geometric body is finally achieved by cyclically executing these steps. The present invention reduces the actual sectioning times of the VGTP unit by determining the effective sectioning plane, reduces the model data volume, and improves the model operation efficiency. At the same time, through a strict sectioning method of the VGTP unit, it is ensured that all units in the model will not mutate into non-VGTP units after sectioning.

Description

Technical Field

[0001] The invention relates to the technical field of three-dimensional geological information models, in particular to a three-dimensional sectioning method of a generalized triangular prism with vertical side edges. Background Art

[0002] The spatial geometry of three-dimensional geological information models is mainly expressed in the form of surfaces or volumes. Volume expression is considered to be a more suitable expression method for the spatial geometry of three-dimensional geological information models because of its characteristics of rapid volume acquisition and convenient spatial calculation.

[0003] The spatial geometry of the three-dimensional geological information model is generally constructed by interpolating the formation interface through drilling, and then the spatial geometry is established through the formation interface. Based on this feature, early scholars recommended the use of quasi-triangular prisms to establish the spatial geometry of the three-dimensional geological information model, that is, first use the triangular network to establish the formation interface, and then connect the triangular networks in the upper and lower adjacent formation interfaces to form a quasi-triangular prism (Analogical Tri-Prism, referred to as ATP), and the spatial geometry of the three-dimensional geological information model is composed of multiple quasi-triangular prisms. Later, in order to better simulate special geological structures such as pinch-out and lens, the generalized tri-prism (Generalize Tri-Prism, referred to as GTP) was derived from the quasi-triangular prism, and it has become a widely recognized way to express the spatial geometry of the three-dimensional geological information model.

[0004] In order to further solve the problem of excavation of three-dimensional geological information models, some scholars have conducted research on the problem of multiple sectioning of arbitrary spatial surfaces of generalized triangular prism models. For example, the invention patent application with publication number CN104318618A, "Three-dimensional Sectioning of Generalized Triangular Prism Spatial Data Model", re-divided the GTP through tetrahedral GTP units after it was sectioned to ensure that all units in the model still meet the characteristics of GTP units, thus realizing multiple sectioning of arbitrary spatial surfaces of the model. However, this method still has the following disadvantages:

[0005] (1) The generalized triangular prism (GTP) unit performs sectioning and segmentation on all section planes, and after each sectioning, the remaining part is further divided into multiple tetrahedral GTP units. As the number of sectioning increases, a large number of GTP units will be generated at the sectioning position, which increases the model storage capacity while reducing the calculation speed of model operation;

[0006] (2) This method uses the GTP unit, which does not require the side edges to be vertical. In comparison, the generalized triangular prism (VGTP) unit with vertical side edges is more efficient in terms of its intersection relationship with the cutting plane and its own volume calculation. Summary of the invention

[0007] In order to overcome the defects in the above-mentioned prior art, the present invention provides a three-dimensional sectioning method of a generalized triangular prism with vertical side edges, which utilizes the positional relationship between the spatial sectioning plane and the generalized triangular prism unit with vertical side edges to reduce the number of times the unit is sectioned, and at the same time, strictly uses the generalized triangular prism unit with vertical side edges to re-divide the units after multiple sectioning, so as to ensure that the model after excavation is still completely composed of generalized triangular prisms with vertical side edges.

[0008] To achieve the above object, the present invention adopts the following technical solutions, including:

[0009] A three-dimensional sectioning method for a generalized triangular prism with vertical side edges, wherein the generalized triangular prism unit with vertical side edges, i.e., a VGTP unit, is surrounded by two upper and lower triangular bases and two / three side faces; the projections of the upper and lower triangular bases along the vertical direction completely overlap, the sides constituting the triangular bases are called base sides, and the vertices of the base sides are called base vertices; the sides of the side faces that intersect and completely overlap in pairs are called side edges, and the sides of the side faces except the side edges completely overlap with the base sides of the upper and lower triangular bases respectively;

[0010] The three-dimensional cutting method comprises the following steps:

[0011] Step 1: Obtain the VGTP unit to be sectioned

[0012] If the VGTP unit has a section plane that intersects with it, it is determined that the VGTP unit needs to be sectioned, and all the section planes that intersect with it are recorded accordingly; otherwise, the VGTP unit does not need to be sectioned;

[0013] Step 2: Determine the effective sectioning plane

[0014] Obtain a VGTP unit that needs to be sectioned and all its corresponding section planes, and determine the judgment method according to the total number n of all the corresponding section planes, specifically: if n=1, the only section plane is a valid section plane; if n>1, let all the section planes independently section the VGTP unit, and obtain the body after all the section planes independently section the VGTP unit, where the body of the VGTP unit that is removed after independent sectioning is recorded as a removed body, and each section plane corresponds to a removed body; compare the removed bodies two by two, if one of the removed bodies completely contains the other removed body, then define the section plane corresponding to the completely contained removed body as an invalid section plane, repeat step 2, eliminate all invalid section planes corresponding to the VGTP units that need to be sectioned, and the remaining section planes are recorded as valid section planes;

[0015] Step 3: VGTP unit sectioning and segmentation

[0016] Step 3.1, obtaining a VGTP unit to be sectioned and all corresponding valid section planes;

[0017] Step 3.2, select a valid cutting plane, calculate the intersection surface of the valid cutting plane and the VGTP unit, record the direction perpendicular to the intersection surface and not pointing to the cut body as the normal of the intersection surface; remove the cut body and obtain the vertical direction of the normal of the intersection surface, record it as the projection direction; determine the projection method of the intersection surface according to the projection direction, specifically: if the projection direction points vertically upward, then project the intersection surface along the vertical direction to the upper bottom surface of the original VGTP unit; if the projection direction points vertically downward, then project the intersection surface along the vertical direction to the lower bottom surface of the original VGTP unit; if the projection direction cannot be obtained, that is, the normal of the intersection surface is completely horizontal, then project the intersection surface along the vertical direction to the upper or lower bottom surface of the original VGTP unit, and define the projection direction as vertically upward or vertically downward accordingly;

[0018] Step 3.3, according to the relationship between the projection surface of the intersecting surface and the remaining part of the upper or lower bottom surface of the VGTP unit after the dissection body is removed, determine the further operation between the two, specifically: if the two completely overlap, then determine this completely overlapping surface as the surface to be divided; if the two do not completely overlap, then the edge of the outer contour of the projection surface that does not overlap with the outer contour of the remaining part of the upper or lower bottom surface is recorded as a trace line, and the remaining part of the upper or lower bottom surface is determined as the surface to be divided; the endpoints of the trace line and the vertices of the surface to be divided are both recorded as segmentation points; the selection of the remaining part of the upper or lower bottom surface needs to be determined according to the projection direction, specifically: if the projection direction is vertically downward, it is the remaining part of the lower bottom surface; otherwise, it is the remaining part of the upper bottom surface; the projection surface of the intersecting surface degenerates into a projection line when the intersecting surface is completely vertical;

[0019] Step 3.4, split the surface to be split into multiple triangular sub-surfaces, the vertices of the triangular sub-surfaces all coincide with the splitting points; if there is a trace line in the surface to be split, then there are also common edges of two triangular sub-surfaces in the triangular sub-surfaces that coincide with the trace line;

[0020] Step 3.5, obtain the common edges between all the triangular sub-faces in the face to be divided, recorded as the dividing line; stretch the dividing line up and down in the vertical direction to form a dividing face; use the dividing face to divide and remove the VGTP unit after the dissected body to form multiple new VGTP units; the upper and lower bottom surfaces of the new VGTP unit need to be determined according to the projection direction in step 3.2, as follows: if the projection direction is vertically downward, then the surface of the new VGTP unit that is located in the same spatial surface as the face to be divided in step 3.3 is determined as the lower bottom surface of the new VGTP unit, and the surface of the new VGTP unit that completely overlaps with the vertical projection of the lower bottom surface of the new VGTP unit is determined as the upper bottom surface of the new VGTP unit; if the projection direction is vertically upward, then the surface of the new VGTP unit that is located in the same spatial surface as the face to be divided in step 3.3 is determined as the upper bottom surface of the new VGTP unit, and the surface of the new VGTP unit that completely overlaps with the vertical projection of the upper bottom surface of the new VGTP unit is determined as the lower bottom surface of the new VGTP unit;

[0021] Step 3.6, select a new effective cutting plane, determine the VGTP unit that needs to be cut by the new effective cutting plane in the new VGTP unit in step 3.5 according to the method in step 1, and cut and split according to the method described in steps 3.2-3.5; repeat this step until all effective cutting planes complete the cutting and splitting of the new VGTP unit generated by the VGTP unit in step 3.1 during the cyclic execution of steps 3.2-3.5;

[0022] Repeat step 3 until all VGTP units that need to be sectioned are sectioned and segmented.

[0023] Preferably, if the two side surfaces of the VGTP unit are two triangular side surfaces, the VGTP unit contains only one vertical side edge, and the VGTP unit is called a VGTP unit with one side edge; if the three side surfaces of the VGTP unit are two triangular side surfaces and one quadrilateral side surface, the VGTP unit contains two vertical side edges, and the VGTP unit is called a VGTP unit with two side edges; if the three side surfaces of the VGTP unit are all quadrilateral side surfaces, the VGTP unit contains three vertical side edges, and the VGTP unit is called a VGTP unit with three side edges.

[0024] Preferably, the specific process of step 1 is as follows:

[0025] Obtain a section plane equation αX+βY+χZ+η=0, wherein any spatial point (X′, Y′, Z′) satisfying the section plane equation is located on the section plane; select a VGTP unit and obtain the spatial coordinates of the six bottom surface vertices of the upper and lower bottom surfaces, wherein the spatial coordinates are composed of plane coordinates and vertical elevations;

[0026] Substitute the plane coordinates of the three vertices of the upper or lower bottom surface into the section plane equation to obtain three vertical elevations, and form three spatial points with the corresponding plane coordinates, which are recorded as the edge intersection points;

[0027] According to the positional relationship between the edge intersection point and the upper and lower bottom surface vertices of the VGTP unit, the edge intersection point is classified, specifically: if the vertical elevation of the edge intersection point is greater than the vertical elevations of the upper and lower bottom surface vertices at the same time, then the edge intersection point is recorded as the upper virtual edge intersection point; if the vertical elevation of the edge intersection point is less than the vertical elevations of the upper and lower bottom surface vertices at the same time, then the edge intersection point is recorded as the lower virtual edge intersection point; otherwise, the edge intersection point is recorded as a real edge intersection point; the edge intersection point and its corresponding upper or lower bottom surface vertex satisfy a positional relationship of coincidence along the vertical projection;

[0028] According to the classification of the edge intersections, it is determined whether the VGTP unit needs to be sectioned. Specifically, if the three edge intersections are all upper virtual edge intersections or all lower virtual edge intersections, it is determined that the VGTP unit does not need to be sectioned; otherwise, it is determined that the VGTP unit needs to be sectioned, and the section surface is recorded accordingly;

[0029] Obtain the equations of the remaining section planes, repeat step 1, and obtain all the VGTP units that need to be sectioned and all the section planes corresponding to each unit.

[0030] Preferably, in step 2, the dissected body that is completely contained is recorded as dissected body B, and the other dissected body compared with it is recorded as dissected body A; the complete containment means that the outer contour of the dissected body B does not exceed the outer contour of the dissected body A.

[0031] Preferably, in step 2, the normal direction of the section plane is any one of the two directions perpendicular to the section plane, and the sectioned body is not in the direction indicated by the normal direction of the section plane.

[0032] A readable storage medium stores a computer program, wherein the computer program includes program instructions, and when the program instructions are executed by a processor, the processor executes the above-mentioned three-dimensional sectioning method of a generalized triangular prism with vertical side edges.

[0033] An electronic device includes a processor, an input device, an output device and a memory, wherein the processor, the input device, the output device and the memory are connected in sequence, the memory is used to store a computer program, the computer program includes program instructions, and the processor is configured to call the program instructions to execute the above-mentioned three-dimensional sectioning method of a generalized triangular prism with vertical side edges.

[0034] The advantages of the present invention are:

[0035] (1) The present invention eliminates invalid section planes by comparing the spatial position relationship between all section planes and VGTP units, thereby reducing the actual number of sectioning times of the VGTP unit, thereby reducing the number of newly generated VGTP units at the sectioning position and saving data storage capacity.

[0036] (2) The present invention proposes a strict VGTP unit sectioning method, by which the sectioned VGTP unit is segmented to generate a new VGTP unit, thereby achieving any number of sectionings of the entire model.

[0037] (3) The present invention strictly uses VGTP units to construct the entire model, which is more efficient in terms of spatial position relationship and volume calculation compared to the model constructed using GTP units.

[0038] (4) The present invention specifically determines the VGTP unit to be sectioned through the intersection relationship between all section planes and the VGTP unit, then determines the effective section plane through the inclusion relationship of the sectioned body, and finally divides the sectioned VGTP unit through a series of vertical sections to form a new VGTP unit, and finally achieves the sectioning of the entire three-dimensional geological information model geometry by cyclically executing these steps. The present invention reduces the actual sectioning times of the VGTP unit through the determination of the effective section plane, reduces the model data volume and improves the model operation efficiency, and at the same time, through the strict sectioning method of the VGTP unit, ensures that all units in the model will not mutate into non-VGTP units after sectioning. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 Three types of diagrams of the generalized triangular prism (VGTP) unit with vertical side edges of the present invention.

[0040] Figure 2 Schematic diagram of effective section plane determination in an embodiment of the present invention.

[0041] Figure 3 Schematic diagram of the cutting and segmentation of a VGTP unit with three side edges in an embodiment of the present invention.

[0042] Figure 4 Schematic diagram of the cutting and segmentation of the VGTP unit of two side edges in an embodiment of the present invention.

[0043] Figure 5 Schematic diagram of the cutting and segmentation of a VGTP unit of a side edge in an embodiment of the present invention.

[0044] Description of reference numerals:

[0045] 1- side edge; 2- upper base; 3- lower base; 4- base vertex; 5- base edge; 6- intersection point of real edges. DETAILED DESCRIPTION

[0046] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0047] Example 1

[0048] A model consisting entirely of generalized triangular prism units with vertical side edges is established by interpolation in the drilling data space. In the embodiment, the generalized triangular prism (VGTP) unit with vertical side edges is recorded in the form of (abd-a'b'd'), wherein abd is the three bottom surface vertices 4 of the lower bottom surface 3, and a'b'd' is the three bottom surface vertices 4 of the upper bottom surface 2. In the embodiment, all the bodies, faces and lines except the VGTP unit are represented by all their vertices.

[0049] Figure 1 Figure 1 is a diagram of three types of VGTP units of the present invention, consisting of Figure 1 It can be seen that the VGTP unit of the present invention is surrounded by two upper and lower triangular bases and two / three side surfaces; the upper and lower triangular bases are completely overlapped when projected in the vertical direction, the sides constituting the triangular base are called base edges 5, and the vertices of the base edges 5 are called base vertex 4; the sides of the two to three vertical triangular or quadrilateral side surfaces that intersect and completely overlap in pairs are called side edges 1, and the sides other than the side edges 1 completely overlap with the base edges 5 of the upper and lower triangular bases respectively.

[0050] The data of the VGTP unit is stored as

[0051] Among them, U1, U2, U3 represent the spatial coordinates of the three bottom vertices 4 of the upper bottom surface 2 of the VGTP unit, and D1, D2, D3 represent the spatial coordinates of the three bottom vertices 4 of the lower bottom surface 3 of the VGTP unit;

[0052] Depend on Figure 1 It can be seen that the VGTP unit can be divided into three types according to the number of side edges 1, specifically: if the two side faces of the VGTP unit are two triangular side faces, then the VGTP unit only contains one vertical side edge 1, and the VGTP unit is called a VGTP unit with one side edge 1; if the three side faces of the VGTP unit are two triangular side faces and one quadrilateral side face, then the VGTP unit contains two vertical side edges 1, and the VGTP unit is called a VGTP unit with two side edges 1; if the three side faces of the VGTP unit are all quadrilateral side faces, then the VGTP unit contains three vertical side edges 1, and the VGTP unit is called a VGTP unit with three side edges 1. Figure 1 (1) shows a VGTP unit with three side edges. Figure 1 (2) shows a VGTP unit with two side edges. Figure 1 (3) shows a VGTP unit with a side edge.

[0053] A three-dimensional sectioning method of a generalized triangular prism unit with vertical side edges of the present invention comprises the following steps:

[0054] Step 1: Obtain the VGTP unit to be sectioned

[0055] Obtain a section plane equation αX+βY+χZ+η=0, wherein any spatial point (X′, Y′, Z′) satisfying the section plane equation is located on the section plane; select a VGTP unit and obtain the spatial coordinates of six bottom vertices 4 of its upper bottom surface and lower bottom surface, wherein the spatial coordinates are composed of plane coordinates and vertical elevations;

[0056] Substitute the plane coordinates of the three bottom surface vertices 4 of the upper bottom surface or the lower bottom surface into the section plane equation to obtain three vertical elevations, and form three spatial points with the corresponding plane coordinates, which are recorded as edge intersection points;

[0057] According to the positional relationship between the edge intersection point and the upper and lower bottom surface vertices 4 of the VGTP unit, the edge intersection point is classified, specifically: if the vertical elevation of the edge intersection point is greater than the vertical elevations of the upper and lower bottom surface vertices 4 at the same time, the edge intersection point is recorded as an upper virtual edge intersection point; if the vertical elevation of the edge intersection point is less than the vertical elevations of the upper and lower bottom surface vertices 4 at the same time, the edge intersection point is recorded as a lower virtual edge intersection point; otherwise, the edge intersection point is recorded as a real edge intersection point 6; the edge intersection point and its corresponding upper and lower bottom surface vertices 4 satisfy a positional relationship of coincidence along the vertical projection;

[0058] According to the classification of the edge intersections, it is determined whether the VGTP unit needs to be sectioned. Specifically, if the three edge intersections are all upper virtual edge intersections or all lower virtual edge intersections, it is determined that the VGTP unit does not need to be sectioned; otherwise, it is determined that the VGTP unit needs to be sectioned, and the section surface is recorded accordingly;

[0059] Obtain the equations of the remaining section planes, repeat step 1, and obtain all the VGTP units that need to be sectioned and all the section planes corresponding to each unit;

[0060] There are multiple VGTP units and multiple section planes in space, and the section plane corresponding to each VGTP unit needs to be determined one by one; step S1 is used to determine whether a section plane intersects with a VGTP unit. If so, it is considered that the section plane is the section plane corresponding to the VGTP unit, that is, the VGTP unit needs to be sectioned by the section plane, so step S1 is repeated for each VGTP unit to determine all the section planes corresponding to each VGTP unit;

[0061] Step 2: Determine the effective sectioning plane

[0062] Figure 2 FIG. 1 is a schematic diagram of determining an effective section plane in an embodiment of the present invention. Figure 2 (3) Take step 2 as an example, as follows:

[0063] Obtain a VGTP unit abd-a'b'd' that needs to be sectioned and all its corresponding section planes, and determine a determination method according to the total number n of section planes; Figure 2 (3) The total number of cutting planes n=3>1, all the cutting planes are used to independently cut the VGTP unit, and the volume after all the cutting planes independently cut the VGTP unit is obtained, wherein the volumes A1B1C1b', A2B2C2b' and A3B3C3b' removed after independent cutting are all recorded as removed volumes, and each cutting plane corresponds to a removed volume; the cutting plane has a normal direction perpendicular to the cutting plane, and the normal direction of the cutting plane is any one of the two directions perpendicular to the cutting plane, and the removed volume is not in the direction indicated by the normal direction of the cutting plane;

[0064] Compare the dissected bodies in pairs, wherein the dissected body A3B3C3b' completely includes the dissected bodies A1B1C1b' and A2B2C2b', and define the section planes corresponding to the completely included dissected bodies A1B1C1b' and A2B2C2b' as invalid section planes, and record the completely included dissected bodies A1B1C1b' and A2B2C2b' as dissected body B, and the other dissected body A3B3C3b' compared therewith as dissected body A; the complete inclusion means that the outer contour of the dissected body B does not exceed the outer contour of the dissected body A;

[0065] Repeat step 2 to eliminate all invalid section planes corresponding to the VGTP units that need to be sectioned, and the remaining section planes are recorded as valid section planes;

[0066] Step 3: Unit cutting and segmentation

[0067] Figure 3-5 Schematic diagram of the cutaway and segmentation of three types of VGTP units in the embodiment of the present invention. Figure 3 (4) Figure 4 (6) and Figure 5 (1) as an example, execute step 3, as follows:

[0068] Step 3.1, obtaining a VGTP unit to be sectioned and all corresponding valid section planes; wherein, Figure 3 The VGTP unit of (4) is abd-a'b'd', Figure 4 The VGTP unit of (6) is abd-a'bd', Figure 5 The VGTP unit of (1) is abd-abd';

[0069] Step 3.2, select an effective section plane, calculate the intersection surface between the effective section plane and the VGTP unit, Figure 3 The intersection of (4) is ABCD, Figure 4 The intersection of (6) is ABCDE, Figure 5 The intersection surface of (1) is ABC, and the direction perpendicular to the intersection surface and not pointing to the cutout body is the normal direction of the intersection surface; remove the cutout body, Figure 3 The dissected body of (4) is ABCDbb', Figure 4 The dissected body of (6) is ABCDEa'b, Figure 5 The cut body of (1) is ABCd', and the vertical direction of the normal of the intersection surface is obtained, which is recorded as the projection direction. Figure 3 (4) Figure 4 (6) Figure 5 (1) are all vertically downward; the projection method of the intersection surface is determined according to the projection direction, specifically: the projection direction is vertically downward, and the intersection surface is projected vertically to the lower bottom surface 3 of the original VGTP unit;

[0070] Step 3.3, determining further operations between the projection surface of the intersection surface and the remaining portion of the bottom surface 3 of the VGTP unit after the cutout body is removed, Figure 3 The projection surface of (4) is A'B'CD, Figure 4 The projection surface of (6) is aB'C'DE, Figure 5 The projection surface of (1) is A'dC', Figure 3 The remaining part of the lower bottom surface 3 of (4) is adCD, Figure 4 The remaining part of the lower bottom surface 3 of (6) is adDE, Figure 5 The remaining part of the lower bottom surface 3 of (1) is abd, specifically: the two do not completely overlap, and the edge of the outer contour of the projection surface that does not overlap with the outer contour of the remaining part of the lower bottom surface 3 is recorded as a trace line, Figure 3 The trace of (4) is A'B', Figure 4 The trace of (6) is B'C', Figure 5 The trace of (1) is A'C', and the remaining part of the lower bottom surface 3 is determined as the surface to be divided; the end points of the trace and the vertices of the surface to be divided are both recorded as dividing points;

[0071] Step 3.4, split the face to be split into multiple triangular sub-faces, Figure 3 The triangular sub-faces of (4) include adB', aA'B', A'B'C and A'DC. Figure 4 The triangular sub-faces of (6) include aB'E, B'EC', C'DE and B'dC', Figure 5The triangular sub-faces of (1) include aA'C', aC'b and A'dC', and the vertices of the triangular sub-faces all coincide with the segmentation points; if there is a trace line in the face to be segmented, then the common edges of two triangular sub-faces in the triangular sub-faces coincide with the trace line, Figure 3 In (4), the two triangular sub-faces whose common edges coincide with the trace are aA'B' and A'B'C. Figure 4 In (6), the two triangular sub-faces whose common edges coincide with the trace are B'EC' and B'dC'. Figure 5 (1) The two triangular sub-faces whose common edges coincide with the trace are aA'C' and A'dC';

[0072] Step 3.5, obtain the common edges between all the triangular sub-faces in the face to be split, and record them as splitting lines; Figure 3 The common edges in (4) include aB', A'B' and A'C. Figure 4 The common edges in (6) include B'E, B'C', and C'E; Figure 5 The common edges in (1) include A'C' and aC'; the dividing line is stretched vertically up and down to form a dividing plane; the VGTP unit after the dissection body is removed is divided by the dividing plane to form multiple new VGTP units Figure 3 (4) The new VGTP units formed include adB'-a'd'B, aA'B'-a'AB, A'B'C-ABC and A'DC-ADC. Figure 4 (6) The new VGTP units formed include aB'E-ABE, B'EC'-BEC, C'DE-CDE and B'dC'-Bd'C. Figure 5 (1) The new VGTP unit formed includes aA'C'-aAC, aC'b-aCb and A'dC'-ABC); the upper and lower bottom surfaces of the new VGTP unit need to be determined according to the projection direction in step 3.2, specifically as follows: the projection direction is vertically downward, and the surface of the new VGTP unit that is located in the same spatial surface as the surface to be divided in step 3.3 is determined as the lower bottom surface 3 of the new VGTP unit, and the surface of the new VGTP unit that completely overlaps with the vertical projection of the lower bottom surface 3 of the new VGTP unit is determined as the upper bottom surface 2 of the new VGTP unit; Figure 3 The lower surface 3 of the new VGTP unit in (4) corresponds to adB', aA'B', A'B'C and A'DC, Figure 4 The lower surface 3 of the new VGTP unit in (6) corresponds to aB'E, B'EC', C'DE and B'dC', Figure 5 The lower base 3 of the new VGTP unit in (1) corresponds to aA'C', aC'b and A'dC'; Figure 3 The upper bottom surface 2 of the new VGTP unit in (4) corresponds to a'd'B, a'AB, ABC and ADC. Figure 4The upper base 2 of the new VGTP unit in (6) corresponds to ABE, BEC, CDE and Bd'C. Figure 5 The upper base 2 of the new VGTP unit in (1) corresponds to aAC, aCb and ABC;

[0073] Step 3.6, select a new effective cutting plane, determine the VGTP unit that needs to be cut by the new effective cutting plane in the new VGTP unit described in step 3.5 according to the method described in step 1, and cut and split according to the method described in steps 3.2-3.5; repeat this step until all effective cutting planes complete the cutting and splitting of the new VGTP unit generated by the VGTP unit described in step 3.1 during the cyclic execution of steps 3.2-3.5;

[0074] Repeat step 3 until all VGTP units that need to be sectioned are sectioned and segmented.

[0075] The present invention specifically determines the VGTP unit to be sectioned through the intersection relationship between all section planes and the VGTP unit, then determines the effective section plane through the inclusion relationship of the sectioned body, and finally divides the sectioned VGTP unit through a series of vertical sections to form a new VGTP unit, and finally achieves the sectioning of the entire three-dimensional geological information model geometric body by cyclically executing these steps. The present invention reduces the actual sectioning times of the VGTP unit through the determination of the effective section plane, reduces the model data volume and improves the model operation efficiency, and at the same time, through the strict sectioning method of the VGTP unit, ensures that all units in the model will not mutate into non-VGTP units after sectioning.

[0076] Example 2

[0077] A readable storage medium stores a computer program, wherein the computer program includes program instructions. When the program instructions are executed by a processor, the processor executes a three-dimensional sectioning method of a generalized triangular prism with vertical side edges according to Example 1.

[0078] Example 3

[0079] An electronic device includes a processor, an input device, an output device and a memory, wherein the processor, the input device, the output device and the memory are connected in sequence, the memory is used to store a computer program, the computer program includes program instructions, and the processor is configured to call the program instructions to execute a three-dimensional sectioning method of a generalized triangular prism with vertical side edges according to Example 1.

[0080] The above are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A three-dimensional sectioning method for a generalized triangular prism with vertical side edges, characterized in that: A model consisting entirely of generalized triangular prism units with vertical side edges is established by interpolation in the borehole data space; the generalized triangular prism unit with vertical side edges, i.e., the VGTP unit, is surrounded by two upper and lower triangular bases and two / three side faces; the projections of the upper and lower triangular bases along the vertical direction completely overlap, the sides constituting the triangular bases are called bases (5), and the vertices of the bases (5) are called base vertices (4); the sides that intersect and completely overlap between the side faces are called side edges (1), and the sides on the side faces other than the side edges (1) completely overlap with the bases (5) of the upper and lower triangular bases respectively; The three-dimensional cutting method comprises the following steps: Step 1: Obtain the VGTP unit to be sectioned If the VGTP unit has a section plane that intersects with it, it is determined that the VGTP unit needs to be sectioned, and all the section planes that intersect with it are recorded accordingly; otherwise, the VGTP unit does not need to be sectioned; Step 2: Determine the effective sectioning plane Obtain a VGTP unit that needs to be sectioned and all its corresponding section planes, and determine the judgment method according to the total number n of all the corresponding section planes, specifically: if n=1, the only section plane is a valid section plane; if n>1, let all the section planes independently section the VGTP unit, and obtain the body after all the section planes independently section the VGTP unit, where the body of the VGTP unit that is removed after independent sectioning is recorded as a removed body, and each section plane corresponds to a removed body; compare the removed bodies two by two, if one of the removed bodies completely contains the other removed body, then define the section plane corresponding to the completely contained removed body as an invalid section plane, repeat step 2, eliminate all invalid section planes corresponding to the VGTP units that need to be sectioned, and the remaining section planes are recorded as valid section planes; Step 3: VGTP unit sectioning and segmentation Step 3.1, obtaining a VGTP unit to be sectioned and all corresponding valid section planes; Step 3.2, select a valid cutting plane, calculate the intersection surface of the valid cutting plane and the VGTP unit, record the direction perpendicular to the intersection surface and not pointing to the cut body as the normal of the intersection surface; remove the cut body and obtain the vertical direction of the normal of the intersection surface, record it as the projection direction; determine the projection method of the intersection surface according to the projection direction, specifically: if the projection direction points vertically upward, then project the intersection surface along the vertical direction to the upper bottom surface (2) of the original VGTP unit; if the projection direction points vertically downward, then project the intersection surface along the vertical direction to the lower bottom surface (3) of the original VGTP unit; if the projection direction cannot be obtained, that is, the normal of the intersection surface is completely horizontal, then project the intersection surface along the vertical direction to the upper or lower bottom surface of the original VGTP unit, and define the projection direction as vertically upward or vertically downward accordingly; Step 3.3, according to the relationship between the projection surface of the intersecting surface and the remaining part of the upper or lower bottom surface of the VGTP unit after the dissection body is removed, further operations between the two are determined, specifically: if the two completely overlap, then the completely overlapping surface is determined as the surface to be segmented; if the two do not completely overlap, then the edge of the outer contour of the projection surface that does not overlap with the outer contour of the remaining part of the upper or lower bottom surface is recorded as a trace line, and the remaining part of the upper or lower bottom surface is determined as the surface to be segmented; the endpoints of the trace line and the vertices of the surface to be segmented are both recorded as segmentation points; the selection of the remaining part of the upper or lower bottom surface needs to be determined according to the projection direction, specifically: if the projection direction is vertically downward, it is the remaining part of the lower bottom surface (3); otherwise, it is the remaining part of the upper bottom surface (2); the projection surface of the intersecting surface degenerates into a projection line when the intersecting surface is completely vertical; Step 3.4, split the surface to be split into multiple triangular sub-surfaces, the vertices of the triangular sub-surfaces all coincide with the splitting points; if there is a trace line in the surface to be split, then there are also common edges of two triangular sub-surfaces in the triangular sub-surfaces that coincide with the trace line; Step 3.5, obtain the common edges between all the triangular sub-faces in the face to be segmented, and record them as segmentation lines; stretch the segmentation lines up and down in the vertical direction to form segmentation planes; use the segmentation planes to segment the VGTP units after the dissection body is removed to form multiple new VGTP units; the upper and lower bottom surfaces of the new VGTP units need to be determined according to the projection direction in step 3.2, as follows: if the projection direction is vertically downward, the surface of the new VGTP unit that is located in the same spatial surface as the face to be segmented in step 3.3 is determined as the new VGTP unit The lower bottom surface (3) of the new VGTP unit, the surface in the new VGTP unit that completely overlaps with the vertical projection of the lower bottom surface (3) of the new VGTP unit is determined as the upper bottom surface (2) of the new VGTP unit; if the projection direction is vertically upward, the surface in the new VGTP unit that is located in the same spatial plane as the surface to be divided in step 3.3 is determined as the upper bottom surface (2) of the new VGTP unit, and the surface in the new VGTP unit that completely overlaps with the vertical projection of the upper bottom surface (2) of the new VGTP unit is determined as the lower bottom surface (3) of the new VGTP unit; Step 3.6, select a new effective cutting plane, determine the VGTP unit that needs to be cut by the new effective cutting plane in the new VGTP unit in step 3.5 according to the method of step 1, and cut and split according to the method described in steps 3.2-3.5; repeat step 3.6 until all effective cutting planes complete the cutting and splitting of the new VGTP unit generated by the VGTP unit in step 3.1 during the cyclic execution of steps 3.2-3.5; Repeat step 3 until all VGTP units that need to be sectioned are sectioned and segmented.

2. A three-dimensional sectioning method for a generalized triangular prism with vertical side edges according to claim 1, characterized in that: If the two side faces of a VGTP unit are two triangular side faces, then the VGTP unit contains only one vertical side edge (1), and the VGTP unit is called a VGTP unit with one side edge (1); if the three side faces of a VGTP unit are two triangular side faces and one quadrilateral side face, then the VGTP unit contains two vertical side edges (1), and the VGTP unit is called a VGTP unit with two side edges (1); if the three side faces of a VGTP unit are all quadrilateral side faces, then the VGTP unit contains three vertical side edges (1), and the VGTP unit is called a VGTP unit with three side edges (1).

3. The three-dimensional sectioning method of a generalized triangular prism with vertical side edges according to claim 1, characterized in that: The specific process of step 1 is as follows: Obtain a section plane equation αX+βY+χZ+η=0, wherein any spatial point (X′, Y′, Z′) satisfying the section plane equation is located on the section plane; select a VGTP unit and obtain the spatial coordinates of the six bottom surface vertices (4) of the upper and lower bottom surfaces, wherein the spatial coordinates are composed of plane coordinates and vertical elevations; Substitute the plane coordinates of the three vertices (4) of the upper or lower bottom surface into the section plane equation to obtain three vertical elevations, and form three spatial points with the corresponding plane coordinates, which are recorded as the edge intersection points; According to the positional relationship between the edge intersection point and the upper and lower bottom surface vertices (4) of the VGTP unit, the edge intersection point is classified, specifically: if the vertical elevation of the edge intersection point is greater than the vertical elevations of the upper and lower bottom surface vertices (4), the edge intersection point is recorded as an upper virtual edge intersection point; if the vertical elevation of the edge intersection point is less than the vertical elevations of the upper and lower bottom surface vertices (4), the edge intersection point is recorded as a lower virtual edge intersection point; otherwise, the edge intersection point is recorded as a real edge intersection point (6); the edge intersection point and its corresponding upper or lower bottom surface vertex (4) satisfy a positional relationship of coincidence along the vertical projection; According to the classification of the edge intersections, it is determined whether the VGTP unit needs to be sectioned. Specifically, if the three edge intersections are all upper virtual edge intersections or all lower virtual edge intersections, it is determined that the VGTP unit does not need to be sectioned; otherwise, it is determined that the VGTP unit needs to be sectioned, and the section surface is recorded accordingly; Obtain the equations of the remaining section planes, repeat step 1, and obtain all the VGTP units that need to be sectioned and all the section planes corresponding to each unit.

4. The three-dimensional sectioning method of a generalized triangular prism with vertical side edges according to claim 1, characterized in that: In step 2, the dissected body that is completely contained is recorded as dissected body B, and the other dissected body compared with it is recorded as dissected body A; the complete containment means that the outer contour of the dissected body B does not exceed the outer contour of the dissected body A.

5. The three-dimensional sectioning method of a generalized triangular prism with vertical side edges according to claim 1, characterized in that: In step 2, the normal direction of the section plane is any one of the two directions perpendicular to the section plane, and the sectioned body is not in the direction indicated by the normal direction of the section plane.

6. A readable storage medium, characterized in that: The storage medium stores a computer program, which includes program instructions. When the program instructions are executed by a processor, the processor executes a three-dimensional sectioning method for a generalized triangular prism with vertical side edges as described in any one of claims 1 to 5.

7. An electronic device, characterized in that: The method comprises a processor, an input device, an output device and a memory, wherein the processor, the input device, the output device and the memory are connected in sequence, the memory is used to store a computer program, the computer program comprises program instructions, and the processor is configured to call the program instructions to execute a three-dimensional sectioning method for a generalized triangular prism with vertical side edges as described in any one of claims 1 to 5.

Images