Earthquake identification method and system for internal structure of small-scale low-sequence fault zone

Abstract

The invention belongs to the technical field of oil and gas field exploration and development, and discloses an earthquake identification method and system for an internal structure of a small-scale low-sequence fault zone, and the method comprises the steps: determining an earthquake frequency division body of a target interval according to geological data; determining a minor fault main component tuning frequency division body and a minor fault background frequency division body according to the seismic frequency division body of the target interval; carrying out breakpoint sharpening processing on the minor fault main component tuning frequency division body to obtain seismic amplitude data after breakpoint sharpening processing; merging the minor fault background frequency division volume and the seismic amplitude data after the breakpoint sharpening processing to obtain a fused data volume; and determining the internal structure of the fault zone according to the fused data volume. According to the method, factors of an internal structure of the small fault zone and seismic reflection characteristics of the small fault zone are considered, breakpoint sharpening processing is introduced for the first time, and the continuous characteristics of small fractures in the longitudinal direction are enhanced. According to the method, the recognition precision of small fault recognition on seismic data can be improved, and the reflected small fault structure is more deterministic.

Description

Technical Field

[0001] The invention belongs to the technical field of oil and gas field exploration and development, and particularly relates to a seismic identification method and system for the internal structure of a small-scale low-order fault zone. Background Art

[0002] At present, the reservoir types in some exploration areas are mainly block oil and gas reservoirs. The oil and gas exploration in the fault lake basin has reached the middle and late stages of exploration. The large-scale structural oil and gas reservoirs have been basically drilled out, and the exploration of the remaining small fault block oil and gas reservoirs is becoming increasingly difficult. The identification of low-order small faults and their internal structures has become a major challenge.

[0003] The existing technology uses the amplitude anomaly of the seismic data body and the morphological misalignment of the travel-time event axis to identify small faults, but does not pay attention to the breakpoint sharpening characteristics based on the frequency response of small faults. Summary of the invention

[0004] In view of the above problems, the present invention provides a seismic identification method and system for the internal structure of a small-scale low-order fault zone, which adopts the following technical solutions:

[0005] The present invention provides a seismic identification method for the internal structure of a small-scale low-order fault zone, comprising the following steps: determining a seismic frequency-fraction body of a target layer segment according to geological data; determining a tuning frequency-fraction body of a main component of a small fault and a background frequency-fraction body of a small fault according to the seismic frequency-fraction body of the target layer segment; performing breakpoint sharpening processing on the tuning frequency-fraction body of the main component of the small fault to obtain seismic amplitude data after the breakpoint sharpening processing; merging the background frequency-fraction body of the small fault and the seismic amplitude data after the breakpoint sharpening processing to obtain a fused data body; and determining the internal structure of the fault zone according to the fused data body.

[0006] Further, determining the seismic frequency division volume of the target layer segment according to the geological data includes the following steps:

[0007] The seismic frequency volume of the target layer is calculated according to the geological data, and the seismic frequency volume is obtained by Fourier transform.

[0008] Further, according to the seismic frequency fraction of the target layer segment, the tuning frequency fraction of the main component of the small fault and the background frequency fraction of the small fault are determined, including the following steps:

[0009] Compare the coherent response characteristics of small faults in each main component tuning frequency division body under the same coherence parameter, and take the frequency range with the best performance of small fault characteristics as the main component tuning frequency division body of small fault;

[0010] The seismic data volume of the tuning frequency volume of the small fault main component is screened out from the seismic amplitude data of the seismic frequency volume corresponding to the target layer segment to obtain the small fault background frequency volume.

[0011] Further, the small fault main component tuning frequency division body is subjected to breakpoint sharpening processing to obtain seismic amplitude data after breakpoint sharpening processing, including the following steps:

[0012] The seismic attributes transformed by Laplacian operator template are used to process the tuning frequency-dividing volume of the main component of the small fault, and the seismic amplitude data after Laplacian sharpening are obtained;

[0013] According to the longitudinal sharpening factor and the seismic amplitude data after Laplacian sharpening, the seismic amplitude data after breakpoint sharpening is determined.

[0014] Further, the small fault background frequency-divided volume and the seismic amplitude data after breakpoint sharpening are merged to obtain a fused data volume, including the following steps:

[0015] The small fault background frequency-dividing body and the seismic amplitude data after breakpoint sharpening are weighted to obtain a fused data body, wherein the weight coefficient of the seismic amplitude data after breakpoint sharpening is not less than 60%.

[0016] Further, the internal structure of the fault zone is determined based on the fused data volume, including the following steps:

[0017] The third generation coherence calculation is carried out on the fused data body, in which the aperture in the coherence parameter takes a 9-point operator, and the internal structure of the fault zone is determined based on the planar combination characteristics of small faults on the coherence results.

[0018] Furthermore, geological data include petroleum seismic data, oil and gas well logging data and small fault interpretation data.

[0019] The present invention also provides a seismic identification system for the internal structure of a small-scale low-order fault zone, comprising:

[0020] The first calculation module is used to determine the seismic frequency division volume of the target layer segment according to the geological data;

[0021] The second calculation module is used to determine the tuning frequency body of the minor fault main component and the minor fault background frequency body according to the seismic frequency body of the target layer segment;

[0022] The third calculation module is used to perform breakpoint sharpening processing on the tuning frequency division body of the main component of the small fault to obtain seismic amplitude data after breakpoint sharpening processing;

[0023] The fourth calculation module is used to merge the small fault background frequency-divided volume and the seismic amplitude data after breakpoint sharpening to obtain a fused data volume;

[0024] The fifth calculation module is used to determine the internal structure of the fault zone based on the fused data volume.

[0025] Furthermore, the second calculation module is specifically used for:

[0026] Compare the coherent response characteristics of small faults in each main component tuning frequency division body under the same coherence parameter, and take the frequency range with the best performance of small fault characteristics as the main component tuning frequency division body of small fault;

[0027] The seismic data volume of the tuning frequency volume of the small fault main component is screened out from the seismic amplitude data of the seismic frequency volume corresponding to the target layer segment to obtain the small fault background frequency volume.

[0028] Furthermore, the third computing module is specifically used for:

[0029] The seismic attributes transformed by Laplacian operator template are used to process the tuning frequency-dividing volume of the main component of the small fault, and the seismic amplitude data after Laplacian sharpening are obtained;

[0030] According to the longitudinal sharpening factor and the seismic amplitude data after Laplacian sharpening, the seismic amplitude data after breakpoint sharpening is determined.

[0031] Furthermore, the fourth calculation module is specifically used for:

[0032] The small fault background frequency-dividing body and the seismic amplitude data after breakpoint sharpening are weighted to obtain a fused data body, wherein the weight coefficient of the seismic amplitude data after breakpoint sharpening is not less than 60%.

[0033] Beneficial effects of the present invention:

[0034] 1. The present invention takes into account the internal structure of the small fault zone and its seismic reflection characteristics, introduces breakpoint sharpening processing for the first time, strengthens the vertical continuity characteristics of the small fault, and fuses it with the background data in a weighted manner to form a thematic data body reflecting the characteristics of the small fault.

[0035] 2. The present invention can improve the recognition accuracy of small faults in seismic data, and the small fault structures reflected are more certain, providing a new method for the interpretation of small faults. The evaluation effect brought by this method is more in line with the actual geological laws, so as to reduce the risk of fault closure drilling, and also provide a geological evaluation reference for optimizing the development plan of block oil and gas reservoirs in the development stage.

[0036] Other features and advantages of the present invention will be described in the following description, and partly become obvious from the description, or be understood by implementing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures pointed out in the description and the drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] 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 use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0038] Figure 1 A schematic flow chart of a method for seismic identification of the internal structure of a small-scale low-order fault zone according to an embodiment of the present invention is shown;

[0039] Figure 2 A schematic structural diagram of a seismic identification system for the internal structure of a small-scale low-order fault zone according to an embodiment of the present invention is shown;

[0040] Figure 3a A schematic diagram of the coherence results of the original processed data of a certain exploration area is shown;

[0041] Figure 3b A schematic diagram of coherent results after the seismic identification method and system for the internal structure of a small-scale low-order fault zone according to an embodiment of the present invention is applied to a certain exploration area is shown. DETAILED DESCRIPTION

[0042] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0043] The present invention provides a seismic identification method and system for the internal structure of a small-scale low-order fault zone, which utilizes the frequency response characteristics and breakpoint sharpening characteristics of small faults to solve the problem of identifying small-scale low-order faults in seismic data.

[0044] like Figure 1 As shown, a seismic identification method for the internal structure of a small-scale low-order fault zone comprises the following steps:

[0045] S1. Obtain geological data, including petroleum seismic data, oil and gas well logging data and fault interpretation data.

[0046] S2. Determining the seismic frequency volume of the target layer segment according to the geological data includes: calculating the seismic frequency volume of the target layer segment according to the geological data, and obtaining the seismic frequency volume by Fourier transform.

[0047] S3, according to the seismic frequency fraction of the target layer segment, determining the tuning frequency fraction of the main component of the small fault and the background frequency fraction of the small fault, including the following steps:

[0048] S31. Compare the coherent response characteristics of the small fault in each main component tuning frequency division body under the same coherence parameter, and take the frequency range with the best performance of the small fault characteristics as the small fault main component tuning frequency division body, denoted as A(x,y).

[0049] S32. Filter out the seismic data volume of the tuning frequency volume of the small fault main component from the seismic amplitude data of the seismic frequency volume corresponding to the target layer segment, and obtain the small fault background frequency volume, recorded as A1(x,y), as the background.

[0050] S4, performing breakpoint sharpening processing on the tuning frequency division body of the main component of the small fault to obtain seismic amplitude data after breakpoint sharpening processing, including the following steps:

[0051] S41. The seismic attributes converted by the Laplacian operator template are used to process the tuning frequency-dividing body of the main component of the small fault to obtain the seismic amplitude data after Laplacian sharpening, as follows:

[0052]

[0053] In the formula, A E (x, y) represents the seismic amplitude data after Laplacian sharpening, A(x, y) is the tuning frequency-dividing body of the main component of the small fault, dimensionless; Represents the seismic attributes transformed by the Laplacian operator template.

[0054] In this step, the image contrast is enhanced by sharpening processing to make the edges, contours and details clearer. Usually, the corresponding sharpening operation is performed for the cause of image blur, which belongs to the content of image restoration. The essence of image blur is that the image is caused by averaging or integration operations, so the image can be restored by operations such as differential operations to make the image clear. From the perspective of spectrum, the essence of image blur is that its high-frequency components are attenuated, so the image can be clarified by high-pass filtering. However, it should be noted that the image that can be sharpened must have a high signal-to-noise ratio, otherwise the signal-to-noise ratio of the image after sharpening will be lower, so that the noise increases more than the signal, so it is generally necessary to remove or reduce the noise before sharpening.

[0055] The basic idea of ​​sharpening by Laplacian algorithm in this step is that when the grayscale of the central pixel of the neighborhood is lower than the average grayscale of other pixels in its neighborhood, the grayscale of this central pixel should be further reduced, and when it is higher, the grayscale of the central pixel should be further increased, thereby achieving image sharpening.

[0056] The specific process of the Laplacian algorithm is: the grayscale value of the image is 9 times its own grayscale value minus the grayscale value of the surrounding pixels. If a bright spot appears in a dark area, the result of the sharpening process is that the bright spot becomes brighter, increasing the noise of the image. This is because the edges in the image are those areas where the grayscale jumps.

[0057] S42, according to the longitudinal sharpening factor and the seismic amplitude data after Laplacian sharpening, determine the seismic amplitude data after breakpoint sharpening, as follows:

[0058] A F (x,y)=K*A E (x, y) (2)

[0059] In the formula, A F (x, y) is the seismic amplitude data after breakpoint sharpening, dimensionless; A E (x, y) is the seismic amplitude data after the commonly used Laplacian sharpening process, which is dimensionless; k is the longitudinal sharpening factor, which is a transformation template related to the sharpening diffusion effect. This parameter aims to enhance the continuous characteristics of small faults in the vertical direction.

[0060] Among them, the vertical sharpening factor k is as follows:

[0061]

[0062] S5, merging the small fault background frequency-dividing volume and the seismic amplitude data after breakpoint sharpening processing to obtain a fused data volume, including the following steps:

[0063] The background frequency-divided volume of the small fault and the seismic amplitude data after breakpoint sharpening are weighted to obtain a fused data volume, wherein the weight coefficient of the seismic amplitude data after breakpoint sharpening is at least 60%, and a fused data volume M(x,y) with the small fault as the theme feature is obtained;

[0064] M(x,y)=w1*A F (x,y)+w2*A1(x,y) (4)

[0065] Among them, W1+W2=1, where W1 represents the weight coefficient of the seismic amplitude data after breakpoint sharpening processing; W2 represents the weight coefficient of the small fault background frequency division body.

[0066] In the formula, A F (x, y) is the seismic amplitude data after breakpoint sharpening; A1(x, y) is the small fault background frequency division body.

[0067] S6. Determine the internal structure of the fault zone according to the fused data volume, including the following steps:

[0068] The publicly available third-generation coherence calculation is performed on the fused data body, in which the aperture in the coherence parameter takes a 9-point operator, and the internal structure of the fault zone is determined based on the planar combination characteristics of small faults on the coherence results.

[0069] The third-generation coherence calculation method was proposed by Gersztenkorn and Marfurt (1999) and applied to the eigenvalue method of seismic coherence analysis. The third-generation coherence calculation method calculates the eigenvalues ​​of the covariance matrix and takes the ratio of the maximum eigenvalue to the matrix trace as the coherence value.

[0070] Usually for the covariance matrix, the eigenvectors and eigenvalues ​​of the matrix are generally arranged in descending order according to their ability to express the energy of the data within the analysis aperture and the main changes in the matrix. The first eigenvalues ​​and their eigenvectors contain most of the information. The third-generation coherent volume technology only needs the trace and the first eigenvalue of the covariance matrix to estimate the characteristic structure coherence.

[0071] In this step, the third-generation coherence volume algorithm is used to convert the seismic data into coherence volume data that is more beneficial for crack and fault detection. It has better detection capabilities for small cracks, making the identification results more accurate.

[0072] Based on the above-mentioned earthquake identification method of the internal structure of small-scale low-order fault zones, such as Figure 2 As shown, the present invention also provides a seismic identification system for the internal structure of a small-scale low-order fault zone, comprising a first calculation module, a second calculation module, a third calculation module, a fourth calculation module and a fifth calculation module.

[0073] Among them, the first calculation module is used to determine the seismic frequency volume of the target layer segment according to the geological data; the second calculation module is used to determine the small fault main component tuning frequency volume and the small fault background frequency volume according to the seismic frequency volume of the target layer segment; the third calculation module is used to perform breakpoint sharpening processing on the small fault main component tuning frequency volume to obtain the seismic amplitude data after breakpoint sharpening processing; the fourth calculation module is used to merge the small fault background frequency volume and the seismic amplitude data after breakpoint sharpening processing to obtain a fused data volume; the fifth calculation module is used to determine the internal structure of the fault zone according to the fused data volume.

[0074] Taking a small fault in a certain exploration area as an example, the implementation process and effect of the present invention are described. Figure 3b The middle figure shows the distribution of small fault planes after applying the identification method of the present invention. It is found that there are obvious changes within the same fault zone, which reflects the difference in its internal structure and the segmented growth of the fault. Figure 3b In the circle shown in the middle, Figure 3aCompared with the middle figure, the number of developed small fractures is greater, reflecting that more small fractures can be identified after being processed by the identification method of the present invention, reflecting the technical advantages of the present invention. Figure 3b The figure shows the characteristics of the common development of small faults in the study area, which is consistent with the geological fact that small faults are common in the Tongbo Temple Formation in this area, indicating that the application effect of the present invention conforms to geological laws and is reliable.

[0075] The present invention takes into account the internal structure of small fault zones and their seismic reflection characteristics, introduces breakpoint sharpening processing for the first time, strengthens the vertical continuity characteristics of small faults, and fuses them with background data in a weighted manner to form a thematic data body reflecting the characteristics of small faults.

[0076] The present invention can improve the recognition accuracy of small faults on seismic data, and the reflected small fault structures are more certain, providing a new method for the interpretation of small faults. The evaluation effect brought by this method is more in line with the actual geological laws, so as to reduce the risk of fault closure drilling, and also provide a geological evaluation reference for optimizing the development plan of block oil and gas reservoirs in the development stage.

[0077] Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent substitutions for some of the technical features therein; and these modifications or substitutions do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A seismic identification method for the internal structure of a small-scale low-order fault zone, characterized in that: The following steps are involved: Determine the seismic frequency distribution of the target layer according to geological data; According to the seismic frequency division body of the target layer segment, the tuning frequency division body of the main component of the small fault and the background frequency division body of the small fault are determined; Perform breakpoint sharpening on the tuning frequency-dividing body of the main component of the small fault to obtain seismic amplitude data after breakpoint sharpening; The small fault background frequency-divided volume and the seismic amplitude data after breakpoint sharpening are merged to obtain a fused data volume; The internal structure of the fault zone is determined based on the fused data volume.

2. The seismic identification method of the internal structure of a small-scale low-order fault zone according to claim 1 is characterized in that: Determining the seismic frequency division volume of the target layer segment according to geological data includes the following steps: The seismic frequency volume of the target layer is calculated according to the geological data, and the seismic frequency volume is obtained by Fourier transform.

3. The seismic identification method of the internal structure of a small-scale low-order fault zone according to claim 1, characterized in that: According to the seismic frequency fraction of the target layer segment, determining the tuning frequency fraction of the main component of the small fault and the background frequency fraction of the small fault includes the following steps: Compare the coherent response characteristics of small faults in each main component tuning frequency division body under the same coherence parameter, and take the frequency range with the best performance of small fault characteristics as the main component tuning frequency division body of small fault; The seismic data volume of the tuning frequency volume of the small fault main component is screened out from the seismic amplitude data of the seismic frequency volume corresponding to the target layer segment to obtain the small fault background frequency volume.

4. The seismic identification method of the internal structure of a small-scale low-order fault zone according to claim 1, characterized in that: Performing breakpoint sharpening processing on the tuning frequency division body of the main component of the small fault to obtain seismic amplitude data after breakpoint sharpening processing includes the following steps: The seismic attributes transformed by Laplacian operator template are used to process the tuning frequency-dividing volume of the main component of the small fault, and the seismic amplitude data after Laplacian sharpening are obtained; According to the longitudinal sharpening factor and the seismic amplitude data after Laplacian sharpening, the seismic amplitude data after breakpoint sharpening is determined.

5. The seismic identification method of the internal structure of a small-scale low-order fault zone according to claim 1, characterized in that: The small fault background frequency-dividing volume and the seismic amplitude data after breakpoint sharpening are merged to obtain a fused data volume, including the following steps: The small fault background frequency-dividing body and the seismic amplitude data after breakpoint sharpening are weighted to obtain a fused data body, wherein the weight coefficient of the seismic amplitude data after breakpoint sharpening is not less than 60%.

6. The seismic identification method of the internal structure of a small-scale low-order fault zone according to claim 1, characterized in that: Determining the internal structure of the fault zone based on the fused data volume includes the following steps: The third generation coherence calculation is carried out on the fused data body, in which the aperture in the coherence parameter takes a 9-point operator, and the internal structure of the fault zone is determined based on the planar combination characteristics of small faults on the coherence results.

7. The seismic identification method for the internal structure of a small-scale low-order fault zone according to any one of claims 1 to 6, characterized in that: Geological data include petroleum seismic data, oil and gas well logging data and small fault interpretation data.

8. A seismic identification system for the internal structure of a small-scale low-order fault zone, characterized in that: include: The first calculation module is used to determine the seismic frequency division volume of the target layer segment according to the geological data; The second calculation module is used to determine the tuning frequency body of the minor fault main component and the minor fault background frequency body according to the seismic frequency body of the target layer segment; The third calculation module is used to perform breakpoint sharpening processing on the tuning frequency division body of the main component of the small fault to obtain seismic amplitude data after breakpoint sharpening processing; The fourth calculation module is used to merge the small fault background frequency-divided volume and the seismic amplitude data after breakpoint sharpening to obtain a fused data volume; The fifth calculation module is used to determine the internal structure of the fault zone based on the fused data volume.

9. The seismic identification system for the internal structure of a small-scale low-order fault zone according to claim 8, characterized in that: The second calculation module is specifically used for: Compare the coherent response characteristics of small faults in each main component tuning frequency division body under the same coherence parameter, and take the frequency range with the best performance of small fault characteristics as the main component tuning frequency division body of small fault; The seismic data volume of the tuning frequency volume of the small fault main component is screened out from the seismic amplitude data of the seismic frequency volume corresponding to the target layer segment to obtain the small fault background frequency volume.

10. The seismic identification system for the internal structure of a small-scale low-order fault zone according to claim 8, characterized in that: The third calculation module is specifically used for: The seismic attributes transformed by Laplacian operator template are used to process the tuning frequency-dividing volume of the main component of the small fault, and the seismic amplitude data after Laplacian sharpening are obtained; According to the longitudinal sharpening factor and the seismic amplitude data after Laplacian sharpening, the seismic amplitude data after breakpoint sharpening is determined.

11. The seismic identification system for the internal structure of a small-scale low-order fault zone according to claim 8, characterized in that: The fourth calculation module is specifically used for: The small fault background frequency-dividing body and the seismic amplitude data after breakpoint sharpening are weighted to obtain a fused data body, wherein the weight coefficient of the seismic amplitude data after breakpoint sharpening is not less than 60%.

Images