我是靠谱客的博主 谦让小丸子,最近开发中收集的这篇文章主要介绍转-超声系统的信号链设计注意事项,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

原文链接:https://kuaibao.qq.com/s/20190802A0K75K00?refer=spider    

       高性能超声成像系统广泛应用于各种医学场景。在过去十年中,超声系统中的分立电路已经被高度集成的芯片(IC)所取代。先进的半导体技术不断推动系统性能优化及尺寸小型化。这些变革都得益于各类芯片技术,如专用低噪声放大器、多通道低功耗 ADC、集成高压发射、优化的硅工艺和多芯片模块封装。随着芯片功耗和尺寸减小至原来的 20%。
         
       此外,得益于低功耗、高性能硅工艺的发展,部分波束合成预处理模块已经集成于通用的模拟或混合信号芯片而非专用的数字处理器。同时,先进的高速串行或是无线接口大大降低了系统布局复杂度,并且能够将尽可能多的 RF 数据转移到系统集成芯片(SOC)、CPU 或 GPU。当前超声技术的应用也从特定的放射学诊断扩展到各类便携式应用,床旁实时监测以及医疗现场就地检查等各个领域。
         
        本应用指南综述了超声系统的架构和原理,分析了系统设计的注意事项,综述了应用于超声芯片的先进技术,最后讲解了医学超声芯片的模拟参数。  
         
1. 医学超声成像
        超声波是一种频率高于 20KHz 的声波。医学超声成像系统常采用 1 MHz 至 20 MHz 的频率,可达到亚毫米级分辨率。第一台商用超声成像系统诞生于 20 世纪 70 年代,可提供实时的 2D 亮度或灰度图像。如今,超声成像凭借安全性、成本效益和实时方面的优势,已经成为重要的医学成像技术。医学超声系统能够有效地监测婴儿发育,也可用于诊断心脏、肝脏、胆囊、脾脏、胰腺、肾脏、膀胱等内脏器官的疾病。
         
        典型的超声系统包括压电换能器、电子电路电路、图像显示单元和 DICOM(医学数字成像和通信)兼容软件。典型超声系统的简化框图如下所示。
         
                                                 
        图 1. 典型超声系统的简化框图
         
2. 声波产生和传播的原理
        超声换能器是超声系统的关键组成部分,由压电元件、连接器和支撑结构组成。压电效应是指某种材料的物理尺寸随施加的电场而变化的现象,反之亦然。如下所示,超声应用中的大多数换能器是双共模式。换能器在发射相(模式)期间将电能转换成机械能。产生的机械波向介质传播,若介质不均匀则会反射。在接收模式中,接收反射的机械波形并由换能器转换成电信号。
         
                                                
        图 2. 换能器振动、声波传播和反射
         
        在换能器被电子激励之后,会产生声波并在介质中传播。在医学超声中,FDA(食品药品管理局)要求所有成像系统满足瞬时、峰值和平均强度的限制。
         
        我们通常将换能器灵敏度或换能器插入损耗(IL)定义为接收(Rx)和发射(Tx)信号幅度之间的比率,如下所示:
         
                                 
                     
        换能器频率由压电材料 L 的厚度和材料中的声速 cm 决定:
         
                                  
         
        如前所述,常用的频率范围为 1MHz 至 20MHz。基于上述方程式,较高频率的换能器需要较薄的材料。因此,构建极高频的换能器具有一定的挑战性。
         
        换能器频率响应或带宽是另一个关键参数。作为一般规则,若换能器被脉冲信号(即短尖峰)激励,则接收回波的持续时间决定了换能器的带宽。具有极快响应(即短回波)的换能器是宽带换能器,反之亦然。在大多数应用中通常优选更宽的带宽。在相同的换能器频率下,宽带换能器可实现更好的轴向分辨率,因为回波长度决定了超声系统的轴向分辨率。与此同时,宽带换能器适用于谐波成像,在该成像模式下超声能量以基频发射,而图像由接收到回波的二次谐波来重建。如没有宽带宽换能器的情况下,换能器灵敏度在其谐波频率点 2f0 处显著下降。因此许多换能器研究人员不断探索新材料、新架构和新制造工艺以进一步改善换能器性能。
         
        在超声成像的早期阶段,用于超声系统的多通道电子电路既昂贵又不成熟。由电机驱动通过机械扫描方式成像单阵元换能器被广泛用于获得二维(2D)图像。由于机械结构的速度和精度限制,早期系统无法实现高帧率或高精度成像。如今,成熟的阵列换能器和多通道电子技术可支持 64 到 512 个阵元的换能器。以电子扫描为基础可获得高达> 100 帧 / 秒的图像。为实现电子扫描,波束合成应用技术于聚焦换能器的声束。波束合成的细节将在下一节中讨论。与光学成像系统类似,超声系统可在聚焦焦点处实现最佳空间分辨率。根据应用,一维(1D)阵列换能器包括线性阵列、弯曲线性阵列和相位阵列。这些换能器之间的主要区别在于光束成形结构、成像范围和图像分辨率。此外,由超过 2000 个元件组成的最新 2D 阵列换能器可支持实时三维(3D)成像。下图所示为单阵元换能器、1D 阵列换能器和 2D 阵列换能器。
         
                                            
        图 3. 典型的换能器。(A)单元件换能器 ;b)1D 阵列换能器;(c)2D 阵列换能器(由 USC、Vermon 和 Philips 提供)。
         
3. 换能器指标与图像质量
        和任何成像系统类似,图像质量是医学超声成像中的重要标准。诸如空间分辨率和成像穿透等共同参数主要通过换能器指标和声波传播理论来决定。超声图像的纵向和横向分辨率与介质中的声波波长成线性关系:
         
                                            
         
        方程式中,c 是介质中的声速,Zf 是焦距,2r 是换能器孔径或直径。当换能器被脉冲信号激励时,τ-6dB 为接收回波的 -6dB 脉冲宽度的持续时间。 τ-6dB 也与波长λ成线性关系。对于宽带阵列换能器,我们可分别比较 5MHz 和 12MHz 的横向分辨率,其工作频率为 5MHz 至 14MHz。成像深度为 5 厘米。在两种情况下,64 个换能器阵元形成有效孔径。元件之间的间距为 0.3mm。介质中的声速为 1540m/s。有效孔径尺寸为 19.2mm。根据超声系统的信号链设计注意事项的公式,对于 5MHz 和 12MHz 的声波,λ分别为 0.31mm 和 0.13mm。根据上述方程式,横向分辨率分别在 5MHz 时为 0.8mm,在 12MHz 时为 0.33mm。因此,更高频率的应用实现更佳的分辨率。
         
        实际上,仅通过增加换能器频率来改善图像质量并非完全可行。一方面,更高频率的换能器需要更薄的压电材料,这需要更精密的制造技术,且成本更高。另一方面,如后面章节所示,较高频率的声波在生物组织中容易衰减。
         
        当介质不均匀时,声波的部分能量可在两个介质的边界处反射。未反射的声波继续传播,直到它在下一个边界被反射,或完全衰减。反射和透射系数由这两种介质的声阻抗(Z=ρc)的差异决定。方程式中,ρ和 c 分别是介质的密度和声速,假设波传播方向垂直于边界。
         
                                          
         
        表 1 所示为所选生物组织、水和空气的特性。在两个声阻抗极其不同的情况下会出现强反射信号。骨骼密度高,声速快;因此它总是超声图像中的强反射组织器。另一方面,血液和肝脏的声阻抗相似,因此这两种组织之间的反射很弱。只有高灵敏度的换能器才能拾取微弱的信号。如表 1 所示,信号在传播过程中会衰减。累积衰减随着传播距离的增加而增加。以方程式 7 计算衰减,其中系数 2 体现了声波双向传播。
         
                                           
         
        在超声波探查体内组织的典型应用中,来自人体表面的回波与来自内部器官的回波之间的动态范围很容易超过 100dB。我们可假设平均衰减系数为 0.7dB/MHz×cm 和 7.5MHz 换能器。在 10cm 的深度处,基于方程式 7,即 7.5×0.7×10×2dB,计算所得 105dB 的衰减。假设表面回波为 1Vpp,体内器官回波的幅度为<10UVpp,非常弱。该示例表明超声信号具有极宽动态范围才能表征皮肤表面至内脏器官的生理结构的差异。因此需要复杂的电子电路以提供足够的动态范围,而这在有限功耗的预算下是不容易实现的。
         
                                                
        表 1:典型组织和介质的声学特性
         
4. 超声成像模式
        当换能器接收到回波时,适当的处理单元需要将这些信号转换成超声波检验师或其他最终用户的可理解的图像信息。超声成像使用几种成像模式来研究组织特征、体液分布及流动、器官功能等。
         
A 模式和 B 模式
        在最早的超声系统中,通过显示回波的幅度及其时域信息来指导临床诊断。即 A 模式(振幅模式)超声成像系统,如下图所示,它以一维逐行扫描为基础,。由于人类视觉对图像更敏感 ,因此开发亮度或灰度成像模式有更积极的临床意义。为构建 2D 图像,需要在特定区域上扫描换能器的声束,且在扫描期间获得多个 A 模式扫描行。这些扫描行构成一帧图像,沿着扫描线的回波幅度以线性或非线性方式映射到像素值。当换能器的声束足够快地进行扫描时,可实现实时图像。这些图像被称为 B 模式(亮度模式)图像,其创建了与扫描方向平行的一个横截面图像。
         
                                                
        图 4. 扫描模式。(a)A 模式扫描行,(b)B 模式图像,(c)3D 声束扫描,以及(d)B 模式(子图 1、2、3)和 3D(子图 4)临床图像(由 Philips 提供)

]         

         
波束合成器
        波束合成器包括发射和接收波束合成器来实现电子聚焦和控制多阵元换能器的声束。,如下图所示从一个换能器阵元到目标的距离与从另一个阵元到目标的距离不同;因此在发送相,针对每个阵元适当地延迟发送的信号,以使得发射器信号同时到达目标并在目标处产生最高声强,也就是获取最强回波。在接收阶段,通过对接收到的回波应用适当的延迟,以实现线性叠加来自多个换能器阵元的回波,以实现最高灵敏度。
         
                                                  
        图 10. 用于在(a)发射相和(b)接收相中聚焦声束的换能器波束合成器
         
        由于发射电路主要是数字型,因此发射延迟实现是通过现场可编程门阵列(FPGA)或数字信号处理器(DSP)等高速计数器完成的。由于接收信号的复杂性,接收波束合成器显然需要更多的算法优化得以实现。早期基于分立晶体管电子电路信号处理能力有限。因此,接收器波束合成波束合成器以基于电感电容组合的模拟延迟线来实现。在 20 世纪 80 年代,接收器波束合成器开始使用多通道模数转化芯片和数字波束合成技术。
         
                                           
        图 11. 数字波束合成器图解
         
        在目前的主流超声系统中,接收波束合成器一般都是数字型的。数字波束合成器通常在具有极高的计算能力的 FPGA、DSP、PC 或 GPU(图形处理单元)中实现。如前所述,较大的换能器孔径可实现更佳的分辨率。因此,在高端超声系统中,256 个换能器元件形成一个聚焦光束,以获得精细分辨率图像。因此高端波束合成器所需的计算能力相当复杂。
         
        生物组织在形状、密度、声速等方面存在异质性。实时延迟计算和校准会基于所涉组织的声学特性和形状。由于波束合成器设计的重要性和复杂性,大多数超声波公司都有自己的 IP。在不影响波束形成性能的情况下,简化波束形成器设计仍然是一个热门话题。相信正在研发新的波束合成器架构将广泛用于未来的超声系统。
         
数字信号处理
        超声信号需要大量的信号处理,以便从原始超声数据中提取各成像模式所需的信息。主要处理模块包括 B 模式图像重建、基于快速傅里叶变换的多普勒频谱信息提取、基于自相关和互相关的彩色多普勒计算、超声图像扫描坐标转换(2D 超声坐标到笛卡尔坐标)、图像增强等。目前,商用处理器,如现场可编程门阵列(FPGA),数字信号处理器(DSP),被广泛使用,。FPGA 使系统设计人员能够对内部逻辑门进行硬连线,并优化其算法的效率。另一方面,DSP 为系统设计人员提供预定义的标准计算模块,可实时更改和优化他们的算法。换言之,FPGA 以硬件效率取胜,而 DSP 以软件灵活性得宠。新的信号处理器,如 PC 和 GPU;它们的计算能力高于 FPGA 和 DSP,而软件开发成本则大大低于 FPGA 和 DSP 然而,由于 PC 和 GPU 的高功耗,它们并不一定适合低功耗便携式系统。
         
6. 模拟前端芯片设计中的工艺选择
        在任何 AFE 设计之前,半导体工艺选择始终是基于设计目标的首要关键考虑注意事项。CMOS 和 BiCMOS 工艺是超声模拟前端设计中最常用的工艺。它们中的每一个都有其自身优点,且适用于相应的电路块。
         
        BiCMOS(双极 CMOS)工艺目前比纯双极工艺更受欢迎,因为它包含用于模拟设计的高性能双极晶体管和用于数字设计的 CMOS 元件。双极晶体管适用于低噪声放大器设计,具有超低 1/f 噪声、宽带宽和良好的功耗 / 噪声效率。双极性工艺还降低了电路电容,以获得良好的总谐波失真。因此,基于双极或 BiCMOS 工艺的放大器可在比基于 CMOS 工艺的放大器小得多的区域和更低的功耗下实现相同的性能。
         
        德州仪器的 0.35um BiCMOS 工艺用于研究双极和 CMOS 器件之间放大器设计的性能影响。下图(a)表明,基于双极晶体管的放大器在相同的偏置电流下可实现更低的噪声;它还说明双极晶体管具有超低 1/f 噪声特性,这对于具有调制和解调电路的多普勒应用至关重要;(b)与类似的 CMOS 设计相比,双极设计显著减小了面积。当然,由于半导体工艺的特征尺寸减小,在 0.35um BiCMOS 工艺和<0.35um CMOS 工艺之间的面积差异变小。然而,一般来讲,由于上述优点,0.35um BiCMOS 工艺仍然极其适于放大器设计。
         
                                                
        图 12. 基于 CMOS 与 BiCMOS 工艺设计的比较
         
        当电路具有更多数字内容和开关元件(如中速 ADC)时,CMOS 工艺更适合。医学超声信号频率处在 1~20MHz 范围内,其 ADC 采样率通常低于 100MSPS,目前大多数 CMOS 工艺都可轻松处理。采用 0.18um~65nm CMOS 工艺,ADC 设计可实现更佳的集成和功耗降低。此外,与可比较的 BiCMOS 工艺相比,CMOS 工艺通常成本更低,且实现更短的制造周期。所有这些都表明 CMOS 工艺适用于超声 AFE 中的 ADC 设计。
         
        总之,当降低噪声 / 功耗是主要目标时,BiCMOS 工艺适于超声 AFE 中的 TGC 放大器设计,即压控放大器(VCA)设计。另一方面,CMOS 工艺是在 ADC 设计中实现低功耗和高集成度的良好选择。特别是在 0.18um 至 65nm 的节点,与 0.35um BiCMOS 工艺相比,具有完整低压数字库的 CMOS 工艺可以具有竞争力的成本获得更高的集成度。
         
        很明显,BiCMOS VCA 和 CMOS ADC 的组合可达到一个噪声<0.8nV/rtHz,功耗<150mW/CH 的出色模拟前端解决方案。这种组合不仅需要专用半导体工艺,还需要先进的封装技术。下图所示为一个模拟前端解决方案,在同一封装中具有两个芯片。实际上,还可以集成两个以上的芯片和多个无源元件。此外,多芯片模块(MCM)可为系统设计提供更大的灵活性。例如,若有更新的 ADC 或 VCA 解决方案,它可取代旧的 AFE 解决方案中的一个解决方案,且仍然保持管脚到管脚的兼容性,以获得更佳性能。
         
                                             
        图 13. 多芯片模块封装
         
        在过去十年中,超声波 AFE 的工艺技术从 0.5um 移至 90nm,从 CMOS 仅移至 BiCMOS 和 CMOS,从单芯片移至封装中具有无源元件的多个芯片。如图所示所有这些技术大大降低了功耗,提升了性能并缩小了芯片尺寸
         
                                                       
        图 14.AFE 集成度的发展

         
7. 超声模拟电路的主要参数
        超声信号有其自身的特点。正如我们在前面部分所讨论的那样,系统中经常会观察到超过 100dB 的动态范围。低频音频电路、高频数字电路、低噪声放大器、低噪声时钟电路存在于同一系统、同一电路板或甚至同一芯片上。AFE 设计和系统设计必须解决这些挑战。
         
过载恢复
        过载信号通常是指高压发射脉冲通过高压收发开关(T/R switch)的泄漏大信号或者是强回波信号。若 AFE 设计中未考虑过载恢复,它们会降低 LNA、PGA、ADC 和 CW 电路的瞬态响应性能。模拟设计人员面临着在有限功耗预算条件下,在大动态范围内实现瞬时恢复响应且响应性能一致的挑战。作为一种较为常见的设计方案,应首先在高压收发开关设计中应用足够的限流限压技术,这可消除对模拟前端的第一级即低噪声放大器的过载影响。在 LNA 设计中,钳位二极管通常可防止 LNA 进一步饱和。  
         
        分析两种常见的过载情况。第一种是由于高压收发开关导通,考虑到超声成像的死区时间通常在 3 到 5us 左右,因此超声模拟前端的过载恢复时间必须达到微秒量级,。目前由基 MOSFET 的高压收发开关处理,仅允许<< 1Vpp 发射泄漏直通;而基于二极管桥构成的高压收发开关,其泄漏电压可达 2Vpp。因此大多数 AFE 设计为可处理~2Vpp 过载信号,以满足各种收发开关的性能。另一种过载情况是由于来自血管壁的大反射信号,超声模拟前端必须立即恢复,以检测血液中的小回波。第二种情况在多普勒应用中极其常见,其性能决定了血流检测灵敏度和准确度。下图所示为模拟血管壁强回波,然后是来自血液的小信号的反应。下述信号具有 60dB 的动态范围,即 5 周期 250mVpp 信号和 5MHz 时的 5 周期 250uVpp 信号;小信号配置为具有 0°或 180°相移。下图所示为超声模拟前端的响应以及 0°和 180°响应之间的差异,即类似于多普勒应用中的相位检测相似。小信号和相位差的提取保证了多普勒应用中的良好性能。
         
                                                   
        图 15. 过载恢复(a)输入信号;(b)输出信号  
         
        除快速过载响应和精确的相位检测之外,多次过载恢复响应的一致性对于频谱多普勒和彩色多普勒应用也是至关重要的。一致的过载恢复可减少系统中的频谱噪声或彩色噪声。大家可通过比较来自多个信号的过载响应差异来评估一致性。
         
        此外,谐波成像是大多数系统中的标准配置。脉冲反转成像被广泛使用。因此系统 AFE 保证对正脉冲和负脉冲有对称的过载响应。最后,在超声系统中经常叠加多个不同类型的图像便于诊断,如双工模式、即 B 模式转换到多普勒模式,甚至三共模式。以。各个工作模式使用到不同的发射电压和占空比的发射波形。因此,AFE 需要快速响应两个或多个图像线内的不同过载信号。当快速切换图像模式时,不同的过载信号不应影响 AFE 过载性能的一致性。   
         
多普勒应用中的信号和噪声调制
        超声系统是复杂的混合信号系统,具有各类数字和模拟电路。数字信号和时钟信号会干扰系统级或芯片级的模拟信号。另外诸如晶体管和二极管之类的非线性元件可调制噪声,也会干扰 RF 信号。
         
        在超声多普勒应用中,系统中的调制效应会影响图像质量和灵敏度。多普勒信号频率范围从 20Hz 到>50Khz 不等。同时,多个系统的时序信号也在此范围内,如帧时钟、成像线时钟等。这些噪声信号可通过接地、电源和控制管脚进入芯片。研究芯片级的调制效应,如:(PSMR)电源调制比就很重要。具有一定频率和幅度的噪声信号可施加在电源管脚上。若存在调制效应,则可找到边带信号。PSMR 表示为载波和边带信号之间的幅度比,如下所示:  
         
                                               
        图 16.PSMR(a)和 IMD3(b)描述
         
        除 PSMR 外,三阶交调干扰(IMD3)是衡量混合信号 IC 性能的关键参数。同时,在超声应用中,用于 IMD3 测量的输入信号具有不同的幅度,其分别代表来自静态组织大回波和来自流动血液的多普勒小信号,其幅度差可达 20 到 30dB 左右。系统设计人员可使用 IMD3 来估计由多普勒镜像频率信号所产生的伪像。多普勒频谱显示中常用 40 到 50dB 的动态范围。因此,优于 50dBc 的 IMD3 不应影响系统性能。   
         
连续波(CW)多普勒参数
        作为中高端系统的关键功能,连续波多普勒已开始成为便携式系统的标配。与 TGC 路径相比,连续波多普勒路径具有处理更大动态范围和更低相位噪声的优点。此外,由于具备这些特点,连续波多普勒波束合成通常在模拟域中实现。多种波束合成方法都在超声系统中应用,包括无源延迟线、有源混频器和无源混频器。在过去几年中,基于混频器的连续波多普勒结构凭借体积小、易于实现且支持多个 CW 频率的良好灵活性而逐渐占据主导地位。此外,连续波多普勒波束合成器已集成在与 TGC 路径相同的芯片上。此外无源混频器不仅降低了功耗和噪声,也满足了连续波多普勒的处理要求,如宽动态范围、低相位噪声、精确 I/Q 通道增益和相位匹配等。
         
        简化的连续波多普勒路径框图如下所示。整个 CW 路径包括 LNA、电压电流转换器、基于开关电路的无源混频器、带低通滤波器的加法器和时钟电路。大多数模块包括性能严格对称的同相和正交通道,以实现良好的镜像频率抑制和波束合成精度。
         
                                                  
        图 17.CW 的简化框图
         
        下述图解和方程式描述了混频器操作的原理。
         
                                                                     
        图 18. 混频器操作的框图
         
        方程式中,Vi(t)、Vo(t) 和 LO(t)分别是混频器的输入、输出和本地振荡器信号。Vi(t)包括高次谐波;LO(t)代表方波,其包含奇数谐波分量,如下式所示:
         
                                        
    
        根据方程式,来自 LO(t) 的 3 阶和 5 阶谐波可与 Vi(t)中这些频带中的 3 阶和 5 阶谐波或宽频噪声相调制。因此,混频器的噪声性能会降低。为了避免这种不良影响,在 LNA 输出或混频器时钟输入上都需要谐波抑制电路,以实现更佳的噪声系数。根据以上方程式,混频器的转换损耗约为 20log2/π 也就是 -4dB 左右。
         
        优于 -46dBc 镜像频率抑制是 CW 成像中的期望参数。CW I/Q 通道匹配也可有助于镜像频率分量。文献表明,0.25°的 I/Q 相位误差可得到 -53dBc 抑制;且 0.05dB 的 I/Q 增益误差可得出 -50dBc 抑制。它们是 CW 路径的设计目标。因此,CW I/Q 路径需要严格的增益和相位匹配。低容差电阻(0.1%)通常用于基于运算放大器的有源滤波器。  
         
        典型的 CW 多普勒移频处在 100Hz 到 20KHz 之间。由于混频特性,CW 信号路径的相位噪声主导低血流速度。因此,大多数 AFE 以载波频偏为 1KHz 时的 CW 相位噪声作为主要性能指标。
         
        最后,CW 路径的动态范围基于输入参考噪声和最大输入信号:
         
                                          
         
        为实现良好的 CW 性能,需要>160dBFS/Hz 的发射器和接收器电路动态范围。
         
8. 总结
        超声成像是一种安全的医学成像模式,具有很大潜力,。越来越多就地检查的床旁应用需要低功耗、低噪声和紧凑的系统。为了充分发挥超声信号的优点,必须选择合适的工艺以实现低功耗、低噪声和小尺寸的目的。BiCMOS 工艺适用于低噪声放大器设计,具有超低 1/f 噪声、宽带宽和良好的功耗 / 噪声效率;而 CMOS 工艺在低功耗下实现了高数字密度。两者结合使用先进的封装技术,可提供最先进的模拟前端解决方案。为达到所需的超声波参数,如快速一致的过载恢复、低 IMD3 和 PSMR、精确的 I/Q 匹配、连续波多普勒混频器中奇次谐波抑制等,需要考虑芯片中的各个参数已到达设计的综合优化。
      

最后

以上就是谦让小丸子为你收集整理的转-超声系统的信号链设计注意事项的全部内容,希望文章能够帮你解决转-超声系统的信号链设计注意事项所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(63)

评论列表共有 0 条评论

立即
投稿
返回
顶部