CN105701848B - A kind of automatic generating method of stratigraphic boundary figure layer - Google Patents
A kind of automatic generating method of stratigraphic boundary figure layer Download PDFInfo
- Publication number
- CN105701848B CN105701848B CN201610023762.2A CN201610023762A CN105701848B CN 105701848 B CN105701848 B CN 105701848B CN 201610023762 A CN201610023762 A CN 201610023762A CN 105701848 B CN105701848 B CN 105701848B
- Authority
- CN
- China
- Prior art keywords
- layer
- boundary
- stratum
- surface element
- contact
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of automatic generating methods of stratigraphic boundary figure layer.This method comprises the following steps:(1) it is based on deposition, volcano stratigraphic boundary figure layer and metamorphic strata unit figure layer, hunts out all of its neighbor stratum in figure layer, builds stratum syntople figure;(2) public boundary that adjacent earth formations are extracted based on stratum syntople figure, generates stratigraphic boundary;(3) it is based on tomography figure layer, identifies fault contact relationship;(4) it is based on occurrence figure layer, identifies basic contact relation, the automation for completing stratigraphic boundary figure layer generates.It is compared with the method for manual manufacture, the stratigraphic boundary figure layer that the method for the present invention generates is not only efficient, quick, and avoids the introducing of human error as far as possible by reducing artificial degree of participation, meets the needs of geological map spatial database stratigraphic boundary figure layer making.
Description
Technical field
The invention belongs to geographical information technologies and application field, and in particular to one kind is based on deposition and volcano stratigraphic unit figure
Layer, metamorphic strata unit figure layer, occurrence figure layer and tomography figure layer, the method that automation generates stratigraphic boundary figure layer.
Background technology
During Construction of Geological Map Space Date-base, the making of stratigraphic boundary figure layer is an essential job.
In professional standard " Digital Geology figure figure layer and property file format [the DZ/T 0197- that China Geological Survey Bureau is formulated
1997] in ", the attribute data of clear ground-to-ground stratum boundary line chart layer is made that regulation.
Currently, the making of stratigraphic boundary figure layer is mainly realized by manual manufacture.Production method includes mainly two aspects
Work:First, scanning dyad, which quantifies existing geological map, obtains stratigraphic boundary, and stratum is completed into edlin to boundary line pel
The extraction in boundary line;Second is that the performance according to Different Strata contact relation on geologic map, manual identified contact relation is simultaneously belonged to
The typing of property information.The production method is easily introduced human error not there is only heavy workload, manual identified efficiency are low, such as
Contact relation identify mistake, attribute information mistake, to related stratigraphic unit figure layer boundary line it is inconsistent the shortcomings of.
During Construction of Geological Map Space Date-base, the deposition and volcano stratigraphic unit figure layer, the change quality that make emphatically
In layer unit figure layer, relatively stratum boundary line position and contact relation information are contained.Wherein, the position of stratigraphic boundary is depositing
With the public affairs for the two stratum elements for showing as that there is neighborhood relationship in volcano stratigraphic unit figure layer and metamorphic strata unit figure layer
Boundary altogether.Contact relation can be based on stratum sequencing, occurrence data, layer data etc. and carry out automatic identification.
Invention content
For this purpose, present invention is generally directed to which the heavy workload occurred when manual manufacture stratigraphic boundary figure layer, recognition efficiency are low, with
And the problems such as being easily introduced human error, it is proposed that one kind based on deposition and volcano stratigraphic unit figure layer, metamorphic strata unit figure layer,
Occurrence figure layer and tomography figure layer carry out the method that automation generates stratigraphic boundary figure layer.
The technical scheme is that:
A kind of automatic generating method of stratigraphic boundary figure layer, includes the following steps:(1) deposition, volcano ground stratum boundary are based on
Line chart layer and metamorphic strata unit figure layer hunt out all of its neighbor stratum in figure layer, build stratum syntople figure (Stratum
Adjacency Relation Graph, abbreviation SARG);(2) common edge of adjacent earth formations is extracted based on stratum syntople figure
Boundary generates stratigraphic boundary;(3) it is based on tomography figure layer, identifies fault contact relationship;(4) it is based on occurrence figure layer, identifies basic contact
Relationship, the automation for completing stratigraphic boundary figure layer generate.
The step (1) is as follows:
Step 11:Load deposition and volcano stratigraphic unit figure layer, metamorphic strata unit figure layer data, and by corresponding stratum
Face element reads in set S={ si| i=1,2 ..., n }, n is the number of stratal surface element;
Step 12:Read each stratal surface element siVertex set vi={ Pj| j=1,2 ..., m }, wherein PjIt indicates to constitute
Stratal surface element siVertex, m be vertex number;And to each stratal surface element siVertex set viCarry out following pre- place
Reason:
D) by vertex set viIn point PjBy up time needle sort;
E) judge each vertex set viIn whether there is P1=Pm, i.e., first node is equal to tail node, and if it exists, then redirects down
One step
C), opposite vertexes collection viIt is handled;Otherwise, jump procedure 13;
F) vertex set v is deletediTail node Pm, ensure that head and the tail do not connect, i.e. vi={ Pj| j=1,2 ..., m-1 };
Step 13:The syntople of stratal surface element in figure layer just sentence:If stratal surface element saVertex set va's
Minimum enclosed rectangle RaWith element sbVertex set vbMinimum enclosed rectangle RbThere is intersection point, then jump procedure 14;Otherwise, it redirects
Step 16;Wherein, a is cyclic variable, a ∈ [1, n-1], initial value 1;B is cyclic variable, b ∈ [a+1, n], initial value a+1;
Step 14:To two stratal surface element saAnd sbSyntople carefully sentenced, detailed process is as follows:
D) vertex set v is calculatedaAnd vbPublic point set T={ p1, p2..., pk, k is the number of common point, and records in T
Point is in vertex set vaMiddle serial number is denoted as N={ n1, n2..., nk};
If e) number k >=2 at the midpoints public point set T, redirect next step c);Otherwise, jump procedure 16;
F) judge to whether there is continuous period n in common point serial number collection Nt, nt+1Man Zu ∣ nt- nt+1∣=1, and if it exists,
That is stratal surface element s1And s2There are syntople, then jump procedure 15;Otherwise, jump procedure 16;
Step 15:Build the side d of stratum syntople figure:Calculate vertex set vaWith vbFocus point CaAnd Cb, point CaAnd Cb
As the vertex of stratum syntople figure, tie point CaAnd CbSide d as stratum syntople figure;
Step 16:If cyclic variable a > n-1, export SARG={ d1, d2... dw, w is stratum syntople figure
The total number on the side of SARG;Otherwise, cyclic variable b adds 1, if b≤n at this time, direct jump procedure 13;Otherwise cyclic variable a
Add 1, cyclic variable b=a+1, then jump procedure 13.
The step (2) is as follows:
Step 21:Read SARG={ d1, d2... dwIn i-th side di, i is cyclic variable, initial value 1;
Step 22:Find out stratal surface element saAnd sbVertex set va={ p1, p2..., pmAnd vb={ q1, q2...,
qn, m, n are respectively stratal surface element sa、sbNumber of vertices, and opposite vertexes collection vaAnd vbCarry out following pretreatment:
D) by vertex set va={ p1, p2..., pmAnd vb={ q1, q2..., qnIn point press up time needle sort;
E) judge vertex set vaAnd vbWhether p is met simultaneously1=pmWith q1=qn, i.e., first node is equal to tail node;If satisfied,
Then jump procedure 23;Otherwise, to ineligible vertex set vaOr vb, execute next step c);
F) opposite vertexes collection vaOr vbContinue to add first point p1Or q1, ensure end to end;
Step 23:Calculate vertex set vaAnd vbPublic point set T={ t1, t2..., tj, j is the number of common point, and remembers
Each point is recorded in public point set T respectively in vertex set vaAnd vbIn serial number, be denoted as Q respectivelyaAnd Qb;
Step 24:The segment processing of public point set T, segmentation processing step are as follows:
D) k is set as cyclic variable, k ∈ [1, j-1], initial value 1;Read serial number QaIn+1 k-th of period i', kth point
Number j';Equally, serial number Q is readbIn+1 k-th of period ii', kth period jj', if the Man Zu ∣ i'-j' Qie of ∣=1 ∣ simultaneously
∣=1 ii'-jj', then 2 points corresponding to period i' and period j' and period ii' and period jj' belong to same section it is public
Boundary;
E) cyclic variable k adds 1, if k < j, continues to execute previous step a) until public point set T segmentations terminate, and successively will
The point for belonging to same segment boundary constitutes side e, obtains public boundary Ea/b={ e1, e2.., em', m' indicates the section of public boundary
Number;If m'>2, execute next step c);Otherwise jump procedure 25;
F) compare two-by-two successively and constitute side etLast point coordinates tL(x, y) and composition side et+1First point coordinates tF
(x, y), t ∈ [1, m'], if there are tL(x, y)=tF(x, y), then by side etWith et+1Merge, side etWith et+1Belong to same section of public affairs
Boundary altogether;
Step 25:Cyclic variable i adds 1, if i≤w, cycle executes step 21-24, obtains the stratigraphic boundary collection of entire figure layer
Close E={ Ea/b, Ec/d... }, wherein Ea/b, Ec/dRespectively represent adjacent stratum saWith sb, adjacent stratum scWith sdBetween ground stratum boundary
Line.
The step (3) is as follows:
Step 31:Tomograph layer data is loaded, and with the fault line F={ f in tomography figure layer1, f2..., fmCentered on
Axis generates tomography buffering area F with distance parameter dB;
Step 32:Extraction interrupts layer line F={ f with layer data1, f2..., fmIntersection stratum syntople figure side
Set D={ d1, d2..., dn, n is the item number that stratum syntople figure intersects with fault line F;
Step 33:Search side d in line set DiCorresponding stratal surface element SaWith stratal surface element SbBetween stratigraphic boundary
Ea/b={ e1, e2..., em, wherein i is cyclic variable, i ∈ [1, n], initial value 1;
Step 34:Judge fault contact:If Ea/b={ e1, e2..., emCertain segment limit line ejIn tomography buffering area FBIt is interior,
Then boundary line ejBoth sides the contact relation of interlayer be fault contact relationship;
Step 35:Cyclic variable i adds 1, if i≤n, repeats step 33,34, until tomography all in figure layer connects
Tactile extraction finishes.
The step (4) is as follows:
Step 41:Based on formation lithology information, intrusive body and non-intruding body are distinguished;The differentiation result on each stratum is recorded
To deposition and volcano stratigraphic unit figure layer, the stratal surface element of metamorphic strata unit figure layer " rock class " attribute in;
Step 42:Based on lithostratic unit's abridged table, formation sequence Table X, and older ordinal number of stratum numerical value is smaller;
Step 43:Read i-th boundary line element E in stratigraphic boundary set Ea/bCorresponding stratal surface element saAnd sb's
Lithostratic unit RocknameaAnd Rocknameb, wherein i is cyclic variable, and i ∈ [1, k], initial value 1, k is stratigraphic boundary
Item number;
Step 44:Read stratal surface element saAnd sb" rock class " attribute;If stratal surface element saAnd sb" rock class " be
Non-intruding body, then jump procedure 45;Otherwise, jump procedure 47;
Step 45:Integrate the judgement of contact and unconformity contact;Read lithostratic unit RocknameaAnd Rocknameb
Serial number x on sequence table Xa、xbIf serial number xa、xbIt is adjacent or identical, i.e. xa=xbHuo person ∣ xa- xb∣=1, then stratal surface want
Plain saAnd sbContact relation be integrate contact;Otherwise, it is that unconformity contacts, jump procedure 46 carries out next step judgement;
Step 46:Disconformity contacts and angular unconformity contact judges, specific judgment method is as follows:
C) occurrence figure layer data are loaded, and read stratal surface element sa、sbAttitude:Move towards StrikeaWith
Strikeb, tendency DipaAnd DipbAnd inclination angle AngleaAnd Angleb;If any attitude information can not obtain, stratum
Face element saAnd sbContact relation be only judged as unconformity contact, jump procedure 48;
D) judge Strikea=StrikebAnd Dipa=DipbQie ∣ Anglea- AnglebWhether ∣≤δ is true, wherein 0 °
10 ° of < δ <;If so, then stratal surface element saAnd sbContact relation be disconformity contact;Otherwise, it is angular unconformity
Contact;
Step 47:The judgement of Deposit contact and intrusive contact;Read lithostratic unit RocknameaAnd Rocknameb
Serial number x on sequence table Xa、xb;As stratal surface element sa" rock class " attribute be intrusive body when, if xa≥xb, then connect for intrusion
It touches;Otherwise, then it is Deposit contact;Equally, as stratal surface element sb" rock class " attribute when being intrusive body, and xa< xb, then to invade
Enter contact;Otherwise, then it is Deposit contact;
Step 48:Record the pel number and contact relation attribute of stratum line feature;And circulation step 43-48, until ground
In stratum boundary line set E all boundary line elements both sides interlayer the equal end of identification of contact relation, complete stratigraphic boundary figure layer system
Make.
Beneficial effects of the present invention:It is compared with the method for manual manufacture, the stratigraphic boundary figure layer that the method for the present invention generates is not
Only efficiently, quick, and the introducing of human error is avoided by the artificial degree of participation of reduction as far as possible, meet geology map space
The needs that database stratigraphic boundary figure layer makes.In addition, the method for the present invention having but also as inspection stratigraphic boundary figure layer quality
Effect tool.
Description of the drawings
Fig. 1 stratigraphic boundary figure layer production flow diagrams of the present invention;
The Mount Lushan Fig. 2 experimental data;
The stratum syntople figure in the Mount Lushan Fig. 3;
The Mount Lushan stratigraphic boundary of Fig. 4 extractions;
Fig. 5 tomography buffering areas;
The stratum syntople figure that Fig. 6 intersects with fault line;
Fig. 7 fault contacts extract result;
The Mount Lushan Fig. 8 stratigraphic boundary thematic map.
Specific implementation mode
The example generated below in conjunction with the accompanying drawings and by describing a stratigraphic boundary figure layer automation, to further illustrate this
The effect of invention.This example selection 1:50000 Mount Lushan deposition and volcano stratigraphic unit figure layer, are broken at metamorphic strata unit figure layer
Layer figure layer and occurrence figure layer are experimental data, as shown in Figure 2.
Specific implementation process is as follows:
(1) SARG is built
Step 11:Mount Lushan deposition and volcano stratigraphic unit figure layer, metamorphic strata unit figure layer data are loaded, as shown in Figure 2.
And corresponding stratal surface element is read in into set S={ si| i=1,2 ..., 419 }, it is wanted altogether comprising 419 stratal surfaces in this example
Element;
Step 12:Read each stratal surface element siVertex set vi, and to each stratal surface element siVertex set carry out
Pretreatment, with stratal surface element s1For (Daye group):
A) by vertex set v1In press up time needle sort, obtain v1=(379043.6582,3277656.1655),
(379018.7804,3277657.9852) ..., (379043.6582,3277656.1655) }, wherein p1=
(379043.6582,3277656.1655), p18=(379043.6582,3277656.1655) include 18 vertex altogether;
B) judge vertex set v1In whether deposit P1=P18。v1Middle P1=(379043.6582,3277656.1655), P18=
(379043.6582,3277656.1655), there are P1=P18(i.e. first node is equal to tail node) then redirects next step c) to top
Point set v1It is handled;
C) vertex set v is deleted1Tail node P18, ensure that head and the tail do not connect.That is v1=(379043.6582,
3277656.1655), (379018.7804,3277657.9852) ..., (379045.2562,3277678.0115) }.
Step 13:To two stratal surface element s in figure layera(a is cyclic variable, a ∈ [Isosorbide-5-Nitrae 18], initial value 1) and sb(b is
Cyclic variable, b ∈ [2,419], initial value a+1) syntople carry out just sentence, with s1(Daye group) and s2(under New Port clay
Section) for carry out step 13-16, i.e. when a=1, b=2:Stratal surface element s1Vertex set v1Minimum enclosed rectangle R1And s2
Vertex set v2Minimum enclosed rectangle R2There are intersection point, jump procedure 14;
Step 14:To two stratal surface element saAnd sbSyntople carefully sentenced, equally, with s1With s2For, specific mistake
Journey is as follows:
A) vertex set v is calculated1With vertex set v2Public point set T=(379043.6582,3277656.1655),
(379018.7804,3277657.9852), (378995.2579,3277644.7475) ... (379043.6582,
3277656.1655), the number k=9 of common point.And the midpoints T are recorded in vertex set v1Middle serial number, be denoted as N=1,2,3,
...18};
B) number k >=2 at the midpoints public point set T, then redirect next step c);
C) there are continuous period 1,2, full ∣=1 foot ∣ 1-2, i.e. stratal surface element s in public point set serial number N1And s2It deposits
In syntople, jump procedure 15.
Step 15:Build the side of stratum syntople figure:Still with s1And s2For, calculate v1With v2Focus point C1
(378978.798955993,3277704.63529242) and C2(379006.309367991,3277476.48297586), point
C1And C2As the vertex of stratum syntople figure, C is connected1And C2A line d as stratum syntople figure1;
Step 16:Cyclic variable a=1<418, then cyclic variable b add 1.B=3≤419, jump procedure 13, judge stratum
Face element s1And s3With the presence or absence of syntople.Cycle executes step 13-15, until stratal surface element s418And s419Adjacent close
System judges to terminate, and obtains the SARG in Mount Lushan, output SARG={ d1, d2... d785, that is, there are 785 stratum syntoples, such as
Shown in Fig. 3.
(2) public boundary between adjacent stratum element, the i.e. extraction of stratigraphic boundary are extracted based on SARG
Step 21:Read SARG={ d1, d2... d785In i-th side di, i is cyclic variable, initial value 1.With d1For
Example, i=1, d1By stratal surface element s1Focus point C1With stratal surface element s2Focus point C2Connection obtains.Extraction stratal surface is wanted
Plain s1With stratal surface element s2Public boundary, execute step 22-26;
Step 22:Find out stratal surface element s1And s2Vertex set v1=(379043.6582,3277656.1655),
(379018.7804,3277657.9852) .., (379043.6582,3277656.1655) } totally 18 vertex and v2=
(379043.6582,3277656.1655), (379048.2663,3277622.0509) .., (379043.6582,
3277656.1655) } totally 23 vertex.And opposite vertexes collection v1And v2Carry out following pretreatment:
A) by vertex set v1And v2In point press up time needle sort, the v obtained in step 221And v2Side clockwise is pressed
To sequence;
B) vertex set v1And v2It is satisfied by p1(379043.6582,3277656.1655)=p18(379043.6582,
3277656.1655) and q1(379043.6582,3277656.1655)=q23(379043.6582,3277656.1655), i.e.,
First node is equal to tail node, jump procedure 23.
Step 23:Calculate vertex set v1And v2Public point set T=(379043.6582,3277656.1655),
(379018.7804,3277657.9852), (378995.2579,3277644.7475) ... (379043.6582,
3277656.1655) }, the number of common point is 9, i.e. j=9.And each point is recorded in T in vertex set v1And v2In serial number, Q1=
{ 1,2,3,4,5,6,7,8,18 } and Q2={ 22,21,20,19,18,17,16,15,0 };
Step 24:The segment processing of public point set T, segmentation processing step are as follows:
A) k is set as cyclic variable, k ∈ [1, j-1], initial value 1.It is segmented by taking k=1 as an example:Cycle reads serial number Q1
In the 1st period 1, the 2nd period 2 equally read serial number Q2In the 1st period 22, the 2nd 21 , ∣ 1-2 Qie of ∣=1 ∣ of period
∣=1 22-21, then serial number Q1Middle period 1 and period 2 or serial number Q2Middle period 22 and period 21, corresponding two points belong to same
Section public boundary;
B) cyclic variable k adds 1, continues to execute k=2,3..., 9, is segmented.Segmentation result shows in v1Middle period 1
Same section of public boundary is belonged to 2 points corresponding to 2,2 and 3,3 and 4,4 and 5,5 and 6,6 and 7,7 and 8.Thus, v1Midpoint
Number 1, the point corresponding to 2,3,4,5,6,7,8 belongs to same segment boundary.And v1And v2Only one section of public boundary, public boundary
E1/2={ e1}.That is m'=1 is not executed c), jump procedure 25.
Step 25:Cyclic variable i adds 1, i=2, i≤785, continues to execute i=2,3..., 785, extracts stratigraphic boundary, obtains
To the stratigraphic boundary set E={ E of entire figure layer1/2, E1/58... E417/418, totally 785 (926 sections) stratigraphic boundaries, and with line
The mode of figure layer shows, as shown in Figure 4.Since there are the situation that stratigraphic boundary includes multistage public boundary, such as E16/388={ e1,
e2, including 2 sections of public boundaries, thus stratigraphic boundary E includes 926 sections of stratigraphic boundaries altogether.
(3) it is based on layer data, identifies fault contact relationship
Step 31:Load Mount Lushan tomograph layer data (Fig. 2).And with the fault line F={ f in tomography figure layer1, f2...,
f51Centered on axis, distance parameter 100m (distance parameter d=100m in this example) generates tomography buffering area FB, such as Fig. 5 institutes
Show;
Step 32:Extraction interrupts layer line F={ f with layer data1, f2..., f51Intersection stratum syntople figure
Line set D={ d1, d2..., d142, 142 sides are shared in SARG in this example and are intersected with fault line F, as shown in Figure 6;
Step 33:Search side d in line set Di(i is cyclic variable, i ∈ [1,142], initial value 1), corresponding stratum
Face element SaWith stratal surface element SbBetween stratigraphic boundary Ea/b={ e1, e2..., em}.With side d129For, it extracts corresponding
Stratal surface element S49With stratal surface element S343Between stratigraphic boundary E49/343=e1;
Step 34:Judge fault contact:If Ea/b={ e1, e2..., emCertain segment limit line ejIn tomography buffering area FBIt is interior,
Then ejBoth sides the contact relation of interlayer be fault contact relationship.Continue with stratigraphic boundary e1For, e1In tomography buffering area FB
It is interior, then e1Both sides stratum, that is, Maoshan Mountain group is with five logical groups in e1The contact relation at place is fault contact relationship;
Step 35:Cyclic variable i adds 1, if i≤142, repeats step 33,34, until tomography all in figure layer
Contact extraction finishes.The final result of fault contact extraction is as shown in Figure 7.
(4) occurrence data are based on, identify basic contact relation, complete the making of stratigraphic boundary figure layer.
Step 41:Based on formation lithology information, intrusive body and non-intruding body are distinguished.It is s that result is distinguished in this example7(arteries and veins
Rock:Granitic pegmatite dike), s9(dyke rock:Quartz vein), s393(Guanyinqiao Gneiss Suite:The long gneiss of the ridges Gao Jia potassium), s394(star
Sub- rock group:Ascent temple granulite) it is intrusive body, remaining stratal surface element is non-intruding body.And will distinguish result be recorded it is heavy
Product and volcano stratigraphic unit figure layer, the stratal surface element of metamorphic strata unit figure layer " rock class " attribute in;
Step 42:Based on lithostratic unit's abridged table, formation sequence Table X, and older ordinal number of stratum numerical value is smaller.This reality
Mount Lushan stratigraphic unit abridged table, formation sequence Table X are based in example, wherein that minimum is Pt1Xq (planetesimal rock group, age are oldest), most
Big is Qhl (even country fair group, age are newest);
Step 43:Read i-th boundary line element E in stratigraphic boundary set Ea/b(i is cyclic variable, i ∈ [1, k], just
Value is item number that 1, k is stratigraphic boundary) corresponding to stratal surface element sa、sbLithostratic unit Rocknamea、
Rocknameb.Respectively with the 227th article of boundary line element E in stratigraphic boundary set E in this example85/275Corresponding stratal surface is wanted
Plain s85(clear water group epimere) and s275The lithostratic unit S of (New Port clay hypomere)1q2And Qp3x1, with the 732nd article of boundary line element
E366/394Corresponding stratal surface element s366(Husband's older sister's mud boulder), s394(planetesimal rock group:Ascent temple granulite) lithostratigraphy it is single
Position Qp2D and Pt1For xg, step 44-48 is executed;
Step 44:Read stratal surface element s85(clear water group epimere) and s275" rock class " attribute of (New Port clay hypomere).Ground
Level element s85(clear water group epimere) and s275" the rock class " of (New Port clay hypomere) is non-intruding body, jump procedure 45;And ground
Level element s366" the rock class " of (Husband's older sister's mud boulder) is non-intruding body, s394(planetesimal rock group:Ascent temple granulite) " rock class " be
Intrusive body, then jump procedure 47;
Step 45:Integrate the judgement of contact and unconformity contact.Read lithostratic unit S1q2And Qp3x1In sequence table X
On serial number be respectively 29,58 ∣=29 ≠ 1 , ∣ 29-58, then stratal surface element s85And s275Contact relation be unconformity contact,
Jump procedure 46 carries out next step judgement;
Step 46:Disconformity contacts and angular unconformity contact judges.Stratal surface element s85(clear water group epimere),
s275The contact relation of (New Port clay hypomere) judge in next step as follows:
A) load Mount Lushan occurrence figure layer data (Fig. 2), and read stratal surface element s85Move towards Strike85=" rock stratum is just
Often trend ", tendency Dip85=" SE " and inclination angle Angle85=25 °;Stratal surface element s275Move towards Strike275=" rock stratum is just
Often trend ", tendency Dip275=" SE " and inclination angle Angle275=20 °;
b)Strike85=Strike275And Dip85=Dip275Qie ∣ Angle85- Angle275∣=5≤δ (δ in this example
Take 5 °) it sets up, then stratal surface element s85(clear water group epimere), s275The contact relation of (New Port clay hypomere) is disconformity
Contact.
Step 47:The judgement of Deposit contact and intrusive contact.Read lithostratic unit Qp2D (Husband's older sister's mud boulder) and Pt1xg
(planetesimal rock group:Ascent temple granulite) serial number on sequence table X is respectively 54,1.Due to stratal surface element s394(planetesimal rock
Group:Ascent temple granulite) " rock class " attribute be intrusive body, and x366> x394, then it is Deposit contact;
Step 48:Record the pel number 85 and 275 and disconformity contact relation attribute and stratum of stratum line feature
The pel number 366 of line feature and 394 and Deposit contact attribute of a relation.Circulation step 43-48, identification stratigraphic boundary set E=
{E1/2, E1/58... E417/418In 926 sections of stratigraphic boundary element both sides interlayer contact relation.Recognition result is to integrate to contact
185 sections, 36 sections of disconformity contact, 6 sections of angular unconformity, 663 sections of unconformity contact, 29 sections of fault contact, intrusive contact 4
Section, 3 sections of Deposit contact, and shown in the form of thematic map, as shown in figure 8, completing the making of stratigraphic boundary figure layer.
Claims (5)
1. a kind of automatic generating method of stratigraphic boundary figure layer, which is characterized in that this method comprises the following steps:
(1) it is based on deposition, volcano stratigraphic boundary figure layer and metamorphic strata unit figure layer, hunts out all of its neighbor stratum in figure layer,
Build stratum syntople figure;
(2) public boundary that adjacent earth formations are extracted based on stratum syntople figure, generates stratigraphic boundary;
(3) it is based on tomography figure layer, identifies fault contact relationship;
(4) it is based on occurrence figure layer, identifies basic contact relation, the automation for completing stratigraphic boundary figure layer generates.
2. a kind of automatic generating method of stratigraphic boundary figure layer according to claim 1, which is characterized in that the step
(1) it is as follows:
Step 11:Load deposition, volcano stratigraphic boundary figure layer and metamorphic strata unit figure layer data, and corresponding stratal surface is wanted
Element reads in set S={ si| i=1,2 ..., n }, n is the number of stratal surface element;
Step 12:Read each stratal surface element siVertex set vi={ Pj| j=1,2 ..., m }, wherein PjIt indicates to constitute stratum
Face element siVertex, m be vertex number;And to each stratal surface element siVertex set viCarry out following pretreatment:
A) by vertex set viIn point PjBy up time needle sort;
B) judge each vertex set viIn whether there is P1=Pm, i.e., first node is equal to tail node, and if it exists, then redirects in next step
C), opposite vertexes collection viIt is handled;Otherwise, jump procedure 13;
C) vertex set v is deletediTail node Pm, ensure that head and the tail do not connect, i.e. vi={ Pj| j=1,2 ..., m-1 };
Step 13:The syntople of stratal surface element in figure layer just sentence:If stratal surface element saVertex set vaMinimum
Boundary rectangle RaWith element sbVertex set vbMinimum enclosed rectangle RbThere is intersection point, then jump procedure 14;Otherwise, jump procedure
16;Wherein, a is cyclic variable, a ∈ [1, n-1], initial value 1;B is cyclic variable, b ∈ [a+1, n], initial value a+1;
Step 14:To two stratal surface element saAnd sbSyntople carefully sentenced, detailed process is as follows:
A) vertex set v is calculatedaAnd vbPublic point set T={ p1, p2..., pk, k is the number of common point, and records the midpoints T and exist
Vertex set vaMiddle serial number is denoted as N={ n1, n2..., nk};
If b) number k >=2 at the midpoints public point set T, redirect next step c);Otherwise, jump procedure 16;
C) judge to whether there is continuous period n in common point serial number collection Nt, nt+1Man Zu ∣ nt- nt+1∣=1, and if it exists, i.e.
Level element s1And s2There are syntople, then jump procedure 15;Otherwise, jump procedure 16;
Step 15:Build the side d of stratum syntople figure:Calculate vertex set vaWith vbFocus point CaAnd Cb, point CaAnd CbAs
The vertex of stratum syntople figure, tie point CaAnd CbSide d as stratum syntople figure;
Step 16:If cyclic variable a > n-1, export SARG={ d1, d2... dw, w is stratum syntople figure SARG's
The total number on side;Otherwise, cyclic variable b adds 1, if b≤n at this time, direct jump procedure 13;Otherwise cyclic variable a adds 1, follows
Ring variable b=a+1, then jump procedure 13.
3. a kind of automatic generating method of stratigraphic boundary figure layer according to claim 2, which is characterized in that the step
(2) it is as follows:
Step 21:Read SARG={ d1, d2... dwIn i-th side di, i is cyclic variable, initial value 1;
Step 22:Find out stratal surface element saAnd sbVertex set va={ p1, p2..., pmAnd vb={ q1, q2..., qn, m, n
Respectively stratal surface element sa、sbNumber of vertices, and opposite vertexes collection vaAnd vbCarry out following pretreatment:
A) by vertex set va={ p1, p2..., pmAnd vb={ q1, q2..., qnIn point press up time needle sort;
B) judge vertex set vaAnd vbWhether p is met simultaneously1=pmWith q1=qn, i.e., first node is equal to tail node;If satisfied, then jumping
Go to step 23;Otherwise, to ineligible vertex set vaOr vb, execute next step c);
C) opposite vertexes collection vaOr vbContinue to add first point p1Or q1, ensure end to end;
Step 23:Calculate vertex set vaAnd vbPublic point set T={ t1, t2..., tj, j is the number of common point, and records public affairs
Each point is respectively in vertex set v in concurrent collection TaAnd vbIn serial number, be denoted as Q respectivelyaAnd Qb;
Step 24:The segment processing of public point set T, segmentation processing step are as follows:
A) k is set as cyclic variable, k ∈ [1, j-1], initial value 1;Read serial number QaIn+1 k-th of period i', kth period j';
Equally, serial number Q is readbIn+1 k-th of period ii', kth period jj', if Man Zu ∣ i'-j' ∣=1 Qie ∣ ii'-jj' simultaneously
∣=1, then 2 points corresponding to period i' and period j' and period ii' and period jj' belong to same section of public boundary;
B) cyclic variable k adds 1, if k < j, continues to execute previous step a) until public point set T segmentations terminate, and will belong to successively
The point of same segment boundary constitutes side e, obtains public boundary Ea/b={ e1, e2.., em', m' indicates the hop count of public boundary;Such as
Fruit m'>2, execute next step c);Otherwise jump procedure 25;
C) compare two-by-two successively and constitute side etLast point coordinates tL(x, y) and composition side et+1First point coordinates tF(x, y),
T ∈ [1, m'], if there are tL(x, y)=tF(x, y), then by side etWith et+1Merge, side etWith et+1Belong to same section of common edge
Boundary;
Step 25:Cyclic variable i adds 1, if i≤w, cycle executes step 21-24, obtains the stratigraphic boundary set E of entire figure layer
={ Ea/b, Ec/d... }, wherein Ea/b, Ec/dRespectively represent adjacent stratum saWith sb, adjacent stratum scWith sdBetween stratigraphic boundary.
4. a kind of automatic generating method of stratigraphic boundary figure layer according to claim 3, which is characterized in that the step
(3) it is as follows:
Step 31:Tomograph layer data is loaded, and with the fault line F={ f in tomography figure layer1, f2..., fmCentered on axis,
Tomography buffering area F is generated with distance parameter dB;
Step 32:Extraction interrupts layer line F={ f with layer data1, f2..., fmIntersection stratum syntople figure line set
D={ d1, d2..., dn, n is the item number that stratum syntople figure intersects with fault line F;
Step 33:Search side d in line set DiCorresponding stratal surface element SaWith stratal surface element SbBetween stratigraphic boundary Ea/b
={ e1, e2..., em, wherein i is cyclic variable, i ∈ [1, n], initial value 1;
Step 34:Judge fault contact:If Ea/b={ e1, e2..., emCertain segment limit line ejIn tomography buffering area FBIt is interior, then boundary
Line ejBoth sides the contact relation of interlayer be fault contact relationship;
Step 35:Cyclic variable i adds 1, if i≤n, repeats step 33,34, until fault contact all in figure layer carries
It takes complete.
5. a kind of automatic generating method of stratigraphic boundary figure layer according to claim 4, which is characterized in that the step
(4) it is as follows:
Step 41:Based on formation lithology information, intrusive body and non-intruding body are distinguished;It is heavy that the differentiation result on each stratum is recorded
Product and volcano stratigraphic unit figure layer, the stratal surface element of metamorphic strata unit figure layer " rock class " attribute in;
Step 42:Based on lithostratic unit's abridged table, formation sequence Table X, and older ordinal number of stratum numerical value is smaller;
Step 43:Read i-th boundary line element E in stratigraphic boundary set Ea/bCorresponding stratal surface element saAnd sbRock
Stratigraphic unit RocknameaAnd Rocknameb, wherein i is cyclic variable, and i ∈ [1, k], initial value 1, k is the item of stratigraphic boundary
Number;
Step 44:Read stratal surface element saAnd sb" rock class " attribute;If stratal surface element saAnd sb" rock class " be non-invade
Enter body, then jump procedure 45;Otherwise, jump procedure 47;
Step 45:Integrate the judgement of contact and unconformity contact;Read lithostratic unit RocknameaAnd RocknamebIn sequence
Serial number x on list Xa、xbIf serial number xa、xbIt is adjacent or identical, i.e. xa=xbHuo person ∣ xa- xb∣=1, then stratal surface element sa
And sbContact relation be integrate contact;Otherwise, it is that unconformity contacts, jump procedure 46 carries out next step judgement;
Step 46:Disconformity contacts and angular unconformity contact judges, specific judgment method is as follows:
A) occurrence figure layer data are loaded, and read stratal surface element sa、sbAttitude:Move towards StrikeaAnd Strikeb, incline
To DipaAnd DipbAnd inclination angle AngleaAnd Angleb;If any attitude information can not obtain, stratal surface element saWith
sbContact relation be only judged as unconformity contact, jump procedure 48;
B) judge Strikea=StrikebAnd Dipa=DipbQie ∣ Anglea- AnglebWhether ∣≤δ is true, wherein 0 ° of < δ <
10°;If so, then stratal surface element saAnd sbContact relation be disconformity contact;Otherwise, it is that angular unconformity contacts;
Step 47:The judgement of Deposit contact and intrusive contact;Read lithostratic unit RocknameaAnd RocknamebIn sequence
Serial number x in Table Xa、xb;As stratal surface element sa" rock class " attribute be intrusive body when, if xa≥xb, then it is intrusive contact;It is no
Then, then it is Deposit contact;Equally, as stratal surface element sb" rock class " attribute when being intrusive body, and xa< xb, then connect for intrusion
It touches;Otherwise, then it is Deposit contact;
Step 48:Record the pel number and contact relation attribute of stratum line feature;And circulation step 43-48, until ground stratum boundary
In line set E all boundary line elements both sides interlayer the equal end of identification of contact relation, complete stratigraphic boundary figure layer making.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023762.2A CN105701848B (en) | 2016-01-14 | 2016-01-14 | A kind of automatic generating method of stratigraphic boundary figure layer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610023762.2A CN105701848B (en) | 2016-01-14 | 2016-01-14 | A kind of automatic generating method of stratigraphic boundary figure layer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105701848A CN105701848A (en) | 2016-06-22 |
CN105701848B true CN105701848B (en) | 2018-07-20 |
Family
ID=56227367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610023762.2A Expired - Fee Related CN105701848B (en) | 2016-01-14 | 2016-01-14 | A kind of automatic generating method of stratigraphic boundary figure layer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105701848B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106384016B (en) * | 2016-10-08 | 2019-04-23 | 南京师范大学 | The automatic inversion method of single dip fault |
CN106934357B (en) * | 2017-02-28 | 2020-04-03 | 南京师范大学 | Automatic parallel fault identification method |
CN106897445A (en) * | 2017-03-01 | 2017-06-27 | 武汉科技大学 | It is a kind of that the novel method that operation delimited on border is improved based on ArcGIS |
CN107065009B (en) * | 2017-04-18 | 2018-12-11 | 中国石油大学(华东) | A kind of buried hill oil gas field underground geological map preparation method based on seismic data |
CN107146283B (en) * | 2017-05-10 | 2020-05-19 | 南京师范大学 | Automatic dividing method for layered rock slope types |
CN111583776B (en) * | 2020-04-28 | 2022-02-11 | 南京师范大学 | Method for acquiring development time sequence of invaded rock mass |
CN111681314B (en) * | 2020-05-29 | 2021-04-06 | 中国地质大学(武汉) | Virtual contour constraint-based fourth-system covering layer three-dimensional modeling method |
CN112231423B (en) * | 2020-09-22 | 2024-02-23 | 南京师范大学 | Geological map-based fault structure automatic recovery method and device |
CN113539051B (en) * | 2021-06-24 | 2022-11-25 | 南京师范大学 | Geological map-based method and device for acquiring stratum boundary point-by-point rock stratum attitude |
CN115544196A (en) * | 2022-10-12 | 2022-12-30 | 兰州天瑞测绘工程有限公司 | Automatic generation method and system for rights boundary line and trend description |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030206655A1 (en) * | 2002-05-06 | 2003-11-06 | Tun-Chi Huang | Method for determining an adjacency relation |
CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
CN102651143A (en) * | 2011-02-25 | 2012-08-29 | 中冶沈勘工程技术有限公司 | Automatically generating method of pinchout geological body three-dimensional complex profile |
CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
-
2016
- 2016-01-14 CN CN201610023762.2A patent/CN105701848B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030206655A1 (en) * | 2002-05-06 | 2003-11-06 | Tun-Chi Huang | Method for determining an adjacency relation |
CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
CN102651143A (en) * | 2011-02-25 | 2012-08-29 | 中冶沈勘工程技术有限公司 | Automatically generating method of pinchout geological body three-dimensional complex profile |
CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
Also Published As
Publication number | Publication date |
---|---|
CN105701848A (en) | 2016-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105701848B (en) | A kind of automatic generating method of stratigraphic boundary figure layer | |
CN103424773B (en) | A kind of palaeogeomorphology based on horizon flattening method | |
CN109903383B (en) | Method for accurately positioning coal mining machine in three-dimensional model of working face coal seam | |
US20120166160A1 (en) | Block model constructing method for complex geological structures | |
CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
CN103969683B (en) | Method for picking position faces in batched mode based on constraint in three-dimensional seismic interpretation | |
CN109658431A (en) | Rock mass point cloud plane extracting method based on region growing | |
CN109763814B (en) | Stratum matching visual analysis method based on multi-dimensional logging data | |
CN104504047A (en) | Estimation system of reserve of solid mineral resources based on two-tier database | |
CN107766978A (en) | The intelligent optimization method of skewed pattern | |
CN109859312A (en) | A kind of fining three-dimensional geological model modeling method based on BIM technology | |
CN108957530B (en) | A kind of crack automatic testing method based on Acceleration Algorithm in Seismic Coherence Cube slice | |
CN115436993A (en) | Seismic signal P wave detection method based on PPNet network | |
CN108986213A (en) | A kind of three dimensional contour line method based on stacking technology | |
CN103824330A (en) | Method for building ore body middle-section layered graph and three-dimensional model | |
CN116721227A (en) | Automatic modeling method for three-dimensional geologic model of complex geologic body | |
CN108873070B (en) | Method for recovering paleogenesis of stratum in basin | |
CN112581558B (en) | Model construction method and system for intrusion structure in cut geological section | |
CN112734926B (en) | Automatic generation method of cut geological section for loose layer coverage area | |
CN112507615B (en) | Intelligent identification and visualization method for lithofacies of continental tight reservoir | |
CN104793263A (en) | Well logging curve automatic tiered value extraction and evaluation method | |
CN113433589A (en) | Weathered crust reservoir bottom interface identification method based on mathematical statistics | |
CN117372643B (en) | Method for constructing three-dimensional stratum model layer by layer from top to bottom based on drilling data | |
CN117251739B (en) | Construction method of sample set for sandstone uranium resource quantitative prediction | |
Li et al. | Automatic Geological Map Generalization for special geobody–A Case Study on Dyke |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180720 Termination date: 20210114 |
|
CF01 | Termination of patent right due to non-payment of annual fee |