Method and device for constructing three-dimensional model of fault structural unit based on geological map

Abstract

The present invention discloses a method and device for constructing a three-dimensional model of a fault tectonic unit based on a geological map, the method comprising: first, reading fault data, stratum data and bedrock surface DEM data; second, segmenting the stratum boundary line based on the stratum adjacency relationship, and calculating the occurrence respectively; then, constructing a three-dimensional model of the upper top surface, side surface and lower bottom surface of the stratum based on the bedrock surface data and the stratum boundary line segmentation; finally, combining the stratum three-dimensional model, binding the material, and exporting the model file. The present invention has higher accuracy, can effectively reduce the modeling complexity of the fault tectonic unit, and improve the degree of automation.

Description

Technical Field

[0001] The present invention relates to the technical field of three-dimensional geological modeling, and in particular to a method and device for constructing a three-dimensional model of a fault structural unit based on a geological map. Background Art

[0002] Three-dimensional geological modeling has very important research significance and application value in many fields such as urban planning, engineering construction, oil and gas storage, and digital mines. Traditional geological modeling methods are mainly based on geological survey data such as boreholes and geological profiles to accurately reconstruct and express three-dimensional geological information. However, this method relies too much on modeling data. When the data is limited, the model accuracy is low or even cannot be constructed. Especially in bedrock areas, borehole and geological profile data are missing or sparse, and traditional geological modeling technology is difficult to use. Compared with borehole and profile data, geological maps are easier to obtain and have a wider coverage. They can provide the geological structure information and stratigraphic information required for modeling in the absence of other geological data in bedrock areas. Summary of the invention

[0003] Purpose of the invention: In view of the problems existing in the prior art, the present invention provides a method and device for constructing a three-dimensional model of a fault structural unit based on a geological map with higher accuracy.

[0004] Technical solution: The method for constructing a three-dimensional model of a fault structure unit based on a geological map of the present invention comprises:

[0005] (1) Read the stratigraphic data, fault data and bedrock surface DEM data into stratigraphic set G, fault set F and pixel set E respectively;

[0006] (2) Read any stratum g from the stratum set G i , in g i The boundary line and g i The intersection of the boundary lines of adjacent strata is the dividing point. i The boundary line segments are obtained to obtain a boundary line segment set S, and the inclination and dip angle of each boundary line segment are calculated;

[0007] (3) Based on the formation g i As the calculation range, a triangulated network is constructed based on the bedrock surface DEM data to obtain the stratum g i The upper triangle face set VT;

[0008] (4) Generate the stratum g based on the pixel set E, the boundary line segment set S, and the inclination and dip angle of the corresponding boundary line segment i The top three-dimensional boundary line and the bottom three-dimensional boundary line of the stratum are constructed according to the bottom three-dimensional boundary line. i The lower base triangle face set VB;

[0009] (5) According to the formation g i The three-dimensional boundary line of the upper top surface and the three-dimensional boundary line of the lower bottom surface are used to construct the side triangle face set VS;

[0010] (6) Combine the top triangle set VT, the bottom triangle set VB and the side triangle set VS to form the stratum g i 3D model of

[0011] (7) Execute steps (2) to (6) repeatedly to complete the construction of the three-dimensional model of all strata;

[0012] (8) Combine the three-dimensional models of all strata to obtain a three-dimensional model of the fault tectonic unit.

[0013] Furthermore, step (1) specifically includes:

[0014] (1-1) Read the stratum data to the stratum set G = {g i |i=1,2,…,gn}; where g i represents the i-th stratum, gn represents the number of strata;

[0015] (1-2) Read the slice data into slice set F = {f k |k=1,2,…,fn}; where f k represents the kth fault, and fn represents the number of faults;

[0016] (1-3) Read the bedrock surface DEM data to the pixel set E = {e pq |p=1,2,…,pn,q=1,2,…,qn}; where e pq Represents the pixel in the pth row and qth column of the bedrock surface DEM data, pn is the number of rows of the bedrock surface DEM data, and qn is the number of columns of the bedrock surface DEM data.

[0017] Furthermore, step (2) specifically includes:

[0018] (2-1) Read any stratum g from the stratum set G i , and get its boundary line gl i ;

[0019] (2-2) Reading the formation g i Any adjacent stratum g j , get g j The boundary line of gl j , and calculate the boundary line gl i With gl j The intersection points are stored in the intersection set P;

[0020] (2-3) Repeat step (2-2) until the formation gi All adjacent layers of are traversed to obtain the intersection point set P;

[0021] (2-4) Take the points in P as the dividing points and divide the boundary line gl i Segment, get the boundary line segment set S = {s u |u=1,2,…,sn}; where sn is the boundary line gl i The number of segments;

[0022] (2-5) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula: u tendency;

[0023]

[0024] u=1,2,...,sn

[0025] Among them, dir u For u The tendency of fdir k is the kth fault f in the fault set F k The tendency of gdir i 、gdir j Formation g i , g j The tendency of j g i any adjacent stratum;

[0026] (2-6) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula u The inclination angle;

[0027]

[0028] u=1,2,...,sn

[0029] Among them, dip u For u The inclination angle, fdip k Fault f k The inclination angle, gdip i 、gdip j Formation g i , g j The inclination angle.

[0030] Furthermore, step (4) specifically includes:

[0031] (4-1) Reading the formation g i Boundary line gl i Any boundary line segment s u , and extract su All vertices form a vertex set s u ={s u,v |v=1,2,…,sun},s u,v Indicates u The vth inflection point in , sun is the number of inflection points;

[0032] (4-2) Based on set s u And pixel set E, calculate the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u ;

[0033] (4-3) According to the three-dimensional boundary line of the upper top surface, segment st u and boundary line segments s u The inclination and dip angle of the bottom surface are deduced to be the three-dimensional boundary line segment bt u ;

[0034] (4-4) Loop through steps (4-1) to (4-3) until the boundary line segment set S is traversed and all bottom three-dimensional boundary line segments bt u Merge to get the three-dimensional boundary line gb of the bottom surface i , segment all the top three-dimensional boundary lines into st u Merge to get the top three-dimensional boundary line gt i ;

[0035] (4-5) Extract the three-dimensional boundary line gb of the bottom surface i The turning points form the set gb i＝ {gb i,a |a=1,2,…,bn}; where gb i,a for gb i The ath vertex of bn is the number of vertex points; extract the three-dimensional boundary line gt of the top surface i The turning points form the set gt i＝ {gt i,a |a=1,2,…,bn}; where gt i,a For gt i The a-th inflection point of ;

[0036] (4-6) follow gb i,1 、gb i,a 、gb i,a+1 The vertex order of constructing the stratum g i The lower base triangle vb a {gb i,1 、gb i,a 、gb i,a+1}, forming the bottom triangle face set VB = {vb a|a=2,…,bn-1}.

[0037] Further, step (4-2) includes:

[0038] (4-2-1) Create an empty three-dimensional point set st u ;

[0039] (4-2-2) Based on the pixel set E, the three-dimensional point st is calculated according to the following formula u,v (x v ,y v ,z v ), and join the set st u ;

[0040]

[0041] Among them, (x u,v ,y u,v ) is s u,v The coordinates of v ,q v Point s u,v The row and column numbers of the corresponding DEM pixels, e yvqv For u,v The corresponding DEM pixel value, (X, Y) is the origin coordinate of the DEM pixel, C is the DEM pixel size, is the floor rounding symbol;

[0042] (4-2-3) Set the st u The line formed by all points is used as the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u .

[0043] Furthermore, step (4-3) specifically includes:

[0044] (4-3-1) Based on boundary line segmentation s u The inclination and dip of the fault plane are calculated according to the following formula:

[0045]

[0046] Among them, dir u For u The tendency of dip u For u The inclination angle;

[0047] (4-3-2) Based on the three-dimensional boundary line segmentation of the upper top surface u , the three-dimensional boundary line segment bt of the lower bottom surface of the fault plane is calculated according to the following formula u ;

[0048]

[0049] bt u ={bt u,v |v=1,2,...,sun}

[0050] In the formula, st u,v for st u The vth inflection point of , D is the deduced depth of the bottom surface of the preset stratum.

[0051] Furthermore, step (5) specifically includes:

[0052] (5-1) Obtain the three-dimensional boundary line gb of the bottom surface i＝ {gb i,a |a=1,2,…,bn} and the top three-dimensional boundary line gt i ={gt i,a |a=1,2,…,bn}, where gb i,a for gb i The a-th inflection point, gt i,a For gt i The a-th inflection point of , bn is the number of inflection points;

[0053] (5-2) Follow gt i,a 、gb i,a ,gt i,a+1 Vertex order to construct the side regular triangle vsl a {gt i,a 、gb i,a ,gt i,a+1}, forming a regular triangle face set VSL = {vsl a |a=2,…,bn-1};

[0054] (5-3) According to gb i,a+1 、gb i,a ,gt i,a+1 Vertex order to construct the inverted triangle surface vsr on the side a {gb i,a+1 、gb i,a ,gt i,a+1}, forming an inverted triangle face set VSR = {vsr a |a=2,…,bn-1};

[0055] (5-4) Integrate the regular triangle face set VSL and the inverted triangle face set VSR as the formation g i The side triangle face set VS.

[0056] Furthermore, step (8) specifically includes:

[0057] (8-1) Create a stratum model material according to the preset texture and bind it to the corresponding stratum 3D model;

[0058] (8-2) Combine the three-dimensional models of all strata to obtain a three-dimensional model of the fault tectonic unit.

[0059] The device for constructing a three-dimensional model of a fault structural unit based on a geological map of the present invention comprises a processor and a computer program stored in a memory and executable on the processor, and the method is implemented when the processor executes the program.

[0060] Beneficial effects: Compared with the prior art, the present invention has the following significant advantages: the present invention has higher accuracy, avoids excessive reliance on geological survey data such as drilling holes and profiles, and has a higher degree of automation. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] Figure 1 It is the DEM data of strata, faults and bedrock surface used in this embodiment;

[0062] Figure 2 is a flow chart of an embodiment of the present invention;

[0063] Figure 3 is the formation occurrence data of this embodiment;

[0064] Figure 4 is the fault occurrence data of this embodiment;

[0065] Figure 5 is an example of a three-dimensional stratum model constructed in this embodiment;

[0066] Figure 6 is the three-dimensional model of the fault unit constructed in this embodiment;

[0067] Figure 7 It is the material binding result of the three-dimensional model of the fault unit in this embodiment. DETAILED DESCRIPTION

[0068] The technical solution of the present invention is further described in detail below. This embodiment selects the bedrock surface DEM data of Tangshan, Nanjing ( Figure 1 ) is used as experimental data, and the projection coordinate system used in the experimental data is Nanjing 92 coordinate system. The following is further explained by describing a specific embodiment in conjunction with the accompanying drawings.

[0069] like Figure 2 As shown, this embodiment provides a method for three-dimensional modeling of fault structural units based on a geological map, which specifically includes the following steps:

[0070] (1) Read the stratigraphic data, fault data and bedrock surface DEM data into stratigraphic set G, fault set F and pixel set E respectively.

[0071] This step specifically includes:

[0072] (1-1) Read the stratum data to the stratum set G = {g i |i=1,2,…,gn}; where g i represents the i-th stratum, gn represents the number of strata; in this embodiment, gn=62;

[0073] (1-2) Read the slice data into slice set F = {f k |k=1,2,…,fn}; where f k represents the kth fault, and fn represents the number of faults; in this embodiment, fn=17;

[0074] (1-3) Read the bedrock surface DEM data to the pixel set E = {e pq |p=1,2,…,pn,q=1,2,…,qn}; where e pq It represents the pixel in the p-th row and q-th column in the bedrock surface DEM data, pn is the number of rows of the bedrock surface DEM data, qn is the number of columns of the bedrock surface DEM data, and in this embodiment, pn=8313, qn=15226.

[0075] (2) Read any stratum g from the stratum set G i , in g i The boundary line and g i The intersection of the boundary lines of adjacent strata is the dividing point. i The boundary line segments are obtained to obtain a set S of boundary line segments, and the inclination and dip angle of each boundary line segment are calculated.

[0076] This step specifically includes:

[0077] (2-1) Read any stratum g from the stratum set G i , and get its boundary line gl i ;

[0078] (2-2) Reading the formation g i Any adjacent stratum g j , get g j The boundary line of gl j , and calculate the boundary line gl i With gl j The intersection points are stored in the intersection set P;

[0079] (2-3) Repeat step (2-2) until the formation g i All adjacent layers of are traversed to obtain the intersection point set P;

[0080] (2-4) Take the points in P as the dividing points and divide the boundary line gl i Segment, get the boundary line segment set S = {s u |u=1,2,…,sn}; where sn is the boundary line gl i The number of segments;

[0081] (2-5) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula: u tendency;

[0082]

[0083] u=1,2,...,sn

[0084] Among them, dir u For u The tendency of fdir k is the kth fault f in the fault set F k The tendency of gdir i 、gdir j Formation g i , g j The tendency of j g i any adjacent stratum;

[0085] (2-6) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula u The inclination angle;

[0086]

[0087] u=1,2,...,sn

[0088] Among them, dip u For u The inclination angle, fdip k Fault f k The inclination angle, gdip i 、gdip j Formation g i , g j In this embodiment, the dip and inclination data of the strata and faults are as follows: Figure 3 and Figure 4 shown.

[0089] (3) Based on the triangulation interface of GIS software, the stratum g i As the calculation range, a triangulated network is constructed based on the bedrock surface DEM data to obtain the stratum g i The top triangular face set VT.

[0090] (4) Generate the stratum g based on the pixel set E, the boundary line segment set S, and the inclination and dip angle of the corresponding boundary line segment i The top three-dimensional boundary line and the bottom three-dimensional boundary line of the stratum are constructed according to the bottom three-dimensional boundary line. i The lower base triangle face set VB.

[0091] This step specifically includes:

[0092] (4-1) Reading the formation g i Boundary line gl i Any boundary line segment s u , and extract s u All vertices form a vertex set s u ={s u,v |v=1,2,…,sun},s u,v Indicates u The vth inflection point in , sun is the number of inflection points;

[0093] (4-2) Based on set s u And pixel set E, calculate the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u ;

[0094] (4-3) According to the three-dimensional boundary line of the upper top surface, segment st u and boundary line segments s u The inclination and dip angle of the bottom surface are deduced to be the three-dimensional boundary line segment bt u ;

[0095] (4-4) Loop through steps (4-1) to (4-3) until the boundary line segment set S is traversed and all bottom three-dimensional boundary line segments bt u Merge to get the three-dimensional boundary line gb of the bottom surface i , segment all the top three-dimensional boundary lines into st u Merge to get the top three-dimensional boundary line gt i ;

[0096] (4-5) Extract the three-dimensional boundary line gb of the bottom surface i The turning points form the set gb i＝ {gb i,a |a=1,2,…,bn}; where gb i,a for gb i The ath vertex of bn is the number of vertex points; extract the three-dimensional boundary line gt of the top surface i The turning points form the set gt i＝ {gt i,a|a=1,2,…,bn}; where gt i,a For gt i The a-th inflection point of ;

[0097] (4-6) follow gb i,1 、gb i,a 、gb i,a+1 The vertex order of constructing the stratum g i The lower base triangle vb a {gb i,1 、gb i,a 、gb i,a+1}, forming the bottom triangle face set VB = {vb a |a=2,…,bn-1}.

[0098] Further, step (4-2) includes:

[0099] (4-2-1) Create an empty three-dimensional point set st u ;

[0100] (4-2-2) Based on the pixel set E, the three-dimensional point st is calculated according to the following formula u,v (x v ,y v ,z v ), and join the set st u ;

[0101]

[0102] Among them, (x u,v ,y u,v ) is s u,v The coordinates of v ,q v Point s u,v The row and column numbers of the corresponding DEM pixels, For u,v The corresponding DEM pixel value, (X, Y) is the origin coordinate of the DEM pixel, C is the DEM pixel size, is a rounding down symbol; in this embodiment, X=172515.64, Y=56491.81, C=10;

[0103] (4-2-3) Set the st u The line formed by all points is used as the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u .

[0104] Wherein, step (4-3) specifically includes:

[0105] (4-3-1) Based on boundary line segmentation s u The inclination and dip of the fault plane are calculated according to the following formula:

[0106]

[0107] Among them, dir u For u The tendency of dip u For u The inclination angle;

[0108] (4-3-2) Based on the three-dimensional boundary line segmentation of the upper top surface u , the three-dimensional boundary line segment bt of the lower bottom surface of the fault plane is calculated according to the following formula u ;

[0109]

[0110] bt u ={bt u,v |v=1,2,...,sun}

[0111] In the formula, st u,v for st u The vth inflection point of , D is the deduced depth of the bottom surface of the preset stratum, and in this embodiment, D=500 meters.

[0112] (5) According to the formation g i The three-dimensional boundary line of the upper top surface and the three-dimensional boundary line of the lower bottom surface are used to construct the side triangle face set VS.

[0113] This step specifically includes:

[0114] (5-1) Obtain the three-dimensional boundary line gb of the bottom surface i = {gb i,a |a=1,2,…,bn} and the top three-dimensional boundary line gt i ={gt i,a |a=1,2,…,bn}, where gb i,a for gb i The a-th inflection point, gt i,a For gt i The a-th inflection point of , bn is the number of inflection points;

[0115] (5-2) Follow gt i,a 、gb i,a ,gt i,a+1 Vertex order to construct the side regular triangle vsl a {gt i,a 、gb i,a ,gti,a+1}, forming a regular triangle face set VSL = {vsl a |a=2,…,bn-1};

[0116] (5-3) According to gb i,a+1 、gb i,a ,gt i,a+1 Vertex order to construct the inverted triangle surface vsr on the side a {gb i,a+1 、gb i,a ,gt i,a+1}, forming an inverted triangle face set VSR = {vsr a |a=2,…,bn-1};

[0117] (5-4) Integrate the regular triangle face set VSL and the inverted triangle face set VSR as the formation g i The side triangle face set VS.

[0118] (6) Combine the top triangle set VT, the bottom triangle set VB and the side triangle set VS to form the stratum g i In this embodiment, the three-dimensional model of the stratum corresponding to the stratum g1 is as follows: Figure 5 shown.

[0119] (7) Execute steps (2) to (6) repeatedly to complete the construction of the three-dimensional models of all strata. Figure 6 shown.

[0120] (8) Combining the three-dimensional models of all strata to obtain a three-dimensional model of the fault structural unit. In this embodiment, the three-dimensional model of the structural unit is as follows: Figure 7 shown.

[0121] This step specifically includes:

[0122] (8-1) Based on the 3D modeling software API, create a stratum model material according to the preset texture and bind it to the corresponding stratum 3D model;

[0123] (8-2) Combine the three-dimensional models of all strata to obtain the three-dimensional model of the fault structural unit and export it as a model file.

[0124] In this embodiment, the triangular facets of the top surface of the stratum are calculated based only on the triangulation network construction interface provided by ArcGIS software. The method can also use the interface of GIS software such as SuperMap and QGIS. In this embodiment, the material is created based only on the three-dimensional model editing software 3ds Max. The method can also use the interface of software such as Maya and AutoCad. In this embodiment, the three-dimensional stratum model is exported only in fbx format. The method can also export the three-dimensional stratum model in other formats such as OBJ.

[0125] This embodiment also provides a device for constructing a three-dimensional model of a fault structural unit based on a geological map, including a processor and a computer program stored in a memory and executable on the processor, wherein the processor implements the above method when executing the program.

Claims

1. A method for constructing a three-dimensional model of a fault tectonic unit based on a geological map. Features The method includes: (1) Read the stratigraphic data, fault data and bedrock surface DEM data into stratigraphic set G, fault set F and pixel set E respectively; (2) Read any stratum g from the stratum set G i , in g i The boundary line and g i The intersection of the boundary lines of adjacent strata is the dividing point. i The boundary line segments are obtained to obtain a boundary line segment set S, and the inclination and dip angle of each boundary line segment are calculated; (3) Based on the formation g i As the calculation range, a triangulated network is constructed based on the bedrock surface DEM data to obtain the stratum g i The upper triangle face set VT; (4) Generate the stratum g based on the pixel set E, the boundary line segment set S, and the inclination and dip angle of the corresponding boundary line segment i The top three-dimensional boundary line and the bottom three-dimensional boundary line of the stratum are constructed according to the bottom three-dimensional boundary line. i The lower base triangle face set VB; (5) According to the formation g i The three-dimensional boundary line of the upper top surface and the three-dimensional boundary line of the lower bottom surface are used to construct the side triangle face set VS; (6) Combine the top triangle set VT, the bottom triangle set VB and the side triangle set VS to form the stratum g i 3D model of (7) Execute steps (2) to (6) repeatedly to complete the construction of the three-dimensional model of all strata; (8) combining the three-dimensional models of all strata to obtain a three-dimensional model of the fault tectonic unit; Step (2) specifically includes: (2-1) Read any stratum g from the stratum set G i , and get its boundary line gl i ; (2-2) Reading the formation g i Any adjacent stratum g j , get g j The boundary line of gl j , and calculate the boundary line gl i With gl j The intersection points are stored in the intersection set P; (2-3) Repeat step (2-2) until the formation g i All adjacent layers of are traversed to obtain the intersection point set P; (2-4) Take the points in P as the dividing points and divide the boundary line gl i Segmentation, get the boundary line segment set S = {s u |u=1,2,…,sn}; where sn is the boundary line gl i The number of segments; (2-5) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula: u tendency; Among them, dir u For u The tendency of fdir k is the kth fault f in the fault set F k The tendency of gdir i 、gdir j Formation g i , g j The tendency of j g i any adjacent stratum; (2-6) Based on the dip data of the strata and faults, the segment s of each boundary line is calculated according to the following formula u The inclination angle; Among them, dip u For u The inclination angle, fdip k Fault f k The inclination angle, gdip i 、gdip j The formations gi and g are j The inclination angle; Step (4) specifically includes: (4-1) Reading the formation g i Boundary line gl i Any boundary line segment s u , and extract s u All vertices form a vertex set s u ={s u,v |v=1,2,…,sun},s u,v Indicates u The vth inflection point in , sun is the number of inflection points; (4-2) Based on set s u And pixel set E, calculate the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u ; (4-3) According to the three-dimensional boundary line of the upper top surface, segment st u and boundary line segments s u The inclination and dip angle of the bottom surface are deduced to be the three-dimensional boundary line segment bt u ; (4-4) Loop through steps (4-1) to (4-3) until the boundary line segment set S is traversed and all bottom surface three-dimensional boundary line segments bt u Merge to get the three-dimensional boundary line gb of the bottom surface i , segment all the top three-dimensional boundary lines into st u Merge to get the top three-dimensional boundary line gt i ; (4-5) Extract the three-dimensional boundary line gb of the bottom surface i The turning points form the set gb i＝ {gb i,a |a=1,2,…,bn}; where gb i,a for gb i The ath vertex of bn is the number of vertex points; extract the three-dimensional boundary line gt of the top surface i The turning points form the set gt i＝ {gt i,a |a=1,2,…,bn}; where gt i,a For gt i The a-th inflection point of ; (4-6) follow gb i,1 、gb i,a 、gb i,a+1 The vertex order of constructing the stratum g i The lower base triangle vb a {gb i,1 、gb i,a 、gb i,a+1 }, forming the lower base triangle face set VB = {vb a |a=2,…,bn-1}; Step (4-2) includes: (4-2-1) Create an empty three-dimensional point set st u ; (4-2-2) Based on the pixel set E, the three-dimensional point st is calculated according to the following formula u,v (x v ,y v ,z v ), and join the set st u ; Among them, (x u,v ,y u,v ) is s u,v The coordinates of v ,q v Point s u,v The row and column numbers of the corresponding DEM pixels, For u,v The corresponding DEM pixel value, (X, Y) is the origin coordinate of the DEM pixel, C is the DEM pixel size, is the floor rounding symbol; (4-2-3) Set the st u The line formed by all points is used as the boundary line segment s u The corresponding upper three-dimensional boundary line segment st u ; Step (4-3) specifically includes: (4-3-1) Based on boundary line segmentation s u The inclination and dip of the fault plane are calculated according to the following formula: Among them, dir u For u The tendency of dip u For u The inclination angle; (4-3-2) Based on the three-dimensional boundary line segmentation of the upper top surface u , the three-dimensional boundary line segment bt of the lower bottom surface of the fault plane is calculated according to the following formula u ; bt u ={bt u,v |v=1,2,...,sun} In the formula, st u,v for st u The vth inflection point of , D is the deduced depth of the bottom surface of the preset stratum.

2. The method for constructing a three-dimensional model of a fault structural unit based on a geological map according to claim 1, Features: Step (1) specifically includes: (1-1) Read the stratum data to the stratum set G = {g i |i=1,2,…,gn}; where g i represents the i-th stratum, gn represents the number of strata; (1-2) Read the slice data into slice set F = {f k |k=1,2,…,fn}; where f k represents the kth fault, and fn represents the number of faults; (1-3) Read the bedrock surface DEM data to the pixel set E = {e pq |p=1,2,…,pn,q=1,2,…,qn}; where e pq Represents the pixel in the pth row and qth column of the bedrock surface DEM data, pn is the number of rows of the bedrock surface DEM data, and qn is the number of columns of the bedrock surface DEM data.

3. The method for constructing a three-dimensional model of a fault structural unit based on a geological map according to claim 1, Features: Step (5) specifically includes: (5-1) Obtain the three-dimensional boundary line gb of the bottom surface i＝ {gb i,a |a=1,2,…,bn} and the top three-dimensional boundary line gt i＝ {gt i,a |a=1,2,…,bn}, where gb i,a for gb i The a-th inflection point, gt i,a For gt i The a-th inflection point of , bn is the number of inflection points; (5-2) Follow gt i,a 、gb i,a ,gt i,a+1 Vertex order to construct the side regular triangle vsl a {gt i,a 、gb i,a ,gt i,a+1 }, forming a regular triangle face set VSL = {vsl a |a=2,…,bn-1}; (5-3) According to gb i,a+1 、gb i,a ,gt i,a+1 Vertex order to construct the inverted triangle surface vsr on the side a {gb i,a+1 、gb i,a ,gt i,a+1 }, forming an inverted triangle face set VSR = {vsr a |a=2,…,bn-1}; (5-4) Integrate the regular triangle face set VSL and the inverted triangle face set VSR as the formation g i The side triangle face set VS.

4. The method for constructing a three-dimensional model of a fault structure unit based on a geological map according to claim 1, Features: Step (8) specifically includes: (8-1) Create a stratum model material according to the preset texture and bind it to the corresponding stratum 3D model; (8-2) Combine the three-dimensional models of all strata to obtain a three-dimensional model of the fault tectonic unit.

5. A device for constructing a three-dimensional model of a fault tectonic unit based on a geological map, comprising a processor and a computer program stored in a memory and executable on the processor, Features: When the processor executes the program, the method according to any one of claims 1 to 4 is implemented.

Images