An anticline structure three-dimensional modeling method based on a plane geological map

Abstract

The application discloses an anticline structure three-dimensional modeling method based on a plane geological map, belongs to the geological field, and is used for solving the problem that the spatial form of an anticline can be presented only by collecting data and conducting field investigation manually for a long time, the method is time-consuming and laborious, and the topographic features on the anticline cannot be presented one by one. The anticline structure three-dimensional modeling method specifically comprises the following steps: grid data processing, intelligent determination of anticline folds, entity space subdivision and modeling, and calculation of denudation amount. The application can quickly and accurately present the spatial form of the anticline structure and the topographic features on the anticline.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of geology and relates to a three-dimensional modeling technology of anticline structure, in particular to a three-dimensional modeling method of anticline structure based on a plane geological map. BACKGROUND

[0002] An anticline refers to a shape that is convex upward when strata are folded, and on a general flat ground, the upper half of the strata of the anticline is eroded and flattened, which forms a stratum arrangement mode in which the middle is old and the two sides are relatively new. The anticline and syncline have important practical significance. Since the anticline strata are convex upward, and the density of oil and gas is smaller than that of water, the anticline is often a good oil and gas storage structure. The anticline is convex upward due to the rock mass, has a supporting effect, and is not easy to seep water, so it is suitable for tunneling under the anticline.

[0003] In the prior art, a large amount of time is spent by people to collect data and conduct field survey, so that the spatial form of the anticline can be presented. Such a method is time-consuming and laborious, and cannot present the denudation and other features on the anticline. Therefore, the three-dimensional modeling method of anticline structure based on a plane geological map is provided. SUMMARY

[0004] In view of the deficiencies in the prior art, the application aims to provide a three-dimensional modeling method of anticline structure based on a plane geological map.

[0005] The technical problem to be solved by the application is how to quickly and accurately present the spatial form of the anticline structure and the topographic features on the anticline.

[0006] The object of the application can be achieved by the following technical scheme.

[0007] The three-dimensional modeling method of anticline structure based on a plane geological map specifically comprises the following steps.

[0008] Step S1: grid data processing.

[0009] Step S2: intelligent determination of anticline folds.

[0010] Step S3: entity space subdivision and modeling.

[0011] Step S4: calculation of denudation amount.

[0012] Further, the grid data processing process specifically comprises the following steps.

[0013] Step S11: extraction of terrain contour lines.

[0014] The terrain contour lines are extracted according to the terrain features on the terrain geological map and are saved in a vector format, such as a wl format of mapgis or a line file of shapefile, and are assigned with an elevation attribute.

[0015] Step S12: Constructing a digital elevation model;

[0016] The topographic contour line vector file obtained by step S11 is combined with the grid parameters given by the user, and is subjected to gridding processing to generate a digital elevation model and is saved in geotiff format;

[0017] Step S13: Extracting rock layer properties and occurrence;

[0018] According to the geological information given by the topographic and geological map, each independent geological unit is vectorized and given a stratum attribute, and is saved as a wp or shapefile surface file of mapgis;

[0019] The occurrence position given in the topographic and geological map is vectorized and saved as a wt or shapefile point file of mapgis, and is given a corresponding occurrence and rock layer attribute;

[0020] Step S14: Extracting geological boundaries;

[0021] According to the geological information given by the topographic and geological map, the rock layer connection is vectorized and given a boundary attribute.

[0022] Further, the grid parameter is a pixel ground resolution;

[0023] The boundary attribute is the rock layer attribute on both sides of the boundary.

[0024] Further, the intelligent determination process of the anticline fold is specifically as follows:

[0025] Step S21: Initialization of spatial elements;

[0026] Modeling range: Bbox = {‘west’: west, ‘east’: east, ‘south’: south, ‘north’: north};

[0027] Digital elevation model: Dem.tif;

[0028] Geological unit: units = {u_1: {‘name’: name_1, ‘era’: era, ‘points’: points_1}, ……};

[0029] Geological boundary: boundaries = {b_1: {‘name’: name_1, ‘points’: points_1, ‘left’: u_m, ‘right’: u_n}, ……};

[0030] Geological occurrence:occurrence={u_1: [[x, y, dip, angle], …], …};

[0031] Step S22: key signs of anticline;

[0032] Step S23: determination of core and two wings;

[0033] According to the definition of anticline structure, according to the obtained geological boundary, as the basic boundary element of the anticline fold, the boundary located at the first position of the set bdy3 is the outermost boundary of the two wings of the anticline ly=[b01, b10], and the geological boundary located at the center of bdy3 is the core of the anticline c=[b0n, bn0].

[0034] Further, the analysis process of step S22 is as follows:

[0035] Step S221: consistency of geological boundary trend;

[0036] Iterate the geological boundary results obtained in step S21, calculate the trend of all geological boundaries, and select the geological boundaries with consistent trend;

[0037] bdy1=[b1, b2, b3……, bn];

[0038] Step S222: uniqueness of geological unit age-new-old-new arrangement;

[0039] Iterate the consistent geological boundary bdy1 obtained in step S21, detect the geological unit age on both sides of the boundary, determine that the geological unit age arrangement follows the new-old-new change rule, remove the geological boundaries that do not meet the constraint, and obtain the result;

[0040] bdy2=[b01, b02, b03, ……b0n];

[0041] Step S223: symmetry of geological unit arrangement;

[0042] Iterate the result bdy2 obtained in step S222, detect whether the geological units on both sides of the geological boundary present spatial arrangement symmetry, and obtain the boundary elements of the structural anticline after removing the geological boundaries that do not meet the constraint condition;

[0043] bdy3=[b01, b02, b03, ……b0n, bn0……, b30, b20, b10].

[0044] Further, the rock layer on the right side of b01 and the left side of b10 is approximately axisymmetric and belongs to the same rock layer unit, and b01 and b10 are the boundary constraints of the anticline fold modeling.

[0045] Further, the entity space profiling and modeling process is specifically as follows:

[0046] Step S31: Anticline constraint boundary calculation;

[0047] According to the geological boundary ly of the two wings, read the occurrence of the b01 and b10 boundary constraint strata, respectively read the coordinates of K01, L01 and K10, L10 two points of b01 and b10 geological boundary;

[0048] And according to the occurrence of b01.left and b10.right strata, construct plane K01L01Q1P1 and plane K10L10Q1P1, K10L10Q1P1 plane is consistent with the dip and dip angle of the rock stratum B on the left side of the geological boundary b01, K01L01Q1P1 plane is consistent with the dip and dip angle of the rock stratum B on the right side of the geological boundary b02;

[0049] The plane K01L01Q1P1 and the plane K10L10Q1P1 intersect at the straight line P1Q1, and the coordinates Q1 (Q1x, Q1y, Q1z) and P1 (P1x, P1y, P1z) of the points Q1P1 are calculated, and the higher value of the elevation of the points Q1 and P1 is taken as the model vertex;

[0050] The face K01L01Q1P1 and the face K10L10Q1P1 constitute the peripheral constraint boundary of the anticline;

[0051] Step S32: stratum constraint interface construction;

[0052] According to the method of step S31 for obtaining the peripheral constraint boundary of the anticline, read the K2, L2 and M2, N2 coordinates of the geological boundary b02 and b20, and read the occurrence of the rock stratum on the right side of b02 and the left side of b20;

[0053] According to the dip and dip angle, respectively, construct the faces K02L02Q2P2 and K20L20Q2P2, and the two faces intersect at the straight line P2Q2, and the constraint interface control point coordinates Q2 (Q2x, Q2y, Q2z) and P2 (P2x, P2y, P2z) of the rock stratum C are calculated respectively;

[0054] Using the same method, according to the parameters and spatial attributes of the remaining geological boundaries in bdy3, until the constraint interface control point coordinates Qn (Qnx, Qny, Qnz) and Pn (Pnx, Pny, Pnz) of the rock stratum N are completed;

[0055] Step S33: cubic profiling and denudation strata three-dimensional reconstruction;

[0056] According to the set minimum subdivision unit cube edge length, taking L01 of the geological boundary b01 as a starting point, searching in horizontal direction and vertical direction according to the digital elevation model dem.tif and the stratum unit units and the geological boundary bdy3 respectively, and taking the B stratum constraint boundary on the right side of b01 as a condition;

[0057] If the current minimum unit center store coordinate (x, y, z) is located in the space surrounded by the constraint interface K01L01Q1P1L10K10 and K02L02Q2P2L20K20, then the attribute of the unit is assigned as B stratum;

[0058] Similarly, until all the cubic units belonging to the stratum N are searched, different colors are given to different cubic subdivision units according to different stratum attributes, that is, the spatial form of the anticline fold of the eroded part is directly presented.

[0059] Further, the working process of the erosion amount calculation is as follows:

[0060] By searching the cubic subdivision units with different stratum attributes, the number thereof is counted, and the approximate volume of the restored eroded stratum is obtained according to the volume of the minimum cubic subdivision unit.

[0061] Compared with the prior art, the present application has the following beneficial effects:

[0062] The anticline structure in the present application does not need to spend a lot of time on data collection and field survey, and is processed through grid data, intelligent determination of anticline folds, entity space subdivision and modeling and erosion amount calculation, so that not only the spatial form of the anticline can be presented, but also the eroded topography and other features on the anticline can be presented one by one. BRIEF DESCRIPTION OF DRAWINGS

[0063] In order to facilitate those skilled in the art to understand, the present application will be further described below with reference to the drawings.

[0064] Figure 1 It is a structural schematic view of the geological boundary in the present application.

[0065] Figure 2 It is a structural schematic view of the anticline peripheral constraint boundary in the present application.

[0066] Figure 3 It is a structural schematic view of the structure surface in the present application.

[0067] Figure 4 It is a structural schematic view of the stratum N constraint interface control point coordinate in the present application.

[0068] Figure 5 The working flowchart of the present application. DETAILED DESCRIPTION

[0069] The technical solutions of the present application will be described clearly and completely in connection with the embodiments. Obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present application.

[0070] The topographic and geological map is a map showing the topographic features, strata, ore bed distribution, rock bed occurrence and structural features of the study area.

[0071] A fold is a deformation of various planes in rock, such as bedding planes, cleavage planes, etc., caused by bending under stress. It is a representation of the change of a nearly flat plane in rock into a curved plane. Most of the deformed planes that form folds are bedding planes. Fold is a common geological structure on the earth's crust. It is most obvious in stratified rocks. Some folds are formed like a flat newspaper being squeezed from both sides to the center by two hands. The newspaper will be raised, and when the raised part is too high, the top will bend and collapse. This illustrates the role of two forces in the formation of folds. One is the horizontal compressive force, and the other is its own gravity. In addition, not all folds are raised upward. The fold plane bending upward is called an anticline, and the fold plane bending downward is called a syncline.

[0072] In an embodiment, referring to FIG. 1, a three-dimensional modeling method of anticline structure based on a planar geological map is shown, and the three-dimensional modeling method of anticline structure specifically includes the following steps. Figures 1-5

[0073] Step S1: grid data processing;

[0074] In this embodiment, the grid data processing is specifically as follows:

[0075] Step S11: extracting topographic contour lines;

[0076] The topographic contour lines are extracted according to the topographic features on the topographic and geological map, and are saved in a vector format, such as a wl format of mapgis or a line file storage of shapefile, and are given an elevation attribute;

[0077] Step S12: constructing a digital elevation model;

[0078] The topographic contour line vector file obtained through step S11 is combined with the grid parameters (pixel ground resolution) given by the user to perform gridding processing, generate a digital elevation model, and save it in a geotiff format;

[0079] Step S13: extracting rock layer attributes and occurrence;

[0080] ​According to the geological information given in the topographic geological map, each independent geological unit is vectorized, and the stratum attribute is given, and saved as the wp or shapefile surface file of mapgis;

[0081] The occurrence position given in the topographic geological map is vectorized, and saved as the wt or shapefile point file of mapgis, and given the corresponding occurrence and stratum attribute;

[0082] Step S14: Extract the geological boundary;

[0083] According to the geological information given in the topographic geological map, the stratum connection is vectorized, and the boundary attribute (stratum attribute on both sides of the boundary) is given;

[0084] Step S2: Intelligent determination of anticline fold;

[0085] In this embodiment, the intelligent determination process of the anticline fold is as follows:

[0086] Anticline refers to the shape of the stratum when it is folded, which is convex upward. On a general flat ground, the upper half of the anticline stratum is eroded and flattened, which forms a middle old and both sides relatively new stratum arrangement;

[0087] Step S21: Initialization of spatial elements;

[0088] Modeling range: Bbox={‘west’:west,’east’:east,’south’:south,’north’:north};

[0089] Digital elevation model: Dem.tif;

[0090] Geological unit: units={u_1:{‘name’:name_1,’era’:era,’points’:points_1},……};

[0091] Geological boundary: boundaries={b_1:{‘name’:name_1,’points’:points_1,’left’:u_m,’right’:u_n},……};

[0092] Geological occurrence: occurence={u_1:[[x,y,dip,angle],……],……};

[0093] Step S22: Key sign of anticline;

[0094] Step S221: Consistency of geological boundary trend;

[0095] Traverse the geological boundary results obtained in step S21, calculate the trend of all geological boundaries, and select the geological boundaries with consistent trends;

[0096] bdy1=[b1, b2, b3……, bn];

[0097] Step S222: Geological unit age-new-old-new arrangement uniqueness;

[0098] Traverse the geological boundaries bdy1 with consistent trends obtained in step S21, detect the geological unit ages on both sides of the boundaries, determine that the geological unit age arrangement follows the new-old-new change rule, remove the geological boundaries that do not meet the constraints, and obtain the results;

[0099] bdy2=[b01, b02, b03,……b0n];

[0100] Step S223: Geological unit arrangement symmetry;

[0101] Traverse the results bdy2 obtained in step S222, detect whether the geological units on both sides of the geological boundaries present spatial arrangement symmetry, and obtain the boundary elements of the structural anticline after removing the geological boundaries that do not meet the constraints;

[0102] bdy3=[b01, b02, b03,……b0n, bn0……, b30, b20, b10];

[0103] Wherein, the rock layer inclination on the right side of b01 and the left side of b10 is approximately axisymmetric, and belongs to the same rock layer unit, b01 and b10 are the boundary constraints of the anticline fold modeling.

[0104] Step S23: Core and two wings determination;

[0105] According to the definition of anticline structure, the geological boundaries obtained in step S223 can be used as the basic boundary elements of the anticline fold, the boundary located at the first position of the set bdy3 is the outermost boundary of the two wings of the anticline ly=[b01, b10], and the geological boundary located at the center of bdy3 is the core of the anticline c=[b0n, bn0];

[0106] Step S3: Entity space subdivision and modeling;

[0107] Step S31: Anticline constraint boundary calculation;

[0108] According to the occurrence of the boundary constrained strata of b01 and b10 of the two wings, the coordinates of K01, L01 and K10, L10 of b01 and b10 geological boundaries are read respectively;

[0109] According to the stratum occurrence of b01.left and b10.right, plane K01L01Q1P1 and plane K10L10Q1P1 are constructed, the plane K10L10Q1P1 is consistent with the dip and the angle of stratum B on the left of geological boundary b01, and the plane K01L01Q1P1 is consistent with the dip and the angle of stratum B on the right of geological boundary b02;

[0110] The plane K01L01Q1P1 and the plane K10L10Q1P1 intersect at straight line P1Q1, the coordinates Q1 (Q1x, Q1y, Q1z) and P1 (P1x, P1y, P1z) of points Q1P1 are calculated, the elevations of points Q1 and P1 are compared, and the point with the higher value is taken as the model vertex;

[0111] The plane K01L01Q1P1 and the plane K10L10Q1P1 constitute the peripheral constraint boundary of the anticline Figure 2 );

[0112] Step S32: stratum constraint interface construction;

[0113] According to the method of step S31 for obtaining the peripheral constraint boundary of the anticline, the K2, L2 and M2, N2 coordinates of geological boundaries b02 and b20 are read, and the stratum occurrence on the right of b02 and on the left of b20 is read;

[0114] According to the dip and the angle, the planes K02L02Q2P2 and K20L20Q2P2 are constructed respectively Figure 3 ), the two planes intersect at straight line P2Q2, and the constraint interface control point coordinates Q2 (Q2x, Q2y, Q2z) and P2 (P2x, P2y, P2z) of stratum C are calculated respectively;

[0115] By using the same method, according to the parameters and spatial attributes of the remaining geological boundaries in bdy3, until the constraint interface control point coordinates Qn (Qnx, Qny, Qnz) and Pn (Pnx, Pny, Pnz) of stratum N are completed, see Figure 4 ;

[0116] Step S33: cubic subdivision and denudation stratum three-dimensional reconstruction;

[0117] According to the set minimum subdivision unit cube edge length, taking L01 of the geological boundary b01 as the starting point, searching in the horizontal direction and the vertical direction according to the digital elevation model dem.tif and the stratum unit units and the geological boundary bdy3 respectively, and taking the B stratum constraint boundary on the right side of b01 as the condition, as long as the current minimum unit center coordinates (x, y, z) are located in the space surrounded by the constraint interfaces K01L01Q1P1L10K10 and K02L02Q2P2L20K20, the attribute of the unit is assigned as the B stratum, and similarly, until all the cubic units belonging to the stratum N are searched. According to different stratum attributes, different colors are given to different cubic subdivision units, so that the spatial form of the denuded anticline fold can be intuitively presented.

[0118] Step S4: denudation amount calculation;

[0119] By searching the cubic subdivision units of different stratum attributes, the number thereof is counted, and according to the volume of the minimum cubic subdivision unit, the approximate volume of the restored denuded stratum can be obtained.

[0120] The above formulas are all de-dimensioned to calculate the numerical values, the formulas are obtained by software simulation of a large amount of data to obtain a formula closest to the real situation, the preset parameters in the formula are set by the person skilled in the art according to the actual situation, and the weight coefficient and the proportion coefficient are specific values obtained by quantifying each parameter, which is convenient for subsequent comparison. As for the size of the weight coefficient and the proportion coefficient, as long as the proportion relationship between the parameters and the quantized values is not affected, it is acceptable.

[0121] The preferred embodiments disclosed above are only used to help explain the present application. The preferred embodiments do not describe all the details and do not limit the application to the specific embodiments. Obviously, according to the content of the specification, many modifications and changes can be made. The specification selects and specifically describes these embodiments in order to better explain the principles and practical applications of the application, so that those skilled in the art can well understand and utilize the application. The application is limited only by the claims and their entire scope and equivalents.

Claims

1. A method for three-dimensional modeling of anticline structure based on a planar geologic map, characterized by, The method for three-dimensional modeling of anticline structure specifically comprises the following steps: Step S1: grid data processing; The grid data processing process specifically comprises the following steps: Step S11: extracting terrain contour lines; According to the terrain features on the topographic and geological map, the terrain contour lines are extracted and saved in a vector format, and the elevation attribute is given; Step S12: constructing a digital elevation model; Through the terrain contour line vector file obtained in step S11, combined with the grid parameters given by the user, the grid processing is performed to generate a digital elevation model, and the digital elevation model is saved in a geotiff format; Step S13: extracting rock layer attributes and occurrence; According to the geological information given by the topographic and geological map, each independent geological unit is vectorized and given with the stratum attribute, and saved as a wp or shapefile surface file of mapgis; The occurrence position given in the topographic and geological map is vectorized and saved as a wt or shapefile point file of mapgis, and given with the corresponding occurrence and rock layer attribute; Step S14: extracting geological boundary lines; According to the geological information given by the topographic and geological map, the rock layer connection lines are vectorized and given with the boundary line attribute; Step S2: intelligent determination of anticline folds; The intelligent determination process of the anticline fold specifically comprises the following steps: Step S21: initialization of spatial elements; Modeling range: Bbox = {'west': west, 'east': east,'south': south, 'north': north}; Digital elevation model: Dem.tif; Geological unit: units = {u_1: {'name': name_1, 'era': era, 'points': points_1}, ……}; Geological boundary line: boundaries = {b_1: {'name': name_1, 'points': points_1, 'left': u_m, 'right': u_n}, ……}; Geological occurrence: occurrence = {u_1: [[x, y, dip, angle], ……], ……}; Step S22: key signs of anticline; Step S23: determination of core and two wings; According to the definition of anticline structure, the obtained geological boundary line can be used as the basic boundary element of the anticline fold, and the boundary line located at the first position of the set bdy3 is the outermost boundary line ly = [b01, b10] of the two wings of the anticline, and the geological boundary line located at the center of bdy3 is the core c = [b0n, bn0] of the anticline; Step S3: entity space subdivision and modeling; The working process of the entity space subdivision and modeling specifically comprises the following steps: Step S31: calculating the anticline constraint boundary; According to the geological boundary line ly of the two wings, the b01 and b10 boundary constraint stratum occurrence is read, and the coordinates of the K01, L01 and K10, L10 two points of the b01 and b10 geological boundary line are read respectively; According to the stratum occurrence of b01.left and b10.right, planes K01L01Q1P1 and K10L10Q1P1 are constructed, the plane K10L10Q1P1 is consistent with the tendency and dip angle of the stratum B on the left side of the geological boundary b01, and the plane K01L01Q1P1 is consistent with the tendency and dip angle of the stratum B on the right side of the geological boundary b02; Planes K01L01Q1P1 and K10L10Q1P1 intersect at a straight line P1Q1, coordinates Q1 (Q1x, Q1y, Q1z) and P1 (P1x, P1y, P1z) of points Q1P1 are calculated, and the higher value of the elevation of the Q1 point and the P1 point is taken as the model vertex; Planes K01L01Q1P1 and K10L10Q1P1 constitute the peripheral constraint boundary of the anticline; Step S32: stratum constraint interface construction; According to the method of step S31 for obtaining the peripheral constraint boundary of the anticline, K2, L2 and M2, N2 coordinates of the geological boundaries b02 and b20 are read, and the stratum occurrence on the right side of b02 and the left side of b20 is read; According to the tendency and dip angle, planes K02L02Q2P2 and K20L20Q2P2 are constructed, two planes intersect at a straight line P2Q2, and constraint interface control point coordinates Q2 (Q2x, Q2y, Q2z) and P2 (P2x, P2y, P2z) of the stratum C are calculated respectively; By using the same method, according to the parameters and spatial attributes of the remaining geological boundaries in bdy3, constraint interface control point coordinates Qn (Qnx, Qny, Qnz) and Pn (Pnx, Pny, Pnz) of the stratum N are completed; Step S33: cubic subdivision and denudation stratum three-dimensional reconstruction; According to the set minimum subdivision unit cube edge length, taking L01 of the geological boundary b01 as the starting point, searching in the horizontal direction and the vertical direction according to the digital elevation model dem.tif, the stratum unit units and the geological boundary bdy3, and taking the B stratum constraint boundary on the right side of b01 as the condition; If the current minimum unit center coordinates (x, y, z) are located in the space surrounded by the constraint interfaces K01L01Q1P1L10K10 and K02L02Q2P2L20K20, the attribute of the unit is assigned as the B stratum; Similarly, until all the cubic units belonging to the stratum N are searched, different colors are given to different cubic subdivision units according to different stratum attributes, that is, the spatial form of the denuded part of the anticline fold is directly presented. Step S4: denudation amount calculation; The working process of the denudation amount calculation is as follows: By searching the cubic subdivision units with different stratum attributes, the number thereof is counted, and the approximate volume of the restored denuded stratum is obtained according to the volume of the minimum cubic subdivision unit.

2. The method according to claim 1, wherein, The grid parameters are pixel ground resolution; The boundary attribute is the stratum attribute on both sides of the boundary.

3. The method according to claim 1, wherein, The analysis process of step S22 is as follows: Step S221: geological boundary trend consistency; The geological boundary results obtained in step S21 are traversed, the trend of all the geological boundaries is calculated, and the geological boundaries with consistent trends are screened. bdy1=[b1, b2, b3……, bn]; Step S222: Geological unit age-new-old-new arrangement uniqueness; Step S21 is traversed to obtain the consistent geological boundary bdy1, the ages of the geological units on both sides of the boundary are detected, it is determined that the arrangement of the geological units follows the change rule of new-old-new, the geological boundaries that do not meet the constraints are removed, and the result is obtained; bdy2=[b01, b02, b03……b0n]; Step S223: Symmetry of geological unit arrangement; Step S222 is traversed to obtain the result bdy2, whether the geological units on both sides of the geological boundary present spatial arrangement symmetry is detected, and the boundary elements of the structural anticline are obtained after removing the geological boundaries that do not meet the constraint conditions; bdy3=[b01, b02, b03……b0n, bn0……, b30, b20, b10].

4. The method according to claim 3, wherein, The rock layer on the right side of b01 and the rock layer on the left side of b10 are approximately axisymmetric and belong to the same rock layer unit, and b01 and b10 are the boundary constraints of the anticline fold modeling.

Images