CN113686735B - 测量血液凝固属性的方法及装置 - Google Patents
测量血液凝固属性的方法及装置 Download PDFInfo
- Publication number
- CN113686735B CN113686735B CN202111088760.9A CN202111088760A CN113686735B CN 113686735 B CN113686735 B CN 113686735B CN 202111088760 A CN202111088760 A CN 202111088760A CN 113686735 B CN113686735 B CN 113686735B
- Authority
- CN
- China
- Prior art keywords
- blood sample
- blood
- sample
- coagulation
- measured
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 80
- 230000023555 blood coagulation Effects 0.000 title claims abstract description 67
- 210000004369 blood Anatomy 0.000 claims abstract description 185
- 239000008280 blood Substances 0.000 claims abstract description 185
- 230000015271 coagulation Effects 0.000 claims abstract description 52
- 238000005345 coagulation Methods 0.000 claims abstract description 52
- 230000008569 process Effects 0.000 claims abstract description 45
- 230000003287 optical effect Effects 0.000 claims abstract description 41
- 238000001514 detection method Methods 0.000 claims abstract description 35
- 238000009792 diffusion process Methods 0.000 claims abstract description 21
- 239000002245 particle Substances 0.000 claims abstract description 20
- 230000001678 irradiating effect Effects 0.000 claims abstract description 10
- 238000012014 optical coherence tomography Methods 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 20
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 14
- 230000035602 clotting Effects 0.000 claims description 11
- 206010053567 Coagulopathies Diseases 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 5
- 230000003197 catalytic effect Effects 0.000 claims description 3
- 239000003153 chemical reaction reagent Substances 0.000 claims description 2
- 238000007711 solidification Methods 0.000 claims description 2
- 230000008023 solidification Effects 0.000 claims description 2
- 210000002381 plasma Anatomy 0.000 abstract description 13
- 238000012544 monitoring process Methods 0.000 abstract description 12
- 238000005259 measurement Methods 0.000 abstract description 10
- 239000000523 sample Substances 0.000 description 163
- 239000007788 liquid Substances 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 9
- 239000007787 solid Substances 0.000 description 7
- 208000007536 Thrombosis Diseases 0.000 description 6
- 230000006835 compression Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 6
- 210000001772 blood platelet Anatomy 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000001727 in vivo Methods 0.000 description 5
- 102000009123 Fibrin Human genes 0.000 description 4
- 108010073385 Fibrin Proteins 0.000 description 4
- BWGVNKXGVNDBDI-UHFFFAOYSA-N Fibrin monomer Chemical compound CNC(=O)CNC(=O)CN BWGVNKXGVNDBDI-UHFFFAOYSA-N 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 4
- 239000011324 bead Substances 0.000 description 4
- 238000002091 elastography Methods 0.000 description 4
- 229950003499 fibrin Drugs 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 238000013169 thromboelastometry Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000006698 induction Effects 0.000 description 3
- 238000004879 turbidimetry Methods 0.000 description 3
- 239000012190 activator Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 239000005995 Aluminium silicate Substances 0.000 description 1
- 230000005653 Brownian motion process Effects 0.000 description 1
- UXVMQQNJUSDDNG-UHFFFAOYSA-L Calcium chloride Chemical compound [Cl-].[Cl-].[Ca+2] UXVMQQNJUSDDNG-UHFFFAOYSA-L 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 235000012211 aluminium silicate Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004159 blood analysis Methods 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 238000005537 brownian motion Methods 0.000 description 1
- 239000001110 calcium chloride Substances 0.000 description 1
- 229910001628 calcium chloride Inorganic materials 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 208000024200 hematopoietic and lymphoid system neoplasm Diseases 0.000 description 1
- 230000002439 hemostatic effect Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000011503 in vivo imaging Methods 0.000 description 1
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
-
- G01N15/01—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N2011/006—Determining flow properties indirectly by measuring other parameters of the system
- G01N2011/008—Determining flow properties indirectly by measuring other parameters of the system optical properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N2013/003—Diffusion; diffusivity between liquids
Abstract
本发明公开了一种测量血液凝固属性的方法及装置,该方法包括:利用光源照射待测血样,通过干涉检测系统采集待测血样随时间变化的光信号;计算光信号在不同延迟时间下的相关系数;计算相关系数随延迟时间的衰减系数,通过衰减系数得到待测血样的扩散系数,利用待测血样的扩散系数与血样粘度之间的关系分析待测血样的凝固过程。该方法通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测,可以用于分析全血或血浆的凝血属性,测量过程中,不需要外力在血样中诱发机械波,只需要光束照射到血样上,即可分析监控血液凝固过程。
Description
技术领域
本发明涉及血液分析技术领域,特别涉及一种测量血液凝固属性的方法及装置。
背景技术
血液凝固是指血液从流动的液体状态变成胶冻状凝块的过程。在凝血过程中,血浆中的纤维蛋白源转变为不溶的血纤维,血纤维交织成网,将很多血细胞网罗在内,形成血凝块。凝血功能的检测可用于诊断凝血功能障碍、血液肿瘤等疾病,评估心脑血管疾病风险,监测药物引起的凝血功能改变。因此,血液凝固属性的测量对于疾病诊断、治疗具有重要意义。
目前,测量血液凝固的方法包括光学比浊法、磁珠检测法、血栓弹力图法、旋转式血栓弹力图法、凝血与血小板功能分析法、光学相干弹性成像方法等。
光学比浊法中,通过测量血样在凝固过程中,光强的吸收或散射的变化,监控血液的凝固过程。向样品中加入凝血激活剂后,随着样品中纤维蛋白凝块的形成过程,透过样品的光强或被样品散射的光强逐步改变,当样品完全凝固以后,透射光或散射光的强度不再变化。这种方法的缺点是样品的光学差异、测试杯的光洁度、样品气泡等都会成为测量的干扰因素。
磁珠检测法中,在样品杯四周设置四个线圈(两个驱动线圈和两个感应线圈),样品杯两侧的驱动线圈产生恒定的交变磁场,导致样品杯中血样里添加的去磁小钢珠受到磁化作用,钢珠产生振荡运动,从而引起样品杯两侧的感应线圈内磁通量变化,产生感应电流,被系统记录。在加入凝血激活剂后,随着纤维蛋白凝块的形成,血浆的粘稠度增加,小钢珠振幅逐渐减小,感应线圈内的电流减小。通过测量感应电流,分析钢珠的振动幅度,监测血液的凝固过程。
由于光学比浊法和磁珠检测法中,被测对象为血浆样品,因此需要对全血进行离心处理,而从血浆对血液凝固问题进行分析时,不能反映血小板在止血过程中的重要作用,因此不能真实完整地评估全血的凝血属性。
在血栓弹力图法(Thromboelastography,TEG)检测时,将待测全血加入到圆柱形样品杯中,将一个悬挂的探针浸泡在杯中的血样里,圆柱形样品杯以一定的角度往复旋转。当血液呈液体状态时,样品杯的往复转动不能带动探针的转动。随着样品凝固过程的进展,因纤维蛋白和血小板产生黏附作用,样品杯的旋转带动悬挂探针转动。探针的运动通过传感器转换为电信号,经数据采集后进行分析和处理。通过对悬挂探针转动开始时间的记录和转动幅度变化分析,评估血液凝固属性。
旋转式血栓弹力图法(Rotational Thromboelastometry,ROTEM)与血栓弹力图法的原理类似,但结构稍有不同。旋转式血栓弹力计测量过程中,悬挂的探针在外力驱动下产生一定角度的往复运动,盛装全血的样品杯固定不动。在血样未发生凝固,呈液态时,悬挂探针能够自由转动。随着血液的凝固,样品杯和悬挂探针之间产生黏附作用,固定的样品杯阻碍悬挂探针的转动,探针的转动幅度与黏附强度呈反比,利用位移传感器测量探针转动幅度,即可生成粘弹力在血凝过程中的变化曲线,分析血液凝固属性。
凝血与血小板功能分析法(Sonoclot,SCT)也是通过检测血凝过程中粘弹性的变化,分析凝血属性。将一个中空的探针插入样品杯中的血样里,中空探针与压电传感器连接。随着压电传感器的振动,中空探针随着产生垂直方向的直线往复振动。当血样处于液体状态时,探针的自由振动不受限制。随着血液凝固的进行,凝块对探针的运动产生阻力。阻力的变化可被探针上的传感器检测,从而反映血液凝块形成过程中,血液粘弹性特征的变化,评估血液凝固属性。
血栓弹力图法、旋转式血栓弹力图法和凝血与血小板功能分析法实现了全血的凝血特性分析,能提供完整的血液凝集发生、发展过程信息。但是,通过探头的运动参数测量,描述血液粘弹性时,它们之间没有准确的关联,所受干扰因素较多,临床检测时缺乏统一的参考标准。
光学相干弹性成像中,利用外力在全血样品中诱发剪切波,通过光学相干层析成像系统检测剪切波,利用剪切波和杨氏模量的定量关系,可以计算血样的弹性属性。液态的血液中无法检测到剪切波,随着血液凝固,血液中出现较慢速度剪切波,而固态的血液中剪切波速度较快。因此,通过剪切波速度测量,可以表征血液凝固过程中弹性的变化,评估血液的凝固属性。
另一种光学相干弹性成像中,利用外力在全血样品中诱发压缩波,通过光学相干层析成像系统检测样品中的压缩波,计算压缩波传播过程中振幅衰减系数,利用衰减系数和样品粘性的关系,分析血样的凝固属性。液态的血液中,样品粘性较小,压缩波衰减系数较小。随着血液凝固,样品粘性增大,压缩波衰减系数增大。而固态的血液中,粘性达到最大值,压缩波衰减系数达到最大。因此,通过压缩波的衰减系数测量,可以表征血液凝固过程中粘度的改变,从而分析血液的凝固属性。
但在光学相干弹性成像方法中,需要一个外部装置在血样中诱发机械波,增加了系统的复杂性。此方法也难以实现多个样品的同时检测,检测通量较低。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种测量血液凝固属性的方法,该方法通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测。
本发明的另一个目的在于提出一种测量血液凝固属性的装置。
为达到上述目的,本发明一方面实施例提出了一种测量血液凝固属性的方法,包括以下步骤:利用光源照射待测血样,通过干涉检测系统采集所述待测血样随时间变化的光信号;计算所述光信号在不同延迟时间下的相关系数;计算所述相关系数随延迟时间的衰减系数,通过所述衰减系数得到所述待测血样的扩散系数,利用所述待测血样的扩散系数与血样粘度之间的关系分析所述待测血样的凝固过程。
本发明实施例的测量血液凝固属性的方法,通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测,可以用于分析全血或血浆的凝血属性,例如开始凝固时间点、凝固过程时间长度和凝固血块粘度等指标。测量过程中,不需要外力在血样中诱发机械波,只需要光束照射到血样上,即可分析监控血液凝固过程,可以实现在体的无接触血凝监控,无需将血样取出放入样品杯中,可以在凝血部位直接检测血凝过程,在不同样品之间移动成像光束时,可以实现多个血液样品的高通量检测。
另外,根据本发明上述实施例的测量血液凝固属性的方法还可以具有以下附加的技术特征:
在本发明的一个实施例中,通过所述干涉检测系统对所述待测血样进行照射和采集所述光信号,所述干涉检测系统包括光学相干层析成像系统(Optical CoherenceTomography,OCT)、马赫-曾德干涉仪(Mach-Zehnder Type Interferometer)、共路干涉仪(Common-Path Interferometer)和迈克尔逊干涉仪(Michelson-Type Interferometer),所述光学相干层析成像系统包括基于扫频光源的光学相干层析成像系统(Swept SourceOCT)和光谱域的光学相干层析成像系统(Spectral Domain OCT)。
在本发明的一个实施例中,所述计算所述光信号在不同延迟时间下的相关系数,包括:计算所述待测血样的同一空间位置点在不同时间段的相关系数,和/或计算所述待测血样的同一个空间范围内在不同时间点的相关系数。
在本发明的一个实施例中,所述计算所述光信号在不同延迟时间下的相关系数为:
其中,g(τ)为相关系数,I(t,z)为t时刻采集的信号序列,I(t+τ,z)为t+τ时刻采集的信号序列,τ为延迟时间,z为所述待测血样的血样深度方向,Cov[I(t,z)I(t+τ,z)]表示I(t,z)序列和I(t+τ,z)序列的协方差,Var[I(t,z)]和Var[I(t+τ,z)]分别表示I(t,z)序列和I(t+τ,z)序列的方差。
在本发明的一个实施例中,所述衰减系数γ与所述待测血样的扩散系数D之间的关系为:
其中,n为折射率,λ为光波长,θ为散射角;
所述待测血样的扩散系数D与血样粘度V之间的关系为:
其中,K表示玻尔兹曼常数,T表示样品温度,R表示散射粒子半径。
在本发明的一个实施例中,所述分析所述待测血样的凝固过程包括:分析所述待测血样的全血凝固开始时间、凝固速率、凝块粘度。
在本发明的一个实施例中,所述待测试样包括:血浆样品、稀释的全血样品、稀释的血浆样品。
在本发明的一个实施例中,还包括:在所述待测血样中加入催化试剂,改变所述待测血样的凝固速度。
在本发明的一个实施例中,还包括:通过扫描光源对所述待测血样的多个位置或多个所述待测血样进行照射,采集所述待测血样的多个位置或多个所述待测血样随时间变化的光信号。
为达到上述目的,本发明另一方面实施例提出了一种测量血液凝固属性的装置,包括:采集模块,用于利用光源照射待测血样,通过干涉检测系统采集所述待测血样随时间变化的光信号;计算模块,用于计算所述光信号在不同延迟时间下的相关系数;分析模块,用于计算所述相关系数随延迟时间的衰减系数,通过所述衰减系数得到所述待测血样的扩散系数,利用所述待测血样的扩散系数与血样粘度之间的关系分析所述待测血样的凝固过程。
本发明实施例的测量血液凝固属性的装置,通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测,可以用于分析全血或血浆的凝血属性,例如开始凝固时间点、凝固过程时间长度和凝固血块粘度等指标。测量过程中,不需要外力在血样中诱发机械波,只需要光束照射到血样上,即可分析监控血液凝固过程,可以实现在体的无接触血凝监控,无需将血样取出放入样品杯中,可以在凝血部位直接检测血凝过程,在不同样品之间移动成像光束时,可以实现多个血液样品的高通量检测。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的测量血液凝固属性的方法流程图;
图2为根据本发明一个实施例的基于扫频光源的光学相干层析成像系统结构示意图;
图3为根据本发明一个实施例的光谱域的光学相干层析成像系统结构示意图;
图4为根据本发明一个实施例的M模式光学相干层析成像图;
图5为根据本发明一个实施例的相关系数随延迟时间的变化曲线示意图;
图6为根据本发明一个实施例的血液凝固过程中的衰减系数变化示意图;
图7为根据本发明一个实施例的血液凝固属性分析流程示意图;
图8为根据本发明一个实施例的测量血液凝固属性的装置结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的测量血液凝固属性的方法及装置。
首先将参照附图描述根据本发明实施例提出的测量血液凝固属性的方法。
图1为根据本发明一个实施例的测量血液凝固属性的方法流程图。
如图1所示,该测量血液凝固属性的方法包括以下步骤:
在步骤S101中,利用光源照射待测血样,通过干涉检测系统采集待测血样随时间变化的光信号。
可以理解的是,在血液凝固过程中,形态从液态到固态的变化时,血液中细胞等粒子随机运动(如布朗运动)幅度存在差异,使用光学成像系统,通过监测血液样品中散射粒子的随机运动,实现血液的凝固属性检测。血液处于液体状态时,散射粒子的随机运动剧烈,成像系统中接收的光信号之间相关性低;血液于固体状态时,散射粒子的随机运动受限,成像系统中接收的光信号之间相关性高,本发明的实施例可以利用该差异对血液凝固过程进行分析。
具体地,通过光源对待测血样进行照射,并通过干涉检测系统采集照射后的待测血样随时间变化的光信号,通过采集的光信号进行后续的血液凝固属性的分析,检测血液中散射粒子的随机运动时,不需要接触血液样品,只需要光束照射到血液样品上,可以实现在体成像,也易于实现多样本同时检测。
在本发明的实施例中,对待测血样进行照射和采集光信号时可以通过干涉检测系统实现,干涉检测系统包括光学相干层析成像系统、马赫-曾德干涉仪、共路干涉仪和迈克尔逊干涉仪,所述光学相干层析成像系统包括基于扫频光源的光学相干层析成像系统和光谱域的光学相干层析成像系统。
作为一种具体的实施方式,本发明实施例可以采用光学相干层析成像系统,可以对血样在一定深度范围内成像,采集X-Y平面内某点位置随时间变化的信号,此信号称为M模式信号。光学相干层析成像系统结构可以是扫频光学相干层析成像系统,如图2所示,也可以是光谱域的光学相干层析成像系统,如图3所示。其中,图2为扫频光学相干层析成像系统。系统主要包括扫频光源,耦合器、参考光路、样品光路和光电探测器等。耦合器实现分光作用,并可以实现参考光路反射的光束与样品光路散射的光束进行干涉。干涉光由光电探测器检测。扫频光源在每一时刻输出单波长(或窄带)激光,光电探测器可以检测不同波长干涉光的强度。图3为光谱域的光学相干层析成像系统。系统主要包括宽谱光源,耦合器、参考光路、样品光路、光栅和相机等。耦合器实现分光作用,并可以实现参考光路反射的光束与样品光路散射的光束进行干涉。宽谱光源在同一时刻发射多波长的光束,不同波长的干涉光由光栅分光后,相机不同位置像素检测不同波长的干涉光强度。
图4中显示的M模式光学相干图像,是在同一X-Y平面点进行采样获得的,Z方向为血样深度方向,横轴为时间,采样时X-Y位置固定不变。因此,图像每一列表示,在同一时间t采集到的随深度z变化的信号序列I(t,z)。利用此M模式图像,即可计算信号的相关性。
可以理解的是,本发明采用的干涉检测系统,包括光学相干层析成像和几类干涉仪,可以获得待测血样中Z方向(光束方向)上不同深度的血液散射属性分布信息。当光束照射到X-Y平面内的某单点时,经过干涉信号分析,可以获得不同深度的血液散射属性,进而通过Z方向空间变化序列或时间变化序列计算相关性,分析血液扩散系数和粘度。改变光束在X-Y平面的照射位置,利用干涉检测系统可以实现三维成像。在测量光信号方面,相比其他的采集系统,通过干涉检测系统采集光信号可以达到更高的分辨率,并具有深度Z方向的成像功能。
在步骤S102中,计算光信号在不同延迟时间下的相关系数。
在本发明的一个实施例中,计算光信号在不同延迟时间下的相关系数,包括:计算待测血样的同一空间位置点在不同时间段的相关系数,和/或计算待测血样的同一个空间范围内在不同时间点的相关系数。
理想情况下,如果血样中颗粒的随机运动被束缚(如固体状态),t时刻采集的信号序列I(t,zi)(i=1,…,n)与t+τ时刻采集的信号序列I(t+τ,zi)(i=1,…,n)数值不变,信号之间的互相关性为1;如果血样中颗粒的随机运动较大(如液体状态),不同时间采集的I(t,zi)和I(t+τ,zi)信号序列值差异较大,信号之间的互相关性小于1。I(t,zi)表示在t时刻深度位置z1至zn这一段空间范围内的信号序列。
计算相关性时,也可以采用另外一种方式。即计算同一深度(或深度范围)z在不同时间点采集的信号I(tj,z)(j=1,…,m)与I(tj+τ,z)(j=1,…,m)信号序列之间的互相关,其中I(tj,z)表示深度z位置在t1时刻至tm时刻这一段时间范围内的信号序列,I(tj+τ,z)表示深度z位置在t1+τ时刻至tm+τ时刻这一段时间范围内的信号序列,τ表示两个序列的时间延迟。理想情况下,如果血样中颗粒的随机运动被束缚(如固体状态),I(tj,z)和I(tj+τ,z)信号序列值不变,信号之间的互相关性为1。如果血样中颗粒的随机运动较大(如液体状态),I(tj,z)和I(tj+τ,z)信号序列值差异较大,信号之间的互相关性小于1。
在本发明的一个实施例中,计算光信号在不同延迟时间下的相关系数为:
其中,g(τ)为相关系数,I(t,z)为t时刻采集的信号序列,I(t+τ,z)为t+τ时刻采集的信号序列,τ为延迟时间,z为待测血样的血样深度方向,Cov[I(t,z)I(t+τ,z)]表示I(t,z)序列和I(t+τ,z)序列的协方差,Var[I(t,z)]和Var[I(t+τ,z)]分别表示I(t,z)序列和I(t+τ,z)序列的方差。
空间变化序列的协方差和方差分别通过以下公式计算:
时间变化序列的协方差和方差分别通过以下公式计算:
利用图4所得的信号,根据上述公式计算,可以获得不同延迟时间τ时,相关系数g(τ)的变化曲线,如图5所示。
在步骤S103中,计算相关系数随延迟时间的衰减系数,通过衰减系数得到待测血样的扩散系数,利用待测血样的扩散系数与血样粘度之间的关系分析待测血样的凝固过程。
在本发明的一个实施例中,可以测量全血凝固的开始时间、凝固速率、凝块粘度等参数。
通过上述计算得到了光信号在不同时延下的相关系数,分析相关系数随时间变化的衰减系数。在图5中,相关系数随延迟时间的变化关系为近似的指数函数,如下所示:
g(τ)=e-γτ≈1-γτ
因此,当τ接近0时,可以使用线性回归的方法计算衰减系数γ。使用图5所示的区间进行线性拟合,计算衰减系数γ。除线性回归的方法计算衰减系数,本发明的实施例还可以通过指数拟合的方法计算衰减系数。
由于衰减系数γ和样品中粒子的扩散系数D有如下关系:
其中n为折射率,λ为光波长,θ为散射角。而样品中粒子的扩散系数D与样品的粘度V存在以下关系:
其中,K表示玻尔兹曼常数,T表示样品温度,R表示散射粒子半径。当散射角为180°时,将上述两个公式合并,即可得到衰减系数γ与样品粘度V之间的关系,如下所示:
因此,当其他参数不变时,可以通过计算g(τ)曲线的衰减系数γ,表征样品粘度。通过分析衰减系数γ的改变,测量样品粘度V的变化,从而监控血样的凝固过程,分析血样的凝固特性。图6显示了血液凝固过程中的衰减系数变化,随着凝固过程中,血样由液态变为固态,衰减系数降低,对应样品的粘度增加。
在图6中,对衰减系数的变化进行拟合后,得到的拟合曲线可以分为初始平台期、快速变化期和终末稳定期三个时间段。初始平台期的结束时间点可以用来表征血液开始凝固的时间点,快速变化期的时间长度可以用来表征凝固过程的时间长度,终末稳定期的衰减系数值可以用来表征凝固后血块的粘度。
在本发明的一个实施例中,在待测血样中可以添加催化试剂,如氯化钙、高岭土等试剂,改变血样凝固的速度。
在本发明的一个实施例中,对凝固曲线进行归一化分析,消除样品之间的测量误差。
在本发明的一个实施例中,待测血样也可以是血浆样品,或者稀释的全血样品,或者稀释的血浆样品。
在本发明的一个实施例中,还包括:通过扫描光源对待测血样的多个位置或多个待测血样进行照射,采集待测血样的多个位置或多个待测血样随时间变化的光信号,本发明的实施例不限于试管中血液样品的分析,同时可以对活体动物或人体上的在体血样进行检测,也可以对体外管道中的血样进行检测。
综上,如图7所示,本发明的实施例可以包括血液样品成像、延迟时间的信号相关性分析、相关系数随延迟时间的衰减系数计算、衰减系数随凝固时间变化计算、血液凝固属性分析五个过程,通过血液样品的光信号的相关性分析得到血液样品的凝固属性。
根据本发明实施例提出的测量血液凝固属性的方法,通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测,可以用于分析全血或血浆的凝血属性,例如开始凝固时间点、凝固过程时间长度和凝固血块粘度等指标。测量过程中,不需要外力在血样中诱发机械波,只需要光束照射到血样上,即可分析监控血液凝固过程,可以实现在体的无接触血凝监控,无需将血样取出放入样品杯中,可以在凝血部位直接检测血凝过程,在不同样品之间移动成像光束时,可以实现多个血液样品的高通量检测。
其次参照附图描述根据本发明实施例提出的测量血液凝固属性的装置。
图8为根据本发明一个实施例的测量血液凝固属性的装置结构示意图。
如图8所示,该测量血液凝固属性的装置包括:采集模块100、计算模块200和分析模块300。
采集模块100,用于利用光源照射待测血样,通过干涉检测系统采集待测血样随时间变化的光信号。计算模块200,用于计算光信号在不同延迟时间下的相关系数。分析模块300,用于计算相关系数随延迟时间的衰减系数,通过衰减系数得到待测血样的扩散系数,利用待测血样的扩散系数与血样粘度之间的关系分析待测血样的凝固过程。
需要说明的是,前述对测量血液凝固属性的方法实施例的解释说明也适用于该实施例的测量血液凝固属性的装置,此处不再赘述。
根据本发明实施例提出的测量血液凝固属性的装置,通过监测血液样品中散射粒子的随机运动,实现血液凝固属性的检测,可以用于分析全血或血浆的凝血属性,例如开始凝固时间点、凝固过程时间长度和凝固血块粘度等指标。测量过程中,不需要外力在血样中诱发机械波,只需要光束照射到血样上,即可分析监控血液凝固过程,可以实现在体的无接触血凝监控,无需将血样取出放入样品杯中,可以在凝血部位直接检测血凝过程,在不同样品之间移动成像光束时,可以实现多个血液样品的高通量检测。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (8)
1.一种测量血液凝固属性的方法,其特征在于,包括以下步骤:
利用光源照射待测血样,通过干涉检测系统采集所述待测血样随时间变化的光信号;
计算所述光信号在不同延迟时间下的相关系数;
计算所述相关系数随延迟时间的衰减系数,通过所述衰减系数得到所述待测血样的扩散系数,利用所述待测血样的扩散系数与血样粘度之间的关系分析所述待测血样的凝固过程;
所述计算所述光信号在不同延迟时间下的相关系数,包括:
计算所述待测血样的同一空间位置点在不同时间段的相关系数,和/或计算所述待测血样的同一个空间范围内在不同时间点的相关系数;
所述计算所述光信号在不同延迟时间下的相关系数为:
其中,为相关系数,/>为/>时刻采集的信号序列,/>为/>时刻采集的信号序列,/>为延迟时间,/>为所述待测血样的血样深度方向,/>表示/>序列和/>序列的协方差,/>和/>分别表示/>序列和/>序列的方差。
2.根据权利要求1所述的方法,其特征在于,通过所述干涉检测系统对所述待测血样进行照射和采集所述光信号,所述干涉检测系统包括光学相干层析成像系统、马赫-曾德干涉仪、共路干涉仪和迈克尔逊干涉仪,所述光学相干层析成像系统包括基于扫频光源的光学相干层析成像系统和光谱域的光学相干层析成像系统。
3.根据权利要求1所述的方法,其特征在于,所述衰减系数与所述待测血样的扩散系数/>之间的关系为:
其中,为折射率,/>为光波长,/>为散射角;
所述待测血样的扩散系数与血样粘度/>之间的关系为:
其中,表示玻尔兹曼常数,/>表示样品温度,/>表示散射粒子半径。
4.根据权利要求1所述的方法,其特征在于,所述分析所述待测血样的凝固过程包括:分析所述待测血样的全血凝固开始时间、凝固速率、凝块粘度。
5.根据权利要求1所述的方法,其特征在于,所述待测试样包括:血浆样品、稀释的全血样品、稀释的血浆样品。
6.根据权利要求1所述的方法,其特征在于,还包括:在所述待测血样中加入催化试剂,改变所述待测血样的凝固速度。
7.根据权利要求1-6任一项所述的方法,其特征在于,还包括:
通过扫描光源对所述待测血样的多个位置或多个所述待测血样进行照射,采集所述待测血样的多个位置或多个所述待测血样随时间变化的光信号。
8.一种测量血液凝固属性的装置,其特征在于,所述装置用于实现如权利要求1所述的一种测量血液凝固属性的方法,所述装置包括:
采集模块,用于利用光源照射待测血样,通过干涉检测系统采集所述待测血样随时间变化的光信号;
计算模块,用于计算所述光信号在不同延迟时间下的相关系数;
分析模块,用于计算所述相关系数随延迟时间的衰减系数,通过所述衰减系数得到所述待测血样的扩散系数,利用所述待测血样的扩散系数与血样粘度之间的关系分析所述待测血样的凝固过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111088760.9A CN113686735B (zh) | 2021-09-16 | 2021-09-16 | 测量血液凝固属性的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111088760.9A CN113686735B (zh) | 2021-09-16 | 2021-09-16 | 测量血液凝固属性的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113686735A CN113686735A (zh) | 2021-11-23 |
CN113686735B true CN113686735B (zh) | 2023-11-07 |
Family
ID=78586535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111088760.9A Active CN113686735B (zh) | 2021-09-16 | 2021-09-16 | 测量血液凝固属性的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113686735B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6019735A (en) * | 1997-08-28 | 2000-02-01 | Visco Technologies, Inc. | Viscosity measuring apparatus and method of use |
CN1540313A (zh) * | 2003-04-21 | 2004-10-27 | 上海兰杰科学仪器有限公司 | 一种用于血凝分析仪的光路系统 |
CN102640004A (zh) * | 2009-12-04 | 2012-08-15 | 株式会社日立高新技术 | 血液凝固分析装置 |
CN111127418A (zh) * | 2019-12-20 | 2020-05-08 | 苏州大学 | 一种基于oct系统检测漆质涂层固化程度的方法 |
CN111307755A (zh) * | 2020-03-20 | 2020-06-19 | 中国科学院上海高等研究院 | 一种基于太赫兹技术检测液体样品的方法 |
CN111436909A (zh) * | 2020-04-28 | 2020-07-24 | 北京信息科技大学 | 一种活体组织的光学相干层析成像系统及方法 |
CN113092369A (zh) * | 2021-03-23 | 2021-07-09 | 华中科技大学 | 一种用于监测血液凝固动态过程的光学装置及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7821620B2 (en) * | 2007-10-03 | 2010-10-26 | University Of Central Florida Research Foundation, Inc. | Optical coagulation monitor and method of use |
WO2012161960A1 (en) * | 2011-05-26 | 2012-11-29 | The General Hospital Corporation | Optical thromboelastography system and method for evaluation of blood coagulation metrics |
EP2936117B1 (en) * | 2012-12-19 | 2020-11-25 | The General Hospital Corporation | Optical blood-coagulation sensor |
WO2015184433A1 (en) * | 2014-05-30 | 2015-12-03 | The General Hospital Corporation | Optical thromboelastography systems and methods |
CN107315094B (zh) * | 2017-07-27 | 2023-12-01 | 深圳传世生物医疗有限公司 | 凝血分析仪及凝血分析方法 |
-
2021
- 2021-09-16 CN CN202111088760.9A patent/CN113686735B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6019735A (en) * | 1997-08-28 | 2000-02-01 | Visco Technologies, Inc. | Viscosity measuring apparatus and method of use |
CN1540313A (zh) * | 2003-04-21 | 2004-10-27 | 上海兰杰科学仪器有限公司 | 一种用于血凝分析仪的光路系统 |
CN102640004A (zh) * | 2009-12-04 | 2012-08-15 | 株式会社日立高新技术 | 血液凝固分析装置 |
CN111127418A (zh) * | 2019-12-20 | 2020-05-08 | 苏州大学 | 一种基于oct系统检测漆质涂层固化程度的方法 |
CN111307755A (zh) * | 2020-03-20 | 2020-06-19 | 中国科学院上海高等研究院 | 一种基于太赫兹技术检测液体样品的方法 |
CN111436909A (zh) * | 2020-04-28 | 2020-07-24 | 北京信息科技大学 | 一种活体组织的光学相干层析成像系统及方法 |
CN113092369A (zh) * | 2021-03-23 | 2021-07-09 | 华中科技大学 | 一种用于监测血液凝固动态过程的光学装置及方法 |
Non-Patent Citations (5)
Title |
---|
P. Grybauskas.Ultrasonic digital interferometer for investgation of blood clotting.Ultrasonics.1978,第Volume 15卷(第Issue 1期),33-36. * |
Study of the blood coagulation by ultrasound;Chih Chung Huang;Journal of Medical and biological Engineering;第31卷(第2期);79-86 * |
基于吸光度法的全自动血凝仪中光电检测系统的研究;付仁凤;中国优秀硕士学位论文全文数据库医药卫生科技(第2017年第3期期);83 * |
杨栋.《时间序列散斑测量信号处理算法的研究》.北京交通大学.2015,第2015.10卷(第2015.10期),第7-24页. * |
血栓弹力图与传统凝血检查监测围关节置换期的凝血状态;耿玮;张志刚;皮斌;程宇;张磊;王雪明;虞秀兰;周立宇;吴贵忠;王治栋;杨同其;徐明;;中国组织工程研究(第48期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113686735A (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8697449B2 (en) | Optical blood coagulation monitor and method | |
JP5134177B2 (ja) | 電場に基づいた光散乱分光法を用いたシステム | |
US10365254B2 (en) | Assessment of blood coagulation using an acoustic radiation force based optical coherence elastography (ARF-OCE) | |
EP2713878B1 (en) | Optical thromboelastography system and method for evaluation of blood coagulation metrics | |
CA2732334C (en) | A system and method for detection, characterization and imaging of heterogeneity using shear wave induced resonance | |
CN105559756A (zh) | 基于全空间调制谱分割角度复合的微血管造影方法与系统 | |
CN208837916U (zh) | 一种血流成像系统 | |
JP4911540B2 (ja) | 超音波診断装置 | |
Xu et al. | Dynamic and quantitative assessment of blood coagulation using optical coherence elastography | |
Moilanen et al. | Photo-acoustic excitation and optical detection of fundamental flexural guided wave in coated bone phantoms | |
Newson et al. | Laser Doppler velocimetry: the problem of fibre movement artefact | |
WO2010143572A1 (ja) | 被検体情報分析装置及び被検体情報分析方法 | |
Li et al. | Evaluation of blood plasma coagulability by laser speckle correlation | |
CN113686735B (zh) | 测量血液凝固属性的方法及装置 | |
Xu et al. | Viscosity monitoring during hemodiluted blood coagulation using optical coherence elastography | |
CN111436910B (zh) | 一种活体组织的光学相干层析多模态成像装置及方法 | |
Garcia-Duitama et al. | Protocol for robust in vivo measurements of erythrocyte aggregation using ultrasound spectroscopy | |
EP1556672A2 (en) | A method of and apparatus for measuring oscillatory motion | |
CN113092369A (zh) | 一种用于监测血液凝固动态过程的光学装置及方法 | |
Hashimoto et al. | Noncontact estimation of stiffness based on optical coherence elastography under acoustic radiation pressure | |
JP2008170390A (ja) | 血液細胞の力学的特性計測システム | |
JPH10221248A (ja) | 光学的検査方法及び光学的検査装置 | |
RU2779368C1 (ru) | Устройство для измерения параметров гемостаза | |
KR102466497B1 (ko) | 굴절률을 이용한 요산염 결정의 검출방법 및 이의 용도 | |
Zhang et al. | Evaluation of blood coagulation by optical vortex tracking |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |