US20060211942A1 - Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements - Google Patents

Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements Download PDF

Info

Publication number
US20060211942A1
US20060211942A1 US11/083,259 US8325905A US2006211942A1 US 20060211942 A1 US20060211942 A1 US 20060211942A1 US 8325905 A US8325905 A US 8325905A US 2006211942 A1 US2006211942 A1 US 2006211942A1
Authority
US
United States
Prior art keywords
blood
blood pressure
blood vessel
flow rate
transmission line
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
US11/083,259
Other versions
US7621876B2 (en
Inventor
Ralph Hoctor
Lawrence Hersh
Bruce Friedman
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Information Technologies Inc
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to US11/083,259 priority Critical patent/US7621876B2/en
Assigned to GE MEDICAL SYSTEMS INFORMATION TECHNOLOGIES, INC. reassignment GE MEDICAL SYSTEMS INFORMATION TECHNOLOGIES, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: FRIEDMAN, BRUCE A., HERSH, LAWRENCE T., HOCTOR, RALPH T.
Publication of US20060211942A1 publication Critical patent/US20060211942A1/en
Priority to US11/669,602 priority patent/US8388544B2/en
Application granted granted Critical
Publication of US7621876B2 publication Critical patent/US7621876B2/en
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/022Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A61B5/02125Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/04Measuring blood pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4209Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
    • A61B8/4236Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames characterised by adhesive patches
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4472Wireless probes

Definitions

  • the invention generally relates to an apparatus and method for measuring blood pressure, and more particularly to a continuous non-invasive blood pressure measurement and monitoring method and apparatus.
  • One technique for providing a continuously measured blood pressure is to insert a saline filled catheter through the patient's vascular system to the point at which the blood pressure measurements are desired.
  • the catheter is connected to a pressure sensor, which measures the pressure in the vessel.
  • a catheter with the pressure sensor at the tip that directly senses the blood pressure can be inserted into the patient's vascular system.
  • the time of the pressure pulse's origin at the heart is determined from an EKG signal and the arrival time of the pulse is measured at another location on the patient. Based on these measurements, the instantaneous blood pressure is determined.
  • initialization data is obtained from a cuff based blood pressure determination and is used along with the time required for the pulse wave to travel between two points to calibrate formulas developed for continuously estimating blood pressure.
  • changes in the pulse transit time can be related to changes in the blood pressure.
  • the measurement of the area of the blood flow passageways, the blood flow, and PWV along with the subsequent blood pressure estimate are determined by the use of an assortment of complex adjustment factors.
  • any type of non-invasive continuous blood pressure monitoring system various factors can affect the accuracy of the measurement. For example, changes in the physiological state of the patient can bring about changes in the arterial wall elasticity. In general, changes in the arterial wall elasticity will affect the measured PWV. If the elastic modulus of the arterial wall changes, the same pressure may then need to be associated with a different cross-sectional area and PWV. If the operating point obtained using a blood pressure cuff calibration, then any changes in the arterial elasticity would require re-calibration, and if no re-calibration were performed, errors in the pressure estimation can clearly occur.
  • U.S. patent application Ser. No. 10/749,181 commonly assigned with the present application and incorporated herein by reference, teaches a method that utilizes pulse wave velocity to create a continuous, non-invasive blood pressure measurement.
  • the '181 application teaches a method of measuring the pulse wave velocity within a patient and a method of relating such pulse transit time to blood pressure.
  • the '181 application teaches a method of accurately measuring the PWV, the relationship to blood pressure can be criticized using some of the same grounds set forth above.
  • the present invention is related to an apparatus and method for continuously and non-invasively monitoring the blood pressure that utilizes a transmission-line model and transcutaneous ultrasound measurements along a blood vessel of a patient.
  • the resultant blood pressure measurements are displayed, in real time, to provide a continuous, non-invasive blood pressure measurement.
  • the apparatus of the invention includes at least a pair of ultrasound patch probes positioned along the skin of a patient.
  • the ultrasound patch probes are separated by a relatively small distance and are positioned above the same blood vessel.
  • the ultrasound patch probes are used by the system and apparatus of the present invention to make measurements relative to the blood vessel on a real time basis. Specifically, the ultrasound patch probes are used to determine the radius of the blood vessel and the blood flow rate at at least two distinct locations along the length of the blood vessel.
  • the ultrasound area and flow measurements are specific to the particular artery being used, unlike pleythysmography.
  • the system of the present invention utilizes a time-domain model of a transmission line where the electrical parameters of capacitance, inductance and resistance of the transmission line model are used in analogy to blood flow characteristics.
  • the capacitance of the electric circuit is analogous to the compliance of the artery
  • the inductance of the electric circuit is analogous to the inertia of the blood
  • the resistance is analogous to the viscosity of blood flow.
  • the input and output current sources in the electrical transmission line model is analogous to the blood flow rate, which are determined at both ends of the modeled arterial segment using ultrasound measurement techniques.
  • a transmission line simulator becomes a circuit representation of the artery, where the voltage across the capacitor of the transmission line model correlates to blood pressure.
  • the system of the present invention uses Kalman filter techniques to provide the best estimates of the state variables, which include the voltage across the capacitor in the transmission line model, which is analogous to a blood pressure measurement.
  • the system and apparatus of the invention utilizes a blood pressure cuff to provide initial calibration measurements for the blood pressure analysis performed in real time.
  • the initiation blood pressure information is used to develop an initial state vector that forms part of the Kalman equation estimates.
  • the system monitors the Kalman filter results, as well as other system variables, and automatically recalibrates the system using measured parameters when the current results of the Kalman filtering or the system variables, as compared to the initiation values, exceed a pre-defined threshold value.
  • the continuous, non-invasive blood pressure monitoring technique utilizes the pair of ultrasound patch probes to make accurate, real time measurements of physical characteristics of the blood vessel through which the blood flow is used to make a blood pressure measurement.
  • the ultrasound patch probes provide accurate data that is used to provide better estimates for the blood flow and size of the blood vessels through which the blood pressure is being determined. In this manner, the continuous, non-invasive blood pressure monitoring method provides greater accuracy as compared to prior methods.
  • FIG. 1 is a illustration of the system of the present invention used to carry out the continuous, non-invasive blood pressure measurement technique
  • FIG. 2 is a drawing showing a cross-sectional view of a typical cMUT cell
  • FIG. 3 is a drawing showing the CNIBP measurement concept in accordance with the disclosed embodiment of the invention.
  • FIG. 4 is a schematic illustration of the transmission line model utilized to develop the blood pressure measurements of the present invention.
  • FIG. 5 a is a graphic illustration of simulated proximal and distal blood flow
  • FIG. 5 b is a graphic illustration of reconstructed pressure and actual pressure in accordance with the present invention.
  • FIG. 6 is a graphic illustration of the error between estimated pressure and the known pressure with a prediction error.
  • FIG. 7 is a high level data flow diagram illustrating the flow of data used to calculate the blood pressure estimates of the present invention.
  • a blood pressure monitor 10 which may be part of a larger patient monitor 12 is shown.
  • the patient monitor 12 may provide medical personnel with information regarding different physiological functions of a patient 14 .
  • the patient monitor 12 may include an electrocardiograph (ECG), a pulse oximeter, and the blood pressure monitor 10 .
  • ECG electrocardiograph
  • pulse oximeter a pulse oximeter
  • the patient monitor 12 is connected to a central processor 16 that is used to carry out the continuous, non-invasive blood pressure monitoring functions to be described in greater detail below.
  • the central processor 16 is coupled to a blood pressure cuff monitor 18 , which in turn controls and monitors the readings from a conventional blood pressure cuff 20 .
  • the blood pressure cuff 20 is shown positioned around the arm of a patient and can be operated in the conventional manner to obtain actual systolic and diastolic blood pressure readings from the patient 14 .
  • the central processor 16 is also coupled to an ultrasound activation circuit 22 , which in turn is coupled to a first ultrasound patch probe 24 and a second ultrasound patch probe 26 . As illustrated in FIG. 1 , the first ultrasound patch probe 24 and the second ultrasound patch probe 26 are spaced from each other by a distance D along the forearm of the patient. Although the first and second ultrasound patch probes 24 , 26 are shown positioned along the forearm of the patient 14 , it should be understood that the ultrasound patch probes 24 , 26 could be positioned at various other locations along the patient.
  • each of the ultrasound patch probes 24 , 26 is positioned on the arm of the patient in such a manner as to enable non-invasive insonification of an artery.
  • Each of the ultrasound patch probes 24 is positioned generally above the radial artery 28 of the patient.
  • the blood pressure cuff is positioned to sense blood pressure in the brachial artery of the patient.
  • non-invasive ultrasound can provide volumetric flow, arterial lumen area and pulse wave velocity information.
  • ultrasound measurements can be taken in such a way that all the data describes a single, uniform arterial segment (e.g., the brachial or common carotid arteries can be used). This means that the measured PWV is directly related to the mechanical properties of the segment of elastic tube (artery) for which the blood volume is being measured.
  • the operator of the ultrasound device is eliminated through the use of software that automatically locates the artery in the ultrasound data, e.g., using known edge detection techniques. Autonomous operation of the ultrasound system allows it to report blood pressure and blood flow traces to the clinical users without those users having to interpret an ultrasound image or operate an ultrasound imaging device.
  • micromachined ultrasonic transducer patch probe technology allows ultrasound data to be taken using a thin, lightweight probe that adheres to the patient's skin.
  • semiconductor processes have been used to manufacture ultrasonic transducers of a type known as micromachined ultrasonic transducers (MUTs), which may be of the capacitive (cMUT) or piezoelectric (pMUT) variety.
  • cMUTs are tiny diaphragm-like devices with electrodes that convert the sound vibration of a received ultrasound signal into a modulated capacitance. For transmission, the capacitive charge is modulated to vibrate the diaphragm of the device and thereby transmit a sound wave.
  • MUTs can be made using semiconductor fabrication processes, such as microfabrication processes grouped under the heading “micromachining”. As explained in U.S. Pat. No. 6,359,367:
  • MUT patch allows a clinician to attach the transducer to the patient's skin.
  • the MUT patch is lightweight and flat.
  • the ultrasound patch probes will be described that are made up of capacitive micromachined ultrasonic transducers (cMUTs).
  • cMUTs capacitive micromachined ultrasonic transducers
  • the patches could instead employ pMUTs.
  • the present invention encompasses not only a device, but also methods for continuous non-invasive blood pressure monitoring, and that the basic method is not limited to devices that employ MUT technology. Instead, the basic method disclosed herein can also be performed by appropriately programmed ultrasound imaging systems having probes comprising conventional piezoceramic transducer elements.
  • the digital signal processor (DSP) 16 Based on acoustic data from the cMUT ultrasound patch probes 24 and 26 and pressure data from the blood pressure cuff 20 , the digital signal processor (DSP) 16 autonomously computes various parameters, including the patient's estimated blood pressure. In addition, the processor software forms and steers ultrasound beams over the volume of space in front of the cMUT ultrasound patch probes 24 and 26 . Using both echo intensity and Doppler blood flow information, the software maps out the location of the vessel, its cross-sectional area, and the area increase at the time of passage of the systolic pressure wave. The pulse wave velocity can also be measured. This will yield useful information about the mechanical properties of the blood vessel wall.
  • the blood pressure cuff 20 is periodically inflated by the central processor 16 to provide the peak and steady state pressures at fixed points in time. More specifically, the processor 16 activates a pump to inflate the blood pressure cuff 20 with air. A blood pressure cuff monitor 18 detects the pressure in the cuff 20 . The analog output of the cuff 20 is converted to a digital value by the cuff monitor 18 , which digital pressure value is used by the processor 16 . Initially, the processor 16 triggers cuff inflation in response to an activation command input by the system operator via an operator interface (not shown). Thereafter, the processor 16 triggers re-inflation of the cuff 20 for re-calibration whenever the estimated compliance value changes by more than some predetermined amount from its last calibration level. The results of the computations are displayed on a screen (not shown) of the patient monitor 10 .
  • the pair of ultrasound patch probes 24 , 26 and the blood pressure cuff 20 pneumatics in the disclosed embodiment are controlled and monitored by electronics coupled to the processor 16 by means of cables, obviously a wireless electrical coupling could be provided.
  • a typical cMUT transducer cell 2 is shown in cross section.
  • An array of such cMUT transducer cells is typically fabricated on a substrate 4 , such as a heavily doped silicon (hence, semiconductive) wafer.
  • a thin membrane or diaphragm 8 which may be made of silicon nitride, is suspended above the substrate 4 .
  • the membrane 8 is supported on its periphery by an insulating support 6 , which may be made of silicon oxide or silicon nitride.
  • the cavity 9 between the membrane 8 and the substrate 4 may be air- or gas-filled or wholly or partially evacuated.
  • a film or layer of conductive material such as aluminum alloy or other suitable conductive material, forms an electrode 11 on the membrane 8 , and another film or layer made of conductive material forms an electrode 13 on the substrate 4 .
  • the bottom electrode can be formed by appropriate doping of the semiconductive substrate 4 .
  • the two electrodes 11 and 13 separated by the cavity 7 , form a capacitance.
  • an impinging acoustic signal causes the membrane 8 to vibrate
  • the variation in the capacitance can be detected using associated electronics (not shown in FIG. 2 ), thereby transducing the acoustic signal into an electrical signal.
  • an AC signal applied to one of the electrodes will modulate the charge on the electrode, which in turn causes a modulation in the capacitive force between the electrodes, the latter causing the diaphragm to move and thereby transmit an acoustic signal.
  • the first data type is a sequence of B-mode images providing several short-axis views of the artery (in planes perpendicular to the direction of blood flow). Two or more short-axis views are taken along the length of the artery.
  • the center of the artery can be found at each measurement location using edge detection techniques.
  • the center of the artery can then be estimated for any location on the patch using interpolation.
  • the B-mode data is used to estimate the vessel diameter for the blood pressure calculation. More specifically, it is assumed that that the cross section of the lumen is a circle, so the estimated diameter can be used to calculate the lumen area A.
  • the second data type is a set of M-mode lines to be used for Doppler processing.
  • the M-mode lines remain locked onto the center of the artery using information from the B-mode data.
  • the multiple M-mode lines can be used with no Doppler steering angle for determining the velocity of the vessel wall or with a Doppler steering angle for determining blood velocities.
  • the M-mode data, either wall or blood velocities, will be used for estimating the pulse wave velocity needed for the blood pressure calculation.
  • the two types of data are acquired by interleaving the acquisitions of the short-axis B-mode beams and the Doppler M-mode beams.
  • a repeating scan sequence consists of a single B-mode beam at each location and all the M-mode beams.
  • the B-mode beams are translated across the patch while the M-mode beams remain in the same position. This is repeated until the desired width for the B-mode image is covered, generating a single frame of the B-mode data.
  • the steps are repeated using any updated information about the vessel center.
  • the result of the acquisition is M-mode data at a high repetition rate ( ⁇ 2 kHz), as required for Doppler processing, and the B-mode images at a slower rate ( ⁇ 100 Hz) determined by the required image size for vessel tracking.
  • the B-mode beams are translated across the cMUT patch by translating an annular array of activated sub-elements across the patch.
  • a uniform translation of the beamforming coefficients produces a new beam at a different location. Repeated frequently, this generates a rectilinear two-dimensional image.
  • An additional bilinear term in the beamforming coefficients produces a beam directed away from the normal. Such beams are transmitted in the M mode to acquire data for Doppler processing that computes the blood velocity based on the Doppler Effect.
  • M-mode beams that are transmitted generally normal to the artery are used to acquire data for Doppler processing of the arterial wall velocity during deflection of the wall as the pulse wave travels down the artery.
  • the maximum deflection occurs when the wall velocity is zero.
  • each of the first and second ultrasound patch probes 24 , 26 can be actuated by the activation circuit 22 to emit an ultrasound pulse 30 .
  • the ultrasound pulse 30 allows the ultrasound activation circuit 22 to make multiple measurements and determinations relative to a uniform stretch of an artery, such as the brachial artery 28 positioned between each of the ultrasound patch probes 24 , 26 .
  • the patch probes 24 , 26 can be spaced on a range of artery in which there would be several bifurcations.
  • the ultrasound patch probes 24 , 26 are shown spaced over a length of the unbifurcated brachial artery.
  • the central processors 16 can determine the internal diameter of the artery 28 at the location beneath the first ultrasound patch probe 24 and at the location beneath the second ultrasound patch probe 26 .
  • the central processor 16 can utilize the Doppler ultrasound flow measurement technique to determine the instantaneous flow of blood at a first location beneath the first ultrasound patch probe 24 and at a second location beneath the second ultrasound patch probe 26 .
  • the central processor 16 through the ultrasound activation circuit 22 and the pair of ultrasound patch probes 24 , 26 can determine the arterial diameter and blood flow rate at two separate locations along the radial artery which are spaced by the distance D.
  • the data collection scheme associated with the present invention requires that at least two ultrasound probes be deployed along the length of a short arterial segment. In the monitoring application, it is anticipated that these will be the ultrasound patch probes 24 and 26 . Each of these probes must have a steerable beam and must be able to operate at an angle with respect to the skin normal so that it can collect Doppler blood flow measurements. Such a probe can also be steered to the normal angle and used to measure the area of the artery's lumen. Optionally, other probes may be deployed between the two end probes, all of them at known locations along the line of the arterial segment.
  • the blood pressure monitoring system and method of the present invention utilizes an algorithm in which a length of an unbifurcated artery is represented by a transmission line model 32 .
  • the transmission line model 32 shown in FIG. 4 is a well known model and technique used for determining the phase velocity of electrical waves propagating in a cable or transmission line.
  • the transmission line model 32 can be analogized to the flow of blood in a distensible tube, and thus used to make blood pressure measurement.
  • the transmission line is an analog for such a tube, with electric charge corresponding to fluid volume, current corresponding to volumetric flow and voltage corresponding to pressure.
  • the specified flow measurements are like direct measurements of current in the analogous electric circuit.
  • the transmission line model 32 will be represented by a ladder network, each stage of which represents a short sub-segment of the arterial segment.
  • the operation of this network can be represented by a set of simultaneous, first-order differential equations, which can be used to write a discrete-time state variable model for the system, given a sample rate.
  • state variable model half of the state variables will be pressures, the other half will be flows.
  • This model can then be driven by ultrasound measurements and the pressure estimate will be given by one of the pressure state variables.
  • the model must have parameters that accurately represent the physical system in order to give meaningful results.
  • the general method of allowing a state variable model to evolve in time is reminiscent of the numerical solution of an initial value problem using Euler's method.
  • the present apparatus approach differs from a typical initial value problem in that the equation to be solved is inhomogeneous, so if the system is stable, the particular solution associated with the driving function dominates the general solution. It is for this reason that the initial values of the state variables are ignored, except as described below in relation to the operating pressure.
  • the transmission line model 32 is represented in the signal processing method as a RLC lossless ladder network shown in FIG. 4 .
  • the term “lossless” indicates that no mass (i.e. blood) is lost from the input to the output of the network, such as in a length of artery without any branching.
  • the analogous quantity to the inductance L of the ladder network is the mass per unit length (per ladder network section) and the analogous quantity to the capacitance C is the elastic compliance of the tube.
  • the resistors R represent energy loss due to the viscosity of blood.
  • the network is driven at both ends by current sources 44 and 46 , representing the ultrasonic blood flow measurements.
  • the RLC ladder network in FIG. 4 is described as being a “lossless” network, it should be understood that if the method of the present invention is used on more complex artery structures, such as the radial artery, or a longer arterial segment, some of the blood mass is diverted in branching. In such an embodiment, a transverse resistance, in parallel with each capacitance, can be included in the ladder network.
  • a transverse resistance in parallel with each capacitance, can be included in the ladder network.
  • the present method does not use ultrasound measurements on the branching arteries, a relationship could be developed between changes in the dimensions and flow of the artery segments being measured and the branches that are not being measured. This relationship can be used to estimate the transverse resistance dynamically. Therefore, the presence of a transverse resistance can be handled within the scope of the described technique.
  • the transmission line model includes multiple RLC ladder network stages that each include an inductance 36 , a capacitance 38 and a series resistance 40 .
  • the ladder network 34 includes n different stages.
  • the order of the transmission line model which is related to the number of serial sections employed, is to be determined by the length of the arterial segment to be modeled, and the parameters represent values of an underlying distributed parameter, expressed per length of the individual section. For example, if the distance between the proximal and distal flow measurements is 5 cm (which is typical), then a 10-section transmission line model would have its compliance value given in units equivalent to “compliance per 5 mm section of tube”. When the length of the individual section is short enough, then the compliance value in the model gives an approximate relationship between area and pressure. This corresponds to the definition of compliance used, for example, in the, Bramwell-Hill equation, discussed below.
  • the state variables are the current through each of the inductors 36 as well as the voltage across each of the capacitors 38 .
  • the ground 48 of the transmission line model 32 is assumed to be zero such that the voltage across capacitor 38 represents the voltage at node 50 .
  • the voltage in the transmission line model corresponds to pressure in the arterial transmission line model. (The ground value should actually be equal to the pressure in the venous system, but this is so low compared to the arterial pressures that it is negligible.)
  • Equation (3) gives a continuous-time expression for the transmission line model.
  • the signal processing operation operates in discrete time, since the proposed measurements are intrinsically discrete time.
  • a new state feedback matrix G is defined, which is the state transition matrix of the continuous system over the sample period.
  • This well-known-approach to discretization of a continuous time linear system represents a delay-free, zero-order hold assumption for the digitization of the driving signals. This is equivalent to representing the driving input signal by a piecewise constant approximation.
  • Other digitization models can also be employed, such as a first-order hold, which represents the input signal as piecewise linear. Higher order hold models may be required at lower sample rates.
  • the basic operation of the ladder network of the transmission line model is to integrate the two flow signals to get a time-varying volume estimate, and to scale that using the compliance (C) to obtain a pressure (P) estimate.
  • the series resistors 40 representing blood viscosity, are very important in implementing the low-pass frequency response characteristic that performs the integration.
  • the network implements a delay for the flow signals. Because of this delay, blood observed entering the system at a given time at one end will be canceled by blood leaving the system at a later time at the other end, with the interval determined by the implemented delay. This is true for the forward flow wave and also for the reflected, or retrograde, flow wave. When the delay does not match the actual delay in the system, the flow is cancelled either too soon or too late, with the result that the pressures are wrong.
  • the phase velocity of the network depends on the product LC.
  • a simulation was performed by taking a smoothed version of an impedance pleythysmography measurement from the femoral artery of a human subject, and taking its first difference to simulate a forward flow at the proximal measurement site.
  • Impedance pleythysmography measures the total blood volume at the measurement site.
  • the simulated proximal forward flow was delayed to get the forward flow at the distal site.
  • a reflected flow wave was then added to both.
  • the reflected wave as modeled as originating at a point 2 centimeters distal to the distal measurement site that reflected the wave with a reflection coefficient of 0.3.
  • FIG. 5 b shows the result of inputting the flow waveforms on FIG. 5 a in to the model above.
  • the pressure variables of the model were initialized to a value of 70.
  • the reconstructed pressure 78 is generally well-estimated, as compared to the actual pressure 80 .
  • the parameters of the model and the scale of the flow measurement were obtained by searching until a correct representation of the flow waveform occurred. Note that the simulation methodology did not enforce a consistent relationship between the delay, frequency response and scaling of the pressure variables. This example has been provided only to show the capability of the technique in producing a continuous non-invasive blood pressure waveform.
  • the Kalman filter algorithm is composed of three main operations, which together recursively update the system state and its covariance matrix, which in turn determines the relative weight given to the observations as opposed to that given to the current state estimates.
  • the state update mechanism is taken to occur with additive noise.
  • process noise is omitted from the model, or equivalently, the covariance matrix is set to zero.
  • an observation matrix H that specifies the observed state variables.
  • H has the same number of columns as the dimension of the state vector and the number of rows is equal to the number of additional observations. In each row, there is only one non-zero entry, which is unity if the measured state variable is a flow and C if the state variable is a pressure and the measurement is an area.
  • the components of z(k) are either noisy flow measurements or noisy area measurements.
  • the innovations sequence is simply the difference between the observation and the observation that would be predicted by the current state vector.
  • the Kalman gain matrix computed in Equation 11 is used to weight the innovations sequence to produce an updated estimate of the state vector at the k th sample time.
  • ⁇ tilde over (x) ⁇ ( k ) ⁇ circumflex over (x) ⁇ ( k ) +K ( k )( z ( k ) ⁇ H ⁇ circumflex over (x) ⁇ ( k )) (12) Note that this equation replaces the current state vector by a weighted sum of itself and the observation.
  • the third and last step in the Kalman algorithm is to project the state vector and its covariance matrix to the next sample time:
  • x ⁇ ⁇ ( k + 1 ) G ⁇ x ⁇ ⁇ ( k ) + F ⁇ [ Q p ⁇ ( k ) Q D ⁇ ( k ) ] ( 13 ⁇ a )
  • P ⁇ ( k + 1 ) G ⁇ [ ( I - K ⁇ ( k ) ⁇ H ) ⁇ P ⁇ ( k ) ] ⁇ G T ( 13 ⁇ b )
  • the effect of the Kalman filter is to replace the observed state variables at every step with a weighted average of their observed and predicted values.
  • the weights depend to a large extent on the value of the noise variance on the main diagonal of R.
  • the artery compliance is the slope of the transfer characteristic between pressure and the cross-sectional area of the vessel lumen. This characteristic is nonlinear because of the collagen fibers in the arterial wall, which are relaxed until a certain stage of the distension of the vessel, and then become taut, decreasing the elasticity. Thus, the compliance can be considered to be dependent on pressure. This nonlinear relationship can be represented by a piecewise linear relationship such that the compliance is different in a number of disjoint pressure ranges. This kind of pressure-dependent compliance can easily be incorporated into the time-domain approach above.
  • Equations (11), (12) and (13) are then executed using that set of matrices that is associated with the estimated pressure at step (k- 1 ).
  • Volumetric flow estimation proceeds in the same way as ultrasound color flow mapping over a limited sector.
  • Color flow mapping is a standard feature of most diagnostic ultrasound devices that estimates blood velocity at a number of locations on the image plane by computing the mean frequency of the spectrum of the returned ultrasound signal. This mean frequency corresponds to a mean Doppler shift induced by the moving blood.
  • This measurement process involves insonifying a volume from which a blood velocity is desired, and this occurs from three to sixteen times for every estimate. This is done along several ultrasound A-lines. Once a set of blood velocity estimates is available at a number of locations covering the vessel lumen, the estimates are all multiplied by the area to which they individually apply and the products are summed. This implements an approximate integration over the lumen and the result is a volumetric flow.
  • ultrasound flow mapping can produce blood velocity estimates that are highly biased down from the true blood velocity if the ultrasound return contains significant energy from stationary tissue reflections. For this reason, the data used to compute the mean frequency estimates (the Doppler shifts) is typically high-pass filtered. This means that very low flows, such as might be encountered during diastole in a peripheral artery, may be missed. Since the majority of the flow undoubtedly occurs when the blood velocity is high, this may not be a major problem, but still, some portion of the flow is bound to be neglected due to this effect.
  • Another effect that can change the performance of the blood pressure measurement is the dependence of the ultrasound velocity measurement on angle.
  • the measured phase change is indicative of the projection of the blood velocity vector on the ultrasound line. If the various flow measurements are based on different Doppler angles, they can be incommensurate.
  • the lossless circuit of FIG. 4 requires that positive current entering the system at the proximal source 44 be canceled by negative current entering the system at the distal source 46 after some delay. If the angles of the two ultrasound lines are different, this cancellation will not take place, and the blood pressure will not be correctly computed. Note that if the ultrasound probes 24 , 26 are mechanically connected to each other by a rigid structure, identical insonification angles are easily obtained.
  • a data-adaptive processing can be utilized to find the undetermined free parameters of the model. This method, however, requires that at least one additional measurement be made, over and above the two required flow measurements described previously.
  • a prediction error is the error between the response of the filter to a known input and the response of the real system to that input. This is a well-known principle in spectrum estimation and system identification. In accordance with the invention, this approach is employed to indicate the degree of agreement between the parameters of the transmission line model and the state of the physical system.
  • Equation (10) For this purpose we use the model of the observations defined in Equation (10) in connection with the Kalman filter approach.
  • the adaptive filter works by applying multiple sets of parameters to the data in parallel over a fixed period of time, and computing E P for each.
  • the set of parameters producing the lowest prediction error energy is taken to be the correct set of parameters.
  • the number of such parallel computations will obviously be limited by computational resources; however, even as few as three parallel computations will allow the adjustment of a single parameter. If one parameter set is chosen to be above the nominal parameter value, one chosen to be below it and another to be exactly at the nominal value, then after each observation period, the new nominal value can be that one of the three with the lowest E P , and the process can be repeated. This allows the system to track a time-varying parameter with minimal computational resources.
  • the parameters that will typically be varied in this adaptation will be the (Ea) product and the blood viscosity, since the L parameter and the other components of the C and R parameters can be measured using ultrasound. Measurement of the characteristic pulse wave velocity of the artery can also be used to help estimate C. It is also possible to limit adaptation of these parameters to a region around the nominal value of ⁇ and the value of C predicted by the Bramwell-Hill equation.
  • the simulation example shown in FIGS. 5 a and 5 b above was repeated with a third, central, flow measurement, and with a range of compliance values around the optimum.
  • the results are depicted in FIG. 6 , in which the squared estimation error between the estimated pressure and the known pressure (line 82 ) is plotted against the prediction error (line 84 ).
  • the two curves are in different units and they have been scaled and shifted so that they could be plotted on the same graph.
  • the minimum value of each is marked by a circle 86 and 88 .
  • the plot shows that the prediction error curve 84 is smooth, and that it can be made to have the same minimum as the actual estimation error 84 , provided the non-adapted parameters are correct.
  • Adaptation of the system parameters using prediction error is consistent with incorporation of the additional measurements using the Kalman filter. All that is required is the execution of the operations given by equations (11), (12) and (13b) in addition to the normal operation given above. There is a potential problem in this approach, however. To the extent that the observation used to generate prediction error updates the state variable to which it is compared, the sensitivity of the prediction error to model errors will be reduced. The degree to which the observations update the state variables is controlled by the noise variance specified in the R matrix; thus, adaptive Kalman filter operation will require that the adaptation occur during fixed periods, during which the noise variance is set to a high level. In periods when the adaptation is not taking place, the noise variance should be set back to a realistic level.
  • the system receives a brachial artery diameter measurement in step 60 .
  • the artery diameter measurement is an ultrasound measurement taken from both of the ultrasound patch probes 24 , 26 illustrated in FIG. 1 .
  • the input blood flow at ultrasound patch probe 24 is measured using ultrasound techniques.
  • the output blood flow is measured at block 64 by the ultrasound patch probe 26 .
  • the three ultrasound input measurements in blocks 60 , 62 and 64 are all transferred to a transmission line simulator shown at block 66 .
  • the transmission line simulator 66 utilizes the input measurements, in connection with additional values to be described below, to carry out the calculation steps previously described.
  • the transmission line simulator 66 includes the Kalman filtering as well as the state variable formulas to generate the blood pressure estimates, which are displayed in block 68 .
  • the data flow also includes a pulse wave velocity estimator, which uses known measurement techniques, such as those taught in U.S. patent application Ser. No. 10/749,181 to calculate the pulse wave velocity.
  • the pulse wave velocity estimator 70 also utilizes the input and output blood flow measurements from blocks 62 , 64 to make an accurate pulse wave velocity estimate.
  • the data flow also includes an elastic estimator 72 .
  • the elastic estimator 72 utilizes the pulse wave velocity measurement from step 70 and an empirical formula to calculate the elasticity of the blood vessel through which the blood pressure is being determined. As illustrated in FIG. 7 , the pulse wave velocity and elastic estimates are sent to the transmission line simulator 66 such that the transmission line simulator 66 can carry out the calculations to determine the blood pressure estimate.
  • the method and apparatus are described as presenting a continuous blood pressure measurement for a patient using ultrasound patch probes placed on the skin of the patient.
  • this continuous measurement will be represented as a waveform that can be displayed.
  • the method of the present invention is equally applicable for estimating the systolic and diastolic blood pressure for each beat of the patient's heart. The beat by beat estimate of the systolic and diastolic pressure would provide a continuous blood pressure measurement as the term is used herein.

Abstract

A method and technique for the continuous, non-invasive measurement of blood pressure. The blood pressure measurement technique of the present invention utilizes ultrasound measurements to determine the diameter of the blood vessel in which the blood pressure is being measured as well as the flow rate of blood at both an input point and an output point along the blood vessel. The system utilizes a transmission line model to relate various blood vessel measurements with electrical components. The transmission line model, in combination with data management techniques including state variable representations and Kalman filtering, is used to develop a blood pressure measurement in real time.

Description

    BACKGROUND OF THE INVENTION
  • The invention generally relates to an apparatus and method for measuring blood pressure, and more particularly to a continuous non-invasive blood pressure measurement and monitoring method and apparatus.
  • In the case of a hospitalized patient, it has long been desired to be able to provide non-invasive beat-by-beat (continuous) systolic and diastolic blood pressure values. Unfortunately, a practical and reliable solution for obtaining this type of information has yet to be developed.
  • One technique for providing a continuously measured blood pressure is to insert a saline filled catheter through the patient's vascular system to the point at which the blood pressure measurements are desired. The catheter is connected to a pressure sensor, which measures the pressure in the vessel. As an alternative method, a catheter with the pressure sensor at the tip that directly senses the blood pressure can be inserted into the patient's vascular system. Although both of these techniques have proven effective and continuously monitor a patient's blood pressure, both techniques involve making an incision into the patient's skin and inserting the catheter into the blood vessel. As a consequence, this invasive procedure entails some risk of complication to the patient and is in most cases undesirable.
  • As yet another alternative, procedures have been developed that favor a tonometric method that does require a blood pressure cuff. However, these methods still require some type of mechanical device that applies pressure to an artery, along with some other type of oscillation (pressure) sensor for the tonometric pressure estimation. Such a device is described in U.S. Pat. No. 6,730,038. Once again, devices of this type have proven unreliable during actual usage.
  • Other available methods for providing a continuous non-invasive blood pressure determination have tried to use formulas derived from the Bramwell-Hill equation. These methods utilize formulas that rely upon measured arterial pulse wave velocity (PWV) and measured arterial blood volume to determine blood pressure. In each of these methods, the arterial pulse wave velocity (PWV) and the arterial area must be measured. The required PWV measurements are typically obtained by observing the pulse transit time (PTT) between two widely separated sites, such as the heart and the finger tip. The pulse arrival times at the measurement sites, such as the finger tip, are typically determined by pleythysmography or pulse oximetry. One known method for determining PWV is described in U.S. Pat. No. 5,857,975. In this patent, the time of the pressure pulse's origin at the heart is determined from an EKG signal and the arrival time of the pulse is measured at another location on the patient. Based on these measurements, the instantaneous blood pressure is determined. In most systems, initialization data is obtained from a cuff based blood pressure determination and is used along with the time required for the pulse wave to travel between two points to calibrate formulas developed for continuously estimating blood pressure. After the initial calibration, changes in the pulse transit time can be related to changes in the blood pressure. In such schemes, the measurement of the area of the blood flow passageways, the blood flow, and PWV along with the subsequent blood pressure estimate are determined by the use of an assortment of complex adjustment factors.
  • In any type of non-invasive continuous blood pressure monitoring system, various factors can affect the accuracy of the measurement. For example, changes in the physiological state of the patient can bring about changes in the arterial wall elasticity. In general, changes in the arterial wall elasticity will affect the measured PWV. If the elastic modulus of the arterial wall changes, the same pressure may then need to be associated with a different cross-sectional area and PWV. If the operating point obtained using a blood pressure cuff calibration, then any changes in the arterial elasticity would require re-calibration, and if no re-calibration were performed, errors in the pressure estimation can clearly occur.
  • Another possible criticism of prior measurement systems concerns the measurement of the arterial blood volume. Specifically, the measurement of the arterial area by using pleythysmography is confounded by the highly elastic nature of the veins in the patient. Since the measurement includes some portion of the venous blood volume, it is difficult to rely upon the measurement to produce a total arterial lumen area. Venous blood volume is strongly affected by the subject position, since the hydrostatic pressure in the patient's body can cause pooling of the blood in highly distensible veins.
  • Yet another criticism of prior methods of making non-invasive blood pressure measurements concerns the manner in which the PWV measurements are made. Typically, any PWV measurement that is made between widely separated sites is measuring the PWV of a collection of branches of the arterial tree further complicating the relationship between PWV and blood pressure.
  • U.S. patent application Ser. No. 10/749,181, commonly assigned with the present application and incorporated herein by reference, teaches a method that utilizes pulse wave velocity to create a continuous, non-invasive blood pressure measurement. The '181 application teaches a method of measuring the pulse wave velocity within a patient and a method of relating such pulse transit time to blood pressure. Although the '181 application teaches a method of accurately measuring the PWV, the relationship to blood pressure can be criticized using some of the same grounds set forth above.
  • Thus, there is a desire to provide a more accurate, continuous, non-invasive blood pressure measurement technique and method.
  • SUMMARY OF THE INVENTION
  • The present invention is related to an apparatus and method for continuously and non-invasively monitoring the blood pressure that utilizes a transmission-line model and transcutaneous ultrasound measurements along a blood vessel of a patient. The resultant blood pressure measurements are displayed, in real time, to provide a continuous, non-invasive blood pressure measurement.
  • The apparatus of the invention includes at least a pair of ultrasound patch probes positioned along the skin of a patient. The ultrasound patch probes are separated by a relatively small distance and are positioned above the same blood vessel. The ultrasound patch probes are used by the system and apparatus of the present invention to make measurements relative to the blood vessel on a real time basis. Specifically, the ultrasound patch probes are used to determine the radius of the blood vessel and the blood flow rate at at least two distinct locations along the length of the blood vessel. The ultrasound area and flow measurements are specific to the particular artery being used, unlike pleythysmography.
  • The system of the present invention utilizes a time-domain model of a transmission line where the electrical parameters of capacitance, inductance and resistance of the transmission line model are used in analogy to blood flow characteristics. Specifically, the capacitance of the electric circuit is analogous to the compliance of the artery, the inductance of the electric circuit is analogous to the inertia of the blood and the resistance is analogous to the viscosity of blood flow. The input and output current sources in the electrical transmission line model is analogous to the blood flow rate, which are determined at both ends of the modeled arterial segment using ultrasound measurement techniques.
  • Based upon the determined values in the model, a transmission line simulator becomes a circuit representation of the artery, where the voltage across the capacitor of the transmission line model correlates to blood pressure. The system of the present invention uses Kalman filter techniques to provide the best estimates of the state variables, which include the voltage across the capacitor in the transmission line model, which is analogous to a blood pressure measurement.
  • The system and apparatus of the invention utilizes a blood pressure cuff to provide initial calibration measurements for the blood pressure analysis performed in real time. The initiation blood pressure information is used to develop an initial state vector that forms part of the Kalman equation estimates. The system monitors the Kalman filter results, as well as other system variables, and automatically recalibrates the system using measured parameters when the current results of the Kalman filtering or the system variables, as compared to the initiation values, exceed a pre-defined threshold value.
  • The continuous, non-invasive blood pressure monitoring technique utilizes the pair of ultrasound patch probes to make accurate, real time measurements of physical characteristics of the blood vessel through which the blood flow is used to make a blood pressure measurement. The ultrasound patch probes provide accurate data that is used to provide better estimates for the blood flow and size of the blood vessels through which the blood pressure is being determined. In this manner, the continuous, non-invasive blood pressure monitoring method provides greater accuracy as compared to prior methods.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The drawings illustrate the best mode presently contemplated for carrying out the invention:
  • In the drawings:
  • FIG. 1 is a illustration of the system of the present invention used to carry out the continuous, non-invasive blood pressure measurement technique;
  • FIG. 2 is a drawing showing a cross-sectional view of a typical cMUT cell;
  • FIG. 3 is a drawing showing the CNIBP measurement concept in accordance with the disclosed embodiment of the invention;
  • FIG. 4 is a schematic illustration of the transmission line model utilized to develop the blood pressure measurements of the present invention;
  • FIG. 5 a is a graphic illustration of simulated proximal and distal blood flow;
  • FIG. 5 b is a graphic illustration of reconstructed pressure and actual pressure in accordance with the present invention;
  • FIG. 6 is a graphic illustration of the error between estimated pressure and the known pressure with a prediction error; and
  • FIG. 7 is a high level data flow diagram illustrating the flow of data used to calculate the blood pressure estimates of the present invention.
  • DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
  • Referring first to FIG. 1, a blood pressure monitor 10, which may be part of a larger patient monitor 12 is shown. The patient monitor 12 may provide medical personnel with information regarding different physiological functions of a patient 14. Specifically, the patient monitor 12 may include an electrocardiograph (ECG), a pulse oximeter, and the blood pressure monitor 10.
  • The patient monitor 12 is connected to a central processor 16 that is used to carry out the continuous, non-invasive blood pressure monitoring functions to be described in greater detail below. The central processor 16 is coupled to a blood pressure cuff monitor 18, which in turn controls and monitors the readings from a conventional blood pressure cuff 20. The blood pressure cuff 20 is shown positioned around the arm of a patient and can be operated in the conventional manner to obtain actual systolic and diastolic blood pressure readings from the patient 14.
  • The central processor 16 is also coupled to an ultrasound activation circuit 22, which in turn is coupled to a first ultrasound patch probe 24 and a second ultrasound patch probe 26. As illustrated in FIG. 1, the first ultrasound patch probe 24 and the second ultrasound patch probe 26 are spaced from each other by a distance D along the forearm of the patient. Although the first and second ultrasound patch probes 24, 26 are shown positioned along the forearm of the patient 14, it should be understood that the ultrasound patch probes 24, 26 could be positioned at various other locations along the patient.
  • As illustrated in FIG. 1, each of the ultrasound patch probes 24, 26 is positioned on the arm of the patient in such a manner as to enable non-invasive insonification of an artery. Each of the ultrasound patch probes 24 is positioned generally above the radial artery 28 of the patient. Likewise, the blood pressure cuff is positioned to sense blood pressure in the brachial artery of the patient.
  • The use of non-invasive ultrasound can provide volumetric flow, arterial lumen area and pulse wave velocity information. In addition, ultrasound measurements can be taken in such a way that all the data describes a single, uniform arterial segment (e.g., the brachial or common carotid arteries can be used). This means that the measured PWV is directly related to the mechanical properties of the segment of elastic tube (artery) for which the blood volume is being measured. In a patient monitoring application, the operator of the ultrasound device is eliminated through the use of software that automatically locates the artery in the ultrasound data, e.g., using known edge detection techniques. Autonomous operation of the ultrasound system allows it to report blood pressure and blood flow traces to the clinical users without those users having to interpret an ultrasound image or operate an ultrasound imaging device.
  • The use of ultrasound measurements in the continuous non-invasive blood pressure monitoring application is enabled by micromachined ultrasonic transducer patch probe technology, which allows ultrasound data to be taken using a thin, lightweight probe that adheres to the patient's skin. Recently, semiconductor processes have been used to manufacture ultrasonic transducers of a type known as micromachined ultrasonic transducers (MUTs), which may be of the capacitive (cMUT) or piezoelectric (pMUT) variety. cMUTs are tiny diaphragm-like devices with electrodes that convert the sound vibration of a received ultrasound signal into a modulated capacitance. For transmission, the capacitive charge is modulated to vibrate the diaphragm of the device and thereby transmit a sound wave.
  • One advantage of MUTs is that they can be made using semiconductor fabrication processes, such as microfabrication processes grouped under the heading “micromachining”. As explained in U.S. Pat. No. 6,359,367:
      • Micromachining is the formation of microscopic structures using a combination or set of (A) Patterning tools (generally lithography such as projection-aligners or wafer-steppers), and (B) Deposition tools such as PVD (physical vapor deposition), CVD (chemical vapor deposition), LPCVD (low-pressure chemical vapor deposition), PECVD (plasma chemical vapor deposition), and (C) Etching tools such as wet-chemical etching, plasma-etching, ion-milling, sputter-etching or laser-etching. Micromachining is typically performed on substrates or wafers made of silicon, glass, sapphire or ceramic. Such substrates or wafers are generally very flat and smooth and have lateral dimensions in inches. They are usually processed as groups in cassettes as they travel from process tool to process tool. Each substrate can advantageously (but not necessarily) incorporate numerous copies of the product. There are two generic types of micromachining . . . 1) Bulk micromachining wherein the wafer or substrate has large portions of its thickness sculptured, and 2) Surface micromachining wherein the sculpturing is generally limited to the surface, and particularly to thin deposited films on the surface. The micromachining definition used herein includes the use of conventional or known micromachinable materials including silicon, sapphire, glass materials of all types, polymers (such as polyimide), polysilicon, silicon nitride, silicon oxynitride, thin film metals such as aluminum alloys, copper alloys and tungsten, spin-on-glasses (SOGs), implantable or diffused dopants and grown films such as silicon oxides and nitrides.
        The same definition of micromachining is adopted herein. The systems resulting from such micromachining processes are typically referred to as “micromachined electromechanical systems” (MEMS).
  • The use of a MUT patch allows a clinician to attach the transducer to the patient's skin. The MUT patch is lightweight and flat. For the purpose of illustration, the ultrasound patch probes will be described that are made up of capacitive micromachined ultrasonic transducers (cMUTs). However, it should be understood that the patches could instead employ pMUTs.
  • An embodiment will now be described that incorporates a MUT patch. However, it should be understood that the present invention encompasses not only a device, but also methods for continuous non-invasive blood pressure monitoring, and that the basic method is not limited to devices that employ MUT technology. Instead, the basic method disclosed herein can also be performed by appropriately programmed ultrasound imaging systems having probes comprising conventional piezoceramic transducer elements.
  • Based on acoustic data from the cMUT ultrasound patch probes 24 and 26 and pressure data from the blood pressure cuff 20, the digital signal processor (DSP) 16 autonomously computes various parameters, including the patient's estimated blood pressure. In addition, the processor software forms and steers ultrasound beams over the volume of space in front of the cMUT ultrasound patch probes 24 and 26. Using both echo intensity and Doppler blood flow information, the software maps out the location of the vessel, its cross-sectional area, and the area increase at the time of passage of the systolic pressure wave. The pulse wave velocity can also be measured. This will yield useful information about the mechanical properties of the blood vessel wall.
  • The blood pressure cuff 20 is periodically inflated by the central processor 16 to provide the peak and steady state pressures at fixed points in time. More specifically, the processor 16 activates a pump to inflate the blood pressure cuff 20 with air. A blood pressure cuff monitor 18 detects the pressure in the cuff 20. The analog output of the cuff 20 is converted to a digital value by the cuff monitor 18, which digital pressure value is used by the processor 16. Initially, the processor 16 triggers cuff inflation in response to an activation command input by the system operator via an operator interface (not shown). Thereafter, the processor 16 triggers re-inflation of the cuff 20 for re-calibration whenever the estimated compliance value changes by more than some predetermined amount from its last calibration level. The results of the computations are displayed on a screen (not shown) of the patient monitor 10.
  • Although the pair of ultrasound patch probes 24, 26 and the blood pressure cuff 20 pneumatics in the disclosed embodiment are controlled and monitored by electronics coupled to the processor 16 by means of cables, obviously a wireless electrical coupling could be provided.
  • Referring to FIG. 2, a typical cMUT transducer cell 2 is shown in cross section. An array of such cMUT transducer cells is typically fabricated on a substrate 4, such as a heavily doped silicon (hence, semiconductive) wafer. For each cMUT transducer cell, a thin membrane or diaphragm 8, which may be made of silicon nitride, is suspended above the substrate 4. The membrane 8 is supported on its periphery by an insulating support 6, which may be made of silicon oxide or silicon nitride. The cavity 9 between the membrane 8 and the substrate 4 may be air- or gas-filled or wholly or partially evacuated. A film or layer of conductive material, such as aluminum alloy or other suitable conductive material, forms an electrode 11 on the membrane 8, and another film or layer made of conductive material forms an electrode 13 on the substrate 4. Alternatively, the bottom electrode can be formed by appropriate doping of the semiconductive substrate 4.
  • The two electrodes 11 and 13, separated by the cavity 7, form a capacitance. When an impinging acoustic signal causes the membrane 8 to vibrate, the variation in the capacitance can be detected using associated electronics (not shown in FIG. 2), thereby transducing the acoustic signal into an electrical signal. Conversely, an AC signal applied to one of the electrodes will modulate the charge on the electrode, which in turn causes a modulation in the capacitive force between the electrodes, the latter causing the diaphragm to move and thereby transmit an acoustic signal.
  • For the CNIBP application, two types of data are acquired simultaneously using the cMUT patch, as shown in FIG. 3. The first data type is a sequence of B-mode images providing several short-axis views of the artery (in planes perpendicular to the direction of blood flow). Two or more short-axis views are taken along the length of the artery. The center of the artery can be found at each measurement location using edge detection techniques. The center of the artery can then be estimated for any location on the patch using interpolation. In addition to finding the artery center, the B-mode data is used to estimate the vessel diameter for the blood pressure calculation. More specifically, it is assumed that that the cross section of the lumen is a circle, so the estimated diameter can be used to calculate the lumen area A.
  • The second data type is a set of M-mode lines to be used for Doppler processing. The M-mode lines remain locked onto the center of the artery using information from the B-mode data. The multiple M-mode lines can be used with no Doppler steering angle for determining the velocity of the vessel wall or with a Doppler steering angle for determining blood velocities. The M-mode data, either wall or blood velocities, will be used for estimating the pulse wave velocity needed for the blood pressure calculation.
  • The two types of data are acquired by interleaving the acquisitions of the short-axis B-mode beams and the Doppler M-mode beams. A repeating scan sequence consists of a single B-mode beam at each location and all the M-mode beams. In the next pass through the sequence, the B-mode beams are translated across the patch while the M-mode beams remain in the same position. This is repeated until the desired width for the B-mode image is covered, generating a single frame of the B-mode data. Then the steps are repeated using any updated information about the vessel center. The result of the acquisition is M-mode data at a high repetition rate (˜2 kHz), as required for Doppler processing, and the B-mode images at a slower rate (˜100 Hz) determined by the required image size for vessel tracking.
  • The B-mode beams are translated across the cMUT patch by translating an annular array of activated sub-elements across the patch. A uniform translation of the beamforming coefficients produces a new beam at a different location. Repeated frequently, this generates a rectilinear two-dimensional image.
  • An additional bilinear term in the beamforming coefficients produces a beam directed away from the normal. Such beams are transmitted in the M mode to acquire data for Doppler processing that computes the blood velocity based on the Doppler Effect.
  • M-mode beams that are transmitted generally normal to the artery are used to acquire data for Doppler processing of the arterial wall velocity during deflection of the wall as the pulse wave travels down the artery. The maximum deflection occurs when the wall velocity is zero.
  • Referring back to FIG. 1, each of the first and second ultrasound patch probes 24, 26 can be actuated by the activation circuit 22 to emit an ultrasound pulse 30. The ultrasound pulse 30 allows the ultrasound activation circuit 22 to make multiple measurements and determinations relative to a uniform stretch of an artery, such as the brachial artery 28 positioned between each of the ultrasound patch probes 24, 26. In an alternate method of the present invention, the patch probes 24, 26 can be spaced on a range of artery in which there would be several bifurcations. In the embodiment of the invention illustrated, the ultrasound patch probes 24, 26 are shown spaced over a length of the unbifurcated brachial artery. Specifically, the central processors 16 can determine the internal diameter of the artery 28 at the location beneath the first ultrasound patch probe 24 and at the location beneath the second ultrasound patch probe 26. In addition to determining the diameter of the artery 28, the central processor 16 can utilize the Doppler ultrasound flow measurement technique to determine the instantaneous flow of blood at a first location beneath the first ultrasound patch probe 24 and at a second location beneath the second ultrasound patch probe 26. Thus, the central processor 16, through the ultrasound activation circuit 22 and the pair of ultrasound patch probes 24, 26 can determine the arterial diameter and blood flow rate at two separate locations along the radial artery which are spaced by the distance D. These two measurements made through the use of the ultrasound patch probes 24, 26 provide important and accurate information to the central processor 16 that aids in the determination of a continuous, non-invasive blood pressure measurement as described in detail below.
  • The data collection scheme associated with the present invention requires that at least two ultrasound probes be deployed along the length of a short arterial segment. In the monitoring application, it is anticipated that these will be the ultrasound patch probes 24 and 26. Each of these probes must have a steerable beam and must be able to operate at an angle with respect to the skin normal so that it can collect Doppler blood flow measurements. Such a probe can also be steered to the normal angle and used to measure the area of the artery's lumen. Optionally, other probes may be deployed between the two end probes, all of them at known locations along the line of the arterial segment.
  • Referring now to FIG. 4, the blood pressure monitoring system and method of the present invention utilizes an algorithm in which a length of an unbifurcated artery is represented by a transmission line model 32. The transmission line model 32 shown in FIG. 4 is a well known model and technique used for determining the phase velocity of electrical waves propagating in a cable or transmission line. In accordance with the present invention, the transmission line model 32 can be analogized to the flow of blood in a distensible tube, and thus used to make blood pressure measurement. In the present approach, the transmission line is an analog for such a tube, with electric charge corresponding to fluid volume, current corresponding to volumetric flow and voltage corresponding to pressure. Thus, the specified flow measurements are like direct measurements of current in the analogous electric circuit.
  • In accordance with the present invention, the transmission line model 32 will be represented by a ladder network, each stage of which represents a short sub-segment of the arterial segment. The operation of this network can be represented by a set of simultaneous, first-order differential equations, which can be used to write a discrete-time state variable model for the system, given a sample rate. In the state variable model, half of the state variables will be pressures, the other half will be flows. This model can then be driven by ultrasound measurements and the pressure estimate will be given by one of the pressure state variables. Naturally, the model must have parameters that accurately represent the physical system in order to give meaningful results.
  • The general method of allowing a state variable model to evolve in time is reminiscent of the numerical solution of an initial value problem using Euler's method. The present apparatus approach differs from a typical initial value problem in that the equation to be solved is inhomogeneous, so if the system is stable, the particular solution associated with the driving function dominates the general solution. It is for this reason that the initial values of the state variables are ignored, except as described below in relation to the operating pressure.
  • The transmission line model 32 is represented in the signal processing method as a RLC lossless ladder network shown in FIG. 4. In the present application, the term “lossless” indicates that no mass (i.e. blood) is lost from the input to the output of the network, such as in a length of artery without any branching. The analogous quantity to the inductance L of the ladder network is the mass per unit length (per ladder network section) and the analogous quantity to the capacitance C is the elastic compliance of the tube. The resistors R represent energy loss due to the viscosity of blood. In the simplest version, the network is driven at both ends by current sources 44 and 46, representing the ultrasonic blood flow measurements.
  • Although the RLC ladder network in FIG. 4 is described as being a “lossless” network, it should be understood that if the method of the present invention is used on more complex artery structures, such as the radial artery, or a longer arterial segment, some of the blood mass is diverted in branching. In such an embodiment, a transverse resistance, in parallel with each capacitance, can be included in the ladder network. Although the present method does not use ultrasound measurements on the branching arteries, a relationship could be developed between changes in the dimensions and flow of the artery segments being measured and the branches that are not being measured. This relationship can be used to estimate the transverse resistance dynamically. Therefore, the presence of a transverse resistance can be handled within the scope of the described technique.
  • In the embodiment illustrated in FIG. 4, the transmission line model includes multiple RLC ladder network stages that each include an inductance 36, a capacitance 38 and a series resistance 40. The ladder network 34 includes n different stages.
  • In practice, the order of the transmission line model, which is related to the number of serial sections employed, is to be determined by the length of the arterial segment to be modeled, and the parameters represent values of an underlying distributed parameter, expressed per length of the individual section. For example, if the distance between the proximal and distal flow measurements is 5 cm (which is typical), then a 10-section transmission line model would have its compliance value given in units equivalent to “compliance per 5 mm section of tube”. When the length of the individual section is short enough, then the compliance value in the model gives an approximate relationship between area and pressure. This corresponds to the definition of compliance used, for example, in the, Bramwell-Hill equation, discussed below.
  • The operation of the network of FIG. 4 is expressed in the following set of simultaneous ordinary differential equations:
    Q P −Q 1 =C dP 1 /dt
    P 1 −P 2 =RQ 1 +L dQ 1 /dt
    Q 1 −Q 2 =C dP 2 /dt
    P 2 −P 3 =RQ 2 +L dQ 2 /dt
    Q 2 −Q D =C dP 3 /dt  (1)
  • In the transmission line model, the state variables are the current through each of the inductors 36 as well as the voltage across each of the capacitors 38. In the present invention, the ground 48 of the transmission line model 32 is assumed to be zero such that the voltage across capacitor 38 represents the voltage at node 50. As previously indicated, the voltage in the transmission line model corresponds to pressure in the arterial transmission line model. (The ground value should actually be equal to the pressure in the venous system, but this is so low compared to the arterial pressures that it is negligible.) Thus, solving for the state variables and specifically the voltage across the capacitors 38, results in a value that is the instantaneous blood pressure estimate desired.
  • In matrix form, these equations are given by [ P . 1 Q . 1 P . 2 Q . 2 P . 3 ] = [ 0 - 1 C 0 0 0 1 L - R L - 1 L 0 0 0 1 C 0 - 1 C 0 0 0 1 L - R L - 1 L 0 0 0 1 C 0 ] [ P 1 Q 1 P 2 Q 2 P 3 ] + [ 1 C 0 0 0 0 0 0 0 0 - 1 C ] [ Q p D ] ( 2 )
    where we use the notation P to denote the first derivative of P with respect to time. Defining the vector of state variables x=[P1 Q1 P2 Q2 P3 ]T, we may write this as x . = Ax + B [ Q P Q D ] ( 3 )
    where the matrices A and B are defined conventionally. Equation (3) gives a continuous-time expression for the transmission line model. The signal processing operation operates in discrete time, since the proposed measurements are intrinsically discrete time. In order to convert Equation (3) into an equivalent discrete-time model, a new state feedback matrix G is defined, which is the state transition matrix of the continuous system over the sample period. An input matrix F is defined that is equivalent to the effect of a constant input signal over the sample period: x ( k + 1 ) = AT x ( k ) + ( 0 T At t ) B [ Q p ( k ) Q D ( k ) ] = Gx ( k ) + F [ Q p ( k ) Q D ( k ) ] ( 4 )
  • This well-known-approach to discretization of a continuous time linear system represents a delay-free, zero-order hold assumption for the digitization of the driving signals. This is equivalent to representing the driving input signal by a piecewise constant approximation. Other digitization models can also be employed, such as a first-order hold, which represents the input signal as piecewise linear. Higher order hold models may be required at lower sample rates. Although the effects of quantization of the flow measurements are not included in this model, these effects can simply be considered as part of the additive observation noise.
  • The basic operation of the ladder network of the transmission line model is to integrate the two flow signals to get a time-varying volume estimate, and to scale that using the compliance (C) to obtain a pressure (P) estimate. The relationship between the cross-sectional area (A) of the lumen and the pressure is simply A=CP, and in a uniform tube, the volume (V) is simply Al, for a segment of length l. The series resistors 40, representing blood viscosity, are very important in implementing the low-pass frequency response characteristic that performs the integration.
  • In addition, the network implements a delay for the flow signals. Because of this delay, blood observed entering the system at a given time at one end will be canceled by blood leaving the system at a later time at the other end, with the interval determined by the implemented delay. This is true for the forward flow wave and also for the reflected, or retrograde, flow wave. When the delay does not match the actual delay in the system, the flow is cancelled either too soon or too late, with the result that the pressures are wrong. The phase velocity of the network depends on the product LC.
  • Most of the parameters of the ladder network 32 can be obtained from ultrasound measurements. The compliance is defined as C = A P = 2 π r 3 Ea ( 5 )
    where E is Young's modulus, a is the thickness of the arterial wall and r is the radius of the lumen. Although the product Ea cannot be measured directly using ultrasound, a quantity that is parameterized by it can be measured. Recently, a method for estimation of the arterial pulse wave velocity (PWV) that compensates for the effects of pulse reflections has been developed, as described in U.S. Patent application Ser. No. 10/749,181. This method can be used to determine the characteristic pulse wave velocity in a short arterial segment. We define the characteristic PWV as that velocity that would be observed in the absence of reflections. For this velocity, we have the Bramwell-Hill equation: v p = π r 2 P ρ A ( 6 )
    where ρ p is the density of blood. Substituting (5) into (6) and solving for C, we get C = π r 2 ρ v p 2 ( 7 )
    so that the C parameter can be measured using PWV, arterial radius and the known density of blood. Note that these equations can also be solved for the Ea product. The blood mass/inertia parameter, L, is given by L = ρ π r 2 ( 8 )
    this can be measured using only the arterial lumen radius. The serial resistance is given by: R = 8 η π r 4 ( 9 )
    where η is the viscosity of blood. Because the viscosity cannot be directly measured by ultrasound, we can use a nominal constant for the viscosity. Alternatively, an adaptive process, described in detail below, can be used to adjust the viscosity until the model is correct.
  • Since there are no pressure sources in the transmission line model 32 of FIG. 4, and since the pressures enter into the equations only by way of first derivatives or of differences between two pressures, any constant could be added to the pressure variables without affecting the operation of the circuit, since the constant will be cancelled in every equation. Further, once a constant is added to the pressures in Equation 1, there is nothing to make it change, since there is no resistor to ground at any stage of the circuit. In practice, the mean pressure level will be set by a blood pressure cuff measurement. In a comparison of the measurements calculated using the present invention to simulated data, it has been found that the pressure state variables should be initialized to be slightly less than the mean pressure. Alternatively, several cycles of pressure observations can be run through the system, and the pressure state variables adjusted so that the mean pressure is that measured by the blood pressure cuff; this adjustment is done by the addition of the same quantity to all the state variables. If, at a later time, the compliance parameter changes, the mean pressure will change in response, since P=A/C, so that the same integrated flow will yield a different pressure. Such a change in mean pressure can also be used to trigger a recalibration using the cuff.
  • As a means of showing the output waveforms that will be provided by this technique, a simulation was performed by taking a smoothed version of an impedance pleythysmography measurement from the femoral artery of a human subject, and taking its first difference to simulate a forward flow at the proximal measurement site. (Impedance pleythysmography measures the total blood volume at the measurement site.) The simulated proximal forward flow was delayed to get the forward flow at the distal site. A reflected flow wave was then added to both. The reflected wave as modeled as originating at a point 2 centimeters distal to the distal measurement site that reflected the wave with a reflection coefficient of 0.3. Since flow waves add algebraically, the reverse flow wave subtracts from the forward wave. Two cycles of the simulated total proximal flow 74 and distal flow 76 are depicted in FIG. 5 a. Note that the distal flow appears more attenuated by the reflected flow because the reflection effect occurs sooner at the distal site.
  • FIG. 5 b shows the result of inputting the flow waveforms on FIG. 5 a in to the model above. The pressure variables of the model were initialized to a value of 70. The reconstructed pressure 78 is generally well-estimated, as compared to the actual pressure 80. The parameters of the model and the scale of the flow measurement were obtained by searching until a correct representation of the flow waveform occurred. Note that the simulation methodology did not enforce a consistent relationship between the delay, frequency response and scaling of the pressure variables. This example has been provided only to show the capability of the technique in producing a continuous non-invasive blood pressure waveform.
  • When additional measurements located between the two ultrasound patch probes 24 and 26 are available, it is not immediately clear how they should be used in the estimation scheme. What is clear is that they cannot be made sources, since that would have the effect of decoupling the sections of the model proximal and distal to the new source.
  • The answer to this dilemma is to turn the estimation procedure into a Kalman filter and to use the new observations as measurements, while retaining the original two measurements as system inputs. Using this approach, both flow and distension, or area, measurements can be incorporated, as available. Each measurement must correspond to a state variable that best represents its distance from the two ends of the model, defined by the two driving flow measurements.
  • The Kalman filter algorithm is composed of three main operations, which together recursively update the system state and its covariance matrix, which in turn determines the relative weight given to the observations as opposed to that given to the current state estimates. In the general form of the Kalman filter, the state update mechanism is taken to occur with additive noise. In accordance with the present invention, process noise is omitted from the model, or equivalently, the covariance matrix is set to zero.
  • In order to accommodate the new, interior observation sites, an observation matrix H is designed that specifies the observed state variables. In general, H has the same number of columns as the dimension of the state vector and the number of rows is equal to the number of additional observations. In each row, there is only one non-zero entry, which is unity if the measured state variable is a flow and C if the state variable is a pressure and the measurement is an area. The observation is then defined as
    z(k)=Hx(k)+v k  (10)
    where vk is an observation noise vector with covariance matrix R=diag(n). The components of z(k) are either noisy flow measurements or noisy area measurements.
  • At the start of every iteration of the Kalman algorithm, it is assumed that an estimate of the current state vector denoted by {circumflex over (x)}(k), is known, and its (time-varying) covariance matrix, denoted P(k), is also known. On the initial iteration, P(1) is set to a diagonal matrix with large entries, to represent the initial uncertainty in the state estimate. From the state and observation covariance matrices, the Kalman gain is computed, given by
    K(k)=P(k)H T(HP(k)H T +R)−1  (11)
  • The innovations sequence is simply the difference between the observation and the observation that would be predicted by the current state vector. The Kalman gain matrix computed in Equation 11 is used to weight the innovations sequence to produce an updated estimate of the state vector at the kth sample time.
    {tilde over (x)}(k)={circumflex over (x)}(k)+K(k)(z(k)−H{circumflex over (x)}(k))  (12)
    Note that this equation replaces the current state vector by a weighted sum of itself and the observation.
  • The third and last step in the Kalman algorithm is to project the state vector and its covariance matrix to the next sample time: x ^ ( k + 1 ) = G x ~ ( k ) + F [ Q p ( k ) Q D ( k ) ] ( 13 a ) P ( k + 1 ) = G [ ( I - K ( k ) H ) P ( k ) ] G T ( 13 b )
  • The effect of the Kalman filter is to replace the observed state variables at every step with a weighted average of their observed and predicted values. In the steady state, the weights depend to a large extent on the value of the noise variance on the main diagonal of R.
  • It is well-known that the actual compliance of arteries is not constant. The artery compliance is the slope of the transfer characteristic between pressure and the cross-sectional area of the vessel lumen. This characteristic is nonlinear because of the collagen fibers in the arterial wall, which are relaxed until a certain stage of the distension of the vessel, and then become taut, decreasing the elasticity. Thus, the compliance can be considered to be dependent on pressure. This nonlinear relationship can be represented by a piecewise linear relationship such that the compliance is different in a number of disjoint pressure ranges. This kind of pressure-dependent compliance can easily be incorporated into the time-domain approach above.
  • The pressure/area relationship has been modeled in the prior literature as an exponential. This transfer function was applied to a Windkessel model of the entire arterial system as seen from the aorta, and was found to have some advantages in terms of predicting the shape of observed pressure waveforms.
  • In the present model, the compliance enters into the definition of G, F and H. In order to implement pressure-dependent compliance, a different set of matrices must be implemented for every distinct pressure interval over which the compliance is constant. Equations (11), (12) and (13) are then executed using that set of matrices that is associated with the estimated pressure at step (k-1).
  • There are several features of flow estimation using ultrasound that must be kept in mind when implementing the method of the present invention. Volumetric flow estimation proceeds in the same way as ultrasound color flow mapping over a limited sector. Color flow mapping is a standard feature of most diagnostic ultrasound devices that estimates blood velocity at a number of locations on the image plane by computing the mean frequency of the spectrum of the returned ultrasound signal. This mean frequency corresponds to a mean Doppler shift induced by the moving blood. This measurement process involves insonifying a volume from which a blood velocity is desired, and this occurs from three to sixteen times for every estimate. This is done along several ultrasound A-lines. Once a set of blood velocity estimates is available at a number of locations covering the vessel lumen, the estimates are all multiplied by the area to which they individually apply and the products are summed. This implements an approximate integration over the lumen and the result is a volumetric flow.
  • However, ultrasound flow mapping can produce blood velocity estimates that are highly biased down from the true blood velocity if the ultrasound return contains significant energy from stationary tissue reflections. For this reason, the data used to compute the mean frequency estimates (the Doppler shifts) is typically high-pass filtered. This means that very low flows, such as might be encountered during diastole in a peripheral artery, may be missed. Since the majority of the flow undoubtedly occurs when the blood velocity is high, this may not be a major problem, but still, some portion of the flow is bound to be neglected due to this effect.
  • Another effect that can change the performance of the blood pressure measurement is the dependence of the ultrasound velocity measurement on angle. The measured phase change is indicative of the projection of the blood velocity vector on the ultrasound line. If the various flow measurements are based on different Doppler angles, they can be incommensurate. The lossless circuit of FIG. 4 requires that positive current entering the system at the proximal source 44 be canceled by negative current entering the system at the distal source 46 after some delay. If the angles of the two ultrasound lines are different, this cancellation will not take place, and the blood pressure will not be correctly computed. Note that if the ultrasound probes 24, 26 are mechanically connected to each other by a rigid structure, identical insonification angles are easily obtained.
  • No matter how many ultrasound lines are used to produce the velocity estimates, a new estimate can be computed for every new line, even though it will be highly correlated with the last such estimate. This means that a new volumetric flow estimate can be obtained for every frame. That is, the frame rate and volumetric flow sample rate are the same. The frame rate should be maximized by using only the minimum number of ultrasound lines required to cover the artery.
  • In the above discussion of the parameters of the transmission line model, it was noted that not all of the important parameters were directly observable using ultrasound. In particular, the product of the modulus of elasticity and the wall thickness (Ea), and the blood viscosity, represent two degrees of freedom not directly observable. We noted that the (Ea) product can be computed from pulse wave velocity, as long as the characteristic velocity rather than the apparent velocity, can be obtained, and that a method described in U.S. patent application Ser. No. 10/749,181 is available to obtain such velocity.
  • Alternatively, a data-adaptive processing can be utilized to find the undetermined free parameters of the model. This method, however, requires that at least one additional measurement be made, over and above the two required flow measurements described previously.
  • Many adaptive filters operate by minimizing a prediction error energy. In general, a prediction error is the error between the response of the filter to a known input and the response of the real system to that input. This is a well-known principle in spectrum estimation and system identification. In accordance with the invention, this approach is employed to indicate the degree of agreement between the parameters of the transmission line model and the state of the physical system.
  • For this purpose we use the model of the observations defined in Equation (10) in connection with the Kalman filter approach. The prediction error power at the kth sample is given by the squared norm of the kth element of the innovations sequence, and the prediction error energy is proportional to the prediction error power summed over some fixed period of time: E p = T s k = k 0 k 1 z ( k ) - H x ^ ( k ) 2 ( 14 )
    where TS is the sample period.
  • The adaptive filter works by applying multiple sets of parameters to the data in parallel over a fixed period of time, and computing EP for each. The set of parameters producing the lowest prediction error energy is taken to be the correct set of parameters. The number of such parallel computations will obviously be limited by computational resources; however, even as few as three parallel computations will allow the adjustment of a single parameter. If one parameter set is chosen to be above the nominal parameter value, one chosen to be below it and another to be exactly at the nominal value, then after each observation period, the new nominal value can be that one of the three with the lowest EP, and the process can be repeated. This allows the system to track a time-varying parameter with minimal computational resources.
  • The parameters that will typically be varied in this adaptation will be the (Ea) product and the blood viscosity, since the L parameter and the other components of the C and R parameters can be measured using ultrasound. Measurement of the characteristic pulse wave velocity of the artery can also be used to help estimate C. It is also possible to limit adaptation of these parameters to a region around the nominal value of η and the value of C predicted by the Bramwell-Hill equation.
  • As an example of the operation of the adaptive filter, the simulation example shown in FIGS. 5 a and 5 b above was repeated with a third, central, flow measurement, and with a range of compliance values around the optimum. The squared error between the known correct pressure and the predicted pressure, given knowledge of the correct mean pressure, was computed, along with the prediction error energy for predicting the third flow. The results are depicted in FIG. 6, in which the squared estimation error between the estimated pressure and the known pressure (line 82) is plotted against the prediction error (line 84). The two curves are in different units and they have been scaled and shifted so that they could be plotted on the same graph. The minimum value of each is marked by a circle 86 and 88. The plot shows that the prediction error curve 84 is smooth, and that it can be made to have the same minimum as the actual estimation error 84, provided the non-adapted parameters are correct.
  • Adaptation of the system parameters using prediction error is consistent with incorporation of the additional measurements using the Kalman filter. All that is required is the execution of the operations given by equations (11), (12) and (13b) in addition to the normal operation given above. There is a potential problem in this approach, however. To the extent that the observation used to generate prediction error updates the state variable to which it is compared, the sensitivity of the prediction error to model errors will be reduced. The degree to which the observations update the state variables is controlled by the noise variance specified in the R matrix; thus, adaptive Kalman filter operation will require that the adaptation occur during fixed periods, during which the noise variance is set to a high level. In periods when the adaptation is not taking place, the noise variance should be set back to a realistic level.
  • Referring now to FIG. 7, there shown is the flow of data utilized for the continuous, non-invasive blood pressure measurement in accordance with the present invention. As shown on the left side of FIG. 7, the system receives a brachial artery diameter measurement in step 60. The artery diameter measurement is an ultrasound measurement taken from both of the ultrasound patch probes 24, 26 illustrated in FIG. 1. In addition to the diameter measurements, the input blood flow at ultrasound patch probe 24 is measured using ultrasound techniques. Likewise, the output blood flow is measured at block 64 by the ultrasound patch probe 26. The three ultrasound input measurements in blocks 60, 62 and 64 are all transferred to a transmission line simulator shown at block 66. The transmission line simulator 66 utilizes the input measurements, in connection with additional values to be described below, to carry out the calculation steps previously described. The transmission line simulator 66 includes the Kalman filtering as well as the state variable formulas to generate the blood pressure estimates, which are displayed in block 68.
  • The data flow also includes a pulse wave velocity estimator, which uses known measurement techniques, such as those taught in U.S. patent application Ser. No. 10/749,181 to calculate the pulse wave velocity. The pulse wave velocity estimator 70 also utilizes the input and output blood flow measurements from blocks 62, 64 to make an accurate pulse wave velocity estimate.
  • The data flow also includes an elastic estimator 72. The elastic estimator 72 utilizes the pulse wave velocity measurement from step 70 and an empirical formula to calculate the elasticity of the blood vessel through which the blood pressure is being determined. As illustrated in FIG. 7, the pulse wave velocity and elastic estimates are sent to the transmission line simulator 66 such that the transmission line simulator 66 can carry out the calculations to determine the blood pressure estimate.
  • In the present embodiment of the invention described above, the method and apparatus are described as presenting a continuous blood pressure measurement for a patient using ultrasound patch probes placed on the skin of the patient. Typically, this continuous measurement will be represented as a waveform that can be displayed. However, the method of the present invention is equally applicable for estimating the systolic and diastolic blood pressure for each beat of the patient's heart. The beat by beat estimate of the systolic and diastolic pressure would provide a continuous blood pressure measurement as the term is used herein.

Claims (27)

1. A method for continuously and non-invasively monitoring blood pressure of a patient, the method comprising the steps of:
positioning at least two non-invasive ultrasound sensors aligned with a blood vessel of the patient, the ultrasound sensors being operable to determine the diameter of the blood vessel and the flow rate of blood through the blood vessel;
utilizing a transmission line model to represent the blood vessel, where the voltage at each node of the transmission line model corresponds to the blood pressure; and
using a transmission-line-based computation to estimate the blood pressure.
2. The method of claim 1 wherein the two ultrasound sensors are each aligned with the same blood vessel, the ultrasound sensors being spaced by a distance.
3. The method of claim 2 wherein the transmission line model includes an input current source and an output current source positioned between an RLC network, wherein the input current source is a first blood flow rate determined by the first ultrasound sensor and the output current source is a second blood flow rate determined by the second ultrasound sensor.
4. The method of claim 3 wherein the RLC ladder network includes an inductance (L), a capacitance (C), and a series resistance (R), wherein the capacitance is analogous to the compliance of the blood vessel, the inductance is analogous to the density of the blood divided by the cross-sectional area of the blood vessel and the series resistance is analogous to eight times the viscosity of blood divided by π times the radius of the lumen of the blood vessel to the fourth power.
5. The method of claim 4 wherein the RLC ladder network further includes a transverse resistance corresponding to loss of blood mass due to branching from the blood vessel.
6. The method of claim 1 further comprising the step of initializing the measurement system by determining a measured blood pressure of the patient by utilizing a blood pressure cuff, the measured blood pressure being utilized to determine the pulse wave velocity.
7. The method of claim 3 wherein the step of using the transmission-line based computation includes the steps of:
developing a series of state variable equations based on the transmission line network;
inputting the measured first blood flow rate, the measured second blood flow rate and a measured blood vessel radius into a Kalman filter; and
utilizing the Kalman filter to estimate the state variables of the transmission line model, wherein the state variables include a capacitor voltage for each node in the transmission line model.
8. The method of claim 7 further comprising the steps of:
determining the pulse wave velocity for the blood flowing through the blood vessel; and
determining the capacitance (C) of the RLC circuit by
C = π r 2 ρ v p 2
where r is the measured radius of the blood vessel, ρ is the density of the blood and vp is the pulse wave velocity of the blood.
9. The method of claim 7 further comprising the step of utilizing the initialization measurement of the blood pressure to define an initial state vector used in the Kalman filter.
10. The method of claim 7 further comprising the steps of:
comparing a current value of the state variables with an initialized value of the state variables; and
re-initializing the measurement system when the current value of the state variables differs from the initialized value of the state variables by more than a pre-determined threshold value.
11. The method of claim 7 further comprising the steps of:
continuously determining the radius of the artery of the patient;
continuously determining the first blood flow rate and the second blood flow rate; and
continuously updating the state variable equation with the blood vessel radius, first blood flow rate and second blood flow rate, such that the Kalman filter provides a continuous, real time blood pressure measurement.
12. The method of claim 7 wherein the state variable equations include at least one undetermined parameter, the method further comprising the steps of:
providing multiple, separate estimates for the value of the undetermined parameter;
calculating a prediction error for each of the estimates, wherein the prediction error is the difference between the response of the Kalman filter for the estimate and a measured blood pressure of the patient;
selecting the estimate of the undetermined parameter that results in the smallest prediction error; and
utilizing the selected estimate with the series of state variable equations and the Kalman filter to determine the continuous blood pressure.
13. A method for continuously and non-invasively monitoring the blood pressure of a patient, the method comprising the steps of:
utilizing an ultrasound measurement to determine the diameter of a blood vessel of a patient;
utilizing an ultrasound measurement to determine a first blood flow rate at a first location along the blood vessel;
utilizing an ultrasound measurement to determine a second blood flow rate at a second location along the blood vessel, the second location being spaced from the first location;
estimating the pulse wave velocity along the blood vessel;
creating a transmission line model for the blood vessel, the transmission line model having multiple RLC ladder stages extending between a first current source and a second current source, where the first current source corresponds to the first blood flow rate and the second current source corresponds to the second blood flow rate; and
using a transmission-line-based computation to estimate a voltage at each node of the transmission line model, wherein the voltage corresponds to the blood pressure.
14. The method of claim 13 wherein each RLC ladder stage of the transmission line model includes a capacitance (C) analogous to the compliance of the blood vessel, an inductance (L) analogous to the density of the blood divided by the area of the blood vessel and a series resistance (R) analogous to eight times the viscosity of blood divided by π times the radius of the lumen of the blood vessel to the fourth power.
15. The method of claim 13 wherein the diameter of the blood vessel, the first blood flow rate and the second blood flow rate are determined by at least two ultrasound sensors positioned on the skin of the patient and aligned with the blood vessel.
16. The method of claim 13 further comprising the step of:
initializing the measurement system by determining an initialized value of the blood pressure of the patient by utilizing a blood pressure cuff, the initialized value of the blood pressure being used to solve the transmission line model.
17. The method of claim 13 wherein the step of using the transmission-line-based computation to estimate the voltage at each node includes the steps of:
developing a series of state variable equations based on the transmission line network;
inputting the measured first blood flow rate, the measured second blood flow rate and the measured blood vessel diameter into a Kalman filter; and
utilizing the Kalman filter to estimate the state variables of the transmission line model, wherein the state variables include a capacitor voltage for each node in the transmission line model.
18. The method of claim 17 further comprising the steps of:
determining the pulse wave velocity for the blood flowing through the blood vessel; and
determining the capacitance of the RLC circuit by
C = π r 2 ρ v p 2
where r is the measured radius of the blood vessel, ρ is the density of the blood and vp is the pulse wave velocity of the blood.
19. The method of claim 17 further comprising the step of utilizing the initialization measurement of the blood pressure to define an initial state vector used in the Kalman filter.
20. The method of claim 17 further comprising the steps of:
comparing a current value of the state variables with an initialized value of the state variables; and
re-initializing the measurement system when the current value of the state variables differs from the initialized value of the state variables by more than a pre-determined threshold value.
21. The method of claim 13 further comprising the steps of:
continuously determining the radius of the blood vessel of the patient;
continuously determining the first blood flow rate and the second blood flow rate; and
continuously updating the state variable equation with the blood vessel radius, first blood flow rate and second blood flow rate, such that the Kalman filter provides a continuous, real time blood pressure measurement.
22. The method of claim 21 wherein the real time blood pressure measurement is presented as a systolic and a diastolic pressure on a beat by beat basis.
23. The method of claim 17 wherein the state variable equations include at least one undetermined parameter, the method further comprising the steps of:
providing multiple, separate estimates for the value of the undetermined parameter;
calculating a prediction error for each of the estimates, wherein the prediction error is the difference between the response of the Kalman filter for the estimate and a measured blood pressure of the patient;
selecting the estimate of the undetermined parameter that results in the smallest prediction error; and
utilizing the selected estimate with the series of state variable equations and the Kalman filter to determine the continuous blood pressure.
24. A continuous, non-invasive blood pressure monitoring system comprising:
at least a pair of non-invasive ultrasound sensors each being operable to determine the diameter of the blood vessel and the flow rate of blood through the blood vessel;
a blood pressure cuff operable to determine an initialization value of the blood pressure of the patient; and
a processor connected to both the ultrasound sensors and the blood pressure cuff for determining the blood pressure of the patient continuously and non-invasively, wherein the processor utilizes a transmission line model to estimate the blood pressure.
25. The apparatus of claim 24 wherein the transmission line model utilized by the processor includes multiple RLC ladder stages extending between a first current source and a second current source, where the first current source corresponds to the first blood flow rate measured by the ultrasound sensor and the second current source corresponds to the second blood flow rate, where the processor develops a series of state variable equations based on the transmission line network and utilizes Kalman filtering to estimate the state variables.
26. The apparatus of claim 24 wherein the capacitance (C) in the transmission line model is analogous to the compliance of the blood vessel, the inductance (L) is analogous to the density of the blood divided by the cross-sectional area of the blood vessel and the resistance (R) is analogous to eight times the viscosity of blood divided by π times the radius of the lumen of the blood vessel to the fourth power.
27. The apparatus of claim 24 wherein the processor compares the current value of the state variables with an initializing value of the state variables and re-initializes the measurement system when the current value of the state variables differs from the initialized value of the state variables by more than a pre-determined threshold value.
US11/083,259 2005-03-17 2005-03-17 Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements Active 2028-09-24 US7621876B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US11/083,259 US7621876B2 (en) 2005-03-17 2005-03-17 Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements
US11/669,602 US8388544B2 (en) 2005-03-17 2007-01-31 System and method for measuring blood viscosity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US11/083,259 US7621876B2 (en) 2005-03-17 2005-03-17 Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/669,602 Continuation-In-Part US8388544B2 (en) 2005-03-17 2007-01-31 System and method for measuring blood viscosity

Publications (2)

Publication Number Publication Date
US20060211942A1 true US20060211942A1 (en) 2006-09-21
US7621876B2 US7621876B2 (en) 2009-11-24

Family

ID=37011307

Family Applications (1)

Application Number Title Priority Date Filing Date
US11/083,259 Active 2028-09-24 US7621876B2 (en) 2005-03-17 2005-03-17 Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements

Country Status (1)

Country Link
US (1) US7621876B2 (en)

Cited By (62)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070123779A1 (en) * 2005-03-17 2007-05-31 General Electric Company System and method for measuring blood viscosity
WO2007134394A1 (en) * 2006-05-22 2007-11-29 Uscom Limited Cardiac measurement system and method
US20080119741A1 (en) * 2006-11-22 2008-05-22 Friedman Bruce A Method and apparatus for automated vascular function testing
US20080119743A1 (en) * 2006-10-23 2008-05-22 Friedman Bruce A Method and apparatus for automated flow mediated dilation
US20080275351A1 (en) * 2007-05-02 2008-11-06 Siemens Corporate Research, Inc. Model-based pulse wave velocity measurement method
WO2008135892A2 (en) * 2007-05-02 2008-11-13 Philips Intellectual Property & Standards Gmbh Apparatus for performing pulse wave velocity measurements
US20080306383A1 (en) * 2007-06-11 2008-12-11 Miyuki Kono Blood flow measuring and evaluating apparatus
US20090093699A1 (en) * 2007-10-04 2009-04-09 Siemens Corporate Research, Inc. Method for monitoring myocardial wall thickness
US20090171213A1 (en) * 2006-03-01 2009-07-02 Koninklijke Philips Electronics, N.V. Linear array ultrasound transducer with variable patch boundaries
US7621876B2 (en) * 2005-03-17 2009-11-24 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements
US20100106016A1 (en) * 2008-10-23 2010-04-29 Skeletal Dynamics Llc Non-Invasive Blood Pressure Monitoring Device and Method
WO2010057495A2 (en) * 2008-11-18 2010-05-27 Sense A/S Methods, apparatus and sensor for measurement of cardiovascular quantities
US20100174202A1 (en) * 2009-01-08 2010-07-08 The General Electric Company Method and system for combining oscillometric blood pressure envelope data obtained from different signal processing paths
US7809420B2 (en) 2003-06-25 2010-10-05 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US7822453B2 (en) 2002-10-01 2010-10-26 Nellcor Puritan Bennett Llc Forehead sensor placement
US20110125034A1 (en) * 2009-11-26 2011-05-26 Nihon Kohden Corporation Vessel wall monitoring apparatus
US20110137182A1 (en) * 2009-12-04 2011-06-09 Anthony Bellezza Methods and devices for assessing intracranial pressure
US20110288420A1 (en) * 2010-05-19 2011-11-24 Seiko Epson Corporation Blood pressure measuring device and blood pressure measuring method
US20110295129A1 (en) * 2008-11-28 2011-12-01 Royal United Hosptial Bath NHS Trust Method of measuring blood pressure and apparatus for performing the same
WO2012009350A1 (en) * 2010-07-12 2012-01-19 Yale University Apparatus, systems and methods analyzing pressure and volume waveforms in the vasculature
US20120179053A1 (en) * 2009-07-31 2012-07-12 Fondazione Toscana Gabriele Monasterio Apparatus for measuring a propagation velocity of a blood pressure wave
US20120197129A1 (en) * 2011-01-31 2012-08-02 Dvx, Llc Flow Measurement Apparatus and Method
WO2012110042A1 (en) * 2011-02-17 2012-08-23 Sense A/S A method of and a system for determining a cardiovascular quantity of a mammal
US8257274B2 (en) 2008-09-25 2012-09-04 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US20120271166A1 (en) * 2009-10-15 2012-10-25 Wuxi Healthcare Medical Technology Co., Ltd. Method and device for detecting elasticity of viscous elastic medium
US8364220B2 (en) 2008-09-25 2013-01-29 Covidien Lp Medical sensor and technique for using the same
US20130079648A1 (en) * 2011-09-27 2013-03-28 Kabushiki Kaisha Toshiba Pulse wave velocity measuring device and pulse wave velocity measurement method
US8412297B2 (en) 2003-10-01 2013-04-02 Covidien Lp Forehead sensor placement
CN103156588A (en) * 2011-12-14 2013-06-19 精工爱普生株式会社 Blood pressure measurement apparatus and blood pressure measurement method
US20130184581A1 (en) * 2005-08-19 2013-07-18 Randall A. Hendriks Systems and methods for capture and display of blood pressure and ultrasound data
US8515515B2 (en) 2009-03-25 2013-08-20 Covidien Lp Medical sensor with compressible light barrier and technique for using the same
WO2014066859A1 (en) * 2012-10-26 2014-05-01 Graham Nichol Systems and methods for real-time assessment of the presence and quantity of carotid blood flow during cardiac arrest
US20140171811A1 (en) * 2012-12-13 2014-06-19 Industrial Technology Research Institute Physiology measuring system and method thereof
US8781548B2 (en) 2009-03-31 2014-07-15 Covidien Lp Medical sensor with flexible components and technique for using the same
US20140330087A1 (en) * 2013-05-01 2014-11-06 Medsensation, Inc. Devices and methods for obtaining physiological data
US20150327786A1 (en) * 2014-05-19 2015-11-19 Qualcomm Incorporated Method of Calibrating a Blood Pressure Measurement Device
WO2016019053A1 (en) * 2014-07-29 2016-02-04 Mayo Foundation For Medical Education And Research System and method for measurement of longitudinal and circumferential wave speeds in cylindrical vessels
US20160228698A1 (en) * 2013-09-19 2016-08-11 Koninklijke Philips N.V. Treatment device for the skin using radio-frequency electric current
US20160345840A1 (en) * 2006-02-21 2016-12-01 Uscom Ltd. Method and apparatus for producing a central pressure waveform in an oscillometric blood pressure system
US20160345930A1 (en) * 2015-05-25 2016-12-01 Seiko Epson Corporation Blood pressure measurement device and blood pressure measurement method
CN106491108A (en) * 2015-09-04 2017-03-15 塔塔咨询服务有限公司 The method and apparatus for monitoring for continuous blood pressure and estimating
CN106510763A (en) * 2016-09-19 2017-03-22 来安中衡物联网设备科技有限公司 Blood relationship spectrum imaging method and device
US9636447B2 (en) 2011-06-30 2017-05-02 Gambro Lundia Ab Filtering of a time-dependent pressure signal
WO2017083343A1 (en) * 2015-11-10 2017-05-18 General Electric Company System and method for estimating arterial pulse wave velocity
WO2017096314A1 (en) * 2015-12-02 2017-06-08 Echo Labs, Inc. Systems and methods for non-invasive blood pressure measurement
RU2648029C2 (en) * 2016-08-10 2018-03-21 Самсунг Электроникс Ко., Лтд. Device and method of blood pressure measurement
WO2018136135A1 (en) * 2017-01-18 2018-07-26 Physio-Control, Inc. Non-invasive blood pressure measurement using ultrasound
WO2018209140A1 (en) * 2017-05-10 2018-11-15 The Regents Of The University Of Michigan Automated ultrasound apparatus and methods to non-invasively monitor fluid responsiveness
US20180369065A1 (en) * 2017-06-23 2018-12-27 Stryker Corporation Patient monitoring and treatment systems and methods
CN110198664A (en) * 2016-11-10 2019-09-03 奥本大学 For assessing the method and system of blood vessel
CN111700638A (en) * 2019-03-18 2020-09-25 通用电气公司 Automated detection and localization of bleeding
US10786161B1 (en) * 2013-11-27 2020-09-29 Bodymatter, Inc. Method for collection of blood pressure measurement
US10959622B2 (en) * 2014-02-24 2021-03-30 Koninklijke Philips N.V. Method for determining pulse wave velocity in an artery
CN112998754A (en) * 2021-02-20 2021-06-22 苏州大学 Noninvasive blood flow detection method, detection system and detection equipment
US11083385B2 (en) * 2015-01-26 2021-08-10 University Of Ulsan Foundation For Industry Cooperation Apparatus for measuring blood circulation disorders, and method therefor
TWI753566B (en) * 2020-09-08 2022-01-21 國立陽明交通大學 Method for continuous beat-to-beat blood pressure monitoring
US11357415B2 (en) 2017-10-27 2022-06-14 Stryker Corporation Light-based non-invasive blood pressure systems and methods
US11391837B2 (en) * 2017-12-27 2022-07-19 Seiko Epson Corporation Ultrasonic device and ultrasonic measurement method
US11413005B2 (en) * 2017-08-14 2022-08-16 Stryker Corporation Constitutive equation for non-invasive blood pressure measurement systems and methods
US11452497B2 (en) * 2020-05-19 2022-09-27 Coravie Medical, Inc. Injectable hemodynamic monitoring devices, systems and methods
US11759168B2 (en) * 2017-11-14 2023-09-19 Koninklijke Philips N.V. Ultrasound vascular navigation devices and methods
EP4275611A1 (en) * 2022-05-09 2023-11-15 Koninklijke Philips N.V. Continuous blood pressure measurement system and method

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102551698B (en) * 2011-12-16 2015-03-18 广州医学院第二附属医院 Neck and brain arterial pulse wave speed measurement system
EP3021758A4 (en) * 2013-07-19 2017-03-15 Volcano Corporation Devices, systems, and methods for assessment of vessels
WO2016081517A2 (en) 2014-11-17 2016-05-26 Borkholder David A Pulse wave velocity, arterial compliance, and blood pressure
WO2018031714A1 (en) 2016-08-11 2018-02-15 Foundry Innovation & Research 1, Ltd. Systems and methods for patient fluid management
EP3725225A1 (en) 2015-02-12 2020-10-21 Foundry Innovation & Research 1, Ltd. Implantable devices for heart failure monitoring
EP3331426A1 (en) 2015-08-03 2018-06-13 Foundry Innovation&Research 1, Ltd. Devices and methods for measurement of vena cava dimensions, pressure, and oxygen saturation
WO2017109080A1 (en) * 2015-12-22 2017-06-29 Koninklijke Philips N.V. Multi-site continuous ultrasound flow measurement for hemodynamic management
JP6971017B2 (en) * 2016-04-12 2021-11-24 ソニーモバイルコミュニケーションズ株式会社 Detection device, detection method, and program
WO2018102435A1 (en) 2016-11-29 2018-06-07 Foundry Innovation & Research 1, Ltd. Wireless resonant circuit and variable inductance vascular implants for monitoring patient vasculature and fluid status and systems and methods employing same
US11206992B2 (en) 2016-08-11 2021-12-28 Foundry Innovation & Research 1, Ltd. Wireless resonant circuit and variable inductance vascular monitoring implants and anchoring structures therefore
US11701018B2 (en) 2016-08-11 2023-07-18 Foundry Innovation & Research 1, Ltd. Wireless resonant circuit and variable inductance vascular monitoring implants and anchoring structures therefore
EP3629921A1 (en) 2017-05-31 2020-04-08 Foundry Innovation & Research 1, Ltd. Implantable sensors for vascular monitoring

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5857975A (en) * 1996-10-11 1999-01-12 Dxtek, Inc. Method and apparatus for non-invasive, cuffless continuous blood pressure determination
US6117087A (en) * 1998-04-01 2000-09-12 Massachusetts Institute Of Technology Method and apparatus for noninvasive assessment of a subject's cardiovascular system
US6261233B1 (en) * 1996-01-05 2001-07-17 Sunlight Medical Ltd. Method and device for a blood velocity determination
US6264611B1 (en) * 1998-11-25 2001-07-24 Ball Semiconductor, Inc. Monitor for interventional procedures
US6647287B1 (en) * 2000-04-14 2003-11-11 Southwest Research Institute Dynamic cardiovascular monitor
US6648828B2 (en) * 2002-03-01 2003-11-18 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for measuring blood pressure using impedance plethysmography
US6676600B1 (en) * 1999-09-03 2004-01-13 Tensys Medical, Inc. Smart physiologic parameter sensor and method
US6730038B2 (en) * 2002-02-05 2004-05-04 Tensys Medical, Inc. Method and apparatus for non-invasively measuring hemodynamic parameters using parametrics
US20050015009A1 (en) * 2000-11-28 2005-01-20 Allez Physionix , Inc. Systems and methods for determining intracranial pressure non-invasively and acoustic transducer assemblies for use in such systems
US20050143640A1 (en) * 2003-12-30 2005-06-30 General Electric Company Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US20050154299A1 (en) * 2003-12-30 2005-07-14 Hoctor Ralph T. Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US20060079773A1 (en) * 2000-11-28 2006-04-13 Allez Physionix Limited Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions
US20060241459A1 (en) * 2005-02-08 2006-10-26 Tai Alan C Automatic signal-optimizing transducer assembly for blood flow measurement

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7621876B2 (en) * 2005-03-17 2009-11-24 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6261233B1 (en) * 1996-01-05 2001-07-17 Sunlight Medical Ltd. Method and device for a blood velocity determination
US5857975A (en) * 1996-10-11 1999-01-12 Dxtek, Inc. Method and apparatus for non-invasive, cuffless continuous blood pressure determination
US6117087A (en) * 1998-04-01 2000-09-12 Massachusetts Institute Of Technology Method and apparatus for noninvasive assessment of a subject's cardiovascular system
US6264611B1 (en) * 1998-11-25 2001-07-24 Ball Semiconductor, Inc. Monitor for interventional procedures
US6676600B1 (en) * 1999-09-03 2004-01-13 Tensys Medical, Inc. Smart physiologic parameter sensor and method
US6647287B1 (en) * 2000-04-14 2003-11-11 Southwest Research Institute Dynamic cardiovascular monitor
US20050015009A1 (en) * 2000-11-28 2005-01-20 Allez Physionix , Inc. Systems and methods for determining intracranial pressure non-invasively and acoustic transducer assemblies for use in such systems
US20060079773A1 (en) * 2000-11-28 2006-04-13 Allez Physionix Limited Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions
US6730038B2 (en) * 2002-02-05 2004-05-04 Tensys Medical, Inc. Method and apparatus for non-invasively measuring hemodynamic parameters using parametrics
US6648828B2 (en) * 2002-03-01 2003-11-18 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for measuring blood pressure using impedance plethysmography
US20050143640A1 (en) * 2003-12-30 2005-06-30 General Electric Company Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US20050154299A1 (en) * 2003-12-30 2005-07-14 Hoctor Ralph T. Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US7125383B2 (en) * 2003-12-30 2006-10-24 General Electric Company Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US7425199B2 (en) * 2003-12-30 2008-09-16 General Electric Company Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US20060241459A1 (en) * 2005-02-08 2006-10-26 Tai Alan C Automatic signal-optimizing transducer assembly for blood flow measurement

Cited By (109)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7822453B2 (en) 2002-10-01 2010-10-26 Nellcor Puritan Bennett Llc Forehead sensor placement
US7809420B2 (en) 2003-06-25 2010-10-05 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US7979102B2 (en) 2003-06-25 2011-07-12 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US7877126B2 (en) 2003-06-25 2011-01-25 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US7877127B2 (en) 2003-06-25 2011-01-25 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US7813779B2 (en) 2003-06-25 2010-10-12 Nellcor Puritan Bennett Llc Hat-based oximeter sensor
US8412297B2 (en) 2003-10-01 2013-04-02 Covidien Lp Forehead sensor placement
US20070123779A1 (en) * 2005-03-17 2007-05-31 General Electric Company System and method for measuring blood viscosity
US8388544B2 (en) 2005-03-17 2013-03-05 General Electric Company System and method for measuring blood viscosity
US7621876B2 (en) * 2005-03-17 2009-11-24 Ge Medical Systems Information Technologies, Inc. Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements
US9445787B2 (en) * 2005-08-19 2016-09-20 Fujifilm Sonosite, Inc. Systems and methods for capture and display of blood pressure and ultrasound data
US20130184581A1 (en) * 2005-08-19 2013-07-18 Randall A. Hendriks Systems and methods for capture and display of blood pressure and ultrasound data
US20160345840A1 (en) * 2006-02-21 2016-12-01 Uscom Ltd. Method and apparatus for producing a central pressure waveform in an oscillometric blood pressure system
US9622668B2 (en) * 2006-02-21 2017-04-18 Uscom Ltd. Method and apparatus for producing a central pressure waveform in an oscillometric blood pressure system
US8177718B2 (en) * 2006-03-01 2012-05-15 Koninklijke Philips Electronics N.V. Linear array ultrasound transducer with variable patch boundaries
US20090171213A1 (en) * 2006-03-01 2009-07-02 Koninklijke Philips Electronics, N.V. Linear array ultrasound transducer with variable patch boundaries
US20100312123A1 (en) * 2006-05-22 2010-12-09 Uscom Limited Cardiac measurement system and method
WO2007134394A1 (en) * 2006-05-22 2007-11-29 Uscom Limited Cardiac measurement system and method
US20080119743A1 (en) * 2006-10-23 2008-05-22 Friedman Bruce A Method and apparatus for automated flow mediated dilation
US8016761B2 (en) * 2006-10-23 2011-09-13 The General Electric Company Method and apparatus for automated flow mediated dilation
US8043223B2 (en) * 2006-11-22 2011-10-25 The General Electric Company Method and apparatus for automated vascular function testing
US20080119741A1 (en) * 2006-11-22 2008-05-22 Friedman Bruce A Method and apparatus for automated vascular function testing
WO2008135892A3 (en) * 2007-05-02 2009-04-30 Philips Intellectual Property Apparatus for performing pulse wave velocity measurements
WO2008135892A2 (en) * 2007-05-02 2008-11-13 Philips Intellectual Property & Standards Gmbh Apparatus for performing pulse wave velocity measurements
US20080275351A1 (en) * 2007-05-02 2008-11-06 Siemens Corporate Research, Inc. Model-based pulse wave velocity measurement method
JP2008302095A (en) * 2007-06-11 2008-12-18 Hitachi Ltd Apparatus for measuring and evaluating blood flow
US20080306383A1 (en) * 2007-06-11 2008-12-11 Miyuki Kono Blood flow measuring and evaluating apparatus
US8323197B2 (en) * 2007-06-11 2012-12-04 Hitachi, Ltd. Blood flow measuring and evaluating apparatus
US8271070B2 (en) * 2007-10-04 2012-09-18 Siemens Aktiengesellschaft Method for monitoring myocardial wall thickness
US20090093699A1 (en) * 2007-10-04 2009-04-09 Siemens Corporate Research, Inc. Method for monitoring myocardial wall thickness
US8364220B2 (en) 2008-09-25 2013-01-29 Covidien Lp Medical sensor and technique for using the same
US8257274B2 (en) 2008-09-25 2012-09-04 Nellcor Puritan Bennett Llc Medical sensor and technique for using the same
US20100106016A1 (en) * 2008-10-23 2010-04-29 Skeletal Dynamics Llc Non-Invasive Blood Pressure Monitoring Device and Method
WO2010048528A2 (en) * 2008-10-23 2010-04-29 Skeletal Dynamics Llc Non-invasive blood pressure monitoring device and method
WO2010048528A3 (en) * 2008-10-23 2010-08-12 Skeletal Dynamics Llc Non-invasive blood pressure monitoring device and method
US20110224529A1 (en) * 2008-11-18 2011-09-15 Sense A/S Methods, apparatus and sensor for measurement of cardiovascular quantities
WO2010057495A2 (en) * 2008-11-18 2010-05-27 Sense A/S Methods, apparatus and sensor for measurement of cardiovascular quantities
US9138161B2 (en) 2008-11-18 2015-09-22 Qualcomm Incorporated Methods, apparatus and sensor for measurement of cardiovascular quantities
WO2010057495A3 (en) * 2008-11-18 2010-08-19 Sense A/S Methods, apparatus and sensor for measurement of cardiovascular quantities
US20110295129A1 (en) * 2008-11-28 2011-12-01 Royal United Hosptial Bath NHS Trust Method of measuring blood pressure and apparatus for performing the same
US20100174202A1 (en) * 2009-01-08 2010-07-08 The General Electric Company Method and system for combining oscillometric blood pressure envelope data obtained from different signal processing paths
US8515515B2 (en) 2009-03-25 2013-08-20 Covidien Lp Medical sensor with compressible light barrier and technique for using the same
US8781548B2 (en) 2009-03-31 2014-07-15 Covidien Lp Medical sensor with flexible components and technique for using the same
US20120179053A1 (en) * 2009-07-31 2012-07-12 Fondazione Toscana Gabriele Monasterio Apparatus for measuring a propagation velocity of a blood pressure wave
US10426432B2 (en) 2009-10-15 2019-10-01 Wuxi Hisky Medial Technologies Co., Ltd. Method and device for detecting elasticity of viscous elastic medium
US10426433B2 (en) 2009-10-15 2019-10-01 Wuxi Hisky Medical Technologies Co., Ltd. Method and device for detecting elasticity of viscous elastic medium
US20120271166A1 (en) * 2009-10-15 2012-10-25 Wuxi Healthcare Medical Technology Co., Ltd. Method and device for detecting elasticity of viscous elastic medium
US9554771B2 (en) * 2009-10-15 2017-01-31 Beijing Sonicexpert Medical Technology Co., Ltd Method and device for detecting elasticity of viscous elastic medium
US10426434B2 (en) 2009-10-15 2019-10-01 Wuxi Hisky Medical Technologies Co., Ltd. Method and device for detecting elasticity of viscous elastic medium
US20110125034A1 (en) * 2009-11-26 2011-05-26 Nihon Kohden Corporation Vessel wall monitoring apparatus
US8814801B2 (en) * 2009-11-26 2014-08-26 Nihon Kohden Corporation Vessel wall monitoring apparatus
US20110137182A1 (en) * 2009-12-04 2011-06-09 Anthony Bellezza Methods and devices for assessing intracranial pressure
US9398861B2 (en) * 2009-12-04 2016-07-26 Third Eye Diagnostics, Inc. Methods and devices for assessing intracranial pressure
US20110288420A1 (en) * 2010-05-19 2011-11-24 Seiko Epson Corporation Blood pressure measuring device and blood pressure measuring method
WO2012009350A1 (en) * 2010-07-12 2012-01-19 Yale University Apparatus, systems and methods analyzing pressure and volume waveforms in the vasculature
US8852110B2 (en) * 2011-01-31 2014-10-07 Dvx, Llc Flow measurement apparatus and method
US20120197129A1 (en) * 2011-01-31 2012-08-02 Dvx, Llc Flow Measurement Apparatus and Method
CN103648374A (en) * 2011-02-17 2014-03-19 森思公司 A method of and a system for determining a cardiovascular quantity of a mammal
WO2012110042A1 (en) * 2011-02-17 2012-08-23 Sense A/S A method of and a system for determining a cardiovascular quantity of a mammal
CN106073731A (en) * 2011-02-17 2016-11-09 高通股份有限公司 For determining the method and system of mammal cardiovascular
US9603533B2 (en) 2011-02-17 2017-03-28 Qualcomm Incorporated Method of and a system for determining a cardiovascular quantity of a mammal
US9636447B2 (en) 2011-06-30 2017-05-02 Gambro Lundia Ab Filtering of a time-dependent pressure signal
US20130079648A1 (en) * 2011-09-27 2013-03-28 Kabushiki Kaisha Toshiba Pulse wave velocity measuring device and pulse wave velocity measurement method
US9445731B2 (en) * 2011-09-27 2016-09-20 Kabushiki Kaisha Toshiba Pulse wave velocity measuring device and pulse wave velocity measurement method
CN103156588A (en) * 2011-12-14 2013-06-19 精工爱普生株式会社 Blood pressure measurement apparatus and blood pressure measurement method
WO2014066859A1 (en) * 2012-10-26 2014-05-01 Graham Nichol Systems and methods for real-time assessment of the presence and quantity of carotid blood flow during cardiac arrest
US20140171811A1 (en) * 2012-12-13 2014-06-19 Industrial Technology Research Institute Physiology measuring system and method thereof
US20170020463A1 (en) * 2012-12-13 2017-01-26 Industrial Technology Research Institute Physiology measuring system and method thereof
US20140330087A1 (en) * 2013-05-01 2014-11-06 Medsensation, Inc. Devices and methods for obtaining physiological data
US20160228698A1 (en) * 2013-09-19 2016-08-11 Koninklijke Philips N.V. Treatment device for the skin using radio-frequency electric current
US10786161B1 (en) * 2013-11-27 2020-09-29 Bodymatter, Inc. Method for collection of blood pressure measurement
US11684270B2 (en) 2013-11-27 2023-06-27 Bodymatter, Inc. Method for collection of blood pressure measurement
US10959622B2 (en) * 2014-02-24 2021-03-30 Koninklijke Philips N.V. Method for determining pulse wave velocity in an artery
US10039455B2 (en) * 2014-05-19 2018-08-07 Qualcomm Incorporated Continuous calibration of a blood pressure measurement device
US10052036B2 (en) 2014-05-19 2018-08-21 Qualcomm Incorporated Non-interfering blood pressure measuring
US20150327786A1 (en) * 2014-05-19 2015-11-19 Qualcomm Incorporated Method of Calibrating a Blood Pressure Measurement Device
US20150327785A1 (en) * 2014-05-19 2015-11-19 Qualcomm Incorporated Continuous Calibration of A Blood Pressure Measurement Device
WO2016019053A1 (en) * 2014-07-29 2016-02-04 Mayo Foundation For Medical Education And Research System and method for measurement of longitudinal and circumferential wave speeds in cylindrical vessels
US11083385B2 (en) * 2015-01-26 2021-08-10 University Of Ulsan Foundation For Industry Cooperation Apparatus for measuring blood circulation disorders, and method therefor
US20160345930A1 (en) * 2015-05-25 2016-12-01 Seiko Epson Corporation Blood pressure measurement device and blood pressure measurement method
JP2016214736A (en) * 2015-05-25 2016-12-22 セイコーエプソン株式会社 Blood pressure measurement device and blood pressure measurement method
CN106491108A (en) * 2015-09-04 2017-03-15 塔塔咨询服务有限公司 The method and apparatus for monitoring for continuous blood pressure and estimating
WO2017083343A1 (en) * 2015-11-10 2017-05-18 General Electric Company System and method for estimating arterial pulse wave velocity
US10438355B2 (en) 2015-11-10 2019-10-08 General Electric Company System and method for estimating arterial pulse wave velocity
US11259753B2 (en) 2015-12-02 2022-03-01 Itamar Medical SPRY 2-21, Limited Partnership Systems and methods for detecting photoplethysmographic device usage
US9943267B2 (en) 2015-12-02 2018-04-17 Spry Health, Inc. Systems and methods for non-invasive respiratory rate measurement
WO2017096314A1 (en) * 2015-12-02 2017-06-08 Echo Labs, Inc. Systems and methods for non-invasive blood pressure measurement
RU2648029C2 (en) * 2016-08-10 2018-03-21 Самсунг Электроникс Ко., Лтд. Device and method of blood pressure measurement
CN106510763A (en) * 2016-09-19 2017-03-22 来安中衡物联网设备科技有限公司 Blood relationship spectrum imaging method and device
CN110198664A (en) * 2016-11-10 2019-09-03 奥本大学 For assessing the method and system of blood vessel
US11166637B2 (en) 2017-01-18 2021-11-09 Stryker Corporation Non-invasive blood pressure measurement using pulse wave velocity
US11877832B2 (en) 2017-01-18 2024-01-23 Stryker Corporation Non-invasive blood pressure measurement using pulse wave velocity
WO2018136135A1 (en) * 2017-01-18 2018-07-26 Physio-Control, Inc. Non-invasive blood pressure measurement using ultrasound
US11701092B2 (en) 2017-05-10 2023-07-18 Regents Of The University Of Michigan Automated ultrasound apparatus and methods to non-invasively monitor fluid responsiveness
WO2018209140A1 (en) * 2017-05-10 2018-11-15 The Regents Of The University Of Michigan Automated ultrasound apparatus and methods to non-invasively monitor fluid responsiveness
US11013488B2 (en) * 2017-06-23 2021-05-25 Stryker Corporation Patient monitoring and treatment systems and methods
US20180369065A1 (en) * 2017-06-23 2018-12-27 Stryker Corporation Patient monitoring and treatment systems and methods
US11413005B2 (en) * 2017-08-14 2022-08-16 Stryker Corporation Constitutive equation for non-invasive blood pressure measurement systems and methods
US11357415B2 (en) 2017-10-27 2022-06-14 Stryker Corporation Light-based non-invasive blood pressure systems and methods
US11759168B2 (en) * 2017-11-14 2023-09-19 Koninklijke Philips N.V. Ultrasound vascular navigation devices and methods
US11391837B2 (en) * 2017-12-27 2022-07-19 Seiko Epson Corporation Ultrasonic device and ultrasonic measurement method
CN111700638A (en) * 2019-03-18 2020-09-25 通用电气公司 Automated detection and localization of bleeding
US11452497B2 (en) * 2020-05-19 2022-09-27 Coravie Medical, Inc. Injectable hemodynamic monitoring devices, systems and methods
US20230000463A1 (en) * 2020-05-19 2023-01-05 Coravie Medical, Inc. Subcutaneous Hemodynamic Monitoring Devices, Systems and Methods
US11826195B2 (en) * 2020-05-19 2023-11-28 Coravie Medical, Inc. Subcutaneous hemodynamic monitoring devices, systems and methods
TWI753566B (en) * 2020-09-08 2022-01-21 國立陽明交通大學 Method for continuous beat-to-beat blood pressure monitoring
CN112998754A (en) * 2021-02-20 2021-06-22 苏州大学 Noninvasive blood flow detection method, detection system and detection equipment
EP4275611A1 (en) * 2022-05-09 2023-11-15 Koninklijke Philips N.V. Continuous blood pressure measurement system and method
WO2023217565A1 (en) 2022-05-09 2023-11-16 Koninklijke Philips N.V. Continuous blood pressure measurement system and method

Also Published As

Publication number Publication date
US7621876B2 (en) 2009-11-24

Similar Documents

Publication Publication Date Title
US7621876B2 (en) Continuous, non-invasive technique for determining blood pressure using a transmission line model and transcutaneous ultrasound measurements
US8388544B2 (en) System and method for measuring blood viscosity
US7125383B2 (en) Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US7425199B2 (en) Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring
US6632181B2 (en) Rapid non-invasive blood pressure measuring device
US6852083B2 (en) System and method of determining whether to recalibrate a blood pressure monitor
US7303530B2 (en) Transducer arrays with an integrated sensor and methods of use
Seo et al. Noninvasive arterial blood pressure waveform monitoring using two-element ultrasound system
WO2018052635A1 (en) Multi-model blood pressure estimation
US20060206032A1 (en) Method and apparatus for assessing hemodynamic parameters within the circulatory system of a living subject
US6662130B1 (en) Systems and methods for calibrating a distorted signal with another signal of known calibration
US20180360412A1 (en) System and Method for Non-Invasive Blood Pressure Measurement
US20160058409A1 (en) Ultrasonic blood pressure measurement apparatus and blood pressure measurement method
EP3784124A2 (en) Methods for the use of inherent frequency shifting mechanisms for sensors response reading with continuous wave excitation
US20180296182A1 (en) Apparatus And Method For Non-Invasive Determination Of Intracranial Pressure
Seo et al. Motion tolerant unfocused imaging of physiological waveforms for blood pressure waveform estimation using ultrasound
WO2016035041A1 (en) Apparatus and method for non-invasive pressure measurement of a fluid confined in a vessel with elastic or rigid walls fitted with an elastic window
EP4221578B1 (en) Estimating central blood pressure
JP2018033663A (en) Pulse wave measurement device and pulse wave measurement method
Seabra Development and Characterization of a Miniaturized Ultrasound Transceiver to Measure Blood Pressure

Legal Events

Date Code Title Description
AS Assignment

Owner name: GE MEDICAL SYSTEMS INFORMATION TECHNOLOGIES, INC.,

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:HOCTOR, RALPH T.;HERSH, LAWRENCE T.;FRIEDMAN, BRUCE A.;REEL/FRAME:015963/0742;SIGNING DATES FROM 20050303 TO 20050308

STCF Information on status: patent grant

Free format text: PATENTED CASE

CC Certificate of correction
FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 12