Positioning method and device for determining interference source based on first arrival information and storage medium

Abstract

The invention discloses a positioning method and device for determining an interference source based on first arrival information, and a storage medium. The method comprises the following steps: picking up the time when the interference source reaches each detector; gridding processing is carried out on the work area, a plane coordinate system is established for each grid, the distances between the grids and the detectors are arranged from small to large to form abscissas, and the first arrival time corresponding to the arranged detectors forms ordinates; and comparing each plane coordinate system, and if each point in the plane coordinate system corresponding to the grid is optimally matched with the sequentially rising trend, determining that the grid is the position of the interference source. The device comprises an interference wave identification module, a first arrival time acquisition module, a grid division module, an offset calculation module and an image output module. The storage medium is used for storing a computer program for executing the method. The method can quickly determine the position of the interference source, greatly reduces the workload, improves the working 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, more and more high-density three-dimensional seismic exploration projects are inevitably entering economically developed towns and surrounding areas. However, in the process of three-dimensional seismic acquisition, due to the developed economy in the region, there are some large-scale equipment or places that generate vibration signals, and often encounter sudden interference with strong energy. These external interference factors are increasing in number, energy is getting stronger, and the range is getting wider, which has a great impact on the quality of the original data of seismic exploration, greatly reducing the signal-to-noise ratio of the acquired data and increasing the difficulty of data processing. In this case, the external interference energy affects hundreds of thousands of original records, and the acquisition construction can only be suspended, and the location of the interference source must be found as soon as possible, and then targeted emergency measures such as coordinated stop and adjustment of the acquisition period can be taken. However, the interference source may come from all directions, and the scope of one-by-one investigation is wide, the speed is slow, and the efficiency is low. Summary of the invention

[0003] The purpose of the present invention is to provide a method for locating an interference source based on first arrival information, so as to achieve the purpose of quickly and accurately determining the location of the interference source; A second object of the present invention is to provide a positioning device for determining interference sources based on first arrival information, which is used to output a visual chart; The third object of the present invention is to provide a computer-readable storage medium; used to store a corresponding computer program for a positioning method for determining an interference source based on first arrival information.

[0004] In order to achieve the above object, the technical solution adopted by the present invention is: A method for locating an interference source based on first arrival information comprises the following steps: S1: Pick up the first arrival information of the interference source When interference source information appears in 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 size of the grid according to the search accuracy requirements, number the divided grids, and assume that each grid is a candidate point for the location of the interference source; S3: Establishing the coordinate system A plane coordinate system is established for each grid, the horizontal axis of which 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 after the arrangement order forms the vertical coordinate; S4: Locate the source of interference 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.

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

[0006] A positioning device for determining an interference source based on first arrival information, comprising: 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 distance 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.

[0007] 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.

[0008] Due to the adoption of the above technical solution, the present invention has the following technical advances compared with the prior art: The present invention grids the work area, and performs statistical analysis in combination with the first arrival information of the interference wave, and utilizes 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, through statistical analysis, as the distance between the detector and the grid increases, the first arrival time tends to increase, the approximate location of the interference source can be quickly determined, and the large-scale investigation can be reduced to a small-scale investigation, which greatly reduces the scope of the investigation, greatly reduces the workload, and improves work efficiency.

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

[0010] Figure 1 It is a signal diagram picked up by the detector after the external interference signal first arrives in an embodiment of the present invention; Figure 2 is a schematic diagram of a work area after grid division processing in an embodiment of the present invention; Figure 3 is a plane coordinate diagram generated for the grid 1 in an embodiment of the present invention; Figure 4 is a plane coordinate diagram generated for grid 2 in an embodiment of the present invention; Figure 5 It is a plane coordinate diagram generated for the grid 3 in the embodiment of the present invention. DETAILED DESCRIPTION

[0011] The present invention is further described in detail below by specific examples. It should be understood that the described examples are only used to explain the present invention, and are not intended to limit the present invention.

[0012] Embodiment 1 Method for determining the location of interference sources based on first arrival information In this embodiment, before determining the location of the interference source, several geophones have been arranged in the work area. The geophones are used to record seismic exploration data, that is, seismic wave information generated by the explosion of the gun point in the work area. When the information picked up by the geophone generates interference information from other vibration sources, it will affect the quality of the seismic recording. 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: S1: Pick up the first arrival information of the interference source When interference waves are generated in the work area, several geophones in the work area can not only detect the seismic waves generated by the explosion of the gun point in the work area, but also detect the interference waves. Therefore, the geophones will also record the interference wave information of the interference source. The interference information is the time when the interference source wave first reaches each geophone, that is, the initial arrival time. At this time, each geophone will correspond to a fixed initial arrival time, that is, the time when each geophone first records the reception of the interference wave. Since the distance between each geophone and the interference source is different, the initial arrival time of the interference wave is also different, but it follows the natural law that the longer the distance, the longer the time to receive the signal. It is assumed that there are 6 geophones in this embodiment, and the numbers of the 6 geophones and the initial 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).

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

[0014] S2: Work area gridding like Figure 2 As shown, the area within the work area and the surrounding area is divided into grids, and the size of the grid is set according to the search accuracy requirement. If the search accuracy is required to be higher, the grid can be divided into smaller ones, otherwise, the grid is divided into larger ones. The divided grids are numbered to form grid 1, grid 2, grid 3, grid 4, etc., assuming that each grid is a candidate point of the interference source location.

[0015] S3: Establishing the 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, which is recorded as the offset distance, and the vertical axis is the time when the interference wave first arrives at the detector, which is recorded as the first arrival time. After the coordinate system is established, it is necessary to determine the position of each point in the coordinate system, that is, the horizontal and vertical coordinates of each point. A coordinate system is established for each grid, and a bar chart or a line chart is generated. In this embodiment, grids 1 to 3 are used as examples to establish a coordinate system and generate a bar chart.

[0016] like Figure 3 As shown, when establishing a coordinate system for grid 1, the coordinates of the center position of grid 1 are obtained, the coordinates of the six detectors are obtained, the distances between grid 1 and the six detectors are calculated respectively, and 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, and the distance values ​​form the horizontal axis. Then, the first arrival times corresponding to the six detectors are used as the vertical axis to establish a bar graph.

[0017] like Figure 4 As shown, when establishing a coordinate system for grid 2, the coordinates of the center position of grid 2 are obtained, the coordinates of the six detectors are obtained, the distances between grid 2 and the six detectors are calculated respectively, and 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, and the distance values ​​form the horizontal axis. Then, the first arrival times corresponding to the six detectors are used as the vertical axis to establish a bar graph.

[0018] like Figure 5 As shown, 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, and the distance values ​​form the horizontal axis. Then, the first arrival times corresponding to the six detectors are used as the vertical axis to establish a bar graph.

[0019] S4: Locate the source of interference 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.

[0020] Figure 3 In the plane coordinate system, the cylinders are in a state of rising successively, 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 close in distance have a longer initial arrival time, which does not conform to the rule that the initial arrival time is short for close distances and long for far distances. Figure 4 In the plane coordinate system, the cylinders fluctuate and do not rise in sequence. The closer the distance, the longer the initial arrival time is. This does not conform to the rule that the closer the initial arrival time is, the longer the initial arrival time is. Therefore, after analyzing each grid graphically, it can be seen through statistical analysis that the location of grid 1 is the location of the interference source. Although grid 1 also covers a range, it greatly reduces the area of ​​investigation.

[0021] Embodiment 2 Positioning device for determining interference source based on first arrival information This embodiment is used to implement the positioning method of Embodiment 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.

[0022] 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 the grid division module.

[0023] First arrival time acquisition module: used to collect the time when the interference wave first arrives at each detector. When the first arrival time acquisition module receives the signal with interference wave sent by the interference wave identification module, the first arrival time acquisition module collects the time when each detector first receives the interference wave, and matches the time with each detector.

[0024] 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 according to the input search accuracy requirements.

[0025] Offset distance calculation module: used to calculate the distance between each grid and each geophone. The grid division module transmits the divided signal to the offset distance calculation module, which identifies the coordinates of the center position of each grid and the coordinates of 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, etc. After the distance between each grid and each geophone is calculated, the distance between grid 1 and each geophone is sorted, the distance between grid 2 and each geophone is sorted, and the distance between grid 3 and each geophone is sorted.

[0026] Image output module: used to convert the coordinate points formed by 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, and the offset calculation module transmits the offset information of each grid after sorting to the image output module. The image output module assigns the first arrival time ordinate corresponding to each detector to the corresponding offset abscissa to form the point value in each coordinate system. Finally, the output forms a bar chart or a line chart.

[0027] Embodiment 3 Computer readable storage medium 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.

[0028] The computer-readable storage medium stores non-transitory computer-readable instructions. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of various embodiments are executed.

[0029] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0030] The above description is only an optional embodiment of the present disclosure and is not intended to limit the present disclosure. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present disclosure shall be included in the protection scope 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 in 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 size of the grid according to the search accuracy requirements, number the divided grids, and assume that each grid is a candidate point for the location of the interference source; S3: Establishing the coordinate system A plane coordinate system is established for each grid, the horizontal axis of which 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 after the arrangement order forms the vertical coordinate; S4: Locate the source of interference 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 locating interference sources 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 detector is: loading the original record of the seismic exploration data containing the interference with the observation system and processing the injection trace information.

3. A positioning device for determining interference sources based on first arrival information, characterized in that: include: 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 distance 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 as described in claim 1 or 2.

Images