CN117876619A

CN117876619A - BIM-based complex three-dimensional geological model automatic modeling method

Info

Publication number: CN117876619A
Application number: CN202410050513.7A
Authority: CN
Inventors: 张志峰; 李翻翻; 吴志刚; 张胜; 赵华宏; 杨大海; 陈修和; 郭城; 陈星�; 刘亦舟
Original assignee: Anhui Transport Consulting and Design Institute Co Ltd
Current assignee: Anhui Transport Consulting and Design Institute Co Ltd
Priority date: 2024-01-12
Filing date: 2024-01-12
Publication date: 2024-04-12

Abstract

The invention discloses a BIM-based complex three-dimensional geological model automatic modeling method, which relates to the technical field of engineering information models and comprises the following steps of: acquiring space coordinates of an orifice point, a stratum demarcation point and a stratum bottom point of a drilled hole and stratum codes among the points to form a coding point; comparing the coding points of all the drill holes from top to bottom to obtain pinch-out and lenticular stratum codes; JT drilling holes with pinch-out and lenticular stratum codes are obtained, and the pinch-out stratum codes and the lenticular stratum codes are distinguished; establishing topography, a conventional stratum interface and a stratum interface to be modified according to conventional stratum codes; establishing a lens stratum body according to the lens stratum body code; establishing an pinch-out stratum body according to the pinch-out stratum code; and correcting the stratum interface to be corrected according to the pinch-out and the lenticular stratum and establishing a conventional stratum. The invention determines the pinch-out and lenticular stratum by combining stratum number comparison, and then generates a geological model by spatial interpolation and curve fitting.

Description

BIM-based complex three-dimensional geological model automatic modeling method

Technical Field

The invention relates to the technical field of engineering information models, in particular to a BIM-based complex three-dimensional geological model automatic modeling method.

Background

With the rise of digital wave of geotechnical engineering, how to build a three-dimensional geological model meeting engineering requirements in a BIM (building information model) platform becomes a problem to be solved urgently. Most of existing BIM software supports basic elements (space points, lines, planes and volumes) of three-dimensional geologic model construction, but when the BIM software is adopted to build a complex stratum model containing pinch-out and lens bodies, technicians usually judge the pinch-out and lens body stratum first, then extract space coordinates of corresponding stratum in a drilling hole to manually outline section curves of the pinch-out and lens body stratum, and then fit curved surfaces by the section curves to generate volumes. The method has the following defects: (1) The difficulty of manually judging the pinch-out and the stratum of the lens body in space is high; (2) The construction of fine pinch-out and lenticular stratum requires a great deal of section curves to be drawn manually, and the workload is high; (3) The formation model is built based on the pinch-out and lenticular formation section curves by hand, which is too subjective.

In order to solve the above-mentioned drawback (1), some scholars judge that the stratum missing in some drilling holes is a pinch-out or lens stratum by establishing a standard layer or two-by-two comparison mode, but in actual engineering, the stratum which is not disclosed in the shallow drilling holes is judged as a pinch-out or lens stratum due to the interference of the shallow drilling holes (the drilling depth is shallower than other drilling holes and all stratum is not disclosed); in order to solve the defects (2) and (3), part of students encrypt through algorithms such as kriging, inverse distance and discrete smooth interpolation by setting a rule for calculating the boundary of a pinch-out or lenticular stratum, for example, chinese patent No. 11 and 26 published in 2021, namely, a method and a device for automatically constructing a three-dimensional model of a geologic body, calculates the spatial position of a stratum tip vanishing point through the virtual thickness obtained by interpolation in a borehole without the pinch-out stratum and the average thickness of the pinch-out stratum, and encrypts through an interpolation algorithm; the Chinese patent (DSI-based three-dimensional stratum surface pinch-out treatment method) published in 12.31.2021 firstly calculates the spatial position of a stratum tip pinch-out point according to the relation between the total thickness of a stratum containing pinch-out and a stratum without pinch-out and the thickness of the pinch-out stratum, adjacent stratum tip pinch-out points, and encrypts the stratum through algorithms such as discrete smooth interpolation; the above approach overcomes disadvantages (2) and (3) to some extent, but ignores the changing characteristics of the pinch-out and lenticular formation interface itself.

Disclosure of Invention

In order to overcome the defects in the prior art, the invention provides a BIM-based complex three-dimensional geologic model automatic modeling method, which comprises the steps of determining pinch-out and lenticular stratum by combining stratum number comparison, and generating a geologic model by spatial interpolation and curve fitting.

In order to achieve the above purpose, the present invention adopts the following technical scheme, including:

a BIM-based complex three-dimensional geological model automatic modeling method comprises the following steps:

step 1: acquiring space coordinates of orifice points, space coordinates of stratum demarcation points and space coordinates of hole bottom points in a drill hole, acquiring stratum codes among the points, and forming corresponding coding points according to the space coordinates of the points and the stratum codes among the points;

the encoding points comprise the spatial coordinates (X, Y, Z) of the points and the upper side encoding Noup and the lower side encoding Nodown in the format of (X, Y, Z, noup, nodown); the upper and lower side codes are the upper and lower side stratum codes of the points, wherein the upper side code of the orifice point corresponding to the code point is the top code, and the lower side code of the hole bottom point corresponding to the code point is the bottom code;

step 2: comparing the coding points of all the drill holes from top to bottom to obtain pinch-out and lenticular stratum codes;

In step 2, extracting upper and lower side codes of all code points to form a code pair-merging code pair group, wherein the format of each code pair in the code pair-merging code pair group is (Noup, nodown); combining the coding pairs in the drill holes, combining the upper coding of the ith coding pair and the lower coding of the (i+j) th coding pair to obtain the combined ith coding pair, combining to make the lower coding of the (i) th coding pair from top to bottom in each drill hole identical, and marking stratum codes which are not included in the final coding pair-combined coding pair group as pinch-out and lenticular stratum codes and marking the rest stratum codes as conventional stratum codes;

step 3: according to the pinch-out and lenticular stratum codes, drilling holes with the pinch-out and lenticular stratum codes are obtained and marked as JT drilling holes, and the JT drilling holes are compared with other drilling holes to distinguish the pinch-out stratum codes and the lenticular stratum codes;

step 4: according to the conventional stratum coding, picking up coding points, and establishing a terrain, a conventional stratum interface and a stratum interface to be modified;

in step 4, the specific process of establishing the terrain, the conventional stratum interface and the stratum interface to be modified is as follows;

step 4.1: selecting a conventional stratum code, traversing all drill holes to obtain coding points above and below the conventional stratum code;

Step 4.2: and (3) establishing topography, a conventional stratum interface and a stratum interface to be modified according to the position of the coding point obtained in the step (4.1), wherein the method comprises the following steps: if the upper coding point position coincides with the orifice point, a topographic surface is directly established by a topographic map data spatial interpolation method; if the space positions of the upper and lower coding points coincide with the stratum demarcation points, respectively establishing an upper stratum interface and a lower stratum interface to be trimmed of the stratum by a space interpolation method on the basis of the space coordinates of the upper and lower coding points; if the lower coding point position coincides with the hole bottom point, no surface is generated;

step 4.3: acquiring the stratum interfaces to be modified of all the stratum generated in the step 4.2, and carrying out modification operation on the stratum interfaces to be modified of the adjacent conventional stratum, wherein the method comprises the following steps: if the geometry of the upper and lower stratum interfaces to be modified of the selected conventional stratum is the same as that of the lower and upper stratum interfaces to be modified of the upper and lower adjacent conventional stratum, merging the upper and lower stratum interfaces to be modified of the selected conventional stratum with the lower and upper stratum interfaces to be modified of the upper and lower adjacent conventional stratum respectively to serve as stratum interfaces; if the geometrical shapes of the upper and lower stratum interfaces to be modified of the selected conventional stratum are different from those of the lower and upper stratum interfaces to be modified of the upper and lower adjacent conventional stratum, the upper stratum interface to be modified of the selected conventional stratum and the lower stratum interface to be modified of the upper adjacent conventional stratum are taken as one group, and the lower stratum interface to be modified of the conventional stratum and the upper stratum interface to be modified of the lower adjacent conventional stratum are taken as the other group to be stored respectively;

Step 5: picking up corresponding stratum demarcation points according to the stratum codes of the lens bodies and establishing stratum bodies of the lens bodies;

step 6: picking up corresponding stratum demarcation points according to the pinch-out stratum codes, and establishing pinch-out stratum bodies;

step 7: and correcting the stratum interface to be corrected according to the pinch-out and the lenticular stratum and establishing a conventional stratum.

Preferably, in step 2, the steps of obtaining the pinch-out and lenticular stratum codes are as follows:

step 2.1: the coding points in the same drilling hole are ordered from high to low according to the elevation Z; extracting upper and lower side codes of all code points to form code pair-merging code pair groups, wherein each code pair in the code pair-merging code pair groups has the format of (Noup, nodown); and setting i=1;

step 2.2: acquiring an ith code pair from top to bottom in each drilling hole, grouping the code pairs according to different lower side codes of the acquired code pairs, and ensuring that each code pair in the same group has the same lower side code; screening out the lower side code in the group into a bottom code group, and eliminating corresponding original code pairs in the code pair-combined code pair group, and marking the rest groups as code pair groups; if the number of the code pairs is more than 1, executing the step 2.3; if the number of the codes is equal to 1, executing the step 2.5; if the number of the codes is equal to 0, executing the step 2.6;

Step 2.3: selecting one of the code pairs to group, carrying out merging code pair operation, repeating the steps until the lower side code of one non-merging code pair group is the same as the lower side code of the merging code pair group, replacing the original code pair by the merged code pair in the code pair-merging code pair group, and simultaneously marking the non-merging code pair group as a main round reference code pair group;

step 2.4: taking the reference code pair group of the round as a reference, carrying out combination code pair operation on the code pair groups except the code pair group selected in the last time in the step 2.3, and replacing the original code pair with the combined code pair in the code pair-combination code pair group;

step 2.5: i, self-increasing 1, entering the next round, and returning to the step 2.2;

step 2.6: the formation codes not included in the final coding pair-merging coding pair group are noted as pinch-out and lenticular formation codes, and the remaining formation codes are noted as regular formation codes;

in steps 2.3 and 2.4, the step of merging the encoding pair operation is as follows:

step 2a: setting j=1, acquiring all code pairs of the drill holes corresponding to the code pair groups needing to be subjected to the combined code pair operation from the code pair-combined code pair group, and establishing a selected code pair expansion group;

Step 2b: extracting the ith code pair and the (i+j) th code pair of the same drill hole in the selected code pair extension group, and combining the ith code pair and the (i+j) th code pair to generate a temporary combined code pair group; the merging rule merges the upper side code of the ith code pair and the lower side code of the (i+j) th code pair;

step 2c: if the temporary combination coding pair group has the same temporary combination coding pair of the upper side coding and the lower side coding, the original coding pair is removed from the selected coding pair extension group and the coding pair-combination coding pair group in the step 2.1, the step 2b is returned, and j is set to be 0;

step 2d: removing the temporary combination coding pairs with the same lower coding and bottom coding, removing all the coding pairs of the drill holes where the temporary combination coding pairs are located from the selected coding pair expansion group, and removing the original coding pairs from the coding pair-combination coding pair group in the step 2.1;

step 2e: grouping the temporary combination coding pairs in the temporary combination coding pair grouping again, wherein the grouping rule is to classify the temporary combination coding pairs with the same lower side coding into the same grouping, and marking the temporary combination coding pairs as temporary combination coding pair sub-groupings; if the number of the sub-packets is greater than or equal to 1, executing the step 2f, otherwise executing the step 2h;

step 2f: comparing the temporary combination coding sub-packet with the appointed coding sub-packet, rejecting all coding pairs in the drill hole where the temporary combination coding sub-packet which is the same as the appointed coding sub-packet is located from the selected coding pair extension packet, rejecting the temporary combination coding sub-packet from the temporary combination coding sub-packet, and simultaneously replacing the original coding pair in the coding pair-combination coding pair group in step 2.1, and returning a matching success signal; the specified code pair group is an uncombined code pair group when the step 2.3 is executed, and is a reference code pair group of the round when the step 2.4 is executed;

Step 2g: if the temporary merging code sub-packet still exists, j is increased by 1, and the step 2b is returned; otherwise, executing the step 2h;

step 2h: if all the code pairs in the selected code pair group are removed and the code pairs are not successfully matched with the appointed code pair group, returning a matching unsuccessful signal; if the matching is successful, a success signal is returned.

Preferably, in step 3, the specific process of distinguishing the pinch-out stratum code from the lenticular stratum code is as follows;

step 3.1: acquiring point plane coordinates (X, Y) of all the drilling holes to generate plane coordinate points, and drawing a closed outer contour line according to the plane coordinate points of all the drilling holes, namely an integral outer contour line;

step 3.2: acquiring the pinch-out and lenticular stratum codes in the step 2, namely JT codes, and selecting one JT code; traversing all the drilling holes, obtaining point plane coordinates (X, Y) of corresponding JT drilling holes to generate plane coordinate points, and drawing a closed outer contour line according to the plane coordinate points of the JT drilling holes, namely a JT outer contour line;

step 3.3: according to the position relation between the JT outer contour line and the integral outer contour line, the JT code corresponding to the JT outer contour line under the following two position relations is determined as the lenticular stratum code:

positional relationship one: the JT outer contour lines are all located within the overall outer contour line;

And the position relation is II: JT outer contour line is partially overlapped with the integral outer contour line, and when the overlapping position is shifted to the non-overlapping position, the tendency of increasing and then decreasing is presented;

step 3.4: repeating steps 3.2 and 3.3 for all JT encodings; determining JT codes not determined as lenticular layer codes as pinch-out layer codes;

preferably, the closed outer contour lines in the steps 3.1 and 3.2 meet that all plane points participating in constructing the closed outer contour line are located in or on the closed outer contour line, and the area is the largest;

the method for determining the trend of increasing and then decreasing of the position relationship II in the step 3.3 is as follows:

step 3a: determining the overlapping part of the JT outer contour line and the integral outer contour line, and drawing the overlapping contour line according to whether the overlapping position is continuous on the integral outer contour line or not, specifically: if continuous, drawing a coincident contour line; if the line is discontinuous, the line is disconnected at the discontinuous position, and a plurality of coincident contour lines are drawn;

step 3b: selecting a coincident contour line, extracting two endpoints of the selected coincident contour line to draw a straight line, and obtaining the length of the straight line cut by the JT outer contour line, namely the cut length;

Step 3c: shifting the straight line to the inner side of the JT outer contour line for a plurality of times, and acquiring a cut-off length once for each shift until the straight line is not intersected with the JT outer contour line;

step 3d: acquiring all offset distances and cut-off lengths, wherein the cut-off lengths show a trend of increasing and then decreasing along with the increase of the offset distances, and returning a success signal;

step 3e: repeating the steps 3b-3d until all the coincident contour lines return a success signal in the step 3d, and determining that the position relationship II is satisfied.

Preferably, in step 4.1, the upper and lower coding points of the conventional formation coding are obtained in the following specific manner:

if two coding points of the conventional stratum coding are contained in one drilling hole, sorting according to the height values, determining an upper coding point with the largest height value as the stratum, and determining a lower coding point with the smallest height value as the stratum;

if the number of the coding points of the conventional stratum coding in one borehole is more than two, the coding points are marked as a coding point group A, the upper coding points with the largest elevation values are determined according to the magnitude sequence of the elevation values, the lower coding points with the smallest elevation values are determined as the upper coding points of the stratum, and the lower coding points are corrected according to the situation, wherein the method comprises the following concrete steps: acquiring a lenticular stratum code in the code point group A, searching code points containing the lenticular stratum code in all drilling holes, and marking the code points as the code point group B, wherein if only the lenticular stratum and the conventional stratum code are in the code point group B, the upper and lower code points of the conventional stratum are not corrected; if the coding point group B contains the adjacent regular stratum coding on the regular stratum, modifying the upper coding point of the regular stratum into the coding point with the lowest coding elevation value of the lens stratum in the drill hole; if the coding point group B contains the coding of the adjacent regular stratum under the regular stratum, the lower coding point of the regular stratum is modified to be the coding point with the highest coding elevation value of the lenticular stratum in the drilling hole.

Preferably, in step 4.3, the adjacent conventional stratum is determined by the coding pair-merging coding pair group obtained in step 2, and the determination method is as follows: acquiring the selected regular stratum code, traversing each code pair of the final code pair-merging code pair group, and taking a stratum corresponding to Noup in the code pair as an upper adjacent regular stratum if the selected regular stratum code appears as Nodown in the code pair; if the regular stratum code is selected to be Noup in the code pair, the stratum corresponding to Nodown in the code pair is taken as the next adjacent regular stratum; the final encoding pair-merging encoding pair group is the encoding pair-merging encoding pair group which completes all merging operations through the step 2.

Preferably, in step 5, the specific process of establishing the lens body stratum is as follows:

step 5.1: selecting a lens stratum code, traversing all drill holes to obtain upper stratum demarcation points and lower stratum demarcation points of the lens stratum; taking the JT outline in the step 3.2 as a boundary, respectively carrying out spatial interpolation on upper and lower stratum demarcation points of the lens stratum, and establishing upper and lower stratum interfaces of the lens stratum with limited range, which are called upper and lower limited lens stratum interfaces;

Step 5.2: taking the centroid of a closed polygon formed by JT outer contour lines as a starting point, making a plurality of radial lines outwards, and simultaneously stretching the radial lines vertically up and down to form a plurality of vertical tangential planes; adopting the vertical section to cut the upper and lower limited lens stratum interfaces, and obtaining upper and lower limited lens stratum interface curves in each vertical section; one end, coinciding with the JT outer contour line, of the vertical projection of the stratum interface curve of the lens body is recorded as a transition end point;

step 5.3: generating complete upper and lower lens stratum interface curves according to the positions of the lens stratum, wherein the curves are as follows:

if the whole lens stratum is positioned in a certain conventional stratum, solving tangent lines of the upper and lower limited lens stratum interface curves at the transition end points of the upper and lower limited lens stratum interface curves in each vertical tangent plane, and acquiring intersection points of the two tangent lines; constructing a triangle enclosing area by using the tangent intersection point and the two transition end points, searching a point with zero thickness of the lens body in the enclosing area, and establishing complete upper and lower lens body stratum interface curves in the vertical tangent plane;

if the lenticular stratum is positioned at the adjacent position of the two conventional strata, generating an intersection line of the vertical tangent plane and the interfaces of the upper stratum and the lower stratum to be trimmed of the two conventional strata at the position in each vertical tangent plane, generating a straight line through two transition end points, searching a point with zero lenticular thickness in an intersection line enclosing area of the straight line and the interfaces of the upper stratum and the lower stratum to be trimmed, and establishing complete upper lenticular stratum interface curves and lower lenticular stratum interface curves in the vertical tangent planes;

Step 5.4: extracting discrete point space coordinates on the complete upper and lower lens stratum interface curves in the step 5.3 according to the set grid size, and complementing the upper and lower stratum interfaces of the lens stratum through a space interpolation algorithm; forming a lens stratum body by complementing the upper stratum interface and the lower stratum interface of the lens stratum;

step 5.5: repeating steps 5.1-5.4 to generate all lens stratum;

in step 5.3, the point with zero thickness of the lens body is searched, and the minimum deviation between the two complete upper and lower lens body stratum interface curves fitted by the point and the discrete points on the upper and lower limited lens body stratum interface curves and the upper and lower limited lens body stratum interface curves is satisfied.

Preferably, in step 6, the concrete process of creating the pinch-out formation body is as follows:

step 6.1: selecting a to-be-built pinch-out stratum code, traversing all drill holes to obtain upper and lower stratum demarcation points of the pinch-out stratum; taking the JT outline in the step 3.2 as a boundary, respectively carrying out spatial interpolation on upper and lower stratum demarcation points of the pinch-out stratum to be built, and establishing upper and lower stratum interfaces of the pinch-out stratum with limited range, which are called upper and lower limited pinch-out stratum interfaces;

step 6.2: extracting discrete points according to the mesh set size at the position where the JT outer contour line and the integral outer contour line are not overlapped in the step 3, and taking the discrete points as vertical points to make a horizontal straight line perpendicular to the JT outer contour line; stretching the straight line up and down along the vertical direction to form a vertical tangent plane, and sectioning upper and lower limited pinch-out stratum interfaces by adopting the vertical tangent plane, wherein upper and lower limited pinch-out stratum interface curves can be obtained in each vertical tangent plane; one end of the vertical projection of the limited pinch-out stratum interface curve, which coincides with the JT outer contour line, is recorded as a transition pinch-out point;

Step 6.3: generating intersection lines of upper and lower stratum interfaces of upper and lower adjacent conventional strata of the to-be-built pinch-out stratum at the position in each vertical tangent plane, and generating complete upper and lower pinch-out stratum interface curves according to whether the upper and lower pinch-out stratum of the to-be-built pinch-out stratum are adjacent to built pinch-out stratum bodies or not, wherein the curves comprise the following concrete steps:

if the built pinch-out stratum exists in the upper and lower adjacent strata of the pinch-out stratum to be built, an upper pinch-out stratum interface curve and a lower pinch-out stratum interface curve of the built pinch-out stratum are generated in a vertical tangent plane, if the built pinch-out stratum interface curve exceeds the transition pinch-out point range, one or two of the upper and lower limited stratum interface curves of the pinch-out stratum to be built are replaced by the built pinch-out stratum interface curve on the premise of keeping the transition pinch-out point unchanged, and the pinch-out point is searched in the intersection area of a straight line formed by connecting the transition pinch-out points and the upper and lower stratum interface to be trimmed, and meanwhile, the complete upper and lower pinch-out stratum interface curves are built; the method for replacing the upper and lower limited stratum interface curves of the stratum to be pinch-out by the established pinch-out stratum interface curves comprises the following steps: if the pinch-out stratum to be built is adjacent to the pinch-out stratum to be built and meets the substitution requirement, substituting an upper limited pinch-out stratum interface curve of the pinch-out stratum to be built by a lower pinch-out stratum interface curve of the pinch-out stratum to be built which is adjacent to the pinch-out stratum to be built; if the pinch-out stratum to be built is adjacent to the built pinch-out stratum, and the substitution requirement is met, substituting the lower limited pinch-out stratum interface curve of the pinch-out stratum to be built with the upper pinch-out stratum interface curve of the built pinch-out stratum which is adjacent to the lower pinch-out stratum;

If the stratum immediately adjacent to the upper stratum and the lower stratum is not the pinch-out stratum or the pinch-out stratum is not established, generating a straight line through two transition pinch-out points, searching the pinch-out points in the intersection line enclosing area of the straight line and the upper stratum interface to be modified and the lower stratum interface to be modified, and establishing a complete upper pinch-out stratum interface curve and a complete lower pinch-out stratum interface curve;

step 6.4: extracting discrete point space coordinates on the complete upper and lower pinch-out stratum interface curves in the step 6.3 according to the set grid size, and complementing the upper and lower stratum interfaces of the pinch-out stratum through a space interpolation algorithm; generating a pinch-out stratum body by complementing the upper stratum interface and the lower stratum interface of the pinch-out stratum;

step 6.5: repeating the steps 6.1-6.4 to generate all the pinch-out stratum bodies;

in step 6.3, the searching point vanishing point meets the requirement that the deviation between two complete upper and lower pinch-out stratum interface curves fitted by the point and discrete points on the upper and lower confined pinch-out stratum interface curves and the upper and lower confined pinch-out stratum interface curves is minimum.

Preferably, in step 7, the step of correcting the stratum interface to be corrected and establishing a stratum body comprises the following specific steps:

step 7.1: acquiring a group of formation interfaces to be modified in the step 4.3, generating all coding points of the formation interfaces to be modified, and judging a modification mode of the group of formation interfaces to be modified according to the coding points, wherein the method comprises the following specific steps:

If lenticular layer coding occurs in the coding points, executing a first trimming mode;

if the pinch-out stratum coding appears in the coding points, executing a second trimming mode;

if the lens body and the pinch-out stratum code appear in the code points at the same time, executing a finishing mode III;

the first trimming mode is as follows: removing coding points containing lenticular stratum codes, only reserving one of the coding points with coincident space coordinates, and generating a conventional stratum interface by adopting a space interpolation algorithm;

the second trimming mode is as follows: according to the stratum pinch-out body established in the pinch-out stratum code searching step 6, extracting a pinch-out boundary of the stratum pinch-out body; extracting discrete points on the pinch-out boundary according to the set size of the grid to serve as pinch-out additionally-arranged stratum demarcation points; merging the code points containing the pinch-out formation code as dummy formation demarcation points within each borehole; reserving one of the coding points with coincident space coordinates, combining the additional coding points with the dummy stratum demarcation point, and generating a conventional stratum interface through a space interpolation algorithm;

the third trimming mode is as follows: removing coding points containing lens body codes by adopting a first trimming mode, and establishing a conventional stratum interface according to a second trimming mode;

Step 7.2: according to the conventional stratum interface, a conventional stratum body to be modified is established, wherein the establishment mode is formed by surrounding upper and lower adjacent conventional stratum interfaces or terrain surfaces;

step 7.3: and checking whether all the regular stratum bodies to be modified are intersected with the pinch-out and lens stratum bodies, and subtracting the pinch-out and lens stratum bodies from the regular stratum bodies through Boolean operation if the regular stratum bodies are intersected, so that the final regular stratum bodies are obtained.

Preferably, in the second modification mode of step 7.1, the step of extracting the pinch-out boundary of the stratum pinch-out body is performed according to the form of the stratum pinch-out body, and specifically includes the following steps: if the pinch-out stratum body does not exist in the stratum pinch-out bodies which are adjacent to each other from top to bottom, the position where the thickness of the pinch-out stratum body is zero is extracted as a pinch-out boundary; if the pinch-out stratum body exists in the upper and lower adjacent pinch-out stratum bodies, extracting the position where the sum of the thicknesses of all the adjacent pinch-out stratum bodies is zero as a pinch-out boundary; in the second trimming mode 7.1, the virtual stratum demarcation points are combined according to the form of stratum pinch-out bodies, and the method specifically comprises the following steps: if the pinch-out stratum body does not exist in the stratum pinch-out bodies which are adjacent to each other from top to bottom, solving a space coordinate average value of a coding point containing the pinch-out stratum code in each drilling hole as a dummy stratum demarcation point; if the pinch-out formation exists in the upper and lower immediate adjacent pinch-out formations, solving a space coordinate average value of coding points comprising the pinch-out formation and the upper and lower immediate adjacent formation codes in each borehole as a dummy formation demarcation point.

The invention has the advantages that:

(1) The method adopts a combined coding mode to automatically judge the pinch-out and the lenticular stratum, and can effectively avoid the interference of shallow drilling (drilling of all stratum is not disclosed) on the pinch-out and lenticular stratum judgment.

(2) The position of the pinch-out and the thickness of the stratum of the lens body is calculated in a mode of spatial interpolation and spatial curve fitting, so that the modeling precision of the pinch-out and the stratum of the lens body is improved.

(3) The formation type judgment, formation boundary and formation body establishment are all automated, so that the modeling efficiency is greatly improved.

Drawings

FIG. 1 is a flow chart of an automated modeling method of a complex three-dimensional geologic model based on BIM.

FIG. 2 is a schematic diagram of a borehole and formation in accordance with an embodiment of the present invention.

FIG. 3 is a schematic representation of formation coding in accordance with an embodiment of the present invention.

FIG. 4 is a table of encoding points according to an embodiment of the present invention.

Fig. 5 is an initial state of the code pair-merging code pair group in the embodiment of the present invention.

FIG. 6 illustrates the final state of the code pair-combining code pair group in an embodiment of the present invention.

Fig. 7 is a schematic diagram of determining a pinchoff formation according to an embodiment of the present invention.

FIG. 8 is a diagram illustrating the positional relationship of lens layers according to an embodiment of the present invention.

FIG. 9 is a diagram illustrating a lens formation positional relationship according to an embodiment of the present invention.

FIG. 10 illustrates a terrain and formation interface to be modified in accordance with an embodiment of the present invention.

FIG. 11 is a schematic diagram illustrating the generation of upper and lower formation interface curves for a lens body formation integrally located within a conventional formation in accordance with an embodiment of the present invention.

FIG. 12 is a schematic representation of the generation of upper and lower formation interface curves for a lens formation integrally located adjacent to two conventional formations in an embodiment of the present invention.

FIG. 13 is a schematic representation of the generation of an upper and lower formation interface curve for an immediate vicinity of a pinch-out formation in an embodiment of the present invention.

The reference numerals are as follows:

1. an orifice point; 2. a formation demarcation point; 3. a hole bottom point; 4. drilling holes; 5. shallow drilling; 6.I formation; a formation ii; a iii formation; iv formation; a formation of 10. V; a vi formation; a VII formation; a viii formation; a ix formation; a formation x; 16. coding a stratum; 17. top encoding; 18. bottom coding; 19. drilling a hole number; 20. coding points; 21. a code pair; 22. an overall outer contour; JT outer contour; 24. cutting off the length; JT drilling; 26. a terrain surface; i, stratum interface to be modified under stratum; IV, the stratum interface to be modified on the stratum; IV, stratum interface to be modified under stratum; VII, forming a stratum interface to be modified on the stratum; the stratum interface to be modified under the VII stratum; a stratum interface to be modified on the stratum X; a limited lens formation interface above the vi formation; vi lower formation limited lens formation interface; a limited pinch-out formation interface curve on a ii formation; a lower limit pinch-out formation interface curve for a ii formation; a limited pinch-out formation interface curve on a iii formation; III, a lower limit pinch-out formation interface curve; a v-formation upper limited lens formation interface; a lower-formation-limited lens formation interface of 40; 41. cutting lines; 42. tangential line intersection points; 43. an enclosing region; 44. vertical cutting surfaces; 45. a transition end point; 46. transition point vanishing point.

Detailed Description

The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.

As shown in FIG. 1, the automatic modeling method of the complex three-dimensional geological model based on BIM comprises the following steps:

step 1: and acquiring the space coordinates of the orifice point 1, the space coordinates of the stratum demarcation point 2, the space coordinates of the hole bottom point 3 and stratum numbers among the points of the drilling hole 4, recoding the stratum numbers, and then forming and storing the stratum codes 16 and the coding points 20 of the drilling hole 4.

In step 1, the stratum numbers are recoded, and the stratum codes 16 which are the same and unique for the same stratum; the encoding points 20 include spatial coordinates and upper and lower side codes, the formats are (X, Y, Z, upper side code Noup, lower side code Nodown), and the encoding points 20 are configured into the following three forms according to the difference of the spatial positions of the encoding points 20:

If the code point 20 is located at the aperture point 1, the spatial coordinates of the code point 20 are the spatial coordinates of the aperture point 1, the upper code Noup of the code point 20 is a code which is different from and fixed to all the stratum codes 16, and is marked as a TOP code 17, the TOP code 17 is TOP, and the lower code Noup of the code point 20 is the stratum code 16 at the lower side of the aperture point 1;

if the coding point 20 is located at the formation boundary point 2, the spatial coordinates of the coding point 20 are the spatial coordinates of the formation boundary point 2, and the upper and lower side codes Noup and Nodown of the coding point 20 are the formation codes 16 on the upper and lower sides of the formation boundary point 2;

if the code point 20 is located at the hole bottom point 3, the spatial coordinates of the code point 20 are those of the hole bottom point 3, the upper code Noup of the code point 20 is the formation code 16 above the hole bottom point 3, the lower code noun of the code point 20 is a code which is different from and fixed to all the formation codes 16 and the top code 17, and is denoted as a bottom code 18, and the bottom code 18 is DOWN.

Fig. 2 and 3 show schematic views of a borehole and a formation according to an embodiment of the present invention, and fig. 3 shows schematic views of formation coding and a borehole numbering according to an embodiment of the present invention. Wherein, 4 is a drilling hole, the drilling hole number 19 is ZK1, ZK2, ZK3, ZK4, ZK5, ZK6, ZK7 and ZK8 in sequence, and the drilling hole 4 in the embodiment of the invention also comprises a shallow drilling hole 5;1 is the orifice point of the borehole 4; 2 is the formation demarcation point in borehole 4; 3 is the bottom point of the drilled hole 4; 6. 7, 8, 9, 10, 11, 12, 13, 14, 15 are the formations between the points, and formation code 16 is I, II, III, IV, V, VI, VII, VIII, IX, X in turn.

The encoding points 20 of all the boreholes in the embodiment of the invention are shown in fig. 4. The storage is to store the coding points 20 in the same drilling hole 4 from high to low according to the height Z and corresponding to the drilling hole number 19.

Step 2: and comparing all the drilled coding points 20 from top to bottom, and acquiring pinch-out and lenticular stratum codes.

In step 2, the coding points 20 of all the boreholes are compared from top to bottom to obtain the pinch-out and lenticular stratum codes, and the flow is as follows:

step 2.1: copying the upper and lower codes of all the code points 20 according to the storage sequence of the code points 20 in the step 1 to form code pair-merging code pair groups, wherein the storage form of each code pair 21 of the code pair-merging code pair groups is Noup, nodown, and the initial state of the code pair-merging code pair groups in the embodiment of the invention is shown in figure 5; setting i=1;

step 2.2: acquiring the ith code pair 21 from top to bottom in each drilling hole in the code pair-merging code pair group, grouping the code pairs 21 according to different lower side codes of the acquired code pairs 21, and ensuring that each code pair 21 in the same group has the same lower side code; screening out the packets of which the lower side is coded as a bottom code 18, and eliminating corresponding original code pairs 21 from the code pair-combined code pair group, and recording the rest of the packets as code pair packets; if the number of the code pairs is more than 1, executing the step 2.3; if the number of the codes is equal to 1, executing the step 2.5; if the number of the codes is equal to 0, executing the step 2.6;

Step 2.3: selecting one of the code pair groups, carrying out merging code pair operation, repeating the steps until the lower side code of one non-merging code pair group is the same as the lower side code of the merging code pair group, replacing the original code pair 21 with the merged code pair 21 in the code pair-merging code pair group, and simultaneously marking the non-merging code pair group as a main round reference code pair group;

step 2.4: taking the reference code pair group of the round as a reference, carrying out combination code pair operation on the code pair groups except the code pair group selected in the last time in the step 2.3, and replacing the original code pair 21 by the combined code pair 21 in the code pair-combination code pair group;

the steps 2.3 and 2.4 combined coding pair operation steps are as follows:

step 2a: setting j=1, acquiring all code pairs 21 of which the operation code pairs are required to be combined and the corresponding drill holes of the code pair group are numbered 19 in the code pair-combined code pair group, and establishing a selected code pair extension group;

step 2b: extracting the ith code pair 21 and the (i+j) th code pair 21 of the same drilling number 19 in the selected code pair extension group, and combining to generate a temporary combined code pair group; the rule of merging merges the upper side code of the i-th code pair 21 and the lower side code of the i+j-th code pair 21;

Step 2c: if the temporary combination coding pair group has the same temporary combination coding pair 21 at the upper side and the lower side, the original coding pair 21 is eliminated from the selected coding pair extension group and the coding pair-combination coding pair group in the step 2.1, the step 2b is returned, and j is set to 0;

step 2d: removing the temporary combination code pairs 21 with the same lower codes as the bottom codes 18, removing all code pairs 21 of the drill holes 4 where the temporary combination code pairs 21 are positioned from the selected code pair extension group, and removing the original code pairs 21 from the code pair-combination code pair group in the step 2.1;

step 2e: reclassifying the temporary combination code pairs 21 in the temporary combination code pair group, wherein the grouping rule is that the temporary combination code pairs 21 with the same lower side code are classified into the same group, and the same group is marked as a temporary combination code pair sub-group; if the number of the sub-packets is greater than or equal to 1, executing the step 2f, otherwise executing the step 2h;

step 2f: comparing the temporary combination coding sub-packet with the appointed coding sub-packet, removing all coding pairs 21 in the drill hole 4 where the temporary combination coding sub-packet which is the same as the appointed coding sub-packet is positioned from the selected coding pair extension packet, removing the temporary combination coding sub-packet from the temporary combination coding sub-packet, and simultaneously replacing the original coding pair 21 in the coding pair-combination coding pair group in step 2.1, and returning a matching success signal; the specified code pair group is an uncombined code pair group when the step 2.3 is executed, and is a reference code pair group of the round when the step 2.4 is executed;

step 2h: if all code pairs 21 in the selected code pair group are removed and the code pair group is not successfully matched with the designated code pair group, returning a matching unsuccessful signal; if the matching is successful, a success signal is returned;

step 2.6: stratigraphic codes 16 that are not included in the final encoding pair-merge encoding pair set are noted as pinch-out and lenticular stratigraphic codes, and the remaining stratigraphic codes 16 are noted as regular stratigraphic codes;

for the embodiment of the invention, the process of acquiring the pinch-out and lenticular stratum codes is as follows:

1) Setting i=1, and acquiring an ith coding pair of each drilling hole in the coding pair-merging coding pair group; forming a code pair packet (TOP, I);

the number of the code pairs is 1, i is increased by 1;

2) Acquiring an ith code pair of each drilling hole in the code pair-combined code pair group; forming code pairs

The number of the code pair packets is 2;

3) Setting j=1, carrying out combined coding pair operation on the selected coding pair group (I, II), obtaining all coding pairs of corresponding drilling holes (ZK 1-ZK 5) of the coding pair group, and establishing a selected coding pair extension group;

4) Combining the ith code pair and the (i+j) th code pair corresponding to the same drilling number in the selected code pair extension group to generate a temporary combined code pair group;

the temporary combination coding pair group does not have the same temporary combination coding pair of the upper side and the lower side; the presence of ZK3 in the temporal combining coding pair packet is the same temporal combining coding pair (I, DOWN) as the bottom coding;

5) Rejecting (I, DOWN) in the temporary combination code pair packet, rejecting all code pairs of the temporary combination code pair corresponding borehole (ZK 3) from the selected code pair extension packet, and rejecting (I, II) and (II, DOWN) code pairs in ZK3 of the code pair-combination code pair group; re-grouping the remaining ones of the temporary combination code pairs in the group to form (I, III) and (I, IV) temporary combination code pairs sub-groups;

6) The temporary combination coding pair group has the same sub-group as the appointed coding pair group (I, IV), all coding pairs of the corresponding drilling holes (ZK 4 and ZK 5) of the sub-group are removed from the selected coding pair extension group, the sub-group is removed from the temporary combination coding pair group, and meanwhile, one combination coding pair (I, IV) is used for replacing two original coding pairs (I, II) and (II, IV) of the ZK4 and the ZK5 in the coding pair-combination coding pair group; returning a matching success signal;

The temporary merging code pair group also has 1 sub-group, j is increased by 1;

7) Combining the ith code pair and the (i+j) th code pair corresponding to the same drilling number in the selected code pair extension group to generate a temporary combined code pair group;

the temporary combination coding pair group does not have the same temporary combination coding pair of the upper side coding and the lower side coding and the same temporary combination coding pair of the lower side coding and the bottom coding;

8) The temporary combination coding pair group has the same sub-group as the designated coding pair group (I, IV), all coding pairs of the corresponding drill holes (ZK 1 and ZK 2) of the sub-group are eliminated from the selected coding pair extension group, the sub-group is eliminated from the temporary combination coding pair group, and simultaneously, one combination coding pair (I, IV) is used for replacing three original coding pairs (I, II), (II, III) and (III, IV) of ZK1 and ZK2 in the coding pair-combination coding pair group; returning a matching success signal;

temporarily combining the code pair group without sub-groups, and determining the code pair group (I, IV) as the reference code pair group of the round according to the success signal; no packet which still needs to be subjected to merging and coding pair operation exists, i is increased by 1;

9) Acquiring an ith code pair of each drilling hole in the code pair-combined code pair group; forming code pairs of packets (IV, VII), (IV, VI) and (IV, V);

The number of the code pair groups is 3;

10 Setting j=1, carrying out combined coding pair operation on the selected coding pair group (IV, VII), obtaining all coding pairs of the corresponding drilling holes (ZK 1-2), and establishing a selected coding pair extension group;

performing merging coding pair operation until the temporary merging coding pair packet is empty and a success signal is not returned;

11 Setting j=1, carrying out combined coding pair operation on the selected coding pair groups (IV, VI), obtaining all coding pairs of corresponding drilling holes (ZK 4-5), and establishing a selected coding pair extension group;

12 Combining the ith code pair and the (i+j) th code pair corresponding to the same drilling number in the selected code pair extension group to generate a temporary combined code pair group;

the temporary combination coding pair has the same temporary combination coding pair of upper and lower side codes in the group;

13 Temporary merge coding pair (IV, IV) is present for ZK4, original coding pairs (IV, VI) and (VI, IV) are culled from ZK4 of the selected coding pair extension group and coding pair-merge coding pair group, and j=0 is set;

14 Combining the ith code pair and the (i+j) th code pair corresponding to the same drilling number in the selected code pair extension group to generate a temporary combined code pair group; re-grouping the temporary combining code pairs in the temporary combining code pair group to form (IV, VII) and (IV, VI) temporary combining code pair sub-groups;

15 The same sub-packet as the specified code pair packet (IV, VII) exists in the temporary merge code pair packet, all code pairs of the sub-packet corresponding borehole (ZK 4) are removed from the selected code pair extension packet, and the sub-packet is removed from the temporary merge code pair packet while replacing (IV, VII) of ZK4 with one merge code pair (IV, VII) in the code pair-merge code pair group; returning a matching success signal;

16 Combining the ith code pair and the (i+j) th code pair corresponding to the same drilling number in the selected code pair extension group to generate a temporary combined code pair group;

17 The same sub-packet as the specified code pair packet (IV, VII) exists in the temporary merge code pair packet, all code pairs of the sub-packet's corresponding borehole (ZK 5) are removed from the selected code pair extension packet, and the sub-packet is removed from the temporary merge code pair packet while the two original code pairs (IV, VI) and (VI, VII) of ZK5 are replaced with one merge code pair (IV, VII) in the code pair-merge code pair group; returning a matching success signal;

Temporarily combining the code gambling groups without sub-groups, and determining the code pair groups (IV, VII) as the round of reference code pair groups according to the success signals;

18 A merging coding pair operation is carried out on the coding pair group (IV, V), and the coding pairs (IV, V), (V, IV), (IV, VI) and (VI, IV) are removed from the ZK6 of the coding pair-merging coding pair group, and the coding pairs (IV, V) and (V, IV) are removed from the ZK7 and the ZK8 after the operation is finished;

19 Continuing to perform step 2 operations on the remaining pairs of pairs in the set of pairs of codes-the (X, DOWN) pairs of codes ZK1 and ZK4-ZK8 being removed, the (VII, DOWN) pairs of codes ZK2 being removed, the (VII, VIII) and (VIII, VII) pairs of codes ZK5-ZK7 being removed, the (VII, IX) and (IX, X) pairs being replaced by a combined pair of codes (VII, X) in ZK6-ZK 8;

20 The codes included in the final code pair-merge code pair group are TOP, I, IV, VII and X; determining that the formations encoded as II, III, V, VI, VIII and IX are the pinch-out and lens formation encodings, and I, IV, VII and X are the conventional formation encodings;

in the embodiment of the present invention, the final states of the coding pair-merging coding pair group are shown in fig. 6, and as can be seen from fig. 6, the pinch-out and lenticular layer coding includes: II. III, V, VI, VIII, IX, conventional formation coding includes: I. IV, VII, X.

Step 3: according to the pinch-out and lenticular stratum codes, obtaining a drilling hole 4 with the pinch-out and lenticular stratum codes, which is marked as JT drilling hole 25, and comparing the JT drilling hole 25 with the rest of the drilling holes 4 to distinguish the pinch-out stratum codes and the lenticular stratum codes; as can be seen from fig. 5, in the embodiment of the present invention, the JT borehole 25 includes: ZK1, ZK2, ZK3, ZK4, ZK5, ZK6, ZK7, ZK8;

in step 3, the process of distinguishing the pinch-out stratum from the lenticular stratum is as follows:

step 3.1: acquiring plane coordinates (X, Y) of all drilling holes 4 to generate plane coordinate points, and drawing a closed outer contour line according to the plane points of all the drilling holes 4, namely an integral outer contour line 22;

step 3.2: acquiring the pinch-out and lenticular layer codes in the step 2, which are called JT codes, wherein the JT codes comprise II, III, V, VIII, IX, and selecting one JT code; traversing all the drilling holes 4, obtaining horizontal coordinates (X, Y) of corresponding JT drilling holes 25 to generate plane coordinate points, and drawing a closed outer contour line, namely JT outer contour line 23;

the closed outer contour line meets the requirement that all plane points participating in constructing the closed outer contour line are positioned in or on the closed outer contour line, and the area is maximum;

step 3.3: according to the positional relationship between the JT outer contour line 23 and the integral outer contour line 22, JT codes corresponding to the JT outer contour line 23 under the following two positional relationships are determined as lenticular layer codes:

Positional relationship one: JT outer contours 23 are all located within the overall outer contour 22;

and the position relation is II: JT outer contour line 23 partially coincides with integral outer contour line 22, and when the offset from the overlapping position to the non-overlapping position is satisfied, the trend of increasing and then decreasing is presented;

the determination method of the trend of increasing and then decreasing in the second positional relationship is as follows:

step 3a: determining the overlapping part of the JT outer contour line 23 and the integral outer contour line 22, and drawing the overlapping contour line according to whether the overlapping part is continuous on the integral outer contour line 22, specifically: if continuous, drawing a coincident contour line; if the line is discontinuous, the line is disconnected at the discontinuous position, and a plurality of coincident contour lines are drawn;

step 3b: selecting a coincident contour line, extracting two endpoints of the selected coincident contour line to draw a straight line, and obtaining the length of the straight line cut by the JT outer contour line 23, which is called cut length 24;

step 3c: shifting the straight line to the inner side of the JT outer contour line 23 for a plurality of times, and acquiring a cut-off length 24 once for each shift until the straight line is stopped when the straight line does not intersect with the JT outer contour line 23;

step 3d: acquiring all offset distances and cut-off lengths 24, wherein the cut-off lengths 24 show a trend of increasing and then decreasing along with the increase of the offset distances, and then a success signal is returned;

Step 3e: repeating the steps 3b-3d until all the coincident contour lines return a success signal in the step 3d, and determining that the second position relation is met;

fig. 7-9 show, fig. 7 is a schematic diagram of determining a pinchoff stratum in an embodiment of the present invention, fig. 8 is a schematic diagram of a lens layer position relationship in an embodiment of the present invention, and fig. 9 is a schematic diagram of a lens layer position relationship in an embodiment of the present invention. The pinch-out formation coding includes: II. III, IX; lenticular formation encoding includes: v, VI, VIII;

step 4: picking up a conventional stratum code, and establishing a terrain, a conventional stratum interface and a stratum interface to be modified by adopting an interpolation algorithm according to the conventional stratum code pick-up code point 20;

in step 4, the process of establishing the terrain, the conventional stratum interface and the stratum interface to be modified is as follows:

step 4.1: selecting a conventional formation code, the conventional formation code comprising: I. IV, VI, VII, X, traversing all the boreholes 4 to obtain the coding points 20 of the conventional stratum coding up and down;

the upper and lower encoding points 20 of the conventional formation code are obtained as follows:

If two encoding points 20 encoded by the conventional stratum are contained in one borehole 4, sorting according to the height values, determining an upper encoding point 20 with the largest height value as the stratum, and determining a lower encoding point 20 with the smallest height value as the stratum;

if the number of the coding points 20 of the conventional stratum coding 16 in one borehole 4 is greater than two, the coding points are marked as a coding point group A, the upper coding points 20 with the largest elevation value and the lower coding points 20 with the smallest elevation value are determined according to the magnitude of the elevation value, and the coding points are corrected according to the situation, and the method is as follows: acquiring the lenticular stratum codes in the code point group A, searching the code points 20 containing the lenticular stratum codes 16 in all the drill holes 4, and marking the code points as a code point group B, wherein if only the lenticular stratum and the conventional stratum codes are in the code point group B, the upper and lower code points 20 of the conventional stratum are not corrected; if the coding point group B contains the coding points 20 of the adjacent regular stratum on the regular stratum, modifying the upper coding points 20 of the regular stratum into the coding points 20 with the lowest coding elevation values of the lens stratum in the drill hole 4; if the coding point group B contains the coding points 20 of the adjacent regular stratum under the regular stratum, modifying the lower coding points 20 of the regular stratum to be the coding points 20 with the highest coding elevation values of the lens stratum in the drilling hole 4;

Step 4.2: the position of the coding point 20 is obtained according to the step 4.1, and the terrain, the conventional stratum interface and the stratum interface to be modified are built, specifically as follows: if the position of the upper coding point 20 coincides with the aperture point 1, the topographic surface 26 is directly established by a topographic map data spatial interpolation method; if the spatial positions of the upper and lower coding points 20 coincide with the stratum demarcation point 2, respectively establishing an upper stratum interface and a lower stratum interface to be trimmed of the stratum by a spatial interpolation method on the basis of the spatial coordinates of the upper and lower coding points 20; if the position of the lower coding point 20 coincides with the hole bottom point 3, no surface is generated;

The adjacent conventional stratum is judged by the coding pair-merging coding pair group obtained in the step 2, and the judging method is as follows: acquiring the selected regular stratum code and traversing each code pair 21 of the final code pair-merging code pair group, and if the selected regular stratum code appears as Nodown in the code pair 21, taking the stratum corresponding to Noup in the code pair 21 as an upper adjacent regular stratum; if the regular formation code is selected to appear as the Noup in the code pair 21, the formation corresponding to the noun in the code pair 21 is taken as the next adjacent regular formation; the final coding pair-merging coding pair group is a coding pair-merging coding pair group which completes all merging operations through the step 2;

FIG. 10 is a schematic view of the terrain and the formation interface to be modified according to an embodiment of the present invention, wherein 26 is a terrain surface; 27 is the stratum interface to be modified under the stratum I; 28 is the stratum interface to be modified on the IV stratum; 29 is the stratum interface to be modified under the IV stratum; 30 is the stratum interface to be modified on the VII stratum; 31 is the stratum interface to be modified under the VII stratum; 32 is the formation boundary to be modified on the X formation.

Step 5: picking up a corresponding stratum demarcation point 2 according to the stratum code of the lens body, and establishing the stratum body of the lens body;

In step 5, the process of establishing the lenticular stratum is as follows:

step 5.1: selecting a lens stratum code, traversing all drill holes 4 to obtain upper stratum demarcation points 2 and lower stratum demarcation points 2 of the lens stratum; taking the JT outer contour line 23 in the step 3.2 as a boundary, respectively performing spatial interpolation on the upper and lower stratum boundary points 2 of the lens stratum, and establishing upper and lower stratum interfaces of the lens stratum with limited range, which are called upper and lower limited lens stratum interfaces;

step 5.2: taking the centroid of a closed polygon formed by the JT outer contour line 23 as a starting point, making a plurality of radial lines outwards, and simultaneously stretching the radial lines vertically up and down to form a plurality of vertical tangential planes 44; sectioning the upper and lower limited lens stratum interfaces by adopting the vertical sections 44, and obtaining upper and lower limited lens stratum interface curves in each vertical section 44; one end, which coincides with the JT outer contour line 23, of the vertical projection of the stratum interface curve of the lens body is recorded as a transition end point 45;

if the whole lens stratum is positioned in a certain conventional stratum, solving tangent lines 41 of the upper and lower limited lens stratum interface curves at transition end points 45 of the tangent lines in each vertical tangent plane 44, and acquiring the intersection point of the two tangent lines 41; constructing a triangle enclosing area 43 by using the tangent intersection point 42 and the two transition end points 45, searching a point with zero thickness of the lens body in the enclosing area 43, and establishing complete upper and lower lens body stratum interface curves in the vertical tangent plane 44;

If the lenticular stratum is positioned at the adjacent position of the two conventional strata, generating an intersection line of the vertical tangent plane 44 and the upper and lower stratum interfaces to be trimmed of the two conventional strata at the position in each vertical tangent plane 44, generating a straight line through two transition end points 45, searching a point with zero lenticular thickness in an intersection line enclosing area 43 of the straight line and the upper and lower stratum interfaces to be trimmed, and establishing complete upper and lower lenticular stratum interface curves in the vertical tangent plane 44;

the point with zero thickness of the searching lens body needs to meet the requirement that two complete upper and lower lens body stratum interface curves fitted by the point and discrete points on the upper and lower limited lens body stratum interface curves have minimum deviation from the upper and lower limited lens body stratum interface curves; the fitting curve is a spline curve or a polynomial curve;

step 5.5: repeating steps 5.1-5.4 to generate all lens stratum;

11-12, FIG. 11 is a schematic diagram illustrating the generation of upper and lower formation interface curves of a lens formation integrally located inside a conventional formation in an embodiment of the present invention, and FIG. 12 is a schematic diagram illustrating the generation of upper and lower formation interface curves of a lens formation integrally located adjacent to two conventional formations in an embodiment of the present invention, wherein 33 is a limited lens formation interface on a VI formation; 34 is the VI formation lower limit lens formation interface; 39 is the upper limit lens body stratum interface of the V stratum; 40 is the lower limit lens stratum interface of the V stratum; 41 is a tangent line; 42 is the tangent intersection; 43 is the enclosed area; 44 is a vertical section; 45 is the transition endpoint;

step 6: picking up a corresponding stratum demarcation point 2 according to the pinch-out stratum code, and establishing a pinch-out stratum body;

in step 6, the process of establishing the pinch-out stratum is as follows:

step 6.1: selecting a to-be-built pinch-out stratum code, traversing all drill holes 4 to obtain upper and lower stratum demarcation points 2 of the pinch-out stratum; taking the JT outer contour line 23 in the step 3.2 as a boundary, respectively carrying out spatial interpolation on an upper stratum demarcation point 2 and a lower stratum demarcation point 2 of a to-be-built pinch-out stratum, and establishing upper stratum interfaces and lower stratum interfaces of the pinch-out stratum with limited range, which are called upper limited pinch-out stratum interfaces and lower limited pinch-out stratum interfaces;

Step 6.2: in the step 3, at the position where the JT outer contour line 23 does not coincide with the integral outer contour line 22, discrete points are extracted according to the mesh set size, and a horizontal straight line perpendicular to the JT outer contour line 23 is made with the discrete points as vertical points; stretching the straight line up and down along the vertical direction to form a vertical tangent plane 44, and sectioning upper and lower limited pinch-out stratum interfaces by adopting the vertical tangent plane 44, wherein upper and lower limited pinch-out stratum interface curves can be obtained in each vertical tangent plane 44; one end of the vertical projection of the limited pinch-out stratum interface curve, which coincides with the JT outer contour line 23, is recorded as a transition tip vanishing point 46;

step 6.3: the intersection line of the upper and lower stratum interfaces of the upper and lower adjacent conventional strata of the to-be-built pinch-out stratum at the position is regenerated in each vertical tangent plane 44, and a complete upper and lower pinch-out stratum interface curve is generated according to whether the upper and lower pinch-out strata of the to-be-built pinch-out stratum have adjacent and built pinch-out stratum bodies or not, specifically as follows:

if the built pinch-out stratum exists in the upper and lower adjacent strata of the pinch-out stratum to be built, an upper pinch-out stratum interface curve and a lower pinch-out stratum interface curve of the built pinch-out stratum are generated in the vertical tangent plane 44, if the built pinch-out stratum interface curve exceeds the range of the transition pinch-out points 46, one or two of the upper and lower limited stratum interface curves of the pinch-out stratum to be built are replaced by the built pinch-out stratum interface curve on the premise of keeping the transition pinch-out points 46 unchanged, and the pinch-out points are searched in the intersection line enclosing area 43 of the straight line formed by the transition pinch-out points 46 and the upper and lower stratum interface to be trimmed, and meanwhile, the complete upper and lower pinch-out stratum interface curves are built; the method for replacing the upper and lower limited stratum interface curves of the stratum to be pinch-out by the established pinch-out stratum interface curves comprises the following steps: if the pinch-out stratum to be built is adjacent to the pinch-out stratum to be built and meets the substitution requirement, substituting an upper limited pinch-out stratum interface curve of the pinch-out stratum to be built by a lower pinch-out stratum interface curve of the pinch-out stratum to be built which is adjacent to the pinch-out stratum to be built; if the pinch-out stratum to be built is adjacent to the built pinch-out stratum, and the substitution requirement is met, substituting the lower limited pinch-out stratum interface curve of the pinch-out stratum to be built with the upper pinch-out stratum interface curve of the built pinch-out stratum which is adjacent to the lower pinch-out stratum;

If the stratum immediately adjacent to the upper stratum and the lower stratum is not the pinch-out stratum or the pinch-out stratum is not established, generating a straight line through two transition pinch-out points 46, searching the pinch-out points in a crossing line enclosing area 43 of the straight line and the upper stratum interface to be modified and the lower stratum interface to be modified, and establishing a complete upper pinch-out stratum interface curve and a complete lower pinch-out stratum interface curve;

the searching point is required to meet the requirement that two complete upper and lower pinch-out stratum interface curves fitted by the point and discrete points on the upper and lower confined pinch-out stratum interface curves have minimum deviation from the upper and lower confined pinch-out stratum interface curves; the fitting curve is a spline curve or a polynomial curve;

FIG. 13 is a schematic diagram of the generation of upper and lower formation interface curves of a formation immediately adjacent to a pinch-out in an embodiment of the present invention, wherein 35 is a limited pinch-out formation interface curve of a formation II; 36 is the lower limit pinch-out formation interface curve of the II formation; and 37 is the limiting pinch-out stratum interface curve on the III stratum; 38 is the lower limit pinch-out formation interface curve for the III formation; 43 is the enclosed area; 44 is a vertical section; 46 is the transition point.

Step 7: correcting a stratum interface to be trimmed according to the pinch-out and the lens stratum, and establishing a conventional stratum;

in step 7, the specific process of correcting the stratum interface to be corrected and establishing the stratum body is as follows:

step 7.1: acquiring a group of formation interfaces to be modified in the step 4.3 and generating all coding points 20 of the formation interfaces to be modified, and judging the modification mode of the group of formation interfaces to be modified according to the coding points 20, wherein the method comprises the following specific steps:

if lenticular layer encoding occurs in the encoding points 20, performing a trimming mode one;

if the pinch-out layer coding occurs in the coding points 20, a second trimming mode is executed;

if the lenticular and pinch-out layer codes occur simultaneously in the code point 20, a finishing mode three is performed;

the first trimming mode is as follows: removing the coding points 20 containing the lenticular stratum codes 16, reserving only one of the coding points 20 with coincident space coordinates, and generating a conventional stratum interface by adopting a space interpolation algorithm;

the second trimming mode is as follows: searching stratum pinch-out bodies established in the step 6 according to the pinch-out stratum codes 16, and extracting pinch-out boundaries of the stratum pinch-out bodies; extracting discrete points on the pinch-out boundary according to the set size of the grid to serve as pinch-out additionally-arranged stratum demarcation points; incorporating within each borehole 4 a code point 20 containing the pinch-out formation code as a dummy formation demarcation point; reserving one of the coding points 20 with the coincident space coordinates, combining the additional formation demarcation points with the dummy formation demarcation points, and generating a conventional formation interface through a space interpolation algorithm;

The method comprises the steps of extracting the pinch-out boundary of a stratum pinch-out body according to the form of the stratum pinch-out body, and specifically comprises the following steps: if the pinch-out stratum body does not exist in the stratum pinch-out bodies which are adjacent to each other from top to bottom, the position where the thickness of the pinch-out stratum body is zero is extracted as a pinch-out boundary; if the pinch-out stratum body exists in the upper and lower adjacent pinch-out stratum bodies, extracting the position where the sum of the thicknesses of all the adjacent pinch-out stratum bodies is zero as a pinch-out boundary;

and in the second finishing mode, the virtual stratum demarcation points are combined according to the stratum pinch-out body mode, and the method specifically comprises the following steps: if the pinch-out formation body does not exist in the formation pinch-out bodies which are immediately adjacent to each other, solving a space coordinate average value of the coding points 20 comprising the pinch-out formation codes 16 in each borehole 4 as a dummy formation demarcation point; if the pinch-out formation exists in the upper and lower adjacent pinch-out formations, solving a space coordinate average value of the coding points 20 comprising the pinch-out formation and the upper and lower adjacent formation codes in each borehole 4 as a dummy formation demarcation point;

the third trimming mode is as follows: removing the coding points 20 containing the lens body codes by adopting a first trimming mode, and establishing a conventional stratum interface according to a second trimming mode;

Step 7.2: the method comprises the steps of establishing a conventional stratum body to be modified according to a conventional stratum interface, wherein the establishment mode is formed by surrounding upper and lower adjacent conventional stratum interfaces or terrain surfaces 26;

step 7.3: checking whether all the regular stratum bodies to be corrected are intersected with the pinch-out and lens stratum bodies, and subtracting the pinch-out and lens stratum bodies from the regular stratum bodies through Boolean operation if the regular stratum bodies are intersected to obtain a final regular stratum body;

the above embodiments are merely preferred embodiments of the present invention and are not intended to limit the present invention, and any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of the present invention.

Claims

1. The automatic modeling method of the complex three-dimensional geological model based on BIM is characterized by comprising the following steps of:

step 1: acquiring space coordinates of an orifice point (1), space coordinates of a stratum demarcation point (2) and space coordinates of a hole bottom point (3) in a drill hole (4), acquiring stratum codes (16) between the points, and forming corresponding coding points (20) according to the space coordinates of the points and the stratum codes (16) between the points;

the encoding points (20) comprise the spatial coordinates (X, Y, Z) of the points and the upper encoding Noup and the lower encoding Nodown in the format (X, Y, Z, noup, nodown); the upper and lower side codes are the upper and lower stratum codes (16) of the point, wherein the upper side code of the orifice point (1) corresponding to the code point (20) is the top code (17), and the lower side code of the orifice bottom point (3) corresponding to the code point (20) is the bottom code (18);

Step 2: comparing the coding points (20) of all the drill holes (4) from top to bottom to obtain pinch-out and lenticular stratum codes;

in the step 2, extracting upper and lower side codes of all the coding points (20) to form a coding pair-merging coding pair group, wherein each coding pair (21) in the coding pair-merging coding pair group has the format of (Noup, nodown); combining the coding pairs (21) in the drilling holes (4), combining the upper coding of the ith coding pair (21) and the lower coding of the (i+j) th coding pair (21) to obtain combined ith coding pair (21), combining to make the lower coding of the (21) th coding pair from top to bottom in each drilling hole (4) identical, and marking stratum codes (16) which are not included in the final coding pair-combined coding pair group as pinch-out and lenticular stratum codes and marking the rest stratum codes (16) as conventional stratum codes;

step 3: according to the pinch-out and lenticular stratum codes, obtaining a drilling hole (4) with the pinch-out and lenticular stratum codes, recording as a JT drilling hole (25), and comparing the JT drilling hole (25) with the rest of the drilling holes (4) to distinguish the pinch-out stratum codes and the lenticular stratum codes;

step 4: the method comprises the steps of (1) picking up coding points (20) according to conventional stratum coding, and establishing a terrain, a conventional stratum interface and a stratum interface to be modified;

step 4.1: selecting a conventional stratum code, traversing all drill holes (4) to obtain coding points (20) above and below the conventional stratum code;

step 4.2: and (3) according to the step 4.1, acquiring the position of the coding point (20), and establishing a terrain, a conventional stratum interface and a stratum interface to be modified, wherein the method comprises the following steps: if the position of the upper coding point (20) coincides with the orifice point (1), a topographic surface (26) is directly established by a topographic map data spatial interpolation method; if the space positions of the upper and lower coding points (20) coincide with the stratum demarcation points (2), respectively establishing an upper stratum interface and a lower stratum interface to be trimmed of the stratum by a space interpolation method on the basis of the space coordinates of the upper and lower coding points (20); if the position of the lower coding point (20) coincides with the hole bottom point (3), no surface is generated;

Step 5: picking up corresponding stratum demarcation points (2) according to the lenticular stratum codes, and establishing lenticular stratum bodies;

step 6: picking up corresponding stratum demarcation points (2) according to the pinch-out stratum codes, and establishing pinch-out stratum bodies;

2. The automated modeling method of a complex three-dimensional geologic model based on BIM according to claim 1, wherein in step 2, the steps of obtaining the pinch-out and lenticular formation codes are as follows:

step 2.1: the coding points (20) in the same drilling hole (4) are ordered from high to low according to the elevation Z; extracting upper and lower side codes of all code points (20) to form code pair-merging code pair groups, wherein each code pair (21) in the code pair-merging code pair groups has the format of (Noup, nodown); and setting i=1;

step 2.2: acquiring an ith code pair (21) from top to bottom in each drilling hole (4), grouping the code pairs (21) according to different lower side codes of the acquired code pairs (21), and ensuring that each code pair (21) in the same group has the same lower side code; screening out the packets of which the lower side codes are bottom codes (18), and eliminating corresponding original code pairs (21) in code pair-combination code pair groups, and recording the rest of the packets as code pair packets; if the number of the code pairs is more than 1, executing the step 2.3; if the number of the codes is equal to 1, executing the step 2.5; if the number of the codes is equal to 0, executing the step 2.6;

Step 2.3: selecting one of the code pair groups, carrying out merging code pair operation, repeating the steps until the lower side code of one non-merging code pair group is the same as the lower side code of the merging code pair group, replacing the original code pair (21) by a merged code pair (21) in the code pair-merging code pair group, and simultaneously marking the non-merging code pair group as a main round reference code pair group;

step 2.4: taking the reference code pair group of the round as a reference, carrying out combination code pair operation on the code pair groups except the code pair group selected in the last time in the step 2.3, and replacing the original code pair (21) by a combined code pair (21) in the code pair-combination code pair group;

step 2.6: the formation codes (16) not included in the final coding pair-merge coding pair group are noted as pinch-out and lenticular formation codes, and the remaining formation codes (16) are noted as regular formation codes;

step 2a: setting j=1, acquiring all code pairs (21) of the drill holes (4) corresponding to the code pair groups needing to be subjected to the combined code pair operation from the code pair-combined code pair group, and establishing a selected code pair expansion group;

Step 2b: extracting and combining an ith code pair (21) and an (i+j) th code pair (21) of the same drilling hole (4) in the selected code pair extension group to generate a temporary combined code pair group; the rule of merging merges the upper side code of the i-th code pair (21) and the lower side code of the i+j-th code pair (21);

step 2c: if the temporary combination coding pair group has the same temporary combination coding pair (21) of the upper side coding and the lower side coding, the original coding pair (21) is eliminated from the selected coding pair extension group and the coding pair-combination coding pair group in the step 2.1, the step 2b is returned, and j is set to 0;

step 2d: removing the temporary combination code pairs (21) with the same lower codes as the bottom codes (18), removing all code pairs (21) of the drill holes (4) where the temporary combination code pairs (21) are positioned from the selected code pair extension group, and removing the original code pairs (21) from the code pair-combination code pair group in the step 2.1;

step 2e: reclassifying the temporary combination code pairs (21) in the temporary combination code pair group, wherein the grouping rule is that the temporary combination code pairs (21) with the same lower side codes are classified into the same group, and marked as temporary combination code sub-groups; if the number of the sub-packets is greater than or equal to 1, executing the step 2f, otherwise executing the step 2h;

Step 2f: comparing the temporary combination coding sub-packet with the appointed coding sub-packet, removing all coding pairs (21) in the drill hole (4) where the same temporary combination coding sub-packet as the appointed coding sub-packet is located from the selected coding pair extension packet, removing the temporary combination coding sub-packet from the temporary combination coding sub-packet, and simultaneously replacing the original coding pair (21) in the coding pair-combination coding pair group in step 2.1, and returning a matching success signal; the specified code pair group is an uncombined code pair group when the step 2.3 is executed, and is a reference code pair group of the round when the step 2.4 is executed;

step 2h: returning a match unsuccessful signal if all code pairs (21) in the selected code pair group are removed and not successfully matched with the designated code pair group; if the matching is successful, a success signal is returned.

3. The automated modeling method of a complex three-dimensional geologic model based on BIM according to claim 1, wherein in step 3, the distinguishing between the pinch-out formation coding and the lenticular formation coding is as follows;

Step 3.1: obtaining point plane coordinates (X, Y) of all the drilling holes (4) to generate plane coordinate points, and drawing a closed outer contour line according to the plane coordinate points of all the drilling holes (4), namely an integral outer contour line (22);

step 3.2: acquiring the pinch-out and lenticular stratum codes in the step 2, namely JT codes, and selecting one JT code; traversing all the drilling holes (4), obtaining point plane coordinates (X, Y) of corresponding JT drilling holes (25) to generate plane coordinate points, and drawing a closed outer contour line according to the plane coordinate points of the JT drilling holes (25), namely JT outer contour line (23);

step 3.3: according to the position relation between the JT outer contour line (23) and the integral outer contour line (22), JT codes corresponding to the JT outer contour line (23) under the following two position relations are determined as lenticular stratum codes:

positional relationship one: JT outer contours (23) are all located within the overall outer contour (22);

and the position relation is II: JT outer contour line (23) is partially overlapped with integral outer contour line (22), and when the overlapping position is shifted to the non-overlapping position, the tendency of increasing and then decreasing is presented;

step 3.4: repeating steps 3.2 and 3.3 for all JT encodings; JT codes that are not determined to be lenticular formation codes are determined to be pinch-out formation codes.

4. The automated modeling method of a complex three-dimensional geological model based on BIM according to claim 1, wherein the closed outer contour in steps 3.1 and 3.2 satisfies that all the planar points involved in constructing the closed outer contour are located in or on the closed outer contour and have the largest area;

step 3a: determining the overlapping part of the JT outer contour line (23) and the integral outer contour line (22), and drawing the overlapping contour line according to whether the overlapping position is continuous on the integral outer contour line (22), specifically: if continuous, drawing a coincident contour line; if the line is discontinuous, the line is disconnected at the discontinuous position, and a plurality of coincident contour lines are drawn;

step 3b: selecting a coincident contour line, extracting two endpoints of the selected coincident contour line to draw a straight line, and obtaining the length of the straight line cut by the JT outer contour line (23), which is called cut length (24);

step 3c: shifting the straight line to the inner side of the JT outer contour line (23) for a plurality of times, and acquiring a cut-off length (24) once for each shift until the straight line is not intersected with the JT outer contour line (23);

step 3d: acquiring all offset distances and cut-off lengths (24), wherein the cut-off lengths (24) show a trend of increasing and then decreasing along with the increase of the offset distances, and returning a success signal;

5. The automated modeling method of a complex three-dimensional geological model based on BIM according to claim 1, wherein in step 4.1, the upper and lower encoding points (20) of the conventional formation code are obtained as follows:

if one borehole (4) contains two coding points (20) coded by the conventional stratum, sorting according to the height values, determining an upper coding point (20) with the largest height value as the stratum, and determining a lower coding point (20) with the smallest height value as the stratum;

if the number of the coding points (20) coded by the conventional stratum in one borehole (4) is more than two, the coding points are marked as a coding point group A, the upper coding points (20) with the largest elevation value are determined according to the magnitude sequence of the elevation value, the lower coding points (20) with the smallest elevation value are determined as the upper coding points (20) of the stratum, and the lower coding points are corrected according to the situation, wherein the method comprises the following steps of: acquiring a lenticular stratum code in the code point group A, searching code points (20) containing the lenticular stratum code in all drill holes (4), and marking the code points as a code point group B, wherein if only the lenticular stratum and the conventional stratum code are in the code point group B, the upper code points and the lower code points (20) of the conventional stratum are not corrected; if the coding point group B contains the adjacent regular stratum coding on the regular stratum, modifying the upper coding point (20) of the regular stratum into the coding point (20) with the lowest coding elevation value of the lens stratum in the drill hole (4); if the code point group B contains the codes of the adjacent regular strata under the regular strata, the lower code point (20) of the regular strata is modified to be the code point (20) with the highest code elevation value of the lens stratum in the drill hole (4).

6. The automated modeling method of a complex three-dimensional geological model based on BIM according to claim 1, wherein in step 4.3, the adjacent conventional strata is determined by the encoding pair-merging encoding pair group obtained in step 2, and the determining method is as follows: acquiring the selected regular stratum code (16) and traversing each code pair (21) of the final code pair-merging code pair group, and taking a stratum corresponding to Noup in the code pair (21) as an upper adjacent regular stratum if the selected regular stratum code (16) appears as Nodown in the code pair (21); if the regular stratum code (16) is selected to be Noup in the code pair (21), the stratum corresponding to Nodown in the code pair (21) is taken as the next adjacent regular stratum; the final encoding pair-merging encoding pair group is the encoding pair-merging encoding pair group which completes all merging operations through the step 2.

7. The automated modeling method of a complex three-dimensional geologic model based on BIM according to claim 1, wherein in step 5, the specific process of creating the lenticular stratum is as follows:

step 5.1: selecting a lens stratum code, traversing all drill holes (4) to obtain upper stratum demarcation points (2) and lower stratum demarcation points of the lens stratum; taking the JT outer contour line (23) in the step 3.2 as a boundary, respectively carrying out spatial interpolation on an upper stratum demarcation point (2) and a lower stratum demarcation point (2) of the lens stratum, and establishing upper stratum interfaces and lower stratum interfaces of the lens stratum with limited range, which are called upper limited lens stratum interfaces and lower limited lens stratum interfaces;

Step 5.2: taking the centroid of a closed polygon formed by JT outer contour lines (23) as a starting point, making a plurality of radial lines outwards, and simultaneously stretching the radial lines vertically up and down to form a plurality of vertical tangential planes (44); adopting the vertical tangential planes (44) to cut the upper and lower limited lens stratum interfaces, and obtaining upper and lower limited lens stratum interface curves in each vertical tangential plane (44); one end, coinciding with the JT outer contour line (23), of the vertical projection of the stratum interface curve of the lens body is marked as a transition end point (45);

if the whole lens stratum is positioned in a certain conventional stratum, solving tangent lines (41) of upper and lower limited lens stratum interface curves at transition end points (45) of the tangent lines in each vertical tangent plane (44), and acquiring an intersection point of the two tangent lines (41); constructing a triangle enclosing area (43) by using the tangent intersection point (42) and the two transition end points (45), searching a point with zero thickness of the lens body in the enclosing area (43), and constructing complete upper and lower lens body stratum interface curves in the vertical tangent plane (44);

if the lenticular stratum is positioned at the adjacent position of two conventional strata, generating an intersection line of the vertical tangent plane (44) and the upper and lower stratum interfaces to be trimmed of the two conventional strata at the position in each vertical tangent plane (44), generating a straight line through two transition end points (45), searching a point with zero lenticular thickness in an intersection line enclosing area (43) of the straight line and the upper and lower stratum interfaces to be trimmed, and establishing complete upper and lower lenticular stratum interface curves in the vertical tangent plane (44);

step 5.5: repeating steps 5.1-5.4 to generate all lens stratum;

8. The automated modeling method of a complex three-dimensional geological model based on BIM according to claim 1, wherein in step 6, the concrete process of creating the pinch-out formation is as follows:

step 6.1: selecting a to-be-built pinch-out stratum code (16), traversing all drill holes (4) to obtain upper stratum demarcation points (2) and lower stratum demarcation points of the pinch-out stratum; taking the JT outer contour line (23) in the step 3.2 as a boundary, respectively carrying out spatial interpolation on an upper stratum demarcation point (2) and a lower stratum demarcation point (2) of a to-be-built pinch-out stratum, and establishing upper stratum interfaces and lower stratum interfaces of the pinch-out stratum with limited range, which are called upper limited pinch-out stratum interfaces and lower limited pinch-out stratum interfaces;

Step 6.2: extracting discrete points according to the mesh set size at the position where the JT outer contour line (23) and the integral outer contour line (22) are not overlapped in the step 3, and taking the discrete points as vertical points to make a horizontal straight line perpendicular to the JT outer contour line (23); stretching the straight line up and down along the vertical direction to form a vertical tangent plane (44), and sectioning upper and lower limited pinch-out stratum interfaces by adopting the vertical tangent plane (44), wherein upper and lower limited pinch-out stratum interface curves can be obtained in each vertical tangent plane (44); one end of the vertical projection of the limited pinch-out stratum interface curve, which coincides with the JT outer contour line (23), is marked as a transition tip vanishing point (46);

step 6.3: generating intersection lines of upper and lower stratum interfaces of upper and lower adjacent conventional stratum to be trimmed at the position of the to-be-built pinch-out stratum in each vertical tangent plane (44), and generating complete upper and lower pinch-out stratum interface curves according to whether the upper and lower pinch-out stratum to be-built has a pinch-out stratum body which is closely adjacent to and built or not, wherein the curves are as follows:

if the built pinch-out stratum exists in the upper and lower adjacent strata of the pinch-out stratum to be built, an upper pinch-out stratum interface curve and a lower pinch-out stratum interface curve of the built pinch-out stratum are generated in a vertical tangent plane (44), if the built pinch-out stratum interface curve exceeds the range of a transition pinch-out point (46), one or two of the upper and lower limited stratum interface curves of the pinch-out stratum to be built are replaced by the built pinch-out stratum interface curve on the premise of keeping the transition pinch-out point (46) unchanged, and the pinch-out points are searched in a line formed by connecting the transition pinch-out points (46) and an intersection line enclosing region (43) of the upper and lower stratum interfaces to be trimmed, and meanwhile, a complete upper and lower pinch-out stratum interface curve is built; the method for replacing the upper and lower limited stratum interface curves of the stratum to be pinch-out by the established pinch-out stratum interface curves comprises the following steps: if the pinch-out stratum to be built is adjacent to the pinch-out stratum to be built and meets the substitution requirement, substituting an upper limited pinch-out stratum interface curve of the pinch-out stratum to be built by a lower pinch-out stratum interface curve of the pinch-out stratum to be built which is adjacent to the pinch-out stratum to be built; if the pinch-out stratum to be built is adjacent to the built pinch-out stratum, and the substitution requirement is met, substituting the lower limited pinch-out stratum interface curve of the pinch-out stratum to be built with the upper pinch-out stratum interface curve of the built pinch-out stratum which is adjacent to the lower pinch-out stratum;

If the stratum immediately adjacent to the upper stratum and the lower stratum is not an pinch-out stratum or is an pinch-out stratum but is not established, generating a straight line through two transition pinch-out points (46), searching the pinch-out points in a crossing line enclosing area (43) of the straight line and the upper stratum interface to be modified and the lower stratum interface to be modified, and establishing a complete upper pinch-out stratum interface curve and a complete lower pinch-out stratum interface curve;

9. The automated modeling method for a complex three-dimensional geological model based on BIM according to claim 1, wherein in step 7, the specific process of correcting the stratum interface to be corrected and establishing the stratum body is as follows:

Step 7.1: acquiring a group of stratum interfaces to be modified in the step 4.3 and generating all coding points (20) of the stratum interfaces to be modified, and judging the modification mode of the group of stratum interfaces to be modified according to the coding points (20), wherein the method comprises the following specific steps:

if lenticular layer coding occurs in the coding points (20), executing a trimming mode I;

if the pinch-out stratum coding occurs in the coding point (20), executing a second trimming mode;

if the lens body and the pinch-out stratum code simultaneously appear in the code points (20), executing a finishing mode III;

the first trimming mode is as follows: removing coding points (20) containing lenticular stratum codes, reserving only one of the coding points (20) with coincident space coordinates, and generating a conventional stratum interface by adopting a spatial interpolation algorithm;

the second trimming mode is as follows: according to the stratum pinch-out body established in the pinch-out stratum code searching step 6, extracting a pinch-out boundary of the stratum pinch-out body; extracting discrete points on the pinch-out boundary according to the set size of the grid to serve as pinch-out additionally-arranged stratum demarcation points; merging within each borehole (4) the coding points (20) containing the pinch-out formation codes as dummy formation demarcation points; reserving one of the coding points (20) with coincident space coordinates, combining the additional formation demarcation points with the dummy formation demarcation points, and generating a conventional formation interface through a space interpolation algorithm;

The third trimming mode is as follows: removing coding points (20) containing the lens body codes by adopting a first trimming mode, and establishing a conventional stratum interface according to a second trimming mode;

step 7.2: according to the conventional stratum interface, a conventional stratum body to be modified is established, wherein the establishment mode is formed by surrounding upper and lower adjacent conventional stratum interfaces or terrain surfaces (26);

10. The automated modeling method of a complex three-dimensional geological model based on BIM according to claim 9, wherein in the second modification mode of step 7.1, the step of extracting the pinch-out boundary of the stratum pinch-out body is performed according to the form of the stratum pinch-out body, specifically as follows: if the pinch-out stratum body does not exist in the stratum pinch-out bodies which are adjacent to each other from top to bottom, the position where the thickness of the pinch-out stratum body is zero is extracted as a pinch-out boundary; if the pinch-out stratum body exists in the upper and lower adjacent pinch-out stratum bodies, extracting the position where the sum of the thicknesses of all the adjacent pinch-out stratum bodies is zero as a pinch-out boundary; in the second trimming mode 7.1, the virtual stratum demarcation points are combined according to the form of stratum pinch-out bodies, and the method specifically comprises the following steps: if the pinch-out stratum body does not exist in the stratum pinch-out bodies which are immediately adjacent to each other, solving a space coordinate average value of the coding points (20) comprising the pinch-out stratum codes in each drilling hole (4) as a dummy stratum demarcation point; if the pinch-out formation exists in the immediate vicinity of the upper and lower pinch-out formations, a spatial coordinate average value of the coding points (20) comprising the pinch-out formation and the immediate vicinity of the upper and lower formations is solved in each borehole (4) as a dummy formation demarcation point.