网站地图
-
随着多核计算时代的来临，人们对信息处理、传输速度的要求不断提高，电互连技术由于串扰、延迟和功耗等缺陷成为系统发展的瓶颈。采用光互连取代电互连可以有效解决这一难题[1]。在光互连的具体实施方案中，硅基光互连以其无可比拟的成本和技术优势成为首选，既能发挥光互连抗干扰、功耗低、速度快、带宽大等优点，又能充分发挥微电子低成本、工艺成熟、高成品率、高集成密度等优势，必将推动新一代光通信系统和高性能计算机的发展，有着广阔的市场应用前景。
光互连的一大优势是可以采用信息复用技术提高带宽密度，增加通信容量。其中波分复用(wavelength-division multiplexing, WDM)是应用最广泛的复用技术，而复用/解复用器是其中最核心的功能器件，是提高系统整体性能的关键。目前硅基光互连中，实现波分复用/解复用功能的波长滤波器件主要有刻蚀衍射光栅(etching diffraction grating, EDG)[2-4]、阵列波导光栅(arrayed waveguide gratings, AWG)[5-7]、微环谐振器(micro-ring resonator, MRR)[8-9]和马赫-曾德尔干涉仪(Mach-Zehnder interferometer, MZI)[10-11]4种。其中EDG和AWG是最典型的平面集成型光学器件，它们均具有集成度高的优点，适合于实现窄通道间隔和大通道数的波分复用器[12]。相比于AWG器件，EDG具有面积小、齿面间距小且数量多、单边输入输出等优点，而且具有成像质量高、集成密度大、通道频谱精度高和易于封装等优势。
本文中在绝缘体上硅(silicon-on-insulator, SOI)上设计一款四通道、20nm通道间隔的EDG，通过在闪耀光栅反射面引入布喇格反射光栅来降低器件插入损耗，在入射波导处引入多模干涉耦合器(multi-mode interference, MMI)实现通道频谱平坦化设计，并通过数值仿真方法对器件最终的传输光谱进行仿真计算。整个器件的制作过程与商用成熟集成电路工艺平台兼容，可以精确地大规模制作。近年来，基于1310nm波段的稀疏波分复用系统由于其在城域网的优势成为研究热点，目前该波段硅基光收发模块[13-14]和器件[15]的研究均有报道，然而用于该波段的硅基EDG却鲜有报道。因此，本文中设计工作于1310nm波段的EDG实现波分复用功能，以满足未来城域光网络的发展需求。
-
EDG设计中通常采用凹面闪耀光栅，利用相邻两个齿面光程差的差异实现色散与聚焦功能。EDG设计通常基于罗兰圆结构，其工作原理如图 1所示。输入输出波导(waveguide, WG)分别位于罗兰圆上，延长线相交于切点P，凹面闪耀光栅位于大圆上并且与小圆相切于点P。对于凹面闪耀光栅，衍射方程为：
$ \sin {\theta _{\rm{i}}} + \sin {\theta _{\rm{d}}} = \frac{{m\lambda }}{{{n_{{\rm{eff}}}}d}} $
(1) 
Figure 1. Schematic of EDG based on Rowland circle
式中，θi和θd分别是入射角和出射角，λ, m, d和neff分别是入射波长、衍射级数、光栅周期和自由传输区有效折射率。
定义线色散度(linear dispersion, LD)dLD为单位波长偏移对应的输出波导位置偏移量，即Δs/Δλ，通过(1)式两边取微分计算角色散关系，可得dLD的表达式为[2]：
$ \begin{array}{*{20}{l}} {{d_{{\rm{LD}}}} = 2R(\frac{{{\rm{d}}{\theta _{\rm{d}}}}}{{{\rm{d}}\lambda }}) = \frac{{2R}}{{\cos {\theta _{\rm{d}}}}}\frac{m}{d}\frac{{{n_{\rm{g}}}}}{{{n_{{\rm{eff}}}}^2}}{\rm{ = }}}\\ {\;\;\;\;\;\;\;\;\frac{{2R}}{\lambda }\frac{{(\sin{\theta _{\rm{i}}} + \sin {\theta _{\rm{d}}})}}{{\cos {\theta _{\rm{d}}}}}\frac{{{n_{\rm{g}}}}}{{{n_{{\rm{eff}}}}}}} \end{array} $
(2) 自由光谱范围(free spectral rang, FSR)ΔλFSR可以表示为[2]：
$ \Delta {\lambda _{{\rm{FSR}}}} = \frac{\lambda }{m}{\left[ {1 - \frac{{m + 1}}{m}\left( {1 - \frac{{{n_{\rm{g}}}}}{{{n_{{\rm{eff}}}}}}} \right)} \right]^{ - 1}} $
(3) 式中，ng是自由传输区的群折射率, ng＝neff－λ·dneff/dλ。
根据以上几式结合设计要求，可以初步确定EDG的基本参量。由(1)式确定的光栅周期d是指各光栅齿面中点在P点切线上的投影间距。这种投影等间距的设计方法忽略了光传输路径上的高阶像差，这些高阶像差在远离切点的光栅齿面上会逐渐增加以至于不能忽略，否则将会影响成像质量。因此，作者采用一点法对各光栅齿面中点位置进行设计[16]。
要满足指定波长零像差输出条件，则要求光信号从罗兰圆上的输入波导处出发，经过任意光栅齿面中点，再汇聚到罗兰圆上的特定输出波导位置处所间隔的光程差为该波长的整数倍。在图 1中，令自由传输区有效折射率为neff，选定波长为λ0时的零像差输出位置为O点，并令P点为光栅第0个槽面的中点，各个光栅齿面中心点P′(x, y)满足如下关系：
$ \begin{array}{l} \sqrt {{{({r_{\rm{i}}}\cos {\theta _{\rm{i}}} - x)}^2} + {{({r_{\rm{i}}}\sin {\theta _{\rm{i}}} - y)}^2}} + \\ \sqrt {{{({r_{\rm{d}}}\cos {\theta _{\rm{d}}} - x)}^2} + {{({r_{\rm{d}}}\sin {\theta _{\rm{d}}} - y)}^2}} - \\ \;\;\;\;\;\;\;\;\;{r_{\rm{d}}} - {r_{\rm{i}}} = \frac{{Km{\lambda _0}}}{{{n_{{\rm{eff}}}}}} \end{array} $
(4) 式中，K为所计算光栅齿面中心点P′与散射中心点P所间隔的光栅周期数，其取值范围为[-M，N]，ri为入射波导中点I到P点之间的距离，rd为零像差输出位置O到P点之间的距离。M和N的大小由输入波导中的出射光衍射传输到光栅处的能量分布决定，它们的选择应该使得衍射光栅伸展到能覆盖绝大部分入射能量，一般应该达到99%以上。
在(4)式中需要解得P′点横纵坐标参量，即衍射光栅各中心点位置坐标。由于P′点位于与罗兰圆相切且半径为2R的大圆上，因此有(5)式成立，与(4)式联立可求得各光栅齿面中点P′的坐标:
$ \sqrt {{{\left( {x - 2R} \right)}^2} + {y^2}} = 2R $
(5) 确定各个光栅齿面中点后，根据输入输出波导位置可以确定每个光栅面的方向。光栅面的反射结构最简单方式就是采用端面反射，如图 2a所示，但这种结构反射效率较差，仅有35%左右; 在此采用全反射光栅结构[17]，如图 2b所示。因此针对1310nm波段进行参量优化，最终反射光栅参量取光栅周期Λ=0.52μm，刻蚀宽度w=0.2μm，仿真的反射光谱如图 3所示。显示其在1270nm~1350nm工作波长范围内反射效率大于0.85，到此整个EDG的参量设计全部完成。
Figure 2. Schematic of SOI grating
Figure 3. Relationship between reflectivity and wavelength of Bragg reflector
作者在SOI上设计一款四通道20nm间隔的EDG，顶层硅厚度H=220nm，具体设计参量如表 1所示。其中入射波导和出射波导采用130nm浅刻蚀脊形波导，同时采用梯形波导将单模波导展宽至2μm，以减小矩形波导与平板波导之间的模式失配，增加耦合效率。
Table 1. Design parameters of 4-channel EDG
λ0 Δλ N θi θd m dLD H 1301nm 20nm 4 42° 36° 10 250 220nm -
设计好EDG基本结构参量之后，需要对光在器件中的传输过程进行仿真计算，以检验所设计器件是否能够满足实际应用要求，同时也可以根据仿真结果对器件的结构和参量进行优化改进。对于EDG工作原理，按照光束的场传输过程，可以将仿真过程分为以下4步。
(1) 利用2维基尔霍夫标量衍射理论[2]计算入射光传输到各光栅齿面上的光场分布，如下所示:
$ \begin{array}{l} {E_{{\rm{inc}}}}\left( {P'} \right) = \frac{1}{2}\sqrt {\frac{{{n_{{\rm{eff}}}}}}{\lambda }} \int_{{\rm{input}}} {{E_{{\rm{input}}}}\left( y \right)} \times \\ \;\;\;\;\;\;\;\frac{{{{\rm{e}}^{ - {\rm{j}}k{r_{\rm{i}}}}}}}{{\sqrt {{r_{\rm{i}}}} }}(1 + \cos \alpha ){\rm{d}}y \end{array} $
(6) 式中，Einc(P′)是各光栅齿面中点处的光场能量，k=2πneff/λ是平板传输区内的波数，Einput是入射波导电场分布。
在低折射率对比的材料系统中，一般采用高斯光束场分布近似代替波导电场分布，此时(6)式可以利用解析方式求解。但在高折射率差对比的材料系统中(例如SOI)，这种高斯近似还不能满足设计的精度要求，主要是由于这类波导电场分布在边界处不连续性比较严重，高斯曲线不能准确表示出这种边界不连续性，而且波导外部的倏逝波衰减与高斯曲线的指数衰减差异较大，因此入射光场需要采用实际的场分布表示。并且EDG是一个偏振敏感器件，作者设计其仅工作在TE模式下，在此利用商业电磁仿真软件采用本征模展开法(eigen-mode expansion, EME)计算入射波导在TE模式时的光传输特性，并得到入射端光场分布，进而提取出入射光场的1维电场分布，如图 4所示。
Figure 4. Electric field intensity distribution of input ridge waveguide and 1-D diagram
仿真计算得到各光栅中点处的场分布如图 5a所示。可以看到入射波导的衍射远场分布情况，由于是正入射，切点P处光强最大，两边逐渐减小，包络外侧有旁瓣。
Figure 5. a—incident electric field at the center point of each grating facet for λ=1.301μm b—output electric field distribution on Rowland circle c—transmission spectrum of 4-channel EDG
(2) 仿真光场在闪耀光栅处发生衍射的过程，在此采用平面波近似法[18]仿真，计算得到衍射后光栅面上电场分布为：
$ {E_{\rm{g}}}^\prime \left( {x', z'} \right) = R'{E_{\rm{g}}}\left( {x', z'} \right) $
(7) 式中，R′为光栅反射率，Eg为入射前光栅面上电场分布，Eg′为衍射后光栅面上电场分布，x′和z′对应的是各个衍射光栅面的坐标。
(3) 仿真光场从衍射光栅到输出面在平板波导中传输的过程，在此仍然采用基尔霍夫标量衍射理论计算。位于罗兰圆上输出点P″处的光场分布Eout等于所有光栅反射回的光场叠加总和。
$ \begin{array}{l} {E_{{\rm{out}}}}\left( {P''} \right) = \frac{\eta }{2}\sqrt {\frac{{{n_{{\rm{eff}}}}}}{\lambda }} \sum {\int_{ - D/2}^{{\rm{ + }}D/2} {{E_{{\rm{inc}}}}\left( {y'} \right) \times } } \\ \frac{{{{\rm{e}}^{ - {\rm{j}}k{r_{\rm{d}}}}}}}{{\sqrt {{r_{\rm{d}}}} }}(cos{\beta _{\rm{i}}} + \cos {\beta _{\rm{d}}}){\rm{d}}y' \end{array} $
(8) 式中，βi和βd是各点光场的入射角和反射角，D是闪耀光栅宽度，η是每个光栅面的反射系数，考虑到工艺偏差与传输损耗，仿真中取η=0.8。得到罗兰圆上光场分布如图 5b所示。可以看到，4个波长对应的光场按照角度均匀聚焦在罗兰圆上的输出波导位置。
(4) 仿真光场耦合进入输出波导的过程，进而得到器件传输光谱曲线，采用叠加积分来仿真计算:
$ I(\lambda ) = \frac{{{{\left| {\int {{E_{{\rm{out}}}}(\lambda , y''){E_{{\rm{output}}}}^*(\lambda , y''){\rm{d}}y''} } \right|}^2}}}{{\int {{{\left| {{E_{{\rm{output}}}}(\lambda , y'')} \right|}^2}{\rm{d}}y'' \cdot \int {{{\left| {{E_{{\rm{input}}}}(\lambda , y)} \right|}^2}{\rm{d}}y} } }} $
(9) 式中，I表示归一化输出光场强度，Eout表示衍射光在罗兰圆上相干叠加后的电场分布，Eoutput表示出射波导电场分布，y″对应的是输出波导面的坐标，上标*表示共轭。
传输光谱仿真结果如图 5c所示。该EDG的插入损耗约为2dB，通道串扰小于-45dB。其中插入损耗主要来源有闪耀光栅、布喇格反射光栅和输出耦合等。首先，由于光束进入自由传输区后会发生衍射发散，为了减小器件面积，闪耀光栅尺寸有限，导致两侧会有部分光能量损耗，如图 5a所示，但其非常小可以忽略。其次，由于位置所限，布喇格反射光栅设计为4个光栅周期，反射效率约为85%，这将引入部分能量损耗。最后，出射面上由于其它干涉级次的存在使得输出波导处聚焦光场有部分能量损耗。这些综合因素导致了整个器件的插入损耗。
-
波分复用器在大容量、超高速波分复用系统中处于核心地位，器件的通道带宽是影响其大规模应用的重要指标之一。按照频谱响应形状可以将其分为高斯型和平坦型。对于高斯型，由于1dB带宽较小，外界因素(环境温度、中心波长和偏振态等漂移)将影响器件在实际系统中的应用。相比之下，频谱平坦型器件具有平坦的输出频谱特性，这可以提高系统的稳定性，增加传输距离和容量，降低系统工作成本。
对于实现EDG频谱平坦化的方法，主要可以通过以下3种结构优化设计实现：衍射光栅、输入波导和输出波导。对于衍射光栅和输入波导的优化，主要是改变了出射面的成像场分布，而输出波导的优化则是改变了输出波导的耦合特性，这些方法都可以通过设计器件的版图来实现。
以上3种方法中，改变输入波导的方法最容易实现，只需在特别小的区域内进行设计调整，而且具有较大的工艺容差。另外，级联不同输入结构时，设计参量可以具有更大的灵活性。考虑到工艺容差和性能方面的要求，作者选择改变输入波导结构对前面的EDG实现通道频谱平坦化设计，在入射波导端引入级联的梯形绝热波导和MMI结构。
本次优化设计中，取MMI长度为器件将一束光均分为两束光自映像所对应的长度，即LMMI=3Lπ/8，通过分别优化MMI输入波导宽度Wi、多模区宽度WMMI和EDG输出波导宽度Wo来实现频谱平坦化设计，如图 6a所示。利用EME方法仿真MMI光传输特性，得到MMI输出场分布，如图 6b所示。对于给定的多模区宽度，MMI拍长随之确定，进而确定LMMI参量，不同的Wi会影响频谱平坦化效果，这主要受MMI的自映像效应影响。此时，EDG的最终结构随之确定，器件整体结构示意图如图 7所示。包含4路不同波长的光信号由基于MMI结构的输入波导进入平板传输区，衍射传输到闪耀光栅面被布喇格光栅反射后聚焦于各个输出波导处，不同波长光信号由相应的输出波导输出，实现波分复用功能。
Figure 6. a—EDG input/output waveguide based on MMI b—electric field intensity of input waveguide based on the cascaded taper and MMI c—1-D electric field intensity of end face
Figure 7. EDG structure based on MMI and Bragg reflector
针对不同宽度的MMI，通过优化输入输出波导参量得到其平坦化频谱输出，结果如图 8a所示。可以看到, 不同WMMI对应不同的平坦化效果，性能差异主要表现在插入损耗和1dB带宽上，其中器件插入损耗与1dB带宽性能成负相关，因此器件参量的确定要在两者之间做折中考虑。通过最终参量优化，选择WMMI=4μm，Wi=3μm，Wo=2μm作为器件频谱平坦化设计参量。仿真结果显示，中心通道1dB带宽大于12nm，插入损耗约为5dB，通道串扰约为40dB。由于MMI具有色散特性，对于非中心通道，首先选择相应的通道波长仿真MMI输出场，再通过输出场计算相应通道传输光谱，结果如图 8b所示。可以看到，器件整体频谱平坦化效果较好，但非中心通道光谱存在轻微倾斜，可以通过调整各通道输出波导参量进行改进。
Figure 8. a—flat-top spectral response with different WMMI (Wo=2μm) b—transmission spectrum of 4-channel flat-top EDG (WMMI=4μm, Wi=3μm)
EDG的一大优势是自由传输区不受刻蚀工艺的影响，且输入输出波导具有较大的制作容差，因此工艺偏差仅会影响闪耀光栅的质量，本次仿真中的参量选取已充分考虑到布喇格光栅的工艺偏差。这为器件设计参量的可靠性提供了保障，可以保证器件的实际应用价值。
-
设计了一款用于稀疏波分复用系统中的刻蚀衍射光栅来实现波分复用/解复用功能。通过在闪耀光栅反射面引入布喇格反射光栅将光栅反射效率提高到85%，减小器件插入损耗至2dB; 在入射波导处引入多模干涉器件实现了各通道平坦化频谱输出，1dB带宽达到12nm。这些为硅基片上刻蚀衍射光栅设计和优化提供了理论基础和设计指导。
硅基片上刻蚀衍射光栅及其平坦化设计
Etching diffraction grating of silicon substrate and design of flatten
-
摘要: 刻蚀衍射光栅（EDG）作为实现波分复用功能的关键器件，对于片上光互连的实现至关重要。为了实现1310nm波段通道间隔为20nm的硅基EDG，采用了基尔霍夫标量衍射理论仿真方法进行理论设计和仿真验证，通过在闪耀光栅反射面引入布喇格反射光栅来提高反射效率、降低器件插入损耗，并在入射波导处引入多模干涉耦合器以实现通道频谱平坦化设计。结果表明，闪耀光栅反射面的反射效率由35%提高到了85%，1dB带宽达到12nm。这对于提高系统稳定性、增大传输距离和容量、降低系统成本具有显著作用，能够满足光互连系统的实际应用需求。Abstract: Etching diffraction grating (EDG), one of the most critical components, can achieve wavelength division multiplex (WDM) function and realize on-chip optical interconnection. In order to realize 4-channel EDG with the wavelength spacing of 20nm at 1310nm wavelength, Kirchhoff scalar diffraction theory was used for theory design and simulation verification. To further improve the reflection efficiency and decrease the insertion loss, Bragg reflection gratings were designed to replace the normal etched facets. Multimode interference (MMI) coupler was also introduced at input waveguide for flat-top frequency response. The simulation results demonstrate the reflection efficiency of grating reflective surface is up to 85% and 1dB bandwidth is up to 12nm. The designed EDG has a significant effect on improving system stability, increasing transmission distance and capacity, and reducing system cost. The design can meet the practical application requirements of optical interconnection system.
-
Figure 4. Electric field intensity distribution of input ridge waveguide and 1-D diagram
Figure 5. a—incident electric field at the center point of each grating facet for λ=1.301μm b—output electric field distribution on Rowland circle c—transmission spectrum of 4-channel EDG
Figure 6. a—EDG input/output waveguide based on MMI b—electric field intensity of input waveguide based on the cascaded taper and MMI c—1-D electric field intensity of end face
Figure 8. a—flat-top spectral response with different WMMI (Wo=2μm) b—transmission spectrum of 4-channel flat-top EDG (WMMI=4μm, Wi=3μm)
Table 1. Design parameters of 4-channel EDG
λ0 Δλ N θi θd m dLD H 1301nm 20nm 4 42° 36° 10 250 220nm 
下载: 导出CSV
-
[1] ZHANG B, SONG Sh H, WANG H J, et al. A novel micro-ring-based 4×4 non-blocking silicon optical router[J]. Laser Technology, 2013, 37(6):731-735 (in Chinese). [2] BROUCKAERT J, BOGAERTS W, DUMON P, et al. Planar concave grating demultiplexer fabricated on a nanophotonic silicon-on-insulator platform [J]. Journal of Lightwave Technology, 2007, 25(5): 1053-1060. [3] CHOWDHURY D. Design of low-loss and polarization-insensitive reflection grating-based planar demultiplexers [J]. IEEE Journal of Selected Topics in Quantum Electronics, 2000, 6(2): 233-239. doi: 10.1109/2944.847758 [4] FENG D, QIAN W, LIANG H, et al. Fabrication intensive echelle grating in silicon-on-insulator platform [J]. IEEE Photonics Technology Letters, 2011, 23(5): 284-287. [5] ZOU J, JIANG X X, XIA X, et al. Ultra-compact birefringence-compensated arrayed waveguide grating triplexer based on silicon-on-insulator [J]. Journal of Lightwave Technology, 2013, 31(12): 1935-1940. doi: 10.1109/JLT.2013.2261857 [6] WANG J, SHENG Z, LI L, et al. Low-loss and low-crosstalk 8×8 silicon nanowire AWG routers fabricated with CMOS technology [J]. Optics Express, 2014, 22(8): 9395-9403. doi: 10.1364/OE.22.009395 [7] YE T, FU Y, QIAO L, et al. Low-crosstalk Si arrayed waveguide grating with parabolic tapers [J]. Optics Express, 2014, 22(26): 31899-31906. doi: 10.1364/OE.22.031899 [8] ZHENG X Z, SHUBIN I, LI G L, et al. A tunable 1×4 silicon CMOS photonic wavelength multiplexer/demultiplexer for dense optical interconnects [J]. Optics Express, 2010, 18(13): 5151-5160. [9] FANG Q, PHANG Y T, TAN C W, et al. Multi-channel silicon photonic receiver based on ring-resonators [J]. Optics Express, 2010, 18(13): 13510-13515. doi: 10.1364/OE.18.013510 [10] HORST F, GREEN W M J, ASSEFA S, et al. Cascaded Mach-Zehnder wavelength filters in silicon photonics for low loss and flat pass-band WDM (de-)multiplexing [J]. Optics Express, 2013, 21(10): 11652-11658. doi: 10.1364/OE.21.011652 [11] JEONG S, TANAKA S, AKIYAMA T, et al. Flat-topped and low loss silicon-nanowire-type optical MUX/DeMUX employing multi-stage microring resonator assisted delayed Mach-Zehnder interferometers [J]. Optics Express, 2012, 20(23): 26000-26011. doi: 10.1364/OE.20.026000 [12] AN J M, ZHANG J S, WANG Y, et al. Study on wavelength division multiplexer for silicon photonics [J]. Laser & Optoelectronics Progress, 2014, 51(11):110006 (in Chinese). [13] DENOYER G, COLE C, SANTIPO A, et al. Hybrid silicon photonic circuits and transceiver for 50Gb/s NRZ transmission over single-mode fiber [J]. Journal of Lightwave Technology, 2015, 33(6): 1247-1254. doi: 10.1109/JLT.2015.2397315 [14] TORRES-FERRERA P, PACHECO-RAMÍREZ L, GUTIÉRREZ-CASTREJÓN R, et al. Next-generation 400Gb/s ethernet PMD over SMF at 1310nm via DD-OFDM with electro-absorption modulator-based transmitters [C]//2015 7th IEEE Latin-American Conference on Communications. New York, USA: IEEE, 2015: 1-6. [15] STRESHINSKY M, DING R, LIU Y, et al. Low power 50Gb/s silicon traveling wave Mach-Zehnder modulator near 1300nm [J]. Optics Express, 2013, 21(25): 30350-30357. doi: 10.1364/OE.21.030350 [16] SHENG Z Y, HE S L, HE J J. 1-Stigma and 2-stigma method for the design of etched diffract ion grating [J]. Journal of Optoelectronics·Laser, 2001, 12(7): 671-674 (in Chinese). [17] BROUCKAERT J, BOGAERTS W, SELVARAJA S, et al. Planar concave grating demultiplexer with high reflective Bragg reflector facets [J]. IEEE Photonics Technology Letters, 2008, 20(4): 309-311. doi: 10.1109/LPT.2007.915585 [18] SHENG Zh Y, HE S L, HE J J. Simulation for etching diffraction grating by use of scalar diffraction theory [J]. Opto-Electronic Engineering, 2001, 28(6): 29-32(in Chinese). -
点击查看大图
图(8) / 表(1)
计量
- 文章访问数: 5424
- HTML全文浏览量: 3603
- PDF下载量: 16
- 被引次数: 0
