張鵬鵬，高雪華，王青敏，呂 振，張新宇

醫(yī)學(xué)超聲關(guān)鍵技術(shù)國(guó)家地方聯(lián)合工程實(shí)驗(yàn)室, 廣東省生物醫(yī)學(xué)信息檢測(cè)與超聲成像重點(diǎn)實(shí)驗(yàn)室,深圳大學(xué)醫(yī)學(xué)部，深圳 518060

【電子與信息科學(xué) / Electronics and Information Science】

角膜黏彈性定量測(cè)量評(píng)估膠原交聯(lián)處理效果

張鵬鵬，高雪華，王青敏，呂 振，張新宇

醫(yī)學(xué)超聲關(guān)鍵技術(shù)國(guó)家地方聯(lián)合工程實(shí)驗(yàn)室, 廣東省生物醫(yī)學(xué)信息檢測(cè)與超聲成像重點(diǎn)實(shí)驗(yàn)室,深圳大學(xué)醫(yī)學(xué)部，深圳 518060

建立基于VerasonicsTMVantage256 超聲開放系統(tǒng)的實(shí)驗(yàn)平臺(tái)，用聲輻射力激勵(lì)角膜振動(dòng)，通過快速平面波成像檢測(cè)A0模式蘭姆波的傳播，根據(jù)蘭姆波的頻散特性估計(jì)牛眼角膜的剪切彈性模量和剪切黏性系數(shù)．通過對(duì)牛眼角膜進(jìn)行離體實(shí)驗(yàn)研究，評(píng)估紫外線-核黃素角膜膠原交聯(lián)術(shù)(corneal collagen cross-linking, CXL)中加入戊二醛和不同紫外線照射時(shí)間的影響．結(jié)果表明，加入質(zhì)量分?jǐn)?shù)為0.1%的戊二醛對(duì)角膜的硬化效果較未加入的更強(qiáng)．加入質(zhì)量分?jǐn)?shù)為0.1%的戊二醛的CXL處理組，紫外線照射5、10、30和60 min對(duì)應(yīng)的楊氏模量比直接CXL處理組的彈性值分別增大了148.78%、335.68%、392.59%和250.47%．

生物醫(yī)學(xué)測(cè)量學(xué)；醫(yī)學(xué)超聲；超聲彈性成像；聲輻射力；角膜；黏彈性；膠原交聯(lián)處理

角膜是眼睛最前面的透明部分，覆蓋虹膜、瞳孔及前房．角膜從外到內(nèi)可分為上皮細(xì)胞層、前彈力層、基質(zhì)層、后彈力層和內(nèi)皮細(xì)胞層，其光學(xué)特性和力學(xué)特性對(duì)維持眼球功能起著重要作用．圓錐角膜疾病，是以角膜中央變薄向前突出，呈圓錐形特征的一種眼病，包括原發(fā)性圓錐角膜和角膜屈光手術(shù)后繼發(fā)性角膜膨?。渲?，以原發(fā)性圓錐角膜為典型代表，發(fā)病率達(dá)0.05%[1]，常造成高度不規(guī)則近視散光和不同程度的視力損害，現(xiàn)是角膜移植的第二常見指征．角膜膠原交聯(lián)術(shù)(corneal collagen cross-linking, CXL)[2]是目前唯一確認(rèn)阻止圓錐角膜病變進(jìn)展的治療手段．CXL經(jīng)過十余年發(fā)展，在歐洲已成為臨床治療圓錐角膜的常規(guī)手段，中國(guó)還在臨床推廣中．中國(guó)作為角膜植片嚴(yán)重匱乏、廣泛大量開展角膜屈光手術(shù)、嚴(yán)重感染性角膜炎仍是主要致盲因素之一的國(guó)家，倡導(dǎo)并推動(dòng)該技術(shù)在眼科的研究和臨床應(yīng)用，具有重要社會(huì)意義[3]．角膜生物力學(xué)特性的在體定量測(cè)量研究可用于圓錐角膜疾病的診斷及CXL術(shù)后效果評(píng)估，為CXL技術(shù)的應(yīng)用和推廣提供了技術(shù)保障．

早期角膜力學(xué)特性測(cè)量主要采用拉伸實(shí)驗(yàn)[4]及眼球膨脹測(cè)試[5]等離體實(shí)驗(yàn)研究．目前美國(guó)Reichert公司的眼反應(yīng)分析儀 (ocular response analyzer，ORA)和德國(guó)Oculus公司推出的Corvis ST角膜生物力學(xué)分析儀都已用于臨床，提供與角膜整體力學(xué)特性有關(guān)的參數(shù)．對(duì)角膜的力學(xué)特性進(jìn)行高分辨率成像的方法也有較快的進(jìn)展，如電子散斑干涉法[6]、布里淵光學(xué)顯微鏡[7]、光學(xué)相干彈性成像法[8]和超聲彈性成像法[9]．這些方法的原理和實(shí)現(xiàn)方式各不相同，目前基本處在離體實(shí)驗(yàn)階段．其中，電子散斑干涉法和光學(xué)相干彈性成像需通過外加機(jī)械力使角膜振動(dòng)，再分別用散斑成像技術(shù)或光學(xué)相干斷層掃描技術(shù)(optical coherence tomography, OCT)技術(shù)追蹤角膜的振動(dòng)，估計(jì)其力學(xué)特性；超聲彈性成像可用同一探頭激勵(lì)角膜振動(dòng)并檢測(cè)，但探頭和角膜間需要水作為耦合劑．布里淵顯微鏡法根據(jù)焦點(diǎn)處的布里淵頻移與硬度的相關(guān)性評(píng)估組織硬度，但是完成一次測(cè)量時(shí)間較長(zhǎng)．Tanter等[9]將剪切波彈性成像(supersonic shear imaging,SSI)技術(shù)應(yīng)用于角膜，根據(jù)簡(jiǎn)化模型估計(jì)角膜的楊氏模量．Zhang等[10]提出基于蘭姆波頻散特性的超聲彈性成像方法定量估計(jì)角膜的黏彈性特性．

本研究建立了基于VerasonicsTMVantage256 超聲開放系統(tǒng)，實(shí)現(xiàn)了角膜黏彈性估計(jì)算法，并設(shè)計(jì)了兩種CXL方案處理角膜，使用所建立的系統(tǒng)平臺(tái)分別測(cè)量?jī)煞N方案處理后角膜的黏彈性變化，定量評(píng)估兩種CXL方案處理后角膜的硬化效果．

1 方法原理

1.1 原 理

本研究基于蘭姆波在角膜組織中傳播的原理，使用超聲方法激勵(lì)并檢測(cè)蘭姆波的傳播特性，估算出角膜的黏彈性．彈性波在無限長(zhǎng)薄板中的傳播模式被命名為蘭姆波，被廣泛用于導(dǎo)波在圓柱殼內(nèi)傳播的研究中[11]．Kanai[12]將漏蘭姆波的頻散特性方程進(jìn)行近似變形，用于估計(jì)心臟血管壁的黏彈性．本研究將超聲輻射力觸發(fā)角膜產(chǎn)生的彈性波近似為蘭姆波．蘭姆波的傳播模式具有高頻散的特性，在低頻范圍(對(duì)角膜低于2 kHz)內(nèi)僅2種模式：零階對(duì)稱模式(zero-order symmetric mode，S0)與零階反對(duì)稱模式(zero-order asymmetric mode，A0)，其中，A0模式振幅大且呈現(xiàn)出明顯的頻散特性．Kanai[12]假設(shè)心室壁為厚度2h的黏彈性薄板且密度與周圍液體相近，并假設(shè)心室壁中縱波的波數(shù)遠(yuǎn)小于蘭姆波的波數(shù)，對(duì)漏蘭姆波的頻散特征方程進(jìn)行近似化簡(jiǎn)，提出蘭姆波在心室壁中傳播的頻散特性方程為

(1)

其中，ω為角頻率；cL(ω)和cS(ω)分別為角頻率為ω時(shí)的A0蘭姆波和對(duì)應(yīng)剪切波的相速度；KL=ω/cL(ω)和Ks=ω/cs(ω)分別為蘭姆波和剪切波的波數(shù)，且

(2)

角膜同樣滿足Kanai提出的假設(shè)．因此本研究采用方程(1)描述角膜內(nèi)A0蘭姆波傳播特性，并引入Voigt模型(由彈簧和阻尼器并聯(lián)組成)描述角膜的黏彈性，則有

μ=μ1+iωμ2

(3)

其中，μ1和μ2分別為剪切彈性和剪切黏性．對(duì)于Voigt材料，剪切波波數(shù)Ks取決于μ1和μ2， 如

(4)

其中，ρ是角膜組織的密度．

將式(4)代入方程(1)可得到蘭姆波速度、頻率、角膜厚度2h、 及μ1和μ2的方程式

f(cL，μ1，μ2，ω，h)=0

(5)

取h等于實(shí)測(cè)角膜厚度的1/2．設(shè)角膜密度ρ≈1.0×103kg/m3，根據(jù)實(shí)測(cè)的蘭姆波相速度，采用非線性優(yōu)化方法解方程(5)的逆問題，得到μ1和μ2的估計(jì)值．算法的基本思想是遍歷<μ1,μ2>的組合，對(duì)于每一組<μ1,μ2>， 代入式(5)得到各頻率下蘭姆波相速度的理論值,找到使估計(jì)的相速度與實(shí)測(cè)相速度cL(ω)之間的誤差平方和最小的<μ1,μ2>組合，作為μ1和μ2的估計(jì)．即

(6)

1.2 超聲測(cè)量平臺(tái)搭建

基于VerasonicsTMVantage256開放性可編程超聲系統(tǒng)構(gòu)建如圖1的實(shí)驗(yàn)平臺(tái)．Vantage256系統(tǒng)可通過Matlab環(huán)境下的腳本程序，獨(dú)立控制各通道發(fā)射和接收，實(shí)時(shí)訪問各通道射頻信號(hào)(radiofrequency，RF)數(shù)據(jù)，并對(duì)數(shù)據(jù)進(jìn)行處理和分析．平臺(tái)采用128陣元的線陣探頭(L11-4)，對(duì)眼球進(jìn)行B-mode超聲成像及超聲彈性測(cè)量．B-mode圖像用于彈性測(cè)量的引導(dǎo)，方便確定測(cè)量位置．在彈性成像模式下，系統(tǒng)控制多個(gè)陣元的超聲波發(fā)射產(chǎn)生聚焦的聲輻射力，打擊測(cè)量點(diǎn)處的角膜，然后使用快速平面波成像方式檢測(cè)組織的振動(dòng)信息和波的傳播．

圖1 超聲實(shí)驗(yàn)過程中眼球與探頭的關(guān)系及角膜組織振動(dòng)位移信息Fig.1 Relationship between eyeball and probe, and vibration displacement information of corneal tissue during ultrasonic experiment

1.3 角膜樣本處理

實(shí)驗(yàn)所用新鮮完整的牛眼球取自深圳市南山肉聯(lián)廠．CXL方法可分為物理類交聯(lián)和化學(xué)類交聯(lián)[13]．物理類交聯(lián)方法主要有核黃素加紫外光誘導(dǎo)和玫瑰紅加綠光誘導(dǎo)交聯(lián)反應(yīng)兩種；化學(xué)交聯(lián)方法主要有利用甘油醛產(chǎn)生非酶糖基化作用硬化角膜和化學(xué)試劑戊二醛硬化作為交聯(lián)劑兩種．本研究使用兩種CXL方案對(duì)離體角膜樣本進(jìn)行硬化，比較處理前后的黏彈性測(cè)量結(jié)果：方案1是紫外光核黃素(UVA-riboflavin)法，角膜表面涂覆質(zhì)量分?jǐn)?shù)為0.1%的核黃素，并給予一定平常的紫外光照射；方案2是在方案1的配方中額外加入質(zhì)量分?jǐn)?shù)為0.1%的戊二醛(glutaraldehyde)增強(qiáng)硬化效果．兩方案中未做CXL處理的角膜均被劃分為空白對(duì)照組(未作任何處理)，去上皮對(duì)照組(僅刮除上皮層)，CXL處理的角膜按照紫外光照射的時(shí)長(zhǎng)分為5、10、30和60min組，共6個(gè)分組．兩個(gè)方案中使用眼球樣本數(shù)(n)如表1，每個(gè)樣本都使用角膜測(cè)厚儀(SP-100，TOMEY，Japan)測(cè)量記錄處理后的角膜厚度來進(jìn)行角膜黏彈性擬合．

紫外光核黃素CXL膠原交聯(lián)溶液配制方法為：

1)備好濃度為0.1mol/L，pH=7.4的磷酸緩沖鹽溶液(phosphatebuffersaline，PBS)；

2)把粉末狀葡聚糖溶入所配置的PBS緩沖液中，得到質(zhì)量分?jǐn)?shù)為20.0%的葡聚糖溶液．(方案2中緩沖液里額外加入質(zhì)量分?jǐn)?shù)0.1%的的戊二醛)；

3)葡聚糖溶液中加入質(zhì)量分?jǐn)?shù)為0.1%的核黃素，配成膠原交聯(lián)溶液．

兩種方案的處理步驟為：

1)刮除角膜上皮，點(diǎn)涂葡聚糖溶液后靜置1h；

2)角膜在膠原交聯(lián)溶液中浸泡5min，紫外光下根據(jù)各小分組的要求定時(shí)照射，間隔5min點(diǎn)涂一次交聯(lián)溶液濕潤(rùn)角膜；

表1 超聲檢測(cè)角膜分組標(biāo)號(hào)及樣本個(gè)數(shù)

3)暗箱內(nèi)靜置45 min，然后用清水清洗角膜表面殘留溶液，用于測(cè)量．

1.4 離體超聲測(cè)量過程及數(shù)據(jù)處理

實(shí)驗(yàn)中將完整眼球固定并浸于純水中，固定探頭，選取角膜頂點(diǎn)為測(cè)量點(diǎn)，一次完整的角膜生物力學(xué)測(cè)量過程為：快速平面波采集10幀靜態(tài)角膜數(shù)據(jù)，聲輻射力激勵(lì)角膜1次，延時(shí)240 μs，快速平面波采集500幀數(shù)據(jù)，處理顯示角膜組織振動(dòng)位移曲線、蘭姆波群速度與蘭姆波頻散分布，保存數(shù)據(jù)．其中，聲輻射力使用了35個(gè)陣元聚焦，發(fā)射頻率為4 MHz的正弦波500個(gè)周期，激勵(lì)電壓最大為58 V；快速平面波采集模式中128個(gè)通道同時(shí)發(fā)射中心頻率為6.25 MHz的短脈沖并接收回波信號(hào)，最大電壓為70 V，脈沖重復(fù)頻率為10 kHz，采集500幀的角膜振動(dòng)數(shù)據(jù)．

對(duì)RF數(shù)據(jù)進(jìn)行正交解調(diào)，參照脈沖多普勒的方法從兩正交信號(hào)中提取角膜振動(dòng)的信息．根據(jù)振動(dòng)傳播路徑上各點(diǎn)的位移，計(jì)算蘭姆波的群速度和相速度．

每個(gè)樣本采集6組數(shù)據(jù)，計(jì)算對(duì)應(yīng)的相速度均值，得到蘭姆波在角膜組織中的傳播頻散特性，結(jié)合式(5)估計(jì)對(duì)應(yīng)的角膜的剪切彈性和剪切黏性值．本研究主要比較各組角膜的蘭姆波頻散曲線及黏彈性的估計(jì)值．用方差分析法(analysis of variance, ANOVA)分析組間均值差異的顯著性．

2 結(jié)果與討論

2.1 兩種方案對(duì)應(yīng)各分組的頻散曲線

方案1各小組的蘭姆波相速度的均值與方差曲線如圖2．圖2顯示R5、R10和R30組的頻散曲線幾乎重合；RC組各頻率相速度明顯低于其他各組，RCE組與其他各組的差異也非常明顯，R60組各頻率相速度均明顯高于其他組．

方案2各組角膜樣本的蘭姆波在200～1 050 Hz上相速度頻散曲線如圖3．由圖3可見，頻散曲線自GC、GCE、G5、G10和G30各組逐漸升高，各組之間差距明顯．

圖2 方案1各分組對(duì)應(yīng)的蘭姆波速度頻散曲線Fig.2 Lamb wave velocity dispersion curves of the groups in CXL scheme 1

圖3 方案2各分組蘭姆波速度頻散曲線Fig.3 Lamb wave velocity dispersion of the groups in CXL scheme 2

2.2 兩種方案的剪切彈性估計(jì)對(duì)比

由每個(gè)樣本6次測(cè)量得到的各頻率下相速度平均值和角膜厚度測(cè)量值，估計(jì)出基于Voigt模型的彈性模量μ1和黏性系數(shù)μ2． 將μ1和μ2代入式(5)，得到蘭姆波相速度的理論值．相對(duì)誤差=(理論值-實(shí)測(cè)值)/實(shí)測(cè)值×100%，對(duì)于所有樣本的相速度，相對(duì)誤差在(7.66±9.09)%，說明理論蘭姆波波速與實(shí)測(cè)結(jié)果吻合較好．

兩種方案各6個(gè)分組的剪切彈性箱線圖如圖4．對(duì)比兩方案分別對(duì)應(yīng)的組別結(jié)果發(fā)現(xiàn)，兩個(gè)方案中正常對(duì)照組與去上表皮組對(duì)照組的彈性結(jié)果非常吻合；方案1中R5、R10、R30和R60組別的彈性值比方案2中對(duì)應(yīng)的G5、G10、G30和G60組別的彈性值分別增了148.78%、335.68%、392.59%和250.47%．方差分析顯示，10、30和60min組兩種方案的彈性具有顯著性差異(P<0.05)， 表明方案2對(duì)角膜硬化程度明顯強(qiáng)于方案1．

圖4 方案1與方案2對(duì)應(yīng)的各分組彈性值箱線圖分布圖Fig.4 Boxplots of estimated shear elasticity for all groups

方案1中各組方差分析結(jié)果顯示，R5與R10，R10與R30組無顯著差異，其余相鄰各組均有顯著差異，說明這6組彈性值呈遞增趨勢(shì)．方案2的彈性方差分析結(jié)果顯示，僅G5組與G10組間無顯著差異，其余各組均有顯著差異．

2.3 兩種方案的剪切黏性估計(jì)對(duì)比

兩種方案對(duì)應(yīng)的黏性結(jié)果顯示，兩方案空白對(duì)照組中的黏性值均明顯大于各處理組中的黏性值，處理后黏性值數(shù)量級(jí)非常小且變化不大．處理后黏性均變小，但各分組間無顯著差異．

由表2兩方案6組的黏性結(jié)果可見，空白對(duì)照組黏性最大，去上皮對(duì)照組黏性降低1～2個(gè)數(shù)量級(jí)，CXL處理后黏性又下降1～2個(gè)數(shù)量級(jí)，但與紫外光照射時(shí)間無明顯相關(guān)．

2.4 結(jié)果的討論

CXL處理中致密的上皮層會(huì)阻礙藥物的滲透，影響交聯(lián)效果，一般都會(huì)刮除去上皮后再進(jìn)行CXL處理，因此本研究在對(duì)照組中設(shè)置了去上皮對(duì)照組．由于角膜是具有多層結(jié)構(gòu)，且每層的力學(xué)特性各不相同．因此，空白對(duì)照組和去上皮對(duì)照組的黏彈性差異主要反映了上皮層對(duì)測(cè)量結(jié)果的影響．而去上皮對(duì)照組與CXL處理后各組之間的差異體現(xiàn)了CXL處理對(duì)角膜力學(xué)特性的改變．

表2 兩方案角膜處理方式中的各分組黏性均值±s)

Tanter等[14]用SSI法測(cè)量到正常豬眼角膜的楊氏模量為(190±32)kPa，經(jīng)UVA-riboflavinCXL處理的角膜楊氏模量為(890±250)kPa，增大了460%．Urs等[15]研究結(jié)果表明，經(jīng)過CXL處理后兔眼角膜硬度明顯增大．Scarcelli等[16]用組織的Brillouin模量來評(píng)估角膜經(jīng)過核黃素-紫外光誘導(dǎo)膠原交聯(lián)后角膜彈性變化，發(fā)現(xiàn)在CXL處理角膜時(shí)，若浸泡時(shí)間不變，紫外光照射時(shí)間越長(zhǎng)，得到處理后角膜Brillouin模量越大．本研究對(duì)牛眼球施以兩種膠原處理，結(jié)果表明，處理后角膜彈性比未處理的角膜彈性顯著增大，并隨紫外光照射時(shí)間的延長(zhǎng)彈性逐漸增大，與文獻(xiàn)[14-16]報(bào)道吻合．

在核黃素溶液中加入戊二醛，相當(dāng)于同時(shí)施加了化學(xué)交聯(lián)和物理交聯(lián)，目的是獲得更強(qiáng)的硬化效果，以評(píng)價(jià)本研究對(duì)角膜硬度的區(qū)分能力．Avila等[17]通過對(duì)離體豬眼球的眼壓變化與體積變化獲取到對(duì)應(yīng)的眼壁硬度系數(shù)，對(duì)應(yīng)的質(zhì)量分?jǐn)?shù)為0.1%的戊二醛法與紫外光核黃素法處理后離體豬眼球眼壁硬度結(jié)果為1 227kPa/L與867kPa/L，即處理相同時(shí)間，戊二醛處理角膜硬化效果比紫外光核黃素法效果更好．本研究結(jié)果顯示，兩種CXL處理方案中，加入戊二醛對(duì)角膜的硬化效果更顯著，CXL處理后的黏性顯著降低．目前尚未見用相似原理進(jìn)行實(shí)驗(yàn)定量估計(jì)角膜黏性系數(shù)的報(bào)道，但有研究發(fā)現(xiàn)膠原交聯(lián)導(dǎo)致硬度增加而黏性降低[18-19]．

上述離體實(shí)驗(yàn)的結(jié)果表明，該超聲測(cè)量方法均可定量區(qū)分出對(duì)照組與處理組間的角膜生物力學(xué)特性變化，不僅能區(qū)分出同一配方不同照射時(shí)間引起彈性的改變，還可區(qū)分出兩方案對(duì)彈性改變的差異．該方法使用同一線陣探頭完成聲輻射力激勵(lì)和檢測(cè)掃描，簡(jiǎn)單便捷，有望成為商用超聲系統(tǒng)的附加功能．

結(jié) 語(yǔ)

本研究基于超聲彈性成像原理，提出一套基于蘭姆波頻散特性估計(jì)角膜黏彈性的方法和系統(tǒng)，評(píng)價(jià)了兩種CXL方案對(duì)角膜處理的交聯(lián)效果，結(jié)果顯示，經(jīng)過CXL處理后彈性增大，黏性降低．方案2處理組紫外光照射5、10、30和60min組的彈性值比方案1中處理組彈性值分別增大了148.78%、335.68%、392.59%和250.47%，表明該超聲方法可很好地定量區(qū)分不同CXL處理方法改變角膜硬度的效果，未來有望用于對(duì)膠原術(shù)治療圓錐角膜的評(píng)估．下一步將著眼于測(cè)量準(zhǔn)確性的驗(yàn)證、黏性估計(jì)的價(jià)值、系統(tǒng)的優(yōu)化、能量安全性評(píng)估和在體動(dòng)物實(shí)驗(yàn)研究．

/ References：

[1] Sykakis E, Karim R, Evans J R, et al. Corneal collagen cross-linking for treating keratoconus [J]. The Cochrane Database of Systematic Reviews, 2015, 3(3): CD10621.

[2] Spoerl E, Huhle M, Seiler T. Induction of cross-links in corneal tissue [J]. Experimental Eye Research, 1998, 66(1): 97-103.

[3] 陳世豪，王勤美．重視角膜膠原交聯(lián)術(shù)再眼科的臨床應(yīng)用和研究[J]．中華眼視光與視覺科學(xué)雜志，2015，17(3)：136-139． Chen Shihao, Wang Qinmei. Attaching importance to the clinical application and research of corneal collagen cross-linking in ophthalmology [J]. Chinese Journal of Optometry Ophthalmology & Visual Science, 2015, 17(3): 136-139.(in Chinese)

[4] Boyce B L, Jones R E, Nguyen T D, et al. Stress-controlled viscoelastic tensile response of bovine cornea.[J]. Journal of Biomechanics, 2007,40(11):2367-2376.

[5] Elsheikh A, Alhasso D, Rama P. Biomechanical properties of human and porcine corneas[J]. Experimental Eye Research, 2008,86(5):783-790.

[6] Gualini M M, Iqbal S, Sixt W, et al. Broader corneal characterization with PulseESPI applied to elasticity measurements[J]. European Journal of Ophthalmology, 2010,20(2):306-309.

[7] Scarcelli G, Pineda R, Yun S H. Brillouin optical microscopy for corneal biomechanics[J]. Invest Ophthalmol & Visual Science, 2012,53(1):185-190.

[8] Qi Wenjuan, Li Rui, Ma Teng, et al. Confocal acoustic radiation force optical coherence Elastography elastography using a ring ultrasonic transducer[J]. Applied Physics Letters, 2014,104(12):123702.

[9] Tanter M, Touboul D, Gennisson J L, et al. High-resolution quantitative imaging of cornea elasticity using supersonic shear imaging[J]. IEEE Transactions on Medical Imaging, 2009,28(12):1881-1893.

[10] Zhang Xinyu, Yin Yin, Guo Yanrong, et al. Measurement of quantitative viscoelasticity of bovine corneas based on Lamb wave dispersion properties[J]. Ultrasound Medicine & Biology, 2015,41(5):1461-1472.

[11] Krauklis P V, Molotkov L A. Low-frequency Lamb waves in cylindrical and spherical layers in an elastic medium[J]. Journal of Mathematical Sciences, 1975,3(1):82-90.

[12] Kanai H. Propagation of spontaneously actuated pulsive vibration in human heart wall and in vivo viscoelasticity estimation propagation of spontaneously actuated pulsive vibration in human heart wall and in vivo viscoelasticity estimation[J]. IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, 2005,52(11):1931-1942.

[13] Greenstein S A, Hersh P S. Corneal collagen cross-linking for keratoconus and corneal ectasia[M]. Berlin: Springer Berlin Heidelberg, 2014: 71-87.

[14] Tanter M, Touboul D, Gennisson J L, et al. High-resolution quantitative imaging of cornea elasticity using supersonic shear imaging[J]. IEEE Transactions on Medical Imaging, 2009,28(12):1881-1893.

[15] Urs R, Lloyd H O, Silverman R H. Acoustic radiation force for noninvasive evaluation of corneal biomechanical changes induced by cross-linking therapy[J]. Journal of Ultrasound in Medicine: Official Journal of the American Institute of Ultrasound in Medicine, 2014,33(8):1417-1426.

[16] Scarcelli G, Kling S, Quijano E, et al. Brillouin microscopy of collagen crosslinking: noncontact depth-dependent analysis of corneal elastic modulus[J]. Investigative Ophthalmology & Visual Science, 2013,54(2):1418-1425.

[17] Avila M Y, Gerena V A, Navia J L. Corneal crosslinking with genipin, comparison with UV-riboflavin in ex-vivo model[J]. Molecular Vision, 2012,18:1068-1073.

[18] Kling S, Remon L, Pérez-Escudero A, et al. Corneal biomechanical changes after collagen cross-linking from porcine eye inflation experiments[J]. Investigative Ophthalmology & Visual Science, 2010,51(8):3961-3968.

[19] Elsheikh A, Wang D, Brown M, et al. Assessment of corneal biomechanical properties and their variation with age[J]. Current Eye Research, 2007,32(1):11-19.

【中文責(zé)編：英 子；英文責(zé)編：子 蘭】

2017-04-12；Accepted：2017-04-29

Associate professor Zhang Xinyu. E-mail: xyzhang9@szu.edu.cn

Assessing the effect of collagen cross-linking by quantitative measurement of corneal viscoelasticity

Zhang Pengpeng, Gao Xuehua, Wang Qingmin, Lyu Zhen, and Zhang Xinyu

National-Regional Key Technology Engineering Laboratory for Medical Ultrasound, Guangdong Key Laboratory for Biomedical Measurements and Ultrasound Imaging, Health Science Center, Shenzhen University, Shenzhen 518060, P.R.China

This study establishes an experimental platform based on VerasonicsTMVantage 256 ultrasonic open system, and uses the acoustic radiation force to excite the corneal vibration and uses plane wave imaging to detect A0 mode Lamb wave. The corneal shear elastic modulus and shear viscosity are estimated based on Lamb wave dispersion characteristics. In order to validate the effectiveness of this system, two kinds of CXL treatments are introduced to treat bovine corneas, and the effect of the treatment with different reagent and UV exposure duration are evaluated. The results show that the CXL with the addition of glutaraldehyde (0.1% mass concentration) in riboflavin solution induces significantly higher increase in elasticity in comparison with those without glutaraldehyde. The elasticity of the group with glutaraldehyde after 5, 10, 30 and 60 min UVA exposure increases by 148.78%、 335.68%, 392.59% and 250.47%, respectively, in comparison with the groups without glutaraldehyde.

biomedical measurement; medical ultrasound; ultrasound elastography; acoustic radiation force; cornea; viscoelasticity; collagen cross-linking

：Zhang Pengpeng, Gao Xuehua, Wang Qingmin, et al. Assessing the effect of collagen cross-linking with quantitative measurement of corneal viscoelasticity [J]. Journal of Shenzhen University Science and Engineering, 2017, 34(4): 428-433.(in Chinese)

R 318.01; R 445.1

A

10.3724/SP.J.1249.2017.04428

國(guó)家自然科學(xué)基金資助項(xiàng)目(61101025，81570552)；深圳市基礎(chǔ)研究基金資助項(xiàng)目 (JCYJ20150625102427087)

張鵬鵬(1991—)，男，深圳大學(xué)碩士研究生.研究方向：醫(yī)學(xué)超聲．E-mail：ZhangPengpengSZ@163.com

Foundation：National Natural Science Foundation of China(61101025，81570552)；Shenzhen Basic Research Project(JCYJ20150625102427087)

引 文：張鵬鵬，高雪華，王青敏，等．角膜黏彈性定量測(cè)量評(píng)估膠原交聯(lián)處理效果[J]. 深圳大學(xué)學(xué)報(bào)理工版，2017，34(4)：428-433.