CN106384016B - The automatic inversion method of single dip fault - Google Patents
The automatic inversion method of single dip fault Download PDFInfo
- Publication number
- CN106384016B CN106384016B CN201610876852.6A CN201610876852A CN106384016B CN 106384016 B CN106384016 B CN 106384016B CN 201610876852 A CN201610876852 A CN 201610876852A CN 106384016 B CN106384016 B CN 106384016B
- Authority
- CN
- China
- Prior art keywords
- point
- intersection
- midpoint
- coordinate value
- interu1
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a kind of automatic inversion methods of single dip fault.Based on geologic map, automatic inversion is carried out to single dip fault.Mainly include the following steps: 1) to calculate separately the intersection that upper and lower disk corresponds to rock stratum and fault line, with this calculating benchmark point;2) upper and lower disk formation coordinate collection is corresponded to decompose;3) it determines datum mark offset, the coordinate set after decomposition is moved;4) merge the upper and lower Pan Duiyingzuobiaoji after movement.The method achieve the automatic inversions of the single dip fault based on geologic map, and algorithm is clear, practical.
Description
Technical field
The invention belongs to geographical information technology application fields, and in particular to a kind of single dip fault based on geologic map
Automatic inversion method.
Background technique
Architectonic inverting refer to based on tectonic landform type along with other conditions for constraint to deduce rock stratum original
The process of pattern.Tomography occurs obviously to be displaced and shape as one of most important construction of the earth's crust, by rock stratum or rock mass along the plane of disruption
At often will form rift valley and cliff.Inverting, the original of rock stratum when simulating and restoring not to be broken are carried out for Fault Evolution process
Beginning pattern has certain research significance for geological hazards prediction, geological resource exploitation, environmental geological protection etc. and applies valence
Value.
The inverting application of tomography is mainly reflected in following three aspects: 1) based on geologic map.Judge that disk body includes in tomography
The consistency of rock stratum, and then the rock stratum being subjected to displacement is merged, ultimately form the original geologic map after removing fault line;2)
Based on log sheet.The discrete point of geological stratum can be restored to the original shape for being deformed to rock stratum according to the backward of fault rupture
Looks;3) it is based on threedimensional model.Main modification tomography entity built-in attribute data, the fault tectonic of deformation is reverted to completely
Matter entity.At present, on the one hand, involved by related application has, but related work is mainly had been manually done by geology expert, is put into
Greatly, inefficient is badly in need of carrying out the automatic inversion technique study of tomography;On the other hand, geology diagram data is due to from a wealth of sources, covering
The features such as range is wide, using less-restrictive, starts with from geologic map, carries out tomographic inversion research, has preferable feasibility and reality
The property used.
Summary of the invention
It is an object of the present invention to overcome defect of the existing technology, a kind of the automatic of single dip fault is proposed
The automatic inversion of single dip fault may be implemented based on geologic map for inversion method.
Main thought is based on the method for the present invention: 1) intersection that upper and lower disk corresponds to rock stratum and fault line is calculated separately,
With this calculating benchmark point;2) upper and lower disk formation coordinate collection is corresponded to decompose;3) datum mark offset is determined, after decomposition
Coordinate set moved;4) merge the upper and lower Pan Duiyingzuobiaoji after movement.
In order to achieve the above-mentioned object of the invention, the technical solution proposed is the present invention: a kind of single dip fault it is automatic
Inversion method, the specific steps are as follows:
(1) data and calculating benchmark point are obtained
Step 1.1: loading the rock stratum face figure layer and fault line line chart layer data of vector format, respectively obtain hanging layer face
Elements combination UpStra={ usi| i=1,2 ..., n } and lower wall rock stratum face elements combination DownStra={ dsi| i=1,2 ...,
N } and the fault line point set FLine={ fl that includesi| i=1,2 ..., m }.What general line feature was made of multiple points,
Referred to herein as fault line as line feature, it includes point set.
Wherein, usiIt is the hanging layer face element of tomography, dsiIt is the lower wall rock stratum face element of tomography, fliIt is the disconnected of tomography
The point that layer line includes, n are the number of rock stratum face element, and m is the number for the point that fault line includes.
Step 1.2: circulation takes the face element us in set UpStra and DownStraiAnd dsi.If usiWith dsiCorresponding rock
The lithostratic unit of layer is consistent, thens follow the steps 1.3;Otherwise, this step is continued to execute.
Step 1.3: obtaining face element us respectivelyiAnd dsiPoint set UsVi={ uij| j=1,2 ..., a } and DsVi={ dij
| j=1,2 ..., b }, remove point set UsViAnd DsViIn it is respective repeat point.
Wherein, uijIt is composition face element usiJ-th point, dijIt is composition face element dsiJ-th point, a be composition face
Element usiPoint number, b be composition face element dsiPoint number.
Step 1.4: calculating intersection.Calculate separately point set UsViWith the intersection UsInter of point set FLinei={ utj| j=1,
2 ..., p }, point set DsViWith the intersection DsInter of point set FLinei={ dtj| j=1,2 ... q }.
Wherein, utjIt is point set UsViWith the common point of point set FLine, dtjIt is point set DsViIt is public with point set FLine
Point, p are intersection UsInteriThe number at midpoint, q are intersection DsInteriThe number at midpoint.
To intersection UsInteriIn point element 1 be ranked up by rule, generate new intersection UsInter'i={ ut'j|j
=1,2 ..., p }.Wherein, Point, SubPoint and SubPoint' respectively represent point set UsVi, intersection UsInteriAnd intersection
UsInter'i;To intersection DsInteriIn point element by rule 1 sort, generate new intersection DsInter'i={ dt'j| j=
1,2,…q}.Wherein, Point, SubPoint and SubPoint' respectively represent point set DsVi, intersection DsInteriAnd intersection
DsInter'i。
Rule 1: it sets up an office and collects Point={ ptj| j=1,2 ... sum }, point set SubPoint={ sptj| j=1,2 ...
sub}.Wherein, sum is the number at the midpoint point set Point, and sub is the number at the midpoint point set SubPoint, andThe step of sequence, is as follows:
1) point set Point is traversed, when finding point ptjMeet condition ptj-1∈ SubPoint andWhen, it jumps
It recycles out, and records point ptj-1;
2) point pt is searchedj-1Position index in point set SubPoint;
3) point set SubPoint is traversed, the point for being wherein designated as index+1 to sub down is added in point set SubPoint';Again
Point set SubPoint' is added in secondary traversal point set SubPoint, the point that 0 to index is designated as by under.
Step 1.5: calculating benchmark point coordinate.For intersection UsInter'i、DsInter'i, respectively according to formula (1), (2)
Calculate two datum mark DaP1 (x1,y1)、DaP2(x2,y2) coordinate value;
Wherein, (ut'1.x,ut'1.y)、(ut'p.x,ut'pIt .y) is set UsInter' respectivelyiThe coordinate of middle head and the tail point
Value, (dt'1.x,dt'1.y)、(dt'q.x,dt'qIt .y) is set DsInter' respectivelyiThe coordinate value of middle head and the tail point.
(2) boundary of rock stratum is decomposed
Step 2.1: to point set UsViIt is decomposed.It is broken down into following several parts:
1) intersection UsInter'i;
2)UsViWith the intersection InterU1={ u1 of left side adjacent strataj| j=1,2 ..., c }, c is in intersection InterU1
The number of point, and InterU1 includes point ut'p;
3)UsViWith the intersection InterU2={ u2 of right side adjacent strataj| j=1,2 ..., d }, d is in intersection InterU2
The number of point, and InterU2 includes point ut'1;
4)UsViWith UsInter'i, InterU1 and InterU2 difference set SubU={ suj| j=1,2 ..., e }, e is poor
Collect the number at the midpoint SubU.
Step 2.2: the point element in set InterU1, InterU2 and SubU being ranked up by rule 1 respectively.Its
In, to three above set, Point represents point set UsVi, SubPoint respectively represents themselves.
Step 2.3: to point set DsViIt is decomposed.It is broken down into following several parts:
1) intersection DsInter'i;
2)DsViWith the intersection InterD1={ d1 of left side adjacent strataj| j=1,2 ..., f }, f is in intersection InterD1
The number of point, and InterD1 includes point dt'q;
3)DsViWith the intersection InterD2={ d2 of right side adjacent strataj| j=1,2 ..., g }, g is in intersection InterD2
The number of point, and InterD2 includes point dt'1;
4)DsViWith DsInter'i, InterD1 and InterD2 difference set SubD={ sdj| j=1,2 ..., h }, h is poor
Collect the number at the midpoint SubD.
Step 2.4: the point element in set InterD1, InterD2 and SubD being ranked up by rule 1 respectively.Its
In, to three above set, Point represents DsVi, SubPoint respectively represents themselves.
(3) mobile rock stratum
Step 3.1: the coordinate at the mobile midpoint set InterU1:
A) according to each point in formula (3) calculating InterU1 in mean deviation the amount aveX and aveY in the direction x, y:
Wherein, (u.x, u.y) is the coordinate value of first point in set InterU1, and (x, y) is the coordinate value of datum mark DaP1,
k1For the number at the midpoint InterU1.
B) set InterU1 is traversed, modifies u1 according to formula (4)jCoordinate value, obtain it is mobile after set
InterU1'={ u1'j| j=1,2 ..., c }, c is the number at the midpoint set InterU1';
Wherein, (vj.x,vjIt .y) is current point u1jCoordinate value, (v'j.x,v'jIt .y) is point u1'jCoordinate value, k2For
The number at the midpoint InterU1.
Step 3.2: the coordinate at the mobile midpoint set InterU2:
C) according to each point in formula (3) calculating InterU2 in mean deviation the amount aveX and aveY in the direction x, y;Herein,
(u.x, u.y) is the coordinate value of tail point in set InterU2, and (x, y) is the coordinate value of datum mark DaP2, k1For in InterU2
The number of point.
D) set InterU2 is traversed, modifies u2 according to formula (5)jCoordinate value, obtain it is mobile after set
InterU2'={ u2'j| j=1,2 ..., d }, d is the number at the midpoint set InterU2';
Wherein, (wj.x,wjIt .y) is current point u2jCoordinate value, (w'j.x,w'jIt .y) is point u2'jCoordinate value.
Step 3.3: merging set.Set InterU1', SubU and InterU2' are successively traversed, the point in set is added
To set NewUsVi={ nuij| j=1,2 ..., s } in, s is set NewUsViThe number at midpoint, and s=c+e+d.
Step 3.4: the coordinate at the mobile midpoint set InterD1:
E) according to each point in formula (3) calculating InterD1 in the mean deviation amount (aveX, aveY) in the direction x, y;This
Place, (u.x, u.y) is the coordinate value of first point in set InterD1, and (x, y) is the coordinate value of datum mark DaP2, k1For InterD1
The number at midpoint.
F) set InterD1 is traversed, modifies d1 according to formula (4)jCoordinate value, obtain it is mobile after set
InterD1'={ d1'j| j=1,2 ..., f }, f is the number at the midpoint set InterD1';Herein, (vj.x,vjIt .y is) current
Point d1jCoordinate value, (v'j.x,v'jIt .y) is point d1'jCoordinate value, k2For the number at the midpoint InterD1.
Step 3.5: the coordinate at the mobile midpoint set InterD2:
G) according to each point in formula (3) calculating InterD2 in the mean deviation amount (aveX, aveY) in the direction x, y;This
Place, (u.x, u.y) is the coordinate value of tail point in set InterD2, and (x, y) is the coordinate value of datum mark DaP1, k1For InterD2
The number at midpoint.
H) set InterD2 is traversed, modifies d2 according to formula (5)jCoordinate value, obtain it is mobile after set
InterD2'={ d2'j| j=1,2 ..., g }, g is the number at the midpoint set InterD2';Herein, (wj.x,wjIt .y is) current
Point d2jCoordinate value, (w'j.x,w'jIt .y) is point d2'jCoordinate value.
Step 3.6: merging set.Set InterD1', SubD and InterD2' are successively traversed, the point in set is added
To set NewDsVi={ ndij| j=1,2 ..., t } in, t is the number at the midpoint set InterD2', and t=f+h+g.
(4) merge rock stratum
Step 4.1: for set NewUsV obtained in step 3.3i, it is traversed, point is added to set
StraViIn.
Step 4.2: to set NewDsV obtained in step 3.6iIt is traversed, in addition to head and the tail are put, successively by other points
It is added to set StraViIn, obtain set StraVi={ svj| j=1,2 ..., r }, r is set StraViThe number at midpoint,
And r=s+t-2.
Step 4.3: circulation executes step 1.2-4.2, until upper and lower disk corresponds to rock stratum merging and finishes, after being merged
Rock stratum face elements combination SumStra={ StraV1,StraV2,…,StraVz, z is the number of rock stratum face element, generates and protects
Face figure layer data after depositing inverting.
The method of the present invention carries out inverting to single dip fault based on geologic map, and realize based on geologic map single inclines
To the automatic inversion of tomography, algorithm is clear, practical.
Detailed description of the invention
Fig. 1 the method for the present invention treatment process schematic diagram
Wherein, scheme (a) calculating benchmark point, figure (b) decomposes the boundary of rock stratum, the mobile rock stratum of figure (c), and figure (d) merges rock stratum.
Fig. 2 the method for the present invention flow chart.
Fig. 3 experimental data source figure.
Fig. 4 experimental data.
Experimental data after Fig. 5 inverting.
Data investigation figure after Fig. 6 initial data and inverting.
Specific implementation method
Below with reference to embodiment and attached drawing, invention is further described in detail.
Embodiment selects the data in 1:100000 regional geologic map (Fig. 3) comprising pregnance tomography.Below with reference to
Attached drawing, rock stratum face figure layer and fault line line chart layer (Fig. 4) based on shp format, the automatic inversion for describing single dip fault are real
Example, to further illustrate effect of the invention.
(1) data and calculating benchmark point are obtained
Step 1.1: loading the rock stratum face figure layer and fault line line chart layer data of vector format, respectively obtaining disk includes rock
Level elements combination UpStra={ usi| i=1,2 ..., 7 and lower wall include rock stratum face elements combination DownStra={ dsi|i
=1,2 ..., 7 and corresponding fault line include point coordinate set FLine=(147692.772656289,
145470.055976739),(147664.719420046,145498.197935516),…,(146785.737324302,
146304.846020359) (totally 40 points).
Step 1.2: circulation takes the face element us in set UpStra and DownStra1And ds1, us1With ds1Corresponding rock
Layer is consistent, executes step 13.
Step 1.3: obtaining face element us respectively1And ds1Vertex set UsV1=(147295.096583363,
146848.889520708),(147318.529740145,146835.520505232),…,(147295.096583363,
146848.889520708) and DsV1=(146873.077847646,146182.502746869),
(146892.673112575,146158.440067434),…,(146873.077847646,146182.502746869)}。
Remove vertex set UsV1And DsV1In it is respective repeat point, obtain UsV1=(147295.096583363,
146848.889520708),(147318.529740145,146835.520505232),…,(147288.294750398,
146823.391549261) (totally 38 points), DsV1=(146873.077847646,146182.502746869),
(146892.673112575,146158.440067434),…,(146859.686997208,146152.173879583)}
(totally 31 points).
Step 1.4: calculating intersection.Calculate separately vertex set UsV1And DsV1With an intersection of coordinate set FLine
UsInter1=(146964.01018912,146093.829101202), (146929.485052078,
..., 146122.412107123) (146866.504650348,146190.574530541) } (totally 5 points), DsInter1=
{(146873.077847646,146182.502746869),(146892.673112575,146158.440067434),…,
(146973.571382636,146087.296876405) } (totally 5 points).Again respectively to intersection UsInter1And DsInter1In
Point element carry out segmentation sequence, generate new intersection UsInter'1=(146964.01018912,
146093.829101202),(146929.485052078,146122.412107123),…,(146866.504650348,
146190.574530541)}、DsInter'1=(146873.007847646,146182.502746869), (146892.6
73112575.146158.440067434),…,(146973.571382636,146087.296876405)}。
Step 1.5: calculating benchmark point coordinate.For intersection UsInter'1、DsInter'1, respectively according to formula (1), (2)
The coordinate value for calculating two datum marks is respectively DaP1 (146869.791248997,146186.538638705), DaP2
(146968.790785878,146090.562988803)。
(2) boundary of rock stratum is decomposed
Step 2.1: opposite vertexes collection UsV1It is decomposed.It is broken down into following several parts:
1) intersection UsInter'1;
2)UsV1With the intersection InterU1=of left side adjacent strata (147295.096583363,
146848.889520708),(146866.504650348,146190.574530541),…,(147288.294750398,
146823.391549261) } (totally 18 points), and InterU1 include point (146866.504650348,
146190.574530541);
3)UsV1With the intersection InterU2=of right side adjacent strata (147400.278282688,
146790.126222549),(147374.129303738,146738.348326925),…,(146964.01018912,
146093.829101202) } (totally 15 points), and InterU2 include point (146964.01018912,
146093.829101202);
4)UsV1With UsInter'1, InterU1 and InterU2 difference set SubU=(147318.529740145,
146835.520505232), (147349.98018661,146817.577489293) (totally 2 points).
Step 2.2: segmentation sequence being carried out to the point element in set InterU1, InterU2 and SubU respectively, is obtained
InterU1=(146866.504650348,146190.574530541), (146882.747013823,
146220.30859269),…,(147295.096583363,146848.889520708)};In set InterU2 and SubU
Point element sequence do not change.
Step 2.3: opposite vertexes collection DsV1It is decomposed.It is broken down into following several parts:
1) intersection DsInter'1;
2)DsV1With the intersection InterD1=of left side adjacent strata (146973.571382636,
146087.296876405),(146955.214877593,146053.760458872),…,(146630.377275767,
145595.372683169) } (totally 11 points), and InterD1 include point (146973.571382636,
146087.296876405);
3)DsV1With the intersection InterD2=of right side adjacent strata (146873.077847646,
146182.502746869),(146554.673697009,145675.404883069),…,(146859.686997208,
146152.173879583) } (totally 13 points), and InterD2 include point (146873.007847646,
146182.502746869);
4)DsV1With DsInter'1, InterD1 and InterD2 difference set SubD=(146611.302175187,
145615.506132667),(146593.977927179,145633.711592918),…,(146562.43399724,
145667.134882314) (totally 4 points).
Step 2.4: segmentation sequence being carried out to the point element in set InterD1, InterD2 and SubD respectively, is obtained
InterD2=(146554.673697009,145675.404883069), (146575.690500073,
145701.901789373),…,(146873.007847646,146182.502746869);In set InterD1 and SubD
Point element sequence do not change.
(3) mobile rock stratum
Step 3.1: the coordinate at the mobile midpoint set InterU1:
A) according to formula (3) calculate InterU1 in each point the direction x, y mean deviation amount (-
0.19332933229366867,0.23740540211849079);
B) set InterU1 is traversed, modifies u1 according to formula (4)jCoordinate value, obtain it is mobile after set
InterU1'=(146869.791248997,146186.538638705), (146885.84028314,
146216.510106256) ..., (147295.096583363,146848.889520708) } (totally 18 points);
Step 3.2: the coordinate at the mobile midpoint set InterU2;
C) according to formula (3) calculate InterU2 in each point the direction x, y mean deviation amount (-
0.3414711970004386,0.23329374275012274);
D) set InterU2 is traversed, modifies u2 according to formula (5)jCoordinate value, obtain it is mobile after set
InterU2'=(147400.278282688,146790.126222549), (147374.470774935,
146738.115033182) ..., (146968.790785878,146090.562988803) } (totally 15 points).
Step 3.3: merging set.Set InterU1', SubU and InterU2' are successively traversed, the point in set is added
To set NewUsV1=(146869.791248997,146186.538638705), (146885.84028314,
146216.510106256) ..., (146968.790785878,146090.562988803) } (totally 35 points);
Step 3.4: the coordinate at the mobile midpoint set InterD1;
E) according to each point in formula (3) calculating InterD1 in the mean deviation amount in the direction x, y
(0.47805967580061404,0.23740540211849079);
F) set InterD1 is traversed, modifies d1 according to formula (4)jCoordinate value, obtain it is mobile after set
InterD1'=(146968.790785878,146090.562988803), (146950.912340511,
146056.699960031) ..., (146630.377275767,145595.372683169) (totally 11 points);
Step 35: the coordinate at the mobile midpoint set InterD2;
G) according to each point in formula (3) calculating InterD2 in the mean deviation amount in the direction x, y
(0.27388322074936394);
H) set InterD2 is traversed, modifies d2 according to formula (5)jCoordinate value, obtain it is mobile after set
InterD2'=(146554.673697009,145675.404883069), (146575.416616852,
145702.238113693) ..., (146869.791248997,146186.538638705) } (totally 13 points).
Step 3.6: merging set.Set InterD1', SubD and InterD2' are successively traversed, the point in set is added
To set
NewDsV1=(146968.790785878,146090.562988803), (146950.912340511,
146056.699960031) ..., (146869.791248997,146186.538638705) } (totally 28 points).
(4) merge rock stratum
Step 4.1: for set NewUsV obtained in step 331, it is traversed, point is added to set
StraV1In.
Step 4.2: to set NewDsV obtained in step 361It is traversed, in addition to head and the tail are put, is successively added other points
It is added to set StraV1In, final StraV1=(146869.791248997,146186.538638705),
(146885.84028314,146216.510106256) ..., (146856.67428178,146155.873447099) } (altogether
61 points).
Step 4.3: circulation executes step 1.2-4.2, until upper and lower disk corresponds to rock stratum merging and finishes, after being merged
Rock stratum face elements combination SumStra={ StraV1,StraV2,…StraV7, it generates and saves the face figure layer data after inverting
(Fig. 5).
Claims (1)
1. the automatic inversion method of single dip fault, the specific steps are as follows:
Step 1 obtains data and calculating benchmark point
Step 1.1: loading the rock stratum face figure layer and fault line line chart layer data of vector format, respectively obtain hanging layer face element
Set UpStra={ usi| i=1,2 ..., n } and lower wall rock stratum face elements combination DownStra={ dsi| i=1,2 ..., n } with
And the point set FLine={ fl that fault line includesj| j=1,2 ..., m };
Wherein, usiIt is the hanging layer face element of tomography, dsiIt is the lower wall rock stratum face element of tomography, fljIt is the fault line of tomography
The point for including, n are the number of rock stratum face element, and m is the number for the point that fault line includes;
Step 1.2: circulation takes the face element us in set UpStra and DownStraiAnd dsi;If usiWith dsiThe rock of corresponding rock stratum
Stone stratigraphic unit is consistent, thens follow the steps 1.3;Otherwise, this step is continued to execute;
Step 1.3: obtaining face element us respectivelyiAnd dsiPoint set UsVi={ uik| k=1,2 ..., a } and DsVi={ dil| l=
1,2,…,b};Remove vertex set UsViAnd DsViIn it is respective repeat point;
Wherein, uikIt is composition face element usiK-th point, dilIt is composition face element dsiFirst point, a be composition face element
usiPoint number, b be composition face element dsiPoint number;
Step 1.4: calculating intersection
Calculate separately point set UsViWith the intersection UsInter of point set FLinei={ utr| r=1,2 ..., p }, point set DsViWith point set
The intersection DsInter of FLinei={ dts| s=1,2 ... q };
Wherein, utrIt is point set UsViWith r-th of common point of point set FLine, dtsIt is point set DsViWith s-th of point set FLine
Common point, p are intersection UsInteriThe number at midpoint, q are intersection DsInterjThe number at midpoint;
Setting rule 1: it sets up an office and collects Point={ pte| e=1,2 ..., sum }, point set SubPoint={ sptf| f=1,2 ...,
Sub }, wherein sum is the number at the midpoint point set Point, and sub is the number at the midpoint point set SubPoint, andThe step of sequence, is as follows:
1) point set Point is traversed, when finding point pteMeet condition pte-1∈ SubPoint andWhen, it jumps out and follows
Ring, and record point pte-1;
2) point pt is searchede-1Position index in point set SubPoint;
3) point set SubPoint is traversed, the point for being wherein designated as index+1 to sub down is added in point set SubPoint';Again time
Point set SubPoint is gone through, point set SubPoint' is added in the point that 0 to index is designated as by under;
To intersection UsInteriIn point element 1 be ranked up by rule, generate new intersection UsInter'i={ ut'r| r=1,
2 ..., p }, wherein Point, SubPoint and SubPoint' respectively represent point set UsVi, intersection UsInteriAnd intersection
UsInter'i;To intersection DsInterjIn point element by rule 1 sort, generate new intersection DsInter'j={ dt's| s=
1,2 ... q }, wherein Point, SubPoint and SubPoint' respectively represent point set DsVj, intersection DsInterjAnd intersection
DsInter'j;
Step 1.5: calculating benchmark point coordinate
For intersection UsInter'i、DsInter'j, two datum mark DaP1 (x are calculated according to formula (1), (2) respectively1,y1)、
DaP2(x2,y2) coordinate value:
Wherein, (ut'1.x,ut'1.y)、(ut'p.x,ut'pIt .y) is set UsInter' respectivelyiThe coordinate value of middle head and the tail point,
(dt'1.x,dt'1.y)、(dt'q.x,dt'qIt .y) is set DsInter' respectivelyjThe coordinate value of middle head and the tail point;
Step 2, the boundary for decomposing rock stratum
Step 2.1: opposite vertexes collection UsViIt is decomposed;It is broken down into following several parts:
1) intersection UsInter'i;
2)UsViWith the intersection InterU1={ u1 of left side adjacent stratag| g=1,2 ..., c }, c is the midpoint intersection InterU1
Number, and InterU1 includes point ut'p;
3)UsViWith the intersection InterU2={ u2 of right side adjacent stratah| h=1,2 ..., d }, d is the midpoint intersection InterU2
Number, and InterU2 includes point ut'1;
4)UsViWith UsInter'i, InterU1 and InterU2 difference set SubU={ suo| o=1,2 ..., ω }, ω is difference set
The number at the midpoint SubU;
Step 2.2: the point element in set InterU1, InterU2 and SubU being ranked up by rule 1 respectively, wherein right
InterU1, InterU2 and SubU tri- set, Point represent point set UsVi, SubPoint respectively represents themselves;
Step 2.3: opposite vertexes collection DsVjIt is decomposed;It is broken down into following several parts:
1) intersection DsInter'j;
2)DsVjWith the intersection of left side adjacent strata For the number at the midpoint intersection InterD1,
And InterD1 includes point dt'q;
3)DsVjWith the intersection of right side adjacent strata For the number at the midpoint intersection InterD2,
And InterD2 includes point dt'1;
4)DsVjWith DsInter'j, InterD1 and InterD2 difference set SubD={ sdw| w=1,2 ..., sw }, sw is difference set
The number at the midpoint SubD;
Step 2.4: the point element in set InterD1, InterD2 and SubD being ranked up by rule 1 respectively, wherein right
Three above set, Point represent DsVj, SubPoint respectively represents themselves;
Step 3, mobile rock stratum
Step 3.1: the coordinate at the mobile midpoint set InterU1:
A) according to each point in formula (3) calculating InterU1 in mean deviation the amount aveX and aveY in the direction x, y:
Wherein, (u.x, u.y) is the coordinate value of first point in set InterU1, and (x, y) is the coordinate value of datum mark DaP1, k1For
The number c at the midpoint InterU1;
B) set InterU1 is traversed, modifies u1 according to formula (4)gCoordinate value, obtain it is mobile after set For the number at the midpoint set InterU1';
Wherein, (vj.x,vjIt .y) is current point u1gCoordinate value, (v'j.x,v'jIt .y) is point u1'gCoordinate value, k2For
The number at the midpoint InterU1;
Step 3.2: the coordinate at the mobile midpoint set InterU2:
C) according to each point in formula (3) calculating InterU2 in mean deviation the amount aveX and aveY in the direction x, y;Herein,
(u.x, u.y) is the coordinate value of tail point in set InterU2, and (x, y) is the coordinate value of datum mark DaP2, k1For in InterU2
The number of point;
D) set InterU2 is traversed, modifies u2 according to formula (5)hCoordinate value, obtain it is mobile after set For the number at the midpoint set InterU2';
Wherein, (wj.x,wjIt .y) is current point u2hCoordinate value, (w'j.x,w'jIt .y) is point u2'ajCoordinate value;
Step 3.3: merging set
Set InterU1', SubU and InterU2' are successively traversed, the point in set is added to set NewUsVi={ nuiw|w
=1,2 ..., s a } in, sa is set NewUsViThe number at midpoint, and
Step 3.4: the coordinate at the mobile midpoint set InterD1:
E) according to each point in formula (3) calculating InterD1 in mean deviation the amount aveX and aveY in the direction x, y;Herein,
(u.x, u.y) is the coordinate value of first point in set InterD1, and (x, y) is the coordinate value of datum mark DaP2, k1For in InterD1
The number of point;
F) set InterD1 is traversed, modifies d1 according to formula (4)uCoordinate value, obtain it is mobile after set InterD1'=
{d1'u| u=1,2 ..., ε }, ε is the number at the midpoint set InterD1';Herein, (vj.x,vjIt .y) is current point d1uCoordinate
Value, (v'j.x,v'jIt .y) is point d1'uCoordinate value, k2For the number at the midpoint InterD1;
Step 3.5: the coordinate at the mobile midpoint set InterD2:
G) according to each point in formula (3) calculating InterD2 in mean deviation the amount aveX and aveY in the direction x, y;Herein,
(u.x, u.y) is the coordinate value of tail point in set InterD2, and (x, y) is the coordinate value of datum mark DaP1, k1For in InterD2
The number of point;
H) set InterD2 is traversed, modifies d2 according to formula (5)vCoordinate value, obtain it is mobile after set InterD2'=
{d2'v| v=1,2 ..., β }, β is the number at the midpoint set InterD2';Herein, (wj.x,wjIt .y) is current point d2vCoordinate
Value, (w'j.x,w'jIt .y) is point d2'vCoordinate value;
Step 3.6: merging set
Set InterD1', SubD and InterD2' are successively traversed, the point in set is added to set NewDsVj={ ndia|a
=1,2 ..., st } in, st is set NewDsVjThe number at midpoint, and st=ε+sw+ β;
Step 4 merges rock stratum
Step 4.1: for set NewUsV obtained in step 3.3i, it is traversed, point is added to set StraVi
In;
Step 4.2: to set NewDsV obtained in step 3.6jIt is traversed, in addition to head and the tail are put, is successively added to other points
Set StraViIn, obtain set StraVi={ svb| b=1,2 ..., sr }, sr is set StraViThe number at midpoint, and sr
=sa+st-2;
Step 4.3: circulation executes step 1.2-4.2, until upper and lower disk corresponds to rock stratum merging and finishes, the rock stratum after being merged
Face elements combination SumStra={ StraV1,StraV2,…,StraVz, z is the number of rock stratum face element, generates and saves anti-
Face figure layer data after drilling.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610876852.6A CN106384016B (en) | 2016-10-08 | 2016-10-08 | The automatic inversion method of single dip fault |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610876852.6A CN106384016B (en) | 2016-10-08 | 2016-10-08 | The automatic inversion method of single dip fault |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106384016A CN106384016A (en) | 2017-02-08 |
CN106384016B true CN106384016B (en) | 2019-04-23 |
Family
ID=57936229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610876852.6A Active CN106384016B (en) | 2016-10-08 | 2016-10-08 | The automatic inversion method of single dip fault |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106384016B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112231423B (en) * | 2020-09-22 | 2024-02-23 | 南京师范大学 | Geological map-based fault structure automatic recovery method and device |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329406A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Broken layer accurate homing method for improving structure diagram precision |
CN104715036A (en) * | 2015-03-17 | 2015-06-17 | 南京师范大学 | GIS surface element geometrical information quantity calculation method based on figure complexity |
CN104915991A (en) * | 2015-06-04 | 2015-09-16 | 南京师范大学 | Spatial subdivision method for intersecting linear geographic elements |
CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
CN105894587A (en) * | 2016-04-01 | 2016-08-24 | 南京师范大学 | Ridge line and valley line filtering method based on rule constraints |
-
2016
- 2016-10-08 CN CN201610876852.6A patent/CN106384016B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329406A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Broken layer accurate homing method for improving structure diagram precision |
CN104715036A (en) * | 2015-03-17 | 2015-06-17 | 南京师范大学 | GIS surface element geometrical information quantity calculation method based on figure complexity |
CN104915991A (en) * | 2015-06-04 | 2015-09-16 | 南京师范大学 | Spatial subdivision method for intersecting linear geographic elements |
CN105225273A (en) * | 2015-09-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Method and the device of FAULT MODEL is set up based on space tomography polygon |
CN105701848A (en) * | 2016-01-14 | 2016-06-22 | 南京师范大学 | Automatic generation method of stratum boundary map layer |
CN105894587A (en) * | 2016-04-01 | 2016-08-24 | 南京师范大学 | Ridge line and valley line filtering method based on rule constraints |
Non-Patent Citations (6)
Title |
---|
Study on the Method of Cartographic Visualization Quality Automated Testing;Anbo Li;《2010 18th International Conference on Geoinformatics》;20100620;第1-6页 |
Three-dimensional structure along the inverted Palaeoproterozoic Fiery Creek Fault System, Mount Isa terrane, Australia;Peter G. Betts et al;《Journal of Structural Geology》;20011231;第23卷(第12期);第1953-1969页 |
反演断层数学模型的建立和应用;施泽进等;《成都理工学院学报》;19940131;第21卷(第1期);第91-98页 |
基于向-位错模型的同震断层运动反演——以龙门山断层为例;张永志等;《地球学报》;20160531;第37卷(第3期);第294-300页 |
基于地壳分层的地震断层同震变形反演分析;李志才等;《武汉大学学报信息科学版》;20080331;第33卷(第3期);第229-232页 |
断层几何形态反演的计算机模拟;陈伟等;《物探化探计算技术》;20081130;第30卷(第6期);第490-493页 |
Also Published As
Publication number | Publication date |
---|---|
CN106384016A (en) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
USRE49507E1 (en) | Faulted geological structures having unconformities | |
Natali et al. | Modeling Terrains and Subsurface Geology. | |
Wu et al. | An approach to computer modeling and visualization of geological faults in 3D | |
US11209560B2 (en) | Assignment of systems tracts | |
CN112381937B (en) | Multi-source geological data coupling modeling method based on drilling and complex geological profile | |
CN103646423B (en) | A kind of method of three-dimensional geological modeling and device | |
Guo et al. | Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm | |
EP3298437B1 (en) | Inversion for tectonic stress | |
NO335854B1 (en) | Automated system for modeling multiple value horizons with faults | |
CN104504754B (en) | A kind of method and device of oil and gas reservoir multi-point statistic modeling | |
CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
Omodeo Salè et al. | Tectono-stratigraphic evolution of an inverted extensional basin: the Cameros Basin (north of Spain) | |
Cosso et al. | Surveying and mapping a cave using 3d laser scanner: the open challenge with free and open source software | |
Wu et al. | On three-dimensional geological modeling and visualization | |
CN108986213B (en) | Three-dimensional stratum modeling method based on stacking technology | |
CN106384016B (en) | The automatic inversion method of single dip fault | |
Sala et al. | The Alpstein in three dimensions: fold-and-thrust belt visualization in the Helvetic zone, eastern Switzerland | |
CN107492130B (en) | Non-continuous equivalent cloud picture generation method based on layer rendering | |
Richards-Rissetto | From mounds to maps to models: visualizing ancient architecture across landscapes | |
CN113887046B (en) | Coal mine tunnel modeling method based on three-dimensional geologic body | |
CN111951394A (en) | Fault structure unit three-dimensional model construction method and device based on geological map | |
Lankford-Bravo | Structural and Stratigraphic History of a Deformed Diapir Margin, Onion Creek Salt Diapir, Paradox Basin, Utah | |
Poggi | Documentation and Analysis Workflow for the On-going Archaeological Excavation with Image-Based 3d Modelling Technique: the Case-study of the Medieval Site of Monteleo, Italy | |
Wang et al. | 3D integral modeling for city surface & subsurface | |
van Dijk | Application of Artificial Intelligence in 3D Modelling of Fault and Fracture Systems: Closing in on the Complexity of Geologic Reality |
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 |