The present status and prospect on exploitotion of tuna and squid fishery resources in South China Sea
-
摘要:
南海中南部深水区蕴藏着丰富的鸢乌贼(Sthenoteuthis oualaniensis)和黄鳍金枪鱼(Thunnus albacares)等大洋性中上层生物资源。文章分别介绍了南海金枪鱼和鸢乌贼的渔业与资源情况,并在此基础上探讨了中国相关渔业的发展前景,并提出南海大洋性生物资源开发应以鸢乌贼灯光罩网渔业为发展重点。
Abstract:There are abundant marine bioresources, such as purpleback flying squid (Sthenoteuthis oualaniensis)and yellowfin tuna (Thunnus albacares), in the deep-sea area of central and southern South China Sea. This paper introduced the status of tuna and squid fisheries and resources in South China Sea respectively.On this basis, the prospect for the development of relevant fisheries in China was discussed and it was proposed that the marine bioresources exploitation in South China Sea should focus on developing flying squid fisheries by using light falling net.
-
水下洗网机采用物理清洗的方法,利用高压水流的能量对网箱网衣附着物进行清洗[1]。在水下洗网机设计中,高压水流常涉及多相混和的介质射流,机理复杂。借计算流体力学的思想和计算机模拟仿真能方便有效地分析流场的动静态特性。这种可视化方法不仅能呈现流场定性,而且能定量地获取流场各处的速度和压力等重要参数[2]。
计算流体力学(computational fluid dynamics,简称CFD),是通过计算机数值计算和图像显示对包含有流体流动相关物理现象的系统所作的分析。通过这种数值模拟,可以得到极其复杂的流场内各个位置上的基本物理量(速度、压力、温度和浓度)的分布,以及这些物理量随时间的变化情况,确定旋涡分布特性、空化特性以及脱流区等。CFD不仅作为一个研究工具,而且还作为设计工具在水利工程、土木工程、环境工程、食品工程、海洋结构工程及工业制造等领域发挥作用[3]。
1. 流场的数学模型
1.1 控制方程
根据计算所得雷诺数大于2 300,流体在洗网机工作腔内的流动是复杂的湍流流动,需要对雷诺时均的质量守恒和动量守恒方程进行求解[4]。这里采用标准k-ε模型作为封闭方程。
质量守恒方程:
$$ \frac{\partial u_i}{\partial x_i}=0 $$ 动量守恒方程:
$$ \frac{\partial u_i}{\partial t}+\frac{\partial\left(u_i u_j\right)}{\partial x_j}=f_i-\frac{1}{\rho} \frac{\partial p}{\rho \partial x_i}+v \frac{\partial^2 u_i}{\partial x_i x_j}-\frac{\partial}{\partial x_i}\left(\overline{u_i^{\prime} u_j^{\prime}}\right) $$ k-ε方程:
$$ \rho \frac{\partial k}{\partial x_i}+\rho \overline{u_j} \frac{\partial k}{\partial x_j}=\frac{\partial}{\partial x_j}\left[\left(\mu+\frac{\mu_t}{\sigma_1}\right) \frac{\partial k}{\partial x_j}\right]+\mu_t\left(\frac{\partial \overline{u_i}}{\partial x_j}\right.\left.+\frac{\partial \overline{u_j}}{\partial x_i}\right) \frac{\partial \overline{u_i}}{\partial x_j}-\rho \varepsilon $$ $$ \rho \frac{\partial \varepsilon}{\partial t}+\rho \overline{u_k} \frac{\partial \varepsilon}{\partial x_k}=\frac{\partial}{\partial x_k}\left[\left(\mu+\frac{\mu_t}{\sigma_2}\right) \frac{\partial k}{\partial x_k}\right]+\frac{C_1 \varepsilon}{k} \mu_t\\\left(\frac{\partial \overline{u_i}}{\partial x_j}+\frac{\partial \overline{u_j}}{\partial x_i}\right) \frac{\partial \overline{u_i}}{\partial x_j}-C_2 \rho \frac{\varepsilon^2}{k} $$ 式中ρ为流体密度(kg·m-3);p为流体压力(Pa);u为流体的速度矢量(m·s-1);μt为流体的动力粘度(Pa·s);ε为耗散率(m2·s-3);k为湍流动能(m2·s-2)。
标准k-ε方程中的相关常数为:Cμ=0.09,C1=1.44,C2=1.92,σ1=1,σ2=1.3。
以上方程组构成求解洗网机内部流场分布规律的封闭方程组,根据实际工况施加相应的边界条件后,构成该方程组的定解问题。
1.2 控制方程的离散
对控制方程的离散常用的数值格式是有限元法(FEM)和有限差分法(FDM) (包括有限体积法)。就其广泛性而言,主要是有限元法。有限元法将连续的结构离散成有限个单元,并在每个单元中设定有限个节点,将连续体离散为一组有限个、且按一定方式相互联结在一起的单元组合体;同时选定场函数的节点值作为基本未知量,并在每一个单元中假设一近似插值函数以表示单元中场函数的分布规律。这样,一个问题的有限元分析中,未知场函数或及其导数在各个节点上的数值就成为新的未知量,即自由度。进而利用力学中的变分原理去建立用以求解节点未知量的有限元方程,于是将一个连续域中的有限自由度问题转化为离散域中的有限自由度问题。随着单元数目的增加,或者随着单元自由度的增加及插值函数精度的提高,解的近似程度也将不断地提高。若单元满足收敛性要求,近似解将趋于收敛为精确解。有限元网格样式多样,大部分问题可以用简单的网格,因而解题自动化程度高。并且能对局部进行网格加密,可以用于复杂边界形状的问题求解[5]。
2. 数值模拟
2.1 有限元模型的建立
流体在洗网机内部的流动实际上是轴对称流动,考虑到流动的特性,将二维流体作为此次模拟的研究对象。根据实际情况,为提高计算效率,简化模型,将工程设计图(图 1-a)转化为ANSYS前处理的几何图(图 1-b)。因结构简单,采用自底向上的方式构建几何模型。
2.2 网格的划分
由于考虑到模拟的精确性,故设置单元尺寸且关键位置处的网格数不少于3个,以保证运算结果的精度。最后的CFD计算模型网格划分如图 1-c所示,共有单元390个,节点474个。在建立模型时,所取单元类型为FLUID141。划分网格时,由于流场在洗网机内腔弯曲处以及洗网机喷嘴出口处比较复杂,因此,在复杂的地方进行了较细密的网格划分,以获得更加精确的模拟结果。
2.3 初始条件的确定
假定流体不可压缩,并且其性质为恒值。选取海水为流动介质,可认为是定常、连续和不可压缩流体。20 ℃常温下海水的平均密度为1.023×103 kg·m-3,动力粘度为1.345×10-3 kg·(m·s)-1。计算设计工况(流量=100 L·min-1)下的流场变化。假设进口速度均匀,初速度=8.3 m·s-1,且垂直于进口流场方向的流体速度为零,出口给定相对压力p=0。由计算可知流体的雷诺数大于2 300,所以流场是湍流模型。为了得到较为精确的结果,使用映射网格划分,能在边界上更好地保持恒定的网格特性[6]。
计算中所使用的洗网机设计工况、流体性质等参数见表 1。
表 1 洗网机工况参数以及流体性质参数Table 1. Work condition parameters and fluid properties of net cleaning machine流量/L·min-1
flow flux海水密度/kg·m-3
sea-water density动力粘度/kg·(m·s)-1
kinematic viscosity初始速度/m·s-1
initial velocity出口相对压力/Pa
relative pressure of outlet100 1.023×103 1.345×10-3 8.3 0 2.4 边界条件
由偏微分方程的理论可知,任何一个方程或方程组具有无数组可能的解,若要获得确切的定解,就必须给定具体流动问题完全的定解条件,即边界条件。在计算流体力学中,边界条件起着重要作用,施加边界条件的准确程度将最终决定洗网机系统流场模拟的可靠性[7]。总的来说,洗网机内部流场的控制方程边界条件的施加有4个部分:进口边界条件、出口边界条件、壁面边界条件和相对压力边界条件。
在洗网机流场的CFD分析中,边界条件包括速度自由度和压力。边界的条件设定为:1)进口采用速度进口,在计算域的进口边界上,给定初始速度。假设入口速度均匀,且垂直于进口流场方向的流体速度为零。入口处施加流体的入口速度,其值为常数;2)出口处施加流体的相对压力边界条件,为常数值零;3)内壁施加无滑移边界条件,所有速度分量为零。施加边界条件后见图 2。
3. 模拟结果分析
仿真分析收敛监测图见图 3。分析中,用变量速度VX,VY、压力(PRES)、湍流动能(ENKE自由度)及动能耗散率(ENDS自由度)等5项参数对收敛情况监控。横坐标表示迭代次数40,纵坐标表示收敛标准值,收敛监测表明求解变量的归一化改变率[8]。从图 3可以看出各项的震荡都不是很明显,在20步的时候已基本完全收敛。
流场整体速度差别不大,在直角拐弯处流速明显减小,出口处的速度为最大(图 4)。因为流体本身具有粘性,所以在主流的带动下产生涡旋,在主流中吸取能量,引起能量的损失,降低能量的利用率;此外涡旋现象也增加了流体的噪声。在直角拐弯处速度的变化率过大,引起很大的稳态及瞬态的液动力,因而较容易产生内壁磨损。流场的总压力在进口处最大,从直角拐弯处开始递减,至出口处最小(图 5)。由于流体离心力作用和直角拐弯处有离心力的作用,从而产生二次涡流效应,在流体流经直角拐弯处的时候改变了流体的速度大小和方向,造成了压力损失。但在出口中间位置,压力明显增大,两侧压力变小。
为了更好地分析洗网机喷嘴出口速度分布情况,进行路径绘制[9]。洗网机喷嘴出口路径上的速度分布曲线(横坐标为出口路径,纵坐标为出口速度)见图 6,洗网机喷嘴出口的压力图见图 7。洗网机出口处形成了明显的速度梯度,由于流场内存在紊流,所以在出口的上下两侧形成了旋涡(图 6)。在出口处产生了较大的流体阻力带来的损失,降低了能量的利用。出口中间出现高压,在喷嘴的出口垂直面的外侧产生很大的负压区,因此,在这些位置会出现旋涡并容易产生气蚀现象(图 7)。
4. 结论
1) 利用ANSYS软件对洗网机内的流场进行了模拟仿真,对洗网机流场的压力分布和速度分布进行了分析,获得了大量的流场信息,较深入地了解了洗网机的结构参数对流动特性的影响,为洗网机结构的合理设计、材料的选用及流体流动方向提供了技术支持,对结构的优化设计有着指导意义。
2) ANSYS可以很好地了解用传统测试和物理试验都很难研究的复杂流体流动现象,有助于提高设计质量和效率,优化整机的开发。
-
表 1 调查渔船鸢乌贼捕捞情况
Table 1 General situation of S.oualaniensis fishing
年份
year渔期/月
fishing season航次
sail作业天数/d
fishing days产量/t
catch日均产量/t·d-1
average catch价格/RMB·kg-1
price产值/×104RMB
value日均产值/RMB·kg-1
average value2006 1~5 9 91 197.1 2.17 4 78.8 8 663 2007 3~5 5 57 127.4 2.24 4 51.0 8 942 2008 3~5 6 60 114.7 1.91 5 57.4 9 560 2009 3~5 3 24 56.7 2.36 2 11.3 4 727 -
[1] 陈新军, 俞存根, 卢伙胜. 渔业资源与渔场学[M]. 北京: 海洋出版社, 2004. https://www.zhangqiaokeyan.com/book-cn/081501459619.html [2] 陈炎, 陈丕茂. 南沙群岛金枪鱼资源初探[J]. 远洋渔业, 2000(2): 7-10. https://www.cqvip.com/QK/82064X/20002/4575068.html [3] 蔡怡卉. 台湾近海鲔延绳钓渔海况变动之研究[D]. 台北: 国立台湾海洋大学, 2005. [4] 张水锴. 台湾鲔延绳钓渔业及重要管理措施简介[J]. 中国水产(台), 2004(621): 15-28. [5] 詹铭慧. 油价上涨对台湾渔业部门之影响评估——渔业部门均衡模型之应用[D]. 台北: 国立台湾海洋大学应用经济研究所, 2007. [6] 黄凌. 台湾近海鲔延绳钓渔获效能之研究[D]. 台北: 国立台湾海洋大学, 2004. [7] 刘佩妤. 以台湾近海鲔延绳钓渔业大目鲔体长数据估计其年龄、成长与死亡率以及单位加入量分析[D]. 台北: 国立台湾大学海洋研究所, 2005. [8] 苏楠杰. 台湾近海金枪鱼延绳钓渔业黄鳍金枪鱼的年龄、成长、死亡率与单位加入生产量[D]. 台北: 国立台湾大学海洋研究所, 2003. [9] LEVIS A D. The tuna fisheries of vietnam-an overview of available information[C/OL]//1st Meeting of the Scientific Committee of the Western and Central Pacific Fisheries Commission Noumea, New Caledonia, August 8-19, 2005[2005-08-02]. http://www.spc.int/OceanFish/Html/WCPFC/sc1/pdf/SC1_ST_IP_5.pdf"> http://www.spc.int/OceanFish/Html/WCPFC/sc1/pdf/SC1_ST_IP_5.pdf
[10] NGAY. Oceanic tuna fishing in Vietnam[EB/OL] [2004-04-26]. http://www.mofi.gov.vn"> http://www.mofi.gov.vn.
[11] LE K L, OLA F, NGUYEN T K A. Economic performance of open-access offshore fisheries the case of Vietnamese longliners in the South China Sea[J]. Fish Res, 2008, 93(2): 296-304. https://www.sciencedirect.com/science/article/pii/S0165783608001616
[12] HOANG T. Ocean tuna fishing and marketing in Viet Nam[J]. Vietfish, 2009, 6(1): 56-59. https://thefishsite.com/articles/ocean-tuna-ranching-in-viet-nam
[13] SIRIPORN P. Tuna longline fishing ground analysis by satellite image [EB/OL] [2009-08-07]. http://mail.vnio.org.vn:8088/traireport/Tuna%20fishing%20ground%20analysis%20by%20satellite%20image_Siriporn.pdf"> http://mail.vnio.org.vn:8088/traireport/Tuna%20fishing%20ground%20analysis%20by%20satellite%20image_Siriporn.pdf.
[14] HA Y. Tuna trawlers look to international waters[EB/OL] [2008-12-11]. http://english.vietnamnet.vn/biz/2006/08/600621/"> http://english.vietnamnet.vn/biz/2006/08/600621/.
[15] 杨吝, 张旭丰, 谭永光, 等. 南海北部灯光罩网渔获组成及其对渔业资源的影响[J]. 南方水产, 2009, 5(4): 41-46. doi: 10.3969/j.issn.1673-2227.2009.04.008 [16] 刘必林. 利用耳石微结构研究印度洋西北海域莺乌贼的年龄和生长[D]. 上海: 上海水产大学, 2006. https://d.wanfangdata.com.cn/thesis/Y918043 [17] LABE L L. Catch rate of oceanic squid by jigging method in the South China Sea area 3: western Philippines: proceedings of the Third Technical Seminar on Marine Fishery Resources Survey in the South China Sea, Area 3: Western Philippines, July 13-15, 1999[C]. Bangkok, Thailand: Southeast Asian Fisheries Development Center, 2000: 19-31.
[18] SIRIRAKSOPHON S, NAKAMURA Y. Ecological aspects of the purpleback flying squid, Sthenoteuthis oualaniensis (Lesson) in the west coast of Philippines: International Conference on the International Oceanographic Data and Information Exchange in the Western Pacific (IODE-WESTPAC) - ICIWP ′99, Langkawi (Malaysia), November 1-4, 1999[C]. Bangkok, Thailand: Southeast Asian Fisheries Development Center, 2001: 187-194.
[19] SIRIRAKSOPHON S, NAKAMURA Y, NATINEE S. Exploration of purpleback flying squid, Sthenoteuthis oualaniensis resources in the South China Sea[M]. Samutprakan, Thailand: Southeast Asian Fisheries Development Center, 2001: 1-81. https://www.semanticscholar.org/paper/Exploration-of-Purpleback-Flying-Squid%2C-Resources-Siriraksophon-Nakamura/07f6dc3a930fda7287c30f433cad23b03e34cec1
[20] 贾晓平, 李永振, 李纯厚, 等. 南海专属经济区和大陆架渔业生态环境与渔业资源[M]. 北京: 科学出版社, 2004: 389-391, 394. https://www.las.ac.cn/front/book/detail?id=42fe2fc1c537854dac87610b13d7d5a0 [21] 张引. Fisheries acoustic studies on the purpleback flying squid resource in the South China Sea[D]. 台北: 国立台湾大学海洋研究所, 2005. [22] DICKSON J O, RAMISCAL R V, ESCOBAR S, Jr. Tuna resource exploration with longline in the South China Sea, area 3: Western Philippines: proceedings of the Third Technical Seminar on Marine Fishery Resources Survey in the South China Sea, Area 3: Western Philippines, July 13-15, 1999 [C]. Bangkok, Thailand: Southeast Asian Fisheries Development Center, 2000: 39-48.
[23] AUSSANEE M, PRATAKPHOL P. Tuna resource exploration with tuna longline in the South China Sea, area 4: Vietnamese waters: proceedings of the Fourth Technical Seminar on Marine Fishery Resources Survey in the South China Sea, Area 4: Vietnamese waters, Hai Phong (Vietnam), September 18-20, 2000[C]. Bangkok, Thailand: Southeast Asian Fisheries Development Center, 2000.
[24] 王健雄. 台湾近海鲔钓渔业资源动态解析[J]. 中国水产(台), 1986(400): 3-15. [25] PHAN H D. The interaction between the tuna fishing operation and sea turtles issue in Viet Nam[C/OL]//Proceedings of the 2nd International Symposium on SEASTAR2000 and Asian Bio-logging Science (The 6th SEASTAR2000 Workshop), December 13-14, 2005, Bangkok, Thailand[2005-12-14]. http://repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/44089/1/6thSEASTAR_79.pdf"> http://repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/44089/1/6thSEASTAR_79.pdf
[26] 王尉任. 西太平洋黄鳍鲔生殖生物学研究[D]. 台北: 国立台湾大学海洋研究所, 2004. [27] 朱淑玲. 西太平洋大目鲔之生殖生物学研究[D]. 台北: 国立台湾大学海洋研究所, 1998. [28] LANGLEY A, HAMPTON J, KLEIBER P. Stock assessment of yellowfin tuna in the western and central Pacific Ocean, including an analysis of management options[C/OL]//Scientific Committee Third Regular Session, August 2007, Honolulu, USA[2007-12-07]. http://www.wcpfc.int/system/files/documents/meetings/scientific-committee/3rd-regular-session/stock-assessment-swg-working-papers/WCPFC-SC3%20SA-SWG%20WP-01.pdf"> http://www.wcpfc.int/system/files/documents/meetings/scientific-committee/3rd-regular-session/stock-assessment-swg-working-papers/WCPFC-SC3%20SA-SWG%20WP-01.pdf.
[29] LANGLEY A, HAMPTONL J, KLEIBER P, et al. Stock assessment of bigeye tuna in the western and central Pacific Ocean, including an analysis of management options[C/OL]//Scientific Committee Fourth Regular Session, August 2008, Papua New Guinea [2008-07-31]. http://www.wcpfc.int/system/files/documents/meetings/scientific-committee/4th-regular-session/stock-assessment-specialist-working-gro-/SC4-SA-WP1-rev.1%20%5B-BET%20Assessment%5D.pdf"> http://www.wcpfc.int/system/files/documents/meetings/scientific-committee/4th-regular-session/stock-assessment-specialist-working-gro-/SC4-SA-WP1-rev.1%20%5B-BET%20Assessment%5D.pdf.
[30] CHEN X, LIU B, CHEN Y. A review of the development of Chinese distant-water squid jigging fisheries[J]. Fish Res, 2008, 89(3): 211-221. https://www.researchgate.net/publication/222570360_A_review_of_the_development_of_Chinese_distant-water_squid_jigging_fisheries
-
期刊类型引用(2)
1. 闫国琦, 倪小辉, 莫嘉嗣. 深远海养殖装备技术研究现状与发展趋势. 大连海洋大学学报. 2018(01): 123-129 . 百度学术
2. 郭少儒, 张尚文, 郭长瑞, 周少斌, 文晓龙, 郭煜, 席良贤, 苏厚德. LNG开架式海水气化器强化传热过程数值模拟. 石油化工设备. 2013(06): 50-53 . 百度学术
其他类型引用(0)