Signal processing method and detector for tunnel engineering

Abstract

The invention belongs to the field of tunnel engineering, and relates to a signal processing method and a detector for tunnel engineering. In each detector: establishing a first data set; clustering the signal amplitude values in the first data set; obtaining a numerical value interval of each cluster, and distributing a priority sequence number for each numerical value interval; respectively executing the following steps on each channel: acquiring an initial signal amplitude value of a current signal and a signal amplitude value transmitted to the output port; acquiring a signal amplitude value of the target signal from other detectors, and establishing a second data set; acquiring a priority sequence number of each signal amplitude value in the second data set, and establishing a third data set; acquiring a priority sequence number of the signal amplitude value; judging whether the priority serial number is the maximum value in the third data set or not; if yes, the channel is kept unblocked; otherwise, closing the channel. According to the invention, low-amplitude signals can be eliminated, high-amplitude signals are ensured to enter an analyzer, and the three-dimensional imaging result is improved.

Description

Technical Field

[0001] The 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, the elastic wave advanced geological prediction instrument is used to detect the surrounding rock structure and geological conditions in front of the tunnel excavation face. The elastic wave advanced geological prediction instrument consists of a vibration source, a detector and an analyzer. Among them, the vibration source excites elastic waves near the tunnel side wall or the face. When the elastic wave encounters an interface with a wave impedance difference, complex reflection, refraction, transmission and scattering phenomena will occur, causing the elastic wave to attenuate its energy; the reflected wave signal is received by multiple detectors arranged in the tunnel, and the detector processes the collected signal and sends the processed signal to the analyzer; the analyzer further processes the signal sent by the detector and finally outputs a three-dimensional geological structure map. The surrounding rock structure and geological conditions in front of the tunnel excavation face are obtained by image analysis of the three-dimensional geological structure map.

[0003] However, in the actual three-dimensional geological structure map, there are often defects such as discontinuous boundaries and blurred images, which makes it difficult to accurately analyze the surrounding rock structure and geological conditions in front of the tunnel excavation face. One of the main reasons is that the amplitude of the signal received by the analyzer is significantly reduced, and these low-amplitude signals will affect the three-dimensional imaging results of the analyzer. Furthermore, the reasons for the reduction in signal amplitude are: on the one hand, elastic waves are attenuated during the propagation process of complex surrounding rock structures; on the other hand, the signal channels in the detector interfere with each other. Summary of the invention

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

[0005] The present invention is achieved through the following technical solutions: In the first aspect, a signal processing method for tunnel engineering is proposed, and the following steps are performed in each detector: extract multiple historically stored signal amplitude values ​​from a 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 according to the median of the numerical interval; perform the following steps on 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 mark to the current signal; S2: use the isolation between the current channel and the other two channels and the initial signal amplitude to obtain the initial signal amplitude value; degree value, obtain the signal amplitude value of the current signal transmitted to the output port through the current channel; S3: obtain the signal amplitude value of the target signal from the output port of the corresponding channel of each other detector, and 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.

[0006] Technical effect: This method realizes the dynamic selection of multiple signal channels (including x channel, y channel and z channel) of multiple detectors in the tunnel through numerical clustering, numerical interval division, isolation analysis and priority comparison, and only allows the signal with the highest amplitude value at the current moment to be sent to the analyzer, thereby eliminating the low-amplitude signal and avoiding the low-amplitude signal from affecting the three-dimensional imaging results of the analyzer. Specifically, the amplitude values ​​of the collected historical signals are clustered first, and the distribution characteristics of the amplitude values ​​of the historical signals can be obtained through clustering, that is, each cluster corresponds to an amplitude value distribution interval, and then the amplitude values ​​in each distribution interval can be graded, and the priority level can be used as the basis for dynamically selecting the signal channel. Furthermore, this method takes into account the influence of the isolation between the channels inside the detector on the signal transmission, and uses the amplitude value of the signal arriving at the output port of the channel as the standard for dynamically selecting the signal channel; the amplitude value at the output port is the amplitude value of the signal entering the analyzer, and the amplitude value will affect the three-dimensional imaging results of the analyzer. 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, the low-amplitude signal is eliminated, thereby ensuring that the high-amplitude signal enters the analyzer and improving the three-dimensional imaging results of the analyzer.

[0007] 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 the 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 for each numerical interval in the numerical interval reference list.

[0008] 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, and 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 value interval; if the upper limit of the previous value interval is less than the extracted signal amplitude value and less than the lower limit of the current value interval, step C2 is executed; if the extracted signal amplitude value is greater than the lower limit of the current value interval, step C3 is executed; step C2: establish two new value intervals; the lower limit of one new value interval is the upper limit of the current value interval, and the upper limit is the extracted signal amplitude value; the lower limit of the other new value interval is the extracted signal amplitude value, and the upper limit is the lower limit of the next value 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.

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

[0010] In the second aspect, a detector for tunnel engineering is proposed, comprising: a data extraction module, used to extract multiple historically stored signal amplitude values ​​from a data storage module to establish a first data set; a data processing module, used to cluster the signal amplitude values ​​in the first data set to obtain multiple clusters; a list construction module, used to obtain the numerical interval of each cluster and establish a reference list of numerical intervals; a first sequence number allocation module, used to allocate a corresponding priority sequence number to each numerical interval according to the median of the numerical interval; a time stamp module, used to obtain the initial signal amplitude value when the current signal arrives at the sensor, and add an arrival time stamp to the current signal; an amplitude calculation module, used to use the isolation between the current channel and the other two channels and the initial signal amplitude value to obtain the signal amplitude value of the current signal transmitted to the output port through the current channel; a first amplitude sampling module, used to obtain 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 collection module is used to obtain 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 mark; a first sequence number extraction module is used to obtain the priority number corresponding to the numerical interval to which each signal amplitude value in the second data set belongs to establish a third data set, and a second sequence number extraction module is used to obtain the priority 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 number in the second sequence number extraction module is the maximum value in the third data set; if so, the first channel control module is controlled to work; otherwise, the second channel control module is controlled 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. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] 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: Figure 1 A flow chart of a signal processing method for tunnel engineering provided in Example 1 of the present invention.

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

[0013] Figure 3 A schematic flow chart of a method for updating a reference list of numerical intervals provided in Embodiment 1 of the present invention. DETAILED DESCRIPTION

[0014] In order to make the purpose, technical scheme and advantages of the present invention clearer, the present invention is further described in detail below in conjunction with the embodiments. The schematic implementation mode of the present invention and its 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, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in the field without creative work are within the scope of protection of the present invention.

[0015] Embodiment 1: In view of the problem of discontinuous boundaries and blurred images in the three-dimensional geological image generated by the analyzer when the elastic wave advanced geological prediction instrument is used to detect the surrounding rock structure and geological conditions in front of the tunnel excavation face in tunnel engineering, this embodiment proposes a signal processing method for tunnel engineering. Starting from the factors affecting the three-dimensional imaging effect of the analyzer, this method shields low-amplitude signals outside the analyzer through means such as numerical clustering, numerical interval division, isolation analysis, priority comparison and dynamic channel selection, so as to avoid the low-amplitude signal from affecting the three-dimensional imaging results of the analyzer.

[0016] This method is performed for each detector set in the tunnel. Figure 1 The following steps are shown: 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.

[0017] The purpose of this method is to select the signal with the highest amplitude value from the signal output by the signal channel (x channel, y channel and z channel) of 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 size 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 interval 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, a4); for distribution area 2, the corresponding amplitude value interval is (a4, a7); for distribution area 3, the corresponding amplitude value interval is (a7, a9). Among them, a1<a4<a7<a9. Combined with Figure 2It can be reasonably analyzed that if the amplitude value of the first signal at the x-channel output port of one detector falls within the amplitude value interval (a1, a4), and the amplitude value of the second signal at the x-channel output port of the other detector falls within the interval (a4, a7), then the amplitude value of the second signal must be greater than the amplitude value of the first signal, and it is more conducive to three-dimensional imaging to select the second signal to be input into the analyzer. In order to facilitate priority comparison, this step uses the median of the amplitude value interval as the priority corresponding to the amplitude value interval, thereby establishing a reference list of numerical intervals, as shown in Table 1.

[0018] Amplitude value range Priority number (a1, a4) 1 (a4, a7) 2 (a7, a9) 3

[0019] Table 1 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.

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

[0021] The three-component geophone records the vibration signals of the elastic wave in the vertical direction and two horizontal directions through three orthogonal sensors, thereby completely capturing the vibration information of the elastic wave in three-dimensional space. Therefore, the method independently performs S1 to S6 on the x channel, y channel and z channel of the geophone.

[0022] 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.

[0023] 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, complex reflection, refraction, transmission and scattering phenomena will occur when the elastic waves encounter the interface with different wave impedance during propagation. Then the reflected waves, refracted waves, transmitted waves and scattered waves propagate along different paths and attenuate to varying degrees due to the 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 arriving in the same time batch. This step adopts the method of adding the 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.

[0024] In addition, elastic waves will attenuate during propagation, resulting in a decrease in the signal amplitude value; further, during the transmission of the collected signal in each channel, the signal amplitude value is further reduced due to mutual interference between channels. In order to obtain the signal amplitude value at the output port of each signal channel, this step first obtains the initial signal amplitude value when the current signal reaches the sensor.

[0025] 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.

[0026] According to the explanation of S1 above, there is mutual interference between the x channel, y channel and z channel in the detector, which causes the amplitude value of the signal to decrease during the transmission process. In order to screen the amplitude value of the signal entering the analyzer, it is necessary to first obtain the amplitude value of the signal when it reaches the output port after passing through the channel transmission.

[0027] In addition, channel isolation refers to the degree of mutual interference between signals between different channels in a multi-channel system. The higher the channel isolation, the smaller the crosstalk between channels, and the higher the accuracy and reliability of signal acquisition. This step calculates the amplitude value of the signal when it reaches the output port after passing through the channel through channel isolation analysis. Prior to this, it is necessary to obtain the isolation between the x channel and the y channel, the isolation between the x channel and the z channel, and the isolation between the y channel and the z channel in advance.

[0028] The calculation model expression of isolation degree is: ;in, Indicates i Channel and j Isolation between channels, Indicates j The average peak-to-peak value of the signal amplitude of each channel, Indicates 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 changes from the maximum value to the minimum value of the signal within a certain period of time. It is a key indicator for measuring the dynamic range and fluctuation amplitude of the signal.

[0029] Taking the x channel as an example, the amplitude value of the signal arriving at the channel output port is calculated as follows: 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 the x channel and the 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 .

[0030] 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.

[0031] Taking the x channel as an example, after obtaining the signal amplitude value at the x channel output port of one of the detectors, 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.

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

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

[0034] 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.

[0035] 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.

[0036] 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 established based on the historical signal amplitude values ​​of the detector, there may be a situation where the amplitude value of the current signal collected by other detectors is not within the range of each value interval established in step 1. In order to make the signal amplitude values ​​collected from other detectors comparable with the signal amplitude value of the detector, it is necessary to avoid the above situation. This embodiment takes the following measures: 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 include Figure 3 The following steps are shown: Step Q1: extract all signal amplitude values ​​that do not fall within any value interval in the value interval reference list.

[0037] 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.

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

[0039] The method for establishing a new numerical interval includes the following steps: 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.

[0040] For example, the three numerical intervals in Table 1 are sorted to obtain [(a1, a4), (a4, a7), (a7, a9)].

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

[0042] 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.

[0043] 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.

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

[0045] 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.

[0046] 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 < the extracted signal amplitude value < the lower limit of the current numerical interval, execute step C2; if the extracted signal amplitude value > the lower limit of the current numerical interval, execute step C3.

[0047] 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.

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

[0049] 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.

[0050] 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.

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

[0052] Taking the three amplitude values ​​a2, a5 and a8 extracted above as an example, after processing in step Q2, the following six new value intervals are newly established, namely: (a1, a2), (a2, a4), (a4, a5), (a5, a7), (a7, a7) and (a8, a9). Among them, a1<a2<a4<a5<a7<a8<a9.

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

[0054] 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.

[0055] Add the above 6 new value intervals to Table 1, and update Table 1 to obtain Table 2.

[0056] Numerical range Priority number (a1, a2) 1 (a2, a4) 2 (a4, a5) 3 (a5, a7) 4 (a7, a8) 5 (a8, a9) 6

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

[0058] 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, S7 is executed, otherwise Table 2 needs to be updated according to the method provided in steps Q1 to Q4.

[0059] For example, if the signal amplitude value of S2 is a6, then its numerical interval is (a5, a7), and the corresponding priority number is 4.

[0060] 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.

[0061] In summary, the signal processing method for tunnel engineering proposed in this embodiment first clusters the amplitude values ​​of the collected historical signals. The distribution characteristics of the amplitude values ​​of the historical signals can be obtained through clustering, that is, each cluster corresponds to an amplitude value distribution interval, and then the amplitude values ​​in each distribution interval can be graded, and the priority can be used as the basis for dynamically selecting the signal channel. Further, this method takes into account the influence of the isolation between the channels inside the detector on the signal transmission, and uses the amplitude value of the signal arriving at the output port of the channel as the standard for dynamically selecting the signal channel; the amplitude value at the output port is the amplitude value of the signal entering the analyzer, and the amplitude value will affect the three-dimensional imaging result of the analyzer. Finally, the priorities corresponding to the amplitude values ​​of the signals arriving at the output port at the same time are compared, and the channels where the highest priority signal is located are kept unblocked, and the channels where the remaining lower priority signals are located are closed, so as to eliminate the low amplitude signals, thereby ensuring that the high amplitude signals enter the analyzer and improving the three-dimensional imaging results of the analyzer.

[0062] Embodiment 2: Corresponding to the above-mentioned embodiment 1, this embodiment proposes a detector for tunnel engineering, including: A data extraction module, used 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, used for clustering the signal amplitude values ​​in the first data set to obtain a plurality of clusters; List building module, used to obtain the value interval of each cluster and build a reference list of value intervals; A first sequence number allocation module, used to allocate a corresponding priority sequence number to each numerical interval according to the median of the numerical interval; A time stamp module, used to obtain the initial signal amplitude value when the current signal arrives at the sensor, and add an arrival time stamp to the current signal; An amplitude calculation module, used to obtain 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 is used to obtain 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 mark; 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 is used to determine whether the priority number in the second 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; A first channel control module, used to keep the channel unobstructed; A second channel control module, used for closing the channel; A first isolation degree acquisition module, used to acquire the isolation degree between the x channel and the y channel; A second isolation degree acquisition module, used to acquire the isolation degree between the x channel and the z channel; A third isolation degree acquisition module, used to acquire the isolation degree between the y channel and the z channel; The calculation model expression of the isolation degree is: ;in, Indicates i Channel and j Isolation between channels, Indicates j The average peak-to-peak value of the signal amplitude of each channel, Indicates i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3; A data transfer module, used for storing the signal amplitude value at the output port into the data storage module; A data storage module, used for storing the data sent by the data transfer module; A second analysis control module is used 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 so, 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, used for adding all new value intervals to the value interval reference list; The second sequence number allocation module is used to reallocate the corresponding priority sequence number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

[0063] The interval component modules include: An interval sorting unit, 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, used for traversing the numerical interval sequence, and controlling the first analysis control unit to work each time a numerical interval is accessed; a first analysis control unit, used to determine whether the currently accessed numerical interval is the first numerical interval; if it is the first numerical interval, control the first numerical comparison unit to work; 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, control the second numerical comparison unit to work; if it is the last numerical interval, control the third numerical comparison unit to work; A first numerical comparison unit is used 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 work; 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 work; 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, used for comparing the extracted signal amplitude value with the upper limit of the current numerical interval, controlling the interval traversal unit to access the next numerical interval, and controlling the operation of the first analysis control unit; The second numerical comparison unit is used 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, the second numerical interval construction unit is controlled to work; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, the third analysis control unit is controlled to work; 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, used to drive the second analysis control unit to work; 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, the second 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 fourth numerical comparison unit is controlled to operate; The fourth numerical comparison unit is used 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, the third numerical interval construction unit is controlled to work; 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.

[0064] 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.

[0065] 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 people familiar with this technology to understand and read, and are not used to limit the limiting conditions for the implementation of the present invention, so they have no substantial technical significance. Any modification of the structure, change of 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 effects and purposes that can be achieved by the present invention. At the same time, the terms such as "upper", "lower", "left", "right", "middle" and the like cited in this specification are only for the convenience of description, and are not used to limit the scope of the implementation of the present invention. The change or adjustment of the relative relationship should also be regarded as the scope of the 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 multiple historically stored signal amplitude values ​​from the data storage module to establish a first data set; clustering the signal amplitude values ​​in the first data set to obtain multiple clusters; obtaining the value interval of each cluster and establishing a reference list of value intervals; assigning a corresponding priority number to each value interval according to the median of the value 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 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: 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 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 in S5 is the maximum value in the third data set; if so, keep the current channel open; otherwise, close the current channel.

2. A signal processing method for tunnel engineering according to claim 1, characterized in that: Before S2, the following steps are also included: obtaining the isolation between the x channel and the y channel; obtaining the isolation between the x channel and the z channel; and obtaining the isolation between the y channel and the z channel.

3. A signal processing method for tunnel engineering according to claim 2, characterized in that: The calculation model expression of the isolation degree is: ;in, Indicates i Channel and j Isolation between channels, Indicates j The average peak-to-peak value of the signal amplitude of each channel, Indicates i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3.

4. A signal processing method for tunnel engineering according to any one of claims 1 to 3, 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.

5. A signal processing method for tunnel engineering according to any one of claims 1 to 3, 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 one of the numerical intervals in the numerical interval reference list; if so, execute S4; otherwise, execute S3.2 to S3.5; S3.2: extracting all signal amplitude values ​​that do not fall within any value interval in the value interval reference list; S3.3: Establish one or more new value intervals using all the extracted signal amplitude values ​​to obtain multiple new value intervals; S3.4: Add all new value intervals to the value interval reference list; S3.5: According to the median of the numerical interval, reallocate a corresponding priority number to each numerical interval in the numerical interval reference list.

6. A signal processing method for tunnel engineering according to claim 5, 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: 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 lower limit of the current numerical interval; if the upper limit of the previous numerical interval < the extracted signal amplitude value < the lower limit of the current numerical interval, execute step C2; if the extracted signal amplitude value > 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.

7. A signal processing method for tunnel engineering according to claim 6, 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.

8. A detector for tunnel engineering, characterized in that: include: A data extraction module, used 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, used for clustering the signal amplitude values ​​in the first data set to obtain a plurality of clusters; List building module, used to obtain the value interval of each cluster and build a reference list of value intervals; A first sequence number allocation module, used to allocate a corresponding priority sequence number to each numerical interval according to the median of the numerical interval; A time stamp module, used to obtain the initial signal amplitude value when the current signal arrives at the sensor, and add an arrival time stamp to the current signal; An amplitude calculation module, used to obtain 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 is used to obtain 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 mark; 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 is used to determine whether the priority number in the second 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; A first channel control module, used to keep the channel unobstructed; The second channel control module is used to close the channel.

9. A detector for tunnel engineering according to claim 8, characterized in that: Also includes: A first isolation degree acquisition module, used to acquire the isolation degree between the x channel and the y channel; A second isolation degree acquisition module, used to acquire the isolation degree between the x channel and the z channel; A third isolation degree acquisition module, used to acquire the isolation degree between the y channel and the z channel; The calculation model expression of the isolation degree is: ;in, Indicates i Channel and j Isolation between channels, Indicates j The average peak-to-peak value of the signal amplitude of each channel, Indicates i The maximum amplitude mean of the channels, i =1,2,3, j =1,2,3.

10. A detector for tunnel engineering according to claim 8 or 9, characterized in that: Also includes: A data transfer module, used for storing 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.

11. A detector for tunnel engineering according to claim 8 or 9, characterized in that: Also includes: A second analysis 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, the first serial number extraction module is controlled to work; otherwise, the second amplitude acquisition module, the interval construction module, the list update module or the second serial number allocation module is controlled 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, used for adding all new value intervals to the value interval reference list; The second sequence number allocation module is used to reallocate the corresponding priority sequence number to each numerical interval in the numerical interval reference list according to the median of the numerical interval.

12. A detector for tunnel engineering according to claim 11, characterized in that: The interval component module includes: An interval sorting unit, 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, used for traversing the numerical interval sequence, and controlling the first analysis control unit to work each time a numerical interval is accessed; a first analysis control unit, used to determine whether the currently accessed numerical interval is the first numerical interval; if it is the first numerical interval, control the first numerical comparison unit to work; 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, control the second numerical comparison unit to work; if it is the last numerical interval, control the third numerical comparison unit to work; A first numerical comparison unit is used 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 work; 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 work; The 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, used for comparing the extracted signal amplitude value with the upper limit of the current numerical interval, controlling the interval traversal unit to access the next numerical interval, and controlling the operation of the first analysis control unit; The second numerical comparison unit is used 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, the second numerical interval construction unit is controlled to work; if the extracted signal amplitude value is greater than the lower limit of the current numerical interval, the third analysis control unit is controlled to work; 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, used to drive the second analysis control unit to work; 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, the second 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 fourth numerical comparison unit is controlled to operate; The fourth numerical comparison unit is used 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, the third numerical interval construction unit is controlled to work; 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