Fracture detection method and device

Abstract

The invention relates to the technical field of oil exploration and development, and particularly discloses a fracture detection method and device, and the method comprises the steps: obtaining a three-dimensional seismic coherent attribute body of a target region; extracting a plurality of two-dimensional profiles from the three-dimensional seismic coherent attribute body; calculating a seismic coherence attribute reciprocal corresponding to each two-dimensional section in the plurality of two-dimensional sections; determining transverse fracture indication factor data corresponding to each two-dimensional section according to the seismic coherence attribute reciprocal corresponding to each two-dimensional section by using a preset transverse judgment rule; constructing a longitudinal judgment rule on the basis of the transverse fracture indication factor data corresponding to each two-dimensional profile; determining longitudinal fracture indication factor data corresponding to each two-dimensional profile by utilizing the longitudinal judgment rule; and based on the longitudinal fracture indication factor data corresponding to each two-dimensional profile, generating a fracture detection result corresponding to the three-dimensional seismic coherent attribute body. According to the scheme, the fracture detection precision can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of oil exploration and development, and in particular to a fracture detection method and device. Background Art

[0002] Fractures are important parameters for describing the petrophysical properties of oil and gas reservoirs. The detailed characterization of fractures is of great significance for oil and gas geophysical exploration, especially for unconventional oil and gas reservoirs. How to effectively identify and characterize fractures is still a major problem. Fracture detection helps to understand the basic oil and gas accumulation elements such as reservoir formation, oil and gas migration and accumulation, and preservation, and is an important part of oil and gas exploration. Conventional post-stack discontinuity attribute extraction is the main technology for fracture detection, and the most representative one is the eigenvalue coherence technology. This series of methods can indicate the more obvious characteristics of the phase axis misalignment, shape change, bifurcation, merging, and distortion caused by the development of fractures. Therefore, it has been widely used, and in practice, a number of derivative technologies such as histogram enhancement, direction enhancement, and frequency enhancement have been generated. In addition, similar to ant tracking technology and maximum likelihood volume attribute technology, secondary optimization based on discontinuous attribute data can also greatly improve the accuracy of fracture detection and be widely used in production practice.

[0003] However, the existing coherent properties are not accurate enough to characterize the fracture system and cannot obtain high-resolution fracture detection results.

[0004] To address the above problems, no effective solution has been proposed yet. Summary of the invention

[0005] The embodiments of this specification provide a fracture detection method and device to solve the problem of insufficient accuracy in describing the fracture system by coherent attributes in the prior art.

[0006] The present invention provides a fracture detection method, including:

[0007] Acquire a three-dimensional seismic coherent attribute volume of the target area; extract multiple two-dimensional sections from the three-dimensional seismic coherent attribute volume; calculate the inverse of the seismic coherent attribute corresponding to each two-dimensional section in the multiple two-dimensional sections;

[0008] Determine the transverse fracture indicator factor data corresponding to each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section by using the preset transverse determination rule;

[0009] Based on the transverse fracture indicator factor data corresponding to each two-dimensional section, a longitudinal determination rule is constructed; and by using the longitudinal determination rule, the longitudinal fracture indicator factor data corresponding to each two-dimensional section is determined;

[0010] Based on the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections, a fracture detection result corresponding to the three-dimensional seismic coherent attribute volume is generated.

[0011] In one embodiment, the preset horizontal determination rule is constructed in the following manner:

[0012] Set up multiple modes according to the lateral distribution morphology of the fault on the inverse of the seismic coherence attribute;

[0013] A corresponding fracture indication factor value is set for each of the multiple modes.

[0014] In one embodiment, using a preset transverse determination rule, according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section, determining the transverse fracture indicator factor data corresponding to each two-dimensional section includes:

[0015] Determine the length of the lateral sliding window according to the horizontal fault distance of the fracture in the target area in the lateral direction;

[0016] For each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the transverse window, so as to obtain the first vector and the second vector corresponding to each point on each of the two-dimensional sections;

[0017] Determine, according to the first vector and the second vector corresponding to each point on each two-dimensional section, a value of a transverse fracture indication factor corresponding to each point on each two-dimensional section;

[0018] According to the transverse fracture indication factor values ​​corresponding to each point on each two-dimensional section, the transverse fracture indication factor data corresponding to each two-dimensional section are obtained.

[0019] In one embodiment, determining the transverse fracture indication factor value corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section includes:

[0020] Calculating the gradients of the first vector and the second vector corresponding to each point on each of the two-dimensional sections to obtain the first vector gradient and the second vector gradient;

[0021] Obtaining a first fracture measure corresponding to each point on each of the two-dimensional sections according to the absolute value of the sum of the positive values ​​and the sum of the negative values ​​included in the first vector gradient and the second vector gradient;

[0022] Obtaining a second fracture measure corresponding to each point on each of the two-dimensional sections according to whether the first vector gradient and the second vector gradient are all positive numbers or all negative numbers;

[0023] Based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a transverse fracture indication factor value corresponding to each point on each two-dimensional section is determined.

[0024] In one embodiment, for each of the multiple points on each of the two-dimensional sections, the length of the transverse sliding window is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window, including:

[0025] Using the length of the transverse sliding window, slide along each two-dimensional section from left to right and from top to bottom, respectively extending to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window;

[0026] Correspondingly, based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a first transverse fracture indication factor value corresponding to each point on each two-dimensional section is determined.

[0027] In one embodiment, the method further comprises:

[0028] Using the length of the transverse sliding window, slide along each two-dimensional section from right to left and from top to bottom, respectively extending to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window;

[0029] Accordingly, based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a second transverse fracture indicator factor value corresponding to each point on each two-dimensional section is determined;

[0030] The method further comprises:

[0031] The larger value of the first transverse fracture indication factor value corresponding to each point on each two-dimensional section and the second transverse fracture indication factor value corresponding to each point on each two-dimensional section is used as the transverse fracture indication factor value corresponding to each point on each two-dimensional section.

[0032] In one embodiment, based on the transverse fracture indicator factor data corresponding to each two-dimensional section, a longitudinal determination rule is constructed, including:

[0033] Determining the longitudinal extension depth of the fracture in the target area according to the transverse fracture indicator factor data corresponding to each two-dimensional section;

[0034] Determining the length of the longitudinal sliding window according to the longitudinal extension depth of the fracture in the target area;

[0035] Based on the length of the vertical sliding window, a vertical decision rule is constructed.

[0036] In one embodiment, using the longitudinal determination rule to determine the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections includes:

[0037] The longitudinal sliding window length is utilized to slide along each two-dimensional section from top to bottom and from left to right. Each time the window slides, the transverse fracture indicator factor data in the longitudinal window is extracted to obtain a longitudinal vector. The longitudinal fracture indicator factor data corresponding to each two-dimensional section is determined based on the number of non-zero elements in two adjacent longitudinal vectors.

[0038] The embodiment of this specification also provides a fracture detection device, including:

[0039] An acquisition module is used to acquire a three-dimensional seismic coherent attribute volume of a target area; extract a plurality of two-dimensional sections from the three-dimensional seismic coherent attribute volume; and calculate the inverse of the seismic coherent attribute corresponding to each of the plurality of two-dimensional sections;

[0040] A first determination module is used to determine the transverse fracture indicator factor data corresponding to each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section by using a preset transverse determination rule;

[0041] The second determination module is used to construct a longitudinal determination rule based on the transverse fracture indication factor data corresponding to each of the two-dimensional sections; and is also used to determine the longitudinal fracture indication factor data corresponding to each of the two-dimensional sections using the longitudinal determination rule;

[0042] A generating module is used to generate a fracture detection result corresponding to the three-dimensional seismic coherent attribute volume based on the longitudinal fracture indicator factor data corresponding to each two-dimensional profile.

[0043] An embodiment of the present specification also provides a computer device, including a processor and a memory for storing processor executable instructions, wherein the processor implements the steps of the fracture detection method described in any of the above embodiments when executing the instructions.

[0044] The embodiments of the present specification also provide a computer-readable storage medium on which computer instructions are stored. When the instructions are executed, the steps of the fracture detection method described in any of the above embodiments are implemented.

[0045] In an embodiment of the present specification, a fracture detection method is provided, which can obtain a three-dimensional seismic coherent attribute body of a target area, extract multiple two-dimensional sections from the three-dimensional seismic coherent attribute body, calculate the inverse of the seismic coherent attribute corresponding to each of the multiple two-dimensional sections, use a preset horizontal judgment rule, and determine the horizontal fracture indicator factor data corresponding to each of the two-dimensional sections according to the inverse of the seismic coherent attribute corresponding to each of the two-dimensional sections. Based on the horizontal fracture indicator factor data corresponding to each of the two-dimensional sections, a vertical judgment rule is constructed, and the vertical fracture indicator factor data corresponding to each of the two-dimensional sections is determined by using the vertical judgment rule. Based on the vertical fracture indicator factor data corresponding to each of the two-dimensional sections, a fracture detection result corresponding to the three-dimensional seismic coherent attribute body is generated. In the above scheme, the typical horizontal morphological features of the fracture in the coherent attribute are used as a reference, the pattern of the discontinuous features of the coherent attribute is enumerated, and a horizontal judgment criterion is formed, so as to match the measured coherent attribute and assign the corresponding fracture indicator factor value, and the accuracy of the fracture position is ensured by taking the best of the forward and reverse judgment results. In addition, considering the ductility of the fracture in the longitudinal direction, a longitudinal judgment criterion is set on the basis of the transverse judgment result to determine the ductility of the fracture in the longitudinal direction. By combining the transverse judgment with the longitudinal judgment, a more accurate fracture detection result of the measured data is finally achieved. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The drawings described herein are used to provide a further understanding of this specification, constitute a part of this specification, and do not constitute a limitation of this specification. In the drawings:

[0047] Figure 1 A flow chart of a fracture detection method in an embodiment of this specification is shown;

[0048] Figure 2 A schematic diagram showing 16 modes of setting according to the lateral distribution morphology of typical fractures on the inverse attribute in an embodiment of this specification is shown;

[0049] Figure 3 A schematic diagram showing a first horizontal determination and a second horizontal determination in an embodiment of the present specification is shown;

[0050] Figure 4 A schematic diagram showing a longitudinal determination in an embodiment of the present specification is shown;

[0051] Figure 5 A flow chart of a fracture detection method in an embodiment of this specification is shown;

[0052] Figure 6 shows a seismic coherent attribute profile image before the implementation of the fracture detection method in an embodiment of this specification;

[0053] Figure 7 The lateral determination result after the fracture detection method in one embodiment of this specification is implemented is shown;

[0054] Figure 8 The longitudinal determination result after the fracture detection method in one embodiment of this specification is implemented is shown;

[0055] Fig. 9 It shows the coherent attribute time slice before the implementation of the fracture detection method in an embodiment of this specification;

[0056] Fig.10 It shows a high-resolution fracture feature slice after the fracture detection method in one embodiment of this specification is implemented;

[0057] Fig.11 A schematic diagram of a fracture detection device in an embodiment of this specification is shown;

[0058] Fig.12 A schematic diagram of a computer device in an embodiment of the present specification is shown. DETAILED DESCRIPTION

[0059] The principles and spirit of this specification will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are provided only to enable those skilled in the art to better understand and implement this specification, and are not intended to limit the scope of this specification in any way. On the contrary, these embodiments are provided to make this specification more thorough and complete, and to fully convey the scope of this disclosure to those skilled in the art.

[0060] Those skilled in the art will appreciate that the embodiments of this specification may be implemented as a system, device, method, or computer program product. Therefore, this specification may be implemented in the following forms: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0061] The embodiments of this specification provide a fracture detection method. Figure 1A flow chart of a fracture detection method in an embodiment of the present specification is shown. Although the present specification provides method operation steps or device structures as shown in the following embodiments or drawings, more or fewer operation steps or module units may be included in the method or device based on routine or without creative labor. In the steps or structures where there is no necessary causal relationship logically, the execution order of these steps or the module structure of the device is not limited to the execution order or module structure described in the embodiments of the present specification and shown in the drawings. When the method or module structure is applied to an actual device or terminal product, it can be connected according to the method or module structure shown in the embodiments or drawings for sequential execution or parallel execution (for example, a parallel processor or a multi-threaded processing environment, or even a distributed processing environment).

[0062] Specifically, Figure 1 As shown, a fracture detection method provided in one embodiment of this specification may include the following steps:

[0063] Step S101, obtaining a three-dimensional seismic coherent attribute volume of a target area; extracting multiple two-dimensional sections from the three-dimensional seismic coherent attribute volume; and calculating the inverse of the seismic coherent attribute corresponding to each of the multiple two-dimensional sections.

[0064] In this embodiment, a three-dimensional seismic coherence attribute body of the target area can be obtained. The three-dimensional seismic coherence attribute body refers to the mutual correlation coefficient of adjacent seismic traces in the three-dimensional seismic data body. The three dimensions of the three-dimensional seismic coherence attribute body are xline offset, inline offset and time. Multiple two-dimensional profiles can be extracted from the three-dimensional seismic coherence attribute body. The horizontal direction of the two-dimensional profile can be the xline offset, and the vertical direction can be the time. The inverse of the seismic coherence attribute corresponding to each two-dimensional profile in the multiple two-dimensional profiles is calculated to obtain the seismic coherence inverse attribute.

[0065] Step S102, using a preset transverse determination rule, according to the inverse of the seismic coherence attribute corresponding to each of the two-dimensional sections, determines the transverse fracture indicator factor data corresponding to each of the two-dimensional sections.

[0066] The preset transverse determination rules can be used to determine the transverse fracture indicator factor data of each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section.

[0067] In some embodiments of the present specification, the preset lateral determination rule may be constructed in the following manner: setting a plurality of modes according to the lateral distribution morphology of the fracture on the inverse of the seismic coherence attribute; and setting a corresponding fracture indication factor value for each of the plurality of modes.

[0068] Specifically, 16 modes can be set according to the lateral distribution of typical fractures on the inverse attribute, and the corresponding fracture indication factor values ​​can be assigned. Figure 2 , a schematic diagram of 16 modes is shown.

[0069] Assume that the length of the one-dimensional sliding window is xstep (xstep is a positive odd number), and at any point (i, j) on the coherent section, extend (xstep-1) / 2 grid points on both sides to extract the coherent attribute value. Take the coherent value from the point (i, j) to the left to get the vector L; take the coherent value from the point (i, j) to the right to get the vector R, calculate the gradients of L and R respectively, and get ZL and ZR.

[0070] Calculate the first fracture metric characterizing the fracture: sum the positive and negative values ​​contained in ZL and take the absolute value to obtain ZL_P_sum and ZL_N_sum; sum the positive and negative values ​​contained in ZR and take the absolute value to obtain ZR_P_sum and ZR_N_sum.

[0071] Calculate the second fracture metric that characterizes the fracture: count whether ZL and ZR are all positive or negative. If they are all positive or all negative, the value is 1; otherwise, the value is 0, that is:

[0072] If ZL are all positive numbers, set the parameter ZLZ1=1;

[0073] If ZL is all negative, set parameter ZLZ2=1;

[0074] If ZR is all positive, set parameter ZRZ1=1;

[0075] If ZR is all negative, set parameter ZRZ2=1.

[0076] Based on the two measures of fracture characterization given above, a numerical value of the indicator factor for fracture is given to the point according to the following scheme:

[0077]

[0078] Among them, “==”, “>”, and “<” are relational operators, which respectively mean “equal to”, “greater than”, and “less than”; “&&” is a logical AND operation, which means that the expressions on both sides must be satisfied at the same time for them to be valid.

[0079] Regarding the value of the indicator factor, the situation in different work areas cannot usually be generalized. The above reference values ​​given in this embodiment are applicable to the situation where the general seismic data conditions are stable.

[0080] In some embodiments of the present specification, using preset horizontal judgment rules, according to the inverse of the seismic coherence attributes corresponding to each two-dimensional section, determining the horizontal fracture indicator factor data corresponding to each two-dimensional section may include: determining the length of the horizontal sliding window according to the horizontal fault distance of the fracture in the target area; for each of the multiple points on each two-dimensional section, using the horizontal sliding window length to extend to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the horizontal window, and obtaining the first vector and the second vector corresponding to each point on each two-dimensional section; determining the horizontal fracture indicator factor value corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section; obtaining the horizontal fracture indicator factor data corresponding to each two-dimensional section according to the horizontal fracture indicator factor value corresponding to each point on each two-dimensional section.

[0081] The length of the transverse sliding window can be determined according to the horizontal fault size of the fault in the target area in the transverse direction. If the horizontal fault is large, then the length of the transverse sliding window is large, otherwise it is small. Afterwards, the transverse judgment can be performed using the transverse sliding window according to the transverse judgment rule. For each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherent attribute in the transverse window, and the first vector L and the second vector R corresponding to each point on each of the two-dimensional sections are obtained. According to the first vector L and the second vector R corresponding to each point on each of the two-dimensional sections, the transverse fracture indicator factor value corresponding to each point on each of the two-dimensional sections is determined. According to the transverse fracture indicator factor value corresponding to each point on each of the two-dimensional sections, the transverse fracture indicator factor data corresponding to each of the two-dimensional sections is obtained.

[0082] In some embodiments of the present specification, determining the numerical value of the transverse fracture indication factor corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section may include: calculating the gradients of the first vector and the second vector corresponding to each point on each two-dimensional section to obtain the first vector gradient and the second vector gradient; obtaining the first fracture measure corresponding to each point on each two-dimensional section according to the absolute value of the sum of the positive values ​​and the sum of the negative values ​​contained in the first vector gradient and the second vector gradient; obtaining the second fracture measure corresponding to each point on each two-dimensional section according to whether the first vector gradient and the second vector gradient are all positive numbers or all negative numbers; determining the numerical value of the transverse fracture indication factor corresponding to each point on each two-dimensional section based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section.

[0083] Specifically, the gradients of the first vector L and the second vector R can be calculated respectively to obtain the gradients ZL and ZR of the first vector and the second vector. The positive and negative values ​​contained in ZL are summed respectively, and the absolute values ​​are taken to obtain ZL_P_sum and ZL_N_sum; the positive and negative values ​​contained in ZR are summed respectively, and the absolute values ​​are taken to obtain ZR_P_sum and ZR_N_sum. The first fracture measure includes four variables: ZL_P_sum, ZL_N_sum, ZR_P_sum and ZR_N_sum. The first fracture measure indicates the coherent numerical fluctuations contained in all points contained in the window passing through the study point, and is used for fracture modes except modes 1, 2, 11, and 12.

[0084] The second fracture measure is calculated as follows: It can be counted whether ZL and ZR are all positive or all negative. If they are all positive or all negative, the value is 1; otherwise, the value is 0, that is:

[0085] If ZL are all positive numbers, set the parameter ZLZ1=1;

[0086] If ZL is all negative, set parameter ZLZ2=1;

[0087] If ZR is all positive, set parameter ZRZ1=1;

[0088] If ZR is all negative, set parameter ZRZ2=1.

[0089] The second fracture measure may include ZLZ1, ZLZ2, ZRZ1, or ZRZ2.

[0090] After obtaining the first fracture measurement and the second fracture measurement, the value of the fracture indication factor of each point on the profile can be obtained according to the above formula (1).

[0091] In some embodiments of the present specification, for each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the transverse window, which can include: using the transverse sliding window length to slide from left to right and from top to bottom along each of the two-dimensional sections, and extending to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the transverse window; accordingly, based on the first fracture measurement and the second fracture measurement corresponding to each point on each of the two-dimensional sections, the first transverse fracture indicator factor value corresponding to each point on each of the two-dimensional sections is determined.

[0092] Specifically, any relevant attribute profile data C can be selected M×N, where M is the number of longitudinal time samples and N is the number of transverse spatial samples. Set the sliding window length to xstep (xstep is an odd number), slide along the coherent attribute profile from left to right and from top to bottom, and extract the coherent value each time you slide to form a decision unit vector c i , where i = 1, 2, ..., xstep. According to the horizontal judgment rule, c i Make a judgment and assign the corresponding fracture indication factor. Iterate the above steps until the entire C M×N The first transverse determination is completed and the result is recorded as the value of the first transverse fracture indication factor.

[0093] In some embodiments of the present specification, the method may further include: utilizing the length of the transverse sliding window to slide along each two-dimensional section from right to left and from top to bottom, and extending to the left and right sides of each point respectively to extract the inverse value of the seismic coherent attribute within the transverse window; correspondingly, based on the first fracture measurement and the second fracture measurement corresponding to each point on each two-dimensional section, determining the second transverse fracture indicator factor value corresponding to each point on each two-dimensional section; the method may further include: using the larger value of the first transverse fracture indicator factor value corresponding to each point on each two-dimensional section and the second transverse fracture indicator factor value corresponding to each point on each two-dimensional section as the transverse fracture indicator factor value corresponding to each point on each two-dimensional section.

[0094] Please refer to Figure 3 , shows a schematic diagram of the first horizontal determination and the second horizontal determination in an embodiment of this specification. Figure 3 As shown, in this embodiment, in the horizontal direction, a sliding window can be used to first make a forward judgment and then a reverse judgment, take the maximum of the two results, assign the fracture indication factor value point by point, and obtain the horizontal judgment result. After obtaining the first horizontal fracture indication factor value, the data matrix C of the profile is M×N Flip left and right, repeat the above steps until the entire section is completed, complete the second transverse determination, and record the result as the second transverse fracture indication factor value corresponding to each point on each two-dimensional section. Extract the maximum value of the first transverse fracture indication factor value and the second transverse fracture indication factor value point by point as the fracture indication factor of the point, and record it as the transverse fracture indication factor value corresponding to each point.

[0095] Step S103, constructing a longitudinal determination rule based on the transverse fracture indication factor data corresponding to each two-dimensional section; and determining the longitudinal fracture indication factor data corresponding to each two-dimensional section by using the longitudinal determination rule.

[0096] After obtaining the transverse fracture indication factor data corresponding to each two-dimensional section, a longitudinal determination rule may be constructed based on the transverse fracture indication factor data corresponding to each two-dimensional section. Then, the longitudinal determination rule and the transverse fracture indication factor data corresponding to each two-dimensional section are used to determine the longitudinal fracture indication factor data corresponding to each two-dimensional section.

[0097] In some embodiments of the present specification, based on the transverse fracture indicator factor data corresponding to each two-dimensional section, constructing a longitudinal determination rule may include: determining the longitudinal extension depth of the fracture in the target area according to the transverse fracture indicator factor data corresponding to each two-dimensional section; determining the longitudinal sliding window length according to the longitudinal extension depth of the fracture in the target area; and constructing a longitudinal determination rule based on the longitudinal sliding window length.

[0098] In this embodiment, based on the transverse fracture indicator factor data corresponding to each two-dimensional section, the sliding window length is set to tstep (tstep is an odd number), and the coherent attribute section is slid from top to bottom and from left to right. Each time it slides, the coherent value is extracted to form a judgment unit vector. i 、F i+1 Two vectors, where i = 1, 2, …, tstep. The setting of tstep needs to be based on the vertical extension depth of the fault in the study area. The longer a fault extends vertically, the farther its impact is, and the larger tstep is. On the contrary, if the faults in a region are small and the impact depth is small, then tstep will be smaller.

[0099] Construct vertical judgment criteria for F i 、F i+1 Make a judgment and assign the corresponding fracture indication factor. Assume that F i 、F i+1 The number of non-zero elements in is N i 、Ni +1 , for any point i, its fracture indicator factor FAULT i It can be expressed as follows:

[0100]

[0101] Among them, max is the maximum value operation of vector elements, and mean is the mean operation of vector orientation.

[0102] In some embodiments of the present specification, using the longitudinal judgment rule to determine the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections may include: using the longitudinal sliding window length to slide along each of the two-dimensional sections from top to bottom and from left to right, each time sliding, extracting the transverse fracture indicator factor data in the longitudinal window to obtain a longitudinal vector, and determining the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections based on the number of non-zero elements in two adjacent longitudinal vectors.

[0103] Please refer to Figure 4 , shows a schematic diagram of longitudinal determination in an embodiment of the present specification. After the longitudinal determination rule is constructed, the longitudinal sliding window length can be used to slide along each two-dimensional section from top to bottom and from left to right. Each time the sliding is performed, the transverse fracture indicator factor data in the longitudinal window is extracted to obtain the longitudinal vector. According to the number of non-zero elements in two adjacent longitudinal vectors, the longitudinal fracture indicator factor data corresponding to each two-dimensional section is determined according to the above formula (2).

[0104] Step S104: generating a fracture detection result corresponding to the three-dimensional seismic coherent attribute volume based on the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections.

[0105] After obtaining the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections, the fracture detection result corresponding to the three-dimensional seismic coherent attribute volume can be generated based on the longitudinal fracture indicator factor.

[0106] In the above embodiment, the typical transverse morphological features of the fracture in the coherent attributes are used as references, the patterns of discontinuous features of the coherent attributes are enumerated, and the transverse determination criteria are formed, so as to match the measured coherent attributes and assign corresponding fracture indication factor values, and the accuracy of the fracture position is ensured by taking the best of the forward and reverse determination results. In addition, considering the longitudinal ductility of the fracture, the longitudinal determination criteria are set on the basis of the transverse determination results to determine the longitudinal ductility of the fracture, and by combining the transverse determination with the longitudinal determination, a more accurate fracture detection result of the measured data is finally achieved.

[0107] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. For details, please refer to the description of the above-mentioned related processing related embodiments, and no further description is given here.

[0108] The above is a description of a specific embodiment of the specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0109] The above method is described below in conjunction with a specific embodiment. However, it should be noted that the specific embodiment is only for better illustrating the present specification and does not constitute an improper limitation on the present specification.

[0110] This specific embodiment discloses a fracture detection method. Figure 5 , shows a flow chart of a fracture detection method in one embodiment of this specification. Figure 5 As shown, the method may include the steps of: obtaining a three-dimensional coherent attribute body, arbitrarily extracting its two-dimensional section, and calculating the inverse attribute; constructing a judgment criterion for fracture in the transverse direction. According to the transverse distribution morphology of typical fractures on the inverse attribute, 16 modes are set and the corresponding fracture indication factor values ​​are assigned; in the transverse direction, a forward judgment is first made in the form of a sliding window, and then a reverse judgment is made, and the maximum value of the two results is taken, and the fracture indication factor value is assigned point by point to obtain the transverse judgment result; based on the transverse judgment result, a judgment criterion for fracture in the longitudinal direction is constructed. In the form of a sliding window, the fracture indication factor value is assigned point by point, thereby strengthening the longitudinal ductility of the fracture; it is iterated cyclically by section, and the above process is repeated to obtain a three-dimensional fracture detection result.

[0111] Construct the judgment criteria of fracture mode. According to the lateral distribution morphology of typical fractures on the inverse attribute, 16 modes are set and the corresponding fracture indication factor values ​​are assigned. Figure 2 shown.

[0112] (1) Assume that the length of the one-dimensional sliding window is xstep (xstep is a positive odd number). At any point (i, j) on the coherent profile, extend (xstep-1) / 2 grid points to both sides to extract the coherent attribute value. The coherent value of the point (i, j) is taken to the left to obtain the vector L; the coherent value of the point (i, j) is taken to the right to obtain the vector R. The gradients of L and R are calculated respectively to obtain ZL and ZR.

[0113] (2) Calculate the measure 1 that characterizes the fracture: sum the positive and negative values ​​contained in ZL and take the absolute value to obtain ZL_P_sum and ZL_N_sum; sum the positive and negative values ​​contained in ZR and take the absolute value to obtain ZR_P_sum and ZR_N_sum;

[0114] (3) Calculate the metric 2 that characterizes the fracture: Count whether ZL and ZR are all positive or all negative. If they are all positive or all negative, the value is 1; otherwise, the value is 0, that is,

[0115] If ZL are all positive numbers, set the parameter ZLZ1=1;

[0116] If ZL is all negative, set parameter ZLZ2=1;

[0117] If ZR is all positive, set parameter ZRZ1=1;

[0118] If ZR is all negative, set parameter ZRZ2=1.

[0119] (4) Based on the two measures of fracture characterization given in (2) and (3), give the point a fracture indicator value according to the following model:

[0120]

[0121] Among them, “==”, “>”, and “<” are relational operators, which respectively mean “equal to”, “greater than”, and “less than”; “&&” is a logical AND operation, which means that the expressions on both sides must be satisfied at the same time for them to be valid.

[0122] Regarding the value of the indicator factor, the situation in different work areas cannot usually be generalized. The above reference values ​​given in this patent application are applicable to the situation where the general seismic data conditions are stable.

[0123] In the horizontal direction, a sliding window is used to make a forward judgment first, and then a reverse judgment. The maximum value of the two results is taken, and the fracture indication factor value is assigned point by point to obtain the horizontal judgment result, which is as follows:

[0124] (1) Take any relevant attribute profile data C M×N , where M is the number of longitudinal time samples and N is the number of transverse spatial samples;

[0125] Please refer to Figure 6 , shows a seismic coherent attribute profile image before the implementation of the fracture detection method in an embodiment of this specification. A horizontal line determination can be performed on the profile.

[0126] (2) Set the sliding window length to xstep (xstep is an odd number), slide along the coherent attribute profile from left to right and from top to bottom, and extract the coherent value each time to form a decision unit vector c i , where i = 1, 2, …, xstep;

[0127] (3) According to the judgment criteria in 2, i Make a judgment and assign the corresponding fracture indication factor.

[0128] (4) Iterate (2) and (3) until the entire C is completed. M×N The first horizontal determination is completed and the result is recorded as F 左 ;

[0129] (5) The data matrix C of the profile M×N Flip left and right, repeat steps (2) and (3) in claim 3 until the entire section is completed and the second lateral determination is completed. The result is recorded as F 右 :

[0130] (6) Extract F point by point 左 、F 右 The maximum value of the fracture indicator factor is taken as the fracture indicator factor of the point, denoted as F 横向 , the formula is as follows:

[0131] F 横向IJ =max I∈[1,M],J∈[1,N] {F 左IJ ,F 右IJ}, where I and J are the horizontal and vertical coordinate values, such as Figure 7 As shown, the lateral determination result after the fracture detection method in one embodiment of this specification is implemented is shown.

[0132] Based on the transverse determination results, the determination criteria for the fracture in the longitudinal direction are constructed. The fracture indication factor value is assigned point by point in a sliding window manner to determine the longitudinal ductility of the fracture, as follows:

[0133] (1) With F 横向 As a basis, set the sliding window length to tstep (tstep is an odd number), slide along the coherent attribute profile from top to bottom and from left to right, and extract the coherent value each time to form a judgment unit vector. i 、F i+1 Two vectors, where i = 1, 2, …, tstep;

[0134] (2) Constructing vertical judgment criteria for F i 、F i+1 Make a judgment and assign the corresponding fracture indication factor. Assume that F i 、F i+1 The number of non-zero elements in is N i 、Ni +1 , for any point i, its fracture indicator factor FAULT i It can be expressed as follows:

[0135]

[0136] Among them, max is the maximum value operation of vector elements, and mean is the mean operation of vector orientation.

[0137] like Figure 8 As shown, it shows the longitudinal determination result after the fracture detection method in one embodiment of this specification is implemented. Figures 9 and 10 , respectively showing the coherent attribute time slice before the implementation of the fracture detection method in an embodiment of this specification and the high-resolution fracture feature slice after the implementation.

[0138] The above process is repeated by iterating the section by section to obtain three-dimensional fracture detection results.

[0139] Based on the same inventive concept, a fracture detection device is also provided in the embodiments of this specification, as described in the following embodiments. Since the principle of solving the problem by the fracture detection device is similar to that of the fracture detection method, the implementation of the fracture detection device can refer to the implementation of the fracture detection method, and the repeated parts will not be repeated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived. Fig.11 is a structural block diagram of a fracture detection device according to an embodiment of this specification, such as Fig.11 As shown, it includes: an acquisition module 1101, a first determination module 1102, a second determination module 1103 and a generation module 1104. The structure is described below.

[0140] The acquisition module 1101 is used to acquire a three-dimensional seismic coherent attribute volume of a target area; extract multiple two-dimensional sections from the three-dimensional seismic coherent attribute volume; and calculate the inverse of the seismic coherent attribute corresponding to each of the multiple two-dimensional sections.

[0141] The first determination module 1102 is used to determine the transverse fracture indicator factor data corresponding to each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section by using the preset transverse determination rule.

[0142] The second determination module 1103 is used to construct a longitudinal determination rule based on the transverse fracture indication factor data corresponding to each two-dimensional section; and is also used to determine the longitudinal fracture indication factor data corresponding to each two-dimensional section using the longitudinal determination rule.

[0143] The generating module 1104 is used to generate the fracture detection result corresponding to the three-dimensional seismic coherent attribute volume based on the longitudinal fracture indicator factor data corresponding to each two-dimensional profile.

[0144] In some embodiments of the present specification, the preset lateral determination rule may be constructed in the following manner: setting a plurality of modes according to the lateral distribution morphology of the fracture on the inverse of the seismic coherence attribute; and setting a corresponding fracture indication factor value for each of the plurality of modes.

[0145] In some embodiments of the present specification, the first determination module can be specifically used to: determine the length of the transverse sliding window according to the horizontal fault distance of the fracture in the target area in the transverse direction; for each of the multiple points on each two-dimensional section, use the transverse sliding window length to extend to the left and right sides of each point to extract the inverse value of the seismic coherent attribute in the transverse window, and obtain the first vector and the second vector corresponding to each point on each two-dimensional section; determine the transverse fracture indicator factor value corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section; obtain the transverse fracture indicator factor data corresponding to each two-dimensional section according to the transverse fracture indicator factor value corresponding to each point on each two-dimensional section.

[0146] In some embodiments of the present specification, determining the numerical value of the transverse fracture indication factor corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section may include: calculating the gradients of the first vector and the second vector corresponding to each point on each two-dimensional section to obtain the first vector gradient and the second vector gradient; obtaining the first fracture measure corresponding to each point on each two-dimensional section according to the absolute value of the sum of the positive values ​​and the sum of the negative values ​​contained in the first vector gradient and the second vector gradient; obtaining the second fracture measure corresponding to each point on each two-dimensional section according to whether the first vector gradient and the second vector gradient are all positive numbers or all negative numbers; determining the numerical value of the transverse fracture indication factor corresponding to each point on each two-dimensional section based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section.

[0147] In some embodiments of the present specification, for each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the transverse window, which can include: using the transverse sliding window length to slide from left to right and from top to bottom along each of the two-dimensional sections, and extending to the left and right sides of each point to extract the inverse value of the seismic coherence attribute in the transverse window; accordingly, based on the first fracture measurement and the second fracture measurement corresponding to each point on each of the two-dimensional sections, the first transverse fracture indicator factor value corresponding to each point on each of the two-dimensional sections is determined.

[0148] In some embodiments of the present specification, the method may further include: utilizing the length of the transverse sliding window to slide along each two-dimensional section from right to left and from top to bottom, and extending to the left and right sides of each point respectively to extract the inverse value of the seismic coherent attribute within the transverse window; correspondingly, based on the first fracture measurement and the second fracture measurement corresponding to each point on each two-dimensional section, determining the second transverse fracture indicator factor value corresponding to each point on each two-dimensional section; the first determination module may further be specifically used to: use the larger value of the first transverse fracture indicator factor value corresponding to each point on each two-dimensional section and the second transverse fracture indicator factor value corresponding to each point on each two-dimensional section as the transverse fracture indicator factor value corresponding to each point on each two-dimensional section.

[0149] In some embodiments of the present specification, the second determination module can be specifically used to: determine the longitudinal extension depth of the fracture in the target area according to the transverse fracture indicator factor data corresponding to each two-dimensional profile; determine the longitudinal sliding window length according to the longitudinal extension depth of the fracture in the target area; and construct a longitudinal judgment rule based on the longitudinal sliding window length.

[0150] In some embodiments of the present specification, the second determination module can be specifically used to: utilize the length of the longitudinal sliding window to slide along each of the two-dimensional sections from top to bottom and from left to right, extract the transverse fracture indicator factor data within the longitudinal window each time sliding, obtain the longitudinal vector, and determine the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections according to the number of non-zero elements in two adjacent longitudinal vectors.

[0151] From the above description, it can be seen that the embodiments of this specification achieve the following technical effects: taking the typical transverse morphological features of the fracture in the coherent attributes as a reference, enumerating the patterns of the discontinuous features of the coherent attributes, forming transverse determination criteria, thereby matching with the measured coherent attributes and assigning corresponding fracture indication factor values, and ensuring the accuracy of the fracture position by taking the best of the forward and reverse determination results. In addition, considering the longitudinal ductility of the fracture, on the basis of the transverse determination results, a longitudinal determination criterion is set to determine the longitudinal ductility of the fracture, and by combining the transverse determination with the longitudinal determination, a more accurate fracture detection result of the measured data is finally achieved.

[0152] This specification also provides a computer device, which can be found in Fig.12The computer device structure diagram of the fracture detection method provided by the embodiment of this specification is shown, and the computer device may specifically include an input device 121, a processor 122, and a memory 123. Among them, the memory 123 is used to store processor executable instructions. When the processor 122 executes the instructions, the steps of the fracture detection method described in any of the above embodiments are implemented.

[0153] In this embodiment, the input device may specifically be one of the main devices for information exchange between the user and the computer system. The input device may include a keyboard, a mouse, a camera, a scanner, a light pen, a handwriting input board, a voice input device, etc.; the input device is used to input the original data and the program for processing these numbers into the computer. The input device can also obtain and receive data transmitted from other modules, units, and devices. The processor can be implemented in any appropriate manner. For example, the processor can take the form of a computer-readable medium, a logic gate, a switch, an application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), a programmable logic controller, and an embedded microcontroller, etc., such as a microprocessor or a processor and a computer-readable program code (such as software or firmware) that can be executed by the (micro) processor. The memory may specifically be a memory device used to store information in modern information technology. The memory may include multiple levels. In a digital system, anything that can store binary data can be a memory; in an integrated circuit, a circuit with a storage function without a physical form is also called a memory, such as a RAM, a FIFO, etc.; in a system, a storage device with a physical form is also called a memory, such as a memory stick, a TF card, etc.

[0154] In this embodiment, the functions and effects specifically realized by the computer device can be explained in comparison with other embodiments and will not be described in detail here.

[0155] A computer storage medium based on the fracture detection method is also provided in the embodiments of the present specification. The computer storage medium stores computer program instructions, and when the computer program instructions are executed, the steps of the fracture detection method described in any of the above embodiments are implemented.

[0156] In this embodiment, the storage medium includes, but is not limited to, a random access memory (RAM), a read-only memory (ROM), a cache, a hard disk (HDD) or a memory card. The memory may be used to store computer program instructions. The network communication unit may be an interface for network connection communication set in accordance with the standard specified by the communication protocol.

[0157] In this embodiment, the functions and effects specifically implemented by the program instructions stored in the computer storage medium can be explained in comparison with other embodiments and will not be repeated here.

[0158] Obviously, those skilled in the art should understand that the modules or steps of the above-mentioned embodiments of this specification can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, and optionally, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order from that here, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. In this way, the embodiments of this specification are not limited to any specific combination of hardware and software.

[0159] It should be understood that the above description is for illustration and not for limitation. Many embodiments and many applications beyond the examples provided will be apparent to those skilled in the art upon reading the above description. Therefore, the scope of this specification should not be determined with reference to the above description, but should be determined with reference to the preceding claims and the full scope of equivalents to which these claims belong.

[0160] The above description is only the preferred embodiment of this specification and is not intended to limit this specification. For those skilled in the art, the embodiments of this specification may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this specification shall be included in the protection scope of this specification.

Claims

1. A fracture detection method, characterized in that: include: Acquire a three-dimensional seismic coherent attribute volume of the target area; extract multiple two-dimensional sections from the three-dimensional seismic coherent attribute volume; calculate the inverse of the seismic coherent attribute corresponding to each two-dimensional section in the multiple two-dimensional sections; Determine the transverse fracture indicator factor data corresponding to each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section by using the preset transverse determination rule; Based on the transverse fracture indicator factor data corresponding to each two-dimensional section, a longitudinal determination rule is constructed; and by using the longitudinal determination rule, the longitudinal fracture indicator factor data corresponding to each two-dimensional section is determined; Based on the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections, a fracture detection result corresponding to the three-dimensional seismic coherent attribute volume is generated.

2. The fracture detection method according to claim 1, characterized in that: The preset horizontal determination rule is constructed in the following manner: Set up multiple modes according to the lateral distribution morphology of the faults on the inverse of the seismic coherence attributes; A corresponding fracture indication factor value is set for each of the multiple modes.

3. The fracture detection method according to claim 1, characterized in that: Using a preset transverse determination rule, according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section, determining the transverse fracture indicator factor data corresponding to each two-dimensional section, including: Determine the length of the lateral sliding window according to the horizontal fault distance of the fracture in the target area in the lateral direction; For each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each of the points to extract the inverse value of the seismic coherence attribute in the transverse window, so as to obtain the first vector and the second vector corresponding to each of the points on each of the two-dimensional sections; Determine, according to the first vector and the second vector corresponding to each point on each two-dimensional section, a value of a transverse fracture indication factor corresponding to each point on each two-dimensional section; According to the transverse fracture indication factor values ​​corresponding to each point on each two-dimensional section, the transverse fracture indication factor data corresponding to each two-dimensional section are obtained.

4. The fracture detection method according to claim 3, characterized in that: Determining the transverse fracture indication factor value corresponding to each point on each two-dimensional section according to the first vector and the second vector corresponding to each point on each two-dimensional section includes: Calculating the gradients of the first vector and the second vector corresponding to each point on each of the two-dimensional sections to obtain the first vector gradient and the second vector gradient; Obtaining a first fracture measure corresponding to each point on each of the two-dimensional sections according to the absolute value of the sum of the positive values ​​and the sum of the negative values ​​included in the first vector gradient and the second vector gradient; Obtaining a second fracture measure corresponding to each point on each of the two-dimensional sections according to whether the first vector gradient and the second vector gradient are all positive numbers or all negative numbers; Based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a transverse fracture indication factor value corresponding to each point on each two-dimensional section is determined.

5. The fracture detection method according to claim 4, characterized in that: For each of the multiple points on each of the two-dimensional sections, the transverse sliding window length is used to extend to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window, including: Using the length of the transverse sliding window, slide along each two-dimensional section from left to right and from top to bottom, respectively extending to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window; Correspondingly, based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a first transverse fracture indication factor value corresponding to each point on each two-dimensional section is determined.

6. The fracture detection method according to claim 5, characterized in that: Also includes: Using the length of the transverse sliding window, slide along each two-dimensional section from right to left and from top to bottom, respectively extending to the left and right sides of each point to extract the inverse value of the seismic coherent attribute within the transverse window; Accordingly, based on the first fracture measure and the second fracture measure corresponding to each point on each two-dimensional section, a second transverse fracture indicator factor value corresponding to each point on each two-dimensional section is determined; The method further comprises: The larger value of the first transverse fracture indication factor value corresponding to each point on each two-dimensional section and the second transverse fracture indication factor value corresponding to each point on each two-dimensional section is used as the transverse fracture indication factor value corresponding to each point on each two-dimensional section.

7. The fracture detection method according to claim 1, characterized in that: Based on the transverse fracture indicator factor data corresponding to each two-dimensional section, a longitudinal determination rule is constructed, including: Determining the longitudinal extension depth of the fracture in the target area according to the transverse fracture indicator factor data corresponding to each two-dimensional section; Determining the length of the longitudinal sliding window according to the longitudinal extension depth of the fracture in the target area; Based on the length of the vertical sliding window, a vertical decision rule is constructed.

8. The fracture detection method according to claim 7, characterized in that: Determining the longitudinal fracture indicator factor data corresponding to each of the two-dimensional sections by using the longitudinal determination rule includes: The longitudinal sliding window length is utilized to slide along each two-dimensional section from top to bottom and from left to right. Each time the window slides, the transverse fracture indicator factor data in the longitudinal window is extracted to obtain a longitudinal vector. The longitudinal fracture indicator factor data corresponding to each two-dimensional section is determined based on the number of non-zero elements in two adjacent longitudinal vectors.

9. A fracture detection device, characterized in that: include: An acquisition module, used for acquiring a three-dimensional seismic coherent attribute volume of a target area; Extracting a plurality of two-dimensional sections from the three-dimensional seismic coherence attribute volume; Calculating the inverse of the seismic coherence attribute corresponding to each of the multiple two-dimensional sections; A first determination module is used to determine the transverse fracture indicator factor data corresponding to each two-dimensional section according to the inverse of the seismic coherence attribute corresponding to each two-dimensional section by using a preset transverse determination rule; The second determination module is used to construct a longitudinal determination rule based on the transverse fracture indication factor data corresponding to each of the two-dimensional sections; and is also used to determine the longitudinal fracture indication factor data corresponding to each of the two-dimensional sections using the longitudinal determination rule; A generating module is used to generate a fracture detection result corresponding to the three-dimensional seismic coherent attribute volume based on the longitudinal fracture indicator factor data corresponding to each two-dimensional profile.

10. A computer-readable storage medium having computer instructions stored thereon, characterized in that: When the instructions are executed by a processor, the steps of the method according to any one of claims 1 to 8 are implemented.

Images