The present application is a U.S. National Stage patent application of International Patent Application No. PCT/US2016/050076, filed on Sep. 2, 2016, the benefit of which is claimed and the disclosure of which is incorporated herein by reference in its entirety.
The present disclosure generally relates to oilfield equipment and, in particular, to downhole tools, drilling and related systems, and techniques for determining a change in an environmental condition in a wellbore. More particularly still, the present disclosure relates to systems and methods for calculating the change in the environmental condition in the wellbore based on a change in an optical frequency from a frequency of a trace measured before the change to a frequency of a trace measured after the change, where the environmental condition can be at least one of temperature, pressure, and strain.
The Joule-Thomson effect causes temperature of a fluid (e.g. a gas or liquid) to change if it is pushed through a throttle, orifice, choke, or similar device, while preventing heat exchange between the fluid and the environment (that is, via an adiabatic process). The strength of the Joule-Thomson effect may depend on the particular fluid, its phase or phases, its composition (for solutions and mixtures), pressure and temperature. The effect is characterized by a coefficient called the Joule-Thomson coefficient and both the sign and magnitude of this coefficient vary with pressure, temperature, and composition of the fluid. In the case of gases in typical ambient conditions a significant cooling occurs when the gas undergoes a pressure change, such as in most refrigerators and air-conditioning systems. In the case of liquids that undergo a pressure change, the Joule-Thomson effect can be much weaker, which can result in a much smaller change in a temperature of the liquid.
Because pressure changes in a gas, when compared to a liquid, may cause more significant temperature changes, these changes can be detected and monitored by systems such as a Distributed Temperature Sensing (DTS) system, which can use Rayleigh, Raman, and/or Brillion backscattering techniques with an optical waveguide positioned in a wellbore to measure environmental conditions (such as temperature, etc.) in the wellbore. However, the change in temperature for liquids due to a pressure change, flow rate, etc. is not readily detected using the DTS systems because the changes in temperature can be below the resolution of these DTS systems.
Therefore, it will be readily appreciated that improvements in the arts of determining changes in environmental conditions in a wellbore are continually needed.
Various embodiments of the present disclosure will be understood more fully from the detailed description given below and from the accompanying drawings of various embodiments of the disclosure. In the drawings, like reference numbers may indicate identical or functionally similar elements. Embodiments are described in detail hereinafter with reference to the accompanying figures, in which:
The disclosure may repeat reference numerals and/or letters in the various examples or Figures. This repetition is for the purpose of simplicity and clarity and does not in itself dictate a relationship between the various embodiments and/or configurations discussed. Further, spatially relative terms, such as beneath, below, lower, above, upper, uphole, downhole, upstream, downstream, and the like, may be used herein for ease of description to describe one element or feature's relationship to another element(s) or feature(s) as illustrated, the upward direction being toward the top of the corresponding figure and the downward direction being toward the bottom of the corresponding figure, the uphole direction being toward the surface of the wellbore, the downhole direction being toward the toe of the wellbore. Unless otherwise stated, the spatially relative terms are intended to encompass different orientations of the apparatus in use or operation in addition to the orientation depicted in the Figures. For example, if an apparatus in the Figures is turned over, elements described as being “below” or “beneath” other elements or features would then be oriented “above” the other elements or features. Thus, the exemplary term “below” can encompass both an orientation of above and below. The apparatus may be otherwise oriented (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein may likewise be interpreted accordingly.
Moreover even though a Figure may depict a horizontal wellbore or a vertical wellbore, unless indicated otherwise, it should be understood by those skilled in the art that the apparatus according to the present disclosure is equally well suited for use in wellbores having other orientations including vertical wellbores, slanted wellbores, multilateral wellbores or the like. Likewise, unless otherwise noted, even though a Figure may depict an offshore operation, it should be understood by those skilled in the art that the method and/or system according to the present disclosure is equally well suited for use in onshore operations and vice-versa. Further, unless otherwise noted, even though a Figure may depict a cased hole, it should be understood by those skilled in the art that the method and/or system according to the present disclosure is equally well suited for use in open hole operations.
As used herein, the words “comprise,” “have,” “include,” and all grammatical variations thereof are each intended to have an open, non-limiting meaning that does not exclude additional elements or steps. While compositions and methods are described in terms of “comprising,” “containing,” or “including” various components or steps, the compositions and methods also can “consist essentially of” or “consist of” the various components and steps. It should also be understood that, as used herein, “first,” “second,” and “third,” are assigned arbitrarily and are merely intended to differentiate between two or more objects, etc., as the case may be, and does not indicate any sequence. Furthermore, it is to be understood that the mere use of the word “first” does not require that there be any “second,” and the mere use of the word “second” does not require that there be any “first” or “third,” etc.
The terms in the claims have their plain, ordinary meaning unless otherwise explicitly and clearly defined by the patentee. Moreover, the indefinite articles “a” or “an,” as used in the claims, are defined herein to mean one or more than one of the element that it introduces. If there is any conflict in the usages of a word or term in this specification and one or more patent(s) or other documents that may be incorporated herein by reference, the definitions that are consistent with this specification should be adopted.
Generally, this disclosure provides a method and system to detect small changes in downhole environmental conditions in a wellbore. For example, temperature changes on the order of 1 to 10 mK (millikelvin) and small strain changes can be detected. The system can include an optical waveguide, a precision frequency optical light source coupled to the waveguide, and an optical receiver coupled to the waveguide. The optical light source can produce an optical signal at a substantially constant frequency which can then be externally tuned using optical modulation techniques to precisely control the optical frequency of the light. The optical waveguide can be installed in a wellbore on a conveyance vehicle (e.g. tubing strings, slickline, wireline, etc.). The method and system of this disclosure can detect small changes in downhole environmental conditions because the optical waveguide's properties are sensitive to changes in the environmental conditions. The optical light source can be coupled to the waveguide for launching light at a substantially constant frequency into the waveguide and the optical receiver can be coupled to the waveguide to receive backscattered light from the waveguide. Because of the response of the waveguide's properties to environmental changes, the backscattered light contains information about how the waveguide's properties were modified by the environment. Therefore, this backscattered light can be used to create one or more measurement traces that represent(s) a profile of an environmental condition along at least a portion of a length of the wellbore. The measurement trace(s) can be created by interference between backscattered light from different inherent scatterers within the waveguide.
When an environmental condition changes, the frequency of the optical light source can be adjusted to compensate for the change in the environmental condition. In other words, the measurement trace obtained when the waveguide “is not” subject to a given environmental condition can be recuperated by using light at a different frequency when the waveguide “is” subject to the given environmental condition. The frequency shift that compensates for the environmental condition change can be determined when a measurement trace at the adjusted frequency substantially correlates with a measurement trace obtained prior to the change of the environmental condition. When the measurement traces, which are taken before and after the environmental change, substantially match each other (i.e. have a high correlation value from 0.8-1.0), then the change in the environmental condition can be determined based on the difference between the optical frequencies used to produce the before and after traces. Furthermore, by determining the change in the environmental condition, other aspects of the fluid flowing in the wellbore may also be determined, such as fluid type, fluid composition, fluid flow rate, pressure drop between the formation and the wellbore, watercut, location(s) of production zone(s) in the wellbore, etc.
It should be understood that the optical waveguide, light source, and backscattered light receiver of this disclosure can also be used to perform additional tasks such as detecting temperature profiles via a Distributed Temperature Sensing (DTS) system using backscattering techniques (e.g. Raman, Rayleigh, and/or Brillouin), detecting absolute pressure measurements downhole, detecting acoustic signals via Distributed Acoustic Sensing (DAS), transmitting command and control data to/from downhole equipment, transmitting collected sensor data and/or telemetry data to/from downhole equipment, etc. These tasks can be performed along with the method and system of the current disclosure.
Turning to
The wellbore production system 10 in
Wellbore 12 may be formed of single or multiple bores, extending into the formation 14, and disposed in any orientation (e.g. vertical, inclined, horizontal, combinations of these, etc.). The wellbore production system 10 can also include multiple wellbores 12 with each wellbore 12 having single or multiple bores. The rig 18 may be spaced apart from a wellhead 40, as shown in
An optical light source/receiver device 52 can be coupled to an optical waveguide 50 installed along the tubing string 30 in the wellbore 12. The optical light source/receiver 52 can be a device known as a Coherent-Optical Time Domain Reflectometry (C-OTDR) device, a Coherent-Optical Frequency Domain Reflectometry (C-OFDR) device, or any other suitable device for launching light into the optical waveguide at a substantially constant frequency and receiving backscattered light from the waveguide (e.g. Rayleigh backscattering). As used herein, the terms “substantially constant frequency” and “stable frequency” refer to a frequency of the optical light source 122 that may vary up to +/−10 MHz long term. Frequencies that remain within this tolerance are seen as being “substantially constant” and/or “stable.” A representative block diagram of a device 52 is shown in
The optical waveguide 50 is shown in
The device 52 and optical waveguide 50 can be used to collect baseline traces at one or more optical signal frequencies. The device 52 can launch an optical light signal 54 at a substantially constant frequency into the waveguide and receive a backscattered light signal 56 from the waveguide 50 (e.g. using Rayleigh backscattering). The device 52 can produce a trace that represents an intensity of the backscattered light signal 56 along a distance of the waveguide 50. One of more baseline traces can be produced that are representative of an environmental condition (e.g. temperature, pressure, strain, etc.) in the wellbore 12 at one or more of the locations of the production zones 60, 62, 64, and/or other locations along the wellbore. The device 52 can be a C-OTDR device.
Polarization controllers 142, 144 can be used to optimize the signal strength at the balanced photo detector 146. The pulsed optical signal 54 can then be amplified (e.g. by an erbium-doped fiber amplifier (EDFA)) and launched into the waveguide 50. As a backscattered light signal 56 is received from the waveguide 50, the signal 56 can be combined at the coupler 148 with the local oscillator signal from the coupler 140 and then detected by the balanced photo diode (PD) receiver 146. The controller 120 can receive a digitized version of the detected signal and produce a representative trace of the signal 56 and transmit the trace to an operator, where the trace can represent the intensity of the backscattered light along at least a portion of the wellbore 12. The C-OTDR device 52 can also produce a trace that is a time resolved intensity of Rayleigh backscattering. Therefore, it is seen that several configurations of the device 52 can be used to support the principles of this disclosure.
If the environmental absolute temperature T is close to some initial temperature T0; then this means that T=T0+ΔT with ΔT<<T0, where ΔT is a change in temperature from the initial temperature T0 to the new temperature T. Similarly, if the optical signal frequency υ of the laser light is close to some initial frequency υ0 (=υlaser+ΔυSSB
A correlation value Corr, can be used to indicate a correlation between two separate traces. If A1(n) is the first trace and A2(n) is the second trace (with n being the index for a specific data point, for example n=1 . . . 10000 if a total of N=10000 data points are measured), then the correlation between the two traces can be defined by equation (1),
where
is the mean of A1(n), and
is its standard deviation, and likewise for μA
In an approximation, for the case of the environmental perturbation affecting only temperature, the trace can depend only on a distance d and a quantity δ, where
Therefore, this indicates that if the temperature changes, the frequency can be fine-tuned afterwards until the quantity δ becomes zero in which case the trace can be substantially equal to the initial trace. In the case where the first and second traces A1(n), A2(n) have a correlation value that is close to “1,” then the first and second traces may be seen as being substantially equal to each other. As used herein, a correlation value “Corr” that is greater than or equal to 0.80 indicates that the correlation is “close to 1” and that the two traces being correlated are substantially equal to each other. With the indication from the correlation value Corr that the two traces are substantially equal to each other, then the quantity δ may be assumed to be zero and the change in temperature can be determined from the frequency change in by solving for ΔT, which results in equation (2):
Therefore, the change in temperature can be calculated based on the change in frequency used to cause the second trace to have a correlation value in the range of 0.8 to 1.0 with the first trace. This example assumes that a strain on the fiber is substantially constant and does not change significantly during the determination of the change in temperature. Therefore, the portion of the change in frequency necessary to compensate for any change in strain of the optical waveguide is considered to be negligible for the purposes of calculating the change in temperature based on the change in frequency between the first and second traces. It should be noted, that the before and after traces refer to two traces that have a correlation Corr from 0.8 and 1.0. A plurality of traces can be collected at the initial and final environmental conditions, with each of the traces collected at the initial conditions correlated to each of the traces collected at the final conditions, thereby providing many independent measurements of the desired frequency shift.
In the example above, traces 90, 94, a temperature shift of −28 mK caused a noticeable difference between the trace 90 produced prior to the temperature change during time period t1, and the trace 94 produced after the temperature change during time period t2. As per the discussion above, the frequency of the optical signal launched into the optical waveguide 50 from the device 52 can be fine-tuned to compensate for the change in the environmental temperature.
Therefore, the change in frequency can be used to calculate the change in the environmental temperature. Accordingly, equation (2) shows a calculation for ΔT based on the change in frequency. As equation (3) shows, the change in frequency of 39 MHz results in a calculated temperature change of 29 mK, which very closely approximates the actual temperature change of −28 mK measured by the temperature probe in this example.
By making similar measurements with the waveguide 50 positioned in a wellbore 12, a change in temperature in a liquid at each of the production zones 60, 62, 64, can be calculated by determining the change in frequency of the optical signal that respectively compensates for the amount of change in the environmental temperature at each of the production zones 60, 62, 64 as fluid flows from the formation 14 into the wellbore 12.
A change in strain of the fiber can be similarly calculated from a change in frequency, assuming the temperature remains constant during the measurements to determine the strain change. Much like the process for determining a change in temperature, the change in strain can have a linear relationship with the change in frequency of the optical signal 54 which is launched into the waveguide 50. The change in the strain can be calculated when two traces taken at different optical signal frequencies (one taken before the strain change and one taken after the strain change) have a correlation value Corr between 0.8 to 1.0. The difference in the optical signal frequencies can be used to calculate the change in the strain. The laser frequency change Δυ that compensates for the strain change Δε is given by equation (4):
Solving for Δε yields equation (5):
In step 206, a backscattered signal 56 can be received from the optical waveguide 50 in response to the introduction of the optical signal 54 during the first time period. A first trace can be produced that represents an intensity of the backscattered light along the waveguide 50 during the first time period. This first trace can also be seen as a baseline trace (i.e., a trace representing the initial conditions) that provides a baseline to which other traces are compared. It should also be noted that this first trace can include multiple traces at multiple different frequencies to provide a wide range of backscattered signals 56 with which to compare against future traces. It is not required that the first, initial, or baseline trace be a singular trace.
After the initial trace or traces are collected in step 206, the environmental condition(s) can be changed in step 208. There are many ways one or more of the environmental conditions downhole can be changed in keeping with the principles of the current disclosure. For example, opening or closing a valve 42 (see
If the valve 42 is then partially (or fully) closed, a backpressure from the closed valve 42 can cause pressure P1 in the wellbore to equalize with pressures P2, P3, P4 in the formation 14 at the respective production zones 60, 62, 64. Equalizing these pressures, reduces (or eliminates) the flow of the fluids 70, 72, 74 into the wellbore 12. With reduced fluid flow, the Joule-Thompson effect in each of the fluids 70, 72, 74 is reduced (or eliminated), thereby allowing the temperatures of the fluids 70, 72, 74 to equalize with the temperatures of the formation 14 at the respective production zones 60, 62, 64. Equalizing the temperatures (and pressures) at the production zones 60, 62, 64, can cause a change in temperature of each of the fluids 70, 72, 74 by reduction and/or elimination of the Joule-Thompson effect in these fluids.
Alternatively, the environmental condition change can also be caused by having the valve 42 initially closed and then opening the valve 42 to permit fluid flow from the wellbore 12. With the valve 42 initially closed, flow of fluid 76 is prevented (or significantly restricted), thereby preventing (or significantly restricting) flow of fluids 70, 72, 74 from the formation 14 into the wellbore 12 at the production zones 60, 62, 64. Without a pressure drop between the production zones 60, 62, 64 and the formation 14, the temperatures and pressures at these zones can equalize with the formation 14 temperatures and pressures to provide a stable environmental condition profile along the wellbore 12. The initial trace(s) can be collected while the environmental conditions (e.g. temperature, pressure, strain, etc.) are stable and the valve 42 is closed. After collection of the initial trace(s), then the environmental condition can be changed by opening the valve 42. Opening the valve 42, can allow fluid 76 to flow through the wellbore 12, thereby creating a pressure differential in the fluids 70, 72, 74 at the respective production zones 60, 62, 64. The pressure differential in each of the fluids 70, 72, 74 can cause the temperature in the fluids to change due to the Joule-Thompson effect.
Once the environmental change has occurred, a next trace can be collected by adjusting the frequency of the light source 122 to a next frequency in step 210 that is different from the first frequency, introducing an optical signal 54 from the light source 122 into the optical waveguide 50 at the next frequency during a next time period in step 212, and receiving a next backscattered signal 56 from the optical waveguide 50 in response to the introduction of the optical signal 54 during the next time period in step 214. The next trace, which represents an intensity of the next backscattered signal 56 along the waveguide 50 during the next time period, can be produced and transmitted to the controller 60 for analysis. The initial (or first) and next (or second) traces can be compared in step 216 to determine if there are any differences between them. This comparison can be performed by calculating the correlation value (i.e. Corr) between the initial and next traces and determining if the resulting correlation value is within an acceptable range. As given above, if the correlation value is within a range from 0.8 to 1.0, then the signals can be deemed as being substantially the same. The steps 210, 212, 214, 216 can be repeated as many times as needed to produce a next (or second) trace that is substantially equal to the initial (or first) trace. In other words, the next frequency can be adjusted (step 210), the next trace collected for the adjusted frequency (steps 212, 214), and the next (or second) trace compared to the first trace (step 216) as many times as desired to produce a next trace that is substantially the same as the first trace (or at least substantially the same as one or more of the first traces). Additionally, the next trace can also include multiple traces, where, in step 216, one of the first traces can be compared to one or more of the next traces until one or more of the first traces is deemed to be substantially equal one or more of the next traces.
When step 216 determines that at least one initial trace is substantially equal to at least one next trace, then the difference between the frequency of the substantially equal initial and next traces can be used to calculate the change in the environmental condition in step 218, such as using equation (2) to calculate a change in temperature, and/or equation (5) above to calculate a change in strain. When the change in the environmental condition has been determined, then one or more fluid characteristics (such as fluid type, fluid composition, fluid flow, fluid pressure differential, etc.) can be determined based on the environmental condition change. For example, if temperature is the environmental condition that changed at a wellbore location during step 208, a strain of the waveguide 50 remains substantially constant, the wellbore volumes at the location are known, and wellbore and formation pressures at the location are known, then the temperature change determined in step 218 can be used to determine an effective Joule-Thompson coefficient for the fluid flowing into the wellbore at the wellbore location. The effective Joule-Thompson coefficient can be used to determine a composition of the fluid flowing into the wellbore 12 at the location, or fluid flow rate (if fluid composition is known), etc. If an actual Joule-Thompson coefficient of the fluid composition is known, then one or more of the other characteristics (such as fluid pressure differential, fluid flow rate, change in strain, etc.) can be determined based on the difference between the frequency of the initial trace and the next trace which has sufficient correlation.
There are several approaches for collecting initial and next traces to identify a change in an environmental condition in the wellbore 12. One example of an approach can be to measure one initial trace at a single frequency, change the flow rate, and then measure multiple traces at many different frequencies. The initial trace can then be compared to each one of the multiple traces to identify which one (or more) of the multiple traces correlate best with the initial trace. The difference in frequency between the initial trace and the one (or more) of the multiple traces, which best correlate, can be used to calculate the change in the environmental condition.
As way of another example, many initial traces at multiple different frequencies can be collected prior to a change in the environmental condition, after which the environmental condition can be changed, and then one final trace at one single frequency can be measured. Each of the initial traces can be compared to the final trace to determine which one or more of the initial traces correlate with the final trace. The difference in frequency between the one or more of the initial traces and the final trace can be used to calculate the change in the environmental condition.
As way of another example, one initial trace at one single frequency can be collected prior to a change in the environmental condition, after which the environmental condition can be changed, and then new traces can be continuously collected and compared to the initial trace as the optical signal frequency is slowly changed. When the comparison between the initial trace and one or more of the new traces identifies the new traces that best correlate to the initial trace, then the frequency difference between these initial and one or more new traces can be used to calculate the change in the environmental condition.
As way of yet another example, many initial traces at multiple different frequencies can be collected prior to the change in the environmental condition, after which the environmental condition can be changed, and many final traces at many different frequencies can be collected. Each of the initial traces can be compared to each one of the final traces to determine which one or more of the initial traces best correlate with one or more of the final traces. By plotting the result as a 2D color map, or a 2D shade level map (as done in
As way of yet another example,
Therefore, as long as the time intervals between the datasets are short enough to ensure that an environmental condition change between two contiguous times can be correlated to determine a frequency shift, the incremental condition change between those contiguous points can be calculated and thus the overall condition change can be calculated. Because each condition increment is determined with a non-zero uncertainty, by adding up the condition increments to obtain the full condition change, a cumulative effect of the uncertainty can occur (following a standard propagation of errors for a sum of measurements). The uncertainties can be minimized by having a device 52 with a large frequency span (2 S) and that datasets be collected at certain condition increments, so as to minimize the number of steps between the current condition shift measurement and the baseline condition.
In step 312, it is determined whether or not the traces in the Current Dataset (C) and the Baseline Dataset (B) can be correlated. Please note that the correlation of these two datasets can be performed as described in detail in this disclosure. If the traces in the Current Dataset (C) and the Baseline Dataset (B) can be correlated (i.e. YES), then the method 300 proceeds to step 314 where an incremental frequency shift Δυ(C,B) can be determined and saved off for future calculations. The incremental shift Δ(C,B) refers to the shift between the Current Dataset (C) and the Baseline Dataset (B), where C and B change per the Current counter (C) and the Baseline counter (B) values. In step 316, it is determined whether or not the environmental condition change event has ended. If YES, then the method 300 proceeds to step 324, where all the incremental frequency shifts Δv(C,B) that were saved off in step 314 are combined (i.e. added together) to produce the total frequency shift between the Baseline Dataset (0) and the last Current Dataset (C). From this total frequency shift, the total environmental condition can be determined, and from that the fluid characteristic(s) can be determined, as described in detail in this disclosure.
If the answer in step 316 is NO, then the method 300 returns back to step 308, were C is incremented by 1 (i.e. C=C+1), and a new Current Dataset (C) is produced in step 310. Steps 308, 310, 312, 314, and 316 are repeated until the Current Dataset (C) cannot correlate to the Baseline Dataset (B) (i.e. NO in step 312), or the change event has ended (i.e. YES in step 316).
If the answer to step 312 is that the Current Dataset (C) cannot be correlated to the Baseline Dataset (B) (i.e. NO in step 312), then in step 318, it is determined whether or not the Current counter (C) is equal to the Baseline counter (B)+1. If YES, then an error has occurred because the environmental condition changes too fast during the change event. In this case, parameters will need to be adjusted to prevent this error and the method 300 reran. If the answer to step 318 is NO, then, in step 320, the Baseline counter (B) is incremented by “1” (i.e. B=B+1), and an incremental Baseline Dataset (B) is set equal to the Current Dataset (C−1) and saved off for future calculations. The method 300 then proceeds to step 312 to verify that the Current Dataset (C) can be correlated to the Baseline Dataset (B). If so, then the steps 314 and 316 can be repeated.
The process of comparing the traces, calculating the correlation values, and identifying the best correlation values can be automated to provide expedited results. The following discussion describes various procedures for determining a pair of traces that have the best correlation out of an array of traces. Baseline data refers to traces collected before an event that changes an environmental condition, and subsequent data refers to traces collected after the event. Please note that the subsequent traces can also be collected while the condition is changing. These procedures and their examples focus on temperature as the environmental condition that changes between time period t1 and time period t2. However, other environmental conditions can be determined in a similar fashion. These procedures are only a few of possible procedures for determining the best correlation values between before and after traces.
A First Procedure:
The first procedure can be used for an approach that has baseline data collected at a single frequency, with subsequent data collected over a range of frequencies. This first procedure can be used to compare a baseline trace A1(fref, n) to a set of M traces A2(fm, n), where the baseline trace A1(fref, n), is collected at frequency f=fref, having N elements (1≤n≤N), and can be representative of a C-OTDR pattern for the measurement that started at time period t=t1. Please note that, in principle, index n is also representative of time, but at a much smaller scale (e.g. intervals less than 1 ns, typically) compared to measurement time intervals Δtmeas of 1 minute or more between the datasets A1 and A2. Also note that baseline trace A1(fref, n) can be obtained from the average of several traces, all taken at f=fref and over a time period from t=t1 to t=t1+Δtmeas.
The set of M traces A2(fm, n), with (1≤m≤M), are each N elements long (1≤n≤N), and are collected starting at t=t2. Please note that the trace for each frequency fm can also consist of an average of many traces, collected at the same frequency, so as to reduce the noise for the trace of that frequency. It is assumed that a range of frequencies can be chosen large enough to cover the effect of the largest temperature shift expected (or desired to be measurable) between t1 and t2. This range of frequencies {fm} can be chosen such that the reference frequency fref is in the middle of the range, but it is not required for the reference frequency fref be in the middle of the range.
For this first procedure, a vector of M cross-correlation can be calculated using equation (6):
with (1≤m≤M) and where
is the mean of A1(fref, n), and
is its standard deviation. Likewise, for the data collected at t=t2, there are M mean values
and M standard deviations
Unless otherwise stated, the correlation function is calculated using the sum
which implies that all data points from 1 to N are part of a physical region of interest (i.e. the optical waveguide 50). If a smaller region of interest is desired, the sum can run only over points which are part of the region of interest.
The set of data points (fm, y2,m) can have a peak centered at frequency f2,peak which can be extracted using a standard peak finding algorithm. A simple way can be to select the frequency fm that corresponds to the highest value y2,m. However, other peak finding algorithms can involve fitting a curve (e.g., a parabola) over the peak portion of the data and taking f2,peak as the location of the peak determined mathematically from the equation of the fitted curve. The offset can be given as Δf1→2=f2,peak−fref. This procedure can be illustrated by the following example for the FIRST PROCEDURE.
An example for the FIRST PROCEDURE:
Optical signals (e.g. laser pulses) can be introduced into the waveguide 50 and a backscattered signal 56 can be measured. During this initial measurement, the SSB modulation frequency can be set to 4.500 GHz.
After 5 minutes, multiple backscattered signals 56 were collected at many different SSB modulation frequencies.
These signals A2(fSSB, n) are then compared to the initial signal A1(4.500 GHz, n) by calculating the correlation function. Equation (7) calculates the correlation value Corr between the region of interest (data point 7500 to data point 10300) of the initial trace measured at the first time and the region of interest of each of the traces recorded at the second time. Here, μ and σ represent the respective median and standard deviation.
As illustrated by the plot of correlation values vs. SSB frequencies shown in
A Second Procedure:
The second procedure can be used for an approach that has baseline data collected over a range of frequencies, with subsequent data also collected over a range of frequencies. This second procedure can be used to compare a set of P traces A1(fp, n), to a set of M traces A2(fm, n), where the P traces A1(fp, n), are collected at a range of frequencies f∈ {fp}, having N elements (1≤n≤N), and can be representative of C-OTDR patterns for the measurements that started at time period t=t1. Please note that the trace for each frequency fp can also consist of an average of many traces, to reduce the noise in the trace for each frequency.
As in the first procedure above, the set of M traces A2(fm, n), with (1≤m≤M), are each N elements long (1≤n≤N), and are collected starting at t=t2. Please note that the trace for each frequency fm can also consist of an average of many traces. Again, it is assumed that a range of frequencies can be chosen large enough to cover the effect of the largest temperature shift expected (or desired to be measurable) between t1 and t2. This range can be chosen such that the reference frequency fref is in the middle of the range, but it is not required for the reference frequency fref be in the middle of the range.
For this second method, P vectors of M cross-correlation values can be regarded as a P×M matrix, and these values can be calculated using equation (8):
For each p there is a corresponding curve made up of the points (fm, y2,p,m) and a location of the peak for each of these curves, expressed as f2,p,peak, can be determined. Thus resulting in a set of P points (fp, f2,p,peak). Those points can describe a line: f2,p,peak=bfp+a. Therefore, the problem can be reduced to a standard linear regression problem with an offset given as Δf1→2=a. (Note that the experiment should yield b=1 and any deviation from this can be interpreted as due to error in the experiment and can be used as a data quality measure.) With the offset being the difference in frequency between a before trace and an after trace, and with the frequency difference, the change in temperature can be calculated as before.
A Third Procedure:
This third procedure incorporates elements that can make it more robust for automated processing because it incorporates quantitative data quality criteria. The third procedure can be used for an approach that has baseline data collected over a range of M frequencies {fm} and subsequent data collected over the same range of frequencies. This third procedure can be used to compare a first set of M traces A1(fm, n) to a second set of M traces A2(fm, n), where the first M traces A1(fm, n) are collected at a range of frequencies f∈{fm}, with (1≤m≤M), and having N elements (1≤n≤N), which can be representative of C-OTDR patterns for the measurements that started at time period t=t1. Additionally, the second set of M traces A2(fm, n) are collected at a range of frequencies f∈{fm}, with (1≤m≤M), and having N elements (1≤n≤N), which are representative of C-OTDR patterns for measurements that started at time period t=t2. Where the first and second set of M traces have a range of frequencies with a span S=max(fm)−min(fm).
For this third procedure, M vectors of M cross-correlation values can be regarded as an M×M matrix, and these values can be calculated using equation (9):
where
is the mean of A1(fi, n) and
is its standard deviation, and likewise for j, μA
A set of K quantities Qk are calculated using Qk=Σm=1M Corrk,m where Corrk,m, is obtained via an interpolation process from the Corri,j defined above. To understand how this interpolation is done, it should be understood that for each k value, there is a fixed offset Bk, such that, for each value of m, is also associated a frequency fk=fm+Bk. Depending on the values of fm and of Bk, fk may fall inside or outside of the range of frequencies collected. If inside the range, fk may not fall on one of the values of {fm} in use, then select fj ∈{fm} such that fj≤fk and fj+1>fk, then Corrk,m is obtained using Equation (10):
If fk falls outside of the range of frequencies collected, then set Corrk,m=0. If other words, Corrk,m=0 if fk<min(fm) or if fk>max(fm).
A normalized quality factor {tilde over (Q)}k can be calculated by dividing the quality factor Qk by the number of correlation values (however, excluding those which were set to zero because fk<min(fm) or fk>max(fm). In other words,
with w(k)=1 for min(fm)<fk<max(fm) and w(k)=0 otherwise. A value {tilde over (k)} can be determined by determining which k value maximizes the quality factor {tilde over (Q)}k. Additionally, the frequency shift which maximizes a correlation between before and after traces is given as B{tilde over (k)}. From this, the temperature change which occurred between t1 and t2 can be expressed as ΔT=aB{tilde over (k)} where the quantity
as stated previously. Operationally, data quality measures can be implemented to verify a validity of the results. For example, the result could be rejected if {tilde over (Q)}k is less than a certain threshold (e.g. 0.7), or if |Bk| is above a threshold (e.g. 0.8 S).
From these three procedures, it has been shown how the before and after traces having the best correlation to each other can be identified, even with each set of before and after traces including multiple traces. It should be understood that these are only a few of the procedures that can be used to determine the before and after traces with the best correlation values, and thereby determine a change in an environmental condition in a wellbore, where the environmental condition can be at least one of temperature, pressure, strain, etc.
Thus, a method for detecting environmental changes in a wellbore (or downhole) has been described. Embodiments of the method may generally include positioning (or installing) an optical waveguide 50 along or within a wellbore 12, the optical waveguide 50 being coupled to an optical laser light source 122, 52 that can launch an optical signal 54 into the waveguide 50, and the waveguide 50 being coupled to a receiver 52 that receives a backscattered light signal 56 from the optical waveguide, where the backscattered light signal 56 represents an intensity of backscattered light along a length of the waveguide 50, and thus along a length (or segment) of the wellbore 12. One or more optical signals can be introduced (or launched) into the waveguide 50 at one or more optical signal frequencies and during one or more time periods, thereby resulting in one or more backscattered signals 56 being received by the receiver 52, which produces a trace for each of the one of more backscattered signals 56, where the trace of each backscattered light signal 56 can represent an environmental condition downhole along a length (or segment) of the wellbore.
By causing a change in the environmental condition (e.g. change in temperature, pressure, strain, etc.), additional backscattered light signals 56 can be obtained at one or more frequencies after the environmental condition change has occurred. Comparing the traces generated before the condition change to those generated after the change, can identify a before trace and an after trace that are substantially equal to each other, with the after trace having a difference in frequency from the before trace. This frequency difference can be used to determine the amount of change in the environmental condition that occurred when the environmental change happened.
Other embodiments of the method of detecting a change in an environmental condition in a wellbore may generally include the features given above, as well as generating one or more traces before the condition change event, and generating one or more traces after the change event. The after event traces can be used to compare to the before event traces to identify a before trace and an after trace that are substantially equal to each other, with the after trace having a difference in frequency from the before trace. This frequency difference can be used to determine the amount of change in the environmental condition that occurred when the environmental change event happened. The change in the environmental condition can also be used to identify a location of a production zone 60, 62, 64 in a wellbore 12, when the production zone 60, 62, 64 locations are unknown, which can be the case even in an uncased wellbore 12 without casing string 34 and perforations 36. The method can produce a profile of the environmental condition along the waveguide and therefore, correspondingly identify a location in the wellbore 12 of variations of the environmental condition with respective production zone locations being identified by these variations.
The environmental change event can be opening and/or closing a valve 42, where the valve 42 can be positioned at the surface or at any point within the wellbore 12, such as interconnected in a tubing string 30, deposited in the wellbore 12, positioned in a casing string 34, positioned within tubing string 30, positioned above and/or below production zones 60, 62, 64. The valve 42 can be a one-time use valve with a degradable material, where the valve 42 is initially closed and can be opened by degrading (e.g. melting, dissolving, eroding, etc.) the material to allow fluid flow through the valve. The valve 42 can include swellable material that actuates the valve 42, where swelling of the swellable material can be used to directly restrict fluid flow, or can be used to actuate the valve 42 to selectively permit and prevent fluid flow through the valve 42.
For the foregoing embodiments, the method may include any one of the following steps, alone or in combination with each other:
Determining the amount of change in the environmental condition can include, calculating a difference in frequency between the first and next frequencies, and calculating the change in the environmental condition based on the calculated difference in frequency. Adjusting the optical signal 54 frequency to a next frequency, introducing the optical signal 54 at the next frequency, receiving the next backscattered signal in response to the introduction of the next frequency optical signal, and identifying the differences can be repeated until the first trace substantially equals the next trace. The intensity of the backscattered signal 56 can represent the environmental condition at locations along the wellbore 12. The environmental condition can be at least one of temperature, pressure, and strain. The first dataset collected prior to the change in the environmental condition can be multiple traces, with the traces being generated in response to multiple different frequencies of the optical signal 54 being launched into the waveguide 50, and multiple resulting backscattered signals 56 being received by the receiver 52 from the waveguide 50, where each of the multiple frequencies is different from the other ones of the multiple frequencies. Identifying differences between each one of the first traces to the next trace, and determining the change in the environmental condition when at least one of the first traces is substantially equal to the next trace. The change in the environmental condition can be caused by opening or closing a valve 42 to respectively increase or decrease fluid flow into the wellbore 12 from one or more production zones 60, 62, 64. Determining the amount of change in the environmental condition can include, at least one of a differential fluid pressure, a fluid flow rate, and a fluid composition based on the determined change in the environmental condition. The wellbore can comprise multiple segments, and the first and second traces can represent the environmental condition along a length of at least one of the multiple segments.
Additionally, another embodiment of the method may generally include positioning (or installing) an optical waveguide 50 along (or within) a wellbore 12, the optical waveguide 50 being coupled to an optical laser light source 122, 52 that can launch (or introduce) a first optical signal 54 into the waveguide 50 at a first frequency during a first time period t1, and the waveguide 50 being coupled to a receiver 52 that receives a first backscattered light signal 56 from the optical waveguide 50 in response to launching the first optical signal 54 in the first time period t1, and producing a first trace 90 that represents an intensity of backscattered light along a length of the waveguide 50, and thus along a length (or segment) of the wellbore 12. Changing the environmental condition and introducing a second optical signal 54 from the light source 122, 52 into the optical waveguide 50 at a second frequency during a second time period t2 while the environmental condition is changing. Receiving a second backscattered signal 56 from the optical waveguide 50 in response to the introduction of the second optical signal 54 and producing a second trace 94. Comparing the first trace 90 to the second trace 94, determining that the first and second traces 90, 94 correlate to each other, and determining a first incremental change in the environmental condition based on differences between the first and second frequencies. Introducing a third optical signal 54 from the light source 122, 52 into the optical waveguide 50 at a third frequency during a third time period t3 which is after the environmental change has occurred and the environmental change is stable, receiving a third backscattered signal 56 from the optical waveguide 50 in response to the introduction of the third optical signal 54 and producing a third trace 96. Comparing the second trace 94 to the third trace 96, determining that the second and third traces 94, 96 correlate to each other, determining a second incremental change in the environmental condition based on differences between the second and third frequencies, and determining the change in the environmental condition by summing the first and second incremental changes.
Additionally, another embodiment of the method may generally include positioning an optical waveguide 50 along the wellbore 12, introducing first optical signals 54 from a light source 52, 122 into the optical waveguide 50 at multiple first frequencies during a first time period, receiving first backscattered signals 56 from the optical waveguide 50 in response to the introduction of the first optical signals 54 and producing a set of first traces, making a set of baseline traces equal to the set of the first traces, initiating a change in the environmental condition, introducing second optical signals 54 from the light source 52, 122 into the optical waveguide 50 at multiple second frequencies after at least a portion of the environmental condition has occurred. Receiving second backscattered signals 56 from the optical waveguide 50 in response to the introduction of the second optical signals 54 and producing a set of second traces, comparing each one of the baseline traces to each one of the second traces, determining that at least one of the baseline traces correlates to at least one of the second traces, determining an incremental change in the environmental condition based on differences between frequencies that are associated with the at least one of the baseline traces and the at least one of the second traces. The method can also include adjusting the multiple second frequencies, to a new set of frequencies to prepare for more sets of traces to be produced, if the environmental change event has not completed (i.e. the environmental condition is still changing).
If the condition is still changing, then the method can repeat 1) the introducing the second optical signals, 2) the receiving the second backscattered signals, 3) the comparing the baseline and the second traces, 4) the determining the correlation, 5) the determining the incremental change and 6) the adjusting the multiple second frequencies until the environmental condition is stable. When the environmental condition is stable, the incremental environmental condition changes can be combined to determine a total environmental change in the wellbore 12.
For the foregoing embodiments, the method may include any one of the following steps, alone or in combination with each other:
The first, second, and third traces 90, 94, 96 can represent an intensity profile of the first, second, and third backscattered signals 54, respectively, along at least a portion of the wellbore 12. Determining at least one of fluid type, fluid composition, fluid flow, and fluid pressure differential based on the determined change in the environmental condition. The method can also include determining at least one of the group consisting of fluid type, fluid composition, fluid flow, and fluid pressure differential based on the determined change in the environmental condition, and/or determining multiple incremental changes in the environmental condition during the changing of the environmental condition, and determining the change in the environmental condition by summing the first, second, and multiple incremental changes.
The methods can also include determining that none of the sensing baseline traces correlate with any of the third traces, making the set of sensing baseline traces equal to a last set of the third traces that included at least one of the third traces that did correlate to at least one of the sensing baseline traces, determining at least one of the group consisting of fluid type, fluid composition, fluid flow, and fluid pressure differential based on the total environmental condition change in the wellbore. After the environmental condition change has occurred and the environmental condition is substantially stable, introducing fourth optical signals from the light source into the optical waveguide at multiple fourth frequencies, receiving fourth backscattered signals from the optical waveguide in response to the introduction of the fourth optical signals, producing a set of fourth traces corresponding to the reference portion, comparing each one of the reference baseline traces to each one of the fourth traces, and determining that at least one of the reference baseline traces correlates to at least one of the fourth traces. Calculating a correction value based on a difference between frequencies that are associated with the at least one of the reference baseline traces and the at least one of the fourth traces, and calculating a compensated value for the change in the environmental condition by removing the correction value from the total environmental condition change.
Additionally, a system for detecting a change in an environmental condition has been described. Embodiments of the system may generally include an optical waveguide 50 installed (or introduced) in a wellbore 12, a light source 122, 52 that introduces (or launches) an optical laser light signal 54 into the waveguide 50, and a receiver 52 that can receive a backscattered light signal 56 which represents an intensity of the environmental condition along at least a length (or segment) of the wellbore 12. The optical laser light signal 54 can be introduced into the waveguide 50 initially at a first frequency during a first time period with the receiver producing a first trace in response to reception of a backscattered light signal 56 from the waveguide 50. The optical laser light signal 54 can be introduced into the waveguide 50 at a second frequency during a second time period with the receiver producing a second trace in response to reception of a backscattered light signal 56 from the waveguide 50. The first and second traces are taken at different environmental conditions in the wellbore and the first and second traces are substantially equal to each other. The second frequency is different than the first frequency and the difference between the first and second frequencies can be used to determine the change in the environmental condition that occurred in the wellbore 12. Other embodiments of the system may generally include the features of the system above, except that at least one of the first and second traces can include multiple traces, and one of each of the multiple first and multiple second traces substantially equal each other. The environmental condition change that can be determined from the difference in frequency between the one of each of the multiple first and multiple second traces can be used to determine the environmental condition change in the wellbore.
For any of the foregoing embodiments, the system may include any one of the following elements, alone or in combination with each other:
The change in the environmental conditions can be due to a valve 42 that is selectively opened and closed, which can variably restrict fluid flow through the wellbore 12. The valve 42 can be closed (or at least partially closed) for the first time period and opened for the second time period. Alternatively, the valve 42 can be opened for the first time period and closed (or at least partially closed) for the second time period. The environmental condition change can also be caused by other suitable events, as well. At least one of a differential fluid pressure, a fluid flow rate, and a fluid composition of a fluid flowing from a production zone 60, 62, 64 into the wellbore 12 can be determined based on the difference between the first and second frequencies. The traces 94, 96, 98 are coherent optical time domain reflectometry traces. The wellbore can include multiple segments, and the first and second traces can represent environmental conditions along a length of at least one of the multiple segments. The first frequency can include multiple frequencies, the first trace can include multiple traces, and the first time period can include multiple time periods, with each of the first traces associated with one of the first time periods and one of the first frequencies, wherein one of the first traces can be substantially equal to the second trace, and the change in the environmental condition along the wellbore can be calculated based on the difference between the frequency associated with the one of the first traces and the second frequency.
Additionally, or alternatively, the second frequency can include multiple frequencies, the second trace can include multiple traces, and the second time period can include multiple time periods, with each of the second traces associated with one of the second time periods and one of the second frequencies, wherein one of the second traces can be substantially equal to the one of the first traces, and the change in the environmental condition along the wellbore can be calculated based on the difference between the frequency associated with the one of the first traces and the frequency associated with the one of the second traces.
Additionally, another embodiment of the system may generally include an optical waveguide positioned along the wellbore, a light source that introduces an optical signal into the optical waveguide at a first frequency during a first time period, a receiver that receives a backscattered signal from the optical waveguide in response to the introduction of the optical signal during the first time period, and processing circuitry that can perform operations comprising, producing a first trace which represents an intensity of the backscattered signal along the waveguide, causing a change in the environmental condition, adjusting a frequency of the light source to a next frequency which is different from the first frequency, introducing an optical signal from the light source into the optical waveguide at the next frequency during a next time period, receiving a next backscattered signal from the optical waveguide in response to the introduction of the optical signal during the next time period, producing a next trace which represents an intensity of the next backscattered signal along the waveguide during the next time period, identifying differences between frequencies associated with the first and next traces, and determining the change in the environmental condition based on the differences when the first and next traces are substantially equal to each other.
The operations can also include calculating a difference in frequency between the first and next frequencies, and calculating the change in the environmental condition based on the calculated difference in frequency.
The operations can also include repeating the adjusting, the introducing the optical signal at the next frequency, the receiving the next backscattered signal in response to the introduction of the optical signal at the next frequency, and the identifying the differences until the first trace substantially equals the next trace. The intensity can represent the environmental condition along the wellbore, and the environmental condition can be temperature.
The first frequency can include multiple frequencies and the first trace can include multiple traces with each one of the first traces corresponding to a separate one of the first frequencies, and with each of the first frequencies being different from other ones of the first frequencies. The operations can also include identifying differences between frequencies associated with each one of the first traces and the next trace, and determining the change in the environmental condition based on the differences when at least one of the first traces is substantially equal to the next trace.
The operations can also include opening a valve to increase fluid flow into the wellbore from a production zone, and/or closing a valve to decrease fluid flow into the wellbore from a production zone. The operations can also include calculating at least one of a differential fluid pressure, a fluid flow rate, and a fluid composition based on the determined change in the environmental condition. Along the wellbore can include multiple segments, where the first and next traces represent the environmental condition along a length of at least one of the multiple segments.
Furthermore, the illustrative methods described herein may be implemented by a system comprising processing circuitry that can include a non-transitory computer readable medium comprising instructions which, when executed by at least one processor of the processing circuitry, causes the processor to perform any of the methods described herein.
Although various embodiments have been shown and described, the disclosure is not limited to such embodiments and will be understood to include all modifications and variations as would be apparent to one skilled in the art. Therefore, it should be understood that the disclosure is not intended to be limited to the particular forms disclosed; rather, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the disclosure as defined by the appended claims.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2016/050076 | 9/2/2016 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2018/044317 | 3/8/2018 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
7499151 | Bao et al. | Mar 2009 | B2 |
7599047 | Zou et al. | Oct 2009 | B2 |
8144314 | Yamamoto | Mar 2012 | B2 |
8520195 | Ramos et al. | Aug 2013 | B2 |
9170149 | Hartog et al. | Oct 2015 | B2 |
20050274513 | Schultz et al. | Dec 2005 | A1 |
20070165691 | Taverner et al. | Jul 2007 | A1 |
20080130707 | Yamamoto | Jun 2008 | A1 |
20110090496 | Samson | Apr 2011 | A1 |
20130257435 | Smithson | Oct 2013 | A1 |
20150146759 | Johnston | May 2015 | A1 |
Number | Date | Country |
---|---|---|
WO 1999045235 | Sep 1999 | WO |
WO 2009158630 | Dec 2009 | WO |
WO 2015076976 | May 2015 | WO |
WO 2016007161 | Jan 2016 | WO |
WO 2016028288 | Feb 2016 | WO |
Entry |
---|
International Search Report and the Written Opinion of the International Search Authority, or the Declaration, dated May 29, 2017, PCT/US2016/050076, 10 pages, ISA/KR. |
Koyamada et al. “Fiber-Optic Distributed Strain and Temperature Sensing with Very High Measurand Resolution Over Long Range Using Coherent OTDR.” Journal of Lightwave Technology. vol. 27, No. 9. pp. 1142-1146. May 2009. |
Lee et al. “Optical Single Sideband Signal Generation Using Phase Modulation of Semiconductor Optical Amplifier.” IEICE Phtonics Technology Letters. vol. 16, No. 5. pp. 1373-1375. May 2004. |
Loranger et al. “Rayleigh Scatter Based Order of Magnitude Increase in Distributed Temperature and Strain Sensing by Simple UV Exposure of Optical Fibre.” Scientific Reports. vol. 5. Jun. 2015. |
Lu et al. “MilliKelvin Resolution in Cryogenic Temperature Distributed Fibre Sensing Based on Coherent Rayleigh Scattering.” Proc. SPIE 23rd International Conference on Optical Fibre Sensors. vol. 9157. |
OzOptics 2012: Fiber Optic Distributed Strain and Temperature Sensors. http://www.amstechnologies.com/fileadmin/amsmedia/downloads/4549_oz_optics_fo_dstssensor.pdf. |
OzOptics 2016: Fiber Optic Distributed Strain and Temperature Sensors. http://www.ozoptics.com/allnew_pdf/dts0115.pdf. |
Sang et al. “One Centimeter Spatial Resolution Temperature Measurements in a Nuclear Reactor Using Rayleigh Scatter in Optical Fiber.” IEEE Sensors Journal. vol. 8, Issue No. 7. Jul. 2008. |
Imahama et al. “Restorability of Rayleigh Backscatter Traces Measured by coherent OTDR with Precisely Frequency-Controlled Light Source.” IEICE Trans. Commun. vol. E91-B. No. 4. pp. 1243-1246. Apr. 2008. |
Koyamada et al. “Novel Fiber-Optic Distributed Strain and Temperature Sensor with Very High Resolution.” IEICE Trans. Commun. vol. E89-B. No. 5. pp. 1722-1725. May 2006. |
Number | Date | Country | |
---|---|---|---|
20210231006 A1 | Jul 2021 | US |