A lumped parameter model for time‐domain inertial soil‐structure interaction analysis of structures on pile foundations

The paper presents a lumped parameter model for the approximation of the frequency‐dependent dynamic stiffness of pile group foundations. The model can be implemented in commercial software to perform linear or nonlinear dynamic analyses of structures founded on piles taking into account the frequency‐dependent coupled roto‐translational, vertical, and torsional behaviour of the soil‐foundation system. Closed‐form formulas for estimating parameters of the model are proposed with reference to pile groups embedded in homogeneous soil deposits. These are calibrated with a nonlinear least square procedure, based on data provided by an extensive non‐dimensional parametric analysis performed with a model previously developed by the authors. Pile groups with square layout and different number of piles embedded in soft and stiff soils are considered. Formulas are overall well capable to reproduce parameters of the proposed lumped system that can be straightforwardly incorporated into inertial structural analyses to account for the dynamic behaviour of the soil‐foundation system. Some applications on typical bridge piers are finally presented to show examples of practical use of the proposed model. Results demonstrate the capability of the proposed lumped system as well as the formulas efficiency in approximating impedances of pile groups and the relevant effect on the response of the superstructure.     

Amplitude variation with incident angle and azimuth inversion for Young's impedance,Poisson's ratio and fracture weaknesses in shale gas reservoirs

By analogy with P- and S-wave impedances, the product of Young's modulus and density can be termed as Young's impedance, which indicates the rock lithology and brittleness of unconventional hydrocarbon reservoirs. Poisson's ratio is also an effective indicator of rock brittleness and fluid property of unconventional reservoirs, and fracture weaknesses indicate the fracture properties (fracturing intensity and fracture fillings) in fracture-induced unconventional reservoirs. We aim to simultaneously estimate the Young's impedance, Poisson's ratio and fracture weaknesses from wide-azimuth surface seismic data in a fracture-induced shale gas reservoir, and use the horizontal transversely isotropic model to characterize the fractures. First, the linearized PP-wave reflection coefficient in terms of Young's impedance, Poisson's ratio, density and fracture weaknesses is derived for the case of a weak-contrast interface separating two weakly horizontal transversely isotropic media. In addition, an orthorhombic anisotropic case is also discussed in this paper. Then a Bayesian amplitude variation with incident angle and azimuth scheme with a model constraint is used to stably estimate Young's impedance, Poisson's ratio and fracture weaknesses with only PP-wave azimuthal seismic data. The proposed approach is finally demonstrated on both synthetic and real data sets with reasonable results.     

Multi-trace basis-pursuit seismic inversion for resolution enhancement

Seismic inversion is an important tool that transfers interface information of seismic data to formation information, which renders the seismic data easily understood by geologists or petroleum engineers. In this study, a novel multi-trace basis-pursuit inversion method based on the Bayesian theory is proposed to enhance the vertical resolution and overcome the lateral instability of inversion results between different traces occasionally seen in the traditional trace-by-trace basis-pursuit inversion method. The Markov process is initially introduced to describe the relationship between adjacent seismic traces and their correlation, which we then close couple in the equation of our new inversion method. A recursive function is further derived to simplify the inversion process by considering the particularity of the coefficient matrix in the multi-trace inversion equation. A series of numerical-analysis and field data examples demonstrates that both the traditional and the new methods for P-wave impedance inversion are helpful in enhancing the resolution of thin beds that are usually difficult to discern from original seismic profiles, thus highlighting the importance of acoustic-impedance inversion for thin bed interpretation. Furthermore, in addition to yielding thin bed inversion results with enhanced lateral continuity and high vertical resolution, our proposed method is robust to noise and cannot be easily contaminated by it, which we verify using both synthetic and field data.     

基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法

针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果.     

非均匀正交各向异性特征参数地震散射反演方法研究

在长波长假设条件下,水平层状地层中发育一组垂直排列的裂缝构成了等效正交各向异性介质.各向异性参数与裂缝弱度参数的估算有助于非均匀各向异性介质的各向异性特征描述,而弹性逆散射理论是非均匀介质参数反演的有效途径.基于地震散射理论,我们首先推导了非均匀正交介质中纵波散射系数方程,并通过引入正交各向异性特征参数,提出了一种新颖的正交各向异性方位弹性阻抗参数化方法.为了提高反演的稳定性与横向连续性,我们发展了贝叶斯框架下的正交各向异性方位弹性阻抗反演方法,同时考虑了柯西稀疏约束正则化和平滑模型约束正则化,最终使用非线性的迭代重加权最小二乘策略实现了各向异性特征参数的稳定估算.模型和实际资料处理表明,反演结果与测井解释数据相吻合,证明了该方法能够稳定可靠地从方位叠前地震资料中获取各向异性特征参数,减小参数估算的不确定性,为非均匀正交介质的各向异性预测提供了一种高可靠性的地震反演方法.     

薄互层干涉对叠前AVO属性的影响分析

常规AVO分析是以单一界面为前提,用来分析薄层的AVO响应特征存在误差。Brekhovski给出的层状传播矩阵方程考虑了层内波的反射和透射以及层厚等影响因素,更适合用于讨论薄互层的地震响应特征。对薄层进行Brekhovski正演,发现不同入射角下振幅与厚度关系不同,与叠后全叠加振幅相比,截距P属性和薄层厚度之间的线性关系更强,说明P属性比全叠加振幅更有利于薄层厚度预测。不同薄互层模型的Brekhovski正演结果表明:净厚度和泊松比是影响AVO曲线特征的两个主要因素,净厚度越小,泊松比对截距P和梯度G属性的影响越小;当泊松比一定时,净厚度和P、G属性呈近似线性关系,P属性与净厚度关系更敏感,可用于薄互层净厚度预测。应用实例证明在薄互层区P属性用于砂地比预测效果比较好,叠后振幅和叠前反演的弹性阻抗预测的净砂岩厚度精度较差。     

岩性反演技术在煤田三维地震勘探中的应用

以开滦集团范各庄矿某采区岩性勘探为例,介绍利用地震反演技术解决煤层顶底板岩性问题。该项技术是以测井资料作为约束条件,对地震资料进行反演,推算出波阻抗资料,进而计算煤层厚度,并对顶底板岩性作出推断。通过对部分资料验证,采用岩性反演技术解释的该采区5煤、7煤、9煤、12煤的深度、厚度及其顶底板岩性,与实际情况基本相符。     

Geological, geophysical and seismological criteria for local response evaluation in Bucharest urban area

Vrancea major intermediate-depth earthquakes produced extreme damage in Bucharest city, located at about 165 km epicenter distance. Our purpose is to investigate the influence of local geological conditions upon the seismic motion in Bucharest in case of large (M>7) Vrancea earthquakes. Two input data sets are used: (a) geological, geotechnical and geophysical information, including in situ measurements, and (b) acceleration recordings of Vrancea earthquakes. Local response evaluation based on first dataset is confirmed by the spectral analysis of the earthquake records. Two main features are outlined: non-stationarity of ground motion dynamic amplification from one event to other and inadequacy of limiting the investigation depth to uppermost 30 m to evaluate ground dynamic characteristics. Consequently (1) we cannot extrapolate the ground motion response determined for moderate and small earthquakes to anticipate the effects of the large Vrancea shocks and (2) the local response is controlled by the entire package of Quaternary deposits which are significantly deeper than 30 m depth beneath Bucharest Area.     

利用全方位P波属性进行裂缝检测的地震方法

胜利油田存在着一种特殊的裂隙性油气藏,即泥岩裂隙油气藏。由于该类裂隙储层的孔隙度很小,岩石物性参数变化不灵敏,并表现出很强的各向异性,因此,其勘探难度很大。到目前为止,国内外还没有一套成熟的地质、物探、测井及钻井等资料综合分析的裂隙方位、分布、密度的识别方法。本文在研究国外裂缝检测方法的基础上,提出了波阻抗随方位变化(IPVA)的新方法,结合罗家地区的泥岩裂隙地震地质特征,研究了利用多方位 P 波资料检测定向垂直裂缝的地震采集、处理和识别方法,对不同共中心点 CMP 位置的全方位 P 波资料在速度随方位变化(VVA)及振幅随方位角变化(RVA)研究的基础上,进行波阻抗随方位变化(IPVA)的研究,用于识别裂缝的分布、走向及密度。经罗家地区实际资料的应用,见到了初步的效果,证明该方法是潜力较大的高角度裂隙的定量检测方法。     

地球电磁响应函数的旋转不变量特征

论述了不变量为整体或全局所共有的约束关系与旋转不变量的地球电性结构研究中的重要性,根据电磁响应函数的旋转特性,提出了几个基本的旋转不变量,也由此导出了其它一些旋转为变量,并指出这些不变量与电性分布的关系,最后,简单给出了大是电磁阻抗张量旋转不变量的Ｍｏｈｒ圆应用。