94 results about "Fault plane" patented technology

The invention discloses a discrete element-based three-dimensional fault / fracture network modeling method and device under a backlash pushing and covering system, and relates to the technical field of three-dimensional fault modeling, and the method comprises the steps: building an initialized fault plane model based on the seismic logging information parameters of a target region; extracting stratum parameters of the target area, and determining a deformation process in combination with an equilibrium profile recovery method; the method comprises the following steps: calibrating mesoscopic parameters of discrete elements according to stratum parameters, depositing a first discrete element model, optimizing and generating a second discrete element model, setting boundary conditions to carry out back-flushing pushing-covering simulation, comparing an output result of the second discrete element model with an actual geological section, adjusting model parameters until the output result and the actual geological section are matched, and establishing a three-dimensional fault and fracture network. And analyzing and establishing spatial distribution characteristics of the target area, and carrying out dynamic evaluation and trend prediction on the activity of the three-dimensional fault. According to the method, the geomechanical model is fused, the subjectivity of manual interpretation is overcome, and the problems of fault tip extension and the like in grid reconstruction are effectively avoided.

The invention discloses a fault slip model smoothing method based on SAR observation data, and the method comprises the steps: firstly, obtaining multi-source observation data (including SAR, GNSS and optical images), and carrying out the joint calculation, and generating a high-resolution earth surface deformation field; the method comprises the following steps: combining a remote sensing image, visually interpreting and drawing a fracture trace, then carrying out data downsampling along an inclination angle direction, and finally constructing an irregular fault geometric model by using a TIN method so as to obtain an initial fault plane model. On this basis, a Laplacian operator suitable for triangular dislocation is defined, a slip inversion model is constructed in combination with an elastic dislocation theory, and a smoothing factor is introduced to balance slip details and calculation efficiency. And finally, optimizing a fault inclination angle and a smoothing factor through a grid search method, evaluating fitting errors and slippage roughness under different combinations, and determining optimal parameters. The final model not only comprises a triangular connection relation, but also provides slip amount distribution and fault dip angle information, and an important basis is provided for seismic activity evaluation and disaster response.

The invention discloses a method and a system for predicting the vertical connectivity of a fault in a petroliferous basin, and relates to the technical field of fault connectivity detection.The method comprises the steps of determining a fault vertical connectivity coefficient of a target position according to influence factors of the fault vertical connectivity of the target position; the influence factors comprise the thickness of the fractured cover layer, the shale content, the normal stress of the fault surface, the formation fluid pressure, the brittleness degree of the cover layer and the fault displacement of the fault; the target position is the target position of a non-drilling area in the petroliferous basin; and predicting the vertical connectivity of the target position according to the fault vertical connectivity coefficient of the target position. According to the method, the identification accuracy of the vertical connectivity of the fault in the petroliferous basin is improved.

The invention discloses an internal space relation speculation method and system based on a rock mass structural surface. The method comprises the steps that based on rock mass structural plane recognition results of a ground plane, a joint plane and a fault plane, correlation of rock mass structural plane surface features and internal space relations of the rock mass structural plane surface features is constructed, and surface constraint conditions of the ground plane, the joint plane and the fault plane are determined; based on key control point setting, thickness constraint and geometric calculation of the triangulated irregular network, realizing thickness-invariant extension of the rock mass structural surface; based on thickness constraint, triangularization curved surface mesh deformation is achieved, and then thickness change extension between different depth layers is achieved; and the triangular irregular net is cut based on the intersecting lines of the intersecting rock mass structural surfaces, and stratum pinching is achieved. According to the method, occurrence information such as the inclination angle of the rock mass internal structure and the internal space relation of the rock mass structural plane are obtained by extending the surface of the structural plane, and then a rock mass internal space relation speculation result serves as an internal constraint to provide data support for subsequent three-dimensional modeling.

The invention provides a shale reservoir hydraulic fracturing wellbore casing deformation risk judgment method which comprises the steps that geological parameters and principal stress of a target area are obtained, and the principal stress comprises horizontal principal stress, maximum principal stress and vertical principal stress. And determining an included angle between the fault plane and the hydraulic fracture plane according to the normal vector of the fault plane and the normal vector of the hydraulic fracture plane. And calculating by substituting the geological parameters, the principal stress and the included angle into an established stress calculation model of the fault plane to obtain shear stress and normal stress of the fault plane in the target area. And calculating by substituting the shear stress and the normal stress into the established displacement component calculation model to obtain the slip component of the fault plane in the target area. And performing safety early warning based on the risk level. The fault slip amount is calculated by accurately identifying the stratum physical property parameters and the fault structure, so that the casing deformation risk is judged.

The invention discloses a method, a device, equipment and a medium for determining magma distribution of an intrusive fault zone, and relates to the technical field of oil and gas field exploration and development, and the method comprises the following steps: loading geological data of a target area; based on the geological data, constructing a fault plane of the magma intrusive fault zone in the target area; calculating a magma distribution indication factor of the magma intrusive fault zone in the target area; and determining a magma distribution position in the magma intrusive fault zone in the target area. According to the method, the value characteristic and the frequency characteristic of the earthquake information amount of the magma invading body are considered, the magma distribution indication factor S of the igneous rock invading fault zone is provided by combining the distribution characteristic of the fracture, the purpose of quantitatively predicting the magma distribution range in the fault zone is achieved, and the method can be used for predicting the closed and oil and gas dredging part in the fault zone. The brought prediction effect is more in line with the actual geological law; and a large number of geological parameters do not need to be counted, research can be directly carried out on seismic data, and the implementation process is more intuitive and rapid.

The invention belongs to the technical field of earthquake simulation tests, and particularly relates to a fault dislocation testing device and method for a three-dimensional earthquake test, and the fault dislocation testing device is composed of two beveling models which are in the same form and are beveled on a cube along diagonal lines of an upper square surface and a lower square surface of the cube. The contact surface of the two beveled models, namely the opposite surface formed by beveling the two beveled models along the diagonal line, is a fault surface, and the fault models are arranged according to different test methods. The experimental device comprises the fault models, a strain observation system, an acoustic emission monitoring system, an active acoustic emission observation system and a laser total reflection fault monitoring system; according to the invention, a laboratory earthquake experiment can be realized, a laser total reflection fault observation system is researched and developed on the basis of a stress-strain observation device and an active and passive acoustic emission device, and synchronous monitoring of fault stress-strain, acoustic emission signals, displacement and fault surface contact conditions can be realized at the same time.

The invention discloses a method and device for constructing a three-dimensional fault plane according to an earthquake catalog, and the method comprises the steps: extracting the position information of an earthquake source according to the earthquake catalog of a target region, obtaining the density information of each data point, and generating a density map; based on the density information and the density map, utilizing an OPTICS algorithm to carry out preliminary clustering analysis on the seismic source position information, outputting an accessibility map, and determining a clustering number N; calculating the local density and the center offset distance of each seismic source point based on the Euclidean distance and the truncation distance between every two seismic source points, generating a decision diagram, and selecting the first N seismic source points with the maximum product of the local density and the center offset distance as a clustering center; and PCA analysis and PCA plane determination are carried out on the DPC classification result, and plane features of the fault plane are obtained when set conditions are satisfied. According to the method, an OPTICS algorithm, a DPC algorithm and PCA analysis are combined, and end-to-end parameterization-free three-dimensional fault plane modeling is achieved only by means of seismic source position information.

A test method and apparatus for the evolution of a neighborhood stress field during fault instability. The test method comprises: step S1: manufacturing a substrate (101); S2, selecting a pre-configured surface (101a) of the substrate to fix a strain sensing member (102); step S3, fixing the substrate (101) to the position of a fault in an experimental model (103), and determining a coordinate position of the substrate (101) on the basis of an inclination angle of the fault; step S4, connecting the strain sensing member (102) to a strain gauge; step S5, applying in-situ triaxial stress to the experimental model (103), and applying disturbance to the experimental model (103) when the experimental model reaches equilibrium; and step S6, calculating a stress value on the basis of a strain component measured by the strain gauge, and capturing the dynamic evolution of a neighborhood stress field during fault plane instability caused by the disturbance.

According to the main shock and aftershock time sequence rapid construction method and device provided by the invention, the influence of similar source generation mechanisms of main shock and aftershock, similar propagation paths, the same field effect and the like is fully considered, and the main shock and aftershock time sequence is rapidly and accurately constructed. The method comprises the steps that 1, the most dangerous fault plane and the spatial position and size of the fault plane are determined according to seismic fault data, sub-events with different circle center positions are generated randomly, and the size of each sub-event is represented by a radius; 2, determining the magnitude of the aftershock according to the set main shock magnitude and the main and aftershock power law relationship; step 3, according to the magnitude of the main earthquake and the aftershock, using rectangles consistent with the size of a main earthquake source and aftershock source to perform frame selection on the fault plane, and using sub-events contained in the rectangles to form an earthquake source model so as to obtain an earthquake source finite fault model of the main earthquake and the aftershock; 4, a deterministic Green function is obtained through decoupling; and step 5, carrying out convolution on the seismic source finite fault model and the deterministic Green function to obtain a final main shock and aftershock sequence.

The invention discloses an ArcGIS-based fault complex structure processing method in geological modeling, which comprises the following steps: collecting related geological data of a research area, laying section lines according to modeling precision and certain density when actual drilling data is lacked, laying virtual drilling points on the section lines, and carrying out modeling on the virtual drilling points; calculating modeling data, including thickness and burial depth, of the underlying stratum of the virtual drill hole through the geological profile map; in combination with geological data and field survey data, finding out basic features of a fault in a modeling area, including a fault angle (alpha) and a stratum thickness (h), and before modeling, arranging a section which has a certain density and passes through the fault according to the step S1; taking the fault plane as a boundary of modeling, and performing independent modeling on geologic bodies on two sides of the fault plane based on ArcGIS; and after the geologic models on the two sides of the fault are built according to the previous steps, the three-dimensional geologic models on the two sides are combined together, so that a complete three-dimensional geologic model capable of reflecting the fault characteristics can be formed.

The invention discloses a data extraction method and system for digital geological mineral exploration, and relates to the technical field of geological mineral exploration, and the method comprises the steps: collecting the original data of a geological mineral exploration site based on a spatial attitude sensing geological detection device; based on the spatial dimension, the angle dimension and the time dimension, preprocessing the collected original data of the geological mineral exploration site to obtain optimized data of the geological mineral exploration site; key geological information is extracted from the optimization data of the geological mineral exploration site to generate geological mineral three-dimensional space features; determining a mineral product space distribution model of the exploration area according to the geological mineral product three-dimensional space characteristics; the spatial information of the geologic structure comprises three-dimensional arrangement information of a stratum, three-dimensional arrangement information of a fracture surface and three-dimensional arrangement information of an ore body.

The invention provides a roadway geology visual modeling method and system under a joint method, and the method comprises the steps: constructing a stratum curved surface model through employing a Kriging interpolation algorithm and an inverse distance weighting method, and employing a TIN-based surface model modeling method; acquiring intersecting lines of the fault and the stratum when each stratum is cut by the same fault so as to determine a control point set of the fault plane, and constructing a three-dimensional geologic model of the fault plane; taking a contour line of a local coordinate system established on each control point along the center line as a section, and carrying out tunneling roadway geometric modeling by using a section expansion modeling method; and performing dynamic prediction and model updating based on the static data of the modeling area and the dynamic data collected in the tunneling process to realize model visualization. According to the method, the form and internal attributes of the geologic body can be taken into consideration, hidden information and feature representation of data can be mined while sample data is learned by means of deep learning, and geological information can be learned and progressively predicted more accurately.

The invention provides a seismic data fault enhancement method based on an optimal local space-time vector, and relates to the technical field of seismic data processing, and the method comprises the steps: calculating and correcting a seismic data event inclination angle, and obtaining a corrected event inclination angle; based on the corrected event inclination angle, obtaining a fault enhancement local space-time vector, and calculating a fault enhancement coefficient and a fault enhancement distance truncation threshold; calculating the distance between the center reference waveform vector and the local space-time vector; extracting an optimal local space-time vector according to a distance truncation threshold value; and finally, carrying out weighted stacking on the seismic data by utilizing the optimal local space-time vector so as to enhance seismic fault features. The method effectively overcomes the defects that an existing reflection seismic data extraction and enhancement method is too strong in dependence on seismic fault information and damages seismic faults, and can remarkably enhance discontinuous characteristics of stratums, so that fault planes are clearer, and oil-gas exploration and development requirements are better met.

The present invention discloses a method for fluid injection-induced earthquake experiments under true triaxial stress conditions. The method first creates a cubic rock model containing a fault, obliquely cuts it to form a fault plane, lays out water injection and water pressure detection holes, and injects water to saturate the fault plane. A true triaxial device is then used to apply three-dimensional unequal stresses and maintain stability. Displacement is then used to drive the fault to spontaneous instability, and the inter-seismic instability period is calculated. Finally, a cyclic or constant-rate water injection scheme is set to simulate industrial fluid activity and explore the mechanism of water injection-induced earthquakes. This method for fluid injection-induced earthquake experiments under true triaxial stress conditions can accurately simulate actual stratum stress and water injection conditions, providing a reliable means for studying water injection-induced earthquakes and having significant implications for earthquake prediction, geological disaster prevention and control, and related industrial safety.

The invention discloses a fault activation monitoring method and system based on an optical fiber sensor. The method comprises the following steps that S1, a fault model is prepared; s2, constructing an experimental sample; s3, preparing and laying fault mud; s4, arranging optical fiber sensors: laying the optical fiber sensors on the surface of the fault surface in a snake shape, a grid shape or a spiral shape, cutting a groove in the outer surface of a sample, and fixing optical fibers on a steel pipe embedded into the groove for monitoring the stress, strain and pore pressure change of the fault surface and shaft deformation in real time; s5, performing experimental simulation and data acquisition; s6, data analysis and early warning: performing fusion analysis on the collected data based on a deep learning algorithm, predicting fault instability risk and generating an early warning signal; by simulating the real stratum condition and the fluid injection process and combining the high-precision optical fiber sensor to monitor the multi-physical field information in the fault zone, technical support is provided for effectively controlling the induced earthquake risk and guaranteeing safe and stable operation of underground engineering.

The present invention provides a method for restoring the slip trajectory of dip-slip faults, including: Step 1, obtaining the structural map of the bottom interface of the target interval, the structural map of the boundary fault plane, and their combined mosaic; Step 2, establishing a profile of the basic interpretation results along the survey line direction; Step 3, obtaining the true dip angle of the fault plane; Step 4, obtaining the spatial position and true length of the fault; Step 5, obtaining the vertical displacement and dip displacement of the fault breakpoint in space; Step 6, obtaining the strike displacement of the fault breakpoint along the fault plane; Step 7, constructing the fault slip trajectory; Step 8, conducting a spatio-temporal displacement difference evaluation. This method for restoring the slip trajectory of dip-slip faults can quantitatively obtain the three-dimensional displacement components of the hanging wall's dip-slip activity along the fault plane, accurately describe the slip trajectory and locate the initial breakpoint position, realize the quantitative characterization of the activity difference characteristics of the spatial target fault, and lay a foundation for the restoration of the prototype basin and its application in oil and gas exploration.

The invention relates to the field of data processing systems suitable for shale gas development management and prediction, in particular to a shale gas-containing fault plane analysis and evaluation system based on seismic data analysis, which comprises a data acquisition module, a fault feature analysis module, a fault analysis module and an enrichment effect evaluation module. According to the method, the seismic data are collected and aggregated, the stratum attribute data set is generated, and the fault features are extracted and clustered, so that the revelation precision of the fault spatial distribution characteristics is improved. And a topological model is further constructed, and the connectivity is analyzed by calculating the connectivity and porosity between the intersection points, so that the understanding of the underground structure is enhanced, and the position and scale evaluation of the shale gas enrichment area is more accurate. In addition, the relation between the fault structure and the gas enrichment area is analyzed, the efficiency of the gas flow channel is calculated, and a scientific basis is provided for gas field development decision making, so that the dominant location is predicted, and the development plan is optimized.

The application provides a seismic motion simulation method based on a stochastic finite fault model. The method comprises the following steps: determining an equivalent magnitude and an equivalent epicentral distance from different exceedance probability levels of a target area, determining geometric parameters of a fault model according to the equivalent magnitude, wherein the geometric parameters comprise a fault burial depth, a fault plane size and a number of sub-faults; dividing the fault plane into a plurality of sub-faults according to the number of sub-faults; calculating rupture propagation time and seismic wave propagation time of each sub-fault to an observation point by using a bilateral rupture mode; generating a seismic motion time history of each sub-fault by using the stochastic finite fault model; introducing a frequency attenuation factor and a path attenuation factor to correct the seismic motion time history; and superimposing corrected seismic motion time histories corresponding to all sub-faults to generate a seismic motion field of the target area. The technical scheme provided by the application realizes fine simulation of a fault rupture process based on a multi-probability scenario, and improves the physical rationality and spatial resolution of target area seismic motion field prediction.

The invention discloses a fault plane total stress monitoring method, which comprises the following steps of: embedding a stress monitoring device into the center of an inclined plane of a to-be-measured body, calculating six independent stress components according to strain values of the to-be-measured body in different directions acquired by the stress monitoring device, and calculating normal stress and shear stress of a fault plane according to the six independent stress components. The fault plane total stress monitoring is realized; according to the invention, the stress monitoring device with six groups of strain rosettes and the base body is embedded into the inclined plane of the object to be measured, and the full stress monitoring of the fault plane is realized through the stress values measured by the six groups of strain rosettes, so that the accurate stress change characteristics during the slippage instability of the fault plane are obtained, and a basis is provided for preventing and controlling disasters induced by the slippage instability of the fault plane.

The invention relates to a method and device for constructing a three-dimensional fracture density model of a breaking joint body and electronic equipment, and belongs to the technical field of oil reservoir exploration and development, and the method comprises the steps: obtaining a field outcrop image of a target area; marking cracks on the field outcrop image of the target area, and obtaining the density of the cracks; obtaining a scatter diagram based on the density of the crack and the first distance between the crack and the fault; fitting data on the scatter diagram to obtain a fitting function; constructing a geologic model of the target area; calculating a second distance from a grid in the geologic model to a fault main shaft surface; and based on the fitting function and the second distance, obtaining a three-dimensional model of the fracture density of the breaking joint body in the target area. According to the method, the fracture density of the breaking joint body can be depicted with higher precision, so that the recovery rate of oil and gas is improved.

The invention provides a method and system for automatically extracting and identifying related tectonic earthquakes under the constraint of a three-dimensional geometric structure of a fault, and relates to the technical field of earthquake activity analysis. According to the method, fault surface trace data, geological occurrence parameters and three-dimensional geometric parameters are obtained, the fault strike angle is determined through principal component analysis, the inclination angle is determined based on orthogonal vector decomposition, and then a three-dimensional fault plane and fault body geometric model is constructed through curved surface modeling; and under the constraint of the three-dimensional geometric structure, judging the spatial relationship between the seismic points in the seismic directory data and the fault, realizing the classification and identification of the fault-related tectonic earthquake, and outputting an analysis result. According to the invention, the spatial integrity of fault geometric expression and the tectonic earthquake identification precision are improved.

The present application relates to a geological modeling method for a geological feature uncertain reservoir. The method comprises the following steps: identifying uncertainty factors of the geological feature uncertain reservoir; parameterizing the uncertainty factors; determining possible distribution modes of the uncertainty parameters; respectively performing random sampling on each uncertainty parameter under each distribution mode to obtain a plurality of certain parameter combinations; each certain parameter combination corresponds to a certain geological model; predicting the geological feature uncertain reservoir through the certain geological model to obtain a plurality of prediction results; and statistically analyzing the prediction results to evaluate the prediction probability and prediction risk of each certain geological model. The present application provides a plurality of certain reservoir geological models varying with key geological features such as stratigraphic plane, fault plane and fracture plane, can be fitted through production history, can invert all the modeled uncertainty parameters, can reduce model uncertainty, and is more suitable for practical use.

According to the method for calculating the fortified length of the cross section expanded section of the strike-slip fault tunnel and the related equipment provided by the invention, through fitting processing of the axes of the tunnels on the two sides after fault dislocation, accurate quantitative characterization of the axis form after dislocation is realized, and basic data support fitting an actual engineering scene is provided for fortified length calculation. The section expanding sections are symmetrically arranged in an extending mode with the section face as the center, full-coverage protection of a fault dislocation influence area is guaranteed, and the problem that local protection is lost possibly due to single-side fortification is avoided. The calculation logic combining the axis fitting and the clearance requirement breaks through the limitation that an existing engineering category method depends on experience judgment, quantitative calculation of the fortification length is achieved, it is guaranteed that the anti-dislocation performance of the tunnel meets the clearance requirement, construction cost waste caused by the too large fortification length can be avoided, and the method is suitable for popularization and application. And a foundation is laid for smooth fitting of a post-earthquake tunnel axis and rapid recovery of functions.

The application relates to an artificial intelligence-based fault interpretation method, which comprises the following steps: converting three-dimensional seismic data into instantaneous phase bodies (S1); intercepting a plurality of instantaneous phase profiles perpendicular to the strike of a fault indication (S2); extracting phase points from each instantaneous phase profile at a certain phase interval, and assigning natural numbers to the points of different phases (S3); determining an interpretation point extraction distance A by using an iterative method according to the distance between two points with the same value (S4); extracting the points with the same value as the condition of meeting the distance A, and connecting the extracted points by lines to obtain a stratum interpretation line (S5); determining a breakpoint extraction distance B by using an iterative method according to the distance between the breakpoints of two different stratum interpretation lines (S6); extracting the breakpoints of each stratum interpretation line as the condition of meeting the distance B, and connecting the extracted breakpoints by lines to obtain a fault plane interpretation line (S7); and splicing the fault plane interpretation lines of each instantaneous phase profile to obtain a three-dimensional model of the fault (S8).

The present application relates to a kind of tension basin extrusion strike-slip fault spatial evolution method, belong to oil and gas field exploration and development, mineral evaluation prediction technical field.It includes the following steps: by prestack, poststack processing of seismic data, using the discontinuity of seismic data properties, using coherence algorithm contrast target point (x, y, z) and with (x, y, z) as center, all data points in the spatial range of seismic bin length radius contrast, similarity is obtained;Inheritance type fault period is distinguished, and the standard for distinguishing is that the point of each change of strike is taken as the distinguishing point;For different period, make fault attribute slice, form fault plane pattern diagram, implement strike-slip fault plane distribution characteristics, judge strike-slip fault relative displacement direction;By implementing associated structure generation period.The problem of fine identification of zone level strike-slip fault is solved, and the identification accuracy of strike-slip fracture zone fault is higher.

The invention discloses a quantitative evaluation method for the damage degree of a fault to a shale gas storage condition. The method comprises the following steps: S1, depicting fault development characteristics of a target layer series of a to-be-detected region; s2, fault levels are divided according to the fault displacement, the extension length on the fault plane and the vertical broken layer; s3, three-dimensional seismic attribute extraction is carried out on the target stratum, natural fracture sensitive attributes are selected based on geological laws and single well core fracture statistical results, and the widths W of fault fracture zones of different levels are quantitatively described; s4, on the basis of the current stress field, the fluid pressure and the fault occurrence, the effective normal stress sigman borne by the three-dimensional space interruption layer is calculated; and S5, establishing a damage index R of the fault to the shale gas storage condition according to the fault fracture zone width W and the size sigma n of the effective normal stress borne by the fault plane, and evaluating the influence of the fault on the shale gas storage condition. Through the method, different scales, different occurrence states, different burial depths and different parts of the same fault can be evaluated.

The invention relates to the technical field of earthquake three-dimensional fault plane construction, and discloses a method and device for constructing a three-dimensional fault plane based on an earthquake catalog, and the method comprises the steps: obtaining the earthquake catalog of a target region, extracting the position information of an earthquake source, and converting the position information to a UTM coordinate system; classifying the seismic source position information by using a DPC algorithm to obtain a plurality of first-order clustering results, performing PCA analysis, and extracting feature values and corresponding feature vectors; if the characteristic value meets a set condition, forming a fault plane based on the first-order clustering result; performing secondary classification on the formed seismic source points in the fault plane by using a DPC algorithm to obtain a plurality of second-order sub-clustering results, performing PCA analysis and PCA plane judgment, and if a set condition is met, obtaining a fault section based on the second-order sub-clustering results; and for the fault plane and the fault section, determining plane parameters through feature vectors obtained through PCA (Principal Component Analysis). According to the method, DPC and PCA are combined, and end-to-end parameterization-free three-dimensional fault plane modeling is realized only by means of seismic source position information.

The present invention discloses a three-dimensional fault / fracture network modeling method and device under a thrust nappe system based on the discrete element method, which relates to the technical field of three-dimensional fault modeling, and includes constructing an initial fault plane model based on the seismic logging data parameters of the target area; extracting the formation parameters of the target area, and combining with the balanced profile restoration method to determine the deformation process; calibrating the mesoscopic parameters of the discrete element according to the formation parameters, depositing the first discrete element model, optimizing and generating the second discrete element model, setting boundary conditions for thrust nappe simulation, and adjusting the model parameters through the comparison between the output results of the second discrete element model and the actual geological profile until the two match, establishing a three-dimensional fault and fracture network, analyzing and establishing the spatial distribution characteristics of the target area, and dynamically evaluating and trend predicting the activity of the three-dimensional fault. The present invention integrates a geomechanical model, overcomes the subjectivity of manual interpretation, and effectively avoids problems such as mesh reconstruction and fault tip extension.

The invention relates to an igneous rock intrusive fault zone structure seismic characterization method and system, and belongs to the technical field of oil and gas field exploration and development. The method comprises the following steps: loading geological data of a target area; constructing a fault plane of the magma intrusive target fault zone; according to the fault plane of the magma intrusive target fault zone, determining the spatial geometry and position of the fault zone in the target area; calculating seismic indication factors and distribution thereof of the fault zone in the target area according to the spatial geometric morphology and position of the fault zone in the target area; and counting the threshold value of the seismic indication factor of the igneous rock intrusive fault zone found in the exploration area, comparing the seismic indication factor of the fault zone to be evaluated with the threshold value of the seismic indication factor, and quantitatively evaluating the structure of the igneous rock intrusive fault zone according to the comparison result. Compared with a traditional SGR method and the like, the method does not need to count a large number of geological parameters, but directly carries out research on seismic data, and the implementation process is more intuitive and rapid.