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

The averaging function p(z) = 5m should optimally be chosen to be compatible with the pulsewidth τ(z) = 50ns, which should be around half the interferometer length L₀ = 10m. With this width of the averaging function, it has no significant effect on the spatial resolution of the DAS. Modeling with and without weighted averaging is presented in Figure 1.9, which demonstrates that significant noise reduction can be achieved. It should be noted that this noise reduction is particularly marked in comparison with the coherent OTDR response, by contrasting with Figure 1.5. Nevertheless, weighted averaging suppresses rather than completely removes the effect of flicker noise, and some channels still demonstrate excessive noise (in addition to shot noise). Hence, the response over all depths at a given time for Figure 1.9 will contain spikes for faded channels.

As is explored in Section 1.3, the problem of flicker noise can be overcome by introducing engineered bright scatter zones along the fiber with constant spatial separation and uniform amplitude. Such scatter zones also reflect more photons, and so improve the shot noise detection limitation. In addition, the use of such engineered fiber allows the use of phase‐detection algorithms with improved sensitivity and extended dynamic range.

Figure 1.9 The left‐hand panel shows modeling of raw DAS acoustic data (Equation 1.12); the right‐hand panel shows the same shot with weighted averaging denoising (Equation 1.13) applied. 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.

      1.1.6. Time Integration of DAS Signal

A DAS interrogator measures, in accordance with Equation 1.13, the speed difference between two sections of fiber that are separated by interferometer length L₀ (referred to also as the gauge length), as presented in Figure 1.10. In pulse‐to‐pulse consideration, the DAS response is linearly proportional to the fiber elongation averaged over the gauge length in the nanometer scale, or strain rate in the nanostrain per second scale. The consideration can also be extended to multiple pulses by time integration of the DAS signal. So, if fiber rests initially and ground displacement equals to zero u(z, t₁) = 0, then:

      (1.22)

      meaning a time integrated DAS signal can be considered as an output of a huge caliper that is measuring fiber elongation between two points with sub‐nanometer precision. This measuring principle is different from that of a geophone but is similar to an electromagnetic linear strain seismograph that can measure changes in distance between two points on the ground (Benioff, 1935).

Figure 1.10 Illustration of two time‐consecutive measurements when DAS output is proportional to fiber elongation between two probe pulses.

1.2. DAS SYSTEM PARAMETERS AND COMPARISON WITH GEOPHONES

In this section, we consider how DAS parameters (such as spatial resolution), gauge length, frequency response, and SNR enable DAS to become an effective tool for seismic measurements. Field data are also presented, with DAS output compared to geophone data.

      1.2.1. DAS Optimization for Seismic Applications

      Distributed fiber sensors measure physical parameters of an external environment continuously through the integration properties of light traveling along a lengthy optical path. This is quite different from point sensors, such as geophones, which make an inertial measurement of ground speed at fixed positions (SEAFOM, 2018). The DAS records a local strain rate, which can be converted into particle velocity to allow direct comparison with geophone data. Following Jousset et al. (2018), we can approximately represent DAS signal A(z, t) via ground displacement u(z, t), where F_S is the DAS sampling frequency and L₀ is the gauge length.

(1.23)

      If F_S → ∞, L₀ → 0, then the DAS signal can be presented in a double differential form:

(1.24)

These simplified expressions (Equations 1.23–1.24) give us a qualitative sense of the DAS algorithm output. For a subsequent quantitative analysis, we shall need the detailed expression that was obtained in the previous section. Namely, for a nonzero interferometer gauge length L₀ and optical pulsewidth τ, averaged over random scattering DAS output, A(z) can be represented by Equation 1.15 in expanded form:

(1.25)

Скачать книгу