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.

Shear wave velocity model using HVSR inversion beneath Bandar Lampung City

The horizontal-to-vertical spectral ratio(HVSR)method has been used to characterize site-effect parameters that are indispensable in seismic hazard and risk-reduction studies in urban areas and rapid land-use planning.This method is widely used because it is the cheapest and simplest geophysical method for the acquisition and processing stages.In subsequent developments,the HVSR method has been widely used to determine elastic rock parameters,particularly shear wave velocity(v_(S)),through the HVSR curve inversion process.Furthermore,the v_(S)structural model can be used to delineate the presence of complex geological structures,particularly faults and sedimentary basins.Bandar Lampung is a city in Lampung Province with many fault structures and groundwater basins to the south.There are 83 HVSR measurement points around Bandar Lampung for delineating the presence of fault structures and groundwater basins.We produced the HVSR curve from the measurement results and then performed an inversion process using the particle swarm optimization algorithm to obtain v_(S)for the depth profile.Subsequently,from this profile,we produced a two-dimensional(2D)lateral and vertical model.The mean v_(S)value was calculated from all the measurement points,and we found stiff soil layers reaching depths of approximately 5 m,with a value of v_(S)<330 m/s.A bedrock layer with a velocity exceeding 1250 m/s was visible at a depth of 100 m.Based on the 2D model,the v_(S)structure shows that the city of Bandar Lampung is divided into two zones,with a NW-SE boundary.The north-middle-eastern part of the city consists of harder rocks.This harder rock is characterized by extremely high v_(S)values,starting from a depth of 50 m.In contrast,the south-middle-west exhibits a low-moderate v_(S)anomaly associated with groundwater basins SW of the city.From the 2D v_(S)structural model,fault structures can be found along the city,characterized by a contrast of v_(S)values from low to medium and from medium to high.

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.

Analysis of the Impact of the Blazhko Effect Both on the Van Hoof Effect and Radial Velocity Amplitude in the Star RR Lyr

The Van Hoof effect is a phase shift existing between the radial velocity curves of hydrogen and metallic lines within the atmosphere of pulsating stars.In this article,we present a study of this phenomenon through the spectra of the brightest pulsating star RR Lyr of RR Lyrae stars recorded for 22 yr.We based ourselves,on the one hand,on 1268 spectra(41 nights of observation)recorded between the years 1994 and 1997 at the Observatory of Haute Provence(OHP,France)previously observed by Chadid and Gillet,and on the other hand on 1569 spectra(46nights of observation)recorded at our Oukaimeden Observatory(Morocco)between 2015 and 2016.Through this study,we have detected information on atmospheric dynamics that had not previously been detected.Indeed,the Van Hoof effect which results in a clear correlation between the radial velocities of hydrogen and those of the metallic lines has been observed and analyzed at different Blazhko phases.A correlation between the radial velocities of different metallic lines located in the lower atmosphere has been observed as well.For the first time,we were able to show that the amplitude of the radial velocity curves deduced from the lines of hydrogen and that of FeⅡ(λ4923.921?)increases toward the minimum of the Blazhko cycle and decreases toward the maximum of the same Blazhko cycle.Furthermore,we found that the Van Hoof effect is also modulated by the Blazhko effect.Thus,toward the minimum of the Blazhko cycle the Van Hoof effect is more visible and at the maximum of the Blazhko cycle,this effect is minimal.We also observed the temporal evolution of the amplitudes of the radial velocities of the lower and upper atmosphere.When observed over a long time,we can interpret it as a function of the Blazhko phases.

Where Is Phase Velocity in Minkowski Space?

In the special theory of relativity, massive particles can travel at neither the speed of light c nor faster. Meanwhile, since the photon was quantized, many have thought of it as a point particle. How pointed? The idea could be a mathematical device or physical simplification. By contrast, the preceding notion of wave-group duality has two velocities: a group velocity vg and a phase velocity vp. In light vp = vg = c;but it follows from special relativity that, in massive particles, vp > c. The phase velocity is the product of the two best measured variables, and so their product constitutes internal motion that travels, verifiably, faster than light. How does vp then appear in Minkowski space? For light, the spatio-temporal Lorentz invariant metric is s2=c2t2−x2−y2−z2, the same in whatever frame it is viewed. The space is divided into 3 parts: firstly a cone, symmetric about the vertical axis ct > 0 that represents the world line of a stationary particle while the conical surface at s = 0 represents the locus for light rays that travel at the speed of light c. Since no real thing travels faster than the speed of light c, the surface is also a horizon for what can be seen by an observer starting from the origin at time t = 0. Secondly, an inverted cone represents, equivalently, time past. Thirdly, outside the cones, inaccessible space. The phase velocity vp, group velocity vg and speed of light are all equal in free space, vp = vg = c, constant. By contrast, for particles, where causality is due to particle interactions having rest mass mo > 0, we have to employ the Klein-Gordon equation with s2=c2t2−x2−y2−z2+mo2c2. Now special relativity requires a complication: vp.vg = c2 where vg c and therefore vp > c. In the volume outside the cones, causality due to light interactions cannot extend beyond the cones. However, since vp > c and even vp >> c when wavelength λ is long, extreme phase velocities are then limited in their causal effects by the particle uncertainty σ, i.e. to vgt ± σ/ω, where ω is the particle angular frequency. This is the first time the phase range has been described for a massive particle.

Intrinsic Wave Velocity Propagation:A Novel Parameter for Assessing the Effect of Anthracycline Chemotherapy Agents on Cardiac Diastolic Function in Breast Cancer Patients

Objective Anthracycline chemotherapeutic agents have significant cardiotoxicity.The present study emphasized the effect of anthracycline chemotherapy drugs on left ventricular(LV)myocardial stiffness in breast cancer patients by measuring the intrinsic wave velocity propagation(IVP),and evaluating the potential clinical value of IVP in detecting early LV diastolic function impairment.Methods A total of 68 newly diagnosed breast cancer patients,who were treated with anthracycline-based chemotherapy,were analyzed.Transthoracic echocardiography was performed at baseline(T0),and after 1,2,3,4 and 8 chemotherapeutic cycles(T1,T2,T3,T4 and T5,respectively).Then,the IVP,LV strain parameters[global longitudinal strain(GLS),longitudinal peak strain rate at systole(LSRs),longitudinal peak strain rate at early diastole(LSRe),longitudinal peak strain rate at late diastole(LSRa),and the E/LSRe ratio],and conventional echocardiographic parameters were obtained and further analyzed.A relative reduction of>15%in GLS was considered a marker of early LV subclinical dysfunction.Results Compared to the T0 stage,IVP significantly increased at the T1 stage.However,there were no significant changes in GLS,LSRs,or LSRe between the T0 and T1 stages.These parameters significantly decreased from the T2 stage.LSRa started to significantly decrease at the T5 stage,and the E/LSRe ratio started to significantly increase at the T3 stage(all P<0.05).At the T0 stage,IVP(AUC=0.752,P<0.001)had a good predictive value for LV subclinical dysfunction after chemotherapy.Conclusions IVP is a potentially sensitive parameter for the early clinical assessment of anthracycline-related cardiac diastolic impairment.

How Classical, Quantum-Mechanical, and Relativistic Wave and Field Equations Are Uniformly Generated by Velocity-Field Divergence Equations

We expand previously established results concerning the uniform representability of classical and relativistic gravitational field equations by means of velocity-field divergence equations by demonstrating that conservation equations for (probability) density functions give rise to velocity-field divergence equations the solutions of which generate—by way of superposition—the totality of solutions of various well-known classical and quantum-mechanical wave equations.

Converted wave AVO inversion for average velocity ratio and shear wave reflection coefficient

Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the relations between the extrema position and amplitude, average velocity ratio across the interface, and shear wave reflection coefficient are derived. The extrema position is a monotonically decreasing function of average velocity ratio, and the extrema amplitude is a function of average velocity ratio and shear wave reflection coefficient. For theoretical models, the average velocity ratio and shear wave reflection coefficient are inverted from the extrema position and amplitude obtained from fitting a power function to converted wave AVO curves. Shear wave reflection coefficient sections have clearer physical meaning than conventional converted wave stacked sections and establish the theoretical foundation for geological structural interpretation and event correlation. ＂The method of inverting average velocity ratio and shear wave reflection coefficient from the extrema position and amplitude obtained from fitting a power function is applied to real CCP gathers. The inverted average velocity ratios are consistent with those computed from compressional and shear wave well logs.

Seismic tomography of Yunnan region using short-period surface wave phase velocity

利用云南数字地震台网23个台站记录的短周期(1～18s)面波资料，用双台窄带通滤波互相关方法，测定了在209条路径上的基阶瑞利波相速度频散曲线，并通过层析成像反演得到各周期的相速度分布图象．给出了沿24°N，25°N，26°N及27°N和100．5°E的相速度频散剖面图．结果表明：云南地区的相速度分布图象具有强烈的横向变化，短周期段的相速度分布与沉积层厚度密切相关，中甸、丽江东部区域在约25km深度存在着明显的低速层，而畹町、云县、思茅一带在该深度范围是一明显的大面积高速区域．小江断裂与红河断裂围成的川滇菱形块体南段为低相速度区域．云南地区强震震中位置主要分布在相速度高速与低速过渡的区域．

Three-dimensional P-wave velocity structure of the crust beneath Hainan Island and its adjacent regions,China

Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999～2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions has been determined. The results show that the pattern of velocity anomalies in the shallower upper crust is somewhat associated with the surface geological tectonics in the region. A relative low-velocity anomaly appears north of the Wangwu-Wenjiao fault zone and a relative high-velocity anomaly appears south of the Wangwu-Wenjiao fault zone,corresponding to the depressed areas in north Hainan Island,where many volcanoes are frequently active and geothermal values are relatively higher,and the uplifted and stable regions in central and south of the Hainan Is-land. In the middle and lower crust velocities are relatively lower in east Hainan than those in west Hainan,possi-bly suggesting the existence of the upwelling of hot materials from the mantle in east Hainan. The pattern of veloc-ity anomalies also indicates that NW faults,i.e.,the Puqian-Qinglan fault,may be shallower,while the E-W Wangwu-Wenjiao fault may be deeper,which perhaps extends down to Moho depth or deeper.

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in cen-tral-western China

采用中国中西部地区(2l°～36°N，98°～112°E)193个地震台在1992～1999年间记录到的9988次地震的Pg和Sg震相走时的读数资料，用Roecker的SPHYPIT90程序，反演了该地区三维地壳P波速度结构，并用SPHREL3D90程序进行了地震的重新定位．反演结果揭示了中国中西部地区地震P波速度结构明显的横向不均匀性，这些不同深度上波速的横向变化多以该地区的活动断裂为分界线．可以看出活动断裂两侧存在明显的速度反差．通过重新定位，得到了6459次地震的震源参数，这些精确定位的地震震中明显沿该区活动断裂呈现条带状分布，其范围和尺度清晰地表示了这一地区地震活动与活动断裂的紧密关系．其中，82％重新精确定位的事件的震源深度在20km以内．这一结果与笔者用双差地震定位法得到的重新定位的震源深度分布相一致．

Three dimensional shear wave velocity structure of crust and upper mantle in South China Sea and its adjacent regions by surface waveform inversion

We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). It partitions the large-scale optimization problem into a number of independent small-scale problems. We adopted surface waveform inversion with an equal block (2((2() discretization in order to acquire the images of shear velocity structure at different depths (from surface to 430 km) in the crust and upper-mantle. The resolution of all these anomalies has been established with (check-board( resolution tests. These results show significant difference in velocity, lithosphere and asthenosphere structure between South China Sea and its adjacent regions.

Digital core based transmitted ultrasonic wave simulation and velocity accuracy analysis

Transmitted ultrasonic wave simulation （TUWS） in a digital core is one of the important elements of digital rock physics and is used to study wave propagation in porous cores and calculate equivalent velocity. When simulating wave propagates in a 3D digital core, two additional layers are attached to its two surfaces vertical to the wave-direction and one planar wave source and two receiver-arrays are properly installed. After source excitation, the two receivers then record incident and transmitted waves of the digital rock. Wave propagating velocity, which is the velocity of the digital core, is computed by the picked peak-time difference between the two recorded waves. To evaluate the accuracy of TUWS, a digital core is fully saturated with gas, oil, and water to calculate the corresponding velocities. The velocities increase with decreasing wave frequencies in the simulation frequency band, and this is considered to be the result of scattering. When the pore fluids are varied from gas to oil and finally to water, the velocity-variation characteristics between the different frequencies are similar, thereby approximately following the variation law of velocities obtained from linear elastic statics simulation （LESS）, although their absolute values are different. However, LESS has been widely used. The results of this paper show that the transmission ultrasonic simulation has high relative precision.

Shear wave velocity is a useful marker for managing nonalcoholic steatohepatitis

AIM:To investigate whether a noninvasive measurement of tissue strain has a potential usefulness for management of nonalcoholic steatohepatitis(NASH).METHODS:In total 26 patients,23 NASHs and 3 normal controls were enrolled in this study.NASH was staged based on Brunt criterion.At a region of interest(ROI),a shear wave was evoked by implementing an acoustic radiation force impulse(ARFI),and the propagation velocity was quantif ied.RESULTS:Shear wave velocity(SWV) could be reproducibly quantified at all ROIs in all subjects except for 4 NASH cases,in which a reliable SWV value was not calculated at several ROIs.An average SWV of 1.34 ± 0.26 m/s in fibrous stage 0-1 was significantly slower than 2.20 ± 0.74 m/s and 2.90 ± 1.01 m/s in stages 3 and 4,respectively,but was not significantly different from 1.79 ± 0.78 m/s in stage 2.When a cutoff value was set at 1.47 m/s,receiver operating characteristic analysis showed significance to dissociate stages 3 and 4 from stage 0-1(P=0.0092) with sensitivity,specificity and area under curve of 100%,75% and 94.2%,respectively.In addition,the correlation between SWV and hyaluronic acid was significant(P<0.0001),while a tendency toward negative correlation was observed with serum albumin(P=0.053).CONCLUSION:The clinical implementation of ARFI provides noninvasive repeated evaluations of liver stiffness at an arbitrary position,which has the potential to shed new light on NASH management.

Shear wave velocity-based liquefaction evaluation in the great Wenchuan earthquake: a preliminary case study

The great Wenchuan earthquake （Ms= 8.0） in 2008 caused severe damage in the western part of the Chengdu Plain. Soil liquefaction was one of the major causes of damage in the plain areas, and proper evaluation of liquefaction potential is important in the definition of the seismic hazard facing a given region and post-earthquake reconstruction. In this paper, a simplified procedure is proposed for liquefaction assessment of sandy deposits using shear wave velocity （V）, and soil liquefaction from the Banqiao School site was preliminarily investigated after the earthquake. Boreholes were made at the site and shear wave velocities were measured both by SASW and down-hole methods. Based on the in-situ soil information and V profiles, the liquefaction potential of this site was evaluated. The results are reasonably consistent with the actual field behavior observed after the earthquake, indicating that the proposed procedure is effective. The possible effects of gravel and fines contents on liquefaction of sandy soils were also briefly discussed.

Seismic tomography of Yunnan region using short-period surface wave phase velocity

The data of short-period (1~18 s) surface waves recorded by 23 stations belonging to the digital seismic network of Yunnan Province of China are used in this paper. From these data, the dispersion curves of phase velocities of the fundamental mode Rayleigh wave along 209 paths are determined by using the two-station narrowband filtering cross-correlation method. Adopting tomography method, the distribution maps of phase velocities at various peri-ods in Yunnan region are inverted. The maps of phase velocities on profiles along 24N, 25N, 26N, 27N and 100.5E and the distribution maps of phase velocities at 3 periods in the study region are given. The results show that the phase velocity distribution in Yunnan region has strong variations in horizontal direction, and the phase velocity distribution in short-period range is closely related to the thickness of sedimentary layers in the shallow crust. The phase velocity in southern part of the Sichuan-Yunnan rhombic block encircled by the Honghe fault and Xiaojiang fault is obviously lower than that in surrounding areas. The epicentral locations of strong earthquakes in Yunnan region are mainly distributed in transitional zones between low and high phase velocities.

Wave-Dependence of Friction Velocity, Roughness Length, and Drag Coefficient over Coastal and Open Water Surfaces by Using Three Databases

The parameterization of friction velocity, roughness length, and the drag coefficient over coastal zones and open water surfaces enables us to better understand the physical processes of air-water interaction. In context of measurements from the Humidity Exchange over the Sea Main Experiment （HEXMAX）, we recently proposed wave-parameter dependent approaches to sea surface friction velocity and the aerodynamic roughness by using the dimensional analysis method. To extend the application of these approaches to a range of natural surface conditions, the present study is to assess this approach by using both coastal shallow （RASEX） and open water surface measurements （Lake Ontario and Grand Banks ERS-1 SAR） where wind speeds were greater than 6.44 m s-1. Friction velocities, the surface aerodynamic roughness, and the neutral drag coefficient estimated by these approaches under moderate wind conditions were compared with the measurements mentioned above. Results showed that the coefficients in these approaches for coastal shallow water surface differ from those for open water surfaces, and that the aerodynamic roughness length in terms of wave age or significant wave height should be treated differently for coastal shallow and open water surfaces.

Fracturing process and effect of fracturing degree onwave velocity of a crystalline rock

The present paper investigates the effect of fracturing degree on P-and S-wave velocities in rock.The deformation of intact brittle rocks under loading conditions is characterized by a microcracking procedure,which occurs due to flaws in their microscopic structure and propagates through the intact rock,leading to shear fracture.This fracturing process is of fundamental significance as it affects the mechanical properties of the rock and hence the wave velocities.In order to determine the fracture mechanism and the effect of fracturing degree,samples were loaded at certain percentages of peak strength and ultrasonic wave velocity was recorded after every test.The fracturing degree was recorded on the outer surface of the sample and quantified by the use of the indices P;（traces of joints/m）,P;（traces of joints/m;） and P;（length of fractures/m;）.It was concluded that the wave velocity decreases exponentially with increasing fracturing degree.Additionally,the fracturing degree is described adequately with the proposed indices.Finally,other parameters concerning the fracture characteristics,rock type and scale influence were found to contribute to the velocity decay and need to be investigated further.

Shear wave velocity-based evaluation and design of stone column improved ground for liquefaction mitigation

The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance （CRR）-shear wave velocity （V）-void ratio （e） of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V of soil is raised high enough （i.e., no less than the critical velocity） to resist the given earthquake loading according to the CRR-V relationship, and then this requirement is transferred to the control of target void ratio （i.e., the critical e） according to the V-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.

Group velocity distribution of Rayleigh waves and crustal and upper mantle velocity structure of the Chinese mainland and its vicinity

Based on the long period digital surface wave data recorded by 11 CDSN stations and 11 IRIS stations, the dispersion curves of the group velocities of fundamental mode Rayleigh waves along 647 paths, with the periods from 10 s to 92 s, were measured by multi-filter. Their distribution at 25 central periods within the region of 18~54N, 70~140E was inverted by Dimtar-Yanovskaya method. Within the period from 10 s to 15.9 s, the group velocity distribution is laterally inhomogeneous and is closely related to geotectonic units, with two low velocity zones located in the Tarim basin and the East China Sea and its north regions, respectively. From 21 s to 33 s, the framework of tectonic blocks is revealed. From 36.6 s to 40 s, the lithospheric subdivision of the Chinese mainland is obviously uncovered, with distinct boundaries among the South-North seismic belt, the Tibetan plateau, the North China, the South China and the Northeast China. Four cross-sections of group velocity distribution with period along 30N, 38N, 90E and 120E, are discussed, respectively, which display the basic features of the crust and upper mantle of the Chinese mainland and its neighboring regions. There are distinguished velocity differences among the different tectonic blocks. There are low-velocity-zones (LVZ) in the middle crust of the eastern Tibetan plateau, high velocity featured as stable platform in the Tarim basin and the Yangtze platform, shallow and thick low-velocity-zone in the upper mantle of the North China. The upper mantle LVZ in the East China Sea and the Japan Sea is related to the frictional heat from the subduction of the Philippine slab and the strong extension since the Himalayan orogenic period.