A method and device for determining the distribution of a fracture system

Abstract

The present invention provides a method and device for determining the distribution of a fault system, the method comprising: obtaining input data; the input data comprising: seismic coherent attribute slices; determining the fracture enhancement result according to the input data; determining the uniform medium background enhancement result according to the input data; determining the fracture system distribution according to the fracture enhancement result and the uniform medium background enhancement result. The present invention uses the seismic coherent attribute slices as basic data, performs enhancement processing on effective fracture features, determines the fracture enhancement result and the uniform medium background enhancement result, and obtains the fracture system distribution; the embodiments of the present invention realize the precise positioning of the main fracture development position, the effective identification of secondary fractures, and the reasonable prediction of the fracture development zone caused by the fracture.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas exploration, and in particular to a method and device for determining the distribution of a fracture system. Background Art

[0002] This section is intended to provide a background or context to embodiments of the invention that are recited in the claims. No statement herein is admitted to be prior art by inclusion in this section.

[0003] Faults, cracks, and caves are important components of fractured oil and gas reservoirs, such as heterogeneous carbonate reservoirs, and they play a crucial role in oil and gas accumulation. Fault development is a key factor controlling the generation, storage, capping, enclosure, migration, and preservation of oil and gas in a reservoir. Consequently, many of the world's largest oil and gas reservoirs are closely related to fault development. Caves, geological bodies formed by tectonic movements, karstification, and other geological processes, serve as oil and gas reservoirs within carbonate reservoirs and are a key exploration target in my country's Tarim Basin. Faults are crucial pathways for oil and gas migration, connecting dispersed reservoir spaces, such as cave systems, to form large-scale reservoirs and providing the necessary permeability within each reservoir space. Therefore, high-precision fracture detection is crucial to support subsequent exploration and development research, including modeling and reservoir simulation.

[0004] A representative example of fracture detection technology is the coherence family. Coherence technology originated in the 1990s and has evolved into three generations of algorithms: the first-generation algorithm based on cross-correlation (C1), the second-generation algorithm utilizing multi-channel similarity (C2), and the third-generation seismic coherence attribute algorithm based on characteristic structure (C3). The C1 algorithm has relatively stringent application requirements. The C2 algorithm suffers from waveform sensitivity and insensitivity to changes in lateral amplitude. The C3 algorithm, by contrast, addresses these deficiencies. Within the migrated 3D data volume, it calculates the coherence of each sample point with the surrounding data, forming a 3D data volume representing the coherence, i.e., the data coherence within the time window. This approach suppresses continuity and highlights discontinuities, while also quantitatively reflecting the lateral variations of seismic characteristics. The resulting results are more intuitive than geological interpretation based on horizontal seismic slices. This method is primarily used for more objective and detailed fracture interpretation, river channel, sand body, and fracture prediction. However, overall, the fracture detection accuracy of the C3 algorithm's coherence attribute technology is significantly affected by factors such as the sliding window size and the signal-to-noise ratio of the seismic signal. Experts, scholars, and engineers have carried out many improvements to this end. For example, they have used methods such as formation dip constraints, structural orientation constraints, frequency division, multi-directional decomposition, and ant tracking to improve the recognition of fault features. However, there are still problems such as inaccurate fault location positioning and low accuracy in identifying secondary faults. In addition, the prediction of fracture development zones caused by faults has always remained at the stage of subjective perception.

[0005] Therefore, how to provide a new solution that can solve the above technical problems is a technical problem that needs to be solved urgently in this field. Summary of the Invention

[0006] An embodiment of the present invention provides a method for determining the distribution of a fracture system, which can achieve accurate characterization of the fracture system distribution. The method includes:

[0007] Acquiring input data; the input data includes: seismic coherent attribute slices;

[0008] Determine the fracture enhancement results based on the input data;

[0009] Determine the uniform medium background enhancement result based on the input data;

[0010] The fracture system distribution is determined based on the fracture enhancement results and the uniform medium background enhancement results.

[0011] An embodiment of the present invention further provides a device for determining a fracture system distribution, comprising:

[0012] An input data acquisition module is used to acquire input data; the input data includes: seismic coherent attribute slices;

[0013] A fracture enhancement result determination module is used to determine the fracture enhancement result based on input data;

[0014] A uniform medium background enhancement result determination module is used to determine the uniform medium background enhancement result based on input data;

[0015] The fracture system distribution determination module is used to determine the fracture system distribution according to the fracture enhancement results and the uniform medium background enhancement results.

[0016] An embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for determining the fracture system distribution when executing the computer program.

[0017] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing the above-mentioned method for determining the distribution of a fracture system.

[0018] An embodiment of the present invention provides a method and device for determining the distribution of a fault system. The method uses seismic coherent attribute slices as basic data, enhances effective fault features, and determines fault enhancement results and uniform medium background enhancement results. The fracture enhancement result enhances the discontinuity characteristics of the fault and obtains as many discontinuity characteristics of suspected faults as possible. The uniform medium background enhancement result obtains regional features where faults are unlikely to occur. By combining the fracture enhancement result and the uniform medium background enhancement result, the uniform medium background enhancement result can be eliminated from the fracture enhancement result to obtain the fault system distribution. The embodiment of the present invention achieves precise positioning of the main fault development position, effective identification of secondary faults, and reasonable prediction of the crack development zone caused by the fault. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0020] Figure 1 Schematic diagram of a method for determining fracture system distribution according to an embodiment of the present invention.

[0021] Figure 2 A schematic diagram of seismic coherence attribute slices for a method for determining fault system distribution according to an embodiment of the present invention.

[0022] Figure 3 A schematic diagram of fracture enhancement results of a method for determining fracture system distribution according to an embodiment of the present invention.

[0023] Figure 4 This is a local comparison diagram of seismic coherence attribute slices and fault enhancement results in a method for determining fault system distribution according to an embodiment of the present invention.

[0024] Figure 5 This is a schematic diagram of the uniform medium background enhancement results of a method for determining the distribution of a fracture system according to an embodiment of the present invention.

[0025] Figure 6 A schematic diagram of fracture system distribution of a method for determining fracture system distribution according to an embodiment of the present invention.

[0026] Figure 7 Schematic diagram of a computer device for executing a method for determining fracture system distribution according to the present invention.

[0027] Figure 8 Schematic diagram of a device for determining fracture system distribution according to an embodiment of the present invention. DETAILED DESCRIPTION

[0028] To make the purpose, technical solutions and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0029] Figure 1 FIG. 1 is a schematic diagram of a method for determining the distribution of a fracture system according to an embodiment of the present invention. Figure 1 As shown, an embodiment of the present invention provides a method for determining the distribution of a fracture system, which can achieve accurate characterization of the distribution of the fracture system. The method includes:

[0030] Step 101: obtaining input data; the input data includes: seismic coherent attribute slices;

[0031] Step 102: Determine the fracture enhancement result based on the input data;

[0032] Step 103: Determine the uniform medium background enhancement result according to the input data;

[0033] Step 104: Determine the fracture system distribution based on the fracture enhancement results and the uniform medium background enhancement results.

[0034] An embodiment of the present invention provides a method for determining the distribution of a fault system. The method uses seismic coherent attribute slices as basic data, enhances effective fault features, and determines fault enhancement results and uniform medium background enhancement results. The fracture enhancement result enhances the discontinuity characteristics of the fault and obtains as many discontinuity characteristics of suspected faults as possible. The uniform medium background enhancement result obtains regional features where faults are unlikely to occur. By combining the fracture enhancement result and the uniform medium background enhancement result, the uniform medium background enhancement result can be eliminated from the fracture enhancement result to obtain the fault system distribution. The embodiment of the present invention achieves accurate positioning of the main fault development position, effective identification of secondary faults, and reasonable prediction of the crack development zone caused by the fault.

[0035] The present invention uses a method for characterizing the distribution of a fault system based on seismic coherent attribute slices (two-dimensional time slices or slices along the horizon) obtained using a third-generation coherence algorithm, combined with adaptive histogram equalization, restricted contrast histogram equalization, and three-level threshold control. The method may include:

[0036] Obtain input data; the input data includes: seismic coherent attribute slices; determine the fracture enhancement result based on the input data; determine the uniform medium background enhancement result based on the input data; determine the fracture system distribution based on the fracture enhancement result and the uniform medium background enhancement result.

[0037] Figure 2 This is a schematic diagram of seismic coherent attribute slices for a method of determining fault system distribution according to an embodiment of the present invention, and relates to seismic coherent attribute slices of the X work area in an example of an embodiment of the present invention; it is known that the distribution trend of the main faults in the X work area is relatively clear, but the boundaries of secondary faults, caves, etc. are unclear, and the ability to characterize the fault system is limited. After the enhanced processing of the embodiment of the present invention, the fault information at all levels is improved to varying degrees, and the final result is superior to other research methods in characterizing the fault system. Figure 2 As shown, when implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, obtaining input data includes:

[0038] Seismic coherent attribute slices are extracted from a three-dimensional seismic data volume; the seismic coherent attribute slices include: two-dimensional time slices or slices along the horizon of the third-generation seismic coherent attributes.

[0039] In the embodiment, two-dimensional time slices of seismic coherent attribute slices or slices along the horizon of third-generation seismic coherent attributes are extracted from the three-dimensional seismic data volume. The seismic coherent attribute slices can be represented as U, and the value range of U is [0,1].

[0040] Figure 3 FIG1 is a schematic diagram of fracture enhancement results of a method for determining fracture system distribution according to an embodiment of the present invention, and relates to the fracture enhancement results of the X work area in an example of an embodiment of the present invention. Figure 4 This is a local comparison diagram of seismic coherence attribute slices and fault enhancement results in a method for determining fault system distribution according to an embodiment of the present invention, involving a local comparison of seismic coherence attribute slices and fault enhancement results in areas A, B, and C of the X work area in an example of an embodiment of the present invention. Figure 3 and Figure 4 As shown, in order to obtain high-precision fracture discontinuity linear features, when implementing a method for determining fracture system distribution according to an embodiment of the present invention, in one embodiment, determining a fracture enhancement result according to input data includes:

[0041] The seismic coherent attribute slices are processed using the adaptive histogram equalization method, the fracture enhancement result and the first threshold are set, and the fracture enhancement result is determined through re-assignment processing.

[0042] To obtain high-precision discontinuous linear features, it is necessary to ensure that the relative differences between the distribution points in the target area remain unchanged while strengthening the numerical differences between the points. Specifically, the seismic coherent attribute slice U is processed using the adaptive histogram equalization method, and the fracture enhancement result U is set. att_en and the first threshold T1, through re-assignment processing, determine the fracture enhancement result U att_en .

[0043] When implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the fracture enhancement result is determined as follows:

[0044] U att_en =AHE(U)

[0045] U att_en =T1(U att_en ≥T1)

[0046] Among them, U att_en is the fracture enhancement result; T1 is the first threshold, and AHE is the adaptive histogram equalization operator.

[0047] The above-mentioned expressions for determining the fracture enhancement results are for illustration only. Those skilled in the art will appreciate that, when implemented, the above formulas may be modified in a certain form and other parameters or data may be added, or other specific formulas may be provided as needed. These variations shall all fall within the scope of protection of the present invention.

[0048] Figure 5 This is a schematic diagram of a uniform medium background enhancement result of a method for determining the distribution of a fracture system according to an embodiment of the present invention, and relates to a uniform medium background enhancement result of an X work area in an example of an embodiment of the present invention. Figure 5 As shown, to obtain a reliable homogeneous medium background, it is necessary to combine the researchers' empirical knowledge of geological characteristics and corresponding discontinuity attributes, and on the other hand, to identify the background with the help of numerical comparative analysis methods. When implementing a method for determining the distribution of a fracture system according to an embodiment of the present invention, in one embodiment, determining the homogeneous medium background enhancement result based on input data includes:

[0049] According to the seismic coherent attribute slices, the seismic coherent attribute slice background area is determined;

[0050] Slice the background area according to the seismic coherence attributes and determine the coherence attribute background removal and equalization processing results;

[0051] The background enhancement result of the uniform medium is determined based on the background equalization processing result of the coherent attribute.

[0052] When implementing a method for determining fault system distribution according to an embodiment of the present invention, in one embodiment, determining a background area of ​​a seismic coherent attribute slice based on seismic coherent attribute slices includes:

[0053] A second threshold is set based on the numerical histogram of the seismic coherence attribute slice and geological knowledge;

[0054] The seismic coherent attribute slice background area is determined according to the second threshold and the seismic coherent attribute slice.

[0055] In the embodiment, the second threshold is set based on the numerical histogram of the seismic coherent attribute slice and the geological knowledge. This embodiment mainly combines the geological background, layer information, and attribute value distribution to determine a rough background area. Next, based on the second threshold T2 and the seismic coherent attribute slice U, the seismic coherent attribute slice background area U is determined. err In the embodiment, the seismic coherence attribute slice background area U is determined as follows: err :

[0056] U err =U-T2

[0057] Among them, U err is the background area of ​​the seismic coherence attribute slice, U is the seismic coherence attribute slice; T2 is the second threshold.

[0058] The above-mentioned expressions of the background areas of seismic coherent attribute slices are for illustration only. Those skilled in the art will appreciate that, when implemented, the above formulas may be deformed in a certain form and other parameters or data may be added as needed, or other specific formulas may be provided. These variations shall all fall within the scope of protection of the present invention.

[0059] When implementing a method for determining fault system distribution according to an embodiment of the present invention, in one embodiment, slicing a background area according to seismic coherent attributes and determining a coherent attribute background removal and equalization processing result include:

[0060] The background area of ​​the seismic coherent attribute slice is processed by restricted contrast histogram equalization to determine.

[0061] When implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the result of the coherent attribute background removal and equalization processing is determined as follows:

[0062] U err_CLAHE =CLAHE(U err )

[0063] Among them, U err_CLAHE is the result of coherent attribute background removal equalization; CLAHE is the limited contrast histogram equalization operator; U err Slice the background area for seismic coherence attributes.

[0064] The above-mentioned expressions for determining the results of the coherent attribute background removal equalization processing are for illustration only. Those skilled in the art will appreciate that, during implementation, the above formulas may be modified in a certain form and other parameters or data may be added as needed, or other specific formulas may be provided. These variations shall all fall within the scope of protection of the present invention.

[0065] In the embodiment, CLAHE is a limited contrast histogram equalization operator, and the relevant parameter settings are determined according to the actual data. err The differences in the attribute interpolation between adjacent points in the image are enhanced, and the scattered points are connected into "slices". The results of this embodiment further polarize the differences between the fracture development zone and the uniform medium background.

[0066] When implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, determining a uniform medium background enhancement result based on a coherent attribute background removal and equalization processing result includes:

[0067] According to the histogram distribution and geological knowledge, the third threshold is set, and the coherent attribute background equalization processing results are threshold processed to determine the uniform medium background enhancement result.

[0068] When implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, a uniform medium background enhancement result is determined as follows:

[0069]

[0070] Among them, U back_en is the background enhancement result of uniform medium; U err_CLAHE is the result of background equalization processing of coherent attributes; T3 is the third threshold.

[0071] The above-mentioned expression for determining the background enhancement result of a uniform medium is for illustration only. Those skilled in the art will appreciate that, when implementing the expression, the above-mentioned formula may be modified in a certain form and other parameters or data may be added as needed, or other specific formulas may be provided. All of these variations shall fall within the scope of protection of the present invention.

[0072] Figure 6 This is a schematic diagram of a fracture system distribution method for determining fracture system distribution according to an embodiment of the present invention, and relates to the fracture system distribution of the X work area in an example of an embodiment of the present invention. Figure 6 As shown, when implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, determining the fracture system distribution according to the fracture enhancement result and the uniform medium background enhancement result includes:

[0073] Combining the fracture enhancement results with the uniform medium background enhancement results, the fracture system distribution is determined by multiplying the corresponding point values.

[0074] When implementing a method for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the fracture system distribution is determined as follows:

[0075] U fault =U att_en ·U back_en

[0076] Among them, U fault is the distribution of the fault system; U att_en is the fracture enhancement result; U back_en is the uniform medium background enhancement result, and is the dot product.

[0077] The above-mentioned expressions for determining the distribution of the fracture system are for illustration only. Those skilled in the art will appreciate that, when implementing the expressions, the above-mentioned formulas may be modified in a certain form and other parameters or data may be added, or other specific formulas may be provided as needed. These variations shall fall within the scope of protection of the present invention.

[0078] Fault identification is an important part of seismic interpretation in oil and gas exploration. Conventional post-stack seismic fault detection methods are subject to the influence of seismic data quality and the smoothing effect of sliding window series algorithms. They can reflect the general trend of the main fault, but the location of the fault is not accurate enough, the accuracy of identifying secondary faults is low, and the prediction of the fracture development zone caused by the fault has always remained at the stage of subjective perception. The present invention uses the time slice (or slice along the layer) of the third-generation coherent attribute results as the basic data, and through the innovative research ideas of combining adaptive histogram equalization, restricted contrast histogram equalization and three-level threshold control, enhances the effective fracture characteristics, strengthens the linear structure of the fracture, eliminates the homogeneous medium background, and achieves the goal of accurately locating the main fault development position, effectively identifying secondary faults and reasonably predicting the fracture development zone caused by the fracture. It is an innovative technology that can comprehensively characterize the distribution of the fracture system. A method for determining the distribution of the fracture system in the embodiment of the present invention is mainly reflected in (1) accurate fracture positioning; (2) strong ability to identify hidden faults; (3) the ability to characterize the fracture system and identify the fracture development zone.

[0079] In another embodiment of the present invention, another implementation method for determining the distribution of the fracture system based on the fracture enhancement results and the uniform medium background enhancement results is provided:

[0080] U fault =U att_en ∩U back_en

[0081] Among them, U fault is the distribution of the fault system; U att_en is the fracture enhancement result; U back_en is the background enhancement result of uniform medium, and (∩) is the intersection calculation.

[0082] In the above embodiment, Ufault is the distribution of the fracture system to be characterized, U att_en is the fracture enhancement result, which means the result of discontinuous property enhancement, U back_en This is the result of uniform medium background enhancement; the discontinuous feature U of the fracture is enhanced by fracture feature enhancement processing. back_en , we obtained as many discontinuous features of suspected fractures as possible; by enhancing the uniform medium background, we obtained the regional features U that are unlikely to produce fractures. back_en , the intersection (∩) of the two is calculated, and the result is equivalent to the fracture enhancement result U att_en The result of eliminating the background enhancement of uniform medium U back_en , the characteristics of the fault system distribution can be obtained.

[0083] Figure 7 A schematic diagram of a computer device for executing a method for determining fracture system distribution implemented by the present invention, as shown in FIG. Figure 7 As shown, an embodiment of the present invention further provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for determining the fracture system distribution when executing the computer program.

[0084] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing the above-mentioned method for determining the distribution of a fracture system.

[0085] The present invention also provides an apparatus for determining the distribution of a fracture system, as described in the following embodiments. Since the principles of the apparatus are similar to those of a method for determining the distribution of a fracture system, the implementation of the apparatus can be referenced to the implementation of a method for determining the distribution of a fracture system, and any repetitions will not be repeated.

[0086] Figure 8 FIG. 1 is a schematic diagram of a device for determining the distribution of a fracture system according to an embodiment of the present invention, as shown in FIG. Figure 8 As shown, an embodiment of the present invention further provides a device for determining the distribution of a fracture system, comprising:

[0087] The input data acquisition module 801 is used to acquire input data; the input data includes: seismic coherent attribute slices;

[0088] The fracture enhancement result determination module 802 is used to determine the fracture enhancement result according to the input data;

[0089] A uniform medium background enhancement result determination module 803 is used to determine a uniform medium background enhancement result based on input data;

[0090] The fracture system distribution determination module 804 is used to determine the fracture system distribution according to the fracture enhancement result and the uniform medium background enhancement result.

[0091] When implementing a device for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, an input data acquisition module is specifically configured to:

[0092] Seismic coherent attribute slices are extracted from a three-dimensional seismic data volume; the seismic coherent attribute slices include: two-dimensional time slices or slices along the horizon of the third-generation seismic coherent attributes.

[0093] When implementing a device for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, a fracture enhancement result determination module is specifically configured to:

[0094] The seismic coherent attribute slices are processed using the adaptive histogram equalization method, the fracture enhancement result and the first threshold are set, and the fracture enhancement result is determined through re-assignment processing.

[0095] When implementing an apparatus for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the fracture enhancement result determination module is further configured to determine the fracture enhancement result in the following manner:

[0096] U att_en =AHE(U)

[0097] U att_en =T1(U att_en ≥T1)

[0098] Among them, U att_en is the fracture enhancement result; T1 is the first threshold, and AHE is the adaptive histogram equalization operator.

[0099] When implementing an apparatus for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, a module for determining a uniform medium background enhancement result is specifically configured to:

[0100] According to the seismic coherent attribute slices, the seismic coherent attribute slice background area is determined;

[0101] Slice the background area according to the seismic coherence attributes and determine the coherence attribute background removal and equalization processing results;

[0102] The background enhancement result of the uniform medium is determined based on the background equalization processing result of the coherent attribute.

[0103] When implementing an apparatus for determining fracture system distribution according to an embodiment of the present invention, in one embodiment, the module for determining a uniform medium background enhancement result is further configured to:

[0104] A second threshold is set based on the numerical histogram of the seismic coherence attribute slice and geological knowledge;

[0105] The seismic coherent attribute slice background area is determined according to the second threshold and the seismic coherent attribute slice.

[0106] When implementing an apparatus for determining fracture system distribution according to an embodiment of the present invention, in one embodiment, the module for determining a uniform medium background enhancement result is further configured to:

[0107] The restricted contrast histogram equalization processing is performed on the background area of ​​the seismic coherent attribute slice to determine the coherent attribute background removal equalization processing result.

[0108] When implementing an apparatus for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the module for determining a uniform medium background enhancement result is further configured to determine a coherent attribute background removal and equalization processing result in the following manner:

[0109] U err_CLAHE =CLAHE(U err )

[0110] Among them, U err_CLAHE is the result of coherent attribute background removal equalization; CLAHE is the limited contrast histogram equalization operator; U err Slice the background area for seismic coherence attributes.

[0111] When implementing an apparatus for determining fracture system distribution according to an embodiment of the present invention, in one embodiment, the module for determining a uniform medium background enhancement result is further configured to:

[0112] According to the histogram distribution and geological knowledge, the third threshold is set, and the coherent attribute background equalization processing results are threshold processed to determine the uniform medium background enhancement result.

[0113] When implementing an apparatus for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the uniform medium background enhancement result determination module is further configured to determine the uniform medium background enhancement result in the following manner:

[0114]

[0115] Among them, U back_en is the background enhancement result of uniform medium; U err_CLAHE is the result of background equalization processing of coherent attributes; T3 is the third threshold.

[0116] When implementing a device for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, a fracture system distribution determination module is specifically configured to:

[0117] Combining the fracture enhancement results with the uniform medium background enhancement results, the fracture system distribution is determined by multiplying the corresponding point values.

[0118] When implementing a device for determining a fracture system distribution according to an embodiment of the present invention, in one embodiment, the fracture system distribution determining module is further configured to determine the fracture system distribution in the following manner:

[0119] U fault =U att_en ·U back_en

[0120] Among them, U fault is the distribution of the fault system; U att_en is the fracture enhancement result; U back_en is the uniform medium background enhancement result, and is the dot product.

[0121] In summary, the embodiment of the present invention provides a method and device for determining the distribution of a fault system. The method uses seismic coherent attribute slices as basic data, enhances the effective fault features, and determines the fault enhancement results and the uniform medium background enhancement results. The fracture enhancement results enhance the discontinuity characteristics of the fault and obtain as many discontinuity characteristics of suspected faults as possible. The uniform medium background enhancement results obtain regional features where faults are unlikely to occur. The fracture enhancement results and the uniform medium background enhancement results are combined to eliminate the uniform medium background enhancement results from the fracture enhancement results to obtain the fault system distribution. The embodiment of the present invention uses an innovative research approach that combines adaptive histogram equalization, restricted contrast histogram equalization, and three-level threshold control to enhance the effective fault features, strengthen the linear structure of the fault, eliminate the uniform medium background, and achieve precise positioning of the main fault development position, effective identification of secondary faults, and reasonable prediction of the crack development zone caused by the fault.

[0122] It will be understood by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0123] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0124] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0125] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0126] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for determining the distribution of a fracture system, characterized in that: include: Get input data; The input data includes: seismic coherent attribute slices; Determine the fracture enhancement results based on the input data; Determine the uniform medium background enhancement result according to the input data; Determine the distribution of the fracture system based on the fracture enhancement results and the uniform medium background enhancement results; Among them, according to the input data, the fracture enhancement results are determined, including: The seismic coherent attribute slices are processed by using an adaptive histogram equalization method, the fracture enhancement result and the first threshold are set, and the fracture enhancement result is determined by re-assignment processing; Based on the input data, determine the uniform medium background enhancement results, including: According to the seismic coherent attribute slices, the seismic coherent attribute slice background area is determined; according to the seismic coherent attribute slice background area, the coherent attribute background removal equalization processing result is determined; according to the coherent attribute background removal equalization processing result, the uniform medium background enhancement result is determined; According to the seismic coherent attribute slices, the seismic coherent attribute slice background area is determined, including: According to the numerical histogram of the seismic coherent attribute slice and the geological knowledge, the second threshold is set; according to the second threshold and the seismic coherent attribute slice, press U err = U-T2, determine the background area of ​​seismic coherent attribute slices, U err is the background area of ​​seismic coherence attribute slice, U is the seismic coherence attribute slice, and T2 is the second threshold; According to the coherent attribute background equalization processing results, the uniform medium background enhancement results are determined, including: According to the histogram distribution and geological knowledge, the third threshold is set, and the coherent attribute background removal equalization processing result is threshold processed to determine the uniform medium background enhancement result; wherein, the uniform medium background enhancement result is determined in the following manner: Among them, U back_en is the background enhancement result of uniform medium; U err_CLAHE is the result of coherent attribute background removal and equalization processing; T3 is the third threshold.

2. The method according to claim 1, characterized in that The obtaining of input data comprises: Seismic coherent attribute slices are extracted from the three-dimensional seismic data body; the seismic coherent attribute slices include: two-dimensional time slices or slices along the horizon of the third-generation seismic coherent attributes.

3. The method according to claim 1, characterized in that Determine the fracture enhancement result as follows: YOU att_en =AHE(U) IN att_en =T1(U att_en ≥T1) Among them, U att_en is the fracture enhancement result; T1 is the first threshold, and AHE is the adaptive histogram equalization operator.

4. The method according to claim 1, characterized in that Slice the background area according to the seismic coherent attributes and determine the coherent attribute background removal equalization processing results, including: The background area of ​​the seismic coherent attribute slice is subjected to restricted contrast histogram equalization processing to determine the coherent attribute background removal equalization processing result.

5. The method according to claim 4, characterized in that The coherent attribute background removal equalization processing result is determined as follows: YOU err_CLAHE =CLAHE(U err ) Among them, U err_CLAHE is the result of coherent attribute background removal equalization; CLAHE is the contrast-constrained histogram equalization operator; U err Slice the background region for seismic coherence attributes.

6. The method according to claim 1, characterized in that According to the fracture enhancement results and the uniform medium background enhancement results, the fracture system distribution is determined, including: Combining the fracture enhancement results with the uniform medium background enhancement results, the fracture system distribution is determined by multiplying the corresponding point values.

7. The method according to claim 6, characterized in that The fracture system distribution is determined as follows: IN fault =U att_en ×U back_en Among them, U fault is the fault system distribution; U att_en is the fracture enhancement result; U back_en is the result of uniform medium background enhancement, and × is the dot product.

8. A device for determining the distribution of a fracture system, characterized in that include: An input data acquisition module, used for acquiring input data; The input data includes: seismic coherent attribute slices; A fracture enhancement result determination module is used to determine the fracture enhancement result according to the input data; A uniform medium background enhancement result determination module is used to determine the uniform medium background enhancement result according to the input data; A fracture system distribution determination module is used to determine the fracture system distribution according to the fracture enhancement result and the uniform medium background enhancement result; Among them, the fracture enhancement result determination module is specifically used for: The seismic coherent attribute slices are processed by using an adaptive histogram equalization method, the fracture enhancement result and the first threshold are set, and the fracture enhancement result is determined by re-assignment processing; The homogeneous medium background enhancement result determination module is specifically used for: According to the seismic coherent attribute slices, the seismic coherent attribute slice background area is determined; according to the seismic coherent attribute slice background area, the coherent attribute background removal equalization processing result is determined; according to the coherent attribute background removal equalization processing result, the uniform medium background enhancement result is determined; The homogeneous medium background enhancement result determination module is also used for: According to the numerical histogram of the seismic coherent attribute slice and the geological knowledge, the second threshold is set; according to the second threshold and the seismic coherent attribute slice, press U err = U-T2, determine the background area of ​​seismic coherent attribute slices, U err is the background area of ​​seismic coherence attribute slice, U is the seismic coherence attribute slice; T2 is the second threshold; The homogeneous medium background enhancement result determination module is also used for: According to the histogram distribution and geological knowledge, the third threshold is set, and the coherent attribute background removal equalization processing result is threshold processed to determine the uniform medium background enhancement result; wherein, the uniform medium background enhancement result is determined in the following manner: Among them, U back_en is the background enhancement result of uniform medium; U err_CLAHE is the result of coherent attribute background removal and equalization processing; T3 is the third threshold.

9. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for determining the fracture system distribution according to any one of claims 1 to 7 when executing the computer program.

10. A computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the method for determining a fracture system distribution according to any one of claims 1 to 7.

Images