A method for accurately positioning node seismic instruments

Abstract

The invention discloses a node seismic instrument precise positioning method, which reads the positioning data of the node seismic instrument, performs preprocessing and data cleaning on it, and obtains the positioning data set of the node seismic instrument; detects the node seismic instrument coordinate points in the neighborhood of any unmarked node seismic instrument coordinate point in the positioning data set of the node seismic instrument, obtains outliers and core points, all outliers form an outlier point set, and all core points form a core point set; searches the core point set until all node seismic instrument coordinate points are searched and marked; removes the outlier point set from the node seismic instrument positioning data set, obtains a new object set, and calculates the geographic center point of the node seismic instrument. Compared with the prior art, the invention can improve the positioning accuracy and operation efficiency of the node instrument during field acquisition.

Description

Technical Field

[0001] The present invention relates to the technical field of geological exploration, and particularly to a precise positioning method for node seismic instruments. Background Art

[0002] With the continuous deepening of oil and gas exploration and development, geophysical exploration construction is facing increasingly complex surface conditions. As a wireless and self - collecting geophysical exploration instrument, the node seismic instrument has solved this problem to a certain extent. After the node instrument is powered on after being deployed in the field and the built - in GPS module completes satellite searching, the node instrument starts data acquisition. After the node instrument starts data acquisition, it will record its coordinate information in real time and store it in the node instrument. Using a dedicated mobile phone APP, through Bluetooth communication, the data stored in the node instrument can be read and sent back to the base through the 4G / 5G network. When the node seismic instrument is collecting data in the field, its built - in GPS module records coordinate information in real time. However, due to the limitations of the accuracy of civilian GPS and the interference of complex field environments, the GPS coordinates collected by the node instrument are scattered. However, when quality control is performed on the node instrument, a unique and accurate position of the node instrument deployed in the field is required.

[0003] Chinese Patent Application No. CN202011201851.4 discloses a method and device for matching stake numbers of on - land wireless node seismic instruments. First, it converts coordinate data from the geodetic coordinate system to the plane rectangular coordinate system; then it calculates the geometric center point CenterPt of all coordinates, the distance dist_i between each point and the center point, and the average distance mean_dist between each point and the center point in turn. If dist_i < mean_dist, it is marked as 1. If dist_i > mean_dist, it is marked as 0. Then, based on the points marked as 1, a new geometric center point newCenterPt is obtained; the new geometric center point is compared with the previous geometric center point. If it is greater than the set threshold, the above steps are repeated until the distance between the new geometric center point and the previous geometric center point is less than the set threshold, and then it ends. Although this method can also obtain the unique and accurate position of the node instrument deployed in the field, this method has the following disadvantages:

[0004] (1) After reading the data, no data cleaning is performed, resulting in too large deviations in the basic data. For null values, 0 values, and positioning data with obvious large deviations, they are not excluded, resulting in too many subsequent operation iterations, affecting the subsequent operation speed and accuracy;

[0005] (2) This algorithm is too sensitive to isolated points with low density, and additional threshold parameters need to be set, resulting in low calculation accuracy, low operation efficiency, and the need for human intervention, etc.;

[0006] (3) Since this algorithm directly converts the geodetic coordinate system into a rectangular coordinate system and then obtains the geometric center point of the rectangular coordinate system, it is not as accurate as obtaining the geographic center point in the original geodetic coordinate system, which affects the final positioning accuracy.

[0007] Based on the above reasons, it is urgent to develop a new node instrument precise positioning method to improve the positioning accuracy and computing efficiency of node instruments during field collection. Summary of the invention

[0008] The purpose of the present invention is to solve the deficiencies in the prior art and to provide a method for accurately positioning a node seismic instrument so as to improve the positioning accuracy and computing efficiency of the node instrument during field acquisition.

[0009] To achieve the above object, the present invention is implemented according to the following technical solutions:

[0010] A node seismic instrument precise positioning method comprises the following steps:

[0011] (1) Read the positioning data of the node seismic instrument, perform preprocessing and data cleaning on it, and obtain the positioning data set of the node seismic instrument;

[0012] (2) Detecting node seismic instrument coordinate points in the neighborhood of any unmarked node seismic instrument coordinate point in the node seismic instrument positioning data set to obtain outlier points and core points, wherein all outlier points constitute an outlier point set, and all core points constitute a core point set; searching the core point set until all node seismic instrument coordinate points have been searched and marked;

[0013] (3) The outlier set is removed from the node seismic instrument location data set to obtain a new object set, and the geographic center point of the node seismic instrument is calculated.

[0014] Furthermore, in the step (1), the preprocessing and data cleaning are specifically as follows: after reading the node seismic instrument positioning data, intercept the positioning data within the blasting time period, extract 1 / 10 of the data in time average, remove abnormal values ​​and null values, and organize the longitude, latitude, and device number into a format suitable for subsequent calculations.

[0015] Furthermore, the step (2) specifically includes:

[0016] 1) Convert all positioning data in the positioning data set from geodetic coordinates to plane rectangular coordinates;

[0017] 2) Set the neighborhood radius Eps and threshold minPts;

[0018] 3) Create empty sets L and M;

[0019] 4) Select any unmarked node seismic instrument coordinate point in the positioning data set and mark it as searched;

[0020] 5) Search for node seismic instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point within a radius of 5; if the number of node seismic instrument coordinate points within the range is less than minPts, mark the point as an outlier and add it to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, mark the point as a core point and add all coordinate points in the neighborhood of the point to the set M;

[0021] 6) In the set M, select any unmarked node seismic instrument coordinate point and mark it as searched;

[0022] 7) Within a radius of 5, search for node instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point; if the number of node seismic instrument coordinate points within the range is less than minPts, mark the point as an outlier and add it to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, mark the point as a core point, and add all node seismic instrument coordinate points in the neighborhood of the point to the set M until all node seismic instrument coordinate points are marked.

[0023] Furthermore, in step (3), calculating the geographic center point of the node seismic instrument specifically includes:

[0024] 1) Convert plane rectangular coordinates into geodetic coordinates;

[0025] 2) Convert the longitude and latitude of the geodetic coordinates into radians, and then further convert them into positions in the coordinate system to obtain the horizontal, vertical and vertical coordinates; add up the horizontal, vertical and vertical coordinates of each node seismic instrument, divide them by the total number of coordinates, and calculate the average values ​​of the horizontal, vertical and vertical coordinates; convert the horizontal, vertical and vertical coordinates into longitude and latitude again to obtain the precise longitude and latitude coordinates of the node seismic instrument.

[0026] Preferably, the neighborhood radius Eps is set to 5, and the threshold minPts is set to 1% of the amount of all positioning data in the positioning data set.

[0027] Compared with the prior art, after reading the coordinate data, the present invention first performs data cleaning to eliminate abnormal values, reduce subsequent operations, and improve operation efficiency; further, the present invention adopts an improved clustering outlier detection algorithm to reduce the influence of isolated points on the initial center point, and no longer needs to set a threshold, thereby improving the overall operation efficiency and positioning accuracy; finally, the precise longitude and latitude coordinates of the node seismic instrument are obtained in the geodetic coordinate system to improve its positioning accuracy; in summary, the present invention can improve the positioning accuracy and operation efficiency of the node instrument during field collection. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 This is a flow chart for the precise positioning of the node seismic instrument of the present invention. DETAILED DESCRIPTION

[0029] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments. The specific embodiments described herein are only used to explain the present invention and are not used to limit the invention.

[0030] like Figure 1 As shown, this embodiment specifically discloses a method for accurately locating a node seismic instrument, and the specific positioning process is as follows:

[0031] First, use a mobile phone APP to read the positioning data of the node seismic instrument, and send the data to the server through the network. The server performs preprocessing and data cleaning. The preprocessing and data cleaning methods in this embodiment are as follows: after reading the positioning data of the node seismic instrument, intercept the positioning data within the blasting time period, extract 1 / 10 of the data by time average, remove abnormal values ​​and null values, and organize the longitude, latitude, and device number into a format suitable for subsequent calculations; obtain the positioning data set of the node seismic instrument; the purpose of data cleaning is to reduce subsequent calculations and improve calculation efficiency. It should be noted that subsequent calculations are all performed on the server.

[0032] Secondly, the server detects node seismic instrument coordinate points in the neighborhood of any unmarked node seismic instrument coordinate point in the node seismic instrument positioning data set, obtains outlier points and core points, all outlier points form an outlier point set, and all core points form a core point set; the core point set is searched until all node seismic instrument coordinate points have been searched and marked;

[0033] 1) Convert all positioning data in the positioning data set from geodetic coordinates to plane rectangular coordinates;

[0034] 2) Set the neighborhood radius Eps and threshold minPts. In actual use, in order to ensure the subsequent positioning accuracy and efficiency, the neighborhood radius Eps is set to 5, and the threshold minPts is set to 1% of all positioning data in the positioning data set;

[0035] 3) Create empty sets L and M;

[0036] 4) Select any unmarked node seismic instrument coordinate point P in the positioning data set and mark it as searched;

[0037] 5) Search for node seismic instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point P within a radius of 5; if the number of node seismic instrument coordinate points within the range is less than minPts, mark the point as an outlier and add it to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, mark the point as a core point and add all coordinate points in the neighborhood of the point to the set M;

[0038] 6) In the set M, select any unmarked node seismic instrument coordinate point Q and mark it as searched;

[0039] 7) Within a radius of 5, search for node instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point Q; if the number of node seismic instrument coordinate points within the range is less than minPts, mark the point as an outlier and add it to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, mark the point as a core point, and add all node seismic instrument coordinate points in the neighborhood of the point to the set M until all node seismic instrument coordinate points are marked.

[0040] Then, the outlier set is removed from the node seismic instrument location data set to obtain a new object set, and the geographic center point of the node seismic instrument is calculated. The calculation of the geographic center point of the node seismic instrument specifically includes:

[0041] 1) Convert plane rectangular coordinates into geodetic coordinates;

[0042] 2) Convert the longitude and latitude of the geodetic coordinates into radians, and then further convert them into positions in the coordinate system to obtain the horizontal, vertical and vertical coordinates; add up the horizontal, vertical and vertical coordinates of each node seismic instrument, divide by the total number of coordinates, and calculate the average of the horizontal, vertical and vertical coordinates; convert the horizontal, vertical and vertical coordinates into longitude and latitude to obtain the precise longitude and latitude coordinates of the node seismic instrument. Obtain the precise longitude and latitude coordinates of the node seismic instrument in the geodetic coordinate system to improve its positioning accuracy.

[0043] The technical solution of the present invention is not limited to the above-mentioned specific embodiments. All technical variations made according to the technical solution of the present invention fall within the protection scope of the present invention.

Claims

1. A method for accurately positioning a node seismic instrument, characterized in that: The following steps are involved: (1) Read the positioning data of the node seismic instrument, and perform preprocessing and data cleaning on it to obtain the positioning data set of the node seismic instrument; the preprocessing and data cleaning are specifically as follows: after reading the positioning data of the node seismic instrument, intercept the positioning data within the blasting time period, extract 1 / 10 of the data in this period according to the time average, remove abnormal values ​​and null values, and organize the longitude, latitude, and device number into a format suitable for subsequent calculations; (2) Detecting node seismic instrument coordinate points in the neighborhood of any unmarked node seismic instrument coordinate point in the node seismic instrument positioning data set to obtain outlier points and core points, wherein all outlier points constitute an outlier point set, and all core points constitute a core point set; searching the core point set until all node seismic instrument coordinate points have been searched and marked; 1) Convert all positioning data in the positioning data set from geodetic coordinates to plane rectangular coordinates; 2) Set the neighborhood radius Eps and threshold minPts; 3) Create empty sets L and M; 4) Select any unmarked node seismic instrument coordinate point in the positioning data set and mark it as searched; 5) Search for node seismic instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point within a radius of 5; if the number of node seismic instrument coordinate points within the range is less than minPts, mark the point as an outlier and add it to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, mark the point as a core point and add all coordinate points in the neighborhood of the point to the set M; 6) In the set M, select any unmarked node seismic instrument coordinate point and mark it as searched; 7) Searching for node instrument coordinate points in the neighborhood of the unmarked node seismic instrument coordinate point within a radius of 5; If the number of node seismic instrument coordinate points within the range is less than minPts, the point is marked as an outlier and added to the set L; if the number of node seismic instrument coordinate points within the range is not less than minPts, the point is marked as a core point, and all node seismic instrument coordinate points in the neighborhood of the point are added to the set M until all node seismic instrument coordinate points have been marked; (3) Remove the outlier set from the node seismic instrument location data set to obtain a new object set, and calculate the geographic center point of the node seismic instrument; The geographic center points of the computing node seismic instruments include: 1) Convert plane rectangular coordinates into geodetic coordinates; 2) Convert the longitude and latitude of the geodetic coordinates into radians, and then further convert them into positions in the coordinate system to obtain the horizontal, vertical and vertical coordinates; add up the horizontal, vertical and vertical coordinates of each node seismic instrument, divide them by the total number of coordinates, and calculate the average values ​​of the horizontal, vertical and vertical coordinates; convert the horizontal, vertical and vertical coordinates into longitude and latitude again to obtain the precise longitude and latitude coordinates of the node seismic instrument.

2. The node seismic instrument precise positioning method according to claim 1, characterized in that: The neighborhood radius Eps is set to 5, and the threshold minPts is set to 1% of the amount of all positioning data in the positioning data set.

Images