CN112231423B - Geological map-based fault structure automatic recovery method and device - Google Patents
Geological map-based fault structure automatic recovery method and device Download PDFInfo
- Publication number
- CN112231423B CN112231423B CN202011000473.3A CN202011000473A CN112231423B CN 112231423 B CN112231423 B CN 112231423B CN 202011000473 A CN202011000473 A CN 202011000473A CN 112231423 B CN112231423 B CN 112231423B
- Authority
- CN
- China
- Prior art keywords
- fault
- rpl
- lpl
- stratum
- boundary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000011084 recovery Methods 0.000 title claims abstract description 24
- 230000011218 segmentation Effects 0.000 claims description 45
- 230000015572 biosynthetic process Effects 0.000 claims description 29
- 238000006073 displacement reaction Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims 1
- 238000002591 computed tomography Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 21
- 239000011435 rock Substances 0.000 description 3
- 238000012876 topography Methods 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/25—Integrating or interfacing systems involving database management systems
- G06F16/252—Integrating or interfacing systems involving database management systems between a Database Management System and a front-end application
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a fault structure automatic recovery method and device based on a geological map, wherein the method comprises the following steps: firstly, reading stratum data and fault data; secondly, stratum and fault sections on two sides of the fault are obtained based on buffer zone neighborhood analysis; then, restoring stratum at two sides of the fault based on stratum corresponding relation; and finally, restoring the fault sections at two sides of the fault based on the fault section corresponding relation, and generating a fault restoration process geological map. The invention realizes an automatic fault structure recovery method based on a geological map. Compared with the prior art, the invention has high efficiency, better fault structure recovery effect and higher automation degree.
Description
Technical Field
The invention relates to the field of application of geographic information technology, in particular to a fault structure automatic recovery method and device based on a geological map.
Background
Recovery of geologic formations refers to the process of extrapolating the original topography of the formation based on the formation topography type plus other conditions as constraints. Faults are one of the most important structures of the crust, and are formed by remarkable displacement of stratum or rock mass along fracture surfaces, so that cracks and scarps are often formed. The method aims at scientific inference of fault evolution process, restores the original shape of rock stratum when fracture does not occur, and has certain research significance and application value for geological disaster prediction, geological resource development, geological environment protection and the like.
The fault recovery application is mainly embodied in the following three aspects: 1) Fault recovery based on geologic maps. Judging the consistency of stratum contained in the tray body in the fault, and then merging the displaced strata to finally form an original geological map after removing the fault line; 2) Fault recovery based on geologic profile. The discrete points of the geological layer can be restored to the original shape of the rock stratum before deformation according to the reverse sequence of fracture of the fault; 3) Fault recovery based on three-dimensional models. And mainly modifying attribute data inside the fault entity to restore the deformed fault structure into a complete geological entity. At present, related applications are related, but related work is mainly completed manually by geological specialists, so that the investment is large, the efficiency is poor, and the research of an automatic recovery method of faults is urgently needed.
The geological map data has the characteristics of rich data, wide sources, wide coverage, less use limit and the like. Starting from the geological map, fault structure recovery research is carried out, and the method has good feasibility and practicality.
Disclosure of Invention
The invention aims to: aiming at the problems existing in the prior art, the invention provides a high-efficiency fault structure automatic recovery method and device based on a geological map.
The technical scheme is as follows: the fault structure automatic recovery method based on the geological map comprises the following steps:
(1) Reading stratum data and fault line data to generate a stratum set P, a fault segmentation set F and a fault segmentation time sequence set T;
(2) Reading a fault section F from the fault section set F according to the fault section timing i ；
(3) Acquiring fault segments f i Left stratum number set LP, left fault section number set LF, right stratum number set RP, right fault section number set RF;
(4) Restoring the boundary line of each stratum in the left stratum number set LP and the right stratum number set RP based on the stratum corresponding relation, and storing the boundary line into the left boundary line set LPL and the right boundary line set RPL respectively;
(5) Restoring fault segmentation f based on the left side layer boundary set LPL and the right side layer boundary set RPL i A fault cutting boundary line of the stratum at two sides, and recovering the stratum based on the fault cutting boundary line and the layer boundary line;
(6) Restoring fault section f based on fault correspondence i Fault segmentation on two sides;
(7) Deriving a geological map based on the recovered strata and faults;
(8) And (3) circularly executing the steps (2) - (7) until the fault segmentation set F is traversed, and recovering all fault structures.
Further, step (1) includes:
(1-1) reading tomographic line vector data to a tomographic segment set f= { F i I = 1,2, …, fn; wherein f i Representing an i-th fault section, fn representing the number of fault sections;
(1-2) acquiring the number of the fault section, and storing the number of the fault section into a fault section number set fid= { FID i I = 1,2, …, fn; wherein, fid i For fault segmentation f i Is the number of (2);
(1-3) acquiring time sequence information of fault segmentation, and storing the time sequence information into a fault segmentation time sequence set T= { T i I = 1,2, …, fn; wherein t is i For fault segmentation f i Is a timing sequence of (2);
(1-4) reading formation data into a formation set p= { P j |j=1, 2, …, pn }; wherein p is j Represents the jth stratum, pn represents the number of strata.
Further, the step (3) includes:
(3-1) buffer zone computing interface based on GIS software, according to preset width W, computing fault section f i Buffer b of (2) i ；
(3-2) creation of a fault segment f i A left stratum number set LP and a right stratum number set RP of (a);
(3-3) reading any one formation P from the formation set P j When p is j And buffer b i When the left boundary of (1) intersects, the stratum number j is stored in the set LP, and when p is j And buffer b i When the right boundary of the (a) is intersected, storing the stratum number j into a set RP;
(3-4) circularly executing the step (3-3) until the stratum set P is traversed to obtain a left stratum number set LP and a right stratum number set RP;
(3-5) creation of a fault segment f i Left side fault section number set LF and right side fault section number set RF;
(3-6) reading any one from the tomosynthesis segment set FFault section f k When f k And buffer b i When the left boundary of (a) is intersected, the fault segment number k is stored in the set LF, and when f k And buffer b i When the right boundary of the fault section is intersected, storing the fault section number k into a set RF;
(3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed to obtain a left fault segmentation number set LF and a right fault segmentation number set RF.
Further, the step (4) includes:
(4-1) calculating the fault section f according to i Is of the displacement vector of (2)
Wherein dir i Representing fault segmentation f i Trend of d i Representing fault segmentation f i Displacement of (2);
(4-2) extracting the boundary line of each stratum in the left stratum number set LP, and storing into the left boundary line set LPL= { LPL u (us, ls, flag) |u=1, 2, …, lpn +1}; therein, lpl u Represents the u-th deck boundary, us being deck boundary lpl u Of the upper stratum type of lpl u In the absence of an upper formation, us takes the form of-1, ls as the bedding line lpl u A type of formation on the underside of (a); when lpl u In the absence of an underlying stratum, ls is-1, and the flag is lpl u The initial value is 0, lpn is the number of elements in LP;
(4-3) extracting the boundary line of each stratum in the right stratum number set RP, and storing into the right boundary line set RPL= { RPL v (us, ls, flag) |v=1, 2, …, rpn +1}; wherein rpl v Represents the v-th layer boundary, us is the layer boundary rpl v Upper stratum type of (1), when rpl v In the absence of an upper strata, us takes the form of-1, ls as the bedding boundary rpl v A type of formation on the underside of (a); when rpl v Without the lower partIn the case of side stratum, ls is taken as-1 and flag is rpl v The initial value is 0, and rpn is the number of elements in RP;
(4-4) matching fault segments f according to the sets LPL and RPL i The boundary line of the left stratum and the boundary line of the right stratum, and the boundary line is restored.
Further, the step (4-4) includes:
(4-4-1) reading the first deck boundary line LPL from the set LPL u Reading a first level boundary RPL from the set RPL v ；
(4-4-2) As lpl u (us) and rpl v (us) the same or lpl u (ls) and rpl v (ls) when the same, performing the step (4-4-3); otherwise, executing the step (4-4-6);
(4-4-3) at lpl u Is to start with the leftmost end point of the layer boundary lpl u The upper points are rearranged in sequence to obtain an arranged point set LPP u ＝{lpl u,g |g=1, 2, …, lgn }; in rpl v Starting from the rightmost end point of (2) and dividing the layer boundary rpl v The upper points are rearranged in sequence to obtain a point set RPP after arrangement v ＝{rpl v,h |h=1, 2, …, lhn }, wherein lpl u,g Representation lpl u Upper g-th point, lgn denotes lpl u Point number, rpl v,h Representation rpl v The upper h point, lhn, represents rpl v Counting points;
(4-4-4) reading LPP u Is lpl of the last point of (5) u,lgn And RPP v Last point rpl of (1) v,lhn Will lpl u,lgn And rpl v,lhn Is used as a reference point and the layer boundary line lpl is adjusted according to the following u 、rpl v Each point is arranged on the upper part;
(4-4-5) the method comprises the steps of lpl u (flag) and rpl v (flag) set to 1, and when u equals lpn +1 or v equals rpn +1, performing step (4-4-7); otherwise, the value of u is incremented by 1, the value of v is incremented by 1, and the step (4) is traced back-4-2)；
(4-4-6) traversing the layer boundary rpl in order v To rpl rpn+1 If rpl is present t Satisfy lpl u (us) and rpl t (us) the same or lpl u (ls) and rpl t (ls) the same, setting v to t, and tracing back to step (4-4-2); otherwise, the value of u is increased by 1, and the step (4-4-2) is traced back;
(4-4-7) setting u to 1 based on the displacement vectorAccording to the following adjustment lpl u And takes the adjusted layer boundary line as a recovered layer boundary line lpl u And will lpl u (flag) set to 2;
(4-4-8) setting v to 1 based on the displacement vectorAdjusting rpl according to v And takes the adjusted layer boundary line as a recovered layer boundary line rpl v And rpl is combined v (flag) set to 2;
(4-4-9) incrementing the value of u by 1 and the value of v by 1, looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is completed.
Further, step (5) includes:
(5-1) extracting the side boundary line of each stratum in the left stratum number set LP and the right stratum number set RP, and adding the side boundary line to the stratum side boundary line set SL;
(5-2) restoring the fault cut boundary based on the left side layer boundary set LPL and the right side layer boundary set RPL, storing the fault cut boundary into the cut boundary set CL;
(5-3) constructing a fault segment f based on the left-side layer boundary line set LPL and the right-side layer boundary line set RPL, the cutting boundary line set CL, and the formation side boundary line set SL i Is stored in a stratum layer set TP;
(5-4) for each floor level TP in the set TP, obtaining an upper level boundary of TP based on the left side level boundary set LPL and the right side level boundary set RPL, and writing a lower side stratum type of the upper level boundary into the TP stratum type attribute; adding tp to the stratum set P;
(5-5) removing the fault section f from the stratigraphic set P according to the left stratigraphic number set LP and the right stratigraphic number set RP i Is formed on both sides of the substrate.
Further, the step (5-2) includes:
(5-2-1) two layer boundary lines LPL in the left layer boundary line set LPL are successively read u 、lpl u+1 When lpl u (flag) =2 or lpl u+1 When (flag) =2, connect lpl u Rightmost end point and lpl u+1 Forming a connecting line at the rightmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-2) continuously reading two boundary lines RPL in the right boundary line set RPL v 、rpl v+1 When rpl v (flag) =2 or rpl v+1 When (flag) =2, rpl is connected v Leftmost endpoint and rpl v+1 Forming a connecting line at the leftmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-3) steps (5-2-1) to (5-2-2) are cyclically performed until the LPL, RPL are traversed.
Further, step (6) includes:
(6-1) reading any unread tomographic section number m in the left tomographic section number set LF to correspond to the tomographic section f m Taking the leftmost end point of the fault section f as a starting point m The upper points are rearranged in sequence to form a point set FM= { f m,p |p＝1,2,…,fpn}；
(6-2) reading any one of the Right fault section number sets RFUnread fault section number n to correspond to fault section f n Starting from the rightmost end point of (a), segment f n The upper points are rearranged to form a point set FN= { f n,q |q＝1,2,…,fqn}；
(6-3) if fault section f m Belonging fault number and fault section f n The fault numbers are the same, and the step (6-4) is executed; otherwise, executing the step (6-2);
(6-4) reading the last endpoint f of the Point set FM m,fpn The last endpoint f of the set of read points FN n,fqn At f m,fpn And f n,fqn Is used as a reference point, and each point of the point set FM and the point set FN is adjusted according to the following steps;
(6-5) generating the adjusted point set FM and the point set FN into the corresponding fault section f m 、f n Merging f m 、f n To one, add to the set of fault segments F, and remove F from the set F m 、f n ；
(6-6) removing the fault section number m from the set LF and the fault section number n from the set RF;
(6-7) circularly executing the steps (6-1) to (6-6) until the LF is traversed;
(6-8) reading any unread fault section number m in the set LF, and adjusting the fault section f according to the following m Each point is located above, and the adjusted fault section is taken as a recovered fault section f m ；
Wherein,representing fault segmentation f i Is a displacement vector of (a);
(6-9) read setAny unread fault section number n in RF, fault section f is adjusted according to the following n Each point is located above, and the adjusted fault section is taken as a recovered fault section f n ；
(6-10) performing steps (6-8) to (6-9) in a loop until the sets LF and RF are traversed;
(6-11) removing the fault section F from the fault section set F m 。
The automatic fault structure recovery device based on the geological map comprises a processor and a computer program which is stored in a memory and can run on the processor, wherein the processor realizes the method when executing the program.
The beneficial effects are that: compared with the prior art, the invention has the remarkable advantages that: the automatic degree is high, the efficiency is high, and the effect is good.
Drawings
FIG. 1 is geologic map data employed in the present embodiment;
fig. 2 is fault attribute data employed in the present embodiment;
FIG. 3 is a schematic flow chart of the present invention;
FIG. 4 is a schematic illustration of various boundaries of a formation;
FIG. 5 is a series of geologic maps generated during a fault recovery process.
Detailed Description
In the following, the technical scheme of the present invention is further described in detail, in this embodiment, geological map data (fig. 1 and 2) of the small Mao Shanzhi mountain area of south Beijing city is selected as experimental data, and a projection coordinate system adopted by the experimental data is a south Beijing 92 coordinate system. Further description will be provided by describing a specific embodiment with reference to the accompanying drawings.
As shown in fig. 3, the embodiment provides a fault structure automatic recovery method based on a geological map, which includes the following steps:
(1) And reading stratum data and fault line data to generate a stratum set P, a fault segmentation set F and a fault segmentation time sequence set T.
The method specifically comprises the following steps:
(1-1) reading tomographic line vector data to a tomographic segment set f= { F i I = 1,2, …, fn; wherein f i Representing an i-th fault section, fn representing the number of fault sections; in this embodiment, fn=6;
(1-2) acquiring the number of the fault section, and storing the number of the fault section into a fault section number set fid= { FID i I = 1,2, …, fn; wherein, fid i For fault segmentation f i Is the number of (2);
(1-3) acquiring time sequence information of fault segmentation, and storing the time sequence information into a fault segmentation time sequence set T= { T i I = 1,2, …, fn; wherein t is i For fault segmentation f i Is a timing sequence of (2);
(1-4) reading formation data into a formation set p= { P j |j=1, 2, …, pn }; wherein p is j The j-th layer is represented, pn represents the number of layers, and pn=45 in this example.
(2) Reading a fault section F from the fault section set F according to the fault section timing i 。
(3) Acquiring fault segments f i Left stratum number set LP, left fault section number set LF, right stratum number set RP, right fault section number set RF.
The method specifically comprises the following steps:
(3-1) buffer zone computing interface based on GIS software, according to preset width W, computing fault section f i Buffer b of (2) i The method comprises the steps of carrying out a first treatment on the surface of the W=20 meters in this example;
(3-2) creation of a fault segment f i A left stratum number set LP and a right stratum number set RP of (a);
(3-3) reading any one formation P from the formation set P j When p is j And buffer b i When the left boundary of (1) intersects, the stratum number j is stored in the set LP, and when p is j And buffer b i When the right boundary of the (a) is intersected, storing the stratum number j into a set RP;
(3-4) circularly executing the step (3-3) until the stratum set P is traversed to obtain a left stratum number set LP and a right stratum number set RP;
(3-5) creation of a fault segment f i Left side fault section number set LF and right side fault section number set RF;
(3-6) reading any one of the fault sections F from the fault section set F k When f k And buffer b i When the left boundary of (a) is intersected, the fault segment number k is stored in the set LF, and when f k And buffer b i When the right boundary of the fault section is intersected, storing the fault section number k into a set RF;
(3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed to obtain a left fault segmentation number set LF and a right fault segmentation number set RF.
(4) And restoring the boundary line of each stratum in the left stratum number set LP and the right stratum number set RP based on the stratum corresponding relation, and storing the boundary lines into the left boundary line set LPL and the right boundary line set RPL respectively.
The method comprises the following steps:
(4-1) calculating the fault section f according to i Is of the displacement vector of (2)
Wherein dir i Representing fault segmentation f i Trend of d i Representing fault segmentation f i Displacement of (2);
(4-2) extracting the boundary line of each stratum in the left stratum number set LP, and storing into the left boundary line set LPL= { LPL u (us, ls, flag) |u=1, 2, …, lpn +1}; therein, lpl u Represents the u-th deck boundary, us being deck boundary lpl u Of the upper stratum type of lpl u In the absence of an upper formation, us takes the form of-1, ls as the bedding line lpl u A type of formation on the underside of (a); when lpl u In the absence of an underlying stratum, ls is-1, and the flag is lpl u The initial value is 0, lpn is the number of elements in LP; in this embodiment, the layer boundary of the formation is shown in FIG. 4;
(4-3) extracting the boundary line of each stratum in the right stratum number set RP, and storing into the right boundary line set RPL= { RPL v (us, ls, flag) |v=1, 2, …, rpn +1}; wherein rpl v Represents the v-th layer boundary, us is the layer boundary rpl v Upper stratum type of (1), when rpl v In the absence of an upper strata, us takes the form of-1, ls as the bedding boundary rpl v A type of formation on the underside of (a); when rpl v In the absence of an underlying stratum, ls is taken to be-1 and flag is rpl v The initial value is 0, and rpn is the number of elements in RP;
(4-4) matching fault segments f according to the sets LPL and RPL i The boundary line of the left stratum and the boundary line of the right stratum, and the boundary line is restored. The method specifically comprises the following steps:
(4-4-1) reading the first deck boundary line LPL from the set LPL u Reading a first level boundary RPL from the set RPL v ；
(4-4-2) As lpl u (us) and rpl v (us) the same or lpl u (ls) and rpl v (ls) when the same, performing the step (4-4-3); otherwise, performing step (4-4-6), wherein (#) represents the attribute# -of the element;
(4-4-3) at lpl u Is to start with the leftmost end point of the layer boundary lpl u The upper points are rearranged in sequence to obtain an arranged point set LPP u ＝{lpl u,g |g=1, 2, …, lgn }; in rpl v Starting from the rightmost end point of (2) and dividing the layer boundary rpl v The upper points are rearranged in sequence to obtain a point set RPP after arrangement v ＝{rpl v,h |h=1, 2, …, lhn }, wherein lpl u,g Representation lpl u Upper g-th point, lgn denotes lpl u Point number, rpl v,h Representation rpl v The upper h point, lhn, represents rpl v Counting points;
(4-4-4) reading LPP u Is lpl of the last point of (5) u,lgn And RPP v Last point in (3)rpl v,lhn Will lpl u,lgn And rpl v,lhn Is used as a reference point and the layer boundary line lpl is adjusted according to the following u 、rpl v Each point is arranged on the upper part;
(4-4-5) the method comprises the steps of lpl u (flag) and rpl v (flag) set to 1, and when u equals lpn +1 or v equals rpn +1, performing step (4-4-7); otherwise, the value of u is increased by 1, the value of v is increased by 1, and the step (4-4-2) is traced back;
(4-4-6) traversing the layer boundary rpl in order v To rpl rpn+1 If rpl is present t Satisfy lpl u (us) and rpl t (us) the same or lpl u (ls) and rpl t (ls) the same, setting v to t, and tracing back to step (4-4-2); otherwise, the value of u is increased by 1, and the step (4-4-2) is traced back;
(4-4-7) setting u to 1 based on the displacement vectorAccording to the following adjustment lpl u And takes the adjusted layer boundary line as a recovered layer boundary line lpl u And will lpl u (flag) set to 2;
(4-4-8) setting v to 1 based on the displacement vectorAdjusting rpl according to v And takes the adjusted layer boundary line as a recovered layer boundary line rpl v And rpl is combined v (flag) set to 2;
(4-4-9) incrementing the value of u by 1 and the value of v by 1, looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is completed.
(5) Restoring fault segmentation f based on the left side layer boundary set LPL and the right side layer boundary set RPL i And recovering the stratum based on the fault cutting boundary line and the layer boundary line.
The method specifically comprises the following steps:
(5-1) extracting the side boundary line of each stratum in the left stratum number set LP and the right stratum number set RP, and adding the side boundary line to the stratum side boundary line set SL;
(5-2) restoring the fault cut boundary based on the left side layer boundary set LPL and the right side layer boundary set RPL, storing the fault cut boundary into the cut boundary set CL;
(5-3) constructing a fault segment f using GIS software based on the left side layer boundary set LPL and the right side layer boundary set RPL, the cutting boundary set CL, and the formation side boundary set SL i Is stored in a stratum layer set TP;
(5-4) for each floor level TP in the set TP, obtaining an upper level boundary of TP based on the left side level boundary set LPL and the right side level boundary set RPL, and writing a lower side stratum type of the upper level boundary into the TP stratum type attribute; adding tp to the stratum set P;
(5-5) removing the fault section f from the stratigraphic set P according to the left stratigraphic number set LP and the right stratigraphic number set RP i Is formed on both sides of the substrate.
Step (5-2) comprises:
(5-2-1) two layer boundary lines LPL in the left layer boundary line set LPL are successively read u 、lpl u+1 When lpl u (flag) =2 or lpl u+1 When (flag) =2, connect lpl u Rightmost end point and lpl u+1 Forming a connecting line at the rightmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-2) continuously reading two boundary lines RPL in the right boundary line set RPL v 、rpl v+1 When rpl v (flag) =2 or rpl v+1 When (flag) =2, rpl is connected v Leftmost endpoint and rpl v+1 Forming a connecting line at the leftmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-3) steps (5-2-1) to (5-2-2) are cyclically performed until the LPL, RPL are traversed.
(6) Restoring fault section f based on fault correspondence i The faults on two sides are segmented.
The method specifically comprises the following steps:
(6-1) reading any unread tomographic section number m in the left tomographic section number set LF to correspond to the tomographic section f m Taking the leftmost end point of the fault section f as a starting point m The upper points are rearranged in sequence to form a point set FM= { f m,p |p＝1,2,…,fpn}；
(6-2) reading any unread tomographic section number n in the right tomographic section number set RF to correspond to the tomographic section f n Starting from the rightmost end point of (a), segment f n The upper points are rearranged to form a point set FN= { f n,q |q＝1,2,…,fqn}；
(6-3) if fault section f m Belonging fault number and fault section f n The fault numbers are the same, and the step (6-4) is executed; otherwise, executing the step (6-2);
(6-4) reading the last endpoint f of the Point set FM m,fpn The last endpoint f of the set of read points FN n,fqn At f m,fpn And f n,fqn Is used as a reference point, and each point of the point set FM and the point set FN is adjusted according to the following steps;
(6-5) generating the adjusted point set FM and the point set FN into the corresponding fault section f m 、f n Merging f m 、f n To one, add to the set of fault segments F, and remove F from the set F m 、f n ；
(6-6) removing the fault section number m from the set LF and the fault section number n from the set RF;
(6-7) circularly executing the steps (6-1) to (6-6) until the LF is traversed;
(6-8) reading any unread fault section number m in the set LF, and adjusting the fault section f according to the following m Each point is located above, and the adjusted fault section is taken as a recovered fault section f m ；
Wherein,representing fault segmentation f i Is a displacement vector of (a);
(6-9) reading any unread fault section number n in the set RF, and adjusting the fault section f according to the following n Each point is located above, and the adjusted fault section is taken as a recovered fault section f n ；
(6-10) performing steps (6-8) to (6-9) in a loop until the sets LF and RF are traversed;
(6-11) removing the fault section F from the fault section set F m 。
(7) Based on the recovered strata and faults, a geologic map is derived.
(8) And (3) circularly executing the steps (2) - (7) until the fault segmentation set F is traversed, and recovering all fault structures. As shown in fig. 5, in this embodiment, 5 fault recovery process geologic maps are co-produced.
The embodiment also provides an automatic fault structure recovery device based on the geological map, which comprises a processor and a computer program stored on a memory and capable of running on the processor, wherein the processor realizes the method when executing the program.
In the embodiment of the invention, the buffer area is created only based on the interface provided by the ArcGIS software, and the method can also use APIs of QGIS, mapGIS and other software.
The above disclosure is only a preferred embodiment of the present invention and should not be construed as limiting the scope of the invention, which is defined by the appended claims.
Claims (7)
1. An automatic fault structure recovery method based on a geological map is characterized by comprising the following steps:
(1) Reading stratum data and fault line data to generate a stratum set P, a fault segmentation set F and a fault segmentation time sequence set T;
(2) Reading a fault section F from the fault section set F according to the fault section timing i ；
(3) Acquiring fault segments f i Left stratum number set LP, left fault section number set LF, right stratum number set RP, right fault section number set RF;
(4) Restoring the boundary line of each stratum in the left stratum number set LP and the right stratum number set RP based on the stratum corresponding relation, and storing the boundary line into the left boundary line set LPL and the right boundary line set RPL respectively;
(5) Restoring fault segmentation f based on the left side layer boundary set LPL and the right side layer boundary set RPL i A fault cutting boundary line of the stratum at two sides, and recovering the stratum based on the fault cutting boundary line and the layer boundary line;
(6) Restoring fault section f based on fault correspondence i Fault segmentation on two sides;
(7) Deriving a geological map based on the recovered strata and faults;
(8) Circularly executing the steps (2) - (7) until the fault segmentation set F is traversed, and recovering all fault structures;
the step (3) specifically comprises:
(3-1) buffer area calculation interface based on GIS software according to preset widthW, computed tomography segmentation f i Buffer b of (2) i ；
(3-2) creation of a fault segment f i A left stratum number set LP and a right stratum number set RP of (a);
(3-3) reading any one formation P from the formation set P j When p is j And buffer b i When the left boundary of (1) intersects, the stratum number j is stored in the set LP, and when p is j And buffer b i When the right boundary of the (a) is intersected, storing the stratum number j into a set RP;
(3-4) circularly executing the step (3-3) until the stratum set P is traversed to obtain a left stratum number set LP and a right stratum number set RP;
(3-5) creation of a fault segment f i Left side fault section number set LF and right side fault section number set RF;
(3-6) reading any one of the fault sections F from the fault section set F k When f k And buffer b i When the left boundary of (a) is intersected, the fault segment number k is stored in the set LF, and when f k And buffer b i When the right boundary of the fault section is intersected, storing the fault section number k into a set RF;
(3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed to obtain a left fault segmentation number set LF and a right fault segmentation number set RF;
the step (4) comprises:
(4-1) calculating the fault section f according to i Is of the displacement vector of (2)
Wherein dir i Representing fault segmentation f i Trend of d i Representing fault segmentation f i Displacement of (2);
(4-2) extracting the layer boundary of each layer in the left layer number set LPThe line is stored in the left layer boundary line set lpl= { LPL u (us, ls, flag) |u=1, 2, …, lpn +1}; therein, lpl u Represents the u-th deck boundary, us being deck boundary lpl u Of the upper stratum type of lpl u In the absence of an upper formation, us takes the form of-1, ls as the bedding line lpl u A type of formation on the underside of (a); when lpl u In the absence of an underlying stratum, ls is-1, and the flag is lpl u The initial value is 0, lpn is the number of elements in LP;
(4-3) extracting the boundary line of each stratum in the right stratum number set RP, and storing into the right boundary line set RPL= { RPL v (us, ls, flag) |v=1, 2, …, rpn +1}; wherein rpl v Represents the v-th layer boundary, us is the layer boundary rpl v Upper stratum type of (1), when rpl v In the absence of an upper strata, us takes the form of-1, ls as the bedding boundary rpl v A type of formation on the underside of (a); when rpl v In the absence of an underlying stratum, ls is taken to be-1 and flag is rpl v The initial value is 0, and rpn is the number of elements in RP;
(4-4) matching fault segments f according to the sets LPL and RPL i The boundary line of the left stratum and the boundary line of the right stratum, and the boundary line is restored.
2. The geological map-based fault structure automatic restoration method of claim 1, wherein: the step (1) comprises:
(1-1) reading tomographic line vector data to a tomographic segment set f= { F i I = 1,2, …, fn; wherein f i Representing an i-th fault section, fn representing the number of fault sections;
(1-2) acquiring the number of the fault section, and storing the number of the fault section into a fault section number set fid= { FID i I = 1,2, …, fn; wherein, fid i For fault segmentation f i Is the number of (2);
(1-3) acquiring time sequence information of fault segmentation, and storing the time sequence information into a fault segmentation time sequence set T= { T i I = 1,2, …, fn; wherein t is i For fault segmentation f i Is a timing sequence of (2);
(1-4) reading formation data to a formation set p= {p j |j=1, 2, …, pn }; wherein p is j Represents the jth stratum, pn represents the number of strata.
3. The geological map-based fault structure automatic restoration method of claim 1, wherein: the step (4-4) comprises:
(4-4-1) reading the first deck boundary line LPL from the set LPL u Reading a first level boundary RPL from the set RPL v ；
(4-4-2) As lpl u (us) and rpl v (us) the same or lpl u (ls) and rpl v (ls) when the same, performing the step (4-4-3); otherwise, executing the step (4-4-6);
(4-4-3) at lpl u Is to start with the leftmost end point of the layer boundary lpl u The upper points are rearranged in sequence to obtain an arranged point set LPP u ＝{lpl u,g |g=1, 2, …, lgn }; in rpl v Starting from the rightmost end point of (2) and dividing the layer boundary rpl v The upper points are rearranged in sequence to obtain a point set RPP after arrangement v ＝{rpl v,h |h=1, 2, …, lhn }, wherein lpl u,g Representation lpl u Upper g-th point, lgn denotes lpl u Point number, rpl v,h Representation rpl v The upper h point, lhn, represents rpl v Counting points;
(4-4-4) reading LPP u Is lpl of the last point of (5) u,lgn And RPP v Last point rpl of (1) v,lhn Will lpl u,lgn And rpl v,lhn Is used as a reference point and the layer boundary line lpl is adjusted according to the following u 、rpl v Each point is arranged on the upper part;
(4-4-5) the method comprises the steps of lpl u (flag) and rpl v (flag) set to 1, and when u equals lpn +1 or v equals rpn +1, performing step (4-4-7); otherwise, the value of u is increased by 1, the value of v is increased by 1, and the step (4-4-2) is traced back;
(4-4-6) traversing the layer boundary rpl in sequence v To rpl rpn+1 If rpl is present t Satisfy lpl u (us) and rpl t (us) the same or lpl u (ls) and rpl t (ls) the same, setting v to t, and tracing back to step (4-4-2); otherwise, the value of u is increased by 1, and the step (4-4-2) is traced back;
(4-4-7) setting u to 1 based on the displacement vectorAccording to the following adjustment lpl u And takes the adjusted layer boundary line as a recovered layer boundary line lpl u And will lpl u (flag) set to 2;
(4-4-8) setting v to 1 based on the displacement vectorAdjusting rpl according to v And takes the adjusted layer boundary line as a recovered layer boundary line rpl v And rpl is combined v (flag) set to 2;
(4-4-9) incrementing the value of u by 1 and the value of v by 1, looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is completed.
4. The geological map-based fault structure automatic restoration method of claim 1, wherein: the step (5) comprises:
(5-1) extracting the side boundary line of each stratum in the left stratum number set LP and the right stratum number set RP, and adding the side boundary line to the stratum side boundary line set SL;
(5-2) restoring the fault cut boundary based on the left side layer boundary set LPL and the right side layer boundary set RPL, storing the fault cut boundary into the cut boundary set CL;
(5-3) constructing a fault segment f based on the left-side layer boundary line set LPL and the right-side layer boundary line set RPL, the cutting boundary line set CL, and the formation side boundary line set SL i Is stored in a stratum layer set TP;
(5-4) for each floor level TP in the set TP, obtaining an upper level boundary of TP based on the left side level boundary set LPL and the right side level boundary set RPL, and writing a lower side stratum type of the upper level boundary into the TP stratum type attribute; adding tp to the stratum set P;
(5-5) removing the fault section f from the stratigraphic set P according to the left stratigraphic number set LP and the right stratigraphic number set RP i Is formed on both sides of the substrate.
5. The geological map based fault structure automatic restoration method of claim 4, wherein: step (5-2) comprises:
(5-2-1) two layer boundary lines LPL in the left layer boundary line set LPL are successively read u 、lpl u+1 When lpl u (flag) =2 or lpl u+1 When (flag) =2, connect lpl u Rightmost end point and lpl u+1 Forming a connecting line at the rightmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-2) continuously reading two boundary lines RPL in the right boundary line set RPL v 、rpl v+1 When rpl v (flag) =2 or rpl v+1 When (flag) =2, rpl is connected v Leftmost endpoint and rpl v+1 Forming a connecting line at the leftmost end point, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
(5-2-3) steps (5-2-1) to (5-2-2) are cyclically performed until the LPL, RPL are traversed.
6. The geological map-based fault structure automatic restoration method of claim 1, wherein: the step (6) comprises:
(6-1) reading any unread tomographic section number m in the left tomographic section number set LF to correspond to the tomographic section f m Taking the leftmost end point of the fault section f as a starting point m The upper points are rearranged in sequence to form a point set FM= { f m,p |p＝1,2,…,fpn}；
(6-2) reading any unread tomographic section number n in the right tomographic section number set RF to correspond to the tomographic section f n Starting from the rightmost end point of (a), segment f n The upper points are rearranged to form a point set FN= { f n,q |q＝1,2,…,fqn}；
(6-3) if fault section f m Belonging fault number and fault section f n The fault numbers are the same, and the step (6-4) is executed; otherwise, executing the step (6-2);
(6-4) reading the last endpoint f of the Point set FM m,fpn The last endpoint f of the set of read points FN n,fqn At f m,fpn And f n,fqn Is used as a reference point, and each point of the point set FM and the point set FN is adjusted according to the following steps;
(6-5) generating the adjusted point set FM and the point set FN into the corresponding fault section f m 、f n Merging f m 、f n To one, add to the set of fault segments F, and remove F from the set F m 、f n ；
(6-6) removing the fault section number m from the set LF and the fault section number n from the set RF;
(6-7) circularly executing the steps (6-1) to (6-6) until the LF is traversed;
(6-8) reading any unread fault section number m in the set LF, and adjusting the fault section f according to the following m Each point is located above, and the adjusted fault section is taken as a recovered fault section f m ；
Wherein,representing fault segmentation f i Is a displacement vector of (a);
(6-9) reading any unread fault section number n in the set RF, and adjusting the fault section f according to the following n Each point is located above, and the adjusted fault section is taken as a recovered fault section f n ；
(6-10) performing steps (6-8) to (6-9) in a loop until the sets LF and RF are traversed;
(6-11) removing the fault section F from the fault section set F m 。
7. An automated fault construction restoration device based on a geological map, comprising a processor and a computer program stored on a memory and executable on the processor, characterized in that: the processor, when executing the program, implements the method of any one of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011000473.3A CN112231423B (en) | 2020-09-22 | 2020-09-22 | Geological map-based fault structure automatic recovery method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011000473.3A CN112231423B (en) | 2020-09-22 | 2020-09-22 | Geological map-based fault structure automatic recovery method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112231423A CN112231423A (en) | 2021-01-15 |
| CN112231423B true CN112231423B (en) | 2024-02-23 |
Family
ID=74107411
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011000473.3A Active CN112231423B (en) | 2020-09-22 | 2020-09-22 | Geological map-based fault structure automatic recovery method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112231423B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113282975A (en) * | 2021-02-25 | 2021-08-20 | 山西焦煤集团有限责任公司 | Method for automatically processing geological profile fault and stratum space morphology |
| CN113539051B (en) * | 2021-06-24 | 2022-11-25 | 南京师范大学 | Geological map-based method and device for acquiring stratum boundary point-by-point rock stratum attitude |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
| CN106384016A (en) * | 2016-10-08 | 2017-02-08 | 南京师范大学 | Automatic inversion method of single dip fault |
| CN106990450A (en) * | 2017-05-13 | 2017-07-28 | 东北石油大学 | A kind of method for recovering fault development history of evolution |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN111443385A (en) * | 2020-04-28 | 2020-07-24 | 南京师范大学 | Method for acquiring fault development time sequence |
-
2020
- 2020-09-22 CN CN202011000473.3A patent/CN112231423B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN106384016A (en) * | 2016-10-08 | 2017-02-08 | 南京师范大学 | Automatic inversion method of single dip fault |
| CN106990450A (en) * | 2017-05-13 | 2017-07-28 | 东北石油大学 | A kind of method for recovering fault development history of evolution |
| CN111443385A (en) * | 2020-04-28 | 2020-07-24 | 南京师范大学 | Method for acquiring fault development time sequence |
Non-Patent Citations (2)
| Title |
|---|
| 几何与语义统一的区域地质构造GIS数据模型;何撼东等;测绘学报;全文 * |
| 地质断层三维构模技术研究;朱良峰;潘信;;岩土力学(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112231423A (en) | 2021-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112231423B (en) | Geological map-based fault structure automatic recovery method and device | |
| Owen | A survey of unstructured mesh generation technology. | |
| CN104635262B (en) | A kind of positive reversed fault isopleth automatic generation method based on extended rectangular grid | |
| CN104851126B (en) | Threedimensional model dividing method and device based on generalized cylinder | |
| Gold et al. | Outside-in: an alternative approach to forest map digitizing | |
| CN104240300A (en) | Large-scale point cloud complex space curved face reconstruction method based on distributed parallelism | |
| CN112464479A (en) | Three-dimensional ore formation prediction method and system based on GOCAD | |
| CN110349262B (en) | BIM-based three-dimensional modeling method and system for advanced geological prediction of tunnel | |
| CN110197013B (en) | Riverbed bedrock surface modeling method based on Morphing | |
| CN112233230B (en) | Three-dimensional model construction method and device for fault structure in cut geological section | |
| CN112116709B (en) | Terrain characteristic line processing method for improving terrain expression precision | |
| CN111951393B (en) | Three-dimensional model construction method and device for horizontally constructing bench | |
| de Oliveira Miranda et al. | Finite element mesh generation for subsurface simulation models | |
| CN104766367B (en) | Three-dimensional gridding topology structure chart construction method in computing three-dimensional model processing | |
| CN115841022A (en) | Data processing method based on face rockfill dam and related equipment | |
| CN106384016B (en) | The automatic inversion method of single dip fault | |
| JP2006009262A (en) | System for creating and displaying geologic structure model, stratum modeling method, and program for creating and displaying geologic structure model | |
| CN112927350B (en) | Multi-domain geologic body model construction method, device, equipment and storage medium | |
| DU et al. | Terrain edge stitching based on least squares generative adversarial networks | |
| Zeng et al. | Developable Surface Segmentation For CAD Models | |
| Altantsetseg et al. | Minimum surface area based complex hole filling algorithm of 3D mesh | |
| CN114565732B (en) | Three-dimensional modeling method and device for generation layer of branch-shaped distribution soil | |
| CN110473270A (en) | A kind of method of Intelligent Reconstruction 2D Stratigraphical sequence model | |
| Song et al. | Integrated Quality Mesh Generation for Poisson Surface Reconstruction in HPC Applications | |
| CN117173357A (en) | Mine three-dimensional geological modeling and layered cutting method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |