51 results about "Fault plane" patented technology

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 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 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.

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 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 invention discloses a fault element constraint-based fault automatic interpretation method. The method comprises the following steps of S1, selecting fault plane sensitive attributes; s2, performing fault plane interpretation on more than two horizon positions according to fault plane sensitive attributes; and S3, projecting a fault plane interpretation result to a seismic section, and performing automatic interpretation in combination with fault elements and fault sensitive attributes. According to the automatic fault interpretation method based on the fault element constraint, the problem of rapid, accurate and automatic interpretation of the fault with the fault displacement smaller than 20 m is solved.

The present invention discloses a method, device, equipment and medium for determining the distribution of magma in an intrusive fault zone, which relates to the technical field of oil and gas field exploration and development, including: loading geological data of a target area; constructing a fault plane of a magma intrusive fault zone in the target area based on the geological data; calculating a magma distribution indicator factor of the magma intrusive fault zone in the target area; and determining the magma distribution position in the magma intrusive fault zone in the target area. The present invention considers the value characteristics and frequency characteristics of the seismic information of the magma intrusion body, and combines the distribution characteristics of the fault to propose a magma distribution indicator factor S of the igneous rock intrusive fault zone, thereby achieving the purpose of quantitatively predicting the distribution range of magma in such fault zones. It can be used to predict the closed and drained oil and gas locations in the fault zone, and the resulting prediction effect is more in line with actual geological laws. It does not require the statistics of a large number of geological parameters, and can be directly studied on seismic data, making the implementation process more intuitive and rapid.

This invention provides a method and system for determining fault planes using stratigraphic gradients and dip curvature, belonging to the field of petroleum exploration. The method includes: determining multiple stratigraphic gradients corresponding to multiple local surface models based on seismic data of a target layer, wherein the target layer includes a fault plane; obtaining cross sections of multiple local surface models based on the multiple stratigraphic gradients corresponding to the multiple local surface models and the dip and strike of each local surface model to determine the true curvature of the multiple cross sections; and determining the fault plane in the target layer based on the multiple stratigraphic gradients corresponding to the multiple local surface models and the true curvature of the multiple cross sections. This invention solves the technical problem of existing technologies being unable to accurately determine the curvature of fault planes, thus failing to accurately determine fault planes, and improves the accuracy of fault plane determination.

This invention provides a method for predicting fault connectivity based on PCA-CRITIC weights, relating to the field of oil and gas exploration technology. It includes determining the influencing factors of fault vertical connectivity, including the thickness of the displaced caprock, clay content, fault plane normal stress, formation fluid pressure, caprock brittleness index, and fault displacement. A comprehensive mathematical formula for the fault vertical connectivity coefficient is established. Based on the oil, gas, and water data from drilled wells and paleoreservoir tracing data near the displaced strata, the vertical connectivity coefficient is calculated using the comprehensive mathematical formula and by well and well depth segment. The vertical connectivity is analyzed based on the fault vertical connectivity coefficient and the oil and gas content of the near displaced strata, predicting the fault vertical connectivity in oil and gas basins. This invention accurately characterizes the oil and gas migration paths within faults, significantly improving prediction accuracy and exploration efficiency.

The invention relates to the technical field of three-dimensional modeling, and discloses a system and method for exporting a stratum section line file adaptive to MapGIS modeling based on Petrel software, and the method comprises the steps: preparing three-dimensional structure planar graph data of a stratum to be modeled; establishing a vertical fault plane in Petrel software along the trend of a target section, and converting the vertical fault plane into a Surface file; utilizing a Surface operation function to generate a space intersection line between the upright fault plane and each structural plane graph; and exporting the intersecting line data as an SHP format file, and converting the intersecting line data into a line file which can be identified by a MapGIS platform. The system comprises a corresponding data preparation module, an upright fault plane construction module, an intersection line generation module and a data export conversion module. According to the method, efficient and accurate conversion of core structure data from a Petrel platform in the field of oil exploration and development to a MapGIS platform in the field of geological survey is realized, and the problem of cross-department and cross-software data barriers is solved.

The invention relates to a fault stress plugging coefficient plane characterization rapid implementation method and device, a medium and equipment. The method includes reading a regional stress direction angle and a fault polygon orientation data file. And calculating the fault stress plugging coefficient of each fault according to the regional stress direction angle and the fault polygonal orientation data file. And forming a plane graph by gridding the fault stress plugging coefficient of each fault. And reading the corresponding fault plane polygon, and externally tangent the plane graph by using the corresponding fault plane polygon to obtain a stress plugging coefficient plane graph reflecting the form of the fault plane. And respectively obtaining stress plugging coefficient plane graphs of different faults and combining the stress plugging coefficient plane graphs to obtain a stress plugging coefficient plane graph describing all the faults. According to the method, calculation is simple, only regional stress direction information needs to be provided, drilled well information does not need to be obtained, the fault plugging performance can be directly analyzed in a well-free and few-well area and an area with fast continental facies stratum phase change, and calculation efficiency is high.

The invention provides a coal face fault modeling method and device. The method comprises the steps that upper and lower wall coordinate string information of a coal face fault is acquired; extracting an elevation value at the intersection point of the fault and the coal seam contour line; acquiring a smooth projection line equation of a fault hanging wall and hanging wall curve in a plane, wherein the smooth projection line equation is obtained by fitting the hanging wall and hanging wall coordinate string information; according to the elevation value and the smooth projection line equation, predicting a plurality of control points of the upper and lower disks; determining a space intersection endpoint of the upper and lower plate curves; performing cubic spline curve fitting on the hanging wall curve and the footwall curve by combining a control point sequence integrated by the elevation value, the plurality of control points and the space intersection end intersection point with geological constraint conditions to obtain a hanging wall line and a footwall line of the fault; and constructing a fault plane model. Therefore, through the three-dimensional modeling method for the fault in the coal seam fusing geological control points, space curve fitting and multiple constraint conditions, the consistency of the precision of the model and geological interpretation is ensured, and the method is suitable for a three-dimensional digital modeling task of a complex fault geologic body.

Embodiments of the present application provide a kind of multi-point seabed seismic generation method and related equipment, can be synthesized multi-point seismic time-domain non-stationary and the spatial variability of complex topography when considering seismic.The method comprises: according to three direction seismic waveform data, random noise and target single-point random seismic model generation one-dimensional seismic waveform data;One-dimensional seismic waveform data is spliced with random data satisfying Dirichlet distribution, to obtain spliced multi-point seismic data;According to the size of semi-space homogeneous isotropic geologic body, the position of fault plane and the size of the fault plane, establish three-dimensional site model containing fault;According to the fault plane parameters corresponding to different moment magnitude and three-dimensional site model, the simulation analysis of seismic is carried out by spectral element method, to obtain the multi-point seismic data of each measuring point in measuring point set;Spliced multi-point seismic data and multi-point seismic data are input into multi-point seismic data generation model, to obtain multi-point seismic.

The embodiment of the invention provides a fracture identification description method based on frequency and maximum likelihood attributes, and belongs to the technical field of oil reservoir development. The method comprises the following steps: processing seismic data corresponding to a target area by using a frequency division phase analysis technology to generate a plurality of frequency division data volumes; performing filtering processing on the plurality of frequency division data volumes to generate corresponding filtering volumes; calculating a maximum likelihood body according to the filtering body, and dividing the maximum likelihood body based on a fault plane direction scale and a fault inclination angle direction scale to generate a plurality of fracture data bodies; and determining a target fracture area based on the fracture data volume. According to the fracture identification description method based on the frequency and the maximum likelihood attribute provided by the embodiment of the invention, the definition of the identified fracture surface can be improved, the technical process is easy to implement in the industry, and the working efficiency can be improved.

The invention discloses a three-dimensional earthquake assessment method for carbonate rock gentle slope co-sedimentary fault activity, and particularly relates to the field of petroleum and natural gas geological exploration, and the method is characterized in that fault plane interpretation is carried out on a co-sedimentary fault, and stratigraphic plane interpretation is carried out on a fault-related stratum; a stratum used for evaluating fault activity is selected from the interpretation result, and a sequence stratigraphic framework is constructed through the stratum; the stratum deposition time of the fault activity period is obtained through calculation according to the time corresponding to the sequence stratum boundary of the sequence stratigraphic framework; selecting a plurality of seismic sections from the fault plane-stratigraphic plane interpretation result, and obtaining a fault plane trajectory and a stratigraphic plane trajectory passing through each seismic section; calculating the average fault activity rate of each seismic section in different periods by combining the stratum deposition time of the fault activity period; and evaluating the activity of the co-sedimentary fault according to the average fault activity rate of each seismic section in different periods to obtain a seismic evaluation result.

This invention discloses a three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes in carbonate rocks, specifically relating to the field of oil and gas geological exploration. Its key points are: interpreting the fault plane of the syn-sedimentary fault and interpreting the stratigraphic plane of the fault-related strata; selecting strata for assessing fault activity from the interpretation results and constructing a sequence stratigraphic framework using these strata; calculating the sedimentary time of the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework; selecting multiple seismic profiles from the fault-stratigraphic plane interpretation results and obtaining the fault plane trajectory and stratigraphic plane trajectory lines passing through each seismic profile; calculating the average fault activity rate of each seismic profile at different periods based on the sedimentary time of the fault activity period; and assessing the activity of the syn-sedimentary fault based on the average fault activity rate of each seismic profile at different periods to obtain the seismic assessment results.

This invention relates to the field of coal mining technology and discloses a Matlab-based system for calculating the critical water pressure for water inrush in non-uniform layered coal seam floor structures. The system includes the following modules: a menu bar module with "Start" and "Help" buttons; clicking the "Start" button allows for data import and export, while clicking the "Help" button displays software-related information; a parameter value module for displaying and modifying user-imported data; a range value module for user input of parameters such as the number of layers, the angle θ1 between the fracture surface and the horizontal plane, the angle β between the fault plane and the horizontal plane, the length of the goaf, and the height of the failure zone; and a display module for showing the relationship between the critical water pressure below layer i and the fracture angle θ1. This invention, through its proposed system for calculating the critical water pressure for water inrush in non-uniform layered coal seam floor structures, provides theoretical and data support for subsequent evaluation of coal seam floor water inrush and is of significant importance for the prevention and control of coal seam floor water inrush.

The present invention relates to a method and system for seismic characterization of the structure of an igneous intrusive fault zone, and belongs to the field of applied oil and gas field exploration and development technology. The method comprises loading geological data of the target area; constructing a fault plane of a magma intrusive target fault zone; determining the spatial geometry and position of the fault zone in the target area based on the fault plane of the magma intrusive target fault zone; calculating the seismic indicator factor and its distribution of the fault zone in the target area based on the spatial geometry and position of the fault zone in the target area; statistically analyzing the threshold value of the seismic indicator factor of the igneous intrusive fault zone that has been discovered in the exploration area, comparing the seismic indicator factor of the fault zone to be evaluated with the threshold value of the seismic indicator factor, and quantitatively evaluating the structure of the igneous intrusive fault zone based on the comparison result. Compared with traditional SGR and other methods, the present invention does not require statistical analysis of a large number of geological parameters, but directly conducts research on seismic data, making the implementation process more intuitive and rapid.

A method for three-dimensional velocity geological modeling with structures and velocities randomly arranged, including determining base points in three-dimensional space, building equation according to the base points to determine planar layered model, complicating a tilt layer of planar layered model, and building a fold layer model of a surface in three-dimensional space; building three-dimensional fault folded model based on the three-dimensional surface fold layer model combined with a fault plane of a random reference point and displacement of each point in a global coordinate system; building a velocity model containing a salt body based on the three-dimensional fault folded model, and simulating salt body intrusion in a geological body of a certain depth; and performing a random velocity amplitude to realize three-dimensional velocity modeling according to the layered type which has been set and according to the set velocity range and the velocity difference range between each layer of geology.

The application provides a fault two-plate sand body connectivity evaluation method based on seismic data, comprising the following steps: step 1, carrying out structural interpretation of a fault plane based on three-dimensional seismic data to form a fault plane planar structure map; step 2, using three-dimensional seismic data to perform layer flattening with the fault plane as a target layer, and establishing a layer flattening seismic data body; step 3, using the layer flattening seismic data body to extract attributes of the upper and lower plates of the fault respectively, and calculating a difference ratio A of the attributes on both sides of the fault; step 4, using the A value of a known open fault to determine the A value range of the sand body connectivity of the two plates of the fault, and defining the A value range as a threshold B, and determining the sand body connectivity of the two plates of the whole fault according to the B and the attribute graph A. The method is based on three-dimensional seismic data, fully utilizes logging data and is relatively simple to apply. The fault two-plate sand body layering connectivity can be effectively carried out in a three-dimensional seismic data coverage area, and a new method is provided for evaluating the feasibility of layering water injection between fault block reservoirs.

The invention relates to the technical field of coal seam transformation engineering, in particular to a complex geologic structure coal seam directional hydraulic slotting and controllable acidification composite transformation method which comprises the steps that the in-situ crustal stress direction of a target coal seam section and the spatial orientation of a fault surface and a fold wing part are obtained, and the arrangement orientation of a directional hydraulic slotting is determined according to the in-situ crustal stress direction and the spatial orientation of the fault surface and the fold wing part; then directional hydraulic slotting is carried out in the drill hole to form initial directional slots which are distributed along the depth of the drill hole in a multi-cluster mode and have different directions, and the initial directional slots are used for limiting the crack initiation position and the crack initiation direction of the crack; on the basis, a controllable acid liquid system is injected into the coal seam section, so that the acid liquid enters the initial directional kerfs and the natural fractures of the coal seam in a low-viscosity state, the pH value of the system is changed by reacting with minerals in coal in the advancing process along the fractures, a gelling unit is triggered to form an in-situ temporary plugging area, and the follow-up acid liquid is guided to flow along different fracture channels; and after the fracture acidizing treatment is completed by continuously injecting acid liquor, flowback or replacement operation is carried out on the coal seam section.

The application provides a palaeogeomorphology recovery method and device, which comprises the following steps: constructing an initial palaeogeomorphology; determining a compaction coefficient based on porosity data and burial depth, and using the compaction coefficient to correct the residual thickness of the initial palaeogeomorphology to obtain a decompacted palaeothickness; determining an absolute palaeowater depth based on uranium element content, and subtracting the absolute palaeowater depth from the decompacted palaeothickness to obtain a palaeotopography of a sediment surface; identifying an unconformity surface, determining an uneroded area and an eroded area according to the unconformity surface, constructing a trend model using the stratigraphic thickness variation trend of the uneroded area, determining the eroded thickness of the eroded area based on the trend model, and obtaining a complete palaeogeomorphology after erosion recovery; taking a selected reference layer as a datum plane, reversely sliding an upper wall of a fault along a fault plane until the upper wall reference layer and a lower wall reference layer are spliced, performing trend interpolation filling on data gaps at the fault to obtain a target palaeogeomorphology after the influence of the fault is removed. The application improves the accuracy and objectivity of palaeogeomorphology recovery.

The invention provides a seismic data fault detection method and device, computer equipment and a storage medium, and the method comprises the steps: obtaining an attribute data body of seismic data, the attribute data body being capable of describing fault features of the seismic data; determining a fault plane direction according to the attribute data volume; establishing a fault coordinate system according to the fault plane direction, and determining vertical parameters and parallel parameters related to the fault plane direction based on the fault coordinate system; based on the vertical parameters and the parallel parameters, utilizing a Laplace Gaussian filtering operator to carry out filtering calculation on the attribute data volume so as to obtain enhanced fault features; the filtering calculation comprises the steps of performing sharpening processing on the attribute data volume in the direction perpendicular to the fault plane and performing smoothing processing on the attribute data volume in the direction parallel to the fault plane. On the basis of a fault enhancement technology of an improved Laplace Gaussian filter operator (LoG), the direction perpendicular to a fault plane is sharpened, the direction parallel to the fault plane is smoother, local non-planar anomaly is smoother, and the tomography quality is improved.