This disclosure relates generally to the detection of low frequency vibrations using laser interferometry. More particularly it pertains to systems, methods, and structures for low frequency vibration (i.e., seismic) detection exhibiting an improved laser phase noise tolerance through the use of fiber turnarounds.
As is known, the detection of low frequency vibrations using ultra-stable laser interferometry is difficult, due in part to laser phase noise, particularly 1/f technical noise, which can swamp (overwhelm) any low frequency vibration(s) that one is trying to detect. A conventional method is to use an ultra-stable laser having a linewidth on the order of 1 Hz, taking special care to isolate the laser chamber mechanically and thermally to reduce the 1/f technical noise. Lasers exhibiting such characteristics are extremely expensive, and mechanical isolation may be operationally difficult depending upon an environment in which it is operating.
An advance in the art is made according to aspects of the present disclosure directed to laser interferometric systems, methods, and structures exhibiting superior laser phase noise tolerance particularly in seismic detection applications.
In sharp contrast to the prior art, systems, methods and structures according to aspects of the present disclosure advantageously relax laser requirements by exhibiting a novel configuration wherein the same laser which generates an outgoing signal is coherently detected using the same laser as local oscillator. According to further aspects of the present disclosure, the use of fiber turnarounds allows for the cancellation and/or mitigation of undesired mechanical vibration.
A more complete understanding of the present disclosure may be realized by reference to the accompanying drawing in which:
The illustrative embodiments are described more fully by the Figures and detailed description. Embodiments according to this disclosure may, however, be embodied in various forms and are not limited to specific or illustrative embodiments described in the drawing and detailed description.
The following merely illustrates the principles of the disclosure. It will thus be appreciated that those skilled in the art will be able to devise various arrangements which, although not explicitly described or shown herein, embody the principles of the disclosure and are included within its spirit and scope.
Furthermore, all examples and conditional language recited herein are intended to be only for pedagogical purposes to aid the reader in understanding the principles of the disclosure and the concepts contributed by the inventor(s) to furthering the art and are to be construed as being without limitation to such specifically recited examples and conditions.
Moreover, all statements herein reciting principles, aspects, and embodiments of the disclosure, as well as specific examples thereof, are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such equivalents include both currently known equivalents as well as equivalents developed in the future, i.e., any elements developed that perform the same function, regardless of structure.
Thus, for example, it will be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative circuitry embodying the principles of the disclosure.
Unless otherwise explicitly specified herein, the FIGs comprising the drawing are not drawn to scale.
Those skilled in the art will understand and appreciate that the complex-valued baseband photocurrent takes the form of:
I(t)=R√{square root over (PrPLO)}exp(j(ϕvib(t)+ϕL(t)−ϕR(t)))+n(t), (1)
where R is the responsivity of the photoreceiver, Pr and PLO are the powers of the received signal and LO at the input of the coherent receiver, respectively. It is assumed that the lasers included in both interrogators—are operating and outputting light at substantially the same wavelength λL=λR=λ, where ϕL(t) and ϕR (t) are phase noises of the left (L) interrogator (transmitter) and right (R) interrogator (receiver-local oscillator), respectively, n(t) is a photocurrent arising from all the noise sources in the system, including amplified spontaneous emission (ASE) of all the inline amplifiers, as well as any shot noise and any thermal noise in the coherent receiver.
The variable of interest is ϕvib(t), which is the optical phase arising from cumulative strain along the fiber, and takes the form of:
where Ltot is the length of the optical cable, and Δϵ(z) is the tensile strain at position z along the fiber optic cable.
A usual method of recovering ϕvib(t) is to take the unwrapped phase of the photocurrent:
ϕ(t)=∠I(t)=ϕvib(t)+ϕL(t)−ϕR(t)+ϕAWGN(t). (3)
where ϕAWGN is the angle formed by the projection of n(t) onto a vector perpendicular to R√{square root over (PrPLO)}exp(j(ϕvib(t)+ϕL(t)−ϕR (t))).
Due to signal propagation delay, the onset of rapid phase variation will only be apparent to the receiver at time tR=t+(L−zvib)/(c/neff), where t is the instance at which the seismic wave first impinges the optical cable at position zvib from the transmitter, which is the closest point to the epicenter. neff is the effective index of the fiber. To determine zvib, a bidirectional link can be implemented: suppose an identical system as above propagates from right to left as shown in
The onset of rapid phase variation will occur for the L interrogator at tL=t+zvib/(c/neff). The time difference tR−tL=(Ltot−2zvib)/(c/neff) can then be used to infer zvio since 4″ is known. Using bidirectional transmission for multiple submarine cables allows determination of the epicenter of an earthquake by triangulation. Ultrastable laser interferometry relies on: (i) ultra-narrow linewidth lasers, and (ii) clock synchronization of all interrogators, which may be achieved by synchronization with respect to global positioning system (GPS).
One key to determining the onset of rapid phase variation is that ϕvib(t) can be distinguished from cumulative noise sources ϕL(t)−ϕR(t)+ϕAWGN(t) Performance can be analyzed in the frequency domain as shown in
Laser frequency noise can be modeled as a two-sided power spectral density (PSD) of:
where Δv is the Lorentzian linewidth (in Hz) arising from spontaneous emission inside the laser cavity. At low frequencies, however, a laser's frequency spectrum is usually dominated by 1/f “technical” noise, which can be interpreted as slow drift of the laser's center frequency. The fitting parameter f0 is the frequency at which frequency noise due to technical noise becomes equal to spontaneous emission. Since laser phase ϕPN(t)=2π ∫−∞t v(t′) dt′ is the integral of frequency noise, the two-sided PSD of laser phase noise is:
while the two-sided PSD of ϕAWGN(t) is:
where η is the signal-to-noise (SNR) ratio of the received signal.
From
We shall now show and describe in greater detail, systems, methods, and structures according to aspects of the present disclosure advantageously employ fiber turnarounds to improve laser phase noise tolerance in low frequency (seismic) detection using ultra-low linewidth interferometry.
As we have previously noted, a conventional, prior-art configuration includes a transmitter and a receiver positioned at opposite ends of a submarine optical fiber link, such that the transmitter laser and local oscillator laser are two different lasers.
In sharp contrast, systems, methods, and structures according to aspects of the present disclosure incorporate a fiber turnaround at the far end of the submarine optical fiber link such that the outgoing (transmitted) signal is routed back to the same transmitter side of the submarine fiber optic cable thereby advantageously allowing a single laser to serve as transmitter and local oscillator in a coherent interferometry arrangement. It will be appreciated by those skilled in the art that such an arrangement helps reduce phase noise at low frequencies where seismic vibration energy is concentrated, thereby allowing increased sensitivity, increased sensing range, and/or the use of less expensive lasers in the interrogators.
Examples of fiber turnarounds incorporated into a submarine cable may be observed in
As shown in that figure are two interrogators positioned at opposite ends of an undersea optical cable. Each of the interrogators L Interrogator and R Interrogator—include a laser light source operating at a specific wavelength for the L Interrogator and R Interrogator that are in optical communication with the undersea optical cable. Along with the laser(s), each of the interrogators include a coherent receiver which is in optical communication with the optical cable. As shown further, the optical cable includes at least a pair of separate optical fibers that are configured directionally—one from the L Interrogator to the R Interrogator and conversely one from the R Interrogator to the L Interrogator.
Advantageously, we improve the performance of detecting ϕvib(t) by using the system architecture as shown illustratively in
The phase of the received signal is thus ϕL,Rx(t)=ϕL,Tx(t−Trt)+(ϕAWGN(t), where Trt=2neffL/c is the round-trip delay, while ϕL,LO(t)=CL,Tx(t). Their coherent beat product produces a phase of ϕ(t)=PL(t−Trt)−ϕL(t)+ϕAWGN(t).
Ignoring the contribution by ASE noise, this phase has a Fourier transform of:
Φ(f)=ΦPN(f)[1−e−j2πfT
Thus PSD of phase noise in this scheme is:
Sϕϕ(f)=4Sϕ
As compared with the original scheme, which produces a phase noise PSD of Sϕϕ(f)=2Sϕ
Phase noise suppression is observed for frequencies below ≈½πTrt=16 Hz. At 1 Hz and 10 Hz, the phase noise suppression ratio (2πfTrt)2 are −24 dB and 4 dB, respectively. Since seismic vibration has its energy concentrated mostly between 1 and 10 Hz, the proposed scheme improves laser phase noise tolerance by the same magnitude.
It is possible to further extend the concept of using a fiber return to reject common-mode mechanical vibration, such as near landing stations in submarine links where mechanical vibration is strongest near the shore. As those skilled in the art will appreciate, to detect a seismic vibration on the ocean floor, it may be necessary to cancel the strong vibrations near the shore.
We note that in the configuration shown in
The optical phase due to mechanical strain at the two frequencies are:
Note that ϕvib,1(t), is the undesired phase due to cumulative mechanical strain between the transmitter and the shore. ϕvib,1(t) has the same form as the first term of ϕvib,2 (t) apart from a scaling factor equal to the ratio of their wavelengths λ1/λ2. It is possible to coherently detect both tones and retrieve their phases ϕ1(t)=ϕvib,1(t)+ϕTx(t−2neffLT/c)−ϕTx(t)+ϕAWGN,1t) and ϕ2=ϕvib,2 (t)+ϕTx (t−2neffL/c)−ϕTx(t)+ϕAWGN,2. The difference ϕ(t)=ϕ2(t)−(λ1/λ2)ϕ1(t) will be free of undesirable landing station vibration ϕvib,1(t). Note this method potentially increases the impact of phase error due to laser phase noise and ASE noise, but if vibration near the landing station, ϕvib,1(t), is the dominant source of interference over seismic frequencies of interest as shown in the configuration of
Note that suppression of near-shore vibration is possible because the two tones traveled in the same fiber cores (both outbound and inbound directions), so they experience identical mechanical strains. Note further that turning the outbound signal around just before the landing station at the remote end will also mitigate mechanical noise the far-end landing station.
As noted previously, such configuration allows the interrogator to perform coherent detection using the same laser as that which generated the outgoing signal. Thought of alternatively, we thus convert laser interferometry using two independent lasers into delay interferometry. As shown by Eq. (7), phase noise is suppressed at low frequencies where seismic vibration energy is concentrated, relaxing the requirement on the phase noise spectrum of the laser. The use of fiber turnarounds in
While we have presented this disclosure using some specific examples, those skilled in the art will recognize that our teachings are not so limited. Accordingly, this disclosure should be only limited by the scope of the claims attached hereto.
This disclosure claims the benefit of U.S. Provisional Patent Application Ser. No. 63/023,283 filed May 12, 2020 the entire contents of which is incorporated by reference as if set forth at length herein.
Number | Name | Date | Kind |
---|---|---|---|
20090059996 | Komeda | Mar 2009 | A1 |
20120039360 | MacDougall | Feb 2012 | A1 |
20170153155 | Hull | Jun 2017 | A1 |
Entry |
---|
6. Marra, Giuseppe, et al. “Seismology with optical links: enabling a global network for submarine earthquake monitoring.” arXiv preprint arXiv:1801.02698 (2017) (Year: 2017). |
7. Marra, Giuseppe, et al. “Ultrastable laser interferometry for earthquake detection with terrestrial and submarine cables.” Science 361.6401 (2018): 486-490 (Year: 2018). |
Number | Date | Country | |
---|---|---|---|
20210356613 A1 | Nov 2021 | US |
Number | Date | Country | |
---|---|---|---|
63023283 | May 2020 | US |