高强度聚焦超声经颅脑肿瘤治疗焦域的仿真研究*

为了研究脑组织和脑肿瘤组织在HIFU治疗时形成焦域的特性及坏死肿瘤组织、不同治疗剂量参数以及多次治疗时焦点间距和时间间隔对HIFU形成焦域的影响,本文以脑胶质瘤患者为例,利用患者头颅CT图像数据建立HIFU经颅治疗的仿真模型,基于Westervelt声波非线性传播方程和Pennes生物热传导方程进行HIFU经颅治疗的仿真研究。结果表明,脑组织和脑肿瘤组织内形成HIFU焦域的差异较小,坏死肿瘤组织对HIFU温度场分布有较小影响;辐照声强越大,焦点温升达到同一温度所用时间越短,焦域长短轴越短,颅骨处温升越低;当焦点间距在一定范围内时,第一次辐照形成的温度场分布对第二个焦点温升达到同一温度所需时间影响较大;两次辐照时间间隔对颅骨处温升和两次聚焦形成60?C以上温度分布影响较小。     

剪切波对HIFU经颅聚焦形成温度场影响的数值仿真研究

在高强度聚焦超声经颅治疗时,既有纵波又有剪切波,为了保障该治疗方法的安全有效性,有必要分析剪切波对HIFU治疗温度场的影响。该文基于人体头颅CT数据和曲率半径为150 mm的256阵元的半球相控换能器建立三维高强度聚焦超声经颅声波传播模型,利用时域有限差分法结合Westervelt声波非线性传播方程、动量方程、质量守恒方程和Pennes生物热传导方程数值仿真其形成温度场,研究在相同输入功率、不同聚焦角度条件下对应阵元数进行激励时,剪切波对换能器形成温度场的影响。结果表明,随换能器聚焦角度减小,在几何焦点处形成的焦域面积逐渐增大,考虑剪切波形成的温度场达到65?C所需时间逐渐延长,焦点前移程度越大;在相同聚焦角度条件下,考虑剪切波的温度场达到65?C所需时间更短,旁瓣更少,在颅骨处的温度更高,对焦点前移几乎没有影响;随换能器聚焦角度减小,考虑剪切波的模型形成的焦域面积变化范围更大;幂指数函数形式对不同聚焦角度下焦域面积大小的拟合优度高,可预测不同聚焦角度换能器形成的焦域面积。     

基于电阻抗层析成像的高强度聚焦超声温度监测技术

作为一种对正常组织无损伤且不易引起癌细胞转移的非入侵肿瘤治疗手段,高强度聚焦超声(HIFU)治疗过程中焦域的温度监测是实现剂量精准控制的关键.本文基于生物组织的温度-电阻抗的关系,将电阻抗层析成像(EIT)和HIFU治疗相结合,提出了一种利用组织焦平面的表面电压实现电阻抗重构的检测技术.建立了HIFU治疗和EIT综合系统模型,在考虑组织的声吸收条件下,对三维Helmholtz方程在柱坐标下的声场计算进行了二维简化,并引入Pennes生物热传导方程来计算HIFU焦域的声压和温升分布特性;引入生物组织的温度-电阻抗关系,基于麦克斯韦电磁场理论,建立了具有温度分布HIFU焦域的电流和电压计算模型,利用恒流注入的边界条件实现电场计算,获得焦平面的表面电压分布.在数值计算中,利用实验聚焦换能器参数,模拟了在固定声功率下组织焦域的声场和温度场分布,以及中心和偏心聚焦条件下不同治疗时刻的电导率分布;然后通过对称电极的循环电流注入,计算了组织模型焦平面内的电流密度和电势分布,获得了焦平面圆周分布的表面电极电压;进一步采用修正的牛顿-拉夫逊算法,利用32×32的表面电极电压实现了焦平面内电导率分布的重建.结果表明,基于温度-电阻抗关系的EIT电导率重建技术不但能准确定位HIFU焦域中心,还能恢复HIFU治疗中焦域的温度分布,证明了EIT用于HIFU治疗中温度监测的可行性,为其疗效评估和剂量控制提供了一种无创电阻抗测量和成像新方法.     

变厚度聚焦换能器对声焦域轴向长度影响的研究*

针对不同厚度的病变组织,改变声焦域轴向长度能提高高强度聚焦超声在临床治疗过程中的安全性和有效性。基于多频超声波叠加原理,该文提出了变厚度(多频)聚焦换能器,并设计了两种类型变厚度聚焦换能器。根据瑞利积分法推导了变厚度聚焦换能器声场,计算和分析了变厚度聚焦换能器的声焦域轴向长度,并与等厚度(单频)聚焦换能器声焦域轴向长度进行对比。结果显示,变厚度聚焦换能器中心到边缘的厚度变化趋势与声焦域轴向长度变化相关,中间薄两边厚换能器声焦域轴向长度缩短,中间厚两边薄换能器声焦域轴向长度变长,且实验验证了理论的正确性。研究结果可为变厚度聚焦换能器声场研究和高强度聚焦超声的临床治疗提供参考。     

基于分数导数研究高强度聚焦超声的非线性声场

在高强度聚焦超声(high intensity focused ultrasound, HIFU) 的研究中, 生物组织的衰减和色散性质会对声能量的空间分布产生影响. 本文提出应用分数导数修正非线性Khokhlov-Zabolotskaya-Kuznetsov (KZK)方程, 研究生物组织中非线性HIFU声场. 对三种生物仿体的衰减和声速色散的理论实验研究表明分数导数应用的可行性, 在此基础上通过数值仿真分析研究了衰减及声速随频率的变化对HIFU焦域分布的影响. 研究结果表明, 在计算强非线性聚焦超声时, 由于高次谐波的强色散作用, 引入分数导数来解决生物组织特殊的衰减以及色散问题可进一步提高HIFU治疗的安全性. 关键词： 分数导数 声衰减 色散 高强度聚焦超声     

微泡对高强度聚焦超声焦域影响的研究

微泡对高强度聚焦超声(HIFU)治疗具有增效作用,而HIFU治疗中不同声学条件下微泡对HIFU治疗焦域的影响尚不清楚。本文基于声传播方程、Yang-Church气泡运动方程、生物热传导方程、时域有限差分法(FDTD)、龙格-库塔(RK)法数值仿真研究输入功率、激励频率和气泡初始半径对HIFU在含气泡体模中形成焦域的影响,并利用含Sono Vue造影剂的仿组织体模研究进行实验验证。结果表明,增大输入功率、气泡初始半径和升高激励频率均可增大焦域,随着输入功率的增大,焦域形状可能发生变化,而随着激励频率升高和气泡初始半径的增大,焦域会向远离换能器的方向移动。     

高强度聚焦超声治疗子宫肌瘤中水囊耦合的影响评估*

水囊作为一种特殊的辅助器件在高强度聚焦超声治疗子宫肌瘤的过程中发挥着重要作用,但它同时也会让声通道变得更加复杂,对超声治疗效率和超声引导影像均造成负面影响。基于实际临床场景,利用k-Wave声学仿真软件建立高强度聚焦超声辐照靶组织及其超声引导过程的仿真模型,在有无水囊两个场景下,对治疗效率和监控影像质量进行定量评估。结果表明,加入水囊后,声波的聚焦性变差,焦域处的最高温度降低;通过分析超声图像的分辨率、对比度、对比噪声以及背景信噪比,发现加入水囊后超声图像的质量降低。建立的数值仿真模型能够初步评估水囊对治疗效率和监控影像质量的影响,可作为一种评估体系用于优化水囊参数,如水囊材料、厚度、形状、内部溶液成分等,为下一步实验探索对治疗效率和影像质量影响最小的水囊耦合方式提供评估手段。     

基于换能器驱动信号特征的高强度聚焦超声焦域损伤实时监测

高强度聚焦超声(high intensity focused ultrasound,HIFU)焦域的实时监测是聚焦超声临床治疗面临的关键问题,目前临床常采用B超图像强回声的变化实现焦域组织损伤的监测,而B超图像出现的强回声大多与焦域处的空化及沸腾气泡相关,无法准确、实时地监测治疗状态.HIFU治疗中焦域组织会伴随温度升高、空化、沸腾和组织特性等变化,换能器表面的声学负载也在不断变化,针对该问题,本文构建了换能器电压电流实时检测平台,通过测量换能器电学参数来感知焦域组织的状态变化.以离体牛肝组织作为HIFU辐照对象,并将相位差变化的结果与离体牛肝组织损伤的结果进行了对照,实验结果表明,在HIFU辐照过程中,换能器电压与电流的相位会出现由相对平稳到大幅波动的过程,此时停止辐照可见焦域出现明显损伤,而此时B超图像灰度无明显变化;此外,当焦域出现空化时,其波动幅度与范围将较之更大.此方法可为HIFU焦域组织损伤监测提供一种新的研究方案和手段.     

高频聚焦超声声场和温度场的仿真研究

为探究临床常用的7 MHz高频聚焦超声在多层生物组织中的声传播以及毫秒级时间内的生物传热规律问题,基于Westervelt方程和Pennes传热方程,使用有限元方法建立高频聚焦超声辐照多层组织的非线性热黏性声传播及传热模型。首先分析了线性模型和非线性模型之间的差异,然后在非线性模型下探究换能器的参数对声场和温度场的影响。仿真结果显示：在7 MHz频率下,当换能器输出声功率超过5 W时,声波传播的非线性效应不可忽视(p <0.05);当声功率从5 W增大到15 W时,非线性模型与线性模型预测的温度偏差从20%增加到34.703%;高频聚焦超声波的非线性行为比低频更加显著,基频能量向高次谐波转移的程度增大,声功率为10 W和15 W时4次谐波与基波之比分别达到7.33%和12.12%;高频换能器参数的改变对组织中声场和温度场分布的影响较大,换能器焦距从12 mm减小到11.2 mm,焦点处最高温度增加了77%。结果表明,7 MHz聚焦超声的非线性声传播需要考虑到4次谐波的影响。该文提出的多层组织非线性仿真模型可为高频聚焦超声换能器参数优化及制定安全、有效的术前治疗方案提供理论参考。     

时间反转聚焦在组织中形成损伤的数值建模

为考察基于时间反转方法的高强度聚焦超声治疗在预设目标点处的组织损伤情况,使用三维有限差分算法求解Westervelt方程,建立非线性声波传播数值模型,采用97阵元相控阵结合虚拟源的时间反转方法进行超声聚焦,分析其形成的声场和热场,并考察目标点偏离轴线时的组织损伤形成规律。结果表明随着目标点偏轴距离的增大,声压旁瓣开始增多。旁瓣的温升较低,不足以形成组织损伤。时间反转方法可用于多点聚焦,在一定的范围内,形成多点目标损伤而不产生额外的周围组织损伤。同时多点聚焦可以形成一个较大的损伤区域,减少超声治疗时间。     

Self-organization of five species in a cyclic competition game

Cyclic competition game models, particularly the “rock–paper–scissors” model, play important roles in exploring the problem of multi-species coexistence in spatially ecological systems. We propose an extended “rock–paper–scissors” game to model cyclic interactions among five species, and find that two of the five can coexistent when biodiversity disappears, which is different from the “rock–paper–scissors” game. As the number of fingers is five, we named the new model the “fingers” game, where the thumb, forefinger, middle finger, ring finger, and little finger cyclically dominate their subsequent species and are dominated by their former species. We investigate the “fingers” model in two ways: direct simulations and nonlinear partial differential equations. An important finding is that the number of species in a cyclic competition game has an influence on the emergence of biodiversity. To be specific, the “rock–paper–scissors” model is in favor of maintaining biodiversity in comparison with the “fingers” model when the variables (population size, reproduction rate, selection rate, and migration rate) are the same. It is also shown that the mobility and reproduction rate can promote or jeopardize biodiversity.     

Contrastive analysis on tone between Putonghua and Taiwan Mandarin

The differences of the pitch and duration of Chinese syllables between Putonghua （PTH） and Taiwan Mandarin （TM） were studied. The speech materials to be used are not only isolated syllables, but also sentences. The results reveal that： For the isolated syllables, T1 and T2 in TM are influenced by Minnan dialect, therefore their pitch are lower than those in PTH. T3 is fall-rise in PTH, while it is fall in TM. Moreover, the syllable duration sequence for different tone is T3〉T2〉T1〉T4 in PTH, while it is T1〉T2〉T3〉T4 in TM. For the syllables in sentences, T2 is mid-rise in PTH, while it is mid-level in TM. And the T3 is longer than T4 but shorter than T1 or T2 in PTH, while it is the shortest in TM. Furthermore the effects of prosodic phrase boundary on duration for different tones are almost the same in PTH, but the lengthening part of T1 or T2 is longer than that of T3 or T4 in TM.     

Experimental study on downstep in Chinese intonation

Downstep in pitch contour of Chinese Putonghua is examined using subtly designed sentences by controlling tone combinations. The results show both automatic and nonautomatic downstep phenomena exist in Chinese. In non-automatic downstep, low tones compress downwards the pitch range of the following syllables. and the main influence of downstep is on topline. Low tone not only lower the topline behind it, but also raise the high tones before it, the effects are compatible with each other. In automatic downstep, the topline of pitch contour in intonational phrase is presented as a linear downtrend, but it differs among speakers due to the effect of personal stress practice. In comparison with downstep phenomenon in other tone or non-tone languages, the downstep ratio in Chinese is not constant, and the domain of downstep is not limited within the adjacent tones.     

Abelian-Higgs phase of SU(2) QCD and glueball energy

It is shown that SU(2) QCD admits an dual Abelian-Higgs phase, with a Higgs vacuum of a type-Ⅱ superconductor. This is done by using a connection decomposition for the gluon field and the random-direction approximation. Using a bag picture with soft wall, we presented a calculational procedure for the glueball energy based on the recent proof for wall-vortices [Nucl. Phys. B 741(2006)1].     

Theoretical analysis and simulation of conjugate heights for dual-conjugate AO system in lidar

A multi-conjugate adaptive optics （MCAO） can offer a possibility of widening field of view （FOV） characterized by the isoplanatic angle, and the choose of conjugate height becomes a basic problem for MCAO, which influences the size of iosplanatic angle. Considering the application of lidar, the isoplanatic angle＇s expressions of two deformable mirrors （DMs） MCAO for uplink and downlink are deduced. The effects of conjugate heights for dual-conjugate AO are thoughtfully discussed, and the isoplanatic angles are further analyzed. The results show that the isopanatic angle varies with the conjugate height and reaches the maximum as the conjugate height is at the optimal altitude. Moreover, the optimal conjugate height changes with the propagation distance.     

Quick single-photon detector with many avalanche photo diodes working on the time division

Due to the limit of response speed of the present single-photon detector, the code rate is still too low to come into practical use for the present quantum key distribution (QKD) system.A new idea is put up to design a quick single-photon detector.This quick single-photon detector is composed of a multi-port optic-fiber splitter and many avalanche photo diodes (APDs).Au of the ports with APDs work on the time division and cooperate with a logic discriminating and deciding unit driven by the clock signal.The operation frequency lies on the number N of ports, and can reach N times of the conventional single-photon detector.The single-photon prompt detection can come true for high repetition-rate pulses.The applying of this detector will largely raise the code rate of the QKD, and boost the commercial use.     

Signal generation in an isotropic medium in scanning electron acoustic microscope

Based on the research in Ref. [5][Materials Science and Engineering, 1989; A122： 57 63], an improved model of heat source is set up, the different modes of Lamb wave in an isotropic sample generated by a chopped electron beam at frequency f are obtained with integral transform and normal function expansion method, and the output signal of PZT coupled at the back surface of the sample is found out. The generation mechanism of SEAM （Scanning Electron Acoustic Microscopy） signal is discussed. It shows that the SEAM is a near field imaging technique with high spatial resolution and its best lateral spatial resolution is about 2√2α （α is the radius of the focused electron beam）. Some of experimental results of SEAM images are presented in the paper and it shows that the spatial resolution of SEAM is better than 0.5 μm and smaller than the thermal diffusion length of the sample. Therefore the character of near field imaging in SEAM is also proved experimentally.     

High efficiency filtering using two-dimensional photonic crystal coupled cavity structures

Transmission spectra of coupled cavity structures (CCSs) in two-dimensional (2D) photonic crystals (PCs) are investigated using a coupled mode theory, and an optical filter based on CCS is proposed. The performance of the filter is investigated using finite-difference time-domain (FDTD) method, and the results show that within a very short coupling distance of about 3λ, where λ is the wavelength of signal in vacuum, the incident signals with different frequencies are separated into different channels with a contrast ratio of 20 dB. The advantages of this kind of filter are small size and easily tunable operation frequencies.