Shot screening method and device, computer device, storage medium and product

Abstract

The application provides a shot point screening method and device, computer equipment, storage medium and product, and belongs to the technical field of earthquakes. The method comprises the following steps: acquiring seismic data; determining a first time range and a second time range corresponding to the seismic data respectively, the first time range being a time range corresponding to zero value data in the seismic data, and the second time range being a time range of missing seismic data; acquiring firing times of a plurality of first shot points; determining a target firing time range corresponding to the seismic data based on the first time range and the second time range; and determining an abnormal shot point in the plurality of first shot points based on the firing times of the plurality of first shot points and the target firing time range. The method reverses the target firing time range corresponding to the first time range and the second time range, and then determines the abnormal shot point in the plurality of shot points by matching the firing times of the shot points and the target firing time range, thereby improving the accuracy of shot point screening.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of seismic technology, and in particular to a shot point screening method and device, computer equipment, storage medium and product. BACKGROUND

[0002] Wireless node seismic acquisition technology has the advantages of low cost, convenient construction, high data acquisition efficiency, etc. In the process of seismic data acquisition, it plays an important role. In the application process of this technology, wireless nodes are used to collect seismic data, which is generated after the shot point excites seismic waves, so it is necessary to screen out the shot point that excites abnormal seismic waves in time to ensure the quality of the seismic data collected by the wireless node. SUMMARY

[0003] The embodiments of the present application provide a shot point screening method, device, computer equipment, storage medium and product, which can improve the accuracy of shot point screening. The technical solution is as follows:

[0004] In one aspect, a shot point screening method is provided, the method comprising:

[0005] obtaining seismic data, the seismic data being generated after a plurality of first shot points excite seismic waves;

[0006] determining a first time range and a second time range corresponding to the seismic data respectively, the first time range being a time range corresponding to zero value data in the seismic data, and the second time range being a time range of missing seismic data;

[0007] obtaining the excitation time of the plurality of first shot points;

[0008] determining a target excitation time range corresponding to the seismic data based on the first time range and the second time range;

[0009] determining an abnormal shot point in the plurality of first shot points based on the excitation time of the plurality of first shot points and the target excitation time range.

[0010] In some embodiments, the seismic data includes a plurality of data channels connected in sequence, and the acquisition time length of each of the plurality of data channels is a target time length. The determination process of the second time range comprises:

[0011] determining the difference between the start times corresponding to any two adjacent data channels in the plurality of data channels to obtain a plurality of difference values;

[0012] If the target difference value exceeding the target time length is included in the plurality of difference values, a sum of a start time of a first track of data and the target time length is taken as a start point of the second time range, and a start time of a second track of data is taken as an end point of the second time range, the first track of data and the second track of data being a previous track of data and a next track of data in the two tracks of data corresponding to the target difference value respectively.

[0013] In some embodiments, the seismic data includes seismic data collected by a plurality of wireless nodes, the target firing time range includes target firing time ranges corresponding to a plurality of abnormal wireless nodes respectively, the abnormal wireless node is a node in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless node has zero value data and loss situation, and the determination of the abnormal shot point in the plurality of first shot points based on the firing time of the plurality of first shot points and the target firing time range includes:

[0014] determining a plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of first shot points;

[0015] determining target shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively based on the target firing time range corresponding to the plurality of abnormal wireless nodes respectively and the firing time of the plurality of second shot points, the firing time of the target shot point matching the target firing time range;

[0016] determining the abnormal shot point from the target shot points corresponding to the plurality of abnormal wireless nodes respectively.

[0017] In some embodiments, the determination of the abnormal shot point from the target shot points corresponding to the plurality of abnormal wireless nodes respectively includes:

[0018] determining an abnormal track proportion of any target shot point by a ratio between a number of abnormal wireless nodes corresponding to the target shot point and a total number of wireless nodes corresponding to the target shot point;

[0019] determining the target shot point as an abnormal shot point in a case where the abnormal track proportion is greater than a target proportion.

[0020] In some embodiments, the seismic data comprises seismic data collected by a plurality of wireless nodes, the first time range and the second time range respectively comprise a first time range and a second time range corresponding to a plurality of abnormal wireless nodes, the target shooting time range comprises a target shooting time range corresponding to the plurality of abnormal wireless nodes, the abnormal wireless nodes are nodes in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless nodes has zero value data and loss; and the method comprises:

[0021] determining a target time range corresponding to the plurality of abnormal wireless nodes from the first time range and the second time range corresponding to the plurality of abnormal wireless nodes, the target time range having a length exceeding a length threshold;

[0022] determining a target shooting time range corresponding to the plurality of abnormal wireless nodes based on the target time range corresponding to the plurality of abnormal wireless nodes.

[0023] In some embodiments, the method comprises:

[0024] obtaining a length adjustment parameter, the length adjustment parameter being used to correct the target time range;

[0025] determining a difference between a start point of the target time range corresponding to the plurality of abnormal wireless nodes and the length adjustment parameter;

[0026] determining a time range composed of the difference and an end point of the target time range corresponding to the plurality of abnormal wireless nodes as the target shooting time range corresponding to the plurality of abnormal wireless nodes.

[0027] In another aspect, a device for shotpoint screening is provided, the device comprising:

[0028] a first obtaining module configured to obtain seismic data, the seismic data being generated after seismic waves are excited by a plurality of first shotpoints;

[0029] a first determining module configured to determine a first time range and a second time range corresponding to the seismic data, the first time range being a time range corresponding to zero value data in the seismic data, and the second time range being a time range of lost seismic data;

[0030] a second obtaining module configured to obtain shooting times of the plurality of first shotpoints;

[0031] a second determining module, configured to determine a target shooting time range corresponding to the seismic data based on the first time range and the second time range;

[0032] a third determining module, configured to determine an abnormal shot point in the multiple first shot points based on shooting times of the multiple first shot points and the target shooting time range.

[0033] In some embodiments, the seismic data comprises multiple channel data connected in sequence, and a collection time length of each of the multiple channel data is a target time length, and the device further comprises:

[0034] a fourth determining module, configured to determine a difference between start times corresponding to any two adjacent channel data in the multiple channel data, to obtain multiple difference values;

[0035] a fifth determining module, configured to, if the multiple difference values include a target difference value exceeding the target time length, take a sum of a start time of a first channel data and the target time length as a starting point of the second time range, take a start time of a second channel data as an ending point of the second time range, the first channel data and the second channel data being a preceding channel data and a subsequent channel data in the two channel data corresponding to the target difference value.

[0036] In some embodiments, the seismic data comprises seismic data collected by multiple wireless nodes, the target shooting time range comprises target shooting time ranges corresponding to multiple abnormal wireless nodes respectively, the abnormal wireless node is a node in the multiple wireless nodes, and seismic data collected by the abnormal wireless node has zero value data and loss, and the third determining module is configured to:

[0037] determine, from the multiple first shot points, multiple second shot points corresponding to the multiple abnormal wireless nodes respectively;

[0038] determine, from the multiple second shot points corresponding to the multiple abnormal wireless nodes respectively, target shot points corresponding to the multiple abnormal wireless nodes respectively based on the target shooting time ranges corresponding to the multiple abnormal wireless nodes respectively and shooting times of the multiple second shot points, the shooting time of the target shot point matching the target shooting time range;

[0039] determine the abnormal shot point from the target shot points corresponding to the multiple abnormal wireless nodes respectively.

[0040] In some embodiments, the third determining module is configured to:

[0041] determine an abnormal channel proportion of any target shot point based on a ratio between a number of abnormal wireless nodes corresponding to the target shot point and a total number of wireless nodes corresponding to the target shot point;

[0042] In a case where the proportion of the abnormal traces is greater than the target proportion, the target shot point is determined as an abnormal shot point.

[0043] In some embodiments, the seismic data comprises seismic data collected by a plurality of wireless nodes, the first time range and the second time range respectively comprise a first time range and a second time range corresponding to a plurality of abnormal wireless nodes respectively, the target firing time range comprises a target firing time range corresponding to the plurality of abnormal wireless nodes respectively, the abnormal wireless nodes are nodes in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless nodes has zero value data and loss, the second determining module is configured to:

[0044] determine, from the first time range and the second time range corresponding to the plurality of abnormal wireless nodes respectively, a target time range corresponding to the plurality of abnormal wireless nodes respectively, the target time range having a time length exceeding a time length threshold;

[0045] determine, based on the target time range corresponding to the plurality of abnormal wireless nodes respectively, a target firing time range corresponding to the plurality of abnormal wireless nodes respectively.

[0046] In some embodiments, the second determining module is configured to:

[0047] obtain a time length adjustment parameter, the time length adjustment parameter being used to correct the target time range;

[0048] determine a difference between a start point of the target time range corresponding to the plurality of abnormal wireless nodes respectively and the time length adjustment parameter;

[0049] compose, as the target firing time range corresponding to the plurality of abnormal wireless nodes respectively, a time range composed of the difference corresponding to the plurality of abnormal wireless nodes respectively and an end point of the target time range.

[0050] In another aspect, a computer device is provided, the computer device comprising one or more processors and one or more memories, the one or more memories having stored therein at least one program code, the at least one program code being loaded and executed by the one or more processors to implement the shot point screening method according to any of the above implementation manners.

[0051] In another aspect, a computer readable storage medium is provided, the computer readable storage medium having stored therein at least one program code, the at least one program code being loaded and executed by a processor to implement the shot point screening method according to any of the above implementation manners.

[0052] In another aspect, a computer program product is provided, including a computer program which, when executed by a processor, implements the shot screening method of any of the above implementation manners.

[0053] The embodiment of the present application provides a shot screening method. Since the first time range is a time range corresponding to zero value data in the seismic data, and the second time range is a time range of missing seismic data, the corresponding target shot time range is reversely obtained based on the first time range and the second time range, the abnormal shot in the multiple shots is determined by matching the shot time of the shot and the target shot time range, and the accuracy of the shot screening is improved. BRIEF DESCRIPTION OF DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0055] Figure 1 is a flowchart of a shot screening method provided by the embodiment of the present application;

[0056] Figure 2 is a flowchart of a method for determining an effective time range provided by the embodiment of the present application;

[0057] Figure 3 is a statistical information diagram of a wireless node provided by the embodiment of the present application;

[0058] Figure 4 is a statistical result diagram of integrity checking of an abnormal wireless node provided by the embodiment of the present application;

[0059] Figure 5 is an information schematic diagram of an abnormal shot provided by the embodiment of the present application;

[0060] Figure 6 is a flowchart of a method for determining a target shot provided by the embodiment of the present application;

[0061] Figure 7 is a flowchart of another shot screening method provided by the embodiment of the present application;

[0062] Figure 8 is a block diagram of a shot screening device provided by the embodiment of the present application;

[0063] Figure 9 is a block diagram of a terminal provided by the embodiment of the present application. DETAILED DESCRIPTION

[0064] For the purposes of the present application, the technical solutions and advantages thereof, the embodiments of the present application will be described in further detail below with reference to the drawings.

[0065] The terms "first", "second", "third", and "fourth" and the like in the description and in the claims of the present application and the drawings are used for distinguishing between similar objects, not necessarily described in a particular order. Also, the terms "comprise", "comprising", "including", and "having" and variations thereof in the description and in the claims of the present application are intended to cover both the case where the stated object is included in the process, method, system, product, or apparatus and the case where the stated object is not included in the process, method, system, product, or apparatus.

[0066] It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data for analysis, stored data, displayed data, etc.), and signals involved in the present application are all authorized by the user or fully authorized by all parties, and the collection, use, and processing of relevant data need to comply with relevant laws, regulations, and standards of relevant countries and regions. For example, the seismic data involved in the present application is obtained under full authorization.

[0067] Figure 1 is a flowchart of a shot point screening method provided by an embodiment of the present application, and the method comprises the following steps.

[0068] 101、The computer device acquires seismic data, which is generated after a plurality of first shot points excite seismic waves.

[0069] The seismic data comprises seismic data collected by a plurality of wireless nodes; one wireless node can collect seismic data generated after a plurality of first shot points excite seismic waves, and seismic data generated after one first shot point excites seismic waves can be collected by a plurality of wireless nodes; each wireless node is marked based on a receiving line number rl i and a point number rp i , and i represents the i th wireless node; each first shot point is marked based on a line number sl h , a point number sp h , and a shot index number si h , and h represents the h th first shot point.

[0070] The computer device obtains the seismic data collected by the wireless node by importing a data file corresponding to the wireless node. Optionally, the data files of the plurality of wireless nodes are represented as f1, f2, …, fi, …, fn, i represents the serial number of the data file, and n represents the number of data files. The seismic data collected by any wireless node includes a plurality of channels of data connected in sequence, and the collection time length of the plurality of channels of data is a target time length, which can be set and changed as needed, and is not specifically limited here.

[0071] 102. The computer device determines a first time range and a second time range corresponding to the seismic data, respectively.

[0072] The first time range is a time range corresponding to zero-value data in the seismic data, and the second time range is a time range of missing seismic data. The first time range and the second time range each include a first time range and a second time range corresponding to a plurality of abnormal wireless nodes, respectively. The abnormal wireless node is a node in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless node has zero-value data and a missing condition; optionally, at least one of the first time range and the second time range corresponds to any abnormal wireless node, and embodiments of the present application are described by taking an example in which the abnormal wireless node corresponds to the first time range and the second time range. Optionally, the first time range and the second time range corresponding to the abnormal wireless node are both multiple, that is, the abnormal wireless node corresponds to multiple time ranges of zero-value data and multiple time ranges of missing seismic data.

[0073] For the ith wireless node, the seismic data of the wireless node includes a plurality of channels of data connected in sequence, and the collection time length of the plurality of channels of data is a target time length. The computer device obtains the first time range, the receiving line number, and the receiving point number of each channel of data from the data channel header of each channel of data, respectively; the first time range corresponding to the abnormal wireless node is represented as ts i1 to te i1 , ts i2 to te i2 , …, ts iu to te iu , …, ts ip to te ip ; optionally, the zero-value data is at least one, u represents the serial number of a continuous zero-value data segment in the continuous channel of data of the ith wireless node, p represents the total number of time ranges of continuous zero-value data segments in the continuous channel of data of the ith wireless node, ts iu represents the start time of the u-th continuous zero-value data segment in the continuous channel of data of the ith wireless node, and te iu represents the end time of the u-th continuous zero-value data segment in the continuous channel of data of the ith wireless node.

[0074] In some embodiments, the determining process of the second time range comprises the following steps: the computer device determines the difference between the start times of any two adjacent data channels in the multi-channel data, obtaining a plurality of difference values. If the plurality of difference values includes a target difference value exceeding a target time length, the computer device takes the sum of the start time of the first data channel and the target time length as the starting point of the second time range, and takes the start time of the second data channel as the ending point of the second time range, the first data channel and the second data channel being the former data channel and the latter data channel of the two data channels corresponding to the target difference value, respectively.

[0075] Optionally, the computer device obtains the start time and the acquisition time length of each data channel from the data channel header of each data channel, respectively. For example, the computer device obtains the start time of the Kth data channel in the continuous data channels of the ith wireless node as t ik , and it should be noted that the channel length of the continuous data channels of any abnormal wireless node, i.e. the acquisition time length, is equal, i.e. the target time length, denoted as L.

[0076] Wherein, the computer device compares the adjacent two data channels t ik -t ik-1 and the target time length L, if t ik -t ik-1 >L, it is determined that the seismic data is missing, and the time range of the missing seismic data is t ik-1 +L to t ik ; and similarly, the computer device can count the time range of the missing seismic data of each abnormal wireless node, which can be represented as t i1-1 +L to t i1 , t i2-1 +L to t i2 , …, t iv-1 +L to t iv , …, t iq-1 +L to t iq . Wherein, k represents the serial number of the data channel of the ith wireless node, v represents the index of the missing seismic data segment in the continuous data channels of the ith wireless node, q represents the total number of missing seismic data segments in the continuous data channels of the ith wireless node, t iv-1 +L represents the start time of the missing seismic data segment in the continuous data channels of the ith wireless node, and t iv represents the end time of the missing seismic data segment in the continuous data channels of the ith wireless node.

[0077] In the embodiment, since the acquisition time lengths of the two adjacent data are both the target time length, and if the difference between the start times of the two adjacent data is greater than the target time length, it indicates that no seismic data is acquired in the time range greater than the target time length, that is, the seismic data loss exists, and therefore the sum of the start time of the first data and the target time length is taken as the start point of the second time range, and the start time of the second data is taken as the end point of the second time range, thereby improving the accuracy of determining the second time range.

[0078] 103. The computer device acquires the firing time of the plurality of first shot points.

[0079] The firing time is the time when the first shot point fires the seismic wave.

[0080] 104. The computer device determines the target firing time range corresponding to the seismic data based on the first time range and the second time range.

[0081] The target firing time range includes a target firing time range corresponding to each of the plurality of abnormal wireless nodes; and the computer device determines the target firing time range corresponding to the seismic data based on the first time range and the second time range, including the following steps (1)-(2).

[0082] (1) The computer device determines the target time range corresponding to each of the plurality of abnormal wireless nodes from the first time range and the second time range corresponding to each of the plurality of abnormal wireless nodes, and the length of the target time range exceeds the length threshold.

[0083] If the first time range and the second time range corresponding to the abnormal wireless node are both multiple, the computer device sorts and merges the first time range and the second time range corresponding to the abnormal wireless node respectively to obtain target time ranges corresponding to multiple abnormal wireless nodes respectively. For each abnormal wireless node, the computer device sorts the multiple first time ranges and the multiple second time ranges according to the start time of each time range; then, the computer device compares the start time and the end time of any two adjacent time ranges, if the sampling points corresponding to the end time of the former time range and the start time of the latter time range are only different by one sampling point, the computer device merges the two time ranges to obtain a third time range, the start time and the end time of the third time range are the start time of the former time range and the end time of the latter time range respectively, and the first time range, the second time range and the third time range are all invalid time ranges; the computer device takes the time range with a time length exceeding a target threshold in the invalid time range as the target time range; it should be noted that since the first time range, the second time range and the third time range are all multiple, the target time range is multiple, that is, multiple target time ranges corresponding to multiple abnormal wireless nodes.

[0084] The target time length can be set and changed as needed; optionally, the target time length is 120s. It should be noted that in the process of collecting seismic data, random errors or system errors will cause the appearance of zero value data or the loss of seismic data for a short time, and in this embodiment, by taking the time range with a time length exceeding the target threshold as the target time range, the accuracy of determining the target time range is ensured.

[0085] In some embodiments, the computer device obtains the total time length of collecting seismic data of any abnormal wireless node based on the data channel head of the abnormal wireless node, and then determines the difference between the total time length and the target time range to obtain an effective time range of the seismic data, which is the range of recording effective seismic data; and then it is convenient to subsequently check the data integrity of the seismic data based on the effective time range.

[0086] Referring to Figure 2 , Figure 2A flowchart for determining an effective time range is provided for the embodiments of the present application. First, the computer device determines a first time range corresponding to zero-value data and a second time range corresponding to missing seismic data for any abnormal wireless node, respectively; then the computer device collates and merges the first time range and the second time range to count the overall invalid time range of the abnormal wireless node; the computer device further obtains the target time range corresponding to the abnormal wireless node by screening out time ranges with a time length greater than a target time length; finally, the computer device obtains the effective time range corresponding to the abnormal wireless node based on the target time range and the total time length of the seismic data collected by the abnormal wireless node.

[0087] Referring to Figure 3 , Figure 3 A statistical information diagram of a wireless node is provided for the embodiments of the present application, which shows the statistical information of the integrity check of the data file corresponding to 25 wireless nodes, including the line number, point number, file name of recorded seismic data, file size (KB), sampling frequency (ms), start time of data recording, end time of data recording, east coordinate, north coordinate, total time length (ms) of effective time range, total time length (ms) of invalid time range, and data status of a plurality of wireless nodes, etc. Referring to Figure 4 , Figure 4 A statistical result diagram of the integrity check of an abnormal wireless node is provided for the embodiments of the present application, which is a statistical report of the wireless node with line number 2878 and point number 3046, which includes the basic information, effective time range situation and invalid time range situation of the wireless node as shown in Figure 3 The effective time range situation includes the total time length of the effective time range, the number of effective time range segments, the start time and end time of each effective time range segment, and the time length of each effective time range segment, etc. The invalid time range situation includes the total time length of the invalid time range, the number of invalid time range segments, the start time and end time of each invalid time range segment, and the time length of each invalid time range segment, etc.

[0088] (2) The computer device determines the target shooting time range corresponding to a plurality of abnormal wireless nodes based on the target time range corresponding to the plurality of abnormal wireless nodes, respectively.

[0089] The computer device determines a target excitation time range corresponding to each of the plurality of abnormal wireless nodes based on a target time range corresponding to each of the plurality of abnormal wireless nodes, including the following steps: the computer device obtains a time length adjustment parameter, the time length adjustment parameter being used to correct the target time range; the computer device determines a difference between a start point of the target time range corresponding to each of the plurality of abnormal wireless nodes and the time length adjustment parameter; and the computer device determines a time range composed of the difference corresponding to each of the plurality of abnormal wireless nodes and an end point of the target time range as the target excitation time range corresponding to each of the plurality of abnormal wireless nodes.

[0090] The time length adjustment parameter includes a length ls of a controllable seismic source scanning signal and a length lc of a seismic data correlation record in a construction parameter of a seismic data acquisition area. c Optionally, if the ith target time range corresponding to the wth abnormal wireless node is represented as Ts iw and Te iw , the computer device determines the target excitation time range as Ts iw -ls-lc to Te iw . Optionally, the length of the controllable seismic source scanning signal is 9000 ms, and the length of the seismic data correlation record is 6000 ms.

[0091] It should be noted that the time range of the seismic data acquisition is the time range after removing the time length adjustment parameter when the wireless node acquires the seismic data; and in this embodiment, the target time range is corrected by the time length adjustment parameter, so that the corrected target excitation time range matches the real invalid time range, and then the subsequent screening of the shot point based on the target excitation time range can improve the accuracy of the shot point screening.

[0092] 105. The computer device determines abnormal shot points in the plurality of first shot points based on the excitation time of the plurality of first shot points and the target excitation time range.

[0093] In some embodiments, the computer device determines abnormal shot points in the plurality of first shot points based on the excitation time of the plurality of first shot points, a first time range and a second time range, including the following steps (1)-(3).

[0094] (1) The computer device determines a plurality of second shot points corresponding to a plurality of abnormal wireless nodes from the plurality of first shot points.

[0095] The computer device determines a plurality of second shot points corresponding to a plurality of abnormal wireless nodes from the plurality of first shot points based on an observation system arrangement relationship of the seismic data; and the observation system arrangement relationship is used to record the correspondence between the wireless nodes and the shot points. Optionally, the computer device determines a receiving line number rl i, whether the point number rpi is located in the line number sl h , the point number is sr h , the gun index number is si h , if yes, the first shot point is taken as the second shot point corresponding to the abnormal wireless node. Similarly, the computer device determines a plurality of second shot points corresponding to a plurality of abnormal wireless nodes.

[0096] (2) The computer device determines a target shot point corresponding to each abnormal wireless node from a plurality of second shot points corresponding to the abnormal wireless node based on a target excitation time range corresponding to the abnormal wireless node and an excitation time of the second shot point, wherein the excitation time of the target shot point matches the target excitation time range.

[0097] wherein, for any abnormal wireless node, the computer device takes the second shot point with an excitation time within the target excitation time range as the target shot point corresponding to the abnormal wireless node; wherein the target shot point corresponding to any abnormal wireless node is multiple.

[0098] Optionally, the computer device obtains a shot point list of the target shot point based on the target shot point corresponding to each abnormal wireless node, and records the line number, point number and index number of the target shot point in the shot point list to mark the target shot point; for example, the target shot points are represented as (sl1, sp1, si1), (sl2, sp2, si2), …, (sl g , sp g , si g ), …, (sl z , sp z , si z ), g represents the serial number of the target shot point, z represents the total number of the target shot point, sl g represents the line number of the target shot point with serial number g, sp g represents the point number of the target shot point with serial number g, si g represents the index number of the target shot point with serial number g, and the index number is used to represent the number of times of exciting the seismic wave by the shot point.

[0099] (3) The computer device determines an abnormal shot point from the target shot point corresponding to each abnormal wireless node.

[0100] In some embodiments, the computer device determines an abnormal shot point from the target shot point corresponding to each abnormal wireless node, including the following steps: the computer device determines the ratio between the number of abnormal wireless nodes corresponding to any target shot point and the total number of wireless nodes corresponding to the target shot point to obtain the abnormal channel proportion of the target shot point; and the computer device determines the target shot point as an abnormal shot point when the abnormal channel proportion is greater than a target proportion.

[0101] Among them, for any line number sl, the computer equipment h The dot number is sp h The index number is si h The target firing point, if any line number among multiple abnormal wireless nodes is rl i The dot number is rp i If an anomalous wireless node is located within the array of the target shot point, meaning the seismic data collected by the anomalous wireless node was generated after the target shot point triggered seismic waves, then the computer equipment will record the number of anomalous wireless nodes corresponding to the target shot point as cn. h Add 1, cn h This indicates the number of abnormal wireless nodes corresponding to the target gun point with serial number h. Furthermore, the computer device uses the number of multiple wireless nodes located within the array range of the target gun point as the total number of wireless nodes corresponding to that target gun point.

[0102] The computer equipment determines the ratio between the number of abnormal wireless nodes corresponding to any target gun point and the number of wireless nodes corresponding to the target gun point, and obtains the abnormal channel ratio of the target gun point, which is achieved by the following formula (1).

[0103] rb g =cb g / ch g (1);

[0104] Among them, rb g cb represents the proportion of abnormal trajectory for target shot point with sequence number g. g This indicates the number of abnormal radio nodes corresponding to the target gun point with index g, ch g This represents the total number of wireless nodes corresponding to the target gun point with serial number g.

[0105] Where, if rb g >rb th The computer equipment then determines that the target firing point is an abnormal firing point. g rb represents the anomaly proportion of the target shot point with sequence number g. th This indicates the target percentage; the target percentage can be set and changed as needed; optionally, the target percentage is 0.02%, which corresponds to 341 abnormal shot points.

[0106] In this embodiment, by determining the second firing points corresponding to multiple abnormal wireless nodes, and then based on the target firing time corresponding to the abnormal wireless nodes and the firing time of the second firing points, the target firing point in the second firing points that matches the target firing time of the abnormal wireless nodes can be determined, thus improving the efficiency of determining the target firing points.

[0107] SeeFigure 5 , Figure 5 An information schematic diagram of an abnormal shot point provided by an embodiment of the present application; the diagram includes information of 25 abnormal shot points, including 7 kinds of information of line number, point number, index number, shot point firing time (TB), east coordinate, north coordinate and corresponding number of abnormal wireless nodes of each of the abnormal shot points.

[0108] In the embodiment, if the abnormal trace ratio exceeds the target ratio, it indicates that the probability of the shot point firing seismic wave anomaly is relatively large, and then the target shot point is taken as the abnormal shot point, thereby improving the accuracy of the shot point screening.

[0109] Referring to Figure 6 , Figure 6 Another flowchart for determining a target shot point provided by an embodiment of the present application. In the flowchart, the computer device determines a third shot point in the plurality of first shot points based on the firing time of the plurality of first shot points and the target firing time range of the plurality of abnormal wireless nodes, and the firing time of the third shot point is located in any target firing time range. Then, the computer device determines whether the abnormal wireless node corresponding to the target firing time range is the wireless node corresponding to the third shot point based on the observation system arrangement relationship of the seismic data, and if yes, the computer device determines that the third shot point is the target shot point. Then, the computer device screens any target shot point based on the ratio between the number of abnormal wireless nodes corresponding to the target shot point and the total number of wireless nodes corresponding to the target shot point and the target ratio.

[0110] Referring to Figure 7 , Figure 7 A flowchart of another shot point screening method provided by an embodiment of the present application. In the flowchart, the computer device imports a data file of a plurality of wireless nodes, and based on the data file of the plurality of wireless nodes, an invalid time range is obtained by statistics. Then, a target shot point affected by the invalid time range is screened. Then, abnormal shot points are screened by statistics of abnormal trace ratios of each target shot point.

[0111] In the embodiment of the present application, the computer device checks the seismic data collected by the wireless nodes to timely find zero value data and missing seismic data in the seismic data, and screens the shot points with abnormal trace ratios based on the invalid time range determined based on the zero value data and the missing seismic data, thereby timely finding the seismic data collected by the wireless nodes with quality problems and the abnormal shot points, and facilitating timely processing and rectification of the seismic data with quality problems and the wireless nodes, and facilitating timely re-shooting, thereby ensuring the quality and efficiency of the collected seismic data.

[0112] The embodiment of the present application provides a shot point screening method, since the first time range is a time range corresponding to zero value data in the seismic data, and the second time range is a time range of missing seismic data, and then the corresponding target excitation time range is obtained by reverse deduction based on the first time range and the second time range, so that the abnormal shot point in the multiple shot points can be determined by matching the excitation time of the shot point and the target excitation time range, and the accuracy of the shot point screening is improved.

[0113] The embodiment of the present application further provides a shot point screening device, referring to Figure 8 , the device comprises:

[0114] The first acquisition module 801 is used for acquiring seismic data, and the seismic data is generated after a plurality of first shot points excite seismic waves;

[0115] The first determination module 802 is used for determining a first time range and a second time range corresponding to the seismic data respectively, the first time range is a time range corresponding to zero value data in the seismic data, and the second time range is a time range of missing seismic data;

[0116] The second acquisition module 803 is used for acquiring excitation time of the multiple first shot points;

[0117] The second determination module 804 is used for determining a target excitation time range corresponding to the seismic data based on the first time range and the second time range;

[0118] The third determination module 805 is used for determining an abnormal shot point in the multiple first shot points based on the excitation time of the multiple first shot points and the target excitation time range.

[0119] In some embodiments, the seismic data comprises a plurality of data channels connected in sequence, the acquisition time length of the plurality of data channels is a target time length, and the device further comprises:

[0120] The fourth determination module is used for determining a difference value of a start time corresponding to any two adjacent data channels in the plurality of data channels, to obtain a plurality of difference values;

[0121] The fifth determination module is used for, if the plurality of difference values include a target difference value exceeding the target time length, taking a sum of the start time of the first data channel and the target time length as a start point of the second time range, taking the start time of the second data channel as an end point of the second time range, and the first data channel and the second data channel are a previous data channel and a next data channel in the two data channels corresponding to the target difference value respectively.

[0122] In some embodiments, the seismic data comprises seismic data collected by a plurality of wireless nodes, the target shot time range comprises target shot time ranges corresponding to a plurality of abnormal wireless nodes respectively, the abnormal wireless node is a node in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless node has zero value data and loss; the third determination module 805 is configured to:

[0123] determine a plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of second shot points based on the target shot time ranges corresponding to the plurality of abnormal wireless nodes respectively and shot times of the plurality of second shot points;

[0124] determine a plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of second shot points based on the target shot time ranges corresponding to the plurality of abnormal wireless nodes respectively and shot times of the plurality of second shot points;

[0125] determine an abnormal shot point from the target shot points corresponding to the plurality of abnormal wireless nodes respectively.

[0126] In some embodiments, the third determination module 805 is configured to:

[0127] determine an abnormal trace ratio of a target shot point by determining a ratio between a number of abnormal wireless nodes corresponding to the target shot point and a total number of wireless nodes corresponding to the target shot point;

[0128] determine the target shot point as an abnormal shot point when the abnormal trace ratio is greater than a target ratio.

[0129] In some embodiments, the seismic data comprises seismic data collected by a plurality of wireless nodes, the first time range and the second time range comprise first time ranges and second time ranges corresponding to a plurality of abnormal wireless nodes respectively, the target shot time range comprises target shot time ranges corresponding to the plurality of abnormal wireless nodes respectively, the abnormal wireless node is a node in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless node has zero value data and loss; the second determination module 804 is configured to:

[0130] determine target time ranges corresponding to the plurality of abnormal wireless nodes respectively from the first time ranges and the second time ranges corresponding to the plurality of abnormal wireless nodes respectively, the target time ranges having a length exceeding a length threshold;

[0131] determine target shot time ranges corresponding to the plurality of abnormal wireless nodes respectively based on the target time ranges corresponding to the plurality of abnormal wireless nodes respectively.

[0132] In some embodiments, the second determination module 804 is configured to:

[0133] obtain a length adjustment parameter, the length adjustment parameter being used to correct the target time range;

[0134] determining a difference between a start point of a target time range corresponding to each of the plurality of abnormal wireless nodes and the time length adjustment parameter;

[0135] composing a time range composed of the difference corresponding to each of the plurality of abnormal wireless nodes and an end point of the target time range as a target firing time range corresponding to each of the plurality of abnormal wireless nodes.

[0136] Embodiments of the present application provide a shot point screening device. Since the first time range is a time range corresponding to zero value data in seismic data, and the second time range is a time range of missing seismic data, and then based on the first time range and the second time range, the target firing time range corresponding thereto is obtained by reverse deduction, and the abnormal shot points in the plurality of shot points can be determined by matching the firing time of the shot points and the target firing time range, and the accuracy of shot point screening is improved.

[0137] In some embodiments, the computer device is configured as a terminal; Figure 9 A structure block diagram of a terminal 900 provided by an exemplary embodiment of the present application is shown. The terminal 900 can be a portable mobile terminal, such as a smart phone, a tablet computer, an MP3 (Moving Picture Experts Group Audio Layer III) player, an MP4 (Moving Picture Experts Group Audio Layer IV) player, a notebook computer or a desktop computer. The terminal 900 can also be referred to as a user equipment, a portable terminal, a laptop terminal, a desktop terminal, and other names.

[0138] Generally, the terminal 900 includes a processor 901 and a memory 902.

[0139] The processor 901 can include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 901 can be implemented in the form of at least one of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), a PLA (Programmable Logic Array), and the like. The processor 901 can also include a main processor and a co-processor. The main processor is a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit). The co-processor is a low-power consumption processor for processing data in a standby state. In some embodiments, the processor 901 can be integrated with a GPU (Graphics Processing Unit) for rendering and drawing content to be displayed by the display screen. In some embodiments, the processor 901 can further include an AI (Artificial Intelligence) processor for processing machine learning-related computing operations.

[0140] The memory 902 can include one or more computer-readable storage media that can be non-transitory. The memory 902 can also include a high-speed random access memory, and a nonvolatile memory such as one or more disk storage devices, flash storage devices. In some embodiments, the non-transitory computer-readable storage medium in the memory 902 is used to store at least one program code for being executed by the processor 901 to implement the shot screening method provided by the method embodiments in the present application.

[0141] In some embodiments, the terminal 900 can also optionally include a peripheral device interface 903 and at least one peripheral device. The processor 901, the memory 902, and the peripheral device interface 903 can be connected through a bus or a signal line. Each peripheral device can be connected to the peripheral device interface 903 through a bus, a signal line, or a circuit board. Specifically, the peripheral device includes at least one of a radio frequency circuit 904, a display screen 905, a camera component 906, an audio circuit 907, a positioning component 908, and a power supply 909.

[0142] The peripheral interface 903 can be used to connect at least one I / O (Input / Output) related peripheral device to the processor 901 and the memory 902. In some embodiments, the processor 901, the memory 902 and the peripheral interface 903 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 901, the memory 902 and the peripheral interface 903 can be implemented on a separate chip or circuit board, and the present embodiments are not limited in this regard.

[0143] The radio frequency circuit 904 is configured to receive and send RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuit 904 communicates with communication networks and other communication devices through electromagnetic signals. The radio frequency circuit 904 converts electrical signals into electromagnetic signals for transmission, or converts electromagnetic signals received into electrical signals. Optionally, the radio frequency circuit 904 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and the like. The radio frequency circuit 904 can communicate with other terminals through at least one wireless communication protocol. The wireless communication protocol includes but is not limited to the World Wide Web, a metropolitan area network, an intranet, various generations of mobile communication networks (2G, 3G, 4G and 5G), a wireless local area network, and / or a WiFi (Wireless Fidelity) network. In some embodiments, the radio frequency circuit 904 can also include NFC (Near Field Communication) related circuitry, and the present application is not limited in this regard.

[0144] The display screen 905 is configured to display a UI (User Interface). The UI can include graphics, text, icons, video, and any combination thereof. When the display screen 905 is a touch display screen, the display screen 905 is further configured to capture touch signals on or above the surface of the display screen 905. The touch signals can be input to the processor 901 as control signals for processing. In this case, the display screen 905 can also be configured to provide virtual buttons and / or virtual keyboard, also known as soft buttons and / or soft keyboard. In some embodiments, the display screen 905 can be one, disposed on the front panel of the terminal 900; in other embodiments, the display screen 905 can be at least two, respectively disposed on different surfaces of the terminal 900 or in a folding design; in other embodiments, the display screen 905 can be a flexible display screen, disposed on a curved surface or a folding surface of the terminal 900. Even, the display screen 905 can also be disposed in an irregular shape other than a rectangle, i.e., a special-shaped screen. The display screen 905 can be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), etc.

[0145] The camera assembly 906 is configured to capture images or videos. Optionally, the camera assembly 906 includes a front camera and a rear camera. Typically, the front camera is disposed on the front panel of the terminal, and the rear camera is disposed on the back of the terminal. In some embodiments, the rear camera is at least two, which are any one of a main camera, a depth-of-field camera, a wide-angle camera, and a telephoto camera, to realize the background blur function by fusing the main camera and the depth-of-field camera, the panoramic shooting and VR (Virtual Reality) shooting function by fusing the main camera and the wide-angle camera, or other fusion shooting functions. In some embodiments, the camera assembly 906 can further include a flash. The flash can be a single-color-temperature flash or a dual-color-temperature flash. The dual-color-temperature flash refers to the combination of a warm light flash and a cold light flash, which can be used for light compensation under different color temperatures.

[0146] The audio circuit 907 can include a microphone and a speaker. The microphone is used to collect sound waves of a user and an environment, and convert the sound waves into an electrical signal input to the processor 901 for processing, or input to the radio frequency circuit 904 to realize voice communication. For the purpose of stereo sound collection or noise reduction, the microphone can be multiple, which are respectively arranged at different parts of the terminal 900. The microphone can also be an array microphone or an omnidirectional collection type microphone. The speaker is used to convert an electrical signal from the processor 901 or the radio frequency circuit 904 into sound waves. The speaker can be a conventional diaphragm speaker, or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only can the electrical signal be converted into a sound wave audible to humans, but also can be converted into a sound wave inaudible to humans for ranging purposes. In some embodiments, the audio circuit 907 can also include a headphone jack.

[0147] The positioning component 908 is used to position the current geographic location of the terminal 900 to realize navigation or LBS (Location Based Service). The positioning component 908 can be a positioning component based on the GPS (Global Positioning System) of the United States, the Beidou system of China or the Galileo system of Russia.

[0148] The power supply 909 is used to supply power to each component in the terminal 900. The power supply 909 can be alternating current, direct current, disposable battery or rechargeable battery. When the power supply 909 includes a rechargeable battery, the rechargeable battery can be a wired charging battery or a wireless charging battery. The wired charging battery is a battery charged through a wired line, and the wireless charging battery is a battery charged through a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0149] In some embodiments, the terminal 900 further includes one or more sensors 190. The one or more sensors 190 include but are not limited to an acceleration sensor 911, a gyroscope sensor 912, a pressure sensor 913, a fingerprint sensor 914, an optical sensor 915 and a proximity sensor 916.

[0150] The acceleration sensor 911 can detect the acceleration size in three coordinate axes of the coordinate system established by the terminal 900. For example, the acceleration sensor 911 can be used to detect the components of the gravitational acceleration in three coordinate axes. The processor 901 can control the display screen 905 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 911. The acceleration sensor 911 can also be used for game or user motion data collection.

[0151] The gyroscope sensor 912 can detect the body direction and rotation angle of the terminal 900, and can collect 3D motions of the user on the terminal 900 in cooperation with the acceleration sensor 911. The processor 901 can implement the following functions according to the data collected by the gyroscope sensor 912: motion sensing (such as changing the UI according to the tilt operation of the user), image stabilization when shooting, game control, and inertial navigation.

[0152] The pressure sensor 913 can be arranged on the side frame of the terminal 900 and / or the lower layer of the display screen 905. When the pressure sensor 913 is arranged on the side frame of the terminal 900, the grip signal of the user on the terminal 900 can be detected, and the left-hand or right-hand recognition or shortcut operation can be performed by the processor 901 according to the grip signal collected by the pressure sensor 913. When the pressure sensor 913 is arranged on the lower layer of the display screen 905, the controllable control on the UI interface can be controlled by the processor 901 according to the pressure operation of the user on the display screen 905. The controllable control includes at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0153] The fingerprint sensor 914 is used to collect the fingerprint of the user, and the identity of the user can be recognized by the processor 901 according to the fingerprint collected by the fingerprint sensor 914, or by the fingerprint sensor 914 according to the collected fingerprint. When the identity of the user is recognized as a trusted identity, the processor 901 authorizes the user to perform related sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, payment, and changing settings. The fingerprint sensor 914 can be arranged on the front, back or side of the terminal 900. When the terminal 900 is provided with a physical button or a manufacturer's logo, the fingerprint sensor 914 can be integrated with the physical button or the manufacturer's logo.

[0154] The optical sensor 915 is used to collect the ambient light intensity. In one embodiment, the processor 901 can control the display brightness of the display screen 905 according to the ambient light intensity collected by the optical sensor 915. Specifically, when the ambient light intensity is high, the display brightness of the display screen 905 is increased; when the ambient light intensity is low, the display brightness of the display screen 905 is decreased. In another embodiment, the processor 901 can also dynamically adjust the shooting parameters of the camera assembly 906 according to the ambient light intensity collected by the optical sensor 915.

[0155] The proximity sensor 916, also referred to as a distance sensor, is usually arranged on the front panel of the terminal 900. The proximity sensor 916 is used to collect the distance between the user and the front of the terminal 900. In an embodiment, when the proximity sensor 916 detects that the distance between the user and the front of the terminal 900 gradually decreases, the display screen 905 is switched from the bright screen state to the screen-off state under the control of the processor 901; when the proximity sensor 916 detects that the distance between the user and the front of the terminal 900 gradually increases, the display screen 905 is switched from the screen-off state to the bright screen state under the control of the processor 901.

[0156] Those skilled in the art can understand that, Figure 9 The structure shown in the figure does not constitute a limitation on the terminal 900, and can include more or fewer components than shown, or combine certain components, or adopt a different arrangement of components.

[0157] The embodiment of the present application also provides a computer readable storage medium, and at least one program code is stored in the computer readable storage medium. The at least one program code is loaded and executed by a processor to implement the shot screening method of any of the above implementation manners.

[0158] The embodiment of the present application also provides a computer program product, comprising a computer program, and the computer program is executed by a processor to implement the shot screening method of any of the above implementation manners.

[0159] In some embodiments, the computer program product related to the embodiment of the present application can be deployed on one computer device for execution, or on multiple computer devices located in one place for execution, or on multiple computer devices distributed in multiple places and interconnected through a communication network for execution. The multiple computer devices distributed in multiple places and interconnected through a communication network can constitute a blockchain system.

[0160] The embodiment of the present application provides a shot screening method. Since the first time range is a time range corresponding to zero value data in the seismic data, and the second time range is a time range of missing seismic data, the corresponding target excitation time range is obtained by reverse deduction based on the first time range and the second time range. Therefore, the abnormal shot point in the multiple shot points can be determined by matching the excitation time of the shot point with the target excitation time range, and the accuracy of shot screening is improved.

[0161] The above is only an optional embodiment of the present application, and does not limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method of shotpoint screening, characterized by, The method comprises: acquiring seismic data generated after seismic waves are excited by a plurality of first shot points, the seismic data comprising seismic data collected by a plurality of wireless nodes; determining a first time range and a second time range corresponding to the seismic data respectively, the first time range being a time range corresponding to zero-value data in the seismic data, and the second time range being a time range of missing seismic data, the first time range and the second time range respectively comprising a first time range and a second time range corresponding to a plurality of abnormal wireless nodes respectively; acquiring excitation times of the plurality of first shot points; based on the first time range and the second time range, determining a target excitation time range corresponding to the seismic data, the target excitation time range comprising a target excitation time range corresponding to the plurality of abnormal wireless nodes respectively, the abnormal wireless nodes being nodes in the plurality of wireless nodes, and the seismic data collected by the abnormal wireless nodes having zero-value data and missing conditions; based on the excitation times of the plurality of first shot points and the target excitation time range, determining abnormal shot points in the plurality of first shot points; wherein, based on the first time range and the second time range, determining a target excitation time range corresponding to the seismic data comprises: determining a target time range corresponding to the plurality of abnormal wireless nodes respectively from the first time range and the second time range corresponding to the plurality of abnormal wireless nodes respectively, the target time range having a duration exceeding a duration threshold; and determining the target excitation time range corresponding to the plurality of abnormal wireless nodes respectively based on the target time range corresponding to the plurality of abnormal wireless nodes respectively; based on the excitation times of the plurality of first shot points and the target excitation time range, determining abnormal shot points in the plurality of first shot points comprises: determining a plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of first shot points; determining target shot points corresponding to the plurality of abnormal wireless nodes respectively from the plurality of second shot points corresponding to the plurality of abnormal wireless nodes respectively based on the target excitation time range corresponding to the plurality of abnormal wireless nodes respectively and the excitation times of the plurality of second shot points, the excitation time of the target shot points matching the target excitation time range; and determining the abnormal shot points from the target shot points corresponding to the plurality of abnormal wireless nodes respectively.

2. The method of claim 1, wherein, The seismic data comprises a plurality of data channels connected in sequence, the collection duration of each of the plurality of data channels being a target duration, and the determination process of the second time range comprises: determining a difference value of the start time corresponding to any two adjacent data channels in the plurality of data channels, to obtain a plurality of difference values; if the plurality of difference values include a target difference value exceeding the target duration, taking the sum of the start time of a first data channel and the target duration as the start point of the second time range, and taking the start time of a second data channel as the end point of the second time range, the first data channel and the second data channel being a preceding data channel and a subsequent data channel in the two data channels corresponding to the target difference value respectively.

3. The method of claim 1, wherein, The step of determining the abnormal gun point from the target gun points corresponding to the plurality of abnormal wireless nodes includes: The ratio of the number of abnormal wireless nodes corresponding to any target gun point to the total number of wireless nodes corresponding to the target gun point is determined to obtain the abnormal trajectory ratio of the target gun point. If the proportion of abnormal paths is greater than the target proportion, the target firing point is determined to be an abnormal firing point.

4. The method of claim 1, wherein, The step of determining the target excitation time range corresponding to each of the multiple abnormal wireless nodes based on the target time range corresponding to each of the multiple abnormal wireless nodes includes: Obtain duration adjustment parameters, which are used to correct the target time range; Determine the difference between the starting point of the target time range corresponding to each of the multiple abnormal wireless nodes and the duration adjustment parameter; The time range formed by the difference between the values ​​corresponding to the multiple abnormal wireless nodes and the end point of the target time range is taken as the target excitation time range corresponding to the multiple abnormal wireless nodes.

5. A shotpoint screening device, characterized in that The device includes: The first acquisition module is used to acquire seismic data, which is generated by seismic waves being excited by multiple first shot points, and includes seismic data collected by multiple wireless nodes. The first determining module is used to determine the first time range and the second time range corresponding to the earthquake data respectively. The first time range is the time range corresponding to the zero value data in the earthquake data, and the second time range is the time range of lost earthquake data. The first time range and the second time range respectively include the first time range and the second time range corresponding to multiple abnormal wireless nodes respectively. The second acquisition module is used to acquire the firing time of the plurality of first shot points; The second determining module is used to determine the target excitation time range corresponding to the seismic data based on the first time range and the second time range. The target excitation time range includes the target excitation time range corresponding to multiple abnormal wireless nodes. The abnormal wireless nodes are nodes among the multiple wireless nodes, and the seismic data collected by the abnormal wireless nodes contains zero-value data and has data loss. The third determining module is used to determine the abnormal shot point among the plurality of first shot points based on the firing time of the plurality of first shot points and the firing time range of the target; The second determining module is used to determine, from the first time range and the second time range corresponding to the plurality of abnormal wireless nodes respectively, the target time range having a duration exceeding a duration threshold; and to determine the target excitation time range corresponding to the plurality of abnormal wireless nodes based on the target time ranges corresponding to the plurality of abnormal wireless nodes respectively. The third determining module is configured to determine, from the multiple first shot points, multiple second shot points corresponding to the multiple abnormal wireless nodes respectively; determine, from the multiple second shot points corresponding to the multiple abnormal wireless nodes respectively, target shot points corresponding to the multiple abnormal wireless nodes respectively based on target firing time ranges corresponding to the multiple abnormal wireless nodes respectively and firing times of the multiple second shot points, the firing time of the target shot point matching the target firing time range; and determine the abnormal shot point from the target shot points corresponding to the multiple abnormal wireless nodes respectively.

6. A computer device, comprising: The computer device comprises one or more processors and one or more memories, and at least one piece of program code is stored in the one or more memories. The at least one piece of program code is loaded and executed by the one or more processors to implement the shot point screening method according to any one of claims 1 to 4.

7. A computer readable storage medium characterized in that, The storage medium stores at least one piece of program code, and the at least one piece of program code is loaded and executed by the processor to implement the shot point screening method according to any one of claims 1 to 4.

8. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the shot point screening method according to any one of claims 1 to 4.

Images