The average upper mantle P wave velocity structure and lateral heterogeneity in the northern part of China and Mongolia are investigated by waveform inversion of broadband body waveform data recorded by CDSN and digit...The average upper mantle P wave velocity structure and lateral heterogeneity in the northern part of China and Mongolia are investigated by waveform inversion of broadband body waveform data recorded by CDSN and digital stations around China. The average model has a low P wave velocity lid (about 7.8~8.0 km·s -1 ) with thickness about 60 km, and two discontinuities with velocity jumps of 0.29 km·s -1 and 0.55 km·s -1 at depth of 410 km and 665 km respectively. In the Jungger basin, the P wave velocity of uppermost mantle is about 7.7 km·s -1 . The lid thickness (90~100 km) and velocity gradient (average gradient is greater than 0.005 5/s) are large. At the depth of 140 km the P wave velocity reaches to 8.2 km·s -1 . Near in Baikal, the lid is about 30 km thick with average P wave velocity of 8.00~8.05 km·s -1 .展开更多
The longitudinal wave term within Faraday’s law of electromagnetic induction (Faraday’s law) underwent recovery to ensure its suitability for theoretical derivation of the equation governing longitudinal electromagn...The longitudinal wave term within Faraday’s law of electromagnetic induction (Faraday’s law) underwent recovery to ensure its suitability for theoretical derivation of the equation governing longitudinal electromagnetic (LEM) waves. The revised Maxwell’s equations include the crucial parameters being the attenuation time constants of magnetic vortex potential and electric vortex potential generated by external electromagnetic field within the propagation medium. Specific expressions for them are obtained through theoretical analysis. Subsequently, a model for propagating magnetic P-wave generated by the superposition of a left-handed photo and a right-handed photon in a vacuum is formulated based on reevaluated total current law and revised Faraday’s law, covering wave equations, energy equation, as well as propagation mode involving mutual induction and conversion between scalar magnetic field and vortex electric field. Furthermore, through theoretical derivations centered around magnetic P-wave, evidence was presented regarding its ability to absorb huge free energy through the entangled interaction between zero-point vacuum energy field and the torsion field produced by the vortex electric field.展开更多
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...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.展开更多
Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal di...Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.展开更多
An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvia...An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.展开更多
Wave reflection and refraction in layered media is a topic closely related to seismology,acoustics,geophysics and earthquake engineering.Analytical solutions for wave reflection and refraction coefficients in multi-la...Wave reflection and refraction in layered media is a topic closely related to seismology,acoustics,geophysics and earthquake engineering.Analytical solutions for wave reflection and refraction coefficients in multi-layered media subjected to P wave incidence from the elastic half-space are derived in terms of displacement potentials.The system is composed of ideal fluid,porous medium,and underlying elastic solid.By numerical examples,the effects of porous medium and the incident wave angle on the dynamic pressures of ideal fluid are analyzed.The results show that the existence of the porous medium,especially in the partially saturated case,may significantly affect the dynamic pressures of the overlying fluid.展开更多
A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property para...A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property parameters, and strong anisotropy, it is very difficult to explore for them. So far, there is no set of mature methods for recognition of direction, distribution, and density of the fractures by an integrated analysis of geologic, geophysical, well log, drilling data, and etc. This paper presents a new method for acoustic impedance variation with azimuth (IPVA), based on existing fracture detection methods. Seismic acquisition, processing, and recognition techniques were developed for detecting directional vertical fractures using multi-azimuth P wave data in combination with the seismic and geological features of shale fractures in the Luojia area. The IPVA research is carried out for recognizing the distribution, strike, and density of fractures based on the study of velocity variation with azimuth (VVA) and amplitude variation with azimuth (AVA) for full azimuth P wave data at different CMP positions. Through practical application in the Luojia area, primary results have been obtained which verifies that the IPVA method provides good potential for quantitative detection of parallel, high angle, shale fractures.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional a...This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.展开更多
The effect of seawater on vertical ground motions is studied via a theoretical method and then actual offshore ground motion records are analyzed using a statistical method. A theoretical analysis of the effect of sea...The effect of seawater on vertical ground motions is studied via a theoretical method and then actual offshore ground motion records are analyzed using a statistical method. A theoretical analysis of the effect of seawater on incident plane P and SV waves at ocean bottom indicate that on one hand, the affected frequency range of vertical ground motions is prominent due to P wave resonance in the water layer if the impedance ratio between the seawater and the underlying medium is large, but it is greatly suppressed if the impedance ratio is small; on the other hand, for the ocean bottom interface model selected herein, vertical ground motions consisting of mostly P waves are more easily affected by seawater than those dominated by SV waves. The statistical analysis of engineering parameters of offshore ground motion records indicate that:(1) Under the infl uence of softer surface soil at the seafl oor, both horizontal and vertical spectral accelerations of offshore motions are exaggerated at long period components, which leads to the peak spectral values moving to a longer period.(2) The spectral ratios(V/H) of offshore ground motions are much smaller than onshore ground motions near the P wave resonant frequencies in the water layer; and as the period becomes larger, the effect of seawater becomes smaller, which leads to a similar V/H at intermediate periods(near 2 s). These results are consistent with the conclusions of Boore and Smith(1999), but the V/H of offshore motion may be smaller than the onshore ground motions at longer periods(more than 5 s).展开更多
This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffract...This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.展开更多
An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The ...An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.展开更多
Although the Zoeppritz equation is suitable for a single interface in a thick deposit, it has some limitations for composite reflection waves from both the floor and the roof of coal seams. Based on the ray model, the...Although the Zoeppritz equation is suitable for a single interface in a thick deposit, it has some limitations for composite reflection waves from both the floor and the roof of coal seams. Based on the ray model, the relationship of the overall reflection coefficient of composite reflection P waves, from coal seam versus incidence angle (AVO), is dis- cussed. The result shows that: 1) the overall reflection coefficient of composite reflection waves from coal seams is a negative value and is determined mainly by the lithology of roof and floor, which is different from the reflection coeffi- cient of a single interface; 2) if the incidence angle ranges from 0° to 6°, the reflection coefficient of composite waves of a coal seam does not change with the incidence angle and 3) if the incidence angle ranges from 6–60° , the reflection coefficient increases monotonically.展开更多
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line sourc...Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.展开更多
A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities signifi...A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities significantly am-plify the surface ground motion nearby. It is suggested that the effect of subways on ground motion should be con-sidered when the subways are planned and designed.展开更多
文摘The average upper mantle P wave velocity structure and lateral heterogeneity in the northern part of China and Mongolia are investigated by waveform inversion of broadband body waveform data recorded by CDSN and digital stations around China. The average model has a low P wave velocity lid (about 7.8~8.0 km·s -1 ) with thickness about 60 km, and two discontinuities with velocity jumps of 0.29 km·s -1 and 0.55 km·s -1 at depth of 410 km and 665 km respectively. In the Jungger basin, the P wave velocity of uppermost mantle is about 7.7 km·s -1 . The lid thickness (90~100 km) and velocity gradient (average gradient is greater than 0.005 5/s) are large. At the depth of 140 km the P wave velocity reaches to 8.2 km·s -1 . Near in Baikal, the lid is about 30 km thick with average P wave velocity of 8.00~8.05 km·s -1 .
文摘The longitudinal wave term within Faraday’s law of electromagnetic induction (Faraday’s law) underwent recovery to ensure its suitability for theoretical derivation of the equation governing longitudinal electromagnetic (LEM) waves. The revised Maxwell’s equations include the crucial parameters being the attenuation time constants of magnetic vortex potential and electric vortex potential generated by external electromagnetic field within the propagation medium. Specific expressions for them are obtained through theoretical analysis. Subsequently, a model for propagating magnetic P-wave generated by the superposition of a left-handed photo and a right-handed photon in a vacuum is formulated based on reevaluated total current law and revised Faraday’s law, covering wave equations, energy equation, as well as propagation mode involving mutual induction and conversion between scalar magnetic field and vortex electric field. Furthermore, through theoretical derivations centered around magnetic P-wave, evidence was presented regarding its ability to absorb huge free energy through the entangled interaction between zero-point vacuum energy field and the torsion field produced by the vortex electric field.
基金Key Project Process Mechanism and Prediction of Geological Hazards (2001CB711005-1-3) and State Key Basic Research Project Mechanism and Prediction of Continental Earthquakes (G1998040702). sponsored by the Ministry of Science and Techno
基金The special project of Detection of Haikou City Earthquake Active Faults from the Tenth Five-year Plan of China Earthquake Administration (0106512)Joint Seismological Science Foundation of China (105086)CAS Key Laboratory of Marginal Sea Geology (MSGL0503).
文摘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.
文摘Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.
基金State Natural Science Foundation of China (No.59878032).
文摘An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.
基金National Natural Science Foundation of China Under Grant No.50309005National Key Basic Research and Development Program Under Grant No.2002CB412709
文摘Wave reflection and refraction in layered media is a topic closely related to seismology,acoustics,geophysics and earthquake engineering.Analytical solutions for wave reflection and refraction coefficients in multi-layered media subjected to P wave incidence from the elastic half-space are derived in terms of displacement potentials.The system is composed of ideal fluid,porous medium,and underlying elastic solid.By numerical examples,the effects of porous medium and the incident wave angle on the dynamic pressures of ideal fluid are analyzed.The results show that the existence of the porous medium,especially in the partially saturated case,may significantly affect the dynamic pressures of the overlying fluid.
文摘A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property parameters, and strong anisotropy, it is very difficult to explore for them. So far, there is no set of mature methods for recognition of direction, distribution, and density of the fractures by an integrated analysis of geologic, geophysical, well log, drilling data, and etc. This paper presents a new method for acoustic impedance variation with azimuth (IPVA), based on existing fracture detection methods. Seismic acquisition, processing, and recognition techniques were developed for detecting directional vertical fractures using multi-azimuth P wave data in combination with the seismic and geological features of shale fractures in the Luojia area. The IPVA research is carried out for recognizing the distribution, strike, and density of fractures based on the study of velocity variation with azimuth (VVA) and amplitude variation with azimuth (AVA) for full azimuth P wave data at different CMP positions. Through practical application in the Luojia area, primary results have been obtained which verifies that the IPVA method provides good potential for quantitative detection of parallel, high angle, shale fractures.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.
基金the 973 Program under Grant No.2011CB013601Basic Science Research Foundation of IEM,CEA under Grant No.2011B02+1 种基金National Natural Science Foundation of China under Grant No.51238012,9121530113,50938006International Science&Technology Cooperation Program of China under Grant No.2012DFA70810
文摘The effect of seawater on vertical ground motions is studied via a theoretical method and then actual offshore ground motion records are analyzed using a statistical method. A theoretical analysis of the effect of seawater on incident plane P and SV waves at ocean bottom indicate that on one hand, the affected frequency range of vertical ground motions is prominent due to P wave resonance in the water layer if the impedance ratio between the seawater and the underlying medium is large, but it is greatly suppressed if the impedance ratio is small; on the other hand, for the ocean bottom interface model selected herein, vertical ground motions consisting of mostly P waves are more easily affected by seawater than those dominated by SV waves. The statistical analysis of engineering parameters of offshore ground motion records indicate that:(1) Under the infl uence of softer surface soil at the seafl oor, both horizontal and vertical spectral accelerations of offshore motions are exaggerated at long period components, which leads to the peak spectral values moving to a longer period.(2) The spectral ratios(V/H) of offshore ground motions are much smaller than onshore ground motions near the P wave resonant frequencies in the water layer; and as the period becomes larger, the effect of seawater becomes smaller, which leads to a similar V/H at intermediate periods(near 2 s). These results are consistent with the conclusions of Boore and Smith(1999), but the V/H of offshore motion may be smaller than the onshore ground motions at longer periods(more than 5 s).
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN), 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.
基金National Natural Science Foundation of China Under Grant No.50378063Excellent Young Teacher Program of Ministry of Education
文摘An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.
基金Project 40574058 Supported by the National Natural Science Foundation of China
文摘Although the Zoeppritz equation is suitable for a single interface in a thick deposit, it has some limitations for composite reflection waves from both the floor and the roof of coal seams. Based on the ray model, the relationship of the overall reflection coefficient of composite reflection P waves, from coal seam versus incidence angle (AVO), is dis- cussed. The result shows that: 1) the overall reflection coefficient of composite reflection waves from coal seams is a negative value and is determined mainly by the lithology of roof and floor, which is different from the reflection coeffi- cient of a single interface; 2) if the incidence angle ranges from 0° to 6°, the reflection coefficient of composite waves of a coal seam does not change with the incidence angle and 3) if the incidence angle ranges from 6–60° , the reflection coefficient increases monotonically.
基金supported by National Natural Science Foundation of China (50978183)
文摘Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green's fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.
基金Supported by National Natural Science Foundation of China (50378063), Excellent Young Teachers Program of MOE and SRF for ROCS, MOE.
文摘A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities significantly am-plify the surface ground motion nearby. It is suggested that the effect of subways on ground motion should be con-sidered when the subways are planned and designed.