Joint inversion of Rayleigh group and phase velocities for S-wave velocity structure of the 2021 M_(S)6.0 Luxian earthquake source area,China

On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.

A trigger-tube tracer dilution technique for determining Darcy and apparent velocities of groundwater in dug wells: A case study on phreatic aquiferous formation in Bamenda-Cameroon

The need to understand flow within aquiferous formations for a complete evaluation of groundwater resource and quality control prompts the determination of groundwater velocities through well dilution technique. Well dilution techniques utilize tracer solutions after establishing an initial homogenous condition to monitor the flow rate of ambient groundwater into the wells. Application of dilution techniques in wells makes it feasible to determine the velocities of groundwater in the aquiferous formation surrounding the well. In this study, a simple trigger-tube tracer dilution technique was employed to determine the Darcy and apparent velocities of groundwater in the phreatic aquiferous formation in Bamenda, Cameroon. Eighteen (18) hand dug-wells at different locations within Bamenda were sampled by utilizing sodium chloride (NaCl) as the conservative tracer. Field estimates of groundwater flow velocities in the phreatic aquiferous formation in Bamenda reveal Darcy's groundwater velocity in the range of 0.39 m/d at Nacho to 130.64 m/d at Foncha Street and apparent velocity in the range of 0.78 m/d at Nacho to 277.86 m/d at Foncha Street. The immense variations in the velocities of groundwater indicate that the groundwater flows at different rates and directions within the aquiferous formation in Bamenda, possibly due to variations in their hydraulic conductivities. Moreover, the spatial variations in the formation types, facies changes, thickness, and layering of the aquiferous formation also contribute to the variation of velocities. Areas with low groundwater velocities are associated with a lower contaminant transport rate when compared to areas with high groundwater velocities. The findings of this study are important for assessing the rates of pollutant movement in the subsurface, as well as the effectiveness and efficacy of the trigger-tube technique in evaluating the hydraulic properties of aquiferous formations.

3D S-wave velocity structure of the Ningdu basin in Jiangxi province inferred from ambient noise tomography with dense array

The Ningdu basin,located in southern Jiangxi province of southwest China,is one of the Mesozoic basin groups which has exploration prospects for geothermal energy.A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration.In this study,we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution.Based on the dense seismic array including 35 short-period(5 s-100 Hz)seismometers with an average interstation distance of~5 km,Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion.Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm.The results revealed obvious low-velocity anomalies in the center of the basin,consistent with the low-velocity Cretaceous sedimentary rocks.The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies.The obvious seismic interface is about 2 km depth in the basin center and decreases to 700 m depth near the basin boundary,suggesting spatial thickness variations of the Cretaceous sediment.The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basincontrolling fault,which may provide possible upwelling channels for geothermal fluid.This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.

Apparent velocity estimation with P-SV ratio method

The apparent velocity of the incident wave is an important parameter for simulating rotational ground motion with theoretical methods, but it is difficult to estimate effectively when there is only a single record. This paper discusses a P-SV ratio method based on elastodynamic theory in a multi-layer isotropic elastic half space. The apparent velocities of four earthquakes in the SMART1 array are calculated with this method. The result is close to a method that uses travel time analysis. Furthermore, the factors that impact the apparent velocity and equivalent incident angle are considered according to records from the Chi-Chi earthquake. There is no obvious relationship between the equivalent incident angle and epicenter distance. However, the equivalent incident angle is obviously dependent on the site conditions.

Numerical Simulation of the Response Features of Apparent Seismic Wave Velocity Ratio in a Horizontal Layered Medium

This paper deals with the response features of AR(apparent ratio of seismic wave velocities to the changes of TR(true ratio of wave velocities)in the horizontal layered model by mathematical modeling.The results show that:(1)the response features of AR are associated with the parameters of the structure and its dynamic changes,and the relative position between the hypocenters and the monitoring networks,showing complicated patterns strongly related to the concrete paths of propagation of seismic waves from the source to the receiver in the observatories of the network;(2)the depth of the seismic source would have important influence on the response features of AR,especially the capacity to carry the anomalous information in the condition of the earth media,being in the anomalous state would be greater for those earthquakes which occur inside the anomalous layers than those underneath the anomalous layers;(3)the response features of AR are clearly related to the changes of TR(true ratio of wave velocities)instead of changes of wave velocities themselves,i.e.the response could be small as the changes in TR is small even in the case of large changes in the wave velocities.It is suggested that more attention must be paid to all these features in combination with detailed investigation of the velocity structure of the earth media in the study region and best fitting of precise hypocenter locations when one wants to obtain the reliable precursors from the changes in AR.

Rock critical porosity inversion and S-wave velocity prediction

A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical porosity value and we can generally take only an empirical critical porosity value which often causes errors. In this paper, we propose a method to obtain the rock critical porosity value by inverting P-wave velocity and applying it to predict S-wave velocity. The applications of experiment and log data both show that the critical porosity inversion method can reduce the uncertainty resulting from using an empirical value in the past and provide the accurate critical porosity value for predicting S-wave velocity which significantly improves the prediction accuracy.

The bound weighted average method(BWAM)for predicting S-wave velocity

The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.

Estimation of Shallow S-Wave Velocity Structure of Two Practical Sites from Microtremors Array Observation in Tangshan Area

Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase velocities of Rayleigh wave are estimated from vertical component records and those of Love wave are estimated from three-component records of microtremors array using modified spatial auto-correlation method.Haskell matrix method is used in calculating Rayleigh and Love wave phase velocities,and the shallow S-wave velocity structure of two practical sites are estimated by means of a hybrid approach of Genetic Algorithm and Simplex.The results are compared with the PS logging data of the two sites,showing it is feasible to estimate the shallow S-wave velocity structure of practical site from the observation of microtremor array.

Estimation of S-wave Velocity for Gas Hydrate Reservoir in the Shenhu Area,North South China Sea

Estimation of S-wave velocity using logging data has mainly been performed for sandstone, mudstone and oil and gas strata, while its application to hydrate reservoirs has been largely overlooked. In this paper we present petxophysical methods to estimate the S-wave velocity of hydrate reservoirs with the P-wave velocity and the density as constraints. The three models used in this paper are an equivalent model （MBGL）, a three-phase model （TPBE）, and a thermo-elasticity model （TEM）. The MBGL model can effectively describe the internal relationship among the components of the rock, and the estimated P-wave velocities are in good agreement with the measured data （2.8% error）. However, in the TPBE model, the solid, liquid and gas phases axe considered to be independent of each other, and the estimation results are relatively low （46.6% error）. The TEM model is based on the sensitivity of the gas hydrate to temperature and pressure, and the accuracy of the estimation results is also high （3.6% error）. Before the estimation, the occurrence patterns of hydrates in the Shenhu area were examined, and occurrence state one （the hydrate is in solid form in the reservoir） was selected for analysis. By using the known P-wave velocity and density as constraints, a reasonable S-wave velocity value （ranging from 400 to 1100 m s 1 and for a hydrate layer of 1100 m s 1） can be obtained through multiple iterations. These methods and results provide new data and technical support for further research on hydrates and other geological features in the Shenhu area.

Mapping crustal S-wave velocity structure with SV-component receiver function method

In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver function, SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings. Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function. Moreover, the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion. The synthetic tests show that the convergence of SV-component receiver function inversion is faster than that of the radial receiver function inversion. As an example, we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method.

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java, Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas(39°N–41.5°N,115.5°E–119.5°E),we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion.We obtained some detailed information about the Tangshan earthquake region and its adjacent areas,including sedimentary thickness,Moho depth,and crustal and upper mantle S-wave velocity.Meanwhile,we also obtained the vP/vS structure along two sections across the Tangshan region.The results show that:(1)the Moho depth ranges from 30 km to 38 km,and it becomes shallower from Yanshan uplift area to North China basin;(2)the thickness of sedimentary layer ranges from 0 km to 3 km,and it thickens from Yanshan uplift region to North China basin;(3)the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure,while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust.Compared with the upper crust,the heterogeneity of the middle and lower crust is more obvious;(4)the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust,and the low vS and high vP/vS beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.

基于分形理论的多孔介质渗透注浆机制

为研究宾汉姆浆液在多孔介质地层的渗透注浆机制,基于分形理论、毛细管模型和宾汉姆浆液流变方程,推导了基于分形理论的宾汉姆浆液渗透扩散表观速度公式和多孔介质渗透注浆球形扩散公式;利用已有的理论模型和室内渗透注浆试验成果对推导的理论公式进行了对比、分析和验证。研究结果表明:与传统的宾汉姆浆液渗透注浆扩散公式相比,基于分形理论的多孔介质渗透注浆扩散公式获得的浆液扩散半径更接近于室内试验成果。该研究成果可为实际多孔介质地层注浆工程提供一定的理论支撑。

固体探空火箭迭代制导方法研究

固体探空火箭在实施某类载荷远程投递任务时,要求接近落点的精度在千米级,由于火箭一般在大气层内的主动段实施有效的制导控制,其制导时间短,气动影响无法忽视,实现落点精度控制难度很大。针对该问题提出了一种适用于固体探空火箭在大气层内的迭代制导方法,通过在线计算剩余视速度和视位置增量,对固体探空火箭落点进行迭代计算,形成制导控制指令。仅在主动段进行制导控制,对再入段进行气动补偿,使得固体火箭落点控制达到了较高的精度。数值仿真分析表明:该方法落点精度较传统的摄动制导方法提高了1个数量级,体现出该迭代制导方法的优越性。

地震超前探测技术在煤矿中的应用研究

采用新一代高速、高分辨矿井地震仪在某矿掘进面进行超前探测研究。由于采空区、断层等地质异常破坏了煤层的完整性,与原煤(未开采煤层)相比会表现出不同的物理性质,主要体现在介质的密度和地震波的传播速度这2方面,针对这一物理现象可利用地震波法进行煤矿地震勘探研究。煤层巷道地震超前探测技术是利用在煤层内部激发的地震波遇到异常反射体(如断层、陷落柱、巷道和采空区等)时产生的反射波,对反射波反射点进行反射归位可得到异常反射体的反射面,根据反射面的走向判断反射体特别是断层的走向等信息。利用地震超前探测技术原理,根据反射时距曲线的负视速度特征,采用了线性Radon变换技术进行上下行波分离,从而来提取反射波信息进行解释。对巷道前方异常的钻探验证表明,钻探结果与物探探测成果吻合。

微动探测方法和高密度电法在地灾调查评价中的综合应用

随着物探技术方法在地质灾害调查评价中应用越来越广泛,对其物探方法的选择即决定了探测效果的好坏。近年来,迅速发展起来的微动探测方法将环境噪声作为信号源,实现了将干扰信号变废为宝,具有抗干扰能力强、施工条件受限少等优势;高密度电法因采样密度高、叠加观测,具有观测精度和分辨率高的特点。本文从实际出发,将微动探测方法和高密度电法综合应用于延安宝塔区地质灾害风险调查评价项目中,工程钻探验证结果显示,该两种方法综合对比研究推断的结果准确可靠,说明了微动探测方法和高密度电法在地灾调查评价中综合应用的可行性和有效性,具有一定的借鉴意义。

The shallow Longmenshan crustal S-velocity structure of the fault zone using ambient noise tomography of a seismic dense array

The Longmenshan fault zone(LMSF),characterized by complex structures and strong seismicity,is located at the junction between the eastern margin of the Tibetan Plateau and the north-western Sichuan basin.Since the Wenchuan earthquake on May 12,2008,abundant studies of the formation mechanism of earthquakes along the LMSF were performed.In this study,a short-period dense seismic array deployed across the LMSF was applied by ambient noise tomography.Fifty-two 3-D seismic instruments were used for data acquisition for 26 days.We calculated the empirical Green's functions(EGFs)between different station-pairs and extracted 776 Rayleigh-wave dispersion curves between 2 and 7 s.And then,we used the direct-inversion method to obtain the fine shallow crustal S-wave velocity structure within 6 km depth in the middle section of the Longmenshan fault zone and nearby areas.Our results show that the sedimentary layer(>5 km)exists in the northwest margin of Sichuan Basin with a low S-wave velocity(~1.5-2.5 km/s)which is much thicker than that beneath the Longmenshan fault zone and the Songpan-Garze block.The high-velocity structures with clear boundaries below the middle of Longmenshan fault zone(~2-4 km)and the Songpan-Garze block(~4.5-6 km)probably reveal the NW-SE distribution patterns of both the Pengguan complex and the high-density belt hidden in the northwest of the Pengguan complex.And the obviously high-velocity anomalies observed at the depth of^1-2 km in the southeastern margin of the Songpan-Garze block can be considered as the Laojungou granites.Our results provide a high-resolution shallow velocity structure for detailed studies of the Longmenshan fault zone.

Interpretation of S-Wave Data from Fanshi-Taipusiqi DSS Profile and Analysis of Correlation between Deep Structural Characteristics and Seismicity

By processing S-wave data from the Fanshi-Huai’an-Taipusiqi DSS profile,which is a three-component,wide-angle reflection/refraction profile,and in the light of the results from P-wave interpretation,two-dimensional(2-D)structures of the crust and upper mantle are presented,including S-wave velocity Vs and the physical parameter of medium-Poisson’s ratio a.Taking other geological and geophysical information into account,and with reference to the results from petrophysical experiments at home and abroad,we carried out interpretation and inference with respect to deep crustal structure,tectonics,and lithologic characters.It has been concluded that in the upper and middle crust,a values are mostly not greater than 0.25,and rocks,which generally assume brittle,are mainly composed of granite; the rocks in the lower layer of the upper crust between Yangyuan-Huai’an containing inorganic CO2 itself releases carbon; for the rocks in the lower crust and crust-mantle transitional zone,which are comparatively

Collimation of High Velocity Molecular Outflows

Statistics of apparent collimation factors for molecular outflows catalogued is presented.Intrinsic collimation factor is calculated with the project effects removed.There is collimation difference between high and low mass outflows.It may concern the stellar sources in the flow regions,the outflow driving process,and the evolution time.

不同地震动输入方向下非对称大跨悬索桥地震响应分析

为了解非对称悬索桥在不同地震动输入方向下的结构地震响应特性,指导该类桥梁的抗震设计,以广东云浮西江特大桥——跨径布置(202+738)m的非对称双塔双跨悬索桥为背景,进行地震响应分析。选取视波速150~3000 m/s的8个工况,采用大质量法模拟多点激励,对桥梁进行非线性时程分析,研究非对称悬索桥在一致激励下和考虑不同地震动输入方向及不同视波速下的结构地震响应规律。结果表明:对于非对称悬索桥,一致激励下,桥跨质量大的一侧构件内力、位移响应较大,位移较内力更为显著,应依据非对称结构地震响应特性,在结构两侧设计不同参数的减隔震装置;行波效应对非对称大跨悬索桥结构关键位置的内力、位移响应影响规律不同,视波速的选取对非对称结构地震响应影响显著,结构不同侧构件在相同视波速下的地震响应存在差异;在不同地震动输入方向下非对称结构的地震响应规律不同,同一构件响应值存在显著差异,建议在非对称桥梁设计时,考虑地震动输入方向的影响,按最不利方向进行设计。