Abstract:
The invention relates to measuring visco-elastic fluid parameters, in particular, in the oil production industry, for defining heavy oil parameters during field development. 
     The method involves the excitation of a hollow resonance device oscillations by sending a continuous variable-frequency signal to two transmitting transducers located on the outer surface the said resonance device. Oscillations are recorded by a receiving transducer. Amplitude-vs-frequency response curve is plotted and effective resonance frequency ω r  is determined. Thereafter, a cavity of the said resonance device is filled with a medium under examination, and oscillations are excited to obtain associated amplitude-vs-frequency response values. An axially symmetric capacity placed in a thermostabilized chamber is used as the resonance device. While implementing the method, the attenuation factor α is also determined for the empty and test medium-filled resonance device as well as the resonance frequency shift Δω r  and the attenuation variance Δα re determined in relation to the hollow resonance device.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates to methods and devices for measuring visco-elastic fluid parameters, specifically, for measuring viscosity and elasticity of a fluid sample, using a resonance rod method while studying acoustic properties of fluids, and may be used, in particular, in the oil production industry, for defining heavy oil parameters during field development. In this regard, studying visco-elastic properties of heavy oils dependent on their temperature, in combination with oil field thermal recovery technologies, is of special interest. 
         [0002]    Therefore, the claimed method and device may be used for all fluid-saturated (derived from the English word fluid—flowing medium) non-consolidated (lack of integral unifying structure; the meaning is opposite to the “consolidated” term (lat.)—consolidatio, derived from the word con (cum)—together, conjointly and solido—compact, strengthen, merge, reinforcement/strengthening of something) media, i.e., essentially, for any visco-elastic medium, which may be placed in a tube, provided that at least satisfactory contact between the medium and the said tube is ensured, in particular, for suspensions (dispersed solutions, jellies), etc. 
       BACKGROUND OF THE INVENTION 
       [0003]    An acoustic logging method, i.e. a method for studying formation properties based on borehole measurements of elastic wave properties in a frequency range of 5-20 kHz, is a commonly used technique for well survey and development. However, visco-elastic properties of heavy oils in a frequency range that is typical for acoustic logging cannot be determined by using standard rheometric devices (laboratory instrumentation), whose operating range is within 0.001-100 Hz. This is why it is necessary to develop adequate devices and methods to perform adequate measurements of visco-elastic properties in the acoustic logging frequency range, in particular, in a broad range of temperatures. 
         [0004]    A resonance method for measuring density, viscosity and visco-elastic properties of fluid is described in U.S. Pat. No. 7,552,619. The above-mentioned method uses a resonance device that is designed as two electromechanically connected resonators and is characterized by a frequency curve with two closely located poles, which correspond to symmetrical and antisymmetrical oscillation modes. Two closely positioned peaks on the resonance curve refer to these poles. The method includes the excitation of a hollow resonance device oscillations by sending a continuous variable-frequency signal to, at least, two transmitting transducers located at the outer surface of the said resonance device, with a recording of oscillations by at least one receiving transducer, located on the outer surface of the resonance device, with a follow-up plotting of an amplitude-vs-frequency response curve, and determination of the effective resonance frequency ω r , the filling of a cavity of the resonance device with a medium under examination; the excitation of medium-filled resonance device oscillations by sending a continuous variable-frequency signal to the transmitting transducers with a follow-up pick-up of amplitude-vs-frequency response curve and determination of viscosity and visco-elastic properties based on the modified functions, which take into account measured data. 
         [0005]    Difficulties in the resonance device fabrication, which increases capital expenditures required for the implementation of the method, is a disadvantage of the above-mentioned method. At the same time, a great number of acoustic surfaces (due to the resonance device complexity) cause uncontrolled spurious resonances, thus reducing the efficiency of the measuring system and the accuracy of the method as a whole. 
       SUMMARY OF THE INVENTION 
       [0006]    This invention aims to to develop a fast-response method for measuring visco-elastic fluid parameters—modulus of elasticity and viscosity—at reduced capital expenditures, improved measurement accuracy in a broad range of temperatures; these goals constitute the claimed technical result to be achieved. 
         [0007]    In addition, an additional technical result is the possibility to implement a method for measuring visco-elastic fluid properties at acoustic logging frequencies, in particular, 5-20 kHz. 
         [0008]    The claimed technical result is reached through the excitation of a hollow resonance device by sending a continuous variable-frequency signal to, at least, two transmitting transducers located at the outer surface of the resonance device, with a recording of oscillation values by at least one receiving transducer, located on the outer surface of the resonance device, with a follow-up determination of the amplitude-vs-frequency response values. Based on the function received, effective resonance frequency ω r  is defined. Thereafter, a cavity of the resonance device is filled with a test medium and the medium-filled resonance device is excited by sending a continuous variable-frequency signal to the transmitting transducers, and associated amplitude-vs-frequency response values are received. While implementing the method, torsional oscillations are excited. Amplitude-vs-frequency response values are determined by the first oscillation mode. An axially symmetric resonance device is placed in a thermostabilized chamber. Dimensions and materials of the resonance device are selected so as to minimize the impact of oscillation modes, which differ from the effective resonance. While implementing the method, the attenuation ratio a of the empty and test medium-filled resonance device, resonance frequency shift Δω r  and the attenuation variance Δα are determined in relation to the hollow resonance device. The below mentioned formula is applied to calculate the real part of the shear modulus μ and dynamic viscosity η: 
         [0000]    
       
         
           
             
               
                 
                   
                     η 
                     f 
                   
                    
                   
                     ( 
                     
                       ω 
                       r 
                     
                     ) 
                   
                 
                  
                 
                   ρ 
                   f 
                 
               
               = 
               
                 
                   
                     
                       Δω 
                       r 
                     
                      
                     Δα 
                   
                   
                     2 
                      
                     
                       ω 
                       r 
                     
                   
                 
                  
                 
                   
                     ( 
                     
                       
                         ρ 
                         s 
                       
                        
                       
                         
                           R 
                           1 
                         
                          
                         
                           ( 
                           
                             
                               
                                 ( 
                                 
                                   
                                     R 
                                     2 
                                   
                                   / 
                                   
                                     R 
                                     1 
                                   
                                 
                                 ) 
                               
                               4 
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                     
                     ) 
                   
                   2 
                 
               
             
             , 
             
               
 
             
              
             
               
                 
                   
                     μ 
                     f 
                   
                    
                   
                     ( 
                     
                       ω 
                       r 
                     
                     ) 
                   
                 
                  
                 
                   ρ 
                   f 
                 
               
               = 
               
                 
                   
                     
                       Δα 
                       2 
                     
                     - 
                     
                       Δω 
                       r 
                       2 
                     
                   
                   4 
                 
                  
                 
                   
                     ( 
                     
                       
                         ρ 
                         s 
                       
                        
                       
                         
                           R 
                           1 
                         
                          
                         
                           ( 
                           
                             
                               
                                 ( 
                                 
                                   
                                     R 
                                     2 
                                   
                                   / 
                                   
                                     R 
                                     1 
                                   
                                 
                                 ) 
                               
                               4 
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                     
                     ) 
                   
                   2 
                 
               
             
             , 
           
         
       
     
         [0000]    where ω r —resonance frequency value for an empty capacity (no fluid), Δω r —resonance frequency shift, Δα—attenuation variance, ρ—density, R 1  &amp; R 2 —inner and outer radiuses of the capacity, respectively; lower indices f and s refer to fluid and capacity, respectively. 
         [0009]    The above-mentioned technical result is also reached through a diametrically opposed positioning of transmitting transducers on the resonance device&#39;s side wall at the bottom level. 
         [0010]    An additional technical result, associated with the possibility of implementation of the method at a required effective resonance frequency, e.g., at a value taken from the acoustic logging frequency range (5-20 kHz), is reached through defining absolute dimensions of the resonance device—length, inner and outer radius—at the established relationship between them and a selected material of the tube (characterized, first of all, by Young&#39;s modulus and shear modulus). 
         [0011]    The claimed technical result is also reached through the implementation of a device for measuring visco-elastic fluid parameters, which comprises a resonance device that comprises at least two transmitting transducers and at least one receiving transducer, both located on the resonance device&#39;s outer surface. The device also includes a computer and a sine wave generator &amp; an analog-to-digital converter, all connected to the computer, wherein an input of the analog-to-digital converter is connected to a receiving converter, whilst the output of the said sine wave generator is connected to the transmitting transducers. The device is provided with a thermostabilized chamber, in which the resonance device is positioned vertically, being hanged along the line of the node of the first mode of torsional oscillations. The resonance device is axially symmetric, whose dimensions and material are selected so as to minimize the impact of oscillation modes, which differ from the effective resonance. A tube may be used as the axially symmetric resonance device. 
         [0012]    Transmitting transducers are positioned diametrically opposed to each other on the resonance device&#39;s side wall at the bottom level. The receiving transducer is positioned on the resonance device&#39;s side wall, at the bottom level, at a distance equally remote from the transmitting transducers. The receiving transducer is connected to the input of the analog-to-digital (A/D) converter via a signal amplifier. 
         [0013]    The thermostabilized chamber is of an air type. The thermostabilized chamber is made with a cylinder-shaped side wall. The side wall is hollow and is filled with a cooling fluid. The thermostabilized chamber is provided with a microfan. The thermostabilized chamber is provided with a temperature transducer connected to a thermostat. 
         [0014]    Absolute dimensions of the resonance device are defined at the established relationship between its length and inner &amp; outer radius for a selected material. 
     
    
     
       BRIEF DESCRIPTION OF THE FIGURES 
         [0015]    The claimed invention is explained by the following drawings: 
           [0016]    FIG.  1 —Experimental unit for measuring shear visco-elastic fluid properties in a temperature range, 
           [0017]    FIG.  2 —Rated amplitude/frequency response vs temperature curve for a 150 mm-long hollow tube. The curves are numbered in increasing order of a temperature: 1° C. ( 1 ), 5° C. ( 2 ), 10° C. ( 3 ), 15° C. ( 4 ), 20° C. ( 5 ), 30° C. ( 6 ), 40° C. ( 7 ) and 60° C. ( 8 )), 
           [0018]    FIG.  3 —Variation dependent on amplitude-and-frequency response temperature for a 150 mm-long tube filled with glycerin. The curves are numbered in increasing order of a temperature: 1° C. ( 1 ), 5° C. ( 2 ), 10° C. ( 3 ), 15° C. ( 4 ), 20° C. ( 5 ), 30° C. ( 6 ), 40° C. ( 7 ) and 60° C. ( 8 ), 
           [0019]    FIG.  4 —Frequency offset and attenuation vs temperature curve for glycerin, 
           [0020]    FIG.  5 —Glycerin shear viscosity vs reverse absolute temperature, 
           [0021]    FIG.  6 —Variation dependent on amplitude-and-frequency response temperature for a 75 mm-long tube filled with oil. The curves are numbered in increasing order of a temperature: −9.3° C. ( 1 ), −4.7° C. ( 2 ), 1° C. ( 3 ), 5° C. ( 4 ), 10° C. ( 5 ), 15° C. ( 6 ), 20° C. ( 7 ), 30° C. ( 8 ), 40° C. ( 9 ) and 60° C. ( 10 ), 
           [0022]    FIG.  7 —Frequency &amp; attenuation ratio variation vs temperature for oil, 
           [0023]    FIG.  8 —Temperature-variant function of the real part and imaginary part of the shear modulus for the Mordovo-Karmalskoye field oil, 
           [0024]    FIG.  9 —Temperature-variant function of the real part and imaginary part of the shear modulus for the Mordovo-Karmalskoye field oil. 
       
    
    
     DETAILED DESCRIPTION 
       [0025]    Methods for measuring acoustic wave velocity and attenuation in high viscosity materials may be split in reflectometry methods, impulse methods, waveguide-impedance methods as well as resonance methods, which are based on measuring a resonance frequency shift and the oscillation factor-of-merit variance for the resonator during its contact with fluid under examination. Hereinafter, within the frames of the claimed engineering solution, any axially symmetric vessel (e.g., a tube) will be assumed as a “resonance device”. 
         [0026]    The developed method is a modification of the resonance rod method that calls for determination of a visco-elastic modulus through measuring a resonance frequency and a factor-of-merit (reciprocal of the attenuation coefficient) of a certain mode of natural oscillations of a sample-filled tube at increased measurement accuracy. 
         [0027]    In regard to the sample factor-of-merit, the following should be mentioned. The factor-of-merit is a parameter of an oscillating system, which defines the resonance sharpness and indicates how many times the energy reserves in the oscillator&#39;s reactive elements are greater than the energy loss at active elements for one oscillation cycle. The factor-of-merit is inversely proportional to a natural oscillation attenuation rate of the system. i.e., the higher the factor-of-merit of a resonance system, the lower the energy losses in each period. Oscillations in a system with a high factor-of-merit attenuate slowly. In other words, the factor-of-merit is defined, first of all, by a specimen material and, in case of a high factor-of-merit, is characterized, first of all, by “narrow peaks” of the effective resonance as well as by the fact that the effective resonance is far away from resonances of other oscillation modes (2 nd , 3 rd  and other modes) and other spurious resonances on the resonance curve. Certainly, this situation is typical, to a great extent, for solid homogenous materials, whose factor-of-merit is typically high. In this regard, for example, in accordance with a resonance ultrasonic spectroscopy (RUS) method which is applicable to solid materials, it&#39;s possible to unambiguously determine elasticity modulus G based on ultrasonic emission. 
         [0028]    While studying visco-elastic fluids, the above-mentioned method cannot be implemented in full, since these media are characterized by a low factor-of-merit—&lt;&lt;overstretched peaks&gt;&gt; in the resonance curve, interference of spurious resonance and other modes of oscillations onto effective resonances. In this case, it is necessary to develop methods which take into account peculiarities of a medium under study. 
         [0029]    It&#39;s common knowledge that a complex shear modulus of fluid is expressed as follows: G f =μ f  −iωη f , where μ—real part of the shear modulus, η—viscosity, ρ—density, ω—frequency, i—imaginary unit equal to √−1, index f refers to a fluid sample under study. Generally, values μ &amp; η represent the “effective” shear modulus and “effective viscosity” and are dependent on the frequency ω. 
         [0030]    As a target case, a problem of a torsional wave propagation in an infinite tube filled up with a visco-elastic material, at a radially symmetric disturbance field, was studied, and the correlation between the above-mentioned problem and the problem of defining resonance properties of torsional oscillations of the finite tube section. On the one side, the tube is a container for a sample, on the other side, the tube is a resonance system whose natural oscillations are defined by its geometry and dimensions. By selecting an optimum ratio between the tube length and its diameter &amp; thickness, at a selected tube material, we attain that a torsional resonance will be as far as possible from spurious resonances (bending resonances, etc.), i.e. we will decrease measurement accuracy and, simultaneously, enhance the factor-of-merit and, consequently, the accuracy of the claimed method for measuring visco-elastic fluid properties. 
         [0031]    On implementation of the method is where torsional oscillations are excited in the tube section by a source of harmonic current. Amplitude-vs-frequency response values are defined for the first mode of oscillations. This is selected due to the fact that torsional waves of the first oscillation mode propagate without dispersion, whilst waves of succeeding modes propagate with dispersion, i.e. the wave propagation velocity depends on frequency. Thereafter, amplitude of oscillations is recorded, dependent on the excitation frequency, i.e. the whole resonance curve is tracked (near the resonance frequency, the curve may be non-symmetrical). The width of the resonance curve allows the attenuation to be determined. Based on the data received, resonance frequency ω r  and attenuation coefficient α are determined. After that, the tube is filled with visco-elastic fluid under study, and the resonance frequency shift Δω r  and attenuation variance Δα is measured. 
         [0032]    The tube section is a major part of a sophisticated resonator, which, in addition to the tube, includes a bottom with transmitting and receiving systems attached thereto. For the tube section, the lowest resonance frequency is ω=πc s /L, where c s —velocity of transverse waves in the tube (c s  is virtually independent from the frequency), L—tube length. The availability of the bottom with the transmitting and receiving systems may be taken into consideration by modeling it like a discrete mass, which give a small correction to the resonance frequency value. 
         [0033]    Within the frames of the suggested theoretical model, a correlation between the real and imaginary parts of a shear modulus of elasticity of the sample G f =μ f −iωη f  vs resonance frequency shift Δω r  and variation of torsional oscillation&#39;s attenuation factor Δα was established: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         η 
                         f 
                       
                        
                       
                         ρ 
                         f 
                       
                     
                     = 
                     
                       
                         - 
                         
                           
                             
                               Δω 
                               r 
                             
                              
                             Δα 
                           
                           
                             2 
                              
                             
                               ω 
                               r 
                             
                           
                         
                       
                        
                       
                         
                           ( 
                           
                             
                               ρ 
                               s 
                             
                              
                             
                               
                                 R 
                                 1 
                               
                                
                               
                                 ( 
                                 
                                   
                                     
                                       ( 
                                       
                                         
                                           R 
                                           2 
                                         
                                         / 
                                         
                                           R 
                                           1 
                                         
                                       
                                       ) 
                                     
                                     4 
                                   
                                   - 
                                   1 
                                 
                                 ) 
                               
                             
                           
                           ) 
                         
                         2 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         μ 
                         f 
                       
                        
                       
                         ρ 
                         f 
                       
                     
                     = 
                     
                       
                         
                           
                             Δα 
                             2 
                           
                           - 
                           
                             Δω 
                             r 
                             2 
                           
                         
                         4 
                       
                        
                       
                         
                           ( 
                           
                             
                               ρ 
                               s 
                             
                              
                             
                               
                                 R 
                                 1 
                               
                                
                               
                                 ( 
                                 
                                   
                                     
                                       ( 
                                       
                                         
                                           R 
                                           2 
                                         
                                         / 
                                         
                                           R 
                                           1 
                                         
                                       
                                       ) 
                                     
                                     4 
                                   
                                   - 
                                   1 
                                 
                                 ) 
                               
                             
                           
                           ) 
                         
                         2 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where R 1,2 —inner and outer diameters of the tube, indices f &amp; s are referred to fluid sample and tube material. 
         [0034]    The validity of the received correlations may be confirmed by studying the Newton fluid, which is known to not have a shear elasticity (μ f =0) i.e., for defining the viscosity, it&#39;s sufficient enough to meter one of the values—frequency shift or attenuation variance. Then, Δα=−Δω, and the function looks like that: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             η 
                             f 
                           
                            
                           
                             ρ 
                             f 
                           
                         
                         = 
                           
                          
                         
                           
                             
                               
                                 ( 
                                 
                                   Δω 
                                   r 
                                 
                                 ) 
                               
                               2 
                             
                             
                               2 
                                
                               
                                 ω 
                                 0 
                               
                             
                           
                            
                           
                             
                               ( 
                               
                                 
                                   ρ 
                                   s 
                                 
                                  
                                 
                                   
                                     R 
                                     1 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         
                                           ( 
                                           
                                             
                                               R 
                                               2 
                                             
                                             / 
                                             
                                               R 
                                               1 
                                             
                                           
                                           ) 
                                         
                                         4 
                                       
                                       - 
                                       1 
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                             2 
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           
                             
                               
                                 ( 
                                 Δα 
                                 ) 
                               
                               2 
                             
                             
                               2 
                                
                               
                                 ω 
                                 0 
                               
                             
                           
                            
                           
                             
                               
                                 ( 
                                 
                                   
                                     ρ 
                                     s 
                                   
                                    
                                   
                                     
                                       R 
                                       1 
                                     
                                      
                                     
                                       ( 
                                       
                                         
                                           
                                             ( 
                                             
                                               
                                                 R 
                                                 2 
                                               
                                               / 
                                               
                                                 R 
                                                 1 
                                               
                                             
                                             ) 
                                           
                                           4 
                                         
                                         - 
                                         1 
                                       
                                       ) 
                                     
                                   
                                 
                                 ) 
                               
                               2 
                             
                             . 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0035]    An experimental unit (see.  FIG. 1 ) was used to implement the suggested method. 
         [0036]    The sample under examination ( 1 ) was placed into an aluminum tube with a bottom. The tube was hanged vertically in the air chamber ( 3 ) with thermostabilization capability to maintain constant temperature of measurements and, therefore, the accuracy of the experiment, since the viscosity of the material is a variable that to a great extent depends on a temperature. The chamber ( 3 ) is designed as a cylinder with outer thermal insulation, with a hollow side wall, in which cooling fluid injected by the thermostat ( 5 ) circulates. The cooling fluid temperature in the thermostat was automatically controlled by the outer transducer ( 4 ) placed in the chamber ( 3 ). To provide a more uniform thermostabilization bahavior, a microfan was installed in the chamber ( 3 ). Within the range of temperatures used in the measuring process (−10÷60° C.), the accuracy of the temperature value maintained in the chamber ( 3 ) was 0.1° C. 
         [0037]    While measuring, it was necessary to take into consideration that temperature variation results in a shift of the resonance frequency and the hollow tube, since its length and shear wave velocity in the tube material changes. Test results for a hollow tube are shown in  FIG. 3 . 
         [0038]    For contact-free acoustic measurements, i.e., to minimize the suspension impact on resonance properties, the tube was banded with a thin fiber in the middle, i.e. along the line of the node of the first mode of torsional oscillations (this line, in fact, is an outer boundary of the tube&#39;s cross-section plane, equally-spaced from the tube ends), and suspension fibers were attached to this fiber. 
         [0039]    Torsional oscillations were excited in the vicinity of the first natural torsional resonance. Oscillations were excited by two point power sources, i.e., the transmitting transducers (piezoelectric accelerometers), which are positioned diametrically opposed to each other on the resonance device&#39;s side wall, at the bottom level (in other words, if a tube is used, on the circumference of the tube bottom). The accelerometer&#39;s axis of sensibility was parallel to a tangent line to the resonance device surface at the accelerometer&#39;s point of contact, whilst the accelerometers themselves were oriented in the opposite direction. 
         [0040]    The third accelerometer (a receiving transducer, located similar to the transmitting transducers, was positioned on the resonance device&#39;s side wall, at the bottom level, at a distance equally remote from the transmitting transducers) was used to record a signal (in other words, in situations where a tube is used, on the tube bottom circumference, in the middle between them). 
         [0041]    For receiving a resonance curve of torsional oscillations, a programmable, frequency-adjustable, sine signal generator ( 6 ), signal preamplifier and analog-to-digital converter (A/D converter) ( 7 ) were used. The generator and A/D converter were controlled from the computing device (computer) ( 8 ), using the Matlab software. Each point of the measuring system&#39;s frequency-response curve was calculated on-line as a relation of spectral components of a given frequency for a signal registered by the receiver as well as for a signal sent to the transmitting transducer. 
         [0042]    The accuracy of the suggested measuring system is defined by a distance of oscillation modes of other types (bending, bell-shaped, second, third modes of torsional oscillations) away from the used first mode of torsional oscillations (which is, in fact, the effective resonance of the system), which, within the frames of the claimed invention, are assumed as the “spurious” resonance. 
         [0043]    For defining a spectrum of natural frequencies and for minimizing the impact of spurious resonances, a finite element method was used for the tube oscillation simulation and Comsol software was applied to select geometry. Comsol software is a software tool applied to perform thermal, acoustic, electrical, chemical and other calculations. Visit website http://www.comsol.com to receive detailed information on the software. 
         [0044]    Therefore, by selecting geometry of the tube, namely, the relationship between the length and inner &amp; outer radius for the selected tube material, we will remove the effective resonance away from spurious resonances as far as possible. At the same time, while defining absolute dimensions of the resonance device and fulfilling the ratio between geometric parameters (length vs inner &amp; outer diameters for the selected tube material), it will be possible to implement the method at a required effective resonance frequency (e.g., a value from acoustic logging frequency range (5-20 kHz)). 
         [0045]    The below-mentioned simulation results are presented as an example that illustrates the possibility of removing spurious resonances from the effective one dependent on the tube geometry, in particular, for a case when the tube length varies while other selected dimensions remain constant. 
         [0046]    The inner and outer radiuses of the tubes were 12.9 mm and 16.7 mm, respectively, the bottom thickness was 1 mm, and the only difference was the length (149.6 mm and 74.8 mm); they were selected in a way to allow measurement of visco-elastic properties at frequencies of 10 and 20 kHz. Density of the aluminum alloy, which was used for tube fabrication, was 2.78 g/cm 3 . 
         [0047]    A rubber-like viscous and elastic material with a shear module within 10 5 ÷10 6  Pa and Poisson ratio of 0.48 was used as Specimen 1. 
         [0048]    For the long tube (149.6 mm), the following natural frequencies were received through simulation: 6.5÷6.6 kHz (spurious resonance), 10 kHz (effective resonance) and 11.3÷14 kHz (spurious resonance). For the short tube (74.8 mm), similar frequencies were virtually independent from the specimen shear modulus, standing at 18.9 kHz (spurious resonance), 19.8 kHz (effective resonance) and 23.3 kHz (spurious resonance). Therefore, the used torsional resonance was removed away from the nearest resonances of other modes by 1.3 kHz for the long tube and by 0.9 kHz for the short tube, i.e. the long tube in the example under consideration was the most acceptable one. 
         [0049]    The above-mentioned example illustrates the spurious-vs-effective resonance position for a case when the tube length is the only variable, at the selected inner &amp; outer diameters for the given tube material. It&#39;s easy to understand that the above-mentioned software simultaneously processes a number of required parameters (parameter selection) with the aim of finding an optimum relationship between them. 
         [0050]    To verify the compliance between the theoretical description and experimental data received, test measurements were conducted for glycerin, i.e., for a fluid that does not have any shear modulus of elasticity within the studied frequency &amp; temperature range. An amplitude-vs-frequency response curve for a tube with a resonance frequency of 10 kHz is shown in  FIG. 3 . Frequency shift and oscillation attenuation variance for a glycerin-filled tube is shown in  FIG. 4 . 
         [0051]    In situations where there is no shear modulus (for Newton fluid), the frequency shift and attenuation variance are equal to each other and, therefore, the curves shown in  FIG. 5  should coincide. However, the difference between the curves received is likely to be explained by an incompletely adequate approximation of the contribution of other modes by a frequency-independent summand and is considered as a systematic measurement error. As follows from  FIG. 5 , the above-mentioned error is higher for a measuring system in which a shorter tube is used. For a 10 kHz measuring system, the error is 6% and reaches 15% for a 20 kHz measuring system. 
         [0052]    The glycerin viscosity calculation as per Formula (3) brought the results, which are close to values obtained during viscosity measurement, using a standard rotary rheometer. The results are presented in  FIG. 5 , where viscosity is plotted in coordinates, which show that Arrhenius law is fulfilled for glycerin, i.e. the viscosity logarithm vs reverse temperature curve is a straight line. The viscous flow activation energy for glycerin, calculated based on experimental data in accordance with Arrhenius law, was 68 kJ/mol, which is close to the spreadsheet value of 63±4 kJ/mol. 
         [0053]    Thus, experimental data received for glycerin (a Newton fluid, i.e. fluid with a zero real part of shear modulus of elasticity)—attenuation curve and resonance frequency curve, which coincide with each other (with adjustment for error) at all temperatures—confirm the validity of using the established relationship (1), (2) for defining parameters of visco-elastic fluids. 
         [0054]    The advantage of the suggested measuring method is that frequency and attenuation shift are independent from the tube length, i.e., from a resonance frequency of torsional oscillations, measurements do not require any data on the tube material shear modulus. Measuring the sample density and tube material is the only thing additionally needed. 
         [0055]    Within the frames of the claimed invention, the known data on the specimen density vs temperature relation was assumed as a prerequisite. 
         [0056]    The density vs temperature curve of the test fluid was defined using Archimedes principle. A weight of silica glass of a calibrated volume (10.000±0.005 cm 3 ) was first weighed in the air and then in the fluid under examination, whose temperature was maintained constant with an accuracy of 0.1° C. in a range of 0÷40° C. A thermocouple, which ensures the specified measurement accuracy, was used to control the test fluid temperature. The scale accuracy was 0.0001 g. A thermostat-connected coil, which embraces a cup filled with the test fluid and which provides a heat insulation coating (foam) from outside, was used for thermostabilization purpose. Within the measurement accuracy limits, the temperature curves for both used fluids (glycerin and oil) may be assumed as a linear function. The density vs temperature curves for glycerin ρ gl  (T) and the Mordovo-Karmalskoye field oil ρ oil (T), which was used as a test fluid, were as follows: 
         [0000]      ρ gl ( T )=127.70±0.3−(0.71±0.02)· T  
 
         [0000]      ρ oil ( T )=95.93±0.2−(0.74±0.007)· T′ 
 
         [0000]    where density is expressed in kg/m 3 , temperature is expressed in ° C. 
         [0057]    Here is an example of the implementation of the claimed resonance method for measuring visco-elastic fluid parameters for a case if the Mordovo-Karmalskoye field oil is used as a test fluid. The results of the experiments are presented in  FIG. 6 ,  7 ,  8 ,  9 . 
         [0058]    At the established constant temperature of 20° C., a continuous time-variant signal at a frequency of 19,800-20,600 Hz was generated in the thermostabilized air chamber ( 3 ) by the sine signal generator ( 6 ) and was sent to the transmitting transducers (point power sources) located on a hollow tube, which, in its turn, excited torsional oscillations. The relationship of geometric parameters of the tube, i.e., length and inner &amp; outer radiuses, as well as the selected tube material were as follows: tube length—75 mm, inner radius—12.9 mm, outer radius—16.7 mm, material—aluminum alloy with a density of 2.78 g/cm 3 . 
         [0059]    The third accelerometer positioned on the bottom circumference in the middle between the transmitting transducers was used to record torsional oscillations. As a result, an amplitude-vs-frequency response curve was received, based on which the resonance frequency value ω r =20,269.5 Hz and the attenuation ratio value α=8.1 s −1  were established; to define the values, a resonance curve simplex approximation method was used. 
         [0060]    Thereafter, the tube was filled up with the Mordovo-Karmalskoye field oil. The oil density calculated as per formula: 
         [0000]      ρ oil ( T )=95.93±0.2−(0.74±0.007)· T,  
 
         [0000]    was 974 kg/m 3 . 
         [0061]    While changing the frequency from 19,800 Hz to 20,600 Hz, which is sent from the sine signal generator ( 6 ) to the transmitting transducers, an amplitude-vs-frequency response curve was also obtained for the oil-filled tube. Based on the relationship obtained, the resonance frequency ω r =20,232.3 Hz and the attenuation factor α=34.1 c −1  were defined. 
         [0062]    Therefore, the resonance frequency shift is Δω r =−37.2 Hz and the attenuation variance is Δα=26 c −1 . 
         [0063]    Formulas (1) &amp; (2) were applied to calculate the real part of the shear modulus (storage modulus) μ and effective viscosity η, which make it possible to determine target “effective” visco-elastic properties of the specimen at the first torsional resonance frequency ω r . These values at a temperature of 20° C. were as follows: μ=14,630 Pa and η=0.364 Pa*s.            
         [0064]    Similar measurements were conduced at other temperatures from a range of −10 . . . +60° C. A frequency-response curve was picked up three times for each temperature, and the repeatability error did not exceed 0.3%. 
         [0065]    A pattern of the curves received ( FIG. 7 ,  FIG. 8 ) evidences that the studied oil showed its visco-elastic properties at temperatures below 20° C. The shear modulus&#39; real part is not presented for temperatures above 30° C., since the modulus values are comparable with the measurement accuracy. 
         [0066]    The viscosity-vs-temperature curve illustrated in  FIG. 9  shows that the studied oil properties deviated from the Arrhenius law that is typical for vitrifying fluids. The viscous flow activation energy is 38 kJ/mol for a high temperature (60° C.) and 78 kJ/mol for a low temperature (−9.3° C.). 
         [0067]    The suggested measuring method may be used in industrial units for on-line study of properties of visco-elastic fluid samples, passing through tubes under increased pressure.