Spatial distribution characteristics analysis and resources assessment of fish in Lugu Lake
-
摘要: 2018年8月15—19日，利用Lowrance elite-7 Ti回声探测仪首次对云南省泸沽湖展开鱼类空间分布及资源量评估。调查采用“之”字形航线对泸沽湖全湖进行水声学探测，累计航程42.71 km，基于深度变化值将探测水体分为若干探测单元进行鱼类资源评估，并运用反距离插值法 (IDW) 基于各单元密度绘制进行空间分布图。水声学探测结果表明，泸沽湖鱼类空间分布不均匀，水平方向上，鱼类主要集群分布于近岸环形水域，其中小渔坝附近水域鱼类密度最高，湖心区域鱼类密度最低；垂直方向上，鱼类在表层水域分布较为集中；2018年8月泸沽湖鱼类密度为0.153尾·(1 000 m3)−1，数量24.39万尾，总资源量50.55 t；其中，鲫 (Carassius auratus) 21.17万尾、17.36 t，鲤 (Cyprinus carpio) 1.74万尾、23.01 t，草鱼 (Ctenopharyngodon idellus) 8.63万尾、8.63 t，大鳞副泥鳅 (Paramisgurnus dabryanus) 0.68万尾、1.33 t，小口裂腹鱼 (Schizothorax microstomus) 0.07万尾、0.22 t。此次调查未发现曾出现记录的宁蒗裂腹鱼 (S. ninglangensis) 和厚唇裂腹鱼 (S. labrosus)。Abstract: From August 15 to 19, 2018, we assessed the spatial distribution and resources of fish in Lugu Lake, Yunnan Province for the first time by using lowrange elite-7 Ti echo detector. Zigzag route was used to detect the whole Lugu Lake, with a total voyage of 42.71 km. Based on the depth variation value, the lake was divided into a number of detection units for fish resources evaluation, and the spatial distribution map was drawn based on the density of each unit by inverse distance weighting (IDW). The results of hydroacoustic survey shows that the spatial distribution of fish in Lugu Lake was uneven. In the horizontal direction, most fish clusters were distributed in the coastal annular waters, among which the fish density in the waters near Xiaoyuba was the highest, while the fish density in the middle of the lake was the lowest; in the vertical direction, the fish distribution in the surface waters was relatively concentrated. In August 2018, the fish density in Lugu Lake was 0.153·(1 000 m3)−1; the total number of fish was 243 900 and the total resource was 50.55 t (Carassius auratus: 211 700 ind., total mass of 17.36 t; Cyprinus carpio: 17 400 ind., total mass of 23.01 t; Ctenopharyngodon idellus: 86 300 ind., total mass of 8.63 t; Paramisgurnus dabryanus: 6 800 ind., total mass of 1.33 t; Schizothorax microstomus: 700 ind., total mass of 0.22 t). In this survey, there was no records of S. ninglangensis and S. labrosus which had been found in previous investigation.
-
Keywords:
- Fish /
- Resource assessment /
- Spatial distribution /
- Hydroacoustic survey /
- Lugu Lake
-
金属卤化物集鱼灯 (简称“金卤灯”) 是光诱渔业最常用的光源类型[1-6]。光诱渔船海面光场分布计算是集鱼灯应用研究的重要课题之一[3,7]。光场分布计算方法包括几何光学法[2,8-9]和蒙特卡罗 (Monte Carlo, MC) 模拟方法[10-12]。几何光学方法中，Bae等[3]将点光源模型 (Point model, PM) 与指数衰减模型结合计算了渔船海面照度，假定灯具空间光场分布各向同性。Choi等[13]提出了线性光源模型 (Line model, LM) 计算灯光船周围海面照度，因其假定灯组为连续发光体，计算值大于实际情况。肖启华和张丽蕊[14]提出了面光源模型，适用于灯具表面积大、计算点与灯具之间小的光场数值的计算。Lai等[15]应用球面模型对LED集鱼灯进行二次透镜设计，并计算船舷两侧的海面光场分布。钱卫国和王飞[9]、叶超[16]、侍炯等[17]提出了配光曲线模型 (Light distribution curve model, LDC) 并做了大量的海上实证研究。由于几何方法主要采用光在水汽界面的折射定律和在空气中直线传播规律，无法求解海面波动状态下的照度分布，学者开始引入蒙特卡罗方法求解光学传输问题[12,18-22]。张涤[12]应用蒙特卡罗方法研究了可见光水下信道问题，指出了不同波长光子信道衰减情况。官文江等[20]结合风浪斜率概率模型[23-24]，进一步应用蒙特卡罗提出了集鱼灯光场分布的计算框架。然而计算框架中的光子数取值问题和光子辐射模型还有待进一步完善，这两个环节决定了集鱼灯光场计算的准确度，进而影响到光场有效诱集范围、渔船之间合理作业间距等方面的估算。因此本文通过数值模拟方法讨论不同光子数取值对模拟结果稳定性的影响，同时根据灯具光度分布提出了新的灯具辐射模型并验证，结合新的辐射模型给出算例，为光诱渔船光场分布计算提供借鉴和参考。
1. 材料与方法
1.1 坐标系
为描述灯具和海面计算点的相对位置，建立
$ xyz $ 三维坐标系 (图1)。其中船尾至船首的纵剖面为$ xoz $ ，船尾中点在海面的投影点为原点$ o $ ，船舯线为$ x $ 轴 (船首为$x $ 正向)，船尾垂直于$ xoz $ 方向为$y$ 轴 (右舷为$ y $ 轴正向)，$ z$ 轴垂直于$ xoy $ 平面 (向上为正)。灯具垂直悬挂于点$ L(x, y, h)$ ，$ h$ 为灯具距离海面的高度，$P$ 点位于$ {xoy} $ 平面内 (图1-a)。定义光束LP在$xoy$ 平面的投影与$x$ 轴顺时针方向的夹角为投影旋转角$φ$ 。定义光束LP与$z$ 轴正方向的夹角为天底角$ {θ} $ 。海面为波动状态，海面与水平面夹角为β (图1-b)。![]()
图 1 灯具与计算点的相对位置及灯具安装示意Figure 1. Relative position between lamp and calculation point and installation diagram of lamp1.2 材料
选取目前光诱渔业中广泛应用的金属卤化物集鱼灯 (2 kW DCJ 2000TT型) 作为研究对象，选取目前光诱渔业中广泛应用的金属卤化物集鱼灯(2 kW DCJ 2000TT 型) 作为研究对象，集鱼灯额定功率2000 W，电压范围210~250 V，额定光通量220000 lm，灯头型号E39/79-A，灯具中心高度295 mm，全长465 mm，最大直径90 mm。
灯具光度分布数据在国家远洋渔业工程技术研究中心进行了测试，灯具配光在垂直方向和水平方向的配光曲线分布见图2。测试设备：远方光电GO-2 000光度分布测试仪，测试模式选系统内置的C−γ。
可以看出光强在水平和垂直剖面内的光强分布曲线存在一定差异，取C90/270和γ90平面内光强分布曲线，借助Matlab (R2015b) 软件用Fourier级数拟合灯具光强分布曲线，公式分别对应如下：
$$\begin{aligned} I(\theta)=& 19\;670-3\;366 \cos (1.851 \theta)+2\;197 \sin (1.851 \theta)-\\ & 1\;201 \cos (3.702 \theta)+2\;245 \sin (3.702 \theta)+\\ & 38.56 \cos (5.553 \theta)+1\;176 \sin (5.553 \theta)\left(R^{2} = 0.96\right) \end{aligned}$$ (1) $$\begin{aligned} g(\varphi)=& 15\;320-9\;206 \cos (1.962 \varphi)-396.3 \sin (1.962 \varphi)-\\ & 3\;534 \cos (3.924 \varphi)-238.8 \sin (3.924 \varphi)-\\ & 1\;305 \cos (5.886 \varphi)-167 \sin (5.886 \varphi)\left(R^{2} = 0.98\right) \end{aligned}$$ (2) 式中：
$ {θ} $ 为灯具垂直剖面γ90配光曲线，决定了灯具不同天底角光束的辐射强度；$ {φ} $ 为灯具水平剖面C90/270配光曲线，决定了灯具不同投影旋转角光束的辐射强度。1.3 模型建立
蒙特卡罗方法是研究光辐射传输特性的经典方法[24-26]，通过将光束能量离散成大量光子，利用随机抽样方法追踪光子路径来解决光束传输问题，计算流程参考官文江等[20]的方法。模型计算假定条件：1) 模拟中所有光子无波长差异；2) 光子在水气界面只考虑反射与入射，无其他形式衰减；3) 不考虑光子在空气中的吸收与散射作用。
1.3.1 光子辐射优化模型
官文江等[20]报道的光子辐射模型中投影旋转角 (
$ {φ} $ ) 使用均匀随机数生成，难以全面反映灯具辐射特征 (图1-a)。对此本文进一步优化，辐射模型考虑光子辐射天底角 ($ {θ} $ ) 和投影旋转角 ($ {φ} $ ) 两个方向的特征，计算公式：$$\left\{\begin{array}{l}F({{{θ}}})=\dfrac{\int_{0}^{{{{θ}}}_{{{i}}}} \mathrm{I}({{{θ}}}) \mathrm{d} {{{θ}}}}{\int_{0}^{2 \pi} \mathrm{I}({{{θ}}}) \mathrm{d} {{{θ}}}} \\ M(\varphi)=\dfrac{\int_{0}^{\varphi_{{{i}}}} \mathrm{~g}({\varphi) \mathrm{d} \varphi}}{\int_{0}^{2 \pi} \mathrm{g}(\varphi) \mathrm{d} \varphi}\end{array}\right.$$ (3) 其中：
$ F $ ($ {θ} $ )、$ {M}{(}{φ}{)} $ 分别为光子天底角$ {θ} $ 和投影旋转角$ {φ} $ 的归一化概率分布函数，考虑到集鱼灯投向天空的光源属于浪费资源，故$ {θ} $ 取值介于[0.5π, π]。$ {θ} $ 和$ {φ} $ 利用反函数法[25]产生。1.3.2 海浪斜率计算
光子从灯具出射后，沿直线投射至水面。而海面受风的影响形成了毛细海浪，因此交互界面有一定的斜率。根据Cox和Munk[27]观测，不同风速下的波浪面倾斜角经验概率密度公式：
$$\rho(\beta)=\frac{2}{\sigma^{2}} {\rm{e}}^{\frac{-\tan ^{2} \beta}{\sigma^{2}}} \tan \beta \sec ^{2} \beta$$ (4) 其中：σ2=0.003+0.005 12V，V为距海面10 m处的风速(m·s−1)；
$ {β} $ 为波浪面与水平面的夹角，同样利用反函数法产生伪随机数[25]确定波浪斜率。确定海浪斜率后即可确定法线方向，光束在海面的入射角
$ {{θ}}_{{i}} $ 调整为：$$ {{θ}}_{{i}}={θ}-{β} $$ (5) 式中：
$ {θ} $ 为光子从灯具出射的天底角；$ {β} $ 同上。1.3.3 海面反射系数
光束经过大气−海水界面发生折射和反射，总反射率Rt计算公式：
$$R_{\mathrm{t}}=R_{\mathrm{f}}+R_{\mathrm{u}}$$ (6) 式中：
$ {{R}}_{\rm{f}} $ 为镜面折射率；$ {{R}}_{\rm{u}} $ 水下反射率。丁明亮[10]指出$ {{R}}_{\rm{u}} $ 泡沫反射率远小于折射率$ {{R}}_{\rm{f}} $ ，因此忽略不计。光子调整角度后在水面反射与入射遵循Fresnel反射定律[28]：$$R_{{f}}=\left[\frac{{n}_{1} \cos {{θ}}_{{i}}-{n}_{2} \sqrt{1-\left(\dfrac{{n}_{1}}{{n}_{2}} \sin {{θ}}_{{i}}\right)^{2}}}{{n}_{1} \cos {{θ}}_{{i}}+{n}_{2} \sqrt{1-\left(\dfrac{{n}_{1}}{{n}_{2}} \sin {{θ}}_{{i}}\right)^{2}}}\right]^{2}$$ (7) 式中：
$ {{θ}}_{{i}} $ 为入射角；$ {{n}}_{{1}} $ 为空气折射率，取1；$ {{n}}_{{2}} $ 为海水折射率，取1.33。若$ {{θ}}_{{i}}{=0} $ ，则$ {{R}}_{{{\rm{f}}}} $ 按下式计算：$$R_{{{\rm{f}}}}=\left(\frac{{n}_{1}-{n}_{2}}{{n}_{1}+{n}_{2}}\right)^{2}$$ (8) 模拟过程中，依据轮盘赌法[26]进行判别：取随机数
$ {R} $ (为0~1内均匀分布随机数)，当$ {R}{ > }{{R}}_{{{\rm{t}}}} $ ，则光子落入水中，否则反射不予统计。对接收面进行栅格化处理后，通过统计每个栅格内的光子数能量与栅格面积之比得到照度值。1.4 模型验证
为验证模型准确性，在光学实验室内完成验证测试。将灯具安装在自动旋转台上，灯高1.6 m，功率2 kW，照度计型号Hobo Mx2202，分辨率为0.01 lx。灯具安装在旋转平台上，照度计放置于地面，与灯具中心水平距离1~5 m，间隔1 m放置一个照度计，见图3-a中的t1—t5位置。测试过程中，转台绕C轴转动，从而测得灯具不同投影旋转角
$ {φ} $ 所对应的地面照度值，实际测试灯具与照度计布置见图3-b。2. 结果
2.1 光子数取值对模拟结果的影响
通过设定不同的光子数值105~1010讨论模拟结果的稳定性，重复模拟1 000次，取模拟结果的差异系数 (Coefficient of variation, CV) 作为稳定性计量指标，差异系数越低稳定性越高。结果显示：1) 对同一计算点，光子数从105增至1010，模拟结果差异系数减少；以与灯具水平距离5 m的计算点为例，当光子数增至108，差异系数为0.07%，光子数取值增大，差异系数减少至0.05%； 2) 相同光子数取值条件下，距灯具水平距离越近的计算点，计算结果稳定性越高；随着光子数取值增大，稳定性差异缩小 (图4)。
![]()
图 4 光子数取值对模拟结果波动性的影响Figure 4. Effect of photon number on volatility of modeling results2.2 光子辐射模型优化
金属卤化物灯由于内部支架等遮挡影响，其空间光场分布在水平切面内并非饶轴对称 (图2-a)。针对新的光子辐射模型不同出射方向 (
$ {φ} $ =90º、60º、30º、0º) 的照度值进行计算并实测，结果显示新模型计算结果与实测相对误差分别为5.12%、4.72%、3.28%、5.03%，平均值为4.53%，未优化模型相对误差均值为5.23% (图5)。不同方向上的计算结果差异系数小于0.08%，模拟稳定性较好 (图6)。![]()
图 5 优化模型不同出射方向计算与实测比较Figure 5. Comparison of different light field calculation models![]()
图 6 优化模型不同出射方向计算的差异系数(光子数取值为108)Figure 6. Coefficient of variation of optimized model simulation results (Photons number of 108)作为对比，基于点光源模型的模拟结果与实际的误差分别为3.57%、4.57%、6.81%、17.1%，平均值8.41%，显然误差增大；基于配光曲线模型计算结果与实际的误差分别为7.07%、8.1%、10.42%、21.04%，平均值11.65%，误差最大。由上可知，本文提出的优化光子输出模型，能够较好地反映灯具光场的空间分布情况。
2.3 优化模型应用实例
结合本文提出的光子辐射优化模型，提供了光诱渔船在实际作业海域的光场分布算例，光诱渔船集鱼灯配置等参数见表1。风速取2 m·s−1，数据来自Noaa Erddap。实测值参照Choi等[13]的研究。
表 1 模型参数表Table 1. Model parameters项目 Item 数值 Value 项目 Item 数值 Value 船舶吨位 Tonnage 19.9 金卤灯功率 Metal halide lamp power/kW 2.0 船长 Length/m 22.0 金卤灯光通量 Luminous flux/lm 200 000 船宽 Width/m 4.8 灯组长度 Lamp length/m 17.0 灯与灯间距 Distance between lamps/m 1.1 灯距水面高度 Height of lamp above water/m 3.0 船尾第一个灯与船尾距离 Fist lamp near stern/m 2.8 渔船集鱼灯总功率 Total power/kW 75 根据上述参数计算了船舷一侧垂直于船舯线不同距离的照度值，并与Choi等[13]的实测值进行对比，结果显示优化后的蒙特卡罗计算结果与实测值平均误差约为5.23%，计算值接近实测值 (图7)。作为对比，应用线光源法和点光源法计算结果与历史研究实测值的平均相对误差分别为52.62%和54.00%，并且随着距离增大，相对误差增加。
![]()
图 7 优化模型应用计算结果与实测值对比(差异系数0~0.23%)Figure 7. Calculated and measured values of horizontal illuminance at specific points beneath water (Difference coefficient 0−0.23%)3. 讨论
蒙特卡罗模拟通过大量随机抽样观察规律，样本量决定模拟结果的准确性。研究发现，光子数取值越大，模拟结果的差异系数越小且趋于极值。究其原因，当光子数取值足够大时，计算机所生成的伪随机数分布情况越接近概率分布函数，当然这还取决于伪随机数的生成方法[25,29]。本文通过预试验发现当光子数取值108时，模拟结果的差异系数为0.07%，此后差异系数趋于稳定值。本试验也发现光子数取值对模拟时长有较大影响。因此光束取值须兼顾模拟结果的稳定性、计算机性能等因素。本试验还发现，距离灯具越近的计算点，模拟的稳定性越高，这种现象在光子数取值较小时尤为明显，当光子数取值增大，远近距离计算点之间结果稳定性的差异减少。这是由于近光源处接受面所对应的灯具立体角较大，相同光子密度条件下，往往造成近光源处接受面能接受到更多光子，光子数越多，进一步造成计算结果越容易趋于稳定值。
针对光子辐射模型方面，本文利用灯具光度分布数据，从水平和垂直两个方向构建光子辐射概率模型并进行验证，发现优化模型的计算结果更接近事实。较官文江等[20]的研究，优化辐射模型充分考虑了灯具不同投影旋转角的光束辐射差异性，可以更加准确地计算灯与灯间光照叠加区域的亮度，进而计算合适的灯间安装距离以达到光源利用最大化。本文提出的光子辐射模型存在局限性，即假定不同的水平和垂直剖面内配光曲线分布趋势一致。因此，本文提出的建模方法无法应用于结构复杂的光源，建议选用二维随机变量分布函数描述[30]。
在算例中，本研究对比分析了几何算法与蒙特卡罗算法结果，发现线光源模型与点光源模型计算结果大于蒙特卡罗算法结果。原因是当海面具有一定斜率，光束在水汽交互界面的透射率发生改变。本文引用Cox和Munk[27]观测模型计算海面倾斜角，可拓展计算不同风速条件下的海面照度情况。
4. 结论
本文利用数值模拟方法研究了蒙特卡罗模拟过程中光子数取值对计算结果稳定性的影响，讨论了造成这一影响的原因。同时，根据光度分布数据提出了光子辐射优化模型并进行验证，发现新的光子辐射模型与实际情况较为接近，能够较为全面地表征灯具不同方向光子的辐射特征。在此基础上，结合实际渔船的灯光配置参数和风速，计算了光诱渔船实际作业环境下的光场分布，结果显示在海面波动条件下，蒙特卡罗计算结果比几何方法计算结果更接近实测值。
-
![]()
图 1 泸沽湖水声学探测航线及渔业捕捞点图
Figure 1. Detection route of hydroacoustic survey and fishing sites in Lugu Lake
![]()
图 3 泸沽湖鱼类水平分布及调查分区鱼类密度分布图
Figure 3. Horizontal distribution of fish in Lugu Lake and fish density in different zones
![]()
图 4 不同水生植物生长情况对鱼类分布的影响示意图
Figure 4. Sketch map of influence of different aquatic plants growth on fish distribution
表 1 各采样点鱼类组成及数量
Table 1 Species composition and number of fish at different sampling sites
种类
Species网具采样区 Area No. and fishing area 合计
Total1 普乐
Pule2 小渔坝
Xiaoyuba3 里格
Lige4 麻风
Mafeng5 大咀
Daju鲫 Carassius auratus 121 99 81 93 99 493 鲤 Cyprinus carpio 5 13 6 9 7 39 草鱼 Ctenopharyngodon idellus 2 8 2 6 0 18 大鳞副泥鳅 Paramisgurnus dabryanus 3 0 5 0 8 16 小口裂腹鱼 Schizothorax microstomus 0 2 0 0 0 2 总数 Amount 131 122 93 108 114 568 占比 Proportion/% 23.0 21.5 16.4 19.0 20.1 100 
下载: 导出CSV
表 2 渔获物组成及体长体质量结果
Table 2 Composition of catches and results of body length and body mass
种类
Species鱼类数量
Number of fish占比
Proportion/%体长
Length/mm体长均值
Average length/mm体质量
Body mass/kg体质量均值
Average body mass/kg鲫 Carassius auratus 493 86.8 101.0~256.0 154.0 0.035~0.305 0.082 鲤 Cyprinus carpio 39 6.9 234.0~501.0 395.0 0.210~3.435 1.367 草鱼 Ctenopharyngodon idellus 18 3.2 263.0~517.0 423.0 0.350~2.420 1.105 大鳞副泥鳅 Paramisgurnus dabryanus 16 2.8 167.0~202.0 172.0 0.182~0.248 0.194 小口裂腹鱼 Schizothorax microstomus 2 0.3 286.0~302.0 294.0 0.291~0.310 0.301
下载: 导出CSV
-
[1] 李英南, 赵晟, 王忠泽. 泸沽湖特有水生生物的保护初探[J]. 云南环境科学, 2000, 19(2): 39-41. [2] 彭徐, 徐大勇, 董艳珍, 等. 泸沽湖鱼类资源现状及保护对策[J]. 西昌学院学报(自然科学版), 2015, 29(2): 1-4. [3] 周杰. 宁蒗泸沽湖自然保护区的评价及其保护管理对策[J]. 林业调查规划, 2004, 29(S1): 89-92. [4] 万晔. 泸沽湖自然生态系统结构研究[J]. 地理与国土研究, 1998, 14(1): 51-54. [5] 孔德平, 陈小勇, 杨君兴. 泸沽湖鱼类区系现状及人为影响成因的初步探讨[J]. 动物学研究, 2006, 27(1): 94-97. [6] 孔德平, 秦涛, 范亦农, 等. 邛海鱼类资源与空间分布的水声学调查[J]. 水生态学杂志, 2019, 40(1): 24-31. [7] AGLEN A. Random errors of acoustic fish abundance estimates in relation to the survey grid density applied[M]. Bergen: Symposium on Fisheries Acoustics, 1983: 293-298.
[8] MISUND O A. Underwater acoustics in marine fisheries and fisheries research[J]. Rev Fish Biol Fish, 1997, 7(1): 1-34. doi: 10.1023/A:1018476523423
[9] 王崇瑞, 张辉, 杜浩, 等. 采用BioSonics DT-X超声波回声仪评估青海湖裸鲤资源量及其空间分布[J]. 淡水渔业, 2011, 41(3): 15-21. doi: 10.3969/j.issn.1000-6907.2011.03.003 [10] 莫伟均, 王从锋, 秦孝辉, 等. 北盘江董箐与光照库区鱼类资源水声学调查[J]. 水生态学杂志, 2015, 36(3): 10-17. [11] 孙明波, 谷孝鸿, 曾庆飞, 等. 基于水声学方法的太湖鱼类空间分布和资源量评估[J]. 湖泊科学, 2013, 25(1): 99-107. doi: 10.3969/j.issn.1003-5427.2013.01.013 [12] 张吉昌, 李显森, 赵宪勇, 等. 中东大西洋中部海域小型中上层鱼类集群形态与时空分布特征[J]. 海洋渔业, 2017, 39(6): 601-610. doi: 10.3969/j.issn.1004-2490.2017.06.001 [13] 凌建忠, 严利平, 李建生. 应用资源密度面积法估算东海区头足类资源量[J]. 海洋渔业, 2008(1): 56-61. doi: 10.3969/j.issn.1004-2490.2008.01.010 [14] 谭细畅, 史建全, 张宏, 等. EY60回声探测仪在青海湖鱼类资源量评估中的应用[J]. 湖泊科学, 2009, 21(6): 865-872. doi: 10.18307/2009.0617 [15] TAMEISHI H, SHINOMIYA H, AOKI I, et al. Understanding Japanese sardine migrations using acoustic and other aids[J]. ICES J Mar Sci, 1996, 53(2): 167-171. doi: 10.1006/jmsc.1996.0018
[16] GODLEWSKA M, WIERZOWSKI A. Hydroacoustical parameters of fish in reservoirs with contrasting levels of eutrophication[J]. Aquat Living Resour, 2003, 16(3): 167-173. doi: 10.1016/S0990-7440(03)00014-7
[17] 褚新洛, 陈银瑞. 云南鱼类志(上册)[M]. 北京: 科学出版社, 1989: 329-353. [18] 宋红桥, 王振华, 单建军, 等. 人为噪音对鱼类影响的研究进展[J]. 农学学报, 2017, 7(11): 46-50. [19] WHEELAND L J, ROSE G A. Quantifying fish avoidance of small acoustic survey vessels in boreal lakes and reservoirs[J]. Ecol Freshw Fish, 2015, 24(1): 67-76. doi: 10.1111/eff.12126
[20] 陈宜瑜, 张卫, 黄顺友. 泸沽湖裂腹鱼类的物种形成[J]. 动物学报, 1982(03): 16-24. [21] 董云仙, 谭志卫, 王俊松. 泸沽湖生态系统问题分析[J]. 环境科学导刊, 2011, 30(1): 39-43. doi: 10.3969/j.issn.1673-9655.2011.01.011 [22] 熊飞, 李文朝, 潘继征, 等. 云南抚仙湖鱼类资源现状与变化[J]. 湖泊科学, 2006, 18(3): 305-311. doi: 10.3321/j.issn:1003-5427.2006.03.018 [23] 费骥慧, 唐涛, 邵晓阳. 洱海渔业资源与渔业发展模式[J]. 湿地科学, 2011, 9(3): 277-283. [24] 陈自明, 杨君兴, 苏瑞凤, 等. 滇池土著鱼类现状[J]. 生物多样性, 2001, 9(4): 407-413. doi: 10.3321/j.issn:1005-0094.2001.04.013 [25] 陈文静, 贺刚, 吴斌, 等. 鄱阳湖通江水道鱼类空间分布特征及资源量评估[J]. 湖泊科学, 2017, 29(4): 923-931. [26] JING D, HAN J, WANG X, et al. A method to estimate the abundance of fish based on dual-frequency identification sonar (DIDSON) imaging[J]. Fish Sci, 2017, 83(5): 685-697. doi: 10.1007/s12562-017-1111-3
[27] MCQUINN I, SIMARD Y, STROUD T, et al. An adaptive, integrated "acoustic trawl" survey design for Atlantic cod (Gadus morhua) with estimation of the acoustic and trawl dead zones[J]. ICES J Mar Sci, 2005, 62: 93-106. doi: 10.1016/j.icesjms.2004.06.023
[28] NERO R W, THOMPSON C H, MICHAEL J J. In situ acoustic estimates of the swimbladder volume of Atlantic herring (Clupea harengus)[J]. J Acoust Soci Am, 2002, 112(3): 323-337.
[29] 孔燕, 余艳红, 苏斌. 云南九大高原湖泊流域现行管理体制及其完善建议[J]. 水生态学杂志, 2018, 39(3): 67-75. [30] 印杰, 田畛, 赵振山. 泥鳅食性的初步研究[J]. 水利渔业, 2000, 20(5): 15-16. doi: 10.3969/j.issn.1003-1278.2000.05.007 [31] 戚文华, 郭延蜀, 李雪芝, 等. 子陵栉鰕虎鱼繁殖特性、胚胎及仔鱼的发育[J]. 动物学杂志, 2008, 43(5): 13-24. doi: 10.3969/j.issn.0250-3263.2008.05.003 [32] 林植华, 雷焕宗, 陈利丽, 等. 棒花鱼形态特征的两性异形和雌性个体生育力[J]. 四川动物, 2007, 26(4): 910-913. doi: 10.3969/j.issn.1000-7083.2007.04.045 [33] 秦玉丽, 李林春, 黄荣静. 麦穗鱼的生物学特性及养殖技术[J]. 江苏农业科学, 2005(3): 114-116. doi: 10.3969/j.issn.1002-1302.2005.03.040 [34] 王权, 王建国, 封琦, 等. 中华鳑鲏的生物学特性及人工养殖技术[J]. 江苏农业科学, 2014, 42(5): 193-194. doi: 10.3969/j.issn.1002-1302.2014.05.064 [35] 连玉喜, 黄耿, GODLEWSKA M, 等. 基于水声学探测的香溪河鱼类资源时空分布特征评估[J]. 水生生物学报, 2015, 39(5): 920-929. doi: 10.7541/2015.121 -
期刊类型引用(2)
1. 邱星宇，刘庆霞，陈作志，蔡研聪，黄洪辉. 2023年春季南沙珊瑚岛礁主要鱼类碳氮稳定同位素研究. 热带海洋学报. 2024(06): 104-113 . 
百度学术
2. 彭谦，王啟芳，宋普庆，黄丁勇，张涵，王建佳，郑新庆. 秋季涠洲岛珊瑚礁主要鱼类营养关系的初步研究. 海洋学报. 2023(09): 91-104 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 4966
- HTML全文浏览量: 2585
- PDF下载量: 85
- 被引次数: 2
粤公网安备 44010502001741号
