CN112581558A - Model construction method and system for intrusion structure in map-cut geological profile - Google Patents
Model construction method and system for intrusion structure in map-cut geological profile Download PDFInfo
- Publication number
- CN112581558A CN112581558A CN202011154215.0A CN202011154215A CN112581558A CN 112581558 A CN112581558 A CN 112581558A CN 202011154215 A CN202011154215 A CN 202011154215A CN 112581558 A CN112581558 A CN 112581558A
- Authority
- CN
- China
- Prior art keywords
- stratum
- invaded
- points
- point
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000010276 construction Methods 0.000 title claims abstract description 6
- 239000011435 rock Substances 0.000 claims abstract description 56
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000004590 computer program Methods 0.000 claims description 3
- 230000009545 invasion Effects 0.000 abstract description 8
- 238000005755 formation reaction Methods 0.000 description 33
- 241000282376 Panthera tigris Species 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
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)
Abstract
The invention discloses a model construction method and a device for an intrusion structure in a map-cutting geological profile, wherein the method comprises the following steps: (1) loading a stratum layer and a stratum boundary layer of the map-cut geological profile; (2) according to the stratum code, extracting invaded stratum data; (3) extracting any invaded rock stratum earth surface intersection point; (4) calculating a bottom intersection point and a control point of the invaded rock stratum; (5) generating a smooth invaded rock stratum inferred formation curve according to the earth surface intersection point, the bottom intersection point and the control point; (6) according to the inferred stratigraphic curve, dividing the original stratigraphic surface elements; (7) and combining to construct the invaded rock stratum. The method can effectively improve the modeling quality of the invasion structure in the map-cut geological profile.
Description
Technical Field
The invention relates to the field of geographic information and geology, in particular to a method and a system for building an intrusion structure model in a map-cut geological profile.
Background
Invasion configuration refers to the configuration created during the process of gradually condensing and consolidating into rock as the temperature decreases during the flow period as part of the crust material invades the overlying surrounding rock in a vertical direction. The invasion structure is not only indispensable content in the research of regional structure, but also has important practical significance in the exploration work of prospecting. The map-cutting geological profile is a geological profile which is compiled by selecting a certain direction on a geological map and using a projection method according to various geological and geographic elements and a certain scale. On a map-cut geological profile, reasonable inference and modeling of an invaded rock mass are carried out based on the stratum attitude and the plane geometric form of an invaded rock layer, which is a necessary premise for accurately modeling the map-cut geological profile in an invaded rock development area and is a requirement for three-dimensional modeling based on a sequence map-cut geological profile. Therefore, the research of the modeling method of the invasion structure in the map-cutting geological profile has important research significance and use value.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems in the prior art, the invention provides a method and a system for constructing a model of an invasion structure in a map-cut geological section with higher modeling quality.
The technical scheme is as follows: the method for constructing the model of the invasion structure in the map-cutting geological profile comprises the following steps:
(1) loading a stratum layer and a stratum line map layer of the map-cutting geological profile to obtain a stratum code set FD and a stratum line set LN of all the stratums;
(2) extracting a set of invaded rock layers FN according to the stratum codes;
(3) reading any invaded rock layer FN from the invaded rock layer set FN to obtain a surface intersection point A, B;
(4) calculating bottom intersection points C 'and D' of the invaded rock layer fn, stratum line core control points H 'and I' and curve control points according to the earth surface intersection point A, B;
(5) generating a pair of smooth invaded rock stratum inferred formation curves according to the earth surface intersection point, the bottom intersection point and the control point;
(6) circularly executing the steps (3) - (5) until the traversal of the invaded rock stratum set FN is completed, obtaining all invaded rock stratum inferred stratum curves, and storing the invaded rock stratum inferred stratum curves into the invaded rock stratum inferred stratum curve set GL;
(7) and deducing a stratum curve set GL based on the invaded stratum, segmenting the original stratum surface, combining and constructing the invaded stratum, and obtaining a model of the invaded structure.
Further, the step (1) specifically comprises:
(1-1) loading a stratum layer and a stratum line layer of the map-cut geological profile, acquiring codes of all stratums from the stratum layer and storing the codes into the groundLayer code set FD ═ fc1,2, …, fc }, where fcRepresenting the c-th stratum code, fc representing the number of stratum codes;
(1-2) all the stratum lines are obtained and stored in the stratum line set LN ═ { LN ═ LN) in sequence along the section line directionb1,2, …, fc-1, where lnbRepresents the b-th stratigraphic line and fc-1 represents the number of stratigraphic lines.
Further, the step (2) specifically comprises:
(2-1) constructing an invaded formation code table FC according to the legend of the geological map of the research area;
(2-2) sorting the stratigraphic codes in the stratigraphic code set FD from left to right along the section line direction;
and (2-3) reading the elements from the stratum code set FD in sequence, and if the elements belong to the stratum codes in the invaded stratum code table FC, storing the elements into the invaded stratum set FN until the reading is finished, thus obtaining the invaded stratum set FN.
Further, the step (3) specifically comprises:
(3-1) reading any invaded formation FN based on the invaded formation set FN;
(3-2) acquiring two stratum lines LN corresponding to the invaded stratum fn according to the stratum line set LNb、lnb+1；
(3-3) separately extracting the horizon lnb、lnb+1Head end point A (x)A,yA)、B(xB,yB) I.e. the surface intersection of the invaded formation fn.
Further, the step (4) specifically comprises:
(4-1) connecting the points A, B, and calculating the difference h between the abscissa and the abscissa of the AB;
(4-2) obtaining the Point C (x) according to the following formulaC,yC) Point D (x)D,yD) And midpoint E (x) of ABE,yE)；
(4-3) respectively making a vertical line through the point C, D, E, wherein the vertical line is respectively intersected with the bottom boundary line of the stratum at the points C ', D' and F, and the points C 'and D' are the bottom intersection points of the invaded stratum fn;
(4-4) obtaining a middle point G of EF, drawing a perpendicular line of EF passing through the point G, and respectively crossing AC 'at a point H and crossing BD' at a point I;
(4-5) obtaining a midpoint H 'of the GH and a midpoint I' of the GI, wherein the points H 'and I' are stratum line core control points of the invaded rock stratum;
and (4-6) calculating a curve control point of the invaded rock stratum according to the inclined direction of the stratigraphic line of the invaded rock stratum.
Further, the step (4-6) specifically comprises:
(4-6-1) connecting CC ', DD', CC 'and DD' with HI at a point S, L respectively;
(4-6-2) calculating midpoint P of point A, S and midpoint Q of point B, L;
(4-6-3) acquiring two formation lines ln corresponding to the invaded formation fnb、lnb+1Calculating the horizon lnbSlope k of1And the formation line lnb+1Slope k of2；
(4-6-4) if k1>0 and k2>0, taking the midpoint P as a curve control point; if k is1<0 and k2<0, taking the middle point Q as a curve control point; if k is1>0 and k2<0, then the midpoint P, Q is simultaneously taken as the curve control point; if k is1<0 and k2>0, there is no curve control point.
Further, the step (5) specifically comprises:
(5-1) if two curve control points P, Q exist, respectively generating a smooth curve between the points A, P, H ', C' and B, Q, I ', D';
(5-2) if only one curve control point P exists, respectively generating a smooth curve between the points A, P, H ', C' and B, I ', D'; if only one curve control point Q exists, a smooth curve is respectively generated between the points A, H ', C' and the points B, Q, I 'and D';
(5-3) if there is no curve control point, a smooth curve is generated between the points A, H ', C' and B, I ', D', respectively.
Further, the method for generating the smooth curve between 4 points in the step (5-1) and the step (5-2) specifically comprises the following steps:
A. sequentially connecting 4 points, and taking the middle two points as vertexes;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 4 points;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
Further, the method for generating the smooth curve between 3 points in the step (5-2) and the step (5-3) specifically comprises the following steps:
A. sequentially connecting 3 points, and taking the middle point as a vertex;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 3 points;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
Further, the step (7) specifically comprises:
(7-1) inferring from the invaded rock a set of formation curves GLReading a pair of invaded rock strata inferred formation curves glpAnd glp+1；
(7-2) inferring a formation curve gl from the invaded rockpAnd glp+1Dividing the original ground level;
(7-3) obtaining the horizon gl based on the graph space relationpAnd glp+1The split horizon in between;
(7-4) combining all the divided stratum surfaces to construct an invaded rock stratum which is used as an invaded structure stratum obtained by modeling;
(7-5) performing steps (7-1) - (7-4) in a loop until the construction of all invasive configurations is completed.
The device for constructing the model of the invasive structure in the map-cutting geological section 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 method can effectively improve the modeling quality of the interrupted invasion structure of the map-cut geological profile, and lays a good data foundation for three-dimensional modeling based on the sequence map-cut geological profile.
Drawings
FIG. 1 is a plot of geological profile data as used in the present embodiment;
FIG. 2 is a flow chart of a method of modeling an intrusive feature in a sliced geological profile provided by the present invention;
FIG. 3 is a schematic diagram of calculating a bottom intersection point and a control point of an invaded formation;
FIG. 4 is a cross-sectional view after generating an inferred formation profile of the invaded formation;
FIG. 5 is a cross-sectional view of an intrusion structure after modeling in the present invention.
Detailed Description
As will be described in further detail below, in this embodiment, a map-cut geological section of a tiger hole in south kyo is selected as experimental data, and as shown in fig. 1, a projection coordinate system adopted by the experimental data is WGS 84. The following further description is provided by describing a specific embodiment in conjunction with the accompanying drawings.
As shown in fig. 2, the method for constructing a model of an intrusive structure in a cutback geological profile provided by the present embodiment includes the following steps:
(1) and loading the stratum layer and the stratum line map layer of the map-cutting geological profile to obtain a stratum code set FD and a stratum line set LN of all the stratums.
The method specifically comprises the following steps:
(1-1) loading a stratum layer and a stratum line layer of the map-cut geological profile, acquiring codes of all the stratums from the stratum layer and storing a stratum code set FD ═ fc1,2, …, fc }, where fcRepresenting the c-th stratum code, fc representing the number of stratum codes; in this embodiment, fc is 15.
(1-2) all the stratum lines are obtained and stored in the stratum line set LN ═ { LN ═ LN) in sequence along the section line directionb1,2, …, fc-1, where lnbRepresents the b-th stratigraphic line and fc-1 represents the number of stratigraphic lines.
(2) According to the stratigraphic code, a set of invaded rock layers FN is extracted.
The method specifically comprises the following steps:
(2-1) constructing an invaded rock stratum code table FC according to the legend of the geological map of the research region, namely constructing the stratum code tables FC of all invaded rock strata according to stratum code coding rules of the legend of the geological map of the research region so as to enable the stratum code types to be consistent; the invaded formation code table FC established in this embodiment is shown in table 1;
table 1 invasion formation code table FC
(2-2) sorting the stratigraphic codes in the stratigraphic code set FD from left to right along the section line direction;
(2-3) reading the elements from the stratum code set FD in sequence, if the elements belong to the stratum codes in the invaded stratum code table FC, storing the elements into the invaded stratum set FN until the reading is finished, and obtaining the invaded stratum set FN.
(3) Any invaded formation FN is read from the invaded formation set FN to obtain its surface intersection A, B.
The method specifically comprises the following steps:
(3-1) reading any invaded formation FN based on the invaded formation set FN;
(3-2) acquiring two stratum lines LN corresponding to the invaded stratum fn according to the stratum line set LNb、lnb+1；
(3-3) separately extracting the horizon lnb、lnb+1Head end point A (x)A,yA)、B(xB,yB) I.e. the surface intersection of the invaded formation fn.
(4) Bottom intersection points C ', D' of the invaded rock layer fn, formation line core control points H ', I' and curve control points are calculated from the surface intersection point A, B.
As shown in fig. 3, the steps specifically include:
(4-1) connecting the points A, B, and calculating the difference h between the abscissa and the abscissa of the AB;
(4-2) obtaining the Point C (x) according to the following formulaC,yC) Point D (x)D,yD) And midpoint E (x) of ABE,yE)；
(4-3) respectively making a vertical line through the point C, D, E, wherein the vertical line is respectively intersected with the bottom boundary line of the stratum at the points C ', D' and F, and the points C 'and D' are the bottom intersection points of the invaded stratum fn;
(4-4) obtaining a middle point G of EF, drawing a perpendicular line of EF passing through the point G, and respectively crossing AC 'at a point H and crossing BD' at a point I;
(4-5) obtaining a midpoint H 'of the GH and a midpoint I' of the GI, wherein the points H 'and I' are stratum line core control points of the invaded rock stratum;
and (4-6) calculating a curve control point of the invaded rock stratum according to the inclined direction of the stratigraphic line of the invaded rock stratum.
The steps (4-6) specifically include:
(4-6-1) connecting CC ', DD', CC 'and DD' with HI at a point S, L respectively;
(4-6-2) calculating midpoint P of point A, S and midpoint Q of point B, L;
(4-6-3) acquiring two formation lines ln corresponding to the invaded formation fnb、lnb+1Calculating the horizon lnbSlope k of1And the formation line lnb+1Slope k of2；
(4-6-4) if k1>0 and k2>0, taking the midpoint P as a curve control point; if k is1<0 and k2<0, taking the middle point Q as a curve control point; if k is1>0 and k2<0, then the midpoint P, Q is simultaneously taken as the curve control point; if k is1<0 and k2>0, there is no curve control point. In this embodiment, only one curve control point Q exists for both invaded formations.
(5) And generating a pair of smooth invaded rock stratum inferred formation curves according to the surface intersection point, the bottom intersection point and the control point.
The method specifically comprises the following steps:
(5-1) if two curve control points P, Q exist, respectively generating a smooth curve between the points A, P, H ', C' and B, Q, I ', D';
(5-2) if only one curve control point P exists, respectively generating a smooth curve between the points A, P, H ', C' and B, I ', D'; if only one curve control point Q exists, a smooth curve is respectively generated between the points A, H ', C' and the points B, Q, I 'and D';
(5-3) if there is no curve control point, a smooth curve is generated between the points A, H ', C' and B, I ', D', respectively.
The method for generating the smooth curve between 4 points in the step (5-1) and the step (5-2) specifically comprises the following steps:
A. sequentially connecting 4 points, and taking the middle two points as vertexes; for example, if a curve is generated between A, P, H ' and C ', A, P, H ' and C ' are connected in sequence, with P, H ' as a vertex;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 4 points; for example, assuming that four points are A, P, H ', C', the 1 st, 2 nd, 3 th, 4 th points are A, P, H ', C', respectively;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
The method for generating the smooth curve between 3 points in the step (5-2) and the step (5-3) specifically comprises the following steps:
A. sequentially connecting 3 points, and taking the middle point as a vertex; for example, if a curve is generated between B, I ' and D ', B, I ' and D ' are connected in order, with I ' as the vertex;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 3 points; assuming that 3 points are B, I 'and D', the 1 st, 2 nd and 3 rd points represent points B, I 'and D', respectively;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
(6) And (5) circularly executing the steps (3) to (5) until the traversal of the invaded rock stratum set FN is completed, obtaining all invaded rock stratum inferred formation curves, and storing the invaded rock stratum inferred formation curves GL.
(7) And deducing a stratum curve set GL based on the invaded stratum, segmenting the original stratum surface, combining and constructing the invaded stratum, and obtaining a model of the invaded structure.
The method specifically comprises the following steps:
(7-1) reading a pair of invaded rock layer inferred formation curves GL from the invaded rock layer inferred formation curve set GLpAnd glp+1；
(7-2) inferring a formation curve gl from the invaded rockpAnd glp+1Segmenting the original ground plane, as shown in fig. 4;
(7-3) obtaining the horizon gl based on the graph space relationpAnd glp+1The split horizon in between;
(7-4) combining all the divided stratum surfaces to construct an invaded rock stratum which is used as an invaded structure stratum obtained by modeling;
(7-5) performing steps (7-1) - (7-4) in a loop until the construction of all invasive configurations is completed. In this embodiment, a constructed cut-map geological profile of the invaded formation is shown in FIG. 5.
The embodiment also provides a model building device for cutting the invasive structure in the geological section, 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 method of modeling an intrusive feature in a cut-to-map geological profile, the method comprising:
(1) loading a stratum layer and a stratum line map layer of the map-cutting geological profile to obtain a stratum code set FD and a stratum line set LN of all the stratums;
(2) extracting a set of invaded rock layers FN according to the stratum codes;
(3) reading any invaded rock layer FN from the invaded rock layer set FN to obtain a surface intersection point A, B;
(4) calculating bottom intersection points C 'and D' of the invaded rock layer fn, stratum line core control points H 'and I' and curve control points according to the earth surface intersection point A, B;
(5) generating a pair of smooth invaded rock stratum inferred formation curves according to the earth surface intersection point, the bottom intersection point and the control point;
(6) circularly executing the steps (3) - (5) until the traversal of the invaded rock stratum set FN is completed, obtaining all invaded rock stratum inferred stratum curves, and storing the invaded rock stratum inferred stratum curves into the invaded rock stratum inferred stratum curve set GL;
(7) and deducing a stratum curve set GL based on the invaded stratum, segmenting the original stratum surface, combining and constructing the invaded stratum, and obtaining a model of the invaded structure.
2. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 1, wherein: the step (1) specifically comprises the following steps:
(1-1) loading a stratum layer and a stratum line layer of the map-cut geological profile, acquiring codes of all the stratums from the stratum layer and storing a stratum code set FD ═ fc1,2, …, fc }, where fcRepresenting the c-th stratum code, fc representing the number of stratum codes;
(1-2) all the stratum lines are obtained and stored in the stratum line set LN ═ { LN ═ LN) in sequence along the section line directionb1,2, …, fc-1, where lnbRepresents the b-th stratigraphic line and fc-1 represents the number of stratigraphic lines.
3. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 1, wherein: the step (2) specifically comprises the following steps:
(2-1) constructing an invaded formation code table FC according to the legend of the geological map of the research area;
(2-2) sorting the stratigraphic codes in the stratigraphic code set FD from left to right along the section line direction;
(2-3) reading elements from the stratum code set FD in sequence, if the elements belong to the stratum codes in the invaded stratum code table FC, storing the elements into the invaded stratum set FN until the elements are read completely, and obtaining the invaded stratum set FN;
the step (3) specifically comprises the following steps:
(3-1) reading any invaded formation FN based on the invaded formation set FN;
(3-2) acquiring two stratum lines LN corresponding to the invaded stratum fn according to the stratum line set LNb、lnb+1；
(3-3) separately extracting the horizon lnb、lnb+1Head end point A (x)A,yA)、B(xB,yB) I.e. the surface intersection of the invaded formation fn.
4. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 1, wherein: the step (4) specifically comprises the following steps:
(4-1) connecting the points A, B, and calculating the difference h between the abscissa and the abscissa of the AB;
(4-2) obtaining the Point C (x) according to the following formulaC,yC) Point D (x)D,yD) And midpoint E (x) of ABE,yE)；
(4-3) respectively making a vertical line through the point C, D, E, wherein the vertical line is respectively intersected with the bottom boundary line of the stratum at the points C ', D' and F, and the points C 'and D' are the bottom intersection points of the invaded stratum fn;
(4-4) obtaining a middle point G of EF, drawing a perpendicular line of EF passing through the point G, and respectively crossing AC 'at a point H and crossing BD' at a point I;
(4-5) obtaining a midpoint H 'of the GH and a midpoint I' of the GI, wherein the points H 'and I' are stratum line core control points of the invaded rock stratum;
and (4-6) calculating a curve control point of the invaded rock stratum according to the inclined direction of the stratigraphic line of the invaded rock stratum.
5. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 4, wherein: the steps (4-6) specifically include:
(4-6-1) connecting CC ', DD', CC 'and DD' with HI at a point S, L respectively;
(4-6-2) calculating midpoint P of point A, S and midpoint Q of point B, L;
(4-6-3) acquiring two formation lines ln corresponding to the invaded formation fnb、lnb+1Calculating the horizon lnbSlope k of1And the formation line lnb+1Slope k of2；
(4-6-4) if k1>0 and k2>0, taking the midpoint P as a curve control point; if k is1<0 and k2<0, taking the middle point Q as a curve control point; if k is1>0 and k2<0, then the midpoint P, Q is simultaneously taken as the curve control point; if k is1<0 and k2>0, there is no curve control point.
6. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 1, wherein: the step (5) specifically comprises the following steps:
(5-1) if two curve control points P, Q exist, respectively generating a smooth curve between the points A, P, H ', C' and B, Q, I ', D';
(5-2) if only one curve control point P exists, respectively generating a smooth curve between the points A, P, H ', C' and B, I ', D'; if only one curve control point Q exists, a smooth curve is respectively generated between the points A, H ', C' and the points B, Q, I 'and D';
(5-3) if there is no curve control point, a smooth curve is generated between the points A, H ', C' and B, I ', D', respectively.
7. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 6, wherein: the method for generating the smooth curve between 4 points in the step (5-1) and the step (5-2) specifically comprises the following steps:
A. sequentially connecting 4 points, and taking the middle two points as vertexes;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 4 points;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
8. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 6, wherein: the method for generating the smooth curve between 3 points in the step (5-2) and the step (5-3) specifically comprises the following steps:
A. sequentially connecting 3 points, and taking the middle point as a vertex;
B. setting the initial value of the parameter t as 0, setting the t to be [0,1], and setting a fixed value Step as a Step length;
C. updating t to t + Step;
D. obtaining the coordinates of the point U according to the following formula;
in the formula, x*、y*Respectively representing the coordinates of the first point in the 3 points;
E. and C, returning to execute the step C until all the U points are obtained, and connecting all the U points into a curve, namely the generated smooth curve.
9. The method of modeling an intrusive feature in a sliced geological profile as defined in claim 1, wherein: the step (7) specifically comprises:
(7-1) reading a pair of invaded rock layer inferred formation curves GL from the invaded rock layer inferred formation curve set GLpAnd glp+1；
(7-2) inferring a formation curve gl from the invaded rockpAnd glp+1Dividing the original ground level;
(7-3) obtaining the horizon gl based on the graph space relationpAnd glp+1The split horizon in between;
(7-4) combining all the divided stratum surfaces to construct an invaded rock stratum which is used as an invaded structure stratum obtained by modeling;
(7-5) performing steps (7-1) - (7-4) in a loop until the construction of all invasive configurations is completed.
10. An apparatus for modeling an invasive structure in a sectioned 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 |
|---|---|---|---|
| CN202011154215.0A CN112581558B (en) | 2020-10-26 | 2020-10-26 | Model construction method and system for intrusion structure in cut geological section |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011154215.0A CN112581558B (en) | 2020-10-26 | 2020-10-26 | Model construction method and system for intrusion structure in cut geological section |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112581558A true CN112581558A (en) | 2021-03-30 |
| CN112581558B CN112581558B (en) | 2024-04-19 |
Family
ID=75119902
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011154215.0A Active CN112581558B (en) | 2020-10-26 | 2020-10-26 | Model construction method and system for intrusion structure in cut geological section |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112581558B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113326606A (en) * | 2021-05-08 | 2021-08-31 | 北京金阳普泰石油技术股份有限公司 | Control line-based stratigraphic structure fluctuation representation processing method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107767429A (en) * | 2016-08-18 | 2018-03-06 | 阿里巴巴集团控股有限公司 | Curve generation method and equipment |
| CN110163965A (en) * | 2019-05-07 | 2019-08-23 | 南京泛在地理信息产业研究院有限公司 | A kind of solid threedimensional model building method of monoclinal strate tectonic landform |
| CN110910495A (en) * | 2019-10-30 | 2020-03-24 | 南京泛在地理信息产业研究院有限公司 | Three-dimensional geometric form recovery method for dome structure |
| CN111174790A (en) * | 2019-12-03 | 2020-05-19 | 河海大学 | Method for forming topographic profile tracking path |
-
2020
- 2020-10-26 CN CN202011154215.0A patent/CN112581558B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107767429A (en) * | 2016-08-18 | 2018-03-06 | 阿里巴巴集团控股有限公司 | Curve generation method and equipment |
| CN110163965A (en) * | 2019-05-07 | 2019-08-23 | 南京泛在地理信息产业研究院有限公司 | A kind of solid threedimensional model building method of monoclinal strate tectonic landform |
| CN110910495A (en) * | 2019-10-30 | 2020-03-24 | 南京泛在地理信息产业研究院有限公司 | Three-dimensional geometric form recovery method for dome structure |
| CN111174790A (en) * | 2019-12-03 | 2020-05-19 | 河海大学 | Method for forming topographic profile tracking path |
Non-Patent Citations (1)
| Title |
|---|
| 陈志军，等: "基于MapGIS的图切地质剖面系统面向对象设计与实现", 《地质学刊》, vol. 36, no. 3, 28 September 2012 (2012-09-28), pages 257 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113326606A (en) * | 2021-05-08 | 2021-08-31 | 北京金阳普泰石油技术股份有限公司 | Control line-based stratigraphic structure fluctuation representation processing method and system |
| CN113326606B (en) * | 2021-05-08 | 2023-08-15 | 北京金阳普泰石油技术股份有限公司 | Control line-based stratum structure fluctuation characterization processing method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112581558B (en) | 2024-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111899341B (en) | Three-dimensional geological rapid modeling method for long and large banded linear engineering | |
| CN109003330B (en) | Three-dimensional stratum modeling method based on bedrock boundary constraint | |
| AU2013296421B2 (en) | Hybrid method of combining multipoint statistic and object-based methods for creating reservoir property models | |
| CN112381937B (en) | Multi-source geological data coupling modeling method based on drilling and complex geological profile | |
| US9026420B2 (en) | Method and system of updating a geocellular model | |
| CN102385067B (en) | Drawing method for isoline containing reverse fault | |
| AU2016346919B2 (en) | Interactive salt model modification | |
| CN106484943A (en) | A kind of online generation method of the log sheet based on nonlinear fitting control | |
| CN108986213B (en) | Three-dimensional stratum modeling method based on stacking technology | |
| CN107180452A (en) | The stratified geological objects modeling method of the multi-source data coupling on the stratum containing standard | |
| Zhang et al. | A case study on integrated modeling of spatial information of a complex geological body | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN110163965A (en) | A kind of solid threedimensional model building method of monoclinal strate tectonic landform | |
| CN114549774A (en) | Three-dimensional stratum modeling method based on drilling data | |
| CN112581558B (en) | Model construction method and system for intrusion structure in cut geological section | |
| CN109472867B (en) | Method for quantifying influence range of drilling position information | |
| CN116721227A (en) | Automatic modeling method for three-dimensional geologic model of complex geologic body | |
| CN111260791B (en) | Method for updating geosteering model | |
| CN112837413B (en) | Geological drilling-oriented virtual stratum deducing method and device | |
| CN110687591B (en) | Method for determining physical property parameters of coal bed and surrounding rock based on waveform matching of prior data | |
| CN114036608A (en) | Method for converting geological profile into ultimate balance calculation model | |
| CN111951394A (en) | Fault structure unit three-dimensional model construction method and device based on geological map | |
| JP2006009262A (en) | System for creating and displaying geologic structure model, stratum modeling method, and program for creating and displaying geologic structure model | |
| CN109658512A (en) | The three-dimensional more volume modeling methods of shallow seismic profile | |
| CN112734926B (en) | Automatic generation method of cut geological section for loose layer coverage area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |