Segmented quantitative division method for strike-slip fault zone and computer equipment

Abstract

The invention relates to a strike-slip fault zone segmentation quantitative division method and device and computer equipment. The method comprises the following steps: creating a target horizon construction contour map according to seismic data and logging data of a target area, and converting the target horizon construction contour map into a digital elevation model (DEM) data map of the target horizon; based on the seismic data, obtaining a fault probability body slice of the target horizon, and determining a strip width and a strip baseline of a strike-slip fault zone of the target horizon; according to a digital elevation model (DEM) data graph of the target horizon, and the strip width and the strip baseline of the strike-slip fault zone, obtaining an elevation-height difference distribution strip graph of the target horizon, drawing a segment difference division baseline, and marking a starting point and an ending point of each translation segment; and performing quantitative fine division and marking on the sections of the strike-slip fault zone by adopting a section difference division scheme. Therefore, fine segmentation division of the strike-slip fault zone in the target horizon can be completed quantitatively, effectively and accurately.

Description

Technical Field

[0001] The present application relates to the technical field of oil and gas geophysical exploration and development, and in particular to a method for quantitative segmentation of strike-slip fault zones, a computer device, and a computer-readable storage medium. Background Art

[0002] Segmental characterization of strike-slip fault zones plays a crucial role in oil and gas exploration and development. Segmental characterization of strike-slip fault zones can reveal geometric changes during their evolution. By analyzing the standard deviation of angles between adjacent faults within a strike-slip fault zone, the extent of strike-slip fault evolution can be quantitatively characterized. Strike-slip fault zones exhibit a certain segmentation along the strike-slip direction, including compressional uplift, tensional depression, and translational segments. These segments often alternate and exhibit distinct characteristics. This segmentation reflects the changes in the geometrical appearance and internal fracture structure of strike-slip faults under varying stresses, which are controlled by strike-slip stress. The mechanism of strike-slip fault segmentation is closely related to their stress state. In some regions, compressive stress decay and a shift in triaxial stress states are the primary drivers of strike-slip fault segmentation. This stress change leads to the formation of compressional and tensional zones in different directions. The segmentation of strike-slip fault zones is of great significance for oil and gas exploration. Analyzing the segmentation characteristics of deep strike-slip fault zones is a key issue in understanding their reservoir control. The segmented differential activity characteristics of strike-slip fault zones not only affect the distribution of oil and gas, but also determine the degree of oil and gas enrichment, resulting in significant differences in the degree of fault-karst development and the scale of oil reservoirs in different fault segments. For example, in the study of the Shunbei strike-slip fault zone, the oil and gas enrichment levels in the uplift and depression segments were higher, while those in the translation segments were lower. It is speculated that the scale and activity strength of the branch faults control their oil and gas enrichment levels. In short, the segmented characterization of strike-slip fault zones not only helps to understand their evolutionary process and stress distribution, but also has important guiding significance for oil and gas exploration and geological research. By conducting in-depth analysis of the segmented characteristics of strike-slip fault zones, we can better reveal their impact on oil and gas distribution, thereby providing a scientific basis for oil and gas exploration.

[0003] In the prior art, methods for segmenting strike-slip fault zones include geological feature analysis, seismicity analysis, and numerical simulation. The geological feature analysis method observes and analyzes the geological features of the strike-slip fault zone (such as rock type, structural characteristics, and fracture surface morphology). For example, strike-slip fault zones are segmented by counting the fault drop of the main fault. This method simplifies the strike-slip fault zone into a single main fault and uses seismic slices perpendicular to the main fault to calculate the fault drop of the main fault, ignoring the influence of secondary faults and resulting in significant statistical errors. Furthermore, the profiles perpendicular to the main fault selected in the statistical method are subjective, and the density of the profiles can significantly affect the final segmentation scheme. Seismicity analysis is based on seismic data, analyzing the seismicity of different regions (such as magnitude, focal depth, and epicenter distribution) and dividing areas with significantly different activity into different segments. This method requires long-term seismic monitoring data and is not accurate in areas where data is scarce, resulting in only large-scale preliminary segmentation. Numerical simulation uses numerical simulation techniques to simulate the behavior and stress state of a fault, and then segment the fault based on the results. This method relies on model construction, and the accuracy of the model depends on the selection of input parameters, and is computationally intensive. Overall, conventional methods for segmenting strike-slip fault zones suffer from the limitations of relying on the researcher's experience and judgment, resulting in high subjectivity, long timescales, and high computational complexity. Consequently, they are unable to provide a quantitative and precise segmentation of strike-slip fault zones, and therefore, cannot form a unified quantitative criterion.

[0004] Strike-slip fault zones are often complicated by a series of secondary faults. With the continuous deepening of oil and gas exploration and development, the requirements for the refinement of exploration results are becoming increasingly higher, which means that there are also higher requirements for the segmentation of strike-slip fault zones. How to achieve the precise quantitative division of each segment of the strike-slip fault zone is an urgent problem that needs to be solved, which affects the subsequent precise exploration, development and utilization of oil and gas resources. Summary of the Invention

[0005] Based on this, it is necessary to provide a strike-slip fault zone segmentation quantitative division method, computer equipment and computer-readable storage medium that can achieve fine quantitative division of each segment of the strike-slip fault zone in order to address the above technical problems.

[0006] A method for quantitatively dividing a strike-slip fault zone into sections, the method comprising:

[0007] Create a structural contour map of the target layer based on the seismic data and well logging data of the target area;

[0008] Performing grid analysis on the structural contour map of the target horizon, and converting the structural contour map of the target horizon into a digital elevation model (DEM) data map of the target horizon;

[0009] Based on the seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability volume analysis to obtain a fault probability volume slice of the target horizon;

[0010] Determining the strip width and strip baseline of the strike-slip fault zone of the target horizon according to the fault probability volume slice of the target horizon;

[0011] Performing data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone to obtain an elevation-elevation distribution strip map of the target layer;

[0012] According to the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmentation formula, a segmentation difference division baseline is drawn in the elevation-elevation difference distribution strip diagram of the target layer to obtain a new elevation-elevation difference distribution strip diagram;

[0013] Determine the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer according to the average elevation strip line in the new elevation-elevation difference distribution strip diagram and the segmented difference division baseline, and mark the starting point and the ending point of each translation segment in the new elevation-elevation difference distribution strip diagram to obtain a marked elevation-elevation difference distribution strip diagram;

[0014] The digital elevation model (DEM) data map of the target layer is overlapped with the fault probability volume slice to form an elevation fault superposition map of the target layer. According to the marked elevation-elevation difference distribution strip map, a segmented difference division scheme is adopted to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

[0015] In one embodiment, the fault probability volume analysis is performed based on the seismic data using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume slice of the target horizon, including:

[0016] performing normalization processing on the seismic data to remove random noise and obtain denoised seismic data;

[0017] Performing fault probability volume analysis on the denoised seismic data using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume;

[0018] performing refinement and enhancement processing on the fault probability volume to obtain an enhanced fault probability volume;

[0019] Extracting a fault probability volume slice of the target horizon along the target horizon from the enhanced fault probability volume.

[0020] In one embodiment, determining the strip width and strip baseline of the strike-slip fault zone in the target horizon based on the fault probability volume slice of the target horizon includes:

[0021] Determining the strike direction of the strike-slip fault zone of the target horizon according to the fault probability volume slice of the target horizon;

[0022] The range covering the width of the strike-slip fault zone is taken as the strip width of the strike-slip fault zone of the target horizon;

[0023] A vector profile line parallel to the strike-slip fault zone is drawn through the symmetry center of the strip width of the strike-slip fault zone of the target horizon as the strip baseline of the strike-slip fault zone of the target horizon.

[0024] In one embodiment, the step of drawing a segmented difference division baseline in the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmentation formula to obtain a new elevation-elevation difference distribution strip diagram includes:

[0025] Obtaining the elevation values of the intersection points of the average elevation strip line and the elevation difference strip line in the elevation-elevation distribution strip diagram of the target layer, and determining the minimum elevation value and the second minimum elevation value from the elevation values of each intersection point;

[0026] According to the minimum elevation value and the second minimum elevation value, a quantitative difference segmentation formula is used to determine the elevation value of the strike-slip fault zone segmentation difference baseline of the target layer;

[0027] According to the elevation value of the segmented difference division baseline of the strike-slip fault zone of the target layer, the segmented difference division baseline of the strike-slip fault zone of the target layer is drawn in the elevation-elevation difference distribution strip diagram of the target layer to obtain a new elevation-elevation difference distribution strip diagram.

[0028] In one embodiment, the quantitative difference segmentation formula is:

[0029]

[0030] Among them, φ H Elevation values of the baseline divided for segment differences; is the minimum elevation value, The second smallest elevation value.

[0031] In one embodiment, the method comprises dividing the baseline of the strike-slip fault zone of the target layer according to the average elevation strip line in the new elevation-elevation difference distribution strip map and the segmented difference, determining the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer, marking the starting point and the ending point of each translation segment in the new elevation-elevation difference distribution strip map, and obtaining the marked elevation-elevation difference distribution strip map, including:

[0032] Draw a tangent line for each point of the average elevation strip line in the new elevation-elevation distribution strip diagram, and calculate the angle between the tangent line of each point and the segment difference division baseline, and determine the point where the angle is less than 20°;

[0033] Traverse the points with an angle less than 20° from left to right along the length direction of the strike-slip fault zone of the target horizon, and determine whether the current point is the first point with an angle less than 20°. If it is the first point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point is traversed, and the next point with an angle less than 20° is taken as the current point, where i∈1, 2, ..., I;

[0034] If it is not the first point with an angle less than 20°, determine whether the current point is continuous with the previous point with an angle less than 20°. If it is not continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point is traversed, and the next point with an angle less than 20° is taken as the current point;

[0035] If it is continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The end point of the segment is traversed, and the next point with an angle less than 20° is traversed as the current point. If it is continuous with the next point with an angle less than 20°, the next point with an angle less than 20° is traversed as the current point until all points with an angle less than 20° are traversed. The continuous segment set X is obtained, X = {X1, X2, ..., X I}, I is the number of continuous segments in the continuous segment set;

[0036] Eliminate the continuous segments whose length is less than a preset threshold value from the continuous segment set to obtain a continuous segment set after elimination;

[0037] The starting point and the ending point of each continuous segment in the eliminated continuous segment set are used as the starting point and the ending point of each translation segment of the strike-slip fault zone of the target layer, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation difference distribution strip diagram to obtain the marked elevation-elevation difference distribution strip diagram.

[0038] In one embodiment, the segment difference division scheme is:

[0039] The area below the segmentation baseline in the average elevation strip line in the marked elevation-elevation distribution strip diagram is divided into a strike-slip fault zone level I depression segment;

[0040] The area above the segmentation baseline in the average elevation strip line in the marked elevation-elevation distribution strip diagram is divided into a strike-slip fault zone level I uplift segment;

[0041] According to the starting point and the ending point of the translation segment in the marked elevation-elevation distribution strip diagram, the area between the starting point and the ending point of each translation segment is divided into a translation segment, and the area where each translation segment is located is further subdivided into a translation segment developed within a level I tension depression segment or a translation segment developed within a level I pressure uplift segment;

[0042] In the case that there is a translation segment in the area where the level I pull-down segment is located, the segments adjacent to each translation segment in the area where the level I pull-down segment is located and lower than the translation segment are finely divided into level II pull-down segments;

[0043] When there is a translation segment in the area where the level I pressure uplift segment is located, the segments adjacent to each translation segment in the area where the level I pressure uplift segment is located and which are higher than the translation segment are finely divided into level II pressure uplift segments.

[0044] A device for quantitatively dividing a strike-slip fault zone into sections, comprising:

[0045] The contour map creation module is used to create a structural contour map of the target layer based on the seismic data and well logging data of the target area;

[0046] A DEM data map conversion module is used to perform grid analysis on the structural contour map of the target layer and convert the structural contour map of the target layer into a digital elevation model DEM data map of the target layer;

[0047] A fault probability volume analysis module is used to perform fault probability volume analysis based on the seismic data using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume slice of the target horizon;

[0048] a strip baseline determination module, configured to determine a strip width and a strip baseline of a strike-slip fault zone of the target horizon according to a fault probability volume slice of the target horizon;

[0049] A data analysis and statistics module is used to perform data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and obtain an elevation-elevation distribution strip map of the target layer;

[0050] A segmented difference division baseline drawing module is used to draw a segmented difference division baseline in the elevation-elevation difference distribution strip map of the target layer according to the elevation-elevation difference distribution strip map of the target layer and the quantitative difference segmentation formula to obtain a new elevation-elevation difference distribution strip map;

[0051] a translation segment marking module, configured to divide the baseline according to the average elevation strip line and the segmented difference in the new elevation-elevation distribution strip diagram, determine the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer, and mark the starting point and the ending point of each translation segment in the new elevation-elevation distribution strip diagram to obtain a marked elevation-elevation distribution strip diagram;

[0052] A segmentation module is used to overlap the digital elevation model (DEM) data map of the target layer with the fault probability volume slice to form an elevation fault superposition map of the target layer. According to the marked elevation-elevation difference distribution strip map, a segmentation difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

[0053] A computer device includes a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the method for quantitative segmentation of a strike-slip fault zone are implemented.

[0054] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the method for quantitative segmentation of a strike-slip fault zone.

[0055] The above-mentioned strike-slip fault zone segmentation quantitative division method, computer equipment and computer-readable storage medium create a target layer structural contour map based on the seismic data and well logging data of the target area; perform grid analysis on the target layer structural contour map, and convert the target layer structural contour map into a digital elevation model DEM data map of the target layer; based on the seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability body analysis to obtain the fault probability body slice of the target layer; according to the fault probability body slice of the target layer, the strip width and strip baseline of the strike-slip fault zone of the target layer are determined; data analysis and statistics are performed based on the digital elevation model DEM data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and the elevation-elevation difference distribution strip map of the target layer is obtained; according to the target layer The target layer's elevation-elevation distribution strip map and quantitative difference segmentation formula are used to draw the segmentation difference division baseline in the elevation-elevation distribution strip map of the target layer to obtain a new elevation-elevation distribution strip map; according to the average elevation strip line and the segmentation difference division baseline in the new elevation-elevation distribution strip map, the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer are determined, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation distribution strip map to obtain a marked elevation-elevation distribution strip map; the digital elevation model DEM data map of the target layer is overlapped with the fault probability volume slice to form an elevation fault superposition map of the target layer; according to the marked elevation-elevation distribution strip map, the segmentation difference division scheme is used to perform quantitative fine division and marking of the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer. Therefore, the fine segmentation of the strike-slip fault zone in the target layer can be completed quantitatively, effectively and accurately, the influence of secondary faults is taken into account, the subjectivity of the existing methods is avoided, the accuracy of the strike-slip fault zone segmentation is improved, and a quantitative and unified segmentation standard is formed. It has great application advantages for work areas with developed strike-slip fault zones and can effectively guide the fine description of strike-slip fault zones and the practice of oil and gas exploration and development. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 A schematic flow chart of a method for quantitatively dividing a strike-slip fault zone into sections in one embodiment;

[0057] Figure 2 is an iso-T0 map of the target horizon T74 in one embodiment;

[0058] Figure 3 This is a structural contour map of the target horizon T74 in one embodiment;

[0059] Figure 4 A digital elevation model DEM data diagram of a target layer T74 in one embodiment;

[0060] Figure 5Schematic diagram of a fault probability volume slice extracted along the target horizon T74 in one embodiment;

[0061] Figure 6 Schematic diagram of a slice of a fault probability volume extracted along the target horizon T74 after refinement and enhancement in one embodiment;

[0062] Figure 7 Schematic diagram of a strip baseline of a strike-slip fault zone of target horizon T74 in one embodiment;

[0063] Figure 8 1 is an elevation-elevation distribution strip diagram of the target layer T74 in one embodiment;

[0064] Figure 9 A new elevation-elevation distribution strip diagram of the target layer T74 in one embodiment;

[0065] Figure 10 This is an elevation fault superposition map of the target horizon T74 in one embodiment;

[0066] Figure 11 Schematic diagram of the final division result of the target layer T74 in one embodiment. DETAILED DESCRIPTION

[0067] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.

[0068] In one embodiment, Figure 1 As shown, a method for quantitative segmentation of strike-slip fault zones is provided. The method is described by taking the application of the method to a terminal as an example, and includes the following steps:

[0069] Step S220: creating a target layer structural contour map based on the seismic data and well logging data of the target area.

[0070] Among them, the target area is the strike-slip fault zone with uplift segments and depression segments.

[0071] Among them, seismic data is a three-dimensional seismic data volume obtained by seismic processing of original seismic data collected through seismic exploration.

[0072] Among them, logging data refers to drilling information including well head, well trajectory, logging curve, and well layer data.

[0073] In one embodiment, creating a target horizon structural contour map based on seismic data and well logging data of a target area includes:

[0074] Based on the seismic data and well logging data of the target area, the location of the fault development in the target area is determined; based on the location of the fault development in the target area, the synthetic seismic record well-seismic fine calibration is carried out in combination with the seismic data and well logging data to determine the position and time-depth relationship of the target layer on the seismic profile; based on the position of the target layer on the seismic profile, the target layer is tracked and interpreted, and then the target layer is interpolated into 1×1, and an iso-T0 map of the target layer is created; based on the time-depth relationship of the target layer, the average velocity method is used for variable velocity mapping, and the iso-T0 map of the target layer is converted into a structural contour map of the target layer.

[0075] Among them, by collecting seismic, well logging and other data, and preprocessing the original well logging data including denoising, outlier removal, and normalization, the position and time-depth relationship of the target layer on the seismic profile are determined through fine calibration of synthetic seismic records according to the location of fault development, and the tracking interpretation work of the target layer is completed (the interpretation accuracy reaches at least 4×4). Then the layer is interpolated into 1×1, and an iso-T0 map of the target layer is created; using the time-depth relationship obtained by fine calibration of well logging, the average velocity method is used for variable velocity mapping to convert the iso-T0 map of the target layer into a structural contour map of the target layer.

[0076] Step S240 , performing grid analysis on the structural contour map of the target layer, and converting the structural contour map of the target layer into a digital elevation model (DEM) data map of the target layer.

[0077] In one embodiment, a target layer structural contour map is subjected to grid analysis, and the target layer structural contour map is converted into a digital elevation model (DEM) data map, including: performing grid analysis on the target layer structural contour map to determine the topological relationship between nodes, triangle edges, and triangle faces that express terrain scatter points; setting a projection coordinate system for target area matching, and converting the target layer structural contour map into a surface TIN data structure using a triangulation algorithm based on the topological relationship; rasterizing the surface TIN data structure, setting the matching output data type, sampling distance, and Z factor, and converting the surface TIN data structure into a digital elevation model (DEM) data map.

[0078] Step S260: Based on the seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability volume analysis to obtain a fault probability volume slice of the target horizon.

[0079] Among them, the hybrid neural network model formed by Resnet (residual network) and U-net in Geoplat software can be used to perform fault probability volume analysis to obtain the fault probability volume slice of the target layer.

[0080] In one embodiment, based on seismic data, a hybrid neural network model formed by ResNet and U-Net is used to perform fault probability volume analysis to obtain a fault probability volume slice of the target horizon, including:

[0081] The seismic data is regularized to remove random noise and obtain denoised seismic data; the denoised seismic data is subjected to fault probability volume analysis using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume; the fault probability volume is refined and enhanced to obtain an enhanced fault probability volume; and the enhanced fault probability volume is used to extract fault probability volume slices of the target layer along the target layer.

[0082] Among them, regularizing the seismic data to remove random noise can improve the signal-to-noise ratio of the seismic data and thus improve the accuracy of analyzing the fault probability body.

[0083] Among them, the refinement and enhancement of the fault probability body can enhance the identification accuracy and continuity of secondary faults and the intersection relationship between secondary faults and main faults.

[0084] Among them, the quality inspection can also be performed on the fault probability body obtained by performing fault probability body analysis on the hybrid neural network model formed by Resnet and U-net.

[0085] Among them, the quality of the fault probability body can be inspected according to the fault structure style, fault combination relationship, secondary fault identification characteristics, etc. When the fault identification effect is not good, the hybrid neural network model can be transferred and learned until the optimal fault probability body is obtained.

[0086] Step S280: determining the strip width and strip baseline of the strike-slip fault zone in the target horizon based on the fault probability volume slice of the target horizon.

[0087] In one embodiment, determining the strip width and strip baseline of the strike-slip fault zone in the target horizon based on the fault probability volume slice of the target horizon includes:

[0088] According to the fault probability volume slice of the target layer, the strike of the strike-slip fault zone of the target layer is determined; the range covering the width of the strike-slip fault zone is used as the strip width of the strike-slip fault zone of the target layer; and a vector profile line parallel to the strike-slip fault zone is drawn through the symmetry center of the strip width of the strike-slip fault zone of the target layer as the strip baseline of the strike-slip fault zone of the target layer.

[0089] Step S300 , performing data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and obtaining an elevation-elevation distribution strip map of the target layer.

[0090] Among them, the digital elevation model DEM data map converted from the linear map of the target layer structure, the strip width and strip baseline of the strike-slip fault zone are used for data analysis and statistics, and the appropriate elevation point step size and the number of parallel profiles in each strip are set to determine the elevation-elevation difference distribution of the strike-slip fault zone of the target layer, and obtain the elevation-elevation difference distribution strip map. The maximum elevation value, minimum elevation value, average elevation value, upper and lower quartile elevation values and elevation difference values of the strip profile area are clarified, and corresponding strip lines are drawn on the elevation-elevation difference distribution strip map, namely the maximum elevation strip line, minimum elevation strip line, average elevation strip line, upper and lower quartile elevation strip line and elevation difference strip line.

[0091] It should be understood that each strip line is connected by a number of sampling points. For example, each sampling point on the average elevation strip line represents the average elevation value of a certain length of the strike-slip fault zone.

[0092] The interval between each sampling point can be set according to actual conditions, such as 5 meters, 10 meters, etc.

[0093] Step S320 , based on the elevation-elevation distribution strip diagram of the target layer and the quantitative difference segmentation formula, a segmentation difference division baseline is drawn in the elevation-elevation distribution strip diagram of the target layer to obtain a new elevation-elevation distribution strip diagram.

[0094] In one embodiment, based on the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmentation formula, a segmentation difference segmentation baseline is drawn in the elevation-elevation difference distribution strip diagram of the target layer to obtain a new elevation-elevation difference distribution strip diagram, including:

[0095] Obtain the elevation values of the intersection points of the average elevation strip line and the elevation difference strip line in the elevation-elevation difference distribution strip map of the target layer, and determine the minimum elevation value and the second minimum elevation value from the elevation values of each intersection point; based on the minimum elevation value and the second minimum elevation value, use the quantitative difference segmentation formula to determine the elevation value of the strike-slip fault zone segmentation difference division baseline of the target layer; based on the elevation value of the strike-slip fault zone segmentation difference division baseline of the target layer, draw the strike-slip fault zone segmentation difference division baseline of the target layer in the elevation-elevation difference distribution strip map of the target layer to obtain a new elevation-elevation difference distribution strip map.

[0096] The average elevation strip line is a strip line formed by connecting the average elevation values of all sampling points.

[0097] Among them, the height difference strip line is a strip line formed by connecting the height difference values of all sampling points, where a certain height difference value is the difference between the maximum elevation value and the minimum elevation value on a certain length of the strike-slip fault zone (that is, the length corresponding to a certain sampling point on the strike-slip fault zone).

[0098] In one embodiment, the quantitative difference segmentation formula is:

[0099]

[0100] Among them, φ H Elevation values of the baseline divided for segment differences; is the minimum elevation value, The second smallest elevation value.

[0101] Step S340, divide the baseline according to the average elevation strip line and the segmented difference in the new elevation-elevation difference distribution strip diagram, determine the starting point and the end point of the translation segment of the strike-slip fault zone of the target layer, and mark the starting point and the end point of each translation segment in the new elevation-elevation difference distribution strip diagram to obtain the marked elevation-elevation difference distribution strip diagram.

[0102] In one embodiment, the baseline is divided according to the average elevation strip line and the segmented difference in the new elevation-elevation distribution strip diagram, the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer are determined, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation distribution strip diagram to obtain the marked elevation-elevation distribution strip diagram, including:

[0103] Draw a tangent to each point of the average elevation strip line in the new elevation-elevation distribution strip map, and calculate the angle between the tangent of each point and the segment difference division baseline, determine the point with an angle less than 20°, traverse the points with an angle less than 20° from left to right along the length direction of the strike-slip fault zone of the target layer, and determine whether the current point is the first point with an angle less than 20°. If it is the first point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as the continuous segment X. i The starting point of the segment X is traversed, and the next point with an angle less than 20° is taken as the current point, where i∈1, 2, ..., I; if it is not the first point with an angle less than 20°, determine whether the current point is continuous with the previous point with an angle less than 20°. If it is not continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point of the segment X is traversed to the next point with an angle less than 20° as the current point. If it is continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X.i The end point of the segment is traversed, and the next point with an angle less than 20° is traversed as the current point. If it is continuous with the next point with an angle less than 20°, the next point with an angle less than 20° is traversed as the current point until all points with an angle less than 20° are traversed. The continuous segment set X is obtained, X = {X1, X2, ..., X I}, I is the number of continuous segments in the continuous segment set; the continuous segments whose length is less than the preset threshold are eliminated from the continuous segment set to obtain the eliminated continuous segment set; the starting point and the ending point of each continuous segment in the eliminated continuous segment set are used as the starting point and the ending point of each translation segment of the strike-slip fault zone of the target layer, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation difference distribution strip diagram to obtain the marked elevation-elevation difference distribution strip diagram.

[0104] Among them, each point of the average elevation strip line refers to a sampling point.

[0105] The preset threshold value can be determined based on the actual situation and the strip width of the strike-slip fault zone in the target layer, such as one-fifth, one-quarter, etc. of the strip width.

[0106] Step S360: Overlap the digital elevation model (DEM) data map of the target layer with the fault probability volume slice to form an elevation fault superposition map of the target layer. Based on the marked elevation-elevation difference distribution strip map, a segmented difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

[0107] It should be understood that the marked elevation-elevation distribution strip map is drawn based on the digital elevation model DEM data map. The elevation fault overlay map of the target layer has a corresponding positional relationship with the strike-slip fault zone represented by the marked elevation-elevation distribution strip map. The characteristics of the average elevation strip line in the marked elevation-elevation distribution strip map are used to divide the various segmented areas of the strike-slip fault zone in the elevation fault overlay map and mark them.

[0108] In one embodiment, the segment difference partitioning scheme is:

[0109] The area below the segmentation baseline in the average elevation strip line of the marked elevation-elevation distribution strip map is divided into the level I tensile depression section of the strike-slip fault zone;

[0110] The area above the segmentation baseline in the average elevation strip line of the marked elevation-elevation distribution strip map is divided into the level I uplift section of the strike-slip fault zone;

[0111] According to the starting and ending points of the translation segments in the marked elevation-elevation distribution strip diagram, the area between the starting and ending points of each translation segment is divided into translation segments, and the area where each translation segment is located is further subdivided into translation segments developed within the first-level tension depression segment or translation segments developed within the first-level pressure uplift segment.

[0112] In the case that there is a translation segment in the area where the level I pull-down segment is located, the segments adjacent to each translation segment in the area where the level I pull-down segment is located and lower than the translation segment are finely divided into level II pull-down segments;

[0113] When there is a translation segment in the area where the level I pressure uplift segment is located, the segments adjacent to each translation segment in the area where the level I pressure uplift segment is located and which are higher than the translation segment are finely divided into level II pressure uplift segments.

[0114] The above-mentioned strike-slip fault zone segmentation quantitative division method creates a target layer structural contour map based on the seismic data and well logging data of the target area; performs grid analysis on the target layer structural contour map, and converts the target layer structural contour map into a digital elevation model (DEM) data map of the target layer; based on the seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability body analysis to obtain the fault probability body slice of the target layer; according to the fault probability body slice of the target layer, the strip width and strip baseline of the strike-slip fault zone of the target layer are determined; data analysis and statistics are performed based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and the elevation-elevation difference distribution strip map of the target layer is obtained; according to the elevation-elevation difference of the target layer, the strip width and strip baseline of the strike-slip fault zone are obtained; A strip map and a quantitative difference segmentation formula are used to draw a segmentation difference division baseline in the elevation-elevation difference distribution strip map of the target layer to obtain a new elevation-elevation difference distribution strip map; according to the average elevation strip line and the segmentation difference division baseline in the new elevation-elevation difference distribution strip map, the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer are determined, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation difference distribution strip map to obtain a marked elevation-elevation difference distribution strip map; the digital elevation model DEM data map of the target layer is overlapped with the fault probability volume slice to form an elevation fault superposition map of the target layer; according to the marked elevation-elevation difference distribution strip map, the segmentation difference division scheme is used to perform quantitative fine division and marking of the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer. In this way, the fine segmentation of the strike-slip fault zone in the target layer can be completed quantitatively, effectively and accurately, the influence of secondary faults is taken into account, the subjectivity of the existing methods is avoided, the accuracy of the strike-slip fault zone segmentation is improved, and a quantitative and unified segmentation standard is formed. It has great application advantages for work areas with developed strike-slip fault zones and can effectively guide the fine description of strike-slip fault zones and the practice of oil and gas exploration and development.

[0115] In one embodiment, taking a typical strike-slip fault zone development area in a certain oil field as an example, the strike-slip fault zone is quantitatively divided into sections using the above-mentioned strike-slip fault zone segmentation method:

[0116] (1) Create a structural contour map of the target layer based on the seismic data and well logging data of the target area:

[0117] Collect seismic, well logging and other data in the target area, and perform denoising and outlier removal on the original seismic data; perform outlier removal and standardization on the well logging curves (such as acoustic time difference (AC), density (DEN), longitudinal wave impedance (P-imp), SP natural potential, etc.); determine the position of the target layer on the seismic profile through fine calibration of synthetic seismic records according to the location of fault development. The target layer selected this time is T74 for example. Complete the tracking interpretation work of the target layer T74 with an interpretation accuracy of 4×4, then interpolate the T74 layer to 1×1, and create the following Figure 2 The iso-T0 map of the target layer T74 is shown in the figure. The time-depth relationship obtained by fine calibration of well seismic data is used to create a variable velocity map using the average velocity method. Figure 3 The structural contour map of the target layer T74 shown can be realized through SMI software.

[0118] (2) Perform grid analysis on the structural contour map of the target layer and convert the structural contour map of the target layer into a digital elevation model (DEM) data map of the target layer:

[0119] The target layer T74 structural contour map is grid analyzed to determine the relationship between nodes, triangle edges and triangle faces, thereby displaying or implicitly expressing the topological relationship of the terrain scattered points. The structural contour map of the target layer T74 is then converted into a surface TIN data structure using a triangulation algorithm. According to the location of the target area, the projection coordinate system matched this time is WGS_1984_UTM_Zone_43N. The target layer T74 surface TIN data structure is rasterized. According to the T74 surface TIN data structure characteristics, the output data type is set to floating point, the sampling distance is set to 5m, and the Z factor is set to 1. Finally, it is converted into the following: Figure 4 The digital elevation model DEM data map shown can be realized through ArcGis software.

[0120] (3) Based on seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability volume analysis to obtain the fault probability volume slice of the target layer:

[0121] Based on the hybrid neural network model formed by Resnet (residual network) and U-net in Geoplat software, the original seismic data is regularized to remove random noise, improve the signal-to-noise ratio of the seismic data volume, and improve the recognition accuracy of faults; the hybrid neural network model is used to calculate the fault probability volume and extract the T74 slice along the target layer in the fault probability volume. Figure 5 As shown in Figure 1, it is found that the fault continuity shown on the fault probability volume along the T74 slice is poor, and some secondary faults cannot be identified. Then the fault probability volume is refined and enhanced to obtain the enhanced fault probability volume. In the enhanced fault probability volume, the following is extracted along the target layer T74: Figure 6 The fault probability volume slice shown in the figure shows that the continuity of the fault has improved, the accuracy of secondary fault identification has increased, and the intersection relationship between secondary faults and main faults has become more obvious.

[0122] (4) According to the fault probability volume slice of the target horizon, the strip width and strip baseline of the strike-slip fault zone of the target horizon are determined:

[0123] According to the probability volume slice of the target layer T74 fault, it can be determined that the strike-slip fault zone is distributed in the NNW-SSE direction. The range that can cover the width of the strike-slip fault zone is taken as the strip width of the strike-slip fault zone. It can be determined that the strip width of the strike-slip fault zone developed in the target area is 1200m. Then, a vector profile line parallel to the strike-slip fault zone is drawn through the symmetry center of the strike-slip fault zone, which is the strip baseline of the strike-slip fault zone, as shown in the following example: Figure 7 shown.

[0124] (5) Data analysis and statistics are performed based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone to obtain the elevation-elevation distribution strip map of the target layer. Specifically:

[0125] The target layer T74 structure isolinear map was converted into a digital elevation model DEM data map, the strike-slip fault zone width (1200m) and the strip baseline were analyzed and statistically analyzed using ArcGis software. The elevation point step size was 1.5 times the DEM unit size, and the number of parallel sections per strip was set to 50, resulting in the following: Figure 8 The elevation-elevation difference distribution strip diagram shown in the figure can be used to clearly define the maximum elevation value, minimum elevation value, average elevation value, upper and lower quartile elevation values and elevation difference value of the strip profile area, and then determine the elevation-elevation difference distribution of the target layer T74 strike-slip fault zone.

[0126] (6) Based on the target layer's elevation-elevation distribution strip map and the quantitative difference segmentation formula, a segmentation difference segmentation baseline is drawn in the target layer's elevation-elevation distribution strip map to obtain a new elevation-elevation distribution strip map. Specifically:

[0127] Obtain the elevation values of the intersection points of the average elevation strip line and the elevation difference strip line in the elevation-elevation distribution strip map of the target layer, and determine the minimum elevation value and the second smallest elevation value from the elevation values of each intersection point; based on the minimum elevation value and the second smallest elevation value, use the quantitative difference segmentation formula to determine the elevation value of the segmentation difference baseline of the strike-slip fault zone of the target layer; based on the elevation value of the segmentation difference baseline of the strike-slip fault zone of the target layer, draw the segmentation difference baseline of the strike-slip fault zone of the target layer in the elevation-elevation distribution strip map of the target layer to obtain a new elevation-elevation distribution strip map; the quantitative difference segmentation formula is:

[0128]

[0129] Among them, φ H Elevation values of the baseline divided for segment differences; is the minimum elevation value, The second smallest elevation value.

[0130] According to the elevation values of the intersection points of the average elevation strip line and the elevation difference strip line of the elevation-elevation distribution strip diagram of the target layer T74 strike-slip fault zone, the minimum elevation value and the second minimum elevation value are determined from the elevation values of each intersection point, that is, -7523m, It is also -7523m, so the segmented difference is divided into baseline elevation values δ H The new elevation-elevation distribution strip map obtained after segmenting the baseline of the strike-slip fault zone of the target layer is -7523m. Figure 9 shown.

[0131] (7) Based on the average elevation strip line and segmented difference in the new elevation-elevation distribution strip map, the baseline is divided, the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer are determined, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation distribution strip map to obtain the marked elevation-elevation distribution strip map. Specifically:

[0132] Draw a tangent to each point of the average elevation strip line in the new elevation-elevation distribution strip map, and calculate the angle between the tangent of each point and the segment difference division baseline, determine the point with an angle less than 20°, traverse the points with an angle less than 20° from left to right along the length direction of the strike-slip fault zone of the target layer, and determine whether the current point is the first point with an angle less than 20°. If it is the first point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as the continuous segment X. iThe starting point of the segment X is traversed, and the next point with an angle less than 20° is taken as the current point, where i∈1, 2, ..., I; if it is not the first point with an angle less than 20°, determine whether the current point is continuous with the previous point with an angle less than 20°. If it is not continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point of the segment X is traversed to the next point with an angle less than 20° as the current point. If it is continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The end point of the segment is traversed, and the next point with an angle less than 20° is traversed as the current point. If it is continuous with the next point with an angle less than 20°, the next point with an angle less than 20° is traversed as the current point until all points with an angle less than 20° are traversed. The continuous segment set X is obtained, X = {X1, X2, ..., X I}, I is the number of continuous segments in the continuous segment set; the continuous segments whose length is less than the preset threshold are removed from the continuous segment set to obtain the removed continuous segment set; the starting point and end point of each continuous segment in the removed continuous segment set are used as the starting point and end point of each translation segment of the strike-slip fault zone in the target layer, and the starting point and end point of each translation segment are marked in the new elevation-elevation distribution strip diagram to obtain the marked elevation-elevation distribution strip diagram. Figure 9 As shown in the figure, it is determined that the target layer T74 strike-slip fault zone has three translation segments, and the elevation values of the starting point and the ending point of each translation segment are

[0133] (8) Overlap the DEM data map of the target layer with the fault probability volume slice to form an elevation fault superposition map of the target layer. Based on the marked elevation-elevation difference distribution strip map, a segmented difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer. Specifically:

[0134] The digital elevation model DEM data map and the fault probability volume slice are superimposed and displayed to form a Figure 10 The target layer T74 elevation fault superposition map is shown in the figure. Based on the marked elevation-elevation distribution strip map of the target layer T74, the segmented difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer T74. The final division result is shown in the figure below. Figure 11As shown in Figure 3, the target layer T74 strike-slip fault zone is divided into three level I uplift segments, two level I tension depression segments, three translation segments and two level II uplift segments.

[0135] The above-mentioned quantitative segmentation method for strike-slip fault zones uses three-dimensional seismic data and actual drilling information to determine the digital elevation model (DEM) data map of the target layer, and determines the intersection relationship of primary and secondary faults and the distribution of strike-slip fault zones through the fault probability body calculated by the hybrid neural network model. The segmentation baseline is determined through the elevation-elevation difference distribution of the strike-slip fault zone and the quantitative segmentation formula of the strike-slip fault zone, achieving fine segmentation of the strike-slip fault zone and improving the accuracy of the strike-slip fault zone segmentation, providing reliable technical support for guiding the exploration and development of oil and gas reservoirs in strike-slip fault-developed areas.

[0136] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. In addition, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0137] In one embodiment, a device for quantitatively dividing a strike-slip fault zone into segments is provided, comprising:

[0138] The contour map creation module is used to create a target layer structural contour map based on the seismic data and well logging data of the target area.

[0139] The DEM data conversion module is used to perform grid analysis on the structural contour map of the target layer and convert the structural contour map of the target layer into a digital elevation model DEM data map of the target layer.

[0140] The fault probability volume analysis module is used to perform fault probability volume analysis based on seismic data using a hybrid neural network model formed by Resnet and U-net to obtain fault probability volume slices of the target layer.

[0141] The strip baseline determination module is used to determine the strip width and strip baseline of the strike-slip fault zone in the target layer according to the fault probability volume slice of the target layer.

[0142] The data analysis and statistics module is used to perform data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and obtain the elevation-elevation difference distribution strip map of the target layer.

[0143] The segmented difference division baseline drawing module is used to draw the segmented difference division baseline in the elevation-elevation difference distribution strip diagram of the target layer according to the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmented division formula to obtain a new elevation-elevation difference distribution strip diagram.

[0144] The translation segment marking module is used to divide the baseline according to the average elevation strip line and segment difference in the new elevation-elevation distribution strip map, determine the starting point and end point of the translation segment of the strike-slip fault zone in the target layer, and mark the starting point and end point of each translation segment in the new elevation-elevation distribution strip map to obtain the marked elevation-elevation distribution strip map.

[0145] The segmentation module is used to overlap the digital elevation model (DEM) data map of the target layer with the fault probability volume slice to form an elevation fault superposition map of the target layer. According to the marked elevation-elevation difference distribution strip map, the segmentation difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

[0146] The specific definitions of the strike-slip fault zone segmentation quantitative division device can be found in the definitions of the strike-slip fault zone segmentation quantitative division method described above and will not be repeated here. Each module in the aforementioned strike-slip fault zone segmentation quantitative division device can be implemented in whole or in part via software, hardware, or a combination thereof. Each of these modules can be embedded in or independent of a processor within a computer device in hardware form, or stored in a computer device memory in software form, allowing the processor to call and execute the corresponding operations of each module.

[0147] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the above-mentioned method for quantitative segmentation of a strike-slip fault zone are implemented.

[0148] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps of the above-mentioned method for quantitative segmentation of a strike-slip fault zone are implemented.

[0149] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0150] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0151] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art could make various modifications and improvements without departing from the spirit of the present application, all of which fall within the scope of protection of the present application. Therefore, the scope of protection of the present patent application shall be determined by the appended claims.

Claims

1. A method for quantitative segmentation of strike-slip fault zones, characterized in that: The strike-slip fault zone segmentation quantitative division method includes: Create structural contour maps of target layers based on seismic data and well logging data in the target area; Performing grid analysis on the structural contour map of the target horizon, and converting the structural contour map of the target horizon into a digital elevation model (DEM) data map of the target horizon; Based on the seismic data, a hybrid neural network model formed by Resnet and U-net is used to perform fault probability volume analysis to obtain a fault probability volume slice of the target horizon; Determining the strip width and strip baseline of the strike-slip fault zone of the target horizon according to the fault probability volume slice of the target horizon; Performing data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone to obtain an elevation-elevation distribution strip map of the target layer; According to the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmentation formula, a segmentation difference division baseline is drawn in the elevation-elevation difference distribution strip diagram of the target layer to obtain a new elevation-elevation difference distribution strip diagram; Determine the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer according to the average elevation strip line in the new elevation-elevation difference distribution strip diagram and the segmented difference division baseline, and mark the starting point and the ending point of each translation segment in the new elevation-elevation difference distribution strip diagram to obtain a marked elevation-elevation difference distribution strip diagram; The digital elevation model (DEM) data map of the target layer is overlapped with the fault probability volume slice to form an elevation fault superposition map of the target layer. According to the marked elevation-elevation difference distribution strip map, a segmented difference division scheme is adopted to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

2. The method for quantitative segmentation of strike-slip fault zones according to claim 1, characterized in that: The method of performing fault probability volume analysis based on the seismic data using a hybrid neural network model formed by ResNet and U-Net to obtain a fault probability volume slice of the target horizon includes: performing normalization processing on the seismic data to remove random noise and obtain denoised seismic data; Performing fault probability volume analysis on the denoised seismic data using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume; performing refinement and enhancement processing on the fault probability volume to obtain an enhanced fault probability volume; Extracting a fault probability volume slice of the target horizon along the target horizon from the enhanced fault probability volume.

3. The method for quantitative segmentation of strike-slip fault zones according to claim 1, characterized in that: The determining of the strip width and strip baseline of the strike-slip fault zone of the target horizon based on the fault probability volume slice of the target horizon includes: Determining the strike direction of the strike-slip fault zone of the target horizon according to the fault probability volume slice of the target horizon; The range covering the width of the strike-slip fault zone is taken as the strip width of the strike-slip fault zone of the target horizon; A vector profile line parallel to the strike-slip fault zone is drawn through the symmetry center of the strip width of the strike-slip fault zone of the target horizon as the strip baseline of the strike-slip fault zone of the target horizon.

4. The method for quantitative segmentation of strike-slip fault zones according to claim 1, characterized in that: The step of drawing a segmented difference division baseline in the elevation-elevation difference distribution strip diagram of the target layer according to the elevation-elevation difference distribution strip diagram of the target layer and the quantitative difference segmented division formula to obtain a new elevation-elevation difference distribution strip diagram includes: Obtaining the elevation values of the intersection points of the average elevation strip line and the elevation difference strip line in the elevation-elevation distribution strip diagram of the target layer, and determining the minimum elevation value and the second minimum elevation value from the elevation values of each intersection point; According to the minimum elevation value and the second minimum elevation value, a quantitative difference segmentation formula is used to determine the elevation value of the strike-slip fault zone segmentation difference baseline of the target layer; According to the elevation value of the segmented difference division baseline of the strike-slip fault zone of the target layer, the segmented difference division baseline of the strike-slip fault zone of the target layer is drawn in the elevation-elevation difference distribution strip diagram of the target layer to obtain a new elevation-elevation difference distribution strip diagram.

5. The method for quantitative segmentation of strike-slip fault zones according to claim 4, characterized in that: The quantitative difference segmentation formula is: Among them, φ H Elevation values of the baseline for segment differences; is the minimum elevation value, The second smallest elevation value.

6. The method for quantitative segmentation of strike-slip fault zones according to claim 1, characterized in that: The method comprises dividing the baseline of the strike-slip fault zone of the target layer according to the average elevation strip line in the new elevation-elevation difference distribution strip map and the segmented difference, determining the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer, marking the starting point and the ending point of each translation segment in the new elevation-elevation difference distribution strip map, and obtaining the marked elevation-elevation difference distribution strip map, including: Draw a tangent line for each point of the average elevation strip line in the new elevation-elevation distribution strip diagram, and calculate the angle between the tangent line of each point and the segment difference division baseline, and determine the point where the angle is less than 20°; Traverse the points with an angle less than 20° from left to right along the length direction of the strike-slip fault zone of the target horizon, and determine whether the current point is the first point with an angle less than 20°. If it is the first point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point is traversed, and the next point with an angle less than 20° is taken as the current point, where i∈1, 2, ..., I; If it is not the first point with an angle less than 20°, determine whether the current point is continuous with the previous point with an angle less than 20°. If it is not continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, continue to traverse the next point with an angle less than 20° as the current point. If it is continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The starting point is traversed, and the next point with an angle less than 20° is taken as the current point; If it is continuous with the previous point with an angle less than 20°, determine whether the current point is continuous with the next point with an angle less than 20°. If it is not continuous with the next point with an angle less than 20°, determine the current point as a continuous segment X. i The end point of the segment is traversed, and the next point with an angle less than 20° is traversed as the current point. If it is continuous with the next point with an angle less than 20°, the next point with an angle less than 20° is traversed as the current point until all points with an angle less than 20° are traversed. The continuous segment set X is obtained, X = {X1, X2, ..., X I }, I is the number of continuous segments in the continuous segment set; Eliminate the continuous segments whose length is less than a preset threshold value from the continuous segment set to obtain a continuous segment set after elimination; The starting point and the ending point of each continuous segment in the eliminated continuous segment set are used as the starting point and the ending point of each translation segment of the strike-slip fault zone of the target layer, and the starting point and the ending point of each translation segment are marked in the new elevation-elevation difference distribution strip diagram to obtain the marked elevation-elevation difference distribution strip diagram.

7. The method for quantitative segmentation of strike-slip fault zones according to claim 1, characterized in that: The segment difference division scheme is: The area below the segmentation baseline in the average elevation strip line in the marked elevation-elevation distribution strip diagram is divided into a strike-slip fault zone level I depression segment; The area above the segmentation baseline in the average elevation strip line in the marked elevation-elevation distribution strip diagram is divided into a strike-slip fault zone level I uplift segment; According to the starting point and the ending point of the translation segment in the marked elevation-elevation distribution strip diagram, the area between the starting point and the ending point of each translation segment is divided into a translation segment, and the area where each translation segment is located is further subdivided into a translation segment developed within a level I tension depression segment or a translation segment developed within a level I pressure uplift segment; In the case that there is a translation segment in the area where the level I pull-down segment is located, the segments adjacent to each translation segment in the area where the level I pull-down segment is located and lower than the translation segment are finely divided into level II pull-down segments; In the case that there is a translation segment in the area where the level I pressure uplift segment is located, the level II pressure uplift segment is finely divided from the segments higher than the translation segment in the adjacent segments of each translation segment in the area where the level I pressure uplift segment is located.

8. A device for quantitatively dividing strike-slip fault zones into segments, characterized in that: include: The contour map creation module is used to create a structural contour map of the target layer based on the seismic data and well logging data of the target area; A DEM data map conversion module is used to perform grid analysis on the structural contour map of the target layer and convert the structural contour map of the target layer into a digital elevation model DEM data map of the target layer; A fault probability volume analysis module is used to perform fault probability volume analysis based on the seismic data using a hybrid neural network model formed by Resnet and U-net to obtain a fault probability volume slice of the target horizon; a strip baseline determination module, configured to determine a strip width and a strip baseline of a strike-slip fault zone of the target horizon according to a fault probability volume slice of the target horizon; A data analysis and statistics module is used to perform data analysis and statistics based on the digital elevation model (DEM) data map of the target layer, the strip width and strip baseline of the strike-slip fault zone, and obtain an elevation-elevation distribution strip map of the target layer; A segmented difference division baseline drawing module is used to draw a segmented difference division baseline in the elevation-elevation difference distribution strip map of the target layer according to the elevation-elevation difference distribution strip map of the target layer and the quantitative difference segmentation formula to obtain a new elevation-elevation difference distribution strip map; a translation segment marking module, configured to divide the baseline according to the average elevation strip line and the segmented difference in the new elevation-elevation distribution strip diagram, determine the starting point and the ending point of the translation segment of the strike-slip fault zone of the target layer, and mark the starting point and the ending point of each translation segment in the new elevation-elevation distribution strip diagram to obtain a marked elevation-elevation distribution strip diagram; A segmentation module is used to overlap the digital elevation model (DEM) data map of the target layer with the fault probability volume slice to form an elevation fault superposition map of the target layer. According to the marked elevation-elevation difference distribution strip map, a segmentation difference division scheme is used to quantitatively and finely divide and mark the segments of the strike-slip fault zone in the elevation fault superposition map of the target layer.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for quantitative segmentation of a strike-slip fault zone according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for quantitative segmentation of a strike-slip fault zone according to any one of claims 1 to 7 are implemented.

Images