Wave intensity analysis

Wave intensity analysis[1] provides a method to calculate the properties of arterial waves that give rise to arterial blood pressure, based on measurements of pressure, P, and velocity, U, waveforms (Figure 1). Wave intensity analysis is applicable to the evaluation of circulatory physiology and quantifying the pathophysiology of disorders such as coronary artery disease. The method is based on discrete, successive wave fronts (wavelets) and is carried out in the time domain. These wavelets travel forward and backwards in the arteries with amplitudes and . The wave intensity, , of a particular wavelet is defined asIt is related to sound intensity in acoustics and describes the power per unit area carried by the wavelet. From the theory discussed below, there is a relationship between the pressure amplitude and the velocity amplitude of a waveletwhere ρ is the density of blood and c is the wave speed of the wavelet. From these equations, generally known as the water hammer equations, it follows that the wave intensity for forward wavelets and for backward wavelets . The ability to determine the direction of a wavelet from its sign is the basis of the practical utility of wave intensity analysis.

P & U
Figure 1. The pressure and velocity measured in the ascending aorta in an elderly man.
Figure 2. The wave intensity calculated from the measured pressure and velocity measured in the ascending aorta of an elderly man.

Net wave intensity

edit

The pressure amplitude of a wavelet can be positive (compression) or negative (decompression) and the velocity amplitude can be positive (acceleration) or negative (deceleration). The measured changes   and   are the sums of the amplitudes of the forward and backward wavelets arriving at the measurement site at the time of the measurement and so the wave intensity   is sometimes called net wave intensity. The   in the Figure 2 shows the normal pattern in the aorta and illustrates four important features:

  1. A large positive peak at the start of systole indicating a dominant forward wave due to the compression of the left ventricle.
  2. A period of relatively small negative wave intensity during mid-systole indicating a small level of reflected wave activity.
  3. A smaller positive peak at the end of systole indicating that the deceleration of blood at the end of systole is predominantly due to a forward deceleration wave instead of backward reflected waves.
  4. The very low level of net wave intensity during diastole

Departures from this pattern of wave intensity is usually indicative of pathology.

 
Figure 3. The pressure (minus the diastolic pressure) measured in the ascending aorta of an elderly man separated into the forward and backward pressures using wave intensity analysis.

Separation of forward and backward waves

edit

The additivity of the forward and backward wavelets coinciding at the site of measurement at a particular time can be combined algebraically with the water-hammer equations to calculate the magnitudes of the two wavelets[2] This method assumes that the wave speed is constant. In general, the wave speed is a function of the pressure. A more complex method of separation involving integrals along the characteristics is available.[3] The forward and backward waveforms follow from summing the magnitudes of the sequential wavelets The pressure shown in Figure 1 is separated into its forward and backward components in Figure 3. This separation is carried out in the time domain and can be applied to irregular, non-periodic data. For periodic heart beats this separation coincides closely with the separation obtained using Fourier analysis methods.[4]

Theoretical basis of wave intensity analysis

edit

Wave intensity analysis is based on the 1-D equations for the conservation of mass and momentum of the blood in the elastic arteries (first published by Euler in 1775) and a relationship between the cross-sectional area of the artery and the pressure within it (generally called a 'tube law'). The resultant equations are hyperbolic in form and can be solved using the method of characteristics. The method is based on the identification of 'characteristic' directions along which the partial differential equations reduce to ordinary differential equations. For the arteries these characteristic directions were ( ) where   is the wave speed where   is the distensibility of the artery wall. These characteristic directions describe the path of the wavelets in  . The derivation of an equation for the local wave speed is a significant advantage of the method of characteristics.

Practical estimation of the local wave speed

edit

A method for estimating local wave speed from clinically measurable data is important in wave intensity analysis. The most common approach is referred to as the 'P-U loop' method. During the initial compression wave at the start of systole when the forward wave is dominant the local wave speed is estimated from the slope of the P-U loop In cases where there is no period during the cardiac cycle when it is obvious that forward waves dominate another method, generally called the 'sum of squares' method[5] is applicable: Here the sums are taken over a complete cardiac period. This estimation has been used extensively in the wave intensity analysis of measurements in the coronary arteries.[6] It should be remembered that both of these estimates are approximations and should be used with caution.

Variants of wave intensity analysis

edit

Wave intensity analysis was developed at an era when intra-arterial pressure and velocity waveforms were measured most commonly in the clinic. Other methods of clinical measurements have emerged (e.g. ultrasound and magnetic resonance imaging) and wave intensity analysis has been recast in terms of the parameters that are measured. No systematic notation has been developed to distinguish the different variants. This is a possible source of confusion.

  1. derivatives vs. differences:[7] Wave intensity is sometime defined as   which gives results that are independent of the sampling time  . This has some advantages but it means that the dimensions of wave intensity are not easily interpreted physically.
  2. area (diameter) vs. pressure[8]: High precision pressure measurements are highly invasive and so some researchers measure the instantaneous area (diameter) of the artery instead. By inverting the tube law to  , it is possible to eliminate   from the mass conservation equation. In this formulation   ( ), which has the dimensions of velocity. This variant retains the advantage that the sign of wave intensity indicates the direction of travel of the dominant wavelets.
  3. volume flow rate vs. mean velocity: Sometimes volume flow rate,  , is most easily measured clinically and wave intensity is defined as  . Since the fractional area change is relatively small, this variant of wave intensity behaves very much like the original form. However, in this variant the sign of wave intensity is not always indicative of its direction of travel.

Applications

edit

Although wave intensity has some clinical applications, its main influence has been on the basic understanding of arterial hemodynamics. The first publication pointed out that the deceleration of blood in the aorta during late systole is mainly due to forward deceleration waves rather than backward reflections from the periphery.[9] Backward wavelets are present as indicated by the negative wave intensity during mid-systole but most of the deceleration is caused by forward decompression waves generated when the contraction of the ventricle can not keep up with the flow of blood established by the initial contraction wave.

Various studies have looked at the magnitudes of the forward and backward peaks of wave intensity in various pathologies.[7] In the 1990s, Aloka introduced clinical ultrasound scanners that measured the instantaneous wave intensity directly from simultaneous measurements of diameter (by B-mode) and blood velocity (by Doppler). Recent studies suggest that wave intensity calculated from catheter measurements in the pulmonary artery can be used to differentiate between pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension.[10]

Wave intensity analysis is frequently used in the study of coronary artery hemodynamics where impedance analysis is seldom, if ever, used.[11] It played an important role in the development of the Instantaneous Wave-free Ratio as a measure of the functional effect of stenoses in the coronary arteries. This method informs an interventional cardiologist if stenting is required by measuring the ratio of downstream and upstream pressures during a period in mid-diastole when the net wave intensity is minimal. Large scale clinical trials have shown that this index is "at least as good as" the Functional Flow Resistance Index which requires the injection of a potent vasodilator.

References

edit
  1. ^ Parker, Kim H. (2009). "An Introduction to Wave Intensity Analysis". Med Biol Eng Comput. 47 (2): 175–88. doi:10.1007/s11517-009-0439-y. PMID 19205773. S2CID 738467.
  2. ^ Parker, K. H.; Jones, C. J. H. (1990). "Forward and Backward Running Waves in the Arteries: Analysis Using the Method of Characteristics". Journal of Biomechanical Engineering. 112 (3): 322–326. doi:10.1115/1.2891191. ISSN 0148-0731. PMID 2214715.
  3. ^ Stergiopulos, N.; Tardy, Y.; Meister, J.-J. (1993). "Nonlinear separation of forward and backward running waves in elastic conduits". Journal of Biomechanics. 26 (2): 201–209. doi:10.1016/0021-9290(93)90049-K. PMID 8429061.
  4. ^ Hughes, Alun D.; Parker, Kim H. (2009). "Forward and backward waves in the arterial system: Impedance or wave intensity analysis?". Med Biol Eng Comput. 47 (2): 207–10. doi:10.1007/s11517-009-0444-1. PMID 19198913. S2CID 9184560.
  5. ^ Aguado-Sierra, J.; Parker, K.H.; Davies, J.E.; Francis, D.; Hughes, A.D.; Mayet, J (2006). "Arterial pulse wave velocity in coronary arteries". 2006 International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE. pp. 867–870. doi:10.1109/IEMBS.2006.259375. ISBN 978-1-4244-0032-4. PMID 17946867. S2CID 14653454.
  6. ^ Broyd, C. J.; Davies, J. E.; Escaned, J. E.; Hughes, A.; Parker, K. (2016-02-02). "Wave intensity analysis and its application to the coronary circulation". Global Cardiology Science and Practice. 2015 (5): 64. doi:10.5339/gcsp.2015.64. hdl:10044/1/52002. ISSN 2305-7823.
  7. ^ a b Sugawara, Motoaki; Niki, Kiyomi; Ohte, Nobuyuki; Okada, Takashi; Harada, Akimitsu (2009). "Clinical usefulness of wave intensity analysis". Med Biol Eng Comput. 47 (2): 197–206. doi:10.1007/s11517-008-0388-x. PMID 18763005.
  8. ^ Feng, J; Khir, A W (2010). "Determination of wave speed and wave separation in the arteries using diameter and velocity". J Biomech. 43 (3): 455–62. doi:10.1016/j.jbiomech.2009.09.046. PMID 19892359.
  9. ^ Parker, Kim H; Jones, Christopher JH; Dawson, J Rex; Gibson, Derek G (1988). "What stops the flow of blood from the heart?". Heart and Vessels. 4 (4): 241–5. doi:10.1007/BF02058593. PMID 3254905. S2CID 26789081.
  10. ^ Su, Junjing; Manisty, Charlotte; Parker, Kim H.; Simonsen, Ulf; Nielsen-Kudsk, Jens Erik; Mellemkjaer, Soren; Connolly, Susan; Lim, P. Boon; Whinnett, Zachary I.; Malik, Iqbal S.; Watson, Geoffrey; Davies, Justin E.; Gibbs, Simon; Hughes, Alun D.; Howard, Luke (2017). "Wave Intensity Analysis Provides Novel Insights Into Pulmonary Arterial Hypertension and Chronic Thromboembolic Pulmonary Hypertension". Journal of the American Heart Association. 6 (11). doi:10.1161/JAHA.117.006679. ISSN 2047-9980. PMC 5721764. PMID 29089339.
  11. ^ Siebes, Maria; Kolyva, Christina; Verhoeff, Bart-Jan; Piek, Jan J.; Spaan, Jos A. (2009). "Potential and limitations of wave intensity analysis in coronary arteries". Medical & Biological Engineering & Computing. 47 (2): 233–239. doi:10.1007/s11517-009-0448-x. ISSN 0140-0118. PMID 19205771.