Micro-motion detection method for pre-generating triangular array position

Abstract

A micro-motion detection method for pre-generating the position of a triangular array comprises the following steps: firstly, determining a reference point in a target area, and then placing a range finder on the reference point; then, reference position data is obtained through measurement of a range finder, a triangular grid is generated, and then numbering is carried out to obtain a triangular point number and a triangular serial number; generating the distance from the range finder to the triangular point number; then enabling the difference between the distance between the seismograph and the distance measuring instrument and the distance of the triangular point number to meet the requirement, namely completing the seismograph layout of a triangular array, and obtaining a micro-motion signal by the seismograph; moving the triangular array until all micro-motion signals are obtained; processing the micro-motion signal to obtain a shear wave velocity three-dimensional structure; the positions of the seismographs are determined through the multiple range finders, so that accurate array position parameters can be obtained, the correlation degree of the parameters and the micro-motion signals is high, and the underground structure inversion precision is high. Therefore, according to the design, array position parameters are accurate, and the underground structure inversion precision is high.

Description

Technical Field

[0001] The invention relates to a micro-motion detection method, belongs to the field of micro-motion detection, and in particular to a micro-motion detection method for pre-generating a triangular array position. Background Art

[0002] Microseismic detection is a geophysical detection method that inverts underground structures based on surface microseismic signals (i.e., weak vibrations generated by nature and human activities).

[0003] A Chinese patent application with application number 202311698530.3 and application date of December 11, 2023 discloses a micro-seismic exploration method and system, the system comprising an acquisition module, a data collection module and a processing module; the acquisition module and the processing module are both communicatively connected to the data collection module; the number of the acquisition modules is multiple; although the design pre-processes the vibration wave data collected in real time by the acquisition module through the data collection module and sends it to the processing module, the processing module then extracts the surface wave dispersion curve based on the vibration wave data, and then inverts the surface wave dispersion curve to obtain the phase velocity and depth map of the underground structure of the area to be measured, it still has the following defects: This design arranges the seismographs into an array through measuring ropes, but the measuring ropes have low accuracy when used, which will lead to inaccurate position parameters of the array, resulting in low accuracy of the underground structure inverted based on the array position parameters and micro-motion signals.

[0004] The information disclosed in this background technology section is only intended to increase the understanding of the overall background of the application, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to ordinary technicians in this field. Summary of the invention

[0005] The purpose of the present invention is to overcome the defects and problems of inaccurate array position parameters and low accuracy of underground structure inversion in the prior art, and to provide a micro-motion detection method for pre-generated triangular array positions with accurate array position parameters and high accuracy of underground structure inversion.

[0006] To achieve the above objectives, the technical solution of the present invention is: A method for detecting micro-motion by pre-generating a triangular array position, the method comprising the following steps: Step 1: First determine the scope of the target area, then select multiple reference points within the scope of the target area, and then place a rangefinder on each reference point in turn; Step 2: First, the rangefinder measures the boundary points of the target area, and then obtains the reference position data, and then sets the coordinate system according to the reference position data, wherein the coordinate system includes the origin, the horizontal axis and the vertical axis, and then generates a triangular mesh within the origin, the horizontal axis and the vertical axis, and then sequentially numbers each intersection position and the midpoint position of the triangle in the triangular mesh within the reference position data to obtain a plurality of triangular point numbers, and then numbers each triangle in the triangular mesh within the reference position data to obtain a triangle serial number, and each triangle serial number corresponds to four triangular point numbers; Step 3: First, generate the distance from each rangefinder to each triangulation point number based on the triangle grid and coordinate system to obtain the point number distance from each triangulation point number to each rangefinder. Then, collect all the point number distances in the order of the triangle number and generate a distance document. Step 4: first place the seismographs to be deployed in the target area, and then start the preliminary deployment process: first arbitrarily select a triangulation point number where no seismograph is deployed, and then use multiple rangefinders to measure the distance from the same seismograph in turn, and then compare the multiple distances with the corresponding point number distances in turn. When the distance is greater than the corresponding point number distance, the seismograph is moved toward the corresponding rangefinder so that the difference between the distance and the corresponding point number distance meets the requirement. When the distance is less than the corresponding point number distance, the seismograph is moved away from the corresponding rangefinder so that the difference between the distance and the corresponding point number distance meets the requirement. When the differences between all distances and the corresponding point number distances meet the requirement, the seismograph is deemed to be located at the selected triangulation point number, and a preliminary deployment process is completed; then select any other triangulation point number in the triangulation sequence corresponding to the aforementioned triangulation point number for the preliminary deployment process, until all triangulation point numbers in the triangulation sequence are deployed with seismographs, and then all seismographs are monitored simultaneously to obtain multiple micro-motion signals; Step 5: Repeat step 4 until all the micro motion signals corresponding to the triangle numbers are obtained; Step 6: First, the micro-motion signal is processed to obtain the three-dimensional structure of the shear wave velocity below the target area, and this method ends.

[0007] In the fourth step, the initial deployment process is to start the initial deployment process for multiple seismographs.

[0008] In the sixth step, the micromotion signal is processed to obtain the three-dimensional structure of the shear wave velocity below the target area by first extracting a dispersion curve from the micromotion signal and then inverting the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area.

[0009] In the sixth step, the dispersion curve is extracted from the micro-motion signal by first processing the micro-motion signal by a spatial autocorrelation method to obtain a spatial autocorrelation coefficient, then fitting the spatial autocorrelation coefficient with a Bessel function, and then extracting the dispersion curve.

[0010] In the sixth step, the processing of the micro-motion signal by the spatial autocorrelation method is to pre-process the micro-motion signal first to extract the effective micro-motion signal from the micro-motion signal, and then process the effective micro-motion signal by the spatial autocorrelation method.

[0011] In the sixth step, the inversion of the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area is to invert the dispersion curve through a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area.

[0012] In the sixth step, the dispersion curve is inverted by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area. The surface wave phase velocity and frequency are first obtained by the dispersion curve, and then the surface wave phase velocity and frequency are substituted into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity. Then, the surface wave phase velocity is converted into the shear wave velocity, and then the relationship between the depth and the shear wave velocity is obtained, that is, the three-dimensional structure of the shear wave velocity below the target area is obtained.

[0013] In the sixth step, the relationship between depth and shear wave velocity is obtained by first processing the relationships between multiple depths and shear wave velocities by interpolation method, and then obtaining the three-dimensional structure of shear wave velocity in the target area.

[0014] In the first step, the selecting of a plurality of reference points within the scope of the target area is selecting greater than or equal to three reference points within the scope of the target area.

[0015] In the first step, selecting three or more reference points within the scope of the target area is selecting three reference points within the scope of the target area.

[0016] Compared with the prior art, the present invention has the following beneficial effects: 1. In a method for detecting micro-motions of pre-generating triangular array positions of the present invention, the method comprises the following steps: the first step: first determine the target area, then select the reference point, and then place a rangefinder on the reference point; the second step: first measure the boundary points of the target area with the rangefinder to obtain reference position data, then generate a coordinate system based on the data, then generate a triangular mesh in the coordinate system, then number each intersection point and midpoint position of the triangles in the target area to obtain a plurality of triangular point numbers, then number each triangle to obtain a plurality of triangular serial numbers, each triangular serial number corresponds to four triangular point numbers; the third step: first generate the point number distance from each rangefinder to the triangular point number, and then obtain a distance file; the fourth step: first place the seismograph In the target area, the layout process is started again: first select a triangulation point number, three rangefinders respectively measure the distance from the point number to the seismograph, and move the seismograph so that the difference between the distance and the point number distance meets the requirement. Meeting the requirement means that the seismograph is located at a triangulation point number, and then repeat the above-mentioned layout process for other triangulation point numbers corresponding to the triangulation point number, that is, the layout of a group of seismographs is completed, and then the seismograph monitors to obtain micro-motion signals; the fifth step: repeat the fourth step until the micro-motion signals of all triangulation numbers are obtained; the sixth step: perform data processing on the micro-motion signals to obtain the three-dimensional structure of the shear wave velocity below the target area, and then judge the three-dimensional structure below the target area according to the three-dimensional structure of the shear wave velocity. The advantages of the present invention also include: First, a triangular grid is generated, and then the triangulation point number, triangulation sequence number and point number distance are obtained. The location of the triangulation point number is the location of the seismograph. Then the distance from the triangulation point number to the seismograph is measured by a distance meter. When the difference between the distance and the point number distance is small, the seismograph is located at the triangulation point number, so the seismograph is located at the location. The location of the seismograph is determined by multiple distance meters, so the location of the seismograph is accurate, and the location parameters of the array are accurate. The array position can be associated with the measured micro-motion signal, so the accuracy of the inverted underground structure is high. Second point: compared with using a measuring rope to lay out a seismograph, the seismograph of the present invention is laid out more accurately, and a distance document can be directly generated without the need for subsequent measurement and recording of the position of the seismograph, thus saving manpower and material costs, and thus having a higher work efficiency; and when laying out a seismograph, it is only necessary to move the seismograph and then use a distance meter to measure, so that position correction can be completed, and thus it is more convenient to use; Third point: multiple groups of seismographs can be deployed at the same time, and multiple groups of seismographs can monitor micro-motion signals at the same time to improve monitoring efficiency; Therefore, the array position parameters of the present invention are accurate and the underground structure inversion accuracy is high.

[0017] 2. In the present invention, in a method for detecting micromotions that pregenerates the position of a triangular array, in the sixth step, the micromotion signals are first preprocessed, and then the micromotion signals are processed by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, and then the spatial autocorrelation coefficient is fitted with the Bessel function to obtain the dispersion curve, and then the surface wave phase velocity and frequency are obtained from the dispersion curve, and then the relationship between the depth and the shear wave velocity is converted by the half-wavelength empirical formula, that is, the three-dimensional structure of the shear wave velocity below the target area is obtained, which is the process of inverting the underground structure; the micromotion signals in a target area are interdependent, so the spatial autocorrelation method is used to process the micromotion signals, and the spatial autocorrelation method can efficiently process a large number of micromotion signals, so the underground structure can be efficiently obtained. Therefore, the present invention is highly efficient.

[0018] 3. In the micro-motion detection method for pre-generating the position of a triangular array of the present invention, in the first step, the number of the reference points is three, and rangefinders are arranged on the three reference points respectively, and then the position of the seismograph is measured by the three rangefinders to determine whether the seismograph is located on the triangulation point number; the number of rangefinders can be set as needed, but three rangefinders can obtain a relatively accurate position of the seismograph. Therefore, the positioning accuracy of the present invention is high. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Figure 1 It is a structural schematic diagram of the present invention.

[0020] Figure 2 It is a schematic diagram of the application of the present invention.

[0021] Figure 3 yes Figure 2 Schematic diagram of the structure of the triangular mesh.

[0022] Figure 4 yes Figure 2 Schematic diagram of the structure of the middle triangle serial number.

[0023] Figure 5 It is a schematic diagram of the case where the distance in Example 1 is smaller than the point distance.

[0024] Figure 6 It is a schematic diagram of the case where the distance in Example 1 is greater than the point distance.

[0025] Figure 7 It is a schematic diagram that the difference between the distance and the point distance in Example 1 meets the requirement.

[0026] Figure 8 yes Figure 2 Schematic diagram of the structure of the medium triangle array.

[0027] Fig. 9 yes Figure 1 Schematic diagram of a mid-range document.

[0028] Fig.10 yes Fig. 9 A magnified view of the midpoint distance.

[0029] Fig.11 yes Figure 1 Schematic diagram of the three-dimensional structure of shear wave velocity in .

[0030] In the figure: target area A, reference point A1, rangefinder B, coordinate system C, origin C1, horizontal axis C2, vertical axis C3, triangular grid C4, triangular point number C5, triangular sequence number C6, seismograph D, triangular array D1. DETAILED DESCRIPTION

[0031] The present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0032] See also Figure 1 — Fig.11 , a method for detecting micro-motion by pre-generating the position of a triangular array, the method comprising the following steps: Step 1: First determine the scope of the target area A, then select multiple reference points A1 within the scope of the target area A, and then place a rangefinder B on each reference point A1 in turn; Step 2: First, the rangefinder B measures the boundary points of the target area A, and then obtains the reference position data, and then sets the coordinate system C according to the reference position data, wherein the coordinate system C includes the origin C1, the horizontal axis C2 and the vertical axis C3, and then generates a triangular mesh C4 within the origin C1, the horizontal axis C2 and the vertical axis C3, and then sequentially numbers each intersection position and the midpoint position of the triangle in the triangular mesh C4 located in the reference position data to obtain a plurality of triangular point numbers C5, and then numbers each triangle in the triangular mesh C4 located in the reference position data to obtain a triangle serial number C6, and each triangle serial number C6 corresponds to four triangle point numbers C5; Step 3: First, based on the triangle grid C4 and the coordinate system C, generate the distance from each rangefinder B to each triangle point number C5 to obtain the point number distance corresponding to each triangle point number C5 to each rangefinder B. Then, collect all the point number distances in the order of the triangle number C6, and then generate a distance document. Step 4: First place the seismograph D to be deployed in the target area A, and then start the preliminary deployment process: first arbitrarily select a triangulation point C5 where no seismograph D is deployed, and then use multiple rangefinders B to measure the distance to the same seismograph D in turn, and then compare the multiple distances with the corresponding point distances in turn. When the distance is greater than the corresponding point distance, move the seismograph D toward the corresponding rangefinder B so that the difference between the distance and the corresponding point distance meets the requirement. When the distance is less than the corresponding point distance, move the seismograph D away from the corresponding point distance. The distance meter B is moved in the direction of the distance meter B so that the difference between the distance and the corresponding point number distance meets the requirement. When the difference between all the distances and the corresponding point number distances meets the requirement, the seismograph D is regarded as being located at the selected triangulation point number C5, that is, a preliminary deployment process is completed; then any other triangulation point number C5 in the triangulation sequence number C6 corresponding to the aforementioned triangulation point number C5 is selected for the preliminary deployment process, until all the triangulation point numbers C5 in the triangulation sequence number C6 are deployed with the seismograph D, and then all the seismographs D are monitored simultaneously to obtain multiple micro-motion signals; Step 5: Repeat step 4 until all micro-motion signals corresponding to triangle number C6 are obtained; Step 6: First, process the micro-motion signal to obtain the three-dimensional structure of the shear wave velocity below the target area A, and this method ends.

[0033] In the fourth step, the initial deployment process is to start the initial deployment process for multiple seismographs D.

[0034] In the sixth step, the micromotion signal is processed to obtain the three-dimensional structure of the shear wave velocity below the target area A by first extracting a dispersion curve from the micromotion signal and then inverting the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area A.

[0035] In the sixth step, the dispersion curve is extracted from the micro-motion signal by first processing the micro-motion signal by a spatial autocorrelation method to obtain a spatial autocorrelation coefficient, then fitting the spatial autocorrelation coefficient with a Bessel function, and then extracting the dispersion curve.

[0036] In the sixth step, the processing of the micro-motion signal by the spatial autocorrelation method is to pre-process the micro-motion signal first to extract the effective micro-motion signal from the micro-motion signal, and then process the effective micro-motion signal by the spatial autocorrelation method.

[0037] In the sixth step, the inversion of the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area A is to invert the dispersion curve through a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area A.

[0038] In the sixth step, the dispersion curve is inverted by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area A. The surface wave phase velocity and frequency are first obtained by the dispersion curve, and then the surface wave phase velocity and frequency are substituted into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity. The surface wave phase velocity is then converted into the shear wave velocity, and the relationship between the depth and the shear wave velocity is obtained, that is, the three-dimensional structure of the shear wave velocity below the target area A is obtained.

[0039] In the sixth step, the relationship between depth and shear wave velocity is obtained by first processing the relationships between multiple depths and shear wave velocities by interpolation method, and then obtaining the three-dimensional structure of shear wave velocity of target area A.

[0040] In the first step, the selecting of a plurality of reference points A1 within the scope of the target area A is selecting greater than or equal to three reference points A1 within the scope of the target area A.

[0041] In the first step, selecting three or more reference points A1 within the target area A is selecting three reference points A1 within the target area A.

[0042] The supplementary description of the present invention is as follows: The measuring rope described in the present invention refers to: a measuring tool for measuring length, generally a steel wire rope with a scale buckle; one end of the steel wire rope is pulled, and the length of the pulled steel wire rope is read through the scale buckle, thereby achieving the effect of measuring distance; during manual measurement, the error of the measuring rope is large, and the measuring rope needs to be pulled once the seismograph D is moved, so the operation is relatively cumbersome.

[0043] The scope of the target area A is determined as follows: the scope of the target area A is determined according to actual detection requirements, the boundary point of the target area A is a point with a large turning point at the edge of the target area A, and then the reference position data is generated according to the position of the boundary point. Figure 2 The purple dotted line in the middle is the virtual boundary line of the target area A automatically generated according to the reference position data; when generating the coordinate system C from the reference position data, the reference position data is placed in the first quadrant as far as possible to facilitate subsequent steps; the triangles in the triangular mesh C4 that are outside the reference position data or on the virtual boundary of the reference position data are not numbered.

[0044] The present invention selects multiple reference points A1 within the scope of the target area A so that there is no physical obstruction between the rangefinder B on the reference point A1 and all the triangulation point numbers C5, so the laser emitted by the rangefinder B can be reflected by the seismograph D, thereby obtaining the distance from the rangefinder B to the seismograph D.

[0045] Embodiment 1: See also Figure 1 — Fig.11, a method for detecting micro-motion by pre-generating the position of a triangular array, the method comprising the following steps: Step 1: First determine the scope of the target area A, then select multiple reference points A1 within the scope of the target area A, and then place a rangefinder B on each reference point A1 in turn; Step 2: First, the rangefinder B measures the boundary points of the target area A, and then obtains the reference position data, and then sets the coordinate system C according to the reference position data, wherein the coordinate system C includes the origin C1, the horizontal axis C2 and the vertical axis C3, and then generates a triangular mesh C4 within the origin C1, the horizontal axis C2 and the vertical axis C3, and then sequentially numbers each intersection position and the midpoint position of the triangle in the triangular mesh C4 located in the reference position data to obtain a plurality of triangular point numbers C5, and then numbers each triangle in the triangular mesh C4 located in the reference position data to obtain a triangle serial number C6, and each triangle serial number C6 corresponds to four triangle point numbers C5; Step 3: First, based on the triangle grid C4 and the coordinate system C, generate the distance from each rangefinder B to each triangle point number C5 to obtain the point number distance corresponding to each triangle point number C5 to each rangefinder B. Then, collect all the point number distances in the order of the triangle number C6, and then generate a distance document. Step 4: First place the seismograph D to be deployed in the target area A, and then start the preliminary deployment process: first arbitrarily select a triangulation point C5 where no seismograph D is deployed, and then use multiple rangefinders B to measure the distance to the same seismograph D in turn, and then compare the multiple distances with the corresponding point distances in turn. When the distance is greater than the corresponding point distance, move the seismograph D toward the corresponding rangefinder B so that the difference between the distance and the corresponding point distance meets the requirement. When the distance is less than the corresponding point distance, move the seismograph D away from the corresponding point distance. The distance meter B is moved in the direction of the distance meter B so that the difference between the distance and the corresponding point number distance meets the requirement. When the difference between all the distances and the corresponding point number distances meets the requirement, the seismograph D is regarded as being located at the selected triangulation point number C5, that is, a preliminary deployment process is completed; then any other triangulation point number C5 in the triangulation sequence number C6 corresponding to the aforementioned triangulation point number C5 is selected for the preliminary deployment process, until all the triangulation point numbers C5 in the triangulation sequence number C6 are deployed with the seismograph D, and then all the seismographs D are monitored simultaneously to obtain multiple micro-motion signals; Step 5: Repeat step 4 until all micro-motion signals corresponding to triangle number C6 are obtained; Step 6: First, process the micro-motion signal to obtain the three-dimensional structure of the shear wave velocity below the target area A, and this method ends.

[0046] In the fourth step, the initial deployment process is to start the initial deployment process for multiple seismographs D.

[0047] Embodiment 2: The basic content is the same as that of Example 1, except that: See also Figure 1 — Fig.11 In the sixth step, the processing of the micromotion signal to obtain the three-dimensional structure of the shear wave velocity below the target area A is to first extract the dispersion curve from the micromotion signal, and then invert the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area A. In the sixth step, the extraction of the dispersion curve from the micromotion signal is to first process the micromotion signal by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, and then fit the spatial autocorrelation coefficient with the Bessel function, and then extract the dispersion curve. In the sixth step, the processing of the micromotion signal by the spatial autocorrelation method is to first pre-process the micromotion signal to extract the effective micromotion signal from the micromotion signal, and then process the effective micromotion signal by the spatial autocorrelation method. In the sixth step, the inversion of the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area A is to invert the dispersion curve by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area A. In the sixth step, the dispersion curve is inverted by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area A. The surface wave phase velocity and frequency are first obtained by the dispersion curve, and then the surface wave phase velocity and frequency are substituted into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity, and then the surface wave phase velocity is converted into the shear wave velocity, and then the relationship between the depth and the shear wave velocity is obtained, that is, the three-dimensional structure of the shear wave velocity below the target area A is obtained. In the sixth step, the relationship between the depth and the shear wave velocity is first processed by the interpolation method for the relationship between multiple depths and the shear wave velocity, and then the three-dimensional structure of the shear wave velocity of the target area A is obtained.

[0048] When applied, in the sixth step, the micromotion signal is first preprocessed to extract the effective micromotion signal in the micromotion signal, and then the effective micromotion signal is processed by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, and then the spatial autocorrelation coefficient is fitted with the Bessel function, and then the dispersion curve is obtained, and then the surface wave phase velocity and frequency are obtained through the dispersion curve, and then the surface wave phase velocity and frequency are introduced into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity, and then the surface wave phase velocity is converted into the shear wave velocity, and then the relationship between the depth and the shear wave velocity is obtained, that is, a group of micromotion signal processing of triangle number C6 is completed; then all the micromotion signals of triangle number C6 are processed by the interpolation method, and then the three-dimensional structure of the shear wave velocity underground in the target area A is obtained, and then the underground structure of the target area A is judged by the three-dimensional structure of the shear wave velocity underground in the target area A.

[0049] Embodiment 3: The basic content is the same as that of Example 1, except that: See also Figure 1 — Fig.10 In the first step, the multiple reference points A1 selected in the target area A are selected at least three reference points A1 in the target area A. In the first step, the three reference points A1 selected in the target area A are selected at least three reference points A1 in the target area A.

[0050] When applied, in the first step, three reference points A1 are selected first, and then a rangefinder B is placed on each reference point A1; measuring the distance to a seismograph D by three rangefinders B is sufficient to determine the position of the seismograph D, even if the seismograph D is located on the triangulation point number C5.

[0051] The above description is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiment. Any equivalent modifications or changes made by ordinary technicians in this field based on the contents disclosed by the present invention should be included in the protection scope recorded in the claims.

Claims

1. A method for detecting micro-motion by pre-generating the position of a triangular array, characterized in that: The method comprises the following steps: Step 1: First determine the scope of the target area (A), then select multiple reference points (A1) within the scope of the target area (A), and then place a rangefinder (B) on each reference point (A1) in turn; Step 2: First, the rangefinder (B) measures the boundary points of the target area (A) to obtain reference position data, and then the coordinate system (C) is set according to the reference position data, wherein the coordinate system (C) includes an origin (C1), a horizontal axis (C2) and a vertical axis (C3), and then a triangular mesh (C4) is generated within the origin (C1), the horizontal axis (C2) and the vertical axis (C3), and then each intersection position and a midpoint position of a triangle in the triangular mesh (C4) located in the reference position data are numbered in sequence to obtain a plurality of triangular point numbers (C5), and then each triangle in the triangular mesh (C4) located in the reference position data is numbered to obtain a triangular serial number (C6), and each triangular serial number (C6) corresponds to four triangular point numbers (C5); Step 3: First, generate the distance from each rangefinder (B) to each triangle point number (C5) based on the triangle grid (C4) and the coordinate system (C) to obtain the point number distance from each triangle point number (C5) to each rangefinder (B). Then, collect all the point number distances in the order of the triangle number (C6) and generate a distance document. Step 4: Place the seismograph (D) to be deployed in the target area (A), and then start the preliminary deployment process: first randomly select a triangulation point number (C5) where no seismograph (D) is deployed, and then use multiple rangefinders (B) to measure the distance to the same seismograph (D) in turn, and then compare the multiple distances with the corresponding point number distances in turn. When the distance is greater than the corresponding point number distance, move the seismograph (D) toward the corresponding rangefinder (B) so that the difference between the distance and the corresponding point number distance meets the requirement. When the distance is less than the corresponding point number distance, move the seismograph (D) away from the corresponding The distance meter (B) moves in the direction so that the difference between the distance and the corresponding point number distance meets the requirement. When the difference between all the distances and the corresponding point number distances meets the requirement, the seismograph (D) is considered to be located at the selected triangulation point number (C5), and a preliminary deployment process is completed; then any other triangulation point number (C5) in the triangulation sequence number (C6) corresponding to the aforementioned triangulation point number (C5) is selected for the preliminary deployment process, until all the triangulation point numbers (C5) in the triangulation sequence number (C6) are deployed on the seismographs (D), and then all the seismographs (D) are monitored simultaneously to obtain multiple micro-motion signals; Step 5: Repeat step 4 until all the micro-motion signals corresponding to the triangle numbers (C6) are obtained; Step 6: First, the micro-motion signal is processed to obtain the three-dimensional structure of the shear wave velocity below the target area (A), and this method ends.

2. The method for detecting micro-motions of pre-generated triangular array positions according to claim 1, characterized in that: In the fourth step, the initial deployment process is to start the initial deployment process for multiple seismographs (D).

3. A method for detecting micro-motions by pre-generating a triangular array position according to claim 1 or 2, characterized in that: In the sixth step, the micromotion signal is processed to obtain the three-dimensional structure of the shear wave velocity below the target area (A) by first extracting a dispersion curve from the micromotion signal and then inverting the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area (A).

4. The method for detecting micro-motions of pre-generated triangular array positions according to claim 3, characterized in that: In the sixth step, the dispersion curve is extracted from the micro-motion signal by first processing the micro-motion signal by a spatial autocorrelation method to obtain a spatial autocorrelation coefficient, then fitting the spatial autocorrelation coefficient with a Bessel function, and then extracting the dispersion curve.

5. The method for detecting micro-motions of pre-generated triangular array positions according to claim 4, characterized in that: In the sixth step, the processing of the micro-motion signal by the spatial autocorrelation method is to pre-process the micro-motion signal first to extract the effective micro-motion signal from the micro-motion signal, and then process the effective micro-motion signal by the spatial autocorrelation method.

6. The method for detecting micro-motions of pre-generated triangular array positions according to claim 5, characterized in that: In the sixth step, the inversion of the dispersion curve to obtain the three-dimensional structure of the shear wave velocity below the target area (A) is to invert the dispersion curve through a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area (A).

7. The method for detecting micro-motions of pre-generated triangular array positions according to claim 6, characterized in that: In the sixth step, the dispersion curve is inverted by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area (A). The surface wave phase velocity and frequency are first obtained by the dispersion curve, and then the surface wave phase velocity and frequency are substituted into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity. The surface wave phase velocity is then converted into the shear wave velocity, and the relationship between the depth and the shear wave velocity is obtained, thereby obtaining the three-dimensional structure of the shear wave velocity below the target area (A).

8. The method for detecting micro-motions of pre-generated triangular array positions according to claim 7, characterized in that: In the sixth step, the relationship between depth and shear wave velocity is obtained by first processing the relationships between multiple depths and shear wave velocities by interpolation method, and then obtaining the three-dimensional structure of shear wave velocity of the target area (A).

9. A method for detecting micro-motions by pre-generating the position of a triangular array according to claim 1 or 2, characterized in that: In the first step, the selecting of a plurality of reference points (A1) within the scope of the target area (A) is selecting greater than or equal to three reference points (A1) within the scope of the target area (A).

10. The method for detecting micro-motion by pre-generating the position of a triangular array according to claim 9, characterized in that: In the first step, selecting three or more reference points (A1) within the target area (A) is selecting three reference points (A1) within the target area (A).

Images