The quantitative assessment of myocardial blood flow using real-time myocardial contrast echocardiography is based on a replenishment of bubble density after bubble destruction by high-power ultrasound exposure (burst). However, all microbubbles in the myocardial vessels are not necessarily completely destroyed, which results in unreliable data of the replenishment curve analysis.

The aim of this study was to propose a corrected equation for the replenishment curve analysis based on the hypothesis in which the initial intensity just after burst should be equivalent to the baseline intensity before contrast infusion, and to examine whether the regional difference of myocardial perfusion parameters could be minimized by the use of corrected equation of replenishment curve.

Myocardial opacification of the left ventricular short-axis view was observed using low mechanical index during infusion of Definity in open-chest dogs. Bubble destruction was set in two ways, either high (0 dB) or low (−11 dB) power burst. The videointensity (VI) of baseline before contrast infusion (f-value) was assumed as an initial intensity after complete bubble destruction. Changes of the VI after burst were fitted to both exponential functions: y = a (1 − e^−βt) + c (conventional equation) and y = a (1 − e^{−β(t − d)}) + f (corrected equation). The c-value was the measured VI just after burst. The d-value was the hypothetical time of onset of the replenishment curve if all bubbles were completely destroyed. The plateau VI was defined as the A-value, which was the sum of a- and c-values or a- and f-values, respectively. The maximal difference of β-value among myocardial regions was calculated by either equation.

The A-value was almost identical in either equation regardless of the acoustic power of burst. The β-value by the conventional equation was higher after the incomplete burst than that after complete burst (0.45 ± 0.12 vs 0.54 ± 0.16). By contrast, the β-value calculated by the corrected equation was almost identical despite complete or incomplete bursts (0.46 ± 0.13 vs 0.48 ± 0.15). The maximal difference of β-value was significantly reduced by the use of corrected equation (conventional 0.24 ± 0.14 vs corrected 0.18 ± 0.10).

Variation of β-value because of the incomplete bubble destruction can be minimized by using the corrected equation: y = a (1 − e^{−β(t − d)}) + f. Further, the corrected equation can improve the regional variation of β-value.