Automatic Analysis Method and Device for Geological Events Based on Geological Profiles

Abstract

This invention discloses an automatic analysis method and apparatus for geological events based on geological profiles. The method includes: (1) reading the geological profile, fault lines, and a list of sedimentary strata in the study area; (2) serializing the surface exposed strata of the profile into a coded sequence; (3) analyzing, reducing, and recording Quaternary sedimentary formation events; (4) analyzing and reducing igneous rock strata and recording intrusion events; (5) detecting faults and determining their types, and then recording the faults and their effects during movement; (6) identifying and recording folds according to the principle of symmetry; and (7) formatting the list of geological events and saving it as a local file. This invention can automatically analyze Quaternary sedimentary events, intrusion events, fault events, and fold deformation events in geological profiles with basic accuracy, and automatically analyze the related strata, with high efficiency and low cost.

Description

Technical Field

[0001] This invention relates to the fields of geographic information and geology, and in particular to an automatic analysis method and apparatus for geological events based on geological profiles. Background Technology

[0002] A geological profile is a map, drawn to a specific scale, depicting geological objects, phenomena, and their interrelationships at a particular location. There are many types of geological profiles. Based on the angle with the strike of strata, they can be divided into longitudinal and transverse profiles. Based on the main content, they can be classified as stratigraphic profiles, intrusive rock profiles, structural profiles, Quaternary profiles, volcanic rock profiles, etc. Based on the drawing method, they can be divided into freehand profiles, measured profiles, and map-cut profiles. Geological profile maps effectively reflect the stratigraphic structure and rock mass characteristics of the study area, facilitating the three-dimensional and intuitive recording and study of geological objects. They are a commonly used type of map for geological experts to record their findings.

[0003] Geological events are geological records that occur suddenly or have special significance within a region during a specific geological period. The analysis of geological events is fundamental to further research into the region's geological history and geological tectonic formation. Based on geoscientific knowledge and expert experience, local geological structures and their corresponding geological events can be inferred from geological profiles and stratigraphic development patterns. Currently, further analysis of geological events based on geological profiles mainly relies on expert experience and is conducted manually. However, with the continuous production and accumulation of measured geological profiles, and the emergence of software that automatically generates map sections from geological maps, manual analysis of geological events is relatively inefficient, costly, and the quality of analysis varies from person to person. Summary of the Invention

[0004] Purpose of the invention: This invention addresses the problems existing in the prior art by providing a highly efficient and low-cost method and apparatus for automatic analysis of geological events based on geological profiles.

[0005] Technical solution: The automatic analysis method for geological events based on geological profiles described in this invention includes the following steps:

[0006] (1) Load the geological profile vector layer, read the strata and store them in the strata set Stra, and at the same time obtain the strata attribute information table FD; load the fault vector layer and store it in the fault set Faults; load all sedimentary strata in the study area and store them in the sedimentary stratigraphic table Sttable and encode each sedimentary strata to form the stratigraphic coding table SF.

[0007] (2) Generate a coding sequence DS consisting of the codes of all exposed strata in the stratigraphic set Stra based on the stratigraphic coding table SF;

[0008] (3) Traverse the coding sequence DS, store the events corresponding to the Quaternary strata into the geological event set GeoEvents in chronological order, store the non-Quaternary strata in the stratigraphic set Stra into the set Stra2, and delete the Quaternary strata from the coding sequence DS.

[0009] (4) Traverse the coding sequence DS, store the intrusion events corresponding to the intrusion strata into the geological event set GeoEvents in chronological order, store the non-intrusion rock strata in Stra2 into the stratum set Stra3, and delete the intrusion strata from the coding sequence DS.

[0010] (5) Traverse the fault set Faults, obtain the adjacent strata on both sides of each fault in Stra3, determine the type of fault and fault event based on the strata on both sides of the fault, and then store it in the geological event set GeoEvents, and add or delete the strata codes affected by the fault in the coding sequence DS accordingly.

[0011] (6) Based on the coding sequence DS, the symmetrical repeating fold strata are searched and stored in the geological event set GeoEvents according to the principle of symmetric repetition;

[0012] (7) Output the GeoEvents collection in a formatted manner.

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

[0014] (1-1) Load the geological profile vector layer and store it in the stratigraphic set Stra = {stra a |a=1,2,…,sn}, where stra a Let 'a' represent the 'a'-th stratum and 'sn' represent the stratum number; sequentially read the stratigraphic codes of all strata in Stra and store them in the stratigraphic code set DC = {dc} a |a=1,2,…,sn} where dc a This represents the stratigraphic designation of the a-th stratum.

[0015] (1-2) Obtain the dip angle and geological age of each stratum in the stratigraphic set Stra, and store them in the stratigraphic attribute table FD;

[0016] (1-3) Load the fault vector layer and store it in the fault set Faults = {f b |b=1,2,…,fn}, where f b Let f_b represent the b-th fault, and f_n represent the number of faults.

[0017] (1-4) Load all sedimentary strata in the study area and store them in the sedimentary strata table Stable={st i |i=0,1,…,m}, where sti This represents the stratigraphic code of the i-th sedimentary stratum from newest to oldest, where m is the total number of all sedimentary strata in the current region.

[0018] (1-5) According to the stratigraphic sequence in Stable, each sedimentary stratum is stratigraphically encoded, that is, a single character is used to represent a sedimentary stratum, forming the stratigraphic encoding table SF.

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

[0020] (2-1) Traverse the strata in the stratigraphic set Stra and check whether they are exposed on the surface. Remove the strata that are not exposed from the stratigraphic set Stra.

[0021] (2-2) Rearrange the strata in the stratigraphic set Stra in order from left to right, and at the same time rearrange the contents of FD in the stratigraphic attribute table in this order;

[0022] (2-3) Based on the stratigraphic coding table SF, map each stratigraphic layer in Stra to a code, and form the coding sequence DS = {ds} in sequence. c |c=1,2,…,cn} where, ds c It is the code of the c-th stratum, and cn is the number of strata.

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

[0024] (3-1) Read any element ds from the encoded sequence DS c Obtain the corresponding stratigraphic set Stra from the stratigraphic set Stra. a and stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ;

[0025] (3-2) If the stratigraphic code is dc a If the stratum begins with Q, then it is a Quaternary stratum, and proceed to step (3-3); otherwise, store the stratum in the set Stra2 which has been reduced from Quaternary strata, and return to step (3-1).

[0026] (3-3) Stra a With dt a Geological objects and their occurrence or generation time are stored in the geological event set GeoEvents;

[0027] (3-4) Repeat steps (3-1)-(3-3) until all elements in the encoded sequence DS have been traversed;

[0028] (3-5) Delete all Quaternary strata codes in the coding sequence DS and merge adjacent consistent strata codes in the coding sequence DS after the Quaternary strata were removed.

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

[0030] (4-1) Read an element ds from the encoded sequence DS c Get the corresponding stratigraphic stra in set Stra2 a and its stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ;

[0031] (4-2) If the stratigraphic code is dc a If the stratum begins with a Greek letter, it is an intrusive stratum, and step (4-3) is executed; otherwise, the stratum is stored in the set Stra3, which has eliminated intrusive strata, and step (4-1) is returned.

[0032] (4-3) Stra a With dt a Geological objects and their occurrence or generation time are stored in the geological event set GeoEvents;

[0033] (4-4) Repeat steps (4-1)-(4-3) until all elements in all encoded sequences DS have been traversed;

[0034] (4-5) Delete all intrusive rock codes in the coding sequence DS and merge adjacent consistent stratigraphic codes in the coding sequence DS after the intrusive rock strata were eliminated.

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

[0036] (5-1) Read a fault f in the fault set Faults in geological age order from newest to oldest. b And calculate the dip angle of the fault. f ;

[0037] (5-2) Find f in set Stra3 a The adjacent left and right strata are denoted as stra l ,stra r ;

[0038] (5-3) Calculate stra l ,stra r The average angle of inclination ave Compare angles ave with angle fIf the values ​​are equal or the difference is within a preset range, proceed to step (5-4); otherwise, proceed to step (5-5).

[0039] (5-4) If stra l ,stra r If the lithology is consistent, then the fault f b Mark the fault as developed within the strata, otherwise mark it as developed on the stratum contact surface, and store the fault in the GeoEvents set, then execute steps (5-7);

[0040] (5-5) Repeatedly search for continuous layers in the adjacent strata on both sides of the fault. If they are found, the type is defined as reverse fault and stored in the GeoEvents set. Then proceed to step (5-7). Otherwise, the type is defined as normal fault and proceed to step (5-6).

[0041] (5-6) Obtain stra l ,stra r The stratigraphic code is used, and the missing sedimentary strata between the two strata are found in the Stable sedimentary stratigraphic table. The missing sedimentary strata code is then added to the coding sequence DS. l ,stra r Between the stratigraphic codes, the fault and the missing sedimentary strata are stored in the GeoEvents set;

[0042] (5-7) Repeat steps (5-1)-(5-6) until all faults have been traversed;

[0043] (5-8) Iterate through the encoded sequence DS sequentially, delete consecutively repeated stratigraphic subsequences, and keep only the unique one among the repeated subsequences.

[0044] Furthermore, step (5-5) specifically includes:

[0045] (5-5-1) Obtain the adjacent strata on both sides of the fault and store them in an ordered manner in the set Stra Near ={stra j |j=1,2,…,k}, where stra j Let j be the j-th nearest stratum, and k be the total number of nearest strata;

[0046] (5-5-2) In the set Stra Near The search for repeating stratigraphic subsequences requires that at least one subsequence is related to the current fault f. b Adjacency;

[0047] (5-5-3) If a repeating subsequence is found, then the current fault f aThe repeated stratigraphic subsequences are stored as reverse fault events in the GeoEvents set; steps (5-7) are executed; otherwise, the current fault f is... a The type is a normal fault; proceed with steps (5-6).

[0048] Steps (5-6) specifically include:

[0049] (5-6-1) Locate and stra in the Stable sedimentary stratigraphic table l ,stra r Get their index number i l i r ;

[0050] (5-6-2) Obtaining the index of the sedimentary stratigraphic table within the Stable in i l i r The codes of the sedimentary strata between them are stored in the set DS. lr ;

[0051] (5-6-3) Set the current fault f a The type is denoted as a normal fault, and is related to the set DS. lr They are stored together in the GeoEvents collection.

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

[0053] (6-1) Obtain a stratigraphic code ds from the coding sequence DS in sequence. c ;

[0054] (6-2) If ds c ds c+1 If they are not the same stratum, then use ds c Centered on the coding sequence DS, the codes of the strata on both sides are compared one by one until the first inconsistent code is encountered. The coding sequence of the compared strata is a palindrome. The strata of the palindrome are stored as folds in the fold set Folds.

[0055] (6-3) If ds c ds c+1 If they are the same stratum, then use ds c ds c+1 Centered on the coding sequence DS, the stratigraphic codes on both sides are compared one by one until the first inconsistent stratigraphic code is encountered. The comparison stops. The sequence of stratigraphic codes that have been compared is a palindrome. The stratigraphic codes of the palindrome are stored as folds in the fold set Folds.

[0056] (6-4) Repeat steps (6-1)-(6-3) until all palindromes are detected;

[0057] (6-5) Traverse all folds in the fold set Folds, and determine whether the palindrome string corresponds to a syncline or anticline based on the arrangement of the old and new strata in the fold.

[0058] (6-6) Store all folds in the fold set Folds into the geological event set GeoEvents.

[0059] Furthermore, step (6-5) specifically includes:

[0060] (6-5-1) Get a fold from the fold set Folds i This fold contains j strata;

[0061] (6-5-2) Obtaining folds i The geological age of the first stratum in the middle is dt1, and the geological age of the second stratum is dt1. Geological age dt of the strata at the location mid ;

[0062] (6-5-3) Compare dt1 and dt mid The geological time sequence, if dt1 is older, then fold i For synclinal; if dt mid Older ones fold i Anticline; record the anticline / syncline type in fold. i In the object;

[0063] (6-5-4) Repeat steps (6-5-1)-(6-5-3) until all wrinkles have been traversed.

[0064] The automatic geological event analysis device based on geological profiles of the present invention includes a processor and a computer program stored in a memory and capable of running on the processor. When the processor executes the program, it implements the above-mentioned method.

[0065] Beneficial effects: Compared with the prior art, the significant advantages of this invention are: This invention matches the characteristics of different geological events based on the sequential adjacency relationships of exposed strata in a geological profile, thus realizing an automatic geological event analysis method based on geological profiles. This invention has a high degree of automation, high efficiency, and low cost. Attached Figure Description

[0066] Figure 1 This is a flowchart of the automatic analysis method for geological events based on geological profiles provided by the present invention;

[0067] Figure 2 This is a flowchart of the analysis of fracture events and fault types in geological profiles in this invention. Detailed Implementation

[0068] The technical solution of the present invention will be further described in detail below. In this embodiment, a geological profile of Qinglong Mountain in Nanjing City was selected as the experimental data, and the projected coordinate system used for the experimental data is WGS84. The following description, in conjunction with the accompanying drawings, will further illustrate the invention through a specific embodiment.

[0069] Example 1

[0070] This embodiment provides an automatic analysis method for geological events based on geological profiles, such as... Figure 1 As shown, it includes the following steps:

[0071] (1) Load the geological profile vector layer, read the strata and store them in the strata set Stra, and at the same time obtain the strata attribute information table FD; load the fault vector layer and store it in the fault set Faults; load all sedimentary strata in the study area and store them in the sedimentary strata table Sttable and encode each sedimentary strata to form the strata coding table SF.

[0072] This step specifically includes:

[0073] (1-1) Load the geological profile vector layer and store it in the stratigraphic set Stra = {stra a |a=1,2,…,sn}, where stra a Let 'a' represent the 'a'-th stratum and 'sn' represent the stratum number; sequentially read the stratigraphic codes of all strata in Stra and store them in the stratigraphic code set DC = {dc} a |a=1,2,…,sn} where dc a This represents the stratigraphic designation of the a-th stratum; in this embodiment, the profile sn = 53;

[0074] (1-2) Obtain the dip angle and geological age of each stratum in the stratigraphic set Stra, and store them in the stratigraphic attribute table FD; the geological age is recorded by the code of sedimentary rocks, while the age of intrusive rocks should be recorded separately in the attribute table; the stratigraphic attribute table FD is shown in Table 1:

[0075] Table 1 Stratigraphic Properties Table FD

[0076]

[0077]

[0078] In this embodiment, there are 42 stratigraphic sets with attributes in Stra. However, the Quaternary strata and igneous strata do not have meaningful dip data, so they are marked as None. There are 11 deep strata in Stra, which are generally considered to be indirectly obtained. Considering that the method takes into account multiple profiles, their location in the geological structure is not considered for the time being.

[0079] (1-3) Load the fault vector layer and store it in the fault set Faults = {f b |b=1,2,…,fn}, where f b Let fn represent the b-th fault, and fn represent the number of faults; in this embodiment, fn = 4.

[0080] (1-4) Load all sedimentary strata in the study area and store them in the sedimentary strata table Stable={st i |i=0,1,…,m}, where st i The stratum code represents the i-th sedimentary stratum from newest to oldest, and m is the total number of sedimentary strata in the current area. The study area in this embodiment belongs to the Ningzhen Mountains, which has a total of 61 sedimentary strata, m = 61 (Jiangsu Provincial Bureau of Geology and Mineral Resources. Geological Records of Ningzhen Mountains [M]. Jiangsu Science and Technology Press, 1989). The specific sedimentary stratum sequence is shown in Table 2.

[0081] Table 2 Correspondence between Sedimentary Stratigraphic Codes and Sub-codes SF

[0082]

[0083] (1-5) According to the stratigraphic sequence in Stable, each sedimentary stratum is stratigraphically encoded, that is, a single character is used to represent a sedimentary stratum, forming the stratigraphic encoding table SF.

[0084] (2) Generate a coding sequence DS consisting of the codes of all exposed strata in the stratigraphic set Stra, based on the stratigraphic coding table SF. This step specifically includes:

[0085] (2-1) Traverse the strata in the stratigraphic set Stra, check whether they are exposed on the surface, remove the unexposed strata from the stratigraphic set Stra and store them in the undetermined spatial location set St deep In this embodiment, there are 11 deep strata that were removed from the set Stra in this step.

[0086] (2-2) Rearrange the strata in the stratigraphic set Stra in left-to-right order, and simultaneously rearrange the contents of FD in the stratigraphic attribute table in this order; in this embodiment, the strata in the set Stra are designated as Q, J. 1-2 xn 1 ,T3f,T2h,Q,T2h,…,P1q,γΠδ 3 ,P1q,P1y,P2l,P2d,Q,K2c,Q;

[0087] (2-3) Based on the stratigraphic coding table SF, map each stratigraphic layer in Stra to a code, and form the coding sequence DS = {ds} in sequence. c |c=1,2,…,cn} where, ds ccn is the code for the c-th stratum, and cn is the number of strata. In this embodiment, cn = 42, and the DS sequence is {0,17,18,19,0,19,…,28,-1,28,27,25,24,0,6,0}.

[0088] (3) Traverse the coding sequence DS, store the events corresponding to the Quaternary strata into the geological event set GeoEvents in chronological order, store the non-Quaternary strata in the stratigraphic set Stra into the set Stra2, and delete the Quaternary strata from the coding sequence DS.

[0089] This step specifically includes:

[0090] (3-1) Read any element ds from the encoded sequence DS c Obtain the corresponding stratigraphic set Stra from the stratigraphic set Stra. a and stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ;

[0091] (3-2) If the stratigraphic code is dc a If the stratum begins with Q, then it is a Quaternary stratum, and proceed to step (3-3); otherwise, store the stratum in the set Stra2 which has been reduced from Quaternary strata, and return to step (3-1).

[0092] (3-3) Stra a With dt a Geological objects and their occurrence or occurrence time are stored in the geological event set GeoEvents; in this embodiment, there are 4 geological events corresponding to the Quaternary strata, and their positions in the DS sequence are 1, 5, 25, 40, and 42;

[0093] (3-4) Repeat steps (3-1)-(3-3) until all elements in the encoded sequence DS have been traversed;

[0094] (3-5) Delete all Quaternary strata codes in the coding sequence DS, and merge adjacent identical strata codes in the coding sequence DS after the removal of Quaternary strata. These strata were originally divided in two due to Quaternary strata coverage. In this embodiment, the strata on both sides of the Quaternary strata at position 5 in the sequence are both designated as T2h, and they are merged when the Quaternary strata are removed.

[0095] (4) Traverse the coding sequence DS, store the intrusion events corresponding to the intrusive strata into the geological event set GeoEvents in chronological order, store the non-intrusive rock strata in Stra2 into the stratum set Stra3, and delete the intrusive strata from the coding sequence DS.

[0096] This step specifically includes:

[0097] (4-1) Read an element ds from the encoded sequence DS c Get the corresponding stratigraphic stra in set Stra2 a and its stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ;

[0098] (4-2) If the stratigraphic code is dc a If the stratum begins with a Greek letter, it is an intrusive stratum, and step (4-3) is executed; otherwise, the stratum is stored in the set Stra3, which has eliminated intrusive strata, and step (4-1) is returned.

[0099] (4-3) Stra a With dt a Geological objects and their occurrence or generation time are stored in the geological event set GeoEvents; in this embodiment, there is an intrusive rock stratum, which is located at 35 in the DS sequence;

[0100] (4-4) Repeat steps (4-1)-(4-3) until all elements in all encoded sequences DS have been traversed;

[0101] (4-5) Delete all intrusive rock codes in the coding sequence DS, and merge adjacent consistent stratigraphic codes in the coding sequence DS that were generated after the removal of intrusive rock strata. These strata were originally divided in two by the intrusion of igneous rocks; in this embodiment, the code is γΠδ. 3 Both sides of the stratum are P1q.

[0102] (5) Traverse the fault set Faults, obtain the adjacent strata on both sides of each fault in Stra3, determine the type of fault and fault event based on the strata on both sides of the fault, and then store them in the geological event set GeoEvents, and add or delete the strata codes affected by the fault in the coding sequence DS accordingly.

[0103] like Figure 2 As shown, this step specifically includes:

[0104] (5-1) Read a fault f in the fault set Faults in geological age order from newest to oldest. b And calculate the dip angle of the fault. f In this embodiment, when the second fault is f b At that time, the tilt angle is 50°;

[0105] (5-2) Find f in set Stra3a The adjacent left and right strata have exactly two exposed strata along the fault line in the cross section, located on both sides of the fault, denoted as stra. l ,stra r In this embodiment, stra l The code is T1x 1 The tilt angle is 79°, stra r The code is P2l with an inclination angle of 62°;

[0106] (5-3) Calculate stra l ,stra r The average angle of inclination ave Compare angles ave with angle f If the values ​​are equal or the difference is within a preset range, it indicates that the fault is likely developed at the interface between the two strata, and step (5-4) is executed; otherwise, step (5-5) is executed. If the dip angles of the fault and the adjacent strata are close, it indicates that the fault develops along the interface between the strata or within the strata joints; otherwise, there may be missing or repeated strata. The threshold for the angle difference needs to be determined according to the study area. In this embodiment, the threshold is set to 45°. The second fault and T1x 1 The faults have opposite dips and large dip angles, therefore the fault is identified as a normal fault.

[0107] (5-4) If the fault dip angle is close to that of the adjacent stratum, it indicates that the fault is located at the interface between the two strata or within a certain stratum. However, the cross-section provides limited information, making it impossible to determine whether the fault is a normal or reverse fault. If stra l ,stra r If the lithology is consistent, then the fault f b Mark the fault as developed within the strata, otherwise mark it as developed on the stratum contact surface, and store the fault in the GeoEvents set, then execute steps (5-7);

[0108] (5-5) If the dip angle of the fault differs significantly from that of the adjacent strata, the fault may be a normal fault or a reverse fault. Repeatedly search for continuous layers in the adjacent strata on both sides of the fault. If found, its type is classified as a reverse fault and stored in the GeoEvents set, then proceed to step (5-7); otherwise, its type is classified as a normal fault and proceed to step (5-6). This step specifically includes: (5-5-1) Obtaining the adjacent strata on both sides of the fault and storing them in an orderly manner in the Set Stra Near ={stra j |j=1,2,…,k}, where stra j Let j be the j-th neighboring stratum, and k be the total number of neighboring strata; (5-5-2) in the set Stra NearThe search for repeating stratigraphic subsequences requires that at least one subsequence is related to the current fault f. b Adjacency; (5-5-3) If a repeating subsequence is found, then the current fault f a The repeated stratigraphic subsequences are stored as reverse fault events in the GeoEvents set; steps (5-7) are executed; otherwise, the current fault f is... a The type is a normal fault; proceed with steps (5-6).

[0109] (5-6) Obtain stra l ,stra r The stratigraphic code is used, and the missing sedimentary strata between the two strata are found in the Stable sedimentary stratigraphic table. The missing sedimentary strata code is then added to the coding sequence DS. l ,stra r The stratigraphic codes are used to identify the fault and the missing sedimentary strata, and the fault and missing sedimentary strata are stored in the GeoEvents set. In this embodiment, the missing strata are designated as P2d. This step specifically includes: (5-6-1) searching for and stra in the Stable sedimentary stratigraphic table. l ,stra r Get their index number i l i r In this embodiment, i l i r 23 and 25 respectively; (5-6-2) Obtain the index of the Stable sedimentary stratigraphic table in i l i r The codes of the sedimentary strata between them are stored in the set DS. lr In this embodiment, DS lr The set is {P2d}; (5-6-3) the current fault f a The type is denoted as a normal fault, and is related to the set DS. lr They are stored together in the GeoEvents collection.

[0110] (5-7) Repeat steps (5-1)-(5-6) until all faults have been traversed;

[0111] (5-8) Iterate through the encoded sequence DS sequentially, delete consecutively repeated stratigraphic subsequences, and keep only the unique one among the repeated subsequences.

[0112] (6) Based on the coding sequence DS, the symmetrical repeating fold strata are searched and stored in the geological event set GeoEvents according to the principle of symmetry and repetition.

[0113] This step specifically includes:

[0114] (6-1) Obtain a stratigraphic code ds from the coding sequence DS in sequence. c ;

[0115] (6-2) If ds c ds c+1 If they are not the same stratum, then use ds c Centered on the coding sequence DS, the codes of the strata on both sides are compared one by one until the first inconsistent code is encountered. The coding sequence of the compared strata is a palindrome. The strata of the palindrome are stored as folds in the fold set Folds.

[0116] (6-3) If ds c ds c+1 If they are the same stratum, then use ds c ds c+1 Centered on the coding sequence DS, the stratigraphic codes on both sides are compared one by one until the first inconsistent stratigraphic code is encountered. The comparison stops. The sequence of stratigraphic codes that have been compared is a palindrome. The stratigraphic codes of the palindrome are stored as folds in the fold set Folds.

[0117] (6-4) Repeat steps (6-1)-(6-3) until all palindromes are detected; in this embodiment, the detected palindrome sequences are {T2x,T1s,T2x}, {T1s,T2x,T1s}, {D3w 1 ,S3m,S2f 2 S2f 1 ,S1g,S2f 1 S2f 2 ,S3m,D3w 1};

[0118] (6-5) Traverse all folds in the fold set Folds, and determine the type of fold corresponding to the palindrome string as syncline or anticline based on the arrangement of the strata in the folds; in this embodiment, {T2x,T1s,T2x} is an anticline, {T1s,T2x,T1s} is a syncline, and {D3w 1 ,S3m,S2f 2 S2f 1 ,S1g,S2f 1 S2f 2 ,S3m,D3w 1} is an anticline; the specific methods for determining this include: (6-5-1) obtaining a fold from the fold set Folds. i This fold contains j strata. In this embodiment, for example, fold i For {D3w 1 ,S3m,S2f 2 S2f 1,S1g,S2f 1 S2f 2 ,S3m,D3w 1}, j = 9; (6-5-2) Get the fold i The geological age of the first stratum in the middle is dt1, and the geological age of the second stratum is dt1. Geological age dt of the strata at the location mid In this embodiment, dt1 is D3w 1 dt mid Let S1g be the value; (6-5-3) Compare dt1 and dt mid The geological time sequence, if dt1 is older, then fold i For synclinal; if dt mid Older ones fold i Anticline; record the anticline / syncline type in fold. i In the object; for example, in this embodiment, dt1 is D3w 1 dt mid For S1g,dt mid Older, so the current fold i It is an anticline; (6-5-4) Repeat steps (6-5-1)-(6-5-3) until all folds have been traversed.

[0119] (6-6) Store all folds in the fold set Folds into the geological event set GeoEvents.

[0120] (7) Format and output the GeoEvents collection. This step specifically includes:

[0121] (7-1) Sort all geological events in GeoEvents in chronological order;

[0122] (7-2) Format and output the geological event objects of different types in GeoEvents as local files. In this embodiment, the list obtained after formatting the GeoEvents collection is shown in Table 3. Mark the geological events on the corresponding profiles:

[0123] Table 3 Geological Events Table

[0124]

[0125] Example 2

[0126] This invention provides a structure for an automatic geological event analysis device 12 based on geological profiles, which provides services for implementing the method of Embodiment 1 above. An exemplary device 12 suitable for implementing the embodiments of this invention is shown. The shown device 12 is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this invention.

[0127] Device 12 is presented in the form of a general-purpose computing device. Components of device 12 may include, but are not limited to: one or more processors or processing units 16, system memory 28, and a bus 18 connecting different system components, including system memory 28 and processing unit 16.

[0128] Bus 18 represents one or more of several bus architectures, including a memory bus or memory controller, a peripheral bus, a graphics acceleration port, a processor, or a local bus using any of the various bus architectures. For example, these architectures include, but are not limited to, the Industry Standard Architecture (ISA) bus, the Micro Channel Architecture (MAC) bus, the Enhanced ISA bus, the Video Electronics Standards Association (VESA) local bus, and the Peripheral Component Interconnect (PCI) bus.

[0129] Device 12 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by device 12, including volatile and non-volatile media, removable and non-removable media.

[0130] System memory 28 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 30 and / or cache memory 32. Device 12 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, storage system 34 may be used to read and write non-removable, non-volatile magnetic media (commonly referred to as "hard disk drives"). Disk drives for reading and writing to removable non-volatile disks (e.g., "floppy disks") and optical disk drives for reading and writing to removable non-volatile optical disks (e.g., CD-ROMs, DVD-ROMs, or other optical media) may be provided. In these cases, each drive may be connected to bus 18 via one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of the embodiments of the present invention.

[0131] A program / utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28. Such program modules 42 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include an implementation of a network environment. Program modules 42 typically perform the functions and / or methods described in the embodiments of the present invention.

[0132] Device 12 can also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), and with one or more devices that enable a user to interact with device 12, and / or with any device that enables device 12 to communicate with one or more other computing devices (e.g., network card, modem, etc.). This communication can be performed via input / output (I / O) interface 22. Furthermore, device 12 can also communicate with one or more networks (e.g., local area network (LAN), wide area network (WAN), and / or public networks, such as the Internet) via network adapter 20. Network adapter 20 communicates with other modules of device 12 via bus 18. It should be understood that, although not shown in the figures, other hardware and / or software modules can be used in conjunction with device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0133] The processing unit 16 executes various functional applications and data processing by running programs stored in the system memory 28, such as implementing the method provided in Embodiment 1 of the present invention.

Claims

1. An automatic analysis method for geological events based on geological profiles, characterized in that... Includes the following steps: (1) Load the geological profile vector layer, read the strata and store them in the strata set Stra, and at the same time obtain the strata attribute information table FD; load the fault vector layer and store it in the fault set Faults; All sedimentary strata in the study area were loaded and stored in the sedimentary stratigraphic table Stable, and each sedimentary strata was coded to form the stratigraphic coding table SF. (2) Generate a coding sequence DS consisting of the codes of all exposed strata in the stratigraphic set Stra based on the stratigraphic coding table SF; (3) Traverse the coding sequence DS, store the events corresponding to the Quaternary strata into the geological event set GeoEvents in chronological order, store the non-Quaternary strata in the stratigraphic set Stra into the set Stra2, and delete the Quaternary strata from the coding sequence DS. (4) Traverse the coding sequence DS, store the intrusion events corresponding to the intrusion strata into the geological event set GeoEvents in chronological order, store the non-intrusion rock strata in Stra2 into the stratum set Stra3, and delete the intrusion strata from the coding sequence DS. (5) Traverse the fault set Faults, obtain the adjacent strata on both sides of each fault in Stra3, determine the type of fault and fault event based on the strata on both sides of the fault, and then store it in the geological event set GeoEvents, and add or delete the strata codes affected by the fault in the coding sequence DS accordingly. (6) Based on the coding sequence DS, the symmetrical repeating fold strata are searched and stored in the geological event set GeoEvents according to the principle of symmetric repetition; (7) Output the GeoEvents collection in a formatted manner; Step (6) specifically includes: (6-1) Obtain a stratigraphic code ds from the coding sequence DS in sequence. c ; (6-2) If ds c ds c+1 If they are not the same stratum, then use ds c Centered on the coding sequence DS, the codes of the strata on both sides are compared one by one until the first inconsistent code is encountered. The coding sequence of the compared strata is a palindrome. The strata of the palindrome are stored as folds in the fold set Folds. (6-3) If ds c ds c+1 If they are the same stratum, then use ds c ds c+1 Centered on the coding sequence DS, the stratigraphic codes on both sides are compared one by one until the first inconsistent stratigraphic code is encountered. The comparison stops. The sequence of stratigraphic codes that have been compared is a palindrome. The stratigraphic codes of the palindrome are stored as folds in the fold set Folds. (6-4) Repeat steps (6-1)-(6-3) until all palindromes are detected; (6-5) Traverse all folds in the fold set Folds, and determine whether the palindrome string corresponds to a syncline or anticline based on the arrangement of the old and new strata in the fold. (6-6) Store all folds in the fold set Folds into the geological event set GeoEvents.

2. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (1) specifically includes: (1-1) Load the geological profile vector layer and store it in the stratigraphic set Stra={stra a |a=1, 2, …, sn}, where stra a Let 'a' represent the 'a'-th stratum and 'sn' represent the stratum number; sequentially read the stratigraphic codes of all strata in Stra and store them in the stratigraphic code set DC={dc a |a=1, 2, …, sn} where dc a This represents the stratigraphic designation of the a-th stratum. (1-2) Obtain the dip angle and geological age of each stratum in the stratigraphic set Stra, and store them in the stratigraphic attribute table FD; (1-3) Load the fault vector layer and store it in the fault set Faults={f b | b=1, 2, …, fn}, where f b Let f_b represent the b-th fault, and f_n represent the number of faults. (1-4) Load all sedimentary strata in the study area and store them in the sedimentary strata table Stable={st i |i=0, 1, …, m}, where st i This represents the stratigraphic code of the i-th sedimentary stratum from newest to oldest, where m is the total number of all sedimentary strata in the current region. (1-5) According to the stratigraphic sequence in Stable, each sedimentary stratum is stratigraphically encoded, that is, a single character is used to represent a sedimentary stratum, forming the stratigraphic encoding table SF.

3. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (2) specifically includes: (2-1) Traverse the strata in the stratigraphic set Stra and check whether they are exposed on the surface. Remove the strata that are not exposed from the stratigraphic set Stra. (2-2) Rearrange the strata in the stratigraphic set Stra in order from left to right, and at the same time rearrange the contents of FD in the stratigraphic attribute table in this order; (2-3) Based on the stratigraphic coding table SF, map each stratigraphic layer in Stra to a code, and form the coding sequence DS={ds} in sequence. c |c=1, 2, …, cn} where ds c It is the code of the c-th stratum, and cn is the number of strata.

4. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (3) specifically includes: (3-1) Read any element ds from the encoded sequence DS c Obtain the corresponding stratigraphic set Stra from the stratigraphic set Stra. a and stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ; (3-2) If the stratigraphic code is dc a If the stratum begins with Q, then it is a Quaternary stratum, and proceed to step (3-3); otherwise, store the stratum in the set Stra2 which has been reduced from Quaternary strata, and return to step (3-1). (3-3) Stra a With dt a Geological objects and their occurrence or generation time are stored in the geological event set GeoEvents; (3-4) Repeat steps (3-1)-(3-3) until all elements in the encoded sequence DS have been traversed; (3-5) Delete all Quaternary strata codes in the coding sequence DS and merge adjacent consistent strata codes in the coding sequence DS after the Quaternary strata were removed.

5. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (4) specifically includes: (4-1) Read an element ds from the encoded sequence DS c Get the corresponding stratigraphic stra in set Stra2 a and its stratigraphic designation DC a And obtain the corresponding geological age dt from the stratigraphic attribute information table FD. a ; (4-2) If the stratigraphic code is dc a If the stratum begins with a Greek letter, it is an intrusive stratum, and step (4-3) is executed; otherwise, the stratum is stored in the set Stra3, which has eliminated intrusive strata, and step (4-1) is returned. (4-3) Stra a With dt a Geological objects and their occurrence or generation time are stored in the geological event set GeoEvents; (4-4) Repeat steps (4-1)-(4-3) until all elements in all encoded sequences DS have been traversed; (4-5) Delete all intrusive rock codes in the coding sequence DS and merge adjacent consistent stratigraphic codes in the coding sequence DS after the intrusive rock strata were eliminated.

6. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (5) specifically includes: (5-1) Read a fault f in the fault set Faults in geological age order from newest to oldest. b And calculate the dip angle of the fault. f ; (5-2) Find f in set Stra3 a The adjacent left and right strata are denoted as stra l ,stra r ; (5-3) Calculate stra l ,stra r The average angle of inclination ave Compare angles ave with angle f If the values ​​are equal or the difference is within a preset range, proceed to step (5-4); otherwise, proceed to step (5-5). (5-4) If stra l ,stra r If the lithology is consistent, then the fault f b Mark the fault as developed within the strata, otherwise mark it as developed on the stratum contact surface, and store the fault in the GeoEvents set, then execute steps (5-7); (5-5) Repeatedly search for continuous layers in the adjacent strata on both sides of the fault. If they are found, the type is defined as reverse fault and stored in the GeoEvents set. Then proceed to step (5-7). Otherwise, the type is defined as normal fault and proceed to step (5-6). (5-6) Obtain stra l ,stra r The stratigraphic code is used, and the missing sedimentary strata between the two strata are found in the Stable sedimentary stratigraphic table. The missing sedimentary strata code is then added to the coding sequence DS. l ,stra r Between the stratigraphic codes, the fault and the missing sedimentary strata are stored in the GeoEvents set; (5-7) Repeat steps (5-1)-(5-6) until all faults have been traversed; (5-8) Iterate through the encoded sequence DS sequentially, delete consecutively repeated stratigraphic subsequences, and keep only the unique one among the repeated subsequences.

7. The automatic analysis method for geological events based on geological profiles according to claim 6, characterized in that: Step (5-5) specifically includes: (5-5-1) Obtain the adjacent strata on both sides of the fault and store them in an ordered manner in the set Stra Near ={stra j |j=1, 2,…, k}, where stra j Let j be the j-th nearest stratum, and k be the total number of nearest strata; (5-5-2) In the set Stra Near The search for repeating stratigraphic subsequences requires that at least one subsequence is related to the current fault f. b Adjacency; (5-5-3) If a repeating subsequence is found, then the current fault f a The repeated stratigraphic subsequences are stored as reverse fault events in the GeoEvents set; steps (5-7) are executed; otherwise, the current fault f is... a The type is a normal fault; proceed with steps (5-6). Steps (5-6) specifically include: (5-6-1) Locate and stra in the Stable sedimentary stratigraphic table l ,stra r Get their index number i l i r ; (5-6-2) Obtaining the index of the sedimentary stratigraphic table within the Stable in i l i r The codes of the sedimentary strata between them are stored in the set DS. lr ; (5-6-3) Set the current fault f a The type is denoted as a normal fault, and is related to the set DS. lr They are stored together in the GeoEvents collection.

8. The automatic analysis method for geological events based on geological profiles according to claim 1, characterized in that: Step (6-5) specifically includes: (6-5-1) Get a fold from the fold set Folds i This fold contains j strata; (6-5-2) Obtaining folds i The geological age dt1 of the first stratum and the geological age dt of the stratum at position ⌊(j+1) / 2⌋ mid ; (6-5-3) Compare dt1 and dt mid The geological time sequence, if dt1 is older, then fold i For synclinal; if dt mid Older ones fold i Anticline; record the anticline / syncline type in fold. i In the object; (6-5-4) Repeat steps (6-5-1)-(6-5-3) until all wrinkles have been traversed.

9. An automatic geological event analysis device based on geological profiles, 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, it implements the method described in any one of claims 1-8.

Images