Forward prediction for tunnel geology and classification of surrounding rock based on seismic wave velocity layered tomography

Excavation under complex geological conditions requires effective and accurate geological forward-prospecting to detect the unfavorable geological structure and estimate the classification of surround-ing rock in front of the tunnel face.In this work,a forward-prediction method for tunnel geology and classification of surrounding rock is developed based on seismic wave velocity layered tomography.In particular,for the problem of strong multi-solution of wave velocity inversion caused by few ray paths in the narrow space of the tunnel,a layered inversion based on regularization is proposed.By reducing the inversion area of each iteration step and applying straight-line interface assumption,the convergence and accuracy of wave velocity inversion are effectively improved.Furthermore,a surrounding rock classification network based on autoencoder is constructed.The mapping relationship between wave velocity and classification of surrounding rock is established with density,Poisson’s ratio and elastic modulus as links.Two numerical examples with geological conditions similar to that in the field tunnel and a field case study in an urban subway tunnel verify the potential of the proposed method for practical application.

Effect of slope gradient on the subsurface water flow velocity of sand layer profile

Subsurface water flow velocity influences the hydrodynamic characteristics of soil seepage and the interaction between subsurface water flow and surface runoff during soil erosion and sediment transport.A visualized method and equipment was adopted in this study to observe the subsurface water flow.Quartz sand was used as the test material of subsurface water flow and fluorescent dye was used as the indicator for tracing subsurface water flow.Water was supplied at the same flow discharge to the three parts at the bottom of the test flume,and the subsurface water flow were determined with four slope gradients(4°,8°,10°,and 12°).The results showed that the seepage velocity gradually increased with increasing slope gradient.The pore water velocity at different depths of sand layer profile increased with increasing slope gradient,whereas the thickness of the flow front gradually decreased.For the same slope gradient,the pore water velocity in the lower layer was the largest,whereas the thickness of the flow front was the smallest.Comparative analysis of the relationship between seepage velocity and pore water velocity at different depths of sand layer profile showed that the maximum relative difference between the measured pore water velocity and the computational pore water velocity at different depths of sand profile in the experiment was 4.38%.Thus,the test method for measuring the subsurface water flow velocity of sand layer profile adopted in this study was effective and feasible.The development of this experiment and the exploration of research methods would lay a good test foundation for future studies on the variation law of subsurface water flow velocity and the determination of flow velocity in purple soils,thus contributing to the improvement of the hydrodynamic mechanism of purple soils.

Experimental studies on flow visualization and velocity field of compression ramp with different incoming boundary layers

Experimental studies which focus on flow visualization and the velocity field of a supersonic laminar/turbulent flow over a compression ramp were carried out in a Mach 3.0 wind tunnel. Fine flow structures and velocity field structures were obtained via NPLS （nanoparticle-tracer planar laser scattering） and PIV （particle image velocimetry） techniques, time- averaged flow structures were researched, and spatiotemporal evolutions of transient flow structures were analyzed. The flow visualization results indicated that when the ramp angles were 25~, a typical separation occurred in the laminar flow, some typical flow structures such as shock induced by the boundary layer, separation shock, reversed flow and reattachment shock were visible clearly. While a certain extent separation occurred in turbulent flow, the separation region was much smaller. When the ramp angles were 28~, laminar flow separated further, and the separation region expanded evidently, flow structures in the separation region were complex. While a typical separation occurred in turbulent flow, reversed flow structures were significant, flow structures in the separation region were relatively simple. The experimental results of velocity field were corresponding to flow visualization, and the velocity field structures of both compression ramp flows agreed with the flow structures well. There were three layered structures in the U component velocity, and the V component velocity appeared like an oblique ＂v＂. Some differences between these two compression ramp flows can be observed in the velocity profiles of the shear layer and the shearing intensity.

Study on the double-logarithmic profile of tidal flow velocity in the near-bed layers

Tidal current velocity profile in the near-bed layers has been widely studied. The results showed that velocity profile in the near-bed layer obviously departure from the traditional logarithmic profile, due to the acceleration or deceleration. Although the logarithmic linear profile can reduce the rate of deviation from this, only it is a lower-order approximate solution. In this paper, considering the unsteady and non-linear features of tidal motion, the double logarithmic profile near-bed layers in estuarine and coastal waters is established on the assumption that the turbulent shear stress along the water depth was parabolic distribution, and on the basis of Prandtl＇s mixing length theory and yon Karman＇s self-similar theory. Having been verified the data observed at the West Solent in the south of England, and comparison of the logarithmic linear profile, it found that the double logarithmic profile is more precious than the latter. At last, the discussed results showed that： （1） The parabolic distribution of the tidal shear stresses verified good by the field data and experimental data, can be better reflected the basic features of the tidal shear stress deviating from linear distribution that is downward when to accelerate, upward when to decelerate. （2） The traditional logarithmic velocity profile is the zero-order approximation solution of the double logarithmic profile, the logarithmic linear profile is the first order, and the logarithmic parabolic profile is the second order. （3） Ignoring the conditions of diffusion and convection in the tida movement, the double logarithmic profile can reflect the tidal properties of acceleration or deceleration, so that the calculation of the friction velocity and roughness length are more reasonable. When the acceleration or the deceleration is about zero, the double logarithmic profile becomes the logarithmic profile.

Lateral Shear Layer and Its Velocity Distribution of Flow in Rectangular Open Channels

The lateral velocity distribution of flow in the shear layer of open channel is required to many problems in river and eco-environment engineering, e.g. distribution of pollutant dispersion, sediment transport and bank erosion, and aquatic habitat. It is not well understood about how the velocity varies laterally in the wall boundary layer. This paper gives an analytical solution of lateral velocity distribution in a rectangular open channel based on the depth-averaged momentum equation proposed by Shiono & Knight. The obtained lateral velocity distributions in the wall shear layer are related to the two hydraulic parameters of lateral eddy viscosity (λ) and depth-averaged secondary flow (Γ) for given roughened channels. Preliminary relationships between the above two parameters and the aspect ratio of channel, B/H, are obtained from two sets of experimental data. The lateral width (δ) of the shear layer was investigated and found to relate to the λ and the bed friction factor (f), as described by Equation (26). This study indicates that the lateral shear layer near the wall can be very wide (δ/H = 14.6) for the extreme case (λ = 0.6 and f = 0.01).

Numerical investigation of velocity distribution of turbulent flow through vertically double-layered vegetation

The velocity structures of flow through vertically double-layered vegetation(VDLV)as well as single-layered rigid vegetation(SLV)were investigated computationally with a three-dimensional(3D)Reynolds stress turbulence model,using the computational fluid dynamics(CFD)code FLUENT.The detailed velocity distribution was explored with a varying initial Froude number(Fr),with consideration of the steady subcritical flow conditions of an inland tsunami.In VDLV flows,the numerical model successfully captured the inflection point in the profiles of mean streamwise velocities in the mixing-layer region around the top of short submerged vegetation.An upward and downward movement of flow occurred at the positions located just behind the tall and short vegetation,respectively.Overall,higher streamwise velocities were observed in the upper vegetation layer due to high porosity,with Pr=98%(sparse vegetation,where Pr is the porosity),as compared to those in the lower vegetation layer,which had comparatively low porosity,with Pr=91%(dense vegetation).A rising trend of velocities was found as the flow passed through the vegetation region,followed by a clear sawtooth distribution,as compared to the regions just upstream and downstream of vegetation where the flow was almost uniform.In VDLV flows,a rising trend in the flow resistance was observed with the increase in the initial Froude number,i.e.,Fr?0.67,0.70,and 0.73.However,the flow resistance in the case of SLV was relatively very low.The numerical results also show the flow structures within the vicinity of short and tall vegetation,which are difficult to attain through experimental measurements.

Ultrasonic method for measuring water holdup of low velocity and high-water-cut oil-water two-phase flow

Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration（of oil droplets） in oil-water two-phase flow,which makes it difficult to measure water holdup in oil wells.In this study,we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in lowvelocity and high water-cut conditions.First,we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling.Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor.Based on the results,we then investigate the effects of oildroplet diameter and distribution on the ultrasonic field.To further understand the measurement characteristics of the ultrasonic sensor,we perform a flow loop test on vertical upward oilwater two-phase flow and measure the responses of the optimized ultrasonic sensor.The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow（D OS/W flow）,but the resolution is favorable for dispersed oil in water flow（D O/W flow） and very fine dispersed oil in water flow（VFD O/W flow）.This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.

Low-velocity Impact Damage Analysis of Composite Laminates Using Self-adapting Delamination Element Method

On the basis ofa 2D 4-node Mindlin shell element method, a novel self-adapting delamination finite element method is presented, which is developed to model the delamination damage of composite laminates. In the method, the sublaminate elements are generated automatically when the delamination damage occurs or extends. Thus, the complex process and state of delamination damage can be simulated practically with high efficiency for both analysis and modeling. Based on the self-adapting delamination method, linear dynamic finite element damage analysis is performed to simulate the low-velocity impact damage process of three types of mixed woven composite laminates. Taking the frictional force among sublaminations during delaminating and the transverse normal stress into account, the analytical results are consistent with those of the experimental data.

VELOCITIES OF PRECURSOR SOLITON GENERATION IN SINGLE LAYER FLOW

In this study, the moving velocitiy of precursor solitons, of the flow in depressed region, and of the zero-crossing of the trailing wavetrain relative to the moving disturbance for single-layer flow over topography were theorecticaly determined in terms of the mass and energy conservation theorems, and were examined with numerical calculations showing good agreement with theoretical results.

A Laboratory Modeling of the Velocity Field in the Convective Boundary Layer with the Particle Image Velocimetry Technique

Based on the research of the convective boundary layer (CBL) temperature field in a convective tank, this paper studies the characteristics of the CBL velocity field in the convective tank. Aluminium powder (400 orders) is used as a tracer particle in the application of the particle image velocimetry (PIV) technique. The experiment demonstrates: the velocity distribution in the mixed layer clearly possesses the characteristics of CBL thermals; the velocity distribution in the top zone of the mixed layer shows entrainment layer characteristics; the vertical distribution of turbulent characteristic variables is reasonable, which is similar to field observations and other tank results; the error analysis demonstrates the validity of aluminium powder, which implies the reliability of the results.

A new analytical model for the low-velocity perforation of thin steel plates by hemispherical-nosed projectiles

Ballistic experiments were conducted on thin steel plates that are normally impacted by hemisphericalnosed projectiles at velocities higher than their ballistic limits. The deformation and failure modes of the thin steel plates were analyzed. A new method was proposed according to the experimental results and the perforation phenomenon of the thin steel plates to determine the radius of the bulging region. In establishing this new method, a dynamic method combined with the plastic wave propagation concept based on the rigid plastic assumption was adopted. The whole perforation process was divided into four consecutive stages, namely, bulging deformation, dishing deformation, ductile hole enlargement, and projectile exit. On the basis of the energy conservation principle, a new model was developed to predict the residual velocities of hemispherical-nosed projectiles that perforate thin steel plates at low velocities.The results obtained from the theoretical calculations by the present model were compared with the experimental results. Theoretical predictions were in good agreement with the experimental results in terms of both the radius of the bulging region and the residual velocity of the projectile when the strain rate effects of the target material during each stage were considered.

GROUP VELOCITY CONTROL SCHEME WITH LOW DISSIPATION

In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.

A Laboratory Study of the Turbulent Velocity Characteristics in the Convective Boundary Layer

Based on the measurement of the velocity field in the convective boundary layer （CBL） in a convection water tank with the particle image velocimetry （PIV） technique, this paper studies the characteristics of the CBL turbulent velocity in a modified convection tank. The experiment results show that the velocity distribution in the mixed layer clearly possesses the characteristics of the CBL thermals, and the turbulent eddies can be seen obviously. The comparison of the vertical distribution of the turbulent velocity variables indicates that the modeling in the new tank is better than in the old one. The experiment data show that the thermal＇s motion in the entrainment zone sometimes fluctuates obviously due to the intermittence of turbulence. Analyses show that this fluctuation can influence the agreement of the measurement data with the parameterization scheme, in which the convective Richardson number is used to characterize the entrainment zone depth. The normalized square velocity wi^2/w＊^2. at the top of the mixed layer seems to be time-dependent, and has a decreasing trend during the experiments. This implies that the vertical turbulent velocity at the top of the mixed layer may not be proportional to the convective velocity （w＊）.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers:A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.

The failure propagation of weakly stable sediment:A reason for the formation of high-velocity turbidity currents in submarine canyons

The long-distance movement of turbidity currents in submarine canyons can transport large amounts of sediment to deep-sea plains.Previous studies show obvious differences in the turbidity current velocities derived from the multiple cables damage events ranging from 5.9 to 28.0 m/s and those of field observations between 0.15 and 7.2 m/s.Therefore,questions remain regarding whether a turbid fluid in an undersea environment can flow through a submarine canyon for a long distance at a high speed.A new model based on weakly stable sediment is proposed(proposed failure propagation model for weakly stable sediments,WS S-PFP model for short)to explain the high-speed and long-range motion of turbidity currents in submarine canyons through the combination of laboratory tests and numerical analogs.The model is based on two mechanisms:1)the original turbidity current triggers the destabilization of the weakly stable sediment bed and promotes the destabilization and transport of the soft sediment in the downstream direction and 2)the excitation wave that forms when the original turbidity current moves into the canyon leads to the destabilization and transport of the weakly stable sediment in the downstream direction.The proposed model will provide dynamic process interpretation for the study of deep-sea deposition,pollutant transport,and optical cable damage.

Tropical Cyclones and Polar Lows: Velocity, Size, and Energy Scales, and Relation to the 26℃ Cyclone Origin Criteria

The goal of this paper is to quantitatively formulate some necessary conditions for the development of intense atmospheric vortices. Specifically, these criteria are discussed for tropical cyclones （TC） and polar lows （PL） by using bulk formulas for fluxes of momentum, sensible heating, and latent heating between the ocean and the atmosphere. The velocity scale is used in two forms： （1） as expressed through the buoyancy flux b and the Coriolis parameter Ic for rotating fluids convection, and （2） as expressed with the cube of velocity times the drag coefficient through the formula for total kinetic energy dissipation in the atmospheric boundary layer. In the quasistationary case the dissipation equals the generation of the energy. In both cases the velocity scale can be expressed through temperature and humidity differences between the ocean and the atmosphere in terms of the reduced gravity, and both forms produce quite comparable velocity scales. Using parameters b and Ic, we can form scales of the area and, by adding the mass of a unit air column, a scale of the total kinetic energy as well. These scales nicely explain the much smaller size of a PL, as compared to a TC, and the total kinetic energy of a TC is of the order 1018 - 1019 J. It will be shown that wind of 33 m s^-1 is produced when the total enthalpy fluxes between the ocean and the atmosphere are about 700 W m-2 for a TC and 1700 W m-2 for a PL, in association with the much larger role of the latent heat in the first case and the stricter geostrophic constraints and larger static stability in the second case. This replaces the mystical role of 26℃ as a criterion for TC origin. The buoyancy flux, a product of the reduced gravity and the wind speed, together with the atmospheric static stability, determines the rate of the penetrating convection. It is known from the observations that the formation time for a PL reaching an altitude of 5-6 km can be only a few hours, and a day, or even half a day, for a TC reaching 15-18 km. These two facts allow us to construct curves on the plane of Ts and ΔT = Ts - Ta to determine possibilities for forming an intense vortex. Here, Ta is the atmospheric temperature at the height z = 10 m. A PL should have AT 〉 20℃ in accordance with the observations and nmnerical simulations. The conditions for a TC are not so straightforward but our diagram shows that the temperature difference of a few degrees, or possibly even a fraction of a degree, might be sufficient for TC development for a range of static stabilities and development times.

Effect of SACMA and QMW Test Methods on Compressive Properties of Composite Laminates after Low Velocity Impact

Experimental studies on the compressive behavior of composite laminates after low velocity impact was carried out with two test methods.One is SACMA Standard,and the other is a small dimensional specimen test method.Impact damage distributions,compressive failure process after impact,quasi static indentation and compression of laminates with a hole were brought into comparison between these two test methods.The results showed that there is a great difference between these two test methods.Compressive behavior of laminates after impact varies with different test methods.Residual compressive strength of laminates after low velocity impact measured with SACMA Standard can reflect stiffness properties of composite resins more wholely than that measured with the other method can do.Small dimensional specimen test method should be improved on as an experimental standard of compressive strength after impact.

MULTI-SCALE COHERENT STRUCTURES IN TURBULENT BOUNDARY LAYER DETECTED BY LOCALLY AVERAGED VELOCITY STRUCTURE FUNCTIONS

The time sequence of longitudinal velocity component at different vertical locations in turbulent boundary layer was finely measured in a wind tunnel. The concept of coarse_grained velocity structure functions, which describes the relative motions of straining and compressing for multi_scale eddy structures in turbulent flows, was put forward based on the theory of locally multi_scale average. Based on the consistency between coarse_grained velocity structure function and Harr wavelet transformation,detecting method was presented, by which the coherent structures and their intermittency was identified by multi_scale flatness factor calculated by locally average structure function. Phase_averaged evolution course for multi_scale coherent eddy structures in wall turbulence were extracted by this conditional sampling to educe scheme. The dynamics course of multi_scale coherent eddy structures and their effects on statistics of turbulent flows were studied.

Effect of Cycling Low Velocity Impact on Mechanical and Wear Properties of CFRP Laminate Composites

The mechanical and wear properties of CFRP laminate were investigated using a method of cycling low velocity impact, to study the trend and mechanism of impact resistance of the CFRP laminate under repeated impact during its service process. The interface responses of CFRP laminate under di erent impact kinetic energy during the cycling impact process were studied were studied experimentally, such as impact contact duration, deformation and energy absorption. The worn surface morphologies were observed through optical microscopy and a 3?D surface profiler and the cross?sectional morphologies were observed through SEM to investigate the mechanism of impact material dam?age. Based on a single?degree?of?freedom damping vibration model, the normal contact sti ness and contact damp?ing of the material in di erent wear stages were calculated. It shows the failure process of CFRP laminate damaged by accumulated absorption energy under the cycling impact of di erent initial kinetic energy. The results indicate that the sti ness and damping coe cients will change at di erent impact velocities or cycle numbers. The damage mechanism of CFRP laminates under cycling low kinetic energy is delamination. After repeated experiments, it was found that there was a threshold value for the accumulated absorption energy before the failure of the CFRP laminate.

Revision of P-wave Velocity and Thickness of Hydrate Layer in Shenhu Area, South China Sea

To revise P-wave velocity and thickness of the hydrate layer in the Shenhu area of the South China Sea, acoustic and resistivity logging curves are reanalyzed. The waterlogging phenomenon is found in the shallow sediments of five drilling wells, which causes P-wave velocity to approximate the propagation velocity of sea water(about 1500 m s-1). This also affects the identification of the hydrate layer and results in the underestimate of its thickness. In addition, because there could be about a 5 m thick velocity ramp above or below the hydrate layer as interpreted by acoustic and resistivity logging curves, the recalibrated thickness of this layer is less than the original estimated thickness. The recalibrated P-wave velocity of the hydrate layer is also higher than the original estimated velocity. For the drilling well with a relatively thin hydrate layer, the velocity ramp plays a more important role in identifying and determining the thickness of the layer.