CN101933812B - 一种尿流动力学检测分析方法 - Google Patents
一种尿流动力学检测分析方法 Download PDFInfo
- Publication number
- CN101933812B CN101933812B CN201010271823XA CN201010271823A CN101933812B CN 101933812 B CN101933812 B CN 101933812B CN 201010271823X A CN201010271823X A CN 201010271823XA CN 201010271823 A CN201010271823 A CN 201010271823A CN 101933812 B CN101933812 B CN 101933812B
- Authority
- CN
- China
- Prior art keywords
- urine
- urethra
- maximum
- urinating
- bladder
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 29
- 238000001514 detection method Methods 0.000 title abstract description 5
- 230000003202 urodynamic effect Effects 0.000 title abstract description 5
- 210000003708 urethra Anatomy 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 42
- 230000008602 contraction Effects 0.000 claims abstract description 30
- 230000001133 acceleration Effects 0.000 claims abstract description 28
- 230000005484 gravity Effects 0.000 claims abstract description 6
- 210000002700 urine Anatomy 0.000 claims description 107
- 230000003187 abdominal effect Effects 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 230000002485 urinary effect Effects 0.000 claims description 7
- 238000000205 computational method Methods 0.000 claims description 4
- QYSXJUFSXHHAJI-YRZJJWOYSA-N vitamin D3 Chemical compound C1(/[C@@H]2CC[C@@H]([C@]2(CCC1)C)[C@H](C)CCCC(C)C)=C\C=C1\C[C@@H](O)CCC1=C QYSXJUFSXHHAJI-YRZJJWOYSA-N 0.000 claims description 4
- 230000001373 regressive effect Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 14
- 208000015181 infectious disease Diseases 0.000 abstract description 5
- 230000027939 micturition Effects 0.000 abstract description 4
- 238000012423 maintenance Methods 0.000 abstract description 2
- 210000003932 urinary bladder Anatomy 0.000 description 34
- 230000006870 function Effects 0.000 description 13
- 238000007689 inspection Methods 0.000 description 10
- 230000000903 blocking effect Effects 0.000 description 8
- 210000003899 penis Anatomy 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 206010046555 Urinary retention Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 201000002327 urinary tract obstruction Diseases 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000036724 intravesical pressure Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000000664 rectum Anatomy 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 206010004446 Benign prostatic hyperplasia Diseases 0.000 description 1
- 206010013554 Diverticulum Diseases 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 206010058914 Hypotonic urinary bladder Diseases 0.000 description 1
- 206010071289 Lower urinary tract symptoms Diseases 0.000 description 1
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 206010034310 Penile pain Diseases 0.000 description 1
- 208000004403 Prostatic Hyperplasia Diseases 0.000 description 1
- 206010055026 Prostatic obstruction Diseases 0.000 description 1
- 206010065584 Urethral stenosis Diseases 0.000 description 1
- 206010046543 Urinary incontinence Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000008828 contractile function Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000002307 prostate Anatomy 0.000 description 1
- 210000003689 pubic bone Anatomy 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 230000001718 repressive effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 210000005070 sphincter Anatomy 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 201000001988 urethral stricture Diseases 0.000 description 1
- 210000001635 urinary tract Anatomy 0.000 description 1
- 238000013022 venting Methods 0.000 description 1
- 230000002747 voluntary effect Effects 0.000 description 1
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/20—Measuring for diagnostic purposes; Identification of persons for measuring urological functions restricted to the evaluation of the urinary system
- A61B5/207—Sensing devices adapted to collect urine
- A61B5/208—Sensing devices adapted to collect urine adapted to determine urine quantity, e.g. flow, volume
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
- G09B23/28—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
- G09B23/30—Anatomical models
- G09B23/303—Anatomical models specially adapted to simulate circulation of bodily fluids
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Business, Economics & Management (AREA)
- Physiology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Medicinal Chemistry (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Urology & Nephrology (AREA)
- Chemical & Material Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种尿流动力学检测分析方法，包括以下步骤：建立膀胱的弹性元模型；建立尿道模型；保持前尿道与重力方向垂直，记录和测量排尿数据，计算排尿参数；计算弹性元的收缩长度，并进一步计算弹性元的收缩速度、收缩加速度和最大收缩加速度；计算尿道模型最大横截面积。本发明可以完全克服传统有创尿流动力学检查方法对患者所造成的痛苦以及感染的可能，整个分析过程都可以在计算机的帮助下自动进行，其结果清晰明了，便于临床记忆和使用，应用此方法所制成的装置结构简单，维护方便，从而降低医疗费用。
Description
技术领域
本发明涉及一种医疗检测方法，特别是一种尿流动力学检测分析方法。
背景技术
尿流动力学是现代医学、生物流体力学、生物流变学的交叉学科，广泛应用于尿路梗阻、尿失禁和其他尿路功能障碍等诸多疾病的基础研究、诊断、治疗和疗效评价，与泌尿外科、妇产科、儿科、内分泌科、神经内科、肛肠科等专业密切相关。
目前，常规尿流动力学直接测量的指标是膀胱和直肠内压力、尿流率、尿量。最常用的是有创尿流动力学检查方法：插入测压管直接测量膀胱内的压力变化，同时测定尿流率；恒速牵引恒流灌注的测压管测量尿道各段压力分布情况；利用AG、PQ图进行分析，诊断疾病。此方法最大限度保证了测量结果与真实情况的拟合度，但检查是在非生理状态下进行，有可能产生假阳性结果，而且侵入性的检查增加患者痛苦以及感染的机率。
除了力学分析外，还有学者从能量的角度进行研究，如排尿过程能量消耗方程F=9.79×10-2P-1.25×10-3Q2+9.8H。由Bernolli粘性液体微元流速方程推导而来的排尿过程能量消耗方程，取膀胱出口、尿道外口为研究平面，流速的计算是用尿流率除以尿道外口的横截面（标准的尿道外口横截面积）。但我们知道Bernolli方程反映的是液流中机械能和其他形式的能量（主要是代表能量损失的热能）间的守恒与转化关系，应用条件是：①定常流动的不可压缩流体，微元流管的形状不随时间而改变；②质量力只有重力，两过水断面间除了水头损失以外，总流没有能量的输入或输出；③所选取的两过水断面必须是缓变流断面，但两过水断面间可以是急变流；④总流的流量沿程不变。从希尔定律可以知道，不同膀胱容量的排尿压力是不同的，也就是说有能量的输入并随时间改变；尿液排泄是以尿道为运动的边界，尿道的形状会随着时间在排尿过程中发生改变，尿流运动是非定常的，尿流动力学问题是流体运动与边界变形运动的耦合，适合使用反映液流与边界上作用力之间关系的动量方程。因此，如果直接将Bernolli粘性液体微元流速方程应用于下尿路尿流动力学的研究，舍弃其他影响因素，将不能准确反映真实的情况。
文献报道了基于光电原理的尿滴分光仪（Drop spectrometer），调制的光线通过水平缝隙，再被聚焦到检测的二极管上，消失的光与尿滴通过光幕的阴影有关。高速摄影机观察可见尿线是由尿滴组成，通过检测脉冲的变化，可以确定每一滴尿液的时间和体积。这项技术可以检测有无梗阻但不能有效判断接近膀胱的梗阻。类似的还有利用CCD作为检测光线变化的发明。
现有技术中，美国专利第5,377,101号，公开了一种用计算机分析尿流率曲线的方法，将检测的尿流率曲线与标准曲线对比，判断是否有异常。T90——排90%尿量的时间，在尿流率曲线上取5%—95%区间测量。QM90——90%尿量的平均尿流率。Tdesc——下降曲线的时间，从最大尿流率开始到95%尿量的时间。Dl/dt40——评价膀胱在40ml时的收缩速度。该方法不能充分描述数量相关、复杂、非线性的流体系统，阻力因素有括约肌、前列腺和远端尿道，而且膀胱和尿道有蠕变性。因为压力和尿流率都是非线性的，所以这种方法还是不能鉴别逼尿肌无力、前列腺增生、尿道狭窄。
现有技术中，文献公开了一种基于流体力学测压管原理的套状尿管检查方法和器械，用阴茎套装的尿管套于阴茎上，夹闭出口后测压管压力上升，当尿流不动时充满尿道的尿液可作为尿管测量膀胱内的压力。虽然不用插入测压管，但仍有逆行感染的可能。此方法不能观察尿流全部过程，而且阴茎套的材料特性难以完全一致，检测结果可重复性差。
现有技术中，美国专利第5,807,278和5,823,972号，公开了一种利用测血压原理测量排尿期的压力，将幼儿的血压计袖套套在阴茎上充气，当患者开始排尿后慢慢放气。阴茎套充气试验：排尿开始后以一定的速率对阴茎套充气。阴茎压迫和放松：排尿时握紧阴茎切断尿流3秒，放松时储存在尿道的尿液建立起伏的尿流率。对比尿流率起伏的大小可以很好地预测流出道的梗阻。充气是会引起阴茎疼痛，而且测量结果受袖套放置位置、大小、形状和材料特性影响，标准化困难。
文献报道通过多普勒估计尿道内的尿流速度。在一小组病例中，发现与侵入性检查有很好的拟合度，如前列腺梗阻和交界部位（通过速度和流率计算），尿道前列腺部和膜部的速率。但实验是在流速大于2ml/s的玻璃管道中检测到的结果，而且处理过程复杂。
文献报道利用尿液通过尿道前列腺部时变成紊流产生声音，在会阴部记录声学特征评价梗阻的程度。在用阴茎套建立的梗阻模型中，发现不同的梗阻程度有不同的声谱。这种方法对检查环境要求高，数据处理复杂，而且放置在会阴部的探头压迫尿道海绵体，影响测量结果。
文献报道通过测量膀胱的重量和厚度预测男性的梗阻。在膀胱容量为250ml时经腹B超检测逼尿肌厚度，预测流出道梗阻。此方法不能鉴别广泛小樑、憩室形成，膀胱壁增厚，膀胱收缩无力的病例。
现有技术中，文献公开了一种近红外光谱法检查方法：在患者耻骨上放置探头，通过红外线探知逼尿肌血红蛋白的变化来推测膀胱压力的变化。此仪器价格贵，临床使用效果有待观察。
发明内容
为了克服现有技术使用时，有创检查方法易引起患者疼痛或感染，无创检查方法不准确、价格贵等不足，本发明提供一种尿流动力学检测分析方法，测量流出尿道外口的尿液质量、速度，根据拓扑学原理建立膀胱弹性元模型、尿道模型，利用能量守恒定律、流体动量方程分析膀胱收缩功能、尿道阻力情况，适用于下尿路症状患者检查以及手术、药物疗效的评价。
本发明为解决其问题所采用的技术方案是：
一种尿流动力学检测分析方法，包括以下步骤：
（1）建立膀胱的弹性元模型：将排尿前膀胱内的尿液拓扑成一个球体，膀胱逼尿肌和腹压的作用看成弹性元的作用，弹性元的模量可随时间、空间改变，根据弹性元模型确定弹性元长度L与膀胱尿量a的函数L =F(a)以及收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)；
（2）建立尿道模型：将尿道视为水平放置的管道，其长度、横截面积随时间发生瞬间改变，尿液在尿道模型中动量、能量的改变等于尿道对尿液的作用；
（3）保持前尿道与重力方向垂直，记录和测量排尿数据，该排尿数据包括排尿原点位置、排尿高度h、不同时间点t i 上排出尿液的重量m i 、水平位移d i ，计算排尿参数，该排尿参数包括不同时间点t i 上的尿流率Q i 、尿流速度vu-i 、单位时间排出尿液的动能E i 和尿道模型横截面积S i ；
（4）根据收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)计算弹性元的收缩长度ΔLi ，并进一步计算弹性元的收缩速度vi 和收缩加速度a d-i ，以弹性元的最大收缩加速度a d-max评价逼尿肌收缩功能；
（5）通过尿道模型横截面积S i 判断尿道阻力情况，通过最大横截面积S max 判断尿道有无梗阻。
本发明的有益效果是：本发明通过外部测量和数学分析的方法对患者尿道的梗阻情况以及逼尿肌的收缩功能进行评价，其可以完全克服传统有创尿流动力学检查方法对患者所造成的痛苦以及感染的可能，整个分析过程都可以在计算机的帮助下自动进行，其结果清晰明了，便于临床记忆和使用，应用此方法所制成的装置结构简单，维护方便，从而降低医疗费用。
附图说明
下面结合附图和实施例对本发明作进一步说明：
图1为膀胱出口梗阻男性患者的尿流率图；
图2为膀胱出口梗阻男性患者的逼尿肌收缩加速度曲线图；
图3为逼尿肌力正常的下尿路梗阻男性患者的最大收缩加速度的频数直方图；
图4为本发明最大动能——收缩长度分析方法的流程方框图；
图5为本发明逼尿肌收缩加速度、尿道横截面积分析方法的流程方框图。
具体实施方式
参照图4和图5，本发明的一种尿流动力学检测分析方法，包括以下步骤：
（1）建立膀胱的弹性元模型：以逼尿肌为研究对象，生物流变学研究表明膀胱逼尿肌有自发收缩能力，无法将张力分解为被动张力、主动张力，逼尿肌弹性元的模量是纤维长度的函数，在给定长度下又是时间的函数，可完全松弛至应力为零，且在静息态下不存在单一的状态，因此不适合使用Hill模型进行研究，而影像尿流动力学显示排尿过程中，在合适尿量下膀胱呈球形，因此可以根据拓扑学将排尿前膀胱内的尿液拓扑成一个球体，膀胱逼尿肌和腹压的作用看成弹性元的作用，弹性元的模量可随时间、空间改变，根据弹性元模型确定弹性元长度L与膀胱尿量a的函数L =F(a)以及收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)，一般地，球体的体积等于膀胱尿量a，弹性元的初始长度Li 等于以球体中心为圆心的周长，相当于膀胱一周逼尿肌纤维的长度，球体体积变化就是尿流率Q，根据上述关系可以根据拓扑函数建立弹性元长度L与膀胱尿量a的函数，为L =F(a) =2*(π2*a)1/3≈8.36084644*a1/3，而收缩长度ΔL与尿流率Q、膀胱尿量a的函数为ΔL=ξ(Q,a)=
Q/{[a1/3+(a－Q) 1/3]2－a1/3*(a－Q) 1/3}；
（2）建立尿道模型：将尿道视为水平放置的管道，其长度、横截面积随时间发生瞬间改变，尿液在尿道模型中动量、能量的改变等于尿道对尿液的作用，生物流变学研究表明尿道具有弹性滞后、应力松弛、蠕变和Boltzmann叠加的特点，不同部位的尿道腔横截面积不相等，而且排尿过程中，尿道腔随时间、尿流动能发生非线性的变化，表现为阻力随时间的非线性改变，当时间无穷小时，阻力只有一个值，那么，可以根据能量守恒定律，建立产生相同阻力的水平管道模型，推导出尿道模型横截面积S i 的计算方法为尿流率除以尿流速度，S i =Q i /vu-i ；
（3）保持前尿道与重力方向垂直，记录和测量排尿数据，该排尿数据包括排尿原点位置、排尿高度h、不同时间点t i 上排出尿液的重量m i 、水平位移d i ，计算排尿参数，该排尿参数包括不同时间点t i 上的尿流率Q i 、尿流速度vu-i 、单位时间排出尿液的动能E i 和尿道模型横截面积S i ，其中某时间点上的尿流率Q i 可以通过计算(m i -m i-1 )/
(t i-t i-1 )求得，尿流速度vu-i 的计算公式为vu-i =d i /t ，t的大小与排尿高度直接相关，t=(2h/g)1/2，单位时间排出尿液的动能E i =(m i -m i-1 )
(vu-i )2*0.5, 尿道模型横截面积S i =Q i /vu-i ，当然，除了上述所提到的排尿数据和排尿参数外，本发明还可以进一步对其他数据和参数进行记录和计算，例如总排尿时间、总尿流时间、平均尿流率、残余尿量（可于排尿完成后通过B超进行测量）、平均尿流速度、膀胱最大收缩速度及时间、总输出功率、平均输出功率、最大输出功率及时间、最大功率加速度、总动量、总动能等等，此外，根据上述得到的各种数据和参数还可以进一步制成尿流率曲线、尿流速度曲线、弹性元收缩速度曲线、输出功率曲线、动量曲线、尿道模型横截面积曲线等等，通过与正常范围曲线的对比进行后续的分析研究；
（4）根据收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)计算弹性元的收缩长度ΔLi ，并进一步计算弹性元的收缩速度vi 和收缩加速度a d-i ，其中收缩速度vi=(ΔLi-ΔLi-1 )/ (t i-t i-1
)，收缩加速度a d-i =(vi-vi-1 )/
(t i-t i-1 )，以弹性元的最大收缩加速度a d-max评价逼尿肌收缩功能，弹性元的最大收缩加速度a d-max的计算方法为：作时间-收缩速度图，取从排尿开始至出现最大尿流速度、最大尿流率和最大尿道模型横截面积的时间段，用回归的方法计算收缩速度的斜率，以该斜率作为最大收缩加速度a d-max，逼尿肌收缩功能的评价方法为：将最大收缩加速度与同年龄组的正常值进行比较，若在正常值范围内，则判断逼尿肌收缩功能正常；若小于正常值，则判断逼尿肌收缩功能有受损可能性，需要进一步进行检查；
（5）通过尿道模型横截面积S i 判断尿道阻力情况，通过最大横截面积S max 判断尿道有无梗阻，由于尿道具有弹性滞后、应力松弛的特点，也可取2～3个相连的横截面积最大值的平均值作为最大横截面积S max ，尿道有无梗阻的判断方法为：
①将最大横截面积与同年龄组的正常值进行比较，若在正常范围内，则判断尿道无梗阻，计算最大尿道半径，进入⑤，若小于正常值，则判断尿道出现梗阻，并进入②；
②绘制时间-尿道横截面积变化曲线，若曲线波幅大，则判断逼尿肌-括约肌失协调或者其他干扰因素；
③ 判断时间-尿道横截面积变化曲线是否出现长平直段：是，则判断为缩窄性梗阻；否，则判断为压迫性梗阻；
④计算最大理论尿道半径；
⑤输出结果。
上述测量过程中水平位移d i 的测量可以采用多种不同的方法，其中一种较优的测量方法为：
①读取排尿的落点图像数据图片；
②将采用的每一帧图片与上一帧图片对比，落点图像重合则忽略；
③如有区别，记录时间点，然后采用最大方差阈值法进行图像二值化，再进行边界搜索，计算落点；
④根据排尿原点位置为圆心的同心圆系统标定，计算出原点位置与落点位置的水平位移d i 。
为了更好地根据所测量的数据对膀胱收缩功能、尿道阻力情况进行分析，本发明还优选包括于步骤（5）之后的下列步骤：
（6）取单位时间排出尿液的动能最大值E imax ，判断排尿时是否出现腹压变化：判断腹压没有变化时，上述动能最大值作为最大动能；判断腹压出现变化时，上述动能最大值减去腹压做功值后作为最大动能，是否出现腹压变化的判断方法为：将收缩加速度a d-i 绝对值与正常值比较，如收缩加速度a d-i 绝对值在开始时间T0和结束时间TEND（T0和TEND的时长都为1S左右）以外的时间段大于正常值的超过2个，且出现于尿流率Q i 的最大值Q max 之前，在排除人为干扰或仪器赝象的情况下，可判断为出现腹压变化；
（7）以最大动能及对应的收缩长度ΔLi 绘制最大动能-收缩长度图，根据最大动能-收缩长度图判断逼尿肌肌力、梗阻等级。
考虑到出现腹压变化会对分析结果产生较大的影响，为了避免分析结果不准确，实际应用时优先选择没有出现腹压变化的数据进行分析计算。
上述分析步骤中所用到的收缩加速度a d-i 绝对值的正常值以及最大收缩加速度a d-max的正常值都通过统计的方式得出，下面以最大收缩加速度a d-max为例进行说明，如图3所示，为逼尿肌力正常的下尿路梗阻男性患者的最大收缩加速度的频数直方图，其样本为361例，通过对直方图的数据进行统计可以得出其最大收缩加速度-0.0022～0.104，算术均数0.013160665±0.0124720779，中位数0.01016，实际应用时，通过检测数据与上述统计结果进行比较即可判断和分析。
下面以一具体实施例对本发明的检测分析方法进行说明：
检查对象为一名60岁的膀胱出口梗阻男性患者，检查时采用以下步骤进行：
（1）保持前尿道与重力方向垂直（女性患者取半坡卧位，男性患者取立位）；
（2）通电，系统初始化，粘贴腹部电极或安置直肠测压管（调零）；
（3）按开始按钮，患者自行排尿；
（3）系统按一定的频率同时采集尿液质量、水平位移图像，经缓存区传入电脑；
（4）结束排尿后，按结束按钮，标定一次排尿过程；
（5）B超测剩余尿量，输入数据；
（6）人工或微机计算、分析；
（7）输出结果、图表、结论。
检查所得的图表如图1、图2所示。
图1是尿流率图。图形显示：曲线正常，尿量为343ml，最大尿流率超过20ml/s。仅从此图判断，患者排尿正常，但残余尿量达1000ml。这是由于膀胱过度扩张，被动张力起主导作用，出现正常的尿流率。
图2可见逼尿肌收缩加速比除了开始的3s其余的均小于0.2，与排尿过程未见腹压波动相吻合，加速度曲线的非线性反映了逼尿肌活动的非线性特征。膀胱尿量1343ml，用回归方法计算收缩速度的回归方程斜率—逼尿肌最大收缩加速度是0.45%，这一数值大大小于最大收缩加速度的统计中位数0.01016，表明该患者存在膀胱收缩功能受损情况，建议进行进一步诊治。
本发明的临床检测分析结果的准确度如下表所示：
其样本为396例，与金标准（有创的尿流动力学分析）对照，本方法的敏感性Se=9/12=75%，特异性Sp=363/384=94.53%，检测准确率T=372/396=93.94%，本发明的检测分析方法具有极高的准确率。
Claims (7)
1. 一种尿流动力学检测分析方法，其特征在于包括以下步骤：
（1）建立膀胱的弹性元模型：将排尿前膀胱内的尿液拓扑成一个球体，膀胱逼尿肌和腹压的作用看成弹性元的作用，弹性元的模量可随时间、空间改变，根据弹性元模型确定弹性元长度L与膀胱尿量a的函数L =ƒ(a)以及收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)；
（2）建立尿道模型：将尿道视为水平放置的管道，其长度、横截面积随时间发生瞬间改变，尿液在尿道模型中动量、能量的改变等于尿道对尿液的作用；
（3）保持前尿道与重力方向垂直，记录和测量排尿数据，该排尿数据包括排尿原点位置、排尿高度h、不同时间点t i 上排出尿液的重量m i 、水平位移d i ，计算排尿参数，该排尿参数包括不同时间点t i 上的尿流率Q i 、尿流速度vu-i 、单位时间排出尿液的动能E i 和尿道模型横截面积S i ；
（4）根据收缩长度ΔL与尿流率Q、膀胱尿量a的函数ΔL=ξ(Q,a)计算弹性元的收缩长度ΔLi ，并进一步计算弹性元的收缩速度vi 和收缩加速度a d-i ，计算弹性元的最大收缩加速度a d-max；
（5）计算尿道模型最大横截面积S max 。
2.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于还包括于步骤（5）之后的下列步骤：
（6）将收缩加速度a d-i 绝对值与正常值比较，如收缩加速度a d-i 绝对值在开始时间T0和结束时间TEND以外的时间段大于正常值的超过2个，且出现于尿流率Q i 的最大值Q max 之前，判断为出现腹压变化；
（7）取单位时间排出尿液的动能最大值E imax ，以最大动能及对应的收缩长度ΔLi 绘制最大动能-收缩长度图，腹压变化时，动能最大值E imax 为计算所得值减去腹压做功值后所得值；
（8）选取最大动能-收缩长度图中没有出现腹压变化时的部分进行分析。
3.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于弹性元长度L与膀胱尿量a的函数为L =ƒ(a) =2*(π2*a)1/3≈8.36084644*a1/3。
4.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于收缩长度ΔL与尿流率Q、膀胱尿量a的函数为ΔL=ξ(Q,a)= Q /{[a1/3+(a－Q)
1/3]2－a1/3*(a－Q) 1/3}。
5.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于水平位移d i 的测量方法为：
①读取排尿的落点图像数据图片；
②将采用的每一帧图片与上一帧图片对比，落点图像重合则忽略；
③如有区别，记录时间点，然后采用最大方差阈值法进行图像二值化，再进行边界搜索，计算落点；
④根据排尿原点位置为圆心的同心圆系统标定，计算出原点位置与落点位置的水平位移d i 。
6.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于尿道模型横截面积S i 的计算方法为尿流率除以尿流速度，S i = Q i/vu-i 。
7.根据权利要求1所述的一种尿流动力学检测分析方法，其特征在于弹性元的最大收缩加速度a d-max的计算方法为：作时间-收缩速度图，取从排尿开始至出现最大尿流速度、最大尿流率和最大尿道模型横截面积的时间段，用回归的方法计算收缩速度的斜率，以该斜率作为最大收缩加速度a d-max。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010271823XA CN101933812B (zh) | 2009-09-16 | 2010-09-03 | 一种尿流动力学检测分析方法 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910190862.4 | 2009-09-16 | ||
| CN200910190862 | 2009-09-16 | ||
| CN201010271823XA CN101933812B (zh) | 2009-09-16 | 2010-09-03 | 一种尿流动力学检测分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101933812A CN101933812A (zh) | 2011-01-05 |
| CN101933812B true CN101933812B (zh) | 2012-08-29 |
Family
ID=43387549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201010271823XA Expired - Fee Related CN101933812B (zh) | 2009-09-16 | 2010-09-03 | 一种尿流动力学检测分析方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20120179387A1 (zh) |
| CN (1) | CN101933812B (zh) |
| WO (1) | WO2011032477A1 (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9909598B1 (en) * | 2014-02-24 | 2018-03-06 | Landtec North America, Inc. | Well monitoring and pressure controlled landfill pump |
| CN111657863B (zh) * | 2014-09-28 | 2023-10-20 | 波特雷罗医疗公司 | 用于感测生理数据以及引流和分析体液的系统、设备和方法 |
| EP3657139A1 (en) | 2014-11-24 | 2020-05-27 | ART Healthcare Ltd. | Drop and drip measurement |
| US10638966B2 (en) * | 2015-08-05 | 2020-05-05 | Art Healthcare Ltd. | Point of care urine analyzer |
| WO2017079380A1 (en) | 2015-11-03 | 2017-05-11 | Franco Intelligent Agent Solutions, Llc | Urodynamic device and procedure |
| CN109805946A (zh) * | 2019-02-27 | 2019-05-28 | 杨小华 | 一种应用尿道阻力测定判别下尿路梗阻程度的方法 |
| CN113440137B (zh) * | 2021-08-02 | 2024-01-26 | 天津市儿童医院 | 一种无创膀胱压力测定装置及其测定方法 |
| CN113712598B (zh) * | 2021-09-09 | 2023-07-25 | 天津理工大学 | 一种便携式膀胱尿量监测系统及方法 |
| CN114659768B (zh) * | 2022-03-16 | 2023-09-19 | 广东工业大学 | 一种尿道阀生物力学模拟实验方法及装置 |
| CN116432564B (zh) * | 2023-06-15 | 2023-08-18 | 天津市第五中心医院 | 一种尿动力学状态分析方法及分析系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19617854A1 (de) * | 1996-04-23 | 1997-11-06 | Viomed Vertriebs Gmbh | Verfahren zur urodynamischen Diagnostik bei Inkontinenz und Meßgerät |
| CN1867294A (zh) * | 2003-10-15 | 2006-11-22 | 不列颠哥伦比亚大学 | 用于尿动力学分析的方法和设备 |
| CN1942140A (zh) * | 2004-02-11 | 2007-04-04 | 伊西康公司 | 利用微电子机械系统进行尿动力学评价的系统和方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060122488A1 (en) * | 2004-11-18 | 2006-06-08 | Abdol-Mohammad Kajbafzadeh | Urodynamic diagnostic method and system |
-
2010
- 2010-09-03 CN CN201010271823XA patent/CN101933812B/zh not_active Expired - Fee Related
- 2010-09-13 WO PCT/CN2010/076835 patent/WO2011032477A1/zh active Application Filing
- 2010-09-13 US US13/496,841 patent/US20120179387A1/en not_active Abandoned
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19617854A1 (de) * | 1996-04-23 | 1997-11-06 | Viomed Vertriebs Gmbh | Verfahren zur urodynamischen Diagnostik bei Inkontinenz und Meßgerät |
| CN1867294A (zh) * | 2003-10-15 | 2006-11-22 | 不列颠哥伦比亚大学 | 用于尿动力学分析的方法和设备 |
| CN1942140A (zh) * | 2004-02-11 | 2007-04-04 | 伊西康公司 | 利用微电子机械系统进行尿动力学评价的系统和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2011032477A1 (zh) | 2011-03-24 |
| US20120179387A1 (en) | 2012-07-12 |
| CN101933812A (zh) | 2011-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101933812B (zh) | 一种尿流动力学检测分析方法 | |
| Scott et al. | Studies on the dynamics of micturition: observations on healthy men | |
| Abrams et al. | Urodynamic techniques | |
| Griffiths | Urodynamics | |
| CN101664313B (zh) | 一种膀胱储尿排尿生理过程的模拟和分析装置 | |
| US20100234758A1 (en) | Valsalva Lung Pressure Monitoring System and Method | |
| Churchill et al. | Urodynamics | |
| CN102599892A (zh) | 腹腔内压力监测系统及其监测方法 | |
| KR20200104827A (ko) | 하부요로 기능이상의 진단 시스템 및 방법 | |
| CN102631191A (zh) | 膀胱压力和腹腔压力的监测系统及方法 | |
| Drach et al. | Disposable peak urinary flowmeter estimates lower urinary tract obstruction | |
| CN202223235U (zh) | 膀胱安全容量检测装置 | |
| Petros et al. | A mathematical model for micturition gives new insights into pressure measurement and function | |
| CN103126666B (zh) | 一种尿流动力学分析仪用的自动检测转换装置 | |
| CN107432747A (zh) | 一种在人体排尿过程中测量尿道阻力能量消耗的方法 | |
| CN111801049A (zh) | 尿动力学检查设备、系统及方法 | |
| CN205083471U (zh) | 一种改进的腹腔压力监测装置 | |
| Lahk et al. | Balloon topography: a simple method of evaluating anal function | |
| JP2007075383A (ja) | 無侵襲超音波膀胱内圧測定器 | |
| CN113440137A (zh) | 一种无创膀胱压力测定装置及其测定方法 | |
| Versi | Relevance of urethral pressure profilometry to date | |
| CN107198518A (zh) | 一种动态膀胱压力容积监测装置 | |
| Lotti | Innovative Technologies In Non-Invasive Urodynamics Diagnostic | |
| Thind et al. | Method for evaluation of the urethral closure mechanism in women during standardised changes of cross-sectional area | |
| Vignoli et al. | Noninvasive urodynamics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120829 |