Distributed Acoustic Sensing in Geophysics. Группа авторов

      where F_S is sampling frequency and A₀ = 115nm is a scale constant (Equation 1.14). So, the velocity field can be recovered by spatial integration starting from a motionless point as:

(1.26)

Then DAS signal (Equation 1.25) can be transformed using shift invariant a(z₁) ⊗ b(z₁ + z₂) = a(z₁ + z₂) ⊗ b(z₁) to:

(1.27)

      where θ(z) is the Heaviside step function, whose value is zero for a negative argument. As expected, the DAS signal is represented (Equation 1.5) as a convolution of a point spread function with v(z).

Spatially integrated signal (Equation 1.27) was modeled for 10 m gauge length and 50 ns pulsewidth, as shown in Figure 1.5 (right panel). The results of modeling (Equation 1.25) are presented in Figure 1.11 (left panel), and the result is converted to geophone‐style data (Equation 1.26) in the right panel. From a practical point of view, low temporal frequencies, out of the range of interest, can be filtered out, and also spatial antialiasing filtering can be used. It is worth mentioning that the right panel of Figure 1.11 demonstrates the real change in polarity of the reflected seismic pulse. Also, spatial integration (Equation 1.26) acts as statistical averaging, which eliminates the randomness of the “staircasing” in Figure 1.5 left panel.

      The most valuable geophysical information is delivered by sound waves with frequencies below F_MAX = 150Hz, as higher frequencies are attenuated by the ground. For a speed of soundC = 3000m/s, this corresponds to an acoustic wavelength C/F_MAX = 20m, so Nyquist’s limit dictates that L_G ≤ C/2F_MAX = 10m is the maximum spacing of conventional sensors. Formally, the linear spline approximation G(z) of conventional antenna velocity v(z) output can be represented using expressions from (Unser, 1999), as:

      (1.28)

      The spatial spectral response of DAS in acoustic angular wavenumber K_z can be represented by Fourier transform ℑ(K_z) following Goodman (2005):

(1.29)

Such spectral responses can be normalized for a constant signal ℑ(K) = 1 (see black line in Figure 1.12). The comb function in (Equation 1.29) is responsible for the repeating of the spatial spectrum with a shift of 2π/Λ, as is shown by the dotted line. To prevent aliasing, the signal spectrum should be inside Nyquist’s limit, which is shown by the gray vertical line.

      Let us compare the conventional velocity sensor with the DAS spectrum, calculated from the spatial resolution expression (Equation 1.25), by Fourier transform as:

(1.30)

Two cases are presented in Figure 1.12: when the optical pulse length is almost equal to the interferometer gauge length τ = L₀, and when it is half the interferometer gauge length τ = L₀/2 (see dashed and solid blue lines, respectively). The absolute value is presented in the figure to aid comparison between curves. In the second case, we have a gain, which is highlighted by the blue filling. This gain can be explained by signal smearing over a long pulse.

Figure 1.11 Acoustic measurements using DAS: The left panel represents strain rate measurement and the right panel displays ground speed measurement, the transform to which comprises filtering and integration. The signals’ cross‐section along the white line is shown in the bottom panels in radians. The modeled source is shown in the right panel of Figure 1.5.

Figure 1.12 Comparison of DAS spectral response with that from a 10 m sensor antenna array.

It