logo资料库

基于3×3迈克耳孙干涉仪的四路检测相位解调φ-OTDR.pdf

发布时间：2022-05-30 发布人：admin 分类：说明书 资料大小：5.50M 资料格式：pdf 举报 版权申诉
第1页 / 共6页
第2页 / 共6页
第3页 / 共6页
第4页 / 共6页
第5页 / 共6页
第6页 / 共6页
资料共6页，全文预览结束
    第３８卷　第８期 ２０１８年８月 光　学　学　报 ACTAOPTICASINICA Vol．３８,No．８ August,２０１８ 基于３×３迈克耳孙干涉仪的四路检测 相位解调ϕＧOTDR 宋牟平∗∗,尹聪∗,陆燕,朱伟基,吴媖,庄守望 浙江大学信息与电子工程学院,浙江 杭州３１００２７ 摘要　相位敏感光时域反射计(ϕＧOTDR)中使用相干性较好 的 窄 线 宽 光 源,因 传 感 光 纤 中 光 脉 冲 产 生 的 散 射 光 发 生干涉,故可利用干涉光的变化检测外界扰动,检测灵 敏 度 较 高.但 传 统 的 幅 度 检 测ϕＧOTDR 的 结 果 难 以 实 现 定 量检测,需要进行相位解调.基于３×３迈克耳孙干涉仪的相位解调系统是一 种 较 好 的 解 调 方 案,但 直 接 解 调 的 结 果存在空间分辨能力小于传统幅度检测ϕＧOTDR 的 问 题.采 用 四 路 检 测 方 法 来 保 持ϕＧOTDR 的 空 间 分 辨 能 力, 只需在系统结构中增加一路检测.同时考虑到窄线宽光源产生的脉冲内干涉对 解 调 造 成 的 不 利 影 响,在 提 出 四 路 检测相位解调ϕＧOTDR 的基础上,从瑞利后向散射光的角度对基于３×３迈克 耳 孙 干 涉 仪 的 相 位 解 调ϕＧOTDR 进 行分析,并通过实验进行验证.完成了对５km 传感光纤上 扰 动 信 号 的 定 量 检 测,检 测 线 性 度 为０．９９５６,解 调 相 位 的幅度达到３１．８５rad,解调系统的空间分辨率与传统的幅度检测ϕＧOTDR 相同. 关键词　光纤光学;分布式光纤传感;相位敏感光时域反射计;相位解调;定量检测;四路检测 中图分类号　 TP２１２　　　文献标识码　A doi:１０．３７８８/AOS２０１８３８．０８０６００１ FourＧChannelDetectingPhaseDemodulationϕＧOTDR Basedon３×３MichelsonInterferometer SongMuping∗∗ YinCong∗ LuYan ZhuWeiji WuYing ZhuangShouwang CollegeofInformationScienceandElectronicEngineering ZhejiangUniversity Hangzhou Zhejiang３１００２７ China Abstract　PhaseＧsensitiveopticaltimedomainreflectometry ϕＧOTDR usesnarrowlineＧwidthlaserwithhigh coherence敭Inthesensingfiber scatteredlightsfromopticalpulseinterfere andtheinterferencelightishighly sensitivetoexternaldisturbance敭ButthetraditionalamplitudeＧdetectionϕＧOTDRisnotquantitativedetection so thatphasedemodulationisnecessary敭Phasedemodulationsystembasedon３×３MichelsoninterferometeriscostＧ effective butisoflowerspatialresolutionthanthetraditionalamplitudeＧdetectionϕＧOTDR敭WeproposeafourＧ channeldetectionmethodtopreservethesamespatialresolutionasamplitudeＧdetectionϕＧOTDRbyaddingonlyone detectionchannel敭Buttheinterferenceofscatteredlightsmayhavecriticalinfluencetophasedemodulation敭AfourＧ channeldetectionsystemofphasedemodulationϕＧOTDRisanalyzedbasedon３×３Michelsoninterferometer from theaspectofRayleighbackscatteringlight敭Theexperimentshowsthatthedisturbancearoundthefiberisdetected quantitatively withalinearityof０敭９９５６ phasesensingrangeof３１敭８５rad andthespatialresolutionissameas thatofamplitudeＧdemodulationϕＧOTDR敭 Keywords　fiberoptics distributedopticalfibersensing phaseＧsensitiveopticaltimedomainreflectometry phase demodulation quantitativedetection fourＧchanneldetecting OCIScodes　０６０敭２３００ ０６０敭２３７０ ０６０敭４０８０ １　引　　言 近年来,分布式光纤传感技术被认为具有广泛 的应用前景,尤其 是 在 国 防、管 道 监 测、大 型 结 构 监 测等方面[１].相位敏感光时域反射计(ϕＧOTDR)因 其灵敏度较高、动态范围大、能够进行多点同时检测 　　收稿日期:２０１８Ｇ０１Ｇ０４;修回日期:２０１８Ｇ０２Ｇ０７;录用日期:２０１８Ｇ０３Ｇ２８ 　 ∗EＧmail:yincong１９９２０３１１＠qq．com;∗∗EＧmail:songmp＠zju．edu．cn ０８０６００１Ｇ１ 等优点而受到 广 大 研 究 人 员 的 关 注[２].目 前,人 们 在提高ϕＧOTDR 传感距离、空间分辨率以及信噪比 等方面进行了深入 研 究,同 时 结 合 光 纤 光 栅 等 辅 助 手段改良其性能,为ϕＧOTDR 的实际应用做出了诸 多努力[３Ｇ５].ϕＧOTDR 利 用 高 相 干 窄 线 宽 激 光 器 为 光源,传感光纤中光 脉 冲 产 生 的 后 向 散 射 光 相 互 干
    光　　　学　　　学　　　报 也有很多研究针对各种检测方法的解调算法进行改 检测方法通过使用频率漂移很小的激光器进行散射 涉,利用干涉光的变化对外界扰动进行检测,具有很 高的灵敏度. 传统的ϕＧOTDR 是 基 于 后 向 瑞 利 散 射 光 幅 度 检测的分布式光纤 传 感 器,但 扰 动 引 起 散 射 光 的 幅 度变化与扰动的强 弱 不 是 线 性 关 系,因 此 难 以 对 扰 动进行定量检 测[６Ｇ７].同 时 幅 度 检 测 在 扰 动 较 强 时 还会产生谐 波,影 响 对 外 界 扰 动 情 况 的 正 常 判 断. 而散射光的相位与 扰 动 的 强 弱 变 化 呈 线 性 关 系,因 此相位解调可以对扰动进行定量测量.已有的外差 光与本征光的相干检测,从而实现相位解调[８];零差 检测方法通过散射 光 的 自 相 干 进 行 解 调,包 括 相 位 生成载波方法及３×３迈克耳孙干涉仪的方法[９Ｇ１０]; 进,以减小解调过程带来的误差以及失真[１１].国内 的王毅等[１２]引进并研究了基于３×３耦合器的解调 方法,指出了该方法的优势与应用价值.基于３×３ 耦合器的迈克耳孙干涉仪,能够对光纤上两处散射光 的相位差进行解调,从而实现对传感光纤上扰动情况 的检测.该方法所用设备结构简单,成本较低,但由 于解调相位是光纤上两处散射光的相位差,需要进一 步处理解调结果才能保持原有的空间分辨率. 本文提出基于３×３ 迈 克 耳 孙 干 涉 仪 的 四 路 检 测相位解调ϕＧOTDR,在保持ϕＧOTDR 系统空间分 辨率和灵敏 度 的 同 时,能 够 对 扰 动 进 行 定 量 检 测. 调结果的影响,通 过 实 验 对 分 析 结 果 进 行 验 证.分 析了干涉仪非平衡 长 度 对 系 统 性 能 的 影 响,并 选 取 适当的非平衡长度实现了对５km 光纤上扰动的定 量检测. ２　原　　理 基于３×３迈克 耳 孙 干 涉 仪 的 四 路 检 测 相 位 解 调ϕＧOTDR 传感原理如图１所示,通过对光纤上两 定量扰动,并通过直 接 检 测 散 射 光 来 对 扰 动 进 行 空 间定位.当有扰动 存 在 时,在 扰 动 位 置 的 散 射 光 幅 度会发生变化,扰动 位 置 及 其 后 向 散 射 光 相 位 也 会 发生变化,从而引起扰动前后散射光的相位差变化, 信息. ２．１　光脉冲内散射光 ϕＧOTDR 中光 脉 冲 在 传 感 光 纤 中 传 播 发 生 散 射,其中后向散射光 携 带 着 光 纤 上 的 扰 动 带 来 的 幅 然后从散射光的角度推导分析脉冲内干涉现象对解 同时检测出幅度变化和相位差即可得到扰动的相关 个相隔一定距离的后向散射光的相位差进行解调来 度和相位变化等 信 息,被 探 测 器 接 收.散 射 模 型 可 以用图２表示. 图１ 基于３×３迈克耳孙干涉仪的四路检测相位 解调ϕＧOTDR 传感原理 Fig．１ SensingprincipleoffourＧchanneldetectingphase demodulationϕＧOTDRbasedon３×３ Michelsoninterferometer 图２ 光脉冲散射模型 Fig．２ Scatteringmodelofopticalpulse i＝１ 涉后的光场可以表示为 Eiexp(jϕi)＝Eexp(jϕ), 光脉冲产生的后向散射光相遇时发生干涉,干 E＝∑n (１) 式中:n为脉冲内散射单元个数;Ei 和ϕi 为第i个 散射点的幅度和相位;E 和ϕ为干涉 后 的 幅 度 和 相 位.当有扰动作用于光纤时,E 和ϕ都 会 随 扰 动 的 变化而变化,ϕＧOTDR 幅度检测 就 是 检 测 扰 动 引 起 的E 的变化. ２．２　散射光相位差解调 基于３×３ 迈 克 耳 孙 干 涉 仪 的 相 位 解 调 ϕＧ OTDR 是对光 纤 上 相 隔 一 定 距 离 的 两 个 后 向 散 射 光之间的相位差进行解调,简化模型如图３所示. 图３ 相位差模型 Fig．３ Modelofphasedifference 向散射光进入３×３耦合器时的光场可以表示为 根据(１)式,没 有 扰 动 时 脉 冲 在 A、B 两 处 的 后 (２) (３) 干涉后从３×３耦合器输出,最终探测器检测到的三 EA ＝EA０exp(jϕA), EB ＝EB０exp(jϕB), 路信号可以表示为 ０８０６００１Ｇ２
    光　　　学　　　学　　　报 l＝lp＋l０, 号平均衰 落 周 期 加 脉 冲 长 度 作 为 干 涉 仪 非 平 衡 长 距离的分布是具有接近但不完全一致的周期性的瑞 一定范围内解调结 果 不 完 全 失 真,就 可 以 确 定 扰 动 引起的相位变化量.由于脉冲产生的散射光幅度沿 利分布,对于 实 际ϕＧOTDR 系 统 来 说,可 以 选 取 信 度,即 (１１) 式中:l为干涉臂 非 平 衡 长 度;lp 为 脉 冲 长 度;l０ 为 ϕＧOTDR 一个平均衰落周期,可通过统 计ϕＧOTDR 曲线衰落周期的平均值来获得[１３]. ３　四路检测系统 基于３×３迈克 耳 孙 干 涉 仪 的 四 路 检 测 相 位 解 调系统如图５ 所 示.激 光 从 中 心 波 长 为１５５０nm、 线宽 为 ２kHz 的 激 光 器 发 出,经 过 声 光 调 制 器 (AOM)调 制 成 光 脉 冲,再 经 掺 铒 光 纤 放 大 器 (EDFA)放大,通 过 环 形 器 进 入 光 纤. 后 向 散 射 光 经环形器进入分 光 比 为１∶１的 耦 合 器,分 出 的 一 路 进入由３×３耦合器与法拉第旋转镜(FRM)构成的 迈克耳孙干涉仪,另一路直接被光电探测器(PD)接 收.四路检测系统中PD１、PD２、PD３接收的信号用 于相位 解 调,PD４ 接 收 的 信 号 用 于 对 扰 动 进 行 定 位.实验中使用直径 为６cm 的 压 电 陶 瓷(PZT)作 用于待测光纤(FUT)产生扰动信号. (４) (５) (６) ö ÷ , ø ö ÷ , ø æ ç è æ ç è E１ ＝M＋Ncosϕ０, E２ ＝M＋Ncosϕ０ ＋２ ３π E３ ＝M＋Ncosϕ０ －２ ３π ３λ(E２A０ ＋E２B０), N＝λEA０EB０, ϕ０ ＝(ϕB －ϕA), M＝１ 其中: (７) (８) (９) 式中:λ为探测器响应 系 数;EA０和EB０为 A、B 两 处 散射光的初始光场幅度;ϕA 和ϕB 为 A、B 两处散射 光的相位;通过差分Ｇ交叉相乘的解调方法得到相位 ϕ０,即为 A、B 两处散射光的相位差.解调框图如图 ４所示[１２],相对于反三角函数解调方法来说,这种解 调方法的优点是即使相位跳变也可以准确解调. 图４ 解调原理图 Fig．４ Diagramofdemodulationprinciple 当光纤上有扰动作用时,假设 A 处未受到扰动 作用,B 处在扰动影响下散射光场变为 (１０) E′B＝E′B０exp(jϕ′B). 此时,M、N 和ϕ０ 均发生变化,其 中 参 数 M 的 变 化 与散射光幅度检测的变化是一致的,ϕ０ 则与扰动强 弱呈线性关系.在相位解调过程中,N 作为与相位 无关的量需要 被 消 去.由 于 脉 冲 内 的 干 涉 效 应,N 在某些位置很小,作 为 分 母 会 放 大 噪 声 对 解 调 结 果 的影响,导致解调相位失真. 当非平衡长度较长时,在 A 处的散射光未受到 扰动影响前,B 处 的 散 射 光 受 到 扰 动 影 响 且 待 相 位 完全发生变化后,这 一 段 范 围 内 两 处 散 射 光 的 相 位 差同时变化,即 可 以 解 调 得 到 连 续 稳 定 的 相 位 差. 响的区域增加,降低了系统的空间分辨率,可以通过 增加一路直接检测 对 扰 动 进 行 定 位,免 去 采 用 数 据 处理手段来保持系统的空间分辨率.但过长的非平 衡长度也没有必要,可以选取适当的非平衡长度,使 增大非平衡长度的方法使采集到的信号受到扰动影 图５ 四路检测相位解调系统 Fig．５ FourＧchanneldetectingphasedemodulationsystem ４　实验结果与分析 首先对系统进行验证,然后在不同干涉仪非平 衡长度下进行解调,比较解调结果,最后选取合适的 非平衡长度在５km 光纤末端对PZT 扰动信号的强 弱进行定量检测.实验中采用１００Hz的 正 弦 信 号 对PZT 进行驱动,脉 冲 的 宽 度 为１００ns,重 复 频 率 为１２kHz. ４．１　结果验证 将PZT 放在距光 纤 起 始 端４２m 的 位 置,PZT 的驱动电压为２Vpp(Vpp为驱动电压的峰值电压),光 ０８０６００１Ｇ３
    光　　　学　　　学　　　报 纤绕制长度为１２m.取５km 长度光纤的散射光幅 度进行统计,得到ϕＧOTDR 衰落的平均周期.计算 得平均衰落周期约为６m.根据(１１)式选取干涉仪 非平衡长度为１６m,对应散射光的光程差为３２m, 使得进行干涉的两 个 位 置 的 散 射 光 分 开 一 定 距 离. 取扰动区域上单点 的 信 号 进 行 分 析,采 集 得 到 该 点 的干涉仪 输 出 信 号 随 时 间 的 变 化 如 图６(a)所 示. 计算得 到 的 M、N 以 及 该 点 的ϕＧOTDR 幅 度 检 测 的信号E０ 如图６(b)所 示,从 图 中 可 以 看 出 M、N 的值随时间改变,其中 M 与E０ 变化一致. 图６ (a)５６m 处单点采集的三路信号E１、E２、E３;(b)计算得到的 M、N 与E０ 的对比 应关系为 lϕ＝l－lp＋lPZT, Fig．６ a AcquiredthreeＧchannelsignalsE１ E２ E３at５６m b comparisonofE０andcalculatedM N ４．２　不同非平衡长度的解调结果及ϕＧOTDR 定位 将PZT 置于光纤起始处,此处干涉衰落位置的 信号仍然较大,解 调 失 真 较 少 以 便 于 分 析 规 律.在 ２Vpp的 PZT 驱 动 电 压 下,干 涉 仪 非 平 衡 长 度 分 别 设置为２０,３０,４０,５０m,对 PZT 扰 动 进 行 解 调,并 用ϕＧOTDR 幅度 检 测 进 行 扰 动 定 位. 光 纤 始 端 相 位解调的结果如图７所示,可以清晰地看到PZT 引 起的相位变化.其中２０m 非平衡长度下解调得到 连续且稳定相位差的范 围 约 为２２m.３０,４０,５０m 非平衡长度对应 的 范 围 分 别 约 为３２,４２,５２m. 对 (１２) 式中:lϕ 为 解 调 连 续 稳 定 相 位 的 范 围;lPZT 为 PZT 作用光纤长度.从图７(b)、(c)中 可 以 看 出 仍 有 一 小段解调结 果 受 到 脉 冲 内 干 涉 的 影 响,产 生 失 真. 因此要选取合适的 干 涉 仪 非 平 衡 长 度,使 得 进 行 干 期,以便对解调结果进行定量分析. 涉的两处 散 射 光 的 距 离 略 大 于 信 号 的 平 均 衰 落 周 图７ 非平衡长度分别为(a)２０m,(b)３０m,(c)４０m,(d)５０m 下的解调结果 Fig．７ Demodulationresultsundertheunbalancedlengthof a ２０m b ３０m c ４０m d ５０m ０８０６００１Ｇ４
    光　　　学　　　学　　　报 　　图８为光 纤 始 端 处PD４ 的 散 射 光 幅 度 检 测 结 果去均值后得到的 信 号,可 以 看 到 振 动 的 起 点 和 终 点,且不受干涉 仪 非 平 衡 长 度 影 响.故 可 以 通 过 四 路检测的相位解调系统,在保留ϕＧOTDR 的定位精 度的同时对扰动的强弱进行定量检测. 图８ 幅度检测去均值后的信号 Fig．８ MeanＧsubtractedsignalofamplitudedetection ４．３　５km 传感光纤末端定量测量结果 将PZT 置 于５km 传 感 光 纤 末 端,采 用１６m 非平衡长度,对不同驱动电压下的PZT 信号进行了 解调,结果如图９ 所 示.可 以 看 到 解 调 相 位 的 幅 度 与PZT 驱动电压呈线性关系,线性度R２＝０．９９５６, 可测量最大相位幅度达到３１．８５rad. 图９ 解调相位与PZT 驱动电压的关系 Fig．９ Relationshipbetweendemodulatedphase anddrivenvoltageofPZT ５　结　　论 为了保持基于３×３ 迈 克 耳 孙 干 涉 仪 的 相 位 解 调ϕＧOTDR 的 空 间 分 辨 率,提 出 了 四 路 检 测 的 系 统,同时考虑到高相 干 窄 线 宽 光 源 给 相 位 解 调 带 来 的影响.从散 射 光 的 角 度 对 四 路 检 测 相 位 解 调ϕＧ OTDR 的原理进行分析,并结合实验进 行 验 证.通 过分析干涉仪非平 衡 长 度 对 解 调 性 能 的 影 响,指 出 了选取适当非平衡 长 度 的 重 要 性,并 在 合 适 的 非 平 衡长度下实现 了５km 末 端 扰 动 的 定 量 检 测,检 测 线性度为０．９９５６,解 调 相 位 幅 度 达 到３１．８５rad,同 时解调系统保持了传统ϕＧOTDR 的空间分辨率. 参 考 文 献 １ 　YuK ZhangYS GuanL etal敭Structuralstate detectionwithtwistedtypefiberopticsensingarray andOTDRtechnology J 敭ActaOpticaSinica １９９９ １９ １１ １５３０Ｇ１５３５敭 　　　郁可 张永胜 管立 等敭基于光 时 域 反 射 技 术 与 缠 绕式光纤传感阵列的结构 状 态 监 测 方 法 J 敭光 学 学 报 １９９９ １９ １１ １５３０Ｇ１５３５敭 ２ 　JuarezJC TaylorH F敭Fieldtestofadistributed fiberＧoptic intrusion sensor system for long perimeters J 敭AppliedOptics ２００７ ４６ １１ １９６８Ｇ １９７１敭 ３ 　PengF Wu H JiaX H etal敭UltraＧlonghighＧ sensitivity ϕＧOTDR for high spatial resolution intrusiondetectionofpipelines J 敭OpticsExpress ２０１４ ２２ １１ １３８０４Ｇ１３８１０敭 ４ 　XieKL RaoYJ RanZL敭Distributedopticalfiber sensingsystembasedofRayleighscatteringlightϕＧ OTDRusingsingleＧmodefiberlaserwithhighpower andnarrowlinewidth J 敭ActaOpticaSinica ２００８ ２８ ３ ５６９Ｇ５７２敭 　　　谢孔利 饶云江 冉曾令敭基于大功率超窄线 宽 单 模 光纤激光器的ϕＧ光时域 反 射 计 光 纤 分 布 式 传 感 系 统 J 敭光学学报 ２００８ ２８ ３ ５６９Ｇ５７２敭 ５ 　LiuS Han X Y Xiong Y C etal敭Distributed vibrationdetectionsystembasedonweakfiberBragg gratingarray J 敭ChineseJournalofLasers ２０１７ ４４ ２ ０２１０００１敭 　　　刘胜 韩新颖 熊玉川 等敭基于 弱 光 纤 光 栅 阵 列 的 分布式振动 探 测 系 统 J 敭中 国 激 光 ２０１７ ４４ ２ ０２１０００１敭 ６ 　Shatalin S V Treschikov V N Rogers A J敭 InterferometricopticaltimeＧdomainreflectometryfor distributedopticalＧfibersensing J 敭AppliedOptics １９９８ ３７ ２４ ５６００Ｇ５６０４敭 ７ 　ParkerT Shatalin S Farhadiroushan M etal敭 Distributedacousticsensing anewtoolforseismic applications J 敭FirstBreak ２０１４ ３２ ２０１０ ６１Ｇ ６９敭 ８ 　LuYL ZhuT ChenL etal敭Distributedvibration sensorbasedoncoherentdetectionofphaseＧOTDR J 敭JournalofLightwave Technology ２０１０ ２８ ２２ ３２４３Ｇ３２４９敭 ９ 　Fang G S Xu T W FengS W etal敭PhaseＧ sensitiveopticaltimedomainreflectometerbasedon phaseＧgenerated carrier algorithm J 敭Journal of LightwaveTechnology ２０１５ ３３ １３ ２８１１Ｇ２８１６敭 １０ 　CaoY L YangF Xu D etal敭Disturbedfiber vibration sensor based on ３ × ３ Michelson ０８０６００１Ｇ５
    光　　　学　　　学　　　报 interferometer J 敭ChineseJournalofLasers ２０１６ ４３ １０ １０１０００１敭 　　　曹玉龙 杨飞 徐丹 等敭基于３×３迈克耳孙干涉仪 的分布式光纤 振 动 传 感 器 J 敭中 国 激 光 ２０１６ ４３ １０ １０１０００１敭 １１ 　WangX Liu M H Yu M etal敭Analysisand improvementfor digital quadrature demodulation algorithm on phaseＧsensitive optical timeＧdomain reflectometricsystem J 敭ChineseJournalofLasers ２０１７ ４４ １２ １２１０００２敭 　　　王旭 刘珉含 于淼 等敭相位敏 感 光 时 域 反 射 系 统 数字正交解调 算 法 分 析 及 改 进 研 究 J 敭中 国 激 光 ２０１７ ４４ １２ １２１０００２敭 １２ 　Jiang Y Lou Y M Wang H W敭 Software demodulationfor３×３ couplerbasedfiberoptic interferometer J 敭ActaPhotonicaSinica １９９８ ２７ ２ １５２Ｇ１５５敭 　　　江毅 娄英 明 王 惠 文敭基 于 对 称３×３ 耦 合 器 的 光 纤干涉 信 号 的 软 件 解 调 技 术 J 敭光 子 学 报 １９９８ ２７ ２ １５２Ｇ１５５敭 １３ 　HealeyP敭FadinginheterodyneOTDR J 敭Electronics Letters １９８４ ２０ １ ３０Ｇ３２敭 ０８０６００１Ｇ６
    分享到：
    赞 收藏

    相关推荐

    资料库

    开发技术

    共收录19306份资料，累计14个分类，关注成员有19位，主要包括：Actionscript,Javascript,PHP,VB,Perl,Delphi,Python,C#,C,Java,C++,Web开发,硬件开发,其它

    热门标签

    Actionscript Javascript PHP VB Perl Delphi Python C# C Java C++ Web开发 硬件开发 其它

    最新资料