CN113393578A - Construction method and device of geological profile three-dimensional model - Google Patents
Construction method and device of geological profile three-dimensional model Download PDFInfo
- Publication number
- CN113393578A CN113393578A CN202110612450.6A CN202110612450A CN113393578A CN 113393578 A CN113393578 A CN 113393578A CN 202110612450 A CN202110612450 A CN 202110612450A CN 113393578 A CN113393578 A CN 113393578A
- Authority
- CN
- China
- Prior art keywords
- section
- node
- line
- dimensional
- stratigraphic
- 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.)
- Pending
Links
- 238000010276 construction Methods 0.000 title claims description 8
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 239000007787 solid Substances 0.000 claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 17
- 230000001131 transforming effect Effects 0.000 claims description 12
- 238000013519 translation Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 description 4
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
Images
Classifications
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a method and a device for constructing a geological profile three-dimensional model, wherein the method comprises the following steps: (1) loading two-dimensional data to form a profile set SP, a stratum profile set S, a section line set SL and a grid set GEoDEM; (2) reading any profile SP from SPoConstructing a ground surface line GL and a reference line RL; (3 reading section spoPerforming space geometric transformation on the stratum section and the surface line GL to obtain a transformed stratum section and a transformed surface line; (4) converting the transformed stratigraphic section into three dimensions to obtain a three-dimensional stratigraphic section; (5) establishing a parallel stratum section according to a preset strip width w based on the three-dimensional stratum section, and constructing a stratum three-dimensional model by adopting the parallel stratum section and the preset strip width w; (6) and circularly executing the steps until all stratums of all the sections are traversed to obtain three-dimensional solid models of all the stratums, and embedding to obtain the three-dimensional geological section models. The invention realizes a method for quickly constructing a three-dimensional section model based on a two-dimensional geological section.
Description
Technical Field
The invention relates to the field of geology, in particular to a method and a device for constructing a geological profile three-dimensional model.
Background
The geological profile is a map reflecting the conditions of the ground layer and the geological structure of the earth surface and the underground at a certain depth in a certain direction, and is basic data for researching the ground layer, the rock mass and the structure. In addition, the geological profile map is matched with the plane geological map, so that the three-dimensional concept of the geological structure can be obtained, and the spatial correlation and the geological evolution relationship of the geological body and the geological structure can be reflected. Based on the important value of geological profiles in research and application, geoscientists have long produced and accumulated a large number of geological profiles.
However, with the continuous maturity of three-dimensional GIS technology and the increasing expansion of three-dimensional geoscience application systems, it is difficult for a traditional two-dimensional geological profile to meet the scientific research and application requirements of users in many aspects such as three-dimensional visual expression, three-dimensional spatial analysis, and the like. How to automatically generate a three-dimensional model of a local area of a geological profile based on a traditional two-dimensional geological profile to support three-dimensional visual expression and three-dimensional space analysis has important research significance and application value.
Disclosure of Invention
The purpose of the invention is as follows: the invention provides a method and a device for constructing a geological section three-dimensional model according to a two-dimensional geological sectional diagram, aiming at the problems in the prior art.
The technical scheme is as follows: the construction method of the geological profile three-dimensional model comprises the following steps:
(1) respectively cutting the section vector data, the section line vector data and the DEM according to the two-dimensional graph to form a section set SP, a stratum section set S, a section line set SL and a grid set GEoDEM;
(2) reading any profile SP from the set of profiles SPoBased on the stratum profile set S to the profile spoConstructing a ground surface line GL and a reference line RL;
(3) reading profiles sp from a set S of stratigraphic profilesoAny stratigraphic section s ofoiSectioning the formation s based on the section line set SL, the surface line GL and the reference line RLoiObtaining a transformed stratigraphic section s 'by performing a spatial geometric transformation with the surface line GL'oiAnd a change surface line GL';
(4) transforming the stratigraphic section s 'based on the grid set GEoDEM, the surface line GL and the transformed surface line GL'oiConverting to three dimensions to obtain a three-dimensional stratum section s'3oi；
(5) Based on three-dimensional stratum section s'3oiCreating a parallel stratigraphic section t from a preset strip width woiConstructing a three-dimensional model of the stratum by adopting the two models;
(6) circularly executing the steps (3) to (5) until the section sp is finishedoTraversing all stratigraphic sections;
(7) and (5) circularly executing the steps (2) to (6) until all the sections are traversed to obtain all the stratum three-dimensional solid models, and inlaying to obtain the three-dimensional geological section models.
Further, the step (1) comprises the following steps:
(1-1) reading two-dimensional geological profile vector data and storing the two-dimensional geological profile vector data into a profile set SP ═ { SP ═o1,2, …, SN, where sp isoDenotes a geological section numbered o, SN denotes the number of geological sections, and stores the strata in each section into a set of stratigraphic sections S ═ S oi1,2, …, SN, i 1,2, …, DNo }; wherein s isoiAn i-th stratigraphic section representing a geological section numbered o, DNo representing the number of stratigraphic sections of the geological section numbered o;
(1-2) reading the section line vector data and storing the section line set SL ═ { SL ═o1,2, …, SN, where sl isoRepresents spoThe cross-hatching of (a), SN represents the number of cross-hatching;
and (1-3) reading DEM data and storing the DEM data into the grid set GEoDEM.
Further, the step (2) comprises the following steps:
(2-1) reading any profile SP from the set of profiles SPoExtracting all stratum sections from the stratum section set S, merging the stratum sections to form an integral geometric surface, and extracting a boundary line BL of the geometric surface;
(2-2) calculating two nodes p on the boundary line BLlAnd prThe specific method comprises the following steps: sorting all nodes on the boundary line BL according to ascending order of the abscissa, obtaining a first node and a node with the same abscissa as the first node to construct a left node set LP, and obtaining a last node and a node with the same abscissa as the last node to construct a right node set RP; then sorting the node sets LP and RP according to the descending order of the ordinate, respectively obtaining the first node in the sets as the node plAnd pr；
(2-3) at node pl、prDividing the boundary BL to obtain two boundary segmentsThe middle origin is plEnd point is prThe boundary segments of (a) are taken as surface lines GL;
(2-4) with prThe abscissa is the abscissa, in plCreation of New node p 'for ordinate'rAnd with plIs a starting point p'rA reference line RL is constructed for the endpoint.
Further, the step (3) comprises the following steps:
(3-1) obtaining reference line RL starting points p respectivelylTerminal p'rAzimuth angle alpha and length LRLAnd spoSection line sl ofoStarting point pfEnd point ptAzimuth angle beta and length Lslo；
(3-2) reading a profile sp from the set S of stratigraphic profilesoAny stratigraphic section s ofoi；
(3-3) sectioning the formation s based on the start point, end point, azimuth angle and length obtained in the step (3-1)oiAffine transformation is carried out on all nodes on the boundary, and a transformed stratum section s 'is generated by adopting the transformed nodes'oi；
(3-4) obtaining a transformed ground line GL' by transforming the ground line GL according to the transformation method of the step (3-3).
Further, the step (3-3) comprises:
(3-3-1) cutting section s from the formationoiObtaining any node p on the boundarysi；
(3-3-2) starting the reference line RL by plTranslating to the origin, and calculating the translation quantity-x (p) in the x directionl) Y-direction translation-y (p)l)；
(3-3-3) calculating a rotation angle θ ═ β - α from the azimuth angle;
(3-3-4) calculating the ratio σ ═ L according to the lengthslo/LRL；
(3-3-5) starting the reference line RL by plShift to pfCalculating the x-direction translation x (p)f) Y-direction translation amount y (p)f)；
(3-3-6) node p is performed based on the following formulasiObtaining a transformed node p'siCoordinate (x (p'si)，y(p′si))：
(3-3-7) returning to execute the step (3-3-1) until all nodes on the boundary are traversed to obtain the transformation nodes of all boundaries;
(3-3-8) generating a transformed stratigraphic section s 'from the transformed nodes'oi。
Further, the step (4) comprises the following steps:
(4-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCalculating the elevation value z (p'si)；
(4-2) circularly executing the step (4-1) to finish transforming the stratum section s'oiThe elevation values of all upper boundary nodes;
(4-3) from the transformed formation section s'oiCreating a three-dimensional stratigraphic section s 'on the three-dimensional coordinate values of all boundary nodes'3oi。
Further, the step (4-1) comprises:
(4-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'si；
(4-1-2) get node p'siTo transformed formation section s'iUpper vertical point pSLAnd connect node p'siPerpendicular point pSLConstructing a line segment PL, extending two ends to obtain an intersection point p of the line segment PL and the ground surface line GLGL；
(4-1-3) respectively acquiring vertical points p based on the grid set GEoDEMSLAn elevation value Z of1And spoSection line sl ofoStarting point pfAn elevation value Z of2；
(4-1-4) if node p'siOn transformed ground surface line GL, then z (p'si)＝Z1Otherwise, executing the step (4-1-5);
(4-1-5) compute node p 'separately'siAnd the vertical point pSLDistance d of1Point of intersection pGLWith the plumbPoint pSLDistance d of2Of node p'siAnd point of intersection pGLDistance d of3And z (p 'is calculated based on the following formula'si)；
Further, the step (5) comprises the following steps:
(5-1) based on preset strip width w and three-dimensional stratum section s'3oiCreating a parallel stratigraphic section toi；
(5-2) transforming the stratigraphic section s 'based on a contour line algorithm'oiParallel stratigraphic section toiBoundary point creation side stoi；
(5-3) mosaic transformation of stratigraphic section s'oiParallel stratigraphic section toiAnd side stoiAnd obtaining a three-dimensional solid model of the stratum.
Further, the step (5-1) comprises:
(5-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCoordinate (x (p'si)，y(p′si)，z(p′si))；
(5-1-2) calculating the parallel section node p by the following equationtiCoordinate (x (p)ti),y(pti)，z(pti))；
Wherein β represents spoSection line sl ofoThe azimuth of (d);
(5-1-3) Return to loop through step (5-1-1) until the formation section s 'is transformed'oiTraversing all nodes on the boundary to obtain all parallel section nodes, and creating a parallel stratigraphic section t based on all the parallel section nodesoi。
The device for constructing the three-dimensional model of the geological profile 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 invention provides a construction method of a geological profile three-dimensional model through links such as data loading, construction under earth surface lines and references, geometric transformation of a stratum profile space, two-dimensional to three-dimensional conversion of a stratum profile, construction of a parallel stratum profile, construction of a stratum three-dimensional model and the like, can complete conversion from a two-dimensional geological sectional view to the three-dimensional geological model, and has high automation degree.
Drawings
FIG. 1 is a geological map (a is a grid map and b is a vector map) used in the present embodiment;
FIG. 2 is a plot of cross-hatching and DEM data as used in the present example;
FIG. 3 is a flow chart of a method for constructing a three-dimensional model of a geological profile according to the present invention;
FIG. 4 is a schematic representation of a surface line and a reference line in accordance with an embodiment of the present invention;
FIG. 5 is a cross-sectional view of an embodiment of the present invention after spatial geometry transformation;
FIG. 6 is a cross-sectional view of the three-dimensionality of an embodiment of the present invention (a is a single schematic view, and b is a cross-sectional view in combination with a DEM);
FIG. 7 is a formation s according to an embodiment of the invention1A front view and a side view;
FIG. 8 is a schematic diagram of a three-dimensional cross-sectional model of an embodiment of the invention.
Detailed Description
To explain the technical solution of the present invention in further detail, the experimental data of this embodiment uses a geological section (fig. 1) and a section line and DEM data (fig. 2) of a certain research area, and the projection coordinate system used by the section line and DEM data is WGS _1984_ UTM _ Zone _ 50N. The following further description is provided by describing a specific embodiment in conjunction with the accompanying drawings.
As shown in fig. 3, the present invention provides a method for constructing a geological section three-dimensional model, comprising:
(1) and respectively cutting the section vector data, the section line vector data and the DEM according to the two-dimensional graph to form a section set SP, a stratum section set S, a section line set SL and a grid set GEoDEM.
The method specifically comprises the following steps:
(1-1) reading two-dimensional geological profile vector data and storing the two-dimensional geological profile vector data into a profile set SP ═ { SP ═o1,2, …, SN, where sp isoDenotes a geological section numbered o, SN denotes the number of geological sections, and stores the strata in each section into a set of stratigraphic sections S ═ S oi1,2, …, SN, i 1,2, …, DNo }; wherein s isoiAn i-th stratigraphic section representing a geological section numbered o, DNo representing the number of stratigraphic sections of the geological section numbered o;
(1-2) reading the section line vector data and storing the section line set SL ═ { SL ═o1,2, …, SN, where sl isoRepresents spoThe cross-hatching of (a), SN represents the number of cross-hatching;
and (1-3) reading DEM data and storing the DEM data into the grid set GEoDEM.
(2) Reading any profile SP from the set of profiles SPoBased on the stratum profile set S to the profile spoA surface line GL and a reference line RL are constructed. In this example, the constructed surface line and the reference line are shown in fig. 4.
The method specifically comprises the following steps:
(2-1) reading any profile SP from the set of profiles SPoExtracting all stratum sections from the stratum section set S, merging the stratum sections to form an integral geometric surface, and extracting a boundary line BL of the geometric surface;
(2-2) calculating two nodes p on the boundary line BLlAnd prThe specific method comprises the following steps: sorting all nodes on the boundary line BL according to ascending order of the abscissa, obtaining a first node and a node with the same abscissa as the first node to construct a left node set LP, and obtaining a last node and a node with the same abscissa as the last node to construct a right node set RP; then sorting the node sets LP and RP according to the descending order of the ordinate, respectively obtaining the first node in the sets as the node plAnd pr；
(2-3) at node pl、prDividing the boundary BL to obtain two boundary segments with a starting point plEnd point is prThe boundary segments of (a) are taken as surface lines GL;
(2-4) with prThe abscissa is the abscissa, in plCreation of New node p 'for ordinate'rAnd with plIs a starting point p'rA reference line RL is constructed for the endpoint.
(3) Reading profiles sp from a set S of stratigraphic profilesoAny stratigraphic section s ofoiSectioning the formation s based on the section line set SL, the surface line GL and the reference line RLoiObtaining a transformed stratigraphic section s 'by performing a spatial geometric transformation with the surface line GL'oiAnd a change surface line GL'; . In the present embodiment, the converted geological profile is shown in fig. 5.
The method specifically comprises the following steps:
(3-1) obtaining reference line RL starting points p respectivelylTerminal p'rAzimuth angle alpha and length LRLAnd section line sloStarting point pfEnd point ptAzimuth angle beta and length Lslo(ii) a The azimuth calculation method is referred to the following documents: [1]Panzhengfeng, Yangzhenyao, Chengxiao, Chengshengjun, Wangtengjun, 2002, digital mapping principle and method, Wuhan university Press, P22.
(3-2) reading a profile sp from the set S of stratigraphic profilesoAny stratigraphic section s ofoi；
(3-3) sectioning the formation s based on the start point, end point, azimuth angle and length obtained in the step (3-1)oiAffine transformation is carried out on all nodes on the boundary, and a transformed stratum section s 'is generated by adopting the transformed nodes'oi；
(3-4) obtaining a transformed ground line GL' by transforming the ground line GL according to the transformation method of the step (3-3).
The step (3-3) comprises:
(3-3-1) cutting section s from the formationoiObtaining any node p on the boundarysi；
(3-3-2) starting the reference line RL by plTranslating to the origin, and calculating the translation quantity-x (p) in the x directionl) In the y directionTranslation-y (p)l)；
(3-3-3) calculating a rotation angle θ ═ β - α from the azimuth angle;
(3-3-4) calculating the ratio σ ═ L according to the lengthslo/LRL；
(3-3-5) starting the reference line RL by plShift to pfCalculating the x-direction translation x (p)f) Y-direction translation amount y (p)f)；
(3-3-6) node p is performed based on the following formulasiObtaining a transformed node p'siCoordinate (x (p'si)，y(p′si))：
(3-3-7) returning to execute the step (3-3-1) until all nodes on the boundary are traversed to obtain the transformation nodes of all boundaries;
(3-3-8) generating a transformed stratigraphic section s 'from the transformed nodes'oi。
(4) Transforming the stratigraphic section s 'based on the grid set GEoDEM, the surface line GL and the transformed surface line GL'oiConverting to three dimensions to obtain a three-dimensional stratum section s'3oi. In this example, a three-dimensional geological profile is constructed as shown in fig. 6.
The method specifically comprises the following steps:
(4-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCalculating the elevation value z (p'si)；
(4-2) circularly executing the step (4-1) to finish transforming the stratum section s'oiThe elevation values of all upper boundary nodes;
(4-3) from the transformed formation section s'oiCreating a three-dimensional stratigraphic section s 'on the three-dimensional coordinate values of all boundary nodes'3oi。
The step (4-1) comprises the following steps:
(4-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'si；
(4-1-2) get node p'siTo transformed formation section s'iUpper vertical point pSLAnd connect node p'siPerpendicular point pSLConstructing a line segment PL, extending two ends to obtain an intersection point p of the line segment PL and the ground surface line GLGL；
(4-1-3) respectively acquiring vertical points p based on the grid set GEoDEMSLAn elevation value Z of1And spoSection line sl ofoStarting point pfAn elevation value Z of2；
(4-1-4) if node p'siOn transformed ground surface line GL, then z (p'si)＝Z1Otherwise, executing the step (4-1-5);
(4-1-5) compute node p 'separately'siAnd the vertical point pSLDistance d of1Point of intersection pGLAnd the vertical point pSLDistance d of2Of node p'siAnd point of intersection pGLDistance d of3And z (p 'is calculated based on the following formula'si)；
(5) Based on three-dimensional stratum section s'3oiCreating a parallel stratigraphic section t from a preset strip width woiAnd constructing a three-dimensional model of the stratum by adopting the two models.
The method specifically comprises the following steps:
(5-1) based on preset strip width w and three-dimensional stratum section s'3oiCreating a parallel stratigraphic section toi；
(5-2) transforming the stratigraphic section s 'based on a contour line algorithm'oiParallel stratigraphic section toiBoundary point creation side sti(ii) a The contour line algorithm is referred to in the following documents: meyers D, Skinner S, Sloan K.1992.surfaces from contacts. ACM Transactions On Graphics,11 (3); 228-258.
(5-3) mosaic transformation of stratigraphic section s'oiParallel stratigraphic section toiKneading dough stoiTo obtain a three-dimensional solid of the stratumA body model. In this embodiment, when i is 1, the formation s1The three-dimensional model is shown in fig. 7.
The step (5-1) comprises the following steps:
(5-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCoordinate (x (p'si)，y(p′si)，z(p′si))；
(5-1-2) calculating the parallel section node p by the following equationtiCoordinate (x (p)ti)，y(pti)，z(pti))；
Wherein β represents spoSection line sl ofoThe azimuth of (d);
(5-1-3) Return to loop through step (5-1-1) until the formation section s 'is transformed'oiTraversing all nodes on the boundary to obtain all parallel section nodes, and creating a parallel stratigraphic section t based on all the parallel section nodesoi。。
(6) Circularly executing the steps (3) to (5) until the section sp is finishedoTraversing all stratigraphic sections;
(7) and (5) circularly executing the steps (2) to (6) until all the sections are traversed to obtain all the stratum three-dimensional solid models, and inlaying to obtain the three-dimensional geological section models.
In this embodiment, a three-dimensional model of a geological map is constructed as shown in FIG. 8. In the embodiment of the invention, partial GIS operation is provided based on the Arcgis Engine API, and related steps can also use the APIs of software such as SuperMap, Arcgis Object and the like to carry out corresponding GIS operation.
The invention also provides a device for constructing the geological profile three-dimensional model, which comprises a processor and a computer program stored on a memory and capable of running on the processor, wherein the processor realizes the method when executing the program.
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 (10)
1. A construction method of a geological profile three-dimensional model is characterized by comprising the following steps:
(1) respectively cutting the section vector data, the section line vector data and the DEM according to the two-dimensional graph to form a section set SP, a stratum section set S, a section line set SL and a grid set GEoDEM;
(2) reading any profile SP from the set of profiles SPoBased on the stratum profile set S to the profile spoConstructing a ground surface line GL and a reference line RL;
(3) reading profiles sp from a set S of stratigraphic profilesoAny stratigraphic section s ofoiSectioning the formation s based on the section line set SL, the surface line GL and the reference line RLoiObtaining a transformed stratigraphic section s 'by performing a spatial geometric transformation with the surface line GL'oiAnd a change surface line GL';
(4) transforming the stratigraphic section s 'based on the grid set GEoDEM, the surface line GL and the transformed surface line GL'oiConverting to three dimensions to obtain a three-dimensional stratum section s'3oi；
(5) Based on three-dimensional stratum section s'3oiCreating a parallel stratigraphic section t from a preset strip width woiConstructing a three-dimensional model of the stratum by adopting the two models;
(6) circularly executing the steps (3) to (5) until the section sp is finishedoTraversing all stratigraphic sections;
(7) and (5) circularly executing the steps (2) to (6) until all the sections are traversed to obtain all the stratum three-dimensional solid models, and inlaying to obtain the three-dimensional geological section models.
2. The method of constructing a three-dimensional model of a geological profile according to claim 1, characterized in that: the step (1) comprises the following steps:
(1-1) reading two-dimensional geological profile vector data and storing the two-dimensional geological profile vector data into a profile set SP ═ { SP ═o1,2, …, SN, where sp isoDenotes the geological section numbered o, SN denotes the number of geological sectionsAnd storing the stratum in each section into a stratum section set S ═ Soi1,2, …, SN, i 1,2, …, DNo }; wherein s isoiAn i-th stratigraphic section representing a geological section numbered o, DNo representing the number of stratigraphic sections of the geological section numbered o;
(1-2) reading the section line vector data and storing the section line set SL ═ { SL ═o1,2, …, SN, where sl isoRepresents spoThe cross-hatching of (a), SN represents the number of cross-hatching;
and (1-3) reading DEM data and storing the DEM data into the grid set GEoDEM.
3. The method of constructing a three-dimensional model of a geological profile according to claim 1, characterized in that: the step (2) comprises the following steps:
(2-1) reading any profile SP from the set of profiles SPoExtracting all stratum sections from the stratum section set S, merging the stratum sections to form an integral geometric surface, and extracting a boundary line BL of the geometric surface;
(2-2) calculating two nodes p on the boundary line BLlAnd prThe specific method comprises the following steps: sorting all nodes on the boundary line BL according to ascending order of the abscissa, obtaining a first node and a node with the same abscissa as the first node to construct a left node set LP, and obtaining a last node and a node with the same abscissa as the last node to construct a right node set RP; then sorting the node sets LP and RP according to the descending order of the ordinate, respectively obtaining the first node in the sets as the node plAnd pr；
(2-3) at node pl、prDividing the boundary BL to obtain two boundary segments with a starting point plEnd point is prThe boundary segments of (a) are taken as surface lines GL;
(2-4) with prThe abscissa is the abscissa, in plCreation of New node p 'for ordinate'rAnd with plIs a starting point p'rA reference line RL is constructed for the endpoint.
4. The method of constructing a three-dimensional model of a geological profile according to claim 1, characterized in that: the step (3) comprises the following steps:
(3-1) obtaining reference line RL starting points p respectivelylTerminal p'rAzimuth angle alpha and length LRLAnd spoSection line sl ofoStarting point pfEnd point ptAzimuth angle beta and length Lslo；
(3-2) reading a profile sp from the set S of stratigraphic profilesoAny stratigraphic section s ofoi；
(3-3) sectioning the formation s based on the start point, end point, azimuth angle and length obtained in the step (3-1)oiAffine transformation is carried out on all nodes on the boundary, and a transformed stratum section s 'is generated by adopting the transformed nodes'oi；
(3-4) obtaining a transformed ground line GL' by transforming the ground line GL according to the transformation method of the step (3-3).
5. The method of constructing a three-dimensional model of a geological profile according to claim 4, wherein: the step (3-3) comprises:
(3-3-1) cutting section s from the formationoiObtaining any node p on the boundarysi；
(3-3-2) starting the reference line RL by plTranslating to the origin, and calculating the translation quantity-x (p) in the x directionl) Y-direction translation-y (p)l)；
(3-3-3) calculating a rotation angle θ ═ β - α from the azimuth angle;
(3-3-4) calculating the ratio σ ═ L according to the lengthslo/LRL；
(3-3-5) starting the reference line RL by plShift to pfCalculating the x-direction translation x (p)f) Y-direction translation amount y (p)f)；
(3-3-6) node p is performed based on the following formulasiObtaining a transformed node p'siCoordinate (x (p'si)，y(p′si))：
(3-3-7) returning to execute the step (3-3-1) until all nodes on the boundary are traversed to obtain the transformation nodes of all boundaries;
(3-3-8) generating a transformed stratigraphic section s 'from the transformed nodes'oi。
6. The method of constructing a three-dimensional model of a geological profile according to claim 1, characterized in that: the step (4) comprises the following steps:
(4-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCalculating the elevation value z (p'si)；
(4-2) circularly executing the step (4-1) to finish transforming the stratum section s'oiThe elevation values of all upper boundary nodes;
(4-3) from the transformed formation section s'oiCreating a three-dimensional stratigraphic section s 'on the three-dimensional coordinate values of all boundary nodes'3oi。
7. The method of constructing a three-dimensional model of a geological profile according to claim 6, wherein: the step (4-1) comprises the following steps:
(4-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'si；
(4-1-2) get node p'siTo transformed formation section s'iUpper vertical point pSLAnd connect node p'siPerpendicular point pSLConstructing a line segment PL, extending two ends to obtain an intersection point p of the line segment PL and the ground surface line GLGL；
(4-1-3) respectively acquiring vertical points p based on the grid set GEoDEMSLAn elevation value Z of1And spoSection line sl ofoStarting point pfAn elevation value Z of2；
(4-1-4) if node p'siOn transformed ground surface line GL, then z (p'si)＝Z1Otherwise, executing the step (4-1-5);
(4-1-5) compute node p 'separately'siAnd the vertical point pSLDistance d of1Point of intersection pGLWith the plumbPoint pSLDistance d of2Of node p'siAnd point of intersection pGLDistance d of3And z (p 'is calculated based on the following formula'si)；
8. The method of constructing a three-dimensional model of a geological profile according to claim 1, characterized in that: the step (5) comprises the following steps:
(5-1) based on preset strip width w and three-dimensional stratum section s'3oiCreating a parallel stratigraphic section toi；
(5-2) transforming the stratigraphic section s 'based on a contour line algorithm'oiParallel stratigraphic section toiBoundary point creation side stoi；
(5-3) mosaic transformation of stratigraphic section s'oiParallel stratigraphic section toiAnd side stoiAnd obtaining a three-dimensional solid model of the stratum.
9. The method of constructing a three-dimensional model of a geological profile according to claim 8, wherein: the step (5-1) comprises the following steps:
(5-1-1) reading converted stratigraphic section s'oiAny node p 'on boundary'siCoordinate (x (p'si)，y(p′si)，z(p′si))；
(5-1-2) calculating the parallel section node p by the following equationtiCoordinate (x (p)ti)，y(pti)，z(pti))；
Wherein β represents spoSection line sl ofoThe azimuth of (d);
(5-1-3) Return to loop through step (5-1-1) until the formation section s 'is transformed'oiAt the boundaryThe nodes are traversed to obtain all parallel profile nodes, and a parallel stratigraphic profile t is created based on all parallel profile nodesoi。
10. An apparatus for constructing a three-dimensional model of a geological profile, 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-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110612450.6A CN113393578A (en) | 2021-06-02 | 2021-06-02 | Construction method and device of geological profile three-dimensional model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110612450.6A CN113393578A (en) | 2021-06-02 | 2021-06-02 | Construction method and device of geological profile three-dimensional model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113393578A true CN113393578A (en) | 2021-09-14 |
Family
ID=77619883
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110612450.6A Pending CN113393578A (en) | 2021-06-02 | 2021-06-02 | Construction method and device of geological profile three-dimensional model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113393578A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115984503A (en) * | 2023-03-20 | 2023-04-18 | 红石恒信(成都)科技有限公司 | Geological profile generation method, geological profile generation system, electronic equipment and medium |
-
2021
- 2021-06-02 CN CN202110612450.6A patent/CN113393578A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115984503A (en) * | 2023-03-20 | 2023-04-18 | 红石恒信(成都)科技有限公司 | Geological profile generation method, geological profile generation system, electronic equipment and medium |
| CN115984503B (en) * | 2023-03-20 | 2023-06-13 | 红石恒信(成都)科技有限公司 | Geological profile generation method, system, electronic equipment and medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106780751A (en) | Three-dimensional point cloud method for reconstructing based on improved shielding Poisson algorithm | |
| CN102831647B (en) | Delaunay triangulation network subdividing method based on space curved surface constraint | |
| CN101246597A (en) | Method for commercial scale vector map data gradual transmission and dynamic reconstruction | |
| CN104240300A (en) | Large-scale point cloud complex space curved face reconstruction method based on distributed parallelism | |
| Demir et al. | Procedural editing of 3d building point clouds | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN104167020A (en) | Spatial reversed curved surface reconstruction method based on restraint Delaunay triangular mesh | |
| CN111415415A (en) | Automatic construction method of three-dimensional map cutting geological profile | |
| CN113393578A (en) | Construction method and device of geological profile three-dimensional model | |
| CN113593038A (en) | Tunnel point cloud central line automatic extraction and triangulation network construction method | |
| Zhu et al. | Propagation strategies for stereo image matching based on the dynamic triangle constraint | |
| CN112231423B (en) | Geological map-based fault structure automatic recovery method and device | |
| CN110058298B (en) | Three-dimensional geologic body spatial interpolation method and system | |
| CN105302951A (en) | Finite element mesh surface triangulation method | |
| WO2012094013A1 (en) | Systems and methods for the construction of closed bodies during 3d modeling | |
| Gonzaga de Oliveira | A review on delaunay refinement techniques | |
| Meng et al. | An IFS-based fractal discrete fracture network for hydraulic fracture behavior of rock mass | |
| Fang et al. | Volumetric data modeling and analysis based on seven-directional box spline | |
| CN111951394A (en) | Fault structure unit three-dimensional model construction method and device based on geological map | |
| CN112233230A (en) | Three-dimensional model construction method and device for fault structure in map-cut geological profile | |
| Hidaka et al. | Polygonization of point cloud of tunnels using lofting operation | |
| KR100397085B1 (en) | Method for the graph representation of 3-dimensional object image for image retrival | |
| Xu et al. | 3D GIS based long-distance project construction progress visualization | |
| Yuan et al. | Filtering and reduction for 3-dimensional surface modeling of laser line scanning point cloud |
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 |