CN112233230B - Three-dimensional model construction method and device for fault structure in cut geological section - Google Patents

Abstract

The invention discloses a three-dimensional model construction method and a device for a cut geological section interrupt layer structure, wherein the method comprises the following steps: (1) Loading a stratum surface layer and a stratum boundary layer of the map cut geological section, and extracting fault line information data; (2) Recoding the stratum code based on stratum code to generate a character string; (3) Judging an unexposed stratum estimated area based on the stratum element morphology; (4) Judging the repetition or deletion condition of the stratum based on the stratum sequence, and obtaining an inferred stratum; (5) calculating the deduced stratum thickness, and dividing fault lines; (6) An inferred formation line is generated, and the original formation-level-like elements are segmented. The invention can effectively improve the modeling quality of fault structures in the cut geological section.

Description

Three-dimensional model construction method and device for fault structure in cut geological section

Technical Field

The invention relates to the field of geographic information and geology, in particular to a three-dimensional model construction method and device for a cut geological section interruption layer structure.

Background

Faults refer to structural forms produced by breaking of rock formations or rock masses when they are subjected to forces exceeding their strength limits, and their integrity is compromised. Faults develop widely in the crust, and are one of the more important structural objects in geological topography research. The map cut geological section is a geological section which is formed by selecting a certain direction on the geological section, and drawing the geological section according to various geological and geographic elements according to a certain scale by a projection method. In geologic maps, the continuity and order of the formation on the left and right sides of the fault line is often disrupted, thereby creating a phenomenon of loss or duplication of the formation sequence. On the map cut geological section map, reasonable inference and modeling of the fault structure related unexposed stratum are carried out based on fault development characteristics and stratum distribution rules, which are necessary preconditions for accurately carrying out map cut geological section modeling and are requirements for carrying out three-dimensional modeling based on sequence map cut geological section. Therefore, the modeling method research of the cut geological section fault layer structure is developed, and the method has important research significance and use value.

Disclosure of Invention

The invention aims to: aiming at the problems existing in the prior art, the invention provides a three-dimensional model construction method and device for a fault structure in a cut geological section with high modeling quality.

The technical scheme is as follows: the three-dimensional model construction method of the fault structure in the cut geological section comprises the following steps:

(1) Acquiring a stratum code set FD, a stratum line set LN and a fault line set FN of a stratum according to a stratum surface layer and a stratum boundary layer of the initial map-cut geological section;

(2) Encoding each stratum code in the stratum code set FD into a single character, and generating a character string Str along the section line direction;

(3) Reading any fault line FN from the fault line set FN _a Judging an unexposed stratum estimated area based on the stratum morphology;

(4) Judging the repetition or deletion condition of the stratum based on the stratum sequence according to the character string Str and the unexposed stratum deducing area to obtain a deduced stratum set GN;

(5) Calculating and deducing the stratum thickness of the stratum according to the set GN, and constructing a stratum thickness set GH;

(6) Generating inferred formation lines according to the formation thickness set GH, and constructing an inferred formation line set GL;

(7) Circularly executing the steps (3) - (6) until all fault lines in the fault line set FN are traversed, so as to obtain an inferred stratum line set corresponding to all fault lines;

(8) The method comprises the steps of dividing an unexposed stratum inferred region based on all inferred stratum line sets, generating inferred stratum levels, and forming a three-dimensional model of a fault structure through stitching.

Further, the step (1) specifically includes:

(1-1) loading stratum layers of the map-cut geological section, reading stratum codes of all geological section strata, and storing stratum codes into a stratum code set FD= { f _c Ic=1, 2, …, FCN }, where f _c Representing a c-th formation code, the FCN representing the number of formation codes;

(1-2) loading stratum boundary line diagram layers of the map-cut geological section, obtaining all stratum lines in the stratum boundary lines, and sequentially storing the stratum lines into a stratum line set LN= { LN along the direction of the section line _b B=1, 2, …, FCN-1}, where ln _b Indicating the b-th formation line, and FCN-1 indicates the number of formation lines;

(1-3) extracting all fault lines from stratum boundary line patterns of the map-cut geological section according to boundary line type attribute of stratum boundary lines, and storing the fault lines into a fault line set FN= { FN _a |a=1, 2, …, FNA }, where fn _a Represents the a-th fault line, and FNA represents the number of fault lines.

Further, the step (2) specifically includes:

(2-1) constructing a stratum symbol encoding table, wherein the stratum symbol encoding table sequences all stratum codes according to the order from new stratum to old stratum, and provides a stratum code for each stratum code, the stratum code is a character, and the new stratum order and the old stratum order can be distinguished through the character;

(2-2) ordering all formation codes in the formation code set FD from left to right along the section line direction;

(2-3) encoding all stratum codes in the sequenced set FD according to a stratum symbol encoding table to generate an encoding character string;

(2-4) if two adjacent characters in the coded character string are equal, merging the two characters into one character to obtain a character string Str.

Further, the step (3) specifically includes:

(3-1) reading any fault line fn from the formation line set LN _a And obtain fault line fn _a Is adjacent to the formation line ln on the left side _l And right adjacent formation line ln _r ；

(3-2) acquiring fault lines fn respectively _a End point p _ta (x _ta ,y _ta )、p _wa (x _wa ,y _wa ) Left adjacent formation line ln _l End point p of (2) _tl (x _tl ,y _tl )、p _wl (x _wl ,y _wl ) And right adjacent formation line ln _r End point p of (2) _tr (x _tr ,y _tr )、p _wr (x _wr ,y _wr )；

(3-3) root calculation of fn _a And ln _l Is the head end point spacing DT of (1) _la And tail point spacing DW _la And fn _a And ln _r Is the head end point spacing DT of (1) _ra And tail point spacing DW _ra ；

(3-4) if DT _la <DW _la And DT (DT) _ra >DW _ra The space shape of the stratum at the left side of the fault line is shown to be narrow at the upper part and wide at the lower part, and the stratum at the left side of the fault line is marked as an unexposed stratum inferred area; if DT _ra <DW _ra And DT (DT) _la >DW _la And indicating that the space shape of the stratum on the right side of the fault line is narrow at the upper part and wide at the lower part, and marking the stratum on the right side of the fault line as an unexposed stratum estimation area.

Further, the step (4) specifically includes:

(4-1) at fault line fn _a Dividing a character string Str into a left character string leftStr and a right character string Str at the position;

(4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, setting the rightStr as a substring to be processed and the lefttstr as a reference substring; if the stratum on the right side of the fault line is an unexposed stratum inferred region, setting the leftStr as a substring to be processed and setting the lightStr as a reference substring;

(4-3) searching whether a public substring exists in the substring to be processed and the reference substring, if so, indicating that the stratum sequence is repeated, and executing the step (4-4); otherwise, indicating that the stratum sequence is missing, and executing the step (4-5);

(4-4) extracting the parts except the common substring in the substring to be processed, and storing the extracted parts in the inferred stratum set GN= { GN _v (e _v ,f _v ) V=1, 2, …, n }, where gn _v Representing the v-th inferred formation, e _v Representing the position of the v-th inferred formation, f _v A formation code representing a v-th inferred formation, n representing a number of inferred formations;

(4-5) if the stratum on the right side of the fault line is an unexposed stratum estimation area, extracting n continuous strata on the left side of the unexposed stratum estimation area, and storing the n continuous strata in an estimated stratum set GN; otherwise, n continuous strata on the right side of the estimated region of the unexposed strata are extracted and stored in the estimated stratum set GN.

Further, the step (4-4) specifically includes:

(4-4-1) searching for a starting position of the occurrence of the common substring in the substring to be processed;

(4-4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, extracting n continuous characters at the right side of the public substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in an inferred stratum set GN;

(4-4-3) if the stratum on the right side of the fault line is an unexposed stratum inferred region, extracting n continuous characters on the left side of the common substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in the inferred stratum set GN.

Further, the step (5) specifically includes:

(5-1) reading the formation lines corresponding to all inferred formations in the set GN from the inferred formation set GN in the formation line set LN, and storing the formation lines into a subset GP= { GP of formation lines _d D=1, 2,..n+1 }, where gp _d Represents the d-th formation line, and n+1 represents the number of formation lines;

(5-2) reading the head end point A of one formation line and the head and tail end points B and C of the adjacent formation line on the right side of the head end point A of the formation line in the GP according to the section direction from left to right;

(5-3) calculating to obtain the area S of the triangle ABC according to the sea-de formula;

(5-4) calculating the distance gh from the point A to the stratum line BC based on the following formula, and correspondingly estimating the thickness of the stratum by taking the gh as the current stratum line;

wherein BC represents the distance between endpoints B and C;

(5-5) cyclically performing steps (5-2) - (5-4) until all inferred formation thicknesses are found and stored in a formation thickness set gh= { GH _t T=1, 2,..n }, where gh _t Represents the thickness of the t th stratum, and n represents the number of strata.

Further, the step (6) specifically includes:

(6-1) calculating a fault line fn _a Is divided into fault lines fn in equal proportion according to the estimated stratum thickness set GH _a ；

(6-2) calculating a fault line fn according to the following formula _a The dividing point p _t (x _t ,y _t )；

Wherein x is _ta 、y _ta Is a fault line fn _a Is the initial point coordinate, gh _t Is a fault line fn _a Correspondingly deducing the stratum thickness;

(6-3) if the formation left of the fault line is an unexposed formation estimation region, the region is estimated based on the division point p _t (x _t ,y _t ) Along the formation line ln _l Parallel lines are arranged in the direction; otherwise based on the partition point p _t (x _t ,y _t ) Along the formation line ln _r The directions are parallel lines, wherein ln _l 、ln _r Respectively fault line fn _a Left side adjacent formation line and right side adjacent formation line;

(6-4) steps (6-2) - (6-3) are performed in a loop until parallel lines on all the division points are completed to be generated and stored in the inferred formation line set GL = { GL _u |u=1, 2, …, m }, where gl _u Represents the ith inferred formation line, and m represents the number of inferred formation lines.

Further, the step (8) specifically includes:

(8-1) dividing the unexposed formation estimation area based on the estimated formation line set GL to generate an estimated formation level;

(8-2) assigning a value to the newly generated inferred formation surface from left to right along the cross-sectional direction based on the formation code f in the inferred formation set GN;

(8-3) stitching all inferred formation faces together to form a three-dimensional model of the fault structure.

The three-dimensional model construction device for the fault structure in the cut geological section comprises a processor and a computer program which is stored in a memory and can run on the processor, wherein the processor realizes the method when executing the program.

The beneficial effects are that: compared with the prior art, the invention has the remarkable advantages that: the invention provides a three-dimensional model construction method and device for a fault layer structure of a cut geological section.

Drawings

FIG. 1 is cut geological profile data employed in the present embodiment;

FIG. 2 is a flow chart of a method for constructing a three-dimensional model of a fault in a cut geological section provided by the invention;

FIG. 3 is a cross-sectional view of an unexposed formation evaluation area after formation lines have been generated;

FIG. 4 is a diagrammatic cross-sectional view of a fault construct modeled in accordance with the present invention.

Detailed Description

In the following, the technical scheme of the present invention is further described in detail, in this embodiment, a map-cut geological section of the mountain of sauvignon in south kyo city is selected as experimental data, and as shown in fig. 1, a projection coordinate system adopted by the experimental data is WGS84. Further description will be provided by describing a specific embodiment with reference to the accompanying drawings.

As shown in fig. 2, the three-dimensional model construction method for a fault structure in a cut geological section provided in this embodiment includes:

(1) And acquiring a stratum code set FD, a stratum line set LN and a fault line set FN of the stratum according to the stratum surface layer and the stratum boundary layer of the initial map cut geological section.

The method specifically comprises the following steps:

(1-1) loading a layer map layer of a map cut geological section, reading all geological sectionsStratum codes of stratum are stored in stratum code set FD= { f _c Ic=1, 2, …, FCN }, where f _c Representing a c-th formation code, fc representing the number of formation codes; in the present embodiment, fcn=9;

(1-2) loading stratum boundary line diagram layers of the map-cut geological section, obtaining all stratum lines in the stratum boundary lines, and sequentially storing the stratum lines into a stratum line set LN= { LN along the direction of the section line _b B=1, 2, …, FCN-1}, where ln _b Indicating the b-th formation line, and FCN-1 indicates the number of formation lines;

(1-3) extracting all fault lines from stratum boundary line patterns of the map-cut geological section according to boundary line type attribute of stratum boundary lines, and storing the fault lines into a fault line set FN= { FN _a |a=1, 2, …, FNA }, where fn _a Represents the a-th fault line, fna represents the number of fault lines. In this embodiment, fna=1, and the 5 th formation line is a fault line.

(2) Each formation code in the formation code set FD is encoded as a single character and a string Str is generated in the cross-hatching direction.

The method specifically comprises the following steps:

(2-1) constructing a stratum symbol encoding table, wherein the stratum symbol encoding table sequences all stratum codes according to the order from new stratum to old stratum, and provides a stratum code for each stratum code, the stratum code is a character, and the new stratum order and the old stratum order can be distinguished through the character; the stratigraphic symbol encoding table is shown in table 1.

Table 1 stratigraphic symbol encoding table

The stratum is sequentially from new to old: q (Q) ₄ 、Q ₃ 、Q ₂ 、...、K ₁ s ⁴ 、K ₁ s ³ 、K ₁ s ² 、...、D ₃ w ¹ 、S ₃ m、S ₂ f ² 、...、Z ₁ l ³ 、Z ₁ l ² 、Z ₁ l ¹ The stratum isThe codes are the stratum codes used in this embodiment, and in other embodiments, the stratum codes may be ordered according to the stratum codes used, and the stratum codes are sequentially given with 0 to 9, a to Z, and a to Z from new to old, and in other embodiments, may be single characters in other orders.

(2-2) ordering all formation codes in the formation code set FD from left to right along the section line direction.

(2-3) encoding all the formation codes in the sorted set FD according to the formation symbol encoding table to generate an encoded string.

(2-4) if two adjacent characters in the coded character string are equal, merging the two characters into one character to obtain a character string Str. The final formation code encoding results in this example are shown in table 2. It can be seen that str=0 fgcavu in this example.

Table 2 formation code encoding

Formation code

Q 4

J 3 l

J 3 x

S 2 f 2

D 3 w 1

S 2 f 2

Formation coding

0

F

G

c

a

c

Formation code

D 3 w 1

C 3 c

P 1 q

Formation coding

a

V

U

(3) Reading any fault line FN from the fault line set FN _a And based on the formation morphology, judging the unexposed formation estimated area.

The method specifically comprises the following steps:

(3-1) reading any fault line fn from the formation line set LN _a And obtain fault line fn _a Is adjacent to the formation line ln on the left side _l And right adjacent formation line ln _r ；

(3-2) acquiring fault lines fn respectively _a End point p _ta (x _ta ,y _ta )、p _wa (x _wa ,y _wa ) Left adjacent formation line ln _l End point p of (2) _tl (x _tl ,y _tl )、p _wl (x _wl ,y _wl ) And right adjacent formation line ln _r End point p of (2) _tr (x _tr ,y _tr )、p _wr (x _wr ,y _wr )；

(3-3) root calculation of fn _a And ln _l Is the head end point spacing DT of (1) _la And tail point spacing DW _la And fn _a And ln _r Is the head end point spacing DT of (1) _ra And tail point spacing DW _ra ；

(3-4) if DT _la <DW _la And DT (DT) _ra >DW _ra The space shape of the stratum at the left side of the fault line is shown to be narrow at the upper part and wide at the lower part, and the stratum at the left side of the fault line is marked as an unexposed stratum inferred area; if DT _ra <DW _ra And DT (DT) _la >DW _la And indicating that the space shape of the stratum on the right side of the fault line is narrow at the upper part and wide at the lower part, and marking the stratum on the right side of the fault line as an unexposed stratum estimation area.

(4) And judging the repetition or deletion condition of the stratum based on the stratum sequence according to the character string Str and the unexposed stratum deducing area, and obtaining a deduced stratum set GN.

The method specifically comprises the following steps:

(4-1) at fault line fn _a Dividing a character string Str into a left character string leftStr and a right character string Str at the position; in this embodiment, leftstr=0fgca, lightstr=cavu;

(4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, setting the rightStr as a substring to be processed and the lefttstr as a reference substring; if the stratum on the right side of the fault line is an unexposed stratum inferred region, setting the leftStr as a substring to be processed and setting the lightStr as a reference substring; in this embodiment, leftStr is the substring to be processed, and lightstr is the reference substring;

(4-3) searching whether a public substring exists in the substring to be processed and the reference substring, if so, indicating that the stratum sequence is repeated, and executing the step (4-4); otherwise, indicating that the stratum sequence is missing, and executing the step (4-5);

(4-4) extracting the parts except the common substring in the substring to be processed, and storing the extracted parts in the inferred stratum set GN= { GN _v (e _v ,f _v ) V=1, 2, …, n }, where gn _v Representing the v-th inferred formation, e _v Representing the position of the v-th inferred formation, f _v A formation code representing a v-th inferred formation, n representing a number of inferred formations;

(4-5) if the stratum on the right side of the fault line is an unexposed stratum estimation area, extracting n continuous strata on the left side of the unexposed stratum estimation area, and storing the n continuous strata in an estimated stratum set GN; otherwise, n continuous strata on the right side of the estimated region of the unexposed strata are extracted and stored in the estimated stratum set GN.

Wherein, the step (4-4) specifically comprises:

(4-4-1) searching for a starting position of the occurrence of the common substring in the substring to be processed;

(4-4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, extracting n continuous characters at the right side of the public substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in an inferred stratum set GN;

(4-4-3) if the stratum on the right side of the fault line is an unexposed stratum inferred region, extracting n continuous characters on the left side of the common substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in the inferred stratum set GN. In this example, the inferred formation set is shown in Table 3.

TABLE 3 inference of stratigraphic sets

k	1	2	3
				e	1	2	3
f	Q 4	J 3 l	J 3 x

(5) And calculating and deducing the stratum thickness of the stratum according to the set GN, and constructing a stratum thickness set GH.

The method specifically comprises the following steps:

(5-1) reading the formation lines corresponding to all inferred formations in the set GN from the inferred formation set GN in the formation line set LN, and storing the formation lines into a subset GP= { GP of formation lines _d D=1, 2,..n+1 }, where gp _d Represents the d-th formation line, and n+1 represents the number of formation lines;

(5-2) reading the head end point A of one formation line and the head and tail end points B and C of the adjacent formation line on the right side of the head end point A of the formation line in the GP according to the section direction from left to right;

(5-3) calculating to obtain the area S of the triangle ABC according to the sea-de formula;

(5-4) calculating the distance gh from the point A to the stratum line BC based on the following formula, and correspondingly estimating the thickness of the stratum by taking the gh as the current stratum line;

wherein BC represents the distance between endpoints B and C;

(5-5) cyclically performing steps (5-2) - (5-4) until all inferred formation thicknesses are found and stored in a formation thickness set gh= { GH _t T=1, 2,..n }, where gh _t Represents the thickness of the t th stratum, and n represents the number of strata.

(6) And generating inferred formation lines according to the formation thickness set GH, and constructing an inferred formation line set GL.

The method specifically comprises the following steps:

(6-1) calculating a fault line fn _a Is divided into fault lines fn in equal proportion according to the estimated stratum thickness set GH _a ；

(6-2) calculating a fault line fn according to the following formula _a The dividing point p _t (x _t ,y _t )；

Wherein x is _ta 、y _ta Is a fault line fn _a Is the initial point coordinate, gh _t Is a fault line fn _a Correspondingly deducing the stratum thickness;

(6-3) if the formation left of the fault line is an unexposed formation estimation region, the region is estimated based on the division point p _t (x _t ,y _t ) Along the formation line ln _l Parallel lines are arranged in the direction; otherwise based on the partition point p _t (x _t ,y _t ) Along the formation line ln _r The directions are parallel lines, wherein ln _l 、ln _r Respectively fault line fn _a Left side adjacent formation line and right side adjacent formation line;

(6-4) steps (6-2) - (6-3) are performed in a loop until parallel lines on all the division points are completed to be generated and stored in the inferred formation line set GL = { GL _u |u=1, 2, …, m }, where gl _u Represents the ith inferred formation line, and m represents the number of inferred formation lines. In this embodiment, since the formation thickness is inferredThe degree exceeds the fault line length, so the number of formation lines m=2 is inferred.

(7) And (3) to (6) are circularly executed until all fault lines in the fault line set FN are traversed, and an inferred stratum line set corresponding to all fault lines is obtained.

(8) The method comprises the steps of dividing an unexposed stratum inferred region based on all inferred stratum line sets, generating inferred stratum levels, and forming a three-dimensional model of a fault structure through stitching.

The method specifically comprises the following steps:

(8-1) dividing the unexposed formation estimation area based on the estimated formation line set GL to generate an estimated formation level; as shown in fig. 3;

(8-2) assigning a value to the newly generated inferred formation surface from left to right along the cross-sectional direction based on the formation code f in the inferred formation set GN;

(8-3) stitching all inferred formation faces together to form a three-dimensional model of the fault structure. In this embodiment, a fault formation in a constructed cut geological section is shown in FIG. 4.

The embodiment also provides a three-dimensional model construction device for the cut geological section interruption layer structure, which comprises a processor and a computer program stored on a memory and capable of running on the processor, wherein the processor realizes the method when executing the program.

The above disclosure is only a preferred embodiment of the present invention and should not be construed as limiting the scope of the invention, which is defined by the appended claims.

Claims (8)

1. A three-dimensional model construction method of a cut geological section interruption layer structure is characterized by comprising the following steps:

(1) Acquiring a stratum code set FD, a stratum line set LN and a fault line set FN of a stratum according to a stratum surface layer and a stratum boundary layer of the initial map-cut geological section;

(2) Encoding each stratum code in the stratum code set FD into a single character, and generating a character string Str along the section line direction;

(3) Reading any fault line FN from the fault line set FN _a Judging an unexposed stratum estimated area based on the stratum morphology;

(4) Judging the repetition or deletion condition of the stratum based on the stratum sequence according to the character string Str and the unexposed stratum deducing area to obtain a deduced stratum set GN;

(5) Calculating and deducing the stratum thickness of the stratum according to the set GN, and constructing a stratum thickness set GH;

(6) Generating inferred formation lines according to the formation thickness set GH, and constructing an inferred formation line set GL;

(7) Circularly executing the steps (3) - (6) until all fault lines in the fault line set FN are traversed, so as to obtain an inferred stratum line set corresponding to all fault lines;

(8) Dividing the unexposed stratum inferred region based on all inferred stratum line sets, generating an inferred stratum level, and forming a three-dimensional model of a fault structure through stitching;

the step (5) specifically comprises:

(5-1) reading the formation lines corresponding to all inferred formations in the set GN from the inferred formation set GN in the formation line set LN, and storing the formation lines into a subset GP= { GP of formation lines _d D=1, 2,..n+1 }, where gp _d Represents the d-th formation line, and n+1 represents the number of formation lines;

(5-2) reading the head end point A of one formation line and the head and tail end points B and C of the adjacent formation line on the right side of the head end point A of the formation line in the GP according to the section direction from left to right;

(5-3) calculating to obtain the area S of the triangle ABC according to the sea-de formula;

(5-4) calculating the distance gh from the point A to the stratum line BC based on the following formula, and correspondingly estimating the thickness of the stratum by taking the gh as the current stratum line;

wherein BC represents the distance between endpoints B and C;

(5-5) cyclically performing steps (5-2) - (5-4) until all inferred formation thicknesses are found and the formation thickness is storedDegree set gh= { GH _t T=1, 2,..n }, where gh _t Represents the thickness of the t stratum, and n represents the number of stratum;

the step (6) specifically comprises:

(6-1) calculating a fault line fn _a Is divided into fault lines fn in equal proportion according to the estimated stratum thickness set GH _a ；

(6-2) calculating a fault line fn according to the following formula _a The dividing point p _t (x _t ,y _t )；

Wherein x is _ta 、y _ta Is a fault line fn _a Is the initial point coordinate, gh _t Is a fault line fn _a Correspondingly deducing the stratum thickness;

(6-3) if the formation left of the fault line is an unexposed formation estimation region, the region is estimated based on the division point p _t (x _t ,y _t ) Along the formation line ln _l Parallel lines are arranged in the direction; otherwise based on the partition point p _t (x _t ,y _t ) Along the formation line ln _r The directions are parallel lines, wherein ln _l 、ln _r Respectively fault line fn _a Left side adjacent formation line and right side adjacent formation line;

(6-4) steps (6-2) - (6-3) are performed in a loop until parallel lines on all the division points are completed to be generated and stored in the inferred formation line set GL = { GL _u |u=1, 2, …, m }, where gl _u Represents the ith inferred formation line, and m represents the number of inferred formation lines.

2. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (1) specifically comprises:

(1-1) Loading cut geologyStratum layers of stratum of the section, stratum codes of all geological section strata are read, and stratum codes are stored into a stratum code set FD= { f _c Ic=1, 2, …, FCN }, where f _c Representing a c-th formation code, the FCN representing the number of formation codes;

(1-2) loading stratum boundary line diagram layers of the map-cut geological section, obtaining all stratum lines in the stratum boundary lines, and sequentially storing the stratum lines into a stratum line set LN= { LN along the direction of the section line _b B=1, 2, …, FCN-1}, where ln _b Indicating the b-th formation line, and FCN-1 indicates the number of formation lines;

(1-3) extracting all fault lines from stratum boundary line patterns of the map-cut geological section according to boundary line type attribute of stratum boundary lines, and storing the fault lines into a fault line set FN= { FN _a |a=1, 2, …, FNA }, where fn _a Represents the a-th fault line, and FNA represents the number of fault lines.

3. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (2) specifically comprises:

(2-1) constructing a stratigraphic symbol encoding table which orders all stratigraphic codes in the order of new stratigraphic codes from old stratigraphic codes, and providing a stratigraphic code for each stratigraphic code, wherein the stratigraphic code is a character, and the new stratigraphic codes and the old stratigraphic codes can be distinguished through the character;

(2-2) ordering all formation codes in the formation code set FD from left to right along the section line direction;

(2-3) encoding all stratum codes in the sequenced set FD according to a stratum symbol encoding table to generate an encoding character string;

(2-4) if two adjacent characters in the coded character string are equal, merging the two characters into one character to obtain a character string Str.

4. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (3) specifically comprises:

(3-1) reading any fault line fn from the formation line set LN _a And obtain fault line fn _a Is adjacent to the formation line ln on the left side _l And right adjacent formation line ln _r ；

(3-2) acquiring fault lines fn respectively _a End point p _ta (x _ta ,y _ta )、p _wa (x _wa ,y _wa ) Left adjacent formation line ln _l End point p of (2) _tl (x _tl ,y _tl )、p _wl (x _wl ,y _wl ) And right adjacent formation line ln _r End point p of (2) _tr (x _tr ,y _tr )、p _wr (x _wr ,y _wr )；

(3-3) root calculation of fn _a And ln _l Is the head end point spacing DT of (1) _la And tail point spacing DW _la And fn _a And ln _r Is the head end point spacing DT of (1) _ra And tail point spacing DW _ra ；

(3-4) if DT _la <DW _la And DT (DT) _ra >DW _ra The space shape of the stratum at the left side of the fault line is shown to be narrow at the upper part and wide at the lower part, and the stratum at the left side of the fault line is marked as an unexposed stratum inferred area; if DT _ra <DW _ra And DT (DT) _la >DW _la And indicating that the space shape of the stratum on the right side of the fault line is narrow at the upper part and wide at the lower part, and marking the stratum on the right side of the fault line as an unexposed stratum estimation area.

5. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (4) specifically comprises:

(4-1) at fault line fn _a Dividing a character string Str into a left character string leftStr and a right character string Str at the position;

(4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, setting the rightStr as a substring to be processed and the lefttstr as a reference substring; if the stratum on the right side of the fault line is an unexposed stratum inferred region, setting the leftStr as a substring to be processed and setting the lightStr as a reference substring;

(4-3) searching whether a public substring exists in the substring to be processed and the reference substring, if so, indicating that the stratum sequence is repeated, and executing the step (4-4); otherwise, indicating that the stratum sequence is missing, and executing the step (4-5);

(4-4) extracting the parts except the common substring in the substring to be processed, and storing the extracted parts in the inferred stratum set GN= { GN _v (e _v ,f _v ) V=1, 2, …, n }, where gn _v Representing the v-th inferred formation, e _v Representing the position of the v-th inferred formation, f _v A formation code representing a v-th inferred formation, n representing a number of inferred formations;

(4-5) if the stratum on the right side of the fault line is an unexposed stratum estimation area, extracting n continuous strata on the left side of the unexposed stratum estimation area, and storing the n continuous strata in an estimated stratum set GN; otherwise, n continuous strata on the right side of the estimated region of the unexposed strata are extracted and stored in the estimated stratum set GN.

6. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (4-4) specifically comprises:

(4-4-1) searching for a starting position of the occurrence of the common substring in the substring to be processed;

(4-4-2) if the stratum at the left side of the fault line is an unexposed stratum inferred region, extracting n continuous characters at the right side of the public substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in an inferred stratum set GN;

(4-4-3) if the stratum on the right side of the fault line is an unexposed stratum inferred region, extracting n continuous characters on the left side of the common substring, searching the corresponding position and stratum code in the Str, and storing the position and stratum code in the inferred stratum set GN.

7. The method for constructing a three-dimensional model of a fault structure in a cut-to-map geological section according to claim 1, wherein: the step (8) specifically comprises:

(8-1) dividing the unexposed formation estimation area based on the estimated formation line set GL to generate an estimated formation level;

(8-2) assigning a value to the newly generated inferred formation surface from left to right along the cross-sectional direction based on the formation code f in the inferred formation set GN;

(8-3) stitching all inferred formation faces together to form a three-dimensional model of the fault structure.

8. A three-dimensional model construction device for a slice geological profile interrupt layer structure, comprising a processor and a computer program stored on a memory and executable on the processor, characterized in that: the processor, when executing the program, implements the method of any one of claims 1-7.