Dynamic Pressure Waves Research Articles

Research and Application of Leakage Monitoring Techniques for Gas Field Water Pipelines

Abstract To study the effectiveness and application effects of leakage monitoring techniques for gas field water pipelines in the Sichuan-Chongqing region, seven kinds of flow limitation holes with different diameters were used to simulate the leakage of the gas field water pipeline to carry out a field test. The monitoring conditions of two leakage monitoring techniques, namely dynamic pressure wave compound method and negative pressure wave and transmission balance method, were compared and analyzed under different working conditions. The test results showed that: (1) When the pipeline was full, the leakage with apertures being smaller than 2 mm could not be effectively monitored. A single-end alarm could be achieved using the dynamic pressure wave compound method, but the positioning was ineffective. An alarm could not be effectively achieved using the transmission balance method. (2) When the pipeline was full, for the leakage with apertures being 2 mm and above, early warning could be well achieved using the both methods, and the monitoring effects were nearly the same, with the average positioning error being about 48 m and the average alarm response time being about 58 s. (3) When there were empty pipe slug flows in the pipeline, only apertures being 5 mm and above could be alarmed; however, positioning could not be effectively achieved, and system technical indicators could not be reached. (4) When there were large elevation differences in the gas field water pipeline and empty pipe slug flows could not be eliminated, the both techniques were not applicable.

Three-dimensional physical experimental study on mechanisms and influencing factors of CO2 huff-n-puff and flooding process in shale reservoirs after fracturing

Carbon dioxide (CO2) capture and geological storage technology is a promising strategy for enhancing oil and gas production and mitigating climate change in unconventional formations. In shale reservoirs, production decline is rapid under volumetric fracturing development. Currently, CO2 injection and huff-n-puff techniques are cost-effective methods to boost shale oil recovery. However, traditional huff-n-puff techniques are limited in capturing reservoir heterogeneity and intricate fracture conditions. This study employed a three-dimensional physical experimental apparatus to explore CO₂ huff-n-puff dynamics across six distinct fracture patterns. By implementing a real-time monitoring system with multiple pressure points, the study elucidated the impacts of varying fracture penetration levels and spacing on dynamic pressure wave patterns, mobilization extent, and development outcomes during CO2 huff-n-puff. The results demonstrate that fractures substantially improve CO₂ huff-n-puff efficiency in shale reservoirs compared to a seamless model, with Model 6 achieving the highest recovery rate of 42.38%. As fracture penetration increases, CO₂ dynamic coverage expands. However, full penetration leads to preferential flow through fractures, reducing matrix contact and diminishing huff-n-puff efficiency. Reduced fracture spacing enlarges the well's control domain, amplifies CO2 huff-n-puff dynamic coverage area, and increases oil production and recovery rate. Conversely, decreasing fracture spacing accelerates oil and gas displacement rates, leading to higher economic expenses. The study concludes that 2/3 fracture penetration in five fractures offers the most efficient recovery method. For a 500 m horizontal well, optimal fracture spacing is recommended at 125 m, with a fracture length of 145.8 m for a well spacing of 417 m.

Experimental Study of the Random Wave-Induced Hydrodynamics and Soil Response in a Porous Seabed Around Double Piles

The evaluation of the wave-induced pore pressures around the offshore piles has attracted great attentions among coastal engineers, because they have been commonly used as foundations of numerous marine infrastructures. This paper presents comparative studies of the random wave-induced transient seabed response around single and double piles in a sandy seabed through a series of wave flume experiments. The influences of relative spacing ratios, wave incidence angles, and front pile diameters under different random wave parameters on oscillatory pore pressures in the vicinity of double piles are examined. In addition, variations in wave profiles and dynamic wave pressures surrounding single and double piles are quantitatively analyzed. Based on the experimental results, the following conclusions can be drawn: (1) under the influence of random waves, the wave profiles around the double piles exhibit obvious irregularity and nonlinearity; (2) the shielding effect existing in the tandem piles results in lower dynamic wave pressures around the rear pile compared to the front pile; (3) the pore pressures on the front surface of the double piles decrease with increasing soil depth, with a decreasing attenuation rate at each layer; (4) when the relative spacing ratio G/D2=3, the group-pile effect weakens, leading to an increase in the pore pressures around the rear pile, approaching the results of a single pile under conditions of lower significant wave heights or periods; (5) the intense disturbance effect caused by large wave incidence angles exacerbates the pore pressure response around the double piles; (6) when the diameter of the front pile in the tandem piles increases, it enhances the shielding effect, thus suppressing the seabed response around the rear pile. In contrast, it causes an increase in the wave surface around the double piles, exacerbating the pore pressure response in the seabed. The latter effect becomes more pronounced when the significant wave height is larger.

Separate effect experimental test and theoretical investigate of pressure wave propagation and pressure rising process upon SGTR accident in advanced LFRs

Steam Generator Tube Rupture accident (SGTR) accident is one of the main safety events to the industrial application of Lead-cooled Fast Reactors (LFRs). This paper presents a separate effect test facility, called Lead-bismuth Eutectic Steam generator tube rupture Test facility (LEST), for testing high-pressure subcooled water jetting into a Lead Bismuth Eutectic (LBE) pool, and analyzes experimental data of maximum initial pressure peak, dynamic pressure attenuation together with the pressure rising process upon SGTR. The experimental analyses put forward a general equation form together with a dimensionless number of the initial pressure peak associated with a high-pressure Newtonian fluid injecting into another low-pressure Newtonian fluid, and the derivation of the general equation form is specifically applied to match a further semi-empirical equation for an initial maximum pressure peak with 5–10 MPa water injection in SGTR. The attenuation curve of dynamic pressure wave in each experiment is matched, and more weakening conditions can change the pressure peak from concentrated and non-period peak group to approximate period 1.4–3.7 s pulse group. The size relationship between pipe length L and nozzle diameter D, which L/D=12, is set to divide SGTR into small wall rupture and tube shear fracture, further three stages of pressure rising are distinguished in the tube shear fracture. This work provides experimental data to help determining the severity of SGTR and assessing the consequences of this event as part of the safety case for the industrial application of LFRs.

A top-mounted, chain-driven, bidirectional point absorber type wave energy convertor

The aims of this paper are to introduce the structure of a point absorber type wave energy convertor (WEC) mounted above sea level, to analyze the motion characteristics and power absorptions of such a WEC and to describe the numerous advantages of installing it in the nearshore. The motion equation is solved by representing the excitation force with the dynamic wave pressure and obtaining an approximate system transfer function. The possibility of the buoy becoming airborne or submerged is considered and the power absorption is increased optimally by reshaping the frequency response using a flywheel. The time, phase and spectral averages of power absorptions are presented for offshore waves, nearshore waves and for concentrated nearshore waves using standard wave spectrums and compared. Approximate gains possible with linear and rectangular arrays are also presented. Considering different phase distributions, it is shown that downgrading WECs based on short term power absorptions is not justifiable. With a buoy of radius 2 m, the power absorption of this kind of a WEC is estimated to be about 114 kW in the deep sea and 117 kW in the nearshore, when the energy flux is about 31 kW/m. By concentrating just 13.6 kW/m nearshore waves, it could be increased to 186 kW. If a 3-element longshore in-line array is placed inside the wave concentrator, the total power absorption could be about 556 kW. A 9-element rectangular array would be capable of absorbing about 1 MW by making use of the remaining energy of the down wave.

Ocena ryzyka przypadkowej inicjacji półfabrykatów zawierających artykuły pirotechniczne o przeznaczeniu technicznym w branży motoryzacyjnej

The identification of the danger of accidental initiation, respectively the establishment of the causes and the possibilities that can generate the triggering of pyrotechnic devices with a technical destination in the automotive field containing pyrotechnic articles of the P1 category, is carried out for each piece of equipment depending on the phases of the technological process of their manufacture. The protection against the accidental initiation of pyrotechnic devices intended for equipping auto vehicles is of particular interest for occupational safety in the manufacturing process of these types of products, because their triggering can endanger the life and health of workers as a result of the uncontrollable effects of functioning, with the generation of thermal and dynamic effects, as well as the emission of toxic reaction products that can affect the human component and/or destroy the work space. The paper highlights the way to analyse the risk of accidental initiation of pyrotechnic devices actuator type from the P1 category, using the latest generation technical-scientific tools in order to computerize the effects following their untimely triggering and establish the safety distances in relation to the amplitude of the degree of damage to the human component, respectively the destruction caused to the work space. In the case of pyrotechnic articles from category P1 type pyrotechnic actuator, the main risks of their accidental initiation are determined by the following technical factors: the production of an electrical discharge greater than 25 KV, the generation of a current with an intensity greater than 0.4 A at the terminals pyrotechnic device, their exposure to a temperature higher than 165°C. The main effects generated after the untimely initiations of these products are determined by: flame, thermal radiation, dynamic pressure waves, projected fragments and hazardous releases of chemical combustion products. Depending on the results obtained following the estimation and assessment of the assessed risks, preventive and countermeasures of a technical and organizational nature are established, in order to secure the predicTab. operations of specific operations with these types of products.

Stability Analysis of Cofferdam with Double-Wall Steel Sheet Piles under Wave Action from Storm Surges

Double-wall steel sheet piles (DSSPs) are widely used in large-span cofferdams for docks due to their good performance against wave action during storm surges. This paper describes a study of the dynamic behavior of a DSSP cofferdam under wave action through flume tests and a numerical simulation that combined computational fluid dynamics (CFD) and the finite element method. The influences of the water level and wave height on the DSSP cofferdam were investigated experimentally and numerically. Tall waves in shallow water broke upon and impacted the seaside pile with large dynamic wave pressure, dramatically increasing the stress and displacement of the seaside pile. The overlap of the traveling and reflected waves increased the excess pore water pressure near the seaside pile due to taller overlapped waves and higher wave frequency. The DSSP cofferdam failed under the combined actions of the dynamic wave pressure and erosion of the landside seabed. The leakage and overflow of the breaking waves resulted in significant erosion of the landside seabed and greatly weakened the support of the seabed. The dynamic wave pressure then pushed the DSSP cofferdam until it failed. The simulation with the combined methods of CFD and FEM resulted in trends that were similar to those of the test measurements. Compared to the quasi-static method and pseudo-dynamic method, the results of the simulation via the present method were much closer to the test results because the simulation included the effects of breaking waves. The reinforced measure worked well to prevent the DSSP cofferdam in a sandy seabed foundation from continuous failures of deformation–leakage–erosion–tilting. However, it failed in a clay interlayer seabed foundation due to the large settlement.

Numerical simulations of wave-soil-structure interactions for a subsea cover

Glass-reinforced plastic (GRP) subsea protection covers are widely used to prevent damage to offshore pipelines placed on the seabed from dropped objects, hydrodynamic wave induced loads, and trawling. The light GRP subsea covers could be stabilized by using skirts which penetrate the seabed soil. The dynamic wave pressure acting on the cover will transfer to the cover bottom and the skirt and further influence the pore pressure and seepage flow inside the soil beneath the cover. The present study performs a numerical analysis for the wave-soil-structure interaction (WSSI) of a subsea cover. Two-dimensional (2D) numerical simulations are carried out using an open-source numerical toolbox for modeling the porous seabed interaction with waves and structures under the framework of the finite-volume-method (FVM) based OpenFOAM. The nonlinear waves are solved to obtain the dynamic wave loadings on the cover and the pressure on the seabed. A soil consolidation model is used to provide the initial effective stress in the soil. Then, a one-way coupling algorithm is applied for the WSSI analysis to obtain the soil response in the vicinity of the cover. The distributions of the wave-induced pore pressure, the soil shear stress, and the seepage flow within the seabed are studied and the influences of the wave heights and the skirt lengths are discussed.

3D numerical investigation of wave-induced seabed response around a tripod pile foundation with different arrangements

This paper developed an integrated numerical model for the dynamic response of the seabed around tripod turbine foundations under dynamic wave pressure. Three-dimensional (3D) numerical analysis is performed by applying an integrated multiphysics model developed in the finite volume method (FVM) based on the OpenFOAM framework. The present work adopts nonlinear wave model, linear elastic structure model, and anisotropic porous seabed model for wave-structure-seabed interactions. Accuracy of the wave-structure-seabed coupling process was ensured by comparing the numerical model with the numerical and experimental results. This study indicates: (1) the existence of tripods will significantly affect the consolidation results of the surrounding seabed; (2) wave force acting on the tripod structure varies with different arrangements, in which the maximum wave force under the wave crest decreases as the angle between wave and structure increases; and (3) the asymmetry 15° arrangement would significantly increase the maximum liquefaction depth of the seabed.

Water Hammer

Dynamic Pressure rise or fall in all kinds of pipelines due to the operation of the valve which are transmitting fluids are very much essential part to consider in the design of pipeline. This paper presents how to calculate these dynamic pressure waves while designing pipelines.

Numerical simulation of interaction between wave-driven currents and revetment on coral reefs

The stability of a revetment breakwater is vital for the safety of islands and reefs formed with reclaimed sand under the action of a complex marine environment. The interactions in wave-reef-revetment systems are complex and difficult to accurately simulate; in addition, it is difficult to describe the dynamic wave pressures distribution by empirical formula, although these parameters are important for the safety of reefs and revetments; hence, more studies are required. In this study, an investigation of the interaction between the wave-driven current on the reef flat and rear revetment breakwater is performed, and the hydrodynamic pressure on the revetment is discussed. A two-dimensional numerical model is established using the Navier–Stokes equation to explore the interaction and calculate the pressure distribution on the revetment, which is verified by an experiment designed and implemented in a wave flume. The simulation illustrates that the wave-driven current is stratified in the vertical direction, and the velocity at the surface is approximately 2–3 times that of the average value, and the maximum flow velocity can reach 13 m/s. By varying the parameters of the geomorphology of the coral reef, the results suggest that for steeper fore-reef slopes, longer reef flats, and larger relative distances, the peak value of the pressure increases consistently. The detailed and sophisticated hydrodynamic pressure distribution on the revetment is obtained, and the numerical results indicate that the fifth-order polynomial can match the regulation of pulsating pressure at the wave surface with the relative submerged water depth.

Experimental Study of Forces Influencing Vertical Breakwater under Extreme Waves

In order to understand the extreme wave acting on the vertical breakwater, a series of experiments were constructed in the wave tank to measure the variations of pressure on the front, rear faces, and below the caisson due to overtopping waves. The front and backward horizontal forces and the uplift forces were estimated by integrating the dynamic wave pressure distributions. The COBRAS numerical model was also used to calculate the wave loads under various overtopping waves. The measured wave pressures and wave forces were compared with the predictions of numerical results and showed good agreement. It was found that the forces acting on the backward side of the vertical structure induced by the wave overtopping should be considered. From the experimental data, new semi-empirical equations for calculating the maximum wave forces are proposed using a least squares approximation.

Assessment of dynamic pressure and wave forces on vertical-caisson type breakwater

The design and construction of coastal structures such as breakwaters, at great water depths is rapidly increasing as a result of the increasing draught of large vessels and off-shore land reclamations. Vertical caisson breakwaters may be the best alternative compared to ordinary rubble mound breakwaters in larger water depths, in terms of performance, total costs, environmental aspects, construction time and maintenance. To fulfilling the functional utility and impact of the structure on the sea environment, it is necessary to study the hydraulic performance. This can be found by field investigation, numerical simulations and by physical modelling. Scale modelling techniques are used to study various coastal engineering problems. This article presents the results obtained by conducting series of experiments in two-dimensional wave flume to assess the hydrodynamic performance of vertical-caisson breakwater, which is made of concrete, with the protection of toe. The dynamic pressure distribution, wave runup, wave reflection, wave forces and stability parameter on the vertical caisson breakwater are discussed. The maximum wave force on the wall breakwater is calculated from measured pressure values and is compared with the forces calculated by Goda’s and Sainflou wave theories.

Research on Gas Drainage Laws of SRV-fractured Horizontal Wells in Different Types of Tight Gas Reservoirs

Horizontal well and large-scale hydraulic fracturing are the main technical means for tight gas development at present. However, tight gas reservoirs are characterized by strong heterogeneity, and different types of gas Wells have different production characteristics and pressure propagation laws. In view of the above problems, based on the physical characteristics and reservoir reconstruction evaluation of tight gas reservoir in GAS Reservoir A in Sichuan Basin, three types of reservoir hydraulic fracture models and numerical simulation models were generated by using numerical simulation method, and the production dynamic prediction and pressure wave propagation law were studied. The simulation results show that: (1) In the beginning of production, the deflated area is mainly hydraulic fracture reconstruction area, and the drainage radius of all kinds of gas Wells expands at similar speed; (2) When the drainage range is expanded in the untransformed area, the dimensionless drainage radius expansion speed of gas well is positively correlated with the reservoir physical property. The dimensionless drainage radius and drainage area of class I and II horizontal Wells are relatively faster, while that of Class III Wells is slightly slower.

A novel design method for wave-induced fatigue of flood gates

This paper presents a novel design method to predict fatigue of flood gates due to dynamic wave loading. The accumulation of fatigue damage is predicted probabilistically over the entire lifetime of the structure rather than with a set of normative events. Load events are defined using a joint probability distribution of historical wind and water level data. The random phase-amplitude model is employed to obtain realisations of the wave state for every combination of environmental conditions. Linear wave theory and pressure-impulse theory are used to predict both quasi-steady and highly dynamic wave pressures. The stress response of the structure is predicted with a hybrid semi-analytical and finite element model. By applying a mode matching technique the fluid-structure interaction is solved in a computationally efficient manner. This facilitates the large number of simulations required for a comprehensive fatigue analysis without making concessions in the physical modelling. The fatigue damage is then evaluated with the linear Palmgren-Miner method by applying a rainflow algorithm. A Monte Carlo analysis is performed to estimate the expected fatigue lifetime of the structure. The modular structure of the model routine allows for easy adaptation to other situations where fatigue due to hydrodynamic loading is of interest. The design method is applied to a case study of a flood gate with an overhang inspired by the situation at the Afsluitdijk. Non-fundamental modes are taken into account without simplification of the fluid-structure interaction process and found to be governing for the fatigue damage for the studied case. Moreover, the interference of vibrations due to consecutive wave impacts is shown to have a significant influence on the outcome of the fatigue assessment. For the case study, the design method leads to a 10-20% reduction of the governing fatigue damage compared to a method commonly used in practice. At specific locations on the flood gate fatigue damage is found to be underestimated by current design methods. The presented design method is therefore found to be a significant improvement.

Upgraded Fiber-Optic Sensor System for Dynamic Strain Measurement in Spallation Neutron Source

We describe an upgraded fiber-optic sensor system and its performance in measuring the dynamic strains in a mercury target of the Spallation Neutron Source (SNS). Strains result from dynamic pressure waves in the stainless-steel mercury target induced by short (~700 ns), intense (up to 23.3 kJ), high-energy (~1 GeV) proton pulses. In the upgraded sensor system, the output of each sensor head is interrogated with a compact, all-fiber based Faraday Michelson interferometer, which generates interference signals with a steady phase shift. Strain waveforms are recovered from the phase-shifted interference signals using a high-speed digital signal processing procedure developed in our previous work. We demonstrate successful measurements of dynamic strain pulses, e.g., <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$400 \mu \varepsilon $ </tex-math></inline-formula> over <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$190 \mu \text{s}$ </tex-math></inline-formula> , on a recently installed SNS target using the upgraded sensor system. The measured strain waveforms are analyzed for more than 20 sensor locations and/or orientations, and provide information regarding the temporal structure of strain profiles and dependence of the strain magnitude on the proton powers of 200 – 1400 kW. The new interrogator also measures the radiation-induced-attenuation (RIA) in the optical fiber, enabling experimental investigations of RIA profiles induced by a 700-ns radiation pulse. The radiation effects on the strain measurement performance are discussed over a radiation dose range of up to <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">${4}\times {10}^{{8}}$ </tex-math></inline-formula> Gy and an RIA compensation method is proposed. The measurements allow insight into the response of this unique piece of equipment and can be used for validation of simulations.

Experimental study on the amplitude characteristics and propagation velocity of dynamic pressure wave for the leakage of gas-liquid two-phase intermittent flow in pipelines

Intermittent flow is one of the most complex flow patterns in gas-liquid two-phase flow in pipelines. The leakage of pipeline intermittent flow poses a threat to the operating safety. In order to improve the application of acoustic method to gas-liquid two-phase leakage detection, the amplitude characteristics and propagation velocity of dynamic pressure wave of intermittent flow for gas-liquid two-phase pipeline leakage were investigated. Experiments with conditions of superficial liquid velocity (vsL) from 0.57 m/s to 1.06 m/s, superficial gas velocity (vsG) from 0.3 m/s to 4.47 m/s, leakage aperture from 4 mm to 9 mm were carried out, and the characteristics were studied by time-domain analysis, frequency-domain analysis and time-frequency analysis. The propagation velocity formula of dynamic pressure wave in two-phase medium was fitted based on Levenberg-Marquardt algorithm, and the propagation velocity of straight pipe and elbow pipe was discussed. The experimental results indicate that the peak-to-peak value of dynamic pressure wave measured by sensor within 4.25 m increases during top or mid sustained leakage when leak aperture is larger than 4 mm in elongated bubble flow and 6 mm in slug flow. The amplitude of sustained leakage dynamic pressure increases as gas flow rate increases. The low frequency energy of signal increases immediately when leakage occurs, the energy of signal mainly concentrates at 0–15 Hz, the frequency amplitude of sustained leakage signal is much higher than normal flow signal and leak occur signal. As the gas volume fraction β increases from 0.23 to 0.75, the dynamic pressure propagation velocity decreases from 33.13 m/s to 26.5 m/s and then increases from 26.5 m/s to 66.25 m/s. The propagation velocity increases when it flows through an elbow pipe section. The results are meaningful for further study of two-phase pipeline slug flow leakage detection and location.

Sound–turbulence interaction model for low mach number flows and its application in natural gas pipeline leak location

When leakage dynamic pressure waves (DPWs) propagate in low Mach number flows, the viscothermal effects are considered the main reason for sound attenuation. However, an experimental analysis conducted in this study shows that the non-equilibrium sound–turbulence interaction process is the main cause. The turbulence effects due to turbulent flows act on the DPWs, and the fluctuations due to the DPWs act on the turbulent flows. Both processes result in the turbulent absorption of the gas to the amplitude of the DPWs, leading to amplitude attenuation at sufficiently low frequencies. To predict the amplitude attenuation, a non-equilibrium sound–turbulence interaction model is established, solved, and verified using analytical and experimental results, which show that attenuation coefficients (ACs) obtained by considering the sound–turbulence interaction effects are 1.6–3.5 times larger than those obtained by only considering the viscothermal effects, even when the Mach number is between 0.0038 and 0.016. The established model can improve leak localization.

Cavitation onset caused by a dynamic pressure wave in liquid pipelines

When liquids flow in the pipelines, the onset of cavitation can be characterized by a variant of the Euler number known as the cavitation number (CN), which is based on the velocity and denoted by C in this paper. Conventionally, cavitation is considered to be induced when C ~ 1. However, experimental observations and several pipe bursts indicate that the CN may incorrectly predict the onset of cavitation. For example, when leakage occurs in the pipeline or a valve in the pipeline is opened, the resultant pressure loss generates a dynamic pressure wave with a small amplitude, which may lead to bubble formation, even though C ~ 1 is not satisfied. Hence, this paper proposes another CN based on the amplitude of the generated dynamic pressure wave, rather than the velocity, for ascertaining the onset of cavitation. The validity of the proposed CN was verified through experiments and a case study. The results indicated that the proposed CN can be effectively used for cavitation prediction induced by pressure fluctuations and for investigating phenomena such as pressure fluctuation, leakage, and corrosion in liquid pipelines, tanks, and pressure vessels, as well as the safety design of liquefied natural gas tanks and tankers.

An integrated numerical model for the stability of artificial submarine slope under wave load

The stability evaluation of artificial submarine slope under wave load is a primary concern in the construction of offshore structures such as submarine pipelines and immersed tunnels. Aiming at this issue, an integrated Finite Element Method (FEM) model of wave-seabed interaction is proposed to investigate the effects of wave action on the wave-induced slope stability of the temporary trench. The Reynolds-Averaged Navier-Stokes (RANS) equations are adopted to describe the wave motion in the fluid domain, while the seabed domain is described using the Mohr-Coulomb constitutive model. The wave-induced seepage pressure near the slope is determined by Darcy's law and considered as an external load in the analysis. The stability of the temporary slope is evaluated by the means of shear strength reduction. The present numerical model is firstly validated through the comparison with the analytical solution and the experimental data available in the literature. The integrated numerical model is then applied to a practical engineering application, including a new evaluation method of stability for the trench slope excavated temporarily for the fully-buried submarine pipelines at an actual site. The results indicate that the stability of submarine slope is closely related to the dynamic wave pressure from the upper seawater. The minimum of FOS (factor of safety) corresponding to the most dangerous situation is introduced as the stability index for the slope with specific slope ratio in the whole process of dynamic wave loading. Additionally, the wave-induced seepage pressure should be considered as external load in the design for the seabed with large permeability such as sandy soils.