A method and system for acquiring spiral angle gathers of dual-azimuth streamer seismic data

Abstract

The present invention belongs to the technical field of seismic data processing, and specifically discloses a method and system for acquiring spiral angle gathers of two-azimuth streamer seismic data. Based on the two-azimuth streamer data, five-dimensional regularization processing is carried out to obtain a more ideal regularization effect of offset distance and azimuth angle distribution between different facets within the effective azimuth sector; by processing into an integrated study of crack prediction, a reasonable unequal sector division method is optimized to ensure the consistency of results of different azimuth sectors; based on the integrated well-seismic analysis, the maximum available offset distance range is quantitatively determined to ensure the reliability of azimuth anisotropic properties; amplitude-preserving signal-to-noise ratio improvement processing is carried out on the spiral incidence angle gather to ensure that the gather has a higher signal-to-noise ratio, meet the demand for pre-stack crack prediction with narrow azimuth data, and provide a data basis for carrying out pre-stack crack prediction.

Description

Technical Field

[0001] The present invention belongs to the technical field of seismic data processing, and in particular relates to a method and system for acquiring spiral angle gathers of dual-azimuth streamer seismic data. Background Art

[0002] Buried-hill reservoirs have become a key area of oil and gas exploration and development in domestic waters, and fracture prediction plays a crucial role in the implementation of these reservoirs. Currently, pre-stack fracture prediction in the industry primarily targets wide-azimuth seismic data. Based on anisotropy theory, and using theories of seismic wave propagation and attenuation in different azimuths, rock physics models of fractured media are established. Methods such as azimuthal anisotropic ellipse fitting and anisotropy inversion are used to describe faults and fracture reservoirs. Marine seismic exploration generally utilizes streamer acquisition. Due to the limitations of this acquisition method, the aspect ratio of streamer acquisition is typically less than 0.05, making it insufficient for pre-stack fracture prediction. In recent years, fracture prediction research for buried-hill reservoirs in offshore exploration areas has primarily focused on post-stack analysis.

[0003] With the deployment of secondary 3D acquisition at sea, many study areas now have two or more sets of streamer seismic data from different acquisition azimuths. These data provide the foundation for prestack fracture prediction. Streamer seismic data acquired at different times each have their own specific research objectives, and the corresponding streamer data acquisition parameters (including source parameters, gun cable depth, and observation system) also vary significantly. This results in inconsistencies in amplitude, frequency, phase, time difference, and bin attributes between blocks of seismic data from different acquisition azimuths. Processing streamer data from different acquisition azimuths requires first effectively eliminating these differences between blocks, then implementing a rational processing pipeline to ultimately obtain bi-directional fusion seismic data. Existing bi-directional fusion processing pipelines can unlock the potential of seismic data from different azimuths and improve imaging quality, but they do not fundamentally improve the study of azimuthal anisotropy and prestack fracture prediction. Summary of the Invention

[0004] The purpose of the present invention is to provide a method and system for acquiring spiral angle gathers of dual-azimuth streamer seismic data to solve the above-mentioned problems existing in the prior art.

[0005] In order to achieve the above object, the present invention adopts the following technical solutions:

[0006] In a first aspect, a method for acquiring spiral angle gathers of two-azimuth streamer seismic data is provided, comprising:

[0007] A. Acquire original streamer seismic data for azimuth 1 and azimuth 2, wherein the valid azimuth sector corresponding to azimuth 1 is α1-α2, where α1 and α2 are the azimuth lower limit and azimuth upper limit, respectively; and the valid azimuth sector corresponding to azimuth 2 is β1-β2, where β1 and β2 are the azimuth lower limit and azimuth upper limit, respectively;

[0008] B. Preprocess the original streamer seismic data of azimuth 1 and azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and perform seismic data regularization on the azimuth 1 seismic shot data A and azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg ;

[0009] C. Orientation-data volume A reg and position data volume B reg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b ;

[0010] D. performing five-dimensional data regularization processing on the azimuth one seismic shot gather data A and the azimuth two seismic shot gather data B to obtain a regularized data set C, wherein the effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2;

[0011] E. Perform azimuth sector division on the effective azimuth sectors α1-α2 and β1-β2, respectively, to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2;

[0012] F. Divide the regularized dataset C according to the small azimuth sectors a1, a2, and a3, and the small azimuth sectors b1, b2, and b3, to obtain data volume A1 corresponding to the small azimuth sector a1, data volume A2 corresponding to the small azimuth sector a2, data volume A3 corresponding to the small azimuth sector a3, data volume B1 corresponding to the small azimuth sector b1, data volume B2 corresponding to the small azimuth sector b2, and data volume B3 corresponding to the small azimuth sector b3;

[0013] G. Based on anisotropic velocity V and anisotropic velocity V a, anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers P b1 、P b2 and P b3 ;

[0014] H. Using the set data consistency judgment conditions to analyze the prestack depth migration CRP gather P a1 、P a2 and P a3 Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain the updated three small azimuth sectors b1, b2 and b3;

[0015] I. Based on the updated three small azimuth sectors a1, a2 and a3, and / or the updated three small azimuth sectors b1, b2 and b3, iteratively execute steps FH until the prestack depth migration CRP gather P is obtained. a1 、P a2 and P a3 The data consistency judgment conditions are met, and the prestack depth migration CRP gather P b1 、P b2 and P b3 The data consistency judgment conditions are met between them, and the final pre-stack depth migration CRP gather P is obtained. a1 、P a2 and P a3 , and the final prestack depth migration CRP gather P b1 、P b2 and P b3 ;

[0016] J. CRP gathers P from final prestack depth migration a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the pre-stack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 ;

[0017] K. The angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

[0018] In one possible design, after obtaining the spiral angle gather data volume G, the method further includes:

[0019] L. Perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G den , and output high signal-to-noise ratio spiral angle gather data volume G den .

[0020] In one possible design, preprocessing the original streamer seismic data at azimuth one and the original streamer seismic data at azimuth two respectively includes:

[0021] The original streamer seismic data at position one and the original streamer seismic data at position two are subjected to fidelity denoising, broadband processing and multiple wave suppression processing respectively, and the original streamer seismic data at position one and the original streamer seismic data at position two are subjected to consistency processing of amplitude, frequency, phase and time difference.

[0022] In a possible design, performing seismic data regularization processing on the azimuth one seismic shot gather data A and the azimuth two seismic shot gather data B respectively includes:

[0023] Seismic data regularization is performed on the seismic shot gather data A at azimuth one and the seismic shot gather data B at azimuth two using a data regularization method based on discrete Fourier transform.

[0024] In a possible design, the five-dimensional data regularization processing of the azimuth one seismic shot gather data A and the azimuth two seismic shot gather data B includes:

[0025] The five-dimensional data regularization processing method based on discrete Fourier transform is used to perform five-dimensional data regularization processing on the seismic shot gather data A of azimuth one and the seismic shot gather data B of azimuth two.

[0026] In one possible design, the data consistency judgment conditions include: the difference in average coverage times is not greater than 20%, the difference in average signal-to-noise ratio of the target layer is not greater than 20%, and the difference in average root mean square amplitude within the statistical time window is not greater than 10%.

[0027] In a possible design, the forward modeling gathers of each well are used to evaluate the AVO characteristics of the corresponding near-well seismic gathers to determine the corresponding resection parameters, including:

[0028] The AVO variation trend of the well forward modeling gathers is statistically analyzed, and the effective offset range is selected according to the AVO variation trend of the well forward modeling gathers. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the effective offset range to determine the available offset range of each target layer. The corresponding resection parameters are determined based on the available offset range of each target layer.

[0029] In a second aspect, a system for acquiring spiral angle gathers of two-directional streamer seismic data is provided, comprising a data acquisition unit, a data sorting unit, a velocity modeling unit, a rule processing unit, a sector division unit, a data division unit, an offset processing unit, a difference determination unit, a resection conversion unit, and a data merging unit, wherein:

[0030] a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 1 is α1-α2, wherein α1 and α2 are the azimuth lower limit and the azimuth upper limit of azimuth 1, respectively; and a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 2 is β1-β2, wherein β1 and β2 are the azimuth lower limit and the azimuth upper limit of azimuth 2, respectively;

[0031] The data sorting unit is used to pre-process the original streamer seismic data of azimuth 1 and the original streamer seismic data of azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and to perform seismic data regularization processing on the azimuth 1 seismic shot data A and the azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg ;

[0032] Velocity modeling unit, used for orientation data volume A reg and position data volume B reg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b ;

[0033] A rule processing unit is used to perform five-dimensional data regularization processing on the azimuth 1 seismic shot gather data A and the azimuth 2 seismic shot gather data B to obtain a regularized data set C, wherein the effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2;

[0034] a sector division unit, configured to divide the effective azimuth sectors α1-α2 and β1-β2 into azimuth sectors, respectively, to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2;

[0035] a data partitioning unit, configured to partition the regularized data set C according to the small azimuth sectors a1, a2, and a3 and the small azimuth sectors b1, b2, and b3, to obtain a data volume A1 corresponding to the small azimuth sector a1, a data volume A2 corresponding to the small azimuth sector a2, a data volume A3 corresponding to the small azimuth sector a3, a data volume B1 corresponding to the small azimuth sector b1, a data volume B2 corresponding to the small azimuth sector b2, and a data volume B3 corresponding to the small azimuth sector b3;

[0036] Migration processing unit for anisotropic velocity V and anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers P b1 、P b2 and P b3 ;

[0037] The difference judgment unit is used to use the set data consistency judgment conditions to judge the prestack depth migration CRP gather P a1 、P a2 and P a3 Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain the updated three small azimuth sectors b1, b2 and b3;

[0038] Excision conversion unit for CRP gathers P from prestack depth migration a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3Perform the resection processing and migrate the pre-stack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 ;

[0039] Data merging unit, used to merge the angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

[0040] In a possible design, the system further includes an amplitude preservation processing unit, which is used to perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G den , and output high signal-to-noise ratio spiral angle gather data volume G den .

[0041] In a third aspect, a system for acquiring spiral angle gathers of two-way streamer seismic data is provided, comprising:

[0042] a memory for storing instructions;

[0043] A processor is used to read the instructions stored in the memory and execute any one of the methods described in the first aspect according to the instructions.

[0044] Beneficial effects:

[0045] 1. The present invention carries out five-dimensional regularization processing on the two-azimuth shot gather data to obtain a more ideal regularization effect of the offset distance and azimuth angle distribution between different bins within the effective azimuth sector.

[0046] 2. The present invention ensures the consistency of results of sectors in different directions by optimizing a reasonable unequal sector division method.

[0047] 3. The present invention quantitatively determines the maximum available offset range based on integrated well-seismic analysis, ensuring the reliability of azimuthal anisotropic properties.

[0048] 4. The present invention performs amplitude-preserving signal-to-noise ratio processing on spiral incidence angle gathers to ensure that the angle gathers have a high signal-to-noise ratio and meet the needs of pre-stack fracture prediction with narrow azimuth data. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] 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 only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0050] Figure 1 Schematic diagram of the steps of the method in Example 1 of the present invention;

[0051] Figure 2 Schematic diagram of unequal sector division for two-way sectors;

[0052] Figure 3 This is a schematic diagram comparing the pre-stack depth migration processing results of unequally divided small azimuth sectors;

[0053] Figure 4 Schematic diagram of AVO gradient attributes corresponding to prestack depth migration CRP gathers of unequally divided small azimuth sectors;

[0054] Figure 5 The final spiral angle gather and the corresponding fracture development layer segment azimuth attribute analysis diagram are shown in Figure 2.

[0055] Figure 6 This is a schematic diagram of the system structure in Example 3 of the present invention. DETAILED DESCRIPTION

[0056] It should be noted that the description of these embodiments is intended to help understand the present invention, but does not constitute a limitation of the present invention. The specific structural and functional details disclosed herein are merely intended to describe exemplary embodiments of the present invention. However, the present invention may be embodied in a variety of alternative forms, and should not be construed as being limited to the embodiments set forth herein.

[0057] It should be understood that, unless otherwise expressly specified or limited, the term "connection" should be understood in a broad sense. For example, it can mean a fixed connection, an electrical connection, a direct connection, an indirect connection through an intermediate medium, or internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the embodiments based on specific circumstances.

[0058] In the following description, certain details are provided to facilitate a thorough understanding of the example embodiments. However, one skilled in the art will appreciate that the example embodiments may be practiced without these specific details. For example, devices may be shown in block diagrams to avoid obscuring the examples with unnecessary detail. In other embodiments, well-known processes, structures, and techniques may be shown without unnecessary detail to avoid obscuring the embodiments.

[0059] Example 1:

[0060] This embodiment provides a method for obtaining spiral angle gathers from dual-azimuth streamer seismic data. The purpose is to obtain azimuth-specific seismic results that can be used for pre-stack fracture prediction by applying a quasi-width azimuth processing flow to dual-azimuth streamer seismic data from offshore seabed. This method can be applied to corresponding seismic data processing servers, such as Figure 1 As shown, the method includes the following steps:

[0061] S1. Obtain the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two, the effective azimuth sector corresponding to azimuth one is α1-α2, wherein α1 and α2 are the lower azimuth limit and the upper azimuth limit of azimuth one, respectively; the effective azimuth sector corresponding to azimuth two is β1-β2, wherein β1 and β2 are the lower azimuth limit and the upper azimuth limit of azimuth two, respectively.

[0062] During specific implementation, it is necessary to first obtain two-azimuth streamer seismic data at sea, including the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two. The effective azimuth sector corresponding to azimuth one is α1-α2, and the effective azimuth sector corresponding to azimuth two is β1-β2. The azimuth angles of the effective azimuth sector of azimuth one and the effective azimuth sector of azimuth two can differ by 60 degrees. For example, the effective azimuth sector corresponding to azimuth one is 90-150 degrees, and the effective azimuth sector corresponding to azimuth two is 150-210 degrees.

[0063] S2. Preprocess the original streamer seismic data of azimuth 1 and azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and perform seismic data regularization on the azimuth 1 seismic shot data A and azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg .

[0064] In specific implementation, fine pre-processing of the two-azimuth streamer seismic data is an important foundation. The preprocessing includes performing fidelity denoising, broadband processing and multiple wave suppression processing on the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two, as well as consistency processing of amplitude, frequency, phase and time difference between the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two. Among them, the consistency processing of amplitude, frequency, phase and time difference can be completed by obtaining the amplitude matching factor, frequency matching factor, phase matching factor and time difference matching factor of the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two, and using the matching factors to perform four matching processing in series. Through pre-processing, two-azimuth shot gather data with good consistency can be obtained.

[0065] After preprocessing to obtain the azimuth one seismic shot data A and the azimuth two seismic shot data B, the data regularization method based on discrete Fourier transform can be used to perform seismic data regularization on the azimuth one seismic shot data A and the azimuth two seismic shot data B respectively to obtain the corresponding azimuth one data volume A. reg and position data volume B reg .

[0066] S3. Orientation-data volume A reg and position data volume B reg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b .

[0067] S4. Perform five-dimensional data regularization processing on the azimuth 1 seismic shot gather data A and the azimuth 2 seismic shot gather data B to obtain a regularized data set C. The effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2.

[0068] In specific implementation, a five-dimensional data regularization processing method based on discrete Fourier transform can be used to perform five-dimensional data regularization processing on the azimuth one seismic shot data A and the azimuth two seismic shot data B to obtain a regularized data set C. The effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2.

[0069] S5. Divide the effective azimuth sectors α1-α2 and β1-β2 into azimuth sectors respectively to obtain three small azimuth sectors a1, a2 and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2 and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2.

[0070] When implementing it specifically, Figure 2 As shown, the effective azimuth sectors α1-α2 and β1-β2 can be divided into azimuth sectors in an unequal manner to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, so that the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal and greater than the azimuth angle range of the small azimuth sector b2. For example, a1 is 90-115 degrees, a2 is 115-125 degrees, a3 is 125-150 degrees, b1 is 150-175 degrees, b2 is 175-185 degrees, and b3 is 185-210 degrees.

[0071] S6. Divide the regularized data set C according to the small orientation sectors a1, a2 and a3 and the small orientation sectors b1, b2 and b3 to obtain the data body A1 corresponding to the small orientation sector a1, the data body A2 corresponding to the small orientation sector a2, the data body A3 corresponding to the small orientation sector a3, the data body B1 corresponding to the small orientation sector b1, the data body B2 corresponding to the small orientation sector b2 and the data body B3 corresponding to the small orientation sector b3.

[0072] In specific implementation, the regularized data set C can be divided according to the small azimuth sectors a1, a2 and a3 and the small azimuth sectors b1, b2 and b3 to obtain data volumes A1, A2, A3, B1, B2 and B3 corresponding to the small azimuth sectors.

[0073] S7. Based on anisotropic velocity V and anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers Pb1 、P b2 and P b3 .

[0074] In specific implementation, the anisotropic velocity V and the anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, carry out anisotropic prestack depth migration processing to obtain the corresponding prestack depth migration CR P gather P b1 、P b2 and P b3 The prestack depth migration processing results corresponding to the six small azimuth sectors are as follows: Figure 3 As shown in the example, it can be seen that the offset results corresponding to the six small azimuth sectors obtained by this method are consistent.

[0075] S8. Use the set data consistency judgment conditions to perform prestack depth migration CRP gather P a1 、P a2 and P a3 Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain three updated small azimuth sectors b1, b2 and b3.

[0076] In specific implementation, the pre-stack depth migration CRP gather P can be compared a1 、P a2 and P a3 If the difference between the three is too large and does not meet the set data consistency judgment conditions, the azimuth sector division of the effective azimuth sector α1-α2 can be changed to obtain the updated three small azimuth sectors a1, a2 and a3. Similarly, the prestack depth migration CRP gather P can be compared.b1 、P b2 and P b3 If the differences between the three are too large and do not meet the set data consistency criteria, the azimuth sector division of the valid azimuth sector β1-β2 can be modified to obtain three updated small azimuth sectors b1, b2, and b3. For example, the data consistency criteria include: 1. The difference in average coverage times is no greater than 20%; 2. The difference in average signal-to-noise ratio of the target layer is no greater than 20%; 3. The difference in average root mean square amplitude within the statistical time window is no greater than 10%. The principle of modification is to reduce the scope of small azimuth sectors with excessive coverage times and correspondingly expand the scope of other small azimuth sectors.

[0077] S9. Based on the updated three small azimuth sectors a1, a2 and a3, and / or the updated three small azimuth sectors b1, b2 and b3, iteratively execute steps S6-S8 until the prestack depth migration CRP gather P is obtained. a1 、P a2 and P a3 The data consistency judgment conditions are met, and the prestack depth migration CRP gather P b1 、P b2 and P b3 The data consistency judgment conditions are met between them, and the final pre-stack depth migration CRP gather P is obtained. a1 、P a2 and P a3 , and the final prestack depth migration CRP gather P b1 、P b2 and P b3 .

[0078] S10. From the final prestack depth migration CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the pre-stack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and Pb3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 .

[0079] In practice, the CRP gathers P can be migrated from the final pre-stack depth a1 、P a2 、P a3 、P b1 、P b2 and P b3 Extract the wellside seismic gathers of the drilled wells and obtain the corresponding well forward gathers. The well forward gathers can be obtained through forward simulation based on the drilled well data. Then, the AVO (Amplitu de variation with offset) trend of the well forward gathers is counted, and the effective offset range with normal AVO trend change is selected. The large offset distance where AVO shows abnormal change is removed. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the effective offset range to determine the available offset range of each target layer. Based on the available offset range of each target layer, the corresponding resection parameters are determined. Figure 4 As shown in the figure, the AVO gradient attribute planes corresponding to the prestack depth migration CRP gathers of six small azimuth sectors obtained by this method are shown. It can be seen that the consistency between different azimuths is good. Finally, the corresponding resection parameters can be used to adjust the final prestack depth migration CRP gathers P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the pre-stack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 .

[0080] S11. The angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and Gb3 Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

[0081] S12. Perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G den , and output high signal-to-noise ratio spiral angle gather data volume G den .

[0082] In specific implementation, the spiral angle gather data volume G can be processed to preserve amplitude to remove random noise and ensure the high signal-to-noise ratio characteristics of the spiral angle gather data volume, and the corresponding high signal-to-noise ratio spiral angle gather data volume G can be obtained. den , and output high signal-to-noise ratio spiral angle gather data volume G den . High signal-to-noise ratio spiral angle gather data volume G den It can meet the demand for pre-stack fracture prediction based on narrow azimuth data and provide data for subsequent fracture prediction work, such as Figure 5 As shown in Figure 3, this is the spiral angle gather finally obtained by this method. In the fracture development layer, analysis of its azimuthal attributes shows obvious azimuthal anisotropy characteristics, and the anisotropy intensity indicates the fracture development situation.

[0083] Example 2:

[0084] This embodiment provides a system for acquiring spiral angle gathers of two-azimuth streamer seismic data, comprising a data acquisition unit, a data sorting unit, a velocity modeling unit, a rule processing unit, a sector division unit, a data division unit, an offset processing unit, a difference determination unit, a resection conversion unit, and a data merging unit, wherein:

[0085] a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 1 is α1-α2, wherein α1 and α2 are the azimuth lower limit and the azimuth upper limit of azimuth 1, respectively; and a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 2 is β1-β2, wherein β1 and β2 are the azimuth lower limit and the azimuth upper limit of azimuth 2, respectively;

[0086] The data sorting unit is used to pre-process the original streamer seismic data of azimuth 1 and the original streamer seismic data of azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and to perform seismic data regularization processing on the azimuth 1 seismic shot data A and the azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg ;

[0087] Velocity modeling unit, used for orientation data volume A reg and position data volume Breg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b ;

[0088] A rule processing unit is used to perform five-dimensional data regularization processing on the azimuth 1 seismic shot gather data A and the azimuth 2 seismic shot gather data B to obtain a regularized data set C, wherein the effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2;

[0089] a sector division unit, configured to divide the effective azimuth sectors α1-α2 and β1-β2 into azimuth sectors, respectively, to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2;

[0090] a data partitioning unit, configured to partition the regularized data set C according to the small azimuth sectors a1, a2, and a3 and the small azimuth sectors b1, b2, and b3, to obtain a data volume A1 corresponding to the small azimuth sector a1, a data volume A2 corresponding to the small azimuth sector a2, a data volume A3 corresponding to the small azimuth sector a3, a data volume B1 corresponding to the small azimuth sector b1, a data volume B2 corresponding to the small azimuth sector b2, and a data volume B3 corresponding to the small azimuth sector b3;

[0091] Migration processing unit for anisotropic velocity V and anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers P b1 、P b2 and P b3 ;

[0092] The difference judgment unit is used to use the set data consistency judgment conditions to judge the prestack depth migration CRP gather P a1 、P a2 and P a3Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain the updated three small azimuth sectors b1, b2 and b3;

[0093] Excision conversion unit for CRP gathers P from prestack depth migration a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the pre-stack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 ;

[0094] Data merging unit, used to merge the angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

[0095] Furthermore, the system further comprises an amplitude preservation processing unit, which is used to perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G. den , and output high signal-to-noise ratio spiral angle gather data volume G den .

[0096] Example 3:

[0097] This embodiment provides a system for acquiring spiral angle gathers of two-way streamer seismic data, such as Figure 6 As shown, including:

[0098] Data interface, used to establish data connection between the processor and the external data terminal;

[0099] a memory for storing instructions;

[0100] The processor is used to read the instructions stored in the memory and execute the method for acquiring spiral angle gathers of dual-azimuth streamer seismic data in embodiment 1 according to the instructions.

[0101] Optionally, the system further includes an internal bus, through which the processor, memory, and data interface can be interconnected. The internal bus can be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, or an EISA (Extended Industry Standard Architecture) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc.

[0102] The memory may include, but is not limited to, random access memory (RAM), read-only memory (ROM), flash memory, first-in-first-out memory (FIFO) and / or first-in-last-out memory (FILO), etc. The processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0103] Example 4:

[0104] This embodiment provides a computer-readable storage medium having instructions stored thereon. When the instructions are executed on a computer, the computer is caused to execute the method for acquiring spiral angle gathers from dual-azimuth streamer seismic data in Example 1. The computer-readable storage medium refers to a data storage medium, which may include, but is not limited to, a floppy disk, an optical disk, a hard disk, a flash memory, a USB flash drive, and / or a memory stick. The computer may be a general-purpose computer, a dedicated computer, a computer network, or other programmable device.

[0105] This embodiment further provides a computer program product, which, when executed on a computer, executes the method for acquiring spiral angle gathers of dual-azimuth streamer seismic data in embodiment 1. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device.

[0106] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention shall be included in the scope of protection of the present invention.

Claims

1. A method for acquiring spiral angle gathers of dual-azimuth streamer seismic data, characterized in that: include: A. Acquire original streamer seismic data for azimuth 1 and azimuth 2, wherein the valid azimuth sector corresponding to azimuth 1 is α1-α2, where α1 and α2 are the azimuth lower limit and azimuth upper limit, respectively; and the valid azimuth sector corresponding to azimuth 2 is β1-β2, where β1 and β2 are the azimuth lower limit and azimuth upper limit, respectively; B. Preprocess the original streamer seismic data of azimuth 1 and azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and perform seismic data regularization on the azimuth 1 seismic shot data A and azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg ; C. Orientation-data volume A reg and position data volume B reg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b ; D. performing five-dimensional data regularization processing on the azimuth one seismic shot gather data A and the azimuth two seismic shot gather data B to obtain a regularized data set C, wherein the effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2; E. Perform azimuth sector division on the effective azimuth sectors α1-α2 and β1-β2, respectively, to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2; F. Divide the regularized dataset C according to the small azimuth sectors a1, a2, and a3, and the small azimuth sectors b1, b2, and b3, to obtain data volume A1 corresponding to the small azimuth sector a1, data volume A2 corresponding to the small azimuth sector a2, data volume A3 corresponding to the small azimuth sector a3, data volume B1 corresponding to the small azimuth sector b1, data volume B2 corresponding to the small azimuth sector b2, and data volume B3 corresponding to the small azimuth sector b3; G. Based on anisotropic velocity V and anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers P b1 、P b2 and P b3 ; H. Using the set data consistency judgment conditions to analyze the prestack depth migration CRP gather P a1 、P a2 and P a3 Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain the updated three small azimuth sectors b1, b2 and b3; I. Based on the updated three small azimuth sectors a1, a2 and a3, and / or the updated three small azimuth sectors b1, b2 and b3, iteratively execute steps FH until the prestack depth migration CRP gather P is obtained. a1 、P a2 and P a3 The data consistency judgment conditions are met, and the prestack depth migration CRP gather P b1 、P b2 and P b3 The data consistency judgment conditions are met between them, and the final pre-stack depth migration CRP gather P is obtained. a1 、P a2 and P a3 , and the final prestack depth migration CRP gather P b1 、P b2 and P b3 ; J. CRP gathers P from final prestack depth migration a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the prestack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 ; K. The angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

2. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: After obtaining the spiral angle gather data volume G, the method further includes: L. Perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G den , and output high signal-to-noise ratio spiral angle gather data volume G den .

3. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: The preprocessing of the original streamer seismic data of azimuth one and the original streamer seismic data of azimuth two respectively includes: The original streamer seismic data at position one and the original streamer seismic data at position two are subjected to fidelity denoising, broadband processing and multiple wave suppression processing respectively, and the original streamer seismic data at position one and the original streamer seismic data at position two are subjected to consistency processing of amplitude, frequency, phase and time difference.

4. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: The seismic data regularization processing is performed on the seismic shot gather data A of azimuth one and the seismic shot gather data B of azimuth two, respectively, including: Seismic data regularization is performed on the seismic shot gather data A at azimuth one and the seismic shot gather data B at azimuth two using a data regularization method based on discrete Fourier transform.

5. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: The five-dimensional data regularization processing of the azimuth one seismic shot gather data A and the azimuth two seismic shot gather data B comprises: The five-dimensional data regularization processing method based on discrete Fourier transform is used to perform five-dimensional data regularization processing on the seismic shot gather data A of azimuth one and the seismic shot gather data B of azimuth two.

6. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: The data consistency determination conditions include: the difference in average coverage times is not greater than 20%, the difference in average signal-to-noise ratio of the target layer is not greater than 20%, and the difference in average root mean square amplitude within the statistical time window is not greater than 10%.

7. The method for obtaining spiral angle gathers of dual-azimuth streamer seismic data according to claim 1, characterized in that: The method of using the forward modeling gathers of each well to evaluate the AVO characteristics of the corresponding near-well seismic gathers and determine the corresponding resection parameters includes: The AVO variation trend of the well forward modeling gathers is statistically analyzed, and the effective offset range is selected according to the AVO variation trend of the well forward modeling gathers. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the effective offset range to determine the available offset range of each target layer. The corresponding resection parameters are determined based on the available offset range of each target layer.

8. A system for acquiring spiral angle gathers of dual-azimuth streamer seismic data, used to implement the method for acquiring spiral angle gathers of dual-azimuth streamer seismic data according to any one of claims 1 to 7, characterized in that: It includes a data acquisition unit, a data sorting unit, a velocity modeling unit, a rule processing unit, a sector division unit, a data division unit, an offset processing unit, a difference determination unit, a resection conversion unit and a data merging unit, wherein: a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 1 is α1-α2, wherein α1 and α2 are the azimuth lower limit and the azimuth upper limit of azimuth 1, respectively; and a data acquisition unit, configured to acquire original streamer seismic data of azimuth 1 and original streamer seismic data of azimuth 2, wherein the effective azimuth sector corresponding to azimuth 2 is β1-β2, wherein β1 and β2 are the azimuth lower limit and the azimuth upper limit of azimuth 2, respectively; The data sorting unit is used to pre-process the original streamer seismic data of azimuth 1 and the original streamer seismic data of azimuth 2 respectively to obtain the corresponding azimuth 1 seismic shot data A and azimuth 2 seismic shot data B, and to perform seismic data regularization processing on the azimuth 1 seismic shot data A and the azimuth 2 seismic shot data B respectively to obtain the corresponding azimuth 1 data volume A. reg and position data volume B reg ; Velocity modeling unit, used for orientation data volume A reg and position data volume B reg Carry out bi-directional anisotropic prestack depth migration velocity modeling to obtain bi-directional fusion anisotropic velocity V and azimuth-corresponding anisotropic velocity V a and the anisotropic velocity V corresponding to orientation 2 b ; A rule processing unit is used to perform five-dimensional data regularization processing on the azimuth 1 seismic shot gather data A and the azimuth 2 seismic shot gather data B to obtain a regularized data set C, wherein the effective azimuth sectors corresponding to the regularized data set C include α1-α2 and β1-β2; a sector division unit, configured to divide the effective azimuth sectors α1-α2 and β1-β2 into azimuth sectors, respectively, to obtain three small azimuth sectors a1, a2, and a3 corresponding to the effective azimuth sectors α1-α2, and three small azimuth sectors b1, b2, and b3 corresponding to the effective azimuth sectors β1-β2, wherein the azimuth angle ranges of the small azimuth sectors a1 and a3 are equal in size and greater than the azimuth angle range of the small azimuth sector a2, and the azimuth angle ranges of the small azimuth sectors b1 and b3 are equal in size and greater than the azimuth angle range of the small azimuth sector b2; a data partitioning unit, configured to partition the regularized data set C according to the small azimuth sectors a1, a2, and a3 and the small azimuth sectors b1, b2, and b3, to obtain a data volume A1 corresponding to the small azimuth sector a1, a data volume A2 corresponding to the small azimuth sector a2, a data volume A3 corresponding to the small azimuth sector a3, a data volume B1 corresponding to the small azimuth sector b1, a data volume B2 corresponding to the small azimuth sector b2, and a data volume B3 corresponding to the small azimuth sector b3; Migration processing unit for anisotropic velocity V and anisotropic velocity V a , anisotropic prestack depth migration processing is performed using data volumes A1, A2, and A3, and the corresponding prestack depth migration CRP gathers P are obtained. a1 、P a2 and P a3 , based on the anisotropic velocity V and the anisotropic velocity V b , respectively, using the data volumes B1, B2 and B3, anisotropic prestack depth migration processing is carried out to obtain the corresponding prestack depth migration CRP gathers P b1 、P b2 and P b3 ; The difference judgment unit is used to use the set data consistency judgment conditions to judge the prestack depth migration CRP gather P a1 、P a2 and P a3 Perform data consistency judgment and use the set data consistency judgment conditions to perform data consistency judgment on the pre-stack depth migration CRP gather P b1 、P b2 and P b3 Perform data consistency judgment. When the prestack depth migration CRP gather P a1 、P a2 and P a3 When the data consistency judgment condition is not met, the azimuth sector division of the effective azimuth sector α1-α2 is changed to obtain the updated three small azimuth sectors a1, a2 and a3. b1 、P b2 and P b3 When the data consistency judgment condition is not satisfied, the azimuth sector division of the effective azimuth sector β1-β2 is changed to obtain the updated three small azimuth sectors b1, b2 and b3; Excision conversion unit for CRP gathers P from prestack depth migration a1 、P a2 、P a3 、P b1 、P b2 and P b3 The wellside seismic gathers are extracted respectively, and the corresponding well forward gathers are obtained. The AVO characteristics of the corresponding wellside seismic gathers are evaluated using the forward gathers of each well, and the corresponding resection parameters are determined. The final prestack depth migration CRP gather P is then a1 、P a2 、P a3 、P b1 、P b2 and P b3 Perform the resection processing and migrate the prestack depth of the resection processed CRP gather P a1 、P a2 、P a3 、P b1 、P b2 and P b3 Convert the incident angle gathers respectively to obtain the corresponding angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 ; Data merging unit, used to merge the angle gather data volume G a1 , G a2 , G a3 , G b1 , G b2 and G b3 Merge, place the trace headers and sort according to the spiral angle gather method to obtain the spiral angle gather data volume G.

9. The system for acquiring spiral angle gathers of dual-azimuth streamer seismic data according to claim 8, characterized in that: The system also includes an amplitude preservation processing unit, which is used to perform amplitude preservation processing on the spiral angle gather data volume G to obtain the corresponding high signal-to-noise ratio spiral angle gather data volume G. den , and output high signal-to-noise ratio spiral angle gather data volume G den .

10. A dual-azimuth streamer seismic data spiral angle gather acquisition system, characterized in that: include: a memory for storing instructions; A processor is used to read the instructions stored in the memory and execute the method for acquiring spiral angle gathers of dual-azimuth streamer seismic data according to any one of claims 1 to 7 according to the instructions.

Images