水声学原理第一章

上传人:壹****1 文档编号:584995817 上传时间:2024-09-01 格式:PPT 页数:47 大小:14.41MB
返回 下载 相关 举报
水声学原理第一章_第1页
第1页 / 共47页
水声学原理第一章_第2页
第2页 / 共47页
水声学原理第一章_第3页
第3页 / 共47页
水声学原理第一章_第4页
第4页 / 共47页
水声学原理第一章_第5页
第5页 / 共47页
点击查看更多>>
资源描述

《水声学原理第一章》由会员分享,可在线阅读,更多相关《水声学原理第一章(47页珍藏版)》请在金锄头文库上搜索。

1、1水声学原理水声学原理上海交通大学振动、冲击、噪声国家重点实验室上海交通大学振动、冲击、噪声国家重点实验室范范 军军2024/9/1第一章第一章 导论导论1.1 海洋与水声技术海洋与水声技术1、为什么用水声技术、为什么用水声技术v 海洋占据地球表面约海洋占据地球表面约70的面积;的面积; 海洋是人类开展交通运输、军事斗争和获取资源的场所。这就必须有观测、通讯、海洋是人类开展交通运输、军事斗争和获取资源的场所。这就必须有观测、通讯、导航、定位的工具。水声技术在其中扮演了重要的角色。导航、定位的工具。水声技术在其中扮演了重要的角色。v声波是迄今为止在水中唯一能有效地远距离传递信息地物理场。声波是迄

2、今为止在水中唯一能有效地远距离传递信息地物理场。 电磁波在水中的衰减:电磁波在水中的衰减: 不能在水中远距离传播不能在水中远距离传播 声波由于介质吸收引起的衰减:声波由于介质吸收引起的衰减: 能远距离传播能远距离传播 声波与电磁波衰减之比:声波与电磁波衰减之比: 10kHz声波水中衰减仅约声波水中衰减仅约 1分贝分贝/公里公里 电磁波为电磁波为4500分贝分贝/公里公里 其它物理场:磁场、水压场、尾流场、温度场,也是可以检测,但可检测距离大其它物理场:磁场、水压场、尾流场、温度场,也是可以检测,但可检测距离大 致与源本身尺度同一量级,不能在水中远距离传递信息。致与源本身尺度同一量级,不能在水中

3、远距离传递信息。21.21.2、声呐与雷达的异同、声呐与雷达的异同声呐与雷达的工作原理相似。但由于信息载体声波与电磁波的差异决定了声呐与雷达的工作原理相似。但由于信息载体声波与电磁波的差异决定了声呐和雷达有重要差别。声呐和雷达有重要差别。a.a.电磁波速度电磁波速度3030万公里万公里/ /秒秒, ,声波在水中声波在水中1.51.5公里公里/ /秒。决定:秒。决定:工作频率差别大。雷达频率约工作频率差别大。雷达频率约GHzGHz( Hz Hz) 声呐频率约声呐频率约kHzkHz( Hz Hz)工作速率差别大。雷达搜速快,声呐搜索慢工作速率差别大。雷达搜速快,声呐搜索慢分辨率差。声图象模糊。分辨

4、率差。声图象模糊。b.b.声呐受海洋信道影响大。声呐环境比雷达环境复杂得多。声呐受海洋信道影响大。声呐环境比雷达环境复杂得多。c.c.声呐的作用距离近。声呐的作用距离近。 31.31.3、水声技术的研究范围、水声技术的研究范围 水声技术呐是研究声波在水中的发射、传输、接收、处理的专门技术。包括:水声技术呐是研究声波在水中的发射、传输、接收、处理的专门技术。包括:a.a.水声换能器和基阵水声换能器和基阵水声传感器系统;水声传感器系统;b.b.水声物理水声物理海洋信道的传播、混响、散射、噪声特性和各种水声目标特性;海洋信道的传播、混响、散射、噪声特性和各种水声目标特性;c.c.水声设备水声设备水声

5、信号处理、水声电子技术。水声信号处理、水声电子技术。水声技术的成果突出反映在两个方面水声技术的成果突出反映在两个方面1 1、声呐性能的不断提高:探测距离原来越远、对目标的定位、跟踪能力越来越强、声呐性能的不断提高:探测距离原来越远、对目标的定位、跟踪能力越来越强2 2、应用声自导或声引信的水中兵器(鱼雷、水雷、深水炸弹等)的作战能力不断、应用声自导或声引信的水中兵器(鱼雷、水雷、深水炸弹等)的作战能力不断提高。提高。 因此,现代舰艇在水下面临的威胁与水声技术的水平有直接的关系。声隐身因此,现代舰艇在水下面临的威胁与水声技术的水平有直接的关系。声隐身性能是潜艇水下隐蔽性的核心。性能是潜艇水下隐蔽

6、性的核心。41.4、水声技术的发展历史、水声技术的发展历史14901490年年 达芬奇就提出声纳的原始概念达芬奇就提出声纳的原始概念泰坦尼克号的沉没,开始最初的声纳设计泰坦尼克号的沉没,开始最初的声纳设计第一次世界大战的爆发促进了一系列军用声纳的发展(值得一提的第一次世界大战的爆发促进了一系列军用声纳的发展(值得一提的是郎之万在换能器上的贡献,并测得了水中是郎之万在换能器上的贡献,并测得了水中15001500米外潜艇回波)米外潜艇回波)一战和二战之间水声工程一直缓慢而稳步发展一战和二战之间水声工程一直缓慢而稳步发展, ,最大的成就是对海最大的成就是对海洋声传播机理的认识。(如洋声传播机理的认识

7、。(如“下午效应下午效应”现象的解释)现象的解释)二战期间为了探测德国潜艇,水声工程有了很大发展,出现了大量二战期间为了探测德国潜艇,水声工程有了很大发展,出现了大量新的理论和技术新的理论和技术战后水声工程随着计算机和电子计算发展,水声工程的应用在军用、战后水声工程随着计算机和电子计算发展,水声工程的应用在军用、民用领域更为广泛。民用领域更为广泛。5 声呐(声纳)声呐(声纳)-SONAR(Sound Navigation and Ranging) 凡是利用水下声信息进行探测、识别、定位、导航和通讯的系统,都通称为声呐系凡是利用水下声信息进行探测、识别、定位、导航和通讯的系统,都通称为声呐系统。

8、声呐的主要应用是军用声呐。按工作方式可以分为:统。声呐的主要应用是军用声呐。按工作方式可以分为:主动声呐主动声呐和和被动声呐。被动声呐。按安装平台分可以分为:按安装平台分可以分为:v潜艇声呐:潜艇上的电子设备是声呐。一般核潜艇装有潜艇声呐:潜艇上的电子设备是声呐。一般核潜艇装有1015部声呐。主要有:艏部声呐。主要有:艏部主、被动综合声呐;被动测距声呐;舷侧阵声呐;拖曳线列阵声呐。部主、被动综合声呐;被动测距声呐;舷侧阵声呐;拖曳线列阵声呐。v水面舰声呐:舰艏声呐;变深拖曳声呐;拖曳线列阵声呐。水面舰声呐:舰艏声呐;变深拖曳声呐;拖曳线列阵声呐。v机载声呐和浮标:吊放声呐;声呐浮标。机载声呐和

9、浮标:吊放声呐;声呐浮标。v海洋水声监视系统:岸站(岸边海底固定式声呐);预警系统海洋水声监视系统:岸站(岸边海底固定式声呐);预警系统v水声对抗器材:鱼雷报警声呐;声诱饵;干扰器;气幕弹水声对抗器材:鱼雷报警声呐;声诱饵;干扰器;气幕弹v水中兵器自导:鱼雷声自导;水雷声引信;水中兵器自导:鱼雷声自导;水雷声引信;v其它:通讯仪、鱼探仪、多普勒测速仪、浅地层剖面仪等。其它:通讯仪、鱼探仪、多普勒测速仪、浅地层剖面仪等。1.5 声呐简介声呐简介6主要声呐图片主要声呐图片7德国德国ATLAS公司研制的公司研制的拖曳线列阵拖曳线列阵英国、法国联合研制的英国、法国联合研制的舷侧阵声呐舷侧阵声呐TSM2

10、253德国德国ATLAS公司研制的公司研制的ASA9225主动拖曳线声呐主动拖曳线声呐英国、法国联合研制的英国、法国联合研制的投吊声呐投吊声呐美国美国Lockheed Martin公司研制的被动测距声公司研制的被动测距声呐呐PUFFS美国美国DTI公司研制的合公司研制的合成孔径声呐成孔径声呐 81.6 声学量的度量、分贝和级声学量的度量、分贝和级声学中采用分贝计量的原因:声学中采用分贝计量的原因:v声学量的变化大到六、七个数量级以上声学量的变化大到六、七个数量级以上 从窃窃私语到大型喷气式飞机起飞的声功率差十个数量级;从窃窃私语到大型喷气式飞机起飞的声功率差十个数量级; 人耳的听阈在频率人耳的

11、听阈在频率1 1kHzkHz时是时是2020PaPa(微帕),痛阈是微帕),痛阈是2020PaPa, 相差六个数量级;相差六个数量级; 在水中,一艘老式潜艇的辐射总声功率达到数瓦,而新型的在水中,一艘老式潜艇的辐射总声功率达到数瓦,而新型的 低噪声潜艇不到低噪声潜艇不到1 1微瓦,相差六、七个数量级。微瓦,相差六、七个数量级。 v人耳(仪器)的响应近似与声压或声强的对数成比例。人耳(仪器)的响应近似与声压或声强的对数成比例。 因此声学中定义一个以对数为基础的因此声学中定义一个以对数为基础的分贝分贝单位,水声也一直沿用。单位,水声也一直沿用。91.6.11.6.1、定义和参考、定义和参考声压、声

12、强和声功率声压、声强和声功率用级和分贝(用级和分贝(dBdB)来量度。他们是:来量度。他们是:参考值 101.6.21.6.2声压级等于声强级声压级等于声强级: :注意参考值不同产生的声级差别注意参考值不同产生的声级差别: :19711971年以前曾用:年以前曾用: =20 =20PaPa2102104 4达因达因/ /厘米厘米2 2,换算到现在标准要加,换算到现在标准要加2626分贝。分贝。 1 1 达因达因/ /厘米厘米2 21 1bb(微巴)微巴)10105 5PaPa,换算到现在的标准换算到现在的标准 要加要加100100分贝。分贝。 俄罗斯标准俄罗斯标准2020PaPa 由于空气声和

13、水声参考值的不同,舱室内声级为由于空气声和水声参考值的不同,舱室内声级为L L分贝的噪声若分贝的噪声若 无损耗地传到水下将变成无损耗地传到水下将变成L L2626分贝的水噪声。分贝的水噪声。11 固体介质中的结构噪声用振动来描述,它的分贝定义实际上就是振固体介质中的结构噪声用振动来描述,它的分贝定义实际上就是振动量的分贝定义。动量的分贝定义。加速度级加速度级 速度级速度级 位移级位移级 加速度、速度和位移参考值是:加速度、速度和位移参考值是: 米/秒2, 米/秒, 米。 应当指出的是,虽然结构噪声级与振动级的定义相同,应当指出的是,虽然结构噪声级与振动级的定义相同,但实际测量和评价方法有区别的

14、。因为结构噪声要反映但实际测量和评价方法有区别的。因为结构噪声要反映连续弹性体的振动特性,所以用一个点的振动级是无法连续弹性体的振动特性,所以用一个点的振动级是无法描述的。通常要用结构的整个辐射面上大量测点的统计描述的。通常要用结构的整个辐射面上大量测点的统计平均来描述。平均来描述。121.6.3 1.6.3 分贝表示的特点分贝表示的特点物理量的乘除运算变成加减运算。物理量的乘除运算变成加减运算。 例如在声学测量中,用灵敏度等于例如在声学测量中,用灵敏度等于S S 伏伏/ /PaPa的水听器接收,经的水听器接收,经过放大倍数等于过放大倍数等于K K的放大器放大后得到电压的放大器放大后得到电压V

15、 V伏。水听器输入端的声压伏。水听器输入端的声压是是 : : (Pa) 声压级声压级: :如果水听器灵敏度、声学测量放大器的放大倍数都用分贝表示,只要简单的如果水听器灵敏度、声学测量放大器的放大倍数都用分贝表示,只要简单的加减运算就可以求出声压级。加减运算就可以求出声压级。13声学中不仅声学量用分贝表示,它们的误差范围也用误差表示,声学中不仅声学量用分贝表示,它们的误差范围也用误差表示,例如例如 。用分贝表示的误差与百分比误差的换算关系:设声压是用分贝表示的误差与百分比误差的换算关系:设声压是分贝表示是分贝表示是 ,则有,图给出其关系曲线:,则有,图给出其关系曲线:14用分贝表示后函数图形发生

16、变化用分贝表示后函数图形发生变化声学中最常见的幂次规律声学中最常见的幂次规律 : :以以 为横坐标声级是一条直线。从直线的斜率可以确定幂次为横坐标声级是一条直线。从直线的斜率可以确定幂次n n。最方最方便的方法是根据频率加倍时声级减小的分贝数便的方法是根据频率加倍时声级减小的分贝数 得到得到 声强随声强随频率衰频率衰减规律减规律声强随声强随距离衰距离衰减规律减规律151.6.4 1.6.4 分贝的基本运算分贝的基本运算相干叠加相干叠加 在讨论分贝运算法则前先要搞清楚声场的叠加原则。因为声压场是标量在讨论分贝运算法则前先要搞清楚声场的叠加原则。因为声压场是标量场,具有可加性。但是,它又是一个波动

17、场,既有振幅又有相位。场,具有可加性。但是,它又是一个波动场,既有振幅又有相位。 相干叠加:相干叠加:当两个以上的有规声波叠加时要同时计及振幅和相位,若是同当两个以上的有规声波叠加时要同时计及振幅和相位,若是同频率的声波叠加会发生干涉现象。若是频率相差不多的两个声波叠加会发频率的声波叠加会发生干涉现象。若是频率相差不多的两个声波叠加会发生生“拍拍”。这些情况称为相干叠加:。这些情况称为相干叠加: 水下声基阵形成指向性就是靠声波到达不同阵元的相位差。水下声基阵形成指向性就是靠声波到达不同阵元的相位差。能能量量叠叠加加:当当两两个个以以上上的的随随机机噪噪声声叠叠加加时时,由由于于相相位位是是随随

18、机机的的,噪噪声声要要按按强强度或能量叠加。因为相位的随机性导致所有的交叉项度或能量叠加。因为相位的随机性导致所有的交叉项 ,所以,所以 。相相当当于于按按强强度度或或能能量量叠叠加加。有有时时候候几几个个噪噪声声之之间间有有一一定定的的相相干干性性,但但是是相相干干性性无无法法估估计计或或测测量量,也也采采用用能能量量叠叠加加的的原原则则处处理理。其其结结果果是是对对相相干干叠叠加加的的一一种种平平均均。所所以以,能能量量叠叠加加是是噪噪声声场场叠叠加加的的基基本本原原则则。因因此此也是分贝计算的指导原则。也是分贝计算的指导原则。161.6.4.1 1.6.4.1 噪声叠加噪声叠加特例:N个

19、声级相同的噪声叠加,总声级是单个噪声声级加 分贝。当N2即两个声级相同的噪声叠加,总声级增加 分贝。当声级为当声级为 的的N个噪声叠加时,按强度叠加得到总声强和个噪声叠加时,按强度叠加得到总声强和总声级总声级:171.6.4.2 1.6.4.2 噪声相减(背景噪声的扣除)噪声相减(背景噪声的扣除) 在实际工作中有时需要从总的噪声中减去某个噪声成分,例如在噪声测量中扣在实际工作中有时需要从总的噪声中减去某个噪声成分,例如在噪声测量中扣除已知的背景干扰。假设总噪声级是除已知的背景干扰。假设总噪声级是 ,背景干扰级是,背景干扰级是 ,扣除后的噪声级,扣除后的噪声级是是 : 图给出修正量图给出修正量

20、与声级差与声级差 的关系的关系曲线。为了保证合理的精度,背景干扰至少曲线。为了保证合理的精度,背景干扰至少要比总声级低要比总声级低3 3dBdB。若背景干扰比总声级低若背景干扰比总声级低1010dBdB以上就可以不修正。以上就可以不修正。 181.6.4.31.6.4.3多个噪声级的平均多个噪声级的平均 对噪声级进行多次测量需要计算其平均值。设对噪声级进行多次测量需要计算其平均值。设N N 次测量的噪声级分别是次测量的噪声级分别是 ,应根据声强的算术平均值计算平均噪声级,应根据声强的算术平均值计算平均噪声级: :若各次测量的噪声级的差值小于若各次测量的噪声级的差值小于3 3dBdB,则可以直接

21、取噪声级的算术平均值代则可以直接取噪声级的算术平均值代替替上上式,得到:式,得到:其其误误差差不不超超过过0.50.5dBdB;若若各各次次测测量量的的噪噪声声级级的的差差值值小小于于5 5dBdB,误误差差不不超超过过0.70.7dBdB。191.6.4.5 1.6.4.5 降噪量的计算降噪量的计算 已已知知总总噪噪声声级级为为140140dBdB,它它由由三三个个声声级级相相同同的的噪噪声声叠叠加加而而成成。可可以以求求出出每每个噪声源的声级是个噪声源的声级是135.2135.2dBdB。三种降噪方案的效果是:三种降噪方案的效果是:方案一、将其中一个噪声源降低方案一、将其中一个噪声源降低1

22、010dBdB,另外两个不变,总噪声级另外两个不变,总噪声级只下降只下降1.61.6dBdB;方案二、将其中的两个噪声源降低方案二、将其中的两个噪声源降低1010dBdB,总噪声级下降了总噪声级下降了4 4dBdB;方案三、将三个噪声源都下降方案三、将三个噪声源都下降1010dBdB,总噪声级也下降总噪声级也下降1010dBdB。 例:已知总噪声级由例:已知总噪声级由N N个声级分别是个声级分别是 的噪声叠加而成的噪声叠加而成。当这。当这N N个噪声分别降低个噪声分别降低 分贝后,总声级降低为分贝后,总声级降低为: :总降噪量是总降噪量是 : :201.7 1.7 频谱和频谱级频谱和频谱级水声

23、信号、特别是噪声信号常常包含有多种频率成分,能量分布在一个频带宽水声信号、特别是噪声信号常常包含有多种频率成分,能量分布在一个频带宽度内。定义单位频带宽度度内。定义单位频带宽度1 1HzHz内的声强度为声强谱密度,用函数内的声强度为声强谱密度,用函数 表示。谱表示。谱密度的分贝表示称为谱密度级。密度的分贝表示称为谱密度级。 将谱密度函数在整个频带内积分就等于总强度:将谱密度函数在整个频带内积分就等于总强度:在频率在频率f f附近带宽内的声强是附近带宽内的声强是:用分贝表示:用分贝表示: 称为频带级,频带级等于谱密度级加称为频带级,频带级等于谱密度级加 。上面的分贝计算法则同样适。上面的分贝计算

24、法则同样适用于频带级和谱密度级的计算。用于频带级和谱密度级的计算。潜艇辐射噪声潜艇辐射噪声21用声压表示时:用声压表示时: 参考值应该理解为:参考值应该理解为: 在声学测量中用到两种滤波器:在声学测量中用到两种滤波器: 恒定带宽滤波器:低频时分析太粗,高频分析太细,无法兼顾。恒定带宽滤波器:低频时分析太粗,高频分析太细,无法兼顾。 恒定百分比或恒定百分比或Q滤波器:用的多,人耳听觉模型是其的组合。声学中的恒定滤波器:用的多,人耳听觉模型是其的组合。声学中的恒定 百分比滤波器称为倍频程滤波器。百分比滤波器称为倍频程滤波器。22设设 是滤波器的下限频率是滤波器的下限频率设设 是滤波器的下限频率是滤

25、波器的下限频率n n倍频程滤波器的数学定义为:倍频程滤波器的数学定义为:或或中心频率为:中心频率为: 带宽:带宽:n=1 1/1倍频程(倍频程(1/1oct) n=1/3 1/3倍频程(倍频程(1/3oct) 对于对于1/31/3octoct来讲,是在间隔为来讲,是在间隔为1 1倍频程的两个频率之间再插入两个频率,使这倍频程的两个频率之间再插入两个频率,使这4 4个频率之间依次相距个频率之间依次相距1/31/3倍频程,即这倍频程,即这4 4个频率值按比例为:个频率值按比例为: ,近似为,近似为 。ISOISO规定规定1 11010HzHz之之间划分间划分1010个频段,中心频率为:个频段,中心

26、频率为:1 1:1.25:1.60:2.0:2.5:3.15:4.0:5.0:6.3:8.0:10.0:1.25:1.60:2.0:2.5:3.15:4.0:5.0:6.3:8.0:10.0。这样的取值。这样的取值得好处是,可使每隔得好处是,可使每隔1010个个1/31/3倍频程频段的两端频率正好相差倍频程频段的两端频率正好相差1010倍,还可以使每隔倍,还可以使每隔3 3个个1/31/3倍频程频倍频程频段为一个段为一个1/11/1倍频程频段。倍频程频段。 231.8 1.8 水声(海军)的一些习惯用量水声(海军)的一些习惯用量 距离单位:英尺(距离单位:英尺(feet),),1英尺英尺12英

27、寸英寸0.3048米米 千码(千码(Kiloyard),1千码千码3000英尺英尺914.4米米 海里(海里(Nautical mile),1海里海里1853米米 链链 ,1链链=1/10海里海里 深度单位:英尺深度单位:英尺 浔(浔(fathom),1浔浔2码码6英尺英尺1.83米米速度单位:速度单位: 节(节(kont),),1节节1海里海里/小时小时0.5米米/秒秒241.9 1.9 声呐方程声呐方程水声探测技术包括水声探测技术包括主动方式主动方式和和被动方式被动方式。相应的声呐也分为。相应的声呐也分为主动声呐主动声呐和和被动声呐被动声呐。主动声呐:主动声呐:由探测设备主动的发射声波,通

28、过接收、分析目标回波实现对目标的由探测设备主动的发射声波,通过接收、分析目标回波实现对目标的探测、定位和识别。探测、定位和识别。 特点:特点:定位和测距精度高,容易暴露自己,非隐蔽探测,今年来又逐渐重要。定位和测距精度高,容易暴露自己,非隐蔽探测,今年来又逐渐重要。被动声呐:被动声呐:由探测设备被动的接收、分析目标发出的噪声对目标进行探测、定位和由探测设备被动的接收、分析目标发出的噪声对目标进行探测、定位和识别。识别。特点:特点:定位和测距精度不如主动声呐高,是隐蔽探测,是目前主要探测方式,被动定位和测距精度不如主动声呐高,是隐蔽探测,是目前主要探测方式,被动声呐包括:水雷引信、鱼雷自导、舰艇

29、被动声呐、拖曳线列阵、海岸预警系统。声呐包括:水雷引信、鱼雷自导、舰艇被动声呐、拖曳线列阵、海岸预警系统。检测阈( DT )指向性指数 ( DT ) 噪声级 ( NL ) 传播损失 ( TL ) 声源级 ( SL )检测阈( DT )指向性指数 ( DI) 声源 声源级 ( DT ) 噪声级 ( NL ) 传播损失 ( TL ) 传播损失 ( TL ) 目标强度 ( TS ) 噪声级 ( NL )被动声呐被动声呐主动声呐主动声呐25 声呐是为了完成特定使命而构成的水声系统。它的性能取决于声呐是为了完成特定使命而构成的水声系统。它的性能取决于目标、传输信目标、传输信道和接收、处理设备道和接收、处

30、理设备的特性。声呐方程是将这三者联系在一起的一个关系式。的特性。声呐方程是将这三者联系在一起的一个关系式。它是声呐系统工作时必须服从的一个关系式。应用这个关系式可以对声呐系统它是声呐系统工作时必须服从的一个关系式。应用这个关系式可以对声呐系统的工作特性作出估计和预报,是声呐设计的基础。的工作特性作出估计和预报,是声呐设计的基础。 无论声呐系统如何复杂,要完成一定的使命必须保证在其输入端满足:无论声呐系统如何复杂,要完成一定的使命必须保证在其输入端满足:或:或:这里检测阈是接收、处理设备对信号作出判决的一个阈值,用这里检测阈是接收、处理设备对信号作出判决的一个阈值,用DTDT表示。表示。DTDT

31、取决取决于:于:1 1、使命的性质;、使命的性质;2 2、完成使命的质量,通常用检测概率和虚惊概率来表示。、完成使命的质量,通常用检测概率和虚惊概率来表示。3 3、接收、处理设备的能力。、接收、处理设备的能力。261.9.1 1.9.1 声呐方程中的各种参数声呐方程中的各种参数在声呐方程中出现的声呐参数可分为三类:在声呐方程中出现的声呐参数可分为三类: 由声呐系统决定的参数,包括:声源级由声呐系统决定的参数,包括:声源级SL:、:、自噪声级自噪声级NL、空间增益空间增益GS(或或DI)、)、时间时间 增益增益GT、检测阈检测阈DT;取决与被探测目标的参数,包括:辐射声源级取决与被探测目标的参数

32、,包括:辐射声源级SL、目标强度目标强度TS;取决于环境的参数,包括:传播损失取决于环境的参数,包括:传播损失TL、海洋环境噪声级海洋环境噪声级NL、等效平面波混等效平面波混 响级响级RL其中两对参数用了相同的符号其中两对参数用了相同的符号(SL与与NL),因为它们本质上是相同的。这些参因为它们本质上是相同的。这些参数都用分贝(数都用分贝(dB)表示,通常还与频率有关。下面我们分别介绍这些参数的定表示,通常还与频率有关。下面我们分别介绍这些参数的定义。义。271.9.2 1.9.2 声源级声源级SLSL 声源级声源级SLSL用来描述主动声呐所发射的声信号强弱,定义为:用来描述主动声呐所发射的声

33、信号强弱,定义为:I I是发射换能器声轴上离声源中心是发射换能器声轴上离声源中心1 1米处的声强,米处的声强,I I0 0是参考声强,约为是参考声强,约为: : 。 为为了了有有效效地地提提高高主主动动声声纳纳的的作作用用距距离离,它它的的发发射射器器总总是是做做成成具具有有一一定定的的发发射射指指向向性性,使使它它所所发发射射的的声声能能主主要要集集中中于于空空间间某某一一方方向向( (通通常常就就是是目目标标所所在在的的方方向向) ),其其余余方方向向上上则则仅仅有有很很少少量量的的发发射射声声能能,下下图图形形象象地地表表示示了了这这种种发发射射指指向性特性。描述发射换能器的这种特性的参

34、数为向性特性。描述发射换能器的这种特性的参数为DIDIT T。281.9.3 1.9.3 发射指向性指数发射指向性指数DIDIT T定义为:设有两个发射相同声功率的发射器,一个具有发射指向性,另一个无发定义为:设有两个发射相同声功率的发射器,一个具有发射指向性,另一个无发射指向性,又设在它们各自的辐射声场的远场测量声强,测量距离相同,测得无射指向性,又设在它们各自的辐射声场的远场测量声强,测量距离相同,测得无指向性发射器辐射声强度为指向性发射器辐射声强度为 I INDND 。在指向性发射器声轴上测得的声强度为在指向性发射器声轴上测得的声强度为I ID D 。见上图见上图则指向性发射器的发射指向

35、性指数则指向性发射器的发射指向性指数DTDTI I 可见,发射指向性指数可见,发射指向性指数DIDIT T实际上就是在相同的距离上,指向性发射器声轴上实际上就是在相同的距离上,指向性发射器声轴上的声级高出无指向性发射器辐射声场声级的分贝数。的声级高出无指向性发射器辐射声场声级的分贝数。DIDIT T值愈大,就表示了声能值愈大,就表示了声能在声轴方向集中的程度愈高,也就是获得了一定的空间增益,也可以表述为在声轴方向集中的程度愈高,也就是获得了一定的空间增益,也可以表述为GSGS,就愈有利于增加设备的作用距离。所以,近代主动声呐的发射器,都在造价、工就愈有利于增加设备的作用距离。所以,近代主动声呐

36、的发射器,都在造价、工程实施等允许的条件下,尽可能地提高发射指向性指数。程实施等允许的条件下,尽可能地提高发射指向性指数。29 发射器的声源级反应了发射器辐射声功率的大小,它们之间有着简单的函发射器的声源级反应了发射器辐射声功率的大小,它们之间有着简单的函数关系。设在无吸收的介质中有一个辐射声功率为数关系。设在无吸收的介质中有一个辐射声功率为P Pa a(w)(w)的点声源,根据声学的点声源,根据声学基础知识可知,距此声源声中心单位距离处的声强度为:基础知识可知,距此声源声中心单位距离处的声强度为:注意到:注意到: ,可以得到发射声源级,可以得到发射声源级SL:对于一个发射声功率为对于一个发射

37、声功率为P Pa a(w) (w) 、指向性指数为指向性指数为DIDIT T的指向性发射器,声源级的指向性发射器,声源级SLSL: 目前,船用声纳的辐射声功率范围为几百瓦到几十干瓦,发射指向性指数目前,船用声纳的辐射声功率范围为几百瓦到几十干瓦,发射指向性指数为为10103030dB,dB,所以其声源级范围约为所以其声源级范围约为210210dBdB到到240240dBdB。 301.9.4 1.9.4 传播损失传播损失TLTL 海水介质是一种不均匀的非理想介质,由于介质本身的吸收、声传播过程中波阵面海水介质是一种不均匀的非理想介质,由于介质本身的吸收、声传播过程中波阵面的扩展及海水中各种不均

38、匀性的散射等原因,声波在传播过程中,声传播方向上的的扩展及海水中各种不均匀性的散射等原因,声波在传播过程中,声传播方向上的声强度将会逐渐减弱,传播损失声强度将会逐渐减弱,传播损失TLTL定量地描述了声波传播一定距离后声强度的衰减定量地描述了声波传播一定距离后声强度的衰减变化,它定义为:变化,它定义为: 式中,式中,I I1 1是离声源声中心是离声源声中心1 1m m处的声强度;处的声强度;I Ir r是距离声源是距离声源r r处的声强度。上式处的声强度。上式定义的传播损失定义的传播损失TLTL值总为正值。值总为正值。311.9.5 1.9.5 目标强度目标强度TSTS 对于主动声呐而言,它是利

39、用目标回波来实现检测的。由声学基础知识可知,对于主动声呐而言,它是利用目标回波来实现检测的。由声学基础知识可知,目标回波的特性除和声波本身的特性如频率、波阵面形状等因素有关外,还与目目标回波的特性除和声波本身的特性如频率、波阵面形状等因素有关外,还与目标的特性如几何形状、组成材料等有关,也就是说,即使是在同样的入射波照射标的特性如几何形状、组成材料等有关,也就是说,即使是在同样的入射波照射下,不同目标的回波也将是不一样的。这一现象反应了目标反射本领的差异。水下,不同目标的回波也将是不一样的。这一现象反应了目标反射本领的差异。水声技术中,通常用目标强度声技术中,通常用目标强度TSTS定量描述目标

40、反射本领的大小,它定义为定量描述目标反射本领的大小,它定义为 式中,式中,I Ii i入是目标处入射声波的强度;入是目标处入射声波的强度;I Ir r=1=1 是在入射声波相反的方向上、离目是在入射声波相反的方向上、离目标声中心标声中心1 1m m处的回波强度。目标强度是目前我们这个研究小组的主要研究方向。处的回波强度。目标强度是目前我们这个研究小组的主要研究方向。321.9.6 1.9.6 海洋环境噪声级海洋环境噪声级NLNL 海水介质中,存在着大量的、各种各样的噪声源。它们各自发出的声波构成海水介质中,存在着大量的、各种各样的噪声源。它们各自发出的声波构成了海洋环境噪声。这种环境噪声,对声

41、纳设备的工作无疑是一种干扰。环境噪了海洋环境噪声。这种环境噪声,对声纳设备的工作无疑是一种干扰。环境噪声级声级NLNL就是用来度量环境噪声强弱的一个量,它定义为:就是用来度量环境噪声强弱的一个量,它定义为: 式中,式中,I IN N是测量带宽内(或是测量带宽内(或lHzlHz频带内频带内) )的噪声强度的噪声强度。 I I0 0是参考声强。是参考声强。331.9.7 1.9.7 等效平面波混响级等效平面波混响级RLRL 对于主动声纳来说,除了环境噪声是背景干扰外,很响也是一种背景干扰。关对于主动声纳来说,除了环境噪声是背景干扰外,很响也是一种背景干扰。关于海水混响的研究指出,混响不同于环境噪声

42、它不是平稳的,也不是各向同性于海水混响的研究指出,混响不同于环境噪声它不是平稳的,也不是各向同性的。为了定量描述混响干扰的强弱,我们引入声纳参数等效平面波混响级的。为了定量描述混响干扰的强弱,我们引入声纳参数等效平面波混响级NLNL。设设有强度为有强度为I I的平面波轴向入射到水听器上,水听器输出某一电压值;如将此水听的平面波轴向入射到水听器上,水听器输出某一电压值;如将此水听器移置于混响场中,使它的声轴指向目标,在混响声的作用下,水听器也输出一器移置于混响场中,使它的声轴指向目标,在混响声的作用下,水听器也输出一个电压。如果这两种情况下水听器的输出相等,那么,就用该平面波的声级来度个电压。如

43、果这两种情况下水听器的输出相等,那么,就用该平面波的声级来度量混响场的强弱,称为等效平面波混响级量混响场的强弱,称为等效平面波混响级RLRL。 具体的计算方法和公式我们将在以后的课程中详细讨论。具体的计算方法和公式我们将在以后的课程中详细讨论。输出电压VI输出电压V平面波平面波混响场混响场341.9.8 1.9.8 接收指向性指数接收指向性指数DIDI接收换能器的接收指向性指数的定义是:接收换能器的接收指向性指数的定义是:设有两个水听器一个无指向性,另一个有指向性,且指向性水听器的轴向灵设有两个水听器一个无指向性,另一个有指向性,且指向性水听器的轴向灵敏度等于无指向性水听器的灵敏度,设为单位值

44、。现将它们置于单位立体角内敏度等于无指向性水听器的灵敏度,设为单位值。现将它们置于单位立体角内的噪声功率为的噪声功率为I Ii i的各向同性噪声场中,此时无指向性水听器产生的均方电压的各向同性噪声场中,此时无指向性水听器产生的均方电压是:是:有指向性水听器产生的均方电压是:有指向性水听器产生的均方电压是:则:则:归一化声束图案函数归一化声束图案函数351.9.9 1.9.9 检测阈检测阈DTDT 声呐设备的接收器工作在噪声环境中,既接收声纳信号,也接收背景噪声,相应声呐设备的接收器工作在噪声环境中,既接收声纳信号,也接收背景噪声,相应地其输出也由这两部分组成。因此这两部分比值的大小将直接影响设

45、备的工作地其输出也由这两部分组成。因此这两部分比值的大小将直接影响设备的工作质量,即如果接收带宽内的信号功率质量,即如果接收带宽内的信号功率( (或均方电压或均方电压) )与与1 1HzHz带宽内带宽内( (或工作带宽内或工作带宽内) )的的噪声功率噪声功率( (或均方电压或均方电压) )的比值较高,则设备就能正常工作,它作的的比值较高,则设备就能正常工作,它作的“判决判决”也是可也是可信的;反之,上述的信噪比值比较低时,设备就不能正常工作,它作出的信的;反之,上述的信噪比值比较低时,设备就不能正常工作,它作出的“判决判决”也就不可信。在水声技术中,习惯上将设备刚好能正常工作所需的处理器输入端

46、也就不可信。在水声技术中,习惯上将设备刚好能正常工作所需的处理器输入端的信噪比值的信噪比值( (用分贝表示用分贝表示) )称作检测阂,它定义为:称作检测阂,它定义为: 由检测阈定义可知,对于完成同样职能的声呐来说,检测阈值较低的设备,其处由检测阈定义可知,对于完成同样职能的声呐来说,检测阈值较低的设备,其处理能力较强,其性能也较好。理能力较强,其性能也较好。361.9.10 1.9.10 时间增益时间增益GTGT声呐信号处理中要完成一系列的过程,这中间时间积累是不可缺少的,因为随声呐信号处理中要完成一系列的过程,这中间时间积累是不可缺少的,因为随机信号的起伏特征只有依靠时间积累才能去除。时间增

47、益定义为:机信号的起伏特征只有依靠时间积累才能去除。时间增益定义为:其中其中(SNR)out是系统输出信噪比,是系统输出信噪比, (SNR)in是系统输入信噪比。对于被动声呐和主是系统输入信噪比。对于被动声呐和主动声呐的最佳检测系统来说,其时间增益动声呐的最佳检测系统来说,其时间增益GT为:为:被动声呐最佳检测系统:被动声呐最佳检测系统:是噪声相关特性的一个度量,可以认为是相关半径。是噪声相关特性的一个度量,可以认为是相关半径。是噪声的自相关函数是噪声的自相关函数主动声呐最佳检测(匹配主动声呐最佳检测(匹配滤波)系统:滤波)系统:信号脉宽信号脉宽或或是输入信号带宽是输入信号带宽信号带宽信号带宽

48、对于方波调制的单频信号对于方波调制的单频信号 2WT1 于是于是GT=0dB371.9.10 1.9.10 主动声呐方程主动声呐方程检测阈( DT )指向性指数 ( DI) 声源 声源级 ( SL) 噪声级 ( NL ) 传播损失 ( TL ) 传播损失 ( TL ) 目标强度 ( TS ) 噪声级 ( NL )对于主动声呐,回声信号级是自身发射的信号传播到目标并经过目标反射回到对于主动声呐,回声信号级是自身发射的信号传播到目标并经过目标反射回到声呐接收端的值,表示为声呐接收端的值,表示为SL-2TL-TSSL-2TL-TS。主动声呐的干扰在多数情况下是混响级主动声呐的干扰在多数情况下是混响级

49、RLRL,当探测距离较远时也可能是噪声背景当探测距离较远时也可能是噪声背景NLNL,这里的噪声包括海洋环境噪声和这里的噪声包括海洋环境噪声和舰艇自噪声。这样单水听器接收的信号信噪(混)比为:舰艇自噪声。这样单水听器接收的信号信噪(混)比为:(SNR)in=SL-2TL+TS-RL混响背景下:混响背景下:噪声背景下:噪声背景下:(SNR)in=SL-2TL+TS-NL加上声呐系统的空间增益加上声呐系统的空间增益GS(DI)和时间增益和时间增益GT可以得到输出信噪(混)比:可以得到输出信噪(混)比:(SNR)out=SL-2TL+TS-NL+GS+GT得到主动声呐方程:得到主动声呐方程:SL-2T

50、L+TS-NL+GS(DI)+GTDT混响背景下:混响背景下:SL-2TL+TS-RL+GS(DI)+GTDT噪声背景下:噪声背景下:38值得注意的是:值得注意的是:1、上述的主动声呐方程适用于收发合置的声呐,对于收发分置的多基地声、上述的主动声呐方程适用于收发合置的声呐,对于收发分置的多基地声 呐,信号的传播损失一般是不同的,不能用简单的呐,信号的传播损失一般是不同的,不能用简单的2TL表示往返的传播损失。表示往返的传播损失。2、主动声呐要经受双层传播损失,这时主动声呐探测距离受到限制的重要原、主动声呐要经受双层传播损失,这时主动声呐探测距离受到限制的重要原因。因。3、以上给出的主动声呐方程

51、与传统的声呐方程形式形式上有差别,主要是在、以上给出的主动声呐方程与传统的声呐方程形式形式上有差别,主要是在声呐发展中基阵处理以及时间处理得到发展,我们考虑了基阵的空间增益声呐发展中基阵处理以及时间处理得到发展,我们考虑了基阵的空间增益GS和时间增益和时间增益GT。一般来说,利用声呐方程的目的是要从中求出声呐的作用距离,因此我们定义:一般来说,利用声呐方程的目的是要从中求出声呐的作用距离,因此我们定义:优质因素FOM(Figure Of Merit)FOM不仅与声呐系统有关,还与海洋环境有关。令不仅与声呐系统有关,还与海洋环境有关。令便可求出便可求出r。391.9.11 1.9.11 被动声呐

52、方程被动声呐方程对于被动声呐信号级是目标辐射噪声经过信道传输后到达声呐基阵的值,可以对于被动声呐信号级是目标辐射噪声经过信道传输后到达声呐基阵的值,可以表示为表示为SL-TL。干扰级是声呐平台自噪声和海洋环境噪声级干扰级是声呐平台自噪声和海洋环境噪声级NL。因此单水听器因此单水听器接收到的信号为:接收到的信号为:SL-TL-NL,信噪比为:信噪比为:(SNR)in=SL-TL-NL经过声呐基阵信号处理,获得空间增益经过声呐基阵信号处理,获得空间增益GS(DI)和时间增益和时间增益GT,得到输出信噪比:得到输出信噪比:(SNR)out=SL-TL-NLGSGT被动声呐方程:被动声呐方程:SL-T

53、L-NLGSGTDT检测阈( DT )指向性指数 ( DI ) 噪声级 ( NL ) 传播损失 ( TL ) 声源级 ( SL )被动声呐方程优质因素:被动声呐方程优质因素:FOMSL-NLGSGTDT401.9.12 1.9.12 声呐方程的基本应用声呐方程的基本应用 例一:例一:设一静止状态工作方式的声呐站,发射声源级为设一静止状态工作方式的声呐站,发射声源级为115分贝,对一个目标强分贝,对一个目标强度为度为15分贝,并位于足够远的目标进行探测,在这一段距离上的总单程传播衰减分贝,并位于足够远的目标进行探测,在这一段距离上的总单程传播衰减为为81分贝,声呐所在的海洋环境噪声和本身总噪声级

54、为分贝,声呐所在的海洋环境噪声和本身总噪声级为40分贝,接收具有指向分贝,接收具有指向性,指向性指数为性,指向性指数为12分贝,求该设备需要多大的检测阈才能可靠探测目标?分贝,求该设备需要多大的检测阈才能可靠探测目标?解:已知SL=115dB, TL=81dB ,NL=-40dB,GS(DI)=12dB,TS=15dB由于在噪声背景下探测目标:使用下列主动声呐方程:由于在噪声背景下探测目标:使用下列主动声呐方程:SL-2TL+TS-NL+DIDT115-2*81+15-(-40)+12=DTDT=20dB结果表明该设备的输入信噪比在结果表明该设备的输入信噪比在20dB时才可能可靠检测到目标。时

55、才可能可靠检测到目标。41例二:例二:设一被动探潜声呐,目标潜艇以设一被动探潜声呐,目标潜艇以60分贝的发射声源级发射分贝的发射声源级发射500Hz的线谱,另的线谱,另一观察艇相距一观察艇相距r处,用无指向性水听器作被动测听,设备的检测阈为处,用无指向性水听器作被动测听,设备的检测阈为8分贝。观察艇分贝。观察艇在在500Hz的噪声级为的噪声级为-34分贝,求在距离分贝,求在距离r上有多大传播衰减时才能可靠测听?上有多大传播衰减时才能可靠测听?设该设备作为被动测听,所用声呐方程为被动声呐方程。设该设备作为被动测听,所用声呐方程为被动声呐方程。SL-TL-NLGSGTDT可以得到:可以得到:TL=

56、-(DT-SL-GS(DI)+NL-GT)这里:这里:SL=60dB,DT=8dB,GS(DI)=0,NL=-34dB,GT=0SL=60dB,DT=8dB,GS(DI)=0,NL=-34dB,GT=0可以得到:可以得到:TL=-(8-60-0-34-0)=86dB结果表明单程传播损失小于结果表明单程传播损失小于86dB才可能可靠测听。才可能可靠测听。42例三:声呐作用距离估计例三:声呐作用距离估计 利用声呐方程可以估算声呐的作用距离,但要注意,声呐方程中有一部分参数利用声呐方程可以估算声呐的作用距离,但要注意,声呐方程中有一部分参数是与海洋环境有关的,同一部声呐(也就是是与海洋环境有关的,同

57、一部声呐(也就是GSGT一样)在不同的海区和不同一样)在不同的海区和不同季节,作用距离可能相差很大。计算作用距离时,有两种结果一样但方法不同的季节,作用距离可能相差很大。计算作用距离时,有两种结果一样但方法不同的估计。估计。第一种是由传播损失计算距离,另一种是把传播损失由近及远画出,用第一种是由传播损失计算距离,另一种是把传播损失由近及远画出,用FOM求声呐距离。求声呐距离。传播损失一般可以通过传播损失一般可以通过理论估计理论估计或或实际测量得到实际测量得到。下图。下图1给出给出同一声呐在不同水文条件下的作用距离,可见与海域、频率有关,图同一声呐在不同水文条件下的作用距离,可见与海域、频率有关

58、,图2给出典型给出典型估计作用距离的图(估计作用距离的图(500Hz)。)。曲线曲线A、B、C由海上测得,由该图可以查处不同由海上测得,由该图可以查处不同FOM下声呐作用距离。下声呐作用距离。500Hz43假设一部舰用主动声呐工作频率为假设一部舰用主动声呐工作频率为500Hz,发射功率发射功率1kW,发射指向性指数发射指向性指数DIT=15dB,接收机空间增益接收机空间增益GS12dB,时间增益时间增益GT=8dB,目标强度目标强度TS=15dB,检检测阈测阈DT=6dB,声呐平台噪声大于环境噪声为声呐平台噪声大于环境噪声为NL=78dB,那么声呐的作用距离?那么声呐的作用距离?首先计算声呐声

59、源级:首先计算声呐声源级:那么得到声呐的优质因素那么得到声呐的优质因素FOM:求距离:求距离:估算方法估算方法1:假设传播是理论上无吸收按球面波扩展,:假设传播是理论上无吸收按球面波扩展,TL20logr,根据根据FOM可以得到作用距离:可以得到作用距离:估算方法估算方法2:假设是在上图的海域中工作,则可以查图得到,:假设是在上图的海域中工作,则可以查图得到,A海域海域r=34kmB海域海域r=22km,C海域海域r=15km441.9.13 1.9.13 使用声呐方程应该注意问题使用声呐方程应该注意问题 声呐方程在声呐设计和对声呐性能预报中起到重要作用,但使用时以下是需声呐方程在声呐设计和对

60、声呐性能预报中起到重要作用,但使用时以下是需要注意的问题。要注意的问题。1、声呐方程是用声强度来描述的,含有能量平均的意思。因此只具有平均意义,声呐方程是用声强度来描述的,含有能量平均的意思。因此只具有平均意义,无法描述声呐的统计特性,而声呐系统恰恰是工作在复杂多变、随机环境中,因无法描述声呐的统计特性,而声呐系统恰恰是工作在复杂多变、随机环境中,因此声呐方程不能精确描述现代声呐的特性。此声呐方程不能精确描述现代声呐的特性。 由声学知识可以知道,声强是声能流在某一时间间隔内的平均值:由声学知识可以知道,声强是声能流在某一时间间隔内的平均值:式中,式中,T是时间间隔,是时间间隔,p,u分别是介质

61、中的声压和介质质点的振速。但是,我们知道,分别是介质中的声压和介质质点的振速。但是,我们知道,当声源发射的声信号是很短的脉冲信号时由于介质的传播效应和目标反射的物理效当声源发射的声信号是很短的脉冲信号时由于介质的传播效应和目标反射的物理效应、接收到的回波声信号波形会产生严重畸变,式所示的平均值会得到不确定的结果应、接收到的回波声信号波形会产生严重畸变,式所示的平均值会得到不确定的结果所以,该式就不再适用。作为一种常用的近似,可以在时间区间所以,该式就不再适用。作为一种常用的近似,可以在时间区间T内对声波的能流内对声波的能流密度求平均而得到声强,即将上式改写成:密度求平均而得到声强,即将上式改写

62、成:45值得注意的是,对于长脉冲声纳,区间值得注意的是,对于长脉冲声纳,区间T即为发射脉冲宽度,而且回波脉冲宽度也大即为发射脉冲宽度,而且回波脉冲宽度也大体等于此值;但对于短的瞬变脉冲声纳信号,区间体等于此值;但对于短的瞬变脉冲声纳信号,区间T一般是不确定的,回声宽度与发一般是不确定的,回声宽度与发射宽度相差甚大。射宽度相差甚大。尤立克研究了短脉冲信号下声纳方程是否还尤立克研究了短脉冲信号下声纳方程是否还适用的问题,指出声纳方程的强度形式仍然适用的问题,指出声纳方程的强度形式仍然适用,只要将声源级改写成如下形式:适用,只要将声源级改写成如下形式:对于时间区间对于时间区间T0发射恒定声源级发射恒定声源级SL的矩形脉的矩形脉冲,脉冲能量密度是平均声强与其脉冲宽度的冲,脉冲能量密度是平均声强与其脉冲宽度的乘积:乘积:供声呐方程使用的有效声源级为:供声呐方程使用的有效声源级为:为回波宽度,为回波宽度, 462、声呐总是工作在有限带宽内的,而且往往是窄带,因为窄带滤波是抗干扰、声呐总是工作在有限带宽内的,而且往往是窄带,因为窄带滤波是抗干扰的最为有效手段。所以声呐参数一般定义在的最为有效手段。所以声呐参数一般定义在1Hz带宽内,当接收系统具有一定带宽内,当接收系统具有一定带带宽时要进行带宽修正,即加宽时要进行带宽修正,即加 。47

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 医学/心理学 > 基础医学

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号