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

Figure 1.12 that DAS low frequency sensitivity is significantly lower than that of a geophone. Practically, however, this is not the case, as the geophone noise rises at low frequencies, and this can be characterized by some high‐pass (HP) filters that limit the range to frequencies around 10 Hz (see dotted line in Figure 1.13). However, DAS has the potential to increase the spectral response at low frequencies by increasing interferometer length—for example, from L₀ = 10m to L₀ = 30m (see Figure 1.13). So, potentially, the DAS response can be synthesized from two measurements: with a short interferometer gauge length to deliver high spatial frequency bandwidth, and a long one to deliver low frequency. As the result, full‐frequency coverage can be as good as from a geophone antenna, or possibly even better, as will be shown in a later SNR comparison. An additional advantage over geophones is the large dynamic range of DAS at low frequencies, which will be discussed later.

Figure 1.13 Low spatial frequency gain in DAS by using long interferometer.

      1.2.2. DAS Directionality in Seismic Measurements

      In the previous section, we analyzed the correspondence between DAS and geophones in the one‐dimensional case and found that “geophone‐style” velocity data can be extracted from DAS signals by spatial integration. However, in 3D analysis, we should consider that DAS is not a velocity sensor but a differential strain sensor. This is a fundamental difference: DAS can measure a component of 3D tensor (strain) but not 3D vector (velocity).

Directionality of the DAS response depends on the fiber optic cable configuration and the cable design, as the device itself is sensitive only to fiber elongation. We will start our consideration where the fiber is placed linearly inside a cable, with no slippage between fiber and cable, nor between the cable and the ground. In this case, fiber displacement will follow ground displacement, and sensitivity will depend on the relative position of fiber and seismic source. A similar mechanical principle was used for the electromagnetic linear strain seismograph to measure variations in the distance between two points of the ground (Benioff, 1935). DAS directional response with respect to incident angle Γ can be found by transformation of the strain tensor components with rotation using geometrical consideration. For a longitudinal (P) apparent wave, it will be cos²Γ, and for transversal (S) wave sin Γ cos Γ, similar to Benioff (1935) (see Figure 1.14). Detailed analysis and diagrams for Rayleigh and Love waves can be found in Martin et al. (2018).

      In vertical seismic profiling (VSP), in the vertical part of the well, both cable and seismic waves are in the same direction for near‐offsets, so the DAS is more sensitive to P‐waves, in which the acoustic displacement vector coincides with the fiber direction. In other applications, such as fracking, the microseismic source is usually on a side of the cable, so shear waves can be effectively detected.

Cable orientation is responsible not only for acoustic amplitude but also for acoustic spatial resolution, even for the same acoustic wavelength. The cable acts as an acoustic antenna where the signal varies rapidly in space if the P‐wave and cable direction coincide, but the signal remains the same over distance if the acoustic wave front is parallel to the cable. To take this effect into consideration, we need to expand the expression for acoustic wavelength K_z along the cable for Equations 1.29–1.30 as:

(1.31)

      For a harmonic wave, directionality will directly affect not only the spatial resolution but also the temporal frequency. After Fourier transfer in the time domain, Equation 1.30, in the absence of aliasing, can be presented as:

(1.32)

Figure 1.14 DAS with linear optical cable is more sensitive to P‐wave in VSP configuration and to S‐wave in microseismic events.

Figure 1.15 2D spectral representation on upgoing and downgoing acoustic waves: The left panel represents original signal and the right panel represents filtered signal, where the spatial filter shape is shown by a blue wavy line.

This expression (Equation 1.32) represents spatial filtering in the 2D Fourier domain as shown in Figure 1.15 for L₀ = 10m for upgoing and downgoing waves, together with white phase noise. As far as harmonic waves can be represented as single lines, the result of spatial filtering is an intensity modulation of these lines. Such modulation is equivalent to temporal frequency modulation, so we can combine Equation 1.31 and Equation 1.32 to get:

(1.33)

It is interesting to mention that, for uniform strain, where C → ∞, we have |ℑ_A(F)| ∝ cos²Γ as expected from Benioff (1935). The result of modeling of Equation 1.33 is presented in Figure 1.16 for different incident angles. An increase in angle expands the measurement frequency range but reduces low‐frequency SNR at the same time.

DAS directivity can be significantly modified through appropriate cable design, which is currently a developing area. For example, field tests have shown that the helical placement of fiber within a squeezable material can deliver omnidirectional sensitivity (Hornman et al., 2013) for P‐waves. The angular dependence