A micro-motion detection method for pre-generating triangular array positions

Abstract

A microtremor detection method for pre-generating triangular array positions, comprising the following steps: first, determining a reference point in a target area, and then placing a range finder on the reference point; then measuring reference position data by the range finder, and then generating a triangular grid, and then numbering to obtain a triangular point number and a triangular serial number; then generating a distance from the range finder to the triangular point number; then making the distance between the seismograph and the range finder and the distance of the triangular point number meet the requirements, that is, completing the seismograph arrangement of a triangular array, and then obtaining microtremor signals by the seismograph; moving the triangular array until all microtremor signals are obtained; and then processing the microtremor signals to obtain a three-dimensional structure of shear wave velocity; the positions of the seismographs are determined by multiple range finders, so that accurate array position parameters can be obtained, the parameters have a high correlation degree with the microtremor signals, and therefore the inversion of the underground structure has a high precision. Therefore, the array position parameters are accurate, and the inversion of the underground structure has a high precision.

Description

TECHNICAL FIELD

[0001] The application relates to a microtremor detection method, and belongs to the microtremor detection field. BACKGROUND

[0002] Microtremor detection is a geophysical detection method based on microtremor signals (i.e. weak vibrations generated by natural phenomena and human activities) to invert underground structures.

[0003] A microtremor exploration method and system are disclosed in Chinese Patent No. 202311698530.3, filed on December 11, 2023. The system includes a collection module, a data collection module, and a processing module. The collection module and the processing module are in communication connection with the data collection module. The number of collection modules is multiple. Although this design sends the vibration wave data collected by the data collection module in real time to the processing module after preprocessing, the processing module extracts the surface wave dispersion curve based on the vibration wave data, then performs surface wave dispersion curve inversion to obtain the phase velocity and depth map of the underground structure of the region to be measured, but still has the following defects:

[0004] This design arranges the seismograph into an array by measuring the rope, but the accuracy of the rope in use is low, which will cause the position parameters of the array to be inaccurate, resulting in low accuracy of the inverted underground structure based on the position parameters of the array and the microtremor signals.

[0005] The information disclosed in this background section is intended only to increase an understanding of the general context of the present application, and it should not be taken as an acknowledgement or any form of suggestion that this information forms prior art that is already widely known in the art. SUMMARY

[0006] The purpose of the present application is to overcome the defects and problems of inaccurate array position parameters and low accuracy of inverted underground structures in the prior art, and to provide a microtremor detection method for pre-generating triangular array positions with accurate array position parameters and high accuracy of inverted underground structures.

[0007] To achieve the above purpose, the technical solution of the present application is:

[0008] A microtremor detection method for pre-generating triangular array positions, the method comprising the following steps:

[0009] Step 1: First, determine the range of the target area, then select multiple reference points in the range of the target area, and then place a range finder on each reference point in turn;

[0010] Second step: first, the range finder 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, the coordinate system includes the origin, the horizontal axis and the vertical axis, and then generates a triangular grid in the range of the origin, the horizontal axis and the vertical axis, and then sequentially numbers each intersection position and the midpoint of the triangle in the triangular grid within the reference position data to obtain a plurality of triangular point numbers, and then numbers each triangle in the triangular grid within the reference position data to obtain a triangular serial number, each triangular serial number corresponds to four triangular point numbers;

[0011] Third step: first, generate the distance from each range finder to each triangular point number based on the triangular grid and the coordinate system to obtain the point number distance corresponding to each triangular point number to each range finder, and then collect all the point number distances in the order of the triangular serial number, and then generate a distance document;

[0012] Fourth step: first, place the seismic instrument to be laid out in the target area, and then start the preliminary layout process: first, randomly select a triangular point number that has not been laid out with a seismic instrument, and then sequentially measure the distance from multiple range finders to the same seismic instrument, and then sequentially compare the multiple distances with the corresponding point number distances, when the distance is greater than the corresponding point number distance, move the seismic instrument towards the corresponding range finder to make the difference between the distance and the corresponding point number distance meet the requirements, when the distance is less than the corresponding point number distance, move the seismic instrument away from the corresponding range finder to make the difference between the distance and the corresponding point number distance meet the requirements, when the difference between all distances and corresponding point number distances meets the requirements, it is considered that the seismic instrument is located at the selected triangular point number, that is, a preliminary layout process is completed; select any other triangular point number in the triangular serial number corresponding to the aforementioned triangular point number for the preliminary layout process, until all triangular point numbers in the triangular serial number are laid out with a seismic instrument, and then all seismic instruments simultaneously monitor to obtain multiple micro-motion signals;

[0013] Fifth step: repeat the fourth step until all triangular serial numbers correspond to the micro-motion signals;

[0014] Sixth step: first, process the micro-motion signals to obtain the three-dimensional structure of the shear wave velocity under the target area, and the method ends.

[0015] In the fourth step, the preliminary layout process starts for multiple seismic instruments.

[0016] In the sixth step, the micro-motion signals are processed to obtain the three-dimensional structure of the shear wave velocity under the target area, which is to extract the dispersion curve from the micro-motion signals, and then invert the dispersion curve to obtain the three-dimensional structure of the shear wave velocity under the target area.

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

[0018] In the sixth step, the microseismic signal is processed by a spatial autocorrelation method by first pre-processing the microseismic signal to extract effective microseismic signals from the microseismic signal, and then processing the effective microseismic signals by a spatial autocorrelation method.

[0019] In the sixth step, the dispersion curve is inverted to obtain the three-dimensional structure of the shear wave velocity under the target area by inverting the dispersion curve by a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity under the target area.

[0020] In the sixth step, the dispersion curve is inverted by a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity under the target area by first obtaining the surface wave phase velocity and frequency from the dispersion curve, then inputting the surface wave phase velocity and frequency into the half-wavelength empirical formula to obtain the relationship between depth and surface wave phase velocity, then converting the surface wave phase velocity to shear wave velocity to obtain the relationship between depth and shear wave velocity, and finally obtaining the three-dimensional structure of the shear wave velocity under the target area.

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

[0022] In the first step, the multiple reference points in the range of the target area are selected by selecting more than or equal to three reference points in the range of the target area.

[0023] In the first step, the more than or equal to three reference points in the range of the target area are selected by selecting three reference points in the range of the target area.

[0024] Compared with the prior art, the beneficial effects of the present application are:

[0025] 1. The micro-motion detection method for pre-generating triangular array position of the application, the method comprises the following steps: first step: first determine the target area, then select the reference point, and then place the range finder on the reference point; second step: first measure the boundary points of the target area by the range finder to obtain the reference position data, then generate the coordinate system according to the data, then generate the triangular grid in the coordinate system, then number the positions of each intersection and midpoint of the triangles in the target area to obtain a plurality of triangle point numbers, and then number each triangle to obtain a plurality of triangle serial numbers, each triangle serial number corresponds to four triangle point numbers; third step: first generate the point number distance from each range finder to the triangle point number, and then obtain the distance document; fourth step: first place the seismograph in the target area, and then start the layout process: first select a triangle point number, and then three range finders measure the distances from them to the seismograph, move the seismograph to make the difference between the distance and the point number distance meet the requirements, and the meeting of the requirements means that the seismograph is located at a triangle point number, then repeat the aforementioned layout process for other triangle point numbers corresponding to the triangle serial number of the triangle point number, that is, complete the layout of a group of seismographs, and then obtain the micro-motion signal by monitoring the seismograph; fifth step: repeat the fourth step until the micro-motion signals of all triangle serial numbers are obtained; sixth step: process the micro-motion signals to obtain the three-dimensional structure of the shear wave velocity under the target area, and then judge the three-dimensional structure under the target area according to the three-dimensional structure of the shear wave velocity, and the advantages of the application further include:

[0026] Firstly, the triangular grid is generated, and then the triangle point number, the triangle serial number and the point number distance are obtained, and the position of the triangle point number is the layout position of the seismograph; then the distance between the range finder and the seismograph is measured, when the distance and the point number distance are small, that is, the seismograph is located at the triangle point number, so the seismograph is located at the layout position; the position of the seismograph is determined by multiple range finders, so the position of the seismograph is accurate, the position parameters of the array are accurate, the array position can be associated with the measured micro-motion signal, and the inversion of the underground structure has high precision;

[0027] Secondly, compared with using a measuring rope to layout the seismograph, the layout position of the seismograph of the application is more accurate, and the distance document can be directly generated without subsequent measurement and recording of the position of the seismograph, so the labor and material resources are saved, and the work efficiency is high; and when the seismograph is laid out, only the seismograph needs to be moved, and then the range finder is measured to complete the position correction, so the use is more convenient;

[0028] Thirdly, multiple groups of seismographs can be laid out at the same time, and multiple groups of seismographs can monitor the micro-motion signal at the same time to improve the monitoring efficiency;

[0029] Therefore, the array position parameters of the application are accurate, and the inversion of the underground structure has high precision.

[0030] 2、The microtremor detection method of the application, in the sixth step, the microtremor signal is pretreated, then the spatial autocorrelation method is used to process the microtremor signal to obtain the spatial autocorrelation coefficient, then the spatial autocorrelation coefficient is fitted with the Bessel function, the dispersion curve is obtained, then the surface wave phase velocity and frequency are obtained from the dispersion curve, the relationship between the depth and the shear wave velocity is obtained by the half wavelength empirical formula, that is, the three-dimensional structure of the shear wave velocity under the target area is obtained, that is, the process of inverting the underground structure; the microtremor signals in a target area are dependent on each other, so the spatial autocorrelation method is used to process the microtremor signals, and the spatial autocorrelation method can efficiently process a large amount of microtremor signals, so the underground structure can be efficiently obtained. Therefore, the application has high efficiency.

[0031] 3、The microtremor detection method of the application, in the first step, the number of the reference points is three, the range finders are arranged on the three reference points respectively, then the positions of the seismographs are measured by the three range finders to determine whether the seismographs are located on the triangle points; the number of the range finders can be set according to the needs, but the positions of the seismographs can be accurately determined by the three range finders. Therefore, the application has high positioning accuracy. BRIEF DESCRIPTION OF DRAWINGS

[0032] Figure 1 is a structural schematic diagram of the application.

[0033] Figure 2 is a schematic diagram of the application.

[0034] Figure 3 is Figure 2 a structural schematic diagram of the triangular grid in the application.

[0035] Figure 4 is Figure 2 a structural schematic diagram of the triangle number in the application.

[0036] Figure 5 is a schematic diagram of the distance less than the point number distance in embodiment 1.

[0037] Figure 6 is a schematic diagram of the distance greater than the point number distance in embodiment 1.

[0038] Figure 7 is a schematic diagram of the distance and the point number distance in embodiment 1.

[0039] Figure 8 is Figure 2 a structural schematic diagram of the triangular array in the application.

[0040] Figure 9 is Figure 1 a schematic diagram of the distance document in the application.

[0041] Figure 10 is Figure 9 Enlarged view of the midpoint distance.

[0042] Figure 11 is Figure 1 Schematic view of the three-dimensional structure of shear wave velocity in

[0043] In the figure: target area A, reference point A1, range finder 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

[0044] The application will be further described in detail in conjunction with the accompanying drawings and specific embodiments.

[0045] Please see Figure 1 — Figure 11 A micro-motion detection method for pre-generating triangular array positions, the method comprising the following steps:

[0046] First step: first determine the range of target area A, then select multiple reference points A1 in the range of target area A, and then place a range finder B on each reference point A1 in turn;

[0047] Second step: first make the range finder B measure the boundary points of the target area A, then obtain the reference position data, then set the coordinate system C from the reference position data, the coordinate system C including the origin C1, the horizontal axis C2 and the vertical axis C3, then generate a triangular grid C4 within the range of the origin C1, the horizontal axis C2 and the vertical axis C3, then sequentially number each intersection position of the triangles within the reference position data in the triangular grid C4 and the midpoint position of the triangles to obtain multiple triangular point numbers C5, and then number each triangle within the reference position data in the triangular grid C4 to obtain a triangular sequence number C6, each triangular sequence number C6 corresponding to four triangular point numbers C5;

[0048] Third step: first generate the distance from each range finder B to each triangular point number C5 based on the triangular grid C4 and the coordinate system C to obtain the point number distance corresponding to each triangular point number C5 to each range finder B, then collect all the point number distances in the order of the triangular sequence number C6, and then generate a distance document;

[0049] Fourth step: first place the seismograph D to be laid in the target area A, and then start the preliminary layout process: first randomly select a triangular point number C5 which has not been laid with a seismograph D, and then measure the distance from multiple range finders B 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 towards the corresponding range finder B to make the difference between the distance and the corresponding point number distance meet the requirements, when the distance is less than the corresponding point number distance, move the seismograph D away from the corresponding range finder B to make the difference between the distance and the corresponding point number distance meet the requirements, when the difference between all distances and corresponding point number distances meets the requirements, it is considered that the seismograph D is located on the selected triangular point number C5, that is, one preliminary layout process is completed; then select any other triangular point number C5 in the triangular number C6 corresponding to the aforementioned triangular point number C5 to perform the preliminary layout process, until all the triangular point numbers C5 in the triangular number C6 are laid with seismographs D, and then all the seismographs D simultaneously monitor to obtain multiple micro-motion signals;

[0050] Fifth step: repeat the fourth step until the micro-motion signals corresponding to all triangular numbers C6 are obtained.

[0051] Sixth step: first process the micro-motion signals to obtain the three-dimensional structure of the shear wave velocity under the target area A, and the method ends.

[0052] In the fourth step, the start of the preliminary layout process is to start the preliminary layout process for multiple seismographs D.

[0053] In the sixth step, the processing of the micro-motion signals to obtain the three-dimensional structure of the shear wave velocity under the target area A is to extract the dispersion curve from the micro-motion signals, and then invert the dispersion curve to obtain the three-dimensional structure of the shear wave velocity under the target area A.

[0054] In the sixth step, the extraction of the dispersion curve from the micro-motion signals is to process the micro-motion signals by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, then fit the spatial autocorrelation coefficient with the Bessel function, and then extract the dispersion curve.

[0055] In the sixth step, the processing of the micro-motion signals by the spatial autocorrelation method is to preprocess the micro-motion signals to extract the effective micro-motion signals in the micro-motion signals, and then process the effective micro-motion signals by the spatial autocorrelation method.

[0056] In the sixth step, the inversion of the dispersion curve to obtain the three-dimensional structure of the shear wave velocity under 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 under the target area A.

[0057] In the sixth step, the dispersion curve is inversed by a half-wavelength empirical formula to obtain the three-dimensional structure of shear wave velocity under the target area A.

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

[0059] In the first step, the plurality of reference points A1 selected in the range of the target area A is three or more reference points A1.

[0060] In the first step, the plurality of reference points A1 selected in the range of the target area A is three or more reference points A1.

[0061] The supplementary description of the present application is as follows:

[0062] The measuring rope refers to a measuring tool for measuring length, which is generally a steel wire rope with a scale buckle; one end of the steel wire rope is pulled, and then the length of the pulled steel wire rope is read through the scale buckle, so that the effect of measuring distance is realized; when manually measuring, the error of the measuring rope is large, and the seismograph D needs to be pulled once every time the seismograph D is moved, so that the operation is relatively complicated.

[0063] The range of the target area A is determined according to the actual detection requirement, the boundary point of the target area A is a point having a large turning at the edge of the target area A, and then the reference position data is generated through the position of the boundary point, Figure 2 The virtual boundary line of the target area A is automatically generated according to the reference position data; when the coordinate system C is generated from the reference position data, the reference position data is preferably located in the first quadrant, so that the subsequent steps are facilitated; the triangles outside the reference position data or on the virtual boundary of the reference position data in the triangular mesh C4 are not numbered.

[0064] The plurality of reference points A1 selected in the range of the target area A is that the distance measuring instrument B on the reference point A1 is not physically blocked from all the triangular point numbers C5, so that the laser emitted by the distance measuring instrument B can be reflected by the seismograph D, so that the distance from the distance measuring instrument B to the seismograph D is obtained.

[0065] Embodiment 1:

[0066] Please refer to Figure 1- Figure 11 A micro-motion detection method for pre-generating triangular array positions, the method comprising the following steps:

[0067] First step: first determine the range of target area A, then select multiple reference points A1 in the range of target area A, and then place a range finder B on each reference point A1 in turn;

[0068] Second step: first make the range finder B measure the boundary points of the target area A, then obtain the reference position data, then set a coordinate system C from the reference position data, the coordinate system C comprising an origin C1, a horizontal axis C2 and a vertical axis C3, then generate a triangular grid C4 within the range of the origin C1, the horizontal axis C2 and the vertical axis C3, then sequentially number each intersection position of the triangles within the reference position data in the triangular grid C4 and the midpoint of the triangles to obtain multiple triangular point numbers C5, and then number each triangle within the reference position data in the triangular grid C4 to obtain a triangular serial number C6, each triangular serial number C6 corresponding to four triangular point numbers C5;

[0069] Third step: first generate the distance from each range finder B to each triangular point number C5 based on the triangular grid C4 and the coordinate system C to obtain the point number distance corresponding to each triangular point number C5 to each range finder B, then collect all the point number distances in the order of the triangular serial number C6, and then generate a distance document;

[0070] Fourth step: first place the seismic instrument D to be deployed within the target area A, then start the preliminary deployment process: first randomly select a triangular point number C5 that has not been deployed with a seismic instrument D, then sequentially measure the distance from multiple range finders B to the same seismic instrument D, then sequentially compare the multiple distances with the corresponding point number distances, when the distance is greater than the corresponding point number distance, move the seismic instrument D towards the corresponding range finder B to make the difference between the distance and the corresponding point number distance meet the requirements, when the distance is less than the corresponding point number distance, move the seismic instrument D away from the corresponding range finder B to make the difference between the distance and the corresponding point number distance meet the requirements, when the difference between all distances and corresponding point number distances meets the requirements, consider that the seismic instrument D is located at the selected triangular point number C5, that is, complete one preliminary deployment process; then select any other triangular point number C5 within the triangular serial number C6 corresponding to the aforementioned triangular point number C5 to perform the preliminary deployment process, until all triangular point numbers C5 within the triangular serial number C6 are deployed with a seismic instrument D, then simultaneously monitor all seismic instruments D to obtain multiple micro-motion signals;

[0071] Fifth step: first repeat the fourth step until all the micro-motion signals corresponding to the triangular serial number C6 are obtained;

[0072] The sixth step is to process the micro-motion signal to obtain the three-dimensional structure of the shear wave velocity under the target area A, and the method ends.

[0073] In the fourth step, the initial layout process is started for the plurality of seismic instruments D.

[0074] Embodiment 2:

[0075] The basic content is the same as that in Embodiment 1, except that:

[0076] See Figure 1 Figure 11 In the sixth step, the processing of the micro-motion signal to obtain the three-dimensional structure of the shear wave velocity under the target area A is to extract the dispersion curve from the micro-motion signal, and then to invert the dispersion curve to obtain the three-dimensional structure of the shear wave velocity under the target area A. In the sixth step, the extraction of the dispersion curve from the micro-motion signal is to process the micro-motion signal by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, to fit the spatial autocorrelation coefficient with the Bessel function, and then to extract the dispersion curve. In the sixth step, the processing of the micro-motion signal by the spatial autocorrelation method is to pre-process the micro-motion signal to extract the effective micro-motion signal in the micro-motion signal, and then to process the effective micro-motion 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 under 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 under the target area A. In the sixth step, the inversion of the dispersion curve by the half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity under the target area A is to obtain the surface wave phase velocity and frequency from the dispersion curve, to bring the surface wave phase velocity and frequency into the half-wavelength empirical formula to obtain the relationship between the depth and the surface wave phase velocity, to convert the surface wave phase velocity into the shear wave velocity, and then to obtain the relationship between the depth and the shear wave velocity, i.e. to obtain the three-dimensional structure of the shear wave velocity under the target area A. In the sixth step, the relationship between the depth and the shear wave velocity is to process a plurality of relationships between the depth and the shear wave velocity by the interpolation method, and then to obtain the three-dimensional structure of the shear wave velocity of the target area A.

[0077] ​In application, in the sixth step, the microtremor signal is preprocessed to extract the effective microtremor signal in the microtremor signal, the effective microtremor signal is processed by the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, the spatial autocorrelation coefficient is fitted with the Bessel function, then the dispersion curve is obtained, the surface wave phase velocity and frequency are obtained from the dispersion curve, then the depth and surface wave phase velocity relationship is obtained by putting the surface wave phase velocity and frequency into the half-wavelength empirical formula, the surface wave phase velocity is converted into the shear wave velocity, then the depth and shear wave velocity relationship is obtained, that is, the microtremor signal processing of a group of triangular serial numbers C6 is completed; all the microtremor signals of the triangular serial numbers C6 are processed by the interpolation method, then the shear wave velocity three-dimensional structure of the underground of the target area A is obtained, and the underground structure of the target area A is judged through the shear wave velocity three-dimensional structure of the underground of the target area A.

[0078] Example 3

[0079] The basic content is the same as that in Example 1, except that:

[0080] See Figure 1 — Figure 10 In the first step, the plurality of reference points A1 selected in the range of the target area A are three or more reference points A1 selected in the range of the target area A. In the first step, the three or more reference points A1 selected in the range of the target area A are three reference points A1 selected in the range of the target area A.

[0081] In application, in the first step, three reference points A1 are selected, and a range finder B is placed on each reference point A1; the distance of one seismograph D is measured by the three range finders B, which is sufficient to determine the position of the seismograph D, that is, the seismograph D is located on the triangular point number C5.

[0082] The above only describes the preferred embodiments of the present application, and the protection scope of the present application is not limited to the above embodiments, but any equivalent modification or change made by those skilled in the art according to the disclosed content of the present application shall be included in the protection scope recorded in the claims.

Claims

1. A method for pre-generating the position of a triangular array for micro-motion detection, characterized in that: The method includes 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, use the rangefinder (B) to measure the boundary points of the target area (A) to obtain the reference position data. Then, set the coordinate system (C) based on the reference position data. The coordinate system (C) includes the origin (C1), the horizontal axis (C2), and the vertical axis (C3). Then, generate a triangular grid (C4) within the range of the origin (C1), the horizontal axis (C2), and the vertical axis (C3). Then, number each intersection point and the midpoint of the triangles in the triangular grid (C4) that are located within the reference position data to obtain multiple triangular point numbers (C5). Then, number each triangle in the triangular grid (C4) that is located within the reference position data to obtain a triangular sequence number (C6). Each triangular sequence number (C6) corresponds to four triangular point numbers (C5). Step 3: First, generate the distance from each rangefinder (B) to each triangular point (C5) based on the triangular grid (C4) and coordinate system (C) to obtain the point distance corresponding to each triangular point (C5) to each rangefinder (B). Then, collect all the point distances in the order of triangular number (C6) and generate a distance document. Step 4: First, place the seismograph (D) to be deployed within the target area (A), and then begin the preliminary deployment process: First, arbitrarily select a triangulation point (C5) where no seismograph (D) has been deployed. Then, have multiple distance measuring instruments (B) measure the distance from this point to the same seismograph (D) in sequence. Compare these distances with the corresponding point distances. If the distance is greater than the corresponding point distance, move the seismograph (D) closer to the corresponding distance measuring instrument (B) so that the difference between the distance and the corresponding point distance meets the requirements. If the distance is less than the corresponding point distance, move the seismograph (D) further away from the corresponding point distance. The rangefinder (B) moves in a direction to ensure that the difference between the distance and the corresponding point number meets the requirements. When the difference between all distances and the corresponding point number meets the requirements, the seismograph (D) is considered to be located on the selected triangulation point (C5), thus completing a preliminary deployment process. Then, any other triangulation point (C5) within the triangulation sequence (C6) corresponding to the aforementioned triangulation point (C5) is selected for the preliminary deployment process, until seismographs (D) are deployed on all triangulation points (C5) within the triangulation sequence (C6). All seismographs (D) are then monitored simultaneously to obtain multiple micro-motion signals. Step 5: Repeat step 4 until you obtain the micro-motion signals corresponding to all the triangular numbers (C6); Step 6: First, process the micro-motion signal to obtain the three-dimensional structure of the shear wave velocity below the target area (A). This method is then complete.

2. The micro-motion detection method for pre-generating the position of a triangular array according to claim 1, characterized in that: In the fourth step, the initial deployment process is to begin the initial deployment process for multiple seismographs (D).

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

4. The micro-motion detection method for pre-generating the position of a triangular array according to claim 3, characterized in that: In the sixth step, the extraction of the dispersion curve from the micro-motion signal involves first processing the micro-motion signal using the spatial autocorrelation method to obtain the spatial autocorrelation coefficient, then fitting the spatial autocorrelation coefficient with the Bessel function, and finally extracting the dispersion curve.

5. The micro-motion detection method for pre-generating the position of a triangular array according to claim 4, characterized in that: In the sixth step, the process of processing the micro-motion signal by spatial autocorrelation involves first preprocessing the micro-motion signal to extract the effective micro-motion signal, and then processing the effective micro-motion signal by spatial autocorrelation.

6. The micro-motion detection method for pre-generating the position of a triangular array 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 performed by inverting the dispersion curve using a half-wavelength empirical formula to obtain the three-dimensional structure of the shear wave velocity below the target area (A).

7. The micro-motion detection method for pre-generating the position of a triangular array according to claim 6, characterized in that: In the sixth step, the dispersion curve is inverted using a half-wavelength empirical formula to obtain the three-dimensional structure of shear wave velocity below the target area (A). First, the surface wave phase velocity and frequency are obtained through the dispersion curve. Then, the surface wave phase velocity and frequency are substituted into the half-wavelength empirical formula to obtain the relationship between depth and surface wave phase velocity. Then, the surface wave phase velocity is converted into shear wave velocity, and the relationship between depth and shear wave velocity is obtained, thus obtaining the three-dimensional structure of shear wave velocity below the target area (A).

8. The micro-motion detection method for pre-generating the position of a triangular array according to claim 7, characterized in that: In the sixth step, the relationship between depth and shear wave velocity is obtained by first processing multiple relationships between depth and shear wave velocity using an interpolation method, and then obtaining the three-dimensional structure of the shear wave velocity of the target area (A).

9. A micro-motion detection method for pre-generating the position of a triangular array according to claim 1 or 2, characterized in that: In the first step, selecting multiple reference points (A1) within the target area (A) means selecting three or more reference points (A1) within the target area (A).

10. The micro-motion detection method for 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) means selecting three reference points (A1) within the target area (A).

Images