The present invention relates to the field of the assessment and condition monitoring of fluid conduits. In particular, aspects of the invention are methods and apparatus which use pressure surges for assessment and monitoring of the internal condition of fluid conduits, including the build-up and deposition of scale, sand, waxes and other materials on the interior surface of conduits. The invention has particular application to pipelines and other conduits in the oil and gas exploration, production and transportation industry.
It is known to use pressure surges for investigating conditions along pipelines. Well known examples in the water industry include the detection and location of leaks (e.g. Taghvaei et al, 2008 [Ref #8], Duan, et al, 2011 [Ref #2]), blockages and restrictions (e.g. Arbon, et al, 2008 [Ref #1], Lee et al, 2008 [Ref #6]). In addition, inverse transient methods have been developed for assessing characteristics of pipeline networks (e.g. Kapelan et al, 2003 [Ref #5], Shamlooa & Haghighi, 2009 [Ref #7]). All of these techniques have enjoyed some success in practical applications.
In the oil industry, surges are sometimes used to determine flowrates or to infer pressure distributions along a line. For example, if there is a valve at the end of a pipe and the area of the pipe is uniform (at least in the region close to the valve), by closing the valve rapidly the resulting pressure rise will satisfy the Joukowski relationship and so can be regarded as a measure of the velocity of the fluid before closure. Assuming that the cross-sectional area of the pipe is known, it is therefore possible to measure flowrates with good accuracy. Although practical issues limit the acceptability of the method, it can be accurate and it requires no special equipment except a pressure sensor. One example is described in international patent publication number WO 96/15427 [Ref #3], in which the estimation of flowrates through the measurement of Joukowski pressures following intentional valve closure is described.
In addition, the subsequent pressure history at the valve could be used to infer the steady-flow pressure gradient along the pipe and hence to assess the flow resistance characteristics, as described in international patent publication number WO 02/25062 [Ref #4]. For a known pipe area, this would be a measure of the surface roughness condition.
However, the methods of WO 96/15427 and WO 02/25062 rely on the rapid closure of a pipeline valve to cause a change in the flowrate from an initial flowing condition to a zero flow condition. The initial flow conditions are inherently noisy, and therefore the collected data has a high level of background noise. The methods of WO 96/15427 and WO 02/25062 also rely on high flowrates, which exacerbate the ambient noise levels. The closure of the valve (which is a valve forming a part of the pipeline system) generates additional noise, and therefore the pressure measurement signal has a poor signal to noise ratio, limiting its usefulness in pipeline assessment applications. Furthermore, the process of closing the valve is difficult to perform quickly, and therefore the resulting pressure surge is generated slowly. This limits the use of the resulting pressure transient for the detection of non-uniformities or gradual changes in the pipeline condition.
In some offshore applications, the operator wishes to know how the diameter/area of a pipe varies along its length. In some cases, this can be achieved by direct geometrical measurement or by the use of an intelligent pig, but these methods are not always suitable, particularly where there is a risk of the pig becoming stuck at a restriction or blockage in the pipeline.
It is amongst the aims and objects of the invention to provide a method of and apparatus for assessing and/or condition monitoring of pipelines which avoid the drawbacks and disadvantages of the prior art.
It is an aim of at least one aspect of the invention to provide a method of using pressure surges to infer area profiles of fluid conduits and a corresponding apparatus.
Other aims and objects will become apparent from reading the following description.
According to first aspect of the invention, there is provided a method of assessing the condition of a fluid conduit, the method comprising:
providing a measurement apparatus comprising a measurement region in fluid communication with the fluid conduit, the measurement apparatus operable to measure pressure and flowrate in the measurement region;
effecting or generating a pressure surge to propagate in the fluid conduit;
measuring a pressure response from the fluid conduit to provide a pressure data set;
measuring a flowrate response from the fluid conduit to provide a flowrate data set;
determining, from the pressure data set and the flowrate data set, a condition data set representative of the cross-sectional area along the fluid conduit.
In the context of the description, the term “the condition of a fluid conduit” is used to generally to refer to the overall state of the flow path defined by a fluid conduit (primarily the effective flow area of the conduit), including its internal condition which may be affected for example by build-up and deposition of scale, sand, waxes and other materials on the interior surface. The term should not be construed as relating exclusively to the inherent physical condition of the conduit or its components (although in some cases there may be a link between the overall condition of the conduit and its inherent physical condition).
The term “pressure surge” is used for convenience and will be understood to refer to a change in pressure conditions in the fluid conduit. Such a change in conditions will of course also change fluid velocity and flowrate. References to pressure surge do not necessarily imply generation of a surge by means of a controlled pressure change; a controlled flowrate change could also effect a “pressure surge”, and indeed the terms “velocity surge” or “flowrate surge” could be used to describe the same effect.
In the present invention, the area is inherently unknown and is inferred from the pressure data set and the flowrate data set. Preferably, the cross-sectional area is determined using the equation:
Q=Q
0
±A·(P−P0)/ρc
where Q and Q0 are flowrates; A is the cross-sectional area of the fluid conduit; P and P0 are pressures; ρ is fluid density; and c is the speed of sound in the fluid.
The measurement region may be a fluid conduit portion of the measurement apparatus, which may be continuous with the fluid conduit of which the condition is to be assessed.
Preferably, the method comprises determining data representative of the cross-sectional area at one or more locations along the fluid conduit. The method may comprise determining data representative of average cross-sectional area of the fluid conduit in a region or section of the fluid conduit.
In preferred embodiments of the invention, one of the pressure or flowrate in the measurement region is controlled, at least during a measurement interval, to render it approximately constant before the pressure surge is effected or generated, and approximately constant after the pressure surge is effected or generated. Preferably, the one of the pressure or flowrate is held substantially constant before and after the pressure surge is generated. By maintaining one of the pressure or flowrate in the measurement region approximately or substantially constant, measurement inaccuracies in the method may be reduced. However, in methods according to alternative embodiments of the invention, slowly varying pressure or flowrate (i.e. variations to pressure or flowrate which are significantly slower than an intentional step-like change) conditions may be used with the contemplated measurement and analysis steps and without a significant effect on the results obtained.
Preferably the method comprises opening a valve to generate a pressure surge in the fluid conduit. Preferably the valve forms a part of the measurement apparatus, and is selected or designed for rapid actuation with low noise generation. This is in contrast to prior art approaches to the application of pressure surges in pipeline systems which use valves that form a part of the pipeline system and which are not designed for rapid or low noise actuation. These pipeline valves are therefore not capable of generating steep fronted pressure surges. The rapid opening of the valve provides improved differentiation between friction and area variations and also enables the cross-sectional area versus distance along the fluid conduit to be analysed in greater detail. This is important for the assessment of gradual changes in cross-sectional area, as may result from the build-up and deposition of scale, sand, waxes and other materials on the interior surface of a pipeline.
The method may comprise causing a flow condition in the measurement region to change from a first flow condition in which there is negligible flow of fluid in the measurement region to a second flow condition in which fluid flows in the measurement region and the fluid conduit, to thereby effect or generate a pressure surge to propagate in the fluid conduit. The method may therefore comprise creating a flow of fluid in the fluid conduit. Preferably the flowrate in the fluid conduit prior to the generation of a pressure surge in the fluid conduit is negligible, and more preferably is zero. Thus the initial flow conditions may provide a relatively low level of background noise, and the collected data may have a good signal to noise ratio. This is a significant advantage over prior art approaches to the use of pressure surges in pipeline applications which use an initial flow condition with a correspondingly high level of background noise. This approach also reduces the influence of friction on the measurements.
Preferably, the method comprises:
effecting or generating a pressure surge in the fluid conduit at a first upstream location or region;
measuring the pressure response from the fluid conduit;
measuring the flowrate response from the fluid conduit;
determining, from a pressure data set and a flowrate data set, a data set representative of the cross-sectional area along the fluid conduit downstream from the first upstream location or region.
Preferably, the pressure response and/or flowrate response are sampled over a measurement interval. This may begin a short time after the step of effecting or generating a pressure surge where the amplitude of the pressure surge is known, but preferably begins before the step of effecting or generating a pressure surge. This has the benefit of identifying the amplitude and character of the pressure surge. The pressure surge is caused to propagate in the fluid conduit. Preferably, the measurement interval finishes before the pressure surge reaches a blockage or complete restriction (such as a fluid conduit termination or a valve) in the fluid conduit.
Preferably, the method comprises sampling the flowrate with high accuracy at a rate greater than 100 Hz, and preferably greater than 1 kHz. The interpretation of the measured data is preferably based in the time domain.
Preferably the fluid conduit is a pipeline flooded with a liquid. Preferably, the method comprises the step of flooding the fluid conduit with a liquid. The method may be performed on a pressurised fluid conduit, and the method may comprise pressurising the fluid conduit to a predetermined pressure. The predetermined pressure may be selected to be sufficient to reduce a volume of gas in the fluid conduit.
The method may include separating the effects of skin friction and local area changes.
Preferably, the method comprises:
effecting or generating a first pressure surge to propagate in the fluid conduit, the first pressure surge having a first amplitude in the fluid conduit;
measuring at least one of a pressure response or a flowrate response from the fluid conduit to provide a first pressure response or a first flowrate response data set;
effecting or generating a second pressure surge to propagate in the fluid conduit, the second pressure surge having a second amplitude in the fluid conduit; and
measuring at least one of a pressure response or a flowrate response from the fluid conduit to provide a second pressure response or a second flowrate response data set.
The method may comprise comparing the first and second data sets to determine a data set representative of the cross sectional area along the fluid conduit.
Preferably the second amplitude is different from the first amplitude.
The method may comprise repeating effecting or generating a plurality or series of pressure surges at a range of different amplitudes, and collecting at least one of a measurement of pressure or flowrate for the plurality of pressure surges.
The method may comprise estimating or determining the cross-sectional area of the fluid conduit in a first interval close to the first upstream location or region, and may comprise estimating or determining the wavespeed of the pressure surge in the fluid conduit at an area close to the first upstream location or region. The first interval may be a first section of the fluid conduit to be assessed.
The method may comprise dividing the fluid conduit into sections, and determining cross sectional area of the sections sequentially.
The method may comprise effecting or generating a pressure surge to propagate in the fluid conduit at a second time, later than the first time;
measuring at least one of a pressure response or a flowrate response from the fluid conduit to a second set of fluid conduit condition data;
comparing the second set of fluid conduit condition data with the first set of fluid conduit condition data.
According to a second aspect of the invention, there is provided a method of monitoring the condition of a fluid conduit, the method comprising:
providing a measurement apparatus comprising a measurement region in fluid communication with the fluid conduit, the measurement apparatus operable to measure at least one of pressure and flowrate in the measurement region;
effecting or generating a pressure surge to propagate in the fluid conduit at a first time;
measuring at least one of a pressure response or a flowrate response from the fluid conduit to provide a first set of fluid conduit condition data;
effecting or generating a pressure surge to propagate in the fluid conduit at a second time, later than the first time;
measuring at least one of a pressure response or a flowrate response from the fluid conduit to a second set of fluid conduit condition data;
comparing the second set of fluid conduit condition data with the first set of fluid conduit condition data.
The method may further comprise effecting or generating a pressure surge in the fluid conduit at a third and/or further time, later than the second time, to determine a third and/or further set of fluid conduit condition data; and comparing the third and/or further set of fluid conduit condition data with the second and/or first set (or other previous set) of fluid conduit condition data.
The method may comprise the step of comparing the second set of fluid conduit condition data with the first set of fluid conduit condition data to determine a change in condition of the fluid conduit. The change of condition may be a change in cross-sectional area of the fluid conduit. The times between acquisition of second and subsequent data sets may be of the order of weeks or months in the context of the monitoring applications envisaged for the present invention.
The method may comprise outputting a signal indicative of a change in condition of the fluid conduit, and may comprise providing a graphical representation of a change in the condition. Alternatively or in addition, the method may comprise providing a graphical representation of the first set of fluid conduit condition data and/or the second or subsequent sets of fluid conduit condition data.
In a preferred embodiment of the invention, the method comprises measuring a pressure response to provide a pressure data set and measuring a flowrate response to provide a flowrate data set. The measured pressure data set and the flowrate data set may then be used to determine the data set representative of the cross-sectional area along the fluid conduit.
In alternative embodiments, one of the pressure response or flowrate response is measured, with the other of the pressure or flowrate in the measurement region being a prescribed or known condition. In other words, one of the pressure data set and flowrate data set may be known a priori, and may be used together with the measured data set to determine the data set representative of the cross sectional area along the fluid conduit. However, due to practical difficulties with prescribing one of the pressure or flowrate in the measurement region with the desired accuracy, it is preferable that both the pressure response and the flowrate response are measured.
Embodiments of the second aspect of the invention may include one or more features of the first aspect of the invention or its embodiments, or vice versa.
According to a third aspect of the invention, there is provided a method of collecting pipeline condition data, the method comprising:
effecting or generating a pressure surge to propagate in a fluid conduit;
measuring a pressure response from the fluid conduit to provide a pressure data set;
measuring a flowrate response from the fluid conduit to provide a flowrate data set.
Embodiments of the third aspect of the invention may include one or more features of the first or second aspects of the invention or their embodiments, or vice versa.
According to a fourth aspect of the invention there is a method of processing pipeline condition monitoring data, the method comprising:
providing a pressure data set and a flowrate data set acquired by effecting or generating a pressure surge to propagate in a fluid conduit and measuring a pressure and flowrate responses from the fluid conduit;
determining from the pressure data set and the flowrate data set, a data set representative of the cross sectional area along the fluid conduit.
In the present invention, the area is inherently unknown and is inferred from the pressure data set and the flowrate data set. Preferably, the cross-sectional area is determined using the equation:
Q=Q
0
±A·(P−P0)/ρc
where Q and Q0 are flowrates; A is the cross-sectional area of the fluid conduit; P and P0 are pressures; ρ is fluid density; and c is the speed of sound in the fluid.
Embodiments of the fourth aspect of the invention may include one or more features of any of the first to third aspects of the invention or their embodiments, or vice versa.
According to a fifth aspect of the invention, there is provided a method of assessing the condition of a pipeline in an oil or gas exploration, production or transportation system, the method comprising:
effecting or generating a pressure surge to propagate in the pipeline at a first tie in point;
measuring a pressure response from the pipeline in a measurement region at or near the tie in point to provide a pressure data set;
measuring a flowrate response from in the pipeline in the measurement region to provide a flowrate data set;
determining, from the pressure data set and the flowrate data set, a condition data set representative of a cross sectional area along the pipeline downstream of the measurement region from the measured pressure and measured flowrate.
Preferably the method comprises detecting the presence and/or location of a partial restriction to the pipeline. Preferably the partial restriction is a distributed restriction on the interior of the pipeline.
Embodiments of the fifth aspect of the invention may include one or more features of any of the first to fourth aspects of the invention or their embodiments, or vice versa.
According a sixth aspect of the invention, there is provided a method of monitoring the condition of a pipeline in an oil or gas exploration, production or transportation system, the method comprising:
effecting or generating a first pressure surge to propagate in the pipeline at a first tie in point at a first time;
measuring a pressure response and a flowrate response from the pipeline in a measurement region at or near the tie in point during a first measurement interval to determine a first set of pipeline condition data;
effecting or generating a second pressure surge in the pipeline at a second time, later than the first time;
measuring a pressure response and a flowrate response from the pipeline at or near the tie in point during a second measurement interval to determine a second set of pipeline condition data;
comparing the second set of pipeline condition data with the first set of pipeline condition data.
Embodiments of the sixth aspect of the invention may include one or more features of any of the first to fifth aspects of the invention or their embodiments, or vice versa.
According to a seventh aspect of the invention there is provided a fluid conduit condition assessment apparatus comprising:
a connection for coupling the apparatus to a tie in point of a fluid conduit;
a pressure source;
a controllable actuator for generating a pressure surge in a fluid conduit in use;
a pressure sensor;
and a flowrate sensor.
Preferably the apparatus comprises a measurement region, which may be a fluid conduit portion of the measurement apparatus, and which may be continuous with the fluid conduit of which the condition is to be assessed.
The apparatus may comprise a valve operable to generate or effect the pressure surge, and the controllable actuator may be configured to operate the valve. Preferably the valve is a fast-acting valve, so that the pressure surge approximates a step change in pressure with a single wavefront.
The valve may be configured to isolate the fluid conduit from a controlled pressure source in a first (closed) condition. The controlled pressure source may comprise an accumulator. Alternatively the valve may be configured to isolate the fluid conduit from a controlled flow source in a first (closed) condition.
In one embodiment, the apparatus comprises a pressure cylinder and a piston which is actuable to deliver a controlled flowrate to the fluid conduit.
The apparatus may be configured to control one of the pressure or flowrate in the measurement region, at least during a measurement interval, to render it approximately constant before the pressure surge is effected or generated, and approximately constant after the pressure surge is effected or generated. Preferably, the one of the pressure or flowrate is held substantially constant before and after the pressure surge is generated. The apparatus may therefore be configured to suddenly generate a step change in pressure to a nominally constant value, and/or may be configured to suddenly generate a step change in flowrate to a nominally constant value. By maintaining one of the pressure or flowrate in the measurement region approximately or substantially constant, measurement inaccuracies in the method may be reduced. Slowly varying pressure or flowrate (i.e. variations to pressure or flowrate which are significantly slower than an intentional step-like change) conditions may be used in methods according to alternative embodiments of the invention.
Preferably the apparatus is configured to perform the methods of any of the first to sixth aspects of the invention, and embodiments of the seventh aspect of the invention may include one or more features of any of the first to sixth aspects of the invention or their embodiments, or vice versa.
According to an eighth aspect of the invention there is provided a pipeline condition assessment system comprising:
a fluid conduit;
an apparatus coupled to a tie in point of the fluid conduit;
a pressure source;
a controllable actuator for generating a pressure surge in the fluid conduit;
a pressure sensor;
and a flowrate sensor.
Preferably the apparatus is configured to perform the methods of any of the first to sixth aspects of the invention, and embodiments of the eighth aspect of the invention may include one or more features of any of the first to seventh aspects of the invention or their embodiments, or vice versa.
According to a ninth aspect of the invention, there is provided a method of assessing the condition of a fluid conduit, the method comprising:
providing a measurement apparatus comprising a measurement region in fluid communication with the fluid conduit, the measurement apparatus operable to measure at least one of pressure and flowrate in the measurement region;
effecting or generating a pressure surge to propagate in the fluid conduit by actuation of a valve in the measurement apparatus;
measuring at least one of a pressure response or a flowrate response from the fluid conduit;
determining, from a pressure data set and a flowrate data set, a data set representative of the cross sectional area along the fluid conduit.
Preferably the valve is a fast-acting valve, so that the pressure surge approximates a step change in pressure with a single wavefront. By generating a pressure surge using a valve in the apparatus, the uncertainties and disadvantages of using a valve which forms a part of the fluid conduit system are avoided. In particular, the valve may be selected to be fast-actuating without generating significant noise levels on opening, to maintain a sufficiently high signal to noise ratio.
Preferably the method comprises causing the flow condition in the measurement region to change from a first flow condition in which there is negligible flow of fluid in the measurement region to a second flow condition in which fluid flows in the measurement region and the fluid conduit, to thereby effect or generate a pressure surge to propagate in the fluid conduit. Thus the initial flow conditions may provide a relatively low level of background noise, and the collected data may have a good signal to noise ratio.
Embodiments of the ninth aspect of the invention may include one or more features of any of the first to eighth aspects of the invention or their embodiments, or vice versa.
According to a tenth aspect of the invention, there is provided a method of assessing the condition of a fluid conduit, the method comprising:
providing a measurement apparatus comprising a measurement region in fluid communication with the fluid conduit, the measurement apparatus operable to measure at least one of pressure and flowrate in the measurement region;
causing the flow condition in the measurement region to change from a first flow condition in which there is negligible flow of fluid in the measurement region to a second flow condition in which fluid flows in the measurement region and the fluid conduit, to thereby effect or generate a pressure surge to propagate in the fluid conduit;
measuring at least one of a pressure response or a flowrate response from the fluid conduit;
determining, from a pressure data set and a flowrate data set, a data set representative of the cross sectional area along the fluid conduit.
Preferably, the first flow condition is substantially zero or zero flow. Thus the initial flow conditions may provide a relatively low level of background noise, and the collected data may have a good signal to noise ratio.
The method may comprise effecting or generating a pressure surge to propagate in the fluid conduit by actuation of a valve in the measurement apparatus.
Embodiments of the tenth aspect of the invention may include one or more features of any of the first to ninth aspects of the invention or their embodiments, or vice versa.
The methods of the various aspects of the invention and/or the critical steps thereof are preferably implemented in software, although it will be understood that the methods or steps thereof may also be implemented in firmware or hardware or in combinations of software, firmware or hardware.
There will now be described, by way of example only, various embodiments of the invention with reference to the drawings, of which:
The basic principles of the invention will be described with reference to
The apparatus comprises a fluid conduit portion 22 which has an outlet 24 coupled to the pipeline 12 via the tie-in point 14. The opposing end of the conduit portion 22 is coupled to a fluid source, which in this example is a liquid vessel 25 connected to an accumulator 26 by a large diameter hose 27. The accumulator 26 is a large air accumulator controlled by a regulator 29, connected to a high pressure air source 31. The conduit portion 22 defines a measurement region for parameters of interest in the inventive method, and includes an accurate flowrate sensor or flowmeter 28, and a pressure sensor 30. Sensors 28 and 30 output measurements to a control module 32, which includes data storage and/or processing means (not shown). The conduit portion 22 also comprises a controllably actuable valve 34, also controlled by the control module 32. The valve 34 is a fast operating valve such as a hydraulically actuated pilot valve.
In use, the apparatus 10 is coupled to the pipeline 12 and the isolation valve 16 is opened. The pipeline 12 is flooded with liquid until it is full and a steady state is reached. The valve 34 is closed and the flowrate in the conduit portion 22 and pipeline 12 is zero. In this zero flowrate condition the ambient noise levels are relatively low compared with a system in which fluid is flowing in the pipeline. The pressure in the accumulator 26 is increased such that there is a pressure differential across the valve 34 of a predetermined amplitude. The valve 34 is controlled to be opened to generate a pressure surge in the conduit and pipeline 12, which propagates through the pipeline system. The valve 34 opens quickly to generate a steep-fronted pressure surge in the system. The valve 34 is designed such that its rapid opening does not generate significant levels of noise, and as such the sensors 28 and 30 are able to collect relatively clean signals. This is an important advantage over previous attempts to analyse data collected from pressure surges created by slamming shut pipeline valves.
Pressure data and flowrate data are sampled using sensors 28, 30 with high accuracy (for example at a frequency of 1 kHz or greater), as the wavefront of the pressure surge propagates in the pipeline 12. The velocity at which the wavefront propagates in the pipeline (i.e. the wavespeed) is known from theoretical calculation or measurement. In some cases, where a good estimate of the cross-sectional area of the pipeline in the region of the pipeline which is adjacent to the apparatus, the wavespeed can be estimated using the equation Q=Q0±A·(P−P0)/ρc; where Q and Q0 are flowrates; A is the cross-sectional area of the fluid conduit; P and P0 are pressures; ρ is fluid density; and c is the wavespeed. Techniques which use the transit time of a pulse propagating in the fluid conduit and reflected from a known position (such as a closed valve) in a conduit may also be used to calculate the wavespeed, and data relating to material properties or the pipeline may be used in wavespeed calculations or estimations.
In the example of
When the wavefront has propagated sufficiently far into the region 12b, it passes through the distributed restriction 40. The local area of flow in region 12b will be significantly smaller than that at the control end and in the first portion 12a. As a consequence of the reduction, the wavefront will behave in a manner that lies between two extreme conditions. At one extreme, the rate of flow remains constant; in this case, the velocity at the wavefront must increase and, hence, so must the pressure. At the opposite extreme, imagine that the pressure remains constant; in this case, the velocity must also be constant and so the rate of flow must have decreased. In reality, neither extreme occurs; instead, the pressure and velocity at the wavefront increase, and the flowrate decreases. The sensors 28, 30 are in the apparatus 10 at the control end of the pipe so the above behaviour will not be observed directly. However, the continuously changing conditions at the wavefront will cause continuous wave reflections that will be detected at the control sensors 28, 30.
Secondary and additional changes in flowrate are detected by the apparatus at later times (from t=0 onwards). These changes are due to interaction with the restriction by secondary and additional wavefronts 42d, 42e, reflected from the end of the pipeline (at the measurement region) back towards the restriction, or wavefront 42f, resulting from a secondary interaction of reflected wavefront 42c with the restriction in the pipeline.
In this embodiment, approximately constant control pressure is prescribed, and therefore the changes to the wavefront pressure are not in fact detected by the pressure sensor 30, as the accumulator adjusts for the change in pressure by regulating the flowrate. However, the changes will be detected in the flowrate measurements: the reduction in flowrate due to the distributed restriction will be detected as a change in the flowrate measurement data; and the increase in pressure will be detected as an additional reduction in the flowrate measurement data as the accumulator 26 adjusts the flowrate to compensate for the pressure increase in the pipeline. The cumulative change in the flowrate is therefore attributable in part to a flowrate change and in part to a pressure change which result from the restriction in the pipeline. In practice, the presence of restriction can be inferred from the decrease in the measured flowrate, and its position can be estimated from the times at which the changes are detected (the speed of the wavefronts in the fluid being known).
In this example, the pressure is prescribed at the measurement region by the accumulator 26. The flow conditions are initially at a steady state, with zero flowrate into the pipeline o and constant pressure (i.e. zero change in pressure, drawn as P=0 in the drawing). Shortly after time t=0, a valve is opened to create an increase in pressure (to 1.0 MPa), which is maintained by the accumulator 26 throughout the measurement interval (shown up to 10 seconds here). The pressure increase causes a flowrate increase to approximately 6.0 litres/second, and the wavefront propagates in the pipeline. At time t1, the wavefront meets a sudden restriction in the pipeline. At this moment, the pressure at the wavefront increases and the flowrate decreases. These changes are detected at the measurement region at time t2, as a decrease in flowrate, with the corresponding change in pressure being compensated for by the accumulator. At a time t4, the detectors measure an increase in flowrate, due to the wavefront encountering the end of the restriction. Again, in practice, the pressure decreases at time t4, but this pressure change is compensated for by the accumulator and is instead seen as a corresponding change in the flowrate. At later times t6 and t8, changes to the flowrate are detected due to the interaction of secondary and additional wavefronts with the restriction (as described with reference to
If further simulations were undertaken with successively shorter restrictions, the reflected pulses would become successively shorter. While this gives rise to practical limitations of the time-domain interpretation method of the invention in drawing reliable inferences about any short-duration features in the reflected signal, there are consequential benefits. For example, because of this characteristic, time-domain methods are potentially able to provide the desired information about extended regions of reduced area even if these regions include one or more local restrictions, and even when such local restrictions cause quite large reductions in area (perhaps as much as 50%). The existence of the local features will not, in itself, detract significantly from the potential of the methods of the invention to achieve the primary objective of assessing the condition of the pipeline over longer pipeline distances. For example, the method may be used to effectively deduce the cross-section of a pipeline in sections beyond the location of a valve in the pipeline.
After completing the assessment of the first section, the analysis is therefore repeated (in slightly modified form) in the next section, namely ΔL1<x<ΔL2. Initially this step will be undertaken on the assumption that the data already inferred in the first section is correct. This is an important choice because the section 0<x<ΔL1 will act as a telescope through which the section ΔL1<x<ΔL2 is viewed.
This process will then be repeated for successively more distant sections until either (a) the whole pipe has been assessed or (b) the inferred sensitivity to plausible alternative pipe profiles becomes too small for conclusions to be drawn with reasonable confidence.
There will now be described, with reference to
In the method of
When the test procedure is ready to begin, the isolation valve 16 is opened and the pipeline 12 is flooded with liquid until it is full and a steady state is reached. At this point, the flow into the pipeline is negligible and is preferably nil. The pressure in the accumulator is regulated to provide a pressure differential across the valve 34 according to the surge amplitude predetermined for the test procedure. The valve 34 is then rapidly opened to create a pressure surge which propagates in the apparatus conduit 22 and pipeline 12 (step 63). The valve 34 is preferably a fast acting hydraulic valve, as it is desirable for the pressure surge to be as close as possible to a square step change in pressures (as shown in
The sensors 28, 30 rapidly sample flowrate and pressure data respectively (step 64), and time series of flowrate and pressure data are recorded by the apparatus (step 65). Typical sampling rates for the flowrate and/or pressure are of the order of tens of Hz, and preferably greater than 100 Hz. It is desirable for sample rates to exceed 1 kHz to provide improved data resolution, particularly in shorter pipeline measurement lengths and/or pipeline systems with higher wave speeds. The data acquisition begins prior to the generation of the surge, so that the time series includes the detection of the wavefront as it passes through the apparatus conduit 22. Data are recorded to digital memory in the apparatus for subsequent analysis. The data may be retrieved by removing digital storage from the apparatus, or by interfacing the apparatus with computer processing equipment or an external storage device. It will also be appreciated that in alternative embodiments, a processing module may be incorporated in apparatus to carry out some or all of the required analysis or interpretation.
The time series flowrate and pressure data is recorded over a measurement interval which enables all useful information relating to the propagation of the wavefront in the fluid conduit to be collected. In practice the useful data measurement interval may finish before the time taken for the wavefront to reach (and be reflected from) a major blockage or complete restriction (such as a fluid conduit termination or closed valve) in the fluid conduit. After this point, the acquired data is dominated by the complete reflection of the pressure surge by the blockage or complete restriction, and little or no useful information about the condition of the pipeline along the length of the pipeline can be derived.
Optionally, the test procedure is repeated at the same pressure amplitude (for example, to provide data redundancy and improved signal to noise ratio) or at a different pressure surge amplitude (the advantages of which will be described later with reference to
When all the relevant data has been collected, it is analysed (step 70) to deduce information about the cross sectional area of the pipeline along its length. The analysis and interpretation of data are described in greater detail with reference to
This embodiment of the invention utilises pressure surge analysis software which is widely used in industry and academic research to model the flow regimes in a range of flow systems. Examples of commercially-available software packages include software available from Flight Software N.I. Limited under the trade mark FluidFlow3; the pressure surge analysis software available from Flowmaster Group and designated by the commercial name FLOWMASTER; the flow assurance and surge analysis software available from Simulation Software Limited and sold under the trade mark VARISIM; and the water hammer analysis software available from KYPipe LLC sold under the product designation SURGE. These software packages (and others) enable a pipeline system to be designed and boundary conditions related to flowrates and pressures to be configured by the user. The software then models the flow regime enabling flowrates and pressures in the modelled pipeline system to be simulated according to the boundary conditions input by the user.
As a first step, data acquired from a previous test procedure is input into the software (step 71). The user also inputs pipeline system data to the software to model the flow system (step 73). The data is input according to the features of the system which are known, which may be approximate topologies and/or pipeline components. In the context of the invention, the true cross-sectional area of the pipeline is unknown, and is the parameter to be deduced, although a reasonable starting value for the analysis may be derived from the original pipeline system design data.
As noted above, to simplify the analysis, cross-sectional areas of sections of the pipeline will be deduced sequentially, section-by-section, starting from a first section closest to the measurement apparatus and working away from the measurement apparatus. The cross-sectional areas attributed to each section will be used to deduce the areas of sequential sections. Therefore as a further preliminary step it is necessary to divide the pipeline into a number of sections (step 72). This may be performed by any of a range of techniques. One automated technique is to calculate the average flowrate over two halves of a pipeline interval extending from a starting point x0 in the pipeline to a distance x1. The distance x1 is extended until the average flowrates in the first and second halves differ by more than a predetermined threshold (say 10%). This indicates a change in pipeline conditions over the length to x1 which is significant enough to take the pipeline section x0 to x1 as the first section for analysis. Other simple techniques include identifying significant changes in the approximate gradient of flowrate with respect to distance along the pipeline, and/or selecting sections in which the gradient is reasonably constant.
With the first pipeline section selected, a starting value for cross-sectional area A is chosen (step 73). The software package is used to calculate, based on its input parameters (including the starting value for A) and the measured pressure data time series, a simulated flowrate data time series (step 74). This simulated flowrate data time series is then compared with the measured flowrate data time series (step 75). Discrepancies between the simulated and measured data indicate one or more inaccuracies or false assumptions in the starting data and/or errors in the measurement signal. The method therefore applies statistical techniques to fit the simulated data to the measured data by optimising the cross-sectional area value A (step 80). This is achieved by making an adjustment to the cross-sectional area value A in the pipeline section being assessed, and then recalculating the simulated data set based on the new value for A. The new simulated data set is compared to the measured data and the process is repeated until the simulated data set is matched to the measured data set to an acceptable degree. The value of A used in the matched simulated data set is taken to be the cross-sectional area for the pipeline in the section analysed.
In either case, the profile of the cross-sectional area is optimised to provide the best fit for the simulated data with the measured data to deduce cross-sectional area information in the pipeline section.
A variety of statistical analysis techniques may be used within the scope of the invention, as would be appreciated by one skilled in the art. Least-squares analysis methods are one example, but equally other techniques may be used. The statistical method will include appropriate allowance for residual noise in the measured signals. This may include the time-averaging of signals to create smoothing and it may include the specification of tolerances with which differences between measured and trial-predictions will be compared.
Returning now to
A pressure wave moving along a pipe is approximately equal the sum of the speed of sound and the fluid speed, U. Thus, relative to the pipe, waves travelling downstream propagate at c+U and waves travelling upstream move slightly more slowly at c−U. The propagating waves are influenced by (a) inertia, (b) variations in area and (c) friction, all of which are experienced simultaneously, and can be described by the following equations:
(satisfying a wave travelling downstream); and
(satisfying a wave travelling upstream).
In both equations: a is cross-sectional area; c is speed of sound; l is the cross-sectional perimeter; p is pressure; r is the time co-ordinate; U is fluid velocity; x is the distance co-ordinate; ρ is fluid density; and τw is shear stress at the pipe wall surface. The three terms on the right-hand side of the equations are measures of the influences of (i) inertia, (ii) variations in area and (iii) friction respectively.
It is intended that the testing will be undertaken with pressure surges that do not cause the measurements to be strongly influenced by friction. Accordingly, the simplified method described above ignores the effect of friction for the purposes of determining the cross-sectional area. However, this might not always be possible. When it is not, the accuracy with which the pipeline geometry can be inferred will depend upon the extent to which the influences of friction and area change can be separated. This will be described below with reference to
In the method 90, data is acquired during a test procedure, for example according to the method 60 of
As with the method of
The user also inputs pipeline system data to the software to model the flow system (step 94). As before, the data is input according to the features of the system which are known, which may be approximate topologies and/or pipeline components. Also as before, the true cross-sectional of the pipeline is unknown, and is the parameter to be deduced, although a reasonable starting value for the analysis may be assumed from the original pipeline system design data. A series of m values of A are chosen for the system input data, and a series of n values of f, representing the frictional interaction between the fluid and the pipe (step 95). The series of A and f each include a range of values which has both smaller and larger values than the most likely values of A and f. The system data input into the system are of course dependent on the pressure surge amplitude, and therefore the software is required to model m×n combinations of area and friction parameters for each pressure surge value used for the pipeline section under analysis.
The pressure surge analysis software is used to calculate (step 96), based on its input parameters (including m×n combinations of the starting values for A and f) and the measured pressure data time series, a simulated flowrate data time series. This is performed for each data set in turn, so that k×m×n simulated flowrate time series are calculated.
The simulated flowrate data time series are then compared (step 97) with the measured flowrate data time series. Discrepancies between the simulated and measured data indicate one or more false assumptions in the starting data (i.e. the choice of A and f) or measured data, and statistical techniques (for example, the mean square difference) are used to derive a statistical value (step 98) and rank the m×n combinations of A and f. An average of the statistical values for each combination of A and f is calculated (i.e. averaged across all k data sets corresponding to pressure surge values) to give an m×n average values (step 99). The m×n average values are then reviewed to select the one that implies the closest match to the measured flowrate time series (step 100). The values of A and f in that gives this selected average value is deemed to be the correct combination of A and f for the pipeline section under analysis.
For the assessment of the next section of the pipeline, the software takes the previously deduced values of A and f as fixed input parameters in the pipeline system data and the method is repeated for the next section (step 101).
The method described with reference to
The above-described embodiments relate to a methodology in which the pressure in the measurement region (or at the tie-in point of the fluid conduit) is the prescribed condition, and is kept approximately constant in the measurement interval. However, it is also within the scope of the invention to perform a similar method in which the flowrate in the measurement region (or at the tie-in point of the fluid conduit) is the prescribed condition, and is kept approximately constant. It is also within the scope of the invention to perform a similar method in which one of the pressure or flowrate is slowly varying or changing (i.e. variations to pressure or flowrate are significantly slower than intentional step-like changes).
The apparatus 210 comprises a pressure cylinder 220 which accommodates a piston 222. The internal volume 223 of the pressure cylinder 220 is in fluid communication with the pipeline, and comprises an inlet 230 to enable the inflow of fluid from a fluid source (not shown). The piston is connected to an actuator 224, which is a threaded rod driven by an electric motor 226 via a mechanical gearbox 228, such that rotation of the electric motor 226 causes the piston to be extended or retracted in the cylinder to reduce or increase the internal volume 223. The gearbox 228 is configured to drive the piston into the cylinder at a preset rate which would not vary during operation. An electric motor is preferred as it rotates at a constant speed, drawing more or less current as the torque load increases, whereas a hydraulically actuated system may vary the rate of progression of the piston as the pressure changes. The apparatus 210 also comprises a control module 232 which controls the electric motor 226, flowmeter 233 and pressure sensor 234 in the measurement region defined by the cylinder during a measurement interval.
In use, the cylinder 220 is filled with a liquid via the inlet 230 with the piston fully retracted to the start position to maximise the volume 223. The isolation valve 216 is opened and the pressure in the cylinder is equalised with the pressure in the pipeline such that conditions are stable. The electric motor drives the piston into the cylinder to displace fluid from the apparatus and into the pipeline at a constant flowrate for a single stroke. The apparatus therefore generates a pressure surge in the pipeline 212. Flowrate and pressure data are measured by flowmeter 233 and pressure sensor 234 in the measurement region. The data acquisition begins prior to the generation of the surge, so that the time series includes the detection of the wavefront as it passes through the measurement region. Data are recorded to digital memory in the apparatus for subsequent analysis.
The data is analogous to the example described with reference to
Thus according to the method described above, pressure and flowrate histories can be derived from a system in which the flowrate is kept approximately constant or slowly varying at one end of the pipeline system. Such an approach may not be preferred in all practical implementations due to the difficulties in achieving a sharply defined incident wavefront (or step-form pressure surge). However, there are also advantages to such an approach, including those associated with the ease of rapid pressure sampling (compared with the flowrate sampling relied upon in the previous embodiments). Data acquired in this way may be analysed according to the embodiments described with reference to
The apparatus 310 differs from the apparatus 210 in that it does not include a pressure cylinder 220 accommodating a piston. Instead, the apparatus 310 has a fluid or slowly varying source 326 which in a first condition, delivers fluid at a controlled and substantially constant or slowly varying flowrate to an inlet conduit portion 327a of the apparatus. A valve 329, which in the example shown is a flapper valve, isolates the inlet conduit portion 327a from the fluid conduit portion 322, and fluid flows into an outlet 327b (in this case a return line is shown to the fluid source 326 although other configurations are possible).
In use, fluid flows at a controlled flowrate in the flow path defined by the inlet conduit portion 327a and the outlet 327b. The isolation valve 316 is opened and the pressure in the cylinder is equalised with the pressure in the pipeline such that conditions are stable. The valve 329 is actuated to move to a position, shown by the dotted line 331, in which it diverts flow from the inlet conduit portion 327a to the conduit portion 322. The apparatus therefore generates a pressure surge in the pipeline 312. Flowrate and pressure data are measured by flowmeter 333 and pressure sensor 334 in the measurement region. The data acquisition begins prior to the generation of the surge, so that the time series includes the detection of the wavefront as it passes through the measurement region. Data are recorded to digital memory in the apparatus for subsequent analysis.
The system 301 has the advantage that the flowrate in the measurement region can be changed from zero to a non-zero flow condition quickly by actuation of the valve, and may therefore produce a well-defined pressure surge in the pipeline, which is useful for detecting gradual changes in the pipeline condition along its length.
A further alternative embodiment of the invention is shown schematically with reference to
The apparatus 410 comprises a pressure cylinder 420 in fluid communication with the pipeline, and comprises an inlet 430 to enable the inflow of fluid from a fluid source (not shown). The piston is connected to an actuator 424, which causes the piston to be extended or retracted in the cylinder to reduce or increase the internal volume 423. The apparatus 410 also comprises a control module 432 which controls the electric motor 426, flowmeter 433 and pressure sensor 434 in the measurement region defined by the cylinder during a measurement interval.
In use, the cylinder 420 is filled with a liquid via the inlet 430 with the piston fully retracted to the start position to maximise the volume 423. The isolation valve 416 is opened and the pressure in the volume of the cylinder 425a on the forward side of the piston 422 is equalised with the pressure in the pipeline such that conditions are stable. In this embodiment the pipeline 412 contained pressurised fluid, at (for example) 4,000 psi (about 27.6 MPa). The volume of the cylinder 425b on the reverse side of the piston 422 is fluidly connected to a pressure source such as a bank of accumulators 427 to balance the pressure across the piston (i.e. in this example to maintain it at 4,000 psi or about 27.6 MPa).
The electric motor drives the piston into the cylinder to displace fluid from the apparatus and into the pipeline at a constant flowrate for a single stroke. The apparatus therefore generates a pressure surge in the pipeline 412, over and above the baseline pressure of 4,000 psi. Flowrate and pressure data are measured by flowmeter 433 and pressure sensor 434 in the measurement region. The data acquisition begins prior to the generation of the surge, so that the time series includes the detection of the wavefront as it passes through the measurement region. Data are recorded to digital memory in the apparatus for subsequent analysis in accordance with the methods described above.
An advantage of the system of
It will be appreciated that the above-described embodiments of the invention may be varied within the scope of the invention. In particular, it will be apparent to one skilled in the art that the method is potentially well-suited to use as the basis of automatic assessment algorithms enabling area profiles to be inferred from the measured signals without detailed human intervention.
According to the invention, there is described a method of and apparatus for assessing the condition of a fluid conduit. The method comprises providing a measurement apparatus comprising a measurement region in fluid communication with the fluid conduit. The measurement apparatus is operable to measure at least one of pressure and flowrate in the measurement region. A pressure surge is effected or generated to propagate in the fluid conduit, and at least one of a pressure response or a flowrate response from the fluid conduit is measured. From a pressure data set and a flowrate data set, a data set representative of the cross sectional area along the fluid conduit is determined. This enables average or overall area profiles along long fluid pipelines to be deduced.
The theoretical basis of the method is a relatively simple application of the pressure surge analysis, and a preferred embodiment utilises pressure and flow measurements at one end of a pipeline during the introduction of a disturbance just upstream of the measurement location. The embodiments are generally described as being implemented in a computer system and it will be appreciated that the invention may be implemented in software, hardware, firmware, or a combination thereof.
It is an advantage of the method that interpretation is based in the time domain that dominates everyday life, not in the frequency domain (in contrast with many prior art methods of assessing measured signals that include rapid changes). Time-domain analysis is well suited to the deduction of overall area profiles even though frequency-domain analysis is usually necessary for detecting and locating local phenomena such as leaks and restrictions.
The method is able to infer cross-sectional areas close to the measurement locations with good accuracy. Allowance can be made for skin friction even though its effects can be closely analogous to those of area variation.
The invention provides measurements with a good signal to noise ratio, which enables small and gradual changes in the pipeline condition to be detected. In particular, the invention enables the use of small flowrates which enable data interpretation to distinguish between noise and intended changes to the flow system.
The method is relatively insensitive to short lengths of restricted area, which enables the technique to be used to infer pipeline information beyond restrictions such a partially closed valve.
Various modifications may be made within the scope of the invention as herein intended, and embodiments of the invention may include combinations of features other than those expressly claimed.
Number | Date | Country | Kind |
---|---|---|---|
1119635.9 | Nov 2011 | GB | national |
1208415.8 | May 2012 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/GB2012/052827 | 11/14/2012 | WO | 00 |