Geological modeling method, device and equipment based on stratigraphic discontinuity pie layering and storage medium

Abstract

The application relates to a stratigraphic modeling method, device and equipment based on stratigraphic interval pie layering and a storage medium, and the method comprises the following steps: obtaining drilling data of a plurality of drill holes, wherein the drilling data of each drill hole comprises drilling position information, inclination information, inclination depth information, elevation information and geological information, the geological information comprises geological era information and main lithology information; determining stratigraphic separation points and stratigraphic interval points of a projection axis of each drill hole based on the drilling position information, the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information; determining stratigraphic separation lines and stratigraphic interval lines between any two adjacent drill holes based on the stratigraphic separation points and the stratigraphic interval points, and generating a two-dimensional stratigraphic pie model between the any two adjacent drill holes; and generating a three-dimensional stratigraphic pie model among any three adjacent drill holes based on the two-dimensional stratigraphic pie model between the any two adjacent drill holes. The application has the effect of improving the drawing rate of geological work.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological exploration, and in particular to a geological modeling method, device and equipment based on stratigraphic discontinuous pie body stratification and a storage medium. BACKGROUND

[0002] In mineral exploration and development, in order to find out the shape of the ore body, thousands of meters of drill holes need to be drilled to explore underground resources, and the stratum geology is inferred according to the data in the drill holes. At present, the geological staff in China mainly analyzes geological problems by drawing geological figures on drawing software. Although some software supports automatic drawing, a large amount of drill hole data needs to be input to generate a geological model. The drill hole data includes drill hole position, hole opening elevation, drill hole depth, depth, thickness, elevation, depth, age, occurrence and lithology of each stratum. The operation amount is large, the modeling period is long, and the efficiency is low. SUMMARY

[0003] In order to improve the modeling efficiency of geological work, the present application provides a geological modeling method, device, equipment and medium based on stratigraphic discontinuous pie body stratification.

[0004] In a first aspect, the present application provides a geological modeling method based on stratigraphic discontinuous pie body stratification, which adopts the following technical scheme: comprising:

[0005] Obtaining drill hole data of a plurality of drill holes, wherein the drill hole data of each drill hole includes drill hole position information, inclination information of each hole section in each drill hole and geological information, and the geological information includes geological age information, inclination depth information, first occurrence information and main lithology information;

[0006] Determining stratigraphic separation points and stratigraphic discontinuous points of a projection axis of each drill hole based on the drill hole position information, the inclination information, the inclination depth information, the elevation information, the geological age information and the main lithology information;

[0007] Determining stratigraphic separation lines and stratigraphic discontinuous lines between any two adjacent drill holes based on the stratigraphic separation points and the stratigraphic discontinuous points, and generating a two-dimensional stratum pie body model between any two adjacent drill holes;

[0008] Based on the two-dimensional stratum pie body model between any two adjacent drill holes, a three-dimensional stratum pie body model between any three adjacent drill holes is generated.

[0009] According to the technical scheme, the electronic device determines the stratigraphic division point and the stratigraphic discontinuity point on the projection axis of each borehole according to the borehole data of the borehole, determines the stratigraphic division line based on the stratigraphic division points between two adjacent boreholes, determines the stratigraphic discontinuity line based on the stratigraphic discontinuity points between two adjacent boreholes, and then generates a two-dimensional stratigraphic pie model between the two adjacent boreholes, and then generates a three-dimensional stratigraphic pie model between any three boreholes according to the two-dimensional stratigraphic pie model between the two adjacent boreholes, which is drawn in a primary and secondary manner and in a zoned and layered manner, improves the drawing rate, generates the three-dimensional stratigraphic pie model according to the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information of the borehole, and does not need the occurrence information, intelligently draws the model by using less borehole data, reduces manual intervention, and then improves the work efficiency and reduces the requirement for computer hardware.

[0010] Further, the determination of the stratigraphic division point and the stratigraphic discontinuity point on the projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information comprises:

[0011] determination of the projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information and the elevation information;

[0012] determination of a plurality of stratigraphic division points based on the geological era information, wherein the geological eras of the stratums on the upper and lower sides of the stratigraphic division point are different;

[0013] determination of a stratigraphic discontinuity point from the plurality of stratigraphic division points based on the geological era information and the main lithology information, wherein the geological eras of the stratums on the upper and lower sides of the stratigraphic discontinuity point are discontinuous;

[0014] determination of the positions of the stratigraphic division point and the stratigraphic discontinuity point on the projection axis based on the inclination information and the inclination depth information.

[0015] According to the technical scheme, the electronic device determines the stratigraphic division point and the stratigraphic discontinuity point according to the characteristics of the geological eras of the stratums, and then determines the positions of the stratigraphic division point and the stratigraphic discontinuity point on the projection axis according to the inclination information and the inclination depth information.

[0016] Further, the stratigraphic discontinuity point comprises a sedimentary discontinuity point, a structural discontinuity point and an intrusive discontinuity point, and the determination of the stratigraphic discontinuity point from the stratigraphic division points based on the geological era information and the main lithology information comprises at least one of the following:

[0017] if the geological eras of the stratums on the upper and lower sides of the stratigraphic division point are discontinuous or discontinuous and reversed, the stratigraphic division point is determined as a sedimentary discontinuity point;

[0018] if the geological era of the strata on both sides of the stratum separation point has a discontinuity and the main lithology information includes structural rock, the stratum separation point is determined as a structural discontinuity point;

[0019] if the geological era of the strata on both sides of the stratum separation point has a discontinuity and the main lithology information includes intrusive rock, the stratum separation point is determined as an intrusive discontinuity point.

[0020] By adopting the above technical solution, the electronic device distinguishes various stratum discontinuity points according to the geological era characteristics and the main lithology characteristics of the strata, facilitates intelligent judgment of the discontinuity points, and improves work efficiency.

[0021] Further, the stratum separation line and the stratum discontinuity line between any two adjacent boreholes are determined based on the stratum separation point and the stratum discontinuity point, and a two-dimensional stratum pie body model between any two adjacent boreholes is generated, including:

[0022] Based on the position information and the depth information of each borehole, a projection axis line for representing the borehole is determined, a homodirectional polyline is connected between the orifices of two boreholes and between the hole bottoms of two boreholes, and a two-dimensional stratum pie body model between two adjacent boreholes is generated, wherein the projection axis line includes the stratum separation point and the stratum discontinuity point;

[0023] Based on the stratum discontinuity points on the two projection axis lines, a stratum discontinuity line is determined, wherein any two adjacent stratum discontinuity lines and two homodirectional polylines form a stratum profile;

[0024] Based on the stratum profile and the stratum separation point, a stratum separation line of each stratum profile is determined.

[0025] By adopting the above technical solution, the electronic device establishes a two-dimensional stratum pie body model between two adjacent boreholes, represents the data of the borehole on the projection axis line, facilitates the electronic device to obtain the stratum discontinuity line according to the corresponding connection of the same stratum discontinuity point on the projection axis line, and further divides the stratum profile according to the stratum discontinuity line, and further determines the stratum separation line in each stratum profile. The electronic device sequentially obtains the stratum discontinuity line and the stratum separation line, improves the drawing speed and accuracy.

[0026] Further, the stratum discontinuity line includes a sedimentary discontinuity line, a structural discontinuity line, and an intrusive discontinuity line, and the determination of the stratum discontinuity line based on the stratum discontinuity points on the two projection axis lines includes at least one of the following:

[0027] determining whether there is a same sedimentary discontinuity point on the two projection axis lines, if yes, connecting the two same sedimentary discontinuity points to obtain a sedimentary discontinuity line; otherwise, determining an estimated point on the nearest same sense fold line or stratigraphic discontinuity line, connecting the estimated point and the sedimentary discontinuity point to obtain a sedimentary discontinuity line, and the sedimentary discontinuity line is parallel to the adjacent stratigraphic layer;

[0028] determining whether there is a same intrusion discontinuity point on the two projection axis lines, if yes, connecting the two same intrusion discontinuity points to obtain an intrusion discontinuity line; otherwise, obtaining first occurrence information at the intrusion discontinuity point, determining a slope of the intrusion discontinuity line based on the first occurrence information, and determining the intrusion discontinuity line between the projection axis line and the nearest stratigraphic discontinuity line based on the intrusion discontinuity point and the slope;

[0029] determining whether there is a same structural discontinuity point on the two projection axis lines, if yes, connecting the two same structural discontinuity points to obtain a structural discontinuity line; otherwise, determining a slope of the intrusion discontinuity line based on occurrence information at the structural discontinuity point, determining an estimated point on the same sense fold line based on the structural discontinuity point and the slope, and connecting the estimated point and the structural discontinuity point to obtain a structural discontinuity line.

[0030] By adopting the above technical solutions, when there is a same stratigraphic discontinuity point on the two projection axis lines, the electronic device directly connects the two stratigraphic discontinuity points to determine a stratigraphic discontinuity line; when there is no stratigraphic discontinuity line on one of the projection axis lines, the electronic device determines an estimated point according to different methods corresponding to various types of stratigraphic discontinuity lines, and determines a stratigraphic discontinuity line according to the stratigraphic discontinuity point and the estimated point, so that the electronic device can generate a stratigraphic discontinuity line of a corresponding type for various situations.

[0031] Further, the determination of the separation line of each stratigraphic section based on the stratigraphic discontinuity surface and the separation point comprises:

[0032] determining whether there is a stratigraphic separation point on a line segment between the two adjacent stratigraphic discontinuity lines; if yes, determining that the stratigraphic section is a multi-layer section; if not, determining that the stratigraphic section is a single-layer section;

[0033] if the stratigraphic section is a multi-layer section, determining a main layer and a interlayer based on the stratigraphic separation point and the geological era information;

[0034] determining the type of each main layer based on the stratigraphic separation point corresponding to each main layer, comprising at least one of the following:

[0035] a connected type stratigraphic layer: there are two stratigraphic separation points of the main layer on the two projection axis lines;

[0036] Half-continuous stratum: both of the projection axis exist the stratum separation point of one side of the main layer, and only one of the projection axis exists the stratum separation point of the other side of the main layer;

[0037] Intermittent stratum: one of the projection axis exists the stratum separation point of both sides of the main layer and does not coincide, and the other of the projection axis exists the stratum separation point of both sides of the main layer and coincides;

[0038] Half-intermittent stratum: only one of the projection axis exists the stratum separation point of both sides of the main layer, and the two stratum separation points do not coincide;

[0039] Based on the type of the main layer, the stratum separation line of both sides of the main layer is determined, comprising:

[0040] Connecting the two stratum separation points of both sides of the main layer of the type of continuous stratum respectively, the stratum separation line of both sides of the main layer is obtained;

[0041] Starting from the continuous stratum, the corresponding stratum separation line of the remaining main layer is determined layer by layer upward and downward respectively, including any one of the following:

[0042] If the main layer is a half-continuous stratum, connecting the two stratum separation points of one side of the main layer, the stratum separation line of one side of the main layer is obtained, determining an estimated point on the stratum discontinuous line located on the other side of the main layer, connecting the stratum separation point of the other side of the main layer and the estimated point, the stratum separation line of the other side of the main layer is obtained;

[0043] If the main layer is an intermittent stratum, connecting the two stratum separation points of both sides of the main layer respectively, the stratum separation line of both sides of the main layer is obtained;

[0044] If the main layer is a half-intermittent stratum and the main layer is not the last stratum, determining an estimated point on the stratum discontinuous line closest to the main layer, connecting the stratum separation point of one side of the main layer and the adjacent estimated point, the stratum separation line of one side of the main layer is obtained, connecting the stratum separation point of the other side of the main layer and the estimated point, the stratum separation line of the other side of the main layer is obtained, and the two stratum separation lines of both sides of the main layer are parallel;

[0045] If the main layer is a half-intermittent stratum and the main layer is the last stratum, connecting the stratum separation point of the main layer with the adjacent estimated point respectively, the stratum separation line of both sides of the main layer is obtained;

[0046] Based on the stratum separation line of both sides of the main layer and the corresponding stratum separation point of the interlayer, the stratum separation line of both sides of the interlayer is determined.

[0047] By adopting the technical scheme, after the electronic device determines that the stratigraphic profile is a multi-layer profile, the type of each main layer is first determined, the stratigraphic separation lines on both sides of the main layer are determined according to the type of the main layer, and then the stratigraphic separation lines on both sides of the interlayer are determined according to the stratigraphic separation lines of the main layer and the stratigraphic separation points of the interlayer. The electronic device analyzes the strata in a primary and secondary manner, achieves the effect of rapid drawing, and reduces the error rate.

[0048] Further, the determination of the stratigraphic separation lines on both sides of the interlayer based on the stratigraphic separation lines on both sides of the main layer and the stratigraphic separation points corresponding to the interlayer comprises at least one of the following:

[0049] If the stratigraphic separation points on both sides of the interlayer exist on both the projection axes, the stratigraphic separation lines on both sides of the interlayer are obtained by connecting the stratigraphic separation points on both sides of the interlayer, respectively;

[0050] If the stratigraphic separation points on one side of the interlayer exist on both the projection axes, and the stratigraphic separation points on the other side of the interlayer exist on only one of the projection axes, the stratigraphic separation lines on both sides of the interlayer are obtained by connecting the two stratigraphic separation points on one side of the interlayer, determining an estimated point on the stratigraphic separation line of the main layer on the other side of the interlayer, and connecting the estimated point and the stratigraphic separation point on the other side of the interlayer.

[0051] If the stratigraphic separation points on both sides of the interlayer exist on one of the projection axes, and the stratigraphic separation points on both sides of the interlayer do not exist on the other projection axis, the stratigraphic separation lines on both sides of the interlayer are obtained by determining an interval point between the two stratigraphic separation lines of the main layer and the two projection axes, and connecting the stratigraphic separation points and the interval point.

[0052] By adopting the technical scheme, the electronic device generates the stratigraphic separation points on both sides of the interlayer by using a corresponding generation method according to the conditions of the stratigraphic separation points corresponding to the interlayer on the two projection axes.

[0053] Further, the generation of the three-dimensional stratigraphic pie model among any three adjacent boreholes based on the two-dimensional stratigraphic pie model between any two adjacent boreholes comprises:

[0054] It is determined whether there are at least two groups of same stratigraphic separation lines or stratigraphic discontinuity lines in the two-dimensional stratigraphic pie model among any three adjacent boreholes. If yes, the same stratigraphic separation lines are connected in a head-to-tail manner, and a stratigraphic separation surface is obtained by point method equation calculation. The same stratigraphic discontinuity lines are connected in a head-to-tail manner, and a stratigraphic discontinuity surface is obtained by point method equation calculation.

[0055] Otherwise, a stratigraphic discontinuity surface is determined based on the stratigraphic discontinuity lines, wherein the stratigraphic discontinuity surface comprises a sedimentation discontinuity surface, a structure discontinuity surface, and an intrusion discontinuity surface, and at least one of the following is included:

[0056] determining a deposition estimation point between any three adjacent boreholes, the deposition estimation point being located on the stratigraphic boundary surface of the nearest main layer, connecting both ends of the deposition discontinuity line with the deposition estimation point, and calculating the deposition discontinuity surface based on the point equation;

[0057] determining a structure discontinuity estimation point at any position between the three adjacent boreholes, connecting both ends of the structure discontinuity line with the structure discontinuity estimation point, and calculating the structure discontinuity surface based on the point equation;

[0058] determining an intrusion discontinuity estimation point at any position between the three adjacent boreholes, connecting both ends of the intrusion discontinuity line with the intrusion discontinuity estimation point, and calculating the intrusion discontinuity surface based on the point equation;

[0059] determining the stratigraphic boundary surface based on the stratigraphic boundary line, including at least one of the following:

[0060] determining at least one main layer estimation point on the stratigraphic boundary surface on both sides of the main layer, connecting both ends of the stratigraphic boundary line on both sides of the main layer with the main layer estimation point on the same side, and calculating the stratigraphic boundary surface on both sides of the main layer based on the point equation;

[0061] determining at least one interlayer estimation point between any three adjacent boreholes, the interlayer estimation point being located between the two stratigraphic boundary surfaces on both sides of the corresponding main layer, connecting both ends of the stratigraphic boundary line on both sides of the interlayer with the interlayer estimation point, and calculating the stratigraphic boundary surface on both sides of the interlayer based on the point equation;

[0062] determining the stratigraphic pie of the discontinuity layer based on the stratigraphic profile of the discontinuity layer corresponding to the two adjacent stratigraphic boundary surfaces and any two adjacent boreholes;

[0063] determining the stratigraphic pie of the main layer based on the stratigraphic profile of the main layer corresponding to any two adjacent boreholes and the stratigraphic boundary surface on both sides of the same main layer;

[0064] determining the stratigraphic pie of the interlayer based on the stratigraphic profile of the interlayer corresponding to any two adjacent boreholes and the stratigraphic boundary surface on both sides of the same interlayer;

[0065] combining the stratigraphic pie of the discontinuity layer, the stratigraphic pie of the main layer, and the stratigraphic pie of the interlayer to generate a three-dimensional stratigraphic pie model between any three adjacent boreholes;

[0066] assigning different colors to any two adjacent stratigraphic pies.

[0067] By adopting the technical scheme, the electronic device connects the same stratum separating lines between any two of the three adjacent boreholes according to the line connection plane method, determines the contour, and determines the equation of the stratum separating plane according to the point method equation. The same principle is used to further obtain the stratum separating planes or the stratum discontinuity planes on both sides of each stratum. If only one group of two adjacent boreholes has a stratum separating line among the three adjacent boreholes, the electronic device estimates the stratum separating lines on both sides of the stratum according to the type of the stratum, and further encloses a three-dimensional stratum pie body according to the stratum separating plane and the stratum profile.

[0068] In another possible implementation, the method further includes:

[0069] obtaining two points determined by a user above the three-dimensional stratum pie body model;

[0070] determining a tangent line connected to the two points;

[0071] based on the tangent line, longitudinally cutting the three-dimensional stratum pie body model in a vertical direction to generate a profile cut by a graph.

[0072] By adopting the technical scheme, when the user needs a profile cut by a graph, the electronic device determines a tangent line according to two points selected by the user, and further determines a profile cut by a graph according to the tangent line, thereby enriching the use function of the three-dimensional stratum pie body model.

[0073] In a second aspect, the application provides a geological modeling device based on stratum discontinuity pie body layering, including:

[0074] an obtaining module, configured to obtain borehole data of a plurality of boreholes, wherein the borehole data of each borehole includes borehole position information, inclination information, inclination depth information, elevation information, and geological information, and the geological information includes geological era information and main lithology information;

[0075] a determining module, configured to determine stratum separating points and stratum discontinuity points of a projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information, the elevation information, the geological era information, and the main lithology information; and a first generating module, configured to determine stratum separating lines and stratum discontinuity lines between any two adjacent boreholes based on the stratum separating points and the stratum discontinuity points, and generate a two-dimensional stratum pie body model between the any two adjacent boreholes;

[0076] a second generating module, configured to generate a three-dimensional stratum pie body model between any three adjacent boreholes based on the two-dimensional stratum pie body model between the any two adjacent boreholes.

[0077] According to the technical scheme, the determining module determines the stratigraphic division points and the stratigraphic discontinuity points on the projection axis of each borehole according to the borehole data information obtained by the obtaining module, the first generating module determines the stratigraphic division line based on the stratigraphic division points between two adjacent boreholes and determines the stratigraphic discontinuity line based on the stratigraphic discontinuity points between the two adjacent boreholes, and then generates a two-dimensional stratigraphic pie body model between the two adjacent boreholes, and then the second generating module generates a three-dimensional stratigraphic pie body model between any three boreholes according to the two-dimensional stratigraphic pie body model between the two adjacent boreholes, which has a main part and a secondary part and is drawn in a zoned and layered manner, improves the drawing rate, generates the three-dimensional stratigraphic pie body model according to the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information of the borehole without the occurrence information, intelligently draws the model by using less borehole data, reduces manual intervention, improves the work efficiency, and reduces the requirement for computer hardware.

[0078] In a third aspect, the present application provides an electronic device, which adopts the following technical scheme:

[0079] at least one processor;

[0080] a memory;

[0081] at least one application program, wherein the at least one application program is stored in the memory and is configured to be executed by the at least one processor, and the at least one application program is configured to execute the method according to any one of the first aspect.

[0082] According to the technical scheme, the processor executes the application program in the memory, determines the stratigraphic division points and the stratigraphic discontinuity points on each borehole according to the borehole data of the borehole, determines the stratigraphic division line based on the stratigraphic division points between two adjacent boreholes, determines the stratigraphic discontinuity line based on the stratigraphic discontinuity points between the two adjacent boreholes, and then generates a two-dimensional stratigraphic pie body model between the two adjacent boreholes, and then generates a three-dimensional stratigraphic pie body model between any three boreholes according to the two-dimensional stratigraphic pie body model between the two adjacent boreholes, which has a main part and a secondary part and is drawn in a zoned and layered manner, improves the drawing rate, generates the three-dimensional stratigraphic pie body model according to the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information of the borehole without the occurrence information, intelligently draws the model by using less borehole data, reduces manual intervention, improves the work efficiency, and reduces the requirement for computer hardware.

[0083] In a fourth aspect, the present application provides a computer readable storage medium, which stores a computer program capable of being loaded and executed by a processor to execute the method according to any one of the first aspect.

[0084] By adopting the technical scheme, the processor executes the computer program in the computer readable storage medium, determines the stratigraphic division point and the stratigraphic discontinuity point on each borehole according to the borehole data of the borehole, determines the stratigraphic division line based on the stratigraphic division points between two adjacent boreholes, determines the stratigraphic discontinuity line based on the stratigraphic discontinuity points between two adjacent boreholes, and then generates a two-dimensional stratigraphic pie body model between the two adjacent boreholes, and then generates a three-dimensional stratigraphic pie body model between any three boreholes according to the two-dimensional stratigraphic pie body model between the two adjacent boreholes, which has a main and a secondary, a partition and a layer, and improves the drawing rate. According to the dip angle information, the oblique depth information, the elevation information, the geological era information and the main lithology information of the borehole, the three-dimensional stratigraphic pie body model can be generated without the occurrence information, the intelligent drawing can be realized by using less borehole data, the manual intervention is reduced, the working efficiency is improved, and the requirement for computer hardware is reduced.

[0085] In summary, the present application has at least one of the following beneficial technical effects:

[0086] 1. Determine the stratigraphic division point and the stratigraphic discontinuity point on each borehole according to the borehole data of the borehole, determine the stratigraphic division line based on the stratigraphic division points between two adjacent boreholes, determine the stratigraphic discontinuity line based on the stratigraphic discontinuity points between two adjacent boreholes, and then generate a two-dimensional stratigraphic pie body model between the two adjacent boreholes, and then generate a three-dimensional stratigraphic pie body model between any three boreholes according to the two-dimensional stratigraphic pie body model between the two adjacent boreholes, which intelligently improves the modeling efficiency of geological work and improves the working efficiency.

[0087] 2. After the electronic device determines that the stratigraphic profile is a multi-layer profile, first determine the type of each main layer, determine the stratigraphic division line on both sides of the main layer according to the type of the main layer, and then determine the stratigraphic division line on both sides of the interlayer according to the stratigraphic division line of the main layer and the stratigraphic division point of the interlayer. The electronic device analyzes the strata in a primary and secondary manner, achieves rapid drawing, and reduces the error rate. BRIEF DESCRIPTION OF DRAWINGS

[0088] Figure 1 is a schematic diagram of a same-direction fold line in an embodiment of the present application.

[0089] Figure 2 is a flowchart of a geological modeling method based on stratigraphic discontinuity pie layering in an embodiment of the present application.

[0090] Figure 3 is a stratigraphic structure diagram between two adjacent boreholes in an embodiment of the present application.

[0091] Figure 4 is a structural schematic diagram of a multi-layer profile in an embodiment of the present application.

[0092] Figure 5 is a structural schematic diagram of a connected stratigraphic layer in an embodiment of the present application.

[0093] Figure 6 is a structural schematic diagram of the continuous formation (upper) and the intermissive formation (lower) in the embodiment of the present application.

[0094] Figure 7 is a structural schematic diagram of the semi-intermissive formation (upper) and the semi-continuous formation (lower) in the embodiment of the present application.

[0095] Figure 8 is a structural schematic diagram of the interlayer in the main layer in the embodiment of the present application.

[0096] Figure 9 is a structural schematic diagram of the formation separation surface in the three-dimensional formation cake model in the embodiment of the present application.

[0097] Figure 10 is a structural schematic diagram of the three-dimensional formation cake model in the embodiment of the present application.

[0098] Figure 11 is a structural block diagram of the geological modeling device based on the formation intermissive cake layering in the embodiment of the present application.

[0099] Figure 12 is a structural block diagram of the electronic device in the embodiment of the present application. DETAILED DESCRIPTION

[0100] In order to make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0101] In addition, the term "and / or" herein is only to describe the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B can mean that A exists alone, A and B exist together, and B exists alone. In addition, the character " / " herein generally represents an "or" relationship between the associated objects unless otherwise specified.

[0102] The geological modeling method based on the formation intermissive cake layering, referred to as "intermissive cake layering method", is a method for establishing two-dimensional and three-dimensional formation cake models based on geological engineering such as drilling.

[0103] Among them, the occurrence is the general term of the spatial output state and the direction of the formation, including the strike, the tendency and the dip angle three elements.

[0104] Geological discontinuity refers to a stratigraphic separation surface composed of sedimentary discontinuity, structural discontinuity and intrusive discontinuity under the joint influence of geological era and main lithology of stratum.

[0105] Cake body refers to a geological structure similar to a cake divided from stratum. In three-dimensional space, the cake body is a solid three-dimensional geological model surrounded by a top surface, a bottom surface and side surfaces. In two-dimensional plane, the vertical section of the cake body is composed of four same-direction fold lines.

[0106] Same-direction fold line: the angle between the vector on any line segment of the fold line and the vector from the start point to the end point of the fold line is less than 90 degrees. For reference Figure 1 , the fold line AB in the figure is a same-direction fold line.

[0107] The embodiment of the present application discloses a geological modeling method based on stratigraphic discontinuity cake body layering. For reference Figure 1 , the electronic device is executed, including (steps S101-S104):

[0108] Step S101: Obtain drilling data of multiple drill holes, wherein the drilling data of each drill hole includes drilling position information, inclination information, inclination depth information, elevation information and geological information, and the geological information includes geological era information and main lithology information.

[0109] Specifically, the drilling data editing entrance and the stratum data editing entrance are set in the operation interface, and the user can input the drilling data in the form of a table through the above two entrances. The electronic device obtains the drilling data of each drill hole input by the user, and can store multiple drilling data.

[0110] Each drill hole includes several hole sections, and the inclination angles of the hole sections are different, and the multiple hole sections constitute a same-direction fold line. The drilling position information is the longitude and latitude of the drill hole, which is convenient for the electronic device to determine the position of the drill hole in the coordinate system; the inclination information is the angle between each hole section and the vertical axis, the inclination depth information is the drilling depth along the hole section direction, and the inclination depth information also includes the inclination depth information of each stratum, and the elevation information is the height difference between the drill hole and the reference surface.

[0111] Step S102: Determine the stratigraphic separation points and stratigraphic discontinuity points of the projection axis of each drill hole based on the drilling position information, the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information, including (steps S1021-S1023):

[0112] Step S1021: Determine multiple stratigraphic separation points based on the geological era information, wherein the geological eras of the strata on the upper and lower sides of the stratigraphic separation points are different.

[0113] Specifically, the stratigraphic separation point is approximately a point where a stratigraphic separation surface intersects the borehole, and the stratigraphic separation point includes a stratigraphic discontinuity point.

[0114] Further, from the direction of deepening into the stratum from the ground, the side of the stratigraphic separation point close to the ground is the upper side, and the side of the stratigraphic separation point away from the ground is the lower side. When the electronic device determines that the geological ages of the two sides of a point are different, the point is determined as a stratigraphic separation point, that is, the two sides of the stratigraphic separation point are two different strata.

[0115] Step S1022: determining a stratigraphic discontinuity point from the plurality of stratigraphic separation points based on the geological age information and the main lithology information, wherein the geological ages of the strata on the upper and lower sides of the stratigraphic discontinuity point are discontinuous.

[0116] Specifically, the stratigraphic discontinuity point is approximately a point where a stratigraphic discontinuity surface intersects the borehole.

[0117] When the electronic device determines that the geological ages of the two sides of a stratigraphic separation point are discontinuous, there is a missing stratum or an increased stratum, and the current stratigraphic separation point is determined as a stratigraphic discontinuity point, and the type of the stratigraphic discontinuity point is determined according to the main lithology.

[0118] The stratigraphic discontinuity point includes a sedimentary discontinuity point, a structural discontinuity point, and an intrusive discontinuity point.

[0119] Step S1022 specifically includes at least one of the following (step S11 to step S13):

[0120] Step S11: If the geological ages of the strata on the upper and lower sides of the stratigraphic separation point have discontinuity or discontinuity and inversion, the stratigraphic separation point is determined as a sedimentary discontinuity point.

[0121] Specifically, the sedimentary discontinuity refers to discontinuous deposition. There is a widespread denudation surface contact relationship between the new and old strata. The lithology of the strata above and below the denudation surface suddenly changes and forms an irregular undulating surface, indicating that after the deposition of the lower stratum, it was subjected to long-term erosion, and then the deposition of the upper new stratum caused the missing of some strata.

[0122] Therefore, when the electronic device determines that the geological ages of the two sides of the stratigraphic separation point have discontinuity, the stratigraphic separation point is determined as a sedimentary discontinuity point.

[0123] In another case, when the stratum is folded due to external force, the ages of the two adjacent strata in the borehole will be discontinuous and inverted. Therefore, when the geological ages of the two sides of the stratigraphic separation point have discontinuity and inversion, the electronic device determines the stratigraphic separation point as a sedimentary discontinuity point.

[0124] Step S12: If the geological ages of the strata on the upper and lower sides of the stratigraphic separation point have discontinuity and the main lithology information includes structural rocks, the stratigraphic separation point is determined as a structural discontinuity point.

[0125] Specifically, the extrusion force or tension of the earth's crust causes the relative displacement of the rock mass on both sides of the fracture surface, so that the geological age of the stratum along the drilling direction of the borehole is reversed, which indicates that the stratum separation point is a structural discontinuity point.

[0126] In addition, tectonic rock, also known as dynamic metamorphic rock, is a kind of rock formed by the crushing, deformation and recrystallization of the original rock in the fracture zone under the action of different stress. Therefore, the electronic device can determine whether the stratum separation point is a structural discontinuity point according to the main lithology information of the stratum.

[0127] Step S13: If there is a discontinuity in the geological age of the stratum on both sides of the stratum separation point and the main lithology information includes intrusive rock, it is determined that the stratum separation point is an intrusive discontinuity point.

[0128] Specifically, intrusive rock refers to the rock formed by the magma in the asthenosphere drilling out when the overlying rock pressure is reduced and condensing at a deep place in the earth's crust.

[0129] In another possible implementation, when the user inputs the borehole data, if it is determined through the geological information of the stratum that the stratum discontinuity point, the user can mark the stratum discontinuity point, and the electronic device automatically determines the type of the stratum discontinuity point according to the artificial marking information. For example, the user marks the current stratum separation point as a sedimentary discontinuity point, and the electronic device determines that the current stratum separation point is a sedimentary separation point.

[0130] Step S1023: Determine the positions of the stratum separation points and the stratum discontinuity points on the projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information and the elevation information.

[0131] Specifically, each borehole is composed of a plurality of hole segments connected end to end, and the inclination of each hole segment is different. Therefore, in order to facilitate the generation of a two-dimensional stratum pie body between two boreholes, it is necessary to convert the borehole into a straight line, preferably a vertical straight line, to facilitate the determination of a plane according to two adjacent boreholes.

[0132] The electronic device determines the position of the borehole mouth in the coordinate system according to the borehole position information and the elevation information, establishes a vertical straight line from the mouth, and determines the projection of the stratum separation points and the stratum discontinuity points on the straight line according to the inclination information and the inclination depth information of the borehole, and further generates a projection axis, i.e. the projection axis includes the stratum separation points and the stratum discontinuity points.

[0133] Step S103: Determine the stratum separation line and the stratum discontinuity line between any two adjacent boreholes based on the stratum separation points and the stratum discontinuity points, and generate a two-dimensional stratum pie body model between any two adjacent boreholes.

[0134] Specifically, refer to Figure 3 and Figure 4, step S103 includes (step S1031~step S1033):

[0135] Step S1031: connecting a same direction fold line between the top of the projection axis of any two adjacent boreholes and the bottom of the two projection axes to generate a two-dimensional strata pie model between the two adjacent boreholes.

[0136] For example Figure 4 In the figure, the projection axis AD and the projection axis BC are the projection axes corresponding to the boreholes, wherein points 3, 4, 5, 6, 7 and 8 are the strata separation points and strata discontinuity points on one of the projection axes, and points 1', 2', 3', 4', 5', 6' and 7' are the strata separation points and strata discontinuity points on the other projection axis.

[0137] Step S1032: determining a strata discontinuity line based on the strata discontinuity points on the two projection axes, wherein any two adjacent strata discontinuity lines and the two same direction fold lines form a strata profile.

[0138] The strata discontinuity line includes a sedimentary discontinuity line, a structural discontinuity line and an intrusive discontinuity line.

[0139] Referring to Figure 3 , in the figure: A: structural discontinuity line; B: main strata layer; C: multi-layer profile; D: strata interlayer; E: single-layer profile; F: sedimentary discontinuity line; G: intrusive discontinuity line.

[0140] The strata pie includes a multi-layer discontinuity body and a single-layer discontinuity body, and when corresponding to the two-dimensional strata pie model, the multi-layer discontinuity body appears as a multi-layer profile and the single-layer discontinuity body appears as a single-layer profile.

[0141] The structural discontinuity surface appears as a structural discontinuity line in the two-dimensional strata pie model; the sedimentary discontinuity surface appears as a sedimentary discontinuity line in the two-dimensional strata pie model; and the intrusive discontinuity surface appears as an intrusive discontinuity line in the two-dimensional strata pie model. The lines on both sides of the main layer and the interlayer in the figure are strata separation lines.

[0142] Step S1032 includes (step Sa~step Sc) when determining the strata discontinuity line based on the strata discontinuity points on the two projection axes.

[0143] Step Sa: determining whether there are the same sedimentary discontinuity points on the two projection axes, if yes, connecting the two same sedimentary discontinuity points to obtain a sedimentary discontinuity line; otherwise, determining an estimated point on the nearest same direction fold line or strata discontinuity line, connecting the estimated point and the sedimentary discontinuity point to obtain a sedimentary discontinuity line, and the sedimentary discontinuity line is parallel to the adjacent strata.

[0144] Specifically, when there is a deposition discontinuity point on one of the projection axes and there is no deposition discontinuity point on the other projection axis, if the closest to the current deposition discontinuity point is a certain same-direction fold line, the estimation point is determined on the closest same-direction fold line, and the generated deposition discontinuity line is made parallel to the adjacent stratum.

[0145] Step Sb: determine whether there is a same intrusion discontinuity point on the two projection axes, if yes, connect the two same intrusion discontinuity points to obtain an intrusion discontinuity line; otherwise, obtain first occurrence information at the intrusion discontinuity point, determine the slope of the intrusion discontinuity line based on the first occurrence information, and determine the intrusion discontinuity line between the projection axis and the nearest stratum discontinuity line based on the intrusion discontinuity point and the slope.

[0146] Specifically, when there is an intrusion discontinuity point on one of the projection axes and there is no same intrusion discontinuity point on the other projection axis, the first occurrence information at the intrusion discontinuity point can represent the occurrence of the intrusive stratum at the borehole, so the electronic device determines the slope of the intrusion discontinuity line according to the first occurrence information, generates the intrusion discontinuity line according to the intrusion discontinuity point and the slope, and makes the intrusion discontinuity line be truncated when it meets the stratum discontinuity line.

[0147] The electronic device obtains the first occurrence information when determining the intrusion discontinuity line, which improves the accuracy of drawing the intrusion discontinuity line. When drawing other positions, the electronic device usually uses intelligent deduction, and only needs the first occurrence information when drawing the intrusion discontinuity line, which reduces the amount of calculation and improves the drawing rate.

[0148] Step Sc: determine whether there is a same structural discontinuity point on the two projection axes, if yes, connect the two same structural discontinuity points to obtain a structural discontinuity line; otherwise, determine the slope of the intrusion discontinuity line based on the occurrence information at the structural discontinuity point, determine the estimation point on the same-direction fold line based on the structural discontinuity point and the slope, and connect the estimation point and the structural discontinuity point to obtain the structural discontinuity line.

[0149] Specifically, when there is a structural discontinuity point on one of the projection axes and there is no same structural discontinuity point on the other projection axis, the other end of the structural discontinuity line may be truncated on the same-direction fold line on one side. The first occurrence information at the structural discontinuity point can represent the occurrence of the structural discontinuity surface at the borehole, so the electronic device determines the slope of the structural discontinuity line according to the first occurrence information, determines the inclination direction of the structural discontinuity line through the slope, and then determines the same-direction fold line on which side to determine the estimation point, and then determines the structural discontinuity line.

[0150] Step S1033: determining a stratum separation line of each stratum profile based on the stratum profile and the stratum separation point, including (step S21 to step S25):

[0151] Step S21: judging whether there is a stratigraphic separation point on the line segment between the two adjacent stratigraphic discontinuity lines of the two projection axes; if yes, determining that the stratigraphic profile is a multi-layer profile; if no, determining that the stratigraphic profile is a single-layer profile.

[0152] Specifically, referring to Figure 3 When the stratigraphic profile is a single-layer profile, the stratigraphic profile is a single layer. When the stratigraphic profile is a multi-layer profile, i.e., the multi-layer profile includes multiple layers, the electronic device determines the stratigraphic separation line in the multi-layer profile, further refines the two-dimensional stratigraphic pie model, and executes steps S22-S25.

[0153] Step S22: if the stratigraphic profile is a multi-layer profile, determining the main layer and the interlayer based on the geological era information and the stratigraphic separation point.

[0154] Specifically, the electronic device pre-stores the geological era information of the layers of different geological eras, and also includes the geological era information of the interlayer in each main layer. Therefore, the electronic device can determine whether the two adjacent stratigraphic separation points are main layers or interlayers based on the geological era information.

[0155] In order to quickly and orderly generate a two-dimensional stratigraphic pie model, the electronic device first generates the stratigraphic separation line on both sides of the main layer, and then generates the stratigraphic separation line on both sides of the interlayer in each main layer.

[0156] Step S23: determining the type of each main layer based on the stratigraphic separation point corresponding to each main layer, including at least one of the following:

[0157] Connected type layer: there are two stratigraphic separation points of the main layer on both projection axes;

[0158] Semi-connected type layer: there are stratigraphic separation points of the main layer on both projection axes on one side, and there are stratigraphic separation points of the main layer on the other side on only one projection axis;

[0159] Intermittent type layer: there are stratigraphic separation points of the main layer on both sides on one projection axis and do not coincide, and there are stratigraphic separation points of the main layer on both sides on the other projection axis and coincide;

[0160] Semi-intermittent type layer: there are stratigraphic separation points of the main layer on both sides on only one projection axis, and the two stratigraphic separation points do not coincide.

[0161] Referring to Figures 5 to 7 , the line segment AB and the line segment CD are two projection axes, the line segment AC and the line segment BD are two same direction polyline, E and F are both stratigraphic separation points, and EF can be a stratigraphic separation line or a stratigraphic discontinuity line.

[0162] Referring to Figure 5 , the two main layers above and below the line segment EF are connected type.

[0163] Referring to Figure 6 , the main layer above the line segment EF is of the connected type; and the main layer below the line segment EF is of the isolated type.

[0164] Referring to Figure 7 , the main layer above the line segment EF is of the semi-isolated type; and the main layer below the line segment EF is of the semi-connected type.

[0165] Step S24: determining the stratigraphic separation line on both sides of the main layer based on the type of the main layer. Referring to Figure 4 , the line segment AD and the line segment BC are respectively a projection axis, and the line segment AB and the line segment CD are a stratigraphic discontinuity line, then step S24 includes (step S241 to step S242):

[0166] Step S241: connecting two stratigraphic separation points on both sides of the main layer of the type of the connected stratigraphic layer to obtain the stratigraphic separation line on both sides of the main layer.

[0167] Specifically, the main layer ⑤⑥⑦ is of the connected type, and both ends are completely connected on both sides of the approximately vertical projection axis. Connecting the stratigraphic separation points 4-4', 5-5', 6-6', and 7-7' on both sides can obtain the second occurrence information of the main layer through the same layer interconnection.

[0168] Step S242: determining the stratigraphic separation line corresponding to the remaining main layer from the connected stratigraphic layer, respectively, and simultaneously upwards and downwards layer by layer.

[0169] Specifically, since the connected stratigraphic layer can be uniquely determined, the accuracy is relatively high, and therefore the electronic device determines upwards and downwards layer by layer from the connected stratigraphic layer, calculates the stratigraphic separation line that cannot be directly obtained according to the trend of the connected stratigraphic layer, and further improves the accuracy of the two-dimensional stratigraphic pie model.

[0170] Step S242 includes any one of the following:

[0171] If the main layer is of the semi-connected type, connecting two stratigraphic separation points on one side of the main layer to obtain the stratigraphic separation line on one side of the main layer, determining an estimated point on the stratigraphic discontinuity line located on the other side of the main layer, and connecting the stratigraphic separation point on the other side of the main layer and the estimated point to obtain the stratigraphic separation line on the other side of the main layer.

[0172] Specifically, the main layer ②⑧ is of the semi-connected type, for example, the main layer ②, connecting 3-3' to obtain the stratigraphic separation line on one side of the main layer, determining the estimated point 2 on the line segment AB, and connecting 2-2' to obtain the stratigraphic separation line on the other side of the main layer.

[0173] Further, when determining the estimated point 2, the electronic device makes the stratigraphic separation line after connecting 2-2' as parallel as possible to the stratigraphic separation line after connecting 3-3'.

[0174] If the main layer is an interdying layer, the two layer separation points on the two sides of the main layer are connected to obtain the layer separation lines on the two sides of the main layer.

[0175] Specifically, the main layer ③④ is an interdying layer. Taking the main layer ④ as an example, the electronic device connects 4-4' and 4-3' to obtain the layer separation lines on the two sides of the main layer ④.

[0176] If the main layer is a semi-interdying layer and the main layer is not the last layer, an estimation point is determined on the layer discontinuity line closest to the main layer, the layer separation point on one side of the main layer is connected to the adjacent estimation point to obtain the layer separation line on one side of the main layer, and the layer separation point on the other side of the main layer is connected to the estimation point to obtain the layer separation line on the other side of the main layer, and the two layer separation lines on the two sides of the main layer are parallel.

[0177] Specifically, the main layer ① is a semi-interdying layer and is not the last layer, and the electronic device determines an estimation point 1 on the line segment AB, wherein 2 is an adjacent estimation point, and the electronic device connects 2-2' and 1-1' to obtain the layer separation lines on the two sides of the main layer ①.

[0178] If the main layer is a semi-interdying layer and the main layer is the last layer, the layer separation points of the main layer are connected to the adjacent estimation points to obtain the layer separation lines on the two sides of the main layer.

[0179] Specifically, the layer above the main layer ① is the last semi-interdying layer, and 1-1' and point 1-B are connected to obtain the layer separation lines on the two sides of the main layer.

[0180] Step S25: determining the layer separation lines on the two sides of the interlayer based on the layer separation lines on the two sides of the main layer and the layer separation points corresponding to the interlayer. Specifically, referring to Figure 8 , the line segment AB and the line segment CD are both projective axes, and the line segment AC and the line segment CD are the layer separation lines on the upper and lower sides of the main layer. Step S25 includes at least one of the following:

[0181] If there are layer separation points on the two sides of the interlayer on both projective axes, the layer separation lines on the two sides of the interlayer are obtained by connecting the layer separation points on the two sides of the interlayer, respectively.

[0182] For example, there are layer separation points 3, 3', 4 and 4' on the two sides of the interlayer ② on both projective axes, and the electronic device connects 3-4' and 4-4' to obtain the layer separation lines on the two sides of the interlayer ②.

[0183] If there are layer separation points on one side of the interlayer on both projective axes, and there are layer separation points on the other side of the interlayer on only one projective axis, the layer separation lines on the two sides of the interlayer are obtained by connecting the two layer separation points on one side of the interlayer and determining an estimation point on the layer separation line of the main layer on the other side of the interlayer, and connecting the estimation point and the layer separation point on the other side of the interlayer.

[0184] For example, the electronic device determines the estimation point 6' on the stratigraphic separation line CD of the main layer, and connects 5-5', 6-6' to obtain the stratigraphic separation lines on both sides of the interlayer ③.

[0185] If there are stratigraphic separation points on both sides of the interlayer on one of the projection axes, and there are no stratigraphic separation points on both sides of the interlayer on the other projection axis, then the interval point is determined between the two stratigraphic separation lines of the main layer and the two projection axes, and the stratigraphic separation points are connected to the interval point to obtain the stratigraphic separation lines on both sides of the interlayer.

[0186] For example, the electronic device determines the interval point 1', and connects 1-1', 2-1' to obtain the stratigraphic separation lines on both sides of the interlayer ①.

[0187] Step S104: generating a three-dimensional stratigraphic pie model among any three adjacent boreholes based on the two-dimensional stratigraphic pie models between any two adjacent boreholes.

[0188] Step S1041: determining whether there are at least two groups of the same stratigraphic separation lines or stratigraphic discontinuity lines in the two-dimensional stratigraphic pie model among any three adjacent boreholes; if so, executing step S1042; if there is only one group of stratigraphic separation lines and stratigraphic discontinuity lines in the two-dimensional stratigraphic pie model among any three adjacent boreholes, executing steps S1043-S1044.

[0189] Step S1042: connecting the same stratigraphic separation lines end to end, and obtaining the stratigraphic separation surface based on the point method equation; connecting the same stratigraphic discontinuity lines end to end, and obtaining the stratigraphic discontinuity surface based on the point method equation.

[0190] Referring to Figure 9 Taking the stratigraphic separation line as an example, there are the same stratigraphic separation lines 1-1', 1'-1", and 1-1" among any three adjacent boreholes AB, CD, and EF.

[0191] The electronic device then connects the same stratigraphic separation lines end to end to obtain the contour of the stratigraphic separation surface, such as 1-1'-1".

[0192] The electronic device presets a plane normal vector And selecting any three points on the stratigraphic separation line, determining two non-collinear vectors according to the three points, and calculating the plane normal vector according to the two non-collinear vectors. For example, selecting the coordinates of points 1, 1', and 1" to obtain two non-collinear vectors and Establishing equation group ① according to the definition of the normal vector ② Solving the equation group and taking one of the solutions.

[0193] Then, the electronic device selects an arbitrary point M(x0, y0, z0) on the stratigraphic separation line, and calculates the equation of the stratigraphic separation surface according to the point-form equation A(x-x0)+B(y-y0)+C=0.

[0194] The electronic device draws the stratigraphic separation surface in the three-dimensional stratigraphic pie model according to the contour of the stratigraphic separation surface and the equation of the stratigraphic separation surface.

[0195] Step S1043: determining a stratigraphic discontinuity surface based on the stratigraphic discontinuity line, wherein the stratigraphic discontinuity surface includes a sedimentary discontinuity surface, a structural discontinuity surface, and an intrusive discontinuity surface, and includes at least one of the following:

[0196] A sedimentary estimation point is determined between any three adjacent boreholes, the sedimentary estimation point is located on the stratigraphic separation surface of the nearest main layer, both ends of the sedimentary discontinuity line are connected to the sedimentary estimation point, and the sedimentary discontinuity surface is calculated based on the point-form equation.

[0197] Specifically, according to the natural generation principle of the stratum, the sedimentary discontinuity surface is connected to the main layer, and then the electronic device determines a sedimentary estimation point on the stratigraphic separation surface of the main layer.

[0198] A structural discontinuity estimation point is determined at any position between any three adjacent boreholes, both ends of the structural discontinuity line are connected to the structural discontinuity estimation point, and the structural discontinuity surface is calculated based on the point-form equation.

[0199] An intrusive discontinuity estimation point is determined at any position between any three adjacent boreholes, both ends of the intrusive discontinuity line are connected to the intrusive discontinuity estimation point, and the intrusive discontinuity surface is calculated based on the point-form equation.

[0200] Specifically, after the electronic device connects both ends of the stratigraphic discontinuity line to the corresponding estimation points, the contour of the stratigraphic discontinuity surface is determined, two non-collinear vectors are determined according to the stratigraphic discontinuity line and the corresponding estimation points, the plane equation of the stratigraphic discontinuity surface is determined according to the point-form equation, and then the stratigraphic discontinuity surface in the three-dimensional stratigraphic pie model is obtained according to the contour of the stratigraphic discontinuity surface and the plane equation.

[0201] Step S1044: determining a stratigraphic separation surface based on a stratigraphic separation line, including at least one of the following:

[0202] At least one main layer estimation point is determined on the stratigraphic separation surface on both sides of the main layer, both ends of the stratigraphic separation line on both sides of the main layer are connected to the main layer estimation point on the same side, and the stratigraphic separation surface on both sides of the main layer is determined based on the point-form equation.

[0203] Determine at least one interlayer estimation point between any three adjacent boreholes, the interlayer estimation point being located between two stratigraphic separation surfaces on both sides of the corresponding main layer, both ends of the stratigraphic separation lines on both sides of the interlayer being connected with the interlayer estimation point at both ends, and the stratigraphic separation surfaces on both sides of the interlayer being determined based on the point equation calculation.

[0204] Specifically, since the main layer is adjacent to the main layer on both sides of the main layer, the stratigraphic separation surfaces on both sides of the main layer are adjacent to the main layer on both sides; the interlayer is located in the main layer, and the interlayer estimation point is located between the stratigraphic separation surfaces on both sides of the main layer.

[0205] Referring to Figure 10 , in the figure: A: tectonic discontinuity surface; B: stratigraphic pie of main layer; C: multi-layer discontinuity; D: stratigraphic pie of interlayer; E: stratigraphic pie of interlayer (single-layer discontinuity); F: sedimentary discontinuity surface; G: intrusive discontinuity surface.

[0206] Step S1045: Determine the stratigraphic pie of the interlayer based on the stratigraphic profile of the interlayer corresponding to any two adjacent boreholes and the two adjacent stratigraphic discontinuity surfaces.

[0207] Step S1046: Determine the stratigraphic pie of the main layer based on the stratigraphic profile of the main layer corresponding to any two adjacent boreholes and the stratigraphic separation surfaces on both sides of the same main layer.

[0208] Step S1047: Determine the stratigraphic pie of the interlayer based on the stratigraphic profile of the interlayer corresponding to any two adjacent boreholes and the stratigraphic separation surfaces on both sides of the same interlayer.

[0209] Specifically, the electronic device can generate the stratigraphic pie of the interlayer by splicing the top surface, the bottom surface and the plurality of side surfaces to form the stratigraphic pie, and thus the same method can be used.

[0210] When the electronic device determines the stratigraphic pie of each layer, the second occurrence information of the layer can be determined according to the stratigraphic separation surface, so as to achieve the purpose of intelligently deriving the occurrence.

[0211] Step S1048: Assign different colors to any two adjacent stratigraphic pies.

[0212] Specifically, after the electronic device assigns different colors to the two adjacent stratigraphic pies, the user can intuitively observe the three-dimensional stratigraphic pie model.

[0213] Therefore, the electronic device can intelligently derive the occurrence of all layers according to the borehole data without the occurrence information of all layers, reduce manual intervention, quickly draw a three-dimensional stratigraphic pie model, and improve work efficiency. And the operation occupies small memory, and the requirement for computer memory is low.

[0214] In another possible implementation, the method further includes: obtaining two points determined by the user above the three-dimensional stratigraphic pie model; determining a tangent line connecting the two points; and generating a profile cut based on the tangent line in the vertical direction.

[0215] Specifically, the user selects two points in the interface of the three-dimensional stratigraphic pie model according to the direction of the profile cut, and the electronic device generates the profile cut after cutting the three-dimensional stratigraphic pie model in the vertical direction based on the tangent line between the two points. This facilitates the user to select different observation angles according to actual needs.

[0216] In another possible implementation, the electronic device can further generate a borehole profile or a histogram in response to a user operation.

[0217] To better implement the method, the embodiments of the present application further provide a geological modeling device based on stratigraphic pie layering, which refers to Figure 11 The geological modeling device based on stratigraphic pie layering 200 includes:

[0218] The obtaining module 201 is configured to obtain borehole data of a plurality of boreholes, wherein the borehole data of each borehole includes borehole position information, inclination information of each hole section in each borehole, and geological information, and the geological information includes geological era information, inclination depth information, first occurrence information, and main lithology information.

[0219] The determining module 202 is configured to determine stratigraphic separation points and stratigraphic discontinuity points of a projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information, the elevation information, the geological era information, and the main lithology information.

[0220] The first generating module 203 is configured to determine stratigraphic separation lines and stratigraphic discontinuity lines between any two adjacent boreholes based on the stratigraphic separation points and the stratigraphic discontinuity points, and generate a two-dimensional stratigraphic pie model between the any two adjacent boreholes.

[0221] The second generating module 204 is configured to generate a three-dimensional stratigraphic pie model between any three adjacent boreholes based on the two-dimensional stratigraphic pie model between the any two adjacent boreholes.

[0222] Further, the determining module 202 is specifically configured to:

[0223] determine a plurality of stratigraphic separation points based on the geological era information, wherein the geological eras of the stratigraphic layers on the upper and lower sides of the stratigraphic separation points are different;

[0224] determine stratigraphic discontinuity points from the plurality of stratigraphic separation points based on the geological era information and the main lithology information, wherein the geological eras of the stratigraphic layers on the upper and lower sides of the stratigraphic discontinuity points are discontinuous;

[0225] Determine the positions of the stratigraphic division points and the stratigraphic discontinuity points on the projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information and the elevation information.

[0226] Further, the determining module 202 is specifically configured to determine the stratigraphic discontinuity points from the stratigraphic division points based on the geological era information and the main lithology information, and is specifically configured to perform at least one of the following:

[0227] If the geological eras of the strata on the upper and lower sides of the stratigraphic division point are discontinuous or discontinuous and reversed, the stratigraphic division point is determined as a sedimentary discontinuity point;

[0228] If the geological eras of the strata on the upper and lower sides of the stratigraphic division point are discontinuous and the main lithology information includes tectonic rocks, the stratigraphic division point is determined as a structural discontinuity point;

[0229] If the geological eras of the strata on the upper and lower sides of the stratigraphic division point are discontinuous and the main lithology information includes intrusive rocks, the stratigraphic division point is determined as an intrusive discontinuity point.

[0230] Further, the first generating module 203 is specifically configured to:

[0231] Connect a homodromous fold line between the top ends of the projection axes of any two adjacent boreholes and between the bottom ends of the two projection axes to generate a two-dimensional stratigraphic pie model between the two adjacent boreholes;

[0232] Determine a stratigraphic discontinuity line based on the stratigraphic discontinuity points on the two projection axes, wherein any two adjacent stratigraphic discontinuity lines and the area surrounded by the two homodromous fold lines form a stratigraphic profile;

[0233] Determine a stratigraphic division line of each stratigraphic profile based on the stratigraphic profile and the stratigraphic division point.

[0234] Further, the first generating module 203 is specifically configured to determine the stratigraphic discontinuity line based on the stratigraphic discontinuity points on the two projection axes, and is specifically configured to perform at least one of the following:

[0235] Determine whether there are the same sedimentary discontinuity points on the two projection axes, if yes, connect the two same sedimentary discontinuity points to obtain a sedimentary discontinuity line; otherwise, determine an estimated point on the nearest homodromous fold line or stratigraphic discontinuity line, connect the estimated point and the sedimentary discontinuity point to obtain the sedimentary discontinuity line, and the sedimentary discontinuity line is parallel to the adjacent strata;

[0236] Determine whether there are the same intrusive discontinuity points on the two projection axes, if yes, connect the two same intrusive discontinuity points to obtain an intrusive discontinuity line; otherwise, obtain first occurrence information at the intrusive discontinuity point, determine the slope of the intrusive discontinuity line based on the first occurrence information, and determine the intrusive discontinuity line between the projection axis and the nearest stratigraphic discontinuity line based on the intrusive discontinuity point and the slope;

[0237] determining whether there is a same structural discontinuity point on the two projection axis lines, if yes, connecting the two same structural discontinuity points to obtain a structural discontinuity line; otherwise, determining a slope of the intrusive discontinuity line based on the occurrence information at the structural discontinuity point, determining an estimated point on the same direction polygonal line based on the structural discontinuity point and the slope, connecting the estimated point and the structural discontinuity point to obtain the structural discontinuity line.

[0238] Further, the first generating module 203, when determining the stratigraphic separation line of each stratigraphic profile based on the stratigraphic profile and the stratigraphic separation point, is specifically used for:

[0239] determining whether there is a stratigraphic separation point on a line segment between the two adjacent stratigraphic discontinuity lines; if yes, determining that the stratigraphic profile is a multi-layer profile; if not, determining that the stratigraphic profile is a single-layer profile;

[0240] if the stratigraphic profile is a multi-layer profile, determining a main layer and a interlayer based on the stratigraphic separation point and the geological era information;

[0241] determining the type of each main layer based on the stratigraphic separation point corresponding to each main layer, including at least one of the following:

[0242] a connected type of stratum: there are two stratigraphic separation points of the main layer on the two projection axis lines;

[0243] a semi-connected type of stratum: there are stratigraphic separation points of the main layer on the two projection axis lines on one side, and there are stratigraphic separation points of the main layer on the other side on only one of the projection axis lines;

[0244] an inter-extinction type of stratum: there are stratigraphic separation points of the main layer on the two projection axis lines on the upper and lower sides and do not coincide, and there are stratigraphic separation points of the main layer on the two projection axis lines on the upper and lower sides and coincide;

[0245] a semi-inter-extinction type of stratum: there are stratigraphic separation points of the main layer on the two projection axis lines on the upper and lower sides, and the two stratigraphic separation points do not coincide;

[0246] determining the stratigraphic separation line on both sides of the main layer based on the type of the main layer, including:

[0247] connecting the two stratigraphic separation points on both sides of the main layer of the connected type of stratum to obtain the stratigraphic separation line on both sides of the main layer;

[0248] starting from the connected type of stratum, determining the stratigraphic separation line corresponding to the remaining main layer layer by layer upwards and downwards at the same time, including any one of the following:

[0249] if the main layer is a semi-connected type of stratum, connecting the two stratigraphic separation points on one side of the main layer to obtain the stratigraphic separation line on one side of the main layer, determining an estimated point on the stratigraphic discontinuity line on the other side of the main layer, and connecting the stratigraphic separation point on the other side of the main layer and the estimated point to obtain the stratigraphic separation line on the other side of the main layer.

[0250] If the main layer is an inter-extinction layer, the two layer separation points on the two sides of the main layer are connected respectively to obtain the layer separation lines on the two sides of the main layer;

[0251] If the main layer is a semi-inter-extinction layer and the main layer is not the last layer, an estimation point is determined on the layer discontinuity line closest to the main layer, the layer separation point on one side of the main layer is connected with the adjacent estimation point to obtain the layer separation line on one side of the main layer, and the layer separation point on the other side of the main layer is connected with the estimation point to obtain the layer separation line on the other side of the main layer, and the two layer separation lines on the two sides of the main layer are parallel;

[0252] If the main layer is a semi-inter-extinction layer and the main layer is the last layer, the layer separation points of the main layer are connected with the adjacent estimation points respectively to obtain the layer separation lines on the two sides of the main layer;

[0253] Based on the layer separation lines on the two sides of the main layer and the layer separation points corresponding to the interlayer, the layer separation lines on the two sides of the interlayer are determined.

[0254] Further, the first determination module 203, when determining the layer separation lines on the two sides of the interlayer based on the layer separation lines on the two sides of the main layer and the layer separation points corresponding to the interlayer, is specifically used for at least one of the following:

[0255] If there are layer separation points on the two projection axes on the two sides of the interlayer, the layer separation lines on the two sides of the interlayer are obtained by connecting the layer separation points on the two sides of the interlayer respectively;

[0256] If there are layer separation points on the two projection axes on one side of the interlayer, and there is only a layer separation point on the other side of the interlayer on one of the projection axes, the layer separation lines on the two sides of the interlayer are obtained by connecting the two layer separation points on one side of the interlayer and determining an estimation point on the layer separation line of the main layer on the other side of the interlayer, and connecting the estimation point and the layer separation point on the other side of the interlayer.

[0257] If there are layer separation points on the two projection axes on the two sides of the interlayer, and there are no layer separation points on the two projection axes on the other side of the interlayer, the layer separation lines on the two sides of the interlayer are obtained by determining an inter-extinction point between the two layer separation lines of the main layer and the two projection axes, and connecting the layer separation points with the inter-extinction point.

[0258] Further, the second generation module 204 is specifically used for:

[0259] Judging whether there are at least two groups of same layer separation lines or layer discontinuity lines in a two-dimensional layer pie model between any three adjacent boreholes, if yes, connecting the same layer separation lines end to end, and obtaining a layer separation surface based on a point method equation; connecting the same layer discontinuity lines end to end, and obtaining a layer discontinuity surface based on a point method equation;

[0260] Otherwise, determining the stratigraphic discontinuity surface based on the stratigraphic discontinuity line, wherein the stratigraphic discontinuity surface comprises the sedimentary discontinuity surface, the structural discontinuity surface and the intrusive discontinuity surface, and comprises at least one of the following:

[0261] Determining a sedimentary estimation point between any three adjacent boreholes, the sedimentary estimation point being located on the stratigraphic separation surface of the nearest main layer, connecting both ends of the sedimentary discontinuity line to the sedimentary estimation point, and calculating the sedimentary discontinuity surface based on the point method equation;

[0262] Determining a structural discontinuity estimation point at any position between any three adjacent boreholes, connecting both ends of the structural discontinuity line to the structural discontinuity estimation point, and calculating the structural discontinuity surface based on the point method equation;

[0263] Determining an intrusive discontinuity estimation point at any position between any three adjacent boreholes, connecting both ends of the intrusive discontinuity line to the intrusive discontinuity estimation point, and calculating the intrusive discontinuity surface based on the point method equation;

[0264] Determining the stratigraphic separation surface based on the stratigraphic separation line, and comprising at least one of the following:

[0265] Determining at least one main layer estimation point on the stratigraphic separation surface on both sides of the main layer, connecting both ends of the stratigraphic separation line on both sides of the main layer to the main layer estimation point on the same side, and calculating the stratigraphic separation surface on both sides of the main layer based on the point method equation;

[0266] Determining at least one interlayer estimation point between any three adjacent boreholes, the interlayer estimation point being located between the two stratigraphic separation surfaces on both sides of the corresponding main layer, connecting both ends of the stratigraphic separation line on both sides of the interlayer to the interlayer estimation point, and calculating the stratigraphic separation surface on both sides of the interlayer based on the point method equation;

[0267] Determining the stratigraphic pie of the discontinuity layer based on the stratigraphic profile corresponding to the discontinuity layer between two adjacent stratigraphic discontinuity surfaces and any two adjacent boreholes;

[0268] Determining the stratigraphic pie of the main layer based on the stratigraphic profile corresponding to the main layer between any two adjacent boreholes and the stratigraphic separation surface on both sides of the same main layer;

[0269] Determining the stratigraphic pie of the interlayer based on the stratigraphic profile corresponding to the interlayer between any two adjacent boreholes and the stratigraphic separation surface on both sides of the same interlayer;

[0270] Combining the stratigraphic pie of the discontinuity layer, the stratigraphic pie of the main layer and the stratigraphic pie of the interlayer to generate a three-dimensional stratigraphic pie model between any three adjacent boreholes;

[0271] Assigning different colors to any two adjacent stratigraphic pies.

[0272] In another possible implementation, the geological modeling apparatus 200 for stratigraphic discontinuity pie layering further comprises: a point acquisition module configured to acquire two points determined by a user above the three-dimensional stratigraphic pie model;

[0273] a tangent line determination module configured to determine a tangent line based on the two points;

[0274] a graph profile generation module configured to generate a graph profile by longitudinally cutting the three-dimensional stratigraphic pie model along a vertical direction based on the tangent line.

[0275] The various changes and specific examples in the method in the foregoing embodiments are also applicable to the geological modeling apparatus 200 for stratigraphic discontinuity pie layering in this embodiment. Based on the detailed description of the geological modeling method for stratigraphic discontinuity pie layering, those skilled in the art can clearly understand the implementation method of the geological modeling apparatus 200 for stratigraphic discontinuity pie layering in this embodiment. Therefore, for the sake of brevity of the specification, the geological modeling apparatus 200 for stratigraphic discontinuity pie layering in this embodiment will not be described in detail.

[0276] To better implement the above method, an electronic device is provided in an embodiment of the present application, which refers to Figure 12 The electronic device 300 comprises a communication control board 301, a memory 303 and a display screen 305. The memory 303 and the display screen 305 are connected to the communication control board 301, for example, through a bus 302. Optionally, the electronic device 300 can further comprise a transceiver 304. It should be noted that the transceiver 304 is not limited to one in actual application, and the structure of the electronic device 300 does not constitute a limitation on the embodiments of the present application.

[0277] The communication control board 301 can be a CPU (Central Processing Unit, central processing unit), a general-purpose processor, a DSP (Digital Signal Processor, digital signal processor), an ASIC (Application Specific Integrated Circuit, application specific integrated circuit), an FPGA (Field Programmable Gate Array, field programmable gate array) or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. It can implement or execute various exemplary logical blocks, modules and circuits described in combination with the disclosure of the present application. The communication control board 301 can also be a combination of computing functions, such as a combination of one or more microprocessors, a combination of DSP and microprocessor, etc.

[0278] The bus 302 can include a path that transmits information between the above-mentioned components. The bus 302 can be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, or the like. The bus 302 can be divided into an address bus, a data bus, a control bus, and the like.

[0279] The memory 303 can be a ROM (Read Only Memory) or other type of static storage device that can store static information and instructions, a RAM (Random Access Memory) or other type of dynamic storage device that can store information and instructions, an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read Only Memory) or other optical disk storage, a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and that can be accessed by a computer, but is not limited thereto.

[0280] The memory 303 is used to store application program codes for implementing the scheme of the present application, and the execution is controlled by the communication control board 301. The communication control board 301 is used to execute the application program codes stored in the memory 303 to realize the content shown in the foregoing method embodiments.

[0281] Figure 12 The electronic device 300 shown is only an example, and should not bring any limitation to the function and use range of the embodiments of the present application.

[0282] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, and the program is executed by a processor to realize the geological modeling method based on stratigraphic discontinuity pie layering provided by the above embodiment. The processor executes the computer program in the computer readable storage medium, determines stratigraphic division points and stratigraphic discontinuity points on each borehole according to borehole data of the borehole, determines stratigraphic division lines based on the stratigraphic division points between two adjacent boreholes, determines stratigraphic discontinuity lines based on the stratigraphic discontinuity points between the two adjacent boreholes, and then generates a two-dimensional stratigraphic pie model between the two adjacent boreholes, and then generates a three-dimensional stratigraphic pie model between any three boreholes according to the two-dimensional stratigraphic pie models between the two adjacent boreholes, so that the drawing efficiency is improved, and the working efficiency is improved.

[0283] In the embodiment, the computer readable storage medium can be a tangible device that maintains and stores instructions for use by an instruction execution device. The computer readable storage medium can be, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any combination thereof. Specifically, the computer readable storage medium can be a portable computer disk, a hard disk, a U disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a platform random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, an optical disk, a magnetic disk, a mechanical coding device, and any combination thereof.

[0284] The computer program in the embodiment includes program codes for executing all the above methods, and the program codes can include instructions corresponding to the execution of the method steps provided by the above embodiment. The computer program can be downloaded from the computer readable storage medium to each computing / processing device, or downloaded to an external computer or an external storage device through a network (such as the Internet, a local area network, a wide area network, and / or a wireless network). The computer program can be completely executed on a user computer, or executed as a separate software package.

[0285] The above are preferred embodiments of the present application, and do not limit the protection scope of the present application, so that: equivalent changes made according to the structure, shape, principle of the present application should be covered within the protection scope of the present application.

[0286] In addition, it should be understood that relational terms such as first and second and the like can be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus.

Claims

1. A method of stratigraphic modeling based on stratigraphic discontinuity slice layering, characterized by, The method comprises: obtaining drilling data of a plurality of drill holes, wherein the drilling data of each drill hole comprises drilling position information, inclination information, depth information, elevation information and geological information, and the geological information comprises geological era information and main lithology information; determining stratigraphic division points and stratigraphic discontinuity points of a projection axis of each drill hole based on the drilling position information, the inclination information, the depth information, the elevation information, the geological era information and the main lithology information; determining stratigraphic division lines and stratigraphic discontinuity lines between any two adjacent drill holes based on the stratigraphic division points and the stratigraphic discontinuity points, and generating a two-dimensional stratigraphic pie model between any two adjacent drill holes; generating a three-dimensional stratigraphic pie model between any three adjacent drill holes based on the two-dimensional stratigraphic pie model between any two adjacent drill holes; the determination of the stratigraphic division points and the stratigraphic discontinuity points of the projection axis of each drill hole based on the drilling position information, the inclination information, the depth information, the elevation information, the geological era information and the main lithology information comprises: determining a plurality of stratigraphic division points based on the geological era information, wherein the stratigraphic division points are on the upper and lower sides of strata with different geological eras; determining stratigraphic discontinuity points from the plurality of stratigraphic division points based on the geological era information and the main lithology information, wherein the stratigraphic discontinuity points are on the upper and lower sides of strata with discontinuous geological eras; the stratigraphic discontinuity points comprise sedimentary discontinuity points, structural discontinuity points and intrusive discontinuity points, and the determination of the stratigraphic discontinuity points from the stratigraphic division points based on the geological era information and the main lithology information comprises at least one of the following: if the geological eras of the strata on the upper and lower sides of the stratigraphic division points are interrupted or interrupted and reversed, the stratigraphic division points are determined as sedimentary discontinuity points; if the geological eras of the strata on the upper and lower sides of the stratigraphic division points are interrupted and the main lithology information includes structural rocks, the stratigraphic division points are determined as structural discontinuity points; if the geological eras of the strata on the upper and lower sides of the stratigraphic division points are interrupted and the main lithology information includes intrusive rocks, the stratigraphic division points are determined as intrusive discontinuity points; determining the positions of the stratigraphic division points and the stratigraphic discontinuity points on the projection axis of each drill hole based on the drilling position information, the inclination information, the depth information and the elevation information; the determination of the stratigraphic division lines and the stratigraphic discontinuity lines between any two adjacent drill holes based on the stratigraphic division points and the stratigraphic discontinuity points, and the generation of the two-dimensional stratigraphic pie model between any two adjacent drill holes, comprises: connecting a same-direction polyline between the top ends of the projection axes of any two adjacent drill holes and between the bottom ends of the two projection axes to generate a two-dimensional stratigraphic pie model between the two adjacent drill holes; determining stratigraphic discontinuity lines based on the stratigraphic discontinuity points on the two projection axes, wherein any two adjacent stratigraphic discontinuity lines and the area enclosed by the two same-direction polylines form a stratigraphic profile; determining stratigraphic division lines of each stratigraphic profile based on the stratigraphic profile and the stratigraphic division points.

2. The method of claim 1, wherein, The stratigraphic discontinuities include sedimentary discontinuities, structural discontinuities, and intrusive discontinuities. Determining stratigraphic discontinuities based on stratigraphic discontinuities on two projective axes includes at least one of the following: Determine whether there are identical sedimentary discontinuities on the two projective axes. If so, connect the two identical sedimentary discontinuities to obtain a sedimentary discontinuity line. Otherwise, determine an estimation point on the nearest unidirectional broken line or stratigraphic discontinuity line, and connect the estimation point and the sedimentary discontinuity to obtain a sedimentary discontinuity line, which is parallel to the adjacent strata. Determine whether there are identical intrusion discontinuities on the two projective axes. If so, connect the two identical intrusion discontinuities to obtain an intrusion discontinuity line. Otherwise, obtain the first attitude information at the intrusion discontinuity, determine the slope of the intrusion discontinuity line based on the first attitude information, and determine the intrusion discontinuity line between the projective axis and the nearest stratigraphic discontinuity line based on the intrusion discontinuity and the slope. Determine whether there are identical structural discontinuities on the two projective axes. If so, connect the two identical structural discontinuities to obtain a structural discontinuity line. Otherwise, determine the slope of the intrusive discontinuity line based on the occurrence information at the structural discontinuity point. Determine an estimation point on the same-direction broken line based on the structural discontinuity point and the slope. Connect the estimation point and the structural discontinuity point to obtain a structural discontinuity line.

3. The method of claim 1, wherein, The determination of the dividing line for each stratigraphic profile based on the stratigraphic discontinuity and the dividing point includes: Determine whether there is a stratigraphic separation point on the line segment between the two adjacent stratigraphic discontinuities of the two projective axes; if yes, then the stratigraphic profile is determined to be a multi-layer profile; if no, then the stratigraphic profile is determined to be a single-layer profile. If the stratigraphic profile is a multi-layer profile, the main layer and interlayers are determined based on the geological age information and stratigraphic separation points; Based on the stratigraphic boundary points corresponding to each principal layer, the type of each principal layer is determined, including at least one of the following: Connected strata: There are two stratigraphic separation points of the main layer on both of the two projective axes; Semi-connected strata: There are stratigraphic separation points on one side of the main layer on both of the two projective axes, and there is a stratigraphic separation point on the other side of the main layer on only one of the projective axes; Intermittent strata: One of the projective axes has stratigraphic separation points on the upper and lower sides of the main layer that do not coincide, and the other projective axis has stratigraphic separation points on the upper and lower sides of the main layer that coincide; Semi-extinction type strata: There are stratigraphic separation points on both sides of the main layer only on any one of the projective axes, and the two stratigraphic separation points do not coincide; Based on the type of the main layer, the stratigraphic boundary lines on both sides of the main layer are determined, including: Connect the two stratigraphic separation points on both sides of the main layer of a connected stratum to obtain the stratigraphic separation lines on both sides of the main layer; Starting with connected strata, simultaneously determine the stratigraphic boundaries of the remaining main strata layer by layer, both upwards and downwards, including any of the following: If the main layer is a semi-continuous layer, two layer separation points on one side of the main layer are connected to obtain a layer separation line on one side of the main layer, an estimated point is determined on the layer interval line on the other side of the main layer, and the estimated point and the layer separation point on the other side of the main layer are connected to obtain a layer separation line on the other side of the main layer; If the main layer is an interlayer, two layer separation points on both sides of the main layer are connected to obtain layer separation lines on both sides of the main layer; If the main layer is a semi-interlayer and the main layer is not the last layer, an estimated point is determined on the layer interval line closest to the main layer, the layer separation point on one side of the main layer is connected to the adjacent estimated point to obtain a layer separation line on one side of the main layer, and the layer separation point on the other side of the main layer is connected to the estimated point to obtain a layer separation line on the other side of the main layer, and the two layer separation lines on both sides of the main layer are parallel; If the main layer is a semi-interlayer and the main layer is the last layer, the layer separation points of the main layer are connected to the adjacent estimated points to obtain layer separation lines on both sides of the main layer; Based on the layer separation lines on both sides of the main layer and the corresponding layer separation points of the interlayer, layer separation lines on both sides of the interlayer are determined.

4. The method of claim 3, wherein, The determination of the layer separation lines on both sides of the interlayer based on the layer separation lines on both sides of the main layer and the corresponding layer separation points of the interlayer includes at least one of the following: If there are layer separation points on both sides of the interlayer on both projection axes, layer separation lines on both sides of the interlayer are obtained by connecting the layer separation points on both sides of the interlayer, respectively; If there are layer separation points on one side of the interlayer on both projection axes, and there is only one layer separation point on the other side of the interlayer on one of the projection axes, the two layer separation points on one side of the interlayer are connected, an estimated point is determined on the layer separation line of the main layer on the other side of the interlayer, and the estimated point and the layer separation point on the other side of the interlayer are connected to obtain layer separation lines on both sides of the interlayer; If there are layer separation points on both sides of the interlayer on one of the projection axes, and there are no layer separation points on both sides of the interlayer on the other projection axis, an interlayer point is determined between the two layer separation lines of the main layer and the two projection axes, and the layer separation points on both sides of the interlayer are connected to the interlayer point to obtain layer separation lines on both sides of the interlayer.

5. The method according to any one of claims 1-4, characterized in that, The generation of a three-dimensional layer pie model among any three adjacent boreholes based on a two-dimensional layer pie model between any two adjacent boreholes includes: determining whether there are at least two groups of identical layer separation lines or layer interval lines in the two-dimensional layer pie model among any three adjacent boreholes, if so, connecting the identical layer separation lines end to end and calculating a layer separation surface based on a point method equation; connecting the identical layer interval lines end to end and calculating a layer interval surface based on a point method equation; otherwise, determining a layer interval surface based on the layer interval line, wherein the layer interval surface includes a deposition interval surface, a structure interval surface and an intrusion interval surface, and includes at least one of the following: determining a deposition estimation point between the three adjacent boreholes, the deposition estimation point being located on a stratigraphic separation surface of a nearest main layer, connecting two ends of the deposition discontinuity line to the deposition estimation point, and calculating a deposition discontinuity surface based on a point equation; determining a structure discontinuity estimation point at any position between the three adjacent boreholes, connecting two ends of the structure discontinuity line to the structure discontinuity estimation point, and calculating a structure discontinuity surface based on a point equation; determining an intrusion discontinuity estimation point at any position between the three adjacent boreholes, connecting two ends of the intrusion discontinuity line to the intrusion discontinuity estimation point, and calculating an intrusion discontinuity surface based on a point equation; determining a stratigraphic separation surface based on the stratigraphic separation line, including at least one of: determining at least one main layer estimation point on the stratigraphic separation surface on both sides of the main layer, connecting two ends of the stratigraphic separation line on both sides of the main layer to the main layer estimation points on the same side, and calculating stratigraphic separation surfaces on both sides of the main layer based on a point equation; determining at least one interlayer estimation point between the three adjacent boreholes, the interlayer estimation point being located between two stratigraphic separation surfaces on both sides of a corresponding main layer, connecting two ends of the stratigraphic separation line on both sides of the interlayer to the interlayer estimation point, and calculating stratigraphic separation surfaces on both sides of the interlayer based on a point equation; determining a stratigraphic pie of an interlayer based on stratigraphic profiles of the interlayer corresponding to two adjacent stratigraphic discontinuity surfaces and any two adjacent boreholes; determining a stratigraphic pie of a main layer based on stratigraphic profiles of the main layer corresponding to stratigraphic separation surfaces on both sides of the main layer and any two adjacent boreholes; determining a stratigraphic pie of an interlayer based on stratigraphic profiles of the interlayer corresponding to stratigraphic separation surfaces on both sides of the interlayer and any two adjacent boreholes; combining the stratigraphic pies of the interlayer, the main layer and the interlayer to generate a three-dimensional stratigraphic pie model between any three adjacent boreholes; assigning different colors to any two adjacent stratigraphic pies.

6. An apparatus for geologically modeling stratigraphic breaks slice-by-slice based on slice-by-slice pie bodies, characterized by, including: an acquisition module configured to acquire borehole data of a plurality of boreholes, wherein the borehole data of each borehole includes borehole position information, inclination information, inclination depth information, elevation information and geological information, and the geological information includes geological era information and main lithology information; a determination module configured to determine stratigraphic separation points and stratigraphic discontinuity points of a projection axis of each borehole based on the borehole position information, the inclination information, the inclination depth information, the elevation information, the geological era information and the main lithology information; a first generation module configured to determine stratigraphic separation lines and stratigraphic discontinuity lines between any two adjacent boreholes based on the stratigraphic separation points and the stratigraphic discontinuity points, and generate a two-dimensional stratigraphic pie model between any two adjacent boreholes; a second generation module configured to generate a three-dimensional stratigraphic pie model between any three adjacent boreholes based on the two-dimensional stratigraphic pie model between any two adjacent boreholes; the determination module is specifically configured to: determine a plurality of stratigraphic separation points based on the geological era information, wherein the geological eras of the stratigraphic separation points on both sides of the stratigraphic separation points are different. determine, based on the geologic age information and the main lithology information, a stratigraphic break point from the plurality of stratigraphic division points, wherein the stratigraphic break point is discontinuous in geologic age of the stratigraphic layers on both sides of the stratigraphic break point; The stratigraphic break point includes a sedimentary break point, a structural break point, and an intrusive break point. Based on the geologic age information and the main lithology information, the stratigraphic break point is determined from the stratigraphic division points, including at least one of the following: If the geologic age of the stratigraphic layers on both sides of the stratigraphic division point is discontinuous or discontinuous and reversed, the stratigraphic division point is determined as a sedimentary break point. If the geologic age of the stratigraphic layers on both sides of the stratigraphic division point is discontinuous and the main lithology information includes tectonic rock, the stratigraphic division point is determined as a structural break point. If the geologic age of the stratigraphic layers on both sides of the stratigraphic division point is discontinuous and the main lithology information includes intrusive rock, the stratigraphic division point is determined as an intrusive break point. Based on the borehole position information, the inclination information, the inclination depth information, and the elevation information, the positions of the stratigraphic division points and the stratigraphic break points on the projection axis of each borehole are determined. The first generation module is specifically configured to: Connect a same-direction fold line between the top ends of the projection axes of any two adjacent boreholes and between the bottom ends of the two projection axes to generate a two-dimensional stratigraphic pie model between the two adjacent boreholes. Determine a stratigraphic break line based on the stratigraphic break points on the two projection axes, wherein the area enclosed by any two adjacent stratigraphic break lines and two same-direction fold lines is a stratigraphic profile. Determine a stratigraphic division line of each stratigraphic profile based on the stratigraphic profile and the stratigraphic division point.

7. An electronic device, comprising: Comprise: At least one processor; Memory; At least one application program, wherein the at least one application program is stored in the memory and configured to be executed by the at least one processor, and the at least one application program is configured to execute the method of any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that: A computer program is stored, which can be loaded and executed by a processor to execute the method of any one of claims 1 to 5.

Images