Inverting the rock mass P-wave velocity field ahead of deep buried tunnel face while borehole drilling

Imaging the wave velocity field surrounding a borehole while drilling is a promising and urgently needed approach for extending the exploration range of the borehole point.This paper develops a drilling process detection(DPD)system consisting of a multifunctional sensor and a pilot geophone installed at the top of the drilling rod,geophones at the tunnel face,a laser rangefinder,and an onsite computer.A weighted adjoint-state first arrival travel time tomography method is used to invert the P-wave velocity field of rock mass while borehole drilling.A field experiment in the ongoing construction of a deep buried tunnel in southwestern China demonstrated the DPD system and the tomography method.Time-frequency analysis of typical borehole drilling detection data shows that the impact drilling source is a pulse-like seismic exploration wavelet.A velocity field of the rock mass in a triangular area defined by the borehole trajectory and geophone receiving line can be obtained.Both the borehole core and optical image validate the inverted P-wave velocity field.A numerical simulation of a checkerboard benchmark model is used to test the tomography method.The rapid convergence of the misfits and consistent agreement between the inverted and observed travel times validate the P-wave velocity imaging.

Velocity distribution of the flow field in the cyclonic zone of cyclone-static micro-bubble flotation column

Laboratory experiments have been conducted to study the flow field in a cyclone static micro-bubble flotation column. The method of Particle Image Velocimetry (PIV) was used. The flow field velocity distribution in both cross section and longitudinal section within cyclonic zone was studied for different circulating volumes. The cross sectional vortex was also analyzed. The results show that in cross section as the circulating volume increases from 0.187 to 0.350 m 3 /h, the flow velocity ranges from 0 to 0.68 m/s. The flow field is mainly a non-vortex potential flow that forms a free vortex without outside energy input. In the cyclonic region the vortex deviates from the center of the flotation column because a single tangential opening introduces circulating fluid into the column. The tangential component of the velocity plays a defining role in the cross section. In the longitudinal section the velocity ranges from 0 to 0.08 m/s. The flow velocity increases as does the circulating volume. Advantageous mineral separation conditions arise from the combined effects of cyclonic flow in cross and longitudinal section.

An analysis on short-wave components of the global stress field

The 10 920 stress indicators collected so far by the WSM (World Stress Map) project represent the observed ori-entations of the maximum horizontal principal stress (sHmax) in a certain region. Assuming that the long-wave component of sHmax is expressed by the absolute direction of plate motions, we can get the relative orientation and the magnitude of the short-wave component resulted from the local tectonic process or other factors with vector analytical technique. The global surface was divided into basic element bins by 2.52.5 dimensions and the WSM indicators were statistically analyzed for each element by weight coefficient method in order to determine the mean orientation of the stress. We calculated the long-wave component of the global stress field using HS2-NUVEL1 model. The relative magnitude or the direction limitation of short-wave component, which reflect the local contribution to the observed stresses, was determined by the angle between the mean sHmax and the orien-tation of the long-wave component. The results of this paper show that the contribution of either the long-wave component or the short-wave component is approximately equal to most of the global plates on the basis of the mean effect of the observed stresses. For some of continental regions, the local active tectonics plays an important role in the observed stresses and controls the generation and occurrence of earthquakes.

Effects of nutrient limitations on the sinking velocity of Thalassiosira weissflogii

The sinking of diatoms is critic al to the formation of oceanic biological pumps and coastal hypoxic zones.However,little is known about the effects of different nutrient restrictions on diatom sinking.In this study,we measured the sinking velocity(SV) of Thalassiosira weissflogii using a new phytoplankton video observation instrument and analyzed major biochemical components under varying nutrient conditions.Our results showed that the SV of T.weissflogii under different nutrient limitation conditions varied substantially.The highest SV of(1.77±0.02) m/d was obtained under nitrate limitation,signific antly surpassing that under phosphate limitation at(0.98±0.13) m/d.As the nutrient limitation was released,the SV steadily decreased to(0.32±0.03) m/d and(0.15±0.05) m/d,respectively.Notably;under conditions with limited nitrate and phosphate concentrations,the SV values of T.weissflogii significantly positively correlated with the lipid content(P <0.001),with R^(2) values of 0.86 and 0.69,respectively.The change of the phytoplankton SV was primarily related to the intracellular compo sition,which is controlled by nutrient conditions but did not significantly correlate with transparent extracellular polymer and biosilica contents.The results of this study help to understand the regulation of the vertical sinking process of diatoms by nutrient restriction and provide new insights into phytoplankton dynamics and their relationship with the marine nutrient structure.

Interseismic slip distribution and locking characteristics of the mid-southern segment of the Tanlu fault zone

We employ the block negative dislocation model to invert the distribution of fault coupling and slip rate deficit on the different segments of the Tanlu(Tancheng-Lujiang) fault zone, according to the GPS horizontal velocity field from 1991 to 2007(the first phase) and 2013 to 2018(the second phase). By comparing the deformation characteristics results, we discuss the relationship between the deformation characteristics with the M earthquake in Japan. The results showed that the fault coupling rate of the northern section of Tancheng in the second phase reduced compared with that in the first phase. However, the results of the two phases showed that the northern section of Juxian still has a high coupling rate, a deep blocking depth, and a dextral compressive deficit, which is the enrapture section of the 1668 Tancheng earthquake. At the same time, the area strain results show that the strain rate of the central and eastern regions of the second phase is obviously enhanced compared with that of the first phase. The occurrence of the great earthquake in Japan has played a specific role in alleviating the strain accumulation in the middle and south sections of the Tanlu fault zone. The results of the maximum shear strain show that the shear strain in the middle section of the Tanlu fault zone in the second phase is weaker than that in the first phase, and the maximum shear strain in the southern section is stronger than that in the first phase. The fault coupling coefficient of the south Sihong to Jiashan section is high, and it is also the unruptured section of historical earthquakes. At the same time, small earthquakes in this area are not active and accumulate stress easily, so the future earthquake risk deserves attention.

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.

Effects of Surface Etch Hole Fault on the Velocity Field in Microchannel Reactors

Microchannel reactors are commonly used in micro-chemical technology. The performance of microreactors is greatly affected by the velocity field in the microchannel. The flow field is disturbed by the cylindrical etch holes caused by air dust on the microchannel surface during its processing procedure. In this approach, a two-dimensional computational fluid dynamics （CFD） model is put forward to study the effect of etch holes on flow field. The influenced area of single or two concave etch holes is studied for the case of laminar flow. The hole diameter, the Reynolds number and the distance between the center of holes are found to have influences on the flow field. Numerical results indicate that the effects of etch hole on the flow field should be evaluated and the way of choosing the economic class of cleanroom for microreactor manufacture is presented.

A Computation on Velocity Field of the Post-flare Loop System of August 17, 1989

The observational line of sight velocity fields of a large post flare loop(PFL) system on August 17,1989 were derived based on “multi cloud model” (MCM) method from the Hβ spectral data obtained with the two dimensional Multi Waveband Spectra Spectro Helio Graph (MW SSHG) at Yunnan Observatory. In order to explain the fundamental features of the velocity fields observed,an assumption and approximation to be taken into account that the matter within a loop system rises up from the loop foot to the loop top along the loop leg in a helical magnetic line under the combined action of solar gravity, magnetic stress and a gradient of atmospheric pressure in the loop system, and a theoretical velocity field of the loop system was calcuated with the aid of MHD theory. From a comparison of both velocity fields, it is found that the computed velocity field is basically similar with the observed one. It indicates that the assumption and approximation mentioned above are suitable basically.

Three-dimensional transformation of magnetization direction and magnetic field component at low latitudes based on vertical relationship

The transformation of the magnetization direction and the magnetic fi eld component is one of the important methods in magnetic data processing and transformation,which can be conducted in both wavenumber and spatial domains.The transformation method in the wavenumber domain has simpler processing expression and higher processing effi ciency than in the spatial domain;however,they are unstable at low latitude.In this paper,the conclusion that the sum is 0 of two vertical magnetic fi eld components(magnetization inclinations are also perpendicular)in 2D is used for the 3D transformation of the magnetization direction and the magnetic field component.In addition,the transformation method at low latitudes based on vertical relationship(VMT)is proposed,which is an iterative algorithm that converts the transformation of the magnetization direction and the magnetic field component at the low latitude into the high latitude.This method restrains the instability of transformation of constant and variable magnetization direction and magnetic fi eld components in low latitudes.The accuracy,stability,and practicality are verifi ed from synthetic models and real data.

Impact of a background velocity field on solidification growth of single-crystal nuclei using the PF-LBM

To promote/inhibit ice formation in the natural environment and industrial systems,the growth and evolution process of ice single-crystal nuclei were simulated using the phase field-lattice Boltzmann method(PF-LBM),and the influence of a background flow field on the growth of single-crystal nucleus dendrites was also analyzed.The results show that the flow field makes dendrite growth asymmetric.The growth of dendrites is more developed on the upstream side than on the downstream side.The dendrite tip growth rate and tip radius are greater on the upstream side than on the downstream side.The solid phase ratio is greater with a background flow field than without one.The higher the flow velocity is,the more developed the dendrites on the upstream side,the faster the dendrites grow,and the higher the dendrite tip growth rate.The dendrites on the backflow side have a lower flow rate and a lower degree of supercooling than those on the upstream side,which inhibits the solidification process,the growth rate is slow,and the dendrites are underdeveloped.

Spatio-temporal variations of shallow seismic velocity changes in Salton Sea Geothermal Field,California in response to large regional earthquakes and long-term geothermal activities

We measure spatio-temporal variations of seismic velocity changes in Salton Sea Geothermal Field,California based on cross correlations of daily seismic traces recorded by a borehole seismic network from December 2007 to January 2014.We find clear co-seismic velocity reductions during the 2010 M 7.2 El Mayor–Cucapah,Mexico earthquake at~100 km further south,followed by long-term recoveries.The co-seismic reductions are larger with longer post-seismic recoveries in higher frequency bands,indicating that material damage and healing process mostly occurred in the shallow depth.In addition,the co-seismic velocity reductions are larger for ray paths outside the active fluid injection/extraction regions.The ray paths inside injection/extraction regions are associated with smaller co-seismic reductions,but subtle long-term velocity increases.We also build 3D transient water flow models based on monthly injection/extraction rates,and find correlations between several water flow parameters and co-seismic velocity reductions.We interpret the relative lack of co-seismic velocity changes within the geothermal region as unclogging of fracture network due to persistent fluid flows of geothermal production.The long-term velocity increase is likely associated with the ground water depletion and subsidence due to net production.

NUMERICAL SIMULATION AND ANALYSIS ON THE VELOCITY FIELD IN BATH OF A ROUND ELECTRIC CLEANING FURNACE

The flowing behavior of liquid slag has an important effect on the heat efficiency, recovery ratio ofvaluable metals and life span of the electric cleaning furnace. The velocity procedure using "κ-ε" model wasdeveloped under the cylindrical coordinate system. The procedure is used to calculate the velocity field ofmolten slag. The forces acting on molten slag were analyzed. The calculating method of electromagnetic forcewas described. The discrete equations, which were solved by using SIMPLE algorithm, were generated bymeans of interlace grids and control volume. The boundary layers near the solid wall were treated as wallfunction. The velocity distribution is obtained. The results show that the velocity of molten slag varies withdifferent horizontal sections. There are six eddies in the surface layer of molten slag. A large circle is formedin the sections below the electrodes, while the longitudinal section includes two large circles around the electrode. The influences of the features of the velocity field on the operation in the electric cleaning furnace arediscussed in detail.

Correction of Observation Interference on Z component of Geomagnetic Field From DC Power Lines

为减小直流输电线路对地磁场观测的干扰,采用数学计算的方法精确求解了线路对地磁Z分量的干扰水平。以观测点为原点建立了地磁仪器观测点坐标系,推导了观测点坐标系各坐标轴在地心坐标系中的单位矢量表达式。根据空间矢量点积的定义,得到了空间电流磁场矢量投影到地磁Z分量方向的求解公式。利用推导的公式和直流线路电流磁场三维计算模型,最终求得直流线路对地磁Z分量观测的干扰水平。由此,开发了＂直流线路对地磁观测干扰校正系统＂,并开展了地磁台站和野外联合地磁观测干扰实验。结果表明,直流线路和地磁观测台站距离不超过10km时,提出的数学模型和计算方法误差不超过5%。

Principal-component estimates of the Kuroshio Current axis and path based on the mathematical verification between satellite altimeter and drifting buoy data

We used satellite altimetry data to investigate the Kuroshio Current because of the higher resolution and wider range of observations. In previous studies, satellite absolute geostrophic velocities were used to study the spatiotemporal variability of the sea surface velocity field along the current, and extraction methods were employed to detect the Kuroshio axes and paths. However, sea surface absolute geostrophic velocity estimated from absolute dynamic topography should be regarded as the geostrophic component of the actual surface velocity, which cannot represent a sea surface current accurately. In this study, mathematical verification between the climatic absolute geostrophic and bin-averaged drifting buoy velocity was established and then adopted to correct the satellite absolute geostrophic velocities. There were some differences in the characteristics between satellite geostrophic and drifting buoy velocities. As a result, the corrected satellite absolute geostrophic velocities were used to detect the Kuroshio axis and path based on a principal-component detection scheme. The results showed that the detection of the Kuroshio axes and paths from corrected absolute geostrophic velocities performed better than those from satellite absolute geostrophic velocities and surface current estimations. The corrected satellite absolute geostrophic velocity may therefore contribute to more precise day-to-day detection of the Kuroshio Current axis and path.

Study on the Lightning-radaited Electric Field Components at Mesospheric Altitudes

In this paper we have extended the semi-analytical model presented by Yashunin et al.to compute both of the lightning vertical and horizontal electric field considering the effect of nonlinear atmospheric electric parameters,and validate its accuracy by using our 2-D FDTD method.It is noted that the semi-analytical model can predict not only the lightning vertical field,but also the horizontal electric field with a satisfied accuracy in a small region with a radius ranging from 5 to 30 km at a height from 80 to 90 km above the lightning strike point.Compared with the results in the free space with a vacuum,the conducting altitude atmosphere has a relatively little effect on the radiated field component,but it has much effect on the electrostatic and induction component because their duration times are longer than the local relaxation time.For example,the induction field component in the conducting nonlinear medium decreases by about 43%,compared with that in the free space with a vacuum.The induction field is larger than the radiation field component within a radius of 25 km at a height of 80 km above the lightning strike point.Therefore,although the induction component may be likely one extra energy source(except for the radiation field)for creating the pancake shape of some elves with a luminous central part,we have to take into account the effect of nonlinear atmospheric electric parameters on the prediction of induced field component peak value that is needed for the production of elves.

Determination of the Velocity of a Magnetic-Field Driven Gliding Arc by Using a Photo-Multiplier

The gliding arc is an important approach to production of non-thermal plasma at atmospheric pressure, it can offer high-energy efficiency and high-electivity for chemical reactions. In this paper, the gliding arc driven by the transverse magnetic field is described and its velocity is measured by using a photo-multiplier. The mean velocity of the gliding arc increases with increasing magnetic induced-intensity, and its value varies from 7.8 m/s to 32 m/s.

A study on the discrete image method for calculation of transient electromagnetic fields in geological media

We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized： the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the ＂ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.

Herrmann Method of Analyzing Structure Design Velocity Field

The shape optimization is studied by adopting the domain integrated method which is based on the calculus of variations during the shape design sensitivity analysis. A new method of improving the efficiency of the design velocity field analysis and the quality of the finite element method (FEM) mesh is put forward. The sensitivity analysis which is based on the calculus of variations is used in the shape optimization. The design velocity field is solved by Herrmann method. An example shows that both the quality of the FEM mesh and the efficiency of the computing of the design velocity field are improved by Herrmann method. So the effect and the efficiency of the shape optimization are guaranteed. If using sensitivity analysis which is based on the calculus of variations in the shape optimization, the sensitivity analysis can be a relatively independent module. The efficiency of computing the design velocity field and the quality of mesh will be improved by using Herrmann method.

3D FEM simulation of flow velocity field for 5052 aluminum alloy multi-row sprocket in cold semi-precision forging process

Based on the design of the multi-row sprocket with a new tooth profile,a cold semi-precision forging process for manufacturing 5052 aluminum alloy multi-row sprocket was presented.Through simulating the forging process of 5052 aluminum alloy sprocket billet with 3D rigid-viscoplastic FEM,both the distributions of flow velocity field in axial(U_Z),radial(U_R) and circumferential(U_θ) directions and the curves of velocity component in different deformation regions were respectively obtained.By comparison and analysis of the velocity varying curves,the velocity component relation conditions for filling the die cavity were clarified.It shows that when the die cavity is almost fully filled,the circumferential velocity U_θ increases sharply,implying that U_θplays a key role in fully filling the die cavity.

ST-LSTM-SA:A New Ocean Sound Velocity Field Prediction Model Based on Deep Learning

The scarcity of in-situ ocean observations poses a challenge for real-time information acquisition in the ocean.Among the crucial hydroacoustic environmental parameters,ocean sound velocity exhibits significant spatial and temporal variability and it is highly relevant to oceanic research.In this study,we propose a new data-driven approach,leveraging deep learning techniques,for the prediction of sound velocity fields(SVFs).Our novel spatiotemporal prediction model,STLSTM-SA,combines Spatiotemporal Long Short-Term Memory(ST-LSTM) with a self-attention mechanism to enable accurate and real-time prediction of SVFs.To circumvent the limited amount of observational data,we employ transfer learning by first training the model using reanalysis datasets,followed by fine-tuning it using in-situ analysis data to obtain the final prediction model.By utilizing the historical 12-month SVFs as input,our model predicts the SVFs for the subsequent three months.We compare the performance of five models:Artificial Neural Networks(ANN),Long ShortTerm Memory(LSTM),Convolutional LSTM(ConvLSTM),ST-LSTM,and our proposed ST-LSTM-SA model in a test experiment spanning 2019 to 2022.Our results demonstrate that the ST-LSTM-SA model significantly improves the prediction accuracy and stability of sound velocity in both temporal and spatial dimensions.The ST-LSTM-SA model not only accurately predicts the ocean sound velocity field(SVF),but also provides valuable insights for spatiotemporal prediction of other oceanic environmental variables.