A parameterized modeling method and device for dynamic symbols of fault structures

Abstract

This invention discloses a parametric modeling method and apparatus for dynamic symbols of fault structures. The method includes: 1) acquiring a user-defined set of stratigraphic parameters T, a set of fault plane parameters F, and a set of fault motion parameters M; 2) generating three-dimensional stratigraphic models of each layer based on T using an attitude inference method; 3) generating three-dimensional fault plane models based on F using an attitude inference method; 4) dividing each three-dimensional fault plane model into left and right halves using a three-dimensional trimming operation method and storing them in a three-dimensional stratigraphic entity model set S′; 5) performing fault slip simulation on each model in S′ using an affine transformation method based on M and storing it in a three-dimensional stratigraphic entity model set S″; 6) using vertex animation technology, generating a three-dimensional dynamic model of the fault symbol through animation interpolation, with S′ and S″ as the start and end frames respectively. This invention is more efficient and more flexible.

Description

Technical Field

[0001] This invention relates to geographic information technology, and more particularly to a parametric modeling method and apparatus for dynamic symbols of fault structures. Background Technology

[0002] Faults are planar structures formed when strata within the Earth undergo significant displacement along a fracture surface under the influence of internal forces. As a common geological structure, the symbolic representation of the dynamic evolution of faults is of great significance.

[0003] Currently, fault symbol representation is mostly limited to the two-dimensional level, which is insufficient to meet the needs of expressing the three-dimensional spatial structure and four-dimensional spatiotemporal evolution of faults in three-dimensional geographic information scenarios. Furthermore, while the three-dimensional models established for fault structures can be used as three-dimensional symbols, they are mostly generated using non-parametric modeling methods based on boreholes or attitude. This method is not only highly dependent on geological survey data but also has low modeling efficiency and limited applicability. Summary of the Invention

[0004] Purpose of the invention: This invention addresses the problems existing in the prior art by providing a more efficient and widely applicable parametric modeling method and apparatus for dynamic symbols of fault structures.

[0005] Technical solution: The parametric modeling method for dynamic symbols of fault structures described in this invention includes:

[0006] (1) Obtain the formation parameters, fault plane parameters and fault motion parameters set by the user, and form the formation parameter set T, the fault plane parameter set F and the fault motion parameter set M;

[0007] (2) Based on the stratigraphic parameter set T, generate three-dimensional stratigraphic models of each layer through the attitude inference method, and store them in the three-dimensional stratigraphic model set S;

[0008] (3) Based on the fault plane parameter set F, a three-dimensional fault plane model is generated by the attitude inference method and stored in the three-dimensional fault plane model set A;

[0009] (4) For each three-dimensional fault plane model in the three-dimensional fault plane model set A, the three-dimensional stratigraphic model in the three-dimensional stratigraphic model set S is divided into left and right blocks by the three-dimensional cutting operation method and stored in the three-dimensional stratigraphic entity model set S′.

[0010] (5) Based on the fault motion parameter set M, the fault sliding simulation is performed on all three-dimensional strata entity models in S′ according to the affine transformation method, and the simulated models are stored in the three-dimensional strata entity model set S″.

[0011] (6) Using vertex animation technology, the three-dimensional stratigraphic model sets S′ and S″ are used as the starting and ending frames respectively, and the three-dimensional dynamic model of the fault symbol is generated by animation interpolation.

[0012] Furthermore, step (1) includes:

[0013] (1-1) Read the formation parameters set by the user, including the reference point location, dip angle, strike, depth, length and width of each formation, and store them in the formation parameter set T;

[0014] (1-2) Read the fault plane parameters set by the user, including the reference point projection distance of each fault plane, the overall direction, the depth of each segment, and the dip angle of each segment, and store them in the fault plane parameter set F.

[0015] (1-3) Read the fault motion parameters set by the user, including the left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, right-side strike-slip component and motion simulation time of the two strata of each fault plane, and store them in the fault motion parameter set M.

[0016] Furthermore, step (2) includes:

[0017] (2-1) Select a stratum i from top to bottom, and read the length, width, dip angle, strike, and depth of stratum i from the stratum parameter set T, and denote them as l, d, and α, respectively. i β i ,dep i ;

[0018] (2-2) If the current stratum is the first stratum, then the coordinates of the reference point in the stratum parameter set T are used as the first feature point p of the three-dimensional stratum model. 0，1 coordinates (x) p，0，1 y p，0，1 , z p，0，1 ), and store it in the first stratigraphic model point set P0; if the current stratigraphy is another stratigraphy, then store it in the first feature point p of the previous stratigraphy. i-1，1 Using the coordinates as a reference, based on the thickness of the previous stratum (dep) i-1 The first feature point p of the current stratum is derived using the following formula. i，1 3D coordinates (x) p，i，1 y p，i，1 , z p，i，1 ), and store it in the current stratigraphic model point set P. i :

[0019]

[0020] In the formula, (x p，i-1，1 y p，i-1，1 , z p，i-1，1 ) represents the first characteristic point p of the upper stratum.i-1，1 The coordinates;

[0021] (2-3) Based on the first feature point p i，1 Length l, Inclination angle α i and towards β i The second characteristic point p of the current stratum is obtained by extrapolating along the dip direction using the following formula. i，2 coordinates (x) p，i，2 y p，i，2 , z p，i，2 ), and store it in the current stratigraphic model point set P. i :

[0022]

[0023] (2-4) Based on the first feature point p i，1 Second feature point p i，2 Based on width d and direction β i The third characteristic point p of the strata is obtained by extrapolating along the strike direction using the following formula. i，3 coordinates (x) p，i，3 y p，i，3 , z p，i，3 ) and the fourth feature point p i，4 coordinates (x) p，i，4 y p，i， 4, z p，i，4 ), and store it in the current stratigraphic model point set P. i :

[0024]

[0025]

[0026] (2-5) Repeat steps (2-1)-(2-4) until all strata have been traversed and all strata model point sets are obtained;

[0027] (2-6) Select any stratigraphic model point set P of stratigraphic i and its underlying stratigraphic strata. i and P i+1 Based on this, a three-dimensional stratigraphic solid model s of stratigraphy i is constructed. i And store it in the three-dimensional stratigraphic model set S;

[0028] (2-7) Repeat step (2-6) until all strata have been traversed, resulting in a three-dimensional stratigraphic model set S = {s i |i=1,2,...,SN}, where SN represents the number of strata.

[0029] Furthermore, step (3) includes:

[0030] (3-1) Select any fault plane j, and read the overall strike of the current fault plane and the distance to the reference point from the fault plane parameter set F, denoted as δ respectively. j f d，j ;

[0031] (3-2) Select a segment k of the current fault plane from top to bottom, and read the depth and dip angle of the current segment from the fault plane parameter set F, denoted as h respectively. j，k and γ j，k ;

[0032] (3-3) If the current segment is the first segment, then the first feature point v of the first segment of the fault plane model is derived along the dip direction of the strata according to the following formula. j，0，1 coordinates (x) v，j，0，1 y v，j，0，1 , z v，j，0，1 ), and store it in the first fault plane point set V0:

[0033]

[0034] In the formula, (x p，0，1 y p，0，1 , z p，0，1 ) is the first feature point p in the first layer of the three-dimensional stratigraphic model. 0，1 The coordinates are β0, which represents the strike of the first stratum, and α0, which represents the dip angle of the first stratum.

[0035] If the current segment is not the first segment, then the third feature point v of the previous segment will be... j，k-1，3 The coordinates are assigned to the first feature point v. j，k，1 The coordinates are stored in the fault plane point set V of the current segment. j，k ;

[0036] (3-4) If the current segment is the first segment, then based on the stratum width d, follow the following formula along the fault strike δ in the fault plane parameters. j The direction was deduced to obtain the second characteristic point v of the first segment of the fault plane. j，0，2 coordinates (x) v，j，0，2 y v，j，0，2 , z v，j，0，2 ), and store it in the fault plane point set V of the first segment. j，0 :

[0037]

[0038] If the current segment is not the first segment, then the fourth feature point v of the previous segment will be... j，k-1，4 The coordinates are assigned to the second feature point v. j，k，2 The coordinates are stored in the fault plane point set V of the current segment. j，k

[0039] (3-5) Based on the first feature point v of the current segment j，k，1 Second feature point v j，k，2 The third characteristic point v of the current segment is obtained by deducing along the dip direction of the fault plane using the following formula. j，k，3 and the fourth feature point v j，k，4 And store it in the fault plane point set V of the current segment. j，k :

[0040]

[0041]

[0042] In the formula, (x v，j，0，3 y v，j，0，3 , z v，j，0，3 ) represents v j，k，3 The coordinates, (x v，j，0，4 y v，j，0，4 , z v，j，0，4 ) represents v j，k，4 The coordinates;

[0043] (3-6) Repeat steps (3-2)-(3-5) until all segments of the current fault plane have been traversed, obtaining the point set V of the current fault plane. j ={V j，k};

[0044] (3-7) Based on the fault plane point set V j The three-dimensional fracture surface model a is obtained by suturing the facets. j And store it in the three-dimensional fault plane model set A;

[0045] (3-8) Repeat steps (3-1)-(3-7) until all fault planes have been traversed, obtaining the three-dimensional fault plane model set A = {a j |j=1,2,...,AN}, where AN represents the number of fault planes.

[0046] Furthermore, step (4) includes:

[0047] (4-1) Copy all three-dimensional stratigraphic models in the three-dimensional stratigraphic model set S to a new three-dimensional stratigraphic entity model set S′;

[0048] (4-2) Select any fault plane model a from the set of three-dimensional fault plane models A. j ;

[0049] (4-3) Select any model s′ from the set of three-dimensional stratigraphic entity models S′ n ;

[0050] (4-4) Based on the fault plane model a jThrough three-dimensional clipping operations, the three-dimensional stratigraphic model s′ n The three-dimensional stratigraphic model is divided into left and right blocks s′ n，j，l and s′ n，j，r The two 3D stratigraphic models are independently encoded and stored in a 3D stratigraphic entity model set S′, and the 3D stratigraphic model s′ is deleted from the set. n ;

[0051] (4-5) Repeat steps (4-3)-(4-4) until all three-dimensional stratigraphic models are fault plane a. j Segmentation complete;

[0052] (4-6) Repeat steps (4-2)-(4-5) until all fault planes have been traversed, and obtain the final updated three-dimensional stratigraphic entity model set S′.

[0053] Furthermore, step (5) includes:

[0054] (5-1) Select any fault plane model a from the set of three-dimensional fault plane models A. j ;

[0055] (5-2) Extracting a from the fault motion parameter set M j The left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, and right-side strike-slip component of the corresponding fault plane j are denoted as ld, respectively. j ls j rd j and rs j ;

[0056] (5-3) Select any element s′ from the set of three-dimensional stratigraphic entity models S′ n And determine s′ based on independent encoding. n Is it part of the left disk obtained after dividing the cross plane j? If so, proceed to step (5-4); otherwise, proceed to step (5-5).

[0057] (5-4) Traverse s′ n All points were used, and new three-dimensional formation solid models s″ were calculated based on the dip and strike-slip components of the left side of the footing. n And store it in the three-dimensional stratigraphic entity model set S″:

[0058] (5-5) Traverse s′ n All points were used, and new three-dimensional formation solid models s″ were calculated based on the right-side dip and strike-slip components to obtain the dip and strike-slip simulation results. n And store it in the three-dimensional stratigraphic entity model set S″:

[0059] (5-6) Repeat steps (5-1)-(5-5) until all fault planes have been traversed and the motion simulation of all fault planes cutting the strata has been completed.

[0060] Furthermore, in step (5-4), the three-dimensional stratigraphic entity model s″ n The simulation method is as follows:

[0061] Perform a tilting simulation using the following formula:

[0062]

[0063] In the formula, (x, y, z) represents s′ n The three-dimensional coordinates of any point in s, (x′, y′, z′) represent s′. n The three-dimensional coordinates of any point after the tilt-slip simulation, γ j，0 Represents s′ n The dip angle corresponding to the first segment of the fault, δ j Represents s′ n The overall strike of the corresponding fault;

[0064] The following formula is used to simulate the skid:

[0065]

[0066] Where (x″, y″, z″) represents s′ n The three-dimensional formation solid model s″ obtained after completing the dip-slip and strike-slip simulation transformation n The three-dimensional coordinates of any point.

[0067] Furthermore, in step (5-5), the three-dimensional stratigraphic solid model s″ n The simulation method is as follows:

[0068] Perform a tilting simulation using the following formula:

[0069]

[0070] In the formula, (x, y, z) represents s′ n The three-dimensional coordinates of any point in s, (x′, y′, z′) represent s′. n The three-dimensional coordinates of any point after the tilt-slip simulation, γ j，0 Represents s′ n The dip angle corresponding to the first segment of the fault, δ j Represents s′ n The overall strike of the corresponding fault;

[0071] The following formula is used to simulate the skid:

[0072]

[0073] Where (x″, y″, z″) represents s′ n The three-dimensional formation solid model s″ obtained after completing the dip-slip and strike-slip simulation transformation n The three-dimensional coordinates of any point.

[0074] Furthermore, step (6) includes:

[0075] (6-1) Select any element s′ from the three-dimensional stratigraphic entity set S′. n Select the corresponding element s″ from the three-dimensional stratigraphic entity set S″. n ;

[0076] (6-2) Select s′ n any point v′ in n，m and s″ n The corresponding point v″ in n，m ;

[0077] (6-3) Using point v″ n，m Using the three-dimensional coordinates as a reference, set point v′ according to the following formula. n，m The deformation channel is (dx, dy, dz), and the motion simulation time m in the fault motion parameter set M is used. time Assign the value to point v′ n，m Deformation time:

[0078]

[0079] In the formula, (x′ n，m y′ n，m , z′ n，m ) represents point v′ n，m The three-dimensional spatial coordinates, (x″) n，m ,y″ n，m , z″ n，m ) represents point v″ n，n Three-dimensional spatial coordinates;

[0080] (6-4) Repeat steps (6-2)-(6-3) until s′ n All points are traversed to obtain the deformation information of all points;

[0081] (6-5) Repeat steps (6-1)-(6-4) until all elements of S′ have been traversed, thus completing the dynamic information generation of all elements;

[0082] (6-6) Export the three-dimensional stratigraphic solid model of S′ and the fault plane model of A, output them as FBX format model files, and bind the corresponding materials.

[0083] The parametric modeling method for dynamic symbols of fault structures according to the present invention includes 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.

[0084] Beneficial Effects: Compared with existing technologies, the significant advantages of this invention are: based on various fault parameters set by the user, this invention generates stratigraphic and fault plane models, and then segments the stratigraphic model through three-dimensional clipping operations. It then performs fault motion simulation and dynamic information representation, achieving parametric modeling of three-dimensional fault dynamic symbols. This invention can construct different types of symbolic fault models based on parameter adjustments, ensuring modeling efficiency while possessing high accuracy and a flexible range of applications, thus having significant research value and application potential. Attached Figure Description

[0085] Figure 1 This is a flowchart of the parametric modeling method for dynamic symbols of fault structures provided by the present invention;

[0086] Figure 2 This is a schematic diagram illustrating the deduction process of various feature points on the ground stratum in the embodiment;

[0087] Figure 3 This is a schematic diagram showing the top feature points of the first segment of the fault plane in this embodiment;

[0088] Figure 4 This is a schematic diagram of the segmented deduction of the fault plane in this embodiment;

[0089] Figure 5 This is a top view of the stratigraphic segmentation and motion simulation process in this embodiment; wherein (a) fault 1_1 segments the stratigraphic model, (b) after fault 1_1 segments the stratigraphic model, motion simulation is performed on the stratigraphic model, (c) fault 2_1 cuts the stratigraphic model after segmentation by fault 1_1, and (d) after fault 2_1 segments the stratigraphic model, motion simulation is performed on the stratigraphic model.

[0090] Figure 6 The results of dynamic fault symbol modeling in this embodiment are as follows: (a) the three-dimensional fault symbol model of the starting frame and (b) the three-dimensional fault symbol model of the ending frame.

[0091] Figure 7 The fault modeling results are divided according to the relative movement of the two sides of the fault; where (a) is a normal fault, (b) is a reverse fault, and (c) is a strike-slip fault / strike-slip fault.

[0092] Figure 8 The fault modeling results are based on the relationship between the fault and the cutting strata; where (a) strike faults / longitudinal faults, (b) dip faults, and (c) oblique faults.

[0093] Figure 9These are modeling results for several typical composite faults; (a) graben, (b) horst, (c) imbricate fault, and (d) step fault. Detailed Implementation

[0094] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0095] Example 1

[0096] This embodiment provides a parametric modeling method for dynamic symbols of fault structures, such as... Figure 1 As shown, it includes the following steps:

[0097] (1) Obtain the formation parameters, fault plane parameters and fault motion parameters set by the user, and form the formation parameter set T, the fault plane parameter set F and the fault motion parameter set M.

[0098] The method for forming the set is as follows: (1-1) Read the stratigraphic parameters set by the user, including the reference point location, dip angle, strike, depth, length and width of each stratum, and store them in the stratigraphic parameter set T; (1-2) Read the fault plane parameters set by the user, including the reference point projection distance, overall strike, depth of each segment, dip angle of each segment, and store them in the fault plane parameter set F; (1-3) Read the fault motion parameters set by the user, including the left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, right-side strike-slip component and motion simulation time of the two sides of each fault plane, and store them in the fault motion parameter set M.

[0099] As shown in Table 1, the relevant parameters of the composite fracture composed of two strike faults were selected as experimental data in this embodiment.

[0100] Table 1 User Settings Parameter Information Table

[0101]

[0102]

[0103] In this embodiment, the reference point is specifically the left vertex.

[0104] (2) Based on the set of stratigraphic parameters T, three-dimensional stratigraphic models of each layer are generated by the attitude inference method and stored in the set of three-dimensional stratigraphic models S.

[0105] This step specifically includes:

[0106] (2-1) Select a stratum i from top to bottom, and read the length, width, dip angle, strike, and depth of stratum i from the stratum parameter set T, and denote them as l, d, and α, respectively. i β i ,dep i In this embodiment, l is 100m, d is 50m, and the inclination angle α of each layer is... i All are 0°, and the orientation β of each layer i All are 0°, and the depth of each layer is dep. i Both are 10m;

[0107] (2-2) If the current stratum is the first stratum, then the coordinates of the reference point in the stratum parameter set T are used as the first feature point p of the three-dimensional stratum model. 0，1 coordinates (x) p，0，1 y p，0，1 , z p，0，1 ), and store it in the first stratigraphic model point set P0; if the current stratigraphy is another stratigraphy, then store it in the first feature point p of the previous stratigraphy. i-1，1 Using the coordinates as a reference, based on the thickness of the previous stratum (dep) i-1 The first feature point p of the current stratum is derived using the following formula. i，1 3D coordinates (x) p，i，1 y p，i，1 , z p，i，1 ), and store it in the current stratigraphic model point set P. i :

[0108]

[0109] In the formula, (x p，i-1，1 y p，i-1，1 , z p，i-1，1 ) represents the first characteristic point p of the upper stratum. i-1，1 The coordinates;

[0110] (2-3) Based on the first feature point p i，1 Length l, Inclination angle α i and towards β i The second characteristic point p of the current stratum is obtained by extrapolating along the dip direction using the following formula. i，2 coordinates (x) p，i，2 y p，i，2 , z p，i，2 ), and store it in the current stratigraphic model point set P. i :

[0111]

[0112] (2-4) Based on the first feature point p i，1 Second feature point p i，2Based on width d and direction β i The third characteristic point p of the strata is obtained by extrapolating along the strike direction using the following formula. i，3 coordinates (x) p，i，3 y p，i，3 , z p，i，3 ) and the fourth feature point p i，4 coordinates (x) p，i，4 y p，i，4 , z p，i，4 ), and store it in the current stratigraphic model point set P. i In this embodiment, four feature points p on the ground plane i，1 p i，2 p i，3 and p i，4 The generation of such as Figure 2 As shown;

[0113]

[0114]

[0115] (2-5) Repeat steps (2-1)-(2-4) until all strata have been traversed and all strata model point sets are obtained;

[0116] (2-6) Select any stratigraphic model point set P of stratigraphic i and its underlying stratigraphic strata. i and P i+1 Based on this, a three-dimensional stratigraphic solid model s of stratigraphy i is constructed. i And store it in the three-dimensional stratigraphic model set S;

[0117] (2-7) Repeat step (2-6) until all strata have been traversed, resulting in a three-dimensional stratigraphic model set S = {s i |i=1,2,...,SN}, where SN represents the number of strata. In this embodiment, SN is 3.

[0118] (3) Based on the fault plane parameter set F, a three-dimensional fault plane model is generated by the attitude inference method and stored in the three-dimensional fault plane model set A.

[0119] This step includes:

[0120] (3-1) Select any fault plane j, and read the overall strike of the current fault plane and the distance to the reference point from the fault plane parameter set F, denoted as δ respectively. j f d，j In this embodiment, there are two fault planes, δ0 and δ1, which are 30° and -60° respectively.

[0121] (3-2) Select a segment k of the current fault plane from top to bottom, and read the depth and dip angle of the current segment from the fault plane parameter set F, denoted as h respectively. j，k and γ j，k In this embodiment, the depth h of any segment of the two faults j，k Both faults are 10m deep, and the dip angle γ of any segment of the two faults is... 0，k and γ 1，k They are -60° and -80° respectively;

[0122] (3-3) If the current segment is the first segment, then the first feature point v of the first segment of the fault plane model is derived along the dip direction of the strata according to the following formula. j，0，1 coordinates (x) v，j，0，1 y v，j，0，1 , z v，j，0，1 ), and store it in the first fault plane point set V0:

[0123]

[0124] In the formula, (x p，0，1 y p，0，1 , z p，0，1 ) is the first feature point p in the first layer of the three-dimensional stratigraphic model. 0，1 The coordinates are given, where β0 represents the strike of the first stratum and α0 represents the dip angle of the first stratum; in this embodiment, the distance f between the reference points of the two faults is calculated. d，0 and f d，1 The heights are 20m and 80m respectively, with the left vertex p of the first floor. 0，1 The coordinates are (0, 0, 0), β0 is 0°, and α0 is 0°;

[0125] If the current segment is not the first segment, then the third feature point v of the previous segment will be... j，k-1，3 The coordinates are assigned to the first feature point v. j，k，1 The coordinates are stored in the fault plane point set V of the current segment. j，k ;

[0126] (3-4) If the current segment is the first segment, then based on the stratum width d, follow the following formula along the fault strike δ in the fault plane parameters. j The direction was deduced to obtain the second characteristic point v of the first segment of the fault plane. j，0，2 coordinates (x) v，j，0，2 y v，j，0，2 , z v，j，0，2 ), and store it in the fault plane point set V of the first segment. j，0 :

[0127]

[0128] If the current segment is not the first segment, then the fourth feature point v of the previous segment will be...j，k-1，4 The coordinates are assigned to the second feature point v. j，k，2 The coordinates are stored in the fault plane point set V of the current segment. j，k In the embodiment, the overall strikes of the two faults, δ0 and δ1, are 30° and -60°, respectively; the characteristic point v of the first segment of fault 0. 0，0，1 and v 0，0，2 The deduction is as follows Figure 3 As shown;

[0129] (3-5) Based on the first feature point v of the current segment j，k，1 Second feature point v j，k，2 The third characteristic point v of the current segment is obtained by deducing along the dip direction of the fault plane using the following formula. j，k，3 and the fourth feature point v j，k，4 And store it in the fault plane point set V of the current segment. j，k :

[0130]

[0131]

[0132] In the formula, (x v，j，0，3 y v，j，0，3 , z v，j，0，3 ) represents v j，k，3 The coordinates, (x v，j，0，4 y v，j，0，4 , z v，j，0，4 ) represents v j，k，4 The coordinates;

[0133] (3-6) Repeat steps (3-2)-(3-5) until all segments of the current fault plane have been traversed, obtaining the point set V of the current fault plane. j ={V j，k};like Figure 4 As shown;

[0134] (3-7) Based on the fault plane point set V j The three-dimensional fracture surface model a is obtained by suturing the facets. j And store it in the three-dimensional fault plane model set A;

[0135] (3-8) Repeat steps (3-1)-(3-7) until all fault planes have been traversed, obtaining the three-dimensional fault plane model set A = {a j |j=1,2,...,AN}, where AN represents the number of fault planes. In this embodiment, AN is 2.

[0136] (4) For each three-dimensional fault plane model in the three-dimensional fault plane model set A, the three-dimensional stratigraphic model in the three-dimensional stratigraphic model set S is divided into left and right plates by the three-dimensional cutting operation method and stored in the three-dimensional stratigraphic entity model set S′.

[0137] This step includes:

[0138] (4-1) Copy all three-dimensional stratigraphic models in the three-dimensional stratigraphic model set S to a new three-dimensional stratigraphic entity model set S′;

[0139] (4-2) Select any fault plane model a from the set of three-dimensional fault plane models A. j ;

[0140] (4-3) Select any model s′ from the set of three-dimensional stratigraphic entity models S′ n ;

[0141] (4-4) Based on the fault plane model a j Through three-dimensional clipping operations, the three-dimensional stratigraphic model s′ n The three-dimensional stratigraphic model is divided into left and right blocks s′ n，j，l and s′ n，j，r Encode the two 3D stratigraphic models independently (any encoding method that can distinguish the left and right sides is acceptable), store them in the 3D stratigraphic entity model set S′, and delete the 3D stratigraphic model s′ from it. n The number of elements in the three-dimensional stratigraphic entity model set S′ is SN′≤SN×AN×2, and it will be continuously updated with the iteration of steps (4-3)-(4-4);

[0142] (4-5) Repeat steps (4-3)-(4-4) until all three-dimensional stratigraphic models are fault plane a. j Segmentation complete;

[0143] (4-6) Repeat steps (4-2)-(4-5) until all fault planes have been traversed, resulting in the final updated three-dimensional stratigraphic entity model set S′. In this embodiment, the final SN′ is 12.

[0144] (5) Based on the fault motion parameter set M, the fault sliding simulation is performed on all three-dimensional strata entity models in S′ according to the affine transformation method, and the simulated models are stored in the three-dimensional strata entity model set S″.

[0145] This step includes:

[0146] (5-1) Select any fault plane model a from the set of three-dimensional fault plane models A. j ;

[0147] (5-2) Extracting a from the fault motion parameter set M j The left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, and right-side strike-slip component of the corresponding fault plane j are denoted as ld, respectively. j ls j rd j and rs j In this embodiment, ld0 and ld1 are both 3m, ls0 and ls1 are both -5m, rd0 and rd1 are both -3m, and rs0 and rs1 are both 5m.

[0148] (5-3) Select any element s′ from the set of three-dimensional stratigraphic entity models S′ n And determine s′ based on independent encoding. n Is it part of the left disk obtained after dividing the cross plane j? If so, proceed to step (5-4); otherwise, proceed to step (5-5).

[0149] (5-4) Traverse s′ n All points are used, and based on the left-side dip and strike-slip components, a new three-dimensional formation solid model s″ after dip and strike-slip simulations is obtained according to the following formula. n And store it in the three-dimensional stratigraphic entity model set S″:

[0150]

[0151]

[0152] In the formula, (x, y, z) represents s′ n The three-dimensional coordinates of any point in s, (x′, y′, z′) represent s′. n The three-dimensional coordinates of any point after the tilt-slip simulation, γ j，0 Represents s′ n The dip angle corresponding to the first segment of the fault, δ j Represents s′ n The overall strike of the corresponding fault; (x″, y″, z″) represents s′ n The three-dimensional formation solid model s″ obtained after completing the dip-slip and strike-slip simulation transformation n The three-dimensional coordinates of any point.

[0153] (5-5) Traverse s′ n All points are used, and based on the right-side dip and strike-slip components, a new three-dimensional formation solid model s″ after dip and strike-slip simulations is obtained according to the following formula. n And store it in the three-dimensional stratigraphic entity model set S″:

[0154]

[0155]

[0156] (5-6) Repeat steps (5-1)-(5-5) until all fault planes have been traversed and the motion simulation of all fault planes cutting the strata is complete. The strata segmentation and motion simulation process in this embodiment is as follows: Figure 5 As shown.

[0157] (6) Using vertex animation technology, the three-dimensional stratigraphic model sets S′ and S″ are used as the starting and ending frames respectively, and the three-dimensional dynamic model of the fault symbol is generated by animation interpolation.

[0158] This step includes:

[0159] (6-1) Select any element s′ from the three-dimensional stratigraphic entity set S′. n Select the corresponding element s″ from the three-dimensional stratigraphic entity set S″. n ;

[0160] (6-2) Select s′ n any point v′ in n，m and s″ n The corresponding point v″ in n，m ;

[0161] (6-3) Using point v″ n，m Using the three-dimensional coordinates as a reference, set point v′ according to the following formula. n，m The deformation channel is (dx, dy, dz), and the motion simulation time m in the fault motion parameter set M is used. time Assign the value to point v′ n，m Deformation time:

[0162]

[0163] In the formula, (x′ n，m y′ n，m , z′ n，m ) represents point v′ n，m The three-dimensional spatial coordinates, (x″) n，m ,y″ n，m , z″ n，m ) represents point v″ n，n Three-dimensional spatial coordinates;

[0164] (6-4) Repeat steps (6-2)-(6-3) until s′ n All points are traversed to obtain the deformation information of all points;

[0165] (6-5) Repeat steps (6-1)-(6-4) until all elements of S′ have been traversed, thus completing the dynamic information generation of all elements;

[0166] (6-6) Export the 3D stratigraphic solid model of S′ and the fault plane model of A, outputting them as FBX format model files, and bind the corresponding materials. For example... Figure 6 As shown.

[0167] This embodiment uses a composite strike fault as an example to construct a three-dimensional fault structure model. Different types of three-dimensional fault models can also be constructed and generated based on different fault plane parameters, such as... Figure 7 , Figure 8 , Figure 9 As shown. In this embodiment, the stratigraphic model segmentation is performed solely based on the clipping operation interface provided by the CGAL open-source code. Other libraries can also be used to complete the 3D clipping operation. This embodiment only exports the fold structure 3D model in FBX format; other formats such as OBJ can also be exported as 3D geological models.

[0168] Example 2

[0169] This embodiment provides a parametric modeling apparatus for dynamic symbols of fault structures, providing services for the implementation of the method in Embodiment 1 of the present invention. It is manifested in the form of a general-purpose computing device, and its components may include, but are not limited to: one or more processors or processing units, system memory, and a bus connecting different system components, including system memory and processing units. Typically, it includes various computer system readable media. These media can be any available media accessible by the device, including volatile and non-volatile media, removable and non-removable media.

[0170] System memory may include computer system readable media in the form of volatile memory, and the apparatus may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, the storage system may be used to read and write non-removable, non-volatile magnetic media. A program / utility having a set (at least one) of program modules may be stored, for example, in memory. Such program modules include, but are not limited to, an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include an implementation of a network environment. Program modules typically perform the functions and / or methods described in the embodiments of the present invention. The processing unit performs various functional applications and data processing, such as implementing the methods provided in Embodiment 1 of the present invention, by running programs stored in system memory.

Claims

1. A parametric modeling method for dynamic symbols of fault structures, characterized in that... The method includes: (1) Obtain the formation parameters, fault plane parameters and fault motion parameters set by the user, and form the formation parameter set T, the fault plane parameter set F and the fault motion parameter set M; (2) Based on the stratigraphic parameter set T, generate three-dimensional stratigraphic models of each layer through the attitude inference method, and store them in the three-dimensional stratigraphic model set S; (3) Based on the fault plane parameter set F, a three-dimensional fault plane model is generated by the attitude inference method and stored in the three-dimensional fault plane model set A; (4) For each three-dimensional fault plane model in the three-dimensional fault plane model set A, the three-dimensional stratigraphic model in the three-dimensional stratigraphic model set S is divided into left and right blocks by the three-dimensional cutting operation method and stored in the three-dimensional stratigraphic entity model set S′. (5) Based on the fault motion parameter set M, the fault sliding simulation is performed on all three-dimensional strata entity models in S′ according to the affine transformation method, and the simulated models are stored in the three-dimensional strata entity model set S″. (6) Using vertex animation technology, with S′ and S″ as the start and end frames respectively, a three-dimensional dynamic model of the fault symbol is generated by animation interpolation.

2. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (1) includes: (1-1) Read the formation parameters set by the user, including the reference point location, dip angle, strike, depth, length and width of each formation, and store them in the formation parameter set T; (1-2) Read the fault plane parameters set by the user, including the reference point projection distance of each fault plane, the overall strike, the depth of each segment, and the dip angle of each segment, and store them in the fault plane parameter set F. (1-3) Read the fault motion parameters set by the user, including the left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, right-side strike-slip component and motion simulation time of the two strata of each fault plane, and store them in the fault motion parameter set M.

3. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (2) includes: (2-1) Select a stratum i from top to bottom, and read the length, width, dip angle, strike, and depth of stratum i from the stratum parameter set T, and denot them as l, d, and α, respectively. i β i ,dep i ; (2-2) If the current stratum is the first stratum, then the coordinates of the reference point in the stratum parameter set T are used as the first feature point p of the three-dimensional stratum model. 0，1 coordinates (x) p，0，1 y p，0，1 , z p，0，1 ), and store it in the first stratigraphic model point set P0; if the current stratigraphy is another stratigraphy, then store it in the first feature point p of the previous stratigraphy. i-1，1 Using the coordinates as a reference, based on the thickness of the previous stratum (dep) i-1 The first feature point p of the current stratum is derived using the following formula. i，1 3D coordinates (x) p，i，1 y p，i，1 , z p，i，1 ), and store it in the current stratigraphic model point set P. i : In the formula, (x p，i-1，1 y p，i-1，1 , z p，i-1，1 ) represents the first characteristic point p of the upper stratum. i-1，1 The coordinates; (2-3) Based on the first feature point p i，1 Length l, Inclination angle α i and towards β i The second characteristic point p of the current stratum is obtained by extrapolating along the dip direction using the following formula. i，2 coordinates (x) p，i，2 y p，i，2 , z p，i，2 ), and store it in the current stratigraphic model point set P. i : (2-4) Based on the first feature point p i，1 Second feature point p i，2 Based on width d and direction β i The third characteristic point p of the strata is obtained by extrapolating along the strike direction using the following formula. i，3 coordinates (x) p，i，3 y p，i，3 , z p，i，3 ) and the fourth feature point p i，4 coordinates (x) p，i，4 y p，i，4 , z p，i，4 ), and store it in the current stratigraphic model point set P. i : (2-5) Repeat steps (2-1)-(2-4) until all strata have been traversed and all strata model point sets are obtained; (2-6) Select any stratigraphic model point set P of stratigraphic i and its underlying stratigraphic strata. i and P i+1 Based on this, a three-dimensional stratigraphic solid model s of stratigraphy i is constructed. i And store it in the three-dimensional stratigraphic model set S; (2-7) Repeat step (2-6) until all strata have been traversed, resulting in a three-dimensional stratigraphic model set S = {s i |i=1,2,...,SN),where SN represents the number of strata.

4. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (3) includes: (3-1) Select any fault plane j, and read the overall strike of the current fault plane and the distance to the reference point from the fault plane parameter set F, denoted as δ respectively. j f d，j ; (3-2) Select a segment k of the current fault plane from top to bottom, and read the depth and dip angle of the current segment from the fault plane parameter set F, denoted as h respectively. j，k and γ j，k ; (3-3) If the current segment is the first segment, then the first feature point v of the first segment of the fault plane model is derived along the dip direction of the strata according to the following formula. j，0，1 coordinates (x) v，j，0，1 y v，j，0，1 , z v，j，0，1 ), and store it in the first fault plane point set V0: In the formula, (x p，0，1 y p，0，1 , z p，0，1 ) is the first feature point p in the first layer of the three-dimensional stratigraphic model. 0，1 The coordinates are β0, which represents the strike of the first stratum, and α0, which represents the dip angle of the first stratum. If the current segment is not the first segment, then the third feature point v of the previous segment will be... j，k-1，3 The coordinates are assigned to the first feature point v. j，k，1 The coordinates are stored in the fault plane point set V of the current segment. j，k ; (3-4) If the current segment is the first segment, then based on the stratum width d, follow the following formula along the fault strike δ in the fault plane parameters. j The direction was deduced to obtain the second characteristic point v of the first segment of the fault plane. j，0，2 coordinates (x) v，j，0，2 y v，j，0，2 , z v，j，0，2 ), and store it in the fault plane point set V of the first segment. j，0 : If the current segment is not the first segment, then the fourth feature point v of the previous segment will be... j，k-1，4 The coordinates are assigned to the second feature point v. j，k，2 The coordinates are stored in the fault plane point set V of the current segment. j，k (3-5) Based on the first feature point v of the current segment j，k，1 Second feature point v j，k，2 The third characteristic point v of the current segment is obtained by deducing along the dip direction of the fault plane using the following formula. j，k，3 and the fourth feature point v j，k，4 And store it in the fault plane point set V of the current segment. j，k : In the formula, (x v，j，0，3 y v，j，0，3 , z v，j，0，3 ) represents v j，k，3 The coordinates, (x v，j，0，4 y v，j，0，4 , z v，j，0，4 ) represents v j，k，4 The coordinates; (3-6) Repeat steps (3-2)-(3-5) until all segments of the current fault plane have been traversed, obtaining the point set V of the current fault plane. j ={V j，k }; (3-7) Based on the fault plane point set V j The three-dimensional fracture surface model a is obtained by stitching the facets. j And store it in the three-dimensional fault plane model set A; (3-8) Repeat steps (3-1)-(3-7) until all fault planes have been traversed, obtaining the three-dimensional fault plane model set A = {a j |j=1,2,...,AN}, where AN represents the number of fault planes.

5. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (4) includes: (4-1) Copy all three-dimensional stratigraphic models in the three-dimensional stratigraphic model set S to a new three-dimensional stratigraphic entity model set S′; (4-2) Select any fault plane model a from the set of three-dimensional fault plane models A. j ; (4-3) Select any model s′ from the set of three-dimensional stratigraphic entity models S′ n ; (4-4) Based on the fault plane model a j Through three-dimensional clipping operations, the three-dimensional stratigraphic model s′ n The three-dimensional stratigraphic model is divided into left and right blocks s′ n，j，l and s′ n，j，r The two 3D stratigraphic models are independently encoded and stored in a 3D stratigraphic entity model set S′, and the 3D stratigraphic model s′ is deleted from the set. n ; (4-5) Repeat steps (4-3)-(4-4) until all three-dimensional stratigraphic models are fault plane a. j Segmentation complete; (4-6) Repeat steps (4-2)-(4-5) until all fault planes have been traversed, and obtain the final updated three-dimensional stratigraphic entity model set S′.

6. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (5) includes: (5-1) Select any fault plane model a from the set of three-dimensional fault plane models A. j ; (5-2) Extracting a from the fault motion parameter set M j The left-side dip-slip component, left-side strike-slip component, right-side dip-slip component, and right-side strike-slip component of the corresponding fault plane j are denoted as ld, respectively. j ls j rd j and rs j ; (5-3) Select any element s′ from the set of three-dimensional stratigraphic entity models S′ n And determine s′ based on independent encoding. n Is it part of the left disk obtained after dividing the cross plane j? If so, proceed to step (5-4); otherwise, proceed to step (5-5). (5-4) Traverse s′ n All points were used, and new three-dimensional formation solid models s″ were calculated based on the dip and strike-slip components of the left side of the footing. n And store it in the three-dimensional stratigraphic entity model set S″: (5-5) Traverse s′ n All points were used, and new three-dimensional formation solid models s″ were calculated based on the right-side dip and strike-slip components to obtain the dip and strike-slip simulation results. n And store it in the three-dimensional stratigraphic entity model set S″: (5-6) Repeat steps (5-1)-(5-5) until all fault planes have been traversed and the motion simulation of all fault planes cutting the strata has been completed.

7. The parametric modeling method for dynamic symbols of fault structures according to claim 6, characterized in that: In step (5-4), the three-dimensional stratigraphic solid model s″ n The simulation method is as follows: Perform a tilting simulation using the following formula: In the formula, (x, y, z) represents s′ n The three-dimensional coordinates of any point in s, (x′, y′, z′) represent s′. n The three-dimensional coordinates of any point after the tilt-slip simulation, γ j，0 Represents s′ n The dip angle corresponding to the first segment of the fault, δ j Represents s′ n The overall strike of the corresponding fault; The following formula is used to simulate the skid: Where (x″, y″, z″) represents s′ n The three-dimensional formation solid model s″ obtained after completing the dip-slip and strike-slip simulation transformation n The three-dimensional coordinates of any point.

8. The parametric modeling method for dynamic symbols of fault structures according to claim 6, characterized in that: In step (5-5), the three-dimensional stratigraphic solid model s″ n The simulation method is as follows: Perform a tilting simulation using the following formula: In the formula, (x, y, z) represents s′ n The three-dimensional coordinates of any point in s, (x′, y′, z′) represent s′. n The three-dimensional coordinates of any point after the tilt-slip simulation, γ j，0 Represents s′ n The dip angle corresponding to the first segment of the fault, δ j Represents s′ n The overall strike of the corresponding fault; The following formula is used to simulate the skid: Where (x″, y″, z″) represents s′ n The three-dimensional formation solid model s″ obtained after completing the dip-slip and strike-slip simulation transformation n The three-dimensional coordinates of any point.

9. The parametric modeling method for dynamic symbols of fault structures according to claim 1, characterized in that: Step (6) includes: (6-1) Select any element s′ from the three-dimensional stratigraphic entity set S′. n Select the corresponding element s″ from the three-dimensional stratigraphic entity set S″. n ; (6-2) Select s′ n any point v′ in n，m and s″ n The corresponding point v″ in n，m ; (6-3) Using point v″ n，m Using the three-dimensional coordinates as a reference, set point v′ according to the following formula. n，m The deformation channel is (dx, dy, dz), and the motion simulation time m in the fault motion parameter set M is used. time Assign the value to point v′ n，m Deformation time: In the formula, (x′ n，m y′ n，m , z′ n，m ) represents point v′ n，m The three-dimensional spatial coordinates, (x″) n，m ,y″ n，m , z″ n，m ) represents point v″ n，m Three-dimensional spatial coordinates; (6-4) Repeat steps (6-2)-(6-3) until s′ n All points are traversed to obtain the deformation information of all points; (6-5) Repeat steps (6-1)-(6-4) until all elements of S′ have been traversed, thus completing the dynamic information generation of all elements; (6-6) Export the three-dimensional stratigraphic solid model of S′ and the fault plane model of A, output them as FBX format model files, and bind the corresponding materials.

10. A parametric modeling method for dynamic symbols of fault structures, comprising a processor and a computer program stored in memory and executable on the processor, characterized in that: When the processor executes the program, it implements the method described in any one of claims 1-9.

Images