李斌，湯勇，孫建富，李佩杰，胡亞東，張劍波

（1.大連海洋大學(xué)海洋科技與環(huán)境學(xué)院，遼寧大連116023；2.中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所，廣東廣州510300；3.大連海洋大學(xué)經(jīng)濟管理學(xué)院，遼寧大連116023；4.河南省水產(chǎn)科學(xué)研究院，河南鄭州450044）

基于聲學(xué)方法的黃河三門峽水庫漁業(yè)資源空間分布研究

李斌1、2，湯勇1，孫建富3，李佩杰1，胡亞東4，張劍波4

（1.大連海洋大學(xué)海洋科技與環(huán)境學(xué)院，遼寧大連116023；2.中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所，廣東廣州510300；3.大連海洋大學(xué)經(jīng)濟管理學(xué)院，遼寧大連116023；4.河南省水產(chǎn)科學(xué)研究院，河南鄭州450044）

為探討聲學(xué)方法在黃河流域生態(tài)環(huán)境研究中的適用性，于2013年11月23日使用分裂式波束科學(xué)魚探儀（Simrad EY60，70 kHz，挪威）對黃河三門峽水庫西灣至庫尾大壩間水域進行了漁業(yè)資源聲學(xué)調(diào)查，并通過現(xiàn)場魚類目標(biāo)強度測定和回波積分方法對庫區(qū)內(nèi)漁業(yè)資源密度、資源現(xiàn)存量和空間分布進行了初步探查與估算。結(jié)果表明：此次調(diào)查庫區(qū)內(nèi)聲學(xué)積分值（Nautical area scattering coefficient，NASC）為5.59 m2/nmi2，漁業(yè)資源密度平均為0.002 ind./m2，資源現(xiàn)存量為42.36 t；在魚類資源密度空間分布統(tǒng)計方面，將調(diào)查區(qū)域庫區(qū)分為上、中、下游3個調(diào)查區(qū)域（A、B、C區(qū)），8個積分區(qū)段，依次為A-Ⅰ、A-Ⅱ、A-Ⅲ，B-Ⅰ、B-Ⅱ，C-Ⅰ、C-Ⅱ、C-Ⅲ區(qū)段，其對應(yīng)的魚類資源密度分別為0.002、0.017、0.013 ind./m2，0.004、0.008 ind./m2，0.013、0.016、0.037 ind./m2，魚類資源主要集中于上游的A-Ⅱ、A-Ⅲ區(qū)段及水庫下游區(qū)域；對應(yīng)魚類的垂直分布，A-Ⅱ、A-Ⅲ、C-Ⅰ、C-Ⅱ區(qū)段魚類主要分布于5 m以淺水層，A-Ⅰ、B-Ⅱ區(qū)段魚類主要分布于10 m以淺水層，而B-Ⅰ區(qū)段魚類均分布于5 m以深水層且數(shù)量較少，C-Ⅲ區(qū)段魚類則主要分布于10～20 m水層；庫區(qū)內(nèi)魚類目標(biāo)強度（Target strength，TS）多分布于-62.5～-47.5 dB之間，且目標(biāo)強度大于-44.5 dB、體長大于25 cm的個體均分布于5 m以深水層。研究表明，聲學(xué)方法在淡水水庫漁業(yè)資源評估中體現(xiàn)出良好的適用性，本研究推動了該方法在黃河流域水域生態(tài)環(huán)境研究中的應(yīng)用，為今后用聲學(xué)方法評估內(nèi)陸河流、湖泊和水庫漁業(yè)資源提供了理論基礎(chǔ)和科學(xué)依據(jù)。

漁業(yè)資源；聲學(xué)評估；空間分布；三門峽水庫

三門峽水利樞紐位于河南省三門峽市的東北部，北鄰山西省，西鄰陜西省，于20世紀(jì)50年代開工興建，是中國在黃河上興建的第一座以防洪為主，兼顧防凌、供水、灌溉、發(fā)電等綜合利用的大型水利工程［1］。三門峽庫區(qū)包括黃河干流龍門禹門口至三門峽大壩河段，全長238.4 km，總面積約84 290 hm2［2］。2011年三門峽庫區(qū)全年入庫水量為260億m3，出庫水量為259.1億m3，桃汛期最大下泄流量達到2660 m3/s［3］。每年10月至翌年6月為水庫的蓄水期，在此期間，沿河淺灘、沼澤等天然濕地為各類水生動植物提供了必需的棲息、索餌和繁殖場所［2］。

自2002年起，為擴大黃河下游主河槽過洪能力，減少小浪底水庫淤積，每年6月中下旬至7月中上旬，黃河水利委員會利用三門峽、萬家寨等水庫與下游的小浪底水利樞紐進行聯(lián)合調(diào)度，實施黃河調(diào)水調(diào)沙［4］。雖然每年一至數(shù)次的調(diào)水調(diào)沙能夠有效沖刷水庫泥沙及下游河槽，對治理黃河起到較為積極的作用，但同時也對黃河流域的魚類種質(zhì)資源、浮游生物和漁業(yè)生態(tài)環(huán)境產(chǎn)生了較大影響，使黃河流域的生態(tài)環(huán)境現(xiàn)狀變動把握困難［5-7］。如何在保護黃河流域生態(tài)資源環(huán)境的前提下發(fā)展水利水電和淡水漁業(yè)，是黃河流域經(jīng)濟發(fā)展的一個重點。目前，三門峽水庫每年11月下旬進行增殖放流，主要放流魚類包括黃河鯉Cyprinus carpio、鰱Hypophthalmichthys molitrix、鳙Aristichthys nobilis、草魚Ctenopharyngodon idella等，苗種規(guī)格大于50 g，每年放流數(shù)量不少于300萬尾，旨在對黃河流域的種質(zhì)資源進行保護，緩解庫區(qū)內(nèi)因調(diào)水調(diào)沙、發(fā)展水利水電引起的生態(tài)結(jié)構(gòu)改變。

由于淡水水域多采用人工增殖放養(yǎng)和定置網(wǎng)具大規(guī)格體長魚類選擇性捕撈，傳統(tǒng)的漁業(yè)資源調(diào)查評估方法較難實施［8］。近年來，漁業(yè)資源聲學(xué)評估方法以其高效性、非接觸性、不損害漁業(yè)資源、回波映像數(shù)據(jù)可保存和回放及對比分析等優(yōu)點，已成為歐美漁業(yè)發(fā)達國家進行水中游泳生物資源現(xiàn)存量調(diào)查的主要方法，是制定海域不同魚類總允許漁獲量（Total allowable catch，TAC）的重要參考依據(jù)，并在漁業(yè)資源評估及生態(tài)環(huán)境調(diào)查等領(lǐng)域得到了廣泛應(yīng)用［9-10］。隨著漁業(yè)聲學(xué)測量技術(shù)在國內(nèi)的推廣，漁業(yè)資源聲學(xué)評估方法在淡水水域（江河、湖泊和水庫等）的應(yīng)用也逐漸增多，均取得了較好的效果［9，11-17］。

本研究中，利用聲學(xué)方法對增殖放流后的黃河三門峽水庫西灣至庫尾大壩間水域的漁業(yè)資源進行探查和評估，對魚類的回波特征及其分布、單體目標(biāo)的大小及其空間結(jié)構(gòu)分布、資源現(xiàn)存量和資源密度空間分布等進行統(tǒng)計分析，旨在探討聲學(xué)方法在黃河流域資源調(diào)查中的適用性，以期為制定增殖放流計劃、改善捕撈方法提供理論基礎(chǔ)，為合理治理黃河、保護黃河流域生態(tài)資源環(huán)境提供科學(xué)依據(jù)。

1 材料與方法

1.1 調(diào)查時間、水域

聲學(xué)調(diào)查時間為2013年三門峽水庫魚類增殖放流后的11月23日，調(diào)查區(qū)域為三門峽水庫西灣至庫尾大壩間（111°8′32.92″～111°20′16.63″E，34°47′23.28″～34°49′59.7″N）約22 km長的水域。

1.2 方法

1.2.1 調(diào)查方法 調(diào)查航線設(shè)計參考中國海洋生物調(diào)查規(guī)范中的《漁業(yè)資源聲學(xué)調(diào)查與評估》［18］，并結(jié)合庫區(qū)狹長的地理地貌特征，將調(diào)查航線設(shè)計為“之”字形。然而受實際駕駛航行條件限制，同時根據(jù)現(xiàn)場情況繞行定置網(wǎng)具，調(diào)查過程中所走實際航線與設(shè)計航線會出現(xiàn)偏差，如圖1所示。

使用分裂式波束科學(xué)魚探儀（Simrad EY60， 70 kHz，挪威）進行聲學(xué)調(diào)查，使用專用軟件ER60進行聲學(xué)數(shù)據(jù)的采集與收錄，動態(tài)經(jīng)緯度位置信息由GPS（Gamin GPS60CSx，美國）獲得，科學(xué)魚探儀主要技術(shù)參數(shù)設(shè)定見表1。調(diào)查前使用直徑為38.1 mm的鎢鋼合金標(biāo)準(zhǔn)球（Tungsten carbide standard sphere），按照國際通用的標(biāo)準(zhǔn)目標(biāo)法對科學(xué)魚探儀的系統(tǒng)收發(fā)增益系數(shù)進行現(xiàn)場校正［19］。聲學(xué)調(diào)查使用漁政快艇（長16 m）進行走航測量，換能器放置于導(dǎo)流罩內(nèi)，利用鋼管懸掛固定于船體左舷外側(cè)，吃水1 m，調(diào)查航速為5～7 kn。

圖1 三門峽水庫聲學(xué)調(diào)查航線及探測區(qū)域劃分Fig.1 Acoustic survey track and the partitions of detecting areas in Sanmenxia Reservoir

表1 EY60科學(xué)魚探儀主要技術(shù)參數(shù)設(shè)定Tab.1 M ain technical parameters in scientific echosounder EY60

1.2.2 魚類資源聲學(xué)評估 使用Echoview 4.6（Myriax，澳大利亞）軟件進行聲學(xué)數(shù)據(jù)后處理及魚類資源評估，參考中國海洋生物調(diào)查規(guī)范聲學(xué)選項中的回波積分法對魚類資源密度及其空間分布特征進行統(tǒng)計分析［18］。針對此次調(diào)查魚類回波特征分布、地理位置特征、水深分布和航行時間，將整個調(diào)查區(qū)域分為3個區(qū)分別進行資源密度和資源量的統(tǒng)計分析，其中，西灣至?xí)d鎮(zhèn)為A區(qū)，會興鎮(zhèn)至棗溝為B區(qū)，棗溝至庫尾大壩為C區(qū)。在魚類資源密度空間分布特征分析方面，對所分區(qū)域按照航行距離進行積分區(qū)段劃分，其中A區(qū)分為A-Ⅰ、A-Ⅱ、A-Ⅲ3個區(qū)段，B區(qū)分為B-Ⅰ、B-Ⅱ兩個區(qū)段，C區(qū)分為C-Ⅰ、C-Ⅱ、C-Ⅲ3個區(qū)段，分別對各個區(qū)段的資源平均密度和以深度方向5 m間隔從表層至庫底不同水層的資源密度進行統(tǒng)計分析。

在對回波映像數(shù)據(jù)進行分析處理的過程中，對不同類型的回波映像進行分類，剔除庫區(qū)表層航行氣泡回波、浮游生物回波、庫區(qū)底部因水庫蓄水淹沒的植被回波和某些干涉信號回波，重新設(shè)置庫區(qū)積分起始水層和積分終止水層，保留用于估算資源密度及資源量的單體魚和魚群回波映像數(shù)據(jù)。

利用Echoview 4.6軟件計算輸出不同調(diào)查區(qū)域的聲學(xué)積分值（Nautical area scattering coefficient，NASC，m2/nmi2）、平均目標(biāo)強度（Target strength，TS，dB）、平均面積后向散射系數(shù)（Area backscattering coefficient，sa）和后向散射截面（Backscattering cross-section，σbs，m2/ind.）等，其計算公式為

在不考慮體長分布的情況下，可獲得各調(diào)查區(qū)域的平均資源密度（ind./m2）為

而在回波積分處理單元內(nèi)，不同體長階段魚類對應(yīng)的平均資源密度為

其中：p為不同目標(biāo)強度階層的頻度分布；i為魚類不同體長階級。

由于本次調(diào)查未對應(yīng)聲學(xué)回波進行網(wǎng)具采樣，且大多數(shù)魚類呈單體形式分布，故采用現(xiàn)場測定法統(tǒng)計魚類目標(biāo)強度頻度分布。利用Echoview 4.6軟件對所獲得的聲學(xué)數(shù)據(jù)進行單體目標(biāo)檢測和單體目標(biāo)軌跡追蹤后，輸出不同目標(biāo)強度階層魚類單體數(shù)量，得到對應(yīng)的平均目標(biāo)強度頻度分布。

同時，因無法對庫區(qū)內(nèi)的魚類進行魚種判別并確定其體長分布，且?guī)靺^(qū)內(nèi)魚類多為管鰾魚類，因此，目標(biāo)強度（TS）與魚類體長（L，cm）的關(guān)系統(tǒng)一使用管鰾魚類模型［18］：

體質(zhì)量（W，g）的計算參考黃河鯽Carassius auratus的體長-體質(zhì)量擬合模型［20］：

庫區(qū)內(nèi)漁業(yè)資源現(xiàn)存量為

其中：Wi為不同體長階級魚類的體質(zhì)量；Aj為不同調(diào)查區(qū)域的面積（km2），利用谷歌軟件進行測算，A、B、C 3個調(diào)查區(qū)域的面積依次為6.46、3.61、8.60 km2。

2 結(jié)果與分析

2.1 各調(diào)查區(qū)域的回波特征

利用Echoview 4.6軟件對庫區(qū)內(nèi)各調(diào)查區(qū)域回波特征進行統(tǒng)計分析。整體上庫區(qū)內(nèi)魚類呈分散的單體形式存在，偶爾會有小型魚群回波出現(xiàn)，庫區(qū)內(nèi)不同調(diào)查區(qū)域的魚類回波特征如圖2所示。在A區(qū)，5 m以淺水層存在回波強度較小且極為分散的單體回波，而在5 m以深水層，存在較多且較為分散的單體回波；在B區(qū)，5 m以淺水層極少出現(xiàn)魚類回波，而在5 m以深水層的近庫區(qū)底部，存在較少且極為分散的單體回波；在C區(qū)，5 m以淺水層魚類分布較少，但存在小型魚群回波且數(shù)量較少，而在5 m以深水層，存在大量強度較大且較為集中的單體回波，同時存在魚群回波（圖2）。

2.2 單體魚類目標(biāo)強度分布

利用Echoview 4.6軟件進行單體目標(biāo)檢測和單體目標(biāo)軌跡追蹤，輸出不同目標(biāo)強度階層魚類單體數(shù)量，對應(yīng)獲得各調(diào)查區(qū)域及全航段的平均目標(biāo)強度頻度分布，結(jié)果如圖3所示。

由圖3可知，此次調(diào)查庫區(qū)內(nèi)魚類目標(biāo)強度大多分布于-62.5～-47.5 dB之間，通過目標(biāo)強度與魚類體長模型進行換算，對應(yīng)3.16～17.78 cm的個體，整體上以小個體魚類分布為主。從西灣至庫尾大壩，A區(qū)目標(biāo)強度大多分布于-62.5～-47.5 dB之間，對應(yīng)3.16～17.78 cm的個體，同庫區(qū)內(nèi)整體分布相吻合；B區(qū)目標(biāo)強度則多分布于-62.5～-53.5 dB之間，對應(yīng)3.16～8.91 cm的個體，較小個體魚類大多分布于該區(qū)域；C區(qū)目標(biāo)強度多分布于-64～-61 dB、-56.5～-47.5 dB之間，對應(yīng)2.66～3.76、6.31～17.78 cm的個體。從西灣至庫尾大壩，目標(biāo)強度分布呈逐漸增大的趨勢，對應(yīng)單體魚類體長分布逐漸增大，但仍以小個體魚類分布占據(jù)主要部分。

圖2 三門峽水庫各調(diào)查區(qū)域的回波特征Fig.2 Features of echo in different areas in Sanmenxia Reservoir

圖3 三門峽水庫不同調(diào)查區(qū)域目標(biāo)強度的頻度分布Fig.3 Frequency profile of target strength in different areas in Sanmenxia Reservoir

根據(jù)庫區(qū)內(nèi)目標(biāo)強度分布以小個體魚類分布為主的特征，且增殖放流多為中下層魚類，為了考察增殖放流效果，結(jié)合三門峽水庫的水深分布，以深度方向5m水深為界限，分別統(tǒng)計5m以淺水層及5 m以深水層目標(biāo)強度的頻度分布，結(jié)果如圖4所示。

由圖4可知，從西灣至庫尾大壩，5 m以淺水層目標(biāo)強度分布均在-50.5 dB以下，而5 m以深水層目標(biāo)強度則分布于-65.5～-29.5 dB之間。其中5 m以淺水層，A區(qū)目標(biāo)強度主要分布于-62.5～-59.5、-53.5～-47.5 dB之間，而B區(qū)目標(biāo)強度主要分布于-65.5～-50.5 dB之間，C區(qū)目標(biāo)強度則主要分布于-65.5～-62.5、-52～-49 dB之間；5 m以深水層，A區(qū)目標(biāo)強度主要分布于-62.5～-47.5 dB間，而B區(qū)目標(biāo)強度主要分布于-64～-61、-59.5～-53.5 dB之間，C區(qū)目標(biāo)強度則主要分布于-64～-61、-56.5～-47.5 dB之間。通過目標(biāo)強度與魚類體長換算模型進行計算，庫區(qū)內(nèi)5 m以淺水層魚類體長均在12.59 cm以下，而目標(biāo)強度大于-44.5 dB、體長大于25 cm的個體均分布于5 m以深水層。

圖4 三門峽水庫不同調(diào)查區(qū)域分水層目標(biāo)強度的頻度分布Fig.4 Layered frequency profile of target strength in different areas in Sanmenxia Reservoir

2.3 資源密度及資源量統(tǒng)計

此次調(diào)查，庫區(qū)內(nèi)聲學(xué)積分值為5.59 m2/nmi2，從西灣至庫尾大壩，A、B、C 3個調(diào)查區(qū)域內(nèi)的聲學(xué)積分值依次為8.25、5.77、8.64 m2/nmi2。在回波積分處理單元內(nèi)，根據(jù)公式（1）、（2）、（4）可獲得各調(diào)查區(qū)域不同目標(biāo)強度階層魚類對應(yīng)的平均資源密度，結(jié)果如表2所示。

表2 各調(diào)查區(qū)域不同目標(biāo)強度階層魚類的平均資源密度Tab.2 Average fisheries resources density of different target strength hierarchies in various surveyed areas ind./m2

在獲得各調(diào)查區(qū)域不同目標(biāo)強度階層對應(yīng)的魚類平均資源密度后，根據(jù)公式（5）、（6）、（7）可計算各調(diào)查區(qū)域漁業(yè)資源量，結(jié)果如表3所示。

表3 各調(diào)查區(qū)域漁業(yè)資源量Tab.3 Fisheries resources in various surveyed areas

2.4 平均資源密度空間分布

2.4.1 水平空間分布 此次調(diào)查庫區(qū)內(nèi)平均資源密度為0.002 ind./m2，從西灣至庫尾大壩，A、B、C 3個調(diào)查區(qū)域平均資源密度分別為0.002、0.004、0.021 ind./m2。為了更為細致地考察資源密度的分布情況，沿西灣至庫尾大壩航線方向，將3個調(diào)查區(qū)域劃分為8個積分區(qū)段，依次為A-Ⅰ、A-Ⅱ、A-Ⅲ區(qū)段，B-Ⅰ、B-Ⅱ區(qū)段，C-Ⅰ、C-Ⅱ、C-Ⅲ區(qū)段，對應(yīng)的平均資源密度分別為0.002、0.017、0.013 ind./m2，0.004、0.008 ind./m2，0.013、0.016、0.037 ind./m2，如圖5所示?？梢钥闯觯瑤靺^(qū)內(nèi)魚類分布以庫尾大壩附近水域（C-Ⅲ區(qū)段）最為密集，A-Ⅱ和C-Ⅱ區(qū)段次之，且兩者差別不大，A-Ⅰ和B-Ⅰ區(qū)段魚類分布則為最少。

圖5 三門峽水庫平均資源密度的水平空間分布Fig.5 Horizontal and spatial distribution of average resources density in Sanmenxia Reservoir

2.4.2 垂直空間分布 對不同水層的魚類回波進行積分，得到三門峽水庫平均資源密度垂直空間分布，如圖6所示。沿西灣至庫尾大壩航線方向，A-Ⅰ區(qū)段魚類主要分布于10 m以淺水層，為0.014～0.020 ind./m2，而A-Ⅱ、A-Ⅲ區(qū)段魚類則主要分布于5 m以淺水層，分別為0.012、0.017 ind./m2。B-Ⅰ區(qū)段魚類主要分布在5 m以深水層且資源密度相差不大，約為0.001 ind./m2，而B-Ⅱ區(qū)段魚類則主要分布于10 m以淺水層，為0.005～0.006 ind./m2。C-Ⅰ區(qū)段魚類主要分布于5 m以淺水層，為0.013 ind./m2，10～15 m水層次之，為0.007 ind./m2，而5～10、15～20 m水層則相差不大，約為0.002 ind./m2；C-Ⅱ區(qū)段魚類集中分布于5 m以淺水層，為0.124 ind./m2，其他水層僅有少量分布；C-Ⅲ區(qū)段魚類主要分布于10～20m水層，為0.008 ind./m2，10m以淺水層次之，為0.004 ind./m2，20 m以深水層最小，為0.002 ind./m2。

圖6 三門峽水庫平均資源密度的垂直空間分布Fig.6 Vertical and spatial distribution of average resources density in Sanmenxia Reservoir

3 討論

3.1 聲學(xué)調(diào)查方法的適用性

針對三門峽庫區(qū)狹長的地理特征，本次聲學(xué)調(diào)查中采用“之”字型走航方式，以規(guī)避庫區(qū)內(nèi)設(shè)置的定置網(wǎng)具。設(shè)計調(diào)查航線覆蓋系數(shù)為6.99，實際調(diào)查航線覆蓋系數(shù)為5.93，可以滿足空間采樣隨機性的聲學(xué)采樣要求，能夠較好地反映魚類在庫區(qū)內(nèi)的空間分布［21］。但由于使用的調(diào)查船只較小，航行穩(wěn)定性較差，同時需繞行定置網(wǎng)具，故實際航線與設(shè)計航線出現(xiàn)了一定的偏差?？紤]到庫區(qū)內(nèi)魚類回波整體上呈分散的單體形式存在，采用了分區(qū)域計算平均聲學(xué)積分值的方法，以降低航線偏差所產(chǎn)生的空間采樣密度不均勻的影響，故認為航線偏差對資源評估結(jié)果影響不大。

在庫區(qū)漁業(yè)資源現(xiàn)存量估算統(tǒng)計方面，《漁業(yè)資源聲學(xué)調(diào)查與評估》［18］規(guī)范中要求使用選擇性較低的專用調(diào)查網(wǎng)具（包括底層拖網(wǎng)和變水層拖網(wǎng)），對產(chǎn)生聲學(xué)回波的生物進行采樣，用以進行魚種識別并確定魚類體長分布。由于本次調(diào)查船只未對產(chǎn)生聲學(xué)回波的生物進行采樣，同時受限于不同魚種和體長的目標(biāo)強度關(guān)系研究尚未開展，魚類體質(zhì)量-體長擬合模型未根據(jù)具體魚種進行測定，捕撈數(shù)據(jù)統(tǒng)計不夠充分等因素影響，故資源量統(tǒng)計的準(zhǔn)確度存在偏差。此外，在調(diào)查過程中科學(xué)魚探儀換能器部分安裝在船體左舷外側(cè)，無法探測到表層1 m以內(nèi)的魚類，而且為排除庫區(qū)底部的反射信號，聲學(xué)評估水深終止于庫底之上0.5 m，形成聲學(xué)探測盲區(qū)，盲區(qū)內(nèi)的魚類視為庫底，不在積分范圍之內(nèi)。因此，資源量的評估結(jié)果與實際資源量存在一定偏差，但對于分析魚類資源平均密度及其空間結(jié)構(gòu)分布影響不大，說明應(yīng)用聲學(xué)方法對該水域進行漁業(yè)資源評估能夠取得較好的效果。

3.2 魚類目標(biāo)強度

本研究中，利用分裂式波束科學(xué)魚探儀進行現(xiàn)場測量并給出庫區(qū)內(nèi)魚類目標(biāo)強度的頻度分布，而現(xiàn)場測量魚類目標(biāo)強度的準(zhǔn)確性會受魚群密度影響，但Sawada等［22］的研究表明，單位采樣體積內(nèi)平均魚尾數(shù)（Nv）在不超過0.04 ind./m3時目標(biāo)強度現(xiàn)場測量有效。本次調(diào)查中，三門峽庫區(qū)Nv值為0.000 05 ind./m3，A、B、C 3個調(diào)查區(qū)域的Nv值分別為0.000 04、0.000 13、0.000 75 ind./m3，說明現(xiàn)場測量目標(biāo)強度方法有效，由此得出的魚類資源密度及其空間分布能夠較為準(zhǔn)確地反映三門峽水庫庫區(qū)內(nèi)魚類分布現(xiàn)狀。

目前，淡水魚類目標(biāo)強度測量的研究相對較少［23-24］，一般參考相近魚類或使用理論模型進行目標(biāo)強度計算，本研究中統(tǒng)一使用一般管鰾魚類目標(biāo)強度與體長關(guān)系來估算庫區(qū)內(nèi)漁業(yè)資源現(xiàn)存量，影響漁業(yè)資源聲學(xué)評估的精度。而對于將來能夠?qū)崿F(xiàn)精確到不同魚種和體長組成的資源量評估，需在今后的研究工作中進行完善。

3.3 調(diào)水調(diào)沙對庫區(qū)漁業(yè)資源的影響及應(yīng)對措施

三門峽水庫的建成為發(fā)展?jié)O業(yè)提供了良好的場所，且發(fā)展?jié)O業(yè)可凈化庫區(qū)水域環(huán)境、維護生態(tài)平衡，對保護庫區(qū)水域生態(tài)安全具有一定的支撐作用。然而每年一至數(shù)次的黃河調(diào)水調(diào)沙對三門峽水庫的漁業(yè)資源會產(chǎn)生較大影響，因此，研究調(diào)水調(diào)沙對三門峽水庫漁業(yè)資源及其生態(tài)環(huán)境的影響具有較高的理論意義和實際價值。

調(diào)水調(diào)沙可導(dǎo)致庫區(qū)內(nèi)魚類大批死亡，大量魚類隨泥沙被沖往下游，產(chǎn)生“流魚”現(xiàn)象，同時庫區(qū)內(nèi)漁業(yè)生態(tài)環(huán)境遭到破壞，干擾水生動植物的自然種群結(jié)構(gòu)、群體補充和繁殖生長周期，使庫區(qū)內(nèi)食物鏈、生態(tài)鏈關(guān)系發(fā)生變化［5，25］。針對以上調(diào)水調(diào)沙產(chǎn)生的問題，應(yīng)積極探索適應(yīng)調(diào)水調(diào)沙的漁業(yè)發(fā)展模式。包括開展可在當(dāng)年捕撈魚類品種的增殖放流，大力發(fā)展放牧式網(wǎng)箱養(yǎng)殖和沿岸休閑漁業(yè)，提高漁業(yè)綜合效益，降低調(diào)水調(diào)沙的影響，促進庫區(qū)漁業(yè)持續(xù)健康發(fā)展。也可嘗試利用聲學(xué)標(biāo)志對庫區(qū)內(nèi)主要經(jīng)濟魚種調(diào)水調(diào)沙期間的行為學(xué)特性進行追蹤，推斷調(diào)水調(diào)沙期間環(huán)境相對適宜魚類的避難所［26-27］。

3.4 回波計數(shù)和回波積分法估算資源密度的比較

針對三門峽庫區(qū)內(nèi)魚類整體低密度分布且大多數(shù)魚類呈單體形式分布的特點，為了更精確地估算庫區(qū)內(nèi)魚類資源密度的空間分布，嘗試利用回波計數(shù)和回波積分法分別估算了魚類的資源體積密度。結(jié)果表明：使用回波計數(shù)法估算整個庫區(qū)內(nèi)魚類資源密度為0.000 23 ind./m3，A、B、C 3個調(diào)查區(qū)域資源密度依次為0.000 15、0.000 16、0.000 72 ind./m3；利用回波積分法估算整個庫區(qū)魚類資源密度為0.000 31 ind./m3，A、B、C 3個調(diào)查區(qū)域資源密度依次為0.000 18、0.000 35、0.001 73 ind./m3。使用回波積分法估算整個庫區(qū)內(nèi)資源密度為用回波計數(shù)法估算資源密度的2.07倍，A、B、C 3個調(diào)查區(qū)域使用回波積分法估算的資源密度依次為用回波計數(shù)法估算資源密度的1.20、2.19和2.40倍?？梢钥闯觯珹區(qū)使用回波計數(shù)和回波積分法估算的資源密度的差別不大，而B區(qū)相差較大，可能是由于該區(qū)域內(nèi)魚類多呈回波強度較小的單體形式存在，軟件計數(shù)設(shè)置的目標(biāo)強度最小閾值限制了目標(biāo)強度較小單體回波的計數(shù)，C區(qū)相差也較大，可能與C區(qū)存在部分密集魚群回波有關(guān)。

3.5 研究不足及今后展望

本研究中存在以下不足之處，首先，聲學(xué)調(diào)查僅進行了一次且缺乏以前三門峽水庫魚類調(diào)查的相關(guān)基礎(chǔ)資料，對于不同季節(jié)魚類的空間分布特征尚不明確，因此，此次調(diào)查僅是對三門峽水庫漁業(yè)資源及其空間分布進行初步探究。今后可在此次調(diào)查的基礎(chǔ)上，對不同季節(jié)魚類隨水溫、溶解氧等水質(zhì)參數(shù)變化，其資源量及其資源密度空間結(jié)構(gòu)分布的變化進行分析，以便更為詳盡地了解庫區(qū)內(nèi)魚類分布特征。其次，本次調(diào)查均在日間進行，且?guī)靺^(qū)內(nèi)部分魚類具有晝夜垂直遷移特性［17，28］，因此，資源聲學(xué)評估存在晝夜誤差。最后，聲學(xué)調(diào)查采用航行斷面生物回波采樣方式進行，少數(shù)庫區(qū)因部分因素未能進行采樣，可能會使資源評估結(jié)果產(chǎn)生一定誤差。

綜上所述，聲學(xué)方法適用于三門峽水庫漁業(yè)資源現(xiàn)存量的調(diào)查與評估，在今后的研究中，可以選擇在調(diào)水調(diào)沙及增殖放流前后的不同季節(jié)對三門峽水庫庫區(qū)內(nèi)漁業(yè)資源現(xiàn)存量及其空間分布進行調(diào)查，對比分析庫區(qū)內(nèi)漁業(yè)資源的生態(tài)變動及調(diào)水調(diào)沙對庫區(qū)漁業(yè)資源的影響，同時，也可以結(jié)合漁獲數(shù)據(jù)將增殖放流對庫區(qū)內(nèi)漁業(yè)資源的恢復(fù)作用進行效果評價，以期為制定增殖放流和捕撈計劃提供科學(xué)依據(jù)，為探索適應(yīng)調(diào)水調(diào)沙的漁業(yè)發(fā)展模式提供新的技術(shù)途徑。

［1］ 王勇澤，李誠，孫樹青，等.黃河三門峽段水環(huán)境健康風(fēng)險評價［J］.水資源保護，2007，23（1）：28-31.

［2］ 陳媛媛.黃河干流禹門口至三門峽段魚類物種多樣性及保護對策［J］.安徽農(nóng)業(yè)科學(xué)，2013，41（15）：6721-6724，6727.

［3］ 劉社強.2011年三門峽庫區(qū)水沙特點及沖淤分析［J］.人民黃河，2014，36（1）：24-26.

［4］ 介子林，胡亞東，張坤，等.減緩小浪底水庫調(diào)水調(diào)沙對魚類資源影響措施研究［J］.河南水產(chǎn)，2012（1）：31-32.［5］ 張軍燕，張建軍，沈紅保，等.小浪底水庫調(diào)水調(diào)沙對壺口至三門峽段魚類群落結(jié)構(gòu)的影響［J］.生態(tài)學(xué)雜志，2012，31（10）：2613-2618.

［6］ 朱國清，趙瑞亮，胡振平，等.小浪底水庫調(diào)水調(diào)沙對黃河中游魚類及生態(tài)敏感區(qū)的影響［J］.水生態(tài)學(xué)雜志，2012，33（5）：7-12.

［7］ 余斌，張建軍，劉一歐，等.小浪底水庫調(diào)水調(diào)沙對黃河壺口至三門峽段浮游生物資源的影響［J］.河北漁業(yè)，2013（1）：15-20.

［8］ 張贊，藺丹清，湯勇，等.基于聲學(xué)測量方法的大伙房水庫魚類資源季節(jié)變動特征［J］.南方水產(chǎn)科學(xué)，2014，10（6）：12-19.

［9］ 李斌，李佩杰，湯勇，等.黃河小浪底水庫主河道水域漁業(yè)資源聲學(xué)評估［J］.水產(chǎn)學(xué)報，2015，39（8）：1134-1143.

［10］ Koslow J A.The role of acoustics in ecosystem-based fishery management［J］.ICES Journal of Marine Science，2009，66（6）：966-973.

［11］ 陶江平，龔昱田，譚細暢，等.長江葛洲壩壩下江段魚類群落變化的時空特征［J］.中國科學(xué)：生命科學(xué)，2012，42（8）：677-688.

［12］ 譚細暢，李新輝，林建志，等.基于水聲學(xué)探測的兩個廣東魴產(chǎn)卵群體繁殖生態(tài)的差異性［J］.生態(tài)學(xué)報，2009，29（4）：1756-1762.

［13］ 張慧杰，楊德國，危起偉，等.葛洲壩至古老背江段魚類的水聲學(xué)調(diào)查［J］.長江流域資源與環(huán)境，2007，16（1）：86-91.

［14］ 譚細暢，史建全，張宏，等.EY60回聲探測儀在青海湖魚類資源量評估中的應(yīng)用［J］.湖泊科學(xué)，2009，21（6）：865-872.

［15］ 譚細暢，夏立啟，立川賢一，等.東湖放養(yǎng)魚類時空分布的水聲學(xué)研究［J］.水生生物學(xué)報，2002，26（6）：585-590.

［16］ 陶江平，陳永柏，喬曄，等.三峽水庫成庫期間魚類空間分布的水聲學(xué)研究［J］.水生態(tài)學(xué)雜志，2008，1（1）：25-33.

［17］ 王靖，張超，王丹，等.清河水庫鰱鳙魚類資源聲學(xué)評估——回波計數(shù)與回波積分法的比較［J］.南方水產(chǎn)，2010，6（5）：50-55.

［18］ 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局.GB/T 12763. 6-2007海洋調(diào)查規(guī)范第6部分：海洋生物調(diào)查［S］.北京：中國標(biāo)準(zhǔn)出版社，2008：90-97.

［19］ Simrad.Simrad ER60，Scientific Echo Sounder Software，Reference Manual，Release2.2.0［M］.Norway：Simrad Maritime AS Kongsberg，2008.

［20］ 邱小琮，王遠吉.黃河鯽魚鱗徑與體長、體長與體重相關(guān)關(guān)系的微機計算分析［J］.寧夏大學(xué)學(xué)報：自然科學(xué)版，2000，21（2）：186-188.

［21］ Parker-Stetter S L，Rudstam L G，Sullivan P J，et al.Standard Operating Procedures for Fisheries Acoustic Surveys in the Great Lakes［M］.Michigan：Special Publication，2009：53-62.

［22］ Sawada K，F(xiàn)urusawa M，Williamson N J.Conditions for the precisemeasurement of fish target strength in situ［J］.The Journal of the Marine Acoustics Society of Japan，1993，20（2）：73-79.

［23］ Warner D M，Rudstam L G，Klumb R A.In situ target strength of alewives in freshwater［J］.Transactions of the American Fisheries Society，2002，131（2）：212-223.

［24］ Rudstam L G，Parker S L，Einhouse DW，et al.Application of in situ target-strength estimations in lakes：examples from rainbow-smelt surveys in Lakes Erie and Champlain［J］.ICES Journal of Marine Science，2003，60（3）：500-507.

［25］ 孫建富，胡亞東，張劍波.小浪底水庫基于調(diào)水調(diào)沙的漁業(yè)措施研究［J］.中國水產(chǎn)，2014（11）：17-21.

［26］ Tang Y，Lin D，Zhang G，et al.Testing the tracking of behaviour of a caged Takifugu rubripes（Temminck＆amp;Schlegel，1850）using acoustic telemetry［J］.Journal of Applied Ichthyology，2013，29（6）：1456-1458.

［27］ Voegeli F A，Lacroix G L，Anderson JM.Development ofminiature pingers for tracking Atlantic salmon smoltsat sea［J］.Hydrobiologia，1998，371-372：35-46.

［28］ 任玉芹，王珂，段辛斌，等.鳙目標(biāo)強度和行為特征的水聲學(xué)研究［J］.淡水漁業(yè)，2011，41（2）：3-9.

Spatial distribution of fisheries resources in Sanmenxia Reservoir in the Yellow River based on an acoustic method

LIBin1，2，TANG Yong1，SUN Jian-fu3，LIPei-jie1，HU Ya-dong4，ZHANG Jian-bo4
（1.College of Marine Science and Environment，Dalian Ocean University，Dalian 116023，China；2.South China Sea Fisheries Research Institute，Chinese Academy of Fishery Sciences，Guangzhou 510300，China；3.College of Economics and Management，Dalian Ocean University，Dalian 116023，China；4.Fisheries Science Research Institute of Henan Province，Zhengzhou 450044，China）

The fisheries resourceswere investigated from Xiwan to the reservoir dam in the Sanmenxia reservoir using split-beam scientific echosounder（Simrad EY60，70 kHz，Norway）on November 23，2013 to evaluate the utilization of an acousticmethod in assessment of ecological environment in the Yellow Sea River basin.The density，standing crop and spatial distribution in differentareas of fishery resourceswere explored and estimated primarily by fish target strength in situ and using of an echo integration method.The results showed that the acoustic integral value（Nautical area scattering coefficient，NASC）in the reservoirwas5.59m2/nmi2，the average density of fisheries resources0.002 ind./m2，and the size of fisheries resources42.36 t.The surveyed area in the reservoirwas divided into upper，middle and lower reaches（A，B and C areas）according to the density spatial distribution of fisheries resources，and the density of fisheries resources in the total of eight integral segmentswere followed by 0.002 ind./m2in A-Ⅰ，0.017 ind./m2in A-Ⅱ，0.013 ind./m2in A-Ⅲ；0.004 ind./m2in B-Ⅰ，0.008 ind./m2in B-Ⅱand 0.013 ind./m2in C-Ⅰ，0.016 ind./m2in C-Ⅱ，and 0.037 ind./m2in C-Ⅲ.The fisheries resourceswere heavily distributed in the upper reaches of A-Ⅱand A-Ⅲsegments and the lower reaches of reservoir.Most fish stock in the segments of A-Ⅱ，A-Ⅲ，C-Ⅱand C-Ⅱwere vertically distributed in 5 m shallow water layer，and fish stock in the segments of A-Ⅰand B-Ⅱwere primarily distributed in 10m shallow water layer.As fish stock in B-Ⅰsegmentwas distributed in 5m deep water layerwith less amount，the fish stock in C-Ⅲsegmentwas largely distributed in the 10 m to 20 m water layer.The fish target strength in the reservoir wasmostly distributed from-62.5 dB to-47.5 dB，and the individuals with target strength of greater than-44.5 dB，and body length of over 25 cm were entirely distributed under 5 m water layer.The findings prove that the acoustic method reflects the favorable applicability at assessment of fisheries resources in freshwater reservoir，promoting the application of acousticmethod in the study of ecological environmentof fishingwaters in the Yellow River，and providing the theoretical foundation and scientific basic for further study on the acoustic assessment of fisheries resources in inland rivers，lakes and reservoirs.

fishery resources；acoustic assessment；spatial distribution；Sanmenxia Reservoir

S932.4

A

10.16535/j.cnki.dlhyxb.2016.05.016

2095-1388（2016）05-0563-09

2016-01-11

公益性行業(yè)（農(nóng)業(yè)）科研專項（201303050）；國家自然科學(xué)基金資助項目（10774021）

李斌（1990—），男，碩士研究生。E-mail：657904475@qq.com

湯勇（1965—），男，博士，教授。E-mail：tang@dlou.edu.cn