Research Article  Open Access
Weiwei Si, Yifan Xue, Yanjun Liu, Zhitong Li, Gang Xue, "Energy Consumption Modeling and Sensitivity Analysis for DeepArgo Otarriinae Profiling Float", Mathematical Problems in Engineering, vol. 2020, Article ID 6035369, 14 pages, 2020. https://doi.org/10.1155/2020/6035369
Energy Consumption Modeling and Sensitivity Analysis for DeepArgo Otarriinae Profiling Float
Abstract
DeepArgo Otarriinae profiling float is a new type of Argo profiling float that has a maximum diving depth of more than 4,000 m. It can collect ocean scientific data allweather and uninterruptedly, which provides reliable data support for the global ocean scientific research. The working time of DeepArgo profiling float is an important indicator of its practicality and economy, and it is clear that the energy consumption is a key factor in determining its working time. In this paper, the single profile energy consumption model with 19 parameters of DeepArgo Otarriinae is established and the main effect indices and total effect indices of the energy consumption parameters to energy consumption are calculated using Sobol’ sensitivity analysis method, aiming to find the parameters that have the greatest impact on energy consumption. The results show that the gliding angle, the diving depth, and the gliding speed have a significant impact on energy consumption of DeepArgo Otarriinae. The results of simulation have a good match with the actual application and have certain reference significance for the determination of the design parameters and the selection of the navigation parameters. This paper also provides a new idea of multiparameter energy consumption modeling for underwater equipment using buoyancy regulation.
1. Introduction
Argo profiling float is a typical type of ocean observation platform, which provides continuous observations of ocean temperature and salinity versus pressure, from the sea surface to the maximum design depth of 2,000 m [1]. It can be equipped with sophisticated sensors and is capable of collecting scientific data under difficult and harsh conditions. It is the key equipment of Array for Realtime Geostrophic Oceanography, which is also named as the ocean observations with the global Argo array [2]. Compared with other ocean observation platforms, Argo profiling float has a simple structure and a long working time. It can periodically float to the sea surface to exchange data with the shore station via satellite. The temporal coverage of Argo data is much better than ship monitoring data [3]. Argo profiling float has become the primary source of nearrealtime data for ocean and atmospheric services and high quality data for climate research [4]. The data from Argo have been used to study climate variability in the Norwegian Sea and interpret the global oceanic heat transport [5, 6].
DeepArgo Otarriinae profiling float is a new type of Argo profiling float that has a maximum diving depth of more than 4,000 m while the conventional one has a maximum diving depth of 2,000 m. DeepArgo extends the observation range of Argo profiling float to the sea floor and accurately measures variation of seawater parameters in the deep sea areas [7, 8]. At present, four models of DeepArgo floats have been developed and are currently in use, including APEX profiling float and DeepSOLO float of the United States, DeepARVOR profiling float of France, and Deep NINJA profiling float of Japan [9–12]. China successfully developed the first set of Argo float prototypes in 2003 and completed the sea trial [13]. Just like other deep ocean exploration equipment, the performance of DeepArgo is also constrained by energy issues. The DeepARVOR profiling float developed by Ifremer extends the current operational depth down to 4,000m, and the results of the atsea experiments showed that the float could work for nine months at a rate of one cycle every two days in the West European Basin [14]. The energy supply of DeepArgo profiling float comes from the disposable or rechargeable batteries it carries. The battery capacity limits the operating time and navigation range of DeepArgo. Therefore, reducing the power consumption of the Argo profiling float is of great significance for improving the performance of DeepArgo. Lowpower technology has become one of the research trends of underwater vehicles. Woithe and Kremer described a new sensor management technique that makes use of published thermocline detection and tracking algorithms to simulate the triggering of sensors while a vehicle is within the thermocline instead of modifying its ﬂight proﬁle. This sensor management technique has also been successfully used in two atsea trials with Slocum Glider and the energy savings of up to 50% and 63%, respectively [15]. Yu established an applied mathematical model of SeaWing underwater glider that relates the gliding range to steady gliding motion parameters and energy consumption. In addition, Yu optimized the steady gliding motion parameters and developed a sensor scheduling strategy to increase the gliding range of SeaWing underwater glider [16]. Zhou presented a new ConstantTime Surfacing A^{∗} algorithmbased iterative path planning algorithm to address underwater glider with variable speeds which can maximize the usage of the favourable currents. From the simulation of hyperbolic currents field, the novel variable speed scheme can get a path with lesser energy cost compared with constant speed [17]. Kelasidi investigated the power consumption of different underwater robotic systems and compared the energy efficiency of different robots depending on the desired motion. In particular, Kelasidi compared the energy efﬁciency of underwater snake robots, which can provide both inspection and intervention capabilities and thus are interesting candidates for the nextgeneration inspection and intervention AUVs and the simulation results showed that the underwater snake robots are more energy efficient for all compared motion modes [18]. Lu proposed and validated a practical path planning algorithm considering both energy efficiency and estimation accuracy for the AUV. Besides, an optimal estimation strategy was developed based on Kalman Filter in order to minimize the ﬁlter uncertainty and energy consumption [19]. Yao proposed an energy consumption reduction method based on MPC (model predictive control) with a quadratic energy consumption term added into the cost function for AUV’s trajectory tracking control. The simulation results showed that the proposed MPC method is feasible and effective for energy consumption reduction of AUV even with great model uncertainty [20].
In this paper, the 4,000 m DeepArgo Otarriinae profiling float developed by Shandong University is regarded as the research object [21]. The shape and the main parameters are shown in Figure 1. Two sea trials with DeepArgo Otarriinae have been conducted successfully near the Mariana Trench since 2017. When the target working depth is 4000m, the number of working profiles of DeepArgo Otarriinae is more than 120. The remaining parts of this paper are organized as follows. In Section 2, the 4,000 m DeepArgo Otarriinae profiling float developed by Shandong University is introduced. In Section 3, the energy consumption model of DeepArgo Otarriinae is established. In Section 4, the parameters in the energy consumption model are determined. The sensitivity analysis of energy consumption parameters is conducted in Section 5. Finally, conclusions are provided in Section 6.
2. DeepArgo Otarriinae Profiling Float
2.1. Workflow of DeepArgo Otarriinae
The DeepArgo Otarriinae profiling float developed by Shandong University is mainly composed of pressure chamber, hydraulic system, communication and positioning unit, embedded control system, and sensors. The pressure chamber provides a sealed and waterproof environment for the internal subsystem. The hydraulic system adjusts the volume of DeepArgo Otarriinae by controlling the hydraulic oil to flow in or flow out the outer oil bladder to achieve buoyancy adjustment. The communication and positioning unit is mainly composed of a global positioning system, Iridium communication module, and wireless communication module. DeepArgo Otarriinae can switch the communication mode according to the distance between DeepArgo Otarriinae and the shore station and the mission demand. The embedded control system manages tasks in real time, performs status monitoring, and data analysis. The sensors are responsible for monitoring the operational status and data acquisition of the DeepArgo Otarriinae profiling float.
The general single profile workflow of DeepArgo Otarriinae is shown in Figure 2, which is mainly composed of communication and waiting stage on the sea surface, diving preparation stage, diving stage, drift stage, ascending preparation stage, and ascending stage [2]. However, in the drift stage, DeepArgo Otarriinae is in standby state and the standby current is very small, about microampere level. So the energy consumption of the drift stage can be simplified in the workflow from the perspective of energy consumption analysis. The simplified single profile workflow of DeepArgo Otarriinae is shown in Figure 3, which excludes drift stage. When DeepArgo Otarriinae is in the communication and waiting stage, the current position of DeepArgo Otarriinae is determined by GPS and the information is exchanged with the deck operator through Iridium or wireless communication. When DeepArgo Otarriinae receives the diving command, DeepArgo Otarriinae adjusts itself to the negative buoyancy through the hydraulic system during the diving preparation stage, and then DeepArgo Otarriinae starts to dive. When DeepArgo Otarriinae reaches the target depth, it adjusts itself to the positive buoyancy through the hydraulic system during the ascending preparation stage, thereby starting to ascend until reaching the water surface. At this point, a singleprofile workflow of DeepArgo Otarriinae ends and then DeepArgo Otarriinae enters the communication and waiting stage of the next workflow. Figure 4 shows the operating state of each subsystem in a single profile, in which the embedded control system always maintains the operating state, and other subsystems just operate when needed.
2.2. The Kinetic Model of DeepArgo Otarriinae
Before building an energy consumption model, we need to find out the relationship between the parameters of DeepArgo Otarriinae. The kinetic parameters and force analysis of DeepArgo Otarriinae during steady gliding are shown in Figure 5, in which denotes the inertial coordinate system; is the body coordinate system; is the buoyancy center of DeepArgo Otarriinae; is the attack angle; is the pitch angle; is the gliding angle; are the drag force and lift force of DeepArgo Otarriinae, respectively; is net buoyancy force; and is the gliding speed of DeepArgo Otarriinae in the body coordinate system.
The geometric relationships and force balance relationships of the steady gliding motion can be expressed as
The drag force and lift force of DeepArgo Otarriinae can be expressed aswhere and are drag coefficients , and are lift coefficients.
Applying the drag force and lift force obtained by equation (4) and (5) into equation (6) and (7), the following equation is obtained:
It can be seen from equation (8) that the attack angle is a function of the gliding angle and the hydrodynamic parameters when DeepArgo Otarriinae is in steady gliding state.
3. Energy Consumption Model of DeepArgo Otarriinae
Combining the single profile workflow of DeepArgo Otarriinae and the operating state of each subsystem, the energy consumption analysis of each subsystem is performed, respectively. As shown in Figure 3, the energy consumption of DeepArgo Otarriinae in a single profile mainly includes three parts: energy consumption of hydraulic system, energy consumption of positioning and communication unit, and the energy consumption of embedded control system and sensors.
3.1. Energy Consumption of Hydraulic System
The energy consumption of hydraulic system mainly includes the energy consumption of the servo motor and the solenoid valve. There are two buoyancy adjustment stages in a single profile workflow, which are diving preparation stage and ascending preparation stage. In the diving preparation stage, the hydraulic oil is driven to the inner oil bladder from outer oil bladder by the plunger pump on the sea surface, the volume of the outer oil bladder will shrink, and buoyancy force on the equipment will decrease, so DeepArgo Otarriinae starts to dive. Conversely, in the ascending preparation stage, the hydraulic oil is driven to the outer oil bladder from inner oil bladder by the plunger pump, the volume of the outer oil bladder will swell, and buoyancy force on the equipment will increase, so DeepArgo Otarriinae starts to ascend. The hydraulic schematic diagram of DeepArgo Otarriinae is shown in Figure 6.
Total energy consumption of the hydraulic system in a single profile can be described aswhere is the energy consumption of hydraulic system, is the power of the hydraulic system when DeepArgo Otarriinae is in the diving preparation stage, is the working time of the hydraulic system in the diving preparation stage, is the target working depth of DeepArgo Otarriinae, is the power of the hydraulic system in the ascending preparation stage, is the working time of the hydraulic system in the ascending preparation stage, is the volume change of the outer oil bladder, is the flow rate of the plunger pump when working in the diving preparation stage, and is the flow rate of the plunger pump in the ascending preparation stage.
Equation (9) can also be described as
From equations (1) and (2), we can get the following equation:
According to the relationship between the buoyancy and the volume change, the following equation can be obtained easily:where is the density of the sea water.
By substituting equations (6), (7), (13), and (14) into equation (12), total energy consumption of the hydraulic system in a single profile can also be described as
3.2. Energy Consumption of Communication and Positioning Unit
Since DeepArgo Otarriinae completes the communication and positioning process on the sea surface, the energy consumption in this process is only related to the power and running time of relevant electronic equipment and has no obvious relationship with the navigation parameters and design parameters of DeepArgo Otarriinae. Energy consumption generated by positioning and communication unit in a single profile can be described aswhere is energy consumption of communication and positioning unit, is the power of the electronic equipment used in communication and waiting stage, and is the time of communication and waiting stage.
3.3. Energy Consumption of Embedded Control System and the Sensors
DeepArgo Otarriinae profiling float can be equipped with a variety of sensors. Sensors can be divided into two types based on different operating modes: fullrange operating sensors and intermittent operating sensors. The sensors embedded in our prototype are shown in Table 1. In addition, our prototype also reserves three sensor interfaces, which can be embedded with different sensors according to different tasks.

The fullrange operating sensors mean that sensors operate throughout the single profile (excluding the communication and waiting stage). For example, the pressure transmitter monitors the change of pressure throughout the whole workflow. The CTD embedded in our prototype is SBE 37SMMicroCAT CT(P) Recorder which is selfcontained. The sampling interval of CTD is 6 s to 6 h, and the sampling time is 1.8 s to 2.6 s. The energy carried by CTD can collect 960,000 data. The energy consumed by CTD is only used to read the data collected by CTD and not to drive the pump and electronic components inside CTD. According to the specifications of the Argo organization, CTD is not open all the way, but only during ascending stage. In the simplified workflow of DeepArgo Otarriinae, the buoyancy adjustment and hydrodynamic coefficient of DeepArgo Otarriinae are basically the same in the diving stage and ascending stage. In other words, the values of the factors that affect the speed of ascent and descent are similar, so the speed of ascent and descent is basically the same and the time of the diving stage and ascending stage is almost identical. In order to facilitate the calculation, we include CTD in the fullrange operating sensors. The energy consumption of the fullrange operating sensors can be described aswhere is the energy consumption of the fullrange operating sensors, is the sum of the average power of the fullrange operating sensors excluding CTD, is the power of CTD, is the time of the diving stage or ascending stage, and is the vertical component of the gliding speed.
The intermittent operating sensors mean that sensors operate at regular depth intervals in order to save energy or to meet measurement requirements. The specific sampling process of intermittent operating altimeter is shown in Figure 7, and the energy consumption of the intermittent operating sensors in the single profile can be described aswhere is the energy consumption of the intermittent operating sensors, is the number of intermittent operating sensors carried by DeepArgo Otarriinae, is the time of single operation of sensor, is the height interval at which the sensor is turned on, and is the power of single operation of sensor.
The embedded control system keeps working during the entire profile, including the communication and waiting stage. The energy consumption of this part can be described aswhere is the energy consumption of the embedded control system, and is the power of the embedded control system.
In summary, the total energy consumption of control system and sensors can be described aswhere is the energy consumption of embedded control system and sensors.
3.4. The Energy Consumption Model in a Single Profile
The total energy consumption model in a single profile can be obtained by adding the energy consumption of each part mentioned above, which can be described as
The energy consumption model of DeepArgo Otarriinae in a single profile has the following characteristics:(1)There are many parameters involved, including hydrodynamic parameters, hydraulic system parameters, navigation parameters, main design parameters, and other 19 parameters.(2)There is a coupling relationship between the energy consumption of each subsystem. The energy consumption of each subsystem in the model is related to multiple parameters, and changes in some parameters will lead to changes of multiple subsystems in energy consumption.(3)The energy consumption model is highly nonlinear.
4. Determination of Energy Consumption Parameters
Some parameters are not determined in the design of DeepArgo Otarriinae. Their values will show a certain change law with the change of working environment, so they need to be obtained through experiment or simulation. Another part of the parameters are determined in the design and their parameter values are shown in Table 2.

4.1. Determination of Hydraulic System Parameters
As shown in Figure 8, in order to obtain the change law of the power and flow of hydraulic system during the diving preparation stage and the ascending preparation stage, we used a hydraulic booster system to simulate the highpressure underwater environment and tested the operation of the hydraulic system under different pressures.
Through the test of DeepArgo Otarriinae hydraulic system, we finally obtained the power of hydraulic system in the diving preparation stage and ascending preparation stage as shown in Figure 9.
(a)
(b)
It can be seen from Figure 9(a) that the power of hydraulic system is about 66W when it is stable in the diving preparation stage. The power of the hydraulic system in the ascending preparation stage at different target depth is fitted by curve, which is shown in Figure 9(b). Fitted curve equation is as follows:
According to the fitting results, the complex correlation coefficient of the fitting curve is 0.9678, which fully meets the accuracy requirement.
The flow of the hydraulic system in the test is shown as Figure 10.
(a)
(b)
From Figure 10(a), we can calculate that the value of flow of hydraulic system is about 0.119 L/min in the diving preparation stage. It can be seen from Figure 10(b) that under different working pressures, the flow of the hydraulic system is stabilized at 0.06 L/min in the ascending preparation stage.
The power and flow of hydraulic system is determined by the performance and design parameters of the system, which can be considered as a fixed value and sensitivity analysis is not required.
4.2. Determination of Hydrodynamic Parameters
There are several methods that will produce results for hydrodynamic parameters based on a given geometry. The methods include analytical, experimental, computational, and semiempirical approaches [22]. In this paper, computational fluid dynamics (CFD) is used to solve the viscous hydrodynamic parameters, and FLUENT is used to calculate the resistance, lift force, and torque of DeepArgo Otarriinae in steady gliding when the angle of attack is 0°, 2°, and 4° and gliding speed is 0.1 m/s, 0.3 m/s, and 0.5 m/s, respectively. The hydrodynamic parameters of DeepArgo Otarriinae in steady gliding are calculated using the least square method, as shown in Table 3.

The hydrodynamic parameters are necessary and important to measure the hydrodynamic performance of DeepArgo Otarriinae. The variation of the hydrodynamic parameters has a very important impact on the energy consumption and voyage of DeepArgo Otarriinae, so sensitivity analysis for hydrodynamic parameters is needed. For the convenience of analysis, we assume that each hydrodynamic coefficient can fluctuate by 20% compared with the prototype and follows the uniform distribution within the respective value range.
4.3. Determination of Navigation Parameters
The navigation parameters are arguments that the operator needs to set directly or indirectly and are also the key to measure the underwater operation state of DeepArgo Otarriinae. Therefore, it is necessary to study the influence of navigation parameters on the energy consumption and voyage of DeepArgo Otarriinae.(1)Gliding speed The gliding speed of the DeepSOLO float is about 0.06 m/s, the gliding speed of the DeepSOLO 2 float is about 0.12 m/s, and the gliding speed of the ArvorC profiling float is about 0.15–0.20 m/s [23, 24]. Therefore, we set the gliding speed of DeepArgo Otarriinae as 0.05–0.20 m/s when it moves steadily, and it follows uniform distribution within this interval.(2)Gliding angle According to equation (1) and equation (2), we can get the equation as follows: By substituting equations (4) and (5) into equation (33), we can get equation as follows: The attack angle α is a certain value when DeepArgo Otarriinae moves steadily, that is to say, equation (34) has a solution and We stipulate that the gliding angle of DeepArgo Otarriinae is within , but in order to ensure the economicality of navigation, generally the gliding angle is within . Therefore, the gliding angle needs to meet the following equation for keeping stable motion: Substituting the hydrodynamic parameters shown in Table 3 into equation (26), we finally determine that the gliding angle of DeepArgo Otarriinae is within . Under the condition that the velocity of ocean current is almost constant, the gliding angle of DeepArgo Otarriinae is set uniformly distributed in the above range.(3)Diving depth The diving depth of DeepArgo Otarriinae is determined by the operator. The diving depth of DeepArgo Otarriinae ranges from 0 to 4,000 m, and can be considered to be uniformly distributed within this range.(4)Communication and waiting time Since the satellite signal intensity will be affected by weather condition, marine environment, satellite operating state, and other factors, and the operator may perform additional data exchange operations with DeepArgo Otarriinae according to the task requirements, the communication and waiting time will have a certain randomness. As a matter of experience, positioning and communication time is within 600–1200 s, and can be considered to be uniformly distributed within this range.(5)Intermittent operating sensors opening depth interval and running time The intermittent operating sensor carried by DeepArgo Otarriinae is only an altimeter that prevents DeepArgo Otarriinae from bottoming out. Too small opening depth interval of altimeter will greatly increase energy consumption while too large opening depth interval will easily cause DeepArgo Otarriinae bottom out. Considering the diving speed of DeepArgo Otarriinae and the buoyancy adjustment time, we set the altimeter opening depth interval to be uniformly distributed between 0m and 50m. Similar to the altimeter opening depth interval, too short running time of the altimeter may cause DeepArgo Otarriinae to bottom out while too long running time will increase the energy consumption. Considering the diving speed of DeepArgo Otarriinae and the installation position of the altimeter, we set the running time to be uniformly distributed between 2 and 10 s.
5. Sensitivity Analysis of Energy Consumption Parameters
5.1. Sensitivity Analysis Based on Sobol’ Method
Sensitivity analysis (SA) aims to identify the key parameters that affect model performance, and it plays important roles in model parameterization, calibration, optimization, and uncertainty quantification [25]. Sensitivity analysis includes local sensitivity analysis (LSA) and global sensitivity analysis (GSA), where global sensitivity analysis is a common technique for evaluating the relative contribution of input parameters to the overall output by considering the change of input factors throughout the change space. Sobol’ method is based on the variance decomposition first proposed by Sobol’ [26], which is suitable for nonlinearity having many input variables and strong stability [27]. In recent years, it has been widely used in many fields including marine engineering [28–30].
Sobol’ method estimates the importance of each input parameter by means of variance decomposition [31]. According to Sobol’ method, the energy consumption model can be decomposed as follows:where is a function expression of the energy consumption model, are M independent input parameters affecting energy consumption, respectively.
Assuming that the square of is integrable, the following equation can be obtained by squaring the left and right sides of equation (30) and then multiple integrals of input variables in space.
The left part of equation (30) is called the total variance, which can be described aswhere is total variance.
Each term on the right side of equation (30) is called partial variance. The single effect of input parameter on the model output can be described aswhere is the single effect of input parameter on the model output.
The joint effect of input parameters on the model output can be described aswhere is the joint effect of input parameters on the model output.
The main effect indices and the total effect indices of the parameter to the model can be described aswhere is the effect of all parameters but .
The difference between the main effect and the total effect indices of the parameter is a measure of the interaction of with all other parameters. The main effect indices reflects effect on the model output when the parameter changes individually while the total effect indices reflects the effect on the model output when the parameter interacts with other parameters. In addition, if the total effect indices of a parameter is close to 0, we can assume that this parameter has no effect on the model output. GSA including the Sobol’ method requires a large number of samples to ensure convergence of the calculation results. In order to get maximisation of the minimum interpoint Euclidean distance between any two sampled points, we used Latin hypercube sampling (LHS) method to design the sample space and the sample space size is . The total number of model evaluations for the approximation of both the main and total effects indices is therefore equal towhere is the total number of model evaluations, n is the number of samples, and is the number of independent parameters obtained per sample.
5.2. Results of Simulation
Since the energy model established in this paper is complex and some parameters are constant, it is impossible and unnecessary to perform sensitivity analysis on each parameter. Therefore, only the parameters most concerned by the researchers are selected for sensitivity analysis. In Section 4, the range and distribution of the 10 parameters for sensitivity analysis have been determined, as shown in Table 4.

The number of independent parameters is 10, and the effect indices of each parameter can be guaranteed to converge when the number of samples is 9600. Therefore, the total number of model evaluations is 115200 according to equation (36). The simulation was conducted on a computer configured with 2.80 GHz i77700HQ processor. The software we used was MATLAB 2019a with SAFE toolbox [32], which took a total time of 150s. We obtain the main effect indices and the total effect indices of the above 10 parameters to the energy consumption model, which is shown in Figures 11 and 12, respectively.
In Figures 11 and 12, the x axis is the sampling times and y axis is the normalized values of main effect or total effect. Besides, the values of main effect or total effect of each parameter tends to be stable and converged with the increase of sampling times, which means our sampling times is suitable and the result of the simulation is convinced. We can see that the sensitivity indices of hydrodynamic parameters and sensor sampling parameters are relatively close and the values are small. The results of Figures 11 and 12 are averaged and ranked by the magnitude of the total sensitivity indices as shown in Figure 13 and Table 5.

5.3. Analysis and Discussion
Analysis and discussion about the above simulation results are given as follows:(1)Within the range of parameters studied in this paper, the gliding angle , the diving depth , and the gliding speed are at a fairly high level and have a significant impact on the energy consumption of DeepArgo Otarriinae in a single profile. The higher values of , , and will increase energy consumption. In fact, these three parameters determine the key parameters such as input power of the hydraulic system and the running time of DeepArgo Otarriinae, which have a good match with the actual application. The main effect indices and total effect indices of other energy consumption parameters are much smaller, indicating that their changes have limited impact on energy consumption.(2)Among the three key parameters that affect the energy consumption, the diving depth is completely set by operators according to the requirements of the working task. There are no posture adjust module and speed adjust module in our prototype. The gliding angle and gliding speed mainly depend on the volume, mass, and shape of DeepArgo Otarriinae which have been determined during the design process. According to Table 5, we can see clearly that smaller hydrodynamic coefficients , , , and will reduce the energy consumption of DeepArgo Otarriinae. But smaller hydrodynamic coefficients will lead to a larger gliding speed which will increase energy consumption greater. Obtaining the optimal solution for hydrodynamic coefficients and gliding speed is a direction for future work.(3)In order to reduce the energy consumption of DeepArgo Otarriinae, for designers, the simulation results above point out the optimization direction in the design process. Designers should focus on reducing the glide angle and the gliding speed through a large number of prototype pool experiments and hydrodynamic simulation. For operators, the simulation results show the influence of diving depth on energy consumption mathematically. When there is no mandatory requirement for the diving depth, operators should select the diving depth reasonably.
6. Conclusions
In this paper, a new modeling method for multiparameter energy consumption model of DeepArgo Otarriinae was proposed. First, the kinematics model of DeepArgo Otarriinae was established to show the relationship between various parameters. Then, the energy consumption of each subsystem of DeepArgo Otarriinae in a single profile was calculated. Finally, an energy consumption model including hydraulic system parameters, navigation parameters, and other 19 parameters was established, which could describe the energy consumption of DeepArgo Otarriinae in a single profile accurately and comprehensively. The energy consumption model we proposed contains many parameters and is highly nonlinear, which has certain research value and provides a new idea of multiparameter energy consumption modeling for underwater equipment using buoyancy regulation because their energy consumption components are similar. To reduce energy consumption of DeepArgo Otarriinae and extend its working time, the Sobol’ method was introduced to analyze the sensitivity of energy consumption parameters and the main effect indices and total effect indices of the selected parameters to energy consumption were calculated. The results of sensitivity analysis show that the gliding angle , the diving depth , and the glide speed have a greater impact on energy consumption. Therefore, reducing the gliding angle and the diving depth in the design and setting the diving depth reasonably in operation can reduce the energy consumption of DeepArgo Otarriinae effectively.
Future work can be done by establishing a more accurate kinematics model which incorporates more parameters and obtaining the optimal solution for hydrodynamic coefficients and gliding speed in the design. Furthermore, we can study the relationship between energy consumption and diving depth or gliding speed to lay the foundation for further research on reducing energy consumption of DeepArgo Otarriinae.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
Weiwei SI and Yifan XUE contributed equally to this article.
Acknowledgments
This paper was supported by the Key Research and Development Program of Shandong Province (NO. 2019GHY112077), the Guangdong Basic and Applied Basic Research Foundation (NO. 2019A1515110040), the Research Fund of State Key Laboratory of Ocean Engineering (Shanghai Jiaotong University) (NO. 1911), the Shandong Province Postdoctoral Innovation Project (NO. 201901001), the Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology (NO. MGQNLM201806), the Shenzhen Science and Technology R&D Foundation (NO. JCYJ20180305164217766), and the Qingdao Postdoctoral Applied Research Project. Gang Xue is the visiting research fellow of State Key Laboratory of Ocean Engineering. Shenzhen Research Institute of Shandong University is the primary supported department.
References
 D. Roemmich, M. H. Alford, H. Claustre et al., “On the future of argo: a global, fulldepth, multidisciplinary array,” Frontiers in Marine Science, vol. 6, 2019. View at: Publisher Site  Google Scholar
 S. C. Riser, H. J. Freeland, D. Roemmich et al., “Fifteen years of ocean observations with the global Argo array,” Nature Climate Change, vol. 6, no. 2, pp. 145–153, 2016. View at: Publisher Site  Google Scholar
 N. Haavisto, L. Tuomi, P. Roiha, S. M. Siiriä, P. Alenius, and T. Purokoski, “Argo floats as a novel part of the monitoring the hydrography of the bothnian sea,” Frontiers in Marine Science, vol. 5, 2018. View at: Publisher Site  Google Scholar
 M. Kuusela and M. L. Stein, “Locally stationary spatiotemporal interpolation of Argo profiling float data,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 474, no. 2220, p. 20180400, 2018. View at: Publisher Site  Google Scholar
 K. A. Mork, Øk Skagseth, and H. Søiland, “Recent warming and freshening of the Norwegian sea observed by argo data,” Journal of Climate, vol. 32, no. 12, pp. 3695–3705, 2019. View at: Publisher Site  Google Scholar
 B. Sun, C. Liu, and F. Wang, “Global meridional eddy heat transport inferred from Argo and altimetry observations,” Scientific Reports, vol. 9, no. 1, 2019. View at: Publisher Site  Google Scholar
 T. Wang, S. T. Gille, M. R. Mazloff, N. V. Zilberman, and Y. Du, “Numerical simulations to project argo float positions in the middepth and deep southwest pacific,” Journal of Atmospheric and Oceanic Technology, vol. 35, no. 7, pp. 1425–1440, 2018. View at: Publisher Site  Google Scholar
 G. C. Johnson, S. G. Purkey, N. V. Zilberman, and D. Roemmich, “Deep Argo quantifies bottom water warming rates in the Southwest Pacific Basin,” Geophysical Research Letters, vol. 46, no. 5, pp. 2662–2669, 2019. View at: Google Scholar
 K. S. Johnson, L. Coletti, H. Jannasch, T. Martz, D. Swift, and S. Riser, “Longterm observations of ocean biogeochemistry with nitrate and oxygen sensors in apex profiling floats,” in Proceedings of the Agu Fall Meeting, vol. 89, San Francisco, CA, USA, December 2008. View at: Google Scholar
 D. Roemmich, J. T. Sherman, R. E. Davis, K. Grindley, M. McClune, and C. J. Parker, “Deep solo: a fulldepth profiling float for the argo Program,” Journal of Atmospheric and Oceanic Technology, vol. 36, no. 10, pp. 1967–1981, 2019. View at: Publisher Site  Google Scholar
 X. André, B. Moreau, and S. Le Reste, “Argos3 satellite communication system: implementation on the arvor oceanographic profiling floats,” Journal of Atmospheric and Oceanic Technology, vol. 32, no. 10, pp. 1902–1914, 2015. View at: Publisher Site  Google Scholar
 T. Kobayashi, K. Amaike, K. Watanabe et al., “Deep NINJA: a new float for deep ocean observation developed in Japan,” in Proceedings of the IEEE Symposium on Underwater Technology and Workshop on Scientific Use of Submarine Cables and Related Technologies, IEEE, Tokyo, Japan, April 2011. View at: Google Scholar
 Z. Liu, X. Wu, J. Xu et al., “China Argo project: progress in China Argo ocean observations and data applications,” Acta Oceanologica Sinica, vol. 36, no. 6, pp. 1–11, 2017. View at: Publisher Site  Google Scholar
 S. Le Reste, V. Dutreuil, X. André et al., “DeepArvor: a new profiling float to extend the argo observations down to 4000m depth,” Journal of Atmospheric and Oceanic Technology, vol. 33, no. 5, pp. 1039–1055, 2016. View at: Publisher Site  Google Scholar
 H. C. Woithe and U. Kremer, “Feature based adaptive energy management of sensors on autonomous underwater vehicles,” Ocean Engineering, vol. 97, pp. 21–29, 2015. View at: Publisher Site  Google Scholar
 J. Yu, F. Zhang, A. Zhang, W. Jin, and Y. Tian, “Motion parameter optimization and sensor scheduling for the seawing underwater glider,” IEEE Journal of Oceanic Engineering, vol. 38, no. 2, pp. 243–254, 2013. View at: Publisher Site  Google Scholar
 Y. Zhou, J. Yu, and X. Wang, “Path planning method of underwater glider based on energy consumption model in current environment,” Lecture Notes in Computer Science, vol. 142–152, 2014. View at: Google Scholar
 E. Kelasidi, K. Y. Pettersen, and J. T. Gravdahl, “Energy efficiency of underwater robots,” IFACPapersOnLine, vol. 48, no. 16, pp. 152–159, 2015. View at: Publisher Site  Google Scholar
 D. Lu, R. Cui, and P. Wang, “Energy efficient path planning of autonomous underwater vehicles for environment modelling,” in Proceedings of the International Conference on Multisensor Fusion and Information Integration for Intelligent Systems, pp. 1–6, IEEE, Beijing, China, September 2014. View at: Google Scholar
 F. Yao, C. Yang, M. Zhang, and Y. Wang, “Optimization of the energy consumption of depth tracking control based on model predictive control for autonomous underwater vehicles,” Sensors, vol. 19, no. 1, p. 162, 2019. View at: Publisher Site  Google Scholar
 Z. Li, Y. Liu, F. Guo, G. Xue, S. Li, and X. Li, “Multiobjective optimization of the shell in autonomous intelligent argo profiling float,” Ocean Engineering, vol. 187, p. 106176, 2019. View at: Publisher Site  Google Scholar
 J. S. Geisbert, Hydrodynamic Modeling for Autonomous Underwater Vehicles Using Computational and Semiempirical Methods, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, 2007, https://hdl.handle.net/10919/33195, Master’s Thesis.
 N. Zilberman and D. Roemmich, “SIO float lab,” 2016, http://www.argo.ucsd.edu/Deep_SOLO_AST16.pdf. View at: Google Scholar
 S. Soylu, B. J. Buckham, and R. P. Podhorodeski, “Redundancy resolution for underwater mobile manipulators,” Ocean Engineering, vol. 37, no. 23, pp. 325–343, 2010. View at: Publisher Site  Google Scholar
 X. Song, J. Zhang, C. Zhan, Y. Xuan, M. Ye, and C. Xu, “Global sensitivity analysis in hydrological modeling: review of concepts, methods, theoretical framework, and applications,” Journal of Hydrology, vol. 523, pp. 739–757, 2015. View at: Publisher Site  Google Scholar
 I. M. Sobol, “On sensitivity estimation for nonlinear mathematical models,” Institute for Mathematical Modelling, vol. 2, no. 1, pp. 112–118, 1990. View at: Google Scholar
 F. Pianosi, K. Beven, J. Freer et al., “Sensitivity analysis of environmental models: a systematic review with practical workflow,” Environmental Modelling & Software, vol. 79, pp. 214–232, 2016. View at: Publisher Site  Google Scholar
 L. J. M. Peeters, G. M. Podger, T. Smith, T. Pickett, R. H. Bark, and S. M. Cuddy, “Robust global sensitivity analysis of a river management model to assess nonlinear and interaction effects,” Hydrology and Earth System Sciences, vol. 18, no. 9, pp. 3777–3785, 2014. View at: Publisher Site  Google Scholar
 K. Menberg, Y. Heo, and R. Choudhary, “Sensitivity analysis methods for building energy models: comparing computational costs and extractable information,” Energy and Buildings, vol. 133, pp. 433–445, 2016. View at: Publisher Site  Google Scholar
 X. Cheng, G. Li, R. Skulstad et al., “Datadriven uncertainty and sensitivity analysis for ship motion modeling in offshore operations,” Ocean Engineering, vol. 179, pp. 261–272, 2019. View at: Publisher Site  Google Scholar
 A. Saltelli and I. M. Sobol, “Sensitivity analysis for nonlinear mathematical models: numerical experience,” Institute for Mathematical Modelling, vol. 7, no. 11, pp. 16–28, 1995, http://mi.mathnet.ru/eng/mm1810. View at: Google Scholar
 F. Pianosi, F. Sarrazin, and T. Wagener, “A matlab toolbox for global sensitivity analysis,” Environmental Modelling & Software, vol. 70, pp. 80–85, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Weiwei Si et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.