Construction of regional and local seismic anisotropic structures from wide-angle seismic data: crustal deformation in the southeast of China

We present a method to obtain spatial distributions of seismic anisotropy associated with regional stress and local faulting in the crust from wide-angle seismic data. The method contains three steps. The first step consists of obtaining radial- and transverse-component seismic sections using a pre-stack depth migration algorithm from the S-wave velocity model determined by conventional interpretation of picked intra-crustal seismic events. In the second step, we compute time delays between split shear-waves and polarizations of fast split shear-waves by minimizing the transverse-component seismic energy. The time delay and polarization in each layer are derived using a layer-stripping method. The final step is to estimate the average splitting parameters along the whole profile. Thus, the average time delay and polarization can be regarded as caused by the effects owing to regional structure and stress fields, whereas the residual values of the splitting parameters are considered to be related to local structures and local faulting. Our method allows us to construct multi-layer anisotropic images, which may later be interpreted in terms of intra-layer coupling/decoupling or deformation. We present results from a set of three-component seismic data acquired by a controlled source experiment in the southeast region of China. The results demonstrate that the average polarizations and time delays are consistent with the direction and strength of the stress field, and their lateral variations related to local anisotropy match the spatial distribution of surface faulting crossing the acquisition seismic profile.     

Ground-motion modelling at regional distances for earthquakes in a continental interior,II. Effect of focal depth,azimuth and attenuation

Multiple-mode surface-wave signals are used to model ground motion at distances of 50 to 500 km for an earthquake source in a continental interior. Motion on a thrust fault is used as the earthquake model. Theoretical ground-motion time histories are generated for this source for various focal depths, receiver azimuths and medium-attenuation models. A shallow source will generate greater values for the ground motion than the same source at a greater depth. Two anelastic attenuation models are considered, one appropriate to the central and eastern United States and the other to southern California. The effects of the difference in the attenuation models are seen at distances greater than 100 km for periods greater than 1.5 sec.     

Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.     

A moving magnetic dipole in a conductive medium

A thorough investigation of the role that the source velocity has in the spatial and temporal variation of the secondary electromagnetic energy scattered by the Earth is necessary because marine controlled-source electromagnetic geophysical surveys employ moving sources. A first step towards this goal is the analysis, for this type of survey, of the difference of the measured value of the secondary electromagnetic energy between two systems: one with a moving source and the other with a fixed source. The model that suffices to stress this kinematic aspect is a vertical magnetic dipole moving at a constant speed along a horizontal line in a homogeneous medium separated from two homogeneous half-spaces by horizontal boundaries both above and below the dipole. The results show that both the velocity and the relative displacement between the source and the medium may cause a measurable variation relative to the static condition. Therefore, it may be necessary to take them into account in geophysical interpretation and to adapt the concepts of time and frequency domain for electromagnetic systems with moving sources.     

甘肃及邻近地区的地壳应力场

本文利用40个地震震源机制解和3个小震综合解研究了甘肃及邻近地区的地壳应力场。结果表明,由甘肃西部到东南部,主压应力方向由北40°东逐渐变为北90°东,P轴仰角大多小于30°,其水平方向垂直于青藏高原东北边缘。反映了该区的力源主要来自于印度板块与欧亚板块的碰撞。一些地区局部应力场的变化对研究断层活动有重要意义。     

一维声波方程联合反演的分步正则化法

本文用一个纵波信息，对一维声波方程的速度和源函数进行联合反演，并考虑到声波方程的反问题是一个不适应问题，对源函数和波速分别和正则化法分步迭代求解，减少反问题的计算工作量，改善该问题的计算稳定性，为计算实际工程和岩性学问题供了一种方法。文中给出只用一个反问题补充条件同时进行多参数反演的公式，并对相应的数值算例进行分析和比较。     

An unconditionally stable hybrid pseudodynamic algorithm

There is a significant motivation to implement an unconditionally stable scheme in the pseudodynamic test method. As more complex experiments with many degrees of freedom are tested, explicit time integration methods limit the size of time step on the basis of the highest natural frequency of the system. This is true even though the response of the structure may be dominated by a few lower frequency modes. The limit on step size is undesirable because it physically increases the duration of a test, but more importantly, because the number of steps to completion increases and error propagation problems increase with the number of steps in a test. In addition, incremental displacements within each step become smaller, introducing the potential for problems associated with stress relaxation. An unconditionally stable algorithm allows the time step to be selected to give accurate response in the modes of interest without regard for higher mode characteristics.     

Comparison of source parameters estimated in the frequency and time domains for seismic events at the Rudna copper mine,Poland

Source parameters estimated in the frequency domain for 100 selected seismic events from the Rudna copper mine, with moment magnitude ranging from 1.4 to 3.6, were collected to study their scaling relations and to compare them with the parameters estimated in the time domain. The apparent stress and static stress drop, corrected for the limited bandwidth recording, increase slightly in a similar manner with increasing seismic moment. The ratio of apparent stress to static stress drop, a measure of radiation efficiency, is practically constant and its mean value is close to 0.1. For 37 seismic events, with moment magnitude between 1.9 and 3.4, source parameters were estimated in the time domain from relative source time functions, that displayed unilateral rupture propagation, and their rupture velocity could be estimated. It ranges from 0.23 to 0.80 of shear wave velocity and is almost independent of seismic moment. The fault length, estimated from the average source pulse width and rupture velocity, is clearly dependent on seismic moment and is smaller than the source radius estimated from the corner frequency on the average by about 25 percent. There is no correlation between the values of static stress drop estimated in the frequency and time domains, but the time domain stress drop is in general similar to that estimated in the frequency domain. The apparent stress increases with increasing rupture velocity, and the ratio of apparent stress to static stress drop seems also to depend on rupture velocity.     

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.     

应变能积累在地震安全性评价中的应用探讨

理论及实验研究显示,震源区的位置和最大震级不仅可能与构造和历史有关,更重要的是与应力状态包括其方向、大小、增加速率、集中位置等以及岩体自身性质,包括其强度、杨氏模量、应变大小及增长速率、裂缝(断层)大小和数量等有关,同时与作用方式有关.综合考虑各种因素,探讨了利用应变能积累确定强震位置和震级的方法:首先对利用形变空间特征变化判断强震震源区位置的方法进行了讨论;接着着重提出估算应变能的初步方法:一是由地表裂缝最大错距和裂缝(断层)长度比估算地震释放能量;二是由新构造运动分区块体体积、年平均应变速率及杨氏模量估算块体应变能积累;进而依据能量估算地震震级;最后以实例论述了应变能积累过程及特点对地震安全性评价的重要作用.     

2008年盈江地震序列的震源参数和震源机制相关系数研究

基于云南数字地震台网记录的2008年盈江地震序列的数字波形资料, 采用波谱分析方法和Brune震源模型, 得到盈江地震序列的震源参数。 利用两次地震事件的相同台站的震源谱参数(零频振幅)计算谱振幅相关系数, 据此对地震的震源机制进行聚类分组, 并收集和对比地震序列已知的震源机制解结果, 发现每组内震源机制解P轴的相关性较好, 且相关系数越大, P轴的方位角就越接近。 将盈江地区划分为三个研究区, 联合震源参数的应力降和聚类分组中每组的平均震源机制解结果研究地震序列发生过程中应力的释放水平和应力场方向变化特征。 结果表明: 不同阶段震源机制解类型的变化和转化特征一定程度上反映了孕震过程中区域应力场随时间的变化特征, 并且地震震源机制解类型在时间段上的集中并向区域构造应力场方向转换的现象可能是发生强震的标志。 震源机制解分组类型和对应类型的地震的应力降有一定依赖关系, 震源机制解类型反映的应力场与区域应力场接近的地震应力降高, 震源机制解类型反映的应力场与区域应力场差距较大的地震应力降普遍较低。     

编码电磁测深

利用逆重复M序列伪随机信号良好的自相关特性,提出相关检测地电传输特性的编码电磁勘探法.分析了编码电磁测深原理及参数提取方法.同时记录发送电流信号和多收发距电磁场响应,对源信号和场信号进行相关运算,解卷积分离接收系统响应后可得到大地的频率特性或时间特性,以此实现地电断面的精细探测.在时间域,大地冲激响应和阶跃响应含有丰富的地电信息.通过冲激响应的峰值时间或阶跃响应的晚期渐近值可估计地电阻率分布.基于层状模型的大地冲激响应和阶跃响应正演计算结果表明,编码电磁测深法对大埋深薄层目标体有精细的分辨能力.可应用于能源、矿产资源、水资源、环境地质及工程地质勘察,有良好的应用前景.     

利用广义叠加反演方法估算长岛震群应力降

地震应力降是表征震源特性的一个重要参数,通常由震源谱计算得到.但如何从大量地震的观测谱中校正路径和台站效应,准确地分离出震源谱并计算应力降,仍是一个难题.本文采用广义叠加反演方法通过分阶段叠加从观测谱中迭代分离出震源项、台站项和路径项,并利用分震级的叠加震源项与理论谱的拟合得到经验格林函数,再利用由经验格林函数校正后的真实震源谱估算出应力降.该方法不需要知道仪器响应、场地响应和传播路径上的衰减系数,并且不用事先假设震源具有自相似的恒应力降特征,计算效率很高,适用于反演大量地震的震源参数.我们将广义叠加反演方法应用于长岛震群,利用S波谱估算了1431个M_L≥1.0地震的应力降.结果表明,广义叠加反演方法获得的震源参数和用传统方法得到的结果一致;长岛震群地震的应力降较小(中值是0.3 MPa),变化范围很大(0.009～3.04MPa之间),属于低应力降事件集;同时应力降中值随地震矩增大而增大,偏离了自相似理论.应力降空间变化不均匀,高应力降事件主要集中在震群的端部(NW端),到最近断层距离的依赖性不强,在9、10km处的应力降中值略高于其他深度.根据应力降的变化特征...     

Estimating the source parameters for the M<Subscript><Emphasis Type="Italic">W</Emphasis></Subscript> 7.6 April 20, 2006 Olyutorskii earthquake from long period <Emphasis Type="Italic">P</Emphasis>-wave records made by the worldwide network

The source parameters of the M _W = 7.6 Olyutorskii earthquake were estimated using the moments of the slip rate function with degrees 1 and 2. The moments were estimated from broadband P-wave records at 52 stations of the worldwide network. The first step was to find a function S(t) for each station; this function is an apparent source time function, i.e., the P-wave slip as radiated by the source toward a station under consideration. The method of empirical Green’s functions was used to estimate S(t). The next step was to calculate the moments of S(t) of degrees 1 and 2 over time and to set up relevant equations to be solved by least squares for the unknown source moments. The horizontal linear source was used as a nonparametric model for calculating the source moments. Haskell’s parametric model was used for further interpretation of the source moments. The resulting estimates are as follows: the source centroid was 13–25 km southwest of the epicenter, the source was 105–120 km long, the source strike was 222°–228°, the rupture velocity was 2.7–3.0 km/s, and the total radiation duration was 24–27 s. These estimates indicate a bilateral rupture dominated by a southwestward sense of rupture propagation. The source characteristics are consistent with the aftershock area geometry and with the focal mechanism, as well as with surface breakage as observed by geologists in the field.     

2019年四川长宁MS6.0地震序列S波分裂变化特征

本文采用质点运动判别与偏振分析相结合的方法对2019年四川长宁MS6.0地震震源区10个台站于2013年4月25日至2019年7月31日记录到的波形数据进行S波分裂参数测定,其中9个台站获取4条以上有效S波分裂参数.结果表明,震源区各台站的S波分裂参数不仅在空间上存在分区特征,还会随时间发生改变.快波偏振方向在空间上的...     

Multiasperity fault model and the nature of short-period subsources

We suggest to consider the breaking of an asperity,i.e., a small contact patch between fault walls, as a typical subsource producing an elementary short-period radiation pulse from a source of a large earthquake. Based on the results ofDas andKostrov we propose formulas to describe amplitudes and spectra of acceleration for a multiasperity fault/source model. The stress drop over an asperity is determined in several ways; the estimates agree to give the average value of several hundred bar. Theoretical acceleration spectral shapes for the case of similar asperities agree with the observed ones, they reproduce such features as lower-frequencyf ¹ slope, peak, and high-frequency cutoff. The statistical stress-drop distribution over the population of asperities, and also the related distribution of peak accelerations are discussed. There distributions are found to be the power-law ones with exponent near to 2. This means that acceleration peaks are formed normally by breaking of individual asperities. We consider small earthquakes as produced by breaking of single asperities, this idea explains the observed correlation between the upper cutoff frequency of acceleration spectrum and the typical characteristic frequency of small earthquakes.     

基于Student′s t分布的不依赖子波最小二乘逆时偏移

最小二乘逆时偏移(Least-Squares Reverse Time Migration,LSRTM)与常规偏移相比具有更高的成像分辨率、振幅保真性及均衡性等优势,是当前研究的热点之一.震源子波的估计直接影响LSRTM结果的好坏,在实际情况下考虑到震源子波的空变特性,其估计十分困难.为了消除子波对LSRTM结果的影响,本文发展了基于卷积目标泛函的不依赖子波LSRTM算法.目标泛函由观测记录卷积模拟记录的参考道以及模拟记录卷积观测记录的参考道组成,由于观测子波和模拟子波在目标泛函的两项中同时存在,从而消除了子波的影响.此外,常用的基于L2范数拟合的LSRTM算法对噪声非常敏感,尤其是当地震数据中含有异常值时,常规LSRTM无法得到满意的结果.Student′s t分布相比L2范数具有更好的稳健性,本文将其推广到不依赖子波LSRTM中,提升了算法的稳健性,最后通过理论模型及实际资料试算验证了算法的有效性和对复杂模型的适应性.     

1995年陡河地震前后小震震源参数及其相互关系

在震源谱的多台多震尾波解法的基础上,应用Brune震源模型,求解1995年10月6日陡河ML5.9地震前后在其周围发生的48次小地震的震源因子和震源波谱,进而估算了地震矩、拐角频率和应力降,并讨论了它们之间的相互关系和在陡河地震前后随时间的变化．结果表明:小震震源因子表现出相当大的稳定性,震源因子的峰值频率和随频率的变化形态都很相似．发现震前一年左右曾多次出现应力降较高的事件．需要指出的是,考虑到应力降由地震矩和拐角频率计算得到以及地震矩和拐角频率的相关关系,实际上,这里的高应力降指的是拐角频率高于正常值的情形．      

大森-宇津定律的一种可能机制以唐山大地震为例

为了探讨大森-宇津定律的物理机制, 本文在余震区等效黏度远低于其外部, 且构造应力场在整个余震活动时间间隔内基本保持不变的假设条件下, 提出了一个开尔文黏弹性地震震源体概念模型． 该模型可用于模拟主震后断层蠕变和震源区应力调整触发的余震序列以及蠕变停止后余震终结、 介质恢复到弹性状态、 断层重新闭锁和积累下一次地震的整个过程． 有限元方法可用来计算非均匀黏弹性地震震源体模型中主震和每次余震所引起的应力场及其随时间的演化过程． 在此基础上, 采用开尔文黏弹性地震震源体概念模型和有限元方法模拟了1976年唐山M_S7.8地震余震序列． 结果表明: 经验的大森-宇津定律可以用开尔文黏弹性震源体模型来解释, 这意味着余震衰减的频度取决于蠕变的速率; 余震序列持续时间受控于震源体的黏度, 即黏度越大, 蠕变时间越长, 余震持续的时间也就越长．      

Using slowness and azimuth fluctuations as new observables for four‐dimensional reservoir seismic monitoring

Four‐dimensional imaging using geophysical data is of increasing interest in the oil and gas industries. While travel‐time and amplitude variations are commonly used to monitor reservoir properties at depth, their interpretation can suffer from a lack of information to decipher the parts played by different parameters. In this context, this study focuses on the slowness and azimuth angle measured at the surface using source and receiver arrays as complementary observables. In the first step, array processing techniques are used to extract both azimuth and incidence angles at the source side (departure angles) and at the receiver side (arrival angles). In the second step, the slowness and angle variations are monitored in a laboratory environment. These new observables are compared with traditional arrival‐time variations when the propagation medium is subject to temperature fluctuations. Finally, field data from a heavy‐oil permanent reservoir monitoring system installed onshore and facing steam injection and temperature variations are investigated. The slowness variations are computed over a period of 152 days. In agreement with Fermat's principle, strong correlations between the slowness and arrival‐time variations are highlighted, as well as good consistency with other techniques and field pressure measurements. Although the temporal variations of slowness and arrival time show the same features, there are still differences that can be considered for further characterization of the physical changes at depth.