Background Right bundle branch block （RBBB） may present as slurred or notched S wave in lead V1. However, slurred or notched S wave may also represent slow conduction in the myocardium. Methods We retrospectively an...Background Right bundle branch block （RBBB） may present as slurred or notched S wave in lead V1. However, slurred or notched S wave may also represent slow conduction in the myocardium. Methods We retrospectively analyzed the QRS patterns in leads VgR to V5R in 7 patients with a slurred or notched S wave in lead V1. Results In the leads V3R to VSR, 6 patients showed incomplete or complete RBBB and 1 patient slurred or notched S wave. Conclusions In the majority of ECGs in a small patient series with slurred or notched S wave in lead V1, QRS morphology indicating incomplete or complete RBBB was present in leads V3R to VSR. A finding of fragmented QRS in these leads may indicate slow conduction in the myocardium.展开更多
Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the ...Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.展开更多
The polarization direction of fast wave and the delay time between fast and slow wave were measured for two earthquake sequences occurred continuously on 21 July （M=6.2） and 16 October （M=6.1） in Dayao, Yunnan in ...The polarization direction of fast wave and the delay time between fast and slow wave were measured for two earthquake sequences occurred continuously on 21 July （M=6.2） and 16 October （M=6.1） in Dayao, Yunnan in 2003 using cross-correlation coefficient method, after determining the high-resolution hypocentral locations of the earthquake sequences using the double-difference earthquake location algorithm. The results indicated that ① The phenomena of S wave splitting are obvious in the two earthquake sequences, and the average polarization directions of fast wave in most stations are almost consistent with regional maximum horizontal compressive stress direction except the station Santai. There are bimodal fast directions in the polarization directions at station Santai and the mean polarization direction is N80°E, indicating an inconsistent phenomenon referred to regional maximum horizontal compressive stress direction. ② There is no apparent relation between delay time and focal depth in the sequences, but the polarization direction show different character in different delay time range. ③ The comparison of S wave splitting results in the two earthquake sequences show that the polarization direction in M=6.2 earthquake sequence is more scattered and its average fast direction is 20° larger than that of M=6.1 sequence, and the delay times between two sequences show a little difference. ④ The spatial variation in S wave splitting polarization direction may be due to the stress disturbance imposed by the M=6.2 and the M=6.1 mainshocks on regional background stress field.展开更多
In this paper, according to the relationship of observation spectrum with source spectrum, travel-path attenuation factor and site effect, we use the digital seismic data of moderate-small earthquakes to invert the so...In this paper, according to the relationship of observation spectrum with source spectrum, travel-path attenuation factor and site effect, we use the digital seismic data of moderate-small earthquakes to invert the soft rock site effect and S wave inelastic attenuation under the constraint of site effect, and determine the parameters of path and site. The parameters obtained by this method exclude the topographic effect of seismic stations, and truly reflect the seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region. The average soft rock site effect is about 1.5 times in the frequency domain of 2-4 Hz and the S wave quality factor is obtained as Qs=278f0.346, which can be directly applied to the stochastic modeling of ground motions in the studied region.展开更多
We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yi...We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yimen, and Luquan, were selected and their attenuations of P and S waves were measured by using the extended coda-normalization method. The six stations were classified into three regions according to their location, that is, Baoshan area, Yongsheng-Lijiang-Heqing (YLH) area and Luquan-Yimen (LY) area. The values of QS-1 and QP-1 are expressed as QS-1=0.00867 f-0.86, QP-1=0.01155 f-0.93, QS-1=0.01824 f-0.92, QP-1=0.02288 f-0.92, and QS-1=0.01647 f-0.91, QP-1=0.02826 f-0.97 in Baoshan, YLH, LY areas respectively. The attenuation of YLH, LY are closer to each other, however, Baoshan area is apparently lower. Comparing attenuation in the three areas with other areas of the world using the same method, it is suggested that the attenuations of P and S waves in YLH and LY areas are close to Kanto of Japan, but much higher than southeast of South Korea. The QS-1 and QP-1 in Baoshan area are slightly higher than southeastern South Korea. Furthermore, the results indicate that our QS-1 in Yunnan area is close to others by analyzing the coda attenuation.展开更多
Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid cha...Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.展开更多
Theory and practice have confirmed that the propagation of converted transverse waves in azimuthally anisotropic media results in shear-wave splitting,causing differences in travel time and reflection amplitude betwee...Theory and practice have confirmed that the propagation of converted transverse waves in azimuthally anisotropic media results in shear-wave splitting,causing differences in travel time and reflection amplitude between S1/S2 waves,which reduces the imaging quality and vertical resolution of converted transverse wave.Therefore,eliminating shear-wave splitting is one of the important steps in the imaging processing of converted transverse wave.Due to the limitations of actual data acquisition and signal-tonoise ratio,it is difficult to directly determine a set of orthogonal four components for shear-wave splitting analysis,which cannot achieve good correction processing results.To this end,starting from Alford rotation,an orthogonal four component equation system is derived,and the least squares method is used to stably solve the high signal-to-noise ratio orthogonal four components.Based on this,the crack orientation and fast/slow wave delay are estimated through azimuth and time-delay scanning,and the shear-wave splitting correction and fast/slow wave separation are completed by layer striping from shallow to deep.The synthetic data and actual data calculations show that this method can achieve good azimuthal anisotropy correction,separation of S1/S2 waves,and accurate imaging.展开更多
Teleseismic datasets at the Shidao Seismographic Station, located in the northwestern South China Sea, are used to determine the earth anisotropy and the vertical distribution pattern of the shear wave velocity by inv...Teleseismic datasets at the Shidao Seismographic Station, located in the northwestern South China Sea, are used to determine the earth anisotropy and the vertical distribution pattern of the shear wave velocity by inversion approaches. The rotated correction function is applied to analyzing high quality SeS records from five earthquakes at distance of 25°-35° to obtain shear wave splitting parameters of the lithosphere. The result from the deepest earthquake among the five events indicates that the polarization of the fast shear wave is N94°E, which means the direction of extensional stress or the moving of the upper mantle mass in Xisha Islands is nearly west to east and confirms that the crust in this region is a transitional one and the driving force beneath the crust is from the moving mass consistent with the Eurasian plate. The anisotropy effective thickness is estimated about 100 km based on the time delay of 1.3 s between the fast and slow shear waves. The receiver function is applied to analyzing high quality P wave records from nine earthquakes at distance of 20°- 60° to obtain the vertical distribution pattern of shear wave velocity beneath the station. The result indicates that the crust could be divided into three layers： the uppermost crust （5 km above） is a velocity gradient zone consisting of several small layers, where the shear wave velocity increases from 1.5 to 3.5 km/s gradually; the 5 - 16 km depth interval also consiss of several small layers of which the mean velocity is about 3.8 km/s; and the lower crust （ 16.0 - 26. 5 km） is an obvious low velocity layer with a velocity of about 3.6 km/s. The buried depth of the Moho discontinuity is 26.5 kin, the mean velocity of the layers beneath the Moho is about 4.7 km/s and there is an obvious low velocity layer just beneath the Moho. Moreover, analysis of the arrival time of converted waves and the swinging variation of velocity around the initial model suggests that smaller layers in the model maybe are not reliable but the low velocity layer between 16 and 26.5 km maybe is the real one that implies the plasticity of the lower crust.展开更多
Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be g...Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.展开更多
The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic ...The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic stations and instrumental response. The path effect of seismic wave propagation, site response of seismic stations and instrumental response must be taken out in the study of source property with seismic data. The path effect of seismic wave propagation （seismic attenuation） involves an important influential factor, the anelastic attenuation of medium, which is measured with quality factor Q, apart from geometric attenuation with the distance. As a basic physical parameter of the Earth medium, Q value is essential for quantitative study of earthquakes and source property （e.g. determination of source parameters）, which is widely used in earthquake source physics and engineering seismology.展开更多
Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are in...Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.展开更多
Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in ...Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.展开更多
On the basis of S wave information from Tai′an Xinzhou DSS profile and with reference to the results from P wave interpretation, the 2 D structures , including S wave velocity V s, ratio γ between V...On the basis of S wave information from Tai′an Xinzhou DSS profile and with reference to the results from P wave interpretation, the 2 D structures , including S wave velocity V s, ratio γ between V p and V s; and Poisson′s ratio σ , are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the above results. Upwarping mantle districts (V s≈4.30 km/s)and sloping mantle districts (V s≈4.50 km/s) of the profile with velocity difference about 4% at the top of upper mantle are divided according to the differences of V s , γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of Niujiaqiao Dongwang high angle ultra crustal fault zone is reaffirmed; the properties of low and high velocity blocks(zones) including the crust mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps, thermal material upwells through the high angle fault, the primary hypocenter was located at the crust mantle juncture 30.0～33.0 km deep, and additional stress excited the M S=6.8 and M S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao Dongwang, the earthquakes took place around the high angle ultra crustal fault and centered in the brittle media and rock strata with low γ and low σ values.展开更多
In order to understand the site soil response of the Xiangtang borehole seismic array under real strong ground motion, reveal the site response, verify the technique of borehole exploration, and improve the precision ...In order to understand the site soil response of the Xiangtang borehole seismic array under real strong ground motion, reveal the site response, verify the technique of borehole exploration, and improve the precision of in-situ test and laboratory test, this paper presents a new approach, which is composed of two methods. One is the layered site seismic response method, whose layer transform matrix is always real. The other is a global-local optimization technique, which uses genetic algorithm （GA）-simplex method. An inversion of multi-component waveforms of P, SV and SH wave is carried out simultaneously. By inverting the records of three moderate and small earthquakes obtained from the Xiangtang borehole array （2^#） site, the soil dynamic characteristic parameters, including P velocity, damping ratio and frequency-dependent coefficient b, which has not been given in previous literatures, are calculated. The results show that the soil S wave velocity of the Xiangtang 2^# borehole is generally greater than that obtained from the 1994 in-situ test, and is close to the velocity of the 3^# borehole, which is more than 200 m away from the 2^# borehole. Meanwhile, perceptible soil nonlinear behavior under peak ground motion of about 60×10^-2 m/s^2 is detected by the inversion analysis. The presented method can be used for studying the soil response of other borehole array sites.展开更多
In this paper, based on the propagation theorics of seismic waves in anisotropic medium and in cracked two-phase medium, the constitutive relations and dynamic equations of the propagation of seismic waves in Cracked ...In this paper, based on the propagation theorics of seismic waves in anisotropic medium and in cracked two-phase medium, the constitutive relations and dynamic equations of the propagation of seismic waves in Cracked twophase anisotropic medium with fourfold rotation symmetry have been derived, and the preliminary theoretical analysis have been made for plane wave as an example.展开更多
In this paper, we coupled the Quantum Mechanics conventional Schrödinger’s equation, for the particles, with the Maxwell’s wave equation, in order to study the potential’s role on the conversion of the ele...In this paper, we coupled the Quantum Mechanics conventional Schrödinger’s equation, for the particles, with the Maxwell’s wave equation, in order to study the potential’s role on the conversion of the electromagnetic field energy to mass and vice versa. We show that the dissipation (“conductivity”) factor and the particle implicit proper frequency are both related to the potential energy. We have also derived a new expression for the Schrödinger’s Equation considering the potential energy into this equation not as an ad hoc term, but also as an operator (Hermitian), which has the scalar potential energy as a natural eigenvalue of this operator.展开更多
The northeastern margin of the South China Sea (SCS), developed from continental rifting and breakup, is usually thought of as a non-volcanic margin. However, post-spreading volcanism is massive and lower crustal high...The northeastern margin of the South China Sea (SCS), developed from continental rifting and breakup, is usually thought of as a non-volcanic margin. However, post-spreading volcanism is massive and lower crustal high-velocity anomalies are widespread, which complicate the nature of the margin here. To better understand crustal seismic velocities, lithology, and geophysical properties, we present an S-wave velocity (VS) model and a VP/VS model for the northeastern margin by using an existing P-wave velocity (VP) model as the starting model for 2-D kinematic S-wave forward ray tracing. The Mesozoic sedimentary sequence has lower VP/VS ratios than the Cenozoic sequence;in between is a main interface of P-S conversion. Two isolated high-velocity zones (HVZ) are found in the lower crust of the continental slope, showing S-wave velocities of 4.0–4.2 km/s and VP/VS ratios of 1.73–1.78. These values indicate a mafic composition, most likely of amphibolite facies. Also, a VP/VS versus VP plot indicates a magnesium-rich gabbro facies from post-spreading mantle melting at temperatures higher than normal. A third high-velocity zone (VP : 7.0–7.8 km/s;VP/VS: 1.85–1.96), 70-km wide and 4-km thick in the continent-ocean transition zone, is most likely to be a consequence of serpentinization of upwelled upper mantle. Seismic velocity structures and also gravity anomalies indicate that mantle upwelling/ serpentinization could be the most severe in the northeasternmost continent-ocean boundary of the SCS. Empirical relationships between seismic velocity and degree of serpentinization suggest that serpentinite content decreases with depth, from 43% in the lower crust to 37% into the mantle.展开更多
In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equ...In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.展开更多
文摘Background Right bundle branch block （RBBB） may present as slurred or notched S wave in lead V1. However, slurred or notched S wave may also represent slow conduction in the myocardium. Methods We retrospectively analyzed the QRS patterns in leads VgR to V5R in 7 patients with a slurred or notched S wave in lead V1. Results In the leads V3R to VSR, 6 patients showed incomplete or complete RBBB and 1 patient slurred or notched S wave. Conclusions In the majority of ECGs in a small patient series with slurred or notched S wave in lead V1, QRS morphology indicating incomplete or complete RBBB was present in leads V3R to VSR. A finding of fragmented QRS in these leads may indicate slow conduction in the myocardium.
基金State Key Basic Development and Programming Project Mechanism and Prediction of Continental Strong Earthquakes (G1998040705).
文摘Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.
基金National Program on Key Basic Projects (2004CB418406), Program for the Tenth Five-Year Plan of China(2004BA601B01-04-03)and Joint Seismological Science Foundation of China (606042).
文摘The polarization direction of fast wave and the delay time between fast and slow wave were measured for two earthquake sequences occurred continuously on 21 July （M=6.2） and 16 October （M=6.1） in Dayao, Yunnan in 2003 using cross-correlation coefficient method, after determining the high-resolution hypocentral locations of the earthquake sequences using the double-difference earthquake location algorithm. The results indicated that ① The phenomena of S wave splitting are obvious in the two earthquake sequences, and the average polarization directions of fast wave in most stations are almost consistent with regional maximum horizontal compressive stress direction except the station Santai. There are bimodal fast directions in the polarization directions at station Santai and the mean polarization direction is N80°E, indicating an inconsistent phenomenon referred to regional maximum horizontal compressive stress direction. ② There is no apparent relation between delay time and focal depth in the sequences, but the polarization direction show different character in different delay time range. ③ The comparison of S wave splitting results in the two earthquake sequences show that the polarization direction in M=6.2 earthquake sequence is more scattered and its average fast direction is 20° larger than that of M=6.1 sequence, and the delay times between two sequences show a little difference. ④ The spatial variation in S wave splitting polarization direction may be due to the stress disturbance imposed by the M=6.2 and the M=6.1 mainshocks on regional background stress field.
基金National Natural Science Foundation of China(50468003)State 973 Project(2002 CB412706).
文摘In this paper, according to the relationship of observation spectrum with source spectrum, travel-path attenuation factor and site effect, we use the digital seismic data of moderate-small earthquakes to invert the soft rock site effect and S wave inelastic attenuation under the constraint of site effect, and determine the parameters of path and site. The parameters obtained by this method exclude the topographic effect of seismic stations, and truly reflect the seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region. The average soft rock site effect is about 1.5 times in the frequency domain of 2-4 Hz and the S wave quality factor is obtained as Qs=278f0.346, which can be directly applied to the stochastic modeling of ground motions in the studied region.
基金The specialized fund for pre-study of national key basic study (2002CCD01700)national key project of science and technology from China Ministry of Science and Technology during the 10th Five-year Plan (2004BA601B01-04-03)
文摘We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yimen, and Luquan, were selected and their attenuations of P and S waves were measured by using the extended coda-normalization method. The six stations were classified into three regions according to their location, that is, Baoshan area, Yongsheng-Lijiang-Heqing (YLH) area and Luquan-Yimen (LY) area. The values of QS-1 and QP-1 are expressed as QS-1=0.00867 f-0.86, QP-1=0.01155 f-0.93, QS-1=0.01824 f-0.92, QP-1=0.02288 f-0.92, and QS-1=0.01647 f-0.91, QP-1=0.02826 f-0.97 in Baoshan, YLH, LY areas respectively. The attenuation of YLH, LY are closer to each other, however, Baoshan area is apparently lower. Comparing attenuation in the three areas with other areas of the world using the same method, it is suggested that the attenuations of P and S waves in YLH and LY areas are close to Kanto of Japan, but much higher than southeast of South Korea. The QS-1 and QP-1 in Baoshan area are slightly higher than southeastern South Korea. Furthermore, the results indicate that our QS-1 in Yunnan area is close to others by analyzing the coda attenuation.
基金supported by the National Science and Technology Major Project(No.2016ZX05047-002-001)
文摘Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.
基金supported by Sinopec Science and Technology Research Project(P23098,P25152).
文摘Theory and practice have confirmed that the propagation of converted transverse waves in azimuthally anisotropic media results in shear-wave splitting,causing differences in travel time and reflection amplitude between S1/S2 waves,which reduces the imaging quality and vertical resolution of converted transverse wave.Therefore,eliminating shear-wave splitting is one of the important steps in the imaging processing of converted transverse wave.Due to the limitations of actual data acquisition and signal-tonoise ratio,it is difficult to directly determine a set of orthogonal four components for shear-wave splitting analysis,which cannot achieve good correction processing results.To this end,starting from Alford rotation,an orthogonal four component equation system is derived,and the least squares method is used to stably solve the high signal-to-noise ratio orthogonal four components.Based on this,the crack orientation and fast/slow wave delay are estimated through azimuth and time-delay scanning,and the shear-wave splitting correction and fast/slow wave separation are completed by layer striping from shallow to deep.The synthetic data and actual data calculations show that this method can achieve good azimuthal anisotropy correction,separation of S1/S2 waves,and accurate imaging.
文摘Teleseismic datasets at the Shidao Seismographic Station, located in the northwestern South China Sea, are used to determine the earth anisotropy and the vertical distribution pattern of the shear wave velocity by inversion approaches. The rotated correction function is applied to analyzing high quality SeS records from five earthquakes at distance of 25°-35° to obtain shear wave splitting parameters of the lithosphere. The result from the deepest earthquake among the five events indicates that the polarization of the fast shear wave is N94°E, which means the direction of extensional stress or the moving of the upper mantle mass in Xisha Islands is nearly west to east and confirms that the crust in this region is a transitional one and the driving force beneath the crust is from the moving mass consistent with the Eurasian plate. The anisotropy effective thickness is estimated about 100 km based on the time delay of 1.3 s between the fast and slow shear waves. The receiver function is applied to analyzing high quality P wave records from nine earthquakes at distance of 20°- 60° to obtain the vertical distribution pattern of shear wave velocity beneath the station. The result indicates that the crust could be divided into three layers： the uppermost crust （5 km above） is a velocity gradient zone consisting of several small layers, where the shear wave velocity increases from 1.5 to 3.5 km/s gradually; the 5 - 16 km depth interval also consiss of several small layers of which the mean velocity is about 3.8 km/s; and the lower crust （ 16.0 - 26. 5 km） is an obvious low velocity layer with a velocity of about 3.6 km/s. The buried depth of the Moho discontinuity is 26.5 kin, the mean velocity of the layers beneath the Moho is about 4.7 km/s and there is an obvious low velocity layer just beneath the Moho. Moreover, analysis of the arrival time of converted waves and the swinging variation of velocity around the initial model suggests that smaller layers in the model maybe are not reliable but the low velocity layer between 16 and 26.5 km maybe is the real one that implies the plasticity of the lower crust.
基金State Natural Science Foundation (49874021).Contribution No. 01FE2002, Institute of Geophysics, China Seismological Bureau.
文摘Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.
基金Early-stage study project of the national key foundation research (2002CCD01700) and the key project of Yunnan Province during the 10th Five-year plan.
文摘The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic stations and instrumental response. The path effect of seismic wave propagation, site response of seismic stations and instrumental response must be taken out in the study of source property with seismic data. The path effect of seismic wave propagation （seismic attenuation） involves an important influential factor, the anelastic attenuation of medium, which is measured with quality factor Q, apart from geometric attenuation with the distance. As a basic physical parameter of the Earth medium, Q value is essential for quantitative study of earthquakes and source property （e.g. determination of source parameters）, which is widely used in earthquake source physics and engineering seismology.
文摘Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.
文摘Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.
文摘On the basis of S wave information from Tai′an Xinzhou DSS profile and with reference to the results from P wave interpretation, the 2 D structures , including S wave velocity V s, ratio γ between V p and V s; and Poisson′s ratio σ , are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the above results. Upwarping mantle districts (V s≈4.30 km/s)and sloping mantle districts (V s≈4.50 km/s) of the profile with velocity difference about 4% at the top of upper mantle are divided according to the differences of V s , γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of Niujiaqiao Dongwang high angle ultra crustal fault zone is reaffirmed; the properties of low and high velocity blocks(zones) including the crust mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps, thermal material upwells through the high angle fault, the primary hypocenter was located at the crust mantle juncture 30.0～33.0 km deep, and additional stress excited the M S=6.8 and M S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao Dongwang, the earthquakes took place around the high angle ultra crustal fault and centered in the brittle media and rock strata with low γ and low σ values.
基金Joint Seismological Science Foundation of China (604036 and 105034)
文摘In order to understand the site soil response of the Xiangtang borehole seismic array under real strong ground motion, reveal the site response, verify the technique of borehole exploration, and improve the precision of in-situ test and laboratory test, this paper presents a new approach, which is composed of two methods. One is the layered site seismic response method, whose layer transform matrix is always real. The other is a global-local optimization technique, which uses genetic algorithm （GA）-simplex method. An inversion of multi-component waveforms of P, SV and SH wave is carried out simultaneously. By inverting the records of three moderate and small earthquakes obtained from the Xiangtang borehole array （2^#） site, the soil dynamic characteristic parameters, including P velocity, damping ratio and frequency-dependent coefficient b, which has not been given in previous literatures, are calculated. The results show that the soil S wave velocity of the Xiangtang 2^# borehole is generally greater than that obtained from the 1994 in-situ test, and is close to the velocity of the 3^# borehole, which is more than 200 m away from the 2^# borehole. Meanwhile, perceptible soil nonlinear behavior under peak ground motion of about 60×10^-2 m/s^2 is detected by the inversion analysis. The presented method can be used for studying the soil response of other borehole array sites.
文摘In this paper, based on the propagation theorics of seismic waves in anisotropic medium and in cracked two-phase medium, the constitutive relations and dynamic equations of the propagation of seismic waves in Cracked twophase anisotropic medium with fourfold rotation symmetry have been derived, and the preliminary theoretical analysis have been made for plane wave as an example.
文摘In this paper, we coupled the Quantum Mechanics conventional Schrödinger’s equation, for the particles, with the Maxwell’s wave equation, in order to study the potential’s role on the conversion of the electromagnetic field energy to mass and vice versa. We show that the dissipation (“conductivity”) factor and the particle implicit proper frequency are both related to the potential energy. We have also derived a new expression for the Schrödinger’s Equation considering the potential energy into this equation not as an ad hoc term, but also as an operator (Hermitian), which has the scalar potential energy as a natural eigenvalue of this operator.
基金South China Sea Institute of Oceanology (SCSIO) for providing R/V Shiyan-2 to carry out this experiment,sponsored by Oceanographic Research Vessel Sharing Plan (NORC2016-08) of National Natural Science Foundation of Chinafunded by National Natural Science Foundation of China (Grant Nos. 41776057, 41761134051, 91858213, 41730532 and 91428039)
文摘The northeastern margin of the South China Sea (SCS), developed from continental rifting and breakup, is usually thought of as a non-volcanic margin. However, post-spreading volcanism is massive and lower crustal high-velocity anomalies are widespread, which complicate the nature of the margin here. To better understand crustal seismic velocities, lithology, and geophysical properties, we present an S-wave velocity (VS) model and a VP/VS model for the northeastern margin by using an existing P-wave velocity (VP) model as the starting model for 2-D kinematic S-wave forward ray tracing. The Mesozoic sedimentary sequence has lower VP/VS ratios than the Cenozoic sequence;in between is a main interface of P-S conversion. Two isolated high-velocity zones (HVZ) are found in the lower crust of the continental slope, showing S-wave velocities of 4.0–4.2 km/s and VP/VS ratios of 1.73–1.78. These values indicate a mafic composition, most likely of amphibolite facies. Also, a VP/VS versus VP plot indicates a magnesium-rich gabbro facies from post-spreading mantle melting at temperatures higher than normal. A third high-velocity zone (VP : 7.0–7.8 km/s;VP/VS: 1.85–1.96), 70-km wide and 4-km thick in the continent-ocean transition zone, is most likely to be a consequence of serpentinization of upwelled upper mantle. Seismic velocity structures and also gravity anomalies indicate that mantle upwelling/ serpentinization could be the most severe in the northeasternmost continent-ocean boundary of the SCS. Empirical relationships between seismic velocity and degree of serpentinization suggest that serpentinite content decreases with depth, from 43% in the lower crust to 37% into the mantle.
基金financially supported by the Natural Science Foundation of China(No.41774133)the Open Funds of SINOPEC Key Laboratory of Geophysics(No.wtyjy-wx2017-01-04)National Science and Technology Major Project of the Ministry of Science and Technology of China(No.2016ZX05024-003-011)
文摘In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.
