A signal processing method and detector for tunnel engineering

Abstract

The present invention belongs to the field of tunnel engineering, and relates to a signal processing method and detector for tunnel engineering. In each detector: establish a first data set; cluster the signal amplitude values ​​in the first data set; obtain the numerical interval of each cluster, and assign a priority number to each numerical interval; perform the following steps on each channel respectively: obtain the initial signal amplitude value of the current signal and the signal amplitude value transmitted to the output port; obtain the signal amplitude value of the target signal from other detectors, and establish a second data set; obtain the priority number of each signal amplitude value in the second data set, and establish a third data set; obtain the priority number of the signal amplitude value; determine whether the priority number is the maximum value in the third data set; if so, keep the current channel unblocked; otherwise, close the current channel. The present invention can eliminate low-amplitude signals, ensure that high-amplitude signals enter the analyzer, and improve three-dimensional imaging results.

Description

Technical Field

[0001] The present invention belongs to the field of tunnel engineering, and in particular relates to a signal processing method and a detector for tunnel engineering. Background Art

[0002] In tunnel engineering, an elastic wave advanced geological prediction instrument is used to detect the surrounding rock structure and geological conditions ahead of the tunnel excavation face. It consists of a vibration source, a detector, and an analyzer. The vibration source excites elastic waves near the tunnel sidewall or working face. When the elastic waves encounter interfaces with different wave impedances, complex reflection, refraction, transmission, and scattering occur, causing the elastic waves to attenuate. The reflected wave signals are received by multiple detectors located within the tunnel. The detectors process the collected signals and transmit them to the analyzer, which further processes the signals sent by the detectors to ultimately output a three-dimensional geological structure map. Image analysis of the three-dimensional geological structure map reveals the surrounding rock structure and geological conditions ahead of the tunnel excavation face.

[0003] However, actual 3D geological structure maps often exhibit defects such as discontinuous boundaries and blurred images, making it difficult to accurately analyze the surrounding rock structure and geological conditions ahead of the tunnel excavation face. One of the main reasons for this is that the amplitude of the signal received by the analyzer is significantly reduced, and this low-amplitude signal affects the analyzer's 3D imaging results. Furthermore, the reduced signal amplitude is caused by, on the one hand, the attenuation of elastic waves during propagation through the complex surrounding rock structure, and, on the other hand, interference between the signal channels within the detector. Summary of the Invention

[0004] The technical problem to be solved by the present invention is: how to eliminate low-amplitude signals that affect three-dimensional imaging results before the analyzer receives the signals.

[0005] The present invention is achieved through the following technical solutions:

[0006] In the first aspect, a signal processing method for tunnel engineering is proposed, and the following steps are performed in each detector: a plurality of historically stored signal amplitude values ​​are extracted from a data storage module to establish a first data set; the signal amplitude values ​​in the first data set are clustered to obtain a plurality of clusters; the numerical interval of each cluster is obtained to establish a numerical interval reference list; a corresponding priority number is assigned to each numerical interval according to the median of the numerical interval; the following steps are performed on the x channel, y channel and z channel respectively: S1: obtaining the initial signal amplitude value when the current signal arrives at the sensor, and adding an arrival time mark to the current signal; S2: using the isolation between the current channel and the other two channels and the initial signal amplitude value, obtaining the signal amplitude value of the current signal transmitted to the output port through the current channel; S3: obtaining the signal amplitude value of the current signal transmitted to the output port through the current channel from each other detector Obtain the signal amplitude value of the target signal at the output port of the corresponding channel to establish a second data set; the target signal and the current signal have the same arrival time mark; S4: obtain the priority number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set; S5: obtain the priority number corresponding to the numerical interval to which the signal amplitude value described in S2 belongs; S6: determine whether the priority number described in S5 is the maximum value in the third data set; if so, keep the current channel unblocked; otherwise, close the current channel; before S2, the following steps are also included: obtain the isolation between the x channel and the y channel; obtain the isolation between the x channel and the z channel; obtain the isolation between the y channel and the z channel; the calculation model expression of the isolation is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3; the S2 includes: S21: according to the isolation between the current channel and the other two channels, obtaining the amplitude value of the crosstalk signal caused by the current channel signal in the other two channels; the amplitude value of the crosstalk signal caused by the current channel signal = , A is the initial signal amplitude value; S22: obtain the sum of the initial signal amplitude value and the amplitude values ​​of the crosstalk signals caused by the current channel signals in the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the current channel.

[0007] Technical Effect: This method dynamically selects multiple signal channels (including x-, y-, and z-channels) from multiple geophones within a tunnel through numerical clustering, numerical interval partitioning, isolation analysis, and priority comparison. Only the signal with the highest amplitude at the current moment is transmitted to the analyzer, thereby eliminating low-amplitude signals and preventing them from affecting the analyzer's 3D imaging results. Specifically, the amplitude values ​​of historical signals are clustered. This clustering reveals the distribution characteristics of the amplitude values ​​of these signals, with each cluster corresponding to an amplitude distribution interval. The amplitude values ​​within each distribution interval are then classified and prioritized, with the priority level used as the basis for dynamic channel selection. Furthermore, this method considers the impact of isolation between channels within the geophone on signal transmission and uses the amplitude value of the signal arriving at the channel output port as the criterion for dynamic channel selection. The amplitude value at the output port is the amplitude value of the signal entering the analyzer, and its high or low amplitude value affects the analyzer's 3D imaging results. Finally, the priorities corresponding to the amplitude values ​​of the signals arriving at the output port at the same time are compared. By keeping the channel where the highest priority signal is located unblocked and closing the channels where the remaining lower priority signals are located, low-amplitude signals are eliminated, thereby ensuring that high-amplitude signals enter the analyzer and improving the analyzer's three-dimensional imaging results.

[0008] Furthermore, before S4, the following steps are also included: S3.1: Determine whether each signal amplitude value in the second data set falls within one of the numerical intervals in the numerical interval reference list; if so, execute S4; otherwise, execute S3.2 to S3.5; S3.2: Extract all signal amplitude values ​​that do not fall within any numerical interval in the numerical interval reference list; S3.3: Use all extracted signal amplitude values ​​to establish one or more new numerical intervals to obtain multiple new numerical intervals; S3.4: Add all new numerical intervals to the numerical interval reference list; S3.5: According to the median of the numerical interval, reallocate the corresponding priority number to each numerical interval in the numerical interval reference list.

[0009] Among them, the method for establishing a new numerical interval is: sort all the numerical intervals in the numerical interval reference list in order of priority number from small to large to obtain a numerical interval sequence; traverse the numerical interval sequence; the traversal includes: step A1: determine whether the currently accessed numerical interval is the first numerical interval; if it is the first numerical interval, execute steps B1 to B3; if it is not the first numerical interval, determine whether the currently accessed numerical interval is the last numerical interval; if it is not the last numerical interval, execute steps C1 to C3; if it is the last numerical interval, execute steps D1 to D3; step B1: compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the extracted signal amplitude value is less than the lower limit of the current numerical interval, execute step B2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step B3; step B2: establish a new numerical interval; the upper limit of the new numerical interval is the extracted signal amplitude value; step B3: compare the extracted signal amplitude value with the upper limit of the current numerical interval, access the next numerical interval, and return to the step A1; step C1: compare the extracted signal amplitude value with the upper limit of the current numerical interval, and access the next numerical interval, and return to the step A1; step C1: compare the extracted signal amplitude value with the lower limit of the current numerical interval. The signal amplitude value of the extracted value is compared with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, step C2 is executed; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, step C3 is executed; step C2: establish two new numerical intervals; the lower limit of one new numerical interval is the upper limit of the current numerical interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new numerical interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next numerical interval; step C3: execute the method described in step B3; step Step D1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, execute the method described in step C2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step D2; Step D2: Compare the extracted signal amplitude value with the upper limit of the current numerical interval; if the extracted signal amplitude value is greater than the upper limit of the current numerical interval, execute step D3: Step D3: Establish a new numerical interval; the lower limit of the new numerical interval is the extracted signal amplitude value.

[0010] Technical effect: This method also takes into account the amplitude value differences of the signals collected by each detector. In order to ensure that the amplitude values ​​of the signals collected from other detectors can find the corresponding numerical intervals in this detector, the numerical interval reference list is dynamically updated by interpolation, thereby improving the accuracy of channel screening.

[0011] In the second aspect, a detector for tunnel engineering is proposed, comprising: a data extraction module for extracting a plurality of historically stored signal amplitude values ​​from a data storage module to establish a first data set; a data processing module for clustering the signal amplitude values ​​in the first data set to obtain a plurality of clusters; a list construction module for obtaining the numerical interval of each cluster and establishing a numerical interval reference list; a first serial number assignment module for assigning a corresponding priority serial number to each numerical interval according to the median of the numerical interval; a time stamp module for respectively obtaining the initial signal amplitude value of the current signal in the x channel, y channel and z channel at the moment when the current signal arrives at the sensor, and adding an arrival time stamp to the current signal; an amplitude calculation module for obtaining the signal amplitude value of the current signal transmitted to the output port through the current channel by using the isolation between the current channel and the other two channels and the initial signal amplitude value; a first amplitude acquisition module for obtaining the signal amplitude value of the target signal from the output port of the corresponding channel of each other detector to establish a second data set; the target The reference signal and the current signal have the same arrival time mark; a first sequence number extraction module is used to obtain the priority sequence number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set; a second sequence number extraction module is used to obtain the priority sequence number corresponding to the numerical interval to which the signal amplitude value obtained by the amplitude calculation module belongs; a first analysis control module is used to determine whether the priority sequence number in the second sequence number extraction module is the maximum value in the third data set; if so, control the first channel control module to work; otherwise, control the second channel control module to work; the first channel control module is used to keep the current channel unobstructed; the second channel control module is used to close the current channel; and further includes: a first isolation acquisition module for acquiring the isolation between the x channel and the y channel; a second isolation acquisition module for acquiring the isolation between the x channel and the z channel; a third isolation acquisition module for acquiring the isolation between the y channel and the z channel; the calculation model expression of the isolation is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3; the amplitude calculation module includes: a first amplitude value calculation unit, configured to obtain the amplitude values ​​of the crosstalk signals caused by the signals of the current channel in the other two channels according to the isolation between the current channel and the other two channels; the amplitude value of the crosstalk signals caused by the signals of the current channel = , A is the initial signal amplitude value; a second amplitude value calculation unit is used to obtain the sum of the initial signal amplitude value and the amplitude values ​​of the crosstalk signals caused by the current channel signals in the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the current channel. BRIEF DESCRIPTION OF THE DRAWINGS

[0012] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings:

[0013] Figure 1 This is a flow chart of a signal processing method for tunnel engineering provided in Example 1 of the present invention.

[0014] Figure 2 A schematic diagram of the amplitude distribution characteristics of historical signals provided in Example 1 of the present invention.

[0015] Figure 3 A flowchart of a method for updating a numerical interval reference list provided in Example 1 of the present invention. DETAILED DESCRIPTION

[0016] In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention is further described in detail below in conjunction with the examples. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention. The embodiments described below are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0017] Example 1: In tunnel engineering, when using an elastic wave advanced geological forecaster to detect the surrounding rock structure and geological conditions ahead of the tunnel excavation face, discontinuous boundaries and blurred images are often observed in the three-dimensional geological images generated by the analyzer. This example proposes a signal processing method for tunnel engineering. This method, based on factors influencing the analyzer's three-dimensional imaging performance, uses numerical clustering, numerical interval partitioning, isolation analysis, priority comparison, and dynamic channel selection to shield low-amplitude signals from the analyzer, preventing them from affecting the analyzer's three-dimensional imaging results.

[0018] This method is performed on each detector set in the tunnel. Figure 1 The following steps are shown:

[0019] Step 1: Extract multiple historically stored signal amplitude values ​​from the data storage module to establish a first data set; cluster the signal amplitude values ​​in the first data set to obtain multiple clusters; obtain the numerical interval of each cluster and establish a numerical interval reference list; assign a corresponding priority number to each numerical interval based on the median of the numerical interval.

[0020] The purpose of this method is to filter out the signal with the highest amplitude value from the signal channels (x channel, y channel and z channel) output by each three-component detector and send it to the analyzer, so as to ensure the three-dimensional imaging effect of the analyzer. Specifically, this step uses the priority corresponding to the signal amplitude value at the output port to filter the signal channel in each direction of each three-component detector. In order to obtain the priority, this step clusters the amplitude values ​​of the historical signals to obtain the distribution characteristics of the amplitude values ​​of the historical signals. Each distribution area corresponds to the amplitude value range of the area. Figure 2 As shown in the figure, after clustering the historically stored signal amplitude values, the distribution of the signal amplitude values ​​is obtained; for distribution area 1, the corresponding amplitude value interval is (a1, a2); for distribution area 2, the corresponding amplitude value interval is (a2, a3); for distribution area 3, the corresponding amplitude value interval is (a3, a4). Among them, a1<a2<a3<a4. Combined with Figure 2 A reasonable analysis concludes that if the amplitude of the first signal at the x-channel output port of one detector falls within the amplitude range (a1, a2), and the amplitude of the second signal at the x-channel output port of the other detector falls within the range (a2, a3), then the amplitude of the second signal is definitely greater than that of the first signal, and inputting the second signal into the analyzer is more conducive to 3D imaging. To facilitate priority comparison, this step uses the median of the amplitude range as the corresponding priority for that range, thereby establishing a reference list of numerical ranges (see Table 1).

[0021] Numerical range Priority number (a1, a2) 1 (a2, a3) 2 (a3, a4) 3

[0022] Table 1

[0023] It should be noted that: First, the historically stored signal amplitude value described in this step refers to the amplitude value of the reflected wave signal collected by the detector facing the same surrounding rock mass in front of the tunnel excavation face, so as to ensure that the established amplitude value range is referenceable; second, the larger the priority number, the higher the priority.

[0024] Step 2: Execute S1 to S6 for the x channel, y channel, and z channel respectively.

[0025] A three-component geophone uses three orthogonal sensors to record elastic wave vibration signals in the vertical direction and two horizontal directions, thereby fully capturing the elastic wave's vibration information in three dimensions. Therefore, this method independently performs S1 to S6 for the geophone's x-, y-, and z-channels.

[0026] S1: Get the initial signal amplitude value when the current signal arrives at the sensor and add an arrival time stamp to the current signal.

[0027] The excited elastic waves propagate in the surrounding rock in front of the tunnel excavation face in the form of spherical waves in three-dimensional space. Due to the complex internal structure of the surrounding rock, when the elastic waves encounter interfaces with different wave impedances during propagation, complex reflection, refraction, transmission, and scattering phenomena will occur. Subsequently, the reflected waves, refracted waves, transmitted waves, and scattered waves propagate along different paths and attenuate to varying degrees due to differences in the characteristics of the propagation medium. Therefore, the detector will receive multiple batches of delayed reflected wave signals. After collecting the signal, the detector processes the signal and then sends it to the analyzer. In order to filter out the signal with the highest amplitude value at the output port of each signal channel, it is necessary to unify the objects, that is, to filter from the signals that arrived in the same time batch. This step adds an arrival time mark to the current signal collected by each channel of the detector to provide a reference standard for object unification during subsequent priority comparison.

[0028] Furthermore, elastic waves attenuate during propagation, reducing the signal amplitude. Furthermore, as the collected signal propagates through each channel, interference between channels further reduces the signal amplitude. To obtain the signal amplitude at the output port of each signal channel, this step first obtains the initial signal amplitude at the moment the current signal reaches the sensor.

[0029] S2: Using the isolation between this channel and the other two channels and the initial signal amplitude value, obtain the signal amplitude value of the current signal transmitted to the output port through this channel.

[0030] As explained in S1 above, interference occurs between the x, y, and z channels within the detector, causing the signal amplitude to decrease during transmission. To filter the amplitude of the signal entering the analyzer, it is necessary to first obtain the amplitude of the signal after it has traveled through the channels and reaches the output port.

[0031] Channel isolation refers to the degree of signal interference between different channels in a multichannel system. Higher channel isolation means less crosstalk between channels, resulting in higher accuracy and reliability in signal acquisition. This step uses channel isolation analysis to calculate the amplitude of the signal when it reaches the output port after traveling through the channels. Prior to this, it's necessary to obtain the isolation between channels x and y, between channels x and z, and between channels y and z.

[0032] The calculation model expression of isolation is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3. It should be noted that the peak-to-peak value describes the total range of variation of a signal from its maximum value to its minimum value within a certain period of time. It is a key indicator for measuring the dynamic range and fluctuation amplitude of a signal.

[0033] S2 includes:

[0034] S2.1: Based on the isolation between the current channel and the other two channels, obtain the amplitude values ​​of the crosstalk signals in the other two channels caused by the current channel signal; the amplitude value of the crosstalk signal caused by the current channel signal = , A is the initial signal amplitude value;

[0035] S2.2: Obtain the sum of the amplitude value of the initial signal and the amplitude values ​​of the crosstalk signals caused by the signals of the current channel in the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the current channel.

[0036] Taking channel x as an example, the amplitude value of the signal arriving at the channel output port is calculated as follows:

[0037] Assume that the initial signal amplitude of the x channel is A, and the isolation between the x channel and the y channel is d 1. The isolation between x channel and z channel is d 2, then the amplitude of the crosstalk signal in the y channel caused by the x channel signal is The amplitude of the crosstalk signal in the z channel caused by the x channel signal is In the x channel, in addition to the original signal, there will be interference signals from the y channel and the z channel. Therefore, the amplitude value of the signal after it is transmitted through the x channel and reaches the output port is A total for .

[0038] S3: Acquire the signal amplitude value of the target signal from the output port of the corresponding channel of each other detector to establish a second data set.

[0039] Taking the x channel as an example, after obtaining the signal amplitude value at the x channel output port of one detector, in order to perform signal screening, it is necessary to collect the corresponding signal amplitude value from the x channel output port of each other detector.

[0040] S4: Store the signal amplitude value at the output port into the data storage module.

[0041] The purpose of this step is to provide data samples for the next processing according to step 1 of this method.

[0042] S5: Obtain the priority number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set.

[0043] The second data set stores the signal amplitude value of the target signal obtained at the output port of the corresponding channel of each detector. In order to filter out the signal with the highest amplitude value through priority comparison, it is necessary to obtain the priority number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs.

[0044] It is worth noting that, since the amplitude value of the current signal collected by the detector is random, and the amplitude value intervals established in step 1 are based on the historical signal amplitude values ​​of the detector, it is possible that the amplitude value of the current signal collected by other detectors is not within the range of the various numerical intervals established in step 1. In order to make the amplitude values ​​of the signals collected from other detectors comparable with the amplitude value of the signal of the detector, and to avoid such a situation, this embodiment adopts the following measures:

[0045] Determine whether each signal amplitude value in the second data set falls within one of the numerical intervals in the numerical interval reference list. If the signal amplitude value in the second data set falls within one of the numerical intervals in the numerical interval reference list, the subsequent steps can be performed smoothly; if the signal amplitude value in the second data set does not fall within any of the numerical intervals in the numerical interval reference list, it is necessary to fill in the missing numerical intervals in the numerical interval reference list. Specifically including Figure 3 The following steps are shown:

[0046] Step Q1: Extract all signal amplitude values ​​that do not fall within any value interval in the value interval reference list.

[0047] For example, three signal amplitude values ​​are extracted that do not fall within any numerical interval in the numerical interval reference list, namely a2, a5, and a8.

[0048] Step Q2: Establish one or more new numerical intervals using all extracted signal amplitude values ​​to obtain multiple new numerical intervals.

[0049] The method for establishing a new numerical interval includes the following steps:

[0050] Step Q2.1: Sort all the numerical intervals in the numerical interval reference list in ascending order of priority numbers to obtain a numerical interval sequence.

[0051] For example, sort the three numerical intervals in Table 1 to obtain [(a1, a2), (a2, a3), (a3, a4)].

[0052] Step Q2.2: Traverse the sequence of numerical intervals.

[0053] The traversal includes:

[0054] Step A1: Determine whether the currently accessed numerical interval is the first numerical interval; if it is the first numerical interval, execute steps B1 to B3; if it is not the first numerical interval, determine whether the currently accessed numerical interval is the last numerical interval; if it is not the last numerical interval, execute steps C1 to C3; if it is the last numerical interval, execute steps D1 to D3.

[0055] Step B1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the extracted signal amplitude value is less than the lower limit of the current numerical interval, execute step B2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step B3.

[0056] Step B2: Establish a new numerical interval; the upper limit of the new numerical interval is the extracted signal amplitude value.

[0057] Step B3: Compare the extracted signal amplitude value with the upper limit of the current numerical interval, access the next numerical interval, and return to the step A1.

[0058] Step C1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, execute step C2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step C3.

[0059] Step C2: Establish two new numerical intervals; the lower limit of one new numerical interval is the upper limit of the current numerical interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new numerical interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next numerical interval.

[0060] Step C3: Execute the method described in step B3.

[0061] Step D1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, execute the method described in step C2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step D2.

[0062] Step D2: Compare the extracted signal amplitude value with the upper limit of the current numerical interval; if the extracted signal amplitude value is greater than the upper limit of the current numerical interval, execute step D3.

[0063] Step D3: Establish a new numerical interval; the lower limit of the new numerical interval is the extracted signal amplitude value.

[0064] Taking the three amplitude values ​​b2, b5, and b8 extracted above as an example, after processing in step Q2, the following six new value intervals are established: (a1, b2), (b2, a2), (a2, b5), (b5, a3), (a3, b8), and (b8, a4). Among them, a1 < b2 < a2 < b5 < a3 < b8 < a4.

[0065] Step Q3: Add all new value intervals to the value interval reference list.

[0066] Step Q4: reallocate a corresponding priority number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

[0067] Add the above 6 new numerical ranges to Table 1, and update Table 1 to obtain Table 2.

[0068]

[0069] Table 2

[0070] S6: Obtain the priority number corresponding to the numerical interval to which the signal amplitude value described in S2 belongs.

[0071] Similarly, before this step, it is also necessary to determine whether the priority number described in S2 falls within one of the numerical intervals in the numerical interval reference list (Table 2). If so, execute S7; otherwise, update Table 2 according to the method provided in steps Q1 to Q4.

[0072] For example, if the signal amplitude value of S2 is a6, then the numerical interval to which it belongs is (b5, a3), and the corresponding priority number is 4.

[0073] S7: Determine whether the priority number described in S6 is the maximum value in the third data set; if so, keep the channel open; otherwise, close the channel.

[0074] In summary, the signal processing method for tunnel engineering proposed in this embodiment first clusters the amplitude values ​​of collected historical signals. This clustering process derives the distribution characteristics of the amplitude values ​​of the historical signals, i.e., each cluster corresponds to an amplitude distribution interval. The amplitude values ​​within each distribution interval can then be classified based on their numerical values, and the priority level can be used as a basis for dynamically selecting signal channels. Furthermore, this method considers the impact of isolation between channels within the detector on signal transmission and uses the amplitude value of the signal arriving at the channel output port as the criterion for dynamically selecting signal channels. The amplitude value at the output port is the amplitude value of the signal entering the analyzer, and the amplitude value will affect the analyzer's three-dimensional imaging results. Finally, the priorities corresponding to the amplitude values ​​of the signals arriving at the output port at the same time are compared. By keeping the channel containing the highest-priority signal open and closing the channels containing the remaining lower-priority signals, low-amplitude signals are eliminated, thereby ensuring that high-amplitude signals enter the analyzer and improving the analyzer's three-dimensional imaging results.

[0075] Example 2: Corresponding to the above-mentioned Example 1, this example proposes a detector for tunnel engineering, including:

[0076] A data extraction module, configured to extract a plurality of historically stored signal amplitude values ​​from the data storage module to establish a first data set;

[0077] a data processing module, configured to cluster the signal amplitude values ​​in the first data set to obtain a plurality of clusters;

[0078] List building module, used to obtain the numerical interval of each cluster and build a numerical interval reference list;

[0079] A first sequence number allocation module is used to allocate a corresponding priority sequence number to each numerical interval according to the median of the numerical interval;

[0080] A time stamp module is used to obtain the initial signal amplitude value of the current signal in the x channel, y channel and z channel when it arrives at the sensor, and add an arrival time stamp to the current signal;

[0081] an amplitude calculation module, configured to obtain a signal amplitude value of the current signal transmitted to the output port through the current channel by using the isolation between the current channel and the other two channels and the initial signal amplitude value;

[0082] The amplitude calculation module includes:

[0083] The first amplitude value calculation unit is configured to obtain the amplitude values ​​of the crosstalk signals in the other two channels caused by the signals of the current channel according to the isolation between the current channel and the other two channels; the amplitude value of the crosstalk signals caused by the signals of the current channel = , A is the initial signal amplitude value;

[0084] a second amplitude value calculation unit, configured to obtain the sum of the amplitude value of the initial signal and the amplitude values ​​of the crosstalk signals caused by the signals of the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the channel;

[0085] A first amplitude acquisition module is configured to acquire a signal amplitude value of a target signal from an output port of a corresponding channel of each other detector to establish a second data set; the target signal and the current signal have the same arrival time stamp;

[0086] A first serial number extraction module is used to obtain the priority serial number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set;

[0087] A second serial number extraction module is used to obtain the priority serial number corresponding to the numerical interval to which the signal amplitude value obtained by the amplitude calculation module belongs;

[0088] a first analysis control module, configured to determine whether the priority number in the second sequence number extraction module is the maximum value in the third data set; if so, control the first channel control module to operate; otherwise, control the second channel control module to operate;

[0089] A first channel control module, configured to keep the channel unobstructed;

[0090] A second channel control module, configured to close the channel;

[0091] A first isolation acquisition module, configured to acquire the isolation between the x channel and the y channel;

[0092] A second isolation acquisition module, configured to acquire the isolation between the x channel and the z channel;

[0093] A third isolation acquisition module is used to acquire the isolation between the y channel and the z channel;

[0094] The calculation model expression of the isolation degree is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3;

[0095] A data transfer module, configured to store the signal amplitude value at the output port into the data storage module;

[0096] A data storage module, used for storing the data sent by the data transfer module;

[0097] a second analysis control module, configured to determine whether each signal amplitude value in the second data set falls within a numerical interval in the numerical interval reference list; if so, control the operation of the first sequence number extraction module; otherwise, control the operation of the second amplitude acquisition module, the interval construction module, the list updating module, or the second sequence number allocation module;

[0098] A second amplitude acquisition module is used to extract all signal amplitude values ​​that do not fall within any numerical interval of the numerical interval reference list;

[0099] An interval construction module is used to establish one or more new numerical intervals using all the extracted signal amplitude values ​​to obtain multiple new numerical intervals;

[0100] a list updating module, configured to add all new value intervals to the value interval reference list;

[0101] The second sequence number allocation module is used to reallocate a corresponding priority sequence number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

[0102] The interval building modules include:

[0103] An interval sorting unit is used to sort all the numerical intervals in the numerical interval reference list in ascending order of priority numbers to obtain a numerical interval sequence;

[0104] An interval traversal unit, configured to traverse the sequence of numerical intervals and control the operation of the first analysis control unit each time a numerical interval is accessed;

[0105] a first analysis control unit, configured to determine whether the currently accessed numerical interval is the first numerical interval; if so, to control the first numerical comparison unit to operate; if not, to determine whether the currently accessed numerical interval is the last numerical interval; if not, to control the second numerical comparison unit to operate; if so, to control the third numerical comparison unit to operate;

[0106] A first numerical comparison unit is configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the extracted signal amplitude value is less than the lower limit of the current numerical interval, the first numerical interval construction unit is controlled to operate; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, the second analysis control unit is controlled to operate;

[0107] A first numerical interval construction unit is used to establish a new numerical interval; the upper limit of the new numerical interval is the extracted signal amplitude value;

[0108] a second analysis control unit, configured to compare the extracted signal amplitude value with the upper limit of the current numerical interval, control the interval traversal unit to access the next numerical interval, and control the operation of the first analysis control unit;

[0109] a second numerical comparison unit, configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, then controlling the second numerical interval construction unit to operate; and if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, then controlling the third analysis control unit to operate;

[0110] The second numerical interval construction unit is used to establish two new numerical intervals; the lower limit of one new numerical interval is the upper limit of the current numerical interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new numerical interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next numerical interval;

[0111] a third analysis control unit, configured to drive the second analysis control unit to operate;

[0112] a third numerical comparison unit, configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, then controlling the second numerical interval construction unit to operate; and if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, then controlling the fourth numerical comparison unit to operate;

[0113] a fourth numerical comparison unit, configured to compare the extracted signal amplitude value with the upper limit of the current numerical interval; if the extracted signal amplitude value is greater than the upper limit of the current numerical interval, controlling the third numerical interval construction unit to operate;

[0114] The third numerical interval construction unit establishes a new numerical interval; the lower limit of the new numerical interval is the extracted signal amplitude value.

[0115] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

[0116] It should be noted that the structures, proportions, sizes, etc. illustrated in the drawings of this specification are only used to match the contents disclosed in the specification for understanding and reading by those familiar with this technology, and are not used to limit the conditions for implementation of the present invention. Therefore, they have no substantial technical significance. Any modification of the structure, change in the proportional relationship, or adjustment of the size should still fall within the scope of the technical content disclosed by the present invention without affecting the efficacy and purpose of the present invention. At the same time, the terms such as "upper", "lower", "left", "right", "middle", etc. quoted in this specification are only for the convenience of description and are not used to limit the scope of implementation of the present invention. Changes or adjustments in their relative relationships should also be regarded as the scope of implementation of the present invention without substantially changing the technical content.

Claims

1. A signal processing method for tunnel engineering, characterized in that: Perform the following steps in each detector: Extracting a plurality of historically stored signal amplitude values ​​from a data storage module to establish a first data set; clustering the signal amplitude values ​​in the first data set to obtain a plurality of clusters; obtaining a numerical interval for each cluster and establishing a numerical interval reference list; and assigning a corresponding priority number to each numerical interval based on the median of the numerical interval; Perform the following steps for the x-channel, y-channel, and z-channel respectively: S1: Obtain the initial signal amplitude value when the current signal arrives at the sensor, and add an arrival time stamp to the current signal; S2: using the isolation between the current channel and the other two channels and the initial signal amplitude value, obtaining the signal amplitude value of the current signal transmitted to the output port through the current channel; S3: Acquire the signal amplitude value of the target signal from the output port of the corresponding channel of each other detector to establish a second data set; the target signal and the current signal have the same arrival time stamp; S4: Obtain the priority number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set; S5: Obtain the priority number corresponding to the numerical interval to which the signal amplitude value described in S2 belongs; S6: Determine whether the priority number in S5 is the maximum value in the third data set; if so, keep the current channel open; otherwise, close the current channel; Before S2, the steps further include: obtaining the isolation between the x channel and the y channel; obtaining the isolation between the x channel and the z channel; obtaining the isolation between the y channel and the z channel; The calculation model expression of the isolation degree is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3; The S2 includes: S21: Obtaining the amplitude values ​​of the crosstalk signals in the other two channels caused by the signals of the current channel according to the isolation between the current channel and the other two channels; the amplitude value of the crosstalk signals caused by the signals of the current channel = , A is the initial signal amplitude value; S22: Acquire the sum of the amplitude value of the initial signal and the amplitude values ​​of the crosstalk signals caused by the signals of the current channel in the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the current channel.

2. A signal processing method for tunnel engineering according to claim 1, characterized in that: Between S2 and S3, the following step is also included: storing the signal amplitude value at the output port into the data storage module.

3. The signal processing method for tunnel engineering according to claim 1, characterized in that: Before S4, the following steps are also included: S3.1: Determine whether each signal amplitude value in the second data set falls within a numerical interval in the numerical interval reference list; if so, execute S4; otherwise, execute S3.2 to S3.5; S3.2: Extract all signal amplitude values ​​that do not fall within any numerical interval in the numerical interval reference list; S3.3: Establish one or more new numerical intervals using all extracted signal amplitude values ​​to obtain multiple new numerical intervals; S3.4: Add all new value intervals to the value interval reference list; S3.5: Reassign a corresponding priority number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

4. A signal processing method for tunnel engineering according to claim 3, characterized in that: The method to create a new value interval is: Sort all the numerical intervals in the numerical interval reference list in ascending order of priority numbers to obtain a numerical interval sequence; Traversing the sequence of numerical intervals; The traversal includes: Step A1: Determine whether the currently accessed numerical interval is the first numerical interval; if it is the first numerical interval, execute steps B1 to B3; if it is not the first numerical interval, determine whether the currently accessed numerical interval is the last numerical interval; if it is not the last numerical interval, execute steps C1 to C3; if it is the last numerical interval, execute steps D1 to D3; Step B1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the extracted signal amplitude value is less than the lower limit of the current numerical interval, execute step B2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step B3; Step B2: Establish a new numerical interval; the upper limit of the new numerical interval is the extracted signal amplitude value; Step B3: comparing the extracted signal amplitude value with the upper limit of the current numerical interval, accessing the next numerical interval, and returning to step A1; Step C1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, execute step C2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step C3; Step C2: establishing two new numerical intervals; the lower limit of one new numerical interval is the upper limit of the current numerical interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new numerical interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next numerical interval; Step C3: executing the method described in step B3; Step D1: Compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, execute the method described in step C2; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, execute step D2; Step D2: Compare the extracted signal amplitude value with the upper limit of the current numerical interval; if the extracted signal amplitude value is greater than the upper limit of the current numerical interval, execute step D3; Step D3: Establish a new numerical interval; the lower limit of the new numerical interval is the extracted signal amplitude value.

5. The signal processing method for tunnel engineering according to claim 4, characterized in that: After S5, the following steps are also included: Determine whether the priority number described in S2 falls within one of the numerical intervals in the numerical interval reference list; if not, execute steps S3.1 to S3.

5.

6. A detector for tunnel engineering, characterized in that: include: A data extraction module, configured to extract a plurality of historically stored signal amplitude values ​​from the data storage module to establish a first data set; a data processing module, configured to cluster the signal amplitude values ​​in the first data set to obtain a plurality of clusters; List building module, used to obtain the numerical interval of each cluster and build a numerical interval reference list; A first sequence number allocation module is used to allocate a corresponding priority sequence number to each numerical interval according to the median of the numerical interval; A time stamp module is used to obtain the initial signal amplitude value of the current signal in the x channel, y channel and z channel when it arrives at the sensor, and add an arrival time stamp to the current signal; an amplitude calculation module, configured to obtain a signal amplitude value of the current signal transmitted to the output port through the current channel by using the isolation between the current channel and the other two channels and the initial signal amplitude value; A first amplitude acquisition module is configured to acquire a signal amplitude value of a target signal from an output port of a corresponding channel of each other detector to establish a second data set; the target signal and the current signal have the same arrival time stamp; A first serial number extraction module is used to obtain the priority serial number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs, and establish a third data set; A second serial number extraction module is used to obtain the priority serial number corresponding to the numerical interval to which the signal amplitude value obtained by the amplitude calculation module belongs; a first analysis control module, configured to determine whether the priority number in the second sequence number extraction module is the maximum value in the third data set; if so, control the first channel control module to operate; otherwise, control the second channel control module to operate; A first channel control module, configured to keep the channel unobstructed; A second channel control module, configured to close the channel; Also includes: A first isolation acquisition module, configured to acquire the isolation between the x channel and the y channel; A second isolation acquisition module, configured to acquire the isolation between the x channel and the z channel; A third isolation acquisition module is used to acquire the isolation between the y channel and the z channel; The calculation model expression of the isolation degree is: ;in, Indicates the i Channel and j Isolation between channels, Indicates the j The average peak-to-peak value of the signal amplitude of each channel, Indicates the i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3; The amplitude calculation module includes: a first amplitude value calculation unit, which is used to obtain the amplitude value of the crosstalk signal caused by the current channel signal in the other two channels according to the isolation between the current channel and the other two channels; the amplitude value of the crosstalk signal caused by the current channel signal = , A is the initial signal amplitude value; The second amplitude value calculation unit is used to obtain the sum of the initial signal amplitude value and the amplitude values ​​of the crosstalk signals caused by the local channel signals in the other two channels, and obtain the signal amplitude value of the current signal transmitted to the output port through the local channel.

7. The detector for tunnel engineering according to claim 6, characterized in that: Also includes: A data transfer module, configured to store the signal amplitude value at the output port into the data storage module; The data storage module is used to store the data sent by the data transfer module.

8. The detector for tunnel engineering according to claim 6, characterized in that: Also includes: a second analysis and control module, configured to determine whether each signal amplitude value in the second data set falls within one of the numerical intervals in the numerical interval reference list; If yes, control the first serial number extraction module to work; otherwise, control the second amplitude acquisition module, the interval construction module, the list update module or the second serial number allocation module to work; A second amplitude acquisition module is used to extract all signal amplitude values ​​that do not fall within any numerical interval of the numerical interval reference list; An interval construction module is used to establish one or more new numerical intervals using all the extracted signal amplitude values ​​to obtain multiple new numerical intervals; a list updating module, configured to add all new value intervals to the value interval reference list; The second sequence number allocation module is used to reallocate a corresponding priority sequence number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

9. The geophone for tunnel engineering according to claim 8, characterized in that: The interval construction module includes: An interval sorting unit is used to sort all the numerical intervals in the numerical interval reference list in ascending order of priority numbers to obtain a numerical interval sequence; An interval traversal unit, configured to traverse the sequence of numerical intervals and control the operation of the first analysis control unit each time a numerical interval is accessed; a first analysis control unit, configured to determine whether the currently accessed numerical interval is the first numerical interval; if so, to control the first numerical comparison unit to operate; if not, to determine whether the currently accessed numerical interval is the last numerical interval; if not, to control the second numerical comparison unit to operate; if so, to control the third numerical comparison unit to operate; A first numerical comparison unit is configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the extracted signal amplitude value is less than the lower limit of the current numerical interval, the first numerical interval construction unit is controlled to operate; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, the second analysis control unit is controlled to operate; A first numerical interval construction unit is used to establish a new numerical interval; the upper limit of the new numerical interval is the extracted signal amplitude value; a second analysis control unit, configured to compare the extracted signal amplitude value with the upper limit of the current numerical interval, control the interval traversal unit to access the next numerical interval, and control the operation of the first analysis control unit; a second numerical comparison unit, configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, then controlling the second numerical interval construction unit to operate; and if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, then controlling the third analysis control unit to operate; The second numerical interval construction unit is used to establish two new numerical intervals; the lower limit of one new numerical interval is the upper limit of the current numerical interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new numerical interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next numerical interval; a third analysis control unit, configured to drive the second analysis control unit to operate; a third numerical comparison unit, configured to compare the extracted signal amplitude value with the lower limit of the current numerical interval; if the upper limit of the previous numerical interval is less than the extracted signal amplitude value and less than the lower limit of the current numerical interval, then controlling the second numerical interval construction unit to operate; and if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, then controlling the fourth numerical comparison unit to operate; a fourth numerical comparison unit, configured to compare the extracted signal amplitude value with the upper limit of the current numerical interval; if the extracted signal amplitude value is greater than the upper limit of the current numerical interval, controlling the third numerical interval construction unit to operate; The third numerical interval construction unit establishes a new numerical interval; the lower limit of the new numerical interval is the extracted signal amplitude value.

Images