CN106446504B - 用于生物医学模拟的系统和方法 - Google Patents

用于生物医学模拟的系统和方法 Download PDF

Info

Publication number
CN106446504B
CN106446504B CN201610608292.6A CN201610608292A CN106446504B CN 106446504 B CN106446504 B CN 106446504B CN 201610608292 A CN201610608292 A CN 201610608292A CN 106446504 B CN106446504 B CN 106446504B
Authority
CN
China
Prior art keywords
simulation
crt
patient
improvement rate
biomedical
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
Application number
CN201610608292.6A
Other languages
English (en)
Other versions
CN106446504A (zh
Inventor
门冈良昌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Publication of CN106446504A publication Critical patent/CN106446504A/zh
Application granted granted Critical
Publication of CN106446504B publication Critical patent/CN106446504B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/10—Office automation; Time management
    • 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
    • G06F19/324—
    • G06F19/328—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/22—Social work or social welfare, e.g. community support activities or counselling services
    • 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
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00—Electrotherapy; Circuits therefor
    • A61N1/18—Applying electric currents by contact electrodes
    • A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
    • A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
    • A61N1/362—Heart stimulators
    • A61N1/3627—Heart stimulators for treating a mechanical deficiency of the heart, e.g. congestive heart failure or cardiomyopathy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00—Electrotherapy; Circuits therefor
    • A61N1/18—Applying electric currents by contact electrodes
    • A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
    • A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
    • A61N1/362—Heart stimulators
    • A61N1/365—Heart stimulators controlled by a physiological parameter, e.g. heart potential
    • A61N1/368—Heart stimulators controlled by a physiological parameter, e.g. heart potential comprising more than one electrode co-operating with different heart regions
    • A61N1/3684—Heart stimulators controlled by a physiological parameter, e.g. heart potential comprising more than one electrode co-operating with different heart regions for stimulating the heart at multiple sites of the ventricle or the atrium
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00—Electrotherapy; Circuits therefor
    • A61N1/18—Applying electric currents by contact electrodes
    • A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
    • A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
    • A61N1/362—Heart stimulators
    • A61N1/365—Heart stimulators controlled by a physiological parameter, e.g. heart potential
    • A61N1/368—Heart stimulators controlled by a physiological parameter, e.g. heart potential comprising more than one electrode co-operating with different heart regions
    • A61N1/3684—Heart stimulators controlled by a physiological parameter, e.g. heart potential comprising more than one electrode co-operating with different heart regions for stimulating the heart at multiple sites of the ventricle or the atrium
    • A61N1/36843—Bi-ventricular stimulation
    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Pathology (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Databases & Information Systems (AREA)
  • Biomedical Technology (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Urology & Nephrology (AREA)
  • Child & Adolescent Psychology (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

本发明公开了一种用于生物医学模拟的系统和方法。生物医学模拟系统基于心脏射血分数和心电图数据来确定是否执行生物医学模拟，并且在确定执行生物医学模拟时，基于层析图像来生成患者的心脏的生物医学模型。接下来，生物医学模拟系统基于生物医学模型来执行对患者执行手术之后患者的心脏的术后状况的生物医学模拟，并且确定生物医学模拟的结果是否满足与症状相关联的预定改善率。生物医学模拟系统在生物医学模拟的结果不满足预定改善率时确定不偿付针对患者的手术的医学治疗成本。当生物医学模拟的结果满足所述预定改善率时，生物医学模拟系统决定对患者执行手术。

Description

用于生物医学模拟的系统和方法
技术领域
本文中讨论的实施例涉及用于生物医学模拟的系统和方法。
背景技术
由于心肌梗死、心衰等，用于心脏激励的传导通路中可能会出现异常。当传导路径中发生异常时，室间隔与心室游离壁之间的收缩时序失去同步，这导致下降的心脏射血分数(EF)和收缩压。射血分数是通过每次心搏过程中由心脏泵出的血液的体积(心搏出量)除以舒张过程中左心室内的血液的体积来计算的。随着患者的症状——如下降的射血分数——的发展，有时会出现患者难保持正常生活的情况。
对于这样的患有严重的心脏疾病的患者，基于可行的医疗准则，在体内植入心脏再同步化治疗(CRT)设备，心脏再同步化治疗(CRT)设备是一种起搏器。准则阐述了例如射血分数≤35％且心电图(ECG)QRS持续时间≥130毫秒的患者作为CRT植入的对象。QRS持续时间是从ECG上的QRS复合波(表示心室激动)中Q波的起始到S波的结束(offset)的时间。
CRT设备具有附接至其上的三条导电线(lead wire)，并且电极被设置在每条导电线的端部上或中间。导电线中的两条导电线被植入到患者的右心房和右心室中，并且第三导电线被布置在冠状窦中。然后，CRT设备以使心室和室间隔与患者中的心脏起搏同步收缩的方式来向导电线上的电极施加电势。
在对通过CRT设备植入产生的双心室起搏做出响应的患者的情况下，心脏运动被重新同步，使得心脏以有节奏的且协调的方式搏动，从而导致改善的心脏泵功能。另一方面，仍有约30％的患者不能对通过CRT设备植入产生的双心室起搏做出响应(称为“不响应者”)。此外，医生通过手术过程中的反复试验和错误来识别最优电极位置，导致外科手术时间延长。这对患者尤其是不响应者造成了相当大的身体上或心理上的负担。鉴于这样的问题，已经提出了与CRT应用确定相关联的各种技术，例如，用于选择适于CRT的患者的技术和用于识别起搏器的电极的最优安装位置的技术。
日本公开专利公报No.2011-98175
国际专利申请的日本国家公开No.2008-534165
国际专利申请的日本国家公开No.2010-528683
Jason Constantino,Yuxuan Hu,Natalia A.Trayanova,"A computationalapproach to understanding the cardiac electromechanical activation sequencein the normal and failing heart,with translation to the clinical practice ofCRT",Progress in Biophysics and Molecular Biology,2012 Oct-Nov,110372-379
Sermesant M,Chabiniok R,Chinchapatnam P,Mansi T,Billet F,Moireau P,Peyrat JM,Wong K,Relan J,Rhode K,Ginks M,Lambiase P,Delingette H,Sorine M,Rinaldi CA,Chapelle D,Razavi R,Ayache N.,"Patient-specific electromechanicalmodels of the heart for the prediction of pacing acute effects in CRT:Apreliminary clinical validation"Medical Image Analysis,2012 Jan；16(1):201-215
Okada J,Sasaki T,Washio T,Yamashita H,Kariya T,Imai Y,Nakagawa M,Kadooka Y,Nagai R,Hisada T,Sugiura S.,"Patient Specific Simulation of BodySurface ECG using the Finite Element Method"Pacing and ClinicalElectrophysiology,2013 Mar；36(3):309-321
CRT设备非常昂贵，并且CRT植入是侵入性的。因此，适于CRT植入的患者中包括不响应者导致浪费的医疗费用和增加的保险成本。例如在美国，医学治疗尤其昂贵，并且有时的情况是即使患者有健康保险，他们仍然要负责一部分的医学治疗费用。如果在CRT植入后患者被证实是不响应者，则医生会面临由被迫承担医疗费用的患者或保险公司提出的寻求巨额资金的损失赔偿的诉讼的风险。医生作为诉讼的败诉方必须承担医学治疗成本。因此，在CRT中，包括患者、医生和保险公司的多方承担支付成本的风险。
另一方面，基于患者自身的准确的心脏模型来运行心脏模拟使得能够在执行CRT植入之前对最优电极位置和CRT的有效性进行准确的估计。CRT有效性的准确估计有助于减少不响应者中的CRT植入的事件，这导致医疗费用的减少。
由承担支付医疗费用的风险的多方做出的适当判断避免执行无效的治疗。但是，不存在使得风险承担方能够共享心脏模拟的结果并且反映实施适当的治疗的决策的结果的现有系统。因此，CRT甚至用于治疗并不期望CRT对其有效的患者，从而增加了医疗费用。
除了CRT之外，存在能够通过使用复制患者的器官的模型的生物医学模拟来预测患者的预后的治疗方法。由于归因于不存在从生物医学模拟获得的结果的适当共享的无效治疗的实践而产生的增加的医疗费用的问题并不仅存在于CRT中，其对于其他医学治疗而言也是普遍的。
发明内容
下文所描述的实施例的一个方面旨在减少医疗费用。
根据一个方面，提供了一种生物医学模拟系统，其包括：生物医学信息存储装置，其用于存储特定于患者的生物医学信息，生物医学信息包括患者的心脏的层析图像、心脏射血分数和心电图数据；模拟执行确定装置，其用于基于心脏射血分数和心电图数据来确定是否执行生物医学模拟；模型生成装置，其用于在模拟执行确定装置确定执行生物医学模拟时，基于层析图像来生成患者的心脏的生物医学模型；模拟执行装置，其用于基于由模型生成装置生成的生物医学模型来执行对患者执行手术之后心脏的术后状况的生物医学模拟；改善率确定装置，其用于确定由模拟执行装置执行的生物医学模拟的结果是否满足与症状相关联的预定改善率；治疗费用偿付装置，其用于在改善率确定装置确定生物医学模拟的结果不满足预定改善率时确定不偿付针对患者的手术的医学治疗成本；以及手术决策装置，其用于在改善率确定装置确定生物医学模拟的结果满足预定改善率时决定对患者执行手术。
附图说明
图1示出了根据第一实施例的生物医学模拟系统的配置示例；
图2示出了CRT设备植入的示例；
图3示出了根据第二实施例的系统配置的示例；
图4是示出用于通过使用CRT模拟来确定是否提供CRT的过程的流程图；
图5示出了医院信息系统的硬件配置的示例；
图6是示出各个系统的功能的框图；
图7示出了根据第二实施例的CRT治疗和医疗费用偿付的过程的示例；
图8是示出确定动态系统模拟参数的处理过程的示例的流程图；
图9是示出心脏模拟的过程的示例的流程图；
图10是示出确定电系统模拟参数的处理过程的示例的流程图；
图11是示出参数调整的处理过程的示例的流程图；
图12是示出计算电势的处理过程的示例的流程图；
图13是示出针对每个CRT设备的心脏模拟的处理过程的示例的流程图；
图14示出了被设计用于双心室起搏CRT设备的心脏模拟的过程的示例；
图15示出了被设计用于三心室(triventricular)起搏CRT设备的心脏模拟的过程的示例；
图16示出了被设计用于四心室(quadriventricular)起搏CRT设备的心脏模拟的过程的示例；
图17是示出基于各种指标来确定最优CRT设备和电极位置的处理过程的示例的流程图；
图18示出了根据第三实施例的CRT治疗和医疗费用偿付的过程的示例；
图19示出了根据第四实施例的处理过程的示例；
图20示出了根据第四实施例的CRT治疗和医疗费用偿付的过程的示例；
图21示出了根据第五实施例的CRT治疗和医疗费用偿付的过程的示例；
图22示出了根据第六实施例的CRT治疗和医疗费用偿付的过程的示例；
图23示出了根据第七实施例的CRT治疗和医疗费用偿付的过程的示例；
图24示出了根据第八实施例的CRT治疗和医疗费用偿付的过程的示例；
图25示出了根据第九实施例的CRT治疗和医疗费用偿付的过程的示例；
图26示出了根据第十实施例的CRT治疗和医疗费用偿付的过程的示例；
图27示出了根据第十一实施例的处理过程的示例；
图28示出了根据第十一实施例的CRT治疗和医疗费用偿付的过程的示例；
图29示出了根据第十二实施例的处理过程的示例；
图30示出了根据第十三实施例的处理过程的示例；
图31示出了根据第十四实施例的处理过程的示例；
图32示出了根据第十五实施例的处理过程的示例；以及
图33示出了根据第十六实施例的处理过程的示例。
具体实施方式
下面将参照附图描述若干实施例，其中在全文中相似的附图标记指代相似的元件。注意，可以以使得不产生冲突的方式来组合下面的实施例中的两个或更多个实施例以用于实现。
(a)第一实施例
接下来描述第一实施例。第一实施例涉及用于通过运行使用患者的器官的生物医学模型的生物医学模拟来调查手术的有效性并且仅在手术预期产生患者的改善时由保险来偿付手术的费用，以从而防止不必要的治疗并且减少医疗费用的技术。
图1示出了根据第一实施例的生物医学模拟系统的配置示例。图1的生物医学模拟系统包括多个计算机11至13。计算机11被安装在例如用于对相关患者执行手术的医疗机构处。计算机12被安装在例如用于运行生物医学模拟的模拟中心处。计算机13被安装在例如用于偿付医学治疗成本的健康保险公司处。
计算机11包括生物医学信息存储单元11a，模拟执行确定单元11b和手术决策单元11c。生物医学信息存储单元11a在其中存储有特定于患者的生物医学信息。生物医学信息包括患者的特定器官(例如心脏)的层析图像、关于心脏中的梗塞区域的信息、射血分数和ECG数据。模拟执行确定单元11b基于生物医学信息来确定是否运行生物医学模拟。例如，在要对患者执行的手术是CRT植入的情况下，模拟执行确定单元11b在射血分数下降至预定值(例如35％)以下、并且QRS持续时间超过预定值(例如130ms)的情况下确定运行生物医学模拟。是否执行生物医学模拟的确定结果被传送至计算机12。手术决策单元11c在计算机12的改善率确定单元12c确定从生物医学模拟获得的结果满足预定改善率时决定对患者执行手术。在决定之后，由医生来执行手术。
计算机12包括模型生成单元12a、模拟执行单元12b和改善率确定单元12c。模型生成单元12a在模拟执行确定单元11b确定执行生物医学模拟时基于生物医学信息来生成患者的特定器官的生物医学模型。例如，模型生成单元12a基于患者的心脏的层析图像来生成患者的心脏的立体模型(心脏模型)。此外，如果患者的心脏具有梗塞区域，则模型生成单元12a设置梗塞区域被表示在心脏模型中的相应位置处。
模拟执行单元12b基于由模型生成单元12a生成的生物医学模型来运行患者的特定器官的治疗后状况的生物医学模拟。在模拟涉及CRT设备的植入的情况下，例如将CRT设备的电极布置模式输入至模拟执行单元12b，以作为手术过程。然后，模拟执行单元12b运行模拟以基于输入的电极布置模式来预测布置CRT设备的电极的手术之后的术后心搏状况。心脏模拟输出例如ECG数据和射血分数值作为其模拟结果。
改善率确定单元12c确定由模拟执行单元12b运行的生物医学模拟的结果是否满足与症状相关联的预定改善率。例如，在已经执行心脏模拟的情况下，将模拟的结果与关于ECG QRS持续时间和射血分数的、从患者获取的实际测量值相比较。如果模拟预测了QRS持续时间和射血分数中的每个的在预定改善率或更高程度上的改善，则认为使用由模拟执行单元12b采用的手术过程的治疗是有效的。另一方面，如果QRS持续时间或射血分数都没有表现出预定改善率或更高程度上的改善，则认为使用所采用的手术过程的治疗不提供益处。模拟结果是否满足预定改善率的确定结果被传送至其它计算机11和13。
计算机13包括治疗费用偿付单元13a。当改善率确定单元12c确定生物医学模拟的结果满足预定改善率时，治疗费用偿付单元13a确定偿付针对患者的手术的治疗成本。另一方面，当改善率确定单元12c确定生物医学模拟结果不满足预定改善率时，治疗费用偿付单元13a确定不偿付针对手术的治疗成本。
根据上述生物医学模拟系统，对于其症状已经变得如此严重使得患者需要手术的患者，确定运行生物医学模拟，并且然后生物医学模拟预测手术之后患者的术后器官状况。患者的器官的生物医学模拟的执行使得能够准确地理解手术的有效性。因此，仅在预期手术产生患者的改善时执行手术，并且然后由保险来偿付手术的费用。因此，根据第一实施例，仅在生物医学模拟预测手术会产生治疗效果时执行手术，这防止无效手术的实践，并且从而实现医疗费用的减少。
注意，各个计算机11至13的每个部件可以被包括在与图1中示出的计算机不同的计算机中。例如，手术决策单元11c可以被包括在计算机13中。此外，改善率确定单元12c可以被包括在计算机13中。此外，手术决策单元11c和改善率确定单元12c可以被包括在计算机13中。
可以通过分开的第一决策单元和第二决策单元来配置手术决策单元11c。在这种情况下，当改善率确定单元12c确定生物医学模拟的结果满足预定改善率时，例如，第一决策单元输出针对保险偿付的申请。当改善率确定单元12c确定生物医学模拟的结果满足预定改善率并且然后第一决策单元已经输出保险偿付申请时，第二决策单元决定对患者执行手术。就此而言，第一决策单元和第二决策单元分别被包括在例如计算机11和计算机13中。由此，当寻求保险偿付时，可以检查保险偿付申请的有效性。
此外，手术决策单元11c可以在确认改善率确定单元12c确定生物医学模拟的结果满足预定改善率、并且治疗费用偿付单元13a确定偿付患者的医学治疗成本时决定对患者执行手术。这防止了其中在手术后发现没有可用的保险偿付的情况。
治疗费用偿付单元13a可以在已经对患者执行了手术时确定偿付针对手术的治疗成本和用于运行生物医学模拟的成本，并且可以在没有执行手术时确定偿付用于运行生物医学模拟的成本。由此，即使在没有对手术成本做出偿付时，患者仍然能够收到针对生物医学模拟服务的保险偿付。
还可以在生物医学模拟系统中设置设备特征数据存储单元以存储表示可能在患者的手术中使用的设备的特征的设备特征数据。在这种情况下，模拟执行单元12b可以参考设备特征数据并且运行其中反映每个设备的特征的模拟。由此，可以提高生物医学模拟的准确性，并且还可以适当地评估设备中的每个设备对患者的适用性。
此外，关于可用于患者的多个手术过程中的每个手术过程，模拟执行单元12b能够运行生物医学模拟以预测使用手术过程的手术之后患者的术后状况。在这种情况下，改善率确定单元12c在可用的手术过程中选择具有表现出最高改善率的模拟结果的一个手术过程，并且确定改善率是否满足预定改善率。这使得能够建议最适合的手术过程。例如，通过运行针对与CRT设备植入相关联的多个电子布置模式中的每个电极布置模式的模拟，最优电极布置模式被确定。然后，改善率确定单元12c将表示具有最高改善率的手术过程的信息传送至被安装在用于对患者执行手术的医疗机构处的信息处理器。这使得能够向将要对患者执行手术的医生呈现最优手术过程(例如，与CRT设备植入相关联的最优电极布置模式)。因此，缩短了手术时间，减轻了患者身体上的负担。
在存在多个电极布置模式、每个电极布置模式产生等于或大于预定改善率的改善率的情况下，改善率确定单元12c可以将这些电极布置模式传送至医疗机构处的信息处理器。例如，改善率确定单元12c传送具有最高的三个左右改善率的电极布置模式。这使得医生能够在多个电极布置模式中选择合适的布置模式并且将其应用到手术中。例如，在实际的手术中，患者的状况可能不允许将电极安装在由最优电极布置模式所表示的位置处。在这样的情况下，医生不是根据最优电极布置模式、而是根据具有第二或第三最高改善率的电极布置模式来安装电极。由此提高了CRT植入的确定性。
模拟执行单元12b可以针对可用于患者的手术的多个设备中的每个设备来运行生物医学模拟。在这种情况下，改善率确定单元12c在可用设备中选择具有表现出最高改善率的模拟结果的一个设备，并且确定改善率是否满足预定改善率。由此，可以适当地预测在使用最适于患者的设备时所产生的治疗有效性。此外，改善率确定单元12c可以将关于产生最高改善率的设备的信息传送至被安装在用于对患者执行手术的医疗机构处的计算机11。这使得能够向医生呈现对于患者最优的设备。此外，改善率确定单元12c可以将关于产生最高改善率的设备的信息也传送至被安装在用于偿付医学治疗成本的健康保险公司处的计算机13。这使得健康保险公司能够从成本效益的角度选择合适的设备。
注意，模拟执行确定单元11b、手术决策单元11c、模型生成单元12a。模拟执行单元12b、改善率确定单元12c和治疗费用偿付单元13a中的每个可以例如由单元与其集成在一起的计算机的处理器来实现。此外，可以例如由存储器或者计算机11的存储设备来实现生物医学信息存储单元11a。在图1中，连接各个部件的每条线表示通信路径的一部分，并且除了图1中示出的通信路径以外的通信路径也是可配置的。
(b)第二实施例
接下来描述第二实施例。第二实施例涉及向确定是否实施涉及CRT设备植入的治疗提供支持的计算机系统。涉及CRT设备植入的治疗在下文简称为“CRT治疗”。
图2示出了CRT设备植入的示例。在CRT中，CRT设备50是一种起搏器，其被植入到患有心脏病的患者40中。CRT设备50具有多个导线(lead)51至53，并且电极51a、52a和53a分别被设置在各个导线51至53上。导线51至53被插入到患者40的身体中，并且电极51a、52a和53a被布置在心脏41中的适当位置处。例如，导线51上的电极51a被植入到右心房中，导线52上的电极52a被植入到右心室中，以及导线53上的电极53a被布置在冠状窦中。CRT设备50以使心室和室间隔与患者40的心脏41的起搏同步收缩的方式向各个电极51a、52a和53a施加电势。
例如，患有心肌梗塞的患者40的梗塞区域不传导电流，并且因此不能够正常地传导脉冲传导系统的电信号。CRT设备50起到支持脉冲传导的作用。注意，CRT治疗之前的患者40的心脏41的脉冲传导的状况根据梗塞的位置而不同。因此，CRT设备50的电极51a、52a和53a的最优布置位置取决于患者40。CRT模拟能够识别适于患者40的CRT设备和电极布置模式。
根据第二实施例的计算机系统使用非侵入性CRT模拟来在CRT治疗之前预测CRT设备的植入之后的效果(ECG QRS持续时间和射血分数的改善)。然后，系统基于预测结果来选择由最适于患者40的制造商生产的CRT设备，并且呈现CRT设备的使用说明(如最优电极安装位置和电极使用配置(如双心室起搏和三心室起搏))。于是，根据第二实施例的系统使得能够在医生、患者和医疗保险公司之间共享基于预测结果的信息。
由此，预先识别出了最适合的制造商的CRT设备和最优电极安装位置，并且因此缩短了外科手术时间，这又减轻了对患者的负担。如果患者可能是不响应者，则选择不同的治疗方法以从而减少不必要的患者的身体上的和经济上的负担，并且还防止浪费医疗费用。此外，通过预测CRT植入之后的效果以及与患者共享关于预测结果的信息，可以降低医生在医疗事故诉讼中被起诉的风险。
图3示出了根据第二实施例的系统配置的示例。医院信息系统100被安装在医院31处。医院信息系统100是用于管理患者的医疗图表和生物医学信息的计算机系统。CRT模拟系统200被安装在用于提供CRT模拟服务的心脏模拟中心(HSC)32处。CRT模拟系统200是用于运行与CRT设备的植入相关联的心脏模拟(CRT模拟)的计算机系统。医疗费用偿付系统300被安装在医疗保险公司33处。医疗费用偿付系统300是用于执行向医院31偿付被保险覆盖的治疗所涉及的成本的过程的计算机系统。
图3的医院信息系统100、CRT模拟系统200和医疗费用偿付系统300协作以顺利地执行CRT模拟，并且使得能够在医生、患者和医疗保险公司33之间共享从模拟获得的结果。具体地，当在医院31处医生确定CRT治疗适用于患者时，将模拟申请从医院信息系统100发送至医疗费用偿付系统300(步骤S101)。在接收到来自医疗保险公司33的批准之后，医院信息系统100向CRT模拟系统200发送模拟请求，医院31在模拟请求上附有模拟所需的患者的数据(步骤S102)。患者数据包括待用于创建特定于患者的心脏模型的生物医学数据。
在HSC 32处，使用模拟器来预测患者中的CRT设备的最优电极安装位置，并且将从模拟获得的结果发送至医院信息系统(步骤S103)。模拟结果包括被预测产生患者的心脏功能的最大改善的电极位置的组合模式以及与所有组合模式中的每个组合模式相关联的心脏功能的数据(如ECG数据、射血分数、dP/dt和使心室运动可视化的数据)。就此而言，HSC32向医疗保险公司33做出针对模拟中产生的成本的索取。作为响应，医疗费用偿付系统300处理支付模拟成本的过程，并且然后向CRT模拟系统200通知过程结果(步骤S104)。支付过程涉及例如向HSC 32的银行账户发送覆盖模拟成本的钱款。
在医院31处，医院信息系统100基于由HSC 32发送的模拟结果来确定是否提供CRT治疗(步骤S105)。当医院信息系统100确定提供CRT治疗时，医生执行CRT植入。在CRT植入之后，医院31将针对医学治疗成本的偿付的请求发送至医疗保险公司33。作为响应，医疗费用偿付系统300处理支付CRT治疗成本的过程并且然后向医院信息系统100通知过程结果(步骤S106)。支付过程涉及例如向医院31的银行账户发送覆盖CRT治疗成本的钱款。以上述方式有效地使用CRT模拟的结果减少了没有效果的治疗。然而，注意，根据患者自己的意愿，他/她仍然能够通过支付全部费用来接受CRT治疗。
图4是示出了通过使用CRT模拟来确定是否提供CRT的过程的流程图。
【步骤S111】医院信息系统100基于患者的生物医学信息来确定以下条件是否被满足：射血分数在35％以下；以及QRS持续时间超过130毫秒。此外，可以增加医学治疗难以治愈的心衰作为另一条件。如果满足条件，则过程移至步骤S113。如果不满足条件，则过程移至步骤S112。
【步骤S112】因为患者的状态没有过于危重而接受CRT治疗，所以医院信息系统100确定不向患者提供CRT并且结束过程。
【步骤S113】医院信息系统100请求CRT模拟系统200运行CRT模拟。作为对请求的响应，CRT模拟系统200运行CRT模拟。
【步骤S114】基于从模拟获得的结果，CRT模拟系统200确定是否已经发现具有正面治疗效果的电极位置。如果已经发现这样的电极位置，则过程移至步骤S115。如果没有发现这样的治疗位置，则过程移至步骤S116。
【步骤S115】CRT模拟系统200将关于具有正面治疗效果的最优电极位置的信息发送至医院信息系统100。然后，医院31处的医生根据作为CRT模拟结果而获得的最优电极位置来执行CRT治疗。在这种情况下，CRT模拟结果还被发送至医疗费用偿付系统300，并且由医疗保险公司33来偿付医学治疗成本。随后过程结束。
【步骤S116】CRT模拟系统200向医院信息系统100通知没有检测到具有正面治疗效果的电极位置。在医院31处，医生询问患者或代表患者行事的人是否患者自付费用来寻求CRT治疗。如果患者决定自付费用来接受CRT治疗，则过程移至步骤S117。如果患者决定不自付费用来接受CRT治疗，则过程移至步骤S118。
【步骤S117】医生执行CRT治疗。在这种情况下，由患者来承担医学治疗成本。随后过程结束。
【步骤S118】医生确定不提供CRT治疗并且尝试不同的治疗方法。这通过不提供无效CRT治疗来实现医疗费用的减少，并且还减轻了对患者的负担。
以上述方式使用CRT模拟使得医生能够基于模拟结果来从手术开始时就将导电线插在最优电极安装位置处。因此，缩短了手术时间。此外，CRT模拟的使用帮助不仅使CRT治疗的有效性清楚，并且还使哪一方要承担医学治疗成本清楚，从而防止了非自愿地负担医学治疗成本的事件。例如，在不使用CRT模拟的情况下，不管患者是否是不响应者，只要满足准则中定义的条件(射血分数＜35％，以及QRS持续时间＞130毫秒)，就实施CRT治疗。另一方面，在使用CRT模拟的情况下，在CRT植入之前进行CRT模拟，并且仅在识别出被预测具有正面治疗效果的电极安装位置时，才执行CRT植入。这使得医疗保险公司33能够仅在客观地确定保证通过CRT植入产生一定程度的效果时才偿付CRT治疗的医学治疗成本，并且从而避免产生用于无效医学治疗的成本。
此外，CRT植入的有效性的高度准确的确定防止了患者在手术后被发现是不响应者的事件，这降低了医生在寻求损失赔偿等的诉讼中被起诉的风险。此外，即使决定不接受CRT治疗的患者之后病情变得危重，医生由于没有植入CRT而被起诉的风险也被降低，因为不接受CRT植入是由患者或代表患者行事的人做出的选择。此外，即使在患者上没有识别出被预测具有正面治疗效果的电极安装位置，患者仍然能够选择通过支付全部费用来接受CRT植入。在这种情况下，患者在清楚地知道他/他可能是不响应者之后接受CRT治疗，并且因此医生免于被患者起诉要求损失赔偿。
接下来更详细地描述第二实施例。首先，叙述了根据第二实施例的系统的硬件配置。图5示出了医院信息系统的硬件配置的示例。由处理器101来进行医院信息系统100的总体控制。存储器102和多个外围设备通过总线109来连接至处理器101。处理器101可以是多处理器。处理器101是例如中央处理单元(CPU)、微处理单元(MPU)，图形处理单元(GPU)，或数字信号处理器(DSP)。通过由处理器101执行程序来实现的功能的至少一部分可以被实现为电子电路，如专用集成电路(ASIC)和可编程逻辑器件(PLD)。存储器102被用作医院信息系统100的主存储器设备。存储器102暂时地存储待由处理器101执行的操作系统(OS)程序和应用程序的至少一部分。存储器102还在其中存储有待被处理器101使用以用于处理器101的处理的各种类型的数据。作为存储器102，可以使用易失性半导体存储设备如随机存取存储器(RAM)。
连接至总线109的外围设备包括硬盘驱动器(HDD)103、图形处理单元104、输入接口105、光盘驱动单元106、设备连接接口107和网络接口108。HDD 103向内置磁盘磁性地写入数据并且从内置磁盘磁性地读取数据，并且用作医院信息系统100的辅助存储设备。HDD103在其中存储有OS程序、应用程序和各种类型的数据。注意，非易失性半导体存储设备(固态驱动器或SSD)如闪速存储器可以代替HDD 103来用作辅助存储设备。
监视器21连接至图形处理单元104。图形处理单元104根据来自处理器101的指令来在监视器21的屏幕上显示图像。例如阴极射线管(CRT)显示器或液晶显示器可以用作监视器21。键盘22和鼠标23连接至输入接口105。输入接口105将发送自键盘22和鼠标23的信号传送至处理器101。注意，鼠标23仅是指向设备的示例，并且可以使用不同的指向设备如触摸面板、平板电脑、触摸垫和跟踪球来代替。光盘驱动单元106使用例如激光来读取记录在光盘24上的数据。光盘24是数据以通过光的反射来被读取的方式被记录在其上的便携式存储介质。光盘24的示例包括数字通用光盘(DVD)、DVD-RAM、致密盘只读存储器(CD-ROM)、CD可记录(CD-R)和CD可再写(CD-RW)。
设备连接接口107是用于将外围设备连接至医院信息系统100的通信接口。例如，存储器设备25和存储器读取器/写入器26可以连接至设备连接接口107。存储器设备25是具有用于与设备连接接口107通信的功能的存储介质。存储器读取器/写入器26是用于向存储卡27写入数据并且从存储卡27读取数据的设备，存储卡27是卡型存储介质。网络接口108连接至网络20。网络接口108通过网络20来向其他计算机和通信设备传送数据并且从其他计算机和通信设备接收数据。
上面描述的硬件配置实现了第二实施例的处理功能。注意，CRT模拟系统200和医疗费用偿付系统300可以被构建成具有与医院信息系统100的硬件配置相同的硬件配置。此外，第一实施例的计算机11至13中的每个也可以被构建成具有与图5的医院信息系统100的硬件配置相同的硬件配置。
医院信息系统100例如通过执行存储在计算机可读存储介质中的程序来实现其根据第二实施例的处理功能。描述待由医院信息系统100实现的处理内容的程序可以被存储在各种类型的存储介质中。例如，待由医院信息系统100执行的程序可以被存储在HDD 103中。处理器101将被存储在HDD 103中的程序的至少一部分加载到存储器102中，并且然后运行程序。此外，待由医院信息系统100执行的程序可以被存储在便携式存储介质如光盘24、存储器设备25和存储卡27中。存储在便携式存储介质中的程序在例如在处理器101的控制下被安装到HDD 103之后变得可执行。替选地，处理器101可以通过从便携式存储介质直接读取程序来运行程序。
根据图5的示例，单个计算机形成计算机系统；但是，计算机系统可以是由多个计算机组成的并行处理系统。例如，如果CRT模拟系统200被配置为具有多个计算机的并行处理系统，则可以在短时间内运行高度准确的CRT模拟。
接下来描述实现第二实施例的各个系统的功能。图6是示出了各个系统的功能的框图。医院信息系统100包括患者数据存储单元110、CRT模拟请求后台程序120和成本管理单元130。患者数据存储单元110在其中存储有通过图片存档和通信系统(PACS)获得的图像数据和针对每个患者的电子健康记录。PACS是用于管理医学图像的系统。例如，通过磁共振成像(MRI)和计算层析成像(CT)产生的图像通过PACS来被存储在患者数据存储单元110中。
CRT模拟请求后台程序120向CRT模拟系统200做出请求以用于进行CRT模拟，并且获取通过CRT模拟获得的模拟结果。此外，CRT模拟请求后台程序120基于模拟结果来确定是否是进行CRT。成本管理单元130管理医学治疗成本。例如，成本管理单元130从医疗费用偿付系统300接收保险覆盖是否可用的通知。然后，如果执行被保险覆盖的CRT植入，则成本管理单元130向医疗费用偿付系统300做出针对手术成本的要求。
CRT模拟系统200包括CRT设备列表存储单元210、CRT模拟管理单元220、患者数据存储单元230、心脏模型创建单元240、参数确定单元250、电极布置模式指定单元260和心脏模拟器270。CRT设备列表存储单元210是由不同的制造商生产和发售的多个CRT设备的列表。例如，CRT设备列表存储单元210包括CRT设备中的每个CRT设备的：制造商的名称、其型号的名称和设备的电极的数目。CRT模拟管理单元220管理CRT模拟的执行。例如，在接收到CRT模拟请求时，CRT模拟管理单元220将附在CRT模拟请求上的患者数据存储在患者数据存储单元230中。然后，CRT模拟管理单元220向心脏模型创建单元240、参数确定单元250和心脏模拟器270发出预定执行指令，以由此实现CRT模拟。
患者数据存储单元230在其中存储有特定于接受CRT模拟的每个患者的数据。被存储在患者数据存储单元230中的数据类似于被存储在医院信息系统100的患者数据存储单元110中的数据。心脏模型创建单元240基于被存储在患者数据存储单元230中的患者的数据来创建患者的心脏模型。心脏模型创建单元240使用例如MRI的层析图像来创建患者的心脏模型。参数确定单元250根据患者的心脏来确定在CRT模拟中使用的参数的值。例如，参数确定单元250以使得通过对没有CRT植入的患者的心脏模型的运动模拟获得的ECG数据与患者的实际ECG数据一致的方式来调节参数值。
电极布置模式指定单元260针对每个CRT设备设置由心脏模型创建单元240创建的患者的心脏模型上的多个电极布置模式。心脏模拟器270关于每个CRT设备来运行针对CRT设备被植入到患者中的情况的CRT模拟。基于目标CRT设备的产品信息来运行CRT模拟。如果已经针对单个CRT设备建议多个电极布置模式，则心脏模拟器270针对所有布置模式运行CRT模拟。重复执行的CRT模拟的结果被发送至CRT模拟管理单元220，并且然后被用于确定对于患者最优的CRT设备和电极位置。
医疗费用偿付系统300包括偿付处理单元310。偿付处理单元310基于CRT模拟结果来确定是否偿付CRT治疗中产生的费用。然后，当执行被保险覆盖的手术时，偿付处理单元310安排成本的支付。
在图6中，连接各个部件的每条线表示通信路径的一部分，并且除了图6中示出的通信路径以外的通信路径也是可配置的。此外，例如通过使计算机执行对应于图6中示出的每个部件的程序模块来实现部件的功能。
接下来详细描述根据第二实施例的CRT治疗和医疗费用偿付处理。图7示出了根据第二实施例的CRT治疗和医疗费用偿付的过程的示例。医院信息系统100的患者数据存储单元110存储待用于CRT模拟的患者数据。患者数据包括例如MRI或CT的层析图像数据、12导线ECG数据、心搏出量、左心室射血分数、血压的变化(dP/dt)和梗塞区域数据。心搏出量是随着心脏的每次收缩被射入主动脉的血液的量。心搏出量以例如毫升(ml)为单位测量。通过以下等式来计算心搏出量：
SV＝EDV–ESV；
其中，SV是心搏出量，EDV是舒张期末容积，以及ESV是收缩期末容积。血压的变化是左心室中血压的变化率。血压的变化是评估整个心脏的收缩功能并且通过心导管检查、多普勒超声心动描记术等测量的常用指数。梗塞区域数据表示患者的心脏的梗塞区域的位置。
医院信息系统100的CRT模拟请求后台程序120确定患者是否是CRT治疗的目标。例如，CRT模拟请求后台程序120从患者数据存储单元110获取患者的射血分数和ECG数据，并且根据ECG数据来计算QRS持续时间。然后，如果射血分数在预定值(例如，35％)以下，并且QRS持续时间超过预定值(例如，130毫秒)，则CRT模拟请求后台程序120确定患者是CRT治疗的目标。当已经将患者确定为CRT治疗目标时，CRT模拟请求后台程序120向CRT模拟系统200发送模拟请求(步骤S121)。患者的患者数据被附在模拟请求上。CRT模拟请求后台程序120可以以加密的形式发送带有患者数据的模拟请求。
在CRT模拟系统200中，CRT模拟管理单元220将接收到的患者数据存储在患者数据存储单元230中。随后，心脏模型创建单元240基于患者的心脏的层析图像来创建表示患者的心脏结构的形态学数据(如心脏形态学模型和胸廓(胸腔)模型)和模拟网格数据(如有限元网格模型)(步骤S122)。此外，心脏模型创建单元240对胸廓(胸腔)进行分割并且设置电导率(对于各个器官、脂肪、肌肉和骨骼)，并且还对心脏进行分割并且设置心脏的电导率。梗塞区域被设置成具有非兴奋性组织。此外，心脏模型创建单元240将表示心脏纤维方向的数据映射到心脏形态学模型上。这个模型创建处理是CRT模拟的准备步骤。
接下来，参数确定单元250基于心搏出量、血压和梗塞区域数据来确定在心脏动态系统模拟中使用的参数的值。(步骤S123)。针对动态系统模拟的参数确定处理的细节将在后面描述(参见图8)。此外，参数确定单元250基于ECG数据来确定在心脏电系统模拟中使用的参数的值(步骤S124)。这使得能够在CRT模拟系统中重建患者的心脏。针对电系统模拟的参数确定处理的细节将在后面描述(参见图10)。
电极布置模式指定单元260从CRT设备列表存储单元210获取CRT设备的列表。然后，电极布置模式指定单元260使用包括在所获取的CRT设备列表中的电极信息来为每个CRT设备指定多个电极布置模式(N个模式，其中N是大于等于1的整数)，多个电极布置模式中的每个电极布置模式具有针对被植入在右心室中和左心室的外侧上的各个电极的位置的不同的组合(步骤S125)。注意，电极布置模式中的每个电极布置模式表示对应的CRT设备的名称和用于CRT设备的植入的电极位置。心脏模拟器270针对每个CRT设备的每个电极布置模式运行心脏行为的模拟(步骤S126)。
通过运行模拟，估计的射血分数、估计的(dP/dt)max等被输出为针对每个电极布置模式的模拟结果。(dP/dt)max是收缩周期过程中变化的dP/dt(心室压力的时间微分)值中的最大值。CRT模拟管理单元220将被计算作为模拟结果的指标(例如，射血分数和(dP/dt)max)中的每个指标与用于计算的对应电极布置模式相关联。然后，CRT模拟管理单元220以降序来对通过重复的心脏模拟计算的指标(例如，射血分数和(dP/dt)max)中的每个指标的值进行排序(步骤S127)。
一旦完成针对所有CRT设备的所有电极布置模式的心脏模拟，CRT模拟管理单元220识别产生最高改善率的CRT设备及其电极布置模式。例如，CRT模拟管理单元220将CRT植入之前的患者的射血分数和(dP/dt)max与作为模拟结果获得的射血分数和(dP/dt)max相比较，并且针对每个电极布置模式来计算改善率。然后，CRT模拟管理单元220确定具有最高改善率的电极布置模式。注意，对应于具有最高改善率的电极布置模式的CRT设备和电极位置被识别为最优CRT设备和电极位置。
CRT模拟管理单元220确定最优CRT设备和电极位置的改善率是否大于或等于指定值(步骤S128)。例如，在预测射血分数改善10％或更多、或者预测(dP/dt)的最大值改善100帕斯卡每秒(Pa/s)或更多时，CRT模拟管理单元220确定改善率大于或等于指定值。如果改善率大于或等于指定值，则CRT模拟管理单元220将关于最优CRT设备和电极位置的信息发送至医院信息系统100(步骤S129)。在医院信息系统100中，成本管理单元130获取关于最优CRT设备和电极位置的信息。在获取到最优CRT设备和电极位置时，成本管理单元130向决策变量X插入(plug in)1(初始值是0)。
如果改善率在指定值以下，则CRT模拟管理单元220向医疗费用偿付系统300通知患者被预测没有经历改善(即，患者可能是不响应者)。在医疗费用偿付系统300中，偿付处理单元310接收来自CRT模拟系统200的通知，并且然后确定对于患者的CRT治疗不被保险覆盖，并且记录决定。随后，偿付处理单元310向医院信息系统100通知针对患者的CRT治疗的保险覆盖不可用(步骤S130)。
在接收到保险覆盖不可用的通知时，医院信息系统100的成本管理单元130向决策变量Y插入0(初始值是0)。然后，成本管理单元130计算通过X+Y获得的值(步骤S131)。当X+Y等于0时，成本管理单元130促使医生确定是否对患者执行CRT植入。例如，成本管理单元130在由医生使用的终端设备上显示指示患者可能是不响应者、并且针对CRT植入的保险覆盖将不可用的信息。医生参照决策(X+Y)的结果来与患者交谈，并且然后确定是否对患者执行CRT植入(步骤S132)。在这种情况下，如果确定执行CRT植入，则医疗保险公司33将不做出针对医疗费用的偿付。当X+Y等于1时，成本管理单元130在医生的终端设备上显示指示CRT治疗预期在患者中有效并且将被保险覆盖的信息。基于模拟结果，医生执行CRT植入(步骤S133)。在这种情况下，将由医疗保险公司33来偿付医学治疗成本。
在上述方式中，CRT模拟结果在医生、患者与医疗保险公司33之间共享，以由此防止不必要的CRT治疗。此外，在执行CRT治疗时，可以决定哪一方将以对于所有这相关的三方可接受的方式来承担成本。
接下来详细地描述用于心脏模拟的参数确定处理。参数确定处理分为图8和图9中示出的用于动态系统模拟的参数确定处理(“动态系统模拟参数确定处理”)以及图10和图11中示出的用于电系统模拟的参数确定处理(“电系统模拟参数确定处理”)。
首先描述动态系统模拟参数确定处理。如图8和图9所示，通过参数确定单元250和心脏模拟器270的协作来确定用于动态系统模拟的参数。图8是示出了确定动态系统模拟参数的处理过程的示例的流程图。根据流程图中的步骤编号来描述图8的处理。这里假设心脏模型创建单元240在开始动态系统模拟参数确定处理之前已经创建表示患者的心脏结构的形态学数据和模拟网格数据。
【步骤S141】参数确定单元250调整心脏模拟中使用的参数。待被调整的参数包括与心搏出量、血压、梗塞区域等相关的参数。例如，参数确定单元250以使得从有限元网格模型获得的心脏模拟结果接近实际从患者获取的生物医学数据(如ECG数据、血压、超声心动描记数据、MRI/CT数据和导管检查数据)的方式来调整参数。例如根据由医生提供的输入、参照之前的模拟结果来做出参数调整。
【步骤S142】心脏模拟器270使用已创建的有限元网格模型和经调整的参数来运行心脏模拟。例如，心脏模拟器270执行高度准确的心脏流体结构交互(FSI)模拟。注意，在FSI模拟中，心脏模拟器270可以使用机械心脏瓣膜或其中整流器电路仿真瓣膜的模拟心脏瓣膜。
【步骤S143】参数确定单元250执行检查术前状态的处理。例如，参数确定单元250显示根据代表此刻患者的活体心脏的有限元网格模型和关于此刻患者的生物医学状态的数据如心搏出量获得的模拟结果，以便对它们进行定量比较。如果模拟结果与患者当前的生物医学状态不同，则基于例如医生的指令来再次对参数进行调整。
【步骤S144】心脏模拟器270显示作为模拟结果而获得的心脏的行为。例如，心脏模拟器270基于模拟结果来以动画方式再现心脏的行为。
接下来详细描述步骤S142的心脏模拟。图9是示出心脏模拟的过程的示例的流程图。根据流程图中的步骤编号来描述图9的处理。
【步骤S151】心脏模拟器270以时间步长来更新值。例如，心脏模拟器270将预定时间增量Δt加到当前时间t上。注意，模拟开始时间是例如t＝0。
【步骤S152】心脏模拟器270执行通过任意拉格朗日-欧拉(ALE)方法和未确定乘数的拉格朗日方法的组合建立的FSI分析。
【步骤153】心脏模拟器270更新ALE网格。
【步骤S154】心脏模拟器270确定计算结果是否已经收敛。如果结果已经收敛，则过程移至步骤S155。如果结果没有收敛，则过程移至步骤S152并且再次执行计算。
【步骤S155】心脏模拟器270确定当前模拟时间t是否已经达到结束时间tend。结束时间tend是例如完成一次心搏的模拟所花费的时间。如果当前时间t已经达到结束时间tend，则心脏模拟结束。如果当前时间t还没有达到结束时间tend，则过程移至步骤S151。
接下来详细描述电系统模拟参数确定处理。图10是示出了确定电系统模拟参数的处理过程的示例的流程图。根据流程图中的步骤编号来描述图10的处理。假设在这里，在开始电系统模拟参数确定处理之前，心脏模型创建单元240已经执行胸廓(胸腔)的分割并且设置电导率(对于各个器官、脂肪、肌肉和骨骼)，并且也已经执行心脏的分割并且设置心脏的电导率。在这样的假设下设置患者的心脏中的梗塞区域的电导率：梗塞区域不是电可兴奋的，并且不传导电流。
【步骤S161】参数确定单元250调整参数x的值以用于调节QRS波。参数x表示例如早期激励起始位点(“早期激励部位”)。例如，在参数x的多个候选值xi(i＝1,2,3,...,imax)中，选择最优值。参数确定单元250调节心内膜表面上的早期激励部位的位置以从而调节QRS波。参数调整处理的细节将在后面描述(参见图11)。
【步骤S162】参数确定单元250确定从使用经调整的参数的心脏模拟获得的ECG上的波形f(t)与从患者获得的ECG上的测量到的波形freal(t)之间的QRS间隔的差异(或者与差异对应的值)是否在阈值以下。例如，参数确定单元250可以将f(t)与freal(t)之间的差异除以freal(t)，并且然后确定相除的结果是否落在阈值以下。这通过以下不等式来表达：
|(f(t)-freal(t))/freal(t)|<阈值。
如果上面的不等式对于QRS复合波中的每个时间t都成立，则过程移至步骤S163。如果存在不等式对于其不成立的时间t，则过程移至步骤S170。
【步骤S163】参数确定单元250将在步骤S161中选择的参数xi确定为用于调节QRS波的参数x的值。
【步骤S164】参数确定单元250调整参数yi的值以用于调节T波。参数yi是例如表示心室心肌中具有长动作电位持续时间(APD)的细胞和具有短APD的细胞的分布的值。APD是心室心肌收缩的时间量，并且对应于ECG波形的一部分，从QRS复合波开始到T波的结束(QT间隔)。例如，参数确定单元250通过调节具有不同的APD的三个细胞模型(心内膜内皮细胞、心外膜细胞和M(中层心肌(mid-myocardial))细胞)的分布来调节T波。M细胞是心肌细胞，被称为心肌细胞(cardiomyocytes)，其具有最长的APD并且分布在心室壁的中间层中。
【步骤S165】参数确定单元250确定从使用经调整的参数的心脏模拟获得的ECG上的波形f(t)与从患者获得的ECG上的测量到的波形freal(t)之间的T波的差异(或者与差异对应的值)是否总是保持在阈值以下。例如，参数确定单元250可以将f(t)与freal(t)之间的差异除以freal(t)，并且然后确定相除的结果是否保持在阈值以下。如果在T波期间T波的波形差异或相除的结果总是保持在阈值以下，则过程移至步骤S166。如果存在T波的波形差异或相除的结果超过或等于阈值的时间段，则过程移至步骤S170。
【步骤S166】参数确定单元250将在步骤S164中选择的参数yi确定为用于调节T波的参数的值。
【步骤S167】参数确定单元250调整参数zi的值以用于调节ECG波形的幅值。参数zi表示例如身体表面上的皮下脂肪的厚度。例如，参数确定单元250通过以下方式来调节幅值：如果ECG幅值总体上低，则将胸腔中的皮下脂肪层设置为薄；如果ECG幅值总体上高，则将胸腔中的皮下脂肪层设置为厚。
【步骤S168】参数确定单元250确定从使用经调整的参数的心脏模拟获得的ECG上的波形f(t)与从患者获得的ECG上的测量到的波形freal(t)之间的幅值的差异(或者与差异对应的值)是否总是保持在阈值以下。例如，参数确定单元250可以将f(t)与freal(t)之间的差异除以freal(t)，并且然后确定相除的结果是否保持在阈值以下。如果幅值的波形差异或相除的结果总是保持在阈值以下，则过程移至步骤S169。如果存在幅值的波形差异或相除的结果超过或等于阈值的时间段，则过程移至步骤S170。
【步骤S169】参数确定单元250将在步骤S167中选择的参数zi确定为用于调节幅值的参数的值。随后，电系统模拟参数确定处理结束。
【步骤S170】参数确定单元250确定没有做出适当的参数调整并且结束过程。在这种情况下，CRT模拟管理单元220向HSC 32的医院信息系统100通知不能够处理相关患者。
接下来描述参数调整处理(步骤S161、S164和S167)的细节。图11是示出参数调整的处理过程的示例的流程图。根据流程图中的步骤编号来描述图11的处理。
【步骤S181】参数确定单元250将变量i设置为初始值1。
【步骤S182】参数确定单元250确定第i个参数值xi。在所准备的候选参数值中，第i个值被确定为参数值xi。
【步骤S183】参数确定单元250请求心脏模拟器270计算对于参数值xi的电势。电势计算的细节将在后面描述(参见图12)。
【步骤S184】参数确定单元250从心脏模拟器270获取作为模拟的结果而获得的ECGfi(t)。
【步骤S185】参数确定单元250将与参数值xi关联的ECG fi(t)存储在存储器中。
【步骤S186】参数确定单元250确定变量i的值是否已经达到最大值imax。如果其没有达到最大值imax，则过程移至步骤S187。如果其已经达到最大值imax，则过程移至步骤S188。
【步骤S187】参数确定单元250使变量i加1，并且然后移至步骤S182。
【步骤S188】参数确定单元250在所有参数值xi中选择具有最高交叉相关性的参数值xi作为参数调整的结果。关于这一点，例如，交叉相关性RNCC(-1≤RNCC≤1)是通过以下等式、基于从使用参数xi的心脏模拟获得的ECG波形A(k,j)和从患者获得的测量到的ECG波形B(k,j)来计算的。
在上面的表达式中，A(k,j)是在从心脏模拟获得的12导线ECG中的第j个ECG中的第k(k＝1,2,3,...,N)模拟时钟时间处获得的模拟结果值(电势)。B(k,j)是在从患者获得的测量到的12导线ECG中的第j个ECG中的与第k模拟时钟时间对应的时间处获得的值。+1的交叉相关性表示完全正相关；0表示不相关；以及-1表示完全负相关。参数确定单元250选择具有最高交叉相关性的xi作为最佳参数值。随后，参数调整处理结束。
在上述方式中，选择了用于调节QRS波的最佳参数值xi。如果从使用所选择的参数xi的心脏模拟获得的ECG满足步骤S162中的条件，则将所选择的参数值xi确定为用于调节QRS波的参数的值。注意，图11示出了调整QRS波调节参数的示例，然而，可以以相同的方式来调整T波调节参数和幅值调节参数中的每个。
接下来详细描述电势计算处理。图12是示出计算电势的处理过程的示例的流程图。根据流程图中的步骤编号来描述图12的处理。
【步骤S191】参数确定单元250将心脏形态学模型和用于做出与测量到的ECG波形非常相似的模拟ECG波形的参数输入至心脏模拟器270。
【步骤S192】心脏模拟器270更新时间步长的值。例如，心脏模拟器270将预定时间增量Δt加到当前时间t上。注意，模拟开始时间是t＝0。
【步骤S193】心脏模拟器270根据细胞模型来计算膜电势Vm的时间演变。
【步骤S194】心脏模拟器270使用有限元方法来使双域(bi-domain)模型的兴奋传播等式离散化，以获得针对细胞外电势φe的线性方程组。
【步骤S195】心脏模拟器270对线性方程组进行求解。
【步骤S196】心脏模拟器270针对离散的下一个时间步长计算因变量即细胞间电势φi和细胞外电势φe的值。
【步骤S197】心脏模拟器270输出计算得到的因变量值。输出值被存储在存储器中。
【步骤S198】心脏模拟器270确定当前模拟时间t是否已经达到结束时间tend。结束时间tend是例如完成一次心搏的模拟所花费的时间。如果当前时间t已经达到结束时间tend，则电势计算处理结束。如果当前时间t还没有达到结束时间tend，则过程移至步骤S192。
通过图8至图12的处理，确定了用于在心脏模拟中准确地再现患者的ECG的参数。使用所确定的参数执行的心脏模拟被认为准确地复制CRT植入之前的患者的心脏的真实的行为。因此，可以通过运行使用所确定的参数对从CRT设备发送至患者的心脏的电信号进行建模的心脏模拟来准确地预测CRT植入之后患者的心脏行为。
接下来详细地描述针对每个CRT设备运行以便确定最优CRT设备和电极位置的心脏模拟。图13是示出针对每个CRT设备的心脏模拟的处理过程的示例的流程图。根据流程图中的步骤编号来描述图13的处理。
【步骤S201】心脏模拟器270选择要进行模拟的一个CRT设备。例如，心脏模拟器270从CRT设备列表存储单元210的CRT设备列表中选择一个CRT设备。
【步骤S202】心脏模拟器270确定所选择的CRT设备的类型。例如，如果输入到其中的电极布置模式包括被布置在心室中的两个电极，则心脏模拟器270确定设备类型为“双心室起搏”。类似地，在三个电极和四个电极被布置在心室中的情况下，设备类型分别被确定为“三心室起搏”和“四心室起搏”。当设备类型是双心室起搏时，过程移至步骤S203。当设备类型是三心室起搏时，过程移至步骤S204。当设备类型是四心室起搏时，过程移至步骤S205。
【步骤S203】心脏模拟器270运行被设计用于双心室起搏CRT设备的心脏模拟。这个处理的细节将在后面描述(参见图14)。随后，过程移至步骤S206。
【步骤S204】心脏模拟器270运行被设计用于三心室起搏CRT设备的心脏模拟。这个处理的细节将在后面描述(参见图15)。随后，过程移至步骤S206。
【步骤S205】心脏模拟器270运行被设计用于四心室起搏CRT设备的心脏模拟。这个处理的细节将在后面描述(参见图16)。随后，过程移至步骤S206。
【步骤S206】心脏模拟器270确定是否存在未被选择的CRT设备。如果存在未被选择的CRT设备，则过程移至步骤S201。如果不存在未被选择的CRT设备，则过程移至步骤S207。
【步骤S207】CRT模拟管理单元220基于从针对每个CRT设备的各个电极布置模式的模拟获得的结果来确定产生最高改善率的CRT设备及其电极位置。例如，CRT模拟管理单元220选择具有射血分数的最高改善率的电极布置模式作为第一候选模式，并且选择具有(dP/dt)max的最高改善率的电极布置模式作为第二候选模式。然后，CRT模拟管理单元220在具有射血分数的改善率的第一候选模式与具有(dP/dt)max的改善率的第二候选模式之间选择具有较高改善率的候选模式，并且将与所选择的候选模式相关联的CRT设备及其电极位置确定为最优CRT设备和电极位置。
假设，例如，在步骤S203至S205中，已经针对N个不同的电极布置模式执行了心脏模拟，并且从1开始对电极布置模式按顺序进行编号。第i个电极布置模式(i＝1,...,N)的射血分数和(dP/dt)max分别被表示为EFi和(dPi/dt)max。最大改善率IRMax被表达为：
IRMax＝{Max((EFi/EFe):i＝0,...,N),Max((dPi/dt)max/(dPe/dt)max:i＝0,...,N)
其中，EFe是治疗前的射血分数，以及(dPe/dt)max是治疗前的(dP/dt)max。与产生最大改善率的电极布置模式相关联的CRT设备和电极位置是最优CRT设备和电极位置。
在上述过程中，确定了最优CRT设备和电极位置。接下来描述图13的各个处理(步骤S203至205)的细节。假设在以下描述中，存在允许电极布置的右心室中的Na个位置(Na是大于或等于1的整数)(在下文中被称为“右心室电极布置允许位置”)和允许电极布置的左心室中的Nb个位置(Nb是大于或等于1的整数)(“左心室电极布置允许位置”)。
首先描述执行双心室起搏的CRT设备的心脏模拟。图14示出了被设计用于双心室起搏CRT设备的心脏模拟的过程的示例。根据流程图中的步骤编号来描述图14的处理。
【步骤S211】心脏模拟器270在Na个右心室电极布置允许位置中选择一个未被选择的位置(即，还没有被选择用于电极布置的位置)，并且将其确定为右心室中的第一电极布置位置。
此时，心脏模拟器270将所有左心室电极布置允许位置的状态重置为“未选择”，并且然后移至步骤S212。
【步骤S212】心脏模拟器270在Nb个左心室电极布置允许位置中选择一个未被选择的位置，并且将其确定为左心室中的第一电极布置位置。
【步骤S213】心脏模拟器270基于所确定的电极布置位置来运行执行双心室起搏的CRT设备的心脏模拟。
【步骤S214】心脏模拟器270基于模拟的结果来计算射血分数和(dP/dt)max，并且然后将与由分别在步骤S211和S212中确定的电极布置位置组成的电极布置模式相关联的射血分数和(dP/dt)max存储在存储器中。
【步骤S215】心脏模拟器270确定是否存在未被选择的左心室电极布置允许位置。如果存在未被选择的左心室电极布置允许位置，则过程移至步骤S212。如果不存在未被选择的左心室电极布置允许位置，则过程移至步骤S216。
【步骤S216】心脏模拟器270确定是否存在未被选择的右心室电极布置允许位置。如果存在未被选择的右心室电极布置允许位置，则过程移至步骤S211。如果不存在未被选择的右心室电极布置允许位置，则心脏模拟结束。
接下来描述执行三心室起搏的CRT设备的心脏模拟。图15示出了被设计用于三心室起搏的CRT设备的心脏模拟的过程的示例。根据流程图中的步骤编号来描述图15的处理。
【步骤S221】心脏模拟器270在Na个右心室电极布置允许位置中选择一个未被选择的位置，并且将其确定为右心室中的第一电极布置位置。
此时，心脏模拟器270将在Nb个左心室电极布置允许位置中选择的所有两个位置的组合(“电极位置对”)的状态重置为“未选择”，并且然后移至步骤S222。
【步骤S222】心脏模拟器270选择左心室中的一个未被选择的电极位置对，并且将包括在所选择的电极位置对中的电极布置允许位置中的一个确定为左心室中的第一电极布置位置。
【步骤S223】心脏模拟器270将在步骤222中选择的电极位置对中的另一个电极布置允许位置确定为左心室中的第二电极布置位置。
【步骤S224】心脏模拟器270基于所确定的电极布置位置来运行执行三心室起搏CRT设备的心脏模拟。
【步骤S225】心脏模拟器270基于模拟的结果来计算射血分数和(dP/dt)max，并且然后将与由分别在步骤S221至S223中确定的电极布置位置组成的电极布置模式相关联的射血分数和(dP/dt)max存储在存储器中。
【步骤S226】心脏模拟器270确定是否存在未被选择的左心室电极位置对。如果存在未被选择的左心室电极位置对，则过程移至步骤S222。如果不存在未被选择的左心室电极位置对，则过程移至步骤S227。
【步骤S227】心脏模拟器270确定是否存在未被选择的右心室电极布置允许位置。如果存在未被选择的右心室电极布置允许位置，则过程移至步骤S221。如果不存在未被选择的右心室电极布置允许位置，则心脏模拟结束。
接下来描述执行四心室起搏的CRT设备的心脏模拟。图16示出了被设计用于四心室起搏CRT设备的心脏模拟的过程的示例。根据流程图中的步骤编号来描述图16的处理。
【步骤S231】心脏模拟器270从在Na个右心室电极布置允许位置中选择的所有两个位置的组合(“电极位置对”)中选择一个未被选择的电极位置对。然后，心脏模拟器270将包括在所选择的电极位置对中的电极布置允许位置中的一个确定为右心室中的第一电极布置位置。
【步骤S232】心脏模拟器270将在步骤231中选择的电极位置对中的另一个电极布置允许位置确定为右心室中的第二电极布置位置。
此时，心脏模拟器270将所有左心室电极位置对的状态重置为“未选择”，并且然后移至步骤S233。
【步骤S233】心脏模拟器270选择左心室中的一个未被选择的电极位置对，并且将包括在所选择的电极位置对中的电极布置允许位置中的一个确定为左心室中的第一电极布置位置。
【步骤S234】心脏模拟器270将在步骤233中选择的电极位置对中的另一个电极布置允许位置确定为左心室中的第二电极布置位置。
【步骤S235】心脏模拟器270基于所确定的电极布置位置来运行执行四心室起搏的CRT设备的心脏模拟。
【步骤S236】心脏模拟器270基于模拟的结果来计算射血分数和(dP/dt)max，并且然后将与由分别在步骤S231至S234中确定的电极布置位置组成的电极布置模式相关联的射血分数和(dP/dt)max存储在存储器中。
【步骤S237】心脏模拟器270确定是否存在未被选择的左心室电极位置对。如果存在未被选择的左心室电极位置对，则过程移至步骤S233。如果不存在未被选择的左心室电极位置对，则过程移至步骤S238。
【步骤S238】心脏模拟器270确定是否存在未被选择的右心室电极位置对。如果存在未被选择的右心室电极位置对，则过程移至步骤S231。如果不存在未被选择的右心室电极位置对，则心脏模拟结束。
通过图13至16的处理，关于每个CRT设备针对各种电极布置模式运行了心脏模拟。然后，将与产生最大改善率的电极布置模式相关联的CRT设备和电极位置确定为最优CRT设备和电极位置。
关于用于比较改善率的指标，可以使用以下指标中的至少一个指标：左心室收缩末期容积(LVESV)、心搏出量(SV)、二尖瓣回流(MR)的改善、脉压、室间(室间(inter-V))不同时性的改善和射血时间(ET)。LVESV是最大心脏收缩的时刻左心室中血液的体积。
二尖瓣回流(MR)是在每次心室收缩时血液向回(回流)通过二尖瓣瓣膜的渗漏。当出现二尖瓣回流时，随着左心室将血液泵入主动脉，一些血液向回渗漏进入左心房，从而增加左心房中的血容量和内部压力。左心房中升高的血压使从肺部延伸向心脏的肺静脉中的血压升高，并且使左心房扩张以容纳从心室向回渗漏的额外的血液。极度扩张的心房通常以不规律的模式快速搏动(称为房颤的病症)，这降低了心脏泵送效率，因为纤维性颤动的心房仅仅是在颤动而不是在泵送。因此，血液不正常流动通过心房，并且在腔室内可能会形成血液凝块。如果凝块松脱并且变成栓，则其被泵出心脏并且可能阻塞动脉，这可能会导致中风或对另一器官造成损伤。严重的回流可能导致心衰，其中心房中升高的压力造成肺中的积液(淤血(congestion))，或者其中从心室到身体的减少的血液的向前流动使得器官不能得到正常量的血液。左心室可能逐渐膨胀和弱化，从而进一步使心衰恶化。因此，如果患有二尖瓣回流的患者在症状上受益于CRT植入的治疗，则有充分理由来执行CRT植入。因此，二尖瓣回流的改善被认为是确定CRT治疗的有效性的重要指标之一。
脉压(pulse pressure)是收缩压与舒张压之间的差异。脉压随着心搏出量的升高而升高，并且随着动脉系统的容积的增加而降低。此外，脉压还由于动脉硬化而升高。脉压与收缩压和舒张压一起用于高血压等的诊断，并且用作确定心脏状况的有效指标。室间不同时性是表示左心室和右心室的不同步搏动的信息。射血时间是从心室射出血液所花费的时间。
所有各种指标的组合可以用于确定最优CRT设备和电极位置。图17是示出基于各种指标来确定最优CRT设备和电极位置的处理过程的示例的流程图。假设在这个示例中，已经针对N个不同的电极布置模式执行了心脏模拟。第i个电极布置模式(i＝1,2,3,...,N)的射血分数和(dP/dt)max分别被表示为EFi和(dPi/dt)max。此外，各个指数的改善率被表示如下：
●第i个电极布置模式的射血分数的改善率：αi,EF＝EFi/EFe
(其中，EFi是第i个电极布置模式的射血分数，以及EFe是治疗前的射血分数)；
●第i个电极布置模式的(dP/dt)max(血压变化率)的改善率：αi,p＝
(dPi/dt)max/(dPe/dt)max
(其中，(dPi/dt)max是第i个电极布置模式的(dP/dt)max，以及(dPe/dt)max是治疗前的(dP/dt)max)；
●第i个电极布置模式的LVESV的改善率：αi,LVESV＝LVESVi/LVESVe(其中LVESVi是第i个电极布置模式的LVESV，以及LVESVe是治疗前的LVESV)；
●第i个电极布置模式的心搏出量的改善率：αi,SV＝SVi/SVe
(其中SVi是第i个电极布置模式的心搏出量，以及SVe是治疗前的心搏出量)；
●第i个电极布置模式的二尖瓣回流的改善率：αi,MR＝MRi/MRe
(其中MRi是第i个电极布置模式的二尖瓣回流，以及MRe是治疗前的二尖瓣回流)；
●第i个电极布置模式的脉压的改善率：αi,PP＝PPi/PPe
(其中PPi是第i个电极布置模式的脉压，以及PPe是治疗前的脉压)；
●第i个电极布置模式的室间不同时性的改善率：αi,IVA＝IVAi/IVAe
(其中IVAi是第i个电极布置模式的室间不同时性，以及IVAe是治疗前的室间不同时性)；以及
·第i个电极布置模式的射血时间的改善率：αi,ET＝ETi/ETe
(其中ETi是第i个电极布置模式的射血时间，以及ETe是治疗前的射血时间)。
接下来根据流程图中的步骤编号来描述图17的处理。
【步骤S241】CRT模拟管理单元220设置变量i，以用于将电极布置模式指定为初始值1。
【步骤S242】CRT模拟管理单元220计算第i个电极布置模式的射血分数的改善率，“αi,EF＝EFi/EFe”。
【步骤S243】CRT模拟管理单元220计算第i个电极布置模式的(dP/dt)max(血压变化率)的改善率，“αi,p＝(dPi/dt)max/(dPe/dt)max”。
【步骤S244】CRT模拟管理单元220计算第i个电极布置模式的LVESV的改善率，“αi,LVESV＝LVESVi/LVESVe”。
【步骤S245】CRT模拟管理单元220计算第i个电极布置模式的心搏出量的改善率，“αi,SV＝SVi/SVe”。
【步骤S246】CRT模拟管理单元220计算第i个电极布置模式的二尖瓣回流的改善率，“αi,MR＝MRi/MRe”。
【步骤S247】CRT模拟管理单元220计算第i个电极布置模式的脉压的改善率，“αi,PP＝PPi/PPe”。
【步骤S248】CRT模拟管理单元220计算第i个电极布置模式的室间不同时性的改善率，“αi,IVA＝IVAi/IVAe”。
【步骤S249】CRT模拟管理单元220计算第i个电极布置模式的射血时间的改善率，“αi,ET＝ETi/ETe”。
【步骤S250】CRT模拟管理单元220使用第i个电极布置模式的各个指标的改善率的优化评估函数F(αi,EF,αi,p,αi,LVESV,αi,SV,αi,MR,αi,PP,αi,IVA,αi,ET)来计算第i个电极布置模式的优化评估值Fi。优化评估函数是例如下面的加权平均函数：
Fi＝βEFαi,EF+βpαi,p+βLVESVαi,LVESV+βSVαi,SV+βMRαi,MR+βPPαi,PP+βIVAαi,IVA+βETαi,ET
其中，βEF、βp、βLVESV、βSV、βMR、βPP、βIVA和βET是针对第i个电极布置模式的各个指标的加权值。加权值中的每个加权值是大于0但小于1的数字。从医疗的角度预先设置这些加权值。
【步骤S251】CRT模拟管理单元220确定变量i是否是N或更大值。如果变量i是N或更大值，则过程移至步骤S253。如果变量i在N以下，则过程移至步骤S252。
【步骤S252】CRT模拟管理单元220将变量i的值加1，并且然后移至步骤S242。
【步骤S253】CRT模拟管理单元220将与产生Fi(i＝1,...,N)的最大值的电极布置模式相关联的CRT设备和电极位置确定为最优CRT设备和电极位置。
在上述方式中，使用预先定义的优化评估函数来针对每个电极布置模式计算优化评估值。然后，将与产生最大优化评估值的电极布置模式相关联的CRT设备和电极位置确定为最优CRT设备和电极位置。
根据第二实施例，CRT模拟使得能够准确地确定对于患者的最优CRT设备和用于植入CRT设备的最优电极位置。然后，如果通过最优CRT设备和电极位置获得的改善率大于或等于指定值，则确定CRT植入对患者是有效的。关于CRT植入的有效性的确定结果在医生、患者和医疗保险公司33之间共享。因此，基于准确的信息来做出关于是否提供CRT植入以及哪一方要承担医学治疗成本的适当的决定。这使得能够防止不必要的CRT植入并且降低医疗费用。
(c)第三实施例
接下来描述第三实施例。第三实施例涉及在决定是否做出保险偿付方面赋予医生的确定以首要重要性。第三实施例与第二实施例的不同之处在于以下点。根据第二实施例，医疗保险公司33基于模拟结果来确定是否偿付医学治疗成本。但是，在模拟结果中加入医生的意见可以给出关于CRT治疗的有效性的更准确的判断。鉴于此，根据第三实施例，仅将医生的确定结果发送至医疗保险公司33。
图18示出了根据第三实施例的CRT治疗和医疗费用偿付的过程的示例。根据第三实施例，如果确定由CRT模拟系统200预测的改善率落在指定值以下(步骤S128)，则将确定结果发送至医院信息系统100。在医院信息系统100中，CRT模拟请求后台程序120基于改善率在指定值以下的确定结果来为变量Y插入0。随后，基于模拟结果来确定CRT治疗的有效性(步骤S131)，这与第二实施例中的情况相同。当预期CRT治疗不提供任何益处时，由医生来确定是否执行CRT治疗(步骤S132-1)。此时，成本管理单元130根据医生关于CRT治疗的确定将处理分成两支(步骤S132a)。如果医生确定执行CRT治疗，则然后医生基于模拟结果来执行CRT治疗(例如，使用最优CRT设备和电极位置)(步骤S133)。随后，成本管理单元130向医疗费用偿付系统300发送针对医学治疗成本和模拟成本的偿付的请求。另一方面，如果医生确定不进行CRT治疗，则成本管理单元130向医疗费用偿付系统300发送针对模拟成本的偿付的请求。医疗费用偿付系统300的偿付处理单元310根据来自医院信息系统100的请求来执行偿付医学治疗成本和/或模拟成本的处理(步骤S134)。
因此，根据第三实施例，模拟结果被发送至医院信息系统100，并且然后由医生来做出关于是否执行CRT治疗的最终决定。医疗保险公司33信任来自医生的反馈，并且做出保险偿付。
(d)第四实施例
接下来描述第四实施例。根据第四实施例，医疗保险公司33基于CRT模拟结果来确定CRT治疗的实施是否提供益处。图19示出了根据第四实施例的处理过程的示例。将模拟请求连同患者数据如医疗图像和ECG数据从医院31的医院信息系统100发送至HSC 32的CRT模拟系统200(步骤S301)。然后，由CRT模拟系统200获得的心脏模拟结果不是被发送至医院信息系统100，而是被发送至医疗保险公司33的医疗费用偿付系统300(步骤S302)。随后，由医疗费用偿付系统300来确定是否做出保险偿付，并且确定结果被发送至医院信息系统100(步骤S303)。
第四实施例与第二实施例的不同之处在于以下点。根据第二实施例，医院信息系统100的CRT模拟请求后台程序120基于模拟结果来决定CRT治疗的实施是否提供益处；然而，根据第四实施例，由医疗费用偿付系统300来做出这个决定。
图20示出了根据第四实施例的CRT治疗和医疗费用偿付的过程的示例。根据第四实施例，如果确定由CRT模拟系统200预测的改善率大于或等于指定值(步骤S128)，则将确定结果发送至医疗费用偿付系统300。在接收到来自CRT模拟系统200的确定结果时，医疗费用偿付系统300的偿付处理单元310为变量X插入1。另一方面，如果偿付处理单元310接收到来自CRT模拟系统200的表示改善率落在指定值以下的确定结果，则偿付处理单元310输出针对CRT治疗的保险覆盖不可用的通知(步骤S130)，并且为变量Y插入0。然后，偿付处理单元310计算通过X+Y获得的值，以从而确定CRT治疗的实施是否提供益处(步骤S131a)。如果X+Y＝1，则偿付处理单元310通知医院信息系统100针对CRT治疗的保险偿付将被做出。如果X+Y＝0，则偿付处理单元310通知医院信息系统100对CRT治疗没有保险偿付可用。
因此，根据第四实施例，医疗保险公司33直接从CRT模拟系统200接收模拟结果和确定结果，并且然后将这些结果呈至医院31。这使得医疗保险公司33能够在执行CRT治疗之前确定是否做出针对CRT治疗的保险偿付，这便利了快速偿付过程。
(e)第五实施例
接下来描述第五实施例。根据第五实施例，医疗保险公司33做出关于是否执行CRT治疗的初步决定；然而，关于是否执行CRT治疗的最终决定是由医生做出的，并且医疗保险公司33尊重医生的决定并且偿付医学治疗成本。
第五实施例与第三实施例的不同之处在于以下点。根据第三实施例，医院信息系统100的CRT模拟请求后台程序120基于模拟结果来决定CRT治疗的实施是否提供益处；然而，根据第五实施例，这个决定是由医疗费用偿付系统300做出的。
图21示出了根据第五实施例的CRT治疗和医疗费用偿付的过程的示例。根据第五实施例，如果确定由CRT模拟系统200预测的改善率大于或等于指定值(步骤S128)，则确定结果被发送至医疗费用偿付系统300。在接收到来自CRT模拟系统200的确定结果时，医疗费用偿付系统300的偿付处理单元310为变量X插入1。另一方面，如果偿付处理单元310接收到来自CRT模拟系统200的表示改善率落在指定值以下的确定结果，则偿付处理单元310输出针对CRT治疗的保险覆盖不可用的通知(步骤S130)，并且为变量Y插入0。然后，偿付处理单元310计算通过X+Y获得的值，以从而确定CRT治疗的实施是否提供益处(步骤S131a)。如果X+Y＝1，则偿付处理单元310通知医院信息系统100针对CRT治疗的保险偿付将被做出。如果X+Y＝0，则偿付处理单元310通知医院信息系统100没有保险偿付将对CRT治疗可用。但是，根据第五实施例，即使偿付处理单元310已经初步地确定不对CRT治疗进行偿付，如果医生决定执行CRT治疗，则偿付处理单元310仍然执行保险偿付处理(步骤S134)。
因此，根据第五实施例，即使从CRT模拟获得的改善率落在预定值以下，如果医院31处的医生确定实施CRT治疗是适当的，则医疗保险公司33信任医生的确定。当接收到来自医院信息系统100的针对医学治疗成本和模拟成本的偿付的请求时，医疗费用偿付系统300针对这些成本做出保险偿付。
(f)第六实施例
接下来描述第六实施例。根据第六实施例，医疗保险公司33进行CRT治疗的有效性的分析，如改善率的计算。第六实施例与第二实施例的不同之处在于以下点。根据第二实施例，CRT模拟系统200确定最优CRT设备和电极位置，以及使用最优CRT设备和电极位置实现的改善率是否大于或等于指定值。然而，根据第六实施例，由医疗费用偿付系统300来执行这些处理。
图22示出了根据第六实施例的CRT治疗和医疗费用偿付的过程的示例。根据第六实施例，CRT模拟系统200的CRT模拟管理单元220将心脏模拟的结果(如射血分数和(dP/dt)max)发送至医疗费用偿付系统300。在医疗费用偿付系统300中，偿付处理单元310参照CRT模拟系统200的患者数据存储单元230并且确定最优CRT设备和电极位置的改善率是否大于或等于指定值(步骤S128a)。然后，如果改善率大于或等于指定值，则偿付处理单元310将关于最优CRT设备和电极位置的信息发送至医院信息系统100(步骤S129a)。
因此，根据第六实施例，关于改善率是否大于或等于指定值的确定不是由CRT模拟系统200做出的，而是由医疗费用偿付系统300做出的。这使得医疗保险公司33能够分析提供保险偿付的指示的改善率，并且从统计学的角度来理解，为实现保险偿付，针对改善率的指定值需要被设置为多高或者多低。因此，对于健康保险公司33而言，确定合理的保险费以及设置针对改善率的适当的指定值变得更容易。
(g)第七实施例
接下来描述第七实施例。根据第七实施例，尽管医疗保险公司33进行CRT治疗的有效性的分析，如改善率的计算(这使得医生能够做出关于是否进行CRT治疗的最终决定)，并且信任来自医生的反馈并且做出保险偿付。第七实施例与第六实施例的不同之处在于以下点。第七实施例包括根据医生关于是否执行CRT治疗的确定来偿付医学治疗成本的功能。
图23示出了根据第七实施例的CRT治疗和医疗费用偿付的过程的示例。根据第七实施例，成本管理单元130根据医生关于是否执行CRT治疗的确定将处理分成两支(步骤S132a)。如果确定执行CRT治疗，则然后医生基于模拟结果来执行CRT治疗(步骤S133)。随后，成本管理单元130向医疗费用偿付系统300发送针对医学治疗成本和模拟成本的偿付的请求。另一方面，如果医生确定不执行CRT治疗，则成本管理单元130向医疗费用偿付系统300发送针对模拟成本的偿付的请求。医疗费用偿付系统300的偿付处理单元310根据来自医院信息系统100的请求来执行偿付医学治疗成本和/或模拟成本的处理(步骤S134)。
因此，根据第七实施例，模拟结果被发送至医疗费用偿付系统300，在医疗费用偿付系统300处，确定CRT治疗的有效性。即使确定CRT治疗不提供益处，如果医生决定执行CRT治疗，则医疗费用偿付系统300仍然执行针对CRT治疗的成本的偿付处理。因此，即使CRT模拟结果不能满足由医疗保险公司33预先设置的改善率，如果根据医生的决定执行了CRT治疗，则仍然由医疗保险公司33来偿付医学治疗成本。
(h)第八实施例
接下来描述第八实施例。根据第八实施例，医疗保险公司33不仅做出通过CRT治疗实现的改善率是否被预测满足指定值的确定，还做出关于是否执行CRT治疗的确定。第八实施例与第六实施例的不同之处在于以下点。根据第六实施例，医疗费用偿付系统300确定改善率是否大于或等于指定值；然而，医院信息系统100确定是否执行CRT治疗。根据第八实施例，医疗费用偿付系统300也做出关于是否执行CRT治疗的确定。
图24示出了根据第八实施例的CRT治疗和医疗费用偿付的过程的示例。根据第八实施例，如果确定改善率大于或等于指定值(步骤S128a)，则偿付处理单元310通知医院信息系统100最优CRT设备和电极位置，并且还为变量X插入1。另一方面，如果确定改善率落在指定值以下，则偿付处理单元310为变量Y插入0。然后，偿付处理单元310计算通过X+Y获得的值，以从而确定CRT治疗的实施是否提供益处(步骤S131a)。如果X+Y＝1，则偿付处理单元310通知医院信息系统100针对CRT治疗的保险偿付将被做出。如果X+Y＝0，则偿付处理单元310通知医院信息系统100没有保险偿付将对CRT治疗可用。
第八实施例使得医疗保险公司33不仅能够预先设置提供保险偿付的指示的改善率，还能够确定是否做出保险偿付，并且告知医院31CRT治疗的实施是否提供益处。这便利了快速偿付过程。
(i)第九实施例
接下来描述第九实施例。根据第九实施例，由医院31和医疗保险公司33二者来确定CRT治疗有效性(即，CRT治疗的实施是否提供益处)(双重检查)。第九实施例与第二实施例的不同之处在于以下点。根据第二实施例，医疗保险公司33直接接受来自医院31的针对医学治疗成本的偿付的申请。然而，根据第九实施例，医疗费用偿付系统300检查保险偿付申请的有效性。
图25示出了根据第九实施例的CRT治疗和医疗费用偿付的过程的示例。根据第九实施例，如果确定由CRT模拟系统200预测的改善率大于或等于指定值(步骤S128)，则确定结果被发送至医院信息系统100和医疗费用偿付系统300二者。当确定CRT治疗的实施提供益处时，医院信息系统100的成本管理单元130将表示保险偿付申请的值1发送至医疗费用偿付系统300(步骤S135)。
另一方面，在从CRT模拟系统200接收到表示改善率大于或等于指定值的确定结果时，医疗费用偿付系统300的偿付处理单元310为变量X插入1。如果从CRT模拟系统200接收到表示改善率落在指定值以下的确定结果，则偿付处理单元310为变量Y插入0。然后，偿付处理单元310计算通过X+Y获得的值，以从而确定CRT治疗的实施是否提供益处(步骤S131a)。如果X+Y＝1，则偿付处理单元310输出值1，其表示针对CRT治疗的保险偿付将被做出。如果X+Y＝0，则偿付处理单元310输出值0，其表示没有保险偿付将对CRT治疗可用。
然后，当由医院信息系统100发出保险偿付申请时，偿付处理单元310确定是否已经输出表示针对CRT治疗的保险偿付将被做出的值1(步骤S136)。例如，偿付处理单元310计算表示保险偿付申请的值“1”与由偿付处理单元310获得的X+Y的结果之间的异或，并且然后确定所计算的结果是否是0(两个输入都是“1”)或1(输入中的一个输入是“0”)。在已经输出表示针对CRT治疗的保险偿付将要被做出的值1(异或的结果是“0”)的情况下，偿付处理单元310通知医院信息系统100针对CRT治疗的保险偿付将要被做出。响应于该通知，在医院31处执行CRT治疗(步骤S133)。
如果还没有输出表示针对CRT治疗的保险偿付将要被做出的值1(异或的结果是“1”)，则偿付处理单元310执行针对保险偿付申请的错误处理(步骤S137)。作为错误处理，例如，偿付处理单元310在管理员的终端上显示表示错误的保险偿付申请已经被提交的消息。
因此，根据第九实施例，由医院31和医疗保险公司33二者来确定CRT治疗的实施是否提供益处，并且当医院31和医疗保险公司33的决定结果彼此一致时，医疗保险公司33对医学治疗成本进行偿付。由此，即使医院31单方面尝试对具有低改善率的患者实施CRT治疗，并且寻求保险偿付，医疗保险公司33仍然能够运行其自己的检查。
(j)第十实施例
接下来描述第十实施例。根据第十实施例，由医院31彼此对照地核查由医院31和医疗保险公司33分别获得的CRT治疗对患者的有效性的确定结果。第十实施例与第九实施例的不同之处在于以下几点。图26示出了根据第十实施例的CRT治疗和医疗费用偿付的过程的示例。根据第十实施例，由医疗费用偿付系统300的偿付处理单元310获得的CRT治疗有效性的确定结果(步骤S131a)被发送至医院信息系统100。医院信息系统100的成本管理单元130对照从医疗费用偿付系统300发送的确定结果来核查其自己的CRT治疗有效性的确定结果(步骤S131)。
例如，成本管理单元310确定两个确定结果是否都是“1”(CRT治疗将是有效的)(步骤S136a)。如果两个确定结果都是“1”，则执行CRT治疗，同时医学治疗成本被保险覆盖(步骤S133)。如果决定结果中的一个决定结果是“0”，则执行错误处理(步骤S137a)。
此外，成本管理单元130确定两个决定结果是否都是“0”(CRT治疗将不产生益处)(步骤S136b)。如果两个决定结果都是“0”，则医生决定是否执行CRT治疗(步骤S132)。如果决定结果中的一个决定结果是“1”，则执行错误处理(步骤S137b)。
因此，根据第十实施例，由医院31和医疗保险公司33分别确定CRT治疗对患者的有效性，并且然后医院31对照彼此来核查决定结果。这使得医院31能够防止医疗保险公司33错误地将具有有利的改善率的患者作为不响应者来对待的事件。
(k)第十一实施例
接下来描述第十一实施例。根据第十一实施例，在具有关于由多个制造商生产的CRT设备的说明书的信息的HSC 32处运行CRT模拟，从而提高确定对于患者最优的CRT设备的准确率。图27示出了根据第十一实施例的处理过程的示例。第十一实施例涉及使用由多个CRT设备制造商34a、34b和34c生产的CRT设备中的一个CRT设备的CRT治疗。CRT设备制造商34a具有制造商信息系统400。制造商信息系统400向CRT模拟系统200提供其CRT设备的说明书(步骤S311)。其余CRT设备制造商34b和34c中的每个也具有制造商信息系统，其向CRT模拟系统200提供其CRT设备的说明书(步骤S312和S313)。
随后，医院信息系统100将模拟请求连同患者数据发送至CRT模拟系统200(步骤S314)。CRT模拟系统200运行CRT模拟。然后，CRT模拟系统200通知医院信息系统100最优CRT设备和电极安装方法(步骤S315)。CRT模拟系统200也向医疗费用偿付系统300通知模拟结果(步骤S316)。随后，医疗费用偿付系统300通知医院信息系统100保险偿付确定的结果(步骤S317)。医院31执行CRT治疗，同时医学治疗成本根据保险偿付确定结果来被覆盖。由此，医院31能够从HSC 32接收对于患者的最有效的CRT设备和CRT设备的最优使用的推荐。
第十一实施例与第二实施例的不同之处在于以下点。图28示出了根据第十一实施例的CRT治疗和医疗成本偿付的过程的示例。根据第十一实施例，提供制造商信息系统400。制造商信息系统400包括设备说明书提供单元410。设备说明书提供单元410将关于CRT设备说明书的信息发送至CRT模拟系统200。CRT模拟系统200设置有设备特征数据存储单元280。设备特征数据存储单元280在其中存储有由制造商信息系统400提供的CRT设备的说明书信息。通过参照设备特征数据存储单元280，心脏模拟器270根据每个模拟目标CRT设备来设置参数并且运行心脏模拟(步骤S126)。
根据第十一实施例，CRT模拟系统200将由各个CRT设备制造商生产的所有CRT设备的说明书保存为数据库。因此，可以向医院31和医疗保险公司33提供选自所注册的CRT设备的最优CRT设备和所选择的CRT设备的最优使用的推荐。
(l)第十二实施例
接下来描述第十二实施例。第十二实施例涉及在CRT设备制造商的发起下进行的CRT治疗过程。图29示出了根据第十二实施例的处理过程的示例。在医生确定CRT治疗在治疗患者的病况方面有效的情况下，医院信息系统100将CRT模拟申请连同患者数据发送至CRT设备制造商43a的制造商信息系统400(步骤S321)。制造商信息系统400将模拟请求连同患者数据发送至HSC 32的CRT模拟系统200(步骤S322)。可以将患者数据从医院信息系统100直接发送至CRT模拟系统200(步骤S323)。
HSC 32的CRT模拟系统200运行CRT模拟然后将模拟的结果发送至医院31(步骤S324)。可以经由CRT设备制造商34a来将模拟结果发送至医院31。随后，医院31使用医院信息系统100来安排给CRT设备制造商34a的模拟成本的支付(步骤S325)。例如，医院信息系统100将覆盖模拟成本的钱款发送至CRT设备制造商34a的银行帐户，并且然后通知制造商信息系统400转账结果。CRT设备制造商34a使用制造商信息系统400来安排给HSC 32的模拟成本的支付(步骤S326)。例如，制造商信息系统400将覆盖模拟成本的钱款发送至HSC 32的银行帐户，并且然后通知CRT模拟系统200转账结果。
在医院31处，医院信息系统100基于模拟结果来确定CRT治疗实施的有效性。如果确定要实施CRT治疗，则医院31接收来自CRT设备制造商34a的CRT设备的递送(步骤S327)。然后，在医院31处，医生执行外科手术过程以将所递送的CRT设备植入到患者中。随后，医院31使用医院信息系统100来向医疗保险公司33做出针对医学治疗成本(包括模拟成本)的要求(步骤S328)。例如，医院信息系统100将医学治疗成本的发票发送至医疗费用偿付系统300。在医疗保险公司33处，医疗费用偿付系统300执行偿付医学治疗成本的处理(步骤S329)。医院31使用医院信息系统100来安排给CRT设备制造商34a的针对CRT设备(包括模拟成本)的支付(步骤S330)。
这样的系统使得CRT设备制造商能够在推荐最优CRT设备方面起到主导作用，同时能够有效地使用CRT模拟。因此，患者能够接受使用对于患者最优的CRT设备的CRT治疗，从而确保对患者的显著治疗效果。
(m)第十三实施例
接下来描述第十三实施例。第十三实施例涉及：当存在多个提供CRT设备的制造商时，运行其中反映每个制造商的CRT设备的特征的CRT模拟，以从而选择对于患者最优的CRT设备。即，不一定仅由一个制造商来提供CRT设备；多个制造商生产全部具有它们的独特特征的CRT设备，并且将它们提供给医院31。有鉴于此，根据第十三实施例，HSC 32基于来自医院31的请求来选择对于患者最优的CRT设备，并且向医院31做出关于CRT设备的最优使用的推荐。
图30示出了根据第十三实施例的处理过程的示例。根据第十三实施例，首先从CRT设备制造商34a和34b中的每个向HSC 32提供对于制造商独特的说明书的足够的公开，并且然后HSC 32准备用于确定制造商的CRT设备的最优使用的模拟并且使其可执行。就此而言，将由CRT设备制造商34a发售的CRT设备“a”的说明书公开给模拟器开发商35a(步骤S341)。类似地，将由CRT设备制造商34b发售的CRT设备“b”的说明书公开给模拟器开发商35b(步骤S342)。例如，将CRT设备说明书的信息从制造商信息系统400a和400b分别发送至各个模拟器开发商35a和35b的模拟器开发系统500a和500b。
模拟器开发商35a开发被设计用于CRT设备“a”的模拟器，并且将其提供给HSC 32(步骤S343)。模拟器开发商35b开发被设计用于CRT设备“b”的模拟器，并且将其提供给HSC32(步骤S344)。例如，模拟器开发商35a和35b的相应的模拟器开发系统500a和500b中的每个将CRT模拟程序发送至HSC 32的CRT模拟系统200。随后，当医院31向HSC 32做出模拟请求时，医院信息系统100将患者数据发送至CRT模拟系统200并且安排给CRT模拟系统200的针对CRT模拟服务的支付(步骤S345)。CRT模拟系统200基于患者数据、使用被设计用于各个CRT设备的模拟器来运行CRT模拟。然后，CRT模拟系统200安排给模拟器开发商35a和35b中的每个的针对模拟器的使用的支付(步骤S346和S347)。
在完成CRT模拟之后，CRT模拟系统200将模拟的结果发送至医院信息系统100(步骤S348)。模拟结果包括关于最优CRT设备以及如何应用CRT设备的最优方法(例如，在四个电极设置在冠状窦导线上的情况下的电极安装位置)的推荐。医院31基于CRT模拟结果来订购最优CRT设备。模拟结果还包括关于每个CRT设备的被预测产生患者的心脏功能的最大改善的电极位置的组合模式，以及与所有组合模式中的每个组合模式相关联的心脏功能的数据(如ECG数据、射血分数、(dP/dt)max和数据可视化心室运动)。已经接收到订单的CRT设备制造商(CRT设备制造商34a或34b)将所订购的CRT设备出售给医院31(步骤S349或S350)。医院31使用医院信息系统100来安排给CRT设备制造商(CRT设备制造商34a或34b)的针对CRT设备的支付(步骤S351或S352)。然后，医院31使用基于CRT模拟结果的方法来实施CRT治疗。随后，医院31使用医院信息系统100来向医疗保险公司33提出针对医学治疗成本包括模拟成本的偿付的要求(步骤S353)。在医疗保险公司33处，医疗费用偿付系统300安排针对医学治疗成本的支付(步骤S354)。
在上述方式中，可以在由多个CRT设备制造商34a和34b生产并发售的CRT设备中选择最优CRT设备，并且由医疗保险来覆盖用于做出选择的模拟的成本。
(n)第十四实施例
接下来描述第十四实施例。根据第十四实施例，不是医院31而是医疗保险公司33预先在多个CRT设备制造商的CRT设备中选择最优CRT设备。医疗保险公司33出于每个CRT设备的成本有效性的考虑来做出选择，这刺激CRT设备制造商之间的价格竞争。
第十四实施例与第十三实施例的不同之处在于以下几点。图31示出了根据第十四实施例的处理过程的示例。根据第十四实施例，医疗保险公司33的医疗费用偿付系统300做出模拟请求并且安排给HSC 32的针对CRT模拟服务的支付(步骤S360)。然后，CRT模拟系统200向医疗费用偿付系统300发送模拟的结果以及关于最优CRT设备和CRT设备的最优使用的信息(步骤S361)。医疗费用偿付系统300然后通过将除了由模拟结果表示的每个CRT设备的例如射血分数的改善率之外的CRT设备的价格考虑在内来做出评估，以从而确定最优CRT设备。然后，医疗费用偿付系统300向医院信息系统100给出最优CRT设备的推荐(步骤S362)。随后，与第十三实施例的步骤S349至S352中相同，做出CRT设备销售和针对CRT设备的支付(步骤S363至S366)。作为响应，医疗费用偿付系统300处理偿付医学治疗成本的过程(步骤S367)。
在上述方式中，通过使得医疗保险公司33能够选择最优CRT设备，也将CRT设备的价格包括在考虑范围内，这导致使用成本有效的CRT设备的CRT治疗，并且由此实现医学治疗成本的下降。
(o)第十五实施例
接下来描述第十五实施例。根据第十五实施例，医生确定最优CRT设备。医院信息系统100预先确定对于患者的最优CRT设备，并且要求所确定的CRT设备的制造商运行CRT模拟。替选地，医院信息系统100可以要求多个CRT设备制造商提交模拟结果并且然后基于结果来确定对于患者最优的CRT设备。
第十五实施例与第十三实施例的不同之处在于以下几点。图32示出了根据第十五实施例的处理过程的示例。根据第十五实施例，医院信息系统100在由多个CRT设备制造商34a和34b生产的CRT设备中确定最优CRT设备。然后，医院信息系统100安排给生产所确定的CRT设备的CRT设备制造商的支付(步骤S381或S382)，并且请求CRT设备制造商递送具有模拟结果的CRT设备。响应于该请求，CRT设备制造商的制造商信息系统通过支付模拟成本来请求CRT模拟系统200运行CRT模拟(步骤S383或步骤S384)。CRT模拟系统200使用被设计用于CRT设备的模拟器来运行CRT模拟，并且然后将模拟的结果发送至制造商信息系统(步骤S385或S386)。制造商信息系统将具有模拟结果的CRT设备递送到医院信息系统100(步骤S387或S388)。在医院31处，执行所递送的CRT设备的植入，并且医院信息系统100向医疗费用偿付系统300做出针对CRT成本和模拟成本的要求(步骤S389)。作为响应，医疗费用偿付系统300安排针对费用的支付(步骤S390)。因此，医生选择适于患者的CRT设备，同时基于医学专业知识来解释模拟结果。
(p)第十六实施例
接下来描述第十六实施例。第十六实施例涉及在医疗保险公司33内部执行CRT模拟。图33示出了根据第十六实施例的处理过程的示例。医院信息系统100将患者数据发送至医疗费用偿付系统300(步骤S401)。医疗费用偿付系统300基于患者数据来运行CRT模拟(步骤S402)并且确定是否由保险来偿付医学治疗成本(步骤S403)。然后，医疗费用偿付系统300向医院信息系统100发送包括关于是否偿付医学治疗成本的信息的模拟的结果(步骤S404)。因此，使用医疗保险公司33的内部系统来执行CRT模拟便利了关于是否做出保险偿付的快速决策。
根据一个方面，可以减少医疗费用。

Claims (13)

1.一种生物医学模拟系统，包括：
生物医学信息存储装置，其用于存储特定于患者的生物医学信息，所述生物医学信息包括所述患者的心脏的层析图像、心脏射血分数和心电图数据；
模拟执行确定装置，其用于基于所述心脏射血分数和所述心电图数据来确定是否执行生物医学模拟；
模型生成装置，其用于在所述模拟执行确定装置确定执行所述生物医学模拟时，基于所述层析图像来生成所述患者的所述心脏的生物医学模型；
模拟执行装置，其用于基于由所述模型生成装置生成的所述生物医学模型来执行对所述患者执行手术之后所述心脏的术后状况的所述生物医学模拟；
改善率确定装置，其用于确定由所述模拟执行装置执行的所述生物医学模拟的结果是否满足与症状相关联的预定改善率；
治疗费用偿付装置，其用于在所述改善率确定装置确定所述生物医学模拟的结果不满足所述预定改善率时确定不偿付用于所述患者的所述手术的医学治疗成本；以及
手术决策装置，其用于在所述改善率确定装置确定所述生物医学模拟的结果满足所述预定改善率时决定对患者执行所述手术。
2.根据权利要求1所述的生物医学模拟系统，其中：
所述生物医学信息存储装置、所述模拟执行确定装置和所述手术决策装置被安装在对所述患者执行所述手术的医疗机构处，
所述模型生成装置、所述模拟执行装置和所述改善率确定装置被安装在模拟中心处，以及
所述治疗费用偿付装置被安装在用于偿付所述医学治疗成本的健康保险公司处。
3.根据权利要求1所述的生物医学模拟系统，其中：
所述生物医学信息存储装置和所述模拟执行确定装置被安装在对所述患者执行所述手术的医疗机构处，
所述模型生成装置、所述模拟执行装置和所述改善率确定装置被安装在模拟中心处，以及
所述手术决策装置和所述治疗费用偿付装置被安装在用于偿付所述医学治疗成本的健康保险公司处。
4.根据权利要求1所述的生物医学模拟系统，其中：
所述生物医学信息存储装置、所述模拟执行确定装置和所述手术决策装置被安装在对所述患者执行所述手术的医疗机构处，
所述模型生成装置和所述模拟执行装置被安装在模拟中心处，以及
所述改善率确定装置和所述治疗费用偿付装置被安装在用于偿付所述医学治疗成本的健康保险公司处。
5.根据权利要求1所述的生物医学模拟系统，其中：
所述生物医学信息存储装置和所述模拟执行确定装置被安装在对所述患者执行所述手术的医疗机构处，
所述模型生成装置和所述模拟执行装置被安装在模拟中心处，以及
所述改善率确定装置、所述手术决策装置和所述治疗费用偿付装置被安装在用于偿付所述医学治疗成本的健康保险公司处。
6.根据权利要求1所述的生物医学模拟系统，其中：
所述手术决策装置包括：用于在所述改善率确定装置确定所述生物医学模拟的结果满足所述预定改善率时输出保险偿付申请的第一决策装置，以及用于在所述改善率确定装置确定所述生物医学模拟的结果满足所述预定改善率并且然后所述第一决策单元输出了所述保险偿付申请时决定对所述患者执行所述手术的第二决策装置，
所述生物医学信息存储装置、所述模拟执行确定装置和所述第一决策装置被安装在对所述患者执行所述手术的医疗机构处，
所述模型生成装置、所述模拟执行装置和所述改善率确定装置被安装在模拟中心处，以及
所述第二决策装置和所述治疗费用偿付装置被安装在用于偿付所述医学治疗成本的健康保险公司处。
7.根据权利要求1所述的生物医学模拟系统，其中：
所述治疗费用偿付装置在所述改善率确定装置确定所述生物医学模拟的结果满足所述预定改善率时确定偿付所述医学治疗成本，以及
所述手术决策装置在所述改善率确定装置确定所述生物医学模拟的结果满足所述预定改善率并且然后所述手术决策装置获得对所述治疗费用偿付装置确定偿付所述医学治疗成本的确认时决定对所述患者执行所述手术。
8.根据权利要求1所述的生物医学模拟系统，其中：
所述治疗费用偿付装置在对所述患者执行所述手术时偿付所述医学治疗成本和用于执行所述生物医学模拟的成本，并且所述治疗费用偿付装置在不对所述患者执行所述手术时偿付用于执行所述生物医学模拟的成本。
9.根据权利要求1所述的生物医学模拟系统，还包括：
设备特征数据存储装置，其用于存储表示待被用在所述患者的所述手术中的设备的特征的设备特征数据，
其中，所述模拟执行装置参照所述设备特征数据来执行其中反映了所述设备的特征的所述生物医学模拟。
10.根据权利要求1所述的生物医学模拟系统，其中：
所述模拟执行装置关于可用于所述患者的多个手术过程中的每个手术过程来执行所述生物医学模拟以预测在使用所述手术过程的手术之后所述患者的术后状况，以及
所述改善率确定装置在可用手术过程中选择具有表现出最高改善率的所述生物医学模拟的结果的手术过程，并且所述改善率确定装置确定所述最高改善率是否满足所述预定改善率。
11.根据权利要求1所述的生物医学模拟系统，其中：
所述模拟执行装置关于可用于所述患者的多个手术过程中的每个手术过程来执行所述生物医学模拟以预测在使用所述手术过程的手术之后所述患者的术后状况，以及
所述改善率确定装置基于评估值来做出与所述预定改善率相关联的确定，其中，所述评估值是值的和，所述值中的每个值通过将针对每个电极布置模式获得的多个指标中的每个指标的改善率乘以对应于所述指标的加权值来获得，所述指标中的每个指标的改善率被包括在从针对各个手术过程执行的所述生物医学模拟获得的结果中。
12.根据权利要求11所述的生物医学模拟系统，其中：
从医学的角度预先设置对应于所述指标中的每个指标的所述加权值。
13.根据权利要求10所述的生物医学模拟系统，其中：
待被用在所述患者的所述手术中的设备是心脏再同步化治疗CRT设备，以及针对所述手术过程中的每个手术过程的所述生物医学模拟是针对与所述CRT设备的植入相关联的多个电极布置模式中的每个电极布置模式的模拟。
CN201610608292.6A 2015-07-31 2016-07-28 用于生物医学模拟的系统和方法 Active CN106446504B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015-151649 2015-07-31
JP2015151649A JP6631072B2 (ja) 2015-07-31 2015-07-31 生体シミュレーションシステムおよび生体シミュレーション方法

Publications (2)

Publication Number Publication Date
CN106446504A CN106446504A (zh) 2017-02-22
CN106446504B true CN106446504B (zh) 2019-04-09

Family

ID=56411491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610608292.6A Active CN106446504B (zh) 2015-07-31 2016-07-28 用于生物医学模拟的系统和方法

Country Status (4)

Country Link
US (1) US20170032095A1 (zh)
EP (1) EP3125141B1 (zh)
JP (1) JP6631072B2 (zh)
CN (1) CN106446504B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016002054A1 (ja) * 2014-07-03 2016-01-07 富士通株式会社 生体シミュレーション装置、生体シミュレーション装置の制御方法、および生体シミュレーション装置の制御プログラム
CN108428477A (zh) * 2018-03-30 2018-08-21 北京航空航天大学 基于数字双胞胎的云医疗仿真平台的构建方法及云医疗系统
BR112020021271A2 (pt) 2018-04-17 2021-01-26 Celldex Therapeutics, Inc. construtos bispecíficos e anticorpos anti- cd27 e anti-pd-l1
WO2020013198A1 (ja) * 2018-07-13 2020-01-16 富士フイルム株式会社 画像診断支援装置、方法及びプログラム、並びに心臓シミュレーションシステム
JP6566409B1 (ja) * 2018-08-24 2019-08-28 株式会社鈴康 情報処理装置、プログラム及び情報処理方法
CN110970133B (zh) * 2019-12-10 2023-03-21 中国医学科学院肿瘤医院 一种crt风险评估方法和风险预测系统
JP7539430B2 (ja) 2022-06-30 2024-08-23 富士通株式会社 立体形状データ生成プログラム、立体形状データ生成方法、および情報処理装置
WO2024038618A1 (ja) 2022-08-19 2024-02-22 株式会社ＵＴ−Ｈｅａｒｔ研究所 インシリコ心疾患データベース活用方法、インシリコ心疾患データベース活用プログラム、および情報処理装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1670755A (zh) * 2004-03-17 2005-09-21 希森美康株式会社 糖尿病诊疗支持系统以及存储介质
CN101077307A (zh) * 2006-05-24 2007-11-28 希森美康株式会社 生物器官功能模拟系统
CN101627411A (zh) * 2007-01-16 2010-01-13 西姆博尼克斯有限公司 手术前的手术模拟
CN102117378A (zh) * 2009-12-31 2011-07-06 苏州瑞派宁科技有限公司 一种基于三维多模影像的肝肿瘤综合手术规划模拟方法及其系统
CN102982238A (zh) * 2012-11-13 2013-03-20 深圳市环球博洋机电科技有限公司 一种基于ct胶片的三维肾肿瘤手术模拟方法及其平台
CN103270513A (zh) * 2010-08-12 2013-08-28 哈特弗罗公司 用于患者特异性血流建模的方法和系统

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001331581A (ja) * 2000-05-22 2001-11-30 Yusuke Tsukamoto 自動診断システム、自動診断方法、治療法自動決定システム、治療法自動決定方法及び記録媒体
JP4108449B2 (ja) * 2002-11-13 2008-06-25 帝人株式会社 酸素療法の治療効果予測方法、酸素療法の実施支援方法
JP2004280807A (ja) * 2003-02-28 2004-10-07 Toshiba Corp サイバーホスピタルシステム
US8036743B2 (en) 2005-03-31 2011-10-11 Proteus Biomedical, Inc. Automated optimization of multi-electrode pacing for cardiac resynchronization
US8936192B2 (en) * 2005-12-20 2015-01-20 United States Postal Service Method and system for interrogating and processing codes
JP4901226B2 (ja) * 2006-01-27 2012-03-21 シスメックス株式会社 医療用シミュレーションシステム及びそのコンピュータプログラム
JP2007312922A (ja) * 2006-05-24 2007-12-06 Sysmex Corp 生体応答予測システム及びそのプログラム
US8224432B2 (en) 2007-05-08 2012-07-17 C.R.Bard, Inc. Rapid 3D mapping using multielectrode position data
JP2008272510A (ja) * 2008-07-14 2008-11-13 Toshiba Corp 治療法の有効性を予測するためのプログラム、データベース、システム及び方法
WO2010132617A2 (en) * 2009-05-12 2010-11-18 Chronicmobile, Inc. Methods and systems for managing, controlling and monitoring medical devices via one or more software applications functioning in a secure environment
JP2011098175A (ja) 2009-11-04 2011-05-19 Naotoshi Maeda 心臓再同期療法（ｃａｒｄｉａｃｒｅｓｙｎｃｈｒｏｎｉｚａｔｉｏｎｔｈｅｒａｐｙ，ＣＲＴ）の適応患者を選択するための検査法
WO2012109618A2 (en) * 2011-02-11 2012-08-16 The Johns Hopkins University System and method for planning a patient-specific cardiac procedure
US10311978B2 (en) * 2012-01-30 2019-06-04 Siemens Healthcare Gmbh Method and system for patient specific planning of cardiac therapies on preoperative clinical data and medical images
US9414776B2 (en) * 2013-03-06 2016-08-16 Navigated Technologies, LLC Patient permission-based mobile health-linked information collection and exchange systems and methods
US9463072B2 (en) * 2013-08-09 2016-10-11 Siemens Aktiengesellschaft System and method for patient specific planning and guidance of electrophysiology interventions
WO2015073927A2 (en) * 2013-11-15 2015-05-21 The Regents Of The University Of California Compositions, devices and methods for diagnosing heart failure and for patient-specific modeling to predict outcomes of cardiac resynchronization therapy
EP3088897B1 (en) * 2013-12-04 2020-11-04 Kazuyuki Yoshizaki Method for predicting therapeutic effect of biological preparation on rheumatoid arthritis

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1670755A (zh) * 2004-03-17 2005-09-21 希森美康株式会社 糖尿病诊疗支持系统以及存储介质
CN101077307A (zh) * 2006-05-24 2007-11-28 希森美康株式会社 生物器官功能模拟系统
CN101627411A (zh) * 2007-01-16 2010-01-13 西姆博尼克斯有限公司 手术前的手术模拟
CN102117378A (zh) * 2009-12-31 2011-07-06 苏州瑞派宁科技有限公司 一种基于三维多模影像的肝肿瘤综合手术规划模拟方法及其系统
CN103270513A (zh) * 2010-08-12 2013-08-28 哈特弗罗公司 用于患者特异性血流建模的方法和系统
CN102982238A (zh) * 2012-11-13 2013-03-20 深圳市环球博洋机电科技有限公司 一种基于ct胶片的三维肾肿瘤手术模拟方法及其平台

Also Published As

Publication number Publication date
EP3125141B1 (en) 2020-05-20
JP2017033227A (ja) 2017-02-09
JP6631072B2 (ja) 2020-01-15
US20170032095A1 (en) 2017-02-02
CN106446504A (zh) 2017-02-22
EP3125141A1 (en) 2017-02-01

Similar Documents

Publication Publication Date Title
CN106446504B (zh) 用于生物医学模拟的系统和方法
Cluitmans et al. Validation and opportunities of electrocardiographic imaging: from technical achievements to clinical applications
Boyle et al. Computationally guided personalized targeted ablation of persistent atrial fibrillation
US11576626B2 (en) Systems and methods for numerically evaluating vasculature
Relan et al. Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia
US12089918B2 (en) System and method for machine learning based mapping function of patient data from one physiological state to another based on computational fluid dynamics simulation of blood via offline training
US11342080B2 (en) Method and system for disease progression modeling and therapy optimizing for individual patients
US20080208068A1 (en) Dynamic positional information constrained heart model
CN108962381B (zh) 用于动脉粥样硬化的个性化评定、长期预测和管理的基于学习的方法
US20110307231A1 (en) Method and arrangement for creating an individualized, computer-aided model of a system, and a corresponding computer program and a corresponding machine-readable storage medium
JP2014534889A (ja) 冠循環のマルチスケールな解剖学的かつ機能的なモデリングの方法およびシステム
JP2020182848A (ja) 心臓アブレーション処置を改善するためのシステム及び方法
JP2017525497A (ja) 血流特性及び病理を心筋血液供給のモデリングを介して決定するためのシステム及び方法
Ge et al. Comparison of Instantaneous Wave‐Free Ratio (iFR) and Fractional Flow Reserve (FFR) with respect to their sensitivities to cardiovascular factors: a computational model‐based study
Pagani et al. Data integration for the numerical simulation of cardiac electrophysiology
Matene et al. Fully automated initiation of simulated episodes of atrial arrhythmias
EP4044927A1 (en) Renal denervation ablation monitoring using perfusion angiography
Sung et al. Evaluation of a deep learning‐enabled automated computational heart modelling workflow for personalized assessment of ventricular arrhythmias
Zhu et al. [Retracted] Artificial Intelligence‐Based Echocardiographic Left Atrial Volume Measurement with Pulmonary Vein Comparison
Loewe et al. Cardiac digital twin modeling
Jia et al. A hybrid catheter localisation framework in echocardiography based on electromagnetic tracking and deep learning segmentation
Falkenberg KIT-Institute of Biomedical Engineering-Research-Publications
Falkenberg KIT-IBT-Forschung-Veröffentlichungen-Publikationen
Huang et al. Quantitative evaluation of the effect of circle of willis structures on cerebral hyperperfusion: A multi-scale model analysis
Galappaththige et al. A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant