KR20150059654A - Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship - Google Patents
Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship Download PDFInfo
- Publication number
- KR20150059654A KR20150059654A KR1020130133301A KR20130133301A KR20150059654A KR 20150059654 A KR20150059654 A KR 20150059654A KR 1020130133301 A KR1020130133301 A KR 1020130133301A KR 20130133301 A KR20130133301 A KR 20130133301A KR 20150059654 A KR20150059654 A KR 20150059654A
- Authority
- KR
- South Korea
- Prior art keywords
- pressure
- flow rate
- coronary
- coronary artery
- flow
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Description
본 발명은 인체 혈류역학에 대한 전산유체역학적 해석 및 이를 이용한 진단기술에 관한 것이다. 보다 상세하게는, 인체 관상 동맥계의 병변에 따른 말초 미세혈관반의 생리학적 변화를 유량-압력 관계식으로 모사하고, 이에 기초하여 관상동맥 혈류의 전산유체해석을 수행하여 관상 동맥 유량과 관상동맥 내 압력을 구하는 방법에 관한 것이다.The present invention relates to computational fluid dynamics analysis of human hemodynamics and diagnostic techniques using the same. More specifically, the physiological changes of peripheral microvasculature according to the lesion of the coronary artery in the human body are simulated by a flow-pressure relationship, and computational fluid analysis of the coronary blood flow is performed based on the flow rate-coronary artery pressure The method comprising:
혈류 전산유체역학기술은 관상동맥 내 협착 병변을 진단하기 위하여 협착병변에 대한 혈관 성형수술을 시행하기 전에 시뮬레이션하는 방향으로 전개되고 있다. 그러나 이러한 기술은 아직 완전히 확립되지 않은 기술로, 보다 많은 임상연구와 함께 실용화 단계로 진입할 것으로 내다보고 있다.Computed tomography (CT) is being developed to simulate intracoronary stenosis before performing angioplasty for stenotic lesions. However, this technology is not yet fully established, and it is expected that it will enter into practical use with more clinical studies.
인체 심장의 각 부위와 연결되는 관상동맥혈관은, 도 1에서 도시된 바와 같이, 장기 전체에 넓게 분포한 모세혈관 및 소동맥(arteriole)을 포함하는 미세혈관반(microvascular bed)을 통해 해당 조직 세포로 혈액을 공급한다. 그런데, 이러한 미세혈관반은 혈관 내부 압력, 대사 산물(이산화탄소 등) 및 자율신경계의 상태에 따라 쉽게 수축 팽창하는 특징을 가지고 있어서 혈관의 직경 변화가 심하며, 이에 따라 혈류저항 및 혈류량도 급변한다. 따라서 만약 도 1과 같이 병변 상류(proximal) 동맥 내의 압력 Pa(aorta pressure)가 협착과 같은 병변에 의해 하류에서 압력 Pd(distal pressure)로 떨어질 경우, 병변 하류의 모든 혈류의 압력은 Pd 이하가 된다. 협착병변 등에 의한 Pd의 변화는 도 1과 같이 collateral vessel (CV)의 형성으로 인하여 억제될 수 있다. 미세혈관반의 압력 또한 기존의 값보다 떨어져서 미세혈관반의 미세혈관들이 수축하게 되면 미세혈관반의 혈류저항(Rvb)이 증가하여 혈류량 Q가 감소하게 된다. 결국, 동맥계 혈류 Q는 Pd 변화에 따라 결정되는 Rvb에 따라 결정되게 된다.As shown in Fig. 1, the coronary artery blood vessels connected to the respective parts of the human heart are connected to the corresponding tissue cells through a microvascular bed including a capillary vessel and an arteriole widely distributed throughout the organs. Blood is supplied. However, these microvascular vasculature easily swells and contracts according to the internal pressure of blood vessels, metabolites (carbon dioxide, etc.) and the state of autonomic nervous system, so that the diameter of blood vessels changes drastically, and blood flow resistance and blood flow also change rapidly. Therefore, if the pressure P a in the proximal artery drops to the pressure P d (distal pressure) from the lesion like stenosis as shown in FIG. 1, the pressure of all bloodstreams downstream of the lesion is P d Or less. Changes in P d due to stenosis or the like can be suppressed by formation of a collateral vessel (CV) as shown in Fig. When the microvessel half pressure is lower than the conventional value, the microvessel microvessels shrink, and the blood flow resistance (R vb ) of the microvascular half increases and the blood flow Q decreases. As a result, the arterial blood flow Q is determined according to R vb determined according to the change in P d .
한편, 현재 임상에서는 협착(stenosis) 병변의 진단을 위해 심근분획 혈류예비력(fractional flow reserve, 이하 FFR)을 측정한다. FFR은 병증의 정도를 파악하는 척도가 된다. FFR은 도 2와 같이 동맥계에 발생한 협착 동맥계(S)의 혈류량을 협착 병변이 없는 동일한 동맥계에 대한 가상의 정상 동맥계(N)의 혈류량으로 나눈 값이다. FFR 값은 정상인 경우 1이 되고 1보다 작을수록 병증의 정도가 심한 것을 나타낸다. 이러한 방법에 대한 구체적인 기술은 미국 특허 제7,775,988호(발명의 명칭: Method for Determining the Blood Flow in a Coronary Artery)에 공개되어 있다. 이 기술에 따르면 협착 동맥계와 동일한 형상을 갖는 가상의 정상 동맥계의 혈류량의 비는 수학식 1과 같이 정의되고, 수학식 2와 같은 압력의 비로 구할 수 있다. 한편 수학식 2의 압력의 비로서 혈류량의 비를 구하기 위하여는 다음과 같은 가정을 한다. 첫째, 심장 근육조직의 미세혈관반의 저항인 Rmyo는 혈관반에 adenosine을 투여할 경우 혈관의 최대 확장을 일으켜 충혈(hyperemia)상태를 만들어 최대 혈류를 만들며 동맥계 상류의 협착병변 등에 의한 압력변화와 상관없이 일정하고 가정한다. 둘째, 심장 근육조직의 정맥 압력 Pv는 항상 일정하고 가정한다.On the other hand, current clinical practice measures fractional flow reserve (FFR) for the diagnosis of stenosis lesions. FFR is a measure of the severity of the disease. FFR is the value obtained by dividing the blood flow of the stenosed arterial system (S) in the arterial system by the blood flow of the virtual steady arterial system (N) for the same arterial system without the stenotic lesion. The FFR value is 1 for normal, and less than 1 indicates severity of pathology. A detailed description of such a method is disclosed in U.S. Patent No. 7,775,988 entitled Method for Determining the Blood Flow in a Coronary Artery. According to this technique, the ratio of the blood flow of the virtual steady arterial system having the same shape as that of the stenosed arterial system is defined by Equation (1) and can be obtained by the ratio of the pressure as shown in Equation (2). In order to obtain the ratio of the blood flow as the ratio of the pressure in Equation (2), the following assumption is made. First, R myo , a resistance of the microvascular layer of cardiac muscle tissue, causes maximal expansion of blood vessels when adenosine is administered to the vascular bed to produce hyperemia, thereby creating a maximum blood flow. Pressure changes due to stenotic lesions upstream of the arterial system Regardless, it is assumed to be constant. Second, the vein pressure P v of the heart muscle tissue is always assumed to be constant.
상기와 같이 정의되는 FFR을 구하기 위해서는 환자에게 아데노신(adenosine)을 투여한 후, 동맥 내부로 미세 압력센서가 장착된 pressure guidewire를 넣어 직접 Pa와 Pd를 측정한다. 이와 같은 침습적인 방법으로 FFR을 구할 경우, 환자에게 적지 않은 고통과 경제적 비용 및 위험성을 초래하게 된다.To obtain the FFR defined above, adenosine is administered to the patient, and a pressure guidewire equipped with a micro-pressure sensor is inserted into the artery to directly measure P a and P d . If such an invasive method is used to obtain the FFR, the patient will suffer a considerable amount of pain and economic costs and risks.
침습적인 방법으로 FFR을 구할 경우 수반되는 환자의 고통, 경제적 부담, 위험성을 줄이기 위하여, 최근에는 CT 나 MRI 기반의 3차원 영상기술과 혈류시뮬레이션 기술이 활용하여 환자의 심혈관의 3차원 구조를 모델링하는 기술이 개발되고 있다. 미국 공개 특허 제2010/0241404호(발명의 명칭, Patient-specific Hemodynamics of the Cardio Vascular System)에는 이러한 전산 유체역학 해석을 활용한 진단기술(computational fluid dynamics diagnosis, CFDD)이 공개되어 있다. 이러한 전산 유체 역학을 활용한 진단 기술은 먼저 환자의 3차원 관상동맥 형상을 가상공간 상에 복원하여 환자의 가상 동맥계에 대한 3차원 모델을 얻고, 여기에 가상 혈류를 흘려보내는 컴퓨터 시뮬레이션을 수행한다. 모델에서 계산된 3차원 혈류 해석 데이터를 처리하여 임상과 동일한 형식으로 FFR을 구한다. In order to reduce the pain, economic burden, and risk associated with the patient's FFR if an invasive method is used, recently, three-dimensional imaging techniques based on CT or MRI and blood flow simulation technology are used to model the three- Technology is being developed. In U.S. Patent Application Publication No. 2010/0241404 (entitled " Patient-specific Hemodynamics of the Cardio Vascular System "), a computational fluid dynamics diagnosis (CFDD) using such computational fluid dynamics analysis is disclosed. This computational fluid dynamics diagnostic technique first reconstructs the patient's 3D coronary artery shape in virtual space, obtains a 3D model of the patient's virtual arterial system, and performs a computer simulation of flowing the virtual blood flow there . The computed three-dimensional blood flow analysis data is processed in the model to obtain the FFR in the same form as the clinical.
그러나 상기 미국 공개 특허 제2010/0241404호에 개시된 기술은 상기 미국 특허 제7,775,988호에서 제시된 가정을 근거로 하고 있으나, 최근에 이러한 가정이 잘못되었다는 것을 지적하는 연구결과가 있다. 예를 들면, Spaan, J.A.E. et al.의 연구에 의하면 도 1에서 Rvb는 충혈상태일지라도 동맥계 상류의 압력 변화(P=Pa-Pd)에 따라 변화되고 있음을 보이고 있다(Physiological Basis of Clinically Used Coronary Hemodynamic Indices. Circulation. 2006;113:446-455, 이하 참고 문헌1). 동맥계의 혈류역학에 대한 컴퓨터 해석 모델에 있어서 Pd, Q, Rvb에 대한 관계는 경계조건으로서 결정적인 역할을 한다. 그러나 이러한 경계조건을 얻기 위해 Rvb를 직접 측정하는 방법은 현재로서는 알려져 있지 않다. 또한, 종래의 경계조건을 구하기 위한 가정이 잘못된 것이라는 연구결과가 있다.However, the technique disclosed in the above-mentioned U.S. Patent Application Publication No. 2010/0241404 is based on the assumption presented in the above-mentioned U.S. Patent No. 7,775,988, but recent research has pointed out that this assumption is wrong. For example, according to a study by Spaan, JAE et al., It is shown that R vb in FIG. 1 is changed according to the pressure change (P = P a -P d ) upstream of the arterial system even if it is in a state of redemption (Physiological Basis of Clinically Used Coronary Hemodynamic Indices. Circulation. 2006; 113: 446-455, hereinafter referred to as reference 1). The relationship between P d , Q, and R vb in the computer analysis model of arterial hemodynamics plays a crucial role as a boundary condition. However, the method of directly measuring R vb to obtain these boundary conditions is not known at present. Further, there is a research result that the assumption for obtaining the conventional boundary condition is wrong.
본 발명은 동맥계 유체역학적 해석에 적용하기 위한 새로운 생리학적 경계조건을 구하는 방법을 제공하고, 이를 이용하여 보다 정확한 생리학적 압력-유량 관계에 기초한 관상동맥 유량과 관상동맥내 압력을 구하는 방법을 제공하는 것을 목적으로 한다.The present invention provides a method for obtaining new physiological boundary conditions for application to arterial hydrodynamic analysis and provides a method for obtaining a coronary flow rate and a pressure in a coronary artery based on a more accurate physiological pressure-flow relationship using the method .
본 발명은 인체 혈류 역학에 대한 전산유체역학적 해석을 통하여 FFR을 구하는 방법을 제공한다. 혈관 내의 유체의 흐름을 해석할 때, 해석하고자 하는 영역에 대한 지배 방정식과 경계조건 및 초기 조건이 필요하다. 혈관내에서의 혈류의 흐름을 지배하는 지배방적식으로 비압축성 나비에 스톡 방정식을 사용하여 해석할 수 있다. 또한, 혈관내에서의 혈류의 흐름은 뉴톤 유체의 흐름으로 가정할 수 있다. 또한, 나비에 스톡 방정식을 삼차원 유한 요소 법으로 해석할 수 있다. 미국 특허 공개 제2010/0241404호에는 유한 요소법으로 혈관내의 유체의 흐름을 해석하는 방법이 공개되어 있다.The present invention provides a method for obtaining FFR through computational fluid dynamics analysis on human hemodynamics. When analyzing the flow of fluid in a vessel, governing equations, boundary conditions and initial conditions are required for the region to be analyzed. It can be interpreted by using the incompressible butterfly equilibrium equation as the dominant equation that governs the flow of blood flow in the blood vessel. In addition, the flow of blood flow in the blood vessel can be assumed as the flow of Newtonian fluid. Also, the stock equations can be interpreted by the three-dimensional finite element method. U.S. Patent Publication No. 2010/0241404 discloses a method for analyzing the flow of a fluid in a blood vessel by a finite element method.
본 발명은, 혈관 내의 유체의 흐름을 해석할 때, 협착 병변이 있는 관상 동맥의 경계 조건으로 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist)의 관계식을 이용하는 것을 특징으로 한다. 본 발명에 따른 컴퓨터를 이용한 관상동맥 유량(MFR)과 관상동맥내의 압력(Pdist)을 구하는 방법은, 생리학적인 관상 동맥 유량(mass flow rate, MFR)과 관상동맥 내의 압력(Pdist) 을 근사하고, 근사된 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist) 을 경계조건으로 관상동맥 혈류역학 지배방정식의 만족하는 수치해를 구하는 단계와, 구해진 수치해가 사전에 정해식 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist)과의 관계식을 만족할 때까지 관상동맥 유량을 조절하는 단계를 포함한다. 또한, 상기 정해진 수학식은 MFR= 0.0186exp(Pdist/44464.308)-0.01982 을 만족하는 것을 특징으로 한다.The present invention is characterized by using a relational expression of a coronary flow rate (MFR) and a pressure (P dist ) in a coronary artery as a boundary condition of a coronary artery having a stenotic lesion when analyzing a flow of fluid in a blood vessel. The method of determining the coronary flow rate (MFR) and the pressure in the coronary artery (P dist ) using the computer according to the present invention approximates the physiological mass flow rate (MFR) and the pressure in the coronary artery (P dist ) (MFR) and the pressure in the coronary artery (P dist ), obtaining a satisfactory numerical solution of the coronary arterial hemodynamic control equation using the boundary condition as the boundary condition, and calculating the numerical solution satisfying the predetermined coronary flow rate (MFR) And adjusting the flow rate of the coronary artery until the relationship between the pressure in the coronary artery and the pressure in the coronary artery (P dist ) is satisfied. Also, the predetermined equation is characterized by satisfying MFR = 0.0186exp (Pdist / 44464.308) -0.01982.
본 발명은 환자의 심혈관 3차원 구조를 3차원 영상기술을 이용하여 가상공간상에 복원된 가상 관상동맥계에 대한 3차원 모델에 대하여 생리학적인 경계조건을 적용하여 가상 혈류에 대한 유체역학적 컴퓨터 시뮬레이션을 방법을 제공한다.The present invention relates to a three-dimensional model of a cardiovascular three-dimensional structure of a patient restored in a virtual space using a three-dimensional imaging technique, by applying a physiological boundary condition to a hydrodynamic computer simulation of a virtual blood flow ≪ / RTI >
본 발명에 따른 방법은 기존의 CFDD 기술에서 적용하고 있는 생리학적 경계조건의 오류 및 한계를 극복할 수 있는 타당한 생리학적 경계조건 수립의 방법을 제안하고 이를 통해 보다 정확한 예측진단 기술을 확보하는 데 기여한다.The method according to the present invention proposes a method of establishing a reasonable physiological boundary condition to overcome the errors and limitations of the physiological boundary conditions applied in the existing CFDD technology, do.
도 1은 병변이 있는 환자의 관상동맥의 모델을 나타내는 개략도
도 2는 FFR을 구하는 방법을 설명하기 위한 관상동맥의 개략도
도 3은 관상 동맥 출구의 a-k 커플된 관상 동맥 모델
도 4는 도 3에 도시된 모델에서 종래의 방법으로 협착 병변이 있는 관상동맥의 유량과 압력을 구하는 방법의 순서도
도 5는 관상 동맥 내의 압력과 유량의 관계를 나타내는 설명도
도 6은 본 발명에 따른 방법으로 협착 병변이 있는 관상동맥의 유량과 압력을 구하는 방법의 순서도
도 7은 본 발명에 따른 방법으로 구한 FFR과 종래의 방법으로 구한 FFR의 비교 그래프1 is a schematic diagram illustrating a model of a coronary artery in a lesioned patient;
2 is a schematic view of a coronary artery for explaining a method of obtaining FFR
Figure 3 shows the AK coupled coronary artery model of the coronary outlet
Figure 4 is a flow chart of a method for determining the flow rate and pressure of a coronary artery with a stenotic lesion according to a conventional method in the model shown in Figure 3
5 is an explanatory diagram showing the relationship between the pressure and the flow rate in the coronary artery
6 is a flowchart of a method for obtaining a flow rate and a pressure of a coronary artery having a stenotic lesion by the method according to the present invention
7 is a graph comparing the FFR obtained by the method according to the present invention and the FFR obtained by the conventional method
미국 공개특허 제2010/0241404호에는, 관상동맥 질병의 진단과 시술을 지원하기 위한 심혈관 시스템의 혈류역학을 예측하고 계량화하기 위한 전산 유체역학적인 방법이 개시되어 있다. 도 3은 상기 특허에서 제시된 coronary outlets a-k coupled to lumped parameter coronary vascular model이다. 상기 특허에서는 관상동맥 모델을 유한요소법으로 해석할 때, 정맥 압력 Pv가 일정한 값을 가지는 것으로 가정하였다. 또한, 미세혈관반의 저항인 Ra-micro, Rv-micro도 일정한 값으로 가정하였다. 일정한 값으로 가정된 미세혈관반 저항들을 바탕으로 관상동맥계 압력-유량관계를 얻어기 위하여는 도 4에 도시된 것과 같은 알고리즘을 이용하여 연산을 한다. 우선 초기 Ra-micro, Rv-micro 등의 말단의 parameter들을 일정한 값으로 정하고 CFD 과정을 수행한다. 그런 후, 관상동맥의 평균 유량 (Qc_mean)을 심박출량의 약 4.0%로 가정하여 Qc_mean이 심박출량의 4.0%가 될 때 까지 도 3의 혈관 컴플라이언스(vascular compliance) Ca와 Cim 등을 조절한다. 따라서 동맥계 말단 미세혈관반의 저항이 고정된 상태에서, 근사화된 Qc_mean 및 정상 동맥압 파형에 기초하여 동맥계의 C 값들을 조절함으로써 동맥계 유량을 결정하고 있다. U.S. Published Patent Application No. 2010/0241404 discloses a computational fluid dynamics method for predicting and quantifying hemodynamics of a cardiovascular system to assist in the diagnosis and treatment of coronary artery disease. FIG. 3 is a coronary outflow ac-coupled to lumped parameter coronary vascular model presented in the patent. When interpreting the coronary model, the finite element method in the patent, was assumed venous pressure P v is has a constant value. In addition, the resistance of the microvessel half, R a-micro and R v-micro , is assumed to be constant. In order to obtain the coronary arterial pressure-flow relationship based on the assumption of a certain value of microvascular half-resistances, the calculation is performed using the algorithm shown in FIG. First, initial parameters such as R a-micro and R v-micro are set to a constant value and CFD process is performed. Then, assuming that the mean flow rate (Q c_mean ) of the coronary artery is about 4.0% of the cardiac output, the vascular compliance C a and C im of FIG. 3 are measured until Q c_mean reaches 4.0% of the cardiac output . Therefore, the arterial flow rate is determined by adjusting the C values of the arterial system based on the approximated Q c_mean and steady - state arterial pressure waveform, with the resistance of the arterial end microvessel half fixed.
그러나 이러한 방법은 동맥계 말단 미세혈관반의 저항변화를 정밀하게 반영하는 방법이 아니며, 근사화된 Qc_mean과 압력 파형을 기준으로 동맥계의 C 값들을 조절하므로 정확한 FFR값 추산에는 여러 한계점을 지니고 있다.However, this method is not a method that accurately reflects the resistance change of the arterial end microvessel, and has several limitations in estimating the accurate FFR value because it controls the C values of the arterial system based on the approximated Q c_mean and the pressure waveform.
최근 연구 결과에 따르면, Rvb는 충혈상태에서도 압력차(ΔP=Pa-Pd)에 따라서 변동하는 것으로 밝혀졌다. 도 4는 참고문헌 1에서 제시된 인체의 관상 동맥 내의 압력(Pdist)과 유량(Q)과의 관계를 나타낸다. 도 4를 참조하면, 관상동맥 내에 발생된 협착병변 등에 의해 ΔP가 증가될 경우 Pdist가 정상적인 수준인 100 mmHg (= Pa,ΔP=0) 보다 떨어지게 된다. 이때 관상동맥 내의 유량은 생리학적 압력 수준 (Pdist > 약 40 mmHg) 에서는 선형적으로 감소하며 이보다 낮은 압력 수준 (Pdist < 약 40 mmHg)에서는 비선형적으로 감소하고 있음을 보여준다. 따라서 이때 최소 Rvb_min.의 값은 ΔP/Q = tanΘ 와 같음을 알 수 있다.Recent studies have shown that R vb varies with the pressure difference (ΔP = P a -P d ) even in the congested state. 4 shows the relationship between the pressure (P dist ) in the coronary artery of the human body and the flow rate (Q) shown in Reference 1. Referring to FIG. 4, when ΔP is increased due to a stenotic lesion or the like occurring in the coronary artery, P dist is lower than a normal level of 100 mmHg (= P a , ΔP = 0). At this time, the flow rate in the coronary arteries decreases linearly at the physiological pressure level (P dist > about 40 mmHg) and non-linearly decreases at the lower pressure level (P dist <about 40 mm Hg). Therefore, the minimum R vb_min. Is equal to? P / Q = tan?.
본 발명은 이러한 관상동맥계 상부 병변의 압력변화에 따른 말초혈관에서의 생리학적 유량변화를 동맥계 말단의 압력 및 유량의 경계조건으로 모델링하는 방법을 제공한다.The present invention provides a method of modeling the change in the physiological flow rate in the peripheral blood vessel according to the pressure change of the upper coronary artery lesion as the boundary condition of the pressure and the flow rate of the arterial end.
우선 도 4에 도시된 인체 동맥계의 생리학적 압력-유량관계를 도 5와 같이 curve fitting하여 아래의 수학식 3과 같은 유량(mass flow rate, MFR)을 압력으로 나타내는 압력-유량 관계식을 유도한다. First, the physiological pressure-flow relationship of the human arterial system shown in FIG. 4 is curve-fitted as shown in FIG. 5 to derive a pressure-flow relation representing a mass flow rate (MFR) as shown in Equation 3 below .
이러한 수학식 3은 도 5의 압력-유량의 관계에 대해 어느 정도의 오차를 가질 수 있다. 또한, 도 5와 같은 압력-유량의 관계는 환자마다 조금씩 다르기 때문에 이에 따른 오차가 어느 정도 존재한다. 그러나 수학식 3의 압력-유량 관계는 관상동맥계 내의 압력변동으로 인한 미세혈관반 혈류저항 Rvb의 변화를 압력에 따른 유량변동으로 대체할 수 있다. 따라서 기존의 고정된 Rvb를 갖는 경계조건을 사용하는 것보다는 좀 더 타당한 FFR 해석결과를 얻을 수 있게 된다. This equation (3) may have some error with respect to the pressure-flow relationship of FIG. In addition, since the relationship between the pressure and the flow rate as shown in Fig. 5 is slightly different for each patient, there is an error to some extent. However, the pressure-flow relationship in equation (3) can replace the change in the microvascular hemodynamic resistance R vb due to the pressure fluctuation in the coronary artery system with the fluctuation in the flow rate according to the pressure. Therefore, it is possible to obtain a more valid FFR analysis result than using the existing boundary condition with fixed R vb .
도 7은 본 발명에 따라서 동맥계 내의 압력변동에 따른 유량변화를 미세혈관반의 혈류저항인 Rvb의 변화에 따라 구하는 방법을 나타낸다. 먼저, 생리학적인 관상 동맥 유량(mass flow rate, MFR)과 관상동맥 내의 압력(Pdist) 을 근사하고, 근사된 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist) 을 경계조건으로 관상동맥 혈류역학 지배방정식의 만족하는 수치해를 구한다. 수치해를 구하기 위하여는 상기 미국 공개특허 제2010/0241404에 공개된 것과 같은 나비에 스토크 방정식을 지배방정식으로 이용하여 혈관내의 혈류를 해석한다. 다음으로 지배 방정식을 만족하는 수치해가 사전에 정해식 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist)과의 관계식을 만족하는지 판단한다. 만약 사전에 정해진 관계식을 만족하지 않으면, 반복적으로 구해진 관상동맥 유량 값을 약간씩 조절하여(섭동 방법, perturbation method) 지배방적식을 풀어서 정해진 오차 범위에 들어올 때까지 반복한다.FIG. 7 shows a method of obtaining a flow rate change according to a pressure fluctuation in an arterial system according to a change of R vb , a blood flow resistance of a microvasculature according to the present invention. First, the physiological approximation of the mass flow rate (MFR) and the pressure in the coronary artery (P dist ), the approximated coronary flow rate (MFR) and the pressure in the coronary artery (P dist ) Obtain the satisfactory numerical solution of the hemodynamic governing equations. To obtain the numerical solution, the Stokes equations are used as a governing equation to interpret the blood flow in the blood vessel as disclosed in the above-mentioned U.S. Patent Publication No. 2010/0241404. Next, it is determined whether the numerical solution satisfying the governing equations satisfies the relation between the predetermined coronary flow (MFR) and the pressure in the coronary artery (P dist ). If the pre-determined relationship is not satisfied, repeat the gauging of the coronary flow values repeatedly (perturbation method), loosen the governance equation and repeat until it reaches the specified error range.
도 8은 협착 병변으로 인해 정상보다 떨어진 Pdist의 경우, 수학식 1을 이용하여 FFR을 계산한 결과이다. 사각형의 점으로 도시된 값은 도 4에 도시된 방법인 constant Rmyo를 적용한 방법으로 구한 FFR 값을 나타낸다. 원형의 점으로 도시된 값은 도 7에 도시된 본 발명에 따른 방법인 variable Rmyo를 적용한 방법으로 구한 FFR 값을 나타낸다. 도면에 도시된 것과 같이, 본 발명에 따른 방법으로 구한 FFR값이 종래의 방법으로 구한 FFR 값보다 낮게 나타난는 것을 알 수 있다. 특히, 병변의 중증을 판단하는 기준이 되는 FFR = 0.8 값 에서 constant Rmyo 방법의 FFR 값은 variable Rmyo 방법의 FFR 값 보다 10 %이상의 오차를 보여주고 있어서 병변의 증상이 심각한 데도 그보다 덜 심각한 것으로 진단할 수 있는 가능성이 있다는 것을 보여 준다.FIG. 8 shows the result of calculating FFR using Equation 1 in the case of P dist , which is distant from the normal due to a stricture lesion. The values indicated by the points of the rectangles represent the FFR values obtained by applying the constant R myo method shown in FIG. The values shown as circular dots represent the FFR values obtained by applying the variable R myo method of the present invention shown in FIG. As shown in the figure, it can be seen that the FFR value obtained by the method according to the present invention is lower than the FFR value obtained by the conventional method. In particular, the FFR value of the constant R myo method at the FFR = 0.8 value, which is the criterion for the severity of the lesion, shows an error of more than 10% from the FFR value of the variable R myo method, It shows that there is a possibility to diagnose.
Claims (2)
생리학적인 관상 동맥 유량(mass flow rate, MFR)과 관상동맥 내의 압력(Pdist) 을 근사하고, 근사된 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist) 을 경계조건으로 관상동맥 혈류역학 지배방정식의 만족하는 수치해를 구하는 단계와,
구해진 수치해가 사전에 정해식 관상동맥 유량(MFR)과 관상동맥 내의 압력(Pdist)과의 관계식을 만족할 때까지 관상동맥 유량을 조절하는 단계를 포함하는 관상동맥 유량과 관상동맥 내의 압력을 구하는 방법.The method of obtaining the coronary flow rate (MFR) and the pressure (P dist ) in the coronary artery using the computer according to the present invention,
The physiologic approximation of the mass flow rate (MFR) and the pressure in the coronary artery (P dist ) and the approximated coronary flow rate (MFR) and the pressure in the coronary arteries (P dist ) Obtaining a satisfactory numerical solution of the governing equations,
A method for determining a coronary flow rate and a pressure in a coronary artery, the method comprising the step of adjusting the flow rate of the coronary artery until the obtained numerical solution satisfies the relationship between the predetermined coronary flow rate (MFR) and the pressure in the coronary artery (P dist ).
상기 정해진 수학식은,
인 관상동맥 유량과 관상동맥 내의 압력을 구하는 방법.The method according to claim 1,
The predetermined equation is:
To determine the coronary flow rate and the pressure in the coronary arteries.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020130133301A KR101753576B1 (en) | 2013-11-05 | 2013-11-05 | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020130133301A KR101753576B1 (en) | 2013-11-05 | 2013-11-05 | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020150102680A Division KR101623186B1 (en) | 2015-07-20 | 2015-07-20 | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| KR20150059654A true KR20150059654A (en) | 2015-06-02 |
| KR101753576B1 KR101753576B1 (en) | 2017-07-04 |
Family
ID=53490781
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020130133301A KR101753576B1 (en) | 2013-11-05 | 2013-11-05 | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship |
Country Status (1)
| Country | Link |
|---|---|
| KR (1) | KR101753576B1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20230029694A (en) | 2020-06-22 | 2023-03-03 | 코로벤티스 리서치 에이비 | Systems and methods for determining microvascular resistance reserve |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8315812B2 (en) | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
-
2013
- 2013-11-05 KR KR1020130133301A patent/KR101753576B1/en active IP Right Grant
Also Published As
| Publication number | Publication date |
|---|---|
| KR101753576B1 (en) | 2017-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106659400B (en) | Device for determining a fractional flow reserve value | |
| US10872698B2 (en) | Method and system for enhancing medical image-based blood flow computations using physiological measurements | |
| KR101677085B1 (en) | A framework for personalization of coronary flow computations during rest and hyperemia | |
| CN108511075B (en) | Method and system for non-invasively acquiring fractional flow reserve | |
| KR101611805B1 (en) | Method and system for patient-specific modeling of blood flow | |
| KR20160038038A (en) | Method and system for modeling blood flow with boundary conditions for optimized diagnostic performance | |
| WO2016030744A1 (en) | Method for determining patient-specific blood vessel information | |
| Cuomo et al. | An experimental–computational study of catheter induced alterations in pulse wave velocity in anesthetized mice | |
| CN113015497B (en) | Method and device for simulating blood flow of blood vessel inherent to object | |
| ES2933276T3 (en) | Patient-specific modeling of hemodynamic parameters in coronary arteries | |
| JP2022500133A (en) | Calculation of boundary conditions for virtual FFR and iFR calculations based on myocardial staining characteristics | |
| EP3457912A1 (en) | Devices and methods for stratification of patients for renal denervation based on intravascular pressure and wall thickness measurements | |
| Ojeda et al. | Sensitivity analysis and parameter estimation of a coronary circulation model for triple-vessel disease | |
| CN114052764A (en) | Method, apparatus, system and computer storage medium for obtaining fractional flow reserve | |
| KR101623186B1 (en) | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship | |
| US20220125324A1 (en) | System and method for determining a blood flow characteristic | |
| CN114947910A (en) | Coronary artery end microvascular resistance calculation method and FFR calculation method and system | |
| CN117438055A (en) | Noninvasive MRR calculation method based on computational fluid dynamics | |
| KR101753576B1 (en) | Method for obtaining mass flow rate and distal pressure of coronary vascular based on physiological pressure-flow relationship | |
| Zhang et al. | Branch flow distribution approach and its application in the calculation of fractional flow reserve in stenotic coronary artery | |
| US20210219924A1 (en) | Noninvasive Diagnostics of Proximal Heart Health Biomarkers | |
| CN110584696B (en) | Fractional flow reserve evaluation method and device and storage medium | |
| KR20180077140A (en) | Method for determining cardiovascular information using mathematical modeling of stenosed lesion | |
| Rabidou | Pulmonary Artery Hemodynamics Using MRI & CFD |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AMND | Amendment | ||
| AMND | Amendment | ||
| A107 | Divisional application of patent | ||
| J201 | Request for trial against refusal decision | ||
| J301 | Trial decision |
Free format text: TRIAL NUMBER: 2015101004113; TRIAL DECISION FOR APPEAL AGAINST DECISION TO DECLINE REFUSAL REQUESTED 20150720 Effective date: 20170220 |
|
| S901 | Examination by remand of revocation | ||
| GRNO | Decision to grant (after opposition) | ||
| GRNT | Written decision to grant |