CN113160921A - 一种基于血流动力学的数字人体心血管***的构建方法和应用 - Google Patents
一种基于血流动力学的数字人体心血管***的构建方法和应用 Download PDFInfo
- Publication number
- CN113160921A CN113160921A CN202110570819.1A CN202110570819A CN113160921A CN 113160921 A CN113160921 A CN 113160921A CN 202110570819 A CN202110570819 A CN 202110570819A CN 113160921 A CN113160921 A CN 113160921A
- Authority
- CN
- China
- Prior art keywords
- signals
- pulse
- human body
- digital human
- distributed
- 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.)
- Pending
Links
- 230000000004 hemodynamic effect Effects 0.000 title claims abstract description 77
- 210000000748 cardiovascular system Anatomy 0.000 title claims abstract description 30
- 238000010276 construction Methods 0.000 title claims abstract description 11
- 238000012544 monitoring process Methods 0.000 claims abstract description 29
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 230000002526 effect on cardiovascular system Effects 0.000 claims abstract description 17
- 230000008321 arterial blood flow Effects 0.000 claims abstract description 9
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 210000001367 artery Anatomy 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 33
- 230000005236 sound signal Effects 0.000 claims description 33
- 238000004458 analytical method Methods 0.000 claims description 20
- 201000010099 disease Diseases 0.000 claims description 19
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 19
- 210000002321 radial artery Anatomy 0.000 claims description 18
- 239000008280 blood Substances 0.000 claims description 15
- 210000004369 blood Anatomy 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 10
- 238000003745 diagnosis Methods 0.000 claims description 9
- 210000002302 brachial artery Anatomy 0.000 claims description 8
- 210000003414 extremity Anatomy 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 7
- 230000036996 cardiovascular health Effects 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 210000001168 carotid artery common Anatomy 0.000 claims description 5
- 238000013186 photoplethysmography Methods 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 5
- 210000002303 tibia Anatomy 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 230000026676 system process Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 11
- 239000003814 drug Substances 0.000 description 11
- 238000011160 research Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000000747 cardiac effect Effects 0.000 description 7
- 230000029058 respiratory gaseous exchange Effects 0.000 description 7
- 210000000709 aorta Anatomy 0.000 description 6
- 230000036772 blood pressure Effects 0.000 description 6
- 210000003141 lower extremity Anatomy 0.000 description 6
- 210000001364 upper extremity Anatomy 0.000 description 6
- 229940079593 drug Drugs 0.000 description 5
- 238000012827 research and development Methods 0.000 description 5
- 230000017531 blood circulation Effects 0.000 description 4
- 230000003205 diastolic effect Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 230000002792 vascular Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000009509 drug development Methods 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 210000004115 mitral valve Anatomy 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 230000036581 peripheral resistance Effects 0.000 description 3
- 230000035488 systolic blood pressure Effects 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 206010019280 Heart failures Diseases 0.000 description 2
- 230000004872 arterial blood pressure Effects 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000003190 augmentative effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000035487 diastolic blood pressure Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 210000002465 tibial artery Anatomy 0.000 description 2
- 210000000591 tricuspid valve Anatomy 0.000 description 2
- 210000000707 wrist Anatomy 0.000 description 2
- 241000905957 Channa melasoma Species 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 101000878457 Macrocallista nimbosa FMRFamide Proteins 0.000 description 1
- 206010031127 Orthostatic hypotension Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 210000001142 back Anatomy 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 210000000269 carotid artery external Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 208000002815 pulmonary hypertension Diseases 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
-
- 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
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- 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
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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
- A61B5/021—Measuring pressure in heart or blood vessels
-
- 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
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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
- A61B5/026—Measuring blood flow
- A61B5/029—Measuring or recording blood output from the heart, e.g. minute volume
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/20—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the management or administration of healthcare resources or facilities, e.g. managing hospital staff or surgery rooms
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
- G16H40/67—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for remote operation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H70/00—ICT specially adapted for the handling or processing of medical references
- G16H70/40—ICT specially adapted for the handling or processing of medical references relating to drugs, e.g. their side effects or intended usage
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Cardiology (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Physiology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- General Business, Economics & Management (AREA)
- Vascular Medicine (AREA)
- Hematology (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medicinal Chemistry (AREA)
- Pharmacology & Pharmacy (AREA)
- Toxicology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开了一种基于血流动力学的数字人体心血管***的构建方法和应用，根据分布式血流动力学监测***划分人体各部位，并确定数字人体模型***的颗粒度；基于动脉血管弹性管模型，构建动脉血流传输函数，并定义函数特征值的生理含义；将分布式血流动力学监测***采集到的信号、特征值、血流动力学参数输入至数字人体模型***中；通过有限元计算获得动态的分布式的数字人体心血管***。本发明通过人体各部位分布式特征值和计算出的血流动力学参数，构建出个体的基于血流动力学的动态数字人体模型。该模型能够实时动态地记录并显示人体各部位的心血管状态，具备丰富的医学价值。
Description
技术领域
本发明涉及数字人体***的构建方法，特别是涉及一种基于血流动力学的数字人体心血管***的构建方法和应用。
背景技术
血流动力学参数与心血管疾病有着密切的关系，其参数的提取较为盲目。无创血流动力学监测技术是包括脉搏波探测、心电信号探测、心音信号探测在内的血流动力学监测技术。各类信号表现出的波形特征(如形态、强度、速度、节律)能够反映出人体心血管***的许多生理和病理特征，如心律不齐、二尖瓣病变、主动脉瓣病变、高血压、肺动脉高压、心力衰竭等。无创血流动力学监测技术有着实时动态、无创、便捷等优势，适合包括临床诊疗、日常健康管理、运动管理在内的多场景心血管健康监测工作。
然而，现有绝大多数的产品、研究工作仅着眼于动脉***的某一节段(如颈动脉、肱动脉、桡动脉、下肢动脉)的血流动力学特征，而对整体动脉***的检测、研究较为缺乏。并且，人体不同部位的动脉的健康状态不尽相同，不同动脉的血流动力学状态也不尽相同。比如，人体的左、右上肢血压可能存在差异，人体上下肢血压也存在差异。人体心血管***不同部位、不同性质的病变，对血流动力学信号的波形特征的影响也不尽相同。现有无创血流动力学监测技术仅限于研究某些节段动脉的血流动力学特征，而对整体动脉***的检测、研究较为缺乏。
数字人体是当今医学科学技术、信息科学、生命科学、人工智能、***科学、计算科学和计算机技术的高度综合技术。现如今数字人体技术停留在基于医学影像的3D层析数字人体，在医学教育起着辅助作用。包括人体心血管***在内的生理参数尚未能构建数字化生理人体模型。
发明内容
发明目的：本发明的目的是提供一种基于血流动力学的数字人体心血管***的构建方法和应用，基于血流动力学的全身分布式监测***，通过脉搏信号、心电信号、心音信号的实时采集、特征信息提取、血流动力学参数计算，实现针对人体动脉***的监测。
技术方案：本发明所述的一种基于血流动力学的数字人体心血管***的构建方法，包括以下步骤：
(1)根据分布式血流动力学监测***划分人体各部位，并确定数字人体模型***的颗粒度；
(2)基于动脉血管弹性管模型，构建动脉血流传输函数，并定义函数特征值的生理含义；
(3)将分布式血流动力学监测***采集到的信号、特征值、血流动力学参数输入至数字人体模型***中；
(4)通过有限元计算获得动态的分布式的数字人体心血管***。
其中，所述分布式血流动力学监测***包括脉搏传感器、心电传感器、心音传感器和信号采集分析***；所述脉搏传感器包括设置于人体表浅动脉位置的动脉脉搏传感器和设置于人体微血管处的光电容积脉搏波传感器，所述心电传感器用于采集不同导联下的心电图信号，所述心音传感器用于采集不同部位下的心音信号；所述信号采集分析***与传感器相连并实时动态接收检测到的人体生理信号，并进行信号分析、特征值提取、血流动力学参数计算，以及信号、特征值、参数的存储。
其中，分布式血流动力学监测***对采集的信号处理，信号处理包括如下步骤：
(1)脉搏、心电、心音信号的预处理，用于提高信号信噪比并排除干扰信号；
(2)提取各类信号的特征值，各特征点之间的强度差与时间差；
(3)多信号协同分析，在各信号同步采集的基础上提取各信号之间的特征值；
(4)计算血流动力学参数，根据脉搏、心电、心音各信号及预设信息计算出多种心血管参数。
优选地，动脉血流传输函数为：
其中，P、Q分别为沿管截面平均的压力和流量；x为坐标，t为时间；L、R、C分别表征流动惯性、黏性阻力和管壁顺应性。
进一步地，动脉血管中，流动惯性
黏性阻力
管壁顺应性
其中，ρ、μ、r0、1、E、h分别是血液密度、血液黏滞系数、无张力状态血管半径、血管段长度、血管壁杨氏模量、管壁厚度。
其中，上述分布式血流动力学监测***监测的人体心血管参数包括：标准肢体导联心电信号、主动脉部位心音信号、左右桡动脉脉搏信号、左右肱动脉脉搏信号、左右脚背动脉脉搏信号、左右胫骨动脉脉搏信号、颈总动脉脉搏信号。
上述步骤(1)中，预处理包括傅里叶滤波、小波变换、自适应滤波和数学形态法。
上述步骤(3)中，多信号协同分析在脉搏、心电、心音信号同步采集的基础上，信号之间的特征值包含着脉搏波传播路线上心血管的健康信息；所述特征值包括：脉搏波到达时间PAT、脉搏波传播时间PTT、射血前期持续时间PEP和脉搏波传播速度PWV；其中，PAT＝PEP+PTT。
本发明还提供了上述方法构建得到的数字人体心血管***。
本发明还提供了上述基于基于血流动力学的数字人体心血管***在非疾病诊断方面的应用，如疾病机理研究、医学教学培训、药物研发、结合增强现实/虚拟现实AR/VR的远程医疗技术等。
(1)人体心血管健康状态预测功能为：基于用户的个体数字人体模型和血流动力学原理，计算估计出用户心血管状态的演变过程，包括已有疾病的演化过程和其他疾病的概率，并提出医疗建议。
(2)疾病机理研究为：通过大数据技术，统计疾病所对应的数字人体模型特征，进而计算估计出疾病的演化过程。该数字人体模型能够动态呈现人体各部位的负载(如血压)，对疾病机理的研究具有积极作用。
(3)医学教学培训，以分布式、动态的数字人体模型和基于大数据的丰富的样本，能够可视化地展示人体心血管***以及相关疾病的演化过程，对相关专业人员的理论技术教学培训起到积极作用。
(4)药物研发仿真功能，药企可将其药物作用数字化，加载至本专利提出的数字人体模型上，计算仿真出人体的心血管状态变化和疾病的演化过程。对药物的研发、给药时间位置的研究起到积极作用。
(5)远程医疗技术，本专利提出的数字人体模型，可结合5G技术、AR/VR技术，实现远程医疗。将用户的数字人体模型实时动态地显示在医生面前，提高诊疗质量。该数字人体模型对远程医疗的实现具有推动作用，对未来医疗模式、医疗资源的再分配起到积极作用。
进一步地，本发明所述脉搏传感器，数量2个及以上，分为动脉脉搏波传感器和光电容积脉搏波传感器两种类型。其中，动脉脉搏波传感器又可分为传感压力、应变的力学传感器或者基于超声原理的超声传感器。动脉脉搏传感器分布放置在人体表浅动脉位置处(如颈总动脉、颈外动脉、肱动脉、桡动脉、胫骨动脉、足背动脉)。光电容积脉搏波传感器则分布放置在人体微血管处(如手脚指尖、耳垂)。
所述心电传感器，导联方式数量1个及以上。其电极分布放置在四肢、胸前、胸后。可按需构建单极导联、双极导联、胸导联等各类导联方式，并采集不同导联下的心电图信号。
所述心音传感器，数量1个及以上。可按需分布放置在主动脉部位、肺部位、三尖瓣部位、二尖瓣部位位置处，采集不同部位下的心音信号。
所述信号采集分析***。包含脉搏/心电/心音各类信号采集、信号存储、信号处理、信号特征值提取、血流动力学参数计算、人机交互(显示、语音)、通讯(与手机、电脑等终端及云端)的功能。
基于一种基于血流动力学的全身分布式监测***，本发明提出一种相应的信号处理算法且构建出一种基于血流动力学的数字人体心血管模型***。
其中，信号处理算法包括对脉搏、心电、心音信号的预处理、各类信号的特征值提取算法、多信号协同分析算法、血流动力学参数计算，具体算法和模型如下：
(1)脉搏、心电、心音信号的预处理算法，包括傅里叶滤波、小波变换、自适应滤波、数学形态法，用于提高信号信噪比并排除干扰信号。
(2)特征值提取算法，提取的特征值包括：脉搏信号收缩期峰强度与时间、潮波强度与时间、舒张期峰强度与时间、波谷强度与时间、始射点强度与时间、收缩峰斜率、舒张期斜率、收缩期波形面积、舒张期波形面积、收缩期持续时间、舒张期持续时间、周期以及各特征点之间的强度差与时间差；心电信号P波、QRS波、T波、U波的强度与时间、周期以及各特征点之间的强度差与时间差；心音信号第一心音峰的强度与时间、第一心音始射点的时间、第二心音峰的强度与时间、第二心音始射点的时间、第一心音时间持续时间、第二心音时间持续时间、周期以及各特征点之间的强度差与时间差。并且，得益于本专利的传感器分布式放置的特性，脉搏信号各特征值具备着空间分布的性质。
(3)信号协同分析算法，旨在在各信号同步采集的基础上，提出各信号之间的特征值，包括：脉搏波到达时间、脉搏波传播时间、射血前期持续时间、脉搏波传播速度。这些特征值根据脉搏、心电、心音各信号特征点的不同有着不同的计算方式。并且，得益于本专利的传感器分布式放置的特性，脉搏波传播速度具备着空间分布的性质。
(4)血流动力学参数计算算法，旨在根据脉搏、心电、心音各信号及预设信息计算出多种心血管参数。预设信息按需包括：年龄、性别、身高、体重、体脂率、血液黏性。心血管参数包括：心率、血氧饱和度、收缩压、舒张压、波形特征量、平均动脉压、心搏出量、心输出量、心搏指数、心脏指数、心功指数、外周阻力、血管顺应性、血流半更新率、血流半更新时间、血流平均滞留时间、总血容量。此外，同时可解调出呼吸率。并且，得益于本专利的传感器分布式放置的特性，相关血流动力学参数具备着空间分布的性质。
本发明的技术难点如下：
(1)分布式血流动力学监测***需具备着高的检测密度和高通量信号采集能力，以实时动态地监测人体整个心血管***的状态，从而保证了数字人体心血管***的高空间颗粒度与时间分辨率。
(2)本发明提出的数字人体模型需具备丰富多样的信号分析技术，以实现精确且全面的生理性数字人体模型构建。
有益效果：
本发明提出的血流动力学监测***，继承了无创血流动力学监测技术时间上实时动态、使用感受上简单便携的特点，并且在空间上实现了全身分布式测量的优势，为数字人体心血管***模型的构建打下硬件基础。
根据本发明的一种基于血流动力学的全身分布式监测***，提出一种相应的信号处理算法且构建出一种基于血流动力学的数字人体心血管***模型。
本发明提出的数字人体模型能够可视化地呈现人体各个部位的心血管参数，并且动态显示人体各个部位动脉的血压、血流的状态和心脏的电-机械行为。
本发明构建的***通过人体各部位分布式特征值和计算出的血流动力学参数，构建出个体的基于血流动力学的动态数字人体模型。该模型能够实时动态地记录并显示人体各部位的心血管状态，具备丰富的医学价值。
附图说明
图1为基于血流动力学的全身分布式监测***的人体佩戴示意图。
图2为***采集到的脉搏信号和脉搏信号特征值示意图；图(a)为脉搏信号，图(b)为脉搏信号特征值示意图。
图3为***采集到的心电信号和心电信号特征值示意图；图(a)为***采集到的心电信号，图(b)心电信号特征值示意图。
图4为***采集到的原始心音信号、呼吸率信号、心音信号和心音信号特征值示意图；图(a)为***采集到的原始心音信号，图(b)为基于原始心音信号获得的呼吸率信号，图(c)为处理后的心音信号，图(d)为心音信号特征值示意图。
图5为心电信号、心音包络信号、肱动脉脉搏信号和桡动脉脉搏信号协同分析特征值示意图。
图6为血流动力学参数计算流程图。
图7为数字人体模型示意图。
图8为智能诊断人体疾病诊断原理示意图。
图9为基于数字人体模型的药物研发仿真流程示意图。
具体实施方式
下面结合实施例对本发明进一步地详细描述。
如图1所示，基于血流动力学的全身分布式监测***通常采用以下部件：信号采集分析***1、上肢脉搏传感器2、下肢脉搏传感器3、颈部脉搏传感器4、心电传感器5、心音传感器6。
所述信号采集分析***1实时动态接收脉搏传感器(上肢脉搏传感器2、下肢脉搏传感器3、颈部脉搏传感器4)、心电传感器5、心音传感器6检测到的人体生理信号，并进行信号分析、特征值提取、血流动力学参数计算，以及信号、特征值、参数的存储。此外，该***具备人机交互能力，能够通过显示图像、语音等方式显示人体生理信号、告知人体生理状态，同时使用者可对该***进行一定控制。该***同时具备与手机、电脑等终端通讯的作用，通过蓝牙、WIFI等无线通讯技术将数据实时显示、存储在各类终端上。并可通过互联网，将数据传输至云端服务器，建立大数据样本库，供人工智能技术学习。
所述脉搏传感器(上肢脉搏传感器2、下肢脉搏传感器3、颈部脉搏传感器4)可按需放置在各表浅动脉或微血管位置处，进行脉搏波信号的采集。在本实施例中，上肢脉搏传感器2选择放置在左右桡动脉、左右肱动脉位置处，下肢脉搏传感器3选择放置在左右脚背动脉、左右胫骨动脉位置处，颈部脉搏传感器放置在颈总动脉位置处。脉搏波信号通过信号采集分析***1采集并存储。脉搏波信号通过15Hz低通滤波滤除高频噪声，再通过数学形态法矫正基线漂移。图2(a)展示了桡动脉位置处处理后的脉搏波信号，其细节示意图如图2(b)所示。波形清晰显示出脉搏波的始射点O、收缩期峰A、潮波B、舒张期峰C和波谷D。通过数学形态法和动态阈值门限法提取发明内容中所述各特征值。凭借分布式脉搏传感器，四肢动脉和主动脉的各自的特征值可分别提取出来，血流动力学参数可分别计算出来。其特征值和参数的差异，给心血管疾病的诊治提供了丰富的分布式信息。
所述心电传感器5可按需构建肢体导联或胸导联方式。在本实施例中，心电电极分别放置在左右手腕和一侧脚腕处，构建标准肢体导联。心电信号通过信号采集分析***1采集并存储。心电信号通过100Hz低通滤波滤除高频噪声，再通过数学形态法矫正基线漂移。图3(a)展示了本***采集的肢体导联心电信号，其细节示意图如图3(b)所示。波形清晰显示出心电波形的P波、QRS波群、T波和U波。通过数学形态法和动态阈值门限法提取发明内容中所述各特征值。
所述心音传感器6可按需分布放置在主动脉部位、肺部位、三尖瓣部位、二尖瓣部位位置处，采集不同部位下的心音信号。在本实施例中，心音传感器放置在主动脉部位。心音信号通过信号采集分析***1采集并存储。图4(a)展示了本***采集的原始心音信号。心音信号的包络蕴含着呼吸带来的胸腔扩张的信息。如图4(b)显示，通过3Hz的低通滤波，提取出呼吸信号，用于呼吸率的计算。对原始心音信号进行20-100Hz带通滤波，实现心音信号的提取，如图4(c)显示。可清晰观测到第一心音事件与第二心音事件。通过归一化香农能量法、数学形态法提取出心音信号的包络(如图4fd))。进而通过数学形态法和动态阈值门限法提取发明内容中所述各特征值。
图5展示了多信号协同分析算法。在图5中，同步采集了心电信号、心音包络信号、肱动脉脉搏信号和桡动脉脉搏信号，对上述各类信号进行协同分析得到特征值，解调出脉搏波传播路线上心血管的健康信息。特征值包括：脉搏波到达时间PAT、脉搏波传播时间PTT、射血前期持续时间PEP、脉搏波传播速度PWV。脉搏波到达时间为心脏产生电信号到远端检测到脉搏波的时间差。脉搏波传播时间为心脏射血到远端检测到脉搏波的时间差。射血前期持续时间为心脏产生电信号到心脏射血的时间差。因此三者通常具有如下关系：PAT＝PEP+PTT。
脉搏波传播速度指脉搏波在动脉血管中传播的速度。得益于本专利的传感器分布式放置的特性，人体不同动脉的脉搏波传播速度可分别计算得到。在本实施例中，实现主动脉段、心脏-肱动脉段、桡动脉段、心脏-胫骨动脉段、胫骨-脚背动脉段的脉搏波传播速度计算。这些特征值根据脉搏、心电、心音各信号特征点的不同有着不同的计算方式。在本实施例中，通过心电信号的R波峰、心音信号S1心音峰、脉搏波信号收缩期峰计算出脉搏波到达时间PATp、脉搏波传播时间PTTp、射血前期持续时间PEP。相关特征值也可通过信号的始射点、中点等特征点计算得出。
血流动力学参数可通过上述脉搏2/3/4、心电5、心音6信号得出。图6展示了人体生理参数的计算流程图。呼吸率RR通过心音6信号的包络求解得出。心率HR和脉搏2/3/4、心电5、心音6信号周期T有关。包括收缩压Ps、舒张压Pd在内的血压和脉搏波到达时间PAT、脉搏波传播时间PTT、射血前期持续时间PEP、脉搏波传播速度PWV以及脉搏波波形有强相关性。在本实施例中，通过公式BP＝a*PWV+b计算血压。血氧饱和度通过光电容积脉搏波的交流直流分量计算得出。其他血流动力学参数计算公式列于下表：
| 波形特征值K | K＝(P<sub>m</sub>-P<sub>d</sub>)/(P<sub>s</sub>-P<sub>d</sub>) |
| 血管顺应性AC | AC＝0.283T/K<sup>2</sup> |
| 心搏出量SV | SV＝0.283T(P<sub>s</sub>-P<sub>d</sub>)/K<sup>2</sup> |
| 心输出量CO | CO＝17(P<sub>s</sub>-P<sub>d</sub>)/K<sup>2</sup> |
| 心功指数W | W＝T(P<sub>s</sub>-P<sub>d</sub>) |
| 外周阻力TPR | TPR＝P<sub>m</sub>/CO |
应理解的是，该公式仅为本发明的具体实施例，并不用于限制本发明。
在本实施方案中，通过本***采集一位26岁心血管健康的男性的血流动力学参数列于下表：
得益于血流动力学***全身分布式的特点，脉搏2/3/4信号、脉搏波传播速度PWV以及相关血流动力学参数(如血压、血氧饱和度、波形特征值K、血管顺应性等)具有空间分布式的差异。以此构建出人体动脉***的健康监测体系。如图7所示，在本实施例实现主动脉段、左右心脏-肱动脉段、左右桡动脉段、左右心脏-胫骨动脉段、左右胫骨-脚背动脉段的分布式监测。此外通过对各部位信号特征值的加权平均，实现了呼吸率RR、心率HR、心搏出量SV、心输出量CO等参量更准确的测量。
下表列出一位26岁心血管健康的男性在不同部位处的部分血流动力学相关参数：
| 部位 | 脉搏波到达/传播时间 | 收缩压 | 舒张压 |
| 主动脉段 | 215.5ms | 132mmHg | 75mmHg |
| 左心脏-肱动脉段 | 253.3ms | 120mmHg | 78mmHg |
| 右心脏-肱动脉段 | 256.7ms | 126mmHg | 78mmHg |
| 左桡动脉段 | 29.0ms | 121mmHg | 80mmHg |
| 右桡动脉段 | 28.2ms | 125mmHg | 79mmHg |
| 左心脏-胫骨动脉段 | 329.3ms | 150mmHg | 101mmHg |
| 右心脏-胫骨动脉段 | 318.7ms | 154mmHg | 99mmHg |
其中，主动脉段、左心脏-肱动脉段、右心脏-肱动脉段、左心脏-胫骨动脉段、右心脏-胫骨动脉段计算的是脉搏波到达时间，即心脏产生心电信号到远端感受到脉搏信号的时间差；左桡动脉段、右桡动脉段计算的是脉搏波传播时间，即脉搏波从近端传播至远端的时间差。
通过人体各部位分布式特征值和计算出的血流动力学参数，构建出个体的基于血流动力学的动态数字人体模型。该模型能够实时动态地记录显示人体各部位的心血管状态，具备丰富的医学价值。
本专利所述数字人体模型，其构建具体步骤如下：
步骤1，根据分布式血流动力学监测***划分人体各部位，决定数字人体模型颗粒度；
步骤2，基于动脉血管弹性管模型，构建动脉血流传输函数，并定义函数特征值的生理意义；
步骤3，将分布式血流动力学监测***采集到的信号、特征值、血流动力学参数输入至数字人体模型中；
步骤4，通过有限元计算，获得动态的分布式的心血管体系数字人体模型。
在本实施例中，分布式血流动力学监测***监测的人体心血管参数有：标准肢体导联心电信号、主动脉部位心音信号、左右桡动脉脉搏信号、左右肱动脉脉搏信号、左右脚背动脉脉搏信号、左右胫骨动脉脉搏信号、颈总动脉脉搏信号。因此，构建的数字人体模型可划分为主动脉段、心脏-肱动脉段、桡动脉段、心脏-胫骨动脉段、胫骨-脚背动脉段。
在本实施例中，动脉血流传输函数为：
其中，P、Q分别为沿管截面平均的压力和流量；x为坐标，t为时间；L、R、C分别表征流动惯性、黏性阻力和管壁顺应性。
在动脉血管中，流动惯性
黏性阻力
管壁顺应性
其中ρ、μ、r0、l、E、h分别是血液密度、血液黏滞系数、无张力状态血管半径、血管段长度、血管壁杨氏模量、管壁厚度。
预设信号(如年龄、性别、身高、体重)能够估计出各血管段长度。在本实施例中，针对一位26岁心血管健康的男性，身高178cm。其主动脉段血管长度估计为20cm，心脏-肱动脉段血管长度估计为45cm，桡动脉段血管长度估计为23cm，心脏-胫骨动脉段血管长度估计为105cm，胫骨-脚背动脉段血管长度估计为33cm。
本***采集的I级信息(脉搏2/3/4、心电5、心音波形信号)反应了各动脉段的压力和流量，以及压力和流量在血管中传播的时间。本实施例中各动脉段的压力和传播时间已列于前表。
本***计算的II级信息(血流动力学参数)，能够计算出各动脉段的血管顺应性和外周阻力。计算方式上述可见。
综合上述3级信息，可通过动脉血流传输函数计算出主动脉段、心脏-肱动脉段、桡动脉段、心脏-胫骨动脉段、胫骨-脚背动脉段流动惯性、黏性阻力和管壁顺应性，进而计算出血液密度、血液黏滞系数、血管壁杨氏模量、管壁厚度等生理参数。
最后通过动脉血流传输函数，细化人体心血管***的划分，降低数字人体模型的颗粒度，获得空间分布率更高的人体心血管***模型。
其功能与用途有：智能诊断人体疾病的精确医疗技术；基于血流动力学预测人体心血管健康状态；疾病机理研究；医学教学培训；药物研发仿真；结合增强现实/虚拟现实AR/VR的远程医疗技术。
在本实施例中，具体描述智能诊断人体疾病的精确医疗技术和药物研发仿真技术。
如图8所示，经用户预设年龄、性别、身高、体重等基本信息并使用本***采集I级信息(脉搏2/3/4、心电5、心音波形信号)、计算的II级信息(血流动力学参数)后，***能够通过疾病与信号的相关性矩阵和神经网络技术综合诊断相关疾病。得益于分布式血流动力学监测***，相应I级信息和II级信息包含人体分布信息，对相关疾病的诊断极大提高了准确度。例如：血栓的诊断和病灶位置的判断、姿态性低血压的诊断、心衰的***性病因溯源等。
如图9所示，企业在药物研发、医疗器械研发的过程中，可先构建相应疾病的数字人体模型。在数字人体模型动态演变的过程中，将作用数字化的药物或医疗器械作为外界激励作用在数字人体模型上，数字人体模型进而对该激励做出反馈。企业可设定人体模型的相关特征参数以评价该药物或医疗器械的效用。此外，企业可罗列药物或医疗器械的相关参数，并定义数字人体模型的评价函数，以快速获得其产品的最优性质。数字人体模型作为一种安全的、高效的、快速的、易使用的研发仿真工具，这对药物或医疗器械剂量和时机的研究起到积极促进作用。
Claims (10)
1.一种基于血流动力学的数字人体心血管***的构建方法，其特征在于，包括以下步骤：
(1)根据分布式血流动力学监测***划分人体各部位，并确定数字人体模型***的颗粒度；
(2)基于动脉血管弹性管模型，构建动脉血流传输函数，并定义函数特征值的生理含义；
(3)将分布式血流动力学监测***采集到的信号、特征值、血流动力学参数输入至数字人体模型***中；
(4)通过有限元计算获得动态的分布式的数字人体心血管***。
2.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，所述分布式血流动力学监测***包括脉搏传感器、心电传感器、心音传感器和信号采集分析***；所述脉搏传感器包括设置于人体表浅动脉位置的动脉脉搏传感器和设置于人体微血管处的光电容积脉搏波传感器，所述心电传感器用于采集不同导联下的心电图信号，所述心音传感器用于采集不同部位下的心音信号；所述信号采集分析***与传感器相连并实时动态接收检测到的人体生理信号，并进行信号分析、特征值提取、血流动力学参数计算，以及信号、特征值、参数的存储。
3.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，分布式血流动力学监测***对采集的信号处理，信号处理包括如下步骤：
(1)脉搏、心电、心音信号的预处理，用于提高信号信噪比并排除干扰信号；
(2)提取各类信号的特征值，各特征点之间的强度差与时间差；
(3)多信号协同分析，在各信号同步采集的基础上提取各信号之间的特征值；
(4)计算血流动力学参数，根据脉搏、心电、心音各信号及预设信息计算出多种心血管参数。
4.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，动脉血流传输函数为：
其中，P、Q分别为沿管截面平均的压力和流量；x为坐标，t为时间；L、R、C分别表征流动惯性、黏性阻力和管壁顺应性。
5.根据权利要求4所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，动脉血管中，流动惯性
黏性阻力
管壁顺应性
其中，ρ、μ、r0、l、E、h分别是血液密度、血液黏滞系数、无张力状态血管半径、血管段长度、血管壁杨氏模量、管壁厚度。
6.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，分布式血流动力学监测***监测的人体心血管参数包括：标准肢体导联心电信号、主动脉部位心音信号、左右桡动脉脉搏信号、左右肱动脉脉搏信号、左右脚背动脉脉搏信号、左右胫骨动脉脉搏信号、颈总动脉脉搏信号。
7.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，步骤(1)中，预处理包括傅里叶滤波、小波变换、自适应滤波和数学形态法。
8.根据权利要求1所述的基于血流动力学的数字人体心血管***的构建方法，其特征在于，步骤(3)中，多信号协同分析在脉搏、心电、心音信号同步采集的基础上，信号之间的特征值包含着脉搏波传播路线上心血管的健康信息；所述特征值包括：脉搏波到达时间PAT、脉搏波传播时间PTT、射血前期持续时间PEP和脉搏波传播速度PWV；其中，PAT＝PEP+PTT。
9.一种权利要求1～8任一项所述方法构建得到的数字人体心血管***。
10.根据权利要求9所述的数字人体心血管***在非疾病诊断方面的应用。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110570819.1A CN113160921A (zh) | 2021-05-26 | 2021-05-26 | 一种基于血流动力学的数字人体心血管***的构建方法和应用 |
| US17/795,881 US20230218178A1 (en) | 2021-05-26 | 2021-07-06 | Construction method and application of digital human cardiovascular system based on hemodynamics |
| PCT/CN2021/106795 WO2022246987A1 (zh) | 2021-05-26 | 2021-07-16 | 一种基于血流动力学的数字人体心血管***的构建方法和应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110570819.1A CN113160921A (zh) | 2021-05-26 | 2021-05-26 | 一种基于血流动力学的数字人体心血管***的构建方法和应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113160921A true CN113160921A (zh) | 2021-07-23 |
Family
ID=76877418
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110570819.1A Pending CN113160921A (zh) | 2021-05-26 | 2021-05-26 | 一种基于血流动力学的数字人体心血管***的构建方法和应用 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20230218178A1 (zh) |
| CN (1) | CN113160921A (zh) |
| WO (1) | WO2022246987A1 (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113854976A (zh) * | 2021-10-26 | 2021-12-31 | 南京大学 | 一种脉搏信号中带潮汐波搏动波形的自动筛选方法 |
| CN114366046A (zh) * | 2022-01-19 | 2022-04-19 | 福州九候生医科技有限公司 | 脉搏波数据处理方法及相关设备 |
| TWI775711B (zh) * | 2022-02-22 | 2022-08-21 | 昌泰科醫股份有限公司 | 基於血流動力學及弦脈分析且與肝火旺盛/心火旺盛相關的特定生理綜合症偵測方法及系統 |
| CN117010299A (zh) * | 2023-08-09 | 2023-11-07 | 山东大学 | 基于血流动力学耦合模型的脑组织血流情况预测*** |
| WO2024082231A1 (zh) * | 2022-10-20 | 2024-04-25 | 深圳迈瑞生物医疗电子股份有限公司 | 监测心血管参数的方法和医疗设备 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104138253B (zh) * | 2013-05-11 | 2016-06-15 | 吴健康 | 一种无创动脉血压连续测量方法和设备 |
| CN104188642A (zh) * | 2014-09-24 | 2014-12-10 | 中国科学院合肥物质科学研究院 | 一种无创获取中心动脉压力波形及其相关参数的装置与获取方法 |
| CN105726000B (zh) * | 2016-01-29 | 2019-03-22 | 北京工业大学 | 一种基于四肢血压脉搏的心脏血管功能参数的装置 |
| CN107928655B (zh) * | 2017-12-14 | 2023-07-11 | 江苏理工学院 | 一种基于人体动脉3d血管网络模型的中心动脉血压测定方法 |
| US20220237864A1 (en) * | 2019-06-26 | 2022-07-28 | The Regents Of The University Of California | Reduced order model for computing blood flow dynamics |
-
2021
- 2021-05-26 CN CN202110570819.1A patent/CN113160921A/zh active Pending
- 2021-07-06 US US17/795,881 patent/US20230218178A1/en active Pending
- 2021-07-16 WO PCT/CN2021/106795 patent/WO2022246987A1/zh active Application Filing
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113854976A (zh) * | 2021-10-26 | 2021-12-31 | 南京大学 | 一种脉搏信号中带潮汐波搏动波形的自动筛选方法 |
| CN114366046A (zh) * | 2022-01-19 | 2022-04-19 | 福州九候生医科技有限公司 | 脉搏波数据处理方法及相关设备 |
| CN114366046B (zh) * | 2022-01-19 | 2024-03-22 | 福州九候生医科技有限公司 | 脉搏波数据处理方法及相关设备 |
| TWI775711B (zh) * | 2022-02-22 | 2022-08-21 | 昌泰科醫股份有限公司 | 基於血流動力學及弦脈分析且與肝火旺盛/心火旺盛相關的特定生理綜合症偵測方法及系統 |
| WO2024082231A1 (zh) * | 2022-10-20 | 2024-04-25 | 深圳迈瑞生物医疗电子股份有限公司 | 监测心血管参数的方法和医疗设备 |
| CN117010299A (zh) * | 2023-08-09 | 2023-11-07 | 山东大学 | 基于血流动力学耦合模型的脑组织血流情况预测*** |
| CN117010299B (zh) * | 2023-08-09 | 2024-02-09 | 山东大学 | 基于血流动力学耦合模型的脑组织血流情况预测*** |
Also Published As
| Publication number | Publication date |
|---|---|
| US20230218178A1 (en) | 2023-07-13 |
| WO2022246987A1 (zh) | 2022-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Miao et al. | Multi-sensor fusion approach for cuff-less blood pressure measurement | |
| CN104873186B (zh) | 一种可穿戴的动脉检测装置及其数据处理方法 | |
| CN104138253B (zh) | 一种无创动脉血压连续测量方法和设备 | |
| CN113160921A (zh) | 一种基于血流动力学的数字人体心血管***的构建方法和应用 | |
| JP5984088B2 (ja) | 非侵襲的連続血圧モニタリング方法及び装置 | |
| CN103892818B (zh) | 一种无创中心动脉血压测量方法和设备 | |
| US10092268B2 (en) | Method and apparatus to monitor physiologic and biometric parameters using a non-invasive set of transducers | |
| CN108366749A (zh) | 动态血压与生命体征监测装置、***和方法 | |
| WO2012149652A1 (en) | Method and apparatus for estimating myocardial contractility using precordial vibration signals | |
| CN102488503A (zh) | 连续血压测量装置 | |
| CN110236508A (zh) | 一种无创血压连续监测方法 | |
| CN105595983A (zh) | 一种血压测量装置及提高血压测量准确度的方法 | |
| CN110881967A (zh) | 一种无创多节段外周动脉血管弹性功能检测方法及其仪器 | |
| JP2008295517A (ja) | 漢方医における脈診の分析システムと方法 | |
| Vakulenkoa et al. | New application of blood pressure monitor with software environment Oranta-AO based on arterial oscillography methods | |
| CN113040738B (zh) | 血压检测装置 | |
| Zhang et al. | Key feature selection and model analysis for blood pressure estimation from electrocardiogram, ballistocardiogram and photoplethysmogram | |
| EP4368105A1 (en) | Method for generating synchronous electrocardiograms on basis of two lead asynchronous electrocardiograms | |
| Shokouhmand et al. | Fingertip Strain Plethysmography: Representation of Pulse Information based on Vascular Vibration | |
| US20220249055A1 (en) | Non-invasive, real-time, beat-to-beat, ambulatory blood pressure monitoring | |
| CN209733969U (zh) | 一种脉搏波传播时间的测量设备 | |
| CN113226170B (zh) | 一种心脏舒张功能评估方法、设备和*** | |
| CN203341726U (zh) | 一种新型智能血压计 | |
| Zhe et al. | Design and research of continuous blood pressure monitoring system based on lora | |
| Koay et al. | Cuffless Blood Pressure Meter with Mobile Application for Home-Care Service |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |