CN112231423A - Fault structure automatic recovery method and device based on geological map - Google Patents

Fault structure automatic recovery method and device based on geological map Download PDF

Info

Publication number
CN112231423A
CN112231423A CN202011000473.3A CN202011000473A CN112231423A CN 112231423 A CN112231423 A CN 112231423A CN 202011000473 A CN202011000473 A CN 202011000473A CN 112231423 A CN112231423 A CN 112231423A
Authority
CN
China
Prior art keywords
fault
rpl
lpl
boundary
stratum
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.)
Granted
Application number
CN202011000473.3A
Other languages
Chinese (zh)
Other versions
CN112231423B (en
Inventor
李安波
徐诗宇
董甜甜
闾国年
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN202011000473.3A priority Critical patent/CN112231423B/en
Publication of CN112231423A publication Critical patent/CN112231423A/en
Application granted granted Critical
Publication of CN112231423B publication Critical patent/CN112231423B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 a device based on a geological map, wherein the method comprises the following steps: firstly, reading formation data and fault data; secondly, acquiring stratums on two sides of the fault and fault segmentation based on buffer zone neighborhood analysis; then, restoring stratums on two sides of the fault based on the corresponding relation of the stratums; and finally, restoring fault sections on two sides of the fault based on the corresponding relation of the fault sections, 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 method has the advantages of high efficiency, good fault structure recovery effect and high automation degree.

Description

Fault structure automatic recovery method and device based on geological map
Technical Field
The invention relates to the field of geographic information technology application, in particular to a fault structure automatic recovery method and device based on a geological map.
Background
The recovery of geological structure refers to the process of deducing the original shape of the rock stratum based on the structural landform type and other conditions as constraints. Faults, one of the most important structures of the earth crust, are formed by the obvious displacement of the stratum or rock along the fracture surface, and often form a valley and a scarp. The method has the advantages that the original morphology of the rock stratum without fracture is restored by aiming at scientific inference of the fault evolution process, and the method has certain research significance and application value for geological disaster prediction, geological resource development, geological environment protection and the like.
The fault recovery application mainly comprises the following three aspects: 1) geological map based fault recovery. Judging the consistency of stratum contained in a disc body in the fault, merging the displaced rock stratums, and finally forming an original geological map after the fault line is removed; 2) fault recovery based on geological 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 order of fracture occurrence of the fault; 3) and fault recovery based on the three-dimensional model. And modifying internal attribute data of the fault entity mainly, and recovering the deformed fault structure into a complete geological entity. At present, related applications are all related, but related work is mainly manually completed by geologists, investment is large, efficiency is poor, and the research of an automatic fault recovery method is urgently needed.
The geological map data has the characteristics of abundant data, wide source, wide coverage range, less use limitation and the like. Starting from the geological map, the fault structure recovery research is carried out, and the feasibility and the practicability are better.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems 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 automatic fault structure recovery method based on the geological map comprises the following steps:
(1) reading stratum data and fault line data, and generating a stratum set P, a fault segmentation set F and a fault segmentation time sequence set T;
(2) reading a fault section F from a fault section set F according to a fault section time sequencei；
(3) Acquiring fault segments fiLeft stratum number set LP and left fault subsection number setLF, right stratum number set RP and right fault segmentation number set RF;
(4) based on the corresponding relationship of the stratums, restoring the layer boundary of each stratum in the left stratum number set LP and the right stratum number set RP, and respectively storing the layer boundary into a left layer boundary set LPL and a right layer boundary set RPL;
(5) restoring fault segments f based on left-side level boundary set LPL and right-side level boundary set RPLiFault cutting boundary lines of the stratums on two sides, and recovering the stratums based on the fault cutting boundary lines and the bedding boundary lines;
(6) restoring fault segments f based on fault correspondencesiSegmenting faults on two sides;
(7) deriving a geological map based on the recovered stratigraphic and fault;
(8) and (5) circularly executing the steps (2) to (7) until the fault segmentation set F is traversed, and all the fault structures are recovered.
Further, the step (1) comprises the following steps:
(1-1) reading fault line vector data into a fault segmentation set F ═ F i1,2, …, fn }; wherein f isiRepresenting the ith fault segmentation, and fn represents the number of fault segmentation;
(1-2) acquiring the serial number of the fault segmentation, and storing the serial number in a fault segmentation serial number set FID ═ FID i1,2, …, fn }; wherein, fidiSegmenting f for faultsiThe 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 }i1,2, …, fn }; wherein, tiSegmenting f for faultsiThe timing of (c);
(1-4) reading formation data into a set of formations P ═ { P ═ P j1,2, …, pn }; wherein p isjDenotes the jth formation and pn denotes the number of formations.
Further, the step (3) comprises the following steps:
(3-1) calculating fault segmentation f according to the preset width W based on a buffer calculation interface of GIS softwareiBuffer b ofi；
(3-2) creating a tomographic section fiThe left stratum number set LP and the right stratum number set RP;
(3-3) reading any stratum P from the stratum set PjWhen p isjAnd a buffer biWhen the left boundary of (a) is intersected, the stratum number j is stored into a set LP, and when p is intersectedjAnd a buffer biWhen the right boundary of the set RP is intersected, the stratum number j is stored 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) creating a tomographic section fiThe left fault subsection number set LF and the right fault subsection number set RF;
(3-6) reading any fault section F from the fault section set FkWhen f iskAnd a buffer biWhen the left boundary of (f) is intersected, the fault segment number k is stored into a set LF, and when f is intersectedkAnd a buffer biWhen the right boundary of the fault is intersected, storing the fault segmentation number k into a set RF;
and (3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed, and obtaining a left fault segmentation number set LF and a right fault segmentation number set RF.
Further, the step (4) comprises the following steps:
(4-1) calculating the tomographic section f according to the following formulaiIs a displacement vector
Figure BDA0002694119260000022
Figure BDA0002694119260000021
In the formula, diriRepresenting fault segments fiRun of (d)iRepresenting fault segments fiThe amount of displacement of (a);
(4-2) extracting the level boundary of each stratum in the left stratum number set LP, and storing the level boundary in the left stratum boundary set LPL ═ LPLu(us, ls, flag) | u ═ 1,2, …, lpn +1 }; wherein, lpluIs shown asu level boundaries, us being level boundary lpluUpper strata type when lpluWithout upper strata, us takes the line-1, ls as the boundary lpluThe type of the lower formation of (a); when lpluWhen there is no lower stratum, ls is-1, and flag is lpluThe initial value of the recovery flag of (1) is 0, and lpn is the number of elements in the LP;
(4-3) extracting the layer boundary of each stratum in the right stratum number set RP and storing the layer boundary in the right side layer boundary set RPL ═ { RPL ═ RPLv(us, ls, flag) | v ═ 1,2, …, rpn +1 }; wherein, rplvDenotes the v-th slice boundary, us being the slice boundary rplvUpper formation type of (2), when rplvWithout upper strata, us takes the value-1, ls as the boundary line rplvThe type of the lower formation of (a); when rpl is reachedvWhen there is no lower stratum, ls is taken as-1, and flag is rplvThe initial value of the recovery flag of (3) is 0, and rpn is the number of elements in RP;
(4-4) matching the fault segments f according to the sets LPL and RPLiThe bedding boundary of the left stratum and the bedding boundary of the right stratum, and restoring the bedding boundary.
Further, the step (4-4) comprises:
(4-4-1) reading LPL the first level boundary from the set LPLuReading the first layer boundary RPL from the RPL setv；
(4-4-2) when lplu(us) and rplv(us) same or lplu(ls) and rplv(ls) when the same, executing the step (4-4-3); otherwise, executing the step (4-4-6);
(4-4-3) at lpluIs the starting point, and the level boundary line lpl is boundeduThe points are rearranged in sequence to obtain an arranged point set LPPu＝{ lpl u,g1,2, …, lgn }; in rplvThe rightmost end point of (a) is a starting point, and the layer boundary lines rplvThe points are rearranged in sequence to obtain an arranged point set RPPv＝{ rpl v,h1,2, …, lhn }, wherein lplu,gRepresentation lpluThe upper g point, lgn for lpluNumber of points, rplv,hDenotes rplvThe upper h point, lhn denotes rplvCounting;
(4-4-4) reading LPPuLplu,lgnAnd RPPvLast point in (rpl)v,lhnLpl will beu,lgnAnd rplv,lhnIs taken as a reference point and the layer boundary line lpl is adjusted according to the following formulau、rplvEach point is arranged;
Figure BDA0002694119260000031
(4-4-5) lplu(flag) and rplv(flag) is set to 1, and when u is equal to lpn +1 or v is equal to rpn +1, step (4-4-7) is performed; otherwise, increasing the value of u by 1, increasing the value of v by 1, and backtracking to the step (4-4-2);
(4-4-6) sequentially traversing the level boundary lines rplvTo rplrpn+1If rpl is presenttSatisfy lplu(us) and rplt(us) same or lplu(ls) and rplt(ls) setting v as t, and backtracking to the step (4-4-2) in the same way; otherwise, increasing the value of u by 1, and backtracking to the step (4-4-2);
(4-4-7) setting u to 1 based on the displacement vector
Figure BDA0002694119260000032
Lpl was adjusted according touAnd the adjusted level boundary is used as the recovered level boundary lpluAnd will lplu(flag) is set to 2;
Figure BDA0002694119260000041
(4-4-8) setting v to 1 based on the displacement vector
Figure BDA0002694119260000042
Adjusting rpl according tovAnd the adjusted level boundary line is used as the recovered level boundary line rplvAnd combining rplv(flag) is set to 2;
Figure BDA0002694119260000043
(4-4-9) incrementing u by 1, incrementing v by 1, and looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is complete.
Further, the step (5) comprises the following steps:
(5-1) extracting a 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 into a stratum side boundary line set SL;
(5-2) restoring the fault cutting boundary line based on the left side layer boundary line set LPL and the right side layer boundary line set RPL, and storing the fault cutting boundary line into the cutting boundary line set CL;
(5-3) constructing fault segments f based on the left-side plane boundary set LPL, the right-side plane boundary set RPL, the cutting boundary set CL and the stratum-side boundary set SLiThe surface of the ground is stored into a surface set TP;
(5-4) for each stratum TP in the set TP, acquiring an upper stratum boundary of TP based on the left side stratum boundary set LPL and the right side stratum boundary set RPL, and writing a lower stratum type of the upper stratum boundary into the TP stratum type attribute; and adding tp to the stratigraphic collection P;
(5-5) removing the fault section f from the stratum set P according to the left stratum number set LP and the right stratum number set RPiThe two side strata.
Further, the step (5-2) comprises:
(5-2-1) continuously reading two slice boundary lines LPL in the left slice boundary set LPLu、lplu+1When lplu(flag) 2 or lplu+1When (flag) is 2, connection lpluRightmost endpoint lplu+1Forming a connecting line at the extreme point at the rightmost side, wherein the connecting line is a fault cutting boundary line and is added to the set CL;
(5-2-2) continuously reading two slice boundary lines RPL in the right slice boundary set RPLv、rplv+1When rpl is usedv(flag) 2 or rplv+1(flag)＝At 2, connecting rplvLeftmost endpoint and rplv+1Forming a connecting line at the extreme point at the leftmost side, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
and (5-2-3) circularly executing the steps (5-2-1) to (5-2-2) until the LPL and the RPL are traversed.
Further, the step (6) comprises the following steps:
(6-1) reading any unread fault segment number m in the left fault segment number set LF to obtain the corresponding fault segment fmIs taken as a starting point, and the fault is segmented into fmThe points are rearranged in sequence to form a point set FM ═ fm,p|p＝1,2,…,fpn}；
(6-2) reading any unread fault segment number n in the right fault segment number set RF to obtain the corresponding fault segment fnIs taken as a starting point, and the fault is segmented into fnRearranging each point to form point set FN ═ fn,q|q＝1,2,…,fqn}；
(6-3) if fault section fmThe fault number and fault section fnThe fault numbers of the faults are the same, and the step (6-4) is executed; otherwise, executing the step (6-2);
(6-4) reading the last end point f of the set of points FMm,fpnLast end point f of set of read points FNn,fqnWith fm,fpnAnd fn,fqnThe central point fp of (2) is used as a reference point, and the points of the point set FM and the point set FN are adjusted according to the following formula;
Figure BDA0002694119260000051
(6-5) generating corresponding fault segmentation f by using the adjusted point set FM and point set FNm、fnMerging fm、fnOne, added to the fault segmentation set F, and F is removed from the set Fm、fn；
(6-6) removing the number m of the fault segments from the set LF and removing the number n of the fault segments 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 segment number m in the set LF, and adjusting the fault segment f according to the following formulamTaking the adjusted fault section as a recovered fault section fm；
Figure BDA0002694119260000052
Wherein the content of the first and second substances,
Figure BDA0002694119260000053
representing fault segments fiThe displacement vector of (2);
(6-9) reading any unread tomographic segment number n in the set RF, and adjusting the tomographic segment f according to the following formulanTaking the adjusted fault section as a recovered fault section fn；
Figure BDA0002694119260000054
(6-10) executing steps (6-8) to (6-9) circularly until the sets LF and RF are traversed;
(6-11) removing the fault segments F from the set of fault segments Fm。
The automatic fault structure recovery device based on the geological map comprises a processor and a computer program which is stored on a memory and can run on the processor, wherein the processor realizes the method when executing the program.
Has the advantages that: compared with the prior art, the invention has the following remarkable advantages: the automation degree is high, the efficiency is high, and the effect is good.
Drawings
FIG. 1 is geological map data used in the present embodiment;
FIG. 2 is fault attribute data employed in the present embodiment;
FIG. 3 is a schematic flow diagram of the present invention;
FIG. 4 is a schematic illustration of various boundaries of a formation;
FIG. 5 is a series of geological maps generated during fault recovery.
Detailed Description
As will be described in further detail below, in this embodiment, geological map data (fig. 1 and 2) from xiaomaoshan to xiangshan in Nanjing is selected as experimental data, and a projection coordinate system adopted by the experimental data is Nanjing 92 coordinate system. The following further description is provided by describing a specific embodiment in conjunction with the accompanying drawings.
As shown in fig. 3, the present embodiment provides a method for automatically recovering a fault structure based on a geological map, which includes the following steps:
(1) and reading the 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 fault line vector data into a fault segmentation set F ═ F i1,2, …, fn }; wherein f isiRepresenting the ith fault segmentation, and fn represents the number of fault segmentation; in this embodiment, fn is 6;
(1-2) acquiring the serial number of the fault segmentation, and storing the serial number in a fault segmentation serial number set FID ═ FID i1,2, …, fn }; wherein, fidiSegmenting f for faultsiThe 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 }i1,2, …, fn }; wherein, tiSegmenting f for faultsiThe timing of (c);
(1-4) reading formation data into a set of formations P ═ { P ═ P j1,2, …, pn }; wherein p isjRepresents the jth stratum, pn represents the number of strata, and pn is 45 in this embodiment.
(2) Reading a fault section F from a fault section set F according to a fault section time sequencei。
(3) Acquiring fault segments fiLeft stratum number set LP, left fault subsection number set LF, right stratum number set RP and right fault subsectionThe set of numbers RF.
The method specifically comprises the following steps:
(3-1) calculating fault segmentation f according to the preset width W based on a buffer calculation interface of GIS softwareiBuffer b ofi(ii) a In the embodiment, W is 20 meters;
(3-2) creating a tomographic section fiThe left stratum number set LP and the right stratum number set RP;
(3-3) reading any stratum P from the stratum set PjWhen p isjAnd a buffer biWhen the left boundary of (a) is intersected, the stratum number j is stored into a set LP, and when p is intersectedjAnd a buffer biWhen the right boundary of the set RP is intersected, the stratum number j is stored 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) creating a tomographic section fiThe left fault subsection number set LF and the right fault subsection number set RF;
(3-6) reading any fault section F from the fault section set FkWhen f iskAnd a buffer biWhen the left boundary of (f) is intersected, the fault segment number k is stored into a set LF, and when f is intersectedkAnd a buffer biWhen the right boundary of the fault is intersected, storing the fault segmentation number k into a set RF;
and (3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed, and obtaining a left fault segmentation number set LF and a right fault segmentation number set RF.
(4) And based on the corresponding relation of the stratums, restoring the layer boundary of each stratum in the left stratum number set LP and the right stratum number set RP, and respectively storing the layer boundary into the left layer boundary set LPL and the right layer boundary set RPL.
The method comprises the following steps:
(4-1) calculating the tomographic section f according to the following formulaiIs a displacement vector
Figure BDA0002694119260000071
Figure BDA0002694119260000072
In the formula, diriRepresenting fault segments fiRun of (d)iRepresenting fault segments fiThe amount of displacement of (a);
(4-2) extracting the level boundary of each stratum in the left stratum number set LP, and storing the level boundary in the left stratum boundary set LPL ═ LPLu(us, ls, flag) | u ═ 1,2, …, lpn +1 }; wherein, lpluDenotes the u-th level boundary, us is level boundary lpluUpper strata type when lpluWithout upper strata, us takes the line-1, ls as the boundary lpluThe type of the lower formation of (a); when lpluWhen there is no lower stratum, ls is-1, and flag is lpluThe initial value of the recovery flag of (1) is 0, and lpn is the number of elements in the LP; in this example, the bedding boundaries of the strata are as shown in FIG. 4;
(4-3) extracting the layer boundary of each stratum in the right stratum number set RP and storing the layer boundary in the right side layer boundary set RPL ═ { RPL ═ RPLv(us, ls, flag) | v ═ 1,2, …, rpn +1 }; wherein, rplvDenotes the v-th slice boundary, us being the slice boundary rplvUpper formation type of (2), when rplvWithout upper strata, us takes the value-1, ls as the boundary line rplvThe type of the lower formation of (a); when rpl is reachedvWhen there is no lower stratum, ls is taken as-1, and flag is rplvThe initial value of the recovery flag of (3) is 0, and rpn is the number of elements in RP;
(4-4) matching the fault segments f according to the sets LPL and RPLiThe bedding boundary of the left stratum and the bedding boundary of the right stratum, and restoring the bedding boundary. The method specifically comprises the following steps:
(4-4-1) reading LPL the first level boundary from the set LPLuReading the first layer boundary RPL from the RPL setv；
(4-4-2) when lplu(us) and rplv(us) same or lplu(ls) and rplv(ls) when the same, executing the step (4-4-3); otherwise, step (4-4-6) is performed, wherein (#) denotes an elementAttribute #;
(4-4-3) at lpluIs the starting point, and the level boundary line lpl is boundeduThe points are rearranged in sequence to obtain an arranged point set LPPu＝{ lpl u,g1,2, …, lgn }; in rplvThe rightmost end point of (a) is a starting point, and the layer boundary lines rplvThe points are rearranged in sequence to obtain an arranged point set RPPv＝{ rpl v,h1,2, …, lhn }, wherein lplu,gRepresentation lpluThe upper g point, lgn for lpluNumber of points, rplv,hDenotes rplvThe upper h point, lhn denotes rplvCounting;
(4-4-4) reading LPPuLplu,lgnAnd RPPvLast point in (rpl)v,lhnLpl will beu,lgnAnd rplv,lhnIs taken as a reference point and the layer boundary line lpl is adjusted according to the following formulau、rplvEach point is arranged;
Figure BDA0002694119260000081
(4-4-5) lplu(flag) and rplv(flag) is set to 1, and when u is equal to lpn +1 or v is equal to rpn +1, step (4-4-7) is performed; otherwise, increasing the value of u by 1, increasing the value of v by 1, and backtracking to the step (4-4-2);
(4-4-6) sequentially traversing the level boundary lines rplvTo rplrpn+1If rpl is presenttSatisfy lplu(us) and rplt(us) same or lplu(ls) and rplt(ls) setting v as t, and backtracking to the step (4-4-2) in the same way; otherwise, increasing the value of u by 1, and backtracking to the step (4-4-2);
(4-4-7) setting u to 1 based on the displacement vector
Figure BDA0002694119260000082
Lpl was adjusted according touAnd the adjusted level boundary is used as the recovered level boundary lpluAnd will lplu(flag) is set to 2;
Figure BDA0002694119260000083
(4-4-8) setting v to 1 based on the displacement vector
Figure BDA0002694119260000084
Adjusting rpl according tovAnd the adjusted level boundary line is used as the recovered level boundary line rplvAnd combining rplv(flag) is set to 2;
Figure BDA0002694119260000085
(4-4-9) incrementing u by 1, incrementing v by 1, and looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is complete.
(5) Restoring fault segments f based on left-side level boundary set LPL and right-side level boundary set RPLiAnd fault cutting boundary lines of the stratums on two sides, and restoring the stratums based on the fault cutting boundary lines and the bedding boundary lines.
The method specifically comprises the following steps:
(5-1) extracting a 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 into a stratum side boundary line set SL;
(5-2) restoring the fault cutting boundary line based on the left side layer boundary line set LPL and the right side layer boundary line set RPL, and storing the fault cutting boundary line into the cutting boundary line set CL;
(5-3) constructing fault segments f by adopting GIS software based on left side layer boundary line set LPL, right side layer boundary line set RPL, cutting boundary line set CL and stratum side boundary line set SLiThe surface of the ground is stored into a surface set TP;
(5-4) for each stratum TP in the set TP, acquiring an upper stratum boundary of TP based on the left side stratum boundary set LPL and the right side stratum boundary set RPL, and writing a lower stratum type of the upper stratum boundary into the TP stratum type attribute; and adding tp to the stratigraphic collection P;
(5-5) removing the fault section f from the stratum set P according to the left stratum number set LP and the right stratum number set RPiThe two side strata.
The step (5-2) comprises the following steps:
(5-2-1) continuously reading two slice boundary lines LPL in the left slice boundary set LPLu、lplu+1When lplu(flag) 2 or lplu+1When (flag) is 2, connection lpluRightmost endpoint lplu+1Forming a connecting line at the extreme point at the rightmost side, wherein the connecting line is a fault cutting boundary line and is added to the set CL;
(5-2-2) continuously reading two slice boundary lines RPL in the right slice boundary set RPLv、rplv+1When rpl is usedv(flag) 2 or rplv+1When (flag) is 2, rpl is connectedvLeftmost endpoint and rplv+1Forming a connecting line at the extreme point at the leftmost side, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
and (5-2-3) circularly executing the steps (5-2-1) to (5-2-2) until the LPL and the RPL are traversed.
(6) Restoring fault segments f based on fault correspondencesiThe fault on both sides is segmented.
The method specifically comprises the following steps:
(6-1) reading any unread fault segment number m in the left fault segment number set LF to obtain the corresponding fault segment fmIs taken as a starting point, and the fault is segmented into fmThe points are rearranged in sequence to form a point set FM ═ fm,p|p＝1,2,…,fpn}；
(6-2) reading any unread fault segment number n in the right fault segment number set RF to obtain the corresponding fault segment fnIs taken as a starting point, and the fault is segmented into fnRearranging each point to form point set FN ═ fn,q|q＝1,2,…,fqn}；
(6-3) if fault section fmThe fault number and fault section fnThe fault numbers of the faults are the same, and the step (6-4) is executed; otherwise, executeStep (6-2);
(6-4) reading the last end point f of the set of points FMm,fpnLast end point f of set of read points FNn,fqnWith fm,fpnAnd fn,fqnThe central point fp of (2) is used as a reference point, and the points of the point set FM and the point set FN are adjusted according to the following formula;
Figure BDA0002694119260000091
(6-5) generating corresponding fault segmentation f by using the adjusted point set FM and point set FNm、fnMerging fm、fnOne, added to the fault segmentation set F, and F is removed from the set Fm、fn；
(6-6) removing the number m of the fault segments from the set LF and removing the number n of the fault segments 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 segment number m in the set LF, and adjusting the fault segment f according to the following formulamTaking the adjusted fault section as a recovered fault section fm；
Figure BDA0002694119260000101
Wherein the content of the first and second substances,
Figure BDA0002694119260000102
representing fault segments fiThe displacement vector of (2);
(6-9) reading any unread tomographic segment number n in the set RF, and adjusting the tomographic segment f according to the following formulanTaking the adjusted fault section as a recovered fault section fn；
Figure BDA0002694119260000103
(6-10) executing steps (6-8) to (6-9) circularly until the sets LF and RF are traversed;
(6-11) removing the fault segments F from the set of fault segments Fm。
(7) A geological map is derived based on the recovered formations and faults.
(8) And (5) circularly executing the steps (2) to (7) until the fault segmentation set F is traversed, and all the fault structures are recovered. As shown in fig. 5, in the present embodiment, 5 fault recovery process geological maps are formed jointly.
The embodiment also provides a geological map-based fault structure automatic recovery device which comprises a processor and a computer program stored on a memory and capable of running on the processor, wherein the processor executes the program to realize the method.
In the embodiment of the invention, the buffer area is only created based on the interface provided by the GIS software ArcGIS, and the method can also use the API of the software such as QGIS, MapGIS and the like.
While the invention has been described in connection with what is presently considered to be the most practical and preferred embodiment, it is to be understood that the invention is not to be limited to the disclosed embodiment, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.

Claims (9)

1. A fault structure automatic recovery method based on geological maps is characterized by comprising the following steps:
(1) reading stratum data and fault line data, and generating a stratum set P, a fault segmentation set F and a fault segmentation time sequence set T;
(2) reading a fault section F from a fault section set F according to a fault section time sequencei；
(3) Acquiring fault segments fiThe left stratum number set LP, the left fault subsection number set LF, the right stratum number set RP and the right fault subsection number set RF;
(4) based on the corresponding relationship of the stratums, restoring the layer boundary of each stratum in the left stratum number set LP and the right stratum number set RP, and respectively storing the layer boundary into a left layer boundary set LPL and a right layer boundary set RPL;
(5) restoring fault segments f based on left-side level boundary set LPL and right-side level boundary set RPLiFault cutting boundary lines of the stratums on two sides, and recovering the stratums based on the fault cutting boundary lines and the bedding boundary lines;
(6) restoring fault segments f based on fault correspondencesiSegmenting faults on two sides;
(7) deriving a geological map based on the recovered stratigraphic and fault;
(8) and (5) circularly executing the steps (2) to (7) until the fault segmentation set F is traversed, and all the fault structures are recovered.
2. The geological map-based fault structure automatic recovery method according to claim 1, characterized in that: the step (1) comprises the following steps:
(1-1) reading fault line vector data into a fault segmentation set F ═ Fi1,2, …, fn }; wherein f isiRepresenting the ith fault segmentation, and fn represents the number of fault segmentation;
(1-2) acquiring the serial number of the fault segmentation, and storing the serial number in a fault segmentation serial number set FID ═ FIDi1,2, …, fn }; wherein, fidiSegmenting f for faultsiThe 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 }i1,2, …, fn }; wherein, tiSegmenting f for faultsiThe timing of (c);
(1-4) reading formation data into a set of formations P ═ { P ═ Pj1,2, …, pn }; wherein p isjDenotes the jth formation and pn denotes the number of formations.
3. The geological map-based fault structure automatic recovery method according to claim 1, characterized in that: the step (3) comprises the following steps:
(3-1) calculating fault segmentation f according to the preset width W based on a buffer calculation interface of GIS softwareiBuffer b ofi；
(3-2) creating a tomographic section fiLeft stratum number set LP and rightA stratum number set RP;
(3-3) reading any stratum P from the stratum set PjWhen p isjAnd a buffer biWhen the left boundary of (a) is intersected, the stratum number j is stored into a set LP, and when p is intersectedjAnd a buffer biWhen the right boundary of the set RP is intersected, the stratum number j is stored 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) creating a tomographic section fiThe left fault subsection number set LF and the right fault subsection number set RF;
(3-6) reading any fault section F from the fault section set FkWhen f iskAnd a buffer biWhen the left boundary of (f) is intersected, the fault segment number k is stored into a set LF, and when f is intersectedkAnd a buffer biWhen the right boundary of the fault is intersected, storing the fault segmentation number k into a set RF;
and (3-7) circularly executing the step (3-6) until the fault segmentation set F is traversed, and obtaining a left fault segmentation number set LF and a right fault segmentation number set RF.
4. The geological map-based fault structure automatic recovery method according to claim 1, characterized in that: the step (4) comprises the following steps:
(4-1) calculating the tomographic section f according to the following formulaiIs a displacement vector
Figure FDA0002694119250000021
Figure FDA0002694119250000022
In the formula, diriRepresenting fault segments fiRun of (d)iRepresenting fault segments fiThe amount of displacement of (a);
(4-2) extracting the layer boundary of each stratum in the left stratum number set LP and storing the layer boundary in the left sideLevel boundary set LPL (LPL) } ═ LPLu(us, ls, flag) | u ═ 1,2, …, lpn +1 }; wherein, lpluDenotes the u-th level boundary, us is level boundary lpluUpper strata type when lpluWithout upper strata, us takes the line-1, ls as the boundary lpluThe type of the lower formation of (a); when lpluWhen there is no lower stratum, ls is-1, and flag is lpluThe initial value of the recovery flag of (1) is 0, and lpn is the number of elements in the LP;
(4-3) extracting the layer boundary of each stratum in the right stratum number set RP and storing the layer boundary in the right side layer boundary set RPL ═ { RPL ═ RPLv(us, ls, flag) | v ═ 1,2, …, rpn +1 }; wherein, rplvDenotes the v-th slice boundary, us being the slice boundary rplvUpper formation type of (2), when rplvWithout upper strata, us takes the value-1, ls as the boundary line rplvThe type of the lower formation of (a); when rpl is reachedvWhen there is no lower stratum, ls is taken as-1, and flag is rplvThe initial value of the recovery flag of (3) is 0, and rpn is the number of elements in RP;
(4-4) matching the fault segments f according to the sets LPL and RPLiThe bedding boundary of the left stratum and the bedding boundary of the right stratum, and restoring the bedding boundary.
5. The geological map-based fault structure automatic recovery method according to claim 4, characterized in that: the step (4-4) comprises the following steps:
(4-4-1) reading LPL the first level boundary from the set LPLuReading the first layer boundary RPL from the RPL setv；
(4-4-2) when lplu(us) and rplv(us) same or lplu(ls) and rplv(ls) when the same, executing the step (4-4-3); otherwise, executing the step (4-4-6);
(4-4-3) at lpluIs the starting point, and the level boundary line lpl is boundeduThe points are rearranged in sequence to obtain an arranged point set LPPu＝{lplu,g1,2, …, lgn }; in rplvThe rightmost end point of (a) is a starting point, and the layer boundary lines rplvThe upper points are rearranged in sequence to obtainArranged point set RPPv＝{rplv,h1,2, …, lhn }, wherein lplu,gRepresentation lpluThe upper g point, lgn for lpluNumber of points, rplv,hDenotes rplvThe upper h point, lhn denotes rplvCounting;
(4-4-4) reading LPPuLplu,lgnAnd RPPvLast point in (rpl)v,lhnLpl will beu,lgnAnd rplv,lhnIs taken as a reference point and the layer boundary line lpl is adjusted according to the following formulau、rplvEach point is arranged;
Figure FDA0002694119250000031
(4-4-5) lplu(flag) and rplv(flag) is set to 1, and when u is equal to lpn +1 or v is equal to rpn +1, step (4-4-7) is performed; otherwise, increasing the value of u by 1, increasing the value of v by 1, and backtracking to the step (4-4-2);
(4-4-6) sequentially traversing the level boundary lines rplvTo rplrpn+1If rpl is presenttSatisfy lplu(us) and rplt(us) same or lplu(ls) and rplt(ls) setting v as t, and backtracking to the step (4-4-2) in the same way; otherwise, increasing the value of u by 1, and backtracking to the step (4-4-2);
(4-4-7) setting u to 1 based on the displacement vector
Figure FDA0002694119250000032
Lpl was adjusted according touAnd the adjusted level boundary is used as the recovered level boundary lpluAnd will lplu(flag) is set to 2;
Figure FDA0002694119250000033
(4-4-8) setting v to 1 based on the displacement vector
Figure FDA0002694119250000034
Adjusting rpl according tovAnd the adjusted level boundary line is used as the recovered level boundary line rplvAnd combining rplv(flag) is set to 2;
Figure FDA0002694119250000035
(4-4-9) incrementing u by 1, incrementing v by 1, and looping through steps (4-4-7) through (4-4-8) until recovery of all level boundaries in the LPL and RPL is complete.
6. The geological map-based fault structure automatic recovery method according to claim 1, characterized in that: the step (5) comprises the following steps:
(5-1) extracting a 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 into a stratum side boundary line set SL;
(5-2) restoring the fault cutting boundary line based on the left side layer boundary line set LPL and the right side layer boundary line set RPL, and storing the fault cutting boundary line into the cutting boundary line set CL;
(5-3) constructing fault segments f based on the left-side plane boundary set LPL, the right-side plane boundary set RPL, the cutting boundary set CL and the stratum-side boundary set SLiThe surface of the ground is stored into a surface set TP;
(5-4) for each stratum TP in the set TP, acquiring an upper stratum boundary of TP based on the left side stratum boundary set LPL and the right side stratum boundary set RPL, and writing a lower stratum type of the upper stratum boundary into the TP stratum type attribute; and adding tp to the stratigraphic collection P;
(5-5) removing the fault section f from the stratum set P according to the left stratum number set LP and the right stratum number set RPiThe two side strata.
7. The geological map-based fault structure automated restoration method according to claim 6, wherein: the step (5-2) comprises the following steps:
(5-2-1) continuously reading two slice boundary lines LPL in the left slice boundary set LPLu、lplu+1When lplu(flag) 2 or lplu+1When (flag) is 2, connection lpluRightmost endpoint lplu+1Forming a connecting line at the extreme point at the rightmost side, wherein the connecting line is a fault cutting boundary line and is added to the set CL;
(5-2-2) continuously reading two slice boundary lines RPL in the right slice boundary set RPLv、rplv+1When rpl is usedv(flag) 2 or rplv+1When (flag) is 2, rpl is connectedvLeftmost endpoint and rplv+1Forming a connecting line at the extreme point at the leftmost side, wherein the connecting line is a fault cutting boundary line, and adding the connecting line to the set CL;
and (5-2-3) circularly executing the steps (5-2-1) to (5-2-2) until the LPL and the RPL are traversed.
8. The geological map-based fault structure automatic recovery method according to claim 1, characterized in that: the step (6) comprises the following steps:
(6-1) reading any unread fault segment number m in the left fault segment number set LF to obtain the corresponding fault segment fmIs taken as a starting point, and the fault is segmented into fmThe points are rearranged in sequence to form a point set FM ═ fm,p|p＝1,2,…,fpn}；
(6-2) reading any unread fault segment number n in the right fault segment number set RF to obtain the corresponding fault segment fnIs taken as a starting point, and the fault is segmented into fnRearranging each point to form point set FN ═ fn,q|q＝1,2,…,fqn}；
(6-3) if fault section fmThe fault number and fault section fnThe fault numbers of the faults are the same, and the step (6-4) is executed; otherwise, executing the step (6-2);
(6-4) reading the last end point f of the set of points FMm,fpnLast end point f of set of read points FNn,fqnWith fm,fpnAnd fn,fqnThe midpoint fp of (a) as a reference point,adjusting each point of the point set FM and the point set FN according to the following formula;
Figure FDA0002694119250000041
(6-5) generating corresponding fault segmentation f by using the adjusted point set FM and point set FNm、fnMerging fm、fnOne, added to the fault segmentation set F, and F is removed from the set Fm、fn；
(6-6) removing the number m of the fault segments from the set LF and removing the number n of the fault segments 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 segment number m in the set LF, and adjusting the fault segment f according to the following formulamTaking the adjusted fault section as a recovered fault section fm；
Figure FDA0002694119250000051
Wherein the content of the first and second substances,
Figure FDA0002694119250000052
representing fault segments fiThe displacement vector of (2);
(6-9) reading any unread tomographic segment number n in the set RF, and adjusting the tomographic segment f according to the following formulanTaking the adjusted fault section as a recovered fault section fn；
Figure FDA0002694119250000053
(6-10) executing steps (6-8) to (6-9) circularly until the sets LF and RF are traversed;
(6-11) removing the fault segments F from the set of fault segments Fm。
9. An apparatus for automated fault structure restoration based on geological maps, comprising a processor and a computer program stored on a memory and executable on the processor, wherein: the processor, when executing the program, implements the method of any of claims 1-8.
CN202011000473.3A 2020-09-22 2020-09-22 Geological map-based fault structure automatic recovery method and device Active CN112231423B (en)

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 true CN112231423A (en) 2021-01-15
CN112231423B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
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
CN113539051A (en) * 2021-06-24 2021-10-22 南京师范大学 Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device

Citations (5)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
何撼东等: "几何与语义统一的区域地质构造GIS数据模型", 测绘学报 *
朱良峰;潘信;: "地质断层三维构模技术研究", 岩土力学, no. 01 *

Cited By (2)

* Cited by examiner, † Cited by third party
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
CN113539051A (en) * 2021-06-24 2021-10-22 南京师范大学 Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device

Also Published As

Publication number Publication date
CN112231423B (en) 2024-02-23

Similar Documents

Publication Publication Date Title
CN112381937B (en) Multi-source geological data coupling modeling method based on drilling and complex geological profile
Owen A survey of unstructured mesh generation technology.
US9507893B2 (en) Tool for rapid configuration of a river model using imagery-based information
CN104635262B (en) A kind of positive reversed fault isopleth automatic generation method based on extended rectangular grid
US8674984B2 (en) Method for generating a hex-dominant mesh of a geometrically complex basin
JP4917764B2 (en) MAP INFORMATION CREATION DEVICE, MAP INFORMATION CREATION METHOD, AND PROGRAM
US10416350B2 (en) Method for exploitation of a subterranean medium in accordance with an exploitation scheme defined by an optimized representation
US9690002B2 (en) Device, system and method for geological-time refinement
CN103500468B (en) Based on the space curved surface method for reconstructing of boundary constraint
CN112231423A (en) Fault structure automatic recovery method and device based on geological map
CN102981182A (en) 2D seismic data all-horizon automatic tracking method based on unsupervised classification
CN110163965B (en) Solid three-dimensional model construction method for monoclinic rock stratum structure landform
KR101764615B1 (en) Spatial knowledge extractor and extraction method
CN104240300A (en) Large-scale point cloud complex space curved face reconstruction method based on distributed parallelism
Wang et al. Automatic identification of the critical slip surface of slopes
CN110197013B (en) Riverbed bedrock surface modeling method based on Morphing
Chentanez et al. Fast grid-free surface tracking
CN112116709B (en) Terrain characteristic line processing method for improving terrain expression precision
CN112233230A (en) Three-dimensional model construction method and device for fault structure in map-cut geological profile
CN104766367A (en) Three-dimensional gridding topology structure chart construction method in computing three-dimensional model processing
Janus et al. MKScal-system for land consolidation project based on CAD platform
JP3748268B1 (en) Method and program for creating fluid three-dimensional moving image
CN112837413A (en) Geological drilling-oriented virtual stratum inference method and device
CN112581558A (en) Model construction method and system for intrusion structure in map-cut geological profile
Tanaka et al. Estimation of fracture effective permeability by upscaling using ensemble Kalman filter and streamline simulation

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