MRI METHOD FOR MEASURING VELOCITY PROFILES IN DRILLING MUD

Abstract

An MRI-based method for determining a velocity profile for a fluid flowing through a pipe, said method comprising: selecting a slice through which said fluid is flowing; selecting a pulse sequence; separating said pulse sequence into a preparation part and a readout part; applying said preparation part to said slice; waiting a predetermined time Rt; and, applying said readout part to said slice.

Claims

1. An MRI-based method for determining a velocity profile for a fluid flowing through a pipe, said method comprising: selecting a slice of said pipe through which said fluid is flowing; selecting a pulse sequence comprising a preparation part and a readout part; applying said preparation part to said slice; waiting a predetermined time R.sub.t; and applying said readout part to said slice, said readout part comprising an imaging sequence wherein said velocity profile is determinable from analysis of said imaging sequence.

2. The method according to claim 1, additionally comprising a step of selecting said imaging sequence selected from a group consisting of: a standard Spin Echo sequence, a very short T.sub.E sequence, a UTE sequence with spiral k-space sampling, and a UTE sequence with spiral k-space sampling and segmentation.

3. The method according to claim 1, additionally comprising a step of determining said velocity for at least one volume within said slice according to $v=\frac{\mathrm{Fl}}{R_t},$ where F is a fraction of a population of spins in said slice that flowed into said slice during said time R.sub.t and l is a thickness of said slice.

4. The method according to claim 1, additionally comprising a step of determining said velocity for at least one volume within said slice according to $.Math.v=\frac{\mathrm{Fl}}{R_t}.Math..Math.\mathrm{and}.Math..Math.F=1-\frac{S}{1-\beta .Math..Math.e.Math.\text{?}},.Math.\text{?}.Math.\text{indicates text missing or illegible when filed}$ where F is a fraction of a population of spins in said slice that flowed into said slice during said time R.sub.t, S is a fraction of a population of spins in said slice that resided in said slice for an entirety of said recovery time R.sub.t and which retain a magnetization imposed upon them during said preparation part, l is a thickness of said slice, β is a magnetization fraction in said slice immediately after said preparation part, and R.sub.1 is a spin-lattice relaxation rate.

5. The method according to claim 1, additionally comprising steps of (i) providing N said predetermined times (R.sub.t).sub.i, 1≦i≦N, and (ii) for each predetermined time (R.sub.t).sub.i, 1≦i≦N, executing said steps of applying said preparation part, waiting said predetermined time (R.sub.t).sub.i, and applying said readout part.

6. The method according to claim 1, wherein said pulse sequence provides a longitudinal magnetization and a transverse magnetization, and said method comprises encoding velocity information into said longitudinal magnetization.

7. The method according to claim 1, additionally comprising a step of applying said method to at least two slices of said fluid in said pipe.

8. The method according to claim 1, additionally comprising a step of selecting said velocity profile from a group consisting of a 2D profile and a 3D profile.

9. The method according to claim 1, wherein said fluid comprises drilling mud.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The invention will now be described with reference to the drawings, wherein

[0017] FIG. 1 presents a schematic flowchart of an exemplary embodiment of an MRI pulse sequence;

[0018] FIG. 2 presents a schematic flowchart of an exemplary embodiment of a process for generating a flow velocity profile in a slice of a flowing fluid; and,

[0019] FIG. 3 presents a graph comparing results for a flow velocity profile made by the method disclosed herein with results obtained from a conventional method.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0020] In the following description, various aspects of the invention will be described. For the purposes of explanation, specific details are set forth in order to provide a thorough understanding of the invention. It will be apparent to one skilled in the art that there are other embodiments of the invention that differ in details without affecting the essential nature thereof. Therefore the invention is not limited by that which is illustrated in the figure and described in the specification, but only as indicated in the accompanying claims, with the proper scope determined only by the broadest interpretation of said claims.

[0021] The present invention discloses an MRI method for measuring velocity profiles in drilling mud by drastically reducing the echo time by encoding velocity information into the longitudinal magnetization (M.sub.z) rather than into the transverse magnetization. This method obviates the need to lower the magnetic field strength in order to make the requisite measurements.

[0022] Reference is now made to FIG. 1, which presents an exemplary embodiment of a flowchart (100) giving a schematic outline of the steps for one pulse sequence of this method. During the “preparation” part (110) of a pulse sequence, at least one slice-selective inversion pulse or at least one saturation pulse is applied to a selected slice of the fluid. A recovery time R.sub.t follows the preparation part (120), after which is the “readout” part, during which an imaging sequence (130) is applied to the selected slice.

[0023] This pulse sequence (100) (preparation part (110)—recovery time (120)—readout part (130)) is repeated at least once, preferably using different R.sub.t values, more preferably for a series of R.sub.t values. In some embodiments, the imaging sequence comprises at least two spatial encoding steps. The results from a set of repeated pulse sequences generate at least one 2D or 3D flow velocity profile of at least once slice of the fluid perpendicular to the flow.

[0024] In less-preferred embodiments of the invention, it can be necessary to wait for full relaxation of the spins between repetitions of the pulse sequence, which can be very time consuming, especially if the spin-lattice relaxation rate is small.

[0025] Reference is now made to FIG. 2, which shows and exemplary embodiment of a flowchart 200 giving a schematic outline of the steps for creating a velocity profile in a slice of the fluid. The number of repeats N (210), where N is greater than 1, and the recovery time (R.sub.t).sub.i for each repeat i (220) are selected. Then, for each repeat I, a pulse sequence is executed (230), using a procedure such as the exemplary procedure (100) outlined in FIG. 1, the pulse sequence having the recovery time (R.sub.t).sub.i. If the number of repeats executed, i, is greater than N (240), a velocity profile is generated (250) and the process terminates. If full relaxation is needed between pulse sequences (260), the system waits for full relaxation (270), then executes the next pulse sequence (230) using the next relaxation time (R.sub.t).sub.i+1. If full relaxation is not needed (260), the next pulse sequence is executed (230) without a waiting time using the next relaxation time (R.sub.t).sub.i+1.

[0026] In some preferred embodiments of the invention, the imaging sequence is a very short T.sub.E standard SE sequence. In other preferred embodiments of the invention, the pulse sequence is an ultra-short T.sub.E (UTE) sequence with spiral k-space sampling. In especially preferred embodiments of the invention, a UTE sequence with spiral k-space sampling is used, with segmentation to speed up image acquisition.

[0027] In the invention herein described, for any recovery time, the spins in the selected slice are divided into three populations, one which comprises spins that flowed into the slice during the recovery time R.sub.t, one which comprises spins that resided in the slice for the whole of the recovery time R.sub.t and which retain the magnetization imposed upon them during the preparation part, and one which comprises spins that resided in the slice for the whole of the recovery time R.sub.t but are not magnetized. The fraction of the total population of spins which is in each magnetization state is F, S and S′, respectively. It should be noted that the population S′ includes both spins that were never magnetized and spins that have relaxed and lost magnetization since the start of the recovery period.

[0028] The time dependence of F is given by equation (1);

[00003]$\begin{array}{cc}.Math.F=\{\begin{array}{cc}\frac{v}{l}.Math.R.Math.\text{?}& R.Math.\text{?}<\frac{l}{v}\\ 1& R.Math.\text{?}\ge \frac{l}{v}\end{array}.Math..Math.\text{?}.Math.\text{indicates text missing or illegible when filed}& \left(1\right)\end{array}$

where v is the velocity of the flow and l is the thickness of the slice. In practice, the fraction F will be slightly less than 1 when R.sub.t=l/v, since some of the fluid is moving slower than v. As R.sub.t increases, it approaches closer to 1.

[0029] Inverting the equation, the flow velocity can be determined from

[00004]$\begin{array}{cc}v=\frac{\mathrm{Fl}}{R_t}& \left(2\right)\end{array}$

[0030] The time dependence of S is given by equation (2) :

S=(1−F)(I−βe.sup.−R.sup.i.sup.R.sup.i)

where (1−F) is the fraction of the spins that resided in the slice for the whole of the time i, (I−β) is the magnetization fraction in the slice immediately after the preparation part (110), and R.sub.1 is the spin-lattice relaxation rate. Both β and R.sub.1 can be measured independently in the absence of a flow.

[0031] From equation (2), F can be calculated from

[00005]$\begin{array}{cc}.Math.F=1-\frac{S}{1-\beta .Math..Math.e.Math.\text{?}}.Math..Math.\text{?}.Math.\text{indicates text missing or illegible when filed}& \left(4\right)\end{array}$

Using the measured results, S, F and v can be found from the measured data using equations 2 and 4.

[0032] Since the model presented here can be applied separately and independently to each voxel in the image, the velocity v of the fluid in the volume of the voxel will be determined independently for each voxel, thereby creating the desired velocity profile.

EXAMPLE

[0033] As a demonstration that the method herein disclosed can reproduce results obtained by standard methods, velocity profiles were obtained for glycerol flowing through a 16 mm diameter pipe by the method disclosed in the present invention and by a conventional method. The results are presented graphically in FIG. 3. Open circles indicate the velocity profile obtained by the method herein disclosed, and the velocity profile obtained by a conventional method is indicated by a solid line. As can be seen in the figure, the method disclosed herein accurately reproduces the flow velocity profile obtained by the conventional method, even at the outer edge of the pipe where the velocity approaches zero.