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.

Assessment of liquefaction potential based on shear wave velocity:Strain energy approach

Liquefaction assessment based on strain energy is significantly superior to conventional stress-based methods.The main purpose of the present study is to investigate the correlation between shear wave velocity and strain energy capacity of silty sands.The dissipated energy until liquefaction occurs was calculated by analyzing the results of three series of comprehensive cyclic direct simple shear and triaxial tests on Ottawa F65,Nevada,and Firoozkuh sands with varying silt content by weight and relative densities.Additionally,the shear wave velocity of each series was obtained using bender element or resonant column tests.Consequently,for the first time,a liquefaction triggering criterion,relating to effective overburden normalized liquefaction capacity energy(WL=s’c)to effective overburden stresscorrected shear wave velocity(eVs1)has been introduced.The accuracy of the proposed criteria was evaluated using in situ data.The results confirm the ability of shear wave velocity as a distinguishing parameter for separating liquefied and non-liquefied soils when it is calculated against liquefaction capacity energy(WL=s’c).However,the proposed WL=s’c-Vs1 curve,similar to previously proposed cyclic resistance ratio(CRR)-Vs1 relationships,should be used conservatively for fields vulnerable to liquefaction-induced lateral spreading.

Shear wave velocity prediction:A review of recent progress and future opportunities

Shear logs,also known as shear velocity logs,are used for various types of seismic analysis,such as determining the relationship between amplitude variation with offset(AVO)and interpreting multiple types of seismic data.This log is an important tool for analyzing the properties of rocks and interpreting seismic data to identify potential areas of oil and gas reserves.However,these logs are often not collected due to cost constraints or poor borehole conditions possibly leading to poor data quality,though there are various approaches in practice for estimating shear wave velocity.In this study,a detailed review of the recent advances in the various techniques used to measure shear wave(S-wave)velocity is carried out.These techniques include direct and indirect measurement,determination of empirical relationships between S-wave velocity and other parameters,machine learning,and rock physics models.Therefore,this study creates a collection of employed techniques,enhancing the existing knowledge of this significant topic and offering a progressive approach for practical implementation in the field.

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.

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.

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.

High-resolution Rayleigh wave phase velocity maps from ambient noise tomography in North China

We presented high-resolution Rayleigh wave phase velocity maps at periods ranging from 5 s to 30 s in the northeast part of the North China Craton (NNCC). Continuous time-series of vertical component between October 2006 and December 2008, recorded by 187 broadband stations temporarily deployed in the NNCC region, have been cross-correlated to obtain estimated fundamental mode Rayleigh wave Green’s functions. Using the frequency and time analysis technique based on continuous wavelet transformation, we measured 3 667 Rayleigh wave phase velocity dispersion curves. High-resolution phase velocity maps at periods of 5, 10, 20 and 30 s were reconstructed with grid size 0.25°× 0.25°, which reveal lateral heterogeneity of shear wave structure in the crust and upper mantle of NNCC. For periods shorter than 10 s, the phase velocity variations are well correlated with the principal geological units in the NNCC, with low-speed anomalies corresponding to the major sedimentary basins and high-speed anomalies coinciding with the main mountain ranges. Within the period range from 20 s to 30 s, high phase velocity observed in eastern NCC is coincident with the thin crust, whereas low phase velocities imaged in central NCC is correlated to the thick crust. However, the low-velocity anomaly in the Beijing-Tianjin-Tangshan region displayed in the 20 s and 30 s phase maps may be associated with fluids.

Tunnel effect of fractal fault and transient S-wave velocity rupture (TSVR) of in-plane shear fault

Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip) fault has been proved, but in 2 dimensional classical model, there are two difficulties in transient S wave velocity rupture, i.e ., initialization difficulty and divergence difficulty in interpreting the realization of TSVR. The initialization difficulty means, when v ↑ v R (Rayleigh wave velocity), the dynamic stress strength factor K 2(t) →+0, and changes from positive into negative in the interval ( v R, β ). How v transit the forbidden of ( v R, β )? The divergence difficulty means K 2(t) →+∞ when v ↓ β . Here we introduce the concept of fractal and tunnel effect that exist everywhere in fault. The structure of all the faults is fractal with multiple cracks. The velocity of fault rupture is differentiate of the length of the fault respect to time, so the rupture velocity is also fractal. The tunnel effect means the dynamic rupture crosses over the interval of the cracks, and the coalescence of the intervals is slower than the propagation of disturbance. Suppose the area of earthquake nucleation is critical or sub critical propagation everywhere, the arriving of disturbance triggers or accelerates the propagation of cracks tip at once, and the observation system cannot distinguish the front of disturbance and the tip of fracture. Then the speed of disturbance may be identified as fracture velocity, and the phenomenon of TSVR appears, which is an apparent velocity. The real reason of apparent velocity is that the mathematics model of shear rupture is simplified of complex process originally. The dual character of rupture velocity means that the apparent velocity of fault and the real velocity of micro crack extending, which are different in physics, but are unified in rupture criterion. Introducing the above mentioned concept to the calculation of K 2 (t) , the difficulty of initialization can be overcome, and the integral equation of triggering the initialization of TSVR is given quantitatively. By solving this integral equation, the lower limit of TSVR is 1.105 3 β , not β , and the divergence difficulty is overcome. TSVR is unstable solution, and may degenerate to sub Rayleigh wave velocity rupture immediately where the non critical condition can be measured. The results of this paper show that the initialization and continuum depends on the condition of earthquake nucleation in seismogenic area.