Method and device for automatically processing the boundary of strata in adjacent geological sections

Abstract

The present invention discloses a method and device for automatically processing the consistency of stratum boundaries in adjacent geological profiles, the method comprising: (1) obtaining adjacent geological profiles A and profile B; (2) extracting a set of stratum boundary lines in the two profiles; (3) creating a merged set of stratum boundary lines formed by any combination of all strata that can be connected in the profiles; (4) comparing the stratum boundary lines created by the same stratum combination in the merged set of stratum boundary lines based on the stratum boundary line set, and finding a unique topologically inconsistent stratum boundary line; (5) creating a virtual stratum for all unique topologically inconsistent stratum boundary lines; (6) adding the virtual stratum to profiles A and B, respectively, to generate adjacent geological profiles with completely consistent stratum boundary line topological relationships. The present invention has high efficiency and high quality.

Description

Technical Field

[0001] The invention relates to the field of geographic information and geology, and in particular to a method and device for processing stratum boundary consistency in geological profile data. Background Art

[0002] 3D geological modeling is an important combination of geographic information science and geology. Through 3D geological modeling, the underground geological structure can be vividly and intuitively reproduced. The basic raw data of 3D geological modeling are mainly drilling data, geological profiles, etc. Adjacent geological profiles are one type of geological profiles. Adjacent geological profiles refer to geological profiles that are adjacent and close in the study area. It is necessary to make full use of the existing adjacent geological profiles for 3D geological modeling.

[0003] In the field of geological profile 3D modeling, the mainstream 3D geological profile modeling method is to establish geological models by using interpolation surface construction and connecting geological profiles (Li Chen. Research on 3D complex geological body modeling technology based on profiles [D]. Chengdu University of Technology, 2018). The interpolation surface construction method is mainly suitable for geological structures with clear hierarchical structures, while the method of connecting geological profiles can be adapted to more common adjacent geological profile data, but this method requires that the same stratum boundary lines between profiles can be well matched. However, due to geological effects such as faults, folds, and intrusions, the contact relationship between strata in adjacent geological profiles has changed significantly. In most cases, the stratum boundary lines in adjacent geological profiles are difficult to match.

[0004] At present, the problem of corresponding stratum boundaries in adjacent geological sections is mainly solved by relying on expert experience to find the corresponding boundary lines through manual interpretation, and artificially constructing transition sections during the change of contact relationship to assist modeling (Qu Honggang, Pan Mao, Wang Yong, et al. Three-dimensional geological modeling based on topological sections [J]. Journal of Peking University: Natural Science Edition, 2006 (06): 24-30). This auxiliary modeling method is inefficient and the quality of the results varies from person to person. Summary of the invention

[0005] Purpose of the invention: In view of the problems existing in the prior art, the present invention provides a highly efficient method and device for processing the consistency of stratum boundaries in geological profile data.

[0006] Technical solution: The method for automatically processing the boundaries of strata in adjacent geological sections according to the present invention comprises:

[0007] (1) Load two adjacent geological sections to be processed, record them as Section A and Section B, obtain the stratigraphic code set FD of all strata, the stratigraphic set SectionA of Section A, and the stratigraphic set SectionB of Section B;

[0008] (2) Extract the boundary lines and boundary line endpoints of Section A and Section B respectively, and store them in boundary line sets ArcsA and ArcsB and boundary line endpoint sets NodesA and NodesB;

[0009] (3) According to the geometric attributes of SectionA, ArcsA, and NodesA, the topological contact relationship table FA1 between the boundary line and the stratum, the topological contact relationship table FA2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FA3 between the strata in Section A are established; for SectionB, ArcsB, and NodesB, the topological contact relationship table FB1 between the boundary line and the stratum, the topological contact relationship table FB2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FB3 between the strata in Section B are obtained by the same processing;

[0010] (4) Generate all combinations of connected strata in profile A and profile B according to adjacency tables FA3 and FB3, and merge the connected strata of each combination respectively. The boundary lines of the merged strata are stored in the merged strata boundary line sets La and Lb respectively;

[0011] (5) Take out any two corresponding merged boundary lines in La and Lb, and check whether the stratum boundary lines before merging correspond to each other according to FA1, FA2, FB1, and FB2, and obtain the set La after the inconsistent lines are grouped. block , Lb block , the corresponding relationship is stored in the relationship table Pair Block ;

[0012] (6) For La block , Lb block For each element in, a virtual stratum is generated and stored in the virtual stratum set SA virtual , SB virtual ;

[0013] (7) Set the virtual stratum SA virtual , SB virtual The virtual stratum in the relationship table Pair Block Renumber and merge with SectionA and SectionB to form Section A with consistent topological relationship fit , Section B fit , Section A fit , Section B fit The results are after consistency processing.

[0014] Furthermore, step (1) specifically includes:

[0015] (1-1) The two adjacent geological sections to be processed are recorded as Section A and Section B, the vector layers of Section A and Section B are read, and the stratigraphic codes of all stratigraphic layers are read from the vector layers and stored in the stratigraphic code set FD = {f i |i=1,2,…,fn}, where f i represents the i-th stratum code, and fn represents the number of stratum codes;

[0016] (1-2) Read all the stratigraphic data in section A, renumber the stratigraphic data, and store them in the stratigraphic set Section A = {StraA a |a=1,2,…,San}, where StraA a represents the ath stratum, and San represents the number of strata in SectionA;

[0017] (1-3) Read all the stratigraphic data in section B, renumber the stratigraphic data, and store them in the stratigraphic set Section B = {StraB b |b=1,2,…,Sbn}, where StraB b represents the b-th stratum, and Sbn represents the number of strata in SectionB.

[0018] Furthermore, step (2) specifically includes:

[0019] (2-1) Read any stratum StraA from the stratum set SectionA i ;

[0020] (2-2) Read all other strata StraA from the stratum set SectionA j , j=1,2,...,San, and j≠i; and analyze StraA i With StraA j If two strata are in contact, the common boundary line arc of the two strata is obtained. t ;

[0021] (2-3) Check arc t Does it exist in the boundary line set ArcsA? If not, set arc t Store in the boundary line set ArcsA;

[0022] (2-4) Steps (2-1) to (2-3) are executed repeatedly until all strata in SectionA are traversed and a complete boundary line set ArcsA is obtained;

[0023] (2-5) Read an arbitrary boundary line arc from the boundary line set ArcsA t , and extract arct The first and last endpoints p 1 、p 2 ;

[0024] (2-6) Check p 1 、p 2 Whether it exists in the boundary endpoint set NodesA, if not, then store it in the boundary endpoint set NodesA;

[0025] (2-7) Loop through steps (2-5) to (2-6) until ArcsA is completely traversed and all boundary endpoints are stored in NodesA;

[0026] (2-8) For set SectionB, perform the same processing as steps (2-1)-(2-7) to obtain the boundary line set ArcsB and the boundary line endpoint set NodesB.

[0027] Furthermore, step (3) specifically includes:

[0028] (3-1) Read any boundary line arc in ArcsA t ;

[0029] (3-2) Arc t Check the topological contact relationship with all the strata in SectionA. If arc t StraA a Contact, then StartA a The number a is stored in the arc in the topological contact relationship table FA1 t The corresponding line;

[0030] (3-3) looping through steps (3-1) to (3-2) until all boundary lines in ArcsA are traversed, and obtaining a topological contact relationship table FA1 that records the contact stratum numbers of all boundary lines;

[0031] (3-4) Read any point p in NodesA k ;

[0032] (3-5) for p k Perform a topological check with all boundary lines in the set ArcsA. If p k On a certain boundary arc t On, the arc t The number t is stored in the point p in the topological contact relationship table FA2 k The corresponding line;

[0033] (3-6) Loop through steps (3-4) to (3-5) until all points in NodesA are traversed, and obtain a topological contact relationship table FA2 that records the contact boundary line numbers of all boundary line endpoints;

[0034] (3-7) Take any stratum StraA from the set SectionA a ;

[0035] (3-8) Take all the remaining strata StraA from the set SectionA a’ , where a≠a'; for StraA a With StraA a’ Perform contact topology check. If there is a contact relationship, (a, a') is stored in the adjacency relationship table FA3 between layers.

[0036] (3-9) Loop through steps (3-7) to (3-8) until all strata in Section A are traversed, and obtain the adjacency table FA3 of all strata in Section A;

[0037] (3-10) For SectionB, ArcsB and NodesB, the same processing is performed according to steps (3-1)-(3-7) to obtain the topological contact relationship table FB1, the topological contact relationship table FB2, and the adjacency relationship table FB3.

[0038] Furthermore, step (4) specifically includes:

[0039] (4-1) Store all the stratum numbers in SectionA and SectionB into the stratum number set Ser = {ser n |n=1,2,…,N}, where ser n is the number of the nth stratum, and N is the total number of strata;

[0040] (4-2) For the stratum number set Ser, generate all combinations of any stratum numbers and store them in the set Comb = {comb m |m=1,2,…,M}, where comb m is the mth stratigraphic number combination, M is the total number of sets;

[0041] (4-3) Read a combination comb in the set Comb m ;

[0042] (4-4) According to the adjacency table FA3, get comb m All strata in comb m The adjacent situation of other strata in the stratum is stored in the set linka m={(u,v)|u,v∈comb m , (u,v)∈FA3}, where u and v are the numbers of adjacent strata respectively, and (u,v) means that strata u and v are adjacent;

[0043] (4-5) Take the strata in sectionA as nodes and linka m The adjacent relationships in are used as edges to construct and query the unionSet a , and linka m All relations in the unionSet are added to the unionSet a ;

[0044] (4-6) For the stratum adjacency table FB3, follow the same process as steps (4-4)-(4-5) to obtain the set linkb m ={(u,v)|u,v∈comb m , (u,v)∈FB3} and unionSet b ;

[0045] (4-7) Check comb m Are all the stratum numbers in the unionSet a belong to the same parent node and are in the unionSet b belong to the same parent node. If they all belong to the same parent node, then comb m All strata corresponding to SectionA are merged into a new stratum StraA m , comb m All the strata corresponding to SectionB are merged into a new stratum StraB m ;

[0046] (4-8) Extraction of StraA m The merged boundary line ea is stored in the merged boundary line set La = {ea m |m≤M}, where ea m Indicates the merged boundary line numbered m, where m is the number of comb m Corresponding combination number in Comb; extract StraB m The merged boundary line eb is stored in the merged boundary line set Lb = {eb m |m≤M}, where eb m Indicates the merged boundary line numbered m;

[0047] (4-9) Loop through steps (4-3) to (4-8) until all combinations in Comb are traversed and a complete set of merged boundary lines La, Lb is generated.

[0048] Furthermore, step (5) specifically includes:

[0049] (5-1) Read any merged boundary line ea in La m , and the merged boundary line eb with the same number in Lb m ;

[0050] (5-2) Traverse all boundary lines in ArcsA and compare them with ea m Completely overlapping boundary lines are stored in the collection arcs ea ; Traverse all the boundary lines in ArcsB and compare them with eb m All completely overlapping boundary lines are stored in the collection arcs eb ;

[0051] (5-3) Set arcs ea All boundary lines in the set arcs are reordered according to the connection relationship in the topological contact relationship table FA2 so that the boundary lines with adjacent numbers are connected; eb All boundary lines in are reordered according to the connection relationship in the topological contact relationship table FB2 so that the boundary lines with adjacent numbers are connected;

[0052] (5-4) Comparison set arcs ea With collection arcs eb The boundary lines in FA1 and FB1 are matched according to the contact relationship between the boundary lines and the strata. If a unique set of topologically inconsistent and continuous boundary lines can be matched, they are stored in the set la. m lb m , and m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship (a m ,b m ) is stored in the relationship set Pair Block ;

[0053] (5-5) Loop through steps (5-1) to (5-4) until all merged boundary lines in sets La and Lb are traversed to obtain a complete set of corresponding relationships Pair Block .

[0054] Further, step (5-4) specifically includes:

[0055] (5-4-1) Traversing arcs ea middle boundary line until a boundary line l is obtained t1, whose contact stratigraphic number is the same as arcs eb A boundary line with a unique number in a contact stratum t2 same;

[0056] (5-4-2) Collect arcs ea The order of the elements in is adjusted to l t1 As the first element; set arcs eb The order of the elements in is adjusted to l t2 is the first element;

[0057] (5-4-3) From the collection arcs ea Read an element l from front to back ea , according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs eb Read an element l from front to back eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ;

[0058] (5-4-4) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2 ) are the same, if they are the same, return to step (5-4-3), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa1, pb1 in;

[0059] (5-4-5) Check whether arcs have been traversed ea or arcs eb If the traversal is not completed, continue to execute step (5-4-6); if the traversal is completed, it means that the actual boundary lines on the merged boundary line can correspond one to one, or there is only a topological inconsistency at the end and the corresponding position at the end of the other merged boundary line is empty. At the same time, it means that the line does not contain a unique continuous boundary line with inconsistent topology, and no new elements are stored in the set La block , set Lb block , set Pair Block , execute step (5-5);

[0060] (5-4-6) From the collection arcs ea Read an element l from back to front ea, according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs eb Read an element l from back to front eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ;

[0061] (5-4-7) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2 ) are the same, if they are the same, return to step (5-4-6), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa2, pb2 in;

[0062] (5-4-8) If pa2-pa1=0 or pb2-pb1=0, it means that the line does not contain a unique topologically inconsistent continuous boundary line, and no new elements are stored in the set La. block , set Lb block , set Pair Block , execute step (5-5); otherwise, continue to execute step (5-4-9);

[0063] (5-4-9) If pa2-pa1>1 and pb2-pb1>1, then traverse arcs ea The elements in the interval [pa1+1,pa2-1] can be eb If the two contacting strata are found to have exactly the same elements in the interval [pb1+1, pb2-1], it means that the interval [pb1+1, pb2-1] contains multiple sets of topologically inconsistent continuous boundary lines. The positions of the elements of the contacting strata found to have exactly the same elements in profile A and profile B are recorded as pat and pbt, and step (5-4-10) is executed; otherwise, step (5-4-11) is executed;

[0064] (5-4-10) will gather arcs ea Divide into sets arcs based on position pat ea1 , arcs ea2 The left and right sets will set arcs eb Divide into sets arcs based on position pbt eb1 , arcs eb2 The left and right sets are respectively for arcs ea1、arcs eb1 , arcs ea2 、arcs eb2 Execute step (5-4) to obtain the topological inconsistency positions in the left and right sets;

[0065] (5-4-11) Arcs ea The elements in the interval [pa1+1,pa2-1] are stored in the set la m ; Arcs eb The elements in the interval [pb1+1,pb2-1] are stored in the set lb m ;

[0066] (5-4-12) m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship is stored in the set Pair Block .

[0067] Furthermore, step (6) specifically includes:

[0068] (6-1) Read La block Any set in la m ;

[0069] (6-2) Create a set la m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SA virtual ;

[0070] (6-3) Read Lb block Any set lb m ;

[0071] (6-4) Create a set lb m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SB virtual .

[0072] Furthermore, step (7) specifically includes:

[0073] (7-1) The strata in Section A and SA virtual The overlapping parts of the virtual strata are erased, and the remaining parts are stored in the collection SectionA new ;

[0074] (7-2) The strata in Section B and SB virtual The overlapping part of the virtual stratum is erased, and the remaining part is stored in the collection SectionB new ;

[0075] (7-3) Read the set Pair Block Any element (a m ,b m ), read the virtual stratum set SA virtual , SB virtual The corresponding stratum sa m 、sb m ;

[0076] (7-4) Give the stratum sa m 、sb m The same new stratigraphic number and the stratigraphic sa m Deposit into Section A new , the formation sb m Deposit into Section B new ;

[0077] (7-5) Loop through steps (7-3)-(7-4) until Pair Block All elements in are traversed;

[0078] (7-6) Section A new Section B new Generate a section A with consistent topological relationship fit , Section B fit , Section A fit , Section B fit The results are after consistency processing.

[0079] The device for automatically processing the consistency of stratum boundaries in adjacent geological profiles described in the present invention includes a processor and a computer program stored in a memory and executable on the processor, and the above method is implemented when the processor executes the program.

[0080] Beneficial effects: Compared with the prior art, the significant advantage of the present invention is that the present invention realizes an automatic processing method for the consistency of stratum boundary lines in adjacent geological sections based on the known correspondence of strata in adjacent geological sections and according to the topological characteristics of the stratum geometric figures in the sections. The method of the present invention realizes the automatic consistency processing of boundary lines, with high automation, high efficiency and high quality, and can facilitate the subsequent three-dimensional modeling of connected geological sections. The method of the present invention has certain application value and research significance for the realization of the three-dimensional modeling method of auxiliary connection address sections and the full use of adjacent geological sections for modeling. BRIEF DESCRIPTION OF THE DRAWINGS

[0081] Figure 1 is the adjacent geological profile vector data map used in this embodiment;

[0082] Figure 2 It is a flow chart of a method for automatically processing the uniformity of stratum boundaries in adjacent geological sections provided by the present invention;

[0083] Figure 3 This is a schematic diagram of the adjacent geological section after extracting the stratum boundary line and endpoints in this embodiment;

[0084] Figure 4 is a diagram of the solution process of the only topologically inconsistent boundary line in this embodiment;

[0085] Figure 5 It is a schematic diagram of the topological relationship on the merged boundary line corresponding to different situations;

[0086] Figure 6 It is the result of this embodiment, and the schematic diagram corresponding to the virtual stratum and the boundary line of the virtual stratum. DETAILED DESCRIPTION

[0087] The technical solution of the present invention is further described in detail below. In this embodiment, the adjacent geological profile of Longtou Mountain in a mining area is selected as experimental data. Figure 1 As shown, the projection coordinate system used by the experimental data is CGCS2000. The following is further explained by describing a specific embodiment in conjunction with the accompanying drawings.

[0088] like Figure 2 As shown, the method for automatically processing the uniformity of stratum boundaries in adjacent geological sections provided in this embodiment includes the following steps:

[0089] (1) Load two adjacent geological sections to be processed, record them as Section A and Section B, and obtain the stratigraphic code set FD of all strata, the stratigraphic set SectionA of Section A, and the stratigraphic set SectionB of Section B.

[0090] This step specifically includes:

[0091] (1-1) The two adjacent geological sections to be processed are recorded as Section A and Section B, the vector layers of Section A and Section B are read, and the stratigraphic codes of all stratigraphic layers are read from the vector layers and stored in the stratigraphic code set FD = {f i |i=1,2,…,fn}, where f i represents the i-th stratum code, and fn represents the number of stratum codes; in this embodiment, fn=5;

[0092] (1-2) Read all the stratigraphic data in section A, renumber the stratigraphic data, and store them in the stratigraphic set Section A = {StraA a |a=1,2,…,San}, where StraA a represents the ath stratum, San represents the number of strata in Section A; in this embodiment, San=5;

[0093] (1-3) Read all the stratigraphic data in section B, renumber the stratigraphic data, and store them in the stratigraphic set Section B = {StraB b |b=1,2,…,Sbn}, where StraB b represents the bth stratum, and Sbn represents the number of strata in Section B. In this embodiment, Sbn=5. Generally speaking, the adjacent stratum sections in the embodiment need to ensure that San=Sbn=fn.

[0094] (2) Extract the boundary lines and boundary line endpoints of SectionA and SectionB respectively, and store them in boundary line sets ArcsA and ArcsB and boundary line endpoint sets NodesA and NodesB.

[0095] This step specifically includes:

[0096] (2-1) Read any stratum StraA from the stratum set SectionA i ;

[0097] (2-2) Read all other strata StraA from the stratum set SectionA j , j=1,2,...,San, and j≠i; and analyze StraA i With StraA j If two strata are in contact, the common boundary line arc of the two strata is obtained. t ;

[0098] (2-3) Check arc t Does it exist in the boundary line set ArcsA? If not, set arc t Store in the boundary line set ArcsA;

[0099] (2-4) Steps (2-1) to (2-3) are executed repeatedly until all strata in SectionA are traversed, and a complete boundary line set ArcsA is obtained. In this embodiment, there are 12 boundary lines in ArcsA.

[0100] (2-5) Read an arbitrary boundary line arc from the boundary line set ArcsA t , and extract arct The first and last endpoints p 1 、p 2 ;

[0101] (2-6) Check p 1 、p 2 Whether it exists in the boundary endpoint set NodesA, if not, then store it in the boundary endpoint set NodesA;

[0102] (2-7) Loop through steps (2-5) to (2-6) until ArcsA is completely traversed, and store all boundary line endpoints in NodesA. In this embodiment, there are 8 midpoints in NodesA.

[0103] (2-8) For set SectionB, the same process is performed according to steps (2-1) to (2-7) to obtain a boundary line set ArcsB and a boundary line endpoint set NodesB. In this embodiment, ArcsA, NodesA, ArcsB, and NodesB contain elements such as Figure 3 shown.

[0104] (3) According to the geometric properties of SectionA, ArcsA, and NodesA, the topological contact relationship table FA1 between the boundary line and the stratum, the topological contact relationship table FA2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FA3 between the strata in Section A are established; for SectionB, ArcsB, and NodesB, the topological contact relationship table FB1 between the boundary line and the stratum, the topological contact relationship table FB2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FB3 between the strata in Section B are obtained by the same processing.

[0105] This step specifically includes:

[0106] (3-1) Read any boundary line arc in ArcsA t ;

[0107] (3-2) Arc t Check the topological contact relationship with all the strata in SectionA. If arc t StraA a Contact, then StartA a The number a is stored in the arc in the topological contact relationship table FA1 t The corresponding line;

[0108] (3-3) Steps (3-1) to (3-2) are executed repeatedly until all boundary lines in ArcsA are traversed, and a topological contact relationship table FA1 recording the contact stratum numbers of all boundary lines is obtained, as shown in Table 1;

[0109] Table 1 Relationship between the boundary line of section A and the topological structure of the strata

[0110]

[0111] The boundary lines are numbered starting from 0 in step (2) according to the cyclic sequence. In this embodiment, the boundary lines of profile A are numbered from 0 to 11. Usually, a boundary line is in contact with two strata. When the boundary line contacts only one stratum, a number -1 is added to the contact stratum sequence to make the number of boundary line strata uniform.

[0112] (3-4) Read any point p in NodesA k ;

[0113] (3-5) for p k Perform a topological check with all boundary lines in the set ArcsA. If p k On a certain boundary arc t On, the arc t The number t is stored in the point p in the topological contact relationship table FA2 k The corresponding line;

[0114] (3-6) Steps (3-4) to (3-5) are executed repeatedly until all points in NodesA are traversed, and a topological contact relationship table FA2 recording the contact boundary line numbers of all boundary line endpoints is obtained, as shown in Table 2 below;

[0115] Table 2 Topological relationship between the boundary line and endpoints of section A

[0116]

[0117] (3-7) Take any stratum StraA from the set SectionA a ;

[0118] (3-8) Take all the remaining strata StraA from the set SectionA a’ , where a≠a'; for StraA a With StraA a’ Perform contact topology check. If there is a contact relationship, store (a, a') in the adjacency relationship table between layers FA3.

[0119] (3-9) Steps (3-7) to (3-8) are executed repeatedly until all strata in Section A are traversed, and an adjacency table FA3 of all strata in Section A is obtained, as shown in Table 3 below;

[0120] Table 3 Stratigraphic adjacency relationship of section A

[0121]

[0122]

[0123] If a stratum has a boundary that does not touch any other stratum, a virtual stratum StraA is added to the adjacency table. -1 , ensuring that the features of the strata that have certain boundaries and do not contact other strata can be recorded.

[0124] (3-10) For SectionB, ArcsB and NodesB, the same processing is performed according to steps (3-1)-(3-7) to obtain the topological contact relationship table FB1, the topological contact relationship table FB2, and the adjacency relationship table FB3.

[0125] (4) Generate combinations of all connected strata in profiles A and B based on adjacency tables FA3 and FB3, and merge the connected strata of each combination respectively. The boundary lines of the merged strata are stored in merged stratum boundary line sets La and Lb respectively.

[0126] This step specifically includes:

[0127] (4-1) Store all the stratum numbers in SectionA and SectionB into the stratum number set Ser = {ser n |n=1,2,…,N}, where ser n is the number of the nth stratum, N is the total number of strata; N = 5

[0128] (4-2) For the stratum number set Ser, generate all combinations of any stratum numbers and store them in the set Comb = {comb m |m=1,2,…,M}, where comb m is the mth stratum number combination, M is the total number of sets; in this embodiment, M=30;

[0129] (4-3) Read a combination comb in the set Comb m ;

[0130] (4-4) According to the adjacency table FA3, get comb m All strata in comb m The adjacent situation of other strata in the stratum is stored in the set linka m ={(u,v)|u,v∈comb m , (u, v)∈FA3}, where u and v are the numbers of adjacent strata respectively, and (u, v) means that the strata u and v are adjacent. In this embodiment, the situation is as follows Figure 4 When shown, comb m ={4,5};

[0131] (4-5) Take the strata in sectionA as nodes and linka m The adjacent relationships in are used as edges to construct and query the unionSet a , and linka m All relations in the unionSet are added to the unionSet a ;

[0132] (4-6) For the stratum adjacency table FB3, follow the same process as steps (4-4)-(4-5) to obtain the set linkb m ={(u,v)|u,v∈comb m , (u,v)∈FB3} and unionSet b ;

[0133] (4-7) Check comb m Do all the stratum numbers in unionSeta belong to the same parent node, and do they belong to the same parent node in unionSeta? b belong to the same parent node. If they all belong to the same parent node, then comb m All strata corresponding to SectionA are merged into a new stratum StraA m , comb m All the strata corresponding to SectionB are merged into a new stratum StraB m In this embodiment, the situation is as follows Figure 4 When shown, StraA m The gray oblique line area in section A, StraB m It is the gray oblique line area in section B;

[0134] (4-8) Extraction of StraA m The merged boundary line ea is stored in the merged boundary line set La = {ea m |m≤M}, where ea m Indicates the merged boundary line numbered m, where m is the number of comb m Corresponding combination number in Comb; extract StraB m The merged boundary line eb is stored in the merged boundary line set Lb = {eb m |m≤M}, where eb m represents the merged boundary line numbered m; in this embodiment, the situation is as follows Figure 4 As shown, the thick black line (and hollow line) in section A shows ea m The thick black line (and hollow line) in section B shows eb m ;

[0135] (4-9) Loop through steps (4-3) to (4-8) until all combinations in Comb are traversed and a complete set of merged boundary lines La, Lb is generated.

[0136] (5) Take out any two corresponding merged boundary lines in La and Lb, and check whether the stratum boundary lines before merging correspond to each other according to FA1, FA2, FB1, and FB2, and obtain the set La after the inconsistent lines are grouped. block , Lb block , the corresponding relationship is stored in the relationship table Pair Block .

[0137] This step specifically includes:

[0138] (5-1) Read any merged boundary line ea in La m , and the merged boundary line eb with the same number in Lb m In this embodiment, the situation is as follows Figure 4 As shown, the thick black line (and hollow line) in section A shows ea m The thick black line (and hollow line) in section B shows eb m ;

[0139] (5-2) Traverse all boundary lines in ArcsA and compare them with ea m Completely overlapping boundary lines are stored in the collection arcs ea ; Traverse all the boundary lines in ArcsB and compare them with eb m All completely overlapping boundary lines are stored in the collection arcs eb In this embodiment, the situation is as follows Figure 4 When shown, arcs ea The set of boundary line numbers is {3, 6, 7, 1, 2}, arcs eb The set of middle boundary line numbers is {1, 6, 7, 2, 3}, which are the actual stratigraphic boundary lines that completely coincide with the merged boundary lines;

[0140] (5-3) Set arcs ea All boundary lines in the set arcs are reordered according to the connection relationship in the topological contact relationship table FA2 so that the boundary lines with adjacent numbers are connected; eb All boundary lines in are reordered according to the connection relationship in the topological contact relationship table FB2 so that the boundary lines with adjacent numbers are connected;

[0141] (5-4) Comparison set arcs ea With collection arcs ebThe boundary lines in FA1 and FB1 are matched according to the contact relationship between the boundary lines and the strata. If a unique set of topologically inconsistent and continuous boundary lines can be matched, they are stored in the set la. m lb m , and m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship (a m ,b m ) is stored in the relationship set Pair Block In this embodiment, the situation is as follows Figure 4 As shown, the hollow line in section A shows la m , the hollow line in section B shows lb m ;

[0142] (5-5) Loop through steps (5-1) to (5-4) until all sets La and L are traversed. b All merged boundary lines in the , get the complete corresponding relationship set Pair Block In this embodiment, Pair Block The number of elements of is 1.

[0143] like Figure 4 and Figure 5 As shown, (5-4) specifically includes:

[0144] (5-4-1) Traversing arcs ea middle boundary line until a boundary line l is obtained t1 , whose contact stratigraphic number is the same as arcs eb A boundary line with a unique number in a contact stratum t2 same;

[0145] (5-4-2) Collect arcs ea The order of the elements in is adjusted to l t1 As the first element; set arcs eb The order of the elements in is adjusted to l t2 is the first element;

[0146] (5-4-3) From the collection arcs ea Read an element l from front to back ea , according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs ebRead an element l from front to back eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ;

[0147] (5-4-4) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2 ) are the same, if they are the same, return to step (5-4-3), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa1, pb1 in;

[0148] (5-4-5) Check whether arcs have been traversed ea or arcs eb If the traversal is not completed, continue to execute step (5-4-6); if the traversal is completed, Figure 5 In case (a) of the above, all the actual boundary lines of the two combined boundary lines correspond to each other completely, or Figure 5 The special case of the case (b) in the above example is that there is only a topological inconsistency at the end and the corresponding position at the end of the other merged boundary line is empty. In short, this indicates that the merged boundary line does not contain a unique topologically inconsistent boundary line, and no new elements are stored in the set La. block , set Lb block , set Pair Block , execute step (5-5);

[0149] (5-4-6) From the collection arcs ea Read an element l from back to front ea , according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs eb Read an element l from back to front eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ;

[0150] (5-4-7) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2) are the same, if they are the same, return to step (5-4-6), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa2, pb2 in;

[0151] (5-4-8) If pa2-pa1=0 or pb2-pb1=0, it means that the line does not contain a unique topologically inconsistent continuous boundary line, and no new elements are stored in the set La. block , set Lb block , set Pair Block , execute step (5-5); otherwise, continue to execute step (5-4-9);

[0152] (5-4-9) If pa2-pa1>1 and pb2-pb1>1, then traverse arcs ea The elements in the interval [pa1+1,pa2-1] can be eb If the two contacting strata are found to have exactly the same elements in the interval [pb1+1, pb2-1], it means that the interval [pb1+1, pb2-1] contains multiple sets of topologically inconsistent continuous boundary lines. The positions of the elements of the contacting strata found to have exactly the same elements in profile A and profile B are recorded as pat and pbt, and step (5-4-10) is executed; otherwise, step (5-4-11) is executed;

[0153] (5-4-10) will gather arcs ea Divide into sets arcs based on position pat ea1 , arcs ea2 The left and right sets will set arcs eb Divide into sets arcs based on position pbt eb1 , arcs eb2 The left and right sets are respectively for arcs ea1 、arcs eb1 , arcs ea2 、arcs eb2 Execute step (5-4) to obtain the topological inconsistency positions in the left and right sets;

[0154] (5-4-11) Arcs ea The elements in the interval [pa1+1,pa2-1] are stored in the set la m ; Arcs eb The elements in the interval [pb1+1,pb2-1] are stored in the set lb m ;like Figure 5In the case of (d), the interval [pa1+1, pa2-1] is a set of topologically inconsistent continuous boundary lines. In this embodiment, Figure 4 In this case, the hollow lines on sections A and B are a set of continuous boundary lines with inconsistent topology.

[0155] (5-4-12) m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship is stored in the set Pair Block .

[0156] (6) For La block , Lb block For each element in, a virtual stratum is generated and stored in the virtual stratum set SA virtual , SB virtual .

[0157] This step specifically includes:

[0158] (6-1) Read La block Any set in la m ;

[0159] (6-2) Create a set la m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SA virtual In this embodiment, SA virtual The number of elements is N v = 1, the virtual stratum in section A is as follows Figure 6 The middle section A is shown as a virtual stratum;

[0160] (6-3) Read Lb block Any set lb m ;

[0161] (6-4) Create a set lb m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SB virtual The virtual strata in section B are as follows: Figure 6 The virtual strata are shown in section B.

[0162] (7) Set the virtual stratum SA virtual , SB virtual The virtual stratum in the relationship table PairBlock Renumber and merge with SectionA and SectionB to form Section A with consistent topological relationship fit , Section B fit , Section A fit , Section B fit The results are after consistency processing.

[0163] This step specifically includes:

[0164] (7-1) The strata in Section A and SA virtual The overlapping parts of the virtual strata are erased, and the remaining parts are stored in the collection SectionA new ;

[0165] (7-2) The strata in Section B and SB virtual The overlapping part of the virtual stratum is erased, and the remaining part is stored in the collection SectionB new ;

[0166] (7-3) Read the set Pair Block Any element (a m ,b m ), read the virtual stratum set SA virtual , SB virtual The corresponding stratum sa m 、sb m ;

[0167] (7-4) Give the stratum sa m 、sb m The same new stratigraphic number and the stratigraphic sa m Deposit into Section A new , the formation sb m Deposit into Section B new In this embodiment, the virtual formation sa m 、sb m like Figure 6 As shown in , it can completely correspond to the boundary line of the adjacent strata one by one, and the topological relationship is consistent;

[0168] (7-5) Loop through steps (7-3)-(7-4) until Pair Block All elements in are traversed;

[0169] (7-6) Section A new Section B new Generate a section A with consistent topological relationship fit , Section B fit , Section A fit , Section Bfit This is the result after the consistency processing. Figure 6 As shown in Section A new Section B new The boundary lines of the virtual strata created in this embodiment correspond to the following: Figure 6 shown.

[0170] The device for automatically processing the consistency of stratum boundaries in adjacent geological profiles described in the present invention includes a processor and a computer program stored in a memory and executable on the processor, and the above method is implemented when the processor executes the program.

[0171] In the embodiment of the present invention, some GIS operations are provided based on the open source library GDAL, and the relevant steps can also use the API of software such as SuperMap and ArcGIS Object to perform corresponding GIS operations.

[0172] The above disclosure is only a preferred embodiment of the present invention, which cannot be used to limit the scope of the present invention. Therefore, equivalent changes made according to the claims of the present invention are still within the scope of the present invention.

Claims

1. A method for automatically processing the boundaries of adjacent strata in geological sections, characterized in that The method includes: (1) Load two adjacent geological sections to be processed, record them as Section A and Section B, obtain the stratigraphic code set FD of all strata, the stratigraphic set SectionA of Section A, and the stratigraphic set SectionB of Section B; (2) Extract the boundary lines and boundary line endpoints of Section A and Section B respectively, and store them in boundary line sets ArcsA and ArcsB and boundary line endpoint sets NodesA and NodesB; (3) According to the geometric attributes of SectionA, ArcsA, and NodesA, the topological contact relationship table FA1 between the boundary line and the stratum, the topological contact relationship table FA2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FA3 between the strata in Section A are established; for SectionB, ArcsB, and NodesB, the topological contact relationship table FB1 between the boundary line and the stratum, the topological contact relationship table FB2 between the boundary line endpoints and the boundary line, and the adjacency relationship table FB3 between the strata in Section B are obtained by the same processing; (4) Generate all combinations of connected strata in profile A and profile B according to adjacency tables FA3 and FB3, and merge the connected strata of each combination respectively. The boundary lines of the merged strata are stored in the merged strata boundary line sets La and Lb respectively; (5) Take out any two corresponding merged boundary lines in La and Lb, and check whether the stratum boundary lines before merging correspond to each other according to FA1, FA2, FB1, and FB2, and obtain the set La after the inconsistent lines are grouped. block , Lb block , the corresponding relationship is stored in the relationship table Pair Block ; (6) For La block , Lb block For each element in, a virtual stratum is generated and stored in the virtual stratum set SA virtual , SB virtual ; (7) Set the virtual stratum SA virtual , SB virtual The virtual stratum in the relationship table Pair Block Renumber and merge with SectionA and SectionB to form Section A with consistent topological relationship fit , Section B fit , Section A fit , Section B fit To make the results consistent; Wherein, step (3) specifically includes: (3-1) Read any boundary line arc in ArcsA t ; (3-2) Arc t Check the topological contact relationship with all the strata in SectionA. If arc t StraA a Contact, then StartA a The number a is stored in the arc in the topological contact relationship table FA1 t The corresponding line; (3-3) looping through steps (3-1) to (3-2) until all boundary lines in ArcsA are traversed, and obtaining a topological contact relationship table FA1 that records the contact stratum numbers of all boundary lines; (3-4) Read any point p in NodesA k ; (3-5) for p k Perform a topological check with all boundary lines in the set ArcsA. If p k On a certain boundary arc t On, the arc t The number t is stored in the point p in the topological contact relationship table FA2 k The corresponding line; (3-6) Loop through steps (3-4) to (3-5) until all points in NodesA are traversed, and obtain a topological contact relationship table FA2 that records the contact boundary line numbers of all boundary line endpoints; (3-7) Take any stratum StraA from the set SectionA a ; (3-8) Take all the remaining strata StraA from the set SectionA a’ , where a≠a'; for StraA a With StraA a’ Perform contact topology check. If there is a contact relationship, (a, a') is stored in the adjacency relationship table FA3 between layers. (3-9) Loop through steps (3-7) to (3-8) until all strata in Section A are traversed, and obtain the adjacency table FA3 of all strata in Section A; (3-10) For SectionB, ArcsB and NodesB, the same processing is performed according to steps (3-1)-(3-7) to obtain the topological contact relationship table FB1, the topological contact relationship table FB2, and the adjacency relationship table FB3.

2. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (1) specifically includes: (1-1) The two adjacent geological sections to be processed are recorded as Section A and Section B, the vector layers of Section A and Section B are read, and the stratigraphic codes of all stratigraphic layers are read from the vector layers and stored in the stratigraphic code set FD = {f i |i=1,2,…,fn}, where f i represents the i-th stratum code, and fn represents the number of stratum codes; (1-2) Read all the stratigraphic data in section A, renumber the stratigraphic data, and store them in the stratigraphic set Section A = {StraA a |a=1,2,…,San}, where StraA a represents the ath stratum, and San represents the number of strata in SectionA; (1-3) Read all the stratigraphic data in section B, renumber the stratigraphic data, and store them in the stratigraphic set Section B = {StraB b |b=1,2,…,Sbn}, where StraB b represents the b-th stratum, and Sbn represents the number of strata in SectionB.

3. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (2) specifically includes: (2-1) Read any stratum StraA from the stratum set SectionA i ; (2-2) Read all other strata StraA from the stratum set SectionA j , j=1,2,...,San, and j≠i; and analyze StraA i With StraA j If two strata are in contact, the common boundary line arc of the two strata is obtained. t ; (2-3) Check arc t Does it exist in the boundary line set ArcsA? If not, set arc t Store in the boundary line set ArcsA; (2-4) Steps (2-1) to (2-3) are executed repeatedly until all strata in SectionA are traversed and a complete boundary line set ArcsA is obtained; (2-5) Read an arbitrary boundary line arc from the boundary line set ArcsA t , and extract arc t The first and last endpoints p1 and p2 on the (2-6) Check whether p1 and p2 exist in the boundary endpoint set NodesA. If not, store them in the boundary endpoint set NodesA. (2-7) Loop through steps (2-5) to (2-6) until ArcsA is completely traversed and all boundary endpoints are stored in NodesA; (2-8) For set SectionB, perform the same processing as steps (2-1)-(2-7) to obtain the boundary line set ArcsB and the boundary line endpoint set NodesB.

4. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (4) specifically includes: (4-1) Store all the stratum numbers in SectionA and SectionB into the stratum number set Ser = {ser n |n=1,2,…,N}, where ser n is the number of the nth stratum, and N is the total number of strata; (4-2) For the stratum number set Ser, generate all combinations of any stratum numbers and store them in the set Comb = {comb m |m=1,2,…,M}, where comb m is the mth stratigraphic number combination, M is the total number of sets; (4-3) Read a combination comb in the set Comb m ; (4-4) According to the adjacency table FA3, get comb m All strata in comb m The adjacent situation of other strata in the stratum is stored in the set linka m ={(u,v)|u,v∈comb m , (u,v)∈FA3}, where u and v are the numbers of adjacent strata respectively, and (u,v) means that strata u and v are adjacent; (4-5) Take the strata in sectionA as nodes and linka m The adjacent relationships in are used as edges to construct and query the unionSet a , and linka m All relations in the unionSet are added to the unionSet a ; (4-6) For the stratum adjacency table FB3, follow the same process as steps (4-4)-(4-5) to obtain the set linkb m ={(u,v)|u,v∈comb m , (u,v)∈FB3} and unionSet b ; (4-7) Check comb m Are all the stratum numbers in the unionSet a belong to the same parent node and are in the unionSet b belong to the same parent node. If they all belong to the same parent node, then comb m All strata corresponding to SectionA are merged into a new stratum StraA m , comb m All the strata corresponding to SectionB are merged into a new stratum StraB m ; (4-8) Extraction of StraA m The merged boundary line ea is stored in the merged boundary line set La = {ea m |m≤M}, where ea m Indicates the merged boundary line numbered m, where m is the number of comb m Corresponding combination number in Comb; extract StraB m The merged boundary line eb is stored in the merged boundary line set Lb = {eb m |m≤M}, where eb m represents the merged boundary line numbered m; (4-9) Loop through steps (4-3) to (4-8) until all combinations in Comb are traversed and a complete set of merged boundary lines La, Lb is generated.

5. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (5) specifically includes: (5-1) Read any merged boundary line ea in La m , and the merged boundary line eb with the same number in Lb m ; (5-2) Traverse all boundary lines in ArcsA and compare them with ea m Completely overlapping boundary lines are stored in the collection arcs ea ; Traverse all the boundary lines in ArcsB and compare them with eb m All completely overlapping boundary lines are stored in the collection arcs eb ; (5-3) Set arcs ea All boundary lines in the set arcs are reordered according to the connection relationship in the topological contact relationship table FA2 so that the boundary lines with adjacent numbers are connected; eb All boundary lines in are reordered according to the connection relationship in the topological contact relationship table FB2 so that the boundary lines with adjacent numbers are connected; (5-4) Comparison set arcs ea With collection arcs eb The boundary lines in FA1 and FB1 are matched according to the contact relationship between the boundary lines and the strata. If a unique set of topologically inconsistent and continuous boundary lines can be matched, they are stored in the set la. m lb m , and m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship (a m ,b m ) is stored in the relationship set Pair Block ; (5-5) Loop through steps (5-1) to (5-4) until all merged boundary lines in sets La and Lb are traversed to obtain a complete set of corresponding relationships Pair Block .

6. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (5-4) specifically includes: (5-4-1) Traversing arcs ea middle boundary line until a boundary line l is obtained t1 , whose contact stratigraphic number is the same as arcs eb A boundary line with a unique number in a contact stratum t2 same; (5-4-2) Collect arcs ea The order of the elements in is adjusted to l t1 As the first element; set arcs eb The order of the elements in is adjusted to l t2 is the first element; (5-4-3) From the collection arcs ea Read an element l from front to back ea , according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs eb Read an element l from front to back eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ; (5-4-4) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2 ) are the same, if they are the same, return to step (5-4-3), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa1, pb1 in; (5-4-5) Check whether arcs have been traversed ea or arcs eb If the traversal is not completed, continue to execute step (5-4-6); if the traversal is completed, it means that the actual boundary lines on the merged boundary line can correspond one to one, or there is only a topological inconsistency at the end and the corresponding position at the end of the other merged boundary line is empty. At the same time, it means that the line does not contain a unique continuous boundary line with inconsistent topology, and no new elements are stored in the set La block , set Lb block , set Pair Block , execute step (5-5); (5-4-6) From the collection arcs ea Read an element l from back to front ea , according to the topological contact relationship table FA1 between the boundary line and the stratum, the ea Contact stratum d ea1 ,d ea2 ; From the collection arcs eb Read an element l from back to front eb According to the topological contact relationship table FB1 between the boundary line and the stratum, the eb Contact stratum d eb1 ,d eb2 ; (5-4-7) Judgment (d ea1 ,d ea2 ) and (d eb1 ,d eb2 ) are the same, if they are the same, return to step (5-4-6), if they are different, stop the loop and record the element read when stopping in arcs ea 、arcs eb Position pa2, pb2 in; (5-4-8) If pa2-pa1=0 or pb2-pb1=0, it means that the line does not contain a unique topologically inconsistent continuous boundary line, and no new elements are stored in the set La. block , set Lb block , set Pair Block , execute step (5-5); otherwise, continue to execute step (5-4-9); (5-4-9) If pa2-pa1>1 and pb2-pb1>1, then traverse arcs ea The elements in the interval [pa1+1,pa2-1] can be eb If the two contacting strata are found to have exactly the same elements in the interval [pb1+1, pb2-1], it means that the interval [pb1+1, pb2-1] contains multiple groups of topologically inconsistent continuous boundary lines. The positions of the elements of the contacting strata found to have exactly the same elements in profile A and profile B are recorded as pat and pbt, and step (5-4-10) is executed; otherwise, step (5-4-11) is executed; (5-4-10) will gather arcs ea Divide into sets arcs based on position pat ea1 , arcs ea2 The left and right sets will set arcs eb Divide into sets arcs based on position pbt eb1 , arcs eb2 The left and right sets are respectively for arcs ea1 、arcs eb1 , arcs ea2 、arcs eb2 Execute step (5-4) to obtain the topological inconsistency positions in the left and right sets; (5-4-11) Arcs ea The elements in the interval [pa1+1,pa2-1] are stored in the set la m ; Arcs eb The elements in the interval [pb1+1,pb2-1] are stored in the set lb m ; (5-4-12) m Deposit into collection La block , lb m Deposit into collection Lb block , will la m lb m The corresponding relationship is stored in the set Pair Block .

7. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (6) specifically includes: (6-1) Read La block Any set in la m ; (6-2) Create a set la m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SA virtual ; (6-3) Read Lb block Any set lb m ; (6-4) Create a set lb m Buffer zone s with distance less than the preset value m , as la m The virtual strata are stored in the virtual strata set SB virtual .

8. The method for automatically processing the uniformity of stratum boundaries in adjacent geological sections according to claim 1, characterized in that: Step (7) specifically includes: (7-1) The strata in Section A and SA virtual The overlapping parts of the virtual strata are erased, and the remaining parts are stored in the collection SectionA new ; (7-2) The strata in Section B and SB virtual The overlapping part of the virtual stratum is erased, and the remaining part is stored in the collection SectionB new ; (7-3) Read the set Pair Block Any element (a m ,b m ), read the virtual stratum set SA virtual , SB virtual The corresponding stratum sa m 、sb m ; (7-4) Give the stratum sa m 、sb m The same new stratigraphic number and the stratigraphic sa m Deposit into Section A new , the formation sb m Deposit into Section B new ; (7-5) Loop through steps (7-3)-(7-4) until Pair Block All elements in are traversed; (7-6) Section A new Section B new Generate a section A with consistent topological relationship fit , Section B fit , Section A fit , Section B fit The results are after consistency processing.

9. An apparatus for automatically processing the uniformity of stratum boundaries in adjacent geological sections, comprising a processor and a computer program stored in a memory and executable on the processor, characterized in that: When the processor executes the program, the method according to any one of claims 1 to 8 is implemented.

Images