Method and device for determining location of interference source based on first arrival information, and storage medium

Abstract

The application discloses a positioning method and device for determining an interference source based on initial information and a storage medium, and steps of the method comprise: picking up the time of the interference source reaching each detector; performing grid processing on a work area, establishing a plane coordinate system for each grid, arranging the distance from the grid to each detector from small to large to form a horizontal coordinate, and arranging the initial time corresponding to the detector after arrangement to form a vertical coordinate; comparing each plane coordinate system, and if each point in the plane coordinate system corresponding to the grid is optimally matched with a rising trend, the grid is the position of the interference source. The device comprises an interference wave identification module, an initial time collection module, a grid division module, a offset distance calculation module and an image output module. The storage medium is used for storing a computer program for executing the method. The application can quickly determine the position of the interference source, greatly reduces the workload, improves the work efficiency, and is suitable for the field of seismic exploration.

Description

Technical Field

[0001] The present invention relates to a method for locating an interference source, in particular to a method for locating an interference source based on first arrival information. Background Art

[0002] With the continuous advancement of seismic exploration technology and the rapid development of the national economy, an increasing number of high-density 3D seismic exploration projects are inevitably being implemented in economically developed cities and surrounding areas. However, due to the presence of large equipment or sites generating seismic signals in these economically developed areas, 3D seismic acquisition often encounters high-energy, sudden interference. These external interference factors are becoming increasingly numerous, powerful, and widespread, significantly impacting the quality of seismic data, significantly reducing the signal-to-noise ratio (SNR) and increasing the difficulty of data processing. In these situations, the interference energy can affect hundreds or even thousands of raw data records, necessitating the suspension of acquisition operations and the rapid identification of the interference source. Targeted emergency measures, such as coordinated suspension and adjustment of acquisition time, can then be implemented. However, interference sources can arise from a wide range of sources, making a one-by-one investigation a long-range, slow, and inefficient process. Summary of the Invention

[0003] The purpose of the present invention is to provide a method for locating interference sources based on first arrival information, so as to achieve the purpose of quickly and accurately determining the location of the interference source;

[0004] A second object of the present invention is to provide a positioning device for determining interference sources based on first arrival information, and outputting a visual chart;

[0005] The third object of the present invention is to provide a computer-readable storage medium for storing a corresponding computer program for a method of locating an interference source based on first arrival information.

[0006] In order to achieve the above object, the technical solution adopted by the present invention is:

[0007] A method for locating an interference source based on first arrival information comprises the following steps:

[0008] S1: Pick up the first arrival information of the interference source

[0009] When interference source information appears on several geophones recording seismic exploration data in the work area, the time when the interference source arrives at each geophone is picked up, so that each geophone corresponds to a fixed first arrival time;

[0010] S2: Work Area Gridding

[0011] Grid the area within and around the work area. Set the grid size based on the search accuracy requirements and number the divided grids. Each grid is assumed to be a candidate point for the interference source location.

[0012] S3: Establishing a coordinate system

[0013] A plane coordinate system is established for each grid. The horizontal axis of the plane coordinate system is the distance between the grid and each detector, and the vertical axis is the time when the interference wave first reaches the detector.

[0014] When determining the coordinates of each point in the plane coordinate system corresponding to the grid, the distance between the grid and each detector is arranged in order from small to large to form the horizontal coordinate, and the first arrival time corresponding to the detector in the arranged order forms the vertical coordinate;

[0015] S4: Determine the location of the interference source

[0016] The plane coordinate systems established for each grid are compared. If each point in the plane coordinate system corresponding to the grid best matches the ascending trend, the grid is the location of the interference source.

[0017] As a limitation of the present invention, in S1, the method for picking up the interference source information in the geophone is: performing observation system loading and injection trace header information processing on the original record of seismic exploration data containing interference.

[0018] A positioning device for determining an interference source based on first arrival information, comprising:

[0019] Interference wave identification module: used to identify whether there are interference waves other than seismic waves;

[0020] First arrival time acquisition module: used to collect the time when the interference wave first arrives at each detector;

[0021] Grid division module: divide the area within and around the work area into grids;

[0022] Offset calculation module: used to calculate the distance between each grid and each detector;

[0023] Image output module: used to convert the coordinate points formed by the corresponding offset distance and first arrival time in each grid into a visual chart.

[0024] A computer-readable storage medium stores a computer program for executing the above-mentioned positioning method for determining an interference source based on first arrival information.

[0025] Due to the adoption of the above technical solution, the present invention has the following technical advancements compared with the prior art:

[0026] The present invention grids the work area and performs statistical analysis based on the first arrival information of the interference wave. By utilizing the propagation law that the information arrives at a closer position in a shorter time and arrives at a farther position in a longer time, as long as statistical analysis shows that the first arrival time tends to increase as the distance between the detector and the grid increases, the approximate location of the interference source can be quickly determined, and a large-scale investigation can be reduced to a small-scale investigation, which greatly narrows the scope of the investigation, greatly reduces the workload, and improves work efficiency.

[0027] In summary, the present invention can quickly determine the location of the interference source, greatly reduces the workload, improves work efficiency, and is suitable for the field of seismic exploration. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 This is a diagram of the signal picked up by the detector after the external interference signal first arrives in an embodiment of the present invention;

[0029] Figure 2 is a schematic diagram of a work area after grid division in an embodiment of the present invention;

[0030] Figure 3 is a plane coordinate diagram generated for grid 1 in an embodiment of the present invention;

[0031] Figure 4 is a plane coordinate diagram generated for grid 2 in an embodiment of the present invention;

[0032] Figure 5 It is a plane coordinate diagram generated for the grid 3 in an embodiment of the present invention. DETAILED DESCRIPTION

[0033] The present invention will be further described in detail below by way of specific examples. It should be understood that the described examples are only used to illustrate the present invention and are not intended to limit the present invention.

[0034] Example 1: Method for determining the location of interference sources based on first arrival information

[0035] In this embodiment, before determining the location of the interference source, several geophones have been deployed in the work area. The geophones are used to record seismic exploration data, that is, seismic wave information generated by the explosion of the blast point in the work area. When the information picked up by the geophones generates interference information from other vibration sources, the quality of the seismic recording will be affected. Figure 1 As shown in the figure, the upper part is the effective seismic wave signal picked up by the detector, and the lower part is the interference wave signal generated by the interference source. In order to remove the interference signal, it is necessary to determine the location of the interference source and avoid the influence of the interference source by coordinating the stop or adjusting the acquisition period. The specific steps to determine the location of the interference source are as follows:

[0036] S1: Pick up the first arrival information of the interference source

[0037] When an interference wave is generated within the work area, the several geophones within the work area can detect not only the seismic waves generated by the blast point explosion within the work area, but also the interference wave. Therefore, the geophones also record the interference wave information of the interference source. This interference information is the time when the interference source wave first reaches each geophone, i.e., the first arrival time. At this time, each geophone will correspond to a fixed first arrival time, i.e., the time when each geophone first records the reception of the interference wave. Since each geophone is at a different distance from the interference source, the first arrival time of the interference wave is also different, but it follows the natural law that the longer the distance, the longer the signal reception time. Assume that there are 6 geophones in this embodiment, and the numbers of the 6 geophones and the first arrival times of the interference waves recorded by each of them are A (8s), B (10s), C (12s), D (9s), E (11s), and F (13s).

[0038] Furthermore, the method for picking up the interference source information in the geophone is: performing observation system loading and injection trace header information processing on the original record of seismic exploration data containing interference.

[0039] S2: Work Area Gridding

[0040] like Figure 2 As shown, the area within and around the work area is divided into grids. The grid size is set according to the search accuracy requirement. If the search accuracy is required, the grid size can be smaller, and vice versa. The grids are numbered as Grid 1, Grid 2, Grid 3, Grid 4, etc., and each grid is assumed to be a candidate point for the interference source location.

[0041] S3: Establishing a coordinate system

[0042] A plane coordinate system is established for each grid. The horizontal axis of the plane coordinate system represents the distance between the grid and each detector, recorded as the offset, and the vertical axis represents the time when the interference wave first reaches the detector, recorded as the first arrival time. After the coordinate system is established, the position of each point in the coordinate system needs to be determined, that is, the horizontal and vertical coordinates of each point. A coordinate system is established for each grid, and a bar chart or line chart is generated. In this example, the coordinate system is established for grids 1 through 3 as an example to generate a bar chart.

[0043] like Figure 3 As shown in the figure, when establishing a coordinate system for grid 1, the coordinates of the center position of grid 1 and the coordinates of the six detectors are obtained, and the distances between grid 1 and the six detectors are calculated respectively. Then, the distances between grid 1 and the six detectors are arranged in order from small to large, namely A, D, B, E, C, and F. The distance values ​​form the horizontal axis, and the first arrival times corresponding to the six detectors are used as the vertical axis to establish a histogram.

[0044] like Figure 4As shown in the figure, when establishing a coordinate system for grid 2, the coordinates of the center position of grid 2 and the coordinates of the six detectors are obtained, and the distances between grid 2 and the six detectors are calculated respectively. Then, the distances between grid 2 and the six detectors are arranged in order from small to large, namely A, B, D, E, C, and F. The distance values ​​form the horizontal axis, and the first arrival times corresponding to the six detectors are used as the vertical axis to establish a histogram.

[0045] like Figure 5 As shown in the figure, when establishing a coordinate system for grid 3, the coordinates of the center position of grid 3 and the coordinates of the six detectors are obtained, and the distances between grid 3 and the six detectors are calculated respectively. Then, the distances between grid 3 and the six detectors are arranged in order from small to large, namely B, A, C, E, D, and F. The distance values ​​form the horizontal axis, and the first arrival times corresponding to the six detectors are used as the vertical axis to establish a histogram.

[0046] S4: Determine the location of the interference source

[0047] The plane coordinate system established for each grid is analyzed and compared. If each point in the plane coordinate system corresponding to the grid best matches the ascending trend, the grid is the location of the interference source.

[0048] Figure 3 In the plane coordinate system, the cylinders are in a rising state, which conforms to the rule that the first arrival time is short for close distances and long for far distances. Figure 3 In the plane coordinate system, the cylinders are in a fluctuating state, not rising in sequence. The cylinders that are closer have longer first arrival times, which does not conform to the rule that the first arrival time is short for close distances and long for far distances. Figure 4 The columns in the plane coordinate system fluctuate, not ascending in sequence. Closer distances have longer initial arrival times, which contradicts the pattern of shorter initial arrival times for closer distances and longer initial arrival times for farther distances. Therefore, after graphically analyzing each grid, statistical analysis reveals that grid 1 is the location of the interference source. Although grid 1 also covers a range, it significantly reduces the area for investigation.

[0049] Example 2 Positioning device for determining interference sources based on first arrival information

[0050] This embodiment is used to implement the positioning method of Example 1, and includes an interference wave identification module, a first arrival time acquisition module, a grid division module, an offset distance calculation module, and an image output module.

[0051] Interference wave identification module: When interference waves generated by interference sources other than seismic waves appear in the work area, the interference wave identification module judges and identifies each detector and sends the judgment and identification results to the first arrival time acquisition module and grid division module.

[0052] The First Arrival Time Acquisition Module collects the time when the interference wave first reaches each detector. Upon receiving the interference wave signal from the Interference Wave Identification Module, the First Arrival Time Acquisition Module collects the time when each detector first receives the interference wave and matches this time to each detector.

[0053] Grid Division Module: Divides the area within and around the work area into grids. When the grid division module receives the signal with interference waves sent by the interference wave identification module, the grid division module divides the area within and around the work area into grids based on the input search accuracy requirements.

[0054] The offset calculation module calculates the distance between each grid and each geophone. The grid division module transmits the divided signal to the offset calculation module. The offset calculation module identifies the coordinates of each grid center and each geophone, and calculates the distance between grid 1 and each geophone, the distance between grid 2 and each geophone, the distance between grid 3 and each geophone, and so on. After calculating the distance between each grid and each geophone, the distances between grid 1 and each geophone are sorted, the distances between grid 2 and each geophone are sorted, and the distances between grid 3 and each geophone are sorted.

[0055] The image output module converts the coordinates of the offset and first arrival time corresponding to each grid into a visual chart. The first arrival time acquisition module transmits the first arrival time information of each detector to the image output module. The offset calculation module transmits the offset information of each grid after sorting to the image output module. The image output module assigns the corresponding first arrival time vertical coordinate of each detector to the corresponding offset horizontal coordinate, forming the point value in each coordinate system. The final output is a bar chart or line chart.

[0056] Example 3 Computer-readable storage medium

[0057] The computer-readable storage medium in this embodiment stores a computer program, and when the computer program is executed by the processor, it is used to implement the positioning method for determining the interference source based on the first arrival information in Example 1.

[0058] The computer-readable storage medium stores non-transitory computer-readable instructions, which, when executed by a processor, execute all or part of the steps of the aforementioned methods.

[0059] The above-mentioned computer-readable storage media include, but are not limited to, optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or mobile hard disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROM (e.g., ROM cartridges).

[0060] The above description is merely an optional embodiment of the present disclosure and is not intended to limit the present disclosure. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present disclosure shall be included in the scope of protection of the present disclosure.

Claims

1. A method for locating interference sources based on first arrival information, characterized in that: The following steps are involved: S1: Pick up the first arrival information of the interference source When interference source information appears on several geophones recording seismic exploration data in the work area, the time when the interference source arrives at each geophone is picked up, so that each geophone corresponds to a fixed first arrival time; S2: Work Area Gridding Grid the area within and around the work area. Set the grid size based on the search accuracy requirements and number the divided grids. Each grid is assumed to be a candidate point for the interference source location. S3: Establishing a coordinate system A plane coordinate system is established for each grid. The horizontal axis of the plane coordinate system is the distance between the grid and each detector, and the vertical axis is the time when the interference wave first reaches the detector. When determining the coordinates of each point in the plane coordinate system corresponding to the grid, the distance between the grid and each detector is arranged in order from small to large to form the horizontal coordinate, and the first arrival time corresponding to the detector in the arranged order forms the vertical coordinate; S4: Determine the location of the interference source The plane coordinate systems established for each grid are compared. If each point in the plane coordinate system corresponding to the grid best matches the ascending trend, the grid is the location of the interference source.

2. The method for determining the interference source positioning based on first arrival information according to claim 1, characterized in that: In S1, the method for picking up the interference source information in the geophone is: performing observation system loading and injection trace header information processing on the original record of seismic exploration data containing interference.

3. A positioning device for determining interference sources based on first arrival information, characterized in that: The method for implementing claim 1 or 2 comprises: Interference wave identification module: used to identify whether there are interference waves other than seismic waves; First arrival time acquisition module: used to collect the time when the interference wave first arrives at each detector; Grid division module: divide the area within and around the work area into grids; Offset calculation module: used to calculate the distance between each grid and each detector; Image output module: used to convert the coordinate points formed by the corresponding offset distance and first arrival time in each grid into a visual chart.

4. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program for executing the positioning method for determining an interference source based on first arrival information according to claim 1 or 2.

Images