• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于臨床數(shù)據(jù)與循環(huán)系統(tǒng)模型融合技術(shù)的心血管功能評(píng)估中的誤差分析

    2016-02-16 07:26:26殷兆芳梁夫友
    關(guān)鍵詞:主要參數(shù)測(cè)量誤差左心室

    李 逸 殷兆芳 梁夫友#*

    1(上海交通大學(xué)船舶海洋與建筑工程學(xué)院,上海 200240)2(上海交通大學(xué)醫(yī)學(xué)院附屬上海第九人民醫(yī)院心內(nèi)科,上海 200011)

    基于臨床數(shù)據(jù)與循環(huán)系統(tǒng)模型融合技術(shù)的心血管功能評(píng)估中的誤差分析

    李 逸1殷兆芳2梁夫友1#*

    1(上海交通大學(xué)船舶海洋與建筑工程學(xué)院,上海 200240)2(上海交通大學(xué)醫(yī)學(xué)院附屬上海第九人民醫(yī)院心內(nèi)科,上海 200011)

    識(shí)別影響患者個(gè)體化心血管功能評(píng)估(基于臨床數(shù)據(jù)與循環(huán)系統(tǒng)模型的融合技術(shù))精度的關(guān)鍵因子,并量化分析各因子的影響程度,為評(píng)估方法的臨床應(yīng)用提供理論支持。通過(guò)參數(shù)敏感度分析并結(jié)合子集選擇方法,識(shí)別與血流動(dòng)力學(xué)計(jì)算和主要模型參數(shù)(與心血管功能的評(píng)估目標(biāo)對(duì)應(yīng))評(píng)估相關(guān)的次主要模型參數(shù)。利用虛擬臨床數(shù)據(jù)開(kāi)展數(shù)值實(shí)驗(yàn),分析臨床數(shù)據(jù)測(cè)量誤差(誤差區(qū)間5%)和次主要參數(shù)變動(dòng)(變動(dòng)率30%)對(duì)主要模型參數(shù)評(píng)估結(jié)果的影響。臨床數(shù)據(jù)的測(cè)量誤差引起評(píng)估結(jié)果的顯著變化(最大變化率達(dá)到16.6%),相對(duì)地,次主要參數(shù)變動(dòng)造成的評(píng)估誤差較小(在10%以內(nèi))。臨床數(shù)據(jù)的精確采集是確保評(píng)估結(jié)果可靠的關(guān)鍵環(huán)節(jié);次主要參數(shù)雖然在患者間存在實(shí)際的個(gè)體差異性,但對(duì)其進(jìn)行固定處理對(duì)評(píng)估結(jié)果的影響較小。

    心血管功能評(píng)估;循環(huán)系統(tǒng)模型;敏感度分析;子集選擇;誤差分析

    引言

    在對(duì)冠心病和心衰患者的治療中,心血管功能指標(biāo)(如心收縮、舒張功能,外周血管阻抗等)是選擇治療時(shí)機(jī)、進(jìn)行恰當(dāng)?shù)男g(shù)后護(hù)理和用藥的重要依據(jù)[1]。臨床上現(xiàn)存多種測(cè)量或評(píng)估心血管功能的方法,例如:心臟超聲可以測(cè)量心臟的體積、射血分?jǐn)?shù)等反映心臟形態(tài)學(xué)特征和心肌收縮力的數(shù)據(jù)[2-3]; MRI(magnetic resonance imaging)可通過(guò)對(duì)心肌運(yùn)動(dòng)特征的區(qū)域化分割處理,識(shí)別潛在的缺血區(qū)域[4];另外,BNP(B-type natriuretic peptide,亦稱為“客觀心功能標(biāo)志物”)是反映心衰程度的重要生化指標(biāo)[5]。

    人體血液循環(huán)系統(tǒng)是由心臟和多種血管構(gòu)成的閉環(huán)系統(tǒng),各部分的血流參數(shù)交互作用,造成在體測(cè)量某一項(xiàng)特定心血管功能指標(biāo)時(shí)通常難以摒除其他因素的影響。例如,心臟超聲測(cè)量的心室射血分?jǐn)?shù)(EF)被認(rèn)為可以反映心室的收縮力,測(cè)量的二尖瓣血流波形頻譜可用于判斷心室的舒張功能。然而,從系統(tǒng)血流動(dòng)力學(xué)的角度分析,EF和二尖瓣血流波形均受心臟前后負(fù)荷(如心室灌注壓、外周血管阻抗)的顯著影響。特別地,在心室舒張功能異常的進(jìn)展期會(huì)出現(xiàn)二尖瓣血流頻譜假性正?;F(xiàn)象,增大了心功能異常的診斷難度。Senzaki等的研究亦表明,分離評(píng)估心功能和心臟的負(fù)荷條件,對(duì)更好地把握患者的心血管功能狀態(tài)和優(yōu)化治療有重要的意義[6]。BNP雖然對(duì)診斷心衰有重要的意義,但無(wú)法有效識(shí)別心收縮或舒張功能不全[7]。迄今為止,關(guān)于心臟的固有功能屬性在臨床上仍沒(méi)有明確的量化評(píng)價(jià)指標(biāo)和對(duì)應(yīng)的測(cè)量手段,但在生理學(xué)領(lǐng)域已有研究者提出量化分析方法(即elastance理論)[8],可實(shí)現(xiàn)對(duì)心臟收縮、舒張功能的量化、分離評(píng)估。然而,對(duì)elastance的評(píng)估通常需要心室容積與壓力的時(shí)變關(guān)系曲線信息[9],臨床可操作性差。

    利用循環(huán)系統(tǒng)的血流動(dòng)力學(xué)模型,通過(guò)逆向工程學(xué)的方法分析和整合離散的臨床數(shù)據(jù),可為上述問(wèn)題的解決提供可行的途徑。前期研究已從數(shù)值實(shí)驗(yàn)[10]、動(dòng)物實(shí)驗(yàn)[11-13]和體外模型實(shí)驗(yàn)[14]等多個(gè)方面對(duì)該方法進(jìn)行了探討,但缺乏對(duì)實(shí)際臨床數(shù)據(jù)的分析和方法學(xué)驗(yàn)證。筆者的研究組在最近的研究中開(kāi)發(fā)了通過(guò)融合循環(huán)系統(tǒng)模型與臨床數(shù)據(jù)來(lái)定量評(píng)估患者心血管功能的新算法,并針對(duì)少量冠脈搭橋和單心室心臟病患者證明了其可行性[15-16]。在上述研究中,研究者為了克服可利用臨床數(shù)據(jù)資源不足這一問(wèn)題,采用參數(shù)敏感度分析法,從模型參數(shù)中識(shí)別出與臨床測(cè)量數(shù)據(jù)(如大動(dòng)脈血壓、肺動(dòng)脈血壓、心容積等)密切相關(guān)的參數(shù)(定義為主要參數(shù),其他參數(shù)為次要參數(shù)),并將其作為變量來(lái)構(gòu)建參數(shù)優(yōu)化的目標(biāo)函數(shù)。在目標(biāo)函數(shù)的求解過(guò)程中,對(duì)次要參數(shù)進(jìn)行固定處理是保證函數(shù)的解存在且唯一的必要條件。然而,一方面,從生理、病理學(xué)的角度考慮,次要參數(shù)所代表的心血管功能在不同患者間存在著客觀差異性;另一方面,血壓、心容積等臨床數(shù)據(jù)包含測(cè)量誤差,這些因素對(duì)求解結(jié)果的影響有待進(jìn)一步研究。

    本研究的目的在于考察臨床數(shù)據(jù)測(cè)量誤差和次要參數(shù)變動(dòng)對(duì)心血管功能參數(shù)評(píng)估結(jié)果的影響,研究結(jié)果將有助于進(jìn)一步檢驗(yàn)評(píng)估算法的可靠性,標(biāo)定評(píng)估結(jié)果的誤差范圍,探明其主導(dǎo)因素,從而為該方法的臨床應(yīng)用奠定理論基礎(chǔ)。

    1 方法

    為了定量地考察評(píng)估結(jié)果受測(cè)量誤差和次要參數(shù)變動(dòng)的影響,采用循環(huán)系統(tǒng)模型進(jìn)行隨機(jī)模擬,生成99組“虛擬患者”數(shù)據(jù),導(dǎo)出肺動(dòng)脈平均壓(Pmpp)、左心房平均壓(Plap)、大動(dòng)脈平均壓(Pmbp)、左心室收縮末期體積(Vesv)和左心室舒張末期體積(Vedv)等“臨床測(cè)量數(shù)據(jù)”,以及對(duì)應(yīng)的主要模型參數(shù)值,將其設(shè)為數(shù)值實(shí)驗(yàn)的參照系。另外,鑒于次要參數(shù)的數(shù)量較大,采用參數(shù)敏感度分析法并結(jié)合參數(shù)間的結(jié)構(gòu)相關(guān)度,分析篩選出對(duì)模型計(jì)算結(jié)果和主要參數(shù)評(píng)估最有影響的數(shù)個(gè)參數(shù),將其定義為次主要參數(shù)。對(duì)有限數(shù)量的次主要參數(shù)進(jìn)行數(shù)值實(shí)驗(yàn),有助于在簡(jiǎn)化誤差分析的同時(shí)降低計(jì)算負(fù)荷。

    1.1 循環(huán)系統(tǒng)的計(jì)算模型

    采用研究組在前期研究中使用的循環(huán)系統(tǒng)模型[17]。該模型使用一系列集中參數(shù),表達(dá)心臟的收縮、舒張功能以及血管的黏性阻抗、順應(yīng)性等特性,其電路類比如圖1所示。該模型涵蓋了循環(huán)系統(tǒng)的各主要部分(如心臟四腔、體循環(huán)和肺循環(huán)),可保障對(duì)心臟與其前后負(fù)荷耦合過(guò)程的合理描述;同時(shí),模型中的每個(gè)參數(shù)均與特定的心血管功能對(duì)應(yīng),有利于定量地分析血壓、血流量等血流動(dòng)力學(xué)參數(shù)與心血管功能的關(guān)系。與更高維度(如一維或三維)的建模方法比較,雖然集中參數(shù)建模法無(wú)法精確描述局部的血流特征和動(dòng)脈系統(tǒng)內(nèi)的血壓波傳播等現(xiàn)象[18],但它是唯一可以描述循環(huán)系統(tǒng)整體血流動(dòng)態(tài)的方法,而且模型的參數(shù)數(shù)量可控,利于參數(shù)篩選和基于有限的臨床數(shù)據(jù)進(jìn)行參數(shù)校準(zhǔn)[19-22]。

    圖1 循環(huán)系統(tǒng)集中參數(shù)模型的電路類比[17]。參數(shù)標(biāo)注:E—倒電容;L—慣性系數(shù);R—黏性阻抗;C—順應(yīng)性;B—伯努利阻抗;S—黏彈性系數(shù);Pit—胸內(nèi)壓;Ppc—心包壓。下標(biāo):ra—右心房;rv—右心室;pua—肺動(dòng)脈;puc—肺毛細(xì)血管;puv—肺靜脈;la—左心房;lv—左心室;av—主動(dòng)脈瓣;mv—二尖瓣;tv—三尖瓣;pv—肺動(dòng)脈瓣;ao—主動(dòng)脈;art—?jiǎng)用};cap—毛細(xì)血管;ven—靜脈;ivc—下腔靜脈;svc—上腔靜脈;_l—下半身;_u—上半身Fig.1 Electric circuit analogy of a lumped parameter model of the cardiovascular system[17]. Notation of parameters: E, elastance; L, inertance; R, viscous resistance; C, compliance; B, Bernoulli’s resistance; S, viscoelasticity coefficient; Pit, intra-thoracic pressure; Ppc, pericardial pressure. Subscripts: ra, right atrium; rv, right ventricle; pua, pulmonary artery; puc, pulmonary capillary; puv, pulmonary vein; la, left atrium; lv, left ventricle; av, aortic vavle; mv, mitral valve; tv, tricuspid valve; pv, pulmonary valve; ao, aorta; art, arteries; cap, capillary bed; ven, veins; ivc, inferior vena cava; svc, superior vena cava; _l, lower body; _u, upper body

    模型支配方程的構(gòu)建遵循基本的血流動(dòng)力學(xué)原理[15-17]。心臟的4個(gè)房室的收縮和舒張功能采用時(shí)變的倒電容(即elastance)來(lái)表達(dá),通過(guò)心臟內(nèi)的血液容積計(jì)算心臟內(nèi)壓,并進(jìn)一步通過(guò)心臟瓣膜模型與外周血管的血流動(dòng)力學(xué)參數(shù)關(guān)聯(lián)。血管內(nèi)血流的支配方程基于流經(jīng)每個(gè)血管區(qū)段血流的質(zhì)量和動(dòng)量守恒原理構(gòu)建,所有方程構(gòu)成一個(gè)非線性的常微分方程組,采用四階龍格-庫(kù)塔法求解。關(guān)于模型支配方程和數(shù)值算法的更詳細(xì)信息,可參考研究組的前期工作[15-17]。

    1.2 臨床數(shù)據(jù)及患者個(gè)體化評(píng)估

    1.2.1 臨床數(shù)據(jù)

    本研究主要利用心臟病診療流程中采集的臨床數(shù)據(jù),包括大動(dòng)脈壓(由動(dòng)脈導(dǎo)管采集)、肺動(dòng)脈壓、肺毛細(xì)血管楔壓(由靜脈導(dǎo)管采集)和左心室體積(由心臟超聲測(cè)量)。為了使模型可以精確地?cái)M合臨床數(shù)據(jù),有必要分析臨床數(shù)據(jù)與模型計(jì)算結(jié)果的對(duì)應(yīng)關(guān)系。例如,考慮集中參數(shù)模型描述血管內(nèi)脈搏波傳播的局限性,取大動(dòng)脈壓、肺動(dòng)脈壓在一個(gè)心動(dòng)周期內(nèi)的平均值而非血壓波形作為擬合目標(biāo);另外,由于本模型未包含肺毛細(xì)血管楔壓的計(jì)算模塊,假定肺毛細(xì)血管楔壓與左心房平均壓一致(可由模型計(jì)算獲得)。最終用于模型參數(shù)評(píng)估的臨床數(shù)據(jù)包括:肺動(dòng)脈平均壓(Pmpp)、左心房平均壓(Plap)、大動(dòng)脈平均壓(Pmbp)、左心室收縮末期體積(Vesv)和左心室舒張末期體積(Vedv)。

    1.2.2 主要參數(shù)識(shí)別與參數(shù)優(yōu)化

    采用的循環(huán)系統(tǒng)模型由58個(gè)參數(shù)構(gòu)成,而可利用的臨床數(shù)據(jù)僅有5項(xiàng),即基于臨床數(shù)據(jù)構(gòu)建的5個(gè)參數(shù)優(yōu)化目標(biāo)函數(shù)包含58個(gè)變量。根據(jù)代數(shù)方程理論,該方程組不封閉,沒(méi)有唯一解存在。本研究所開(kāi)展的數(shù)值實(shí)驗(yàn)也證明,目標(biāo)函數(shù)的解隨模型初始條件的改變而變化,無(wú)法唯一地確定模型參數(shù)。為了解決該問(wèn)題,對(duì)模型系統(tǒng)進(jìn)行了參數(shù)敏感度分析,從中選取與臨床數(shù)據(jù)的模型計(jì)算結(jié)果最相關(guān)的數(shù)個(gè)模型參數(shù),設(shè)定為主要參數(shù),即體循環(huán)血管順應(yīng)性(Csys)、左心室基準(zhǔn)elastance(Elvb)、左心室收縮期elastance振幅(Elva)、體循環(huán)阻抗(Rsys)、肺循環(huán)阻抗(Rpul),其他參數(shù)設(shè)為次要參數(shù)。根據(jù)臨床數(shù)據(jù)和主要參數(shù),構(gòu)建的參數(shù)優(yōu)化目標(biāo)函數(shù)為

    (1)

    式中,Xe(xe,1,xe,2…xe,5)和X0(x0,1,x0,2…x0,53)分別表示主要參數(shù)向量和次要參數(shù)向量,在此,Xe為目標(biāo)函數(shù)的變量,X0為常數(shù)(默認(rèn)值),Pj和Yj分別表示模型計(jì)算結(jié)果和臨床實(shí)測(cè)值,對(duì)應(yīng)著肺動(dòng)脈平均壓和左心室體積等。

    由于目標(biāo)函數(shù)是包含多個(gè)變量的非線性方程組,本研究采用改良的Nelder-Mead法進(jìn)行求解[16],求解所得Xe值即為與臨床數(shù)據(jù)(Yj)對(duì)應(yīng)的患者個(gè)體化心血管功能參數(shù)。

    1.3 虛擬臨床數(shù)據(jù)的生成

    如前所述,筆者的主要目的是研究臨床數(shù)據(jù)測(cè)量誤差和次要參數(shù)變動(dòng)對(duì)主要參數(shù)估測(cè)值的影響,為此目的,需要有主要參數(shù)值作為已知條件和參照系,以便開(kāi)展定量分析。然而,在實(shí)測(cè)的臨床數(shù)據(jù)中缺乏心血管功能參數(shù)信息,無(wú)法開(kāi)展上述研究。為了解決這一問(wèn)題,采用循環(huán)系統(tǒng)模型生成一系列虛擬的臨床數(shù)據(jù),并記錄其所對(duì)應(yīng)的心血管功能參數(shù),構(gòu)建分析估測(cè)誤差的參照系。具體來(lái)說(shuō),在模型計(jì)算中,對(duì)主要模型參數(shù)在其生理范圍內(nèi)進(jìn)行隨機(jī)變動(dòng)(相對(duì)其默認(rèn)值在20%~200%間變動(dòng)),并求解肺動(dòng)脈平均壓、左心房壓和左心室體積等,將其作為虛擬的臨床數(shù)據(jù)。利用該方法,共生成了符合統(tǒng)計(jì)學(xué)要求的99組虛擬臨床數(shù)據(jù),每組數(shù)據(jù)均包含病人標(biāo)號(hào)(n)、性別(sex)、體重(wight)、身高(hight)、心動(dòng)周期(Tduration)、肺動(dòng)脈平均壓(Pmpp)、左心房平均壓(Plap)、大動(dòng)脈平均壓(Pmbp)、左心室收縮末期體積(Vesv)、左心室舒張末期體積(Vedv)等信息。各項(xiàng)虛擬數(shù)據(jù)的分布區(qū)間和統(tǒng)計(jì)結(jié)果(表示為mean±SD)如下:肺動(dòng)脈平均壓[12.1,37.2](22.7±6.7)[mmHg],左心房壓[7.9,30.9](15.3±5.5)[mmHg],大動(dòng)脈平均壓[33.2,190.8](98.3±38.7)[mmHg],左心室收縮末期體積[12.6,217.7](57.9±36.1)[ml],左心室舒張末期體積 [109.3,463.2](217.7±74.7)[ml]。上述數(shù)據(jù)覆蓋了臨床所見(jiàn)血壓和心室容積的主要變動(dòng)范圍,并且不同組間的數(shù)值變異顯著,因此可以較好地體現(xiàn)實(shí)際患者間的個(gè)體差異性。

    1.4 次主要參數(shù)的篩選

    本研究采用參數(shù)敏感度分析和子集選擇相結(jié)合的方法,識(shí)別對(duì)臨床數(shù)據(jù)的模型仿真結(jié)果影響較大同時(shí)與主要參數(shù)相關(guān)的次主要參數(shù)。

    1.4.1 參數(shù)敏感度分析

    由于循環(huán)系統(tǒng)是復(fù)雜的非線性系統(tǒng),模型參數(shù)的敏感性與系統(tǒng)狀態(tài)有關(guān),因此在單一心血管功能狀態(tài)下開(kāi)展的參數(shù)敏感度分析有漏選某些重要參數(shù)的風(fēng)險(xiǎn)[16]。為了解決上述問(wèn)題,首先采用Latin-Hypercube全域參數(shù)取樣技術(shù)[23],對(duì)關(guān)鍵模型參數(shù)(5個(gè)主要參數(shù)+20個(gè)次要參數(shù)(根據(jù)其生理意義選擇))在其生理范圍內(nèi)進(jìn)行隨機(jī)變動(dòng),生成100組模型參數(shù)子集來(lái)表達(dá)多種心血管功能狀態(tài),從而構(gòu)成進(jìn)行全域參數(shù)敏感度分析的參數(shù)空間。需要指出,雖然進(jìn)一步增加參數(shù)組合可以提高數(shù)據(jù)的分割精度(研究組在前期工作中采用了3000組數(shù)據(jù)[16]),但會(huì)增大敏感度矩陣的復(fù)雜性及后續(xù)的數(shù)據(jù)分析難度。Latin-Hypercube取樣技術(shù)的優(yōu)勢(shì)在于可以利用有限數(shù)量的數(shù)據(jù)構(gòu)成反映參數(shù)全域變動(dòng)屬性的數(shù)據(jù)空間,本研究開(kāi)展的數(shù)值實(shí)驗(yàn)亦表明,在參數(shù)分布區(qū)間不變的條件下,將參數(shù)組合由100組增大至1 000組識(shí)別出的敏感參數(shù)不變。在上述參數(shù)空間內(nèi),通過(guò)模型仿真考察每個(gè)參數(shù)變動(dòng)(本研究采用的變化率為Δ=5%)而引起的肺動(dòng)脈平均壓(Pmpp)、左心房平均壓(Plap)等的變化率。最后,通過(guò)統(tǒng)計(jì)學(xué)分析計(jì)算Pmpp、Plap等對(duì)每個(gè)參數(shù)的敏感度指數(shù)。

    對(duì)每項(xiàng)臨床數(shù)據(jù),其參數(shù)敏感度指數(shù)定義為

    (2)

    式中,y表示與臨床數(shù)據(jù)對(duì)應(yīng)的模型計(jì)算值,θ表示模型參數(shù),下標(biāo)i(=1,M),j(=1,N)分別表示模型參數(shù)子集的序列號(hào)和各參數(shù)在模型中的標(biāo)號(hào),得到敏感度矩陣表示為

    (3)

    式中,M=100,N=58。

    表1 次主要參數(shù)篩選列表Table 1 List of secondary major parameters

    1.4.2 子集選擇

    雖然準(zhǔn)次主要參數(shù)對(duì)臨床數(shù)據(jù)的模型計(jì)算結(jié)果影響較大,但對(duì)主要參數(shù)估測(cè)結(jié)果的影響程度未知。為了明確其對(duì)主要參數(shù)估測(cè)結(jié)果的影響,同時(shí)進(jìn)一步減少擬分析次要參數(shù)的數(shù)量,本研究引入了子集選擇方法。子集選擇有結(jié)構(gòu)相關(guān)度分析、特征值分解以及特征向量子空間分析等方法[24]??紤]本研究的目的,采用了結(jié)構(gòu)相關(guān)度分析法計(jì)算所有模型參數(shù)與主要參數(shù)的相關(guān)度γ,γ值越大,則兩參數(shù)間的相關(guān)度越大。

    (4)

    1.5 對(duì)主要參數(shù)評(píng)估結(jié)果的影響

    1.5.1 臨床測(cè)量誤差

    將虛擬臨床數(shù)據(jù)的每項(xiàng)數(shù)據(jù)(即肺動(dòng)脈平均壓(Pmpp)、左心房平均壓(Plap)、大動(dòng)脈平均壓(Pmbp)、左心室收縮末期體積(Vesv)和左心室舒張末期體積(Vedv))分別施加5%的測(cè)量誤差(臨床測(cè)量誤差一般處于0~5%之間[12]),利用參數(shù)優(yōu)化算法重新評(píng)估主要參數(shù)(即體循環(huán)血管順應(yīng)性(Csys)、左心室基準(zhǔn)elastance(Elvb,反映左心室的舒張功能)、左心室收縮期elastance振幅(Elva,反映左心室的收縮功能)、體循環(huán)阻抗(Rsys)、肺循環(huán)阻抗(Rpul),并與主要參數(shù)的參照值比較,計(jì)算評(píng)估結(jié)果的變化率(以百分比計(jì))。最后,對(duì)基于99組虛擬臨床數(shù)據(jù)的求解結(jié)果進(jìn)行統(tǒng)計(jì)學(xué)分析,計(jì)算變化率的平均值及標(biāo)準(zhǔn)差(mean±SD)。

    1.5.2 次主要參數(shù)變動(dòng)

    對(duì)表1中的7個(gè)次主要參數(shù)(即右心室基準(zhǔn)elastance(Ervb)、右心房基準(zhǔn)elastance(Erab)、右心室收縮期elastance振幅(Erva)、左心房收縮期elastance振幅(Elaa)、下腔靜脈阻抗(Rvca_l)、左心房基準(zhǔn)elastance(Elab)、右心房收縮期elastance振幅(Eraa)),在其默認(rèn)值的基礎(chǔ)上分別施加30%的變化,同時(shí)保持虛擬臨床數(shù)據(jù)的值不變,重新評(píng)估主要參數(shù)并計(jì)算評(píng)估結(jié)果的變化率。最后,對(duì)基于99組虛擬臨床數(shù)據(jù)的求解結(jié)果進(jìn)行統(tǒng)計(jì)學(xué)分析,計(jì)算其平均值及標(biāo)準(zhǔn)差(mean±SD)。

    2 結(jié)果

    2.1 臨床測(cè)量誤差的影響

    圖2表示各項(xiàng)臨床數(shù)據(jù)存在5%的測(cè)量誤差時(shí)主要參數(shù)評(píng)估結(jié)果的變化率。由圖可知,測(cè)量誤差可引起主要參數(shù)評(píng)估結(jié)果的顯著變化。其中,肺動(dòng)脈平均壓(Pmpp)的測(cè)量誤差對(duì)肺循環(huán)阻抗估測(cè)值的影響最大,達(dá)到16.6%;左心房平均壓(Plap)的測(cè)量誤差對(duì)肺循環(huán)阻抗和左心室基準(zhǔn)elastance評(píng)估結(jié)果的影響較大,分別為11.4%和8.8%;大動(dòng)脈平均壓(Pmbp)的測(cè)量誤差對(duì)體循環(huán)阻抗的評(píng)估結(jié)果影響顯著,為9.7%;左心室舒張末期體積(Vedv)的測(cè)量誤差對(duì)各參數(shù)評(píng)估結(jié)果的影響均較明顯,變化率在3.5%~10.3%之間。相對(duì)來(lái)說(shuō),左心室收縮末期體積(Vesv)的測(cè)量誤差影響較小,僅顯著影響左心室收縮期elastance振幅(Elva)的評(píng)估結(jié)果。

    值得注意的是,測(cè)量誤差引起的評(píng)估結(jié)果變化率與所使用的臨床數(shù)據(jù)有關(guān),在特定情況下(如Pmpp和Plap的測(cè)量誤差對(duì)肺循環(huán)阻抗的影響),若標(biāo)準(zhǔn)差大于平均值,則提示測(cè)量誤差對(duì)參數(shù)評(píng)估結(jié)果的影響程度與心血管功能狀態(tài)有關(guān),即存在個(gè)體差異性。

    2.2 次主要參數(shù)變動(dòng)的影響

    圖3顯示了將次主要參數(shù)的值變動(dòng)30%后主要參數(shù)的估測(cè)結(jié)果相對(duì)其參照值的變化率。整體而言,次主要參數(shù)變動(dòng)所引起的評(píng)估結(jié)果變化率較小,最大值為7.4%。體循環(huán)血管順應(yīng)性(Csys)的評(píng)估結(jié)果受右心室功能(Ervb,Erva)的影響相對(duì)明顯,而左心房基準(zhǔn)elastance(Elab)的變動(dòng)對(duì)左心室基準(zhǔn)elastance(Elvb)的評(píng)估結(jié)果也有較為顯著的影響。

    圖2 臨床測(cè)量數(shù)據(jù)存在5%誤差時(shí)主要參數(shù)估測(cè)值的變化率Fig.2 Changes in estimated main parameters as 5% errors are introduced to the clinical data

    圖3 次主要參數(shù)變化30%時(shí)主要參數(shù)估測(cè)值的變化率Fig.3 Changes in estimated main parameters as secondary main parameters are varied by 30%

    3 討論

    基于臨床數(shù)據(jù)對(duì)循環(huán)系統(tǒng)的計(jì)算模型進(jìn)行患者個(gè)體化校準(zhǔn)或評(píng)估患者的心血管功能,已成為心血管生物力學(xué)領(lǐng)域的一個(gè)研究熱點(diǎn)[27-28]。由于可利用的臨床數(shù)據(jù)有限,研究人員通常需要對(duì)模型進(jìn)行簡(jiǎn)化處理,或構(gòu)建僅包含有限數(shù)量模型參數(shù)的目標(biāo)函數(shù),以使參數(shù)優(yōu)化算法收斂于唯一的解[24,29]。在最近的研究中,筆者的研究組發(fā)展了通過(guò)融合臨床數(shù)據(jù)和循環(huán)系統(tǒng)理論模型來(lái)評(píng)估患者心血管功能的方法[15-16]。在上述研究中,參數(shù)優(yōu)化僅針對(duì)特定的模型參數(shù)(即主要參數(shù)),忽略了臨床數(shù)據(jù)的測(cè)量誤差和其他參數(shù)(即次要參數(shù))對(duì)評(píng)估結(jié)果的影響。實(shí)際上,任何臨床數(shù)據(jù)的測(cè)量均無(wú)法避免不確定誤差的存在,而次要參數(shù)所表達(dá)的心血管功能在患者間存在差異,這些因素對(duì)主要參數(shù)評(píng)估結(jié)果的影響有待量化分析,以明確該評(píng)估方法的臨床應(yīng)用價(jià)值或發(fā)現(xiàn)其潛在問(wèn)題。為了研究上述問(wèn)題,筆者進(jìn)一步發(fā)展了針對(duì)測(cè)量誤差和參數(shù)變動(dòng)的誤差分析方法,其間開(kāi)發(fā)了通過(guò)參數(shù)敏感度分析與子集選擇相結(jié)合來(lái)挑選次主要參數(shù)的方法,從而有效降低了計(jì)算負(fù)擔(dān),簡(jiǎn)化了誤差分析過(guò)程。

    本研究的數(shù)值實(shí)驗(yàn)表明,測(cè)量誤差對(duì)主要參數(shù)的評(píng)估結(jié)果有顯著的影響,并且研究結(jié)果與循環(huán)系統(tǒng)血流動(dòng)力學(xué)知識(shí)吻合,部分驗(yàn)證了所開(kāi)發(fā)的誤差分析方法的有效性。例如,肺動(dòng)脈壓和左心房壓共同決定肺循環(huán)的壓力梯度,因此對(duì)肺循環(huán)阻抗的評(píng)估結(jié)果影響顯著,左心房壓和左心室舒張末期體積是反映左心室舒張功能的間接指標(biāo),本研究的結(jié)果亦表明兩者的測(cè)量誤差對(duì)左心室基準(zhǔn)elastance(反映左心室的僵硬度,指示舒張功能)的評(píng)估結(jié)果影響最大?;谏鲜霭l(fā)現(xiàn),筆者認(rèn)為,在實(shí)際應(yīng)用中,對(duì)可能存在較大測(cè)量誤差的臨床數(shù)據(jù),需要在恰當(dāng)評(píng)估測(cè)量誤差影響的基礎(chǔ)上,解釋參數(shù)優(yōu)化結(jié)果。與臨床數(shù)據(jù)的測(cè)量誤差比較,次主要參數(shù)變動(dòng)對(duì)評(píng)估結(jié)果的影響相對(duì)較小,即使變動(dòng)了30%,所引起的最大評(píng)估誤差也在10%以內(nèi)。該發(fā)現(xiàn)表明,雖然次要參數(shù)可能在不同患者間差異明顯,但在主要參數(shù)的評(píng)估計(jì)算中固定這些參數(shù)對(duì)評(píng)估結(jié)果的影響不大,證明課題組開(kāi)發(fā)的以主要參數(shù)為目標(biāo)的參數(shù)評(píng)估方法有較好的魯棒性。

    本研究的局限性主要體現(xiàn)在兩個(gè)方面。首先,基于有限的虛擬臨床數(shù)據(jù)開(kāi)展誤差分析研究,無(wú)法涵蓋所有患者的心血管功能狀態(tài);其次,在進(jìn)行誤差分析時(shí),每組數(shù)值實(shí)驗(yàn)僅針對(duì)某一項(xiàng)臨床數(shù)據(jù)的測(cè)量誤差或某個(gè)模型參數(shù)的變動(dòng),未研究多項(xiàng)臨床數(shù)據(jù)同時(shí)存在誤差或多個(gè)模型參數(shù)同時(shí)變動(dòng)時(shí)評(píng)估結(jié)果的變化規(guī)律。由于人體血液循環(huán)系統(tǒng)呈現(xiàn)高度非線性特征,各心血管功能參數(shù)間相互影響,無(wú)法通過(guò)線性疊加的方法來(lái)拓展現(xiàn)有的研究結(jié)果。為了彌補(bǔ)上述不足,需要在后續(xù)的研究中,在進(jìn)一步擴(kuò)大樣本數(shù)據(jù)的基礎(chǔ)上,引入多項(xiàng)測(cè)量誤差和多參數(shù)變動(dòng),以提高研究結(jié)果的統(tǒng)計(jì)學(xué)意義。

    4 結(jié)論

    為了彌補(bǔ)前期患者個(gè)體化心血管功能評(píng)估研究在參數(shù)依存性和誤差分析方面的不足,本研究進(jìn)一步發(fā)展了次主要參數(shù)的識(shí)別方法,并量化分析了由臨床數(shù)據(jù)的測(cè)量誤差和模型參數(shù)變動(dòng)引起的評(píng)估誤差。研究結(jié)果表明,臨床數(shù)據(jù)的測(cè)量誤差顯著影響評(píng)估結(jié)果,相對(duì)地,評(píng)估結(jié)果受次主要參數(shù)變動(dòng)(或患者個(gè)體差異性)的影響程度較低。因此,確保臨床數(shù)據(jù)的采集精度,將是決定患者個(gè)體化心血管功能評(píng)估的臨床應(yīng)用效果的關(guān)鍵環(huán)節(jié),而本研究的發(fā)現(xiàn)也為解釋評(píng)估結(jié)果和分析誤差提供了理論依據(jù)。在后續(xù)的研究中,將進(jìn)一步考察多項(xiàng)測(cè)量誤差共存和多參數(shù)變動(dòng)對(duì)評(píng)估結(jié)果的綜合影響。

    (致謝 上海交通大學(xué)劉浩教授、上海第九人民醫(yī)院心內(nèi)科范虞琪醫(yī)師為本研究的開(kāi)展和論文撰寫(xiě)提供了寶貴的意見(jiàn)和建議,在此表示衷心的感謝。)

    [1] Evonich RF, Stephens JC, Merhi W, et al. The role of temporary biventricular pacing in the cardiac surgical patient with severely reduced left ventricular systolic function[J]. Journal of Thoracic and Cardiovascular Surgery, 2008, 136(4):915-921.

    [2] Pavlicek M, Wahl A, Rutz T, et al. Right ventricular systolic function assessment: rank of echocardiographic methods vs cardiac magnetic resonance imaging [J]. European Journal of Echocardiography, 2011, 12(11):871-880.

    [3] Vitarelli A, Conde Y, Cimino E, et al. Quantitative assessment of systolic and diastolic ventricular function with tissue Doppler imaging after Fontan type of operation [J]. International Journal of Cardiology, 2005, 102(1):61-69.

    [4] Natsume T, Amano T, Takehara Y, et al. Quantitative assessment of regional systolic and diastolic functions and temporal heterogeneity of myocardial contraction in patients with myocardial infarction using cine magnetic resonance imaging and Fourier fitting [J]. Magnetic Resonance Imaging, 2009, 27(10):1440-1446.

    [5] Vodovar N, Séronde MF, Laribi S, et al. Elevated plasma B-type natriuretic peptide concentrations directly inhibit circulating neprilysin activity in heart failure [J]. JACC: Heart Failure, 2015, 3(8):629-636.

    [6] Senzaki H, Chen C, Masutani S, et al. Assessment of cardiovascular dynamics by pressure-area relations in pediatric patients with congenital heart disease [J]. Journal of Thoracic and Cardiovascular Surgery, 2001, 122(3):535-547.

    [7] Troughton RW, Prior DL, Pereira JJ, et al. Plasma B-type natriuretic peptide levels in systolic heart failure: importance of left ventricular diastolic function and right ventricular systolic function [J]. J Am Coll Cardiol, 2004, 43(3):416-422.

    [8] Suga H, Sagawa K, Shoukas AA. Load independence of the instantaneous pressure-volume ratio of the canine left ventricle and effects of epinephrine and heart rate on the ratio [J]. Circulation Research, 1973,32:314-322.

    [9] Senzaki H, Chen CH, Kass DA. Single-beat estimation of end-systolic pressure-volume relation in humans: A new method with the potential for noninvasive application [J]. Circulation, 1996, 94(10):2497-2506.

    [10] Xiao X, Ozawa ET, Huang Y, et al. Model-based assessment of cardiovascular health from noninvasive measurements [J]. Ann Biomed Eng, 2002, 30(5):612-623.

    [11] Hann CE, Revie J, Stevenson D, et al. Patient specific identification of the cardiac driver function in a cardiovascular system model[J]. Computer Methods and Programs in Biomedicine, 2011, 101(2):201-207.

    [12] Revie J, Stevenson D, Chase JG, et al. Validation of subject-specific cardiovascular system models from porcine measurements[J]. Computer Methods and Programs in Biomedicine, 2013, 109(2):197-210.

    [13] Chase JG, Lambermont B, Starfinger C, et al. Subject-specific cardiovascular system model-based identification and diagnosis of septic shock with a minimally invasive data set: animal experiments and proof of concept[J]. Physiological Measurement, 2011, 32(1):65-82.

    [14] Yu YC, Boston JR, Simaan MA, et al. Minimally invasive estimation of systemic vascular parameters[J]. Ann Biomed Eng, 2001, 29(7):595-606.

    [15] Sughimoto K, Liang F, Takahara Y, et al. Assessment of cardiovascular function by combining clinical data with a computational model of the cardiovascular system[J]. Journal of Thoracic and Cardiovascular Surgery, 2013, 145(5):1367-1372.

    [16] Liang F, Sughimoto K, Matsuo K, et al. Patient-specific assessment of cardiovascular function by combination of clinical data and computational model with applications to patients undergoing Fontan operation[J]. International Journal for Numerical Methods in Biomedical Engineering, 2014, 30(10):1000-1018.

    [17] Liang F, Senzaki H, Kurishima C, et al. Hemodynamic performance of the Fontan circulation compared with a normal biventricular circulation: a computational model study[J]. American Journal of Physiology-Heart and Circulatory Physiology, 2014,307(7), H1056-H1072.

    [18] Liang F, Oshima M, Huang H, et al. Numerical study of cerebroarterial hemodynamic changes following carotid artery operation: A comparison between multiscale modeling and stand-alone three-dimensional modeling [J]. Journal of Biomechanical Engineering, 2015, 137(10):101011.

    [19] 樊瑜波,陳君楷,袁支潤(rùn),等. 全分支的體循環(huán)混合參數(shù)計(jì)算機(jī)模擬研究[J]. 生物醫(yī)學(xué)工程學(xué)雜志,1995,12(3):128-223.

    [20] 丁光宏,魏崗之. 腦循環(huán)脈動(dòng)流的集中參數(shù)模型[J]. 力學(xué)學(xué)報(bào),1996,28(3):336-346.

    [21] 王盛章,姚偉,丁光宏. 腦血流自動(dòng)調(diào)節(jié)功能的三單元模型[J]. 醫(yī)用生物力學(xué),2007,22(3):268-272.

    [22] 谷凱云,高斌,劉有軍,等. 基于主動(dòng)脈內(nèi)血泵的集中參數(shù)模型研究[J]. 醫(yī)用生物力學(xué), 2011, 26(4):367-372.

    [23] Van Griensven A, Meixner T, Grunwald S, et al. A global sensitivity analysis tool for the parameters of multi-variable catchment models [J]. Journal of Hydrology, 2006, 324(1):10-23.

    [24] Olufsen MS, Ottesen JT. A practical approach to parameter estimation applied to model predicting heart rate regulation[J]. Journal of Mathematical Biology, 2013, 67(1):39-68.

    [25] Cintron-Arias A, Banks HT, Capaldi A, et al. A sensitivity matrix based methodology for inverse problem formulation[J]. Journal of Inverse and Ill-posed Problems, 2009, 17(6):545-564.

    [26] Yue H, Brown M, Knowles J, et al. Insights into the behavior of systems biology models from dynamic sensitivity and identifiability analysis: a case study of an NF-kB signaling pathway[J]. Mol Biosyst, 2006, 2(12):640-649.

    [27] Pennati G, Corsini C, Cosentino D, et al. Boundary conditions of patient-specific fluid dynamics modeling of cavopulmonary connections: possible adaptation of pulmonary resistances results in a critical issue for a virtual surgical planning [J]. Interface Focus, 2011, 1(3):297-307.

    [28] Sun Qi, Liu Jinlong, Qian Yi, et al. Computational haemodynamic analysis of patient-specific virtual operations for total cavopulmonary connection with dual superior venae cavae [J]. European Journal of Cardio-Thoracic Surgery, 2014, 45(3):564-569.

    [29] Pope SR, Ellwein LM, Zapata CL, et al. Estimation and identification of parameters in a lumped cerebrovascular model[J]. Math Biosci Eng, 2009, 6(1):93-115.

    Error Analysis on the Assessment of Cardiovascular Function Based on Integration of Clinical Data and Cardiovascular Model

    Li Yi1Yin Zhaofang2Liang Fuyou1#*

    1(SchoolofNavalArchitecture,OceanandCivilEngineering,ShanghaiJiaoTongUniversity,Shanghai200240,China)2(DepartmentofCardiology,ShanghaiNinthPeople’sHospital,ShanghaiJiaoTongUniversity,Shanghai200011,China)

    The purpose of this study is to identify the key factors that affect the accuracy of patient-specific assessment of cardiovascular function (by means of integration of clinical data and cardiovascular model) and quantify their effects so as to provide theoretical evidence for guiding the clinical application of the assessment method. Parameter sensitivity analysis was performed in combination with parameter subset selection to identify the secondary main parameters that are related closely to the model-based prediction of hemodynamic variables and the assessment of the main parameters (corresponding to the assessed cardiovascular function). Numerical experiments were carried out based on a series of virtual clinical data to quantify the changes in assessment results induced by measurement errors (in a range of 5%) of clinical data and variations in the secondary main parameters (rate of change being 30%). Measurement errors of clinical data were found to induce pronounced changes in assessment results (up to 16.6%), relatively, variations in the secondary main parameters had less influence on the assessment (rates of changes being within 10%). Accurate clinical data measurement is a key step to guaranteeing the reliability of cardiovascular function assessment. The secondary main parameters only have limited influence on the assessment although they may vary significantly among patients.

    Assessment of cardiovascular function; cardiovascular model; sensitivity analysis; subset selection; error analysis

    10.3969/j.issn.0258-8021. 2016. 01.006

    2015-09-13, 錄用日期:2015-11-29

    國(guó)家自然科學(xué)基金項(xiàng)目(81370438);上海交通大學(xué)醫(yī)工交叉研究基金項(xiàng)目(YG2012MS24)

    R318

    A

    0258-8021(2016) 01-0047-08

    # 中國(guó)生物醫(yī)學(xué)工程學(xué)會(huì)會(huì)員(Member, Chinese Society of Biomedical Engineering)

    *通信作者(Corresponding author), E-mail: fuyouliang@sjtu.edu.cn

    猜你喜歡
    主要參數(shù)測(cè)量誤差左心室
    密度測(cè)量誤差分析
    心電向量圖診斷高血壓病左心室異常的臨床應(yīng)用
    生活垃圾分類對(duì)垃圾主要參數(shù)的影響分析
    有機(jī)硅流化床氣體分布板主要參數(shù)設(shè)計(jì)
    縱向數(shù)據(jù)下變系數(shù)測(cè)量誤差模型的漸近估計(jì)
    碎邊剪剪切特性分析與主要參數(shù)確定
    初診狼瘡腎炎患者左心室肥厚的相關(guān)因素
    牽引變壓器功率測(cè)量誤差分析
    影響輪軌粘滑振動(dòng)的主要參數(shù)分析
    IMU/GPS測(cè)量誤差對(duì)斜視條件下機(jī)載重軌干涉
    午夜爱爱视频在线播放| 欧美97在线视频| 色5月婷婷丁香| 日本免费在线观看一区| 亚洲真实伦在线观看| 五月伊人婷婷丁香| 在线观看三级黄色| 久久久久久久精品精品| 97精品久久久久久久久久精品| 尤物成人国产欧美一区二区三区| www.av在线官网国产| 日韩av在线免费看完整版不卡| 国产精品久久久久久精品电影| 久久ye,这里只有精品| 一本久久精品| 中文字幕人妻熟人妻熟丝袜美| 日韩,欧美,国产一区二区三区| 91午夜精品亚洲一区二区三区| av又黄又爽大尺度在线免费看| 69av精品久久久久久| 在现免费观看毛片| 日日啪夜夜爽| 99精国产麻豆久久婷婷| 最新中文字幕久久久久| 免费大片18禁| 亚洲精品第二区| 新久久久久国产一级毛片| 免费看光身美女| 亚洲精品aⅴ在线观看| 高清日韩中文字幕在线| 成人无遮挡网站| 欧美3d第一页| 国产一区二区三区综合在线观看 | 免费播放大片免费观看视频在线观看| 少妇裸体淫交视频免费看高清| 免费在线观看成人毛片| 欧美高清成人免费视频www| 寂寞人妻少妇视频99o| 久久久成人免费电影| 久久6这里有精品| 美女高潮的动态| 另类亚洲欧美激情| 精品人妻偷拍中文字幕| 亚洲欧洲日产国产| 水蜜桃什么品种好| 国产欧美日韩精品一区二区| 男的添女的下面高潮视频| 免费播放大片免费观看视频在线观看| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 国产毛片在线视频| 日韩欧美精品v在线| 亚洲最大成人av| 黄色视频在线播放观看不卡| 亚洲欧美成人精品一区二区| 黄片无遮挡物在线观看| 欧美一区二区亚洲| 久久精品夜色国产| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 国产亚洲午夜精品一区二区久久 | 久久精品人妻少妇| 激情五月婷婷亚洲| 熟女人妻精品中文字幕| 久久精品熟女亚洲av麻豆精品| 日本免费在线观看一区| 国产精品人妻久久久影院| 国语对白做爰xxxⅹ性视频网站| 99久久人妻综合| 亚洲成色77777| 亚洲欧美成人综合另类久久久| 在线观看美女被高潮喷水网站| 99热这里只有是精品在线观看| 人人妻人人看人人澡| 久久精品人妻少妇| 99九九线精品视频在线观看视频| 色视频www国产| 国产伦精品一区二区三区四那| av.在线天堂| 3wmmmm亚洲av在线观看| 一个人观看的视频www高清免费观看| 啦啦啦啦在线视频资源| 久久精品国产亚洲av涩爱| 亚洲欧美精品专区久久| 一区二区三区四区激情视频| 国产乱人偷精品视频| av天堂中文字幕网| 亚洲av免费在线观看| 免费观看的影片在线观看| 精品亚洲乱码少妇综合久久| 欧美日韩综合久久久久久| 91午夜精品亚洲一区二区三区| 国产亚洲5aaaaa淫片| 自拍偷自拍亚洲精品老妇| 亚洲精品日韩在线中文字幕| 最新中文字幕久久久久| 色哟哟·www| 一区二区三区四区激情视频| 午夜爱爱视频在线播放| 狠狠精品人妻久久久久久综合| 欧美潮喷喷水| 日韩,欧美,国产一区二区三区| 日韩欧美 国产精品| 黄片wwwwww| 国产av不卡久久| 一区二区三区免费毛片| 少妇熟女欧美另类| 亚洲av成人精品一区久久| 国产精品.久久久| 日韩欧美 国产精品| 高清av免费在线| 亚洲精品日本国产第一区| 精品久久久噜噜| 欧美日韩综合久久久久久| 各种免费的搞黄视频| 久久久久网色| 久久久久久久久久成人| 制服丝袜香蕉在线| 日韩大片免费观看网站| 插阴视频在线观看视频| 免费黄频网站在线观看国产| 日韩av不卡免费在线播放| 国产高清国产精品国产三级 | 成人美女网站在线观看视频| 日产精品乱码卡一卡2卡三| 亚洲av一区综合| 女的被弄到高潮叫床怎么办| 久久久久久久久久人人人人人人| 少妇的逼好多水| 一级黄片播放器| 色5月婷婷丁香| 少妇人妻一区二区三区视频| 国模一区二区三区四区视频| 在线观看一区二区三区激情| 亚洲精品久久久久久婷婷小说| 综合色丁香网| 亚洲图色成人| 亚洲精品乱码久久久久久按摩| 国产精品精品国产色婷婷| 成年免费大片在线观看| 内射极品少妇av片p| 日本色播在线视频| 王馨瑶露胸无遮挡在线观看| 精品一区二区三卡| 国产黄频视频在线观看| 国产精品久久久久久av不卡| 可以在线观看毛片的网站| freevideosex欧美| 69av精品久久久久久| 老女人水多毛片| 大陆偷拍与自拍| 亚州av有码| 男人爽女人下面视频在线观看| 小蜜桃在线观看免费完整版高清| 亚洲美女视频黄频| 国内精品宾馆在线| 看黄色毛片网站| 国产一区二区三区av在线| 在线观看美女被高潮喷水网站| 欧美成人精品欧美一级黄| 黄色日韩在线| av一本久久久久| 成人二区视频| 1000部很黄的大片| 亚洲最大成人av| 极品教师在线视频| 三级经典国产精品| 国产男人的电影天堂91| 国产精品一二三区在线看| 日韩国内少妇激情av| 少妇丰满av| 国产成人a∨麻豆精品| tube8黄色片| 国产爱豆传媒在线观看| 色视频在线一区二区三区| 97在线人人人人妻| 天天躁夜夜躁狠狠久久av| 亚洲av电影在线观看一区二区三区 | 直男gayav资源| 久久久欧美国产精品| 人妻制服诱惑在线中文字幕| 日韩欧美 国产精品| 久久ye,这里只有精品| 国产精品蜜桃在线观看| 大话2 男鬼变身卡| 青春草亚洲视频在线观看| 一级av片app| 黄色日韩在线| 搡老乐熟女国产| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 亚洲av.av天堂| 午夜爱爱视频在线播放| 精品酒店卫生间| 日韩三级伦理在线观看| 97精品久久久久久久久久精品| 伦理电影大哥的女人| 大码成人一级视频| 青春草国产在线视频| 国产精品一区二区三区四区免费观看| 亚州av有码| 久久精品久久久久久久性| 日韩国内少妇激情av| 久久久久久伊人网av| 丝袜脚勾引网站| 欧美+日韩+精品| 精品久久久久久电影网| 老司机影院成人| 亚洲精华国产精华液的使用体验| 亚洲美女视频黄频| 在线免费十八禁| 伦精品一区二区三区| 久久久久精品久久久久真实原创| 午夜爱爱视频在线播放| 成人亚洲欧美一区二区av| 97在线视频观看| 一级av片app| 精品国产露脸久久av麻豆| 国产黄a三级三级三级人| 人妻系列 视频| 成人高潮视频无遮挡免费网站| 我的女老师完整版在线观看| 超碰97精品在线观看| 国产男女超爽视频在线观看| 国产男女内射视频| 亚洲人成网站在线播| 午夜福利在线观看免费完整高清在| 亚洲最大成人av| 在线亚洲精品国产二区图片欧美 | 色综合色国产| 精品人妻偷拍中文字幕| 又粗又硬又长又爽又黄的视频| 美女cb高潮喷水在线观看| 国产大屁股一区二区在线视频| 老司机影院毛片| 久久久久精品久久久久真实原创| 少妇猛男粗大的猛烈进出视频 | 久久国产乱子免费精品| 激情 狠狠 欧美| 久久精品国产亚洲网站| 欧美变态另类bdsm刘玥| 大话2 男鬼变身卡| 精品午夜福利在线看| 日本欧美国产在线视频| 国产亚洲av片在线观看秒播厂| 亚洲欧美日韩东京热| 搡女人真爽免费视频火全软件| 久久久色成人| 三级经典国产精品| 国产高清国产精品国产三级 | 成人免费观看视频高清| 五月玫瑰六月丁香| 中文字幕制服av| 永久网站在线| 性插视频无遮挡在线免费观看| 亚洲真实伦在线观看| 69av精品久久久久久| 国产免费又黄又爽又色| 蜜臀久久99精品久久宅男| 国产免费一区二区三区四区乱码| 看免费成人av毛片| 搡老乐熟女国产| 别揉我奶头 嗯啊视频| 国产亚洲5aaaaa淫片| 亚洲最大成人手机在线| 色综合色国产| 欧美另类一区| 男女那种视频在线观看| eeuss影院久久| 老司机影院成人| 内射极品少妇av片p| av又黄又爽大尺度在线免费看| 成人鲁丝片一二三区免费| 久久久久久久久久人人人人人人| 日本av手机在线免费观看| 免费人成在线观看视频色| 男女啪啪激烈高潮av片| 天天躁日日操中文字幕| 欧美日本视频| 精品午夜福利在线看| 国产一级毛片在线| 制服丝袜香蕉在线| 久久久久久久午夜电影| 国产高清三级在线| 日韩中字成人| 免费看av在线观看网站| 国产免费一级a男人的天堂| 18禁在线无遮挡免费观看视频| 成人午夜精彩视频在线观看| 精品午夜福利在线看| 国产精品蜜桃在线观看| av在线观看视频网站免费| av在线播放精品| 国产欧美日韩精品一区二区| 午夜激情久久久久久久| 免费黄网站久久成人精品| 男的添女的下面高潮视频| 卡戴珊不雅视频在线播放| 18+在线观看网站| 亚洲精品乱久久久久久| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 久久精品国产亚洲网站| 身体一侧抽搐| 免费观看无遮挡的男女| 色5月婷婷丁香| 国产av不卡久久| 别揉我奶头 嗯啊视频| 日韩欧美精品v在线| 麻豆成人午夜福利视频| 国产高清不卡午夜福利| 综合色av麻豆| 可以在线观看毛片的网站| 久久久成人免费电影| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人 | 国产精品精品国产色婷婷| 新久久久久国产一级毛片| 91aial.com中文字幕在线观看| 毛片女人毛片| 九草在线视频观看| 一级二级三级毛片免费看| 新久久久久国产一级毛片| 免费观看a级毛片全部| 秋霞在线观看毛片| 亚洲成人av在线免费| 在线观看国产h片| 亚洲在线观看片| 久久精品国产自在天天线| 国产视频内射| 蜜臀久久99精品久久宅男| 色婷婷久久久亚洲欧美| 一级毛片黄色毛片免费观看视频| 一个人观看的视频www高清免费观看| 国产精品人妻久久久久久| 亚洲精品色激情综合| 久久精品国产a三级三级三级| 哪个播放器可以免费观看大片| 韩国av在线不卡| 欧美三级亚洲精品| 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品| 日日撸夜夜添| 亚洲不卡免费看| 免费看a级黄色片| 日韩免费高清中文字幕av| 亚洲成色77777| 美女xxoo啪啪120秒动态图| 国产精品麻豆人妻色哟哟久久| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 少妇人妻久久综合中文| 精品人妻一区二区三区麻豆| 三级国产精品欧美在线观看| 中国三级夫妇交换| 亚洲高清免费不卡视频| 久久精品国产a三级三级三级| 一级爰片在线观看| 三级国产精品欧美在线观看| 97在线视频观看| 欧美高清性xxxxhd video| 最近中文字幕2019免费版| 国产一区二区亚洲精品在线观看| 如何舔出高潮| 大陆偷拍与自拍| 激情 狠狠 欧美| 老师上课跳d突然被开到最大视频| 亚洲精品乱码久久久久久按摩| 18禁在线无遮挡免费观看视频| 人人妻人人看人人澡| 99久久中文字幕三级久久日本| 成人鲁丝片一二三区免费| 大话2 男鬼变身卡| 欧美激情在线99| 人妻 亚洲 视频| 国产成人91sexporn| 美女高潮的动态| 91久久精品电影网| 91精品伊人久久大香线蕉| 日韩视频在线欧美| 亚洲欧美日韩无卡精品| 国产男女内射视频| 亚洲欧洲日产国产| 亚洲国产精品国产精品| 国产真实伦视频高清在线观看| 在线免费十八禁| 大陆偷拍与自拍| 听说在线观看完整版免费高清| 日日啪夜夜爽| 婷婷色综合大香蕉| 最近2019中文字幕mv第一页| 在线a可以看的网站| 久久久久久久精品精品| tube8黄色片| 日韩不卡一区二区三区视频在线| 黄色视频在线播放观看不卡| 看十八女毛片水多多多| 美女xxoo啪啪120秒动态图| 精品久久久久久久人妻蜜臀av| 各种免费的搞黄视频| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 人妻少妇偷人精品九色| 免费观看av网站的网址| 亚洲婷婷狠狠爱综合网| 内射极品少妇av片p| 2018国产大陆天天弄谢| 亚洲美女搞黄在线观看| 日韩中字成人| 国产亚洲一区二区精品| 国内少妇人妻偷人精品xxx网站| 丰满少妇做爰视频| 女人久久www免费人成看片| 久久精品久久精品一区二区三区| 国产伦精品一区二区三区四那| 97超视频在线观看视频| av线在线观看网站| 日本wwww免费看| 久久久久久久大尺度免费视频| 身体一侧抽搐| 精品久久久久久久末码| 亚洲最大成人中文| 国产高清国产精品国产三级 | 久久久久久久国产电影| 亚洲天堂av无毛| 午夜激情福利司机影院| 三级国产精品欧美在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲性久久影院| 97在线视频观看| 97在线人人人人妻| 久久99蜜桃精品久久| 日本一二三区视频观看| 国产精品.久久久| 美女xxoo啪啪120秒动态图| av专区在线播放| 亚洲精品色激情综合| 久久精品国产a三级三级三级| 欧美日韩在线观看h| 熟女人妻精品中文字幕| 免费观看在线日韩| 午夜福利视频精品| 日日摸夜夜添夜夜添av毛片| 大片电影免费在线观看免费| 国产在线一区二区三区精| 国产午夜精品久久久久久一区二区三区| 午夜爱爱视频在线播放| 欧美日韩综合久久久久久| 丝袜美腿在线中文| 欧美激情国产日韩精品一区| 在线观看国产h片| 亚洲人成网站在线观看播放| 高清av免费在线| 免费少妇av软件| 一区二区av电影网| 中文在线观看免费www的网站| 少妇的逼好多水| 另类亚洲欧美激情| 一二三四中文在线观看免费高清| 你懂的网址亚洲精品在线观看| 黄色视频在线播放观看不卡| av免费观看日本| 中文在线观看免费www的网站| 日本猛色少妇xxxxx猛交久久| 中国美白少妇内射xxxbb| 亚洲国产精品成人综合色| av.在线天堂| 少妇人妻精品综合一区二区| 少妇丰满av| 国产精品麻豆人妻色哟哟久久| 国产视频内射| 日韩欧美精品v在线| 久久精品国产鲁丝片午夜精品| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲人成网站高清观看| 一级毛片电影观看| 国产成人精品福利久久| 亚洲精品日本国产第一区| 2021天堂中文幕一二区在线观| 亚洲美女搞黄在线观看| 国产男女内射视频| 久久鲁丝午夜福利片| 男插女下体视频免费在线播放| 成人毛片a级毛片在线播放| 啦啦啦啦在线视频资源| 我的女老师完整版在线观看| 人人妻人人看人人澡| 91久久精品国产一区二区成人| 寂寞人妻少妇视频99o| 韩国高清视频一区二区三区| 日韩制服骚丝袜av| 成人欧美大片| 最近中文字幕高清免费大全6| 全区人妻精品视频| 少妇 在线观看| 在线观看一区二区三区| 大又大粗又爽又黄少妇毛片口| av在线天堂中文字幕| 三级经典国产精品| 亚洲欧美一区二区三区国产| 最新中文字幕久久久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久午夜乱码| 永久网站在线| 嘟嘟电影网在线观看| 亚洲av.av天堂| 看黄色毛片网站| 国产成人a∨麻豆精品| 日韩欧美 国产精品| 香蕉精品网在线| 最近最新中文字幕免费大全7| 97超碰精品成人国产| 国产免费一区二区三区四区乱码| 亚洲精品乱码久久久久久按摩| 免费黄网站久久成人精品| 少妇的逼好多水| 人妻制服诱惑在线中文字幕| 亚洲精品日本国产第一区| 亚洲精品日韩av片在线观看| 欧美一级a爱片免费观看看| 18禁在线播放成人免费| 你懂的网址亚洲精品在线观看| 麻豆精品久久久久久蜜桃| 精品久久久久久电影网| 欧美高清成人免费视频www| 王馨瑶露胸无遮挡在线观看| 婷婷色麻豆天堂久久| 亚洲图色成人| 亚洲最大成人av| 亚洲成人av在线免费| 韩国av在线不卡| 在线观看一区二区三区| av线在线观看网站| 99久久中文字幕三级久久日本| 嫩草影院入口| 亚洲av成人精品一区久久| 亚洲精品乱久久久久久| 又黄又爽又刺激的免费视频.| 亚洲av男天堂| 日本一本二区三区精品| 国产男人的电影天堂91| 免费人成在线观看视频色| 精品久久国产蜜桃| 人人妻人人看人人澡| 一级毛片久久久久久久久女| 免费看光身美女| 毛片一级片免费看久久久久| 欧美日韩综合久久久久久| 亚洲四区av| 国产爱豆传媒在线观看| 建设人人有责人人尽责人人享有的 | 亚洲av免费高清在线观看| 久久热精品热| 美女高潮的动态| 男女下面进入的视频免费午夜| 涩涩av久久男人的天堂| 国产成人精品婷婷| 免费av观看视频| 国产高清有码在线观看视频| 在线a可以看的网站| 老司机影院成人| 蜜臀久久99精品久久宅男| 亚洲国产成人一精品久久久| 我的女老师完整版在线观看| 久久这里有精品视频免费| 欧美一区二区亚洲| 亚洲av二区三区四区| 国产亚洲av嫩草精品影院| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 少妇丰满av| 久久久欧美国产精品| 男插女下体视频免费在线播放| 亚洲熟女精品中文字幕| 国产欧美日韩精品一区二区| av国产久精品久网站免费入址| 自拍偷自拍亚洲精品老妇| 国产成人a区在线观看| 国产探花在线观看一区二区| 91精品伊人久久大香线蕉| 亚洲精品久久久久久婷婷小说| tube8黄色片| 女人十人毛片免费观看3o分钟| 男女边摸边吃奶| 国产 一区 欧美 日韩| 最近中文字幕高清免费大全6| 国产爱豆传媒在线观看| 国产日韩欧美亚洲二区| av在线天堂中文字幕| 成人亚洲欧美一区二区av| 亚洲人与动物交配视频| 国产精品无大码| 99久久人妻综合| 一级片'在线观看视频| kizo精华| 91久久精品国产一区二区成人| 国产av码专区亚洲av| 久热久热在线精品观看| 久久久久久久午夜电影| 亚洲国产高清在线一区二区三| 中文字幕制服av| 我要看日韩黄色一级片| 人妻夜夜爽99麻豆av| 亚洲国产成人一精品久久久| 久久精品综合一区二区三区| 亚洲国产高清在线一区二区三| 午夜精品国产一区二区电影 | 国产成人精品一,二区| 91精品伊人久久大香线蕉| 久久99蜜桃精品久久| 成人一区二区视频在线观看| 一区二区三区免费毛片| 日本色播在线视频|