Method for detecting underground pipe gallery

An underground pipe gallery and detection method technology, applied in the field of engineering geophysical exploration, can solve the problems of being easily affected by surrounding human activities, high-density electrical methods being easily affected by the underground environment, poor detection effect, etc.

Active Publication Date: 2019-01-11
WUHAN ENG SCI & TECH RESINST
10 Cites 1 Cited by

AI-Extracted Technical Summary

Problems solved by technology

However, these methods have their own shortcomings: the electromagnetic induction method and the ground radar method have better detection effects on metal pipelines, and poorer detection effects on non-metallic pipelines, especially PE pipes, which are widely used now. The seismic ...
View more

Abstract

The invention relates to a method for detecting an urban underground pipe gallery, comprising the following steps: Step 1, uniformly arranging a plurality of seismic wave collecting device groups along the pipe gallery on the ground above a to-be-detected area of the urban underground pipe gallery, wherein each seismic wave collecting device in the seismic wave collecting device groups is a sampling point in the collecting area; step 2, collecting seismic wave signals received by seismic wave collecting device groups at each collecting point of each collecting area; step 3, performing cross-correlation calculation by using the seismic wave signals of each collecting point in each collecting area; Step 4, calculating the corresponding empirical Green's function; Step 5, extracting the surface wave dispersion curve; Step 6, decomposing the fundamental mode and the higher-order mode of the dispersion curve; Step 7, obtaining a corresponding speed profile of the to-be-detected area; Step 8, determining the location of the pipe gallery. According to the method for detecting an urban underground pipe gallery, the urban underground pipe gallery can be conveniently and efficiently detected.

Application Domain

Seismic signal processing

Technology Topic

Surface waveSeismic wave +3

Image

  • Method for detecting underground pipe gallery
  • Method for detecting underground pipe gallery
  • Method for detecting underground pipe gallery

Examples

  • Experimental program(1)

Example Embodiment

[0041] The present invention will be further described in detail below in conjunction with the drawings and specific embodiments:
[0042] An urban underground pipe gallery detection method, such as Figure 1~5 As shown, it includes the following steps:
[0043] Step 1: According to the existing urban underground pipe gallery information, determine the direction of the urban underground pipe gallery to be detected, and evenly arrange multiple seismic wave acquisition device groups along the pipe gallery on the ground above the area to be detected of the urban underground pipe gallery. The setting position of each seismic wave acquisition device group is used as an acquisition area. The multiple seismic wave acquisition devices in each seismic wave acquisition device group are evenly arranged on the ground along the vertical direction of the pipe gallery. Each seismic wave acquisition device in the seismic wave acquisition device group collects The device is a sampling point in the collection area;
[0044] Step 2: Turn on the seismic wave acquisition device, and collect the seismic wave signals received by the seismic wave acquisition device group at each acquisition point in each acquisition area;
[0045] Step 3: Use the seismic wave signals of each acquisition point in each acquisition area to perform cross-correlation calculations, and calculate the seismic cross-correlation function between the most lateral acquisition point and the nth acquisition point in each acquisition area (equivalent to acquiring at the most lateral A seismic source is arranged in the device, and receivers are arranged in other acquisition devices to improve the accuracy of subsequent calculations. The specific formula is as follows:
[0046]
[0047] Among them, A 1 (t) represents the seismic wave signal of the most lateral acquisition point, A n (t) represents the seismic wave signal of the nth acquisition point, n is counted from the second acquisition point outside the most lateral acquisition point, and N is A 1 (t) and A n (t) The sum of the number of discrete data points of two seismic wave signals, m represents the variable, Indicates that m is taken from 1 to N-1 in the formula, and then summed, t represents the time, the seismic wave signal is a discrete signal, A 1 (m) represents the value of the m-th discrete point of the seismic wave signal at the most lateral acquisition point, A n (m+t) represents the value of the m+t discrete point of the seismic wave signal at the nth acquisition point;
[0048] Step 4: Use the seismic wave cross-correlation function obtained in step 3 to calculate the corresponding empirical Green's function (the real Green's function is difficult to obtain directly, here you can use the empirical Green's function to replace the real Green's function for calculation and processing), the specific operations are as follows :
[0049] The seismic wave cross-correlation function A between a most lateral acquisition point and the nth acquisition point in each acquisition area 1n (t), empirical Green's function And the real Green function G 1n The relationship between (t) is expressed as:
[0050]
[0051] Among them, Indicates the empirical Green's function when the geophone at the n-th collection point receives the signal as the seismic source, G 1n (t) represents the true Green's function when the geophone at the most lateral collection point is used as the seismic source and the geophone at the nth collection point receives the signal, Indicates the empirical Green's function when the detector at the nth acquisition point is used as the seismic source, and the detector at the most lateral acquisition point receives the signal, G n1 (t) represents the true Green's function when the detector at the nth acquisition point is used as the seismic source and the detector at the most lateral acquisition point receives the signal;
[0052] The empirical Green's function can be obtained by the time differentiation of the seismic cross-correlation function;
[0053]
[0054] A n1 (-t) represents the seismic wave cross-correlation function between the nth acquisition point and one of the most lateral acquisition points in each acquisition area;
[0055]
[0056] A n (m) represents the value of the m-th discrete point of the seismic wave signal at the n-th acquisition point, A 1 (m+t) represents the value of the m+t discrete point of the seismic wave signal of a most lateral acquisition point;
[0057] Step 5: Use the empirical Green's function obtained in Step 4 to extract the surface wave dispersion curve. The operation is as follows:
[0058] The result obtained in step 4 Combine it into:
[0059]
[0060] And Perform a two-dimensional Fourier transform to obtain a frequency wavenumber spectrum F(f,k), where f represents frequency and k represents wavenumber;
[0061] According to the definition of wave number
[0062]
[0063] Where v s Is the transverse wave velocity;
[0064] Transform the frequency wave number domain record to the frequency wave speed domain: use the frequency wave number domain record:
[0065]
[0066] Resample (i.e. by formula Convert the frequency wave number domain to the frequency wave speed domain) to obtain the frequency wave speed spectrum. At any frequency in the frequency wave speed domain, from the extreme point obtained from the frequency wave speed spectrum, the corresponding phase speed value of the frequency can be obtained, thereby obtaining the surface Dispersion curve of wave recording;
[0067] Step 6: Because there is a pipe gallery underground, the dispersion curve obtained in step 5 has multi-mode characteristics. Decompose the fundamental mode and higher-order mode of the dispersion curve to obtain the fundamental mode part and each higher-order mode part of the dispersion curve;
[0068] Step 7: Due to the characteristics of surface waves, the velocities of waves of different wavelengths correspond to the velocities of the corresponding depth positions underground, so based on the several sets of dispersion curves obtained in step 6, the corresponding velocity profile of the area to be detected can be obtained;
[0069] Step 8: Detect the area to be detected in the urban underground pipe gallery, and compare and analyze with the sections obtained in step 7, discard inconsistent sections, leave consistent sections, and find out from the sections Underground wave speed abnormal area, the underground wave speed abnormal area corresponds to the location of the underground pipe gallery.
[0070] In the step 8, the area to be detected of the urban underground pipe gallery is detected by the geological radar or the transient Rayleigh surface wave method.
[0071] In step 8, the anomalous area of ​​underground wave speed indicates an area where the wave speed is higher than the surrounding wave speed.
[0072] The most lateral collection point is the leftmost collection point or the rightmost collection point.
[0073] The distance between two adjacent seismic wave acquisition device groups ranges from 50 to 200 meters (in the straight-line traveling area, there are two accurate end points within this distance, and the specific location of the pipe gallery can be determined).
[0074] The distance between two adjacent seismic wave acquisition devices in each seismic wave acquisition device group ranges from 0.15 to 2 meters (the smaller the distance, the higher the measurement accuracy, but the use cost and efficiency issues need to be considered).
[0075] The content not described in detail in this specification belongs to the prior art known to those skilled in the art.

PUM

no PUM

Description & Claims & Application Information

We can also present the details of the Description, Claims and Application information to help users get a comprehensive understanding of the technical details of the patent, such as background art, summary of invention, brief description of drawings, description of embodiments, and other original content. On the other hand, users can also determine the specific scope of protection of the technology through the list of claims; as well as understand the changes in the life cycle of the technology with the presentation of the patent timeline. Login to view more.
Who we serve
  • R&D Engineer
  • R&D Manager
  • IP Professional
Why Eureka
  • Industry Leading Data Capabilities
  • Powerful AI technology
  • Patent DNA Extraction
Social media
Try Eureka
PatSnap group products