A method and device for removing abnormal first arrivals based on azimuth strip

Abstract

The application discloses a kind of based on azimuth strip's abnormal first arrival elimination method and device, including to all shot points according to geographical location is automatically divided into area, obtains multiple groups of shot points;All shot points are automatically divided according to 360° azimuth angle and obtain 360 rectangular strip areas;To the shot point data corresponding to each group, automatically extract the data on the receiving point of each shot point same azimuth angle in the rectangular strip area range, obtain multiple group data sets;Draw multiple group first arrival time-offset distance intersection diagram through the multiple group set data obtained;Through the multiple group first arrival time-offset distance intersection diagram drawn, automatically identify abnormal first arrival and eliminate abnormal data.The application is automatically divided into area through shot point, automatically divided into azimuth strip and the like means, excludes the case that first arrival time does not increase with offset distance in different directions possibly caused by the case that lateral velocity changes greatly, to improve the quality and efficiency of seismic wave first arrival picking.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of pre-stack seismic data preprocessing, and more particularly relates to an abnormal first arrival elimination method and device based on azimuth strips. BACKGROUND

[0002] Static correction is one of the key steps in seismic data processing in complex surface areas such as mountains, loess tableland and desert, and is of great significance to improving imaging results. However, the application effect of the current refraction static correction and tomographic static correction methods mainly depends on the quality of first arrival picking. With the increasing data volume of seismic exploration data processing, the workload of first arrival picking is also greatly increased. Therefore, it is more important to improve the precision and speed of automatic first arrival picking. At present, more and more software companies have launched various first arrival automatic picking software or algorithms. However, due to the complex near-surface conditions, the first arrival is affected by multiple factors, resulting in that the quality of the first arrival automatic picking is often not very ideal. Therefore, a large amount of manual modification is required to meet the production requirements, which greatly increases the workload and prolongs the processing cycle.

[0003] The current conventional abnormal value elimination method for first arrival automatic picking includes griding the shot-receiver distance and first arrival data in all seismic data, removing the abnormal first arrival data according to the threshold value of the first arrival number set in the grid; or obtaining a singular value threshold range according to the average value and standard deviation of the first arrival difference sequence after sorting the multiple first arrival data and the multiple coefficient values of the singular value threshold range, and iteratively removing the abnormal first arrival; or extracting multiple shot-receiver pairs from the seismic data, applying a basic static correction amount to the first arrival of the shot-receiver pairs to obtain a first time of the shot-receiver pairs, and setting a threshold value of the first time to identify the abnormal first arrival. However, the above methods do not consider the case that the lateral velocity changes greatly in a complex near-surface area, such as the case that the lateral velocity changes greatly in a complex near-surface area such as a mountain or a loess tableland. The above methods cannot exclude the case that the lateral velocity changes greatly in different directions, which may cause the first arrival time not to increase with the increase of offset distance, and cannot improve the accuracy of identifying abnormal first arrival, thereby cannot improve the quality of first arrival automatic picking and cannot reduce the workload of manual interactive modification of first arrival. SUMMARY

[0004] Therefore, the embodiments of the present application provide an abnormal first arrival elimination method and device based on azimuth strips, which at least solve the technical problem in the prior art that the case that the lateral velocity changes greatly in different directions, which may cause the first arrival time not to increase with the increase of offset distance, cannot be excluded, thereby the quality and efficiency of first arrival picking of seismic waves cannot be improved.

[0005] In a first aspect, the embodiments of the present application provide an abnormal first arrival elimination method based on azimuth strips, which includes:

[0006] The automatic area division is performed on all the shooting points according to geographical positions to obtain multiple groups of shooting points.

[0007] The 360° azimuth automatic division is performed on all the shooting points to obtain 360 rectangular strip areas.

[0008] For the shot gather data corresponding to each group of shooting points, the data of the receiving points within the rectangular strip area at the same azimuth of each shooting point is automatically extracted to obtain multiple data sets.

[0009] The multiple first arrival time-offset distance intersection diagrams are drawn through the multiple data sets.

[0010] The abnormal first arrivals are automatically identified and the abnormal data are removed through the multiple first arrival time-offset distance intersection diagrams.

[0011] Optionally, the automatic area division is performed on all the shooting points according to geographical positions to obtain multiple groups of shooting points, which comprises the following steps.

[0012] The automatic area division is performed on all the shooting points according to the geographical position proximity principle to obtain multiple groups of shooting points.

[0013] Optionally, the 360° azimuth automatic division is performed on all the shooting points to obtain 360 rectangular strip areas, which comprises the following steps.

[0014] The position of each shooting point is taken as a short edge midpoint, 2 times of the receiving point interval is taken as a short edge, and the maximum offset distance is taken as a long edge, and the 360° azimuth automatic division is performed to obtain 360 rectangular strip areas.

[0015] Optionally, the rectangular strip area is divided every 1° azimuth, the number of the obtained rectangular strip areas is 360, and the angle of the azimuth is 1°.

[0016] Optionally, for the shot gather data corresponding to each group of shooting points, the data of the receiving points within the rectangular strip area at the same azimuth of each shooting point is automatically extracted to obtain multiple data sets, which comprises the following steps: for the shot gather data corresponding to each group of shooting points, the seismic data of the receiving points within the rectangular strip area at the same azimuth is extracted every 1° azimuth, and then each group of shot gather data in each geographical area is automatically divided into 360 data sets according to the azimuth.

[0017] In a second aspect, the embodiment of the present application further provides an abnormal first arrival removal device based on an azimuth strip, which comprises:

[0018] The area division module is used for automatically dividing areas according to geographical positions of all the shooting points to obtain multiple groups of shooting points.

[0019] An azimuth angle division module is configured to automatically divide and obtain 360 rectangular strip regions according to 360° azimuth angles for all the shooting points.

[0020] A set module is configured to automatically extract data of receiving points in the rectangular strip region at the same azimuth angle for each shooting point for each group of shot gather data, and obtain multiple data sets.

[0021] A drawing module is configured to draw multiple first arrival time-offset intersection diagrams by using the obtained multiple data sets.

[0022] A rejection module is configured to automatically identify abnormal first arrivals and reject abnormal data by using the drawn multiple first arrival time-offset intersection diagrams.

[0023] Optionally, the region division module comprises:

[0024] The region division module is configured to automatically divide all the shooting points according to a geographical position proximity principle, thereby obtaining multiple shooting points.

[0025] Optionally, the azimuth angle division module comprises:

[0026] The azimuth angle division module is configured to take the position of each shooting point as a short edge midpoint, take 2 times of the receiving point interval as a short edge, take the maximum offset as a long edge, and automatically divide and obtain 360 rectangular strip regions along the direction of 360° azimuth angles.

[0027] Optionally, the angle of the azimuth angle is 1°.

[0028] Optionally, the set module comprises:

[0029] The set module is configured to extract seismic data of receiving points in the rectangular strip at the same azimuth angle every 1° azimuth angle for each group of shot gather data corresponding to each shooting point, and automatically divide each group of shot gather data in each geographical region into 360 data sets according to azimuth angles.

[0030] The present application has the following beneficial effects:

[0031] The present application can finely divide the corresponding rectangular strip region at each azimuth angle by automatically dividing all the shooting points and automatically dividing the azimuth strip, strictly follows the principle that the first arrival time of a seismic wave in the same direction increases with the increase of offset, eliminates the influence of unreasonable first arrival time caused by the change of lateral velocity in different directions, and makes the region division in the far offset region more fine and the first arrival quality control at the far offset end more accurate, thereby improving the accuracy of the first arrival automatic picking quality control as a whole and providing a good data basis for the later static correction link.

[0032] Further, the whole process of removing abnormal first arrivals is automated, reducing the workload of manual interaction and greatly improving the efficiency of quality control of automatic first arrival picking, saving the project time and cost of seismic exploration data processing.

[0033] Other features and advantages of the present application will be made clear to those skilled in the art from the following detailed description of the application. BRIEF DESCRIPTION OF DRAWINGS

[0034] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings in which like reference characters refer to like parts throughout the figures, and in which:

[0035] Figure 1 A first arrival time-offset crossplot of the prior art is shown;

[0036] Figure 2 A schematic diagram of the prior art of dividing fan-shaped regions according to azimuth angle is shown;

[0037] Figure 3 An effect diagram of the prior art of abnormal first arrival removal of fan-shaped region division is shown;

[0038] Figure 4 A schematic diagram of the present application of dividing rectangular strip regions according to azimuth angle is shown;

[0039] Figure 5 A flowchart of the present application of abnormal first arrival removal based on azimuth strip is shown;

[0040] Figure 6 An effect diagram of the present application of abnormal first arrival removal based on azimuth strip is shown. DETAILED DESCRIPTION

[0041] Preferred embodiments of the present application will be described in more detail below. Although the preferred embodiments of the present application are described below, it is understood that the present application can be implemented in various forms and should not be limited by the embodiments set forth herein.

[0042] Currently, the basic quality control method of automatic first arrival picking is to form a time-offset crossplot using the first arrival time and offset information of all seismic data in the work area, such as Figure 1As shown in the figure, it can be seen that the scattered points far away from the main axis of the time-offset data are abnormal first arrivals. Directly removing these abnormal first arrivals from the seismic data is a relatively rough quality control method, because this method superimposes all the seismic data for comparison, without considering factors such as the geographical location of the work area, the offset size, and the lateral velocity change, and may also remove the correct first arrival data; for example, in a complex near-surface situation with large lateral velocity changes, the first arrival times of receiving points in the same offset segment in different azimuth directions may vary greatly, and even the situation where the larger the offset, the smaller the first arrival time may occur. In this case, it is inappropriate to put all the first arrivals at different azimuths together for abnormal removal.

[0043] There is also a relative Figure 1 A more accurate method is to use the coordinates of each excitation point in the work area as the center of the circle and the maximum offset as the radius, and automatically divide the circle into 4-8 fan-shaped areas according to the 45° or 90° azimuth angle, such as Figure 2 As shown in the figure, the coordinate origin represents the excitation point, and the scattered points distributed along the line segments represent the receiving points. The seismic data of the receiving points in the fan-shaped area corresponding to the same azimuth of each excitation point are extracted to draw the first arrival time-offset intersection diagram, and the abnormal first arrivals are eliminated based on these data sets. This follows the law of seismic wave propagation along rays to a certain extent. However, this method of dividing the fan-shaped area by azimuth is also relatively rough, especially in areas with far offset distances. The fan-shaped area covers a large area, and the quality control of far-offset seismic data is still insufficient.

[0044] An embodiment of the present invention provides a method for eliminating abnormal first arrivals based on azimuth strips, comprising:

[0045] All excitation points are automatically divided into regions according to their geographical locations to obtain multiple groups of excitation points;

[0046] Automatically dividing all the excitation points into 360 rectangular strip areas according to a 360° azimuth angle;

[0047] For the shot gather data corresponding to each group of excitation points, the data on the receiving points within the rectangular strip area at the same azimuth angle of each excitation point are automatically extracted to obtain multiple sets of data;

[0048] Draw multiple sets of first arrival time-offset intersection diagrams using multiple sets of acquired data sets;

[0049] Abnormal first arrivals are automatically identified and abnormal data are eliminated by drawing multiple sets of first arrival time-offset intersection diagrams.

[0050] The application can finely divide the corresponding rectangular strip region on each azimuth angle by automatic regional division and azimuth strip automatic division of all shooting points, strictly follows the principle that the first arrival time of seismic wave in the same direction increases with the increase of offset, eliminates the unreasonable influence of lateral velocity variation on the first arrival time in different directions, so that the regional division of the far offset area is more fine, the first arrival quality control of the far offset end is more accurate, thereby the precision of the first arrival automatic picking quality control is improved as a whole, and a good data basis is provided for the later static correction link.

[0051] Optionally, all the shooting points are automatically regionally divided according to geographical positions to obtain a plurality of groups of shooting points, including:

[0052] The all shooting points are automatically regionally divided according to geographical position proximity principle to obtain a plurality of groups of shooting points.

[0053] Optionally, the all shooting points are automatically divided according to 360° azimuth and 360 rectangular strip regions are obtained, including:

[0054] Each position of the shooting point is taken as a short edge midpoint, 2 times of the receiving point interval is taken as a short edge, and the maximum offset is taken as a long edge, so that 360 rectangular strip regions are automatically divided and obtained along the direction of 360° azimuth.

[0055] Optionally, the rectangular strip regions are divided every 1° azimuth, the number of the obtained rectangular strip regions is 360, and the angle of the azimuth is 1°.

[0056] Optionally, for the shot gather data corresponding to each group of shooting points, the data of the receiving points in the rectangular strip region range on the same azimuth of each shooting point are automatically extracted to obtain a plurality of data sets, including: for the shot gather data corresponding to each group of shooting points, the seismic data of the receiving points in the rectangular strip region range on the same azimuth are extracted every 1° azimuth, and then each group of shot gather data in each geographical region is automatically divided into 360 data sets according to the azimuth.

[0057] The embodiment of the application further provides an abnormal first arrival removing device based on an azimuth strip, including:

[0058] The regional division module is used for automatically regionally dividing all the shooting points according to geographical positions to obtain a plurality of groups of shooting points.

[0059] The azimuth division module is used for automatically dividing the all shooting points according to 360° azimuth and obtaining 360 rectangular strip regions.

[0060] The collection module automatically extracts data on the receiving points within the range of the rectangular strip region at the same azimuth angle of each excitation point from the shot gather data corresponding to each group of excitation points, and obtains multiple sets of data collections.

[0061] The drawing module draws multiple sets of first arrival time-offset intersection diagrams through the obtained multiple sets of data collections.

[0062] The elimination module automatically identifies abnormal first arrivals and eliminates abnormal data through the multiple sets of first arrival time-offset intersection diagrams drawn.

[0063] Optionally, the region division module comprises:

[0064] The region division module automatically divides the regions according to the principle of geographical position proximity, thereby obtaining multiple groups of excitation points.

[0065] Optionally, the azimuth angle division module comprises:

[0066] The azimuth angle division module takes the position of each excitation point as a short edge midpoint, takes 2 times the receiving point interval as a short edge, and takes the maximum offset as a long edge, and automatically divides and obtains 360 rectangular strip regions along the direction of 360° azimuth angle.

[0067] Optionally, the angle of the azimuth angle is 1°.

[0068] Optionally, the collection module comprises:

[0069] The collection module automatically divides each set of shot gather data in each geographical region into 360 data collections according to the azimuth angle after extracting the seismic data on the receiving points in the rectangular strip at the same azimuth angle every 1°.

[0070] Embodiment one:

[0071] As shown in Figure 5 , a method for removing abnormal first arrivals based on azimuth strips, the method comprising the following steps:

[0072] S1, automatic region division

[0073] Since the area covered by all the excitation points is large, the region division module first needs to automatically divide the regions according to the principle of geographical position proximity, thereby reducing the influence of lateral velocity variation in a large range on the identification of abnormal first arrivals.

[0074] S2, automatic azimuth strip division

[0075] The azimuth strip division module takes the position of each excitation point as a short edge midpoint, takes 2 times the receiving point interval as a short edge, and takes the maximum offset as a long edge, and automatically divides 360 rectangular strip regions along the direction of 360° azimuth angle, as shown in Figure 4as shown.

[0076] S3, automatically extracting azimuth strip data

[0077] For the shot gather data corresponding to each group of shot points in step S1, a set of data groups of the receiving points in the rectangular strip range at the same azimuth angle of each shot point is automatically extracted.

[0078] The first arrival time-offset intersection diagram is drawn by extracting the seismic data on the receiving points in the rectangular strip range corresponding to the same azimuth angle of each shot point, and the abnormal first arrival is removed on the basis of the data set, so that the law of seismic wave propagation along the ray is followed to a certain extent.

[0079] S4, drawing a time-offset intersection diagram

[0080] The first arrival time-offset intersection diagram is drawn according to the azimuth rectangular strip data set extracted in step S3.

[0081] S5, automatic removal of abnormal first arrival

[0082] The abnormal first arrival is automatically identified and removed according to the first arrival time-offset intersection diagram drawn in step S4, such as Figure 6 as shown.

[0083] The method can finely divide the corresponding rectangular strip area at each azimuth angle, strictly follows the principle that the first arrival time of seismic wave in the same direction increases with the increase of offset, eliminates the unreasonable influence of lateral velocity variation on the first arrival time in different directions, and especially the coverage area of the rectangular strip in the far offset area is much smaller than the commonly used fan-shaped area. Therefore, the regional division in the far offset area is more fine, and the first arrival quality control at the far offset end is more accurate.

[0084] In order to verify the effect of the present application, the automatic first arrival picking result of the actual field acquisition seismic data is selected to verify the quality control effect.

[0085] Figure 3 is the result after the abnormal first arrival is removed by the method of fan-shaped azimuth area division. As can be seen from the figure, in the area with large offset, part of the abnormal first arrival is not identified, so the applicability of the method in the far offset area needs to be improved.

[0086] Figure 6 is the result obtained after the abnormal first arrival is removed by the method of the present application based on the azimuth rectangular strip, and then the normal first arrival result of the adjacent receiving points is interpolated and compensated at the receiving point where the abnormal first arrival is removed. It can be seen that in the far offset area, Figure 3After all the abnormal first arrivals of the convex are identified and removed, the missing data are compensated by interpolation method, which greatly improves the effect of automatic picking of the first arrival and lays a good foundation for the later static correction.

[0087] Embodiment two:

[0088] The embodiment of the application provides an abnormal first arrival removing device based on azimuth strips, which comprises:

[0089] A region division module is configured to automatically divide all the shot points into regions according to geographical positions and obtain multiple groups of shot points.

[0090] Specifically, the region division module comprises:

[0091] The region division module automatically divides all the shot points into regions according to geographical positions, thereby obtaining multiple groups of shot points.

[0092] An azimuth angle division module is configured to automatically divide all the shot points into 360 rectangular strip regions according to 360° azimuth angles and obtain the 360 rectangular strip regions.

[0093] Specifically, the angle of the azimuth angle is 1°.

[0094] Specifically, the azimuth angle division module comprises:

[0095] The azimuth angle division module automatically divides all the shot points into 360 rectangular strip regions according to 360° azimuth angles and obtains the 360 rectangular strip regions, with the position of each shot point as the midpoint of the short side, 2 times the receiving point interval as the short side, and the maximum offset distance as the long side.

[0096] A set module is configured to automatically extract the data of the receiving points in the rectangular strip region range of the same azimuth angle of each shot point from the shot gather data corresponding to each group of shot points, and obtain multiple data sets.

[0097] Specifically, the set module comprises:

[0098] The set module automatically divides the shot gather data corresponding to each group of shot points into 360 data sets according to azimuth angles, after extracting the seismic data of the receiving points in the rectangular strip region of the same azimuth angle every 1° azimuth angle.

[0099] A drawing module is configured to draw multiple first arrival time-offset distance intersection diagrams from the obtained multiple data sets.

[0100] A removing module is configured to automatically identify abnormal first arrivals and remove abnormal data from the drawn multiple first arrival time-offset distance intersection diagrams.

[0101] The device can finely divide the corresponding rectangular strip area in each azimuth angle, strictly follows the principle that the first arrival time of seismic wave in the same direction increases with the increase of offset, eliminates the unreasonable influence of lateral velocity variation on the first arrival time in different directions, and especially the coverage area of the rectangular strip in the far offset area is far less than the commonly used fan-shaped area, so that the area division in the far offset area is more fine, and the first arrival quality control at the far offset end is more accurate.

[0102] In order to verify the technical effect of the present application, the automatic picking result of the first arrival of the actual field collected seismic data is selected to verify the quality control effect.

[0103] Figure 3 is the result after the abnormal first arrival is removed by the fan-shaped azimuth area division method, and from the figure, it can be seen that in the area with large offset, part of the abnormal first arrival is not identified, so the applicability of this method in the far offset area needs to be improved.

[0104] Figure 6 is the result after the abnormal first arrival is removed based on the azimuth rectangular strip in the present application, and after the normal first arrival result on the adjacent receiving point is interpolated and compensated at the receiving point where the abnormal first arrival is removed, it can be seen that in the far offset area, Figure 3 the protruding abnormal first arrival is identified and removed, and then the complete result is obtained by interpolation method, which greatly improves the effect of automatic picking of the first arrival, and lays a good foundation for the later static correction link.

[0105] In summary, through automatic area division of all shooting points, automatic azimuth strip division and other means, the present application can finely divide the corresponding rectangular strip area in each azimuth angle, strictly follows the principle that the first arrival time of seismic wave in the same direction increases with the increase of offset, eliminates the unreasonable influence of lateral velocity variation on the first arrival time in different directions, so that the area division in the far offset area is more fine, the first arrival quality control at the far offset end is more accurate, thereby the precision of the automatic picking quality control of the first arrival is improved, and a good data basis is provided for the later static correction link.

[0106] Further, the whole abnormal first arrival removal process is fully automated, reduces the workload of manual interaction and modification of the first arrival, greatly improves the efficiency of the automatic picking quality control of the first arrival, and saves the construction period and cost for the seismic exploration data processing project.

[0107] The above has described the embodiments of the present application, the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for eliminating abnormal first arrivals based on azimuth strips, characterized in that: include: All excitation points are automatically divided into regions according to their geographical locations to obtain multiple groups of excitation points; Automatically dividing all the excitation points into 360 rectangular strip areas according to a 360° azimuth angle; For the shot gather data corresponding to each group of excitation points, the data on the receiving points within the rectangular strip area at the same azimuth angle of each excitation point are automatically extracted to obtain multiple sets of data; Draw multiple sets of first arrival time-offset intersection diagrams using multiple sets of acquired data sets; Abnormal first arrivals are automatically identified and abnormal data are eliminated by drawing multiple sets of first arrival time-offset intersection diagrams.

2. The method for eliminating abnormal first arrivals based on azimuth strips according to claim 1 is characterized in that: All the excitation points are automatically divided into regions according to their geographical locations to obtain multiple groups of excitation points, including: All the excitation points are automatically divided into regions according to the principle of geographical proximity, thereby obtaining multiple groups of excitation points.

3. The method for eliminating abnormal first arrivals based on azimuth strips according to claim 1 is characterized in that: All the excitation points are automatically divided into 360 rectangular strip areas according to the 360° azimuth angle, including: With the position of each excitation point as the midpoint of the short side, twice the receiving point spacing as the short side, and the maximum offset as the long side, 360 rectangular strip areas are automatically divided and acquired along the 360° azimuth direction.

4. The method for eliminating abnormal first arrivals based on azimuth strips according to claim 3 is characterized in that: The rectangular strip area is divided every 1° azimuth angle, the number of the obtained rectangular strip areas is 360, and the division interval of the azimuth angle is 1°.

5. The method for eliminating abnormal first arrivals based on azimuth strips according to claim 1 is characterized in that: For the shot gather data corresponding to each group of excitation points, data on the receiving points within the rectangular strip area at the same azimuth of each excitation point are automatically extracted to obtain multiple data sets, including: for the shot gather data corresponding to each group of excitation points, seismic data on the receiving points within the rectangular strip at the same azimuth every 1° azimuth are extracted, and each group of shot gather data in each geographical area is automatically divided into 360 data sets according to azimuth.

6. A device for rejecting abnormal first arrivals based on azimuth strips, characterized in that: include: The regional division module is used to automatically divide all the excitation points into regions according to their geographical locations to obtain multiple groups of excitation points; An azimuth division module is used to automatically divide all the excitation points according to the 360° azimuth angle and obtain 360 rectangular strip areas; The collection module automatically extracts the data of the receiving points within the rectangular strip area at the same azimuth angle of each excitation point for the shot gather data corresponding to each group of excitation points, and obtains multiple groups of data sets; The drawing module draws multiple sets of first arrival time-offset intersection diagrams using multiple sets of acquired data sets; The elimination module automatically identifies abnormal first arrivals and eliminates abnormal data by drawing multiple sets of first arrival time-offset intersection diagrams.

7. The device for rejecting abnormal first arrivals based on azimuth strips according to claim 6, characterized in that: The area division module includes: All the excitation points are automatically divided into regions according to the principle of geographical proximity, thereby obtaining multiple groups of excitation points.

8. The device for rejecting abnormal first arrivals based on azimuth strips according to claim 6, characterized in that: The azimuth angle division module includes: With the position of each excitation point as the midpoint of the short side, twice the receiving point spacing as the short side, and the maximum offset as the long side, 360 rectangular strip areas are automatically divided and acquired along the 360° azimuth direction.

9. The device for rejecting abnormal first arrivals based on azimuth strips according to claim 8, characterized in that: The division interval of the azimuth angle is 1°.

10. The device for rejecting abnormal first arrivals based on orientation strips according to claim 6, characterized in that: The collection module includes: For the shot gather data corresponding to each group of excitation points, seismic data of receiving points within a rectangular strip at the same azimuth are extracted every 1° azimuth, and then each group of shot gather data in each geographical area is automatically divided into 360 data sets according to azimuth.

Images