频率域2.5-D井间波形层析成像及其实际应用

从频率域3-D声波波动方程出发,结合井间观测方式的特点,基于Tarantola广义反演理论,提出了一种频率域2.5-D井间波形层析成像方法.数值模型试验结果表明：该方法对薄层厚度的分辨能力能达到约主频波长的1/4,且分辨率显著高于走时层析成像,尤其垂直分辨率有实质改善.模拟资料的抗噪试验表明：在信噪比为0.8的情况下,随机噪声对波形层析成像的影响较小;而相干噪声对全波形层析成像的影响显著,特别是初至波附近的强振幅干扰影响更为严重.井间实际资料的试处理结果表明：波形层析成像能很好地刻画井间介质的分布情况与储层连通性,对于油藏开发阶段的方案实施具有指导意义.     

马边-永善地震带构造形式及地震特征

根据该地褶皱形态、断裂特征、应力场分布等资料研究,指出马边～永善及其两侧的构造形迹自西而东有一个相应强～弱变化过程。推测其力源主要来自西部喜马拉雅构造域的侧向推挤,所形成的马边-永善地震带是处于马边-昭通-曲靖断裂构造带的一部分,从而成为川西经向构造带的东部边界,也是青藏高原东南边缘与四川盆地西南边缘的浅部界线     

具有低频超宽衰减域的部分埋入式表面波屏障

针对低频Rayleigh表面波,设计了部分埋入式工字形截面周期波屏障。利用有限元方法计算了结构的频散曲线,分析了带隙的形成机理,讨论了屏障埋入土体深度和截面参数对带隙的影响,在此基础上设计了具有低频超宽衰减域的梯度及分段梯度波屏障并计算了其传输谱。结果表明:周期波屏障存在较宽带隙,板埋入深度和端部尺寸是影响带隙的关键参数,通过参数调节可实现不同频段Rayleigh波的调控。工字形变截面波屏障比等截面具有更优越的隔震性能且节省材料。梯度及分段梯度波屏障显著拓宽了衰减域的频率范围,对1.5~20 Hz范围内的Rayleigh表面波实现了全覆盖,用小尺寸控制了大波长。     

接地导线源电磁场全域有效趋肤深度

为了更好地推进广域电磁法的发展和应用,本文以接地长导线源为例,研究了可控源电磁场全场域的有效趋肤深度.利用频域电偶极源在均匀半空间产生电磁场的闭合表达式,计算了不同电磁场分量定义的有效趋肤深度,并讨论了在不同频率、不同偏移距、不同电导率情况下,有效趋肤深度的变化特性.根据有效趋肤深度随偏移距的变化特征以及与平面波趋肤深度之间的关系,利用多项式拟合的办法在五个不同频率范围内给出了适用于全场域的有效趋肤深度快速估算公式.研究结果表明：不同分量定义的有效趋肤深度是不同的,但是它们随偏移距、频率、电导率等参数的变化趋势是类似的,而且在一定范围内都趋近于平面波趋肤深度.理论模型的研究表明,有效趋肤深度可以作为测量参数选取和数据解释工作的参考依据.     

Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.     

基于神经网络的时间域直升机电磁数据电导率深度成像

采用G-S变换以及高斯数值积分法,形成了时间域直升机的航空电磁响应正演样本集,分析了飞机测量过程中吊舱高度变化对电磁响应的影响,并将吊舱高度的变化等效成电导率为零的假层厚度的变化,以去除高度计等的影响.以假层半空间模型为基础,研究了基于人工神经网络的电导率深度成像算法,通过分析两个三层模型的电导率深度成像结果得出,神经网络方法计算时间域航空电磁探测的视电导率精度较高,特别是对高阻层的视电导率计算.     

Application of perfectly matched layers in the transient analysis of dam–reservoir systems

The transient analysis of dam–reservoir systems by employing perfectly matched layers has been investigated. In previous studies, boundary conditions of the PML region in the reservoir have been neglected. In this paper, they are incorporated completely in the formulation. Moreover, a technique is introduced to involve the effect of incident waves caused by vertical ground motions at the reservoir bottom in the analysis. Performing several numerical experiments indicates that applying boundary conditions of the PML domain and utilizing the proposed method for vertical excitation cases reduce the computational cost significantly and make the PML method a very efficient approach for the transient analysis of dam–reservoir systems.     

On a modal damping identification model for non-classically damped linear building structures subjected to earthquakes

A simple modal damping identification model developed by the present authors for classically damped linear building frames is extended here to the non-classically damped case. The modal damping values are obtained with the aid of the frequency domain modulus of the roof-to-basement transfer function and the resonant frequencies of the structure (peaks of the transfer function) as well as the modal participation factors and mode shapes of the undamped structure. The assumption is made that the modulus of the transfer function of the non-classically damped structure matches the one of the classically damped structure in a discrete manner, i.e., at the resonant frequencies of that function modulus. This proposed approximate identification method is applied to a number of plane building frames with and without pronounced non-classical damping under different with respect to their frequency content earthquakes and its limitations and range of applicability are assessed with respect to the accuracy of both the identified damping ratios and that of the seismic structural response obtained by classical mode superposition and use of those identified modal damping ratios.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.