Advanced azimuth angle optimization method based on binary coupling

By using a binary coupling method of seismic amplitude and travel time change rate of OVT gathers, the range of dominant azimuth angles is optimized, solving the problem of quantitatively selecting dominant azimuth angles in existing technologies and improving the accuracy of seismic imaging and reservoir prediction.

CN122072364APending Publication Date: 2026-05-22CHINA PETROLEUM & CHEMICAL CORP +1
View PDF 3 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-11-22
Publication Date
2026-05-22

AI Technical Summary

Technical Problem

Existing technologies lack quantitative methods to optimize the advantageous azimuth range of the OVT domain with a wide azimuth angle, resulting in poor seismic imaging performance and affecting the accuracy of structural interpretation and reservoir prediction.

Method used

Based on well logging and other data, the range of dominant azimuth angles is calculated by binary coupling of seismic amplitude and travel time rate of OVT gathers. The dominant azimuth angle is then selected by utilizing the variation characteristics of OVT gather amplitude and two-way travel time with azimuth.

Benefits of technology

It improves the accuracy of seismic interpretation and reservoir prediction, provides reliable guidance on the range of advantageous azimuth angles, and enhances the imaging effect of geological bodies.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN122072364A_ABST
    Figure CN122072364A_ABST
Patent Text Reader

Abstract

The invention provides a dominant azimuth angle optimization method based on binary coupling. The dominant azimuth angle optimization method based on binary coupling comprises the following steps: step 1, determining a seismic reflection axis corresponding to a target layer; 2, anisotropy analysis of an OVT gather target layer is carried out; step 3, extracting travel time of a target layer and a seismic amplitude value of the OVT gather; 4, calculating the seismic amplitude and travel time change rate of a target layer of the OVT gather; and step 5, coupling the OVT gather target layer seismic amplitude with the travel time change rate, and determining the dominant azimuth angle range. According to the dominant azimuth angle optimization method based on binary coupling, optimization of the dominant azimuth angle of the target stratum can be conveniently carried out by utilizing the change characteristics of the amplitude and travel time of the OVT gather data along with the azimuth, the azimuth-based superposition of the pre-stack seismic data of the exploratory area is guided, and high-quality seismic data is provided for the comprehensive research of seismic geology of the exploratory area.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of exploration geophysics, and in particular to a method for optimizing the preferred azimuth angle based on binary coupling. Background Technology

[0002] In recent years, with the continuous advancement of exploration technology and the development of OVT (Offset Vector Tile) technology, more and more OVT domain pre-stack gathering data with wide azimuth and large offset have emerged and been applied in production. Compared with conventional pre-stack gathering seismic data, OVT domain pre-stack gatherings retain more azimuth and offset information.

[0003] Information about subsurface geological bodies recorded at different azimuths and offsets varies during earthquake propagation. Therefore, high-density, wide-azimuth, and long-offset seismic data provide a solid data foundation for structural interpretation and reservoir prediction. Their wide-azimuth and wide-bandwidth characteristics offer new insights for detailed structural interpretation and reservoir prediction. Utilizing CRP gathers (OVGs) after OVT migration for seismic interpretation is an important direction for future seismic exploration.

[0004] By utilizing the travel time and velocity variation information obtained from the five-dimensional gathers in the OVT data domain, and selecting seismic attribute variations within a limited range of azimuth and shot-receiver distance, it is beneficial to focus on geological target layers, overcome problems such as energy instability, density imbalance, and low signal-to-noise ratio in the original pre-stack gathers, improve the accuracy of target layer tracking, determine specific geological structural features and their spatial distribution, display structural and reservoir development details, and improve the accuracy of structural interpretation and reservoir prediction.

[0005] When using OVT gathers for detailed fault characterization and reservoir description, it was found that different azimuth angles yielded different results, with a particular azimuth providing better seismic imaging and thus enabling better structural interpretation and reservoir prediction. This azimuth angle that enhances the imaging effect of geological bodies is usually referred to as the dominant azimuth. Qualitatively, the azimuth perpendicular to the geological body is considered the dominant azimuth when characterizing geological bodies, but there is still no quantitative method to optimize the range of specific dominant azimuth angles.

[0006] Chinese patent application CN202010529547.6 discloses a method and system for converting OVT domain data into azimuth domain imaging gathers. The method first calculates the Di x layer velocity and uses it to set the initial layer velocity in the depth domain space. An optimized objective function is then used to determine the optimized layer velocity in the depth domain space. Based on the optimized layer velocity in the depth domain space, the maximum reflection angle is calculated, and the reflection angle range is determined. The regions within the azimuth and reflection angle ranges are then meshed to obtain the azimuth domain mesh. Finally, based on the optimized layer velocity in the depth domain space and the root mean square velocity of the pre-stack time offset in the OVT domain imaging gather data, the reflection angle corresponding to each azimuth is calculated, resulting in the azimuth domain imaging gather. This invention, by solving an objective function, achieves the conversion of OVT domain data in laterally variable speed media into azimuth domain imaging gathers to improve the accuracy of AVA / AVAZ inversion.

[0007] Chinese patent application CN202311115158.9 discloses a method and system for accurately interpreting small faults based on OVT domain gathers, including: optimizing azimuth angles using fault points as constraints to obtain first, second, and third optimized sub-azimuth angles; determining the starting azimuth angle based on the survey line direction of the study area; performing sub-azimuth angle stacking starting from the starting azimuth angle to obtain first, second, and third optimized sub-azimuth angle stacked profiles, and extracting composite attribute slices from the first, second, and third optimized sub-azimuth angle stacked profiles; denoising the optimized sub-azimuth angle stacked profiles and performing comparative analysis based on the denoising results; performing comparative analysis on the optimized sub-azimuth angle stacked composite attribute slices; obtaining the first, second, and third optimized sub-azimuth angle stacked comprehensive interpretation based on the comparative analysis results of the stacked profiles and stacked comprehensive attribute slices; and accurately interpreting small faults based on the optimized sub-azimuth angle stacked comprehensive interpretation. This invention achieves the goal of accurately interpreting small faults in coal seams.

[0008] Chinese patent application CN202310930152.0 discloses a method for determining an azimuth gather partitioning scheme, comprising the following steps: determining multiple initial partitioning angles of seismic data according to a preset angle partitioning step size; for each initial partitioning angle, dividing the seismic data into multiple azimuth gather data according to the initial partitioning angle and preset partitioning rules, thereby obtaining multiple azimuth gather partitioning schemes to be evaluated; for each azimuth gather partitioning scheme to be evaluated, acquiring seismic characteristic data corresponding to each azimuth gather data, and determining an azimuth anisotropy fracture prediction accuracy evaluation factor based on the seismic characteristic data; and determining the target azimuth gather partitioning scheme based on the azimuth anisotropy fracture prediction accuracy evaluation factor corresponding to each azimuth gather partitioning scheme to be evaluated. This invention provides a quantitative evaluation method for seismic data azimuth gather partitioning schemes, achieving rapid and efficient determination of the optimal seismic data azimuth gather partitioning scheme, thereby improving the accuracy of fracture prediction.

[0009] The above-mentioned existing technologies are all quite different from the present invention and have failed to solve the technical problem we want to solve. Therefore, we have invented a new method for optimizing the advantageous azimuth angle based on binary coupling. Summary of the Invention

[0010] The purpose of this invention is to provide a method for optimizing the azimuth angle based on binary coupling, which is based on data such as well logging and well logging, makes full use of the differences in azimuth anisotropy of OVT domain gathers, optimizes the reliable range of azimuth angles, and lays a solid foundation for the superposition of advantageous azimuth angles.

[0011] The objective of this invention can be achieved through the following technical measures: a method for selecting the dominant azimuth angle based on binary coupling, which includes:

[0012] Step 1: Identify the seismic reflection axis corresponding to the target layer;

[0013] Step 2: Conduct anisotropy analysis of the target layer of the OVT gather;

[0014] Step 3: Extract the travel time and seismic amplitude values ​​of the target layer from the OVT gather as a function of azimuth;

[0015] Step 4: Calculate the seismic amplitude and travel-time variation rate of the target layer in the OVT gather;

[0016] Step 5: Couple the seismic amplitude and travel time variation rate of the target layer of the OVT gather to determine the range of dominant azimuth angles.

[0017] The objective of this invention can also be achieved through the following technical measures:

[0018] In step 1, synthetic record calibration of typical wells is carried out based on comprehensive seismic geological analysis to identify the seismic reflection axis corresponding to the target layer.

[0019] In step 1, the comprehensive seismic geological analysis data includes: well logging, well logging and drilling geological data, and seismic data.

[0020] In step 2, an azimuth anisotropy analysis of the seismic amplitude and travel time of the target layer in the actual work area OVT gather is carried out; it is confirmed that the target layer has azimuth anisotropy, and an azimuth anisotropy study is conducted.

[0021] In step 3, after arranging the OVT gathers by azimuth, the seismic amplitude and travel time curves of the target layer at typical well locations are extracted as a function of azimuth.

[0022] In step 4, when calculating the rate of change of seismic amplitude, the seismic amplitude calibrated by the single-well synthetic record or the average seismic amplitude of each direction of the target layer in the OVT gather can be selected as the reference value for calculation; when calculating the rate of change of travel time, the travel time calibrated by the single-well synthetic record or the average travel time of each direction of the target layer in the OVT gather can be selected as the reference value for calculation.

[0023] In step 4, the formula for calculating the rate of change of seismic amplitude and travel time with azimuth for the target layer of the OVT gather is as follows:

[0024] ΔA = (A′ - A) / A * 100

[0025] Where A′ is the amplitude value of the target layer extracted by the OVT gather, A is the seismic amplitude value at that point on the conventional seismic profile or the average amplitude of the target layer in each direction of the OVT gather, and ΔA is the seismic amplitude variation rate;

[0026] ΔT=(T′-T) / T*100,

[0027] Where T′ is the travel time of the target layer extracted by the OVT gather, T is the travel time of that point on the conventional seismic profile or the average travel time of each direction of the target layer in the OVT gather, and ΔT is the rate of change of travel time.

[0028] In step 5, the coupling relationship value F between the seismic amplitude and travel time rate of change of the OVT gather is calculated. The range of azimuth angles with small and stable travel time and amplitude rate of change, and large amplitude value is identified as the dominant azimuth angle range.

[0029] In step 5, the formula for calculating the coupling relationship value F between the seismic amplitude and the rate of change of travel time of the target layer in the OVT gather is as follows:

[0030] F = ΔA / |ΔT|,

[0031] Where ΔA is the rate of change of earthquake amplitude and ΔT is the rate of change of travel time. The range of dominant F values ​​is defined according to the reflection characteristics of the target layer. The corresponding azimuth range is selected based on the dominant F value range, which is the dominant azimuth range.

[0032] The objective of this invention can also be achieved through the following technical measures: a binary coupling-based dominant azimuth angle optimization system, which uses a binary coupling-based dominant azimuth angle optimization method to optimize the dominant azimuth angle of the target layer by utilizing the amplitude and two-way travel time variation characteristics of OVT gather data with azimuth.

[0033] The dominant azimuth angle optimization method based on binary coupling in this invention is based on geological data such as well logging and well logging, and fully utilizes the anisotropy analysis of OVT gather data. The optimized dominant azimuth angle range is reliable and can broadly guide the selection of dominant azimuth angles when stacking OVT domain gather data by azimuth. Compared with the prior art, the dominant azimuth angle optimization method based on binary coupling in this invention can conveniently utilize the amplitude and azimuth variation characteristics of OVT gather data to optimize the dominant azimuth angle of the target layer, guide the azimuth stacking of pre-stack seismic data in the exploration area, and provide high-quality seismic data for comprehensive seismic geological research in the exploration area. Attached Figure Description

[0034] Figure 1 This is a flowchart of a specific embodiment of the preferred azimuth angle optimization method based on binary coupling of the present invention;

[0035] Figure 2 This is a schematic diagram of a typical single-well synthetic record in the Shengli-1 work area according to a specific embodiment of the present invention. The target layer is Tg, and the corresponding seismic reflection axis is shown below. Figure 2 ;

[0036] Figure 3 This is a seismic amplitude anisotropy analysis diagram of the target layer Tg of the OVT gather in the Shengli No. 1 work area according to a specific embodiment of the present invention;

[0037] Figure 4 This is a travel-time anisotropy analysis diagram of the target layer Tg of the OVT gather in Shengli No. 1 work area according to a specific embodiment of the present invention;

[0038] Figure 5 This is a graph showing the variation of seismic amplitude with azimuth angle in a typical well of the OVT collection in the Shengli No. 1 work area according to a specific embodiment of the present invention.

[0039] Figure 6 This is a graph showing the travel time of a typical well in the OVT collection area of ​​Shengli No. 1 working zone as a function of azimuth angle in a specific embodiment of the present invention.

[0040] Figure 7 This is a graph showing the variation of seismic amplitude with azimuth angle in a typical well of the OVT collection in the Shengli No. 1 work area according to a specific embodiment of the present invention.

[0041] Figure 8This is a graph showing the travel time of a typical well in the OVT collection area of ​​Shengli No. 1 working zone as a function of azimuth angle in a specific embodiment of the present invention.

[0042] Figure 9 This is a coupled graph of the seismic amplitude and travel time versus azimuth rate of change curves of a typical well in the OVT collection of the Shengli No. 1 work area, according to a specific embodiment of the present invention. Detailed Implementation

[0043] It should be noted that the following detailed descriptions are exemplary and intended to provide further illustration of the invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0044] It should be noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the exemplary embodiments of the present invention. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Furthermore, it should be understood that when the terms "comprising" and / or "including" are used in this specification, they indicate the presence of features, steps, operations, and / or combinations thereof.

[0045] The preferred azimuth angle optimization method based on binary coupling of the present invention specifically includes the following steps:

[0046] Step 1: Conduct comprehensive seismic geological analysis and single-well synthetic record calibration to determine the seismic reflection axis corresponding to the target layer. This includes:

[0047] Synthetic record calibration of typical wells was carried out based on comprehensive seismic geological analysis;

[0048] The comprehensive seismic geological analysis data includes well logging, well logging and drilling geological data, and seismic data.

[0049] Step 2: Conduct anisotropy analysis of the target layer of the OVT gather, specifically including:

[0050] By conducting azimuth anisotropy analysis of seismic amplitude and travel time of the target layer in the actual work area OVT gather, it is confirmed that the target layer has azimuth anisotropy, which allows for azimuth anisotropy research.

[0051] Step 3: Extract the travel time and seismic amplitude values ​​of the target layer from the OVT gather as a function of azimuth, specifically including:

[0052] After arranging the OVT gathers by azimuth, extract the seismic amplitude and travel time variation curves of the target layer at typical well locations as a function of azimuth.

[0053] Step 4: The formula for calculating the rate of change of seismic amplitude and travel time with azimuth for the target layer of the OVT gather is as follows:

[0054] ΔA = (AA) / A * 100

[0055] Where A is the amplitude value of the target layer extracted by the OVT gather, A is the seismic amplitude value at that point on the conventional seismic profile or the average amplitude of the target layer in each direction of the OVT gather, and ΔA is the seismic amplitude variation rate.

[0056] ΔT = (TT) / T * 100,

[0057] Where T is the travel time of the target layer extracted by the OVT gather, T is the travel time of that point on the conventional seismic profile or the average travel time of each direction of the target layer in the OVT gather, and ΔT is the rate of change of travel time.

[0058] Step 5: Calculate the coupling relationship value F between seismic amplitude and travel time rate of change of the target layer in the OVT gather, and find the azimuth range where both the travel time and seismic amplitude rate of change are small and stable, and the amplitude value is large.

[0059] F = ΔA / |ΔT|,

[0060] Where ΔA is the rate of change of earthquake amplitude and ΔT is the rate of change of travel time. The range of dominant F values ​​is defined according to the reflection characteristics of the target layer. The corresponding azimuth range is selected based on the dominant F value range, which is the dominant azimuth range.

[0061] The following are several specific embodiments of the application of the present invention.

[0062] Example 1

[0063] Applying a specific embodiment 1 of the present invention, the preferred azimuth angle optimization method based on binary coupling of the present invention specifically includes:

[0064] Step 1: Conduct comprehensive seismic geological analysis and perform single-well synthetic record calibration to determine the seismic reflection axis corresponding to the target layer. This specifically includes:

[0065] Synthetic record calibration of typical wells was carried out based on comprehensive seismic geological analysis;

[0066] The comprehensive seismic geological analysis data includes well logging, well logging and drilling geological data, and seismic data.

[0067] Step 2: Conduct anisotropy analysis of the target layer of the OVT gather, specifically including:

[0068] By conducting azimuth anisotropy analysis of seismic amplitude and travel time of the target layer in the actual work area OVT gather, it is confirmed that the target layer has azimuth anisotropy, which allows for azimuth anisotropy research.

[0069] Step 3: Extract the travel time and seismic amplitude values ​​of the target layer from the OVT gather as a function of azimuth. Specifically, this includes:

[0070] After arranging the OVT gathers by azimuth, extract the seismic amplitude and travel time variation curves of the target layer at typical well locations as a function of azimuth.

[0071] Step 4: Calculate the seismic amplitude and travel time variation rate of the target layer in the OVT gather.

[0072] When calculating the rate of change of seismic amplitude, the seismic amplitude calibrated by the single-well synthetic record or the average seismic amplitude of each direction of the target layer in the OVT gather can be selected as the reference value for calculation; when calculating the rate of change of travel time, the travel time calibrated by the single-well synthetic record or the average travel time of each direction of the target layer in the OVT gather can be selected as the reference value for calculation.

[0073] The formula for calculating the rate of change of seismic amplitude and travel time with azimuth for the target layer of the OVT gather is as follows:

[0074] ΔA = (A′ - A) / A * 100

[0075] Where A′ is the amplitude value of the target layer extracted by the OVT gather, A is the seismic amplitude value at that point on the conventional seismic profile or the average amplitude of the target layer in each direction of the OVT gather, and ΔA is the seismic amplitude variation rate;

[0076] ΔT=(T′-T) / T*100,

[0077] Where T′ is the travel time of the target layer extracted by the OVT gather, T is the travel time of that point on the conventional seismic profile or the average travel time of each direction of the target layer in the OVT gather, and ΔT is the rate of change of travel time.

[0078] Step 5: Calculate the coupling relationship value F between seismic amplitude and travel time rate of change of the target layer in the OVT gather, and find the azimuth range where both the travel time and seismic amplitude rate of change are small and stable, and the amplitude value is large.

[0079] The formula for calculating the coupling relationship F between seismic amplitude and travel-time rate of change of the target layer in the OVT gather is as follows:

[0080] F = ΔA / |ΔT|,

[0081] Where ΔA is the rate of change of earthquake amplitude and ΔT is the rate of change of travel time. The range of dominant F values ​​is defined according to the reflection characteristics of the target layer. The corresponding azimuth range is selected based on the dominant F value range, which is the dominant azimuth range.

[0082] Example 2

[0083] A specific embodiment 2 of the present invention is applied. Figure 1 This is a flowchart of a method for selecting the optimal azimuth angle based on binary coupling according to the present invention. The method for selecting the optimal azimuth angle based on binary coupling according to the present invention includes the following steps:

[0084] Step 101: Based on well logging, well logging, and drilling data, conduct comprehensive geological analysis, select typical wells for single-well seismic synthetic record calibration, and clarify the seismic reflection axis corresponding to the target layer. Figure 2 This is a typical single-well synthetic record from the Shengli No. 1 work area in a specific embodiment of the present invention. The process proceeds to step 102.

[0085] Step 102: Conduct anisotropic analysis of seismic amplitude and travel time of the target layer in the actual OVT gather of the work area to clarify that the target layer has azimuth anisotropy and can be studied. Figure 3 Anisotropy analysis diagram of amplitude of a specific layer in the OVT gather of Shengli No. 1 work area. Figure 4 The OVT gather for the Shengli No. 1 work area is shown in the travel time anisotropy analysis diagram for a target layer. The amplitude and travel time slices both indicate that the amplitude and travel time are different in different azimuths of the target layer, that is, there is azimuth anisotropy. The process proceeds to step 103.

[0086] Step 103: Extract the travel time and seismic amplitude variation curves of the target layer at typical well locations in the OVT gather, as well as the variation curves with azimuth. Figure 5 This is a curve showing the seismic amplitude versus azimuth at a typical well in the OVT collection of the Shengli No. 1 work area, according to a specific embodiment of the present invention. Figure 6 This is a travel time versus azimuth curve of a typical well at the OVT collection site in the Shengli-1 work area, as shown in a specific embodiment of the present invention. The process proceeds to step 104.

[0087] In step 104, using the average seismic amplitude and travel time of the target layer in the typical well composite record as A and T, the rate of change curves of seismic amplitude as a function of azimuth and travel time as a function of azimuth for the target layer at the typical well of the OVT gather are calculated. Figure 7 This is a specific embodiment of the present invention, showing the rate of change of seismic amplitude with azimuth at a typical well in the OVT collection of the Shengli No. 1 work area; Figure 8 This is a travel time versus azimuth rate curve of a typical well at the OVT collection point in the Shengli-1 work area, as shown in a specific embodiment of the present invention. The process proceeds to step 105.

[0088] In step 105, the coupling relationship between the amplitude and the rate of change of travel time of the target layer at a typical well in the OVT gather is calculated as F, where F = ΔA / |ΔT|. Figure 9 This is a coupling curve of seismic amplitude and travel time at a typical well in the Shengli-1 work area OVT gather, according to a specific embodiment of the present invention. In this area, an F value greater than 80 represents the dominant coupling range, corresponding to an azimuth of 55-115°. End of process.

[0089] This embodiment provides a method for optimizing the dominant azimuth angle based on binary coupling. By optimizing the dominant azimuth angle of the target layer based on the variation characteristics of the amplitude and travel time of the target layer in a typical single well of the OVT gather within the work area, the most suitable range of dominant azimuth angles for seismic interpretation or description of the target layer can be selected, thus establishing a foundation for azimuth-based superposition of OVT gathers.

[0090] Finally, it should be noted that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

[0091] Except for the technical features described in the specification, all other technologies are known to those skilled in the art.

Claims

1. A method for optimizing the preferred azimuth angle based on binary coupling, characterized in that, The method for optimizing the advantageous azimuth angle based on binary coupling includes: Step 1: Identify the seismic reflection axis corresponding to the target layer; Step 2: Conduct anisotropy analysis of the target layer of the OVT gather; Step 3: Extract the travel time and seismic amplitude values ​​of the target layer from the OVT gather as a function of azimuth; Step 4: Calculate the seismic amplitude and travel-time variation rate of the target layer in the OVT gather; Step 5: Couple the seismic amplitude and travel time variation rate of the target layer of the OVT gather to determine the range of dominant azimuth angles.

2. The method for selecting the optimal azimuth angle based on binary coupling according to claim 1, characterized in that, In step 1, synthetic record calibration of typical wells is carried out based on comprehensive seismic geological analysis to identify the seismic reflection axis corresponding to the target layer.

3. The method for selecting the optimal azimuth angle based on binary coupling according to claim 2, characterized in that, In step 1, the comprehensive seismic geological analysis data includes: well logging, well logging and drilling geological data, and seismic data.

4. The method for selecting the optimal azimuth angle based on binary coupling according to claim 1, characterized in that, In step 2, an azimuth anisotropy analysis of the seismic amplitude and travel time of the target layer in the actual work area OVT gather is carried out; it is confirmed that the target layer has azimuth anisotropy, and an azimuth anisotropy study is conducted.

5. The method for selecting the optimal azimuth angle based on binary coupling according to claim 1, characterized in that, In step 3, after arranging the OVT gathers by azimuth, the seismic amplitude and travel time curves of the target layer at typical well locations are extracted as a function of azimuth.

6. The method for selecting the optimal azimuth angle based on binary coupling according to claim 1, characterized in that, In step 4, when calculating the rate of change of seismic amplitude, the seismic amplitude calibrated by the single-well synthetic record or the average seismic amplitude of each direction of the target layer in the OVT gather can be selected as the reference value for calculation; when calculating the rate of change of travel time, the travel time calibrated by the single-well synthetic record or the average travel time of each direction of the target layer in the OVT gather can be selected as the reference value for calculation.

7. The method for selecting the optimal azimuth angle based on binary coupling according to claim 6, characterized in that, In step 4, the rate of change of seismic amplitude and travel time with azimuth for the target layer of the OVT gather is calculated using the following formula: ΔA = (A′ - A) / A * 100 Where A′ is the seismic amplitude value of the target layer extracted by the OVT gather, A is the seismic amplitude value at that point on the conventional seismic profile, and ΔA is the seismic amplitude variation rate; ΔT=(T′-T) / T*100, Where T′ is the travel time of the target layer extracted by the OVT gather, T is the travel time of that point on the conventional seismic profile, and ΔT is the rate of change of travel time.

8. The method for selecting the optimal azimuth angle based on binary coupling according to claim 1, characterized in that, In step 5, the coupling relationship value F between the seismic amplitude and travel time rate of change of the OVT gather is calculated. The range of azimuth angles with small and stable travel time and amplitude rate of change, and large amplitude value is identified as the dominant azimuth angle range.

9. The method for selecting the optimal azimuth angle based on binary coupling according to claim 8, characterized in that, In step 5, the formula for calculating the coupling relationship value F between the seismic amplitude and the rate of change of travel time of the target layer in the OVT gather is as follows: F = ΔA / |ΔT|, Where ΔA is the rate of change of earthquake amplitude and ΔT is the rate of change of travel time. The range of dominant F values ​​is defined according to the reflection characteristics of the target layer. The corresponding azimuth range is selected based on the dominant F value range, which is the dominant azimuth range.

10. A system for optimizing the preferred azimuth angle based on binary coupling, characterized in that, The dominant azimuth angle optimization system based on binary coupling adopts the dominant azimuth angle optimization method based on binary coupling as described in any one of claims 1-9, and uses the variation characteristics of seismic amplitude and travel time of the target layer with azimuth in OVT gather data to optimize the dominant azimuth angle of the target layer.

Citation Information

Patent Citations

  • Method and System for Converting OVT Domain Data into Azimuth Angle Domain Imaging Gathers

    CN111538084B

  • Method for determining azimuth gather division scheme

    CN116931083A

  • Method and system for accurately explaining minor fault based on azimuth angle optimization and attribute extraction of OVT domain gather

    CN117348073A