断续节理岩体中地震波响应的数值与理论分析

为研究地震产生的应力波在断续节理岩体中的传播规律和应力分布趋势,首先,采用数值模拟方法分析应力波通过贯通节理的传播规律,并与已有理论研究结果进行对比,验证数值分析的准确性和适用性;然后,对应力波在断续节理岩体中的传播进行数值模拟,分析透射系数在水平方向的分布趋势以及不同节理连续性对波传播的影响,并结合波的衍射原理,给出定性的理论解释。结果表明:应力波通过断续节理时,节理的透射作用会使应力波振幅减小,引起波的衰减,岩桥的衍射作用则会使波阵面由平面变为曲面,波的传播方向发生改变,从而导致应力波振幅在水平方向的分布发生变化;应力波通过断续节理的透射系数与岩桥尺寸L_r和衍射角μ相关,当衍射角比较小时,透射系数主要受岩桥尺寸L_r的影响,当衍射角较大时,岩桥尺寸L_r和衍射角μ共同影响应力波在岩体中的传播。      

Seismic physical modeling and quality factor

Accurate Q parameter is hard to be obtained, but there is great difference between Q measurements from different measurement methods in seismic physical modelling. The influence factors, stability and accuracy of different methods are analyzed through standard sample experiment and the seismic physical modelling. Based on this, we proposed an improved method for improving accuracy of pulse transmission method, in which the samples with similar acoustic properties to the test sample are selected as the reference samples. We assess the stability and accuracy of the pulse transmission, pulse transmission insertion, and reflection wave methods for obtaining the quality factor Q using standard and reference samples and seismic physical modeling. The results suggest that the Q-values obtained by the pulse transmission method are strongly affected by diffraction and the error is 50% or greater, whereas the relative error of the improved pulse transmission method is about 10%. By using a theoretical diffraction correction method and the improved measurement method, the differences among the Q-measuring methods can be limited to within 10%.     

Simulation and analysis of GPR signal based on stochastic media model with an ellipsoidal autocorrelation function

Target detection using ground penetrating radar (GPR) is based on the contrast between the electrical parameters of the target and the background medium, such as dielectric permittivity, conductivity and permeability. The application mainly concentrates on the detection of the medium interface and the target shape. In any theoretical study, a simulation model is built with a homogeneous medium. However, real detection encounters heterogeneous media which might produce scattering and diffraction at electrical interfaces and distort the radar pulse shape and affect the detection resolution. In this paper, we build multi-scale random media model with an ellipsoidal autocorrelation function and use FDTD method to simulate the GPR signal response. We then estimate and analyze the arrival time, layer thickness, permittivity and the physics relation in different scale random models according to the S transform method and the transmission wave method. The results demonstrate that we can use GPR to obtain geophysical information of multi-scale heterogeneous media, and provide a foundation for real media detection and complex media inversion.     

One dimensional elastic wave propagation in a non-homogeneous conical rod

Summary A one dimensional problem concerning the wave propagation in a non-homogeneous conical rod having its cross-section directly proportional to the square of the distance from a chosen point is considered in this paper. One end of the rod is subjected to a pressure step while the other end is stress free and the pulse shape is examined when the propagation speed of the wave has a particular continuous distribution over the entire rod.     

Direct and indirect inversions

A bridge is highlighted between the direct inversion and the indirect inversion. They are based on fundamental different approaches: one is looking after a projection from the data space to the model space while the other one is reducing a misfit between observed data and synthetic data obtained from a given model. However, it is possible to obtain similar structures for model perturbation, and we shall focus on P-wave velocity reconstruction. This bridge is built up through the Born approximation linearizing the forward problem with respect to model perturbation and through asymptotic approximations of the Green functions of the wave propagation equation. We first describe the direct inversion and its ingredients and then we focus on a specific misfit function design leading to a indirect inversion. Finally, we shall compare this indirect inversion with more standard least-squares inversion as the FWI, enabling the focus on small weak velocity perturbations on one side and the speed-up of the velocity perturbation reconstruction on the other side. This bridge has been proposed by the group led by Raul Madariaga in the early nineties, emphasizing his leading role in efficient imaging workflows for seismic velocity reconstruction, a drastic requirement at that time.     

Numerical modeling of Gaussian beam propagation and diffraction in inhomogeneous media based on the complex eikonal equation

Gaussian beam is an important complex geometrical optical technology for modeling seismic wave propagation and diffraction in the subsurface with complex geological structure. Current methods for Gaussian beam modeling rely on the dynamic ray tracing and the evanescent wave tracking. However, the dynamic ray tracing method is based on the paraxial ray approximation and the evanescent wave tracking method cannot describe strongly evanescent fields. This leads to inaccuracy of the computed wave fields in the region with a strong inhomogeneous medium. To address this problem, we compute Gaussian beam wave fields using the complex phase by directly solving the complex eikonal equation. In this method, the fast marching method, which is widely used for phase calculation, is combined with Gauss–Newton optimization algorithm to obtain the complex phase at the regular grid points. The main theoretical challenge in combination of this method with Gaussian beam modeling is to address the irregular boundary near the curved central ray. To cope with this challenge, we present the non-uniform finite difference operator and a modified fast marching method. The numerical results confirm the proposed approach.     

Ground penetrating radar sounding of a temperate glacier; modelling of a multilayered medium1

A numerical simulation of electromagnetic propagation through a multilayered medium is performed in order to explain and interpret the signal received from the radar sounding of a temperate glacier. During the winter of 1990, several radar profiles were obtained on the Mont-de-Lans glacier in the French Alps with a ground penetrating radar which uses a phase modulation of the transmitted pulse by coded sequences. The pulse compression is obtained by applying the matched filter to the received signal, which provides a range-resolution of about 8 m in the ice. The profiles recorded on the temperate glacier do not show a single clear reflection from the ice-bedrock interface, but they exhibit widely distributed energy decreasing with depth. This may be due to the inhomogeneous inner structure of the temperate glacier and we use a simple model of a layered medium to compute a simulation of the propagation. Thus, partial reflection at each layer and scattering from a rough basal interface may explain the observed signal. A computer-based technique is used to locate on the data the bottom of the glacier in order to estimate the ice thickness. The results from the different radar profiles are consistent and are a good fit to the thickness which has been determined by other geophysical methods.     

哨声谱与出口点的关系

一、引言 哨声波是一种在空间等离子体中传播的色散电磁波,一般以自然界中的闪电作为激励源。利用地面哨声资料可推算出赤道面内哨声路径顶点处的电子浓度N_(eq)和管电子含量N_T等电离层参数。要准确提取这些信息就需知道哨声传播的路径参数,即入口点和出口点位置。由于部分哨声透出电离层后,在地一电离层波导中传播了相当长的距离才被     

WHERE IS ZERO TIME?*

Mis-ties are all-too-common results of seismic surveys made at the same place but at different times with different equipment or by different organizations. Even after removal of positioning or polarity errors, reflection times often appear to differ by several tens of milliseconds. Zero time appears to fluctuate. How can zero time differ on surveys with only minor differences in acquisition or processing? What can be done to identify the true zero time for each survey? The first step toward establishing zero time is to record the source pulse. It is well-known that the different sources currently used in reflection seismic prospecting (propane-oxygen explosions, compressed-air discharges, explosives, steam bubbles, mechanical implosions, vibrations, etc.) yield different pressure wavefronts as the input to the seismic reflection system. By recording this wavefront we capture the basic pulse shape and we establish the initial time delay. The second step is to process the recorded source pulse as if it were reflection data to establish the additional time and shape changes introduced by data processing. Then, display the recorded and processed source pulse as an auxiliary variable at the ends of the seismic section. From this display the interpreter can systematically establish the time shifts appropriate to each picked event. He can determine also whether the pick should be a peak or a trough. He can see why surveys which appear to tie for shallow reflections appear to mis-tie for deep reflections. The display of the processed source pulse constitutes a major interpretation aid which, in a readily useable form, increases the information content of the basic seismic section.     

A tale of two beams: an elementary overview of Gaussian beams and Bessel beams

An overview of two types of beam solutions is presented, Gaussian beams and Bessel beams. Gaussian beams are examples of non-localized or diffracting beam solutions, and Bessel beams are example of localized, non-diffracting beam solutions. Gaussian beams stay bounded over a certain propagation range after which they diverge. Bessel beams are among a class of solutions to the wave equation that are ideally diffraction-free and do not diverge when they propagate. They can be described by plane waves with normal vectors along a cone with a fixed angle from the beam propagation direction. X-waves are an example of pulsed beams that propagate in an undistorted fashion. For realizable localized beam solutions, Bessel beams must ultimately be windowed by an aperture, and for a Gaussian tapered window function this results in Bessel-Gauss beams. Bessel-Gauss beams can also be realized by a combination of Gaussian beams propagating along a cone with a fixed opening angle. Depending on the beam parameters, Bessel-Gauss beams can be used to describe a range of beams solutions with Gaussian beams and Bessel beams as end-members. Both Gaussian beams, as well as limited diffraction beams, can be used as building blocks for the modeling and synthesis of other types of wave fields. In seismology and geophysics, limited diffraction beams have the potential of providing improved controllability of the beam solutions and a large depth of focus in the subsurface for seismic imaging.     

Modified ray-mode (phase) theory: Understanding counter-wind propagation effects from atmospheric explosions

We have incorporated horizontal winds into ray-mode theory including the full spectrum of acoustic–gravity waves for a perfectly stratified, range-independent, steady-state model atmosphere for frequencies from 10^?4 to ～10 Hz. This approach has also been applied to a specific atmospheric propagation problem that has long defied a solution, namely counter-wind propagation arrivals at a location ～300 km up-wind of the source. Our modified ray-mode theory predicts reliable up-wind solutions, but only if small-scale sound speed fluctuations were added to the mean seasonal sound speed profiles. Since full-waveguide theory and modified ray-mode mode theory incorporate diffraction and scattering propagation effects, we have performed additional analyses to determine the mechanism through which these fluctuations produce the up-wind signals. We have concluded that the dominant mechanism is through diffraction due to the presence of semi-permanent turbulence and internal gravity waves located near the stratopause.     

Representation of near-fault pulse-type ground motions

Near-fault ground motions with long-period pulses have been identified as critical in the design of structures. To aid in the representation of this special type of motion, eight simple pulses that characterize the effects of either the fling-step or forward-directivity are considered. Relationships between pulse amplitudes and velocity pulse period for different pulses are discussed. Representative ratios and peak acceleration amplification can exhibit distinctive features depending on variations in pulse duration, amplitude and the selected acceleration pulse shape. Additionally, response spectral characteristics for the equivalent pulses are identified and compared in terms of fixed PGA and PGV, respectively. Response spectra are strongly affected by the duration of pulses and the shape of the basic pulses. Finally, dynamic time history response features of a damped SDOF system subjected to pulse excitations are examined. These special aspects of pulse waveforms and their response spectra should be taken into account in the estimation of ground motions for a project site close to a fault.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.     

Effective elastic properties of randomly fractured soils: 3D numerical experiments

This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results.     

Numerical simulation of the impact of sediment supply and streamflow variations on channel grain sizes and Chinook salmon habitat in mountain drainage networks

Climatically driven changes in streamflow and hillslope sediment supply could potentially alter stream surface grain size distribution patterns and thereby impact habitat for a number of threatened and endangered in‐stream fish species. Relatively little is known about hydrograph (shape, peak flow) influence or the relative importance of chronic and episodic hillslope inputs on channel conditions. To better understand these external drivers, we calculated sediment routing through a gravel‐bedded river network using a one‐dimensional (1D) bedload transport model. We calculated changes in grain sizes and estimated Chinook salmon habitat suitability caused by a dry year and an extreme flood hydrograph, and chronic (diffusive, overland flow) or pulse (landslide, debris flow) hillslope sediment supplies. To obtain accurate channel conditions, a relatively high reference Shields stress, representative of steep mountain streams, was needed. An extreme event flood without any hillslope sediment inputs caused widespread bed coarsening and a decrease in aquatic habitat. Chronic sediment input combined with this hydrograph eliminated any changes in grain size and habitat, although when combined with a dry year flow, caused systematic bed fining. The influence of a given hydrograph therefore highly depends on the hillslope sediment supply. Regardless of the flow hydrograph or sediment pulse timing, grain size distribution or location, pulse sediment inputs did not cause widespread grain size changes despite being 100 times the total chronic input volume. Widespread and continuous hillslope sediment inputs may influence channel grain sizes and aquatic habitat more than a single discrete sediment pulse. Depending on the magnitudes of flow hydrograph and sediment supply alterations, climate change may induce no differences in grain sizes or very dramatic changes with significant consequences for long‐term sustainability. Copyright © 2013 John Wiley & Sons, Ltd.     

龙门山高倾角逆断层结构与孕震机制

针对汶川8.0级地震的主破裂面是否以陡立倾角延伸至地壳深部的争议,我们的研究旨在明确发震断层形态随深度分布的二维结构特征, 即断层的倾角及相应段落的深度, 并在此基础上分析发震断层的孕震和发震机理. 利用子空间置信域非线性反演方法, 通过拟合近场的同震水准变形, 获得了分别对应于清平-北川和南坝-青川发震断层的二维弹性位错模型. 结果显示, 此次龙门山中央断裂带的发震断层系统存在明显的南北分区特征. 以北川-南坝为过渡带, 北川以南至清平的发震构造为二元结构, 包括两部分:一是浅部高倾角的逆断层, 倾角在70°~80°以上, 底部深度可达10~15 km, 同震位错主要发生在10 km以上深度, 平均位错超过6 m;二是深部缓倾角的逆断层, 反演得其倾角约25°, 底部可达30 km深度, 位错主要分布于断层的顶部和底部, 平均位错约4 m. 南坝以北的发震断层为单一结构的逆断层, 倾角约60°~70°, 逆断层位错分布于10 km深度以上, 平均位错小于2 m.余震分布和主震震源机制也支持清平以南发震断层为二元结构的推论. 有限元模拟显示,在二元结构的逆断层系统中, 通过提高断层面上正的库仑应力, 深部缓倾角的逆断层活动对浅部高倾角逆断层有明显的促震作用. 模拟还显示地壳缩短不是现今松潘地块地表垂直变形的主要原因, 垂直变形更可能反映了青藏高原东缘相对四川盆地的差异抬升. 对高倾角逆断层的库仑应力作用显示, 差异抬升对龙门山逆断层活动有重要的促进作用.     

Pulse-like versus non-pulse-like ground motion records: Spectral shape comparisons and record selection strategies

Pulse-like records are well recognized for their potential to impose higher demands on structures when compared with ordinary records. The increased severity of the structural response usually caused by pulse-like records is commonly attributed to the spectral increment around the pulse period. By comparing the building response to sets of spectrally equivalent pulse-like and ordinary records, we show that there are characteristics of pulse-like records beyond the shape of the acceleration response spectrum that affect the results of nonlinear dynamic analysis. Nevertheless, spectral shape together with the ratio of pulse period to the first-mode structural period, T_p/T₁, are confirmed as “sufficient” predictors for deformation and acceleration response metrics in a building, conditioned on the seismic intensity. Furthermore, the average spectral acceleration over a period range, AvgSA, is shown to incorporate to a good proxy for spectral shape, and together with T_p/T₁, form an efficient and sufficient intensity measure for response prediction to pulse-like ground motions. Following this latter route, we propose a record selection scheme that maintains the consistency of T_p with the hazard of the site but uses AvgSA to account for the response sensitivity to spectral shape.     

Diffraction imaging using an adaptive phase filter

As an ideal carrier of high-resolution information, seismic diffraction can be used to clarify and locate small-scale discontinuities or inhomogeneities in the subsurface. However, a diffraction is weak and thus be suppressed by the specular reflection. Furthermore, a diffraction would be destroyed by the conventional imaging method due to the polarity reversal of diffraction. In this paper, we analyse the behaviour of diffraction and reflection. For the image point on a horizontal or oblique reflector, the zone on both sides of the stationary point has the same energy after using a cosine weight function. Based on the behaviour, we propose the adaptive phase filter to adjust the polarity of the energy on both sides, and calculate it through the illumination angle and the reflector dip angle. This method avoids the calculation of the Fresnel zones and can further suppress residual reflection that disturb the diffraction images. Synthetic and field data applications show that the desired imaging results can be obtained by the proposed method. The test results demonstrate that the method is efficient in detecting small-scale discontinuities or inhomogeneities in the subsurface and can provide high-resolution information for seismic interpretation.     

基于振动台实验的结构损伤识别研究

选择美国加州大学圣地亚哥分校7层钢筋混凝土剪力墙足尺结构振动台实验,开展结构损伤识别研究,实验采用白噪声、环境振动和不同强度的地震动交替激发,记录地震动激发实验前后的结构反应。基于该记录计算和对比自振频率和振型曲率的变化、剪切波走时及其变化和结构层间位移角,分析发现一层和二层振型曲率较大,走时较长,走时变化也较大,现场检查发现一层和二层的破坏也较为严重,这些参数可用于识别结构损伤程度和定位损伤位置,而自振频率和层间位移角变化仅可反映出结构损伤程度,难以揭示结构损伤位置。