A method for iterative modeling of three-dimensional structural models using massive focal mechanism solution data.

Abstract

This invention discloses a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data, belonging to the field of geophysical exploration and geological modeling technology. The method includes the following steps: acquiring and preprocessing triangular meshes and source data; extracting normal vectors, strike angles, and dip angles; constructing the geostress tensor; calculating structural surface mechanical parameters, slip angles, and average parameters of the modeling trend surface; defining a buffer zone to filter source data and extract source trend surface parameters; comparing the difference between the normal vectors and slip angles of the two data sources; setting a threshold for iterative model correction; generating a three-dimensional comparison chart to annotate the difference areas and iteration progress; and achieving dynamic matching between the model and seismic activity characteristics. This invention constructs a quantitative correlation between a three-dimensional structural model and source parameters using massive source data, and dynamically iteratively corrects the model based on trend surface differences. This improves the utilization rate of source data and the accuracy of the model, providing high-precision geomechanical support for shale gas development and reservoir safety monitoring projects.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration and geological modeling technology, specifically to a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data. Background Technology

[0002] Three-dimensional structural models are stereoscopic visualization models created using computer graphics and geographic information system (GIS) technologies. They can accurately describe the shape, size, and spatial relationships of structures. Utilizing X, Y, and Z coordinate data, a discrete smooth interpolation algorithm is used to construct triangular meshes to create structural surfaces that conform to geological characteristics, achieving a digital representation of structures. These models are widely used in oil exploration, geothermal development, geological hazard assessment, and scientific research, becoming an indispensable stereoscopic language in the digital age. The core value of three-dimensional structural models lies in their intuitiveness and accuracy. Through 3D structural modeling, designers and engineers can clearly observe the three-dimensional spatial distribution characteristics of structures, construct spatial relationships between structures, optimize drilling plans, and define hazard assessments. Furthermore, the models support virtual simulation and real-time interaction, effectively saving time and costs.

[0003] To address the issues of mismatched accuracy between the 3D structural model and source parameters, and low efficiency in source data utilization during iterative modeling of 3D structural models, existing techniques employ manual qualitative analysis using a limited amount of source mechanism solution data. However, this approach often fails to effectively integrate massive amounts of source mechanism solution data into the modeling process, resulting in a lack of quantitative basis for iterative and progressive approximation of the 3D structural model, significant randomness in the modeling results, and a disconnect from actual seismic activity characteristics. To resolve these issues, a method for iterative modeling of 3D structural models using massive amounts of source mechanism solution data is proposed. Summary of the Invention

[0004] The purpose of this invention is to provide a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data, so as to solve the problems mentioned in the background art.

[0005] To solve the above-mentioned technical problems, the technical solution adopted by this invention is: a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data, comprising the following steps:

[0006] S1. Acquire and preprocess triangular mesh data of structural surfaces and focal mechanism solution data in the three-dimensional structural model;

[0007] S2. Extract the normal vector, orientation angle, and dip angle of each triangle in the modeling trend surface from the preprocessed triangular mesh data;

[0008] S3. Based on the magnitude and direction of the principal stress obtained from the source mechanism solution inversion, construct the geostress tensor, calculate the traction force, normal stress and shear stress on the structural surface, decompose the components of shear stress in the strike vector and dip vector directions, and calculate the slip angle.

[0009] S4. Based on the normal vector and sliding angle of each triangular face, calculate the average normal vector and average sliding angle of the modeling trend surface, define the structural surface buffer, filter the focal mechanism solution data falling into the structural surface buffer, and calculate the average normal vector and average sliding angle of the focal trend surface.

[0010] S5. Compare the average normal vector angle and average slip angle difference between the modeling trend surface and the source trend surface, set the angle threshold, and iteratively update the three-dimensional construction model.

[0011] S6. Generate a spatial comparison map of the three-dimensional structural model and the earthquake source trend surface, and mark the areas of difference and the iteration progress.

[0012] A further improvement to the technical solution of this invention lies in the following: In step S1, the acquisition and preprocessing of the triangular mesh data of the structural surface and the focal mechanism solution data in the three-dimensional structural model includes:

[0013] Three-dimensional geological modeling software is deployed on high-performance computing servers and workstations, and geological exploration data, including borehole data, seismic profile data and ground-penetrating radar scan data, is imported. Professional technicians use the three-dimensional geological modeling software to integrate the geological exploration data and perform three-dimensional geological modeling based on their understanding of geological structures, generate triangular mesh models of structural surfaces, output the node number and three-dimensional coordinates of each triangular surface, and obtain triangular mesh data.

[0014] Broadband seismographs are deployed to form a regional seismic network. The seismographs record the arrival times and amplitudes of P-waves and S-waves during seismic events. Focal mechanism solution inversion software is deployed. This software receives the arrival times and amplitudes of P-waves and S-waves generated during natural earthquakes and inverts the focal parameters from the seismic waveforms. These focal parameters include the focal coordinates P... s =(x s ,y s ,z s ), direction angle Inclination angle δ s and sliding angle r s Obtain focal mechanism solution data;

[0015] The triangular mesh data is topologically checked to correct mesh breaks and overlaps, ensuring that node coordinates are continuous and conform to geological patterns. The cleaned triangular mesh data is then stored in partitions according to structural surfaces.

[0016] The source parameters were unified into an east-north-vertical downward coordinate system. Data with source depths exceeding the range of the three-dimensional structural model were removed. Data records with missing strike angles, dip angles, and slip angles were deleted. Abnormal source parameters deviating more than three times the mean were removed using the standard deviation method. The cleaned source mechanism solution data were stored with spatiotemporal labels and associated with triangular mesh data.

[0017] A further improvement to the technical solution of this invention lies in the following: In step S2, the process of extracting the normal vector, orientation angle, and dip angle of each triangular facet from the preprocessed triangular mesh data includes:

[0018] Based on the coordinates of the three nodes A(x1,y1,z1), B(x2,y2,z2), and C(x3,y3,z3) of each triangular face in the preprocessed triangular mesh data, the eastern component of the normal vector of that triangular face is calculated by vector cross product. North component and vertical component The calculation process is as follows:

[0019]

[0020] Based on the east and north components of the normal vector, calculate the dip azimuth angle θ and the strike angle. The calculation process is as follows:

[0021]

[0022] Where θ is the direction of the inclination of the structural surface;

[0023] Based on the vertical component of the normal vector and normal vector magnitude The inclination angle δ is calculated as follows:

[0024]

[0025] Wherein, normal vector

[0026] A further improvement to the technical solution of this invention lies in the following: In step S3, the process of constructing the geostress tensor, calculating the traction force, normal stress, and shear stress on the structural surface, decomposing the components of the shear stress in the strike vector and dip vector directions, and calculating the slip angle includes:

[0027] Based on the principal stress magnitudes S1≥S2≥S3 and principal stress directions obtained from focal mechanism inversion, each principal stress direction is represented by unit vectors v1, v2, and v3, and the geostress tensor σ is constructed. The calculation process is as follows:

[0028] v x =sinδ·sinθ,vy =sinδ·cosθ,v y =cosδ;

[0029] σ=S1v1v1+S2v2v2+S3v3v3;

[0030] Based on the triangular face normal vector Calculate the traction force on the structural surface And decompose the normal stress σ n And the shear stress τ, the calculation process is as follows:

[0031] t = σ·n;

[0032] σ n =t·n;

[0033] τ=t-σ n n;

[0034] The shear stress is decomposed into a strike vector s representing a unit vector in the horizontal direction and a dip vector d pointing towards the dipping direction of the structural surface. The shear stress is then projected onto the strike and dip directions respectively, and the angle between the shear stress direction and the strike vector is calculated as the slip angle r. The calculation process is as follows:

[0035]

[0036] d = s × n;

[0037]

[0038] Where 0°≤r≤180°, and τ δ These are the components of the shear stress in the strike and dip directions, respectively. Then r = 180° - r.

[0039] A further improvement to the technical solution of this invention lies in the following: In step S4, the process of calculating the average normal vector and average sliding angle of the modeling trend surface based on the normal vector and sliding angle of each triangular face includes:

[0040] In the 3D construction model, a construction surface F is defined, which contains M triangular faces, each of which has a face T. i From its vertex coordinates A i (x1,y1,z1), B i (x2,y2,z2) and C i Composed of (x3, y3, z3), the arithmetic mean of the normal vectors of the triangular faces within the construction surface and the triangular sliding angles is used to obtain the average normal vector of the modeling trend surface. With average sliding angle The calculation process is as follows:

[0041]

[0042] in, and These are the normal vector and sliding angle of the i-th triangle, respectively.

[0043] A further improvement to the technical solution of this invention lies in the following: In step S4, the process of defining a structural surface buffer and filtering focal mechanism solution data falling into the structural surface buffer includes:

[0044] Triangle T i Extending a distance R along the direction of the normal vector forms a three-dimensional spatial region, such that point P = (x p ,y p ,z p If the distance d ≤ R between the triangular faces, then the buffer zone of the constructed face is... Wherein d(p,T) i () indicates the distance from point P to triangle T i The distance;

[0045] Filter the focal mechanism solution data that fall within the structural buffer zone, and calculate P for each focal point. s to the construction surface buffer zone Γ F The minimum distance d(P) s ,Γ F The calculation process is as follows:

[0046]

[0047] If d(P) s ,Γ F If R ≤ R, then the earthquake source is determined to fall within the tectonic buffer zone Γ. F If d(P) s ,Γ F If R > R, then the earthquake source is determined not to have fallen into the tectonic buffer zone Γ. F .

[0048] A further improvement to the technical solution of this invention lies in the following: In step S4, the process of calculating the average normal vector and average slip angle of the focal trend surface based on the focal mechanism solution data of the structural surface buffer includes:

[0049] If N focal mechanism solutions fall within the tectonic surface buffer, then the average normal vector of its focal trend surface is... and average sliding angle The calculation process is as follows:

[0050]

[0051] Among them, based on the strike angle of the epicenter and tilt angle δ s The normal vector of each earthquake source

[0052] A further improvement to the technical solution of this invention lies in the following: In step S5, the process of comparing the difference between the average normal vector angle and the average slip angle between the modeling trend surface and the source trend surface includes:

[0053] Based on the average normal vector and average slip angle of the modeling trend surface and the source trend surface, the angle between their normal vectors is calculated using the dot product formula of the normal vectors. and the absolute difference between the average sliding angles of the two The calculation process is as follows:

[0054]

[0055] in, This is the dot product of the mean normal vector of the earthquake source trend surface and the mean normal vector of the modeling trend surface. and These are the magnitudes of the average normal vector of the source trend surface and the average normal vector of the modeling trend surface, respectively.

[0056] A further improvement to the technical solution of this invention lies in the following: In step S5, the process of setting an angle threshold and iteratively updating the three-dimensional construction model includes:

[0057] A preset angle threshold η is used to determine whether the modeling trend surface matches the earthquake source trend surface. At that time, it was considered that the structural surface obtained by modeling did not match the focal mechanism solution data. It is necessary to select the focal mechanism solution data characteristics based on the angle of the modeling trend surface and carry out iterative modeling control.

[0058] When the difference between the average normal angle and the average sliding angle of two trend surfaces is less than one of the angle thresholds, the average normal angle of the trend surface takes priority. When the average normal angle is greater than the threshold angle, it is necessary to consider whether to iterate. When the average normal angle is 4 to 5 times the threshold angle, iterative modeling is required.

[0059] The iterative modeling process includes adjusting the node coordinates of the triangular mesh data in the 3D structural model based on the directional difference of the average normal vector angle, so that the normal vectors of the modeling trend surface and the source trend surface are close. If the difference in the average slip angle is greater than the angle threshold, the slip angle of the structural surface is adjusted by adjusting the ratio of the shear stress components in the strike and dip directions, so that it is close to the slip angle of the source trend surface. The updated triangular mesh data iteratively generates new modeling trend surface normal vectors and slip angles until... At this point, it is assumed that the structural surface obtained from the modeling is consistent with the source mechanism solution data, and no iteration is required. The three-dimensional structural model at this time is then output.

[0060] A further improvement to the technical solution of this invention lies in the following: In step S6, the process of generating a spatial comparison map of the three-dimensional structural model and the earthquake source trend surface, and marking the difference areas and iteration progress, includes:

[0061] Based on the triangular mesh data of the 3D structural model and the parameters of the source trend surface, the geometry of the structural surface is rendered by a 3D visualization engine, and the source mechanism solution point cloud data is superimposed. Color coding is used to represent the slip angle and the source depth. The normal vectors of the modeling trend surface and the source trend surface are projected into the same space, and the difference areas are marked with color levels. According to the difference between the angle between the local average normal vector and the average slip angle, the areas exceeding the angle threshold are marked as highlighted blocks, and the global iteration progress and the proportion of areas that have not met the standards are displayed, providing a visual criterion for model correction.

[0062] After each iteration of the 3D model construction, the coordinates of the triangular mesh nodes are automatically updated and the 3D comparison image is re-rendered. The annotations of the difference areas and the convergence curve are refreshed in real time. Users can manually adjust the angle threshold and update the highlighted area range instantly. Clicking on a single triangular face allows you to view its average normal vector angle, average sliding angle difference, and iteration history data.

[0063] Due to the adoption of the above technical solution, the technical progress achieved by this invention compared to the prior art is as follows:

[0064] 1. This invention provides a method for iterative modeling of a three-dimensional structural model using massive source mechanism solution data, which solves the problem of mismatch between the accuracy of the three-dimensional model and the source parameters. By constructing the geostress tensor and calculating the mechanical parameters of the structural surface, the non-uniform triangular mesh model and the source parameters are uniformly mapped to the same platform, realizing the consistency of geometric and mechanical parameters and improving the accuracy of the model.

[0065] 2. This invention provides a method for iterative modeling of a three-dimensional structural model using massive source mechanism solution data, which significantly improves the efficiency of source data utilization. By defining a structural surface buffer to filter spatially correlated source data and calculating source trend surface parameters, massive source mechanism solution data is automatically integrated into the modeling process, thereby improving data utilization.

[0066] 3. This invention provides a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data. It provides quantitative iteration criteria and dynamic closed-loop correction capabilities. An angle threshold is set based on the angle between the normal vectors. When the difference exceeds the threshold, the coordinates of the triangular mesh nodes are automatically adjusted to drive the model to approximate the actual seismic activity characteristics, reduce the number of iterations, and improve the model convergence speed. It can be applied to engineering scenarios such as shale gas development and reservoir safety monitoring. Attached Figure Description

[0067] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this invention. For those skilled in the art, other drawings can be obtained based on these drawings.

[0068] Figure 1 This is a flowchart of the present invention. Detailed Implementation

[0069] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, 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, 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.

[0070] Examples, such as Figure 1 As shown, this invention provides a method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data, comprising the following steps:

[0071] S1. Acquire and preprocess triangular mesh data and focal mechanism solution data of structural surfaces in the 3D structural model. Deploy 3D geological modeling software on high-performance computing servers and workstations, import geological exploration data, including borehole data, seismic profile data, and ground-penetrating radar scan data. Professional technicians use the 3D geological modeling software to integrate the geological exploration data and perform 3D geological modeling based on their understanding of geological structures, generating triangular mesh models of structural surfaces, outputting the node number and 3D coordinates of each triangular facet, acquiring triangular mesh data, deploying broadband seismometers to form a regional seismic network, recording the arrival times and amplitudes of P-waves and S-waves of seismic events, deploying focal mechanism solution inversion software, receiving the arrival times and amplitudes of P-waves and S-waves generated during natural earthquakes, and inverting the focal parameters from the seismic waveforms. The focal parameters include the focal coordinates P... s =(x s ,y s ,z s ), direction angle Inclination angle δ s and sliding angle r sThe process involves acquiring focal mechanism solution data, performing topological checks on the triangular mesh data to correct mesh breaks and overlaps, ensuring continuous node coordinates that conform to geological patterns, storing the cleaned triangular mesh data in partitions according to structural surfaces, unifying the focal parameters into an east-north-vertical downward coordinate system, removing data whose focal depth exceeds the range of the three-dimensional structural model, deleting missing strike angle, dip angle, and slip angle data records, and removing abnormal focal parameters that deviate from the mean by more than 3 times using the standard deviation method. The cleaned focal mechanism solution data is stored with spatiotemporal labels and associated with the triangular mesh data.

[0072] S2. Extract the normal vector, orientation angle, and dip angle of each triangle in the modeling trend surface from the preprocessed triangular mesh data. Based on the three node coordinates A(x1,y1,z1), B(x2,y2,z2), and C(x3,y3,z3) of each triangle in the preprocessed triangular mesh data, calculate the eastward component of the normal vector of the triangle by vector cross product. North component and vertical component The calculation process is as follows:

[0073]

[0074] Based on the east and north components of the normal vector, calculate the dip azimuth θ and strike angle. The calculation process is as follows:

[0075]

[0076] Where θ is the direction of the tilt of the construction surface, based on the vertical component of the normal vector. and normal vector magnitude The inclination angle δ is calculated as follows:

[0077]

[0078] Wherein, normal vector

[0079] S3. Based on the magnitude and direction of the principal stresses obtained from the focal mechanism inversion, construct the geostress tensor, calculate the traction force, normal stress, and shear stress on the structural surface, decompose the shear stress components in the strike and dip vectors, and calculate the slip angle. Based on the magnitudes S1≥S2≥S3 and directions of the principal stresses obtained from the focal mechanism inversion, represent each principal stress direction with unit vectors v1, v2, and v3, and construct the geostress tensor σ. The calculation process is as follows:

[0080] v x =sinδ·sinθ,v y =sinδ·cosθ,v y=cosδ;

[0081] σ=S1v1v1+S2v2v2+S3v3v3;

[0082] Based on the triangular face normal vector Calculate the traction force on the structural surface And decompose the normal stress σ n And the shear stress τ, the calculation process is as follows:

[0083] t = σ·n;

[0084] σ n =t·n;

[0085] τ=t-σ n n;

[0086] The shear stress is decomposed into a strike vector s representing a unit vector in the horizontal direction and a dip vector d pointing towards the dipping direction of the structural surface. The shear stress is then projected onto the strike and dip directions respectively, and the angle between the shear stress direction and the strike vector is calculated as the slip angle r. The calculation process is as follows:

[0087]

[0088] d = s × n;

[0089]

[0090] Where 0°≤r≤180°, and τ δ These are the components of the shear stress in the strike and dip directions, respectively. Then r = 180° - r;

[0091] S4. Based on the normal vector and slip angle of each triangular facet, calculate the average normal vector and average slip angle of the modeling trend surface. Define a structural surface buffer, filter the focal mechanism solution data falling into the structural surface buffer, and calculate the average normal vector and average slip angle of the focal trend surface. In the three-dimensional structural model, define a structural surface F, which contains M triangular faces, each triangular facet T... i From its vertex coordinates A i (x1,y1,z1), B i (x2,y2,z2) and C i Composed of (x3, y3, z3), the arithmetic mean of the normal vectors of the triangular faces within the construction surface and the triangular sliding angles is used to obtain the average normal vector of the modeling trend surface. With average sliding angle The calculation process is as follows:

[0092]

[0093] in, and Let T be the normal vector and the sliding angle of the i-th triangle face, respectively. i Extending a distance R along the direction of the normal vector forms a three-dimensional spatial region, such that point P = (x p ,y p ,z p If the distance d ≤ R between the triangular faces, then the buffer zone of the constructed face is... Wherein d(p,T) i () indicates the distance from point P to triangle T i The distance is used to filter the focal mechanism solution data that fall within the structural buffer zone, and the distance to each focal point P is calculated. s to the construction surface buffer zone Γ F The minimum distance d(P) s ,Γ F The calculation process is as follows:

[0094]

[0095] If d(P) s ,Γ F If R ≤ R, then the earthquake source is determined to fall within the tectonic buffer zone Γ. F If d(P) s ,Γ F If R > R, then the earthquake source is determined not to have fallen into the tectonic buffer zone Γ. F If N focal mechanism solutions fall within the tectonic surface buffer, then the average normal vector of its focal trend surface is... and average sliding angle The calculation process is as follows:

[0096]

[0097] Among them, based on the strike angle of the epicenter and tilt angle δ s The normal vector of each earthquake source

[0098] S5. Compare the difference between the average normal vector angle and the average slip angle between the modeling trend surface and the source trend surface, set an angle threshold, iteratively update the 3D construction model, and calculate the angle between the normal vectors based on the average normal vectors and average slip angles of the modeling trend surface and the source trend surface using the dot product formula of the normal vectors. and the absolute difference between the average sliding angles of the two The calculation process is as follows:

[0099]

[0100] in, This is the dot product of the mean normal vector of the earthquake source trend surface and the mean normal vector of the modeling trend surface. and These are the magnitudes of the average normal vector of the focal trend surface and the average normal vector of the modeling trend surface, respectively. A preset angle threshold η is used to determine whether the modeling trend surface and the focal trend surface match. If the modeled structural surface is found to be inconsistent with the focal mechanism solution data, iterative modeling control is required based on the angle between the modeling trend surfaces and the focal mechanism solution data characteristics. When the difference between the average normal angle and the average slip angle of the two trend surfaces is less than one of the angle thresholds, the average normal angle of the trend surface is prioritized. When the average normal angle is greater than the threshold angle, iteration needs to be considered. When the average normal angle is 4 to 5 times the threshold angle, iterative modeling is required. The iterative modeling process includes adjusting the node coordinates of the triangular mesh data in the 3D structural model based on the directional difference of the average normal angle to make the normal vectors of the modeling trend surface and the focal trend surface approximate each other. If the difference in the average slip angle is greater than the angle threshold, the slip angle of the structural surface is adjusted by adjusting the ratio of the shear stress components in the strike and dip directions to make it closer to the slip angle of the focal trend surface. The updated triangular mesh data iteratively generates new modeling trend surface normal vectors and slip angles until... At this point, it is assumed that the structural surface obtained from the modeling is consistent with the focal mechanism solution data, and no iteration is required; the three-dimensional structural model at this time is then output.

[0101] S6. Generate a spatial comparison map of the 3D structural model and the source trend surface, marking the difference areas and iteration progress. Based on the triangular mesh data of the 3D structural model and the parameters of the source trend surface, render the geometry of the structural surface through a 3D visualization engine, overlay the source mechanism solution point cloud data, and use color coding to represent the slip angle and source depth. Project the normal vector of the modeling trend surface and the normal vector of the source trend surface into the same space, and mark the difference areas with color levels. According to the difference between the local average normal vector angle and the average slip angle, the area exceeding the angle threshold is marked as a highlight block, and the global iteration progress and the proportion of the non-compliant area are displayed, providing a visual criterion for model correction. After each iteration of the 3D structural model, the coordinates of the triangular mesh nodes are automatically updated and the 3D comparison map is re-rendered, and the marking of the difference areas and the convergence curve are refreshed in real time. Users can manually adjust the angle threshold and update the highlighted area range in real time. Clicking on a single triangular face can view its average normal vector angle, average slip angle difference, and iteration history data.

[0102] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data, characterized in that, Includes the following steps: S1. Acquire and preprocess triangular mesh data of structural surfaces and focal mechanism solution data in the three-dimensional structural model; S2. Extract the normal vector, orientation angle, and dip angle of each triangle in the modeling trend surface from the preprocessed triangular mesh data; S3. Based on the magnitude and direction of the principal stress obtained from the source mechanism solution inversion, construct the geostress tensor, calculate the traction force, normal stress and shear stress on the structural surface, decompose the components of shear stress in the strike vector and dip vector directions, and calculate the slip angle. S4. Based on the normal vector and sliding angle of each triangular face, calculate the average normal vector and average sliding angle of the modeling trend surface, define the structural surface buffer, filter the focal mechanism solution data falling into the structural surface buffer, and calculate the average normal vector and average sliding angle of the focal trend surface. S5. Compare the average normal angle and average slip angle difference between the modeling trend surface and the source trend surface, set an angle threshold, and iteratively update the 3D construction model. The process of comparing the average normal angle and average slip angle difference between the modeling trend surface and the source trend surface includes: Based on the average normal vector and average slip angle of the modeling trend surface and the source trend surface, the angle between their normal vectors is calculated using the dot product formula of the normal vectors. and the absolute difference between the average sliding angles of the two The calculation process is as follows: in, This is the dot product of the mean normal vector of the earthquake source trend surface and the mean normal vector of the modeling trend surface. and These are the magnitudes of the mean normal vector of the epicenter trend surface and the mean normal vector of the modeling trend surface, respectively. The process of setting an angle threshold and iteratively updating the 3D construction model includes: A preset angle threshold η is used to determine whether the modeling trend surface matches the earthquake source trend surface. At that time, it was considered that the structural surface obtained by modeling did not match the focal mechanism solution data. It is necessary to select the focal mechanism solution data characteristics based on the angle of the modeling trend surface and carry out iterative modeling control. When the difference between the average normal angle and the average sliding angle of two trend surfaces is less than one of the angle thresholds, the average normal angle of the trend surface takes priority. When the average normal angle is greater than the threshold angle, it is necessary to consider whether to iterate. When the average normal angle is 4 to 5 times the threshold angle, iterative modeling is required. The iterative modeling process includes adjusting the node coordinates of the triangular mesh data in the 3D structural model based on the directional difference of the average normal vector angle, so that the normal vectors of the modeling trend surface and the source trend surface are close. If the difference in the average slip angle is greater than the angle threshold, the slip angle of the structural surface is adjusted by adjusting the ratio of the shear stress components in the strike and dip directions, so that it is close to the slip angle of the source trend surface. The updated triangular mesh data iteratively generates new modeling trend surface normal vectors and slip angles until... At this point, it is assumed that the structural surface obtained from the modeling is consistent with the focal mechanism solution data, and no iteration is required; the three-dimensional structural model at this time is then output. S6. Generate a spatial comparison map of the three-dimensional structural model and the earthquake source trend surface, and mark the areas of difference and the iteration progress.

2. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 1, characterized in that: In S1, the acquisition and preprocessing of triangular mesh data of structural surfaces and focal mechanism solution data in the three-dimensional structural model includes: Three-dimensional geological modeling software is deployed on high-performance computing servers and workstations. Geological exploration data, including borehole data, seismic profile data and ground-penetrating radar scan data, is imported. Using the three-dimensional geological modeling software, the geological exploration data is integrated and three-dimensional geological modeling is performed based on the understanding of geological structures. A triangular mesh model of the structural surface is generated, and the node number and three-dimensional coordinates of each triangular surface are output to obtain the triangular mesh data. Broadband seismographs are deployed to form a regional seismic network. The seismographs record the arrival times and amplitudes of P-waves and S-waves during seismic events. Focal mechanism solution inversion software is deployed. This software receives the arrival times and amplitudes of P-waves and S-waves generated during natural earthquakes and inverts the focal parameters from the seismic waveforms. These focal parameters include the focal coordinates P... s =(x s ,y s ,z s ), direction angle Inclination angle δ s and sliding angle r s Obtain focal mechanism solution data; The triangular mesh data is topologically checked to correct mesh breaks and overlaps, ensuring that node coordinates are continuous and conform to geological patterns. The cleaned triangular mesh data is then stored in partitions according to structural surfaces. The source parameters were unified into an east-north-vertical downward coordinate system. Data with source depths exceeding the range of the three-dimensional structural model were removed. Data records with missing strike angles, dip angles, and slip angles were deleted. Abnormal source parameters deviating more than three times the mean were removed using the standard deviation method. The cleaned source mechanism solution data were stored with spatiotemporal labels and associated with triangular mesh data.

3. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 2, characterized in that: In step S2, the process of extracting the normal vector, orientation angle, and dip angle of each triangle facet in the modeling trend surface from the preprocessed triangular mesh data includes: Based on the coordinates of the three nodes A(x1,y1,z1), B(x2,y2,z2), and C(x3,y3,z3) of each triangular face in the preprocessed triangular mesh data, the eastern component of the normal vector of that triangular face is calculated by vector cross product. North component and vertical component The calculation process is as follows: Based on the east and north components of the normal vector, calculate the dip azimuth θ and strike angle. The calculation process is as follows: Where θ is the direction of the inclination of the structural surface; Based on the vertical component of the normal vector and normal vector magnitude The inclination angle δ is calculated as follows: Wherein, normal vector 4. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 3, characterized in that: In S3, the process of constructing the geostress tensor, calculating the traction force, normal stress, and shear stress on the structural surface, decomposing the components of the shear stress in the strike and dip vectors, and calculating the slip angle includes: Based on the principal stress magnitudes S1≥S2≥S3 and principal stress directions obtained from focal mechanism inversion, each principal stress direction is represented by unit vectors v1, v2, and v3, and the geostress tensor σ is constructed. The calculation process is as follows: v x =sinδ·sinθ,v y =sinδ·cosθ,v y =cosδ; σ=S1v1v1+S2v2v2+S3v3v3; Based on the triangular face normal vector Calculate the traction force t on the structural surface and decompose the normal stress σ. n And the shear stress τ, the calculation process is as follows: t = σ·n; s n =t·n; τ=t-σ n n; The shear stress is decomposed into a strike vector s representing a unit vector in the horizontal direction and a dip vector d pointing towards the dipping direction of the structural surface. The shear stress is then projected onto the strike and dip directions respectively, and the angle between the shear stress direction and the strike vector is calculated as the slip angle r. The calculation process is as follows: d = s × n; Where 0°≤r≤180°, and τ δ These are the components of the shear stress in the strike and dip directions, respectively. Then r = 180° - r.

5. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 4, characterized in that: In step S4, the process of calculating the average normal vector and average sliding angle of the modeling trend surface based on the normal vector and sliding angle of each triangular face includes: In the 3D construction model, a construction surface F is defined, which contains M triangular faces, each of which has a face T. i From its vertex coordinates A i (x1,y1,z1), B i (x2,y2,z2) and C i Composed of (x3, y3, z3), the arithmetic mean of the normal vectors of the triangular faces within the construction surface and the triangular sliding angles is used to obtain the average normal vector of the modeling trend surface. With average sliding angle The calculation process is as follows: in, and These are the normal vector and sliding angle of the i-th triangle, respectively.

6. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 5, characterized in that: In S4, the process of defining a structural surface buffer and filtering focal mechanism solution data falling into the structural surface buffer includes: Triangle T i Extending a distance R along the direction of the normal vector forms a three-dimensional spatial region, such that point P = (x p ,y p ,z p If the distance d ≤ R between the triangular faces, then the buffer zone of the constructed face is... Wherein d(p,T) i () indicates the distance from point P to triangle T i The distance; Filter the focal mechanism solution data that fall within the structural buffer zone, and calculate P for each focal point. s to the construction surface buffer zone Γ F The minimum distance d(P) s ,Γ F The calculation process is as follows: If d(P) s ,Γ F If R ≤ R, then the earthquake source is determined to fall within the tectonic buffer zone Γ. F If d(P) s ,Γ F If R > R, then the earthquake source is determined not to have fallen into the tectonic buffer zone Γ. F .

7. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 6, characterized in that: In step S4, the process of calculating the average normal vector and average slip angle of the focal trend surface based on the focal mechanism solution data of the structural surface buffer includes: If N focal mechanism solutions fall within the tectonic surface buffer, then the average normal vector of its focal trend surface is... and average sliding angle The calculation process is as follows: Among them, based on the strike angle of the epicenter and tilt angle δ s The normal vector of each earthquake source 8. The method for iteratively modeling a three-dimensional structural model using massive source mechanism solution data according to claim 7, characterized in that: In step S6, the process of generating a spatial comparison map of the three-dimensional structural model and the seismic source trend surface, and marking the difference areas and iteration progress, includes: Based on the triangular mesh data of the 3D structural model and the parameters of the source trend surface, the geometry of the structural surface is rendered by a 3D visualization engine, and the source mechanism solution point cloud data is superimposed. Color coding is used to represent the slip angle and the source depth. The normal vectors of the modeling trend surface and the source trend surface are projected into the same space, and the difference areas are marked with color levels. According to the difference between the angle between the local average normal vector and the average slip angle, the areas exceeding the angle threshold are marked as highlighted blocks, and the global iteration progress and the proportion of areas that have not met the standards are displayed, providing a visual criterion for model correction. After each iteration of the 3D model construction, the coordinates of the triangular mesh nodes are automatically updated and the 3D comparison image is re-rendered. The annotations of the difference areas and the convergence curve are refreshed in real time. Users can manually adjust the angle threshold and update the highlighted area range instantly. Clicking on a single triangular face allows you to view its average normal vector angle, average sliding angle difference, and iteration history data.

Images