Figure 1.
Technical framework of the monitoring and early warning system for ISWs.
| China Geological Survey Chinese Academy of Geological Sciences | Host |
| Citation: |
Dan-yi Su, Bin-bin Guo, Qian-yong Liang, Chu-jin Liang, Fei-long Lin, Su-meng Jiang, Yi-fei Dong, Xue-min Wu, 2023. Application of the monitoring and early warning system for internal solitary waves: Take the second natural gas hydrates production test in the South China Sea as an example, China Geology, 6, 676-684. doi: 10.31035/cg2022043
|
Application of the monitoring and early warning system for internal solitary waves: Take the second natural gas hydrates production test in the South China Sea as an example
-
Abstract
Internal solitary waves (ISWs) contain great energy and have the characteristics of emergency and concealment. To avoid their damage to offshore engineering, a new generation of monitoring and early warning system for ISWs was developed using technologies of double buoys monitoring, intelligent real-time data transmission, and automatic software identification. The system was applied to the second natural gas hydrates (NGHs) production test in the Shenhu Area, South China Sea (SCS) and successfully provided the early warning of ISWs for 173 days (from October 2019 to April 2020). The abrupt changes in the thrust force of the drilling platform under the attack of ISWs were consistent with the early warning information, proving the reliability of this system. A total of 93 ISWs were detected around the drilling platform. Most of them occurred during the spring tides in October–December 2019 and April 2020, while few of them occurred in winter. As suggested by the theoretical model, the full-depth structure of ISWs was a typical current profile of mode-1, and the velocities of wave-induced currents can reach 80 cm/s and 30 cm/s, respectively, in the upper ocean and near the seabed. The ISWs may be primarily generated from the interactions between the topography and semidiurnal tides in the Luzon Strait, and then propagate westward to the drilling platform. This study could serve as an important reference for the early warning of ISWs for offshore engineering construction in the future.
-
1. Introduction
Oceanic internal waves occur in the ocean with stable stratification and their maximum amplitudes are up to tens to hundreds of meters. Low-frequency internal waves have wavelengths of tens or even hundreds of kilometers and propagation speeds of up to several meters per second, while high-frequency internal waves have wavelengths of hundreds to tens of hundreds of meters, propagation speeds of tens of centimeters per second, and periods ranging from 5‒10 mins to 2‒5 hours (Lee CY and Beardsley RC, 1974; Maxworthy T, 1979). The generation of internal waves is closely related to topography, the density stratification of water masses, and the disturbance of tidal and ocean currents (Holloway PE et al., 1997; Gerkema T, 2001; Du T et al., 2008; VanGastel P et al., 2009).
The South China Sea (SCS), which is the largest semi-closed marginal sea in China, boasts abundant oil and natural gas hydrates resources (Xiong LL et al., 2020; He JX et al., 2020). Internal waves are active in the SCS due to the complex geomorphic units, dramatic topographic changes, and pronounced seasonal thermoclines (Hsu MK et al., 2000; Fang XH and Du T, 2005). As a special type of internal waves, internal solitary waves (ISWs) have concentrated energy and large scales. Moreover, ISWs can result in the convergence of seawater and cause abrupt strong currents on the sea surface during their propagation. As a result, sea surface flow fields and density pycnoclines will significantly change. Without early warning of ISWs, abrupt strong wave-induced currents could cause vessels and platform on the sea to be out of control. Furthermore, the vessels may suddenly change course and even collide with offshore platforms. Offshore platforms, which are equipped with a dynamic positioning system nonetheless, respond to strong current attacks with a certain delay. Therefore, abrupt wave-induced currents will make offshore platforms lose balance, thus affecting their ongoing operations and threatening the safety of offshore engineering (You YX et al., 2010; Zhao JR et al., 2018).
The offshore engineering construction and drilling operations in the SCS have been repeatedly affected by ISWs since the 1980s. For instance, the semi-submersible drilling platform deployed near the Dongsha Islands by the China National Offshore Oil Corporation (CNOOC) once swayed by 110° in less than five minutes due to strong ISWs. During the drilling of Well LW13-1-1, the West Aquarius platform experienced an ISW at a speed of 3.8 knots, which moved the platform 37 m, broke the transmission pipelines, and made them wind around the tug propellers. Besides, the Nanhai VIII semi-submersible drilling platform moved hundreds of meters due to multiple attacks by ISWs when it drilled Well LH29-4-1. This resulted in damage to the platform facilities and the collision between the platform and a standby vessel (Ebbesmeyer CC et al., 1991). Therefore, it is of great practical significance and application value to build a monitoring and early warning system to achieve accurate detection and tracking of ISWs and reduce their potential threat to offshore engineering.
Many studies have been conducted on the observation (Ramp SR et al., 2004; Yang YJ et al., 2009; Huang XD et al., 2016) and numerical simulation (Cai SQ et al., 2000, 2012; Zhang Z et al., 2011) of ISWs in the SCS. Owing to the randomness of ISWs, it is difficult to employ existing numerical simulation technologies to achieve the accurate early warning of a single ISW without in-situ observation. Therefore, the monitoring and early warning of ISWs for offshore engineering are mainly implemented using the following three methods: (1) long-distance visual inspection and radar detection from standby vessels; (2) monitoring the temperature, salinity, and ocean currents using standby vessels; and (3) using buoys to monitor ISWs (Goff M et al., 2010; Wang HP et al., 2021; Chun MH et al., 2021). The first method has risks of artificial omissions or errors in identification. Moreover, this method tends to yield large errors in the arrival time of ISWs since it cannot obtain the characteristic parameters such as the propagation speed of ISWs. The second method suffers the disadvantages such as errors induced by the drifting of standby vessels, high shipping fees, and high labor costs. The third method can effectively improve the early warning accuracy of ISWs. However, the performance of this method was often degraded in the presence of data transmission delay and low robustness of buoy systems in the face of severe sea conditions (Hu WJ et al., 2015). Therefore, there is an urgent need to develop a new generation of early warning system using monitoring buoys for ISWs.
This study upgraded the data transmission mode and built a highly robust monitoring and early warning system for ISWs (also referred to as the ISW System), which used the real-time monitoring of double buoys and the technology of automatic software identification. This system was applied in the second natural gas hydrates (NGHs) production test in the Shenhu Area, SCS in 2020 (Ye Jl et al., 2020). In addition, to deepen the understanding of ISWs and provide an important basis for future engineering construction in the SCS, this study analyzed the characteristics and generation mechanisms of ISWs in the study area using observation data.
This paper is organized as follows. Section 2 describes the data and methods used in this study. Section 3 introduces the characteristics of ISWs and the application effects of the ISW System. Section 4 includes the discussion of the generation mechanisms of ISWs in the study area, and the final section provides the conclusions.
2. Data and method
The ISW System was developed using the technologies of double buoys monitoring, intelligent real-time data transmission, and automatic software identification. In this system, two buoys deployed at different locations along the arrival direction of ISWs were used to monitor and send the oceanic data to an onshore server through a communication satellite. The server was adopted to receive, process, and store data from buoys and then sent data through the Internet to the ISW monitoring software installed on an offshore platform and in an onshore monitoring laboratory. The ISW monitoring software would automatically identify and output signals of ISWs and issue warnings, thus providing the early warning of ISWs for the offshore platform. Moreover, the onshore monitoring laboratory was equipped with a 24-hour artificial service, where the personnel would immediately inform operators on the offshore platform and proposed suggestions once they discovered the signals of ISWs (Fig. 1).
2.1 Real-time monitoring using buoys
As shown by the historical observations near the drilling platform for the second NGHs production test, the ISWs in this area propagated at approximately 2.1‒2.5 m/s (average: 2.2 m/s) in the directions of approximately 270° (western) (Jackson CR et al., 2009). Therefore, buoy No. 1 of the ISW System was deployed 27 km in the direction of 95° to the drilling platform, and buoy No. 2 was placed 13.5 km in the direction of 85° to the drilling platform (Fig. 2). They were used to provide the early warning of ISWs for the drilling platform 1.5–3 hours ahead of time.
Figure 2. Locations of the drilling platform and monitoring buoys in the second NGHs production test. The red box denotes the study area (i.e., the target area for the monitoring and early warning of ISWs); the yellow box denotes the site of tide forecast data; the red triangles indicate the buoy stations; the black square is the position of the drilling platform; and the black lines indicate the water depth (unit: m).The BETA-3M-C buoys were used in the ISW System. Each of them consisted of a buoy body, a sensor subsystem, an anchoring subsystem, and subsystems for data acquisition, data transmission, power supply, and tracking. Moreover, they were equipped with positioning and warning devices such as GPS and anchor lights (Fig. 3). Their buoy bodies had maximum swing angles of less than 25°, less than 35°, and less than 45°, respectively in the case of the wave height of 5 m, 10 m, and 15 m. The buoy can automatically restore itself to the normal attitude even when it was inclined at an angle of 90° and can be placed in an area with a wind speed of less than 100 m/s, or with a wave height of less than 20 m, or with a flow velocity of less than 5 m/s. The sensors equipped in the BETA-3M-C buoys mainly included conductivity-temperature-depth sensors (CTDs, SBE 37-SM) and Acoustic Doppler Current Profilers sensors (ADCPs; TDRI Workhorse Quartermaster 150 kHz and TDRI Workhorse Sentinel 300 kHz). These sensors can monitor the pressure, temperature, and salinity of surface seawater and the velocity profiles of upper ocean currents (Table 1). The monitoring data and the operating state of buoys were transmitted to the onshore server in real time through an Iridium 9522B satellite transceiver.
Figure 3. Structure of a buoy used in the ISW System. Physical picture (a); overall structure (b).Table 1. Monitoring indices and data collection settings of buoy sensors.Sensor Element Range Precision Sampling frequency Data transmission interval Depth CTD Temperature −5°C–45°C ±0.002°C 10 mins 1 hour Sea surface Conductivity 0–70 mS/cm ±0.003 mS/cm ADCP Velocity ±5 m/s ±0.5%×
measured value or
±5 mm/s ±0.5% × measured value or
±5 mm/s3 mins 12 mins 0–124 m Direction 0–360° ±2° | Show Table
DownLoad:
CSV
The ISW System was designed to comprise double buoys with dual ADCPs each and an intelligent data transmission mode, which further improved its robustness. Specifically, each of the buoys was equipped with a 150 kHz and a 300 kHz ADCP, thus successfully combining a high measuring range and a high accuracy. Moreover, the ISW System can provide effective early warning of ISWs as long as one of the ADCPs normally worked. In addition, given the massive amounts of monitoring data, the ISW System adopted a new intelligent data transmission mode to avoid data transmission failure caused by factors such as bad weather and sea conditions. In detail, the latest monitoring data would be preferentially transmitted when large amounts of historical monitoring data were waiting to be sent, thus ensuring the timeliness of early warning. Furthermore, the users can remotely control the data transmission switches, the priority orders of the four ADCPs and other auxiliary sensors in real time according to the states of data transmission and the Iridium 9522B satellite transceiver. The ISW System operation and real-time data transmission performed well in the second NGHs production test. The data from CTDs and ADCPs were transmitted at an interval of 1 hour and 12 mins, respectively. Owing to such high-frequency data transmission, the delay in the early warning of ISWs was less than 12 mins.
2.2 Data processing and software identification of ISWs
The ISW System adopted automatic software identification technology to detect ISWs. After receiving the observation data transmitted from buoys in real time, the ISW monitoring software automatically conducted data processing, identified the signals of ISWs, and issued warnings. Firstly, based on the finite depth model (Joseph RI et al., 1978) and the temperature vs. salinity profiles, the authors calculated the parameters of various modes of ISWs in the study area (including the velocity of wave-induced current and the depth, acceleration, and duration of ISWs) and built up a database. During the application, the ISW monitoring software processed the observation data transmitted from the buoys in real time and obtained the characteristics of upper ocean currents. Afterward, it compared these characteristics with the database and calculated the similarity. Once the similarity exceeded a certain threshold, the software automatically sent out warning signals of ISWs. Moreover, the propagation speed of an ISW can be calculated based on the distance between the two buoys and the time interval of the ISW passing through them. Subsequently, the time when the ISW arrived at the drilling platform was predicted according to the distance from buoy No. 2 to the drilling platform, and previous information was correspondingly updated. In sum, the ISW information issued by the software included the time when an ISW passed through two buoys; the velocity and depth of the maximum wave-induced current observed; the maximum velocities of the upwelling and downwelling currents; and the predicted time when an ISW arrived at the drilling platform.
To reduce potential risks such as software failures and omissions caused by the automatic identification algorithm, a 24-hour artificial service was arranged for the onshore monitoring laboratory, where manual auxiliary diagnoses were made on the following bases: (1) ISWs may occur if the horizontal velocity of upper ocean currents surges to more than 60 cm/s in 6–20 mins; (2) ISWs may occur if strong upwelling (>5 cm/s) and downwelling (<−5 cm/s) simultaneously take place in the upper ocean in 10–30 mins; (3) westward or northwestward ISWs may occur if the upper ocean currents are greater than 45 cm/s in the west or northwest direction. Besides, the personnel further confirmed signals of ISWs by comparing the propagation velocities measured by a single buoy with those calculated using the observation data from the two buoys. They would immediately report the information of ISWs whose modes were not registered in the database to the operators on the drilling platform, and confirmed the impacts of these ISWs. If any impact was caused, it can be determined that an ISW with a special mode had occurred and would be saved into the database, thus constantly updating the database. In addition, the software can monitor the status of buoys in real time and send out information on their abnormal working states, displacement, and failures, which assist in maintaining the buoys in time to keep them operating normally and stably.
To fulfill the safety construction requirements of the second NGHs production test, the study defined the intensity grading of ISWs based on the velocities of wave-induced currents and proposed corresponding measures (Table 2).
Table 2. Measures against ISWs of different intensity gradings.ISW intensity Velocity/(m/s) Measures Low u<0.3 Recording the information of ISWs; the drilling platform need not take special measures. Medium 0.3 
u<0.45The ISW System issues warnings; the drilling platform need not take special measures. High 0.45 
u<0.6The ISW System sends out emergency responses; the drilling platform adjusts its heading direction to the arrival direction of ISWs to reduce their lateral impacts on the platform. Very high u 
0.6The ISW System sends out emergency responses; the drilling platform stops all submarine operations. | Show TableDownLoad:
CSV
2.3 Simulation of ISWs
The ADCPs equipped in the buoys cannot obtain the full-depth vertical structure of ISWs since their maximum measuring range was only 124 m. However, the impacts of ISWs on the seabed also need to be considered in offshore engineering. Given this, the authors simulated the vertical structures of ISWs using the finite depth theory (Grimshaw R, 1985; Pierini S, 1989). The amplitude of an ISW in a finite depth fluid satisfies the following equation:
∂η∂t+c∂η∂x+αη∂η∂x−βFc2H0∂2∂x2∫+∞−∞η(x′,t)[cothπ(x−x′)2H0−sgn(x−x′)]dx′=0, (1) where η and c are the amplitude and the linear phase speed of ISW, respectively; α is a nonlinear parameter; H0 is the water depth; βF is the dispersion parameter in the finite depth theory. Solving this equation yields the following horizontal velocity component (Cai SQ et al., 2015):
u=−cη0cosh2(ϕ)+sinh2(ϕ)/a2b2∙dWdz, (2) where
ϕ=a(x−Vt) , and a is a parameter similar to the wavenumber and satisfies equationsabtan(aH0)=1 andb=−4βFcαη0 . Moreover, the nonlinear phase speed V can be calculated based on the distance between the two buoys and the time interval of an ISW passing through them. Among them, the amplitude and nonlinear phase speed are two key parameters to be determined, from one of which other parameters of ISWs involved in the finite depth theory can be derived. Therefore, the authors repeatedly assumed the amplitude and calculated the horizontal velocity in the upper ocean. Until the calculated horizontal velocity was the closest to the measured velocity, the finally assumed amplitude was considered as the actual amplitude. Then, the full-depth flow field was calculated using this amplitude.3. Application results
During the second NGHs production test, two buoys were deployed to the east of the drilling platform to ensure the safety of the engineering. They successfully provided real-time monitoring and early warning of ISWs for 173 days (from October 12, 2019 to April 17, 2020).
3.1 Monitoring effects and statistical characteristics of ISWs
A total of 93 ISWs were detected in the study area during the second NGHs production test, including 8 with very high intensity, 11 with high intensity, 25 with medium intensity, and 49 with low intensity. These ISWs mainly occurred at 8–96 m in the upper ocean. Moreover, the velocity, directions, and durations of the wave-induced currents were 19.1–79.8 cm/s, 250°–303°, and 3–27 mins, respectively (Table 3). As for the 44 ISWs with medium or higher intensity, the ISW System issued effective early warnings to the drilling platform and provided accurate information of ISWs in time, thus ensuring the safety of the second NGHs production test.
Table 3. Characteristics of ISWs (from October 12, 2019 to April 17, 2020).Date ISW number Calculation
methodVelocity/
(cm/s)Direction/° Depth/
mDuration/
min2019.10 16 range 20.9‒66.9 250‒298 8‒80 6‒24 average 46.7 268 53 14 2019.11 17 range 20.2‒68.0 256‒286 8‒92 6‒21 average 44.9 269 58 11 2019.12 5 range 28.3‒68.7 275‒295 12‒96 6‒18 average 47.5 285 47 13 2020.01 9 range 20.3‒55.0 282‒290 42‒57 6‒9 average 43.5 286 50 8 2020.02 0 range / / / / average / / / / 2020.03 8 range 20.4‒37.6 273‒277 28‒20 17‒18 average 37.1 275 24 18 2020.04 38 range 19.1‒79.8 255‒303 32‒92 3‒27 average 47.0 276 65 11 | Show TableDownLoad:
CSV
Fig. 4 shows the strongest ISW detected during the second NGHs production test. This ISW, which propagated from the east toward the drilling platform, arrived at buoy No. 2 at 12:42 on April 11, 2020. The velocities of currents at 60 m surged from 40 cm/s to 108 cm/s (Fig. 4a) in the direction of approximately 262°. The currents with velocity of greater than 60 cm/s lasted for 30 mins (Fig. 4c). Then this ISW passed by buoy No. 1 at 13:57 that day. Its peaks detected by both buoys had the time interval of 75 mins. Since the distance between the two buoys was 13.5 km, the propagation speed of the strongest ISW was 10.8 km/h. The ISW System sent out an early warning to the drilling platform and predicted that the strongest ISW would arrive at 15:12 according to its propagation speed and the distance between buoy No. 1 and the drilling platform (13.5 km). When an ISW arrived at the drilling platform, the dynamic positioning system of the platform would exert a reaction force (i.e., the thrust force) to prevent the platform from moving. The thrust force of the drilling platform showed that the ISW attacked the platform at 15:12 and that the enhanced thrust force lasted for approximately 30 mins (Fig. 5), which was consistent with the early warning information. It proved that the new generation of the monitoring and early warning system can provide accurate early warning of ISWs for offshore platforms.
Figure 4. Velocity profiles of upper ocean currents (a, b) and average velocity at 0–100 m (c, d) when the strongest ISW passed through buoys No. 2 (left) and No. 1 (right) on April 11, 2020.
Figure 5. Changes in the thrust force of the drilling platform when the strongest ISW passed on April 11, 2020 (unit: tons).3.2 Temporal variations of ISWs
Based on the 6-month observation data, the authors analyzed and summarized the temporal variations of ISWs in the study area. Fig. 6 shows the statistical characteristics of ISWs and the time series of semidiurnal tides near the Luzon Strait (Fig. 2; yellow box). The tide forecast data in this figure originated from the global tide model TPXO (Egbert GD et al., 1994; Egbert GD and Erofeeva SY, 2002). It can be inferred that the ISWs in the study area mainly occur during spring tides from October to December, 2019 and in April, 2020. For instance, ISWs were generally detected every day from October 17 to October 22 and from October 30 to November 4. These two periods had an interval of approximately 15 days, which exactly corresponded to the spring tide period.
Figure 6. Time series of ISWs in the study area. Velocity of wave-induced currents (a); depth of the maximum velocity (b); duration of the ISWs (c); semidiurnal tides near the Luzon Strait (d).As shown by the monthly statistics, the ISW number, the velocities and durations of wave-induced currents gradually decreased from October 2019 to March 2020 (Table 3, Figs. 6a and 6c), and the average direction progressively shifted from southwest to west. However, the ISWs significantly intensified and occurred more frequently in April 2020. That is to say, the ISWs in the study area had apparent seasonality and occurred the least frequently in winter. This phenomenon may be related to the seasonal changes in the upper density pycnoclines in the SCS and the westward energy flux from the Luzon Strait (Li MJ, 2012; Chen L et al., 2019). In addition, as shown by the daily statistics, the ISWs, especially those with very high intensity, were concentrated at 05:00–06:00 and 13:00–14:00 (Fig. 7). This occurred possibly because the energy of the ISWs was primarily sourced from semidiurnal tides (Duda TF et al., 2004; Li Q and Farmer DM, 2011).
Figure 7. Daily distribution of the velocity of wave-induced currents during the monitoring.3.3 Vertical structure of ISWs
Fig. 8 shows the profile and the vertical shear profile of the horizontal velocity of wave-induced currents when the strongest ISW passed by buoy No. 2 on April 11, 2020, as well as the temperature vs. salinity profile near the buoy in April 2021. It can be found that the strongest shear occurred in the thermoclines (approximately 80–90 m, Figs. 8b and 8c) when the ISW approached and left. Besides, the authors simulated the vertical structure of the ISWs using the observation data from buoy No. 2 and the finite depth theory (Fig. 9). The simulation results at 0‒100 m were consistent with the observation results. Specifically, both results showed that the maximum velocities of wave-induced currents were up to 100 cm/s, the cores of the currents were roughly located at 100 m, and the wave-induced currents with velocities of greater than 60 cm/s lasted for approximately 30 mins (Figs. 9a and 9b). The simulated velocity profile indicated that the velocity of wave-induced currents at 1500 m was up to 30 cm/s and its direction was opposite to that in the upper ocean (Fig. 9c). These results suggested that the vertical structure of ISWs in the study area was a typical current profile of mode-1.
Figure 8. The horizontal velocity profile (a) and the vertical shear profile (b) when the strongest ISW passed by buoy No. 2 on April 11, 2020, and the temperature vs. salinity profile (c) near buoy No. 2 in April 2021.
Figure 9. Vertical structure of the strongest ISW passed by buoy No. 2 on April 11, 2020. Simulation results at 0‒100 m (a); monitoring results at 0‒100 m (b); simulated full-depth structure (c).4. Discussion
According to previous studies, the ISWs in the SCS are mostly generated from the Luzon Strait due to lee wave and the interactions between tides and local topography (Liu AK et al., 1998, 2004; Zhao Z et al., 2004; Buijsman MC et al., 2010). Based on the empirical formula proposed by Jackson CR et al. (2009), the authors calculated the time it takes for the ISWs originating in the Luzon Strait to propagate to different locations in the SCS.
(∂T(x,y)∂x)2+(∂T(x,y)∂y)2=1C(x,y)2 (3) C(x,y)=Cmax√tanh(B1+H(x,y)B2) (4) In equation (3), T(x,y) is the time it takes for waves to propagate from the start position (x0,y0) to other position (x,y), and C(x,y) is the empirical formula of the phase speed of waves (Jackson CR et al., 2009), which is applicable to the SCS. In equation (4), H(x,y) is the water depth; Cmax is the maximum phase speed (~2.971 m/s); B1 and B2 are terms that determine how rapidly variations in depth will affect the slowing and refracting of the waves. In this study, B1 = 0.003 and B2 = 1390.758.
As suggested by the temporal distribution of wave-induced currents and tidal currents in the Luzon Strait (Fig. 6d), the time when the ISWs intensively occurred in the study area lagged the period of semidiurnal tides in the Luzon Strait for 2–3 days. This delay exactly matches the time it takes for the ISWs originating in the Luzon Strait to propagate to the study area (approximately 2.7 days, Fig. 10). Therefore, it can be inferred that some of the ISWs in the study area mainly originated from the interactions between the topography and the semidiurnal tides in the Luzon Strait.
Figure 10. Time it takes for ISWs to propagate from the Luzon Strait to the study area (black solid lines; unit: day). The colorful solid lines denote water depth (unit: m), and the red box denotes the study area.5. Conclusions
This study developed a new generation of monitoring and early warning system for ISWs in order to provide accurate early warning of ISWs for offshore engineering and better ensure its safe implementation. The ISW System was successfully applied to the second NGHs production test in the Shenhu Area, SCS from October 12, 2019 to April 17, 2020.
Based on two buoys deployed near the drilling platform, the ISW System obtained the observation data in real time. Afterward, the data were transmitted to an onshore server through a satellite, and sent through the Internet to the ISW monitoring software. Then, the independently developed ISW monitoring software identified the signals of ISWs and issued warnings in time. The ISW System adopted double buoys with dual ADCPs each, an intelligent data transmission mode, and automatic software identification technology, which greatly improved its robustness in the face of severe sea conditions.
During the second NGHs production test, 93 ISWs were detected in total, including 8 with very high intensity and 11 with high intensity. As shown by the thrust force of the drilling platform, the time at which the strongest ISW reached the platform and the duration of the wave-induced current were consistent with the early warning information. This consistency suggested that the ISW System can provide accurate and effective warning signals of ISWs for offshore platforms. The statistical results showed that the ISWs in the study area occurred the least frequently in winter and that most of them occurred during the spring tides from October to December, 2019 and in April, 2020. The ISWs mainly appeared at 8–96 m, with durations of 3–27 mins. The velocities of wave-induced currents were at 19.1–79.8 cm/s in directions of 250°–303°.
As revealed by the theoretical model and the observation data, the vertical structure of ISWs in the study area was a typical current profile of mode-1, and the wave-induced currents near the seabed can reach 30 cm/s. The ISWs were primarily generated from the interactions between the semidiurnal tides and local topography in the Luzon Strait. The characteristics of ISWs revealed by this study can provide important information for the design of other offshore engineering in the Shenhu Area, South China Sea. Moreover, the ISW System is expected to be widely used in offshore engineering construction.
CRediT authorship contribution statement
Bin-bin Guo, Qian-yong Liang and Chu-jin Liang conceived of the presented idea. Yi-fei Dong and Xue-min Wu carried out the implementation, Dan-yi Su, Chu-jin Liang, Fei-long Lin and Su-meng Jiang developed the theoretical formalism, performed the analytic calculations, and contributed to the interpretation of the results. Dan-yi Su and Bin-bin Guo wrote the manuscript with input from all authors.
Declaration of competing interest
The authors declare no conflicts of interest.
Acknowledgment
This study was funded by the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0307), the Marine Geological Survey Program of China Geological Survey (DD20190218, DD20221706), the Key Program of Marine Economy Development Special Foundation of Department of Natural Resources of Guangdong Province (GDNRC [2020] 043), and the National Natural Science Foundation of China (41806074, 41730528). The authors would like to extend their gratitude to the staff of all expeditions of R/V Haiyangdizhi 6 and R/V Haiyangdizhi 4 for their assistance in the setting, picking up, and maintenance of the buoys used in this study.
-
References
Buijsman MC, Kanarska Y, McWilliams JC. 2010. On the generation and evolution of nonlinear internal waves in the South China Sea. Journal of Geophysical Research: Oceans, 115(C2). doi: 10.1029/2009jc005275. Cai SX, Huang QZ, Qiu Z, Mai BQ. 2000. A numerical study on dynamic mechanism of internal tide. Tropic Oceanology, 19(3), 27–32 (in Chinese). Cai SQ, Xie JS, He JL. 2012. An Overview of internal solitary waves in the South China Sea. Surveys in Geophysics, 33(5), 927–943. doi: 10.1007/s10712-012-9176-0. Cai SQ, et al. 2015. Numerical model of internal isolated wave and its application in the South China Sea. Beijing. China Ocean Press. 1–158 (in Chinese). Chen L, Zheng QA, Xiong XJ, Yuan YL, Xie HR, Guo YL, Yu L, Yun SJ. 2019. Dynamic and statistical features of internal solitary waves on the continental slope in the northern South China Sea derived from mooring observations. Journal of Geophysical Research: Oceans, 124, 4078–4097. doi: 10.1029/2018jc014843. Chun MH, Leng XD, Liang SS, Wang Z, Luo XQ, Yang XD, Xu S. 2021. Research and application of internal wave monitoring and warning technology in South China Sea. Petroleum Engineering Construction, 47(4), 70–74 (in Chinese). Du T, Tseng YH, Yan XH. 2008. Impacts of tidal currents and Kuroshio intrusion on the generation of nonlinear internal waves in Luzon Strait. Journal of Geophysical Research:Atmospheres, 113, C08015. doi: 10.1029/2007jc004294. Duda TF, Lynch JF, Irish JD, Beardsley RC, Ramp SR, Chiu CS, Tang TY, Yang YJ. 2004. Internal tide and nonlinear internal wave behavior at the continental slope in the Northern South China Sea. IEEE Journal of Oceanic Engineering, 29(4), 1105–1130. doi: 10.1109/joe.2004.836998. Ebbesmeyer CC, Coomes CA, Hamiton RC. 1991. New observations on internal waves (solitons) in the South China Sea using an acoustic doppler current profiler. Marine Technology Society Journal, 91, 165–175. Egbert GD, Bennett AF, Foreman MGG. 1994. TOPEX/POSEIDON tides estimated using a global inverse model. Journal of Geophysical Research: Atmospheres, 99(C12), 24821–24852. doi: 10.1029/94jc01894. Egbert GD, Erofeeva SY. 2002. Efficient inverse modeling of barotropic ocean tides. Journal of Atmospheric and Oceanic Technology, 19(2), 183–204. doi: 10.1175/1520-0426(2002)0192.0.CO;2. Fang XH, Du T. 2005. Fundamentals of Oceanic Internal Waves and Internal Waves in the China Seas. China Qingdao, Ocean University of China Press, 337 (in Chinese). Gerkema T. 2001. Internal and interfacial tides: beam scattering and local generation of solitary waves. Journal of Marine Research, 59, 227–255. doi: 10.1357/002224001762882646. Goff M, Jeans G, and Harrington-Missin L, Baschenis C. 2010. Soliton early warning system for offshore applications. Oceanology International (OI), London, United Kingdom. Grimshaw R, 1985. Evolution equations for weakly nonlinear long internal waves in a rotating fluid. Stud. Studies in Applied Mathematics, 73(1), 1–33. doi: 10.1002/sapm19857311. He JX, Zhong CM, Yao YJ, Yan P, Wang YL, Wan ZF, Guan J, Zhang JF. The exploration and production test of gas hydrate and its research progress and exploration prospect in the northern South China Sea. Marine Geology Frontiers, 2020, 36(12): 1-14. doi: 10.16028/j.1009-2722.2020.127 (in Chinese). Holloway PE, Pelinovsky E, Talipova T, Barnes B. 1997. A nonlinear model of internal tide transformation on the Australian North West Shelf. Journal of Physical Oceanography, 27(6), 871–896. doi: 10.1175/1520-0485(1997)027<0871:ANMOIT>2.0.CO;2. Hsu MK, Liu AK, Liu C. 2000. A study of internal waves in the China Seas and Yellow Sea using SAR. Continental Shelf Research, 20(4-5), 389–410. doi: 10.1016/s0278-4343(99)00078-3. Huang XD, Chen ZH, Zhao W, Zhang ZW, Zhou C, Yang QX, Tian JW. 2016. An extreme internal solitary wave event observed in the northern South China Sea. Scientific Reports, 6(1), 30041. doi: 10.1038/srep30041. Hu WJ, Liu ZL, Chen B. 2015. Impacts of internal waves in the South China Sea on deepwater drilling safety and corresponding countermeasures. Oil Drilling & Production Technology, 37(1), 160–162 (in Chinese). Jackson CR. 2009. An empirical model for estimating the geographic location of nonlinear internal solitary waves. Journal of Atmospheric and Oceanic Technology. 26(10), 2243–2255. doi: 10.1175/2009jtecho638.1. Joseph RI, Egri R. 1978. Multi-soliton solutions in a finite depth fluid. Journal of Physics A:Mathematical and General, 11(5), 97–102. doi: 10.1088/0305-4470/11/5/002. Lee CY, Beardsley RC. 1974. The generation of long nonlinear internal waves in a weakly stratified shear flow. Journal of Geophysical Research: Atmospheres, 79(3), 453–462. doi: 10.1029/jc079i003p00453. Li MJ. 2012. Three-dimensional Numerical Simulation and Analysis of Internal Tides in the Northern South China Sea. Beijing, University of Chinese Academy of Sciences, Ph. D thesis, 1–116 (in Chinese). Li Q, Farmer DM. 2011. The generation and evolution of nonlinear internal waves in the deep basin of the South China Sea. Journal of Physical Oceanography, 41(7), 1345–1363. doi: 10.1175/2011jpo4587.1. Liu AK, Hsu MK, Liang NK. 1998. Evolution of nonlinear internal waves in China Seas, Journal of Geophysical Research: Atmospheres, 103 (C4), 7995–8008. Liu AK, Ramp SR, Zhao YH, Tang TY. 2004. A case study of internal solitary wave propagation during ASIAEX 2001. IEEE Journal of Oceanic Engineering, 29(4), 1144–1156. doi: 10.1109/joe.2004.841392. Maxworthy T. 1979. A note on the internal solitary waves produced by tidal flow over a three-dimensional ridge. Journal of Geophysical Research:Atmospheres, 84(C1), 338. doi: 10.1029/jc084ic01p00338. Pierini S. 1989. A model for the Alboran Sea internal solitary wave. Journal of Physical Oceanography, 19, 755–772. doi: 10.1175/1520-0485(1989)019<0755:AMFTAS>2.0.CO;2. Ramp SR, Tang TY, Duda TF, Lynch JF, Liu AK, Chiu CS, Bahr FL, Kim HR, Yang YJ. 2004. Internal solitons in the northeastern South China Sea Part I: sources and deep water propagation. Journal of Oceanic Engineering, 29(4), 1157–1181. doi: 10.1109/JOE.2004.840839. VanGastel P, Ivey GN, Meuleners MJ, Antenucci JP, Fringer O. 2009. The variability of the large-amplitude internal wave field on the Australian North West Shelf. Continental Shelf Research, 29(11–12), 1373–1383. doi: 10.1016/j.csr.2009.02.006. Wang HP, Chen L, Guo YL, Cao LZ, Hu XM, Xiong XJ. 2021. Observing, identification and early warning technology of internal solitary wave and its application in Liuhua 16-2 oilfield group development project. The Ocean Engineering, 39(2), 162–170 (in Chinese). Xiong LL, Yang CP, Wu JQ, Gao HF, Yao YJ, Li XJ, Zhu XY, Cheng ZH. 2020. Characteristics of stratigraphic sedimentary in the south-north continental margin basin of the South China Sea and its differential control on hydrocarbon accumulation. Geology in China, 47(5), 1407–1425. doi: 10.12029/gc20200508. Yang YJ, Fang YC, Chang MH, Ramp SR, Kao CC, Tang TY. 2009. Observations of second baroclinic mode internal solitary waves on the continental slope of the northern South China Sea. Journal of Geophysical Research:Oceans, 114, C10003. doi: 10.1029/2009jc005318. Ye JL, Qin XW, Xie WW, Lu HL, Ma BJ, Qiu HJ, Liang JQ, Lu LA, Kuang ZG, Lu C, Liang QY, Wei SP, Yu YJ, Liu CS, Li B, Shen KX, Shi HX, Lu QP, Li J, Kou BB, Song G, Li B, Zhang HE, Lu HF, Ma C, Dong YF, Bian H. 2020. The second natural gas hydrate production test in the South China Sea. China Geology, 3(2), 197–209. doi: 10.31035/cg2020043. You YX, Hu TQ, Xu H, Yu Z, Fang YL, Shi ZM, Qu Y, Xiao Y. 2010. Experiments on interaction of internal waves with semi-submersible platform in a stratified fluid. Chinese Journal of Theoretical and Applied Mechanics, 42(3), 400–406 (in Chinese). Zhang Z, Fringer OB, Ramp SR. 2011. Three-dimensional, nonhydrostatic numerical simulation of nonlinear internal wave generation and propagation in the South China Sea. Journal of Geophysical Research:Oceans, 116, C05022. doi: 10.1029/2010jc006424. Zhao Z, Klemas V, Zheng Q, Yan XH. 2004. Remote sensing evidence for baroclinic tide origin of internal solitary waves in the northeastern South China Sea. Geophysical Research Letters, 31(6), L06302. doi: 10.1029/2003gl019077. Zhao JR, Xie B, Su J, Wang SS, Li B. 2018. Transient mooring analysis of a semi-submersible drilling platform under internal wave condition. Ship Science and Technology, 40(21), 75–79 (in Chinese). -
Access History
Figures(10)
Tables(3)
Export File
Citation
Dan-yi Su, Bin-bin Guo, Qian-yong Liang, Chu-jin Liang, Fei-long Lin, Su-meng Jiang, Yi-fei Dong, Xue-min Wu, 2023. Application of the monitoring and early warning system for internal solitary waves: Take the second natural gas hydrates production test in the South China Sea as an example, China Geology, 6, 676-684. doi: 10.31035/cg2022043
Format
Content
- Table 1. Monitoring indices and data collection settings of buoy sensors.
Sensor Element Range Precision Sampling frequency Data transmission interval Depth CTD Temperature −5°C–45°C ±0.002°C 10 mins 1 hour Sea surface Conductivity 0–70 mS/cm ±0.003 mS/cm ADCP Velocity ±5 m/s ±0.5%×
measured value or
±5 mm/s ±0.5% × measured value or
±5 mm/s3 mins 12 mins 0–124 m Direction 0–360° ±2° | Show TableDownLoad:
CSV
- Table 2. Measures against ISWs of different intensity gradings.
ISW intensity Velocity/(m/s) Measures Low u<0.3 Recording the information of ISWs; the drilling platform need not take special measures. Medium 0.3 u<0.45The ISW System issues warnings; the drilling platform need not take special measures. High 0.45 u<0.6The ISW System sends out emergency responses; the drilling platform adjusts its heading direction to the arrival direction of ISWs to reduce their lateral impacts on the platform. Very high u 0.6The ISW System sends out emergency responses; the drilling platform stops all submarine operations. | Show TableDownLoad:
CSV
- Table 3. Characteristics of ISWs (from October 12, 2019 to April 17, 2020).
Date ISW number Calculation
methodVelocity/
(cm/s)Direction/° Depth/
mDuration/
min2019.10 16 range 20.9‒66.9 250‒298 8‒80 6‒24 average 46.7 268 53 14 2019.11 17 range 20.2‒68.0 256‒286 8‒92 6‒21 average 44.9 269 58 11 2019.12 5 range 28.3‒68.7 275‒295 12‒96 6‒18 average 47.5 285 47 13 2020.01 9 range 20.3‒55.0 282‒290 42‒57 6‒9 average 43.5 286 50 8 2020.02 0 range / / / / average / / / / 2020.03 8 range 20.4‒37.6 273‒277 28‒20 17‒18 average 37.1 275 24 18 2020.04 38 range 19.1‒79.8 255‒303 32‒92 3‒27 average 47.0 276 65 11 | Show TableDownLoad:
CSV
Sensor Element Range Precision Sampling frequency Data transmission interval Depth CTD Temperature −5°C–45°C ±0.002°C 10 mins 1 hour Sea surface Conductivity 0–70 mS/cm ±0.003 mS/cm ADCP Velocity ±5 m/s ±0.5%×
measured value or
±5 mm/s ±0.5% × measured value or
±5 mm/s3 mins 12 mins 0–124 m Direction 0–360° ±2° ISW intensity Velocity/(m/s) Measures Low u<0.3 Recording the information of ISWs; the drilling platform need not take special measures. Medium 0.3 u<0.45The ISW System issues warnings; the drilling platform need not take special measures. High 0.45 u<0.6The ISW System sends out emergency responses; the drilling platform adjusts its heading direction to the arrival direction of ISWs to reduce their lateral impacts on the platform. Very high u 0.6The ISW System sends out emergency responses; the drilling platform stops all submarine operations. Date ISW number Calculation
methodVelocity/
(cm/s)Direction/° Depth/
mDuration/
min2019.10 16 range 20.9‒66.9 250‒298 8‒80 6‒24 average 46.7 268 53 14 2019.11 17 range 20.2‒68.0 256‒286 8‒92 6‒21 average 44.9 269 58 11 2019.12 5 range 28.3‒68.7 275‒295 12‒96 6‒18 average 47.5 285 47 13 2020.01 9 range 20.3‒55.0 282‒290 42‒57 6‒9 average 43.5 286 50 8 2020.02 0 range / / / / average / / / / 2020.03 8 range 20.4‒37.6 273‒277 28‒20 17‒18 average 37.1 275 24 18 2020.04 38 range 19.1‒79.8 255‒303 32‒92 3‒27 average 47.0 276 65 11
-
Figure 1.
Technical framework of the monitoring and early warning system for ISWs.
-
Figure 2.
Locations of the drilling platform and monitoring buoys in the second NGHs production test. The red box denotes the study area (i.e., the target area for the monitoring and early warning of ISWs); the yellow box denotes the site of tide forecast data; the red triangles indicate the buoy stations; the black square is the position of the drilling platform; and the black lines indicate the water depth (unit: m).
-
Figure 3.
Structure of a buoy used in the ISW System. Physical picture (a); overall structure (b).
-
Figure 4.
Velocity profiles of upper ocean currents (a, b) and average velocity at 0–100 m (c, d) when the strongest ISW passed through buoys No. 2 (left) and No. 1 (right) on April 11, 2020.
-
Figure 5.
Changes in the thrust force of the drilling platform when the strongest ISW passed on April 11, 2020 (unit: tons).
-
Figure 6.
Time series of ISWs in the study area. Velocity of wave-induced currents (a); depth of the maximum velocity (b); duration of the ISWs (c); semidiurnal tides near the Luzon Strait (d).
-
Figure 7.
Daily distribution of the velocity of wave-induced currents during the monitoring.
-
Figure 8.
The horizontal velocity profile (a) and the vertical shear profile (b) when the strongest ISW passed by buoy No. 2 on April 11, 2020, and the temperature vs. salinity profile (c) near buoy No. 2 in April 2021.
-
Figure 9.
Vertical structure of the strongest ISW passed by buoy No. 2 on April 11, 2020. Simulation results at 0‒100 m (a); monitoring results at 0‒100 m (b); simulated full-depth structure (c).
-
Figure 10.
Time it takes for ISWs to propagate from the Luzon Strait to the study area (black solid lines; unit: day). The colorful solid lines denote water depth (unit: m), and the red box denotes the study area.