A critical safety wave pressure model of typical fishes under the action of underwater blasting shock waves
-
摘要: 为了探究水下爆破冲击波对鱼类损伤的影响,首先通过理论分析探讨了水下冲击波在鱼体的传播过程和对鱼鳔的损伤机理,构建了典型有鳔鱼类临界安全波压模型,并结合实测数据、参考相关文献,确定了模型参数。结果表明,典型有鳔鱼类临界安全波压与鱼体长呈线性关系。其中,水介质和鱼鳔壁介质的波阻抗比、鱼鳔宽度因数、鱼鳔壁厚因数、鱼鳔径向临界拉应力因数和鱼鳔形状因数分别为0.3~2.0、0.04~0.09、0.002、60和0.6~1.1。然后,根据典型有鳔鱼类临界安全波压模型参数,得到最大和最小鱼类临界安全波压模型分别为pic,max=30L和pic,min=3L(pic,max的单位为kPa,L的单位为cm)。最后,通过鱼类损伤情况的实测数据和文献数据,验证了鱼类临界安全波压模型。结果表明,不同体长的鱼类在不同冲击波压力时的受损状况分布,基本与鱼类所能承受的最大和最小的临界安全波压范围符合。并根据鱼类临界安全波压与体长的关系,划分了死亡区、存活有影响区和存活无影响区。Abstract: To investigate the propagation process of the underwater blasting shock wave in a fish body and its effect on typical swim bladder fishes, a critical safety wave pressure model for typical fishes was established and verified through theoretical analysis and field tests. According to the transverse reflection pattern of one-dimensional elastic compression wave between different media, the relationship between the critical safety wave pressure and the body length of typical swim bladder fishes was established. The length and mechanical properties of the swim bladder and fish body were measured using a vernier caliper, a digital micrometer, and a microcomputer tensile tester. Based on the measured data, the positive correlations of the length, width, wall thickness, and radial critical tensile stress of the swim bladder with the fish body length were determined, and the parameters in the fish critical safety wave pressure model were calibrated. The wave impedance ratio of the water medium and swim bladder wall medium was 0.3–2.0. The width, wall thickness, shape, and radial critical tensile stress coefficients of the swimming bladder were 0.04–0.09, 0.002, 0.6–1.1 and 60, respectively. The underwater blasting shock wave pressure and its effect on the fishes were measured using a blast wave tester, and the damage to the fishes was divided into three types: death, survival with influence, and survival without influence. The fish critical safety wave pressure model was verified by the statistical results of fish damage. The results show that the damage states of fishes with different body lengths at different shock pressures are conformed with the maximum and minimum critical safety wave pressure that the fishes can withstand. The proposed fish critical safety wave pressure model can be used to describe the relationship between the critical safety wave pressure and body length of the swim bladder fishes under the action of the underwater blasting shock waves. The research achievement can provide a theoretical basis for ecological protection of the fishes in the waterway regulation project of the upper reaches of the Yangtze River.
-
长江上游黄金航道是长江经济带综合交通运输体系的重要组成部分。在西部大开发和成渝双城经济圈等国家战略背景下,目前,正论证和实施朝天门至涪陵段航道、宜宾至重庆段航道的等级提升工程。其中,向家坝水电站坝轴线下1.8 km至重庆长江马桑溪大桥江段属于长江上游珍稀特有鱼类国家级自然保护区,重庆广阳镇至涪陵南沱镇江段属于长江重庆段四大家鱼国家级水产种质资源保护区[1-3]。在长江上游航道整治工程中,多采用水下钻孔爆破方式对巨量碍航礁石予以清除。而在航道整治工程水下爆破清礁时,传统爆破施工技术和工艺等会对水域生态环境造成一定影响,尤其爆破产生的冲击波可能会对鱼类造成损伤[4-6]。因此,在长江大保护战略下,开展水下爆破冲击波对长江上游典型鱼类损伤的研究具有重要意义。
对于有鳔鱼类,在水下爆破瞬间形成的冲击波从水体传播到鱼体时,因鱼体密度与周围水体密度近似,可认为冲击波在到达鱼体与水交界面时,以全透射形式直接通过鱼体向前传播[7]。当冲击波从鱼体传播至鱼鳔壁时,因鱼体与鱼鳔壁的波阻抗不同,在鱼体与鱼鳔壁的交界面瞬间会发生一次反射及透射,透射后的冲击波在鱼鳔壁中传播[8-10]。当透射后的冲击波从鱼鳔壁传播至鱼鳔内部时,在鱼鳔壁与鱼鳔内部的交界面瞬间发生二次反射及透射,因鱼鳔壁与鱼鳔内空气的波阻抗相差较大,鱼鳔内部二次透射后的冲击波峰值压力相对较小。第一次与第二次透射后的冲击波会在鱼鳔内、外形成较大压强差,在冲击波传播方向上对鱼鳔壁产生挤压,使鱼鳔变形,并在鱼鳔壁面形成影响鱼鳔的径向拉应力,一旦径向拉应力超过了鱼鳔自身的临界拉应力,鱼鳔就会受损而导致鱼类死亡。此外,鱼体中肝、脾、肾等器官也易受冲击波震动损伤,导致鱼类死亡[11-16]。
对于无鳔鱼类,水下爆破冲击波主要对其肝脏器官有损伤,相比于有鳔鱼类,冲击波对无鳔鱼类的损伤较轻[7]。国家标准GB 6722—2014《爆破安全规程》[17]中,给出了水下爆破冲击波作用下不同敏感程度鱼类的安全控制标准,然而,对不同体长有鳔鱼类的临界安全波压控制标准还缺乏系统研究。
本文中,针对长江上游典型有鳔鱼类,通过理论分析阐述水下爆破冲击波在鱼体的传播过程及对鱼鳔的损伤机理,进一步构建典型有鳔鱼类的临界安全波压模型。通过模型参数的测定和试验,验证所建模型能较好地反映水下爆破冲击波作用下鱼类临界安全波压与鱼体长的相互关系。
1. 模型构建
鱼体介质和鱼鳔壁介质在交界面处始终保持接触,根据介质连续条件和牛顿第三定律,水下爆破冲击波作用在鱼体与鱼鳔壁交界面上的瞬间,经反射、透射后,界面两侧质点速度相等、波压相等:
vi+vr=vt (1) pi+pr=pt (2) 式中:p为质点波压,v为质点速度,下标i、r和t分别表示入射波、反射波和透射波。
根据冲击波阵面动量守恒,水下爆破冲击波在传播过程中,质点波压p与介质密度ρ0、波速c0和质点速度v的关系为:
p=−ρ0c0v,c0=√E/ρ0 (3) 式中:E为介质弹性模量,ρ0c0为介质波阻抗。
根据波阵面动力学关系,联立式(1)~(3),可得:
piρwcw−prρwcw=ptρfcf (4) 式中:下标w、f分别表示水介质和鱼鳔壁介质。
结合式 (2)、(4),整理得到透射波压与入射波压的关系:
pt=21+npi (5) 式中:n为水介质和鱼鳔壁介质的波阻抗比。同样,若取n为鱼鳔壁介质和空气介质的波阻抗比,可得一次透射后冲击波从鱼鳔壁传播到鱼鳔内部的二次透射波压ptt。
构建典型有鳔鱼类临界安全波压模型的核心在于,影响鱼鳔的径向拉应力刚好达到鱼鳔自身所能承受的临界拉应力,以此时传播到鱼体表面的水下冲击波压pi为鱼类的临界安全波压。
鱼鳔因鱼的种类不同而形态各异,为使模型计算简便,假设鱼鳔形状为椭球体。为减小鱼鳔形状简化后产生的偏差,引入鱼鳔形状因数α(鱼鳔高度为直径的圆周长与椭球体竖向截面周长的比),如图1所示。鱼鳔临界拉应力为:
σ=α(pt−ptt)d4δ,α=πhπh+2(d−h) (6) 式中: d、δ和h为鱼鳔的宽度、壁厚和高度。
水下冲击波在传播过程中会逐渐衰减形成声波,因此介质的波阻抗可用声阻抗表示,根据文献[18]的声阻抗表得到水介质和空气介质的波阻抗为1.48和0.4 g/(m2·s)。鱼鳔壁介质的平均密度为1.24 kg/m3[19],根据鲫鱼、草鱼和鲢鱼的力学拉伸试验数据,鱼鳔壁介质的弹性模量范围为0.51~24.5 MPa,得到鱼鳔壁介质的波阻抗范围为0.78~5.51。则可得水介质和鱼鳔壁介质的波阻抗比范围、鱼鳔壁介质和空气介质的波阻抗比范围分别为0.3~2.0、2500~13775。鱼鳔壁介质和空气介质的波阻抗比较大,致使二次透射后的冲击波峰值压力较小,在构建鱼类临界安全波压模型时,可不考虑二次透射后的冲击波对鱼鳔临界拉应力的影响。
结合式(5)~(6),可得:
σ=αpid2δ(1+n) (7) 鱼鳔宽度、鱼鳔壁厚和鱼鳔临界拉应力与鱼体长分别存在对应关系[20-23],引入鱼鳔宽度因数β1及指数m1、鱼鳔壁厚因数β2及指数m2和鱼鳔临界拉应力因数β3及指数m3,鱼鳔宽度、鱼鳔壁厚和鱼鳔临界拉应力与鱼体长L的关系分别为:
d=β1Lm1,δ=β2Lm2,σ=β3Lm3 (8) 将式(8)代入式(7),可建立水下爆破冲击波作用下典型有鳔鱼类的临界安全波压模型:
pic=2(1+n)β2β3αβ1Lm2+m3−m1 (9) 2. 模型参数
2.1 参数的测定方法
选用游标卡尺测量鱼鳔长度、鱼鳔宽度和鱼鳔高度,选用卷尺测量鱼体长,选用电子数显千分尺测量鱼鳔壁厚,每组数据共测3次后取均值。
鱼鳔临界拉应力主要反映胶原纤维性质,与胶原纤维排列方向一致时力学强度大,而鱼鳔胶原纤维主要沿周向排列,因此鱼鳔周向的临界拉应力大于鱼鳔径向的临界拉应力。水下爆破冲击波对鱼鳔的破坏主要为径向拉应力,因此本文中分析鱼鳔的径向临界拉应力。采用力学拉伸试验测量鱼鳔的临界拉应力和弹性模量等[20],试验在电脑控制拉力试验机(STS5000N)上进行。
2.2 参数的测量结果
通过对鲫鱼(Carassius auratus)、鲢鱼(Hypophthalmichthys molitrix)和草鱼(Ctenopharyngodon idella)的体长及其鱼鳔的长度、宽度、厚度和径向临界拉应力的测量,确定水下爆破冲击波作用下典型有鳔鱼类临界安全波压模型中的4组相关参数(鱼鳔宽度因数及指数、鱼鳔壁厚因数及指数、鱼鳔临界拉应力因数及指数和鱼鳔形状因数),结果见表1。
表 1 临界安全波压模型相关参数的实测数据Table 1. Measured parameters related to the critical safety wave pressure model鱼种类 鱼体长/cm 鱼鳔长/cm 鱼鳔宽/cm 鱼鳔高/cm 鱼鳔壁厚/cm 鱼鳔径向临界拉应力/kPa 鱼鳔形状因数 鲫鱼 12.6 2.7 1.2 1.0 0.0141 591 0.89 14.5 2.8 1.3 1.2 0.0189 423 0.95 16.2 3.4 1.3 1.3 0.0164 457 1.00 19.4 3.4 1.3 1.1 0.0187 523 0.90 22.9 3.5 2.0 1.2 0.0191 654 0.70 24.6 3.6 2.1 1.7 0.0193 628 0.87 28.6 4.0 2.3 1.5 0.0232 811 0.75 31.5 4.6 2.4 2.2 0.0256 1021 0.95 31.7 4.7 2.4 2.3 0.0260 953 0.97 35.5 5.6 3.2 1.8 0.0287 1625 0.67 鲢鱼 42.6 6.4 3.5 3.7 0.0518 1100 1.04 46.3 7.7 3.8 3.3 0.0591 1160 0.91 49.0 8.5 4.2 3.6 0.0660 1240 0.90 51.5 8.8 3.8 2.4 0.0709 2210 0.73 51.7 9.3 4.4 4.7 0.0718 2230 1.04 56.4 10.2 4.2 2.1 0.0756 2440 0.61 58.0 10.8 4.6 2.6 0.0771 2500 0.67 59.3 11.2 4.7 3.6 0.0812 2660 0.84 61.0 11.4 5.8 3.5 0.0832 3450 0.70 63.1 11.9 5.9 3.2 0.0826 3210 0.65 草鱼 36.6 7.7 1.8 1.7 0.0476 830 0.96 38.0 8.2 2.0 1.1 0.0501 1960 0.66 39.3 8.5 1.9 1.6 0.0514 1390 0.89 39.5 9.9 2.2 2.0 0.0585 1340 0.94 40.6 10.4 2.6 2.2 0.0614 1440 0.90 41.0 11.0 2.2 1.7 0.0537 1960 0.84 41.4 8.4 1.9 2.2 0.0621 1890 1.10 50.2 12.1 2.3 2.3 0.0776 2430 1.00 59.1 14.0 3.0 2.6 0.0842 2560 0.91 60.1 14.4 3.2 2.1 0.0714 2840 0.75 2.3 模型参数的确定
Sobradillo等[23]对58条穆氏暗光鱼的形态研究显示,鱼鳔长、鱼鳔体积、鱼鳔等效球体半径与鱼体长呈显著正相关性。本文中,综合前人关于鱼鳔长与鱼体长的研究结论,通过测量鲫鱼、鲢鱼和草鱼总计30组鱼鳔长与鱼体长(见表1),拟合得到鲫鱼、鲢鱼和草鱼3种鱼类鱼鳔长与鱼体长的线性关系,如图2所示,鱼鳔长与鱼体长的因数为0.212。
Sobradillo等[23]也指出,鱼鳔宽与鱼鳔长呈正相关性。草鱼等鱼鳔形状偏细长,鲫鱼和鲢鱼等鱼鳔形状偏椭圆[24],本文中综合前人关于鱼鳔宽与鱼鳔长的研究结论,通过测量鲫鱼、鲢鱼和草鱼总计30组鱼鳔长与鱼鳔宽(见表1),拟合得到鱼鳔宽与鱼鳔长的两种线性关系,鱼鳔宽与鱼鳔长的因数为0.186~0.441,如图3所示。结合鱼鳔长度和鱼体长的关系,建议取鱼鳔宽度因数β1=0.04~0.09,指数m1=1。通过测量鲫鱼、鲢鱼和草鱼总计30组鱼鳔宽与鱼鳔高数据(见表1),代入鱼鳔形状因数公式,得到鱼鳔形状因数α=0.6~1.1。
Fine等[22]对43条蟾鱼的鱼鳔壁厚结构研究显示,成年蟾鱼鱼鳔壁厚与鱼体长呈正线性相关。本文中,综合前人关于鱼鳔壁厚与鱼体长的研究结论,通过测量鲫鱼、鲢鱼和草鱼总计30组鱼鳔壁厚(见表1),拟合得到鱼鳔壁厚与鱼体长的线性关系,如图4所示。建议取鱼鳔壁厚因数β2=0.002,指数m2=1。
通过测定鲫鱼、鲢鱼和草鱼总计30组鱼鳔径向临界拉应力(见表1),拟合得到鱼鳔径向临界拉应力与鱼体长的线性关系y=60x−570,如图5所示。建议取鱼鳔径向临界拉应力因数β3=60,指数m3=1。
3. 模型验证
3.1 试验
基于水下爆破冲击波作用下鱼类损伤试验,依托重庆市果园港重大件泊位港池开挖爆破工程,验证典型有鳔鱼类临界安全波压模型。受试鱼种为长6~48 cm的青鱼(Mylopharyngodon piceus)、草鱼、鲢鱼和鳙鱼(Aristichthys nobilis),将鱼类放置于鱼笼中,使用受力绳将鱼笼固定于爆破工作船的船头、船中部和船尾的船舷处,每个鱼笼各放置12~20条鱼,分别距爆源50、70和90 m,现场爆破药量为22~138 kg。采用爆破冲击波测试仪(Blast PRO)、压力传感器(TP-SJB-5M)及水下爆破专用信号线缆,监测水下爆破冲击波压力,测点布置如图6所示。
3.2 结果和讨论
将鱼类损伤情况分为3种(鱼类存活、无影响,鱼类存活、有影响,鱼类死亡),并按照鱼群死亡率对有鳔鱼类损伤结果进行判定。当鱼群死亡率为0,判定为鱼类存活、无影响;当鱼群死亡率低于10%,判定为鱼类存活、有影响;当鱼群死亡率高于10%,判定为鱼类死亡。鱼类损伤情况见表2,表中也收集了文献[11]中的鲫鱼损伤数据。
表 2 有鳔鱼类的损伤Table 2. Damages for swim bladder fishes鱼种类 鱼体长/cm 压力/kPa 损伤结果 判定结果 草鱼 30~48 81 鱼类存活,无变化 鱼类存活,无影响 6~15 514 鱼类存活,但行动迟缓 鱼类存活,有影响 鲢鱼 15~30 150 鱼类存活,但行动迟缓 鱼类存活,有影响 16~35 2000 60%受到严重损伤当即死亡,
40%受到轻微损伤在爆破结束后一小时内逐渐死亡鱼类死亡 青鱼、草鱼、鲢鱼、鳙鱼 6~18 706 青鱼全部死亡,草鱼约75%死亡,鲢鱼和鳙鱼约50%死亡 鱼类死亡 鲫鱼[11] 15~40 927 鱼类存活,但行动迟缓 鱼类存活,有影响 1193 33%鱼类死亡 鱼类死亡 2611 鱼类全部死亡 鱼类死亡 将已确定的模型参数代入式(9)。取最大波阻抗比为2,最小鱼鳔宽度因数为0.04,最小鱼鳔形状因数为0.6,得到典型有鳔鱼类最大临界安全波压模型为:pic,max=30L。取最小波阻抗比为0.3,最大鱼鳔宽度因数为0.09,最大鱼鳔形状因数为1.1,得到典型有鳔鱼类最小临界安全波压模型为:pic,min=3L。其中,pic,max的单位为kPa,L的单位为cm。
典型有鳔鱼类最大和最小临界安全波压模型曲线将冲击波对鱼类的影响划分为3个区域,用鱼类损伤结果(见表2)来验证典型有鳔鱼类的临界安全波压模型。结果表明,不同体长的鱼类在不同冲击波压力时的受损状况分布基本与鱼类所能承受的最大和最小的临界安全波压范围符合。根据典型有鳔鱼类临界安全波压模型结果与鱼类损伤结果对比(见图7),典型有鳔鱼类临界安全波压模型曲线较好地描述了水下爆破冲击波作用下鱼类临界安全波压与体长的关系,可划分为:鱼类死亡区,鱼类存活、有影响区,鱼类存活、无影响区。
然而,典型有鳔鱼类最大临界安全波压模型的结果与试验结果仍有稍许偏差。由图7可见,体长范围较小且存活的部分鲢鱼和鲫鱼被划入了死亡区域,体长范围较大且死亡的部分鲫鱼被划入了存活区域。这主要原因是,因本文中以鱼群死亡率10%作为鱼类存活或死亡的临界判定依据,致使模型结果出现部分差异。另外,本文中通过理论研究将复杂的水下爆破过程抽象化,最终归纳为鱼类临界安全波压与体长的数学模型,在建模过程中忽略了鱼类的其余器官损伤情况、身体结构等的影响,造成了一定的模型误差。对图7偏差较大的数据,采取方差分析方法进行有效性检验,计算的P为0.03 (小于0.05),表明个别存在偏差的数据对模型整体的准确性影响较小。
综上所述,本文中通过理论分析得到鱼类临界安全波压与体长的线性关系,以此构建了典型有鳔鱼类临界安全波压模型,并通过不同体长的鱼类在不同冲击波压力时的受损情况,验证了模型的准确性。研究结果可为长江上游航道整治工程中有鳔鱼类的生态保护工作提供理论参考。
4. 结 论
基于水下爆破冲击波的传播过程及其作用下的鱼类损伤机理,构建了鱼类临界安全波压模型。
(1) 由水下爆破冲击波对鱼类损伤过程的理论分析,获得了典型有鳔鱼类临界安全波压与体长的关系模型。由模型参数实测数据,分析了鱼鳔长、鱼鳔宽、鱼鳔壁厚、鱼体长与鱼鳔径向临界拉应力的相互正相关性,并率定了模型参数n、β1、β2、β3和α分别为0.3~2.0、0.04~0.09、0.002、60和0.6~1.1。
(2) 基于不同体长的鱼类在不同冲击波压力下的损伤,将鱼类的损伤划分为死亡、存活有影响、存活无影响等3种情况。进一步对鱼类临界安全波压模型进行了验证,鱼类的受损状况基本符合鱼类所能承受的最大和最小临界安全波压范围。
构建了典型有鳔鱼类临界安全波压模型,测定了模型相关参数,并通过现场试验对该模型进行了验证。但不同鱼类的模型参数有一定差异,在应用于水下爆破的鱼类安全防护实践中,需要针对不同鱼种进一步细化相关模型参数,增强模型的实际应用效果。
-
表 1 临界安全波压模型相关参数的实测数据
Table 1. Measured parameters related to the critical safety wave pressure model
鱼种类 鱼体长/cm 鱼鳔长/cm 鱼鳔宽/cm 鱼鳔高/cm 鱼鳔壁厚/cm 鱼鳔径向临界拉应力/kPa 鱼鳔形状因数 鲫鱼 12.6 2.7 1.2 1.0 0.0141 591 0.89 14.5 2.8 1.3 1.2 0.0189 423 0.95 16.2 3.4 1.3 1.3 0.0164 457 1.00 19.4 3.4 1.3 1.1 0.0187 523 0.90 22.9 3.5 2.0 1.2 0.0191 654 0.70 24.6 3.6 2.1 1.7 0.0193 628 0.87 28.6 4.0 2.3 1.5 0.0232 811 0.75 31.5 4.6 2.4 2.2 0.0256 1021 0.95 31.7 4.7 2.4 2.3 0.0260 953 0.97 35.5 5.6 3.2 1.8 0.0287 1625 0.67 鲢鱼 42.6 6.4 3.5 3.7 0.0518 1100 1.04 46.3 7.7 3.8 3.3 0.0591 1160 0.91 49.0 8.5 4.2 3.6 0.0660 1240 0.90 51.5 8.8 3.8 2.4 0.0709 2210 0.73 51.7 9.3 4.4 4.7 0.0718 2230 1.04 56.4 10.2 4.2 2.1 0.0756 2440 0.61 58.0 10.8 4.6 2.6 0.0771 2500 0.67 59.3 11.2 4.7 3.6 0.0812 2660 0.84 61.0 11.4 5.8 3.5 0.0832 3450 0.70 63.1 11.9 5.9 3.2 0.0826 3210 0.65 草鱼 36.6 7.7 1.8 1.7 0.0476 830 0.96 38.0 8.2 2.0 1.1 0.0501 1960 0.66 39.3 8.5 1.9 1.6 0.0514 1390 0.89 39.5 9.9 2.2 2.0 0.0585 1340 0.94 40.6 10.4 2.6 2.2 0.0614 1440 0.90 41.0 11.0 2.2 1.7 0.0537 1960 0.84 41.4 8.4 1.9 2.2 0.0621 1890 1.10 50.2 12.1 2.3 2.3 0.0776 2430 1.00 59.1 14.0 3.0 2.6 0.0842 2560 0.91 60.1 14.4 3.2 2.1 0.0714 2840 0.75 表 2 有鳔鱼类的损伤
Table 2. Damages for swim bladder fishes
鱼种类 鱼体长/cm 压力/kPa 损伤结果 判定结果 草鱼 30~48 81 鱼类存活,无变化 鱼类存活,无影响 6~15 514 鱼类存活,但行动迟缓 鱼类存活,有影响 鲢鱼 15~30 150 鱼类存活,但行动迟缓 鱼类存活,有影响 16~35 2000 60%受到严重损伤当即死亡,
40%受到轻微损伤在爆破结束后一小时内逐渐死亡鱼类死亡 青鱼、草鱼、鲢鱼、鳙鱼 6~18 706 青鱼全部死亡,草鱼约75%死亡,鲢鱼和鳙鱼约50%死亡 鱼类死亡 鲫鱼[11] 15~40 927 鱼类存活,但行动迟缓 鱼类存活,有影响 1193 33%鱼类死亡 鱼类死亡 2611 鱼类全部死亡 鱼类死亡 -
[1] 张先炳, 杨胜发, 杨威, 等. 长江上游宜宾-江津与涪陵-丰都江段鱼类早期资源分布研究 [J]. 淡水渔业, 2021, 51(5): 51–59. DOI: 10.3969/j.issn.1000-6907.2021.05.007.ZHANG X B, YANG S F, YANG W, et al. The distribution of the early-stage fish resources between Yibin-Jiangjin and Fuling-Fengdu in the upper reaches of the Yangtze River [J]. Freshwater Fisheries, 2021, 51(5): 51–59. DOI: 10.3969/j.issn.1000-6907.2021.05.007. [2] 高天珩, 田辉伍, 叶超, 等. 长江上游珍稀特有鱼类国家级自然保护区干流段鱼类组成及其多样性 [J]. 淡水渔业, 2013, 43(2): 36–42. DOI: 10.3969/j.issn.1000-6907.2013.02.007.GAO T H, TIAN H W, YE C, et al. Diversity and composition of fish in the mainstream of national nature reserve of rare and endemic fish in the upper Yangtze River [J]. Freshwater Fisheries, 2013, 43(2): 36–42. DOI: 10.3969/j.issn.1000-6907.2013.02.007. [3] 段辛斌. 长江上游鱼类资源现状及早期资源调查研究 [D]. 武汉: 华中农业大学, 2008: 1−6. DOI: 10.7666/d.y1598376. [4] 喻灿星, 覃国杰, 曾丽. 内河航道水下炸礁钻孔爆破技术研究 [J]. 工程爆破, 2021, 27(4): 58–63. DOI: 10.19931/j.EB.20200273.YU C X, QIN G J, ZENG L. Research on technology of underwater reef drilling blasting in an inland waterway [J]. Engineering Blasting, 2021, 27(4): 58–63. DOI: 10.19931/j.EB.20200273. [5] 李金河, 赵继波, 谭多望, 等. 炸药水中爆炸的冲击波性能 [J]. 爆炸与冲击, 2009, 29(2): 172–176. DOI: 10.11883/1001-1455(2009)02-0172-05.LI J H, ZHAO J B, TAN D W, et al. Underwater shock wave performances of explosives [J]. Explosion and Shock Waves, 2009, 29(2): 172–176. DOI: 10.11883/1001-1455(2009)02-0172-05. [6] 赵根, 吴从清, 王文辉. 爆破水中冲击波对鱼类损伤研究 [J]. 工程爆破, 2011, 17(4): 103–105, 93. DOI: 10.3969/j.issn.1006-7051.2011.04.025.ZHAO G, WU C Q, WANG W H. Research on blasting shock wave in water to damage of fish [J]. Engineering Blasting, 2011, 17(4): 103–105, 93. DOI: 10.3969/j.issn.1006-7051.2011.04.025. [7] 李文涛, 张秀梅. 水下爆破施工对鱼类影响的估算及预防措施 [J]. 海洋科学, 2003, 27(11): 20–23. DOI: 10.3969/j.issn.1000-3096.2003.11.005.LI W T, ZHANG X M. Impact and mitigation measures for fish communities exposed to underwater explosion [J]. Marine Sciences, 2003, 27(11): 20–23. DOI: 10.3969/j.issn.1000-3096.2003.11.005. [8] 贾虎, 沈兆武. 空气隔层对水中冲击波的衰减特性 [J]. 爆炸与冲击, 2012, 32(1): 61–66. DOI: 10.11883/1001-1455(2012)01-0061-06.JIA H, SHEN Z W. An investigation into attenuation of underwater shockwave by air interlayer [J]. Explosion and Shock Waves, 2012, 32(1): 61–66. DOI: 10.11883/1001-1455(2012)01-0061-06. [9] 樊自建, 沈兆武, 马宏昊. 水中空气隔层对冲击波传播衰减作用的初步探讨 [J]. 工程爆破, 2007, 13(2): 7–10. DOI: 10.3969/j.issn.1006-7051.2007.02.002.FAN Z J, SHEN Z W, MA H H. Primary study on attenuation of underwater shock wave by using air insulation [J]. Engineering Blasting, 2007, 13(2): 7–10. DOI: 10.3969/j.issn.1006-7051.2007.02.002. [10] 徐爽, 赵宁, 王春武, 等. 水/气多介质问题的界面处理方法 [J]. 爆炸与冲击, 2015, 35(3): 326–334. DOI: 10.11883/1001-1455-(2015)03-0326-09.XU S, ZHAO N, WANG C W, et al. Interface treating methods for the gas-water multi-phase flows [J]. Explosion and Shock Waves, 2015, 35(3): 326–334. DOI: 10.11883/1001-1455-(2015)03-0326-09. [11] 尚龙生, 戴云丛, 刘现明, 等. 水中爆破对双台子河口渔场的影响 [J]. 海洋环境科学, 1994, 13(3): 23–26,32. [12] KEEVIN T M. A review of natural resource agency recommendations for mitigating the impacts of underwater blasting [J]. Reviews in Fisheries Science, 1998, 6(4): 281–313. DOI: 10.1080/10641269891314302. [13] 周杰, 陶钢, 王健. 爆炸冲击波对肺损伤的数值模拟 [J]. 爆炸与冲击, 2012, 32(4): 418–422. DOI: 10.11883/1001-1455(2012)04-0418-05.ZHOU J, TAO G, WANG J. Numerical simulation of lung injury induced by shock wave [J]. Explosion and Shock Waves, 2012, 32(4): 418–422. DOI: 10.11883/1001-1455(2012)04-0418-05. [14] WIERNICKI C J, LIANG D, BAILEY H, et al. The effect of swim bladder presence and morphology on sound frequency detection for fishes [J]. Reviews in Fisheries Science & Aquaculture, 2020, 28(4): 459–477. DOI: 10.1080/23308249.2020.1762536. [15] HALVORSEN M B, CASPER B M, WOODLEY C M, et al. Threshold for onset of injury in Chinook salmon from exposure to impulsive pile driving sounds [J]. PLoS One, 2012, 7(6): e38968. DOI: 10.1371/journal.pone.0038968. [16] 杨志焕, 朱佩芳, 蒋建新, 等. 水下冲击波的生物效应 [J]. 爆炸与冲击, 2003, 23(2): 134–139.YANG Z H, ZHU P F, JIANG J X, et al. Bio-effects of underwater blast waves [J]. Explosion and Shock Waves, 2003, 23(2): 134–139. [17] 中华人民共和国国家标准编写组. 爆破安全规程: GB 6722—2014 [S]. 北京: 中国标准出版社, 2015. [18] 孙玉. 声阻抗梯度变化材料中声传播特性研究 [D]. 黑龙江: 哈尔滨工程大学, 2015: 91–107. DOI: 10.7666/d.D749650. [19] 杜伟东, 李海森, 陈宝伟, 等. 一种基于声散射特性的有鳔鱼特征获取方法 [J]. 应用声学, 2014, 33(6): 505–511. DOI: 10.11684/j.issn.1000-310X.2014.06.005.DU W D, LI H S, CHEN B W, et al. Features acquisition of fish with swim bladder based on acoustic scattering characteristics [J]. Journal of Applied Acoustics, 2014, 33(6): 505–511. DOI: 10.11684/j.issn.1000-310X.2014.06.005. [20] 李鑫. 鱼鳔作为新型心血管外科手术生物材料的试验研究 [D]. 上海: 第二军医大学, 2013: 1–32. DOI: 10.7666/d.Y2339973. [21] 张更申, 张庆俊, 孙国柱, 等. 应用鲤鱼鳔进行家兔硬脑膜修补术实验研究 [J]. 河北医科大学学报, 2000, 21(6): 337–340. DOI: 10.3969/j.issn.1007-3205.2000.06.006.ZHANG Q S, ZHANG Q J, SUN G Z, et al. Experimental duraplasty with carp swim-bladder in rabbits [J]. Journal of Hebei Medical University, 2000, 21(6): 337–340. DOI: 10.3969/j.issn.1007-3205.2000.06.006. [22] FINE M L, KING T L, ALI H, et al. Wall structure and material properties cause viscous damping of swimbladder sounds in the oyster toadfish Opsanus tau [J]. Proceedings of the Royal Society B: Biological Sciences, 2016, 283(1841): 20161094. DOI: 10.1098/rspb.2016.1094. [23] SOBRADILLO B, BOYRA G, MARTINEZ U, et al. Target strength and swimbladder morphology of Mueller’s pearlside (Maurolicus muelleri) [J]. Scientific Reports, 2019, 9(1): 17311. DOI: 10.1038/s41598-019-53819-6. [24] 李彬寒. 鱼鳔源抗钙化心血管生物材料的研究 [D]. 北京: 北京协和医学院, 2020: 48–59. DOI: 10.27648/d.cnki.gzxhu.2020.000519. -