Energy Consumption Modeling and Sensitivity Analysis for Deep-Argo Otarriinae Profiling Float

Deep-Argo Otarriinae profiling float is a new type of Argo profiling float that has a maximum diving depth of more than 4,000m. It can collect ocean scientific data all-weather and uninterruptedly, which provides reliable data support for the global ocean scientific research.2e working time of Deep-Argo 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 Deep-Argo 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.2e results show that the gliding angle, the diving depth, and the gliding speed have a significant impact on energy consumption of Deep-Argo Otarriinae. 2e 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. 2is paper also provides a new idea of multiparameter energy consumption modeling for underwater equipment using buoyancy regulation.


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 Real-time 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. e temporal coverage of Argo data is much better than ship monitoring data [3]. Argo profiling float has become the primary source of near-real-time data for ocean and atmospheric services and high quality data for climate research [4]. e data from Argo have been used to study climate variability in the Norwegian Sea and interpret the global oceanic heat transport [5,6].
Deep-Argo 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. Deep-Argo 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 Deep-Argo floats have been developed and are currently in use, including APEX profiling float and Deep-SOLO float of the United States, Deep-ARVOR profiling float of France, and Deep NINJA profiling float of Japan [9][10][11][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 Deep-Argo is also constrained by energy issues. e Deep-ARVOR 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]. e energy supply of Deep-Argo profiling float comes from the disposable or rechargeable batteries it carries. e battery capacity limits the operating time and navigation range of Deep-Argo.
erefore, reducing the power consumption of the Argo profiling float is of great significance for improving the performance of Deep-Argo. Low-power 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 flight profile.
is sensor management technique has also been successfully used in two at-sea trials with Slocum Glider and the energy savings of up to 50% and 63%, respectively [15]. Yu established an applied mathematical model of Sea-Wing 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 Sea-Wing underwater glider [16]. Zhou presented a new Constant-Time 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 efficiency 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 filter 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. e 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 Deep-Argo Otarriinae profiling float developed by Shandong University is regarded as the research object [21]. e shape and the main parameters are shown in Figure 1. Two sea trials with Deep-Argo Otarriinae have been conducted successfully near the Mariana Trench since 2017. When the target working depth is 4000m, the number of working profiles of Deep-Argo Otarriinae is more than 120. e remaining parts of this paper are organized as follows. In Section 2, the 4,000 m Deep-Argo Otarriinae profiling float developed by Shandong University is introduced. In Section 3, the energy consumption model of Deep-Argo Otarriinae is established. In Section 4, the parameters in the energy consumption model are determined. e sensitivity analysis of energy consumption parameters is conducted in Section 5. Finally, conclusions are provided in Section 6.

Workflow of Deep-Argo Otarriinae.
e Deep-Argo Otarriinae profiling float developed by Shandong University is mainly composed of pressure chamber, hydraulic system, communication and positioning unit, embedded control system, and sensors. e pressure chamber provides a sealed and waterproof environment for the internal subsystem. e hydraulic system adjusts the volume of Deep-Argo Otarriinae by controlling the hydraulic oil to flow in or flow out the outer oil bladder to achieve buoyancy adjustment. e communication and positioning unit is mainly composed of a global positioning system, Iridium communication module, and wireless communication module. Deep-Argo Otarriinae can switch the communication mode according to the distance between Deep-Argo Otarriinae and the shore station and the mission demand. e embedded control system manages tasks in real time, performs status monitoring, and data analysis. e sensors are responsible for monitoring the operational status and data acquisition of the Deep-Argo Otarriinae profiling float. e general single profile workflow of Deep-Argo 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, Deep-Argo 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. e simplified single profile workflow of Deep-Argo Otarriinae is shown in Figure 3, which excludes drift stage. When Deep-Argo Otarriinae is in the communication and waiting stage, the current position of Deep-Argo Otarriinae is determined by GPS and the information is exchanged with the deck operator through Iridium or wireless communication. When Deep-Argo Otarriinae receives the diving command, Deep-Argo Otarriinae adjusts itself to the negative buoyancy through the hydraulic system during the diving preparation stage, and then Deep-Argo Otarriinae starts to dive. When Deep-Argo 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 single-profile workflow of Deep-Argo Otarriinae ends and then Deep- Argo 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.

e Kinetic Model of Deep-Argo Otarriinae.
Before building an energy consumption model, we need to find out the relationship between the parameters of Deep-Argo Otarriinae.
e kinetic parameters and force analysis of Deep-Argo Otarriinae during steady gliding are shown in Figure 5, in which O − XZ denotes the inertial coordinate system; B − e 1 e 3 is the body coordinate system; B is the buoyancy center of Deep-Argo Otarriinae; α is the attack angle; θ is the pitch angle; ξ is the gliding angle; D, L are the drag force and lift force of Deep-Argo Otarriinae, respectively; mg is net buoyancy force; and v is the gliding speed of Deep-Argo Otarriinae in the body coordinate system. e geometric relationships and force balance relationships of the steady gliding motion can be expressed as e drag force and lift force of Deep-Argo Otarriinae can be expressed as where K D0 and K D are drag coefficients , K L0 and K L 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 Deep-Argo Otarriinae is in steady gliding state.

Energy Consumption Model of Deep-Argo Otarriinae
Combining the single profile workflow of Deep-Argo 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 Deep-Argo 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.

Energy Consumption of Hydraulic
System. e energy consumption of hydraulic system mainly includes the energy consumption of the servo motor and the solenoid valve.
ere 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 Deep-Argo 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 Deep-Argo Otarriinae starts to ascend. e hydraulic schematic diagram of Deep-Argo Otarriinae is shown in Figure 6.
Total energy consumption of the hydraulic system in a single profile can be described as where Q n is the energy consumption of hydraulic system, P hydra0 is the power of the hydraulic system when Deep-Argo Otarriinae is in the diving preparation stage, t 1 is the working time of the hydraulic system in the diving preparation stage, H is the target working depth of Deep-Argo Otarriinae, P hydra (H) is the power of the hydraulic system in the ascending preparation stage, t 2 is the working time of the hydraulic system in the ascending preparation stage, ΔV is the volume change of the outer oil bladder, q pump0 is the flow rate of the plunger pump when working in the diving preparation stage, and q pumpH 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 ρ water 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

Energy Consumption of Communication and Positioning
Unit. Since Deep-Argo 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 Deep-Argo Otarriinae. Energy consumption generated by positioning and communication unit in a single profile can be described as where Q f is energy consumption of communication and positioning unit, P f is the power of the electronic equipment used in communication and waiting stage, and t f is the time of communication and waiting stage.

Energy Consumption of Embedded Control System and the Sensors.
Deep-Argo 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. e 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.  Mathematical Problems in Engineering e full-range 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. e CTD embedded in our prototype is SBE 37-SMMicroCAT CT(P) Recorder which is self-contained. e sampling interval of CTD is 6 s to 6 h, and the sampling time is 1.8 s to 2.6 s. e energy carried by CTD can collect 960,000 data. e 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 Deep-Argo Otarriinae, the buoyancy adjustment and hydrodynamic coefficient of Deep-Argo 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 full-range operating sensors. e energy consumption of the full-range operating sensors can be described as where Q t1 is the energy consumption of the full-range operating sensors, P t1 is the sum of the average power of the full-range operating sensors excluding CTD, P CTD is the power of CTD, t 3 is the time of the diving stage or ascending stage, and v Z is the vertical component of the gliding speed. e intermittent operating sensors mean that sensors operate at regular depth intervals in order to save energy or to meet measurement requirements. e 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 as where Q t2 is the energy consumption of the intermittent operating sensors, j is the number of intermittent operating sensors carried by Deep-Argo Otarriinae, Δt i is the time of single operation of i − th sensor, Δh i is the height interval at which the sensor is turned on, and P si is the power of single operation of i − th sensor. e embedded control system keeps working during the entire profile, including the communication and waiting stage. e energy consumption of this part can be described as where Q t3 is the energy consumption of the embedded control system, and P t3 is the power of the embedded control system. In summary, the total energy consumption of control system and sensors can be described as    6 Mathematical Problems in Engineering where Q t is the energy consumption of embedded control system and sensors.

e Energy Consumption Model in a Single
Profile. e 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 e energy consumption model of Deep-Argo Otarriinae in a single profile has the following characteristics: (1) ere are many parameters involved, including hydrodynamic parameters, hydraulic system parameters, navigation parameters, main design parameters, and other 19 parameters. (2) ere is a coupling relationship between the energy consumption of each subsystem. e 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) e energy consumption model is highly nonlinear.

Determination of Energy Consumption Parameters
Some parameters are not determined in the design of Deep-Argo Otarriinae. eir 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.

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. rough the test of Deep-Argo 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.
It can be seen from Figure 9(a) that the power of hydraulic system P hydra0 is about 66W when it is stable in the diving preparation stage. e 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 R 2 of the fitting curve is 0.9678, which fully meets the accuracy requirement. e flow of the hydraulic system in the test is shown as Figure 10.
From Figure 10(a), we can calculate that the value of flow q pump0 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 q pumpH of the hydraulic system is stabilized at 0.06 L/min in the ascending preparation stage. e 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.

Determination of Hydrodynamic Parameters.
ere are several methods that will produce results for hydrodynamic parameters based on a given geometry. e 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 Deep-Argo 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. e hydrodynamic parameters of Deep-Argo Otarriinae in steady gliding are calculated using the least square method, as shown in Table 3. e hydrodynamic parameters are necessary and important to measure the hydrodynamic performance of Deep-Argo Otarriinae. e variation of the hydrodynamic parameters has a very important impact on the energy Mathematical Problems in Engineering 7 consumption and voyage of Deep-Argo 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.

Determination of Navigation
Parameters. e 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 Deep-Argo Otarriinae. erefore, it is necessary to study the influence of navigation parameters on the energy consumption and voyage of Deep-Argo Otarriinae.
(1) Gliding speed v e gliding speed of the Deep-SOLO float is about 0.06 m/s, the gliding speed of the Deep-SOLO 2 float is about 0.12 m/s, and the gliding speed of the Arvor-C profiling float is about 0.15-0.20 m/s [23,24]. erefore, we set the gliding speed of Deep-Argo 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: e attack angle α is a certain value when Deep-Argo Otarriinae moves steadily, that is to say, equation (34) has a solution and We stipulate that the gliding angle of Deep-Argo Otarriinae is within (−π/2, 0) ∪(0, π/2), but in order to ensure the economicality of navigation, generally the gliding angle is within (−π/4, 0) ∪(0, π/4). erefore, 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 Deep-Argo Otarriinae is within (−7.76°, 0) ∪(0, 7.51°). Under the condition that the velocity of ocean current is almost constant, the gliding angle of Deep-Argo Otarriinae is set uniformly distributed in the above range. 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 Deep-Argo 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 t f is   Mathematical Problems in Engineering interval Δh i to be uniformly distributed between 0m and 50m. Similar to the altimeter opening depth interval, too short running time Δt i of the altimeter may cause Deep-Argo Otarriinae to bottom out while too long running time Δt i will increase the energy consumption. Considering the diving speed of Deep-Argo Otarriinae and the installation position of the altimeter, we set the running time Δt i to be uniformly distributed between 2 and 10 s.

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][29][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 f(X) is a function expression of the energy consumption model, x 1, x 2 , . . . , x M are M independent input parameters affecting energy consumption, respectively. Assuming that the square of f(X) 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.
e left part of equation (30) is called the total variance, which can be described as where D is total variance. Each term on the right side of equation (30) is called partial variance. e single effect of i − th input parameter on the model output can be described as where D i is the single effect of i − th input parameter on the model output. e joint effect of input parameters on the model output can be described as where D ij is the joint effect of input parameters on the model output. e main effect indices S i and the total effect indices S Ti of the parameter x i to the model can be described as where D ∼i is the effect of all parameters but x i . e difference between the main effect and the total effect indices of the parameter x i is a measure of the interaction of x i with all other parameters. e main effect indices S i reflects effect on the model output when the parameter x i changes individually while the total effect indices S Ti reflects the effect on the model output when the parameter x i 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 (n, M).
e total number of model evaluations for the approximation of both the main and total effects indices is therefore equal to where N is the total number of model evaluations, n is the number of samples, and M is the number of independent parameters obtained per sample.

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. erefore, 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. e number of independent parameters M is 10, and the effect indices of each parameter can be guaranteed to converge when the number of samples n is 9600. erefore, the total number of model evaluations N is 115200 according to equation (36). e simulation was conducted on a computer configured with 2.80 GHz i7-7700HQ processor. e 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. e 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.

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 H, and the gliding speed v are at a fairly high level and have a significant impact on the energy consumption of Deep-Argo Otarriinae in a single profile. e higher values of ξ, H, and v 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 Deep-Argo Otarriinae, which have a good match with the actual application. e 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 H is completely set by operators according to the requirements of the working task. ere are no posture adjust module and speed adjust module in our  prototype. e gliding angle ξ and gliding speed v mainly depend on the volume, mass, and shape of Deep-Argo Otarriinae which have been determined during the design process. According to Table 5, we can see clearly that smaller hydrodynamic coefficients K L , K L0 , K D , and K D0 will reduce the energy consumption of Deep-Argo 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.

Conclusions
In this paper, a new modeling method for multiparameter energy consumption model of Deep-Argo Otarriinae was proposed. First, the kinematics model of Deep-Argo Otarriinae was established to show the relationship between various parameters. en, the energy consumption of each subsystem of Deep-Argo 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 Deep-Argo Otarriinae in a single profile accurately and comprehensively. e 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 Deep-Argo 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. e results of sensitivity analysis show that the gliding angle ξ , the diving depth H, and the glide speed v have a greater impact on energy consumption. erefore, 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 Deep-Argo 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 Deep-Argo Otarriinae.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e 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.