大亚湾水交换的数值模拟研究

王聪, 林军, 陈丕茂, 章守宇

王聪, 林军, 陈丕茂, 章守宇. 大亚湾水交换的数值模拟研究[J]. 南方水产科学, 2008, 4(4): 8-15.
引用本文: 王聪, 林军, 陈丕茂, 章守宇. 大亚湾水交换的数值模拟研究[J]. 南方水产科学, 2008, 4(4): 8-15.
WANG Cong, LIN Jun, CHEN Pi-mao, ZHANG Shou-yu. Numerical simulation on water exchange in Daya Bay[J]. South China Fisheries Science, 2008, 4(4): 8-15.
Citation: WANG Cong, LIN Jun, CHEN Pi-mao, ZHANG Shou-yu. Numerical simulation on water exchange in Daya Bay[J]. South China Fisheries Science, 2008, 4(4): 8-15.

大亚湾水交换的数值模拟研究

基金项目: 

国家高技术研究发展计划(863计划)重大项目 2006AA100303

上海市教委重点基金项目 05ZZ51

国家自然科学基金项目 30471332

农业部渔业生态环境重点开放实验室开放基金项目 2006-4

详细信息
    作者简介:

    王聪(1982-),女,硕士研究生,从事海洋生态系统工程研究。E-mail: cwang@stmail.shfu.edu.cn

    通讯作者:

    林军,E-mail: jlin@shfu.edu.cn

  • 中图分类号: P731

Numerical simulation on water exchange in Daya Bay

  • 摘要:

    对大亚湾的潮流场进行三维数值模拟, 分别用Lagrange质点追踪法和保守物质输运扩散2种方法研究大亚湾的水交换能力, 并把大亚湾分成7个子区域, 计算整个海湾和各子区域的水体平均存留时间和更新时间, 比较各子区域的水交换能力的大小, 分析2种方法结果差异产生的原因。大亚湾的水交换能力受潮致余流流场结构支配。质点追踪模拟结果显示, 大亚湾南部及大鹏澳区域质点迁移出海湾的速度较快, 哑铃湾、大亚湾北部及范和港区域较慢; 保守物质输运扩散模型计算结果显示, 大亚湾东南部水体的更新时间最快, 水交换能力最强, 其次分别为大亚湾中东部、西南部、范和港、大亚湾北部、大鹏澳、哑铃湾。

    Abstract:

    The tide and currents in Daya Bay was simulated by a 3D numerical estuarine and coastal ocean model. Lagrangian particle tracking and matter transportation and diffusion were applied in water exchange study. Daya Bay was divided into seven parts to calculate the averaged water residence time and refresh time in the whole bay as well as in each part. The structure of the tidal residual currents, as well as landform and depth play decisive roles in water exchange in Daya Bay. The result showed that the discharge of particles was faster in south of Daya Bay and Dapeng′ao and slower in Yalin gulf, north of Daya Bay and Fanhe gulf as indicated by lagrangian particle tracking. Water refresh time was shorter in southeast of Daya Bay than in middle east, southwest of Daya Bay, Fanhe gulf, north of Daya Bay, Dapeng′ao and Yalin gulf.

  • 占地球表面积71%的海洋蕴藏着丰富的生物、矿物和能源等资源,是人类繁衍生息和持续发展的重要基础[1-2]。我国是海洋大国,丰富的海洋资源为水产养殖提供了天然的资源基础[3]。树立大食物观,建设海洋牧场,是推进我国粮食安全战略的重要保障。近年来,我国重视发展深远海绿色养殖,鼓励深远海大型智能化养殖平台建设,使深远海养殖逐渐朝机械化、自动化、信息化方向发展[4-5]。其中,最具代表性的是养殖工船养殖装备的设计与研发。养殖舱内的水质检测是养殖工船作业的关键环节。水质环境中温度、pH、溶解氧 (DO)、氨氮等是影响鱼类生长的重要因素[6-8],实时水质监测可控制养殖舱的水环境,降低水产养殖风险,提高养殖效益。

    水下机器人又称无人遥控潜水器,是一种在水下作业的机器人,它可在常规技术不能到达的深度或区域进行综合考察和研究,能自主完成多种作业使命,成为探索、开发和利用海洋的重要手段之一[9-10]。目前,水下机器人已应用于多场景海洋作业中,如水下矿物、矿石和水下考古观测[11],海底电缆铺设[12],石油和天然气工业[13],在军事、渔业等领域也占据举足轻重的地位[14-16]。在水产养殖中,利用水下机器人进行水质监测,可明显提升监测的灵活性、效率和范围。王懿偲等[17]设计了一款监测水质与鱼群活动的机械鱼机器人,搭建了机器鱼三维路径跟踪控制系统和误差模型,以及控制精度高和鲁棒性强的模糊滑模控制器,但鉴于实际试验环境中干扰更加复杂,需进一步优化算法。刘宁等[18]开发了一款用于定期水质监测的新型蛇形机器人,采用仿生蛇形和沉浮结构设计,分析和模拟了蛇形机器人的稳定性,为水下机器人的结构设计提供了全新的思路,但其无法满足定点悬停采样以及实时监测。李翔飞和戴逸飞[19]构建了一款用于监测水下环境的六足水下机器人,其控制系统使用Python程序对步态进行控制,采用仿生结构将水底行走与水中游动相结合,可代替人工水下作业,但其作业深度受限。Ravalli等[20]开发了一种用于水质监测和预警的多传感器仿生鱼机器人,在结构上将pH传感器与改良后的聚苯胺传感器同水下机器人相结合,将采集的信号转换为电子信号控制水下机器人,利用多传感器平台对环境中的氢离子 (H+) 浓度进行检测,但检测的参数较少。Jo等[21]提出了一种新的、完全开源、低成本、小型无人水面飞行器 (Unmanned Surface Vessel, USV),用于实时测量近地表水质,USV设计的硬件和软件均为开源,开发Android程序通过蓝牙控制USV,可实现自主导航并显示传感器和GPS数据,但其使用场景受限。Huang等[22]设计了一种用于水产养殖水质检测的机器人,采用开放式框架结构设计、三螺旋桨推进方式,机器人配备水质检测模块,采集定点位置的淡水水温、pH、DO等参数,但实际作业中无法实现实时监测。

    目前,水质检测技术正向多元化、智能化、信息化方向发展[23-24],但存在技术不成熟、运营成本高、普及范围小、水质污染严重等问题[25]。现今水质监测主要靠人工采集,再进行仪器检测,虽然准确度高,但效率低且周期长;另一种方式是采用传感器网络技术,在检测区域布置各类水质监测传感器,可以实现实时检测,却无法实现水质分层监测,而且传感器费用较高,作业区域也受限。基于上述问题,本研究设计了一款流线型作业型缆控水下机器人 (Remotely operated vehicle, ROV),将ROV和多参数水质监测传感器组结合,搭建了一套完整的水质检测系统,通过抽水方式将采集水样输送到采样桶,然后由多参数传感器组和上位机软件实时监测水温、pH、DO及盐度参数,可解决水产养殖中水质检测结果的片面性、局部性等问题,且能实现5自由度灵活运动,可在不同深度下完成分层水质抽样检测;利用多参数传感器监测系统可以实时监测,同时具有观察鱼类生长环境和监控作业的功能,可提升养殖成效,助力海洋牧场建设。

    ROV整体水质检测系统由多参数水质监测传感器组和水下机器人组成,多参数水质监测传感器组检测采样的水质,水下机器人 (ROV) 为采样的载体。该水质监测系统作业过程为:在养殖舱甲板或岸上安装多参数水质监测传感器组,用水管搭载到机器人载体上,依靠机器人的高机动性通过手柄操控在养殖舱内到达不同的水层,经自吸泵吸水输送到传感器组,对水质的DO、温度、pH、盐度进行监测,利用上位机软件实现水质参数分析、记录和储存功能。整体水质检测系统流程见图1-a。根据3种传感器的结构尺寸设计对应的固定结构,将传感器组固定在同一水平线上,保证水质参数采集数据的一致性和准确性,将3个传感器通过螺栓和螺母固定在空隙位置,可以根据水位高低调节传感器组的位置,可以更换不同类型水质传感器获取不同的水质参数。整体传感器组见图1-b。

    图  1  ROV水质检测系统设计
    Figure  1.  Design of ROV water quality detection system

    ROV整体采用流线型外壳结构设计,由6个推进器、2k分辨率的摄像头和可调节亮度的照明灯等组成。根据设计指标,要求ROV可以满足多自由度运动包括升潜、转舵、进退、回转和横摇等5个自由度;可以实现在指定位置悬停并配合水管进行抽水作业;能够拍摄水下视频并传输记录水底场景,同时记录ROV水下姿态角、水下深度等信息。具体设计指标见表1

    表  1  ROV整体设计参数指标
    Table  1.  Parameter indicators of ROV overall design
    项目 Project参数 Parameter
    本体结构 Body structure 流线型 Streamline
    结构尺寸 Structure size 545 mm×508 mm×172 mm
    作业深度 Operation depth/m 0~50
    工作航速 Operation speed/kn 4
    空气中质量 Mass in air/kg 6.92
    控制方式 Control method 手动(遥控)控制
    Manual (Remote) control
    推进器布置 Thruster arrangement 水平2个,竖直4个
    Horizontal 2, Vertical 4
    工作时间 Operation time/min 45
    LED照明灯 LED lighting 4档可调节亮度
    Adjustable brightness in 4 levels
    下载: 导出CSV 
    | 显示表格

    ROV通过内部蓄电池为其供电,通过线缆为上位机发送信号指令并将收集的图像、传感器信号以及位置等信息传送给上位机,通过人机交互界面直观、清楚、准确地反映水下环境。ROV主要由水下机器人本体、线缆以及PC端上位机组成。ROV系统主要由动力系统、观测及照明系统、控制系统、外壳结构和抽水检测结构等组成。ROV具体系统组成及系统设计流程图分别见图2图3

    图  2  ROV系统组成
    Figure  2.  ROV system composition
    图  3  ROV系统设计流程图
    Figure  3.  Flowchart of ROV system design

    小型水下机器人具有运动速度敏捷、机动性高等特点,外壳作为水下机器人的重要载体与水直接接触,流线型外壳设计更好地满足了运动要求。由于本研究设计的外壳模型呈现不规则形状较多,且对水下压力要求较小,结合实际需求,选择采用高强度ABS[26-27]作为外壳材料,参数见表2。采用Solidworks软件从美观、性能、实用性等角度出发设计外壳,外壳的主体结构包含摄像头和照明灯外壳、抽水外壳、推进器的外壳和固定结构。采用光固化3D打印技术制造外壳 (图4)。

    表  2  ABS物理参数
    Table  2.  ABS physical parameters
    参数 Parameter数值 Value
    弹性模量 Elastic modulus E/MPa 2×103
    抗剪模量 Shear modulus/MPa 318.9
    质量密度 Mass density ρ/(kg·m–3) 1.02×103
    泊松比 Poisson's ratio μ 0.394
    屈服强度 Yield strength/MPa 46
    下载: 导出CSV 
    | 显示表格
    图  4  ROV外壳
    Figure  4.  Housing of ROV

    为验证外壳设计是否满足强度要求,对ROV外壳进行应力仿真分析。利用Solidworks软件进行建模并保存为IGES格式导入Ansys中;将上、下外壳定义为刚体,并定义材料属性为ABS材料,参数设置见表2;采用自动网格划分方法划分网格;施加压力载荷,需满足水下机器人极限作业深度50 m,经计算施加载荷压力为0.5 MPa;将上、下盖两端的接口处添加固定支撑进行求解。求解后上、下外壳等效应力和总变形见图5

    图  5  ROV上、下外壳数值模拟
    Figure  5.  Numerical simulation of upper and lower housing of ROV

    仿真分析表明,水下机器人的上、下外壳最大总变形分别为0.061 382和0.022 448 mm (图5-a, 5-b),属于微小变形,可忽略不计;上、下外壳最大等效应力 (MPa) 分别为 11.01$ \leqslant {{\sigma }}_{\mathrm{s}}=46 $$ 4.192\,3\leqslant {\sigma }_{\mathrm{s}}=46 $ (图5-c, 5-d),可满足ABS材料屈服强度要求。根据数值模拟可知,外壳设计满足应力要求,说明整体外壳能满足水下50 m的强度和压力要求。

    ROV动力系统组主要由电机、推进器、供电系统组成,通过蓄电池为电机提供动能带动推进器运动,为ROV提供在水中推进和移动的动力支持。ROV整体布局采用6个推进器总体设计,4个在竖直方向和2个在水平方向,在水平和竖直方向上的2个推进器螺旋桨转速不同,还能够产生一定的力矩[28],竖直方向4个推进器提供浮潜、纵倾、俯仰和转舵4个自由度和水平方向2个推进器为机器人提供进退自由度。为满足ROV水下工作条件,采用一正一反6个推进器,6个推进器布置在同一轴线上,利用板材将6个推进器固定在耐压舱两端,推进器的工作的平均功率在270 W,最大推力为700 N。

    在竖直方向推进器的推力和扭矩表达式为:

    $$ \left[\begin{array}{ll} T & M \end{array}\right]^T{\text{=}}\left[\begin{array}{llllll} T_1{\text{+}}T_2 & 0 & 0 & 0 & a\left(T_1{\text{−}}T_2\right) & 0 \end{array}\right]^T$$ (1)

    在水平方向推进器的推力和扭矩表达式为:

    $$ \begin{array}{c} \left[\begin{array}{ll} T & M \end{array}\right]^T{\text{=}} \left[\begin{array}{llllll} 0 & 0 & T_3{\text{+}}T_4{\text{+}}T_5{\text{+}}T_6 & b\left(T_3{\text{−}}T_4+T_5{\text{−}}T_6\right) & 0 & a\left(T_1{\text{−}}T_2\right) \end{array}\right]^T \end{array} $$ (2)

    由公式(1)、(2)可知,整体ROV的推力和扭矩表达式为:

    $$ \left[\begin{array}{ll} T & M \end{array}\right]^T{\text{=}}\left[\begin{array}{llllll} T_1+T_2 & 0 & T_3{\text{+}}T_4{\text{+}}T_5{\text{+}}T_6 & b\left(T_3{\text{−}}T_4{\text{+}}T_5{\text{−}}T_6\right) & a\left(T_1{\text{−}}T_2\right) & a\left(T_1{\text{−}}T_2\right) \end{array}\right]^T$$ (3)

    式中:$T_1 $$T_2 $x轴上推进器空间推力;$T_3 $$T_4 $$T_5 $$T_6 $z轴上推进器空间推力;T为表示空间推力;M为力矩。

    ROV进退自由的推力为$ {T}_{1}{\text{+}}{T}_{2} $,浮潜推力为$ {T}_{3}{\text{+}}{T}_{4}{\text{+}}{T}_{5}{{\text{+}}T}_{6} $,纵倾力矩为 $ b\;({T}_{3}{\text{−}}{T}_{4}{\text{+}}{T}_{5}{{\text{−}}T}_{6}) $,俯仰力矩为 $ a\;({T}_{1}{\text{−}}{T}_{2}) $ 和转舵力矩为 $ a\;({T}_{1}{\text{−}}{T}_{2}) $。总体的分布见图6

    图  6  推进器分布图
    Figure  6.  Thruster setup diagram

    ROV控制系统是整个系统的心脏,控制ROV的执行作业。控制硬件主要有主控制器Pixhawk、树莓派、电力载波模块、线缆和地面基站等组成。ROV整体控制系统需要使用3.3、5、12 V 3种电压,通过内置的24 V可充电锂电池,经由多路DC-DC电压转换模块分压即可获得相应电源。其中,5 V电压供电给主控制器、树莓派和电力载波模块,推进器电机总功率需要在额定电压24 V下,3.3 V为整体的航姿参考系统供电;主控制器Pixhawk由两片芯片STM32F427和STM32F100组成,STM32F427为主芯片负责各类传感器数据采集和处理,负责视频传输和通信;ROV的加速度计负责采集各坐标轴角度、角速度及角加速度;深度信息利用深度计通过压力信息解析获取,同时可以获取海水温度,实现温度计的功能。推进器采用差速P75推进器,使用直流无刷电机并集成了电子调速器,经由控制器STM32F100产生PWM信号对水下6个推进器进行控制,实现推进器正反转的控制和调速功能;主控制器和树莓派通过USB连接进行通信并传递相应控制信息,树莓派3B (Raspberry Pi 3B)微型计算机搭载有1.2 Ghz 64位ARM Cortex-A53四核处理器,通过线缆连接到水下电力载波模块Fathom X,Fathom X以太网通讯板可利用双绞线完成长距离的以太网通讯,其中一个布置于水下机器人本体,另一个布置于地面操作系统内,二者均用于标准以太网通讯与双绞线以太网通讯的转化,水下电力载波模块将树莓派的视频信号转换到两根电力线上进行通信,水面电力载波模块将从电力线上传输过来的视频信号转换回网线接口连接至路由器上,以完成高清视频数据的采集发送以及以太网接口下的数据传输。整体的控制原理见图7

    图  7  ROV控制系统原理图
    Figure  7.  Schematic diagram of ROV control system

    姿态和深度控制是ROV运动控制系统的重点和难点,其控制效果直接影响机器人的性能。ROV的控制从全局上可以理解为对电机进行控制,控制电机力矩的输入和输出。本研究采用非集中控制方式,通过4个独立的控制器来对ROV的姿态和深度进行控制,其由艏向角控制器、纵倾角控制器、横倾角控制器和深度控制器组成 (图8)。其中,艏向角控制器用于控制机器人在水中艏向角的平衡,它的输出量分别等量且相反控制机器人两侧2个垂直推进电机;纵倾角控制器用于控制机器人在水中纵倾角的平衡;横倾角控制器用于控制机器人的偏航角;深度控制器用于控制下潜深度。4个姿态控制器均采用PD控制器。4个姿态和深度控制器的输入分别为姿态角和深度的期望值,由解算得到的姿态角和深度实际值的比较偏差,通过上位机软件显示再进行调整。当需要定深时,以深度期望值与压力传感器测量计算得到的深度实际值比较得到的偏差,作为深度控制器的输入量。深度控制器输出深度控制量为PWM信号控制电调增量,电调控制水平方向的4个推进器进而控制ROV的下潜深度。

    图  8  ROV姿态和深度控制原理图
    Figure  8.  Schematic diagram of ROV attitude and depth control

    在ROV外壳中端位置上采用抛物线线型设计3口转换吸水管口,通过转换头1连接水管,2和3为吸水口。整体抽水采用抛物线线型设计有利于减轻流体阻力,管口整体尺寸偏小长度为30 mm,宽度为60 mm,高度为45 mm,3个管口半径为10 mm见 (图9)。

    图  9  抽水外壳
    Figure  9.  Drawing of pumping housing

    ROV的固定结构主要包括整体耐压舱的固定、6个推进器同轴线固定、摄像机和照明灯固定 (图10)。其中摄像机和照明灯采用环形夹具和卡扣配合螺母固定;耐压舱舱体采用密封舱安装座、安装座固定环和安装座支撑杆组成的固定环形柱固定,防止机器人在运动时耐压舱移位,顶端位置与外壳通过螺母固定;推进器采用板材固定,板材一端和固定环形柱通过螺丝固定,另一端等距对称固定6个推进器。

    图  10  ROV的主要的固定结构
    Figure  10.  Main fixed structures of ROV

    对上述固定部件进行有限元分析,判断是否满足设计应力要求。耐压舱的固定部件材料为铝合金,摄像机和照明灯的固定部件材料为结构钢,推进器固定板材采用胶木材料。3种材料的参数见表3

    表  3  3种材料参数
    Table  3.  Parameters of three materials
    参数
    Parameter
    结构钢
    Structural steel
    铝合金
    Aluminium
    胶木
    Plywood
    密度 Density/(kg·m–3) 7 850 2 770 1 450
    泊松比 Poisson's ratio μ 0.3 0.33 0
    杨氏模量 Young's modulus/MPa 2.0×105 7.1×104 0
    剪切模量 Shear modulus/MPa 7.69×104 2.67×104 0
    拉伸屈服强度 Tensile yield strength/MPa 250 280 ≥100
    抗压极限强度 Ultimate compressive strength/MPa 460 310 0
    垂直层向弯曲强度 Vertical laminar bending strength/MPa 0 0 ≥120
    垂直层向压缩强度 Vertical laminar compression strength/MPa 0 0 ≥250
    下载: 导出CSV 
    | 显示表格

    对所有固定部件进行前处理,网格划分及外部载荷施加,施加压力为0.5 MPa,进行求解。求解后数值模拟云图见图11

    图  11  各个固定部件数值模拟云图
    Figure  11.  Numerical simulation of each fixed part

    由固定部件的仿真结果可知,图11-a—11-b照明灯固定部件的最大变形度为0.019 449 mm,变形度较小,可忽略不计,等效应力为171.44 MPa,小于结构钢的屈服强度 (250 MPa),满足50 m的压力要求;图11-c—11-d摄像机的固定部件的最大变形度为0.057 932 mm,变形度较小,等效应力为228.71 MPa,小于结构钢的屈服强度 (250 MPa),满足50 m的压力要求;图11-e—11-f耐压舱的固定部件的最大变形度为0.268 74 mm,变形度较小,等效应力为235.28 MPa,小于铝合金的屈服强度 (280 MPa),满足50 m的压力要求;图11-g—11-h 6个螺旋桨固定部件的最大变形度为0.009 344 6,变形度较小,等效应力为44.274 MPa,小于胶木的垂直层向弯曲强度和屈服强度(100 MPa),满足50 m的压力要求。综上所述,ROV所有的固定部件的强度和设计均满足水下50 m的要求。

    耐压舱内装载了ROV核心部件,包括控制板、电源板、电子罗盘、水下摄像机、照明器(LED等)、深度传感器、温度传感器和各类连接线等。耐压舱需具备密封保护电子元件的作用,基于流线型的外壳将耐压舱设计为圆柱形状,采用了亚克力材料(PMMA),具体参数见表4。首先对整个耐压舱舱体强度、安全系数等进行理论计算[29-30],然后通过Ansys软件的静应力分析判断舱体应力、应变及安全系数是否满足在水中的强度和刚度要求。

    表  4  亚克力材料属性
    Table  4.  Properties of acrylic material
    参数 Parameter数值 Value
    密度 Density/(kg·m–3) 1 190
    泊松比 Poisson's ratio μ 0.35
    抗剪模量 Shear modulus G/MPa 2051
    张力强度 Tension strength σt/MPa 61
    弹性模量 Elastic modulus E/MPa 3 160
    屈服强度 Yielding strength σs/MPa 45
    下载: 导出CSV 
    | 显示表格

    在标准大气压下,压强的计算公式为:

    $$ \begin{array}{c}P{\text{=}}\rho gh{\text{=}}0.5\left({\rm{MPa}}\right)\end{array} $$ (4)

    式中:ρ为海水密度 (kg·m–3);g为重力加速度 (9.8 N·kg–1);h为作业深度 (m)。

    由劳氏规范安全系数知,PMMA材料的安全系数最小为2.5[31],目前设计的耐压舱的厚度为5 mm,长度L = 280 mm,直径D = 100 mm,ROV潜水深度介于1~50 m,取ξ = 5 mm,计算临界压力 $ {P}_{\mathrm{C}} $ 为:

    $$ \begin{array}{c}{P}_{\mathrm{C}}{\text{=}}\dfrac{2.59E{\xi }^{2}}{LD\sqrt{D/\xi }}{\text{=}}\dfrac{2.59·3\;160·{5}^{2}}{280·100·\sqrt{\dfrac{100}{5}}}{\text{=}}1.63\left({\rm{MPa}}\right)\end{array} $$ (5)

    安全系数为:

    $$ \begin{array}{c}n{\text{=}}\dfrac{{P}_{\rm{C}}}{p}{\text{=}}\dfrac{1.63}{0.5}{\text{=}}3.268\end{array} $$ (6)

    由于n = 3.268 > 2.5,故安全系数能满足水下50 m的设计要求。

    强度计算要求判断耐压舱所受外力分类,主要有长圆筒、短圆筒两种,由公式判断为:

    $$ \begin{array}{c}{L}_{\rm{C}}{\text{=}}1.17D\sqrt{D/\xi }{\text{=}}1.17·100·\sqrt{\dfrac{100}{5}}{\text{=}}523.224\;({\rm{mm}})\end{array} $$ (7)

    耐压舱的长度L = 270 (mm) <${L}_{\rm{C}}{\text{=}} 523.224 \; (\mathrm{m}\mathrm{m} )$,故判断为短圆筒。利用米塞斯 (Mises) 公式简化后,计算最小临界应力的公式为:

    $$ \begin{array}{c}{\sigma }_{c}{\text{=}}\dfrac{{P}_{\rm{C}}D}{2\xi }{\text{=}}\dfrac{1.3E}{\raisebox{1ex}{$L$} \left/ \raisebox{-1ex}{$D$}\right.}{\left(\dfrac{\xi }{D}\right)}^{1.5}{\text{=}}\\ \dfrac{1.3·3\;160}{\raisebox{1ex}{$280$} \left/ \raisebox{-1ex}{$100$}\right.}·{\left(\dfrac{5}{100}\right)}^{1.5}{\text{=}}16.403\;2\;({\rm{MPa}})\end{array} $$ (8)

    取安全系数n=2.5,则:

    $$ \begin{array}{c}{\sigma }_{n}{\text{=}}\dfrac{{\sigma }_{c}}{n}{\text{=}}\dfrac{16.403\;2}{2.5}{\text{=}}6.561\;28\end{array} $$ (9)

    耐压舱壳体的压力 ${\sigma }_{n}$为6.561 28 MPa,小于PMMA许用应力的值 (22.5 $ \mathrm{M}\mathrm{P}\mathrm{a} $),即耐压舱壳体 能满足50 m 的强度要求。

    在理论计算的基础上,利用Ansys有限元仿真软件对耐压舱壳体进行静应力仿真分析,求解后的等效应力、合位移、等效应变、安全系数云图等仿真结果如图12所示。

    图  12  耐压舱数值模拟
    Figure  12.  Numerical simulation of pressure-resistant chamber

    耐压舱外壳的最大总变形为0.108 73 mm,最大等效弹性应为0.001 51,属于微小变形,可忽略不计;最大等效应力为 $ 9.24\mathrm{ }\mathrm{M}\mathrm{P}\mathrm{a}\leqslant {\sigma }_{s}= 45\mathrm{ }\mathrm{M}\mathrm{P}\mathrm{a} $,满足PMMA材料的物理特性,在理论计算和静应力仿真上均满足设计要求,说明耐压舱舱体可以满足水深50 m的设计要求。

    对ROV试验测试包括水下性能测试和水质检测试验。在浙江海洋大学拖曳水池对ROV的自由度、稳定性、密封性以及上位机软件控制性能等进行试验,验证ROV设计的合理性与科学性。试验器材包括:试验水池,ROV整体,传感器组 (DO传感器、盐度传感器、pH传感器),传感器固定支架,采样桶,自吸泵,水管及浮力材料等 (图13)。

    图  13  试验场地和器材
    Figure  13.  Experimental site and equipment

    上位机软件调试包括ROV控制系统软件和多参数传感器组数据采集软件。ROV上位机软件采用开源的QGround Control (QGC),通过QGC调试ROV推进器推力、各类传感器作业参数、照明和观测系统、数据传输等功能,为水下试验提供作业前提。多参数传感器组数据采集软件由pH传感器、DO传感器、盐度传感器组成。调试方法为将传感器组调平后放到水中采集10 min,采集速率为1 s,通过软件观察曲线变化。两者调试试验见图14

    图  14  ROV控制系统软件和传感器数据采集软件调试
    Figure  14.  ROV control system software and sensor data acquisition software commissioning

    通过X-BOX控制ROV到达指定的深度,采用定深悬停命令在指定的位置悬停,通过自吸泵吸水输送到传感器组,对水质DO、温度、pH和盐度进行监测;在光线暗淡或水质较差的作业环境下,可调节ROV的照明系统,以满足水下高清视频、图片要求。现场抽水及观测试验见图15

    图  15  抽水及观测试验
    Figure  15.  Pumping and observation experiment

    ROV是一个复杂的系统,其中水下自由度为其重要组成部分。为进一步验证ROV在不同水深下正常作业的可行性,于水池中进行了潜水试验。测试中,ROV的目标下潜深度设定为4 m,通过X-BOX控制,验证ROV是否满足设计要求,ROV自由度测试见图16。经过试验验证,ROV的水下自由度能满足设计要求。

    图  16  ROV自由度测试
    Figure  16.  Freedom test of ROV degrees

    在本次测试中,ROV测试的下潜深度为0~4 m,并在池中保持定深状态,通过自吸泵吸入采样桶中,30 min内定点采集水池的水样,经传感器组检测,通过上位机软件实时监测水温、pH、DO和盐度4个水质指标随时间的变化情况。ROV在指定1 m位置悬停的4个参数变化见图17,在1、2、3 m的变化见图18

    图  17  ROV在水深1 m的温度、溶解氧、pH和盐度变化
    Figure  17.  Changes of temperature, DO, pH and salinity of ROV at water depth of 1 m
    图  18  ROV在不同水深的温度、溶解氧、pH和盐度变化
    Figure  18.  Changes of temperature, DO, pH and salinity of ROV at different water depths

    目前国内外针对水质检测的ROV设计与研究报道较少,水产养殖中的水质检测仍存在监测片面性、局部性、作业区域受限以及监测周期长等问题,尚未得到有效解决。基于国内外的研究现状,本研究设计了一种小型作业型ROV,用于监测水产养殖中的水质。样机的性能测试和水质参数测试结果显示,ROV的性能可满足工作要求;在水质参数测试中,ROV可以达到指定位置进行分层和定点水质采样。此外,还可以更换不同传感器,如氨氮传感器、氧化还原电位传感器(ORP)、浊度传感器等进行抽水测试,以满足不同的水质检测需要。

  • 图  1   大亚湾各分区

    Ⅰ.范和港湾;Ⅱ.大亚湾顶部;Ⅲ.港口列岛西部及哑铃湾;Ⅳ.大鹏澳;Ⅴ.大亚湾西南部;Ⅵ.大亚湾中东部;Ⅶ.大亚湾东南部;C1、C2.实测站点

    Figure  1.   Partition of Daya Bay

    Ⅰ.Fanhe gulf; Ⅱ.top of Daya Bay; Ⅲ.west of Gangkou islands and Yalin gulf; Ⅳ.Dapeng′ao; Ⅴ.southwest of Daya Bay; Ⅵ.east and middle of Daya Bay; Ⅶ.southeast of Daya Bay; C1, C2.survey stations

    图  2   大亚湾网格图

    Figure  2.   Model grids of Daya Bay

    图  3   潮位、流向、流速验证

    Figure  3.   Validation of tidal elevation and currents

    图  4   大亚湾欧拉余流图

    Figure  4.   Euler tidal residual currents in Daya Bay

    图  5   各时刻流出湾外质点的百分数

    Figure  5.   Percent of discharged particles flowing out the bay

    图  6   大亚湾保守物质输运扩散浓度分布

    Figure  6.   Concentration distributions of the passive tracer by the force of advection and diffusion

  • [1]

    DAVIES A. Meteorologically-induced circulation on the North-West European continental shelf: from a three-dimensional numerical model[J]. Oceanological Acta, 1982, 5(3): 269-280. https://www.semanticscholar.org/paper/Meteorologically-induced-circulation-on-the-from-a-Davies/6f8e3e2ddd0309326365a13585f5bf7a144edf2c

    [2]

    BACKHAUS J. Estimates of the variability of low frequency currents and flushing-times of the North Sea[P]. ICES Hydrography Committee, C.M. 1984, c: 24.

    [3] 匡国瑞, 杨殿荣, 喻祖祥, 等. 海湾水交换的研究——乳山东湾环境容量初步探讨[J]. 海洋环境科学, 1987, 6(1): 13-23.
    [4] 高抒, 谢钦春. 狭长形海湾与外海水体交换的一个物理模型[J]. 海洋通报, 1991, 10(3): 1-9. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsQ0hJTmV3UzIwMjQxMTA1MTcxMzA0Eg5RSzAwMDAwMTQ5NTI5NhoIczZkbWlib28%3D
    [5]

    TAKEOKA H. Water exchange in a time varying transport field[J]. J Oceanogr Soc Jap, 1987, 43(1): 21-27. doi: 10.1007/BF02110630

    [6]

    SIGNELL R, BUTMAN B. Modeling tidal exchange and dispersion in Boston Harbor[J]. J Geophysical Res, 1992, 97(c10): 15 591-15 606. doi: 10.1029/92JC01429

    [7] 袁柱瀚. 渤海海峡水交换的研究[D]. 青岛: 青岛海洋大学, 1997: 22. https://xueshu.baidu.com/usercenter/paper/show?paperid=5f6983a79ac239319afd8447cba8dfd5
    [8] 赵亮, 魏皓, 赵建中. 胶州湾水交换的数值研究[J]. 海洋与湖沼, 2002, 33(1): 23-29. doi: 10.11693/hyhz200201004004
    [9]

    ROGER L, POHLMANN T. Calculation of water exchange times in the ICES Boxes with a Eulerian model using a half-life time approach[J]. Dtsch Hydrogr Z, 1995, 47(4): 287-299. doi: 10.1007/BF02737789

    [10] 魏皓, 田怡, 周锋, 等. 渤海水交换的数值研究-水质模型对半交换时间的模拟[J]. 青岛海洋大学学报: 自然科学版, 2002, 32(4): 519-525. doi: 10.3969/j.issn.1672-5174.2002.04.003
    [11] 陈伟, 苏纪兰. 狭窄海湾潮交换的分段模式Ⅱ在象山港的应用[J]. 海洋环境科学, 1999, 18(3): 7-10. doi: 10.3969/j.issn.1007-6336.1999.03.002
    [12] 董礼先, 苏纪兰. 象山港水交换数值研究Ⅰ. 对流-扩散型的水交换模式[J]. 海洋与湖沼, 1999, 30(4): 410-415. https://www.semanticscholar.org/paper/%E8%B1%A1%E5%B1%B1%E6%B8%AF%E6%B0%B4%E4%BA%A4%E6%8D%A2%E6%95%B0%E5%80%BC%E7%A0%94%E7%A9%B6%EF%BC%9AI.%E5%AF%B9%E6%B5%81%E2%80%94%E6%89%A9%E6%95%A3%E5%9E%8B%E7%9A%84%E6%B0%B4%E4%BA%A4%E6%8D%A2%E6%A8%A1%E5%BC%8F-%E8%91%A3%E7%A4%BC%E5%85%88-%E8%8B%8F%E7%BA%AA%E5%85%B0/5994ea7a841ec72e203fb4fa64c55d693854e032
    [13] 董礼先, 苏纪兰. 象山港水交换数值研究Ⅱ. 模型的应用和水交换研究[J]. 海洋与湖沼, 1999, 30(5): 465-470. doi: 10.1088/0256-307X/15/12/024
    [14] 董礼先, 苏纪兰. 象山港盐度分布和水体混合Ⅰ. 盐度分布和环流结构[J]. 海洋与湖沼, 2000, 31(2): 151-158. doi: 10.3321/j.issn:0029-814X.2000.02.007
    [15] 董礼先, 苏纪兰. 象山港盐度分布和水体混合Ⅱ. 混合分析[J]. 海洋与湖沼, 2000, 31(3): 322-326. doi: 10.11693/hyhz200003014014
    [16] 中国海湾志编委员会. 中国海湾志第九分册[M]. 北京: 海洋出版社, 1993: 221-298.
    [17] 苏柯. 大亚湾核电液态排放的潮流污染扩散及海洋生态转移研究[D]. 广州: 暨南大学, 2005: 6-13. https://d.wanfangdata.com.cn/thesis/D464188
    [18] 刘云旭, 温伟英, 王文介. 大亚湾海域物理自净能力的时空差异性研究[J]. 热带海洋, 1999, 18(4): 61-68. http://www.irgrid.ac.cn/handle/1471x/2562021
    [19]

    BLUMBERG A F, MELLOR G L. A description of a three-dimensional coastal ocean circulation model in three-dimensional Coastal Models[M]. [S. l. ]: America Geophysical Union: Coastal and Estuarine Sciences, 1987: 1-161. doi: 10.1029/CO004p0001

    [20]

    CHEN C, ZHU J, RALPH E, et al. Prognostic modeling studies of the Keweenaw Current in Lake Superior: Part Ⅰ, formation and evolution[J]. J Phys Oceanogr, 2001, 31(2): 379-395. doi: 10.1175/1520-0485(2001)031<0379:PMSOTK>2.0.CO;2

    [21] 朱建荣, 朱首贤. ECOM模式的改进及在长江河口、杭州湾及邻近海区的应用[J]. 海洋与湖沼, 2003, 34(4): 364-388. doi: 10.3321/j.issn:0029-814X.2003.04.003
    [22] 杨国标. 大亚湾海区潮流运动特征[N]. 人民珠江, 2001(1): 30-32. 10.3969/j.issn.1001-9235.2001.01.010
    [23] 吴岩, 陶建华. 大亚湾潮流数学模型的校准及验证[J]. 港工技术, 1998(1): 6-11. https://xueshu.baidu.com/usercenter/paper/show?paperid=e9cd7b7b26aee768c47d569096f8d0cf&site=xueshu_se&hitarticle=1
图(6)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2008-01-13
  • 修回日期:  2008-04-06
  • 刊出日期:  2008-08-04

目录

/

返回文章
返回